地源熱泵地埋管換熱器換熱性能試驗(yàn)研究

沈德安

(中鐵第一勘察設(shè)計(jì)院集團(tuán)有限公司，陜西 西安 710043)

地源熱泵地埋管換熱器換熱性能試驗(yàn)研究

沈德安

(中鐵第一勘察設(shè)計(jì)院集團(tuán)有限公司，陜西 西安 710043)

影響地源熱泵系統(tǒng)效果的關(guān)鍵因素是室外地埋管換熱器的設(shè)置情況，影響地埋管換熱器換熱性能的主要因素是巖土的綜合導(dǎo)熱系數(shù)、綜合比熱容以及鉆孔的總熱阻等熱物性參數(shù)。各工程所在地地質(zhì)條件差異很大，如何得到地埋管換熱器換熱性能參數(shù)是地源熱泵系統(tǒng)設(shè)計(jì)必須解決的問題。本文基于線熱源模型得出換熱器內(nèi)流體平均溫度與時(shí)間的對(duì)數(shù)關(guān)系曲線，通過測試地埋管換熱器進(jìn)出口水溫，進(jìn)而推算出巖土的綜合導(dǎo)熱系數(shù)及熱阻，大大簡化了計(jì)算過程。此外，通過對(duì)地埋管測試井長期連續(xù)監(jiān)測發(fā)現(xiàn)，在寒冷或嚴(yán)寒地區(qū)對(duì)于冬夏季冷熱負(fù)荷不平衡的系統(tǒng)，室外地埋管換熱區(qū)巖土溫度呈下降趨勢(shì)，設(shè)計(jì)時(shí)應(yīng)考慮解決措施。

地源熱泵 地埋管換熱器 換熱特性 線熱源

地源熱泵技術(shù)是利用電能等高位能源從地表淺層提取地?zé)豳Y源的一種新能源運(yùn)用技術(shù)，地埋管地源熱泵系統(tǒng)由換熱器、熱泵主機(jī)和末端用戶三部分組成，其中影響系統(tǒng)的關(guān)鍵因素是地埋管換熱器的設(shè)計(jì)與施工［1］。目前，工程上普遍采用的方法是通過冷熱響應(yīng)模擬地源熱泵系統(tǒng)實(shí)際運(yùn)行工況，在施工安裝前進(jìn)行現(xiàn)場測試［2-3］。對(duì)地源熱泵U型管地下?lián)Q熱器運(yùn)行工況的模擬，當(dāng)前國內(nèi)外工程中使用較多的是線熱源模型和柱熱源模型［4-5］。影響地埋管換熱器性能的直接因素是巖土的綜合導(dǎo)熱系數(shù)、綜合比熱容和鉆孔總熱阻等熱物性參數(shù)。

近年來，很多地源熱泵工程實(shí)際應(yīng)用中出現(xiàn)了問題，主要是由于地埋管換熱器的選型設(shè)計(jì)或施工不當(dāng)，導(dǎo)致實(shí)際運(yùn)行效果不好，亦或是工程造價(jià)過高等，這主要是因?yàn)楦鞯貛r土的熱物性差異很大［6］。據(jù)文獻(xiàn)［7］，地下巖土導(dǎo)熱系數(shù)發(fā)生10%的偏差便會(huì)造成設(shè)計(jì)的地下埋管長度產(chǎn)生4.5% ～5.8%的偏差，最終導(dǎo)致鉆孔和地埋管長度發(fā)生偏差，影響工程的運(yùn)行效果和工程造價(jià)。在工程測試中，為推算出巖土的熱物性參數(shù)，首先應(yīng)建立合理的地埋管換熱器的傳熱計(jì)算模型?，F(xiàn)在有關(guān)換熱器的計(jì)算模型有很多文獻(xiàn)，目前已經(jīng)提出并公開發(fā)表的已經(jīng)有30余種模型［8-10］。研究表明［11］，在換熱初期階段，柱熱源和線熱源模型的計(jì)算結(jié)果會(huì)有顯著差別，但隨著換熱時(shí)間的延長，兩種模型計(jì)算結(jié)果的相對(duì)誤差會(huì)逐步縮小，而且時(shí)間越長誤差越小，因此在工程研究中完全可以采用線熱源模型進(jìn)行分析計(jì)算，而應(yīng)用圓柱面模型會(huì)使計(jì)算變得復(fù)雜許多。本文結(jié)合具體的太陽能耦合地源熱泵系統(tǒng)工程來研究地源熱泵地埋管換熱器換熱性能。

