水熱耦合模擬中的凍土滲透系數(shù)模型選擇

王青林， 陳 磊， 明 鋒

（1. 青海省湟源公路工程建設(shè)有限公司， 青海 西寧 810000； 2. 青島農(nóng)業(yè)大學(xué) 建筑工程學(xué)院，山東 青島 266000；3. 中國科學(xué)院 西北生態(tài)環(huán)境資源研究院 凍土工程國家重點(diǎn)實(shí)驗(yàn)室， 甘肅 蘭州 730000）

0 引言

我國是第三凍土大國，凍土區(qū)面積約占全國面積的53.5%［1］。受凍土區(qū)氣溫變化的影響，不僅土中原位水發(fā)生相變產(chǎn)生約9%的體積膨脹，而且未凍結(jié)土層中的孔隙水在溫度梯度作用下將向凍結(jié)鋒面遷移并凍結(jié)成冰。由此將引起地表出現(xiàn)大幅度隆起，局部變形量甚至超過40 cm，這將導(dǎo)致修建在寒區(qū)環(huán)境中的建筑物產(chǎn)生凍害［2-4］。為解決因凍脹引起的凍害問題，研究者提出了諸多水熱耦合模型來評(píng)估凍脹［5-7］。在這些凍脹模型中，凍土滲透系數(shù)是控制水分遷移速率及遷移量的關(guān)鍵參數(shù)。因此，選擇合理的滲透系數(shù)模型就成為準(zhǔn)確預(yù)測(cè)凍脹的前提條件。

獲取凍土滲透系數(shù)的方法有兩種：試驗(yàn)測(cè)試和理論推導(dǎo)。盡管飽和凍土滲透系數(shù)試驗(yàn)存在費(fèi)時(shí)、費(fèi)力、離散性大等缺陷，但諸多研究人員仍得到了凍土滲透系數(shù)隨溫度的變化規(guī)律［8-12］。基于已有試驗(yàn)結(jié)果，國內(nèi)外學(xué)者提出了一系列經(jīng)驗(yàn)?zāi)Ｐ秃屠碚撃Ｐ蛠韽浹a(bǔ)試驗(yàn)方法的不足。經(jīng)驗(yàn)?zāi)Ｐ偷奶攸c(diǎn)是形式簡單，但沒有嚴(yán)格的推導(dǎo)過程，僅僅是根據(jù)試驗(yàn)結(jié)果所給出擬合公式，通常具有一個(gè)或多個(gè)經(jīng)驗(yàn)參數(shù)［13］。而且，經(jīng)驗(yàn)?zāi)Ｐ椭袃鐾翝B透性系數(shù)是未凍水含量或者溫度的函數(shù)。在水熱耦合模型中，大多采用經(jīng)驗(yàn)?zāi)Ｐ蛠砻枋鰞鐾翝B透系數(shù)變化規(guī)律，如應(yīng)用較多的冰阻模型［7］。該模型認(rèn)為飽和凍土的滲透系數(shù)由于受到了冰的阻力作用，致使其降低至相同液態(tài)水含量條件下非飽和未凍土滲透系數(shù)的1/I（I是冰阻系數(shù)，其常為含冰量或溫度的函數(shù)）。經(jīng)驗(yàn)?zāi)Ｐ椭械膮?shù)是通過試驗(yàn)結(jié)果擬合得到，所以這些參數(shù)僅適用于特定土體。當(dāng)將其用于其他土體時(shí)，模型有效性將有所降低［14］。所以，在不同水熱模型中，通常會(huì)選擇不同的經(jīng)驗(yàn)?zāi)Ｐ蛠硖岣哳A(yù)測(cè)精度。相對(duì)而言，理論模型具有嚴(yán)格的推導(dǎo)過程，其模型參數(shù)具有明確的物理意義。由于理論模型減少了對(duì)經(jīng)驗(yàn)系數(shù)的依賴，其普適性更強(qiáng)，更值得在水熱模型中推廣使用。因此，本文對(duì)凍土滲透系數(shù)理論模型在水熱耦合模型中的有效性進(jìn)行討論。

