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

    煤礦區(qū)重構(gòu)土壤剖面水氣變化及其對溫度梯度的響應(yīng)

    2021-05-21 09:09:18陳孝楊王校剛張凌霄朱建明
    煤炭學(xué)報 2021年4期
    關(guān)鍵詞:土壤層煤矸石矸石

    陳 敏,陳孝楊,王校剛,張凌霄,張 迪,朱建明

    (1.安徽理工大學(xué) 深部煤礦采動響應(yīng)與災(zāi)害防控國家重點實驗室,安徽 淮南 232001; 2.安徽理工大學(xué) 地球與環(huán)境學(xué)院,安徽 淮南 232001)

    煤礦開采為國內(nèi)經(jīng)濟(jì)帶來巨大貢獻(xiàn)的同時,也帶來了嚴(yán)重的生態(tài)環(huán)境問題及地質(zhì)災(zāi)害[1-2]。另外,煤炭生產(chǎn)過程也會產(chǎn)生大量的固體廢棄物,其中煤矸石作為主要廢棄物,占原煤總量的15%~20%[3]。在中國,煤矸石總量已達(dá)38億t,大量煤矸石的堆積也會造成嚴(yán)重的環(huán)境問題[4-6]。盡管煤矸石資源化利用已成為一個重要的研究領(lǐng)域,但煤矸石的利用率仍未超過15%[7]。經(jīng)過長期的理論研究與實踐,利用煤矸石充填重構(gòu)土壤來恢復(fù)或改善礦區(qū)生態(tài)環(huán)境是經(jīng)濟(jì)可行的,這不僅修復(fù)了受損土地,也在一定程度上解決了煤矸石的堆積問題[8-11]。因此,煤礦區(qū)用煤矸石作為基質(zhì)進(jìn)行土壤重構(gòu),恢復(fù)生態(tài)環(huán)境和土地復(fù)墾的立地條件工程大量存在。

    但是,煤矸石的水氣熱力學(xué)性質(zhì)與土壤存在顯著差異,其飽和含水量遠(yuǎn)低于土壤,而飽和導(dǎo)水率、熱擴(kuò)散系數(shù)和導(dǎo)氣率均大于土壤[12-14],這勢必會對重構(gòu)土壤水氣熱分布狀況及其變化產(chǎn)生影響。WANG等[15]研究發(fā)現(xiàn),重構(gòu)土壤表層土壤(0~30 cm)和填充基質(zhì)層(50~80 cm)的含水量存在顯著差異。宋楊睿等[16]研究發(fā)現(xiàn),重構(gòu)土壤含水量比自然土壤低,并且由于煤矸石持水能力極弱,無法實現(xiàn)地下水對表層土壤的有效補(bǔ)給。筆者所在研究團(tuán)隊研究發(fā)現(xiàn),重構(gòu)土壤導(dǎo)氣率受覆土類型和煤矸石的共同影響,覆土土壤導(dǎo)氣率決定了重構(gòu)土壤的導(dǎo)氣率大小,煤矸石導(dǎo)氣率決定了煤矸石對重構(gòu)土壤導(dǎo)氣率的影響系數(shù)[13];矸石層的存在會影響表層土壤溫度,當(dāng)覆土厚度較薄時,表層土壤溫度的晝夜變化幅度相對較大[17],并且矸石層對土壤層溫度的影響隨著太陽輻射的增加而增強(qiáng)[14]。特別地,當(dāng)煤矸石中黃鐵礦(FeS2)含量較高時,會發(fā)生氧化反應(yīng)并伴有釋熱過程[18]。盡管氧化過程緩慢,但依舊會對重構(gòu)土壤的溫度狀況產(chǎn)生影響,導(dǎo)致土壤層與矸石層溫度的差異。鄭國強(qiáng)等[19]對覆土矸石山溫度進(jìn)行調(diào)查,發(fā)現(xiàn)不同地塊覆土矸石山垂直方向80 cm深度溫度在28.69~52.9 ℃,溫度隨著深度的增加而升高,剖面存在顯著的溫度梯度,并且隨著熱量的持續(xù)上升,地表的溫度也隨之升高。溫度作為水氣運動的主要影響因素[20-21],會對土壤水分和養(yǎng)分的有效性、土壤的呼吸過程以及微生物的活性等產(chǎn)生影響,進(jìn)而影響植物的生長發(fā)育。那么,煤矸石氧化放熱過程對重構(gòu)土壤剖面的水氣變化的影響,勢必會影響矸石山及礦區(qū)的生態(tài)修復(fù)效果。

    然而,目前國內(nèi)對重構(gòu)土壤剖面水氣變化的研究還比較欠缺,特別是在煤矸石氧化放熱對重構(gòu)土壤剖面溫度分布及水氣變化影響這方面。因此,筆者設(shè)計了一種重構(gòu)土壤水氣熱耦合運移室內(nèi)模擬裝置,分析重構(gòu)土壤剖面水氣變化,并通過底部加熱裝置模擬煤矸石氧化放熱過程,分析其對重構(gòu)土壤剖面溫度分布及水氣變化的影響,旨在為矸石山及礦區(qū)生態(tài)修復(fù)提供理論基礎(chǔ)和技術(shù)支持。

    1 材料與方法

    1.1 試驗材料

    煤矸石樣品采自安徽省淮南市潘集區(qū)某矸石山。潘集區(qū)屬于淮南典型的煤炭開采區(qū),該矸石山現(xiàn)已停止煤矸石的輸入,呈“錐形”,表層有風(fēng)化跡象。從矸石山的山底、山腰和山頂處隨機(jī)采集等質(zhì)量風(fēng)化樣品充分混合,帶回實驗室。在實驗室內(nèi),通過篩分法測定煤矸石樣品的質(zhì)地,將粒徑大于5 cm的煤矸石挑選出來,再分別用孔徑為2 mm和10 mm的篩網(wǎng)進(jìn)行篩分,并對其飽和含水量進(jìn)行測定[12],文中含水量均為體積含水量。煤矸石樣品的機(jī)械組成及飽和含水量見表1。

    表1 煤矸石樣品的機(jī)械組成及飽和含水量Table 1 Mechanical composition and water content of coal gangue %

    土壤樣品采自淮南市大通濕地公園。取樣時,去掉表層5 cm土壤,分別使用環(huán)刀和鋁盒取樣,用于測定樣品的容重與含水量;再采集5~20 cm深度的土壤若干袋,用于室內(nèi)土柱填充。在實驗室內(nèi),用烘干法測定樣品的初始容重和含水量。從土壤樣品中隨機(jī)取樣3次,風(fēng)干后去除其中的碎石、枯枝等,研磨過2 mm篩,用Rise-2006激光粒度分析儀測定土壤顆粒組成,并依據(jù)國際分級標(biāo)準(zhǔn)對土壤質(zhì)地進(jìn)行命名,其結(jié)果見表2。

    1.2 研究方法

    1.2.1試驗裝置

    筆者設(shè)計了一種重構(gòu)土壤水氣熱耦合運移室內(nèi)模擬裝置(圖1),該裝置主要由土樣箱、加熱系統(tǒng)、傳感器、進(jìn)氣系統(tǒng)和計算機(jī)等組成。土柱箱長60 cm,寬60 cm,高120 cm,在土柱箱內(nèi)部周圍放置2 cm厚的隔熱棉,減少外界溫度對重構(gòu)土壤的影響,并在裝置底部用密封膠泥密封保證裝置的密封性。先將隨機(jī)采集的煤矸石混合樣均勻地填充在土柱底部并壓實,然后按初始容重將土樣分層均勻的填入土柱內(nèi),在填充煤矸石和土樣的過程中,每20 cm放置一對溫濕度傳感器。其中0~60 cm為土壤層,60~100 cm為矸石層[14,17],20,40,60,80和100 cm深度處為監(jiān)測點,埋有溫濕度傳感器。

    表2 土壤樣品的機(jī)械組成及物理性質(zhì)Table 2 Mechanical composition and physical properties of soil

    圖1 重構(gòu)土壤水氣熱耦合運移室內(nèi)模擬裝置Fig.1 Coupling simulation device of water,air and heat for reconstruction soil

    1.2.2試驗設(shè)置

    在土柱上方灌溉60 mm的水,通過濕度傳感器監(jiān)測水分的動態(tài)變化,待水分入滲過程穩(wěn)定后,通過CO2氣瓶以3 L/h的速率通氣,改變底部的初始CO2體積分?jǐn)?shù),使用泵吸式CO2檢測儀每隔2 h測定一次監(jiān)測點的CO2體積分?jǐn)?shù),再通過底部加熱板,分別設(shè)置30,40和50 ℃[19],監(jiān)測重構(gòu)土壤溫度的變化,待溫度穩(wěn)定后持續(xù)加熱,記錄此刻剖面的穩(wěn)態(tài)溫度和含水量(Ⅰ-Ⅰ,Ⅰ-Ⅱ,Ⅰ-Ⅲ),監(jiān)測重構(gòu)土壤水分的動態(tài)變化的同時并進(jìn)行相同的通氣試驗,最后再次在土柱上方灌溉60 mm量的水,重復(fù)上一階段試驗(Ⅱ-Ⅰ,Ⅱ-Ⅱ,Ⅱ-Ⅲ),研究加熱條件下含水量對重構(gòu)土壤剖面氣體變化的影響。不同試驗階段重構(gòu)土壤剖面溫度及含水量狀況見表3。

    1.2.3水分特征參數(shù)

    土壤水分特征曲線作為土壤基本水力特性的重要參數(shù)之一,對土壤水分運移研究具有重要意義。一般土壤水分特征曲線用經(jīng)驗公式來描述,常用的有Brooks-Corey(BC)[22]和Van Genuchten(VG)[23]經(jīng)驗公式。

    BC經(jīng)驗公式:

    (1)

    式中,θ為土壤吸力h(cm)下的土壤含水量,%;θr為滯留含水量,%;θs為飽和含水量,%;hd為土壤進(jìn)氣吸力,cm;N為擬合參數(shù)。

    表3 不同試驗階段重構(gòu)土壤剖面溫度及含水量Table 3 Temperature and water content of reconstruction soil in different test stages

    VG經(jīng)驗公式:

    (2)

    式中,α,n,m為擬合參數(shù),m=1-1/n。

    重構(gòu)土壤水流模型采用VG模型,選擇無滯后效應(yīng)。樣品水分特征參數(shù)α,n,飽和含水量(θs)及殘余含水量(θr)通過VG模型擬合土壤含水量與水勢的關(guān)系進(jìn)行逆向求解獲得,主要水分特征參數(shù)數(shù)值見表4(Ks為飽和導(dǎo)水率,cm/h)。

    表4 樣品水分特征參數(shù)Table 4 Hydraulic properties of samples

    1.3 數(shù)據(jù)處理與計算

    利用Excel 2010和Origin Pro 2016軟件統(tǒng)計分析數(shù)據(jù)及作圖,SPSS 22.0軟件對數(shù)據(jù)進(jìn)行相關(guān)性分析并建立回歸方程,Hydrus-1D軟件對重構(gòu)土壤水分分布及入滲進(jìn)行數(shù)值模擬[24]。利用軟件模擬時,重構(gòu)土壤水流模型采用VG模型,選擇無滯后效應(yīng);將重構(gòu)土壤分為兩層,兩層的參數(shù)分別輸入上覆土壤和煤矸石的土壤非飽和參數(shù);上邊界選擇大氣邊界(無徑流),下邊界選擇滲漏面邊界;上邊界在開始時,瞬時灌水量設(shè)為60 mm。

    2 結(jié)果與分析

    2.1 重構(gòu)土壤剖面水分變化及其數(shù)值模擬

    2.1.1重構(gòu)土壤含水量分布及其數(shù)值模擬

    通過對重構(gòu)土壤剖面含水量的監(jiān)測,發(fā)現(xiàn)土壤層和矸石層之間含水量存在顯著差異(圖2)。矸石層含水量明顯低于土壤層,并且在層間界面(土壤層-矸石層)含水量發(fā)生突變,層間界面土壤層的含水量顯著高于煤矸石,這與自然土壤含水量分布存在明顯區(qū)別[25]。為了驗證現(xiàn)有模型能否準(zhǔn)確模擬重構(gòu)土壤的含水量分布,根據(jù)煤矸石和土壤的水分特征參數(shù),不考慮滯后作用,對重構(gòu)土壤剖面含水量分布進(jìn)行模擬,其結(jié)果如圖2所示。通過對比,重構(gòu)土壤剖面含水量的模擬值與實測值之間雖然有些差別,但總體是相近的,特別是層間界面土壤層和矸石層的含水量。這說明根據(jù)水分特征參數(shù)模擬重構(gòu)土壤含水量時,能夠準(zhǔn)確表征不同基質(zhì)層的基質(zhì)勢差異。層間界面含水量突變主要是受基質(zhì)勢差異的影響,在同一吸力下,煤矸石的含水量要顯著低于土壤[26]。因此,當(dāng)層間界面水勢達(dá)到動態(tài)平衡時,層間界面矸石層的含水量要遠(yuǎn)低于土壤層。另外,土壤層含水量的模擬值要略低于實測值,這可能是受矸石層的影響,土壤層和矸石層存在顯著的基質(zhì)勢差異,而現(xiàn)有的方程尚未根據(jù)這種特殊的層狀土壤進(jìn)行修正。所以對重構(gòu)土壤含水量進(jìn)行擬合時,與實測值略有出入,但總體上能較好的描述重構(gòu)土壤的含水量分布狀況。

    圖2 不同深度重構(gòu)土壤含水量監(jiān)測值與模擬值Fig.2 Monitored and simulated values of water content at different depths in reconstruction soil

    2.1.2重構(gòu)土壤水分入滲及數(shù)值模擬

    重構(gòu)土壤水分入滲過程可以概化為特殊的層狀土壤水分運動問題,但與均質(zhì)土壤存在差異,其層間孔隙差異明顯。為了分析現(xiàn)有的土壤水分運動方程是否適用于重構(gòu)土壤,現(xiàn)仍然應(yīng)用Richards方程來描述重構(gòu)土壤水分運動。設(shè)置60 mm的灌溉量,模擬水分入滲過程,記錄灌溉后5,10,15和20 h不同深度的含水量變化,其結(jié)果如圖3所示。在實際水分入滲過程中,0~20 h時間段內(nèi),20~40 cm土壤層含水量先增加后減少,最大變幅為10.76%。這主要是基質(zhì)勢的影響,基質(zhì)勢作為非飽和土壤水運動過程的主要驅(qū)動力[27],水分總是從基質(zhì)勢大的地方向基質(zhì)勢小的方向運動,最后達(dá)到能量的動態(tài)平衡。在灌溉的條件下,水分由上向下運動,使20~40 cm的土壤層含水量增加,而后隨著表層土壤水分的蒸發(fā),表層土壤含水量逐漸降低,20~40 cm土壤層開始對表層土壤進(jìn)行水分的補(bǔ)給,使得20~40 cm土壤層的含水量開始降低。在入滲過程中,20~40 cm土壤層的含水量有明顯的波動,而40 cm深度以下的土壤層及矸石層含水量沒有明顯變化,并且40 cm土壤層含水量要顯著高于60 cm土壤層含水量。這說明水分在入滲過程中,在20~40 cm土壤層有一個累積的過程。而在模擬過程中,20~40 cm土壤層的含水量同樣是先增加后減小,但是其變幅較小僅有3.47%,并且其含水量降幅要小于實際值。這表明重構(gòu)土壤實際的水分蒸發(fā)速率要大于模擬值,這可能是由于煤矸石的結(jié)構(gòu)孔隙與土壤存在顯著差異,其含水量偏低,毛細(xì)管作用小,不利于對表層土壤水分的補(bǔ)給[28],并且矸石層的存在會降低了重構(gòu)土壤土壤層保水和蓄水能力[29]。另外,重構(gòu)土壤水分入滲時,水分主要在20~40 cm土壤層累積,而模擬值主要是在層間界面累積,與模擬值相比,該過程存在明顯的滯后性。這主要是受層間孔隙差異的影響,當(dāng)細(xì)質(zhì)土覆蓋粗質(zhì)土?xí)r,濕潤鋒穿過層間界面時,水分會滯留在層間界面[30],并由于層間孔隙差異導(dǎo)致毛管障礙的存在,也會阻礙水分的入滲[31]。重構(gòu)土壤與自然層狀土壤相比,其層間孔隙差異更加明顯,這可能會增強(qiáng)層面界面的滯留效果,水分未達(dá)到層間界面就開始滯留,導(dǎo)致水分運動的滯后性,使得現(xiàn)有的水分運動方程不能很好地描述重構(gòu)土壤的水分入滲問題,而準(zhǔn)確模擬水分運動過程對重構(gòu)土壤水分管理具有重要意義[32]。

    圖3 不同深度重構(gòu)土壤水分再分布監(jiān)測值與模擬值Fig.3 Monitored and simulated values of water redistribution at different depths in reconstruction soil

    2.2 重構(gòu)土壤剖面氣體變化

    重構(gòu)土壤不同深度CO2體積分?jǐn)?shù)的變化趨勢基本一致。由于底部通氣,底部CO2初始體積分?jǐn)?shù)最高,其體積分?jǐn)?shù)快速下降;80 cm層CO2體積分?jǐn)?shù)先升高后降低;60 cm層CO2體積分?jǐn)?shù)呈緩慢上升趨勢;20 cm和40 cm層CO2體積分?jǐn)?shù)沒有明顯變化(圖4)。在整個監(jiān)測過程中,矸石層的CO2體積分?jǐn)?shù)擴(kuò)散速率要顯著高于土壤層。這主要是由于煤矸石的孔隙較大,在相同含水量下,其導(dǎo)氣率要顯著高于土壤,有利于氣體在矸石層中的運動[13],并且矸石層的含水量要遠(yuǎn)低于土壤層,導(dǎo)氣率隨含水量的增加而降低[33]。同時,在80 cm層監(jiān)測到了體積分?jǐn)?shù)峰,這說明CO2在擴(kuò)散過程中在80 cm層有一個緩慢的累積過程。這可能是兩者導(dǎo)氣率和層間孔隙差異的影響,首先煤矸石的導(dǎo)氣率遠(yuǎn)高于土壤,氣體更易在矸石層中擴(kuò)散,當(dāng)氣體向土壤層擴(kuò)散時,氣體擴(kuò)散受到限制,CO2開始在矸石層中累積,并且受層間孔隙差異的影響,氣體在重構(gòu)土壤擴(kuò)散過程中表現(xiàn)出了明顯的滯后性,氣體擴(kuò)散至80 cm層后有一個受阻過程,導(dǎo)致氣體在80 cm層累積,該過程與水分入滲過程相似。

    圖4 重構(gòu)土壤剖面CO2體積分?jǐn)?shù)隨時間的變化Fig.4 Variation of CO2 concentration with time in recons-truction soil profile

    2.3 底部供熱對重構(gòu)土壤溫度的影響

    當(dāng)?shù)撞繙囟葹?0 ℃時,隨著加熱時間的增加,矸石層溫度迅速升高,土壤層溫度緩慢升高,后逐漸達(dá)到平衡(圖5)。越靠近加熱板,在開始階段溫度升高越快,溫度變化越明顯,而后溫度變化逐漸平緩,溫度穩(wěn)定所需的時間越短。在加熱 6 h時,矸石層的溫度迅速升高,而土壤層卻無明顯變化,在6 h之后土壤層溫度才緩慢上升。在整個加熱過程中,20,40,60,80和100 cm層溫度分別升高了2.06,2.33,4.04,7.99和13.47 ℃,矸石層的溫度明顯高于土壤層。這主要是因為熱量通過加熱裝置由矸石層傳向土壤層,隨著距離的增加,重構(gòu)土壤溫度的波動受加熱板的影響逐漸減小。但是,筆者所在課題組在潘一礦生態(tài)修復(fù)區(qū)對重構(gòu)土壤剖面溫度監(jiān)測時發(fā)現(xiàn),即使表層土壤溫度升高至30 ℃,對80~100 cm深度的矸石層溫度無顯著影響[14],這與重構(gòu)土壤底部加熱對表層土壤溫度影響存在差異。這主要是受煤矸石熱力學(xué)性質(zhì)影響,煤矸石是一種灰色巖石,其熱容遠(yuǎn)小于土壤,但導(dǎo)熱率卻又遠(yuǎn)高于土壤,其熱擴(kuò)散系數(shù)要遠(yuǎn)大于土壤。另外,熱擴(kuò)散系數(shù)主要受含水量、礦物組成和粒徑分布等影響[34]。煤矸石的孔隙結(jié)構(gòu)較大,并且矸石層的含水量遠(yuǎn)低于土壤層,這就有利于熱量由底部向土壤層中的傳輸,從而影響表層土壤的溫度。當(dāng)溫度達(dá)到平衡時,重構(gòu)土壤存在著顯著的溫度梯度。因此,通過改變底部加熱裝置的溫度,進(jìn)一步研究溫度對重構(gòu)土壤剖面溫度梯度的影響。

    圖5 表層土壤30 ℃時剖面溫度隨時間的變化Fig.5 Change of profile temperature with time at topsoil 30 ℃

    由圖6可知,隨著溫度的升高,矸石層的穩(wěn)態(tài)溫度和升幅要遠(yuǎn)大于土壤層,并且兩者的溫差也越來越明顯。重構(gòu)土壤20,40,60,80和100 cm的穩(wěn)態(tài)溫度分別為10.12~14.04,11.86~16.98,14.99~21.76,19.37~30.55和24.66~39.89 ℃。利用SPSS軟件對矸石層和土壤層的溫度與深度之間進(jìn)行線性擬合,發(fā)現(xiàn)兩者存在顯著的線性關(guān)系(P<0.05),根據(jù)回歸方程的斜率計算出土壤層和矸石層的溫度梯度。矸石層的溫度梯度要高于土壤層,并且隨著溫度的增加,矸石層的溫度梯度增幅要大于高土壤層。隨著溫度的增加,矸石層的溫度梯度從0.25 ℃/cm增加到0.45 ℃/cm,而土壤層的溫度梯度僅從0.11 ℃/cm增加到0.19 ℃/cm。這可能是由于土壤含水量要高于矸石層,土壤的熱容隨著含水量的增加而提高[35],使得溫度保持相對穩(wěn)定,溫度梯度相對較小。通過上述研究,可以說明當(dāng)煤矸石氧化放熱時不僅會使重構(gòu)土壤內(nèi)部形成溫度梯度,而且會影響表層土壤的溫度。溫度作為土壤水分運移和氣體擴(kuò)散的主要驅(qū)動力,重構(gòu)土壤內(nèi)部溫度梯度的存在將會改變重構(gòu)土壤內(nèi)部能量,驅(qū)動水分及氣體的運動[36]。

    圖6 不同加熱溫度下重構(gòu)土壤剖面溫度梯度Fig.6 Temperature gradient of reconstruction soil profile at different heating temperatures

    2.4 溫度梯度對重構(gòu)土壤水氣運動的影響

    2.4.1溫度梯度對重構(gòu)土壤水分運動的影響

    在溫度梯度的作用下,重構(gòu)土壤水分由底部向上遷移,100 cm層含水量顯著下降,其他層含水量逐漸增加(圖7)。

    圖7 溫度梯度下重構(gòu)土壤剖面含水量的動態(tài)變化Fig.7 Dynamic change of water content in reconstruction soil profile under temperature gradient

    在溫度梯度的作用下,重構(gòu)土壤內(nèi)部基質(zhì)勢發(fā)生改變,導(dǎo)致水分由溫度較高的區(qū)域向較低的區(qū)域遷移[37]。利用SPSS軟件對含水量與時間之間進(jìn)行相關(guān)性分析并建立回歸方程,發(fā)現(xiàn)兩者存在顯著的線性關(guān)系(P<0.05)。對回歸方程進(jìn)行分析,100 cm層含水量以0.142%/h的速度下降,20,40,60和80 cm層含水量分別以0.039,0.052,0.057和0.040%/h速度增加??梢园l(fā)現(xiàn),土壤層(40 cm和60 cm層)的含水量增速要大于矸石層(80 cm層),并且層間界面(60 cm層)的含水量增速最大,這與重構(gòu)土壤氣體變化規(guī)律不一致。這可能是由于矸石層的溫度梯度高于土壤層,在溫度梯度的作用下,水分在矸石層中的運動速度要大于土壤層,這也是100 cm層含水率迅速降低的原因。另外,在加熱的條件下,加速了水分的蒸發(fā),由于煤矸石充氣孔隙豐富,水汽更易在矸石層中擴(kuò)散[38],當(dāng)水汽遇到溫度相對較低的土壤層時,水汽容易液化,導(dǎo)致水分在層間界面累積。隨著水分在層間界面的累積,其基質(zhì)勢也隨之增加,促使水分向40 cm土壤層運動,進(jìn)一步影響含水量的變化。

    2.4.2溫度梯度及含水量對重構(gòu)土壤氣體運動的影響

    與不加熱只通氣相比,溫度梯度的存在明顯改變了氣體在重構(gòu)土壤中的運動(圖4,圖8)。在溫度梯度的作用下,重構(gòu)土壤底部CO2初始體積分?jǐn)?shù)明顯升高,均達(dá)到5.0%以上;體積分?jǐn)?shù)峰出現(xiàn)的時間縮短,其峰值更高,CO2達(dá)到平衡時間縮短,并且在60 cm層也監(jiān)測到了體積分?jǐn)?shù)峰。而在未加熱條件下,底部CO2初始體積分?jǐn)?shù)僅為4.28%,而且在60 cm層未監(jiān)測到體積分?jǐn)?shù)峰。另外,不同加熱溫度條件下,氣體變化規(guī)律也有所不同。當(dāng)?shù)撞繙囟葹?0 ℃時,60 cm和80 cm層的CO2體積分?jǐn)?shù)峰值分別為1.17%和2.58%;當(dāng)?shù)撞繙囟葹?0 ℃時,60 cm的CO2體積分?jǐn)?shù)峰值為1.51%,80 cm層的峰值達(dá)到5.0%以上;當(dāng)?shù)撞考訜釡囟葹?0 ℃時,60 cm的CO2體積分?jǐn)?shù)峰值為3.15%,80 cm層的峰值達(dá)到5.0%以上。隨著底部溫度的升高,60 cm和80 cm層的峰值逐漸升高,并且當(dāng)加熱溫度大于40 ℃時,80 cm層的CO2體積分?jǐn)?shù)峰值均高于4.99%??梢园l(fā)現(xiàn),氣體容易在80 cm層及層間界面累積,并且累積效應(yīng)隨著加熱溫度的升高而升高,而土壤層CO2體積分?jǐn)?shù)沒有明顯的波動。這主要受土壤層的阻礙作用和層間孔隙差異的影響,并且在溫度梯度的作用下,溫度差會引起介質(zhì)的遷移,隨著溫度的升高會加快氣體分子的運動速度,進(jìn)而有利于氣體的擴(kuò)散[39],而且煤矸石大孔隙發(fā)育良好,氣體會優(yōu)先通過大孔隙[40],使得煤矸石的氣體擴(kuò)散系數(shù)明顯高于土壤層。

    圖8 溫度梯度下重構(gòu)土壤剖面CO2體積分?jǐn)?shù)動態(tài)變化Fig.8 Dynamic change of CO2 concentration in reconstruction soil under temperature gradient

    另外,在溫度梯度相近的情況下(表3),含水量的增加改變了重構(gòu)土壤氣體的運動規(guī)律,特別是當(dāng)?shù)撞繙囟葹?0 ℃和40 ℃時。當(dāng)?shù)撞繙囟葹?0 ℃時,含水量的增加不僅降低了重構(gòu)土壤底部的初始體積分?jǐn)?shù),其體積分?jǐn)?shù)僅為3.83%,而且降低了60 cm和80 cm層CO2的體積分?jǐn)?shù)峰值,其體積分?jǐn)?shù)分別為0.91%和1.19%;當(dāng)?shù)撞繙囟葹?0 ℃時,含水量的增加降低了60 cm和80 cm層CO2的體積分?jǐn)?shù)峰值,其體積分?jǐn)?shù)分別為1.25%和3.75%;當(dāng)?shù)撞繙囟葹?0 ℃時,60 cm的CO2體積分?jǐn)?shù)峰值略微降低,為2.92%,80 cm層的峰值在5.0%以上。另外,當(dāng)?shù)撞繙囟葹?0 ℃和40 ℃時,體積分?jǐn)?shù)峰出現(xiàn)的時間也明顯延遲,只有當(dāng)?shù)撞繙囟葹?0 ℃時,體積分?jǐn)?shù)峰出現(xiàn)時間基本一致。這是由于含水量的增加降低了氣體在重構(gòu)土壤中的擴(kuò)散速率,土壤氣體擴(kuò)散主要是在土壤中相互連接的充氣孔隙進(jìn)行,而土壤中充氣孔隙與含水量密切相關(guān),含水量的增加會使充氣孔隙減少,從而降低氣體的擴(kuò)散[41]。此外,當(dāng)?shù)撞繙囟葹?0 ℃時,雖然60 cm層的CO2體積分?jǐn)?shù)峰值略微降低,但是總體上I-III和II-III的試驗結(jié)果基本一致。通過對比,可以說明含水率對氣體擴(kuò)散的影響隨溫度的升高而減小,當(dāng)溫度較低時,含水量作為重構(gòu)土壤氣體擴(kuò)散的主要影響因素;當(dāng)溫度較高時,溫度作為氣體擴(kuò)散的主要影響因素。

    3 討 論

    在層狀土壤中,當(dāng)細(xì)質(zhì)土覆蓋粗質(zhì)土?xí)r,由于毛管障礙的存在,濕潤鋒穿過層間界面時會在層間界面滯留[42]。而在煤礦區(qū)重構(gòu)土壤中,水分入滲過程存在明顯的滯后性,在土壤層(40 cm)累積,而非層間界面(60 cm)。同樣的,氣體擴(kuò)散過程也存在明顯的滯后性,氣體在擴(kuò)散過程中在矸石層(80 cm)累積。在底部加熱的情況下,發(fā)現(xiàn)氣體在矸石層中和層間界面均有累積的過程。氣體在擴(kuò)散過程中,首先在矸石層中受阻,當(dāng)氣體穿過矸石層后,又受到層間界面的阻礙作用。這說明在重構(gòu)土壤中可能存在一個層間界面“障礙帶”,當(dāng)水分和氣體接觸到層間界面“障礙帶”時開始受阻,而非層間界面。這可能是由于矸石層孔隙較大,層間孔隙差異明顯,導(dǎo)致毛管障礙加強(qiáng)[43],形成層間界面“障礙帶”。因此,在重構(gòu)土壤中,層間界面對水氣運動過程的影響持續(xù)時間更長,滯后性更顯著,從而導(dǎo)致現(xiàn)有的水分運動方程不能很好地描述重構(gòu)土壤的水分入滲過程。另外,卜囡等[30]在層狀土入滲試驗中發(fā)現(xiàn),當(dāng)細(xì)質(zhì)土覆蓋粗質(zhì)土?xí)r,由于粗質(zhì)地土壤的斥水性會導(dǎo)致指流的形成,并且下層介質(zhì)粒徑越大越易出現(xiàn)指流,且指流行進(jìn)速度越快。汪順生等[44]在層狀土毛管水運動特性研究中發(fā)現(xiàn),層狀土層間界面存在較大的毛細(xì)屏障,毛管水穿過層間界面需要一定的時間,導(dǎo)致層狀土中毛管水上升速率明顯低于均質(zhì)土。因此,由于“障礙帶”的存在,增加了水分在土壤層中累積時間,短期時間更有利于土壤水分的保持,阻礙水分的流失。但是,這也使得重構(gòu)土壤中更易形成指流,當(dāng)水分突破“障礙帶”時,導(dǎo)致水分更易流失,并且受“障礙帶”的影響,下層毛管水上升速度更慢,不利于對土壤層水分的補(bǔ)給,使得表層土壤含水量偏低,從而影響礦區(qū)生態(tài)修復(fù)的效果。因此,在煤矸石充填復(fù)墾區(qū)域,特別是覆土厚度較薄時,應(yīng)通過少量多次灌溉來防止水分突破“障礙帶”,保障水分在土壤層中的累積以減少水分的流失。另外,由于重構(gòu)土壤層間孔隙差異顯著導(dǎo)致“障礙帶”的形成,那么在利用煤矸石進(jìn)行充填重構(gòu)時,應(yīng)先將煤矸石破碎,以減小層間孔隙的差異,防止“障礙帶”的形成,并適當(dāng)增加覆土厚度,增加“障礙帶”的深度,削弱對表層土壤的影響。

    特別地,當(dāng)煤矸石氧化放熱時,重構(gòu)土壤中會形成穩(wěn)定的溫度梯度,甚至影響表層土壤的溫度。在溫度梯度的作用下,重構(gòu)土壤水分會向上遷移,并在層間界面“障礙帶”累積。這也會縮短毛管水穿過層間界面“障礙帶”需要的時間,增加毛管水的上升速率。在這個過程中,可能會伴隨著矸石層污染物的遷移,導(dǎo)致污染物在層間界面“障礙帶”累積并影響土壤層。鄭永紅等[45]發(fā)現(xiàn)重構(gòu)土壤層間界面土壤層重金屬元素Cu,Ni,Pb和Cd 的含量要明顯高于表層土壤,表明重金屬在層間界面富集,也證實了這一點。不管是水分入滲還是向上遷移過程中,受層間孔隙差異的影響,水分易在層間界面“障礙帶”累積。水分作為污染物遷移的重要媒介,其運動過程會對污染物遷移過程產(chǎn)生重要影響,所以水分在層間界面“障礙帶”的累積必然會伴隨著污染物的累積。在重構(gòu)土壤中,當(dāng)水分入滲時,易形成指流,可能會攜帶著污染物對地下水造成污染;當(dāng)毛管水向上遷移時,毛管水穿過層間界面時可能會攜帶著污染物,對土壤層造成污染。另外,在溫度梯度的作用下,氣體易在層間界面“障礙帶”累積。那么,當(dāng)煤矸石中黃鐵礦(FeS2)含量較高時,矸石層的大孔隙結(jié)構(gòu)有利于黃鐵礦與空氣和水接觸,在微生物的催化作用下,會釋放多種酸性氣體[46],酸性氣體易在層間界面“障礙帶”累積,則會導(dǎo)致重構(gòu)土壤內(nèi)部局部酸化[47],并隨著水分運動過程遷移。因此,在利用煤矸石進(jìn)行充填重構(gòu)時,應(yīng)先測定煤矸石中有毒微量元素及黃鐵礦的含量,考慮其有毒微量元素及酸性氣體在累積遷移過程中對土壤及地下水環(huán)境的影響。

    綜上,由于層間界面“障礙帶”的存在,使重構(gòu)土壤中水分和氣體的運動規(guī)律與均質(zhì)土壤具有很大的差異,變得更加復(fù)雜。在重構(gòu)土壤中,層間界面“障礙帶”對水氣運動的阻礙作用更大,持續(xù)時間更長,但目前對層間界面“障礙帶”對水氣流障礙發(fā)展與突變機(jī)理的研究鮮有報道。而這方面的研究,可以了解水氣流在重構(gòu)土壤中的滯留空間、障礙時間、突變條件及突變流通量變化,修正現(xiàn)有土壤水氣運動方程或建立重構(gòu)土壤數(shù)值模型,這對準(zhǔn)確預(yù)報污染物在土壤中運移及土壤水分、養(yǎng)分遷移具有重要理論和實際意義[48]。另外,重構(gòu)土壤截面孔隙結(jié)構(gòu)和剖面土層結(jié)構(gòu)是人為因素引起的[49],揭示其水氣流微觀機(jī)理及其與自然土壤的差異,可以探究原因進(jìn)行有效調(diào)控。因此,課題組下一步將開展這方面的研究,這對矸石山管理及礦區(qū)生態(tài)修復(fù)具有重要意義。

    4 結(jié) 論

    (1)Hydrus-1D軟件能夠較好地模擬重構(gòu)土壤含水量分布,但受矸石層的影響,水分入滲受阻并在土壤層累積,導(dǎo)致重構(gòu)土壤水分入滲實際過程比模擬過程存在明顯的滯后性,使得現(xiàn)有的水分運動方程不能很好地描述重構(gòu)土壤的水分入滲過程。

    (2)底部加熱時,重構(gòu)土壤剖面會形成穩(wěn)定的溫度梯度,并且矸石層的溫度明顯高于土壤層。隨著距離的增加,重構(gòu)土壤溫度的波動受加熱板的影響逐漸減小。另外,矸石層的溫度梯度要高于土壤層,并且隨著溫度的增加,矸石層的溫度梯度增幅要大于土壤層。

    (3)在溫度梯度的作用下,重構(gòu)土壤剖面水分向上遷移并在層間界面累積,層間界面的含水量增速最大,為0.057%/h;氣體擴(kuò)散速率增加,在80 cm層及層間界面累積,并且累積效應(yīng)隨著加熱溫度的升高而上升,而土壤層CO2體積分?jǐn)?shù)沒有明顯的波動。溫度和含水量均會影響重構(gòu)土壤氣體的擴(kuò)散,但溫度梯度的增加會削弱含水量對氣體擴(kuò)散的影響。

    (4)由于矸石層孔隙較大,層間孔隙差異明顯,導(dǎo)致毛管障礙加強(qiáng),可能在層間界面形成“障礙帶”。當(dāng)水分和氣體接觸到層間界面“障礙帶”時開始受阻,從而導(dǎo)致水分在土壤層開始累積,而氣體在擴(kuò)散過程中,首先在矸石層中受阻,當(dāng)氣體穿過矸石層后,又在層間界面累積。

    猜你喜歡
    土壤層煤矸石矸石
    淺談煤矸石在固廢混凝土中的應(yīng)用
    礦井矸石山環(huán)境危害與防治措施分析
    山西冶金(2022年3期)2022-08-03 08:40:28
    雄安新區(qū)上游油松林土壤層物理性質(zhì)研究
    東江中下游流域森林土壤有機(jī)碳空間分布特征研究*
    礦山矸石綠色充填系統(tǒng)設(shè)計及參數(shù)研究*
    陜西煤炭(2021年6期)2021-11-22 09:12:26
    邢東礦ZC5160/30/50型復(fù)合型充填液壓支架的研究與應(yīng)用
    煤炭與化工(2021年5期)2021-07-04 02:52:12
    煤矸石的綜合利用
    上海建材(2018年2期)2018-06-26 08:50:56
    甘肅省夏河地區(qū)影響冬蟲夏草種群分布的土壤理化因子調(diào)查
    灤河典型林分枯落物層與土壤層的水文效應(yīng)
    改性煤矸石對污水中油的吸附作用
    亚洲欧美精品专区久久| 国语对白做爰xxxⅹ性视频网站| 国产伦在线观看视频一区| 婷婷色av中文字幕| 免费观看性生交大片5| 亚洲av福利一区| 三级国产精品欧美在线观看| 内射极品少妇av片p| 日韩一区二区视频免费看| 又爽又黄a免费视频| 午夜老司机福利剧场| 欧美国产精品一级二级三级 | 国产在线视频一区二区| 97精品久久久久久久久久精品| 亚洲精品久久午夜乱码| 久久影院123| 七月丁香在线播放| 美女脱内裤让男人舔精品视频| 欧美xxxx性猛交bbbb| 国国产精品蜜臀av免费| 欧美精品一区二区免费开放| 你懂的网址亚洲精品在线观看| av专区在线播放| 中文字幕亚洲精品专区| 少妇熟女欧美另类| 男人和女人高潮做爰伦理| 午夜免费鲁丝| 国产精品一区二区在线不卡| 在线观看三级黄色| 最黄视频免费看| 久久精品人妻少妇| 网址你懂的国产日韩在线| 日韩,欧美,国产一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 久久精品国产亚洲网站| 国产精品99久久99久久久不卡 | 菩萨蛮人人尽说江南好唐韦庄| 国产极品天堂在线| 精品熟女少妇av免费看| 伦精品一区二区三区| 久久久久久久亚洲中文字幕| 特大巨黑吊av在线直播| 老司机影院成人| 成人漫画全彩无遮挡| 99久久精品热视频| 亚洲欧美日韩无卡精品| 欧美三级亚洲精品| 夫妻性生交免费视频一级片| 成人一区二区视频在线观看| 狠狠精品人妻久久久久久综合| 国产精品国产三级国产专区5o| 午夜激情福利司机影院| 国产精品久久久久久精品电影小说 | 久久99精品国语久久久| 亚洲精品自拍成人| 亚洲精品aⅴ在线观看| 国产av码专区亚洲av| 日韩,欧美,国产一区二区三区| 国产精品国产av在线观看| 国产 一区 欧美 日韩| 国产一级毛片在线| 国产在视频线精品| 黑人猛操日本美女一级片| 亚洲av男天堂| 久久久久精品久久久久真实原创| av在线观看视频网站免费| 久久久午夜欧美精品| 男女边摸边吃奶| 亚洲久久久国产精品| 大又大粗又爽又黄少妇毛片口| 精品亚洲成国产av| 蜜桃久久精品国产亚洲av| 嘟嘟电影网在线观看| 简卡轻食公司| 少妇人妻久久综合中文| 多毛熟女@视频| 精品午夜福利在线看| av女优亚洲男人天堂| 91久久精品电影网| 亚洲成人av在线免费| 国产av码专区亚洲av| 身体一侧抽搐| 亚洲自偷自拍三级| 成人免费观看视频高清| 久久久久久久久久久免费av| 99热这里只有是精品50| 看非洲黑人一级黄片| 免费高清在线观看视频在线观看| 日韩人妻高清精品专区| 91aial.com中文字幕在线观看| 精品人妻偷拍中文字幕| 午夜福利在线观看免费完整高清在| 我的老师免费观看完整版| 欧美精品一区二区大全| 欧美亚洲 丝袜 人妻 在线| 日产精品乱码卡一卡2卡三| 日韩av在线免费看完整版不卡| av在线蜜桃| 国产女主播在线喷水免费视频网站| 久久国产精品大桥未久av | 尾随美女入室| 少妇高潮的动态图| 久久鲁丝午夜福利片| 国产探花极品一区二区| 亚洲性久久影院| 日韩,欧美,国产一区二区三区| 亚洲最大成人中文| 久久亚洲国产成人精品v| 国产伦理片在线播放av一区| 熟女av电影| 久久久欧美国产精品| 亚洲av中文字字幕乱码综合| 另类亚洲欧美激情| 在线观看免费视频网站a站| 精品一区二区三卡| 成人午夜精彩视频在线观看| 97在线人人人人妻| videos熟女内射| 精品亚洲乱码少妇综合久久| 联通29元200g的流量卡| 精品一品国产午夜福利视频| 人人妻人人看人人澡| 亚洲无线观看免费| 国产免费福利视频在线观看| 国产黄色免费在线视频| 看非洲黑人一级黄片| 成人国产av品久久久| 久久精品国产鲁丝片午夜精品| 内地一区二区视频在线| 丝袜喷水一区| 秋霞伦理黄片| 日韩电影二区| 小蜜桃在线观看免费完整版高清| 亚洲,欧美,日韩| 久久女婷五月综合色啪小说| 成人毛片60女人毛片免费| 制服丝袜香蕉在线| 视频区图区小说| 男女边吃奶边做爰视频| av一本久久久久| 综合色丁香网| videossex国产| 热re99久久精品国产66热6| 亚洲精品国产av蜜桃| 老熟女久久久| 成人美女网站在线观看视频| 免费看不卡的av| 成人亚洲精品一区在线观看 | 成人国产麻豆网| 日本猛色少妇xxxxx猛交久久| 成人漫画全彩无遮挡| 男男h啪啪无遮挡| 我要看黄色一级片免费的| 最近的中文字幕免费完整| 香蕉精品网在线| 日本欧美视频一区| 人妻一区二区av| 三级国产精品欧美在线观看| 欧美bdsm另类| 高清av免费在线| 视频中文字幕在线观看| 亚洲aⅴ乱码一区二区在线播放| 国精品久久久久久国模美| 99久久精品国产国产毛片| 久久精品国产亚洲网站| 成人美女网站在线观看视频| 久久久久国产网址| 91久久精品电影网| 3wmmmm亚洲av在线观看| 亚洲精品乱久久久久久| 日本av手机在线免费观看| 亚洲精品国产av成人精品| 99热全是精品| 欧美高清性xxxxhd video| 日韩电影二区| 日韩国内少妇激情av| 1000部很黄的大片| 2022亚洲国产成人精品| 韩国高清视频一区二区三区| 极品少妇高潮喷水抽搐| 国产av一区二区精品久久 | 纵有疾风起免费观看全集完整版| 91久久精品电影网| 久久久久视频综合| 天美传媒精品一区二区| 激情五月婷婷亚洲| 六月丁香七月| 国产精品伦人一区二区| 一级爰片在线观看| 赤兔流量卡办理| 五月开心婷婷网| 国产高潮美女av| 视频区图区小说| 多毛熟女@视频| 婷婷色综合www| 免费看av在线观看网站| 夜夜爽夜夜爽视频| 国产精品.久久久| 这个男人来自地球电影免费观看 | 一区二区三区四区激情视频| 老司机影院成人| 观看免费一级毛片| 国产在线男女| 久久久久国产网址| 交换朋友夫妻互换小说| 亚洲成人中文字幕在线播放| 黑人高潮一二区| 卡戴珊不雅视频在线播放| 国产精品一二三区在线看| 久久久a久久爽久久v久久| 国产精品不卡视频一区二区| 肉色欧美久久久久久久蜜桃| 亚洲天堂av无毛| 国产亚洲91精品色在线| 欧美xxxx黑人xx丫x性爽| 在线免费十八禁| 久久精品国产鲁丝片午夜精品| 亚洲性久久影院| 日韩人妻高清精品专区| 日韩制服骚丝袜av| 欧美日韩视频高清一区二区三区二| 2021少妇久久久久久久久久久| 欧美极品一区二区三区四区| 久久亚洲国产成人精品v| 男的添女的下面高潮视频| 亚洲伊人久久精品综合| 国产高清国产精品国产三级 | 国产有黄有色有爽视频| 欧美精品一区二区免费开放| 国产精品欧美亚洲77777| 国内精品宾馆在线| av黄色大香蕉| 亚洲性久久影院| 这个男人来自地球电影免费观看 | 国语对白做爰xxxⅹ性视频网站| 免费黄网站久久成人精品| 97超碰精品成人国产| 肉色欧美久久久久久久蜜桃| 精品久久久久久电影网| 人人妻人人添人人爽欧美一区卜 | 亚洲性久久影院| 国产黄色免费在线视频| 国国产精品蜜臀av免费| 亚洲国产精品成人久久小说| 尤物成人国产欧美一区二区三区| 日本黄色日本黄色录像| 22中文网久久字幕| 日日撸夜夜添| 美女内射精品一级片tv| 一级av片app| 99re6热这里在线精品视频| 亚洲成人中文字幕在线播放| 欧美日韩视频高清一区二区三区二| 久久久久性生活片| 欧美高清成人免费视频www| 人妻夜夜爽99麻豆av| 我的女老师完整版在线观看| 久久99热这里只频精品6学生| 日本黄色片子视频| 赤兔流量卡办理| 午夜福利高清视频| 免费久久久久久久精品成人欧美视频 | 久久精品熟女亚洲av麻豆精品| 性高湖久久久久久久久免费观看| 男的添女的下面高潮视频| av免费在线看不卡| 大片免费播放器 马上看| 亚洲人成网站在线播| 国产在线视频一区二区| 午夜日本视频在线| 亚洲自偷自拍三级| 能在线免费看毛片的网站| 欧美高清性xxxxhd video| 女人久久www免费人成看片| 久久99精品国语久久久| 一级毛片电影观看| 中文字幕精品免费在线观看视频 | 国产成人精品一,二区| 三级经典国产精品| 精品久久久久久久久亚洲| 久久久精品94久久精品| 我要看黄色一级片免费的| 亚洲av在线观看美女高潮| 岛国毛片在线播放| 欧美成人午夜免费资源| 99久久精品一区二区三区| 日韩av在线免费看完整版不卡| 亚洲伊人久久精品综合| 女人十人毛片免费观看3o分钟| 欧美成人精品欧美一级黄| 国产精品一区二区三区四区免费观看| freevideosex欧美| 激情五月婷婷亚洲| a 毛片基地| 国产在线男女| 两个人的视频大全免费| 亚州av有码| 午夜福利影视在线免费观看| 国产 一区 欧美 日韩| 夜夜看夜夜爽夜夜摸| 久久久久视频综合| 亚洲精品日本国产第一区| 日韩一区二区视频免费看| 国产又色又爽无遮挡免| 交换朋友夫妻互换小说| 午夜激情久久久久久久| 夜夜骑夜夜射夜夜干| 熟妇人妻不卡中文字幕| 日韩在线高清观看一区二区三区| 欧美三级亚洲精品| 高清av免费在线| 中文乱码字字幕精品一区二区三区| 国产欧美日韩精品一区二区| 高清黄色对白视频在线免费看 | 最近最新中文字幕免费大全7| 欧美成人a在线观看| 免费黄网站久久成人精品| .国产精品久久| 99久久精品热视频| 嫩草影院入口| 涩涩av久久男人的天堂| .国产精品久久| 日本vs欧美在线观看视频 | 天天躁日日操中文字幕| 纯流量卡能插随身wifi吗| 欧美人与善性xxx| av播播在线观看一区| 舔av片在线| 欧美丝袜亚洲另类| 亚洲在久久综合| 哪个播放器可以免费观看大片| 一二三四中文在线观看免费高清| 亚洲av不卡在线观看| 日本黄色日本黄色录像| 成人国产麻豆网| 王馨瑶露胸无遮挡在线观看| 高清在线视频一区二区三区| 高清午夜精品一区二区三区| 精品久久久久久久久av| 欧美xxxx性猛交bbbb| 熟女av电影| 欧美日韩精品成人综合77777| 在线观看免费高清a一片| 亚洲不卡免费看| 身体一侧抽搐| 成人影院久久| 国产一级毛片在线| 一个人看的www免费观看视频| 一区二区av电影网| 亚洲综合精品二区| 免费少妇av软件| 亚洲精品aⅴ在线观看| 夫妻午夜视频| 国产无遮挡羞羞视频在线观看| 亚洲成色77777| 国产男人的电影天堂91| 色婷婷久久久亚洲欧美| 一个人免费看片子| 国产高潮美女av| 青春草亚洲视频在线观看| 纯流量卡能插随身wifi吗| h日本视频在线播放| 成人国产av品久久久| 插逼视频在线观看| 日韩三级伦理在线观看| 亚洲欧美一区二区三区国产| 新久久久久国产一级毛片| 女性生殖器流出的白浆| 日韩一本色道免费dvd| 国内揄拍国产精品人妻在线| 午夜福利影视在线免费观看| 日韩欧美一区视频在线观看 | 交换朋友夫妻互换小说| 国产高清国产精品国产三级 | 男女啪啪激烈高潮av片| 亚洲精品乱久久久久久| 国产人妻一区二区三区在| 另类亚洲欧美激情| 少妇人妻一区二区三区视频| 亚洲高清免费不卡视频| 97在线人人人人妻| 大话2 男鬼变身卡| 91久久精品国产一区二区三区| 国产精品伦人一区二区| 一个人看的www免费观看视频| 在线观看三级黄色| 水蜜桃什么品种好| 亚洲性久久影院| www.色视频.com| videossex国产| 青春草视频在线免费观看| 久久99精品国语久久久| 色哟哟·www| 欧美xxxx性猛交bbbb| 黄色配什么色好看| 国产伦精品一区二区三区视频9| 一二三四中文在线观看免费高清| 成人美女网站在线观看视频| 欧美bdsm另类| 天堂俺去俺来也www色官网| 久久久久久久久久人人人人人人| 国产精品女同一区二区软件| 久久国产精品男人的天堂亚洲 | 人妻一区二区av| 欧美日韩在线观看h| 欧美日韩国产mv在线观看视频 | 久久久久久伊人网av| 观看美女的网站| 九九爱精品视频在线观看| 免费观看av网站的网址| 久久人人爽人人爽人人片va| 少妇 在线观看| 日日摸夜夜添夜夜添av毛片| 日韩中字成人| 美女cb高潮喷水在线观看| 免费看不卡的av| 日韩精品有码人妻一区| 久久久久网色| 99国产精品免费福利视频| 亚洲av男天堂| 51国产日韩欧美| 日韩在线高清观看一区二区三区| 成人国产麻豆网| 少妇熟女欧美另类| 国产在线视频一区二区| 一本—道久久a久久精品蜜桃钙片| 欧美 日韩 精品 国产| 中文天堂在线官网| 97超视频在线观看视频| 一级av片app| 欧美日韩一区二区视频在线观看视频在线| 国产亚洲欧美精品永久| 免费看av在线观看网站| 你懂的网址亚洲精品在线观看| 一级二级三级毛片免费看| 在线观看人妻少妇| 久久99蜜桃精品久久| 日韩亚洲欧美综合| 黄色配什么色好看| 国国产精品蜜臀av免费| 久久精品国产亚洲网站| h日本视频在线播放| 国产精品欧美亚洲77777| av在线app专区| 在线看a的网站| 久久 成人 亚洲| 欧美一级a爱片免费观看看| 亚洲精品日韩av片在线观看| 在现免费观看毛片| 国产在线视频一区二区| 亚洲国产欧美在线一区| 狂野欧美激情性bbbbbb| 精品国产一区二区三区久久久樱花 | 午夜福利在线在线| 性色avwww在线观看| 黄色配什么色好看| 成人影院久久| 午夜免费观看性视频| 欧美日韩精品成人综合77777| 国产美女午夜福利| 国产熟女欧美一区二区| 成人漫画全彩无遮挡| 亚洲人成网站在线观看播放| 亚洲精品乱久久久久久| 成年av动漫网址| 午夜福利高清视频| 久久国内精品自在自线图片| 日韩成人伦理影院| 久久久国产一区二区| videos熟女内射| 国产一区二区在线观看日韩| 日日摸夜夜添夜夜爱| 欧美亚洲 丝袜 人妻 在线| 26uuu在线亚洲综合色| 我要看日韩黄色一级片| 色视频在线一区二区三区| 免费人成在线观看视频色| 国产高清三级在线| 男人爽女人下面视频在线观看| 午夜视频国产福利| 国产 一区精品| 王馨瑶露胸无遮挡在线观看| 中国美白少妇内射xxxbb| 久久精品久久久久久噜噜老黄| 久久精品久久久久久久性| 插逼视频在线观看| 成人18禁高潮啪啪吃奶动态图 | 亚洲美女搞黄在线观看| 汤姆久久久久久久影院中文字幕| 毛片一级片免费看久久久久| 联通29元200g的流量卡| 这个男人来自地球电影免费观看 | 日本wwww免费看| 国模一区二区三区四区视频| 国产大屁股一区二区在线视频| 精品熟女少妇av免费看| 啦啦啦在线观看免费高清www| 韩国av在线不卡| 久久久久久久久久成人| 有码 亚洲区| 欧美高清性xxxxhd video| 亚洲av中文av极速乱| 高清不卡的av网站| 久久久久久久亚洲中文字幕| 亚洲一区二区三区欧美精品| 一二三四中文在线观看免费高清| 午夜免费男女啪啪视频观看| 亚洲国产av新网站| 欧美3d第一页| 一级毛片黄色毛片免费观看视频| 国产亚洲最大av| 国产色爽女视频免费观看| 午夜日本视频在线| 热99国产精品久久久久久7| 人妻 亚洲 视频| 久久久久精品久久久久真实原创| 亚洲精品国产av成人精品| h日本视频在线播放| 免费久久久久久久精品成人欧美视频 | 久久精品夜色国产| 国产女主播在线喷水免费视频网站| 99九九线精品视频在线观看视频| 啦啦啦啦在线视频资源| 精品人妻偷拍中文字幕| 久久午夜福利片| 日韩,欧美,国产一区二区三区| 免费久久久久久久精品成人欧美视频 | 亚洲精品亚洲一区二区| 自拍欧美九色日韩亚洲蝌蚪91 | 国产综合精华液| 成人二区视频| 五月开心婷婷网| 男女免费视频国产| 99九九线精品视频在线观看视频| 男女无遮挡免费网站观看| 精品国产乱码久久久久久小说| 国产精品伦人一区二区| 国产爱豆传媒在线观看| 欧美精品亚洲一区二区| 国产精品久久久久久久电影| 久久韩国三级中文字幕| 国产精品人妻久久久久久| 亚洲国产欧美人成| 在线亚洲精品国产二区图片欧美 | 午夜福利在线在线| 免费观看a级毛片全部| 日韩免费高清中文字幕av| 久久久欧美国产精品| 欧美亚洲 丝袜 人妻 在线| 久久综合国产亚洲精品| 精品亚洲乱码少妇综合久久| 一区二区av电影网| 国产精品一二三区在线看| 日韩一本色道免费dvd| av卡一久久| 人人妻人人爽人人添夜夜欢视频 | 国产精品嫩草影院av在线观看| 18+在线观看网站| 久久久精品94久久精品| av线在线观看网站| 久久99蜜桃精品久久| 久久久色成人| 亚洲精品aⅴ在线观看| 欧美成人a在线观看| av天堂中文字幕网| 免费久久久久久久精品成人欧美视频 | 我要看黄色一级片免费的| av在线老鸭窝| 国产白丝娇喘喷水9色精品| 国产精品三级大全| 秋霞在线观看毛片| 久久久亚洲精品成人影院| 午夜老司机福利剧场| 欧美最新免费一区二区三区| 嫩草影院新地址| 精品酒店卫生间| 99九九线精品视频在线观看视频| 91精品一卡2卡3卡4卡| 菩萨蛮人人尽说江南好唐韦庄| 香蕉精品网在线| 中文字幕av成人在线电影| 久久精品夜色国产| 久久99精品国语久久久| 久久久久网色| 欧美精品人与动牲交sv欧美| 最后的刺客免费高清国语| 99视频精品全部免费 在线| 成人影院久久| 国产精品不卡视频一区二区| 亚洲国产欧美人成| 国产成人午夜福利电影在线观看| 国产成人免费观看mmmm| 中文在线观看免费www的网站| 免费观看在线日韩| 久久 成人 亚洲| 黄色日韩在线| 欧美一区二区亚洲| 国产成人一区二区在线| 草草在线视频免费看| 狂野欧美白嫩少妇大欣赏| 我的老师免费观看完整版| 国产人妻一区二区三区在| 国产精品一区二区性色av| 欧美成人a在线观看| 日韩成人伦理影院| 亚洲精品国产色婷婷电影| 又爽又黄a免费视频| 爱豆传媒免费全集在线观看| 国产黄色视频一区二区在线观看| 亚洲精品亚洲一区二区| 最后的刺客免费高清国语| 人妻系列 视频| 国产精品偷伦视频观看了| 久久久久国产精品人妻一区二区| 亚洲精品,欧美精品|