張映梅 夏瓊 王旭

摘 要：為進一步掌握西北干旱地區(qū)非飽和黃土水分蒸發(fā)規(guī)律，進行不同環(huán)境溫度下非飽和黃土水分室內(nèi)模擬蒸發(fā)試驗，探討蒸發(fā)過程中土體含水率變化、溫度對蒸發(fā)速率及累計蒸發(fā)量的影響。結(jié)果表明：隨蒸發(fā)時間的延長，土體含水率日改變率逐漸減小;環(huán)境溫度對土體含水率的影響與土體深度范圍、蒸發(fā)時段均有關(guān)，相同的升溫值導(dǎo)致土體含水率的變化量是不相等的，不同的溫度范圍內(nèi)土體含水率敏感變化深度不同;環(huán)境溫度與室溫接近時，土體每天的蒸發(fā)速率隨時間延長沒有明顯變化，但當(dāng)環(huán)境溫度高于室溫一定值時，土體每天蒸發(fā)速率隨時間呈遞減狀態(tài);環(huán)境溫度一定時，氣候模擬箱自身的冷凝作用使得非飽和土體累計蒸發(fā)量隨蒸發(fā)時長呈波動增加趨勢。

關(guān)鍵詞：非飽和黃土：溫度;含水率;蒸發(fā)速率;累計算蒸發(fā)量

中圖分類號：TU444 文獻標(biāo)志碼：A

doi：10.3969/j.issn.1000-1379.2020.05.030

Abstract： In order to further understand the law of water evaporation of unsaturated loess in the arid area of Northwest China. Based on the indoor model test， the water evaporation experiments of unsaturated loess under different ambient temperatures were carried out. The effects of soil moisture content and temperature on evaporation rate and cumulative evaporation were discussed. The results show that with the extension of evaporation time， the diurnal change of soil moisture content is decreased gradually. The influence of environmental temperature to the soil moisture content is related with the depth range of soil sample， evaporation stage and reference temperature. The same temperature increase results in different moisture changes， and different temperature ranges correspond to different depth ranges in which the soil moisture content is changed sensitively. When the environmental temperature is close to the indoor temperature， the evaporation rate develops with time without obvious fluctuation， but the evaporation rate is fluctuated and decreased with time obviously when the environmental temperature is higher than the indoor temperature to a certain extent.

Key words： unsaturated loess; temperature; moisture content; evaporation rate; cumulative evaporation

我國黃土主要分布在弱降雨、強蒸發(fā)的干旱與半干旱地區(qū)，地下水位埋藏較深，因此該地區(qū)黃土大多處于非飽和狀態(tài)。降雨入滲和蒸發(fā)作用使得土體內(nèi)水分發(fā)生動態(tài)變化，進而使得土體結(jié)構(gòu)、強度、承載力及滲透性等發(fā)生改變，經(jīng)常引發(fā)一系列巖土工程問題。工程實踐發(fā)現(xiàn)，反復(fù)的干濕循環(huán)會引起邊坡土體開裂、地基收縮，進而導(dǎo)致邊坡失穩(wěn)、地基沉降、建筑物發(fā)生傾斜開裂等[1-2]?，F(xiàn)有的關(guān)于土體水分蒸發(fā)方面的研究側(cè)重于土壤學(xué)、農(nóng)學(xué)和氣象學(xué)等學(xué)科，但由于研究背景和目的的不同，其相關(guān)成果很難直接應(yīng)用于巖土工程領(lǐng)域。隨著“一帶一路”倡議的持續(xù)推進，西北地區(qū)的工程建設(shè)將進一步加快，因此針對西北干旱與半干旱地區(qū)黃土進行特定環(huán)境下的土體蒸發(fā)研究，可為該地區(qū)黃土土工結(jié)構(gòu)物水分變化及其引起的應(yīng)力應(yīng)變等分析和預(yù)測提供依據(jù)。