由于水熱耦合模擬結(jié)果對(duì)凍土滲透系數(shù)極為敏感，所以合理選擇凍土滲透模型就變得十分重要。為評(píng)價(jià)理論模型在水熱耦合模型中的適用性，本文以四組飽和凍土的滲透系數(shù)理論模型為例，首先將滲透系數(shù)理論模型與經(jīng)驗(yàn)?zāi)Ｐ皖A(yù)測(cè)結(jié)果進(jìn)行對(duì)比，然后選擇預(yù)測(cè)效果較好的滲透系數(shù)模型帶入水熱耦合模型，以凍脹量變化為判據(jù)評(píng)價(jià)其適用性。

1 凍土滲透系數(shù)理論模型

在凍土滲透實(shí)驗(yàn)的基礎(chǔ)上，人們提出了諸多經(jīng)驗(yàn)?zāi)Ｐ图袄碚撃Ｐ停? 列舉了四組飽和凍土的滲透系數(shù)理論模型。根據(jù)模型特點(diǎn)，將四組模型分別命名為分形模型［13］、離散統(tǒng)計(jì)模型［14］、連續(xù)統(tǒng)計(jì)模型［15］和等價(jià)模型［16］，各模型的數(shù)學(xué)表達(dá)式及主要計(jì)算參數(shù)如表1所示。

表1 凍土滲透系數(shù)模型Table 1 Hydraulic conductivity models for frozen soil

上述四組模型的詳細(xì)推導(dǎo)過程，可參閱文獻(xiàn)［13-16］。為更好體現(xiàn)四組滲透系數(shù)模型的優(yōu)缺點(diǎn)，下面對(duì)四組理論模型進(jìn)行簡單介紹。分形模型是將分形理論應(yīng)用到凍土中，可更好地表征土中孔隙真實(shí)無序的狀態(tài)。該模型結(jié)果主要取決于凍土中的最大孔徑，孔隙率、孔徑分維數(shù)和迂曲度，其中最大孔徑對(duì)凍土的滲透系數(shù)影響最大［13］。離散統(tǒng)計(jì)模型是基于土體凍結(jié)特征曲線與土體孔隙分布函數(shù)的相似性，以及非飽和土的滲透系數(shù)與土水特征曲線之間的關(guān)系而提出的，該模型的計(jì)算結(jié)果主要取決于其中有效孔隙分布［14］。連續(xù)統(tǒng)計(jì)模型則是通過建立凍結(jié)特征曲線和孔徑分布函數(shù)的關(guān)系函數(shù)，提出了滲透系數(shù)模型，該模型只需要對(duì)不同未凍水含量條件下的土體凍結(jié)特征曲線進(jìn)行積分即可得到滲透系數(shù)［15］。等價(jià)模型借鑒非飽和土的滲透系數(shù)模型，認(rèn)為孔隙冰與空氣對(duì)水分遷移的阻礙作用是相當(dāng)?shù)模J(rèn)為飽和凍土的滲透系數(shù)等于相同液態(tài)水含量條件下非飽和未凍土的滲透系數(shù)［16］。由于部分文獻(xiàn)中并沒有給出飽和未凍土的滲透系數(shù)，所以只能計(jì)算其相對(duì)滲透系數(shù)進(jìn)行比較。當(dāng)采用相對(duì)滲透系數(shù)表述預(yù)測(cè)結(jié)果時(shí)，只需要土體凍結(jié)特征曲線和土體顆粒級(jí)配曲線。

為量化滲透系數(shù)預(yù)測(cè)值和實(shí)測(cè)值之間的差異，采用均方根誤差（RMSE）進(jìn)行分析，其計(jì)算公式如下：

式中：N是實(shí)測(cè)數(shù)據(jù)數(shù)量；Krti是第i個(gè)實(shí)測(cè)相對(duì)滲透系數(shù)（由經(jīng)驗(yàn)?zāi)Ｐ痛_定）；Krpi是第i個(gè)預(yù)測(cè)相對(duì)滲透系數(shù)（由理論模型確定）。

2 結(jié)果與分析

