王麗慧 張 嫄 吳喜平 劉 俊

上海理工大學(xué)

隨著地鐵運行年限的增加，隧道區(qū)間溫度逐年上升的現(xiàn)象日益引起關(guān)注，部分區(qū)間隧道達到37℃[1]，溫度遠遠高出了設(shè)計值。區(qū)間溫升過高一方面將影響夏季列車車廂冷凝器散熱，導(dǎo)致車廂空調(diào)能耗上升；另一方面，區(qū)間的高溫空氣會隨活塞風(fēng)進入站臺（非屏蔽門系統(tǒng)），或者區(qū)間熱量也會通過門熱傳導(dǎo)或者在門開啟關(guān)閉過程中進入站臺（屏蔽門系統(tǒng)），這些都將引起站臺空調(diào)能耗的上升。區(qū)間隧道多深埋在距離地表15～20 m的位置，初始溫度一般在15℃左右，列車在行駛過程中放熱量是區(qū)間溫升的主要熱源，而區(qū)間半無限大土體的蓄放熱（隧道水平側(cè)和豎直向下側(cè)為半無限大土體）是區(qū)間溫升逐年變化的根本原因。本文則從機理上探討半無限大土體在周期性熱擾作用下的土體熱庫蓄放熱特性，為后續(xù)分析和解決區(qū)間隧道溫升過高的問題服務(wù)。

既有研究采用了現(xiàn)場實測，實驗臺縮尺模型模擬，理論分析和數(shù)值模擬等方法開展了相關(guān)研究工作。南京大學(xué)的李曉昭老師等[2]現(xiàn)場實測了隧道內(nèi)空氣溫度，壁面溫度，壁面熱流等變化規(guī)律，估算了圍巖傳熱量占總產(chǎn)熱量的比例最高可達26.6%；并以水代替空氣作為熱媒介質(zhì)采用縮尺模型實驗臺研究了穩(wěn)態(tài)土體溫度熱量傳遞。肖琳等[3]利用模型實驗研究了地鐵圍巖內(nèi)土體熱導(dǎo)率值；而絕大部分學(xué)者選擇了理論分析和數(shù)值模擬，王海彥和胡增輝等[4,5]通過對隧道內(nèi)有限范圍的圍巖體溫度場的理論分析，得出了隧道土體內(nèi)各個位置處的溫度變化規(guī)律。同濟大學(xué)于連廣等[6,7]提出了考慮氣固耦合傳熱、地下水滲流等因素的隧道土體溫度三維預(yù)測模型，指出土體蓄熱作用在地鐵環(huán)控系統(tǒng)的設(shè)計中應(yīng)予考慮；既有相關(guān)研究中，多以水代替空氣模擬區(qū)間氣體與隧道圍巖土體之間的熱傳遞，而忽略了熱濕耦合對傳熱的影響；隧道圍巖大多以穩(wěn)態(tài)傳熱為主，未考慮非穩(wěn)態(tài)傳熱實際工況。本文運用軟件模擬分析上海地區(qū)地鐵隧道圍巖溫度場和含水率變化規(guī)律，為以后地鐵環(huán)控系統(tǒng)的設(shè)計提供參考。

圖1 土壤熱濕傳遞實驗臺原理圖

1 縮尺模型實驗

本實驗主要是研究地鐵區(qū)間隧道內(nèi)活塞風(fēng)的溫度對隧道土壤溫度場的影響，因此采用空氣作為熱量傳遞介質(zhì)，這樣既與實際模型相吻合又彌補了以水作為熱媒介質(zhì)不能體現(xiàn)濕傳遞的缺陷；變工況的空氣溫度以上海地鐵一號線衡山路站全年溫度的實測值為依據(jù)，取全年每月的平均溫度作為實驗臺上空氣的溫度，這樣實驗臺上的空氣溫度既方便調(diào)節(jié)，又能夠較好的接近隧道內(nèi)的實測溫度。依照上述需求，設(shè)計和搭建了一套用于研究隧道圍巖土壤溫度場率變化規(guī)律的實驗裝置如下。

1.1 實驗臺組成

本次實驗設(shè)計的土壤熱濕傳遞實驗臺原理圖如圖1所示。

