楊文榮， 馬曉燕， 徐茂林， 邊鑫磊

(省部共建電工裝備可靠性與智能化國(guó)家重點(diǎn)實(shí)驗(yàn)室， 河北工業(yè)大學(xué)， 天津 300130)

1 引言

因全球能源危機(jī)和環(huán)境污染問(wèn)題日益突出，太陽(yáng)能、風(fēng)能等可再生清潔能源得到迅速發(fā)展。風(fēng)光出力的間歇性和隨機(jī)性以及負(fù)荷的波動(dòng)性等不確定因素威脅著電力系統(tǒng)的安全穩(wěn)定運(yùn)行，同時(shí)在孤島、偏遠(yuǎn)區(qū)域和微電源滲透率高的一些地區(qū)，難以實(shí)現(xiàn)供電要求。微電網(wǎng)是一種將分布式電源、儲(chǔ)能裝置、負(fù)荷和保護(hù)裝置等結(jié)合起來(lái)的小型發(fā)配電系統(tǒng)，可有效解決上述問(wèn)題[1-3]。

微電網(wǎng)的優(yōu)化調(diào)度即在滿足負(fù)荷需求和微電源約束條件等的前提下，使微電源的發(fā)電成本、停電損失、污染排放等達(dá)到最小。目前大部分研究方向集中在模型的優(yōu)化、調(diào)度策略的完善或有效的智能算法[4,5]。文獻(xiàn)[6]利用監(jiān)控和數(shù)據(jù)采集系統(tǒng)(SCADA)控制孤島型智能微電網(wǎng)，采用混合線性整數(shù)規(guī)劃算法(MILP)計(jì)算各分布式電源的出力，以達(dá)到最優(yōu)的方案。但由于目標(biāo)函數(shù)和約束條件及微電網(wǎng)系統(tǒng)的復(fù)雜性，MILP有一定局限性。文獻(xiàn)[7]提出了一個(gè)基于凈現(xiàn)值的主從式目標(biāo)函數(shù)，通過(guò)二次規(guī)劃算法和粒子群算法(PSO)最小化運(yùn)行成本和停電損失，采用模糊邏輯工具箱在Pareto曲線中尋找“從目標(biāo)函數(shù)”中最優(yōu)策略，然后運(yùn)用PSO的不斷迭代最終解決“主目標(biāo)函數(shù)”中的優(yōu)化問(wèn)題，但PSO易陷入局部最優(yōu)且收斂速度慢。文獻(xiàn)[8]考慮光伏發(fā)電、風(fēng)力發(fā)電和內(nèi)燃機(jī)發(fā)電，用蓄電池來(lái)儲(chǔ)存發(fā)電單元的剩余電量。仿真證明蓄電池在運(yùn)行過(guò)程中根據(jù)可再生能源和負(fù)荷的需要自由充放電比限制蓄電池充放電策略運(yùn)行成本要低。

文獻(xiàn)[9]為一種分級(jí)式控制策略，針對(duì)PSO可能不收斂以及參數(shù)過(guò)多難以確定等問(wèn)題，利用量子機(jī)制和波函數(shù)來(lái)表示解的情況，通過(guò)求解薛定諤函數(shù)和蒙特卡洛模型找出粒子的相應(yīng)位置，即采用量子粒子群算法(QPSO)優(yōu)化微電源的出力，使系統(tǒng)的經(jīng)濟(jì)成本最小。仿真證明QPSO算法優(yōu)于PSO，但此模型未考慮污染物的排放。文獻(xiàn)[10]以系統(tǒng)經(jīng)濟(jì)和環(huán)保成本最小為優(yōu)化目標(biāo)，從臨時(shí)Pareto解中篩選出N個(gè)擁擠距離最大的個(gè)體作為精英保留，采用此改進(jìn)遺傳算法對(duì)微電源的出力進(jìn)行優(yōu)化，證明了改進(jìn)算法可以跳出局部最優(yōu)解，且多目標(biāo)函數(shù)的優(yōu)化結(jié)果比單目標(biāo)函數(shù)更經(jīng)濟(jì)環(huán)保。文獻(xiàn)[6-10]的算法對(duì)于復(fù)雜超高維的問(wèn)題難以求得理想效果，且考慮單因素不利于微電網(wǎng)的調(diào)度研究。

