高校供熱系統(tǒng)調(diào)控策略優(yōu)化研究

王睿鑫， 周志剛， 荊 強(qiáng)， 諶薛蛟

(哈爾濱工業(yè)大學(xué) 建筑學(xué)院，黑龍江哈爾濱150006)

1 概述

目前，我國(guó)北方城鎮(zhèn)集中供熱面積達(dá)到91.4×108m2，建筑供暖能耗占建筑總能耗的24.1%[1]。高校建筑由于建筑類型雜，體量大，在建筑供暖能耗中占有相當(dāng)大的比重。據(jù)統(tǒng)計(jì)，高等院校的人均能源消耗是普通居民的4倍，北方高校供暖能耗已達(dá)學(xué)校整體能耗的50%以上，位居高校能耗首位[2]。因此，供熱系統(tǒng)的節(jié)能優(yōu)化是高等院校節(jié)能工作的重點(diǎn)。高校供熱系統(tǒng)的運(yùn)行控制是節(jié)能的關(guān)鍵，在尋找合理的調(diào)控策略過(guò)程中，許多學(xué)者以負(fù)荷優(yōu)化為目標(biāo)研究了供熱系統(tǒng)的部分負(fù)荷特性、影響因素[3-4]以及負(fù)荷預(yù)測(cè)方法[5-6]。部分學(xué)者則以多種能源供給為前提，研究了供熱參數(shù)調(diào)控優(yōu)化算法[7-8]。

本文針對(duì)北方高校，以穩(wěn)態(tài)傳熱模型為基礎(chǔ)，開(kāi)發(fā)包含熱源、二級(jí)管網(wǎng)、建筑在內(nèi)的供熱系統(tǒng)仿真模型。結(jié)合校園建筑分時(shí)段使用的特點(diǎn)，提出供熱系統(tǒng)啟停、建筑使用期間室溫調(diào)控優(yōu)化方法。采用仿真模型，對(duì)優(yōu)化后的調(diào)控策略效果進(jìn)行模擬分析。

2 高校供熱系統(tǒng)仿真模型

2.1 數(shù)學(xué)模型

供熱系統(tǒng)數(shù)學(xué)模型架構(gòu)見(jiàn)圖1。由圖1可知，供熱系統(tǒng)數(shù)學(xué)模型涵蓋了熱力站模型、二級(jí)管網(wǎng)模型、用戶模型。

圖1 供熱系統(tǒng)數(shù)學(xué)模型架構(gòu)

① 熱力站模型

換熱器一級(jí)側(cè)熱平衡關(guān)系式為：

式中m1——換熱器一級(jí)側(cè)水的質(zhì)量，kg

cp,1——換熱器一級(jí)側(cè)水的比定壓熱容，J/(kg·K)

θ1,out——換熱器一級(jí)側(cè)出水溫度，℃

t——時(shí)間，s

qm,1——換熱器一級(jí)水的質(zhì)量流量，kg/s

θ1,in——換熱器一級(jí)側(cè)進(jìn)水溫度，℃

K——換熱器的傳熱系數(shù)，W/(m2·K)

A——換熱器的換熱面積，m2

Δtm——換熱對(duì)數(shù)平均溫差，℃

換熱器二級(jí)側(cè)熱平衡關(guān)系為：

式中m2——換熱器二級(jí)側(cè)水的質(zhì)量，kg

cp,2——換熱器二級(jí)側(cè)水的比定壓熱容，J/(kg·K)

θ2,out——換熱器二級(jí)側(cè)出水溫度，℃

qm,2——換熱器二級(jí)水的質(zhì)量流量，kg/s

θ2,in——換熱器二級(jí)側(cè)進(jìn)水溫度，℃

② 二級(jí)管網(wǎng)模型

由熱平衡關(guān)系，得到二級(jí)管網(wǎng)模型為[9]：

式中C——管網(wǎng)的熱容，J/K

θout——管網(wǎng)的出水溫度，℃

