許舒榮 XU Shu-rong

（北京中煤礦山工程有限公司，北京 100013）

0 引言

人工地層凍結(jié)是指在預(yù)先埋置在地層中的凍結(jié)管內(nèi)循環(huán)低溫冷媒，通過凍結(jié)管與土體進(jìn)行熱交換，從而形成完整性好、隔水效果效果佳、支護(hù)性能強(qiáng)的凍結(jié)帷幕。凍結(jié)法以其適用地層廣、環(huán)境污染小等特點(diǎn)，現(xiàn)已廣泛應(yīng)用于地鐵聯(lián)絡(luò)通道施工。根據(jù)凍結(jié)工程所處地層、凍結(jié)體量等因素，所采用的凍結(jié)設(shè)計(jì)參數(shù)有明顯差異[1-5]。為了保證安全的凍結(jié)過程以及良好的凍結(jié)效果，考慮土體初始溫度對(duì)凍結(jié)溫度場(chǎng)影響規(guī)律研究顯得意義重大。

黃建華等[6]建立了三維有限元模型，系統(tǒng)研究了水泥改良后的地鐵聯(lián)絡(luò)通道凍結(jié)溫度場(chǎng)發(fā)展特性，得出“水泥改良+人工凍結(jié)”對(duì)淤泥質(zhì)土地層具有顯著的加固效果；劉先觀[7]采用數(shù)值模擬與模型試驗(yàn)相結(jié)合的方法，考慮了凍結(jié)管直徑、間距、凍結(jié)管端部距管片的距離、等因素對(duì)聯(lián)絡(luò)通道凍結(jié)溫度場(chǎng)的影響規(guī)律；鄭立夫等[8]基于熱力耦合理論，通過有限差分軟件FLAC 3D建立了有限元模型，對(duì)凍結(jié)施工全過程進(jìn)行了模擬，研究了不同凍結(jié)壁厚度條件下地表凍脹、融沉變形及隧道管片變形規(guī)律，最終實(shí)現(xiàn)了凍結(jié)壁厚度的優(yōu)化設(shè)計(jì)。但目前國(guó)內(nèi)許多學(xué)者對(duì)凍結(jié)溫度場(chǎng)數(shù)值模擬研究多集中于聯(lián)絡(luò)通道，對(duì)地鐵出入口凍結(jié)溫度場(chǎng)發(fā)展規(guī)律的研究較少。本文利用有限元軟件，以上海地鐵某淺覆土地鐵出入口凍結(jié)加固工程為背景，建立了三維數(shù)值計(jì)算模型，對(duì)該出入口凍結(jié)溫度場(chǎng)發(fā)展與分布規(guī)律展開了數(shù)值模擬研究；并獲得了不同初始地溫條件下地鐵出入口凍結(jié)溫度場(chǎng)發(fā)展規(guī)律，為今后不同土層中類似工程的凍結(jié)設(shè)計(jì)與施工提供合理化建議。

1 工程概況

上海地鐵18號(hào)線某地鐵出入口凍結(jié)工程，該出入口周邊環(huán)境復(fù)雜，有酒店、寫字樓等構(gòu)筑物。出入口頂板覆土厚度約2.8m，且上方有電力、污水、信息等管線；出入口結(jié)構(gòu)呈“L”形，開挖斷面寬6.4m、高5.32m，暗挖段展開長(zhǎng)度為13.112m，開挖斷面大，若施工不當(dāng)易引起塌方冒頂?shù)群蠊?；根?jù)地勘報(bào)告顯示，出入口主要處于灰色淤泥質(zhì)粉質(zhì)粘土，該土層具有高含水量、高壓縮性、高靈敏度、低強(qiáng)度的特性。

