翟國兵，趙大軍，張金寶，劉玉民，孫梓航，劉華南

(吉林大學建設工程學院，吉林長春130026)

0 引言

在油頁巖原位高溫開采過程中，需建立安全可靠的地下冷凍墻以隔絕周圍的地下水滲入，防止油氣泄漏［1］。地下冷凍墻還可用于礦山豎井開挖、基坑臨時支護等領域，目前在豎井開挖中應用最為廣泛。冷凍墻形成凍結的過程為:凍結管周圍的地層凍結，形成圓柱體，凍結柱體沿徑向擴展，在相鄰凍結體之間形成凍結交圈［2］。最后，凍結交圈發(fā)展成為冷凍墻屏障。

地下冷凍墻的制冷液溫度影響凍結交圈時間、制冷效率及成本。因此，本文基于有限元模擬，分析制冷液溫度對地下溫度場的影響規(guī)律，優(yōu)化在不同制冷方式下的制冷液溫度，有利于提高制冷效率、節(jié)約能源、減少施工成本和運行費用，為地下工程建設尤其是油頁巖原位開采提供技術支撐，具有非常重要的理論及現(xiàn)實意義。

在應用地下冷凍墻于油頁巖原位開采方面，國內相關研究報道極少，國外的殼牌公司對此研究相對成熟［3－4］。關于人工凍結，國內外學者主要集中于凍結井結構設計、凍結土物理力學性質及凍土帷幕溫度等方面研究。李功洲對現(xiàn)場凍結壁溫度場和位移進行實測研究，認為降低冷凍液的溫度可提高凍結壁的擴展速度，提高強度，同時可抑制深井壁的位移過大［5］。李述訓等基于數(shù)學建模描述了巖土凍結過程中溫度場變化規(guī)律，并結合實測資料研究了相鄰凍結管之間溫度場擴展變化及外側凍結壁的形成變化規(guī)律［6－7］。

但國內研究僅限于分析模擬凍結過程中整個凍結井溫度場的變化和凍結交圈規(guī)律，從理論上預測凍結井溫度場的影響因素，并未在明確控制其它變量的條件下研究制冷液溫度對制冷效果的影響規(guī)律。本文基于有限元分析方法，對制冷液溫度對地下冷凍墻制冷效果的影響規(guī)律進行了分析研究，并以提高效率、降低成本為目標，對制冷液溫度進行優(yōu)化。

1 基本理論及模型建立

1．1 假設條件與理論基礎

為便于分析求解地下?lián)Q熱器復雜的換熱過程，作出以下假設:

(1)假設在傳熱過程中，土體物理成分、熱物性參數(shù)保持不變;

(2)假設換熱器內同一截面流體速度、溫度分布均勻;

(3)假設熱量傳遞在土體中沿徑向和垂直方向，忽略圓周方向的導熱過程;

(4)假設回填土與孔壁、換熱器與回填土接觸完好，忽略接觸熱阻;

(5)假設土體與換熱器之間熱量傳遞為純導熱的傳熱過程，忽略水分遷移的影響;

(6)在所建求解區(qū)域物理模型中，假設求解區(qū)域半徑為遠邊界半徑，由遠邊界理論可知，為確保計算精度，避免受熱流擴散影響，應將遠邊界半徑選取的足夠大［8］，此處選取為10 m。

根據以上假設條件，建立地面以下徑向二維瞬態(tài)的傳熱模型如圖1，管外側土體熱傳導是徑向一維不穩(wěn)定傳熱方式，其在平面直角坐標系內滿足下列熱傳導微分方程:

圖1 求解區(qū)物理模型

式中，r——土體半徑，m;a——土體導溫系數(shù)，m2/s;k——土體導熱系數(shù)，W/(m·℃);θ——過余溫度(θ=t－ t0)，℃;t0——初始溫度，℃;θf——管內載冷劑的過余溫度(θf=tf－ t0)，℃;tf——載冷劑在換熱器內平均溫度，℃;α——載冷劑與孔壁間對流換熱系數(shù)，W/(m·℃);τ——時間，s。

1．2 傳熱初始條件與邊界條件

