魏家柱

(中國(guó)大唐集團(tuán)科學(xué)技術(shù)研究院有限公司華東電力試驗(yàn)研究院,安徽 合肥 230088)

0 引言

根據(jù)國(guó)務(wù)院2016年發(fā)布的《“十三五”節(jié)能減排綜合工作方案》要求,2020年國(guó)內(nèi)萬(wàn)元生產(chǎn)總值所消耗能源和氮氧化物排放總量與2015年相比都需降低15%以上。在不同特性的機(jī)組之間合理分配負(fù)荷,研究電力系統(tǒng)的環(huán)保經(jīng)濟(jì)發(fā)電調(diào)度(environmental-economic dispatch,EED)問(wèn)題,對(duì)促進(jìn)節(jié)能減排具有重要意義。

現(xiàn)階段在研究EED問(wèn)題時(shí),一些文獻(xiàn)采用了多目標(biāo)優(yōu)化算法[1-4]。如文獻(xiàn)[1]在經(jīng)濟(jì)成本和環(huán)境成本中考慮了實(shí)際的負(fù)荷需求和多種燃料資源,并在多目標(biāo)粒子群優(yōu)化算法中引入了多目標(biāo)解集更新和變鄰域搜索2種策略,很好地降低了煤耗和發(fā)電成本。文獻(xiàn)[2]將AC-DC混合微電網(wǎng)的操作屬性以及負(fù)載和可再生能源發(fā)電的不確定性納入最佳調(diào)度問(wèn)題,同基于容量的下垂設(shè)置相比調(diào)度過(guò)程中的運(yùn)營(yíng)成本和排放量都明顯減少。文獻(xiàn)[3]提出了一種基于安全約束的最優(yōu)潮流多目標(biāo)發(fā)電調(diào)度模型,所獲取的帕累托最優(yōu)解引起不安全風(fēng)險(xiǎn)的可能性較低,對(duì)于系統(tǒng)的不確定變化表現(xiàn)出較強(qiáng)的魯棒性。文獻(xiàn)[4]為了解決電力系統(tǒng)中凸和非凸相結(jié)合的EED問(wèn)題提出了1種改進(jìn)的基于生物地理的優(yōu)化算法,并在3種不同的測(cè)試系統(tǒng)上仿真,仿真結(jié)果表明所提算法具有較好的收斂性。

相對(duì)于多目標(biāo)優(yōu)化算法,單目標(biāo)優(yōu)化算法較為簡(jiǎn)單且在運(yùn)算時(shí)間方面具有顯著的優(yōu)勢(shì),因此在負(fù)荷多目標(biāo)優(yōu)化分配的問(wèn)題中得到了廣泛的應(yīng)用[5-8]。文獻(xiàn)[5-6]分別利用加權(quán)和法和模糊集理論將對(duì)經(jīng)濟(jì)環(huán)保2個(gè)性能進(jìn)行研究的多目標(biāo)優(yōu)化問(wèn)題轉(zhuǎn)化為單目標(biāo)優(yōu)化問(wèn)題,在固定權(quán)重組合下均提高了優(yōu)化效果。但是只在固定權(quán)重下研究?jī)?yōu)化問(wèn)題,需要多次計(jì)算不同目標(biāo)權(quán)重組合下的優(yōu)化結(jié)果,不能自動(dòng)生成最優(yōu)的權(quán)重組合。

