水-熱耦合作用下嚴(yán)寒地區(qū)高速鐵路隧道溫度場分布規(guī)律

趙希望, 張 恒

(1.廣東省隧道結(jié)構(gòu)智能監(jiān)控與維護(hù)企業(yè)重點(diǎn)實(shí)驗(yàn)室， 廣州 511440； 2.中鐵隧道局集團(tuán)有限公司勘察設(shè)計(jì)研究院， 廣州 511440； 3.華南理工大學(xué)土木與交通學(xué)院， 廣州 510641)

寒冷地區(qū)隧道的修筑打破了圍巖的原始熱平衡，圍巖中存在的水分會(huì)在寒季凍結(jié)成冰[1]。由于液相轉(zhuǎn)為固相體積增大，進(jìn)而使隧道襯砌背后圍巖產(chǎn)生凍脹力；在暖季，凍結(jié)的冰融化成水，在凍結(jié)時(shí)產(chǎn)生的位移不能完全恢復(fù)，造成的空隙又會(huì)增加滲漏，嚴(yán)重時(shí)可能發(fā)生融陷，這種襯砌背后圍巖的反復(fù)凍融作用使襯砌防水層發(fā)生破損，進(jìn)而產(chǎn)生更大的凍害，造成襯砌結(jié)構(gòu)變形，使隧道運(yùn)營安全受到威脅[2-4]。

眾多專家和學(xué)者對新疆、青海等嚴(yán)寒地區(qū)的隧道進(jìn)行了研究。賴遠(yuǎn)明等[5]運(yùn)用伽遼金法對青海大坂山隧道溫度場和滲流場耦合非線性問題進(jìn)行了分析和研究，通過無量綱分析法和攝動(dòng)技術(shù)求出寒區(qū)圓形隧道溫度場的解析解。高焱等[6]對寒區(qū)隧道洞內(nèi)空氣溫度及圍巖溫度場進(jìn)行了分析，給出了保溫層的鋪設(shè)范圍。陳建勛等[7]首先結(jié)合寒區(qū)公路隧道防凍保溫層的設(shè)計(jì)，對保溫材料現(xiàn)場的保溫性能進(jìn)行了測試，確定了最佳的防凍保溫材料組合和施工，然后運(yùn)用正弦函數(shù)回歸法對隧道拱頂、拱腰、邊墻和路面四個(gè)部位實(shí)測的溫度數(shù)據(jù)進(jìn)行了分析。張學(xué)富等[8]根據(jù)青藏鐵路昆侖山隧道內(nèi)實(shí)測大氣溫度和圍巖溫度，考慮冰-水相變影響，應(yīng)用伽遼金法推導(dǎo)了有限元公式，對昆侖山隧道施工期間融化的圍巖進(jìn)行了回凍預(yù)測研究，分析了隧道鋪設(shè)保溫材料和無保溫材料的溫度場，得出保溫材料對多年凍土區(qū)隧道融化回凍過程所起的作用。郭春香等[9]以昆侖山隧道為例，利用有限元法分別計(jì)算了考慮混凝土水化熱和不考慮混凝土水化熱的隧道圍巖融化回凍過程，分析了混凝土水化熱對圍巖融化、回凍過程的影響。楊旭等[10]基于圍巖溫度場現(xiàn)場測試的結(jié)果，同時(shí)考慮水文地質(zhì)條件、混凝土水化熱、大氣溫度和地溫等因素，比較了季節(jié)凍土區(qū)隧道施加保溫層和未施加保溫層的凍融循環(huán)圈的差異。周元輔等[11-12]推導(dǎo)了多年凍土隧道隔熱層的數(shù)學(xué)優(yōu)化模型，分析了隔熱層用于多年凍土隧道的合理性和經(jīng)濟(jì)性。從以上對寒區(qū)隧道的研究中可以看出，防凍保溫層是寒區(qū)隧道防治凍害的有效措施，確定保溫層鋪設(shè)的合理性和經(jīng)濟(jì)性對寒區(qū)隧道工程有著重要的意義。但是，以上分析中或是以多年凍土區(qū)隧道施工期間造成的圍巖融化為對象，研究襯砌背后圍巖的回凍過程，并不適用于初始融土地區(qū)隧道修筑后，由于圍巖熱平衡破壞造成的襯砌背后圍巖凍結(jié)的問題，如文獻(xiàn)[8-9]；或是在隧道溫度場的分析中只考慮了襯砌、圍巖熱傳導(dǎo)，沒有考慮水分遷移或冰水相變，如文獻(xiàn)[10-12]。

