穆 玄 鄒 行 裴 鵬 林華穎 郝定溢

（1.貴州大學(xué)礦業(yè)學(xué)院 貴陽 550025;2.中國礦業(yè)大學(xué)礦業(yè)工程學(xué)院 徐州 221116）

0 引言

地源熱泵是以淺層地表熱能為熱源，通過輸入少量的高品位能源，實(shí)現(xiàn)低品位熱能向高品位熱能轉(zhuǎn)移的熱泵系統(tǒng)。地源熱泵屬于可再生能源利用技術(shù)，對(duì)能量的利用具有經(jīng)濟(jì)高效、無污染、運(yùn)行費(fèi)用低，不受地質(zhì)條件限制等優(yōu)點(diǎn)，被認(rèn)為是一種極具潛力的綠色能源技術(shù)[1-4]。

地埋管地源熱泵是向土壤或把土壤作為熱交換器來傳輸熱量，地下?lián)Q熱器是系統(tǒng)中進(jìn)行熱量交換的主要設(shè)備[5]，其主要是通過埋管內(nèi)流體與埋管周圍土壤之間進(jìn)行熱量傳遞來實(shí)現(xiàn)熱量的轉(zhuǎn)移，因此其傳熱問題一直被諸多學(xué)者關(guān)注研究[6-8]，

水平埋管式地源熱泵系統(tǒng)其管體布置于淺層土壤之中，其換熱性能受土壤熱參數(shù)影響較大[9]?；诖藝鴥?nèi)外學(xué)者做了大量研究，Leong W H 等[10]通過計(jì)算機(jī)模擬使用了五種不同飽和度（0、12.5、25、50 和100%）的土壤（砂、粉質(zhì)壤土和粉質(zhì)粘土），指出土壤熱泵系統(tǒng)的換熱性能受土壤含水率及土壤類型影響。Hikari 等[11,12]使用FEELOW 模擬軟件對(duì)水平“slinky”型螺旋管的優(yōu)化設(shè)計(jì)進(jìn)行了數(shù)值模擬研究，結(jié)果顯示，模擬值與實(shí)驗(yàn)值有較好的一致性，證明了該數(shù)值模型可用于模擬研究水平“slinky”型螺旋管的換熱問題以及用于其優(yōu)化設(shè)計(jì)。

然而，之前的研究中，多將非飽和區(qū)域土壤不同深度含水飽和度視作常數(shù)，往往忽略了地下水位線以上區(qū)域毛管水分布特征，從而使得其模擬結(jié)果與實(shí)際情況偏差較大，這在一定程度上不利于實(shí)際工程指導(dǎo)。因此該文基于土壤毛管水力特征曲線對(duì)土壤中非飽和區(qū)域水分的分布情況做了定量描述，并討論了不同因素對(duì)水平管換熱性能的影響，其研究結(jié)果可用于更準(zhǔn)確地計(jì)算水平埋管換熱器的換熱性能和估算地源熱泵機(jī)組的技術(shù)經(jīng)濟(jì)指標(biāo)。

1 理論及控制方程

1.1 多孔介質(zhì)中的水力傳導(dǎo)

達(dá)西定律作為描述飽和巖土體中水分流動(dòng)的經(jīng)典定律，指出巖土體中的水分在壓力梯度下由高水勢(shì)流向低水勢(shì)流動(dòng)，方程如下[13,14]：

式中：v為達(dá)西速度，m/s；κ為多孔介質(zhì)的滲透系數(shù)，m2；μ為流體的動(dòng)力粘度，Pa·s；▽?duì)诪槿我鈨牲c(diǎn)水勢(shì)差；ψg為重力勢(shì)；ψp為壓力勢(shì)。

同樣的，將裂隙看做為多孔介質(zhì)，裂隙流方程使用了達(dá)西定律的切向?qū)?shù)來定義沿內(nèi)部邊界（多孔介質(zhì)中的裂隙）的流動(dòng)。

式中：κf為是裂隙的滲透率，m2；df為裂隙張開度或厚度，m；u為裂隙中的達(dá)西流速，m/s。

