聶永輝，馮浩然，于永利，王中杰

（1.東北電力大學 教務處，吉林 吉林 132012；2.東北電力大學 電氣工程學院，吉林 吉林 132012；3.吉林石化分公司精細化學品廠，吉林 吉林 132022）

0 引言

目前，以風電為代表的可再生能源得到了快速發(fā)展［1-2］，研究含風電并網(wǎng)的大規(guī)模電力系統(tǒng)優(yōu)化運行問題是一個基礎而又重要的課題［3］。各國學者主要采用2種方法應對風電本身具有的不確定性：一種方法是基于風電功率預測的大規(guī)模風電并網(wǎng)電力系統(tǒng)優(yōu)化求解，為減少風電不確定性，已有多種方法用來提高對未來時間段內風電機組功率出力預測的準確性［4-6］；另一種是在最優(yōu)模型中考慮到這種不確定性的影響，即采用期望值模型、機會約束規(guī)劃和相關機會規(guī)劃等不確定性的理論和方法來研究大規(guī)模風電接入下的優(yōu)化運行問題［7-9］。

風力發(fā)電不需要消耗石油、煤炭等非可再生能源，可以節(jié)省燃料成本，具有保護環(huán)境價值。如果優(yōu)化潮流模型沒有考慮傳統(tǒng)火電機組環(huán)境成本，將造成部分電源的環(huán)境價值被無償占有，這不僅阻礙了可再生能源的發(fā)展，不能合理反映以風電為代表的新能源環(huán)境價值，而且使得污染比較嚴重的傳統(tǒng)能源過度發(fā)展，加重了環(huán)境污染程度。另一方面，由于風能具有波動性、間歇性和隨機性等特點，為保證電網(wǎng)頻率的變化在允許的范圍內，調度部門常以此為借口決策風電場棄風運行，降低了風電的利用效率，造成風能資源浪費［2］。

為使風電環(huán)境價值得到合理的反映，國內外研究人員對含風電的優(yōu)化計算問題進行了大量的研究。文獻［3］考慮了風電功率的波動性和難以預測性對動態(tài)經(jīng)濟調度旋轉備用的新要求，引入了正、負旋轉備用約束，將火電機組閥點效應引入到發(fā)電總成本中，使得優(yōu)化模型更趨合理。文獻［10］從發(fā)電收益、發(fā)電成本、輔助服務成本、環(huán)保效益等各方面因素考慮風力并網(wǎng)發(fā)電的價值。文獻［11］引入和可利用的風能與實際發(fā)出的風能之差成正比的懲罰成本，然后轉化成和風能概率密度相關的積分形式來表示風能的隨機性，仿真表明引入懲罰成本確保了風能得到充分利用。文獻［12］采用隨機模擬技術分析風電并網(wǎng)對電力系統(tǒng)運行方式和發(fā)電成本的影響規(guī)律。文獻［13］通過計算下一調度日風電場實際出力的條件期望與計劃出力的差值，確定了風電對系統(tǒng)正、負旋轉備用的需求。

可以看出，上述文獻并未考慮棄風損失。本文風電場采用籠型異步風力發(fā)電機組，在風電場功率輸出預估基礎上，引入基于極限穿透功率的風電場棄風運行懲罰成本（WAOPC）來度量由于風電場棄風行為而造成的功率損失；引入風電備用成本（RCCC）來處理隨著系統(tǒng)風電穿透水平增加而引起的系統(tǒng)備用上升；同時考慮火電名義環(huán)境補償成本（NECC）來度量常規(guī)能源在發(fā)電時污染環(huán)境所帶來的成本。因此，在電力市場環(huán)境下，考慮風電特性、能源利用效率、整個系統(tǒng)運行及環(huán)境保護等因素，建立大規(guī)模風電并網(wǎng)動態(tài)優(yōu)化潮流新模型。該模型考慮電能生產(chǎn)中的環(huán)境污染和資源消耗等外部成本，引入了棄風成本，使以風電為代表的新能源價值得到合理的體現(xiàn)。

1 異步風電機組建模