針對融土地區(qū)隧道修筑后，襯砌背后圍巖凍結(jié)的問題，以吉圖琿客運(yùn)專線榆樹川隧道為依托，根據(jù)瞬態(tài)溫度場問題的能量平衡方程和水分遷移方程，考慮水分遷移、冰-水相變等因素，建立了寒區(qū)隧道圍巖水-熱耦合模型，并結(jié)合現(xiàn)場實(shí)測大氣溫度，利用有限元法對嚴(yán)寒地區(qū)隧道溫度場進(jìn)行計(jì)算，比較有、無保溫層對隧道溫度的影響，確定隧道凍害發(fā)生的最不利位置和時(shí)間，分析保溫材料對隧道凍融圈的影響，給出保溫層厚度和導(dǎo)熱系數(shù)的合理值。

1 水-熱耦合作用下傳熱數(shù)學(xué)模型

1.1 溫度場方程

土體在凍結(jié)時(shí)水蒸氣蒸發(fā)耗熱很少，考慮水分遷移和冰水相變，土體內(nèi)發(fā)生的熱傳導(dǎo)[13]可以表示如下。

融化區(qū):

(1)

凍結(jié)區(qū):

(2)

式中：f、u為圍巖凍結(jié)和融化狀態(tài)；λ為熱傳導(dǎo)系數(shù)，W/(m·℃)；T為溫度，℃；t為時(shí)間；c為圍巖的比熱容，J/(kg·℃)；ρ和ρi為分別為圍巖的密度和冰的密度，kg/m3；L為單位質(zhì)量冰水相變潛熱，kJ/kg；θi為體積含冰量。

在凍結(jié)鋒面位置處s(t)，應(yīng)滿足連續(xù)性條件

Tf[s(t),t]=Tu[s(t),t]=Tm

(3)

式(3)中：Tm為土體的凍結(jié)溫度，℃。

1.2 水分場方程

根據(jù)文獻(xiàn)[14-15]，土體內(nèi)的空氣和水蒸氣在孔隙內(nèi)的流動(dòng)帶動(dòng)的水分遷移量相當(dāng)微弱，因此，忽略空氣和水蒸氣遷移對水分的遷移的影響，凍結(jié)和融化過程中圍巖內(nèi)水分遷移滿足Darcy定律，可表達(dá)為

(4)

式(4)中：K為導(dǎo)水系數(shù)，m/s；φ為基質(zhì)勢，Pa；θu為未凍水體積含量；ρw為水的密度，kg/m3。

根據(jù)徐斅祖等[16]的實(shí)驗(yàn)結(jié)果，未凍水含量θu與負(fù)溫始終保持動(dòng)態(tài)平衡的關(guān)系，即

θu=a1|T|-b1

(5)

式(5)中：a1和b1為與土質(zhì)因素有關(guān)的經(jīng)驗(yàn)常數(shù)。

1.3 溫度場方程和水分場方程的耦合

土體中的微分水容量C(m-1)是等溫條件下基質(zhì)勢與體積含水量θu之間的定量關(guān)系[16]，表達(dá)式為

(6)

土壤的水分?jǐn)U散系數(shù)D可以表示為土體微分水容量C與導(dǎo)水系數(shù)K的關(guān)系[16]，即

(7)

Taylor等[17]認(rèn)為未凍結(jié)土的土-水特征曲線也適用于凍結(jié)土的情況。因此，根據(jù)凍土中未凍水含量與溫度的關(guān)系以及土-水特征曲線可確定不同溫度的未凍水基質(zhì)勢，表達(dá)式為

(8)