鳥(niǎo)群算法(Bird Swarm Algorithm，BSA)是Xian-Bing Meng等人在2015年提出的一種新型生物啟發(fā)式全局優(yōu)化算法，來(lái)源于對(duì)鳥(niǎo)群覓食、警惕和飛行行為的模仿以及對(duì)鳥(niǎo)群覓食過(guò)程中共享信息的研究，BSA具備調(diào)節(jié)參數(shù)少、收斂精度高和魯棒性能好等特點(diǎn)，已驗(yàn)證BSA在優(yōu)化函數(shù)方面優(yōu)于PSO和微分進(jìn)化算法(DE)[11]。因BSA的認(rèn)知行為和群體行為調(diào)節(jié)參數(shù)易導(dǎo)致進(jìn)化種群收斂精度和迭代次數(shù)偏大，通過(guò)在BSA基礎(chǔ)上引入線性微分遞減策略的慣性權(quán)重和線性調(diào)整的學(xué)習(xí)系數(shù)，來(lái)提高收斂速度和搜索精度。本文采用改進(jìn)的鳥(niǎo)群算法(IBSA)對(duì)微電網(wǎng)并網(wǎng)運(yùn)行模式下的夏季典型日進(jìn)行調(diào)度分析，建立了考慮微電源的經(jīng)濟(jì)性和環(huán)保性的優(yōu)化調(diào)度模型，最后通過(guò)各目標(biāo)函數(shù)下微電網(wǎng)的調(diào)度情況和IBSA、BSA算法的對(duì)比，驗(yàn)證了本文調(diào)度模型和IBSA算法的正確性。

2 并網(wǎng)運(yùn)行方式下微電網(wǎng)優(yōu)化調(diào)度數(shù)學(xué)模型

本文采用的分布式電源為可再生清潔不可控機(jī)組光伏電池(Photovoltaic Cells，PV)和風(fēng)力發(fā)電機(jī)(Wind Turbine，WT)、效率高污染小的可控機(jī)組燃料電池(Fuel Cell，F(xiàn)C)和微型燃?xì)廨啓C(jī)(Micro Turbine，MT)、平抑可再生能源發(fā)電功率波動(dòng)的儲(chǔ)能蓄電池(Storage Battery，BT)，微電源的出力模型如文獻(xiàn)[12-14]中所述。微電網(wǎng)優(yōu)化調(diào)度的目標(biāo)是在各種約束條件下優(yōu)化微電源的出力使得系統(tǒng)的經(jīng)濟(jì)性和環(huán)保性達(dá)到最好。

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

2.1.1 微電網(wǎng)系統(tǒng)的經(jīng)濟(jì)性成本最低

目標(biāo)函數(shù)為：

minC1,COST=min(CFC&MT+COM+CDP+Cgrid)

(1)

(2)

ηFC(t)=-0.0023PFC(t)+0.6735

(3)

(4)

(5)

(6)

(7)

(8)

式中，CFC&MT、COM、CDP和Cgrid分別為分布式電源的燃料、運(yùn)行維護(hù)、折舊和與主網(wǎng)電能互購(gòu)成本[14,15]；N為微電源種類(lèi)數(shù)；C為天然氣價(jià)格；LHV為天然氣低熱值；PFC(t)、PMT(t)分別為t時(shí)刻FC和MT的輸出功率；ηFC(t)、ηMT(t)分別為t時(shí)刻FC和MT的發(fā)電效率；KOM,i、Pi(t)分別為第i種分布式電源的運(yùn)行維護(hù)系數(shù)和在t時(shí)刻的輸出功率；ADCCi為微電源i的年均折舊資本；Ccost,i為單位容量初始安裝成本；ki、ri、Li、Pfci分別為微電源的容量因數(shù)、年折舊利率(6%)、折舊年限和額定功率；cbuy(t)、ccell(t)分別為t時(shí)刻的購(gòu)售電電價(jià)；Pbuy(t)、Psell(t)分別為t時(shí)刻的購(gòu)售電功率。

2.1.2 微電網(wǎng)系統(tǒng)的污染物處理費(fèi)用最小

本文污染物指CO2、SO2和NOx，目標(biāo)函數(shù)為：

minC2,COST=CDG,i+Cgrid,i

(9)

(10)

(11)

式中，CDG,i、Cgrid,i分別為微電源i和主網(wǎng)的排污治理費(fèi)用(微電源中產(chǎn)生污染的為FC和MT)；αj為第j種污染物單位排污處理成本；βi,j、βgrid,j分別為微電源i和主網(wǎng)的第j種污染物排放系數(shù)；Pgrid(t)為t時(shí)刻微電網(wǎng)與主網(wǎng)的交互功率；Δt為時(shí)間間隔。

2.1.3 微電網(wǎng)系統(tǒng)的綜合效益最好

