基于熱平衡分析的地埋管地源熱泵換熱方案模擬優化

駱祖江，杜菁菁

?

基于熱平衡分析的地埋管地源熱泵換熱方案模擬優化

駱祖江，杜菁菁

（河海大學地球科學與工程學院，南京 211100）

為了確保丹陽中心城區淺層地熱能可持續開發利用，避免丹陽中心城區地埋管地源熱泵運行期間出現熱堆積問題，基于地下水滲流和熱量運移原理，建立了丹陽中心城區地下水非穩定滲流與熱量運移三維耦合數值模型，結合未來地埋管地源熱泵系統的運行工況，預測丹陽中心城區地下溫度場的熱平衡發展趨勢。在此基礎上，規劃設計了3種優化方案：①按行政區劃Ⅰ～Ⅴ個開發利用分區，調整地埋管間距分別為18，23，17，20，20 m；②地埋管間距5 m，加熱秋季生活用水Ⅰ區157 798 m3/d、Ⅱ區413 235 m3/d、Ⅲ區339 322 m3/d、Ⅳ區261 266 m3/d、Ⅴ區27 6205 m3/d，加熱春季生活用水Ⅰ區473 394 m3/d、Ⅱ區123 9705 m3/d、Ⅲ區1 017 966 m3/d、Ⅳ區783 798 m3/d、Ⅴ區828 615 m3/d；③地埋管間距5 m，增設冷卻塔輔助地埋管換熱孔進行夏季排熱，冷卻塔夏季冷卻溫度Ⅰ區為4.35 ℃、Ⅱ區為6.12 ℃、Ⅲ區為4.87 ℃、Ⅳ區為5.29 ℃、Ⅴ區為4.80 ℃。3種方案均可以有效減緩和避免丹陽中心城區地埋管地源熱泵運行期內的熱堆積問題。地下水非穩定滲流與熱量運移三維耦合數值模型是優化確定淺層地熱能地埋管地源熱泵可持續開發利用方案的有效方法。

熱泵系統；模型；熱平衡；熱堆積；淺層地熱能

0 引 言

隨著城市基礎建設的發展，建筑能耗不斷增加，淺層地熱能憑借其清潔、易開采、儲量大等優點受到人們的重視[1-2]。20世紀40、50年代，歐美國家首先興起了對地源熱泵技術的研究；20世紀70年代，第一次能源危機席卷全球，地源熱泵技術重新引起人們的關注[3]。地埋管地源熱泵（ground-dource heat pump）是地源熱泵系統[4-7]的一個分支，是以地表淺層土壤為低品位熱源，通過由地下埋管組成的地熱換熱器與水循環系統輸入少量高品位能源，實現熱量由低到高的轉變，進而實現城市夏季供冷冬季供熱的需求[8-11]。進入21世紀以來，中國陸續開展了淺層地熱能地埋管地源熱泵的開發利用，目前中國地埋管地源熱泵系統的應用面積已超過2 000萬m2，且每年以超過15%的速度遞增[12]。但由于目前尚處于起步階段，普遍缺乏科學的設計和規劃，地源熱泵系統在應用過程中相繼出現了熱堆積問題，嚴重影響了淺層地熱能的可持續開發利用[13-16]。

本文以丹陽中心城區地埋管地源熱泵為例，通過建立江蘇丹陽中心城區地下水非穩定滲流與熱量運移三維耦合數值模型，結合未來地源熱泵系統的運行工況，模擬預測了丹陽地下溫度場的熱平衡發展趨勢，在此基礎上優化模擬設計了3種不同開發利用方案，即調整地埋管間距[17-19]，保持地埋管間距不變增加春秋季生活用水，和保持地埋管間距不變增設冷卻塔輔助換熱孔排熱[20-22]，有效解決了地源熱泵系統運行過程中的熱堆積問題，確保了淺層地熱能的可持續開發利用[23]。

1 地質概況

丹陽地處江蘇南部，長江南岸，屬于北亞熱帶南部季風氣候區，地表水系發達。全市地勢西北高、東南低，地處低山丘陵與平原過渡帶，地貌以低山丘陵和平原為主，前第四系地層屬揚子地層區。本次研究范圍為丹陽中心城區，詳見圖1，圖中DW1，DW4，DW6為熱響應試驗孔，限制深度為100 m。基巖面埋深為30～90 m，包含山前崗地第四系全新世地層、下蜀組、柏山組，太湖沖積平原第四系全新世如東組、更新世滆湖組、昆山組、啟東組，巖性以黏土、粉砂為主，少有鐵錳質結核。丹陽中心城區缺失第四系上更新統下部的第Ⅰ承壓含水層，第四系中、下更新統第Ⅱ、Ⅲ承壓含水層。潛水含水層接受大氣降水及農田灌溉的補給，由全新統和上更新統上部的褐黃、灰黃色亞黏土和灰、青灰色亞砂土夾粉砂組成，厚度為5～20 m，水位埋深為0.50～4.31 m。

計算區域恒溫帶深度為10～30 m，鉆孔埋深100 m范圍內，初始平均溫度為18.7～19.26 ℃，巖土體的平均導熱系數為1.68～2.07 W/(m·K)。

