不同鹽水溫度下單管凍結(jié)溫度場數(shù)值分析

吳雨薇，胡俊，汪樹成

(海南大學(xué) 土木建筑工程學(xué)院，海口 570228)

不同鹽水溫度下單管凍結(jié)溫度場數(shù)值分析

吳雨薇，胡俊*，汪樹成

(海南大學(xué) 土木建筑工程學(xué)院，海口 570228)

為了得到在不同鹽水溫度下單管凍結(jié)的溫度場分布規(guī)律，運用有限元分析軟件，采用單因素分析法進行數(shù)值模擬，得到單根凍結(jié)管在不同鹽水溫度下-1℃和-10℃等溫線的變化規(guī)律，為工程實踐提出優(yōu)化方案。主要得出：在一定的范圍內(nèi)降低鹽水溫度可以有效地擴大凍結(jié)范圍；單根凍結(jié)管的降溫過程、有效的凍結(jié)范圍(距凍結(jié)管軸心2 m處)以及-1℃與-10℃等溫線之間的距離不隨鹽水降溫計劃的改變而改變；通過擬合鹽水溫度與凍結(jié)半徑的關(guān)系得到最低鹽水溫度每降溫5℃，其-1℃等溫線的半徑增大0.109 m，而-10℃等溫線的半徑增大0.085 7 m。所得結(jié)果可為日后的相關(guān)實際工程提供參考依據(jù)。

凍結(jié)法；溫度場；數(shù)值模擬；單管凍結(jié)

0 引言

在我國的東南沿海及珠江三角洲地區(qū)，覆蓋著大量的軟土層，而近年來，隨著國民經(jīng)濟的日益發(fā)展，地下空間工程開發(fā)與利用逐漸成為與人們?nèi)粘Ｉ盥?lián)系日益緊密的項目。隧道開挖是地下工程中十分重要的一個工程分支，而在軟土地基中進行隧道開挖，會遇到很多工程問題，如土的強度不足，開挖過程中隧道滲水等，為了解決這些問題，人工凍結(jié)技術(shù)應(yīng)運而生。人工凍結(jié)技術(shù)是對軟土進行臨時加固且非常有效的一種常用加固技術(shù)手段，它是一種可逆的環(huán)境友好型的軟土加固技術(shù)，這一技術(shù)手段可以有效加強土體強度并降低土體中滲流的產(chǎn)生。人工凍結(jié)技術(shù)被廣泛運用于地下空間的開挖支護、滑坡穩(wěn)定性控制及地下水滲流的控制等。在復(fù)雜水文及工程地質(zhì)情況下的隧道開挖中，人工凍結(jié)技術(shù)常用于含水軟土中地下水的滲透控制[1-5]。因此，對人工凍結(jié)技術(shù)的研究也從未止步，由前蘇聯(lián)學(xué)者提出的經(jīng)典人工地層凍結(jié)解析解一直被沿用至今，如單管凍結(jié)溫度場公式[6]，單排凍結(jié)管溫度場公式以及雙排凍結(jié)管溫度場公式等[7]是當(dāng)今工程實踐中有意義的參照標(biāo)準(zhǔn)，國內(nèi)的學(xué)者也提出了一些相關(guān)的理論解析解答，如胡向東提出的環(huán)形單圈管凍結(jié)穩(wěn)態(tài)溫度場解析解等[8]，這些解析解公式雖然精確，但是其理論解的假設(shè)為穩(wěn)態(tài)溫度場，這樣的假設(shè)與工程實際有較大的出入，為了能夠探究凍結(jié)過程中單根凍結(jié)管周圍溫度的變化過程，本文運用ADINA有限元軟件模擬單根凍結(jié)管在不同鹽水降溫計劃下凍結(jié)半徑的發(fā)展過程，數(shù)值模型為帶相變的瞬態(tài)模型，這樣的假設(shè)更接近工程實際，因此對實際工況而言，有更強的參考意義。以此為依據(jù)應(yīng)用于工程實際中，不僅可以優(yōu)化人工凍結(jié)技術(shù)的凍結(jié)管布置方案，還可以更加準(zhǔn)確的預(yù)測凍土帷幕溫度，避免工程事故的發(fā)生。

1 建立溫度場數(shù)值模型

假定土層為單一均質(zhì)土，且符合連續(xù)、均勻、各向同性；土層凍結(jié)溫度取-1℃；忽略水分遷移的影響[9-14]。

