考慮風(fēng)電功率預(yù)測不確定性的火電機(jī)組調(diào)度策略研究

趙征，于悅波

考慮風(fēng)電功率預(yù)測不確定性的火電機(jī)組調(diào)度策略研究

趙征，于悅波

（華北電力大學(xué) 控制與計(jì)算機(jī)工程學(xué)院，河北 保定 071003）

針對大規(guī)模風(fēng)電并網(wǎng)導(dǎo)致電力系統(tǒng)調(diào)度不確定性增大的問題，提出了一種考慮風(fēng)電功率不確定性的火電機(jī)組優(yōu)化調(diào)度方法。首先，采用均值聚類對歷史風(fēng)功率預(yù)測值和預(yù)測誤差進(jìn)行區(qū)間劃分，分區(qū)間采用高斯混合模型對預(yù)測誤差的概率密度函數(shù)進(jìn)行擬合。其次，建立考慮風(fēng)電預(yù)測不確定性的火電機(jī)組經(jīng)濟(jì)調(diào)度模型。最后，在麻雀算法的基礎(chǔ)上，結(jié)合天牛須尋優(yōu)算法進(jìn)行火電機(jī)組經(jīng)濟(jì)調(diào)度的求解。

大規(guī)模風(fēng)電并網(wǎng)；火電機(jī)組調(diào)度策略；麻雀優(yōu)化算法；均值聚類算法；高斯混合模型

0 引言

為實(shí)現(xiàn)節(jié)能減排，風(fēng)力發(fā)電技術(shù)得到了廣泛應(yīng)用，風(fēng)電裝機(jī)容量快速增加[1]。具有間歇性、隨機(jī)性的大規(guī)模風(fēng)電并網(wǎng)勢必會(huì)對電網(wǎng)的供電品質(zhì)甚至安全造成沖擊[2]。因此，研究考慮風(fēng)電不確定性的電力系統(tǒng)經(jīng)濟(jì)調(diào)度問題尤為重要。

目前，針對風(fēng)電功率預(yù)測不確定性的電力系統(tǒng)經(jīng)濟(jì)調(diào)度方法主要分為3類，分別是魯棒優(yōu)化法[3]、場景生成法[4]和以概率密度函數(shù)為基礎(chǔ)的隨機(jī)優(yōu)化法[5]。魯棒優(yōu)化法特點(diǎn)是不確定誤差分布，僅考慮誤差波動(dòng)范圍，要求約束條件在區(qū)間范圍內(nèi)最惡劣的情形下都要滿足，優(yōu)化結(jié)果過于保守。場景生成法雖然考慮了誤差的分布情況，但是需要生成大量場景后再對場景進(jìn)行縮減來選取典型場景，計(jì)算量巨大。基于概率密度函數(shù)的隨機(jī)優(yōu)化法是目前研究熱點(diǎn)。

文獻(xiàn)[5]用Beta分布對風(fēng)電功率預(yù)測誤差不確定性進(jìn)行描述，建立考慮風(fēng)電功率預(yù)測不確定性的電力系統(tǒng)經(jīng)濟(jì)調(diào)度模型。文獻(xiàn)[6]采用通用分布來描述風(fēng)電功率預(yù)測誤差的不確定性，通過對目標(biāo)函數(shù)和約束條件的轉(zhuǎn)化與分析，將隨機(jī)優(yōu)化模型轉(zhuǎn)化為一個(gè)非線性凸優(yōu)化問題，結(jié)合二次規(guī)劃算法和內(nèi)點(diǎn)法求解經(jīng)濟(jì)調(diào)度問題。文獻(xiàn)[7]將風(fēng)功率看作服從正態(tài)分布的隨機(jī)變量進(jìn)行經(jīng)濟(jì)調(diào)度，采用非支配排序遺傳算法實(shí)現(xiàn)混合整數(shù)多目標(biāo)優(yōu)化問題求解。以上單一概率分布模型無法同時(shí)描述風(fēng)電功率的非對稱、尖峰、多峰和厚尾特性。

文獻(xiàn)[8]針對短期風(fēng)電功率預(yù)測誤差不對稱甚至多峰的分布特性，提出采用高斯混合模型對預(yù)測誤差概率密度分布進(jìn)行擬合，取得了極好的擬合效果并應(yīng)用到了概率潮流的計(jì)算中。但是文獻(xiàn)[8]沒有考慮不同預(yù)測誤差區(qū)間下分布函數(shù)不同的情況。文獻(xiàn)[9]指出風(fēng)電功率預(yù)測值處于不同區(qū)間時(shí)誤差分布不同，并分區(qū)間進(jìn)行了概率密度函數(shù)擬合。文獻(xiàn)[10]分區(qū)間采用高斯混合模型進(jìn)行風(fēng)電功率不確定性建模并進(jìn)行火電機(jī)組經(jīng)濟(jì)調(diào)度研究，但是文獻(xiàn)[10]的區(qū)間劃分依然采用傳統(tǒng)經(jīng)驗(yàn)劃分方法，以預(yù)測值的0.3、0.6進(jìn)行劃分，區(qū)間劃分不準(zhǔn)確會(huì)對高斯混合模型參數(shù)產(chǎn)生影響進(jìn)而影響到風(fēng)電功率區(qū)間預(yù)測的準(zhǔn)確性；同時(shí)文獻(xiàn)[10]僅在約束條件中考慮風(fēng)電功率不確定性，目標(biāo)函數(shù)中未建立風(fēng)電功率不確定性成本。

