深季節(jié)凍土區(qū)涵洞熱狀況分析

汪江紅，牛永紅，慕青松，牛富俊，李蒙蒙

( 1.蘭州大學 西部災害與環(huán)境力學教育部重點實驗室；土木工程與力學學院，蘭州 730000;2.中國科學院寒區(qū)旱區(qū)環(huán)境與工程研究所 凍土工程國家重點實驗室，蘭州 730000)

深季節(jié)凍土區(qū)涵洞熱狀況分析

汪江紅1，牛永紅2，慕青松1，牛富俊2，李蒙蒙1

( 1.蘭州大學 西部災害與環(huán)境力學教育部重點實驗室；土木工程與力學學院，蘭州 730000;2.中國科學院寒區(qū)旱區(qū)環(huán)境與工程研究所 凍土工程國家重點實驗室，蘭州 730000)

為分析深季節(jié)凍土區(qū)涵洞對路基及其下地基熱狀況的影響，建立了涵洞溫度場計算模型，基于已有的試驗涵洞監(jiān)測數(shù)據(jù)，對比分析了所建模型的可靠性，再在該模型基礎上，通過改變材料參數(shù)和模型尺寸，進行不同情況下的數(shù)值計算。結果表明：模擬季節(jié)凍土區(qū)涵洞溫度場時，需要充分考慮熱對流效果，忽略填土熱對流效果是不合適的，考慮熱對流效果的模型在定性和定量上均與現(xiàn)場監(jiān)測情況一致；在季節(jié)性凍土區(qū)涵洞設計工作中，應該采取一些積極措施避免或減少對流傳熱對涵底凍深發(fā)展的影響；涵底最大凍結深度隨著含水量呈三段式的變化規(guī)律：隨其增大先減小，后減小速度變緩甚至凍深略有增加，最后又趨于減?。缓畢^(qū)涵洞內(nèi)徑的尺寸對涵底凍結深度方面的影響較大，而凈高對涵底凍結深度方面的影響很小，可以不考慮。

深季節(jié)凍土；涵洞；溫度分布；模擬計算

寒區(qū)涵洞建設是路基工程中必不可少的一部分，它的設置改善了凍土地區(qū)路基橫向排水條件，優(yōu)化了交通路線，但同時也改變了路基原有的水、熱平衡狀態(tài)，對路基的熱穩(wěn)定性有一定的影響。到目前為止，針對多年凍土區(qū)涵洞對路基熱穩(wěn)定性的影響，很多學者展開了研究，也取得了很多成果：俄羅斯專家Passek等綜合運用數(shù)值模擬和現(xiàn)場觀測方法，研究了大口徑波紋管涵在道路工程中的應用；Kunitake等[1]采用三維溫度場的有限元法對一混凝土箱涵進行了數(shù)值模擬；王凱等[2]針對青藏鐵路自身所處的特殊自然條件，對涵節(jié)預制、涵節(jié)運輸、涵洞基礎施工和涵節(jié)現(xiàn)場拼裝等技術作了技術改進，克服了現(xiàn)澆涵洞中所遇到的一些困難缺點；吳少海等[3]根據(jù)青藏公路涵洞的經(jīng)驗，結合青藏鐵路的特點，提出了適合多年凍土區(qū)的橋涵結構、支座類型、墩臺結構形式、基礎類型、混凝土材料及其施工方案；張學富等[4]通過對青藏鐵路某處涵洞現(xiàn)澆混凝土基礎水化熱進行數(shù)值分析，發(fā)現(xiàn)現(xiàn)澆混凝土水化熱對涵洞周圍凍土的熱狀況有較大影響，并計算分析了在涵底基礎下鋪設保溫層時凍土熱狀況；張旭芝[5]通過對青藏鐵路沱沱河試驗段兩座試驗涵洞進行地溫現(xiàn)場監(jiān)測和對觀測數(shù)據(jù)的分析，研究了適合青藏高原特殊施工環(huán)境的涵洞施工工藝、最佳施工季節(jié)、施工對多年凍土的影響以及沿涵洞軸向多年凍土上限的變化特征；孫文等[6]運用滲流理論和傳熱學理論，計算并討論了滲流對涵洞地基溫度場的影響以及涵洞對路基熱穩(wěn)定性的影響。由于早期季節(jié)凍土區(qū)工程經(jīng)驗不足，因此有關季節(jié)凍土區(qū)涵洞對路基熱穩(wěn)定性的影響研究不多。