籠型異步發(fā)電機在超同步速運行情況下以發(fā)電方式運行，此時吸收風力機提供的旋轉機械能，輸出有功功率，但其自身沒有勵磁裝置，因此沒有電壓調節(jié)功能，需從電網(wǎng)吸收無功功率提供其建立磁場所需的勵磁電流，為了減小風電場并網(wǎng)運行對系統(tǒng)電壓分布的影響，一般在籠型異步風電機組機端安裝并聯(lián)電容器組以補償無功損耗，其簡化等值電路如圖1所示。圖中，忽略定子電阻，x1為定子漏抗，xm為激磁電抗，x2為轉子漏抗，s為轉差率，r2為轉子電阻。當圖中功率參考方向采用發(fā)電方式時，可以推出下面的功率表達式：

其中，x=x1+x2。

由上述兩式可推出無功功率為：

圖1 異步風力發(fā)電機的簡化等值電路Fig.1 Simplified equivalent circuit of asynchronous wind-power induction generator

風力發(fā)電機輸出功率和風速的函數(shù)關系采用分段函數(shù)表示為［7］：

其中，vci、vr、vco分別為風機切入風速、額定風速和切出風速；PW為風機有功功率輸出；Pr為風機額定功率輸出；f（v）一般采用一次函數(shù)表達式來近似表示。

在某一時段，風機輸出功率可用其期望值PWF來近似，即將典型日風速分布曲線轉化成典型日功率分布曲線，利用下式求解得出：

其中，t0和t1分別為某時段的起始時刻和終止時刻。

2 大規(guī)模異步風電并網(wǎng)動態(tài)優(yōu)化潮流建模

2.1 基于風電場出力預估的風電備用容量補償成本

風電并網(wǎng)所增加的備用容量是由于風能的隨機性引起的，從而導致風電場預報出力與實際出力存在偏差，因此風電并網(wǎng)所增加的備用容量費用可表示為：

其中，CRCCCj為風電備用容量補償成本系數(shù)，其值與提供備用容量的火電機組的排污特性有關；PWFj為第j個風電場的輸出功率的期望值；PWj為第j個風電場計劃調度出力；nW為風電場個數(shù)。

2.2 基于極限穿透功率的風電場棄風運行懲罰成本

風電極限穿透功率是指風電場裝機容量占系統(tǒng)總負荷的比例。其是從整個電力系統(tǒng)出發(fā)，表明一個給定電網(wǎng)可以承受的最大風電容量大小，旨在考慮在風電投入和退出運行的2種情況下，電力系統(tǒng)的可調節(jié)容量應能保證電網(wǎng)頻率的變化在允許的范圍內［2］。這就使得調度部門常以此為借口決策風電場棄風運行，降低了風電的利用效率，造成風能資源浪費。定義由于調度原因造成的風電場棄風行為為風電場棄風運行，在定義的基礎上，提出了基于風電場極限穿透功率的風電場棄風運行懲罰成本，即：

其中，F(xiàn)WAOPC為總的風電場棄風運行懲罰成本；PWPPWj為第j個風電場的裝機容量，雖然風電極限穿越功率為一個全網(wǎng)的概念，但其是指所有風電場裝機容量之和占系統(tǒng)總負荷的比例，本文中將定義為第j個風電場的極限穿越功率；CWAOPCj為第j個風電場棄風運行懲罰價格系數(shù)。

2.3 火電名義環(huán)境補償成本

常規(guī)能源在生產(chǎn)電能時會釋放各種污染氣體及顆粒，對環(huán)境造成了巨大的危害。然而電力系統(tǒng)常規(guī)最優(yōu)潮流模型沒有體現(xiàn)火電的環(huán)境成本，從而制約了潔凈能源的發(fā)展。在體現(xiàn)火電購電成本的傳統(tǒng)優(yōu)化模型基礎上，根據(jù)火力發(fā)電機組的排污特性，增加火電機組的環(huán)境補償成本，這個附加的成本在機組的調度優(yōu)先權中體現(xiàn)，并不附加在消費者的購電費用中，因此將該成本定義為火電機組的名義環(huán)境補償成本。

火電機組的名義環(huán)境補償成本表示為：

其中，F(xiàn)NECC為總的火電機組的名義環(huán)境補償成本；nG為系統(tǒng)中的火電機組數(shù)目；CNECCi為火電名義環(huán)境補償成本價格系數(shù)；fei為排污特性方程，其是通過將單位時間內火電機組排放的各種污染物折算成NOx的重量來表示，fei的表達式為：

其中，PGi為第i臺火電機組的出力；αi、βi、γi、ζi、λi為第i臺火電機組的排污特性系數(shù)，可通過測量-擬合得到。

