某220kv電網(wǎng)潮流計算及輸電線路繼電保護配置

目錄

一、任務(wù)書......................................................3?4

二、網(wǎng)絡(luò)參數(shù)的計算..............................................5?9

三、電力系統(tǒng)潮流的計算..........................................10?*

附：PQ分解法潮流計算源程序及結(jié)果...........................18?30

四、短路電流的計算..............................................31?46

附：1、對稱短路程序.........................................36?37

2、不對稱短路程序.......................................43?46

五、線路保護的整定...............................................47?64

六、附錄：

a：對稱短路計算的系統(tǒng)接線圖.................................65

b：零序網(wǎng)絡(luò)圖...............................................66

c；潮流分布圖...............................................67

d：線路保護配置圖..........................................68

某220kv電網(wǎng)潮流計算及輸電線路繼電保護配置

(-)：已知(1)系統(tǒng)最大運行方式為四臺發(fā)電機滿發(fā)和系統(tǒng)投入運行;系統(tǒng)最小運行方式為停兩臺發(fā)電機(Fl,F3),

各負荷減半.

(2)系統(tǒng)各負荷及線路參數(shù)如圖所示，各變壓器及發(fā)電機型號分別為：

F1?F4:SF-100-40/8540Pc=1OOMWUC=13.8KVCOS^=0.9n=98.16%

Bl?B3:SFP7-15(X)(X)/220220±2x2.5%/13.8Ud=I4%YN,d11

B4:SFP7-150000/100110±2x2.5%/l3.8Ud=13%YN,d11

B5-B8:SFPSZ7-12(X)00/220220±8x1.25%/115/10.5YN,ynO,dll

UI2=23.5%U23=13.3UIJ=7.7%

B9:RT=3Q,Xr=110Q

(3)線路參數(shù)如圖中所示.

(-)：設(shè)計任務(wù)

(I)計算各元件阻抗標(biāo)么值(Sj=IO()MVA,Vj二VN),并畫出正序，負序，零序等效網(wǎng)絡(luò)圖；

(2)對系統(tǒng)進行潮流計算;(用C語言)

(3)對5、6、7、8點進行各種類型的短路電流計算;(QB)

(4)為線路6-7選擇保護方式；

(5)對所選保護進行整定計算；

(6)對保護進行選型；

(7)畫出保護原理圖；

(8)書寫設(shè)計說明書及準(zhǔn)備答辯.

(三)：設(shè)計成果

(1)系統(tǒng)潮流分布圖一張及短路電流表一份；

(2)線路保護配置圖一張；

(3)保護原理圖一張；

(4)設(shè)計說明書一份。

△Ps2=l/2（APs（l-2）+APS（2-3）-APS（3-I））

APS3=1/2（APS（2-3）+APS（3-1）-APs（I-2））

求出各繞組的短路損耗后，便可導(dǎo)出雙繞組變壓器計算電阻相同形式的算式，

232

KP：{Ri}Q={PSi}Kw{VN}KVX10/{SN}KVA

2.電抗X1,X2,X3：和雙繞組變壓器一樣，近似地認為電抗上的電壓降就等于短路電壓.在給出短路電壓力后,

與電阻的計算公式相似,各繞組的短路電壓為

Vsi%=l/2（VS（i-2）%+Vs（3-l）%-Vs（2-3）%）

VS2%=1/2（Vs（l-2）%+vS（2-3）%-VS（3-I）%）

Vs3%=1/2（VS（2-3）%+Vs（3-1）%-VS（1-2）%）

各繞組的等值電抗為：（乂山尸\^%乂{\^2}|＜丫'103/100/{51\1}1＜丫人

3.導(dǎo)納GTTBT和變比ki2、k13＞k23：計算與雙繞組變壓器相同.

三.線路參數(shù)的計算.

輸電線路的參數(shù)有四個:反映線路通過電流時產(chǎn)生有功功率損失效應(yīng)的電阻;反映截流導(dǎo)線周圍產(chǎn)生磁場效應(yīng)

的電感;反映線路帶電時絕緣介質(zhì)中產(chǎn)生泄漏電流及導(dǎo)線附近空氣游離而產(chǎn)生有功功率損失的電導(dǎo);反映帶電導(dǎo)線

周用電場效應(yīng)的電容.輸電線路的這些參數(shù)通常可以認為是沿全長均勻分布的,每單位長度的參數(shù)為r、x、g及b.

當(dāng)線路長為l（km）時，