本文采用均值聚類算法，首先對風(fēng)電功率預(yù)測區(qū)間進(jìn)行劃分；分區(qū)間采用高斯混合模型描述預(yù)測誤差的分布；建立考慮風(fēng)電功率不確定性的調(diào)度模型，在目標(biāo)函數(shù)中建立風(fēng)電功率不確定性成本；在高維度尋優(yōu)能力突出的麻雀算法的基礎(chǔ)上，提出結(jié)合天牛須尋優(yōu)算法彌補(bǔ)麻雀算法后期種群單一的問題；應(yīng)用改進(jìn)麻雀算法進(jìn)行電力系統(tǒng)經(jīng)濟(jì)調(diào)度的求解。

1 風(fēng)電功率不確定性建模

1.1 k均值聚類

均值聚類(-means)是以歐氏距離為基礎(chǔ)的聚類算法[11]，具有理論可靠、處理數(shù)據(jù)效率高、便于實(shí)現(xiàn)等優(yōu)點(diǎn)。

該方法能把數(shù)據(jù)集={1,2,···,x}?R×d劃分為個(gè)子集={c,=1,2,···,}，每個(gè)子集是一個(gè)類，聚類中心u為均值，即：

均值算法隨機(jī)選取聚類中心，再計(jì)算樣本與中心的歐氏距離，然后將樣本劃分到最近的中心所屬類別中。按照式（1）確定新的聚類中心，多次迭代后直至滿足收斂條件，最后得到聚類結(jié)果。

1.2 高斯混合模型

高斯混合模型結(jié)合參數(shù)和非參數(shù)方法的特點(diǎn)，可對任意隨機(jī)變量的概率分布函數(shù)進(jìn)行擬合，目前已經(jīng)在負(fù)荷和風(fēng)電功率預(yù)測上得到了很好的應(yīng)用。

高斯混合模型的核心思想是：任意隨機(jī)變量的概率分布都可以用一組正態(tài)分布來表示：

式中：whole為整體的概率密度函數(shù)；part,(,)為第個(gè)均值為、方差為的正態(tài)分布；為權(quán)重系數(shù)，其中?[1,]；為part,(,)的個(gè)數(shù)。受下式約束：

式（2）中的參數(shù)、通常通過EM算法迭代得到，EM算法具體步驟可見文獻(xiàn)[12]。

1.3 誤差分布擬合及評價(jià)

本文選取西北某風(fēng)電場的預(yù)測數(shù)據(jù)進(jìn)行誤差分析，分辨率為15 min。為對風(fēng)電預(yù)測誤差做進(jìn)一步研究，對風(fēng)電場風(fēng)功率預(yù)測誤差與風(fēng)功率預(yù)測值分布進(jìn)行統(tǒng)計(jì)；采用均值聚類方法按照聚類原則對分布情況進(jìn)行劃分，具體的劃分結(jié)果如圖1所示。

圖1 風(fēng)速–風(fēng)功率誤差聚類結(jié)果

由圖1可見，以預(yù)測值大小進(jìn)行劃分，將區(qū)間劃分為3份，2條分割線分別為1=83.43 MW、2=168.80 MW，對應(yīng)占比為0.304 6和0.637 8，區(qū)別于傳統(tǒng)經(jīng)驗(yàn)觀察法的0.3和0.6。在每一段區(qū)間上采用高斯混合模型GMM-3計(jì)算其相應(yīng)概率密度函數(shù)，不同區(qū)間的概率分布結(jié)果如圖2所示。誤差分布擬合參數(shù)如表1所示。限于文章篇幅，本文只列出預(yù)測誤差小于83.43 MW的誤差分布擬合參數(shù)。

觀察圖2可知，不同區(qū)間下的風(fēng)電功率預(yù)測誤差的分布函數(shù)不同。當(dāng)風(fēng)電功率預(yù)測值小于1時(shí)，誤差的概率密度函數(shù)是多峰瘦尾的；當(dāng)風(fēng)電功率預(yù)測值大于1且小于2時(shí)，風(fēng)電功率預(yù)測誤差的概率密度函數(shù)近似多峰厚尾；當(dāng)風(fēng)電功率預(yù)測值大于2時(shí)，風(fēng)電功率預(yù)測誤差的概率密度函數(shù)是單峰的。由此可見，分區(qū)間考慮誤差分布情況能更全面具體地描述誤差分布。

