基于多源協(xié)同的主動(dòng)配電網(wǎng)優(yōu)化控制策略

李慶波， 羅 銳， 翟旭京， 辛超山， 邊家瑜

（國(guó)網(wǎng)新疆電力有限公司 經(jīng)濟(jì)技術(shù)研究院， 烏魯木齊 830000）

隨著我國(guó)社會(huì)經(jīng)濟(jì)的發(fā)展，電力需求持續(xù)增長(zhǎng)，清潔低碳成為我國(guó)能源轉(zhuǎn)型發(fā)展的重要特征.如何在保證電力系統(tǒng)安全、穩(wěn)定運(yùn)行的同時(shí)，最大程度地消納各種清潔能源成為當(dāng)前電網(wǎng)轉(zhuǎn)型發(fā)展亟待解決的問(wèn)題[1].與傳統(tǒng)的被動(dòng)式運(yùn)營(yíng)管理配電網(wǎng)相比，主動(dòng)配電網(wǎng)可以充分利用各種先進(jìn)信息與管理技術(shù)，實(shí)現(xiàn)對(duì)分布式能源靈活有效的控制[2-3].

目前關(guān)于主動(dòng)配電網(wǎng)的優(yōu)化控制研究多是基于能量管理系統(tǒng)的角度，較少涉及主動(dòng)配電網(wǎng)供需兩側(cè)的組合優(yōu)化問(wèn)題[12-13].本文針對(duì)主動(dòng)配電網(wǎng)的多源協(xié)同問(wèn)題，充分考慮可再生能源以及負(fù)荷需求，構(gòu)建了以節(jié)約系統(tǒng)調(diào)度成本、提升可再生能源利用率及提高用戶(hù)滿(mǎn)意度為目標(biāo)的多源協(xié)同主動(dòng)配電網(wǎng)優(yōu)化調(diào)度模型.利用改進(jìn)微分進(jìn)化算法對(duì)模型進(jìn)行求解，實(shí)現(xiàn)了對(duì)各類(lèi)資源的有效利用.

1 多源協(xié)同優(yōu)化模型

1.1 配電網(wǎng)模型

主動(dòng)配電網(wǎng)中分布式電源的出力預(yù)測(cè)與負(fù)荷預(yù)測(cè)是實(shí)現(xiàn)對(duì)其優(yōu)化控制的基礎(chǔ)，因此首先建立了包含風(fēng)力發(fā)電、光伏發(fā)電、電動(dòng)汽車(chē)動(dòng)態(tài)充放電等多種電源的配電網(wǎng)模型.

由統(tǒng)計(jì)數(shù)據(jù)可知，風(fēng)力發(fā)電機(jī)的輸出功率與風(fēng)速之間關(guān)系式[14-15]為

（1）

式中：PWT為風(fēng)力發(fā)電機(jī)的實(shí)際輸出功率；PWR為風(fēng)機(jī)的最大功率；v、vci、vco、vr分別為風(fēng)機(jī)所處位置的風(fēng)速、切入風(fēng)速、切出風(fēng)速與額定風(fēng)速.

光伏發(fā)電輸出功率表達(dá)式[16]為

PM=ARη

（2）

式中：PM為光伏發(fā)電輸出功率；R為光照強(qiáng)度；A為方陣面積；η為光電轉(zhuǎn)換效率.

電動(dòng)汽車(chē)的日行駛里程md服從對(duì)數(shù)正態(tài)分布，其概率密度函數(shù)表達(dá)式為

（3）

電動(dòng)汽車(chē)充電時(shí)間與離網(wǎng)時(shí)間的概率密度函數(shù)表達(dá)式為

（4）

（5）

式中，t0、t1分別為電動(dòng)汽車(chē)接入電網(wǎng)時(shí)間與脫離電網(wǎng)時(shí)間.計(jì)算過(guò)程中，通過(guò)計(jì)算機(jī)產(chǎn)生隨機(jī)數(shù)的方法模擬電動(dòng)汽車(chē)的充放電行為.

負(fù)荷概率密度函數(shù)可用正態(tài)分布表示為

（6）

式中：PL為負(fù)荷功率；μL為負(fù)荷功率期望；σL為負(fù)荷功率標(biāo)準(zhǔn)差.則負(fù)荷的累計(jì)分布函數(shù)為

（7）

1.2 優(yōu)化調(diào)度模型

文中從配網(wǎng)調(diào)度的經(jīng)濟(jì)性、清潔能源的利用率及用戶(hù)滿(mǎn)意度三個(gè)方面，建立了多源協(xié)同的主動(dòng)配電網(wǎng)優(yōu)化控制模型，在滿(mǎn)足系統(tǒng)功率平衡約束條件下，采用改進(jìn)的微分進(jìn)化算法對(duì)其進(jìn)行求解.

