冉洪伍，范繼輝，黃 菁,3

（1.中國科學(xué)院水利部成都山地災(zāi)害與環(huán)境研究所 山地表生過程與生態(tài)調(diào)控重點實驗室，四川 成都 610041；2.中國科學(xué)院大學(xué)，北京 100049；3.武漢理工大學(xué)資源與環(huán)境工程學(xué)院，湖北 武漢 430070）

凍融過程是指在多年凍土和季節(jié)凍土區(qū)，隨著周期性的氣溫變化導(dǎo)致地表土層反復(fù)凍結(jié)和融化的過程。凍融過程中土壤水熱發(fā)生著復(fù)雜的變化，對于寒區(qū)的農(nóng)業(yè)生產(chǎn)和建筑設(shè)施都有著很大的影響[1]。早在20世紀(jì)40年代，人們就發(fā)現(xiàn)了在農(nóng)業(yè)生產(chǎn)中土壤凍結(jié)后產(chǎn)生的冰晶阻塞土體孔隙，影響根系生長，土壤凍融造成土壤體積的變化導(dǎo)致房屋受損，一些學(xué)者開始針對凍融過程的水熱問題展開研究，但只停留在對于過程的定性描述上。直至60年代土水勢被定義之后才開始對凍土的水熱問題開始定量研究。Hoekstra[2]與Dirksen和Miller[3]通過室內(nèi)實驗發(fā)現(xiàn)了凍融過程水分運動變化的基本規(guī)律。隨著土壤水勢理論的不斷完善和土壤參數(shù)實時監(jiān)測技術(shù)的實現(xiàn)，凍土水熱問題定量化研究得以進(jìn)一步發(fā)展。20世紀(jì)90年代，徐學(xué)祖和鄧友生[4]對正凍土和已凍土進(jìn)行室內(nèi)土柱模擬實驗和野外開放系統(tǒng)現(xiàn)場實驗，為后來凍土野外實驗積累了寶貴經(jīng)驗。相比于凍土的水熱問題，凍土的凍脹問題早在17世紀(jì)就被人們關(guān)注[5]，但并不能解釋凍脹的產(chǎn)生原因，直到20世紀(jì)對凍土水熱問題的研究發(fā)現(xiàn)，凍脹產(chǎn)生的根本原因是水分遷移[6]，凍土的水熱力研究從單獨研究力學(xué)性質(zhì)、水熱性質(zhì)到水熱力耦合的研究。近年來，隨著人類在寒區(qū)活動的增加以及全球變化問題的嚴(yán)峻，凍土水熱力耦合的研究成為了凍土研究的熱點，眾多的水熱力耦合模型也被提出并廣泛運用于實際問題。

凍土的水熱力耦合是溫度、水分和應(yīng)力相互作用的一個極其復(fù)雜的熱力學(xué)、物理化學(xué)和力學(xué)的綜合問題。凍土中由于冰的存在使得土壤水熱特性和未凍土差異巨大，在土壤孔隙中的冰對土壤水流運動產(chǎn)生阻礙，使得土壤導(dǎo)水率降低，進(jìn)而影響導(dǎo)熱率等其他土壤水熱特性，造成熱量和水量的重新分配。加之凍土中水的相變，使得水熱運移有很大的耦合性。土壤在凍結(jié)過程中，液態(tài)水凍結(jié)成冰，阻礙水分的運移，在水勢梯度的作用下，未凍水向凍結(jié)區(qū)轉(zhuǎn)移，致使凍結(jié)區(qū)的含冰量增加；隨著溫度的上升，土壤表層開始逐漸融化，下部由于溫度較低，依然有冰的存在，由于冰的阻礙作用，上層融化的水無法下滲，使得上層土壤含水率增加；同時在這些過程中，土體由于凍脹作用發(fā)生變形[7]。凍土水熱力耦合問題作為水循環(huán)的重要內(nèi)容，在農(nóng)業(yè)、水資源、環(huán)境、建筑工程等領(lǐng)域受到高度重視。

1 凍土凍融過程研究

1.1 凍融過程土壤水分變化研究

