辦公建筑運行能耗的混沌時間序列復(fù)合預(yù)測

于軍琪，段佳音，趙安軍，井文強，王佳麗

(西安建筑科技大學(xué) 建筑設(shè)備科學(xué)與工程學(xué)院，西安 710055)

辦公建筑在全生命周期中普遍存在能耗大、能效低等問題，導(dǎo)致嚴(yán)重的能源浪費。其節(jié)能潛力巨大，降耗幅度可以達(dá)到30%～50%[1]。其中，辦公建筑運行能耗占比最大。因此，研究辦公建筑的運行能耗情況具有重要意義，對其進(jìn)行實時準(zhǔn)確地預(yù)測可以為優(yōu)化運行效率提供數(shù)據(jù)決策，從而實現(xiàn)節(jié)能目標(biāo)。

建筑能耗的預(yù)測方法主要分為兩大類：正演模型和數(shù)據(jù)驅(qū)動模型。筆者根據(jù)獲得的辦公建筑歷史能耗數(shù)據(jù)特點，利用數(shù)據(jù)驅(qū)動模型中的時間序列法進(jìn)行能耗短期預(yù)測。機器學(xué)習(xí)預(yù)測方法是最常用的短期能耗預(yù)測方法，已經(jīng)廣泛應(yīng)用于風(fēng)速預(yù)測[2]、電力需求預(yù)測[3]、建筑能耗預(yù)測[4]、建筑冷負(fù)荷預(yù)測[5]等能源領(lǐng)域。目前國內(nèi)外建筑能耗時間序列預(yù)測方法應(yīng)用較為廣泛的有：人工神經(jīng)網(wǎng)絡(luò)[6](ANN, artificial neural network)、差分整合移動平均自回歸[7](ARIMA, autoregressive integrated moving average)、支持向量回歸[8](SVR, support vector regression)、多元線性回歸[9](MLR, multiple linear regression)等。辦公建筑除了室外氣象參數(shù)和圍護(hù)結(jié)構(gòu)等本體能耗外，還有機電設(shè)備、節(jié)假日、人員占用率等復(fù)雜的運行能耗，造成其能耗時間序列具有非線性特征。為了深入挖掘其時間序列內(nèi)部的非線性規(guī)律，提高預(yù)測精度，研究者們將不同的算法與人工神經(jīng)網(wǎng)絡(luò)和支持向量機結(jié)合，建立了多種辦公建筑混合預(yù)測模型。文獻(xiàn)[9]根據(jù)天氣、人員等多種因素利用MLR方法對辦公建筑的逐日冷負(fù)荷進(jìn)行了回歸預(yù)測，與實際負(fù)荷相比平均絕對百分誤差小于8%。文獻(xiàn)[10]利用反饋神經(jīng)網(wǎng)絡(luò)(BPNN, back propagation neural network)，徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)(RBFNN, radial basis function neural network)，廣義回歸神經(jīng)網(wǎng)絡(luò)(GRNN, generalized regression neural network)，支持向量機(SVM, support vector machine)方法分別對辦公建筑的逐時能耗進(jìn)行分析，以中國廣州某辦公樓為實例，取得了較好的預(yù)測效果。文獻(xiàn)[11]建立了小波-支持向量機(wavelet-SVM)，小波-偏最小二乘-支持向量機(Wavelet-PLS-SVM, wavelet-partial least squares-support vector machine)模型對辦公建筑的逐時能耗進(jìn)行了分析，最終得出1 h，2 h，3 h和24 h的預(yù)測結(jié)果。在其他應(yīng)用領(lǐng)域[12]將平均速度、平均占有率和平均交通流時間序列利用貝葉斯理論集成到一個時間序列中，將混沌理論和SVR相結(jié)合用于交通流的時間序列預(yù)測，達(dá)到了較高的預(yù)測精度。上述研究都是基于多變量輸入的，相較于單變量時間序列需要考慮的因素復(fù)雜，文獻(xiàn)[13]使用歷史時刻的辦公建筑冷負(fù)荷根據(jù)貝葉斯理論映射到同一時間維度作為機器學(xué)習(xí)預(yù)測模型的單變量輸入，利用Chaos-SVR和WD-SVR(wavelet decomposition-support vector regression)分別對建筑冷負(fù)荷時間序列進(jìn)行預(yù)測，結(jié)果表明，Chaos-SVR效果更優(yōu)。因此，選擇恰當(dāng)?shù)念A(yù)測模型不僅可以避免單變量輸入的局限性，而且會取得較好的預(yù)測效果，但由于組合模型在迭代過程中存在誤差累積的問題，導(dǎo)致不能很好地滿足實際優(yōu)化運營的需求。

