諧波狀態(tài)估計(jì)諧波狀態(tài)估計(jì)技術(shù)

國外對(duì)于諧波狀態(tài)估計(jì)問題研究較早，1989年著名學(xué)者Heydt就提出了諧波狀態(tài)估計(jì)問題，認(rèn)為諧波狀態(tài)估計(jì)是諧波潮流的逆問題，并提出了一種利用最小方差估計(jì)器的諧波源識(shí)別算法。作者利用關(guān)聯(lián)矩陣建立了諧波量測(cè)量與狀態(tài)變量之間的數(shù)學(xué)模型，選用注入視在功率和線路視在功率作量測(cè)量，并將節(jié)點(diǎn)分為非諧波源和可疑諧波源兩種類型，以減少未知狀態(tài)變量的數(shù)目。但是在波形畸變的情況下，無功功率的定義尚未得到統(tǒng)一認(rèn)識(shí)，因此采用視在功率的方法欠缺說服力，但研究開創(chuàng)了諧波狀態(tài)估計(jì)研究的先河，具有重要的意義。

Meliopoulos 和張帆等人的研究成果中將諧波狀態(tài)估計(jì)問題看作為優(yōu)化問題，并給出了一種最小方差估計(jì)算法。

Ma Haili和Girgis在1996年提出了一種應(yīng)用卡爾曼濾波器識(shí)別諧波源的新算法，適用于非平衡三相電力系統(tǒng)中諧波測(cè)量?jī)x表的優(yōu)化配置，以及諧波源位置及其注入電流大小的最優(yōu)動(dòng)態(tài)估計(jì)。以諧波電流為狀態(tài)變量，諧波電壓為量測(cè)量，建立狀態(tài)方程和量測(cè)方程。對(duì)于確定數(shù)目的諧波測(cè)量?jī)x表，通過計(jì)算不同配置條件時(shí)誤差協(xié)方差矩陣的跡，得到諧波測(cè)量?jī)x表的最佳配置方案和諧波注入的最優(yōu)估計(jì)值。

由于電網(wǎng)中非諧波源母線的數(shù)量可能遠(yuǎn)大于諧波源母線數(shù)量，為減少未知狀態(tài)變量的數(shù)目，杜振平和Arrillaga提出了一種電力系統(tǒng)連續(xù)諧波的狀態(tài)估計(jì)算法。利用關(guān)聯(lián)矩陣的概念建立起諧波量測(cè)量與狀態(tài)變量的數(shù)學(xué)模型，并且將系統(tǒng)母線分為非諧波源母線和可能的諧波源母線兩種類型；此外，還將可能的諧波源母線分為測(cè)量母線和未測(cè)母線兩類。采用上述方法可極大減少未知狀態(tài)變量的數(shù)目，從而極大減少計(jì)算工作量，同時(shí)還可使諧波估計(jì)方程由欠定變?yōu)槌?，增加了估?jì)結(jié)果的可信度。

2000 年， S.S.Matair 和Watson 提出將奇異值分解（Singular Value Decomposition，簡(jiǎn)稱SVD）算法用于電力系統(tǒng)諧波狀態(tài)估計(jì)，該算法能夠在系統(tǒng)非完全可觀即部分可觀、估計(jì)方程欠定時(shí)的情況下進(jìn)行有效估計(jì)， 降低了對(duì)測(cè)量冗余的要求。當(dāng)系統(tǒng)完全可觀，估計(jì)方程正定或超定時(shí)，SVD 算法能給出一個(gè)唯一解，并以新西蘭南島220 kV電網(wǎng)為例，分別給出系統(tǒng)完全可觀、部分可觀時(shí)的狀態(tài)估計(jì)結(jié)果，并且與實(shí)際值進(jìn)行對(duì)比，對(duì)比結(jié)果表明奇異值分解法能夠在系統(tǒng)可觀、部分可觀的情況下給出有效估計(jì)值。

選擇節(jié)點(diǎn)電壓作為狀態(tài)量，母線注入電流、母線電壓、支路電流同步量測(cè)作為量測(cè)量進(jìn)行狀態(tài)估計(jì)。對(duì)于有足夠測(cè)量（超定）的方程且測(cè)量方程無病態(tài)時(shí)，通過節(jié)點(diǎn)編號(hào)優(yōu)化，運(yùn)用分層算法對(duì)測(cè)量矩陣進(jìn)行預(yù)處理后再進(jìn)行矩陣求解；對(duì)于測(cè)量方程病態(tài)、欠定時(shí)，采用SVD算法進(jìn)行求解諧波狀態(tài)估計(jì)問題，求得估計(jì)方程的最小二乘解。以IEEE14節(jié)點(diǎn)系統(tǒng)為例，建立系統(tǒng)模型，運(yùn)用MATLAB編程仿真驗(yàn)證了算法的可靠性。而且，還在SVD 算法的基礎(chǔ)上分析了部分可觀系統(tǒng)的測(cè)量問題，進(jìn)而對(duì)測(cè)量配置進(jìn)行了優(yōu)化。