土壤在凍融過程中水分發(fā)生著液態(tài)水和水汽的遷移過程[6]。20世紀(jì)80年代，Jame等[8]針對非飽和土壤中的水分遷移做了初步研究，發(fā)現(xiàn)在液態(tài)水凍結(jié)成冰的過程中，隨著未凍水含量的降低，土壤基質(zhì)勢也較小，未凍水不斷向凍結(jié)鋒面遷移。在融化過程中，隨著溫度的升高，土壤上層的冰開始融化，液態(tài)水在重力作用下向下遷移，使得融化鋒面出現(xiàn)含水率峰值。水汽遷移由土壤中的水汽密度的梯度引起，基質(zhì)勢和土壤的溫度決定了水汽密度，但在研究中這部分水分遷移量較少，通常忽略不計。

由于凍土的特殊性，凍土中水分運動研究受到許多學(xué)者的重視。Singh等[9]通過水熱傳遞實驗對比了凍結(jié)土與未凍結(jié)土的水分運移特性，得到了在不同溫度下水分運移的函數(shù)；Shoop等[10]根據(jù)野外監(jiān)測數(shù)據(jù)建立了非飽和土壤水分遷移的數(shù)學(xué)模型并通過田間實驗的驗證。在國內(nèi)也針對凍融土壤水分運移做了大量研究，周德源[11]對內(nèi)蒙古河套地區(qū)季節(jié)性凍土凍脹規(guī)律展開了研究，揭示了河套地區(qū)季節(jié)凍土的水分運移規(guī)律；李偉強等[12]通過室內(nèi)土柱實驗和大田試驗，得到了硬殼覆蓋條件下季節(jié)性凍土的水分運移規(guī)律。目前，全球氣候變化問題越來越嚴(yán)峻，許多國家對凍土的研究也很重視，美國、加拿大、瑞典等國家多次召開有關(guān)專題會議，推動了凍土水分運動定量研究的發(fā)展[13]。

1.2 凍融過程土壤熱量變化

在凍融過程中，水分的凍結(jié)和融化必將伴隨著熱量的交換。發(fā)生的熱量變化主要由下面3個過程引起[6]：一是溫度梯度引起的熱量由高溫區(qū)向低溫區(qū)的傳導(dǎo)過程；二是液態(tài)水凍結(jié)或冰融化產(chǎn)生的熱量釋放或吸收，這個過程是土壤凍融過程所特有的現(xiàn)象；三是土壤與外界環(huán)境之間的熱量交換。同時，凍融過程也會對土壤的熱傳導(dǎo)率產(chǎn)生影響。Campbell[14]的研究發(fā)現(xiàn)，土壤熱傳導(dǎo)率和液態(tài)含水率呈現(xiàn)顯著正相關(guān)關(guān)系，土壤在凍結(jié)過程中，液態(tài)水相變一定程度上降低了土壤的熱傳導(dǎo)率。同時，凍融時的熱傳導(dǎo)率還受到溫度、土壤質(zhì)地的影響。譚賢君等[15]針對不同含水率、質(zhì)地的巖土導(dǎo)熱系數(shù)在凍融狀態(tài)下的變化進(jìn)行了研究。另外，凍融過程中熱量的改變和水分的變化密不可分，許多學(xué)者將水分遷移和熱量遷移相結(jié)合展開了研究。Jame等[8]基于Harlan模型對非飽和季節(jié)凍土的熱量遷移做了數(shù)值模擬實驗并通過田間試驗的驗證，對凍融過程的熱量遷移進(jìn)行了定量的描述；Wollschlaeger等[16]通過探地雷達(dá) (ground penetrating radar, GPR)監(jiān)測了不同地表狀況下的凍土水熱特征，全面描述了區(qū)域土壤內(nèi)部水熱分布和動態(tài)變化。

1.3 凍融過程土壤力學(xué)性質(zhì)變化

