宗曉英，張志鵬，黃大為

(東北電力大學(xué)電氣工程學(xué)院，吉林吉林132012)

隨著化石能源的日益減少和環(huán)境污染的加劇，傳統(tǒng)常規(guī)機(jī)組燃料燃燒所引起的環(huán)境污染問題已經(jīng)引起了大家的廣泛關(guān)注。由于風(fēng)能等新能源具有可再生、環(huán)保等特點(diǎn)，越來越受到人們的重視。其中風(fēng)力發(fā)電是發(fā)展最快、裝機(jī)容量最大的新能源，因此風(fēng)電成為緩解化石能源枯竭和環(huán)境污染的必然選擇，由于受外界環(huán)境的影響較大，風(fēng)力發(fā)電的不確定性給大規(guī)模風(fēng)電并網(wǎng)的電力系統(tǒng)調(diào)度帶來一系列新的問題［1－3］。

關(guān)于含有風(fēng)電的電力系統(tǒng)調(diào)度目前主要存在兩大問題:一是如何提高風(fēng)電出力預(yù)測(cè)的精確度;二是減少發(fā)電成本，協(xié)調(diào)風(fēng)電出力與煤耗量之間的關(guān)系。文獻(xiàn)［4］引入新的衡量電網(wǎng)安全的電網(wǎng)安全條件風(fēng)險(xiǎn)價(jià)值(CVaR)指標(biāo)，在不同置信水平下得到不同的條件風(fēng)險(xiǎn)調(diào)度方案。文獻(xiàn)［5］首先通過人工神經(jīng)網(wǎng)絡(luò)對(duì)風(fēng)向、風(fēng)速進(jìn)行了預(yù)測(cè)，然后根據(jù)功率曲線得到風(fēng)電出力預(yù)測(cè)。提出了一種基于預(yù)測(cè)誤差的非參數(shù)置信區(qū)間估計(jì)方法，建立了預(yù)測(cè)誤差概率模型，并在該模型的基礎(chǔ)上得到了風(fēng)電功率的預(yù)測(cè)值。文獻(xiàn)［6］引入風(fēng)速預(yù)測(cè)誤差的概率密度函數(shù)來處理風(fēng)電預(yù)測(cè)存在的誤差，在目標(biāo)函數(shù)中加入了誤差懲罰因子，并考慮了旋轉(zhuǎn)備用約束，建立了含風(fēng)電場(chǎng)的動(dòng)態(tài)預(yù)測(cè)模型。文獻(xiàn)［7］對(duì)于風(fēng)電出力的不確定性采用區(qū)間估計(jì)的方法，給出了有效備用響應(yīng)能力與區(qū)間滿足度約束的定義，但文中只局限于考慮單時(shí)段的調(diào)度。文獻(xiàn)［8］以區(qū)間數(shù)的形式對(duì)風(fēng)電出力不確定性進(jìn)行刻畫和建模，采用雙線性規(guī)劃模型得到區(qū)間優(yōu)化的兩個(gè)精確的上下邊界，基于空間分支定界法找到該模型的可行解也呈現(xiàn)區(qū)間形式。文獻(xiàn)［9］提出了一種考慮系統(tǒng)有功調(diào)節(jié)能力的風(fēng)電功率場(chǎng)景選取方法，將風(fēng)電出力按場(chǎng)景劃分，并對(duì)場(chǎng)景的選取進(jìn)行修正，考慮了多個(gè)時(shí)段內(nèi)的調(diào)度，建立了動(dòng)態(tài)經(jīng)濟(jì)調(diào)度模型。文獻(xiàn)［10］把風(fēng)功率視作負(fù)的負(fù)荷對(duì)含有風(fēng)電的系統(tǒng)進(jìn)行調(diào)度，采用區(qū)間估計(jì)的方法對(duì)風(fēng)電出力的不確定性進(jìn)行刻畫，并考慮了多時(shí)段下的動(dòng)態(tài)經(jīng)濟(jì)調(diào)度。傳統(tǒng)的電力系統(tǒng)調(diào)度方法已日臻完善，針對(duì)含有風(fēng)電出力的電力系統(tǒng)動(dòng)態(tài)經(jīng)濟(jì)調(diào)度的研究逐漸成為熱點(diǎn)。因此提高風(fēng)功率預(yù)測(cè)的精度，減少發(fā)電成本是所有含風(fēng)電機(jī)組的電力系統(tǒng)調(diào)度的主要研究問題。

