中深層地埋管換熱器取熱穩(wěn)定性及熱影響半徑

黃 帥，朱 科，董建鍇，李 驥，姜益強(qiáng)，方肇洪

(1.哈爾濱工業(yè)大學(xué) 建筑學(xué)院，哈爾濱 150006；2.寒地城鄉(xiāng)人居環(huán)境科學(xué)與技術(shù)工業(yè)和信息化部重點(diǎn)實(shí)驗(yàn)室(哈爾濱工業(yè)大學(xué))，哈爾濱 150090；3.山東建筑大學(xué) 熱能工程學(xué)院，濟(jì)南 250101；4.中國(guó)建筑科學(xué)研究院有限公司，北京 100101)

地?zé)崮苜Y源豐富，具有分布廣、清潔環(huán)保、儲(chǔ)量大、穩(wěn)定可靠等特點(diǎn)[1]。在建筑行業(yè)地?zé)崮荛_(kāi)發(fā)利用過(guò)程中，淺層地源熱泵系統(tǒng)應(yīng)用較為廣泛[2]，但其對(duì)空調(diào)全年冷熱負(fù)荷要求較高，長(zhǎng)期運(yùn)行后系統(tǒng)性能會(huì)顯著下降[3]，并且淺層地埋管換熱器埋管占地面積較大。近年來(lái)，中深層地源熱泵系統(tǒng)作為新興的地?zé)崮芾眯问剑ǔＦ渎窆苌疃仍? 000～3 000 m，具有地埋管換熱器占地面積小、對(duì)全年冷熱負(fù)荷平衡要求低、系統(tǒng)取熱性能穩(wěn)定等優(yōu)勢(shì)[4]，受到了廣泛關(guān)注。

與傳統(tǒng)淺層地埋管換熱器相比，中深層地埋管供熱技術(shù)只在近年來(lái)得以迅速發(fā)展[5]，并且其發(fā)展現(xiàn)狀是工程實(shí)踐先于理論研究，有很多問(wèn)題尚待于研究[6]。其中，在中深層地埋管換熱器取熱特性方面，Li等[7]對(duì)影響中深層地埋管換熱器取熱特性的因素進(jìn)行分析，研究結(jié)果表明增大外管管徑和流速有利于提升中深層地埋管換熱器的取熱量；Liu等[8]對(duì)中深層地埋管換熱器管內(nèi)循環(huán)工質(zhì)的流速進(jìn)行了優(yōu)化，并提出最佳流速為0.7 m/s；Fang等[9]發(fā)現(xiàn)增大中深層地埋管換熱器的內(nèi)套管熱阻能夠提升其換熱能力。而對(duì)中深層地埋管換熱器周?chē)膸r土熱恢復(fù)特性等問(wèn)題研究較少，其中，賈林瑞等[10]研究了單取熱工況地下巖土溫度熱響應(yīng)，發(fā)現(xiàn)1 400～1 800 m處的巖土溫度變化率最快；Cai等[11]對(duì)間歇運(yùn)行模式下中深層地埋管換熱器取熱穩(wěn)定性進(jìn)行了研究，發(fā)現(xiàn)不同運(yùn)停比運(yùn)行時(shí)的巖土溫度熱恢復(fù)效果較好；Song等[12]研究結(jié)果表明巖土經(jīng)過(guò)為期8個(gè)月的熱恢復(fù)期，其溫度基本能恢復(fù)至初始狀態(tài)。但上述研究主要是針對(duì)單個(gè)供熱周期，其目的是量化設(shè)計(jì)參數(shù)對(duì)中深層套管式地埋管換熱器取熱性能的影響以及揭示巖土的熱恢復(fù)規(guī)律。考慮到中深層地?zé)岬芈窆芄嵯到y(tǒng)需長(zhǎng)期運(yùn)行，因此，為明確利用該系統(tǒng)供熱的可持續(xù)性，需要對(duì)中深層套管式地埋管換熱器長(zhǎng)期取熱的穩(wěn)定性進(jìn)行相關(guān)研究。

