基于廣域時(shí)空隨機(jī)響應(yīng)的低頻振蕩模態(tài)辨識(shí)

賈 勇，何正友，廖 凱

（西南交通大學(xué) 電氣工程學(xué)院，四川 成都 610031）

0 引言

隨著互聯(lián)電網(wǎng)規(guī)模和復(fù)雜度不斷增加，電力系統(tǒng)低頻振蕩問題已成為制約輸電聯(lián)絡(luò)線路傳輸能力的瓶頸，并嚴(yán)重威脅系統(tǒng)安全穩(wěn)定運(yùn)行［1-2］。近年來，以相量測量單元PMU（Phasor Measurement Unit）為基礎(chǔ)的廣域測量系統(tǒng)WAMS（Wide Area Measurement System）在系統(tǒng)調(diào)度中心得到了廣泛應(yīng)用，實(shí)現(xiàn)了大電網(wǎng)運(yùn)行狀態(tài)信息的廣域時(shí)空統(tǒng)一測量。系統(tǒng)低頻振蕩監(jiān)測與分析是目前WAMS主要高級應(yīng)用之一［3-4］。系統(tǒng)正常運(yùn)行中負(fù)荷隨機(jī)投切等引起小幅波動(dòng)的隨機(jī)響應(yīng)，能夠反映當(dāng)前系統(tǒng)振蕩模態(tài)特征［5-7］。由于電網(wǎng)客觀存在的拓?fù)溥B接及其電磁作用關(guān)系，每個(gè)PMU量測信息具有自相關(guān)性，且多個(gè)PMU量測信息之間具有直接或間接的關(guān)聯(lián)性［8］。為便于系統(tǒng)運(yùn)行人員及時(shí)全面掌握系統(tǒng)運(yùn)行狀態(tài)，研究利用系統(tǒng)廣域時(shí)空隨機(jī)響應(yīng)準(zhǔn)確辨識(shí)振蕩模態(tài)有著重要意義。

低頻振蕩模態(tài)信息包括振蕩頻率、阻尼比和振蕩振型，通常又將振蕩頻率、阻尼比統(tǒng)稱為模式參數(shù)。模式參數(shù)是振蕩模態(tài)在時(shí)間維度上的特征，其中阻尼比參數(shù)是系統(tǒng)小信號穩(wěn)定性的重要刻畫指標(biāo)；振蕩振型則是振蕩模態(tài)在空間維度上的特征，能夠反映振蕩界面及振蕩在電網(wǎng)中的分布情況［9］?，F(xiàn)有基于廣域測量隨機(jī)響應(yīng)的低頻振蕩分析方法大多對振蕩模式參數(shù)或振型分別展開辨識(shí)研究。在模式參數(shù)辨識(shí)方面，針對系統(tǒng)正常運(yùn)行下隨機(jī)負(fù)荷投切等激勵(lì)所產(chǎn)生的系統(tǒng)隨機(jī)響應(yīng)，集中于利用單一通道信號估計(jì)振蕩頻率和阻尼比，如自回歸AR（Auto-Regressive）模型［10-11］、自回歸滑動(dòng)平均 ARMA（Auto-Regressive and Moving Average）模型［12-15］等辨識(shí)方法。由于某一通道信號不能保證對系統(tǒng)的多個(gè)主導(dǎo)模態(tài)同時(shí)都具有較好的可觀性，因而利用單一信號無法同時(shí)準(zhǔn)確辨識(shí)系統(tǒng)的多個(gè)主導(dǎo)振蕩模態(tài)；此外，對某個(gè)振蕩模態(tài)，各個(gè)通道信號的辨識(shí)結(jié)果之間會(huì)存在一定偏差，如何選取或綜合各通道信號的辨識(shí)結(jié)果也將是一個(gè)難題。在振型辨識(shí)研究方面，與振蕩模式參數(shù)估計(jì)相比文獻(xiàn)則相對較少。振型辨識(shí)需要利用電網(wǎng)多個(gè)監(jiān)測點(diǎn)的同步數(shù)據(jù)，現(xiàn)有文獻(xiàn)中主要有 Green 函數(shù)法［16］、功率譜法［17］、傳遞函數(shù)法［18］及頻率分解法［19］等。文獻(xiàn)［18］推導(dǎo)了兩通道系統(tǒng)輸出間的傳遞函數(shù)與振型的聯(lián)系，通過估計(jì)一組系統(tǒng)輸出相互間的有源自回歸（ARX）模型并建立傳遞函數(shù)，進(jìn)而估計(jì)出振蕩振型。文獻(xiàn)［20］對比分析了利用系統(tǒng)隨機(jī)響應(yīng)的振型辨識(shí)方法，指出功率譜法、頻率分解法是廣義傳遞函數(shù)的特殊形式。由于傳遞函數(shù)法不能直接辨識(shí)振蕩模式參數(shù)，需要結(jié)合其他方法估計(jì)得到系統(tǒng)振蕩模式參數(shù)，然后才能估計(jì)對應(yīng)于該振蕩模式的振型。

