張 貞 張 穎 李安均 熊腡成

（智能電網(wǎng)運(yùn)行與控制湖南省重點(diǎn)實(shí)驗(yàn)室（長(zhǎng)沙理工大學(xué)），長(zhǎng)沙 410114）

微電網(wǎng)由分布式電源和電力負(fù)荷構(gòu)成，具有高度的可靠性和穩(wěn)定性[1-2]。隨著微電網(wǎng)控制技術(shù)的不斷發(fā)展和風(fēng)能、太陽能等可再生能源發(fā)電成本的持續(xù)降低，微電網(wǎng)系統(tǒng)可以實(shí)現(xiàn)延緩電網(wǎng)投資、降低網(wǎng)損、提高電能質(zhì)量、節(jié)能減排[3-6]。

目前，國(guó)內(nèi)外一些學(xué)者對(duì)微電網(wǎng)規(guī)劃的理論與實(shí)踐研究，取得了一定的學(xué)術(shù)成果。文獻(xiàn)[7]針對(duì)微電源組合方式多樣、運(yùn)行方式靈活等特點(diǎn)，將微電網(wǎng)評(píng)價(jià)指標(biāo)分為可靠性、經(jīng)濟(jì)性、市場(chǎng)運(yùn)營(yíng)和環(huán)保性4個(gè)方面，對(duì)各指標(biāo)進(jìn)行展開說明，并構(gòu)建了樹狀結(jié)構(gòu)；文獻(xiàn)[8]提出一種含多種可再生能源的分布式發(fā)電系統(tǒng)成本最優(yōu)組合模型；文獻(xiàn)[9]提出微電網(wǎng)發(fā)電備用配置優(yōu)化問題，并從模型建立與求解兩個(gè)方面對(duì)該問題進(jìn)行較為全面的評(píng)述。但是，文獻(xiàn)[7-9]均未考慮分布式電源的輸出功率受風(fēng)能和太陽能隨機(jī)性和波動(dòng)性的影響，加之儲(chǔ)能裝置充放電具有隨機(jī)性，因此，還不能有效、充分、合理地利用微電網(wǎng)來解決一些供電負(fù)荷。

本文對(duì)風(fēng)能、太陽能以及儲(chǔ)能蓄電池中存儲(chǔ)能量的隨機(jī)性進(jìn)行了分析，建立了微電網(wǎng)隨機(jī)優(yōu)化的數(shù)學(xué)模型，并利用SAA近似及粒子群算法對(duì)微電網(wǎng)隨機(jī)優(yōu)化模型進(jìn)行求解，實(shí)現(xiàn)微電網(wǎng)的最優(yōu)經(jīng)濟(jì)調(diào)度分配。算例分析證明了該隨機(jī)優(yōu)化模型的合理性和有效性。

1 風(fēng)光儲(chǔ)混合供電的微電網(wǎng)隨機(jī)性分析

本文研究的微電網(wǎng)系統(tǒng)由風(fēng)電機(jī)組、光伏發(fā)電機(jī)組、儲(chǔ)能裝置等分布式電源及中央控制系統(tǒng)組成，由公共連接點(diǎn)與外網(wǎng)相連。其結(jié)構(gòu)如圖1所示。

圖1 微電網(wǎng)系統(tǒng)結(jié)構(gòu)模型

圖1中，風(fēng)電機(jī)組和光伏發(fā)電機(jī)組保持輸出功率為最大；通信及中央控制系統(tǒng)，實(shí)現(xiàn)對(duì)系統(tǒng)有效管理和控制，并且通過其與各個(gè)系統(tǒng)的相連來控制供電系統(tǒng)的功率輸出。

1.1 風(fēng)力發(fā)電系統(tǒng)的隨機(jī)特性

風(fēng)電場(chǎng)風(fēng)速服從兩參數(shù)Weibull分布[10]，其概率密度函數(shù)為