綜合式(2)、式(4)、式(6)～式(8)，土體凍結(jié)區(qū)的熱傳導(dǎo)方程可以表示為

(9)

用等效熱傳導(dǎo)系數(shù)λe和等效比熱容Ce代替圍巖融化和凍結(jié)狀態(tài)的熱傳導(dǎo)系數(shù)λu、λf和比熱容cu、cf，則考慮冰-水相變和水分遷移的圍巖熱傳導(dǎo)方程可以表示為

(10)

式(1)～式(10)組成了隧道圍巖凍結(jié)過程中的水-熱耦合模型，采用COMSOL Multiphysics 軟件的PDE (partial differential equations)模塊進(jìn)行求解。以T和θ為因變量，得到兩個(gè)廣義型PDE標(biāo)準(zhǔn)形式為

(11)

式(11)中：da為質(zhì)量系數(shù)；U為微分方程因變量，U={T,θ}；Γ為通量向量；f為源項(xiàng)。

將節(jié)點(diǎn)溫度T和含水量θ作為基本未知量，運(yùn)用隱式歐拉向后差分(BDF)法對控制方程中的時(shí)間項(xiàng)進(jìn)行離散，采用非線性迭代修正阻尼牛頓法求解方程。利用有限差分時(shí)間步長的方法，對該微分方程進(jìn)行時(shí)間差分，并通過連續(xù)線性循環(huán)解法求解。

1.4 模型驗(yàn)證

為驗(yàn)證提出的水-熱耦合模型的正確性，采用胡和平等[18]對于張掖壤土的凍結(jié)試驗(yàn)和結(jié)果進(jìn)行數(shù)值模擬。胡和平等采用張掖壤土，在試驗(yàn)室內(nèi)進(jìn)行了一組單向凍結(jié)試驗(yàn)。試樣的高度為13.68 cm，直徑為11.36 cm。試樣的初始溫度如表1所示，試樣的初始體積含水量為22.08%。凍結(jié)實(shí)驗(yàn)時(shí)，試樣的側(cè)面為熱絕緣條件，頂部與底部溫度按照如圖1所示變化，試驗(yàn)在凍結(jié)2 830 min后結(jié)束，并將試樣取出，然后將試樣切片，在烤箱中烘干得到總含水量。

圖1 試樣頂部和底部溫度變化Fig.1 Temperature variations at the top and bottom of the sample

圖2 計(jì)算模型Fig.2 Calculation model

表1 試樣初始溫度分布

計(jì)算時(shí)，取軸對稱計(jì)算模型(圖2)，邊界條件和初始條件如下：

溫度邊界條件：AC為對稱邊界；BD為絕緣邊界；初始溫度場沿深度方向按照表1分布。AB和CD邊界溫度隨時(shí)間變化情況如圖1所示。

滲流邊界條件：AC為對稱邊界；AB、BD和CD都是不透水邊界，即試樣處于封閉環(huán)境中。

試驗(yàn)所用土為張掖壤土，熱物理參數(shù)如表2所示。

表2 張掖壤土熱物理參數(shù)

在土體凍結(jié)過程中，土體孔隙中生成的冰晶體會(huì)阻礙孔隙內(nèi)水分的遷移孔道，從而導(dǎo)致水分?jǐn)U散系數(shù)D和滲透系數(shù)K減小。受土體孔隙內(nèi)冰晶體的影響，土體在凍結(jié)過程中的分水?dāng)U散系數(shù)D和滲透系數(shù)K[19]可以表示為

(12)

(13)

式中：I為阻抗系數(shù)；a2、b2、a3和b3分別為物理試驗(yàn)常數(shù)，其值如表3所示。

表3 張掖壤土物理參數(shù)