為了解決以上問題，筆者探究了混沌現(xiàn)象背后隱藏的決定性規(guī)律使得將混沌理論應(yīng)用于辦公建筑運行能耗預(yù)測領(lǐng)域的可能性，展開了以下研究工作：

1)分析辦公建筑運行能耗的特點及構(gòu)成，進(jìn)行數(shù)據(jù)預(yù)處理和時間序列的混沌特性判別。

2)將混沌理論(Chaos theory)和支持向量理論結(jié)合，構(gòu)建Chaos-SVR辦公建筑能耗預(yù)測模型。

3)利用Markov鏈修正Chaos-SVR組合模型中由于參數(shù)傳遞存在的累積誤差。

4)以西安某辦公建筑監(jiān)測數(shù)據(jù)為例進(jìn)行實例分析，并與其他預(yù)測方法進(jìn)行比較，驗證了模型的有效性。

文中主要研究辦公建筑運行能耗的混沌特性，其他大型公共建筑比如商業(yè)建筑、醫(yī)療建筑等能耗是否具有混沌特性，模型的可推廣性將是下一步的研究重點。

1 能耗數(shù)據(jù)特征分析

辦公建筑能耗混沌時間序列預(yù)測研究的前提是確定該時間序列具有混沌特征[14]。文中通過西安市建筑能耗監(jiān)測平臺，采集某辦公建筑2019年8月1日0點到9月30日23點各項分類能耗數(shù)據(jù)，如表1所示。

表1 西安某辦公建筑能耗量統(tǒng)計

從表1可以看出，辦公建筑運行能耗包括水耗量、電耗量、燃?xì)饬?、集中供?熱量和其他能源使用量。為了對辦公建筑運行能耗時間序列進(jìn)行混沌特性分析，由于其組成結(jié)構(gòu)復(fù)雜，需要對多項數(shù)據(jù)統(tǒng)一量綱。因此，將所有運行能耗的量綱轉(zhuǎn)化為千克標(biāo)準(zhǔn)煤(kgce)的形式，從而得到該辦公建筑運行能耗的逐時序列。然后對該逐時序列數(shù)據(jù)進(jìn)行處理，利用最大 Lyapunov指數(shù)法來判斷其是否具有混沌特性。最大Lyapunov指數(shù)法是基于相空間重構(gòu)的，因此對時間序列進(jìn)行相空間重構(gòu)是研究的前提。

1.1 相空間重構(gòu)

相空間是一種工具，用于描述由確定自由度的確定性混沌系統(tǒng)生成的單變量或多變量時間序列重構(gòu)的動態(tài)系統(tǒng)[15]。相空間重構(gòu)的主要原理是在拓?fù)涞葍r意義下從高維空間中恢復(fù)混沌吸引子的動力學(xué)特性[15]，對于總能耗的逐時時間序列x1,x2,x3,…,xn-1,xn，適當(dāng)選取嵌入維數(shù)m和延遲時間τ，得到重構(gòu)相空間

Ym(i)=[x(i),x(i+τ),x(i+2τ),…,x(i+(m-1)τ)]，

(1)

式中，i=1,2,…,n;N=n-(m-1)τ。N表示由時間序列構(gòu)成m維后，由n個相點變成了N個相點。

