太陽能熱泵地板供暖系統(tǒng)測(cè)試與經(jīng)濟(jì)性研究

王　琛，趙　麒

(1. 長春工程學(xué)院 設(shè)計(jì)研究院 長春 130012；2. 長春工程學(xué)院能源動(dòng)力工程學(xué)院 長春 130012)

?

太陽能熱泵地板供暖系統(tǒng)測(cè)試與經(jīng)濟(jì)性研究

王琛1，趙麒2

(1. 長春工程學(xué)院 設(shè)計(jì)研究院 長春 130012；2. 長春工程學(xué)院能源動(dòng)力工程學(xué)院 長春 130012)

在中國北方城市一住宅內(nèi)設(shè)置了太陽能熱泵地板供暖裝置，測(cè)試了系統(tǒng)的冬季運(yùn)行情況，結(jié)果顯示系統(tǒng)滿足該用戶的冬季供暖要求.并基于有限時(shí)間熱力學(xué)理論和集熱器熱損失模型，建立了太陽能熱泵地板供暖系統(tǒng)的熱力學(xué)模型，并對(duì)該系統(tǒng)進(jìn)行了熱經(jīng)濟(jì)分析.研究在給定供熱率和初投資的條件下，以系統(tǒng)的供熱系數(shù)作為熱經(jīng)濟(jì)性目標(biāo)函數(shù)，得出了在目標(biāo)函數(shù)取最大時(shí)系統(tǒng)最佳的運(yùn)行性能系數(shù)和設(shè)計(jì)參數(shù).同時(shí)還研究了初投資對(duì)系統(tǒng)運(yùn)行以及設(shè)計(jì)參數(shù)的影響，得出了對(duì)應(yīng)給定供熱率系統(tǒng)的最佳初投資及其相應(yīng)的設(shè)計(jì)參數(shù).

有限時(shí)間熱力學(xué)；太陽能熱泵；熱經(jīng)濟(jì)性；目標(biāo)函數(shù)

在能源和環(huán)境問題日益嚴(yán)峻的今天，太陽能技術(shù)因其具有顯著的節(jié)能性和環(huán)境友好性，得到了越來越廣泛的關(guān)注[1-3].其發(fā)展前景廣闊，在嚴(yán)寒地區(qū)應(yīng)用該系統(tǒng)進(jìn)行供暖獲得了很好的運(yùn)行效果[4].該系統(tǒng)以太陽熱作為熱泵的低位熱源，熱泵的蒸發(fā)溫度較高，供熱系數(shù)較其它熱泵供暖系統(tǒng)高.集熱器溫度升高會(huì)使熱泵系能系數(shù)提高.但是在一定的負(fù)荷下，只增加集熱器面積會(huì)使集熱器的工作溫度升高，導(dǎo)致集熱器的熱損失增加，集熱效率降低，使得整個(gè)系統(tǒng)的運(yùn)行性能系數(shù)增加很小，系統(tǒng)的經(jīng)濟(jì)性下降.有限時(shí)間熱力學(xué)為探索熱力過程的性能界限、熱力學(xué)優(yōu)化問題提供了有效的方法.Bejan[5]研究了傳熱不可逆情況下制冷裝置冷熱換熱面積的最佳分配方法.陳林根等[6]研究了帶有變溫?zé)嵩吹目稍偕豢赡骈]式Joule-Brayton循環(huán)性能.

目前的研究主要集中在熱機(jī)和熱泵裝置本身的最佳工作溫度和面積的最佳優(yōu)化分配上，本文應(yīng)用有限時(shí)間熱力學(xué)理論，考慮熱泵內(nèi)部不可逆因素的影響，同時(shí)考慮集熱器線性熱損失情況下，建立了太陽能熱泵地板供暖系統(tǒng)熱力學(xué)模型，進(jìn)行熱經(jīng)濟(jì)性分析，為該系統(tǒng)的優(yōu)化設(shè)計(jì)提供理論依據(jù).

