朱艷峰，吳亞平，李文英

(1.廣州番禺職業(yè)技術(shù)學院，廣東廣州 511483;2.蘭州交通大學土木工程學院，甘肅蘭州 730070;3.天津鐵道職業(yè)技術(shù)學院，天津 300240)

0 引言

全球多年凍土約占陸地總面積的1/4，寒區(qū)則分布更為廣泛。我國陸地面積的43.5%為寒區(qū)［1］，其中多年凍土面積達215×104 km2［2］。凍土最大的特點就是水、熱、力性質(zhì)不穩(wěn)定，對溫度和外界因素影響極為敏感。寒區(qū)長大隧道建成運營后，高速列車在隧道中行駛，將會引起相應的空氣動力學效應，行車阻力增大，導致相當多的熱量散發(fā)到隧道環(huán)境中，且隨著行車速度的提高及行車密度的增大而加劇。此外，運營期間空調(diào)冷凝器散熱，制動電阻、閘瓦、牽引電機和隧道照明等設(shè)備的發(fā)熱，也會使隧道溫度不斷上升，相對濕度不斷下降。隨著隧道內(nèi)熱環(huán)境的日益惡化，寒區(qū)多年凍土隧道結(jié)構(gòu)周圍土體溫度升高，凍土天然的水熱收支平衡將遭到破壞，多年凍土的熱融沉降會引起土中應力的重新分布，從而導致地基的整體或不均勻下沉以及穩(wěn)定性下降［3－4］。襯砌漏水、路面(隧底)冒水等現(xiàn)象，將嚴重影響運營期間隧道結(jié)構(gòu)的安全性及壽命。

隨著長大隧道的建設(shè)以及列車速度的不斷提高，從20世紀80年代開始，國內(nèi)外學者開始關(guān)注特長隧道的溫升問題。H.Barrow［5］假設(shè)隧道壁面與列車間環(huán)狀區(qū)域為穩(wěn)態(tài)充分發(fā)展流，所有壁面的熱流密度為0，研究了英法海底隧道內(nèi)空氣溫度以正弦規(guī)律變化時對隧道巖土層溫度分布的影響;Arturo Baron等［6］應用非穩(wěn)態(tài)的一維可壓縮黏性Naviser-Stokes方程，通過有限體積法給出了高速列車在阻塞比為0.52的60 km特長隧道中運行時車周溫度的變化;麥繼婷等［7］將隧道內(nèi)空氣溫度滿足的方程與巖體溫度方程聯(lián)立求解，采用隧道壁面處涌水量以及隧道內(nèi)發(fā)熱量的估計值預測了秦嶺隧道運營較長時間后隧道內(nèi)溫度的變化情況;針對寒區(qū)隧道，賴遠明［8］等運用Galerkin法對寒區(qū)隧道溫度場、滲流場和應力場耦合的非線性問題進行了研究，求出了寒區(qū)圓形截面隧道溫度場的解析解;白國權(quán)［9］對高海拔嚴寒地區(qū)公路隧道溫度場分布規(guī)律及襯砌凍脹力進行了數(shù)值模擬研究;此外，2010年，朱艷峰等［10］建立了溫度應力耦合模型，對高速列車通過長隧道引起的熱、力效應進行了研究;2013年，龍垚［11］對昆侖山隧道溫度場進行了三維有限元分析和預測，以探討寒區(qū)隧道保溫隔熱層的設(shè)防長度?？傊壳拔覈鴮ｉT針對長期運營下的高速鐵路長大隧道進行隧道內(nèi)模擬分析和溫度預測的研究文獻還較少，關(guān)于隧道內(nèi)部濕熱環(huán)境的研究還基本上停留在對已修建完成的隧道進行現(xiàn)場觀測的程度上，在理論分析上仍顯不足。本文針對寒區(qū)多年凍土隧道圍巖溫度場具有導熱與對流換熱耦合邊界并伴有相變的非穩(wěn)態(tài)溫度場的特點，通過建立隧道內(nèi)流體及凍土體非穩(wěn)態(tài)溫度場的控制方程，將有限體積法與有限單元法綜合應用，首次系統(tǒng)研究了寒區(qū)高速鐵路長大隧道在長期運營后隧道空間以及隧道周圍凍土區(qū)的溫度變化，為寒區(qū)高速鐵路長大隧道隔熱保溫層的設(shè)置效果研究以及將來增加行車密度、提高列車運營速度提供理論依據(jù)。