此外，中深層地埋管換熱器取熱對(duì)周?chē)鷰r土的熱擾與諸多參數(shù)相關(guān)，如運(yùn)行時(shí)間、管徑、循環(huán)流量、回填材料導(dǎo)熱系數(shù)及巖土熱物性等。在工程實(shí)踐過(guò)程中，為解決中深層地埋管換熱器鉆孔間距不當(dāng)造成的熱短路以及土地資源浪費(fèi)等問(wèn)題，首先要確定中深層地埋管換熱器取熱對(duì)周?chē)鷰r土的熱擾與上述哪些參數(shù)相關(guān)性較大，因此，需要對(duì)上述參數(shù)做敏感性分析。

基于此，本文建立了中深層地埋管換熱器數(shù)值傳熱模型，基于有限差分法將控制方程離散求解，并利用項(xiàng)目實(shí)測(cè)數(shù)據(jù)進(jìn)行驗(yàn)證，在此基礎(chǔ)上研究長(zhǎng)期運(yùn)行下中深層套管式地埋管換熱器的出水溫度以及巖土的溫度變化。此外，對(duì)影響中深層套管式地埋管換熱器熱擾的運(yùn)行參數(shù)進(jìn)行了敏感性分析，探尋中深層地埋管換熱器的熱影響半徑與相關(guān)參數(shù)的表達(dá)式，為日后的工程實(shí)踐提供理論指導(dǎo)。

1 數(shù)值模型建立

1.1 物理模型

中深層同軸套管式地埋管換熱器的取熱原理如圖1所示，在循環(huán)水泵的驅(qū)動(dòng)下循環(huán)工質(zhì)(水)由套管外部環(huán)腔注入，經(jīng)地埋管換熱器周?chē)邷貛r土層加熱后通過(guò)內(nèi)套管反向流動(dòng)，流出地面后進(jìn)入熱泵機(jī)組。

圖1 中深層地埋管換熱器的換熱原理Fig.1 Heat exchange principle of MDBHE

中深層地埋管換熱器的取熱過(guò)程主要包括內(nèi)套管流體與外部環(huán)腔流體的對(duì)流換熱、環(huán)腔流體與外管壁的對(duì)流換熱以及外套管與周?chē)鷰r土的導(dǎo)熱，其換熱過(guò)程較為復(fù)雜，并且地下巖土環(huán)境具有不確定性，因此做出以下假設(shè)：

1)忽略地下滲流作用的影響，將巖土中傳熱視為單純的導(dǎo)熱問(wèn)題[13]；

2)忽略地表溫度波動(dòng)的影響[14]；

3)假定數(shù)值模擬區(qū)域的徑向邊界處的溫度分布不受地埋管換熱器的影響；

4)假定同軸套管式換熱器內(nèi)循環(huán)體主要以對(duì)流傳熱，忽略其軸向的熱傳導(dǎo)。

1.2 控制方程

基于以上假設(shè)，每層巖土的導(dǎo)熱方程可寫(xiě)為[15]

(1)

式中：ak為不同巖土層的熱擴(kuò)散率，m2/s；Ts為巖土層溫度，℃；Δτ為時(shí)間步長(zhǎng)，s；r、z分別為徑向長(zhǎng)度與豎直方向長(zhǎng)度，m。

當(dāng)管內(nèi)流體工質(zhì)以圖1所示的方式流動(dòng)時(shí)，外部環(huán)腔內(nèi)流體的能量方程為

(2)

內(nèi)管流體的能量方程為

(3)

式(2)、(3)描述的是在中深層地埋管換熱器正常取熱階段下的能量方程，當(dāng)系統(tǒng)停運(yùn)時(shí)，外部環(huán)腔內(nèi)流體及內(nèi)管流體的流速為零，此時(shí)地埋管周?chē)鷰r土處于熱恢復(fù)狀態(tài)。與正常取熱階段下的能量方程相比，熱恢復(fù)階段能量方程中的系數(shù)C為零。

其中：

式中：C=MCf為循環(huán)液的熱容流量,kJ/(s·K)；db為鉆孔半徑,m；di、do分別為套管的內(nèi)徑和外徑，m；ρc、ρ1c1、ρ2c2、ρgcg分別為循環(huán)水比定容熱容，外管比定容熱容、內(nèi)管比定容熱容、回填材料比定容熱容，J/(m3·K)。此外，外管內(nèi)循環(huán)水和鉆孔壁之間的熱阻R1以及內(nèi)外管循環(huán)水之間的熱阻R2由下式求得：

