鄭偉業(yè)，吳文傳，張伯明，孫宏斌

?

基于內點法的交直流混聯(lián)系統(tǒng)抗差狀態(tài)估計

鄭偉業(yè)，吳文傳，張伯明，孫宏斌

（電力系統(tǒng)及發(fā)電設備控制和仿真國家重點實驗室(清華大學電機系)，北京 100084）

針對交直流混聯(lián)系統(tǒng)的狀態(tài)估計問題，分析了其數學模型，對其求解方法進行了研究，并開發(fā)了軟件。介紹了交直流混聯(lián)系統(tǒng)的抗差狀態(tài)估計模型，根據該模型的特點進行變量等價轉換，以降低其復雜度。在利用內點法進行求解時，通過對原問題的海森矩陣進行近似，確保其正定且提升了算法的效率。通過IEEE 9節(jié)點、14節(jié)點、30節(jié)點、118節(jié)點、145節(jié)點系統(tǒng)這5個算例，分析了算法效率以及收斂域。結果表明，與牛頓法相比，基于內點法的狀態(tài)估計的收斂域更寬廣，可以應對更惡劣的現場狀況，具備現場應用潛力。

電力系統(tǒng)；交直流混聯(lián)系統(tǒng)；內點法；抗差狀態(tài)估計

0 引言

高壓直流輸電在我國發(fā)展迅速，例如南方電網已建設成世界上最復雜的交直流混聯(lián)電網。交直流的狀態(tài)估計比純交流的問題，要更為復雜，更難求解。從1980年開始，出現了有關交直流互聯(lián)電網狀態(tài)估計的研究。文獻[1]對于20世紀狀態(tài)估計問題的研究狀況作出了系統(tǒng)性的總結。按照求解思路可分為整體求解法和解耦求解法這兩類。

（1）整體求解法

整體求解法把交直流部分統(tǒng)一納入狀態(tài)估計的數學模型中，并且不考慮他們之間的解耦。文獻[2]提出交直流互聯(lián)系統(tǒng)狀態(tài)估計的概念并給出了實際算例，但是，該文并沒有為換流器建模。文獻[3]對量測方程進行泰勒展開，并且只保留至二階，使狀態(tài)估計在迭代過程中的信息矩陣變?yōu)槎ǔｊ?，對交直流系統(tǒng)進行整體建模，然后一起進行狀態(tài)估計。

總體來說，整體求解法收斂性好，而且建模簡單，缺點是求解起來難度較大，計算量大。

（2）解耦求解法

解耦求解法則在數學建模的時候，把原系統(tǒng)分成相互獨立的部分，以減小問題的求解規(guī)模和復雜度，同時又可以沿用原有交流系統(tǒng)的狀態(tài)估計程序。文獻[4]首先提出基于雙端直流系統(tǒng)的有功-無功-直流順序迭代的交直流估計算法。文獻[5]把交直流混聯(lián)系統(tǒng)解耦為直流網絡、交流網絡以及交直流網絡互聯(lián)部分，并對這三部分分別進行計算，求出狀態(tài)估計的次優(yōu)解。文獻[6]采用了主從分裂的算法，考慮了交直流之間的耦合關系，但是，并沒有提供理論上的收斂性證明。文獻[7]利用高斯消去法，把直流部分對交流部分的影響在信息矩陣中消去，而對于直流部分，則把交流的耦合放到量測方程右側，從而實現交直流的解耦。

簡而言之，解耦求解法求解方法簡單，可以部分沿用先前的程序，但是建模技巧性要求高，不恰當的簡化與假設會給結果帶來誤差，并且，在交替迭代的過程中經常會遇到不收斂的問題。

傳統(tǒng)狀態(tài)估計一般采用牛頓法，但是，在現場應用中發(fā)現，在重載等等惡劣情況下，牛頓法存在收斂性問題。牛頓法存在對初值敏感、收斂域窄的不足，以及無法考慮多種復雜不等式約束和等式約束的缺點，內點法以及其變種算法在電力系統(tǒng)中得到了日益廣泛的應用。其中，文獻[14-17]是內點法在電力系統(tǒng)狀態(tài)估計中的應用。

