回填料對地埋管換熱器性能的影響

李隆鍵，鮑建鎮，廖 全，崔文智

(重慶大學 低品位能源利用技術及系統教育部重點實驗室，重慶400044)

現場巖土熱物性測試是合理設計埋管換熱器的基礎[1－2]，現在國際通用的巖土熱物性測試方法是恒熱流法，恒熱流法最先由C.Eklof等[3]于1996年提出，是IGSHPA的標準和美國采暖制冷與空調工程師學會(ASHRAE)手冊所推薦的方法。恒熱流法熱響應測試(Thermal Response Test，TRT)方法中的數據處理是基于線熱源模型[4]，它是地埋管換熱器傳熱分析的基礎，但是現有的線源法模型將管子周圍的大地土壤連同回填料部分看作是一個無限大的實體，不能考慮回填料熱物性對土壤當量導熱系數的影響。

另一方面，地埋管換熱器鉆孔埋管后回填是地源熱泵工程的重要環節，回填料的選擇和施工對地下巖土導熱系數有重要影響[5－6]。即使同一種巖石或地質成分，如果回填料的物性差別較大，測得的土壤當量導熱系數相差也比較大。S.P.Kavanaugh[7]研究表明，地下巖土導熱系數10%的偏差會造成地下埋管換熱器長度(或鉆孔深度)4.5%～5.8%的偏差，直接影響地埋管換熱器性能及系統的運行費用與初投資。然而，回填料對巖土當量導熱系數影響的定量分析至今還很少，缺乏對工程施工的理論指導。

數值模擬是研究地埋管換熱器性能的有效方法，通過實驗的驗證具有較高的精度[8－9]。本文針對實際地埋管換熱器建立土壤和冷卻介質耦合傳熱的三維瞬態模型，通過數值分析考察回填料物性參數對巖土當量導熱系數的影響。

1 物理數學模型

1.1 假設條件

由于地埋管換熱器傳熱過程的復雜性，還有地下水的影響，難以建立能反映所有實際情況的模型并求解。為簡化起見，須作以下假設：

1)土壤是均勻的，而且在整個傳熱過程中土壤的熱物性不變。

2)忽略土壤中水分遷移的影響。

3)忽略U型管管壁與回填材料、回填材料與土壤之間的接觸熱阻。

4)忽略地表溫度波動以及埋管深度對土壤溫度的影響，認為土壤初始溫度均勻一致，初始階段為當地的年平均氣溫。

1.2 物理模型

對于單U形管結構，假定地埋管換熱器內流動與傳熱都是相對于上升和下降管軸心線所在豎直剖面對稱的，因此只取整個換熱器結構一半作為計算區域。地埋管鉆孔深度100m，土壤半徑3m，U型管是外徑32mm，壁厚3mm，內徑是26mm的高密度聚乙烯管，U型管上升和下降管中心距為64mm，鉆孔直徑為130mm，U型彎底端再延伸3m作為外層邊界。簡化模型示意圖如圖1所示。很顯然，在TRT測試中地埋管換熱器內是一個三維非穩態對流－導熱耦合傳熱過程。

1.3 控制方程

由于地埋管換熱器中流體的流動一般都處于湍流狀態，因此在進行CFD模擬時，選用了κ－ε雙方程模型[10]來解得U型管中水的湍流流動和傳熱，并將管內對流換熱與回填土和土壤中的導熱過程耦合起來。管內湍流流動的控制微分方程，包括連續性方程、動量方程、能量方程、湍動能方程以及湍流耗散方程依次為：

圖1 幾何模型示意圖

回填土和土壤傳熱的控制方程：

式中：ρ為流體的密度，kg·m－3；u為速度矢量：μeff為有效黏性系數，Pa·s；p為流體的壓力，Pa；h為流體的焓，J·kg－1；k為湍流脈動動能，μ為流體的黏性系數，Pa·s；μt為湍流黏性系數，Pa·s；ε為湍流脈動動能耗散率；σT、σk、σε、σε1、σε2為常數；ρs為回填土或土壤密度，kg·m－3；c為回填土或土壤的比熱容，J·kg－1·K－1；T為回填土或土壤溫度，K；λ為回填土或土壤的熱導率，W·m－1·K－1。

