3排管凍結(jié)梯形-拋物弓疊合等效溫度場模型和平均溫度

胡向東，任 輝

(1.同濟(jì)大學(xué) 地下建筑與工程系，上海 200092；2.同濟(jì)大學(xué) 巖土及地下工程教育部重點(diǎn)實(shí)驗(yàn)室，上海 200092)

3排管凍結(jié)梯形-拋物弓疊合等效溫度場模型和平均溫度

胡向東1,2，任 輝1,2

(1.同濟(jì)大學(xué) 地下建筑與工程系，上海 200092；2.同濟(jì)大學(xué) 巖土及地下工程教育部重點(diǎn)實(shí)驗(yàn)室，上海 200092)

由于3排管凍結(jié)溫度場的復(fù)雜性，不便于對凍土帷幕溫度場進(jìn)行簡便描述。故以現(xiàn)有的3排管凍結(jié)溫度場的勢函數(shù)疊加法解析解為理論基礎(chǔ)，根據(jù)溫度場的周期性，選擇主面與界面之間的合適位置作為等效截面，以等效截面的溫度場等效整體溫度場。以平均溫度為檢驗(yàn)標(biāo)準(zhǔn)，簡化等效截面的溫度場曲線，得出方便應(yīng)用的梯形拋物弓疊合等效溫度場模型?；诘刃囟葓瞿Ｐ痛_定平均溫度計(jì)算方法，得出免于積分運(yùn)算的平均溫度簡便計(jì)算公式。該公式計(jì)算精度滿足工程要求，表明等效溫度場及其平均溫度計(jì)算方法具有良好的實(shí)用性。

3排管凍結(jié)；等效溫度場；平均溫度；凍土帷幕；人工地層凍結(jié)法

凍結(jié)溫度場的研究是人工地層凍結(jié)法理論研究的基礎(chǔ)，也是凍結(jié)設(shè)計(jì)施工的理論依據(jù)。隨著人工地層凍結(jié)法在礦山建設(shè)、城市地下工程和河底隧道工程中的廣泛應(yīng)用，單排凍結(jié)管溫度場以及雙排凍結(jié)管溫度場的研究取得了一系列成果[1-3]，筆者也進(jìn)行了單、雙排管凍結(jié)溫度場及其平均溫度的完善和應(yīng)用性研究[4-10]。

由于工程實(shí)踐的需要，為了實(shí)現(xiàn)厚度較大和平均溫度較低的凍土帷幕，多排(3排及3排以上)管凍結(jié)方案被廣泛地應(yīng)用于大水流地層礦山井筒和深井凍結(jié)、盾構(gòu)進(jìn)出洞地基加固等工程中。但是對多排管凍土帷幕溫度場及其平均溫度計(jì)算方法的研究極少，有學(xué)者通過數(shù)值模擬獲得了經(jīng)驗(yàn)公式[11]，筆者基于巴霍爾金解析解的雙排管凍結(jié)溫度場公式通過作圖的方法應(yīng)用于描述多排管凍土帷幕溫度場，并得到多排管凍土帷幕厚度及平均溫度的計(jì)算公式[12]。也有學(xué)者在上海地鐵某一區(qū)間隧道修復(fù)采用四排局部垂直凍結(jié)時(shí)，對凍結(jié)壁溫度場發(fā)展進(jìn)行實(shí)時(shí)監(jiān)測，分析多排管局部凍結(jié)凍土壁溫度場特性[13]，但沒有對多排管凍土壁進(jìn)行深入的理論研究。工程界也有采用計(jì)算單排管凍土帷幕平均溫度的“成冰”公式[14]近似計(jì)算多排管凍土帷幕平均溫度的情況。故筆者針對3排管凍結(jié)溫度場進(jìn)行研究，得出3排管凍結(jié)溫度場的勢函數(shù)疊加法解析解，并驗(yàn)證了其解析解結(jié)果的可靠性[15]?？紤]3排管凍結(jié)溫度場的復(fù)雜性，為方便研究凍土帷幕力學(xué)性能和指導(dǎo)凍結(jié)設(shè)計(jì)施工，筆者將以3排管凍結(jié)溫度場解析解為理論基礎(chǔ)，探求3排管凍土帷幕的在二維平面的等效溫度場模型。

1 3排管凍結(jié)溫度場

