基于深度學(xué)習(xí)的房間冷負荷預(yù)測模型

林 越 劉廷章

(1.上海大學(xué)機電工程與自動化學(xué)院, 上海 200444;2.海南熱帶海洋學(xué)院理學(xué)院, 海南三亞 572022)

民用建筑能耗約占社會總能耗的20%, 有巨大的節(jié)能潛力[1].針對建筑冷負荷節(jié)能, 不同領(lǐng)域的學(xué)者從各自專業(yè)角度提出了不同的建筑節(jié)能方案.在空調(diào)系統(tǒng)運行過程節(jié)能與控制研究中[2-3], 除了在設(shè)定溫度不變情況下對供能側(cè)的調(diào)度管理外, 還有學(xué)者提出適用于空調(diào)房間的溫濕度比例-積分-微分(propotional-integral-differential, PID) 解耦控制系統(tǒng).基于PCTT(principal component threshold technique) 矩陣算法的溫濕度PID 解耦控制系統(tǒng)不但消除了溫、濕度調(diào)節(jié)通道之間存在的相互影響, 且其控制性能指標更好.到目前為止, 很多控制策略基本通過控制房間溫度變化來實現(xiàn)建筑節(jié)能.但這些研究尚處于起始階段, 缺乏溫度波動工況下的普適、定量的負荷預(yù)測模型.

估計空調(diào)的冷負荷, 就是在房間的熱量動態(tài)變化的狀態(tài)下, 尋找使得房間溫度維持在設(shè)定值的能量需求.現(xiàn)有負荷計算方法主要有2 種: 基于機理分析和基于數(shù)據(jù)學(xué)習(xí).基于機理建模的方法需要精確的參數(shù)測定且計算過程復(fù)雜, 適用于離線計算, 通常用于建筑設(shè)計階段以及空調(diào)安裝過程.基于數(shù)據(jù)學(xué)習(xí)的方法則在冷負荷計算過程中引入先進的控制算法, 例如自回歸(auto regression, AR)[4]、數(shù)據(jù)挖掘(data mining)[5]、灰色系統(tǒng)(grey system)[6]和人工神經(jīng)網(wǎng)絡(luò)(artificial neural network, ANN)[7]等.基于數(shù)據(jù)學(xué)習(xí)的方法往往并不區(qū)分冷負荷與供冷量, 通常用于供能側(cè)調(diào)控[8-9].這2 種方法均不適用于調(diào)溫模式下的房間冷負荷預(yù)測, 無法滿足基于空調(diào)運行過程中的負荷側(cè)調(diào)控需求.基于此, 本工作提出了基于深度循環(huán)神經(jīng)網(wǎng)絡(luò)的針對調(diào)溫模式下房間冷負荷預(yù)測方法.

神經(jīng)網(wǎng)絡(luò)憑借其超強的非線性擬合與泛化能力得到了廣泛的關(guān)注與應(yīng)用.深度學(xué)習(xí)[10]自2006 年作為一種快速訓(xùn)練方法被提出以來, 廣泛應(yīng)用到語音、圖像、視頻等領(lǐng)域[11].通過深度架構(gòu)的開發(fā), 深度學(xué)習(xí)方法能夠從數(shù)據(jù)中發(fā)現(xiàn)更深層次的隱藏結(jié)構(gòu)與特征的抽象關(guān)系.深度學(xué)習(xí)基本架構(gòu)包括深度信念網(wǎng)絡(luò)(deep belief network, DBN)、堆疊自動編碼器(stacked auto encoder, SAE)、卷積神經(jīng)網(wǎng)絡(luò)(convolutional neural network, CNN)、循環(huán)神經(jīng)網(wǎng)絡(luò)(recurrent neural networks, RNN).RNN 考慮了樣本之間的關(guān)聯(lián)關(guān)系, 這種關(guān)聯(lián)關(guān)系以神經(jīng)網(wǎng)絡(luò)之間的鏈接來體現(xiàn).