初始溫度值為在實驗系統(tǒng)未運行時，在換熱器中注入制冷液后，待地下溫度穩(wěn)定時所測溫度值［9－10］。由實驗現(xiàn)場數(shù)據得，地下4 m以深初始溫度平均值為8.25℃，所以取該值為模擬的初始溫度值［11］:

熱流密度Q在地下?lián)Q熱器內壁均勻加載，忽略循環(huán)介質影響，低溫循環(huán)液的平均溫度即為凍結管內壁平均溫度，即:

式中:Q——熱流密度，J/(m2·s);ks——換熱器導熱系數(shù)，W/(m·℃);n——經過該點等溫線上的法向單位矢量;?！叵侣窆軗Q熱器的內壁。

2 仿真模擬與分析

2．1 模型建立與網格劃分

本文采用ANSYS分析軟件模擬不同制冷工況下的地下土體溫度場變化情況，采用 ANSYS PLANE55單元，單元中每個節(jié)點只有一個自由度，即為溫度。采用自動劃分網格，凍結管及整個求解區(qū)域的網格劃分情況如圖2所示。對于16個同軸管換熱器區(qū)域，由于溫度變化明顯，溫度梯度較大，因此此區(qū)域網格劃分細密，網格最大邊長5 mm，網格精度為5;對于其它區(qū)域，劃分網格時采用按比例增大，控制網格最大邊長50 mm。

2．2 模擬參數(shù)設計

在仿真分析溫度對凍結效果的影響時模擬相關參數(shù)如表1所示。

圖2 地層凍結區(qū)域網格劃分

表1 有限元分析模擬參數(shù)

獲取數(shù)據時，選取5個溫度采集孔。如圖3所示，沿凍結管排布圓周，距凍結管內側1、2 m處分別為2號、1號采集點，距凍結管外側1、2 m處分別為3號、4號采集點，凍結孔中間為5號采集點。

2．3 制冷液溫度對凍結效果的影響分析

分別在制冷液溫度為 －5、－10、－15、－20℃時，模擬地下制冷系統(tǒng)運行180 d后的土體凍結情況，并得到1～5號數(shù)據采集點的溫度云圖和溫度變化曲線，如圖4～7所示。

圖3 測溫孔分布示意圖

圖4 溫度為－5℃的溫度云圖和各個測溫點的溫度變化曲線

由圖4～7可知，在凍結時間15.6×106s(4320 h)后:制冷液溫度為－5℃時，1號、2號測溫孔溫度下降速度分別快于3號、4號，但均未達到0℃，凍結圈沿內擴展約0.5 m，外側約0.3 m，均未擴展到1 m，位于凍結管中間的5號測溫孔溫度降低速度最快，最終下降到－0.4℃，最先進入凍結階段，但并未達到完全凍結溫度，仍需一段時間的凍結才能完全交圈。制冷液為－10℃時，除4號外其他各孔均進入凍結階段，凍結管內側凍結圈已擴展至1.6 m，外側不足1 m，約為0.7 m。同時，5號測溫孔曲線末溫度達到－3.67℃，說明相鄰凍結管中間地層已形成有效凍結交圈，即已經完全凍結。制冷液為－15℃時，內側擴展寬度約3 m，外側約1.2 m。凍結圈均擴展1 m以上，內側寬度大于外側。5號測溫點在凍結4.7×106s(1319 h)后，溫度下降至－2℃，完成交圈，并繼續(xù)凍結，進一步擴展交圈寬度。制冷液為－20℃時，內側凍結圈擴展至5 m，外側擴展至1.8 m。5號測溫點在凍結3.2×106s(880 h)時，溫度已低于－2℃，完成交圈。

圖5 溫度為－10℃的溫度云圖和各個測溫點的溫度變化曲線

圖7 溫度為－20℃的溫度云圖和各個測溫點的溫度變化曲線

在制冷液溫度為－5～－20℃條件下，模擬得到各測溫孔凍結15.6×106s(4320 h)后的溫度值和溫度變化幅度，統(tǒng)計如圖8、圖9所示。

由圖8、圖9可知，溫度曲線基本呈直線下降，1～5號測溫孔末溫度值隨制冷液溫度降低而降低，而1～5號測溫孔的降溫幅度隨著制冷液溫度降低基本呈增大趨勢，由此說明在－5～－20℃，隨制冷液溫度降低，地層制冷溫度逐漸降低，制冷效率提高。