2.4 大規(guī)模異步風電并網(wǎng)動態(tài)優(yōu)化潮流新模型

2.4.1 目標函數(shù)

為體現(xiàn)傳統(tǒng)火電機組環(huán)境成本和風電機組的環(huán)境價值，本文采用火電名義環(huán)境補償成本來度量常規(guī)能源在發(fā)電時污染環(huán)境所帶來的成本；在風電場功率輸出預估基礎上，引入基于極限穿透功率的風電場棄風運行懲罰成本來度量由于風電場棄風行為而造成的功率損失；同時，引入風電備用成本來處理隨著系統(tǒng)風電穿透水平增加而引起的系統(tǒng)備用上升。因此，在電力市場環(huán)境下，考慮風電特性、能源利用效率、整個系統(tǒng)運行及環(huán)境保護等因素，以火力發(fā)電機有功出力和無功出力、異步風電機組的有功出力、并聯(lián)無功補償裝置的無功輸出、節(jié)點電壓實部和虛部以及有載調壓變壓器變比為優(yōu)化變量，建立了大規(guī)模風電并網(wǎng)動態(tài)優(yōu)化潮流新模型。該模型的目標函數(shù)可表示為：

其中，t代表調度時段；F（t）為系統(tǒng)總的成本費用，由風電備用容量補償成本、基于極限穿透功率的風電場棄風運行懲罰成本、風電場營運成本、火電機組的名義環(huán)境補償成本、火電營運成本5個部分組成。

風電機組的營運成本包括風電的投資成本和運行成本。在整個壽命周期內，可表示為與風電發(fā)電出力呈線性關系［2］：

其中，CWOPj為第j個風電場營運成本系數(shù)。

火電營運成本考慮火電投資、能量消耗和運行維護等成本，可表示為［2］：

其中，ai、bi、ci為與火力發(fā)電機組的耗量特性有關的火電機組的購電成本特性系數(shù)，可根據(jù)火電機組的報價曲線來確定。

2.4.2 等式約束

基于風電場期望值和異步發(fā)電機電壓與無功方程（3），動態(tài)優(yōu)化潮流必須滿足每個調度時段的第i個交流潮流平衡方程，可表示為：

其中，為籠型異步風電機組的有功計劃出力；為由式（3）表示的籠型異步風電機組的無功需求，是關于籠型異步風電機組機端電壓和有功出力的函數(shù)；分別為火力發(fā)電機有功和無功出力；為并聯(lián)補償裝置無功出力，包括風電場的無功補償出力；分別為節(jié)點電壓的實部和虛部；為有載調壓變壓器變比。對于通過風電匯集升壓站接入公共電網(wǎng)的風電場，配置容性無功補償容量能補償風電場滿發(fā)時送出線路上的無功損耗；配置感性無功補償容量能補償風電場空載時送出線路上的充電功率。風電場無功補償裝置容量總和一般不小于風電裝機容量的30%～50%。

2.4.3 不等式約束

不等式約束包括節(jié)點電壓幅值、火電發(fā)電機有功出力和無功輸出、風電機組有功輸出、并聯(lián)無功補償輸出（包括風電接入點的并聯(lián)無功補償出力）、有載調壓變壓器變比和支路負載的上下限約束限制及火電機組的爬坡速率約束。分別表示為：

其中，分別為第i臺火電機組的上升和下降速率；ΔT（t）為調度時段的時間間隔；n、nG、nW、nC、nT、nB分別為節(jié)點、常規(guī)發(fā)電機、風電場、并聯(lián)無功補償裝置和支路數(shù)量；rij、xij分別為支路電阻、電抗；為支路電流的上限。

在不等式約束中，只包含了籠型異步風力發(fā)電機組有功出力的上下限約束，而沒有包含其無功需求的上下限約束，這是因為其無功需求通過式（1）、（2）可表示為端電壓和轉速的函數(shù)，當轉速固定時，無功需求隨著端電壓的升高而升高。因此，當風電機組穩(wěn)定運行時，可利用端電壓的約束限制無功需求的范圍。

本文所論述的風火電運行成本建模中，各成本系數(shù)的取值反映了各個成本因素在總成本費用中的權重關系。該模型既考慮了傳統(tǒng)電力生產(chǎn)所帶來的環(huán)境污染和資源消耗等外部成本，同時也計及了因棄風行為造成的風電損失，使以風電為代表的新能源形式的電能價值得到合理的反映。本文利用預測校正原對偶內點法實現(xiàn)上述優(yōu)化問題，關于該算法的介紹詳見文獻［14-15］，這里不再贅述。