(4)

(5)

式中：kg、kp1、kp2分別為回填材料、外管和內(nèi)管的導(dǎo)熱系數(shù)，W/(m·K)；h1、h2分別為套管外管和內(nèi)管的對(duì)流換熱系數(shù)，W/(m2·K)。

1.3 初始及邊界條件

在任一深度的地層的初始地溫可以表達(dá)為下式[9]：

(6)

式中：Hj為第j層地層底部的坐標(biāo)；qg為大地?zé)崃鳎琖/m2；ta為地表溫度，℃；ha為空氣與地表的對(duì)流換熱系數(shù)，W/(m2·K)；km為巖土的導(dǎo)熱系數(shù)，W/(m·K)。

外部環(huán)腔流體與內(nèi)管流體能量方程的邊界體條件為：

(7)

Tf1=Tf2,Z=H

(8)

式中：Q為取熱功率，kW；C=MCf為循環(huán)液的熱容流量,kJ/(s·K)；Z為埋管深度，m。

2 數(shù)值模型求解及驗(yàn)證

2.1 模型求解

在非穩(wěn)態(tài)數(shù)值傳熱模擬中，時(shí)間步長(zhǎng)Δτ、縱向步長(zhǎng)Δz及徑向步長(zhǎng)對(duì)模擬計(jì)算速度及結(jié)果的精確度均有一定影響。為了保證結(jié)果的準(zhǔn)確性需要對(duì)其做無(wú)關(guān)性驗(yàn)證，由圖2可知Δτ可取900 s，此外，考慮徑向較遠(yuǎn)處熱流趨近于零，鉆孔附近的熱流及溫度梯度較大，故以鉆孔壁為界采用變步長(zhǎng)，即徑向步長(zhǎng)按照等比級(jí)數(shù)變化(等比系數(shù)為1.2)，徑向離鉆孔第1個(gè)節(jié)點(diǎn)距離為0.168 m,徑向共取40個(gè)節(jié)點(diǎn)，縱向步長(zhǎng)Δz取10 m。由圖3可知節(jié)點(diǎn)數(shù)量可取7 700個(gè)，滿足計(jì)算精度要求。其余模擬所需的基礎(chǔ)物理參數(shù)見(jiàn)表1，巖土熱物性參數(shù)除特別說(shuō)明外取下文中的案例1。

表1 模擬所需的物理參數(shù)Tab.1 Physical parameters required for simulation

圖2 時(shí)間步長(zhǎng)無(wú)關(guān)性驗(yàn)證Fig.2 Verification of time step independence

圖3 節(jié)點(diǎn)數(shù)量無(wú)關(guān)性驗(yàn)證Fig.3 Verification of node number independence

在中深層地埋管換熱器與周?chē)鷰r土層的耦合傳熱問(wèn)題中，可將其看作圓柱體的二維瞬態(tài)導(dǎo)熱。對(duì)于二維瞬態(tài)耦合傳熱問(wèn)題，若采用隱式格式差分方法得到非邊界處節(jié)點(diǎn)方程中的未知數(shù)通常有5個(gè)，采用迭代法求解耗時(shí)較長(zhǎng)[16]。因此，本文基于微元體熱平衡和交替時(shí)間步長(zhǎng)法[17]建立節(jié)點(diǎn)方程，離散網(wǎng)格如圖4所示。

圖4 模型離散網(wǎng)格示意Fig.4 Schematic of discretization mesh for numerical model

2.2 模型驗(yàn)證