1.1 模型的幾何尺寸

如圖1所示，基本模型為30 m(長度)×12 m(深度)的土體，在土體中設(shè)置一根圓柱形凍結(jié)管，埋深6 m，采用不同的鹽水降溫凍結(jié)計劃，分析當(dāng)最低鹽水溫度發(fā)生改變時，單根凍結(jié)管凍結(jié)范圍的變化規(guī)律。

圖1 數(shù)值模型幾何尺寸及網(wǎng)格劃分示意圖(mm)Fig.1 Schematic of geometric size and meshing diagram of numerical model

1.2 計算模型參數(shù)選取

采用帶相變的瞬態(tài)導(dǎo)熱模型作為溫度場的計算模型。網(wǎng)格劃分選取九節(jié)點網(wǎng)格劃分模式，規(guī)定單個網(wǎng)格的邊長為0.5 m。網(wǎng)格劃分后的計算模型如圖1所示。凍結(jié)的初始地溫設(shè)置為18℃(一般地層埋深10 m以內(nèi)的土體溫度為15～20℃)，以不同的鹽水降溫計劃作為溫度荷載加載到單根凍結(jié)管管壁上，不同鹽水降溫計劃見表1(凍結(jié)前5d降溫計劃一致，從凍結(jié)第10天一直到第40天，鹽水溫度一致維持在最低鹽水溫度，分別為-15、-20、-25、-30、-35、-40、-45℃這7種情況)。設(shè)置凍結(jié)時間步為40 d，24 h為一個時間步步長。模型土體采用熱導(dǎo)單元，具體參數(shù)詳見表2。

表1 不同鹽水溫度降溫計劃表Tab.1 Cooling plan for brine temperature

表2 土體材料參數(shù)Tab.2 Parameters of soils material

1.3 研究路徑劃分

為了更加翔實的研究不同鹽水溫度下，單根凍結(jié)管周圍的凍土帷幕發(fā)展情況，設(shè)置如圖1所示的觀察路徑1，在距地面以下3 m深處垂直于凍結(jié)管設(shè)置平行于水平地表的路徑1，自凍結(jié)管軸心算起，每隔0.5 m設(shè)置一個觀察點分別為1、2、3、4、5、6號，通過研究不同時間各個觀測點的溫度變化情況，分析單管凍結(jié)40 d的溫度變化全過程。

2 溫度場計算結(jié)果與分析

2.1 不同鹽水溫度下單根凍結(jié)管的凍結(jié)效果

見表1，選取的最低鹽水溫度分別為-15、-20、-25、-30、-35、-40、-45 ℃這7種情況，著重對比在這7種情況下單根凍結(jié)管的凍結(jié)半徑情況，從而優(yōu)化工程應(yīng)用中凍結(jié)管的排布問題。圖2為不同最低鹽水溫度下單根凍結(jié)管凍結(jié)40 d時的凍土帷幕發(fā)展情況。

(a)最低鹽水溫度-15 ℃

(b)最低鹽水溫度-20 ℃

(c)最低鹽水溫度-25 ℃

(d)最低鹽水溫度-30 ℃

(e)最低鹽水溫度-35 ℃

(f)最低鹽水溫度-40 ℃

(g)最低鹽水溫度-45 ℃圖2 不同鹽水溫度下凍結(jié)40 d時溫度場計算云圖與溫度等值線圖Fig.2 Band plot and contours of temperature field of various final temperature of freezing pale with the 40-day freezing

由圖2可以知道，經(jīng)過40 d的積極凍結(jié)后，凍結(jié)管周圍都形成了穩(wěn)定的凍土，凍結(jié)溫度在沿凍結(jié)管方向的分布較為一致，而在垂直于凍結(jié)管方向上的溫度分布規(guī)律隨著最低鹽水溫度的改變產(chǎn)生了差異性的凍結(jié)效果，即隨著最低鹽水溫度的降低，單根凍結(jié)管的最大凍結(jié)半徑也在擴大。觀察凍結(jié)40 d后的等值線及云圖，可以得到：在不同鹽水溫度下，從-15 ℃到-45 ℃，單根凍結(jié)管的凍結(jié)半徑(土層凍結(jié)溫度取-1 ℃)分別為0.70、0.90、1.00、1.10、1.20、1.20、1.25 m，而在不同最低鹽水溫度下得到的-10℃的等溫線半徑分別為0.25、0.45、0.60、0.75、0.80、0.90、0.90 m。