注：DW1、DW4、DW6為孔號，100為孔深度。

2 地下水非穩定滲流與熱量運移模型

2.1 地下水非穩定滲流與熱量運移概念模型

以丹陽中心城區為計算區域，垂向上按地層巖性結構由上到下進行概化，共4層，依次為厚度為3～20.4 m的全新世地層、厚度為5.9～44.6 m的晚更新世地層、厚度為6.7～29.4 m的中更新世及厚度為2.2～16.3 m的部分早更新世地層和基巖層，計算深度按換熱器長度取100 m。各熱儲層水平方向上均概化為非均質各向同性，各層水平方向和垂直方向的滲透性相差較大，為水平和垂直2個方向上的各向異性，各層之間均發生水力聯系，地下水流態為三維非穩定流。系統頂部為補給邊界接受大氣降雨的補給，排泄邊界為地下水通過系統頂部進行蒸發之處，隔水邊界為系統底部。各熱儲層邊界上水力交換密切且其四周均概化為通用水頭邊界。設定計算范圍內地下水和含水介質骨架的熱動平衡瞬時完成，忽略由于溫度差引起水的密度不同而引起的上下自然對流的影響，并將熱儲層四周概化為第二類熱量邊界條件。

2.2 地下水非穩定滲流與熱量運移數學模型

取坐標軸方向與熱儲層各向異性主滲透系數方向一致，可以建立工作區地下水非穩定滲流與熱量運移三維耦合數值模型[24-30]。

地下水非穩定滲流數學模型：

式中k、k、k為各向異性主方向滲透系數，m/d；為點（,,）在0時刻的水頭值，m；為源匯項，d-1；為時間，d；為計算域；0(,,,0)為點（,,）處的初始水頭值，m；(,,)為第二類邊界上單位面積的補給量；cos(,)，cos(,)，cos(,)為流量邊界外法線方向與坐標軸方向夾角的余弦；為飽和差（自由面上升）或給水度（自由面下降）；2、3為第二類邊界和自由面邊界。

地下水熱量運移數學模型

式中、、為各方向水的熱動力彌散系數，W/(m·K)；c為水的熱容量，J/(m3·K)；為含水介質的熱容量，J/(m3·K)；v、v、v為地下水滲流速度分量，m/d；0(,,)為點(,,)處初始溫度值，K；(,,,)為邊界溫度函數；Q為熱源匯項Q=cW(T?)；T為源匯項的溫度，℃；(,,,)為第二類邊界上已知的熱量或熱流函數（J/(m3·d)）。

地下水流運動方程：

將式式（1）、（2）、（3）耦合，共同構成地下水非穩定滲流與熱量運移三維耦合數值模型。

2.3 模型識別與校正

上述模型采用伽遼金有限元方法進行求解，并采用預處理共軛梯度法求解線性代數方程組。整個計算過程采用FDRTRAN 90語言，按結構化矩陣的設計要求，編制成計算機程序，在PC機加以完成。

平面上將計算域剖分成5 324個三角形單元，每層節點2 764個，垂向上考慮到地層結構和巖性，共剖分為4個地層，5個計算層面，共計21 296個單元，13 820個節點。計算域平面與空間單元剖分見圖2。

模型的校正、識別采用熱響應試驗孔地溫監測數據進行溫度擬合，其中水文地質參數根據土樣滲透系數室內測試統計賦值。根據熱響應測試運行期間所取得的地溫監測資料，反演土體熱物性參數，選取時段2016年3月22日14：00至2016年4月19日5：00作為模型的校正識別時段，共663個應力期，每個應力期一個時間步長。

識別驗證階段采用2016年3月22日14：00各熱儲層實測溫度，作為模型各層的初始溫度值；根據地下水位統測資料，確定2016年3月22日14：00各熱儲層的水位，作為模型各層的初始水位值。

圖2 計算域

通過地溫監測資料，不斷調整參數分區，并進行反演求參，最終將整個計算區分成了30個參數分區，并得到了各分區的參數值，各熱儲層淺層地熱能熱物性參數分區見圖3，各參數分區的參數值見表1。圖4以DW4、DW6為例說明了溫度觀測值和計算值擬合精度，從圖4中可以看出DW4孔溫度觀測值與計算值的差不超過6 ℃，DW6孔溫度觀測值與計算值的差不超過2 ℃。

圖3 第1～4熱儲層淺層地熱能熱物性參數分區

表1 各層位各參數分區表

a. DW4

b. DW6

圖4 DW6孔溫度觀測值與計算值擬合圖

Fig.4 Temperature observation value and calculated value fitting diagram in DW6

3 熱平衡模擬計算

3.1 熱平衡發展趨勢模擬分析

按行政區將整個丹陽中心城區劃分成Ⅰ～Ⅴ個開發利用分區，根據丹陽市實際情況，采取南方典型的冷暖運行模式：冷暖聯供模式：每年的6－9月為制冷期（122 d），10－11月為停采期（61 d），12至次年2月為供暖期（90 d），3－5月為停采期（92 d），此為一個運行周期。模擬假設系統運行10 a，且換熱系統全天運行10 h，Ⅰ～Ⅴ區井間距均為5 m。

