基于數(shù)值模擬的隧道溫度場分布規(guī)律研究

鐘小春

(中鐵五局集團(tuán)第一工程有限責(zé)任公司， 湖南長沙 410117)

季凍區(qū)隧道工程由于其所處的特殊環(huán)境，會受到外界環(huán)境低溫的影響，較容易發(fā)生凍害[1]，防寒抗凍問題一直是寒區(qū)隧道工程建設(shè)的重點(diǎn)問題。因此國內(nèi)外學(xué)者對季凍區(qū)隧道進(jìn)行了大量的研究，研究成果主要集中在幾個(gè)方面。

郝飛等[2]通過現(xiàn)場實(shí)測和數(shù)值模擬手段研究了寒區(qū)公路隧道的溫度場分布，并對保溫層的保溫效果進(jìn)行了評定；李鐵根[3]基于溫度場的控制方程和有限元公式，對實(shí)際隧道工程進(jìn)行理論計(jì)算和數(shù)值仿真，從而研究了氣溫對寒區(qū)隧道溫度場的影響；朱艷峰等[4]將有限體積法和有限單元法聯(lián)合使用，對圍巖溫度場進(jìn)行了有限元分析，探究了高速列車對溫度場的影響；孫克國等[5]依托實(shí)際寒區(qū)隧道工程，借助理論分析和數(shù)值計(jì)算手段得到了隧道橫縱斷面的分布規(guī)律；李思等[6]對寒區(qū)隧道溫度場分布規(guī)律進(jìn)行了研究，并探究了圍巖熱學(xué)參數(shù)對溫度場的影響；王仁遠(yuǎn)等[7]根據(jù)隧道現(xiàn)場實(shí)測數(shù)據(jù)，通過有限元軟件計(jì)算研究了不同因素對隧道溫度場的影響情況；李昊波等[8]根據(jù)實(shí)際隧道工程建立數(shù)值傳熱模型，研究隧道溫度場時(shí)空分布規(guī)律，并分析了風(fēng)速和風(fēng)溫的影響。

圍巖的熱學(xué)參數(shù)對溫度場有著較大影響，而在關(guān)于溫度場數(shù)值計(jì)算中，大多將圍巖簡化為單一巖性，計(jì)算結(jié)果可能與實(shí)際情況有較大差別。因此，本文依托青海關(guān)角山隧道，根據(jù)山體實(shí)際地形地貌，考慮圍巖分層情況，建立三維傳熱模型，研究隧道溫度場的時(shí)空分布規(guī)律。相關(guān)成果能夠?yàn)轭愃乒こ痰姆纼鲈O(shè)計(jì)起到一定指導(dǎo)作用。

1 工程概況

隧址區(qū)位于天峻縣關(guān)角山，屬構(gòu)造剝蝕中高山地貌。山嶺走向與構(gòu)造線走向大體一致，總體呈東北—西南向展布，隧道穿越處最高點(diǎn)標(biāo)高4 126.64 m，最低點(diǎn)標(biāo)高3 682.70 m，相對高差443.94 m。隧道頂板最大埋深410.53 m。

天峻縣屬大陸性高原氣候，具寒長暑短，四季不分明，無霜期短，日溫差大，多風(fēng)少雨，蒸發(fā)量大等特點(diǎn)。縣內(nèi)平均氣溫在0 ℃以下，為-1.1 ℃，各地區(qū)氣溫的分布，主要受地勢影響，海拔愈高氣溫愈低。極端最高氣溫28 ℃，極端最低氣溫-35.8 ℃，無霜期32天。天峻縣大風(fēng)日數(shù)較多，沙塵暴日數(shù)5.4天，全年平均大風(fēng)日數(shù)為70天，平均風(fēng)速為3.6 m/s，最大風(fēng)速為24 m/s。