cp——管網(wǎng)中水的比定壓熱容，J/(kg·K)

qm——管網(wǎng)中水的質(zhì)量流量，kg/s

θin——管網(wǎng)的進(jìn)水溫度，℃

Kp——管中水與土壤間的傳熱系數(shù)，W/(m2·K)

Ap——管子的散熱面積，m2

θsoil——土壤溫度，℃

③ 熱用戶模型

由熱平衡關(guān)系，得到散熱器模型為：

式中Cs——散熱器熱容，J/K

θs,out——散熱器的出水溫度，℃

cp,s——散熱器內(nèi)水的比定壓熱容，J/(kg·K)

qm,s——散熱器中水的質(zhì)量流量，kg/s

θs,in——散熱器的進(jìn)水溫度，℃

Ks——散熱器的傳熱系數(shù)，W/(m2·K)

As——散熱器的傳熱面積，m2

θr——室內(nèi)溫度，℃

圍護(hù)結(jié)構(gòu)的傳熱可以近似為一維導(dǎo)熱，圍護(hù)結(jié)構(gòu)模型為：

hw,outAw,out(θw,out-θo)

式中Cw——墻體的熱容，J/K

θw,in——墻體內(nèi)表面溫度，℃

hw,in——墻體內(nèi)壁表面?zhèn)鳠嵯禂?shù)，W/(m2·K)

Aw,in——墻體內(nèi)壁傳熱面積，m2

hw,out——墻體外壁表面?zhèn)鳠嵯禂?shù)，W/(m2·K)

Aw,out——墻體外壁傳熱面積，m2

θw,out——墻體外表面溫度，℃

θo——室外溫度，℃

忽略太陽(yáng)輻射得熱量、室內(nèi)人員及設(shè)備散熱，室溫模型為[10]：

hw,inAw,in(θr-θw,in)-KcAc(θr-θo)-

0.278ρa(bǔ)cp,aVr(θr-θo)

式中Cr——房間的熱容，J/K

Kc——外窗傳熱系數(shù)，W/(m2·K)

Ac——外窗面積，m2

ρa(bǔ)——空氣密度，kg/m3

cp,a——空氣的比定壓熱容，kJ/(kg·K)

Vr——房間的容積，m3

2.2 仿真模型

根據(jù)數(shù)學(xué)模型，采用TRNSYS構(gòu)建校園供熱系統(tǒng)仿真模型。模擬對(duì)象為哈爾濱某高校供熱系統(tǒng)，室內(nèi)采用散熱器供暖系統(tǒng)，總供熱面積為7.9×104m2，總設(shè)計(jì)熱負(fù)荷為7 745 kW。供暖室內(nèi)設(shè)計(jì)溫度為18 ℃，供暖室外計(jì)算溫度為-24.2 ℃，室外的平均風(fēng)速為3.2 m/s。建筑熱力入口安裝分時(shí)供熱控制器、電動(dòng)調(diào)節(jié)閥、熱量表，供回水管安裝溫度、壓力傳感器。熱力站采用質(zhì)-量并調(diào)方式，二級(jí)循環(huán)泵為變頻泵。供熱系統(tǒng)平面布置見(jiàn)圖2。

圖2 供熱系統(tǒng)平面布置

TRNSYS仿真系統(tǒng)見(jiàn)圖3。室外氣象參數(shù)來(lái)自國(guó)家氣象局網(wǎng)站，由室外氣象模塊調(diào)用對(duì)應(yīng)模擬時(shí)間段內(nèi)的氣象參數(shù)。

圖3 TRNSYS仿真系統(tǒng)(軟件截圖)

