案例4：電力系統(tǒng)的電壓控制方法

案例4：電力系統(tǒng)的電壓控制方法

一、案例正文

1.1電壓穩(wěn)定的概念

電力工業(yè)的市場化發(fā)展和各級電網(wǎng)的互聯(lián)互通在優(yōu)化資源配置

和提高經(jīng)濟(jì)效益的同時(shí)，對電力系統(tǒng)的電壓穩(wěn)定運(yùn)行提出了更大的挑

戰(zhàn)。電力系統(tǒng)的電壓穩(wěn)定運(yùn)行是保證電力系統(tǒng)發(fā)展的基礎(chǔ)，對電力系

統(tǒng)的電壓穩(wěn)定性研究是電力系統(tǒng)規(guī)劃設(shè)計(jì)和運(yùn)行研究領(lǐng)域中非常重

要的課題。國際大可網(wǎng)會議（CIGRE）和國際電氣與電子工程師學(xué)會

電力工程分會（IEEE/PES）認(rèn)為：電壓穩(wěn)定性是指電力系統(tǒng)遭受擾動

后母線維持穩(wěn)定電壓的能力。我國在《電力系統(tǒng)安全穩(wěn)定導(dǎo)則》中將

電壓穩(wěn)定定義為：電力系統(tǒng)受到小擾動后，系統(tǒng)電壓能夠保持或恢復(fù)

到允許的范圍內(nèi)，不發(fā)生電壓崩潰的能力。

1.2電壓穩(wěn)定分析方法

潮流方程是對電力系統(tǒng)電壓穩(wěn)定進(jìn)行靜態(tài)分析的核心，以擴(kuò)展潮

流方程為基本模型，通過分析臨界運(yùn)行狀態(tài)下系統(tǒng)的各類特征判斷當(dāng)

前狀態(tài)是否穩(wěn)定，根據(jù)不同特征選擇不同的分析方法計(jì)算臨界點(diǎn)，評

判電壓穩(wěn)定的靜態(tài)分析方法主要有:連續(xù)潮流法靈敏度分析法、非線

性規(guī)劃法、奇異值分解法和分岔分析法。

1.2.1連續(xù)潮流法

通過連續(xù)潮流法可有效分析電力系統(tǒng)的靜態(tài)電壓穩(wěn)定性，連續(xù)潮

流法的計(jì)算步驟主要包括預(yù)測、校正、參數(shù)化和補(bǔ)償控制。該方法常

用于電壓穩(wěn)定及負(fù)荷能力的分析,能夠比較準(zhǔn)確地計(jì)算出負(fù)荷裕度。

通過連續(xù)潮流計(jì)算不僅能夠得到系統(tǒng)的電壓穩(wěn)定臨界點(diǎn)，而且還可以

得到完整的負(fù)荷節(jié)點(diǎn)PV曲線?；诜蔷€性潮流方程的連續(xù)潮濟(jì)法計(jì)

算速度比較慢，一種改進(jìn)的連續(xù)潮流法，通過結(jié)合電力網(wǎng)絡(luò)的結(jié)構(gòu)特

性引入新的狀態(tài)變量，在每次迭代時(shí)不需要重新計(jì)算雅可比矩陣，從

而減少了大量計(jì)算，提高了連續(xù)潮流法計(jì)算的整體速度；還有基于特

勒根定理的改進(jìn)連續(xù)潮流模型，通過對步長進(jìn)行控制，有效提高了計(jì)

算速度。

1.2.2靈敏度分析法

靈敏度分析法在電力系統(tǒng)靜態(tài)電壓穩(wěn)定分析中應(yīng)用廣泛，靈敏度

分析法基于電力系統(tǒng)潮流方程計(jì)算各物理量之間的微分關(guān)系，以此來

研究系統(tǒng)的電壓穩(wěn)定性，如dQi/d%、dPQi/dVI。通過靈敏度分析有

助于理解系統(tǒng)元件中冬控制變量對系統(tǒng)電壓穩(wěn)定性的影響，以研究負(fù)

荷裕度與各控制變量之間的關(guān)系。靈敏度方法求解方便，物理概念清

晰，計(jì)算量小，在簡單電力系統(tǒng)中應(yīng)用取得了不錯(cuò)的效果。有文獻(xiàn)提

出靜態(tài)電壓穩(wěn)定預(yù)防控制在線決策的核心環(huán)節(jié)是負(fù)荷裕度及其與控

制變量間靈敏度關(guān)系的快速計(jì)算，通過常規(guī)潮流方程總結(jié)了調(diào)控量與