根據(jù)工程特點(diǎn)，設(shè)計(jì)采用管棚凍結(jié)法加固地層后形成具有一定強(qiáng)度的凍結(jié)止水帷幕，再采用礦山法暗挖。共設(shè)計(jì)57個(gè)φ89×8mm的凍結(jié)管，19個(gè)φ273×10mm的管棚，同時(shí)在管棚中下放的φ89×8mm的凍結(jié)管；設(shè)計(jì)頂部?jī)鼋Y(jié)壁厚度為1m，側(cè)墻凍結(jié)壁厚度為1.9m，底部?jī)鼋Y(jié)壁厚度設(shè)置為2.5m。

圖1 凍結(jié)孔布置圖

2 數(shù)值模型建立

數(shù)值模型進(jìn)行以下簡(jiǎn)化，出入口所處地層較為復(fù)雜，為簡(jiǎn)化計(jì)算，將土層視為單一均質(zhì)的灰色淤泥質(zhì)粉質(zhì)粘土，溫度場(chǎng)計(jì)算過程中忽略地下水滲流與水分遷移的影響；出入口周圍的土體溫度變化是均勻的，凍結(jié)時(shí)潛熱集中在凍結(jié)界面連續(xù)放出。

根據(jù)勘察資料以及室內(nèi)試驗(yàn)結(jié)果，土體的計(jì)算參數(shù)如下，密度1860kg/m3，凍結(jié)溫度-0.6℃，凍土的導(dǎo)熱系數(shù)2.1w/（m·℃），未凍土的導(dǎo)熱系數(shù)1.4w/（m·℃），凍土比熱容為1kJ/（kg·℃），未凍土比熱容為1.3kJ/（kg·℃），相變潛熱取47.9kJ/kg。工程設(shè)計(jì)積極凍結(jié)時(shí)間為40d，鹽水降溫計(jì)劃為，7d內(nèi)鹽水溫度降至-18℃，15d內(nèi)鹽水溫度降到-24℃，開挖構(gòu)筑前鹽水溫度降至-28℃，凍結(jié)管外壁溫度荷載與鹽水降溫曲線一致。

依據(jù)實(shí)際工程圖紙中凍結(jié)孔設(shè)計(jì)文件，利用有限元軟件ANSYS建立含土體與凍結(jié)管的計(jì)算模型，采用熱力學(xué)計(jì)算模塊；出入口和凍結(jié)管依照實(shí)際工程1∶1建立，共建立72根φ89*8mm的凍結(jié)管，土體與凍結(jié)管均選用8節(jié)點(diǎn)的SOLID 70單元?？紤]到凍結(jié)法施工對(duì)周邊環(huán)境的影響，整體模型尺寸取35m（長(zhǎng)度）×25m（寬度）×13m（縱向）；靠近凍結(jié)管與設(shè)計(jì)凍結(jié)壁的土體的區(qū)域進(jìn)行網(wǎng)格加密，凍結(jié)管區(qū)域附近網(wǎng)格密度設(shè)置0.1m，以提高計(jì)算精度，遠(yuǎn)離凍結(jié)管與設(shè)計(jì)凍結(jié)壁的區(qū)域，減小網(wǎng)格密度，遠(yuǎn)離凍結(jié)管網(wǎng)格密度設(shè)置為0.3m，共計(jì)劃分214014個(gè)單元。

圖3 三維有限元模型圖

3 計(jì)算結(jié)果分析

3.1 凍結(jié)溫度場(chǎng)結(jié)果分析

按照建立的有限元模型進(jìn)行計(jì)算，選取10d、20d、30d、40d的溫度場(chǎng)云圖，如圖4所示。對(duì)計(jì)算結(jié)果分析如下。

圖4 凍結(jié)溫度場(chǎng)云圖

①凍結(jié)初期，凍結(jié)管與外界土體發(fā)生復(fù)雜的熱交換，凍結(jié)管周圍逐漸形成環(huán)狀凍土體，環(huán)狀凍土隨時(shí)間不斷向外擴(kuò)展，距離凍結(jié)管越近土體溫度下降越快；凍結(jié)第20d后，由于凍結(jié)管四周土體已降至冰點(diǎn)以下，未凍土轉(zhuǎn)變凍土?xí)r將釋放熱量，此時(shí)土體降溫速度減緩。