為評(píng)價(jià)滲透系數(shù)理論模型的適用性，此處選取文獻(xiàn)［17-19］中的三個(gè)水熱耦合模型為例進(jìn)行分析。在這三組水熱耦合模型中，不僅都采用的是飽和凍土的滲透系數(shù)經(jīng)驗(yàn)?zāi)Ｐ停椅墨I(xiàn)中都給出了詳細(xì)的模型參數(shù)。模型驗(yàn)證從兩個(gè)方面開展：滲透系數(shù)變化和凍脹變化。為方便比較，將案例一、案例二和案例三中的滲透系數(shù)經(jīng)驗(yàn)?zāi)Ｐ头謩e命名為經(jīng)驗(yàn)?zāi)Ｐ鸵唬?jīng)驗(yàn)?zāi)Ｐ投徒?jīng)驗(yàn)?zāi)Ｐ腿?/p>

2.1 滲透系數(shù)變化

2.1.1 案例一

案例一為文獻(xiàn)［17］中的水熱耦合模型，該模型采用經(jīng)驗(yàn)?zāi)Ｐ鸵粊砻枋鰞鐾翝B透系數(shù)變化。經(jīng)驗(yàn)?zāi)Ｐ鸵豢梢员硎緸椋?0］：

式中：k為飽和凍土的滲透系數(shù)；ks為飽和未凍土的滲透系數(shù)；θu為體積未凍水含量；θ為初始飽和體積含水量；b為與土顆粒級(jí)配相關(guān)的擬合參數(shù)，可以表示為：

式中：dg為幾何平均顆粒直徑；σg為幾何標(biāo)準(zhǔn)差；dg與σg可由土顆粒級(jí)配曲線獲得。

基于文獻(xiàn)［17］中給出的土體凍結(jié)特征曲線、土顆粒級(jí)配曲線和飽和滲透系數(shù)，利用表1中的模型，可計(jì)算得到四組模型計(jì)算的相對(duì)滲透系數(shù)（圖1）。從圖1可以看出本文所列舉的四組理論模型可以有效反映滲透系數(shù)隨溫度的變化規(guī)律，且預(yù)測(cè)值與文獻(xiàn)［17］中的滲透系數(shù)值比較接近。應(yīng)該指出的是，分形模型存在有效的臨界溫度，對(duì)應(yīng)凍結(jié)過程中的粉質(zhì)黏土、粉土和粉砂分別是-18.0 ℃、-15.0 ℃和-13.7 ℃。當(dāng)溫度低于該臨界溫度時(shí)，分形模型將不再適用。從圖1 中的計(jì)算結(jié)果來看，該臨界溫度較低，可滿足凍土滲透系數(shù)計(jì)算要求。離散統(tǒng)計(jì)模型高估了粉質(zhì)黏土的滲透系數(shù)，但低估了粉砂土的滲透系數(shù)。整體而言，分形模型、等價(jià)模型和連續(xù)統(tǒng)計(jì)模型的表現(xiàn)較好。

圖1 不同土體的相對(duì)滲透系數(shù)隨溫度變化Fig. 1 Variation of the relative hydraulic conductivity with the temperature

圖2 進(jìn)一步在雙對(duì)數(shù)坐標(biāo)中對(duì)比了滲透系數(shù)模型計(jì)算值和文獻(xiàn)［17］中所用的滲透系數(shù)值。在《土工試驗(yàn)方法標(biāo)準(zhǔn)》（GB/T 50123—2019）中，土體滲透系數(shù)的最大試驗(yàn)允許差值為±2×10-n（n為自然數(shù)）。參考此誤差標(biāo)準(zhǔn)，本文確定相對(duì)滲透系數(shù)預(yù)測(cè)值的可接受試驗(yàn)誤差值為±2×10-n。如圖2 所示，大多數(shù)的數(shù)據(jù)點(diǎn)都落在可接受試驗(yàn)誤差區(qū)域內(nèi)，僅少部分?jǐn)?shù)據(jù)點(diǎn)落在可接受實(shí)驗(yàn)誤差區(qū)域外。需要注意的是，有三組滲透系數(shù)數(shù)據(jù)點(diǎn)甚至落在一個(gè)數(shù)量級(jí)誤差區(qū)域外，如粉砂的分形模型，粉土和粉砂的離散統(tǒng)計(jì)模型。總體而言，本文所列四組滲透系數(shù)模型預(yù)測(cè)值和文獻(xiàn)［17］水熱耦合模型中所用的滲透系數(shù)值較接近。這表明本文所列四組滲透系數(shù)理論模型可代替文獻(xiàn)［17］中的經(jīng)驗(yàn)?zāi)Ｐ汀?/p>

