周榮昱

(福州大學(xué)電氣工程與自動(dòng)化學(xué)院，福建 福州 350108)

1 引言

隨著化石能源危機(jī)及環(huán)境保護(hù)問(wèn)題的日益突出，分布式電源(Distributed Generation，DG)作為一種可再生清潔能源在電力工業(yè)中發(fā)展迅速。然而，DG具有強(qiáng)烈的隨機(jī)性、間歇性和難預(yù)測(cè)性導(dǎo)致其難以持續(xù)穩(wěn)定的輸出功率，大規(guī)模DG并網(wǎng)將給電力系統(tǒng)的安全可靠運(yùn)行帶來(lái)不利影響[1-2]，從而對(duì)電網(wǎng)的調(diào)度有了更高的要求。儲(chǔ)能技術(shù)具有對(duì)功率和能量的時(shí)間遷移能力，能有效改善DG的輸出特性，增強(qiáng)DG的可調(diào)度性，是解決DG大規(guī)模并網(wǎng)問(wèn)題及促進(jìn)資源利用的關(guān)鍵技術(shù)[3-6]。

目前，國(guó)內(nèi)外關(guān)于儲(chǔ)能容量配置的研究文獻(xiàn)[7]主要以綜合成本最低、電壓偏移量最小為優(yōu)化目標(biāo)建立數(shù)學(xué)模型。根據(jù)目標(biāo)函數(shù)和決策變量的不同，配電網(wǎng)儲(chǔ)能容量配置模型主要分為線性規(guī)劃模型和非線性規(guī)劃模型[8]。常用的儲(chǔ)能容量?jī)?yōu)化配置方法有內(nèi)點(diǎn)法[9]、離散傅立葉變換、加權(quán)移動(dòng)平均控制法和人工智能算法[10]等。在儲(chǔ)能優(yōu)化配置儲(chǔ)能過(guò)程中，文獻(xiàn)[11]將容量配置與充放電策略分別選作上層決策與下層決策,并以投資費(fèi)用最低和失負(fù)荷概率最小為子目標(biāo)，建立多目標(biāo)優(yōu)化配置模型。文獻(xiàn)[12]采用基于全壽命周期成本的優(yōu)化模型，并對(duì)不同類型儲(chǔ)能電池的年平均成本及循環(huán)效率進(jìn)行對(duì)比。文獻(xiàn)[13]對(duì)配電網(wǎng)的節(jié)點(diǎn)電壓偏離值和線路損耗最小為目標(biāo)函數(shù)，采用線性加權(quán)和法將多目標(biāo)問(wèn)題轉(zhuǎn)化為含權(quán)重系數(shù)的單目標(biāo)問(wèn)題。

本文建立以系統(tǒng)綜合成本最低、網(wǎng)損率最低、電壓穩(wěn)定性最高作為儲(chǔ)能配置優(yōu)化目標(biāo)，綜合考慮DG電源、負(fù)荷和儲(chǔ)能裝置相應(yīng)的模型約束條件，實(shí)現(xiàn)了將配電網(wǎng)穩(wěn)定經(jīng)濟(jì)運(yùn)行與儲(chǔ)能容量的合理配置建立相關(guān)模型，對(duì)維持配網(wǎng)功率平衡與降低綜合成本具有顯著意義。

2 儲(chǔ)能的經(jīng)濟(jì)效益模型

配電網(wǎng)工程項(xiàng)目投入使用后，其工作周期較長(zhǎng)，因此在配置儲(chǔ)能容量時(shí)，需要考慮不同時(shí)期的投入費(fèi)用。儲(chǔ)能系統(tǒng)的綜合投入費(fèi)用由建設(shè)階段、運(yùn)營(yíng)維護(hù)階段構(gòu)成，以下對(duì)各階段費(fèi)用做詳細(xì)分析。

2.1 一次投資成本

建設(shè)初期費(fèi)用即一次成本C1，包含超級(jí)電容與蓄電池購(gòu)置成本、功率轉(zhuǎn)換設(shè)備的技術(shù)投資成本、電池能量管理系統(tǒng)投資成本，計(jì)算公式為：

