山東建筑大學(xué) 關(guān)春敏 山東中瑞新能源科技有限公司 趙樹旺 山東建筑大學(xué) 張文科 孫炳巖 王建華 孫知曉 山東中瑞新能源科技有限公司 方肇洪

0 引言

在能源緊張及氣候變暖、霧霾天氣頻頻出現(xiàn)等環(huán)境問題的壓力下，熱泵清潔供暖技術(shù)成為業(yè)界關(guān)注和研究的熱點[1]。地源熱泵技術(shù)作為一種新能源技術(shù)，具有環(huán)保節(jié)能、能效比高等優(yōu)勢[2]。

淺層地源熱泵技術(shù)主要利用地下深度范圍在50～200 m的淺層地?zé)崮埽枰紦?jù)一定的土地面積布置地埋管，而且在冷熱負(fù)荷不平衡地區(qū)或單供暖地區(qū)，地下介質(zhì)的溫度逐年變化會導(dǎo)致地下傳熱性能越來越差。

鑒于淺層地?zé)崮艿卦礋岜眉夹g(shù)的缺點，為提高地?zé)崮艿睦眯?，國?nèi)外的專家學(xué)者近些年開展了對中深層地?zé)崮艿难芯?，把鉆孔埋管的深度加大到1 500～3 000 m，稱之為中深層地埋管換熱器技術(shù)。一個鉆孔地埋管換熱器可滿足幾千甚至上萬m2的建筑物供暖，不僅減少了占地面積，而且可利用的土壤溫度顯著提高，地下溫度基本穩(wěn)定，適合于單供暖的情況。相對于水熱型的中深層地?zé)嵯到y(tǒng)，中深層地埋管換熱器技術(shù)采用閉式埋管循環(huán)系統(tǒng)，故對地質(zhì)條件的限制比較少，可以靈活應(yīng)用于多種地質(zhì)條件。避免了開式系統(tǒng)與地下水或巖層的直接接觸或質(zhì)傳遞，從而避免了對地下水生物化學(xué)方面的影響。

在從地下1 500～3 000 m深、溫度在70～90 ℃甚至更高范圍的巖石中提取地?zé)崮茏鳛闊岜孟到y(tǒng)低溫?zé)嵩吹倪^程中，為了保證地下?lián)Q熱裝置的結(jié)構(gòu)穩(wěn)定性，通常采用套管結(jié)構(gòu)，即循環(huán)水在內(nèi)外管嵌套的地埋管內(nèi)流動，實現(xiàn)與地下介質(zhì)的換熱。外管通常采用鋼管起固定及強化傳熱的作用，而內(nèi)管則用塑料管以對循環(huán)水進行一定程度的保溫。從提高換熱效果的角度考慮，循環(huán)水在水泵的驅(qū)動下從內(nèi)外管的空隙向下流動，到達豎直管的底部后，再由內(nèi)管向上流出換熱裝置，地埋管換熱器的結(jié)構(gòu)及循環(huán)水流動方式如圖1所示。

圖1 中深層套管式地埋管換熱器示意圖

中深層地源熱泵作為一種對地?zé)崮芾幂^為高效的方式，受到業(yè)內(nèi)研究人士的重視且已逐漸推廣應(yīng)用。但是，由于中深孔地源熱泵打井和埋設(shè)管道的成本相對較高且一旦施工結(jié)束后無法進行維修，所以必須對工程項目進行分析和預(yù)測，探索不同因素對套管換熱器取熱能力的影響，研究循環(huán)水在系統(tǒng)運行過程中的溫度變化，以及揭示單取熱工況對地下土壤介質(zhì)的影響。本文根據(jù)中深層套管式地埋管換熱器的傳熱模型，以實際供暖工程項目為研究對象，對地埋管換熱器的運行工況進行傳熱分析，為實際的技術(shù)應(yīng)用提供理論指導(dǎo)，進一步促進中深層地?zé)崮艿耐茝V應(yīng)用。

1 地?zé)峋こ谈艣r

該工程位于河北省邯鄲市，地處寒冷地區(qū)，當(dāng)?shù)卮髿饽昶骄鶞囟葹?3 ℃，地表對流換熱系數(shù)為15 W/(m2·K)。單井的地下介質(zhì)包括5個地層，各個地層介質(zhì)的信息見表1。

表1 邯鄲某工程中深層鉆井地層信息

