袁 斌，方芩璐，羅滇生，王 娟

（1.湖南大學電氣與信息工程學院，長沙 410082；2.四川省成都電業(yè)局，成都 610016；3.63663部隊技術部，烏魯木齊 841700）

電力負荷預測，尤其是短期電力負荷預測在實時控制和保證電力系統(tǒng)經(jīng)濟、安全和可靠運行方面起著重要作用，對于系統(tǒng)運行和生產消費具有非常重大的意義［1，2］。近年來，大多數(shù)大中城市的負荷呈現(xiàn)出居民類負荷和第三產業(yè)負荷在社會總負荷中比重逐年上升的趨勢。而這些負荷與天氣狀況尤其是氣溫存在密切關系，為氣溫敏感型負荷。因此要提高預測精度，關鍵是如何更加合理地考慮氣象因素，尤其是氣溫對負荷的影響［3～5］。

目前的預測方法中，模糊線性回歸法、人工神經(jīng)網(wǎng)絡法、專家系統(tǒng)法等方法對氣象因素進行了很好的量化和估測，對傳統(tǒng)預測方法進行了改進。文獻［6］提出了一種利用多地區(qū)氣溫因素來進行電力負荷預測的方法，通過提高預測模型中輸入（最高氣溫）與輸出（負荷）的關聯(lián)度來提高負荷預測精度。

一些系統(tǒng)事實上已經(jīng)把這些預測模型用到了負荷調度中心。湖南省電力公司已經(jīng)安裝了用人工神經(jīng)網(wǎng)絡方法進行96點短期負荷預測的負荷預測系統(tǒng)，并在超短期和短期負荷預測中取得了可喜的成果。不過該系統(tǒng)僅用長沙一個地區(qū)氣象因素代替整個湖南省的氣象因素，當湖南省其它地區(qū)氣象因素不同于長沙市的氣象因素時可能會產生很大的負荷預測偏差。另外，在系統(tǒng)輸入端氣溫因素只采用了最高氣溫，而輸出端的預測值是96個時間斷面點的負荷，因此在某些時間點會出現(xiàn)負荷與氣溫關聯(lián)度不高的情況。為解決這些問題，本文提出了利用粒子群算法動態(tài)地優(yōu)化湖南省各個地區(qū)氣溫比率權值，提高預測系統(tǒng)輸入輸出關聯(lián)度的方法來提高負荷預測的精度。

1 負荷預測模型

預測系統(tǒng)的神經(jīng)網(wǎng)絡是三層前向網(wǎng)絡如圖1

圖1 人工神經(jīng)網(wǎng)絡結構Fig.1 Structure of ANN

所示。目標負荷是下一天的96個時間斷面點的負荷。預測系統(tǒng)的輸入端讀取歷史負荷數(shù)據(jù)庫數(shù)據(jù)，同時輸入天氣、氣溫、空氣濕度等信息。湖南省氣候四季分明，在不同的季節(jié)電力負荷和天氣條件的關聯(lián)度不同。根據(jù)聚類分析思想，在四個季節(jié)分別建立網(wǎng)絡模型，并且區(qū)分工作日和節(jié)假日模型［7］。

按照平均負荷大小和負荷曲線形狀按來劃分季節(jié)類型，劃分如下：3月1日到6月30日為春季；7月1日到9月15日為夏季；9月16日到10月30日為秋季；11月1日到3月31日為冬季。每個星期利用一次歷史實際負荷數(shù)據(jù)對神經(jīng)網(wǎng)絡進行訓練。

2 對輸入輸出關聯(lián)度的改進

2.1 全局氣溫的選取

湖南省幅員遼闊，從最東部到最西部大約有530 km，從最北部到最南部大約有610 km，這個地域有平原和高山。不同地區(qū)的氣象因素都不一樣。因此在選取氣溫因素時，應綜合考慮全省各個地區(qū)的氣溫信息。這里有多種方法，如將各個地區(qū)氣溫簡單的算術平均和根據(jù)各個地區(qū)負荷需求比率來加權平均等。然而這些都不是最好和最有效的方法。本文提出的利用粒子群算法優(yōu)化幾個地區(qū)的氣溫權值比率，從而得出最佳的全局氣溫。按湖南省電力公司行政區(qū)劃為14個地區(qū)，各個地區(qū)電力負荷需求比率如表1所示。

表1 各個地區(qū)的負荷需求比率Tab.1 Load demand rate of different areas

2.2 動態(tài)氣溫的采用

負荷預測系統(tǒng)將一天24 h分為96個時間斷面，每隔15 min為一個時間斷面點，分別進行預測，從而得到下一天的96個時間斷面點的負荷。因此，取動態(tài)的氣溫信息能更準確的反應負荷與氣溫的關聯(lián)度?？紤]到氣溫的變化的特點，以及氣溫數(shù)據(jù)的采集成本等原因，取每個整點的氣溫共24點（或者每兩個小時的氣溫共12個點）。每個時間斷面點的負荷對應這個時間段的氣溫。

