王忠勇 ,李嘯驄
1.廣東石油化工學院計算機與電子信息學院,廣東 茂名 525000;
2.廣西電力系統最優化與節能技術重點實驗室,廣西 南寧 530004
風能的利用近年來得到了巨大的發展,但風電大規模并網給電力系統的穩定運行帶來了新的挑戰,并網運行后電網的擾動或風速的變化會導致機網扭振,對風電機組的穩定運行和設備安全產生很大的影響[1-5]。在大容量風電場中,嚴重的扭振現象會造成風電機組損耗加劇甚至損壞,影響電力系統的安全運行[6-7]。因此,深入研究風電機組的動態特性及其與電網的相互作用,以及風機軸系的扭振問題是保證機組安全運行的重要課題。
直驅永磁風電傳動系統具有“柔性”,受到電網擾動時軸系出現的低頻振蕩,使傳動軸出現應力疲勞[8-12]。國內學者已經探討了風機的穩定性運行和軸系扭振的問題[13-16],已有學者對非線性控制理論應用于電力系統穩定分析進行了深入的探討[17-19]。
本文基于直驅永磁風機(direct-driven permanent magnet synchronous generators,D-PMSG )軸系扭振,建立了風機發電機組的模型,采用非線性控制策略實現協調控制來抑制軸系間扭矩振蕩,仿真結果驗證了該控制策略抑制軸系扭振的有效性。
風電機組系統采用如圖1所示的系統模型,由空氣動力學可知,風力機產生的機械功率為風力機從空氣中捕獲的風能,由(1)式表示:

式(1)中 ρ為空氣的密度,r為風力機葉片半徑,Vw為風速,Cp為風能利用系數,是γ和 β的非線性函數,β為葉片的槳距角,γ為風力機葉尖速比,ωw為風力機的轉速,P=Tωw。

圖1 直驅永磁風力發電系統結構圖Fig.1 Structure diagram of PMSG system
1.2.1 機械軸系模型 風電機的軸系通常等效為兩個或三個質量塊[4],這里采用兩質塊,數學模型為式(2):

式(2)中 D為阻尼系數,K為剛性系數,Hw、Hg為風機、發電機慣性時間常數,θw,θg分別為相對角位移,ωw、ωg為轉子角速度,Tw,Te分別為機械轉矩、電磁轉矩
1.2.2 槳距角模型 風速通過葉片槳距角β的控制來實現,如下方程表示槳距角模型[8]:

式(3)中Ts是槳距控制系統的慣性時間常數,β0是槳距角初始值。
1.2.3 D-PMSG模型 設dq坐標系以同步速度旋轉,且q軸超前于d軸,則在dq軸坐標下的D-PMSG模型為式(4):

式(4)中 u1d、u1q分別為d、q軸電壓,i1d、i1q分別為d、q軸電流;L1d、L1q分別為定子直軸電感和交軸電感;λ0為永磁體磁鏈;np為轉子的極對數;ωm及Ra分別為定子電角速度和電阻。
根據能量守恒,不計變流器的能量損失,直流環節表示為式(5):

式(5)中,m1,θ1,m2,θ2分別為變換器直流環節兩側變比和導通角。
1.2.4 電網側數學模型 由結構圖,取d軸與a相電壓矢量重合,d軸超前q軸90°電角度,則網側模型為[13]式(6):

式(6)中 u2d、u2q分別為d、q軸電壓,i2d、i2q分別為 d、q軸電流;為電網電壓,ω 為電網角速g度,L2、R分別為連接電感及等值電阻。
綜合式(2)~式(6)構成了分析直驅永磁風電機軸系扭振的數學模型,可以寫出式(7)形式[18]:

其中

根據對軸系扭振模型和控制目標的分析,確定3個輸入量分別為:

根據對風機控制做分析,將槳距角,機側、網側變流器作為控制對象。
考慮非線性仿射控制系統[16]:

式(8)中 x∈Rn為狀態量,f(x)=[f1(x)…fn(x)]T∈Rn、gj(x)=[g1j(x)…gnj(x)]T∈Rn為光滑向量,u1…um為控制標量;h1(x)…hm(x)為輸出函數;y1…ym為輸出標量。
為保證協調控制的效果,輸出函數選為多個狀態量的線性組合[19]:
槳距角 h1(x)=c11ωw+c12ωg+c13θwg+c14i1d,
機側變流器h2(x)=c23θwg+c24β+c26i1q+c27udc,
網側變流器h3(x)=c32ωg+c37udc+c39i2q。
將三個輸出函數寫成矩陣的形式:

根據非線性控制理論,先計算矩陣[22]:

經過計算,系統的相對階分別為r1+r2+r3=5<n=9,故可進行部分精確線性化設計控制律,將x空間轉化為z空間的第二標準型,

反饋控制律u可通過下式求取