2004 年，吳篤貴、徐政提出了一種基于相量測(cè)量裝置PMU（Phasor Measurement Unit）的狀態(tài)估計(jì)方法。選取節(jié)點(diǎn)電壓相量作為狀態(tài)變量，節(jié)點(diǎn)電壓、支路電流和注入電流相量作為量測(cè)量，采用加權(quán)最小二乘法進(jìn)行狀態(tài)估計(jì)。

上述的諧波狀態(tài)估計(jì)方法都有自己的特點(diǎn)，在某種特定的條件下可在一定程度上實(shí)現(xiàn)諧波狀態(tài)估計(jì)，但也均存在一定的缺點(diǎn)，精度高、速度快與可觀性好的諧波狀態(tài)估計(jì)方法的研究還需進(jìn)一步深化。

早期的諧波狀態(tài)估計(jì)技術(shù)選取功率作為量測(cè)量是受到傳統(tǒng)電力系統(tǒng)狀態(tài)估計(jì)影響的結(jié)果，但因諧波無功功率的定義存在爭(zhēng)議且其測(cè)量裝置沒有統(tǒng)一的標(biāo)準(zhǔn)，由此，采用無功功率的方法依據(jù)不足。隨著基于全球定位系統(tǒng)（GPS）的相量測(cè)量裝置（PMU）的出現(xiàn)以及相應(yīng)的廣域測(cè)量系統(tǒng)（WAMS）在電力系統(tǒng)的逐步建設(shè)，調(diào)度中心可以獲得一種新型的實(shí)時(shí)測(cè)量數(shù)據(jù)———相量量測(cè)，它主要包括節(jié)點(diǎn)電壓的幅值、相角量測(cè)和支路電流的幅值、相角量測(cè)， 有的系統(tǒng)還包括發(fā)電機(jī)節(jié)點(diǎn)的功角量測(cè)。

與SCADA系統(tǒng)相比，相量量測(cè)主要有以下幾個(gè)特點(diǎn)：增加了一種測(cè)量類型，即相角量測(cè)；數(shù)據(jù)測(cè)量精度高，這主要源于測(cè)量裝置的等級(jí)高以及數(shù)據(jù)測(cè)量同步性好，即測(cè)量同步誤差小兩方面；數(shù)據(jù)采集、傳輸速度非?？?。由于技術(shù)的發(fā)展，諧波量的同步測(cè)量成為可能，為此，近年來在諧波狀態(tài)估計(jì)研究中基本采用諧波同步相量作為量測(cè)量，而不采用功率作為量測(cè)量，這樣可使得狀態(tài)估計(jì)方程線性化，計(jì)算量減小并可提高估計(jì)的精度。

目前，諧波狀態(tài)估計(jì)研究主要集中在諧波狀態(tài)估計(jì)建模、諧波狀態(tài)估計(jì)的求解算法、諧波狀態(tài)估計(jì)的可觀性和誤差、PMU 量測(cè)優(yōu)化配置等方面。

隨著大量非線性負(fù)荷的投入使用，越來越多的諧波注入到電網(wǎng)中，使得電網(wǎng)諧波污染引起的電能質(zhì)量問題越來越嚴(yán)重。為有效地治理諧波，必須首先確定電網(wǎng)諧波的分布與狀態(tài)。但在感興趣的所有節(jié)點(diǎn)上都安裝電能質(zhì)量監(jiān)測(cè)裝置是不可能的，為此，需根據(jù)有限節(jié)點(diǎn)上測(cè)量的各次諧波電壓和電流值，通過網(wǎng)絡(luò)分析估算出未知節(jié)點(diǎn)的各次諧波電壓以及各支路的諧波電流，以掌握整個(gè)系統(tǒng)的諧波狀況，為諧波分析和治理提供科學(xué)的依據(jù)，這就是電力系統(tǒng)諧波狀態(tài)估計(jì)問題。

