地下近鄰平行精密恒溫隧道間傳熱的模擬研究*

何 玲 霍鏡濤 朱學錦 朱 喆 葉 蔚△ 張 旭

(1.同濟大學,上海;2.上海建筑設計研究院有限公司,上海)

0 引言

為減少地面溫濕度變化、太陽輻射、振動等環境因素對科學實驗的干擾,以北京正負電子對撞機、江門中微子實驗站等為代表的國家大科學裝置常設計于地下。借助地下相對穩定、無強熱擾的“室外”環境,這類特殊地下空間通過加裝恒溫恒濕空調以確保“室內”滿足精密恒溫的要求。與常規地下恒溫室不同,上海某在建大科學工程內含有3條近鄰平行精密恒溫隧道,隧道頂位于地下38 m。當3條隧道中的1條或1條以上隧道長期運行時,相鄰隧道互相成為恒溫空調的熱擾,影響程度需進一步研究。

本質上,相鄰恒溫空間會通過隧道圍護結構和周圍土壤進行熱量傳遞。Zhang等人發現隧道內部空間溫度除了受內部熱源影響外,圍護結構和土壤的傳熱也起著重要作用,熱量的傳遞與土壤及隧道圍護結構的特性有關[1]。胡政等人通過實際隧道巖體探測,結合地溫熱量傳遞理論,對隧道與周圍土壤的熱量傳遞進行了理論預測分析[2]。Song等人對土壤溫度場和隧道壁溫隨時間的演變進行了分析,求解熱平衡方程,提供了快速確定土壤蓄熱能力的方法[3]。此外,不同土壤的熱物性參數會影響恒溫空間的傳熱性能和溫度分布,很多學者也對此展開了研究[4-5]。

通常認為埋深大于12 m的地下建筑為深埋地下建筑。針對深埋地下建筑,學者們發現長通道式深埋地下建筑圍護結構沿長度方向上的傳熱過程幾乎可忽略[6];周期性內熱源作用下,巖壁1個月后傳熱進入穩定期[7]。就研究方式而言,Zeng等人通過縮尺實驗研究了對流傳導與隧道入口低溫氣流的耦合效應[8];Liu等人通過建立計算模型并結合現場實測,得出了隧道內含軸向坐標和時間的復雜邊界條件[9];王開運等人采用數值模擬研究了保溫層厚度對圍巖溫度的影響[10];王麗慧等人結合縮尺實驗和模擬,探究了土壤的蓄熱規律[11]。

對于地下空間和土壤傳熱的研究已有很多,但是針對地下多空間、多熱源之間相互影響的研究卻不多。本文針對某大科學實際工程存在的3條平行隧道,采用數值模擬方法,對地下隧道圍護結構的傳熱過程及不同運行工況下隧道圍護結構傳熱的影響進行研究。本研究可為地下空間多熱源的傳熱影響研究提供參考。

1 工程概況

本文研究對象為某大科學工程,含有3條隧道,隧道頂位于地下38 m,由文獻[12]可知,此處土壤處于恒溫層,年平均溫度約為18 ℃。隧道內左右兩側存在發熱電纜,熱量傳遞到隧道內的空氣中,且隧道內空氣溫度約為25 ℃,高于土壤溫度,因此隧道會通過圍護結構向周圍土壤散熱。對于該工程而言,3條隧道之間相距僅4 m,不但存在隧道與土壤的傳熱,還存在隧道之間的相互影響。

2 模型建立

2.1 模型簡介

為了研究隧道圍護結構的傳熱及3條隧道的相互影響,建立如圖1、2所示的2個模型。模型1僅含有1條隧道,模型2有3條隧道,將模型2中3條隧道分別命名為L、M、R。模型包括隧道、隧道外混凝土結構和土壤層三部分。根據文獻[12],土壤層為灰色砂質粉土和粉砂。

圖1 隧道模型1示意圖

圖2 隧道模型2示意圖

隧道的直徑為5.9 m,混凝土結構外徑為6.6 m,隧道之間最短距離為4 m,模型厚度為1 m。由于本文的研究對象為隧道圍護結構與土壤層的傳熱,故隧道內部簡化其他部分,僅有發熱電纜。電纜發熱簡化為面熱源,發熱量為35.25 W/m2;隧道內縱向風速約為0.42 m/s,空氣溫度約為25 ℃。

模型采用ICEM進行網格劃分,通過網格獨立性驗證后,最終選取的網格總數為1 169 564,土壤區域采用六面體結構化網格,隧道區域和混凝土結構采用非結構化網格,且對混凝土區域網格進行加密處理。計算求解采用分離式求解器,壓力速度耦合采用SIMPLE算法,采用二階迎風格式離散、RealizableK-ε模型,隧道內壁采用對流邊界。