圖3 凍結(jié)40 d時單根凍結(jié)管凍結(jié)半徑 隨最低鹽水溫度變化規(guī)律示意圖Fig.3 Scatter plot of the freezing radius of the single freezing tube with the different brine temperature in the 40th days of freezing

圖3為凍結(jié)40 d時單根凍結(jié)管凍結(jié)半徑隨最低鹽水溫度變化規(guī)律，由凍結(jié)土的凍結(jié)半徑變化規(guī)律擬合出線性關(guān)系。可以看出：當(dāng)鹽水降溫計劃的最低鹽水溫度分別設(shè)定為-15、-20、-25、-30、-35、-40、-45℃時，單根凍結(jié)管在凍結(jié)40 d時-1℃與-10℃等溫線之間的距離不隨鹽水降溫計劃的改變而改變，而等溫線的半徑變化規(guī)律表明：最低鹽水溫度每降溫5℃，其-1℃等溫線的半徑增大0.109 m，而-10℃的等溫線半徑增大0.085 7 m。

2.2 路徑分析

選取如圖1所示的觀察路徑1，當(dāng)最低鹽水溫度設(shè)置為-15、-20、-25、-30、-35、-40、-45℃時，分別測量各觀測點溫度的發(fā)展情況，觀察點1～ 6平行布置，間隔0.5 m，自距凍結(jié)管軸心0.5 m處起設(shè)置。圖4為不同最低鹽水溫度下各觀測點溫度變化空間分布曲線，可知：①距凍結(jié)管2 m處開始向更遠處延伸，得到的最低鹽水溫度對土體的凍結(jié)作用開始不顯著。②若-10℃為能保證凍結(jié)質(zhì)量的最高溫度，則隨著設(shè)置的最低鹽水溫度逐漸降低，可以保證凍結(jié)質(zhì)量的凍結(jié)半徑逐漸擴大，但擴大的效果并不顯著。③在凍結(jié)位置上看，觀察點1和觀察點2上溫度的變化最為劇烈，可以考慮在重要的凍結(jié)工程作業(yè)或較為軟弱的地基土凍結(jié)過程中，在觀察點2處采取加強措施以保證凍結(jié)質(zhì)量(如適當(dāng)加密凍結(jié)管，或是在大的凍結(jié)管之間夾插小型凍結(jié)管)，在凍結(jié)期間加強溫度和壓力的監(jiān)測，以保證凍結(jié)質(zhì)量[15]。④從凍結(jié)時間上看，從凍結(jié)開始至20 d的范圍內(nèi)，鹽水降溫的效果較為顯著；而從20 d開始到40 d的范圍內(nèi)，鹽水降溫的效果并不如前20 d顯著。

(a)最低鹽水溫度為-15 ℃

(b)最低鹽水溫度為-20 ℃

(c)最低鹽水溫度為-25 ℃

(d)最低鹽水溫度為-30 ℃

(e)最低鹽水溫度為-35 ℃

(f)最低鹽水溫度為-40 ℃

(g)最低鹽水溫度為-45 ℃圖4 不同最低鹽水溫度下各觀測點 溫度變化空間分布曲線Fig.4 Scatter plot of temperature variation of observation point at different minimum brine temperatures

考慮經(jīng)濟效益，凍結(jié)時長取40 d時，建議采取-30℃的最低鹽水溫度，此時單根凍結(jié)管的有效凍結(jié)半徑(-10℃等溫線半徑)控制在0.75 m左右，可保證凍結(jié)質(zhì)量與效果。

圖5為凍結(jié)40 d時不同最低鹽水溫度下各觀測點溫度空間分布情況圖，可知：有效凍結(jié)范圍隨著鹽水降溫計劃的變化而改變，從-15℃到-35℃范圍每降低5℃時凍結(jié)有效范圍擴大較快，-35℃到-45℃范圍內(nèi)，凍結(jié)管在鹽水降溫計劃下凍結(jié)范圍擴大的速度變慢。但隨著鹽水降溫計劃的改變，Ⅱ區(qū)的土體溫度都幾乎沒有受到凍結(jié)管溫度變化的影響，說明在此區(qū)域內(nèi)，凍結(jié)管鹽水溫度的改變不會對此區(qū)域的土體強度產(chǎn)生影響，即凍結(jié)管溫度變化的影響范圍在Ⅰ區(qū)(距凍結(jié)管2 m范圍內(nèi))。