對 式 (12) 的 反 饋 控 制 律 u1=Tw,u2=m1cosθ1, u3=m2sin θ2, 即可通 過 控 制 風 機 輸入,變換器的 mi(i=1,2),θi(i=1,2)來實現對扭矩的控制。
根據軸系扭振的簡化模型,在MATLAB/Simu?link環境下做仿真分析軸系扭振抑制效果,用于仿真永磁風力機組和變換器的參數采用文獻[14]數據。
設電網側發生三相短路0.5 s時發生三相短路,在0.15 s后切除故障,分析軸系扭振的兩個主要狀態量ΔδSG、ΔωSG的響應曲線如圖2~圖4所示。
曲線表明,非線性控制方法抑制軸系扭振的效果較好,偏差振蕩更小。
設定0.5 s時,風機有功功率增加10%,分析軸系扭振相關量ΔδSG、ΔωSG的變化。從曲線看出非線性控制下的振蕩更小,表明軸系扭振的幅度更小。
設定0.5 s時,機端電壓增大5%,分析扭振相關狀態量ΔδSG、ΔωSG的變化,曲線表明非線性方法可用更短的時間恢復到穩定值。

圖2 短路擾動時系統的動態響應曲線:(a)ΔωSG;(b)ΔδSGFig.2 Dynamic response curves of system during short-circuit disturbance:(a) ΔωSG ;(b)ΔδSG

圖3 功率變化時系統的動態響應曲線:(a)ΔδSG;(b)ΔωSGFig.3 Dynamic response curves of system at variable powers:(a)ΔδSG ;(b)ΔωSG

圖4 調壓擾動時系統的動態響應曲線:(a)ΔδSG;(b)ΔωSGFig.4 Dynamic response curves of system under voltage regulation disturbance:(a)ΔδSG ;(b)ΔωSG
針對直驅永磁風電機組的軸系剛體模型,結合變流器和風機功率控制建立了該軸系的數學模型,采用協調控制的方法進行多指標控制策略的設計。仿真表明,該方法可使風力發電機組受到擾動時,有效消弱軸系間的扭振,提高系統穩定運行的能力。
[1] LALOR G,MULLANE A,O’MALLEY M.Frequency controland wind turbine technologies[J]. IEEE Transactions on Power Systems,2005,20(4):1905-1913.
[2] 王曉敏,李自成,剪浩杰.基于前饋和PI控制的光伏離網逆變器的研究[J].武漢工程大學學報,2016,38(3):293-294.
[3] 楊敏霞,解大,婁宇成,等.定速異步機型風電場機網扭 振 建 模 及 仿 真[J].電 網 技 術 ,2014,38(9):2422-2428.
[4] 解大,王瑞琳,王西田,等.多機型風電機組機網扭振的 模 型 與 機 理[J]. 太 陽 能 學 報 ,2011,32(9):1281-1283.
[5] 楊秀媛,肖洋,陳樹勇.風電場風速和發電功率預測研究[J].中國電機工程學報,2005,25(11):1-5.
[6] 陳樹勇,戴慧珠,白曉民.風電場的發電可靠性模型及其應用[J].中國電機工程學報,2000,20(3):26-29.
[7] 王忠勇.基于微分幾何的電池儲能系統輸出穩定控制器 設 計[J].武 漢 工 程 大 學 學 報 ,2011,33(8):103-106.
[8] 張紅光,張粒子,陳樹勇.大容量風電場對電力系統小干擾穩定和阻尼特性的影響[J].電網技術,2007,31(13):76-80.
[9] 王毅,張祥宇,李和明,等.永磁直驅風電機組對系統功率振蕩的阻尼控制[J].電工技術學報,2012,27(12):163-165.
[10] KIA S H,HENAO H.Torsional vibration effects on induction machine current and torque signatures in gearbox-based electromechanical system[J].IEEE Transactions on Industrial Electronics,2009,56(11):4689-4699.
[11] 耿華,許德偉,吳斌,等.永磁直驅變速風電系統的控制及穩定性分析[J].中國電機工程報,2009,29(30):68-75.
[12] 郝正航,余貽鑫.勵磁控制引起的雙饋風電機組軸系扭振機理[J].電力系統自動化,2010,34(21):81-86.
[13] 魯玉普,解大,孫俊博,等.風電場機網扭振的小信號建模及仿真[J].電網技術,2016,40(2):3-6.
[14] 蔣長江,劉俊,劉友波,等.計及風電隨機激勵的電力系統暫態穩定分析[J].電力自動化設備,2016,36(3):100-103.
[15] 張友鵬,李萌,高鋒陽,等.直驅型機組風電場并網靜態電壓穩定性研究[J].電源技術,2013,137(5):811-813.
[16] 徐光虎,王杰,陳陳,等.基于微分代數模型的AC/DC系統非線性控制器設計[J].中國電機工程學報,2005,25(7):52-57.
[17] 李嘯驄,程時杰,韋化,等.一種高性能的非線性勵磁控制[J]. 中國電機工程學報,2003,23(12):37-42.
[18] 李嘯驄,程時杰,韋化,等.輸出函數在單輸入單輸出非線性控制系統設計中的重要作用[J].中國電機工程學報,2004,24(10):50-56.
[19] 李嘯驄,程時杰,韋化.輸出函數在多輸入多輸出非線性控制系統設計中的重要作用[J].中國電機工程學報,2006,26(9):87-93.