為了驗(yàn)證建立的數(shù)學(xué)模型及所編程序的準(zhǔn)確性，作者所在團(tuán)隊(duì)于2019年11月至2020年1月期間對(duì)西安市某住宅工程進(jìn)行了現(xiàn)場(chǎng)實(shí)測(cè)，測(cè)試參數(shù)主要包括中深層地埋管換熱器進(jìn)出水溫度、機(jī)組及循環(huán)水泵耗電量、用戶側(cè)進(jìn)出水溫度及循環(huán)水泵耗電量，測(cè)試期間機(jī)組的平均COP達(dá)4.7，限于篇幅，其余內(nèi)容詳見(jiàn)文獻(xiàn)[7]。其中，測(cè)試期間單孔地?zé)峋骄鶕Q熱量為270 kW，循環(huán)水平均流量為25.5 m3/h，中深層地埋管換熱器進(jìn)出水溫差約為9.2 ℃，將該地?zé)峋\(yùn)行參數(shù)導(dǎo)入所編制的程序，進(jìn)而得出模擬地源側(cè)進(jìn)出口溫度，選取650 h后的測(cè)試數(shù)據(jù)與模擬值進(jìn)行對(duì)比，地源側(cè)進(jìn)出口溫度的模擬值與實(shí)測(cè)值對(duì)比如圖5所示。從圖5中出水溫度可知，除個(gè)別波動(dòng)點(diǎn)外，最大誤差為6.01%，進(jìn)而驗(yàn)證了數(shù)值傳熱模型及所編程序的正確性，為下文相關(guān)分析奠定了基礎(chǔ)。

圖5 地源側(cè)進(jìn)、出口溫度實(shí)際值與模擬值比對(duì)Fig.5 Comparison of actual and simulated temperature at inlet and outlet of ground source side

3 數(shù)值模擬結(jié)果與分析

考慮到中深層地?zé)岬芈窆芄嵯到y(tǒng)需長(zhǎng)期運(yùn)行，因此，為明確利用該系統(tǒng)供熱的可持續(xù)性，需要對(duì)中深層套管式地埋管換熱器長(zhǎng)期取熱的穩(wěn)定性以及對(duì)其周?chē)鷰r土的熱擾進(jìn)行相關(guān)研究。

3.1 中深層套管式地埋管換熱器的穩(wěn)定性

3.1.1 長(zhǎng)期出水溫度

對(duì)中深層套管式地埋管換熱器的長(zhǎng)期出水溫度進(jìn)行為期15 a的模擬計(jì)算，此時(shí)的延米換熱量為100 W/m，流量為28 m3/h，其余參數(shù)見(jiàn)表1。

圖6所示為中深層套管式地埋管換熱器在15 a運(yùn)行期間的出水溫度。從圖6中可以看出，在每個(gè)供熱周期內(nèi)的中深層套管式地埋管換熱器的出水溫度變化較大，以第1個(gè)供暖期為例，初始階段的出水溫度為35.6 ℃，采暖季結(jié)束后的出水溫度為25.9 ℃。從整體看，中深層套管式地埋管換熱器的出水溫度隨著運(yùn)行年限的增加呈遞減趨勢(shì)，并且遞減程度逐漸減小，第1 a、第5 a、第10 a和第15 a采暖季末期的中深層套管式地埋管換熱器的出水溫度分別為25.9、24.4、23.9、23.5 ℃。分別以第1 a、第5 a和第10 a的出水溫度為基準(zhǔn)，第5 a、第10 a和第15 a的出水溫度在此基礎(chǔ)上分別下降5.8%、2.0%和1.8%。此外，從圖6中也可以看出，從第11 a開(kāi)始出水溫度基本處于準(zhǔn)穩(wěn)態(tài)階段，第11 a、第12 a、第13 a、第14 a和第15 a采暖季末期的中深層套管式地埋管換熱器的出水溫度分別為23.79、23.72、23.66、23.60、23.54 ℃，表明中深層套管式地埋管換熱器長(zhǎng)期運(yùn)行下的出水溫度較為穩(wěn)定。

圖6 15 a內(nèi)的中深層套管式地埋管換熱器出水溫度Fig.6 The oulet temperature of MDBHE within 15 a

3.1.2 長(zhǎng)期巖土溫度變化

對(duì)距地埋管換熱器外壁1 m，不同埋管深度的巖土溫度進(jìn)行為期15 a的模擬計(jì)算，此時(shí)的延米換熱量為100 W/m，流量為28 m3/h，其余參數(shù)見(jiàn)表1。