1.4 數值求解

對上述非穩態流固耦合傳熱問題，采用控制容積有限元方法(CVFEM)進行求解，控制方程中的對流項全部采用二階迎風差分離散。為了保證計算良好的收斂性，壓力和動量采用欠松弛的迭代法。

在劃分網格時，在流體的近壁面處加邊界層，在U形管下部的彎管處，由于流場變化劇烈，且曲度較大，所以要對該部分網格加密，避免網格有較大的傾斜角，整個模型均用六面體網格。在U型管進口處定義為VELOCITY＿INLET邊界類型，出口處為OUTFLOW邊界類型。用恒熱流法做TRT熱響應測試時，水是在一個封閉的系統中流動，當水從U型管的出口出來時，會經過一個加熱器，使水有一個ΔT的溫升，為了實現此功能在U型管的進出口用UDF編一個周期性邊界條件。由于模型的長寬比較大，所以選擇雙精度的求解器。

2 模擬結果及分析

2.1 模型的驗證

目前，土壤導熱系數最主要的測試方式是現場熱響應實驗測試法[11]。該實驗法要求測試的鉆孔結構盡量與工程實際一致，循環工質、導管尺寸、回填料、井孔的深度與半徑等都要按照工程設計要求確定。近年來，瑞典、德國和美國等分別研制出了移動的地下熱響應實驗測量設備[12－13]。筆者參照熱響應測試系統和方法，開發出巖土熱物性現場測試儀[14]。測量儀主要部件由加熱器、循環水泵、流量測量裝置、溫度測量裝置、信號變送裝置、微機控制與處理裝置等構成。測量儀中的管路與地埋管換熱器地下回路相接，循環水泵驅動流體在回路中循環流動，流體經過加熱器加熱后流經地下回路與地下巖土進行換熱。為了檢驗數值分析模型及其數值計算方法的正確性，針對某地埋管換熱器實驗測試工況進行了相應的計算，計算所需材料物性參數見表1。

表1 材料物性參數表

管內流體、管壁、回填土、土壤的初始溫度均為18℃。由TRT熱響應測試實驗測得，在加熱功率為6kW，流速為0.48m／s下，測得土壤的導熱系數為2.7W·m－1·K－1，由此在模型中設土壤和回填土的導熱系數均為2.7W·m－1·K－1，在定加熱功率6kW、流速為0.48m／s的條件下，流體進入加熱器后會有一個ΔT的溫升，由計算得到為5.57℃，可以通過FLUENT的UDF來實現此功能。在實驗和模擬的情況下，地埋管換熱器連續運行76h，U型管進出口算術平均水溫隨時間的變化曲線如圖2所示。

由于實驗時每10s記錄一次數據，數據量較大，所以在圖上呈上下震蕩的曲線，而模擬的數據400s記錄一下數據，呈現出比較規律的曲線。由圖2可以看出，模擬結果與實驗結果基本吻合，進出口平均溫度最大相對誤差為2.6%。由此可見，數值計算可以再現TRT測試的整個過程，證明了所建立的模型及其數值計算方法是正確的。

圖2 U型管進出口平均水溫模擬值與實驗值的比較

2.2 回填土導熱系數的影響

描述回填材料熱物理性質的基本參數有導熱系數、密度、比熱容等，除地質條件中地下水對換熱器性能影響重大外，其中導熱系數是最關鍵的，也是決定系統效率高低的主要因素。圖3為地埋管換熱器連續運行70h，土壤導熱系數為2.7W／(m·K)，回填土導熱系數為1W／(m·K)和4W／(m·K)時鉆孔內部50m深度的溫度分布。從圖中可以看出，在不同導熱系數下，管內流體的溫度相差4℃，可見回填土導熱系數對地埋管換熱器的影響還是比較大。

圖3 回填土導熱系數為1W／(m·K)和4W／(m·K)時50m深度的溫度分布