本文在分析負(fù)荷和風(fēng)電預(yù)測(cè)的不確定性對(duì)含風(fēng)電機(jī)組的電力系統(tǒng)經(jīng)濟(jì)調(diào)度影響的基礎(chǔ)上，提出用三元區(qū)間數(shù)來刻畫負(fù)荷預(yù)測(cè)和風(fēng)電預(yù)測(cè)的不確定性，并建立基于三元區(qū)間數(shù)的電力系統(tǒng)動(dòng)態(tài)經(jīng)濟(jì)調(diào)度模型，針對(duì)凈負(fù)荷的三元區(qū)間對(duì)電力系統(tǒng)調(diào)度的經(jīng)濟(jì)性的影響進(jìn)行深入研究。

1 負(fù)荷與風(fēng)電功率的不確定性模型

負(fù)荷和風(fēng)電功率的不確定性是影響電力系統(tǒng)調(diào)度準(zhǔn)確性的最大影響因素，如何恰當(dāng)?shù)目坍嬤@些不確定量是提高電力系統(tǒng)調(diào)度準(zhǔn)確性的關(guān)鍵因素，風(fēng)速的變化可能不只服從一種概率分布［7］，本文假設(shè)風(fēng)電功率的變化服從正態(tài)分布，為了準(zhǔn)確的描述風(fēng)電功率的變化，采用三元區(qū)間數(shù)［11］的概念來刻畫負(fù)荷預(yù)測(cè)和風(fēng)功率預(yù)測(cè)的不確定性。

我們將負(fù)荷預(yù)測(cè)的三元區(qū)間定義為［Ll，Lm，Lu］，風(fēng)電預(yù)測(cè)的三元區(qū)間定義為［Wl，Wm，Wu］。風(fēng)力發(fā)電功率與系統(tǒng)負(fù)荷在很多方面有很強(qiáng)的正相關(guān)性，因此把風(fēng)力發(fā)電作為一個(gè)“負(fù)”的負(fù)荷，于是在風(fēng)電接入研究中用凈負(fù)荷當(dāng)作系統(tǒng)的總的負(fù)荷值［10］。凈負(fù)荷就是負(fù)荷預(yù)測(cè)值與風(fēng)電預(yù)測(cè)值的代數(shù)和，此時(shí)風(fēng)力發(fā)電出力當(dāng)作負(fù)的負(fù)荷。根據(jù)三元區(qū)間數(shù)的算術(shù)運(yùn)算法則［12］可求得凈負(fù)荷區(qū)間:

式中:Du為預(yù)測(cè)的凈負(fù)荷三元區(qū)間數(shù)的上限;Dl為預(yù)測(cè)的凈負(fù)荷三元區(qū)間數(shù)的下限;Dm為預(yù)測(cè)的凈負(fù)荷三元區(qū)間數(shù)的中間值。

負(fù)荷預(yù)測(cè)和風(fēng)電功率預(yù)測(cè)是在一定的區(qū)間內(nèi)上下波動(dòng)，因此我們給定一個(gè)置信度α:那么負(fù)荷預(yù)測(cè)和風(fēng)電功率預(yù)測(cè)的區(qū)間屬于那些包含三元區(qū)間數(shù)中間值Dm的區(qū)間的可信程度為1－α，如圖1。則:凈負(fù)荷預(yù)測(cè)區(qū)間的置信度滿足