3 算例仿真

本文對IEEE30節(jié)點標準測試系統(tǒng)進行修改，將風電場通過升壓變壓器和1號節(jié)點相連，風電場容量為20×2 MW，風電場無功補償裝置容量為12 Mvar，風力機的切入風速、額定風速和切出風速分別為3、14、20 m/s，異步風力發(fā)電機的定子阻抗為j0.09985 p.u.，轉子阻抗為 0.00373+j0.10906 p.u.，勵磁電抗為3.54708p.u.，調度時段間隔為30 min，設火電機組上升和下降速率相等，對于2、22、27、23、13號節(jié)點所連接的火電機組的爬坡速率分別為0.30、0.24、0.18、0.24、0.30 MW/min。為方便起見，本文用文獻［7］中式（4）、（5）得出的風電場有功功率期望值計算結果，共包含5個時段，各時段的計算結果分別是 40、28、17、0、10 MW。首先對第2 個時段進行靜態(tài)優(yōu)化求解，以研究各個成本系數(shù)對機組調度權的影響，然后對所建動態(tài)優(yōu)化建模型進行求解。另外，本文各成本系數(shù)按比例采用相對價格進行仿真。

3.1 風電購電成本系數(shù)變化對機組出力的影響

設風電備用容量成本、棄風運行懲罰成本和火電名義環(huán)境補償成本系數(shù)為0。風電購電成本系數(shù)變化時，各機組出力（標幺值，后同）的仿真優(yōu)化結果如圖2（a）所示。由圖可知，在傳統(tǒng)的價格機制下，不考慮各機組的名義環(huán)境補償成本時，對環(huán)境破壞大的機組反而擁有更大的調度優(yōu)先權，火電機組2一直處于滿負荷狀態(tài)運行（該機組額定出力為0.8p.u.）；風電購電成本系數(shù)大于3.1時，除火電機組4外，其他所有火電機組都獲得不同容量的有功出力，而風電機組卻處于停運狀態(tài)，只有當購電成本系數(shù)降到足夠低的水平時才能獲得調度權。圖2（b）顯示了風電購電成本系數(shù)變化對優(yōu)化后總成本的影響，風電購電成本系數(shù)小于2.9左右時，優(yōu)化后總成本隨著風電購電成本系數(shù)變化緩慢地呈線性增加，但風電購電成本系數(shù)大于2.9左右時，優(yōu)化后總成本隨著風電購電成本系數(shù)變化大幅度地增加，最后穩(wěn)定在 673.8239$/h不變。

圖2 風電購電成本系數(shù)變化對機組出力的影響Fig.2 Influence of wind-power purchase cost coefficient on unit output

3.2 火電機組名義環(huán)境補償成本系數(shù)變化對機組出力的影響

設風電備用容量補償成本和棄風運行懲罰成本系數(shù)為0，風電運營成本系數(shù)為1?；痣姍C組名義環(huán)境補償成本系數(shù)變化時，各機組出力仿真優(yōu)化結果如圖3（a）所示。由圖可知，機組3、5和6的調度出力隨名義環(huán)境補償成本系數(shù)的提高而降低，對環(huán)境更友好的火電機組2、3一直處于滿發(fā)狀態(tài)，當名義環(huán)境補償成本系數(shù)為2.6時，風電機組由于火電機組名義環(huán)境補償成本的增加而獲得了調度權?？傊?，一方面，名義環(huán)境補償成本改變了火電機組的調度權；另一方面，名義環(huán)境補償成本賦予風電等環(huán)境友好的可再生能源調度優(yōu)先權，確保風電的優(yōu)先發(fā)展。圖3（b）顯示了名義環(huán)境補償成本系數(shù)對優(yōu)化后總成本的影響，可以看出優(yōu)化后總成本曲線是一條不光滑曲線，這是排污特性方程非線性項作用的結果，當火電機組名義環(huán)境補償成本系數(shù)達到2.6左右時，優(yōu)化后總成本大幅度下降。

圖3 名義環(huán)境補償成本系數(shù)變化對機組出力的影響Fig.3 Influence of NECC coefficient on unit output

3.3 風電備用容量補償成本系數(shù)變化對機組出力的影響