C1=mfba+nfuc+KpPn+KeEn

(1)

式中：C1一次投入成本，m、n分別為蓄電池與超級(jí)電容配置數(shù)量，fba為蓄電池購(gòu)置成本，fnc為超級(jí)電容購(gòu)置成本，Kp為池能量管理系統(tǒng)投資成本；功率轉(zhuǎn)換設(shè)備的單位成本，Pn為混合儲(chǔ)能額定功率，Ke為能量管理系統(tǒng)的單位成本，En為混合儲(chǔ)能額定容量。

2.2 運(yùn)營(yíng)維護(hù)成本

儲(chǔ)能系統(tǒng)在全壽命周期內(nèi)的運(yùn)營(yíng)成本C2包括對(duì)蓄電池和超級(jí)電容的維護(hù)成本Cmc與更換成本Cc：

C2=Cmc+Cc

(2)

Cmc=N(mKmcbafba+nKmcucfuc)

(3)

Cc=mQbafba+nQucfuc

(4)

式中，N為混合儲(chǔ)能系統(tǒng)設(shè)計(jì)服役時(shí)間；Kmcba、Kmcuc分別為蓄電池和超級(jí)電容單位運(yùn)維成本；Qba、Quc分別為需要更換的蓄電池和超級(jí)電容數(shù)量，考慮到超級(jí)電容具有較高的循環(huán)次數(shù)和使用壽命，因此在分析時(shí)可不考慮超級(jí)電容的更換與維護(hù)費(fèi)用。

2.3 退役處置成本

蓄電池在達(dá)到設(shè)計(jì)使用期限后，進(jìn)入退役處置階段，這一階段所產(chǎn)生的費(fèi)用統(tǒng)稱為退役處置成本，其費(fèi)用主要為微網(wǎng)項(xiàng)目終止時(shí)儲(chǔ)能設(shè)備的報(bào)廢成本，退役階段所產(chǎn)生的處置成本為：

C3=mKdcbafba+nKdcucfuc

(5)

式中：Kdcba為蓄電池處置成本系數(shù)；Kdcuc為超級(jí)電容處置成本系數(shù)。

綜上，本文中混合儲(chǔ)能的綜合成本C可表示為：

C=C1+C2+C3

(6)

3 配電網(wǎng)儲(chǔ)能系統(tǒng)配置

3.1 網(wǎng)損指標(biāo)

電網(wǎng)網(wǎng)損作為配電網(wǎng)經(jīng)濟(jì)運(yùn)行的重要指標(biāo)，儲(chǔ)能系統(tǒng)接入電網(wǎng)后，因此降低電力系統(tǒng)有功損耗，從而提高經(jīng)濟(jì)運(yùn)行指標(biāo)。圖1表示配電網(wǎng)中接入儲(chǔ)能系統(tǒng)后的簡(jiǎn)易饋線模型，其中Pm+jQm表示負(fù)荷吸收功率，Pm表示有功負(fù)荷，Qm表示無(wú)功負(fù)荷。

圖1 儲(chǔ)能系統(tǒng)接入的饋線模型

為改善潮流分布，有效控制系統(tǒng)的網(wǎng)損，建立系統(tǒng)網(wǎng)損最小的目標(biāo)函數(shù)，即儲(chǔ)能裝置安裝前后網(wǎng)損變化最大，其數(shù)學(xué)表達(dá)式為：

(7)

式中：PLoss為接入儲(chǔ)能裝置前后系統(tǒng)網(wǎng)損的變化量,PLoss>0時(shí)表明接入儲(chǔ)能裝置后總網(wǎng)損降低,PLoss<0時(shí)表明接入儲(chǔ)能裝置后總網(wǎng)損增加。由上式分析可知，系統(tǒng)網(wǎng)損與儲(chǔ)能容量大小及接入位置均有關(guān)。

3.2 電壓偏移指標(biāo)

配電網(wǎng)系統(tǒng)電壓穩(wěn)定性以降低電壓偏移作為優(yōu)化目標(biāo)為f2，以日均電壓偏差為基準(zhǔn)，一天中取樣間隔時(shí)間為1h，計(jì)算每個(gè)節(jié)點(diǎn)電壓偏差平方和，電壓指標(biāo)函數(shù)計(jì)算模型為：

