鄭曉菁

（安徽理工大學電氣與信息工程學院，安徽 淮南 232001）

基于人工蜂群優(yōu)化法的多區(qū)域電力系統(tǒng)經(jīng)濟調(diào)度
*

鄭曉菁

（安徽理工大學電氣與信息工程學院，安徽 淮南 232001）

針對多區(qū)域電力系統(tǒng)經(jīng)濟調(diào)度問題，在滿足聯(lián)絡線傳輸限制、多種燃料特征、閥點效應和禁止運轉(zhuǎn)區(qū)的約束條件下，綜合考慮多區(qū)域電力負載成本最小的要求，建立數(shù)學計算模型，利用人工蜂群優(yōu)化法快速地尋找全局最優(yōu)解。通過兩個不同規(guī)模、不同程度復雜性的仿真測試系統(tǒng)進行計算，結果驗證了所提算法的可行性?？紤]獲得解的質(zhì)量，將人工蜂群優(yōu)化算法與DE、EP、RCGA算法進行對比分析，結果表明所提算法在實際電力系統(tǒng)中解決多區(qū)域經(jīng)濟分配問題具有有效性和優(yōu)越性。

多區(qū)域經(jīng)濟分配；人工蜂群優(yōu)化算法；聯(lián)絡線約束；電力系統(tǒng)；功率平衡約束

1 引言

在電力系統(tǒng)運行中，經(jīng)濟調(diào)度［1，2］是一個重要的優(yōu)化問題，其目標是在滿足單區(qū)域范圍內(nèi)受到各種自身限制因素的影響前提下，使發(fā)電總成本最小化。而多區(qū)域電力系統(tǒng)經(jīng)濟調(diào)度通常是將發(fā)電機組劃分為幾個相互連接的發(fā)電區(qū)域，在滿足電力需求、電機特征等約束下，尋求系統(tǒng)的發(fā)電能力和各區(qū)域之間的電力交換，從而使總體發(fā)電量成本最小化。

國內(nèi)外學者對于經(jīng)濟調(diào)度問題的研究較多，但并未考慮電力傳輸約束的影響。文獻［3］闡述了區(qū)域間電力輸入輸出限制的經(jīng)濟調(diào)度問題，研究提出了多區(qū)域發(fā)電計劃的完整公式和框架。Romano R等［4］對于多區(qū)域電力系統(tǒng)約束的經(jīng)濟調(diào)度提出了丹沃爾夫分解原理。Doty K W和McEntire P L［5］運用空間動態(tài)規(guī)劃法研究了多區(qū)域經(jīng)濟調(diào)度問題，獲得全局優(yōu)化最優(yōu)解。Streiffert D［6］提出了用網(wǎng)路流模型解決受傳輸限制的多區(qū)域經(jīng)濟調(diào)度問題。Yalcinoz T和Short M J［7］使用霍普菲爾神經(jīng)網(wǎng)路法解決多區(qū)域經(jīng)濟調(diào)度問題。Jayabarathi T等［8］使用進化規(guī)劃法解決受聯(lián)絡線限制的多區(qū)域經(jīng)濟調(diào)度問題。文獻［9］研究了多區(qū)域經(jīng)濟環(huán)境調(diào)度問題。文獻［10］采用隨時間變化的變異差分進化法處理儲備約束多區(qū)域經(jīng)濟調(diào)度問題。

由于智能優(yōu)化算法對問題特征沒有明確要求，并具有很強的全局搜索能力，近年來，在經(jīng)濟調(diào)度問題上得到廣泛應用，如實數(shù)編碼遺傳算法［11］、粒子 群 算 法［12］、差 分 進 化 算法［13］和 進 化 規(guī) 劃法［14］等。

人工蜂群優(yōu)化算法［15］是 一種新 型的 群智 能優(yōu)化算法，該算法模擬蜜蜂的采蜜行為，通過獨特的角色分配，可簡便、有效地解決復雜的組合優(yōu)化問題。

本文針對多區(qū)域經(jīng)濟調(diào)度問題，提出了人工蜂群優(yōu)化算法，并應用于考慮運行區(qū)域、傳輸損失、閥點效應和多種燃料特征的經(jīng)濟調(diào)度問題。

2 問題描述

多區(qū)域電力系統(tǒng)經(jīng)濟調(diào)度的目標是在滿足功率平衡約束、發(fā)電限制及傳送線容量限制等約束條件下，使得供應給整個區(qū)域內(nèi)電力負載的總費用最小。