設火電機組名義環(huán)境補償成本系數(shù)為0，風電棄風運行懲罰成本系數(shù)為2，風電運營成本系數(shù)為2.2。風電備用容量補償成本系數(shù)變化時各機組優(yōu)化分配仿真結果如圖4（a）所示。由圖可知，當風電備用容量補償成本系數(shù)比較小時，對環(huán)境更友好的機組2、3和風電機組獲得優(yōu)先調度權而處于滿發(fā)狀態(tài)，由于沒有考慮火電機組名義環(huán)境補償成本，所有機組出力固定不變，但當風電備用容量補償成本系數(shù)達2.7左右時，由于風電備用容量補償成本的影響，風電出力下降，其他未滿發(fā)的火電機組出力上升。圖4（b）顯示了風電備用容量補償成本系數(shù)變化對優(yōu)化后總成本的影響，風電備用容量補償成本系數(shù)小于2.7左右時，優(yōu)化后總成本固定不變，但風電備用容量補償成本系數(shù)大于2.7左右時，優(yōu)化后總成本隨著風電購電成本系數(shù)變化大幅度地增加。

圖4 風電備用容量補償成本系數(shù)變化對機組出力的影響Fig.4 Influence of wind-power RCCC coefficient on unit output

3.4 風電棄風運行懲罰成本系數(shù)變化對機組出力的影響

火電名義環(huán)境補償成本、備用容量補償成本和風電運營成本系數(shù)分別為0、2和2.2時，風電棄風運行懲罰成本系數(shù)變化時各機組優(yōu)化仿真結果如圖5（a）所示。由于風電備用容量補償成本和風電運營成本的作用，風電機組在風電棄風運行懲罰成本系數(shù)比較小時未獲得調度權，當風電棄風運行懲罰成本系數(shù)達到1.1時，風電機組獲得調度權并處于滿發(fā)狀態(tài)，可防止調度部門以穿透功率為借口導致棄風行為；由于沒有考慮火電機組名義環(huán)境補償成本，對環(huán)境更友好的機組2、3獲得優(yōu)先調度權而處于滿發(fā)狀態(tài)，而其他能耗大的機組處于降額運行。圖5（b）顯示了風電棄風運行懲罰成本系數(shù)變化對優(yōu)化后總成本的影響，可以看出，優(yōu)化后總成本隨著風電棄風運行懲罰成本系數(shù)的增加而降低，這充分體現(xiàn)了風電機組的環(huán)境價值，確保風電機組的優(yōu)先發(fā)展。

圖5 風電棄風運行懲罰成本系數(shù)變化對機組出力的影響Fig.5 Influence of WAOPC coefficient on unit output

3.5 所建動態(tài)優(yōu)化模型仿真分析

本文采用文獻［7］中式（4）、（5）求取的風電場功率期望值對所建動態(tài)優(yōu)化模型進行求解，風電場功率期望值計算結果共包含5個時段，各時段的計算結果分別是 40、28、17、0、10 MW。設火電機組名義環(huán)境補償成本系數(shù)為1，風電備用容量補償成本系數(shù)為2，風電運營成本系數(shù)為2.2，仿真風電棄風運行懲罰成本系數(shù)對機組出力的影響，仿真結果如表1所示，表中有功功率為標幺值，后同。為使風電機組的優(yōu)化出力接近于期望功率，應調整風電棄風運行懲罰成本系數(shù)，當風電期望功率的值比較小時，應降低風電棄風運行懲罰成本系數(shù)以減少風電機組的優(yōu)化出力；當風電期望功率的值比較大時，應提高風電棄風運行懲罰成本系數(shù)以增加風電機組的優(yōu)化出力。另外，當風電期望功率的值比較大并且調整風電棄風運行懲罰成本系數(shù)以使風電機組按期望功率發(fā)電時，優(yōu)化后的發(fā)電總成本有明顯的降低，反之，優(yōu)化后的發(fā)電總成本大幅增加，這充分體現(xiàn)了風電機組的環(huán)境價值。

表1 考慮風電棄風運行懲罰成本系數(shù)的動態(tài)潮流仿真結果Table 1 Simulative results of dynamic OPF with WAOPC