R=rl;X=xl;G=gl;B=bl由于沿絕緣子的泄漏很小，可設(shè)G=0.

四.標(biāo)么值的折算.

建立電力網(wǎng)絡(luò)和電力系統(tǒng)的數(shù)學(xué)模型,需解決標(biāo)么值的折算問題.

進行電力系統(tǒng)計算時，除采用有單位的阻抗、導(dǎo)納、電壓、電流、功率等進行運算外,還可采用沒有單位的阻抗、

導(dǎo)綱、電壓、電流、功率等的相對值進行運算.前者稱有名制,后者稱標(biāo)么制.標(biāo)么制之所以能在相當(dāng)寬廣的范圍內(nèi)

取代有名制,是由于標(biāo)么制具有計算結(jié)果清晰、便于迅速判斷計算結(jié)果的正確性、可大量簡化計算等優(yōu)點。

標(biāo)么值=實際有名值（任意單位）/基準(zhǔn)值（與有名值同單位）

對于直接電氣聯(lián)系的網(wǎng)絡(luò),在制訂標(biāo)么值的等值電路時，各元件的參數(shù)必須按統(tǒng)一的基準(zhǔn)值進行歸算.由「各元件

的額定值可能不同，因此，必須把不同基準(zhǔn)值的標(biāo)么阻抗換算成統(tǒng)一基準(zhǔn)值的標(biāo)么值.現(xiàn)統(tǒng)一選定的基準(zhǔn)電壓和基

準(zhǔn)功率分別為V和S,則電抗的實際有名值換算成標(biāo)么值，即

在工程計算中規(guī)定,各個電壓等級都以其平均額定電壓V作為基準(zhǔn)電壓.根據(jù)我國現(xiàn)行的電壓等級，各級平均額定

電壓規(guī)定為3.15,6.3,10.5,15.75,37,115,230,345,525KV

五.具體計算及結(jié)果

根據(jù)原始數(shù)據(jù)及資料顯示，通過公式計算可得:（選擇基準(zhǔn)值:SB=100MVA,VB=Vav）

1、雙繞組變壓器:B1,B2,B3(型號SFP7----150(X)0/220)

(參數(shù)2^0=14()1<\乂2^*=4501<\￥25%=14%,10%=().8,丫1\1—dl1)

RT=APSXUN2X10?/SN2=450X2202X103/l500002=0.968Q

RT^RTXSB/VBM.968X100/2302=0.0018

2323

XT=Vx%xVNX10/1()()/SN=14X220X10/100/150000=45.1730

XT*=XTXSB/VB2=45.173X100/2302=0.0854

BT=IO%SNX10-3/100/VN2=0.8X150000x10-3/l00/2302=0.0227x103S

BT*=BTXVB2/SB=0.0227X10-3x2302/l00=0.012

B4(型號SFP7----150000/110)

(參數(shù)2^-1071<鞏2^\二5471<\￥、5%=13%,10%=06丫1\1—<111)

RT=547XH()2XIO3/150(X)02=0.2942Q

RT*=0.2942X100/1152=0.0022

Xi=13xll02xlOV100/150000=10.487。

XT-=10.487x100/1152=0.0793

Br=0.6x150000x103/100/1102=0.744x104S

2

BT*=0.744XIO^XII5/l00=0.0098

2、三繞組變壓器:B5?B8(型號SFPSZ7-120(X)0/220220±8x1.25%/!15/10.5

YN,ynO,dllt/12=23.5%=13.3%f/I3=7.7%APv=425▽4)=dO)

.%=-x(23.5+13.3-7.7)=14.55%

s}2

