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

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

1 概述

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

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

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

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

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

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

① 熱力站模型

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

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

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

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

t——時間，s

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

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

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

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

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

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

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

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

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

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

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

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

由熱平衡關(guān)系，得到二級管網(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——室外溫度，℃

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

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

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

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

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

圖4 建筑幾何模型

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

2.3 模型驗證

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

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

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

表3 各建筑模擬日供暖能耗的仿真值與實測值的相對誤差 %

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

3.1 啟停優(yōu)化

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

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

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

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

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

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

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

① 室內(nèi)溫升

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

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

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

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

供熱開始時間為6:00時：

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

供熱開始時間為7:00時：

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

② 室內(nèi)溫降

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

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

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

停熱時間為20:00時：

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

停熱時間為21:00時:

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

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

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

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

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

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

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

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

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

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

a——室內(nèi)溫度升高1 ℃閥門相對開度減小的步長

θr,ac——實際室溫，℃

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

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

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

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

4 結(jié)果與討論

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

4.1 調(diào)控策略1

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

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

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

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

4.2 調(diào)控策略2

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

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

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

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

4.3 調(diào)控策略3

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

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

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

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

4.4 優(yōu)化效果

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

5 結(jié)論

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

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

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

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

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