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

    包氣帶巖性結(jié)構(gòu)對地下水生態(tài)功能影響特征

    2022-09-21 02:33:58崔浩浩張光輝劉鵬飛王金哲田言亮
    水文地質(zhì)工程地質(zhì) 2022年5期
    關(guān)鍵詞:包氣毛細(xì)巖性

    崔浩浩,張光輝,劉鵬飛,王金哲,田言亮,王 茜

    (1.中國地質(zhì)科學(xué)院水文地質(zhì)環(huán)境地質(zhì)研究所, 河北 石家莊 050061;2.中國地質(zhì)大學(xué)(北京), 北京100083;3.自然資源部地下水科學(xué)與工程重點(diǎn)實(shí)驗(yàn)室, 河北 石家莊 050061)

    我國西北內(nèi)陸流域下游的天然綠洲區(qū),氣候干旱、少雨,天然植被生存、發(fā)育和生態(tài)維系對地下水水位埋深十分敏感[1]。地下水(潛水)直接或通過支持毛細(xì)作用向包氣帶表層輸供水分維系天然植被的生態(tài)能力或效應(yīng),稱之為地下水生態(tài)功能。隨著人類活動(dòng)對地下水開發(fā)利用程度的不斷提高,一旦地下水天然平衡狀態(tài)被打破,地下水失去生態(tài)功能,旱區(qū)天然植被將會(huì)嚴(yán)重退化,甚至出現(xiàn)天然綠洲荒漠化[2-5]。而包氣帶是連接植被與地下水的重要紐帶,其巖性結(jié)構(gòu)是影響潛水向天然植被根系層輸供水分的主要因子之一,它對地下水生態(tài)功能具有重要影響[6-7]。因此,包氣帶巖性結(jié)構(gòu)如何影響地下水生態(tài)功能,不僅是值得探討的學(xué)科前沿問題,也是西北干旱區(qū)地下水生態(tài)保護(hù)亟待解決的重要問題[8-9]。

    干旱區(qū)包氣帶對生態(tài)環(huán)境具有重要的意義[10-11],這一點(diǎn)眾多學(xué)者已經(jīng)達(dá)成共識(shí)。喬曉英等[12]分析了地下水生態(tài)功能在西北內(nèi)陸盆地的具體表現(xiàn),指出包氣帶巖性結(jié)構(gòu)不同,其陸表植被的生態(tài)效應(yīng)也不同,相對于單一的砂土結(jié)構(gòu),粉土和黏土的互層結(jié)構(gòu)更有利于植被的生長。陶正平等[13]研究了鄂爾多斯盆地風(fēng)積沙覆基巖型包氣帶結(jié)構(gòu)及特征,指出該結(jié)構(gòu)利于土壤水分的富集與保持,對維持當(dāng)?shù)氐纳鷳B(tài)系統(tǒng)具有重要作用;陳敏建等[14]指出地下水毛細(xì)上升高度是計(jì)算極限生態(tài)水位埋深的關(guān)鍵,而影響毛細(xì)水上升高度的主要因子是包氣帶巖性結(jié)構(gòu)(圖1)。在相同潛水位埋深下,包氣帶巖性結(jié)構(gòu)不同,直接影響地下水通過支持毛細(xì)作用向天然植被根系層土壤輸供水分的能力[15],導(dǎo)致陸表植被具有不同的生態(tài)響應(yīng)。前人關(guān)于包氣帶巖性結(jié)構(gòu)對地下水生態(tài)功能的影響研究取得了一定的成果,定性分析了旱區(qū)包氣帶巖性結(jié)構(gòu)的生態(tài)效應(yīng),但是對包氣帶不同巖性及結(jié)構(gòu)組合的研究尚不全面,也缺乏對地下水生態(tài)功能強(qiáng)弱影響的定量分析。

    圖1 旱區(qū)地下水生態(tài)功能維系天然植被生態(tài)機(jī)理簡圖Fig.1 Ecological mechanism of groundwater ecological function maintaining natural vegetation in an arid area

    針對上述問題,本文以甘肅省石羊河流域下游民勤盆地天然綠洲區(qū)為重點(diǎn)研究區(qū),依托國家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目開展的野外調(diào)查、室內(nèi)土柱試驗(yàn)和數(shù)值模擬工作,研究不同包氣帶巖性結(jié)構(gòu)與地下水耦合作用的生態(tài)效應(yīng),分析不同巖性結(jié)構(gòu)包氣帶獲取地下水供給水分和持水能力差異,定量對比不同包氣帶巖性結(jié)構(gòu)對地下水生態(tài)功能的影響特征。研究結(jié)果加深了對包氣帶在地下水生態(tài)功能中調(diào)節(jié)作用的認(rèn)識(shí),可以為旱區(qū)水資源的精細(xì)化管理及生態(tài)保護(hù)提供科學(xué)依據(jù)。

    1 研究區(qū)概況

    石羊河流域下游民勤盆地面積為1.59×104km2,其中綠洲面積僅占6%。區(qū)內(nèi)年均降水量為113.2 mm,年均蒸發(fā)量為2 675.6 mm,年均氣溫為8.8 °C,年均日照時(shí)間為3 142.2 h,是我國西北內(nèi)陸典型的干旱區(qū),被騰格里沙漠和巴丹吉林沙漠包圍,生態(tài)環(huán)境極為脆弱[16-18]。研究區(qū)位于民勤盆地農(nóng)田與騰格里沙漠之間的天然綠洲區(qū),調(diào)查點(diǎn)位布置見圖2。

    圖2 石羊河流域地形及研究區(qū)圖Fig.2 Topography of Shiyang River Basin and investigation work deployment of study area

    研究區(qū)包氣帶巖性結(jié)構(gòu)主要為“上粗下細(xì)”、“上細(xì)下粗”的多元結(jié)構(gòu)和單一“粗”結(jié)構(gòu)(圖3)。其中“上粗下細(xì)”多元結(jié)構(gòu)是:包氣帶上部地層為粗粒的粉細(xì)砂、中部為細(xì)粒的亞砂土和亞黏土、下部為粗粒的細(xì)砂地層,包括細(xì)粒土夾層結(jié)構(gòu)(研究區(qū)大部分地層為該結(jié)構(gòu),如圖3 中ZB01 鉆孔)和粗細(xì)互層結(jié)構(gòu)(如圖3中ZB12 鉆孔)?!吧霞?xì)下粗”多元結(jié)構(gòu)是:包氣帶上部地層為細(xì)粒的亞砂土或亞黏土、下部為粗粒的粉細(xì)砂或細(xì)砂地層;包括粗粒土夾層結(jié)構(gòu)或者細(xì)粗互層結(jié)構(gòu)(如圖3 中ZB13 鉆孔)。單一“粗”結(jié)構(gòu)是:包氣帶上部地層為粉細(xì)砂、下部為細(xì)砂的粗粒地層。

    從圖3 可見,單一“粗”結(jié)構(gòu)包氣帶主要分布在沙漠邊緣,遠(yuǎn)離沙漠2 km 之后,出現(xiàn)“上粗下細(xì)”的多元結(jié)構(gòu)包氣帶。在調(diào)查點(diǎn)分布方向上,自東南至西北,潛水位埋深逐漸變深;沙漠邊緣潛水位埋深約2 m,鄰近農(nóng)田區(qū)潛水位埋深增大至約為3.5 m。隨著潛水位埋深由淺至深,天然植被類型由鹽爪爪、黑枸杞和蘆葦,過渡為鹽爪爪、黑枸杞、檉柳和白刺。

    2 研究方法

    2.1 包氣帶巖性結(jié)構(gòu)調(diào)查

    在國家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目“我國西部特殊地貌區(qū)地下水開發(fā)利用與生態(tài)功能保護(hù)”成果—石羊河流域地下水功能評價(jià)與區(qū)劃的基礎(chǔ)上[19-20],選擇荒漠—綠洲過渡帶的天然植被區(qū)布設(shè)2 條包氣帶巖性結(jié)構(gòu)詳查路線,每條路線長度約10 km。采用原位探測方法,每個(gè)調(diào)查點(diǎn)自地表至潛水水面,進(jìn)行垂向探測,垂向間隔為20 cm,包括采集土壤樣和記錄包氣帶巖性結(jié)構(gòu),調(diào)查點(diǎn)位分布見圖2,調(diào)查結(jié)果見圖3。然后,以典型點(diǎn)位包氣帶剖面為基礎(chǔ)開展室內(nèi)土柱釋水試驗(yàn),研究地下水排空條件下,土壤剖面中有效持水量繼續(xù)維持典型旱區(qū)植被的生態(tài)效應(yīng);以土柱釋水試驗(yàn)為原型應(yīng)用Hydrus-1D 建立非飽和帶水分運(yùn)移數(shù)值模型,適當(dāng)調(diào)整參數(shù)并開展了模型的識(shí)別和驗(yàn)證,利用校正識(shí)別后的模型模擬評估不同包氣帶巖性結(jié)構(gòu)與地下水耦合作用對陸表植被的生態(tài)效應(yīng),進(jìn)而分析包氣帶巖性結(jié)構(gòu)對地下水生態(tài)功能的影響。

    圖3 包氣帶巖性結(jié)構(gòu)分布特征剖面圖Fig.3 Profile of lithologic structure distribution in the vadose zone

    2.2 室內(nèi)土柱釋水試驗(yàn)

    根據(jù)包氣帶巖性結(jié)構(gòu)調(diào)查結(jié)果,選擇巖性全面、結(jié)構(gòu)完整的點(diǎn)位為代表采集土樣,進(jìn)行室內(nèi)土柱釋水試驗(yàn)。土柱剖面巖性結(jié)構(gòu)是根據(jù)野外探查結(jié)果適當(dāng)調(diào)整確定,自上而下包氣帶地層巖性為亞砂土(40 cm)、亞黏土(40 cm)、粉細(xì)砂(20 cm)、細(xì)砂(50 cm),底部為25 cm 厚的粗砂作為含水層(圖4),按照原干容重每隔5 cm 進(jìn)行回填,層與層間打毛使得接觸更緊密。試驗(yàn)土柱高度200 cm,內(nèi)徑20 cm,自距離土柱頂部5 cm 位置開始,向下每隔10 cm 布設(shè)一個(gè)取樣點(diǎn),共20 個(gè),取樣點(diǎn)孔徑為2 cm,用橡膠塞封堵。試驗(yàn)土柱中心安裝Trime 土壤含水率測管,用于實(shí)時(shí)監(jiān)測土柱不同深度土壤含水率變化,測點(diǎn)間距10 cm。試驗(yàn)結(jié)束后,在各取樣點(diǎn)進(jìn)行取樣,并通過烘干法和Trime 測試各測點(diǎn)深度的土壤含水率,利用烘干法和Trime 測試的土壤含水率擬合關(guān)系,對試驗(yàn)過程中Trime 實(shí)時(shí)測試的土壤含水率進(jìn)行校正。在毛細(xì)水上升過程中,由馬氏瓶向土柱定水頭供水;在釋水過程中,馬氏瓶中的水放空,充當(dāng)排水裝置。

    圖4 驗(yàn)證試驗(yàn)(土柱)裝置示意圖Fig.4 Soil column test equipment

    2.3 數(shù)值模型

    根據(jù)天然植被區(qū)包氣帶巖性結(jié)構(gòu)的原位調(diào)查結(jié)果,確定多元結(jié)構(gòu)包氣帶從上到下巖性為“粉細(xì)砂、亞砂土、亞黏土和細(xì)砂”,單一結(jié)構(gòu)包氣帶巖性為“粉細(xì)砂和細(xì)砂”。求解土壤水分運(yùn)移模型為Hydrus-1D 模型基礎(chǔ)上增加源匯項(xiàng)的Richards 方程[21]:

    式中:θ—土壤體積含水率/(cm3·cm-3);

    t—時(shí)間/d;

    z—垂向空間坐標(biāo),向上為正;

    φ—壓力水頭/cm;

    K(θ )—水力傳導(dǎo)系數(shù)/(cm·d-1);

    S—根系吸水速率/(cm·d-1)。

    (1)初始條件及邊界條件

    本研究模擬中,暫不考慮根系吸水,由此S為0。在潛水支持毛細(xì)水上升的模擬中,設(shè)定初始壓力水頭為-1.5 MPa,上邊界為定流量邊界、流量為零;下邊界為定水頭邊界,地下水水位埋深設(shè)定為5 m,模擬時(shí)間為100 d。在土柱釋水模擬中,設(shè)定土壤剖面含水率的初始條件為30%,初始地下水水位埋深為5 m,上邊界為定流量邊界、流量為零;下邊界為自由排水邊界,模擬地下水水位埋深大幅增大的情況,模擬時(shí)間為250 d。

    (2)土壤模型結(jié)構(gòu)和水力特性參數(shù)

    土壤模型結(jié)構(gòu)以土柱釋水試驗(yàn)為依據(jù)進(jìn)行識(shí)別與校正。土壤水力參數(shù)應(yīng)用Hydrus-1D 模型的內(nèi)置模塊進(jìn)行擬合確定。其中殘余含水率是利用烘干法實(shí)測數(shù)據(jù),細(xì)砂參數(shù)采用內(nèi)置土壤分類的“sand”參數(shù)設(shè)定,其它巖性是根據(jù)土壤顆粒組成(表1),通過Rosetta模塊預(yù)測確定。預(yù)測土壤水分運(yùn)動(dòng)特征參數(shù)[22],見表2。

    表1 模擬模型巖性確定依據(jù)Table 1 Lithology of the simulation model

    表2 校正識(shí)別后不同巖性土壤水力特性參數(shù)Table 2 Corrected hydraulic characteristic parameters of different lithology soils

    利用校正識(shí)別后的參數(shù),在Hydrus-1D 中分別建立包氣帶不同巖性結(jié)構(gòu)模型,模擬其毛細(xì)水上升過程的供水能力和土柱釋水過程的持水能力。在包氣帶不同巖性單一結(jié)構(gòu)模擬中,建立粉細(xì)砂、亞砂土、亞黏土和細(xì)砂4 種模型,設(shè)定包氣帶厚度均為5 m;在包氣帶不同組合結(jié)構(gòu)模擬中,建立上粗下細(xì)、上細(xì)下粗、細(xì)粒夾層和粗粒夾層4 種模型。在上粗下細(xì)結(jié)構(gòu)中,上部細(xì)砂和下部亞黏土厚度各設(shè)為2.5 m;在上細(xì)下粗結(jié)構(gòu)中,上部亞黏土和下部細(xì)砂厚度各設(shè)為2.5 m;在夾層結(jié)構(gòu)中,設(shè)置夾層厚度(亞黏土或細(xì)砂)為1 m,上部地層厚度(細(xì)砂或亞黏土)為3 m,下部地層厚度(細(xì)砂或亞黏土)為1 m。

    (3)有效持水量

    當(dāng)植物開始發(fā)生永久凋萎時(shí)的土壤含水率,稱為凋萎含水率或萎蔫點(diǎn)[23-25];當(dāng)土壤基質(zhì)勢為-1.5 MPa時(shí)[26-27],對應(yīng)的土壤含水率為凋萎含水率。從保障天然植被生態(tài)需水角度,凋萎含水率與土壤含水率之間的土壤水量為有效持水量[28]。本研究中,基于表2 中土壤水力特性參數(shù),求解各巖性土壤水分特征曲線,進(jìn)而獲得各巖性土壤凋萎含水率,由此計(jì)算和獲得各巖性土壤有效持水量。

    3 結(jié)果

    3.1 包氣帶延滯釋水的生態(tài)效應(yīng)

    選取研究區(qū)多元結(jié)構(gòu)包氣帶不同深度土壤,包括亞砂土、亞黏土、粉細(xì)砂和細(xì)砂,采樣并建立試驗(yàn)土柱。試驗(yàn)開始之前,在定水頭下通過馬氏瓶向試驗(yàn)土柱的含水層供水,促使試驗(yàn)土體能夠充分獲得含水層(潛水)支持毛細(xì)水補(bǔ)給。317.45 h 后,支持毛細(xì)上升濕潤鋒不再繼續(xù)上移,此后開始釋水試驗(yàn),試驗(yàn)結(jié)果如圖5 所示。

    圖5 釋水后包氣帶含水率減降變化特征Fig.5 Variation characteristics of water content in the vadose zone after water release

    17日21:00 開始釋水試驗(yàn),至次日8:30 細(xì)砂和粉細(xì)砂的含水率明顯變小,而亞黏土和亞砂土的含水率變化微小。釋水過程歷時(shí)156 h(至24日9:00),再次測定試驗(yàn)土柱各深度土壤含水率,并與釋水前的土壤含水量比較,結(jié)果是試驗(yàn)土柱持水量由340.2 mm減少至305.1 mm,可供天然植被吸用的有效持水量為182.54 mm[29]。張陽陽等[30]研究表明適宜生態(tài)水位下,天然植被梭梭生長季內(nèi)(4—9月)耗水量為279 mm;極限生態(tài)水位下,梭梭耗水量為131 mm。因此,在地下水水位大幅下降之后,研究區(qū)包氣帶的有效持水量仍然能夠維持天然植被梭梭等生存一個(gè)生長季,確保旱區(qū)地下水生態(tài)功能有效。

    3.2 不同包氣帶巖性結(jié)構(gòu)對地下水生態(tài)功能影響特征

    為了研究包氣帶延滯釋水的生態(tài)效應(yīng)是否具有普遍性,以土柱試驗(yàn)的包氣帶結(jié)構(gòu)為原型,通過Hydrus-1D 建立數(shù)值模型,利用實(shí)測值對模擬值進(jìn)行識(shí)別和驗(yàn)證,開展不同包氣帶巖性結(jié)構(gòu)對地下水生態(tài)功能影響研究。通過圖6 可以看出模擬值和實(shí)測值吻合比較一致,因此,可以利用校正識(shí)別后的模型參數(shù)(表2)擬合各巖性的土壤水分特征曲線(圖7),并建立不同包氣帶結(jié)構(gòu)的數(shù)值模型,研究不同包氣帶巖性結(jié)構(gòu)對地下水生態(tài)功能影響特征,包括不同巖性的單一結(jié)構(gòu)、上粗下細(xì)結(jié)構(gòu)、上細(xì)下粗結(jié)構(gòu)、以及不同質(zhì)地的夾層結(jié)構(gòu)。

    圖6 土柱剖面各監(jiān)測點(diǎn)土壤含水率模擬值與實(shí)測值Fig.6 Simulated and measured values of soil moisture content in the soil column profile at each monitoring point

    圖7 不同巖性土壤水分特征曲線Fig.7 Soil water characteristic curve of different lithologies

    (1)包氣帶不同巖性單一結(jié)構(gòu)

    植被的極限生態(tài)水位埋深為植被根系發(fā)育帶厚度和支持毛細(xì)水上升高度之和[14]。因此,在相同植被和地下水水位埋深條件下,支持毛細(xì)水上升高度越高的巖性結(jié)構(gòu),其地下水生態(tài)功能越強(qiáng)。通過對比4 種不同包氣帶巖性單一結(jié)構(gòu)的最大支持毛細(xì)水上升高度、累計(jì)獲得地下水補(bǔ)給量和土壤含水率分布(圖8)可以看出,包氣帶巖性顆粒越細(xì),其支持毛細(xì)水上升高度和速度就越大,并且土壤獲得地下水的補(bǔ)給水分越快越多,對地表植被的生長越有利。因此,從地下水支持毛細(xì)水上升高度角度,在4 種包氣帶單一巖性結(jié)構(gòu)中,地下水生態(tài)功能強(qiáng)弱為:亞黏土>亞砂土>粉細(xì)砂>細(xì)砂。

    圖8 不同巖性包氣帶獲取地下水供給水分能力差異Fig.8 Difference characteristics of groundwater supply capacity in the vadose zone with different lithologies

    在初始條件下,設(shè)定包氣帶不同巖性土壤總持水量相同。從圖9 可以看出,地下水水位大幅下降后,包氣帶不同巖性剖面持水性能不同,各巖性剖面在前20 d 內(nèi)總持水量迅速減少,然后進(jìn)行緩慢釋水,在50 d之后逐漸趨于穩(wěn)定。結(jié)合各巖性的凋萎含水率,計(jì)算了包氣帶不同巖性的土壤凋萎含水量和釋水50 d 后的持水量、有效持水量,見表3。當(dāng)包氣帶中的有效持水量越大時(shí),其地下水生態(tài)功能越強(qiáng),從表3 中可以看出,地下水水位大幅下降時(shí),不同巖性的包氣帶土壤持水量也大幅下降,但是其有效持水量依然可觀,可以繼續(xù)維持地表植被存活。以梭梭適宜生態(tài)水位埋深下耗水量279 mm 為例,各巖性的土壤有效持水量可以維持其一個(gè)生長季的存活,具有明顯的正生態(tài)效應(yīng)。雖然釋水后細(xì)粒亞黏土的土壤持水量最大,但是有效持水量卻是中等顆粒的亞砂土最大。

    表3 釋水50 d 后不同巖性包氣帶有效持水量Table 3 Effective water holding capacity in the vadose zone of different lithologies after 50 days of water release/cm

    圖9 不同巖性包氣帶持水能力差異Fig.9 Difference of water holding capacity in the vadose zone of different lithologies

    在地下水不同水位動(dòng)態(tài)條件下,包氣帶不同巖性單一結(jié)構(gòu)對地下水生態(tài)功能的影響不同。在地下水水位上升時(shí),細(xì)粒巖性毛細(xì)水上升速度和高度都優(yōu)于粗粒巖性,表現(xiàn)出較強(qiáng)的生態(tài)功能;在地下水水位大幅下降時(shí),盡管顆粒較細(xì)的亞黏土持水量最大,但是顆粒中等的亞砂土有效持水量最高,表現(xiàn)出較強(qiáng)的生態(tài)功能。

    (2)包氣帶不同組合結(jié)構(gòu)

    包氣帶常見的結(jié)構(gòu)為層狀結(jié)構(gòu),主要包括上粗下細(xì)、上細(xì)下粗、細(xì)粒夾層和粗粒夾層等結(jié)構(gòu)。為方便討論不同結(jié)構(gòu)對地下水生態(tài)功能影響特征,細(xì)粒巖性選取亞黏土,粗粒巖性選取細(xì)砂。細(xì)粒夾層和粗粒夾層為上粗下細(xì)、上細(xì)下粗結(jié)構(gòu)的特殊情況,當(dāng)支持毛細(xì)水不能穿透中間夾層時(shí),與上述結(jié)構(gòu)情況相同。因此,在夾層結(jié)構(gòu)模擬中設(shè)置夾層厚度為1 m,以便毛細(xì)水能夠穿透夾層達(dá)到上部地層。利用Hydrus-1D 設(shè)置不同的結(jié)構(gòu)模型進(jìn)行地下水毛細(xì)上升過程模擬和土壤釋水過程模擬。除了巖性結(jié)構(gòu)外,其他條件與單一巖性結(jié)構(gòu)的模擬相同。

    通過圖10 發(fā)現(xiàn),相同巖性不同結(jié)構(gòu)的包氣帶,其對地下水生態(tài)功能影響不同。上粗下細(xì)結(jié)構(gòu)的毛細(xì)水上升高度和地下水累計(jì)供水量要遠(yuǎn)大于上細(xì)下粗結(jié)構(gòu),下部細(xì)粒的亞黏土能夠通過毛細(xì)作用將地下水輸送到上部的細(xì)砂中,利于地表植被吸水,具有明顯的正生態(tài)效應(yīng)。上細(xì)下粗結(jié)構(gòu)中,雖然支持毛細(xì)水沒有達(dá)到上部亞黏土,但是下部粗粒細(xì)砂支持毛細(xì)水上升高度(200 cm)比單一巖性的細(xì)砂(185 cm)要高15 cm,表明上部細(xì)粒的亞黏土對下伏細(xì)砂產(chǎn)生了一定的影響,使得地下水生態(tài)功能得到一定程度的提升。

    圖10 不同結(jié)構(gòu)包氣帶獲取地下水供給水分能力差異Fig.10 Difference characteristics of groundwater supply capacity in the vadose zone with different structures

    從圖10 可以看出,細(xì)粒夾層(250 cm)和粗粒夾層(260 cm)的毛細(xì)水上升高度基本相同,但是地下水向二者的剖面累計(jì)供水量不同,細(xì)粒夾層的累計(jì)供水量(38.31 cm)大于粗粒夾層的累計(jì)供水量(29.75 cm)。因此,細(xì)粒夾層結(jié)構(gòu)具有較強(qiáng)的供水能力,地下水生態(tài)功能較強(qiáng)。

    在土壤釋水過程中,包氣帶不同巖性結(jié)構(gòu)的土壤持水量隨釋水時(shí)間的動(dòng)態(tài)變化與不同巖性單一結(jié)構(gòu)具有相似的特征(圖9 和圖11)。通過土壤有效持水量計(jì)算結(jié)果發(fā)現(xiàn)表4,包氣帶結(jié)構(gòu)的組合方式對土壤持水性能具有顯著影響,即使巖性完全相同,僅組合結(jié)構(gòu)不同,其有效持水量也會(huì)有較大差異。

    圖11 不同結(jié)構(gòu)包氣帶持水能力差異Fig.11 Difference in water holding capacity in the vadose zone of different structures

    表4 釋水50 d 后不同結(jié)構(gòu)包氣帶有效持水量Table 4 Effective water holding capacity in the vadose zone of different structures after 50 days of water release/cm

    當(dāng)包氣帶中的有效持水量越大時(shí),其地下水生態(tài)功能越強(qiáng)。通過上粗下細(xì)和上細(xì)下粗結(jié)構(gòu)對比,二者初始持水量和凋萎含水量相同,但是在自由釋水50 d 后,上細(xì)下粗結(jié)構(gòu)釋水量較少,有效持水量較多;粗粒夾層和細(xì)粒夾層對比,在自由釋水50 d 后,粗粒夾層有效持水量較多,表明上細(xì)下粗和粗粒夾層結(jié)構(gòu)其生態(tài)功能較強(qiáng)。綜合對比包氣帶單一巖性和組合結(jié)構(gòu),在釋水后,多種巖性的組合結(jié)構(gòu)有效持水量較大。尤其是通過對比單一亞黏土結(jié)構(gòu)和粗粒夾層,單一細(xì)砂結(jié)構(gòu)和細(xì)粒夾層,不管夾層的巖性如何,夾層的存在都使得包氣帶剖面的有效持水量增加。因此,當(dāng)?shù)叵滤宦裆畲蠓龃蠛螅喾N巖性的組合結(jié)構(gòu)更有利于持有較大的有效持水量,生態(tài)效應(yīng)更強(qiáng)。

    4 討論

    在同一研究區(qū)域,天然植被類型、地下水水位埋深和氣象影響因素基本相同,唯有包氣帶巖性結(jié)構(gòu)不同。不同的包氣帶巖性結(jié)構(gòu)對地下水生態(tài)功能的影響主要表現(xiàn)在支持毛細(xì)水高度不同和潛水位大幅下降過程中植被主根系層土壤釋水和含水量減少程度方面。理論上,相同潛水水位埋深下,地下水支持毛細(xì)水上升高度越大,地下水供給天然植被根系層土壤水分能力就越強(qiáng),地下水生態(tài)功能則越強(qiáng);在潛水水位大幅下降過程中,天然植被主根系層土壤釋水越慢、含水率減小越少,越有利于維系天然植被生態(tài),地下水生態(tài)功能則越強(qiáng)。在西北內(nèi)陸干旱區(qū),不同植被的適宜生態(tài)水位埋深均值為2.9 m,極限生態(tài)水位埋深均值為5.5 m[31]。根據(jù)前期在青土湖區(qū)開展的地下水水位埋深動(dòng)態(tài)監(jiān)測數(shù)據(jù)可知,由于農(nóng)田灌溉或者未得到及時(shí)的人工輸水,天然植被區(qū)地下水水位埋深會(huì)在一段時(shí)間內(nèi)大幅增大,甚至大于極限生態(tài)水位(圖12)。從圖12 可知,在青土湖核心區(qū),地下水水位埋深在6—8月時(shí)段大于適宜生態(tài)水位;在外圍區(qū),地下水水位埋深在8—10月時(shí)段大于極限生態(tài)水位,雖然此時(shí)植被無法通過支持毛細(xì)作用持續(xù)獲得地下水的補(bǔ)給,但是包氣帶的有效持水量會(huì)繼續(xù)維持植被存活一段時(shí)間。包氣帶有效持水量越大,其地下水生態(tài)功能越強(qiáng);待地下水水位埋深逐漸減小時(shí),支持毛細(xì)水上升速度越快越有利于對土壤水的補(bǔ)給,其地下水生態(tài)功能越強(qiáng)。

    圖12 天然植被區(qū)地下水水位埋深動(dòng)態(tài)變化規(guī)律Fig.12 Dynamic response of groundwater depth to artificial water conveyance in the natural vegetation area

    在土柱釋水試驗(yàn)中驗(yàn)證了包氣帶中有效持水量具有繼續(xù)維持典型旱區(qū)植被的生態(tài)效應(yīng);數(shù)值模擬中分析了包氣帶不同巖性和結(jié)構(gòu)與地下水耦合作用的生態(tài)效應(yīng),因此,包氣帶巖性結(jié)構(gòu)對地下水生態(tài)功能具有重要的影響。在包氣帶單一巖性結(jié)構(gòu)中,顆粒的粗細(xì)決定著地下水支持毛細(xì)水上升高度和土壤有效持水量,從而影響地下水生態(tài)功能。顆粒越細(xì),其支持毛細(xì)水上升高度越大,持水量也越大,但是由于其較高的凋萎含水率,細(xì)粒巖性的有效持水量不是最高;而中等巖性顆粒的持水量雖然不是最大,但是由于其較低的凋萎含水率,有效持水量反而最高[32]。根據(jù)土壤水分特征曲線可以看出(圖7),在獲得相同含水率情況下,巖性顆粒越細(xì),其需要的吸持力越大。因此,細(xì)粒巖性土壤對水分具有較高的吸持力,水分不容易被植被利用,而中等顆粒的亞砂土中砂粒和粉粒含量較高,顆粒較粗,雖然總持水量相對較少,但是土壤對水分的吸持力也較小,可供植被吸收利用的水分較多。而粉細(xì)砂和細(xì)砂以砂粒為主的粗顆粒含量較高,土壤顆粒大不容易吸持水分,持水性能最低。

    在包氣帶不同結(jié)構(gòu)中,不管是上粗下細(xì)結(jié)構(gòu)、上細(xì)下粗結(jié)構(gòu)還是不同質(zhì)地的夾層結(jié)構(gòu),導(dǎo)致其不同地下水生態(tài)功能差異的原因就是水分在穿過粗、細(xì)巖性界面時(shí)發(fā)生的不同吸持力的變化。當(dāng)?shù)叵滤上孪蛏线\(yùn)移時(shí),由于細(xì)粒對水分的吸持力較大,水分由細(xì)粒向粗粒運(yùn)移時(shí)速度會(huì)減緩,而水分由粗粒向細(xì)粒運(yùn)移時(shí)速度會(huì)加快;當(dāng)?shù)叵滤缮舷蛳逻\(yùn)移時(shí),由于細(xì)粒對水分的吸持力較大,以及細(xì)粒的滲透性能較差,水分由細(xì)粒向粗粒運(yùn)移或者由粗粒向細(xì)粒運(yùn)移時(shí)速度都會(huì)減緩,從而加強(qiáng)了包氣帶延滯釋水效應(yīng),使得包氣帶持有的有效水分能夠在一段時(shí)期內(nèi)繼續(xù)供給地表植被,發(fā)揮地下水的生態(tài)功能。

    需要說明的是,上述不同包氣帶不同組合結(jié)構(gòu)的模擬結(jié)果只是針對設(shè)定的結(jié)構(gòu)情景。由于不同包氣帶結(jié)構(gòu)的巖層厚度不同,夾層的巖性、厚度、位置以及層數(shù)不同,其毛細(xì)上升高度和持水量也不同[33],這些因素均會(huì)影響地下水生態(tài)功能。

    5 結(jié)論

    (1)當(dāng)?shù)叵滤宦裆钪饾u減小時(shí),在相同植被條件下,包氣帶巖性結(jié)構(gòu)不同,其支持毛細(xì)水上升高度和速度不同,對地下水生態(tài)功能的影響也不同。單一巖性中,顆粒越細(xì),其支持毛細(xì)水上升高度和速度越大,土壤獲得地下水的補(bǔ)給水分越快越多,對地表植被的生長越有利;其中,亞黏土的毛細(xì)水上升高度為5 m,地下水累計(jì)供水量為92.17 cm,表現(xiàn)出較強(qiáng)的生態(tài)效應(yīng)。不同巖性結(jié)構(gòu)組合的支持毛細(xì)水上升高度比較復(fù)雜,因巖性厚度、位置、層數(shù)等而異;其中,上粗下細(xì)結(jié)構(gòu)的毛細(xì)水上升高度為3.20 m,地下水累計(jì)供水量為54.80 cm,表現(xiàn)出較強(qiáng)的生態(tài)效應(yīng)。

    (2)當(dāng)?shù)叵滤宦裆畲蠓龃蠛?,旱區(qū)包氣帶具有繼續(xù)維持陸表植被存活的生態(tài)效應(yīng),不同巖性和結(jié)構(gòu)的包氣帶其有效持水量均可以維持旱區(qū)典型植被(梭梭)存活一個(gè)生長季(4—9月),從而拓展了“生態(tài)水位”的含義。單一巖性中,包氣帶巖性顆粒越細(xì),其持水量越大,但是其凋萎含水量也越高;其中,中等巖性顆粒的亞砂土有效持水量最大,為69.68 cm,對維持植被的生態(tài)效應(yīng)最明顯。與單一巖性相比,多種巖性組合結(jié)構(gòu)的有效持水量更大,生態(tài)效應(yīng)更強(qiáng);其中,粗粒夾層結(jié)構(gòu)有效持水量最大,為74.96 cm,上細(xì)下粗性結(jié)構(gòu)有效持水量次之,為73.57 cm,表現(xiàn)出較強(qiáng)的生態(tài)效應(yīng)。

    因此,充分認(rèn)識(shí)旱區(qū)包氣帶不同巖性和結(jié)構(gòu)在地下水水位埋深增大和減小情況下對地下水生態(tài)功能的影響,利用包氣帶在地下水生態(tài)功能中的調(diào)節(jié)作用,可以為當(dāng)?shù)厮Y源的精細(xì)化管理和生態(tài)保護(hù)提供科學(xué)依據(jù)。

    猜你喜歡
    包氣毛細(xì)巖性
    “毛細(xì)”政務(wù)號(hào)關(guān)停背后
    廉政瞭望(2024年5期)2024-05-26 13:21:07
    基于氯離子示蹤法深厚包氣帶地區(qū)地下水補(bǔ)給特征
    層狀非均質(zhì)結(jié)構(gòu)包氣帶入滲過程單相流與兩相流數(shù)值模擬對比研究
    多孔建筑材料毛細(xì)吸水系數(shù)連續(xù)測量的方法
    能源工程(2020年6期)2021-01-26 00:55:16
    一種識(shí)別薄巖性氣藏的地震反射特征分析
    云南化工(2020年11期)2021-01-14 00:50:58
    K 近鄰分類法在巖屑數(shù)字圖像巖性分析中的應(yīng)用
    錄井工程(2017年1期)2017-07-31 17:44:42
    出現(xiàn)憋喘 可能是毛細(xì)支氣管炎!
    媽媽寶寶(2017年3期)2017-02-21 01:22:16
    高滲鹽水霧化吸入治療毛細(xì)支氣管炎的療效觀察
    包氣帶淺層地?zé)崛萘坑?jì)算方法商榷
    鉛和鉻污染包氣帶及再釋放規(guī)律的實(shí)驗(yàn)研究
    欧美午夜高清在线| 制服人妻中文乱码| 美国免费a级毛片| 精品国产超薄肉色丝袜足j| 美女午夜性视频免费| 国产精品日韩av在线免费观看| 国产亚洲av嫩草精品影院| 成年人黄色毛片网站| 精品国产乱码久久久久久男人| 欧美精品亚洲一区二区| 亚洲中文字幕一区二区三区有码在线看 | 中文字幕久久专区| 午夜成年电影在线免费观看| 欧美黄色淫秽网站| 一卡2卡三卡四卡精品乱码亚洲| 国产亚洲欧美精品永久| 亚洲天堂国产精品一区在线| 中出人妻视频一区二区| 中文字幕精品亚洲无线码一区 | 色哟哟哟哟哟哟| 香蕉av资源在线| 最近最新中文字幕大全免费视频| 美女扒开内裤让男人捅视频| 少妇 在线观看| 亚洲国产看品久久| 亚洲精华国产精华精| 人妻久久中文字幕网| 欧美乱色亚洲激情| 国产一区在线观看成人免费| 最近最新免费中文字幕在线| 欧美久久黑人一区二区| 国产爱豆传媒在线观看 | 99精品久久久久人妻精品| 亚洲欧美日韩高清在线视频| 日韩免费av在线播放| 亚洲精品国产区一区二| 精品久久久久久久毛片微露脸| 制服诱惑二区| 男男h啪啪无遮挡| 身体一侧抽搐| 久久久国产精品麻豆| 中文字幕人妻熟女乱码| 99久久99久久久精品蜜桃| 好男人电影高清在线观看| АⅤ资源中文在线天堂| 免费在线观看亚洲国产| 女性被躁到高潮视频| 国产野战对白在线观看| 欧美激情高清一区二区三区| 97超级碰碰碰精品色视频在线观看| 波多野结衣巨乳人妻| 亚洲 国产 在线| 99在线人妻在线中文字幕| 欧美成狂野欧美在线观看| 亚洲人成网站高清观看| 久久精品人妻少妇| 久久久久久国产a免费观看| 中文字幕人妻熟女乱码| 国产成+人综合+亚洲专区| 人人妻,人人澡人人爽秒播| 91成年电影在线观看| 热99re8久久精品国产| 久久 成人 亚洲| 12—13女人毛片做爰片一| 亚洲国产欧洲综合997久久, | 久久欧美精品欧美久久欧美| 亚洲欧美一区二区三区黑人| 亚洲avbb在线观看| 精品不卡国产一区二区三区| 亚洲成av片中文字幕在线观看| 国产亚洲av嫩草精品影院| 50天的宝宝边吃奶边哭怎么回事| 久久久久亚洲av毛片大全| 国产精品爽爽va在线观看网站 | 色哟哟哟哟哟哟| 变态另类丝袜制服| bbb黄色大片| 淫妇啪啪啪对白视频| 亚洲午夜理论影院| 免费观看精品视频网站| 亚洲精品国产一区二区精华液| 午夜久久久久精精品| 精品国产乱子伦一区二区三区| 看片在线看免费视频| 成人欧美大片| av电影中文网址| 色哟哟哟哟哟哟| 久热这里只有精品99| 欧美日韩中文字幕国产精品一区二区三区| 欧美zozozo另类| 黄色片一级片一级黄色片| 99热这里只有精品一区 | 18禁国产床啪视频网站| 亚洲av日韩精品久久久久久密| 色精品久久人妻99蜜桃| 可以免费在线观看a视频的电影网站| 日韩欧美国产在线观看| 亚洲国产高清在线一区二区三 | 91麻豆av在线| 久热爱精品视频在线9| 亚洲精品久久成人aⅴ小说| av超薄肉色丝袜交足视频| 级片在线观看| 身体一侧抽搐| 国产亚洲精品久久久久5区| 精品欧美一区二区三区在线| 欧美不卡视频在线免费观看 | 嫁个100分男人电影在线观看| 级片在线观看| 神马国产精品三级电影在线观看 | 可以免费在线观看a视频的电影网站| 日韩精品中文字幕看吧| 中文字幕人妻丝袜一区二区| 国内精品久久久久久久电影| 又黄又爽又免费观看的视频| 国产成年人精品一区二区| 久久久国产成人免费| 麻豆一二三区av精品| 人人妻人人澡欧美一区二区| 久久久久国内视频| 色综合婷婷激情| 在线视频色国产色| 18禁裸乳无遮挡免费网站照片 | 国产高清激情床上av| 视频在线观看一区二区三区| 麻豆一二三区av精品| 无遮挡黄片免费观看| 人成视频在线观看免费观看| 精品国产乱子伦一区二区三区| 97碰自拍视频| 午夜福利在线观看吧| 人妻丰满熟妇av一区二区三区| 91字幕亚洲| 免费观看人在逋| svipshipincom国产片| 动漫黄色视频在线观看| 国内毛片毛片毛片毛片毛片| 国产片内射在线| 欧美亚洲日本最大视频资源| 亚洲av第一区精品v没综合| 国产精品永久免费网站| 不卡一级毛片| 精品免费久久久久久久清纯| 我的亚洲天堂| 97人妻精品一区二区三区麻豆 | 日本熟妇午夜| 国产欧美日韩精品亚洲av| 97超级碰碰碰精品色视频在线观看| 婷婷亚洲欧美| 女性生殖器流出的白浆| 国产久久久一区二区三区| 在线观看www视频免费| 成人免费观看视频高清| 久久久久精品国产欧美久久久| 欧美亚洲日本最大视频资源| 亚洲av熟女| 级片在线观看| 波多野结衣av一区二区av| 中亚洲国语对白在线视频| 一级毛片女人18水好多| 国产不卡一卡二| 在线观看免费日韩欧美大片| 天天一区二区日本电影三级| 欧美激情极品国产一区二区三区| 欧美成人免费av一区二区三区| 日本成人三级电影网站| 成在线人永久免费视频| 欧美中文综合在线视频| 不卡一级毛片| 大型av网站在线播放| 最新美女视频免费是黄的| 日韩免费av在线播放| 51午夜福利影视在线观看| 成人18禁高潮啪啪吃奶动态图| 老熟妇仑乱视频hdxx| 中文字幕另类日韩欧美亚洲嫩草| 欧美国产日韩亚洲一区| 亚洲精品一区av在线观看| 欧美zozozo另类| xxx96com| 亚洲国产看品久久| 人人妻人人澡人人看| 不卡av一区二区三区| 亚洲成人国产一区在线观看| 可以免费在线观看a视频的电影网站| 欧美成人午夜精品| 黄频高清免费视频| 亚洲在线自拍视频| 国产成人系列免费观看| 久99久视频精品免费| 亚洲最大成人中文| 淫秽高清视频在线观看| 国产又黄又爽又无遮挡在线| 亚洲国产精品sss在线观看| 欧美成人一区二区免费高清观看 | 亚洲人成网站高清观看| 男女下面进入的视频免费午夜 | 亚洲人成网站在线播放欧美日韩| 久久久国产精品麻豆| 日韩欧美在线二视频| 老熟妇仑乱视频hdxx| 成人午夜高清在线视频 | 香蕉丝袜av| 国产人伦9x9x在线观看| 久久国产乱子伦精品免费另类| 高清毛片免费观看视频网站| 好看av亚洲va欧美ⅴa在| 久久天堂一区二区三区四区| 精品久久蜜臀av无| 天天躁夜夜躁狠狠躁躁| 亚洲一码二码三码区别大吗| 看片在线看免费视频| 久久精品夜夜夜夜夜久久蜜豆 | 中出人妻视频一区二区| 日韩精品青青久久久久久| 中国美女看黄片| 色综合婷婷激情| 精品国产美女av久久久久小说| 国产国语露脸激情在线看| 国产黄片美女视频| 淫秽高清视频在线观看| 免费观看精品视频网站| 国产蜜桃级精品一区二区三区| 国产成人一区二区三区免费视频网站| 中文亚洲av片在线观看爽| av在线天堂中文字幕| 国产在线精品亚洲第一网站| 久久精品夜夜夜夜夜久久蜜豆 | 国产激情久久老熟女| 精品电影一区二区在线| 久久久久九九精品影院| 一本大道久久a久久精品| 亚洲久久久国产精品| 久久草成人影院| 999精品在线视频| 99国产精品一区二区蜜桃av| 精品国产国语对白av| 黄片小视频在线播放| 最新在线观看一区二区三区| 老司机靠b影院| 啪啪无遮挡十八禁网站| 亚洲一区二区三区不卡视频| 久久久久久久午夜电影| 夜夜躁狠狠躁天天躁| 18禁黄网站禁片免费观看直播| 欧美日韩一级在线毛片| 国产午夜福利久久久久久| 黄频高清免费视频| 亚洲国产日韩欧美精品在线观看 | 国产伦一二天堂av在线观看| 亚洲av日韩精品久久久久久密| 日韩欧美国产在线观看| 亚洲精品久久成人aⅴ小说| 美女 人体艺术 gogo| 国产亚洲av高清不卡| 精品国产乱码久久久久久男人| 日韩大尺度精品在线看网址| 久久中文看片网| 婷婷丁香在线五月| 亚洲真实伦在线观看| 青草久久国产| 午夜激情av网站| 国产精品精品国产色婷婷| 精品久久久久久久久久免费视频| 他把我摸到了高潮在线观看| 亚洲男人的天堂狠狠| 午夜福利在线在线| 可以在线观看的亚洲视频| 一级黄色大片毛片| 老鸭窝网址在线观看| 看黄色毛片网站| 亚洲成人国产一区在线观看| 在线播放国产精品三级| 在线观看一区二区三区| 精品久久蜜臀av无| 亚洲精品在线观看二区| 51午夜福利影视在线观看| 亚洲欧洲精品一区二区精品久久久| 天天一区二区日本电影三级| www日本黄色视频网| 久久久久久久精品吃奶| √禁漫天堂资源中文www| 国产精品亚洲一级av第二区| 老汉色av国产亚洲站长工具| 99久久国产精品久久久| 亚洲国产日韩欧美精品在线观看 | 欧美在线一区亚洲| 久久久久精品国产欧美久久久| 国产一级毛片七仙女欲春2 | 亚洲中文日韩欧美视频| 久久久久亚洲av毛片大全| 欧美zozozo另类| 精品久久久久久久久久免费视频| 亚洲精品在线观看二区| 母亲3免费完整高清在线观看| 成人一区二区视频在线观看| 国产精品一区二区三区四区久久 | 精品久久久久久久毛片微露脸| 美女扒开内裤让男人捅视频| 色哟哟哟哟哟哟| 久久人妻福利社区极品人妻图片| 男女那种视频在线观看| 精品久久久久久,| 亚洲国产精品成人综合色| 国产aⅴ精品一区二区三区波| 精品电影一区二区在线| 免费搜索国产男女视频| 国产精品一区二区免费欧美| 国产亚洲av高清不卡| 波多野结衣高清作品| 亚洲欧美激情综合另类| 淫秽高清视频在线观看| 美女扒开内裤让男人捅视频| 一区二区三区激情视频| 色综合站精品国产| 欧美一级毛片孕妇| 国产精品98久久久久久宅男小说| 中文字幕人妻丝袜一区二区| 中文字幕最新亚洲高清| 宅男免费午夜| 制服人妻中文乱码| 精品国产美女av久久久久小说| 欧美日韩亚洲国产一区二区在线观看| 国产不卡一卡二| 久久草成人影院| 人人妻人人澡欧美一区二区| 国产97色在线日韩免费| 天天躁狠狠躁夜夜躁狠狠躁| 久久精品国产99精品国产亚洲性色| 亚洲av成人不卡在线观看播放网| 无人区码免费观看不卡| 悠悠久久av| 午夜精品在线福利| 亚洲国产欧美日韩在线播放| av福利片在线| 日韩免费av在线播放| 麻豆国产av国片精品| 亚洲 欧美 日韩 在线 免费| 青草久久国产| 久久婷婷成人综合色麻豆| 久久久久国内视频| 国产精品亚洲一级av第二区| 精品久久久久久久末码| 亚洲人成网站在线播放欧美日韩| 国产精华一区二区三区| 动漫黄色视频在线观看| 欧美色视频一区免费| 夜夜看夜夜爽夜夜摸| 日韩大码丰满熟妇| 成人特级黄色片久久久久久久| 视频在线观看一区二区三区| 熟女电影av网| 99riav亚洲国产免费| 侵犯人妻中文字幕一二三四区| 69av精品久久久久久| 搞女人的毛片| cao死你这个sao货| 久热爱精品视频在线9| 久久天躁狠狠躁夜夜2o2o| 老司机在亚洲福利影院| 久久精品成人免费网站| 麻豆成人av在线观看| 1024视频免费在线观看| 日韩中文字幕欧美一区二区| 国产亚洲av嫩草精品影院| 一区二区三区激情视频| 麻豆国产av国片精品| 丝袜在线中文字幕| 国产成人av教育| 久久久久久九九精品二区国产 | 亚洲 欧美一区二区三区| 亚洲国产欧洲综合997久久, | 成在线人永久免费视频| 欧美绝顶高潮抽搐喷水| 亚洲欧美一区二区三区黑人| 女性被躁到高潮视频| 午夜福利免费观看在线| 欧美黄色淫秽网站| 99国产精品一区二区三区| 一区二区日韩欧美中文字幕| 最近在线观看免费完整版| 窝窝影院91人妻| 无遮挡黄片免费观看| 在线观看午夜福利视频| 亚洲国产日韩欧美精品在线观看 | 丁香六月欧美| 听说在线观看完整版免费高清| 亚洲国产日韩欧美精品在线观看 | 成人特级黄色片久久久久久久| 桃红色精品国产亚洲av| 一区二区三区高清视频在线| 中文字幕精品亚洲无线码一区 | av福利片在线| 亚洲真实伦在线观看| 欧美黄色淫秽网站| 一区福利在线观看| 精品第一国产精品| 黄色视频不卡| 国产亚洲精品一区二区www| 欧美激情高清一区二区三区| 白带黄色成豆腐渣| 一本久久中文字幕| 日韩免费av在线播放| 国产成人av教育| 国产男靠女视频免费网站| 久久久国产成人精品二区| 国产成人系列免费观看| 亚洲第一电影网av| 天堂动漫精品| 看片在线看免费视频| 亚洲狠狠婷婷综合久久图片| 国产伦一二天堂av在线观看| 一个人观看的视频www高清免费观看 | 欧美乱妇无乱码| 每晚都被弄得嗷嗷叫到高潮| 国产国语露脸激情在线看| 国产av在哪里看| 成年女人毛片免费观看观看9| 丁香欧美五月| 美女扒开内裤让男人捅视频| 宅男免费午夜| 中文字幕av电影在线播放| 亚洲狠狠婷婷综合久久图片| 欧美日韩瑟瑟在线播放| 麻豆av在线久日| 久久精品91蜜桃| 别揉我奶头~嗯~啊~动态视频| 大型黄色视频在线免费观看| 久久久国产成人精品二区| 一a级毛片在线观看| 亚洲欧洲精品一区二区精品久久久| 亚洲男人天堂网一区| 丁香六月欧美| 国产伦人伦偷精品视频| 少妇 在线观看| 制服诱惑二区| 嫩草影视91久久| 亚洲成国产人片在线观看| 久久精品夜夜夜夜夜久久蜜豆 | netflix在线观看网站| 精品第一国产精品| 黄片小视频在线播放| 精品国内亚洲2022精品成人| 99国产综合亚洲精品| 欧美激情久久久久久爽电影| 在线观看免费日韩欧美大片| 久久精品人妻少妇| 欧美成人免费av一区二区三区| 国产免费av片在线观看野外av| 黄色a级毛片大全视频| 成人免费观看视频高清| 色婷婷久久久亚洲欧美| 亚洲精品久久国产高清桃花| 国产高清视频在线播放一区| 性欧美人与动物交配| 久久精品国产亚洲av香蕉五月| 日本三级黄在线观看| 深夜精品福利| 热re99久久国产66热| 91av网站免费观看| 女同久久另类99精品国产91| 99精品久久久久人妻精品| 国产男靠女视频免费网站| 亚洲精品国产一区二区精华液| 黄片小视频在线播放| 波多野结衣巨乳人妻| 日本a在线网址| 精品不卡国产一区二区三区| 18禁黄网站禁片午夜丰满| 国产野战对白在线观看| 国产av一区在线观看免费| 十八禁网站免费在线| 久久久久久久久中文| 男女之事视频高清在线观看| 亚洲中文字幕日韩| 日韩一卡2卡3卡4卡2021年| 美女免费视频网站| 日本一本二区三区精品| 欧美日韩中文字幕国产精品一区二区三区| 女性被躁到高潮视频| 午夜福利视频1000在线观看| 日韩免费av在线播放| 老司机午夜福利在线观看视频| 亚洲激情在线av| 成人精品一区二区免费| cao死你这个sao货| 99久久久亚洲精品蜜臀av| e午夜精品久久久久久久| 久久精品aⅴ一区二区三区四区| 亚洲人成网站高清观看| 欧美日韩亚洲综合一区二区三区_| 777久久人妻少妇嫩草av网站| 香蕉国产在线看| 国产亚洲精品第一综合不卡| 美女大奶头视频| 最近最新中文字幕大全电影3 | 九色国产91popny在线| 91成年电影在线观看| 亚洲色图av天堂| 久久午夜亚洲精品久久| 亚洲欧美日韩无卡精品| 色尼玛亚洲综合影院| 国产一区二区激情短视频| 老司机深夜福利视频在线观看| a在线观看视频网站| 91字幕亚洲| 国产激情久久老熟女| 99在线人妻在线中文字幕| 欧美 亚洲 国产 日韩一| 在线观看午夜福利视频| 人人妻,人人澡人人爽秒播| 日本三级黄在线观看| 女生性感内裤真人,穿戴方法视频| 欧美绝顶高潮抽搐喷水| 国产色视频综合| 91成人精品电影| 91大片在线观看| 成人亚洲精品一区在线观看| 欧美精品亚洲一区二区| 国产麻豆成人av免费视频| 久久精品夜夜夜夜夜久久蜜豆 | 久久久久国产精品人妻aⅴ院| 午夜精品在线福利| 18禁黄网站禁片免费观看直播| 精品一区二区三区av网在线观看| 免费女性裸体啪啪无遮挡网站| 久久精品91蜜桃| 首页视频小说图片口味搜索| 美女高潮喷水抽搐中文字幕| 国产成+人综合+亚洲专区| 一区二区三区精品91| 免费看a级黄色片| 一本精品99久久精品77| 国内毛片毛片毛片毛片毛片| 欧美三级亚洲精品| 成年女人毛片免费观看观看9| 欧美黑人欧美精品刺激| 欧美性长视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 国产又色又爽无遮挡免费看| 50天的宝宝边吃奶边哭怎么回事| 亚洲七黄色美女视频| 欧美三级亚洲精品| 老司机午夜福利在线观看视频| 波多野结衣高清作品| 日本熟妇午夜| 别揉我奶头~嗯~啊~动态视频| 国产又色又爽无遮挡免费看| 久久 成人 亚洲| 亚洲精华国产精华精| 欧美又色又爽又黄视频| 国产人伦9x9x在线观看| 丝袜美腿诱惑在线| 成人特级黄色片久久久久久久| 性欧美人与动物交配| 亚洲国产欧洲综合997久久, | 久久午夜亚洲精品久久| 欧美在线黄色| 可以在线观看毛片的网站| 久久亚洲真实| 国产精品影院久久| 久久国产精品人妻蜜桃| 国产伦人伦偷精品视频| 午夜视频精品福利| 99热这里只有精品一区 | 国产亚洲精品第一综合不卡| 精品一区二区三区av网在线观看| av视频在线观看入口| 欧美亚洲日本最大视频资源| 熟女少妇亚洲综合色aaa.| 国产午夜精品久久久久久| 精品福利观看| 很黄的视频免费| 国产精品免费一区二区三区在线| 熟妇人妻久久中文字幕3abv| 别揉我奶头~嗯~啊~动态视频| 久热这里只有精品99| 黑人操中国人逼视频| 国产色视频综合| 午夜福利欧美成人| 波多野结衣巨乳人妻| 国产成人影院久久av| 亚洲国产精品999在线| 欧美久久黑人一区二区| 国产高清有码在线观看视频 | av欧美777| 成人精品一区二区免费| 嫁个100分男人电影在线观看| 亚洲第一av免费看| 一级毛片高清免费大全| 俺也久久电影网| 啦啦啦免费观看视频1| 波多野结衣巨乳人妻| 桃红色精品国产亚洲av| 国产一区在线观看成人免费| 热re99久久国产66热| 91国产中文字幕| 国产精品香港三级国产av潘金莲| 亚洲av片天天在线观看| 一本一本综合久久| 美女免费视频网站| 最新美女视频免费是黄的| 亚洲av日韩精品久久久久久密| 两性午夜刺激爽爽歪歪视频在线观看 | 12—13女人毛片做爰片一| av超薄肉色丝袜交足视频| 丝袜在线中文字幕| 日韩一卡2卡3卡4卡2021年| 国产一区在线观看成人免费| 亚洲天堂国产精品一区在线| 国产黄a三级三级三级人| 国产伦在线观看视频一区| 欧美成人一区二区免费高清观看 | 亚洲av电影不卡..在线观看| 免费电影在线观看免费观看|