1 計算模型及基本假設(shè)

寒區(qū)高速鐵路隧道是指所在地區(qū)常年平均氣溫低于0℃的高速鐵路隧道。本文以2種不同的寒區(qū)高速鐵路長大隧道作為分析模型，隧道埋深為100 m，長6 km，隧道截面面積分別為105.1 m2和65.1 m2，列車橫截面積At=13.664 m2，車頭形狀為45°傾角，阻塞比分別為0.13和0.21。對于特長隧道而言，隧道中部受外界氣溫變化的影響并不大，因此，本文在計算寒區(qū)高速鐵路長大隧道中部一點溫度場時，未考慮外界氣溫的影響。不考慮滲流的影響、由于地表溫度變化而產(chǎn)生的向隧道巖土層傳遞的熱量、涌水與圍巖間的對流換熱以及壁面處涌水的蒸發(fā)、凝結(jié)引起的熱量變化，只考慮通過隧道內(nèi)壁面向凍土區(qū)的徑向熱傳遞?？諝馀c隧道壁面的對流換熱系數(shù)α取15.0 W/(m2·℃)，混凝土的材料參數(shù)及多年凍土區(qū)隧道周圍巖土層的熱力學參數(shù)如表1和表2所示。

表1 混凝土的材料參數(shù)Table 1 Parameters of concrete

表2 凍土區(qū)圍巖熱力學參數(shù)Table 2 Thermodynamic parameters of surrounding rock in permafrost region

隧道內(nèi)空氣溫度的變化主要和空氣與隧道壁面間的對流換熱及空氣與設(shè)備均勻、局部發(fā)熱間熱量交換的2個熱傳遞有關(guān)。隧道中熱量的主要來源是高速列車運行時產(chǎn)生的列車阻力熱和空調(diào)熱［12］。

列車阻力熱是列車在行駛過程中全部機械能損失產(chǎn)生的熱量，行車阻力由機械阻力和氣動阻力2部分組成。對時速為250～300 km/h的流線型列車來說，氣動阻力占全部阻力的90%以上，因而在熱量計算中，忽略機械阻力的影響。氣動阻力主要依賴于列車在隧道中的阻塞率、隧道和列車的長度、車頭和車尾的形狀、空氣柱的存在、隧道橫向連接結(jié)構(gòu)、隧道壁面的粗糙度以及在隧道中是否有其他列車存在等［13］。對于隧道內(nèi)流體區(qū)域，本文采用了以下的基本假定:1)高速列車在隧道中運行時，將隧道內(nèi)的空氣看作是理想氣體，不計入空氣浮力;2)沿著隧道長度方向的隧道橫斷面不變;3)考慮列車為一個有熱源的實心體;4)除了沿著隧道軸向方向存在著風速外，隧道橫斷面上亦存在著水平方向和垂直方向的速度分布;5)列車在隧道內(nèi)以恒定的速度運行。

1.1 寒區(qū)高速鐵路隧道內(nèi)流體溫度場基本方程

高速列車在寒區(qū)高速鐵路隧道中運行時，隧道內(nèi)流場是黏性可壓縮的、存在邊界層分離的三維不定常紊流流動，且為充分發(fā)展的湍流。流體溫度場方程包括:

1)連續(xù)方程。

式中:ρ為空氣的密度;t為時間;u為速度矢量。

2)動量方程。

3)能量方程。

式中:e為氣體的內(nèi)能;T為氣體的溫度;λ為氣體的導熱系數(shù);τij為黏性流體應力張量。

