郭占全 苑中顯 張 琦

跨季節(jié)土壤蓄熱太陽能數(shù)值模擬研究

郭占全 苑中顯 張 琦

（北京工業(yè)大學環(huán)境與能源工程學院 北京 100124）

為減少太陽能利用所受的季節(jié)性限制，以太陽能跨季節(jié)蓄熱為背景，土壤為蓄熱介質，利用ANSYS軟件對水平地埋管與土壤的非穩(wěn)態(tài)換熱過程進行三維數(shù)值模擬。在流速、管徑、土壤初始溫度以及運行時間等條件相同的情況下，分別討論埋管深度、管內流體溫度、土壤熱物性參數(shù)以及埋管間距對土壤蓄熱效果的影響。結果表明：對于埋深為1.5m、2.0m、2.5m這三種工況，埋深2.5m的土壤蓄熱量最大，埋深1.5m的蓄熱量最小。對于流體溫度為60℃、70℃、80℃這三種工況，溫度為80℃的土壤蓄熱量最大，60℃的蓄熱量最小。對于沙土、粘土、褐土和花崗巖這四種土壤，褐土是最有利于蓄熱的，沙土是最不利的。對于埋管間距為2.8m、3.2m、3.6m、4m這四種工況，土壤蓄熱量隨埋管間距的增大而增大，管間距為4m的土壤蓄熱量最大，2.8m的蓄熱量最小。

太陽能；水平地埋管；數(shù)值模擬；土壤蓄熱

0 引言

隨著能源危機的出現(xiàn)，人們對新型能源的需要更加迫切，太陽能作為一種新型能源受到越來越多的關注，有關太陽能的研究利用范圍也在日益擴大。1956年，美國伊利諾伊大學的Penrod等人首次提出使用太陽能集熱器與地埋管換熱器協(xié)同運行的想法[1]，并提出了太陽能集熱器和地埋管換熱器的設計方法，此后，各國的學者相繼開始了對太陽能跨季節(jié)蓄熱系統(tǒng)的研究。太陽能跨季節(jié)蓄熱系統(tǒng)能有效克服太陽能利用易受時間、氣候等影響的缺點，可以將太陽資源豐富季節(jié)的熱量儲存到太陽資源不足的季節(jié)使用。我國太陽能資源豐富,年平均日照時數(shù)在2200h以上，輻射值為1500～2300kWh/m2，氣象條件適合太陽能的推廣[2]。將太陽能作為冬季供暖熱源有其獨特優(yōu)點。對于這種系統(tǒng)應采用何種蓄熱方式以及如何提高蓄熱效率成為目前研究的一個重點。

對于土壤地埋管蓄熱而言，由于其成本低，就地利用土壤，適用范圍廣泛，可大規(guī)模適用等優(yōu)點，通常被人們大量采用。但是，國內外關于垂直地埋管蓄熱的研究較多，而對于水平地埋管的研究相對較少。水平地埋管蓄熱的研究主要集中在以下幾個方面：水平地埋管與土壤耦合換熱的傳熱模型，管內流體的溫度，水平地埋管的鋪設形式，水平地埋管的埋深，土壤的熱物性參數(shù)以及回填材料的分析等對蓄熱的影響。

本文建立了地下埋管換熱器與周圍土壤的計算流體力學數(shù)值模型。通過數(shù)值模擬分析地埋管的埋深、管內流體溫度、土壤的熱物性參數(shù)以及地埋管間距對地埋管周圍土壤溫度分布的影響，其次，由土壤溫度的變化進一步得到土壤蓄熱量的大小，最后討論以上因素對太陽能跨季節(jié)土壤蓄熱效果的影響。

1 物理模型的建立

圖1 太陽能儲熱系統(tǒng)示意圖

圖1為太陽能蓄熱系統(tǒng)的示意圖。在太陽光的照射下，低溫水被太陽能集熱器加熱后成為高溫水。由于高溫水的溫度遠高于土壤的溫度，因此高溫水在水平埋管內流動與土壤進行換熱，土壤的溫度不斷升高，從而太陽熱能被蓄集在地下土壤。