本文把文獻[18]提出的抗差狀態(tài)估計模型推廣到交直流混聯(lián)系統(tǒng)中，并采用內點法對交直流混聯(lián)系統(tǒng)進行整體求解。本文主要工作有以下三點：

（1）提出了變量等價轉換的方法，對原抗差狀態(tài)估計模型進行轉換，降低其復雜度；

（2）開發(fā)了基于內點法的交直流混聯(lián)系統(tǒng)抗差狀態(tài)估計的程序，用算例驗證了其寬廣的收斂域；

（3）提出了近似海森矩陣的方法，算例結果表明，這可以有效提升算法效率。

1 抗差狀態(tài)估計模型及其等價變量轉換

在文獻[18]中提出的抗差狀態(tài)估計模型，可描述成以下形式。

因此需要對原抗差狀態(tài)估計模型進行等價轉換，引入殘差變量

（3）

則式（1）所描述的模型可以轉換為以下模型。

其中：的狀態(tài)量；的定義如式（3）所示。

即可滿足零注入約束。

做了這樣的變換之后，把量測方程放在等式約束中，而在目標函數中只保留殘差變量，那樣雅可比矩陣各元素的表達式被極大地簡化了。

（7）

（8）

海森矩陣各元素的表達式也變得簡潔不少，且大量元素變?yōu)?。

（10）

（12）

（13）

（15）

這樣變換可以大大降低程序編寫難度。

此外，只需要少量的代碼改動，即可實現狀態(tài)估計模型的轉換。因為在式（4）的模型下，不同狀態(tài)估計器之間的區(qū)別只在于目標函數，量測方程的等式約束都是一樣的。例如，實現最常見的加權最小二乘，只需把式（4）的目標函數改為

等式約束部分無需改動。

這時雅可比矩陣的元素只修改式（7）。

海森矩陣的元素也只需要修改式（12）。

（18）

而雅可比矩陣以及海森矩陣其余元素均無需改動。

2 交直流混聯(lián)電網的整體建模

雙端直流輸電系統(tǒng)如圖1所示。

圖1 雙端直流輸電系統(tǒng)

雙端高壓直流輸電系統(tǒng)本質上可以看成是雙端口元件，潮流的參考方向和電流的參考方向均取流入高壓直流系統(tǒng)為正方向。

在此列出直流內部的量測方程（上標m表示量測）為

（20）

（21）

（23）

（24）

交直流邊界處的量測方程為

（26）

式（20）和式（26）所示的無功功率量測方程中含有根號項，在內點法在迭代的過程其并不總在該函數的定義域內，因此會造成程序錯誤。

（28）

這樣，可以避免在迭代過程中引起計算錯誤。與此同時，原來的直流支路無功功率方程則變?yōu)?/p>

3 內點法求解過程以及海森矩陣的近似

采用原始-對偶內點法來對式（4）的模型進行求解，算法流程如下。

（1）為了簡明起見，把式（4）的模型改寫成標準型為

其中：是狀態(tài)量和殘差變量合并后的復合向量，；函數()包含所有等式約束；()包含所有不等式約束。

其中，為函數()的維數。

（3）構建增廣拉格朗日函數為

（4）設置迭代初始值，將原始變量中，狀態(tài)量中的電壓幅值初值設為1，電壓相角初值設為0，將殘差變量的初始值設為0，將對偶變量的初始值設為0，設置迭代次數，設置最大迭代次數，的取值范圍為：10～100，對于零注入母線，設定殘差變量的第個分量。

（6）如果=1，則直接進行步驟（7）。

其中：和為松弛變量；和為對偶變量；為函數()的維數。

（8）按照以下KKT方程組，計算上述增廣拉格朗日函數對原-對偶內點法中各原始變量和對偶變量的一階偏導數為

（35）

（36）

其中：

以保持海森矩陣正定進行迭代。

此時一般t較小，海森矩陣正定性可以得到滿足，沿下降方向尋找最優(yōu)解。

（11）根據步驟（10）的修正步長，更新上述原始變量和對偶變量， 使

（40）

（12）重復步驟（5）。

4 算例分析

（1）內點法和牛頓法收斂域的對比

