向 山,向迪輝
(張家界電業局,湖南 張家界 427000)
基于量子粒子群混合算法的電網經濟調度
向 山,向迪輝
(張家界電業局,湖南 張家界 427000)
電網經濟調度 (ED)是電力系統中一個重要的經濟問題。針對ED問題目標函數非線性、不連續的特點,提出了一種量子粒子群混合算法 (QPSO)。該算法用量子疊加態表示單個粒子,增加了粒子的多樣性;使用量子旋轉門更新粒子,提高了粒子的收斂速度。通過實際算例對算法進行仿真驗證,結果表明,QPSO算法能夠有效解決ED問題,對降低發電成本、減少燃料消耗有一定實用價值。
電網;經濟調度;粒子群算法;量子疊加態;QPSO
電網經濟調度 (ED)[1]問題的目標是在滿足系統運行以及負荷約束的條件下,給各臺發電機合理分配負荷,使得發電機組總的發電成本最低。ED問題是電力系統運行中必須面對的一類優化問題。合理解決ED問題能降低發電成本,具有重大的經濟意義,是實現電網生產、管理、運行最優的基本保障。
實際電力系統中建立的ED數學模型呈現不連續、非線性、不可微、不可導的特點,使得像拉格朗日和等微增率法這類傳統的算法,在處理ED問題時不能取得理想的結果。在這個背景下,一系列人工智能算法被用來解決ED問題,如遺傳算法[2]、粒子群優化算法 (PSO)[3]、混沌優化算法[4]。其中,PSO算法因為具有容易理解和實現的特點,被廣大學者進行研究,并取得了一定成果。
PSO算法的思想起源是鳥類覓食的過程,通過種群團體合作來尋找最優解,具有計算速度快、流程簡單易實現的特點。因為PSO算法存在后期收斂速度變慢和容易陷入局部最優的問題,所以又出現了一些改進的粒子群算法。本文所采用的量子粒子群混合算法 (QPSO)[5],用量子疊加態表示單個粒子,增加了粒子的多樣性;使用量子旋轉門更新粒子,提高了粒子的收斂速度。

式中:F表示發電機總的發電費用;Fi(Pi)表示第i臺發電機的耗量特性;Ei為閥點效應引起的耗量特性變化;Ng為發電機總臺數;αi、βi、γi、hi、gi表示第i臺發電機的耗量特性曲線系數;Pi表示第i臺發電機的有功功率;Pi,min表示第i臺發電機的有功功率下限。
a. 發電機運行約束

式中:Pi為第 i臺發電機的出力;Pi,min、Pi,max分別表示第i臺發電機的有功功率下限和上限。
b. 功率平衡約束

式中:PL為系統總負荷;PLoss為系統總網損。網損一般采用潮流計算或B系數法求得,B系數法網損可按照下式計算:

式中:P為Ng維發電機有功功率列向量;PT為P的轉置;B為Ng×Ng維對稱方陣;B0為Ng維列向量;B00為常數。
設初始種群具有N個粒子,每個粒子維數是m,每個粒子都用一個相位θ來描述,即每個粒子都可以按下面的公式編碼:

式中:j=1,2,…,N;d=1,2,…,m;Pjd(α)、Pjd(β)分別表示第j個粒子第d維處于0態和1態的概率分別為|α|2、|β|2;θjd表示j個粒子第d維的量子相位。
設量子粒子群算法的實際搜索空間是[a,b],而粒子的編碼概率都在[0-1]內,因而需要將編碼概率映射到實際的參數空間,空間變換公式如下:

式中:P為狀態表達的選擇概率;R為 (0,1)之間的隨機數;Pjd為第j個粒子第d維的實際位置。
本文采用的算法采用標準粒子群算法更新公式來更新量子相位,更新公式如下:

式中:w為慣性權重;c1、c2是大于0的學習因子(一般取1.5~2.5);r1、r2為0到1 之間均勻分布的隨機數,t是迭代次數;Δθt+1jd、Δθtjd分別為第t+1、t次迭代中第j個粒子第d維的相位修正量;θjd表示j個粒子第d維的量子相位;θjdb為第j個粒子歷史最優值第d維的相位;θdg為全局最優粒子第d維的相位。
在量子粒子群混合算法中,粒子是以疊加態存在的,因此,需要通過量子旋轉門把相位轉換為以疊加態表示的粒子,轉化公式為

