• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    耕作后土壤水分入滲參數(shù)BP預報模型

    2018-09-10 04:22:37武雯昱樊貴盛
    人民黃河 2018年2期
    關鍵詞:耕作土壤水分含水率

    武雯昱 樊貴盛

    摘要:根據(jù)黃土高原區(qū)耕作后土壤水分入滲資料所建立的基于Kostiakov三參數(shù)經(jīng)驗模型公式的BP神經(jīng)網(wǎng)絡預報模型,對土壤水分入滲參數(shù)進行預測。預測結(jié)果表明:利用土壤基本理化參數(shù)耕作層(0~20cm)土壤含水率及其密度、砂粒含量、粉粒含量、有機質(zhì)含量預報入滲模型參數(shù)具有可行性。入滲模型參數(shù)入滲指數(shù)。、入滲系數(shù)k的平均相對誤差約為1%;穩(wěn)滲率.fo的平均相對誤差約為8%,在可接受范圍內(nèi);90min累計入滲量與I90平均相對誤差小于3%。模型整體預測精度較高,實現(xiàn)了利用土壤基本理化參數(shù)預報土壤水分入滲參數(shù),可為黃土高原區(qū)耕作后土壤的灌水技術參數(shù)的確定提供理論與技術支撐。

    關鍵詞:拼作土壤;Kostiakov三參數(shù)經(jīng)驗模型;入滲參數(shù);BP模型;耕作層

    中圖分類號:S152.7;TV93 文獻標志碼:A doi ;10.3969/j.issn.1000-1379.2018.02.032

    土壤耕作對農(nóng)業(yè)種植業(yè)的產(chǎn)量具有決定性影響。土壤在耕作之后,疏松程度大大提高,其團粒結(jié)構(gòu)得到良好的恢復,耕作過程還創(chuàng)造了適度的孔隙比例,改善了土壤肥力,從而可以形成高產(chǎn)土壤。同時,耕作對土壤的入滲有顯著的促進作用,Abrisqueta等[1]在地中海杏林地區(qū)進行試驗研究后發(fā)現(xiàn),耕作后的土壤能夠顯著增加入滲,減少徑流,水蝕得到有效緩解。鄭子成等[2]研究發(fā)現(xiàn)耕作能夠使地表粗糙度增大,從而使徑流減小,入滲量加大。入滲是土壤中水分循環(huán)的重要環(huán)節(jié),目前,關于土壤水分入滲已有大量的研究成果[3~6],然而針對耕作后土壤水分入滲的研究卻鮮有報道,尤其是在有關耕作后土壤水分入滲參數(shù)預測方面的研究。土壤水分入滲參數(shù)是確定灌溉技術參數(shù)和指導農(nóng)業(yè)節(jié)水灌溉的重要依據(jù)。目前,土壤水分入滲參數(shù)有入滲儀法和大田灌水法兩種獲取方法,而這兩種方法相對復雜,因此筆者試圖通過易獲得的土壤常規(guī)理化參數(shù)預測土壤水分入滲參數(shù)值。

    常見的土壤水分入滲理論經(jīng)驗模型有Green-Ampt入滲模型、Kostiakov二參數(shù)模型、Kostiakov三參數(shù)模型Horton入滲模型,Philip入滲模型等。對于耕作后的土壤,大量試驗表明Kostiakov三參數(shù)模型的擬合效果較好[7],可以此建立入滲參數(shù)與土壤常規(guī)參數(shù)間的非線性關系。BP模型在表達非線性模糊關系時有很好的非線性映射能力,可以很好地逼近任意非線性函數(shù)。因此,本文擬建立耕作土壤的基于Kostiakov三參數(shù)模型的土壤水分入滲參數(shù)BP預報模型,利用已獲得的土壤常規(guī)理化參數(shù)(土壤密度、土壤質(zhì)地、土壤含水量及土壤有機質(zhì)含量)預測土壤水分入滲參數(shù),從而確定農(nóng)業(yè)灌水技術參數(shù),為耕作土壤的節(jié)水灌溉提供理論支撐。

    1 試驗區(qū)土壤條件與方法

    1.1 試驗區(qū)土壤條件

    試驗區(qū)位于黃土高原區(qū),土壤以草甸土與褐土為主,從地形地貌上涵蓋山區(qū)、丘陵、沖擊平原以及洪積平原。本次土壤水分入滲的試驗選取樣本的土質(zhì)為翻松土,耕作層土壤的砂粒含量變化范圍為25.99%~74.08%,粉粒含量變化范圍為13.41%~56.09%,黏粒含量變化范圍為0.89%~22.90%。耕作層(0~20cm)含水率的變化范圍為8.60%~31.99%,其干密度變化范圍為1.0991~1.5870g/cm3;犁底層(20~40cm)土壤含水率變化范圍為3%~35%,其土壤干密度變化范圍為1.38~1.78g/cm3。代表性樣本土壤基本理化參數(shù)見表1。

    1.2 試驗方法

    本研究試驗旨在通過大田試驗了解耕作后土壤水分入滲特性及其主導影響因素,建立土壤入滲參數(shù)預報模型。試驗針對影響土壤水分入滲特性的主要因素而提出設計,分為大田試驗和室內(nèi)試驗。

    1.2.1 大田試驗

    運用積水入滲方法測定土壤水分入滲過程,在試驗田形成地表含水率約等于飽和含水率的一維垂直人滲條件,其邊界條件屬于灌溉入滲模型。試驗設備是內(nèi)環(huán)直徑26.0cm、外環(huán)直徑64.4cm、內(nèi)外環(huán)高度均為25cm的雙套環(huán)入滲儀。將該設備預埋于試驗點,人滲過程中分時段記錄土壤水分入滲的時間及對應的人滲量,當入滲達到穩(wěn)定狀態(tài)時停止試驗。土壤入滲試驗一般在60min時達到穩(wěn)定狀態(tài),為了增強數(shù)據(jù)的可靠性以及考慮到數(shù)據(jù)的統(tǒng)一性,將試驗入滲時間定為90min。

    1.2.2 室內(nèi)試驗

    室內(nèi)試驗包括土壤密度、土壤質(zhì)地、土壤含水率及有機質(zhì)的測定。

    土壤密度的測定采用環(huán)刀法,用標準規(guī)格的環(huán)刀在土壤剖面取土、稱量,烘干后計算單位體積土的質(zhì)量,通常用相對體積質(zhì)量來表示[8]。

    土壤含水率的測定主要采用烘干法,人工鉆土取得相應土層的土壤,稱取待測土壤,烘箱烘干,稱量干土質(zhì)量,計算其含水率通常用質(zhì)量含水率或體積含水率表示[9]。

    土壤質(zhì)地的測定利用粒度分析儀。土壤風干、研磨過篩后,用粒度分析儀測定。

    土壤有機質(zhì)的測定采用重鉻酸鉀容量法。加熱并有硫酸存在的條件下,用過量的一定濃度的重鉻酸鉀溶液氧化土壤中的有機質(zhì)(碳),過剩的重鉻酸鉀溶液用標準硫酸亞鐵溶液滴定,進而依據(jù)消耗的重鉻酸鉀量計算有機碳的量。

    1.3 試驗樣本數(shù)據(jù)的建立

    Kostiakov三參數(shù)模型的表達式[10]為

    I=ktα+f0t(1)式中:I為土壤水分入滲時間為t時的累計入滲量,cm;t為土壤水分入滲時間,min;k為入滲系數(shù),cm/min;f0為穩(wěn)定入滲率,cm/min;α為入滲指數(shù)。

    利用入滲試驗得到入滲時間及其累計入滲量,進行非線性擬合,得到其相應的耕作土壤入滲參數(shù)值:人滲系數(shù)k、穩(wěn)定入滲率f0及入滲指數(shù)α。隨機選取3組耕作后土壤的Kostiakov三參數(shù)模型參數(shù)值,見表2。

    建立基本理化參數(shù)(見表1)與模型參數(shù)(見表2)間一一對應關系,構(gòu)成樣本序列,剔除個別有明顯錯誤的數(shù)據(jù),選取100組有代表性的試驗數(shù)據(jù)建立模型。

    2 BP神經(jīng)網(wǎng)絡模型的建立

    2.1 輸入、輸出參數(shù)的確定

    預報模型的輸入?yún)?shù)是根據(jù)耕作后土壤水分入滲特性的影響因素確定的。但若把全部影響因素作為輸人變量,必將使模型應用變得繁雜,因此在選擇模型的輸入?yún)?shù)時,僅考慮土壤水分入滲特性的主要影響因素[11]。經(jīng)過試驗分析,將土壤密度、土壤含水率、土壤質(zhì)地(砂粒、粉粒含量)、土壤有機質(zhì)含量作為輸入變量。由于水分通過地表進入土壤,地表對土壤水分入滲能力起控制作用[12],因此各影響因素均取0~20cm土層的平均值作為輸入?yún)?shù)。

    預報模型的輸出參數(shù)即模型預報量。反映土壤水分入滲能力大小的指標之一是土壤水分累計入滲量,為了確保土壤水分入滲達到穩(wěn)定入滲階段,選取90min土壤水分累計入滲量I90作為土壤水分入滲指標,即模型的輸出參數(shù);本文基于kostiakov經(jīng)驗模型公式建立BP模型,故選取入滲系數(shù)k、穩(wěn)定入滲率f0、入滲指數(shù)。作為預報模型的輸出變量。因此,本文所確定的模型輸出參數(shù)為90min土壤水分累計入滲量編、入滲系數(shù)k、穩(wěn)定入滲率f0、入滲指數(shù)α。

    2.2 BP神經(jīng)網(wǎng)絡的結(jié)構(gòu)與方法

    建立的BP神經(jīng)網(wǎng)絡模型由輸入層、隱含層、輸出層構(gòu)成。其中:輸入層有5個輸入?yún)?shù),即5個神經(jīng)元;輸出層有4個輸出參數(shù),即4個神經(jīng)元;而隱含層的節(jié)點數(shù)需要多次運算確定,通過逐漸增加隱含層節(jié)點數(shù),反復訓練樣本,達到模型精度時停止訓練[13],經(jīng)運算分析,本文建立的模型隱含層的神經(jīng)元數(shù)為18時可以達到模型精度,故建立的網(wǎng)絡拓撲結(jié)構(gòu)為5:18:4。

    網(wǎng)絡訓練函數(shù)為trainlm函數(shù),為加快參數(shù)的收斂,需對輸入、輸出樣本進行歸一化處理。Matlab7.0提供了歸一化函數(shù)premnmx以及還原函數(shù)post-mamx[14]。根據(jù)歸一化后數(shù)據(jù)的特點,建立的網(wǎng)絡模型選擇正切函數(shù)tansig函數(shù)作為隱含層的激活函數(shù),線性函數(shù)purelin函數(shù)作為輸出層的激活函數(shù)。模型訓練參數(shù)設為學習率0.01、最大學習迭代次數(shù)2000、訓練精度0.0001。BP模型是一種以誤差逆?zhèn)鞑ピ頌榛A的多層前饋神經(jīng)網(wǎng)絡,通過誤差反向傳播對輸入樣本反復訓練,不斷改變網(wǎng)絡鏈接權(quán)值,致使網(wǎng)絡輸出不斷接近期望值。

    2.3 入滲參數(shù)的預報模型

    2.3.1 模型結(jié)構(gòu)

    建立的BP網(wǎng)絡模型結(jié)構(gòu)為

    net=newff(min max(traininput),[18,4],

    {'tansig','purelin'},'trainlm')(2)式中:net為創(chuàng)建的BP神經(jīng)網(wǎng)絡;newff為在Matlab7.0中生成BP神經(jīng)網(wǎng)絡的函數(shù);min max()為輸入向量取值范圍的矩陣;[18,4]為隱含層和輸出層神經(jīng)元個數(shù);{'tansig',purelin'}為隱含層和輸出層的傳輸函數(shù);'trainlm'為網(wǎng)絡的訓練函數(shù)。

    BP神經(jīng)網(wǎng)絡的訓練結(jié)果為

    [α,k,f0,I90]=purelin(lw2×(tansig(iw1×p+b1))+b2)(3)式中:iw1為模型輸入層到隱含層的權(quán)值;lw2為模型隱含層到輸出層的權(quán)值;b1為模型輸入層到隱含層的閾值;b2為模型隱含層到輸出層的閾值;p=[δ,θ,ω,ψ,β],其中δ為耕層(0~20cm)土壤密度,θ為耕層(0~20cm)土壤質(zhì)量含水率,ω為耕層土壤砂粒含量,ψ為耕層土壤粉粒含量,β為耕層有機質(zhì)含量。

    2.3.2 預測結(jié)果分析

    基于Kostiakov三參數(shù)模型建立的BP神經(jīng)網(wǎng)絡對100組數(shù)據(jù)資料進行運算后,得到入滲指數(shù)α、入滲系數(shù)k、穩(wěn)定入滲率f0、90min累計入滲量I90的結(jié)果分別見表3、表4、表5、表6。

    預報模型對入滲指數(shù)α運算的結(jié)果顯示:其相對誤差最大值為4.3531%,最小值為0.0355%,平均值為0.9552%;其絕對誤差最大值為0.0108,最小值為0,平均值為0.0019。相對誤差均小于5%,模型預報精度較高。

    預報模型對入滲系數(shù)k運算結(jié)果顯示:其相對誤差最大值為5.8779%,最小值為0.0042%,平均值為1.0081%;其絕對誤差最大值為0.0757,最小值為0.0002,平均值為0.0190??梢娤鄬φ`差均小于6%,模型預報精度較高。

    預報模型對穩(wěn)滲率f0運算結(jié)果顯示:其相對誤差最大值為15.0228%,最小值為0.0032%,平均值為8.3113%;其絕對誤差最大值為0.0169,最小值為0.0004,平均值為0.0055??梢?,預報模型對穩(wěn)滲率f0預測的相對誤差大于人滲指數(shù)α與人滲系數(shù)k的,其預報精度較二者略低,但在可接受范圍內(nèi)。

    預報模型對90min累計入滲量I90運算結(jié)果顯示:其相對誤差最大值為9.3280%,最小值為0.0343%,平均值為2.9348%;其絕對誤差最大值為0.9695cm,最小值為0.0050cm,平均值為0.3155cm??梢?,其相對誤差均小于10%,預報精度較高。

    分析預報模型的訓練結(jié)果,說明耕作后的土壤基于Kostiakov三參數(shù)模型建立的BP預報模型對土壤水分入滲參數(shù)預測精度較高,預測結(jié)果可靠,可行性較強。

    3 結(jié)語

    本文建立的土壤水分入滲參數(shù)BP預報模型能夠客觀地反映耕作后土壤基本理化參數(shù)與Kostiakov三參數(shù)模型中的模型參數(shù)之間復雜的非線性關系,其預測結(jié)果有很高的精度。綜合分析認為,預報模型預測的入滲指數(shù)α、人滲系數(shù)k、穩(wěn)滲率f0及累計入滲量I90平均相對誤差分別為0. 9552%、1. 0081%、8.3113%、2.9348%,穩(wěn)滲率f0預測精度略低于其他參數(shù),但在可接受范圍內(nèi),分析認為:對于入滲指數(shù)α、入滲系數(shù)k,耕層(0~20cm)土壤理化參數(shù)對其影響很大,其關聯(lián)性較強,故預測精度較高;而對于穩(wěn)滲率f0,不僅與耕層(0~20cm)土壤理化參數(shù)有關聯(lián),而且與犁底層(20~40cm)土壤理化性質(zhì)也有一些關系,故其預測精度相對較低??傊?,本文所建立的BP預報模型整體預測精度較高,對耕作后土壤水分入滲參數(shù)的預測是可行的,這將為優(yōu)化耕作后農(nóng)田灌水技術參數(shù)提供理論與技術支撐。

    基于Kostiakov三參數(shù)經(jīng)驗模型公式建立的BP模型中,把耕層土壤密度、質(zhì)量含水率、土壤砂粒含量、粉粒含量及有機質(zhì)含量作為模型的輸入變量,預測其模型參數(shù)α、k、f0及I90,均可獲得比較可靠的結(jié)果。

    由于BP模型結(jié)構(gòu)的設計主要依賴于實踐經(jīng)驗,所需模型的樣本范圍及模型訓練頻數(shù)均在嘗試中,而且樣本也有一定的誤差,因此有待開發(fā)更為理想的方法對網(wǎng)絡模型結(jié)構(gòu)進行優(yōu)化,從而提高模型預測精度。

    參考文獻:

    [1] ABRISQUETA J M,PLANA V,MOUNZER 0 H,et al. Effectsof Soil Tillage on Runoff Generation in a Mediterranean Apricotor Chard[J].純ricultural Water Management,2007,93:11-18.

    [2]鄭子成,何淑勤,吳發(fā)啟.人為管理措施對產(chǎn)流產(chǎn)沙作用的研究[J].水土保持學報,2003,17(6):150-152.

    [3]雷廷武,毛麗麗,李鑫,等.土壤入滲性能的線源入流測量方法研究[J].農(nóng)業(yè)工程學報,2007,23(1):1-5.

    [4]曹崇文,樊貴盛.耕作土壤入滲能力衰減機理[J].太原理工大學學報,2007,38(2):116-118.

    [5]史曉楠,王全九,巨龍.微咸水人滲條件下Philip與Green-Ampt模型參數(shù)的對比分析[J].土壤學報,2007,44(2):360-363.

    [6]陳洪松,邵明安,王克林.土壤初始含水率對坡面降雨入滲及土壤水分再分布的影響[J].農(nóng)業(yè)工程學報,2006,22(1):44-47.

    [7]蘇正安,張建輝.耕作導致的土壤再分布對土壤水分入滲的影響[J].水土保持學報,2010,24(3):194-198.

    [8]李雪轉(zhuǎn),樊貴盛.土壤有機質(zhì)含量對土壤入滲能力及參數(shù)影響的試驗研究[J].農(nóng)業(yè)工程學報,2006,22(3):188-190.

    [9]李雪轉(zhuǎn),樊貴盛,王雪.影響土壤人滲積水時間因素的試驗[J].農(nóng)業(yè)工程學報,2009,25(9):27-30.

    [10]雷志棟,楊詩秀,謝森傳[M].北京:清華大學出版社,1988:35-41.

    [11]楊素宜,渾水入滲規(guī)律及畦灌補源技術研究[D].太原:太原理工大學,2006:11-23.

    [12]MAGNUS Porsson,RONNY Bemdtsson. Estimating TransferParameters in an Undisturbed Soil Column Using Time DomainReflectometry and Transfer Function Theory[J].Journal ofHydrology,1998,205:232-247.

    [13]伍春香,劉琳,王葆元.三層BP網(wǎng)隱層節(jié)點數(shù)確定方法的研究[J].武漢測繪科技大學學報,1999,24(2):85-87.

    [14]尹京川,馬孝義,孫永勝,等.基于BP神經(jīng)網(wǎng)絡與GIS可視化的作物需水量預測[J].中國農(nóng)村水利水電,2012(2):13-15.

    [15]周青云,孫西歡,康紹忠.蓄水坑灌條件下土壤水分運動的數(shù)值模擬[J].水利學報,2006,37(3):342-347.

    猜你喜歡
    耕作土壤水分含水率
    昆明森林可燃物燃燒機理研究
    基于表面能的濕煤顆粒含水率表征方法
    弱膨脹土增濕變形量試驗及路堤填筑分析
    西藏高原土壤水分遙感監(jiān)測方法研究
    原油含水率在線測量技術研究
    電子制作(2016年1期)2016-11-07 08:42:56
    耕作深度對紫色土坡地旋耕機耕作侵蝕的影響
    玉米保護性耕作的技術要領
    草地耕作技術在澳大利亞的應用
    土壤與作物(2015年3期)2015-12-08 00:46:58
    不同覆蓋措施對棗園土壤水分和溫度的影響
    西洞庭湖區(qū)免耕耕作模式及其配套技術
    作物研究(2014年6期)2014-03-01 03:39:12
    精品久久久精品久久久| 毛片一级片免费看久久久久| 欧美丝袜亚洲另类| 欧美精品亚洲一区二区| 中文字幕制服av| 国产在视频线精品| 亚洲国产色片| 99久久中文字幕三级久久日本| 亚洲国产欧美人成| 国产精品精品国产色婷婷| 色吧在线观看| 亚洲精品色激情综合| av在线播放精品| 女性被躁到高潮视频| 成年免费大片在线观看| 极品教师在线视频| 国产亚洲欧美精品永久| 免费人成在线观看视频色| 在线亚洲精品国产二区图片欧美 | 欧美性感艳星| 成人二区视频| 久久久久国产网址| 日韩精品有码人妻一区| 视频中文字幕在线观看| 亚洲欧美日韩另类电影网站 | 亚洲天堂av无毛| 亚洲av电影在线观看一区二区三区| 丝袜脚勾引网站| 又黄又爽又刺激的免费视频.| 成人亚洲精品一区在线观看 | 中文字幕精品免费在线观看视频 | 国产成人a∨麻豆精品| 噜噜噜噜噜久久久久久91| 成人国产av品久久久| 在线免费观看不下载黄p国产| 亚洲精品456在线播放app| 久热这里只有精品99| 99久久精品热视频| 一个人看视频在线观看www免费| 久久久久久久亚洲中文字幕| 欧美一级a爱片免费观看看| 成人综合一区亚洲| 国产免费一级a男人的天堂| 国产精品一及| 毛片女人毛片| 久久久精品免费免费高清| 我的老师免费观看完整版| 亚洲伊人久久精品综合| 18禁在线无遮挡免费观看视频| 极品教师在线视频| 夫妻性生交免费视频一级片| 丰满人妻一区二区三区视频av| 王馨瑶露胸无遮挡在线观看| 狠狠精品人妻久久久久久综合| 日韩在线高清观看一区二区三区| 成年美女黄网站色视频大全免费 | 午夜福利影视在线免费观看| 91在线精品国自产拍蜜月| 日韩亚洲欧美综合| 亚洲第一av免费看| 久久热精品热| 天堂8中文在线网| 美女高潮的动态| 蜜桃在线观看..| 国产乱人视频| 国产亚洲91精品色在线| 哪个播放器可以免费观看大片| 亚洲国产精品国产精品| 欧美精品人与动牲交sv欧美| 国产精品一区二区三区四区免费观看| 深夜a级毛片| 偷拍熟女少妇极品色| 久久精品久久久久久久性| 高清视频免费观看一区二区| 免费观看a级毛片全部| 欧美日韩视频高清一区二区三区二| 99热这里只有是精品在线观看| 中国美白少妇内射xxxbb| 这个男人来自地球电影免费观看 | 尾随美女入室| 亚洲天堂av无毛| 一区二区三区四区激情视频| 久久午夜福利片| 精品人妻熟女av久视频| 99热国产这里只有精品6| 一本色道久久久久久精品综合| 丰满迷人的少妇在线观看| 美女cb高潮喷水在线观看| 久久久久国产精品人妻一区二区| 亚洲av在线观看美女高潮| 赤兔流量卡办理| 精品酒店卫生间| 91狼人影院| 欧美激情极品国产一区二区三区 | 亚洲精华国产精华液的使用体验| 中文字幕人妻熟人妻熟丝袜美| a 毛片基地| 91精品国产九色| 国产淫片久久久久久久久| 欧美成人一区二区免费高清观看| 免费观看无遮挡的男女| 国产黄色视频一区二区在线观看| 日韩 亚洲 欧美在线| 中文字幕亚洲精品专区| 高清毛片免费看| 国产欧美日韩精品一区二区| 免费看光身美女| 久久热精品热| 亚洲国产精品专区欧美| 一本色道久久久久久精品综合| 亚洲av电影在线观看一区二区三区| 亚洲高清免费不卡视频| 美女高潮的动态| 秋霞伦理黄片| 人人妻人人澡人人爽人人夜夜| 欧美亚洲 丝袜 人妻 在线| 99热国产这里只有精品6| 大码成人一级视频| 18禁在线播放成人免费| av国产久精品久网站免费入址| 国产精品一区www在线观看| 熟妇人妻不卡中文字幕| 黄色欧美视频在线观看| 女人久久www免费人成看片| 性高湖久久久久久久久免费观看| 91午夜精品亚洲一区二区三区| 久久久午夜欧美精品| 五月伊人婷婷丁香| 日本午夜av视频| 色吧在线观看| 久久人人爽av亚洲精品天堂 | 欧美老熟妇乱子伦牲交| 寂寞人妻少妇视频99o| 丰满人妻一区二区三区视频av| 综合色丁香网| 国产视频内射| 美女中出高潮动态图| 九九在线视频观看精品| 亚洲美女视频黄频| 免费av中文字幕在线| 一区二区三区免费毛片| 亚洲成人一二三区av| 国产在线一区二区三区精| 成年av动漫网址| 日韩一区二区视频免费看| kizo精华| av.在线天堂| 亚洲美女视频黄频| 国产成人91sexporn| 久久青草综合色| 一级av片app| 免费看不卡的av| 你懂的网址亚洲精品在线观看| 国产永久视频网站| 亚洲熟女精品中文字幕| 日本猛色少妇xxxxx猛交久久| 国产av国产精品国产| 国产成人aa在线观看| 高清毛片免费看| 欧美老熟妇乱子伦牲交| 日韩欧美 国产精品| 日韩中字成人| 狂野欧美激情性xxxx在线观看| 国产免费视频播放在线视频| 美女视频免费永久观看网站| 狠狠精品人妻久久久久久综合| 成人亚洲精品一区在线观看 | 成人国产av品久久久| 国产久久久一区二区三区| av在线老鸭窝| 亚洲欧美成人精品一区二区| 久久精品国产亚洲av天美| 成人一区二区视频在线观看| 十八禁网站网址无遮挡 | 一级爰片在线观看| 精品酒店卫生间| 日韩欧美 国产精品| 亚洲欧美日韩无卡精品| 五月伊人婷婷丁香| 97超视频在线观看视频| 一区二区三区四区激情视频| 国内精品宾馆在线| av线在线观看网站| 尾随美女入室| 日韩中字成人| 亚洲成人手机| 国产免费一区二区三区四区乱码| 久久久久性生活片| 97精品久久久久久久久久精品| 99久久人妻综合| 亚洲最大成人中文| 永久网站在线| 久久久成人免费电影| 免费观看无遮挡的男女| 国产爽快片一区二区三区| 熟女av电影| 欧美精品人与动牲交sv欧美| 搡女人真爽免费视频火全软件| 最新中文字幕久久久久| 久久 成人 亚洲| a级毛色黄片| 99久久人妻综合| 久久精品国产亚洲av涩爱| 亚洲欧美一区二区三区黑人 | 国产成人午夜福利电影在线观看| 伊人久久精品亚洲午夜| 亚洲av日韩在线播放| 国产精品国产三级国产av玫瑰| 国产男人的电影天堂91| 一区二区三区四区激情视频| 国产精品久久久久成人av| 老师上课跳d突然被开到最大视频| 国产白丝娇喘喷水9色精品| 精品熟女少妇av免费看| 日韩精品有码人妻一区| 91久久精品国产一区二区三区| 嫩草影院入口| 色视频www国产| 久久久a久久爽久久v久久| 九草在线视频观看| 男女边摸边吃奶| 久久99蜜桃精品久久| 亚洲伊人久久精品综合| 国产精品女同一区二区软件| 妹子高潮喷水视频| 国产成人精品一,二区| 大香蕉久久网| 亚洲欧美精品自产自拍| 免费人妻精品一区二区三区视频| 又粗又硬又长又爽又黄的视频| 精品久久久久久久末码| 女的被弄到高潮叫床怎么办| 成人一区二区视频在线观看| 熟女电影av网| 欧美日韩视频精品一区| 男女国产视频网站| 麻豆成人午夜福利视频| 亚洲av成人精品一区久久| 国产在线免费精品| 国产女主播在线喷水免费视频网站| 免费观看a级毛片全部| 啦啦啦在线观看免费高清www| 一级二级三级毛片免费看| 18禁裸乳无遮挡免费网站照片| 国产视频内射| 观看免费一级毛片| 国产久久久一区二区三区| 一级爰片在线观看| h日本视频在线播放| 男女边摸边吃奶| 九色成人免费人妻av| 国产精品国产三级国产专区5o| 亚洲丝袜综合中文字幕| 亚洲国产高清在线一区二区三| 久久久国产一区二区| 男人和女人高潮做爰伦理| av在线老鸭窝| 久久久久久久久久久免费av| 寂寞人妻少妇视频99o| 大片电影免费在线观看免费| 春色校园在线视频观看| 直男gayav资源| 久久久久久久久久久免费av| 97超碰精品成人国产| 在线观看av片永久免费下载| 亚洲国产精品一区三区| 最近的中文字幕免费完整| 狂野欧美激情性bbbbbb| 国产乱人偷精品视频| 99久久精品热视频| 各种免费的搞黄视频| 久久久久精品性色| 久久精品国产a三级三级三级| 九九在线视频观看精品| 亚洲国产精品专区欧美| 国产亚洲精品久久久com| 国产日韩欧美在线精品| 夜夜看夜夜爽夜夜摸| 多毛熟女@视频| 亚洲欧美日韩另类电影网站 | 久久久久视频综合| 久久久久久久久久久丰满| 一区二区三区精品91| 国产精品久久久久久精品电影小说 | 亚洲国产精品999| 午夜激情久久久久久久| 又爽又黄a免费视频| 观看av在线不卡| 国产精品一区www在线观看| 熟女电影av网| 插阴视频在线观看视频| 成人亚洲精品一区在线观看 | 少妇的逼水好多| 久久久久久久久大av| 欧美成人精品欧美一级黄| 日韩人妻高清精品专区| 在现免费观看毛片| 肉色欧美久久久久久久蜜桃| 美女内射精品一级片tv| 欧美激情极品国产一区二区三区 | 波野结衣二区三区在线| 女人十人毛片免费观看3o分钟| 亚洲精品一二三| 精品午夜福利在线看| 九九爱精品视频在线观看| 一区二区三区精品91| 熟女人妻精品中文字幕| 国产有黄有色有爽视频| 一本久久精品| 夫妻性生交免费视频一级片| 王馨瑶露胸无遮挡在线观看| 国产精品一区二区在线观看99| 欧美性感艳星| 久久久久久久亚洲中文字幕| 亚洲精品视频女| 在线观看免费高清a一片| 新久久久久国产一级毛片| 在线看a的网站| 一本色道久久久久久精品综合| 亚洲欧美中文字幕日韩二区| 久久精品国产亚洲av涩爱| 丰满人妻一区二区三区视频av| 国产69精品久久久久777片| 国产成人freesex在线| 久久99精品国语久久久| 一区二区三区四区激情视频| 国产在线一区二区三区精| 亚洲成人av在线免费| 九色成人免费人妻av| 超碰97精品在线观看| 日韩欧美 国产精品| 免费观看的影片在线观看| 亚洲色图av天堂| av在线播放精品| 99re6热这里在线精品视频| 成人特级av手机在线观看| 国产免费又黄又爽又色| 日韩三级伦理在线观看| 最近最新中文字幕大全电影3| www.色视频.com| 久久久久久久大尺度免费视频| 国产亚洲5aaaaa淫片| 狂野欧美白嫩少妇大欣赏| 亚洲精品aⅴ在线观看| 99久国产av精品国产电影| 久久国产精品大桥未久av | 晚上一个人看的免费电影| 国产精品国产三级专区第一集| 高清欧美精品videossex| 国产男女超爽视频在线观看| 18禁在线无遮挡免费观看视频| 少妇高潮的动态图| 视频中文字幕在线观看| 成人高潮视频无遮挡免费网站| 精品久久久久久久久av| 五月天丁香电影| 丝袜喷水一区| 九九爱精品视频在线观看| 亚洲av二区三区四区| 日本猛色少妇xxxxx猛交久久| 91精品国产九色| 人人妻人人澡人人爽人人夜夜| 成人漫画全彩无遮挡| 直男gayav资源| 精品国产三级普通话版| av在线app专区| 一级毛片 在线播放| 涩涩av久久男人的天堂| www.色视频.com| 国产乱人偷精品视频| 色婷婷av一区二区三区视频| 国产一区有黄有色的免费视频| 日韩三级伦理在线观看| h视频一区二区三区| 国产一区亚洲一区在线观看| 成人美女网站在线观看视频| 有码 亚洲区| 久久久精品免费免费高清| 欧美另类一区| 亚洲精品一二三| 18禁在线播放成人免费| 欧美激情极品国产一区二区三区 | 麻豆成人av视频| 人人妻人人添人人爽欧美一区卜 | 欧美国产精品一级二级三级 | 欧美亚洲 丝袜 人妻 在线| 国产视频内射| 欧美bdsm另类| 国产91av在线免费观看| 七月丁香在线播放| 超碰av人人做人人爽久久| h视频一区二区三区| 久久久国产一区二区| 99久久精品一区二区三区| 精品久久久久久久久av| 欧美老熟妇乱子伦牲交| 亚洲,一卡二卡三卡| 久久久久久久久久久免费av| 中文字幕精品免费在线观看视频 | 极品教师在线视频| 菩萨蛮人人尽说江南好唐韦庄| 精品亚洲成a人片在线观看 | 国产一区二区三区综合在线观看 | 欧美日韩一区二区视频在线观看视频在线| 久久久久网色| 日产精品乱码卡一卡2卡三| 97精品久久久久久久久久精品| 只有这里有精品99| 蜜桃亚洲精品一区二区三区| 亚洲欧美成人综合另类久久久| 高清在线视频一区二区三区| 一级片'在线观看视频| 午夜视频国产福利| 校园人妻丝袜中文字幕| 黄色一级大片看看| 国语对白做爰xxxⅹ性视频网站| 五月开心婷婷网| 亚洲国产最新在线播放| 亚洲欧美一区二区三区黑人 | 一级爰片在线观看| 在线观看三级黄色| 亚洲精华国产精华液的使用体验| 有码 亚洲区| 欧美高清性xxxxhd video| 欧美xxⅹ黑人| 亚洲精品色激情综合| 一级毛片电影观看| 一边亲一边摸免费视频| 97精品久久久久久久久久精品| 人妻一区二区av| 久久久久国产网址| 美女内射精品一级片tv| 国产在线一区二区三区精| 精华霜和精华液先用哪个| 99热全是精品| 日本免费在线观看一区| 国产高清国产精品国产三级 | 日本与韩国留学比较| 欧美精品亚洲一区二区| 国产女主播在线喷水免费视频网站| 99久国产av精品国产电影| 99视频精品全部免费 在线| 日本猛色少妇xxxxx猛交久久| 大片免费播放器 马上看| 少妇的逼水好多| 精品久久久久久久久av| 国产亚洲av片在线观看秒播厂| 99久国产av精品国产电影| 国产69精品久久久久777片| 久久久亚洲精品成人影院| 大香蕉97超碰在线| 成人亚洲精品一区在线观看 | 观看美女的网站| 午夜免费男女啪啪视频观看| 亚洲内射少妇av| 精品亚洲乱码少妇综合久久| av在线观看视频网站免费| 最近中文字幕高清免费大全6| 欧美人与善性xxx| 中国国产av一级| 精品一区二区三卡| 成人黄色视频免费在线看| 毛片女人毛片| 黄色日韩在线| 超碰av人人做人人爽久久| 在线观看免费视频网站a站| 亚洲av.av天堂| 国产精品一二三区在线看| 亚洲第一区二区三区不卡| 国模一区二区三区四区视频| 又黄又爽又刺激的免费视频.| 日本av免费视频播放| 高清欧美精品videossex| 伊人久久精品亚洲午夜| 久久久精品94久久精品| 免费人成在线观看视频色| 成人国产av品久久久| 日本一二三区视频观看| 精品国产一区二区三区久久久樱花 | 久久综合国产亚洲精品| 欧美日本视频| 夫妻午夜视频| 午夜免费观看性视频| 自拍欧美九色日韩亚洲蝌蚪91 | av视频免费观看在线观看| 丰满人妻一区二区三区视频av| 老熟女久久久| 少妇人妻 视频| 妹子高潮喷水视频| 免费人妻精品一区二区三区视频| 国产黄色视频一区二区在线观看| 永久免费av网站大全| 99久久精品国产国产毛片| 岛国毛片在线播放| 精品酒店卫生间| 成人漫画全彩无遮挡| 亚洲国产精品999| 精品一区在线观看国产| 小蜜桃在线观看免费完整版高清| av免费在线看不卡| 高清毛片免费看| 久久精品国产自在天天线| 日韩 亚洲 欧美在线| 性高湖久久久久久久久免费观看| xxx大片免费视频| 一级黄片播放器| 亚洲成人av在线免费| 91久久精品国产一区二区三区| 国产乱人偷精品视频| 国产极品天堂在线| 18禁裸乳无遮挡免费网站照片| 中国美白少妇内射xxxbb| 精品一区二区免费观看| 国产中年淑女户外野战色| 亚洲中文av在线| 黑人猛操日本美女一级片| 欧美丝袜亚洲另类| 国产精品久久久久久久电影| 观看美女的网站| 成人综合一区亚洲| av国产免费在线观看| 在线观看免费高清a一片| 午夜免费观看性视频| av免费观看日本| 在线免费十八禁| 大陆偷拍与自拍| 偷拍熟女少妇极品色| 精品亚洲成国产av| 自拍偷自拍亚洲精品老妇| 欧美日韩国产mv在线观看视频 | 亚洲欧洲国产日韩| 亚洲国产成人一精品久久久| 在线观看一区二区三区| 97在线人人人人妻| 一本—道久久a久久精品蜜桃钙片| 欧美97在线视频| 观看美女的网站| 亚洲欧美一区二区三区黑人 | 国产白丝娇喘喷水9色精品| 乱系列少妇在线播放| 欧美日韩视频精品一区| 男人添女人高潮全过程视频| 校园人妻丝袜中文字幕| 两个人的视频大全免费| 如何舔出高潮| 精品久久国产蜜桃| 亚洲av二区三区四区| 亚洲精品第二区| 国产成人免费观看mmmm| 麻豆成人av视频| 男人和女人高潮做爰伦理| 美女主播在线视频| 大又大粗又爽又黄少妇毛片口| 少妇裸体淫交视频免费看高清| 亚洲欧美精品自产自拍| 18禁动态无遮挡网站| 亚洲欧洲日产国产| 亚洲色图av天堂| 国产中年淑女户外野战色| 国产成人免费观看mmmm| 99re6热这里在线精品视频| 亚洲精品日本国产第一区| 少妇丰满av| 一级毛片aaaaaa免费看小| 青春草亚洲视频在线观看| 肉色欧美久久久久久久蜜桃| 91精品国产九色| 麻豆国产97在线/欧美| 国产精品一区二区在线观看99| 少妇熟女欧美另类| 久久久久视频综合| 岛国毛片在线播放| 少妇人妻精品综合一区二区| 我的女老师完整版在线观看| 日日摸夜夜添夜夜添av毛片| 欧美成人一区二区免费高清观看| 亚洲美女黄色视频免费看| 亚洲三级黄色毛片| 青青草视频在线视频观看| 看免费成人av毛片| 精品久久久久久久末码| 亚洲av综合色区一区| av福利片在线观看| 毛片女人毛片| freevideosex欧美| 极品教师在线视频| 中国三级夫妇交换| 人妻 亚洲 视频| 这个男人来自地球电影免费观看 | 国产91av在线免费观看| 久久久久久伊人网av| 美女中出高潮动态图| 国国产精品蜜臀av免费| 精品少妇黑人巨大在线播放| 国模一区二区三区四区视频| 交换朋友夫妻互换小说| 性色av一级| 国产一区二区在线观看日韩| 国产成人午夜福利电影在线观看| 97热精品久久久久久| 亚洲熟女精品中文字幕| 80岁老熟妇乱子伦牲交| 97热精品久久久久久| 日本色播在线视频| 联通29元200g的流量卡| 国产精品久久久久久av不卡| 日日啪夜夜撸| 一个人免费看片子| 在线观看免费高清a一片| 久久99蜜桃精品久久| 在线观看免费视频网站a站| 性色avwww在线观看| 国产 一区精品| 另类亚洲欧美激情| 午夜福利影视在线免费观看|