1 工程概況

工程位于山西省五寨縣三岔站站區(qū)，項(xiàng)目于2008年8月開展設(shè)計(jì)，2010年7月投入運(yùn)行，至今系統(tǒng)運(yùn)行狀況良好。整個(gè)太陽能耦合地源熱泵系統(tǒng)為站區(qū)綜合樓提供空調(diào)、采暖及生活熱水供應(yīng)，綜合樓冬季采暖設(shè)計(jì)熱負(fù)荷為1 420 kW，夏季空調(diào)設(shè)計(jì)冷負(fù)荷為1 016 kW，生活熱水最大耗熱量157 kW。機(jī)房內(nèi)設(shè)置2臺(tái)螺桿式地源熱泵機(jī)組，壓縮機(jī)為半封閉式雙螺桿壓縮機(jī)，采用4段容量調(diào)節(jié)方式，即100% ～75% ～50%～25%，其中25%為起動(dòng)過程。通過人性化的液晶PLC中文界面操作微電腦控制器，能隨時(shí)監(jiān)控機(jī)組運(yùn)行狀態(tài);空調(diào)系統(tǒng)末端設(shè)備為風(fēng)機(jī)盤管。室外共設(shè)地埋管換熱井536口，井深100 m，其中有494口雙 U型地埋管換熱器井和42口單U型地埋管換熱器井。通過水平集水管和分支管，將4個(gè)并聯(lián)環(huán)路連接成一個(gè)完整的換熱系統(tǒng)。當(dāng)?shù)夭膳谧援?dāng)年10月底至次年的4月初，約5個(gè)月?？照{(diào)期自6月至8月，約3個(gè)月。

2 地埋管換熱器傳熱線熱源模型

建立線熱源模型，首先應(yīng)作以下假設(shè):

1)將鉆孔內(nèi)當(dāng)作一個(gè)整體，將地埋管換熱器當(dāng)作一個(gè)當(dāng)量直徑的線熱源;

2)地埋管換熱器按照輻射狀向周圍巖土傳熱;

3)將管子周圍的巖土和回填材料看成是無限大的實(shí)體。

基于以上假設(shè)，換熱井周圍的傳熱實(shí)際上可以簡化為一維軸對(duì)稱問題。由于鉆孔半徑較小，該模型的數(shù)學(xué)表達(dá)式可以簡化為［12］

式中:T=T(r，t)，是t時(shí)刻r處的巖土溫度，℃;λZ為巖土綜合導(dǎo)熱系數(shù)，W/(m·K);ρs為巖土密度，kg/m3;Cs為巖土的比熱容，kJ/(kg·k);T0為未受擾動(dòng)時(shí)的巖土原始溫度，℃;rb為鉆孔半徑，m;ql為每延米長線熱源的熱流強(qiáng)度，W/m。

根據(jù)上述線熱源模型，可以得出其理論解析解，流體的平均溫度Tf(t)為

根據(jù)地埋管換熱器的線熱源模型理論，換熱量恒定時(shí)，系統(tǒng)循環(huán)流體的平均溫度與時(shí)間存在對(duì)數(shù)關(guān)系，即

巖土綜合導(dǎo)熱系數(shù)λZ和鉆孔內(nèi)的總熱阻R0的計(jì)算公式如下

式中，k為擬合曲線的斜率值。

3 試驗(yàn)測試

試驗(yàn)測試過程中，在地埋管換熱區(qū)中間和邊緣共選取2口有代表性的換熱井進(jìn)行溫度的測量，溫度傳感器分別安裝在距地面5，35，65和95 m的位置。溫度測點(diǎn)布置如圖1所示。

