李小玲，馬貴陽(yáng)，趙 鵬

（1. 遼寧石油化工大學(xué)石油天然氣工程學(xué)院，遼寧 撫順113001； 2. 中石化管道儲(chǔ)運(yùn)公司京唐輸油處，江蘇 徐州 221000）

地源熱泵埋管換熱器數(shù)值模擬研究

李小玲1，馬貴陽(yáng)1，趙 鵬2

（1. 遼寧石油化工大學(xué)石油天然氣工程學(xué)院，遼寧 撫順113001； 2. 中石化管道儲(chǔ)運(yùn)公司京唐輸油處，江蘇 徐州 221000）

針對(duì)地源熱泵垂直U型埋地?fù)Q熱器的傳熱特性，建立了換熱器及周?chē)寥罍囟葓?chǎng)的數(shù)學(xué)物理模型，模擬計(jì)算分析了單鉆孔垂直U型埋地?fù)Q熱器短期運(yùn)行工況下系統(tǒng)運(yùn)行時(shí)間和不同灌漿材料對(duì)埋管周?chē)寥罍囟葓?chǎng)的影響。研究結(jié)論可為合理設(shè)計(jì)埋管換熱器提供參考。

地源熱泵；運(yùn)行時(shí)間；灌漿材料；溫度場(chǎng)

地源熱泵系統(tǒng)具有運(yùn)行穩(wěn)定、溫度變化范圍小等優(yōu)點(diǎn)。用土壤作熱泵的熱源，夏季土壤溫度比室外氣溫低，冬季比室外氣溫高，這樣更有利于提高熱泵系統(tǒng)的運(yùn)行性能[1,2]。地源熱泵系統(tǒng)的性能很大程度上取決于埋管換熱器的性能，而且熱泵系統(tǒng)的負(fù)荷、初投資造價(jià)均與埋管換熱器的尺寸相關(guān)聯(lián)。

采用有限元計(jì)算方法，編制程序，對(duì)垂直U型埋管換熱器運(yùn)行工況進(jìn)行模擬計(jì)算。研究分析了單鉆孔垂直U型埋地?fù)Q熱器短期運(yùn)行工況下系統(tǒng)運(yùn)行時(shí)間和不同灌漿材料對(duì)埋管周?chē)寥罍囟葓?chǎng)的影響。

1 傳熱模型

垂直單U型埋管換熱器，管子埋在地表下的豎直鉆孔中，用回填材料填實(shí)，流體在U型埋管內(nèi)流動(dòng)并與土壤進(jìn)行換熱。

模型的建立基于如下假設(shè)：忽略地表面溫度波動(dòng)對(duì)土壤溫度的影響，認(rèn)為土壤初始溫度均勻一致，恒定不變，且等于邊界土壤的溫度；忽略埋管和回填土接觸熱阻；不考慮管壁熱阻；不考慮水分遷移對(duì)熱量傳遞的影響。

控制方程：

式中：λ—土壤的導(dǎo)熱系數(shù)，W/(m·K)

ρ—土壤的密度，kg/m3

C—土壤的比熱，J/(kgC)?o

邊界條件：

（1）管內(nèi)流體與埋管壁交界施加第三類(lèi)邊界條件

式中：α—管璧與防凍液換熱系數(shù)，W/（m2·℃）

T—管璧的溫度，℃

Tw—管內(nèi)流體的溫度，℃

（2）計(jì)算土壤遠(yuǎn)邊界施加絕熱邊界條件

式中：n為邊界外法線(xiàn)方向。

（3）邊界土壤溫度保持土壤初始溫度不變

式中：T∞—土壤初始溫度，℃

初始條件：

網(wǎng)格劃分：用三節(jié)點(diǎn)單元對(duì)U型管周?chē)靥钔良巴寥肋M(jìn)行單元?jiǎng)澐?。由于地下?lián)Q熱器傳熱過(guò)程中，溫度沿徑向方向變化較大，因此在水平方向上對(duì)U型管周?chē)W(wǎng)格進(jìn)行了局部加密。如圖1所示。

