基于細(xì)菌菌落算法的含分布式電源多目標(biāo)無功優(yōu)化

黃沂平, 馬雷鵬, 張永超, 葉 莘

(1.長樂市電力公司, 福建 長樂 350200; 2.上海電力學(xué)院, 上海 200090)

基于細(xì)菌菌落算法的含分布式電源多目標(biāo)無功優(yōu)化

黃沂平1, 馬雷鵬2, 張永超2, 葉 莘2

(1.長樂市電力公司, 福建 長樂 350200; 2.上海電力學(xué)院, 上海 200090)

以含風(fēng)電和光伏的配電網(wǎng)為研究對象,提出了含風(fēng)電和光伏的配電網(wǎng)多目標(biāo)無功優(yōu)化模型.用場景分析法對雙饋異步風(fēng)力發(fā)電機(jī)(DFIG)進(jìn)行了場景分析,用Beta分布函數(shù)來模擬光伏出力,通過蒙特卡洛方法對配電網(wǎng)進(jìn)行了無功補(bǔ)償選址,最后通過細(xì)菌菌落優(yōu)化算法對IEEE33節(jié)點配電系統(tǒng)進(jìn)行了算例分析,驗證了該模型及方法的有效性.

多目標(biāo)無功優(yōu)化; 分布式電源; 細(xì)菌菌落算法; 蒙特卡洛法

分布式發(fā)電(Distributed Generation,DG)具有靈活、高效、環(huán)保、削峰填谷等優(yōu)勢,同時能提高供電質(zhì)量,發(fā)展DG意義尤為重大[1].在清潔替代和電能替代的要求下,風(fēng)電和光伏的發(fā)展成為必然的選擇,但是風(fēng)電和光伏出力的隨機(jī)性及不確定性給配電網(wǎng)無功優(yōu)化帶來了困難.

文獻(xiàn)[2]和文獻(xiàn)[3]探討了不同出力模式下,分布式電源對配電網(wǎng)的電壓、規(guī)劃、損耗、供電可靠性等方面帶來的影響;文獻(xiàn)[4]和文獻(xiàn)[5]探討了含有風(fēng)電的配電網(wǎng)無功優(yōu)化,用遺傳算法控制無功補(bǔ)償?shù)娜萘考巴肚须娙萜鞯慕M數(shù);文獻(xiàn)[6]計及光伏電站并網(wǎng)逆變器的無功調(diào)節(jié)能力,采用改進(jìn)粒子群算法對規(guī)劃模型進(jìn)行優(yōu)化;文獻(xiàn)[7]采用基本螢火蟲算法驗證分布式電源接入能夠提高電壓穩(wěn)定性,然而基本螢火蟲算法易陷入局部最優(yōu)解,且算法收斂速度不快;文獻(xiàn)[8]至文獻(xiàn)[10]根據(jù)細(xì)菌菌落算法的特性,驗證了該算法具有較強(qiáng)的全局尋優(yōu)能力,且收斂速度快,魯棒性好,可以作為求解配電網(wǎng)優(yōu)化問題的一種新途徑.本文首先對雙饋異步風(fēng)力發(fā)電機(jī)組(Double-Fed Induction Generator,DFIG)的出力進(jìn)行了場景分析,并用Beta函數(shù)對光伏進(jìn)行了模擬,建立包括網(wǎng)絡(luò)損耗、節(jié)點電壓偏移及靜態(tài)穩(wěn)定性在內(nèi)的多目標(biāo)無功優(yōu)化模型,然后采用細(xì)菌菌落算法對該函數(shù)模型進(jìn)行了優(yōu)化.

1 DG模型

1.1 基于場景分析的DFIG模型

DFIG出力的隨機(jī)性使得配電網(wǎng)無功優(yōu)化具有不確定性,因此采用場景理論對其輸出功率進(jìn)行場景分析,可將DFIG出力的隨機(jī)性轉(zhuǎn)化為多個階段內(nèi)的確定性來進(jìn)行求解[9].可從DFIG的功率特性曲線上得出在某一風(fēng)速下的有功功率Pw.圖1為1.5 MW DFIG的功率特性曲線.

圖1 DFIG的功率特性曲線

由圖1可知,DFIG有功輸出功率Pw與風(fēng)速v的函數(shù)關(guān)系式為:

(1)