嵌入維數(shù)和延遲時間的合理選取決定了辦公建筑能耗時間序列相空間重構(gòu)的質(zhì)量，同時也保證了在實際預(yù)測中有足夠高的預(yù)測精度。研究表明[16]，將τ和m聯(lián)合起來的嵌入窗寬的確定也會影響重構(gòu)相空間的質(zhì)量。C-C[17]法通過序列的關(guān)聯(lián)積分來構(gòu)成統(tǒng)計量，通過統(tǒng)計量和延遲時間的關(guān)系圖同時計算出延遲時間τ和最佳窗寬[16]τw，再根據(jù)τω=(m-1)τ求出嵌入維數(shù)。該方法既能有效地減少互信息量法的計算量，又能保持時間序列的非線性特征。故本文采用C-C法求取延遲時間和嵌入維數(shù)，具體步驟如下：

(2)

式中，dij=‖Xi-Xj‖，Xi為相空間的點，θ(·)為階躍函數(shù)：

(3)

2)統(tǒng)計量代表了非線性時間序列的相關(guān)性。通過時間序列的關(guān)聯(lián)積分來構(gòu)成每個不相交的時間子序列的檢驗統(tǒng)計量：

(4)

令n→∞，則：

(5)

3)如果時間序列是獨立同分布的，對于固定的m和t，當(dāng)N→∞時，r取任意值時，S(m,r,t)恒等于零。但實際的時間序列是有限的，因此一般S(m,r,t)不等于零。那么，局部最大時間間隔可以取S(m,r,t)的零點或?qū)λ械陌霃絩相互差別最小的時間點，這意味著這些點幾乎是均勻分布的。選擇對應(yīng)值最大和最小的2個半徑r，定義關(guān)于r的最大偏差：

ΔS(m,t)=max{S(m,rj,t)}-min{S(m,rj,t)}。

(6)

局部最大時間t應(yīng)該是S(m,r,t)的零點和ΔS(m,t)的最小值。但是S(m,r,t)的零點對所有m,r應(yīng)幾乎相等；ΔS(m,t)的最小值對所有m應(yīng)幾乎相等。延遲時間τ對應(yīng)著這些局部最大時間t中的第一個。

4)通過統(tǒng)計量和延遲時間的關(guān)系圖來確定τ和嵌入窗寬τω，從而確定出嵌入維數(shù)[16]。

幾種重要漸進(jìn)分布的數(shù)學(xué)統(tǒng)計結(jié)果表明：當(dāng)2≤m≤5,σ/2≤r≤2σ,N≥500時(σ指時間序列的均方差或標(biāo)準(zhǔn)差)，漸近分布可以通過有限序列很好的近似，S(m,n,r,1)能代表序列的相關(guān)性。具體計算時，根據(jù)辦公建筑能耗時間序列特性，取m=2,3,4,5,ri=iσ/2,i=1,2,3,4,相點數(shù)n=3 000，計算下列3個檢驗統(tǒng)計量[17]：

(7)

(8)

(9)

根據(jù)上述計算結(jié)果在MATLAB中繪制圖像，得到C-C法重構(gòu)辦公建筑能耗時間序列統(tǒng)計量曲線如圖1所示。

圖1 C-C法重構(gòu)辦公建筑總能耗時間序列統(tǒng)計量曲線圖

因此，該辦公建筑總能耗時間序列重構(gòu)的相空間為Y2(i)=[x(i),x(i+2)]。

1.2 混沌特性判別

一個系統(tǒng)的Lyapunov指數(shù)譜中有一個或一個以上的正值就可以肯定其存在混沌特性[18]。在相空間重構(gòu)后，相空間中的兩條相鄰軌線將隨時間逐漸發(fā)散或聚合，而Lyapunov指數(shù)是其軌線的收斂或發(fā)散率，正的最大Lyapunov指數(shù)反映時間序列具有混純特征，負(fù)的最大Lyapunov指數(shù)表明時間序列具有隨機性或周期性。最大Lyapunov指數(shù)越大，時間序列的非線性越強，對初值越敏感[18]。

小數(shù)據(jù)量法[19]可以充分利用所有能夠利用的數(shù)據(jù)，對小數(shù)據(jù)組來說比較可靠，計算量不大，相對容易操作，計算結(jié)果精度高。因此采用小數(shù)據(jù)量方法計算最大Lyapunov指數(shù)。小數(shù)據(jù)量法求取最大Lyapunov指數(shù)的步驟如下：