圖1 土壤溫度場(chǎng)水平面網(wǎng)格Fig.1 Horizontal plane of underground temperature field

土壤計(jì)算面積為3 m×3 m矩形區(qū)域，鉆孔孔徑300 mm,U型管外徑32 mm，埋管兩支管中心間距180 mm土壤導(dǎo)熱系數(shù)1.2 W/(m·K)，土壤密度1 925 kg/m3，土壤比熱1 039 J/（kg·℃），土壤遠(yuǎn)邊界溫度15 ℃，回填材料導(dǎo)熱系數(shù)0.75 W/(m·K)，回填材料比熱2 026 J/（kg·℃），熱泵在制冷模式下，防凍液進(jìn)出口溫度分別為37、33 ℃（進(jìn)出口液體最大溫差大約為5℃），管內(nèi)璧與防凍液換熱系數(shù)2 300 W/(m2·℃)，土壤初始溫度15 ℃。由于一般換熱器的水平管都設(shè)在地下1~2 m處，因此本節(jié)研究垂直深度為2 m處換熱器周?chē)寥浪椒较驕囟葓?chǎng)的分布情況。

2 系統(tǒng)運(yùn)行時(shí)間的影響

模擬了系統(tǒng)運(yùn)行240 h內(nèi)換熱器周?chē)寥罍囟葓?chǎng)。從圖2可以看出，隨著系統(tǒng)運(yùn)行時(shí)間的增加，土壤溫度升高，但各點(diǎn)的溫度變化規(guī)律不盡相同。

圖2 土壤溫度隨系統(tǒng)運(yùn)行時(shí)間的變化Fig.2 The change of soil temperature with the system operating time

土壤換熱器附近的點(diǎn)，如左孔壁（x=1.35 m）、鉆孔中心（x=1.5 m）及右孔壁（x=1.65 m），系統(tǒng)運(yùn)行初期土壤溫度急劇上升，但隨著運(yùn)行時(shí)間的增加，這些點(diǎn)的溫度逐漸趨于穩(wěn)定，溫度緩慢增加。而離土壤換熱器較遠(yuǎn)的點(diǎn)，如x=0.5 m、x=2.5 m，系統(tǒng)運(yùn)行初期土壤溫度基本不變，隨著運(yùn)行時(shí)間的增加，土壤溫度穩(wěn)步上升，上升的幅度越來(lái)越大。

為了更好的說(shuō)明問(wèn)題，在此定義Tg為土壤溫度隨系統(tǒng)運(yùn)行時(shí)間的變化率：

式中:τ—系統(tǒng)運(yùn)行時(shí)間，h。

將系統(tǒng)運(yùn)行時(shí)間分別由10, 40, 80, 120, 160, 200 h增加到40, 80, 120, 160, 200，240 h，分別稱(chēng)為變化a,b,c,d,e,f。圖3給出了土壤溫度變化率隨x的變化。由該圖可以看出，系統(tǒng)運(yùn)行的前40 h（變化a）土壤換熱器附近的溫度變化率較大，最大值出現(xiàn)在左孔壁（x=1.35 m）處，Tg達(dá)到了0.092 7，而離換熱器較遠(yuǎn)處土壤溫度變化率很小，x=2.5 m處Tg為0.001 2，Tg的最大和最小值相差0.091 5。系統(tǒng)運(yùn)行40 h后便較快的趨于穩(wěn)定，Tg的最大值急劇降低且向孔壁外轉(zhuǎn)移，如變化b，土壤溫度變化率的最大值出現(xiàn)在x=1 m處，Tg為0.028 4，Tg的最大和最小值相差0.022 5。這是因?yàn)殡S系統(tǒng)運(yùn)行時(shí)間的增加熱量開(kāi)始向外傳遞導(dǎo)致土壤溫度的升高。直到系統(tǒng)運(yùn)行240 h（變化f），Tg的最大和最小值僅相差0.006，系統(tǒng)幾乎處于穩(wěn)定狀態(tài)，熱量近似達(dá)到平衡，土壤溫度雖仍在升高但上升的幅度逐漸減小，最終趨于零。