綜合考慮經(jīng)濟(jì)成本和環(huán)境污染因素，由于兩目標(biāo)函數(shù)量綱相同，故用線性加權(quán)求和法將多目標(biāo)轉(zhuǎn)化為單目標(biāo)，再用二元對(duì)比定權(quán)法確定權(quán)系數(shù)[16]。本文取經(jīng)濟(jì)性重要性等級(jí)“稍稍”在環(huán)保性之上，通過(guò)建立重要性定性排序標(biāo)度矩陣E和經(jīng)過(guò)一致性檢驗(yàn)，將模糊語(yǔ)氣算子轉(zhuǎn)換為隸屬度指標(biāo)的非歸一化權(quán)重向量 [1，0.739]，歸一化后得到的指標(biāo)的權(quán)重向量W=[0.575，0.425]，目標(biāo)函數(shù)為：

minC=0.575C1,COST+0.425C2,COST

(12)

2.2 約束條件

(1)功率平衡約束

(13)

式中，PBT(t)為蓄電池充放電功率，放電為正，充電為負(fù)；PPV(t)、PWT(t)分別為光伏電池、風(fēng)機(jī)的出力；PL(t)為微電網(wǎng)負(fù)荷。

(2)微電源輸出的功率約束

Pi.min≤Pi≤Pi.max

(14)

式中，Pi.max、Pi.min分別為微電源的有功出力上、下限。

(3)微電網(wǎng)和主網(wǎng)間的交互功率約束

Pgrid,min≤Pgrid≤Pgrid,max

(15)

式中，Pgrid,max、Pgrid,min分別為傳輸容量上、下限。

(4)儲(chǔ)能蓄電池的運(yùn)行約束

PBT,i,min≤PBT,i≤PBT,i,max

(16)

EBT(0)=EBT(T)

(17)

SOCmin≤SOC≤SOCmax

(18)

式中，PBT,i,max、PBT,i,min分別為蓄電池充放電功率最大值、最小值；SOCmax、SOCmin分別為蓄電池荷電狀態(tài)允許的最大值、最小值；EBT(0)、EBT(T)分別為蓄電池一個(gè)調(diào)度周期的初始儲(chǔ)能和最終儲(chǔ)能。

3 改進(jìn)鳥(niǎo)群算法及微電網(wǎng)優(yōu)化調(diào)度求解步驟

3.1 鳥(niǎo)群算法仿生學(xué)基本原理

BSA算法是根據(jù)鳥(niǎo)類(lèi)覓食、警惕和飛行行為衍生出的生物啟發(fā)式算法，這些社會(huì)行為和活動(dòng)的本質(zhì)是群體智能，鳥(niǎo)群從社會(huì)活動(dòng)中獲取最佳食物的位置即對(duì)應(yīng)于目標(biāo)函數(shù)中尋找最優(yōu)解。鳥(niǎo)類(lèi)的行為簡(jiǎn)化規(guī)則如下：

(1)規(guī)則1：每只鳥(niǎo)隨機(jī)決定在警惕和覓食行為之間轉(zhuǎn)換身份。

(2)規(guī)則2：當(dāng)選擇覓食時(shí)，鳥(niǎo)類(lèi)及時(shí)記錄和更新該個(gè)體和種群先前最好的覓食位置，并以信息的方式共享于整個(gè)種群，以此來(lái)更好地尋找食物。

(3)規(guī)則3：當(dāng)保持警惕狀態(tài)時(shí)，每只鳥(niǎo)都試圖向種群的中心位置飛去，此行為受到種群間競(jìng)爭(zhēng)的干擾，高儲(chǔ)備量的鳥(niǎo)類(lèi)更容易飛到種群的中心位置。

(4)規(guī)則4：鳥(niǎo)類(lèi)會(huì)定期飛到另一個(gè)區(qū)域，當(dāng)?shù)竭_(dá)終點(diǎn)時(shí)，鳥(niǎo)類(lèi)將會(huì)在生產(chǎn)者和乞食者之間轉(zhuǎn)換，食物儲(chǔ)存量多的為生產(chǎn)者，反之為乞食者。其他鳥(niǎo)類(lèi)則隨機(jī)選擇成為生產(chǎn)者或乞食者。

(5)規(guī)則5：生產(chǎn)者積極尋找食物，乞食者則隨機(jī)跟隨生產(chǎn)者尋找食物。

鳥(niǎo)群規(guī)模設(shè)為N，飛行空間維度為D，xit(i∈[1,2,…,N])代表第i只鳥(niǎo)在t時(shí)刻的覓食或飛行位置。

3.1.1 覓食行為

規(guī)則1為隨機(jī)決策，若rand(0,1)