為了有效利用大電網(wǎng)廣域時(shí)空數(shù)據(jù)，深度挖掘電網(wǎng)固有的時(shí)空關(guān)聯(lián)特性信息，全面給出系統(tǒng)振蕩模態(tài)特征，通過分析廣域多通道信號的向量自回歸VAR（Vector AutoRegressive）模型與系統(tǒng)振蕩模態(tài)之間的聯(lián)系，提出了基于VAR模型的振蕩模態(tài)辨識(shí)方法。通過新英格蘭系統(tǒng)的蒙特卡羅仿真信號及WECC系統(tǒng)實(shí)測信號的測試，與子空間辨識(shí)方法對比分析，結(jié)果表明所提方法計(jì)算更加簡單，能有效辨識(shí)多個(gè)主導(dǎo)振蕩模態(tài)的模式參數(shù)和振蕩振型。

1 基于廣域時(shí)空隨機(jī)響應(yīng)的振蕩模態(tài)辨識(shí)

1.1 振蕩模態(tài)辨識(shí)理論基礎(chǔ)

在研究系統(tǒng)的小信號穩(wěn)定性時(shí)，通常將系統(tǒng)在某一工作點(diǎn)線性化，從而系統(tǒng)可用線性時(shí)不變的狀態(tài)空間模型表示。系統(tǒng)穩(wěn)態(tài)運(yùn)行時(shí)，假定各個(gè)負(fù)荷的隨機(jī)投切所引起的擾動(dòng)為高斯白噪聲，且無法測量，系統(tǒng)的隨機(jī)響應(yīng)是同步等間隔測量的時(shí)間序列，則系統(tǒng)可簡化為如式（1）所示的離散隨機(jī)狀態(tài)空間模型。

其中，xt∈Rn×1為系統(tǒng)狀態(tài)向量；yt∈Rd×1為系統(tǒng)輸出向量；wt和vt分別為系統(tǒng)噪聲和測量噪聲；A為系統(tǒng)狀態(tài)矩陣；C為輸出矩陣。

式（1）中的狀態(tài)空間模型可等效為一個(gè)多變量ARMA 模型［21］：

其中，Λi（i=1，2，…，p）∈Rd×d和 Qj（j=1，2，…，q）∈Rd×d分別為自回歸系數(shù)矩陣和滑動(dòng)平均系數(shù)矩陣；ek（k=t，t-1，…，t-q）為白噪聲。

利用廣域時(shí)空系統(tǒng)隨機(jī)響應(yīng)辨識(shí)振蕩模態(tài)，系統(tǒng)的振蕩模態(tài)特征僅與自回歸系數(shù)矩陣相關(guān)?？梢宰C明，式（2）也可用一個(gè)高階自回歸模型近似［22］。 由于自回歸模型系數(shù)可通過簡單的最小二乘法計(jì)算，文中考慮利用d通道隨機(jī)響應(yīng)向量序列yt，yt=［y1t，y2t，…，ydt］T∈Rd×1，構(gòu)建 VAR 模型［23］：

估計(jì)得到VAR參數(shù)Λi后，系統(tǒng)振蕩模態(tài)的計(jì)算將非常簡單。首先，利用Λi構(gòu)造系統(tǒng)狀態(tài)矩陣H ∈Rdp×dp：

然后，對狀態(tài)矩陣H進(jìn)行特征值分解，得到離散特征值 λr及其對應(yīng)的右特征向量 φr，φr∈Rdp×1。連續(xù)系統(tǒng)模型的特征值可通過式（5）計(jì)算：

其中，Ts為采樣周期。

最后，系統(tǒng)振蕩模態(tài)sr的頻率及阻尼比為：

輸出向量yt的各通道信號對模態(tài)sr的振型為：

其中，I、zi（i=1，2，…，p-1）分別為 d×d 階的單位矩陣和零矩陣。因此，φr的前d個(gè)元素即為模態(tài)sr的振型 ψr。

1.2 VAR模型參數(shù)的最小二乘估計(jì)

假定有連續(xù)等間隔測量的系統(tǒng)隨機(jī)響應(yīng)向量y1、…、yN，VAR 模型階數(shù)為 p，為方便起見，定義以下向量和矩陣：

