程浩忠,顧穎中,熊 寧

(上海交通大學(xué)電氣工程系,上海 200240)

1 引言

無功規(guī)劃與傳統(tǒng)的無功優(yōu)化不同。無功優(yōu)化RPO(Reactive Pow erOptim ization)是在系統(tǒng)現(xiàn)有的無功補償元件及調(diào)節(jié)設(shè)備范圍內(nèi),在滿足系統(tǒng)各項約束條件下,使系統(tǒng)運行費用(通常是使網(wǎng)絡(luò)損耗)最小。而無功規(guī)劃RPP(Reactive Pow er Programming)是為達到某種目標(通常是使系統(tǒng)靜態(tài)電壓穩(wěn)定裕度高于某個閥值)而對需投入無功備用設(shè)備進行規(guī)劃的問題。

常規(guī)無功規(guī)劃問題的求解方法,主要有試驗誤差法、非線性規(guī)劃法、動態(tài)規(guī)劃法、智能進化算法以及線性規(guī)劃法[1,2]。其中,內(nèi)點法以其計算效率的高效性和對約束處理的方便性得到了廣泛的應(yīng)用。Karm arkar在1984年提出內(nèi)點法后[3],許多學(xué)者對其進行了擴展研究,相繼提出了一些完善方法,并在電力系統(tǒng)中得到了廣泛地應(yīng)用,如將對數(shù)壁壘函數(shù)與牛頓法結(jié)合起來應(yīng)用到非線性規(guī)劃問題[4,5]中等。在內(nèi)點法的基礎(chǔ)上,M ehrotra提出了基于預(yù)測校正的原對偶內(nèi)點法[6],由于該方法在計算中引入預(yù)測、校正環(huán)節(jié),故可以計及互補松弛條件的二次性和動態(tài)地確定向心參數(shù)的取值,從而能夠較好地協(xié)調(diào)解的最優(yōu)性及可行性之間的關(guān)系,改善算法的收斂性能。

2 負荷裕度指標的定義

負荷裕度指標是評估電力系統(tǒng)靜態(tài)電壓穩(wěn)定性的有效手段之一。它是指從系統(tǒng)當前運行點到電壓崩潰點的“負荷距離”,具有較好的線性和直觀性。它可以使運行人員清晰地了解系統(tǒng)目前所處的狀態(tài),并提前采取有效措施避免電壓崩潰的發(fā)生。

如圖1所示的P-U曲線,當前系統(tǒng)的運行點是PL0,即表示當前系統(tǒng)的負荷。Pcr對應(yīng)著系統(tǒng)電壓崩潰的臨界點,即代表系統(tǒng)所能承受的極限負荷量。Pcr-PL0即為該系統(tǒng)的靜態(tài)電壓穩(wěn)定負荷裕度。由此可以定義系統(tǒng)的相對負荷裕度λ為:

由P-U曲線可以看出,當系統(tǒng)的負荷裕度較大時,系統(tǒng)的電壓水平也較高;同樣,當系統(tǒng)接近電壓崩潰臨界點時,系統(tǒng)裕度很小,系統(tǒng)電壓水平也相應(yīng)下降。利用此電壓穩(wěn)定評價指標,可以清晰地表示出系統(tǒng)電壓的穩(wěn)定水平。

圖1 P-U曲線

3 無功規(guī)劃模型

以新增無功容量總和最小作為目標函數(shù),無功規(guī)劃模型可用式(2)表示:

(1)等式約束方程h(x)

式中:Ui——節(jié)點 i的電壓幅值;θi——節(jié)點 i的相角; Bij——節(jié)點i與節(jié)點 j之間的互電納;Gij——節(jié)點 i與 j之間的互電導(dǎo);Qc——補償容量;λcr——臨界點的負荷增量;Qci——節(jié)點i補償無功功率的容量;SB——所有節(jié)點的集合;bp——有功負荷增長方向;bq——無功負荷增長方向;SG——無功源集合;PGi——第i個節(jié)點發(fā)電機的有功;QGi——第i個節(jié)點發(fā)電機的無功;PLi——第i個節(jié)點負荷的有功;QLi——第i個節(jié)點負荷的無功;Qci——第i個節(jié)點的無功補償量。