式中，v為風(fēng)速，k為形狀參數(shù)（無量綱）；c為尺度參數(shù)（m/s）。風(fēng)速的大小決定風(fēng)電機(jī)組機(jī)出力，在不考慮風(fēng)電機(jī)組中發(fā)電機(jī)與傳動(dòng)裝置的特性以及風(fēng)機(jī)間相互聯(lián)系的前提下，以此假設(shè)多臺(tái)風(fēng)機(jī)的風(fēng)向和風(fēng)速相同，可以用一臺(tái)等效的風(fēng)機(jī)代替整個(gè)風(fēng)電機(jī)組[11]，則風(fēng)速-風(fēng)機(jī)出力的函數(shù)關(guān)系表示如下[12]

式中，p( v)表示風(fēng)電機(jī)組的出力；vin為切入風(fēng)速，vrate為額定風(fēng)速，vout為切出風(fēng)速，prate為額定功率，a, b為風(fēng)力機(jī)的功率特性參數(shù)可表示為

1.2 光伏發(fā)電系統(tǒng)的隨機(jī)特性

根據(jù)實(shí)測(cè)經(jīng)驗(yàn)，考慮干擾因素（如陰影、云層等）時(shí)，一天內(nèi)光電機(jī)組接收的光照強(qiáng)度服從Beta分布[13]，其概率密度函數(shù)為

Beta分布是連續(xù)分布，取值在0～1之間。其中α，β均為取值為正數(shù)的形狀參數(shù)，其取值時(shí)要因地制宜；rmax(w/m2)是該段時(shí)間內(nèi)的最大光照強(qiáng)度。

光伏電池的輸出功率在環(huán)境溫度、負(fù)載阻抗等參數(shù)保持不變的情況下與光照強(qiáng)度成正比，則光電機(jī)組的總輸出功率為

式中，η為光電機(jī)組的總效率，A為光伏陣列的總面積，S為光照強(qiáng)度，w/m2。

1.3 儲(chǔ)能蓄電池的隨機(jī)特性

荷電狀態(tài)（state of charge, SOC）是指蓄電池的剩余容量與其完全充電狀態(tài)的容量的比值，常用百分?jǐn)?shù)表示。蓄電池初始荷電狀態(tài)SOC0成正態(tài)分布，其概率密度函數(shù)為

蓄電池的充放電曲線，如圖2所示。

根據(jù)圖2所示曲線，對(duì)蓄電池進(jìn)行充放電擬合實(shí)驗(yàn)表明，其功率成正態(tài)分布。

圖2 蓄電池充放電曲線

2 微電網(wǎng)的隨機(jī)優(yōu)化模型

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

受天氣和環(huán)境等因素的影響，風(fēng)電機(jī)組和光電機(jī)組出力是不確定的。同時(shí)，由于蓄電池中能量的不確定，使得蓄電池的輸出功率具有隨機(jī)性。在滿足負(fù)荷需求條件下盡量避免或減少外電網(wǎng)供電的功率，即保證系統(tǒng)經(jīng)濟(jì)穩(wěn)定運(yùn)行。為此，本文采用了概率的形式來描述風(fēng)電機(jī)組、光電機(jī)組出力和蓄電池供電的隨機(jī)性，提出微電網(wǎng)供電的隨機(jī)優(yōu)化模型?；诖耍㈦娋W(wǎng)以用電費(fèi)用最小為目標(biāo)的隨機(jī)優(yōu)化模型的目標(biāo)函數(shù)可表示為

式中，cWT為風(fēng)電系統(tǒng)出力的成本系數(shù)、cPV為光伏發(fā)電系統(tǒng)出力的成本系數(shù)、cbat為蓄電池出力的成本系數(shù)、cG為外網(wǎng)系統(tǒng)出力的成本系數(shù)。為風(fēng)電機(jī)組出力、為光電機(jī)組出力、為蓄電池出力、為常規(guī)機(jī)組出力。T=1表示白天光電機(jī)組有功率輸出的情況，T=2表示晚上光電機(jī)組無功率輸出的情況。

2.2 約束條件

依據(jù)微電網(wǎng)系統(tǒng)運(yùn)行的具體要求，應(yīng)滿足如下約束條件。

（1）功率平衡約束

（2）運(yùn)行約束

式中，PWT.max表示風(fēng)力發(fā)電系統(tǒng)可輸出的最大功率；PPV.max表示光伏發(fā)電系統(tǒng)可輸出的最大功率；表示儲(chǔ)能蓄電池T時(shí)段可輸出的最大功率；γPV，γWT，γbat分別為事先給定的置信水平；P{·}為不等式約束成立的概率。當(dāng)T=1時(shí)，蓄電池出力由初始荷電狀態(tài)決定，具有不確定性；當(dāng)T=2時(shí)，