則式（3）可表示為如式（9）所示的矩陣形式。

VAR模型參數(shù)矩陣Λ可通過如式（10）所示的最小二乘估計(jì)得到。

直接計(jì)算式（10）的協(xié)方差矩陣ZZT存在條件數(shù)變大的問題，由于矩陣Z的QR分解能保持原問題的條件數(shù)不變，因此本文采用QR分解來實(shí)現(xiàn)VAR模型參數(shù)的最小二乘估計(jì)［24］。利用測量數(shù)據(jù)構(gòu)造矩陣 K，K∈RN×（dp+d）。

對矩陣K進(jìn)行QR分解，即K=QR，得到一正交矩陣 Q∈RN×N和一上三角矩陣 R∈RN×（dp+d）。

其中，R11∈Rdp×dp，R12∈Rdp×d，R22∈Rd×d。

VAR模型參數(shù)矩陣及誤差協(xié)方差矩陣的計(jì)算分別如式（13）和式（14）所示。

1.3 隨機(jī)響應(yīng)的功率譜計(jì)算

對于低頻振蕩監(jiān)測等工程應(yīng)用，系統(tǒng)運(yùn)行人員重點(diǎn)關(guān)注的是系統(tǒng)中的主導(dǎo)模態(tài)?；趶V域多通道隨機(jī)響應(yīng)的VAR模型分析辨識(shí)得到的系統(tǒng)振蕩模態(tài)包含許多虛假模態(tài)，因此，從辨識(shí)的模態(tài)中選擇出系統(tǒng)主導(dǎo)模態(tài)也是一個(gè)重要步驟。由于系統(tǒng)隨機(jī)響應(yīng)的能量主要集中在主導(dǎo)模態(tài)的頻率分量上，可通過計(jì)算系統(tǒng)輸出隨機(jī)響應(yīng)的功率譜，并根據(jù)功率譜的峰值分布來確定主導(dǎo)模態(tài)。

VAR 模型的功率譜可由式（13）、（14）計(jì)算的系數(shù)矩陣Λ及誤差協(xié)方差矩陣Δ計(jì)算得到。

1.4 振蕩模態(tài)的辨識(shí)流程

基于廣域時(shí)空隨機(jī)響應(yīng)的振蕩模態(tài)辨識(shí)示意圖如圖1所示，主要包含以下步驟。

（1）數(shù)據(jù)讀取及預(yù)處理：利用發(fā)電機(jī)母線電壓相角信號辨識(shí)系統(tǒng)的振蕩模態(tài)。在辨識(shí)振蕩模態(tài)前，需要對相角信號進(jìn)行預(yù)處理。采用一階差分濾波及高通濾波去除穩(wěn)態(tài)趨勢，然后以10 Hz的頻率進(jìn)行采樣，最后采用有限脈沖響應(yīng)（FIR）零相移低通濾波器（截止頻率為 2.5 Hz）進(jìn)行濾波。

（2）VAR建模及參數(shù)估計(jì)：通過最終預(yù)測誤差準(zhǔn)則和反復(fù)試驗(yàn)相結(jié)合的方法確定VAR模型階數(shù)，并采用QR分解實(shí)現(xiàn)VAR模型參數(shù)的最小二乘估計(jì)，分別由式（13）、（14）計(jì)算模型系數(shù)矩陣和誤差協(xié)方差矩陣。

（3）振蕩模態(tài)參數(shù)計(jì)算：在估計(jì)得到VAR模型參數(shù)后，通過式（4）—（7）計(jì)算振蕩模態(tài)的頻率、阻尼比及振型。

（4）主導(dǎo)模態(tài)確定：利用式（15）計(jì)算低頻振蕩頻率范圍（0.2～2 Hz）的功率譜，并將 d通道隨機(jī)響應(yīng)的自功率譜疊加，然后尋找幅值占絕對優(yōu)勢的局部極大值所對應(yīng)的頻率值，將其與步驟（3）計(jì)算的模態(tài)頻率比較，確定系統(tǒng)的主導(dǎo)模態(tài)。

圖1 振蕩模態(tài)辨識(shí)示意圖Fig.1 Schematic diagram of oscillation mode identification

2 仿真信號分析

通過新英格蘭系統(tǒng)的蒙特卡羅仿真，測試分析所提模態(tài)辨識(shí)方法的有效性。新英格蘭系統(tǒng)包含10臺(tái)發(fā)電機(jī)、39條母線，接線圖如圖2所示。系統(tǒng)中發(fā)電機(jī)采用4階模型，除發(fā)電機(jī)G1外，其余發(fā)電機(jī)均配備了勵(lì)磁器，負(fù)荷設(shè)置為恒阻抗負(fù)荷。對該測試系統(tǒng)進(jìn)行特征值分析，得到3個(gè)區(qū)間振蕩模態(tài)，各模態(tài)的特征值、模態(tài)頻率、阻尼比及主要參與機(jī)組如表1所示。