本實驗裝置主要由三部分組成，分別為：空氣加熱及送風(fēng)裝置（A）、土體及溫濕度和熱流密度采集裝置（B）、實驗臺支撐框架（C）?？諝饧訜峒八惋L(fēng)裝置（A）主要由變頻風(fēng)機、空氣加熱器、控制柜、風(fēng)管、送風(fēng)參數(shù)測試孔構(gòu)成，；土體及溫濕度采集裝置（B）主要由土壤表面熱流密度探頭、土壤濕度探頭、土壤溫度探頭、土體、數(shù)據(jù)采集板構(gòu)成，熱流密度板布置在空氣與土壤接觸的表面；實驗臺支撐框架（C）主要由風(fēng)管支撐架、加熱器支撐架、可移動整體支撐架構(gòu)成，所有的裝置放置在可移動整體支撐架上面，且可拆卸移動。實驗過程中土壤熱濕傳遞實驗臺置于恒溫恒濕室，空氣從恒溫恒濕室送風(fēng)系統(tǒng)進入，依次經(jīng)過變頻風(fēng)機、空氣加熱器、風(fēng)管，風(fēng)管出口與恒溫恒濕室排風(fēng)口相連；利用變頻風(fēng)機調(diào)節(jié)風(fēng)速，通過空氣加熱器調(diào)節(jié)空氣溫度，恒溫恒濕室空調(diào)箱調(diào)節(jié)溫度的降低，同時也可以調(diào)節(jié)空氣的濕度，空氣經(jīng)過風(fēng)管時掠過土壤表面，反應(yīng)實際隧道模型土壤熱濕傳遞過程。

2 CHAMPS軟件介紹

2.1 CHAMPS-BES概況

CHAMPS-BES（Coupled Heat、Air、Moisture and Pollutant Simulation in Building Envelope Systems）由德國德累斯頓工業(yè)大學(xué)和美國雪城大學(xué)等合作研發(fā)，獲得美國環(huán)保局和美國能源部的科研資金支持，專門用于多孔介質(zhì)熱濕傳遞模擬的軟件，可對一維、二維及旋轉(zhuǎn)對稱的三維模型進行模擬。針對地鐵隧道內(nèi)空氣和土壤熱濕耦合傳遞模型，CHAMPS-BES可以出色的完成幾何模型的幾何建模、網(wǎng)格離散劃分；可以添加包括模型介質(zhì)溫度、相對濕度、材料參數(shù)，邊界空氣溫度、風(fēng)速、對流換熱系數(shù)等條件，支持逐時標準年氣象數(shù)據(jù)鏈接，使模擬條件盡可能與實際情況相符合。

2.2 模型的建立及定解條件的確定

給出隧道內(nèi)全年空氣溫度和空氣流速、土壤初始溫度和含水率、隧道壁面與空氣的對流換熱系數(shù)對空氣和土壤的熱濕傳遞過程求解采用第三類邊界條件。模型建立時做如下假設(shè)：（1）不考慮圍巖土壤內(nèi)地下水滲流的影響；（2）假設(shè)隧道的埋深足夠深，處于地下恒溫層內(nèi)，不考慮空氣溫度的年變化對土壤溫度的影響；（3）假設(shè)隧道壁面溫度分布在X方向上是均勻的；（4）假設(shè)隧道壁面和土壤接觸良好，對土壤的傳熱特性沒有影響。

邊界條件的取值如表1所示。

表1 各項邊界條件參數(shù)

圖2 隧道圍巖結(jié)構(gòu)

圖3 實驗和軟件模擬在各年的溫度對比曲線

2.3 地鐵隧道圍巖物理模型

地鐵隧道圍巖物理結(jié)構(gòu)如圖2所示，分析可知地鐵傳熱計算模型采用多層平壁計算模型，溫度的分布為深度和時間的函數(shù)。結(jié)構(gòu)由空氣、隧道混凝土圍巖和土壤構(gòu)成，空氣-土壤耦合傳熱過程簡化為沿土壤深度方向（Y方向）的傳熱；參照地鐵站與站之間的隧道長度大概為1 000 m，此物理結(jié)構(gòu)中X方向的長度取為1 000 m，為了使傳熱充分且不受邊界的影響，Y方向深度取為50 m，同時在空氣與土壤之間有一層厚度為20 cm的混凝土結(jié)構(gòu)；空氣溫度為變工況條件，作用于圍巖1表面，與土壤進行熱濕耦合的傳熱、傳濕。

2.4 實驗結(jié)果與軟件模擬的對比驗證

為了驗證實驗結(jié)果的準確性，同時為了說明軟件模擬準確可靠，能夠用于后續(xù)的研究中，采用軟件模擬和實驗結(jié)果相對比的方式互相驗證。

表2 土壤表面實測溫度與設(shè)定溫度偏離程度