3排管凍土帷幕凍結(jié)管布置方式一般可分為兩大類：一為3排凍結(jié)管直線型對齊排布，二為相鄰兩排凍結(jié)管管間距錯(cuò)開l/2排列，成為3排凍結(jié)管直線型交錯(cuò)排布，如圖1所示。

根據(jù)3排管凍結(jié)溫度場的勢函數(shù)疊加法解析解理論[15]，以圖1所示坐標(biāo)系(xoy)為準(zhǔn)，3排管凍土帷幕穩(wěn)態(tài)溫度場分布的解析解表達(dá)式為

(1)

式中，t(x,y)為任意點(diǎn)(x,y)的溫度，℃；tct為凍結(jié)管表面的溫度，℃；l為每排兩根凍結(jié)管的間距，m；L為3排凍結(jié)管的排距，m；ξ為每排凍結(jié)管至凍土墻邊緣的凍土厚度，m；r0為凍結(jié)管外半徑，m；ω為區(qū)別3排管兩種排列方式的符號：對于3排凍結(jié)管對齊排列，ω=1；對于3排凍結(jié)管交錯(cuò)排列，ω=-1。

圖1 3排管凍土帷幕凍結(jié)管布置方式Fig.1 Arrangement of three-row-pipe freezing

2 3排管等效溫度場

2.1 等效溫度場的確定

不論從式(1)還是圖2都可以直觀地看出3排管凍土帷幕的溫度場的復(fù)雜，不利于理論研究和工程應(yīng)用，為此我們希望找出一個(gè)等效截面x=xe，并以此截面上的溫度t(xe,y)來替代真實(shí)溫度場。

圖2 3排管凍土帷幕溫度場示意Fig.2 Temperature field distribution of three-row-pipe freezing

考慮凍結(jié)管布置的對稱性，及溫度場在x軸上的周期性，整個(gè)凍土帷幕只需要研究其中0≤x

分析式(1)可知，當(dāng)凍結(jié)管外半徑和凍結(jié)管分布形式確定時(shí)，整個(gè)凍土帷幕的溫度場主要受凍結(jié)管間距l(xiāng)、凍結(jié)管排距L、凍結(jié)管至凍土墻邊緣的凍土厚度ξ的影響，針對上述3個(gè)參數(shù)的不同取值，利用式(1)計(jì)算結(jié)果考察0≤x

選取x=0,l/6,l/4,l/3,l/2這5個(gè)特征截面，觀察其凍土帷幕溫度曲線t(x,y)的分布情況。根據(jù)凍結(jié)管兩種不同排布方式，選擇不同參數(shù)取值得到352組凍土帷幕的5個(gè)特征截面溫度分布情況，參數(shù)具體取值詳見表1。觀察每組數(shù)據(jù)的溫度曲線可知，在相同的凍結(jié)管排布方式下，曲線的形狀都是類似的。故針對凍結(jié)管兩類不同排布方式，選取1組參數(shù)(L=1.0m，l/L=1.0，ξ/L=0.5)為例進(jìn)行具體說明，描述這5個(gè)特征截面的溫度分布情況，如圖3所示。

表1參數(shù)取值
Table1Parametervalues

參數(shù)取值L1 0，1 2，1 4，1 6ml/L0 7，0 8，0 9，1 0ξ/L0 5，0 6，0 7，0 8，0 9，1 0，1 1，1 2，1 3，1 4，1 5

圖3 特征截面凍土帷幕溫度場Fig.3 Temperature field distribution of the characteristic cross-sections

根據(jù)圖3顯示可知，3排管在主面(x=0)與界面(x=l/2)之間不同截面的溫度場主要是在凍結(jié)管位置溫度變化較大，尤其是3排管外面兩排凍結(jié)管位置溫度變化更明顯，并且在主面(x=0)到截面(x=l/6)溫度變化梯度最大，隨后逐漸減小，這就表明在凍結(jié)管位置從主面(x=0)到界面(x=l/2)溫度是迅速升高，使得截面溫度曲線趨于一致。由此可知，僅主面附近小部分范圍因靠近凍結(jié)管溫度較低，使得此部分截面溫度曲線在凍結(jié)管位置有較大變化，而其他大部分截面溫度曲線型狀基本類似。因此，考慮用其中某一特征截面的溫度曲線來近似描述整個(gè)溫度場的溫度變化情況，此特征截面即為等效截面，筆者擬采用主面(x=0)與界面(x=l/2)之間的中間位置(xe=l/4)作為等效截面。