1)平均周期是從相空間的意義上度量一個吸引子全局屬性的最小單元，利用快速傅里葉變換得到頻譜信息，計算出能耗的平均頻率根據(jù)其倒數(shù)估算出時間序列的平均周期τp。

經(jīng)過傅里葉變換后，平均周期可表示為

(10)

式中：fi(i=1,2,…L)是0 Hz和奈奎斯特頻率之間均勻分布的頻率點；Pi是在頻率點fi處的幅值。用C-C法計算出延遲時間τ和嵌入維數(shù)m。對于辦公建筑能耗時間序列{xi}，重構(gòu)后的相空間為式(1)。

(11)

3)根據(jù)相空間中每個點的最鄰近點的平均發(fā)散速率估計最大Lyapunov指數(shù)為

(12)

式中：Δt為樣本周期。對相空間中的每個點Ym(t)，計算出該鄰域點經(jīng)過i個離散事件步長后的距離

(13)

估計得到：

Lt(i)=Cjeλ(iΔt)，Cj=Lt(0)。

(14)

最大Lyapunov指數(shù)的幾何意義是量化初始閉軌道的指數(shù)發(fā)散和估計系統(tǒng)的總體混沌水平的量[11]，對式(14)兩邊取對數(shù)得到

lnLt(i)=lnCj+λ(iΔt)(i=1,2,…N)。

(15)

4)對每個i求出所有t的lnLt(i)平均值y(i)，即：

(16)

式中：q是非零Lt(i)的數(shù)目。最大Lyapunov指數(shù)相當(dāng)于式(15)直線的斜率，可以通過最小二乘法逼近這條直線得到。

根據(jù)相空間重構(gòu)C-C法求出的最佳延遲時間τ=2及嵌入維數(shù)m=2，利用小數(shù)據(jù)量法計算得出如圖2所示辦公建筑能耗時間序列的最小二乘法擬合直線。

圖2 最小二乘法擬合直線

根據(jù)直線斜率可以得到最大Lyapunov指數(shù)值為0.007 7，此為正值，由此可得出辦公建筑總能耗量時間序列具有混沌特性。

2 Chaos-SVR組合模型預(yù)測

辦公建筑的Chaos-SVR的預(yù)測方法主要運用混沌理論，構(gòu)造非線性映射利用相空間重構(gòu)技術(shù)將原始時間序列數(shù)據(jù)所隱含的信息進(jìn)行還原重構(gòu)，并將其作為SVR模型的輸入，進(jìn)行非線性訓(xùn)練和預(yù)測?；谵k公建筑能耗的Chaos-SVR時間序列模型建立及預(yù)測過程如下：

1)構(gòu)建預(yù)測模型結(jié)構(gòu)，主要包括3部分：輸入向量、支持向量的非線性變化和輸出。西安某辦公建筑運行能耗的時間序列進(jìn)行重構(gòu)后的相空間為Y2(i)=[x(i),x(i+2)]，相點數(shù)為N=n-2，每個相點包含了每個序列的主要特征，能夠逼近建筑能耗的真實狀態(tài)。向量xi=(xi,l)是模型的輸入向量，β=(β1,β2,…βl)是輸出的權(quán)重向量，yi=f(xi)是模型的輸出。Chaos-SVR模型結(jié)構(gòu)如圖3所示。

圖3 Chaos-SVR模型結(jié)構(gòu)圖

3)進(jìn)行網(wǎng)絡(luò)循環(huán)學(xué)習(xí)，修正SVR的輸出及權(quán)值，直到誤差控制在允許范圍或迭代次數(shù)達(dá)到上限，訓(xùn)練結(jié)束。

4)取樣本前n個數(shù)據(jù)之后連續(xù)的Q個數(shù)據(jù)點作為預(yù)測數(shù)據(jù)。對模型進(jìn)行預(yù)測，網(wǎng)絡(luò)的輸出值即為預(yù)測值，將預(yù)測值進(jìn)行相應(yīng)的反歸一化處理，得到預(yù)測的實際值。