縣區(qū)降水量季節(jié)性顯著，暖季降水充沛，各地5～9月的5個(gè)月中，集中了全年降水量的90%左右，區(qū)內(nèi)地質(zhì)災(zāi)害受降雨較明顯，全年以7～8月份是地質(zhì)災(zāi)害發(fā)生頻率最高，也是地質(zhì)災(zāi)害的易發(fā)時(shí)段。

2 數(shù)值模型及參數(shù)選取

本次建模采用ICEM-CFD建立網(wǎng)格模型(圖1)并用FLUENT軟件進(jìn)行數(shù)值模擬分析。為了保證數(shù)值模擬計(jì)算結(jié)果的準(zhǔn)確性，根據(jù)關(guān)角山隧道縱向地質(zhì)剖面圖和平面圖，考慮了實(shí)際的地形和巖石巖性來建立三維數(shù)值模型，模型長寬為3 600 m×70 m，高度的最大值為470 m，網(wǎng)格單元數(shù)量一共為776 266個(gè)，節(jié)點(diǎn)數(shù)一共為675 000個(gè)。

圖1 關(guān)角山隧道模型示意

本模型所考慮的巖性包括砂巖(強(qiáng)風(fēng)化、中風(fēng)化)、礫巖(強(qiáng)風(fēng)化、中風(fēng)化)、白云巖(強(qiáng)風(fēng)化、中風(fēng)化)、頁巖(強(qiáng)風(fēng)化、中風(fēng)化)，材料熱工參數(shù)如表1。流體為理想空氣，比熱容Ca=1006.43 J/(kg·k)，導(dǎo)熱系數(shù)λa=2.42×10-2W/(m·k)，黏度μ=1.789×10-5 kg/(m·s)，工作環(huán)境為標(biāo)準(zhǔn)大氣壓。洞內(nèi)風(fēng)流動狀態(tài)為紊流，水力直徑為11.242 m，湍流強(qiáng)度2.794。

依據(jù)隧址氣象站監(jiān)測數(shù)據(jù)進(jìn)行分析，對其日氣溫隨時(shí)間變化進(jìn)行正弦曲線擬合(以12月21日為起始點(diǎn)，見圖2)，氣溫?cái)M合函數(shù)為：

表1 材料熱工參數(shù)

T(d)=1-8.418×sin(2π(d+90)/365)

式中：T為洞口氣溫；d為天數(shù)。

圖2 日均溫度擬合曲線

由于建模時(shí)采用分塊組裝的方法，在考慮傳熱時(shí)，不同的組采用不同耦合面進(jìn)行熱的耦合，這里不同巖性間的耦合面采用自動生成的shadow耦合面，而巖石和初支采用手動設(shè)置的interface耦合面，排水管和巖體之間也采用手動生成的interface耦合面。

還原初始溫度場時(shí)，首先進(jìn)行穩(wěn)態(tài)計(jì)算，模型的上邊界取年平均溫度1 ℃(wall)；模型前、后、左、右邊界均為絕熱壁面(wall)；模型下邊界設(shè)定為熱流密度邊界(flux-wall)，其中距洞口500 m段地?zé)釤崃髅芏劝磓=0.04 W/m2取值，其他地段q=0.08 W/m2取值。然后再進(jìn)行瞬態(tài)計(jì)算，模型的上邊界取周期氣溫函數(shù)T(d)計(jì)算100年，模型前、后、左、右仍然設(shè)為絕熱壁面，下表面保持原來的設(shè)置不變。

計(jì)算隧道挖通后通風(fēng)的溫度場時(shí)，模型頂部仍然施加溫度為T(d)的溫度荷載，底部保持原來heat-flux設(shè)置不變。隧道進(jìn)口設(shè)置為velocity-inlet，其中速度為2 m/s，溫度取氣溫函數(shù)T(d)，隧道出口設(shè)置為pressure-outlet，計(jì)算的時(shí)間步長取為86 400 s。

3 隧道溫度場分布規(guī)律