圖2 文獻(xiàn)［17］所用相對(duì)滲透系數(shù)與四組理論模型計(jì)算值對(duì)比Fig. 2 Comparison of the relative hydraulic conductivity used in Reference ［17］ and the values calculated by the four hydraulic conductivity models

表2匯總了這四組模型的均方根誤差值。可以看出，四組模型的均方根誤差值均較小，說明了四組模型具有良好的穩(wěn)定性。四種模型均對(duì)粉砂的預(yù)測(cè)效果較差，對(duì)粉土和粉質(zhì)黏土的效果較為準(zhǔn)確。等價(jià)模型的均方根誤差最小，其次是連續(xù)統(tǒng)計(jì)模型、分形模型和離散統(tǒng)計(jì)模型。對(duì)于文獻(xiàn)［17］中的土樣，等價(jià)模型的預(yù)測(cè)能力最好。

表2 模型的均方根Table 2 RMSE of the model

2.1.2 案例二

案例二以文獻(xiàn)［18］中的水熱耦合模型為例進(jìn)行分析，并以粉質(zhì)黏土為研究對(duì)象。文獻(xiàn)［18］中使用的滲透系數(shù)為經(jīng)驗(yàn)?zāi)Ｐ投浔磉_(dá)式為：

式中：模型所需參數(shù)可查詢參考文獻(xiàn)［18］。

圖3對(duì)比了四組理論模型預(yù)測(cè)值與文獻(xiàn)［18］所用的相對(duì)滲透系數(shù)。如圖所示，相對(duì)滲透系數(shù)經(jīng)驗(yàn)?zāi)Ｐ椭蹬c分形模型和離散統(tǒng)計(jì)模型的預(yù)測(cè)值較為接近，但是與連續(xù)統(tǒng)計(jì)模型和等價(jià)模型的預(yù)測(cè)值稍有偏差。鑒于該經(jīng)驗(yàn)?zāi)Ｐ皖A(yù)測(cè)結(jié)果位于本文四組模型預(yù)測(cè)范圍內(nèi)，且被O’Neill 等［21］、周國慶等［22］、周揚(yáng)等［23］多次采用，表明該模型具有一定的合理性。需要指出的是，該模型的指數(shù)對(duì)不同土并不是一個(gè)常數(shù)。對(duì)于等價(jià)模型，根據(jù)顆粒級(jí)配曲線計(jì)算得到其指數(shù)n=12.08，與文獻(xiàn)［18］所用的n=9 有所差別，所以導(dǎo)致相對(duì)滲透系數(shù)預(yù)測(cè)值偏低。對(duì)分形模型，該粉質(zhì)黏土分形模型有效的臨界溫度是-22 ℃，低于該溫度則分形模型不再適用。可以看出，土體顆粒級(jí)配曲線不僅嚴(yán)重影響等價(jià)模型中的指數(shù)n的取值，而且對(duì)土體凍結(jié)特征曲線的數(shù)值和形態(tài)也有較大影響，進(jìn)而影響分形模型、離散統(tǒng)計(jì)模型和連續(xù)統(tǒng)計(jì)模型的預(yù)測(cè)精度。

圖4 在雙對(duì)數(shù)坐標(biāo)中對(duì)比了四組理論模型與經(jīng)驗(yàn)?zāi)Ｐ投臐B透系數(shù)預(yù)測(cè)結(jié)果。從圖4 中可以看出，只有分形模型和離散統(tǒng)計(jì)模型的數(shù)據(jù)點(diǎn)落在了可接受試驗(yàn)誤差區(qū)域內(nèi)，等價(jià)模型的數(shù)據(jù)點(diǎn)則大多落在了可接受實(shí)驗(yàn)誤差區(qū)域的邊緣。連續(xù)統(tǒng)計(jì)模型的數(shù)據(jù)點(diǎn)雖然大多落在了可接受實(shí)驗(yàn)誤差區(qū)域外，但也基本分布在一個(gè)量級(jí)的誤差區(qū)域內(nèi)。因此，可以采用分形模型和離散統(tǒng)計(jì)模型來代替經(jīng)驗(yàn)?zāi)Ｐ投?/p>