為進(jìn)行諧波狀態(tài)估計(jì)，需按照一定的估計(jì)準(zhǔn)則，對(duì)量測(cè)值進(jìn)行處理，以得到目標(biāo)函數(shù)的最優(yōu)狀態(tài)值。雖然有很多方法均可應(yīng)用于諧波狀態(tài)估計(jì)領(lǐng)域，如最小方差估計(jì)、極大驗(yàn)后估計(jì)、極大似然估計(jì)和最小二乘估計(jì)以及由其衍生出來的加權(quán)最小二乘估計(jì)算法等，然而由于最小二乘法具有所需的先驗(yàn)統(tǒng)計(jì)知識(shí)少、算法簡(jiǎn)單、計(jì)算量小與收斂性好等特點(diǎn)，其在諧波狀態(tài)估計(jì)領(lǐng)域得到了更廣泛的應(yīng)用。

在給定網(wǎng)絡(luò)結(jié)線、支路參數(shù)和量測(cè)系統(tǒng)的條件下，諧波狀態(tài)估計(jì)的量測(cè)方程可寫為：

式中：

z ———量測(cè)量；

x ———狀態(tài)變量，一般是節(jié)點(diǎn)電壓幅值和相位；

v ———量測(cè)誤差。

對(duì)電力系統(tǒng)諧波狀態(tài)估計(jì)問題，引入PMU量測(cè)裝置以后，存在其節(jié)點(diǎn)電壓、支路電流和注入電流相量量測(cè)與狀態(tài)量節(jié)點(diǎn)電壓相量的線性關(guān)系，從而，諧波狀態(tài)估計(jì)數(shù)學(xué)模型大大簡(jiǎn)化。在測(cè)量噪聲可以忽略不計(jì)的情況下，即有如下的線性狀態(tài)方程:

式中：

H_T———量測(cè)矩陣，與系統(tǒng)具體的拓?fù)浣Y(jié)構(gòu)以及量測(cè)點(diǎn)的配置有關(guān)；

X^*———系統(tǒng)狀態(tài)變量；

Z_m———系統(tǒng)量測(cè)量；

W———量測(cè)加權(quán)矩陣。

求解該方程，則得到諧波狀態(tài)問題的最優(yōu)估計(jì)值。

在此基礎(chǔ)上，對(duì)于同一個(gè)系統(tǒng)的狀態(tài)方程，應(yīng)根據(jù)量測(cè)的不同情況采用不同的求解算法。

諧波狀態(tài)估計(jì)的質(zhì)量是量測(cè)數(shù)量和量測(cè)裝置安裝地點(diǎn)的函數(shù)。在實(shí)際電力系統(tǒng)中，具有諧波量測(cè)功能的裝置的配置畢竟是有限的，從經(jīng)濟(jì)角度考慮不可能在感興趣的所有節(jié)點(diǎn)都裝設(shè)。因此，在有限數(shù)量的量測(cè)裝置情況下，為達(dá)到狀態(tài)估計(jì)的目的，需要進(jìn)行量測(cè)點(diǎn)最優(yōu)配置問題的研究，以得到最優(yōu)量測(cè)裝置數(shù)量及量測(cè)點(diǎn)位置。

Heydt 首次提出采用系數(shù)矩陣的最小條件數(shù)法進(jìn)行諧波量測(cè)點(diǎn)配置。針對(duì)靜態(tài)諧波狀態(tài)估計(jì)提出了一種混合非線性最小二乘法，然而，它研究的是選擇一個(gè)最優(yōu)位置放置一個(gè)量測(cè)裝置來確定諧波源的問題，而不是選擇量測(cè)裝置的最優(yōu)數(shù)量和對(duì)諧波數(shù)量的準(zhǔn)確估計(jì)的評(píng)價(jià)。

通過算例分析比較了粒子群算法與遺傳算法、枚舉法在量測(cè)配置方案的計(jì)算收斂速度，對(duì)安裝不同數(shù)量量測(cè)設(shè)備的配置問題，應(yīng)用粒子群算法能夠得到最優(yōu)解，并且收斂速度相對(duì)較快。

a）進(jìn)行諧波狀態(tài)估計(jì)中發(fā)電機(jī)、變壓器、輸電線路以及諧波源等模型建立方面的進(jìn)一步研究。

b）對(duì)諧波狀態(tài)估計(jì)算法進(jìn)一步研究，可以借鑒現(xiàn)有電力系統(tǒng)基波狀態(tài)估計(jì)方法，提出更具普遍和實(shí)用性的新算法。

c）對(duì)量測(cè)點(diǎn)的優(yōu)化配置原則進(jìn)一步研究，在保證系統(tǒng)可觀性的條件下合理布局量測(cè)設(shè)備的安置地點(diǎn)，以更經(jīng)濟(jì)的方式使?fàn)顟B(tài)估計(jì)的結(jié)果質(zhì)量更高。