(8)

式中：NL為電力系統(tǒng)中負(fù)荷節(jié)點(diǎn)數(shù)；Ui和Ui*分別為第i節(jié)點(diǎn)的測(cè)量電壓和標(biāo)準(zhǔn)參考電壓，一般取Ui*=1；ΔUi,max為第i節(jié)點(diǎn)處電壓的最大允許偏差值，即允許范圍內(nèi)最大值與最小值之差，ΔUi,max=Ui,max-Ui,min。

3.3 混合儲(chǔ)能多目標(biāo)優(yōu)化配置

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

綜合考慮成本和電壓偏移這兩個(gè)指標(biāo)，混合儲(chǔ)能優(yōu)化配置的目標(biāo)函數(shù)如下：

minF1=C-f1

(9)

minF2=f2

(10)

式中，F(xiàn)1是以綜合成本為優(yōu)化的子目標(biāo)一，F(xiàn)2是以電壓偏移為優(yōu)化的子目標(biāo)二。

3.3.2 約束條件

(1)儲(chǔ)能系統(tǒng)相關(guān)約束

相較于功率型儲(chǔ)能器件，蓄電池的循環(huán)充放電次數(shù)較低，因此為延長(zhǎng)其循環(huán)使用壽命從而降低更換成本，需要考慮避免對(duì)蓄電池過(guò)度充放電，因此對(duì)蓄電池充放電功率和能量建立約束條件，具體公式如下：

充放電功率約束：

(11)

充放電狀態(tài)轉(zhuǎn)換約束：

(12)

初始能量約束：

(13)

儲(chǔ)能容量約束：

(14)

相較于蓄電池而言，超級(jí)電容功率密度高且循環(huán)壽命長(zhǎng)，因此不考慮放電對(duì)其使用壽命的影響，超級(jí)電容組運(yùn)行約束如下：

(15)

其中，

(16)

(17)

式中：Ecmin、Ecmax分別為超級(jí)電容器組的最小和最大儲(chǔ)能容量；Pc,t為超級(jí)電容在t時(shí)刻的輸出功率；ucmax、ucmin分別為超級(jí)電容的最大、最小工作電壓值。

(2)DG運(yùn)行約束

DG電源發(fā)電機(jī)組的輸出功率大小主要與設(shè)備額定功率和自然環(huán)境因素決定，因此需對(duì)風(fēng)光電源的輸出功率做約束，約束表達(dá)式為：

Pimin≤Pi,t≤Pimax

(18)

式中：Pimin為配電網(wǎng)中各類型分布式電源的最小輸出功率；Pimax為配網(wǎng)中各類型分布式電源的最大輸出功率。

對(duì)儲(chǔ)能容量配置時(shí)，需要考慮對(duì)DG中風(fēng)電和光伏出力的消納率，即在優(yōu)化配置后降低棄風(fēng)棄光率，約束表達(dá)式為：

(19)

式中：rPV為優(yōu)化配置后分布式電源中光伏的實(shí)際棄光率；rPV,max為設(shè)置的最大允許棄光率；rWT為優(yōu)化配置后風(fēng)力發(fā)電機(jī)組的實(shí)際棄光率；rWT,max為風(fēng)電機(jī)組最大棄風(fēng)率。

(3)功率平衡約束

分布式電源與負(fù)荷間的正相關(guān)性較弱，通過(guò)儲(chǔ)能設(shè)備進(jìn)行削峰填谷，從而減少棄負(fù)荷概率，因此要求在任一時(shí)刻，儲(chǔ)能系統(tǒng)的剩余容量均大于等于分布式電源發(fā)電量與負(fù)荷需求量的差值，約束表達(dá)式為：

NESSESmax≥η(Pw,t+PV,t-PL,t)·Δt+Nes·PES,t·Δt

(20)

式中：NES為實(shí)際配備的儲(chǔ)能設(shè)備數(shù)量；Pw,t、Pv,t分別為t時(shí)刻，風(fēng)電機(jī)組和光伏機(jī)組的輸出總功率；PL,t為t時(shí)刻負(fù)荷的功率需求；PES,t為t時(shí)刻儲(chǔ)能系統(tǒng)的充電功率；SESmax為混合儲(chǔ)能系統(tǒng)的最大容量。