2.3 目標函數(shù)和約束

有著高的輸入輸出關聯(lián)度的輸入數(shù)據(jù)一般更加能實現(xiàn)精確地負荷預測。因此，最高關聯(lián)度的最好加權比率就能實現(xiàn)最精確地負荷預測。本文中，收集到湖南省14個地區(qū)的歷史每日96點時間斷面的負荷數(shù)據(jù)以及歷史氣溫信息。目標函數(shù)的關聯(lián)度r2是在每個時間段，14個地區(qū)的氣溫的權重平均和負荷的近似二次函數(shù)曲線。目標函數(shù)和約束如下所示。

其中r2為關聯(lián)度系數(shù)；xn為氣溫權值比率。

2.4 粒子群算法優(yōu)化氣溫權值比率

粒子群優(yōu)化PSO（particle swarm optimization）算法［8，9］是一種基于集群智能的隨機優(yōu)化算法，其基本思想是：優(yōu)化問題的每一個解稱為一個粒子。定義一個符合度函數(shù)來衡量每個粒子解的優(yōu)越程度。每個粒子根據(jù)自己和其它粒子的“飛行經(jīng)驗”群游，從而達到從全空間搜索最優(yōu)解的目的。具體搜索過程如下。

每個粒子在解空間中同時向兩個點接近，第一個點是整個粒子群中所有粒子在歷代搜索過程中所達到的最優(yōu)解，被稱為全局最優(yōu)解g＊；另一個點則是每個粒子在歷代搜索過程中自身所達到的最優(yōu)解，這個解被稱為個體最優(yōu)解p＊。每個粒子表示在n維空間中的一個點。

式（3）表示第i個粒子的氣溫權值，其個體最優(yōu)解（第i個粒子最小適應值所對應的解）為

全局最優(yōu)解（整個粒子群在歷代搜索過程中最小適應值所對應的解）為

而xi的第k次迭代的修正量（粒子移動的速度）表示為

其計算公式為

其中i＝1，2，…，M，M為種群數(shù)；d＝1，2，…，N，N為粒子大??；c1、c2是大于0的權重因子，r1、r2［0，1］之間的隨機數(shù)；w是慣性權重系數(shù)；k是迭代次數(shù)。

粒子群優(yōu)化算法步驟如下：

步驟1設定種群大小為10，粒子大小為14。初始化粒子向量的粒子（各個地區(qū)氣溫權值和速度，設定迭代次數(shù)80次；

步驟2計算每個粒子在當前狀態(tài)下的適應函數(shù)值（氣溫和負荷之間的關聯(lián)度）pi；

步驟3將步驟2中計算的適應函數(shù)值pi與自身的優(yōu)化解進行比較，如果｜pi｜＜｜｜，則用新的適應函數(shù)值取代前一輪的優(yōu)化解，用新的粒子取代前一輪粒子，即←pi←xi；

步驟5完成以上的計算后，再進行新一輪的計算，按式（5）將粒子進行移動，從而產生新的粒子（即新解），返回步驟2。直至完成設定的迭代次數(shù)或滿足事先給定的精度要求為止。

2.5 優(yōu)化方法的流程

優(yōu)化方法的流程圖如圖2所示。所有的解決方案通過下面的步驟進行評估。

列出相應時刻14個地區(qū)氣溫，14個地區(qū)的氣溫和各自的權重比率相乘后相加得出一個全局氣溫，作為預測該個時間斷面負荷的輸入氣溫。

全局氣溫和斷面負荷之間的關系可以擬合成一條二次曲線。

計算權值氣溫和二次曲線函數(shù)（目標函數(shù)）之間的關聯(lián)系數(shù)。

這時，更新權值比率。

圖2 氣溫權值比率優(yōu)化流程Fig.2 Flow of temperature weight rate improvement

3 仿真結果

對本文所提出的方法和傳統(tǒng)的方法進行了比較，體現(xiàn)出該方法有效性。仿真的數(shù)據(jù)利用湖南省2008年的日負荷96點數(shù)據(jù)，取11∶00時刻的歷史負荷值，驗證氣溫與負荷的關聯(lián)度。仿真結果如表2所示。方案1中用的是長沙市的日最高氣溫，沒有取動態(tài)氣溫，也不是各地區(qū)加權平均氣溫（方案1曾是湖南省電力公司一種傳統(tǒng)負荷預測方法）。在方案2中，權值比率等于每個地區(qū)的電力負荷需求比率。本文所提的方法是方案3，用粒子群算法以最大化目標函數(shù)來優(yōu)化氣溫權值比率。方案3所計算出來的關聯(lián)度r2是所有方案中最高的。