式中:Pr——DFIG的額定輸出功率;vc,vr,vo——切入風(fēng)速、額定風(fēng)速和切出風(fēng)速.

由圖1可知,DFIG有3種運(yùn)行狀態(tài),分別與3種場景相對應(yīng),即:停機(jī)狀態(tài)-零輸出場景;欠額定運(yùn)行狀態(tài)-欠額定場景;額定運(yùn)行狀態(tài)-額定場景.零輸出和額定場景下的場景功率取0和Pr,欠額定場景下的場景功率Pw2為:

(2)

本文采用Weibull模型來模擬風(fēng)速情況,其概率的密度及分布函數(shù)為:

(3)

式中:v——DFIG的風(fēng)速;k,c——Weibull兩參數(shù)模型中的形狀和尺度參數(shù),分別取1.93和9.19[11].

由此可以求出每一種場景下的功率.因此,本文選用了一臺風(fēng)機(jī),在實際工程應(yīng)用中,情況將更為復(fù)雜,場景可以根據(jù)風(fēng)電場的規(guī)模劃分得更細(xì).

1.2 基于Beta分布的光伏發(fā)電模型

光伏發(fā)電的有功功率輸出主要隨光照強(qiáng)度的變化而變化.據(jù)分析,光照強(qiáng)度服從Beta分布[12-13],因此可以用Beta分布來描述光伏發(fā)電的功率.其概率的密度函數(shù)為:

(4)

(5)

式中:PM,Pmax——有功功率輸出和最大有功功率輸出;

α,β——Beta分布模型的形狀參數(shù);

μ,δ——一段時間內(nèi)光照強(qiáng)度的平均值和方差.

光伏電站會同時向配電網(wǎng)提供有功功率和無功功率,說明其也參與配電網(wǎng)的無功優(yōu)化.提供的無功功率為:

(6)

式中:Smax,Pa——光伏并網(wǎng)逆變器提供的最大視在功率和最大有功功率.

2 配電網(wǎng)無功優(yōu)化模型

2.1 目標(biāo)函數(shù)

本文所建立的無功優(yōu)化模型同時考慮了配電網(wǎng)的經(jīng)濟(jì)性、可靠性及穩(wěn)定性.所建立的多目標(biāo)無功優(yōu)化模型函數(shù)為:

(7)

式中:F1——配電網(wǎng)網(wǎng)損;F2,F3——節(jié)點電壓偏差及穩(wěn)定性指標(biāo);Na,Nb——配電網(wǎng)的支路數(shù)和節(jié)點數(shù);Gij——節(jié)點i與j之間的電導(dǎo);Ui,Uj——節(jié)點i和j的電壓幅值;δij——節(jié)點i和j之間的電壓相角差;UiN,ΔUi,max——節(jié)點i的電壓期望值和最大電壓偏差值;

Usta.i——節(jié)點i的電壓穩(wěn)定性指標(biāo)值;

Rij,Xij——節(jié)點i與j之間的電阻和電抗;

Pj,Qj——節(jié)點i與j注入的有功和無功功率.

由于目標(biāo)函數(shù)的量綱不同,因此對上述目標(biāo)函數(shù)作歸一化處理,可得:

(8)

最終目標(biāo)函數(shù)為:

(9)

式中:F10,F20,F30——各指標(biāo)的最優(yōu)值;F12,F22,F32——各指標(biāo)優(yōu)化前的初始值;λ1,λ2,λ3——最終目標(biāo)函數(shù)對網(wǎng)損、節(jié)點電壓偏差、電壓穩(wěn)定性指標(biāo)的偏好程度,滿足λ1+λ2+λ3=1且λ1,λ2,λ3均大于零.

2.2 約束條件

2.2.1 等式約束條件

負(fù)荷節(jié)點必須滿足有功功率及無功功率的平衡,即:

(10)

式中:Pi,PGi,PLi——節(jié)點i注入的有功功率、發(fā)電機(jī)節(jié)點的有功出力和負(fù)荷節(jié)點消耗的有功功率;

Qi,QGi,QCi,QLi——節(jié)點i注入的無功功率、發(fā)電機(jī)節(jié)點的無功出力、無功補(bǔ)償容量和負(fù)荷節(jié)點消耗的無功功率;

Gij,Bij,θij——節(jié)點i與j之間的電導(dǎo)、電納和相角差.