在其他成本系數(shù)保持不變、風電棄風運行懲罰成本系數(shù)取為0的情況下，對所建動態(tài)優(yōu)化模型進行仿真求解，如表2所示。由表2可知，所有時段的風電機組的優(yōu)化出力低于期望功率，除了時段1，其余時段的風電機組的優(yōu)化出力為0，這是備用容量補償成本系數(shù)和風電運營成本系數(shù)共同作用的結果，時段1的調度出力為0.2586p.u.，從而使該時段優(yōu)化后總成本最低。更重要的是，在忽略風電棄風運行懲罰成本系數(shù)的情況下，在所有5個時段優(yōu)化后的總成本比考慮風電棄風運行懲罰成本系數(shù)時有明顯的增加，可以看出，風電棄風運行懲罰成本的引入使風電獲得了優(yōu)先調度權，充分利用了風電機組的環(huán)境價值。

表2 忽略風電棄風運行懲罰成本系數(shù)的動態(tài)潮流仿真結果Table 2 Simulative results of dynamic OPF without WAOPC

為了進一步驗證所建模型的有效性，在將火電機組名義環(huán)境補償成本和風電棄風運行懲罰成本系數(shù)設為0、其他成本系數(shù)保持不變的情況下，對所建模型進行靜態(tài)仿真求解，所有時段的各臺機組優(yōu)化結果相同，為 0、0.8、0.4、0、0.3、0.4196p.u.。對于不同時段，風電機組的優(yōu)化出力都為0，低于期望功率，這是備用容量補償成本系數(shù)和風電運營成本系數(shù)共同作用導致風電不參與優(yōu)化調度，從而火電機組在不同時段具有相同的仿真結果。更重要的是，5個時段分別進行靜態(tài)優(yōu)化后的總成本比上述2種動態(tài)仿真情況下有明顯的增加，為730.3768$/h，這進一步驗證了所建模型的有效性，體現(xiàn)了風電的環(huán)境價值。

4 結論

由于風速具有隨機變化的特點，本文首先計算每個時段的風電場有功輸出的期望值；在此基礎上，引入火電名義環(huán)境補償成本和極限穿透功率的風電場棄風運行懲罰成本來考慮環(huán)境效益；同時，引入風電備用成本來處理隨著系統(tǒng)風電穿透水平增加而引起的系統(tǒng)備用上升，由此建立了大規(guī)模異步風電并網(wǎng)的動態(tài)優(yōu)化潮流新模型。對IEEE30節(jié)點系統(tǒng)進行修改測試，通過仿真得出以下結論。

a.由于籠型異步風力發(fā)電機組具有特殊的電壓-無功特性，易于建立動態(tài)優(yōu)化潮流新模型。

b.對環(huán)境污染大的機組調度權隨名義環(huán)境補償成本系數(shù)的提高而降低，對環(huán)境更友好的機組隨名義環(huán)境補償成本系數(shù)的提高獲得優(yōu)先調度權而處于滿發(fā)狀態(tài)，由此可以看出，名義環(huán)境補償成本不僅改變了火電機組的調度權，而且賦予風電等環(huán)境友好的可再生能源調度優(yōu)先權，確保風電的優(yōu)先發(fā)展。

c.對不同的時段，為使風電機組的實際出力按預測值進行調度，應調節(jié)風電棄風運行懲罰成本系數(shù)，當風電期望功率的值比較大時，應提高風電棄風運行懲罰成本系數(shù)以增加風電機組的優(yōu)化出力，使得優(yōu)化后總成本減少；反之，降低風電棄風運行懲罰成本系數(shù)以減少風電機組的優(yōu)化出力，使得優(yōu)化后總成本增加。

該模型既考慮傳統(tǒng)電力生產(chǎn)所帶來的環(huán)境污染和資源消耗等外部成本，同時也計及因棄風行為造成的風電損失，可防止調度部門以穿透功率為借口導致棄風行為的發(fā)生，使風電一直按期望值處于滿發(fā)狀態(tài)，而其他能耗大的機組一直處于降額運行狀態(tài)，使以風電為代表的新能源環(huán)境價值得到合理的反映。

［1］李生虎.風力電力系統(tǒng)分析［M］.北京：科學出版社，2012：1-5.

［2］周雙喜，魯宗相.風力發(fā)電與電力系統(tǒng)［M］.北京：中國電力出版社，2011：1-12，390-424.

［3］周瑋，彭昱，孫輝，等.含風電場的電力系統(tǒng)動態(tài)經(jīng)濟調度［J］.中國電機工程學報，2009，29（25）：13-18.ZHOU Wei，PENG Yu，SUN Hui，etal.Dynamiceconomic dispatch in wind power integrated system［J］.Proceedings of the CSEE，2009，29（25）：13-18.

