徐天碩，王 劍，金佳煜

(杭州電子科技大學自動化院，浙江 杭州 310018)

0 引 言

太陽能-熱泵復合熱水系統(tǒng)(簡稱S-H系統(tǒng))既發(fā)揮了太陽能集熱系統(tǒng)低成本、低污染的優(yōu)點，又可在太陽能不足的情況下利用空氣源熱泵系統(tǒng)制取熱水，具有結(jié)構(gòu)簡單、熱效率高、環(huán)境適應能力強等優(yōu)點。根據(jù)運行環(huán)境參數(shù)和用戶需求，對參數(shù)進行優(yōu)化匹配實現(xiàn)系統(tǒng)高效穩(wěn)定運行，達到性能和費用的最優(yōu)，是當前S-H系統(tǒng)的研究熱點。Z.M.Amin等[1]對新加坡地區(qū)S-H系統(tǒng)進行了應用性研究分析，通過建立數(shù)學模型，得出在太陽能系統(tǒng)與熱泵系統(tǒng)合理配合時，系統(tǒng)集熱效率可達到80%～90%，系統(tǒng)能效比(Coefficient of Performance，COP)達到8.0以上。J.V.Anderson等[2]采用f-chart方法，確定太陽能保證率和太陽能系統(tǒng)承擔的負荷量，從而確定熱泵系統(tǒng)應承擔的熱水負荷量，實現(xiàn)了S-H系統(tǒng)的最優(yōu)化設計。G.Panaras等[3]研究太陽能-熱泵熱水系統(tǒng)的性能，發(fā)現(xiàn)熱泵溫度設定點是系統(tǒng)性能的一個重要參數(shù)，設定點的值越高，熱泵的效率就越低。劉雨曦[4]以系統(tǒng)全生命周期成本為評價函數(shù)，分析各參數(shù)變化對系統(tǒng)的影響程度，通過對S-H系統(tǒng)多個參數(shù)進行同步優(yōu)化匹配，得出適用于西昌地區(qū)的S-H系統(tǒng)優(yōu)化參數(shù)模型。王崗等[5]通過研究各運行模式下系統(tǒng)運行特性受環(huán)境參數(shù)的影響，提出基于環(huán)境參數(shù)、氣象參數(shù)、硬件參數(shù)的太陽能熱源與輔助熱源最優(yōu)匹配的研究方向。丁勇等[6]對重慶地區(qū)S-H系統(tǒng)的關鍵因素進行模型化研究分析，得出適合重慶的熱水系統(tǒng)性能參數(shù)，提高了能源利用率。目前對S-H系統(tǒng)的研究主要是對系統(tǒng)中單個參數(shù)進行優(yōu)化，或者對系統(tǒng)多個參數(shù)采用單因素分析法進行研究，沒有深入分析各參數(shù)之間的耦合現(xiàn)象以及多參數(shù)交互作用對系統(tǒng)性能的影響。另外，大部分研究主要通過參數(shù)優(yōu)化匹配降低系統(tǒng)的成本和運行費用，對系統(tǒng)本身運行性能指標的研究較少。本文以某高校學生公寓太陽能-熱泵熱水工程為對象，分析主要設計參數(shù)對系統(tǒng)性能及費用的影響，針對系統(tǒng)優(yōu)化匹配問題展開研究。

1 S-H系統(tǒng)的基本結(jié)構(gòu)

S-H系統(tǒng)由太陽能集熱器、空氣源熱泵機組、儲熱水箱、用戶供給側(cè)組成。本文中太陽能集熱器和空氣源熱泵采取并聯(lián)式結(jié)構(gòu)，其制造的熱水共同匯集于儲熱水箱后，統(tǒng)一通過供水管路供給用戶熱水。S-H系統(tǒng)的基本結(jié)構(gòu)如圖1所示。

圖1 S-H系統(tǒng)的基本結(jié)構(gòu)

2 S-H系統(tǒng)的主要技術指標

2.1 S-H系統(tǒng)的能效指標

S-H系統(tǒng)的評價指標主要包括太陽能系統(tǒng)集熱系數(shù)、熱泵系統(tǒng)制熱系數(shù)、熱泵能源轉(zhuǎn)換效率之比(簡稱能效比或性能系數(shù))COPH、系統(tǒng)能效比COPS-H等。其中系統(tǒng)能效比COPS-H能夠較全面地評價S-H系統(tǒng)中輸入-產(chǎn)出的綜合能效，本文選用COPS-H作為研究S-H系統(tǒng)的運行性能指標。

