非飽和黃土水分蒸發(fā)溫度效應(yīng)試驗(yàn)研究

張映梅 夏瓊 王旭

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

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

中圖分類號(hào)：TU444 文獻(xiàn)標(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

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

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

1 試驗(yàn)材料與方法

1.1 試驗(yàn)土樣

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

1.2 試驗(yàn)儀器

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

1.3 試驗(yàn)步驟

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

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

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

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

1.4 試驗(yàn)方案與過(guò)程控制

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

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

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

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

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

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

土體含水率隨蒸發(fā)時(shí)長(zhǎng)的變化規(guī)律曲線見(jiàn)圖3（圖3中Δω/t表示土體不同深度處，某一時(shí)段內(nèi)土體含水率的平均日減少率）。從圖3可看出，環(huán)境溫度一定時(shí)，土體含水率隨蒸發(fā)時(shí)間的日減少率不僅與蒸發(fā)階段有關(guān)，還沿其深度方向不一樣;蒸發(fā)初期土體含水率日減少率最大，隨著蒸發(fā)時(shí)間的延長(zhǎng)，土體含水率日改變率逐漸減小，直到蒸發(fā)后期蒸發(fā)時(shí)長(zhǎng)對(duì)含水率日減少率的影響不再明顯，如在距離土體表層1 cm處0～5 d、5～8 d、8～11 d、11～13 d、13～15 d時(shí)段內(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ā)進(jìn)行到11 d后，土體含水率的日減少率沿深度的分布曲線接近重合，說(shuō)明在此后一定蒸發(fā)時(shí)間內(nèi)，同一深度處土體含水率的日減少率接近，即蒸發(fā)速率接近。

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

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

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

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

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

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

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

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

3 結(jié) 論

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

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

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

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

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

參考文獻(xiàn)：

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

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

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

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

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

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

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

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

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

[10] 歐陽(yáng)斌強(qiáng)，唐朝生，王德銀，等.土體水分蒸發(fā)研究進(jìn)展[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é)任編輯 趙宏偉】