軟件模擬建立了實驗臺土體的幾何模型，模擬了土壤從初始恒溫14.3℃到17個作用周期以后溫度場的分布情況，將實驗結(jié)果和軟件模擬結(jié)果的對比如下：

從圖3可見，實驗結(jié)果和軟件模擬值存在一定的差異，主要表現(xiàn)在實驗溫度值在熱庫峰值之前比軟件模擬值高，而在熱庫峰值之后又比模擬值低，但是整體分布曲線吻合較好，隨土壤深度的增加，溫度場分布規(guī)律一致，熱庫峰值大小和熱庫峰值位置都相差不大，實測溫度值和軟件模擬溫度值之差保持在0.8℃之內(nèi)。由上述對比分析，一方面說明了實驗結(jié)果準確可信，另一方面驗證了軟件分析結(jié)果和實際情況有很好的吻合，證明了用軟件進行后續(xù)研究的可行性。

3 上海隧道圍巖土壤溫度和含水率演化規(guī)律研究

空氣側(cè)溫度以上海地鐵衡山路站全年空氣溫度的實測值作用于土壤表面，土壤的土質(zhì)參數(shù)主要包括密度、導(dǎo)熱系數(shù)、比熱容、容積含水率。土壤參數(shù)的測試實驗及上海地區(qū)的土質(zhì)為黏性黑土，取得上海地區(qū)地鐵隧道圍巖土壤的參數(shù)如表3所示。

表3 上海地區(qū)土壤參數(shù)

按照上述氣象參數(shù)及土壤參數(shù)利用CHAMPSBES軟件模擬上海地區(qū)地鐵隧道圍巖從新建到運行20年以后土壤溫度和含水率變化情況，結(jié)果分析如下。

3.1 地鐵隧道圍巖溫度場變化特性

軟件模擬是一種理想化的計算，只要計算一開始模型內(nèi)各個位置處的溫度都會發(fā)生變化，只是變化的大小隨深度會有不同，離壁面越近則變化幅度越大，而實際地鐵隧道模型傳熱要比軟件模擬復(fù)雜，會受到諸多因素的干擾，溫度場的變化相比軟件模擬不會那樣靈敏，從嚴格意義上講，軟件模擬中溫度恒定的位置是不存在的，因此將軟件模擬中溫度變化小于0.5℃忽略不計，只有在溫度變化大于0.5℃時才認為溫度場受到了波動。

由圖4可見，通過軟件模擬的上海地鐵圍巖土壤溫度場隨年限的變化規(guī)律和實驗結(jié)果是相似的。溫度場的變化主要發(fā)生在1～10年，在此期間，隨年限的增加，熱庫峰值不斷升高，峰值位置和熱庫厚度不斷加深；第10年以后溫度場趨于穩(wěn)定，圍巖土壤和隧道空氣達到吸放熱的動態(tài)平衡，軟件模擬每年末的熱庫曲線幾乎保持不變，第10年以后在0.0～0.6 m范圍內(nèi)溫度基本不變，而在0.6～1.0 m的范圍內(nèi)溫度繼續(xù)升高。

由圖4可見，整體溫度場在初始時刻溫度為17℃，傳熱進行1年以后形成熱庫曲線，該曲線為類似于隧道內(nèi)空氣溫度變化曲線的拋物線，對應(yīng)存在熱庫曲線的熱庫峰值和熱庫厚度，每1年末熱庫峰值逐漸增高，熱庫厚度逐漸擴大，最后逐漸趨向穩(wěn)定；1～5年溫度場分布變化明顯，熱庫峰值由第1年末的23.5℃升高到第5年末的25.2℃，波峰位置由距壁面1.4 m移動到2 m，熱庫厚度則由6m加深到16 m左右，熱庫峰值逐漸變大，峰值位置沿土壤深度方向移動，熱庫厚度加深；6～10年熱庫峰值有較小幅度的增高，熱庫范圍由第6年末的16 m左右變化到第10年末的22 m左右，波峰位置維持在2 m左右基本不變，溫度場的變化主要集中在5～22 m的范圍之內(nèi)，而在5 m之前的位置溫度場變化已經(jīng)趨于穩(wěn)定；11～15年期間，熱庫峰值基本穩(wěn)定在25.7℃左右，熱庫位置在距壁面距離2.2 m，熱庫范圍在25 m左右，6～22 m的范圍之內(nèi)有較小的溫度升高；16～20年期間，熱庫峰值、波峰位置、熱庫范圍的變化都較小，變化基本可以忽略，認為保持不變。

