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

    新疆干旱區(qū)膜下滴灌土壤水平衡模擬

    2017-03-21 01:20:42陸寶宏王惠鳳程增輝陸建宇
    節(jié)水灌溉 2017年3期
    關鍵詞:區(qū)域模型

    胡 越,陸寶宏,2,熊 絲, 王惠鳳, 程增輝, 陸建宇

    (1.河海大學水文水資源學院,南京 210098;2.河海大學水文水資源與水利工程科學國家重點實驗室,南京 210098;3.中水北方勘測設計研究有限責任公司,天津 300222)

    0 引 言

    我國是農業(yè)大國,農業(yè)用水占到全國用水總量60%以上[1],尤其新疆地區(qū)水資源更加貧乏,其水資源短缺問題已嚴重制約了區(qū)域農業(yè)的發(fā)展。因此,研究區(qū)域農業(yè)灌溉技術對保障區(qū)域農業(yè)生產安全,推動新疆地區(qū)農業(yè)生產進步及經濟持續(xù)繁榮具有重要意義。

    近年來,新疆生產建設兵團一直不斷探索節(jié)水灌溉種植方式,其中膜下滴灌將薄膜覆蓋技術與滴灌技術相結合,能夠在保證作物正常生長發(fā)育同時,最大限度節(jié)約水資源,是多種灌溉技術中最為有效的一種新型農田灌溉技術[2]。迄今為止,在模擬不同灌溉條件下的土壤水分運動方面,眾多專家學者做了大量針對性研究并取得一系列成果。張建軍[3]較好地模擬了滴灌條件下的土壤水運動狀況,并給出了不同灌水量和滴頭流量對土壤含水率的影響;劉群昌[4]等模擬了隔畦灌溉在不同深度的土壤水運動情勢;李亮[5]研究了內蒙古河套灌區(qū)荒地的土壤水分模擬問題;VB Bufon[6]等模擬并比較了德克薩斯高平原地下滴灌種植棉花灌溉前后土壤水變化;劉梅先[7]等模擬了新疆地區(qū)膜下滴灌膜中、管下和膜間等位置處土水情勢,并運用實測數據反演模型參數;李坡[8]等針對不同土壤類型,分別作了適宜數值模擬,并推求了不同條件下需水量。上述研究成果在一定程度上推動了我國農業(yè)灌溉技術的發(fā)展。但是,關于土壤水平衡問題的研究仍然較少,尤其缺乏干旱區(qū)以1膜4行覆膜方式棉花種植研究。目前,1膜4行覆膜方式具有提高地膜覆蓋率、增加種植密度等優(yōu)勢,十分有利于實現(xiàn)作物的保產、穩(wěn)產和增產,可在全國較大范圍內推廣應用。

    鑒于此,本文以石河子試驗站覆膜棉花種植觀測資料為基礎,運用HYDRUS-2D模型模擬棉花在生育期土壤水運動狀況,并分析了試驗區(qū)土壤水平衡情勢,以期研究干旱灌區(qū)的土壤水平衡問題,為掌握干旱灌區(qū)土壤水平衡,促進區(qū)域農業(yè)生產提供一定的參考。

    1 研究區(qū)概況

    瑪納斯河流域位于新疆天山北麓中段,地處準噶爾盆地南緣,歐亞大陸腹地,遠離海洋,屬于典型的大陸性氣候干旱區(qū)。該流域介于東經85°01′~86°32′,北緯43°27′~45°21′之間,流域面積19 800 km2,其中平原面積14 650 km2,是北疆重要的經濟核心帶[9]。

    本文即以瑪納斯河流域山前平原區(qū)石河子試驗站為研究區(qū)域。試驗地點位于新疆石河子大學節(jié)水灌溉試驗站,屬瑪納斯河流域洪積沖積扇中部,年平均氣溫年平均7.4 ℃,日照時數為2 705 h,多年平均降雨量207 mm,年蒸發(fā)量1 506 mm(小型蒸發(fā)皿),無霜期平均160 d,灌溉用水取自30 m淺井水,井水平均礦化度0.87 g/L,潛水位常年在7~11 m之間變化。

    2 材料來源與方法

    2.1 田間試驗

    本試驗區(qū)為試驗站B5小區(qū),面積54 m2(5.4 m×10 m)。試驗時間為2009年6-9月,即棉花的生育中期(6-8月)和生育末期(9月)。試驗區(qū)原始土壤類型為荒漠灰鈣土,經過多年耕作改良熟化為耕作灰漠土,土壤質地為砂壤土, 肥力中等。土壤屬性見表1。

    表1 試驗區(qū)土壤物理特性Tab.1 Soil physical characteristics in the test area

    膜下滴灌采用1膜4行栽培種植模式(見圖1):即1膜內鋪設2條毛管,1 條毛管灌溉2行棉花,行距配置為30 cm +60 cm + 30 cm,即作物間距30 cm和60 cm,毛管間距90 cm,滴頭間距30 cm,膜間距40 cm,膜寬140 cm?;跀祿?,本次模擬從2009年7月27日至9月3日共計39 d。在模擬期內,從7月27日起,每隔6~8 d灌水1次,共灌水5次,滴頭流量1~2 L/h。

    圖1 1膜4行布置圖(單位:cm)Fig.1 The layout of one mulch with four rows

    觀測內容包括試驗區(qū)水文氣象要素,如降水、蒸發(fā)、氣溫和日照等;以及土壤含水率。由于1膜2管4行種植具有高度對稱性,故從地膜中心位置處一分為二直到兩個地膜中間處截止為一個單元,水平向即為90 cm。土壤含水率觀測在豎直向上,地面以下5 cm為第一層次,地面以下10 cm是第二層次,往下每間隔10 cm觀測一次,直到距土壤表面100 cm處截止,觀測點分別布設在膜中、膜內、管下、膜外和膜間5個位置處(如圖2),共計46個觀測點,觀測頻度是每10 d一次,降雨后與灌溉前后加測。

    圖2 土壤含水率觀測點布置Fig.2 The disposal of observation points of soil water content

    2.2 數學模型

    HYDRUS-2D軟件由計算機程序和交互式圖形界面組成,基于有限元方法求解修改過的Richards方程來模擬二維飽和-非飽和土壤水運動問題。本次模擬運用HYDRUS-2D模型,建立新疆石河子地區(qū)基于1膜4行覆膜方式的棉花土壤水運動模型,采用實測土壤含水率數據對模型進行參數率定和模型檢驗,最后用檢驗好的模型分析膜下滴灌土壤水平衡情況。

    2.2.1 模型構建

    1膜4行下的土壤水運動是三維問題,當假設土壤均質各向同性,且由于覆膜方式又具有對稱特性時,可模擬一個覆膜單元膜中間分開右側(左側),膜下滴灌問題即可簡化為對稱線源二維問題。在模擬期內,根據實測數據依次輸入日降水量、日蒸散發(fā)量以及灌溉水量,模擬區(qū)域水平方向90 cm,豎直方向100 cm,生成三角形有限元網格,形成節(jié)點1 878個,為考慮水分運動和根系吸水兩個模擬過程,將土壤分層設置,其中0~40 cm為耕作層,40~100 cm為基層,考慮到棉花根系位于豎直方向20~60 cm范圍內,水平向滴管處兩側各15 cm范圍內,最低和最高時間步長乘法因子取1.3和0.7,同時,以含水率方式指定初始條件,并依據實際土壤含水率觀測點設置模擬觀測點,以便提高模型檢驗效率。

    2.2.2 基本方程

    (1)土壤水運動方程?;诰€性的Galerkin有限元方法,模擬非飽和水流求解Richards方程:

    (1)

    式中:θ為土壤體積含水率,cm3/cm3;Kh為水力傳導度,cm/d;t為時間,d;h為土壤負壓水頭,cm;x代表水平方向,cm;z代表豎直方向,cm;S為源匯項,d-1。

    (2)根系吸水?;贔eddes模型,求解源匯項S:

    S(h)=α(h)β(x,z)Tp

    (2)

    式中:α(h)為水分脅迫影響函數;β(x,z)為根系吸水密度分布;Tp為潛在蒸騰速率,cm/d。

    (3)土壤水分特征方程?;趘an Genuchten方程,求解土壤水分特征曲線θ(h)以及非飽和水力傳導度函數K(h):

    (3)

    θ(h)=θs(h≥0)

    (4)

    K(h)=KsS0.5e[1-(1-S1/me)m]2(h≤0)

    (5)

    (6)

    式中:θs為土壤飽和含水率,cm3/cm3;θr為土壤殘余含水率,cm3/cm3;Ks為飽和土壤導水率,cm/d;Se為土壤相對飽和度,cm3/cm3;α、n與m(m=1-1/n)是形狀參數,θ是土壤含水率,cm3/cm3。

    2.2.3 初始條件和邊界條件

    取2009年7月26日實測土壤含水率為初始土壤含水率,模擬區(qū)域初始含水率,見圖3。由模型邊界條件(見圖4)可知,由于模擬區(qū)域上邊界水平方向0~70 cm覆膜,故采用不透水邊界,但在滴管處考慮入滲半徑,此處采用變流量邊界,水平方向70~90 cm膜間裸露在大氣下,采用大氣邊界;不考慮地下水對模擬區(qū)影響,下邊界取自由排水邊界;模擬區(qū)左右對稱,因此左右邊界取不透水邊界。

    圖3 初始含水率Fig.3 Initial water content

    圖4 邊界條件Fig.4 Boundary condition

    2.2.4 模型參數

    (1)土壤水力參數。模擬區(qū)豎直方向(0~100 cm)土壤分為耕作層(0~40 cm)與基層(40~100 cm),相應土壤水力參數也要分層設置。以2009年7月27日至8月15日土壤含水率模擬值與實測值進行擬合,率定后利用8月20日至9月3日的模擬值與實測值進行驗證,經均方根誤差(RMSE)評價,率定后的土壤水力參數滿足試驗精度,具體數值見表2。

    表2 土壤水力參數Tab.2 Soil hydraulic parameters

    (2)根系吸水參數。 結合本次模擬實際情況,并查閱相關資料[10],確定根系吸水的7個參數值,見表3。

    表3 根系吸水參數Tab.3 Root water uptake parameters

    3 結果分析與討論

    3.1 不同位置處土壤含水率模擬

    基于HYDRUS-2D模型,得到模擬期內不同位置處土壤含水率值,并將其與實測土壤含水率進行比較,見圖5。同時,以均方根誤差(RMSE)和平均相對誤差(RE)[11],分析土壤含水率模擬精度,見表4。

    圖5 不同位置處含水率比較Fig.5 The compassion of water content in different positions

    評價指標膜中膜內管下膜外膜間RMSE(cm3·cm-3)0.0430.0580.0660.0610.051RE/%10.218.517.718.713.9

    由表4可知,5個位置處RMSE均小于0.1,且RE在20% 以內,說明誤差在可接受精度范圍內,整體模擬效果較好。但膜外處RE相對較大,可能因此處上邊界是大氣邊界,易受降水、蒸發(fā)以及左側灌溉水等因素綜合影響,因而土壤含水率變化較大,再結合圖5(d),土壤含水率模擬值略小于實測值,這可能與模擬過程中蒸散發(fā)量分配有關;管下處RMSE相對較大,可能因此處上邊界條件是變流量邊界,當有灌溉時,此處土壤含水率變化最大,結合圖5(c),管下處的實測含水率相對較大可能與根系吸水參數以及土壤水力參數設置以及實測資料精確度有關;膜中處,RMSE與RE均較小,可能是此處離滴管以及大氣邊界距離較遠,土壤含水率所受干擾較少。總體上圖5中的實測值與模擬值繪出的點與1:1線十分接近,表明HYDRUS-2D可較好模擬干旱區(qū)土壤水運動狀況。

    3.2 不同深度處土壤含水率模擬

    模擬期內土層不同深度含水率實測值與模擬值擬合情況以及隨時間的變化規(guī)律,見圖6。從圖6可看出:①不同深度處(-5、-10、-20、-50、-80、-100 cm)土壤含水率模擬值與實測值偏差較小,整體模擬效果較好;②-5、-10、-20 cm深度處模擬效果最佳, -50、-80、-100 cm處模擬效果相對一般,可能與模型中土壤水力參數取值,尤其是飽和滲透系數(Ks)以及土壤的空間變異性有關;③模擬期內,淺層土(-0~-50 cm)土壤含水率在0.25~0.30 cm3/cm3波動幅度較大,較深層土(-50~-100 cm)土壤含水率維持0.30 cm3/cm3上下,這可能因淺層土易受灌溉水、地表蒸發(fā)以及降雨等影響。較深層土距滴管以及裸地表面較遠,土壤含水率變化趨勢相對平緩,再通過圖6(a)、(b)、(c)易知,每次灌水土壤含水率均先升后降,第10~20 d灌溉2次,時間間隔較小,故淺層土壤含水率相對較高,且沒有明顯降低趨勢。

    3.3 土壤水平衡分析

    根據新疆石河子棉田覆膜灌溉實際情況,不難得出:模擬區(qū)域入水項為從大氣邊界流入的降雨量以及通過滴管進入土壤的灌溉量;模擬區(qū)域出水項為根系吸水量,大氣邊界土壤的無效蒸發(fā)量以及自由下滲量。

    圖6 不同深度處含水率比較Fig.6 The compassion of water content in different depths

    本次忽略膜下土壤的蒸發(fā)作用,同時不考慮區(qū)域地表徑流。HYDRUS-2D模型模擬區(qū)域水平衡公式如下:

    ΔW=W2-W1=P+I-E-T-Rc-Rs-Rg

    (7)

    式中:ΔW為模擬區(qū)域水量的變化量,mm;W1和W2分別為初始時和最終時土壤含水量,mm;P為降雨量,mm;I為灌溉量,mm;E、T為蒸發(fā)和蒸騰量,mm;Rc為植株水分儲變量,mm;Rs為地表徑流,mm;Rg為土壤水與地下水交換量,mm。

    結合本試驗模擬情景,可將式(7)簡化為:

    ΔW=W2-W1=P+I-E-R-Rg

    (8)

    式中:R代表根系吸水量,mm;其余字母含義與式(7)同。

    新疆石河子模擬區(qū)模擬期內入水量及出水量累計變化情況,見圖7和圖8。由圖7、圖8可知:①模擬期內,降雨量值很小,僅在灌溉當天區(qū)域累計入水量會有明顯增加,整體呈階梯性增長變化趨勢;②根系吸水累計量占出水累計量的絕大部分,并遠大于下滲量以及無效蒸發(fā)量,而下滲量和無效蒸發(fā)量相對較小,且下滲量大體呈線性增加而無效蒸發(fā)量在29 d后比29 d前增加的要快,這是由于在29 d后棉田進入生育末期,其作物的蒸騰量有所降低,相應的無效蒸發(fā)量會比之前生育中期大,區(qū)域出水累計量大體上呈現(xiàn)線性增長。

    新疆石河子模擬區(qū)域模擬期內的各出水項速率變化情況,見圖9。由圖9可知,根系吸水速率在前25 d穩(wěn)步下降,29 d后下降趨勢顯著增加;蒸發(fā)速率在前29 d穩(wěn)定在0.33 mm/d,29 d后在1 mm/d左右,有起伏是受大氣邊界降雨的影響;而下滲速率穩(wěn)定在0.5 mm/d左右,變化幅度不明顯。

    圖7 模擬區(qū)域入水量累計Fig.7 Accumulated water inflow in the simulation

    圖8 模擬區(qū)域出水量累計Fig.8 Accumulated water outflow in the simulation region

    圖9 模擬區(qū)域出水各項速率Fig.9 The outflow rates of each items in the simulation region

    基于HYDRUS-2D模型,給出模擬區(qū)域39 d出入水量結果,見表5。由表5可知,模擬期內,區(qū)域入水量為289.23 mm,其中灌溉量占到入水量的97.6%,而降雨量只有總入水量的2.34%,可見新疆瑪納斯河流域種植棉花所消耗大部分水是由灌溉來水提供的,而降雨補給相對較少;出水量為272.8 mm,其中根系吸水量占總出水量的85.7%,無效蒸發(fā)與下滲各占總出水量的7.28%和6.99%,可見1膜4行覆膜方式可使土壤無效蒸發(fā)大大減少,下滲量也有一定程度減少,80%以上的水量均由作物吸收,這大幅提高了區(qū)域水資源利用效率;區(qū)域總入水量比出水量多16.43 mm,可見作物生育中、末期以1膜4行覆膜方式下的土壤水量在得到灌水補給后不會降低,甚至會略有升高。

    表5 模擬區(qū)域39 d累計出入水量表 mm

    4 結 語

    基于HYDRUS-2D模型模擬了新疆瑪納斯河流域石河子試驗站種植棉田土壤水運動狀況,得到以下結論。

    (1)模擬期內, 5個位置處(膜中、膜內、管下、膜外、膜間)RMSE均小于0.1,RE在20% 以內,不同深度處(-5、-10、-20、-50、-80、-100 cm)模擬效果較好,說明HYDRUS-2D模型可用來模擬新疆干旱區(qū)棉田1膜4行覆膜方式的土壤水運動狀況。

    (2)模擬區(qū)域水量基本保持平衡狀態(tài),其中灌溉水量是區(qū)域入水的主要部分,而作物生長發(fā)育根系吸水占到了模擬區(qū)域耗水量80%以上,滲漏和裸地無效蒸發(fā)均不到耗水總量的10%。因此,1膜4行覆膜方式可以有效利用灌溉水,也將減少地下深層滲漏以及裸地無效蒸發(fā)。

    (3)因實測資料有限,未涉及新疆棉田全生育期土壤水平衡模擬研究。未來在擴充研究資料、研究對象的基礎上,進一步探討和論證HYDRUS-2D模型在新疆地區(qū)土壤水平衡情勢分析中的適用性。

    致謝:感謝楊玫老師提供的數據資料,使本文順利完成。

    [1] 佚 名. 節(jié)水灌溉篇[J]. 中國農村水利水電,2015,(12).

    [2] 康 靜,黃興法. 膜下滴灌的研究及發(fā)展[J]. 節(jié)水灌溉, 2013,(9):71-74.

    [3] 張建君. 滴灌施肥灌溉土壤水氮分布規(guī)律的試驗研究及數學模擬[D]. 北京:中國農業(yè)科學研究院, 2002.

    [4] 劉群昌,王韶華,蘇軼醒. 隔畦灌溉適宜的土壤及灌水量研究[J]. 節(jié)水灌溉, 2007,(6):39-42.

    [5] 李 亮,史海濱,賈錦鳳,等. 內蒙古河套灌區(qū)荒地水鹽運移規(guī)律模擬[J]. 農業(yè)工程學報, 2010,26(1):31-35.

    [6] Bufon V B, Lascano R J, Bednarz C, et al. Soil water content on drip irrigated cotton: Comparison of measured and simulated values obtained with the Hydrus 2-D model[J]. Irrigation Science, 2012,30(4):259-273.

    [7] Liu M X, Yang J S, Xiao-Ming L I, et al. Numerical simulation of soil water dynamics in a drip irrigated cotton field under plastic mulch[J]. Pedosphere, 2013,23(5):620-635.

    [8] 李 坡, 吳飛青, 胡正峰,等. 節(jié)水灌溉中土壤水運動模擬研究[J]. 節(jié)水灌溉, 2015,(3):46-51.

    [9] 王振華. 典型綠洲區(qū)長期膜下滴灌棉田土壤鹽分運移規(guī)律與灌溉調控研究[D].北京:中國農業(yè)大學, 2014.

    [10] 吳爭光. 棉花膜下滴灌土壤水鹽運移規(guī)律研究及數值模擬[D]. 烏魯木齊:新疆農業(yè)大學, 2011.

    [11] 郝遠遠,徐 旭,任東陽,等. 河套灌區(qū)土壤水鹽和作物生長的HYDRUS-EPIC模型分布式模擬[J]. 農業(yè)工程學報, 2015,(11):110-116,315.

    猜你喜歡
    區(qū)域模型
    一半模型
    永久基本農田集中區(qū)域“禁廢”
    分割區(qū)域
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權M-估計的漸近分布
    3D打印中的模型分割與打包
    關于四色猜想
    分區(qū)域
    FLUKA幾何模型到CAD幾何模型轉換方法初步研究
    基于嚴重區(qū)域的多PCC點暫降頻次估計
    電測與儀表(2015年5期)2015-04-09 11:30:52
    精品一品国产午夜福利视频| 亚洲熟女精品中文字幕| 亚洲欧洲日产国产| av女优亚洲男人天堂| 中文字幕免费在线视频6| 成人午夜精彩视频在线观看| 亚洲电影在线观看av| 久久精品久久久久久噜噜老黄| 欧美少妇被猛烈插入视频| 亚洲精品国产av成人精品| 亚洲人成网站在线播| 高清视频免费观看一区二区| 男人舔奶头视频| 少妇的逼好多水| 国产亚洲精品久久久com| 亚洲精品乱码久久久v下载方式| 亚洲久久久国产精品| 国内少妇人妻偷人精品xxx网站| 久久久成人免费电影| 啦啦啦中文免费视频观看日本| 国产av一区二区精品久久 | 亚洲自偷自拍三级| 欧美xxxx黑人xx丫x性爽| 国产91av在线免费观看| 成人亚洲精品一区在线观看 | 99热6这里只有精品| 成人综合一区亚洲| 人妻夜夜爽99麻豆av| 亚州av有码| 能在线免费看毛片的网站| 久久韩国三级中文字幕| 亚洲欧洲国产日韩| 欧美成人一区二区免费高清观看| 直男gayav资源| 少妇高潮的动态图| av免费观看日本| 99九九线精品视频在线观看视频| 蜜桃亚洲精品一区二区三区| 伦理电影大哥的女人| 欧美日韩精品成人综合77777| 青青草视频在线视频观看| 日韩一区二区三区影片| 伦精品一区二区三区| 亚洲av国产av综合av卡| 丰满迷人的少妇在线观看| 午夜激情福利司机影院| 国产伦理片在线播放av一区| 中文欧美无线码| 国产亚洲最大av| 国产 精品1| 亚洲三级黄色毛片| 国产黄色免费在线视频| 天堂8中文在线网| 国产精品一二三区在线看| 久久国产精品男人的天堂亚洲 | 少妇人妻精品综合一区二区| 国产在线男女| 国产午夜精品久久久久久一区二区三区| 舔av片在线| 久久精品国产亚洲网站| 国产毛片在线视频| 欧美日韩精品成人综合77777| av国产精品久久久久影院| 久热久热在线精品观看| 久久久久久久久久人人人人人人| 18禁裸乳无遮挡动漫免费视频| 婷婷色综合www| 国产久久久一区二区三区| 亚洲成色77777| 国产精品秋霞免费鲁丝片| 久久6这里有精品| 欧美成人精品欧美一级黄| 日韩国内少妇激情av| 我的老师免费观看完整版| av福利片在线观看| 亚洲av福利一区| 老司机影院成人| 一级爰片在线观看| 亚洲一区二区三区欧美精品| 久久综合国产亚洲精品| a级毛色黄片| 尾随美女入室| 免费播放大片免费观看视频在线观看| 插阴视频在线观看视频| 久热这里只有精品99| 精品人妻偷拍中文字幕| 久久ye,这里只有精品| 狂野欧美激情性bbbbbb| 国产成人freesex在线| av免费观看日本| 卡戴珊不雅视频在线播放| 精品亚洲乱码少妇综合久久| 我的老师免费观看完整版| 亚洲国产欧美人成| 亚洲av不卡在线观看| 欧美激情极品国产一区二区三区 | 最近最新中文字幕大全电影3| 亚洲欧美一区二区三区黑人 | 久久精品夜色国产| 国产成人freesex在线| 国内揄拍国产精品人妻在线| 最近手机中文字幕大全| 夜夜爽夜夜爽视频| 亚洲欧美精品专区久久| 亚洲va在线va天堂va国产| 又黄又爽又刺激的免费视频.| 寂寞人妻少妇视频99o| 国产成人精品婷婷| 一级av片app| 免费大片黄手机在线观看| 在线观看一区二区三区| 日韩亚洲欧美综合| 亚洲,欧美,日韩| 日韩精品有码人妻一区| 亚洲精品国产av成人精品| 大又大粗又爽又黄少妇毛片口| 交换朋友夫妻互换小说| 校园人妻丝袜中文字幕| 涩涩av久久男人的天堂| 亚洲精品亚洲一区二区| 欧美成人a在线观看| 极品少妇高潮喷水抽搐| 亚洲精品自拍成人| 日本爱情动作片www.在线观看| 亚洲不卡免费看| 一级毛片黄色毛片免费观看视频| a级毛片免费高清观看在线播放| 免费在线观看成人毛片| 亚洲经典国产精华液单| 亚洲国产欧美人成| 18禁动态无遮挡网站| 黄色配什么色好看| 一级片'在线观看视频| 日韩一区二区三区影片| 精品国产露脸久久av麻豆| xxx大片免费视频| 免费黄色在线免费观看| 欧美+日韩+精品| 日韩精品有码人妻一区| 黄色视频在线播放观看不卡| 久久精品国产亚洲av涩爱| 男女下面进入的视频免费午夜| 97热精品久久久久久| 夜夜爽夜夜爽视频| 26uuu在线亚洲综合色| 丰满少妇做爰视频| 欧美日本视频| 一级片'在线观看视频| 国产 精品1| 高清视频免费观看一区二区| 精品午夜福利在线看| 国产中年淑女户外野战色| 亚洲久久久国产精品| 精品一区二区免费观看| 国产乱来视频区| 亚洲中文av在线| 黑人高潮一二区| 51国产日韩欧美| 黄色配什么色好看| 爱豆传媒免费全集在线观看| 嘟嘟电影网在线观看| av天堂中文字幕网| 噜噜噜噜噜久久久久久91| 一级毛片我不卡| av黄色大香蕉| 人妻夜夜爽99麻豆av| 国产伦精品一区二区三区视频9| 在线 av 中文字幕| 日日摸夜夜添夜夜添av毛片| 美女国产视频在线观看| 最近的中文字幕免费完整| 国产一区二区三区综合在线观看 | 蜜桃在线观看..| 成人黄色视频免费在线看| 97热精品久久久久久| 午夜免费观看性视频| 亚洲天堂av无毛| 免费观看的影片在线观看| 亚洲人成网站在线观看播放| 亚洲国产精品专区欧美| 天堂俺去俺来也www色官网| 亚洲中文av在线| 午夜福利网站1000一区二区三区| 免费观看a级毛片全部| 卡戴珊不雅视频在线播放| 中文字幕久久专区| 国内揄拍国产精品人妻在线| 日韩一区二区三区影片| 午夜福利在线观看免费完整高清在| 亚洲精品国产色婷婷电影| 午夜精品国产一区二区电影| 91午夜精品亚洲一区二区三区| 日本爱情动作片www.在线观看| 国产色爽女视频免费观看| 日日啪夜夜撸| 中文字幕精品免费在线观看视频 | 久久97久久精品| 99热这里只有是精品50| 深夜a级毛片| 性色avwww在线观看| a级一级毛片免费在线观看| 精华霜和精华液先用哪个| 在现免费观看毛片| 内射极品少妇av片p| 色吧在线观看| 三级国产精品片| 久久女婷五月综合色啪小说| 一级二级三级毛片免费看| 老师上课跳d突然被开到最大视频| 亚洲精品亚洲一区二区| 热99国产精品久久久久久7| 国产伦在线观看视频一区| 各种免费的搞黄视频| 色视频www国产| 九九在线视频观看精品| 天美传媒精品一区二区| 成人午夜精彩视频在线观看| 黑人高潮一二区| 一级av片app| 在线看a的网站| 久久99热这里只频精品6学生| 丝袜喷水一区| 国产精品av视频在线免费观看| 欧美xxxx性猛交bbbb| 国产白丝娇喘喷水9色精品| 欧美另类一区| 又黄又爽又刺激的免费视频.| 午夜日本视频在线| 一本一本综合久久| 女性生殖器流出的白浆| 青春草亚洲视频在线观看| 又爽又黄a免费视频| 自拍偷自拍亚洲精品老妇| 最近手机中文字幕大全| 99国产精品免费福利视频| 黑人高潮一二区| 国产精品蜜桃在线观看| 亚洲精品久久久久久婷婷小说| 国模一区二区三区四区视频| 国产成人freesex在线| 免费久久久久久久精品成人欧美视频 | 五月天丁香电影| 热re99久久精品国产66热6| 欧美日韩一区二区视频在线观看视频在线| 校园人妻丝袜中文字幕| 日韩在线高清观看一区二区三区| 中文精品一卡2卡3卡4更新| 在线亚洲精品国产二区图片欧美 | 最近中文字幕2019免费版| 激情 狠狠 欧美| 国产日韩欧美亚洲二区| 久久精品国产自在天天线| 国产黄色视频一区二区在线观看| 在线看a的网站| 少妇猛男粗大的猛烈进出视频| 午夜视频国产福利| 亚洲av福利一区| 欧美成人午夜免费资源| 欧美日韩精品成人综合77777| 国产美女午夜福利| 日本黄大片高清| 国产av国产精品国产| 久久毛片免费看一区二区三区| 色吧在线观看| 日日啪夜夜爽| 亚洲av成人精品一二三区| 色视频www国产| 国产女主播在线喷水免费视频网站| 天天躁日日操中文字幕| 国产在线视频一区二区| 亚洲国产精品一区三区| 小蜜桃在线观看免费完整版高清| 99久久中文字幕三级久久日本| 超碰97精品在线观看| 天天躁夜夜躁狠狠久久av| 丰满少妇做爰视频| 欧美最新免费一区二区三区| 女人十人毛片免费观看3o分钟| 国产又色又爽无遮挡免| 亚洲欧洲日产国产| 国产 一区 欧美 日韩| 一区在线观看完整版| 国产精品一区www在线观看| 欧美区成人在线视频| 国精品久久久久久国模美| 亚洲精品一二三| 一级片'在线观看视频| 又爽又黄a免费视频| 美女主播在线视频| 三级国产精品片| 91久久精品电影网| 我要看黄色一级片免费的| 亚洲中文av在线| 视频中文字幕在线观看| 国产欧美日韩精品一区二区| 国产精品一及| 国产大屁股一区二区在线视频| 欧美亚洲 丝袜 人妻 在线| 五月天丁香电影| 欧美成人一区二区免费高清观看| 亚洲美女搞黄在线观看| 中文欧美无线码| 久久久久国产精品人妻一区二区| 久久人妻熟女aⅴ| 国产精品久久久久久av不卡| 成人免费观看视频高清| 少妇裸体淫交视频免费看高清| 少妇的逼水好多| 18禁裸乳无遮挡动漫免费视频| 亚洲美女搞黄在线观看| 免费播放大片免费观看视频在线观看| 中文字幕亚洲精品专区| 亚洲精品日韩av片在线观看| 一级av片app| 日日啪夜夜撸| a级毛色黄片| 亚洲激情五月婷婷啪啪| 国产人妻一区二区三区在| 亚洲第一区二区三区不卡| av国产久精品久网站免费入址| 国产在线视频一区二区| av在线播放精品| 欧美+日韩+精品| 日韩强制内射视频| 国产在线男女| 日韩强制内射视频| 久久亚洲国产成人精品v| 国国产精品蜜臀av免费| 亚洲av中文字字幕乱码综合| 男女国产视频网站| 蜜桃在线观看..| 丰满乱子伦码专区| 国产av一区二区精品久久 | 亚洲美女搞黄在线观看| 日本wwww免费看| 国产精品久久久久久久久免| 亚洲成人中文字幕在线播放| 成人二区视频| 中文欧美无线码| 99热这里只有是精品在线观看| 日韩人妻高清精品专区| 最近的中文字幕免费完整| 成人毛片a级毛片在线播放| 一本一本综合久久| 边亲边吃奶的免费视频| 国内精品宾馆在线| 国产精品一区二区性色av| 国产亚洲精品久久久com| 熟女电影av网| 欧美成人一区二区免费高清观看| 国产熟女欧美一区二区| 日本欧美国产在线视频| 亚洲精品aⅴ在线观看| av不卡在线播放| 亚洲国产高清在线一区二区三| 日韩欧美精品免费久久| 亚洲经典国产精华液单| av卡一久久| www.av在线官网国产| 人人妻人人添人人爽欧美一区卜 | 免费av中文字幕在线| 久久精品熟女亚洲av麻豆精品| 中文精品一卡2卡3卡4更新| 国产成人a区在线观看| 丰满人妻一区二区三区视频av| 高清在线视频一区二区三区| 色吧在线观看| 少妇精品久久久久久久| 大片电影免费在线观看免费| 精品久久久噜噜| 国产一区亚洲一区在线观看| 亚洲丝袜综合中文字幕| 欧美国产精品一级二级三级 | 亚洲av国产av综合av卡| 国产亚洲一区二区精品| 亚洲欧美日韩无卡精品| 国产 精品1| 男人添女人高潮全过程视频| 色哟哟·www| 成人综合一区亚洲| 欧美国产精品一级二级三级 | 亚洲精品亚洲一区二区| 天堂俺去俺来也www色官网| 在线观看免费日韩欧美大片 | 国产精品99久久99久久久不卡 | 三级国产精品欧美在线观看| 简卡轻食公司| 亚洲真实伦在线观看| 久久99热6这里只有精品| av播播在线观看一区| 99热这里只有是精品在线观看| 国产黄频视频在线观看| 欧美 日韩 精品 国产| 精品国产乱码久久久久久小说| 18禁动态无遮挡网站| 日韩精品有码人妻一区| 精品亚洲成a人片在线观看 | 欧美+日韩+精品| 精品少妇久久久久久888优播| 国产一区二区三区综合在线观看 | 国产女主播在线喷水免费视频网站| 亚洲欧美中文字幕日韩二区| 亚洲av国产av综合av卡| 亚洲精品456在线播放app| 少妇被粗大猛烈的视频| 九九爱精品视频在线观看| 精品久久久久久电影网| 精品久久久噜噜| 亚洲国产欧美人成| 国产成人一区二区在线| 欧美zozozo另类| h日本视频在线播放| 亚洲婷婷狠狠爱综合网| 成年美女黄网站色视频大全免费 | 最近手机中文字幕大全| 成人毛片a级毛片在线播放| 少妇人妻精品综合一区二区| 久久国产乱子免费精品| 日日啪夜夜爽| 九色成人免费人妻av| 小蜜桃在线观看免费完整版高清| 色视频在线一区二区三区| 国产成人精品福利久久| 美女cb高潮喷水在线观看| 永久免费av网站大全| 国产黄色视频一区二区在线观看| 亚洲三级黄色毛片| 欧美日韩在线观看h| 国产v大片淫在线免费观看| 欧美bdsm另类| 一个人免费看片子| 97在线视频观看| 日韩成人伦理影院| 日日摸夜夜添夜夜爱| .国产精品久久| av不卡在线播放| 日韩人妻高清精品专区| 国产精品伦人一区二区| 国产成人91sexporn| 极品少妇高潮喷水抽搐| 国产精品三级大全| 亚洲在久久综合| 啦啦啦视频在线资源免费观看| 久久久久国产网址| 好男人视频免费观看在线| 久久人妻熟女aⅴ| 熟女人妻精品中文字幕| 亚洲欧美日韩无卡精品| 尤物成人国产欧美一区二区三区| 国产精品免费大片| 99re6热这里在线精品视频| 欧美成人a在线观看| 久久人妻熟女aⅴ| 狂野欧美白嫩少妇大欣赏| 日本av免费视频播放| 高清毛片免费看| 精品午夜福利在线看| 国产又色又爽无遮挡免| 久久国产亚洲av麻豆专区| 99精国产麻豆久久婷婷| 国产 精品1| 国产精品无大码| 久久久久久久久久久丰满| 免费不卡的大黄色大毛片视频在线观看| 国产成人精品一,二区| 麻豆国产97在线/欧美| 97超碰精品成人国产| 日本爱情动作片www.在线观看| av在线蜜桃| 日韩一区二区三区影片| 亚洲美女视频黄频| 大又大粗又爽又黄少妇毛片口| 久久影院123| 久久久久国产网址| 日本欧美国产在线视频| 久久99精品国语久久久| 久久精品国产亚洲av天美| 99久久精品热视频| 午夜免费观看性视频| 久久久久久久精品精品| 成人免费观看视频高清| 国产片特级美女逼逼视频| 午夜福利在线在线| 99久久人妻综合| 又粗又硬又长又爽又黄的视频| 亚洲精华国产精华液的使用体验| 嫩草影院入口| 日韩人妻高清精品专区| 亚洲精品国产av蜜桃| 黄片wwwwww| 最近最新中文字幕大全电影3| 激情五月婷婷亚洲| 九草在线视频观看| 亚洲av免费高清在线观看| 只有这里有精品99| 尤物成人国产欧美一区二区三区| av一本久久久久| 亚洲色图av天堂| 亚洲国产精品成人久久小说| 国产精品福利在线免费观看| 乱码一卡2卡4卡精品| 国产成人精品婷婷| 免费看光身美女| 夜夜骑夜夜射夜夜干| 免费高清在线观看视频在线观看| 直男gayav资源| 男女边吃奶边做爰视频| 伦理电影大哥的女人| 啦啦啦中文免费视频观看日本| 免费黄网站久久成人精品| 亚洲一级一片aⅴ在线观看| 小蜜桃在线观看免费完整版高清| 毛片一级片免费看久久久久| 精品人妻偷拍中文字幕| 欧美高清性xxxxhd video| 亚洲av福利一区| 99re6热这里在线精品视频| 春色校园在线视频观看| 只有这里有精品99| 国产 一区 欧美 日韩| 狂野欧美激情性xxxx在线观看| 性高湖久久久久久久久免费观看| 国产 精品1| 欧美日韩视频精品一区| 中文字幕人妻熟人妻熟丝袜美| 国产视频首页在线观看| 大陆偷拍与自拍| 韩国av在线不卡| 菩萨蛮人人尽说江南好唐韦庄| 80岁老熟妇乱子伦牲交| 欧美zozozo另类| 中文天堂在线官网| 熟妇人妻不卡中文字幕| 国产精品三级大全| 亚洲经典国产精华液单| 久久久久国产精品人妻一区二区| 国产伦精品一区二区三区视频9| 国产在线男女| 欧美成人一区二区免费高清观看| 亚洲精品久久久久久婷婷小说| 极品少妇高潮喷水抽搐| 美女xxoo啪啪120秒动态图| 全区人妻精品视频| 久久精品久久久久久久性| 爱豆传媒免费全集在线观看| 国产探花极品一区二区| 亚洲av免费高清在线观看| 免费观看a级毛片全部| 青春草视频在线免费观看| 国产综合精华液| 亚洲国产高清在线一区二区三| 欧美3d第一页| 午夜激情久久久久久久| 大香蕉久久网| 五月天丁香电影| 看非洲黑人一级黄片| 美女中出高潮动态图| 校园人妻丝袜中文字幕| 寂寞人妻少妇视频99o| 国产精品蜜桃在线观看| 亚洲国产欧美人成| 久久久久久久久久人人人人人人| 天堂俺去俺来也www色官网| 黑人高潮一二区| 日韩,欧美,国产一区二区三区| 男人和女人高潮做爰伦理| 国产一区二区在线观看日韩| a级一级毛片免费在线观看| 国产免费一区二区三区四区乱码| 免费人妻精品一区二区三区视频| 国产色爽女视频免费观看| 校园人妻丝袜中文字幕| 中文天堂在线官网| 99久久人妻综合| 一级二级三级毛片免费看| 美女主播在线视频| 精品亚洲成a人片在线观看 | 少妇的逼好多水| 欧美3d第一页| 制服丝袜香蕉在线| 国产视频内射| 精品久久久久久久久亚洲| 国产日韩欧美亚洲二区| 少妇的逼好多水| 人人妻人人爽人人添夜夜欢视频 | 精品午夜福利在线看| 日韩电影二区| 免费大片18禁| 国产av精品麻豆| 亚洲欧洲日产国产| 亚洲,欧美,日韩| 国语对白做爰xxxⅹ性视频网站| 亚洲精品aⅴ在线观看| av女优亚洲男人天堂| 身体一侧抽搐| 国内精品宾馆在线| 亚洲av国产av综合av卡| 国产精品国产三级国产专区5o| 男女啪啪激烈高潮av片| 在现免费观看毛片| 亚洲av欧美aⅴ国产| 99久久精品国产国产毛片| 亚洲精品自拍成人| 99久久人妻综合| 亚洲精品国产色婷婷电影| 一个人看的www免费观看视频| 麻豆成人午夜福利视频| 日本av免费视频播放| 老师上课跳d突然被开到最大视频| 男女国产视频网站| 亚洲中文av在线| 黄色视频在线播放观看不卡| 在线观看一区二区三区| 亚洲三级黄色毛片| 在线观看美女被高潮喷水网站|