狀態(tài)量間的靈敏度關(guān)系，基于戴維南等值參數(shù)辨識方法并結(jié)合等值參

數(shù)與負(fù)荷裕度之間的數(shù)學(xué)關(guān)系，推演出調(diào)控量與負(fù)荷裕度間的解析關(guān)

聯(lián)靈敏度，算例分析證明了靈敏度方法的有效性，通過靈敏度分析對

冬控制變量進(jìn)行調(diào)整可有效提高系統(tǒng)的靜態(tài)電壓穩(wěn)定性。

1.2.3非線性規(guī)劃法

非線性規(guī)劃法是求解最優(yōu)化問題的堂用方法，屬于運(yùn)籌學(xué)的重要

部分，其中包括目標(biāo)函數(shù)和非線性約束條件。非線性規(guī)劃法在電力系

統(tǒng)中的應(yīng)用相當(dāng)廣泛，也是求取靜態(tài)電壓穩(wěn)定臨界點(diǎn)的一種方法，常

規(guī)非線性規(guī)劃法求解靜態(tài)電壓穩(wěn)定臨界點(diǎn)時(shí)沒有考慮負(fù)荷增長時(shí)冬

發(fā)電機(jī)組有功出力的增長方式，導(dǎo)致計(jì)算結(jié)果過于理想。有一種結(jié)合

動態(tài)潮流思想，提出一種考慮發(fā)電機(jī)有功出力增長方式約束的動態(tài)非

線性規(guī)劃法模型，結(jié)合實(shí)際選定的發(fā)電機(jī)有功出力增長方式來分擔(dān)系

統(tǒng)全部有功不平衡功率，使得利用動態(tài)非線性規(guī)劃法求得的電力系統(tǒng)

負(fù)荷裕度結(jié)果更加貼近實(shí)際。非線性規(guī)劃法也可有效避免臨近電壓穩(wěn)

定極限時(shí)潮流雅可比矩陣奇異及潮流不收斂的情形。隨著系統(tǒng)規(guī)模的

擴(kuò)大，約束方程數(shù)的急劇增加使非線性規(guī)劃求解的難度增加，使該方

法的計(jì)算規(guī)模受到限制。

1.2.4奇異值分解法

奇異值分解是有關(guān)矩陣問題的強(qiáng)有力的計(jì)算工具，已在眾多領(lǐng)域

中得到應(yīng)用。該方法由潮流雅可比矩陣的行列式符號決定所研究系統(tǒng)

是穩(wěn)定的或不穩(wěn)定的，可以將潮流方程的雅可比矩陣進(jìn)行奇異值分解,

潮流雅可比矩陣的最小奇異值作為靜態(tài)電壓穩(wěn)定性的指標(biāo)，將其加入

目標(biāo)函數(shù)構(gòu)造優(yōu)化模型。以最小奇異值作為靜態(tài)電壓穩(wěn)定性的指標(biāo)計(jì)

算量小，實(shí)現(xiàn)簡單。其中最小奇異值大小為當(dāng)前運(yùn)行點(diǎn)和電壓穩(wěn)定極

限之間的距離，左奇異向量和右奇異向量中的最大元素分別指示功率

注入最靈敏的方向和電壓幅值。

1.2.5分岔分析法

從數(shù)學(xué)的角度來看，分岔分析法研究的主要是非線性方程中參數(shù)

統(tǒng)運(yùn)行點(diǎn)接近鞍結(jié)分岔點(diǎn)時(shí)數(shù)值發(fā)散的問題，該方法通過在常規(guī)潮流

方程中添加連續(xù)性參數(shù)，可以沿著負(fù)荷增長的方向逐步追蹤P-V曲線

的軌跡達(dá)到鼻點(diǎn)為止（如圖1所示），電力系統(tǒng)的負(fù)荷裕度可定義為當(dāng)

負(fù)荷按照某種增長方式逼近電壓崩潰點(diǎn)時(shí)：系統(tǒng)當(dāng)前運(yùn)行點(diǎn)到電壓崩

潰點(diǎn)的距離。

圖IP-V曲線示意圖

參數(shù)化的連續(xù)潮流方程可以表示為：

f（x,a,p）=o（1）

式中，X為狀態(tài)變量，a為負(fù)荷裕度，p表示系統(tǒng)指定參數(shù)。節(jié)點(diǎn)i的連

續(xù)潮流方程展開式為：