圖3 土壤溫度變化率隨徑向x的變化Fig.3 The change rate of soil temperature along the radial x

3 鉆孔中灌漿材料的影響

相比于土壤，回填材料充當(dāng)換熱器和土壤之間的橋梁，使兩者之間可以更流暢的進(jìn)行熱交換，它的作用更多的應(yīng)該體現(xiàn)在它的熱量傳遞方面，所以它的導(dǎo)熱系數(shù)就顯得尤為重要。對(duì)導(dǎo)熱系數(shù)分別為0.75 W/(m·K)、1.10 W/(m·K)、1.50 W/(m·K)、2.08 W/(m·K)、2.42 W/(m·K)[3]這五種回填材料對(duì)土壤溫度場(chǎng)的影響進(jìn)行了模擬分析。圖4為回填材料不同的系統(tǒng)運(yùn)行240 h后土壤沿徑向溫度分布對(duì)比圖。

圖4 回填材料導(dǎo)熱系數(shù)不同的系統(tǒng)運(yùn)行240 h時(shí)后土壤沿徑向溫度分布Fig.4 The soil temperature distribution along the radial after 240 hours of system operation when thermal conductivity of backfill materials is varied

對(duì)比圖4中5種導(dǎo)熱系數(shù)不同的回填材料在徑向的溫度分布也不難看出，隨著回填材料導(dǎo)熱系數(shù)的增大，雖然增大了整個(gè)傳熱效果，但同時(shí)也增大了兩管間的熱短路[4,5]。如：以回填土導(dǎo)熱系數(shù)分別為0.75 W/(m·K)和2.42 W/(m·K)的兩條曲線(xiàn)為例（進(jìn)液管中心在徑向1.41 m處，出液管中心在徑向1.59 m處），前者在徑向1.461 m和1.563 m處土壤溫度為31 ℃，而后者在徑向1.513 m和1.548 m處土壤溫度為31 ℃，說(shuō)明回填土導(dǎo)熱系數(shù)越大，兩管間的熱短路現(xiàn)象也就越明顯。

下面通過(guò)幾個(gè)點(diǎn)的分析說(shuō)明問(wèn)題：取y=1.5 m，x分別為1.00 m、1.35 m（左孔壁）、1.5 m（鉆孔中心）、1.65 m（右孔壁）及2.00 m的5個(gè)點(diǎn)研究土壤水平方向的溫度場(chǎng)隨回填材料導(dǎo)熱系數(shù)變化的分布情況。

圖5給出了各點(diǎn)土壤溫度隨回填材料導(dǎo)熱系數(shù)的變化。

圖5 土壤溫度隨回填材料導(dǎo)熱系數(shù)的變化Fig.5 The change of soil temperature with thermal conductivity of backfill materials

由圖可以看出，如前所述，隨著回填材料導(dǎo)熱系數(shù)的增加，土壤溫度增加。但各點(diǎn)溫度增加的變化規(guī)律不盡相同。在靠近換熱器處，土壤溫度增幅較小，而遠(yuǎn)離換熱器處，土壤溫度增幅較大。如，當(dāng)回填材料導(dǎo)熱系數(shù)由0.75 W/(m·K)增加到2.42 W/(m·K)時(shí)，x=1.35 m（左孔壁）處土壤溫度由28.575℃增加到了31.307 ℃，增加了2.732 ℃，而在x=1.00 m處土壤溫度由19.253 ℃增加到了23.312 ℃，增加了4.059 ℃。

為了更好的說(shuō)明問(wèn)題，在此定義Tg為土壤溫度隨回填材料導(dǎo)熱系數(shù)的變化率：

式中:λ1—回填材料導(dǎo)熱系數(shù)，W/(m·K)。