圖2 新英格蘭系統(tǒng)單線圖Fig.2 Single-line diagram of New England system

表1 新英格蘭系統(tǒng)的區(qū)間振蕩模態(tài)Table 1 Inter-area oscillation modes of New England system

為了模擬實(shí)際系統(tǒng)中小幅度變化的隨機(jī)響應(yīng)，將系統(tǒng)中各個(gè)負(fù)荷容量的1%設(shè)置為高斯白噪聲，等效實(shí)際電網(wǎng)的負(fù)荷隨機(jī)變化。對新格蘭系統(tǒng)進(jìn)行蒙特卡羅仿真200次，每次仿真時(shí)長為10 min，收集各發(fā)電機(jī)母線的電壓相角信號，然后疊加20 dB高斯白噪聲來模擬測量噪聲的影響。

通過特征值分析可知，3個(gè)區(qū)間振蕩模態(tài)的主要參與機(jī)組為 G1、G2、G3、G5和 G9，因此選取這 5 臺(tái)發(fā)電機(jī)的母線電壓相角信號來進(jìn)行系統(tǒng)振蕩模態(tài)辨識(shí)。圖3為母線1的電壓相角預(yù)處理后得到的零均值信號。 對 G1、G2、G3、G5和 G9的母線相角信號分別分析，采用最終預(yù)測誤差準(zhǔn)則確定的模型階數(shù)分別為18、20、20、24和 24。 綜合200次仿真的模式參數(shù)辨識(shí)結(jié)果，計(jì)算振蕩頻率、阻尼比的平均值及標(biāo)準(zhǔn)差，結(jié)果如表2所示。表中“—”表示利用某通道隨機(jī)響應(yīng)未能有效辨識(shí)出該振蕩模態(tài)。

圖3 預(yù)處理后母線1的相角Fig.3 Phase angle of Bus 1 after pre-processing

從表2可以看出，盡管G5與G9是模態(tài) 3的主要參與機(jī)組，在模型階數(shù)p=24時(shí)未能有效辨識(shí)出模態(tài)3。考慮到G5與G9也是模態(tài)2的主要參與機(jī)組，同時(shí)模態(tài)2與模態(tài)3的頻率相近，可能是模態(tài)間的相互影響導(dǎo)致未能從G5和G9的母線電壓相角信號中辨識(shí)出模態(tài)3。多次試驗(yàn)發(fā)現(xiàn)增加自回歸階數(shù)至一定數(shù)值后，如p=100，從G5和G9的母線電壓相角信號中辨識(shí)出了3個(gè)振蕩模態(tài)。限于篇幅，僅給出了G9母線電壓相角信號的估計(jì)結(jié)果，見表2。對比階數(shù)p=24和p=100時(shí)的估計(jì)結(jié)果，盡管增加自回歸階數(shù)能夠估計(jì)出3個(gè)振蕩模態(tài)，但3個(gè)模態(tài)的頻率和阻尼比的估計(jì)結(jié)果的方差大幅增加，估計(jì)精度較差。

表2 單通道信號的模式參數(shù)估計(jì)結(jié)果Table 2 Results of modal parameter estimation based on single-channel signals

從以上單通道隨機(jī)響應(yīng)的測試分析可以看出，利用單通道信號無法準(zhǔn)確辨識(shí)多個(gè)區(qū)間振蕩模態(tài)，且不同通道信號的辨識(shí)結(jié)果存在一定偏差。此外，利用單通道信號不能給出振蕩模態(tài)的振型信息。因此，有必要深度挖掘電網(wǎng)固有的時(shí)空關(guān)聯(lián)特性，利用電網(wǎng)廣域時(shí)空隨機(jī)響應(yīng)進(jìn)行系統(tǒng)模態(tài)辨識(shí)。

利用本文方法對5臺(tái)發(fā)電機(jī)的母線電壓相角信號進(jìn)行VAR建模分析，先采用最終預(yù)測誤差方法來初步確定自回歸階數(shù)，并通過反復(fù)試驗(yàn)法確定的模型階數(shù)為20。5臺(tái)發(fā)電機(jī)的電壓相角信號的疊加功率譜如圖4所示。從圖中可以看出，功率譜中包含3個(gè)局部峰值，分別對應(yīng)表1中的3個(gè)區(qū)間振蕩模態(tài)，因而可以通過功率譜峰值所對應(yīng)的頻率值來確定系統(tǒng)主導(dǎo)模態(tài)。