每只鳥(niǎo)由自身和群體的經(jīng)驗(yàn)來(lái)尋找食物。規(guī)則2用數(shù)學(xué)公式表示如下：

(19)

式中，j∈[1,2,…,D]；rand(0,1)代表(0，1)之間的獨(dú)立均勻分布數(shù)；C和S(C,S>0)分別為認(rèn)知系數(shù)和社會(huì)系數(shù)；pi,j為第i只鳥(niǎo)的先前最優(yōu)位置；gj為種群的先前最優(yōu)位置。

3.1.2 警惕行為

根據(jù)規(guī)則3，鳥(niǎo)類(lèi)在試圖飛向種群的中心時(shí)會(huì)不可避免地與其他鳥(niǎo)類(lèi)發(fā)生競(jìng)爭(zhēng)。因此每只鳥(niǎo)不會(huì)直接飛向種群的中心，這些行為可表示為：

(20)

(21)

(22)

式中，k(k≠i)為[1,N]之間的隨機(jī)正整數(shù)；a1，a2∈[0,2]；pFiti、pFitk分別為第i、k只鳥(niǎo)的最佳適應(yīng)度值；sumFit為種群最佳適應(yīng)度值之和；ε為用來(lái)避免零分割的計(jì)算機(jī)中很小的常數(shù)；meanj為整個(gè)種群中第j維的平均位置。

3.1.3 飛行行為

因捕食者的威脅、覓食或其他原因，鳥(niǎo)類(lèi)將會(huì)飛向另一個(gè)區(qū)域，此時(shí)他們將會(huì)重新覓食。一些鳥(niǎo)作為生產(chǎn)者來(lái)尋找食物，另一些則根據(jù)這些生產(chǎn)者找到的食物信息來(lái)尋找食物，規(guī)則4可以區(qū)分出鳥(niǎo)類(lèi)生產(chǎn)者和乞食者，“生產(chǎn)者”的數(shù)學(xué)表達(dá)式為：

(23)

式中，randn(0,1)為服從均值為0、標(biāo)準(zhǔn)差為1的Gaussian中的一個(gè)隨機(jī)數(shù)。

“乞食者”的數(shù)學(xué)表達(dá)式為：

(24)

式中，

3.2 改進(jìn)鳥(niǎo)群算法(IBSA)

3.2.1 線性微分遞減策略的慣性權(quán)重的改進(jìn)

慣性權(quán)重是鳥(niǎo)群算法的重要改進(jìn)參數(shù)，通過(guò)設(shè)置慣性權(quán)重可以改變上一時(shí)刻覓食位置對(duì)當(dāng)前時(shí)刻位置的影響。迭代初期較大的慣性權(quán)重能使BSA算法保持較強(qiáng)的全局搜索能力，迭代后期較小的慣性權(quán)重有利于算法更精確的局部搜索。為了更好地平衡BSA算法的全局與局部搜索能力，本文采用線性微分遞減策略[17]，計(jì)算公式為：

(25)

式(25)可轉(zhuǎn)化為：

(26)

式中，wmax、wmin分別為慣性權(quán)重的最大值和最小值，本文取wmax=0.9，wmin=0.4；tmax為最大迭代次數(shù)；t為當(dāng)前迭代次數(shù)。隨著迭代的進(jìn)行，慣性權(quán)重由0.9線性遞減到0.4。

3.2.2 線性調(diào)整學(xué)習(xí)系數(shù)

BSA算法中C和S分別代表認(rèn)知系數(shù)和社會(huì)系數(shù)，學(xué)習(xí)因子可以調(diào)整自身認(rèn)知經(jīng)驗(yàn)和社會(huì)經(jīng)驗(yàn)的比重值。鳥(niǎo)群搜索覓食的前期C取較大值，S取較小值，增加認(rèn)知經(jīng)驗(yàn)的比重，使鳥(niǎo)類(lèi)的全局搜索能力增強(qiáng)；搜尋覓食的后期C取較小值，S取較大值，增加搜尋食物的社會(huì)經(jīng)驗(yàn)的比重，使得局部搜索增強(qiáng)，線性調(diào)整的學(xué)習(xí)系數(shù)更新公式為：

(27)

(28)

式中，Ce=Ss=0.5；Cs=Se=2.5。

改進(jìn)后的覓食公式為:

(29)

3.3 IBSA算法求解微電網(wǎng)并網(wǎng)優(yōu)化調(diào)度步驟