我國山區村鎮居民防災知識相對比較缺乏，要通過防災宣傳教育，使人們掌握應對各種地質災害的防災知識，提高自救能力。對泥石流來說，應懂得：不要停留或躲在屋里；露宿時避開陡峭的懸崖和溝壑；逃生時往旁側面跑，跑到開闊地帶或堅固的高地上；山谷中出現巨大沉悶的響聲、溝槽斷流和溝水變渾等現象，預示可能發生泥石流。

在上述工況下，預測地埋管換熱方式淺層地熱能開采運行10 a地下水滲流與熱量運移的變化趨勢。預測時間段為2015年6月初至2025年5月底，每年12個應力期，每個應力期一個時間步長。初始溫度場等值線圖及運行10 a溫度場等值線圖見圖5、圖6，以Ⅲ區DW6孔為例，運行10 a該監測孔處地溫歷時曲線見圖7，各監控點處的溫度變化值見表2。

表2 監控點溫度變化值 Table 2 Temperature change value 計算分區Partition地埋管間距Spacing of buried pipe/m監控點編號Monitory point no.地層溫度增幅Formation temperature increase/℃均值Mean value/℃ Ⅰ區512.776.90 211.03 Ⅱ區5312.2510.24 48.22 Ⅲ區5511.8411.64 611.45 Ⅳ區575.397.44 89.50 Ⅴ區598.386.71 105.03

圖5 2015年6月1日模型第1～4層初始溫度場圖 Fig.5 Natural temperature field map of 1-4 layer in Jun 1, 2015

圖6 2025年5月31日第1～4層預測溫度場圖 Fig.6 Forecast temperature field map of 1-4layer in May 31, 2025

圖7 DW6監測孔處運行10 a的地溫歷時曲線 Fig.7 Earth temperature diachronic curve of DW6 after operating 10 years

由于地埋管換熱方式淺層地熱能的開發，隨著換熱孔夏季放熱冬季取熱，溫度在局部范圍內將發生較大變化，熱量逐年堆積。由此將對周圍水環境產生一系列的影響，也將造成地埋管換熱系統在第一個運行季節的中期或是后期換熱效果減弱，換熱系統能效降低甚至無法正常運行。

3.2 地埋管換熱方式優化模擬分析

由運行方案預測結果可知，由于地埋管間距過近，淺層地熱能開發引起熱量堆積在地下難以消耗，隨時間推移，熱堆積越來越嚴重。

為實現淺層地熱能的可持續開發利用，應用所建立的研究區地下水非穩定滲流與熱量運移三維耦合數值模型，分3種方案進行優化模擬分析，討論3種方案減緩熱量堆積的程度：方案1當各計算分區單孔年內取放熱量確定時，運行工況不變，針對熱堆積問題，模擬計算各分區淺層地熱能開發利用的最小地埋管間距；方案2當各計算分區地埋管間距為5 m、單孔年內取放熱量確定時，適當改變運行工況，針對熱堆積問題，模擬計算春秋季需加熱生活用水的水量；方案3當各計算分區地埋管間距為5 m、單孔年內取放熱量確定時，增加冷卻塔輔助設備，針對熱堆積問題，模擬計算夏季循環介質經冷卻塔需冷卻的溫度。

3.2.1 方案1模擬結果與分析

研究數據選取科睿唯安公司的德溫特數據庫；檢索時間為2018 年7月；檢索方式利用“關鍵詞”+“IPC分類號”進行組合檢索；檢索范圍為1963年至今與自動駕駛技術相關的所有專利數據；檢索初步結果為相關專利共計22 184條，經查重、清洗得到最終有效專利22 177條。

趙新宇等(2013)利用吉林大學公眾幸福指數課題組關于2012年中國公眾主觀幸福感問卷調查數據，運用有序概率模型考察了絕對收入、相對收入和預期對公眾主觀幸福感的影響。發現我國存在“幸福悖論”現象，相對收入對公眾主觀幸福感有顯著促進作用，其效果強于絕對收入；預期對于中、低收入群體的主觀幸福感具有顯著正向作用。

方案1通過模擬計算，將Ⅰ～Ⅴ區的最小地埋管間距改變為18、23、17、20、20 m，各監控點10 a后100 m以上淺地層平均溫度增幅及可開采資源量見表3，DW6孔處運行10 a的溫度歷時曲線見圖8。

表3 方案1監控點溫度變化及可開采資源量 Table 3 Temperature change value and recoverable resources in scheme 1 計算分區Partition地埋管間距Spacing of buried pipe/m監控點編號Monitory point No.地層溫度增幅Formation temperature increase/℃均值Mean value/℃夏季總排熱量Total summer heat /kW冬季總取熱量Total winter heat/kW Ⅰ區1810.490.4837 355.2933 187.23 20.47 Ⅱ區2330.480.4445 829.0134 421.35 40.40 Ⅲ區1750.460.4755 537.8446 883.83 60.49 Ⅳ區2070.410.4354 658.9544 335.95 80.45 Ⅴ區2090.460.4261 665.9252 484.28 100.38