圖4 電壓相角的疊加功率譜Fig.4 Superimposed power spectrum of voltage phase-angle

對200次仿真的辨識(shí)結(jié)果進(jìn)行統(tǒng)計(jì)分析，并與子空間辨識(shí)方法［25］的辨識(shí)結(jié)果比較，如表3所示。對比表1和表3可以發(fā)現(xiàn)，本文方法與子空間辨識(shí)方法的辨識(shí)結(jié)果非常接近，3個(gè)模態(tài)的頻率、阻尼比的平均值與模態(tài)真實(shí)參數(shù)的偏差較小，且辨識(shí)結(jié)果有較小的標(biāo)準(zhǔn)差，表明2種方法都能有效辨識(shí)系統(tǒng)多個(gè)主導(dǎo)模態(tài)的振蕩頻率和阻尼比。與阻尼比相比，模態(tài)頻率的估計(jì)結(jié)果更為精確。

表3 2種方法的模式參數(shù)估計(jì)結(jié)果對比Table 3 Comparison of estimated modalparameters between two methods

振蕩振型的辨識(shí)有助于直觀地顯示某模態(tài)在電網(wǎng)中的分布情況，能夠?yàn)橄到y(tǒng)同調(diào)機(jī)組識(shí)別及快速失步解列等提供重要信息。為方便振蕩振型辨識(shí)結(jié)果的比較分析，以振型幅值最大的通道為基準(zhǔn)，對幅值進(jìn)行歸一化；同時(shí)對振型相位進(jìn)行移相處理，使得幅值最大的通道的相位為0°。類似地，計(jì)算200次仿真的各通道信號的振型幅值和相位的平均值A(chǔ)、θ及標(biāo)準(zhǔn)差σA、σθ，結(jié)果如表4所示，表4中同時(shí)給出了3個(gè)振蕩模態(tài)振型的特征值分析結(jié)果和子空間辨識(shí)方法的估計(jì)結(jié)果，其中振型幅值的平均值和標(biāo)準(zhǔn)差均為標(biāo)幺值。

從表4的估計(jì)結(jié)果可以發(fā)現(xiàn)，當(dāng)振型相位有較大的標(biāo)準(zhǔn)差時(shí)其幅值都很?。ㄈ缒B(tài)2的G1、G2、G3，模態(tài)3的G1）；此外，由于相隔較近的模態(tài)2和模態(tài)3的相互影響，同時(shí)參與這2個(gè)模態(tài)的機(jī)組G5和G9的振型估計(jì)方差相對較大，但估計(jì)結(jié)果能夠滿足工程應(yīng)用要求。本文方法與子空間辨識(shí)方法的辨識(shí)結(jié)果非常接近，且與特征值分析結(jié)果相符，都能有效估計(jì)出振蕩模態(tài)的振型。表4中振蕩振型歸一化幅值量化了各通道信號對振蕩模態(tài)的相對參與程度，振型相角則顯示了各通道信號的相互搖擺關(guān)系。

以上新英格蘭系統(tǒng)蒙特卡羅仿真信號的測試分析表明，本文方法及子空間辨識(shí)方法都能利用廣域時(shí)空隨機(jī)響應(yīng)有效辨識(shí)系統(tǒng)振蕩模態(tài)，2種方法的辨識(shí)性能非常接近。最后，在2.4 GHz i3處理器、4 G內(nèi)存的電腦上測試了2種方法200次仿真的辨識(shí)效率。本文方法和子空間辨識(shí)方法的辨識(shí)時(shí)間分別是24.41 s和231.96 s，可見，與子空間辨識(shí)方法相比，本文方法計(jì)算更加簡單有效。

3 實(shí)測信號分析

利用北美廣域頻率監(jiān)測系統(tǒng)記錄的一組WECC系統(tǒng)的實(shí)測信號，分析本文方法利用實(shí)際電網(wǎng)廣域隨機(jī)響應(yīng)辨識(shí)系統(tǒng)模態(tài)的適應(yīng)性。北美廣域頻率監(jiān)測系統(tǒng)的測量終端頻率擾動(dòng)記錄儀FDR（Frequency Disturbance Recorder）通過GPS時(shí)鐘同步采集配電網(wǎng)單相電壓信號，計(jì)算頻率和電壓相量［26］。