哈爾濱-大連高速鐵路(簡稱“哈大高鐵”)是我國在嚴寒地區(qū)設計、建造的第一條鐵路客運專線，該線穿越東北地區(qū)中部山前平原，處于中-深季節(jié)凍土區(qū)，隨著它的建成通車，為鐵路建設提供了大量寶貴的資料和經(jīng)驗，基于此，許健等[7]運用數(shù)值方法模擬分析了季節(jié)凍土區(qū)保溫路基設計參數(shù)；Liu等[8]分析了季節(jié)性凍融路堤和地基不同的溫度變化過程，牛富俊等[9]對哈大高鐵某涵洞的路涵過渡段路基土體的凍結特征展開監(jiān)測，并分析了路涵過渡段路基及其下地基土的凍結特征和時空變化過程及其對路基穩(wěn)定性的影響，取得了一定的成果，但對于季節(jié)凍土區(qū)涵洞的研究還可以更加深入。

本文以哈大高鐵典型涵洞為原型，以熱學理論為基礎，建立了季節(jié)凍土區(qū)涵洞的溫度場計算模型，并與監(jiān)測數(shù)據(jù)進行對比分析該模型的可靠性，通過改變參數(shù)和尺寸，對不同模型進行數(shù)值計算，分析了冷季涵底基礎的對流效果、地基土含水量和涵洞尺寸對最大凍結深度發(fā)展的影響，對深季節(jié)凍土區(qū)涵洞工程建設和維護具有實際意義。

1 試驗涵洞概況

哈大客運專線涵洞防凍脹試驗研究試驗場地位于吉林省中西部的公主嶺市西北郊范家屯附近，海拔204 m，屬于溫帶半濕潤季風氣候，其特點為：冬冷夏熱，溫度、雨量、光照等季節(jié)性變化顯著，春季干旱多大風，回暖迅速；夏季熱而多雨；秋季溫暖多晴朗天氣；冬季漫長而寒冷。年平均氣溫5.6 ℃。據(jù)氣象資料及現(xiàn)場監(jiān)測的天然場地凍結深度變化數(shù)據(jù)顯示，最大凍結深度為1.56 m。

現(xiàn)場試驗涵洞設計里程為DK673+820，是1孔2.0 m的鋼筋混凝土箱型涵。該涵為排地表水而設，與線路正交。上下游均采用標準鋪砌，翼墻采用直徑φ12 mm、間距15 cm的護面鋼筋，基礎底地基處理方法同相鄰路基處理。上游側路基邊坡上設置簡易檢查梯一座。

試驗涵洞處的底層巖性特征為：

黏質黃土：厚5～16 m，分布于地表，軟塑為主，局部硬塑，Ⅲ級硬土。

粗砂：厚1.0～3.0 m，透鏡狀分布于黏質黃土之下，中密，飽和，Ⅰ級松土。

泥巖夾砂巖：成巖作用差，全風化泥巖巖芯呈泥柱狀、土狀，強風化泥巖巖芯呈塊狀、碎石狀、短柱狀，錘擊易碎，風化層厚度大于10 m，Ⅲ級硬土，弱風化，泥巖為主，泥質膠結，Ⅳ級軟石，為極軟巖，泥巖具有弱膨脹性。

地下水位埋深1～4 m，地下水對混凝土不具有侵蝕性，監(jiān)測資料顯示，地下水位在3到4月份保持在最低值，之后隨降水量的增加逐漸上升，到次年7～8月降水量最大時達到最高位，然后隨降水量的減少而緩慢下降，并保持水位在1 m以內(nèi)到11月底，之后，降水量很小，地下水位迅速下降，并于3月初達到最低值。

試驗涵洞段的路堤表層上是厚0.1 m的強化基床表層混凝土、其下是0.65 m的級配碎石砂礫土、1.00 m的非凍脹性A/B組填料、1.30 m的普通A/B組填料、基床以下路堤是1.66 m的普通A/B/C組填料，再下面是0.50 m的褥墊層和黏質黃土地基。路涵過渡段路基采用的是一段式過渡，涵洞高2.50 m，涵側混凝土壁厚0.22 m，涵頂和涵底厚0.28 m，涵身以下設置了厚1.02 m的涵底地基，總厚度達1.30 m，各構件組成圖如圖1所示。

圖1 涵洞各組成構件位置示意圖

2 計算模型

目前有關凍土區(qū)溫度場數(shù)值模型的研究比較多[10-14]，本文采用了一種較成熟的數(shù)值模擬方法，旨在研究深季節(jié)凍土區(qū)涵洞模型參數(shù)的選取以及分析各參數(shù)對其溫度場的影響。

2.1 理論模型

控制微分方程及其有限元公式

只考慮熱傳導和相變的情況下，溫度場的數(shù)學模型方程為[15]

式中：C、λ、ρd分別為模型中各材料的容積熱容量、導熱系數(shù)和干容重，Wi為含冰量；L為冰、水的相變潛熱。