圖2 k均值聚類后不同區(qū)間的概率密度分布

表1 預(yù)測誤差小于83.43 MW誤差分布擬合參數(shù)

為驗(yàn)證采用均值聚類劃分區(qū)間相較于傳統(tǒng)經(jīng)驗(yàn)劃分區(qū)間的優(yōu)越性，分別采用均值聚類劃分區(qū)間和傳統(tǒng)經(jīng)驗(yàn)劃分區(qū)間得到兩組高斯混合模型，并基于兩組高斯混合模型分別進(jìn)行區(qū)間預(yù)測，評價(jià)指標(biāo)采用區(qū)間覆蓋率（PICP）及區(qū)間帶寬（PINAW），結(jié)果如表2和圖3所示。

表2 不同劃分方法對應(yīng)的區(qū)間預(yù)測評價(jià)指標(biāo)

由表2和圖3可知：采用均值聚類對區(qū)間進(jìn)行劃分，高斯混合模型區(qū)間預(yù)測PICP為0.927 1，區(qū)間帶寬為0.167 0；而采用傳統(tǒng)的經(jīng)驗(yàn)劃分，高斯混合模型區(qū)間預(yù)測PICP為0.916 7，區(qū)間帶寬為0.172 9。兩者比較可以發(fā)現(xiàn)，利用均值聚類對區(qū)間進(jìn)行劃分可以在區(qū)間覆蓋率提高1.04%的同時(shí)，區(qū)間帶寬也能降低3.41%。這表明采用均值聚類劃分區(qū)間能夠有效提高風(fēng)電功率區(qū)間預(yù)測效果。

圖3 兩種區(qū)間劃分方式下區(qū)間預(yù)測結(jié)果

2 考慮風(fēng)電功率不確定性的火電機(jī)組經(jīng)濟(jì)調(diào)度模型

2.1 模型構(gòu)建

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

本文的調(diào)度模型中，將風(fēng)功率預(yù)測誤差作為隨機(jī)變量，以風(fēng)電功率預(yù)測值為確定變量，以每臺火電機(jī)組的出力情況為決策變量，以總發(fā)電成本最低為目標(biāo)函數(shù)，同時(shí)考慮滿足系統(tǒng)運(yùn)行的相關(guān)約束。

目標(biāo)函數(shù)包括3部分：火電機(jī)組燃料成本、火電機(jī)組閥點(diǎn)效應(yīng)成本、風(fēng)電功率預(yù)測不確定性成本；不考慮風(fēng)電場運(yùn)行成本，認(rèn)為風(fēng)電場歸電網(wǎng)所有，除固定的投資成本和維修成本外無燃料成本[10,13]。具體表達(dá)式為：

2.1.2 約束條件

（1）功率平衡約束

式中：W為風(fēng)電功率預(yù)測值；L為負(fù)荷需求預(yù)測值。

（2）火電機(jī)組爬坡速率約束

（3）火電機(jī)組出力上下限約束

（4）火電機(jī)組旋轉(zhuǎn)備用容量約束

同理，向下旋轉(zhuǎn)備用為：

2.2 不確定性模型轉(zhuǎn)化為確定性模型

在考慮風(fēng)電功率預(yù)測的火電機(jī)組經(jīng)濟(jì)調(diào)度模型中，風(fēng)電功率預(yù)測誤差為隨機(jī)變量，因此在調(diào)度過程中需要將不確定的隨機(jī)變量轉(zhuǎn)換為確定的量。通過風(fēng)電功率預(yù)測誤差的概率密度函數(shù)，可以經(jīng)過如下方式，將包含不確定隨機(jī)變量的式（10）（11）進(jìn)行轉(zhuǎn)換。

式中：(·)為風(fēng)功率誤差的概率分布函數(shù)。

2.3 可行化調(diào)度方法

應(yīng)用尋優(yōu)方法求解多約束、高維度的火電機(jī)組經(jīng)濟(jì)調(diào)度模型時(shí)，如果將所有的約束條件都作為懲罰項(xiàng)加入到目標(biāo)函數(shù)中，會(huì)導(dǎo)致算法容易陷入局部極值，而且無法保證滿足所有約束。本文在傳統(tǒng)懲罰函數(shù)法[15]的基礎(chǔ)上做“減法”，利用文獻(xiàn)[16]中的可行化調(diào)整方法，將目標(biāo)函數(shù)中的功率平衡約束、火電機(jī)組出力上下限約束移到目標(biāo)函數(shù)之外進(jìn)行處理，在目標(biāo)函數(shù)中只保留旋轉(zhuǎn)備用約束和爬坡速率約束。這樣避免了算法陷入局部極值的情況，有效緩解了高維度、多約束的火電機(jī)組經(jīng)濟(jì)調(diào)度模型難以獲得有效調(diào)度結(jié)果的困難。具體的可行化調(diào)整方法見文獻(xiàn)[16]。

