趙洪山，薛 寧，時 寧

（華北電力大學 新能源電力系統(tǒng)國家重點實驗室，河北 保定 071003）

0 引言

目前，在電力系統(tǒng)計算、仿真、分析與控制過程中，經(jīng)常采用非線性高階動態(tài)模型來提高描述系統(tǒng)動態(tài)行為的能力，但同時帶來了系統(tǒng)模型動態(tài)方程維數(shù)高、模型計算的復雜度（包括計算時間、速度以及存儲空間等）增加等問題，給工程研究人員帶來了巨大的挑戰(zhàn)。尤其對于多機電力系統(tǒng)，隨著系統(tǒng)規(guī)模不斷增大，其動態(tài)模型的階數(shù)也隨之增加，這使得通過計算系統(tǒng)模型來快速完成各種功能變得非常困難[1-3]。

模型降階是解決上述非線性電力系統(tǒng)模型高維動態(tài)行為計算復雜、難以分析等問題的有效方法之一。學者們也已進行了大量研究，如奇異值攝動法[4-7]、解耦算法[8-10]、同調等值理論[11-14]、模態(tài)分析方法[15-16]、基于相關性的模型降階方法[17]以及基于Krylov子空間的模型降階方法[18]。但上述系統(tǒng)模型降階多用于分析研究，較少關注系統(tǒng)的輸入控制，或將輸入控制作恒定處理；而應用于控制設計的系統(tǒng)模型降階，則需要考慮降階模型是否能夠保持原系統(tǒng)的動態(tài)行為和輸入/輸出特性的問題。

在眾多模型降階方法中，Gramian方法已在線性電力系統(tǒng)廣泛應用。文獻[19]為了簡化區(qū)域電網(wǎng)振蕩阻尼控制器的設計，利用Gramian方法研究了電力系統(tǒng)模型降階，并驗證基于降階系統(tǒng)設計的控制器沒有降低原系統(tǒng)閉環(huán)特性的結論；文獻[20]也在抑制低頻振蕩問題上利用Gramian降階方法研究了更大規(guī)模電力系統(tǒng)的模型降階問題。但對于一般非線性動態(tài)電力系統(tǒng)模型降階，由于缺少有效的理論方法，其模型降階仍是研究的難點。本文利用經(jīng)驗Gramian方法和平衡實現(xiàn)理論，探討非線性電力系統(tǒng)動態(tài)模型的降階過程，并以某實際多機非線性電力系統(tǒng)為例，驗證經(jīng)驗Gramian降階方法對非線性電力系統(tǒng)模型降階的有效性。

1 非線性系統(tǒng)模型平衡降階基本原理

1.1 模型降階基本思想

具有輸入和輸出的非線性系統(tǒng)模型一般可以表示為：

其中，f（x，u）描述非線性系統(tǒng)的動態(tài)行為；x（t）?Rn為非線性系統(tǒng)的狀態(tài)變量；u（t）?Rp為非線性系統(tǒng)的控制輸入；y（t）?Rq為非線性系統(tǒng)的輸出。

模型降階的基本思想就是對形如式（1）的非線性動態(tài)系統(tǒng)，通過已知標準列正交轉換矩陣，將原系統(tǒng)變換到一個低維空間Rk（k≤n）上，并可用低維空間中的降階系統(tǒng)來表示原系統(tǒng)，從而降低了原系統(tǒng)理論分析難度和數(shù)據(jù)運算量。假設已知轉換矩陣P，由原始非線性系統(tǒng)（1）可以得到降階系統(tǒng)：

1.2 Galerkin投影

目前，模型降階中最重要和最基本的方法就是Krylov子空間類方法，經(jīng)驗Gramian方法也基于Krylov子空間原理，即通過構建轉換矩陣P來進行模型降階。其中，矩陣P是由Galerkin投影獲得的。

Galerkin投影：一旦得到子空間的特征函數(shù)，就可以利用Galerkin投影將系統(tǒng)（1）投影到該子空間。Galerkin投影矩陣P通?？杀硎緸椋?/p>