文獻(xiàn)[7]選取煤耗、變負(fù)荷時(shí)間和污染物排放做為優(yōu)化目標(biāo),通過(guò)決策理論和層次分析法簡(jiǎn)化了研究問(wèn)題并獲得了3個(gè)目標(biāo)的權(quán)重,和單個(gè)目標(biāo)優(yōu)化的結(jié)果相比整體優(yōu)化效果得到了提升。文獻(xiàn)[8]在進(jìn)行無(wú)功功率優(yōu)化時(shí)選取網(wǎng)損和節(jié)點(diǎn)電壓偏移量作為優(yōu)化指標(biāo),并利用熵權(quán)法確定了各指標(biāo)的權(quán)重,通過(guò)和遺傳算法的對(duì)比說(shuō)明了所提算法在計(jì)算時(shí)間上的優(yōu)勢(shì)。但是,不管是文獻(xiàn)[7]的層次分析法還是文獻(xiàn)[8]的熵權(quán)法得到的權(quán)重組合都沒(méi)有進(jìn)一步和其他權(quán)重組合進(jìn)行對(duì)比,難以說(shuō)明所得權(quán)重組合的先進(jìn)性。

基于上述研究現(xiàn)狀,本文在對(duì)EED模型求解時(shí)提出1種利用發(fā)電機(jī)組數(shù)量和優(yōu)化目標(biāo)個(gè)數(shù)共同確定粒子維數(shù)的自適應(yīng)多目標(biāo)權(quán)重粒子群優(yōu)化(adaptive multi-objective weights particle swarm optimization algorithm,AMWPSO)算法,在得到各機(jī)組出力的同時(shí)給出將多目標(biāo)優(yōu)化問(wèn)題轉(zhuǎn)化為單目標(biāo)優(yōu)化問(wèn)題時(shí)的唯一權(quán)重組合方案,最后通過(guò)和2組相鄰近的權(quán)重組合下的優(yōu)化結(jié)果進(jìn)行對(duì)比說(shuō)明所提優(yōu)化算法的先進(jìn)性和有效性。

1 調(diào)節(jié)容量分配的多目標(biāo)優(yōu)化模型

火電機(jī)組運(yùn)行時(shí)產(chǎn)生的海量運(yùn)行數(shù)據(jù)由邊界數(shù)據(jù)集和運(yùn)行數(shù)據(jù)集構(gòu)成[9]。在研究EED問(wèn)題時(shí),需要收集機(jī)組的成本系數(shù)和排放系數(shù)。

1.1 經(jīng)濟(jì)性單目標(biāo)函數(shù)

綜合考慮某一區(qū)域所有機(jī)組的燃煤成本,可以得到調(diào)頻成本為

式中:P i為機(jī)組i所承擔(dān)的調(diào)節(jié)功率;N為協(xié)同參與頻率調(diào)節(jié)機(jī)組個(gè)數(shù);a i、b i、ci為機(jī)組輸出有功和發(fā)電成本之間的關(guān)系系數(shù)。

1.2 環(huán)保性單目標(biāo)函數(shù)

燃煤過(guò)程中產(chǎn)生的SO2、CO2和NOx這3種氣體中,從治污成本的角度可以?xún)H僅考慮NOx的排放[10],其排放表達(dá)式為

式中:αi、βi、γi、εi和λi為機(jī)組輸出有功和 NOx排放量之間的關(guān)系系數(shù)。

1.3 權(quán)重法建立的容量分配模型

電網(wǎng)運(yùn)行的經(jīng)濟(jì)性與環(huán)保性無(wú)法同時(shí)達(dá)到最優(yōu),體現(xiàn)為某一指標(biāo)的優(yōu)化需以惡化另一指標(biāo)為代價(jià)[11]。考慮采用權(quán)重法將EED問(wèn)題轉(zhuǎn)化為單目標(biāo)優(yōu)化問(wèn)題時(shí),不同目標(biāo)之間疊加的前提條件是各目標(biāo)的量綱相同。因此在建立模型時(shí)首先應(yīng)將模型去量綱化,每個(gè)目標(biāo)函數(shù)再乘以權(quán)重因子組成新的單目標(biāo)函數(shù)[12]。根據(jù)以上思路,可以得到整合2個(gè)目標(biāo)函數(shù)后新的負(fù)荷優(yōu)化分配數(shù)學(xué)模型:

式中:we和wg分別為發(fā)電成本、排放量的權(quán)重系數(shù);{l1,l2,…,l N}和{m1,m2,…,m N}分別為經(jīng)濟(jì)指標(biāo)最優(yōu)下和環(huán)保指標(biāo)最優(yōu)下各機(jī)組負(fù)荷分配結(jié)果。

1.4 約束條件

1)考慮網(wǎng)絡(luò)損耗時(shí),其對(duì)應(yīng)的約束條件為

式中:PD為該區(qū)域所承擔(dān)的總負(fù)荷需求;PL表示該區(qū)域產(chǎn)生的網(wǎng)絡(luò)損耗,其計(jì)算表達(dá)式為

式中B,B0,B00為網(wǎng)損系數(shù)。

2)機(jī)組出力需限制在一定的區(qū)間內(nèi),可將其以不等式的形式表達(dá)為

式中Pmaxi和Pmini分別為該區(qū)域機(jī)組i的出力上限和下限。

2 負(fù)荷分配模型求解

2.1 改進(jìn)措施

在粒子群優(yōu)化算法中,每一個(gè)粒子都代表著優(yōu)化問(wèn)題的一個(gè)可能解[13],且都對(duì)應(yīng)一個(gè)適應(yīng)度函數(shù)值,以判斷當(dāng)前粒子位置距離目標(biāo)最優(yōu)位置的遠(yuǎn)近[14]。由公式(3)可知,適應(yīng)度函數(shù)由各機(jī)組出力大小和目標(biāo)的權(quán)重大小共同決定。采用粒子群優(yōu)化算法進(jìn)行求解時(shí),可以利用發(fā)電機(jī)組的數(shù)量和目標(biāo)函數(shù)的數(shù)量共同決定單個(gè)粒子的維度,并對(duì)機(jī)組出力大小和目標(biāo)權(quán)重大小一起進(jìn)行搜索。在第t次迭代時(shí),群中單個(gè)粒子的位置X t可以表示為

式中:表示第t次迭代發(fā)電機(jī)i的出力;表示第t次迭代目標(biāo)i的權(quán)重大小;N為發(fā)電機(jī)組個(gè)數(shù);M為所優(yōu)化的目標(biāo)個(gè)數(shù)。

2.2 約束條件的處理

本文通過(guò)設(shè)置平衡機(jī)組和跟隨權(quán)重的方式使求解結(jié)果滿(mǎn)足等式約束條件。選取N個(gè)機(jī)組中的機(jī)組1作為出力平衡機(jī)組,M個(gè)優(yōu)化目標(biāo)函數(shù)中目標(biāo)1的權(quán)重作為跟隨權(quán)重。

機(jī)組1的發(fā)電出力可由公式(4)(5)得到:

2.3 求解流程

根據(jù)以上思路,所提粒子群優(yōu)化算法在求解機(jī)組調(diào)節(jié)容量分配問(wèn)題時(shí)的步驟如下文所述。

1)在機(jī)組2,…,N承擔(dān)的負(fù)荷界限內(nèi)隨機(jī)選取Q個(gè)點(diǎn)作為粒子初始位置,目標(biāo)2,…,M的權(quán)重在[0,1]之間選取Q個(gè)點(diǎn)作為粒子初始位置。為了加速尋優(yōu)過(guò)程,初始化公式[15]為

式中r為0～1的隨機(jī)數(shù)。

機(jī)組1的出力初始大小和目標(biāo)函數(shù)1權(quán)重的初始大小分別由公式(8)(9)確定,粒子初始化為Q行N+M列的矩陣。

2)將初始化后的粒子代入公式(3)中的目標(biāo)函數(shù),當(dāng)前的粒子位置設(shè)定為個(gè)體極值pbest,則pbest初始化為Q行N+M列的矩陣。目標(biāo)函數(shù)最小的粒子存儲(chǔ)在全局極值gbest中,則gbest初始化為1行N+M列的矩陣。

3)逐行更新粒子速度和位置,更新公式為