為了驗證性能，開發(fā)了基于上述方法的程序，程序中調用了基于原始-對偶內點法的開源包IPOPT（內點法優(yōu)化包）來求解交直流混聯(lián)抗差狀態(tài)估計模型。測試系統(tǒng)的硬件環(huán)境為英特爾雙核i5-3210 M CPU 2.50 GHz，8 GB內存，操作系統(tǒng)為Win7 64 bit，開發(fā)環(huán)境為Visual Studio 2008，IPOPT版本為3.5.5。

以IEEE 30節(jié)點標準測試系統(tǒng)為算例，核估計法的窗寬按經驗設為（也可設為別的值），電壓相角初值均選為0，基于同樣的抗差狀態(tài)估計模型，分別對內點法和牛頓法的收斂域作對比。針對相同數學模型，這兩種求解算法，在收斂的情況下可以得到相同的結果。在此處關心的是兩種算法的收斂域的區(qū)別。

從表1可見內點法的收斂域比牛頓法更寬廣，也就意味著內點法更為魯棒，能在質量惡劣的初值下迅速收斂。這是因為內點法引入了罰函數，當狀態(tài)量在可行域之外時，罰函數能把狀態(tài)量壓回可行域內。而牛頓法遇到這種情況，可能就會發(fā)散。

表1 內點法和牛頓法收斂域對比

內點法優(yōu)良的收斂性能是具有現實意義的。在極端的現場情況下，牛頓法的程序有時會不收斂。例如，負荷水平過高或過低時，節(jié)點電壓幅值會偏離標準電壓較遠，但狀態(tài)估計設初值一般采用平啟動，由于牛頓法收斂域較窄，這時可能就不收斂了。而如果采用內點法來計算，就有可能應對這些極端情況。

依然以IEEE 30節(jié)點標準測試系統(tǒng)為算例，核估計法的窗寬不變，調節(jié)負荷水平，采用平啟動，對內點法和牛頓法兩種算法的收斂域做對比，結果如表2所示。

表2 不同負荷水平下內點法和牛頓法收斂域對比

（2）近似海森矩陣的效果

雖然基于內點法的抗差狀態(tài)估計程序具有更寬廣的收斂域，但是在程序調試時發(fā)現，如果一直采用精確的海森矩陣，當系統(tǒng)規(guī)模較大時（例如對于IEEE 145節(jié)點系統(tǒng)），程序的迭代次數過多，以至于計算的效率較低。為了提高效率，使得內點法更加可能滿足現場對狀態(tài)估計的實時性的需求，本文提出了加速算法的一種方法——在迭代前期對海森矩陣采用式（38）的近似。

為了驗證這種方法的效果，以145節(jié)點標準測試系統(tǒng)為算例，對兩種情況作對比：

1）不近似：一直采用精確的海森矩陣進行計算；

2）近似：在迭代前期對海森矩陣作式（38）的近似處理，轉換閾值設為計算精度的2倍。當互補間隙時，就采用精確的海森矩陣。

得到結果如表3和圖2所示。

表3 迭代次數對比

圖2 估計精度對比

可見，在迭代初期對海森矩陣作數值近似，能有效地減少迭代次數，提高了求解效率，同時不影響狀態(tài)估計的求解精度。

下面對IEEE 9節(jié)點系統(tǒng)、14節(jié)點系統(tǒng)、30節(jié)點系統(tǒng)、118節(jié)點系統(tǒng)、145節(jié)點系統(tǒng)共5個標準系統(tǒng)進行測試，窗寬按經驗取，采用平啟動方式，再對近似和不近似的兩種情況進行測試。計算時間如圖3所示。

圖3 不同規(guī)模系統(tǒng)下的計算時間

5個系統(tǒng)的結果均表明，近似海森矩陣是一種加速算法的有效方法。

（3）抗差性能

在此處驗證一下該程序在交直流混聯(lián)系統(tǒng)中的抗差性能。請注意，抗差性能只與數學模型有關，而與算法的實現無關。如果計算精度設定一致，而且結果收斂的話，牛頓法和內點法的計算結果是一樣的。跟牛頓法相比，內點法的優(yōu)點是收斂域更寬，如前文所述。

把IEEE 30節(jié)點標準測試系統(tǒng)的其中一條線路換成直流線路，得到交直流混聯(lián)系統(tǒng)，其結構圖如圖4所示。其中，節(jié)點10、20之間的直流線路參數如表4~表6。

