風電系統參數對Hopf分岔影響的仿真

李 季, 周雪松, 馬幼捷, 王海洋

(1. 天津大學 電氣與自動化工程學院, 天津 300072；2. 天津理工大學 天津市復雜系統控制理論及應用重點實驗室, 天津 300384; 3. 白城師范學院 機械工程學院, 吉林 白城 137000)

隨著風力發電風電容量的迅速增加, 風電并網對系統的影響日益顯著, 但風電大規模接入電網會導致電能質量下降、 線路傳輸功率越限、 系統短路容量增加和暫態穩定性變化等問題[1-3]. 文獻[4]指出風速變化會影響并網風電系統的電能質量, 引起系統電壓的波動； 文獻[5]指出風電機組對風電系統電壓的影響程度決定于電網結構強度和風機容量大小.

風電的大規模接入必然會使系統運行點向穩定極限點靠近, 一對共軛特征根出現于虛軸而產生的Hopf分岔是該特征的一種現象. Hopf分岔是常見的分岔現象[6], 是導致電壓失穩的3種分岔形式之一[7]. 分岔理論已成功解釋了DAE(differential algebraic equation)方程所表示的非線性動態電力系統電壓失穩機理[8-10]. 文獻[11]研究了含風電場的電力系統靜態電壓穩定問題； 文獻[12]基于延拓法追蹤了基于風電系統靜態模型的二維參數分岔邊界. 但已有研究很少涉及風電系統DAE方程的Hopf分岔. 本文應用延拓算法追蹤風電系統的平衡解曲線及兩參數Hopf分岔邊界, 分析參數對風電系統電壓穩定性的影響, 并研究靜止無功補償器(static var compensation, SVC)對Hopf分岔的控制作用.

1 延拓算法

對于動力學系統[13]

(1)

其中：x表示狀態變量；λ表示分岔參數.

設式(1)在平衡點(x,λ)出現Hopf分岔, 則該系統的Jacobi矩陣A=?f/?x在平衡點處有一對共軛純虛根, 記為jω[13]. 對Hopf分岔, 有

f(x,λ)=0,Aq=jωq, 〈q,q0〉=1,

(2)

式(2)為復數, 即向量q和q0,x,λ,q,ω為方程的未知數, 方程的個數與未知數個數相等.

設式(1)的平衡點方程為f(x,λ)=0, 則該方程在n+1維的空間y=(x,λ)∈Rn×R1中定義了一個一維流形M, 稱為平衡解曲線或平衡解流形, 延拓算法可追蹤該平衡解流形, 即用一系列離散點(y1,y2,…)逼近曲線M, 該離散點滿足平衡點的要求. 追蹤該平衡流形的過程中, 通過檢測局部分岔條件判斷在該流形上是否存在分岔點. 但一般的延拓法僅能追蹤到一維流形, 因此只能計算單參數分岔值. 為了實現微分代數方程表示的風電系統二維參數分岔邊界的追蹤. 可從單參數分岔點出發, 假設已追蹤得到式(1)的單參數分岔點, 據此可探索將延拓算法應用到求解直接算法所能求解的數學方程中. 即等價于應用延拓算法追蹤局部分岔滿足的流形, 并非上述意義下的平衡解流形, 這時可求得式(1)的二維參數局部分岔邊界.

2 風電系統模型

風電系統模型如圖1所示. 該模型由2個發電機和1個風電場組成, 發電機G1為無窮大電源.

圖1 風電系統模型Fig.1 Model of wind power system

2.1 等值發電機模型

等值發電機G2采用下述二階模型:

其中：Em,M,Tm和Dm分別表示發電機機端電壓、 發電機的慣量、 輸入轉矩及阻尼系數；ω和δm分別表示發電機角頻率和發電機功角；ym和am表示網絡參數.

2.2 風電場模型

異步發電機技術較成熟、 性價比高、 運行可靠, 是目前多數風電場的主流機型[14], 異步發電機運行時發出有功功率, 同時吸收無功功率, 使整個風電場對無功功率需求較大, 降低了風電接入地區電網的電壓穩定性.

