謝國輝,張粒子,舒 雋,王家興

(1.華北電力大學(xué)電氣與電子工程學(xué)院,北京 102206;2.福建省福州超高壓輸變電局,福州 350013)

節(jié)能發(fā)電調(diào)度以節(jié)能、環(huán)保為目標(biāo),以全電力系統(tǒng)內(nèi)發(fā)、輸、供電設(shè)備為調(diào)度對象,優(yōu)先調(diào)度可再生和清潔發(fā)電資源,按能耗和污染物排放水平,由低到高依次調(diào)用化石類發(fā)電資源,最大限度地減少能源、資源消耗和污染物排放[1]。研究節(jié)能發(fā)電調(diào)度日前機(jī)組組合優(yōu)化方法,是制定日前節(jié)能發(fā)電調(diào)度計劃的關(guān)鍵問題。

目前,關(guān)于制定機(jī)組組合方法已有研究成果。針對傳統(tǒng)機(jī)組組合模型,國內(nèi)外眾多專家學(xué)者已試驗了大量的數(shù)學(xué)規(guī)劃方法,如:優(yōu)先次序法、動態(tài)規(guī)劃法[2]、拉格朗日松弛法[3]、智能方法[4]和混合整數(shù)規(guī)劃法等。針對電力市場環(huán)境下機(jī)組組合方法研究,文獻(xiàn)[5]提出了在電力市場機(jī)制下機(jī)組組合問題的數(shù)學(xué)模型,并運用遺傳算法求解。文獻(xiàn)[6]提出一種組合優(yōu)化方法解決電力市場環(huán)境下機(jī)組組合問題;文獻(xiàn)[7]提出一種用于求解中期電力市場機(jī)組組合問題啟發(fā)式方法。針對節(jié)能發(fā)電調(diào)度環(huán)境,文獻(xiàn)[8]提出節(jié)能發(fā)電調(diào)度發(fā)電計劃的編制算法,首先根據(jù)火電機(jī)組發(fā)電序位表確定機(jī)組開機(jī)方法,然后采用等微增率方法進(jìn)行功率分配。

上述研究是針對傳統(tǒng)經(jīng)濟(jì)調(diào)度和電力市場,目前針對節(jié)能發(fā)電調(diào)度環(huán)境制定火電機(jī)組組合方法的研究文獻(xiàn)較少。本文針對節(jié)能發(fā)電調(diào)度模式,提出分步優(yōu)化和聯(lián)合優(yōu)化的兩種制定火電機(jī)組日前節(jié)能發(fā)電調(diào)度機(jī)組組合方法。分步優(yōu)化方法基于機(jī)組排序表確定機(jī)組的啟停狀態(tài)再進(jìn)行經(jīng)濟(jì)負(fù)荷分配;聯(lián)合優(yōu)化方法則直接建立以能耗最小為目標(biāo)、考慮多種約束的日前節(jié)能發(fā)電調(diào)度機(jī)組組合模型,再基于啟發(fā)式動態(tài)規(guī)劃結(jié)合等微增率的混合算法求解。最后通過實例對兩類方法進(jìn)行比對分析。

1 機(jī)組組合的分步優(yōu)化方法

1.1 機(jī)組啟停狀態(tài)確定

根據(jù)機(jī)組的排序表確定每臺機(jī)組發(fā)電的優(yōu)先順序。機(jī)組排序表根據(jù)火電機(jī)組的設(shè)計標(biāo)準(zhǔn)煤耗水平由低到高排序形成,設(shè)計煤耗相同時,按機(jī)組污染物排放水平確定優(yōu)先順序。在機(jī)組發(fā)電優(yōu)先順序基礎(chǔ)上確定機(jī)組啟停狀態(tài),具體步驟如下:

1)根據(jù)機(jī)組優(yōu)先順序,滿足t時段負(fù)荷和備用需求,確定機(jī)組t時段啟停狀態(tài);

2)進(jìn)行機(jī)組t-1時段啟停狀態(tài)向t時段狀態(tài)轉(zhuǎn)移判斷;

3)若滿足機(jī)組啟停約束,轉(zhuǎn)入步驟4);若不滿足,進(jìn)行t時段啟停狀態(tài)調(diào)整,然后轉(zhuǎn)入步驟4);