圖8 方案1下DW6孔運行10 a的溫度歷時曲線 Fig.8 Earth temperature diachronic curve of DW6 in scheme 1 after operating 10 years

由計算結果可知，地埋管淺層地熱能開采10 a后，由于地埋管間距增加，年內地層溫度升幅減小；由監控點多年歷時曲線圖可知，制冷期內監控點處土壤溫度逐漸上升，供熱期內監控點處土壤溫度逐漸下降，隨著淺層地熱能的開發利用，土壤溫度呈緩慢上升的趨勢，且10 a內各點處升溫均不超過0.5 ℃，熱堆積問題得到有效緩解。

3.2.2 方案2模擬結果與分析

已知Ⅰ～Ⅴ區需加熱生活用水的總面積為173.23 km2，對各計算分區單孔換熱量分別按5%損失計算，將生活用水從15 ℃加熱到50 ℃，獲得各區加熱生活用水量詳見表4。通過模擬計算，保持地埋管地埋管間距為5 m，增加春秋季加熱生活用水1 447 826.19 m3/d的運行方案下，熱堆積問題得到緩解。各監控點10 a后100 m以上淺地層平均溫度增幅見表5，淺層地熱能可開采資源量見表6，DW6孔處運行10 a的溫度歷時曲線見圖9。

表4 方案2春秋季各區加熱生活用水量 Table 4 Heating water consumption in each district in spring and autumn in scheme 2 運行時間Running time加熱水量Heating water consumption /(m3·d-1) Ⅰ區Ⅱ區Ⅲ區Ⅳ區Ⅴ區合計Total 秋季Autumn11月157 798413 235339 322261 266276 2051 447 826 春季Spring3月157 798413 235339 322261 266276 2051 447 826 4月157 798413 235339 322261 266276 2051 447 826 5月157 798413 235339 322261 266276 2051 447 826

表5 方案2監控點溫度變化值 Table 5 Temperature change values of observe points in sheme 2 計算分區Partition地埋管間距Spacing of buried pipe/m監控點編號Monitory point No.地層溫度增幅Formation temperature increase/℃均值Mean value/℃ Ⅰ區510.140.26 20.39 Ⅱ區530.360.31 40.25 Ⅲ區550.330.29 60.25 Ⅳ區570.220.27 80.32 Ⅴ區590.280.23 100.17

表6 方案2淺層地熱能可開采資源量 Table 6 Shallow geothermal energy recoverable resources in scheme 2 計算分區Partition面積Proportion/km2夏季總排熱量Total summer heat discharge/kW秋季總取熱量Total autumn heat absorption/kW冬季總取熱量Total winter heat absorption/ kW春季總取熱量Total spring heat absorption/ kW Ⅰ區24.34484 311.6760 255.45430 174.1260 255.45 Ⅱ區41.78809 232.13159 761.57605 509.36159 761.57 Ⅲ區36.33889 268.92130 489.79750 238.68130 489.79 Ⅳ區36.53631 064.91100 121.91512 168.62100 121.91 Ⅴ區34.25712 711.50105 606.96606 147.68105 606.96 合計Total173.233 526 589.13556 235.692 904 238.45556 235.69

圖9 DW6孔運行10 a的溫度歷時曲線 Fig.9 Earth temperature diachronic curve of DW6 after operating 10 years

從計算結果可知，地埋管換熱方式淺層地熱能開采10 a后，由于增加春秋季加熱生活用水，年內地層溫度幾乎無增幅；由監控點溫度多年歷時曲線圖9可知，夏季制冷期內監控點處土壤溫度逐年上升，秋季、冬季、春季供熱期內監控點處土壤溫度逐漸下降，隨著淺層地熱能的開發利用，土壤溫度呈緩慢上升的趨勢，且10 a內各點處的溫度上升平均不超過0.5 ℃，熱堆積問題得到有效緩解。

測試設備使用了超寬帶厘米波設備，具有成本低、輻射小和便于推廣的特點.實驗結果驗證可以準確識別介電常數差別較大的液體危險品，如汽油和水.算法為微波檢測液體危險品提供了一個新的研究方向.算法存在對介電常數相似的樣品識別率低、散射曲線平滑程度特征散布較大這兩個難點，今后將從分析樣本的分布、改進訓練算法、隔離外部干擾源等方面做進一步研究.

本文選取國際原油市場比較有代表性的WTI原油現貨價格與1月、2月、3月和4月期的期貨價格的月度數據為研究對象，為消除時間序列可能存在的異方差，對現貨價格和期貨價格取自然對數，分別用lnS、lnF(1)、lnF(2)、lnF(3)和lnF(4)表示，數據來源于國際能源信息署，所使用的樣本數據的時間跨度為1986年1月到2013年5月，共329個觀測。對相應合約在某個給定日并無交易的情形，遵循Alquist和Kilian(2010)的方法，將最接近月度最后一個交易日的相應期貨價格作為其月度值。表1給出WTI現貨與期貨的描述性統計量和單位根檢驗結果。