為了驗證模型的有效性，選取2019年9月4日～9月24日共計20天480組數(shù)據(jù)作為訓(xùn)練數(shù)據(jù)，以9月25日和26日2天共計48組數(shù)據(jù)作為預(yù)測數(shù)據(jù)來進(jìn)行模型驗證。并與目前在時間序列預(yù)測領(lǐng)域應(yīng)用廣泛的2種機器學(xué)習(xí)方法——非線性自回歸神經(jīng)網(wǎng)絡(luò)(nonlinear auto regressive neural network)和支持向量回歸(SVR)方法的預(yù)測值和實際值進(jìn)行比較，利用均方根誤差RMSE(root mean squared error)平均絕對百分比誤差MAPE(mean absolute percent error)作為評價指標(biāo)。預(yù)測對比結(jié)果如圖4所示。

圖4 預(yù)測結(jié)果對比圖

預(yù)測結(jié)果誤差對比如圖5所示。

圖5 預(yù)測結(jié)果誤差對比圖

預(yù)測結(jié)果誤差如表2所示。

表2 預(yù)測模型誤差值

從圖4中可以看出，Narnet預(yù)測曲線波動最大，SVR預(yù)測曲線較為平緩，Chaos-SVR預(yù)測曲線距離真實值曲線最接近，部分值幾乎完全擬合。由圖5得，Narnet預(yù)測方法誤差最大，SVR預(yù)測方法誤差較大，Chaos-SVR預(yù)測方法誤差在0值附近小范圍浮動。由表2可以看出，Chaos-SVR預(yù)測模型的RMSE為6.326 1，MAPE為0.287 5，相較于其他2種方法均最小。實驗是在同樣的環(huán)境和同一組數(shù)據(jù)的基礎(chǔ)上進(jìn)行的，Chaos-SVR預(yù)測方法相較于Narnet和SVR方法，預(yù)測效果最好，預(yù)測精度明顯提高。但Chaos-SVR組合模型在預(yù)測的過程中，還有部分值存在較大的相對誤差，為了更加精確地進(jìn)行辦公建筑運行能耗預(yù)測，需要對這部分的誤差進(jìn)行修正。導(dǎo)致這種誤差存在的原因可能是由于組合模型在迭代過程中，由于參數(shù)傳遞而出現(xiàn)的過程誤差積累，這樣的誤差積累會造成預(yù)測效果存在一定的缺陷。為了進(jìn)一步提高預(yù)測精度，文中利用Markov鏈進(jìn)行修正。

3 Markov鏈修正誤差

Markov鏈?zhǔn)且环N具有無后效性的隨機過程[20]，Markov鏈適用于對存在波動性的預(yù)測問題進(jìn)行修正描述[21]。它可以根據(jù)某一已知條件的時刻條件下推求出下一時刻的概率分布，得到下一時刻的狀態(tài)而與其他時刻無關(guān)。因此，在傳統(tǒng)Chaos-SVR預(yù)測模型基礎(chǔ)之上利用Markov鏈對累積誤差進(jìn)行修正，從而實現(xiàn)辦公建筑能耗時間序列的精確預(yù)測。

1)誤差狀態(tài)區(qū)間[22]劃分。

(17)

(18)

2)狀態(tài)轉(zhuǎn)移概率矩陣。

(19)

(20)

使用C-K方程(切普曼-柯爾莫哥洛夫方程)計算k步的轉(zhuǎn)移概率矩陣為

P(k)=(P(1))k。

(21)

3)預(yù)測模型。

設(shè)初始向量為P(0)，經(jīng)k步轉(zhuǎn)移后的狀態(tài)向量為

P(k)=P(0)×P(k)=P(0)×(P(1))k，

(22)

據(jù)此可以判斷k步所處的狀態(tài)空間。

4)根據(jù)步驟2)所得概率轉(zhuǎn)移矩陣和步驟3)所得狀態(tài)向量，得出誤差修正后的值。