式中:c1和c2為學(xué)習(xí)因子;表示第j個(gè)粒子當(dāng)前的位置;表示第j個(gè)粒子當(dāng)前的速度,其中選取的飛行粒子的總維度為N+M-2,即

采用公式(8)(9)更新機(jī)組1的出力大小和目標(biāo)函數(shù)1權(quán)重的大小。

4)將更新后的結(jié)果代入公式(3)中的目標(biāo)函數(shù),評(píng)價(jià)粒子適應(yīng)度。

5)更新后的粒子適應(yīng)度與pbest和gbest的適應(yīng)度作比較,較好的一方成為新的pbest和gbest;

6)若迭代次數(shù)不滿(mǎn)足則回到步驟3);迭代次數(shù)滿(mǎn)足后輸出gbest,得到1行N+M列的矩陣,即為優(yōu)化后各機(jī)組承擔(dān)的功率和各目標(biāo)權(quán)重的大小。

上述過(guò)程的實(shí)施流程如圖1所示,整個(gè)過(guò)程的實(shí)現(xiàn)均依靠機(jī)組的本地控制器。

3 實(shí)例分析

為了驗(yàn)證本文所提粒子群優(yōu)化算法的正確性和有效性,設(shè)定研究對(duì)象為標(biāo)準(zhǔn)的包括6臺(tái)機(jī)組的IEEE 30總線(xiàn)系統(tǒng)[16],發(fā)電機(jī)組特性參數(shù)詳見(jiàn)表1,總負(fù)荷PD=2.834 pu。

對(duì)應(yīng)的網(wǎng)損系數(shù)如下:

圖1 粒子群優(yōu)化算法的求解流程Fig.1 Solution process of particle swarm optimization algorithm

為了給決策人員提供更好的決策依據(jù),將權(quán)重和機(jī)組出力同時(shí)放入粒子中優(yōu)化,即此時(shí)的粒子維度為6+2=8,設(shè)置Q=1 000,Tmax=100。本文參考文獻(xiàn)[15]給出了計(jì)及和忽略網(wǎng)損2種情況下的計(jì)算結(jié)果。為進(jìn)一步論證本文所提優(yōu)化算法的先進(jìn)性,參考文獻(xiàn)[5]的權(quán)重設(shè)置方法同時(shí)計(jì)算了2組相鄰的權(quán)重組合下的優(yōu)化結(jié)果。

其中,忽略和計(jì)及網(wǎng)損時(shí)的單一目標(biāo)下的優(yōu)化結(jié)果{l1,l2,…,l N}和{m1,m2,…,m N}見(jiàn)表2。忽略網(wǎng)損時(shí)其適應(yīng)度函數(shù)迭代曲線(xiàn)和對(duì)比結(jié)果分別如圖2和表3所示;計(jì)及網(wǎng)損時(shí)其適應(yīng)度函數(shù)迭代曲線(xiàn)和對(duì)比結(jié)果分別如圖3和表4所示。

圖2 忽略網(wǎng)損時(shí)的適應(yīng)度函數(shù)曲線(xiàn)Fig.2 Fitness function curve when network loss is ignored

由表2中的結(jié)果能夠得出,選取不同優(yōu)化目標(biāo)時(shí)對(duì)優(yōu)化結(jié)果的影響較大。以經(jīng)濟(jì)性指標(biāo)為例,忽略網(wǎng)損時(shí)選取不同優(yōu)化目標(biāo)發(fā)電成本的差值為34.130 9$/h,計(jì)及網(wǎng)損時(shí)選取不同優(yōu)化目標(biāo)發(fā)電成本的差值為39.607 4$/h。尋找2個(gè)目標(biāo)的最優(yōu)權(quán)重對(duì)獲取經(jīng)濟(jì)和環(huán)保2個(gè)指標(biāo)的最優(yōu)綜合表現(xiàn)具有決定性的作用,這也進(jìn)一步論證了本文將各目標(biāo)權(quán)重作為粒子的一部分進(jìn)行尋優(yōu)的必要性和有效性。