隧道圍巖土壤30 m范圍內(nèi)各個位置處1～10年、11～20年溫度變化情況如表4和圖5所示。由圖可見，隧道壁面到18m距離范圍之內(nèi)土壤溫度場在1～10年之內(nèi)的變化要大于11～20年之內(nèi)的變化；1～10年的變化隨深度方向呈拋物線的形式，1 m處的溫升為1.51℃，沿深度方向逐漸增大，到6 m位置處達到最大值4.61℃，之后逐漸降低到1.22℃；11～20年的變化隨深度方向呈對數(shù)曲線的形式，距壁面越遠溫升越大，說明在此期間離壁面越遠傳熱波動越強烈，由0.17℃逐漸升高到1.22℃。18～30 m范圍內(nèi)土壤溫度場在11～20年之內(nèi)的溫升幅度要大于1～10年之內(nèi)的溫升幅度；主要是因為土壤傳熱存在時間上的延遲，11年以后土壤熱量傳遞到了土壤深處，同時1～10、11～20年18 m之后沿深度方向溫升幅度逐漸減小。

由此可見，上海地區(qū)新建地鐵隧道圍巖土壤溫度場的變化主要集中在從新建到投入運營前10年之內(nèi)，在此期間隧道圍巖與空氣進行熱量交換，吸熱大于放熱，土壤內(nèi)溫度場逐年升高；運營10年以后近處土壤達到吸放熱平衡，遠處土壤溫度仍逐年有少許上升，可以認為整體溫度場分布趨于穩(wěn)定，隧道圍巖土壤和空氣之間達到了吸放熱的平衡，土壤吸熱能力下降。

圖4 整體溫度場每年分布曲線

表4 各個位置處1～10年、11～20年溫度變化值

圖5 各個位置處1～10年、11～20年溫度變化

圖6 土壤容積含水率分布曲線

圖7 不同位置處容積含水率變化曲線

圖8 20年內(nèi)不同位置處容積含水率的變化曲線

3.2 土壤容積含水率的變化過程

隧道圍巖土壤容積含水率的分布曲線如圖6所示，土壤的初始容積含水率為0.426 m3/m3，容積含水率的變化主要發(fā)生在3.5 m的范圍之內(nèi)，隧道壁面到0.2 m范圍內(nèi)容積含水率變化較大，且第1年末到第20年末基本保持不變；0.2 m～3.0 m范圍內(nèi)為容積含水率主要發(fā)生變化的區(qū)域，隨著年限的推移，容積含水率逐年減少，但是每一年變化的幅度越來越小，容積含水率的分布逐漸趨向于穩(wěn)定；超過3.0 m范圍的土壤容積含水率基本保持不變。

隧道圍巖不同位置處容積含水率隨年限的變化過程如圖7所示，隨著離壁面距離的增加，容積含水率的變化越來越??；0.5 m處容積含水率的變化最大，第1年到第20年減小了0.052 m3/m3，4.0 m處容積含水率保持不變。隨著距壁面距離的增加，容積含水率的變化同樣存在時間上的延遲，離壁面距離越近，變化越快，例如0.5 m和1.0 m處的容積含水率呈現(xiàn)指數(shù)分布的形式，開始幾年變化大，之后變化逐年減少減慢；離壁面越遠，開始變化得越慢，如1.5 m處從第3年開始變化，2.5 m處從第8年左右開始變化，4.0 m處幾乎沒有變化。

20年內(nèi)不同位置處容積含水率變化的分布情況如圖8所示，由此可見容積含水率的變化區(qū)域大約在土壤深處3 m左右。

3.3 土體熱庫蓄放熱量演化特性

土體熱庫演化本質(zhì)在于土體蓄放熱能力。對于實驗臺中土體蓄放熱量大小，一方面可根據(jù)土體表面熱流密度板采集的熱流密度得到，另一方面可通過土體溫度沿深度方向的分布曲線積分而得。兩種方法相互驗證，同時也證明了該實驗臺的熱平衡。

表5給出了實驗臺土體表面各周期凈熱流密度和，其考慮了土體吸放熱過程熱流密度的代數(shù)和（正值表示熱量從空氣傳入土體，負值表示熱量由土體傳入空氣）。從實驗結(jié)果可見，土體逐年均從空氣側(cè)吸熱，但吸熱量逐年減少；同時可知新建地鐵區(qū)間隧道的年蓄熱量為遠期年蓄熱量的11.6倍。