圖1 凈負(fù)荷預(yù)測(cè)區(qū)間滿足的置信度

2 區(qū)間動(dòng)態(tài)經(jīng)濟(jì)調(diào)度的數(shù)學(xué)模型

如圖2，我們假設(shè)系統(tǒng)在t1時(shí)刻運(yùn)行在A點(diǎn)，由于凈負(fù)荷的不確定性影響，系統(tǒng)內(nèi)的某一機(jī)組可能在t2時(shí)刻運(yùn)行在B點(diǎn)，這時(shí)機(jī)組的發(fā)電出力需要在Δt的時(shí)間內(nèi)向上調(diào)節(jié)ΔPup;機(jī)組也可能運(yùn)行在C點(diǎn)，這時(shí)機(jī)組的發(fā)電出力需要在Δt的時(shí)間內(nèi)向下調(diào)節(jié)ΔPdown。在預(yù)測(cè)負(fù)荷的三元區(qū)間內(nèi)機(jī)組總是能夠快速反應(yīng)，保證了系統(tǒng)安全穩(wěn)定運(yùn)行，避免了因過負(fù)荷而導(dǎo)致的切機(jī)切負(fù)荷事件。

圖2 三元區(qū)間動(dòng)態(tài)經(jīng)濟(jì)調(diào)度方式

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

本文的主要研究對(duì)象是如何合理的安排常規(guī)機(jī)組的出力，以達(dá)到期望發(fā)電成本最小的目的。在三北地區(qū)，夜間用電低谷時(shí)期，往往是風(fēng)力大發(fā)時(shí)期，這時(shí)風(fēng)電機(jī)組出力較多，相應(yīng)的常規(guī)機(jī)組出力可以減小，體現(xiàn)在常規(guī)機(jī)組出力的權(quán)重下降;而在白天用電高峰期，風(fēng)速減小，風(fēng)電出力遠(yuǎn)不如夜間，此時(shí)需增加常規(guī)機(jī)組出力權(quán)重以滿足負(fù)荷變化的需求。因此目標(biāo)函數(shù)如下:

2.2 約束條件

(1)系統(tǒng)的供需平衡約束:

(2)常規(guī)機(jī)組的出力上下限約束:

式中:Rdi為常規(guī)機(jī)組i運(yùn)行時(shí)的最大爬坡速率;Δt為同一負(fù)荷水平下機(jī)組爬坡速率的反應(yīng)時(shí)間與不同負(fù)荷水平下機(jī)組爬坡速率的反應(yīng)時(shí)間。公式(10－12)表示在同一負(fù)荷水平下，即常規(guī)機(jī)組分別處于凈負(fù)荷預(yù)測(cè)的中間值、上限、下限三種情況下機(jī)組要滿足爬坡速率的約束;式(13－14)表示在不同的負(fù)荷水平下，即凈負(fù)荷在三元區(qū)間數(shù)內(nèi)變化時(shí)，機(jī)組的爬坡速率要滿足的約束。

(4)線路上的最大流動(dòng)功率約束:

輸電元件ij上允許流過的有功潮流分量Pij也是約束電力系統(tǒng)調(diào)度的一個(gè)因素，本文采用直流潮流這一近似計(jì)算來計(jì)算線路上實(shí)際流過的功率，并按照給定的裕度系數(shù)來確定線路上允許的最大流動(dòng)功率Pmaxij，并將Pmaxij作為已知條件應(yīng)用到常規(guī)機(jī)組出力上下限的約束中。

式中:δi為節(jié)點(diǎn)i的相角;Xij為線路i－j的電抗;Kij為線路上允許流過的功率的裕度系數(shù);Pij線路上實(shí)際流過的功率線路上允許流過的功率的最大值。

2.3 模型求解策略

3 算例分析

3.1 仿真參數(shù)