4)t=t+1,如果t不大于總時段數(shù)T,轉(zhuǎn)入步驟1),確定機(jī)組t時段啟停狀態(tài);反之停止,得到機(jī)組全時段的啟停狀態(tài)。

1.2 經(jīng)濟(jì)負(fù)荷分配

在上述確定的機(jī)組啟停狀態(tài)基礎(chǔ)上,進(jìn)行每臺機(jī)組全時段經(jīng)濟(jì)負(fù)荷分配,優(yōu)化數(shù)學(xué)模型為

約束條件:

功率平衡約束

機(jī)組功率約束

機(jī)組爬坡速度約束

直流潮流約束

式(1)～(6)中:i為機(jī)組號(i=1,2,…,It),It為每個時段入圍機(jī)組數(shù);t為時段號(t=1,2,…,T),T為時段數(shù)為機(jī)組i在t時段的有功出力,單位MW為機(jī)組i的煤耗特性曲線,可以通過參數(shù)實測或配套的煤耗在線監(jiān)測設(shè)備實時獲得。PD,t為第t時段的負(fù)荷分別為機(jī)組i的最小和最大出力;Pupi、Pdowni分別為機(jī)組i的升降出力速率;Sj,k為第j節(jié)點注入功率對第k條線路潮流的靈敏度為第j節(jié)點的凈注入功率;為第k條線路的穩(wěn)定極限。

1.3 求解算法

直流潮流約束下,模型(1)中所有約束條件都為線性約束,目標(biāo)函數(shù)中的機(jī)組耗量特性曲線函數(shù)可以近似用二次函數(shù)進(jìn)行擬合,故本文采用二次規(guī)劃算法求解模型(1)。

2 機(jī)組組合的聯(lián)合優(yōu)化方法

分步優(yōu)化方法簡易快速,不過,由于機(jī)組啟停狀態(tài)和經(jīng)濟(jì)負(fù)荷分配二者相互依存,分步優(yōu)化方法對二者進(jìn)行了簡單解耦優(yōu)化,可能不能獲得最佳的優(yōu)化效果。為了精細(xì)化制定日前發(fā)電調(diào)度機(jī)組組合,本文提出聯(lián)合優(yōu)化方法。

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

以降低系統(tǒng)煤耗為目標(biāo),并考慮啟動煤耗,建立日前節(jié)能發(fā)電調(diào)度的機(jī)組組合模型,通過求解較完備的機(jī)組組合模型實現(xiàn)確定機(jī)組啟停狀態(tài)和經(jīng)濟(jì)負(fù)荷分配的聯(lián)合優(yōu)化。

2.2 約束條件

啟停變量約束

功率平衡約束

機(jī)組功率約束

最小啟停時間約束

負(fù)荷備用約束

直流潮流約束

式(7)～(16)中:i為機(jī)組號(i=1,2,…,I),I為所有機(jī)組總數(shù);U i,t為機(jī)組i在t時段運行狀態(tài),Ui,t=1運行,U i,t=0停機(jī);SCi為機(jī)組i的啟動煤耗,為停機(jī)時間的函數(shù)分別為機(jī)組i允許最小連續(xù)停機(jī)時間和允許最小連續(xù)開機(jī)時間;Rt為t時段系統(tǒng)備用需求。

2.3 求解算法

上述建立的聯(lián)合優(yōu)化模型是復(fù)雜多約束的機(jī)組組合模型。本文采用啟發(fā)式動態(tài)規(guī)劃結(jié)合等微增率法的混合優(yōu)化算法求解模型(7),整個算法流程如圖1所示。

圖1 聯(lián)合優(yōu)化算法流程Fig.1 Algorithm o f joint optim ization

1)確定各時段可能的啟停狀態(tài)組合

首先,根據(jù)每個時段的負(fù)荷需求,按上述確定的機(jī)組發(fā)電優(yōu)先順序排序確定每個時段的邊際機(jī)組。然后,在邊際機(jī)組上下范圍內(nèi)適當(dāng)增加或減少機(jī)組數(shù)目,同時滿足系統(tǒng)的備用需求,形成每個時段可能的機(jī)組組合狀態(tài)空間集。