圖4 文獻(xiàn)［18］中所用的相對(duì)滲透系數(shù)與四組理論模型計(jì)算值對(duì)比Fig. 4 Comparison of the relative hydraulic conductivity used in Reference ［18］ and the values calculated by the four hydraulic conductivity models

表3 匯總了這四組模型的均方根誤差值，可以看到四組模型的均方根誤差值均較小，說明四組模型的穩(wěn)定性較好。相對(duì)而言，分形模型和離散統(tǒng)計(jì)模型的均方根誤差較小，表明這兩組模型可以代替經(jīng)驗(yàn)?zāi)Ｐ投?/p>

表3 模型的均方根Table 3 RMSE of the model

2.1.3 案例三

案例三以文獻(xiàn)［19］的水熱耦合模型為例進(jìn)行分析。該耦合模型采用的是經(jīng)驗(yàn)?zāi)Ｐ腿浔磉_(dá)式為：

從公式（5）可以看出，經(jīng)驗(yàn)?zāi)Ｐ腿谋磉_(dá)式簡單，只需要溫度和飽和滲透系數(shù)就可以預(yù)測(cè)凍土滲透系數(shù)。圖5給出了本文所列舉的滲透系數(shù)模型計(jì)算的相對(duì)滲透系數(shù)。如圖所示，除等價(jià)模型外，其余三組理論模型與經(jīng)驗(yàn)?zāi)Ｐ腿臐B透系數(shù)預(yù)測(cè)值較為接近。需要指出的是，該粉質(zhì)黏土對(duì)應(yīng)的分形模型有效臨界溫度為-1.40 ℃。對(duì)比前兩個(gè)案例，此粉質(zhì)黏土的分形有效臨界溫度偏高，這是因?yàn)樵撏翗拥膬鼋Y(jié)特征曲線斜率過大，使得大部分水分被迅速凍結(jié)。等價(jià)模型與經(jīng)驗(yàn)?zāi)Ｐ腿念A(yù)測(cè)結(jié)果相差最大，可能是由于土體顆粒級(jí)配曲線不合理，導(dǎo)致計(jì)算得到的指數(shù)n偏離極大所致。盡管經(jīng)驗(yàn)?zāi)Ｐ腿嬖诮?jīng)驗(yàn)參數(shù)，但仍適用于描述此粉質(zhì)黏土樣品滲透系數(shù)變化。

圖5 相對(duì)滲透系數(shù)隨溫度的變化曲線Fig. 5 Variation of the relative hydraulic conductivity with the temperature

圖6 進(jìn)一步在雙對(duì)數(shù)坐標(biāo)中對(duì)比了本文所列舉的四組滲透系數(shù)模型和經(jīng)驗(yàn)?zāi)Ｐ腿?jì)算的相對(duì)滲透系數(shù)數(shù)值。如圖6 所示，連續(xù)統(tǒng)計(jì)模型的數(shù)據(jù)點(diǎn)基本都落在可接受試驗(yàn)誤差范圍內(nèi)；分形模型的數(shù)據(jù)點(diǎn)也大多落在可接受試驗(yàn)誤差區(qū)域內(nèi)；離散統(tǒng)計(jì)模型只有少部分?jǐn)?shù)據(jù)點(diǎn)落在可接受試驗(yàn)誤差內(nèi)，但是大部分都落在了一個(gè)量級(jí)的誤差區(qū)域內(nèi)；等價(jià)模型預(yù)測(cè)結(jié)果則幾乎全部不在可接受試驗(yàn)誤差區(qū)域內(nèi)。這說明等價(jià)模型并不適用于此土樣，而其余三組理論模型可以代替經(jīng)驗(yàn)?zāi)Ｐ腿鐾翝B透系數(shù)變化趨勢(shì)。