2.2 模擬工況簡介

為比較不同隧道運行模式下圍護結構傳熱過程的差異,設置如表1所示工況。

表1 模擬工況

對于高精度溫控隧道而言,圍護結構單位面積傳熱量越小,隧道內溫度精度要求越容易實現。因此,分別對上述幾個工況進行模擬計算,比較不同工況下隧道外混凝土圍護結構與土壤的熱流密度。

3 結果分析

3.1 工況0(單隧道運行)

為了研究圍護結構傳熱過程和比較隧道間傳熱影響,單設1組空白對照組。不同運行時間隧道及土壤的溫度場分布如圖3所示。

圖3 工況0不同運行時間土壤溫度分布

由圖3可以看出:整體而言,隨著時間增長,隧道的發熱影響范圍變大;剛開始運行時,隧道與周圍土壤的溫差很大,隨著時間的增長,溫差和溫度梯度都減小;當僅有1條隧道時,土壤初始溫度一致且物理性質不變,各向同性,故土壤中熱量沿各個方向的擴散速率基本相等,等溫線為圓形。

表2給出了不同運行時間工況0隧道圍護結構的平均熱流密度。可以看出:在開始運行的第1個月內,平均熱流密度急劇減小,減小了92.45%;而后趨于平緩,維持一個較低的水平。

表2 工況0不同運行時間隧道圍護結構平均熱流密度

3.2 工況1(M隧道單獨運行10年)

圖4為M隧道不同運行時間隧道及土壤溫度分布云圖。

圖4 工況1不同運行時間隧道及土壤溫度分布

由圖4可以看出:隨著時間的增長,隧道的發熱影響范圍變大;同工況0,剛開始運行時,隧道與周圍土壤的溫度梯度很大,而時間越長,溫度梯度越小;兩側的空置隧道會對熱量擴散造成影響,在豎直方向上,傳熱幾乎相等,而水平方向上,受兩側隧道的影響,隧道間的溫度梯度更大;土壤等溫線表現為扁平的近橢圓形。

表3給出了不同運行時間M隧道圍護結構的平均熱流密度。可以看出,在M隧道開始運行的1個月內,熱流密度就降低了92.6%,從第1個月到第6個月,熱流密度減小了約1.75 W/m2,而后的每一年,圍護結構單位面積的傳熱量越來越低;在初始運行時間段,溫度梯度最大,平均熱流密度最大。

表3 工況1不同運行時間M隧道圍護結構平均熱流密度

同時,結合模型1空白對照組可發現:首先,土壤中僅1條隧道運行和M隧道單獨運行,圍護結構熱流密度變化規律一致,均在開始使用的第1個月內急劇減小,而后趨于平緩;其次,對比工況0和工況1的熱流密度數據,兩者變化基本一致,隨著運行時間的增長,兩者之間的偏差雖會增大,但是差值仍很小。

因此,對于研究結果作出推測,研究的3條隧道在地下土壤中平行存在,當只有1條隧道運行時,其圍護結構傳熱過程可簡化為單條隧道圍護結構的傳熱過程。

3.3 工況2(M隧道單獨使用10年后L、R隧道開始使用)

圖5為工況2不同運行時間隧道及土壤溫度分布云圖。

圖5 工況2不同運行時間隧道及土壤溫度分布

由圖5可以看出:M隧道使用10年后開始使用L、R隧道,相當于在原有溫度場的基礎上,左右兩邊各增加了新的穩定熱源,隨著時間的增長,溫度梯度逐漸減小,土壤中等溫線逐漸形成橢圓形;此外,新熱源增加也會影響原來隧道M的溫度分布和熱量擴散,M隧道水平方向溫度梯度減小,熱量擴散速率大大減小,而豎直方向上熱量擴散速率增大。

圖6顯示了開始使用兩側隧道后,M隧道周圍混凝土圍護結構熱流密度隨時間的變化,虛線對應的時間為10年后L和R開始運行的時間。可以看出:當開始使用左右兩側隧道后,M隧道的圍護結構單位面積傳熱量迅速減小,而后逐漸趨于平穩;兩側隧道的運行對M隧道產生了影響,使得M隧道由圍護結構向周圍土壤傳熱的熱流密度減小。

圖6 M隧道運行10年后圍護結構熱流密度的變化

表4給出了不同運行時間L、R隧道圍護結構的平均熱流密度。

表4 工況2不同運行時間L、R隧道圍護結構的平均熱流密度

2條新增運行隧道熱流密度隨時間的變化趨勢與M隧道單獨開始運行時一致,但數值上,L、R運行后的平均熱流密度遠小于M單獨開始運行時的熱流密度。由于兩側隧道開始運行時,中間的隧道已經運行了10年,周圍土壤的溫度已經升高,故兩側隧道圍護結構的熱流密度要小很多。