(/“1I1%=-2x(7.7+13.3-23.5)/=-1.25%=0

。仙%=lx(23.5+7.7-13.3)=8.95%

{RT)n={APs}KW{VN2}KVX103/(SN2}KVA=425X2302X103/l200002=1.5613

2

/?r.=1.5613x100/230=0.002951

2

14.552203

100X120000xlO=58.685Q

8.952202

Xni=-----x----------xlO3=36.098Q

100120000

58.685x100

X[*==O.llO9=X

2302op

36.098x100

=0.068=X。川.

2302

X[.+X][.=0.1109

〃242/121,

K=-----------=1

230/115

3、聯(lián)絡(luò)變壓器B9：=更吧=0.00567=0.2079K=231/110=0.957

23(f23O2230/115

4、線路L1?L7:(L1,L5,L6為雙回路,L2,L3,L4,L7為單回路)

L1：型號為2XLGJ-400/75km,參數(shù)為:r=0.08C/km,x=0.397Q/km,b=2.92X10-6s/km

Ri=l/2Lxr=75x0.08/2=3QXi=l/2Lx=0.397><75/2=14.888Q

Bi=2xLb=2x2.92x10-6x75=4.38x104S

RP=RXSB/VB2=3X]00/2302=0.00567XI*=XIXSB/VB2=14.888X100/2302=0.0281

Bi*=B?XVB2/SB=4.38X10-4x2302/l00=0.2317Bi“2=0.1158

L2：型號為LGJ-400/50km.參數(shù)為:r=0.08Q/km,x=0.397Q/km,b=2.92X10-6s/km

64

R2=Lxr=50x().08=4QX2=Lx=0.397x50=19.85QB2=Lb=2.92x1()-X5()=1.46x1()S

22

R2,=RXSB/VB2=4X100/115J00302X2*=X2XSB/VB=19.85X100/115=0.1501

BL=B，KV『/SB=1.46x104x1152/l00=0.0193B-/2=0.00965

L3:型號為LGJ-185/70km.參數(shù)為:r=0.17C/km,x=0.395Q/km,b=2.79X10-6s/kin

6

R3=Lxr=70x0.17=11.9QX3=Lx=0.395x70=27.65CB3=Lb=2.79x10x70=1.953x109S

222

R3*=RXSB/VB2=11.9x100/115=0.09X3?=X3XSB/VB=27.65X100/115=0.2091

242

B3*=B3XVB/SB=1.953X]0-X|15/l00=0.0258B3*/2=0.0129

L4:型號為LGJ-150/74km.參數(shù)為：r=0.21Q/km,x=0.403Q/km,b=2.74X10-6s/km

64

R4=Lxr=74x0.21=15.54。X4=Lx=0.403x74=29.822QB4=Lb=2.74x10-x74=2.03x10S

222

R4*=RXSB/VB2=15.54x100/115=0.1175X4*=X4XSB/VB=29.822X100/115=0.2255

242

B4*=B4XVB/SB=2.03X10-X115/1()0=0.0268B4*/2=0.0134

L5：型號為2XLGJ-120/70km.參數(shù)為：r=0.27Q/km,x=0.409C/km,b=2.69X10-6szkm

R5=l/2Lxr=70x0.27/2=9.45QX5=l/2Lx=0.409x70/2=14.315Q

64

B5=2xLb=2x2.69x10x70=3.766x10-S

J22

R5*=R5XSB/VB'9.45X100/115=0.0714X5*=X5XSB/VB=14.315X100/115=0.1082

242

B5*=B5XVB/SB=3.766X10-x115/100=0.0498B5*/2=0.0249

L6：型號為2XLGJTOO/73km.參數(shù)為:r=0.08Q/km,x=0.397Q/km,b=2.92X10-6s/km

R6=l/2Lxr=73x0.08/2=2.92QX6=l/2Lx=0.397x73/2=14.491Q

64

B6=2xLb=2x2.92x10-x73=4.26x10S

2222

R6*=ROXSB/VB=2.92X100/230=0.0055X6*=X6XSB/VB=14.491x100/230=0.0274

242

B6*=B6XVB/SB=4.26X10-x230/l00=0.2254B6*/2=0.1127

L7：型號為LGJ-400/135km.參數(shù)為:r=0.08Q/km,x=0.397Q/km,b=2.92X10-6s/km

R7=Lxr=135x0.08=10.8QX7=Lx=0.397x135=53.595QB7=Lb=2.92x1()%135=7.844x1(尸S

222

R7*=R7XSB/VB=10.8x100/2302=0.0204X7*=X7XSB/VB=53.595X100/230=0.1013

242

B7*=B7XVB/SB=7.844X10-x230/l00=0.4149B7*/2=0.20745

5、發(fā)電機FLF4：(型號：SF-100/40-854Pc=l()()MWcos。=0.9U,=13.8KV

X".=0.2034X(r=0.0927

6、PQ節(jié)點PQ標(biāo)么值:

4：S4=180+jl00S4*=S4/SB=1.8+jl.O

5：S5=72+j40S5*=S5/SB=0.72+j0.4

8：S8=120+j63S8*=S8/SB=1.2+j0.63

9：S9=36+j20S9*=S9/SB=0.36+J0.2

10：Sio=4O+j25Sio*=Sio/SB=O.4+jO.25

7、零序電抗標(biāo)么值：

X向LJJ=4.L5IXZ1=4.5x0.0281=0.12645

=3X「=3x0.1501=0.4503

Xg=3XA=3x0.2091=0.6273

XL=3X/4=3X0.2255=0.6765

=3X/s=3x().1082=0,3246

Xi=4.5X/6=4.5x0.0274=0.1233

乂…=4.5Xl=4.5x0.1013=0.45585

電力系統(tǒng)潮流計算

P-Q分解法潮流計算

主要步驟

1.導(dǎo)納矩陣的形成

2.因子表的形成（三角分解法）

3.給定電壓初值

4.計算不平衡功率△P*除以V2

5.計算不平衡功率AQ*除以丫產(chǎn)

6.判斷收斂性

7.回代，修正△盤k）,AV*

8.計算全線路功率

電力系統(tǒng)的潮流計算機算法

潮流計算的任務(wù)：對給定的運行條件確定運行狀態(tài)，如果各母線上的電壓（幅值及相角）網(wǎng)絡(luò)

中的功率分布以及功率損耗等，幾個節(jié)點電力系統(tǒng)的潮流方程的一般形式

根據(jù)電力系統(tǒng)的實際運行條件，按給定變量的不同，一般節(jié)點可分三種：

1.PQ節(jié)點：有功功率P和無功功率Q給定的，節(jié)點電壓（V,8）是待求量。

2.PV節(jié)點：有功功率P和電壓幅值V給定的，節(jié)點的無功功率Q電壓的相位角6是待求量。

3.平衡節(jié)點：網(wǎng)絡(luò)中至少有一個節(jié)點的有功功率P不能給定，此節(jié)點承擔(dān)系統(tǒng)的有功功率平

衡。

網(wǎng)絡(luò)方程是線性方程，但由于在定解條件中不能給定節(jié)點電流，只能給出節(jié)點功率，從而使潮

流方程變?yōu)榉蔷€性方程，由于平衡節(jié)點的電壓已經(jīng)給定，假設(shè)系統(tǒng)中有n個節(jié)點，其中有m個PQ

節(jié)點，n-（m+l）個PV節(jié)點和一個平衡節(jié)點，平衡節(jié)點不參加求解。

形成導(dǎo)納矩陣的方法及框圖

1.形成導(dǎo)納矩陣的程序框圖

為了形成導(dǎo)納矩陣，必須個計算機輸入必要的原始數(shù)據(jù)，條支路的原始數(shù)據(jù)應(yīng)包括兩端的節(jié)點號和支路阻

抗，對變壓器支路的原始數(shù)據(jù)應(yīng)包括他的變比，對于電力線路還應(yīng)包括它的容納，即：I,J,R,X,K（-B/2）e

對以上數(shù)據(jù)說明如下：

（1）當(dāng)支路為接地支路時，規(guī)定節(jié)點I處置零，J處填寫接地支路所在的節(jié)點號，支路參數(shù)用阻抗填寫,

對于K處置零，但必須填寫。

（2）對于純阻抗支路，K處置零。

（3）對于具有容納的線路，電容電納不作單獨支路處理，而且把總?cè)菁{的一半負值填寫在K處，其符

號用以區(qū)別支路的性質(zhì)。

（4）對于變壓器支路，采用下圖所示等值電路，即非標(biāo)準(zhǔn)變比在J側(cè)，變壓器阻抗在I側(cè)，K處填寫

實際變化

YB/K

*

ic———----------13-----。J

工上j

□K2

ZB

—0J

-------(Uh-—

2.框圖

注：Bi不記充電電容和非標(biāo)準(zhǔn)變比時，導(dǎo)納矩陣的虛部

導(dǎo)納矩陣的虛舒

B2

二、形成因子表的方法及框圖

1.用行消去過程形成因子表的程序框圖

2.因子表的形成（三角分解法）

消去法求非線性方程組的一種常用算法是對方程式的系數(shù)矩陣A進行三角分解，在本次電力系統(tǒng)潮流計算中

采用的三角分解是將非奇方陣A分解為單位下三角矩陣L和上三角矩陣R的乘積。

A=LR,非奇方陣A被表示為矩陣L和R的乘積：這兩個三角矩陣稱為A的因子矩陣，兩個因子矩陣的元

素計算公式：

4-1

0=2,3,???n、

j=i,2,…r

將人=1^代入線性方程組，便得LRX=B,這個方程又可以分解為以下兩個方程:

①LF=B②RX=F

先由方程組①自上而下地依次算出fl、f2、fn，其計算通式為

Z=2一乞Lufj(，=1,2,???n)

7=1

這一步演算相當(dāng)于消元過程中對原方程式右端常數(shù)向量所作的變換，只順用到下三角因子矩陣。方程組②的

求解屬于回代過程，只順用到上三角因子矩陣以及經(jīng)過消元變浜的右端常數(shù)向量，方程組可以自下而上地逐步

算出待求量，其計算通式為：

一之芍…M)

j=M

X.i=n,n-1,??-1)

(Z-1)

三、節(jié)點電壓的表示：

1.極坐標(biāo)表示法：

匕=匕/&=K(sing.)

節(jié)點功率方程表示為：

「〃

e=K￡V/g.cos4+B)sin%)

