某220kv电网潮流计算及输电线路继电保护配置

目录

一、任务书......................................................3〜4

二、网络参数的计算..............................................5〜9

三、电力系统潮流的计算..........................................10〜*

附：PQ分解法潮流计算源程序及结果...........................18〜30

四、短路电流的计算..............................................31〜46

附：1、对称短路程序.........................................36〜37

2、不对称短路程序.......................................43〜46

五、线路保护的整定...............................................47〜64

六、附录：

a：对称短路计算的系统接线图.................................65

b：零序网络图...............................................66

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

d：线路保护配置图..........................................68

某220kv电网潮流计算及输电线路继电保护配置

(-)：已知(1)系统最大运行方式为四台发电机满发和系统投入运行;系统最小运行方式为停两台发电机(Fl,F3),

各负荷减半.

(2)系统各负荷及线路参数如图所示，各变压器及发电机型号分别为：

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)线路参数如图中所示.

(-)：设计任务

(I)计算各元件阻抗标么值(Sj=IO()MVA,Vj二VN),并画出正序，负序，零序等效网络图；

(2)对系统进行潮流计算;(用C语言)

(3)对5、6、7、8点进行各种类型的短路电流计算;(QB)

(4)为线路6-7选择保护方式；

(5)对所选保护进行整定计算；

(6)对保护进行选型；

(7)画出保护原理图；

(8)书写设计说明书及准备答辩.

(三)：设计成果

(1)系统潮流分布图一张及短路电流表一份；

(2)线路保护配置图一张；

(3)保护原理图一张；

(4)设计说明书一份。

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

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

求出各绕组的短路损耗后，便可导出双绕组变压器计算电阻相同形式的算式，

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.导纳GTTBT和变比ki2、k13＞k23：计算与双绕组变压器相同.

三.线路参数的计算.

输电线路的参数有四个:反映线路通过电流时产生有功功率损失效应的电阻;反映截流导线周围产生磁场效应

的电感;反映线路带电时绝缘介质中产生泄漏电流及导线附近空气游离而产生有功功率损失的电导;反映带电导线

周用电场效应的电容.输电线路的这些参数通常可以认为是沿全长均匀分布的,每单位长度的参数为r、x、g及b.

当线路长为l（km）时，

R=rl;X=xl;G=gl;B=bl由于沿绝缘子的泄漏很小，可设G=0.

四.标么值的折算.

建立电力网络和电力系统的数学模型,需解决标么值的折算问题.

进行电力系统计算时，除采用有单位的阻抗、导纳、电压、电流、功率等进行运算外,还可采用没有单位的阻抗、

导纲、电压、电流、功率等的相对值进行运算.前者称有名制,后者称标么制.标么制之所以能在相当宽广的范围内

取代有名制,是由于标么制具有计算结果清晰、便于迅速判断计算结果的正确性、可大量简化计算等优点。

标么值=实际有名值（任意单位）/基准值（与有名值同单位）

对于直接电气联系的网络,在制订标么值的等值电路时，各元件的参数必须按统一的基准值进行归算.由「各元件

的额定值可能不同，因此，必须把不同基准值的标么阻抗换算成统一基准值的标么值.现统一选定的基准电压和基

准功率分别为V和S,则电抗的实际有名值换算成标么值，即

在工程计算中规定,各个电压等级都以其平均额定电压V作为基准电压.根据我国现行的电压等级，各级平均额定

电压规定为3.15,6.3,10.5,15.75,37,115,230,345,525KV

五.具体计算及结果

根据原始数据及资料显示，通过公式计算可得:（选择基准值:SB=100MVA,VB=Vav）

1、双绕组变压器:B1,B2,B3(型号SFP7----150(X)0/220)

(参数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)