針對房間溫度動態(tài)變化工況下的房間冷負荷預(yù)測, 本工作首先根據(jù)房間能量平衡方程, 通過分析冷負荷、供冷量、蓄熱量之間的關(guān)系, 提出房間冷負荷預(yù)測模型; 然后, 利用頻域分解法實現(xiàn)蓄熱放熱系數(shù)的辨識以及冷負荷軟測量, 結(jié)合深度循環(huán)神經(jīng)網(wǎng)絡(luò)的數(shù)據(jù)學(xué)習(xí)能力進行室溫恒定條件下的冷負荷預(yù)測; 最后, 將室溫恒定條件下的負荷預(yù)測值與溫度波動所帶來的蓄熱放熱相加, 即可得到調(diào)溫模式下的房間冷負荷預(yù)測值.本方法的有效性通過基于建筑能耗仿真軟件的仿真實驗和基于某實驗樓的實測實驗進行驗證.

1 模 型

冷負荷是指房間溫度維持在設(shè)定溫度時所需要的能量排出率.在室內(nèi)溫度非恒定的工況下, 房間熱力系統(tǒng)存在能量平衡方程[12]:

式中:cl(n) 為以恒定室溫作為參考溫度的當前時刻逐時冷負荷;he(n) 為當前時刻空調(diào)系統(tǒng)逐時供冷量, 供冷量的一部分用于抵消室內(nèi)負荷; 差值部分hs(n) 為房間逐時蓄熱量.

空調(diào)供冷量可以通過計算空調(diào)出風(fēng)口與進風(fēng)口之間焓差的方式得到.

式中:G(n) 為當前時刻空氣逐時質(zhì)量流量;hout(n) 為出風(fēng)口空氣的焓;hin(n) 為進風(fēng)口空氣的焓;cg為干空氣定壓比熱;cq為水蒸氣定壓比熱;dout為出風(fēng)口單位質(zhì)量干空氣所含水蒸氣質(zhì)量;din(n) 為進風(fēng)口單位質(zhì)量干空氣所含水蒸氣質(zhì)量;tout(n) 為出風(fēng)口空氣溫度;tin(n) 進風(fēng)口空氣溫度.

蓄熱量定義為房間所需冷負荷與空調(diào)系統(tǒng)實際供冷量之間的能量差, 存儲于建筑圍護結(jié)構(gòu)、室內(nèi)家具等蓄熱體.蓄熱量與房間實際溫度相對于設(shè)定溫度偏離量的關(guān)系為

式中: ?t(n-i) 為i小時前室內(nèi)實際溫度相對于參考溫度的偏離量,i=0 為當前時刻;ki為蓄熱放熱系數(shù), 反應(yīng)了蓄熱與該時刻溫度差的數(shù)值關(guān)系.

以某一恒定室溫作為參考, 溫度動態(tài)調(diào)控模式下房間冷負荷計算包含2 部分: ①室溫恒定這一理想假設(shè)下的房間冷負荷; ②房間的蓄熱放熱[13].

式中:Load(n+1) 為下一時刻房間冷負荷;he(n+1) 為下一時刻讓室內(nèi)溫度維持在某一設(shè)定溫度應(yīng)該有的空調(diào)供冷量;cl(n+1) 為下一時刻房間溫度恒定在參考溫度下的冷負荷;hs(n+1) 為下一時刻房間的蓄熱放熱.

綜上可知, 如果能實現(xiàn)恒定溫度下冷負荷與房間蓄熱量的計算, 就可實現(xiàn)冷負荷的預(yù)測.

2 方 法

根據(jù)式(4), 要實現(xiàn)溫度波動工況下的房間冷負荷預(yù)測, 需要將其分解為2 部分: 跟溫度波動相關(guān)的蓄熱量和房間溫度恒定工況下的冷負荷預(yù)測值.

2.1 頻域分解

