劉逸 ，陳培強(qiáng)

（1. 哈爾濱商業(yè)大學(xué)能源與建筑工程學(xué)院，黑龍江哈爾濱 150028）

隨著建筑能耗的日益增加，建筑節(jié)能技術(shù)日漸被我國所了解與重視，利用可再生能源，采用土壤源熱泵系統(tǒng)取代以往的一次能源對建筑進(jìn)行供熱供冷的方式為節(jié)能方面提供了一個新的思路［1］。近年來，針對地埋管鉆井回填材料進(jìn)行的研究，得出了回填材料的導(dǎo)熱系數(shù)大可以促進(jìn)換熱效果的結(jié)論［2］。通過加快管內(nèi)流速，管道與周邊土壤的傳熱增加，但進(jìn)出口水溫溫差減少［3］。Koyun等人對埋管材質(zhì)進(jìn)行研究，認(rèn)為埋管材料選用鋁管時，較普通塑料管可大幅度提升埋管換熱量［4］。Hikmet Esen等人對垂直埋管埋深（30米、60米、90米）與周邊巖土溫度變化進(jìn)行實驗，證明井深對巖土體冷熱堆積無明顯影響［5］。G.Florides等人［6］進(jìn)行巖土豎直方向分區(qū)，對水平及垂直埋管方式進(jìn)行比較，得出垂直埋管方式換熱效果更好的結(jié)論。由于土壤源熱泵的熱源是土壤，所以長期運行勢必帶來土壤熱失衡問題［7］，進(jìn)而在地面上生長的植物也會因土壤溫度的改變而受到一定的影響［8］。

綜上所述，針對土壤源熱泵系統(tǒng)中地埋管換熱器的各項參數(shù)對換熱器換熱性能影響方面的研究較為豐富，但對于考慮有無地下水滲流對換熱器換熱性能影響的研究還很少。本文以大連某一別墅用土壤源熱泵實際項目為研究對象，依據(jù)實際工程建立數(shù)值計算三維模型，模型尺寸、模擬過程中所涉及到的各項物性參數(shù)及循環(huán)流體入口溫度、流速均由實際工程熱響應(yīng)實驗測得，通過將埋管換熱器出口溫度實測值與模擬值的對比驗證所建模型可靠性。當(dāng)?shù)芈窆茉诠┡r下連續(xù)運行一周時，研究其有無滲流對系統(tǒng)換熱及周邊巖土溫度變化的影響。得出有地下水滲流時可有效減少土壤冷堆積問題并能使埋管換熱效率提升的結(jié)論，有利于為實際工程項目減少不必要的初投資，同時為后續(xù)在大連開展小型建筑土壤源熱泵采暖方式的廣泛性研究奠定基礎(chǔ)。

1 地埋管換熱器三維傳熱數(shù)值模型

1.1 數(shù)學(xué)模型

1.1.1 鉆孔外數(shù)學(xué)模型

將巖土體近似看做物性參數(shù)符合各向同性的理想多孔介質(zhì)，土壤最初溫度T0，豎直方向上溫度變化量不計，假設(shè)地下水滲流沿水平方向流動，因其彎管處尺寸很小所以可視地埋管為有限長線熱源，循環(huán)流體的物性參數(shù)決定了埋管在巖土體中的熱流密度。則地埋管傳熱方程為［9］,［10］:

其初始條件及邊界條件為；

其中，T為溫度，℃；U為當(dāng)量滲流速度，m·s-1；a為熱擴(kuò)散系數(shù)，m2·s-1；τ為時間，s；r′為徑向距離；r為坐標(biāo)系空間距離；cp為 流體比熱容，J·kg-1·K-1；T為埋管進(jìn)口水溫，℃；T為埋管inout出口水溫，℃；M為埋管內(nèi)流體的質(zhì)量流速，kg/s；H為鉆孔深度，m。

基于虛擬熱源法，巖土體的瞬時溫度解析解為［11］：

其中

當(dāng)時間趨于零時，地下巖土體溫度解析解為

1.1.2 鉆孔內(nèi)數(shù)學(xué)模型