如上所述，建議采用-30℃的最低鹽水溫度積極凍結(jié)40 d，模擬實際工況，在地下插入5根間距為1.5 m的凍結(jié)管，凍結(jié)管管長為6 m，凍結(jié)效果如圖6所示。自左向右，以模擬工況中第一根凍結(jié)管向左1.5 m遠距離處為坐標(biāo)原點，得到該工況下各個時段的凍結(jié)溫度空間分布規(guī)律如圖7所示。

圖5 凍結(jié)40 d時不同最低鹽水溫度下 各觀測點溫度空間分布情況Fig.5 The spatial distribution of the temperature at each point at different minimum brine temperatures during freezing 40 days

圖6 實際工況凍結(jié)模擬云圖Fig.6 Band plot of simulate actual working conditions freezing

圖7 凍結(jié)過程面溫度分布Fig.7 Temperature distribution along different planes

由圖7模擬實驗結(jié)果可知，設(shè)置最低鹽水溫度為-30℃時，可以得到較好的凍結(jié)效果，分區(qū)對凍結(jié)效果做出測評，則Ⅱ區(qū)的凍結(jié)效果最好，凍結(jié)管周圍的溫度最低，各凍結(jié)管之間的溫度雖然有一些回升，但是總體來看，凍結(jié)管之間能夠形成完整的凍結(jié)帷幕，保證凍結(jié)土體的凍結(jié)強度。前15 d時，各個凍結(jié)管，及凍結(jié)管周圍的溫差較大，凍結(jié)管之間的溫度回升也較快。在凍結(jié)管的外側(cè)較遠的地方，溫度較高，且降溫效果并不顯著，15 d至凍結(jié)結(jié)束的時間段內(nèi)，凍結(jié)的效果較好，溫度回升較慢，凍結(jié)管之間的溫差開始變小，這對實際工程的指導(dǎo)意義為：凍結(jié)工程中為了保證凍結(jié)效果，需保證足夠的凍結(jié)時間，否則會影響完整凍結(jié)帷幕的形成。凍結(jié)管周圍的波浪形邊界也在凍結(jié)管周圍的凍土交圈完成后變得平整。針對外側(cè)凍結(jié)效果欠缺的現(xiàn)象，可以選擇加密凍結(jié)管，或者對外側(cè)凍結(jié)管實施更低的鹽水溫度的措施。

3 結(jié)束語

本文采用ADINA有限元模擬軟件對單根凍結(jié)管在不同鹽水溫度下的凍結(jié)效果做了數(shù)值模擬，主要研究了降溫計劃作為變量時，單根凍結(jié)管周圍溫度場發(fā)展及分布規(guī)律，得出：

(1)凍結(jié)管的有效凍結(jié)半徑(-10℃等溫線半徑)隨著鹽水降溫計劃中最低鹽水溫度的不斷降低而增大，但鹽水降溫計劃與有效凍結(jié)半徑并不是線性的反比關(guān)系，在-15～-35 ℃范圍內(nèi)，凍結(jié)發(fā)展較快，而到-35℃以后，對凍結(jié)有效凍結(jié)半徑的影響不大。

(2)在一定范圍內(nèi)降低鹽水溫度可以有效增大凍結(jié)管的凍結(jié)范圍，但距單根凍結(jié)管2 m處得到的最低鹽水溫度對土體的凍結(jié)作用開始不顯著。

(3)-1℃等溫線與-10℃等溫線之間的距離不隨鹽水降溫計劃的改變而改變。

(4)最低鹽水溫度每降溫5℃，其-1℃等溫線的半徑增大0.109 m，而-10℃等溫線的半徑增大0.085 7 m。

[1] 胡俊，王效賓，袁云輝.盾構(gòu)隧道端頭杯型凍結(jié)壁溫度場發(fā)展與分布規(guī)律研究[M].北京：中國水利水電出版社，2015.

[2] 胡俊.高水壓砂性土層地鐵大直徑盾構(gòu)始發(fā)端頭加固方式研究[D].南京：南京林業(yè)大學(xué)，2012.