圖1 換熱井內(nèi)溫度測點(diǎn)分布(單位:m)

測試結(jié)果表明，原始無熱源干擾時(shí)循環(huán)溫度在10.5℃左右波動(dòng)，巖土初始溫度可視為10.5℃。試驗(yàn)測試參數(shù)如表1所示。

表1 試驗(yàn)測試參數(shù)

4 地埋管換熱器性能分析

4.1 巖土熱物性分析

圖2是2011年8月地埋管換熱區(qū)巖土溫度變化曲線?？梢钥闯觯瑴\層巖土溫度波動(dòng)十分明顯，地下5 m處的巖土溫度波動(dòng)較大，隨著巖土深度增加，溫度波動(dòng)越來越小。這是因?yàn)闇\層巖土受地面和氣候環(huán)境的影響比較大，而深層的巖土受到的影響較小。此外巖土溫度的分布呈現(xiàn)一定的梯度，隨著深度的增加有一定的溫升，但從地下5 m到95 m之間，巖土溫度升高在1℃以內(nèi)，這說明地埋管區(qū)溫度的分布還是比較均勻的，設(shè)計(jì)中可以把它看成一個(gè)較為均勻的溫度場。

圖2 8月巖土溫度變化曲線

圖3 測試井中7月底至次年2月初巖土溫度變化曲線

圖3是從2011年7月底至2012年2月初的兩個(gè)測試井內(nèi)巖土溫度變化曲線。由圖可以看出，巖土溫度剛開始緩慢上升，這是因?yàn)楫?dāng)?shù)叵募究照{(diào)期比較短，而且夏季空調(diào)冷負(fù)荷也比較小。從10月底開始供暖后，巖土溫度出現(xiàn)持續(xù)下降，這是因?yàn)楸卷?xiàng)目處于寒冷地區(qū)，冬季供暖期較長，熱負(fù)荷較大。2012年2月初的巖土溫度與巖土初始溫度(10.5℃)相比下降了將近3℃，因此本項(xiàng)目所研究的地源熱泵系統(tǒng)在冬夏季運(yùn)行過程中，其吸排熱量未能夠達(dá)到平衡。

圖4所示為地源熱泵系統(tǒng)連續(xù)穩(wěn)定運(yùn)行情況下，根據(jù)實(shí)測數(shù)據(jù)所擬合的地埋管進(jìn)出水平均溫度與運(yùn)行時(shí)間關(guān)系曲線。

圖4 地埋管進(jìn)出水平均溫度與運(yùn)行時(shí)間關(guān)系曲線

由圖4可知，擬合曲線斜率 k等于1.538 3，根據(jù)實(shí)測數(shù)據(jù)，ql為15.3 W/m，故由式(4)可計(jì)算出 λZ等于0.792 W/(m·℃)。本工程當(dāng)?shù)赝临|(zhì)以砂土為主，查有關(guān)資料得巖土密度ρs=1 374 kg/m3，巖土的比熱容Cs=1.51 kJ/(kg·k)［13］，由式(5)計(jì)算可得 R0=0.314 m·k/W。

4.2 地埋管換熱器換熱特性分析

1)夏季工況

圖5所示為2011年7月26日、8月9日和8月17日這三天中某一時(shí)段地埋管換熱器的進(jìn)出水溫度及換熱量變化曲線。由圖5可以看到，在不同的進(jìn)水溫度、相近的換熱時(shí)間內(nèi)，換熱器的出水溫度比較接近。在8月17日，進(jìn)水溫度有上升趨勢(shì)，但經(jīng)過一段時(shí)間后，出水溫度最終接近7月26日的溫度值。這是因?yàn)樵谠撓到y(tǒng)中，夏季地埋管換熱器的釋熱能力充分，巖土的換熱能力充足，完全滿足夏季制冷需求。夏季地埋管換熱器平均換熱性能測試結(jié)果如表2所示。

