朱雪凌，張翠影，趙臣鵬，劉林飛

(華北水利水電大學(xué)，河南 鄭州450045)

牛頓法具有較好的收斂性，在解最優(yōu)潮流時(shí)必須用到Hessian 矩陣［1］的逆矩陣，其存儲(chǔ)量及計(jì)算量大，使問(wèn)題變得復(fù)雜，因而如何簡(jiǎn)化成為首要問(wèn)題.

1984年，臺(tái)灣學(xué)者Sun D I 等［2］提出應(yīng)用二次罰函數(shù)的牛頓法處理該問(wèn)題. 該算法不用區(qū)分狀態(tài)變量和控制變量，充分利用電力網(wǎng)絡(luò)的物理特征，運(yùn)用Hessian 矩陣的導(dǎo)納稀疏結(jié)構(gòu)，把等式約束條件和不等式約束條件［3］用Lagrange 乘子引入到目標(biāo)函數(shù)中，直接對(duì)拉格朗日函數(shù)的Karush-Kuhn-Tucker 條件［4］(簡(jiǎn)稱KKT 條件)進(jìn)行牛頓法迭代求解，不等式約束用二次罰函數(shù)來(lái)處理. 文中采用二次罰函數(shù)的牛頓法來(lái)求解最優(yōu)潮流，并經(jīng)試驗(yàn)驗(yàn)證了該方法具有很強(qiáng)的實(shí)用性及經(jīng)濟(jì)性.

1 牛頓法的數(shù)學(xué)模型

1.1 非線性規(guī)劃的數(shù)學(xué)模型

典型的非線性規(guī)劃問(wèn)題［5］就是求解目標(biāo)函數(shù)的極大值或極小值問(wèn)題，文中所求的是極小值，數(shù)學(xué)模型可表示為:

1.2 牛頓法的描述

只考慮等式約束g(x)= 0 時(shí)，Lagrange 函數(shù)可表示為

其中λ 是Lagrange 乘子［6］，

根據(jù)庫(kù)恩－ 塔克條件［7］，在極小值點(diǎn)(x*，λ*)進(jìn)行Taylor 展開(kāi):

將二次項(xiàng)及高次項(xiàng)忽略，式(5)變?yōu)?/p>

式中H 和J 分別為Hessian 和Jacobian 矩陣.

將等式約束g(x)= 0 在變量初始值x0處進(jìn)行Taylor 展開(kāi):

忽略二次項(xiàng)與高次項(xiàng)得:

由式(6)和式(7)得:

式(8)則為求等式約束非線性規(guī)劃問(wèn)題的牛頓修正方程式［8］.而不等式約束條件h(x)≥0，用二次罰函數(shù)［9］來(lái)處理，擴(kuò)展后的Lagrange 函數(shù)表示為

式中:Ci為罰因子［10］;i 為不等式約束的個(gè)數(shù).把

作為擴(kuò)展目標(biāo)函數(shù)，考慮不等式約束后的牛頓修正方程為

可見(jiàn)，不等式約束只影響Hessian 矩陣系數(shù)和等式的右側(cè).

2 最優(yōu)潮流的數(shù)學(xué)模型

2.1 最優(yōu)潮流

最優(yōu)潮流(OPF)問(wèn)題［11］是一個(gè)典型的帶約束條件的非線性優(yōu)化問(wèn)題，進(jìn)行最優(yōu)潮流計(jì)算時(shí)，一般以系統(tǒng)發(fā)出有功、無(wú)功成本最小為目標(biāo)函數(shù)，其數(shù)學(xué)模型為

式中fpi(Pgi)，fqi(Qgi)為機(jī)組i 的燃料耗費(fèi).

2.2 約束條件

等式約束條件為

式中:Pgi，PLi分別為機(jī)組i 有功出力和有功負(fù)荷;Qgi，QLi分別為機(jī)組i 無(wú)功出力和無(wú)功負(fù)荷;P(V，θ)，Q(V，θ)分別為有功和無(wú)功網(wǎng)損. 式(12)和式(13)也是節(jié)點(diǎn)潮流方程［12］.

2.3 不等式約束

不等式約束條件為

3 算法步驟

算法步驟如下:

1)輸入原始數(shù)據(jù)，給出初始值θ0，V0，λP0，λQ0;

2)節(jié)點(diǎn)進(jìn)行優(yōu)化排隊(duì);

3)形成Hessian 矩陣和Jacobian 矩陣，進(jìn)行懲罰修正;