近年來，國內(nèi)學(xué)者王鐵行等[3-5]通過室內(nèi)試驗和數(shù)值模擬，研究了不同干密度、含水率和溫度影響下黃土土水特征曲線;劉小軍等[6]在等溫條件下研究了非飽和黃土中的氣態(tài)、液態(tài)水遷移規(guī)律;騰繼東等[7]通過理論推導(dǎo)和試驗分析探討了風(fēng)速影響下土體水分運移規(guī)律;唐朝生等[8]通過室內(nèi)試驗探討了土體水分蒸發(fā)運動過程的變化規(guī)律;吳珺華等[9]通過數(shù)值模擬計算，獲得了蒸發(fā)條件下膨脹土內(nèi)部水分分布特征和變化規(guī)律。但是，這些研究僅在試驗、模擬角度考慮單因素條件下土體水分遷移規(guī)律，沒有考慮晝夜溫差較大時，大氣中的水氣在溫度較低的夜晚會冷凝形成冷凝水，滴落在土體表面，從而對土體水分遷移的影響，因此很難指導(dǎo)大氣環(huán)境中土體水分運移規(guī)律研究。此外，近些年來在西北黃土高原旱區(qū)綠化和邊坡土體植被防護，以及工程防護中，與土體水分遷移及其有關(guān)的巖土工程問題比較嚴峻[10]。為進一步掌握西北干旱地區(qū)非飽和黃土水分蒸發(fā)規(guī)律，開展大氣作用下非飽和黃土水分蒸發(fā)試驗很有必要。鑒于大氣環(huán)境中影響因素較多，筆者采用控制變量法研究環(huán)境溫度對土體水分蒸發(fā)的影響，以期為黃土地區(qū)蒸發(fā)引起的工程病害和機制分析提供理論支撐。

1 試驗材料與方法

1.1 試驗土樣

試驗用土為蘭州地區(qū)非飽和黃土，其基本物理力學(xué)性質(zhì)見表1。

1.2 試驗儀器

試驗儀器為智能人工氣候模擬箱，見圖1。試驗采用控制變量法將濕度、風(fēng)速和光照等級設(shè)為定值（濕度50%，風(fēng)速2 m/s，光照度22 000 lx），通過控制不同的環(huán)境溫度（控溫精度0.1 ℃），研究大氣作用下環(huán)境溫度對土體水分蒸發(fā)的影響。

1.3 試驗步驟

（1）將原狀黃土試樣風(fēng)干粉碎，過2 mm篩孔，按相關(guān)試驗要求將初始含水率配置為15%，并用塑料布覆蓋靜置48 h以上，保證土體中水分擴散均勻。

（2）試驗土體在直徑16 cm、高21 cm的PVC桶中制備，桶內(nèi)每層填土厚度為2 cm，壓實度控制為0.8。相鄰兩層填土的界面用小刀刮毛，使上下層土體充分接觸。

（3）土體裝完后，在其底部和側(cè)壁用保溫材料包裹，頂面敞開，使土體在試驗過程中只有頂面可與外界接觸，其他部位與外界沒有水分、熱量交換，保證桶內(nèi)土體處于一維蒸發(fā)狀態(tài)。

（4）將制備好的試樣桶和電子稱（精度0.01 g）一起放入智能人工氣候模擬箱，研究不同蒸發(fā)時長下土體水分變化規(guī)律。

1.4 試驗方案與過程控制

試驗共配置4組試樣，分別置于25、30、35、40 ℃的氣候模擬箱中進行蒸發(fā)試驗。其他控制條件相同，即濕度為50%，風(fēng)速為2 m/s，光照度為22 000 lx。試驗期間每天9：00打開機器，17：00關(guān)閉，每2 h通過電子稱讀取土體質(zhì)量及水桶的質(zhì)量，研究水分的蒸發(fā)規(guī)律。每組試樣配置5個平行試樣，在試驗進行的第5、8、11、13、15 d的17：00依次取一個土樣，然后立即沿深度方向剖開，采用烘干法測定不同深度處土體的含水率。試驗過程中在氣候模擬箱中放一個與蒸發(fā)試樣等蒸發(fā)界面的水桶，用來測定相同環(huán)境條件下自由水面的潛在蒸發(fā)量。

2 試驗結(jié)果及分析

2.1 土體含水率變化規(guī)律

2.1.1 含水率沿深度分布規(guī)律

圖2為不同環(huán)境溫度下土體蒸發(fā)過程中含水率沿深度方向變化曲線。從圖2可看出，土體表面含水率自試驗開始一直在降低，試驗結(jié)束時土樣表面的含水率最小值在5%以上（40 ℃下15 d土樣），說明非飽和土體蒸發(fā)速率自試驗開始就由土體中水分的多少及其特性決定，不像飽和土體蒸發(fā)開始階段，蒸發(fā)速率是由外界大氣蒸發(fā)強度決定的;當(dāng)蒸發(fā)時長一定時，距離土體表層越近土體含水率變化越大，距離土體表層越遠土體含水率變化越小，原因主要是土體蒸發(fā)從表面開始，下層土體中的水分通過毛細管水勢傳輸?shù)秸舭l(fā)面，對非飽和土來說，蒸發(fā)時長一定，土體表面的蒸發(fā)速率遠大于毛細管勢作用下水分的傳輸速率，導(dǎo)致土體表面含水率減小最顯著，越接近土體底部其含水率減小越少，且含水率沿深度變化量越小;在蒸發(fā)初期（5 d），土體底部的含水率已減小，說明在試驗條件下，蒸發(fā)影響深度已達到土體底部。