2 水平地埋管與土壤耦合換熱的數(shù)學模型

地埋管與土壤之間的換熱是一個非穩(wěn)態(tài)的過程，其傳熱過程復雜且影響因素眾多。這些影響因素包括管內水的流動速度、材料的物性參數(shù)、土壤的熱物性參數(shù)、回填密實度以及地下水分遷移等問題。下面首先討論單根地埋管工況。圖2為單根地埋管的數(shù)值分析模型。

圖2 水平地埋管蓄熱體模型圖

為了便于計算分析，本文做以下簡化[3,4]：

（1）土壤為各項同性、均質的剛性含濕多孔介質體，在整個過程中土壤的熱物性參數(shù)保持不變；

（2）不考慮直管管壁與回填材料、回填材料與土壤之間的接觸熱阻；

（3）忽略土壤的熱濕遷移耦合作用的影響；

（4）管內流體為不可壓縮流體；

在上述簡化條件下的控制方程如下：

連續(xù)性方程：

動量守恒方程：

能量守恒方程：

3 數(shù)值求解方法

3.1 網(wǎng)格劃分

本文以北京地區(qū)土壤為模擬對象，通過查閱文獻得到當?shù)赝寥赖臒釋省⒚芏?、比定壓熱容、熱擴散率以及不同深度土壤的大致溫度。表1為北京地區(qū)部分土壤及室外空氣主要物性表。

表1 各土壤及空氣的物性參數(shù)

模型中整個蓄熱體由土壤和埋管兩部分構成其網(wǎng)格劃分，見圖3，全部采用結構性網(wǎng)格劃分。由于地埋管尺寸相對于周圍土壤蓄熱體的尺寸小很多，為平穩(wěn)過渡，需要在兩者交界處進行網(wǎng)格加密處理。遠離埋管的土壤受流體加熱的影響較小，那里的網(wǎng)格劃分相對稀疏。

此外，有研究表明，水平地埋管換熱器的管長存在最大長度[5]，一旦超過該長度極限，長度的再增加對蓄熱性能的改善作用將變小。一般情況下，地埋管長度大約在100m范圍內。本文取管長為40m，采用PE管，管直徑為40mm，壁厚3mm。

圖3 蓄熱體三維網(wǎng)格圖及局部放大圖

3.2 求解方法

本文設置流體的流速為0.3m/s，通過計算流體的Re數(shù)約為47000，遠遠大于2300。采用ANSYS FLUENT進行求解，管內傳熱介質的流動模型選擇?湍流模型，采用非穩(wěn)態(tài)計算方式，動量方程與能量方程耦合求解。

3.3 計算條件

北京市供暖時間為11月15日至3月15日，運行時間為整個非供暖季，即3月16號到11月15號。選擇3月16日0點時刻作為計算初始時間，土壤初始溫度[6]按照下式計算：

式中：為土壤深度，m；為土壤初始溫度，K。

由于地表面長時間暴露在隨時間不斷變化的環(huán)境條件下，地表面與大氣之間的對流換熱過程受到室外溫度、相對濕度、風速、風向以及太陽輻射等各種因素的影響，是一個復雜的動態(tài)過程。綜合上述因素，張文武[7]等在研究風速對土壤溫度分布影響時，指出土壤表面的對流換熱系數(shù)與風速變化之間的統(tǒng)計關系可以用下式表述。

其中，的單位為W/m2·K，為當?shù)販y點表面的風速，m/s。

本文取北京地區(qū)地表年均風速為2.3m/s[8]，經(jīng)過式（7）計算，得到地表面與大氣的年平均對流換熱系數(shù)20.9W/m2·K。