3 隨機(jī)優(yōu)化模型的SAA近似及粒子群算法優(yōu)化求解

3.1 隨機(jī)優(yōu)化模型的SAA近似

根據(jù)隨機(jī)優(yōu)化模型的復(fù)雜性，需要將優(yōu)化模型中的不確定因素轉(zhuǎn)化為確定性因素，這里采用抽樣平均近似（SAA）的方法。

計(jì)算得出，來評(píng)價(jià)其對(duì)應(yīng)的約束條件是否得到滿足。

可計(jì)算的確定性條件由不確定性條件經(jīng)以上步驟轉(zhuǎn)化而來：

將不確定性的優(yōu)化問題經(jīng)SAA處理轉(zhuǎn)化為確定性的優(yōu)化問題：

式（17）即為微電網(wǎng)隨機(jī)優(yōu)化調(diào)度模型，它考慮了風(fēng)電、光伏發(fā)電以及蓄電池隨機(jī)性的。

3.2 粒子群算法及其實(shí)現(xiàn)

1）基本粒子群算法

PSO是由Kennedy博士和Eberhart博士于l995年提出的一種基于群體智能演化的計(jì)算技術(shù)[14]。其基本機(jī)理是通過群體之間的信息共享和個(gè)體自身經(jīng)驗(yàn)總結(jié)來修正個(gè)體行動(dòng)策略，最終求取優(yōu)化問題的解[15]。

現(xiàn)假設(shè)Xi＝(xi1,xi2,…,xin)為粒子i的當(dāng)前位置，Vi＝(vi1,vi2,…,vin)為粒子i的當(dāng)前飛行速度，Pi＝(pi1,pi2,…,pim)為粒子i的歷史最優(yōu)位置，Pg＝(pg1,pg2,…,pgm)為全局最優(yōu)位置。本粒子群算法的進(jìn)化方程可描述為

式中，下標(biāo)i表示粒子i，j表示粒子的第j維，t表示第t代，w慣性權(quán)重系數(shù)，c1∈U(0,2),c2∈U (0,2)為加速常數(shù)，r1∈U(0,1),r2∈U (0,1)為兩個(gè)相互獨(dú)立的隨機(jī)函數(shù)。

2）粒子群算法的實(shí)現(xiàn)

基本PSO算法的實(shí)現(xiàn)步驟如下[16]：

（1）各粒子的飛行速度與位置隨機(jī)初始化。

（2）評(píng)價(jià)每個(gè)粒子的適應(yīng)度（包括約束適應(yīng)度和目標(biāo)函數(shù)適應(yīng)度）。

（3）比較當(dāng)前最優(yōu)位置Pi的適應(yīng)度與各個(gè)粒子的適應(yīng)度，選擇較優(yōu)，更新為當(dāng)前最優(yōu)位置。

（4）比較全局歷史最優(yōu)位置Pg的適應(yīng)度與各個(gè)粒子的適應(yīng)度，選擇較優(yōu)，更新當(dāng)前的全局最優(yōu)位置。

（5）對(duì)粒子的速度和位置根據(jù)式（18）式和（19）進(jìn)行進(jìn)化。

（6）如未達(dá)到預(yù)先設(shè)定的停止準(zhǔn)則（通常設(shè)置為最大迭代次數(shù)Gmax），則返回步驟（2），若達(dá)到則停止。

4 算例分析