2.1.2 含水率隨時間變化規(guī)律

土體含水率隨蒸發(fā)時長的變化規(guī)律曲線見圖3（圖3中Δω/t表示土體不同深度處，某一時段內(nèi)土體含水率的平均日減少率）。從圖3可看出，環(huán)境溫度一定時，土體含水率隨蒸發(fā)時間的日減少率不僅與蒸發(fā)階段有關(guān)，還沿其深度方向不一樣;蒸發(fā)初期土體含水率日減少率最大，隨著蒸發(fā)時間的延長，土體含水率日改變率逐漸減小，直到蒸發(fā)后期蒸發(fā)時長對含水率日減少率的影響不再明顯，如在距離土體表層1 cm處0～5 d、5～8 d、8～11 d、11～13 d、13～15 d時段內(nèi)土體含水率的日減少率依次為1.186%/d、0.974%/d、0.748%/d、0.694%/d、0.670%/d，距離土體表層19 cm處含水率的日減少率依次為0.520%/d、0.456%/d、0.432%/d、0.426%/d、0.420%/d;蒸發(fā)進行到11 d后，土體含水率的日減少率沿深度的分布曲線接近重合，說明在此后一定蒸發(fā)時間內(nèi)，同一深度處土體含水率的日減少率接近，即蒸發(fā)速率接近。

2.2 環(huán)境溫度對土體蒸發(fā)的影響規(guī)律

2.2.1 環(huán)境溫度對土體含水率的影響

環(huán)境溫度是大氣蒸發(fā)強度的主要影響因素之一，對非飽和土體蒸發(fā)來說，環(huán)境溫度對其蒸發(fā)后含水率的分布會產(chǎn)生直接影響。圖4為不同溫度下，距離土體表層1 cm和19 cm處的含水率隨時間變化曲線，實線表示距離土體表層1 cm處，虛線表示距離土體表層19 cm處。從圖4可看出大部分時段環(huán)境溫度對土體表層的含水率影響較土體底部更明顯。從蒸發(fā)開始到第8 d，距離土體表層1 cm處，土體含水率的減小與蒸發(fā)時間幾乎呈線性關(guān)系，溫度為25、30、35、40 ℃時，所對應(yīng)的斜率依次為0.53、0.68、0.77、0.97;蒸發(fā)到第8 d以后，40 ℃下，距離土體表層1 cm處含水率隨時間的變化斜率明顯較之前減小，但含水率仍大于5%，仍在繼續(xù)減小，而在25、30 ℃和35 ℃下，與各自8 d之前相比沒有明顯減小;對比4個溫度條件下8 d后的含水率隨時間的變化斜率還可發(fā)現(xiàn)，此時25 ℃下對應(yīng)的斜率反而較其他3種溫度下更大一些，這是由于25 ℃下土體在0～8 d內(nèi)蒸發(fā)量最少，土體中剩余的水分較多，在8 d后土體表層的水分較多，且土體底部的水由毛細管作用傳送至表層也較多，導(dǎo)致其表層蒸發(fā)更快，含水率隨時間的變化更明顯。在距離土體表層19 cm處，0～5 d內(nèi)溫度越高含水率隨時間的變化斜率越大;5～11 d內(nèi)，4個溫度條件下含水率隨時間的變化斜率基本相等;11 d后，40 ℃下土體底部含水率隨時間的變化斜率明顯減小，且小于其他3個溫度下土體底部含水率隨時間的變化斜率。

