徐 虹 蘆晶晶 孫宇斌 艾 欣 張建華

（1. 華北電力大學(xué)新能源電力系統(tǒng)國家重點實驗室，北京 102206；2.中國電力科學(xué)研究院，北京 100192）

近年來，人們對經(jīng)濟、清潔、高效的分布式電源的需求不斷提高，促使傳統(tǒng)電網(wǎng)急需向含有大量分布式電源的智能電網(wǎng)轉(zhuǎn)變[1]。但是，各種分布式電源的特性不同，運行方式相差較大，接入后會給配電網(wǎng)的規(guī)劃、預(yù)測和運行帶來很大的不確定性[2-3]。如柴油機發(fā)電雖然成本低，但污染嚴重；風(fēng)力發(fā)電和光伏發(fā)電無污染，但一次設(shè)備成本高，易受天氣影響，功率輸出不易控制；燃料電池污染小、易控制，但反應(yīng)速度慢、成本高；儲能設(shè)備雖然反應(yīng)速度快，但不易與其他電源協(xié)調(diào)運行[4]。因此對配電網(wǎng)內(nèi)各個分布式電源的優(yōu)化調(diào)度是十分必要的，還需要同時考慮經(jīng)濟性和環(huán)保因素，也就是說是一個多目標(biāo)尋優(yōu)問題[5]。

本文綜合考慮運行費用和環(huán)境效益，建立了含柴油機、風(fēng)力發(fā)電機、光伏發(fā)電單元和燃料電池等分布式電源的智能配電網(wǎng)優(yōu)化調(diào)度模型，并在滿足各種約束條件的基礎(chǔ)上，使用遺傳算法求解了該模型。

1 多目標(biāo)優(yōu)化調(diào)度模型

對于含分布式電源的配電網(wǎng)優(yōu)化調(diào)度，既要考慮安全經(jīng)濟性又要考慮保因性，因此是一個多目標(biāo)優(yōu)化問題，并且各個子目標(biāo)相互制約、相互矛盾，模型求解過程必須準(zhǔn)確評估各種因素的影響。

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

式中，0≤ω1≤1，0≤ω2≤1為權(quán)系數(shù)，且 w1+w2=1。Cp為最小污染物處理費用，Cy為最小運行費用。該目標(biāo)函數(shù)綜合考慮了運行費用和環(huán)境效益，采用線性加權(quán)法表示，為多目標(biāo)優(yōu)化問題。由于運行費用最小和污染物處理費用最小兩者之間存在一定的矛盾關(guān)系，并且帕雷托最優(yōu)解通常也很多，因而希望優(yōu)化算法能夠給出最優(yōu)解集以供選擇，從而折中考慮或側(cè)重某一方面。

Cp和Cy又可以分別表示為

式中，nDG接入配電網(wǎng)的分布式電源個數(shù) ；m為所排放的污染物類型（SO2、CO2、NOX等）；Cm為處理每 kg 污染物的費用（元/kg）；函數(shù)f為輸出電能P時所排放的污染物（kg/kW·h）。CDGi為第i個分布式電源的固定投資費用（元）；WDGi為第i個分布式電源的檢修、維護費用（元）)；WDN配電網(wǎng)的檢修、維護費用（元）。

1.2 約束條件

1）等式約束條件：

式中，PLt為t時刻系統(tǒng)中的總有功負荷（kW）；PDGt為分布式電源容量；PDNt配電網(wǎng)和可中斷負荷；p為可調(diào)度發(fā)電單元數(shù)目；q為不可調(diào)度單元數(shù)目；Pit可調(diào)度型發(fā)電單元 t時刻的功率輸出（kW）；Pkt不可調(diào)度型發(fā)電單元t時刻的功率輸出（kW）。

2）不等式約束條件：

式（6）為分布式電源功率輸出限值約束，Pitmin是第 i個發(fā)電單元的最小輸出功率；Pitmax是第 i個發(fā)電單元的最大輸出功率。式（7）為儲能單元存儲容量約束，Smin是t時刻的容量下限，Smax是t時刻的容量上限。式（8）為安全裕度約束，Ps為安全裕度預(yù)留功率。

2 模型求解

遺傳算法[6-7]具有效率高、搜索的空間大、魯棒性好等特點，并且由于含分布式電源的智能配電網(wǎng)優(yōu)化調(diào)度模型是多階離散非線性優(yōu)化問題，因此采用遺傳算法[8]。計算流程如圖1所示。

圖1 計算流程圖

為了實現(xiàn)含分布式電源的配電網(wǎng)的綜合效益最大，這里使用隨機加權(quán)法對目標(biāo)函數(shù)進行優(yōu)化，按式（1）計算目標(biāo)函數(shù)時，每一個體在迭代過程中均從區(qū)間[0，1]中產(chǎn)生隨機數(shù)ω1和ω2作為權(quán)系數(shù)，并且 w1+w2=1。通過這種遺傳算法，搜索方向隨機，因此只要群體規(guī)模足夠大，產(chǎn)生的隨機數(shù)足夠均勻，就能得到所有帕雷托最優(yōu)解。本模型為雙目標(biāo)函數(shù)，圖2所示是其遺傳算法的搜索方向。

