預(yù)估–校正原對偶內(nèi)點(diǎn)具體實現(xiàn)步驟如下：

1）優(yōu)化變量初始化。迭代從第0次開始，設(shè)置k?0，定義碰壁參數(shù)初值?0，在滿足非負(fù)變量嚴(yán)格為正的前提下對優(yōu)化變量進(jìn)行初始化，給出初始運(yùn)行點(diǎn)y0。

式中：x為優(yōu)化變量向量；f(x)為優(yōu)化目標(biāo)函數(shù)；g(x) 和h(x)分別為等式及不等式約束函數(shù)向量；h和h 分別為不等式約束的下限和上限向量。

2）預(yù)估環(huán)節(jié)。在預(yù)估環(huán)節(jié)中，在式(37)中省去

52 中 國 電 機(jī) 工 程 學(xué) 報 第32卷

了碰壁參數(shù)?k和?項，可計算獲得當(dāng)前運(yùn)行點(diǎn)仿射方向。

45 MW。風(fēng)電場中風(fēng)速均值取為各個時段的風(fēng)速預(yù)測值。3個時段風(fēng)速預(yù)測值分別為8.更多內(nèi)容請訪問久久建筑網(wǎng)
5、10.5和

3）校正環(huán)節(jié)。在當(dāng)前運(yùn)行點(diǎn)計算完整的牛頓方向。

13 m/s，風(fēng)速預(yù)測誤差的方差取預(yù)測值的8%。系統(tǒng)風(fēng)險檻值?設(shè)置為0.05。

4）更新優(yōu)化變量。沿牛頓方向計算迭代步長，更新原對偶變量。

1）模式1。

在這種運(yùn)行模式下，假設(shè)風(fēng)電場歸電力系統(tǒng)所有，那么總成本中的風(fēng)電成本和負(fù)旋轉(zhuǎn)備用成本均不計。優(yōu)化結(jié)果為：總發(fā)電成本為16.157萬$，其中正備用成本為5 061.8 $，常規(guī)發(fā)電機(jī)組燃料成本為15.651萬$。表3為所得的調(diào)度分配結(jié)果，包括常規(guī)發(fā)電機(jī)組計劃出力。各時段風(fēng)電場計劃出力及系統(tǒng)正旋轉(zhuǎn)備用需求如表4所示。

首先分析設(shè)定不同風(fēng)險檻值?對動態(tài)調(diào)度結(jié)果的影響。正如圖1所示，動態(tài)經(jīng)濟(jì)調(diào)度的目標(biāo)函數(shù)值隨著?的增加而減小。說明系統(tǒng)對可靠性要求的

表3 系統(tǒng)常規(guī)發(fā)電機(jī)輸出的有功功率 Tab. 3 Scheduled power outputs of

conventional generators

機(jī)組號

時段1

1 2 3 4 5 6 7 8 9 10

453.461 2 443.462 4 323.465 7 275.156 0 227.306 8 159.994 0 129.996 1 100.100 0 23.070 6 55.000 0

機(jī)組各時段出力/MW

時段2 452.498 1 442.499 2 322.502 2 244.321 8 226.773 2 159.993 9 129.996 1 84.192 7 21.446 0 55.000 0

時段3 452.266 9 442.269 0 322.272 7 175.638 5 226.708 2 159.993 9 129.996 1 75.970 7 21.222 2 55.000 0

5）判斷收斂條件。如果新的運(yùn)行點(diǎn)滿足收斂條件，迭代停止；否則，置k?k?1，更新碰壁參數(shù)，返回步驟2）。

限于篇幅，關(guān)于迭代步長、初始化及判斷收斂條件部分不再闡述，具體過程參見本文作者已發(fā)表的文獻(xiàn)[23]。

4 算例分析

將本文所提的基于風(fēng)險備用約束的含風(fēng)電場動態(tài)經(jīng)濟(jì)調(diào)度優(yōu)化模型在具有一個并網(wǎng)風(fēng)場的10機(jī)組系統(tǒng)中進(jìn)行仿真計算。該風(fēng)場共有100臺額定功率為2 MW的風(fēng)力發(fā)電機(jī)。假設(shè)風(fēng)力發(fā)電機(jī)組不提供旋轉(zhuǎn)備用且不考慮其強(qiáng)迫停運(yùn)的可能性。常規(guī)發(fā)電機(jī)組強(qiáng)迫停運(yùn)的概率如表1所示。風(fēng)場中使用如下風(fēng)速參數(shù)：切入風(fēng)速vi?4 m/s、額定風(fēng)速vr??

12.5 m/s以及切出風(fēng)速vo?20 m/s。各個常規(guī)發(fā)電機(jī)組的燃料特性以及有功出力極限見文獻(xiàn)[24]，常規(guī)發(fā)電機(jī)組爬坡率數(shù)據(jù)如表2所示。系統(tǒng)正旋轉(zhuǎn)備用成本系數(shù)為20 $/(MW?h)。為方便地對優(yōu)化調(diào)度結(jié)果及其影響因素進(jìn)行分析，系統(tǒng)研究周期設(shè)定為 3