以IEEE－30節(jié)點(diǎn)系統(tǒng)為例，節(jié)點(diǎn)1為平衡節(jié)點(diǎn)，6臺(tái)常規(guī)機(jī)組。假設(shè)負(fù)荷和風(fēng)電功率服從N(35，152)的正態(tài)分布［8］，相應(yīng)的取α=0.05，Δt=1 h，T=6 h，Kij=2。風(fēng)電功率預(yù)測(cè)的三元區(qū)間數(shù)為［0.056，0.35，0.644］，負(fù)荷預(yù)測(cè)的三元區(qū)間數(shù)為［2.559，3，3.441］，凈負(fù)荷預(yù)測(cè)的三元區(qū)間數(shù)為［2.503，2.65，2.797］。發(fā)電機(jī)的參數(shù)見表1，權(quán)重系數(shù)的取值見表2。

表1 系統(tǒng)發(fā)電機(jī)參數(shù)

以上功率基準(zhǔn)值均為SB=100 MVA。

權(quán)重的物理意義是表述常規(guī)機(jī)組在凈負(fù)荷預(yù)測(cè)的三元區(qū)間數(shù)下的出力可能性的大小。當(dāng)風(fēng)電出力較大時(shí)，常規(guī)機(jī)組的出力需要下調(diào)，給風(fēng)電提供上網(wǎng)空間，常規(guī)機(jī)組出力由A調(diào)節(jié)到C的可能性增大(如圖2)，因此常規(guī)機(jī)組出力運(yùn)行在凈負(fù)荷預(yù)測(cè)三元區(qū)間數(shù)的下限的可能性增大，故β的值偏大。權(quán)重的取值考慮了風(fēng)電的出力特性，基于上述理論，按照逐一驗(yàn)證的方法，我們選取了以上三組權(quán)重系數(shù)。

表2 權(quán)重β的取值

3.2 仿真結(jié)果及分析

不同權(quán)重系數(shù)下，各機(jī)組的出力分配見圖3。

圖3 不同權(quán)重下三元區(qū)間數(shù)動(dòng)態(tài)經(jīng)濟(jì)調(diào)度

表3給出了動(dòng)態(tài)經(jīng)濟(jì)調(diào)度在不同權(quán)重下的期望成本，可以看出考慮了權(quán)重系數(shù)的調(diào)度方案期望成本較小。

通過圖表數(shù)據(jù)可以發(fā)現(xiàn)，在不同調(diào)度方式下，部分機(jī)組的出力是不變的，而有些機(jī)組出力變化較大。機(jī)組1，2，3在不同的權(quán)重和凈負(fù)荷下的出力變化很小，這些機(jī)組的經(jīng)濟(jì)性較差，可以考慮不參與含有風(fēng)電的經(jīng)濟(jì)調(diào)度，只承擔(dān)一些波動(dòng)較小的基礎(chǔ)負(fù)荷;機(jī)組4，5，6的爬坡能力強(qiáng)，機(jī)組出力變化較大，可以承擔(dān)波動(dòng)較大的負(fù)荷。同時(shí)考慮風(fēng)電出力的時(shí)間特性，在調(diào)度過程中增加權(quán)重系數(shù)可以減少發(fā)電費(fèi)用。通過對(duì)仿真數(shù)據(jù)的總結(jié)分析將不同出力特性的機(jī)組進(jìn)行分群，根據(jù)每個(gè)機(jī)組的出力特性安排機(jī)組的啟停和運(yùn)行時(shí)間，以達(dá)到減少發(fā)電費(fèi)用的目的。

表3 不同權(quán)重下三元區(qū)間數(shù)動(dòng)態(tài)經(jīng)濟(jì)調(diào)度的期望成本

4 結(jié)語