通過計(jì)算，可得到隧道區(qū)域的初始溫度場分布云圖，具體如圖3所示。從中可以看出隧址區(qū)域的初始地溫場分布較為規(guī)律，在隧道埋深淺的區(qū)域地溫較低，埋深大的區(qū)域地溫較高。觀察溫度云圖的等溫線可以發(fā)現(xiàn)，隧址初始地溫的等溫線走勢與地形走勢基本一致，地溫場溫度梯度較為均勻。

圖3 初始溫度場溫度云圖

在初始溫度場基礎(chǔ)上進(jìn)行隧道挖通后的溫度場計(jì)算，這時(shí)考慮空氣流動的影響，風(fēng)速設(shè)為2 m/s，溫度按氣溫函數(shù)T(d)取，先進(jìn)行3年通風(fēng)計(jì)算來減少初始溫度場的影響，之后再進(jìn)行1年的通風(fēng)計(jì)算作為最終結(jié)果，以代表隧道運(yùn)營期的溫度分布狀態(tài)。

3.1 溫度場時(shí)間分布規(guī)律

提取距進(jìn)口100 m處橫斷面隧道拱底的不同埋深測點(diǎn)溫度數(shù)據(jù)，并采用正弦函數(shù)回歸法對不同測點(diǎn)溫度年變化進(jìn)行擬合，其一般形式：

式中:tm為年平均氣溫；Ay為年溫度振幅；d0為相位。

隧道橫斷面拱底各測點(diǎn)的溫度時(shí)間擬合曲線見圖4，其擬合的參數(shù)如表2所示。

圖4 拱底測點(diǎn)溫度擬合曲線

由圖4可知：

(1)距隧道進(jìn)口100 m處氣溫及隧道地溫溫度變化都符合正弦函數(shù)分布，但是存在明顯相位變化，氣溫和地溫間的相位差最大可達(dá)到78天，而不同埋深測點(diǎn)間的相位差為7～9天。

表2 拱底溫度擬合參數(shù)

(2)不同地溫測點(diǎn)的溫度隨時(shí)間的變化最低溫度都在50天左右，最高溫度都在220天左右。

(3)測點(diǎn)1最低溫度為-3.2 ℃，最高溫為5.4 ℃，正溫的時(shí)間段是從120天到320天，共200天，在剩下的165天為負(fù)溫，容易發(fā)生凍害應(yīng)該采取相應(yīng)工程措施。

3.2 隧道橫向溫度場分布規(guī)律

提取2月21日隧道進(jìn)口處、距隧道進(jìn)口500 m、1 000 m、1 500 m、2 000 m、2 500 m、3 000 m以及隧道出口處的隧道橫斷面上溫度分布云圖，具體如圖5所示。

圖5 沿隧道縱向不同橫斷面的溫度分布云圖

從圖5可以看出，縱向上靠近洞口的截面整體溫度相對于隧道中部的截面低，符合溫度分布兩邊低中間高的規(guī)律。同時(shí)可以看到徑向上洞身周圍及埋深較淺的區(qū)域的溫度低，埋深深的且離隧洞距離遠(yuǎn)的區(qū)域的溫度高，符合寒區(qū)隧道徑向上的溫度分布規(guī)律。

提取距隧道進(jìn)口100 m處不同時(shí)間的橫斷面溫度分布云圖見圖6。

圖6 距進(jìn)口100m不同時(shí)間橫斷面溫度云圖

由圖6可知：

(1)該斷面上不同時(shí)間的最低氣溫分布情況：在2月21日(-4.05 ℃)、10月21日(-1.92 ℃)、12月21日(-7.36 ℃)均為負(fù)溫，其中12月21日溫度最低；空氣溫度在4月21日(3.92 ℃)、6月21日(9.4 ℃)、8月21日(6.34 ℃)均為正溫，其中6月21日為最高溫，這和洞口外氣溫規(guī)律一致，且其隨時(shí)間呈明顯的三角函數(shù)變化規(guī)律。