從微電網(wǎng)系統(tǒng)的隨機(jī)優(yōu)化調(diào)度模型的合理性和有效性出發(fā)來，假設(shè)某一地區(qū)內(nèi)有微電網(wǎng)，它和外網(wǎng)同時(shí)為這地區(qū)提供電能。風(fēng)電機(jī)組采用金風(fēng)82/1500kW機(jī)組，其額定容量為1500kW，切入風(fēng)速、額定風(fēng)速、切出風(fēng)速分別是3、10.3、22m/s；光伏發(fā)電系統(tǒng)由晶硅太陽能電池板組成，基本技術(shù)參數(shù)為：標(biāo)稱功率230wp，開路電壓36.6V，峰值點(diǎn)電壓30.6V，短路電流7.98A，峰值點(diǎn)電流7.51A，光伏并網(wǎng)逆變器轉(zhuǎn)換效率為90%，系統(tǒng)連接效率為95%，光電轉(zhuǎn)換效率為15%，在光照強(qiáng)度S=1000w/m2、T=25℃時(shí)最大輸出功率為2MW；采用GEFC公司125V 200A VRB釩電池電堆，其額定容量為50kW·h，其額定功率為25kW，放電效率認(rèn)為是100%，充電效率采用雙程效率為72%，每小時(shí)自放電率為0.01%。成本系數(shù)分別為cWT=0.21，cPV=0.35，cbat=-0.55，cG=0.50；在以下算例分析中均采用PSO對(duì)目標(biāo)函數(shù)進(jìn)行求解。

4.1 隨機(jī)優(yōu)化模型對(duì)調(diào)度的影響

在隨機(jī)優(yōu)化調(diào)度模型中風(fēng)電系統(tǒng)、光伏發(fā)電系統(tǒng)和蓄電池按相等原則設(shè)定0.98為置信水平，設(shè)定為1500為M onte Carlo的抽樣次數(shù)，最優(yōu)調(diào)度方案和用電成本經(jīng)PSO算法求解得到，如表1所示；同時(shí)與確定性模型做了比較。

表1 不同決策方式下的調(diào)度方案和用電成本

從表1中可以看出，與確定性模型相比，隨機(jī)優(yōu)化模型的風(fēng)電、光伏發(fā)電和蓄電池的輸出功率均有所提高。確定性方案比隨機(jī)優(yōu)化方案的成本高169.46元，表明隨機(jī)優(yōu)化模型經(jīng)濟(jì)優(yōu)勢(shì)明顯。

4.2 置信度水平和抽樣次數(shù)對(duì)隨機(jī)優(yōu)化調(diào)度的影響

針對(duì)本文提出的隨機(jī)優(yōu)化調(diào)度模型，設(shè)置機(jī)會(huì)約束的不同的置信水平和SAA中不同的Monte Carlo抽樣次數(shù)，假設(shè)γWT= γPV=γbat，再利用PSO算法求解目標(biāo)函數(shù)得到優(yōu)化方案。表2和圖3為不同置信水平下和不同抽樣次數(shù)下，某個(gè)微電網(wǎng)系統(tǒng)白天時(shí)與外網(wǎng)同時(shí)為負(fù)載提供電能時(shí)的調(diào)度方案。

表2 不同置信水平下和不同抽樣次數(shù)下的調(diào)度方案和用電成本

從表2中可以看出，各機(jī)組的輸出功率在不同的置信水平下各不相同，適當(dāng)降低置信水平γWT/γPV/γbat，增加了風(fēng)電、光伏發(fā)電和蓄電池的輸出功率，減少了外網(wǎng)輸出功率，使負(fù)載的用電成本降低，證明機(jī)會(huì)約束是有效的、合理的。從圖3中可以看出各機(jī)組輸出功率在不同的抽樣次數(shù)下改變較大，增多Monte Carlo抽樣次數(shù)，將減小SAA的近似誤差。當(dāng)抽樣次數(shù)超過1500，優(yōu)化計(jì)算出的發(fā)電成本值基本趨于穩(wěn)定。因此，在計(jì)算中取1500為合理的抽樣次數(shù)，微電網(wǎng)系統(tǒng)輸出功率增加，外網(wǎng)輸出功率減小，達(dá)到減小負(fù)載用電成本的目的。所以在風(fēng)光儲(chǔ)微電網(wǎng)系統(tǒng)的實(shí)際運(yùn)行中，為減少負(fù)載用電費(fèi)用，應(yīng)合理的設(shè)置置信度水平和抽樣次數(shù)，以便更合理的分配微電網(wǎng)系統(tǒng)的輸出功率。

圖3 不同置信水平下和不同抽樣次數(shù)下的優(yōu)化調(diào)度結(jié)果

5 結(jié)論

