張洋，程海峰，張舉，楊葉萌，王庚

（1.安徽建筑大學(xué)環(huán)境與能源工程學(xué)院，安徽 合肥 230601；2.建筑能效控制與評估教育部工程研究中心，安徽 合肥 230601；3.智能建筑與建筑節(jié)能安徽省重點(diǎn)實(shí)驗(yàn)室，安徽 合肥 230601；4.安徽建筑大學(xué)建筑室內(nèi)熱濕環(huán)境實(shí)驗(yàn)室，安徽 合肥 230601）

1 概述

本文通過對宿州地區(qū)單U深井式地埋管換熱器建立地埋管換熱器的管內(nèi)外耦合全尺寸三維動態(tài)數(shù)值計(jì)算模型，對該地埋管換熱器的換熱性能等進(jìn)行數(shù)值仿真模擬并與實(shí)測結(jié)果進(jìn)行比較分析。項(xiàng)目地點(diǎn)位于宿州市埇橋區(qū)（其地理坐標(biāo)為：東經(jīng) 11656'28"～11701'09"，北緯 3337'52"～3340'42"）。鉆井深度100m，回填材料為原土回填，其中單U地埋管采用高密度聚乙烯材質(zhì)（HDPE），管徑為20mm，地下鉆井直徑140mm。單U地埋管換熱器截面和剖面結(jié)構(gòu)如圖1、圖 2。

圖1 單U換熱器橫截面圖

圖2 單U換熱器剖面圖

2 模型的建立

2.1 基本假設(shè)

適當(dāng)?shù)募僭O(shè)能夠在確保研究數(shù)據(jù)準(zhǔn)確度為前提的條件下加快模型運(yùn)算速度，合理的假設(shè)對于解決問題而言十分重要。為了更好地分析，按照實(shí)際情況做出下列的假設(shè)：

①假設(shè)地表溫度變化忽略不計(jì)，土壤的初始溫度在無窮遠(yuǎn)處也不發(fā)生變化且均勻一致；

②假設(shè)換熱器周圍的復(fù)雜固液混合物為均勻固體，忽略周圍巖土體滲流帶來的影響，熱量交換僅以導(dǎo)熱形式進(jìn)行；

③假設(shè)換熱管內(nèi)流體為理想工況即管內(nèi)同一截面的流體具有相同的溫度與流速；

④假設(shè)換熱裝置周圍為均勻的巖石、土層和回填料，往徑向和豎向延伸，其熱物性參數(shù)不變；

⑤假設(shè)忽略回填材料的接觸熱阻，地埋管與周圍巖土以線熱源方式進(jìn)行換熱。

2.2 三維模型的建立

根據(jù)以上假設(shè)，利用COMSOL軟件建立模型，如圖3、圖4。

圖3 頂部模型圖

圖4 底部模型圖

2.3 三維模型網(wǎng)格分析

本模型外部土壤模擬為半徑5m內(nèi)土壤換熱，井內(nèi)PE管直徑僅為20mm，構(gòu)建網(wǎng)格較為復(fù)雜，為使模擬計(jì)算盡可能精準(zhǔn)快捷，所以采用分區(qū)構(gòu)建的方法。構(gòu)建如圖5、圖6。

圖5 換熱系統(tǒng)網(wǎng)格

圖6 底部單U部分網(wǎng)格

3 模型的驗(yàn)證

為了驗(yàn)證模型，根據(jù)地源熱泵運(yùn)轉(zhuǎn)情況實(shí)測分析，得到包括進(jìn)出口水溫在內(nèi)數(shù)據(jù)，查閱資料及對宿州地區(qū)進(jìn)行巖土熱響應(yīng)實(shí)驗(yàn)，可以得到COMSOL模型各部分熱物性。宿州地區(qū)地屬華北平原，地層大多屬于第四系，土壤類型以粉質(zhì)黏土、粉土和粉質(zhì)沙土為主，土壤平均熱導(dǎo)率為 1.8W/（m·K），平均比熱容為2.8MJ/（m·K），初始地溫為18℃。地層平均初始溫度隨時間變化曲線圖如圖7。

圖7 地層平均初始溫度隨時間變化