h，每個時段為1 h。假設(shè)各個時段負(fù)荷不變，負(fù)荷預(yù)測值和預(yù)測誤差的標(biāo)準(zhǔn)差分別為2 250和

表1 常規(guī)發(fā)電機(jī)組強(qiáng)迫停運(yùn)的概率

Tab. 1 Probability of convetional generator trip

發(fā)電機(jī) 1 2 3 4 5

強(qiáng)迫停運(yùn)概率 0.000 8 0.002 0 0.000 3 0.002 0 0.002 0

發(fā)電機(jī) 6 7 8 9 10

強(qiáng)迫停運(yùn)概率 0.000 5 0.000 8 0.002 0 0.000 3 0.002 0

表4 風(fēng)電場計劃出力及系統(tǒng)正旋轉(zhuǎn)備用需求 Tab. 4 Scheduled power outputs of wind farm and

up spinning reserve demand

時段號

風(fēng)電場出力/MW

正旋轉(zhuǎn)備用需求/MW

1 58.991 9 82.495

8 2

110.781 5

85.066 1

3 188.666 7 85.530 5

總發(fā)電成本/105 $

表2 常規(guī)發(fā)電機(jī)組的爬坡率

Tab. 2 Ramp rate data of conventional generator

發(fā)電機(jī)

ru/

rd/

發(fā)電機(jī)

ru/ (MW/min)

rd/ (MW/min)

(MW/min) (MW/min)

1 2.67 6 1.67 1.67 2 2.67 7 1.00 1.00 3 2.67 8 1.00 1.00 4 1.67 9 1.00 1.00 5 1.67 1.00

??

圖1 發(fā)電總成本關(guān)于風(fēng)險檻值?的函數(shù)曲線 Fig. 1 Total generation cost as a function of ?

第1期 周瑋等：計及風(fēng)險備用約束的含風(fēng)電場電力系統(tǒng)動態(tài)經(jīng)濟(jì)調(diào)度 53

降低有利于經(jīng)濟(jì)性的提升。

根據(jù)圖2－4所描述的曲線來分析不同正備用成本系數(shù)對三時段調(diào)度結(jié)果的影響。此時，保持風(fēng)發(fā)電總成本會隨著它的增加而不斷22 $/(MW?h)時，

增加。超過22 $/(MW?h)以后，發(fā)電總成本則幾乎不變。這一結(jié)果可以通過圖3和圖4來解釋。為了在滿足系統(tǒng)風(fēng)險要求的前提下防止不斷增長的備 用成本系數(shù)導(dǎo)致發(fā)電總成本的大幅度增加，系統(tǒng)必須限制風(fēng)電場的計劃出力，從而降低了正備用需求量。當(dāng)正備用成本系數(shù)增大到一定程度的時候，風(fēng)電場的計劃出力和系統(tǒng)需求的正旋轉(zhuǎn)備用容量不再降低，幾乎保持恒定的數(shù)值。

kr/($/(MW?h))?

2）模式2。

當(dāng)風(fēng)電場成為獨(dú)立的運(yùn)營機(jī)構(gòu)后，電網(wǎng)必須向風(fēng)電供應(yīng)商支付相應(yīng)的發(fā)電費(fèi)用。調(diào)度時需要考慮價格取為10 $/(MW?h)，風(fēng)險檻值?也取0.05，其他參數(shù)與模式1相同。表5和表6給出了負(fù)旋轉(zhuǎn)備用成本系數(shù)為11 $/(MW?h)時的調(diào)度結(jié)果。優(yōu)化結(jié)果表 明，總發(fā)電成本為16.795 6萬$，其中常規(guī)發(fā)電機(jī)組燃料成本為15.597 3萬$，風(fēng)電功率成本為 3 857.5 $，正備用成本為5 608.7 $，負(fù)備用成本為2 517.3 $。進(jìn)而研究負(fù)旋轉(zhuǎn)備用成本系數(shù)(即風(fēng)能浪費(fèi)懲罰程度)改變對調(diào)度結(jié)果的影響。

表5 系統(tǒng)常規(guī)發(fā)電機(jī)輸出的有功功率 Tab. 5 Scheduled power outputs of

conventional generator

機(jī)組號

時段1

機(jī)組各時段出力/MW

時段2

時段3

險檻值?為0.05。當(dāng)正備用成本系數(shù)小于 風(fēng)電成本以及負(fù)旋轉(zhuǎn)備用成本。此時，將風(fēng)電成本

總發(fā)電成本/105 $

1 449.37 447.20 447.31 2 439.37 437.21 437.28 3 319.37 317.20 317.34

4 276.63 251.18 182.33 5 226.37 225.40 225.38 6 159.98 159.91 159.61 7 129.99 129.98 129.86 8 102.48 84.21 75.40

圖2 發(fā)電總成本關(guān)于正備用成本系數(shù)的函數(shù)曲線 Fig. 2 Total generation cost as a function of

up reserve cost coefficient

旋轉(zhuǎn)備用需求/MW

kr/($/(MW?h))?