由于微電網(wǎng)中新能源發(fā)電的隨機(jī)波動(dòng)性，使得微電網(wǎng)的調(diào)度與運(yùn)行非常復(fù)雜。為了解決微電網(wǎng)中分布式電源出力的優(yōu)化配置問題，本文在滿足系統(tǒng)功率平衡約束和運(yùn)行約束的條件下建立基于機(jī)會(huì)約束理論的微電網(wǎng)經(jīng)濟(jì)優(yōu)化數(shù)學(xué)模型，并采用抽樣平均近似（SAA）及粒子群優(yōu)化算法（PSO）進(jìn)行簡(jiǎn)化和尋優(yōu)求解，實(shí)現(xiàn)了微電網(wǎng)中微電源出力的優(yōu)化，使得系統(tǒng)用電總費(fèi)用最小。該方法不但滿足微電網(wǎng)的供電要求，還可有效地評(píng)估微電網(wǎng)的經(jīng)濟(jì)性。算例分析表明了所提出方法的有效性。

[1] 趙宏偉,吳濤濤.基于分布式電源的微網(wǎng)技術(shù)[J].電力系統(tǒng)及其自動(dòng)化學(xué)報(bào), 2008, 20(1): 121-128.

[2] 王凌,李培強(qiáng),李欣然,等.微電源建模及其在微電網(wǎng)仿真中的應(yīng)用[J].電力系統(tǒng)及其自動(dòng)化學(xué)報(bào), 2010,22(3): 32-38.

[3] 錢科軍,袁越,石曉丹,等.分布式發(fā)電的環(huán)境效益分析[J].中國(guó)電機(jī)工程學(xué)報(bào), 2008, 28(29): 11-15.

[4] HATZIARGYRIOU N, ASANO H, IRAVANI R, et al.Microgrids[J]. IEEE Power and Energy Magazine,2007, 5(4): 78-94.

[5] 戴上,張焰,祝達(dá)康.含有微電網(wǎng)的配電網(wǎng)規(guī)劃方法[J].電力系統(tǒng)動(dòng)化, 2010, 34(22): 41-45.

[6] 解翔,袁越,李振杰.含微電網(wǎng)的新型配電網(wǎng)供電可靠性分析[J].電力系統(tǒng)自動(dòng)化, 2011, 35(9): 67-72.

[7] 楊琦,馬世英,唐曉駿,等.微電網(wǎng)規(guī)劃評(píng)價(jià)指標(biāo)體系構(gòu)建與應(yīng)用[J].電力系統(tǒng)自動(dòng)化, 2012, 36(9): 13-17.

[8] 丁明,包敏,吳紅斌,等.復(fù)合能源分布式發(fā)電系統(tǒng)的機(jī)組組合問題[J].電力系統(tǒng)自動(dòng)化, 2008, 32(6): 46- 50.

[9] 李磊,陳海牛,呂亞洲.微電網(wǎng)發(fā)電備用配置優(yōu)化問題[J].華東電力, 2011, 39(8): 1153-1156.

[10] BILLINTON Y G R. Adequacy assessment of generating systems containing wind power considering wind speed correlation[J]. IET Renewable Power Generation, 2008, 217-226.

[11] WANG L, SINGH C.Balancing risk and cost in fuzzy economic dispa-tch including wind power penetration based on Particle Swarm Optimization[J]. Electric Power Systems Research, 2008, 78(8): 1361-1368.

[12] JOHNSON G L. Wind Energy System, 2001.

[13] 王敏.布式電源的概率建模及其對(duì)電力系統(tǒng)的影響[D]. [合肥工業(yè)大學(xué)博士學(xué)位論文],合肥,電力系統(tǒng)及其自動(dòng)化, 2010.

[14] KENNEDY J. EBERHART R. Particle swarm optimization [C]//Proceedings of the 1995 IEEE International Conference on Neural Networks, Perth,1995.

[15] 劉自發(fā),張建華.基于自適應(yīng)小生境粒子群優(yōu)化算法的電系統(tǒng)無功優(yōu)化.電力自動(dòng)化設(shè)備, 2009, 29(11):27-30.

[16] 劉偉,梁新蘭,安曉龍.基于蜜蜂進(jìn)化型粒子群算法的電力系統(tǒng)無功優(yōu)化.電力系統(tǒng)保護(hù)與控捌, 2010,38(7): 16-21.