Y汩

①

?

、。,=匕￡匕(GysinM-Bi.cos^)

j=i

M=5,一S,.(兩節(jié)點電壓的相位角)

方程式①把節(jié)點功率表示為節(jié)點電壓的幅值和相角的函數(shù)

在有n個節(jié)點的系統(tǒng)中，假定第1?m號節(jié)點為PQ節(jié)點，第m+l~n-l號節(jié)點為平衡節(jié)點。Vn和50是給定

的，PV節(jié)點的電壓幅值Vm+i~Vz也是給定的。因此，只剩下n-1個節(jié)點的電壓相角火…61H和m個節(jié)點的電

壓幅值V]…Vm是未知的量。

對于每一個PQ節(jié)點或每一個PV節(jié)點都可以列寫一個有功功率不平衡方程：

n

△4=4一片=4-匕2匕(％35%+82嗝)=0(i=l,2,??.n-l)對于每

j=i

一個PQ節(jié)點可以列寫一個無功功率不平衡方程式：

△0=Qis-Qi=2s-匕2匕(G“sin%-%cos%)=0(i=l,2,…m)

y=i

所以可以寫出修正方程式

②

A。二

其中②

AV匕

AV=252

匕

四、P?Q分解法

1.在交流高壓電網(wǎng)中，輸電線路的電抗要比電阻大得多，系統(tǒng)中母線有功功率的變化主要受電壓相位的影