1　系統(tǒng)模型

如圖1所示，系統(tǒng)主要由太陽集熱器、蓄熱裝置、熱泵蒸發(fā)器、熱泵冷凝器和地板輻射供暖裝置構(gòu)成.蓄熱裝置的作用是緩解在一天內(nèi)太陽輻射的不均勻性，提高集熱器集熱量，保證系統(tǒng)持續(xù)高效運(yùn)行，增加了系統(tǒng)的供暖穩(wěn)定性.

1-太陽能集熱器；2-蓄熱水箱；3-蒸發(fā)器；4-冷凝器；5-熱用戶；6-壓縮機(jī)；7-采暖循環(huán)泵；8-節(jié)流機(jī)構(gòu)；9，10-循環(huán)泵圖1　太陽能熱泵地板供暖系統(tǒng)原理圖

2　運(yùn)行測(cè)試

集熱器水溫直接影響熱泵的性能,而集熱器水溫又依賴太陽的輻射強(qiáng)度，所以有必要對(duì)太陽輻射強(qiáng)度及低溫側(cè)水溫進(jìn)行測(cè)試研究.

2.1測(cè)試條件

試驗(yàn)地點(diǎn)選擇長春市世紀(jì)大街與武漢路交匯處某小區(qū)的一間頂層民房內(nèi).北向外墻490mm厚紅磚外墻，外窗為1 500×1 500mm雙玻塑鋼窗；西側(cè)間墻為240mm厚紅磚墻，隔壁房間未供暖；東側(cè)間墻為240mm厚紅磚墻，相鄰房間為廚房，未供暖；南側(cè)間墻為240mm厚紅磚墻，門口尺寸為900×2 000mm；頂棚為120mm鋼筋混凝土樓板外加100mm厚容重10kg/m3苯板保溫層；地面為低溫地板輻射供熱面.

低溫地板輻射供熱結(jié)構(gòu)為樓板上加設(shè)一層厚度為20mm的苯板保溫層，減少熱量向樓下房間傳送；在苯板上鋪設(shè)反光層，減少輻射向下傳送；以200mm間距在反光層上盤地?zé)峁?，共分兩環(huán)，每環(huán)長不超過65m(房間面積為3.3×5.1=16.83m2)共約120m.管與管之間用小粒徑河卵石填充，最上部為20mm厚水泥砂漿保護(hù)層.集熱裝置選擇天普全玻璃真空管太陽熱水器，玻璃真空管管外徑為58mm，每根長為1 800mm，共20根，外形尺寸為1 480mm×1 610mm(左右寬×前后寬).太陽熱水器的集熱面積為1.8m2(0.05×1.8×20).安裝位置為樓頂坡屋面，朝向正南，傾角約為45°.貯熱水箱保溫材料使用硬質(zhì)無氟聚氨酯泡沫塑料，性能見表1.

表1　硬質(zhì)無氟聚氨酯泡沫塑料性能

2.2測(cè)試結(jié)果

表2 為當(dāng)月室外氣溫最低日即2013年12月22日和2014年1月4日的太陽輻射強(qiáng)度值，可以看出太陽輻射強(qiáng)度最大值分別出現(xiàn)在上午12點(diǎn)和11點(diǎn)，最小值出現(xiàn)在早晨8：00(太陽剛升起)和下午15：00～16：00(即將日落).

表2　太陽輻射強(qiáng)度測(cè)試

圖2顯示了采暖期監(jiān)測(cè)的太陽輻射強(qiáng)度變換.太陽輻射強(qiáng)度隨天氣和時(shí)間變化較大，輻射強(qiáng)度越大，可提供集熱器內(nèi)低溫水的熱量越多，由于蓄熱水箱的存在，可使提供熱泵低溫?zé)崃康乃疁夭恢伦兓瘎×遥岣吡藷岜醚h(huán)的穩(wěn)定性.圖3為采暖期熱泵高低溫側(cè)進(jìn)出水溫度隨時(shí)間的變化.需要注意的是3.16日起，隨著天氣的轉(zhuǎn)暖，室內(nèi)供熱負(fù)荷減少，太陽輻射強(qiáng)度增加，集熱器中水溫升高，此時(shí)可以用集熱器的水直接供熱，熱泵可以停用.熱泵的啟?？捎捎脩糇孕袥Q定.

