金碧瑤，李占培，劉廷章，張修紅，閆 斌，張穎婍

（上海大學機電工程與自動化學院，上海市電站自動化技術重點實驗室，上海 200072）

基于傳熱機理的實時冷負荷預測

金碧瑤，李占培，劉廷章，張修紅，閆 斌，張穎婍

（上海大學機電工程與自動化學院，上海市電站自動化技術重點實驗室，上海 200072）

針對現(xiàn)有冷負荷預測方法不能有效解決動態(tài)實時性的問題，提出一種基于傳熱機理的實時冷負荷預測方法。根據(jù)諧波反應法建立針對室外氣象因子的實時冷負荷預測簡化模型，用最小二乘法辨識模型參數(shù)，從而降低計算量。實驗結果表明該模型能較好地反映室外氣象因子的實時動態(tài)性，具有較高準確度。

冷負荷預測；動態(tài)實時性；室外氣象因子；模型

0 引 言

能耗的迅速增加是造成我國電力供應緊張的一個重要因素。空調(diào)的能耗約占建筑能耗的40%～55%[1]，可見空調(diào)系統(tǒng)的節(jié)能是建筑節(jié)能的重點之一。準確計算空調(diào)的實時負荷有助于調(diào)整空調(diào)的實際運行，同時盡可能地減少能源消耗。

空調(diào)冷負荷[2-3]形成過程中，室內(nèi)熱量的產(chǎn)生遵循一定散熱規(guī)律，可通過簡化的工程方法進行計算。而室外環(huán)境對房間冷負荷形成的影響有圍護結構、窗戶傳熱兩個方面。室外環(huán)境包括室外溫度、室外空氣濕度、太陽輻射、室外風速等，這些因素具有不確定的變化規(guī)律，是研究冷負荷計算的難點。現(xiàn)有的空調(diào)負荷計算方法[4-8]可以分為兩類：面向空調(diào)工程設計的負荷計算和面向空調(diào)運行的負荷預測。前者準確性較高，但計算量大；后者建模時間較長，模型建立后可以實時快速地計算，但無法涵蓋室內(nèi)外動態(tài)變化的因素。

本文提出一種針對室外氣象因子的基于傳熱機理的實時冷負荷預測方法。該方法融合了兩類方法的優(yōu)點，能有效地解決動態(tài)實時性問題，同時具有較好的準確性。

1 冷負荷模型的建立

在不考慮內(nèi)擾的情況下，基于傳熱機理的冷負荷形成過程如圖1所示。房間冷負荷主要由通過圍護結構得到的冷負荷（CLQ1）和經(jīng)窗戶得熱形成的冷負荷（CLQ2）組成。因此，某時刻的冷負荷為

圖1 建模總體方案

1.1 通過圍護結構得熱量形成冷負荷

圍護結構主要是指墻體、屋頂?shù)?。在諧波反應法[9-10]中，通過圍護結構形成的冷負荷由對流熱成分和輻射熱成分組成。只考慮傅里葉的低階級數(shù)時，某時刻通過圍護結構得熱量形成冷負荷CLQ1為

式中：K——圍護結構傳熱系數(shù)，W/（m2·K）；

F——墻體或屋頂?shù)拿娣e，m2；

tN——室內(nèi)空氣溫度，℃；

βd——對流成分占總得熱量的比例；

βf——輻射成分占總得熱量的比例；

αN——圍護結構內(nèi)表面放熱系數(shù)，W/（m2·K）；

ΔtZ·1——圍護結構的一階外擾；

ν1——圍護結構對一階綜合溫度擾量的衰減度；

μ1——房間對一階墻體或屋頂傳導得熱中輻射擾量的衰減度。

室外空氣綜合溫度tz受室外空氣溫度tW和太陽總輻射I的影響而室外空氣綜合溫度是以24h為周期的，因此其平均值以24h為單位取平均。將式（2）中的室外平均溫度用室外空氣溫度和太陽總輻射來表示并簡化可得：

式中：ρ——圍護結構外表面對太陽輻射熱的吸收系數(shù)；

αW——圍護結構內(nèi)外表面的總放熱系數(shù)，W/（m2·K）；

外空氣溫度和總輻射，單位分別為℃，W/m2。

1.2 通過窗戶得熱量形成冷負荷

通過窗戶得熱量形成的冷負荷包括瞬變傳導形成的冷負荷（CLQA）和日射形成的冷負荷（CLQB）即：

某時刻瞬變傳導形成的冷負荷為

式中A1為一階室外空氣溫度的波幅，℃。

相應的日射形成的冷負荷為

式中：Cs——窗玻璃的遮陽系數(shù)；

C1——窗內(nèi)遮陽設施的遮陽系數(shù)；

B1——一階日射得熱因數(shù)諧波的波幅。

日射得熱取決于很多因素，從太陽輻射方面來說，輻射強度和入射角均依緯度、月份、日期、時間的不同而不同。從窗戶本身的結構來說，它隨玻璃的光學性能，窗戶結構、特性、是否有遮陽裝置等而異。此外，還與內(nèi)外放熱系數(shù)有關。為了計算方便，這里主要選擇太陽直射輻射強度IC·Z和太陽散射輻射強度IS這兩個因素。因此式（7）和式（8）可表示為

αD——玻璃窗直射吸收率；

αS——玻璃窗散射吸收率。

1.3 房間總的冷負荷

綜合上述分析，某時刻房間的冷負荷可以表示為與室外空氣溫度tW、太陽總輻射I、太陽直射輻射強度IC·Z及太陽散射輻射強度IS有關的關系式。而太陽的總輻射、散射和直射是存在一定聯(lián)系的，知道其中兩個量，必可推倒出第3個量。因此，將室外氣象因子減為室外空氣溫度、總輻射和散射3個量，可用矩陣表示：