響，無功功率的變化規(guī)則主要受母線電壓幅值變化的影響，在修正方程式的系數(shù)矩陣中，偏導(dǎo)數(shù)立竺和婆

dV

的數(shù)值相當(dāng)于偏導(dǎo)數(shù)之絲和2半是相當(dāng)小的，作為簡化第一步，可將方程②中N、k略去不計，既認為它

65dV

們的元素都等于零，這樣，便可分解為n-1階和m階的兩個方程：

③"=-HN④\Q=-LV^AV

所謂P-Q分解法，節(jié)點有功功率不平衡量又用于修正電壓的相位。節(jié)點的無功功率不平衡量值用于修正電

壓幅值，方程③、④分別輪流迭代。

2.H、L元素都是節(jié)點電壓幅值和相角差函數(shù)，其數(shù)值在迭代過程中是不斷變化的，從而，最關(guān)鍵的一

步是把系數(shù)矩陣H、L簡化為常數(shù)矩陣。

方法：在一般情況下，在線路兩端電壓的相角差是不大的（不超過0。?20。），因此，以為

cos4=1,Gijsg的WBij,此外，與系統(tǒng)各節(jié)點無功功率柜適應(yīng)的導(dǎo)納Bui必須小于該節(jié)點自導(dǎo)納虛部，即

%*=3也或2工匕％

??.矩陣H、L可簡化為⑤H=B'Vj?L=VD2-B'-VD2將⑤⑥分別帶入③④

得:V-'AP=-B'V/y

r——、＜,'

△Q=-vB^^vV-\\Q=-BnAV

D2X.

■雙-

B”B]2B|3…Bg

v,1

…B2,n-lV2Ab2

△匕

?????

(1)—??????*

v??????

2k-1

B〃_]2…

…Bun-l.n-\_AM.

Mi

_K..J

-AV/

■AS8B

HnBI3…%一

匕

B……BmAV

2222

—

△02??????

(2)?????*

匕??????

BB...BAV

\_Dm\Dm2…mmJ

△Q,“

利用計算節(jié)點功率的不平衡量，用修正方程⑴⑵解出修正量A3及AV,并換下述條件:

max|）|}<8Pmaxb。，）}<氣

校驗收斂.這就完成分解法的計算了.

3.P-Q分解法計算潮流的程序框圖

PQ分解法源程序

include"maih.h"