②距離凍結(jié)管位置近的土體降溫速率快，溫度較低；距離凍結(jié)管較遠(yuǎn)位置的土體，降溫速率慢，溫度較高。凍結(jié)至40d時(shí)，此時(shí)凍結(jié)管內(nèi)側(cè)的低溫區(qū)域面積明顯大于凍結(jié)管外側(cè)的低溫區(qū)域面積。

③若40d后繼續(xù)凍結(jié)將導(dǎo)致凍結(jié)孔內(nèi)側(cè)的凍土平均溫度不斷降低，強(qiáng)度隨之提高，對(duì)后續(xù)的通道開挖造成一定的困難，影響工程施工效率。

選擇右側(cè)墻凍結(jié)壁右邊界處測(cè)溫點(diǎn)T1，其降溫曲線如圖5可以看出。T1測(cè)點(diǎn)在凍結(jié)第40d時(shí)，其溫度由20℃下降到-3.26℃，已達(dá)到冰點(diǎn)，滿足設(shè)計(jì)要求。凍結(jié)第20d，T1測(cè)點(diǎn)溫度由20℃下降到4.76℃，凍結(jié)前20d平均降溫速率為0.76℃/d；凍結(jié)第30d，T1測(cè)點(diǎn)溫度由4.76℃下降到0.35℃，平均降溫速率為0.44℃/d；凍結(jié)第40d，T1測(cè)點(diǎn)溫度由0.35℃下降到-3.26，平均降溫速率為0.36℃/d。

圖5 T1測(cè)點(diǎn)降溫曲線

3.2 初始地溫對(duì)凍結(jié)溫度場(chǎng)的影響規(guī)律

為了解不同初始地溫條件下出入口凍結(jié)溫度場(chǎng)發(fā)展規(guī)律，選取測(cè)溫點(diǎn)T1，研究初始地溫15℃、20℃、25℃、30℃時(shí)測(cè)溫點(diǎn)T1的降溫規(guī)律，其降溫曲線如圖6所示?？梢钥闯觯煌跏嫉販貤l件下，測(cè)點(diǎn)降溫趨勢(shì)基本一致；凍結(jié)40d，初始地溫為15℃、20℃、25℃、30℃時(shí)，測(cè)溫點(diǎn)T1溫度分別為-4.17℃、-3.27℃、-1.90℃、-1.14℃，可以看出在本工程條件下初始地溫每升高5℃，T1測(cè)點(diǎn)的溫度約下降1℃；初始地溫為15℃、20℃、25℃、30℃時(shí)，測(cè)點(diǎn)降至冰點(diǎn)所需時(shí)間分別為30.0d、31.8d、33.8d、35.9d；可以看出在本工程條件下，初始地溫每升高5℃，T1測(cè)點(diǎn)降至冰點(diǎn)時(shí)間延長(zhǎng)約2d。

圖6 不同初始地溫條件下T1測(cè)點(diǎn)降溫曲線

4 結(jié)語

本文依托上海軌道交通18號(hào)線某出入口凍結(jié)工程，采用數(shù)值模擬分析了淺覆土地鐵出入口凍結(jié)溫度場(chǎng)發(fā)展規(guī)律，并研究了不同初始地溫對(duì)凍結(jié)溫度場(chǎng)的影響規(guī)律，主要得出以下結(jié)論：

①淺覆土地鐵出入口凍結(jié)溫度場(chǎng)發(fā)展規(guī)律表現(xiàn)為，凍結(jié)初期快速降溫速度較快，隨著凍結(jié)時(shí)間增長(zhǎng)，凍結(jié)降溫速率逐漸減緩；

②隨著初始溫度升高，達(dá)到設(shè)計(jì)凍結(jié)壁所需時(shí)間增長(zhǎng)；初始地溫每升高5℃，達(dá)到設(shè)計(jì)凍結(jié)壁所需時(shí)間延長(zhǎng)2d。