其中，Ik×k為 k×k 階的正交矩陣；0k×（n-k）為 k×（n-k）階的全零矩陣。 因此，非線性降階系統(tǒng)（2）中，有 class="picture_character" src="images/2bdd96436f03e63f1ee327b703d8a6f2.png" />（t）=Px（t）。

1.3 經(jīng)驗Gramian矩陣

對于一般非線性系統(tǒng)（1），其輸入與輸出特性通常由系統(tǒng)的可控Gramian矩陣Wc和可觀Gramian矩陣Wo來描述，可控Gramian矩陣描述控制輸入變化對系統(tǒng)狀態(tài)的影響，而可觀Gramian矩陣描述輸出對系統(tǒng)狀態(tài)變化的反應。

對于線性系統(tǒng)，其模型降階可通過求解Lyapunov方程獲得可控Gramian矩陣和可觀Gramian矩陣，然后通過分析Gramian矩陣的特征值，來確定線性系統(tǒng)模型的降階階數(shù)[21]。而對于非線性系統(tǒng)，由于無法獲得Lyapunov方程，目前仍沒有直接計算可控Gramian矩陣和可觀Gramian矩陣的解析方法。Sirovich在文獻[22]中提出：通過仿真或實驗獲得系統(tǒng)狀態(tài)變量和輸出變量樣本，再根據(jù)采樣的經(jīng)驗數(shù)據(jù)樣本，計算得到非線性系統(tǒng)模型的經(jīng)驗可控Gramian矩陣和經(jīng)驗可觀Gramian矩陣。

經(jīng)驗Gramian方法得到的系統(tǒng)數(shù)據(jù)樣本應滿足在其穩(wěn)定域內這一條件，并且要盡可能反映各種可能的擾動行為。為此，定義如下集合：

其中，Tn為一組r個正交的n×n階矩陣，r為激勵或者擾動方向的矩陣數(shù)量；M為一組s個正常數(shù)，s為在每個方向不同擾動大小的數(shù)量；Ep為Rp中標準單位矢量。

定義1設Tn、Ep和M為以上所給定集合，對于非線性系統(tǒng)（1），經(jīng)驗可控Gramian矩陣Wc定義如下：

對于具有輸入和輸出的可控非線性系統(tǒng)，只研究輸入與狀態(tài)的行為是不全面的，同時也應考慮輸出反映系統(tǒng)狀態(tài)的行為。下面定義給出了經(jīng)驗可觀Gramian矩陣的構建方式。

定義2設Tn、Ep和M為以上所給定集合，對于非線性系統(tǒng)（1），經(jīng)驗可觀Gramian矩陣Wo定義如下：

對經(jīng)驗可控和可觀Gramian矩陣分別進行奇異值分解，其非零特征值所對應的特征向量就可生成一個新子空間，然后可以將系統(tǒng)（1）投影到該子空間，得到降階系統(tǒng)（2）。由此，實現(xiàn)了反映系統(tǒng)可觀性與可控性非線性系統(tǒng)的模型降階。

1.4 基于經(jīng)驗Gramian平衡降階

由上文可知，一旦得到非線性系統(tǒng)（1）的可控與可觀Gramian矩陣，便可對其進行奇異值分解，并利用可控或可觀Gramian矩陣的奇異值大小來確定Galerkin投影矩陣P，進而得到降階模型的階數(shù)。但對于一般非線性系統(tǒng)，直接得到的可控或可觀Gramian矩陣往往是不相同的，如果以可控Gramian矩陣的奇異值進行降階，降階系統(tǒng)的可控性獲得較好的性能，但系統(tǒng)的可觀性稍差；反之，以可觀Gramian矩陣的奇異值進行降階，降階系統(tǒng)可得到較好的可觀性，可控性能稍差。

為了使降階系統(tǒng)的可控性和可觀性都能獲得較好的性能，利用平衡實現(xiàn)理論對非線性系統(tǒng)（1）進行處理，將經(jīng)驗可控與可觀Gramian矩陣變換成平衡的可控和可觀Gramian矩陣，即將系統(tǒng)（1）變換成平衡系統(tǒng)。所謂平衡系統(tǒng)，是指系統(tǒng)的可控Gramian矩陣與可觀Gramian矩陣是相等的，并且是對角的。下面給出了利用經(jīng)驗可控Gramian矩陣Wc和可觀Gramian矩陣Wo計算系統(tǒng)（1）的平衡系統(tǒng)的步驟。