2.2.2 不等式約束條件

(1) 控制變量不等式約束條件為:

(11)

式中:Ti——變壓器第i個分接頭的位置;QCi——第i組并聯(lián)電容器的補(bǔ)償容量;QGi——分布式電源發(fā)出的無功功率;

Ti,max,Ti,min——變壓器變比的上下限;

QCi,max,QCi,min——并聯(lián)電容器補(bǔ)償容量的上下限;

QGi,max,QGi,min——分布式電源無功功率的上下限.

(2) 狀態(tài)變量不等式約束條件為:

(12)

式中:Ui,max,Ui,min——節(jié)點i的電壓上下限值.

3 無功優(yōu)化模型的求解算法及算例分析

3.1 基于蒙特卡洛仿真的無功補(bǔ)償選址

配電網(wǎng)未來運(yùn)行方式是不確定的,因此對其在特定方式下進(jìn)行無功補(bǔ)償存在一定的局限性.基于蒙特卡洛仿真的無功補(bǔ)償選址充分考慮了各種不確定性因素,可對配電網(wǎng)運(yùn)行進(jìn)行隨機(jī)性模擬[14].其基本思路為:對含分布式電源的配電網(wǎng),建立各種運(yùn)行環(huán)境下的概率模型,對各種電氣元件進(jìn)行仿真建模分析,從而確定系統(tǒng)的運(yùn)行狀態(tài),并確定該運(yùn)行狀態(tài)下的潮流分布情況.若某些負(fù)荷節(jié)點出現(xiàn)電壓越下限等情況,則對這些節(jié)點進(jìn)行無功補(bǔ)償,其統(tǒng)計指標(biāo)如下[15].

(1) 節(jié)點電壓平均值為:

(2) 節(jié)點無功補(bǔ)償概率為:

基于蒙特卡洛仿真的無功補(bǔ)償選址流程如圖2所示.

圖2 基于蒙特卡洛仿真的無功補(bǔ)償選址

3.2 基于細(xì)菌菌落優(yōu)化算法的無功優(yōu)化

細(xì)菌菌落算法通過模擬單個及群體細(xì)菌的生長過程來尋找最優(yōu)解[16].其規(guī)定個體的生長繁殖不能超出種群的規(guī)模N,而且個體在滿足規(guī)定的死亡條件時死亡.細(xì)菌培養(yǎng)液及培養(yǎng)液中營養(yǎng)物質(zhì)的濃度分別表示待優(yōu)化模型的解空間及單個目標(biāo)函數(shù)的優(yōu)化程度,培養(yǎng)液濃度和外界條件可以制約細(xì)菌個體的生長繁殖.

細(xì)菌個體性能的優(yōu)劣程度受待優(yōu)化模型目標(biāo)函數(shù)適應(yīng)值的影響.在每一次迭代計算后,細(xì)菌不停地進(jìn)行更新來改變其位置及速度,直至找到最優(yōu)點.細(xì)菌個體具有游走和翻滾兩個屬性.游走即沿著上次的運(yùn)動方向向有利的環(huán)境方向繼續(xù)運(yùn)動;翻滾就是換個方向運(yùn)動,避開不利環(huán)境.經(jīng)過不停的游走和翻滾,可以實現(xiàn)目標(biāo)函數(shù)的優(yōu)化.細(xì)菌個體游走的更新公式為:

(13)

式中:Vk——第k次迭代時細(xì)菌的前進(jìn)方向;

Xk——第k次迭代時細(xì)菌的位置;

fbest——上一次迭代后細(xì)菌的最優(yōu)位置;

gbest——目前菌落全局的最優(yōu)位置;

α,r1,r2——系數(shù);

R——0與1之間的隨機(jī)數(shù).

當(dāng)個體目標(biāo)模型適應(yīng)值沒有上一次優(yōu)越時,表示細(xì)菌個體沒有到達(dá)營養(yǎng)物濃度高的點.這時細(xì)菌個體會通過翻滾的運(yùn)動方式來改變運(yùn)動的方向及位置,其更新公式為:

(14)

式中:r——搜索半徑;Rn——-1與1之間的隨機(jī)數(shù).

細(xì)菌菌落優(yōu)化流程如圖3所示.

圖3 細(xì)菌菌落優(yōu)化流程

3.3 算例分析