根據(jù)上述分析結(jié)果，xe=l/4特征截面作為等效截面來描述整個(gè)凍土帷幕溫度場，但xe=l/4代入式(1)后所得溫度場公式仍然很復(fù)雜，故對特征截面(xe=l/4)進(jìn)行處理，觀察xe=l/4的溫度曲線，可將溫度曲線簡化為如圖4所示的等效溫度曲線(在參數(shù)相同的情況下，xe=l/4截面的對齊排布和交錯(cuò)排布溫度曲線一致)，曲線分上下兩部分，上部分為拋物弓型狀，下部分為等腰梯形型狀，拋物弓頂點(diǎn)和梯形頂點(diǎn)位置溫度取特征截面xe=l/4所對應(yīng)的溫度，分別為t(l/4,0)和t(l/4,L)，此時(shí)等效溫度曲線可表示為

(2)

其中：

式(2)即為等效溫度場分布的表達(dá)式。

圖4 等效溫度場示意Fig.4 Temperature field distribution of the equivalent model

2.2 等效溫度場的驗(yàn)證

等效溫度場與真實(shí)溫度場等效的標(biāo)準(zhǔn)有2條：① 在等效截面上等效溫度場曲線與真實(shí)溫度場曲線吻合；② 等效溫度場的平均溫度與真實(shí)溫度場的平均溫度相等。第①條標(biāo)準(zhǔn)由等效溫度場曲線的定義(式(2))已經(jīng)滿足，由圖4已可見。以下對第②條標(biāo)準(zhǔn)進(jìn)行驗(yàn)證。

(1)以3排管凍結(jié)溫度場的勢函數(shù)疊加法解析解為理論基礎(chǔ)，根據(jù)凍結(jié)管布置的對稱性及溫度場在x軸上的周期性，整個(gè)凍土帷幕的平均溫度等于0≤x

(3)

(4)

圖5 凍土帷幕平均溫度計(jì)算結(jié)果比較Fig.5 Comparison of calculated average temperature

從圖6可以看出，相對誤差位于2.61%～4.95%，誤差控制在5.0%以內(nèi)，且隨著l/L的增大，誤差增大，隨著ξ/L的增大，誤差減小。

圖6 等效溫度場平均溫度的誤差Fig.6 Accuracy of average temperatures of the equivalent temperature field

3 利用等效溫度場求解平均溫度

3.1 平均溫度計(jì)算公式

上一節(jié)采用平均溫度這一指標(biāo)來驗(yàn)證等效溫度場的準(zhǔn)確性，據(jù)此，可以導(dǎo)出凍土帷幕平均溫度的近似簡單求解方法，就是利用等效溫度場來替代真實(shí)的溫度場，聯(lián)合式(2)，(4)得

(5)

其中

3.2 平均溫度計(jì)算公式的適用性

人工地層凍結(jié)過程是一個(gè)瞬態(tài)熱傳導(dǎo)過程，是個(gè)具有相變并且相變界面移動的問題。然而本文計(jì)算公式由于是基于3排管凍結(jié)穩(wěn)態(tài)溫度場解析解[15]，其適用性值得注意。對于采用穩(wěn)態(tài)溫度場近似瞬態(tài)溫度場的適用性問題，學(xué)術(shù)界和工程界普遍接受的觀點(diǎn)是：由于人工地層凍結(jié)是個(gè)發(fā)展相對緩慢的過程，尤其是在凍結(jié)的后期，其溫度場非常接近穩(wěn)態(tài)導(dǎo)熱溫度場，在人工地層凍結(jié)法的溫度場理論與工程中，對此狀態(tài)可按穩(wěn)態(tài)導(dǎo)熱近似求解人工凍結(jié)溫度場。世界上最流行和實(shí)用的人工凍結(jié)溫度場解析解有前蘇聯(lián)[1]、美國[2]和日本[3]公式，他們均為穩(wěn)態(tài)溫度場公式。對此問題，筆者曾針對文獻(xiàn)[1]的單排管和雙排管凍結(jié)巴霍爾金公式的準(zhǔn)確性進(jìn)行了分析[9]，結(jié)果表明，在凍結(jié)的中后期(ξ/l> 0.7)，溫度場任何一點(diǎn)的計(jì)算誤差不超過1℃。作為本文基礎(chǔ)的3排管凍結(jié)溫度場公式[15]與巴霍爾金公式基本特征相同，因此可以推測，在凍結(jié)的中后期，本文計(jì)算公式具有足夠的準(zhǔn)確性。