從圖7可以看出，在15 a運(yùn)行期間，距鉆孔壁1 m處不同深度巖土層的溫度呈周期性變化，并且?guī)r土經(jīng)過(guò)熱恢復(fù)期不能完全恢復(fù)到初始狀態(tài)，與初始溫度相比，巖土溫度隨運(yùn)行年限的增加逐年下降。以z=1 800 m為例進(jìn)行分析，與巖土初始溫度相比，第1 a溫降比例為2.03%，第5 a溫降比例為4.30%，第10 a溫降比例為5.37%，第13～15 a溫降比例分別為5.78%，5.90%和6.01%，溫降增幅逐年減小。分析原因：在運(yùn)行參數(shù)不變的條件下，被取熱后的巖土不能完全恢復(fù)到初始溫度，會(huì)造成下一年地埋管換熱器取熱能力略有下降，進(jìn)而對(duì)其周?chē)鷰r土的熱干擾作用減弱，經(jīng)過(guò)長(zhǎng)時(shí)間的周期循環(huán)這種現(xiàn)象會(huì)逐漸減弱，最后會(huì)趨于相對(duì)平衡，此外，經(jīng)過(guò)為期15 a運(yùn)行，最大溫降比例僅為6.50%，說(shuō)明中深層地埋管換熱器取熱穩(wěn)定性較好，巖土具有較強(qiáng)的熱恢復(fù)能力。

圖7 不同埋管深度下徑向1 m處的巖土溫度變化Fig.7 Temperature variation of rock and soil at 1 m in radial direction under different buried depths

3.2 中深層套管式地埋管換熱器的熱擾

Δtr,z,τ=tbnd,z,0-tr*,z,τ

(9)

式中：tbnd,z,0為τ=0時(shí)刻徑向邊界溫度，℃；tr,z,τ為τ時(shí)刻半徑r處的溫度,℃。

3.2.1 不同運(yùn)行時(shí)間下熱影響半徑

圖8所示為中深層地埋管換熱器運(yùn)行5、10、15 a后的巖土溫度場(chǎng)以及其不同埋管深度處徑向巖土溫度的變化曲線。以圖8(f)為例進(jìn)行分析，當(dāng)運(yùn)行15 a后，與徑向邊界171.47 m處的巖土溫度相比，當(dāng)中深層地埋管換熱器的埋管深度分別為1 200、1 400、1 600、1 800、2 000 m時(shí)，距離其0.17 m處的巖土溫降分別為20.79、23.07、22.91、23.81、24.51 ℃，相比而言z=1 200 m處的溫降最小，表明中深層地埋管換熱器的取熱區(qū)域主要集中在1 400～2 000 m處的巖土。此外，隨著運(yùn)行時(shí)間的增加中深層地埋管換熱器的熱影響半徑逐漸增大，第5、10、15 a的影響半徑分別為39.87、57.42、82.69 m，以第10 a對(duì)應(yīng)的熱影響半徑為基準(zhǔn)，第15 a的熱影響半徑在其基礎(chǔ)上增加44%。分析原因：由前文分析可知經(jīng)取熱后的巖土不能完全恢復(fù)到初始狀態(tài)，并且，隨著運(yùn)行年限的增加，巖土溫度逐漸降低，由于熱傳導(dǎo)的作用，周?chē)词苡绊懙膸r土?xí)虻芈窆苤行倪M(jìn)行熱量傳遞，故熱影響半徑會(huì)隨運(yùn)行時(shí)間的增加而逐漸增大。

圖8 不同運(yùn)行年限下溫度場(chǎng)分布和不同埋管深度徑向溫度的變化Fig.8 Distribution of temperature field and variation of radial temperature at different buried depths under different operation periods

3.2.2 主要設(shè)計(jì)參數(shù)與熱影響半徑的關(guān)系