在凍結過程中，未凍水含量是溫度的函數(shù)[16]，不考慮水分遷移，則

將式(2)帶入式(1)，簡化為

假設劇烈相變發(fā)生在(Tm-ΔT,Tm)范圍內(nèi)，且在這范圍內(nèi)，體積熱容量C及導熱系數(shù)λ均隨著溫度線性變化[9-10]。則

凍融過程的熱傳導是一個二維變系數(shù)非線性問題，采用具有迭代的有限元法，并通過減少時間步長提高精度，計算模型的有限元公式導出為

式中：Δt為時間步長；T0、T1分別表示每個計算步長開始之前和結束時溫度列向量；K表示導熱系數(shù)矩陣；M表示總變溫矩陣；Q1與邊界條件有關的熱流量矢量列陣。

采用Geo-Studio軟件中TEMP/W(地下熱傳遞分析軟件)進行模擬計算，Geo-Studio是一套適用于巖土工程和巖土環(huán)境模擬計算的仿真軟件。在TEMP/W軟件中，不再需要編制繁雜的求解器，它含有全面的地熱模型，用戶可以根據(jù)自己的需要設定，能對水轉化成冰和冰轉化成水的過程中釋放和吸收的相變潛熱問題進行計算分析，并且將其表示為未凍水含量的函數(shù)。TEMP/W是一種對各種不同地面條件進行分析的強大軟件。

2.2 模型計算區(qū)域及土體熱物理參數(shù)

試驗涵洞在幾何形狀上是關于橫截面和中截面對稱的，忽略陰陽坡效應等一些不對稱因素，將涵洞及其過渡段簡化為一半來計算。根據(jù)《高速鐵路設計規(guī)范(試行)》(2009年)，綜合考慮現(xiàn)場實際情況(圖1)，將計算區(qū)域延伸到天然地表以下14 m，沿路基方向取25 m，模型示意簡圖如圖2。

圖2 計算模型示意簡圖

2.3 初始條件及邊界條件的選擇

由于涵頂在2009年底做了堆載預壓處理，于2010年6月底卸除完畢，并從2010年7月開始以每周一次的頻率進行正式的地溫監(jiān)測，所以以2010年7月10日這一天的實測溫度場數(shù)值作為初始溫度場。

涵洞及其過渡段溫度場初始溫度條件為：T(t)|t=0=T0。

監(jiān)測結果顯示涵洞處深部地溫為8 ℃左右，且年變幅只有0.1 ℃，因此取模型的下部邊界條件為：T=8 ℃。

對于上邊界和涵洞內(nèi)壁溫度邊界條件，為了使得計算模型更加符合實際情況，減小由于邊界條件的選擇而產(chǎn)生的誤差，擬采取實測溫度(如圖3)進行計算。

圖3 溫度邊界條件

涵洞內(nèi)壁沒有太陽直射，取室外大氣溫度作為計算溫度，室外溫度理論上應該是一天內(nèi)各個時間段氣溫的加權平均值，采用平均溫度之半加上最低溫度之半所得的值作為涵洞內(nèi)壁溫度邊界條件的計算數(shù)值，對于上邊界混凝土路面的溫度邊界條件，采用地表的實測溫度值作為上邊界的溫度邊界條件。

3 參數(shù)選取及計算結果分析

3.1 參數(shù)的選取

為了方便模型的建立，引入如下基本假定：

1) 假定在每次的模擬過程中涵洞及路基填料無水分補給和排水作用，即忽略水分遷移，這樣可以將含水量固定為一個常量，同時也是考慮到試驗涵洞現(xiàn)行的排水措施的有效性。

2) 假定填料及地基土介質為剛性體，在凍結過程中土體不發(fā)生變形。

3) 各區(qū)域的材料均勻且各向同性。

4)考慮熱傳導和相變作用，忽略填料內(nèi)部熱對流的影響。

5)土層熱物理參數(shù)的選擇，由于溫度對填料的體積熱容量影響比較小，為了簡便起見，計算過程中，土體在未相變區(qū)和凍結相變區(qū)的體積熱容量只取其融、凍兩種狀態(tài)的數(shù)值，即融土和凍土，同樣，忽略溫度對導熱系數(shù)的影響，在未相變區(qū)和凍結相變區(qū)土體的導熱系數(shù)只考慮凍融狀態(tài)[17-18]。

2011年春季以后，該處涵洞進行了修補操作，在涵底加鋪了20 cm的混凝土層等，考慮到水化熱等影響較大，本文所要研究的主要內(nèi)容是圍繞最大凍結深度等內(nèi)容展開，而4月份以后氣溫已經(jīng)升到0 ℃以上，所以所模擬的時間區(qū)間為2010-07-10到2011-04-02。