圖2　采暖期太陽能輻射強(qiáng)度變化

圖3　熱泵高低溫側(cè)供回水溫度變化

3　熱經(jīng)濟(jì)優(yōu)化

3.1熱泵模型

考慮熱泵為不可逆變溫?zé)嵩礋岜醚h(huán)，工質(zhì)與熱源換熱為逆流式，由傳熱學(xué)效能-傳熱單元數(shù)法(ε-NTU)可得冷凝器、蒸發(fā)器側(cè)的換熱量分別為：

Qc=Ccεc(Tc-Tci)

(1)

Qe=Ceεe(Tei-Te)

(2)

εc=1-exp(-NTUc)

(3)

εe=1-exp(-NTUe)

(4)

NTUc=KcAc/Cc

(5)

NTUe=KeAe/Ce

(6)

其中：Cc、Ce為高、低溫側(cè)熱源的熱容(質(zhì)量流率與定壓比熱的乘積)；εc、εe分別為冷凝器、蒸發(fā)器側(cè)換熱器的效能；Tci、Tco、Tei、Teo分別為冷凝器、蒸發(fā)器的進(jìn)出口溫度；NTUc、NTUe分別為冷凝器、蒸發(fā)器側(cè)傳熱單元數(shù)；Kc、Ke、Ac、Ae分別為冷凝器、蒸發(fā)器側(cè)的傳熱系數(shù)和傳熱面積.考慮到熱泵循環(huán)中存在耗散效應(yīng)等不可逆因素，引入不可逆因子：

I=ΔSc/ΔSe

(7)

ΔSc、ΔSe分別為高溫側(cè)流出的熵流量和低溫側(cè)流入的熵流量，由熱力學(xué)第二定律和式(7)：

Qc/Qe=ITc/Te

(8)

忽略熱漏損失，由式(1)～(8)得：

(9)

3.2集熱器模型

太陽能集熱器主要分為平板型和真空管型兩大類.本文以真空管型太陽能集熱器為研究對(duì)象.某時(shí)間從集熱器得到的有用能量是吸熱板吸收的太陽能量與散失到周圍環(huán)境的能量差可由公式[7]得

(10)

其中：Qj為集熱器收集的有用能量(W)；Aj為集熱器面積(m2)；αs集熱器表面吸收率；Ij為集熱器單位面積上的太陽能輻射率(W/m2)；α集熱器表面對(duì)流換熱系數(shù)，W/(m2·K)；Tin、Tout分別為集熱器進(jìn)水與出水溫度，K；Ta為周圍空氣溫度，K；σb黑體輻射常數(shù)，5.67×10-8W/(m2·K4)；Tsky天空溫度，K；αsky集熱器表面對(duì)天空輻射的吸收率；εs集熱器表面的發(fā)射率.

3.3地板供暖模型

地板傳熱量包括兩部分：對(duì)流傳熱量和輻射傳熱量.對(duì)流傳熱量與輻射傳熱量分別由式(11)和(12)確定[8].

qd=2.17(tp-tn)1.31

(11)

(12)

其中：qd為地板對(duì)流換熱量(W/m2)；qf為輻射換熱量(W/m2)；Tp為地板輻射面平均溫度(K)；Tf為非加熱表面平均溫度(K)；tn為采暖室內(nèi)設(shè)計(jì)溫度(℃).

由式(11)可知，影響地板傳熱量的參數(shù)有Tp和Tf.其中是tn是設(shè)計(jì)參數(shù)，由設(shè)計(jì)確定，在計(jì)算中可認(rèn)為是定值.則實(shí)際計(jì)算中影響地板傳熱量的參數(shù)只有Tp(本文近似為冷凝器進(jìn)、出口的算數(shù)平均溫度)和Tf，根據(jù)文獻(xiàn)[8]有：

