基于rbf神經(jīng)網(wǎng)絡(luò)的建筑逐時(shí)空調(diào)負(fù)荷智能預(yù)測(cè)系統(tǒng)

基于rbf神經(jīng)網(wǎng)絡(luò)的建筑逐時(shí)空調(diào)負(fù)荷智能預(yù)測(cè)系統(tǒng)

正確預(yù)測(cè)建筑物空調(diào)負(fù)荷是實(shí)現(xiàn)空調(diào)系統(tǒng)優(yōu)化運(yùn)行的先決條件。目前，預(yù)測(cè)空調(diào)負(fù)荷預(yù)測(cè)方法主要包括線(xiàn)性回歸法、指數(shù)平滑法、神經(jīng)網(wǎng)絡(luò)法和層析法。其中，人工神經(jīng)網(wǎng)絡(luò)具有強(qiáng)大的非線(xiàn)性映射能力，因此在預(yù)測(cè)空調(diào)負(fù)荷方面得到了廣泛應(yīng)用。然而，目前采用的人工神經(jīng)網(wǎng)絡(luò)算法基本上是bp算法或其改進(jìn)版本。當(dāng)bp網(wǎng)絡(luò)用于函數(shù)接近時(shí)，加權(quán)值的調(diào)整采用負(fù)梯度下降法。這種調(diào)整權(quán)值的方法存在收斂慢、局部小的缺點(diǎn)。徑向函數(shù)（rbf）的神經(jīng)網(wǎng)絡(luò)是一個(gè)比p網(wǎng)絡(luò)更快的神經(jīng)網(wǎng)絡(luò)。因此，作者使用rbf神經(jīng)網(wǎng)絡(luò)建立了智能預(yù)測(cè)模型，并對(duì)該模型進(jìn)行了應(yīng)用于辦公樓和圖書(shū)館的空調(diào)負(fù)荷進(jìn)行預(yù)測(cè)，并將其與bp-s算法的預(yù)測(cè)結(jié)果進(jìn)行了比較。我們將驗(yàn)證rbf神經(jīng)網(wǎng)絡(luò)對(duì)建筑物空調(diào)負(fù)荷預(yù)測(cè)的可行性。1rbf神經(jīng)網(wǎng)絡(luò)用于預(yù)測(cè)空調(diào)負(fù)荷1.1rbf神經(jīng)網(wǎng)絡(luò)RBF神經(jīng)網(wǎng)絡(luò)是一種具有全局收斂特性的線(xiàn)性學(xué)習(xí)算法的前饋網(wǎng)絡(luò).它已被廣泛應(yīng)用于數(shù)據(jù)的分類(lèi)和時(shí)間序列的預(yù)測(cè)等方面.RBF神經(jīng)網(wǎng)絡(luò)由3層組成,其結(jié)構(gòu)如圖1所示.輸入層節(jié)點(diǎn)只傳遞輸入信號(hào)到隱層,隱層節(jié)點(diǎn)一般由輻射狀的作用函數(shù)構(gòu)成,通常采用高斯函數(shù).輸入層實(shí)現(xiàn)從xi→Rj(xi)的非線(xiàn)性映射;而輸出層節(jié)點(diǎn)通常選取線(xiàn)性函數(shù),從而使輸出層實(shí)現(xiàn)從Rj(xi)→yk的線(xiàn)性映射.1.2rbf基于神經(jīng)網(wǎng)絡(luò)的空調(diào)負(fù)荷預(yù)測(cè)模型1.2.1空調(diào)負(fù)荷預(yù)測(cè)模型建筑物空調(diào)負(fù)荷受到諸多參數(shù)的影響,關(guān)系較為復(fù)雜.一般而言,建筑物空調(diào)負(fù)荷是由建筑物外擾和內(nèi)擾引起的,室外氣象參數(shù)(溫度、濕度、太陽(yáng)輻射和風(fēng)速等)和室內(nèi)人員變動(dòng)與設(shè)備啟停是引起建筑物負(fù)荷波動(dòng)的主要原因.對(duì)于某些人員變動(dòng)和設(shè)備啟停較為規(guī)律的建筑,如辦公樓、圖書(shū)館等,內(nèi)擾變化規(guī)律相對(duì)穩(wěn)定,對(duì)空調(diào)負(fù)荷的影響不大.而在室外氣象參數(shù)中,溫度、濕度及太陽(yáng)輻射對(duì)空調(diào)負(fù)荷的影響最大.因此,文中將溫度、濕度和太陽(yáng)輻射作為空調(diào)負(fù)荷預(yù)測(cè)模型的輸入變量.此外,考慮到溫度、太陽(yáng)輻射強(qiáng)度對(duì)空調(diào)動(dòng)態(tài)負(fù)荷的延遲作用,將它們?cè)谇耙粋€(gè)或兩個(gè)時(shí)刻的歷史值也納入輸入層.令當(dāng)前時(shí)刻為ue84aSymboltAp,則RBF模型的輸入量可選擇以下6個(gè)參數(shù):ue84aSymboltAp、ue84aSymboltAp-1、ue84aSymboltAp-2時(shí)刻室外空氣干球溫度tue84aSymboltAp、tue84aSymboltAp-1、tue84aSymboltAp-2,ue84aSymboltAp時(shí)刻室外空氣含濕量due84aSymboltAp,ue84aSymboltAp、ue84aSymboltAp-1時(shí)刻太陽(yáng)輻射照度lue84aSymboltAp、lue84aSymboltAp-1.建筑物的逐時(shí)空調(diào)負(fù)荷P作為模型的輸出.1.2.2室外空氣干球溫度為了提高計(jì)算效率,防止計(jì)算過(guò)程中出現(xiàn)個(gè)別數(shù)據(jù)溢出的不良現(xiàn)象,需要將輸入和輸出參數(shù)預(yù)先進(jìn)行歸一化處理,即把數(shù)據(jù)處理成0～1之間的小數(shù)值.文中的參數(shù)采取如下歸一化處理方法:溫度t′=t-tmintmax-tmin(1)t′=t?tmintmax?tmin(1)含濕量d′=d-dmindmax-dmin(2)d′=d?dmindmax?dmin(2)太陽(yáng)輻射照度l′=l-lminlmax-lmin(3)l′=l?lminlmax?lmin(3)空調(diào)負(fù)荷Ρ′=Ρ-ΡminΡmax-Ρmin(4)P′=P?PminPmax?Pmin(4)式(1)～(4)中:t、d、l、P分別為樣本中室外空氣干球溫度、含濕量、太陽(yáng)輻射照度及空調(diào)負(fù)荷;下標(biāo)max、min分別表示最大、最小值;t′、d′、l′、P′分別為t、d、l、P經(jīng)過(guò)歸一化處理后的值.1.2.3結(jié)構(gòu)神經(jīng)網(wǎng)絡(luò)建模將經(jīng)過(guò)歸一化處理后的ue84aSymboltAp、ue84aSymboltAp-1和ue84aSymboltAp-2時(shí)刻室外空氣干球溫度t′ue84aSymboltAp、t′ue84aSymboltAp-1、t′ue84aSymboltAp-2,ue84aSymboltAp時(shí)刻含濕量d′ue84aSymboltAp,ue84aSymboltAp和ue84aSymboltAp-1時(shí)刻太陽(yáng)輻射照度l′ue84aSymboltAp、l′ue84aSymboltAp-1等輸入量,用一個(gè)向量xue84aSymboltAp表示,xTue84aSymboltAp={t′ue84aSymboltAp,t′ue84aSymboltAp-1,t′ue84aSymboltAp-2,d′ue84aSymboltAp,l′ue84aSymboltAp,l′ue84aSymboltAp-1}.ue84aSymboltAp時(shí)刻的建筑物空調(diào)負(fù)荷P′ue84aSymboltAp用yi表示.對(duì)于N個(gè)時(shí)刻,每個(gè)時(shí)刻的相關(guān)參數(shù)作為一個(gè)樣本,則樣本總數(shù)為N,樣本集為{(xi,yi)}Νi=1Ni=1.對(duì)于已確定的樣本集{(xi,yi)}Νi=1Ni=1,參照文獻(xiàn),可進(jìn)行RBF神經(jīng)網(wǎng)絡(luò)建模過(guò)程.將中間層神經(jīng)元的轉(zhuǎn)換函數(shù)取為高斯泛函數(shù),第k個(gè)神經(jīng)元的轉(zhuǎn)換函數(shù)為xi,k=G(|Qi-Τk|)=exp-12|Qi-Τk|σk)2](5)式中:xi,k為第i個(gè)樣本經(jīng)過(guò)第k個(gè)神經(jīng)元的輸出值;Qi為第i個(gè)樣本矢量;Tk為第k個(gè)中間層神經(jīng)元的中心;σk為控制神經(jīng)元相互間重疊程度的參數(shù);|Qi-Τk|=√(Qi-Τk)Τ(Qi-Τk)(6)輸出層的神經(jīng)元實(shí)現(xiàn)了如下的插值過(guò)程:在n維空間Rn中,有N個(gè)不同樣本點(diǎn)的集合{xiue84aSymbolNCpRni=1,2,…,N},并在一維空間R1中相應(yīng)給定N個(gè)實(shí)數(shù)的集合{yiue84aSymbolNCpR1i=1,2,…,N},尋求一函數(shù)f:Rn→R1,使之滿(mǎn)足插值條件f(xi)=yi(i=1,2,…,N).上述基函數(shù)采用了最小二乘法來(lái)實(shí)現(xiàn)輸出層的線(xiàn)性插值,于是,整個(gè)網(wǎng)絡(luò)的映射關(guān)系為f(xi)=n∑k=1ue84aSymbollApkexp-12|Qi-Τk|σk)2](7)按下式把預(yù)測(cè)值進(jìn)行反歸一化,還原為所求的預(yù)測(cè)值,則建立了基于RBF神經(jīng)網(wǎng)絡(luò)的空調(diào)負(fù)荷預(yù)測(cè)模型:?Ρi=Ρmin+f(xi)(Ρmax-Ρmin)(8)2實(shí)例分析2.1示例12.1.1主要建筑構(gòu)件參數(shù)用于驗(yàn)證文中所建立的基于RBF神經(jīng)網(wǎng)絡(luò)的空調(diào)負(fù)荷預(yù)測(cè)模型的建筑物是廣州地區(qū)一幢長(zhǎng)56m、寬20m、高14.4m的4層辦公樓,其外形如圖2所示.該建筑層高3.6m,集中供冷,建筑空調(diào)面積為3581m2.外窗尺寸為2.5m×2.0m,窗臺(tái)離地高度為0.8m.每天空調(diào)運(yùn)行時(shí)間為8:00～21:00,總計(jì)每天運(yùn)行13h.主要建筑構(gòu)件參數(shù)如下:外墻:水泥砂漿+24磚墻+石灰砂漿,熱阻0.3409m2·K/W;內(nèi)墻:水泥砂漿+12磚墻+石灰砂漿,熱阻0.1904m2·K/W;屋頂:水泥砂漿+加氣混凝土+鋼筋混凝土+水泥砂漿,熱阻0.607m2·K/W;地面:水泥砂漿+混凝土+水泥砂漿,熱阻0.026m2·K/W;樓板:水泥砂漿+鋼筋混凝土+水泥砂漿,熱阻0.0975m2·K/W;門(mén):單層實(shí)體木制外門(mén),熱阻0.372m2·K/W;外窗:普通中空玻璃+淺色窗簾,無(wú)外遮陽(yáng),熱阻0.256m2·K/W.2.1.2評(píng)估指標(biāo)采用均方根誤差ΔRMSE和平均相對(duì)誤差ΔMRE兩個(gè)統(tǒng)計(jì)量來(lái)評(píng)價(jià)模型的擬合程度和預(yù)測(cè)效果.式中:?Ρi為Pi的預(yù)測(cè)值.2.1.3模型輸入和計(jì)算方法取7月份(共31天)的數(shù)據(jù)作為訓(xùn)練樣本,分別取5月份(共31天)、6月份(共30天)、8月份(共31天)和10月份(共31天)的數(shù)據(jù)作為檢驗(yàn)樣本,構(gòu)成模型樣本空間.相關(guān)的輸入?yún)?shù)采用廣州典型氣象年的氣象數(shù)據(jù).空調(diào)負(fù)荷的基準(zhǔn)值采用動(dòng)態(tài)負(fù)荷計(jì)算軟件DeST-c的計(jì)算值.2.1.4逐時(shí)負(fù)荷預(yù)測(cè)結(jié)果該辦公樓空調(diào)負(fù)荷的RBF及BP神經(jīng)網(wǎng)絡(luò)模型的訓(xùn)練學(xué)習(xí)過(guò)程、預(yù)測(cè)輸出過(guò)程及誤差曲線(xiàn)如圖3所示.由于兩種模型的訓(xùn)練樣本輸出和實(shí)測(cè)值幾乎完全一致,所以圖3中“o”和“+”基本重疊在基準(zhǔn)線(xiàn)上.用7月份訓(xùn)練后的RBF神經(jīng)網(wǎng)絡(luò)模型和BP網(wǎng)絡(luò)模型對(duì)5月全月的逐時(shí)負(fù)荷進(jìn)行預(yù)測(cè),以驗(yàn)證模型的正確性,其檢驗(yàn)過(guò)程與誤差曲線(xiàn)如圖4所示.同樣,用7月份訓(xùn)練后的RBF和BP網(wǎng)絡(luò)模型對(duì)6月、8月與10月份的逐時(shí)負(fù)荷進(jìn)行預(yù)測(cè),以進(jìn)一步驗(yàn)證模型的正確性,其檢驗(yàn)過(guò)程與誤差曲線(xiàn)如圖5～7所示.訓(xùn)練樣本和檢驗(yàn)樣本預(yù)測(cè)輸出的均方根誤差ΔRMSE和平均相對(duì)誤差ΔMRE見(jiàn)表1.由以上結(jié)果可知,訓(xùn)練過(guò)程輸出誤差均很小,說(shuō)明均獲得了正確的輸入-輸出行為,建立了可表示建筑物空調(diào)負(fù)荷和相關(guān)運(yùn)行參數(shù)關(guān)系的BP和RBF神經(jīng)網(wǎng)絡(luò)模型.在驗(yàn)證模型正確性的過(guò)程中,RBF模型的預(yù)測(cè)誤差均比BP網(wǎng)絡(luò)方法的小.對(duì)于5、6、8、10共4個(gè)月的檢驗(yàn)樣本集,RBF模型預(yù)測(cè)的均方根誤差ΔRMSE分別是BP網(wǎng)絡(luò)方法的63.47%、63.60%、64.28%和63.72%,平均相對(duì)誤差ΔMRE分別是BP網(wǎng)絡(luò)方法的63.60%、63.58%、64.10%和63.98%.這一結(jié)果表明用RBF方法來(lái)預(yù)測(cè)建筑物空調(diào)負(fù)荷,其數(shù)據(jù)吻合程度明顯優(yōu)于BP網(wǎng)絡(luò)方法,均方根誤差ΔRMSE和平均相對(duì)誤差ΔMRE都僅為BP網(wǎng)絡(luò)方法預(yù)測(cè)結(jié)果的64%左右,精度能夠滿(mǎn)足實(shí)際使用的需要.這主要是因?yàn)镽BF網(wǎng)絡(luò)采用局部激勵(lì)函數(shù),具有良好的泛化能力,容易適應(yīng)新數(shù)據(jù),學(xué)習(xí)速度快,其隱層節(jié)點(diǎn)的數(shù)目也在訓(xùn)練過(guò)程中確定,保證其收斂性.2.2示例22.2.1空調(diào)運(yùn)行時(shí)間及主要建筑構(gòu)件參數(shù)為了進(jìn)一步驗(yàn)證已建立的預(yù)測(cè)模型,作者進(jìn)一步將模型應(yīng)用于廣州一棟長(zhǎng)約57m、寬約15m、高26.2m的圖書(shū)館(如圖8所示).該圖書(shū)館共7層,底層層高4.6m,二層以上樓層層高3.6m.采用中央空調(diào),建筑空調(diào)面積3936m2,每天空調(diào)運(yùn)行時(shí)間為8:00～18:00,總計(jì)每天運(yùn)行10h.主要建筑構(gòu)件參數(shù)如下:外墻:水泥砂漿+24磚墻+石灰砂漿,熱阻0.3409m2·K/W;內(nèi)墻:水泥砂漿+陶粒混凝土+水泥砂漿,熱阻0.430m2·K/W;屋頂:種植屋面,熱阻1.146m2·K/W;地面:水泥砂漿+混凝土+水泥砂漿,熱阻0.026m2·K/W;樓板:水泥砂漿+鋼筋混凝土+水泥砂漿,熱阻0.0975m2·K/W;門(mén):圖書(shū)館一層外門(mén)為厚度為18mm的玻璃外門(mén),熱阻0.031m2·K/W,二層以上外門(mén)及內(nèi)門(mén)均為單層實(shí)體木制門(mén),熱阻0.372m2·K/W;外窗:普通單層玻璃+淺色窗簾,無(wú)外遮陽(yáng),熱阻0.156m2·K/W.2.2.2rbf模型預(yù)測(cè)的均方根誤差提升取6月份(共30天)的數(shù)據(jù)作為訓(xùn)練樣本,分別取5月份(共31天)、7月份(共31天)、8月份(共31天)和10月份(共31天)的數(shù)據(jù)作為檢驗(yàn)樣本,構(gòu)成模型樣本空間.與實(shí)例1相同,采用均方根誤差和平均相對(duì)誤差兩個(gè)統(tǒng)計(jì)量來(lái)評(píng)價(jià)模型的擬合程度和預(yù)測(cè)效果.限于篇幅,文中僅給出訓(xùn)練樣本和檢驗(yàn)樣本預(yù)測(cè)輸出的均方根誤差ΔRMSE和平均相對(duì)誤差ΔMRE,見(jiàn)表2.由表2可見(jiàn),對(duì)于5、7、8、10月份的檢驗(yàn)樣本集,RBF模型預(yù)測(cè)的均方根誤差ΔRMSE分別為BP網(wǎng)絡(luò)方法的63.99%、63.86%、64.02%和64.54%,平均相對(duì)誤差ΔMRE分別為BP網(wǎng)絡(luò)方法的63.75%、63.46%、64.18%和64.46%.這一結(jié)果進(jìn)一步表明RBF方法可更為有效地預(yù)測(cè)建筑物空調(diào)負(fù)荷,精度能夠滿(mǎn)足實(shí)際使用的需要.3空調(diào)負(fù)荷預(yù)測(cè)系統(tǒng)的組成為了使建筑物空調(diào)負(fù)荷智能預(yù)測(cè)模型能夠在現(xiàn)實(shí)中得到應(yīng)用,筆者構(gòu)建了一套智能預(yù)測(cè)軟件構(gòu)架.目前,由于時(shí)間以及技術(shù)條件等限制,預(yù)期的現(xiàn)場(chǎng)預(yù)測(cè)系統(tǒng)還沒(méi)有完全建立起來(lái),只能根據(jù)目前的研究結(jié)果對(duì)建筑物空調(diào)負(fù)荷預(yù)測(cè)系統(tǒng)做一個(gè)大致的構(gòu)建.該系統(tǒng)由室外氣象參數(shù)、建筑物參數(shù)設(shè)置、空調(diào)運(yùn)行時(shí)間設(shè)置、智能預(yù)測(cè)模型選擇和計(jì)算及結(jié)果顯示5個(gè)模塊構(gòu)成,如圖9所示.其中,顯示逐時(shí)空調(diào)負(fù)荷的預(yù)測(cè)曲線(xiàn)的界面如圖10所示.4空調(diào)負(fù)荷預(yù)測(cè)由于RBF神經(jīng)網(wǎng)絡(luò)具有學(xué)習(xí)速度快、精度高、全局最優(yōu)及泛化性好等特點(diǎn),將其應(yīng)用于建筑物空調(diào)負(fù)荷的預(yù)測(cè)中,建立基于RBF神經(jīng)網(wǎng)絡(luò)的空調(diào)負(fù)荷預(yù)測(cè)模型,可將建筑物空調(diào)負(fù)荷特性與各個(gè)相關(guān)參數(shù)之間的復(fù)雜關(guān)系及參數(shù)之間的耦合關(guān)系簡(jiǎn)單地表示出來(lái),準(zhǔn)確有效地對(duì)空調(diào)負(fù)荷進(jìn)行預(yù)測(cè).將RBF建筑物空調(diào)負(fù)荷預(yù)測(cè)方法應(yīng)用于廣州地區(qū)一棟辦公樓和一棟圖書(shū)館的空調(diào)負(fù)荷預(yù)測(cè)中,RBF模型預(yù)測(cè)的均方根誤差ΔRMSE和平均相對(duì)誤差ΔMRE都僅是BP網(wǎng)絡(luò)方法的64%左右,表