3.2.3 方案3模擬結果與分析

將冷卻塔與換熱孔共同作為冷卻源供機組運行，經過冷卻塔冷卻的水，進入地埋管換熱器，降低了換熱器的冷凝溫度，實現兩者的互補，提高整個系統的運行效率。保持地埋管間距5 m、冷卻塔與地埋管換熱孔數量一致，冷卻塔在6～8月開機，并且與地埋管換熱器每日同時運行10 h的工況下，忽略換熱器與地下埋管之間管路的熱量損失。根據巖土體現場熱響應試驗及地溫監測資料，夏季的換熱溫差一般為13.24～14.27 ℃，冬季的換熱溫差一般為10.73～11.76 ℃；根據《江蘇省地源熱泵工程技術規程》（DGJ32/TJ89－2009），地源熱泵系統的制冷工況：空調側出水溫度為7～12 ℃，冷卻水溫差為5 ℃，因此Ⅰ～Ⅴ區夏季各區地埋管及冷卻塔冷卻水溫詳見表7。各監控點10 a后100 m以上淺地層平均溫度增幅及可開采資源量見表8，DW6孔處運行10 a的溫度歷時曲線見圖10。

表7 方案3下Ⅰ～Ⅴ區地埋管及冷卻塔冷卻水溫 Table7 Cooling water temperature of buried pipe and cooling tower in Ⅰ～Ⅴarea in scheme 3 分區編號Number分區面積Proportion/km2換熱溫差夏季Summer heat exchange temperature difference/℃換熱溫差冬季Winter heat exchange temperature difference /℃冷卻塔夏季冷卻溫度Cooling tower summer cooling temperature/℃ Ⅰ區24.3413.2411.764.35 Ⅱ區41.7814.2710.736.12 Ⅲ區36.3313.5611.444.87 Ⅳ區36.5313.8011.205.29 Ⅴ區34.2513.5111.494.80 合計Total173.23———

圖10 方案3下DW6孔運行10 a溫度歷時曲線 Fig.10 Earth temperature diachronic curve of DW6 after operating 10 years in scheme 3

表8 方案3監控點溫度變化及可開采資源量 Table 8 Temperature change value and recoverable resources in scheme 3 計算分區Partition地埋管間距Spacing of buried pipe/m監控點編號Monitory point no.地層溫度增幅Formation temperature increase/℃均值Mean value/℃夏季總排熱量Total summer heat discharge/kW冷卻塔夏季總散熱量Summer total heat dissipation of cooling tower/kW冬季總取熱量Total winter heat absorption /kW Ⅰ區510.230.32484 311.67159 200.40430 237.27 20.41 Ⅱ區530.320.26809 232.13348 168.73607 176.35 40.20 Ⅲ區550.310.31889 268.92319 667.45750 197.07 60.31 Ⅳ區570.220.26631 064.91241 990.93512 153.66 80.30 Ⅴ區590.310.24712 711.50253 016.90606 176.37 100.18

由計算結果可知，地埋管換熱方式淺層地熱能開采10 a后，由于采用冷卻塔與地埋管混合式換熱系統，彌補了地埋管單獨排熱導致熱堆積的缺陷，保證溫度有較好的恢復，10 a內地層溫度幾乎無增幅；由監控點溫度多年歷時曲線圖10可知，夏季制冷期內監控點處土壤溫度逐漸上升，冬季供暖期內監控點處土壤溫度逐漸下降，隨著淺層地熱能的開發利用，土壤溫度呈緩慢上升的趨勢，且10 a內各點處溫度上升平均不超過0.5 ℃，熱堆積問題得到有效緩解。

這類現象還有很多，比如 “閉月羞花”中的“閉”和“羞”、“打草驚蛇”中的“驚”、“楚楚動人”中的“動” 、“驚心動魄”中的“驚”和“動”、“感人肺腑”中的“感”、“翻江倒海”中的“翻”和“倒”、“除舊布新”中的“除”和“布”都用作使動。

3.3 結果分析

從上述模擬計算結果可以看出，研究區如不加輔助設施，采用5 m的地埋管間距開采淺層地熱能會導致熱堆積問題，嚴重影響淺層地熱能的可持續開發利用，必須增大地埋管的間距，Ⅰ～Ⅴ開發利用分區的地埋管間距分別要增至18，23，17，20，20 m，才能避免熱堆積問題，實現淺層地熱能的可持續開發利用。

但在實際開發利用過程中會受到開發利用場地條件的限制，無法采取大間距的地埋管方式來實現淺層地熱能的可持續開發利用。目前江蘇南部地區常用的地埋管間距一般為5 m左右，因此研究區要實行淺層地熱能的可持續開發利用必須在地埋管間距5 m的基礎上，增加輔助設施，如方案2、3采用加熱春秋季生活用水或增設地面冷卻塔的辦法來消除熱堆積問題，實現淺層地熱能的可持續開發利用。

4 結 論

1）基于地下水滲流與熱泵運移原理建立的地下水非穩定滲流和熱量運移模型，可以準確模擬預測地埋管地源熱泵三維耦合數值模型系統不同工況運移過程中的熱平衡發展趨勢，確定出淺層地熱能的可持續性開發利用方案。