tf=tn-1.1

(13)

為方便計(jì)算，借助局部阻力當(dāng)量長度的轉(zhuǎn)換思想，將輻射換熱量折算成對(duì)流換熱量，擬合公式得：

qz=4.298(tp-tn)1.31

(14)

擬合后的計(jì)算值與原值變化(表3).由表3可見，擬合后的總換熱量計(jì)算值與未擬合的對(duì)流、輻射換熱量之和的偏差≤3.6%.

則地板供暖系統(tǒng)的散熱可以用式(15)計(jì)算，

Qz=4.298(tp-tn)1.31Ad

(15)

Ad為地板供暖的散熱面積(m2).

表3　不同供回水平均溫度下?lián)Q熱量擬合計(jì)算表

3.4約束條件

在供熱負(fù)荷一定的情況下，copmax的表達(dá)式f：

copmax=f(Ce、Cc、Ae、Aj、Ad、I、Ij、T0、UL)

(16)

若以冷凝器與蒸發(fā)器面積來表征熱泵的初投資，以集熱器面積、地板散熱面積分別表征集熱器和地板采暖系統(tǒng)的初投資，則總初投資M可以表示為：

M=mjAj+mcAc+meAe+mdAd

(17)

其中：mj、mc、me、md分別為單位面積集熱器、冷凝器、蒸發(fā)器和地板采暖的價(jià)格.

3.5求解方法

由于cop的表達(dá)式比較繁瑣，同時(shí)還存在關(guān)于Ac、Ae、Af的非線性問題，因此本文采用數(shù)值計(jì)算的方法，編制計(jì)算機(jī)程序，使As、Ac、Ae、Af在滿足式(17)的情況下，分別以一定的步長變化，并求出對(duì)應(yīng)的cop值，從而找出當(dāng)cop值達(dá)到最大值時(shí)，對(duì)應(yīng)的As、Ac、Ae、Af值.

上述理論模型可以用于研究在一定的初投資M條件下，系統(tǒng)內(nèi)部組件的參數(shù)匹配問題.而當(dāng)系統(tǒng)的初投資M發(fā)生變化時(shí)，其系統(tǒng)最大copmax值及對(duì)應(yīng)的As、Ac、Ae、Af值也相應(yīng)發(fā)生變化.易知，當(dāng)初投資M增加時(shí)，As、Ac、As、Af及copmax也相應(yīng)增加，使得系統(tǒng)的運(yùn)行費(fèi)用降低.因此，為了求出系統(tǒng)的最佳設(shè)計(jì)參數(shù)，使得系統(tǒng)在運(yùn)行年限內(nèi)每年折合初投資與運(yùn)行費(fèi)用之和L(M)為目標(biāo)函數(shù)，求出當(dāng)L(M)取得最優(yōu)(最小)時(shí)，系統(tǒng)的最佳設(shè)計(jì)參數(shù)As、Ac、Ae、Af值.若以熱泵的耗電量來表征系統(tǒng)的運(yùn)行費(fèi)用，可得：

(18)

其中：nx為系統(tǒng)運(yùn)行年限；γ為系統(tǒng)運(yùn)行其他費(fèi)用(包括循環(huán)水泵、風(fēng)機(jī)及維修費(fèi)用等)占熱泵耗電量的百分比；D為采暖天數(shù)；E為電費(fèi).

4　結(jié)果與分析