表1 發(fā)電機(jī)組特性參數(shù)Table 1 Generator characteristic parameters

表2 單一目標(biāo)下的優(yōu)化結(jié)果Table 2 Optimization result when only a single goal is considered

表3 忽略網(wǎng)損時(shí)的優(yōu)化結(jié)果Table 3 Optimization result when network loss is ignored

圖3 計(jì)及網(wǎng)損時(shí)的適應(yīng)度函數(shù)曲線(xiàn)Fig.3 Fitness function curve when network loss is considered

忽略網(wǎng)損時(shí),分析圖2和表3中的結(jié)果可以發(fā)現(xiàn),80次迭代后適應(yīng)度函數(shù)的大小達(dá)到穩(wěn)定,得到的權(quán)重組合為[0.4527 0.5473]。參照文獻(xiàn)[5]中的權(quán)重設(shè)置方法和設(shè)置精度,對(duì)比權(quán)重組合[0.4000 0.6000]而言,發(fā)電成本下降2.395 4$/h,下降幅度可達(dá)0.387 3%;污染物排放升高0.000 4 t/h,升高幅度僅僅為0.198 1%。對(duì)比權(quán)重組合[0.500 0 0.500 0]而言,發(fā)電成本上升0.467 8$/h,上升幅度僅僅為0.076 0%;污染物排放降低0.001 7 t/h,降低幅度可達(dá)0.833 3%。因此,忽略網(wǎng)損時(shí)本文所提方法和文獻(xiàn)[5]中的手動(dòng)設(shè)置權(quán)重組合的方法相比綜合表現(xiàn)更佳。

計(jì)及網(wǎng)損時(shí),分析圖3和表4中的結(jié)果可以發(fā)現(xiàn),28次迭代后適應(yīng)度函數(shù)的大小達(dá)到穩(wěn)定,得到的權(quán)重組合為[0.4761 0.5239]。參照文獻(xiàn)[5]中的權(quán)重設(shè)置方法和設(shè)置精度,對(duì)比權(quán)重組合[0.4000 0.6000]而言,發(fā)電成本下降10.177 7$/h,下降幅度可達(dá)1.067 6%;污染物排放升高0.0003 t/h,升高幅度僅僅為0.152 1%。對(duì)比權(quán)重組合[0.500 0 0.500 0]而言,發(fā)電成本上升2.604 4$/h,上升幅度僅僅為0.419 9%;污染物排放降低0.004 9 t/h,降低幅度可達(dá)2.419 8%。因此,計(jì)及網(wǎng)損時(shí)利用本文所提方法求解模型仍然能夠得到更好的綜合表現(xiàn),且優(yōu)勢(shì)相對(duì)于文獻(xiàn)[5]中的手動(dòng)設(shè)置方法更加明顯。

表4 計(jì)及網(wǎng)損時(shí)的優(yōu)化結(jié)果Table 4 Optimization result when network loss is considered

綜合上述,不管是在忽略還是在計(jì)及網(wǎng)損的條件下,本文所提方法都能較好地適用于負(fù)荷多目標(biāo)優(yōu)化分配問(wèn)題的求解中。

4 結(jié)論

在將EED多目標(biāo)優(yōu)化問(wèn)題轉(zhuǎn)化為單目標(biāo)優(yōu)化問(wèn)題時(shí),利用本文所提的依據(jù)機(jī)組數(shù)量和優(yōu)化目標(biāo)數(shù)量共同決定粒子維度的AMWPSO算法,有效避免了繁瑣的權(quán)重篩選過(guò)程。且相對(duì)于相鄰的權(quán)重組合,其經(jīng)濟(jì)性和環(huán)保性的綜合性能表現(xiàn)較好。此外,本文所提方法也可以應(yīng)用于3個(gè)及以上目標(biāo)的負(fù)荷優(yōu)化分配問(wèn)題中。