當(dāng)環(huán)境溫度改變，試驗結(jié)束后，土體不同深度處溫度變化對土體含水率的影響見表2。分析表2中數(shù)據(jù)發(fā)現(xiàn)，除距離土體表層1 cm外，在土體其他深度處，溫度由25 ℃升至30 ℃時含水率減少量明顯大于相應(yīng)深度處溫度由30 ℃升至35 ℃、35 ℃升至40 ℃時的，說明環(huán)境溫度每升溫1 ℃所導(dǎo)致的含水率減少量不一定完全相等，且與升溫前的參考溫度有關(guān);當(dāng)溫度由25 ℃升至30 ℃時，土體含水率減少最明顯的深度范圍是土樣的中間部位，溫度由30 ℃升至35 ℃、35 ℃升至40 ℃時，土體含水率減少最明顯的部位是土體表層。由于試驗數(shù)據(jù)有限，不能得出更具體的結(jié)論，但可以看出溫度對土體含水率的影響是復(fù)雜的，相同的升溫值導(dǎo)致土體含水率的變化量是不相等的，不同的溫度范圍內(nèi)土體含水率敏感變化深度不同。

2.2.2 溫度對蒸發(fā)速率的影響

不同環(huán)境溫度下非飽和土體每天的蒸發(fā)速率變化見圖5。從圖5可看出，在蒸發(fā)前期，隨環(huán)境溫度升高，蒸發(fā)速率增大，但蒸發(fā)到第5 d之后，隨環(huán)境溫度升高，蒸發(fā)速率反而降低，這是因為環(huán)境溫度越高后期土體中的水分越少所致，如環(huán)境溫度為25、30、35、40 ℃時，蒸發(fā)到第5 d土體表層含水率依次為12.48%、11.07%、10.85%、9.07%;當(dāng)環(huán)境溫度為30、35、40 ℃時，隨蒸發(fā)時間的延長，蒸發(fā)速率呈遞減趨勢，環(huán)境溫度為25 ℃時，蒸發(fā)速率隨時間發(fā)展較平緩，由于夜晚環(huán)境溫度降低，水的飽和蒸汽壓減小，當(dāng)模擬箱內(nèi)空氣中水蒸汽壓接近飽和蒸汽壓時，水蒸汽發(fā)生冷凝變成液態(tài)水，根據(jù)水的飽和蒸汽壓與溫度之間的關(guān)系式[11]，計算得到環(huán)境溫度由30 ℃降低到18.5 ℃，或由40 ℃降低到27.5 ℃時，箱體內(nèi)空氣中的水分會發(fā)生冷凝成液態(tài)水，進而滴落至土樣表面，滲入土體表層，使得夜晚土體含水率增大，當(dāng)環(huán)境溫度為25 ℃時，箱體內(nèi)夜晚溫度需要降低到14 ℃左右時才會發(fā)生冷凝，此時室內(nèi)夜晚溫度（18 ℃左右）高于此溫度，因此箱體內(nèi)空氣中水蒸氣不會發(fā)生冷凝，夜間土體蒸發(fā)速率的減小不明顯，使得土體水分每天的蒸發(fā)速率沒有明顯變化。另外，此環(huán)境溫度下試驗期間土體中水分減少得較少，也是每天的蒸發(fā)速率變化平緩的原因。

2.3 累計蒸發(fā)量的變化規(guī)律

不同環(huán)境溫度下非飽和土體的累計蒸發(fā)量隨時間的變化曲線見圖6。圖7為潛在累計蒸發(fā)量隨時間的變化曲線，作為外界環(huán)境蒸發(fā)強度參考。從圖6、圖7可看出，環(huán)境溫度一定時，非飽和土體累計蒸發(fā)量均隨蒸發(fā)時間的延長呈波動遞增，潛在累計蒸發(fā)量也呈現(xiàn)類似規(guī)律，這主要是因為在試驗條件下，夜晚環(huán)境溫度處于降溫情況，在密閉箱體內(nèi)，空氣中的水分不能散發(fā)到箱體外且氣候模擬箱本身具有冷凝作用，使得箱內(nèi)空氣中的水分冷凝成液態(tài)水，滴落至土體表面，進而滲入土體表層，使土體的蒸發(fā)量在夜晚不增加反而減小;隨著環(huán)境溫度的升高，土體在蒸發(fā)初期蒸發(fā)量顯著增加，當(dāng)蒸發(fā)到某一時刻，40 ℃溫度下土體波動曲線逐漸趨于平緩，35 ℃溫度下土體波動曲線在試驗后期開始轉(zhuǎn)入平緩，由于試驗時間較短，環(huán)境溫度為25 ℃和30 ℃下沒有體現(xiàn)出來，這主要是因為一方面溫度的升高促使水分子的運動加劇，使水分子的運動動能增大，縮短了水分子逃逸到空氣中所需要的時間;另一方面溫度的升高使土體的持水能力減小、降低水的黏滯性及減小顆粒間毛細水彎液面上的表面張力，使得土體的滲透性增強[12]，進而使得累計蒸發(fā)量趨于最后的平衡所需要的時間縮短。