2.1 目標函數(shù)

目標函數(shù)為：

其中，Eij（Pij）為在區(qū)域i中第j個發(fā)電機組的成本函數(shù)，通常用二次多項式表示；aij、bij和cij為在區(qū)域i中第j個發(fā)電機組的成本系數(shù)；N為區(qū)域數(shù)量；Mi為在區(qū)域i中的發(fā)電機數(shù)量；Pij為在區(qū)域i中第j個發(fā)電機的實際輸出功率。然而，在汽輪機進氣閥突然開啟時出現(xiàn)的拔絲現(xiàn)象會在機組耗量曲線上疊加一個脈沖效應，產(chǎn)生閥點效應［16］。當計及發(fā)電機組的閥點效應時，目標函數(shù)（1）應變?yōu)椋?/p>

其中，dij、eij為在區(qū)域i內(nèi)與第j個發(fā)電機閥點效應的價值系數(shù)。

由 于發(fā) 電 機 燃 料 來源是 多 樣 化 的［17］，因 此 通過幾個分段二次函數(shù)的正弦項來反映燃料類型的變化，以確定最經(jīng)濟燃燒燃料。對于同時考慮 NF種燃料和閥點效應，第i個發(fā)電機的發(fā)電成本為：

2.2 約束條件

函數(shù)Et要達到最小值受到以下條件約束。

2.2.1 實際功率平衡約束

其中：

其中，PDi為區(qū)域i的實際功率需求；PLi為區(qū)域i的傳輸損失；Tik為從區(qū)域i傳到區(qū)域k的聯(lián)絡線實際功率，當功率從i傳輸?shù)絢時，Tik取正值；反之取負值，Bilj為傳輸損失系數(shù)。

2.2.2 聯(lián)絡線性能約束

基于安全考慮，從區(qū)域i到區(qū)域k的聯(lián)絡線的實際傳輸功率Tik不應超過聯(lián)絡線的傳輸能力，應滿足：

2.2.3 發(fā)電機組約束

發(fā)電機組的功率上下限制為：

2.2.4 禁止運轉(zhuǎn)區(qū)

對于存在禁止運轉(zhuǎn)區(qū)的發(fā)電機組i的輸出功率限制約束為：

3 平衡發(fā)電機功率確定

假設有（Mi—1）個發(fā)電機的電力負荷是已知的，第Mi個發(fā)電機（平衡發(fā)電機）的電力負荷可表示為：

此時，傳輸損失PLi可以用一個包括平衡發(fā)電機的所有發(fā)電機輸出電力的函數(shù)來表示：

將式（10）代入式（9），式（9）變?yōu)椋?/p>

利用標準代數(shù)的方法，求解式（11）可得到平衡發(fā)電機的負荷。

4 人工蜂群優(yōu)化算法在多區(qū)域經(jīng)濟調(diào)度問題中的應用

4.1 人工蜂群優(yōu)化法

人工蜂群算法是Karaboga D等人于2005年提出的一種模擬蜜蜂采蜜行為的隨機搜索優(yōu)化算

法。此算法中，模擬三種蜜蜂的搜索行為：采蜜蜂、觀察蜂和偵察蜂。采蜜蜂在記憶中選擇一個鄰近蜜源；通過蜂巢內(nèi)的觀察蜂共享蜜源信息，再由觀察蜂在鄰近的蜜源內(nèi)選擇一個蜜源，此時采蜜蜂放棄蜜源，變成偵查蜂，并開始全局隨機搜索質(zhì)量更高的新蜜源。每只蜜蜂對應了一個解，采蜜蜂代表構成當前種群的現(xiàn)有解；觀察蜂代表潛在的鄰域搜索解，有機會進入種群成為現(xiàn)有解；偵查蜂則代表全局隨機搜索解，可以代替廢棄的現(xiàn)有解。

4.2 計算流程圖

人工蜂群優(yōu)化算法的流程如圖1所示。4.3 多區(qū)域經(jīng)濟調(diào)度問題的實現(xiàn)

Figure 1 Flowchart of the artificial bee colony optimization圖1 人工蜂群優(yōu)化法流程圖

4.3.1 初始化