[3] 胡俊，楊平，董朝文，等.盾構(gòu)始發(fā)端頭化學(xué)加固范圍及加固工藝研究[J].鐵道建筑，2010，15(2)：47-51.

[4] Hu J，Wang X B，Jiang B R.Numerical Analysis of Temperature Field of Vertical Frozen Soil Wall Reinforcement at Shield Shaft[C] .2014 3rd International Conference on Micro Nano Devices，Structure and Computing Systems(MNDSCS 2014)，Singapore，2014，3：4-10.

[5] 胡俊，楊平.大直徑杯型凍土壁溫度場數(shù)值分析[J].巖土力學(xué)，2015，36(2)：523-531.

[6] Трупак Н Г.Замораживание горных пород при про-ходке стволов[M].Москва：Углетехиздат，1954.

[7] Бахолдин Б В.Выбор оптимального режима замора-живания грунтов в строительных целях[M].Москва：Госстройиздат，1963.

[8] 胡向東，陳錦，汪洋，等.環(huán)形單圈管凍結(jié)穩(wěn)態(tài)溫度場解析解[J] 巖土力學(xué)，2013，34(3)：875-877

[9] 胡俊.盾構(gòu)隧道端頭垂直凍結(jié)加固不同凍結(jié)管直徑的溫度場數(shù)值分析[J].鐵道建筑，2014，19(9)：57-60.

[10] 胡俊，劉勇，曾暉.新型管幕凍結(jié)法不同管幕填充形式的溫度場數(shù)值對比分析[J].森林工程，2015，31(6)：135-141.

[11] 胡俊，衛(wèi)宏，劉勇.凍土帷幕設(shè)置加熱限位管時溫度場數(shù)值分析[J].隧道建設(shè)，2016，36(6)：688-694.

[12] 馮明揚，劉欣益，伍明強.某預(yù)應(yīng)力混凝土箱梁橋火災(zāi)下溫度場數(shù)值模擬[J].公路工程，2016，41(2)：208-212.

[13] 胡俊，劉勇.加熱限位管對凍土帷幕厚度的影響研究[J].森林工程，2016，32(2)：69-74.

[14] 胡俊.水泥改良前后土體凍結(jié)溫度及力學(xué)特性試驗研究[J].鐵道建筑，2013，18(4)：156-159.

[15] 李巖，劉波，張建新.豎向直排凍結(jié)條件水平凍脹力試驗研究[J] 巖土力學(xué)，2014，35(11)：3200-3206.

NumericalAnalysisofSingleFreezingTubeinFreezingTemperatureFieldatDifferentMinimumBrineTemperature

Wu Yuwei，Hu Jun*，Wang Shucheng

(College of Civil Engineering and Architecture，Hainan University，Haikou 570228)

In order to obtain the distribution of temperature field at different minimum brine temperature，the finite element analysis software was used to simulate the changing law of single freezing tube at different temperatures of -1 ℃ and -10 ℃ isotherm line，which could put forward optimization scheme for engineering practice.The results are as follows.In a certain range，lowering the brine temperature can effectively enlarge the freezing range.The cooling process of the single freezing tube，the effective freezing range(from frozen tube axis at 2 m)and the distance between -1 ℃ and -10 ℃ isotherms do not change with the change of the brine cooling plan.By fitting the relationship between brine temperature and freezing radius，it can be shown that when the minimum brine temperature was reduced by 5 ℃，the radius of the -1 ℃ isotherms increased by 0.109 m and the -10 ℃ isotherm increased by 0.085 7 m.The results can be used as a reference for future practical projects.

Freezing method；temperature field；numerical simulation；single tube freezing

U 455.43

A

1001-005X(2017)06-0060-07

2017-05-10

海南省教育廳高等學(xué)校科研項目(Hnky2016ZD-7，Hnky2015-10)；人社部留學(xué)人員科技活動擇優(yōu)資助項目(人社廳函[2016]176-3)；海南大學(xué)科研團隊培育專項資助項目(hdkytg201708)

吳雨薇，碩士研究生。研究方向：隧道及地下工程。

*通信作者：胡俊，博士，副教授。研究方向：隧道及地下工程。E-mail：183633299@qq.com

吳雨薇，胡俊，張友良，等.不同鹽水溫度下單管凍結(jié)溫度場數(shù)值分析[J].森林工程，2017，33(6)：60-66.