#inckide"stdio.h"

#include"stdlib.h"

man()

(

inta=l.k.i.j,d=l.m.n.ls,nl,kk;

floatz[100|[5].pqf!001[3],p.qxx,kb.b.g:

floatyg[100][100],yb[100][100],bri00][100],b2[100][100],v[100],w[10(?];

floatpp[l001.qq[100].e[l00].fll001.ggn00][2].tt[100][2j,pw|10()].qw[100].wwll00];

floatt,t1,tp,tq,p1,p2,qI,q2,t2;

printf("\n");

prin(f("電力系統(tǒng)潮流計算\n")；

printfCAn");

printfC指導(dǎo)老師：設(shè)計：2002年12月\n");

printfCAn");

printf("\n");

printfC注：各參數(shù)用標(biāo)么值表示\n")；

printf("\n");

begii).piinlfT請選擇.1-建立電網(wǎng),2-直看數(shù)據(jù),3-修改數(shù)據(jù),4-運行計算,5-退出程序

if(a==l)

{

gotocrcat;

)

elseif(a==2)

(

gotoview;

}

elseif(a==3)

(

gotomodify;

}

elseif(a==4)

I

gotorun;

}

else

(

exit(O);

}

creat:printf("請輸入系統(tǒng)的節(jié)點數(shù)N:");scanf("%d",&n);

printf("請輸入系統(tǒng)的支路數(shù)LS:");scanf("%d",&ls);

printf("請輸入系統(tǒng)的PQ節(jié)點數(shù)M:");scanf("%d",&m);

for(k=l;k<=ls;k++)