此外，為了使數(shù)值模擬的結果更接近實際，本文考慮了北京地區(qū)氣溫的年變化以及太陽日照時長的年變化對土壤蓄熱效果的影響。利用FLUENT處理流動或傳熱問題時，模型的速度和溫度等邊界條件是不隨時間變化的，難以實現(xiàn)動態(tài)模擬。而本文中地埋管的進口流速是動態(tài)變化的，白天有太陽光照時，管內流體速度為0.3m/s，對土壤進行流動蓄熱。夜間無太陽光照時，管內流體速度為0，不對土壤進行蓄熱。此外，地表的氣溫一年四季內也是動態(tài)變化的。因此，需要借助FLUENT擴展功能中的用戶自定義函數(shù)（User Defined Function）來解決。UDF是用戶使用C語言編寫的擴展FLUENT程序代碼，可加載到FLUENT求解器中，實現(xiàn)邊界條件、物性參數(shù)和源項等條件的動態(tài)輸入。

北京地區(qū)氣溫的年變化用余弦函數(shù)表示：

北京地區(qū)的一天中的晝長可用下式表示：

其中太陽赤緯可由Cooper方程近似的進行計算：

計算機在進行模擬計算時得到的流體入口處平均速度隨時間的變化如圖4。速度隨時間呈現(xiàn)出周期性變化。

本文采用的數(shù)值算法是基于原始變量的SIMPLE方法，壓力項采用Second Order方案離散，動量方程和能量方程采用二階迎風格式進行離 散[9]。土壤與埋管內水體采取耦合換熱方式，避免了將壁面設置為恒壁溫或恒熱流邊界條件的問題。

圖4 模擬計算時流體入口平均速度的變化

4 模擬結果及分析

4.1 不同埋管深度對土壤蓄熱效果的影響

在工程實踐中，水平埋管換熱器的埋設深度一般取1.5～3.0m，處于大氣影響范圍內[10]，但無論任何情況均應埋設在當?shù)貎鐾翆右韵?。北京地區(qū)的最大凍土深度為0.85m[11]。為了研究埋管深度對土壤蓄熱效果的影響，本算例取流體的入口溫度為353K，土壤為褐土，對埋管深度為1.5m、2.0m、2.5m這三種設計工況的單根埋管土壤蓄熱系統(tǒng)進行模擬計算。系統(tǒng)蓄熱結束時，分析上述圖2中部所示截面的土壤溫度場，繪制不同埋深下褐土土壤的溫度在水平方向（）的分布情況，如圖5所示。

圖5 不同埋深下土壤在管長中段處x方向上的溫度變化（y=20m）

由圖中地埋管周圍土壤的溫度變化情況可以看出：不同埋深下土壤的溫度存在差距，分別對埋深為1.5m、2.0m、2.5m工況時管長中段處截面的土壤溫度場用下式進行積分。

式中：為蓄熱結束時土壤的蓄熱量；T為截面處某點蓄熱結束后的土壤溫度；為截面處對應點的初始土壤溫度；C分別為土壤蓄熱體的密度、定壓比熱。

求得埋管深度為1.5m、2.0m、2.5m時土壤蓄熱量分別為14.28GJ、19.07GJ、22.59GJ?？梢园l(fā)現(xiàn)，隨著埋深增加，土壤蓄熱量逐漸增大，蓄熱效果越來越好。但隨著埋管深度的增加，施工和材料成本也會大幅增加。

4.2 不同管內流體溫度對土壤蓄熱效果的影響

圖6 不同流體溫度對應的截面處土壤溫度分布云圖

分別對上述三個截面的土壤溫度場用公式（11）進行積分求解，得到流體進口溫度為333K、343K和353K時單根埋管的土壤蓄熱量分別為13.77GJ、16.08GJ、19.07GJ。從以上數(shù)據(jù)可以發(fā)現(xiàn)，不同管內流體溫度對太陽能跨季節(jié)土壤蓄熱效果存在影響，隨著流體溫度的升高，換熱溫差增大，流體的熱流密度增大。到蓄熱結束時，土壤的蓄熱量也逐漸增大，土壤蓄熱效果越來越好。