圖3所示為試樣在不同凍結(jié)時(shí)刻(34、528、2 830 min)，不同深度處溫度實(shí)驗(yàn)值與計(jì)算值的比較。可以看出，各凍結(jié)時(shí)刻試樣內(nèi)溫度的計(jì)算值和實(shí)驗(yàn)值能夠很好地相吻合，這說明該水-熱模型能夠很好地模擬土體凍結(jié)過程中的溫度分布規(guī)律。同時(shí)，可以看出：凍結(jié)初始階段，試樣處于正溫狀態(tài)，但隨著凍結(jié)的進(jìn)行，試樣發(fā)生自上而下的凍結(jié)，凍結(jié)鋒面不斷下移。凍結(jié)初期溫度傳遞速率較快，隨凍結(jié)時(shí)間的增加溫度傳遞速率逐漸減小，最終溫度趨于穩(wěn)定。

圖3 不同時(shí)刻土樣溫度計(jì)算值與試驗(yàn)值的比較Fig.3 Comparison between the experimental and simulated temperatures of soil sample at different time

圖4為實(shí)驗(yàn)結(jié)束時(shí)試樣內(nèi)含水量的計(jì)算值和實(shí)驗(yàn)值的比較，含水量的實(shí)驗(yàn)值和計(jì)算值在形狀和分布特征上都是一致的。試驗(yàn)開始時(shí)，試樣的初始含水量為22.08%，但試驗(yàn)結(jié)束時(shí)，試樣內(nèi)的最大含水量的實(shí)驗(yàn)值為29.2%，計(jì)算值為32.2%；最小值含水量的實(shí)驗(yàn)值為18.2%，計(jì)算值為18.4%，計(jì)算值與實(shí)驗(yàn)值非常接近，并且最大含水量位置的計(jì)算值和實(shí)驗(yàn)值也相當(dāng)吻合。這說明水-熱模型可以很好地模擬土體凍結(jié)過程中的水分遷移現(xiàn)象。

圖5為在不同凍結(jié)時(shí)刻(100、528、2 830 min)，試樣不同深度處的含冰量分布?？梢钥闯?，在凍結(jié)過程中試樣發(fā)生自上而下的凍結(jié)，試樣含冰量隨時(shí)間的變化也表現(xiàn)出凍結(jié)初期主要分布在土體淺層，隨著凍結(jié)的進(jìn)行，冰向土體內(nèi)部不斷增長，在凍結(jié)鋒面位置處含冰量達(dá)到最大。由于在實(shí)驗(yàn)中土體含冰量的變化無法進(jìn)行測量，在實(shí)驗(yàn)結(jié)束后測得的都是液態(tài)水的含量，所以對土柱在凍結(jié)過程中冰含量的分布和變化的分析是一個(gè)難題。水-熱耦合模型可以實(shí)現(xiàn)對土體在凍結(jié)過程中冰含量變化的監(jiān)測，為分析土體內(nèi)冰的分布和含量提供了方便。

圖4 實(shí)驗(yàn)結(jié)束時(shí)土樣含水量的計(jì)算值和實(shí)驗(yàn)值的比較Fig.4 Comparison between the experimental and simulated water content of soil sample at the end of experiment

圖5 不同時(shí)刻土體含冰量的計(jì)算值Fig.5 The calculated ice content of soil sample at different time

2 數(shù)值模擬

2.1 工程概況

新建吉圖琿客運(yùn)專線位于吉林省東部吉林市和延邊朝鮮族自治州境內(nèi)，線路全長360 km，設(shè)計(jì)時(shí)速250 km/h，其中，榆樹川隧道進(jìn)口位于延吉市安圖縣境內(nèi)附近，隧道全長2 210 m，最大埋深約160 m，如圖6所示。根據(jù)延吉市地區(qū)氣象資料和水文地質(zhì)條件，隧址年平均氣溫5.32 ℃，極端最高氣溫37.7 ℃，極端最低氣溫-29.2 ℃，最冷月平均氣溫-16.5 ℃，最大凍結(jié)深度為1.68 m。隧道進(jìn)口處大氣月平均溫度如圖7所示。

圖6 榆樹川隧道照片F(xiàn)ig.6 Photo of the Yushuchuan tunnel

圖7 隧道外大氣溫度Fig.7 Air temperature outside the tunnel

