電力系統(tǒng)潮流計(jì)算高斯

一、高斯一一塞德?tīng)柗ǔ绷饔?jì)算

以導(dǎo)納矩陣為基礎(chǔ)的潮流計(jì)算。

設(shè)系統(tǒng)中有n個(gè)節(jié)點(diǎn)，其中有m個(gè)PQ點(diǎn)、n-（m+l）個(gè)PV節(jié)點(diǎn)和一

個(gè)平衡節(jié)點(diǎn)。平衡節(jié)點(diǎn)不參加迭代。

從方程式可以解出：

匕=^^^-￡亞］。（12-14）

Yuv>?=i

7*1

將上式改寫(xiě)成高斯一一塞德?tīng)柗ǖ碌袷剑?/p>

W=白/或-WzM”-￡叫。（12-15）

在用這個(gè)迭代公式時(shí)，PQ節(jié)點(diǎn)的功率是給定的，因此只要給出節(jié)

點(diǎn)電壓的初值匕?,可以進(jìn)行迭代計(jì)算。

對(duì)于PV節(jié)點(diǎn)，節(jié)點(diǎn)有功功率6和電壓幅值匕是給定的。但是節(jié)點(diǎn)

的無(wú)功功率只在迭代開(kāi)始時(shí)給出初值0?,此后的迭代值必須在迭代過(guò)

程中依次的算出。因此，在每一次迭代中，對(duì)于PV節(jié)點(diǎn)，必須作以下

幾項(xiàng)計(jì)算。

1、修正節(jié)點(diǎn)電壓

在迭代計(jì)算中，由公式（12-15）求得的節(jié)點(diǎn)電壓，其幅值不一定

等于給定的電壓幅值匕。為滿足這個(gè)條件，我們只保留節(jié)點(diǎn)電壓

的相位s*,而把其幅值直接取為給定值匕，即令

印二匕o（12-16）

2、計(jì)算節(jié)點(diǎn)無(wú)功功率

其計(jì)算公式為：

Q?=ImM⑹優(yōu)）］=1皿匕（幻（*匕邛川+￡＞尸）］（12-17）

7=1j=i

3、無(wú)功功率越線檢查

由上式算出的無(wú)功功率須按以下的不等式進(jìn)行檢驗(yàn)：

2mm＜中＜2"。（12-18）

如果。尸＞22，則令*=%；如果&&）＜%，則令中=“。

做完上述三項(xiàng)計(jì)算后，才應(yīng)用公式（12-15）計(jì)算節(jié)點(diǎn)電壓的新值。

平衡節(jié)點(diǎn)的電壓幅值和相位都是給定的，不必進(jìn)行迭代。

迭代收斂的判據(jù)是：

max{|"M）-片幻|}＜￡o（12-19）

迭代結(jié)束后，還要算出平衡節(jié)點(diǎn)的功率和網(wǎng)絡(luò)中的功率分布。輸

電線路功率的計(jì)算公式如下：（見(jiàn)圖12-2）

尾+閥=v,hj=匕2%+匕（匕-匕）%（12-20）

圖12-2支路功率計(jì)算

圖12-3是高斯-塞德?tīng)柗ǔ绷饔?jì)算的流程框圖。

輸入原始數(shù)據(jù)

形式節(jié)點(diǎn)‘導(dǎo)納矩陣

送節(jié)點(diǎn)電壓初值；PV節(jié)點(diǎn)無(wú)功初值

置迭代t十算k=0

，置節(jié)點(diǎn)計(jì)數(shù)i=1及AVniax=0

/是否平衡機(jī)節(jié)是

八占、、

▼否

＜是否PV節(jié)點(diǎn)＞----------

__________\A___________

用匕代v陽(yáng)由公式(12-17)計(jì)算°廣

。：—是＞。*223？是

否"I否_,_________\

用。.代H|用一代r用匕5代r

「LT

用公式(12?15)計(jì)算曠”卜

▼

A匕小”＞AJ?,置△匕s=AV；|U,)

，否

，是

用(12*4)求平衡節(jié)點(diǎn)功率，用(12-20)求全部線路功率

輸出‘結(jié)果