采用含風(fēng)電和光伏的IEEE33節(jié)點配電系統(tǒng)作為算例來驗證本文所提的多目標(biāo)無功優(yōu)化模型及算法的適用性.設(shè)系統(tǒng)的基準(zhǔn)電壓值為12.66 kV,基準(zhǔn)功率值為10 MW,有載調(diào)壓變壓器的調(diào)節(jié)范圍為0.95～1.05,共9個檔位,步進(jìn)量為1.25%.在節(jié)點3處并入光伏電站,用PV表示,其服從Beta分布函數(shù)的形狀參數(shù)α=β=0.85.在節(jié)點18及33處接入DFIG,用W1,W2表示.

DFIG的參數(shù):額定功率為1.5 MW,vc,vr,vo分別取3 m/s,12 m/s,25 m/s;場景1,場景2,場景3發(fā)生的概率為0.109 9,0.703 5,0.186 6;由式(3)計算出各場景功率分別為0 MW,0.283 9 MW,1.5 MW,則可以求出DFIG的期望輸出功率為0.479 6 MW.細(xì)菌種群規(guī)模N=20,生長周期S=4,繁殖和死亡條件Np=2,系數(shù)α=2,r1=r2=1.5,搜索半徑R=0.1.

改進(jìn)的IEEE33節(jié)點配電系統(tǒng)如圖4所示.

圖4 改進(jìn)的IEEE33節(jié)點配電系統(tǒng)

采用蒙特卡洛算法對改進(jìn)的系統(tǒng)進(jìn)行抽樣,所得出的系統(tǒng)節(jié)點的補(bǔ)償次數(shù)和補(bǔ)償概率見表1.

表1 含DG的IEEE 33節(jié)點配電系統(tǒng)無功補(bǔ)償結(jié)果

由表1可知,此系統(tǒng)中節(jié)點30和14電壓損耗最大,需要進(jìn)行無功補(bǔ)償?shù)拇螖?shù)和概率也最大,即在這兩節(jié)點處進(jìn)行補(bǔ)償能達(dá)到最佳效果.因此,在節(jié)點14和30處分別并入4組和7組電容器,每組并聯(lián)電容器的補(bǔ)償容量為100 kW.

對改進(jìn)的配電系統(tǒng)及加入DG后的配電系統(tǒng)的無功優(yōu)化結(jié)果進(jìn)行對比,結(jié)果如表2所示.

表2 改進(jìn)的IEEE 33節(jié)點配電系統(tǒng)優(yōu)化結(jié)果對比 (p.u.)

由表2可以看出,在加入風(fēng)電機(jī)組和光伏后,通過DG的無功調(diào)節(jié)能力,系統(tǒng)的網(wǎng)損、電壓偏差、電壓穩(wěn)定指標(biāo)都得到了有效改善.

在配電系統(tǒng)經(jīng)過優(yōu)化之后,可以得到總的目標(biāo)函數(shù)值F,進(jìn)而可以得到各個分目標(biāo)函數(shù)值,即網(wǎng)損Ploss,電壓偏差ΔU,電壓穩(wěn)定性指標(biāo)Usta,因此可以得到各節(jié)點電壓值Ui.優(yōu)化前后各節(jié)點電壓值如圖5所示.由圖5可以看出,本文的無功優(yōu)化方法能有效降低配電網(wǎng)的網(wǎng)損、保持節(jié)點電壓的穩(wěn)定性等.

圖5 優(yōu)化前后的各節(jié)點電壓

4 結(jié) 論

(1) 采用場景分析法和Beta函數(shù)有效解決了風(fēng)機(jī)及光伏的隨機(jī)性問題.

(2) 采用蒙特卡羅法對不同運(yùn)行方式下的配電網(wǎng)進(jìn)行無功補(bǔ)償選址,得出節(jié)點的無功補(bǔ)償次數(shù)及概率;用細(xì)菌菌落算法對模型進(jìn)行優(yōu)化,并用算例驗證二者具有優(yōu)越性,由此表明本文所建立的多目標(biāo)無功優(yōu)化模型具有較好的適用性.

[1] 呂忠,周強(qiáng),蔡雨昌.含分布式電源的DEIWO算法配電網(wǎng)無功優(yōu)化[J].電力系統(tǒng)保護(hù)與控制,2015(4):69-73.