1） 調(diào)度的經(jīng)濟(jì)性.主動(dòng)配電網(wǎng)運(yùn)行總成本主要由調(diào)度成本C1、購(gòu)電成本C2、分布式電源運(yùn)維成本C3，清潔能源發(fā)電補(bǔ)貼C4與電動(dòng)汽車(chē)充電收入C5構(gòu)成，其表達(dá)式為

C0=C1+C2+C3-（C4+C5）

（8）

除去成本外，配電網(wǎng)向用戶(hù)售電的總收益F表示為

（9）

式中：PL，t為主動(dòng)配電網(wǎng)在t時(shí)刻的負(fù)荷；CEX，t為當(dāng)前時(shí)刻的電價(jià)；Sum為統(tǒng)計(jì)周期的總數(shù).則關(guān)于調(diào)度成本的目標(biāo)函數(shù)表達(dá)式為

（10）

2） 清潔能源利用率.定義清潔能源利用率為調(diào)度周期內(nèi)清潔能源的消納量占其發(fā)電出力的比例，則清潔能源利用率表達(dá)式為

（11）

式中，PRES，t為t時(shí)刻清潔能源的總發(fā)電功率.則目標(biāo)函數(shù)可表示為

maxf2=γ

（12）

3） 用戶(hù)滿(mǎn)意度.用戶(hù)滿(mǎn)意度屬于用戶(hù)的一種用電體驗(yàn)，為對(duì)其進(jìn)行量化衡量，本文從用戶(hù)用電時(shí)間β1、用電量β2及電動(dòng)汽車(chē)的荷電狀態(tài)β3構(gòu)建的滿(mǎn)意度函數(shù)為

maxf3=β1+β2+β3

（13）

其中，β1、β2、β3分別表示為

（14）

式中：S1、S2分別為削減負(fù)荷、可轉(zhuǎn)移負(fù)荷參與調(diào)度的時(shí)間段；ΔQD，t、QDT，t為可削減負(fù)荷在t時(shí)刻的削減量和總用電量；ΔQS，t、QST，t為可轉(zhuǎn)移負(fù)荷在t時(shí)刻內(nèi)的轉(zhuǎn)移量與總用電量.

（15）

式中：Us為用戶(hù)總數(shù)；ΔTD，y、TDT，y為可削減負(fù)荷用戶(hù)y的負(fù)荷削減時(shí)間和可削減總時(shí)間；ΔTS，y、TST，y為可轉(zhuǎn)移負(fù)荷用戶(hù)y的負(fù)荷轉(zhuǎn)移時(shí)間與可轉(zhuǎn)移總時(shí)間.

（16）

式中：NEV為電動(dòng)汽車(chē)總數(shù)；SocE，y、SocB，y分別為汽車(chē)現(xiàn)有電量與滿(mǎn)電量.

綜合目標(biāo)函數(shù)f1、f2、f3，則總的目標(biāo)函數(shù)可表示為

（17）

1.3 微分進(jìn)化算法

本文建立的優(yōu)化調(diào)度模型屬于非線性規(guī)劃問(wèn)題，因此采用微分進(jìn)化算法對(duì)其進(jìn)行求解.

1） 初始化.假設(shè)共有N個(gè)種群，自變量的個(gè)數(shù)為D，第n代的種群為Xn.

（18）

潮州電廠引進(jìn)的1000MW#3、#4超超臨界機(jī)組，分別于2009年11月9日、2010年1月18日通過(guò)168小時(shí)試運(yùn)行，于2010年10月18 -22日進(jìn)行了現(xiàn)場(chǎng)監(jiān)測(cè)及調(diào)查。電廠廢水為循環(huán)利用，實(shí)現(xiàn)零排放；大氣顆粒污染物主要采用靜電除塵器進(jìn)行處理，脫硫采用石灰石-石膏濕法煙氣脫硫工藝，采用低NOx燃燒器、分級(jí)配風(fēng)、降低燃燒溫度水平等方式控制NOx產(chǎn)生；固體廢物由專(zhuān)門(mén)的公司負(fù)責(zé)清理。

3） 交叉操作.交叉操作通過(guò)變異量與目標(biāo)量之間的隨機(jī)重組，將變異量引入自變量的解集中.交叉過(guò)程表示為

（19）

式中：i=1，2，…，N；j=1，2，…，D；φ為0～1之間的隨機(jī)數(shù)；cR為交叉因子，用于指定下一代種群的變化率；r為交叉概率，是一個(gè)在0～D中的隨機(jī)整數(shù).

4） 選擇操作.假設(shè)待求優(yōu)化函數(shù)為minf（x），則選擇操作表達(dá)式為

（20）

根據(jù)實(shí)際運(yùn)行經(jīng)驗(yàn)，本文在計(jì)算過(guò)程中選擇種群個(gè)數(shù)為500，縮放因子為0.6，cR取0.7.