4 結(jié) 論

(1)以3排管溫度場解析解為基礎(chǔ)，建立了等效溫度場模型，即找到一個(gè)等效截面，以等效截面的溫度場替代整體溫度場。將此截面溫度場分布圖形簡化為由上下兩部分構(gòu)成的圖形，上部分為拋物弓，下部分為等腰梯形。

(2)利用平均溫度驗(yàn)證等效溫度場的準(zhǔn)確性。等效溫度場的平均溫度和解析解溫度場的平均溫度兩者的對比分析表明，等效溫度場平均溫度計(jì)算結(jié)果的相對誤差不超過5.0%，精度滿足工程要求，表明等效溫度場具有良好的實(shí)用性。

(3)基于等效溫度場，建立了3排管凍土帷幕平均溫度計(jì)算公式。該公式為免于積分計(jì)算的簡便公式，計(jì)算準(zhǔn)確，且結(jié)果偏于安全。

[1] Bakholdin B V.Selection of optimized mode of ground freezing for construction purpose[M].Moscow.State Construction Press，1963.

[2] Sanger F J.Ground freezing in construction[J].Journal of the Soil Mechanics and Foundations Division,1968,94(SM1):131-158.

[3] 戸部暢,秋元攻.凍土內(nèi)溫度分布計(jì)算式と,その応用[J].冷凍,1979,54(622):3-11. Tobe N,Akimoto O.Temperature distribution formula in frozen soil and its application[J].Refrigeration,1979,54(622):3-11.

[4] 胡向東,黃 峰,白 楠.考慮土層凍結(jié)溫度時(shí)人工凍結(jié)溫度場模型[J].中國礦業(yè)大學(xué)學(xué)報(bào),2008,37(4):550-555. Hu Xiangdong,Huang Feng,Bai Nan.Models of artificial frozen temperature field considering soil freezing point[J].Journal of China University of Mining & Technology,2008,37(4):550-555.

[5] 胡向東,白 楠,余 鋒.單排管凍結(jié)溫度場ТРУПАК和БАХОЛДИН公式的適用性[J].同濟(jì)大學(xué)學(xué)報(bào)(自然科學(xué)版),2008,36(7):906-910. Hu Xiangdong,Bai Nan,Yu Feng.Analysis of Trupak and Bakholdinformulas for temperature field of single-row-pipe frozen soil wall[J].Journal of Tongji University(Natural Science),2008,36(7):906-910.

[6] 胡向東.直線形單排管凍土帷幕平均溫度計(jì)算方法[J].冰川凍土,2010,32(4):778-785. Hu Xiangdong.Average temperature calculation for the straight single-row-pipe frozen soil wall[J].Journal of Glaciology and Geocryology,2010,32(4):778-785.

[7] 胡傳鵬,胡向東,朱合華.單排管凍結(jié)巴霍爾金溫度場控制參數(shù)敏感度分析[J].煤炭學(xué)報(bào),2011,36(8):938-944. Hu Chuanpeng,Hu Xiangdong,Zhu Hehua.Sensitivity analysis of control parameters of Bakholdin solution in single-row-pipe freezing[J].Journal of China Coal Society,2011,36(8):938-944.

[8] Hu Xiangdong.Average temperature model of double-row-pipe frozen soil wall by equivalent trapezoid method[A].API Conference Proceedings[C].2010,1233:1333-1338.

[9] 胡向東,趙俊杰.人工凍結(jié)溫度場巴霍爾金模型準(zhǔn)確性研究[J].地下空間與工程學(xué)報(bào),2010,6(1):96-101. Hu Xiangdong,Zhao Junjie.Research on precision of Bakholdin model for temperature field of artificial ground freezing[J].Chinese Journal of Underground Space and Engineering,2010,6(1):96-101.

[10] 胡向東,趙 飛,佘思源,等.直線雙排管凍結(jié)壁平均溫度的等效拋物弓形模型[J].煤炭學(xué)報(bào),2012,37(1):28-32. Hu Xiangdong,Zhao Fei,She Siyuan,et al.Equivalent parabolic arch method of average temperature calculation for straight double-row-pipe frozen soil wall[J].Journal of China Coal Society,2012,37(1):28-32.