圖12-3高斯-塞德?tīng)柗ǔ绷饔?jì)算的流程框圖

二、牛頓拉夫遜法潮流計(jì)算

把牛頓法用于潮流計(jì)算。為此，先應(yīng)將潮流方程(12-4)的右端

展開(kāi)，并且分開(kāi)實(shí)部和虛部?！痉匠?12-4)匕+川=喧；皿(；1，2，?-，〃)】

>=i

采用直角坐標(biāo)，節(jié)點(diǎn)電壓表示為：

%=%+jfi

導(dǎo)納矩陣元數(shù)則表示為：

力=%+jf；

將上述表示式代入(12-4)式的右端，展開(kāi)并分出實(shí)部和虛部，

便得：

pi=eiX+fjt(G/j+B/p,

:：(12-39)

Qi=￡Z(。溝-8//),2電力+8汽)。

j=i》

按照節(jié)點(diǎn)的分類，PQ節(jié)點(diǎn)的有功功率和無(wú)功功率是給定的，第i

個(gè)節(jié)點(diǎn)的給定功率設(shè)為巳和&。假定系統(tǒng)中的第L2,…，m號(hào)節(jié)點(diǎn)

為PQ節(jié)點(diǎn)，對(duì)其中每一個(gè)節(jié)點(diǎn)可列寫(xiě)方程：

然嗎-ei￡(Gaj-B㈤芭力(G/j+B曲)=0,

j=ij=i(12-40)

△2=QIX-Q,=Qijf；X(G/廠B/j)+c,.￡(G應(yīng).+B/j)=0,

j=1J=l

(i=l,2,???,m),

PV節(jié)點(diǎn)的有功功率和節(jié)點(diǎn)電壓幅值是給定的。假設(shè)系統(tǒng)中的第

m+l,m+2,…n?l號(hào)節(jié)點(diǎn)為PV節(jié)點(diǎn)，則對(duì)其中每一個(gè)節(jié)點(diǎn)可以列寫(xiě)方

程：

岫北弋及-e$(G/j-B乩)-f；汽(G/,+B四)=0,

j=ij=i

△匕2=4—匕2=匕：一@2+/2)=0,(12-41)

(i=m+1,m+2,…，〃-l)0

第n號(hào)節(jié)點(diǎn)為平衡節(jié)點(diǎn)，其電壓匕,=儲(chǔ)+〃時(shí)給定的，故不參加迭

代。

式(12-40)和(12-41)總共包含了2(n-l)個(gè)方程，帶待求變

量有色,工勺,人,...明,裔也是2(n-l)個(gè)。我們還看到，方程(12-40)

和(12-41)已經(jīng)具備了方程組(12-8)的形式。因此，不難寫(xiě)出如下

的修正方程式：

AW=-JAV,

式中，

%

陰

△2.以

■■

AW=9AV=9

M向

M3

△加

M3

J為雅可比矩陣的各元素，可以對(duì)(12-40)和(12-41)式求偏

導(dǎo)數(shù)獲得。

當(dāng)，?打?時(shí)，

當(dāng)二警二一5工),

22=口-f,(12-43)

%%-以

處國(guó)=0。

陽(yáng)帶

當(dāng)J=i時(shí)，

r\\Pn

堂=_§(GE-%/)-G消-B",

G

~T7~=一￡(ijfj~旦汽)+瓦￡一GJ1,

0Ji/-I

辦2

=（GJj-4臼）+5咨—GJ,

de

i尸(12-44)

~t（G次一斗力）—G擊一用辦

然

7=1

oei

抽匕2

二一2/。

F

修正方程式（12-42）還可以編寫(xiě)成分塊矩陣的形式:

-AW14112…

△嗎，21^22???A.W-1

--,(12-45)

??

*?

△此1A-1.2???△匕L

LJ

式中，△叱和△匕都是二維列向量；々是2x2階方陣。

AV.=(12-46)

醞

對(duì)于PQ節(jié)點(diǎn),

M

△叱=(12-47)

△2

d\P,d\Pt

de

j二j斫(12-48)

廠-2M2

對(duì)于PV節(jié)點(diǎn),

△叱二(12-49)

AV.2

33M

Jdf

=s(12-50)

ij~cAV；2d\V^

百百

用牛頓一一拉夫遜法計(jì)算潮流的流程框圖示與圖12-5.首先要

輸入網(wǎng)絡(luò)的原始數(shù)據(jù)以及各節(jié)點(diǎn)的給定值并形成節(jié)點(diǎn)導(dǎo)納矩陣。輸入

節(jié)點(diǎn)電壓初值靖）和尸）），置迭代計(jì)數(shù)k=0。然后開(kāi)始進(jìn)入牛頓法的迭代

過(guò)程。在進(jìn)行第k+1次迭代時(shí)，其計(jì)算步驟如下：

（1）按上一次迭代算出的節(jié)點(diǎn)電壓值"和/的（當(dāng)k=0時(shí)即

為給定的初值），利用式（12-40）和（12-41）計(jì)算各

類節(jié)點(diǎn)的不平衡量△票二和△匕9）。

（2）按條件（12-37）校驗(yàn)收斂，即

max{|M"),2XQ?}<￡。

如果收斂，迭代到此結(jié)束，轉(zhuǎn)入計(jì)算各線路潮流和平衡節(jié)點(diǎn)的功

率，并打印輸出結(jié)果。不收斂則繼續(xù)計(jì)算。

（3）利用式（12-43）和（12-44）計(jì)算雅可比矩陣各元素。

（4）解修正方程（12-42）求節(jié)點(diǎn)電壓的修正量

（5）修正各節(jié)點(diǎn)的電壓，

6，向）=彳⑹+△4腦力"7=/⑹+M⑹

（6）迭代計(jì)數(shù)加1,返回第一步繼續(xù)迭代過(guò)程。

三、牛頓法潮流計(jì)算的極坐標(biāo)形式

當(dāng)節(jié)點(diǎn)電壓用極坐標(biāo)形式表示時(shí)，

甘=匕exp范=V（cos4+JsinR）o

節(jié)點(diǎn)功率方程（12-4）將寫(xiě)成：

6=匕￡V（G）cos%+B.sin%），

，：（12-53）

2=匕力匕（，sin%+與cos%）

式中，是i,j兩節(jié)點(diǎn)電壓的相角差。

方程式（12-53）把節(jié)點(diǎn)功率表示為節(jié)點(diǎn)電壓的幅值和相角的函數(shù)。

在有n個(gè)節(jié)點(diǎn)的系統(tǒng)中，假定第Pm號(hào)節(jié)點(diǎn)為PQ節(jié)點(diǎn)，第m+Pn-1

號(hào)節(jié)點(diǎn)為PV節(jié)點(diǎn)，第n號(hào)節(jié)點(diǎn)為平衡節(jié)點(diǎn)。匕和，是給定的，PV節(jié)

點(diǎn)的電壓幅值匕.「匕川也是給的。因此，只剩下n—l個(gè)節(jié)點(diǎn)的電壓相

角心,…&和m個(gè)節(jié)點(diǎn)的電壓幅值匕是未知量。

實(shí)際上，對(duì)于每一個(gè)PQ節(jié)點(diǎn)或每一個(gè)PV節(jié)點(diǎn)都可以列寫(xiě)一個(gè)有

功功率不平衡量方程式：

kPH=匕一匕￡V（Ggcos%+與sin琢）=0（i=1,2,…，〃-1）。（12-54）

J=i

而對(duì)于每一個(gè)PQ節(jié)點(diǎn)還可以再列寫(xiě)一個(gè)無(wú)功功率不平衡量方程式：

△Q,=Qi,~Qi=Qis-匕￡匕（Gjsin%+%cos%）=0（i=1,2,…,m）。（12-55）

j=i

式（12-54）和式（12-55）一共包含了nT+m個(gè)方程式，正好同未知

量的數(shù)目相等，而比直角坐標(biāo)形式的方程式少了一個(gè)n-l-m個(gè)。

對(duì)于方程式（12-54）和（12-55）可以寫(xiě)出修正方程式如下：

APHN

△QKLV￡AV

式中,

△a

△P=卜3=

地一】

心」(12-57)

AV.

AV

AV=2