圖5 夏季地埋管換熱器進(jìn)出水溫度及換熱量變化曲線

表2 夏季地埋管換熱器的換熱性能測試結(jié)果

2)冬季工況

圖6所示為系統(tǒng)在冬季供暖過程中，地埋管換熱器進(jìn)出水溫度的變化曲線?？梢钥闯?，在系統(tǒng)剛啟動(dòng)階段，地埋管換熱器進(jìn)出水溫度比較穩(wěn)定，隨著換熱持續(xù)，換熱器進(jìn)出水溫度先呈現(xiàn)緩慢下降的趨勢(shì)，但隨著換熱過程的延長，換熱器又逐漸進(jìn)入平穩(wěn)換熱狀態(tài)。

圖6 地埋管換熱器進(jìn)出水溫度曲線

圖7和圖8分別為2011年1月14日、1月25日和2月10日某一時(shí)段地埋管換熱器進(jìn)出口水溫、每延米井深的換熱量變化曲線。

圖7 冬季地埋管換熱器進(jìn)出水溫度變化曲線

圖8 冬季換熱器每延米井深換熱量變化曲線

從圖7中可以看出，每天的進(jìn)出水溫度基本穩(wěn)定，但系統(tǒng)持續(xù)運(yùn)行一段時(shí)間后，地埋管換熱器出水溫度有所下降。從圖8可以看出，1月14日每延米井深換熱量約為11 W/m，到1月25日后略有下降，到2月10日每延米井深換熱量僅約為10 W/m。這說明該系統(tǒng)在冬季連續(xù)運(yùn)行后，地埋管換熱器的換熱能力有所減弱，但并未出現(xiàn)明顯的惡化。這主要是因?yàn)殚L期的取熱過程中，地埋管換熱區(qū)的巖土溫度沒有及時(shí)恢復(fù)，故地埋管換熱器的換熱能力有所下降。根據(jù)現(xiàn)場測試來看，該供暖季巖土溫度下降對(duì)系統(tǒng)沒有造成明顯的影響，系統(tǒng)可以穩(wěn)定運(yùn)行。須說明的是，地源熱泵系統(tǒng)是否能常年維持穩(wěn)定高效運(yùn)行，取決于非采暖季節(jié)是否有足夠的熱量補(bǔ)償至室外地埋管區(qū)域。根據(jù)前面關(guān)于巖土溫度變化的分析，可以看出系統(tǒng)常年運(yùn)行后會(huì)出現(xiàn)換熱量減小、效率降低的趨勢(shì)。本工程設(shè)計(jì)中預(yù)留有太陽能跨季節(jié)蓄熱措施，可以解決此問題。冬季地埋管換熱器平均換熱性能測試結(jié)果如表3所示。

表3 冬季地埋管換熱器換熱性能測試結(jié)果

5 結(jié)論

1)在實(shí)際工程測試中，基于線熱源模型研究地源熱泵室外地埋管換熱器換熱性能及巖土熱物性的方法切實(shí)可行，根據(jù)線熱源模型理論解析解，可以計(jì)算出巖土綜合導(dǎo)熱系數(shù)λZ和鉆孔內(nèi)的總熱阻R0等熱物性參數(shù)。

2)對(duì)于寒冷或嚴(yán)寒地區(qū)，冬夏季冷熱負(fù)荷嚴(yán)重不平衡的地源熱泵系統(tǒng)，根據(jù)一個(gè)空調(diào)季和一個(gè)采暖季對(duì)地埋管區(qū)巖土溫度監(jiān)測結(jié)果顯示，埋管區(qū)巖土溫度會(huì)呈現(xiàn)持續(xù)下降趨勢(shì)。這會(huì)影響地埋管換熱器換熱量大小以及地源熱泵系統(tǒng)性能，需要采取措施防止系統(tǒng)運(yùn)行性能惡化。

3)淺層巖土溫度波動(dòng)較大，深層巖土溫度波動(dòng)較小，且隨著巖土深度的增加巖土溫度有所升高，但從地下5 m到95 m，巖土溫度升高在1℃以內(nèi)。在地源熱泵工程設(shè)計(jì)中，可以將地埋管換熱區(qū)視為一個(gè)均勻溫度場進(jìn)行設(shè)計(jì)分析。