TRNSYS仿真系統(tǒng)主要組件包括建筑物模塊、室外氣象模塊、時(shí)間控制模塊、計(jì)算模塊。供熱系統(tǒng)包含7幢建筑(見(jiàn)圖2)，建筑幾何模型(見(jiàn)圖4)按實(shí)際結(jié)構(gòu)尺寸構(gòu)建，利用SketchUp軟件建立。建筑幾何模型以特定文件形式儲(chǔ)存，可被TRNSYS軟件中的建筑物模塊調(diào)用。圍護(hù)結(jié)構(gòu)參數(shù)在建筑物模塊中進(jìn)行設(shè)置，圍護(hù)結(jié)構(gòu)參數(shù)見(jiàn)表1。外窗為雙層玻璃塑鋼窗，傳熱系數(shù)為2.0 W/(m2·K)。

圖4 建筑幾何模型

表1 建筑物圍護(hù)結(jié)構(gòu)參數(shù)

2.3 模型驗(yàn)證

利用TRNSYS仿真模型模擬各建筑2020年1月20日—23日的室內(nèi)溫度，以相同時(shí)段的實(shí)測(cè)數(shù)據(jù)進(jìn)行驗(yàn)證。各建筑室內(nèi)溫度仿真值與實(shí)測(cè)值的最大誤差、標(biāo)準(zhǔn)差見(jiàn)表2，各建筑模擬日供暖能耗的仿真值與實(shí)測(cè)值的相對(duì)誤差見(jiàn)表3。

表2 各建筑室內(nèi)溫度仿真值與實(shí)測(cè)值的最大誤差、 標(biāo)準(zhǔn)差

由表2可知，對(duì)于室內(nèi)溫度，仿真值與實(shí)測(cè)值的最大誤差不超過(guò)1.3 ℃，標(biāo)準(zhǔn)差不超過(guò)0.6 ℃。由表3可知，對(duì)于供暖能耗，仿真值與實(shí)測(cè)值的最大相對(duì)誤差為5.19%，平均相對(duì)誤差不超過(guò)3%。這表明仿真結(jié)果可信。

表3 各建筑模擬日供暖能耗的仿真值與實(shí)測(cè)值的相對(duì)誤差 %

3 調(diào)控策略優(yōu)化方法

3.1 啟停優(yōu)化

高校建筑的運(yùn)行時(shí)間為8:00—22:00，由此制定的室內(nèi)目標(biāo)溫度分時(shí)控制原則見(jiàn)圖5。由圖5可知，分時(shí)控制原則是在夜間僅滿足室內(nèi)供暖系統(tǒng)防凍需要(室內(nèi)最低溫度5 ℃)，并保證建筑在開(kāi)始使用時(shí)室內(nèi)溫度滿足要求。實(shí)際上，得益于建筑熱惰性的作用，停熱期間室內(nèi)溫度為緩慢下降過(guò)程，這也為供熱開(kāi)始后室內(nèi)溫度迅速回到設(shè)計(jì)溫度提供了條件。

圖5 室內(nèi)目標(biāo)溫度分時(shí)控制原則

常見(jiàn)的啟停策略往往根據(jù)經(jīng)驗(yàn)，在建筑開(kāi)始使用前的3～5 h供熱，在建筑結(jié)束使用時(shí)停熱。但這種啟停策略缺少理論依據(jù)，導(dǎo)致在建筑使用時(shí)間外，也保持了較高的室內(nèi)溫度。

筆者利用仿真模型，模擬了2019—2020年供暖期，不同室外平均溫度(指供熱開(kāi)始至建筑開(kāi)始使用期間的室外平均溫度)下，不同供熱開(kāi)始時(shí)間的室內(nèi)溫升(指供熱開(kāi)始時(shí)間至8:00的室內(nèi)溫升)，以及不同室外平均溫度(指停熱至建筑停止使用期間的室外平均溫度)下，不同停熱時(shí)間的室內(nèi)溫降(指停熱時(shí)間至22:00的室內(nèi)溫降)，分別見(jiàn)圖6、7。由圖6可知，相同室外平均溫度下，供熱開(kāi)始時(shí)間越早，室內(nèi)溫升越大。由圖7可知，相同室外平均溫度下，停熱時(shí)間越晚，室內(nèi)溫降越小。因此，對(duì)于具體項(xiàng)目，可根據(jù)仿真結(jié)果選取供熱開(kāi)始時(shí)間、停熱時(shí)間，對(duì)啟停策略進(jìn)行優(yōu)化。