該工程鉆井深度為2 500 m，鉆井口徑為?311 mm，使用閉式套管式換熱器，以水為循環(huán)介質(zhì)，套管的內(nèi)外管均為鋼管。其中，內(nèi)管經(jīng)過保溫處理后的導(dǎo)熱系數(shù)為0.02 W/(m·K)，單位體積熱容為3.627 MJ/(m3·K)；外管的導(dǎo)熱系數(shù)為45 W/(m·K)，單位體積熱容為3.634 MJ/(m3·K)。在0～420 m處設(shè)置隔熱段，隔熱段使用導(dǎo)熱系數(shù)為1.48 W/(m·K)的普通水泥回填，而420～2 445 m處使用導(dǎo)熱系數(shù)為2.68 W/(m·K)的普通水泥回填，經(jīng)過計算，回填材料的平均導(dǎo)熱系數(shù)為2.476 W/(m·K)，平均單位體積熱容為2.632 MJ/(m3·K)。

2 傳熱模型

求解套管式地埋管換熱器的傳熱模型需考慮換熱器及周圍巖土地層的實際特點。由于鉆孔埋管的深度很大，整個鉆孔可穿越若干具有不同介質(zhì)的地層，每個地層介質(zhì)的參數(shù)應(yīng)分別統(tǒng)計，鉆孔底部的溫度可能比地表溫度高50 ℃或者更高，軸向地溫梯度不可忽略。由于地下介質(zhì)的溫度沿鉆孔深度方向逐漸變化，熱量自地下至地面以導(dǎo)熱的方式進行傳遞，成為傳熱分析中必須注意的一個因素。循環(huán)液的流動方式已在圖1中說明，向下和向上的流動要區(qū)別對待。此外，整個傳熱研究較為復(fù)雜，為便于探索，可對一些情形進行簡化。因此，可總結(jié)建立傳熱模型的前提條件，并對循環(huán)液和巖土層分別建立能量控制方程，進而對傳熱模型求解分析。

2.1 前提條件

建立中深層地埋管換熱器的傳熱模型，需要基于一定的前提條件，對其進行總結(jié)如下：

1) 在所考察的范圍內(nèi)將巖土層看作是均勻的介質(zhì)[3]，或分層均勻的介質(zhì)，介質(zhì)的物性不隨溫度改變，傳熱機理僅考慮導(dǎo)熱。

2) 不考慮地下水流動的影響[4]，將單一套管式地埋管的換熱問題簡化成二維非穩(wěn)態(tài)的導(dǎo)熱問題。

3) 套管埋管內(nèi)外管中流體的流動和傳熱都按一維問題考慮，即不考慮管內(nèi)流體在橫截面上的速度和溫度分布，而僅考慮在橫截面上的平均溫度和平均流速。

4) 根據(jù)所考慮的時間尺度，在徑向限定一個足夠大而有限的范圍，認(rèn)為在該徑向邊界上的溫度保持原來的溫度而沒有受到套管傳熱的影響。

5) 在所考察的范圍內(nèi)有均勻向上的大地?zé)崃?。大地?zé)崃魇侵赣傻匦南蛲鈧鬟f的能量，即單位時間內(nèi)流經(jīng)地表單位面積的熱量[5],其計算公式如下：

(1)

式中q為大地?zé)崃?，W/m2；λd為地層的導(dǎo)熱系數(shù)，W/(m·K)；t為地溫，℃；z為軸向距離，m。

6) 初始時刻巖土在徑向的溫度分布是均勻的，在軸向則存在溫度梯度。因為大地?zé)崃鞅徽J(rèn)為是均勻的，所以不同水平地層中的地溫梯度也不同。

7) 地表以上的空氣與地面進行對流換熱，空氣的溫度及表面對流換熱系數(shù)保持不變。

2.2 控制方程

2.2.1套管內(nèi)外管流體溫度的控制方程

來自土壤的熱量在徑向傳遞過程中受到的阻礙包括外管內(nèi)循環(huán)液和鉆孔壁之間單位長度的熱阻R1和內(nèi)管流體與外管流體之間單位長度的熱阻R2，其計算公式如下:

(2)

(3)

式(2)、(3)中d2i和d2o分別為外管的內(nèi)徑與外徑，m；h2和h1分別為套管外管與內(nèi)管的對流換熱系數(shù)，W/(m2·K)；λp2、λg、λp1分別為套管外管、回填漿料及套管內(nèi)管的導(dǎo)熱系數(shù)，W/(m·K)；db為鉆孔直徑，m；d1i、d1o分別為內(nèi)管的內(nèi)徑與外徑，m。