1.2 多孔介質(zhì)中的熱量傳遞

多孔介質(zhì)中的能量守恒方程為[15]：

式中：ρ為土壤基質(zhì)各相態(tài)的密度，kg/m3；Cp為土壤基質(zhì)各相態(tài)質(zhì)量比熱容，J/kg·K；T為溫度，K；t為時(shí)間，s；ρw為地下水的密度，kg/m3；Cp,w為地下水的質(zhì)量比熱容，J/kg·K；v為土壤基質(zhì)中流體速度，m/s（該速度場(chǎng)耦合于達(dá)西速度）；keff為多孔介質(zhì)有效導(dǎo)熱系數(shù)，W/m·K；(ρCp)eff為多孔介質(zhì)的有效體積比熱容，J/m3·K；Qwall為與埋管換熱器的換熱量，W/m3。

對(duì)裂隙中的熱傳遞的控制方程為：

式中：(ρCp)’eff為裂隙體積熱容，J/m3·K；qfr為裂隙流體的傳導(dǎo)熱通量，W/m2；u為裂隙中達(dá)西速度，m/s。k’eff為裂隙有效導(dǎo)熱系數(shù)，W/m·K；Q為可能存在的熱源，W/m3；q0為流體和固體邊界上的熱通量，W/m2；?t是切向梯度因子。

在此將巖體以及裂隙視作為多相介質(zhì)，主要由固相、氣相、液相等組成，表示為：

固相體積分?jǐn)?shù)：

氣相體積分?jǐn)?shù)：

液相體積分?jǐn)?shù)：

式中：η為對(duì)應(yīng)巖體或裂隙的孔隙度。θ為水分體積分?jǐn)?shù)。

多孔介質(zhì)有效導(dǎo)熱系數(shù)keff與有效體積比熱容(ρCp)eff可基于不同相態(tài)表示為[17,18]：

式中：Ki為為多孔介質(zhì)或裂隙中固相、氣相、液相的導(dǎo)熱系數(shù)，W/m·K；ρi為為多孔介質(zhì)或裂隙中固相、氣相、液相的密度，J/m3·K；Ci為多孔介質(zhì)或裂隙中固相、氣相、液相的質(zhì)量熱容，J/Kg·K。

1.3 管道中的熱傳遞

U 型管中不可壓縮流體流動(dòng)的能量方程如下[16]：

式中，ρf為管內(nèi)流體密度，kg/m3；A為流體流過的斷面面積，m2；Cp,f為管中流體的質(zhì)量比熱容，J/kg·K；Tf為循環(huán)流體的溫度，K；t為時(shí)間，s；u’為水平管中換熱工質(zhì)的速度，m/s；kf為循環(huán)流體導(dǎo)熱系數(shù)，W/m·K；fD為達(dá)西摩擦系數(shù)；dh為平均水力直徑，mm。

U 型管與周圍巖體熱量交換方程如下：

式中：(hZ)eff為熱傳遞系數(shù)有效值（其中h為管壁等效對(duì)流傳熱系數(shù)，W/m2·K；Z為管壁周長，m）；Text為管壁外多孔介質(zhì)的溫度，K。

在巖溶管道內(nèi)的地下水傳熱方程類似于工質(zhì)在U 型管內(nèi)的能量控制方程，其不同點(diǎn)主要在于巖溶管道將巖體視為管壁。

2 不同地下水位線深度對(duì)水平管換熱性能的影響

2.1 仿真模型構(gòu)建

基于土壤毛管水理論，該節(jié)主要考慮了三種水力模型。

地下水位線位于埋管上方，水位線距地表0m，埋管水平孔隙含水飽和度（水所占的孔隙的體積與巖石孔隙體積之比）100%，含水體積分?jǐn)?shù)0.434（見圖1）。

圖1 100%含水飽和度模型X-Z 切面圖Fig.1 X-Z cross-sectional view of 100%saturation model

