范劉洋, 汪可友, 李國杰, 吳 巍, 葛維春
(1. 電力傳輸與功率變換控制教育部重點(diǎn)實(shí)驗(yàn)室(上海交通大學(xué)), 上海市 200240;2. 國網(wǎng)遼寧省電力有限公司, 遼寧省沈陽市 110006)
近年來,隨著技術(shù)的成熟和成本的降低,風(fēng)電、太陽能等可再生能源裝機(jī)容量得到了快速的提高[1]。國家風(fēng)電發(fā)展“十二五”規(guī)劃中提出,2020年中國風(fēng)電裝機(jī)規(guī)模將達(dá)到2×108kW,風(fēng)電裝機(jī)容量將占總裝機(jī)容量的11%[2]。風(fēng)電具有間歇性、波動(dòng)性等特點(diǎn),日前預(yù)測誤差可達(dá)在20%,遠(yuǎn)高于負(fù)荷預(yù)測誤差[3]。大規(guī)模風(fēng)電的接入,給電力系統(tǒng)調(diào)度帶來了巨大挑戰(zhàn)。
近年來,魯棒優(yōu)化理論被應(yīng)用于這一領(lǐng)域,并取得良好效果[4-9]。魯棒優(yōu)化使用閉凸集合描述參數(shù)的不確定性,并計(jì)算目標(biāo)函數(shù)在“最壞情況”下的最優(yōu)問題[10-13]。相比于隨機(jī)規(guī)劃,魯棒優(yōu)化不要求獲得隨機(jī)變量分布函數(shù),不需要生成典型場景,更易于實(shí)施。但另一方面,魯棒優(yōu)化結(jié)果一般保守性較強(qiáng)。
魯棒優(yōu)化結(jié)果的保守性直接受到不確定參數(shù)集合的影響。最早,文獻(xiàn)[10]采用每個(gè)不確定變量在一定置信水平下的上下界作為不確定集。近年來,文獻(xiàn)[11-12]采用橢球集合描述不確定集,文獻(xiàn)[4-5,7]采用盒式集合描述不確定集,進(jìn)一步降低了魯棒優(yōu)化的保守性,由于盒式集合的線性性質(zhì),這種不確定集合在電力系統(tǒng)中應(yīng)用更為廣泛。
但是,以上不確定集合均假定每個(gè)時(shí)刻不確定參數(shù)服從的是獨(dú)立分布律,這并不符合實(shí)際情況。以含風(fēng)電的調(diào)度問題為例,當(dāng)風(fēng)電場之間距離較近時(shí),兩個(gè)或多個(gè)風(fēng)電場之間的出力會(huì)具有相關(guān)性。這種現(xiàn)象已經(jīng)在風(fēng)電場建設(shè)規(guī)劃、風(fēng)險(xiǎn)機(jī)組組合、隨機(jī)規(guī)劃場景生成等問題中被注意到[14-17]。
然而對(duì)風(fēng)電預(yù)測誤差時(shí)間相關(guān)性(以下簡稱“時(shí)間相關(guān)性”)的研究仍然較少。電力系統(tǒng)在制定機(jī)組組合計(jì)劃時(shí),需要的是提前一天到一周內(nèi)的風(fēng)速預(yù)測信息。這一時(shí)間尺度的風(fēng)速預(yù)測往往使用數(shù)值天氣預(yù)報(bào)(numeric weather predictor,NWP)模型[18]。該方法根據(jù)風(fēng)場周圍的氣壓、氣溫、等高線以及障礙物等物理信息構(gòu)建風(fēng)速的微分方程表達(dá)式,對(duì)風(fēng)速進(jìn)行預(yù)測。文獻(xiàn)[18-19]根據(jù)大量的風(fēng)電場運(yùn)行歷史數(shù)據(jù),建立了1 h步長的一階風(fēng)速馬爾可夫鏈模型。其中,狀態(tài)轉(zhuǎn)移概率矩陣呈近似對(duì)角分布,對(duì)角線及其附近元素?cái)?shù)值較大,這從側(cè)面驗(yàn)證了風(fēng)速在時(shí)間上的相關(guān)特點(diǎn)。文獻(xiàn)[20]分析了愛爾蘭風(fēng)電場的風(fēng)電數(shù)據(jù),直接指出風(fēng)電預(yù)測誤差序列具有明確的自相關(guān)性特性,并提出利用該特性降低風(fēng)電預(yù)測誤差的統(tǒng)計(jì)方法。文獻(xiàn)[21]將風(fēng)電預(yù)測誤差的協(xié)方差矩陣加入高斯分布中,并驗(yàn)證考慮相關(guān)性的模型可以提高魯棒最優(yōu)潮流模型的經(jīng)濟(jì)性。在魯棒機(jī)組組合(RUC)問題中,理論上同樣可以利用時(shí)間相關(guān)性特點(diǎn),剔除一些概率密度較低的場景,降低不確定集合的保守性。但是魯棒優(yōu)化本質(zhì)為多層優(yōu)化問題,對(duì)約束條件要求較高,統(tǒng)計(jì)方法的分析結(jié)果需要表示為線性約束才能應(yīng)用于魯棒優(yōu)化算法中。
本文首先以德國巴登—符騰堡州風(fēng)電接入真實(shí)數(shù)據(jù)為基礎(chǔ),驗(yàn)證了時(shí)間相關(guān)性的存在,并以Pearson相關(guān)系數(shù)(PCC)為相關(guān)性測度,提出了不確定變量的相關(guān)性約束。接著,根據(jù)利用不確定集合的離散型特點(diǎn),將相關(guān)性約束近似簡化為線性約束。本文將簡化后的相關(guān)性約束加入RUC中,并在列與限制生成(C&CG)算法的基礎(chǔ)上改進(jìn)Bender’s分解后子問題的求解方法,用于求解含有相關(guān)性約束的RUC。最后在真實(shí)風(fēng)電預(yù)測數(shù)據(jù)的基礎(chǔ)上,大量仿真實(shí)驗(yàn)驗(yàn)證了所提相關(guān)性約束的有效性。
在進(jìn)行風(fēng)電預(yù)測誤差的時(shí)間相關(guān)性建模之前,首先需要確認(rèn)其存在性。對(duì)于某一個(gè)時(shí)間序列x(t),自相關(guān)系數(shù)用于表征序列本身在不同時(shí)間的相關(guān)程序,該數(shù)值越高,表示歷史數(shù)據(jù)對(duì)當(dāng)前時(shí)刻的影響程度越大。自相關(guān)系數(shù)計(jì)算公式為:
(1)
t時(shí)刻的相對(duì)預(yù)測誤差et可以表示為:
(2)
記相對(duì)誤差序列為:
s0=[e1e2…eN-1]
(3)
s1=[e2e3…eN]
(4)
式中:N為風(fēng)電預(yù)測的樣本數(shù)量。顯然,s1是s0向后平移一個(gè)時(shí)間間隔的序列。
以德國巴登—符騰堡州的風(fēng)電接入數(shù)據(jù)為預(yù)測誤差序列樣本[22]。圖1是以2015年1月樣本數(shù)據(jù)序列生成s0和s1向量所繪制的散點(diǎn)圖。圖中:下標(biāo)i=1,2,…,N-1;j=2,3,…,N。從圖1中可以直觀地看到,s0和s1向量具有明顯的正相關(guān)性。為了進(jìn)一步確認(rèn)相關(guān)性水平,計(jì)算s0和s1的PCC。相關(guān)系數(shù)越大,則說明歷史數(shù)據(jù)對(duì)當(dāng)前數(shù)據(jù)的影響越強(qiáng),預(yù)測誤差的時(shí)間相關(guān)性越高。若相關(guān)系數(shù)大于0.6,則認(rèn)為預(yù)測誤差具有顯著的時(shí)間相關(guān)性??紤]到該地區(qū)氣候的季節(jié)特性,對(duì)風(fēng)電預(yù)測數(shù)據(jù)以月為單位進(jìn)行檢驗(yàn)。預(yù)測數(shù)據(jù)的時(shí)間間隔為1 h。
圖1 預(yù)測誤差序列散點(diǎn)圖Fig.1 Scatter plot of predictive error sequence
圖2是對(duì)2015年每個(gè)月內(nèi)的預(yù)測數(shù)據(jù)進(jìn)行相關(guān)性計(jì)算的結(jié)果??梢钥吹?各個(gè)月份都體現(xiàn)出了正相關(guān)性,平均相關(guān)系數(shù)在0.69左右,表明之前的預(yù)測誤差具有時(shí)間相關(guān)性假設(shè)是成立的。
圖2 以1 h為間隔的預(yù)測誤差樣本Pearson相關(guān)性系數(shù)Fig.2 Pearson correlation coefficient of 1 hour interval predictive error sequence
在確認(rèn)時(shí)間相關(guān)性后,可以利用這一點(diǎn),進(jìn)一步縮小不確定集合規(guī)模,將相關(guān)性較低的預(yù)測誤差序列剔除,降低魯棒優(yōu)化的保守程度。
圖3 某月內(nèi)每天預(yù)測誤差相關(guān)系數(shù)分布Fig.3 Distribution of daily correlation coefficient for predictive error within a month
圖3是某典型月內(nèi),每天的預(yù)測誤差相關(guān)系數(shù)分布直方圖。其中的紅色曲線是對(duì)預(yù)測誤差分布進(jìn)行Kernel擬合的曲線??梢钥吹?由于樣本數(shù)量較少,預(yù)測誤差分布規(guī)律不明顯,并非典型的正態(tài)分布。但是可以確定的是,置信水平為95%時(shí),樣本的相關(guān)系數(shù)均在[0.36,1]范圍內(nèi)。基于以上分析,可以對(duì)不確定集合加入時(shí)間相關(guān)性約束。
(5)
式中:pe為風(fēng)電出力的日前預(yù)測場景;p為風(fēng)電出力的可能場景;C(·)為PCC計(jì)算函數(shù);γ為可能場景與預(yù)測場景的時(shí)間相關(guān)性下限,與預(yù)設(shè)的置信水平有關(guān);ρ為置信水平。
式(5)本質(zhì)上是預(yù)測誤差時(shí)間序列自相關(guān)性的體現(xiàn)。由于相關(guān)性約束的存在,獨(dú)立分布的時(shí)間序列(比如白噪聲)將不被優(yōu)化模型考慮,因此魯棒優(yōu)化可以得到更加經(jīng)濟(jì)的結(jié)果。γ可以作為魯棒優(yōu)化保守程度的調(diào)節(jié)參數(shù),置信水平越高,γ越小,不確定集合越大,優(yōu)化結(jié)果越保守。
式(5)是從時(shí)間相關(guān)性角度建立的約束條件。但是相關(guān)系數(shù)計(jì)算函數(shù)C(·)計(jì)算過程復(fù)雜,是一個(gè)高度非線性的函數(shù)。式(5)難以直接加入魯棒優(yōu)化中,必須進(jìn)行簡化處理。
盒式不確定集合一般形式可以表示為:
(6)
Δ不等于0時(shí),式(6)描述的不確定集合是一個(gè)多維空間中的多面體。文獻(xiàn)[4]指出,魯棒優(yōu)化過程中,最差情景時(shí)不確定變量d取值必定為該多面體的某個(gè)極點(diǎn)。進(jìn)一步地,文獻(xiàn)[23]指出,當(dāng)Δ為整數(shù)時(shí),對(duì)于魯棒優(yōu)化問題,式(6)描述的集合可以改寫為以下形式:
(7)
容易理解,在式(7)表示方式下,原本連續(xù)的預(yù)測誤差序列e=d+-d-。所以式(5)中對(duì)從e中截取而生成的s0和s1的相關(guān)系數(shù)的約束,可以轉(zhuǎn)化成對(duì)d+和d-的約束??紤]到d+和d-的在式(7)完全對(duì)稱,下文的討論將以d+為例。
(8)
(9)
記d+序列之和為Δ+,這里定義的Δ+可以類比式(6)中的不確定預(yù)算Δ。并定義變化標(biāo)志變量v+為:
(10)
相應(yīng)的,記v+求和為Λ+,類比Δ+,稱Λ+為該序列的 “變化量預(yù)算”。因?yàn)棣?代表著原序列出現(xiàn)前后不一致的次數(shù),并給出以下定理。
證明1PCC的計(jì)算公式如下:
(11)
式中:cov(·)為求協(xié)方差函數(shù)。
(12)
(13)
這4種組合中,出現(xiàn){0,1}和{1,0}次數(shù)之和為Λ+。出現(xiàn){1,1}的次數(shù)為(2Δ+-Λ+)/2。出現(xiàn){0,0}的次數(shù)為N-1-(2Δ+-Λ+)/2-Λ+。
定理2是式(5)化簡的重要依據(jù)。預(yù)測誤差在使用連續(xù)集合式(6)描述時(shí),本身具有無窮多種可能場景。即使使用離散集合式(7)進(jìn)行描述,也存在多達(dá)2N種可能。由定理2可知,并不需要對(duì)如此多種可能場景的相關(guān)系數(shù)一一計(jì)算。所有Λ+和Δ+相同序列的相關(guān)系數(shù)都相同。為了更形象地解釋定理2,利用表1進(jìn)行說明。
表1 典型序列的PCC與不確定預(yù)算、變化量預(yù)算的關(guān)系Table 1 Relationship between PCC and uncertainty budget/variation budget in typical sequences
表1展示了6組典型序列的相關(guān)系數(shù)與不確定預(yù)算、變化量預(yù)算的關(guān)系。比較前三組和后三組序列可以看到,在變化量預(yù)算Λ+一定的情況下,即使序列具體排列不同,相關(guān)系數(shù)也不會(huì)變化。這樣一來,僅需要研究Λ+,Δ+和相關(guān)系數(shù)的關(guān)系即可。對(duì)相系數(shù)的約束可以轉(zhuǎn)化為對(duì)Λ+和Δ+參數(shù)的線性約束。
至此,考慮風(fēng)電相關(guān)性的不確定集合可表示為:
(14)
式中:Λ為總的不確定預(yù)算,即d+和d-序列的不確定預(yù)算Λ+與Λ-之和;上標(biāo)u表示啟動(dòng),d表示關(guān)停。
對(duì)于某個(gè)風(fēng)電場在N個(gè)時(shí)間間隔內(nèi)出力的不確定集合,式(7)需要2N個(gè)整數(shù)變量,而式(14)需要6N個(gè)整數(shù)變量,故而考慮時(shí)間相關(guān)性約束后,魯棒優(yōu)化的求解時(shí)間可能增加。
在第1節(jié)中,通過對(duì)相關(guān)系數(shù)的簡化,將非線性的相關(guān)性約束式(5)簡化為線性約束,并用式(14)表示改進(jìn)后的不確定集合。本節(jié)將討論含有式(14)的魯棒優(yōu)化求解算法。
在電力系統(tǒng)中,RUC的一般形式為(考慮到篇幅限制,這里只列出簡潔模型,詳細(xì)模型見附錄A):
(15)
s.t.Fx≤f
(16)
Ax+By≤g
(17)
Iyy+Idd=k
(18)
My+Nd≤l
(19)
式中:x為機(jī)組狀態(tài)變量,是二進(jìn)制離散變量;y為機(jī)組參數(shù)變量,是連續(xù)型變量;其余變量含義見文獻(xiàn)[4]。式(16)表示機(jī)組對(duì)機(jī)組狀態(tài)的約束,包括啟停邏輯約束、最小啟停時(shí)間約束、備用率約束等;式(17)表示機(jī)組狀態(tài)和機(jī)組出力關(guān)系約束,包括容量約束,爬坡率約束等;式(18)表示系統(tǒng)的功率平衡約束;式(19)表示傳輸線安全約束。
本文文獻(xiàn)[13]提出的C&CG算法進(jìn)行求解。由于式(14)為離散型集合,本文中對(duì)Bender’s分解后產(chǎn)生的子問題的求解不再采用文獻(xiàn)[13]的方式。Bender’s分解后子問題為:
(20)
寫出min問題部分對(duì)偶形式,式(20)變?yōu)殡p線性規(guī)劃問題:
(21)
s.t.b+λTB+ηTIy+θTM=0
(22)
(23)
式中:λ,η,θ分別為式(17)至式(19)的拉格朗日乘子。
顯然,對(duì)雙線性項(xiàng)(ηTIdd,θTNd)的處理是求解的關(guān)鍵。將式(14)同樣寫為簡潔形式:
{d|d=Tz+T0,Rz≤r}
(24)
式中:z為包括v,d+,d-的二進(jìn)制變量。則雙線性項(xiàng)變?yōu)棣荰Cz,θTGz,其中C=IdT,G=NT。以ηTCd為例:
(25)
顯然,只有cij不為0時(shí),對(duì)式(25)的計(jì)算才是有意義的??紤]到zj的二進(jìn)制性質(zhì),引入附加變量qk:
-zjQ≤qk≤zjQ
(26)
ηi-(1-zj)Q≤qk≤ηi+(1-zj)Q
(27)
其中,q的下標(biāo)k與i,j的關(guān)系為:矩陣C中第k個(gè)不為零的元素為cij。Q為q的上界。容易理解,式(26)和式(27)等價(jià)于當(dāng)zj等于0時(shí),qk等于0;當(dāng)zj不等于0時(shí),qk等于ηi。因此,qk=ηizj。
至此,子問題由雙線性規(guī)劃轉(zhuǎn)化為混合整數(shù)線性規(guī)劃,可以被CPLEX等成熟的商業(yè)求解器求解。但是,注意到該問題的求解效率與矩陣C和矩陣G的稀疏性密切相關(guān)。矩陣C由功率平衡約束轉(zhuǎn)化而來,規(guī)模較小;而矩陣G由傳輸線安全約束轉(zhuǎn)化而來,規(guī)模較大。在實(shí)際計(jì)算時(shí),可以先忽略傳輸線安全約束,待計(jì)算完成后,再校驗(yàn)各條線路上功率是否越限。如越限,則將該線路加入傳輸線安全約束中,再次計(jì)算。
為了驗(yàn)證本文提出的考慮時(shí)間相關(guān)性的不確定集合的有效性,對(duì)以下三種機(jī)組組合模型進(jìn)行對(duì)比。
1)傳統(tǒng)機(jī)組組合(TUC)。采用確定性優(yōu)化算法,不考慮風(fēng)電預(yù)測誤差。
2)RUC。采用第3節(jié)所述的魯棒優(yōu)化算法,考慮風(fēng)電預(yù)測誤差,但不考慮時(shí)間相關(guān)性,對(duì)不確定集的描述如式(7)所示。
3)計(jì)及時(shí)間相關(guān)性的魯棒機(jī)組組合(TRUC)。采用第3節(jié)所述的魯棒優(yōu)化算法,對(duì)不確定集的描述如式(14)所示。
電力系統(tǒng)模型采用IEEE 118節(jié)點(diǎn)標(biāo)準(zhǔn)算例,包括118條母線、54臺(tái)發(fā)電機(jī)以及186條傳輸線[24]。風(fēng)電接入和系統(tǒng)負(fù)荷數(shù)據(jù)來自德國巴登—符騰堡州電網(wǎng)運(yùn)營商Transnet BW的公開數(shù)據(jù)[22]。對(duì)風(fēng)電及負(fù)荷數(shù)據(jù)稍作調(diào)整,使其符合IEEE 118節(jié)點(diǎn)電力系統(tǒng)正常運(yùn)行的功率范圍。調(diào)整后系統(tǒng)最大負(fù)荷大約為6 500 MW,最大風(fēng)電滲透率大約為40%,平均風(fēng)電滲透率大約為8%。以該地區(qū)2015年運(yùn)行數(shù)據(jù)為歷史數(shù)據(jù),對(duì)不確定集合參數(shù)進(jìn)行估計(jì)。
三種機(jī)組組合算法在MATLAB 2014b及YALMIP上完成[25]。混合整數(shù)規(guī)劃采用CPLEX 12.6進(jìn)行求解。魯棒優(yōu)化結(jié)果的間隙誤差(即上下界之差除以下界)上限為0.1%。所有計(jì)算在配備處理器型號(hào)為Intel Core i7-4790 4核3.6 GHz以及16 GB RAM的個(gè)人計(jì)算機(jī)上完成。
為了驗(yàn)證本文提出的TRUC的有效性,首先以某典型日風(fēng)電預(yù)測情況為輸入數(shù)據(jù),計(jì)算三種機(jī)組組合模型的求解結(jié)果,并對(duì)機(jī)組組合結(jié)果進(jìn)行安全校驗(yàn)。安全校驗(yàn)指的是在風(fēng)電場真實(shí)出力情況下,日前制定的機(jī)組組合計(jì)劃是否能夠滿足電力系統(tǒng)安全性要求[5-7]。其數(shù)學(xué)模型為:
(28)
式中:s1至s4為保證式(17)和式(19)滿足而引入的松弛變量。
若R(x)=0則說明在s1至s4均為零的情況下,系統(tǒng)功率平衡約束、傳輸線安全約束等可以滿足。反之,R(x)>0則代表機(jī)組組合無法滿足系統(tǒng)安全性要求,此時(shí)需要實(shí)施切負(fù)荷或棄風(fēng)等緊急措施。在Δ等于8,Λ等于2時(shí),三種模型求解結(jié)果如表2及圖4所示。
表2 三種模型典型日優(yōu)化結(jié)果對(duì)比(Δ=8,Λ=2)Table 2 Comparison of optimization results of three models in a typical day (Δ=8, Λ=2)
圖4 典型日風(fēng)電預(yù)測及“最壞情況”求解結(jié)果曲線圖(Δ=8,Λ=2)Fig.4 Wind power prediction and calculation results of “the worst case scenario” in a typical day (Δ=8, Λ=2)
該算例的優(yōu)化結(jié)果顯示,相比于RUC算法,TRUC算法在保證安全性的基礎(chǔ)上,提高了機(jī)組組合經(jīng)濟(jì)性,降低了魯棒優(yōu)化的保守程度。為進(jìn)一步解釋TRUC機(jī)理,將RUC與TRUC計(jì)算出的“最壞情境”在圖4中標(biāo)出。
RUC算法得到的“最壞情境”如圖4中紅線標(biāo)注,其中出現(xiàn)多次風(fēng)電爬坡事件,風(fēng)電出力在預(yù)測值及最小值之間多次跳動(dòng)。根據(jù)前文分析,該情境出現(xiàn)概率較低。而TRUC算法中,在相同的不確定預(yù)算時(shí),由于變化量預(yù)算的存在,限制了預(yù)測誤差變化次數(shù),剔除了概率較低的情境,提高了優(yōu)化經(jīng)濟(jì)性。
TRUC優(yōu)化結(jié)果受到兩個(gè)參數(shù)的影響,不確定預(yù)算Δ和變化量預(yù)算Λ。下面比較不同的Δ和Λ取值對(duì)優(yōu)化結(jié)果的影響,如圖5所示。
圖5 不同Δ和Λ情況下的求解結(jié)果Fig.5 Optimization results under different Δ and Λ
從圖5中可以看到,對(duì)優(yōu)化結(jié)果影響最大的是不確定預(yù)算的值。比較Δ等于6,8,10的三條曲線,優(yōu)化結(jié)果大約變化了2%。而變化量預(yù)算對(duì)優(yōu)化結(jié)果的影響要相對(duì)小一些。以Δ曲線為例,比較Λ等于1和Λ等于6,成本變化了大約0.4%。并且,在Λ大于4時(shí),成本變化并不明顯。由此,可以將Λ看作一個(gè)輔助參數(shù),是在決策者對(duì)Δ進(jìn)行估計(jì)后,對(duì)調(diào)度保守性的更精細(xì)的調(diào)節(jié)手段。
3.2節(jié)中,已經(jīng)看到變化量預(yù)算對(duì)優(yōu)化結(jié)果的調(diào)節(jié)作用。但是擔(dān)心加入變化量預(yù)算的約束,是否會(huì)對(duì)系統(tǒng)的安全性造成影響。從圖5中可以看到,這種可能性是存在的。因此,對(duì)一段較長時(shí)間內(nèi)TRUC算法的安全性進(jìn)行評(píng)估。以德國巴登—符騰堡州2016年1月到7月間風(fēng)電接入為輸入數(shù)據(jù),計(jì)算三種算法的日前機(jī)組組合結(jié)果,并進(jìn)行安全評(píng)估。
Δ的取值一般由決策者估計(jì),Δ對(duì)優(yōu)化結(jié)果的影響已經(jīng)在多篇文章中有討論[4-7]。這里固定Δ取值為8,改變?chǔ)娜≈?來考察Λ對(duì)優(yōu)化結(jié)果及安全性的影響。
表3是三種機(jī)組組合模型計(jì)算結(jié)果的對(duì)比。其中,通過安全校驗(yàn)比例代表了不出現(xiàn)切負(fù)荷或棄風(fēng)等措施的可能性,該數(shù)值越大,代表算法的可靠性越高。從表3中可以看到,相比于RUC算法,TRUC算法提高了機(jī)組組合的經(jīng)濟(jì)性。在Λ=2時(shí),TRUC可靠性與RUC接近,機(jī)組組合成本降低了48萬美元。在Λ=3時(shí),可靠性與RUC相同,但成本降低了約8萬美元。
表3 三種模型長期優(yōu)化結(jié)果對(duì)比Table 3 Comparison of optimization results of three models in a long term
本文研究了風(fēng)電預(yù)測誤差的時(shí)間相關(guān)性特點(diǎn),并提出了TRUC算法。該算法中在“不確定預(yù)算”外,加入了“變化量預(yù)算”參數(shù),作為魯棒優(yōu)化保守性的另一個(gè)控制參數(shù)。由本文的理論研究及仿真實(shí)驗(yàn)得出如下結(jié)論。
1)某些預(yù)測方法得到的風(fēng)電預(yù)測誤差存在強(qiáng)烈的時(shí)間相關(guān)性。利用時(shí)間相關(guān)特點(diǎn),可以進(jìn)一步縮小不確定集合范圍,降低魯棒優(yōu)化保守性。
2)利用不確定集合離散型特點(diǎn),預(yù)測誤差序列的自相關(guān)系數(shù)計(jì)算公式可以得到簡化。相關(guān)性約束可以近似簡化為線性形式。
3)在真實(shí)風(fēng)電數(shù)據(jù)進(jìn)行的仿真實(shí)驗(yàn)上,驗(yàn)證了該算法能夠在不降低機(jī)組組合魯棒性的前提下,提高調(diào)度計(jì)劃的經(jīng)濟(jì)性。
在含風(fēng)電的RUC問題中,如何將“棄風(fēng)”加入模型研究中,通過魯棒優(yōu)化手段,進(jìn)一步降低“棄風(fēng)”比例,是本文下一步研究的方向。
附錄見本刊網(wǎng)絡(luò)版(http://www.aeps-info.com/aeps/ch/index.aspx)。