式中:左邊和右邊的列向量分別表示第j個粒子第d維經過t+1、t迭代后的概率幅。
步驟1,隨機產生N個m維的以量子疊加態表示的粒子 (包括它們初始位置和自身最優位置)。
步驟2,通過式 (8)進行解空間變換,并通過目標函數計算適應度,粒子的開始局部最優位置設為它們的初始位置,整個粒子群中適應度最高的粒子位置設為全局最優位置。
步驟3,由式 (9)計算量子相位修正量 (相當于更新速度),再通過式 (10)更新粒子的位置。
步驟4,進行解空間變換,計算適應度,并統計新的局部最優位置和新的全局最優位置。檢查是否達到優化條件,如果達到誤差精度或迭代次數,則結束,否則重復步驟3。
為了驗證QPSO算法在解決ED問題時的有效性,對文獻[6]中的3機和13機系統分別進行仿真測試,仿真中,考慮發電機的閥點效應,忽略網損。最后,本文的仿真結果將與文獻 [6]中采用的三種改進算法 (分別為FEP、MFEP、IFEP)以及文獻 [7]的IGA遺傳算法結果進行比較。發電機的相關參數見文獻 [6]。
算例1:本例采用的是3機6母線系統,系統總負荷為850 MW。仿真中,粒子數d取50,空間維數m取3,學習因子c1、c2都取為1.5,選擇概率p取為0.95,慣性權重w取為0.7,迭代次數為50。發電費用的收斂結果如圖1所示。表1為本算例各種改進算法的結果比較。


表1 算例1各種算法結果的比較
由圖1可知,本文所采用的QPSO算法具有快速收斂的特性,由表1可知,QPSO算法比其他算法的總花費低。可見QPSO算法有一定的優越性。
算例2:本例采用13機組系統,系統承擔總負荷為1 800 MW,空間維數取m為13,迭代次數取200,其他參數和算例1中一致。發電費用的收斂結果如圖2所示,各發電機出力仿真結果如表2所示,表3為本算例各種改進算法的結果比較。


表2 算例2 QPSO發電機出力仿真結果

表3 算例2各種算法結果的比較
由圖2可知,在維數較高情況下,QPSO算法依然能夠較快收斂,由表3數據可知,QPSO依然比其他算法的總花費低,更好地獲得了最優解。
本文首次將量子粒子群混合算法應用到電網經濟調度問題中來,這種算法不僅收斂速度快,同時能夠更好地獲得最優解,克服了粒子群算法容易陷入局部最優的缺陷,能夠有效解決經濟調度的問題,具有一定的實用價值。
[1] 黎 輝.多元投資環境下電網經濟調度研究 [J].東北電力技術,1998,19(4):28-31.
[2] 邢維健,張國立.改進遺傳算法在有功經濟調度中的應用[J].東北電力技術,2004,25(9):10-12.
[3] 白江斌,金慰剛,張建華.基于粒子群的模糊神經網絡[J].東北電力技術,2007,28(3):16-19.
[4] 唐 巍,李殿璞.電力系統負荷經濟分配的混沌優化算法[J].中國電機工程學報,2000,20(10):36-40.
[5] 鄭玲峰,林 輝,樂尚利,等.基于量子粒子群混合算法的電力系統無功優化[J].華中電力,2011,24(2):16-19.
[6] Sinha N,Chakrabarti R,Chattopadhyay PK.Evolutionary Programming techniques for economic load dispatch [J].IEEE Transaction on Evolutionary computation,2003,7(1):83-94.
[7] Ling S H,Lam H K,Leung F H F,et al.Improved genetic algorithm for economic load dispatch with valve-point loadings[C]/The 29 Annual conference of the IEEE Industrial Electronics Society,Piscataway.NJ,2003,1:442-447.
Economic Dispatch Based On Quantum Particle Swarm Optimization In Power System
XIANG Shan,XIANG Di-hui
(Zhangjiajie Electric Power Bureau,Zhangjiajie,Hunan 427000,China)
Economic dispatch(ED)is an important economic problem in power system,as well as an environmental issue.This paper proposes a quantum particle swarm optimization(QPSO)algorithm,for ED problem has the features of nonlinear and discontinuous.The proposed algorithm is demonstrated on practical examples,which reveal that the algorithm can effectively solve the ED problem.
Grid;Economic dispatch;Particle swarm optimization;Quantum superposition state;QPSO
TM73
A
1004-7913(2014)01-0005-03
向 山 (1988—),女,學士,助理工程師,從事電力系統計量工作。
2013-09-30)