根據(jù)現(xiàn)場實測資料和有關參考資料的研究[7,18]，土層的熱物理參數(shù)見表1。

表1 模擬計算中各土層熱物理參數(shù)

圖4 涵下不同深度探頭處實測值與計算值的比較

為了驗證計算模型的可行性，取涵底下0.3、0.8、1.3、1.8和2.3 m各位置上的溫度變化情況(如圖4中的模擬值1)，與實測結果分別進行對比。

對涵底以下各探頭處實測值與計算值進行對比，可以看到，按表1所示的熱物理參數(shù)計算，暖季涵下溫度一致性比較好，但是冷季溫度遠高于監(jiān)測數(shù)據(jù)，從而由此得出的凍結深度也遠小于監(jiān)測數(shù)據(jù)。考慮到涵底基礎是由片石混凝土或漿砌片石鋪砌，涵底在暖季熱傳導起主流作用，在冷季熱對流會對其產(chǎn)生明顯影響，所以在模擬過程中忽略冷季熱對流效果會導致很大誤差，假定4)并不適用。因此假設涵底基礎冷季當量導熱系數(shù)為暖季的2倍(具體熱物理參數(shù)如表2所示)，并進行計算對比，對于不同涵底基礎熱對流效果對涵底最大凍深的影響在后面進行討論。

表2 涵底基礎熱物理參數(shù)

圖4中模擬值2是假設涵底基礎冷季當量導熱系數(shù)為暖季的2倍時涵底下不同深度探頭的溫度計算值，對比實測值可以看到，在模擬時間段內(nèi)，各探頭點的模擬計算結果與實測數(shù)據(jù)吻合的相當好，誤差都比較小，該模擬結果不僅在溫度上升下降的整體趨勢是與實測值相符合，在對涵洞進行定量計算都有一定的可靠性。

為了進一步分析模擬結果與實測結果之間的差值隨深度的變化規(guī)律，取2010-12-04和2011-02-26涵底不同位置溫度示意圖(圖5)，可以看出進行參數(shù)修正以后，模擬精度得到提高，尤其在凍結層附近的溫度擬合精度較高。

圖5 涵底地溫隨深度分布對比示意圖

3.2 涵洞熱狀況分析

根據(jù)第二次模擬結果，對涵洞的熱狀況進行分析。

如圖6所示，取路基下0.65和3.25 m距涵洞中心不同距離處的溫度進行對比，可以看出涵洞對涵周路基溫度分布有顯著影響，加大了路基溫度隨氣溫變化而改變的幅度，對氣溫改變的響應更為迅速，而且越接近涵壁，這種影響更為劇烈。

圖6 路基下距涵洞中心不同距離處的溫度變化圖

如圖7所示的涵頂填料中不同位置溫度時程曲線圖，可以看到在高溫段和低溫段，涵頂溫度發(fā)展較均衡，溫差較小，而且由于太陽輻射的作用，從上往下溫度降低，另外在降溫和升溫過程中，分別表現(xiàn)出上下低于中心和上下高于中心的特點，雙向凍結和雙向融化明顯。

為便于對涵洞下凍結深度的發(fā)展過程有個直觀了解，根據(jù)有限元模型的計算結果，繪制了涵下基礎中不同深度溫度時程曲線(如圖8)。其中橫坐標表示時間歷程，縱坐標表示從涵洞底向下的深度。

圖7 涵頂填料中不同高度溫度時程曲線

圖8 涵洞基礎中心下不同深度溫度時程圖

計算模型溫度邊界條件采用的是實測溫度，在凍結初期，由于溫度的擾動產(chǎn)生了一段波動，但隨著溫度整體趨勢趨于下降，涵洞下土體正式進入凍結階段，正式凍結開始后，涵底1 m范圍內(nèi)溫度下降迅速，在1周內(nèi)即整體下降到0 ℃以下，然后趨于平緩，以一定的速率向下發(fā)展，經(jīng)歷大概100 d，達到最大凍結深度；之后進入融化階段，凍深慢慢減小，仔細觀察包絡線，會發(fā)現(xiàn)融化結束前0 ℃線幾乎是垂直的，表明涵洞中心下溫度梯度很小。

如圖9所示的涵周路基凍結圈發(fā)展歷程，在凍結期內(nèi)，涵洞周圍路基內(nèi)形成的凍結圈呈“倒鐘”形。涵周填土首先形成凍結圈，然后涵頂雙向凍結，涵底單向凍結，第4 d涵頂填料就整體凍結，形成“倒鐘”形凍結圈，再經(jīng)過大約90 d，于2011-02-22形成最大凍結圈；再經(jīng)過40 d左右，涵頂和大部分涵下凍層消失，一個月以后，于2011-05-06凍層徹底消失，這些與現(xiàn)場實測數(shù)據(jù)體現(xiàn)的凍結圈發(fā)展歷程基本吻合。涵底中心凍結層最大厚度為2.26 m。