a.對經(jīng)驗Wc和Wo，進行Cholesky因式分解，得到矩陣X和Y。

b.構建積矩陣YTX，并對積矩陣進行奇異值分解，得到對角矩陣∑、正交矩陣U和V。

其中，∑ ?Rn×n為 Hankel奇異值矩陣。

c.根據(jù)式（8）可獲得平衡轉換矩陣T。

至此，利用該平衡變換矩陣T就可獲得系統(tǒng)（1）的平衡系統(tǒng)：

平衡系統(tǒng)與其原系統(tǒng)的動態(tài)行為是完全相同的，只是平衡系統(tǒng)的狀態(tài)變量變?yōu)橐唤M新的狀態(tài)變量。從線性系統(tǒng)理論而言，平衡系統(tǒng)任一狀態(tài)都是原系統(tǒng)所有狀態(tài)的線性組合，即。 因此，對非線性系統(tǒng)（1）進行模型降階，若可控和可觀Gramian矩陣不等，就可先把系統(tǒng)轉換成平衡系統(tǒng)，再對平衡系統(tǒng)進行模型降階，該過程稱為平衡實現(xiàn)。

針對非線性系統(tǒng)對應的平衡系統(tǒng)，計算其經(jīng)驗Gramian矩陣，并對經(jīng)驗 Gramian矩陣和進行奇異值分解。根據(jù)Hankel矩陣∑奇異值σ的大小，當滿足σk?σk+1或者當前k個奇異值能量所占所有奇異值能量絕大部分時，就認為k后面的奇異值所對應的狀態(tài)變量對系統(tǒng)的輸入影響和輸出影響較小，將其略去以實現(xiàn)模型降階的目的。因此，可以將分為主要和非主要兩部分，如下：

2 非線性電力系統(tǒng)經(jīng)驗Gramian模型降階

2.1 非線性電力系統(tǒng)模型

假設多機電力系統(tǒng)有Z臺發(fā)電機，同步發(fā)電機采用計及轉子超瞬變過程且轉子q軸要考慮阻尼繞組的6階動態(tài)模型。第i臺發(fā)電機模型描述如下：

其中，i=1，2，…，Z；δi為第 i臺發(fā)電機轉子角；ωi為第i臺發(fā)電機角頻率；E′di、E′qi、E″di、E″qi分別為第i臺發(fā)電機d軸和q軸暫態(tài)電勢和次暫態(tài)電勢；ωs為發(fā)電機額定角頻率；ωref為發(fā)電機參考角頻率；Pmi為第i臺發(fā)電機的原動機輸出機械功率；Pei為第i臺發(fā)電機的電磁功率；Hi為第i臺發(fā)電機機組慣性時間常數(shù)；T′d0i、T′q0i、T″d0i、T″q0i分別為第i臺發(fā)電機d軸和q軸開路暫態(tài)、次暫態(tài)時間常數(shù)；Efi為第i臺發(fā)電機機組的勵磁電壓；Di為第i臺發(fā)電機定常阻尼系數(shù)；rai、xdi、xqi、x′di、x′qi、x″di、x″qi分別為第i臺發(fā)電機定子電阻和d 軸和 q 軸暫同步、瞬變、超瞬變電抗；Idi、Iqi、Udi、Uqi分別為第i臺發(fā)電機d軸和q軸暫態(tài)電流和機端電壓，其關系為Udi=E″di-raiIdi-x″qiIqi，Uqi=E″qi+x″diIdi-raiIqi。

發(fā)電機與網(wǎng)絡的接口方程為：

其中，Ixy=CSIdq；Uxy=CUdq；C 為 dq－xy坐標系的變化矩陣；S為網(wǎng)絡容量和發(fā)電機容量之間的變化矩陣；Ys為電力系統(tǒng)節(jié)點導納矩陣。

將xy坐標下的Ixy和Uxy轉換到dq坐標系下，代入式（12），消去Idi和Iqi，得到如下形式的具有輸入和輸出多機非線性電力系統(tǒng)模型：