針對(duì)風(fēng)電功率不確定性對(duì)電網(wǎng)有功優(yōu)化調(diào)度的影響，建立了區(qū)間動(dòng)態(tài)經(jīng)濟(jì)調(diào)度模型。通過算例分析表明:由本文方法得到的區(qū)間動(dòng)態(tài)經(jīng)濟(jì)調(diào)度方案能確保在凈負(fù)荷預(yù)測(cè)的三元區(qū)間數(shù)內(nèi)滿足風(fēng)電和負(fù)荷的隨機(jī)變化，具有較好的適應(yīng)性;所得調(diào)度方案杜絕出現(xiàn)當(dāng)系統(tǒng)承受高負(fù)荷時(shí)因無法及時(shí)反應(yīng)而切機(jī)切負(fù)荷等極端情況，保證了供電可靠性，同時(shí)也確保機(jī)組安全平穩(wěn)的運(yùn)行。

［1］楊明，朱思萌，韓學(xué)山，等.風(fēng)電場(chǎng)輸出功率的多時(shí)段聯(lián)合概率密度預(yù)測(cè)［J］.電力系統(tǒng)自動(dòng)化，2013，10(7):23－28.

［2］高凱，朱加明，葛延峰，等.聯(lián)網(wǎng)風(fēng)電場(chǎng)集群運(yùn)行特性分析［J］.東北電力大學(xué)學(xué)報(bào)，2014，34(4):11－16.

［3］楊茂，賈云彭，錢為，等.基于動(dòng)態(tài)權(quán)重的風(fēng)電功率組合預(yù)測(cè)方法研究［J］.東北電力大學(xué)學(xué)報(bào)，2013，33(1/2):131－136.

［4］周任軍，姚龍華，童小嬌，等.采用條件風(fēng)險(xiǎn)方法的含風(fēng)電系統(tǒng)安全經(jīng)濟(jì)調(diào)度［J］.中國電機(jī)工程學(xué)報(bào)，2012，36(1):56－63.

［5］范高鋒，王偉勝，劉純，等.基于人工神經(jīng)網(wǎng)絡(luò)的風(fēng)電功率預(yù)測(cè)［J］.中國電機(jī)工程學(xué)報(bào)，2008，34(5):118－123.

［6］周瑋，孫輝，顧宏，馬千，等.計(jì)及風(fēng)險(xiǎn)備用約束的含風(fēng)電場(chǎng)電力系統(tǒng)動(dòng)態(tài)經(jīng)濟(jì)調(diào)度［J］.中國電機(jī)工程學(xué)報(bào)，2012，22(01):47－55.

［7］王士柏，韓學(xué)山，楊明.計(jì)及機(jī)組備用響應(yīng)能力的電力系統(tǒng)區(qū)間經(jīng)濟(jì)調(diào)度［J］.中國電機(jī)工程學(xué)報(bào)，2013，21(70):99－108.

［8］丁濤，郭慶來，柏瑞，等.考慮風(fēng)電不確定性的區(qū)間經(jīng)濟(jì)調(diào)度模型及空間分支定界法［J］.中國電機(jī)工程學(xué)報(bào)，2014，22(3):3707－3714.

［9］黃大為，張偉，韓學(xué)山.一種考慮系統(tǒng)有功調(diào)節(jié)能力的風(fēng)電功率場(chǎng)景選取方法［J］.電力系統(tǒng)保護(hù)與控制，2013，34(1):33－39.

［10］Kwok W.Cheung，Senior Member，Ricardo Rios－Zalapa.Smart Dispatch for Large Grid Operations with Integrated Renewable Resources［J］.IEEE Transactions on Energy Conversion，2011.

［11］胡啟洲，張衛(wèi)華.一種基于三元區(qū)間數(shù)的多指標(biāo)決策法［J］.統(tǒng)計(jì)與決策，2007，24(1):26－28.

［12］黃華，宋艷萍，德娜·吐熱汗.三元區(qū)間數(shù)線性規(guī)劃及其解法［J］.數(shù)學(xué)的實(shí)踐與認(rèn)識(shí)，2011，19(5):134－141.