圖9 涵周路基凍結圈的發(fā)展歷程圖

綜上對比分析計算結果和實測數(shù)據(jù)，可以看出，這種根據(jù)試驗涵洞實際鋪砌情況及實測溫度建立的涵洞段溫度場模型，不僅在定性分析涵洞熱狀態(tài)變化趨勢上與現(xiàn)場實際情況有較好的吻合，而且在定量計算上都有一定的可靠性。這說明在溫度場模型建立過程中，應該根據(jù)填料熱對流實際情況，確定其冷季當量導熱系數(shù)，由此建立模型計算結果才能夠準確的反應溫度狀況，對其它涵洞工程的設計具有一定的參考價值。

3.3 涵底基礎熱對流效果對涵底最大凍深的影響

涵底基礎是由片石混凝土或漿砌片石鋪砌的，在前面的分析討論中看到，涵底基礎的冷季對流效果明顯，加深了涵底最大凍結深度。

在多年凍土區(qū)，碎石通風路基是一種積極保護多年凍土的工程措施，工作原理是：在寒冷季節(jié)，冷空氣有較大的密度，在自重和風的作用下使片石間隙中的熱空氣上升，冷空氣下降并進入地基；而在溫暖的季節(jié)，熱空氣密度小，很難進入地基，類似于熱開關效應。碎石路基在多年凍土區(qū)對于保護地基土凍結狀態(tài)具有優(yōu)勢，被國內(nèi)外廣泛應用于凍土區(qū)房屋建筑，相關資料指出，某試驗中塊石層的有效導熱系數(shù)在冬季是夏季的12.2倍[19]。

但在季節(jié)性凍土區(qū)，由于防止凍脹災害關鍵在于減小凍結深度，減小凍脹量，而非多年凍土區(qū)的保護凍土層，防止凍土層的退化。在上面建立的涵洞溫度場計算模型中，考慮到對流作用，對涵底基礎取2倍于暖季導熱系數(shù)的當量導熱系數(shù)為其冷季熱參數(shù)值，與實測數(shù)據(jù)的比較也能看出，這種處理辦法是符合客觀規(guī)律的。

為了討論季節(jié)性凍土區(qū)涵底碎石對流效果對涵底最大凍結深度的影響，采用上面建立的溫度場模型，取冷季當量導熱系數(shù)分別為暖季的1～12倍，其它參數(shù)不變進行計算。

如圖10所示最大凍結深度與冷、暖季涵底基礎當量導熱系數(shù)比值之間的關系，可以用下面的函數(shù)公式進行較好的擬合：

隨著涵底碎石冷、暖季當量導熱系數(shù)比值的增加，涵底凍結深度呈增函數(shù)形式變化，而且增長先快后慢，當λc/λw=2時，最大凍結深度較沒有熱對流時增加了0.38 m，從而看出，對于深季節(jié)凍土區(qū)，碎石層的對流傳熱對防止凍脹災害是極其有害的，所以在季節(jié)性寒區(qū)涵洞設計工作中，應該避免一般涵洞設計中，涵底基礎統(tǒng)一用片石混凝土或漿砌片石，應該采用新工藝，減少碎石中的孔隙，或者減少孔隙連通性，減少對流傳熱對涵底凍深發(fā)展的影響，使得涵底碎石基礎在穩(wěn)固路基的同時，還能夠起到很好的保溫效果。

圖10 當W=25.0%時，最大凍結深度隨λc/λw的變化趨勢

3.4 地基土含水量對涵底最大凍深的影響

為了研究地基土含水量對最大凍深的影響，分別取不同含水量的地基土所對應的熱學參數(shù)(如圖11、12所示)，帶到上面有限元模型中進行計算。

圖11 λu、λf隨含水量的變化規(guī)律

圖12 Cu、Cf隨含水量的變化規(guī)律

不同體積含水量地基土涵底最大凍結深度變化趨勢如圖13所示，隨著含水量的增加，最大凍結深度的變化規(guī)律大致可分為3段：1)當含水量W<17.5%，凍深隨著含水量的增大而減少；2)當17.5%27.5%時凍深隨著含水量的增加而減小。

圖13 最大凍結深度隨體積含水量變化圖