異步發電機采用忽略定子磁通暫態過程, 考慮轉子繞組暫態過程的暫態模型.

1) 定子電壓方程為

u=(r1+jx′)Is+E′.

(3)

2) 電磁暫態方程為

(4)

3) 轉子運動方程為

(5)

3 風電系統參數對Hopf分岔的影響

3.1 有功功率和無功功率對Hopf分岔的影響

Hopf分岔是最具代表性的動態分岔, 也是其他動態分岔的基礎. Hopf分岔和電力系統的參數密切相關, 不同參數對Hopf分岔的影響也不同. 下面針對幾個風電系統的典型參數對Hopf分岔影響進行分析.

由于風速的隨機性和不確定性, 感應發電機消耗或吸收的無功功率也會不斷變化, 導致電網電壓的波動, 引起風電發電機組機端電壓下降, 因而極大降低了系統電壓的穩定水平. 以風電場吸收的無功功率Q1作為分岔參數時, 風電場機端電壓u的Hopf分岔曲線如圖2所示, 本文僅分析穩定平衡態的上半分支. 當無功功率Q1逐漸增加時, 風電場電壓u快速降落, 在系統運行到極限點(limit point, LP)前的某時刻,Q1約達到1.499 884, 風電系統發生了Hopf分岔, 此時電壓u=0.908 584, 如圖2中的H1, Hopf分岔理論指出, 當系統發生Hopf分岔時, 系統出現振蕩, 進而維持周期性振蕩或振幅不斷增大而導致最終失穩, 因此風電場消耗的無功功率將影響風電系統動態電壓的穩定性.

風速的不斷變化導致感應風電機組在運行中有功功率和無功功率同時變化, 二者相互影響, 因此針對風電場的兩參數Hopf分岔分析更有意義. 如圖3所示,P1和Q1為分岔參數風電系統電壓u的二維Hopf分岔邊界, 隨著風電機組發出有功功率P1的增加, 其所吸收的無功功率Q1也會相應增加, 此時系統電壓u下降. 由圖2可見, 風電系統無功功率的增加易導致Hopf分岔發生, 所以圖3給出了避免Hopf分岔發生, 提高電壓穩定性的規律, 即隨著風電場發出有功功率的增大, 必須減小風電場吸收的無功功率. 相應地, 如果風電場的無功功率增加, 則必須減少風電場的有功功率, 這樣雖有效避免了系統運行在Hopf的二維分岔邊界點上, 但Hopf分岔的消失以犧牲風電場的有功功率為代價. 所以控制風電場所吸收的無功功率是提高系統電壓穩定性的一種有效方式.

圖2 Q1變化時風電系統的Hopf分岔曲線Fig.2 Hopf bifurcation curve of wind power system with Q1 variation

圖3 P1-Q1的二維Hopf分岔邊界Fig.3 Two-dimensional Hopf bifurcation boundary of P1-Q1

3.2 傳輸線路導納對Hopf分岔的影響

圖4為y0變化時的風電系統Hopf分岔曲線. 由圖4可見, 線路的導納參數對系統電壓水平影響較小. 圖5給出了系統導納參數y0和ym共同作用下二維Hopf分岔邊界. 由圖4和圖5可見, 隨著線路導納的增加, 系統分岔邊界電壓水平變化不大. 當y0=5.626 611時, 系統發生Hopf分岔, 如圖4的H1, 當增大傳輸線路導納y0時, 即縮短了線路的電氣距離, 降低了線路所需的無功功率, 相應地, 系統供給風電場的無功功率也會增加, 因而使系統的Hopf分岔延后, 提高了系統電壓的穩定性.

圖4 y0變化時風電系統的Hopf分岔曲線Fig.4 Hopf bifurcation curve of wind power system with y0 variation

圖5 y0-ym的二維Hopf分岔邊界Fig.5 Two-dimensional Hopf bifurcation boundary of y0-ym

4 靜止無功補償器對Hopf分岔的影響