S-H系統(tǒng)的能效比COPS-H為：

(1)

式中，QW為熱水需求量，單位為kg；CW為水的定壓比熱容，單位為kJ/(kg·℃)；Tend為供熱水溫，單位為℃；TL為冷水溫度，單位為℃；NS為太陽能集熱系統(tǒng)耗電量，單位為kJ；NH為熱泵制熱系統(tǒng)耗電量，單位為kJ。

2.2 S-H系統(tǒng)的經(jīng)濟指標

S-H系統(tǒng)的經(jīng)濟指標[7]中包含系統(tǒng)投資費用和系統(tǒng)運行費用。本文采用動態(tài)費用法將系統(tǒng)初投資折算成月費用，并與月運行總費用相加得出月費用總值：

(2)

式中，Zm為按動態(tài)法計算的月費用，單位為元；i為利率，一般取8%；m為系統(tǒng)使用年限，一般取15；K為初投資，S-H系統(tǒng)中集熱器按工程經(jīng)驗取1 500元/m2；C為月運行總費用，主要為S-H系統(tǒng)的運行費用和維護費用，運行費用可通過計算得出系統(tǒng)消耗的電能，電價按0.82元/(kW·h)計算，年維護費按系統(tǒng)投資5%計算。

3 設計參數(shù)對S-H系統(tǒng)性能的影響

本文以浙江地區(qū)某學生公寓S-H系統(tǒng)為例，基于逐月日平均氣象數(shù)據(jù)，研究集熱器傾角、供熱溫度、設計太陽能保證率(本文不涉及系統(tǒng)運行時的實際太陽能保證率，下文簡稱太陽能保證率)等設計參數(shù)變動時，S-H系統(tǒng)制熱性能和月費用的變化情況。該系統(tǒng)用水定額QW=4.8×104L/d。由于集熱面積有限，故采用熱泵為主、太陽能為輔的熱水系統(tǒng)，希望在運行中力求最大化地利用太陽能，實現(xiàn)節(jié)能降耗。

3.1 集熱器傾角β對S-H系統(tǒng)性能的影響

在研究集熱器傾角對S-H系統(tǒng)影響時，太陽能集熱面積Ac為：

(3)

式中，Ti為年平均冷水溫度，根據(jù)實際情況取10 ℃；f為太陽能保證率；JT為年平均太陽輻射量，單位為kJ/m2；ηc為年平均集熱效率，一般取0.4～0.6；ηL為管路熱損率，一般取0.2～0.3。

太陽能集熱量Qs為：

QS=QWCW(Tlend-TL)=JTdAcηcd(1-ηL)

(4)

式中，Tlend為太陽能集熱器單獨工作時可提升的集熱水箱的水溫，單位為℃；ηcd為逐月的日平均集熱效率；JTd為逐月氣象條件下日平均太陽輻射量，單位為kJ/m2，JTd與集熱器傾角β之間的關系可由Hay模型[8]求得：

(5)

式中，Hb為水平面上月平均日直射輻射量，單位為kJ/m2；Rb為傾斜面與水平面月平均日直接輻射總量的比值；Hd為水平面上月平均日散射輻射量，單位為kJ/m2；Rd為傾斜面與水平面月平均日散射輻射量的比值；ρ為地表反射率，一般取0.2。

在已知環(huán)境溫度情況下，擬合環(huán)境平均溫度Ta和冷水溫度TL的關系如下：

TL=4.717e0.041Ta

(6)

在本文的熱水系統(tǒng)中，熱泵額定制熱量NH為40.5 kW/h，其運行功率除了跟熱泵機組的系統(tǒng)匹配性相關，還受到環(huán)境溫度和熱泵進口溫度的影響，在不考慮水箱熱損的情況下，根據(jù)經(jīng)驗公式[9]可得：

P=b1Ta+b2Tend+b3

(7)

在環(huán)境平均溫度Ta為24 ℃、供熱溫度Tend為55 ℃、太陽能保證率f為20%時，集熱器傾角變化對傾斜表面月平均日輻射量和太陽能集熱停止時水箱溫度的影響如圖2(a)所示；不同傾斜角下，系統(tǒng)性能系數(shù)COPS-H和月費用Zm如圖2(b)所示。

圖2 集熱器傾角β對S-H系統(tǒng)的影響