土壤中的孔隙水在溫度降至臨界溫度時凍結(jié)成冰，凍結(jié)過程中水的相變造成體積增大，對周圍顆粒產(chǎn)生擠壓作用，使得土壤大顆粒不斷破碎成小團(tuán)聚體，小團(tuán)聚體向相對大的顆粒聚集，使得土壤顆粒變得更加緊實，土體縫隙增多增大。同時，未凍水不斷向凍結(jié)鋒面遷移，使得土體被進(jìn)一步分割[6]。由此可見，凍融可以改變土壤物理結(jié)構(gòu)和穩(wěn)定性。

凍土力學(xué)性質(zhì)的變化是寒區(qū)工程建設(shè)中的重要研究內(nèi)容之一，受到眾多國內(nèi)外學(xué)者的關(guān)注。Lehrsch等[17-18]的研究顯示，當(dāng)土壤黏粒含量在15%～41%內(nèi)，團(tuán)聚體的穩(wěn)定性隨土壤中黏粒含量增加而增加，當(dāng)土壤黏粒含量大于41%，團(tuán)聚體穩(wěn)定性幾乎不受到含水量和凍融交替次數(shù)的影響，當(dāng)土壤黏粒含量小于15%，團(tuán)聚體穩(wěn)定性隨土壤凍融次數(shù)的增加和含水量的增加而降低。同時，凍融過程中土壤的結(jié)構(gòu)改變，將直接導(dǎo)致土壤特性參數(shù)的改變。鄧西民等[19]對壤質(zhì)粘土進(jìn)行室內(nèi)模擬凍融試驗發(fā)現(xiàn)，凍融處理后土壤容重降低3.5%～9.3%，孔隙度增大6.1%～16.3%，飽和導(dǎo)水率提高1.4～7.7倍。

2 凍融過程水熱力耦合作用機理研究

2.1 水分遷移的驅(qū)動力

凍土中的未凍水在溫度梯度或壓力梯度的驅(qū)動下產(chǎn)生運移，從而誘發(fā)鹽分的擴散機制，同時造成土顆粒的位移[6]，另外，土體中的水分遷移也會改變凍土的冰水比例和結(jié)構(gòu)[20]，因此，確定凍土中水分遷移的驅(qū)動力是解決凍土水熱耦合問題的關(guān)鍵，也是防治土壤鹽漬化、凍脹、下沉等問題的理論基礎(chǔ)。目前凍土中的水分遷移機理主要有兩種理論[21]：第一種是土水勢理論，認(rèn)為水分遷移的主要驅(qū)動力來自土水勢。但試驗中土壤冰水相壓力難以測定，通常采用Clausius-Clapeyron方程[22]來描述某一時刻冰水壓之間的關(guān)系。此類方法只能描述單物質(zhì)平衡狀態(tài)下物理量的變化，對于動態(tài)問題并不是特別適用。第二種為水動力理論，認(rèn)為土壤凍融過程中發(fā)生了水分遷移、熱量傳遞、水分相變3個主要的物理過程，且和水分場、溫度場相互作用、相互影響。

為確定水分遷移驅(qū)動力，學(xué)者曾提出毛細(xì)力、結(jié)晶勢、化學(xué)勢等多種假說，但由于當(dāng)時技術(shù)條件的限制，并沒有對這些假說進(jìn)行驗證[6]。1935年Beskow[23]提出細(xì)顆粒薄膜水遷移理論和吸附-薄膜理論，認(rèn)為未凍水薄膜厚度是溫度的函數(shù)，水從水薄膜較薄、水分子較活躍處遷移。之后便提出了水勢的概念，Harlan[24]認(rèn)為水勢梯度是水分遷移的動力，利用Clausius-Clapeyron方程求得土水勢。之后，Taylor和Luthin[25]在Harlan理論的基礎(chǔ)上進(jìn)行了改進(jìn)，認(rèn)為未凍水含量梯度是水分遷移驅(qū)動力，根據(jù)未凍水含量和土水特征曲線來確定土水勢。Horiguchi[26]認(rèn)為溫度梯度和壓力梯度共同作用于水分遷移。但是，水分遷移是力學(xué)、物理和物理化學(xué)三者作用之和，上述假說也只能說明特定條件下的水分遷移驅(qū)動力[6]。20世紀(jì)60年代初期，國際土壤學(xué)會提出了土壤水勢的定義及劃分，指出土壤中的總水勢是壓力、重力、溫度等分勢的總和，水分遷移是其中各個分勢梯度作用的結(jié)果[4]。經(jīng)過長期試驗研究，證明了影響凍土水分遷移的重要因素有地下水埋深、溫度、土質(zhì)條件等。但是由于凍土結(jié)構(gòu)和水分遷移過程的復(fù)雜性，至今任何一個模型都不能完全模擬現(xiàn)實條件下的凍土水分遷移。因此，還需對這些水分遷移驅(qū)動力模型進(jìn)行完善。