圖4 30節(jié)點交直流混聯(lián)系統(tǒng)結構圖

為了測試，在量測值加上一個服從(0,0.01)正態(tài)分布的隨機誤差，核估計法的窗寬也取為0.01，則當壞數據比例從0逐漸增至10%時，可得估計結果如表7。

表4 30節(jié)點交直流混聯(lián)系統(tǒng)的直流線路參數(A)

表5 30節(jié)點交直流混聯(lián)系統(tǒng)的直流線路參數(B)

表6 30節(jié)點交直流混聯(lián)系統(tǒng)的直流線路參數(C)

表7 30節(jié)點算例計算結果

可見式（4）的模型應用到交直流混聯(lián)系統(tǒng)中，依然有較強的抗差性。

5 結論

傳統(tǒng)的牛頓法存在收斂域窄的問題，為了克服這個缺點，本文從研究對象、狀態(tài)估計模型、求解算法三方面著手，以指數型目標函數為抗差估計模型，對交直流混聯(lián)電網進行整體建模，并基于內點法開發(fā)了一套抗差狀態(tài)估計程序。

本文提出的等價轉換的技巧，簡化了程序編寫的難度。此外，算例表明，采用內點法來做抗差狀態(tài)估計，具有比牛頓法更為寬廣的收斂域，在質量極為惡劣的初值下也能迅速收斂，能應對現場出現的極端情況。

內點法的抗差狀態(tài)估計程序目前存在效率較低的問題，為此，本文提出一種加速算法的方法——近似海森矩陣，并從理論和仿真實驗兩方面說明其有效性。

本文工作還需要從理論上對算法效率進行研究，對Parzen窗寬進行自適應調整，進一步加速算法，以供現場大規(guī)?；ヂ?lián)電網使用。

[1] 李碧君, 薛禹勝, 顧錦汶, 等. 電力系統(tǒng)狀態(tài)估計問題的研究現狀和展望[J]. 電力系統(tǒng)自動化, 1998, 22(11): 53-60.

LI Bi-jun, XUE Yu-sheng, GU Jin-wen, et al. Status quo and prospect of power system state estimation[J]. Automation of Electric Power Systems, 1998, 22(11): 53-60.

[2] GLOVER J D, SHEIKOLESLAMI M. State estimation of interconnected HVDC/AC systems[J]. IEEE Power Engineering Review, 1983, 3(6): 55-56.

[3] HAQUE M Z, KALAM A. Second-order AC-DC power system state estimator[J]. International Journal of Power & Energy Systems, 1997, 17(2): 138-144.

[4] SIRISENA H R, BROWN E P M. Inclusion of HVDC links in AC power-system state estimation[J]. IEE Proceedings on Generation, Transmission and Distribution, 1981, 128(3): 147-154.

[5] LEITE DA SILVA A M, PERROTTA G, PRADA R B, et al. State estimation for integrated multi-terminal DC/AC systems[J]. IEEE Trans on Power Apparatus and Systems, 1985, 104(9): 2349-2355.

[6] 李大志. 南方電網交直流混合狀態(tài)估計的研究[D]. 北京: 清華大學, 2008.

LI Da-zhi. Study on state estimation of the AC/DC hybrid power system in China Southern Power Grid[D]. Beijing: Tsinghua University, 2008.

[7] DING Qi-feng, CHUNG T S, ZHANG Bo-ming. An improved sequential method for AC/MTDC power system state estimation[J]. IEEE Trans on Power Systems, 2001, 16(3): 506-512.

[8] KARMARKAR N. A new polynomial-time algorithm for linear programming[C] // Proceedings of the Sixteenth Annual ACM Symposium on Theory of Computing, 1984: 302-311.

[9] MONTEIRO R D C, ADLER I, RESENDE M G C. A polynomial-time primal-dual affine scaling algorithm for linear and convex quadratic programming and its power series extension[J]. Mathematics of Operations Research, 1990, 15(2): 191-214.

[10]鐘世霞, 袁榮湘. 內點法在電力系統(tǒng)中的應用述評[J]. 高電壓技術, 2005, 31(12): 76-79.

ZHONG Shi-xia, YUAN Rong-xiang. An application of interior point method in optimization of power systems [J]. High Voltage Engineering, 2005, 31(12): 76-79.