根據(jù)水文地質(zhì)條件，榆樹川隧道通過地層為強(qiáng)風(fēng)化、全風(fēng)化砂巖。隧道開挖后，地層之間的熱平衡被打破，隧道的進(jìn)口和出口處因受到大氣溫度的作用，在冷季圍巖中的水分發(fā)生凍結(jié)，從而使隧道襯砌發(fā)生凍脹破壞。

2.2 有限元模型

考慮到隧道縱向長度較長，計(jì)算斷面位于榆樹川隧道里程DK273+704處，該位置距離隧道進(jìn)口30 m，隧道內(nèi)空氣溫度可取為大氣溫度[15]。計(jì)算模型寬AD=BC=60 m，高AB=CD=80 m，隧道埋深15 m，計(jì)算模型如圖8所示，圍巖和襯砌的熱力學(xué)參數(shù)見表4。

圖8 隧道計(jì)算模型Fig.8 Numerical model of the tunnel

表4 圍巖和襯砌熱力學(xué)參數(shù)

2.3 邊界條件

溫度邊界：模型AB和CD為絕熱邊界；隧道襯砌邊界與隧道內(nèi)空氣對流換熱，對流換熱系數(shù)為12.5 W/(m2·℃)。模型底部BC的熱流密度為0.01 W/m2。

水分邊界：考慮到隧道頂部的植被會(huì)減弱降雨和降雪對圍巖含水量的影響，同時(shí)在初期支護(hù)與二次襯砌之間施做防水層，因此，模型AB、BC、CD和AD邊界均為不透水邊界。

3 計(jì)算結(jié)果與分析

3.1 無保溫材料隧道溫度場分布

在隧道運(yùn)營期間，隧道進(jìn)口暴露在空氣之中，進(jìn)口處圍巖的凍融狀態(tài)受空氣溫度的直接影響[20]。圖9為隧道運(yùn)營20年后，襯砌支護(hù)背面最低溫度的分布。可以看出：隧道圍巖溫度沿隧道軸線呈對稱分布，最低溫度的最小值出現(xiàn)在隧道墻頂DL和DR處，最低溫度的最大值出現(xiàn)在仰拱下部G處。因此，將隧道的墻頂處DL和DR處作為隧道發(fā)生凍結(jié)的最不利位置。圖10為隧道運(yùn)營20年內(nèi)最不利位置處溫度隨時(shí)間變化。由圖可以看出，在隧道運(yùn)營期，圍巖發(fā)生凍結(jié)的最不利位置處溫度隨時(shí)間呈周期性波動(dòng)，并且最不利位置處的年最低溫度隨時(shí)間不斷降低，在隧道運(yùn)營后第8年2月2日溫度降低至最低，并在隨后的年份內(nèi)年最低溫度無明顯變化。因此，將隧道運(yùn)營后第8年2月2日作為隧道內(nèi)最不利位置發(fā)生凍結(jié)的最不利時(shí)間。

圖11和圖12分別給出了最不利凍結(jié)時(shí)間時(shí)，無保溫層隧道圍巖溫度和水分分布。由圖可以看出，在最不利凍結(jié)時(shí)間，圍巖0 ℃等溫線呈環(huán)向分布，并且在隧道圍巖處徑向形成厚度約為2.845 m的凍結(jié)圈。圍巖水分發(fā)生了明顯地重分布。沿著凍結(jié)圈，圍巖的未凍水含量和含冰量變化劇烈，在凍結(jié)圈內(nèi)，靠近隧道襯砌支護(hù)的位置，圍巖含冰量急劇增大，并達(dá)到30%，遠(yuǎn)遠(yuǎn)超過了圍巖的初始含水量18%，這就是隧道發(fā)生凍脹破壞的原因。而在靠近0 ℃等溫線的位置處，圍巖未凍水含量為18%，含冰量為12%，圍巖的總含水量為30%，也超過了圍巖的初始含水量18%，這主要是由于在凍融循環(huán)作用下，圍巖內(nèi)的水分向凍結(jié)鋒面發(fā)生了遷移。因此，對于嚴(yán)寒地區(qū)隧道需施做相應(yīng)的保溫措施，使隧道襯砌不因圍巖發(fā)生凍結(jié)而破壞。