設Pn=［（P11，P12，…，P1M1），…，（Pi1，Pi2，…，PiMi），…，（PN1，PN 2，…，PNM N），（T12，T13，…，T1N），（T23，T24，…，T2N），…，（T（N—1）N）］為種群進化的第n個向量，n=1，2，…，NP。Pn為所有區(qū)域內(nèi)發(fā)電機的實際功率輸出和聯(lián)絡線實際功率流。在區(qū)域i內(nèi)的第j個實際功率輸出通過設定Pij～U（，）來確定，式中i=1，2，…，N；j=1，2，…，Mi。聯(lián)絡線實際功率流由Tik～U（—，）確定。U（a，b）表示均勻分布的隨機變量的范圍［a，b］。每個向量應滿足式（4）、式（6）～式（8）的約束條件。

4.3.2 適應度函數(shù)

初始種群的適應度函數(shù)如下：

4.3.3 初始種群的選擇

基于最佳適應度的鄰域搜索，選擇m個最優(yōu)解，由此確定每個最優(yōu)解的鄰域搜索的范圍。

4.3.4 鄰域解

在鄰域搜索范圍內(nèi)，圍繞每個被選擇的解，使用式（13）和式（14）產(chǎn)生nb個鄰域解。

其中，mulG和mulT分別表示實際發(fā)電功率和聯(lián)絡線功率傳輸?shù)谋壤蜃?，N（0，1）表示標準正態(tài)分布。每個鄰域解應滿足式（4）、式（6）～式（8）的約束條件。

4.3.5 選擇

利用式（12）對m×nb個解的適應度進行評價，并從中選擇m個最優(yōu)解。

4.3.6 終止

若循環(huán)達到規(guī)定的最大迭代次數(shù)（Nmax），搜索過程結束，當前狀況下的解為最優(yōu)解；否則選擇由鄰域解產(chǎn)生的m個最優(yōu)解。

5 仿真實驗

為了驗證人工蜂群優(yōu)化算法ABCO（Artificial Bee Colony Optimization）的有效性，將相同的兩個測試系統(tǒng)分別用差分進化法DE（Differential E-volution）、進化規(guī)劃法EP（Evolutionary Programming）和實數(shù)編碼遺傳算法RCGA（Real Coded Genetic Algorithm）進行運算并進行結果對比。實驗以MATLAB7.0為仿真環(huán)境，在Intel Core i5 3.4 GHz/8 GB/Windows 7的PC上完成。

5.1 算例1

該算例系統(tǒng)由兩個區(qū)域組成，每個區(qū)域由三臺存在禁止運轉(zhuǎn)區(qū)的發(fā)電機組成，并考慮傳輸損失，其具體數(shù)據(jù)見文獻 ［18］。系統(tǒng)總功率要求是1 263 MW，功率流的限制為100 MW，區(qū)域1和區(qū)域2分別占總功率要求的60%和40%。對于此問題，設置ABCO算法的參數(shù)為：

ns=50，m=30，nb=10，mulG=0.1，mulT= 0.01，Nmax=100。對于該算例系統(tǒng)，分別采用DE、EP和RCGA算法進行對比驗證。設定差分進化

法中種群大小、比例因子和交叉常數(shù)分別為200、1.0和1.0。進化規(guī)劃法中種群大小和比例因子分別選擇為100和0.1；實數(shù)編碼遺傳算法中種群大小、交叉和變異概率分別選取100、0.9和0.2。DE、EP和RCGA三種方法的最大迭代次數(shù)取100次。該測試系統(tǒng)計算結果見表1，成本收斂特性如圖2所示。

Table 1 Simulation results for test system 1表1 測試系統(tǒng)1仿真結果

Figure 2 Cost convergence characteristic of test system 1圖2 測試系統(tǒng)1成本收斂特性圖

5.2 算例2

該算例系統(tǒng)由具有閥點負載的10個發(fā)電機組和3個燃料選項的多燃料源組成，并且考慮傳輸損失，其具體數(shù)據(jù)見文獻［17］。系統(tǒng)總功率要求是2 700 MW。10個發(fā)電機組分成3個區(qū)域，區(qū)域1到區(qū)域3的發(fā)電機組數(shù)量按順序依次是4、3、3，區(qū)域1到區(qū)域3的功率要求分別占總功率要求的50%、25%和25%。各區(qū)域之間的功率流限制為100 MW。采用ABCO算法選取的參數(shù)為ns=50，m=30，nb=10，mulG=0.1，mulT=0.01，Nmax= 300。對于該算例系統(tǒng)，分別采用DE、EP和RCGA算法進行對比驗證。設定差分進化法中種群大小、比例因子和交叉常數(shù)分別為200、1.0和1.0。進化規(guī)劃法中種群大小和比例因子分別選擇為100和0.1；實數(shù)編碼遺傳算法中種群大小、交叉和變異概率分別選取100、0.9和0.2。DE、EP和RCGA三種方法的最大迭代次數(shù)取100次。該測試系統(tǒng)計算結果見表2，成本收斂特性如圖3所示。采用DE、EP和RCGA算法進行對比驗證。設定差分進化法中種群大小、比例因子和交叉常數(shù)分別為400、1.0和1.0。進化規(guī)劃法中種群大小和比例因子分別選擇為200和0.1；實數(shù)編碼遺傳算法中種群大小、交叉和變異概率分別選取200、0.9和0.2。DE、EP和RCGA三種方法的最大迭代次數(shù)取500次。該測試系統(tǒng)計算結果見表4，成本收斂特性如圖4所示。