套管內(nèi)流體的流動方式為外進內(nèi)出，內(nèi)管和外管中流體流動的能量方程式分別為

(4)

(5)

式(4)、(5)中C1、C2分別為內(nèi)、外管道熱流通道單位長度的熱容量，J/(m·K)；C為流體的熱容量，J/(s·K)；tf1、tf2分別為內(nèi)、外管中流體的溫度，℃；τ為時間，s；tb為鉆孔壁的溫度，℃。

C1、C2、C的計算式如下：

(6)

(7)

C=mc

(8)

式(6)、(7)中ρ、ρ1、ρ2、ρg分別為水、內(nèi)管、外管和回填漿料的密度，kg/m3；c、c1、c2、cg分別為水、內(nèi)管、外管和回填漿料的比熱容，J/(kg·K)；m為套管內(nèi)水的質(zhì)量流量，kg/s。

2.2.2巖土中的能量控制方程

每個巖土層的導(dǎo)熱方程如下：

(9)

式中a為地層的熱擴散率，m2/s；r為徑向坐標(biāo)，m。

由于柱坐標(biāo)中徑向溫度梯度分布受地埋管換熱器的影響，分布不均勻，其中心部分的溫度梯度大而邊緣部分的溫度梯度小，在徑向可以采用不均等的差分步長。為此引進新坐標(biāo)σ：

(10)

式中r0為鉆孔的半徑，m。

對式(9)進行坐標(biāo)變換得:

(11)

如果σ坐標(biāo)采用等步長Δσ，則對應(yīng)的r坐標(biāo)就成為一個等比級數(shù):

(12)

式中ri+1為第i+1個徑向節(jié)點的徑向坐標(biāo)，m；ri為第i個徑向節(jié)點的徑向坐標(biāo)，m；r1為第1個徑向節(jié)點的徑向坐標(biāo)，m。

計算中，取β=1.2，Δσ=0.182 32。

3 項目的傳熱分析

3.1 名義取熱量的研究

3.1.1名義取熱量的定義

名義取熱量是一個通過模擬來確定的虛擬的量化指標(biāo)，便于工程師與非專業(yè)人士進行溝通。其定義為在特定運行工況下，一個鉆孔換熱器可以提供的最大取熱量。這種特定工況指的是[6-8]：

1) 取熱量在3個月內(nèi)(90 d)是恒定的。

2) 巖土層的初始溫度分布值是預(yù)先確定好的，在這里認(rèn)為取熱開始時巖土層未受中深層地埋管換熱器系統(tǒng)的擾動。

3) 中深層地埋管換熱器系統(tǒng)的進口溫度在整個取熱期間不得低于5 ℃。

3.1.2鉆井深度對名義取熱量的影響

為了研究鉆井深度對名義取熱量的影響，在該工程的基礎(chǔ)上，其他條件不變只改變鉆井深度，分別對鉆井深度為1 500、1 800、2 100、2 400、2 700、3 000 m的情況進行了模擬研究，得到了名義取熱量隨鉆井深度的變化，如圖2所示。

圖2 名義取熱量隨鉆井深度的變化

從圖2可以看出，隨著鉆井深度的增大，名義取熱量也明顯增大，這是因為鉆井越深，越接近地心，受到地心溫度的影響就越大，鉆孔底部的溫度也越高，地埋管換熱器可以通過與土壤的換熱獲得更多的熱量。但是，鉆孔越深也就意味著工程成本越高，因此應(yīng)以此工程模擬為依據(jù)并綜合考慮經(jīng)濟效益與取熱量，選擇最優(yōu)的鉆孔深度以實現(xiàn)效益最大化。

3.1.3循環(huán)液流量不同時的名義取熱量

根據(jù)Q=cmΔt(其中Q為換熱量；Δt為溫差)可知，循環(huán)液流量越大，地埋管換熱器取得的熱量越多。為了對此進行驗證，分別對循環(huán)液流量為20、30、40、50、60、70、80 m3/h時的運行工況進行模擬，得到循環(huán)液流量對鉆孔地埋管換熱器名義取熱量的影響，結(jié)果如圖3所示。

圖3 名義取熱量隨循環(huán)液流量的變化