4.3 不同土壤的蓄熱效果

通過改變土壤的種類，對地埋管埋深=2m，流體進口溫度為353K的土壤蓄熱體的蓄熱效果進行研究，當系統(tǒng)運行整個非供暖季結束時，得到管長中段處截面土壤的溫度分布云圖，如圖7所示。

圖7 不同土壤種類對應的截面處土壤溫度分布云圖（D=2m）

分別對上述三個管長中段處截面的土壤溫度場用公式（11）進行積分求解。得到土壤類型分別為沙土、粘土和花崗巖時單根埋管的土壤蓄熱體的蓄熱量分別為4.41GJ、6.97GJ、8.92GJ。雖然從圖7看不出三種蓄熱體溫度上有多大差別，但是由于花崗巖的密度大大高于其他土壤（見表1），導致在相同蓄熱期之后，其蓄熱量也明顯高于其他土壤。即使如此，與褐色土壤相比，其蓄熱能力還是大為遜色。連同上述計算的褐土土壤，四種土壤蓄熱量的對比圖如圖8所示。

圖8 不同土壤種類的土壤蓄熱量

由上述數(shù)據(jù)可知土壤類型對埋管換熱性能及土壤溫度變化有顯著影響，導溫系數(shù)高有利于土壤中熱量的擴散，比熱容高有利于降低土壤溫度上升速率和幅度，密度高有利于土壤蓄熱，褐土是上述四種土壤中最有利于蓄熱的土壤，而沙土是最不利的。

4.4 雙埋管間距對土壤蓄熱效果的影響

為了研究雙埋管的太陽能蓄熱效果，本文建立了雙根水平地埋管土壤蓄熱體的數(shù)值模型，分別設置埋管間距為2.8m、3.2m、3.6m、4.0m四種情況，至蓄熱結束時得到如下土壤溫度分布云圖。

圖9 不同埋管間距對應的土壤溫度分布云圖

同樣對上述四個土壤溫度場用公式（11）進行積分。得到埋管間距分別為2.8m、3.2m、3.6m、4.0m時雙根埋管的土壤蓄熱體的蓄熱量分別為48.16GJ、50.06GJ、51.66GJ、52.20GJ，如圖10所示。

圖10 不同埋管間距的土壤蓄熱量

以上數(shù)據(jù)可以充分證明埋管間距對土壤的蓄熱量存在影響，隨著埋管間距的減小，土壤地埋管換熱器的附加溫變熱阻增大，溫度場之間的相互干擾加強，從而導致蓄熱量減少。

5 結論

通過對地埋管的數(shù)值模擬分析，得到如下幾點結論：

（1）埋管深度對太陽能跨季節(jié)土壤蓄熱效果存在影響。埋管深度增加，使得蓄熱體上表面受環(huán)境熱干擾減小，蓄熱結束時土壤蓄熱體的平均溫度升高，土壤的蓄熱量增大。但隨著深度的增加，其蓄熱量的增量越來越小，此外，深度增大會大幅增加系統(tǒng)的初投資，因此，埋管的深度不易取過大。

（2）不同管內流體溫度對太陽能跨季節(jié)土壤蓄熱效果存在影響。隨著流體溫度的升高，換熱溫差增大，流體的熱流密度增大，到蓄熱結束時，土壤的蓄熱量也增大。

（3）土壤類型對埋管換熱性能及土壤溫度變化有顯著影響。導溫系數(shù)高有利于土壤中熱量的擴散，比熱容高有利于降低土壤溫度上升速率和幅度，密度高有利于土壤蓄熱，褐土是上述四種土壤中最有利于蓄熱的土壤，沙土是最不利的，花崗巖和粘土介于兩者之間。

（4）埋管間距對太陽能跨季節(jié)土壤蓄熱效果存在影響。隨著埋管間距的逐漸減小，溫度場之間的相互干擾加強，導致蓄熱量的減少。適當增加埋管間距可以減小附加溫變熱阻，使得傳熱性能提高，從而使得土壤的蓄熱量增加，同時較大埋管間距蓄熱的散失量要少，土壤對蓄熱熱量的保持性好。