printf("請輸入支路的%d的數(shù)據(jù)1(接地支路為0):",k);scanf("%d'\&i);

prinlf("請輸入支路的%d的數(shù)據(jù)J:",k);scanf(',%d,,,&j);

printfC請輸入支路的%d的數(shù)據(jù)R:".k):scanf("%f,.&r);

prinlf("請輸入支路的%d的數(shù)據(jù)X:".k);scanf(M%f;&x);

prinlf("請輸入支路的%d的數(shù)據(jù)KB:';k);scanf("%f,.&kb);

z[k][l]=i;z[k][2]=j;z[k][3]=r;z[kH4]=x;z[k][5]=kb;

I

ibr(k=l：k<=m;k++)

(

pq[kl[I]=k;

prinlf("請輸入PQ節(jié)點%d的注入有功P:",k);scanf("%f',&p);

printf("請輸入PQ節(jié)點泡的注入無功Q:';k);scanf("%f;&q);

pq[k][2]=p;pq[k][3]=q;

}

for(k=m+1;k<=n-l;k++)

(

pq[k][l]=k;

prinlfT請輸入PV節(jié)點％d的注入有功P:",k);scanf("％r',&p);

prinlf("請輸入PV節(jié)點％d的電壓模二k);scanf("%T,&q);

pq[k][2J=p;pq(k)[3]f;

I

pq(k][l]=n；

printf("請輸入平衡節(jié)點％d的電壓模二n);scanf("%F.&p);

prinlf("請輸入平衡節(jié)點颯的電壓角:”,n);scanf("%F,&q);

pq[n][2]=p;pq[n][3]=q;

gotobegin:

view:printf("\n");prin(f("N=%d'dLS=%d\tM=%d",n,ls,ni);prinlf(',\n");

printfC支路數(shù)據(jù)”);printfC\ir);

printt("\t節(jié)點l\t節(jié)點八電阻V電抗\l變比或電納");prinlf("\n");

fbr(i=l;i<=ls;i++)

(

printfC支路％d",i);

for(j=l；j<=5:j++)

{

if(j<=2)

(

printf("\t%1.0r,z[iJ[j])；

)

else

(

printf("%14.6r,z(i]g]);

}

)

printf("\n");

printf("\tPQ節(jié)點\l注入有功V注入無功)prinlf("\n");

Ibr(i=l;i<=m;i++)

for(j=1;j<=3;j++)

(

if(j<=D

(

prinlfC\t%1.0f'.pq「皿)；

)

else

(

printf("V%14.6r',pq[i][j]);

)

)

printf("\n");

)

printf("\tPV節(jié)點\l注入有辦l電壓模)prinlf("\n");

fbr(i=rn+1;i<=n-l;i++)

{

for(j=lj<=3J++)

{

if(j<=D

{

pr!ntfC'\t%I.Or.pq[iHj]);

)

else

{

printf("\t%14.6F.pq[i皿)；

}

)

printf("\n");

}

printf("\n");

printf("\t平衡節(jié)點\l電壓模\l電壓角)prinlf("\n")；

for0=l；j<=3;j++)

I

if(J<=D

{

pr!ntfC'\t%1.0r.pq[n][j]);

}

else

{

printf("\t%14.6f'.pqln]Ul)；

)

printf("\n");

gotobegin;

[11。由1丫邛面可”1-修改支路數(shù)據(jù);2-修改節(jié)點注入功率;3-返回主菜單:");2】值"％(1”.&(1);

if(d==l)

(

gotom1;

}

elseif(d==2)

(

gotom2:

}

elseif(d==3)

{

gotobegin;

I

m1:printf("請輸入要修改的支路號:");scanf("%d",&k);

printf("請輸入支路的％d的新數(shù)據(jù)I:",k);scanf("%d",&i);

pinifT請輸入支路的%d的新數(shù)捷J.",k),scanf("%d",&j),

printf("請輸入支路的％d的新數(shù)據(jù)R：k);scanf("%f',&r);

prinlf(”請輸入支路的％d的新數(shù)據(jù)X:",k);scanf("%f,,&x);

printf("請輸入支路的％d的新數(shù)扼KB:",k);scanf("%r,&kb);

z[kHI]=i;z[k][2]=j;z[k][31=r;z[k][4；=x;z[k][5]=kb;

gotomodify;

m2:k=0;

printf("請輸入需修改的節(jié)點號:");$canf("%d",&k);pq[k][l]=k;

prinlf("請重新輸入節(jié)點的注入有功或電壓參數(shù):");scanf("%口&pq[k]⑵);

printf("請重新輸入節(jié)點的注入無功或電壓參數(shù)巧;scanf("%匚&四出]閭)；

gotomodify;

run:for(i=l;i<=n:i++)

(

fbr(j=l;j<=n;j++)

{

yb[iJUl=o；

}

}

/*形成Bl,不計線路電容和變壓器非標(biāo)準(zhǔn)變比*/

for(k=l;k<=ls;k++)

(

i=z[k][l]U=z[k][2];r=z[k][31；x=z[k](4];

b=-x/(r*r+x*x);

if(i==0)

I

gotopO;

)

yb(il[i]=ybIi][i)+b;yb|i]Ul=yt[i|[jl-b:yb[j||il=ybniUl；

pO：yb(jl(jl=ybrjlfjl+b;

|

for(i=l;i<n;i++)

(

for(j=l;j<=n;j++)

(

bl[illjl=yb[i]Ul；

}

}

/*形成導(dǎo)綱矩陣*/

for(i=l;i<=n;i++)

(

for(j=l；j<=n;j++)

{

yg[i]U]=o；yb[i]U]=o；

I

}

for(k=l;k<=ls;k++)

(

i=z[k][l]戶z[kH2J;尸z[k][3];x=z岡⑷;kb=z[k][5];

b=-x/(r*r+x*x);

g=r/(r*r+x*x);

if(i==O)

(

gotopl;

}

it(kb>0)

(

gotop2;

}

yb[i][i]=yb[i][i]-kb;

ybU]U]=yb[j皿-kb;

p=l;q=l;gotop3;

pi:p=l;gotop4;

p2:p=kb*kb;q=kb;

p3：yg[il[i]=yg[i](il+g;ybni[i]=yb[i][il+b；

yg[i]U]=yg[i]lj]-gzq^b[i]|j]=yb[i]U)-b/q；

p4：yg[j][j]=yg[j][j]+g/p;yb[j][j]=ybU][j]+b/p；

}

/*形成b2*/

for(i=l;i<=m;i++)

fbr(j=kj<=m;j++)

b2[iHj]=ybmij];

}

prinif("矩陣B:\n");

for(i=l;i<=n;i++)

(

for(j=l；j<=n；j++)

{

printf("%14.7f',yb[i][j])；

}

printfCNn");

}

prinif("矩陣G:\n");

for(i=l;i<=n;i++)

(

for(j=l;j<=n;j++)

I

piintf("%14.7f',yg[i]U])；

}

printf("\ii");

}

/*形成bl因子表*/

for(i=2;i<=n-l;i++)

(

for(k=1;k<=i-l;k++)

bHiJLkJ=-blliJlkJ/bHkJlkj;

for(j=k+ly<=n-ly++)

{

bl[i]U)=bl[i]U]+bl[i][k]*bl[k][j];

)

}

}

printfC'Bl因子表:\n");

for(i=l;i<=n-l;i++)

(

for(j=l;j<=n-l;j++)

(

prinlf("%14.7f',bl[i]UD;

I

printfCVn");

}

/*形成b2因子表:*/

for(i=2;i<=m;i++)

for(k=l;k<=i-l;k++)

b2[i][k]=-b2[i][k]/b2[k][k];

fbr(j=k+1;j<=m:j++)

(

b2[iJU]=b2[i]UJ+b2[i]lk]*b2[k]UJ;

)

}

}

printf("B2因子表:\n");

for(i=l;i<=m;i++)

(

ibr(j=l;j<=m;j++)

{

printf("%14.7f',b2[i]U]);

I

Printf("\n")；

}

printf("\n");

/*給定電壓初值*/

for(i=l;i<=m;i++)

(

v(i]=l;

w[i]=O;

)

for(i=m+1;i<=n-l;i++)

(

v[i]=pq[i]⑶；

w[i]=O：

)

v[n]=pq[n][2];

w[n]=pq[n][3];

kk=O;

/*計算n-1個節(jié)點不平衡有功并除其電壓*/

al:for(i=l;i<=n-l;i++)

(

l=0;

for(j=1;j<=n;j++)

(

if(yb[i]UJ!=0)

t=t+v[j]*(yg[i]|j]*ccs(w[i]-w(j])+yb[i][j]*sin(w[i]-w[j]));

}

pp|i]=pq[i][2]/v[i]-t;

)

/*計算m個節(jié)點不平衡無功并除其電壓*/

for(i=l;i<=m;i++)

t=0;

for(j=l;j<=n:j++)

(

if(yb[iJ[jJ!=O)

t=t+v[jl*(yg[i]|j]*sin(w[il-wUl)-yb[iHjl*cos(w[i]-w|j]));

)

qq[i]=pq[i][3]/v[i]-t；

)

for(i=l;i<=n-l;i++)

(

if(fabs(pp[i])>0.00001)

(

gotoa4;

}

)

for(i=l;i<=m;i++)

{

if(fabs(qq|i])>0.00001)

(

gotoa2;

}

)

gotoa3;

產(chǎn)利用Bl因子表對右端項變換*/

a4:for(i=l；i<=n-2;i++)

(

t=0；

for(j=l;j<=ij++)

(

t=t+bl[i+l]g]?ppUl；

}

pp[i+l]=pp[i+l]+t;

)

nl=n-l;

/*回代求解*/

pp[nl]=pp[nl]/bl[nl][nl];

for(k=l;k<=nl-l;k++)

(

t=0:

for(j=nl-k+1:j<=nly++)

t=t+bl[nl-k]|j]*pp|j];

1

pp|nl-kl=(pp|nl-k]-t)/bHnl-k][nl-k];

)

/*修正電壓角*/

for(i=l;i<=nl;i++)

(

w(i]=w[i)-pp[i];

)

/*利用b2因子表對右端頂變換為

a2:fbr(i=l;i<=m-l;i++)

(

t=0:

for(j=l;j<=ij++)

(

t=t+b2[i+ll|jl*qqlj]；

)

qq[i+l]=qq[i+l]+t;

)

qq[m]=qq[m]/b2[m]|m];

蘆回代求解*/

for(k=l;k<=m-l;k++)

{

1=0;

for(j=n-k+1:j<=m:j++)

(

t=t+qq[j]*b2[m-kHj];

qqlm-k]=(qqlm-k]-t)/b2lni-kj[ni-kj;

)

/*修正電壓模*/

for(i=l：i<=ni:i++)

(

)

kk=kk++:

gotoa1;

a3:fbr(i=l；i<=n;i++)

{

e(i]=v[i]*cos(w[i]);

f[i]=v[i]*sin(w[i]);

ww[i]=w(i]*180/3.1415926;

)

/*計算支路功率*/

for(k=l;k<=ls;k++)

i=zfklfl];j=zfk](21;kb=zfk][51;

tl=-yg[i]|j]*(e[i]-e|j])+yb[i]U]*[f[i]-fU]);

t2=-yb[i][j]*(e[i]-e[j])-yg[i][j]*(f[i]-flj])：

節(jié)點間有功的部分Vi(Ui-Uj)Yij*/

ql=f[i]*tl-e[i]*t2;

p2=e|j]*(-tl)+f|j]*(-t2);

q2=f[j]*(-tl)-e|jl*(-t2);

if(kb>0)

I