分析圖11、12中融、凍地基土的導熱系數(shù)、體積熱容量與含水量之間的關系，可以看到，融土導熱系數(shù)隨含水量增加的變化規(guī)律是前面增加的快，后面增加的慢，凍土的變化曲線分3個階段，兩頭緩中間陡；融土的容積熱容量隨含水量呈直線增大，凍土的容積熱容量隨含水量呈折線增大，整體上凍土容積熱容量隨含水量增大的斜率較融土的緩[18]。據(jù)此可以分析最大凍結深度隨含水量變化規(guī)律。

第1階段，干燥即結合水未達到最大值時，這段時間內(nèi)，凍、融土的導熱系數(shù)與體積熱容量基本上都呈線性變化趨勢，同時，隨著含水量的增加，潛熱變大，從而導致涵底最大凍結深度不斷減??；第2階段，凍、融土的體積熱容量依然呈線性變化趨勢，融土的導熱系數(shù)隨含水量的變化依然保持第1階段的增速，但凍土的導熱系數(shù)隨著含水量的變化快速增大，從而導致在冷季，地基土中能量大量流失，從而導致涵下最大凍結深度減小速度減緩，隨著含水量的持續(xù)增加，凍深稍稍有些加深；第3階段，隨著含水量的增加，凍、融土的導熱系數(shù)和體積熱容量的變化規(guī)律又趨于一致，由于潛熱的增大，涵底最大凍結深度又開始不斷減少。這與凍土導熱系數(shù)隨含水量變化規(guī)律分3個階段是一致的。整體上看，涵底最大凍結深度隨著含水量的增加呈減小趨勢。

在實測中，涵底最大凍深為天然地層的兩倍左右，究其原因，通過上面的計算討論，最有可能是有兩個方面：1)冷季涵底基礎熱對流效果明顯，當量導熱系數(shù)遠大于天然地層，散熱嚴重；2)由于天然地層的含水量較涵底碎石基礎的要高，根據(jù)最大凍結深度隨含水量變化規(guī)律，天然地層的凍深較淺顯而易見。涵底凍深較大，一方面體現(xiàn)出涵洞及路基內(nèi)采取的隔水措施是有效的；另一方面，也表明當區(qū)域地基土是凍脹敏感性土，且隔、排水措施完善，人為凍深也需要考慮的相對較深。

3.5 涵洞尺寸對涵底凍深的影響

試驗涵洞為鋼筋混凝土箱型涵，為了排地表水而設，根據(jù)《鐵路橋涵設計基本規(guī)范》[20]，泄水涵洞的孔徑不宜小于2.0 m，涵洞頂至軌底的填方厚度不小于1.2 m，并且對涵洞凈高與涵洞長度做出了規(guī)定。為了研究涵洞幾何尺寸對季節(jié)凍土區(qū)涵底最大凍結深度的影響，利用上面的模型，針對不同尺寸的涵洞計算其涵底最大凍結深度，為了滿足排水需求，而且在規(guī)范要求內(nèi)，將涵洞孔徑和凈高均控制在2.00～2.95 m之間，同時計算過程中不考慮強度、穩(wěn)定性等要求，僅僅討論凈高、內(nèi)徑的改變對涵底凍深的影響，求得最合適的涵洞尺寸。計算結果如圖14所示。

圖14 不同尺寸下涵底最大凍結深度值的對比

從圖14可以看出，箱涵的凈高對涵底最大凍結深度的影響很小，當內(nèi)徑恒定為2 m，箱涵的凈高從2 m增加到2.95 m時，最大凍結深度僅增加了4.5 cm，增加率僅為2.22%；相較于凈高，內(nèi)徑對涵底最大凍結深度的影響相當明顯，隨著內(nèi)徑的增大，最大凍結深度呈線性加深，當凈高恒定為2.95 m，內(nèi)徑從2 m增加到2.95 m時，最大凍結深度增加了0.251 m，增加率為12.23%。所以在寒區(qū)涵洞的設計過程中，在滿足強度、穩(wěn)定性、功能性等一系列要求外，尺寸方面優(yōu)先考慮內(nèi)徑，取允許的最小內(nèi)徑，由于凈高對涵底凍結深度方面的影響很小，可以不考慮。

4 結 論

1) 根據(jù)試驗涵洞實際鋪砌情況及現(xiàn)場氣象資料，建立了考慮相變潛熱的涵洞段溫度場模型，對比分析計算結果和實測數(shù)據(jù)，表明考慮冷季熱對流效果后，計算結果不僅在定性分析涵洞熱狀態(tài)變化趨勢上與現(xiàn)場實際情況有較好的吻合，而且在定量計算上都有一定的可靠性，能夠準確的反應溫度狀況，對季節(jié)凍土區(qū)其它涵洞工程設計有一定的參考價值。