5） 算法改進(jìn).微分進(jìn)化算法在進(jìn)化過(guò)程中，變異與交叉過(guò)程是隨機(jī)的，在尋找最優(yōu)解的過(guò)程中會(huì)導(dǎo)致收斂速度較慢.因此本文在保持種群多樣性的同時(shí)，通過(guò)將種群的最優(yōu)個(gè)體引入變異量中，以提高算法的收斂速度，算法改進(jìn)表達(dá)式為

（21）

變異量中的縮放因子與交叉因子用于控制算法的搜索速度與范圍.尋優(yōu)過(guò)程中，較大的縮放因子有利于完成全局快速搜索，但不利用算法的局部尋優(yōu).本文對(duì)縮放因子進(jìn)行動(dòng)態(tài)控制，令其值在0.5～0.75范圍內(nèi)動(dòng)態(tài)變化，控制表達(dá)式為

（22）

交叉因子用于控制算法的局部搜索能力，交叉因子越小局部搜索能力越好，在初期可采用較小的交叉因子以提高局部搜索能力；在后期通過(guò)增大交叉因子防止算法陷入局部最優(yōu).本文對(duì)交叉因子進(jìn)行動(dòng)態(tài)控制，令其值在0.3～0.6內(nèi)變化，控制表達(dá)式為

（23）

2 算例分析

文中以IEEE 11節(jié)點(diǎn)網(wǎng)絡(luò)為例，對(duì)算法進(jìn)行驗(yàn)證.圖1為電網(wǎng)結(jié)構(gòu)，表1為線路基本參數(shù)，其值為標(biāo)幺值.

圖1中，L1～L5為配電網(wǎng)負(fù)荷.其中固定負(fù)荷所占比例為75%，可轉(zhuǎn)移負(fù)荷與可削減負(fù)荷所占比例分別為15%和10%.各負(fù)荷用戶(hù)參與轉(zhuǎn)移或削減調(diào)度過(guò)程中實(shí)行階梯型補(bǔ)償機(jī)制.補(bǔ)償費(fèi)用如表2所示.

假設(shè)該配電網(wǎng)所有負(fù)荷的用戶(hù)數(shù)為20，某一時(shí)刻參與負(fù)荷削減調(diào)度的用戶(hù)數(shù)量服從均值為12，方差為1.25的正態(tài)分布，負(fù)荷削減量服從均值為10，方差為2.5的正態(tài)分布；某一時(shí)刻參與轉(zhuǎn)移負(fù)荷調(diào)度的用戶(hù)數(shù)量服從均值為12，方差為1.4的正態(tài)分布，轉(zhuǎn)移負(fù)荷量服從均值為6，方差為2.5的正態(tài)分布，轉(zhuǎn)移負(fù)荷目標(biāo)時(shí)間服從均值為3，方差為1.2的正態(tài)分布.

圖1 配電網(wǎng)結(jié)構(gòu)Fig.1 Distribution network structure

表1 配電網(wǎng)線路阻抗參數(shù)Tab.1 Line impedance parameters of distribution network

表2 負(fù)荷階梯型補(bǔ)償費(fèi)用Tab.2 Compensation cost of load stages

參與調(diào)度的電動(dòng)汽車(chē)為30輛，補(bǔ)償價(jià)格為0.6元/kW·h，每臺(tái)電動(dòng)汽車(chē)的充電容量為80 kW·h，額定充電功率為10 kW，期望行駛里程500 km.假設(shè)電動(dòng)汽車(chē)每天充、放電各一次，開(kāi)始充電時(shí)間服從均值為18，方差為3.5的正態(tài)分布；結(jié)束充電時(shí)間服從均值為8，方差為3.5的正態(tài)分布.

電網(wǎng)電價(jià)采用三費(fèi)率分時(shí)電價(jià).峰時(shí)段（10∶00～15∶00；18∶00～21∶00）電價(jià)為0.8元/kW·h；平時(shí)段（7∶00～10∶00；15∶00～18∶00；21∶00～23∶00）電價(jià)為0.45元/kW·h；低谷時(shí)段（23∶00～7∶00）電價(jià)為0.25元/kW·h.

利用本文提出的風(fēng)力、光伏及負(fù)荷預(yù)測(cè)模型得到可再生能源日內(nèi)典型出力與負(fù)荷曲線如圖2所示.

圖2 可再生能源出力及負(fù)荷預(yù)測(cè)模型Fig.2 Prediction model of renewable energy output and loads

為驗(yàn)證調(diào)度策略的有效性，本文采用靜態(tài)優(yōu)化調(diào)度與動(dòng)態(tài)優(yōu)化調(diào)度方法，即常規(guī)微分進(jìn)化算法與改進(jìn)的微分進(jìn)化算法對(duì)配電網(wǎng)進(jìn)行優(yōu)化控制，仿真結(jié)果如圖3所示.