巖土體平均初始溫度測試歷時5.8h，在未向地埋管提供冷、熱量的情況下保持地埋管內(nèi)水以1.5m/h的流量循環(huán)流動，測得的循環(huán)水溫度即為巖土的平均初始溫度。經(jīng)測定，試驗(yàn)孔巖土平均初始溫度為18.2℃，即初始地溫。

根據(jù)實(shí)際的進(jìn)水溫度，經(jīng)過模擬，可得實(shí)際與模擬所得的出水溫度曲線如圖8。

圖8 出口溫度模擬計(jì)算結(jié)果與實(shí)測結(jié)果對比圖

本次模擬為從當(dāng)日下午8點(diǎn)到次日下午7點(diǎn)，每隔1小時記錄一次出口溫度，共模擬了機(jī)組在夏季穩(wěn)定運(yùn)行24小時的工作狀況，與相同工況下，機(jī)組實(shí)際穩(wěn)定運(yùn)行24h所測得的出口溫度進(jìn)行對比。由圖8可知，實(shí)測出口溫度比模擬值較低，是因?yàn)橥寥乐杏械叵滤拇嬖冢酝寥乐写嬖跓釋α?，而本模擬把土壤的傳熱看作為純導(dǎo)熱處理，因此實(shí)際的傳熱效果要好于模擬效果。通過查閱宿州地區(qū)水文地質(zhì)資料，宿州地區(qū)100m深度土壤溫度大概在18℃，與實(shí)驗(yàn)?zāi)M設(shè)定地溫一致。而試驗(yàn)出口水溫值取穩(wěn)定運(yùn)行狀態(tài)下24h時間段的數(shù)據(jù)，可以看出實(shí)驗(yàn)與模擬所對應(yīng)的趨勢一致，出口水溫值較為吻合，誤差不足實(shí)際數(shù)值的5%，二者之間最大差值僅為0.6℃，即驗(yàn)證了本文所采用的的數(shù)值模擬和計(jì)算方法的合理性。

4 換熱效率影響因素的分析

4.1 循環(huán)流體流速

實(shí)際測算過程中發(fā)現(xiàn)，循環(huán)流體的流速對系統(tǒng)換熱性能有較大影響。為了探究循環(huán)流體流速對地?zé)峋畵Q熱性能的具體影響，本文展開具體的分析研究。本文第2、3章所建傳熱模型模擬流速與實(shí)際機(jī)組運(yùn)行循環(huán)流體流速一致，為1.5m/h，且經(jīng)過驗(yàn)證符合實(shí)際情況。所以在原模型基礎(chǔ)上改變流速，得到不同流速下相同時刻的進(jìn)、出口溫度模擬如圖9～圖13。

圖9 流速為0.5m3/h

圖10 流速為2.5m3/h

圖11 流速為3.5m3/h

圖12 流速為4.5m3/h

圖13 流速1.5m3/h，進(jìn)口水溫26℃

由圖9～圖13可知，當(dāng)循環(huán)流體流速分別為0.5 m/h，1..5 m/h，2.5 m/h，3.5 m/h，4.5 m/h逐漸增加時，系統(tǒng)進(jìn)水溫度一定，設(shè)定為26℃，出水口溫度逐漸升高。假設(shè)以穩(wěn)定進(jìn)口流速工作，不同流速各自工作一小時，出口溫度及換熱量模擬數(shù)據(jù)繪制如圖14、圖15。

圖14 出口溫度隨循環(huán)體流速變化圖

圖15 換熱量隨循環(huán)體流速變化圖

由圖14、圖15可知，當(dāng)循環(huán)流體流速逐漸增加時，出口溫度逐漸增大，溫差減小，但因?yàn)榱魉俅?，循環(huán)水體流量大，所以整體上系統(tǒng)的換熱量仍處于增加趨勢。而且可以看出，當(dāng)流速較小時，流速的增加，系統(tǒng)換熱能力顯著增加，當(dāng)流量增加到一定程度，換熱量增加的趨勢逐漸變緩并趨于穩(wěn)定。而對于地源熱泵系統(tǒng)，流量越大，循環(huán)水泵功耗也就越大，系統(tǒng)功耗也越大。所以循環(huán)體流速對地埋管換熱器換熱效率影響很大，但并非流速越大，系統(tǒng)就越節(jié)能，要根據(jù)實(shí)際需求及系統(tǒng)需要的出水溫度，結(jié)合系統(tǒng)功耗及換熱量，比較換熱量和水泵功耗，確定最合理的循環(huán)體流速，以保證系統(tǒng)在最節(jié)能的方式下正常穩(wěn)定運(yùn)行。由分析可知，進(jìn)水流速達(dá)到3.5m/h，系統(tǒng)換熱量基本穩(wěn)定，最適宜的入口流速區(qū)推薦在1.3m/h～2.3m/h。

