唐俊

（四川省電力公司資陽公司，四川 資陽 641300）

引言

電力系統(tǒng)狀態(tài)估計是能量管理系統(tǒng)(energy management system，EMS)中的重要組成部分，其結(jié)果直接影響電網(wǎng)調(diào)度的智能化分析與決策。狀態(tài)估計是利用實時量測系統(tǒng)的冗余度來提高數(shù)據(jù)精度，自動排除隨機(jī)干擾所引起的錯誤信息，估計或預(yù)報系統(tǒng)的狀態(tài)。

1969年，F(xiàn).c.schweppe等人針對電力系統(tǒng)安全監(jiān)控問題把靜態(tài)狀態(tài)估計理論引入電力系統(tǒng)；同期，A.Debs等人提出了Kalman濾波算法，并建立了動態(tài)狀態(tài)估計理論。隨著電力系統(tǒng)的發(fā)展和電網(wǎng)規(guī)模的壯大，對狀態(tài)估計有了新要求，涌現(xiàn)出新的理論和新的技術(shù)，為解決狀態(tài)估計的某些問題提供了可能。本文對電力系統(tǒng)經(jīng)典狀態(tài)估計理論的發(fā)展和狀態(tài)估計新方法進(jìn)行了綜合闡述。

1 靜態(tài)狀態(tài)估計

靜態(tài)狀態(tài)估計不考慮系統(tǒng)狀態(tài)隨時間的變化，僅利用某一時間斷面的量測信息估計電網(wǎng)的狀態(tài)。靜態(tài)狀態(tài)估計算法主要包括加權(quán)最小二乘估計方法和抗差估計方法。

1.1 加權(quán)最小二乘估計

電力系統(tǒng)狀態(tài)估計器采用的估計準(zhǔn)則大多是極大似然估計，其中加權(quán)最小二乘(Weight Least Square，WLS)估計是當(dāng)前普遍采用的經(jīng)典法方程(Normal Equations，NE)估計方法。WLS估計器的優(yōu)點是模型簡單，計算量小，收斂性好，對理想正態(tài)分布的量測量，估計具有一致且無偏等優(yōu)良統(tǒng)計特性。然而，由于該算法的計算量和使用內(nèi)存量較大，難以應(yīng)用于大型電力系統(tǒng)的實時計算；而且在權(quán)重因子相差較大、節(jié)點注入型量測較多或長線路和短線路相連等情況下可能出現(xiàn)病態(tài)，引起數(shù)值計算穩(wěn)定性問題。為了克服法方程法的缺陷，出現(xiàn)了各種改進(jìn)算法。

結(jié)合電力系統(tǒng)的物理特性，吸取潮流計算經(jīng)驗而建立了快速解耦狀態(tài)估計算法，對有功和無功進(jìn)行分解計算，同時使雅可比矩陣常數(shù)化，兼顧了計算速度、收斂穩(wěn)定性、內(nèi)存使用量和對各種類型量測量的適應(yīng)性等方面的優(yōu)點。量測變換狀態(tài)估計算法(又稱唯支路法)僅采用支路潮流量測，計算速度快，使用內(nèi)存少，編程簡單，但該法難以處理節(jié)點注入型量測。

為提高法方程法的數(shù)值穩(wěn)定性，主要的改進(jìn)算法有正交變換法、Peter and wilkinson法、混合法、帶等式約束的法方程 (NE/C)法以及Hachtel增廣矩陣法等。正交變換法對雅可比矩陣進(jìn)行QR正交變換以提高迭代表達(dá)式的抗病態(tài)條件能力，具有良好的數(shù)值穩(wěn)定性，但計算量和內(nèi)存使用量大，且無法進(jìn)行快速解耦計算；Peter and wilkinson法則對雅可比矩陣進(jìn)行LU分解，在數(shù)值穩(wěn)定性和矩陣稀疏方面是比較好的折中；混合法是法方程法與正交變換法的混合算法，該法首先對雅可比矩陣進(jìn)行正交變換，然后用前推、回代來求解法方程；NE/C法把零注入功率虛擬量測作為等式約束引入估計中，解決虛擬量測權(quán)重過大問題，在數(shù)值穩(wěn)定性和計算效率方面是比較好的折中；Hachtel法將殘差作為獨立變量進(jìn)行處理，把殘差方程增廣到等式約束法中，兼顧了計算效率和數(shù)值穩(wěn)定性，但NE/C法和Hachtel法的迭代方程系數(shù)矩陣不再是正定矩陣。在實際應(yīng)用中，應(yīng)根據(jù)系統(tǒng)規(guī)模和量測配置情況來選擇適當(dāng)?shù)姆椒ā?/p>