3 結(jié) 論

通過對非飽和黃土進行不同溫度下水分蒸發(fā)試驗，探討了不同環(huán)境溫度對蒸發(fā)過程中土體含水率、蒸發(fā)速率的影響規(guī)律，分析非飽和土體蒸發(fā)過程中含水率分布規(guī)律和累計蒸發(fā)量隨時間的變化規(guī)律，主要得出以下結(jié)論：

（1）環(huán)境溫度對土體含水率的影響規(guī)律與土體深度、蒸發(fā)時段均有關(guān);蒸發(fā)初期，環(huán)境溫度越高表層的含水率減小越明顯，到蒸發(fā)后期，環(huán)境溫度較低土體含水率變化率反而更明顯;在蒸發(fā)試驗的大部分時段土體底部含水率受溫度的影響沒有表層明顯，底部含水率在不同時段受溫度的影響規(guī)律基本與土體表層類似，只是時間轉(zhuǎn)變點相對滯后于表層。

（2）其他條件等同前提下，相同的升溫值導(dǎo)致土體含水率的變化量是不相等的，不同的溫度范圍內(nèi)土體含水率敏感變化深度不同。

（3）環(huán)境溫度與室溫接近時，土體夜晚蒸發(fā)速率減小不多，每天的蒸發(fā)速率隨時間沒有明顯變化，但當(dāng)環(huán)境溫度高于室溫一定值時，受夜晚空氣中水蒸汽冷凝影響，夜晚土體中含水率會增加，夜晚蒸發(fā)速率減小明顯，土體每天蒸發(fā)速率隨時間呈遞減趨勢，在蒸發(fā)初期，環(huán)境溫度越高遞減越明顯。

（4）隨著環(huán)境溫度的升高，土體在蒸發(fā)初期蒸發(fā)量顯著增加，進而使得累計蒸發(fā)量趨于平衡所需要的時間縮短，且由于冷凝作用使得土體累計蒸發(fā)量呈現(xiàn)波動增加趨勢。

參考文獻：

[1] 陳守義.考慮入滲和蒸發(fā)影響的土坡穩(wěn)定性分析方法[J].巖土力學(xué)，1997，18（2）：8-12，22.

[2] 楊洋，姚海林，盧正.蒸發(fā)條件下路基對氣候變化的響應(yīng)模型及影響因素分析[J].巖土力學(xué)，2009，30（5）：1209-1214，1220.

[3] 王鐵行，陸海紅.溫度影響下的非飽和黃土水分遷移問題探討[J].巖土力學(xué)，2004，25（7）：1081-1084.

[4] 王鐵行，陳晶晶，李彥龍.非飽和黃土地表蒸發(fā)的試驗研究[J].干旱區(qū)研究，2014，31（6）：985-990.

[5] 王鐵行.非飽和黃土路基水分場的數(shù)值分析[J].巖土工程學(xué)報，2008，30（1）：41-45.

[6] 劉小軍，閆玉闖，羅揚.非飽和黃土水分遷移規(guī)律及遷移量的試驗研究[J].西安建筑科技大學(xué)學(xué)報（自然科學(xué)版），2015，47（1）：56-61.

[7] 滕繼東，單鋒，張升，等.考慮風(fēng)速影響的土表蒸發(fā)計算方法[J].巖土工程學(xué)報，2018，40（1）：100-107.

[8] 唐朝生，施斌，顧凱.土中水分的蒸發(fā)過程試驗研究[J].工程地質(zhì)學(xué)報，2011，19（6）：875-881.

[9] 吳珺華，袁俊平，何利軍，等.蒸發(fā)條件下土體內(nèi)部水分分布及變化規(guī)律研究[J].防災(zāi)減災(zāi)工程學(xué)報，2013，33（3）：290-294.

[10] 歐陽斌強，唐朝生，王德銀，等.土體水分蒸發(fā)研究進展[J].巖土力學(xué)，2016，27（3）：625-634.

[11] 盧寧，力科思.非飽和土力學(xué)[M].北京：高等教育出版社，2012：29-30.

[12] ROMERO E， GENS A， LLORET A. Temperature Effects on the Hydraulic Behaviour of an Unsaturated Clay[J].Geotechnical and Geological Engineering，2001，19： 311-332.

【責(zé)任編輯 趙宏偉】