另一方面，用MATLAB將圖5中的土體溫度擬合為深度的函數(shù)，然后根據(jù)各深度溫度差值積分求得土體蓄熱量的年變化量，并將蓄熱量折算到原型隧道（長1 000 m，直徑為6 m）。最后將兩種方法得到的土體逐年蓄熱量進行對比。

從圖9可見，兩種方法得到的區(qū)間隧道土體蓄熱量隨年限的變化規(guī)律基本一致，尤其在1～4年及10～17年吻合較好?！盁釒靹討B(tài)拓展期”內(nèi)土體每年的凈吸熱量逐年減少，在未考慮隧道盾構(gòu)和未控制土體含水率的情況下，對于一個直徑6 m，長度1 000m的地鐵隧道，第1年的圍巖吸熱量可達80 000 MJ。雖然土體的蓄熱量在逐年減小，但前10年總蓄熱量可達300 000 MJ；在10～17年的“熱庫動態(tài)穩(wěn)定期”內(nèi)，吸熱量曲線變化較小，保持在5 000 MJ（265.26 kJ/m2·年）左右，說明各年的吸熱量稍微大于放熱量，土體溫度有較小的變化，因此每個作用周期末的土體熱庫曲線基本相同，熱庫厚度和熱庫峰值變化較小，土體的蓄熱量與吸熱量基本實現(xiàn)動態(tài)平衡。

表5 實驗對應(yīng)實際模型每年的換熱量

4 結(jié)論

本文主要采用CHAMPS-BES軟件對實際隧道圍巖土壤傳熱過程，土壤不同參數(shù)對圍巖土壤的傳熱影響進行了分析，確定了軟件模擬的各個邊界條件；其次，通過軟件模擬了實驗臺的傳熱過程，對比軟件模擬和實驗結(jié)果發(fā)現(xiàn)兩者偏差較小，證明了軟件模擬準確可靠；最后以上海衡山路站隧道內(nèi)全年空氣溫度為依據(jù)和上海地區(qū)土質(zhì)參數(shù)為邊界條件，模擬了實際地鐵隧道圍巖土壤的熱濕耦合過程，分析了熱庫的形成、不同位置處溫度的波動、土壤容積含水率的變化過程。利用實驗采集的熱流密度與溫度分布曲線分別計算出一個直徑6 m，長度1 000 m的地鐵隧道圍巖總的蓄熱量。在“熱庫動態(tài)拓展期”，第1年的吸熱量可達80 000 MJ，前10年總的年蓄熱量可達300000MJ?？梢婋S著運行年限地增加，隧道圍巖吸熱能力逐年下降，導(dǎo)致區(qū)間熱量不斷累積，引起區(qū)間溫升，增大車廂站臺能耗。因此在今后的研究工作中，宜考慮采用合理的技術(shù)措施以可持續(xù)的利用區(qū)間盾構(gòu)蓄存熱量的能力，減緩區(qū)間溫度的上升。

圖9 17年每年土體的總蓄熱量對比圖（對應(yīng)直徑6 m，長1 000 m的地鐵隧道）

[1] 龔蓂杰. 地鐵區(qū)間隧道溫度特性及其對站臺環(huán)境的影響[D].重慶：重慶大學(xué), 2014.

[2] 李曉昭，熊志勇，喬恒君，馬娟，張學(xué)華，杜茂金. 地鐵圍巖傳熱規(guī)律的監(jiān)測分析[J].地下空間與工程學(xué)報，2012，8（1）：105-110.

[3] 肖琳，楊成奎，胡增輝，李曉昭，李茉. 地鐵隧道圍巖內(nèi)溫度分布規(guī)律的模型實驗及其熱導(dǎo)率反算研究[J].巖土力學(xué)，2010，31（2）：86-91.

[4] 王海彥，駱憲龍，楊石柱. 隧道圍巖溫度場變化規(guī)律理論分析[J].石家莊鐵路職業(yè)技術(shù)學(xué)院學(xué)報，2006，5（2）：25-29.

[5] 胡增輝，李曉昭，趙曉豹，肖琳，武偉. 隧道圍巖溫度場分布的數(shù)值分析及預(yù)測[J].地下空間與工程學(xué)報，2009，5（5）：867-872.

[6] 于連廣，吳喜平，駱澤彬. 圓形腔體內(nèi)熱桿件軸向同心移動湍流換熱模型[J].同濟大學(xué)學(xué)報：自然科學(xué)版，2010，10：1501-1505.

[7] 于連廣，吳喜平，李昊翔.三維地鐵隧道土體溫度預(yù)測模型[J].沈陽工業(yè)大學(xué)學(xué)報，2011，33（2）：234-240.