(参数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、联络变压器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,参数为: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.参数为: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.参数为: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.参数为：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.参数为：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.参数为: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.参数为: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、发电机FLF4：(型号：SF-100/40-854Pc=l()()MWcos。=0.9U,=13.8KV

X".=0.2034X(r=0.0927

6、PQ节点PQ标么值:

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、零序电抗标么值：

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

电力系统潮流计算

P-Q分解法潮流计算

主要步骤

1.导纳矩阵的形成

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

3.给定电压初值

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

5.计算不平衡功率AQ*除以丫产

6.判断收敛性

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

8.计算全线路功率

电力系统的潮流计算机算法

潮流计算的任务：对给定的运行条件确定运行状态，如果各母线上的电压（幅值及相角）网络

中的功率分布以及功率损耗等，几个节点电力系统的潮流方程的一般形式

根据电力系统的实际运行条件，按给定变量的不同，一般节点可分三种：

1.PQ节点：有功功率P和无功功率Q给定的，节点电压（V,8）是待求量。

2.PV节点：有功功率P和电压幅值V给定的，节点的无功功率Q电压的相位角6是待求量。

3.平衡节点：网络中至少有一个节点的有功功率P不能给定，此节点承担系统的有功功率平

衡。

网络方程是线性方程，但由于在定解条件中不能给定节点电流，只能给出节点功率，从而使潮

流方程变为非线性方程，由于平衡节点的电压已经给定，假设系统中有n个节点，其中有m个PQ

节点，n-（m+l）个PV节点和一个平衡节点，平衡节点不参加求解。

形成导纳矩阵的方法及框图

1.形成导纳矩阵的程序框图

为了形成导纳矩阵，必须个计算机输入必要的原始数据，条支路的原始数据应包括两端的节点号和支路阻

抗，对变压器支路的原始数据应包括他的变比，对于电力线路还应包括它的容纳，即：I,J,R,X,K（-B/2）e

对以上数据说明如下：

（1）当支路为接地支路时，规定节点I处置零，J处填写接地支路所在的节点号，支路参数用阻抗填写,

对于K处置零，但必须填写。

（2）对于纯阻抗支路，K处置零。

（3）对于具有容纳的线路，电容电纳不作单独支路处理，而且把总容纳的一半负值填写在K处，其符

号用以区别支路的性质。

（4）对于变压器支路，采用下图所示等值电路，即非标准变比在J侧，变压器阻抗在I侧，K处填写

实际变化

YB/K

*

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

工上j

□K2

ZB

—0J

-------(Uh-—

2.框图

注：Bi不记充电电容和非标准变比时，导纳矩阵的虚部

导纳矩阵的虚舒

B2

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

1.用行消去过程形成因子表的程序框图

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

消去法求非线性方程组的一种常用算法是对方程式的系数矩阵A进行三角分解，在本次电力系统潮流计算中

采用的三角分解是将非奇方阵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

这一步演算相当于消元过程中对原方程式右端常数向量所作的变换，只顺用到下三角因子矩阵。方程组②的

求解属于回代过程，只顺用到上三角因子矩阵以及经过消元变浜的右端常数向量，方程组可以自下而上地逐步

算出待求量，其计算通式为：

一之芍…M)

j=M

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

(Z-1)

三、节点电压的表示：

1.极坐标表示法：

匕=匕/&=K(sing.)

节点功率方程表示为：

「〃

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

Y汩

①

・

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

j=i

M=5,一S,.(两节点电压的相位角)

方程式①把节点功率表示为节点电压的幅值和相角的函数

在有n个节点的系统中，假定第1〜m号节点为PQ节点，第m+l~n-l号节点为平衡节点。Vn和50是给定

的，PV节点的电压幅值Vm+i~Vz也是给定的。因此，只剩下n-1个节点的电压相角火…61H和m个节点的电

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

对于每一个PQ节点或每一个PV节点都可以列写一个有功功率不平衡方程：

n

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

j=i

一个PQ节点可以列写一个无功功率不平衡方程式：

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

y=i

所以可以写出修正方程式

②

A。二

其中②

AV匕

AV=252

匕

四、P・Q分解法

1.在交流高压电网中，输电线路的电抗要比电阻大得多，系统中母线有功功率的变化主要受电压相位的影

响，无功功率的变化规则主要受母线电压幅值变化的影响，在修正方程式的系数矩阵中，偏导数立竺和婆

dV

的数值相当于偏导数之丝和2半是相当小的，作为简化第一步，可将方程②中N、k略去不计，既认为它

65dV

们的元素都等于零，这样，便可分解为n-1阶和m阶的两个方程：

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

所谓P-Q分解法，节点有功功率不平衡量又用于修正电压的相位。节点的无功功率不平衡量值用于修正电

压幅值，方程③、④分别轮流迭代。

2.H、L元素都是节点电压幅值和相角差函数，其数值在迭代过程中是不断变化的，从而，最关键的一

步是把系数矩阵H、L简化为常数矩阵。

方法：在一般情况下，在线路两端电压的相角差是不大的（不超过0。〜20。），因此，以为

cos4=1,Gijsg的WBij,此外，与系统各节点无功功率柜适应的导纳Bui必须小于该节点自导纳虚部，即

%*=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,“

利用计算节点功率的不平衡量，用修正方程⑴⑵解出修正量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("电力系统潮流计算\n")；

printfCAn");

printfC指导老师：设计：2002年12月\n");

printfCAn");

printf("\n");

printfC注：各参数用标么值表示\n")；

printf("\n");

begii).piinlfT请选择.1-建立电网,2-直看数据,3-修改数据,4-运行计算,5-退出程序

if(a==l)

{

gotocrcat;

)

elseif(a==2)

(

gotoview;

}

elseif(a==3)

(

gotomodify;

}

elseif(a==4)

I

gotorun;

}