圖9 隧道運(yùn)營第20年初期支護(hù)背面最低溫度分布Fig.9 The minimum temperature at the boundary of preliminary lining of the selected section in the 20th year

圖10 隧道運(yùn)營20年內(nèi)最不利位置溫度隨時(shí)間變化Fig.10 Temperature variation at the unfavorable position of tunnel during 20 years

圖11 最不利時(shí)間時(shí)無保溫層隧道圍巖溫度分布Fig.11 Temperature distribution of the selected section of tunnel without thermal insulation layer at the unfavorable time

圖12 最不利時(shí)間時(shí)無保溫隧道圍巖水分分布Fig.12 Moisture distribution of the selected section of tunnel without thermal insulation layer at the unfavorable time

3.2 大氣溫度對凍結(jié)圈的影響

大氣溫度處于負(fù)溫是隧道圍巖內(nèi)水分發(fā)生凍結(jié)的原因，進(jìn)而引發(fā)隧道襯砌的凍脹破壞。所以，對不同溫度條件下，凍融圈變化規(guī)律的研究必不可少，而溫度的變化主要體現(xiàn)在空氣年平均溫度的變化方面。圖13是空氣年平均溫度分別為1.82、2.32、3.32、4.32、5.32 ℃時(shí)隧道的溫度分布?？梢钥闯?，空氣溫度對隧道溫度場的分布影響巨大，隨著空氣年平均溫度的增大，隧道各位置處的凍結(jié)深度顯著降低，以隧道墻頂為例，空氣年平均溫度為1.82 ℃時(shí)，墻頂處的最大凍結(jié)深度為4.13 m；隨著空氣年平均溫度升高至2.32 ℃，墻頂處的最大凍結(jié)深度減小為3.575 m；當(dāng)空氣年平均溫度升高至5.32 ℃，墻頂處的最大凍結(jié)深度減小為2.048 m。隧道最大凍結(jié)深度隨空氣年平均溫度的變化如圖14所示，對隧道最大凍結(jié)深度隨年平均溫度的變化擬合函數(shù)關(guān)系，得

圖13 不同空氣溫度條件下無保溫層隧道溫度分布Fig.13 Temperature distribution of the selected section of tunnel without thermal insulation layer under different atmospheric temperature

d=6.132 4Ta-0.651 4

(14)

式(14)中：d為隧道最大凍結(jié)深度，m；Ta為空氣年平均溫度，℃。

隧道內(nèi)不同位置處的空氣溫度往往與隧道外空氣溫度不同，如果按隧道外空氣溫度確定隧道內(nèi)不同深度處的最大凍結(jié)深度會(huì)存在誤差。運(yùn)用擬合得到的隧道最大凍結(jié)深度與空氣年平均溫度的變化關(guān)系[式(14)]，只需知道隧道內(nèi)不同深度處的年平均溫度，即可確定隧道內(nèi)不同深度處的最大凍結(jié)深度，可以很好地可以解決上述存在的問題，為隧道內(nèi)不同深度處最大凍結(jié)深度的確定提供依據(jù)。

圖14 隧道最大凍結(jié)深度隨空氣年平均溫度的變化Fig.14 Variation in the maximum frozen depth of tunnel under different atmospheric temperature

3.3 保溫材料對凍融圈的影響

在寒區(qū)隧道工程中，為了防止隧道發(fā)生凍脹破壞，常采用圍巖注漿，鋪設(shè)防水層，設(shè)置止水帶，設(shè)置保溫層等綜合措施，而作為嚴(yán)寒地區(qū)隧道防治凍害的有效措施，保溫層常常鋪設(shè)于初期支護(hù)與二次襯砌之間以達(dá)到防止圍巖發(fā)生凍脹的目的?，F(xiàn)采用上述水-熱耦合模型，對鋪設(shè)導(dǎo)熱系數(shù)為0.025 W/(m·℃)，厚度為5 cm保溫材料的隧道運(yùn)營8年后的溫度場進(jìn)行分析。