4)此外，對于完全氣體，熱力學狀態(tài)方程為

式中R為摩爾氣體常數(shù)。

5)隧道內(nèi)氣流與隧道壁面之間的對流熱交換是一個包含自然對流和強迫對流的復雜對流過程。通常對于空氣的自然對流系數(shù)一般在5.50 W/(m2·K)以內(nèi)。考慮到高速鐵路隧道一般為深埋，而地表大氣溫度、濕度及氣壓的日變化和季節(jié)性變化對深層土壤溫度影響微小，因此深層土壤溫度按恒溫處理。此時，隧道內(nèi)壁面的邊界條件為第三類邊界條件，熱交換方程為

列車壁面設(shè)為絕熱邊界條件，?T/?n=0。

式中:Tc為巖石溫度;α為隧道內(nèi)空氣與隧道壁面的對流換熱系數(shù);T為隧道壁面處空氣的溫度。

1.2 寒區(qū)非穩(wěn)態(tài)溫度場的控制方程

寒區(qū)凍土中的傳熱過程是一種伴隨著冰水相變的復雜傳熱過程，包括:1)冰與水的相變換熱;2)固體骨架(顆粒)之間由于相互接觸的導熱過程;3)水的導熱和對流換熱過程;4)固體顆粒與未凍水之間的對流換熱過程;5)固體顆粒與空隙中氣體以及顆粒之間的輻射換熱過程。凍土隧道周圍土體溫度場的控制方程是伴有相變的非穩(wěn)態(tài)導熱微分方程。相變界面把計算區(qū)域G分成2個部分，即正凍區(qū)域Gf和融化區(qū)域Gu。2個區(qū)域有各自的溫度場和物性參數(shù)，分別滿足穩(wěn)態(tài)導熱方程及各自的邊界和初始條件，其控制方程為

在Gf內(nèi)

在Gu內(nèi)

2個區(qū)域的溫度場在相變界面s(t)上耦合，滿足溫度連續(xù)條件和能量守恒條件。

初始條件為Tu|t=0=T0，固定邊界上可能的第三類邊界條件為

式中:Tf，Cf，λf和Tu，Cu，λu分別為凍結(jié)區(qū)域及融化區(qū)域內(nèi)巖土的溫度(℃)、熱容量(j·m－3·℃－1)和導熱系數(shù)(W·m－1·℃ －1);Lγd(W－Wu)為含水巖土的相變潛熱(J/kg)。

相變潛熱為等溫等壓條件下系統(tǒng)由一種相態(tài)變?yōu)榱硪环N相態(tài)的焓變H。系統(tǒng)的內(nèi)能與系統(tǒng)內(nèi)的焓，比容和壓力的關(guān)系為

用焓法求解伴有相變的凍土體非穩(wěn)態(tài)溫度場控制微分方程為

2 寒區(qū)隧道溫度場

本文采用數(shù)值積分方法進行求解，對流體區(qū)域應用有限體積法，固體區(qū)域應用有限單元法。流體區(qū)域網(wǎng)格數(shù)量為365 976;固體區(qū)域內(nèi)邊緣為混凝土襯砌，厚度為0.4 m，劃分為49 200個四節(jié)點單元;襯砌外為凍土區(qū)域，劃分為457 200個四節(jié)點單元。計算模型及網(wǎng)格劃分如圖1所示。

圖1 計算模型及網(wǎng)格劃分Fig.1 Calculation model and grid division

2.1 隧道空間溫度場

本文寒區(qū)高速鐵路隧道算例參考昆侖山隧道，多年凍土上限為2.1～2.7 m，隧道最大埋深為106 m，隧道年平均氣溫設(shè)為－2.5℃。對于長大隧道而言，隧道中部通過活塞效應引起的自然空氣交換來通風的效果并不理想。在高速列車長期連續(xù)不斷的運營下，列車能量的消耗以及附屬設(shè)備的散熱疊加會越來越嚴重，長大隧道內(nèi)熱量的積聚會使隧道內(nèi)的氣溫逐年升高。