4 基于小生境的多目標(biāo)改進(jìn)粒子群算法

多目標(biāo)粒子群優(yōu)化(Multi-objective particle swarm Optimization，MOPSO)算法是一種基于鳥類群體智能搜索的隨機(jī)演化計(jì)算方法，對(duì)求解大規(guī)模優(yōu)化問(wèn)題有很快的收斂速度和全局尋優(yōu)能力。

4.1 小生鏡技術(shù)原理

小生境技術(shù)源于自然界中小生境的概念，是指在某種特定的地域中，存在一群相同或相似的物種，各物種間存在交流和競(jìng)爭(zhēng)關(guān)系，通過(guò)“優(yōu)勝劣汰”的生存機(jī)制，淘汰適應(yīng)能力較弱的物種，將生存能力強(qiáng)的物種保留，從而使物種得到進(jìn)化。

粒子間的距離可表示為：

(21)

式中，dij為粒子i與粒子j間的距離；xi、xj分別為D維空間中第i個(gè)和第j個(gè)粒子；N為小生境內(nèi)全部粒子數(shù)。

小生境半徑可表示為：

di=min(‖xi-xj‖)(i,j=1,2,…,m)

(22)

(23)

式中，di為粒子i的最小歐式距離；m為種群內(nèi)粒子數(shù)；C為初始值常數(shù)，m=1時(shí)，小生境半徑取常數(shù)C，否則取di的平均值；當(dāng)粒子間距離小于小生境半徑，即dij

粒子當(dāng)前適應(yīng)度可表示為：

(24)

式中，fsh(dij)代表共享函數(shù)，當(dāng)粒子間距離dij小于小生鏡半徑Rch時(shí)，則共享函數(shù)的計(jì)算公式為：

(25)

式中，λ為控制共享函數(shù)形狀的參數(shù)。

4.2 小生境多目標(biāo)粒子群優(yōu)化算法的計(jì)算流程

4.2.1 算法流程

小生境多目標(biāo)粒子群優(yōu)化算法(Niche Multi-objective Particle Swarm Optimization,NMOPSO)在原有算法的基礎(chǔ)上，引入小生境處理機(jī)制，將這種多樣性和自我調(diào)節(jié)能力特性引入到MOPSO算法中，提出了基于小生境技術(shù)的多目標(biāo)粒子群算法對(duì)問(wèn)題求解。利用混沌變異技術(shù)，有效保證粒子全局搜索，避免算法后期的“早熟”問(wèn)題。小生境多目標(biāo)粒子群算法步驟如下：

(1)參數(shù)初始化：根據(jù)控制變量與約束條件隨機(jī)生成初始種群，設(shè)置種群規(guī)模與最大迭代次數(shù)。計(jì)算小生境半徑和個(gè)體適應(yīng)度，并生成小生境種群。計(jì)算粒子適應(yīng)度時(shí)，引入小生境競(jìng)爭(zhēng)機(jī)制保證了種群的多樣性?；谛∩掣?jìng)爭(zhēng)機(jī)制的概率選擇公式為：

(26)

式中，xi和Si表示第i個(gè)粒子及其適應(yīng)度函數(shù)，粒子i的相似度越高，則該粒子被選中的概率越小，通過(guò)概率機(jī)制使適應(yīng)度值低的個(gè)體能繼續(xù)參與進(jìn)化，因此提高了種群多樣性。

(2)更新粒子的速度和位置，將外部檔案集初始化，并定義為全局最優(yōu)解集。

(27)

式中，粒子標(biāo)號(hào)可表示成i=1,2,3,…,m；k為迭代代數(shù)；c1、c2代表學(xué)習(xí)因子，默認(rèn)取值范圍為(1.4,2.1)；r1、r2為介于0到1的隨機(jī)數(shù)；P表示粒子群的位置；V表示迭代速度。