3 改進(jìn)麻雀算法

3.1 麻雀優(yōu)化算法基本原理

麻雀算法（SSA）是文獻(xiàn)[17]根據(jù)麻雀覓食和反捕食行為提出的一種新型優(yōu)化算法。算法將麻雀覓食過程抽象為發(fā)現(xiàn)者–加入者模型并在這一模型中加入偵察者機(jī)制。

算法中，每只麻雀的位置對應(yīng)一個(gè)解。在維空間中，麻雀數(shù)量為，則第（≤）只麻雀在此搜索空間中的位置可以為=[x1,x2,···,x,···,x]。其中，?[1,]；?[1,]；x表示第只麻雀在第維的位置。發(fā)現(xiàn)者位置更新如式（12）所示：

式中：為當(dāng)前迭代次數(shù)；為最大迭代次數(shù)；是(0, 1)之間的隨機(jī)數(shù)；2?[0,1]為預(yù)警值，?[0.5, 1]為安全值。是一個(gè)服從標(biāo)準(zhǔn)正態(tài)分布的隨機(jī)數(shù)；是一個(gè)1×的矩陣，每個(gè)元素均為1。

加入者位置更新如式（13）所示：

式中：為一個(gè)隨機(jī)數(shù)，其服從均值為0、方差為1的標(biāo)準(zhǔn)正態(tài)分布；為一個(gè)隨機(jī)數(shù)，?[–1,1]。i、g、w分別表示麻雀的適應(yīng)度、麻雀種群的最優(yōu)適應(yīng)度、麻雀種群的最差適應(yīng)度。

3.2 改進(jìn)麻雀算法

麻雀算法的加入者在全維度向最佳位置靠近時(shí)，雖然能達(dá)到快速收斂的效果，但是也減少了種群多樣性，使算法容易陷入局部最優(yōu)[18]。

天牛須算法（BAS）是一種個(gè)體啟發(fā)類算法。由于個(gè)體位置更新不受群體影響，因此將天牛須算法與麻雀算法結(jié)合可以克服加入者后期向最佳位置靠近時(shí)種群多樣性減少的缺點(diǎn)。種群多樣性的增加，增大了算法跳出局部最優(yōu)解的機(jī)會(huì)。同時(shí)，為了保持麻雀算法本身的特性，做出以下改進(jìn)，如圖4所示。

圖4 改進(jìn)麻雀算法流程圖

在加入者中選擇1只加入者。對于這1只加入者，將其位置更新公式（12）修改為天牛須[19]的更新方式：

其中，由式（16）確定：

3.3 改進(jìn)麻雀算法有效性驗(yàn)證

為驗(yàn)證改進(jìn)麻雀算法在高維尋優(yōu)上的有效性，采用標(biāo)準(zhǔn)測試函數(shù)對改進(jìn)前后的麻雀算法進(jìn)行測試：麻雀算法的種群個(gè)數(shù)為50；迭代次數(shù)為500；發(fā)現(xiàn)者、加入者的比例為7:3。

采用不同類型的測試函數(shù)來驗(yàn)證算法的改進(jìn)效果。測試函數(shù)如表3所示[13]。表中，高維單峰函數(shù)為1、2，高維多峰函數(shù)為3、4。對每個(gè)測試函數(shù)測試10次，尋優(yōu)結(jié)果的均值和方差如表4所示。

表3 測試函數(shù)表

注：表中表示維度，=30。

表4 測試結(jié)果

由表3可知，對于高維單峰函數(shù)1、2而言，改進(jìn)麻雀算法能夠準(zhǔn)確找到最優(yōu)值，同時(shí)方差也為0：這表征改進(jìn)麻雀算法能找到最優(yōu)值，同時(shí)也具有極高的穩(wěn)定性。對于高峰多維函數(shù)4，其在改進(jìn)前后都有良好的準(zhǔn)確性和穩(wěn)定性，這表征麻雀算法在高維尋優(yōu)問題上具有突出尋優(yōu)能力。對于函數(shù)3，傳統(tǒng)麻雀算法的尋優(yōu)結(jié)果與最優(yōu)值相差甚大；而改進(jìn)麻雀算法的尋優(yōu)結(jié)果已經(jīng)與最優(yōu)值非常接近，且改進(jìn)麻雀算法的方差也得到了大幅度降低：這表征改進(jìn)麻雀算法在尋優(yōu)精度和準(zhǔn)確性上有大幅度提升。