圖2表示阻塞比為0.13和0.21時，每h通過6列速度為500 km/h的高速列車在寒區(qū)長大隧道中連續(xù)運行1 d時，隧道空間中部一點溫度隨時間的變化。從圖2可以看出，隧道內(nèi)溫度在1 d中緩慢增長，24 h內(nèi)溫度變化約為1.7℃和0.62℃。

圖2 列車連續(xù)運行1 d后寒區(qū)隧道空間中部一點溫度變化Fig.2 Change of temperature in tunnels in permafrost regions after one-day train operation

圖3 表示高速列車連續(xù)運行1年后寒區(qū)長大隧道空間中部一點溫度的變化情況。顯然隨著時間的延續(xù)，長大隧道內(nèi)的溫度也在不斷升高，一年的連續(xù)運營后，阻塞比為0.13的特大斷面隧道內(nèi)溫度由初始的－2.5℃變化為接近0℃;阻塞比為0.21的大斷面隧道內(nèi)溫度由初始的0.62℃變化為接近0.64℃，變化幅度只有0.02℃。而且可以看到，特大斷面隧道內(nèi)的溫度在開始運營的前幾個月上升較快，半年以后，上升速度明顯減緩，但在1年內(nèi)溫度仍維持在0℃以下。

圖4表示高速列車連續(xù)運行25年后寒區(qū)長大隧道中部一點溫度的變化情況。從圖4可以直觀地看出，特大斷面隧道大約在連續(xù)運營2年后，溫度由負變正。10年間，溫度的變化幅度約為2.7℃。雖然數(shù)值不大，但對于寒區(qū)高溫凍土隧道而言，這將是非常重要的。10年間，阻塞比為0.21的大斷面隧道的溫度變化速度持續(xù)減緩，總變化幅度約為0.015℃，幾乎可以忽略不計。隧道運營至25年，特大斷面和大斷面隧道內(nèi)溫度將分別為0.28℃和0.66℃，50年后，隧道內(nèi)環(huán)境溫度預計分別為0.36℃和0.67℃。由圖1還可看出，特大斷面隧道在開始運營階段熱量變化的梯度小于大斷面隧道，隧道空間溫度場達到穩(wěn)態(tài)的過程也將會滯后于大斷面隧道。25年后，隧道溫度基本保持不變，但將長期維持在0℃以上，隧道內(nèi)溫度場以近似穩(wěn)定的狀態(tài)存在。

圖3 列車連續(xù)運行1年后寒區(qū)隧道空間中部一點溫度變化Fig.3 Change of temperature in tunnels in permafrost regions after one-year train operation

圖4 列車連續(xù)運行25年后寒區(qū)隧道空間中部一點溫度變化Fig.4 Change of temperature in tunnels in permafrost regions after long-term train operation

2.2 隧道周圍凍土區(qū)溫度場預測

對于寒區(qū)隧道工程而言，最為關(guān)注的就是隧道周圍凍土區(qū)的溫度問題。例如:青藏鐵路穿越550 km的多年凍土區(qū)，其中，昆侖山隧道處于高原腹地，多年凍土下限為100～110 m，隧道洞身基本完全處于多年凍土范圍內(nèi)。高溫凍土區(qū)接近0℃的地溫及持續(xù)不斷的熱積累是引起下伏多年凍土不斷融化的主要原因。低溫凍土區(qū)進入多年凍土的熱積累暫時以增高地溫耗熱為主，隨著地溫的增高，低溫凍土區(qū)也可能發(fā)生強烈的凍土融化［14－17］。

圖5和圖6表示阻塞比為0.13的特大斷面寒區(qū)隧道，在速度為500 km/h的高速列車連續(xù)不斷運營下，隧道周圍凍土區(qū)溫度隨時間的變化情況及不同時間凍土的融化范圍。