等我從衛(wèi)生間里出來，竟然和那位小說家不期而遇。讓我想不到的是那個寫童話的女孩一臉笑容，緊緊地跟在小說家的身邊。我剛想和他們打招呼，但小說家卻神色尷尬地轉(zhuǎn)過身去……

2.2 土壤凍結(jié)凍脹問題研究

土體凍脹可以分為原位凍脹和分凝凍脹?？紫端l(fā)生凍結(jié)時位置不發(fā)生改變稱作原位凍脹，原位凍脹造成體積增大9%。由外界水分補給并在土中發(fā)生凍結(jié)稱作分凝凍脹，體積增加1.09倍[6]，形成的冰透鏡體成為分凝冰。對于開放系統(tǒng)來說，分凝凍脹是構(gòu)成土體凍脹的主要分量。因此，研究分凝冰的形成機制以及發(fā)育過程是解決土體凍脹問題的關(guān)鍵。

國內(nèi)外的學(xué)者在基于長期的試驗觀察下，對分凝冰的生成和凍脹過程的解釋提出了不同的理論觀點。在20世紀(jì)，人們逐漸認(rèn)識到導(dǎo)致土體凍脹的根本原因是土壤水分遷移。20世紀(jì)60年代，Everett[27]根據(jù)毛細(xì)理論分別對凍脹和凍脹力進(jìn)行定量的估算和解釋，在當(dāng)時這一理論被廣為接受并得以發(fā)展，并被稱作第一凍脹理論。Williams[28]根據(jù)這一理論對凍土的性質(zhì)以及變化特征進(jìn)行了總結(jié)；Jones等[29]根據(jù)這一理論提出了一種測定巖石吸力特性的方法。但是，毛細(xì)理論卻并不能解釋分凝冰的形成原因[30]，Miller等[31]認(rèn)識到這些不足從而提出了第二凍脹理論，也被稱作凍結(jié)緣理論。他認(rèn)為在冰透鏡底面和凍結(jié)鋒面之間存在一個低導(dǎo)濕率、低含水率、無凍脹的凍結(jié)緣帶。第二凍脹理論克服了毛細(xì)理論的不足，得到很多學(xué)者的認(rèn)可，國內(nèi)外在此基礎(chǔ)上對凍結(jié)緣進(jìn)行深入研究并取得了大量的成果[32-34]。目前，隨著國內(nèi)外工程實踐的發(fā)展需求，對凍結(jié)緣的研究已經(jīng)成為揭示凍脹機理、解決融沉災(zāi)害問題的基礎(chǔ)和重點[5]。

2.3 凍融過程水熱力耦合作用

從宏觀上來認(rèn)識溫度梯度下土壤水分遷移分為3個過程[35]：首先，當(dāng)溫度降至凍結(jié)的臨界溫度時，由于土體和水的相互作用導(dǎo)致一部分水處于凍結(jié)狀態(tài)，一部分水處于未凍結(jié)狀態(tài)，此時水分遷移以水頭勢作用下的滲流運動為主，稱作直接遷移；隨著溫度的降低，未凍水含量和水壓力梯度均降低，使得未凍水向低溫處轉(zhuǎn)移，在凍結(jié)緣鋒面聚集；當(dāng)溫度遠(yuǎn)低于凍結(jié)溫度時，水分運移表現(xiàn)為從暖段向冷端的蒸汽運移，此時的水分遷移稱作間接(耦合)遷移。土體凍融過程是溫度場、水分場、應(yīng)力場三者相互作用的過程，當(dāng)凍融過程中未凍水向凍結(jié)鋒面轉(zhuǎn)移的同時，產(chǎn)生相變潛熱，影響土體的溫度場；溫度場的變化導(dǎo)致土壤水分特征參數(shù)改變，造成土體凍脹，進(jìn)而導(dǎo)致土體應(yīng)力場的改變；同時，應(yīng)力場又作用于溫度場和水分場。凍土水熱耦合過程中，這種相互動態(tài)作用一直持續(xù)[6]。