采用控制變量法分析主要設(shè)計(jì)參數(shù)和中深層地埋管換熱器熱影響半徑之間的關(guān)系，其中，主要設(shè)計(jì)參數(shù)包括管徑、工質(zhì)循環(huán)流量、回填材料導(dǎo)熱系數(shù)以及延米換熱量。在不同設(shè)計(jì)參數(shù)下，系統(tǒng)運(yùn)行10 a后z=2 000 m處徑向巖土溫度的分布如圖9所示。從圖9可以看出，管徑、工質(zhì)循環(huán)流量以及回填材料導(dǎo)熱系數(shù)的變化對(duì)中深層地埋管換熱器周?chē)鷰r土層溫度分布影響較小，相對(duì)于其他3種設(shè)計(jì)參數(shù)，延米換熱量對(duì)中深層地埋管換熱器周?chē)鷰r土層溫度分布影響較大，隨著延米換熱換熱量的增大周?chē)鷰r土層的溫度越低。因此，在進(jìn)行工程設(shè)計(jì)時(shí)應(yīng)采用合理的延米換熱量，文獻(xiàn)[18]對(duì)延米換熱量做了詳細(xì)研究，并指出延米換熱量不宜大于150 W/m。此外，從圖中可以看出在一定范圍內(nèi)，上述4種設(shè)計(jì)參數(shù)雖對(duì)中深層地埋管換熱器周?chē)鷰r土具有一定的影響，但對(duì)其熱半徑的熱影響范圍有限，即中深層地埋管換熱器的熱影響半徑趨于一致。

圖9 不同設(shè)計(jì)參數(shù)下運(yùn)行10 a后徑向巖土溫度的分布Fig.9 Radial temperature distribution of rock and soil under different design parameters after 10 a operation

3.2.3 巖土熱物性參數(shù)和熱影響半徑的關(guān)系

本文主要分析巖土熱物性參數(shù)和熱影響半徑之間的關(guān)系。巖土不同熱物性參數(shù)見(jiàn)表2，此時(shí)的延米換熱量為100 W/m，流量為28 m3/h，其余參數(shù)見(jiàn)表1。

在表2所示的巖土熱物性參數(shù)下，圖10所示為z=500 m處巖土的徑向溫度分布。從圖10中可以看出，與徑向邊界171.47 m處的巖土溫度相比，當(dāng)巖土熱物性參數(shù)分別取案例1～5時(shí)，距離中深層地埋管換熱器0.17 m處的巖土溫降分別為10.55、11.61、16.13、14.94、5.58 ℃，由此可知，中深層地埋管換熱器取熱對(duì)其周?chē)鷰r土溫度場(chǎng)的熱擾主要與巖土本身熱物性參數(shù)相關(guān)。

表2 巖土物理參數(shù)Tab.2 Physical parameters of rock and soil

圖10 不同巖土熱物性參數(shù)下徑向巖土溫度分布Fig.10 Radial temperature distribution along depth of DBHE under different geotechnical thermal physical parameters

案例1～5的巖土單位體積平均熱容分別為2.410×106、1.390×106、2.040×106、2.210×106、2.680×106J/(m3·K)，其中案例2的單位體積平均熱容最小，案例5的單位體積平均熱容最大。當(dāng)z=2 000 m時(shí)，案例2的熱影響半徑最大，而案例5的熱影響半徑最小，其中，案例1～5的巖土平均熱擴(kuò)散率分別為1.43×10-6、2.01×10-6、8.60×10-7、6.19×10-7、4.13×10-7m2/s，所對(duì)應(yīng)的中深層地埋管換熱器熱影響半徑分別為82.69、99.23、68.91、57.43、27.69 m。因此，中深層地埋管換熱器熱影響半徑與巖土單位體積熱容及熱擴(kuò)散率等物性參數(shù)相關(guān)性較大，巖土單位體積熱容越小、熱擴(kuò)散率越大，熱影響半徑也越大。但隨著熱擴(kuò)散率的增大，中深層地埋管換熱器的半徑增加的速率逐漸減小，當(dāng)熱擴(kuò)散率由4.13×10-7m2/s增加到6.19×10-7m2/s時(shí)，熱影響半徑由27.69 m增加到57.43 m；而當(dāng)熱擴(kuò)散率由1.43×10-6m2/s增加到2.01×10-6m2/s時(shí)，熱影響半徑由82.69 m增加到99.23 m，增加率為20%。