圖15為最不利凍結(jié)時(shí)間時(shí)有保溫層隧道圍巖溫度分布。由圖可以看出，鋪設(shè)保溫層后，在最不利凍結(jié)時(shí)間，圍巖已無0 ℃等溫線存在，這說明保溫層的鋪設(shè)使嚴(yán)寒地區(qū)隧道圍巖凍結(jié)狀況得到了明顯地改善，降低了圍巖水分凍結(jié)的可能性，避免了襯砌凍害的產(chǎn)生。

圖15 最不利時(shí)間時(shí)有保溫層隧道圍巖溫度分布Fig.15 Temperature distribution of the selected section of tunnel with thermal insulation layer at the unfavorable time

以隧道墻頂位置處，取二次襯砌表面、初期支護(hù)表面和圍巖表面作為關(guān)鍵點(diǎn)，分析保溫層對隧道保溫效果的影響。圖16為隧道運(yùn)營第8年，有、無保溫層隧道各關(guān)鍵點(diǎn)溫度隨時(shí)間的變化規(guī)律。由圖可以看出：隧道襯砌在無保溫層時(shí)，襯砌表面溫度受隧道內(nèi)空氣溫度的直接影響，呈正、負(fù)溫度交替變化，在冷季形成的低溫凍結(jié)圈在暖季全部融化，如此寒暖交替，凍融反復(fù)。而鋪設(shè)保溫層后，雖隧道二次襯砌溫度還隨時(shí)間呈正、負(fù)溫度交替變化，但初期支護(hù)和圍巖表面溫度全年都為正溫，從而避免了圍巖水分的凍結(jié)?？梢?，保溫層的鋪設(shè)對嚴(yán)寒地區(qū)隧道防凍有重要的意義。

圖16 運(yùn)營第8年無保溫層和有保溫層隧道關(guān)鍵點(diǎn)溫度變化Fig.16 Temperature variation at the classical position of tunnel without and with thermal insulation layer in the 8th year

影響保溫效果的一個(gè)重要因素是保溫材料的導(dǎo)熱系數(shù)?？紤]到目前寒區(qū)隧道采用的保溫材料的導(dǎo)熱系數(shù)一般在0.03 W/(m·℃)以下[21]，下面選取五個(gè)典型的導(dǎo)熱系數(shù)λ作趨勢分析。圖17是隧道分別鋪設(shè)導(dǎo)熱系數(shù)為0.008、0.010、0.012、0.015、0.017 W/(m·℃)的保溫材料(鋪設(shè)厚度均為10 cm)時(shí)，隧道運(yùn)行第7～10年時(shí)間內(nèi)，初期支護(hù)與保溫層接觸位置溫度變化規(guī)律?？梢钥闯?，該位置溫度變化幅度隨保溫材料導(dǎo)熱系數(shù)λ的增加而增加，在隧道表面敷設(shè)厚度為10 cm，導(dǎo)熱系數(shù)為0.012 W/(m·℃)的保溫材料后，已經(jīng)可以保證隧道不出現(xiàn)負(fù)溫的不利工況。

圖17 不同導(dǎo)熱系數(shù)條件下初期支護(hù)表面溫度變化規(guī)律Fig.17 Variation of temperature at the unfavorable position of preliminary lining in the condition of different thermal conductivity

表5為表面鋪設(shè)不同導(dǎo)熱系數(shù)的保溫材料，隧道運(yùn)行第7～10年時(shí)間內(nèi)各典型位置的最大凍深。由表可以看出，隧道墻頂?shù)膬錾罹S著保溫材料導(dǎo)熱系數(shù)的減小而減小，當(dāng)采用厚度為10 cm，導(dǎo)熱系數(shù)為0.017 W/(m·℃)的保溫材料時(shí)，隧道墻頂最大凍深為53.45 cm；當(dāng)采用導(dǎo)熱系數(shù)為0.015 W/(m·℃)的保溫材料后，墻頂處的凍深減小到了31.04 cm；減小了42%，而當(dāng)采用導(dǎo)熱系數(shù)為0.012 W/(m·℃)的保溫材料時(shí)，墻頂已經(jīng)不會(huì)出現(xiàn)負(fù)溫的情況，說明隧道防凍與保溫層導(dǎo)熱系數(shù)密切相關(guān)。