圖3 優(yōu)化控制仿真結(jié)果Fig.3 Simulation results of optimal control

由圖3可知，通過(guò)對(duì)供給側(cè)與需求側(cè)資源的優(yōu)化調(diào)度，負(fù)荷曲線相比于優(yōu)化前有明顯的改善.優(yōu)化前負(fù)荷曲線的極差為615.6，標(biāo)準(zhǔn)差為175.5；采用靜態(tài)優(yōu)化后的負(fù)荷曲線極差為238.6，標(biāo)準(zhǔn)差為59.0；采用動(dòng)態(tài)優(yōu)化后的負(fù)荷曲線極差為195.8，標(biāo)準(zhǔn)差為51.8.由此可見(jiàn)，基于本文所提優(yōu)化算法得出的負(fù)荷曲線具有較為明顯的削峰填谷特性，且動(dòng)態(tài)優(yōu)化的效果優(yōu)于靜態(tài)優(yōu)化.

圖4為兩種優(yōu)化算法的結(jié)果比較，由圖4可知，相比于優(yōu)化前，靜態(tài)優(yōu)化與動(dòng)態(tài)優(yōu)化模式下分別節(jié)約調(diào)度成本15.3%、18.9%.優(yōu)化前為了滿(mǎn)足用戶(hù)的用電需求，不僅需要額外的調(diào)度成本，同時(shí)存在大量的棄風(fēng)棄電現(xiàn)象.優(yōu)化后，可再生能源的利用率由優(yōu)化前的85%提高到優(yōu)化后的90%以上，動(dòng)態(tài)優(yōu)化清潔能源利用率較靜態(tài)優(yōu)化提高3.12%.通過(guò)對(duì)主動(dòng)配電網(wǎng)的優(yōu)化控制，減小了棄風(fēng)棄光成本與購(gòu)電成本，具有顯著的環(huán)保效益與經(jīng)濟(jì)效益.

圖4 兩種模式參數(shù)比較Fig.4 Comparison of parameters between two modes

圖5為優(yōu)化后配網(wǎng)各微電源的輸出功率.由圖5可知，光伏出力主要集中在6∶00～16∶00區(qū)間內(nèi)，基本為滿(mǎn)負(fù)荷出力.風(fēng)力發(fā)電機(jī)輸出功率主要集中在18點(diǎn)后，覆蓋了晚間用電的高峰時(shí)段，充分實(shí)現(xiàn)了清潔能源的有效利用.18∶00～22∶00為一天的用電高峰期，僅依靠風(fēng)力發(fā)電難以滿(mǎn)足負(fù)荷需求，此時(shí)也有部分儲(chǔ)能電源開(kāi)始向電網(wǎng)輸出功率.此外，按照本文的調(diào)度策略以及負(fù)荷的削減與轉(zhuǎn)移原則，23∶00～10∶00為用電低谷與平時(shí)段，電動(dòng)汽車(chē)和儲(chǔ)能電池陸續(xù)充電，既可以節(jié)約成本，同時(shí)又消納了多余的風(fēng)電.

圖5 各電源輸出功率Fig.5 Output power of each power supply

綜上所述，采用基于微分進(jìn)化算法的多源協(xié)同主動(dòng)配電網(wǎng)優(yōu)化控制策略能夠在不降低用戶(hù)滿(mǎn)意度的前提下，有效節(jié)約調(diào)度成本，達(dá)到削峰填谷目的，同時(shí)有力促進(jìn)了清潔能源的消納及綜合能源的利用效率.

3 結(jié) 論

本文研究了一種基于微分進(jìn)化算法的多源協(xié)同主動(dòng)配電網(wǎng)優(yōu)化控制策略.在微分進(jìn)化算法的基礎(chǔ)上，通過(guò)引入動(dòng)態(tài)縮放因子與交叉因子，對(duì)多源協(xié)同的主動(dòng)配電網(wǎng)進(jìn)行負(fù)荷預(yù)測(cè)和控制.仿真結(jié)果表明，采用動(dòng)態(tài)優(yōu)化比傳統(tǒng)微分進(jìn)化算法節(jié)約成本3.6%，用戶(hù)滿(mǎn)意度提高4.11%，清潔能源利用率提高3.12%.以此為基礎(chǔ)，通過(guò)協(xié)調(diào)需求側(cè)與供給側(cè)資源的互動(dòng)機(jī)制，建立了多源協(xié)同的主動(dòng)配電網(wǎng)優(yōu)化控制模型.優(yōu)化后配網(wǎng)調(diào)度成本可節(jié)約15%以上，清潔能源利用率可達(dá)90%以上.