d）研究與開發(fā)具有同步時(shí)標(biāo)的電能質(zhì)量監(jiān)測(cè)設(shè)備，綜合監(jiān)控電能質(zhì)量各個(gè)指標(biāo)，降低PMU的投資成本。

e）開展諧波狀態(tài)估計(jì)的應(yīng)用研究，在應(yīng)用過程中不斷提高諧波狀態(tài)估計(jì)的精度，使該項(xiàng)研究盡早實(shí)用化。 2100433B

1、傳統(tǒng)狀態(tài)估計(jì)方法

狀態(tài)估計(jì)(state estimation，SE)技術(shù)在電力系統(tǒng)中有較早的應(yīng)用，在諧波應(yīng)用領(lǐng)域中主要用來進(jìn)行諧波源識(shí)別。Heydt 等學(xué)者首先提出了諧波的SE問題并給出了一種利用最小方差估計(jì)器的諧波源識(shí)別算法：選用注入視在功率和線路視在功率作量測(cè)量，利用廣義逆求解欠定方程組，從而獲得對(duì)狀態(tài)變量的最優(yōu)估計(jì)，并以此來獲得負(fù)荷注入系統(tǒng)的諧波功率：當(dāng)注入諧波功率為正時(shí)，則判定該負(fù)荷為諧波源，從而進(jìn)行諧波源識(shí)別。傳統(tǒng)SE方法的估計(jì)結(jié)果是檢測(cè)設(shè)備數(shù)量和安裝位置的函數(shù)，需要較多的諧波檢測(cè)設(shè)備，而冗余的檢測(cè)設(shè)備有時(shí)還是必需的，這使得該方法實(shí)現(xiàn)成本較高。

此外，該方法需要用傅里葉變換(Fourier transform，F(xiàn)T)技術(shù)(如FFT、DFT 等)獲取諧波分量來作為SE 方法的輸入量以估計(jì)諧波水平，但采樣率不足、頻譜泄漏等問題會(huì)影響該方法的準(zhǔn)確性，特別是非正弦非周期信號(hào)、間諧波信號(hào)等影響尤大。這也使得傳統(tǒng)的SE方法無法處理間諧波問題。為改進(jìn)傳統(tǒng)SE 方法的缺陷，后繼的研究人員、學(xué)者提出了不少的改良SE 方法。

2、 結(jié)合最小方差的SE 方法

Meliopoulo 等人將諧波的SE 問題看作優(yōu)化問題，給出了相應(yīng)的最小方差估計(jì)算法。選用諧波電壓作為狀態(tài)變量，電壓和電流作量測(cè)量，使其估計(jì)算法更具普遍性。同時(shí)，描述諧波量測(cè)系統(tǒng)(harmonic measurement system，HMS)的基本結(jié)構(gòu)，并就不對(duì)稱和不平衡狀態(tài)估計(jì)特性的靈敏度分析及HMS 的可觀性分析進(jìn)行了說明。

3、連續(xù)狀態(tài)估計(jì)方法

Du Z P和Arrillag等提出了電力系統(tǒng)連續(xù)諧波的SE算法，以使所裝設(shè)的檢測(cè)儀器盡可能少。該算法不僅利用關(guān)聯(lián)矩陣的概念建立起諧波量測(cè)量與狀態(tài)變量的數(shù)學(xué)模型，而且考慮了電力系統(tǒng)中諧波分布的實(shí)際情況，將母線分為非諧波源母線和可能的諧波源母線(通常非諧波源母線的數(shù)量遠(yuǎn)大于可能的諧波源母線數(shù)量)，這令算法的計(jì)算資源需求大幅減小(未知狀態(tài)變量的數(shù)目大大減少)。通過上述處理，不僅大大減少了算法的計(jì)算量，而且使諧波估計(jì)方程由欠定變?yōu)槌?，估?jì)結(jié)果的可信度也大為增加。因此該算法效率較高，也更經(jīng)濟(jì)實(shí)用。

4、引入小波的間諧波源處理方法

Van 等在Du Z P等學(xué)者的工作基礎(chǔ)上，將SE算法和小波變化技術(shù)結(jié)合起來進(jìn)行間諧波估計(jì)：首先利用小波變化檢測(cè)并獲取間諧波分量，克服FT技術(shù)的缺陷，然后利用符號(hào)觀察分析方法(symbolicobservability analysis，SOA)從一系列的測(cè)量結(jié)果中決定系統(tǒng)是否具有可觀察性，最后再用SE 方法來估計(jì)間諧波的水平，并從估計(jì)結(jié)果中得出間諧波源的類型辨識(shí)和定位結(jié)果。