圖2 遺傳算法的兩個搜索方向

同時采取排擠小生境技術(shù)，這樣就可以維持種群的多樣性[9]，保證帕雷托解的多樣性。遺傳操作過程時，在某代群體中隨機抽取 1/τ（τ為排擠因子，一般取 2～3）的個體，由這些個體組成排擠成員，如果排擠成員和新產(chǎn)生的個體相似性較大（個體編碼之間的歐氏距離較?。?，就會被來排擠掉。歐氏距離的計算式：

若干次迭代運算和排擠過程一并進行，逐步分類種群中個體，形成小生境，從而維持了群體的多樣性。

為了加快算法的效率，需要對初始種群進行處理，在編碼過程中對于問題的決策變量 Pit采用浮點數(shù)編碼方式，優(yōu)先使用分布式電源 Pit=Pimax。群體初始化時，優(yōu)先檢查分布式電源是否滿足約束條件，如滿足，就在可調(diào)度功率區(qū)間[Pitmin，Pitmax]內(nèi)產(chǎn)生一個隨機數(shù)。當(dāng)產(chǎn)生不可行解時，可以在目標(biāo)函數(shù)的基礎(chǔ)上利用加法方式增加懲罰項作為新目標(biāo)函數(shù)，將不可行解的懲罰約束問題轉(zhuǎn)換為無約束問題。然后采用錦標(biāo)賽選擇方式，不需要設(shè)計專門的適應(yīng)度函數(shù)，只需從上一代群體中隨機選擇并挑選出適應(yīng)度最優(yōu)的個體，多次按照此種方式操作，構(gòu)成新群體。接下來在個體編碼中隨機設(shè)置交叉點并隨機選擇基因，采用高斯變異方式，增加高斯變異項[10]。在更新群體時，還可以采用精英保留策略，來保證迭代過程中最優(yōu)解的不退化。

3 算例分析

假設(shè)配電網(wǎng)中電源主要有柴油機（DE）、風(fēng)力發(fā)電（WT）、光伏發(fā)電（PV）、燃料電池（FC）以及配電網(wǎng)常規(guī)機組（DN）等，其相關(guān)參數(shù)見表1。

表1 電源參數(shù)表

假設(shè)依據(jù)風(fēng)速和輻射強度預(yù)測各時段分布式電源發(fā)電量，則根據(jù)目標(biāo)函數(shù)，各電源各時段優(yōu)化調(diào)度出力，如圖3所示。

圖3 各電源優(yōu)化調(diào)度出力

表2還給出了目標(biāo)函數(shù)的計算結(jié)果，驗證了模型及算法可行性。

由表2可知，在目標(biāo)函數(shù)下，運行費用、污染物處理費用以及總費用分別節(jié)省了39.34%、28.31%和37.35%，既體現(xiàn)了經(jīng)濟效益，又體現(xiàn)了經(jīng)濟效益。

表2 目標(biāo)函數(shù)下的費用/元

4 結(jié)論

通過計算可以看出，綜合考慮運行費用和環(huán)境效益的多目標(biāo)優(yōu)化調(diào)度模型具有良好的經(jīng)濟效益和環(huán)境效益。在解決靜態(tài)調(diào)度問題過程中，遺傳算法能夠很好地收斂并且得出精確度較高的結(jié)果，但是，將各種動態(tài)約束條件考慮進去后，該算法的收斂性和計算速度都下降了。因此，以后還應(yīng)該加強多目標(biāo)優(yōu)化算法的理論和實踐研究。

[1]丁明,王敏.分布式發(fā)電技術(shù)[J].電力自動化設(shè)備,2004,24(7):31-36.

[2]PENG F Z.Editorial special issue on distributed power generatio[J]. IEEE Transactions on Power Electronic,2004, 19(5):1157-1158.

[3]鄭劍,劉克,申金平.分布式供電系統(tǒng)對電網(wǎng)的影響[J].東北電力技術(shù),2006(4):4-6.

[4]梁振鋒,楊曉萍,張娉.分布式發(fā)電技術(shù)及其在中國的發(fā)展[J].西北水電,2006(1):51-53.

[5]王民量,張伯明,夏清.考慮多種約束條件的機組組合新算法[J].電力系統(tǒng)自動化,2000, 24(12): 29-35.

[6]陳國良.遺傳算法及其應(yīng)用[M].北京:人民郵電出版社,2001.

[7]王小平, 曹立明. 遺傳算法－理論、應(yīng)用與軟件實現(xiàn)[M].西安交通大學(xué)出版社, 2002.

[8]EI-KHATTAM WALID, HEGAZY, Y G, SALAMA M M A.An integrated distributed generation optimization model for distribution system planning[J].IEEE Trans on Power Systems, 2005, 20 (2):1158-1165.

[9]陳皓勇,王錫凡.機組組合問題的優(yōu)化方法綜述[J].電力系統(tǒng)自動化,1999,23(4):51-56.

[10]KAZARLIS S A, BAKIRTZIS A G, PETRIDIS V.A genetic algorithm solution to the unit commitment problem[J]. IEEE Trans on Power Systems,1996,11(1):83-92.