將回填材料導(dǎo)熱系數(shù)分別由0.75, 1.1, 1.5, 2.08 W/(m·K)增加到1.1, 1.5, 2.08, 2.42 W/(m·K), 分別稱(chēng)為變化a,b,c,d。圖6給出了土壤溫度變化率隨x的變化。由該圖可知，隨著回填材料導(dǎo)熱系數(shù)的增加土壤溫度變化率在鉆孔中心左側(cè)是減小的，而在鉆孔中心右側(cè)是增加的。而且，在鉆孔內(nèi)Tg急劇變化，而在鉆孔外Tg變化率減小。例如，在x分別為1.00, 1.35, 1.5 m時(shí)，變化a的Tg分別為3.640, 2.517, 1.626，是逐漸減小的，而且x由1.00 m到1.35 m（土壤內(nèi)），Tg的變化率為3.21，由1.35 m到1.5 m（鉆孔內(nèi)），Tg的變化率為5.94。這主要是因?yàn)榛靥畈牧系膶?dǎo)熱系數(shù)主要影響的是回填材料的熱傳遞速度，而土壤的熱物性不變，所以離回填材料越遠(yuǎn)土壤溫度變化所受的影響越小。

圖6 土壤溫度變化率隨徑向x的變化Fig.6 The change rate of soil temperature along the radial x

此外，比較變化a到變化d的4條曲線(xiàn)亦可看出，離鉆孔中心距離相同處，隨著回填材料導(dǎo)熱系數(shù)的增加，Tg值是減小的。如，當(dāng)x=1.35 m時(shí)，變化a的Tg為2.517，而變化b的Tg為1.855，變化c的Tg為1.326，而變化d的Tg為1.000。說(shuō)明回填物導(dǎo)熱系數(shù)增加可改善熱傳導(dǎo)性能，但是隨著回填物導(dǎo)熱系數(shù)進(jìn)一步增加，導(dǎo)熱的增加率卻遞減。這主要是因?yàn)殡S著回填材料導(dǎo)熱系數(shù)的增加，系統(tǒng)在同一運(yùn)行時(shí)間內(nèi)進(jìn)入土壤的熱量傳遞的更遠(yuǎn)，如圖，當(dāng)Tg=2.0時(shí)，對(duì)應(yīng)變化a的x約為1.43 m，而對(duì)應(yīng)變化b的x約為1.30 m，對(duì)應(yīng)變化c的x為1.00 m（鉆孔中心處x為1.5 m）。

綜上所述，盡管回填在埋管附近的回填材料比土壤的體積小得多，但對(duì)埋管和土壤的換熱影響很大。因此，地埋管地源熱泵換熱器鉆孔回填材料應(yīng)該在滿(mǎn)足現(xiàn)場(chǎng)施工可行的基礎(chǔ)上具有較高的導(dǎo)熱性，以提高地下?lián)Q熱能力。高性能的回填材料可以通過(guò)減小地埋管鉆孔深度等方法降低地源熱泵的初裝成本以及地源熱泵后期的運(yùn)行維護(hù)費(fèi)用。

4 結(jié) 論

通過(guò)模擬計(jì)算分析單鉆孔垂直U型埋地?fù)Q熱器短期運(yùn)行工況下系統(tǒng)運(yùn)行時(shí)間和不同灌漿材料對(duì)埋管周?chē)寥罍囟葓?chǎng)的變化情況，得出如下結(jié)論：

（1）隨著系統(tǒng)運(yùn)行時(shí)間的增加，土壤溫度升高，但各點(diǎn)的溫度變化規(guī)律不盡相同。土壤換熱器附近的點(diǎn)，系統(tǒng)運(yùn)行初期土壤溫度急劇上升，但隨著運(yùn)行時(shí)間的增加，這些點(diǎn)的溫度逐漸趨于穩(wěn)定，溫度緩慢增加。而離土壤換熱器較遠(yuǎn)的點(diǎn)，系統(tǒng)運(yùn)行初期土壤溫度基本不變，隨著運(yùn)行時(shí)間的增加，土壤溫度穩(wěn)步上升，上升的幅度越來(lái)越大。本文模擬結(jié)果表明，系統(tǒng)運(yùn)行40 h后便較快的趨于穩(wěn)定，Tg的最大值急劇降低且向孔壁外轉(zhuǎn)移。直到系統(tǒng)運(yùn)行240 h時(shí)Tg的最大和最小值僅相差0.006，系統(tǒng)幾乎處于穩(wěn)定狀態(tài)，熱量近似達(dá)到平衡，土壤溫度雖仍在升高但上升的幅度逐漸減小，最終將趨于零。