1.2 抗差估計

量測誤差服從正態(tài)分布是WLS估計獲得最優(yōu)、一致、無偏估計結(jié)果的前提，然而實際量測并不完全服從正態(tài)分布，因此許多學(xué)者研究抗差估計器。所謂抗差估計是指在粗差不可避免的情況下，選擇適當(dāng)?shù)墓烙嫹椒ㄊ勾植顚烙嬛档挠绊懽钚??？共罟烙嬈鞔笾路譃槿悾篗估計器、L估計器、R估計器。M估計是經(jīng)典的極大似然估計的推廣，通常稱為廣義極大似然估計；L估計是基于排序統(tǒng)計量的線性組合估計，計算簡單，適用于計算機(jī)模擬不同類型的抗差估計并對其抗差能力和效率進(jìn)行檢驗；R估計是通過秩檢驗導(dǎo)出估計量，又稱秩檢驗估計。

三類抗差估計中，M估計應(yīng)用最多，它最接近經(jīng)典的WLS估計。M估計最常用的方法是加權(quán)最小絕對值估計(weighed least absolute value，WLAV)。WLAV估計一般采用線性規(guī)劃方法來求解，計算量大，易受杠桿點的影響。抗差估計還包括非二次準(zhǔn)則估計、最小平方中值法估計和最小修剪平方法估計，二次準(zhǔn)則估計器假設(shè)量測服從Huber分布，最小平方中值法估計和最小平方中值法估計假設(shè)量測服從拉普拉斯分布。

抗差估計的顯著優(yōu)點是能夠較真實地反映量測的實際分布模式，但其具有模型復(fù)雜和計算量大等缺點，并未得到廣泛的應(yīng)用。

2 動態(tài)狀態(tài)估計

動態(tài)狀態(tài)估計是相對靜態(tài)狀態(tài)估計而言的。靜態(tài)狀態(tài)估計僅根據(jù)某時刻的量測數(shù)據(jù)確定該時刻狀態(tài)量；動態(tài)狀態(tài)估計則是按運動方程與某一時刻的量測數(shù)據(jù)作為初值進(jìn)行下一時刻狀態(tài)量的估計，兼有狀態(tài)估計和狀態(tài)預(yù)報的功能。動態(tài)狀態(tài)估計具有預(yù)報功能，在處理系統(tǒng)異常情況及電力系統(tǒng)安全分析評估、經(jīng)濟(jì)調(diào)度、電壓穩(wěn)定、預(yù)防控制等方面有較長的決策時間，具有重要的實用價值。

動態(tài)狀態(tài)估計主要是基于擴(kuò)展卡爾曼濾波(Extended Kalman Filter，EKF)理論所建立。該理論是經(jīng)由系統(tǒng)上的量測值及狀態(tài)向量的預(yù)測值，輔以加權(quán)最小二乘(WLS)法的模式進(jìn)行濾波運算，從而推導(dǎo)出系統(tǒng)的最佳狀態(tài)值。

20世紀(jì)70年代初，DebS等人就提出了實用Kalman濾波算法，并用最簡單的系統(tǒng)模型建立動態(tài)狀態(tài)估計。隨后，Schweppe等人提出了跟蹤狀態(tài)估計器，把前一時刻的狀態(tài)估計值作為下一時刻狀態(tài)估計的初值，是一種沒有考慮任何系統(tǒng)動態(tài)模型但卻跟蹤系統(tǒng)狀態(tài)的有效狀態(tài)估計算法。在考慮負(fù)荷變化的情況下，Nishiya等人在該系統(tǒng)模型中引入傾斜因子。該方法把傾斜因子看成固定不變的，但該因子實際上隨系統(tǒng)運行狀態(tài)變化，即是隨機(jī)改變的。為了克服這一缺點，Silva等人用常參數(shù)指數(shù)平滑法建立系統(tǒng)模型，以提高估計精度。Mallieu等人提出了一種更為實際的動態(tài)狀態(tài)估計模型，該模型基于節(jié)點負(fù)荷預(yù)報，預(yù)報環(huán)節(jié)是在注入節(jié)點上進(jìn)行，而不是傳統(tǒng)的狀態(tài)量處理模型。這種方法的前提是：l)系統(tǒng)動態(tài)改變是由負(fù)荷引起的；2)負(fù)荷之間是相互獨立的；3)節(jié)點負(fù)荷的動態(tài)改變可以通過簡單、準(zhǔn)確的負(fù)荷預(yù)報模型來判斷。隨后，Silva等人提出了一種基于在線處理系統(tǒng)狀態(tài)的動態(tài)模型。這個模型具有快速的動態(tài)反應(yīng)能力和良好的狀態(tài)預(yù)測能力。