針對(duì)第2節(jié)建立的微電網(wǎng)優(yōu)化調(diào)度模型，每個(gè)微電源的出力為待優(yōu)化的參數(shù)，鳥(niǎo)群中每個(gè)個(gè)體所處的空間位置為優(yōu)化變量，適應(yīng)度函數(shù)的大小表示該位置的優(yōu)劣，根據(jù)鳥(niǎo)群覓食、警惕和飛行行為來(lái)不斷更新鳥(niǎo)類(lèi)個(gè)體的空間位置，直至尋找到群體最佳位置，即微電源的最優(yōu)出力[18,19]。具體步驟如下：

(1)初始化BSA算法中各參數(shù)和不同種類(lèi)微電源的出力上下限，輸入夏季典型日24h的用電情況。利用式(30)初始化微電源在空間的位置。

xit=lb+(ub-lb)·rand(1,D)

(30)

式中，xit為微電源i在t時(shí)刻的位置；ub、lb分別為微電源的出力上、下限。

(2)計(jì)算各目標(biāo)函數(shù)的適應(yīng)度值，并從中隨機(jī)選取個(gè)體最優(yōu)和全局最優(yōu)位置。

(3)采用線性微分遞減策略改進(jìn)慣性權(quán)重，并線性調(diào)整認(rèn)知系數(shù)和社會(huì)系數(shù)。

(4)依據(jù)鳥(niǎo)類(lèi)覓食、警惕和飛行行為更新種群，限定微電源的出力滿足式(13)～式(18)約束，將微電源的功率平衡等式約束轉(zhuǎn)換為不等式判定，在超過(guò)范圍之后加入不同程度的懲罰因子。

(5)計(jì)算新種群適應(yīng)度值，與前一時(shí)刻適應(yīng)度值作比較，剔除劣解，留下優(yōu)解，更新歷史最優(yōu)解。

(6)判斷是否滿足終止條件，若達(dá)到要求，輸出最終最優(yōu)解，否則t+1返回至步驟(4)繼續(xù)迭代，直至滿足終止條件或者達(dá)到最優(yōu)解[20]。

4 算例分析

本文以工作在并網(wǎng)運(yùn)行模式下的包含PV、WT、FC、MT、BT和負(fù)荷的微電網(wǎng)系統(tǒng)為例，其中PV和WT工作在最大功率點(diǎn)跟蹤模式(MPPT)，本文以最大化利用可再生能源為前提，基于改進(jìn)的鳥(niǎo)群算法，采用MATLAB R2014a編寫(xiě)了適用于微電網(wǎng)算例的優(yōu)化調(diào)度程序。

4.1 微電網(wǎng)基礎(chǔ)數(shù)據(jù)

本文計(jì)算周期為1天，24個(gè)時(shí)段，每小時(shí)代表1個(gè)時(shí)段，圖1為夏季典型日負(fù)荷出力曲線圖，全天溫度、風(fēng)速、光照強(qiáng)度數(shù)據(jù)曲線分別如圖2～圖4所示，由此計(jì)算出的PV和WT出力如圖5所示。微電網(wǎng)系統(tǒng)分布式電源參數(shù)、污染物處理排放系數(shù)以及全天分時(shí)電價(jià)分別如表1～表3所示。

圖1 24h負(fù)荷出力曲線圖Fig.1 Load output curve in 24h

圖2 24h溫度曲線圖Fig.2 Temperature curve in 24h

圖3 24h風(fēng)速曲線圖Fig.3 Wind speed curve in 24h

圖4 24h光照強(qiáng)度曲線圖Fig.4 Light intensity curve in 24h

圖5 24h風(fēng)光出力曲線圖Fig.5 WT and PV output curve in 24h

數(shù)學(xué)模型主要參數(shù)為：PV參數(shù)：標(biāo)準(zhǔn)測(cè)試條件下光照強(qiáng)度為1000W/m2，參考溫度為25℃，功率溫度系數(shù)為-0.0047/℃；WT參數(shù)：切入、切出和額定風(fēng)速分別為2.5m/s、25m/s和12m/s；FC、MT參數(shù)：天然氣價(jià)格取2.5元/m3，天然氣低熱值取9.7(kW·h)/m3；BT參數(shù)：最大儲(chǔ)能量150kW·h，每小時(shí)充放電功率限值為-30kW、30kW，初始荷電狀態(tài)為0.2，荷電狀態(tài)上、下限為1和0.2，充電為負(fù)， 放電為正；微電網(wǎng)與主網(wǎng)的最大交互功率為50kW， 從主網(wǎng)購(gòu)電為正，售電為負(fù)。

表2 污染物處理排放系數(shù)Tab.2 Pollutant discharge coefficient