［4］江岳文，溫步瀛.結合風電功率超短期預測值偏差的實時市場調度［J］.電力自動化設備，2015，35（3）：12-17.JIANG Yuewen，WEN Buying.Real-time market dispatch based on ultra-short-term forecast error of wind power［J］.Electric Power Automation Equipment，2015，35（3）：12-17.

［5］王成福，梁軍，張利，等.基于機會約束規(guī)劃的風電預測功率分級處理［J］.電力系統(tǒng)自動化，2011，35（17）：14-19.WANG Chengfu，LIANG Jun，ZHANG Li，et al.Classified treatment of wind power predictive power based on chance constrained programming［J］.Automation of Electric Power Systems，2011，35（17）：14-19.

［6］孔波利，崔麗艷，丁釗，等.基于風光混合模型的短期功率預測方法研究［J］.電力系統(tǒng)保護與控制，2015，43（18）：62-66.KONG Boli，CUILiyan，DING Zhao，etal.Shortterm power prediction based on hybrid wind-PV forecasting model［J］.Power System Protection and Control，2015，43（18）：62-66.

［7］陳金富，陳海焱，段獻忠.含大型風電場的電力系統(tǒng)多時段動態(tài)優(yōu)化潮流［J］.中國電機程學報，2006，26（3）：31-35.CHEN Jinfu，CHEN Haiyan，DUAN Xianzhong.Multi-period dynamic optimal power flow in wind power integrated system［J］.Proceedings of the CSEE，2006，26（3）：31-35.

［8］劉學，李暉，周明，等.考慮風速相關性的多目標電網(wǎng)規(guī)劃［J］.電力自動化設備，2015，35（10）：87-94.LIU Xue，LI Hui，ZHOU Ming，et al.Multiple-objective power grid planning considering wind speed correlation ［J］.Electric Power Automation Equipment，2015，35（10）：87-94.

［9］FRANCOIS B，GALIANA F D.Stochastic security for operations planning with significant wind power generation［J］.IEEE Transactions on Power Systems，2008，23（2）：306-316.

［10］周雙喜，王海超，陳壽孫.風力發(fā)電運行價值分析［J］.電網(wǎng)技術，2006，30（14）：98-102.ZHOU Shuangxi，WANG Haichao，CHEN Shousun.Analysis on operation value of wind power resources ［J］.Power System Technology，2006，30（14）：98-102.

［11］HETZER J，YU D C，BHATTARAI K.An economic dispatch model incorporating wind power［J］.IEEE Transactions on Energy Conversion，2008，23（2）：603-611.

［12］雷亞洲，王偉勝，印永華，等.風電對電力系統(tǒng)運行的價值分析［J］.電網(wǎng)技術，2002，26（5）：10-14.LEI Yazhou，WANG Weisheng，YIN Yonghua，etal.Analysis of wind powervaluetopowersystem operation［J］.Power System Technology，2002，26（5）：10-14.

［13］張海峰，高峰，吳江，等.含風電的電力系統(tǒng)動態(tài)經(jīng)濟調度模型［J］.電網(wǎng)技術，2013，37（5）：1128-1303.ZHANG Haifeng，GAO Feng，WU Jiang，etal.A dynamic economic dispatching model for power grid containing wind power generation system［J］.Power System Technology，2013，37（5）：1128-1303.

［14］聶永輝，李天云，劉鳳蘭，等.基于非線性內點理論的交直流系統(tǒng)無功優(yōu)化新模型［J］.電力系統(tǒng)保護與控制，2014，42（10）：35-40.NIE Yonghui，LI Tianyun，LIU Fenglan，et al.A new AC-DC reactive optimal power flow model based on nonlinear interior point method［J］.Power System Protection and Control，2014，42（10）：35-40.

［15］孫國強，任賓，衛(wèi)志農(nóng)，等.基于自適應加權預測-校正內點法的含VSC-HVDC電力系統(tǒng)最優(yōu)潮流［J］.電力自動化設備，2015，35（3）：54-59，67.SUN Guoqiang，REN Bin，WEI Zhinong，et al.Optimal power flow of power system with VSC-HVDC based on AWPC-IPM［J］.Electric PowerAutomation Equipment，2015，35 （3）：54-59，67.