室內(nèi)溫度恒定是一個理想條件, 與實際情況不相符.文獻[14]利用頻域分解的方法, 在室溫波動的實際工況下實現(xiàn)蓄熱量的估計, 并用于實現(xiàn)室溫恒定條件房間冷負荷的軟測量.仿真實驗和實測實驗均顯示該方法有較好的計算效果.

恒定室溫下的冷負荷在頻域分解之后, 能量集中在以24 h 為周期的頻率點及其諧波上.經(jīng)過傅里葉變換, 能量平衡方程重寫為

式中:cl(ω) 為室溫恒定條件下冷負荷的頻域變換;he(ω) 為空調(diào)供冷量頻域變換;dt(ω) 為逐時溫度差的頻域變換.

根據(jù)室溫恒定條件下冷負荷的頻域特性, 可得

式中:K為蓄熱計算截斷階數(shù);ω1、ω2、ω3分別為以24 h 為周期的1、2、3 次諧波頻率.

綜上, 可求得式(3) 中的蓄熱放熱系數(shù), 并計算蓄熱量.

2.2 深度循環(huán)神經(jīng)網(wǎng)絡(luò)

在得到蓄熱放熱系數(shù)后, 可以計算蓄熱量, 并根據(jù)式(1) 得到溫度恒定條件下的冷負荷的軟測量值.當前時刻負荷計算是為了之后的供能提供參考.神經(jīng)網(wǎng)絡(luò)能實現(xiàn)非常復(fù)雜的非線性映射而被廣泛應(yīng)用.根據(jù)冷負荷的軟測量值, 采用神經(jīng)網(wǎng)絡(luò)可以預(yù)測

式中:NN(·) 為神經(jīng)網(wǎng)絡(luò);CL(n) 為截至當前時刻的冷負荷時間序列.

循環(huán)神經(jīng)網(wǎng)絡(luò)是適合于序列數(shù)據(jù)建模的常用模型.一般來講, RNN 包含有神經(jīng)元輸出再到輸入的前饋鏈接, 潛在的缺點是其針對缺少當前時間步長輸入的分層處理.深度循環(huán)神經(jīng)網(wǎng)絡(luò)(deep recurrent neural network, DRNN) 正是為克服RNN 的缺點而開發(fā), 能處理多時間步分層信息, 并成功應(yīng)用于語音分離.

圖1 為DRNN 信息流向示意圖, 其中I代表輸入,O為輸出,L(1)、L(2)、L(3) 為神經(jīng)網(wǎng)絡(luò)層數(shù),I(i) 為第i時刻輸入,O(i) 為第i時刻輸出.圖1(b) 為圖1(a) 沿時間方向展開, 箭頭代表信息傳遞方向.因為同層之間的鏈接不需要都有, 所以其中一些箭頭可以省略.神經(jīng)網(wǎng)絡(luò)層數(shù)、各層神經(jīng)元個數(shù)可以根據(jù)實際需要調(diào)節(jié), 神經(jīng)元與其前面時刻的輸出之間的鏈接也可以根據(jù)情況選擇連接或者不連接.

圖1 DRNN 信息流向示意圖Fig.1 Information flow in deep recurrent neural networks

圖2 高斯-牛頓-LM 自適應(yīng)切換Fig.2 Adaptive switching between GN and LM

權(quán)重更新過程為

式中:wd+1為更新后的權(quán)重;wd為更新前權(quán)重; ?w為權(quán)重修正值.

DRNN 權(quán)重更新方法為沿時間反向傳播(back propagation through time, BPTT).在BPTT 過程中, 梯度下降(gradient descent, GD) 法收斂速度較慢, 不適用于處理大數(shù)據(jù)的深度學(xué)習(xí), 而高斯-牛頓(Gauss-Newton, GN) 法收斂速度快, 但在當前值遠離最優(yōu)值時不適用.萊文伯格-馬夸特(Levenberg-Marquardt, LM) 算法是一個求解非線性方程的重要方法, 相當于梯度下降法與高斯-牛頓法的折中方案, 但是如果參數(shù)選取不當, 該方法可能會導(dǎo)致過早收斂于局部最優(yōu).基于此, 本工作選取一種高斯-牛頓-LM 法自適應(yīng)切換的參數(shù)調(diào)整算法.