AE=PGi-PLi-Vi^=i匕（Gijcos%+Bijsin辦）+入（△七，一△

P〃）=0（2）

n

△Qi=Qa~QLI一匕W，（Gi，sin%—cos8訂）+入（一△QLi）=0

7=i

（3）

式中，匕和,分別表示節(jié)點(diǎn)i和節(jié)點(diǎn)/的電壓幅值，PGi和QGi分別表示節(jié)

點(diǎn)i上常規(guī)電源的有功功率和無功功率，Qij表示節(jié)點(diǎn)i和節(jié)點(diǎn)，的電壓

相角差值，2和Q〃，分別表示節(jié)點(diǎn)i上的負(fù)荷功率，Gu和%分別表示

節(jié)點(diǎn)i和節(jié)禹之間的網(wǎng)絡(luò)導(dǎo)納的實(shí)部和虛部，APGI表示節(jié)點(diǎn)i的發(fā)電機(jī)

注入變化趨勢，AP〃和表示節(jié)點(diǎn)i的負(fù)荷變化趨勢。

L3.3鞍結(jié)點(diǎn)分岔點(diǎn)的靈敏度

在鞍結(jié)點(diǎn)分岔點(diǎn)(%*〃*),系統(tǒng)的雅可比知陣奇異，因此系統(tǒng)的雅

可比矩陣存在零特征值，且存在一個(gè)相對于該零特征值的左特

征向量3(%*,不)滿足：

①此人*)瞿|=0(4)

■I-

在分岔點(diǎn)對于不同的參數(shù)p均滿足上式，將公式(1)進(jìn)行線性化處

理可得：

+察△入+及Ap=0(5)

dxdXdpl

式中，d/ya%為系統(tǒng)的雅可比矩陣，為潮流方程對負(fù)荷增長系

數(shù)a的導(dǎo)數(shù)a//3p為潮流方程對任意參數(shù)P的導(dǎo)數(shù)。將零特征值對應(yīng)的

左特征向量3(%*"*)左乘上式，可得：

38，入*)修△入喘”)=。⑹

從而得出鞍結(jié)點(diǎn)分岔點(diǎn)處負(fù)荷裕度a對任意參數(shù)P的靈敏度公式

為為:

加=益=一3(”入埸/3(X*,入嚏(7)

1.3.3.1線路投切的靈敏度計(jì)算

在線路投切的靈敏度計(jì)算中，任意參數(shù)p表示開斷的線路參數(shù)，

包括線路電納、對地容納以及線路電導(dǎo)，AP表示開斷線路的負(fù)參數(shù)向

量，0〃即表示潮流方程對開斷線路參數(shù)的導(dǎo)數(shù)矩陣。根據(jù)連續(xù)潮流

方程可得:如果開斷線路i-只有節(jié)點(diǎn)i和節(jié)點(diǎn)/的功率平衡方程與開

斷線路參數(shù)有關(guān)系，因此矩陣。△Pi/OPij中只有?！魈?。句、MQi/

dPij、d^Qj/dPtj.d△號/。均非零。這4個(gè)非零元素分別為：

n

△Pl=Pot—匕*W，*(G〃c。s%+%sine；j)+r^PGi-APQ

;=i

=0

(8)

n

AQi=Qoi-ViW4(G〃sin%—BijCos魴)+矛(-AQQ=0

J=I

⑼

n

%=Poj-Vt24(%c。s%+Bjtsin穌)+矛(APGJ-"切)=0

；=i

(io)

n

△Qj=Qoj-ViW4(%sine-t-BjiCOs所)+a*(-AQu)=0

i=l

(II)

式中，Poi、Qoi是初始狀態(tài)下節(jié)點(diǎn)i注入的有功功率和無功功率，Poj、

Qoj是初始狀態(tài)下節(jié)點(diǎn)j注入的有功功率和無功功率。

線路i-j斷開后導(dǎo)致負(fù)荷裕度有所變化，其變化量△入叮由式(12)計(jì)

算：

pOAPtcQdAQ]pdAPj?八dAQ,

△電=3西邸ij+3兩5/+3'麗即。+3，標(biāo)啊i

（12）

式中，M和的分別表示節(jié)點(diǎn)i和節(jié)點(diǎn)/的有功功率平衡方程，AQi和

△Q,分別表示節(jié)點(diǎn)i和節(jié)點(diǎn)/的無功功率平衡方程,表示線路開斷前

后輸電線路參數(shù)變化量，s=-co（x*,r）/（co（x*,r）（a//dA））,砂為

向量小對應(yīng)節(jié)點(diǎn)i尢功功率平衡方程中的元素，為向量小對應(yīng)節(jié)點(diǎn)i

有功功率平衡方程中的元素，方Q/為向量G對應(yīng)節(jié)點(diǎn)/無功功率平衡方

程中的元素，5JP/為向量G對應(yīng)節(jié)點(diǎn)/有功功率平衡方程中的元素。

1.3.3.2變壓器分接頭的靈敏度計(jì)算

通過計(jì)算負(fù)荷裕度對變壓器分接頭的靈敏度，可有效篩選出提高

負(fù)荷裕度效果較強(qiáng)的變壓器，對篩選出的變壓器進(jìn)行調(diào)整可有效提高

系統(tǒng)的靜態(tài)電壓穩(wěn)定性。關(guān)于負(fù)荷裕度對變壓器分接頭的靈敏度計(jì)算,

計(jì)算方法與線路投切類似，將變壓器等效成7T型電路如圖2所示。

圖2變壓器丸型等效電跖

變壓器分接頭位置可以表示為：

k=1+kp?t(13)

式中,k為變壓器分接頭；芍為變壓器調(diào)整時(shí)單位檔位的大小，其數(shù)值通

常為0.0125;t為變壓器分接頭所處的檔位，通常情況下

加[—8,—7,…7,8],若升高一檔則t+1,若降低一檔則t—1。假設(shè)變

壓器可以連續(xù)調(diào)節(jié)，則dk/dt=0.0125o

變壓器節(jié)點(diǎn)i處的自導(dǎo)納為:

GCI;n_(kT)y/工如_”

ii+必t-----yT(14)

由式(14)可得，節(jié)點(diǎn)i處的自導(dǎo)納與分接頭位置無關(guān)，則其對t的

偏導(dǎo)為:

㈣=0(15)

九

嗎=0(16)

和

變壓器支路的互導(dǎo)納為-笠，其對珀勺偏導(dǎo)為：

ik

dGijd(G)dkkpG

--=-------Lr----(17)

UAtLLdLikv(dXLA/"2

dBU=*)dkR=M

dtdkdtk2(18)

變壓器節(jié)點(diǎn)i處的自導(dǎo)納為：

cIjD_(l-k)yr,yr_yr

Gjj+jBjj+工—記(19)

對t的偏導(dǎo)為:

dGJJ_d信)dkc=-2kpG

dtdkdtUk3(20)

dBJj=壯忌)dkB="B

dtdkdtUk3(21)

根據(jù)上式關(guān)于節(jié)點(diǎn)導(dǎo)納陣的偏導(dǎo)計(jì)算結(jié)果，結(jié)合公式(8)?(11)推導(dǎo)

功率平衡方程對變壓器檔位的偏導(dǎo)為：

誓=一碧匕匕(Geos勿+BsinB^(22)

Uviv

嘿=-%VM(Gcos9ij_8sin%)(23)

寫2=察-*V"j(Gcoseij+Bsin*j)(24)

=一魯叫2-今匕匕.(Gsin%?-BcosBtj)(25)

1.4滿足靜態(tài)電壓穩(wěn)定增強(qiáng)要求的發(fā)電機(jī)最佳調(diào)整方法

1.4.1發(fā)電機(jī)最佳調(diào)整模型

優(yōu)化發(fā)電機(jī)出力能夠提升電力系統(tǒng)的靜態(tài)電壓穩(wěn)定性，欲將系統(tǒng)

的負(fù)荷裕度提高指定百分比％%,計(jì)算需要調(diào)整的發(fā)電機(jī)并確定其有

功出力。

1.4.1.1目標(biāo)函數(shù)

①參與出力調(diào)整的發(fā)電機(jī)臺數(shù)最少

mi咯日I(26)

式中，尤表示節(jié)點(diǎn)i處的發(fā)電機(jī)是否參與出力調(diào)整，心￡{0,1},若該發(fā)

電機(jī)參與出力調(diào)整，則叫=1(反之，匕=0)；4為全部發(fā)電機(jī)所在

節(jié)點(diǎn)組成的集合。

②發(fā)電成本最小

mink?(PGi)(27)

式中，C《)為節(jié)點(diǎn)i處發(fā)電機(jī)的發(fā)電成本函數(shù)；h為節(jié)點(diǎn)i處發(fā)電機(jī)調(diào)

整后的有功出力，即PG?=P&i+k/PGi；其中，P&為初始狀態(tài)下節(jié)點(diǎn)

i處發(fā)電機(jī)的有功出力，APG?為所求的有功出力調(diào)整量。

1.4.1.2約束條件

①負(fù)荷裕度約束

目標(biāo)函數(shù)(式(26)?式(27))需在負(fù)荷裕度滿足指定條件下，

即

A>(1+x%)〃(28)

式中，大丁和九分別為優(yōu)化前和優(yōu)化后系統(tǒng)的負(fù)荷裕度；x%為負(fù)荷裕度

提升需求，由運(yùn)行人員根據(jù)系統(tǒng)運(yùn)行情況斃前指定。

②連續(xù)潮流平衡約束

對節(jié)點(diǎn)i(im=1,23…m)，應(yīng)滿足連續(xù)潮流平衡方程：

PGI+^PGi~PLi~匕￡為Vj^Gijeos%+B^sin%)+