3.4 工況3(M、R隧道同時開始運行)

圖7為當M、R 2條距離較近的隧道同時運行時不同運行時間隧道及土壤溫度分布云圖。

圖7 工況3不同運行時間隧道及土壤溫度分布

當M、R隧道同時運行時,由于它們相距較近,圍護結構的傳熱會互相影響。與1條隧道單獨運行相比,2條隧道同時運行時,中間的土壤溫升很快。此外,由于左側空隧道的影響,熱量在兩側的傳播速率不同,土壤中等溫線為類橢圓形。

表5給出了不同運行時間M和R隧道圍護結構的平均熱流密度。

表5 M、R隧道圍護結構不同運行時間的平均熱流密度

同時開始運行2條隧道時,圍護結構傳熱變化規律與1條隧道基本一致。但數值上,與僅運行M隧道相比,運行2條隧道的圍護結構熱流密度衰減速率更快。由于2條隧道圍護結構傳熱過程的疊加效應和相互影響,2條隧道附近的土壤溫升比單獨1條隧道大,故隧道圍護結構的熱流密度也更小。

3.5 工況4(L、R隧道同時開始運行)

圖8為2條相距較遠的隧道L和R同時運行時不同運行時間隧道及土壤溫度分布云圖。

圖8 工況4不同運行時間隧道及土壤溫度分布

從圖8可以看出:由于L、R隧道相距較遠,相較工況3來說,L、R隧道的傳熱過程要經過較長時間才能互相影響;剛開始運行時類似僅1條隧道運行的情況,隨著時間的增長,等溫線分布近似于啞鈴狀。

表6給出了不同運行時間L、R隧道圍護結構的平均熱流密度。

表6 L、R隧道圍護結構不同運行時間的平均熱流密度

L、R隧道同時開始運行時,圍護結構傳熱變化趨勢與其他工況基本一致。與工況3相比,工況4中2條隧道的圍護結構熱流密度較大,原因是2條隧道相距較遠,隧道之間圍護結構傳熱影響較小。即2條隧道距離越遠,隨著時間的增長,隧道的平均熱流密度受彼此影響程度越小。

3.6 工況5(3條隧道同時開始運行)

圖9為3條隧道同時運行時不同時間隧道及土壤溫度分布云圖。

圖9 工況5不同運行時間隧道及土壤溫度分布

由圖9可以看出:當3條隧道同時開始運行時,土壤的溫升很快,傳熱過程也很快互相影響;隨著時間的增長,土壤中等溫線表現為類橢圓形。

表7給出了不同運行時間3條隧道圍護結構的平均熱流密度。

表7 3條隧道同時運行圍護結構不同運行時間的平均熱流密度

當3條隧道同時運行時,隧道圍護結構的傳熱變化趨勢也與其他工況相似。數值上,M隧道同時受到左右兩側隧道的影響,隨著時間的增長,M隧道的圍護結構熱流密度比兩側的小。

3.7 小結

各運行模式下熱流密度變化趨勢的比較如圖10所示。圖中紅色虛線為6個月運行時間的分界線。

圖10 各運行模式熱流密度隨時間的變化

由圖10可以看出:不同運行模式下,熱流密度在6個月后出現明顯差異;隧道運行數量越多、距離越近,相互影響就越大,此時熱流密度也越小,如M隧道在3條隧道同時運行時,受兩邊隧道的影響,熱流密度最小,而在單獨運行時熱流密度最大。

4 結論

本文采用數值模擬方法,對某大科學裝置的地下3條近鄰隧道運行時的隧道間傳熱進行了研究,比較了不同運行工況下隧道圍護結構熱流密度和周圍土壤溫度分布,得出以下結論:

1) 僅中間隧道使用時,圍護結構熱流密度在開始的1個月內下降速度最快,而隨著時間的增長逐漸變小。實際隧道運行過程中可在初始時間段對隧道進行通風處理,防止初始時間段的不穩定。

2) 兩側空隧道對中間隧道圍護結構傳熱的影響較小。當中間隧道運行后再開始運行另外2條隧道,因土壤溫度的升高,兩側隧道開始運行時的熱流密度僅為中間隧道不運行時的28%;且兩側隧道運行后,中間隧道圍護結構熱流密度也會衰減。

3) 當同時運行2條或3條隧道時,各隧道的圍護結構傳熱過程基本一致。隧道運行數量越多、距離越近,相互影響就越大,圍護結構熱流密度也越小;僅運行1條隧道的工況下,圍護結構熱流密度最大。

4) 各運行模式均在初始第1個月內熱流密度下降最快,可下降90%以上。不同運行模式下,熱流密度在6個月后出現明顯差異,最大值是最小值的115%;而10年后差異可達到268%。