式中: Hessi 為所構(gòu)造近似海森矩陣的特征值最小值; Flag 是一個標志位; GN 是高斯-牛頓法程序; LM 是LM 法程序.在每次權(quán)重成功更新之后, 置Flag=0; 如果每次權(quán)重更新之后誤差沒有沿著下降的方向取值, 則該次權(quán)重更新不被接受, 置Flag=1, 轉(zhuǎn)向LM 法程序.

綜上可知, 通過式(7) 可以實現(xiàn)恒溫模式下空調(diào)房間冷負荷預(yù)測.

2.3 冷負荷預(yù)測

根據(jù)式(4), 綜合溫度波動情況下的蓄熱量計算和溫度恒定條件下的冷負荷預(yù)測值, 可以得到調(diào)溫模式下房間冷負荷預(yù)測.

圖3 給出了本方法的房間冷負荷預(yù)測的實現(xiàn)過程, 其中輸入量為逐時空調(diào)供冷量序列、逐時室內(nèi)溫度序列, 預(yù)測時刻的房間溫度為期望達到的房間溫度設(shè)定.

圖3 房間冷負荷預(yù)測Fig.3 Prediction for room cooling load

3 實 驗

根據(jù)式(4) 和圖3, 采用本方法實現(xiàn)調(diào)溫模式下房間冷負荷預(yù)測, 需要2 組逐時數(shù)據(jù), 即供冷量序列和溫度差序列, 且這2 組序列均可測或有公認方法可計算.進一步通過仿真實驗與實測實驗驗證所提方法的有效性, 其中仿真實驗基于DeST 能耗仿真軟件, 實測實驗選在上海市寶山區(qū)某實驗樓中的一個房間內(nèi)進行.

3.1 仿真實驗

DeST 需要依賴外部畫圖軟件給出建筑結(jié)構(gòu)模型.軟件所需相關(guān)選項設(shè)置好之后, 輸入房間逐時溫度序列, 即可得到空調(diào)系統(tǒng)逐時供冷量.本實驗輸入逐時溫度為(25±5)?C, 該序列采用MATLAB 隨機數(shù)生成函數(shù)生成, 然后手動輸入DeST.

訓(xùn)練數(shù)據(jù)選擇全年第4 000～5 000 h (典型氣象年下全年共8 760 h).驗證數(shù)據(jù)為第5 151～5 200 h.預(yù)測結(jié)果如圖4 所示, 其中參考值為DeST 的負荷計算值.

圖4 仿真實驗Fig.4 Simulation experiment

DRNN 有4 層: 第一層為輸入層, 2 個神經(jīng)元; 第二層為隱含層, 5 個神經(jīng)元, 連接至前3時刻的輸出; 第3 層為隱含層, 3 個神經(jīng)元, 連接至前2 時刻輸出; 第四層為輸出層.

3.2 實測實驗

選取上海市寶山區(qū)某實驗樓的一個房間進行實測實驗和數(shù)據(jù)采集.溫度采集用Pt100溫度傳感器, 校準后連接至ICP DASR PISO-813 數(shù)據(jù)采集卡, 上位機選用ADVANTECHR 610H 工控機, 組態(tài)軟件為力控ForceControlR6.1.濕度采集使用福祿克FLUKER 54II, 因為影響較小且室內(nèi)溫度波動不大, 所以采樣周期較長, 但保持在每天采樣5 次以上, 用插值法估計逐時濕度.

選取2015 年1 月5 日—1 月13 日的數(shù)據(jù)作為訓(xùn)練數(shù)據(jù), 選取1 月14 日—1 月15 日的數(shù)據(jù)作為驗證數(shù)據(jù).實驗結(jié)果如圖5 所示, 其中參考值為下一時刻的空調(diào)供冷量.