FDR在WECC系統(tǒng)的分布如圖5所示。WECC系統(tǒng)主要包含2個(gè)區(qū)間振蕩模態(tài)：N-S模態(tài)（約0.24 Hz）和 Alberta模態(tài)（約 0.36 Hz）。N-S 模態(tài)是WECC系統(tǒng)北部電網(wǎng)對南部電網(wǎng)的振蕩，而Alberta模態(tài)是加拿大Alberta地區(qū)、AZ、NW、CO及CA南部對 WECC 中部電網(wǎng)的振蕩［27］。

表4 2種方法的振蕩振型估計(jì)結(jié)果對比Table 4 Comparison of estimated modal shapes between two methods

圖5 FDR分布圖Fig.5 Distribution of FDR

以2012年WECC系統(tǒng)發(fā)生的某一大擾動(dòng)事件為例，選取含有效測量數(shù)據(jù)的6個(gè)FDR，以字母A—F對其編號。讀取大擾動(dòng)事件發(fā)生時(shí)刻前后各60 min的系統(tǒng)響應(yīng)，圖6為FDR A記錄的大擾動(dòng)發(fā)生前后的電壓相角信號，其中0 min為事件發(fā)生時(shí)刻。對擾動(dòng)后包含明顯振蕩的暫態(tài)響應(yīng)信號采用Matrix Pencil方法進(jìn)行模態(tài)辨識(shí)，僅得到N-S振蕩模態(tài)（頻率0.251 Hz，阻尼比6.7%），表明擾動(dòng)事件未能充分激勵(lì)系統(tǒng)的Alberta模態(tài)。由于系統(tǒng)的真實(shí)模態(tài)未知，以下分析將暫態(tài)響應(yīng)信號的模態(tài)辨識(shí)結(jié)果作為評估系統(tǒng)隨機(jī)響應(yīng)的模態(tài)辨識(shí)結(jié)果的參考值。

圖6 WECC系統(tǒng)實(shí)測相角Fig.6 Measured phase-angle of WECC system

圖7 N-S模態(tài)的模式參數(shù)估計(jì)結(jié)果Fig.7 Estimated parameters of N-S mode

采用本文方法對切機(jī)前后截取的6通道隨機(jī)響應(yīng)進(jìn)行模態(tài)辨識(shí)，設(shè)置時(shí)間窗為10 min，滑動(dòng)步長為5 min，得到10組模態(tài)估計(jì)結(jié)果。利用系統(tǒng)廣域時(shí)空隨機(jī)響應(yīng)有效地辨識(shí)出系統(tǒng)的2個(gè)主要區(qū)間振蕩模態(tài)，由于擾動(dòng)事件分析結(jié)果只得到N-S振蕩模態(tài)，限于篇幅，僅給出系統(tǒng)N-S模態(tài)辨識(shí)結(jié)果。圖7為N-S模態(tài)的振蕩頻率和阻尼比的估計(jì)結(jié)果，可以發(fā)現(xiàn)隨機(jī)響應(yīng)的模式參數(shù)估計(jì)值分布在參考值上下，且與參考值的偏差較小。圖8為N-S模態(tài)的振型估計(jì)結(jié)果，從圖8（a）的歸一化振型幅值圖可知，F(xiàn)DR A在N-S模態(tài)的參與度最大；從圖8（b）的振型相角關(guān)系可知，F(xiàn)DR A、B、C 與 FDR D、E、F 之間相差約180°，顯示出WECC系統(tǒng)北部電網(wǎng)和南部電網(wǎng)相互搖擺關(guān)系，這與文獻(xiàn)［27］中的分析結(jié)果是一致的。因此本文方法利用廣域時(shí)空隨機(jī)響應(yīng)，不僅能準(zhǔn)確估計(jì)振蕩頻率及阻尼比，還能有效辨識(shí)出振蕩振型。

圖8 N-S模態(tài)的振型估計(jì)結(jié)果Fig.8 Estimated shapes of N-S mode

4 結(jié)論

為了有效利用系統(tǒng)廣域時(shí)空隨機(jī)響應(yīng)，提出一種基于VAR建模分析的振蕩模態(tài)辨識(shí)方法。通過新英格蘭系統(tǒng)的蒙特卡羅仿真信號和WECC系統(tǒng)實(shí)測信號的測試分析，表明所提方法能準(zhǔn)確估計(jì)多個(gè)主導(dǎo)模態(tài)的模式參數(shù)，同時(shí)能辨識(shí)出相應(yīng)模態(tài)的振型信息。與子空間辨識(shí)方法相比，廣域多通道隨機(jī)響應(yīng)的VAR建模分析方法計(jì)算更加簡單有效。因此，本文方法能充分利用廣域多通道隨機(jī)響應(yīng)的時(shí)空特性，全面給出系統(tǒng)的主要振蕩模態(tài)特征信息，有利于系統(tǒng)運(yùn)行人員準(zhǔn)確掌握系統(tǒng)的運(yùn)行狀態(tài)。