其中，u（t）=[Pm1… PmZEf1… EfZ]T?Rp；x（t）=[δ1…δZω1…ωZE′q1…E′qZE′d1…E′dZE″q1…E″qZE″d1…E″dZ]T?Rn；f（x，u）描述電力系統(tǒng)動態(tài)行為；y（x）為系統(tǒng)的輸出函數(shù)，由關于狀態(tài)x（t）的非線性函數(shù)來表示。

其中，Ui為第 i臺發(fā)電機機端電壓；y（t）?Rq。

2.2 非線性電力系統(tǒng)經(jīng)驗Gramian矩陣

對于非線性電力系統(tǒng)（13），其經(jīng)驗可控Gramian矩陣和可觀Gramian矩陣的計算，由于涉及積分計算，計算量非常龐大。假設系統(tǒng)的樣本是在離散時間t1、…、tm得到的，那么狀態(tài)或輸出的相關矩陣可以寫成如下離散形式：

其中，xk（t）為系統(tǒng)的樣本值；xss為系統(tǒng)狀態(tài)的穩(wěn)態(tài)值。

下面分別給出了經(jīng)驗可控和可觀Gramian矩陣的離散定義。

設Tn、Ep和M為前文所給定的集合，p為系統(tǒng)的控制輸入個數(shù)，則系統(tǒng)經(jīng)驗可控Gramian矩陣Wc定義為：

設Tn、Ep和M為前文所給定的集合，則系統(tǒng)經(jīng)驗可觀Gramian矩陣Wo定義為：

利用經(jīng)驗樣本數(shù)據(jù)，通過式（15）和（16）便可計算出非線性電力系統(tǒng)（13）的經(jīng)驗可控Gramian矩陣Wc和經(jīng)驗可觀Gramian矩陣Wo，其計算精度取決于樣本個數(shù)，樣本數(shù)量越大，精度越高；此外，還受外加激勵大小的選擇cj、外加激勵方向的選擇Tlei以及外加激勵的個數(shù)p等因素的影響。

2.3 基于經(jīng)驗Gramian的非線性電力系統(tǒng)平衡降階步驟

對于非線性電力系統(tǒng)動態(tài)模型（13），利用經(jīng)驗Gramian進行平衡降階方法的步驟如下。

a.根據(jù)定義（15）和（16），計算非線性多機電力系統(tǒng)（13）的經(jīng)驗可控Gramian矩陣Wc和經(jīng)驗可觀Gramian 矩陣 Wo。

b.利用式（6）—（8）計算平衡變換矩陣 T，得到非線性多機電力系統(tǒng)（13）的平衡系統(tǒng)模型。

c.利用平衡變換矩陣T，計算平衡系統(tǒng)的經(jīng)驗可控Gramian矩陣和經(jīng)驗可觀Gramian矩陣：

e.確定非線性多機電力系統(tǒng)降階模型的階數(shù)。如果滿足

則平衡系統(tǒng)的降階模型的階數(shù)為r。其中，ε為接近于1的值。利用階數(shù)r構造Galerkin投影矩陣P。

f.利用平衡變換矩陣T和Galerkin投影矩陣P，計算得到形如式（10）的電力系統(tǒng)平衡降階模型。

3 算例分析

對某實際20機電力系統(tǒng)模型進行降階仿真分析，系統(tǒng)接線圖見圖1。

首先，形成多機電力系統(tǒng)分析模型，每臺同步發(fā)電機采用6階動態(tài)模型，負荷為恒阻抗模型。由于發(fā)電機模型采用dq坐標，網(wǎng)絡模型采用xy坐標，因此，需要xy-dq坐標變換實現(xiàn)發(fā)電機與網(wǎng)絡的互聯(lián)，其過程見式（11）—（13），最后形成一個 120階的形如式（13）的非線性電力系統(tǒng)動態(tài)模型。

圖1 某實際20機系統(tǒng)接線圖Fig.1 Connection diagram of an 20-generator power system

圖1 某實際20機系統(tǒng)接線圖Fig.1 Connection diagram of an 20-generator power system