實測得到的實驗數(shù)據(jù)相對較少, 為避免過擬合, 選擇的DRNN 結(jié)構(gòu)層數(shù)和神經(jīng)元個數(shù)都相對較少(相對于仿真實驗).DRNN 的第一層為輸入層, 2 個神經(jīng)元; 第二層為隱含層, 5 個神經(jīng)元, 連接至前2 個時刻的該層輸入; 第三層為輸出層.

3.3 不同方法的比較

為了驗證本方法的有效性,選取了常見的BP 神經(jīng)網(wǎng)絡(luò)(back-propagation neural network,BPNN) 模型、自回歸(auto-regressive, AR) 模型和自回歸滑動平均(auto-regressive moving average, ARMA) 模型作為對比實驗, 同時選取均方誤差(mean-square error, MSE)、平均絕對百分誤差(mean absolute percentage error, MAPE) 和標準差(standard deviation, SD) 作為精度評價指標.

式中:N為測試樣本容量;yi和xi分別為測試樣本的第i個預(yù)測值和真實值.

4 種模型對應(yīng)的仿真實驗和實測實驗結(jié)果如表1 所示.

表1 4 種模型的預(yù)測性能Table 1 The forecasting performance of 4 kinds of models

由表1 可知:

(1) 無論在仿真實驗階段或者在實測實驗階段, 本方法的各種精度評價指標均優(yōu)于其他3種傳統(tǒng)機器學(xué)習(xí)方法.因此, 本方法取得了優(yōu)于傳統(tǒng)機器學(xué)習(xí)的良好表現(xiàn), 特別是平均絕對百分誤差僅為1.44%, 這個結(jié)果能滿足工程要求.

(2) 仿真實驗階段的精度整體優(yōu)于實測實驗階段.這是因為一方面建筑能耗仿真軟件具有線性假設(shè), 另一方面實測實驗受實驗條件限制, 采集的數(shù)據(jù)相對較少.這導(dǎo)致該階段的測試精度、模型誤差、擬合誤差等比仿真實驗階段總體要差, 但與其他機器學(xué)習(xí)算法比較, 本方法依然達到了較好的精度.

(3) 本方法在冷負荷預(yù)測所需要的計算時間反映了其復(fù)雜度, 側(cè)面反應(yīng)了其經(jīng)濟性.當學(xué)習(xí)率與判斷誤差在小范圍內(nèi)搭配時, 達到所規(guī)定的判斷誤差時所需的時間最短, 這可從計算次數(shù)中得到反映.

4 結(jié) 論

暖通空調(diào)系統(tǒng)在調(diào)溫模式下房間負荷預(yù)測傳統(tǒng)方法缺少基于模型的定量研究.針對房間冷負荷建模與預(yù)測的需求, 本工作首先在房間溫度動態(tài)變化工況下分析冷負荷、供冷量、蓄熱量之間的關(guān)系, 并在此基礎(chǔ)上搭建一個普適模型, 再利用頻域分解法辨識蓄熱放熱系數(shù), 結(jié)合深度循環(huán)神經(jīng)網(wǎng)絡(luò)實現(xiàn)房間溫度恒定工況下的冷負荷預(yù)測.綜合溫度波動工況下的蓄熱量和溫度恒定情況下的負荷預(yù)測值, 可得到調(diào)溫模式下的房間冷負荷預(yù)測值.仿真和實測實驗表明, 本方法能快速有效地實現(xiàn)房間逐時冷負荷預(yù)測.本方法快速精確地實現(xiàn)了調(diào)溫模式下房間逐時負荷需求預(yù)測, 可用于實現(xiàn)建筑被動熱儲能以及負荷側(cè)調(diào)控的定量計算, 也為整個電網(wǎng)需求側(cè)直接負荷控制提供了可借鑒的思路.