為細(xì)化對函數(shù)3的尋優(yōu)過程，分別采用傳統(tǒng)麻雀算法和改進(jìn)麻雀算法對3尋優(yōu)并畫出目標(biāo)值與迭代次數(shù)的關(guān)系圖，結(jié)果如圖5和圖6所示。

對于3，傳統(tǒng)麻雀算法已經(jīng)完全陷入了局部極值，無法找到準(zhǔn)確的最小值；而改進(jìn)后的算法則可以準(zhǔn)確找到全局最小值點(diǎn)。尋優(yōu)過程如圖5、圖6所示。

圖5 原始麻雀算法對F3尋優(yōu)過程

圖6 改進(jìn)麻雀算法對F3尋優(yōu)過程

觀察圖5、圖6可以發(fā)現(xiàn)，對于高維多峰的函數(shù)3，改進(jìn)的麻雀算法雖然在153次迭代前陷入了極值，但是在第154次迭代時(shí)能夠跳出局部極值找到真實(shí)的最小值。傳統(tǒng)麻雀算法在迭代到500次的時(shí)候依然無法找到最小值，已陷入局部極值。這表明改進(jìn)后的麻雀算法具備良好的種群多樣性，在算法后期依然有跳出局部極值的能力。

4 算例分析

本文算例的風(fēng)電數(shù)據(jù)來自于某風(fēng)電場實(shí)際運(yùn)行數(shù)據(jù)，火電機(jī)組為10臺，具體的火電機(jī)組參數(shù)可參考文獻(xiàn)[20]。算例中：風(fēng)電容量與火電容量之比為2:8，滿足大規(guī)模風(fēng)電并網(wǎng)條件。調(diào)度周期為1 h；調(diào)度模型中置信度1、2均取值0.95；旋轉(zhuǎn)備用響應(yīng)時(shí)間10為10 min；改進(jìn)麻雀算法的發(fā)現(xiàn)者、加入者的比例為7:3。應(yīng)用改進(jìn)麻雀算法進(jìn)行火電機(jī)組調(diào)度，計(jì)算結(jié)果如圖7所示。

觀察圖7可以發(fā)現(xiàn)，通過本文的調(diào)度策略安排火電機(jī)組的未來出力，在大規(guī)模風(fēng)電入網(wǎng)的前提下，仍能滿足功率平衡約束及其他各項(xiàng)約束，調(diào)度結(jié)果具有有效性。為研究改進(jìn)麻雀算法（SSA-BAS）及均值區(qū)間劃分對功率分配（ED）中發(fā)電成本的影響，與粒子群算法（PSO）、灰狼算法（GWO）進(jìn)行比較，分6個(gè)場景研究了火電機(jī)組的發(fā)電成本情況。具體的場景及發(fā)電成本如表5所示。

圖7 SSA-ED-k-means劃分區(qū)間調(diào)度結(jié)果

表5 不同場景下預(yù)測效果評價(jià)

注：發(fā)電成本為0:00—24:00內(nèi)火電機(jī)組燃料成本、火電機(jī)組閥點(diǎn)效應(yīng)成本、風(fēng)電功率預(yù)測不確定性成本之和

通過表5，對比場景1和場景2可以發(fā)現(xiàn)，改進(jìn)后的麻雀算法可以有效降低發(fā)電成本；對比場景2和場景3可知，對預(yù)測區(qū)間的劃分可以有效降低發(fā)電成本；對比場景3和場景4可以發(fā)現(xiàn)，使用-means聚類對區(qū)間劃分比經(jīng)驗(yàn)劃分能更多地降低發(fā)電成本。對比場景4～6可以發(fā)現(xiàn)，此實(shí)驗(yàn)背景下，本文所提的改進(jìn)麻雀算法在求解高維度多約束的火電機(jī)組經(jīng)濟(jì)調(diào)度問題上，較灰狼算法與粒子群算法而言具有優(yōu)越性，發(fā)電成本最低。

為進(jìn)一步研究置信度1、2對火電機(jī)組調(diào)度的影響，表6給出了不同置信度下火電機(jī)組旋轉(zhuǎn)備用容量的變化情況。

表6 不同置信度下火電機(jī)組旋轉(zhuǎn)備用容量的變化情況

觀察表6可以發(fā)現(xiàn)：當(dāng)置信度1自0.80增加到0.95而置信度2保持0.95不變的情況下，火電機(jī)組提供的上升旋轉(zhuǎn)備用隨置信度1的增加而增大，從1 155.7 MW增加到1 264.4 MW。當(dāng)置信度1保持0.95不變而置信度2自0.95減小到0.80時(shí)，火電機(jī)組的下降旋轉(zhuǎn)備用自1 079.3 MW減少到了861.571 2 MW。這表征，置信度水平的增加雖然降低系統(tǒng)風(fēng)險(xiǎn)，但是會(huì)導(dǎo)致火電機(jī)組的旋轉(zhuǎn)備用容量增加。