3 凍融過程水熱力耦合模型研究

凍土水熱力耦合模型的發(fā)展始于最開始的經(jīng)驗半經(jīng)驗?zāi)Ｐ?。Takashi等[36]建立了凍脹量的半經(jīng)驗?zāi)Ｐ停@些半經(jīng)驗?zāi)Ｐ投际峭ㄟ^現(xiàn)場或室內(nèi)模擬實驗得出的數(shù)據(jù)基于凍脹的物理本質(zhì)建立的，為后來凍脹量定量化和凍土水熱力耦合模型的發(fā)展提供了理論依據(jù)，但是經(jīng)驗和半經(jīng)驗?zāi)Ｐ蛢H僅適用于某一特定土質(zhì)和特定條件，難以全面推廣。后來，在這些經(jīng)驗半經(jīng)驗?zāi)Ｐ偷幕A(chǔ)之上，結(jié)合物理定律，提出了一系列的凍土水熱力耦合模型并在實際研究中應(yīng)用。Masters等[37]研究了土體形變、孔隙壓力和溫度三者間的耦合關(guān)系，Lai等[38]對凍土滲流場和溫度場的耦合問題做了研究，Michalowski等[39]基于熱力學(xué)原理模擬土體中溫度場和應(yīng)力場的變化。發(fā)展到目前，較成熟的水熱力耦合模型主要有水動力學(xué)模型、剛性冰模型和熱力學(xué)模型3類，以及基于這3類基本模型根據(jù)實際需求演變而成的模型。本文對當(dāng)前主流的模型進(jìn)行對比(表1)，并選取當(dāng)前最具代表性的水動力學(xué)模型、剛性冰模型和熱力學(xué)模型進(jìn)行簡要介紹。

3.1 水動力學(xué)模型

水動力力學(xué)模型是基于非飽和土體水分遷移，將凍脹機理歸結(jié)為未凍水含量隨溫度變化的函數(shù)。其中，Harlan[24]提出的水熱耦合模型應(yīng)用最為廣泛，其方程如下：

表1 水熱力耦合過程的代表性模型Table 1 Typical models of the hydrothermal coupling process

3.2 剛性冰模型

剛性冰模型的理論基礎(chǔ)是第二凍脹理論[6]，即在冰透鏡體與凍結(jié)鋒面之間存在著一個無凍脹、低導(dǎo)水率、低含水率的凍結(jié)緣帶。該模型將活動冰透鏡體和位于其下的孔隙冰看作一個剛性的整體，以相同的速度移動[42]，這也被稱作重結(jié)冰假設(shè)[49]。剛性冰模型對分凝冰的形成過程給出了描述：孔隙壓力大于外荷載造成土體骨架的斷裂，形成分凝冰。剛性冰模型的基本方程如下：

相較其他模型，剛性冰模型考慮了冰晶的移動，對凍結(jié)緣內(nèi)參數(shù)、分凝冰產(chǎn)生、凍脹給出了描述，但該模型參數(shù)復(fù)雜且難以得到，在實際問題的研究中應(yīng)用并不廣泛[5]。為了解決模型參數(shù)復(fù)雜的問題，Black[50]發(fā)展了針對該模型的參數(shù)輸入系統(tǒng)，并根據(jù)試驗進(jìn)行了驗證。也有學(xué)者為解決實際的問題對剛性冰模型進(jìn)行改進(jìn)和簡化，Sheng[51]為了解決野外條件下分層土、非飽和土等的凍脹預(yù)報、冰透鏡位置等問題，在模型中考慮應(yīng)力和應(yīng)變的特征，而對模型中的熱質(zhì)遷移方程不考慮，從而在實際工作中得以更好的使用。

