李 昱,程 強(qiáng),孫高龍,俞 雷,李海濤
(1.中國電建集團(tuán)貴陽勘測設(shè)計(jì)研究院有限公司,貴陽 550000;2.國網(wǎng)湖北省電力有限公司超高壓公司,湖北宜昌 443000)
在眾多可再生能源中,風(fēng)電是目前發(fā)展最快最成熟的發(fā)電技術(shù)[1]。我國幅員遼闊,海岸線狹長,海上風(fēng)電具有廣闊的開發(fā)前景[2]。隨著海上風(fēng)電規(guī)模的不斷增長,海上風(fēng)電并網(wǎng)量越來越大,風(fēng)電的隨機(jī)性和波動(dòng)性對(duì)電力系統(tǒng)的安全性和穩(wěn)定性帶來了巨大挑戰(zhàn)[3]。因此,開展海上風(fēng)電并網(wǎng)調(diào)度研究,建立更合理的海上風(fēng)電并網(wǎng)調(diào)度模型,對(duì)于提高風(fēng)電利用率、減小棄風(fēng)電量具有重要意義。
近年來,海上風(fēng)電并網(wǎng)調(diào)度問題成為研究的熱點(diǎn),國內(nèi)外專家學(xué)者對(duì)此進(jìn)行了大量研究。文獻(xiàn)[4]為了降低海上風(fēng)電的調(diào)度成本,通過增加約束條件的方式完善海上風(fēng)電優(yōu)化調(diào)度模型,并采用實(shí)際風(fēng)電場數(shù)據(jù)對(duì)模型的正確性進(jìn)行驗(yàn)證。文獻(xiàn)[5]為了減小海上風(fēng)電并網(wǎng)對(duì)電力系統(tǒng)的影響,提出一種海水抽水蓄能電站與海上風(fēng)電場聯(lián)合運(yùn)行的優(yōu)化調(diào)度模型,并用實(shí)際算例驗(yàn)證了模型的正確性和有效性。文獻(xiàn)[6]為了滿足海上風(fēng)電調(diào)度的經(jīng)濟(jì)性要求,分析了經(jīng)濟(jì)調(diào)度的相關(guān)特征和管理流程,建立了基于機(jī)會(huì)約束規(guī)劃理論的海上風(fēng)電經(jīng)濟(jì)調(diào)度模型。文獻(xiàn)[7]建立了海上風(fēng)電智慧安全調(diào)度平臺(tái),利用該平臺(tái)的監(jiān)測和預(yù)警功能,有效減少了海上風(fēng)電非計(jì)劃停運(yùn)時(shí)間,提高了海上風(fēng)電并網(wǎng)率。綜上所述,目前海上風(fēng)電調(diào)度模型和平臺(tái)雖然很多,但仍存在一些不足,其目標(biāo)函數(shù)和約束條件有待進(jìn)一步完善。
基于此,本研究以海上風(fēng)電運(yùn)行成本、火力發(fā)電成本、備用容量成本和棄風(fēng)懲罰成本等組成的綜合成本為目標(biāo)函數(shù),綜合考慮系統(tǒng)平衡、電網(wǎng)安全、機(jī)組運(yùn)行、旋轉(zhuǎn)備用、機(jī)組爬坡和支路潮流等各類約束條件,建立基于Petri網(wǎng)理論的海上風(fēng)電并網(wǎng)調(diào)度模型,采用實(shí)際算法驗(yàn)證模型的正確性和實(shí)用性。
Petri網(wǎng)理論是數(shù)學(xué)家Carl提出的,它是一種系統(tǒng)建模分析工具,能夠?qū)﹄x散事件系統(tǒng)中各類關(guān)系進(jìn)行描述,應(yīng)用廣泛[8]。Petri網(wǎng)的原理及分析步驟可參考文獻(xiàn)[9]。
根據(jù)Petri網(wǎng)理論,以增加約束條件的方式優(yōu)化海上風(fēng)電并網(wǎng)調(diào)度流程,減少棄風(fēng)電量,提高海上風(fēng)電利用率,盡可能降低海上風(fēng)力發(fā)電產(chǎn)生的各項(xiàng)成本。圖1給出了基于Petri網(wǎng)理論的海上風(fēng)電并網(wǎng)調(diào)度流程。
圖1 海上風(fēng)電并網(wǎng)調(diào)度流程
在圖1中,海上風(fēng)電并網(wǎng)調(diào)度模型共有四層,即任務(wù)層、任務(wù)描述層、資源層和任務(wù)結(jié)束層。各層功能如下:
(1)任務(wù)層:包括常規(guī)機(jī)組、備用機(jī)組、風(fēng)電場輸出功率和負(fù)荷需求等數(shù)據(jù),其作用是根據(jù)相關(guān)數(shù)據(jù)明確調(diào)度目標(biāo),以確定機(jī)組發(fā)電順序;
(2)任務(wù)描述層:結(jié)合調(diào)度任務(wù)制定相應(yīng)的調(diào)度計(jì)劃,當(dāng)調(diào)度計(jì)劃變遷后,執(zhí)行完本次調(diào)度計(jì)劃后轉(zhuǎn)移至下一個(gè),任務(wù)描述層中的機(jī)組運(yùn)行順序表、運(yùn)行總成本和污染物排放的懲罰成本均可觸發(fā)該條件;
(3)資源層:代表系統(tǒng)中的所有約束,比如系統(tǒng)平衡約束、機(jī)組運(yùn)行約束、電網(wǎng)安全約束、旋轉(zhuǎn)備用約束、機(jī)組爬坡約束、支路潮流約束等,建模時(shí)需要充分考慮海上風(fēng)電并網(wǎng)調(diào)度模型的各種約束;
(4)任務(wù)結(jié)束層:調(diào)度任務(wù)結(jié)束后,將該任務(wù)從資源層中移除,本次調(diào)度任務(wù)結(jié)束。
根據(jù)Petri網(wǎng)理論建立的海上風(fēng)電并網(wǎng)調(diào)度模型,構(gòu)建以綜合成本最小為調(diào)度目標(biāo)函數(shù),綜合成本的表達(dá)式如下:
式中:CZ為綜合成本;CY為海上風(fēng)電運(yùn)行成本;CF為系統(tǒng)發(fā)電成本;CB為備用容量成本;CQ為棄風(fēng)懲罰成本。
(1)海上風(fēng)電運(yùn)行成本
海上風(fēng)電運(yùn)行成本指海上風(fēng)電場在運(yùn)行過程中產(chǎn)生的成本,其表達(dá)式為:
式中:h為時(shí)段,h=1,2,…,H,H為總時(shí)段;i為機(jī)組編號(hào),i=1,2,…,N,N為機(jī)組總數(shù);Pi,h為h時(shí)段內(nèi)第i臺(tái)機(jī)組的風(fēng)電輸出功率;Ii,h為h時(shí)段內(nèi)第i臺(tái)機(jī)組的狀態(tài);Fci(Pi,h,Ii,h)為h時(shí)段內(nèi)第i臺(tái)機(jī)組的運(yùn)行成本,也稱第i臺(tái)機(jī)組的特性函數(shù)。
特性函數(shù)通常由多項(xiàng)式構(gòu)成,本研究中特性函數(shù)選取如下二次函數(shù),其表達(dá)式為:
式中:ai、bi、ci分別為發(fā)電機(jī)組的費(fèi)用系數(shù);Pi,min為第i臺(tái)機(jī)組最小輸出功率;Pi,max為第i臺(tái)機(jī)組最大輸出功率。
(2)系統(tǒng)發(fā)電成本
系統(tǒng)發(fā)電成本指系統(tǒng)內(nèi)機(jī)組發(fā)電產(chǎn)生的成本,其表達(dá)式為:
式中:kY為機(jī)組發(fā)電成本系數(shù)。
(3)備用容量成本
考慮到風(fēng)電的隨機(jī)性和不確定性,需要設(shè)置備用容量,該成本的表達(dá)式為:
式中:kB為備用容量成本系數(shù);Pact為風(fēng)電實(shí)際輸出功率。
(4)棄風(fēng)懲罰成本
風(fēng)電機(jī)組輸出功率過剩時(shí),會(huì)產(chǎn)生棄風(fēng)懲罰成本,其表達(dá)式為:
式中:kQ為棄風(fēng)懲罰成本系數(shù)。
綜上所述,海上風(fēng)電并網(wǎng)調(diào)度模型最終的目標(biāo)函數(shù)為:
海上風(fēng)電并網(wǎng)調(diào)度模型通??紤]的約束條件有系統(tǒng)平衡約束、電網(wǎng)安全約束和機(jī)組運(yùn)行約束[10]。為了使模型更接近實(shí)際情況,本研究還考慮了備用容量、機(jī)組爬坡和支路潮流等因素的影響,并分別增加相應(yīng)的約束條件,因此,考慮的約束條件共有6個(gè),具體如下。
(1)系統(tǒng)平衡約束
式中:PD,h為h時(shí)段內(nèi)系統(tǒng)負(fù)載值;m為風(fēng)電機(jī)組序號(hào);W為風(fēng)電機(jī)組總數(shù);Pf,m,h為h時(shí)段內(nèi)機(jī)組m的計(jì)劃輸出功率。
(2)電網(wǎng)安全約束
(3)機(jī)組運(yùn)行約束
式中:Pi,min、Pi,max分別為第i臺(tái)機(jī)組輸出功率的上、下限。
(4)旋轉(zhuǎn)備用約束
式中:Rup,h為h時(shí)段內(nèi)上調(diào)的備用容量;Rdown,h為h時(shí)段內(nèi)下調(diào)的備用容量。
(5)機(jī)組爬坡約束
式中:URi為機(jī)組爬坡速率。
(6)支路潮流約束
采用某海上風(fēng)電場2019年實(shí)際運(yùn)行數(shù)據(jù)進(jìn)行算例分析,以驗(yàn)證海上風(fēng)電并網(wǎng)調(diào)度模型的正確性和有效性,模型中相關(guān)參數(shù)取值可參考文獻(xiàn)[11]。根據(jù)風(fēng)電場實(shí)際運(yùn)行數(shù)據(jù),結(jié)合系統(tǒng)中各機(jī)組實(shí)際出力曲線,可以得到風(fēng)功率概率密度函數(shù);根據(jù)已構(gòu)建模型的目標(biāo)函數(shù)和約束條件,在MATLAB中仿真計(jì)算,對(duì)模型求解,同時(shí)計(jì)算出不同置信區(qū)間下的風(fēng)功率曲線,具體如圖2所示。
圖2 不同置信區(qū)間下的風(fēng)功率曲線
在圖2中,Q10~Q90分別表示不同置信區(qū)間下的風(fēng)功率曲線。從圖2可以看出,以1 h~4 h為例,風(fēng)功率實(shí)際值高于預(yù)測值,如果根據(jù)預(yù)測功率調(diào)度,則會(huì)增加棄風(fēng)電量,增大棄風(fēng)懲罰成本,同時(shí)為了滿足負(fù)荷需要,機(jī)組出力增加,發(fā)電成本也會(huì)增加,此時(shí)應(yīng)增加風(fēng)電出力,使風(fēng)電輸出功率處于Q60~Q70之間,盡可能降低棄風(fēng)懲罰成本和發(fā)電成本。在6 h~8 h,風(fēng)功率實(shí)際值低于預(yù)測值,如果直接根據(jù)預(yù)測功率調(diào)度,則需要啟動(dòng)系統(tǒng)備用容量,此時(shí)應(yīng)減小風(fēng)電出力,使風(fēng)電輸出功率調(diào)整至Q30~Q40之間,盡量減少使用系統(tǒng)備用容量,降低備用容量成本。因此,調(diào)度部門應(yīng)根據(jù)海上風(fēng)電預(yù)測功率偏差,及時(shí)調(diào)整海上風(fēng)力發(fā)電計(jì)劃,在提高風(fēng)電使用率的同時(shí)盡可能減少使用系統(tǒng)備用容量,降低系統(tǒng)發(fā)電成本和備用容量成本,從而降低綜合成本。
表1給出了海上風(fēng)電并網(wǎng)調(diào)度模型計(jì)算的不同置信區(qū)間下風(fēng)功率曲線各成本情況。由表1可知,Q90的海上風(fēng)電運(yùn)行成本、系統(tǒng)發(fā)電成本、備用容量成本和棄風(fēng)懲罰成本均最小,風(fēng)電輸出功率最大,風(fēng)電利用率最高,此時(shí)對(duì)應(yīng)的綜合成本為2694.1萬元。因此,電力部門在對(duì)含海上風(fēng)電的系統(tǒng)進(jìn)行并網(wǎng)調(diào)度時(shí),應(yīng)綜合考慮海上風(fēng)電運(yùn)行成本、系統(tǒng)發(fā)電成本、備用容量成本和棄風(fēng)懲罰成本等各項(xiàng)成本,在保障系統(tǒng)安全性和穩(wěn)定性的同時(shí),盡可能提高風(fēng)電利用率。
表1 不同置信區(qū)間下的風(fēng)功率曲線 單位:萬元
以海上風(fēng)電運(yùn)行成本、系統(tǒng)發(fā)電成本、備用容量成本和棄風(fēng)懲罰成本組成的綜合成本為目標(biāo)函數(shù),綜合考慮系統(tǒng)平衡、電網(wǎng)安全、機(jī)組運(yùn)行、旋轉(zhuǎn)備用、機(jī)組爬坡和支路潮流等各類約束條件,建立基于Petri網(wǎng)理論的海上風(fēng)電并網(wǎng)調(diào)度模型,并采用實(shí)際海上風(fēng)電場數(shù)據(jù)進(jìn)行算例分析,結(jié)果表明,Q90置信區(qū)間下調(diào)度模型的綜合成本為2694.1萬元,驗(yàn)證了模型的正確性和實(shí)用性。