針對(duì)慣性權(quán)重問(wèn)題，傳統(tǒng)線性遞減的慣性權(quán)重策略無(wú)法同時(shí)兼顧全局與局部尋優(yōu)能力，本文采用一種新的非線性慣性權(quán)重遞減策略，ω隨迭代次數(shù)增加呈非線性遞減：

ω=ωstart-k(ωstart-ωend)sin(πt/2tmax)

(28)

式中：ωstart和ωend分別為起始和終止慣性權(quán)重；k為調(diào)節(jié)系數(shù)，與當(dāng)前迭代次數(shù)t有關(guān)，當(dāng)t=0，k=1時(shí)，此時(shí)ω為初始權(quán)重值，全局搜索能力較強(qiáng)，不易陷入局部最優(yōu)。ω隨迭代次數(shù)增加呈非線性遞減趨勢(shì)，使算法由全局搜索向局部搜過(guò)過(guò)渡。當(dāng)達(dá)到最大迭代次數(shù)tmax時(shí)，調(diào)節(jié)系數(shù)k較大，ω接近終止慣性權(quán)重值，增強(qiáng)了局部搜索能力，有效解決算法后期的“早熟”問(wèn)題。

(3)計(jì)算適應(yīng)度函數(shù)并更新外部檔案，判斷新產(chǎn)生粒子i與外部檔案中粒子是否存在支配關(guān)系，若粒子i與其他粒子間非支配，則將粒子i添加到外部檔案中，否則繼續(xù)迭代。

(4)利用小生境的選擇機(jī)制，在不同生境內(nèi)淘汰個(gè)體適應(yīng)度較差的粒子，保證Pareto前沿解始終最優(yōu)。

(5)采用輪盤賭注策略，隨機(jī)選取N個(gè)粒子變異為新的群體。

(6)對(duì)外部檔案的優(yōu)良個(gè)體進(jìn)行選擇，選取新的最優(yōu)解并繼續(xù)迭代，產(chǎn)生新種群。

(7)計(jì)算適應(yīng)度函數(shù)，獲取全局最優(yōu)值。

(8)判斷程序是否滿足終止，當(dāng)達(dá)到最大迭代次數(shù)時(shí)循環(huán)結(jié)束并輸出最優(yōu)解集，否則返回步驟3。

小生境多目標(biāo)粒子群優(yōu)化算法流程圖如圖2所示。

圖2 小生境多目標(biāo)粒子群優(yōu)化算法流程圖

5 算例仿真與分析

為驗(yàn)證本文所建混合儲(chǔ)能優(yōu)化配置模型的有效性，以修改后的IEEE33配網(wǎng)為例進(jìn)行仿真分析，采用改進(jìn)后的多目標(biāo)粒子群算法進(jìn)行求解。分析儲(chǔ)能接入的位置與容量對(duì)系統(tǒng)電壓穩(wěn)定性、網(wǎng)損和綜合費(fèi)用的影響。

5.1 基礎(chǔ)參數(shù)

以IEEE33節(jié)點(diǎn)標(biāo)準(zhǔn)算例進(jìn)行混合儲(chǔ)能的混合配置計(jì)算，該配電網(wǎng)絡(luò)有32條支路，系統(tǒng)電壓為12.66kV，該網(wǎng)絡(luò)結(jié)構(gòu)圖如圖3所示。儲(chǔ)能待選類型有超級(jí)電容和蓄電池，儲(chǔ)能可接入配電網(wǎng)任意位置，其參數(shù)如表1所示。

圖3 IEEE33節(jié)點(diǎn)系統(tǒng)連接圖

表1 待選儲(chǔ)能基本參數(shù)

本文研究的混合儲(chǔ)能在變電站10kV母線側(cè)接入，作為負(fù)荷的供電補(bǔ)充，減少負(fù)荷的波動(dòng)對(duì)配電網(wǎng)設(shè)備動(dòng)作的影響，調(diào)整網(wǎng)絡(luò)有功潮流的波動(dòng)，減少有功功率從變電站節(jié)點(diǎn)向配電網(wǎng)末端流動(dòng)，從而減少負(fù)荷損失和電壓波動(dòng)。

圖4 儲(chǔ)能接入配電網(wǎng)方式