4 模型的求解

采用預(yù)測-校正內(nèi)點法對模型式(2)至式(8)進行求解,構(gòu)造其拉格朗日函數(shù),如式(9):

式(9)中,變量μ＞0為擾動因子(或稱障礙常數(shù));l,u為松弛變量,且u＞0。

應(yīng)用牛頓法得到待求問題的修正方程為:

將式(12)至式(17)簡化后寫成矩陣,如式(18):

只需求解式(18)矩陣右下角分塊中的函數(shù)值,然后逐一回代,即可求出所有變量。求解方程組如式(19):

將所得解Δx,Δy回代到式(20),即可求出所有變量的解:

變量更新的迭代方程組為:

式中αd,αp分別為步長:

引入預(yù)測及校正環(huán)節(jié)動態(tài)確定向心參數(shù)的取值,并且利用互補松弛條件的二次性,預(yù)估計算結(jié)果,以較好地協(xié)調(diào)解的最優(yōu)性及可行性之間的關(guān)系,改善算法的收斂性能。

(1)預(yù)測環(huán)節(jié) 對式(12)至式(17)的修正方程,取 μ=0,可求得 ΔxAff,ΔyAff,ΔzAff,ΔlAff, ΔuAff,ΔwAff。

式中:Aff表示仿射方向,即擾動因子取0時的牛頓方向。

根據(jù)預(yù)測過程,可確定仿射方向下的修正步長,如式(28)及相應(yīng)的對偶間隙如式(29):

預(yù)測-校正原對偶內(nèi)點法根據(jù)當前點的對偶間隙與仿射方向下的對偶間隙間的關(guān)系動態(tài)給定向心參數(shù)值,根據(jù)Mehrotra的經(jīng)驗[6],σ一般為:

(2)校正環(huán)節(jié) 考慮到互補松弛條件式(16)及式(17)的二次性,將二次泰勒方程式展開,可表示為:

將 ΔZAff、ΔlAff、ΔWAff、ΔuAff分別代入,對ΔZΔl及ΔWΔu進行近似,對互補松弛條件進行校正,即:

預(yù)測-校正原對偶內(nèi)點法引入了預(yù)測-校正機制,充分利用了互補松弛條件的二次性,避開了向心參數(shù)的人為選取,同時也提高了收斂性。

5 IEEE標準算例仿真分析與應(yīng)用

將無功規(guī)劃模型和求解方法用于IEEE 14、30、57、118標準算例,求解各系統(tǒng)初始負荷裕度結(jié)果如表1所示。其中,負荷增長方式設(shè)定為根據(jù)初始負荷大小全網(wǎng)同比例增長;發(fā)電機出力方式設(shè)定為系統(tǒng)中所有發(fā)電機(平衡機除外)根據(jù)各自有功備用的大小按比例分攤負荷的增量;優(yōu)化目標為將系統(tǒng)原有的負荷裕度提升10%。

表1 IEEE標準算例負荷裕度求解結(jié)果

(1)IEEE 14節(jié)點無功補償方案為:10號節(jié)點的無功補償量為4.57 MW,13號節(jié)點的無功補償量為0.33 MW,14號節(jié)點的無功補償量為11.53MW,總計為16.43 MW。

(2)IEEE 30節(jié)點無功補償方案為:24號節(jié)點的無功補償量為3.85 MW,26號節(jié)點的無功補償量為4.22 MW,29號節(jié)點的無功補償量為1.01 MW,30號節(jié)點的無功補償量為6.28 MW,總計為15.36MW。

(3)IEEE 57節(jié)點無功補償方案為:30號的節(jié)點補償量為4.35 MW,33號節(jié)點補償量為2.91 MW,31號節(jié)點補償量為5.75 MW,總計為13.01MW。

(4)IEEE 118節(jié)點無功補償方案為:1號的節(jié)點補償量為28.28 MW,9號的節(jié)點補償量為86.67MW,總計為114.93MW。