Hopf分岔可導致風電系統電壓不穩定, 所以控制系統中的Hopf分岔能避免電壓崩潰. 由上面分析可知, 系統Hopf分岔和無功功率關系密切, 而異步風電機組通常只能吸收無功, 無功控制主要依靠電容器就地補償, 無法實現快速的無功調節. 靜止無功補償器(SVC)作為靈活的動態無功補償裝置, 在實際應用中可有效控制分岔[12].

本文采用在風電系統的風電場端安裝靜止無功補償器有效控制無功功率, 如圖6所示.

圖6 加裝SVC的風電系統模型Fig.6 Wind power system with SVC

數學模型為

(6)

其中：u表示補償點電壓；uref表示參考電壓；Tr表示時間常數；Kr表示控制器增益；bSVC表示輸出電壓； 無功功率Q=bSVCu2.

圖7 加裝SVC的風電系統Hopf分岔曲線Fig.7 Hopf bifurcation curve of wind power system with SVC

以Q1對風電系統的Hopf分岔影響為例, 分析靜止無功補償器對系統Hopf分岔及電壓穩定性的影響. 當Kr=1.5,Tr=0.02時, 求得平衡點處風電場端電壓數值隨分岔參數Q1的變化情況如圖7所示. 與圖2比較可見, 含SVC的風電系統,Q1增加時, 提高了風電場的電壓水平, 且系統的Hopf分岔點H1無限接近系統的LP點, 對應的無功功率也提高到Q1=11.531 989. 表明靜止無功補償器補償了風電場吸收的無功功率, 同時延遲了Hopf分岔, 擴大了風電系統的電壓穩定域.

綜上所述, 本文應用延拓算法求取了風電系統的平衡解流形, 并針對幾種分岔參數和靜止無功補償器對風電系統的Hopf分岔影響進行了研究, 進而分析了其對風電系統電壓穩定性的影響： 1) 隨著風電場發出有功功率的增加, 需減小風電場吸收的無功功率, 既避免了系統運行在Hopf分岔邊界上, 又提高了系統電壓水平; 2) 風電系統線路導納的增加, 使風電系統的Hopf分岔延后, 增加了系統穩定域, 提高了風電系統的電壓穩定性; 3) 靜止無功補償器為風電場有效提供了無功功率, 可延遲系統的Hopf分岔點, 增加負荷極限, 提高風電系統的電壓穩定性.

[1] LIN Li, SUN Cai-xin, WANG Yong-ping, et al. Calculation Analysis and Control Strategy for Voltage Stability of Power Grid with Large Capacity Wind Farm Interconnected [J]. Power System Technology, 2008, 32(3)： 41-46. (林莉, 孫才新, 王永平, 等. 大容量風電場接入后電網電壓穩定性的計算分析與控制策略 [J]. 電網技術, 2008, 32(3)： 41-46.)

[2] Vittal E, O’Malley M, Keane A. A Steady-State Voltage Stability Analysis of Power Systems with High Penetrations of Wind [J]. IEEE Transactions on Power Systems, 2010, 25(1)： 433-442.

[3] CHEN Ning, ZHU Ling-zhi, WANG Wei. Strategy for Reactive Power Control of Wind Farm for Improving Voltage Stability in Wind Power Integrated Region [J]. Proceedings of the CSEE, 2009, 29(10)： 102-108. (陳寧, 朱凌志, 王偉. 改善接入地區電壓穩定性的風電場無功控制策略 [J]. 中國電機工程學報, 2009, 29(10)： 102-108.)

[4] ZHANG Feng, CHAO Qin. Research on Improving Transient Voltage Stability of Wind Farm by STATCOM [J]. Power System Technology, 2008, 32(9)： 70-73. (張鋒, 晁勤. STATCOM改善風電場暫態電壓穩定性的研究 [J]. 電網技術, 2008, 32(9)： 70-73.)