2)最優(yōu)路徑轉(zhuǎn)移策略

最優(yōu)路徑轉(zhuǎn)移策略流程如圖2。

綜上所述：腦卒中篩查中采納頸動脈超聲，可有效診斷出頸動脈狹窄程度及其血流參數(shù)，并分析其超聲特征，值得臨床信賴并進(jìn)一步推廣。

圖2 最優(yōu)路徑轉(zhuǎn)移策略Fig.2 Strategy of optimal shift path

3 算例分析

采用某電網(wǎng)中10臺火電機(jī)組數(shù)據(jù),分別采用分步優(yōu)化和聯(lián)合優(yōu)化兩種方法進(jìn)行日前24時段節(jié)能發(fā)電調(diào)度機(jī)組組合模擬分析。為簡化分析,暫不考慮電網(wǎng)安全約束。日前24時段負(fù)荷預(yù)測見圖3。

圖3 日前24時段負(fù)荷預(yù)測Fig.3 Load forecast of day-ahead 24 periods

每個時段系統(tǒng)的備用需求取為該時段負(fù)荷的10%。本文算例根據(jù)機(jī)組出力點對應(yīng)的煤耗值(已折算為標(biāo)準(zhǔn)煤),采用二次曲線擬合機(jī)組煤耗特性曲線,具體參數(shù)見表1,表中a、b、c分別為擬合二次曲線的二次項、一次項和常數(shù)項系數(shù)。

表1 二次曲線擬合的煤耗參數(shù)Tab.1 Coal consump tion parameters using conicfitting

3.1 分步優(yōu)化模擬分析

(1)機(jī)組優(yōu)先發(fā)電順序

根據(jù)機(jī)組的設(shè)計煤耗水平進(jìn)行排序,4號機(jī)組和6號機(jī)組的設(shè)計煤耗水平相同,由于4號機(jī)組尚沒有裝脫硫設(shè)備,故考慮環(huán)保因素6號機(jī)組優(yōu)先順序高于4號機(jī)組,那么最終確定的機(jī)組優(yōu)先發(fā)電順序為:1、2、3、6、4、5、7、8、9和10。

(2)優(yōu)化結(jié)果分析

分步優(yōu)化的機(jī)組組合情況見表2。

表2 分步優(yōu)化機(jī)組組合Tab.2 Unit comm itmentby step optim ization

3.2 聯(lián)合優(yōu)化模擬分析

(1)聯(lián)合優(yōu)化結(jié)果見表3。

表3 聯(lián)合優(yōu)化機(jī)組組合Tab.3 Unit comm itment by joint op tim ization

(2)與分布優(yōu)化對比分析

分步優(yōu)化與聯(lián)合優(yōu)化二者最明顯差異是3、6號和5號機(jī)組,9號和10號機(jī)組的啟停狀態(tài)上。3、6號發(fā)電優(yōu)先順序都高于5號機(jī)組,故分步優(yōu)化方法下3、6號機(jī)組全時段基本不停機(jī);而聯(lián)合優(yōu)化方式下5號機(jī)組全時段啟機(jī)運行,相反3、6號機(jī)組停機(jī)的時段則較多。雖然5號機(jī)組發(fā)電優(yōu)先順序較低,但5號機(jī)組煤耗曲線的c值相比3、6號機(jī)組要小得多,由于聯(lián)合優(yōu)化方法擴(kuò)展了機(jī)組啟停狀態(tài)的尋優(yōu)空間,故權(quán)衡多方面因素后選擇5號機(jī)組是更優(yōu)結(jié)果。9號和10號機(jī)組情況與此類似?？梢?分步優(yōu)化方法能夠嚴(yán)格體現(xiàn)機(jī)組的發(fā)電優(yōu)先順序,而聯(lián)合優(yōu)化方法則更體現(xiàn)了算法的尋優(yōu)能力。

進(jìn)一步對比兩種方法下每個時段的運行煤耗。根據(jù)經(jīng)濟(jì)負(fù)荷分配優(yōu)化得到的每臺機(jī)組在各個時段的出力,計算各時段兩種方法運行煤耗的差值,正值表明分步優(yōu)化方法運行煤耗高于聯(lián)合優(yōu)化方法,如圖4所示。