由圖3可以看出，隨著地埋管換熱器中循環(huán)液流量的增大，名義取熱量也呈現(xiàn)增大的趨勢，但是增加幅度越來越小。這是由于當(dāng)循環(huán)液流量較大時，地埋管換熱器持續(xù)從土壤中吸收較多的熱量，使得土壤溫度下降后不能較快地恢復(fù)。如果要維持取熱量的穩(wěn)定，則循環(huán)水進出口溫度會降低，甚至?xí)陀诮Y(jié)冰點，從而影響整個系統(tǒng)的正常運行。同時，循環(huán)液流量增大時，對循環(huán)水泵的要求也隨之提高。因此，在建設(shè)實際工程時應(yīng)選取合適的循環(huán)液流量，以達到更好的運行效果，建議以35～50 m3/h為宜。

3.1.4大地?zé)崃鞑煌瑫r的名義取熱量

大地?zé)崃魇悄軌蚍从惩寥罒釥顟B(tài)的物理量[9]，它反映了在單位時間、單位面積內(nèi)由地心向地表傳輸能量的多少。本文模擬了大地?zé)崃鞣謩e為0.04、0.05、0.06、0.07、0.08 W/m2時名義取熱量的變化，結(jié)果如圖4所示。

圖4 名義取熱量隨大地?zé)崃鞯淖兓?/p>

從圖4可以看出，隨著大地?zé)崃鞯脑龃?，地埋管換熱器的名義取熱量呈現(xiàn)明顯的增大趨勢。這說明當(dāng)?shù)叵陆橘|(zhì)的導(dǎo)熱系數(shù)不發(fā)生變化時，大地?zé)崃髟酱笠馕吨叵陆橘|(zhì)的溫度梯度越大，地埋管換熱器可以從周圍土壤中獲得更多的熱量，越有利于地埋管換熱器的運行。

3.2 模擬數(shù)據(jù)與實驗數(shù)據(jù)的對比

在本次模擬中，除了對名義取熱量的研究外，其余的研究均在額定取熱量490 kW下進行，為了探究該取熱量的選取是否合理，將通過實驗獲得的取熱量數(shù)據(jù)與其進行對比，實驗運行時的循環(huán)液流量為45 m3/h，該工況的實驗時間為2019年7月31日09：00至8月10日09：00，共240 h，其中實驗在進行過程中是每隔一段時間讀取一個數(shù)據(jù)，對比結(jié)果如圖5所示。

圖5 模擬數(shù)據(jù)與實驗數(shù)據(jù)的比較

從圖5可以看出，實驗數(shù)據(jù)與選定的模擬數(shù)據(jù)中的取熱量有一定的誤差，其原因可能是由于地下地質(zhì)結(jié)構(gòu)較為復(fù)雜，可能存在地下水滲流的影響，使得地埋管換熱器在真實的運行工況下無法保證其取熱量為一個恒定的數(shù)值，但實驗數(shù)據(jù)仍然圍繞模擬選定的取熱量數(shù)據(jù)波動，存在的誤差較小，可以忽略。因此，在接下來的模擬中選擇取熱量為490 kW是合理的。

3.3 進出口溫度隨時間的變化

循環(huán)水進出口溫度的變化也是衡量土壤受地埋管換熱器影響的一個標(biāo)準(zhǔn)，此處通過模擬考察了系統(tǒng)在連續(xù)運行工況下持續(xù)運行5 a進出口溫度的變化，如圖6所示。

圖6 系統(tǒng)連續(xù)運行5 a進出口溫度隨時間的變化

由圖6可以看出，循環(huán)液的進出口溫度逐年緩慢下降，但在第5年仍然能保持在60 ℃左右，說明雖然地埋管換熱器的取熱會對土壤造成一些影響，但是土壤溫度在非供暖時間段可以得到較好的恢復(fù)，以維持其供熱能力。

3.4 換熱器對土壤的熱影響半徑

熱影響半徑定義為：在鉆孔深度范圍內(nèi)，某深度處某個半徑坐標(biāo)的土壤溫度與該土層土壤初始溫度之差的絕對值小于0.5 ℃時，則認(rèn)為這個半徑為該深度的熱影響半徑[10]。

通過模擬得出地埋管換熱器運行10 a后鉆孔500、1 000、1 500、2 000、2 500 m深處土壤的徑向溫度分布，如圖7所示。由圖7可以看出：隨著徑向距離的增大，土壤溫度先逐漸升高后基本保持不變；徑向距離相同時，隨著深度增大，土壤溫度逐漸升高。

圖7 不同深度土壤溫度沿徑向的變化