3 狀態(tài)估計新方法

3.1 新息圖狀態(tài)估計

新息圖狀態(tài)估計是利用新息圖的方法，把新息和網(wǎng)絡(luò)圖結(jié)合起來進(jìn)行計算和推理，以處理壞數(shù)據(jù)、拓?fù)渥兓拓?fù)荷突變等異常情況，進(jìn)而計算系統(tǒng)的運行狀態(tài)。新息圖狀態(tài)估計的主要特點是建立了新息向量元素之間的空間聯(lián)系。既利用預(yù)報狀態(tài)來檢測不同時刻每個量測自身在時間方面的縱向變化，又運用傳統(tǒng)靜態(tài)估計方法中所利用的不同量測之間的符合電路定律的空間的橫向聯(lián)系。通過縱橫兩方面的結(jié)合，使得新息圖兼顧了靜態(tài)狀態(tài)估計和動態(tài)狀態(tài)估計的優(yōu)點，識別功能大大提高，同時給出了狀態(tài)估計和狀態(tài)預(yù)估的結(jié)果，為電網(wǎng)調(diào)度和安全監(jiān)控提供快速可靠的實時數(shù)據(jù)。

新息圖法和傳統(tǒng)狀態(tài)估計方法相比，識別拓?fù)浣Y(jié)構(gòu)變化能力強(qiáng)，計算速度快，所需冗余度小。但它在建立新息網(wǎng)絡(luò)圖時，采用的是直流潮流模型，即忽略了有功功率的損耗，因而也就不可能考慮線路的充電電容支路，因而新息圖模型并不是一個很精確的模型。

3.2 諧波狀態(tài)估計

電力諧波會引起通訊干擾和電磁污染，諧波的傳播嚴(yán)重影響電網(wǎng)的電能質(zhì)量，對電力系統(tǒng)本身的安全經(jīng)濟(jì)運行也會造成危害。對電力系統(tǒng)的諧波進(jìn)行監(jiān)測和抑制必須首先確定諧波源的分布和狀態(tài)，只有實時地掌握電網(wǎng)中的諧波狀態(tài)才能經(jīng)濟(jì)有效地解決諧波治理問題。

諧波狀態(tài)估計是一項有效的諧波監(jiān)控與分析技術(shù)，是傳統(tǒng)潮流狀態(tài)估計在諧波領(lǐng)域的延伸。它利用有限的諧波測量數(shù)據(jù)，根據(jù)一定的估計準(zhǔn)則來估計整個電網(wǎng)的諧波分布，推斷電網(wǎng)的諧波電壓與諧波電流狀態(tài)，并分析找出系統(tǒng)中的諧波源，達(dá)到對整個系統(tǒng)進(jìn)行諧波監(jiān)測和諧波管理的目的。諧波狀態(tài)估計的根本目標(biāo)在于建立諧波情況下三相電力系統(tǒng)的數(shù)學(xué)模型，進(jìn)而提出一套適合于電力系統(tǒng)諧波狀態(tài)估計的計算方法，根據(jù)有限點上的量測數(shù)據(jù)，通過所選擇的估計器實現(xiàn)對諧波源位置、類型和注入電流大小的識別，從而為電力系統(tǒng)諧波的管理和抑制提供依據(jù)。

3.3 配電網(wǎng)狀態(tài)估計