圖4 兩種優(yōu)化方法各時段運行煤耗差值Fig.4 Period coal-consump tion dif ferences by two op timalmethods

由于每個時段機(jī)組的啟停狀態(tài)不同,故兩種優(yōu)化方法下每個時段的運行煤耗存在明顯差異。除7、13和24三個時段外,聯(lián)合優(yōu)化方法每個時段的運行煤耗都低于分步優(yōu)化方法,從全時段運行煤耗來看,聯(lián)合優(yōu)化方法相比分步優(yōu)化方法減少標(biāo)準(zhǔn)煤108540 kg。計及啟動煤耗后,聯(lián)合優(yōu)化總標(biāo)準(zhǔn)煤耗量為13428750 kg,而分步優(yōu)化為13530450 kg,減小了101700 kg。可見,聯(lián)合優(yōu)化方法要優(yōu)于分步優(yōu)化方法。

4 結(jié)語

本文提出了制定火電機(jī)組日前節(jié)能發(fā)電調(diào)度機(jī)組組合的兩種優(yōu)化方法。同一算例對比分析表明,分步優(yōu)化方法嚴(yán)格遵循機(jī)組排序表,而聯(lián)合優(yōu)化方法則更體現(xiàn)了算法的尋優(yōu)能力,故聯(lián)合優(yōu)化方法的優(yōu)化效果更佳。其機(jī)組組合制定方案計及啟動煤耗后全天共節(jié)省標(biāo)準(zhǔn)煤101700 kg,節(jié)省近0.8個百分點。不過,聯(lián)合優(yōu)化需要獲得較精確的煤耗曲線,并且需要提高算法執(zhí)行效率從而提高計算速度,這將做進(jìn)一步深入研究。

[1] 國務(wù)院辦公廳.國務(wù)院辦公廳關(guān)于轉(zhuǎn)發(fā)發(fā)展改革委等部門節(jié)能發(fā)電調(diào)度辦法(試行)的通知[EB/OL].http://www.gov.cn/gongbao/content/2007/content_744115.htm,2007.

[2] Pang Chok Khon,Sheble Gerald B,A lbuyyeh Farrokh.Evaluation of dynamic programm ing based method and mu ltiple area representation for thermal unit comm itment[J].IEEE Trans on Power Appara-tus and Systems,1981,100(3):1212-1218.

[3] Ongsakul Weerakorn,Petcharaks Nit.Unit comm itmentby enhanced adap tive Lagrangian relaxation[J].IEEE Trans on Pow er System s,2004,19(1):620-628.

[4] Juste K A,Kita H,Tanaka E,et a l.Evolutionary programm ing solution to the unit comm itment problem[J].IEEE T rans on Pow er System s,1999,14(4):1452-1459.

[5] 劉繼春,李興源(Liu Jichun,LiXingyuan).電力市場中機(jī)組組合的智能優(yōu)化算法(Intelligent optimization app roach to unit comm itment in powermarket)[J].電力系統(tǒng)及其自動化學(xué)報(Proceedings of the CSUEPSA),2002,14(5):48-51.

[6] 王承民,郭志忠(W ang Chengm in,Guo Zhizhong).電力市場環(huán)境下解決機(jī)組組合問題的新方法(A new method of unit comm itment in electrical pow er market)[J].電力自動化設(shè)備(Electric Power Automation Equipment),2001,21(11):14-17.

[7] 楊秀媛,張南,王海寧,等(Yang Xiuyuan,Zhang Nan,Wang Haining,et a l).用于中期電力市場的啟發(fā)式機(jī)組組合算法(A heuristic unit comm itment algorithm formedium-term electricitymarket)[J].電網(wǎng)技術(shù)(Pow er System Techno logy),2004,28(24):21-24.

[8] 施建華,譚素梅(Shi Jianhua,Tan Sumei).節(jié)能發(fā)電調(diào)度發(fā)電計劃編制算法(A lgorithm of energy saving generation dispatch scheduling)[J].電力系統(tǒng)自動化(Automation o f Electric Pow er Systems),2008,32(24):48-51.