4)求解修正方程，得出:

5)求得步驟4 中結(jié)果看是否符合庫(kù)恩－塔克條件，若符合則結(jié)束運(yùn)算，否則返回到第2 步;

6)停止運(yùn)算.

4 實(shí)例計(jì)算

以IEEE14［14］節(jié)點(diǎn)標(biāo)準(zhǔn)系統(tǒng)為例，運(yùn)用MATLAB 編程進(jìn)行最優(yōu)潮流計(jì)算，所得支路節(jié)點(diǎn)和母線最優(yōu)潮流結(jié)果見(jiàn)表1和表2. 該算法求得的最優(yōu)潮流收斂時(shí)間在5. 52 s 以內(nèi)，系統(tǒng)的發(fā)電成本為8 081.53 |S/h. 由此可以看出，該算法收斂速度較快，求得的發(fā)電成本較低.

表1 支路節(jié)點(diǎn)最優(yōu)潮流計(jì)算結(jié)果

表2 母線最優(yōu)潮流計(jì)算結(jié)果

5 結(jié) 語(yǔ)

由試驗(yàn)數(shù)據(jù)看出:對(duì)于復(fù)雜的電力系統(tǒng)最優(yōu)潮流問(wèn)題，牛頓法可以較為精確地求出計(jì)算的結(jié)果;用二次罰函數(shù)處理不等式約束條件，使復(fù)雜問(wèn)題簡(jiǎn)單化.同時(shí)，二次罰函數(shù)的牛頓法的收斂性較好，運(yùn)算速度較快，求得的發(fā)電成本較低，具有很強(qiáng)的經(jīng)濟(jì)性與實(shí)用性，適合求解大系統(tǒng)的最優(yōu)潮流問(wèn)題.

［1］李勝淵.牛頓法最優(yōu)潮流算法的研究及實(shí)踐［D］.北京:中國(guó)電力科學(xué)研究院，1996.

［2］Sun D T，Ashley B. Optimal power flow by Newton approach［J］.IEEE PAS，1984，103(10):401－410.

［3］劉明波.大電網(wǎng)最優(yōu)潮流計(jì)算［M］.北京:科學(xué)出版社，2010:385－390.

［4］范宏，韋化.基于擾動(dòng)KKT 條件的原始－對(duì)偶內(nèi)點(diǎn)法和分界點(diǎn)法的最優(yōu)潮流研究［J］.電力自動(dòng)化設(shè)備，2004，24(5):5－9.

［5］胡運(yùn)洪.非線性最優(yōu)化問(wèn)題及其算法研究［J］.運(yùn)城學(xué)院學(xué)報(bào)，2003，21(3):9－10.

［6］王永剛，彭世康，靳現(xiàn)林. 牛頓法最優(yōu)潮流的改進(jìn)［J］.繼電器，2003，31(3):1－5.

［7］葛結(jié)根.庫(kù)恩－塔克定理在微觀經(jīng)濟(jì)學(xué)中的應(yīng)用［J］.統(tǒng)計(jì)與決策，2004(7):47－53.

［8］顏偉，黃正波，余娟，等.牛頓法潮流計(jì)算的高效綜合稀疏技術(shù)［J］.中國(guó)電力，2010，43(7):19－23.

［9］趙普泉，侯志儉，吳際舜. 牛頓法最優(yōu)潮流算法中離散控制量的新處理方法［J］. 電力系統(tǒng)自動(dòng)化，1999，23(23):37－58.

［10］張力平，何大愚，朱太秀.牛頓法最優(yōu)潮流與最優(yōu)無(wú)功補(bǔ)償［J］.中國(guó)電機(jī)工程學(xué)報(bào)，1987，7(1):11－18.

［11］林睦綱.電力系統(tǒng)最優(yōu)潮流新算法的研究［D］.長(zhǎng)沙:長(zhǎng)沙理工大學(xué)，2005.

［12］陳珩.電力系統(tǒng)穩(wěn)態(tài)分析［M］. 北京:中國(guó)電力出版社，2007:123－135.

［13］邱道尹，張凌云，顧波，等. 基于徑向基神經(jīng)網(wǎng)絡(luò)的雙饋風(fēng)力發(fā)電機(jī)低電壓穿越控制研究［J］. 華北水利水電學(xué)院學(xué)報(bào)，2013，34(6):100－105.

［14］IEEE Committee Report. IEEE reliability test system［J］.IEEE Trans PWRS，1979，98(6):2047－2054.