表3 購(gòu)電和售電電價(jià)Tab.3 Purchase and sale of electricity price

4.2 仿真結(jié)果

4.2.1 微電網(wǎng)系統(tǒng)不同目標(biāo)函數(shù)下出力曲線分析

采用改進(jìn)的鳥(niǎo)群算法，分別對(duì)微電網(wǎng)的經(jīng)濟(jì)性、環(huán)保性和采用二元對(duì)比定權(quán)法的綜合性目標(biāo)函數(shù)進(jìn)行分析，PV和WT始終以MPPT模式最先出力，24h的優(yōu)化出力曲線如圖6～圖8所示，不同目標(biāo)函數(shù)下的成本對(duì)比如表4所示。

圖6 運(yùn)行成本最低情況下優(yōu)化調(diào)度曲線Fig.6 Optimal scheduling curve under lowest operating cost

圖7 污染物處理費(fèi)用最低情況下優(yōu)化調(diào)度曲線Fig.7 Optimal scheduling curve under lowest cost of pollutant treatment

圖8 綜合效益最大情況下優(yōu)化調(diào)度曲線Fig.8 Optimal scheduling curve under condition of comprehensive benefit

目標(biāo)函數(shù)運(yùn)行成本/元環(huán)境費(fèi)用/元運(yùn)行成本最小環(huán)境費(fèi)用最小綜合成本最低214055254000225724324313088531212

由圖6可知，在經(jīng)濟(jì)成本最低情況下由于FC的運(yùn)行成本低于MT，除電網(wǎng)電價(jià)谷時(shí)段，其余時(shí)刻優(yōu)化調(diào)度時(shí)優(yōu)先考慮FC。1∶00～7∶00時(shí)段，電網(wǎng)負(fù)荷較輕，電網(wǎng)電價(jià)處于谷時(shí)段，優(yōu)先考慮從主網(wǎng)購(gòu)電，主網(wǎng)發(fā)出的電量、FC和MT向負(fù)荷供電的剩余電量給BT充電。夏季典型日11∶00～14∶00、18∶00～21∶00時(shí)段處于用電高峰期，F(xiàn)C接近滿發(fā)，不足的功率由MT、BT放電和從主網(wǎng)購(gòu)電來(lái)提供。

由圖7可知，污染物處理費(fèi)用最低的模型下，在FC、MT和主網(wǎng)同時(shí)進(jìn)行調(diào)峰時(shí)，環(huán)境污染度最小的MT發(fā)揮了主要作用，其次是FC，最后是主網(wǎng)。夏季典型日10∶00～15∶00光照非常充足，此時(shí)光伏出力增加，8∶00、15∶00和晚間風(fēng)速很大，風(fēng)機(jī)出力增加，雖然此時(shí)蓄電池吸收能量，但是由于用電負(fù)荷很高， 8∶00以后蓄電池始終處于放電狀態(tài)。

由圖8可知，綜合考慮經(jīng)濟(jì)成本和環(huán)境污染處理費(fèi)用，在1∶00～7∶00用電谷時(shí)段，電價(jià)較低，此時(shí)優(yōu)先從主網(wǎng)購(gòu)電來(lái)滿足負(fù)荷需求，晚間風(fēng)速大，風(fēng)機(jī)出力增加，因此使用從電網(wǎng)購(gòu)電和微電源出力剩余電量給BT充電。夏季用電負(fù)荷普遍偏高，故在8∶00后BT一直保持放電狀態(tài)。11∶00～15∶00、18∶00～21∶00時(shí)段為用電高峰期，F(xiàn)C接近滿發(fā)來(lái)使運(yùn)行成本降低，其余時(shí)刻利用污染度小的MT優(yōu)先出力。

由表4可知，夏季典型日綜合成本最低目標(biāo)函數(shù)中，運(yùn)行成本與環(huán)境費(fèi)用均介于運(yùn)行成本最小與環(huán)境費(fèi)用最低目標(biāo)函數(shù)下的運(yùn)行成本與環(huán)境費(fèi)用之間，綜合考慮經(jīng)濟(jì)性與環(huán)境污染處理費(fèi)用的多目標(biāo)優(yōu)化調(diào)度很好地協(xié)調(diào)了經(jīng)濟(jì)與環(huán)境之間的平衡。

4.2.2 BSA、IBSA兩種智能算法對(duì)比分析