4.2 U型管進(jìn)口水溫

因?yàn)閱蜺型地埋管換熱器進(jìn)口水溫不同會直接影響系統(tǒng)在地下的換熱量，所以在確定系統(tǒng)的進(jìn)水口流速情況下，分別設(shè)置不同的進(jìn)水溫度，根據(jù)已驗(yàn)證的COMSOL模型，得到各個不同進(jìn)口水溫條件下，進(jìn)出口溫度分布圖17～圖20。

由圖13，圖16～圖19可知，當(dāng)循環(huán)水流速一定時，換熱系統(tǒng)分別以22℃、26℃、30℃、34℃、38℃的進(jìn)口水溫工作，出口水溫也有較大差別，為逐漸升高趨勢。假設(shè)系統(tǒng)以1.5m/h的進(jìn)口流速下穩(wěn)定運(yùn)行一小時，由模型計(jì)算，可得出口溫度及換熱量在不同進(jìn)口水溫下的變化如圖20、圖21。

圖16 進(jìn)口水溫22℃

圖17 進(jìn)口水溫30℃

圖18 進(jìn)口水溫34℃

圖19 進(jìn)口水溫38℃

圖20 出口溫度隨循環(huán)水進(jìn)口溫度變化圖

圖21 換熱量隨循環(huán)水進(jìn)口溫度變化圖

由圖20、圖21可知，當(dāng)循環(huán)水進(jìn)口水溫逐漸增大時，出口水溫隨之增大，換熱量也逐漸增大，但增大趨勢明顯減小，由此可知進(jìn)口溫度的高低對地下?lián)Q熱系統(tǒng)的性能具有較大影響。而且系統(tǒng)的換熱能力是受到實(shí)際條件限制的，夏季制冷工況，過高的進(jìn)水溫度，會使出水溫度更高，系統(tǒng)的制冷效率降低，熱泵機(jī)組效率下降。當(dāng)?shù)芈窆艹隹谒疁爻^熱泵機(jī)組進(jìn)口水溫時，則會導(dǎo)致機(jī)組停機(jī)保護(hù)，甚至可能會使地下土壤層造成“熱堆積”現(xiàn)象，對環(huán)境和周圍生態(tài)平衡造成損害。因此，換熱器內(nèi)循環(huán)水進(jìn)口溫度的選擇應(yīng)該綜合考慮系統(tǒng)換熱量及熱泵機(jī)組的性能等因素，在滿足系統(tǒng)正常穩(wěn)定運(yùn)行的前提下確定最佳的循環(huán)流體進(jìn)口溫度。由分析可知，進(jìn)水溫度達(dá)到37℃，換熱量基本不變，最適宜的進(jìn)水溫度區(qū)推薦在26℃～30℃。

5 結(jié)論

①進(jìn)口流速和進(jìn)水溫度對地埋管換熱器換熱量都有很大影響。隨著流速增大，進(jìn)水溫度升高，換熱量總體上呈現(xiàn)上升趨勢，最后分別在進(jìn)口流速3.5 m/h，進(jìn)水溫度37℃時達(dá)到穩(wěn)定。

②當(dāng)循環(huán)水進(jìn)口流速不同時，流速越大，系統(tǒng)的出水溫度越大，換熱量越大，但加大流速會增加水泵的功耗，所以進(jìn)口流速在 1.3m/h～2.3m/h范圍內(nèi)時，換熱量增幅最大，換熱效率最高。

③當(dāng)進(jìn)水溫度不同時，進(jìn)水溫度越高，系統(tǒng)的出水溫度也越高，換熱量隨之增大，但過高溫度會使熱泵制冷效率下降，所以進(jìn)水溫度應(yīng)控制在26℃～30℃，此時系統(tǒng)換熱效率最高。