地下水位線位于埋管下方，水位線距地表8m，埋管水平孔隙含水飽和度50%，含水體積分?jǐn)?shù)0.218（見圖2）。

圖2 50%含水飽和度模型X-Z 切面圖Fig.2 X-Z cross-sectional view of 50%saturation model

地下水位線位于埋管下方，水位線距地表200m，埋管水平孔隙含水飽和度12%，含水體積分?jǐn)?shù)0.052（視為殘余水量）（見圖3）。

表2 三種土壤的熱力學(xué)性能Table 2 Thermodynamic properties of three soils

圖3 12%含水飽和度模型X-Z 切面圖Fig.3 X-Z cross-sectional view of 12%saturation model

根據(jù)三種水力模型，模型的基本參數(shù)如表1所示。

表1 模擬參數(shù)Table 1 Simulation parameters

2.2 模擬結(jié)果與分析

模擬一個(gè)制冷期，后文中截取的溫度場(chǎng)選取時(shí)間點(diǎn)為制冷末期（第100 天），X-Y 切面位于埋管水平，將同時(shí)展示3 種不同含水飽和度分別在不同時(shí)間點(diǎn)的俯視溫度云圖，具體俯視切面位置如圖4所示。

圖4 X-Y 切面展示圖Fig.4 X-Y section display

2.2.1 制冷末期（第80 天）X-Y(Z=-2)切面溫度場(chǎng)分布情況

隨著水平管內(nèi)高溫工質(zhì)持續(xù)向土壤中排入熱量，埋管周圍土壤溫度均呈升高趨勢(shì)。對(duì)比圖5（a）、（b）、（c）可知，土壤含水飽和度越高，受排熱影響范圍越小，隨著土壤含水飽和度的降低，埋管區(qū)域土壤“熱堆積”現(xiàn)象加劇，表現(xiàn)為更多的高溫黃色區(qū)域。數(shù)據(jù)表明，土壤含水飽和度100%時(shí)，埋管附近最高溫度升至302.5K，土壤含水飽和度50%時(shí)，埋管附近最高溫度升至302.9K，土壤含水飽和度12%時(shí)，埋管附近最高溫度升至303.3K，比含水飽和度100%時(shí)多升高了0.8K，如圖5所示。

圖5 制冷末期溫度場(chǎng)分布圖Fig.5 Temperature field distribution diagram at the end of refrigeration

2.2.2 水平管出口水溫

制冷期內(nèi)隨著土壤含水飽和度的增加，其出水溫度越來越低，意味著埋管向土壤中提取了更多的“冷量”。數(shù)據(jù)表明，在制冷末期，土壤含水飽和度100%時(shí)，出口水溫301.7K，埋管進(jìn)出水溫差高達(dá)4.3K。土壤含水飽和度50%時(shí)，出水口溫302.3K，埋管進(jìn)出水溫差3.7K。當(dāng)土壤含水飽和度降低至12%時(shí)，其出口水溫升高至302.7K，較含水飽和度100%時(shí)升高了23.2%，埋管進(jìn)出口水溫差3.3K，如圖6所示。

圖6 地埋管出口水溫Fig.6 Outlet water temperature of buried pipe

2.2.3 延米換熱量

埋管與周圍土壤熱量交換能力在制冷末期趨于平穩(wěn)，此時(shí)，埋管延米換熱量基于土壤含水飽和度的不同而存在差異。土壤含水飽和度100%時(shí)，水平管延米換熱量35W/m，土壤含水飽和度50%時(shí)，水平管延米換熱量30W/m，隨著含水飽和度降低至12%時(shí)，水平管延米換熱量下降至27W/m。土壤含水飽和度100%較含水飽和度12%每延米換熱量提高了30%，收益明顯，其原因主要在于土壤含水飽和度越高，比熱容更大，容納熱量的能力較強(qiáng)，因此在排入相同的熱量時(shí)，含水飽和度高的土壤溫升幅度會(huì)更小，如圖7所示。

圖7 埋管延米換熱量Fig.7 Heat transfer per meter of buried pipe

