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

    庫布齊沙地土壤水分運(yùn)移規(guī)律數(shù)值模擬研究

    2023-06-12 10:02:12楊建軍孫貴榮丁俊峰王涵孫凱阿拉騰蘇和
    林業(yè)科技 2023年2期
    關(guān)鍵詞:設(shè)置研究

    楊建軍 孫貴榮 丁俊峰 王涵 孫凱 阿拉騰蘇和

    摘要: ?以砂土下進(jìn)行地表滴灌的沙柳、速生楊、旱柳和榆樹人工林為研究對(duì)象,采用粒徑檢測(cè)激光分析儀進(jìn)行土壤質(zhì)地分析得到各土樣粒徑分布情況,通過HYDRUS-2D/3D數(shù)值模擬研究灌水及水分再分布過程中壓力水頭和含水量的動(dòng)態(tài)變化規(guī)律。結(jié)果表明:庫布齊沙漠北部邊緣地區(qū)人工林立地土壤中砂粒含量高于70%且黏粒含量低于2.5%,屬典型砂質(zhì)土壤。滴灌流量為3.0 L/h,時(shí)長為6 h時(shí),各類人工林根深范圍內(nèi)土壤均能達(dá)到飽和含水量,能滿足人工林用水需求。

    關(guān)鍵詞: ?滴灌; ?水分運(yùn)移; ?壓力水頭; ?含水量; ?HYDRUS-2D/3D

    中圖分類號(hào): ? S 152. 7 ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)識(shí)碼: ? A ? ? ? ? ? ? ? ?文章編號(hào):1001 - 9499(2023)02 - 0028 - 05

    Numerical Simulation of Soil Water Movement

    in the Kubuqi Sandy Land

    YANG Jianjun1 SUN Guirong1 WANG Han2 DING Junfeng2 SUN Kai2 Alatengsuhe1

    (1. ?Erdos afforestation general plant, Inner Mongolia Shulinzhao 014300; ?2. ?College of energy and transportation engineering, ?Inner Mongolia Agricultural University, ?Inner Mongolia Hohhot 010018)

    Abstract Taking Salix psammophila, Populus popular, Pinus sylvestris and Ulmus pumila for surface drip irrigation under sandy soil as the research object, the particle size distribution of each soil sample was obtained by soil texture analysis using particle size detection laser analyzer. The dynamic change laws of pressure head, soil water content and infiltration rate during irrigation and water redistribution were studied through HYDRUS-2D/3D numerical simulation. The results showed that the sand content in the soil of artificial forest sites in the northern edge of Kubuqi Desert was higher than 70% and the clay content was lower than 2.5%, which was typical sandy soil. When the flow rate of drip irrigation is 3.0 L/h and the duration is 6 h, the soil within the root depth of various plantations can reach the saturated water content, which can meet the water demand of plantations.

    Key words drip irrigation; water movement; pressure head; water content; HYDRUS-2D/3D

    滴灌是可根據(jù)植物需水規(guī)律及時(shí)補(bǔ)充根部水分的現(xiàn)代節(jié)水灌溉技術(shù),也是人工林維持生產(chǎn)力的重要技術(shù)措施[ 1 - 3 ]。眾多研究表明,與其他灌溉方式相比,滴灌能顯著提高植物水分利用效率達(dá)14%~35%[ 4 -6 ]。滴灌通過布置的管道系統(tǒng)使水流緩緩滴到植物體的根部和土壤中,形成土壤濕潤體供作物吸收大量的水分和養(yǎng)分,因而土壤濕潤體的大小和形狀能對(duì)人工林生長和發(fā)育產(chǎn)生重要影響[ 7 - 8 ]。因此,研究滴灌下不同人工林土壤水分運(yùn)移規(guī)律,揭示不同林地土壤濕潤體時(shí)空變化特征,是制定高效灌溉策略的重要依據(jù)。

    通??赏ㄟ^研究土壤壓力水頭和含水量的時(shí)空變化規(guī)律開展土壤水分運(yùn)移規(guī)律的田間試驗(yàn)研究,但田間試驗(yàn)具有耗時(shí)長、投入人力物力大、可研究情形有限等缺點(diǎn)。因此,基于田間試驗(yàn)資料,進(jìn)行土壤水分運(yùn)動(dòng)建模,通過數(shù)值模擬計(jì)算得到較優(yōu)的滴灌模式和參數(shù),是較為經(jīng)濟(jì)、高效的研究方法。近年來,有學(xué)者采用軟件數(shù)值模擬法研究滴灌條件下的土壤水分運(yùn)移規(guī)律,如Hiba Ghazouani等[ 9 ]利用HYDRUS-2D建立馬鈴薯(Potato)田研究模型,對(duì)比了灌溉量、灌溉水質(zhì)對(duì)土壤含水量和電導(dǎo)率的影響;李豆豆等[ 1 ]利用HYDRUS-2D/3D對(duì)毛白楊(Populus tomentosa)人工林砂壤土建模,研究了不同初始土壤含水量下濕潤鋒水平和垂直運(yùn)移距離隨灌溉時(shí)長的變化規(guī)律;黃凱等[ 10 ]利用HYDRUS 對(duì)赤紅壤甘蔗( Saccharum officinarum) 田建模分析了滴頭間距、流量對(duì)灌水均勻性的影響;Naglic等[ 11 ]利用HYDRUS-2D/3D 針對(duì)不同質(zhì)地土壤模擬了滴頭流量和初始含水量對(duì)濕潤體的影響。然而,利用HYDRUS 系列軟件模擬滴灌下土壤水分運(yùn)動(dòng)的研究多集中在毛白楊、甘蔗、小麥、葡萄等產(chǎn)值較高或單一作物田地及不同質(zhì)地裸地上[ 12 - 15 ],而在沙柳(Salix psammophila)、速生楊(Populus popu- ?lar)、樟子松(Pinus sylvestris)、榆樹(Ulmus pumila)等多種人工林上的對(duì)比研究卻鮮有涉及。

    人工林可通過改變土壤容重、土壤孔隙度,從而增強(qiáng)土壤入滲性能,使土壤儲(chǔ)水能力增加。但另一方面,庫布齊沙漠北部邊緣地區(qū)人工造林成活后,由于蒸騰耗水作用強(qiáng)烈,沙地水分被大量消耗,沙地水分與植物生長供需嚴(yán)重失調(diào),水分成為限制人工林生長的主要影響因子,人工林出現(xiàn)衰退和死亡現(xiàn)象,使得精準(zhǔn)滴灌顯得尤為重要。因此,在內(nèi)蒙古庫布齊沙漠北部邊緣地區(qū)人工林典型栽植立地砂壤土上,開展地表滴灌下不同人工林林地土壤水分運(yùn)移的數(shù)值模擬研究,進(jìn)而提出灌溉策略建議,對(duì)于沙柳等人工林生產(chǎn)力的快速提升具有重要意義。本研究對(duì)不同人工林下滴灌土壤水分運(yùn)移規(guī)律開展研究,在庫布齊沙漠北部邊緣地區(qū)砂壤土立地下,探討地表滴灌人工林的土壤水分運(yùn)移規(guī)律,通過HYDRUS-2D/3D模擬試驗(yàn)情境,得到不同林地土壤壓力水頭和含水量的時(shí)空變化特征,以期明確滴灌下各類林地土壤水分運(yùn)移規(guī)律,為不同林分制訂精準(zhǔn)滴灌策略提供數(shù)據(jù)支撐和理論依據(jù)。

    1 研究區(qū)概況

    試驗(yàn)地位于內(nèi)蒙古鄂爾多斯市達(dá)拉特旗鄂爾多斯林場(chǎng)曹四灘護(hù)林站(110°39'14″E,40°14'24″N)(圖1)。該地屬溫帶大陸性季風(fēng)氣候,其顯著特征是常年干燥多風(fēng)少雨,多年平均氣溫6.1 ℃;全年無霜期135天,年平均日照時(shí)數(shù)3 100 h,年平均降水量246.9~310.8 mm,降雨集中在每年7~9月,7月最多。人工林主要樹種為沙柳(Salix psammophila)、速生楊(Populus popular)、旱柳(Pinus sylvestris)及榆樹(Ulmus pumila),人工林下草本植被以沙生植被為主。

    于2021年8月初在研究區(qū)采集土樣,研究區(qū)內(nèi)包含沙柳林、速生楊林、旱柳和榆樹,針對(duì)每種林木各選擇一塊具有代表性的試驗(yàn)樣地,其林木均勻配置。取土樣方設(shè)置在各試驗(yàn)小區(qū)相鄰的林木林帶之間,長2.0 m,寬1.0 m,深1.2 m,對(duì)樣方內(nèi)的土壤進(jìn)行剖面分層取土。由于研究區(qū)的地下水位較低,且速生楊、榆樹、旱柳等喬木根系較深,故取土深度統(tǒng)一設(shè)置為100 cm,取樣間隔為20 cm,即0~20、20~40、40~60、60~80、80~100 cm,分別用直徑為5.046 cm、高為5cm、體積為100 cm3環(huán)刀分層取土,每層隨機(jī)取3個(gè)重復(fù)土樣,采用德國SYMPATEC GmbH公司的RODOS/M粒徑檢測(cè)激光分析儀進(jìn)行土壤質(zhì)地分析得到各土樣粒徑分布情況,并計(jì)算3個(gè)土樣均值(表1)。

    2 研究方法

    2. 1 數(shù)值模擬基礎(chǔ)條件設(shè)計(jì)

    (1)模擬基本環(huán)境

    本次模擬的是距樹30 cm的干旱未灌溉條件,滴灌裝置選取為滴頭流量3 L/h的滴灌帶,所有土壤的模擬滴灌時(shí)長均為6 h。根據(jù)土樣采集基本數(shù)據(jù),滴灌模擬土壤的總深度設(shè)置為地下0~100 cm,以0~20、20~40、40~60、60~80和80~100 cm分別設(shè)置為單一層級(jí)的土壤,共設(shè)置5層。

    (2)時(shí)間離散設(shè)置

    在進(jìn)行模型的迭代運(yùn)算時(shí),離散時(shí)間單位取小時(shí)(h)。其中所選用的最短時(shí)間步長0.001 h;在允許偏差下,土壤含水量迭代精度設(shè)置為0.001,壓力水頭的迭代精度設(shè)置為1;迭代的最大次數(shù)為10,若迭代次數(shù)超過10,則系統(tǒng)將用新時(shí)間步長取代原有時(shí)間步長的1/3重新開始有關(guān)迭代計(jì)算。除此之外,本次模擬增大迭代的比例設(shè)置為1.3,增大步長的迭代信號(hào)設(shè)置為3,縮小迭代的比例設(shè)置為0.7,縮小步長的迭代信號(hào)設(shè)置為7。在網(wǎng)格劃分和離散設(shè)置中四類人工林土壤結(jié)構(gòu)選取范圍相同。

    2. 2 土壤環(huán)境設(shè)置

    (1)土壤粒徑分布參數(shù)

    各類人工林中土壤砂粒、粉粒和黏粒含量的體積百分比如表1所示,按該參數(shù)對(duì)數(shù)值模擬土壤進(jìn)行分層設(shè)置。

    (2)土壤水動(dòng)力參數(shù)設(shè)置

    根據(jù)表1中5種土壤各土層顆粒組成,利用Rosetta網(wǎng)絡(luò)神經(jīng)模塊,初步擬合并逐步計(jì)算推導(dǎo)出各種土壤水動(dòng)力參數(shù)如表2。

    (3)初始條件及邊界條件設(shè)置

    模型區(qū)域:本次模擬將模擬區(qū)域設(shè)置成寬為100 cm、深為100 cm的正方形區(qū)域。

    初始條件:以進(jìn)水前測(cè)量結(jié)果擬定出來的含水量剖面圖作為進(jìn)行本次模擬計(jì)算的初始條件之一,并以此假設(shè)含水量能在同一厚度土層內(nèi)均勻分布。

    邊界條件:在模擬中的滴頭處預(yù)先設(shè)置變流量邊界,在開始進(jìn)水后以其周圍半徑為W的區(qū)域?yàn)轱柡蛥^(qū),且在進(jìn)行模擬試驗(yàn)時(shí),可以忽略飽和區(qū)半徑和進(jìn)水時(shí)間的關(guān)系。

    下排水邊界條件和其余邊界條件:根據(jù)研究區(qū)地下水位平均埋深約為10 m,距地面最近地下水埋深為5.26 m。而由于模擬深度為1 m,因此模擬區(qū)域下排水邊界應(yīng)被設(shè)定為自由排水邊界;其余所有未預(yù)先設(shè)置的邊界均默認(rèn)保存為零通量邊界。

    (4)觀測(cè)點(diǎn)選取

    依據(jù)土壤分層設(shè)置的特點(diǎn),所選觀測(cè)點(diǎn)深度分別為20、40、60、80和100 cm,自上而下依次為觀測(cè)點(diǎn)1~5。

    3 結(jié)果與分析

    3. 1 觀測(cè)點(diǎn)壓力水頭變化

    根據(jù)HYDRUS-2D/3D模擬結(jié)果,依據(jù)圖2,灌溉開始前各人工林任意觀測(cè)點(diǎn)的起始?jí)毫λ^均為-100 cm,滴灌6 h后其數(shù)值均超過0 cm,且在之后穩(wěn)定不變。可見此次滴灌的結(jié)果是100 cm土壤深度范圍內(nèi)不同土層水分均達(dá)到穩(wěn)定壓力水頭值,且該值大于0 cm,隨土層加深該壓力水頭值逐漸減小,同時(shí)不同人工林水分運(yùn)移的時(shí)空變化規(guī)律存在一定差異。

    對(duì)比不同人工林,速生楊林各觀測(cè)點(diǎn)均達(dá)到穩(wěn)定狀態(tài)用時(shí)最短,約為1.5 h,且各土層最終形成的穩(wěn)定壓力水頭值最大,為97~152 cm,沙柳林和旱柳林次之,榆樹林各土層達(dá)到穩(wěn)定狀態(tài)用時(shí)最長,約為5.7 h,且穩(wěn)定壓力水頭值最小,為0.2~8.8 cm。上述結(jié)果與該四類人工林土壤質(zhì)地有關(guān),由表2可知,速生楊林各土層中砂粒和粉粒含量之比(以下簡稱“砂粉比”)均值最大,而沙柳林和旱柳林次之,榆樹林土壤砂粉比最小。可見土壤顆粒分布對(duì)其滴灌過程水分運(yùn)移規(guī)律影響較大,土壤質(zhì)地越粗則水分運(yùn)移中土壤中各土層水分越能快速趨于穩(wěn)定,且最終各土層壓力水頭也越大,反之質(zhì)地越細(xì)則速度越慢且最終壓力水頭也越小。

    3. 2 觀測(cè)點(diǎn)含水量變化

    依據(jù)圖3,灌溉開始前各人工林任意觀測(cè)點(diǎn)的起始含水量均為0.05~0.20 cm3/cm3,高于其殘余含水量(表2),滴灌6 h后其數(shù)值均接近0.40 cm3/cm3,且在之后穩(wěn)定不變。由表2可知,滴灌結(jié)果為100 cm土壤深度范圍內(nèi)各土層水分均達(dá)到各自飽和含水量,該值與土壤質(zhì)地有關(guān)且非常接近,該特點(diǎn)與壓力水頭有明顯差別,可見相同水分運(yùn)移過程中土壤壓力水頭和含水量的時(shí)空變化規(guī)律存在差異。

    對(duì)比圖2和圖3,不同人工林各觀測(cè)點(diǎn)含水量均達(dá)到飽和含水量與壓力水頭達(dá)到穩(wěn)定狀態(tài)用時(shí)相同,仍呈現(xiàn)速生楊<沙柳<旱柳<榆樹的特點(diǎn),同時(shí)各土層均隨其初始含水量的不同而有微小差異,但這對(duì)不同土層含水量達(dá)到飽和含水量的時(shí)間影響較小。以上規(guī)律進(jìn)一步驗(yàn)證了土壤砂粉比對(duì)水分運(yùn)移過程中含水量時(shí)空變化規(guī)律的影響,土壤質(zhì)地越粗則水分運(yùn)移中土壤中各土層含水量越能快速趨于飽和,反之質(zhì)地越細(xì)則速度越慢,而最終各土層含水量主要由其飽和含水量決定。

    4 結(jié) 論

    本文采用數(shù)值模擬的方法對(duì)庫布齊沙地不同人工林下滴灌土壤水分運(yùn)移規(guī)律開展研究,主要得出如下結(jié)論:

    (1)庫布齊沙漠北部邊緣地區(qū)人工林立地土壤中砂粒含量高于70%且黏粒含量低于2.5%,屬于典型砂質(zhì)土壤。滴灌流量為3.0 L/h,時(shí)長為6 h時(shí),該地區(qū)各類人工林100 cm土層深度范圍內(nèi)土壤均能達(dá)到飽和含水量,能滿足人工林用水需求。

    (2)庫布齊沙地土壤質(zhì)地下,粒徑分布尤其是砂粒和粉粒含量之比對(duì)滴灌過程中水分運(yùn)移結(jié)果影響較大,具體表現(xiàn)為砂粉比越大則人工林各土層達(dá)到壓力水頭最大值和土壤飽和含水量的時(shí)間越短。上述規(guī)律與已有文獻(xiàn)結(jié)論類似,表明采用數(shù)值模擬的手段研究土壤多孔介質(zhì)水分運(yùn)移規(guī)律切實(shí)可行。

    (3)當(dāng)通過滴灌使一定深度范圍土壤均達(dá)到飽和含水量時(shí),各土層土壤含水量值較為接近,但壓力水頭值卻與土壤砂粉比和土層深度有關(guān)。隨土層加深,水分運(yùn)移至深土層處水勢(shì)梯度相應(yīng)減小,滴灌后穩(wěn)定的壓力水頭值也減小。

    參考文獻(xiàn)

    [1] 李豆豆, ?席本野, ?唐連峰, ?等. ?砂壤土下滴灌毛白楊幼林土壤水分運(yùn)移規(guī)律與模擬[J]. ?林業(yè)科學(xué), ?2018, 54(12): 157 - 168.

    [2] 秘洪雷, ?蘭再平, ?彭晶晶, ?等. ?滴灌與溝灌對(duì)楊樹細(xì)根空間分布的影響[J]. ?浙江農(nóng)林大學(xué)學(xué)報(bào): ?2021, 38(X): 1 - 8.

    [3] 張玉豪, ?姚素梅, ?孟麗, ? 等. ?滴灌土壤基質(zhì)勢(shì)調(diào)控對(duì)紅豆杉生理特性及次生代謝物含量的影響[J]. ?干旱地區(qū)農(nóng)業(yè)研究,2021, 39(5): 129 - 137.

    [4] Bhunia S R, Verma I M, Sahu M P, et al. Effect of drip irrigation and bioregulators on yield, economics and water use of fenugreek ( Trigonella foenum-graecum) [J]. Journal of Spices and Aromatic Crops, 2015, 24(2): 102 - 105

    [5] 李久生, ?栗巖峰, ?王軍, ?等. ?微灌在中國: ?歷史,現(xiàn)狀和未來[J]. ?水利學(xué)報(bào),2016, ?47(3): 372 - 381.

    [6] 賈詠霖, ?屈忠義, ?丁艷宏, ?等. ?不同灌溉方式下施用生物炭對(duì)土壤水鹽運(yùn)移規(guī)律及玉米水分利用效率的影響[J]. ?灌溉排水學(xué)報(bào), ?2020, 39(8): 44 - 51.

    [7] Barragn J, Bralts V, Wu I P.Assessment of emission uniformity for micro-irrigation design. Biosystems Engineering[J]. 2006, 93(1) : 89 - 97.

    [8] 傅建平, ?蘭再平, ?孫尚偉, ?等. ?滴灌條件下楊樹人工林土壤的水分運(yùn)移[J]. ?林業(yè)科學(xué), ?2013, 49(6): 25 - 29.

    [9] Hiba Ghazouani, Giovanni Rallo, Amel Mguidiche, et al.Assessing Hydrus-2D Model to Investigate the Effects of Different On-Farm Irrigation Strategies on Potato Crop under Subsurface Drip Irrigation[J]. Water, 2019, 11(3): 85 - 89.

    [10] 黃凱, ?蔡德所, ?潘偉, ?等. ?廣西赤紅壤甘蔗田間滴灌帶合理布設(shè)參數(shù)確定[J]. ?農(nóng)業(yè)工程學(xué)報(bào), ?2015, 31(11): 136 - 143.

    [11] Naglic B, Kechavarzi C, Coulon F, et al. Numerical investigation of the influence of texture,surface drip emitter discharge rate and initial soil moisture condition on wetting pattern size[J]. Irrigation Science, 2014, 32(6): 421 - 436.

    [12] 李明思, ?康紹忠, ?孫海燕. ?點(diǎn)源滴灌滴頭流量與濕潤體關(guān)系研究[J]. ?農(nóng)業(yè)工程學(xué)報(bào), ?2006, 22(4): 32 - 35.

    [13] Wang F X, Kang Y, Liu S P. Effects of drip irrigation frequency on soil wetting pattern and potato growth in North China Plain[J].Agricultural Water Management, 2006, 79(3): 248 - 264.

    [14] 王振華, ?鄭旭榮, ?任杰, ?等. ?灌水頻率對(duì)地下滴灌線源入滲土壤水分運(yùn)動(dòng)規(guī)律影響的試驗(yàn)研究[J]. 水利學(xué)報(bào), ?2007(s1): 299 - 302.

    [15] McClymont L, Goodwin I. Effects of wetted soil volume on young pear trees[J]. Acta Horticulturae, 2014, 11(30): 479 - 484.

    猜你喜歡
    設(shè)置研究
    FMS與YBT相關(guān)性的實(shí)證研究
    2020年國內(nèi)翻譯研究述評(píng)
    遼代千人邑研究述論
    中隊(duì)崗位該如何設(shè)置
    視錯(cuò)覺在平面設(shè)計(jì)中的應(yīng)用與研究
    科技傳播(2019年22期)2020-01-14 03:06:54
    EMA伺服控制系統(tǒng)研究
    7招教你手動(dòng)設(shè)置參數(shù)
    新版C-NCAP側(cè)面碰撞假人損傷研究
    本刊欄目設(shè)置說明
    中俄臨床醫(yī)學(xué)專業(yè)課程設(shè)置的比較與思考
    国产午夜精品久久久久久一区二区三区 | 亚洲精品国产成人久久av| 免费电影在线观看免费观看| 一夜夜www| av在线播放精品| 毛片女人毛片| 黄色欧美视频在线观看| 麻豆久久精品国产亚洲av| 婷婷六月久久综合丁香| 一区二区三区免费毛片| 日本在线视频免费播放| 国产老妇女一区| 美女大奶头视频| 久久久精品94久久精品| 春色校园在线视频观看| 亚洲av成人精品一区久久| 亚洲专区国产一区二区| 男女做爰动态图高潮gif福利片| 最近中文字幕高清免费大全6| 精品无人区乱码1区二区| 精品国内亚洲2022精品成人| av.在线天堂| 久久精品国产亚洲av香蕉五月| а√天堂www在线а√下载| 秋霞在线观看毛片| 大又大粗又爽又黄少妇毛片口| 在线天堂最新版资源| 亚洲美女视频黄频| 日韩欧美精品免费久久| 99热这里只有是精品50| 亚洲成人久久性| 免费不卡的大黄色大毛片视频在线观看 | 午夜精品一区二区三区免费看| 小说图片视频综合网站| 国产成年人精品一区二区| 午夜福利在线观看吧| 日本三级黄在线观看| 婷婷亚洲欧美| 亚洲av.av天堂| av卡一久久| 蜜桃亚洲精品一区二区三区| 欧美日韩国产亚洲二区| 婷婷色综合大香蕉| av在线播放精品| 欧美成人精品欧美一级黄| 国产午夜精品论理片| 一进一出抽搐gif免费好疼| 国产男人的电影天堂91| 狂野欧美白嫩少妇大欣赏| 国产精品电影一区二区三区| 精品一区二区三区人妻视频| 亚洲国产欧美人成| 美女高潮的动态| 真人做人爱边吃奶动态| 婷婷亚洲欧美| 夜夜夜夜夜久久久久| 国产成人a∨麻豆精品| 美女内射精品一级片tv| 一级毛片我不卡| 成人无遮挡网站| 国产亚洲精品久久久com| 精品一区二区三区视频在线| 日韩 亚洲 欧美在线| 一级av片app| 特大巨黑吊av在线直播| av在线老鸭窝| 日韩在线高清观看一区二区三区| 成人永久免费在线观看视频| 尾随美女入室| 婷婷精品国产亚洲av| 亚洲一区二区三区色噜噜| 日本五十路高清| 男女下面进入的视频免费午夜| 免费不卡的大黄色大毛片视频在线观看 | 色噜噜av男人的天堂激情| 丰满乱子伦码专区| 国模一区二区三区四区视频| 免费看av在线观看网站| 国产精品久久久久久久久免| 亚洲七黄色美女视频| 国产精品女同一区二区软件| 1024手机看黄色片| 99热这里只有精品一区| 波多野结衣巨乳人妻| 美女被艹到高潮喷水动态| 亚洲电影在线观看av| 国产伦精品一区二区三区视频9| 嫩草影院精品99| 男人的好看免费观看在线视频| 简卡轻食公司| 美女高潮的动态| 亚洲天堂国产精品一区在线| 免费无遮挡裸体视频| 少妇人妻一区二区三区视频| 精品久久国产蜜桃| 嫩草影院新地址| 精品人妻熟女av久视频| 国产又黄又爽又无遮挡在线| 色综合色国产| av国产免费在线观看| 麻豆国产97在线/欧美| 嫩草影院入口| 欧美xxxx性猛交bbbb| 国产日本99.免费观看| 精品久久久久久久久久久久久| 国产在线男女| 久久久久免费精品人妻一区二区| 成人av一区二区三区在线看| 香蕉av资源在线| 最近中文字幕高清免费大全6| 亚洲人与动物交配视频| 欧洲精品卡2卡3卡4卡5卡区| 日日摸夜夜添夜夜添av毛片| 亚洲国产高清在线一区二区三| 91av网一区二区| 成人漫画全彩无遮挡| 最近视频中文字幕2019在线8| 亚洲精品影视一区二区三区av| 99在线人妻在线中文字幕| 九九热线精品视视频播放| 深爱激情五月婷婷| 国产淫片久久久久久久久| 久久久国产成人精品二区| 国产精品伦人一区二区| 国产色爽女视频免费观看| 99热精品在线国产| 亚洲精品影视一区二区三区av| 一区二区三区四区激情视频 | 三级毛片av免费| 国产在视频线在精品| 波野结衣二区三区在线| 欧美日韩精品成人综合77777| 99九九线精品视频在线观看视频| 三级毛片av免费| 真实男女啪啪啪动态图| 国产老妇女一区| 久久久久久久久大av| 我要搜黄色片| 亚洲成人久久爱视频| 乱码一卡2卡4卡精品| 欧美激情久久久久久爽电影| 深夜a级毛片| 搡老妇女老女人老熟妇| 桃色一区二区三区在线观看| 波多野结衣巨乳人妻| 亚洲中文字幕日韩| 22中文网久久字幕| 国产三级在线视频| 草草在线视频免费看| 听说在线观看完整版免费高清| 99久久精品国产国产毛片| 亚洲欧美精品综合久久99| av专区在线播放| 三级男女做爰猛烈吃奶摸视频| 亚洲自拍偷在线| 亚洲熟妇中文字幕五十中出| 免费搜索国产男女视频| 乱系列少妇在线播放| 欧美潮喷喷水| 亚洲成人av在线免费| 精品乱码久久久久久99久播| 久久99热6这里只有精品| 亚洲欧美日韩东京热| 91麻豆精品激情在线观看国产| 熟女人妻精品中文字幕| 亚洲av二区三区四区| 亚洲18禁久久av| 亚洲欧美中文字幕日韩二区| 网址你懂的国产日韩在线| 国产欧美日韩精品亚洲av| 麻豆乱淫一区二区| 国产精品亚洲一级av第二区| 久久精品国产清高在天天线| 给我免费播放毛片高清在线观看| 一个人看视频在线观看www免费| 久久久久久伊人网av| 久久久色成人| 国产白丝娇喘喷水9色精品| 国产精品人妻久久久影院| 国产麻豆成人av免费视频| 不卡视频在线观看欧美| 国产精品,欧美在线| 成人av在线播放网站| 国产成人一区二区在线| 天堂动漫精品| 久久99热6这里只有精品| 三级经典国产精品| 长腿黑丝高跟| 中出人妻视频一区二区| 成人特级黄色片久久久久久久| 成人特级av手机在线观看| 女的被弄到高潮叫床怎么办| 1000部很黄的大片| 少妇猛男粗大的猛烈进出视频 | 三级男女做爰猛烈吃奶摸视频| 国产av不卡久久| 久久精品国产亚洲网站| 亚洲人成网站在线播放欧美日韩| 色综合站精品国产| 老师上课跳d突然被开到最大视频| 一区二区三区高清视频在线| 欧美在线一区亚洲| 欧洲精品卡2卡3卡4卡5卡区| 国产老妇女一区| a级毛片免费高清观看在线播放| 桃色一区二区三区在线观看| 变态另类成人亚洲欧美熟女| 国产蜜桃级精品一区二区三区| 免费看美女性在线毛片视频| 日日啪夜夜撸| 搡女人真爽免费视频火全软件 | 最新中文字幕久久久久| 欧美性猛交╳xxx乱大交人| 六月丁香七月| aaaaa片日本免费| 大又大粗又爽又黄少妇毛片口| 人人妻人人澡人人爽人人夜夜 | 22中文网久久字幕| 免费不卡的大黄色大毛片视频在线观看 | 久久久久久大精品| 欧美另类亚洲清纯唯美| 亚洲精品日韩在线中文字幕 | 老司机福利观看| 亚洲欧美精品自产自拍| 成人午夜高清在线视频| 国产在线男女| 亚洲av免费在线观看| 51国产日韩欧美| 内射极品少妇av片p| 99热只有精品国产| 成年av动漫网址| 一进一出抽搐gif免费好疼| 亚洲av熟女| 淫妇啪啪啪对白视频| 国产精品久久久久久亚洲av鲁大| 久久国内精品自在自线图片| 亚洲av五月六月丁香网| 丰满乱子伦码专区| 日韩欧美国产在线观看| 3wmmmm亚洲av在线观看| 日本撒尿小便嘘嘘汇集6| 国产精品1区2区在线观看.| 国产精品三级大全| 亚洲精品乱码久久久v下载方式| 1000部很黄的大片| 欧美激情在线99| 日韩一区二区视频免费看| 看非洲黑人一级黄片| 亚洲欧美日韩高清专用| 日韩精品中文字幕看吧| 91久久精品国产一区二区三区| 天堂影院成人在线观看| 久久国产乱子免费精品| 51国产日韩欧美| 变态另类成人亚洲欧美熟女| а√天堂www在线а√下载| av专区在线播放| 国产伦精品一区二区三区视频9| 国产成人aa在线观看| 免费av不卡在线播放| 九九热线精品视视频播放| 精品人妻一区二区三区麻豆 | 好男人在线观看高清免费视频| 少妇熟女aⅴ在线视频| 精品乱码久久久久久99久播| 欧洲精品卡2卡3卡4卡5卡区| 亚洲乱码一区二区免费版| 99久久精品国产国产毛片| 久久久久国产网址| 亚洲无线在线观看| 亚洲在线观看片| 色噜噜av男人的天堂激情| 国产麻豆成人av免费视频| 亚洲成人av在线免费| 校园人妻丝袜中文字幕| 免费av观看视频| 欧美精品国产亚洲| 国产又黄又爽又无遮挡在线| www.色视频.com| 欧美丝袜亚洲另类| 老熟妇仑乱视频hdxx| 久久午夜福利片| 中文资源天堂在线| 久久精品综合一区二区三区| 亚洲在线观看片| 九九久久精品国产亚洲av麻豆| 秋霞在线观看毛片| 国产三级中文精品| 国产淫片久久久久久久久| 小说图片视频综合网站| 少妇猛男粗大的猛烈进出视频 | 国产精品福利在线免费观看| 国产精品嫩草影院av在线观看| 国产精品国产高清国产av| 99九九线精品视频在线观看视频| 国产欧美日韩精品一区二区| 美女黄网站色视频| 中文在线观看免费www的网站| 久久人人精品亚洲av| 亚洲欧美中文字幕日韩二区| 亚洲一区二区三区色噜噜| 97热精品久久久久久| 欧美日韩国产亚洲二区| 久久久精品欧美日韩精品| 一级毛片aaaaaa免费看小| 久久精品国产亚洲av香蕉五月| 啦啦啦啦在线视频资源| 97碰自拍视频| 能在线免费观看的黄片| 1000部很黄的大片| av女优亚洲男人天堂| 变态另类成人亚洲欧美熟女| 狠狠狠狠99中文字幕| 午夜免费男女啪啪视频观看 | 夜夜爽天天搞| 亚洲中文日韩欧美视频| 国产一区二区亚洲精品在线观看| 可以在线观看毛片的网站| 国产精品,欧美在线| 国产91av在线免费观看| 18禁黄网站禁片免费观看直播| 午夜激情福利司机影院| 中文字幕av成人在线电影| 一本久久中文字幕| .国产精品久久| 日日摸夜夜添夜夜爱| 国产成人一区二区在线| a级毛片a级免费在线| 国产成人a区在线观看| 最新中文字幕久久久久| 亚洲国产精品久久男人天堂| 成人三级黄色视频| 亚洲精品456在线播放app| 欧美又色又爽又黄视频| 久久精品夜色国产| 91午夜精品亚洲一区二区三区| 99久国产av精品国产电影| 一个人看的www免费观看视频| 欧美日韩乱码在线| 亚洲无线观看免费| or卡值多少钱| 国产成人一区二区在线| 欧美一区二区亚洲| 搞女人的毛片| 深爱激情五月婷婷| 国产精品亚洲美女久久久| 久久精品综合一区二区三区| 看片在线看免费视频| 亚洲18禁久久av| 久久久a久久爽久久v久久| 九九久久精品国产亚洲av麻豆| 精品午夜福利视频在线观看一区| 少妇丰满av| 成年女人毛片免费观看观看9| 国产男靠女视频免费网站| 午夜爱爱视频在线播放| 日本免费一区二区三区高清不卡| 午夜激情欧美在线| 波多野结衣高清无吗| 国产人妻一区二区三区在| 国产免费一级a男人的天堂| 美女cb高潮喷水在线观看| 亚洲av二区三区四区| 国产成人a∨麻豆精品| 黄色一级大片看看| 亚洲国产精品成人综合色| 黄色一级大片看看| 国产成人a∨麻豆精品| 久久久久国产精品人妻aⅴ院| 日本免费a在线| 尾随美女入室| 免费观看人在逋| 熟妇人妻久久中文字幕3abv| 日韩欧美免费精品| 亚洲最大成人中文| 女人十人毛片免费观看3o分钟| 一个人看视频在线观看www免费| 午夜亚洲福利在线播放| 国产成人a区在线观看| 搡老岳熟女国产| 日日摸夜夜添夜夜爱| 精品久久久久久久久久久久久| 国产高潮美女av| 亚洲国产精品成人综合色| 一级毛片aaaaaa免费看小| 丰满的人妻完整版| 国产91av在线免费观看| 国产男靠女视频免费网站| 久久韩国三级中文字幕| 成年av动漫网址| 国产精品亚洲一级av第二区| 久久精品影院6| 高清毛片免费观看视频网站| 亚洲人成网站在线观看播放| 在线观看一区二区三区| 高清毛片免费观看视频网站| 免费av观看视频| 久久99热这里只有精品18| 久久久色成人| 久久久久免费精品人妻一区二区| 日韩三级伦理在线观看| 麻豆乱淫一区二区| 日韩国内少妇激情av| 国产精品国产三级国产av玫瑰| aaaaa片日本免费| 久久精品91蜜桃| av视频在线观看入口| 亚洲最大成人中文| 中国美女看黄片| 国产高清不卡午夜福利| 性色avwww在线观看| av.在线天堂| 久久九九热精品免费| 永久网站在线| 国产精品嫩草影院av在线观看| 身体一侧抽搐| 狂野欧美白嫩少妇大欣赏| 男人舔奶头视频| 久久人人爽人人片av| 国产男靠女视频免费网站| 欧美丝袜亚洲另类| 少妇丰满av| 亚洲国产精品久久男人天堂| 国产精品国产三级国产av玫瑰| 日韩大尺度精品在线看网址| a级毛片免费高清观看在线播放| 日韩在线高清观看一区二区三区| 欧美一区二区亚洲| 看免费成人av毛片| 午夜福利在线观看吧| 99视频精品全部免费 在线| 久久国内精品自在自线图片| 三级国产精品欧美在线观看| 最近的中文字幕免费完整| 人人妻,人人澡人人爽秒播| 夜夜夜夜夜久久久久| 亚洲国产精品sss在线观看| 成人亚洲精品av一区二区| 国产伦精品一区二区三区视频9| 最近最新中文字幕大全电影3| 精品乱码久久久久久99久播| 亚洲丝袜综合中文字幕| 免费看美女性在线毛片视频| 少妇高潮的动态图| 国产亚洲精品久久久com| 国产欧美日韩精品亚洲av| 亚洲在线自拍视频| 哪里可以看免费的av片| 国产成人一区二区在线| 卡戴珊不雅视频在线播放| 在线a可以看的网站| 1024手机看黄色片| 99国产极品粉嫩在线观看| 特大巨黑吊av在线直播| 日日撸夜夜添| 亚洲三级黄色毛片| 综合色av麻豆| 久久久久国产网址| 国产精品久久久久久av不卡| 身体一侧抽搐| 免费在线观看成人毛片| 久久国产乱子免费精品| 国产又黄又爽又无遮挡在线| 亚洲精品国产成人久久av| 国产精品乱码一区二三区的特点| 高清毛片免费观看视频网站| 日韩精品青青久久久久久| 最新在线观看一区二区三区| www日本黄色视频网| 国产久久久一区二区三区| 国产高清视频在线观看网站| 亚洲av中文字字幕乱码综合| 91在线精品国自产拍蜜月| 国产精品1区2区在线观看.| 日韩高清综合在线| 嫩草影视91久久| 国产单亲对白刺激| 99热网站在线观看| 99热精品在线国产| 亚洲av第一区精品v没综合| 免费av不卡在线播放| 日本欧美国产在线视频| 日产精品乱码卡一卡2卡三| 男人狂女人下面高潮的视频| 国产精品久久久久久精品电影| 日韩大尺度精品在线看网址| 成人漫画全彩无遮挡| 国产黄片美女视频| 欧美色视频一区免费| 国产精品福利在线免费观看| 亚洲图色成人| 变态另类丝袜制服| 伊人久久精品亚洲午夜| 在线观看午夜福利视频| 欧美+亚洲+日韩+国产| 91精品国产九色| 国产精品国产三级国产av玫瑰| 日日摸夜夜添夜夜添小说| av在线天堂中文字幕| 精品久久久久久久久亚洲| 直男gayav资源| 天堂动漫精品| 天堂av国产一区二区熟女人妻| 精品免费久久久久久久清纯| 中文亚洲av片在线观看爽| 国产美女午夜福利| 成人国产麻豆网| 黄色欧美视频在线观看| 日韩欧美免费精品| 亚洲第一区二区三区不卡| 一级黄片播放器| 久久这里只有精品中国| 99热精品在线国产| 99国产精品一区二区蜜桃av| av.在线天堂| 久久久久性生活片| 久久精品夜夜夜夜夜久久蜜豆| 国产大屁股一区二区在线视频| 欧美一区二区国产精品久久精品| 又黄又爽又免费观看的视频| 国产精华一区二区三区| 亚洲自拍偷在线| 国产在视频线在精品| 国产精品爽爽va在线观看网站| 婷婷亚洲欧美| 国产精品美女特级片免费视频播放器| 久久国产乱子免费精品| 国内精品久久久久精免费| 亚洲性夜色夜夜综合| 精品福利观看| 三级毛片av免费| 久久精品夜夜夜夜夜久久蜜豆| 美女被艹到高潮喷水动态| 男女之事视频高清在线观看| 精品熟女少妇av免费看| 最新在线观看一区二区三区| 久久精品国产清高在天天线| 亚洲无线观看免费| 在线观看美女被高潮喷水网站| 波野结衣二区三区在线| 好男人在线观看高清免费视频| 亚洲激情五月婷婷啪啪| 欧美中文日本在线观看视频| 久久国产乱子免费精品| 精品久久久久久久久久免费视频| 色综合站精品国产| 久久精品国产亚洲av涩爱 | 69人妻影院| 黄色配什么色好看| а√天堂www在线а√下载| 高清午夜精品一区二区三区 | 亚洲欧美中文字幕日韩二区| 亚洲乱码一区二区免费版| 亚洲第一电影网av| 九色成人免费人妻av| a级一级毛片免费在线观看| 日本a在线网址| 精品不卡国产一区二区三区| 精品欧美国产一区二区三| 国产乱人视频| 日本黄色片子视频| 最新中文字幕久久久久| 淫妇啪啪啪对白视频| 久久综合国产亚洲精品| 久久久久国产网址| 久久久色成人| 此物有八面人人有两片| 久久草成人影院| 亚洲av.av天堂| 99热精品在线国产| 国产精品女同一区二区软件| 国内精品一区二区在线观看| 一本一本综合久久| 亚洲av中文av极速乱| 国产精品一二三区在线看| 丰满的人妻完整版| 一级毛片aaaaaa免费看小| 日韩av在线大香蕉| 成人特级黄色片久久久久久久| 欧美丝袜亚洲另类| 久久精品久久久久久噜噜老黄 | 麻豆国产97在线/欧美| 亚洲在线观看片| 变态另类丝袜制服| 99热这里只有是精品50| 一a级毛片在线观看| 一级av片app| av在线观看视频网站免费| 成人av在线播放网站| 欧美最新免费一区二区三区| 精品免费久久久久久久清纯| 日韩国内少妇激情av| 我要看日韩黄色一级片| 小蜜桃在线观看免费完整版高清| 麻豆国产97在线/欧美| 欧美日韩国产亚洲二区| 国产高清有码在线观看视频| 国语自产精品视频在线第100页| 国产亚洲av嫩草精品影院| 国产高潮美女av| 在线看三级毛片| 午夜a级毛片| 国产色爽女视频免费观看| 国产一区二区在线av高清观看| 国产精品亚洲美女久久久| 日日摸夜夜添夜夜爱| 亚洲人成网站在线播放欧美日韩| 国产真实乱freesex| 日产精品乱码卡一卡2卡三| 国产v大片淫在线免费观看| 国产私拍福利视频在线观看| 成人综合一区亚洲| 嫩草影院精品99| 看十八女毛片水多多多| 高清毛片免费观看视频网站|