風(fēng)電光伏分別接入節(jié)點(diǎn)8和20，其預(yù)測(cè)出力如圖5所示。負(fù)荷預(yù)測(cè)值如圖6所示。

圖5 風(fēng)光預(yù)測(cè)出力

圖6 電網(wǎng)負(fù)荷出力

5.2 仿真結(jié)果分析

5.2.1 MOPSO與NMOPSO仿真結(jié)果

基于小生境技術(shù)的改進(jìn)粒子群算法能最大程度的保證粒子的豐富性以及種群的多樣性，避免的算法過(guò)早收斂而陷入局部最優(yōu)。從圖7的計(jì)算顯示，小生境粒子群算法(NMOPSO)的優(yōu)化效果優(yōu)于多目標(biāo)粒子群算法(MOPSO)。

圖7 MOPSO與NMOPSO的優(yōu)化結(jié)果

5.3 混合儲(chǔ)能優(yōu)化配置分析

在系統(tǒng)接入兩套混合儲(chǔ)能裝置后，一方面可以在主變低壓側(cè)實(shí)現(xiàn)對(duì)系統(tǒng)的有功潮流的改善，避免潮流從變壓器高壓側(cè)往低壓側(cè)流動(dòng)而造成網(wǎng)損的增加，同時(shí)降低電網(wǎng)電壓的波動(dòng)。

由圖8可以看出，在節(jié)點(diǎn)6和節(jié)點(diǎn)30接入兩套混合儲(chǔ)能裝置后，可以很大的程度改善電網(wǎng)網(wǎng)損。圖10的仿真結(jié)果表明，加入儲(chǔ)能裝置后電網(wǎng)網(wǎng)損降低了大約0.05MW。

圖8 配網(wǎng)加入混合儲(chǔ)能后網(wǎng)損對(duì)比

圖9 配網(wǎng)加入混合儲(chǔ)能后電壓對(duì)比

圖10 配網(wǎng)加入混合儲(chǔ)能后電網(wǎng)電壓穩(wěn)定性LCPI指標(biāo)

根據(jù)圖9的仿真結(jié)果表明，加入混合儲(chǔ)能后對(duì)不同時(shí)間段均提升了系統(tǒng)最低電壓值。在節(jié)點(diǎn)6和節(jié)點(diǎn)30并入兩套混合儲(chǔ)能裝置后，可以很大的程度改善電網(wǎng)電壓。從數(shù)量關(guān)系上可以看出，加入儲(chǔ)能裝置后最低電壓值提升了0.01pu左右，對(duì)于提升電網(wǎng)的電壓有一定的提升。

LCPI代表電壓穩(wěn)定指標(biāo)[14]，該值越大代表電壓穩(wěn)定性越好。由圖10仿真結(jié)果表明，配置混合儲(chǔ)能后對(duì)于提升配電網(wǎng)電壓穩(wěn)定性效果明顯，這是因?yàn)閮?chǔ)能參與到電壓調(diào)節(jié)與優(yōu)化模型中，通過(guò)儲(chǔ)能的充放電策略來(lái)平滑風(fēng)光出力，減少電壓波動(dòng)，實(shí)現(xiàn)降低網(wǎng)絡(luò)損耗，最終使得電網(wǎng)的電壓穩(wěn)定性得到提高。

6 結(jié)論

本文綜合考慮分布式電源、儲(chǔ)能裝置及配網(wǎng)安全運(yùn)行約束條件，建立以儲(chǔ)能成本、網(wǎng)損率以及電壓穩(wěn)定性作為優(yōu)化目標(biāo)的配網(wǎng)儲(chǔ)能優(yōu)化配置模型，并以IEEE33節(jié)點(diǎn)配網(wǎng)為算例，運(yùn)用小生境多目標(biāo)粒子群優(yōu)化算法NMOPSO對(duì)該模型進(jìn)行求解，驗(yàn)證所提模型及算法的有效性。結(jié)果表明小生境多目標(biāo)粒子群優(yōu)化算法適合求解多目標(biāo)儲(chǔ)能優(yōu)化配置模型，可兼顧經(jīng)濟(jì)性和安全性指標(biāo)。