3 不同土壤類型對(duì)水平管換熱性能的影響

3.1 仿真模型構(gòu)建

該節(jié)主要考慮了三種土壤類型，即砂土、壤土、黏土，具體模型如圖8所示。

圖8 三種土壤類型及模型Fig.8 Three soil types and their models

根據(jù)三種土壤類型，模擬參數(shù)設(shè)置如表3所示。

表3 模擬參數(shù)Table 3 Simulation parameters

3.2 模擬結(jié)果與分析

3.2.1 溫度場(chǎng)分布

在制冷末期，壤土的溫度場(chǎng)分布情況優(yōu)于比熱容最大的黏土及導(dǎo)熱系數(shù)最大的砂土。數(shù)據(jù)表明，在制冷期第90 天，砂土的土壤最高溫度升至304.7 K，壤土的土壤溫度升至304.2 K，黏土的土壤溫度升至304.8 K，相較于砂土及黏土而言，壤土的土壤溫度場(chǎng)分布情況最佳，熱堆積程度最輕，如圖9所示。

圖9 制冷期末期溫度場(chǎng)分布：（a 黏土、b 壤土、c 砂土）Fig.9 Distribution of temperature field at the end of the refrigeration period:(a clay,b loam,c sandy soil)

3.2.2 水平埋管出口水溫

在制冷期第90 天，砂土出水溫度303.8 K，壤土303.4 K，黏土304 K。水平管在壤土中有更低的出口水溫，制冷效率最佳。對(duì)于土壤導(dǎo)熱系數(shù)最大的砂土而言，其出水溫度相較于黏土有一定的下降趨勢(shì)，在一定程度上可認(rèn)為增加土壤的導(dǎo)熱系數(shù)比增加土壤中水分能夠獲得更直觀的效益。所列三種土壤中，導(dǎo)熱系數(shù)及比熱容適中的壤土具有更低的出水溫度，意味著其與周圍土壤的熱交換過程最為利。因此，影響埋管與土壤的換熱過程并非完全取決于土壤中液相及固相所占比例，而是在兩者共同的影響機(jī)理作用下進(jìn)行熱量交換，只有在土壤固相所占比例較大的情況下，增加土壤中水分才能獲得更大的換熱效率，如圖10所示。

圖10 地埋管出口水溫Fig.10 Outlet water temperature of buried pipe

4 不同換熱器構(gòu)造對(duì)水平管換熱性能的影響

4.1 仿真模型構(gòu)架

基于目前常用的幾種水平管換熱器構(gòu)造，該節(jié)主要考慮了三種水平管模型，如圖11所示。

圖11 三種水平管埋管模型Fig.11 Three models of buried horizontal pipes

根據(jù)三種水平管模型，模型參數(shù)設(shè)置如表4所示。

表4 模擬參數(shù)Table 4 Simulation parameters

4.2 模擬結(jié)果與分析

4.2.1 溫度場(chǎng)分布情況

采用水平直管布置方式，在經(jīng)歷90 天的制冷期后，埋管附近土壤溫度最高升至304.4 K，水平“spiral”型螺旋埋管附近土壤溫度升高至306.9 K，水平“slinky”型螺旋埋管附近土壤溫度升高至307.6 K，比水平直管溫度多升高了3.2 K。同時(shí)可明顯觀測(cè)到水平“slinky”型螺旋埋管附近土壤熱堆積現(xiàn)象嚴(yán)重，溫度場(chǎng)分布情況較差，換熱性能受之影響較大。相較而言，水平直管周圍土壤溫度場(chǎng)分布情況最佳，水平“spiral”型螺旋埋管附近土壤溫度分布情況欠佳。

圖12 三種埋管方式溫度場(chǎng)變化情況Fig.12 Variation of temperature field in three buried pipe modes

4.2.2 水平管出口水溫