圖5 不同時間隧道周圍凍土區(qū)溫度場Fig.5 Temperature field of permafrost soil around tunnel at different time

從圖5和圖6可以看出，在給定的邊界及初始條件下，凍土區(qū)域在最初的1年內(nèi)，由于隧道空間溫度升高導致凍土融化的范圍非常小，僅限于襯砌周圍很薄一層。隨著時間的增長，融化速率穩(wěn)定增長，融化區(qū)域逐步擴大。5年以后，凍土融化速率雖逐步減慢，但融化帶寬度卻隨時間的推移不斷增加。本文定義了－0.2℃為融化溫度，從計算結(jié)果來看，隧道經(jīng)過10年運營后融化區(qū)域距襯砌深度約在2.5 m范圍內(nèi)。

圖7為2種不同阻塞比下的凍土區(qū)高速鐵路長大隧道由于列車運行導致隧道周圍凍土融化深度的比較。從圖7可以看出，大斷面隧道引起的凍土融化深度遠大于特大斷面隧道;因此，對于寒區(qū)多年凍土區(qū)長大隧道而言，采用特大斷面隧道可以降低隧道周圍凍土的熱融沉降。

圖6 不同時間隧道周圍凍土區(qū)融化范圍Fig.6 Thawing scope of permafrost soil around tunnel at different time

圖7 不同阻塞比下圍巖凍土融化深度Fig.7 Thawing depth of surrounding rock under different congestion ratio

3 結(jié)論與建議

通過算例可知，在寒區(qū)建設(shè)高速鐵路長大隧道，采用特大斷面隧道將是比較有利的;但無論阻塞比如何變化，運營10年后寒區(qū)高速鐵路長大隧道中部的溫度都已在0℃以上。50年后，寒區(qū)隧道內(nèi)環(huán)境溫度預計將分別達到0.36℃和0.64℃，隧道中部溫度將以近似穩(wěn)態(tài)的形式存在，對于寒區(qū)高溫凍土隧道而言，多年凍土的融化必將導致隧道基底結(jié)構(gòu)變形，對隧道結(jié)構(gòu)的穩(wěn)定性及高速列車的安全運行造成不利影響。實際上，隧道內(nèi)溫度的影響因素很多，寒區(qū)隧道內(nèi)溫度一年四季差別較大，每天早晚溫度差異也很大，保溫層的設(shè)置以及隨著時間的增長保溫功效的衰減，對長大隧道空間以及多年凍土的影響將是非常重要的，還有待于進一步研究。

［1］ 陳仁升，康爾泗，吳立宗，等.中國寒區(qū)分布探討［J］.冰川凍土，2005，27(4):469－475.(CHEN Rensheng，KANG Ersi，WU Lizong，et al.Cold regions in China［J］.Journal of Glaciology and Geocryology，2005，27(4):469 －475.(in Chinese))

［2］ 楊針娘.冰川水文學［M］.重慶:重慶出版社，2001.(YANG Zhenniang.Glacier hydrology［M］.Chongqing:Chongqing Press，2001.(in Chinese))

［3］ 程國棟，楊成松.青藏鐵路建設(shè)中的凍土力學問題［J］.力學與實踐，2006，28(3):1－8(CHENG Guodong，YANG Chengsong.Mechanics related with frozen ground in construction of Qinghai-Tibet railway［J］.Mechanics in Engineering，2006，28(3):1 －8.(in Chinese))

［4］ Henson D A，Bradbury WM S，Atkins WS.The aerodynamics of Channel tunnel trains［M］.［s.l.］:Elsevier Science Publishers Ltd.，1991:927 －956.

［5］ Barrow H.A simple analysis of flow and heat transfer in railway tunnels［J］.Heat and Fluid Flow，1987，8(2):119 －123.

［6］ Baron Arturo，Mossi Michele，Sibilla Stefano.The alleviation of the aerodynamic drag and wave effects of high-speed train in very long tunnels［J］.Journal of Wind Engineering and Industrial Aerodynamics，2001(89):365－401