-APQ=0(29)

m

QQLI-匕分，(Gi/sin%-Bjjcos%)+乂―AQQ=0

J=i

(30)

式中，P8和盤i表示初始狀態(tài)下節(jié)點(diǎn)i處發(fā)電機(jī)的有功出力與無功

出力；和表示節(jié)點(diǎn)i處的負(fù)荷功率；7幾為系統(tǒng)的總節(jié)點(diǎn)個(gè)數(shù)。

③發(fā)電機(jī)出力調(diào)整約束

調(diào)整過程中，忽略網(wǎng)絡(luò)功率損耗，則參與調(diào)節(jié)的發(fā)電機(jī)有功功率

調(diào)整量總和為零，即

￡畛33=0(31)

同時(shí)，調(diào)整過程應(yīng)滿足有功出力上下限約束，即

<PGI+KGI<2i,max，i€A(32)

式中，PGi.max和PGLznin分別為節(jié)點(diǎn)i處的發(fā)電機(jī)有功出力允許的上下

限值。

④靜態(tài)安全運(yùn)行約束

各節(jié)點(diǎn)電壓幅值應(yīng)滿足約束：

匕,min——匕,max,，=1,2,3,…6(33)

式中，匕、％max和匕,m譏分別為節(jié)點(diǎn)t的電壓幅值及允許的電壓上下限

值。

輸電線路〃?上視在功率的取值范圍:

\Sij\<Sjmax，i，j=1,2,3,??.zn（34）

式中，片和Sjm”分別為連接在節(jié)點(diǎn)i和/之間的線路容量和允許的線

路容量限值。

142模型轉(zhuǎn)化

根據(jù)第二章靈敏度計(jì)算方法可以得到鞍結(jié)分岔點(diǎn)（％*,不）處的負(fù)

荷裕度入對各臺發(fā)電機(jī)有功出力p的靈敏度：

.=.=-3（%*"*境/3。/堵（35）

ApUpUA

式中，力為負(fù)荷裕度對節(jié)點(diǎn)i處發(fā)電機(jī)有功出力的靈敏度。力>0表示

若增加節(jié)點(diǎn)i處發(fā)電機(jī)有功出力可以提高負(fù)荷裕度；反之，/<0表示

節(jié)點(diǎn)i處發(fā)電機(jī)減少出力可以提高負(fù)荷裕度；而當(dāng)歷《0,則節(jié)點(diǎn)i處

發(fā)電機(jī)有功出力調(diào)整對負(fù)荷裕度的影響不明顯，因此本文不考慮靈敏

度接近于0的發(fā)電機(jī)參與調(diào)整。發(fā)電機(jī)有功出力最大調(diào)整量為：

ADfPGi,max—^Gi>。.A

△PG「=髭"-PG…V。L"（36）

盡管部分發(fā)電機(jī)靈敏度數(shù)值較大，但其可調(diào)有功區(qū)間較小，因此

本章采用式（37）來估算負(fù)荷裕度的最大變化量，并以此為依據(jù)將發(fā)

電機(jī)進(jìn)行排序。

△%=7］心PGmax，1EA附7）

式中，△々是通過靈敏度數(shù)值與發(fā)電機(jī)有功出力最大可調(diào)整量估算出

的負(fù)荷裕度變化量。

本文所提模型中發(fā)電機(jī)是否參與出力調(diào)整可表示為（M整數(shù)變量,

發(fā)電機(jī)有功功率調(diào)整量有界且連續(xù)，因此該問題是一個(gè)多目標(biāo)混合整

數(shù)非線性規(guī)劃問題。為降低該問題的計(jì)算復(fù)雜性，本文將多目標(biāo)規(guī)劃