為了考察回填料導熱系數對TRT當量導熱系數的影響，取不同的回填料導熱系數進行熱響應測試過程的數值模擬，然后根據數值模擬計算結果按照TRT測試的數據處理方法得到土壤當量導熱系數。圖4為土壤導熱系數分別為2.2W／(m·K)和2.7W／(m·K)時，土壤當量導熱系數隨回填土導熱系數的變化曲線。由圖可見，當回填料導熱系數小于土壤導熱系數時，TRT得到的當量導熱系數受回填料導熱系數的影響非常顯著。此時，隨著回填料導熱系數的減小，當量導熱系數呈指數規律遞減。當回填土導熱系數為0.4W／(m·K)時，計算所得土壤當量導熱系數比土壤的真實導熱系數2.2W／(m·K)和2.7W／(m·K)分別要小10.5%和8%。

圖4 土壤當量導熱系數隨回填土導熱系數的變化曲線

當回填料導熱系數大于土壤導熱系數時，TRT測試的當量導熱系數已非常接近土壤導熱系數。理論上分析，隨著回填土導熱系數的增大，計算所得結果應該分別逐漸趨近于土壤的導熱系數。這是因為當回填土導熱系數增大時，鉆孔熱阻減小，當導熱系數無限大時，鉆孔熱阻就會無限小，此時就可以認為整個地埋管換熱器的熱阻只有土壤部分起作用，而導熱系數就趨近于土壤的導熱系數。為了在TRT熱響應測試中能夠得到更真實的土壤的導熱系數，應該盡可能增大回填土的導熱系數，但是當回填土的導熱系數超過土壤導熱系數的時候，得到的土壤當量導熱系數基本不變，所以也不能一味的增大其導熱系數。

2.3 回填土比熱的影響

回填土比熱也是影響TRT熱響應測試結果的重要因素，比熱是單位質量物質的熱容量，比熱越大，對熱的容量越大，熱慣性越大，對土壤和冷卻介質之間傳熱的阻滯作用越大。圖5為地埋管換熱器連續運行70h，土壤導熱系數為2.7W／(m·K)回填土比熱為620J／(kg·K)和3 620J／(kg·K)時50m深度的溫度分布，從圖中可以看出回填土比熱較大時管內水的溫度相對較高。

圖5 回填土比熱為620J／(kg·K)和3 620J／(kg·K)時50m深度的溫度分布

從文獻[15]可知，回填土比熱與鉆孔熱阻有很強的正相關性，即回填土比熱越大，相當于鉆孔熱阻越大，反之，回填土比熱越小，鉆孔熱阻越小。從圖6可以看出，土壤導熱系數分別取2.2W／(m·K)和2.7W／(m·K)，回填土比熱大于1 000J／(kg·K)時，計算得到的土壤的當量導熱系數是隨著回填土比熱的增大線性減小的，且下降幅度基本一致，當回填土比熱小于1 000J／(kg·K)時，隨著回填土比熱的減小，土壤的當量導熱系數逐漸趨近于土壤的真實導熱系數2.2W／(m·K)和2.7W／(m·K)。在土壤的2種導熱系數下，當回填土比熱為3 620J／(kg·K)時，計算所得土壤當量導熱系數比土壤的真實導熱系數分別要小8%和7%，可見回填土比熱對TRT熱響應測試結果還是有一定影響的。

圖6 土壤當量導熱系數隨回填土比熱的變化曲線

因為當回填土比熱增大時，熱慣性增大，相當于鉆孔熱阻增加，相應的土壤的當量導熱系數減小，所以為了在TRT熱響應測試中能夠得到更真實的土壤的導熱系數，應該盡可能減小回填土的比熱。

3 結 論

應用地埋管換熱器三維非穩態耦合傳熱的物理數學模型，分析了回填土導熱系數和比熱對TRT熱響應測試得到的土壤當量導熱系數的影響，得出了如下結論：

1)地埋管換熱器中當回填料導熱系數小于土壤導熱系數時，土壤當量導熱系數隨回填料導熱系數的減小呈指數規律減小。當回填料導熱系數大于土壤導熱系數時，土壤當量導熱系數幾乎不再變化；設計或施工過程中選擇回填料的導熱系數至少不得低于土壤導熱系數。