5、結(jié)合卡爾曼濾波器的方法

Husam等較早將卡爾曼濾波器用于諧波源的狀態(tài)估計(jì)中，但該方法仍無法解決傳統(tǒng)SE 方法需要較多冗余檢測(cè)設(shè)備和依賴FT技術(shù)的問題。

綜上，Heydt 將諧波源識(shí)別問題作為諧波潮流的逆問題利用SE 方法求解，而后繼的工作以減少測(cè)量裝置數(shù)量和脫離對(duì)FT 方法的依賴為目的來改進(jìn)的。

但從諧波功率原理上看，SE 方法仍存在著缺陷：由于負(fù)荷的諧波功率并不能完全表征負(fù)荷的特性，線性負(fù)荷雖然不會(huì)發(fā)出諧波功率，但在復(fù)雜的多諧波源網(wǎng)絡(luò)中，諧波源有可能吸收諧波功率。

因此，這種識(shí)別方法難免會(huì)得出錯(cuò)誤的結(jié)果。并且，在波形存在畸變的情況下，關(guān)于無功功率目前尚無一種能普遍接受的定義。此外，參數(shù)分散性(包括網(wǎng)絡(luò)元件的諧波參數(shù)、網(wǎng)絡(luò)結(jié)構(gòu)的變動(dòng)等)也會(huì)影響SE 算法的性能。

狀態(tài)估計(jì)的數(shù)學(xué)模型是基于反映網(wǎng)絡(luò)結(jié)構(gòu)、線路參數(shù)、狀態(tài)變量和實(shí)時(shí)量測(cè)之間相互關(guān)系的量測(cè)方程:

z=h(x) v

其中z是量測(cè)量;h(x)是狀態(tài)變量,一般是節(jié)點(diǎn)電壓幅值和相位角;v是量測(cè)誤差;它們都是隨機(jī)變量。

狀態(tài)估計(jì)器的估計(jì)準(zhǔn)則是指求解狀態(tài)變量二的原則,電力系統(tǒng)狀態(tài)估計(jì)器采用的估計(jì)準(zhǔn)則大多是極大似然估計(jì),即求解的狀態(tài)變量二`使量測(cè)值z(mì)被觀測(cè)到的可能性最大,用數(shù)學(xué)語言描述,即:

其中f(z)是量測(cè)z概率分布密度函數(shù)。

顯然,具體的目標(biāo)函數(shù)表達(dá)式與量測(cè)z的分布模式密切相關(guān),對(duì)每個(gè)f(幼都有相應(yīng)的極大似然估計(jì)函數(shù)。對(duì)同一系統(tǒng)的相同實(shí)時(shí)量測(cè),若假定的量測(cè)分布模式不同,則得到的估計(jì)結(jié)果不完全相同,因此有不同估計(jì)準(zhǔn)則的估計(jì)器 。2100433B

狀態(tài)估計(jì)是當(dāng)代電力系統(tǒng)能量管理系統(tǒng)(EMS)的重要組成部分,尤其在電力市場(chǎng)環(huán)境中發(fā)揮更重要的作用。狀態(tài)估計(jì)問題的提出激發(fā)了許多學(xué)者的研究興趣,他們以數(shù)學(xué)、控制理論和其它新理論為指導(dǎo),根據(jù)當(dāng)時(shí)的計(jì)算機(jī)軟件和硬件條件,結(jié)合電力系統(tǒng)的特點(diǎn),在理論方面進(jìn)行廣大量研究。同時(shí),以狀態(tài)估計(jì)軟件實(shí)用為目標(biāo),針對(duì)實(shí)際工程面臨的問題,探索和總結(jié)出許多可行的寶貴經(jīng)驗(yàn)。狀態(tài)估計(jì)的理論研究促進(jìn)了工程應(yīng)用,而狀態(tài)估計(jì)軟件的工程應(yīng)用也推動(dòng)了狀態(tài)估計(jì)理論的研究和發(fā)展。迄今為止,這兩方面都取得了大量成果。然而,狀態(tài)估計(jì)領(lǐng)域仍有不少問題未得到妥善解決,隨著電力系統(tǒng)規(guī)模的不斷擴(kuò)大,電力工業(yè)管理體制向市場(chǎng)化邁進(jìn),對(duì)狀態(tài)估計(jì)有了新要求,各種新技術(shù)和新理論不斷涌現(xiàn),為解決狀態(tài)估計(jì)的某些問題提供了可能 。