Table 2 Simulation results for test system 2表2 測試系統(tǒng)2仿真結果

Figure 3 Cost convergence characteristic of test system 2圖3 測試系統(tǒng)2成本收斂特性圖

Table 3 Power flow limit between areas表3 各區(qū)域間的功率流限制

5.3 算例3

該算例系統(tǒng)由具有閥點負載的40個發(fā)電機組組成。各個發(fā)電機參數(shù)值見文獻［14］，系統(tǒng)總功率要求是10 500 MW。40個發(fā)電機組按機組數(shù)量依順序平均分成4個區(qū)域，每個區(qū)域機組數(shù)量是10個，區(qū)域1到區(qū)域4的功率要求分別占總功率要求的15%、40%、30%和15%。各區(qū)域之間的功率流限制見表3。

采用ABCO算法選取的參數(shù)為ns=100，m= 50，nb=20，mulG=0.1，mulT=0.01，Nmax=500。本算例中不考慮傳輸損失。對于該算例系統(tǒng)，分別

采用同樣實驗環(huán)境，分別將本文 ABCO算法與DE、EP、RCGA算法求解此問題的計算結果進行比較。通過對比分析表1、表2和表4中的數(shù)據(jù)可知，ABCO算法對多區(qū)域電力系統(tǒng)經(jīng)濟調(diào)度問題是有效的，而且算法性能也較其他算法更優(yōu)，改善了計算結果，使系統(tǒng)總運行費用有所下降，節(jié)約了燃料的運行成本，減少了運算時間。從圖2、圖3和圖4可見，ABCO算法的收斂速度更快，且算法的優(yōu)化過程更趨于連續(xù)穩(wěn)定。

Table 4 Simulation results for test system 3表4 測試系統(tǒng)3仿真結果

6 結束語

Figure 4 Cost convergence characteristic of test system 3圖4 測試系統(tǒng)3成本收斂特性圖

理論分析和仿真計算表明，本文所提人工蜂群優(yōu)化算法在解決多區(qū)域電力系統(tǒng)經(jīng)濟調(diào)度分配問題時，能在復雜的約束條件下，找到全局最優(yōu)解，具有很強的尋優(yōu)能力，在收斂速度和解的質(zhì)量上都優(yōu)于DE、EP、RCGA算法，驗證了本文算法的可行性和有效性。

［1］ Liu Zhuo.Optimization principle and its application in power system［M］.Harbin：Publishing House of Harbin Institute of Technology，1988.（in Chinese）

［2］ Chowdhury B H，Rahman S.A review of recent advances in economic dispatch［J］.IEEE Transactions on Power Systems，1990，5（4）：1248-1259.

［3］ Shoults R R，Chang S K，Helmick S，et al.A practical approach to unit commitment economic dispatch and savings allocation for multiple-area pool operation with import/export constraints［J］.IEEE Transactions on Power Apparatus and Systems，1980，99（2）：625-635.

［4］ Romano R，Quintana V H，Lopez R，et la.Constrained economic dispatch of multi-area systems using the Dantzig-Wolfe decomposition principle［J］.IEEE Transactions on Power Apparatus and Systems，1981，100（4）：2127-2137.

［5］ Doty K W，McEntire P L.An analysis of electric power brokerage systems［J］.IEEE Transactions on Power Apparatus and Systems，1982，101（2）：389-396.

［6］ Streiffert D.Multi-area economic dispatch with tie line constraints［J］.IEEE Transactions on Power Systems，1995，10 （4）：1946-1951.

［7］ Yalcinoz T，Short M J.Neural networks approach for solving economic dispatch problem with transmission capacity constraints［J］.IEEE Transactions on Power Systems，1998，13 （2）：307-313.