為了有效地監(jiān)控配電網(wǎng)的安全性和經(jīng)濟(jì)性，并進(jìn)行正確的分析與決策，配電調(diào)度中心需要了解配電網(wǎng)絡(luò)的實時信息，進(jìn)行配電網(wǎng)狀態(tài)估計，進(jìn)而正確、全面地掌握配電系統(tǒng)的運行狀態(tài)。

配電網(wǎng)狀態(tài)估計就是利用配電網(wǎng)量測信息估計出高精度的、完整的、可靠的配電系統(tǒng)狀態(tài)。配電網(wǎng)狀態(tài)估計大都基于量測或狀態(tài)量變換技術(shù)，或者把功率量測轉(zhuǎn)化為電流量測，或者以支路電流、支路功率、負(fù)荷電流為狀態(tài)向量，通過變換使量測函數(shù)線性定常，雅可比矩陣常數(shù)化，狀態(tài)向量實部和虛部解耦，三相解耦，提高計算速度。配電網(wǎng)狀態(tài)估計的結(jié)果用于配電自動化的很多方面，例如短路電流計算、網(wǎng)絡(luò)重構(gòu)、電壓/無功控制等。而且，由于配電網(wǎng)狀態(tài)估計減小了數(shù)據(jù)誤差，修正了數(shù)據(jù)錯誤，提高了數(shù)據(jù)精度，因此還可用于檢查負(fù)荷數(shù)據(jù)，并對負(fù)荷模型和負(fù)荷數(shù)據(jù)加以校正，提高配電系統(tǒng)的安全與經(jīng)濟(jì)運行水平。

4 研究展望

隨著電力系統(tǒng)規(guī)模的不斷擴(kuò)大，各種新理論、新技術(shù)的不斷涌現(xiàn)，無論從理論方面還是從實際應(yīng)用需求方面，狀態(tài)估計算法仍有許多問題需要深入研究。狀態(tài)估計算法在以下方面有重要的研究價值，還需進(jìn)一步研究：

1）WAMS系統(tǒng)和PMU技術(shù)逐漸成熟，研究基于PMU量測的實時狀態(tài)估計算法。2）電網(wǎng)規(guī)模不斷擴(kuò)大，研究計算速度快和數(shù)值穩(wěn)定性好的算法，縮短面向大系統(tǒng)的狀態(tài)估計執(zhí)行周期。3）研究模型簡潔、計算量少、應(yīng)用于實際系統(tǒng)的新型抗差狀態(tài)估計算法。4）研究計算速度快、計算精度高、冗余數(shù)據(jù)少、對不良數(shù)據(jù)免疫性能強(qiáng)的多功能結(jié)合的動態(tài)狀態(tài)估計算法。5）其它新理論、新技術(shù)應(yīng)用于電力系統(tǒng)狀態(tài)估計算法的探討和研究。

[1]于爾鏗.電力系統(tǒng)狀態(tài)估計[M].北京：水利電力出版社，1985.

[2]Schweppe F，Wildes J.Power system statiestate estimation，Partl：exactmode[J].IEEE Trans on power Apparatus and Systems，1970.

[3]李可文，張步涵，曲偉君.一種基于量測變換的快速解藕狀態(tài)估計方法[J].華中理工大學(xué)學(xué)報，2000，28(2).

[4]Holten L，Gjelsvik A，Aam S，et al.Comparison of different methods for state estimation[J].IEEE Trans on Power Systems，1988.

[5]周江文，黃幼才，楊元喜，等[M].抗差最小二乘法.武漢：華中理工大學(xué)出版社，1997.

[6]李碧君，薛禹勝，顧錦汶，等.抗差估計理論及其在電力系統(tǒng)中的應(yīng)用 [J].電力系統(tǒng)自動化，1999.

[7]Mili L，Cheniae M，Viehare N，et al.Robust state estimation based on projeetion statisties[J].IEEE Trans on Power Systems，1996，11(2).

[8]周蘇荃，柳悼.新息圖的智能特征[J].電力系統(tǒng)自動化，2000，24(13).

[9]周蘇荃，柳悼.新息圖法識別多重網(wǎng)絡(luò)結(jié)構(gòu)動態(tài)變化[J].中國電機(jī)工程學(xué)報，2001，21(10).

[10]Baran M，Kelley A.A braneh-current-based state estimation method for distribution systems[J].IEEE Trans on Power Systems，1995，10(l).