為進(jìn)一步探究中深層地埋管換熱器熱影響半徑與巖土熱擴(kuò)散率的關(guān)系，圖11所示為兩者擬合關(guān)系式。由圖11(a)可知，兩者并非簡(jiǎn)單的線性關(guān)系，具體表現(xiàn)為隨巖土熱擴(kuò)散率的增加，中深層地埋管換熱器的熱影響半徑在前期增加速率較快，而后期較慢，滿足指數(shù)關(guān)系式的變化規(guī)律。從圖11(b)中可以看出，當(dāng)熱影響半徑r和熱擴(kuò)散率a的關(guān)系式滿足r=-151.99×exp[-a/(5.14×10-7)]+98.14時(shí)，此時(shí)的R2為0.97，能夠較好的表述兩者之間的關(guān)系。

圖11 巖土熱擴(kuò)散率與中深層地埋管熱影響半徑的關(guān)系Fig.11 Relationship between thermal diffusivity of rock and soil and heat-affected radius of MDBHE

但應(yīng)當(dāng)指出：利用中深層地埋管換熱器熱影響半徑確定鉆孔間距并不是最好的方法，熱影響半徑和間距問(wèn)題并不完全等價(jià)，最好的方法是直接求解多根中深層地埋管換熱器之間耦合，而且要求解反問(wèn)題，即算出取熱量隨間距的變化，當(dāng)然難度和計(jì)算工作量較大?，F(xiàn)有文獻(xiàn)的報(bào)道多關(guān)于單根中深層地埋管換熱器，但在實(shí)際工程中，延米換熱量取值一般要小于150 W/m，單根中深層地埋管換熱器不能滿足現(xiàn)有公共建筑及住宅建筑的供暖需求，往往需要更多根中深層地埋管換熱器來(lái)滿足建筑總供熱負(fù)荷。因此，為使研究成果更好的服務(wù)實(shí)際工程，應(yīng)研究中深層地埋管換熱器管群，分析管群的布管形式以及地埋管換熱器之間的相互影響，但在現(xiàn)有研究背景下，本研究成果仍有一定的參考意義。

4 結(jié) 論

1)中深層套管式地埋管換熱器長(zhǎng)期運(yùn)行下的出水溫度較為穩(wěn)定。從第11 a開(kāi)始出水溫度基本處于準(zhǔn)穩(wěn)態(tài)階段，第11 a、第12 a、第13 a、第14 a和第15 a采暖季末期的中深層套管式地埋管換熱器的出水溫度分別為23.79、23.72、23.66、23.60、23.54 ℃。

2)中深層套管式地埋管換熱器周?chē)鷰r土熱恢復(fù)能力較強(qiáng)。對(duì)距地埋管換熱器外壁1 m位置處不同埋管深度的巖土溫度進(jìn)行為期15 a計(jì)算，結(jié)果顯示經(jīng)過(guò)15 a的運(yùn)行，與初始巖土溫度相比最大溫降比例僅為6.50%。

3)中深層套管式地埋管換熱器的熱影響半徑隨運(yùn)行時(shí)間的增加而逐漸增大。運(yùn)行5、10、15 a的影響半徑分別為39.87、57.42、82.69 m，以第10 a對(duì)應(yīng)的熱影響半徑為基準(zhǔn)，5 a后的熱影響半徑在其基礎(chǔ)上增加44%。

4)相比于中深層套管式地埋管換熱器外管管徑、回填材料導(dǎo)熱系數(shù)及工質(zhì)循環(huán)流量，延米換熱量對(duì)中深層換熱器周?chē)鷰r土的溫度分布影響較大。但在一定范圍內(nèi)，上述設(shè)計(jì)參數(shù)和中深層套管式地埋管換熱器熱影響半徑之間無(wú)明顯關(guān)系。

5)隨巖土熱擴(kuò)散率的增加，中深層套管式地埋管換熱器的熱影響半徑在前期增加速率較快，而后期較慢，當(dāng)熱擴(kuò)散率由4.13×10-7m2/s增加到6.19×10-7m2/s時(shí)，熱影響半徑由27.69 m增加到57.43 m；而當(dāng)熱擴(kuò)散率由1.43×10-6m2/s增加到2.01×10-6m2/s時(shí)，熱影響半徑由82.69 m增加到99.23 m，增長(zhǎng)率為20%，滿足指數(shù)關(guān)系式的變化規(guī)律。并且，熱影響半徑r和熱擴(kuò)散率a的關(guān)系式滿足r=-151.99×exp[-a/(5.14×10-7)]+98.14，其中，R2=0.97。