數(shù)值計(jì)算參數(shù)如下：Qc=10kW，Cc=1 260W/K，Ce=1 050W/K，Ke=Kc=1 500W/ (m2·K)，I=1.1，Ij=200W/m2，αs=0.92，α=2.5W/(m2·K)，εs=0.07，Ts=278K，Ta=258K，Tsky=0K，mj=400 ￥/m2，mc= me=1 200 ￥/m2，md=90 ￥/m2，E=0.55￥/kWh，nx=20a，γ=10﹪，D=174d.計(jì)算得到冷凝器、蒸發(fā)器、地板供暖最優(yōu)面積、集熱器最優(yōu)面積、最佳冷凝溫度、最佳蒸發(fā)溫度、系統(tǒng)最大供熱系數(shù)、最小耗功率、年折合初投資及運(yùn)行費(fèi)用隨初投資的變化關(guān)系，進(jìn)而通過系統(tǒng)年總費(fèi)用最低的優(yōu)化目標(biāo)確定最佳初投資及各換熱裝置的最佳面積. 由圖4～6可以看出隨初投資M的增加，集熱器、蒸發(fā)器、冷凝器與地板供暖的面積都相應(yīng)增加，而且基本呈線性增加的趨勢(shì)，只是各自的增長速度有所不同，這主要是由于各個(gè)部分的傳熱機(jī)理、傳熱系數(shù)、熱損失系數(shù)以及單位面積的價(jià)格不同等因素造成.

圖4　冷凝器與蒸發(fā)器最優(yōu)面積與初投資M的關(guān)系

圖5　地板供暖最優(yōu)面積與初投資M的關(guān)系

圖6　集熱器最優(yōu)面積與初投資M的關(guān)系

圖7顯示隨著初投資M的增加，系統(tǒng)的最佳蒸發(fā)溫度逐漸上升，而冷凝溫度逐漸下降.因?yàn)榈責(zé)岜P管、蒸發(fā)器與冷凝器面積隨M變大，使得換熱

器的換熱效率增加，在供熱率和室內(nèi)設(shè)計(jì)溫度一定的條件下，冷凝器進(jìn)口溫度下降，冷凝換熱的最大溫差變小，從而使得冷凝溫度變小.蒸發(fā)溫度升高的主要原因是太陽集熱器面積的增加，使得系統(tǒng)的集熱量增加，進(jìn)而提高了蒸發(fā)器側(cè)的平均水溫.圖8反映了系統(tǒng)最大供熱性能系數(shù)與最小耗功率隨初投資M的變化情況，由圖中可以看出，隨M的增加系統(tǒng)消耗的功率減小，進(jìn)而使得系統(tǒng)的供熱系數(shù)增加，這主要是隨著M的增加使得蒸發(fā)溫度升高，而冷凝溫度降低共同作用的結(jié)果. 由圖9、10可以看出，隨著系統(tǒng)初投資M的增加，每年折和的初投資增加，而系統(tǒng)的運(yùn)行費(fèi)用逐漸下降，這主要因?yàn)橄到y(tǒng)的耗功率隨M下降，由圖10可知，在上述的算例計(jì)算參數(shù)下，系統(tǒng)的最佳投資為M=4.6萬元，此時(shí)的折算費(fèi)用之和最小，相應(yīng)的copmax=5.2，最佳集熱器面積Aj=55.275m2，冷凝器面積Ac=1.831m2，蒸發(fā)器面積Ae=1.715m2，地板供暖面積Ad=204.72m2，熱泵功率P=1 922.13W.

圖7　最佳冷凝溫度和蒸發(fā)溫度值與初投資M的關(guān)系

圖8　系統(tǒng)最大供熱系數(shù)、最小耗功率隨初投資變化

圖9　系統(tǒng)每年折合初投資、運(yùn)行費(fèi)用隨初投資變化

圖10　系統(tǒng)每年總費(fèi)用與初投資M的關(guān)系

5　結(jié)　語

太陽能熱泵地板供暖系統(tǒng)是一種節(jié)能、經(jīng)濟(jì)、環(huán)保的新型供暖系統(tǒng)，為了合理有效的利用該系統(tǒng)，使其特點(diǎn)能夠得到充分發(fā)揮.本文測(cè)試了中國北方城市某住宅內(nèi)冬季使用該系統(tǒng)供熱的運(yùn)行情