由于鉆孔內(nèi)流體與回填材料熱容比土壤小很多，所以鉆孔內(nèi)對流換熱及導(dǎo)熱可近似看做穩(wěn)態(tài)傳熱問題，系統(tǒng)運行時管內(nèi)流體是時刻都在運動的，溫度也在不斷地降低或提高，根據(jù)能量守恒，得出流體與管壁換熱可以采用以下兩個方程［12］：

其中定解條件為：

Tf1為 下行管內(nèi)流體溫度，℃；Tf2為上行管內(nèi)流體溫度，℃。

鉆孔內(nèi)熱阻表達(dá)式為［13］

式中R11為流體進(jìn)口所在管路與流體的熱阻、R22為 流體出口所在管路與流體的熱阻、R12為兩支管間熱阻，m·K·W-1；λg為回填導(dǎo)熱系數(shù)、λs為巖土導(dǎo)熱系數(shù)、λp為管的導(dǎo)熱系數(shù)，W·m-1·K-1；rb為鉆孔半徑、ri為聚乙烯管的內(nèi)半徑、ro為 聚乙烯管的外半徑，m；D為聚乙烯管支管管間距，m；根據(jù)以下公式確定流體對流換熱系數(shù)hi［14］：

其中，供熱時n=0.4，制冷時n=0.3。

以上方程組用數(shù)學(xué)方法求解，簡化得出進(jìn)出口溫度與壁溫的函數(shù)關(guān)系式：

1.1.3 多孔介質(zhì)數(shù)學(xué)模型

1）連續(xù)性方程［15］

2）能量方程

3）運動方程

式中，d為定性尺寸，m；u為流體流速，m/s；v為流體的運動粘度，m2/s。Re＜2300，屬于層流，符合達(dá)西定律。

1.2 物理模型

豎直地埋管換熱器與周圍土壤的傳熱過程極其復(fù)雜，模擬過程中計算量太大，需要對模型中影響換熱狀況較小的因素進(jìn)行簡化。故對本文中所建立的模型做出如下假設(shè)：

（1）實際室外溫度是隨時間變化的，為保證地表溫度不受其周期變化的影響，忽略地表空氣溫度對地下溫度的影響；

（2）假設(shè)土壤為各向同性、均質(zhì)的多孔介質(zhì)且回填、埋管、土壤在換熱過程中物性參數(shù)固定不變；

（3）假設(shè)含地下水滲流區(qū)域滲流方向只沿水平方向；

（4）假設(shè)管壁與回填、回填與土壤之間充分接觸。

（5）為方便確定模擬時滲流方向，依據(jù)實際工程項目對模型整體簡化成令土壤橫截面為對角線6m的正方形，高為實際埋深55m，整體呈長方體規(guī)則結(jié)構(gòu)。

（6）建立多孔介質(zhì)區(qū)域模型時，并不是實際建立帶有孔隙的幾何體，而是在模擬求解時定義孔隙率反映出該區(qū)域為多孔介質(zhì)區(qū)域。

1.3 幾何模型及網(wǎng)格劃分

圖1 幾何模型

幾何模型如圖1所示，鉆孔半徑為100mm，井深55m，鉆孔處與土壤模型外邊界最短距離約2.12m，鉆孔內(nèi)U型管外徑32mm、支管中心間距70mm。

利用模擬軟件畫出網(wǎng)格總數(shù)達(dá)到130多萬,考慮重力作用，重力方向沿z軸負(fù)方向。由于彎管處流場梯度變化很大，故采用非結(jié)構(gòu)網(wǎng)格，兩支管及回填和土壤模型部分使用結(jié)構(gòu)性網(wǎng)格。網(wǎng)格扭曲率0.7的占100%，在0.95以內(nèi)，符合模擬標(biāo)準(zhǔn)。

1.4 初始條件及邊界條件

（1）初始條件

認(rèn)為系統(tǒng)初始時達(dá)到熱平衡，土壤、回填、地埋管、循環(huán)流體的溫度都為土壤初始溫度。

（2）邊界條件