2 模型參數(shù)的辨識

在實驗過程中，由于測量儀器有限，無法實時采集太陽總輻射、散射等數(shù)值，并且冷負荷又難以測量，而DeST[11]是一款可以實現(xiàn)房間冷負荷模擬的仿真軟件。因此，本文以上海某建筑為原型建立實際房間的模型，選取特定房間，通過仿真得到7月的冷負荷數(shù)據(jù)。將DeST中的典型氣象值和仿真得到的冷負荷數(shù)據(jù)中前25日的逐時值作為模型辨識數(shù)據(jù)，剩余數(shù)據(jù)作為模型驗證數(shù)據(jù)。

本文采用最小二乘法在Matlab平臺上進行辨識。將室外空氣溫度、太陽總輻射、散射以及模擬房間冷負荷數(shù)據(jù)的24 h樣本序列輸入Matlab程序中進行計算。7月份的部分樣本數(shù)據(jù)如表1所示。辨識的結果為

3 模型驗證

由辨識得到的模型參數(shù)結合7月26日至7月31日的氣象數(shù)據(jù)，可以得到一組計算負荷，并與DeST的仿真負荷進行比較驗證，其結果如圖2所示。橫軸為時間，每天從0時刻到23時刻反復進行；縱軸為冷負荷值，單位是kW。藍色為DeST的仿真負荷曲線，紅色為利用模型計算得到的冷負荷曲線。雖然在局部點上，計算負荷與仿真負荷存在差距，但二者在總體趨勢上是一致的。

圖2 模型計算結果與DeST仿真結果對比

圖3為相對誤差統(tǒng)計圖，模型的計算與仿真結果的最大相對誤差為0.084%，說明該模型計算結果在合理的范圍內(nèi)，從而驗證了模型的準確性。

圖4為冷負荷與氣象因子的趨勢圖。從圖中可以看出，影響冷負荷的最主要因子是室外空氣溫度，同時也驗證了模型能較好地反映室外氣象因子的動態(tài)實時性。

表1 7月份的部分樣本數(shù)據(jù)

圖3 相對誤差統(tǒng)計

圖4 冷負荷與氣象因子的趨勢圖

4 結束語

本文提出的基于傳熱機理的實時冷負荷模型方法能夠在空調(diào)運行工況下，進行房間冷負荷計算，為優(yōu)化房間空調(diào)供冷量控制策略提供了理論依據(jù)。進一步分析了房間冷負荷形成的重要影響因素，實現(xiàn)空調(diào)冷負荷簡化計算。該模型能有效地解決動態(tài)實時性問題，同時又保證其準確度。

[1]清華大學建筑節(jié)能研究中心.中國建筑節(jié)能年度發(fā)展研究報告2009[M].北京：中國建筑工業(yè)出版社，2009：92-127.

[2]Fouda A，Melikyan Z，Mohamed M A，et al.A modified method of calculating the heating load for residential buildings[J].Energy and Buildings，2014（75）：170-175.

[3]Yao Y，Chen J. Global optimization of a central air-conditioning system using decomposition coordination method[J].Energy Build，2010，42（5）：570-583.

[4]吳杰.冰蓄冷空調(diào)系統(tǒng)負荷預測模型和系統(tǒng)優(yōu)化控制研究[D].杭州：浙江大學，2002.

[5]Kusiak A，Li M，Zhang Z.A data-driven approach for steam load prediction in building[J].Applied Energy，2010，87（3）：925-933.

[6]Yun K，Luck R，Mago P J，et al.Building hourly thermal load prediction using an indexed ARX model[J]. Energy Build，2012（54）：225-233.

[7]Yun K，Cho H，Luck R，et al.Real-time combined heat and power operational strategy using a hierarchical optimization [J]. Proceedings of the Institution of Mechanical Engineers，Part A：Journal of Power and Energy，2011，225（4）：403-412.

[8]Yool P，Jong I P.Analysis of energy performance by window directions in apartment[C]∥ Proceedings of Summer Conference of Air-Conditioning Refrigeration Engineering，2012：1041-1044.

[9]趙榮義，范存養(yǎng)，薛殿華，等.空氣調(diào)節(jié)[M].北京：中國建筑工業(yè)出版社，2008：20-50.

[10]蒲凈，劉桂蓮，閻哲，等.建筑空調(diào)逐時負荷的確定[J].華北電力大學學報，2009，36（3）：20-22.

[11]江億.建筑環(huán)境系統(tǒng)模擬分析方法——DeST[M].北京：中國建筑工業(yè)出版社，2006：117-141.

Forecasting of real-time cooling load based on heat transfer mechanism

JIN Biyao，LI Zhanpei，LIU Tingzhang，ZHANG Xiuhong，YAN Bin，ZHANG Yingqi
（Shanghai Key Laboratory of Power Station Automation Technology，School of Mechanical Engineering and Automation，Shanghai University，Shanghai 200072，China）

This paper proposes a methodology in predicting real-time cooling load based on heattransfer mechanism，which solves the problem that the current calculation and prediction on cooling load cannot work out the dynamic real-time performance efficiently.Considering outdoor meteorologicalfactors，the cooling load simplified modelis builtvia harmonic reaction and identified by least square method，which reduces the amount of computation.The result indicates that the model improves dynamic real-time performance of outdoor meteorological factors as well as accuracy.

cooling load forecasting；dynamic real-time；outdoor meteorological factors；model

A

：1674-5124（2015）07-0112-04

10.11857/j.issn.1674-5124.2015.07.026

2014-11-12；

：2015-01-16

國家自然科學基金項目（61273190）

金碧瑤（1989-），女，浙江安吉市人，碩士研究生，專業(yè)方向為空調(diào)節(jié)能控制。