【正文】 后，通過仿真，得到um(w(k))、fm(w(k))，通過（）～（）式求得Xu(w(k)) 、Xf(w(k))和w(k+1)；3．若Xu(w(k))＞1或Xf(w(k))＞1，則令： ()如果△W(k)≤1MW,停止迭代，取上次迭代計(jì)算結(jié)果為最終解；4．如果Xu(w(k))=1且Xf(w(k))＜1，或Xf(w(k))=1且Xu(w(k))＜1，則W(k)為最終解；5．重復(fù)2～4，直到求得最終解。算例采用IEEE(New England)39節(jié)點(diǎn)系統(tǒng)，對(duì)上述的風(fēng)電場極限功率求取及優(yōu)化做了驗(yàn)證。算例系統(tǒng)中，風(fēng)電場通過節(jié)點(diǎn)8接入系統(tǒng)。系統(tǒng)中除風(fēng)電場外其他發(fā)電廠總裝機(jī)容量為7050MW，系統(tǒng)最大負(fù)荷為6097MW。假設(shè)風(fēng)電場初始容量為300MW，應(yīng)用本文提出的優(yōu)化算法，計(jì)算三種運(yùn)行方式下，風(fēng)電場的穿透功率極限。 算例IEEE(New England)39節(jié)點(diǎn)系統(tǒng)主接線圖1．故障運(yùn)行方式。在風(fēng)電場運(yùn)行時(shí)，由于故障導(dǎo)致風(fēng)電場與系統(tǒng)解列。仿真次數(shù)（K）風(fēng)電場容量（MW）節(jié)點(diǎn)電壓最小值（pu）系統(tǒng)頻率最小值（pu）風(fēng)電場容量修正值（MW）10100300230034512．風(fēng)電場輸出功率短時(shí)間內(nèi)迅速上升?？紤]風(fēng)電場在短時(shí)間(4秒鐘)內(nèi)，輸出功率由零增至額定值。仿真次數(shù)（K）風(fēng)電場容量（MW）節(jié)點(diǎn)電壓最大值（pu）系統(tǒng)頻率最大值（pu）風(fēng)電場容量修正值（MW）101003002300345613．風(fēng)電場輸出功率短時(shí)間內(nèi)迅速下降。考慮風(fēng)電場在短時(shí)間(4秒鐘)內(nèi)，輸出功率由額定值降至零。仿真次數(shù)（K）風(fēng)電場容量（MW）節(jié)點(diǎn)電壓最小值（pu）系統(tǒng)頻率最小值（pu）風(fēng)電場容量修正值（MW）101003002300345671綜合上述仿真結(jié)果可知，風(fēng)電場通過節(jié)點(diǎn)8接入系統(tǒng)時(shí)，%。通過上述三種運(yùn)行方式的仿真計(jì)算結(jié)果可以看出，在第3種運(yùn)行方式下，系統(tǒng)因失去大量有功導(dǎo)致系統(tǒng)節(jié)點(diǎn)電壓、頻率降低。在第2種運(yùn)行方式下，系統(tǒng)因短時(shí)間內(nèi)有功激增，導(dǎo)致系統(tǒng)節(jié)點(diǎn)電壓、頻率上升。改變風(fēng)電場接入系統(tǒng)的節(jié)點(diǎn)，通過上述方法計(jì)算風(fēng)電場最大容量。風(fēng)電場接入點(diǎn)風(fēng)電場最大容量（MW）仿真迭代次數(shù)平均值風(fēng)電場穿透功率極限故障運(yùn)行方式風(fēng)電場輸出功率劇減風(fēng)電場輸出功率激增2%3%4%5%7%當(dāng)風(fēng)電場接入節(jié)點(diǎn)選擇不同時(shí)，風(fēng)電場的最大容量、風(fēng)電場穿透功率極限也不同。一般來說，風(fēng)電場接入系統(tǒng)節(jié)點(diǎn)離發(fā)電廠越近時(shí)，它的穿透功率極限會(huì)越小，因?yàn)榇藭r(shí)它對(duì)系統(tǒng)頻率影響較大。上述仿真結(jié)果驗(yàn)證了優(yōu)化算法的真實(shí)準(zhǔn)確性。應(yīng)用此方法免除了人為的去估算風(fēng)電場最大容量的修正值，使得風(fēng)電場極限穿透功率求取簡便而快速。本節(jié)在應(yīng)用動(dòng)態(tài)仿真求取風(fēng)電場極限功率的過程中，通過近似線性規(guī)劃優(yōu)化方法對(duì)風(fēng)電場容量值的修正進(jìn)行優(yōu)化。驗(yàn)證結(jié)果表明，此種優(yōu)化方法計(jì)算準(zhǔn)確而迭代快速，它可與不同的仿真軟件結(jié)合，快速準(zhǔn)確的求取系統(tǒng)風(fēng)電場穿透功率極限。5總結(jié)與展望本文比較全面和深入地分析了大型風(fēng)電場并網(wǎng)運(yùn)行對(duì)系統(tǒng)正常運(yùn)行的影響，了解到風(fēng)電系統(tǒng)并網(wǎng)對(duì)電力系統(tǒng)電壓、頻率等的影響。進(jìn)而簡單介紹了一系列計(jì)算風(fēng)電場極限穿透功率的方法，由于動(dòng)態(tài)仿真求取風(fēng)電場穿透功率極限時(shí)，如若考慮到系統(tǒng)各種運(yùn)行方式的話，計(jì)算量將會(huì)很大，為了能夠減少計(jì)算量，本文提出了一種基于近似線性規(guī)劃的風(fēng)電場穿透功率極限優(yōu)化的改進(jìn)算法，通過近似線性規(guī)劃方法，將原計(jì)算風(fēng)電場穿透功率極限的非線性目標(biāo)約束函數(shù)作線性化再應(yīng)用線性近似解去逼近非線性真實(shí)解。應(yīng)用該方法可以快速準(zhǔn)確地求取風(fēng)電場穿透功率極限。通過在New England 39節(jié)點(diǎn)系統(tǒng)上仿真計(jì)算，驗(yàn)證了所提方法的有效性及其快速準(zhǔn)確的特點(diǎn)，得到了風(fēng)電場接入節(jié)點(diǎn)選擇不同時(shí)，風(fēng)電場的最大容量、風(fēng)電場穿透功率極限也不同的結(jié)論。一般來說，風(fēng)電場接入系統(tǒng)節(jié)點(diǎn)離發(fā)電廠越近時(shí)，它的穿透功率極限會(huì)越小，因?yàn)榇藭r(shí)它對(duì)系統(tǒng)頻率影響較大。通過對(duì)風(fēng)電場穿透功率極限的分析研究的方法優(yōu)化，來進(jìn)一步研究解決以下問題：、電源規(guī)劃情況，研究風(fēng)電穿透功率對(duì)電網(wǎng)調(diào)度（發(fā)電）計(jì)劃、經(jīng)濟(jì)運(yùn)行、系統(tǒng)調(diào)節(jié)能力及控制手段的影響；．對(duì)電網(wǎng)風(fēng)電容量可信度進(jìn)行正確評(píng)價(jià)，在對(duì)風(fēng)電容信度正確評(píng)價(jià)的基礎(chǔ)上深入研究風(fēng)電穿透功率極限值；。從電力系統(tǒng)運(yùn)行的角度出發(fā)，可以在上述約束條件中考慮風(fēng)電價(jià)值與風(fēng)電場所在地的風(fēng)資源特性、風(fēng)電機(jī)組技術(shù)參數(shù)、風(fēng)電負(fù)荷在系統(tǒng)負(fù)荷中所占比例之間的關(guān)系等；，解決風(fēng)電并網(wǎng)時(shí)系統(tǒng)的無功功率控制，頻率控制，電壓波動(dòng)等問題。致 謝在論文完成之際，我首先要向陳功貴老師表示最真摯的謝意。在論文寫作期間，陳老師給我悉心的指導(dǎo)，讓我得到了多方面的關(guān)懷和幫助，老師淵博的知識(shí)和嚴(yán)謹(jǐn)?shù)目茖W(xué)態(tài)度及崇高的敬業(yè)精神使我受益匪淺，在此謹(jǐn)向他表示我最衷心的感謝和敬意。向所有關(guān)心和幫助過我的老師、同學(xué)和朋友表示真誠的感謝。最后，向評(píng)閱論文的各位老師表示最誠摯的謝意。由于本人學(xué)識(shí)有限，文中不免有錯(cuò)誤和待改進(jìn)之處，真誠歡迎各位師長、同行提出寶貴意見。參考文獻(xiàn)[1] 趙旺初.“十五”風(fēng)力發(fā)電目標(biāo)已確定[J].汽輪機(jī)技術(shù)，8(13)：615.[2] , . Alternate Energy Sources and Energy Conversation[C] CIGRE, Paris, 2002, 54207.[3] Christensen J. Falck. Methods and models for evaluating the impact of decentralized generation[C]. CIGRE, Paris, 1998, 38301.[4] ，新疆風(fēng)能網(wǎng)站，.[5] Falck C. Methods and models for evaluating the impact of decentralized generation [C].CIGRE, Paris, 1998, 38301.[6] Henrik Bindner, Per Lundsager. Integration of wind power in the power system. IEEE 2002 28th Annual Conference, Volume: 4, 2002, pp: 33093316.[7] P. Christiansen. Wind farm controller and dynamic stability. Proc. EWEC2001, Copenhagen Denmark, July 2001.[8] 邰圣平，趙天洪，陳思寧，何云，[J].電力系統(tǒng)自動(dòng)化，1998,25(41)：2534.[9] P. Sorensen et al. Models for simulation of transient events in a wind farm[C]. Global Wind Energy Conference, Paris France, March 2002.[10] , , . Wind models for prediction of power fluctuations from wind farms [J]. JWE , Kyoto Japan, 2001，14(11)：518.[11] 吳學(xué)光，張學(xué)成，印永華，[J].電網(wǎng)技術(shù)，,31(22)：6278.[12] 陳樹勇，戴慧珠，白曉民，[J].電機(jī)工程學(xué)報(bào)，21(7)：5674.[13] 吳學(xué)光，王偉勝，[J].風(fēng)力發(fā)電，,16(43)：2328.[14] 著，施鵬飛 [M].. [15] 申洪，梁軍，[J].電網(wǎng)技術(shù)，2002,26(8),pp:811.[16] 雷亞洲，王偉勝，印永華，等．基于機(jī)會(huì)約束規(guī)劃的風(fēng)電穿透功率極限計(jì)算[J]．中國電機(jī)工程學(xué)報(bào),2002,22(5):3235.[17] 吳俊，李國杰，孫元章．基于隨機(jī)規(guī)劃的并網(wǎng)風(fēng)電場最大注入功率計(jì)算．電網(wǎng)技術(shù)，2007，14(31)：1519.[18] 劉建新，牛永志，[J].華北電力大學(xué)學(xué)報(bào)，2008，35(6)：5256.[19] 陳衍．電力系統(tǒng)穩(wěn)態(tài)分析（第三版）[M]．北京：中國電力出版社，2007.[20] Papadopoulos M, et al. Simulation and Analysis of small and medium size power systems containing wind turbines [J]. IEEE Transactions on Power Apparatus and Systems, 1991, 6(4): 14531458.[21] Zaininger H W. Potential dynamic impacts of wind turbines on unity systems [J]. IEEE Transactions on Power Apparatus and Systems, 1981, 100(12): 48214829.[22] 雷亞洲，王偉勝，印永華，[J].中國電機(jī)工程學(xué)報(bào)，,: 3235.[23] 雷亞洲，王偉勝，[J].中國電機(jī)工程學(xué)報(bào)，2001,21(6):2528.[24] 陳衍．電力系統(tǒng)穩(wěn)態(tài)分析[M]．北京：中國電力出版社，1995.[25] 陳樹勇，戴慧珠，白曉民，[J].中國電力，2000，3 (11):3031.[26] Wu Xueguang, etc. Application of models of the wind energy conversion system to wind power dynamic analysis[C]. International Conference on Power System Technology, 1998, Volume: 2, pp: 14061411.[27] 王子若，陳永昌．優(yōu)化計(jì)算方法[M]．北京：機(jī)械工業(yè)出版社，1989. 附 錄IEEE (New England) 39節(jié)點(diǎn)系統(tǒng)參數(shù)：表1 節(jié)點(diǎn)負(fù)荷數(shù)據(jù)所在節(jié)點(diǎn)名恒有功（pu）恒無功（pu）所在節(jié)點(diǎn)名恒有功（pu）恒無功（pu）3322234500184247842522485221762613917128827281153201532820616329291815830312068010339110425021274115表2 發(fā)電機(jī)數(shù)據(jù)所在節(jié)點(diǎn)名有功（有名）無功（有名）電壓（pu）30 31 32 33 34 35 36 37 38 39 表3線路數(shù)據(jù)節(jié)點(diǎn)名1節(jié)點(diǎn)名2電阻（pu）電抗（pu）B/2（pu）2139132252431835414465857611687983991110131014131514161517161916211624161817271722212322242326252726282629262928表4變壓器數(shù)據(jù)節(jié)點(diǎn)名1節(jié)點(diǎn)名2電阻（pu）電抗（pu）變比1112131231632103319342035223623372530238292019