問題轉(zhuǎn)化成兩個(gè)子問題進(jìn)行求解，算法步驟如圖3所示。

子問題1：采用靈敏度方法快速估算不同發(fā)電機(jī)調(diào)整臺數(shù)的負(fù)荷

裕度提升區(qū)間，求解所提問題的整數(shù)變量；

子問題2：應(yīng)用靈敏度法將原問題的非線性等式約束進(jìn)行線性化

處理，轉(zhuǎn)化成線性規(guī)劃問題來計(jì)算最優(yōu)解。

圖3滿足負(fù)荷裕度要求的發(fā)電機(jī)調(diào)整流程圖

1.4.3發(fā)電機(jī)最佳調(diào)整問題的求解

1.發(fā)電機(jī)調(diào)整臺數(shù)計(jì)算

該子問題采用基于靈敏度的方法估算發(fā)電機(jī)調(diào)整臺數(shù)，前文已經(jīng)

介紹了靈敏度方法及裕度估算，因此，為實(shí)現(xiàn)參與調(diào)整發(fā)電機(jī)的計(jì)算,

本文構(gòu)造了一個(gè)基于靈敏度的發(fā)電機(jī)調(diào)整臺數(shù)估算子問題的數(shù)學(xué)模

型，根據(jù)負(fù)荷裕度提升需求確定需要調(diào)整的發(fā)電機(jī)集合。

max2g匕/APGi(38)

stEE如/IVE-kg<￡3此忻H(39)

Eiwa履=n,nG{2,3,4,}(40)

式中，九為參與調(diào)整的發(fā)電機(jī)臺數(shù)；式(39)保證了在選取發(fā)電機(jī)增

強(qiáng)負(fù)荷裕度時(shí)其靈敏度數(shù)值要有止有負(fù)，式(31)保證了參與調(diào)整的

發(fā)電機(jī)臺數(shù)最少為2。

上述子問題解決的是：當(dāng)采取兩臺發(fā)電機(jī)(即n=2)參與出力調(diào)

整時(shí)，通過上述模型可得到：①估算出系統(tǒng)最大負(fù)荷裕度變化量；②

參與調(diào)整的兩臺發(fā)電機(jī)節(jié)點(diǎn)組成的集合B；③各發(fā)電機(jī)有功出力調(diào)節(jié)

量APGi(iC8)。根據(jù)計(jì)算出的發(fā)電機(jī)節(jié)點(diǎn)集合及有功出力調(diào)節(jié)量，運(yùn)

用連續(xù)潮流法計(jì)算出此時(shí)具體的負(fù)荷裕度變化量△鼠以此類推，當(dāng)采

取多臺發(fā)電機(jī)(即n>2)參與調(diào)整時(shí)，均可用該方法求解負(fù)荷裕度最

大變化量。

根據(jù)上述子問題的解可以得到不同臺數(shù)發(fā)電機(jī)調(diào)整與負(fù)荷裕度

提升范圍的對應(yīng)關(guān)系，以此作為依據(jù)，運(yùn)行人員可根據(jù)提升需求選擇

對應(yīng)的發(fā)電機(jī)，求解出所提問題的整數(shù)變量，將對應(yīng)的發(fā)電機(jī)節(jié)點(diǎn)集

合8傳送至下一階段進(jìn)行具體調(diào)節(jié)量的計(jì)算。

2.發(fā)電機(jī)調(diào)節(jié)量計(jì)算

該子問題以子問題1的解作為給定值，求解發(fā)電機(jī)有功出力的具

體調(diào)節(jié)量。本文采用靈敏度方法，將負(fù)荷裕度約束線性化處理，針對

所提問題構(gòu)造如下數(shù)學(xué)模型：

s?t.SiGB/APGI>ax%AT(41)

式中，a是補(bǔ)償因子，一般取為(1.01?1.03),它主要是為補(bǔ)償線性靈敏

度的不足而設(shè)計(jì)，后續(xù)仿真中a取1.02o

該子問題在數(shù)學(xué)上屬于線性混合整數(shù)規(guī)劃問題，利用現(xiàn)有的線性

規(guī)劃求解器，求解出發(fā)電機(jī)有功出力調(diào)整量APG&eB),并確定新的

發(fā)電機(jī)出力。

在本文中由于負(fù)荷裕度的計(jì)算是非線性的，因此需要檢查電力系

統(tǒng)在新的出力下是否滿足系統(tǒng)負(fù)荷裕