2）進行地埋管換熱方式淺層地熱能開發時，換熱系統夏季開機122 d、冬季開機90 d、每天運行10 h，選取Ⅰ～Ⅴ區地埋管間距分別為18，23，17，20，20 m，可以有效緩解地溫場熱堆積問題。

2）裂縫寬度大于0.3mm的環向裂縫，屬結構受力裂縫。但隧道結構為拱形截面受力，此類裂縫不影響結構的安全性，采用環氧樹脂漿抹面封閉，滲漏嚴重的，先排水后封閉處理即可。

3）進行地埋管換熱方式淺層地熱能開發時，保持地埋管間距為5 m，春季、秋季分別對生活用水加熱92、31 d，每天運行10 h，將生活用水從15 ℃加熱到50 ℃，可以有效緩解地溫場熱堆積問題。

4）進行淺層地熱能開發時，地埋管換熱方式為保持地埋管間距5 m，換熱系統夏季開機122 d、冬季開機90 d、每天運行10 h，冷卻塔夏季開機122 d、每天運行10 h，可以有效緩解地溫場熱堆積問題。

[參 考 文 獻]

[1] 張甫仁，彭清元，朱方圓，等. 重慶主城區淺層地溫能資源量評價研究[J]. 中國地質，2013，40(3)：974－980.

Zhang Furen, Peng Qingyuan, Zhu Fangyuan, et al. The evaluation of shallow geothermal energy resources in Chongqing[J]. Geology in China, 2013, 40(3): 974－980. (in Chinese with English abstract)

[2] 柴立龍，馬承偉，張義. 北京地區溫室地源熱泵供暖能耗及經濟性分析[J]. 農業工程學報，2010，26(3)：249－254.