根據(jù)劃分的狀態(tài)區(qū)間對Chaos-SVR模型預(yù)測值的誤差值進(jìn)行狀態(tài)劃分，劃分結(jié)果如表3所示。

表3 Chaos-SVR模型預(yù)測誤差分類結(jié)果

根據(jù)上述誤差狀態(tài)區(qū)間的確定可以得到Chaos-SVR預(yù)測結(jié)果的Markov鏈的狀態(tài)轉(zhuǎn)移情況，進(jìn)一步計算出狀態(tài)轉(zhuǎn)移概率矩陣P為

(23)

根據(jù)Markov鏈預(yù)測模型得到辦公建筑9月25日的能耗預(yù)測狀態(tài)向量，針對Chao-SVR預(yù)測模型存在的累積誤差進(jìn)行修正，得到該辦公建筑能耗的最終預(yù)測值。為了驗證基于Markov鏈修正后的Chaos-SVR預(yù)測效果，對修正前后的預(yù)測結(jié)果進(jìn)行比較。對比結(jié)果如圖6所示。

圖6 Markov修正前后預(yù)測結(jié)果對比圖

Markov鏈修正前后預(yù)測結(jié)果誤差對比如圖7所示。

圖7 Markov鏈修正前后預(yù)測結(jié)果誤差對比圖

Markov鏈修正前后的誤差計算結(jié)果如表4所示。

表4 Markov鏈修正前后誤差值

由圖6得，經(jīng)Markov修正后的模型預(yù)測曲線和真實值曲線擬合度很高，直觀地可以看到主要對2個時間點的運行能耗值進(jìn)行了修正。從圖7預(yù)測結(jié)果誤差對比中可以得到證實，經(jīng)過Markov鏈修正之后得誤差非常小，不僅對部分累積誤差進(jìn)行修正接近于0值，其他誤差也進(jìn)行了相應(yīng)修正，總體預(yù)測精度顯著提高。由表4可以看到，經(jīng)Markov修正之后的RMSE和MAPE均降低，這說明經(jīng)過Markov鏈對組合模型中參數(shù)傳遞存在的累積誤差進(jìn)行了有效改進(jìn)，它的修正過程適用于辦公建筑運行能耗這一非線性系統(tǒng)。Markov鏈的無效性特征在組合模型的誤差修正中得到了充分體現(xiàn)，只對最后的預(yù)測結(jié)果做處理而不影響傳統(tǒng)預(yù)測模型的運算過程，更符合實際辦公建筑運行能耗的變化規(guī)律。因此，與Chaos-SVR時間序列預(yù)測模型相比，基于Markov鏈修正的Chaos-SVR時間序列預(yù)測模型更優(yōu)良。

4 結(jié) 論

針對辦公建筑運行能耗預(yù)測問題，本文根據(jù)混沌理論構(gòu)造非線性映射對辦公建筑運行能耗時間序列進(jìn)行相空間重構(gòu)，提出Chaos-SVR的組合預(yù)測模型，并利用Markov鏈對組合模型中累積誤差進(jìn)行修正的方法。并將預(yù)測模型應(yīng)用于西安某實際的辦公建筑，得出以下結(jié)論：

1)用C-C法求解辦公建筑總能耗時間序列嵌入維數(shù)和延遲時間，對其進(jìn)行了相空間重構(gòu)。同時根據(jù)小數(shù)據(jù)量法求解最大Lyapunov指數(shù)值為0.007 7，由此判定辦公建筑能耗時間序列存在混沌特性。

2)將常用的預(yù)測方法比如Narnet,SVR與本文所建立方法Chaos-SVR預(yù)測結(jié)果進(jìn)行比較，本文所建立方法預(yù)測誤差最小，RMSE為6.326 1，MAPE為0.287 5，利用Markov鏈對組合模型累積誤差修正之后的RMSE和MAPE分別降為0.565 5和0.044 1。結(jié)果表明，Markov修正之后的預(yù)測精度顯著提高，為辦公建筑的優(yōu)化運營與節(jié)能管理提供充分的決策依據(jù)，符合工程實際應(yīng)用。