圖8 不同制冷液溫度下各個測點溫度值曲線

圖9 不同制冷液溫度下各個測點溫度降低值曲線

完全凍結溫度為－2℃，當制冷液溫度為－5℃時，1～5號測溫孔均未達到完全凍結。制冷液－10℃時，只有5號和1號達到完全凍結階段。制冷液－15℃時，除4號測溫孔仍在0℃以上，其他各孔均達到完全凍結。制冷液－20℃時，除4號以外各測溫孔仍均處于完全凍結狀態(tài)，且各自溫度均低于－15℃時所測，4號溫度雖降低至－0.45℃，但未完全凍結。從－5～－15℃，達到完全凍結狀態(tài)的測溫孔逐漸增多，說明完全凍結區(qū)域面積(－2℃以下)隨制冷液溫度降低而擴大。而5號、1號、2號、3號依次達到完全凍結狀態(tài)，說明凍結圈在凍結管內側擴展速度快于外側，并且由于地層凍結交圈作用，位于兩相鄰凍結管中間地層凍結速度最快，凍結效率最高。

在同一制冷液溫度下，5號溫度值均為最低，1號、2號測溫曲線溫度值低于3號、4號溫度值，且1號溫度均低于2號，3號均低于4號;5號溫度降幅最大，1號、2號降幅大于3號、4號，1號降幅大于2號，3號大于4號。說明沿凍結管所排布的圓周，其內側的制冷效果和速率強于外側，并隨距凍結管中心距離增大而減弱。制冷受凍結交圈作用影響，相鄰凍結管之間地層制冷效果最好，效率最高。

同時可以看出，制冷液溫度由－15℃降低為－20℃時，雖然各個測溫孔溫度均降低，但達到完全凍結狀態(tài)的測溫孔個數(shù)并未增多。而圖9中4號測溫孔溫度降低幅度在－5～－10℃增長緩慢，在－10～－15℃增長迅速，在－15～－20℃溫度降低幅度有輕微下降，說明地層制冷效果隨制冷液溫度的降低而增強，但過低的制冷液溫度對凍結管外圍的地層的制冷效果影響不大。

3 冷凍液溫度的優(yōu)化

為優(yōu)化地下冷凍墻的制冷液溫度，將制冷液溫度分別設定為－5、－7.5、－10、－15、－20、－25℃進行模擬凍結，得到不同制冷液溫度下的凍結交圈時間。制定以制冷液溫度為橫坐標，交圈時間為縱坐標的坐標系，將以上點繪制在坐標系中，并將其擬合成一條曲線，如圖10所示。

圖10 冷凍液溫度與交圈時間趨勢

由圖10可知，總體上，從－5～－25℃，交圈時間隨制冷液溫度的降低而減少，但在不同的溫度段降低的趨勢有所不同。

溫度由－5～－10℃，曲線基本呈直線狀下降趨勢，交圈時間隨溫度降低而迅速減少;溫度由－10～－15℃，曲線呈凹形，交圈時間隨溫度降低而減少的趨勢放緩;溫度由－15～－25℃，曲線更加平緩呈近似水平直線狀態(tài)，交圈時間隨溫度的降低變化很小，基本趨于穩(wěn)定。

分析產生以上趨勢的原因為:在制冷液溫度較高的條件下，地層與制冷液溫度梯度較小，使得地層溫度降低速度較慢。隨著制冷系統(tǒng)的運行，地層溫度逐漸下降，當其降低至0℃以下而開始凍結時，溫度梯度更小，溫度下降速率更慢，凍結時間延長，因此交圈時間較長。而在制冷液溫度過低的條件下，地層與制冷液溫度梯度雖較大，但由于凍結過程受地層熱傳導系數(shù)限制，溫度下降速率達到最大值并保持穩(wěn)定，因此凍結交圈時間值雖小卻不能無限縮短，最終趨于穩(wěn)定值。