else

(

exit(O);

}

creat:printf("请输入系统的节点数N:");scanf("%d",&n);

printf("请输入系统的支路数LS:");scanf("%d",&ls);

printf("请输入系统的PQ节点数M:");scanf("%d",&m);

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

printf("请输入支路的%d的数据1(接地支路为0):",k);scanf("%d'\&i);

prinlf("请输入支路的%d的数据J:",k);scanf(',%d,,,&j);

printfC请输入支路的%d的数据R:".k):scanf("%f,.&r);

prinlf("请输入支路的%d的数据X:".k);scanf(M%f;&x);

prinlf("请输入支路的%d的数据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节点%d的注入有功P:",k);scanf("%f',&p);

printf("请输入PQ节点泡的注入无功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节点％d的注入有功P:",k);scanf("％r',&p);

prinlf("请输入PV节点％d的电压模二k);scanf("%T,&q);

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

I

pq(k][l]=n；

printf("请输入平衡节点％d的电压模二n);scanf("%F.&p);

prinlf("请输入平衡节点飒的电压角:”,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支路数据”);printfC\ir);

printt("\t节点l\t节点八电阻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节点\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节点\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平衡节点\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-修改支路数据;2-修改节点注入功率;3-返回主菜单:");2】值"％(1”.&(1);

if(d==l)

(

gotom1;

}

elseif(d==2)

(

gotom2:

}

elseif(d==3)

{

gotobegin;

I

m1:printf("请输入要修改的支路号:");scanf("%d",&k);

printf("请输入支路的％d的新数据I:",k);scanf("%d",&i);

pinifT请输入支路的%d的新数捷J.",k),scanf("%d",&j),

printf("请输入支路的％d的新数据R：k);scanf("%f',&r);

prinlf(”请输入支路的％d的新数据X:",k);scanf("%f,,&x);

printf("请输入支路的％d的新数扼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("请输入需修改的节点号:");$canf("%d",&k);pq[k][l]=k;

prinlf("请重新输入节点的注入有功或电压参数:");scanf("%口&pq[k]⑵);

printf("请重新输入节点的注入无功或电压参数巧;scanf("%匚&四出]闾)；

gotomodify;

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

(

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

{

yb[iJUl=o；

}

}

/*形成Bl,不计线路电容和变压器非标准变比*/

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；

}

}

/*形成导纲矩阵*/

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个节点不平衡有功并除其电压*/

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个节点不平衡无功并除其电压*/

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;

产利用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])：

节点间有功的部分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