圖6 不同室外平均溫度下不同供熱開(kāi)始時(shí)間的 室內(nèi)溫升

圖7 不同室外平均溫度下不同停熱時(shí)間的室內(nèi)溫降

根據(jù)仿真結(jié)果，對(duì)室內(nèi)溫升、溫降與室外平均溫度的關(guān)系式分別進(jìn)行擬合。

① 室內(nèi)溫升

供熱開(kāi)始時(shí)間為5:00時(shí)，室內(nèi)溫升與室外平均溫度的關(guān)系式為：

Δθr,r=0.062θo,av+3.2

式中 Δθr,r——室內(nèi)溫升，℃

θo,av——室外平均溫度，℃

供熱開(kāi)始時(shí)間為6:00時(shí)：

Δθr,r=0.043θo,av+2.4

供熱開(kāi)始時(shí)間為7:00時(shí)：

Δθr,r=0.041θo,av+1.88

② 室內(nèi)溫降

停熱時(shí)間為19:00時(shí)，室內(nèi)溫降與室外平均溫度的關(guān)系式：

Δθr,d=-0.057θo,av+0.91

式中 Δθr,d——室內(nèi)溫降，℃

停熱時(shí)間為20:00時(shí)：

Δθr,d=-0.051θo,av+0.825

停熱時(shí)間為21:00時(shí):

Δθr,d=-0.048θo,av+0.27

3.2 建筑使用期間的室溫調(diào)控策略優(yōu)化

建筑使用期間的室溫調(diào)控策略優(yōu)化主要是根據(jù)室外溫度、供暖室內(nèi)計(jì)算溫度預(yù)測(cè)熱負(fù)荷，調(diào)節(jié)熱力入口電動(dòng)調(diào)節(jié)閥相對(duì)開(kāi)度(本文將建筑使用期間的室溫調(diào)控優(yōu)化方法稱為熱力入口電動(dòng)閥室溫調(diào)控)。實(shí)現(xiàn)按需供熱，提高供熱質(zhì)量，保持良好的經(jīng)濟(jì)性。熱力入口電動(dòng)閥室溫調(diào)控流程見(jiàn)圖8。每15 min讀取1次當(dāng)前室內(nèi)溫度，是否有閥門開(kāi)度調(diào)節(jié)操作的判定時(shí)間t為15 min。建筑使用期間的室內(nèi)溫度控制范圍取18～20 ℃。

圖8 熱力入口電動(dòng)閥室溫調(diào)控流程

利用機(jī)器學(xué)習(xí)中的監(jiān)督學(xué)習(xí)方式計(jì)算熱力入口電動(dòng)調(diào)節(jié)閥相對(duì)開(kāi)度，并建立相應(yīng)的回歸模型?？紤]到供熱數(shù)據(jù)規(guī)模比較小，執(zhí)行機(jī)構(gòu)(熱力入口電動(dòng)調(diào)節(jié)閥)相對(duì)簡(jiǎn)單，復(fù)雜、多權(quán)重、多參數(shù)的算法模型并不適用。線性回歸算法作為機(jī)器學(xué)習(xí)中較為基礎(chǔ)的一種，具有簡(jiǎn)單、建模迅速、容易實(shí)現(xiàn)的特點(diǎn)，對(duì)小數(shù)據(jù)量、簡(jiǎn)單的關(guān)系很有效。而且訓(xùn)練和預(yù)測(cè)速度比較快，結(jié)果具有很好的可解釋性，因此本文采用線性回歸算法。優(yōu)化訓(xùn)練數(shù)據(jù)為340組，驗(yàn)證集數(shù)據(jù)為120組。

熱力入口電動(dòng)調(diào)節(jié)閥目標(biāo)相對(duì)開(kāi)度ka的回歸模型為：

ka=β0+β1k+a(θr,ac-θr,ai)