本文采用的BSA參數(shù)如下：遷移周期FQ=10，a1=a2=1.5,C=S=1.5；IBSA參數(shù)如下：FQ=10，a1=a2=1.5，其余參數(shù)3.2節(jié)已經(jīng)給出；兩種算法種群規(guī)模均為300，迭代300次。針對(duì)本文提出的微電網(wǎng)優(yōu)化模型，兩種算法均獨(dú)立運(yùn)行30次，求取平均值，對(duì)比結(jié)果如圖9和表5所示。

圖9 綜合目標(biāo)函數(shù)下的適應(yīng)度收斂曲線對(duì)比圖Fig.9 Contrast curve of fitness convergence curve under comprehensive objective function

算法運(yùn)行時(shí)間/s收斂次數(shù)運(yùn)行成本/元BSA3587182264797IBSA3582141225724

綜上可得，在同等種群規(guī)模和迭代次數(shù)條件下，IBSA算法相比BSA算法其收斂速度、收斂次數(shù)和尋優(yōu)精度更好，由此證明了IBSA算法的優(yōu)越性。

5 結(jié)論

本文建立了微電網(wǎng)并網(wǎng)運(yùn)行模式下考慮運(yùn)行成本和環(huán)境處理費(fèi)用的優(yōu)化調(diào)度模型，采用新型生物啟發(fā)式算法——改進(jìn)的鳥(niǎo)群算法分別優(yōu)化經(jīng)濟(jì)性、環(huán)保性和綜合性目標(biāo)函數(shù)下各分布式電源的出力，并比較了IBSA和BSA兩種算法的優(yōu)劣。結(jié)果表明，多目標(biāo)優(yōu)化調(diào)度有效協(xié)調(diào)了經(jīng)濟(jì)與環(huán)境之間的平衡，比單目標(biāo)情況下的結(jié)果更加經(jīng)濟(jì)環(huán)保和高效，同時(shí)IBSA算法的收斂速度和尋優(yōu)精度均優(yōu)于BSA算法。算例仿真結(jié)果驗(yàn)證了模型和IBSA算法的正確性。

[1] 宋璇坤， 韓柳， 鞠黃培， 等(Song Xuankun，Han Liu，Ju Huangpei，et al．)． 中國(guó)智能電網(wǎng)技術(shù)發(fā)展實(shí)踐綜述(A review on development practice of smart grid technology in China)[J]． 電力建設(shè)(Electric Power Construction)，2016，37(7)：1-11．

[2] 李鵬， 竇鵬沖， 李雨薇， 等(Li Peng，Dou Pengchong，Li Yuwei，et al．)． 微電網(wǎng)技術(shù)在主動(dòng)配電網(wǎng)中的應(yīng)用(Application of microgrid technology in active distribution network)[J]． 電力自動(dòng)化設(shè)備(Electric Power Automation Equipment)，2015，35(4)：8-16．

[3] 尤毅， 劉東， 鐘清， 等(You Yi，Liu Dong，Zhong Qing，et al．)． 主動(dòng)配電網(wǎng)優(yōu)化調(diào)度策略研究(Research on optimal schedule strategy for active distribution network)[J]． 電力系統(tǒng)自動(dòng)化(Automation of Electric Power Systems)，2014，38(9)：177-183．

[4] 王成山， 武震， 李鵬(Wang Chengshan，Wu Zhen，Li Peng)． 微電網(wǎng)關(guān)鍵技術(shù)研究(Research on key technologies of microgrid)[J]． 電工技術(shù)學(xué)報(bào)(Transactions of China Electrotechnical Society)，2014, 29(2)：1-12．

[5] 劉文勝， 毛曉明， 陳志勇(Liu Wensheng，Mao Xiaoming，Chen Zhiyong)． 微網(wǎng)優(yōu)化運(yùn)行研究進(jìn)展與展望(A review on microgrid optimal operation)[J]． 電工電能新技術(shù)(Advanced Technology of Electrical Engineering and Energy)，2012，31(2)：64-68．

[6] Lazar E， Petreus D， Etz R， et al． Optimal scheduling of an islanded microgrid based on minimum cost [A]. International Spring Seminar on Electronics Technology[C]. 2016. 290-295．

[7] Moradi M H， Eskandari M， Mahdi Hosseinian S． Operational strategy optimization in an optimal sized smart microgrid [J]． IEEE Transactions on Smart Grid，2015，6(3)：1087-1095．

[8] Toma L， Tristiu I， Bulac C， et al． Optimal generation scheduling strategy in a microgrid[A]. IEEE Transportation Electrification Conference and Expo, Asia-Pacific [C]. 2016. 491-496．

[9] 黃晨， 張靠社(Huang Chen, Zhang Kaoshe)． 基于量子粒子群算法的微電網(wǎng)優(yōu)化調(diào)度(Micro grid scheduling optimization based on quantum particle swarm optimization(PSO)algorithm)[J]． 電網(wǎng)與清潔能源(Power System and Clean Energy)，2016, 32(2)：72-76．