為研究火電機(jī)組變負(fù)荷速率對火電機(jī)組發(fā)電成本以及對風(fēng)電消納容量的影響，分別做火電機(jī)組在原爬坡速率下、增加10%爬坡速率、增加20%爬坡速率、增加30%爬坡速率場景下的火電機(jī)組調(diào)度策略安排，結(jié)果如表7和表8所示。

表7 火電機(jī)組爬坡速率提升對發(fā)電成本的影響

注：見表5注。

根據(jù)表7可以明顯發(fā)現(xiàn)，火電機(jī)組爬坡速率的提高會(huì)降低火電機(jī)組的發(fā)電成本。分析其原因?yàn)椋号榔滤俾实奶岣邥?huì)導(dǎo)致受爬坡速率約束式（6）約束的搜索空間變大。

在火電機(jī)組原爬坡速率下、增加10%爬坡速率、增加20%爬坡速率、增加30%爬坡速率4種場景下，不斷增加風(fēng)電并網(wǎng)容量，直到調(diào)度安排不滿足某一項(xiàng)約束條件為止。以能滿足各項(xiàng)約束的最大風(fēng)電容量來研究火電機(jī)組變負(fù)荷速率對風(fēng)電消納容量的影響。不同爬坡速率增長情況下，對應(yīng)的多消納的風(fēng)電容量結(jié)果如表8所示。

表8 火電機(jī)組爬坡速率提升對風(fēng)電消納容量的影響

由表8可知，當(dāng)火電機(jī)組爬坡速率增長時(shí)，系統(tǒng)可以有效增加風(fēng)電消納。對于本文討論的系統(tǒng)而言，當(dāng)爬坡速率增加30%的時(shí)候，可以多消納風(fēng)電49.86 MW。因此，在實(shí)際應(yīng)用中，應(yīng)增加火電機(jī)組爬坡速率以有效消納更多風(fēng)電資源。

5 結(jié)論

（1）較之傳統(tǒng)劃分方法，采用均值聚類算法對預(yù)測值進(jìn)行劃分，能夠更準(zhǔn)確地描述誤差的分布，使風(fēng)電功率區(qū)間預(yù)測結(jié)果更可靠。

（2）改進(jìn)的麻雀算法能夠提高傳統(tǒng)麻雀算法在高維度的尋優(yōu)能力。算法在后期依然有較高的種群豐富度，進(jìn)而增大了后期跳出局部最優(yōu)值的能力。同時(shí)，本文所提的改進(jìn)麻雀算法在求解高維度多約束的火電機(jī)組經(jīng)濟(jì)調(diào)度問題上，較灰狼算法與粒子群算法更具有優(yōu)越性。

火電機(jī)組爬坡速率的提高，能夠有效降低發(fā)電成本、增加對風(fēng)電資源的消納。

[1] 趙征, 于悅波, 孫昊天. 基于凝結(jié)水節(jié)流的新型協(xié)調(diào)優(yōu)化控制策略[J]. 動(dòng)力工程學(xué)報(bào), 2021, 41(2): 107-112.

ZHAO ZHENG, YU YUEBO, SUN HAOTIAN. Optimization of a new coordinated control strategy based on condensate throttling[J]. Journal of Chinese Society of Power Engineering, 2021, 41(2): 107-112(in Chinese).

[2] 楊錫運(yùn), 邢國通, 馬雪, 等. 一種核極限學(xué)習(xí)機(jī)分位數(shù)回歸模型及風(fēng)電功率區(qū)間預(yù)測[J]. 太陽能學(xué)報(bào), 2020, 41(11): 300-306.

YANG XIYUN, XING GUOTONG, MA XUE, et al. A model of quantile regression with kernel extreme learning machine and wind power interval prediction[J]. Acta Energiae Solaris Sinica, 2020, 41(11): 300-306(in Chinese).

[3] 崔明勇, 靳賀, 何良策, 等. 基于風(fēng)電消納評估的電力系統(tǒng)魯棒低碳經(jīng)濟(jì)調(diào)度[J]. 太陽能學(xué)報(bào), 2020, 41(8): 270-280.

CUI MINGYONG, JIN HE, HE LIANGCE, et al. Robust low-carbon economic dispatch of power system based on wind power accommodation evaluation[J]. Acta Energiae Solaris Sinica, 2020, 41(8): 270-280(in Chinese).

[4] 伍棟文, 于艾清. 壞場景集下含風(fēng)電多源系統(tǒng)的改進(jìn)魯棒調(diào)度模型[J]. 太陽能學(xué)報(bào), 2018, 39(10): 2929-2936.