保溫材料對隧道的保溫效果不僅受到保溫材料類型的影響，還和保溫材料的厚度密切相關(guān)。圖18是隧道鋪設(shè)導(dǎo)熱系數(shù)為0.015 W/(m·℃)，厚度分別為8、10、13、14、15 cm的保溫材料時(shí)，隧道運(yùn)行第7～10年時(shí)間內(nèi)，初期支護(hù)與保溫層接觸位置溫度變化規(guī)律?？梢钥闯?，該位置處溫度變化振幅隨保溫層厚度的增加而減小，在隧道表面敷設(shè)導(dǎo)熱系數(shù)為0.015 W/(m·℃)厚度為13 cm的保溫材料已經(jīng)能保證隧道不出現(xiàn)負(fù)溫的不利工況。

圖18 不同厚度條件下初期支護(hù)表面溫度變化規(guī)律Fig.18 Variation of temperature at the unfavorable position of preliminary lining in the condition of different thickness

表6為表面鋪設(shè)不同厚度的保溫材料，隧道運(yùn)行第7～10年時(shí)間內(nèi)各典型位置的最大凍深。由表可以看出，隧道墻頂?shù)膬錾铍S著保溫材料厚度的增大而減小，當(dāng)保溫層厚度為8 cm時(shí)，隧道墻頂最大凍深為46.7 cm；當(dāng)保溫層厚度為10 cm時(shí)，凍深則減小到了31 cm，減小了33.6%；而當(dāng)保溫層厚度為13 cm時(shí)，該位置則已經(jīng)不會(huì)出現(xiàn)負(fù)溫的情況，說明隧道防凍與保溫層厚度密切相關(guān)。

表6 最大凍深與保溫層厚度的關(guān)系

4 結(jié)論

(1)本文建立的土體凍結(jié)過程中的水-熱耦合模型，考慮了水分遷移和冰-水相變潛熱對土體凍結(jié)的影響，并借助前人的實(shí)驗(yàn)結(jié)果，通過數(shù)值計(jì)算的方法驗(yàn)證耦合模型的準(zhǔn)確性，計(jì)算結(jié)果與前人實(shí)驗(yàn)結(jié)果十分吻合，驗(yàn)證了本文建立的水-熱耦合模型在研究土體凍結(jié)過程中的有效性。

(2)隧道的墻頂處是隧道發(fā)生凍結(jié)的最不利位置，最不利位置處溫度隨時(shí)間呈周期性波動(dòng)，并且最不利位置處的年最低溫度隨時(shí)間不斷降低，在隧道運(yùn)營后第8年2月2日溫度降低至最低，以此作為隧道內(nèi)最不利位置發(fā)生凍結(jié)的最不利時(shí)間。

(3)大氣溫度是隧道溫度場的分布決定因素，隨著大氣年平均溫度的增大，隧道各位置處的凍結(jié)深度顯著降低，隧道最大凍結(jié)深度與年平均溫度之間存在一定的關(guān)系。因此，可根據(jù)隧道內(nèi)不同深度處的年平均溫度確定隧道內(nèi)不同深度處的最大凍結(jié)深度。

(4)保溫層的鋪設(shè)使嚴(yán)寒地區(qū)隧道圍巖凍結(jié)狀況得到了明顯地改善，避免了襯砌凍害的產(chǎn)生。隧道的凍深隨保溫材料導(dǎo)熱系數(shù)的減小而減小，隨保溫層厚度的增加而減小。對于榆樹川隧道，敷設(shè)厚度為10 cm，導(dǎo)熱系數(shù)為0.012 W/(m·℃)的保溫材料，則墻頂已經(jīng)不會(huì)出現(xiàn)0 ℃以下的情況；敷設(shè)導(dǎo)熱系數(shù)為0.015 W/(m·℃)，厚度為13 cm的保溫材料己經(jīng)能保證該位置不出現(xiàn)負(fù)溫的情況。