將地埋管換熱器運行第1 h的土壤溫度分布作為土壤的初始溫度，然后用已經(jīng)得出的10 a后土壤在不同深度處的徑向溫度分布與其作差，并找出絕對值小于0.5 ℃時的徑向距離，即熱泵運行到該時刻在該深度處的熱影響半徑，將其統(tǒng)計列入表2中。

表2 不同深度處的熱影響半徑

由表2可以看出，隨著深度的增加，地埋管換熱器的熱影響半徑呈現(xiàn)先增大后減小的趨勢，這可能是由于當(dāng)深度增加到某一數(shù)值時，地埋管換熱器中的循環(huán)液與地下介質(zhì)之間的換熱密度增大，土壤溫度受影響的范圍也會變大；而在鉆孔的底部，由于此時循環(huán)液已經(jīng)到達了地埋管換熱器的邊界處，與土壤之間的換熱程度減弱，對土壤中溫度的影響范圍也隨之減小。同時，可以看出，該地埋管換熱器在運行10 a后的最大熱影響半徑為44.14 m，因此，為了防止熱干擾現(xiàn)象的發(fā)生，在距離該地埋管換熱器中心88.28 m的區(qū)域內(nèi)不宜設(shè)置其他地埋管換熱器。

3.5 鉆孔壁溫度沿深度方向的變化

在整個換熱過程中，鉆孔壁的溫度受土壤溫度與地埋管換熱器內(nèi)循環(huán)液溫度的雙重影響，圖8顯示了地埋管換熱器運行的第720 h(1個月)、1 440 h(2個月)、2 160 h(3個月)、2 880 h(4個月)鉆孔壁溫度分布情況。

圖8 不同時刻鉆孔壁溫度分布

由圖8可以看出，鉆孔壁溫度先升高，在480 m左右處出現(xiàn)一個小幅波動，然后繼續(xù)升高。這是由于在距地面0～480 m左右的深度處土壤溫度低于循環(huán)液進口溫度，鉆孔壁會吸收循環(huán)液中的熱量，而在480 m以下的位置，土壤溫度高于循環(huán)液溫度，循環(huán)液開始從鉆孔壁處吸收熱量，同時鉆孔壁也會吸收來自土壤的熱量，在二者的雙重影響下，鉆孔壁溫度繼續(xù)升高。同時，隨著地埋管換熱器運行時間的增長，鉆孔壁溫度呈現(xiàn)降低的趨勢，這說明，地埋管換熱器的運行會使土壤溫度有一定程度的降低。

4 結(jié)論

1) 地埋管換熱器的名義取熱量隨著鉆井深度、循環(huán)液流量及大地?zé)崃鞯脑龃蠖饾u增大。鉆井深度增大意味著工程項目的初投資增加，應(yīng)結(jié)合經(jīng)濟成本考慮鉆井深度。循環(huán)液流量的增加對水泵的要求提高，且名義取熱量隨流量的增大幅度越來越小，應(yīng)根據(jù)不同項目設(shè)定合理的流量值。大地?zé)崃鞯脑龃罂梢允姑x取熱量呈現(xiàn)明顯的增大趨勢，說明當(dāng)?shù)叵陆橘|(zhì)的導(dǎo)熱系數(shù)不發(fā)生變化時，大地?zé)崃髟酱螅叵陆橘|(zhì)的溫度梯度越大，地埋管換熱器可以從周圍土壤介質(zhì)獲得更多的熱量，有利于地埋管換熱器的運行。

2) 隨著運行時間的增加，循環(huán)液進出口溫度逐年緩慢下降，但在第5年仍然可以維持在60 ℃左右，這說明地埋管換熱器仍然能夠維持良好的取熱能力。

3) 地埋管換熱器運行10 a后，隨著深度的增加，熱影響半徑先增大后減小，這可能是由于隨著深度的增大，地埋管換熱器中的循環(huán)液與地下介質(zhì)之間的換熱密度增大，土壤溫度受影響的范圍也隨之增大；而在鉆孔的底部，由于循環(huán)液已經(jīng)到達了地埋管換熱器的邊界處，與土壤介質(zhì)之間的換熱程度減弱，對土壤溫度的影響范圍也隨之減小。

4) 地埋管換熱器鉆孔壁的溫度隨著深度的增大而升高，隨著時間的增長而降低，這是因為當(dāng)深度越大時，土壤溫度越高；而隨著地埋管換熱器運行時間的增長，周圍土壤的溫度也會受到一定程度的影響而有所降低。