式中ka——熱力入口電動(dòng)調(diào)節(jié)閥目標(biāo)相對(duì)開(kāi)度

β0、β1——參數(shù)

k——當(dāng)前相對(duì)開(kāi)度

a——室內(nèi)溫度升高1 ℃閥門相對(duì)開(kāi)度減小的步長(zhǎng)

θr,ac——實(shí)際室溫，℃

θr,ai——目標(biāo)室溫，℃

采用最小二乘法的訓(xùn)練結(jié)果為：

ka=0.066+0.919k-0.14(θr,ac-θr,ai)

采用平均絕對(duì)誤差MAE、均方根誤差RMSE、相關(guān)指數(shù)R2作為驗(yàn)證集的評(píng)價(jià)指標(biāo)，計(jì)算結(jié)果分別為0.016、0.02、0.979。由計(jì)算結(jié)果可知，熱力入口電動(dòng)調(diào)節(jié)閥目標(biāo)相對(duì)開(kāi)度回歸模型效果精確度比較高。

4 結(jié)果與討論

以教學(xué)樓1～4作為分析對(duì)象，選擇2019—2020年供暖期的嚴(yán)寒期、末寒期各1個(gè)典型日作為模擬時(shí)段。嚴(yán)寒期選擇1月8日8:00至1月9日8:00(以下簡(jiǎn)稱嚴(yán)寒日)，末寒期選擇3月10日8:00至3月11日8:00(以下簡(jiǎn)稱末寒日)。設(shè)定1月8日8:00、3月10日8:00初始室溫為18 ℃，熱力入口電動(dòng)調(diào)節(jié)閥初始相對(duì)開(kāi)度為0.71。

4.1 調(diào)控策略1

調(diào)控策略1：3:00供熱，22:00停熱，在建筑使用期間，二級(jí)管網(wǎng)流量固定，熱力入口電動(dòng)調(diào)節(jié)閥保持初始相對(duì)開(kāi)度不變。

調(diào)控策略1下嚴(yán)寒日、末寒日室內(nèi)溫度、室外溫度隨時(shí)間的變化分別見(jiàn)圖9、10。由圖9、10可知，調(diào)控策略1下，嚴(yán)寒日、末寒日教學(xué)樓使用期間的室內(nèi)溫度均滿足甚至高于要求。22:00停熱后室內(nèi)溫度經(jīng)過(guò)快速下降后變得緩慢，這主要得益于建筑熱惰性。由于供熱開(kāi)始時(shí)間過(guò)早，建筑非使用期間室內(nèi)溫度偏高，1月9日8:00、3月11日8:00室內(nèi)溫度分別達(dá)到了19.5、21.8 ℃，存在著較為嚴(yán)重的過(guò)量供熱問(wèn)題。

圖9 調(diào)控策略1下嚴(yán)寒日室內(nèi)溫度、室外溫度 隨時(shí)間的變化

圖10 調(diào)控策略1下末寒日室內(nèi)溫度、室外溫度 隨時(shí)間的變化

4.2 調(diào)控策略2

調(diào)控策略2(啟停策略優(yōu)化后)：嚴(yán)寒日停熱時(shí)間為21:00，供熱時(shí)間為5:00。末寒日停熱時(shí)間為19:00，供熱時(shí)間為7:00。在建筑使用期間，二級(jí)管網(wǎng)流量固定，熱力入口電動(dòng)調(diào)節(jié)閥相對(duì)開(kāi)度保持初始相對(duì)開(kāi)度不變。

調(diào)控策略1、2下嚴(yán)寒日、末寒日室內(nèi)溫度隨時(shí)間的變化分別見(jiàn)圖11、12。由圖11、12可知，與調(diào)控策略1相比，由于調(diào)控策略2進(jìn)行了啟停策略優(yōu)化，進(jìn)一 步提前停熱、延后供熱，仍能滿足建筑使用期間的室內(nèi)溫度要求，建筑非使用時(shí)段的室內(nèi)溫度也明顯合理。建筑熱惰性確保了停熱后室內(nèi)溫度的有限下降，嚴(yán)寒日供熱開(kāi)始時(shí)，室內(nèi)溫度為16.6 ℃；末寒日供熱開(kāi)始時(shí)，室內(nèi)溫度為17.9 ℃。