[5] WEI Xiao-guang, CHI Yong-ning, WEI Xiao-yun, et al. Transient Voltage Stability Enhancement of Induction Generator-Based Wind Farm by VSC-HVDC [J]. Transactions of China Electrotechnical Society, 2008, 23(5)： 100-107. (魏曉光, 遲永寧, 魏曉云, 等. 電壓源換相高壓直流輸電改善異步發電機風電場暫態電壓穩定性的研究 [J]. 電工技術學報, 2008, 23(5)： 100-107.)

[6] CHEN Yu-dong, PEI Chun-yan, CAO Zong-jie. Stability Analysis of Nonlinear System at Critical Point of Hopf Bifurcation [J]. Journal of Jilin University: Engineering and Technology Edition, 2009, 39(4): 993-996. (陳宇東, 裴春艷, 曹宗杰. 非線性系統Hopf分叉臨界點的穩定性 [J]. 吉林大學學報: 工學版, 2009, 39(4): 993-996.)

[7] Radunskaya A, Williamson R, Yinger R. A Dynamic Analysis of the Stability of a Network of Induction Generators [J]. IEEE Transactions on Power Systems, 2008, 23(2)： 657-663.

[8] Kataoka Yoshihiko, Shinoda Yukio. Voltage Stability Limit of Electric Power Systems with Generator Reactive Power Constraints Considered [J]. IEEE Transactions on Power Systems, 2005, 20(2)： 951-962.

[9] LI Hong-zhong, CHENG Hao-zhong, TENG Le-tian, et al. A Direct Method for Computing Hopf Bifurcation Point in Power System Dynamic Voltage Stability [J]. Proceedings of the CSEE, 2006, 26(8)： 28-32. (李宏仲, 程浩忠, 滕樂天, 等. 以簡化直接法求解電力系統動態電壓穩定Hopf分岔點 [J]. 中國電機工程學報, 2006, 26(8)： 28-32.)

[10] JIANG Wei, WANG Cheng-shan, YU Yi-xin, et al. A New Method for Direct Calculation the Critical Point of Static Voltage Stability [J]. Proceedings of the CSEE, 2006, 26(10): 1-6. (江偉, 王成山, 余貽鑫, 等. 直接計算靜態電壓穩定臨界點的新方法 [J]. 中國電機工程學報, 2006, 26(10): 1-6.)

[11] MA You-jie, ZHANG Ji-dong, ZHOU Xue-song, et al. Study on Steady State Voltage Stability of Power System Containing Wind Farm Based on Bifurcation Theory [J]. Power System Technology, 2008, 32(9)： 74-79. (馬幼捷, 張繼東, 周雪松, 等. 基于分岔理論的含風電場電力系統靜態電壓穩定問題研究 [J]. 電網技術, 2008, 32(9)： 74-79.)

[12] MA You-jie, WEN Hu-long, ZHOU Xue-song, et al. Calculation and Research of Two-Dimensional Parameter Bifurcation Boundary of Power System Voltage Stability Model with Wind Farm by Continuation Method [J]. Proceedings of the CSEE, 2011, 30(19)： 26-30. (馬幼捷, 問虎龍, 周雪松, 等. 采用延拓法的風電系統穩定模型二維參數分岔邊界的計算與研究 [J]. 中國電機工程學報, 2011, 30(19)： 26-30.)

[13] GUAN Qing-guang, ZANG Lin. Numerical Computation of Periodic Solutions Generated by Generalized Hopf Bifurcation of 3-Dimensional Piecewise Smooth Filippov-Type Equation [J]. Journal of Jilin University: Science Edition, 2010, 48(3): 371-374. (管慶光, 臧林. 三維分片光滑Filippov-型方程Hopf分支周期解的數值計算 [J]. 吉林大學學報: 理學版, 2010, 48(3): 371-374.)

[14] ZHANG Rui-ming, ZHANG Xin-yan. Study on Static Voltage Stability of Power System Based on Static Models of Common and Doubly Fed Asynchronous Generators [J]. Power System Technology, 2011, 35(1)： 175-179. (張瑞明, 張新燕. 基于普通異步發電機和雙饋風力發電機靜態數學模型的系統靜態電壓穩定性研究 [J]. 電網技術, 2011, 35(1)： 175-179.)