上述無功補償量計算中所需時間隨節(jié)點數(shù)量的增多而線性增大,但是其增加的速率基本是多項式級的,由此證明內(nèi)點法適合解算大規(guī)模系統(tǒng)問題。

利用上述無功規(guī)劃模型,在國內(nèi)某個大型電力系統(tǒng)上驗證所提方法的有效性。負荷增長方式和發(fā)電機出力方式與上述定義相同,補償候選點集合為220 kV主變的低壓側(cè)節(jié)點,各節(jié)點無功補償下限為0,上限為主變?nèi)萘康?0%。

采用預(yù)測-校正原對偶內(nèi)點法對該電網(wǎng)進行無功規(guī)劃,補償方案如表2所示。

如表2所示,該系統(tǒng)共有38個節(jié)點需要安裝無功補償設(shè)備,補償容量總計為354.05 MW。

表2 國內(nèi)某系統(tǒng)無功補償方案 MW

算法收斂過程如圖2所示。

圖2 國內(nèi)某電網(wǎng)無功規(guī)劃互補間隙收斂過程

預(yù)測-校正原對偶內(nèi)點法經(jīng)過27步迭代,耗時118 s收斂。在保證系統(tǒng)負荷裕度達到設(shè)定閥值的前提下,圓滿地完成了以補償量最小為目標的無功規(guī)劃。

6 結(jié)語

通過建立無功規(guī)劃的最優(yōu)化模型,應(yīng)用預(yù)測-校正原對偶內(nèi)點法對其進行求解。在IEEE 14,IEEE 30,IEEE 57和IEEE 118系統(tǒng)和國內(nèi)某個大型電網(wǎng)上,驗證了上述模型和求解方法的有效性。

[1] 文學(xué)鴻.基于內(nèi)點優(yōu)化算法的電力系統(tǒng)靜態(tài)電壓穩(wěn)定分析與最優(yōu)控制[碩士].河海大學(xué),2007.

[2] 陳宏偉,張興凱,王 寬.電力系統(tǒng)無功優(yōu)化的研究現(xiàn)狀和展望[J].電氣應(yīng)用,2006(12).

[3] N.KARMARKAR.A new polynom ial-tim ealgorithm for linear programm ing[C],STOC'84 Proceedings of six teen th annual ACM symposium on Theory of computing,1984.

[4] 王錫凡.現(xiàn)代電力系統(tǒng)分析[M].北京:科學(xué)出版社,2003.

[5] 張伯明.高等電力網(wǎng)絡(luò)分析[M].北京:清華大學(xué)出版社,2007.

[6] S.MEHROTRA.On the Imp lementation of a Primal-Dual Interior Point Method[J],SIAM Jou rnal on Optim ization,Nov.1984.

[7] 劉沛津,谷立臣,韓行.基于內(nèi)點法與改進遺傳法的無功規(guī)劃優(yōu)化混合算法[J].電力系統(tǒng)保護與控制,2008(17).

[8] 鐘世霞,袁榮湘.內(nèi)點法在電力系統(tǒng)中的應(yīng)用述評[J].高電壓技術(shù),2005(12).

[9] 文學(xué)鴻,袁 越,鞠 平.靜態(tài)電壓穩(wěn)定負荷裕度分析方法比較[J].電力自動化設(shè)備,2008,28(5).

[10] 王佳賢,程浩忠,胡澤春.多負荷水平下的配電網(wǎng)無功優(yōu)化規(guī)劃方法[J].電網(wǎng)技術(shù),2008,(19).

[11] 吳 浩.電力系統(tǒng)電壓穩(wěn)定研究[博士].浙江大學(xué),2002.

[12] 韋 化,丁曉鶯.基于現(xiàn)代內(nèi)點理論的電壓穩(wěn)定臨界點算法[J].中國電機工程學(xué)報,2002,22(3):27-31.

[13] 程浩忠,吳浩著.電力系統(tǒng)無功與電壓穩(wěn)定性[M].北京:中國電力出版社.2004.