其次，利用經(jīng)驗Gramian平衡降階算法，計算上述非線性電力系統(tǒng)的平衡變換矩陣T，得到該系統(tǒng)的平衡系統(tǒng)。對平衡系統(tǒng)的經(jīng)驗Gramian矩陣進行奇異值分解，利用奇異值分解得到Hankel奇異值矩陣∑，奇異值大小分布如圖2所示。

圖2 Hankel奇異值的分布圖Fig.2 Distribution of Hankel singular values

最后，利用Hankel奇異值大小，根據(jù)模型降階階數(shù)判定條件確定降階系統(tǒng)的維數(shù)，得到Galerkin投影矩陣P。由此，利用平衡變換矩陣T和投影矩陣P，得到20機非線性電力系統(tǒng)的降階模型。

從Hankel奇異值角度分析，當ε取0.99，可以將原120階非線性電力系統(tǒng)投影到50維的子空間中；當ε取0.95，可將原系統(tǒng)降為32階的非線性動態(tài)系統(tǒng)。理論上，ε取值越小，降階系統(tǒng)的階數(shù)就越小，但投影到子空間后的動態(tài)行為失真就越大。因此，需要合理選擇降階模型的階數(shù)。

對于該20機實際電力系統(tǒng)，每臺發(fā)電機均考慮調速控制和勵磁控制功能，這些控制器的參數(shù)在各降階系統(tǒng)中都保持不變。在考慮三相短路故障（故障點在線路33－34之間50%處）情況下，對各降階模型的暫態(tài)動態(tài)變化過程進行仿真分析。為了節(jié)省篇幅，以發(fā)電機7為例，給出了其相應的仿真曲線如圖3—5所示。圖中，功角δ、角速度ω、輸出電壓Ut均為標幺值。

圖3 降階系統(tǒng)與原系統(tǒng)功角仿真曲線Fig.3 Simulative angle curves of reduced system and original system

圖4 降階系統(tǒng)與原系統(tǒng)角速度仿真曲線Fig.4 Simulative angular speed curves of reduced system and original system

圖5 降階系統(tǒng)與原系統(tǒng)輸出電壓仿真曲線Fig.5 Simulative output voltage curves of reduced system and original system

從仿真曲線圖3—5中可以看出：對于50階的降階系統(tǒng)，除功角幅值誤差稍大些外，角速率和電壓幅值與原系統(tǒng)輸出相比誤差都很小，而且3個輸出響應的頻率特性（即暫態(tài)過程）與原系統(tǒng)的幾乎相同；當系統(tǒng)降到49階時，功角、角速度和電壓幅值3個輸出響應的幅值變化不大，但較50階的降階系統(tǒng)，其與原系統(tǒng)的誤差要比50階系統(tǒng)大，并有逐步增大的趨勢。尤其，當系統(tǒng)降到47階或低于47階時，降階系統(tǒng)輸出響應的幅值誤差增大，觀察功角和角速度曲線出現(xiàn)頻率偏移，即產(chǎn)生了頻率失真。從Hankel奇異值角度分析，投影到50維子空間的能量約為原系統(tǒng)能量的99.45%，其能保留原系統(tǒng)的大部分能量，即保留了原系統(tǒng)的各種特性和動態(tài)行為。

通過對所有節(jié)點仿真曲線進行統(tǒng)計分析可知：當ε選擇大于0.99時，降階模型的誤差都比較小，能夠將原系統(tǒng)的輸入和輸出動態(tài)行為保留下來。為了使降階系統(tǒng)得到較為精確的幅頻特性，通常選擇ε為0.99。本例中，降階系統(tǒng)模型的階數(shù)可降到50階，其所有節(jié)點狀態(tài)和輸出與原系統(tǒng)相比，誤差都非常小。

4 結論

本文提出了一種針對多機非線性電力系統(tǒng)模型的經(jīng)驗Gramian平衡降階方法，并通過某實際20機非線性多機電力系統(tǒng)模型進行仿真驗證。仿真結果表明：降階系統(tǒng)能夠很好地保留原非線性電力系統(tǒng)輸入和輸出的動態(tài)行為以及原非線性電力系統(tǒng)的穩(wěn)態(tài)值。對于該20機非線性電力系統(tǒng)，可以將120階非線性動態(tài)模型降到50階，仿真結果證實了所提方法在非線性電力系統(tǒng)模型降階中的有效性。