• <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)研究
    在线观看www视频免费| 精品国产乱码久久久久久小说| 成人av一区二区三区在线看 | a在线观看视频网站| 日韩一区二区三区影片| 亚洲精品av麻豆狂野| 亚洲精华国产精华精| 我要看黄色一级片免费的| 最近最新免费中文字幕在线| 欧美日韩av久久| 在线永久观看黄色视频| 黄色片一级片一级黄色片| 国产视频一区二区在线看| 老司机影院毛片| 亚洲欧美精品综合一区二区三区| bbb黄色大片| 精品国产一区二区三区四区第35| 男女国产视频网站| 在线精品无人区一区二区三| 多毛熟女@视频| 精品人妻1区二区| av网站在线播放免费| 亚洲国产欧美日韩在线播放| 美女大奶头黄色视频| 成人国产一区最新在线观看| 日本猛色少妇xxxxx猛交久久| 国产男女超爽视频在线观看| 一级a爱视频在线免费观看| 最近最新免费中文字幕在线| 亚洲 国产 在线| 国产区一区二久久| 久久国产精品男人的天堂亚洲| 汤姆久久久久久久影院中文字幕| av免费在线观看网站| 亚洲国产欧美一区二区综合| 亚洲精品国产精品久久久不卡| 色视频在线一区二区三区| 久久影院123| 国产成人av激情在线播放| 精品国产乱子伦一区二区三区 | 亚洲国产欧美网| 国产99久久九九免费精品| 夜夜夜夜夜久久久久| 久久人妻福利社区极品人妻图片| 日韩视频在线欧美| 欧美日韩精品网址| netflix在线观看网站| 欧美一级毛片孕妇| 菩萨蛮人人尽说江南好唐韦庄| 桃花免费在线播放| 久久久精品国产亚洲av高清涩受| 夜夜夜夜夜久久久久| 日韩一卡2卡3卡4卡2021年| 精品国产一区二区久久| 亚洲天堂av无毛| 可以免费在线观看a视频的电影网站| 精品久久蜜臀av无| 欧美中文综合在线视频| 多毛熟女@视频| 欧美日韩亚洲国产一区二区在线观看 | 国产视频一区二区在线看| 真人做人爱边吃奶动态| av网站在线播放免费| 黄色片一级片一级黄色片| 一级片免费观看大全| 成年人免费黄色播放视频| 他把我摸到了高潮在线观看 | 人人妻人人爽人人添夜夜欢视频| 日本a在线网址| 亚洲五月婷婷丁香| 亚洲熟女精品中文字幕| 精品一区二区三区四区五区乱码| 成人黄色视频免费在线看| 1024视频免费在线观看| cao死你这个sao货| 亚洲色图 男人天堂 中文字幕| 午夜精品久久久久久毛片777| 国产av又大| 一区在线观看完整版| 女警被强在线播放| 亚洲精品乱久久久久久| 亚洲欧美成人综合另类久久久| 久久女婷五月综合色啪小说| 国产一区二区三区综合在线观看| 女性生殖器流出的白浆| 成年美女黄网站色视频大全免费| 后天国语完整版免费观看| 国产欧美日韩综合在线一区二区| 亚洲精品国产一区二区精华液| 99久久99久久久精品蜜桃| 中文字幕制服av| 久久av网站| 亚洲欧美精品自产自拍| 欧美在线一区亚洲| 视频区图区小说| 久9热在线精品视频| 国产成人欧美在线观看 | 精品视频人人做人人爽| 巨乳人妻的诱惑在线观看| 老熟女久久久| 悠悠久久av| 天堂俺去俺来也www色官网| 老汉色∧v一级毛片| 午夜视频精品福利| 午夜精品国产一区二区电影| 97精品久久久久久久久久精品| 亚洲精品一卡2卡三卡4卡5卡 | 日韩制服丝袜自拍偷拍| 亚洲全国av大片| 国产精品国产三级国产专区5o| www日本在线高清视频| 色播在线永久视频| 久久国产亚洲av麻豆专区| 中文精品一卡2卡3卡4更新| 老司机午夜十八禁免费视频| 成人黄色视频免费在线看| 国产极品粉嫩免费观看在线| 美女大奶头黄色视频| cao死你这个sao货| 亚洲av日韩在线播放| av有码第一页| 99re6热这里在线精品视频| 亚洲久久久国产精品| 黄片小视频在线播放| 午夜免费观看性视频| 精品人妻在线不人妻| 欧美午夜高清在线| 国产一区二区激情短视频 | 国产97色在线日韩免费| 国产区一区二久久| 国产精品 欧美亚洲| 少妇猛男粗大的猛烈进出视频| 搡老熟女国产l中国老女人| 下体分泌物呈黄色| 亚洲精品成人av观看孕妇| 国产成人欧美| 99国产极品粉嫩在线观看| 一区二区三区激情视频| 亚洲欧美色中文字幕在线| 国产真人三级小视频在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲少妇的诱惑av| 久久久久网色| 亚洲欧美精品综合一区二区三区| 九色亚洲精品在线播放| 99久久综合免费| 精品久久久精品久久久| 国产成+人综合+亚洲专区| 午夜日韩欧美国产| 人妻人人澡人人爽人人| av电影中文网址| 在线十欧美十亚洲十日本专区| 久久久久久久精品精品| 免费在线观看视频国产中文字幕亚洲 | 窝窝影院91人妻| netflix在线观看网站| 少妇猛男粗大的猛烈进出视频| 真人做人爱边吃奶动态| 欧美大码av| 精品一区在线观看国产| 正在播放国产对白刺激| 色综合欧美亚洲国产小说| 久久久久国产一级毛片高清牌| 国产国语露脸激情在线看| 国产精品秋霞免费鲁丝片| 人人妻人人澡人人看| 操美女的视频在线观看| 日本撒尿小便嘘嘘汇集6| 国产欧美日韩一区二区精品| 精品福利永久在线观看| 最近最新中文字幕大全免费视频| 久久性视频一级片| 久久久精品94久久精品| 热re99久久国产66热| 每晚都被弄得嗷嗷叫到高潮| 欧美变态另类bdsm刘玥| 国产精品一区二区在线观看99| 香蕉国产在线看| 一级片'在线观看视频| 精品少妇久久久久久888优播| 女人爽到高潮嗷嗷叫在线视频| 亚洲欧美精品综合一区二区三区| 久久久久久人人人人人| 免费在线观看完整版高清| 亚洲性夜色夜夜综合| 久久免费观看电影| 色婷婷av一区二区三区视频| 亚洲 欧美一区二区三区| netflix在线观看网站| 日本av手机在线免费观看| 日日爽夜夜爽网站| 一级毛片精品| 97人妻天天添夜夜摸| 亚洲午夜精品一区,二区,三区| 欧美性长视频在线观看| 午夜91福利影院| 欧美人与性动交α欧美软件| 亚洲自偷自拍图片 自拍| 国产成人a∨麻豆精品| 国精品久久久久久国模美| 国产日韩一区二区三区精品不卡| 成人18禁高潮啪啪吃奶动态图| 搡老乐熟女国产| 日本av手机在线免费观看| 天堂8中文在线网| 国产精品一区二区免费欧美 | 国产亚洲精品一区二区www | 国产不卡av网站在线观看| 亚洲中文av在线| 一区二区三区精品91| 日韩制服骚丝袜av| 午夜老司机福利片| 美女扒开内裤让男人捅视频| 色视频在线一区二区三区| 热re99久久国产66热| 免费在线观看视频国产中文字幕亚洲 | 曰老女人黄片| 欧美激情久久久久久爽电影 | av视频免费观看在线观看| 91精品三级在线观看| 19禁男女啪啪无遮挡网站| 国产日韩欧美在线精品| 久久久久久久大尺度免费视频| 十八禁高潮呻吟视频| 一个人免费看片子| 精品人妻熟女毛片av久久网站| 国产有黄有色有爽视频| 母亲3免费完整高清在线观看| 久久人妻福利社区极品人妻图片| 最近最新免费中文字幕在线| 国产成人影院久久av| 国产av又大| 黄色片一级片一级黄色片| 日本wwww免费看| 一区二区av电影网| 精品久久久久久电影网| 最新在线观看一区二区三区| 成人国产av品久久久| 好男人电影高清在线观看| 美女中出高潮动态图| 国产野战对白在线观看| 久久精品国产亚洲av香蕉五月 | 高潮久久久久久久久久久不卡| 亚洲欧美日韩高清在线视频 | 成人影院久久| 热re99久久国产66热| a 毛片基地| 999久久久精品免费观看国产| 亚洲av成人不卡在线观看播放网 | 免费高清在线观看视频在线观看| 欧美成狂野欧美在线观看| 免费不卡黄色视频| 一进一出抽搐动态| 国产国语露脸激情在线看| 亚洲成国产人片在线观看| 亚洲精品国产色婷婷电影| 国产精品 欧美亚洲| 中国国产av一级| 精品少妇内射三级| 中国国产av一级| 亚洲国产中文字幕在线视频| 搡老乐熟女国产| 欧美中文综合在线视频| 精品高清国产在线一区| 精品久久久精品久久久| 欧美日韩黄片免| 亚洲第一av免费看| 久久久久精品人妻al黑| 亚洲一码二码三码区别大吗| 99精品久久久久人妻精品| 少妇的丰满在线观看| 正在播放国产对白刺激| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产欧美日韩在线播放| 国产伦理片在线播放av一区| 国产免费现黄频在线看| 久久久久久久久久久久大奶| 久久精品久久久久久噜噜老黄| 一级片'在线观看视频| 久久青草综合色| 亚洲国产精品成人久久小说| 每晚都被弄得嗷嗷叫到高潮| 精品国产一区二区三区四区第35| 啦啦啦 在线观看视频| 欧美在线一区亚洲| 欧美日韩福利视频一区二区| 狂野欧美激情性xxxx| 国产一区二区 视频在线| 国产亚洲精品一区二区www | 精品人妻熟女毛片av久久网站| 国产av精品麻豆| 亚洲国产欧美一区二区综合| 三级毛片av免费| 男女高潮啪啪啪动态图| 丰满少妇做爰视频| 国产视频一区二区在线看| 亚洲欧美色中文字幕在线| 欧美激情极品国产一区二区三区| av网站免费在线观看视频| 欧美黄色淫秽网站| 精品第一国产精品| 手机成人av网站| 国产成人系列免费观看| 一本一本久久a久久精品综合妖精| 永久免费av网站大全| 三级毛片av免费| 亚洲精品国产色婷婷电影| 午夜福利影视在线免费观看| 国产精品久久久久久精品古装| 免费黄频网站在线观看国产| 久热这里只有精品99| 国产一区二区三区综合在线观看| 男女免费视频国产| 免费观看a级毛片全部| 成人影院久久| 99国产精品一区二区三区| 丰满饥渴人妻一区二区三| 黑人巨大精品欧美一区二区蜜桃| 亚洲精品美女久久av网站| 91国产中文字幕| 日本av手机在线免费观看| 夜夜骑夜夜射夜夜干| 涩涩av久久男人的天堂| 十八禁网站免费在线| 久久免费观看电影| 亚洲午夜精品一区,二区,三区| 无遮挡黄片免费观看| 天天添夜夜摸| 久久久久久亚洲精品国产蜜桃av| 亚洲精品久久午夜乱码| www日本在线高清视频| 在线观看一区二区三区激情| 久久亚洲国产成人精品v| 成年人黄色毛片网站| 国产1区2区3区精品| 在线亚洲精品国产二区图片欧美| 日本黄色日本黄色录像| 亚洲天堂av无毛| 亚洲av电影在线观看一区二区三区| 国产免费视频播放在线视频| 国产亚洲一区二区精品| 波多野结衣av一区二区av| 岛国毛片在线播放| 久久天躁狠狠躁夜夜2o2o| 日本wwww免费看| 丰满人妻熟妇乱又伦精品不卡| 99国产精品99久久久久| 久久久久久免费高清国产稀缺| 精品视频人人做人人爽| 一个人免费在线观看的高清视频 | 国产成人av教育| 我的亚洲天堂| 精品国产超薄肉色丝袜足j| 久久精品aⅴ一区二区三区四区| 俄罗斯特黄特色一大片| 别揉我奶头~嗯~啊~动态视频 | 日本欧美视频一区| 国产三级黄色录像| 汤姆久久久久久久影院中文字幕| 2018国产大陆天天弄谢| 免费一级毛片在线播放高清视频 | 狠狠狠狠99中文字幕| 视频区图区小说| xxxhd国产人妻xxx| 日韩欧美国产一区二区入口| 国产精品.久久久| 一级片'在线观看视频| 亚洲中文字幕日韩| 黄片小视频在线播放| 亚洲美女黄色视频免费看| 久久ye,这里只有精品| 欧美日韩精品网址| 国产又爽黄色视频| 亚洲国产欧美一区二区综合| 国产人伦9x9x在线观看| 亚洲专区国产一区二区| 久久亚洲国产成人精品v| 日韩人妻精品一区2区三区| 免费不卡黄色视频| 中文字幕高清在线视频| 999久久久精品免费观看国产| 欧美变态另类bdsm刘玥| 亚洲av电影在线进入| 日本av免费视频播放| 黑人猛操日本美女一级片| 日韩一卡2卡3卡4卡2021年| 日韩免费高清中文字幕av| 不卡一级毛片| 久久中文看片网| 老司机午夜十八禁免费视频| 天堂8中文在线网| 色婷婷久久久亚洲欧美| 每晚都被弄得嗷嗷叫到高潮| 一区福利在线观看| 亚洲精品粉嫩美女一区| 日本av免费视频播放| 久久综合国产亚洲精品| 亚洲欧美日韩另类电影网站| 宅男免费午夜| 亚洲一区二区三区欧美精品| 男女国产视频网站| 天天躁夜夜躁狠狠躁躁| 最近最新免费中文字幕在线| 人人澡人人妻人| 亚洲专区中文字幕在线| 黄频高清免费视频| 男女国产视频网站| 精品人妻1区二区| 男女之事视频高清在线观看| 老司机福利观看| 一边摸一边做爽爽视频免费| 久久亚洲国产成人精品v| 亚洲熟女精品中文字幕| 美女福利国产在线| 亚洲av成人不卡在线观看播放网 | 国产一区二区在线观看av| 亚洲欧美一区二区三区久久| 精品人妻一区二区三区麻豆| 亚洲av日韩精品久久久久久密| 国产麻豆69| 一本综合久久免费| 亚洲欧美色中文字幕在线| a级毛片在线看网站| 久久午夜综合久久蜜桃| 12—13女人毛片做爰片一| 久久中文字幕一级| 大码成人一级视频| 超碰成人久久| 久久久久久久国产电影| 免费黄频网站在线观看国产| 欧美激情高清一区二区三区| 欧美日韩一级在线毛片| 在线观看免费午夜福利视频| 成人18禁高潮啪啪吃奶动态图| 国产亚洲精品久久久久5区| 91精品三级在线观看| 一本—道久久a久久精品蜜桃钙片| 婷婷色av中文字幕| 午夜福利免费观看在线| 首页视频小说图片口味搜索| 一级黄色大片毛片| 亚洲av日韩在线播放| 亚洲精品av麻豆狂野| 青春草视频在线免费观看| 欧美精品一区二区免费开放| 一本久久精品| 国产精品国产av在线观看| 99热国产这里只有精品6| 亚洲欧美清纯卡通| 欧美老熟妇乱子伦牲交| 国产精品欧美亚洲77777| 十八禁高潮呻吟视频| 丝袜人妻中文字幕| 欧美精品啪啪一区二区三区 | 免费观看人在逋| 丁香六月欧美| 久久狼人影院| 亚洲熟女精品中文字幕| 50天的宝宝边吃奶边哭怎么回事| 国产亚洲精品第一综合不卡| 欧美亚洲 丝袜 人妻 在线| 亚洲欧美日韩高清在线视频 | 丁香六月天网| 嫁个100分男人电影在线观看| 免费一级毛片在线播放高清视频 | 夜夜夜夜夜久久久久| 老熟妇乱子伦视频在线观看 | 在线看a的网站| 国产高清国产精品国产三级| 飞空精品影院首页| 国产男女超爽视频在线观看| 色综合欧美亚洲国产小说| 亚洲成人免费电影在线观看| www.av在线官网国产| 日韩大片免费观看网站| 97精品久久久久久久久久精品| 日韩三级视频一区二区三区| 日韩中文字幕视频在线看片| 在线亚洲精品国产二区图片欧美| 中国美女看黄片| 免费日韩欧美在线观看| 国产欧美亚洲国产| svipshipincom国产片| 成年av动漫网址| 欧美精品人与动牲交sv欧美| 精品少妇内射三级| 日韩欧美一区二区三区在线观看 | 亚洲午夜精品一区,二区,三区| 久久精品人人爽人人爽视色| 2018国产大陆天天弄谢| 999久久久国产精品视频| 欧美亚洲 丝袜 人妻 在线| 超碰成人久久| 大香蕉久久网| 久久影院123| 国产精品影院久久| 免费高清在线观看视频在线观看| 久久久久网色| av有码第一页| 青春草亚洲视频在线观看| www.自偷自拍.com| 欧美日韩成人在线一区二区| 久久精品亚洲熟妇少妇任你| 性高湖久久久久久久久免费观看| 欧美精品av麻豆av| 精品一区二区三卡| 日韩中文字幕视频在线看片| 亚洲av成人不卡在线观看播放网 | 久久精品国产综合久久久| 午夜福利视频在线观看免费| 操出白浆在线播放| 欧美日韩精品网址| 国产精品一区二区在线观看99| 国产成人啪精品午夜网站| 伊人久久大香线蕉亚洲五| 色精品久久人妻99蜜桃| 色播在线永久视频| 亚洲成人国产一区在线观看| 国产精品亚洲av一区麻豆| 不卡一级毛片| 国产免费福利视频在线观看| 黄色怎么调成土黄色| 国产又爽黄色视频| 一区在线观看完整版| 女性被躁到高潮视频| 深夜精品福利| av一本久久久久| 精品国产一区二区三区四区第35| 欧美日韩福利视频一区二区| av免费在线观看网站| 99国产综合亚洲精品| 亚洲情色 制服丝袜| 男人添女人高潮全过程视频| 日韩有码中文字幕| 秋霞在线观看毛片| 性色av乱码一区二区三区2| 精品国产国语对白av| 99re6热这里在线精品视频| 正在播放国产对白刺激| 午夜福利视频在线观看免费| 18禁黄网站禁片午夜丰满| 狂野欧美激情性bbbbbb| 肉色欧美久久久久久久蜜桃| 91精品国产国语对白视频| 亚洲视频免费观看视频| 国产精品亚洲av一区麻豆| 久久久久久人人人人人| 18在线观看网站| 岛国在线观看网站| 天天添夜夜摸| 免费一级毛片在线播放高清视频 | 一进一出抽搐动态| 99精品欧美一区二区三区四区| 91成年电影在线观看| a级毛片在线看网站| 一区二区三区精品91| 国产精品一区二区免费欧美 | 99久久国产精品久久久| 搡老熟女国产l中国老女人| avwww免费| 久久久久久亚洲精品国产蜜桃av| 一边摸一边抽搐一进一出视频| 亚洲国产看品久久| 亚洲av男天堂| 中文字幕精品免费在线观看视频| 国产xxxxx性猛交| 国产97色在线日韩免费| 亚洲中文字幕日韩| 精品久久久久久久毛片微露脸 | 色婷婷av一区二区三区视频| 12—13女人毛片做爰片一| 亚洲欧美一区二区三区黑人| 亚洲情色 制服丝袜| 三上悠亚av全集在线观看| 黄色视频,在线免费观看| 日本一区二区免费在线视频| 午夜久久久在线观看| a 毛片基地| 一个人免费看片子| 亚洲欧美一区二区三区黑人| 人妻人人澡人人爽人人| 亚洲精品美女久久久久99蜜臀| 各种免费的搞黄视频| 国产日韩欧美亚洲二区| 亚洲av国产av综合av卡| 色播在线永久视频| 日韩三级视频一区二区三区| 熟女少妇亚洲综合色aaa.| 欧美 日韩 精品 国产| 香蕉国产在线看| 久久精品亚洲熟妇少妇任你| 另类亚洲欧美激情| 一级毛片电影观看| 国产在线观看jvid| 18在线观看网站| 亚洲熟女毛片儿| 99精国产麻豆久久婷婷| 欧美精品av麻豆av| 久久亚洲精品不卡| 日本wwww免费看| 秋霞在线观看毛片| 精品免费久久久久久久清纯 | 国产成人免费无遮挡视频| 青春草视频在线免费观看| 首页视频小说图片口味搜索| 热re99久久精品国产66热6| 狠狠婷婷综合久久久久久88av| 国产99久久九九免费精品| 精品国产一区二区三区四区第35| 日韩有码中文字幕| 妹子高潮喷水视频| 天堂中文最新版在线下载| 亚洲国产精品一区三区| 国产精品香港三级国产av潘金莲| 如日韩欧美国产精品一区二区三区|