[2] 李斌,劉天琪,李興源.分布式電源接入對系統(tǒng)電壓穩(wěn)定性的影響[J].電網(wǎng)技術(shù),2009,33(3):84-88.

[3] 韋鋼,吳偉力,胡丹云,等.分布式電源及其并網(wǎng)時對電網(wǎng)的影響[J].高電壓技術(shù),2007,33(1):36-40.

[4] 陳樹勇,申洪,張洋,等.基于遺傳算法的風(fēng)電場無功補(bǔ)償及控制方法的研究[J].中國電機(jī)工程學(xué)報,2005,25(8):1-6.

[5] 劉學(xué)平,劉天琪,李興源.含風(fēng)電機(jī)組的配電網(wǎng)無功優(yōu)化補(bǔ)償[J].電力系統(tǒng)保護(hù)與控制,2010,38(20):130-135.

[6] 孫卓新,朱永強(qiáng),倪一峰,等.基于粒子群算法的含光伏電站的配電網(wǎng)無功優(yōu)化[J].電力建設(shè),2014,35(4):25-30.

[7] HERBADJI O,NADHIR K,SLIMANI L,etal.Optimal power flow with emission controlled using firefly algorithm[C].International Conference on Modeling,IEEE,2013:1-6.

[8] ABDI Sh,AFSHAR K.Application of IPSO—monte carlo for optimal distributed generation allocation and sizing[J].Int J. Elect. Power Energy System,2013,44(1):786-797.

[9] 任新偉,徐建政,趙斌.含光伏電站的配電網(wǎng)無功優(yōu)化[J].電力電容器與無功補(bǔ)償,2014,35(1):12-15.

[10] 簡獻(xiàn)忠,李瑩,周海,等.基于細(xì)菌菌落算法的電力系統(tǒng)無功優(yōu)化[J].控制工程,2014,21(6):935-938.

[11] 李闖,陳民鈾,付昂,等,含風(fēng)電場的配電網(wǎng)無功優(yōu)化策略研究[J].電力系統(tǒng)保護(hù)與控制,2013(9):100-105.

[12] ABOUZAHR I,RAMAKUMAR R.Loss of power supply probability of stand-alone photovoltaic systems:a closed form solution approach[J].IEEE Transactions on Energy Conversion,1991,6(1):1-11.

[13] KARAKI S H,CHEDID R B,RAMADAN R.Probabilistic performance assessment of autonomous solar-wind energy conversion systems[J].IEEE Transactions on Energy Conversion,1999,14(3):766-772.

[14] 張瑞華,宋云亭.基于蒙特卡羅仿真和電壓安全約束的無功優(yōu)化算法[J].電力系統(tǒng)自動化,2002,26(7):23-27.

[15] 黃澤霞,俞攸紅,黃德才.慣性權(quán)自適應(yīng)調(diào)整的量子粒子群優(yōu)化算法[J].上海交通大學(xué)學(xué)報,2012,46(2):228-232.

[16] 李明,楊成梧.細(xì)菌菌落優(yōu)化算法[J].控制理論與應(yīng)用,2011,28(2):223-22.

(編輯 胡小萍)

Multi-objective Reactive Power Optimization with DistributedGeneration Based on Bacterial Colony Optimization Algorithm

HUANG Yiping1, MA Leipeng2, ZHANG Yongchao2, YE Shen2

(1.ChanglePowerSupplyCorporation,Changle350200,China; 2.ShanghaiUniversityofElectricPower,Shanghai200090,China)

Distribution network with wind power and photovoltaic is studied,and the multi-objective reactive power optimization model of the distribution network is proposed.The scene analysis method is used to analyze the scene of Double Fed Induction Generator(DFIG) and the Beta distribution function is used to simulate the PV output,the location of reactive power compensation device for distribution network is determined by Monte Carlo simulation.Finally,the IEEE 33 node system is used as a case by the bacterial colony optimization algorithm,verifying the validity of the model and method.

multi-objective reactive power optimization; distributed generation; bacterial colony optimization algorithm; Monte Carlo simulation

10.3969/j.issn.1006-4729.2016.06.007

2016-03-24

簡介：馬雷鵬(1989-),男,在讀碩士,山東昌邑人.主要研究方向為配電網(wǎng)無功優(yōu)化.E-mail:maleipeng1990@126.com.

TM714.3

A

1006-4729(2016)06-0538-05