圖3、圖4分別顯示了方案1中負荷與日最高氣溫的二次曲線函數(shù)擬合和3次曲線函數(shù)擬合的關聯(lián)度；圖5、圖6分別顯示了方案3中負荷與加權優(yōu)化氣溫的二次曲線函數(shù)擬合和3次曲線函數(shù)擬合的關聯(lián)度。由圖可見，方案1的數(shù)據(jù)和方案3比較起來分散得更加大一些；二次曲線函數(shù)擬合和3次曲線函數(shù)擬合的結果都證明了方案3的關聯(lián)度更高。

表2 三種方案的關聯(lián)度對比Tab.2 Contrast of three methods

圖3 方案1二次曲線函數(shù)擬合結果（A點坐標（26.1，9270））Fig.3 Quadratic curve fitting results of scheme 1（Point A＇s coordinate is（26.1，9270））

圖4 方案1的三次曲線函數(shù)擬合結果Fig.4 Cubic curve fitting results of scheme 1

圖3和圖5中的A數(shù)據(jù)為同一時刻采用方案1和方案3分別驗證的結果。圖3中長沙市當日最高氣溫是26.1℃；當日湖南省其它地區(qū)的氣溫比較高，圖5中加權優(yōu)化后的氣溫為28.5℃。優(yōu)化后的氣溫與負荷的關聯(lián)度更高，以這個氣溫預測的負荷結果更加準確。

圖5 方案3二次曲線函數(shù)擬合結果（A點坐標（28.5，9270））Fig.5 Quadratic curve fitting results of scheme 3（Point A＇s coordinate is（28.5，9270））

圖6 方案3的三次曲線函數(shù)擬合結果Fig.6 Cubic curve fitting results of scheme 3

4 結論

將用粒子群算法以最大化目標函數(shù)來優(yōu)化加權比率得到的全局氣溫作為負荷預測模型的氣溫輸入信息，用改進后的模型返回去驗證2008年的負荷預測結果，事實證明本文方法氣溫與負荷的關聯(lián)度更高，預測更加準確，更加合理。

［1］ 吳杰康，陳明華，陳國通，等（Wu Jiekang，Chen Minghua，Chen Guotong，et al）.基于PSO的模糊神經(jīng)網(wǎng)絡短期負荷預測（Fuzzy neural network model based on particle swarm optimization for short－term load forecasting）［J］.電力系統(tǒng)及其自動化學報（Proceedings of the CSU－EPSA），2007，19（1）：63－67.

［2］ 邵偉，楊敏虹，徐政，等（Shao Wei，Yang Minhong，Xu Zheng，et al）.發(fā)電市場中負荷預測的功能及其方法（Functions and methods of load forecasting in power generation market）［J］.中國電力（Electric Power），2000，33（11）：54－57.

［3］ 王波，邰能靈，翟海青，等（Wang Bo，Tai Nengling，Zhai Haiqing，et al）.基于混合粒子群算法的短期負荷預測模型（Hybrid optimization method based on evolutionary algorithm and particle swarm optimization for short－term load forecasting）［J］.電力系統(tǒng)及其自動化學報（Proceedings of the CSU－EPSA），2008，20（3）：50－55.

［4］ 羅鳳章，王成山，肖峻，等（Luo Fengzhang，Wang Chengshan，Xiao Jun，et al）.計及氣溫因素的年度負荷預測修正方法（Exploration on calibration method for yearly load forecasting considering temperature）［J］.電力系統(tǒng)及其自動化學報（Proceedings of the CSU－EPSA），2009，21（3）：75－79.

［5］ 楊旭，劉自發(fā)，張建華，等（Yang Xu，Liu Zifa，Zhang Jianhua，et al）.考慮氣候因素的城市電力負荷預測（Urban power load forecasting considering climate factor）［J］.現(xiàn)代電力（Modern Electric Power），2008，25（3）：28－32.

［6］ Mizutani Akihiro，Yukawa Tetsuya，Numa Kazuyuki，et al.Improvement of input－output correlations of electric power load forecasting by scatter search［C］∥The 13th International Conference on Intelligent Systems Application to Power Systems，Arlington，USA：2005.

［7］ 楊爭林，唐國慶，宋燕敏，等（Yang Zhenglin，Tang Guoqing，Song Yanmin，et al）.改進的基于聚類分析的超短期負荷預測方法（Improved cluster analysis based ultra－short term load forecasting method）［J］.電力系統(tǒng)自動化（Automation of Electric Power Sys－tems），2005，29（24）：83－86，97.

［8］ 吳昌友，王福林，董志貴，等（Wu Changyou，Wang Fulin，Dong Zhigui，et al）.改進粒子群優(yōu)化算法在電力負荷組合預測模型中的應用（Application of improved particle swarm optimization in power load combination forecasting model）［J］.電網(wǎng)技術（Power System Technology），2009，33（2）：27－30.

［9］ 許傳敏，楊洪耕（Xu Chuanmin，Yang Honggeng）.電壓波動估計的混沌簡化粒子群算法（Chaotic simple particle swarm optimization algorithm for voltage fluctuation measurement）［J］.電力系統(tǒng)及其自動化學報（Proceedings of the CSU－EPSA），2009，21（3）：18－22.