[1] Pinel P, Cruickshank C A, Beausoleil-Morrison I, et al. A review of available methods for seasonal storage of solar thermal energy in residential applications[J]. Renewable & Sustainable Energy Reviews, 2011,15(7):3341-3359.

[2] 沈義.我國太陽能的空間分布及地區(qū)開發(fā)利用綜合潛力評價[D].蘭州:蘭州大學,2014.

[3] 孫純武,張素云.水平埋管換熱器地熱源熱泵實驗研究及傳熱模型[J].土木建筑與環(huán)境工程,2001,23(6):49-55.

[4] 卿菁,王勇.不同埋深水平蛇形地埋管換熱器換熱性能比較分析[J].制冷與空調,2015,(5):496-501.

[5] W S Johnson. Seasonal Performance Evaluation of Two Horizontal-coil Ground-coupled Heat Pump Systems[J]. ASHRAE Trans, NT-87-19-4.

[6] 王麗.跨季節(jié)太陽能水箱蓄熱供暖系統(tǒng)的計算與研究[D].北京:華北電力大學,2012.

[7] 張文武.城市不透水表面對流換熱系數(shù)的實測和模擬研究[D].哈爾濱:哈爾濱工業(yè)大學,2008.

[8] 杜吳鵬,房小怡.北京近年地表風速和大氣混合層厚度變化特征研究[J].環(huán)境科學與技術,2017,(6):149-156.

[9] 熊莉芳,林源,李世武,等.-湍流模型及其在FLUENT軟件中的應用[J].工業(yè)加熱,2007,36(4):13-15.

[10] 曾召田,趙艷林,呂海波,等.氣象波動對水平埋管換熱器傳熱影響的數(shù)值模擬[J].太陽能學報,2018,(5):1179- 1186.

[11] 陳博,李建平.中國凍土的時空變化特征[C].中國氣象學會年會“氣候變化及其機理和模擬”分會場,2006.

Numerical Simulation of Solar Energy Storage in Trans-seasonal Soils

Guo Zhanquan Yuan Zhongxian Zhang Qi

( College of Environmental and Energy Engineering, Beijing University of Technology, Beijing, 100124 )

In order to reduce the seasonal limitation of solar energy utilization, this paper takes solar energy as a background of seasonal heat storage, and the soil is a heat storage medium. The ANSYS software is used to simulate the unsteady heat exchange process of horizontal buried pipe and soil. Under the conditions of the same flow rate, pipe diameter, initial soil temperature and running time, the effects of buried pipe depth, pipe fluid temperature, soil thermal property parameters and buried pipe spacing on soil heat storage were discussed. The results show that for the three working conditions of buried depth of 1.5m, 2.0m and 2.5m, the soil with a buried depth of 2.5m has the largest heat storage capacity and the minimum buried heat of 1.5m. For the three working conditions of fluid temperature of 60℃, 70℃, and 80℃, the temperature of the soil at 80℃ is the largest, and the heat storage at 60℃ is the smallest. For the four soils of sand, clay, cinnamon and granite, cinnamon is the most favorable for heat storage, and sand is the most unfavorable. For the four working conditions of buried pipe spacing of 2.8m, 3.2m, 3.6m and 4m, the soil heat storage increases with the increase of the buried pipe spacing. When the pipe spacing is 4m, the soil heat storage is the largest while 2.8m is the smallest.

Solar energy; horizontal buried pipe; numerical simulation; soil heat storage

TK02

A

1671-6612（2020）04-399-07

國家“973”項目（NO.2015CB251303）

郭占全（1993.10-），男，碩士，E-mail：835884151@qq.com.

苑中顯（1962.11-），男，博士，教授，E-mail：zxyuan@bjut.edu.cn

2019-11-06