通過對可變成本中的制冷功率及電費、設備折舊費、材料費、鉆進費、檢查費等費用的綜合分析，將不同制冷液溫度方案的成本整理并繪制曲線如表2及圖11所示?？梢钥闯?，成本隨制冷液溫度的降低而增大，在－5～－7.5℃，曲線較平緩，成本隨制冷液溫度降低變化幅度不大;在－7.5～－15℃，曲線變化率有所增大，成本隨制冷液溫度降低呈線性增長;由－15℃降至－20℃時，曲線變化幅度陡然變大，成本隨溫度的降低而增長的幅度最大;由－20～－25℃，曲線最為平緩，隨溫度的降低，成本變化很小。

表2 不同制冷液溫度的成本分析

圖11 不同制冷液溫度的成本曲線

因此從交圈時間曲線可知，當制冷液溫度高于－15℃時，隨制冷液溫度降低，交圈時間有較大的降低空間，凍結速率提高較快，而當制冷溫度低于－15℃時，隨制冷溫度降低，交圈時間基本趨于穩(wěn)定，變化幅度極小，凍結速率提高空間極小。同時由成本曲線可看出，若采用制冷機組制冷，制冷液溫度越低，所需制冷機組功率越大，當制冷液溫度低于－15℃時，設備的投資與運行成本會大幅增加;若采用風冷機組制冷，制冷液溫度受制于環(huán)境溫度影響，利用東北地區(qū)冬季自然冷源，環(huán)境溫度平均溫度為－15℃，制冷液平均溫度可達－10℃。因此綜合考慮各種因素條件，雖然由前文可知隨制冷液溫度下降，地層制冷溫度降低，但過低的制冷液溫度對凍結交圈速度的提高影響不大，并且一味追求更低的制冷溫度還會受到成本和環(huán)境等因素的制約，因此，采用制冷機組制冷或風冷制冷方式，選取制冷液溫度在－10～－15℃較為合適。

4 結論

(1)以地下冷凍墻實驗平臺為模擬對象，使用ANSYS數(shù)值模擬軟件，進行了制冷液溫度對凍結效果影響的相關模擬分析，得到了不同制冷液溫度下的溫度場。

(2)在綜合考慮地下冷凍墻交圈時間、制冷效率、工程成本等因素，對制冷液溫度進行優(yōu)化，得出在制冷機組或風冷制冷方式下，制冷液溫度應在－10～－15℃較為合適。

雖已進行了很多研究工作，但受限于軟件模擬精度、網格劃分細密程度等因素，仍存在不足之處。在之后的研究工作中，為提高結論精確性，進行不同制冷液溫度的模擬時，要減小模擬溫度間距，并著重進行制冷溫度在－10～－15℃區(qū)間的凍結模擬，進一步提高優(yōu)化值的精度。

［1］ 李淑敏，趙大軍，劉玉民，等．油頁巖原位開采地下冷凍墻緩沖距離研究［J］．探礦工程(巖土鉆掘工程)，2013，40(12):9－12．

［2］ 張金寶．油頁巖原位開采地下冷凍墻聯(lián)合制冷系統(tǒng)的實驗研究［D］．吉林長春:吉林大學，2014．

［3］ 段會軍．地下冷凍墻換熱器的研究［D］．吉林長春:吉林大學，2012．

［4］ 于磊．油頁巖原位開采地下冷凍墻制冷系統(tǒng)的實驗研究［D］．吉林長春:吉林大學，2012．

［5］ 李功洲．深井凍結壁位移實測研究［J］．煤炭學報，1995，20(1):99－104．

［6］ 李述訓．立井凍結法鑿井工程中的熱工計算［J］．冰川凍土，1994，16(1):21 －30．

［7］ 南卓銅，李述訓，程國棟，等．地面凍結數(shù)模型及其在青藏高原的應用［J］．冰川凍土，2012，(1):89 －95．

［8］ 趙研．太陽能地下混凝土存取熱試驗研究與數(shù)值模擬［D］．吉林長春:吉林大學，2011．

［9］ 張金寶，趙大軍，于磊，等．嚴寒地區(qū)自然冷源用于冷凍墻制冷系統(tǒng)［J］．吉林大學學報(地球科學版)，2012，(S3):329 －336．

［10］ 李淑敏．基于自然冷源的地下冷凍墻制冷系統(tǒng)研究［D］．吉林長春:吉林大學，2014．

［11］ 王?。晚搸r原位開采溫度場的數(shù)值模擬［D］．吉林長春:吉林大學，2011．