3.3 熱力學(xué)模型

20世紀(jì)90年代，F(xiàn)remond等[52]和Duquennoi等[43]基于能量、動量、質(zhì)量以及熵增平衡定律，提出了熱力學(xué)模型。模型考慮水熱遷移和相態(tài)變化產(chǎn)生的孔隙力，從微觀角度對凍土中的水熱耦合進(jìn)行描述[39]。

基于質(zhì)量守恒定律，得出如下的方程：

熱力學(xué)模型可以很好地描述由于孔隙水凍結(jié)、水分和熱量遷移引起的吸力，模型基于一般的力學(xué)公理，有著很高的準(zhǔn)確性，但是其參數(shù)眾多，且大多數(shù)為微觀物性參數(shù)，難以獲得，使用時需對模型簡化才能得到有實際意義的結(jié)果，因此，熱力學(xué)模型后續(xù)應(yīng)用并不廣泛，其意義也主要體現(xiàn)在理解水熱力耦合過程影響因素這一層面[49]。

4 討論

本研究闡述了凍融過程水分遷移、熱量、力學(xué)性質(zhì)的變化，對水熱力耦合理論及模型的發(fā)展、應(yīng)用做了簡要的介紹。雖然國內(nèi)外學(xué)者對凍土水熱耦合問題進(jìn)行了大量的研究，也取得較大的進(jìn)展，但是目前仍然存在如下重要問題有待解決。

1)凍土中的土水特征的研究工作目前大多只停留在未凍含水量和溫度的經(jīng)驗關(guān)系之上，極少有研究涉及其內(nèi)在機理的探究。以Harlan模型對凍融水熱耦合的模擬為例，目前的研究主要是通過室內(nèi)模擬實驗得到未凍水含量和溫度的耦合關(guān)系，再根據(jù)此關(guān)系對凍融過程的水熱過程進(jìn)行模擬，這種方式存在精度不可控、實驗周期長的缺點，因此，提出在土壤凍結(jié)狀態(tài)下基質(zhì)吸力和未動含水量的通用模型也應(yīng)該成為土水特征的研究重點。

2)大多針對凍融耦合的研究和模型都是假設(shè)土體是飽和的，而對于不飽和的土體研究較少。在大多數(shù)多年凍土區(qū)都是非飽和土體，因此，非飽和土體水熱耦合過程及耦合模型還有待完善。針對這個問題，筆者認(rèn)為今后應(yīng)針對不同區(qū)域、不同土壤類型進(jìn)行長時間尺度的野外觀測實驗，并結(jié)合室內(nèi)模擬實驗，深入研究在不同土壤含水率、不同土壤質(zhì)地下的凍融水熱耦合關(guān)系。另外，研究的焦點不應(yīng)局限在凍融過程本身，應(yīng)建立大氣環(huán)境和土壤凍融之間的相互作用關(guān)系，并擴展至不同環(huán)境氣候條件對凍融過程的影響乃至大時空尺度上凍融耦合研究，例如，不同地表覆蓋層對凍融過程的影響，以及在區(qū)域尺度上氣候變化背景下凍融過程的變化等內(nèi)容。

3)分凝冰的形成機制作為揭示土壤凍脹機理和解決凍土地區(qū)沉降問題的基礎(chǔ)和重點，其形成規(guī)律已成為凍土凍脹研究的前沿問題，但目前實驗設(shè)備的技術(shù)水平有限，無法從微觀上對凍結(jié)過程中分凝冰的形成進(jìn)行觀察，因此，引入可視化微型設(shè)備對凍結(jié)緣內(nèi)部的結(jié)構(gòu)變化和分凝冰的形成是很有必要的。另外，目前針對凍結(jié)的實驗大多是在恒溫的條件下進(jìn)行的，而在現(xiàn)實情況下溫度的變化情況復(fù)雜，在研究中可以將溫度的變化看成線性降溫模式，揭示在線性降溫模式下的土壤凍結(jié)規(guī)律，以更好地了解分凝冰的形成規(guī)律和凍脹并更好的解決寒區(qū)建筑工程中的凍脹問題。