參考文獻(xiàn)：

［1］洪峰，陳金富，段獻(xiàn)忠，等.弱互聯(lián)大區(qū)電網(wǎng)聯(lián)絡(luò)線功率振蕩研究［J］.中國電機(jī)工程學(xué)報(bào)，2011，31（10）：46-51.HONG Feng，CHEN Jinfu，DUAN Xianzhong，et al.Research on tie line power flow oscillation in weak interconnection power grids［J］.Proceedings of the CSEE，2011，31（10）：46-51.

［2］郝思鵬，楚成彪，張仰飛，等.基于多類型受擾軌跡的電力系統(tǒng)低頻振蕩分析［J］.電力自動(dòng)化設(shè)備，2014，34（12）：52-57.HAO Sipeng，CHU Chengbiao，ZHANG Yangfei，et al.Analysis of power system low-frequency oscillation based on multi-type perturbation trajectory［J］.Electric Power Automation Equipment，2014，34（12）：52-57.

［3］段剛，嚴(yán)亞勤，謝曉冬，等.廣域相量測量技術(shù)發(fā)展現(xiàn)狀與展望［J］.電力系統(tǒng)自動(dòng)化，2015，39（1）：73-80.DUAN Gang，YAN Yaqin，XIE Xiaodong，etal.Development status quo and tendency of wide area phasor measuring technology［J］.Automation of Electric Power Systems，2015，39（1）：73-80.

［4］HAUER J，TRUDNOWSKI D，DESTEESE J.A perspective on WAMS analysis tools for tracking of oscillation dynamics［C］∥IEEE PowerSocietyGeneralMeeting.Tampa，F(xiàn)L，USA：IEEE,2007：1-10.

［5］柳勇軍，時(shí)伯年.基于類噪聲數(shù)據(jù)的電力系統(tǒng)低頻振蕩模式辨識(shí)方法及應(yīng)用［J］.南方電網(wǎng)技術(shù)，2013，7（5）：74-77.LIU Yongjun，SHIBonian.Ambientdatabased identification method for the low frequency oscillation modes of power system and its utilization［J］.Southern Power System Technology，2013，7（5）：74-77.

［6］楊德友，蔡國偉，陳家榮.基于動(dòng)態(tài)響應(yīng)的機(jī)電振蕩特征識(shí)別與能控性評估方法［J］.電力自動(dòng)化設(shè)備，2014，34（2）：54-60.YANG Deyou，CAI Guowei，CHAN K W.Characteristics extraction and controllabilityevaluation based on dynamic response of electro-mechanical oscillation［J］.Electric Power Automation Equipment，2014，34（2）：54-60.

［7］王祥超，張鵬，甄威，等.基于自然激勵(lì)技術(shù)和TLS-ESPRIT方法的低頻振蕩模式辨識(shí)［J］.電力系統(tǒng)自動(dòng)化，2015，39（10）：75-80，130.WANG Xiangchao，ZHANG Peng，ZHEN Wei，et al.Identification of low frequency oscillation modes based on NExT and TLSESPRIT algorithm［J］.Automation of Electric Power Systems，2015，39（10）：75-80,130.

［8］劉道偉，張東霞，孫華東，等.時(shí)空大數(shù)據(jù)環(huán)境下的大電網(wǎng)穩(wěn)定態(tài)勢量化評估與自適應(yīng)防控體系構(gòu)建［J］.中國電機(jī)工程學(xué)報(bào)，2015，35（2）：268-276.LIU Daowei，ZHANG Dongxia，SUN Huadong，et al.Construction of stability situation quantitative assessment and adaptive control system for large-scale power grid in the spatio-temporal big data environment［J］.Proceedings of the CSEE，2015，35（2）：268-276.

［9］薛禹勝，郝思鵬，劉俊勇，等.關(guān)于低頻振蕩分析方法的評述［J］.電力系統(tǒng)自動(dòng)化，2009，33（3）：1-8.XUE Yusheng，HAO Sipeng，LIU Junyong，et al.A review of analysis methods for low-frequency oscillations［J］.Automation of Electric Power Systems，2009，33（3）：1-8.

［10］PIERRE J W，TRUDNOWSKI D J，DONNELLY M K.Initial results in electromechanical mode identification from ambient data［J］.IEEE Transactions on Power Systems，1997，12（3）：1245-1251.