［8］ Jayabarathi T，Sadasivam G，Ramachandran V.Evolutionary programming based multi-area economic dispatch with tie line constraints［J］.Electrical Machines and Power Systems，2000，28：1165-1176.

［9］ Wang Ling-feng，Singh Chan-an.Reserve-constrained multiarea environmental/economic dispatch based on particle swarm optimization with local search［J］.Engineering Applications of Artificial Intelligence，2009，22（2）：298-307.

［10］ Kalpana J，Manjaree P.Discussion of reserve constrained multiarea economic dispatch employing differential evolution with time-varying mutation［J］.International Journal of Electrical Power and Energy Systems，2012，39（1）：68-69.

［11］ He Da-kuo，Wang Fu-li，Mao Zhi-zhong.Economic dispatch of power systems based on an improved genetic algorithm ［J］.Control and Decision，2007，22（2）：230-232.（in Chinese）

［12］ Wang Yu，Li Bin，Yuan Bo.Hybrid of comprehensive learning particle swarm optimization and SQP algorithm for large scale economic load dispatch optimization of power system ［J］.Science China，2010，40（3）：403-411.（in Chinese）

［13］ Manisha S，Manjaree P，Laxmi S.Reserve constrained multi-area economic dispatch employing differential evolution with time-varying mutation［J］.International Journal of Electrical Power and Energy Systems，2011，33（3）：753-766.

［14］ Sinha N，Chakrabarti R，Chattopadhyay P K.Evolutionary programming techniques for economic load dispatch［J］. IEEE Transactions on Evolutionary Computation，2003，7 （1）：83-94.

［15］ Karaboga D，Basturk B.A powerful and efficient algorithm for numerical function optimization：Artificial bee colony （ABC）algorithm［J］.Journal of Global Optimization，2007，39（3）：459-471.

［16］ Walter D C，Sheble G B.Genetic algorithm solution of economic dispatch with valve point loading［J］.IEEE Transactions on Power Systems 1993，8（3）：1325-1332.

［17］ Chiang C-L.Improved genetic algorithm for power economic dispatch of units with valve-point effects and multiple fuels ［J］.IEEE Transactions on Power Systems，2005，20（4）：1690-1699.

［18］ Gaing Z-L.Particle swarm optimization to solving the economic dispatch considering the generator constraints［J］. IEEE Transactions on Power Systems，2003，18（3）：1187-1195.

附中文參考文獻：

［1］ 柳焯.最優(yōu)化原理及其在電力系統(tǒng)中的應用［M］.哈爾濱：哈爾濱工業(yè)大學出版社，1988.

［11］ 何大闊，王福利，毛志忠.基于改進遺傳算法的電力系統(tǒng)經(jīng)濟負荷分配［J］.控制與決策，2007，22（2）：230-232.

［12］ 王瑜，李斌，袁博.混合SQP的基于完全學習的粒子群優(yōu)化算法在電力系統(tǒng)中經(jīng) 濟分配問題 的應用［J］.中國科學，2010，40（3）：403-411.

鄭曉菁（1978），女，安徽滁州人，碩士，實驗師，研究方向為電氣信息監(jiān)測和計算機仿真。E-mail：liht@aust.edu.cn

ZHENG Xiao-jing，born in 1978，MS，experimentalist，her research interests include electric information monitoring，and simulation.

Multi-area economic dispatch of power system based on artificial bee colony optimization

ZHENG Xiao-jing
（School of Electrical and Information Engineering，Anhui University of Science&Technology，Huainan 232001，China）

Aiming at the problems of multi-area economic dispatch（MAED）of the power system，such as tie line transmission losses，multiple fuels，valve-point loading and prohibited operating zones，we design a mathematical model in which the requirement of the minimum cost of multi-area power load is taken into account and the artificial bee colony optimization is utilized to quickly search for the global optimal solution.The effectiveness and feasibility of the proposed algorithm have been verified on two different test systems，both small and large，involving varying degrees of complexity.Compared with algorithms including differential evolution，evolutionary programming and real coded genetic algorithm，the proposed algorithm is a promising alternative approach for solving the MAED problems in practical power system.

multi-area economic dispatch；artificial bee colony optimization；tie line constraints；power system；power balance constraints

TM721

A

10.3969/j.issn.1007-130X.2015.08.017

1007-130X（2015）08-1533-07

2014-10-14；

2014-12-16

通信地址：232001安徽省淮南市安徽理工大學電氣與信息工程學院

Address：School of Electrical and Information Engineering，Anhui University of Science&Technology，Huainan 232001，Anhui，P.R.

China