| <script
var cpro_id = "u1457042"; <iframe id="iframeu1457042_0" scbm?rdid=1457042&dc=2&di=u1457042&dri=0&dis=0&dai=3&ps=345x362&dcb=BAIDU_SSP_define&dtm=BAIDU_DUP_SETJSONADSLOT&dvi=0.0&dci=-1&dpt=none&tsr=0&tpr=1459711225295&ti=%E5%A4%A7%E5%B9%B2%E6%89%B0%E7%A8%B3%E5%AE%9A%E4%B8%AD%E4%BD%8E%E9%A2%91%E6%8C%AF%E8%8D%A1%E6%A8%A1%E5%BC%8F%E7%9A%84%E4%BD%9C%E7%94%A8%E7%A0%94%E7%A9%B6_%E7%94%B5%E6%B0%94%E8%87%AA%E5%8A%A8%E5%8C%96%E6%8A%80%E6%9C%AF%E7%BD%91&ari=1&dbv=0&drs=1&pcs=645x335&pss=970x346&cfv=0&cpl=22&chi=50&cce=true&cec=gbk&tlm=1402381903<u=http%3A%2F%2Fwww.dqjsw.com.cn%2Fdiangongdianzi%2Fdianlidiangong%2F3075.html&ecd=1&psr=1366x768&par=1366x728&pis=-1x-1&ccd=24&cja=false&cmi=34&col=zh-CN&cdo=-1&tcn=1459711225&qn=d95b10249b0250f8&tt=1459711225266.159.257.258" vspace="0" hspace="0" marginwidth="0" marginheight="0" scrolling="no" style="border:0; vertical-align:bottom;margin:0;" allowtransparency="true" align="center,center" width="200" height="200" frameborder="0"> |
1 引言
現(xiàn)代電力系統(tǒng)中存在的低頻振蕩現(xiàn)象是增幅性低頻振蕩小干擾下系統(tǒng)失穩(wěn)的主要原因之一。弱阻尼低頻振蕩模式是增幅性低頻振蕩發(fā)生的內(nèi)在因素。所以，長期以來人們一直都把小干擾穩(wěn)定分析的重點放在對低頻振蕩模式的研究上, 這個時期，低頻振蕩模式是與大干擾強非線性無關(guān)的。
向量場正則形理論作為分析非線性系統(tǒng)的一個新的有效工具，已被用來研究大干擾下stress系統(tǒng)的動態(tài)特性[1~4]。文[1]論證了2階解在大干擾模式間的非線性相關(guān)作用的有效性，文[2] 在2階解的基礎(chǔ)上研究了模式間非線性相關(guān)作用對控制器性能的影響，文[3]利用模式間的非線性相關(guān)作用提出了確定經(jīng)典電力系統(tǒng)模型臨界切除時間的新方法, 文[4]推出了在諧振與準(zhǔn)諧振條件下的2階解析解，找到了大干擾下易失穩(wěn)的參數(shù)域。向量場正則形理論突破了傳統(tǒng)穩(wěn)定分析的局限，把模式和大干擾下系統(tǒng)的動態(tài)特性聯(lián)在了一起，而在電力系統(tǒng)中應(yīng)用向量場正則形理論的關(guān)鍵是求解出非線性正則變換系數(shù)。
本文提出了數(shù)值求解向量場非線性正則變換系數(shù)的算法（ND算法），該算法簡單、方便、實用、有效，適用于任何復(fù)雜的電力系統(tǒng)，給出了實用的鑒別主導(dǎo)低頻振蕩模式的方法。在此基礎(chǔ)上，通過研究低頻振蕩模式與其它模式以及狀態(tài)變量間的非線性相關(guān)作用，把低頻振蕩模式與系統(tǒng)大干擾穩(wěn)定聯(lián)系在一起，探索了低頻振蕩模式在大干擾穩(wěn)定中充當(dāng)?shù)慕巧妥饔?，從另一個側(cè)面揭示了以往大干擾穩(wěn)定分析中所無法涉及的一些新現(xiàn)象，得到了一些新的觀點和新的見解。
2 向量場的正則變換
移平衡點到原點，對n機系統(tǒng)，消去非發(fā)電機節(jié)點的N維狀態(tài)方程為
式中 x為狀態(tài)向量，Y 為電壓和電流組成的中間變量。
若U為系統(tǒng)的右特征向量陣，取線性變換x=UY, 式(1)變?yōu)榧s當(dāng)形系統(tǒng)
式中 Y為約當(dāng)形變量；J為由系統(tǒng)特征根組成的對角矩陣；Y2(Y)為系統(tǒng)的二階項；Yh為高階項。
向量場正則形理論指出，通過正則變換，約當(dāng)形系統(tǒng)中的高階項可以消掉，系統(tǒng)變?yōu)檎齽t形[5]即線性系統(tǒng)如下：
式中 zj為正則形變量；lj為系統(tǒng)第j個特征根。
式(4)是非線性正則變換系數(shù)或稱為非線性相關(guān)系數(shù)，它代表模式間非線性相關(guān)作用的大小。
為矩陣C j的第K行第l列,
式中 V為規(guī)格化的左特征向量；VT=U-1；H P為原系統(tǒng)海森矩陣H的第P個子陣。
式(3)的解析解很容易寫出，在式(3)解析解的基礎(chǔ)上，再利用上述正則變換矩陣的反變換陣Z=Y-h2(Y)，可得式(2)的2階解析解[5]為
從上式可看出，模式仍是2階解的主要成分，而非線性正則變換系數(shù)則是構(gòu)成2階解的一個基本參數(shù)，它包含著一系列重要的非線性信息，是展現(xiàn)系統(tǒng)非線性特性的源泉 [1,2,4]。要求得2階解，其最重要、最關(guān)鍵的一步是求出非線性正則變換系數(shù)。
3 求解非線性正則變換系數(shù)的ND算法
由于電力系統(tǒng)的狀態(tài)方程是由狀態(tài)變量和中間變量y共同組成的，只有從狀態(tài)方程、非線性網(wǎng)絡(luò)方程和機端電壓方程中消掉y，才能得到僅含狀態(tài)變量的封閉的狀態(tài)方程，進而解析求得系統(tǒng)的海森矩陣，再求得非線性正則變換系數(shù)。但由于非線性網(wǎng)絡(luò)方程的存在，要從機端電壓方程和網(wǎng)絡(luò)方程中消掉中間變量y，得到封閉的不含y 的狀態(tài)方程，是根本辦不到的。按以往線性化方法得到的雅克比矩陣，也只是系統(tǒng)一階偏導(dǎo)數(shù)在平衡點處的值，是無法繼續(xù)求海森矩陣的。即使是走別的解析求導(dǎo)的渠道，也將是一件困難和繁雜的事情。所以，求取系統(tǒng)狀態(tài)方程的高階偏導(dǎo)數(shù)就成了研究大干擾下低頻振蕩模式作用的第一道難題。
本文避開了解析求導(dǎo)，在文[7]、[9]的基礎(chǔ)上，提出了數(shù)值求解狀態(tài)方程海森矩陣的算法(ND算法)，不但成功地解決了第一道難題，同時因為ND算法具有簡單方便，適應(yīng)性強的特點，所以也為其它復(fù)雜系統(tǒng)求取高階偏導(dǎo)數(shù)提供了一個有效的工具。
這里，特別要提及的是，f對狀態(tài)變量的一階偏導(dǎo)數(shù)是用解析方法求得的，以往小干擾穩(wěn)定分析中任何一種解析求f雅可比矩陣A的算法都可用，數(shù)值微商處理的僅僅是2階偏導(dǎo)，這是本文所提用數(shù)值微商求海森矩陣的精髓所在。
數(shù)值微商中另一個重要問題就是增量Dx的選取和誤差分析。在式(8)中，x*是系統(tǒng)的原平衡點，一般情況下它不會是零，所以可以給定數(shù)值微商中第j個狀態(tài)變量在平衡點處的增量Dxj為
則數(shù)值微商的截斷誤差應(yīng)是O(10-8Dxj)=O(10-10)階; 若取Dxj=10-5，則數(shù)值微商的截斷誤差應(yīng)是O(10-10)階, 而舍入誤差是O(u/Dxj)=O(10-9)階。這完全可以滿足實際需要。
階, 而舍入誤差是O(u/ 求得系統(tǒng)的海森矩陣，再求非線性正則變換系數(shù)就只是簡單的矩陣乘法運算的問題了。
4 大干擾下主導(dǎo)低頻振蕩模式的鑒別
小干擾下，弱阻尼低頻振蕩模式是影響系統(tǒng)穩(wěn)定的主要因素。文[1]、[2]、[4]指出，大干擾下，振蕩模式的非線性相關(guān)作用是主導(dǎo)系統(tǒng)動態(tài)特性的主要因素。對一個n機系統(tǒng)，應(yīng)有n-1對低頻振蕩模式，究竟哪一個模式對系統(tǒng)動態(tài)特性的影響更大，可以通過約當(dāng)形系統(tǒng)的一階解和2階解的比較看出。
一階解主要提供了模式和狀態(tài)變量間線性相關(guān)的信息，2階解含有模式和模式、模式和狀態(tài)變量間非線性相關(guān)作用的信息，所以通過一階解和2階解的比較可以看出在大干擾下哪一個低頻振蕩模式將被更強烈的激勵，表現(xiàn)出更強烈的非線性相關(guān)作用[2,4]。整個解的比較是很繁瑣的事，文[2]、[4]提出了一種簡化的鑒別所有主導(dǎo)振蕩模式(包括主導(dǎo)低頻振蕩模式和主導(dǎo)控制模式)的公式為
一般情況下，這一條件是能滿足的。但在系統(tǒng)有多個頻率相近的低頻振蕩模式時，會出現(xiàn)兩項或多項
的模近似相等的情況，再僅用一項
就有可能給主導(dǎo)低頻振蕩模式的鑒別帶來誤差。為此本文對式(11) 作了如下修正：
5 算例分析
以中國電力科學(xué)研究院研制的綜合穩(wěn)定程序中的8機系統(tǒng)為算例，見圖1。所有發(fā)電機都采用3階模型(d, w, Eq)，除1號機外，其它發(fā)電機的勵磁系統(tǒng)也都用3階模型，系統(tǒng)總階數(shù)為45階，求得的低頻振蕩模式列于表1。選定大干擾的形式為在節(jié)點30處發(fā)生的三相瞬時短路，0.165s切除故障。
為說明本文提出的ND算法的有效性，對上述系統(tǒng)計算海森矩陣，再求出非線性正則變換系數(shù)，在此列出其前6個元素：
再用解析方法求得海森矩陣，進而求得非線性正則變換系數(shù)，其對應(yīng)的前6個元素為：
-0.00942869069496-j0.000360394117736
0.00340272976985+j0.00030146649439
兩者各元素的前10位數(shù)字是一樣的，僅后3位數(shù)字有誤差，個別的后4位數(shù)字有誤差，與第3節(jié)中的分析是相同的。用非線性正則變換系數(shù)的其它元素來比較，結(jié)果也是如此，14位數(shù)字中也僅后4位數(shù)字有誤差。對文[7]所示的3機系統(tǒng)，文[8] 所示的單機系統(tǒng)，計算非線性正則變換系數(shù)的結(jié)果也同樣。由此可見本文所提ND算法的有效性。
按本文所提算法算得主導(dǎo)低頻振蕩模式是λ21 和λ22 , 若僅按式(11)計算，主導(dǎo)低頻振蕩模式應(yīng)是λ27 和λ28 。由特征根和狀態(tài)變量的線性相關(guān)因子計算知，這對模式與δ1 線性強相關(guān)，也就是與第7臺發(fā)電機線性強相關(guān)。由特征根和狀態(tài)變量的非線性相關(guān)因子[2,4]計算知，這對模式與w1 非線性強相關(guān)，同樣也是與第7臺發(fā)電機非線性強相關(guān)。同時按文[2]、[4]提出的方法求得模式21與模式22，40(或模式22與模式21，40)間的非線性相互作用最大，這里,
算知，模式40與e'q,6非線性強相關(guān)。也就是說，在當(dāng)前的大干擾下，由于主導(dǎo)低頻振蕩模式21和模式40間強烈的非線性相互作用，使得7號發(fā)電機和6號發(fā)電機間也將發(fā)生強烈的非線性相互作用。換句話說，節(jié)點30處發(fā)生短路故障，7號機受到的影響最大，遠離故障的6號發(fā)電機受到的影響應(yīng)該較小。但由于主導(dǎo)低頻振蕩模式21與模式40間強烈的非線性相關(guān)作用，進而與6號發(fā)電機狀態(tài)變量間強烈的非線性相關(guān)作用的結(jié)果，使得遠離故障的6號發(fā)電機也將受到較大的擾動。
用文[2]中提出的非線性正則變換系數(shù)對系統(tǒng)參數(shù)的靈敏度思想同樣可以說明這一特性。按文[2]中的式(16)，用本文提出的數(shù)值微商算法求取
對所有勵磁系統(tǒng)參數(shù)的靈敏度，得模最大的是相對第8號機勵磁系統(tǒng)放大系數(shù)Kt8 的靈敏度，為1.18∠68.4 ，次之為相對于第6號機勵磁系統(tǒng)放大系數(shù)Kt6 的靈敏度，為1.09∠46.7 。這表明兩者的微增將使隨之加大，非線性增強。或者說，從勵磁系統(tǒng)這個角度看，Kt8 和Kt6 在較大程度上決定了代表的非線性的強弱，也就是模式21、模式22和模式40 間非線性相關(guān)作用的強弱。這種相關(guān)作用將使得第8號機﹑第6號機和模式21，進而第7號機間非線性相關(guān)聯(lián)系緊密。
為驗證上述分析的正確性，使用中國電力科學(xué)研究院研制的綜合穩(wěn)定程序做時域仿真，所得以1號機為參考機的相對功角曲線如圖2所示。由圖2可見，對上述大干擾，7號機和8號機在1s時δ7.1 和δ8.1 已超過180°，失去同步；而遠離故障的6號發(fā)電機在1s時δ6.1 也超過180°，失去同步，而且頭三擺呈增幅狀態(tài)。δ7.1 在4擺衰減振蕩后及δ8.1 在2擺衰減振蕩后又都出現(xiàn)增幅振蕩，這也是由于低頻振蕩模式強烈的非線性相關(guān)作用造成的典型的非線性現(xiàn)象[4]。由此驗證了本文提出的這一新的觀點：低頻振蕩模式強烈的非線性相關(guān)作用是大干擾下影響系統(tǒng)動態(tài)特性的重要因素，極端情況時，主導(dǎo)低頻振蕩模式與其它模式、發(fā)電機狀態(tài)變量間強烈的非線性相關(guān)作用是造成遠離故障發(fā)電機失穩(wěn)的主要原因之一。
6 結(jié)論
（1）向量場正則形理論把大干擾強非線性系統(tǒng)的動態(tài)特性研究與系統(tǒng)內(nèi)部的結(jié)構(gòu)特性聯(lián)系在一起，從另一側(cè)面為分析強非線性下系統(tǒng)的穩(wěn)定性以及動態(tài)特性提供了一個新的有效途徑。
（2）本文提出的數(shù)值求解非線性向量場正則變換系數(shù)的算法(ND算法)簡單、方便、實用、有效，適用于任何復(fù)雜的電力系統(tǒng)，解決了應(yīng)用向量場正則形理論的最基本的問題，也為復(fù)雜系統(tǒng)求解高階偏導(dǎo)數(shù)提供了一個極其有效的工具。
（3）本文提出的鑒別大干擾下主導(dǎo)低頻振蕩模式的方法，定量評價模式非線性相關(guān)作用大小的方法，為探索低頻振蕩模式在大干擾穩(wěn)定中的作用，從另一個側(cè)面揭示了以往大干擾穩(wěn)定分析中所難以解釋的一些現(xiàn)象。
（4）實例計算驗證了結(jié)論（2），（3）的有效性，同時得到了一新的觀點：模式間，特別是主導(dǎo)低頻振蕩模式與其它模式、發(fā)電機狀態(tài)變量間非線性相關(guān)作用是大干擾下影響系統(tǒng)動態(tài)特性和穩(wěn)定性的主要因素。