(2)隧道洞身附近的溫度場受洞內(nèi)氣流影響導(dǎo)致一定范圍內(nèi)發(fā)生變異，在冬季出現(xiàn)凍結(jié)圈(其中凍結(jié)圈12月份為2.5 m，2月份為4.5 m)，徑向一定深度范圍內(nèi)出現(xiàn)變溫圈，其深度約為10 m，徑向深度超過10 m后，地層溫度基本不受洞內(nèi)氣流的影響，呈現(xiàn)恒溫規(guī)律。

3.3 隧道縱向溫度場分布規(guī)律

為進(jìn)一步研究隧道溫度場縱向上的時(shí)空分布規(guī)律，這里提取隧道底部深度分別為0.5 m、1.0 m、1.5 m、2.0 m、3.5 m的縱向測線，繪制溫度沿縱向的變化曲線，具體見圖7、圖8。

圖7 隧道底部溫度沿縱向距離變化曲線(埋深0.5m)

由圖7知：

(1)隧道縱向溫度場分布初始溫度場和通風(fēng)后溫度場分布規(guī)律差異較大，最大區(qū)別為溫度沿隧道縱向的溫度差，初始溫度場最低溫(1 ℃)發(fā)生在埋深最淺的進(jìn)口處溫度，最高溫度(14 ℃)發(fā)生在埋深最深處，其溫差為13 ℃，而通風(fēng)后隧道2月縱向溫度場最大溫差是3 ℃。

(2)隧道縱向溫度總體上看是滿足中間高兩邊低的規(guī)律，且溫度和地形起伏有明顯關(guān)聯(lián)。

(3)從2月21日到12月21日，隧道縱向上溫度先升高再降低，其中2月21日為溫度最低的月份，8月21日溫度為最高月份。這和氣溫的最低溫度月份12月21日最高溫度月份6月21日都具有2個(gè)月的滯后。

(4)此監(jiān)測線2 500 m之前最高溫度發(fā)生在6月21日，最低溫發(fā)生在12月21日，這和空氣溫度規(guī)律一致，但是2 500 m之后最高溫發(fā)生在8月21日，溫度最低發(fā)生在2月21日，這說明埋深為0.5 m時(shí)隧道縱向上溫度滯后規(guī)律不一致，縱向距離越大滯后現(xiàn)象越明顯。

圖8 隧道底部溫度沿縱向距離變化曲線

由圖7和圖8可知:

(1)埋深大于1 m后，溫度沿縱向上溫度分布規(guī)律基本一致且溫度隨時(shí)間變化規(guī)律也一致，同時(shí)溫度滯后的規(guī)律基本保持一致。

(2)隨著埋深從1 m到3 m的增加最低溫也從-3.5 ℃增加到-1.5 ℃，這說明埋深越大隧道溫度場溫度越高。

4 結(jié)論

本文通過建立三維數(shù)值傳熱模型研究季凍區(qū)隧道溫度場分布規(guī)律，依據(jù)數(shù)值計(jì)算結(jié)果，分析得到結(jié)論：

(1)隧址原始地溫的等溫線走勢與地形走勢基本一致，地溫場溫度梯度較為均勻。隧道開通后，隧道縱向溫度總體上呈現(xiàn)中間高兩邊低的規(guī)律，其分布與地形起伏有明顯關(guān)聯(lián)。

(2)與洞外氣溫變化規(guī)律一致，隧道內(nèi)氣溫和地溫隨時(shí)間發(fā)生正弦函數(shù)型變化，二者存在明顯的相位變化。

(3)隧道圍巖溫度隨徑向深度的增加而升高，當(dāng)埋深達(dá)到10 m后，地層溫度則基本不受洞內(nèi)氣流的影響，呈現(xiàn)恒溫規(guī)律。