圖11 調(diào)控策略1、2下嚴(yán)寒日室內(nèi)溫度隨時(shí)間的變化

圖12 調(diào)控策略1、2下末寒日室內(nèi)溫度隨時(shí)間的變化

4.3 調(diào)控策略3

調(diào)控策略3(啟停優(yōu)化策略+熱力入口電動(dòng)閥室溫調(diào)控)：嚴(yán)寒日停熱時(shí)間為21:00，供熱時(shí)間為5:00。末寒日停熱時(shí)間為19:00，供熱時(shí)間為7:00。在建筑使用期間，調(diào)節(jié)熱力入口電動(dòng)調(diào)節(jié)閥相對(duì)開(kāi)度，使室內(nèi)溫度保持在18～20 ℃。

調(diào)控策略2、3下嚴(yán)寒日、末寒日室內(nèi)溫度隨時(shí)間的變化分別見(jiàn)圖13、14。由圖13、14可知，與調(diào)控策略2相比，運(yùn)行策略3可將建筑使用期間的室內(nèi)溫度控制在要求范圍內(nèi)，進(jìn)一步降低供熱能耗。

圖13 調(diào)控策略2、3下嚴(yán)寒日室內(nèi)溫度隨時(shí)間的變化

圖14 調(diào)控策略2、3下末寒日室內(nèi)溫度隨時(shí)間的變化

4.4 優(yōu)化效果

采取啟停優(yōu)化策略+熱力入口電動(dòng)閥室溫調(diào)控，不僅使建筑使用期間的室內(nèi)溫度控制在要求范圍內(nèi)，還利用建筑熱惰性確保停熱后室內(nèi)溫度的有限下降，有利于降低供熱能耗。

5 結(jié)論

① 以穩(wěn)態(tài)傳熱模型為基礎(chǔ)，開(kāi)發(fā)的供熱系統(tǒng)(包含熱力站、二級(jí)管網(wǎng)、建筑)仿真模型的仿真結(jié)果可信。

② 調(diào)控策略1下，嚴(yán)寒日、末寒日教學(xué)樓使用期間的室內(nèi)溫度均滿足甚至高于要求。22:00停熱后室內(nèi)溫度經(jīng)過(guò)快速下降后變得緩慢，這主要得益于建筑熱惰性。由于供熱開(kāi)始時(shí)間過(guò)早，建筑非使用期間室內(nèi)溫度偏高，1月9日8:00、3月11日8:00室內(nèi)溫度分別達(dá)到了19.5、21.8 ℃，存在著較為嚴(yán)重的過(guò)量供熱問(wèn)題。

③ 與調(diào)控策略1相比，調(diào)控策略2進(jìn)行了啟停策略優(yōu)化，進(jìn)一步提前停熱、延后供熱，仍能滿足建筑使用期間的室內(nèi)溫度要求，建筑非使用時(shí)段的室內(nèi)溫度也明顯合理。建筑熱惰性確保了停熱后室內(nèi)溫度的有限下降，嚴(yán)寒日供熱開(kāi)始時(shí)，室內(nèi)溫度為16.6 ℃；末寒日供熱開(kāi)始時(shí)，室內(nèi)溫度為17.9 ℃。

④ 與調(diào)控策略2相比，調(diào)控策略3可將建筑使用期間的室內(nèi)溫度控制在要求范圍內(nèi)，進(jìn)一步降低供熱能耗。

⑤ 采取啟停優(yōu)化策略+熱力入口電動(dòng)閥室溫調(diào)控，不僅使建筑使用期間的室內(nèi)溫度控制在要求范圍內(nèi)，還利用建筑熱惰性確保停熱后室內(nèi)溫度的有限下降，有利于降低供熱能耗。