2)土壤的當量導熱系數隨著回填土比熱的增大而線性減小。在其它條件相同的情況下，當回填土比熱范圍從620J／(kg·K)增大到3 620J／(kg·K)時，土壤當量導熱系數減小7%～8%；實際TRT測試和工程施工中測試中，應盡量減小回填料比熱。

[1]NINA MATTSSON，GILBERT STEINMANN，LYESSE LALOUI.Advanced compact device for the in situ determination of geothermal characteristics of soils[J].Energy and Buildings，2008，40：1344－1352.

[2]方亮，張方方.關于地埋管換熱器熱響應試驗的討論[J].建筑熱能通風空調，2009，28(4)：48－51.FANG LIANG，ZHANG FANG－FANG.Comments on the thermal response test of the ground heat exchangers[J].Building Energy ＆ Environment，2009，28(4)：48－51.

[3]ESEN H，INALLI M.In－situ thermal response test for ground source heat pump system in elazi｛dotless｝g，Turkey[J].Energy and Buildings，2009，41(4)：395－401.

[4]KAVANAUGH S P.Field tests for ground thermal properties Inethods and impact on ground source heat pump design[J].ASHRAE Transactions，1992，98(2)：607－615.

[5]高青，余傳輝，馬純強，等.地下土壤導熱系數確定中影響因素分析[J].太陽能學報，2008 29(5)：581－585.GAO QING， YU CHUAN－HUI， MA CHUANQIANG，et al.Analysis of influence factors on determining the ground thermal conductivity[J].Acta Energiae Solaris Sinica，2008，29(5)：581－585.

[6]趙進，王景剛.地源熱泵土壤熱物性測試與分析[J].河北工程大學學報：自然科學版，2010，27(1)：58－60.ZHAO JIN，WANG JING－GANG.Testing and analysis on the soil thermal properties by ground source heat pump[J].Journal of Hebei University of Engineering：Natural Science Edition，2010，27(1)：58－60.

[7]KAVANAUGH S P.Field test s for ground thermal properties methods and impact on ground－source heat pumps[J].ASHRAE Trans，1998，104(2)：347－355.

[8]HAKAN DEMIR，AHMET KOYUN，GALIP TEMIR.Heat transfer of horizontal parallel pipe ground heat exchanger and experimental verification[J].Applied Thermal Engineering，2009，29(2／3)：224－233

[9]SHARQAWY M H，MOKHEIMER E M，HABIB M A.Energy，exergy and uncertainty analyses of the thermal response test for a ground heat exchanger[J].International Journal of Energy Research，2009，33(6)：582－592.

[10]陶文銓.數值傳熱學[M].2版.西安交通大學出版社，2001：333－385.

[11]朱瑞琪，張健明.熱響應試驗方法與機理[J].西部制冷空調與供暖，2004，24(2)：1－3.ZHU RUI－QI，ZHANG JIAN－MING.Method and mechanism of thermal response test[J].Western Refrigeration，Air Conditioning and Heating，2004，24(2)：1－3.

[12]MOGENSEN P.Fluid to duct wall heat transfer in duct system heat storages[C]／／Proceedings of the International Conference on Subsurface Heat Storage in Theory and Practice.Swedish Council for Building Research，1983.

[13]SIGNHILD GEHLIN.Thermal response test－in situ measurements of thermal properties in hard tock[D].Hdea University of Technology，Sweden，1998.

[14]彭清元，常桂欽，崔文智，等.基于線熱源理論的巖土熱響應測試研究現狀[J].節能，2010，(8)：15－19.PENG QING－YUAN，CHANG GUI－QIN，CUI WENZHI，et al.Recent status of in－situ thermal response test based on line－source theory[J].Energy Conservation，2010(8)：15－19.

[15]MARCOTTE D，PASQUIER P.On the estimation of thermal resistance in borehole thermal conductivity test[J].Renewable Energy，2008(33)：2407－2415.