［11］ZHOU N，TRUDNOWSKI D J，PIERRE J W，et al.Electromechanicalmode onlineestimation using regularized robust RLS methods［J］.IEEE Transactions on Power Systems，2008，23（4）：1670-1680.

［12］TRUDNOWSKI D J，PIERRE J W，ZHOU N， et al.Performance ofthree mode-meterblock-processing algorithms for automated dynamic stability assessment［J］.IEEE Transactions on Power Systems，2008，23（2）：680-690.

［13］陳剛，段曉，張繼紅，等.基于ARMA模型的低頻振蕩模式在線辨識(shí)技術(shù)研究［J］.電網(wǎng)技術(shù)，2010，34（11）：48-54.CHEN Gang，DUAN Xiao，ZHANG Jihong，et al.A new approach foronline identification oflow frequency oscillation modes based on auto-regressive moving-average model［J］.Power System Technology，2010，34（11）：48-54.

［14］吳超，陸超，韓英鐸，等.Prony方法和ARMA法在低頻振蕩模式辨識(shí)中的適用性比較［J］.電力自動(dòng)化設(shè)備，2010，30（3）：30-34.WU Chao，LU Chao，HAN Yingduo，etal. Comparison of applicability in low frequency oscillation mode identification between Prony and ARMA［J］.Electric Power Automation Equipment，2010，30（3）：30-34.

［15］賈勇，楊健維，孫斌，等.計(jì)及不確定性的電力系統(tǒng)隨機(jī)響應(yīng)模式估計(jì)方法［J］. 中國電機(jī)工程學(xué)報(bào)，2015，35（11）：2735-2742.JIA Yong，YANG Jianwei，SUN Bin，et al.A mode estimation method for power system stochastic response considering the uncertainty［J］.Proceedings of the CSEE，2015，35（11）：2735-2742.

［16］吳超，陸超，韓英鐸，等.基于廣域測量類噪聲信號的節(jié)點(diǎn)間振蕩相位辨識(shí)［J］.清華大學(xué)學(xué)報(bào)（自然科學(xué)版），2010，50（10）：1613-1618.WU Chao，LU Chao，HAN Yingduo，etal. Identification of oscillation phase relationship among nodes based on wide area measured ambient signals［J］.Journal of Tsinghua University（Science and Technology），2010，50（10）：1613-1618.

［17］TRUDNOWSKI D J.Estimating electromechanical mode shape from synchrophasormeasurements ［J］.IEEE Transactions on Power Systems，2008，23（3）：1188-1195.

［18］ZHOU N，HUANG Z，DOSIEK L，et al.Electromechanical mode shape estimation based on transfer function identification using PMU measurements［C］∥Power& EnergySociety General Meeting.Calgary，AB,Canada：IEEE，2009：1-7.

［19］NING Jiawei，PAN Xueping，VENKATASUBRAMANIAN V.Oscillation modal analysis from ambient synchrophasor data using distributed frequency domain optimization［J］.IEEE Transactions on Power Systems，2013,28（2）：1960-1968.

［20］DOSIEK L，ZHOU N，PIERRE J W，et al.Mode shape estimation algorithmsunderambientconditions：a comparative review［J］.IEEE Transactions on Power Systems，2013，28（2）：779-787.

［21］HANNAN E J，DEISTLER M.The statistical theory of linear systems［M］.Philadelphia，PA，USA：SIAM，1988：9-20.

［22］LJUNG L.System identification：theory for the user［M］.Upper Saddle River，NJ，USA：Prentice Hall，1999：81-93.

［23］VU V H，THOMAS M，LAKIS A A，et al.Operational modal analysis by updating autoregressive model［J］.Mechanical Systems and Signal Processing，2011，25（3）：1028-1044.

［24］NEUMAIER A，SCHNEIDER T.Estimation of parameters and eigenmodes of multivariate autoregressive models［J］.ACM Transactions on Mathematical Software（TOMS），2001，27（1）：27-57.

［25］LARSSON M，LAILA D S.Monitoring of inter-area oscillations under ambient conditions using subspace identification［C］∥IEEE Power&Energy Society General Meeting.Calgary，AB，USA：IEEE，2009：1-6.

［26］ZHANG Yingchen，MARKHAM P，XIA Tao，etal.Wide-area frequency monitoring network（FNET） architecture and applications［J］.IEEE Transactions on Smart Grid，2010，1（2）：159-167.

［27］HAUER J F，MITTELSTADT W A，MARTIN K E，et al.Use of the WECC WAMS in wide-area probing tests for validation of system performance and modeling ［J］.IEEE Transactions on Power Systems，2009，24（1）：250-257.