圖6 文獻(xiàn)［19］中所用的相對(duì)滲透系數(shù)與四組模型計(jì)算值對(duì)比Fig. 6 Comparison of the relative hydraulic conductivity used in Reference ［19］ and the values calculated by the four hydraulic conductivity models

表4匯總了這四組模型的均方根誤差值。由表4 可見，除了等價(jià)模型，其余三組模型的均方根誤差值均較小，說明了這三組模型的可靠性。其中，連續(xù)統(tǒng)計(jì)模型的均方根誤差最小，其次為分形模型和離散統(tǒng)計(jì)模型，而等價(jià)模型的均方根誤差最大。因此，可以采用分形模型和連續(xù)統(tǒng)計(jì)模型預(yù)測(cè)該凍結(jié)粉質(zhì)黏土樣品的滲透系數(shù)。

表4 模型的均方根Table 4 RMSE of the model

2.2 凍脹變化規(guī)律

凍土滲透系數(shù)模型大多應(yīng)用于預(yù)測(cè)凍土中的水熱分布及凍脹、融沉變形。從前一節(jié)的分析結(jié)果來看，分形模型在三個(gè)案例中的適用性更強(qiáng)。因此，本節(jié)采用分形模型代替文獻(xiàn)［18］中的經(jīng)驗(yàn)?zāi)Ｐ投⒗迷闹械南嚓P(guān)參數(shù)對(duì)凍脹試驗(yàn)結(jié)果進(jìn)行預(yù)測(cè)（圖7）。如圖7 所示，在凍結(jié)前期（0～50 h），無論是經(jīng)驗(yàn)?zāi)Ｐ投€是分形模型，其預(yù)測(cè)的凍脹量與實(shí)測(cè)凍脹量相差很小。在凍結(jié)后期（50～90 h），經(jīng)驗(yàn)?zāi)Ｐ投念A(yù)測(cè)結(jié)果偏小，而采用分形模型的預(yù)測(cè)結(jié)果偏大。由于經(jīng)驗(yàn)?zāi)Ｐ投臐B透系數(shù)預(yù)測(cè)值略小于分形模型的預(yù)測(cè)值，所以采用分形模型的凍脹預(yù)測(cè)結(jié)果偏大。但從凍脹預(yù)測(cè)誤差來看，分形模型的預(yù)測(cè)結(jié)果接近實(shí)測(cè)值（最大相對(duì)誤差小于10%），可證明分形模型的有效性。

圖7 不同滲透系數(shù)模型下的凍脹量預(yù)測(cè)結(jié)果Fig. 7 Prediction of the frost heave with different hydraulic conductivity models

3 討論

3.1 理論模型的優(yōu)缺點(diǎn)

從三個(gè)案例來看，本文所列四組理論模型可以有效的反映文獻(xiàn)中滲透系數(shù)的變化規(guī)律，而且凍脹算例也證實(shí)了理論模型的可靠性。由于選擇不合理的凍土滲透系數(shù)模型，將導(dǎo)致水熱分布預(yù)測(cè)結(jié)果不準(zhǔn)確，進(jìn)而影響凍脹預(yù)測(cè)精度。因此，若在水熱耦合模型中，采用滲透系數(shù)理論模型，不僅可有效避免經(jīng)驗(yàn)公式參數(shù)的選取，而且還能保證預(yù)測(cè)精度。但需要注意的是，滲透系數(shù)理論模型也存在一定的局限性：（1）最大孔徑對(duì)分形模型的預(yù)測(cè)結(jié)果影響較大，且該模型存在溫度下限，但該下限可滿足凍脹模型的要求；（2）離散統(tǒng)計(jì)模型雖然解決了分形模型的缺點(diǎn)，但其計(jì)算結(jié)果并不是一條平滑曲線，難以通過一個(gè)公式擬合，因此很難應(yīng)用到數(shù)值計(jì)算中。（3）連續(xù)統(tǒng)計(jì)模型所需參數(shù)較少，但對(duì)部分土體的適用性較差。（4）等價(jià)模型不僅沒有嚴(yán)格的理論推導(dǎo)，而且對(duì)部分土體的預(yù)測(cè)效果較差。從四組模型對(duì)三個(gè)案例的預(yù)測(cè)結(jié)果來看，不同理論模型表現(xiàn)出不同的預(yù)測(cè)效果，而分形模型在三個(gè)案例中表現(xiàn)較為均衡，其適用性最強(qiáng)。因此，建議在水熱耦合模型中采用分形模型來描述凍土滲透系數(shù)變化。

