劉仁孝 ,張曉華 ,胡 騰,孫 悅
(1.東北電力大學(xué),吉林 吉林 132012; 2.國網(wǎng)吉林省電力有限公司 ,長春 130028;3.中國電力工程顧問集團(tuán)東北電力設(shè)計院有限公司,長春 130021)
隨著化石資源短缺和環(huán)境污染等問題日益嚴(yán)峻,以風(fēng)力發(fā)電為代表的可再生能源在我國能源占比日益增加,截至2019年底,我國風(fēng)電累計裝機(jī)2.1×108kW,風(fēng)電裝機(jī)占全部發(fā)電裝機(jī)的10.4%,位居全球第一[1-2]。然而,由于風(fēng)電的間歇性和隨機(jī)性,其大規(guī)模并網(wǎng)給電網(wǎng)的安全穩(wěn)定運(yùn)行帶來了極大的不確定性。傳統(tǒng)的確定性調(diào)度理論難以滿足日益增加的風(fēng)電并網(wǎng)規(guī)模的要求。
風(fēng)電的不確定性導(dǎo)致電網(wǎng)難以準(zhǔn)確進(jìn)行日前調(diào)度,預(yù)留足夠的旋轉(zhuǎn)備用容量有利于應(yīng)對風(fēng)電的波動性,保證電網(wǎng)的安全穩(wěn)定運(yùn)行,文獻(xiàn)[3]引入抽水蓄能電站作為系統(tǒng)備用降低風(fēng)電出力波動對系統(tǒng)的影響,文獻(xiàn)[4]利用直接搜索法求解了風(fēng)電和熱電相協(xié)調(diào)的最優(yōu)組合備用模型。文獻(xiàn)[3-4]均沒有對系統(tǒng)運(yùn)行風(fēng)險的評估,只是以一定的比例作為系統(tǒng)的備用容量,容易導(dǎo)致由于風(fēng)險評估不足而出現(xiàn)在備用不足或預(yù)留過多。文獻(xiàn)[5]以機(jī)會約束規(guī)劃代替?zhèn)鹘y(tǒng)確定性約束條件,通過置信度對旋轉(zhuǎn)備用容量進(jìn)行了優(yōu)化求解。文獻(xiàn)[6]提出了基于區(qū)間優(yōu)化的考慮風(fēng)電不確定性的調(diào)度模型,并將模型轉(zhuǎn)化為線性規(guī)劃進(jìn)行求解。文獻(xiàn)[7]采用魯棒優(yōu)化方法,提出了衡量潛在風(fēng)險成本的方法并對多類型能源進(jìn)行了優(yōu)化。文獻(xiàn)[8]基于條件風(fēng)險價值方法構(gòu)建了電力系統(tǒng)條件風(fēng)險調(diào)度模型。上述文獻(xiàn)模型的可靠性和置信水平都是由調(diào)度員設(shè)定,置信水平的高低反應(yīng)了調(diào)度員對電力系統(tǒng)運(yùn)行可靠性與經(jīng)濟(jì)性的權(quán)衡,但對于電力系統(tǒng)的可靠性與經(jīng)濟(jì)性的平衡考慮的不夠全面,為有效處理這一問題,本文綜合考慮運(yùn)行風(fēng)險與消納成本兩方面因素,對風(fēng)-火-蓄系統(tǒng)進(jìn)行聯(lián)合優(yōu)化,在機(jī)會約束的基礎(chǔ)上引入風(fēng)險理論建立相應(yīng)的風(fēng)險量化模型進(jìn)行研究。最后利用算例驗(yàn)證了模型的有效性。
(1)
式中Wt為風(fēng)電總裝機(jī)容量。
圖1 風(fēng)功率分布概率曲線
為應(yīng)對消納區(qū)間可能出現(xiàn)的風(fēng)電預(yù)測誤差,系統(tǒng)需預(yù)留相應(yīng)的備用容量,保證其安全可靠運(yùn)行。根據(jù)風(fēng)電消納區(qū)間對系統(tǒng)備用容量的要求,針對系統(tǒng)備用,基于機(jī)會約束規(guī)劃(chance constrained programming,CCP)方法可建立相應(yīng)約束如下:
Pr{P+Ru+Pw+ew≥Pd}≥β
(2)
Pr{P-Rd+Pw+ew≤Pd}≥β
(3)
(4)
綜合考慮系統(tǒng)運(yùn)行風(fēng)險與消納成本,優(yōu)化風(fēng)電消納區(qū)間,合理配置系統(tǒng)備用容量。可建立風(fēng)險模型如下:
(5)
本文綜合考慮風(fēng)-火-蓄系統(tǒng)的發(fā)電特性與運(yùn)行約束,合理地考慮運(yùn)行風(fēng)險和消納成本使整個聯(lián)合系統(tǒng)的運(yùn)行成本F達(dá)到最優(yōu)。
minF=Fr+FG+FE
(6)
式中:FG為常規(guī)機(jī)組運(yùn)行成本,包括機(jī)組啟停成本、煤耗成本、備用成本;FE為抽水儲能成本,包括儲能運(yùn)行成本、備用成本。本文暫不考慮抽水蓄能機(jī)組啟動時的成本,只考慮抽水蓄能機(jī)組管理運(yùn)行維護(hù)費(fèi)用和基建成本平攤到各協(xié)調(diào)周期。
2.2.1 火電機(jī)組約束
火電機(jī)組的特性約束包括功率約束、備用約束、爬坡約束、邏輯約束和運(yùn)行時間約束:
(7)
(8)
(9)
Vi,t-Vi,t-1=yi,t-zi,t
(10)
2.2.2 抽水蓄能電站約束
抽水蓄能電站系統(tǒng)約束包括抽水蓄能電站功率約束、庫容運(yùn)行約束與備用約束:
(11)
(12)
(13)
(14)
(15)
(16)
2.2.3 系統(tǒng)運(yùn)行約束
(17)
由于概率運(yùn)算與積分運(yùn)算的存在,因此模型為復(fù)雜非線性模型,求解難度較大,需要利用計算性能較好的智能算法進(jìn)行求解。本文將通過數(shù)學(xué)轉(zhuǎn)化將模型進(jìn)行化簡,降低求解難度。
a.機(jī)會約束轉(zhuǎn)化。根據(jù)文獻(xiàn)[10] 中定理指出,機(jī)會約束內(nèi)不等式中的隨機(jī)變量可以從多項(xiàng)式中分離時,可實(shí)現(xiàn)將機(jī)會約束由不確定形式向確定形式的轉(zhuǎn)化,假設(shè)隨機(jī)變量的累計分布函數(shù)為φ,則可以轉(zhuǎn)化為如下確定形式:
Pd-P+Ru+Pw≤φ-1(1-β)
(18)
P-Rd+Pw-Pd≤-φ-1(β)
(19)
(20)
圖2 分段線性化示意圖
經(jīng)過以上處理,模型最終轉(zhuǎn)化為混合整數(shù)二次規(guī)劃模型問題。這樣就可以借助商用求解器CPLEX高效求解。具體步驟為:讀入風(fēng)電、負(fù)荷和儲能等數(shù)據(jù),讀入系統(tǒng)模型基本參數(shù);對輸入的相關(guān)數(shù)據(jù)進(jìn)行計算處理;構(gòu)建待求解問題的數(shù)學(xué)模型,包括目標(biāo)函數(shù)和各類等式約束;調(diào)用CPLEX求解器對模型進(jìn)行求解;輸出最優(yōu)解作為最優(yōu)運(yùn)行方案。
為驗(yàn)證本文所提方法的有效性和適用性,以簡單四元單母線為例進(jìn)行測試(見圖3),包含2個常規(guī)機(jī)組、1座抽水蓄能電站、1個風(fēng)電場、1個負(fù)荷。
圖3 四單元母線源系統(tǒng)
考慮到切負(fù)荷風(fēng)險對系統(tǒng)安全問題的影響遠(yuǎn)大于棄風(fēng)風(fēng)險,因此,系統(tǒng)的切負(fù)荷風(fēng)險成本系數(shù)取300元/MW,棄風(fēng)風(fēng)險成本因子取4.2元/MW。在利用分段線性化處理積分函數(shù)時,考慮到模型精度,分段數(shù)k取50?;痣姍C(jī)組參數(shù)見表1,可調(diào)備用成本為其煤耗成本的15%,抽水蓄能機(jī)組參數(shù)為最大放電/抽水功率25/25 MW,最大/最小庫容100/10(MW·h),放電/抽水效率(85/95)%,運(yùn)行成本825元/(MW·h)。
利用本文所提方法,對系統(tǒng)24 h日前調(diào)度進(jìn)行建模求解,機(jī)會約束條件置信水平 ,典型日風(fēng)電預(yù)測出力曲線、負(fù)荷預(yù)測曲線和凈負(fù)荷曲線見圖 4。
圖4 風(fēng)電預(yù)測出力曲線、負(fù)荷預(yù)測曲線和凈負(fù)荷曲線
根據(jù)本文所提日前調(diào)度模型,求解所得調(diào)度計劃及風(fēng)電每小時的波動預(yù)測概率見圖5。系統(tǒng)向下的消納區(qū)間會比設(shè)定區(qū)間大一些,這是由于本文所提模型會在機(jī)會約束設(shè)定置信水平的基礎(chǔ)上判斷系統(tǒng)所承擔(dān)的風(fēng)險來確定上、下可調(diào)備用容量??梢源龠M(jìn)系統(tǒng)備用合理配置,更好地調(diào)整風(fēng)電的消納區(qū)間。
圖6為系統(tǒng)正備用與負(fù)備用的風(fēng)險概率和風(fēng)險期望。風(fēng)電消納區(qū)間的下邊界與系統(tǒng)正備用相互對應(yīng),系統(tǒng)正備用在整個研究周期都出現(xiàn)了調(diào)整,只有13 h沒有調(diào)整,此時系統(tǒng)正備用不足的風(fēng)險概率少量增加,此時正備用風(fēng)險增至 6.13%,與相對應(yīng)的切負(fù)荷風(fēng)險比原來增加了1.01 MW。未出現(xiàn)調(diào)整的原因是由于凈負(fù)荷曲線在13 h出現(xiàn)了短時波峰,若此時使風(fēng)險降對最低則需啟動額外的機(jī)組且短時間運(yùn)行后又需停機(jī),這對整個系統(tǒng)來說是不經(jīng)濟(jì)的,對火電機(jī)組來說也是不利的。若參照優(yōu)化結(jié)果,則可以避免額外的啟停計劃,減少整個系統(tǒng)的成本費(fèi)用,改善系統(tǒng)運(yùn)行風(fēng)電消納區(qū)間的上邊界與系統(tǒng)負(fù)備用相互對應(yīng),由圖6可知,風(fēng)電消納區(qū)間在上邊界幾乎與設(shè)定值完全重合,這是由于棄風(fēng)風(fēng)險因子設(shè)定較低,模型將棄風(fēng)風(fēng)險在初始設(shè)定的基礎(chǔ)上最大化棄風(fēng)風(fēng)險,此時可以適當(dāng)修改風(fēng)險系數(shù)來調(diào)整棄風(fēng)風(fēng)險,有助于優(yōu)化系統(tǒng)調(diào)度計劃,提高系統(tǒng)調(diào)度經(jīng)濟(jì)性。此外,通過對比優(yōu)化后的系統(tǒng)正負(fù)備用和風(fēng)險概率可知,模型可以根據(jù)風(fēng)險情況,優(yōu)先避免切負(fù)荷風(fēng)險的增加,降低電力系統(tǒng)的運(yùn)行風(fēng)險。
圖5 預(yù)測概率
分析同一置信水平下,不同風(fēng)險成本系數(shù)對模型優(yōu)化結(jié)果的影響,以風(fēng)險成本系數(shù)為例(見圖7)。
向下可調(diào)備用容量隨著風(fēng)險成本系數(shù)的增加而增加,風(fēng)險成本系統(tǒng)增大表示系統(tǒng)的可靠性增加,風(fēng)險成本在整個目標(biāo)函數(shù)的比重增加。此時,需增大系統(tǒng)的備用容量使系統(tǒng)風(fēng)險降低,減少風(fēng)險成本,即增加系統(tǒng)的可調(diào)備用容量,當(dāng)風(fēng)險成本系數(shù)增加到某一數(shù)值時,此時系統(tǒng)備用容量將達(dá)到最大值。
表1 火電機(jī)組參數(shù)
圖6 系統(tǒng)正備用與負(fù)備用
圖7 優(yōu)化后的風(fēng)電消納區(qū)間
不同方法的成本結(jié)果對比,將本文所提方法與固定比例法(一般取負(fù)荷的10%)、傳統(tǒng)CCP備用方法進(jìn)行對比結(jié)果見表2。由表2可知,傳統(tǒng)固定比例法配置的旋轉(zhuǎn)備用容量沒有考慮風(fēng)電波動的差異性,在某些時段配置的可調(diào)備用容量不足,導(dǎo)致切負(fù)荷量及棄風(fēng)量增多,因此在該備用策略下系統(tǒng)的風(fēng)險成本和綜合成本較高。而傳統(tǒng)CCP雖然在備用成本上有所增加,但風(fēng)險成本降低了很多,總成本上下降了27 467.2元;而本文所提方法在傳統(tǒng)CCP的基礎(chǔ)上雖然增加了少許備用成本,但系統(tǒng)總成本仍然降低了6 853.79元,本文所建的備用模型在風(fēng)電出力較大的時段合理地增加向上、向下可調(diào)備用容量,降低了系統(tǒng)應(yīng)對風(fēng)電出力不確性的風(fēng)險成本,而在風(fēng)電出力較小時段合理地減少向上、向下可調(diào)備用容量,使發(fā)電機(jī)組可以運(yùn)行在效率更高的發(fā)電狀態(tài),降低了整個系統(tǒng)的運(yùn)行成本。
針對風(fēng)電出力隨機(jī)性和波動性給電力系統(tǒng)優(yōu)化運(yùn)行帶來的新難題,本文在機(jī)會約束的基礎(chǔ)上引入風(fēng)險理論,構(gòu)建了考慮風(fēng)電出力誤差概率的棄風(fēng)風(fēng)險、切負(fù)荷風(fēng)險模型,并以系統(tǒng)運(yùn)行成本最小化作為目標(biāo),建立了風(fēng)-火-蓄聯(lián)合運(yùn)行優(yōu)化模型,合理地權(quán)衡了系統(tǒng)風(fēng)電消納區(qū)間調(diào)整的成本和風(fēng)險。算例結(jié)果表明:
a.本文所提模型可以在人為規(guī)定置信度的基礎(chǔ)上識別由風(fēng)力發(fā)電引起的調(diào)度不確定性合理權(quán)衡風(fēng)險成本與發(fā)電運(yùn)行成本的關(guān)系合理地配置電力系統(tǒng)的備用容量;
b.與傳統(tǒng)CCP方法和固定比例方法相比,本文所提方法將風(fēng)險表示為調(diào)度決策的函數(shù)。有助于制定調(diào)度計劃,更直觀地了解發(fā)電成本和風(fēng)險成本的限制;
c.將本文的模型中機(jī)會約束的不確定部分轉(zhuǎn)化為確定性問題、積分問題進(jìn)行分段線性化,將原模型重新表述為混合整數(shù)型規(guī)劃問題,不增加額外的變量。這種緊湊的建模有助于提高本文模式的計算效率。
表2 不同方案對比 元