WU DONGWEN, YU AIQING. Improved robust scheduling model for multi-source power system with wind power based on bad-scenario set[J]. Acta Energiae Solaris Sinica, 2018, 39(10): 2929-2936(in Chinese).

[5] 吳棟梁, 王揚(yáng), 郭創(chuàng)新, 等. 電力市場環(huán)境下考慮風(fēng)電預(yù)測誤差的經(jīng)濟(jì)調(diào)度模型[J]. 電力系統(tǒng)自動(dòng)化, 2012, 36(6): 23-28.

WU DONGLIANG, WANG YANG, GUO CHUANGXIN, et al. An economic dispatching model considering wind power forecast errors in electricity market environment[J]. Automation of Electric Power Systems, 2012, 36(6): 23-28(in Chinese).

[6] 王豹, 徐箭, 孫元章, 等. 基于通用分布的含風(fēng)電電力系統(tǒng)隨機(jī)動(dòng)態(tài)經(jīng)濟(jì)調(diào)度[J]. 電力系統(tǒng)自動(dòng)化, 2016, 40(6): 17-24.

WANG BAO, XU JIAN, SUN YUANZHANG, et al. Stochastic dynamic economic dispatch of power systems considering wind power based on versatile probability distribution[J]. Automation of Electric Power Systems, 2016, 40(6): 17-24(in Chinese).

[7] 馬燕峰, 范振亞, 劉偉東, 等. 考慮碳權(quán)交易和風(fēng)荷預(yù)測誤差隨機(jī)性的環(huán)境經(jīng)濟(jì)調(diào)度[J]. 電網(wǎng)技術(shù), 2016, 40(2): 412-418.

MA YANFENG, FAN ZHENYA, LIU WEIDONG, et al. Environmental and economic dispatch considering carbon trading credit and randomicity of wind power and load forecast error[J]. Power System Technology, 2016, 40(2): 412-418(in Chinese).

[8] 葉林, 張亞麗, 巨云濤, 等. 用于含風(fēng)電場的電力系統(tǒng)概率潮流計(jì)算的高斯混合模型[J]. 中國電機(jī)工程學(xué)報(bào), 2017, 37(15): 4379-4387.

YE LIN, ZHANG YALI, JU YUNTAO, et al. Gaussian mixture model for probabilistic power flow calculation of system integrated wind farm[J]. Proceedings of the CSEE, 2017, 37(15): 4379-4387(in Chinese).

[9] 趙書強(qiáng), 王揚(yáng), 徐巖. 基于風(fēng)電預(yù)測誤差隨機(jī)性的火儲聯(lián)合相關(guān)機(jī)會(huì)規(guī)劃調(diào)度[J]. 中國電機(jī)工程學(xué)報(bào), 2014, 34(S1): 9-16.

ZHAO SHUQIANG, WANG YANG, XU YAN. Dependent chance programming dispatching of integrated thermal power generation and energy storage system based on wind power forecasting error[J]. Proceedings of the CSEE, 2014, 34(S1): 9-16(in Chinese).

[10] 易明月, 童曉陽. 考慮風(fēng)荷預(yù)測誤差不確定性的動(dòng)態(tài)經(jīng)濟(jì)調(diào)度[J]. 電網(wǎng)技術(shù), 2019, 43(11): 4050-4056.

YI MINGYUE, TONG XIAOYANG. Dynamic economic dispatch considering uncertainties of wind power and load forecast error[J]. Power System Technology, 2019, 43(11): 4050-4056(in Chinese).

[11] 黃雨薇, 彭道剛, 姚峻, 等. 基于SSA和K均值的TD-BP神經(jīng)網(wǎng)絡(luò)超短期光伏功率預(yù)測[J]. 太陽能學(xué)報(bào), 2021, 42(4): 229-238.

HUANG YUWEI, PENG DAOGANG, YAO JUN, et al. Ultra-short-team photovoltaic power forecast of TD-BP neural network based on SSA and K-means[J]. Acta Energiae Solaris Sinica, 2021, 42(4): 229-238(in Chinese).

[12] 楊茂, 董駿城. 基于混合高斯分布的風(fēng)電功率實(shí)時(shí)預(yù)測誤差分析[J]. 太陽能學(xué)報(bào), 2016, 37(6): 1594-1602.

YANG MAO, DONG JUNCHENG. Real-time prediction error analysis of wind power based on mixed Gaussian distribution model[J]. Acta Energiae Solaris Sinica, 2016, 37(6): 1594-1602(in Chinese).

[13] 周瑋, 孫輝, 顧宏, 等. 計(jì)及風(fēng)險(xiǎn)備用約束的含風(fēng)電場電力系統(tǒng)動(dòng)態(tài)經(jīng)濟(jì)調(diào)度[J]. 中國電機(jī)工程學(xué)報(bào), 2012, 32(1): 47-55.