[10] 趙磊， 曾芬鈺， 王霜， 等(Zhao Lei，Zeng Fenyu， Wang Shuang， et al．)． 基于經(jīng)濟(jì)性與環(huán)保性的微電網(wǎng)多目標(biāo)優(yōu)化調(diào)度研究(Research on multi-objective optimal operation of microgrid based on economic and environmental protection)[J]． 高壓電器(High Voltage Apparatus)，2015, 51(6)：127-132．

[11] Meng X B， Gao X Z， Lu L H， et al． A new bio-inspired optimisation algorithm：Bird swarm algorithm [J]． Journal of Experimental & Theoretical Artificial Intelligence，2015, 28(4)：673-687．

[12] 張麗娟(Zhang Lijuan)． 含多種分布式電源的微網(wǎng)系統(tǒng)環(huán)保經(jīng)濟(jì)優(yōu)化調(diào)度研究(Research on environmental and economic optimal dispatch of microgrid with multi-distributed genertion)[D]．北京：華北電力大學(xué)(Beijing：North China Electric Power University)，2015．

[13] 肖浩， 裴瑋， 孔力， 等(Xiao Hao，Pei Wei，Kong Li，et al．)． 綜合能源微網(wǎng)運(yùn)行優(yōu)化調(diào)度方法研究(Research on optimal operation method of integrated energy microgrid system)[J]． 電工電能新技術(shù)(Advanced Technology of Electrical Engineering and Energy)，2016，35(12):1-11．

[14] 段程煜(Duan Chengyu)．基于粒子群算法的微電網(wǎng)優(yōu)化調(diào)度研究(The research on optimal scheduling of microgrid based on particle swarms optimization(PSO))[D]．北京：華北電力大學(xué)(Beijing：North China Electric Power University)，2014．

[15] 周曉燕， 劉天琪， 沈浩東， 等(Zhou Xiaoyan, Liu Tianqi, Shen Haodong, et al．)． 含多種分布式電源的微電網(wǎng)經(jīng)濟(jì)調(diào)度研究(Research on economic dispatch of microgrid with multi-distributed generation)[J] ． 電工電能新技術(shù)(Advanced Technology of Electrical Engineering and Energy)，2013，32(1)：5-8．

[16] 鐘宇峰， 黃民翔， 葉承晉(Zhong Yufeng，Huang Minxiang，Ye Chengjin，et al．)． 基于電池儲(chǔ)能系統(tǒng)動(dòng)態(tài)調(diào)度的微電網(wǎng)多目標(biāo)運(yùn)行優(yōu)化(Multi-objective optimization of microgrid operation based on dynamic dispatch of battery energy storage system)[J]． 電力自動(dòng)化設(shè)備(Electric Power Automation Equipment)，2014，34(6)：114-121．

[17] 梅真(Mei Zhen)． 基于改進(jìn)粒子群算法的微電網(wǎng)優(yōu)化運(yùn)行策略研究(Research on the optimal operation strategy of micro-grid based on improved PSO)[D]． 武漢：湖北工業(yè)大學(xué)(Wuhan: Hubei University of Technology)，2015．

[18] 崔東文， 金波(Cui Dongwen，Jin Bo)． 改進(jìn)鳥(niǎo)群算法及其在梯級(jí)水庫(kù)優(yōu)化調(diào)度中的應(yīng)用(Improved bird swarm algorithm and its application to reservoir optimal operation)[J]． 三峽大學(xué)學(xué)報(bào)(自然科學(xué)版)(Journal of China Gorges University(Natural Sciences))，2016，38(6): 7-14.

[19] 崔東文， 金波(Cui Dongwen，Jin Bo)． 鳥(niǎo)群算法-投影尋蹤回歸模型在多元變量年徑流預(yù)測(cè)中的應(yīng)用(Application of the bird swarm algorithm - projection pursuit regression model to prediction of multivariate annual runoff)[J]． 人民珠江(Pearl River)，2016，37(11)：26-30．

[20] 曾嶒， 彭春華， 王奎， 等(Zeng Ceng， Peng Chunhua, Wang Kui，et al．)． 基于鳥(niǎo)群算法的微電網(wǎng)多目標(biāo)運(yùn)行優(yōu)化(Multi-objective operation optimization of micro grid based on bird swarm algorithm)[J]． 電力系統(tǒng)保護(hù)與控制(Power System Protection and Control)，2016，44(13)：117-122．