2) 在多年凍土區(qū)，忽略對流傳熱的計算是基于一種不利的情況進行計算，對于工程安全穩(wěn)定趨于保守的，但在季節(jié)凍土區(qū)，忽略對流傳熱的影響相當于是基于一種有利的情況進行計算，得到的人為凍結深度要淺于實際人為凍結深度，由此指導設計的涵洞工程就不能達到預定的安全系數(shù)。

3) 在深季節(jié)性凍土區(qū)涵洞設計工作中，冷季涵底基礎對流傳熱的影響不容忽視，應該避免一般的涵洞設計中，涵底基礎統(tǒng)一用片石混凝土或漿砌片石；應該采用新工藝，減少碎石中的孔隙，或者減少孔隙連通性，從而減小對流傳熱對涵底凍深發(fā)展的影響，使得涵底碎石基礎在穩(wěn)固路基的同時，還能夠起到很好的保溫效果。

4) 涵底最大凍結深度隨著含水量變化規(guī)律呈三段：當含水量W<17.5%，凍深隨著含水量的增大而減少；當17.5%27.5%時凍深隨著含水量的增加而減小。涵底最大凍結深度整體上隨著含水量的增加呈減小趨勢。據(jù)此，在隔、排水措施完善的區(qū)域，人為凍深相較更深，需要給予更多注意。

5) 研究了寒區(qū)涵洞尺寸對涵底凍深的影響，分別計算了不同凈高和不同內(nèi)徑下涵底最大凍結深度值，通過對比可以看出其內(nèi)徑對涵底凍結深度方面的影響很大，在滿足強度、穩(wěn)定性和功能性等一系列要求后應該優(yōu)先考慮，而凈高對涵底凍結深度方面的影響很小，可以不考慮。

[1] Kunitake M,Ogata H,Nakazawa T,et al.Temperature distribution in a concrete box culvert using 3-D FEM [C]//Proceedings of the 2nd Asian-Pacific Conference on Computational Mechanics，1993:215-228.

[2] 王凱,羅廣宇,羅英,等.高含冰量凍土地區(qū)拼裝涵洞施工技術[J].長安大學學報:建筑與環(huán)境科學版,2003,20(1):30-32.

Wang K,Luo G Y,Luo Y,et al.The construction technique of precast installed culvert at high-ice content frozen soil [J].Journal of Chang’an University:Architecture and Environment Science Edition,2003,20(1):30-32.

[3] 吳少海.青藏鐵路多年凍土區(qū)橋涵設計與施工[J].冰川凍土,2003,25(Sup1):59-68.

Wu S H.Design and construction of bridges and culverts of the Qinghai-Tibet Railway in permafrost regions [J].Journal of Glaciology and Geocryology,2003,25(Sup1):59-68.

[4] 張學富,賴遠明,喻文兵,等.青藏鐵路涵洞現(xiàn)澆混凝土基礎水化熱的數(shù)值分析[J].冰川凍土,2004,26(1):106-111.

Zhang X F,Lai Y M,Yu W B,et al.Numerical analysis for the hydration heat of cast-in-situ concrete foundation of the culvert on the Qinghai-Tibet Railway [J].Journal of Glaciology and Geocryology,2004,26(1):106-111.

[5] 張旭芝,王星華.凍土鐵路涵洞施工對地基土地溫的影響[J].中國鐵道科學,2007,28(4):19-24.

Zhang X Z,Wang X H.Impact of railway culvert construction on the ground temperature of foundation soil in permafrost regions [J].China Railway Science,2007,28(4):19-24.

[6] 孫文.多年凍土區(qū)滲流對路基熱穩(wěn)定性的影響分析[J]. 城市道橋與防洪,2012,29(2):108-111.

Sun W.Analysis on influence of seepage flow in permafrost area on thermal stability of roadbed [J].Urban Roads Bridges & Flood Control,2012,29(2):108-111.

[7] 許健,牛富俊,牛永紅,等.季節(jié)凍土區(qū)保溫路基設計參數(shù)[J].土木建筑與環(huán)境工程,2009,31(3):83-89.

Xu J,Niu F J,Niu Y H,et al.The design parameters of roadbed with insulation in seasonal frozen ground [J].Journal of Civil，ArchitectvraI & Environmental Engineering,2009,31(3):83-89.

[8] Liu H,Niu F J,Niu Y H,et al.Experimental and numerical investigation on temperature characteristics of high-speed railway’s embankment in seasonal frozen regions [J].Cold Regions Science and Technology,2012,81:55-64.

[9] 牛富俊,劉華,牛永紅.高速鐵路路涵過渡段路基凍結特征試驗研究[J].巖石力學與工程學報,2014,33(3):639-646.

Niu F J,Liu H,Niu Y H.Study on the freeze characteristic in roadbed-culvert transition section along a high-speed railway [J].Chinese Journal of Rock Mechanics and Engineering,2014,33(3):639-646.