（2）就本文研究的5種灌漿材料而言，其導(dǎo)熱系數(shù)對(duì)埋管和土壤的換熱影響很大。當(dāng)回填材料導(dǎo)熱系數(shù)由0.75 W/(m·K)增加到2.42 W/(m·K)時(shí)，x=1.35 m（左孔壁）處土壤溫度增加了2.732 ℃，而在x=1.00 m（左孔壁外）處土壤溫度增加了4.059℃。這充分說(shuō)明：①隨著回填材料導(dǎo)熱系數(shù)的增加，土壤溫度增加，Tg（土壤溫度變化率）為正值。但各點(diǎn)溫度增加的變化規(guī)律不盡相同。在靠近換熱器處，土壤溫度增幅較小，而遠(yuǎn)離換熱器處，土壤溫度增幅較大。②此外，模擬結(jié)果也表明，對(duì)同一支管間距，隨著回填材料導(dǎo)熱系數(shù)的增大，雖然增大了整個(gè)傳熱效果，但同時(shí)也增大了兩管間的熱短路。③同時(shí)，回填物導(dǎo)熱系數(shù)增加雖然可改善熱傳導(dǎo)性能,但是隨著回填物導(dǎo)熱系數(shù)進(jìn)一步增加,導(dǎo)熱的增加率卻遞減。因此，回填材料應(yīng)該在滿(mǎn)足現(xiàn)場(chǎng)施工可行的基礎(chǔ)上具有較高的導(dǎo)熱性，同時(shí)盡可能增大埋管間距。

［1］裴俠風(fēng).地源熱泵技術(shù)的應(yīng)用現(xiàn)狀及展望[J].制冷與空調(diào)，2004，3：76-78.

［2］李高建，胡玉葉.地源熱泵技術(shù)的研究與應(yīng)用現(xiàn)狀[J].節(jié)能技術(shù), 2007，25（142）：176-178.

［3］高青，李明,等.地下群井換熱強(qiáng)化與運(yùn)行模式影響規(guī)律[J].太陽(yáng)能學(xué)報(bào), 2006,27（1）：83-89.

［4］李新國(guó)，汪洪軍，趙軍，等.不同回填材料對(duì)U型垂直埋管換熱性能的影響[J].太陽(yáng)能學(xué)報(bào)，2003，24（6）：810-813.

［5］楊衛(wèi)波，施明恒.地源熱泵中U型埋管傳熱過(guò)程的數(shù)值模擬[J].東南大學(xué)學(xué)報(bào),2007,37（1）：78-83.

Numerical Simulation Study on Heat Exchanger for Ground-source Heat Pump

LI Xiao-ling1，MA Gui-yang1，ZHAO Peng2
（1. College of Petroleum and Natural Gas Engineering，Liaoning Shihua University, Liaoning Fushun 113001，China；2. Sinopec Pipeline Storage and Transportation Branch，Jiangsu Xushou 221000，China）

The heat transfer model of underground vertical buried tube heat exchanger was built with considering the heat transfer characteristics of U-tubes ground heat exchanger of ground-source heat pump. Then influences of the system operating time and borehole’s backfill material on the soil temperature field were analyzed. The research results may provide a guidance in the design of ground coupled heat pump systems.

Ground-source heat pump; Running time; Grouting material; Temperature distribution

TQ 019

A

1671-0460（2011）11-1198-04

2011-08-30

李小玲（1982-），女，青海西寧人，講師，博士在讀，研究方向：從事地源熱泵技術(shù)工作。E-mail：lingzi_1015@163.com。