由圖2(a)可以看出：太陽能集熱器傾斜表面的月平均日輻射量和月平均日太陽能集熱終止時的水箱水溫隨集熱器傾角變化的趨勢基本相同，傾角為30°左右時，兩者都達到最大。由圖2(b)可以看出：傾角為30°左右時，性能評價指標COPS-H和經(jīng)濟性評價指標Zm也達到最優(yōu)。因此，在進行S-H系統(tǒng)優(yōu)化匹配問題后續(xù)研究中，集熱器傾角選取30°。

3.2 供熱溫度Tend對S-H系統(tǒng)的影響

由式(3)可知，集熱面積主要和系統(tǒng)供熱溫度Tend與太陽能保證率f有關。在研究供熱溫度Tend對太陽能熱泵系統(tǒng)運行性能指標和經(jīng)濟性指標的影響時，集熱器傾角為30°，其它相關參數(shù)和3.1節(jié)一致。供熱溫度變化對S-H系統(tǒng)的影響如圖3所示。

圖3 供熱溫度Tend對S-H系統(tǒng)的影響

由圖3(a)可以看出：當供熱溫度從50 ℃～60 ℃逐漸變化時，太陽能集熱水箱終溫逐漸升高，其與熱泵運行初始溫度的溫差也逐漸增大，溫差決定了熱泵需提供的熱量，因此熱泵運行的時間也逐漸增加。圖3(b)可以看出：在S-H系統(tǒng)運行中，供熱溫度的改變與系統(tǒng)性能系數(shù)呈反比關系。就經(jīng)濟性而言，隨著供熱溫度的升高，集熱器集熱面積和所需熱泵制熱量也隨之增大，月費用與供熱溫度成正比關系。

3.3 太陽能保證率f對S-H系統(tǒng)的影響

由于受集熱面積的約束，太陽能保證率不能取值過高，本文取20%～30%變化區(qū)間。集熱器傾角為30°，供熱溫度為Tend為55 ℃，其余參數(shù)與前文一致。仿真結(jié)果如圖4所示。

圖4 太陽能保證率f對S-H系統(tǒng)的影響

由圖4(a)可以看出：當太陽能保證率f增大時，集熱系統(tǒng)終溫Tlend也增大。由于供熱溫度不變，熱泵制熱水終溫與集熱終溫溫差逐漸減小，熱泵月平均日制熱量逐漸下降。圖4(b)可以看出：隨著熱泵制熱量減小，熱泵系統(tǒng)運行時間降低，性能指標COPS-H隨之增大。另一方面，系統(tǒng)月費用也隨保證率的增加呈上升趨勢，這是因為隨著保證率的增大，系統(tǒng)設計集熱面積隨之增大，從而使系統(tǒng)初投資增加，由上述分析可知：隨著保證率增大，熱泵系統(tǒng)運行能耗降低，運行費用減少。但分析發(fā)現(xiàn)，初投資和維護費用的增長對月費用的影響大于運行費用減少的影響，從而造成動態(tài)月計算費用逐漸增大。

4 S-H系統(tǒng)中設計參數(shù)的優(yōu)化匹配問題

4.1 優(yōu)化匹配問題

本文所討論的優(yōu)化匹配問題可描述為：在滿足用戶用水需求和集熱面積約束的前提下，通過合理選取設計參數(shù)，使系統(tǒng)的運行性能和費用達到最優(yōu)。選取系統(tǒng)性能系數(shù)COPS-H與月費用Zm作為目標函數(shù)，分別研究在單一優(yōu)化目標和復合優(yōu)化目標下的參數(shù)匹配問題，通過MATLAB軟件運用遺傳算法分別求解。

4.2 遺傳算法優(yōu)化

遺傳優(yōu)化算法以二進制編碼的若干個個體為初始點，其中每個個體由優(yōu)化變量(集熱器傾角、供熱溫度和太陽能保證率)來表示其特征。將目標函數(shù)(COPS-H和Zm)轉(zhuǎn)化為適應性函數(shù)，通過計算得到其適應度,并判斷是否達到優(yōu)化標準，若符合則輸出其對應的最優(yōu)解結(jié)束優(yōu)化；否則根據(jù)適應度依次進行選擇、交叉與變異運算。首先選取初始群體中適應度高的個體保留，接著通過交叉與變異操作得到新一代個體,使個體組成的群體得到更新,然后用適應性函數(shù)來計算新個體的適應度,直至滿足優(yōu)化準則或者達到最大迭代次數(shù)。

4.3 單目標函數(shù)優(yōu)化結(jié)果

4.3.1 系統(tǒng)能效比COPS-H