埋管入口采用定溫度、定流速邊界條件，鉆孔內(nèi)埋管進(jìn)出口邊界設(shè)置為velocity inlet和pressure outlet，進(jìn)出口位置沿x軸正方向布置，管壁及回填與土壤界面處thermal conditios設(shè)置為coupled［16］。無限遠(yuǎn)處設(shè)置為絕熱壁面temperature。

（3）土壤溫度分析

地面以下土壤溫度分布得計算，假想巖土體為無限大均質(zhì)物體，根據(jù)傅里葉定律求得地下巖土溫度計算公式為［17］：

利用函數(shù)繪圖軟件繪制出如圖大連地區(qū)不同深度土壤溫度分布，得出土壤初始溫度為284.3k。

圖2 大連全年土壤溫度分布

2 評價指標(biāo)與模擬驗證

2.1 模擬結(jié)果評價參數(shù)的定義

（1）地埋管換熱效率

地埋管換熱效率為φ，關(guān)系式如下：

實際上換熱器換熱效率的大小最終取決于埋管進(jìn)出口溫差，即換熱效率為進(jìn)出口溫差與進(jìn)口與地面初始溫度之差的比的絕對值，當(dāng)系統(tǒng)工作滿足建筑負(fù)荷要求時，盡可能將φ利用到最大。從理論的角度出發(fā)，實際換熱量和理想中最大換熱量相等，即φ=1；若埋管換熱無效，則φ=0；若實際運行時當(dāng)?shù)芈窆軗Q熱效率φ低于0.05時就近似認(rèn)為其換熱失效。

（2）單位井深換熱量

式中：ql為 單井換熱量，W/m；A—地埋管橫截面面積，m2；u—流體流速，m/s；

2.2 地埋管模型可靠性驗證

對大連某一以土壤源熱泵系統(tǒng)為采暖方式的獨棟別墅進(jìn)行研究，并嚴(yán)格按照土壤源熱泵施工規(guī)范［18］要求進(jìn)行熱響應(yīng)實驗。在實際施工過程中發(fā)現(xiàn)地下40m左右有滲流的情況，因此對模型進(jìn)行分層設(shè)置，40m以上為無滲流區(qū)，40m以下為有滲流區(qū)且流動狀態(tài)為層流。interval設(shè)置為6s，number of iterations為300，autosave every為200，材料均采用自定義模式，在區(qū)域設(shè)置中勾選Laminar zone和Porous zone，孔隙率為0.385。埋管入口流體溫度277K，流速0.68m·s-1。利用數(shù)值模擬軟件進(jìn)行計算，可以得出埋管出口溫度的數(shù)值解，通過將模型計算出的數(shù)值解與自行研制的TRT裝置進(jìn)行熱響應(yīng)實驗［19］現(xiàn)場實測得的埋管出口溫度值進(jìn)行對比，來驗證所建模型的可靠性。

在模擬過程中所涉及到的各部分物性參數(shù)均由實際工程中熱響應(yīng)實驗測得，具體數(shù)據(jù)見表1。

表1 豎直埋管換熱器模型各項參數(shù)

將地埋管換熱系統(tǒng)運行一周的現(xiàn)場實測出口溫度與本文所建模型得出的埋管出口溫度繪制成如圖3所示的折線圖，可明顯看出埋管出口溫度隨時間的變化情況。由本文所建立的地埋管三維傳熱模型計算出的埋管出口溫度與實驗值偏差1.03%，故該模型符合實際運行工況，可以進(jìn)行下一步分析研究。

圖3 地埋管出口溫度對比圖

3 模擬結(jié)果及分析

3.1 有無滲流時埋管出口溫度

模擬中對無滲流、滲流流速為30m/a、200m/a分別進(jìn)行研究對比。從圖4可以看出系統(tǒng)初運行時出口溫度急劇下降，這是由于模擬時認(rèn)為系統(tǒng)運行前已處于熱平衡狀態(tài)，故管內(nèi)流體初始溫度為土壤初始溫度，因管內(nèi)流體在系統(tǒng)運行時處于流動狀態(tài)，埋管入口低溫流體不斷與管內(nèi)流體混合換熱，從而導(dǎo)致埋管出口溫度急劇下降。