[11] 王劍, 劉天琪, 劉學平. 基于優(yōu)先級排序和內點法的機組優(yōu)化組合[J]. 電力系統(tǒng)保護與控制, 2010, 38(13): 55-60.

WANG Jian, LIU Tian-qi, LIU Xue-ping. Unit commitment optimization based on priority list and interior point methods[J]. Power System Protection and Control, 2010, 38(13): 55-60.

[12] 姚煜, 蔡燕春. 離散粒子群與內點法結合的電力系統(tǒng)無功優(yōu)化[J]. 電力系統(tǒng)保護與控制, 2010, 38(3): 48-52.

YAO Yu, CAI Yan-chun. A hybrid strategy based on DPSO and IPM for optimal reactive power flow[J]. Power System Protection and Control, 2010, 38(3): 48-52.

[13] 張鋒, 段余平, 邱軍, 等. 基于粒子群算法與內點算法的無功優(yōu)化研究[J]. 電力系統(tǒng)保護與控制, 2010, 38(13): 11-16.

ZHANG Feng, DUAN Yu-ping, QIU Jun, et al. Research on reactive power flow based on particle swarm optimization and interior point method[J]. Power System Protection and Control, 2010, 38(13): 11-16.

[14] SINGH H, ALVARADO F L. Weighted least absolute value state estimation using interior point methods[J]. IEEE Trans on Power Systems, 1994, 9(3): 1478-1484.

[15] JABR R A. Primal-dual interior-point approach to compute the solution of the state estimation problem[J]. IEE Proceedings on Generation, Transmission and Distribution, 2005, 152(3): 313-320.

[16] 王雅婷, 何光宇, 劉鎧誠, 等. 求解狀態(tài)估計問題的內點半定規(guī)劃法[J]. 電網技術, 2012, 36(10): 209-215.

WANG Ya-ting, HE Guang-yu, LIU Kai-cheng, et al. An semi-definite programming method for state estimation problem[J]. Power System Technology, 2012, 36(10): 209-215.

[17] RAMIREZ-ARREDONDO J M, BAROCIO E E, CHACON O L. The affine-scaling dual algorithm as an alternative to solve the power system state estimation problem[J]. IEEE Power Engineering Review, 1999, 19(6): 50-52.

[18] WU Wen-chuan, GUO Ye, ZHANG Bo-ming, et al. A robust state estimation method based on maximum exponential square (MES)[J]. IET Generation, Transmission & Distribution, 2011, 5(11): 1165-1172.

Robust state estimator for AC/DC hybrid power system based on an interior point method

ZHENG Wei-ye, WU Wen-chuan, ZHANG Bo-ming, SUN Hong-bin

(State Key Lab of Power Systems, Dept of Electrical Engineering, Tsinghua University, Beijing 100084, China)

A robust SE model is introduced for AC/DC hybrid system. First, to simplify its formulation, the SE model is transformed into an equivalent form by introducing equality constraints and new variables. In this model, the hybrid system is formulated as a whole. Then, an interior point method is used to solve this problem. And an approximating technique is proposed to guarantee the Hessian matrix positive definite, which can significantly improve the efficiency of the solution. Five numerical tests with IEEE 9 node system, IEEE 14 node system, IEEE 30 node system, IEEE 118 node system, and IEEE 145 node system show that the proposed algorithm is efficient, robust and reliable. It has the potential to be applied in real application.

This work is supported by National Key Basic Research Program of China (973 Program) (No. 2013CB228203).

power system; AC/DC hybrid system; interior point method; robust state estimation

TM71

A

1674-3415(2014)21-0001-08

2014-02-19；

2014-03-12

鄭偉業(yè)(1990-)，男，博士生，研究方向為電力系統(tǒng)狀態(tài)估計、分布式自律優(yōu)化與控制；E-mail: neo_adonis@139.com

吳文傳(1973-)，男，博士，教授，博士生導師，主要從事控制中心自動化系統(tǒng)的研究和教學工作；

張伯明(1948-)，男，博士，教授，博士生導師，IEEE Fellow，主要從事電力系統(tǒng)運行、分析和控制的研究工作。

國家重點基礎研究發(fā)展計劃(973 計劃)（2013CB228203）