ZHOU WEI, SUN HUI, GU HONG, et al. Dynamic economic dispatch of wind integrated power systems based on risk reserve constraints[J]. Proceedings of the CSEE, 2012, 32(1): 47-55(in Chinese).

[14] 國家能源局華北監(jiān)管局. 華北區(qū)域發(fā)電廠并網(wǎng)運(yùn)行管理實(shí)施細(xì)則(試行)[Z]. 北京: 華北電監(jiān)局, 2019.

[15] 翁振星, 石立寶, 徐政, 等. 計(jì)及風(fēng)電成本的電力系統(tǒng)動(dòng)態(tài)經(jīng)濟(jì)調(diào)度[J]. 中國電機(jī)工程學(xué)報(bào), 2014, 34(4): 514-523.

WENG ZHENXING, SHI LIBAO, XU ZHENG, et al. Power system dynamic economic dispatch incorporating wind power cost[J]. Proceedings of the CSEE, 2014, 34(4): 514-523(in Chinese).

[16] 張海峰, 高峰, 吳江, 等. 含風(fēng)電的電力系統(tǒng)動(dòng)態(tài)經(jīng)濟(jì)調(diào)度模型[J]. 電網(wǎng)技術(shù), 2013, 37(5): 1298-1303.

ZHANG HAIFENG, GAO FENG, WU JIANG, et al. A dynamic economic dispatching model for power grid containing wind power generation system[J]. Power Systems Technology, 2013, 37(5): 1298-1303(in Chinese).

[17] XUE J, SHEN B. A novel swarm intelligence optimization approach: sparrow search algorithm[J]. Systems Science & Control Engineering, 2020, 8(1): 22-34.

[18] 呂鑫, 慕曉冬, 張鈞. 基于改進(jìn)麻雀搜索算法的多閾值圖像分割[J]. 系統(tǒng)工程與電子技術(shù), 2021, 43(2): 318-327.

LV XIN, MU XIAODONG, ZHANG JUN, et al. Multi-threshold image segmentation based on improved sparrow search algorithm[J]. Systems Engineering and Electronics, 2021, 43(2): 318-327(in Chinese).

[19] 喬貴方, 呂仲艷, 張穎, 等. 基于BAS-PSO算法的機(jī)器人定位精度提升[J]. 光學(xué)精密工程, 2021, 29(4): 763-771.

QIAO GUIFANG, LV ZHONGYAN, ZHANG YING, et al. Improvement of robot kinematic accuracy based on BAS-PSO algorithm[J]. Optics and Precision Engineering, 2021, 29(4): 763-771(in Chinese).

[20] 梁吉, 左藝, 張玉琢, 等. 基于可再生能源配額制的風(fēng)電并網(wǎng)節(jié)能經(jīng)濟(jì)調(diào)度[J]. 電網(wǎng)技術(shù), 2019, 43(7): 2528-2534.

LIANG JI, ZUO YI, ZHANG YUZHUO, et al. Energy-saving and economic dispatch of power system containing wind power integration under renewable portfolio standard[J]. Power System Technology, 2019, 43(7): 2528-2534(in Chinese).

Research on Dispatching Strategy of Thermal Power Units Considering the Uncertainty of Wind Power Forecast

ZHAO Zheng, YU Yuebo

(School of Control and Computer Engineering, North China Electric Power University, Baoding 071003, China)

To solve the problem of due to the increased uncertainty in power system dispatch, an optimal dispatch method of thermal power units considering the uncertainty of wind power is proposed. First, this paper uses-means clustering to divide the historical wind power prediction value from the prediction error, and the Gaussian mixture model to fit the probability density function of the prediction error between the partitions. Secondly, it establishes an economic dispatch model for thermal power units that takes into account the forecast errors of wind power. Finally, based on the basis of the Sparrow algorithm, it gives the optimal solution of economic dispatch of thermal power units with the Tianniu optimization algorithm.

large-scale wind power grid connection; dispatching strategy of thermal power unit; sparrow optimization algorithm;-means clustering algorithm; Gaussian mixture model

10.3969/j.ISSN.1672-0792.2021.11.004

TM614; TM734

A

1672-0792(2021)11-0032-10

2021-08-14

中央高校基本科研業(yè)務(wù)費(fèi)（2017MS133）；北京市自然科學(xué)基金（3202027）

趙 征（1978—），女，副教授，主要研究方向?yàn)轱L(fēng)速、風(fēng)功率預(yù)測，熱工過程建模與控制優(yōu)化；

于悅波（1996—），男，碩士研究生，主要研究方向?yàn)轱L(fēng)速、風(fēng)功率預(yù)測，電力系統(tǒng)優(yōu)化調(diào)度。

于悅波