運用遺傳算法對性能指標系統(tǒng)能效比COPS-H的優(yōu)化結(jié)果如表1所示。

表1 系統(tǒng)能效比COPS--H優(yōu)化結(jié)果

由表1可以看出：集熱器傾角為33.1°，供熱溫度為50.6 ℃，太陽能保證率為28.1%時，系統(tǒng)的運行性能系數(shù)達到最大，為8.3。與優(yōu)化前的系統(tǒng)性能系數(shù)相比，COPS-H提升了13.70%。通過計算可知：優(yōu)化后的設計方案使月運行費用提高約4.24%。

4.3.2 月費用Zm

運用遺傳算法對月折算費用Zm優(yōu)化結(jié)果如表2所示。

由表2可以看出：集熱器傾角為31.5°，供熱溫度為50.1 ℃，太陽能保證率為20.4%時，系統(tǒng)月費用最低，為14 913元。與優(yōu)化前相比，優(yōu)化后的系統(tǒng)月費用Zm降低了12.35%；通過計算可知：此時COPS-H提高約2.39%?？梢妰?yōu)化后的方案使月費用降低較大，但對系統(tǒng)性能的優(yōu)化并不明顯。

對單一目標函數(shù)的優(yōu)化結(jié)果表明：在對某一目標函數(shù)的優(yōu)化時，往往造成另一目標函數(shù)的結(jié)果不理想。為平衡系統(tǒng)性能系數(shù)函數(shù)和月費用函數(shù)之間的矛盾，本文引入雙目標函數(shù)優(yōu)化方法。

4.4 雙目標函數(shù)優(yōu)化

對系統(tǒng)性能、費用雙目標函數(shù)同時進行優(yōu)化，并基于遺傳算法進行求解。目標函數(shù)如下：

(8)

式中，a1，a2分別為雙目標函數(shù)的加權系數(shù)，且a1+a2=1，本設計采用權重法中的容限分析法計算得出a1=0.460 33，a2=0.539 67；自變量x1，x2，x3分別對應本設計中的集熱器傾角、供熱溫度和太陽能保證率；遺傳算法默認求取函數(shù)最小值，此時F1(x)為系統(tǒng)性能系數(shù)函數(shù)的倒數(shù)；F2(x)為月費用函數(shù)。

在遺傳算法中，確定目標函數(shù)后，設置種群規(guī)模為300，交叉概率為0.4，變異概率為0.2，最大迭代次數(shù)為500。由于求解多目標函數(shù)優(yōu)化問題的過程就是尋找Pareto最優(yōu)解(非劣解、有效解)，而算法所獲得的非劣最優(yōu)解與實際的Pareto前端會存在一定的距離，因此需要多次調(diào)用遺傳算法求解，計算平均值，使算法所找到的Pareto前端更接近實際Pareto前端。運用遺傳算法求解出30組非劣解，其中2次遺傳算法迭代優(yōu)化結(jié)果如圖5所示。

圖5中，橫坐標表示迭代次數(shù)，縱坐標表示經(jīng)過平衡后的目標函數(shù)值，目標函數(shù)值越接近1表示優(yōu)化結(jié)果越好。30次優(yōu)化結(jié)果的平均值如表3所示。

表3 雙目標函數(shù)優(yōu)化結(jié)果

由表3可以看出：集熱器傾角為31.8°，供熱溫度為50.1 ℃，太陽能保證率為22.4%時，目標函數(shù)達到最優(yōu)。在此參數(shù)下，系統(tǒng)運行性能系數(shù)達到7.67，比原系統(tǒng)提高了5.48%，月折算費用為15 554元，降低了8.24%。

5 結(jié)束語

本文研究了集熱器傾角、供熱溫度和太陽能保證率等設計參數(shù)對S-H系統(tǒng)的系統(tǒng)性能和經(jīng)濟性的影響，并基于遺傳算法對系統(tǒng)優(yōu)化匹配問題進行求解。在多目標函數(shù)的優(yōu)化中，通過合理分配加權因子得出S-H系統(tǒng)的最佳設計參數(shù)，解決了系統(tǒng)性能指標和經(jīng)濟性指標之間的矛盾。后續(xù)研究主要集中在建立全流程動態(tài)S-H系統(tǒng)模型，針對遺傳算法收斂速度較慢等問題對算法進行改進，并運用改進后的遺傳算法對S-H系統(tǒng)優(yōu)化匹配問題求解，驗證優(yōu)化方案的可行性。