9 10

24.77 55.00

21.99 55.00

22.12 55.00

表6 風(fēng)電場計劃出力及系統(tǒng)正、負(fù)旋轉(zhuǎn)備用需求 Tab. 6 Scheduled power outputs of wind farm,

up and down spinning reserve demand

時段號

風(fēng)電場 出力/MW

正旋轉(zhuǎn)備用 需求/MW

負(fù)旋轉(zhuǎn)備用 需求/MW

1 66.67 90.19 71.51

圖3

風(fēng)電場計劃出力關(guān)于正備用成本系數(shù)函數(shù)曲線 Fig. 3 Scheduled wind power output as a function of

up reserve cost coefficient

旋轉(zhuǎn)備用需求/MW

2 120.72 95.01 86.77 3 198.36 95.24 70.56

從圖5可以看出，系統(tǒng)總發(fā)電成本隨著負(fù)備用成本系數(shù)的增大而逐漸增加。由于加大了對浪費(fèi)風(fēng)電資源的懲罰力度，因而會刺激電力系統(tǒng)調(diào)度運(yùn)行人員增加風(fēng)電計劃出力值，以降低在同一風(fēng)險水平下所需求的負(fù)備用容量值，從而避免調(diào)度過程中產(chǎn)生的負(fù)備用成本的進(jìn)一步提升；相反地，風(fēng)電功率增加必然會提高系統(tǒng)失負(fù)荷的風(fēng)險。因此，系統(tǒng)需

kr/($/(MW?h))?

要增加正旋轉(zhuǎn)備用以確保失負(fù)荷風(fēng)險保持在要求的范圍內(nèi)。圖6—8分別展示了各個時段風(fēng)電計劃出力、系統(tǒng)正旋轉(zhuǎn)備用需求和負(fù)旋轉(zhuǎn)備用需求隨著

圖4 正旋轉(zhuǎn)備用需求關(guān)于正備用成本系數(shù)函數(shù)曲線 Fig. 4 Spinning reserve demand as a function of

up reserve cost coefficient

kp的變化關(guān)系。

54

總發(fā)電成本/105 $

中 國 電 機(jī) 工 程 學(xué) 報 第32卷

型，在模型目標(biāo)函數(shù)中考慮了正、負(fù)旋轉(zhuǎn)備用成本。通過引入基于失負(fù)荷風(fēng)險和風(fēng)能浪費(fèi)風(fēng)險指標(biāo)的備用約束條件，確保了系統(tǒng)能夠在風(fēng)險允許的范圍內(nèi)安排各個常規(guī)發(fā)電機(jī)組的出力計劃，同時可獲取風(fēng)電場的計劃出力值和正、負(fù)旋轉(zhuǎn)備用需求量，最終實現(xiàn)經(jīng)濟(jì)性最優(yōu)。研究表明，含風(fēng)電場的電力系

kp/($/(MW?h))?

統(tǒng)動態(tài)經(jīng)濟(jì)調(diào)度結(jié)果與模型中風(fēng)險檻值、正/負(fù)備用成本系數(shù)等因素密切相關(guān)。這一模型的提出可以為電力系統(tǒng)調(diào)度人員根據(jù)系統(tǒng)要求合理地安排風(fēng)電和其他常規(guī)機(jī)組的出力計劃提供理論依據(jù)。

圖5 發(fā)電總成本關(guān)于負(fù)備用成本系數(shù)的函數(shù)曲線 Fig. 5 Total generation cost as a function of

down reserve cost coefficient

風(fēng)電計劃出力/MW

參考文獻(xiàn)

[1] Potter C W，Negnevitsky M．Very short-term wind

forecasting for tasmanian power generation[J]．IEEE Transactions on Power Systems，2006，21(2)：965- 972． [2] 楊秀媛，肖洋，陳樹勇．風(fēng)電場風(fēng)速和發(fā)電功率預(yù)測研

究[J]．中國電機(jī)工程學(xué)報，2005，25(11)：1-5．

kp/($/(MW?h))?

Yang Xiuyuan，Xiao Yang，Chen Shuyong．Wind speed and generated power forecasting in wind farm [J]．Proceedings of the CSEE，2005，25(11)：1-5(in Chinese)．

[3] 袁鐵江，晁勤，李義巖，等．大規(guī)模風(fēng)電并網(wǎng)電力系統(tǒng)

經(jīng)濟(jì)調(diào)度中風(fēng)電場出力的短期預(yù)測模型[J]．中國電機(jī)工

圖6 風(fēng)電計劃出力關(guān)于負(fù)備用成本系數(shù)函數(shù)曲線

Fig. 6 Scheduld wind power output as a function of down reserve cost coefficient

正旋轉(zhuǎn)備用需求/MW

程學(xué)報，2010，30(13)：23-27．

Yuan Tiejiang，Chao Qin，Li Yiyan，et al．Short-term wind power output forecasting model for economic dispatch of power system incorporating large-scale wind farm [J]．Proceedings of the CSEE，2010，30(13)：23-27(in

kp/($/(MW?h))?