Chai Lilong, Ma Chengwei, Zhang Yi, et al. Energy consumption and economic analysis of ground source heat pump used in greenhouse in Beijing[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010, 26(3): 249－254. (in Chinese with English abstract)

[3] 駱祖江，李偉，王琰，等. 地下水源熱泵系統熱平衡模擬三維數值模型[J]. 農業工程學報，2014，30(2)：198－204.

Luo Zujiang, Li Wei, Wang Yan, et al. Three-dimensional numerical model for heat balance simulation of groundwater heat pump[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(2): 198－204. (in Chinese with English abstract)

[4] 胡金強. 地源熱泵系統熱平衡分析及其在大型公共建筑中的應用[J]. 制冷技術，2015，35(2)：63－67.

Hu Jinqiang. Analysis of heat balance in ground source heat pump system and its application in large scale public buildings[J]. Chinese Journal of Refrigeration Technology, 2015, 35(2): 63－67. (in Chinese with English abstract)

[5] John W Lund, Derek H Freeston, Tonya L Boyd. Direct utilization of geothermal energy 2010 worldwide review[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2011, 40(3): 159－180. (in Chinese with English abstract)

[6] 田良河，閆震鵬，劉新號. 鄭州市地下水源熱泵適宜區淺層地熱能資源量評價[J]. 城市地質，2011，6(3)：12－16.

Tian Lianghe, Yan Zhenpeng, Liu Xinhao. The shallow geothermal energy resource evaluation on the suitable area of groundwater source heat pump in Zhengzhou city[J]. City Geology, 2011, 6(3): 12－16. (in Chinese with English abstract)

[7] Pike Christopher, Whitney Erin. Heat pump technology: An Alaska case study[J].Journal of Renewable and Sustainable Energy, 2017, 9(6)

[8] 岳高凡，趙紫威，王婉麗，等. 地源熱泵系統對土壤溫度場影響的數值模擬研究——以北京某場地為例[J]. 中國科技論文，2017，9(12)：1049－1053.

Yue Gaofan, Zhao Ziwei, Wang Wanli, et al. Numerical simulation study on soil temperature field influenced by GSHP: A case study in Beijing[J]. China Science Paper, 2017, 9(12): 1049－1053. (in Chinese with English abstract)

[9] 徐偉，劉志堅. 中國地源熱泵技術發展與展望[J]. 建筑科學，2013，29 (10)：26－33.

Xu Wei, Liu Zhijian. Development and prospect of ground source heat pump technology in China[J]. Building Science, 2013, 29 (10): 26－33. (in Chinese with English abstract)

[10] 駱祖江，徐曉，常曉軍. 基于Monte-Carlo參數模擬的體積法評價淺層地熱能資源[J]. 勘查科學技術，2017(6)：15－22.

Luo Zujiang, Xu Xiao, Chang Xiaojun. Evaluation of shallowgeothermal energy resources based on Monte-Carlo parameter simulation in volume method[J]. Site Investigation Science and Technology, 2017(6): 15－22. (in Chinese with English abstract)

[11] 岳麗燕，孟令軍，趙蘇民，等. 天津市淺層地熱能存在的熱堆積問題及解決方法探討[J]. 地質調查與研究，2017，40 (1)：76－80.

Yue Liyan, Meng Lingjun, Zhao Sumin, et al. The cold & heat accumulation problem and solution methods of the shallow geothermal resource in Tianjin[J]. Geological Survey and Research, 2017, 40(1): 76－80. (in Chinese with English abstract)

[12] 方慧，楊其長，孫驥. 地源熱泵—地板散熱系統在溫室冬季供暖中的應用[J]. 農業工程學報，2008，24(12)：145－149.

Fang Hui, Yang Qichang, Sun Ji. Application of ground-source heat pump and floor heating system to greenhouse heating in winter[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2008, 24(12): 145－149. (in Chinese with English abstract)

[13] 龍西亭，袁瑞強，皮建高，等. 長沙淺層地熱能資源調查與評價[J]. 自然資源學報，2016，31(1)：163－176.

Long Xiting, Yuan Ruiqiang, Pi Jiangao, et al. Survey and evaluation of shallow geothermal energy in Changsha[J]. Journal of Natural Resources, 2016, 31(1): 163－176. (in Chinese with English abstract)

[14] Barbier E. Geothermal energy technology and current status: Aan overview[J]. Renewable & Sustainable Energy Reviews, 2002, 6(1/2): 3－65.

[15] Hettiarachchia H D Madhawa, Golubovica Mihajlo, Worek William M, et al. Optimum design criteria for an organic Rankine cycle using low-temperature geothermal heat source[J]. Energy, 2007, 32(9): 1698－1706.

[16] 張遠東，萬育生. 我國地下水源熱泵應用現狀和監管措施探討[J]. 水資源管理，2009(21)：39－41.

Zhang Yuandong, Wan Yusheng. Current status of groundwater source heat pump utilization in China and supervisory suggestions[J]. China Water Resources, 2009(21): 39－41. (in Chinese with English abstract)

[17] 胡繼華，張延軍，于子望，等. 水源熱泵系統中地下水流貫通及其對溫度場的影響[J]. 吉林大學學報：地球科學版，2008，38(6)：992－998.

Hu Jihua, Zhang Yanjun, Yu Ziwang, et al. Groundwater flow transfixion of groundwater source heatpump system and its influence on temperature field[J]. Journal of Jilin University: Earth Science Edition, 2008, 38(6): 992－998. (in Chinese with English abstract)

[18] 徐貴來，劉紅衛，張晴，等. 地下水源熱泵系統運行期間巖土層溫度變化規律研究[J]. 資源環境與工程，2010，24(2)：180－184.

Xu Guilai, Liu Hongwei, Zhang Qing, et al. Geotemperature changes in rock and soil stratum during groundwater source heat pump system operation[J]. Resources Environment and Engineering, 2010, 24(2): 180－184. (in Chinese with English abstract)

[19] Hepbasli A. Energy and exergy analysis of a ground source (geothermal) heat pump system[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2004, 45(5): 737－753. (in Chinese with English abstract)

[20] Lim K B, Lee S H, Lee C H. Application of ground-heat- source water thermal-storage with a heat-pump cooling and heating system[J]. Proceedings of the Institution of Mechanical E-journal of Process Mechanical Engineering, 2010, 224(E1): 73－78.

[21] Paksoy H O, Anderson O, Abaci S, et al. Heating and cooling of a hospital using solar energy coupled with seasonal thermal energy storage in an aquifer[J]. Renewable Energy, 2000, 19(1): 117－122.

[22] Camdali U, Tuncel E. An Economic Analysis of Horizontal Ground Source Heat Pumps (GSHPs) For Use in Heating and Cooling in Bolu, Turkey[J].Energy Sources Part B-Economics Planning and Policy, 2013, 8(3): 290－303.

[23] Kavanaugh S. Simple, efficient, reliable - Ground source heat pumps[J]. Ashrae journal-American Society of Heating Refrigerating and Air-conditioning Engineers, 1998, 40(10): 31.

[24] 楊衛波，施明恒. 地源熱泵中U型埋管傳熱過程的數值模擬[J]. 東南大學學報：自然科學版，2007，37(1)：78－83.

Yang Weibo, Shi Mingheng. Numerical simulation on heat transfer process in U-bend ground heat exchanger of ground source heat pump[J]. Journal of Southeast University: Natural Science Edition, 2007, 37(1): 78－83. (in Chinese with English abstract)

[25] Luo Z, Wang Y, Zhou S, et al. Simulation and prediction of conditions for effective development of shallow geothermal energy[J]. Applied Thermal Engineering, 2015, 91: 370－376.

[26] Ren Xiaohong, Sun Chunwu, Hu Yanhui. A three- dimensional numerical model and experimental research for U-tube ground heat exchanger in heating operation[J]. Journal of Chongqing Jianzhu University, 2004, 26(5): 90－95.

[27] 周彥章，周志芳，傅志敏. 地源熱泵系統可持續運行地下水熱量運移模擬[J]. 東南大學學報：自然科學版，2009，39(增刊)：297－303.

Zhou Yanzhang, Zhou Zhifang, Fu Zhimin. Simulation of thermal transport in aquifer for the sustainable running of GWHP system[J]. Journal of Southeast University: Natural Science Edition, 2009, 39(Supp.): 297－303. (in Chinese with English abstract)

[28] 朱小波，李幽錚，周彥章. 地源熱泵系統冷負荷設計方案變化地下水熱量運移模擬[J]. 勘察科學技術，2010(5)：13－17.

Zhu Xiaobo, Li Youzheng, Zhou Yanzhang. Simulation of groundwater thermal migration for GWHP system[J]. Site Investigation Science and Technology, 2010(5): 13－17. (in Chinese with English abstract)

[29] 靳萍，邵景力，李長青，等. 基于T-PROGS的地下水三維數值模擬及應用[J]. 水文地質工程地質，2009，36(4)：21－26.

Jin Ping, Shao Jingli, Li Changqing, et al. Application of T-PROGS to a 3-D numerical simulation of groundwater flow[J]. Hydrogeology & Engineering Geology, 2009, 36(4): 21－26. (in Chinese with English abstract)

[30] Berthier L, Garrahan J P. Numerical study of a fragile three-dimensional kinetically constrained model[J]. Journal of Physical Chemistry B, 2005, 109(14): 6916.

Heat exchange scheme simulation optimization for ground source heat pump system with buried pipes by thermal equilibrium analysis

Luo Zujiang, Du Jingjing

(School of Earth Science and Engineering, Hohai University, Nanjing 211100, China)

Abstract: Since the 21st century, the development and utilization of ground source heat pump in shallow geothermal energy is carried out in China successively. But it is still in its infancy, and there is a general lack of scientific design and planning, the problem of heat accumulation occurred in the process of application of ground source heat pump system, seriously affecting the sustainable development and utilization of shallow geothermal energy. In order to guarantee shallow geothermal energy sustainable development and utilization and avoid the problem of heat accumulation during the running period of ground-source heat pump in the downtown of Danyang, a three-dimensional coupled numerical model of groundwater unsteady seepage and heat transfer was established and applied to the demonstration project of the buried pipe heat pump system in Danyang, Jiangsu Province. The model was based on the principle of groundwater seepage and heat transport, combined with future ground-source heat pump system operating conditions, forecasting the thermal equilibrium development trend of underground temperature field. Three schemes were designed on this basis. Scheme 1 is to adjust heat transfer hole spacing. By simulation and calculation, we changed the minimum heat transfer 5 m hole spacing to 18, 23,17, 20, 20 m inⅠ-Ⅴarea respectively, and results showed that due to the increase of heat transfer hole spacing through mining buried pipe shallow geothermal energy for ten years, the rise in formation temperature decreased. Scheme 2 is to keep the heat exchange hole spacing constant, and heating spring and autumn living water. The total area of heating domestic water in study area was 173.23 km2, single hole heat exchange was calculated with 5% heat loss and heating water from 15 to 50 ℃. By the simulation and calculation, keeping 5 m transfer hole spacing, the heating water amount reach 1 447 826.19 m3/d in spring and autumn with the development and utilization of shallow geothermal energy. Mining 10 years, due to increasing heating water in spring and autumn, there is little increase in formation temperature during the year. The temperature of soil in winter, autumn and spring increases slowly, the problem of heat accumulation is alleviated. Scheme 3 is to keep the heat exchange hole spacing constant and add the cooling tower auxiliary heat exchange hole for summer heat discharge. The cooling tower is in accordance with the number of heat exchange holes in buried pipes and keeping 5 m transfer hole spacing, running 10 hours a day, let water through cooling towers of different cooling temperatures. Mining 10 years, the mixed heat transfer system with cooling tower and buried pipe can eliminate the heat accumulation caused by separate buried pipe system, ensure a good temperature recovery. Due to the limit of ground source heat pump engineering site, Scheme 1 can’t be used, heating water in spring and autumn or adding cooling tower auxiliary heat exchange hole for summer heat discharge can be adopted in the actual project. The three-dimensional coupled numerical model of groundwater unsteady seepage and heat pump is an effective way to optimize and determine the sustainable development and utilization plan of ground source heat pump for shallow geothermal energy.

Keywords: heat pump systems; models; heat balance; heat accumulation; shallow geothermal energy

駱祖江，杜菁菁. 基于熱平衡分析的地埋管地源熱泵換熱方案模擬優化[J]. 農業工程學報，2018，34(13)：246－254. doi：10.11975/j.issn.1002-6819.2018.13.030 http://www.tcsae.org

Luo Zujiang, Du Jingjing. Heat exchange scheme simulation optimization for ground source heat pump system with buried pipes by thermal equilibrium analysis[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(13): 246－254. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2018.13.030 http://www.tcsae.org

收稿日期：2018-02-27

修訂日期：2018-05-31

基金項目：中國地質調查局江蘇省合作項目（12120114023101）

作者簡介：駱祖江，博士，教授，博士生導師，主要從事水文地質等方面研究與教學工作。Email：luozujiang@ sina.com

doi：10.11975/j.issn.1002-6819.2018.13.030

中圖分類號：TU111; TK52

文獻標志碼：A

文章編號：1002-6819(2018)-13-0246-09