況，結(jié)果顯示系統(tǒng)滿足該用戶的冬季供暖要求.最后在供熱率一定的條件下，利用有限時(shí)間熱力學(xué)的方法，建立了太陽能熱泵地板供暖系統(tǒng)的數(shù)學(xué)模型，用數(shù)值方法對(duì)系統(tǒng)進(jìn)行了熱經(jīng)濟(jì)性分析，得出了對(duì)應(yīng)供熱率的最佳初投資及其相應(yīng)優(yōu)化的設(shè)計(jì)參數(shù)，為太陽能熱泵地板供暖系統(tǒng)的設(shè)計(jì)起到一定的參考作用.

[1]宋偉. 大連市某污水源熱泵-太陽能熱水供應(yīng)系統(tǒng)運(yùn)行特性分析[D]. 大連: 大連理工大學(xué), 2014.

[2]張龍. 青海地區(qū)太陽能地下水源熱泵復(fù)合系統(tǒng)供熱特性研究[D]. 邯鄲. 河北工程大學(xué), 2014.

[3]蘇婉昀. 太陽能跨季節(jié)貯熱-地源熱泵熱水系統(tǒng)在住宅中的應(yīng)用研究[D]. 西安: 長安大學(xué), 2014.

[4]韓宗偉, 鄭茂余, 劉威, 等. 嚴(yán)寒地區(qū)太陽能-土壤源熱泵相變蓄熱供暖系統(tǒng)[J]. 太陽能學(xué)報(bào), 2006, 27(12): 1214-1218.

[5]BEJANA.Powergenerationandrefrigerati-onmodelswithheattransferirreversibilities[J].HeatTransferSocJpn, 1994, 33(128): 68-75.

[6]CHENLG,SUNFR,WUC, et al.TheoreticalanalysisoftheperformanceofaregenerativeclosedBraytoncyclewithinternalirreversibilities[J].EnergyConversM-gmt, 1997, 3(9) :871-877.

[7]羅運(yùn)俊, 李元哲, 趙承龍. 太陽熱水器原理、制造與施工[M]. 北京: 化學(xué)工業(yè)出版社, 2005. 56-58.

[8]WANGSK．Handbookofconditionandrefrigeration[M].NewYork:McGraw-HillCompanies, 2000.

Test and thermo-economic analysis for solar source heat pump and floor heating system

WANG Chen1, ZHAO Qi2

(1. Design and Research Institute, Changchun Institute of Technology, Changchun 130012, China;2.SchoolofEnergyandPower,ChangchunInstituteofTechnology,Changchun130012,China)

AsourcesolarheatpumpandfloorheatingdeviceissetinaresidentialroomtotesttheperformanceofthesysteminwinterinnortherncityofChina.Theresultsshowedthesystemwellmeetthewinterheatingrequirementsoftheuser.Inthispaper,amathematicalmodelofasolarsourceheatpumpandfloorheatingsystemwasestablishedandthethermo-economicperformanceofthesystemwascarriedoutbasedonfinite-timethermodynamics.Underthegivenheatingloadandthefirstinvestmentcondition,thethermo-economicobjectivefunctionforthissystemwasobtained.Theoptimalcoefficientofheatingperformanceanddesignparametersformaximizedobjectivefunctionwerealsoobtained.Theinfluenceofthefirstinvestmentforcoefficientofheatingperformanceanddesignparameterswasinvestigated,andtheoptimalfirstinvestmentandcorrespondingdesignparametersforagivenheatingloadwereobtained.

finitetimethermodynamics;solarsourceheatpumpandfloorheatingsystem;thermo-economic;objectivefunction

2015-11-26.

吉林省教育廳科學(xué)技術(shù)研究項(xiàng)目(120150044)

王琛(1982-)，女，工程師，研究方向：建筑環(huán)境與設(shè)備工程.

TU832.1

A

1672-0946(2016)04-0492-06