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

吳雨薇，胡俊，汪樹成

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

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

吳雨薇，胡俊*，汪樹成

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

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

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

0 引言

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

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

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

1.1 模型的幾何尺寸

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

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

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

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

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

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

1.3 研究路徑劃分

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

2 溫度場(chǎng)計(jì)算結(jié)果與分析

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

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

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

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

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

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

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

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

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

由圖2可以知道，經(jīng)過(guò)40 d的積極凍結(jié)后，凍結(jié)管周圍都形成了穩(wěn)定的凍土，凍結(jié)溫度在沿凍結(jié)管方向的分布較為一致，而在垂直于凍結(jié)管方向上的溫度分布規(guī)律隨著最低鹽水溫度的改變產(chǎn)生了差異性的凍結(jié)效果，即隨著最低鹽水溫度的降低，單根凍結(jié)管的最大凍結(jié)半徑也在擴(kuò)大。觀察凍結(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時(shí)單根凍結(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時(shí)單根凍結(jié)管凍結(jié)半徑隨最低鹽水溫度變化規(guī)律，由凍結(jié)土的凍結(jié)半徑變化規(guī)律擬合出線性關(guān)系?？梢钥闯觯寒?dāng)鹽水降溫計(jì)劃的最低鹽水溫度分別設(shè)定為-15、-20、-25、-30、-35、-40、-45℃時(shí)，單根凍結(jié)管在凍結(jié)40 d時(shí)-1℃與-10℃等溫線之間的距離不隨鹽水降溫計(jì)劃的改變而改變，而等溫線的半徑變化規(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℃時(shí)，分別測(cè)量各觀測(cè)點(diǎn)溫度的發(fā)展情況，觀察點(diǎn)1～ 6平行布置，間隔0.5 m，自距凍結(jié)管軸心0.5 m處起設(shè)置。圖4為不同最低鹽水溫度下各觀測(cè)點(diǎn)溫度變化空間分布曲線，可知：①距凍結(jié)管2 m處開始向更遠(yuǎn)處延伸，得到的最低鹽水溫度對(duì)土體的凍結(jié)作用開始不顯著。②若-10℃為能保證凍結(jié)質(zhì)量的最高溫度，則隨著設(shè)置的最低鹽水溫度逐漸降低，可以保證凍結(jié)質(zhì)量的凍結(jié)半徑逐漸擴(kuò)大，但擴(kuò)大的效果并不顯著。③在凍結(jié)位置上看，觀察點(diǎn)1和觀察點(diǎn)2上溫度的變化最為劇烈，可以考慮在重要的凍結(jié)工程作業(yè)或較為軟弱的地基土凍結(jié)過(guò)程中，在觀察點(diǎn)2處采取加強(qiáng)措施以保證凍結(jié)質(zhì)量(如適當(dāng)加密凍結(jié)管，或是在大的凍結(jié)管之間夾插小型凍結(jié)管)，在凍結(jié)期間加強(qiáng)溫度和壓力的監(jiān)測(cè)，以保證凍結(jié)質(zhì)量[15]。④從凍結(jié)時(shí)間上看，從凍結(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 不同最低鹽水溫度下各觀測(cè)點(diǎn) 溫度變化空間分布曲線Fig.4 Scatter plot of temperature variation of observation point at different minimum brine temperatures

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

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

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

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

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

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

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

3 結(jié)束語(yǔ)

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

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

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

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

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

[1] 胡俊，王效賓，袁云輝.盾構(gòu)隧道端頭杯型凍結(jié)壁溫度場(chǎng)發(fā)展與分布規(guī)律研究[M].北京：中國(guó)水利水電出版社，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] 胡俊，楊平.大直徑杯型凍土壁溫度場(chǎng)數(shù)值分析[J].巖土力學(xué)，2015，36(2)：523-531.

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

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

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

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

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

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

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

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

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

[15] 李巖，劉波，張建新.豎向直排凍結(jié)條件水平凍脹力試驗(yàn)研究[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é)?？蒲许?xiàng)目(Hnky2016ZD-7，Hnky2015-10)；人社部留學(xué)人員科技活動(dòng)擇優(yōu)資助項(xiàng)目(人社廳函[2016]176-3)；海南大學(xué)科研團(tuán)隊(duì)培育專項(xiàng)資助項(xiàng)目(hdkytg201708)

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

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

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