［1］程宗科，柯軍，劉秋新，等．地源熱泵地埋管換熱設(shè)計(jì)計(jì)算［J］．制冷與空調(diào)，2008，22(2):58-62．

［2］RYBACH L，SANNER B．Ground-Source Heat Pump Systems:The European Experience［J］．GHC Bulletin，2000(1):16-26．

［3］RAFFERTY K．Design Issues in the Commercial Application of GSHP Systems in the U．S．［J］．GHC Bulletin，2000(1):6-10．

［4］GU Yian，O’NEAL D L．An Analytical Solution to Transient Heat Conduction in a Composite Region with a Cylindrical Heat Source［J］．Journal of Solar Energy Engineering，1995，117(3):242-248．

［5］柳曉雷，王德林，方肇洪．垂直埋管地源熱泵的圓柱面?zhèn)鳠崮Ｐ图昂喕?jì)算［J］．山東建筑工程學(xué)院學(xué)報(bào)，2001，16(1):47-51．

［6］李敏．地埋管地源熱泵系統(tǒng)的區(qū)域適應(yīng)性研究［D］.南京:南京理工大學(xué)，2010．

［7］KAVANSUGH S P．Field Tests for Gound Thermal Properties Methods and Impact on Ground-Source Heat Pumps［J］．ASHRAE Transactions，1998，104(2):347-355．

［8］仇君，李莉葉，朱曉慧．地源熱泵埋管換熱器傳熱模型的研究進(jìn)展［J］．廣西大學(xué)學(xué)報(bào):自然科學(xué)版，2007，32(增):1-4，7．

［9］潘昌榮，白辛民，王永波，等．利用地?zé)崽岣呤彝獾叵陆o水井內(nèi)部溫度的新方法［J］．鐵道建筑，2012(9):132-134．

［10］荊志東，王春雷，謝強(qiáng)．凍區(qū)鹽漬土導(dǎo)熱系數(shù)測試方法的研究［J］．鐵道建筑，2008(3):69-70．

［11］張長興，胡松濤，李緒泉．基于線熱源模型的地源熱泵系統(tǒng)運(yùn)行模擬［J］．建筑科學(xué)，2008，24(12):71-74．

［12］于明志，彭曉峰，方肇洪，等．基于線熱源模型的地下巖土熱物性測試方法［J］．太陽能學(xué)報(bào)，2006，27(3):279-281．

［13］張熙民，任澤霈，梅飛鳴．傳熱學(xué)［M］．北京:中國建筑工業(yè)出版社，2001．

Experimental study on heat transfer performance of ground source heat pump buried tube heat exchanger

SHEN De'an

(China Railway First Survey and Design Institute Group Co．，Ltd．，Xi'an Shaanxi 710043，China)

The key of the ground source heat pump system is the design of the ground heat exchanger．Some soil thermal properties including soil heat conductivity，specific heat capacity and entire thermal resistance of the well drilling affect heat transfer characteristics of ground heat exchanger significantly．However，geological conditions are different in different area，resulting in the difficulty in acquiring ground heat exchanger heat transfer capacity in every project designing．In this paper，the log curve of the relationship between average temperature and time was obtained based on the line heat resource theory．According to the test of the inlet and outlet water temperature，the soil heat conductivity and thermal resistance were calculated，leading to a simpler calculation．Additionally，it was revealed that when summer cooling load and winter heating load was unbalanced in cold area，the temperature dropped in buried pipes outdoors，requiring solutions at the design stage．

Ground source heat pump;Ground heat exchanger;Heat transfer characteristics;Linear heat resource

TK523

A

10．3969/j．issn．1003-1995．2014．10．28

1003-1995(2014)10-0111-05

2014-06-20;

2014-07-25

沈德安(1981— )，男，云南騰沖人，工程師，碩士。

(責(zé)任審編 李付軍)