[10] Lunardini V J.Heat transfer with freezing and thawing [M].Amsterdam:Elsevier,1991.

[11] Miao T D,Guo L,Niu Y H.et al.Modeling on coupled heat and moisture transfer in freezing soil using mixture theory [J].Science in China:Series D,1999,42(Sup1):9-16.

[12] Wang Z L,Fu Q,Jiang Q X,et al.Numerical simulation of water-heat coupled movements in seasonal frozen soil [J].Mathematical and Computer Modelling,2011,54(3/4):970-975.

[13] Zhang K,Li D Q,Niu F J,et al.Cooling effects study on ventilated embankments under the influence of the temperature differences between the sunny slopes and the shady slopes [J].Cold Regions Science and Technology,2011,65(2):226-233.

[14] Wu D,Jin L,Peng J B,et al.The thermal budget evaluation of the two-phase closed thermosyphon embankment of the Qinghai-Tibet Highway in permafrost regions [J].Cold Regions Science and Technology,2014,103:115-122.

[15] 安維東. 凍土的溫度水分應力及其相互作用[M].蘭州:蘭州大學出版社,1990:21-54.

[16] Lai Y M,Zhang S J,Mi L.Effect of climatic warming on the temperature fields of embankments in cold regions and a countermeasure [J].Numerical Heat Transfer Part A:Application,2004,45(2):191-210.

[17] 丑亞玲,盛煜,馬巍.多年凍土區(qū)道路邊坡熱狀況差異對多年凍土融化形態(tài)的影響[J].冰川凍土,2007,29(6):977-985.

Chou Y L,Sheng Y,Ma W.Effect of the differences of thermal regime within roadbed slopes on thawing morphology in permafrost regions [J].Journal of Glaciology and Geocryology,2007,29(6):977-985.

[18] 徐斅祖,王家澄,張立新.凍土物理學[M].2版.北京:科學出版社,2001:75-86.

[19] 葛建軍.青藏鐵路高原多年凍土區(qū)片石通風路基[J].冰川凍土,2003,25(Sup1):52-53.

Ge J J.The ventilative sheet-stone roadbed of the Qinghai-Tibet Railway in plateau permafrost regions [J].Journal of Glaciology and Geocryology,2003,25(Sup1):52-53.

[20] 中華人民共和國鐵道部.TB 10002.1—2005 鐵路橋涵設計基本規(guī)范[S].北京：中國鐵道出版社，2005.

(編輯 王秀玲)

ThermalConditionsofCulvertinDeepSeasonallyFrozenRegion

WangJianghong1,NiuYonghong2,MuQingsong1,NiuFujun2,LiMengmeng1

( 1. Key Laboratory of Mechanics on Disaster and Environment in Western China, the Ministry of Education of China;College of Civil Engineer and Mechanics, Lanzhou University, Lanzhou 730000, P.R.China;2. State Key Laboratory of Frozen Soil Engineering, Cold and Arid Regions Environmental and Engineering Research Institute, CAS, Lanzhou 730000, P.R.China)

In order to explore the impact of culvert on the roadbed-culvert transition section in deep seasonally frozen region, a finite element model was established to simulate culvert’s temperature fields considering latent heat, and judged the reliability through comparing the simulation results with the existing monitoring data. Models of different parameters or different sizes were stimulated to investigate the influence on the maximum frozen depth. The results indicate that: 1) The measured temperatures and simulated temperatures agree well which demonstrated the accuracy of the numerical model. 2) In the design of culvert in seasonally frozen region, foundation convection under the culvert in cold season could not be ignored. Some active measures should be taken to protect the subsoil from the effect of convective heat transfer. 3) The culvert maximum frozen depth changing with the amount of water content is divided into three steps. At first, it decreased with increasing water content. And then, the reduction got slow, even the frozen depth increased slightly. At last, it decreased again; 4) The diameter of culvert greatly affected the maximum frozen depth, while the size of clear height had little effect on it. Diameter is priority for the requirements of strength, stability, functionality, etc. and the clear height could be considered at last.

seasonally frozen ground; culverts; temperature distribution; simulated calculation

10.11835/j.issn.1674-4764.2014.05.001

2014-05-06

國家自然科學基金(41071045)；鐵道部科技研究開發(fā)計劃項目(Z2012-062)；凍土工程國家重點實驗室自主立項課題(SKLFSE-ZY-15)

汪江紅(1988-)，男，主要從事熱力學、巖土工程方面研究，(E-mail)wjh-2008-ok@163.com。 牛永紅(通信作者)，男，博士，副研究員，(E-mail)yhniu@lzb.ac.cn。

TU445；U24

A

1674-4764(2014)05-0001-10