3.2 凍融過程中的滲透系數(shù)滯后性

從表1 中可以看出，凍土滲透系數(shù)不僅可以采用溫度為變量，而且也可以采用未凍水含量為變量。圖8 以經(jīng)驗(yàn)?zāi)Ｐ鸵唬═arnawski 和Wagner 模型）為基礎(chǔ)，計(jì)算得到凍結(jié)與融化過程中的土體相對(duì)滲透系數(shù)變化規(guī)律［20］。如圖8 所示，在相同溫度條件下，計(jì)算得到的粉質(zhì)黏土的相對(duì)滲透系數(shù)最大，其次是粉土，粉砂的相對(duì)滲透系數(shù)最小；而在相同體積未凍水含量下，正好相反，粉砂的相對(duì)滲透系數(shù)最大，其次是粉土，而粉質(zhì)黏土的相對(duì)滲透系數(shù)最小。這是因?yàn)門arnawski 和Wagner 模型中，土體凍結(jié)特征曲線和土體顆粒級(jí)配曲線共同決定了滲透系數(shù)的大小。土體顆粒級(jí)配曲線不僅決定了模型參數(shù)2b+ 3 的大小，而且也影響土體中的未凍水含量θu。依據(jù)土體顆粒級(jí)配曲線，計(jì)算得到粉質(zhì)黏土，粉土和粉砂的模型參數(shù)2b+ 3 分別是12.22、10.83和8.52。如圖8（a）所示，計(jì)算出來的融化過程的滲透系數(shù)要小于凍結(jié)過程的滲透系數(shù)，這是因?yàn)橥馏w在融化過程中的未凍水含量要小于其在凍結(jié)過程中的未凍水含量。但是隨著溫度的降低，滯后效應(yīng)對(duì)滲透系數(shù)的影響也在逐漸減小。但是采用未凍水含量為變量時(shí)，滲透系數(shù)在凍融過程并沒有體現(xiàn)出滯后效應(yīng)［圖8（b）］，這和以往的研究結(jié)果一致［24］。因此，建議以未凍水含量為變量表示凍土滲透系數(shù)模型。

圖8 凍結(jié)和融化過程中相對(duì)滲透系數(shù)隨溫度（a）和體積未凍水含量（b）的變化曲線Fig. 8 Variation of the relative hydraulic conductivity with the temperature （a） and unfrozen water content （b） during the freezing and thawing process

4 結(jié)論

本文通過對(duì)四組飽和凍土的滲透系數(shù)模型和三個(gè)水熱耦合模型中選取的經(jīng)驗(yàn)?zāi)Ｐ偷挠?jì)算結(jié)果進(jìn)行對(duì)比，討論了四組飽和凍土的滲透系數(shù)模型的預(yù)測(cè)能力，得出以下主要結(jié)論：

（1）水熱耦合模型中所用的滲透系數(shù)經(jīng)驗(yàn)?zāi)Ｐ蛿?shù)值和本文所列舉的四組理論模型數(shù)值大多都非常接近，建議使用滲透系數(shù)理論模型來描述凍土滲透系數(shù)變化規(guī)律。

（2）四組理論模型均能有效描述凍土滲透系數(shù)變化，但理論模型預(yù)測(cè)結(jié)果受土體顆粒級(jí)配曲線影響明顯，推薦使用分形模型計(jì)算凍土滲透系數(shù)。

（3）用溫度表示的滲透系數(shù)表現(xiàn)出滯后效應(yīng)，而用未凍水含量表示的滲透系數(shù)并沒有表現(xiàn)出滯后效應(yīng)，建議在凍土滲透系數(shù)模型中以未凍水含量為變量。