水平“spiral”型螺旋埋管的出水溫度303.8K，進(jìn)出水溫差4.3K，埋管與周圍土體的熱量交換能力較強(qiáng)，制冷效果顯著。針對(duì)水平“slinky”型螺旋埋管與直埋管，在制冷季前期，水平“slinky”型螺旋埋管出水溫度稍低于直埋管，因?yàn)樗健皊linky”型螺旋埋管總鋪設(shè)長度基數(shù)大，管路與周圍土壤接觸換熱面積大，因此其換熱效率相對(duì)于直管更加明顯。但隨著熱量的不斷投入，周圍土壤溫度持續(xù)走高，加上水平“slinky”型螺旋埋管本身管路間距較小，管間熱干擾現(xiàn)象明顯，從而導(dǎo)致埋管與周圍土體溫度交換能力急劇下滑，盡管其埋設(shè)管路較長，但其后期的制冷效率不及直管，如圖13所示。

圖13 三種埋管方式出水口溫度變化情況Fig.13 Variation of outlet temperature in three buried pipe modes

5 結(jié)論

（1）地源熱泵系統(tǒng)水平埋管換熱器埋深較淺，其換熱性能易受到土壤水力學(xué)特征的影響。改變地下水位線的埋設(shè)深度，將直接導(dǎo)致了水平埋管管體周圍土壤含水量的變化，從而間接改變了管體周圍土壤的熱物理性質(zhì)。土壤中水分含量越高，則水平管的換熱性能越強(qiáng)，工作效率越高。數(shù)據(jù)表明，在306K 的制冷工況下，土壤含水飽和度100%時(shí)，出口水溫301.7K，埋管進(jìn)出水溫差高達(dá)4.3K；土壤含水飽和度50%時(shí)，出水口溫302.3K，埋管進(jìn)出水溫差3.7K；當(dāng)含水飽和度降低至12%時(shí)，出口水溫升至302.7K，比含水飽和度100%時(shí)升高了23.2%，埋管進(jìn)出口水溫差3.3K。土壤含水飽和度越高其制冷效率越高，同時(shí)基于制冷末期溫度場(chǎng)分布情況可觀測(cè)到，當(dāng)土壤含水飽和度較低時(shí)，土壤熱失衡風(fēng)險(xiǎn)較大，容易引起“熱堆積”現(xiàn)象而使得系統(tǒng)無法正常運(yùn)行。由此可知，土壤的水力學(xué)特性對(duì)水平管的換熱過程有著極其顯著的影響。

（2）在相同的地下水位深度下，由于土壤類型不同同樣會(huì)引起土壤中水分分布的差異性，進(jìn)而影響水平管的換熱性能。本文對(duì)比了水平管換熱器在砂土、壤土、黏土中的換熱性能，其中水平管在壤土中的換熱性能最佳，在黏土中的換熱效果最差。土壤的導(dǎo)熱系數(shù)與比熱容同時(shí)決定了水平管的換熱能力，且增加土壤導(dǎo)熱系數(shù)比增加土壤比熱容更能促進(jìn)水平管與土壤的換熱性能。

（3）水平管的換熱性能受管體自身構(gòu)造的影響較大，連續(xù)模擬制冷90 天后，采用“Slinky”型螺旋埋管時(shí)土壤熱堆積現(xiàn)象嚴(yán)重，采用水平直埋管時(shí)土壤僅有輕微的熱堆積現(xiàn)象，更加利于相同長時(shí)間穩(wěn)定運(yùn)行。基于總換熱量指標(biāo)評(píng)價(jià)，采用“spiral”型螺旋埋管時(shí)末期制冷功率可達(dá)到4005W，采用水平直埋管時(shí)降低至3010W，當(dāng)采用水平“slinky”型螺旋埋管時(shí)僅為2611W。由此可知，“spiral”型螺旋埋管能夠獲得更佳的換熱效率，但缺點(diǎn)在于溫度場(chǎng)分布情況并不樂觀，當(dāng)采用水平直管時(shí)，盡管其溫度場(chǎng)分布情況有利，但由于其土壤利用率較為低下，導(dǎo)致其總換熱量數(shù)值較小難以滿足換熱量要求。