[11] 彭益成,胡向東.丁文其.多排管凍土帷幕平均溫度特征與計(jì)算方法[J].低溫建筑技術(shù),2009,36(10):89-91. Peng Yicheng,Hu Xiangdong,Ding Wenqi.Characteritics and calculation method of average temperature for multi-row-pipe frozen soil wall[J].Low Temperature Architecture Technology,2009,36(10):89-91.

[12] 胡向東,何挺秀.多排管直線凍土墻平均溫度的等效梯形計(jì)算方法[J].煤炭學(xué)報(bào),2009,34(11):1465-1469. Hu Xiangdong,He Tingxiu.Equivalent-trapezoid method of average temperature calculation for multi-row-pipe straight frozen soil wall[J].Journal of China Coal Society,2009,34(11):1465-1469.

[13] 肖朝昀，胡向東，張慶賀.多排管局部凍結(jié)凍土壁溫度場特征[J].巖石力學(xué)與工程學(xué)報(bào)，2007，26(S1)：2694-2700. Xiao Zhaoyun,Hu Xiangdong,Zhang Qinghe.Characters of temperature field in frozen soil wall with multirow freeze-tubes and limited depth freezing[J].Chinese Journal of Rock Mechanics and Engineering,2007,26(S1):2694-2700.

[14] 陳文豹,湯志斌.潘集礦區(qū)凍結(jié)壁平均溫度及凍結(jié)孔布置圈徑的探討[J].煤炭學(xué)報(bào),1982(1):46-52. Chen Wenbao,Tang Zhibin.The average temperature in ice wall and the diameter of frozen circle in Panji coal field[J].Journal of China Coal Society,1982(1):46-52.

[15] 胡向東,汪 洋.三排管凍結(jié)溫度場的勢函數(shù)疊加法解析解[J].巖石力學(xué)與工程學(xué)報(bào),2012,31(5):1071-1080. Hu Xiangdong,Wang Yang.Analytical solution of three-row-piped frozen temperature field by means of superposition of potential[J].Chinese Journal of Rock Mechanics and Engineering,2012,31(5):1071-1080.

Atrapezoidal-parabolicarchsuperimposedequivalenttemperaturefieldmodelanditsaveragetemperatureformulafortriple-row-pipefreezing

HU Xiang-dong1,2，REN Hui1,2

(1.DepartmentofGeotechnicalEngineering,TongjiUniversity,Shanghai200092,China.2.KeyLaboratoryofGeotechnicalandUndergroundEngineeringofMinistryofEducation,TongjiUniversity,Shanghai200092,China)

The temperature field of triple-row-pipe freezing is difficult and inconvenient to be described simply due to its complexity.Based on the existing analytical solution of triple-row-pipe frozen temperature field,according to the periodicity of the temperature field,an appropriate position between the main section and the intersection was chosen to be the equivalent section,so that the true temperature field could be replaced by the temperature field of the equivalent section.Using the average temperature as equivalence criteria,a trapezoidal-parabolic arch superimposed equivalent temperature field model was proposed by simplifying the temperature field curve of the equivalent section.In addition,asimplified average temperature formula,which was noneintegral operation,was deduced from the equivalent temperature field model.The calculation results show that the equivalent temperature field and its average temperature formula are accurate enough.

triple-row-pipe freezing;equivalent temperature field;average temperature;frozen soil wall;artificial ground freezing

10.13225/j.cnki.jccs.2013.0166

國家自然科學(xué)基金資助項(xiàng)目(50578120,51178336)

胡向東(1961—)，男，浙江衢州人，副教授，博士生導(dǎo)師。Tel：021-65988771，E-mail：anton.geotech@#edu.cn。通訊作者：任 輝(1986—)，男，湖南常德人，碩士研究生。Tel：021-65985249，E-mail：ren445@126.com

TD265.3

A

0253-9993(2014)01-0078-06

胡向東，任 輝.3排管凍結(jié)梯形-拋物弓疊合等效溫度場模型和平均溫度[J].煤炭學(xué)報(bào),2014,39(1):78-83.

Hu Xiangdong，Ren Hui.A trapezoidal-parabolic arch superimposed equivalent temperature field model and its average temperature formula for triple-row-pipe freezing[J].Journal of China Coal Society,2014,39(1):78-83.doi:10.13225/j.cnki.jccs.2013.0166