［7］ 麥繼婷，陳春光.秦嶺特長隧道內(nèi)溫度預測［J］.西南交通大學學報，1998，33(2):153 －157.(MAI Jiting，CHEN Chunguang.Prediction of temperatures in the Qinling extra-long tunnel［J］. Journal of Southwest Jiaotong University，1998，33(2):153 －157.(in Chinese))

［8］ 賴遠明，吳紫汪，朱元林，等.寒區(qū)隧道溫度場、滲流場和應力場耦合問題的非線性分析［J］.巖土工程學報，1999，21(5):529 － 533.(LAI Yuanming，WU Ziwang，ZHU Yuanlin，et al.Nonlinear analyses for the couple problem of temperature，seepage and stress fields in cold region tunnels［J］. Chinese Journal of Geotechnical Engineering，1999，21(5):529 －533.(in Chinese))

［9］ 白國權(quán).高海拔嚴寒地區(qū)隧道溫度場分布規(guī)律及襯砌凍脹力數(shù)值模擬研究［D］.成都:西南交通大學，2006.(BAI Guoquan.Numerical simulation study on distribution of temperature field and frost force of lining for tunnel in high altitude and severe cold region［D］.Chengdu:Southwest Jiaotong University，2006.(in Chinese))

［10］ 朱艷峰，吳亞平，楊昌宇，等.高速列車通過長隧道引起的熱、力效應研究［J］.鐵道工程學報，2010(4):66－70.(ZHU Yanfeng，WU Yaping，YANG Changyu，et al.Study on thermal and mechanics effects of high-speed train running through long tunnel［J］.Journal of Railway Engineering Society，2010(4):66 －70.(in Chinese))

［11］ 龍垚.寒區(qū)長大隧道溫度場的三維分析及保溫措施研究［D］，蘭州:蘭州交通大學，2013.(LONG Yao.The three-dimensional analysis of temperature field of long tunnel in cold regions and study of heat preservation［D］.Lanzhou:Lanzhou Jiaotong University， 2013. (in Chinese))

［12］ Dappas Yves.Le tunnel sous la Manehe［J］.LeRail，1990，18(4):14－20.

［13］ Raghunathan Raghu S，Kim H-D，Setoguchi T.Aerodynamics of high speed railway train［J］.Progress in Aerospace Sciences，2002(38):469 －514.

［14］ 張旭芝.高原多年凍土涵洞溫度場及地基土凍融變形規(guī)律研究［D］.長沙:中南大學，2007.(ZHANG Xuzhi.Study on temperature fields and frost heave and thaw settlement of railway culvert in permafrost region［D］.Changsha:Central South University，2007.(in Chinese))

［15］ 朱正剛，錢江，李南生.寒區(qū)路基工程與多年凍土間相互作用問題研究進展［J］.結(jié)構(gòu)工程師，2004，20(6):43 － 48.(ZHU Zhenggang，QIAN Jiang，LI Nansheng.Recent progress in research on mutual-affections between roadbed and perennial frozen soils［J］. Structural Engineers，2004，20(6):43 －48.(in Chinese))

［16］ 盛煜，劉永智，張建明，等.青藏公路下伏多年凍土的融化分析［J］.冰川凍土，2003，25(1):43－48.(SHENG Yu，LIU Yongzhi，ZHANG Jianming，et al.Analysis on the thawing of permafrost underlying Qinghai-Tibet highway［J］.Journal of Glaciology and Geocryology，2003，25(1):43 －48.(in Chinese))

［17］ 胡元芳，王建宇.青藏鐵路昆侖山隧道凍脹壓力計算［J］.現(xiàn)代隧道枝術(shù)，2002，39(2):28 － 32.(HU Yuanfang， WANG Jianyu. Freezing caused swelling pressure on the lining of Kunlunshan tunnel on Qinghai-Tibet railway ［J］.Modern Tunnelling Technology，2002，39(2):28 －32.(in Chinese))