從圖中還可以看出有滲流比無滲流時埋管出口溫度高，滲流流速越大地埋管出口溫度越高。這是由于水的熱容比土壤大很多且地下水處于流動狀態(tài)進(jìn)而促進(jìn)其換熱。冬季埋管進(jìn)出口溫差依次為2.9k、3.47k、3.9k。

圖4 冬季換熱器出口溫度

3.2 有無滲流時換熱性能比較

從圖5中看出有滲流比無滲流的換熱效率高，且滲流速度越大地埋管換熱效率越高。比較發(fā)現(xiàn)冬季單井換熱量、換熱效率分別為33.5W·m-1、φ=0.4，39.7W·m-1、φ=0.48，44.7W·m-1、φ=0.53。

圖5 冬季換熱器換熱性能

3.3 鉆孔側(cè)不同位置土壤溫度波動

從圖6中看出，滲流速度為30m·a-1時，靠近鉆孔側(cè)的土壤溫度下降緩慢，距地下水滲流下游鉆孔側(cè)一米處的土壤溫度波動影響可近似忽略，土壤平均溫度波動不明顯。流速為200m·a-1時，靠近鉆孔的土壤溫度在3h內(nèi)迅速下降，對距地下水滲流下游鉆孔側(cè)一米處的土壤溫度在48h后有明顯下降趨勢，土壤平均溫度呈現(xiàn)出稍許下降趨勢。這是由于滲流速度大時埋管換熱器與土壤換熱增強(qiáng)，土壤對埋管傳遞的熱量增加，進(jìn)而靠近鉆孔側(cè)土壤溫度迅速下降。滲流的流速越大，由地下水帶入到下游鉆孔側(cè)一米處的冷量越多，使得下游鉆孔側(cè)一米處土壤溫度下降得越快。

整體來看，地下水滲流流速快，鉆孔周圍土壤溫度下降趨勢越明顯，達(dá)到平衡態(tài)溫度的時間越短。但不管滲流流速大還是小，運行一段時間以后土壤任一處溫度都將趨于平穩(wěn)。

圖6 鉆孔側(cè)土壤溫度

3.4 有無滲流時土壤溫度分布云圖

從圖7可以看出，無滲流時，鉆孔周圍土壤溫度波動以埋管為圓心呈圓形對稱分布，溫度梯度由內(nèi)向外均勻減小。在滲流流速為30m·a-1時，鉆孔周圍土壤溫度波動沿著滲流方向偏移。在滲流流速為200m·a-1時，鉆孔周圍土壤溫度呈上下對稱的偏橢圓形分布，且沿著滲流方向偏移量更大。這是由于無滲流時土壤與換熱器之間換熱過程是純導(dǎo)熱，溫度均勻變化。有滲流時地下水處于流動狀態(tài)，勢必會帶走土壤中的熱量，進(jìn)而影響土壤溫度波動范圍。

有滲流與無滲流相比時，土壤溫度波動范圍有較大的不同，同時在不同滲流速度下土壤溫度波動也大不相同，沿著滲流方向影響范圍變大，垂直滲流方向影響范圍變小，且鉆孔上游側(cè)土壤溫度比下游高，故埋管換熱器進(jìn)口在上游時換熱效果好，有利于提高埋管換熱器的換熱性能。

圖7 有無滲流作用下的溫度云圖

4 結(jié)論

（1）地埋管換熱器在有地下水滲流的土壤中換熱時，可提高其換熱效率。在200m·a-1的滲流速度下，冬季埋管進(jìn)出口溫差3.9k，單位鉆井深度取熱量44.7W·m-1。與在無滲流土壤中換熱相比，取熱量增加約33.4%，冬季φ提升32.5%，能夠大幅度減少系統(tǒng)的初投資。

（2）有滲流時，鉆孔側(cè)上游土壤溫度比下游高，地下水將地埋管換熱產(chǎn)生的冷量帶到了下游，這使得土壤冷量不容易集中在埋管測，解決了由于地埋管長期運行造成的土壤冷堆積問題。