电力系统毕业设计220kv电网潮流计算及输电线路继电保护配置

电力系统及其自动化专业设计人：学号： 指导老师目录一、 任务书 34二、 网络参数的计算 59三、 电力系统潮流的计算 1017附： PQ 分解法潮流计算源程序及结果 .18? 30四、 短路电流的计算 31 ? 46附：1、对称短路程序 36? 372 、不对称短路程序 4346五、 线路保护的整定 .47? 64六、 附录：a ：对称短路计算的系统接线图 .65b: 零序网络图 66c ：潮流分布图 67d: 线路保护配置图 68某220kv电网潮流计算及输电线路继电保护配置H4nuSMVA屈飞OKV)TnTB9STH3Q-nuH5KVMJ J4W、名 Ml(F1,F3),(一) : 已知 ( 1)系统最大运行方式为四台发电机满发和系统投入运行； 系统最小运行方式为停两台发电机各负荷减半?(2) 系统各负荷及线路参数如图所示， 各变压器及发电机型号分别为 ：F1F4:SF-100-40/8540P e=100MWU e=13.8KV COS P =0.9 n =98.16%B1B3:SFP7-150000/220220 _ 2 2.5% /13.8 Ud=14%YN,d11B4:SFP7-150000/100 110 _ 2 2.5%/13.8 Ud=13%YN,d11B5B8:SFPSZ7-120000/220220 8 1.25% /115/10.5 YN,yn0,d11U12=23.5%U 23=13.3U 13=7.7%b9: R = 3 1 , Xt =110 1(3) 线路参数如图中所示 .(二) : 设计任务(1) 计算各元件阻抗标么值(Sj=100MVA,V j=VN) ,并画出正序， 负序 ， 零序等效网络图 ；(2) 对系统进行潮流计算;(用 C 语言 )(3) 对 5 、 6、 7 、 8 点进行各种类型的短路电流计算； (QB)(4) 为线路 6-7 选择保护方式 ；(5) 对所选保护进行整定计算；(6) 对保护进行选型；(7) 画出保护原理图；(8) 书写设计说明书及准备答辩 .(三) ：设计成果(1) 系统潮流分布图一张及短路电流表一份；(2) 线路保护配置图一张；(3) 保护原理图一张；(4) 设计说明书一份。网络参数计算一、双绕组变压器的参数计算变压器的参数一般是指其等值电路中的电阻Rt,电抗电导Gt和电纳Bt ,变压器的变比Ko根据铭牌上所给的短路损耗Ps,短路电压Vs%,空载损耗Po,空载电流Io%。前两个数据由短路试验得到，用以确定Rt和Xt;后两个数据由空载试验得到，用以确定Gt和Bt。1. 电阻RT ： 变压器作短路试验时 ， 将一侧绕组短接， 在另一侧绕组施加电压， 使短路绕组的电流达到额定值?由于此时外加电压较小，相应的铁耗也小，可以认为短路损耗即等于变压器通过额定电流时原、昌U方绕组电阻的总损耗?在电力系统计算中， 常用变压器三相额定容量和额定线电压进行参数计算 ， 则公式232为:R T Q= APS KWV N KV X10 /S N KVA2. 电抗XT: 当变压器通过额定电流时， 在电抗上产生的电压降的大小 ， 可以用额定电压的百分数表示对于大容量变压器， 其绕组电阻比电抗小得多 ， 则公式 ：X t沪 Vs% X Vn2kv X 103/100/S N KVA3. 电导Gt：变压器的电导是用来表示铁芯损耗的.由于空载电流相对额定电流来说是很小的，绕组中的铜耗也很小 ， 所以近似认为变压器的铁耗就等于空载损耗， 则公式为 ：-32G T S= AP0 KW X10 /V N KV4. 电纳Bt：变压器的电纳代表变压器的励磁功率.变压器空载电流包含有功分量和无功分量，与励磁功率对应的是无功分量. 由于有功分量很小 ， 无功分量和空载电流在数值上几乎相等.B T S=I 0% XS N KVA5. 变压比Kt：在三相电力系统计算中，变压器的变压比通常是指两侧绕组空载线电压的比值.对于星形和三角形接法的变压器， 变压比与原副方绕组匝数比相等； 对于星三角形接法的变压器， 变压比为原副方绕组匝数比的 .3 倍.根据电力系统运行调节的要求， 变压器不一定工作在主抽头上 ， 因此 ,变压器运行中的实际变比 ， 应是工作时两侧绕组实际抽头的空载线电压之比三绕组变压器的参数计算三绕组变压器等值电路中的参数计算原则与双绕组变压器的相同 ， 下面分别确定各参数的计算公式 .1. 电阻R1,R2 ,R3 ： 为了确定三个绕组的等值阻抗， 要有三个方程 ， 为此 ， 需要有三种短路试验的数据.三绕组变压器的短路试验是依次让一个绕组开路， 按双绕组变压器来作 .通过查手册可得短路损耗分别为 ， 则有A S1 = 1/2( APS(1-2)+ APS(3-1)- APS(2-3)APS2=1/2( APS( 1-2)+ 护 S(2-3) - APS(3-1)APS3=1/2( APs(2-3)+ APs(3-1)-护 S(1-2)求 出 各 绕 组 的 短 路 损 耗 后 ,便 可 导 出 双 绕 组 变 压 器 计 算 电 阻 相 同 形 式 的 算 式 即 :R i Q=Psi KWV J KV X10/S KVA2. 电抗 X 1,X 2,X 3:和双绕组变压器一样,近似地认为电抗上的电压降就等于短路电压.在给出短路电压力 后,与电阻的计算公式相似,各绕组的短路电压为Vs1%=1/2(V s(1-2) %+Vs(3-1) %-Vs(2-3) %)Vs2%=1/2(V s(1-2) %+Vs(2-3) %-Vs(3-1) %)Vs3%=1/2(V s(2-3) %+Vs(3-1) %-Vs(1-2) %)各绕组的等值电抗为 ： X i Q=VSi%cXV N?KV Vb2/Sb=0.0227 10 -3X2302/100=0.012d11) Rt=547 1102B4 (型号 SFP7-150000/110 )(参数 AP=107KW, AP=547KW,Vs%=13%,I %=0.6,YNX103/150000 2=0.2942 QRt*=0.2942 100/115 2=0.0022Xt=13 X1102 X103/100/150000=10.487 QXt*=10.487 100/115 2=0.0793Bt=0.6 150000 X0 3/100/11C2=0.744 10 -4SBt*=0.744 10 -4 X1152/100=0.00982、三绕组变压器:B5B8 (型号 SFPSZ7-120000/220220_8 1.25%/115/10.5YN,y n0,d11=23.5% U 23 =13.3%53=7.7%P，= 425： P90)1Usi%(23.5 13.3 -7.7) =14.55%1Usn %(7.7 13.3 - 23.5) - -1.25% : 01Usm%(23.5 7.7 -13.3) =8.95% R tq= APskwVn2kv X103/S n2kva=425 X 2302 X 103/120000 2=1.5613Rt* =1.5613 100/230 2= 0.002951214.552202Xi220103 = 58.685100120000it远!C3=36.098. 1100 120000 58.685 100xi*-0.1109 = Xci230236.098 100 0.068 = X c 皿2302xi* Xn* =0.1109242/121K1230/115 300100X10231 /1103、 联络变压器B9: RT*20.00567 X T*20.2079 K0.95723022302230/1154、 线路L1L7: （L1,L5,L6为双同路 2,L3,L4,L7为单回路）L1 ： 型号为2X LGJ F00/75km . 参数为 ： r=0.08Q/km,x=0.397 Q/km,b=2.92 X 10-6s/kmR1=1/2L X r=75 X 0.08/2=3 QX1=1/2Lx=0.397 X 75/2=14.888 Q_6_4B1=2XLb=2 X2.92 X06X75=4.38 X04S2 2 2 2R1*=RXSB/VB =3X100/230 =0.00567X1*=X1 XSB/VB =14.888 X00/230 =0.02812-42B1*=B1XVB2/SB=4.38 为 042302/100=0.2317 B 1*/2=0.1158L2：型号为 LGJ F00/50km.参数为：r=0.08Q/km,x=0.397 Q/km,b=2.92 X 10-6s/km-6-4R2=LXr=50 X0.08=4 Q X2=Lx=0.397 50=19.85 QB2=Lb=2.92 10 X50=1.46 10 SR2*=RXSB/VB2=4X100/1152=0.0302X2*=X2XSB/VB2=19.85 X00/115 2=0.1501B2*=B2 XVB2/SB=1.46 为 0-4 %15 2/100=0.0193B2*/2=0.00965L3：型号为 LGJ 85/70km.参数为：r=0.17 Q/km,x=0.395 Q/km,b=2.79 X 10-6s/km-6-4R3=LXr=70 X0.17=11.9 Q X3=Lx=0.395 X70=27.65 QB3=Lb=2.79 X106 *0=1.953 10 4SR3*=RXSB/VB2=11.9 X00/11#=0.09X3*=X3XSB/VB2=27.65 X00/115 2=0.2091B3*=B3XVB2/SB=1.953 10-4X115/100=0.0258B3*/2=0.0129L4：型号为 LGJ 50/74km.参数为：r=0.21 Q/km,x=0.403 Q/km,b=2.74 X 10-6s/km-6-4R4=LXr=74 X0.21=15.54 Q X4=Lx=0.403 74=29.822 QB4=Lb=2.74 X10 X74=2.03 10 SR4*=RXSb/Vb2=15.54 100/115 A=0.1175X4*=X4 0/Vb2=29.822 100/115 2=0.2255L5:型号为 2XLGJ -20/70km.B4*=B4XVB2/SB=2.03 为 0-4%152/100=0.0268B4*/2=0.0134参数为 ： r=0.27 Q/km,x=0.409 Q/km,b=2.69 X 10-6s/kmR5=1/2L X r=70 X 0.27/2=9.45 Q-6-4B5=2XLb=2X2.69 10 6 八0=3.766 10 4sX5=1/2Lx=0.409 X 70/2=14.315 QR5*=R5 XSb/Vb2=9.45 八00/115 八=0.0714X5*=X5 XSb/Vb2=14.315 100/115 2=0.1082B5*=B5 XVB2/SB=3.766 10 -4 X115/100=0.0498 B5*/2=0.0249L6：型号为 2 XLGJ Y00/73km .参数为：r=0.08 Q/km,x=0.397 Q/km,b=2.92 X 10-6s/kmR6=1/2L X r=73 X 0.08/2=2.92 QX6=1/2Lx=0.397 X 73/2=14.491 Q-6-4B6=2 XLb=2 X2.92 10 6 X73=4.26 10 4S2 2 2 2R6*=R6 XSB/VB =2.92 X00/230 =0.0055X6*=X6 XSB/VB =14.491 100/230 =0.02742-42B6*=B6 XVB2/SB=4.26 X0 4 X2302/100=0.2254B6*/2=0.1127L7：型号为 LGJ F00/135km.参数为：r=0.08 Q /km,x=0.397Q/km,b=2.92 X10-6s/kmR7=LXr=135 X 0.08=10.8 CX 7=Lx=0.397 X35=53.595 QB7=Lb=2.92 1O -6X135=7.844 采 0-4S2 2 2 2R7*=R7XSb/Vb =10.8 X00/230 =0.0204X7*=X7Xsb/Vb =53.595 X00/230 =0.10132-42B7*=B7XVb2/Sb=7.844 X0 4X2302/100=0.4149B7*/2=0.207455、发电机 F1-F4:(型号：SF-100/40-854Pe=100MWcos =0.9 Ue=13.8KVXd* =0.2034Xo* -0.0 9 2 76 PQ节点PQ标么值:4S4=180+j100S4*=S4/Sb = 1.8+j1.0S5=72+j408 S8=120+j639 _S9=36+j2010S10=40+j25S5*=S5/SB=0.72+j0.4S8*=S8/SB=1.2+j0.63Sg*=S9/SB=0.36+j0.2S10*=S10/SB=0.4+j0.27、零序电抗标么值Xo (l1)=4.5Xl1 =4.5 0.0281 =0.12645Xo (l2) - 3Xl2 =3 0.1501 = 0.4503X (L3)=3Xl3 =3 0.2091 = 0.6273Xo (l4)=3Xl4 -3 0.2255 =0.6765Xo (l5)=3Xl5 =3 0.1082 =0.3246Xo (l6)=4.5Xl6 =4.5 0.0274=0.1233Xo (l7)=4.5Xl7 =4.5 0.1013 = 0.45585电力系统潮流计算P Q 分解法潮流计算主要步骤1. 导纳矩阵的形成2. 因子表的形成（三角分解法）3. 给定电压初值4. 计算不平衡功率4Pi（k）除以V严）5. 计算不平衡功率4 Qi的除以V严）6. 判断收敛性7. 回代，修正V , Vi（k）8. 计算全线路功率电力系统的潮流计算机算法潮流计算的任务：对给定的运行条件确定运行状态， 如果各母线上的电压（ 幅值及相角 ） 网络 中的功率分布以及功率损耗等，几个节点电力系统的潮流方程的一般形式根据电力系统的实际运行条件，按给定变量的不同，一般节点可分三种：1. PQ节点：有功功率P和无功功率Q给定的，节点电压（V, S）是待求量。2. PV节点：有功功率P和电压幅值V给定的，节点的无功功率Q电压的相位角S是待求量3. 平衡节点：网络中至少有一个节点的有功功率 P 不能给定，此节点承担系统的有功功率平衡。网络方程是线性方程， 但由于在定解条件中不能给定节点电流， 只能给出节点功率， 从而使潮 流方程变为非线性方程，由于平衡节点的电压已经给定，假设系统中有n 个节点，其中有m 个 PQ 节点， n-（m+1）个PV节点和一个平衡节点，平衡节点不参加求解。形成导纳矩阵的方法及框图1 .形成导纳矩阵的程序框图为了形成导纳矩阵，必须个计算机输入必要的原始数据，一条支路的原始数据应包括两端的节点号和支路阻抗，对变压器支路的原始数据应包括他的变比，对于电力线路还应包括它的容纳，即：I , J , R , X , K( -B/2 )，对以上数据说明如下：(1) 当支路为接地支路时，规定节点I处置零，J处填写接地支路所在的节点号，支路参数用阻抗填写,对于K处置零，但必须填写。(2) 对于纯阻抗支路，K处置零。(3) 对于具有容纳的线路，电容电纳不作单独支路处理，而且把总容纳的一半负值填写在K处，其符号用以区别支路的性质。(4) 对于变压器支路，采用下图所示等值电路，即非标准变比在J侧，变压器阻抗在I侧，K处填写实际变化YB/Kl:k2 .框图P=1I开始I .导纳矩阵清零n k=1 , L取支路K数福I二J，RX，K (-B/2 )是接地 路吗？是变压器支路吗？P=k、 Q=k线路容纳追加到两节点自导纳中P=1 , Q=1支路导纳追加到节点的自导纳中卜J之间互导纳中将支路导纳追加至J节点- 的自导纳中输出结果注：Bi 不记充电电容和非标准变比时，与纳矩阵的虚部B2与纳矩阵的虚部、形成因子表的方法及框图1用行消去过程形成因子表的程序框图2.因子表的形成（三角分解法）消去法求非线性方程组的一种常用算法是对方程式的系数矩阵A进行三角分解，在本次电力系统潮流计算中采用的三角分解是将非奇方阵A分解为单位下三角矩阵L和上三角矩阵R的乘积。A=LR ,非奇方阵 A被表示为矩阵 L和R的乘积：这两个三角矩阵称为A的因子矩阵，两个因子矩阵的元素计算公式：j -1Lj(aj - xlik rkj )k =1r1 ijrja iji =1昱1 ik r kj k =1 1, 2 , -nj i , i+ 1 ,n将A=LR代入线性方程组，便得 LRX=B ,这个方程又可以分解为以下两个方程 LF=BRX=F展开为1L 21L 31L 41L 32l42L43 1LniLn2Ln312P 22M，a ar 13r 2nnnfn先由方程组自上而下地依次算出f2、fn,其计算通式为bii - 1送L i.j =这一步演算相当于消元过程中对原方程式右端常数向量所作的变只顺用到下三角因子矩阵。方程组的求解属于回代过程，只顺用到上三角因子矩阵以及经过消元变换的右端常数向量，方程组可以自下而上地逐步 算出待求量，其计算通式为：(fi(TXin一(i 1)ij()Xj )j =i(i - 1 )三、节点电压的表示：i.极坐标表示法：Vi 二丫 二 Vii)节点功率方程表示(sin为：Vj (Gj cos 川Bij sin j)QiVj(Gjsin j-Bijcosij)H二i-j (两节点电压的相位角)方程式把节点功率表示为节点电压的幅值和相角的函数在有n个节点的系统中，假定第 1m号节点为PQ节点，第m+1n-1号节点为平衡节点。 V和3 n是给定 的。PV 节点的电压幅值 Vm+iVn-i也是给定的。因此，只剩下 n-1个节点的电压相角3仁？ 3 n-i和m个节点的电 压幅值V1 Vm 是未知的量。Bj sin jh0 (i 2n-1) 对于每对于每一个PQ节点或每一个 PV节点都可以列写一个有功功率不平衡方程：p = Ps - p = Ps - ViVj 9aos 可jT一个PQ节点可以列写一个无功功率不平衡方程式：g =Q is-Q =Q is-VQ Vj (Gj sin% jyBj C0Sj)= 0 (i = 1,2 )m)所以可以写出修正方程式其中gp 111Ap211 :111r ii i2111Ad11:111d1 111I1 Ax,V2YD?匕1 1113m.V3JP和二还eVN、k略去不计，既认为它四、P-Q分解法1 ?在交流高压电网中，输电线路的电抗要比电阻大得多，系统中母线有功功率的变化主要受电压相位的影响，无功功率的变化规则主要受母线电压幅值变化的影响，在修正方程式的系数矩阵中，偏导数 的数值相当于偏导数二二 E和：Q是相当小的，作为简化第一步，可将方程中V们的元素都等于零，这样，便可分解为n-1阶和m阶的两个方程：P - -H Q = -LVDA: V所谓P-Q分解法，节点有功功率不平衡量又用于修正电压的相位。节点的无功功率不平衡量值用于修正电压幅值，方程、分别轮流迭代。2. H、L元素都是节点电压幅值和相角差函数，其数值在迭代过程中是不断变化的，从而，最关键的一步是把系数矩阵H、L简化为常数矩阵。方法：在一般情况下，在线路两端电压的相角差是不大的（不超过0。20*）,因此，以为cos ij =1 , Gij sin : ij - Bij,止匕外，与系统各节点无功功率相适应的导纳BLdi必须小于该节点自导纳虚部， 即Qi2BLDi = V2 - Bij 或 Qi - Vi BijVi?矩阵H、L可简化为H AVDi BVdiL - VD2 B VD2将分别带入 得：AP = VdQVd,心Vd： 1AP = -BVd,仝=VD2B“迪 V I VDAQ = BAV即：PiVi*V2 mAPnAVn1 B12 BI3 Bl,nJ VU1B22 B2,nV V2 也合 2：I | ：I IU,1 B.4,2 昭 2必异% _1 一校验收敛应就完成分解法的计算了Vi| V2 aAQm.B11B12 B13Blm1,DD D|B21B22B2m1 aaaaa(D D D-A11Dm2韦 m-3 11 1FV2 11: 1II3m-Vm利用计算节点功率的不平衡量,用修正方程解出修正量 3及4 V,弁换下述条件max 也 Pj( kJm aXAQ# 讥名 q3. P-Q分解法计算潮流的程序框图输入原始数据形成矩阵ET及形成B：闪子表设P（电点电压初值，各节点电冰相由初值笆迭代次数K=0| i| 算 AR/Vi (N-1)个+_想/V. ( m )个迭代收敛了吗？. 计算支路功率一利用B1因子表对后端项变换回代求解计算支路功率损耗修正电压角输出.果利用B对右端项变换_ +J结.回代求解 I修正电压模k=k+1PQ 分解法源程序#include math.h#include stdio.h#include stdlib.hmain ()int a=1,k,i,j,d=1,m,n,ls,nl,kk;float z1005,pq1003,p,q,r,x,kb,b,g;float yg100100,yb100100,b1100100,b2100100,v100,w100;float pp100,qq100,e100,f100,gg1002,tt1002,pw100,qw100,ww100;float t,t1,tp,tq,p1,p2,q1,q2,t2;printf(n);printf( 电力系统潮流计算n);printf(n);printf(指导老师：设计：2002年12月n);printf(n);printf(n);printf(注：各参数用标么值表示n);printf(n);begin:printf(n); printf( 请选择 : 1 - 建立电网 ;2-查看数据 ;3-修改数据 ;4-运行计算;5-退出程序 :);scanf(%d,&a);if (a=1)goto creat;else if (a=2)goto view;else if (a=3) goto modify;else if (a=4)goto run; elseexit(0);creat:printf( 请输入系统的节点数N:);scanf(%d,&n);printf( 请输入系统的支路数LS:);scanf(%d,&ls);printf( 请输入系统的 PQ 节点数 M:);scanf(%d,&m);for(k=1;k=ls;k+)printf( 请输入支路的 %d 的数据 I( 接地支路为 0) :,k);scanf(%d,&i);printf( 请输入支路的 %d 的数据 J:,k);scanf(%d,&j);printf( 请输入支路的 %d 的数据 R:,k);scanf(%f,&r); printf( 请输入支路的 %d 的数据X:,k);scanf(%f,&x);printf( 请输入支路的 %d 的数据 KB:,k);scanf(%f,&kb);zk1=i;zk2=j;zk3=r;zk4=x;zk5=kb;for(k=1;k=m;k+)pqk1=k;printf(请输入 PQ节点 d的注入有功P:,k);scanf(%f,&p); printf( 请输入 PQ节点 d的注入无功Q:,k);scanf(%f,&q);pqk2=p;pqk3=q;for(k=m+1;k=n-1;k+)pqk1=k;printf( 请输入 PV 节点 %d 的注入有功 P:,k);scanf(%f,&p);printf( 请输入 PV 节点 %d 的电压模 :,k);scanf(%f,&q); pqk2=p;pqk3=q;pqk1=n;printf( 请输入平衡节点%d 的电压模:,n);scanf(%f,&p);printf( 请输入平衡节点%d 的电压角:,n);scanf(%f,&q); pqn2=p;pqn3=q;goto begin; view:printf(n);printf(N=%dtLS=%dtM=%d,n,ls,m);printf(n);printf( 支路数据 );printf(n);printf(t 节点 It 节点 J 电阻 t 电抗 t 变比或电纳 );printf(n);for(i=1;i=ls;i+)printf( 支路 %d,i);for(j=1;j=5;j+)if (j=2)printf(t %1.0f,zij);elseprintf(%14.6f,zij);printf(n);printf(tPQ 节点 t 注入有功 t 注入无功 );printf(n); for(i=1;i=m;i+)for(j=1;j=3;j+)if (j=1)printf(t %1.0f,pqij);else printf(t%14.6f,pqij); printf(n);printf(tPV 节点 t 注入有功 t 电压模 ” p);rintf(n);for(i=m+1;i=n-1;i+) for(j=1;j=3;j+) if (j=1)printf(t %1.0f,pqij);else printf(t%14.6f,pqij); printf(n);printf(n);printf(t 平衡节点 t 电压模 t 电压角 );printf(n);for(j=1;j=3;j+) if (j=1)printf(t %1.0f,pqnj);else printf(t%14.6f,pqnj);printf(n); goto begin;modify: printf(1- 修改支路数据;2- 修改节点注入功率;3- 返回主菜单:);scanf(%d,&d); if (d=1)goto m1;else if (d=2) goto m2; else if (d=3)goto begin;m1: printf( 请输入要修改的支路号:);scanf(%d,&k);printf( 请输入支路的printf( 请输入支路的printf( 请输入支路的%d 的新数据I:,k);scanf(%d,&i);%d 的新数据J:,k);scanf(%d,&j);%d 的新数据R:,k);scanf(%f,&r);printf( 请输入支路的printf( 请输入支路的%d 的新数据X:,k);scanf(%f,&x);%d 的新数据KB:,k);scanf(%f,&kb);zk1=i;zk2=j;zk3=r;zk4=x;zk5=kb;goto modify;m2: k=0;printf( 请输入需修改的节点号:);scanf(%d,&k);pqk1=k;printf( 请重新输入节点的注入有功或电压参数printf( 请 :);scanf(%f,&pqk2);重新输入节点的注入无功或电压参数goto modify;:);scanf(%f,&pqk3);run : for(i=1;i=n;i+) for(j=1;j=n;j+)ybij=0;/*形成 B1 ，不计线路电容和变压器非标准变比 */ for(k=1;k=ls;k+) i=zk1;j=zk2;r=zk3;x=zk4;b=-x/(r*r+x*x);if(i=0) goto p0;ybii=ybii+b;ybij=ybij-b;ybji=ybij; p0:ybjj=ybjj+b; for(i=1;in;i+) for(j=1;j=n;j+) b1ij=ybij; /* 形成导纳矩阵*/for (i=1;i=n;i+)for (j=1;j=n;j+)ygij=0; ybij=0;for (k=1;k0)goto p2;ybii=ybii-kb;ybjj=ybjj-kb;p=1;q=1;goto p3;p1: p=1;goto p4; p2:p=kb*kb;q=kb;p3:ygii=ygii+g;ybii=ybii+b; ygij=ygij-g/q;ybij=ybij-b/q; ygji=ygij;ybji=ybij;p4:ygjj=ygjj+g/p;ybjj=ybjj+b/p; /* 形成 b2*/for(i=1;i=m;i+)for(j=1;j=m;j+)b2ij=ybij;printf ( 矩阵 B:n);for(i=1;i=n;i+)for(j=1;j=n;j+)printf(%14.7f,ybij);printf(n);printf ( 矩阵 G:n);for(i=1;i=n;i+)for(j=1;j=n;j+) printf(%14.7f,ygij); printf(n);/*形成bl 因子表 */for(i=2;i=n-1;i+)(for(k=1;k=i-1;k+) b1ik=-b1ik/b1kk; for (j=k+1;j=n-1;j+)b1ij=b1ij+b1ik*b1kj;)printf(B1 因子表:n);for(i=1;i=n-1;i+)for(j=1;j=n-1;j+) printf(%14.7f,b1ij);)printf(n);)/*形成b2因子表:*/for(i=2;i=m;i+)for(k=1;k=i-1;k+)b2ik=-b2ik/b2kk;for (j=k+1;j=m;j+)b2ij=b2ij+b2ik*b2kj; )printf(B2 因子表:n);for(i=1;i=m;i+)for(j=1;j=m;j+)printf(%14.7f,b2ij);)printf(n);)printf(n);/*给定电压初值*/for(i=1;i=m;i+)vi=1;wi=0;)for(i=m+1;i=n-1;i+)vi=pqi3;wi=0;)vn=pqn2;wn=pqn3;kk=0;/*计算 n-1 个节点不平衡有功并除其电压*/a1:for(i=1;i=n-1;i+) t=0;for(j=1;j=n;j+) if (ybij!=0)t=t+vj*(ygij*cos(wi-wj)+ybij*sin(wi-wj);ppi=pqi2/vi-t;/* 计算m 个节点不平衡无功并除其电压 */for(i=1;i=m;i+) t=0;for(j=1;j=n;j+) if (ybij!=0) t=t+vj*(ygij*sin(wi-wj)-ybij*cos(wi-wj); qqi=pqi3/vi-t;for(i=1;i0.00001)goto a4 ;for(i=1;i0.00001)goto a2;goto a3;/* 利用 B1 因子表对右端项变换*/a4:for(i=1;i=n-2;i+) t=0;for(j=1;j=i;j+) t=t+b1i+1j*ppj; ppi+1=ppi+1+t;nl=n-1;/* 回代求解 */ppnl=ppnl/b1nlnl;for(k=1;k=nl-1;k+)t=0;for(j=nl-k+1;j=nl;j+) t=t+b1nl-kj*ppj;ppnl-k=(ppnl-k-t)/b1nl-knl-k;/* 修正电压角 */for(i=1;i=nl;i+)wi=wi-ppi;/* 利用 b2 因子表对右端项变换*/a2:for(i=1;i=m-1;i+)t=0;for(j=1;j=i;j+) t=t+b2i+1j*qqj;qqi+1=qqi+1+t;qqm=qqm/b2mm;/* 回代求解 */for(k=1;k=m-1;k+)t=0;for(j=n-k+1;j=m;j+) t=t+qqj*b2m-kj;qqm-k=(qqm-k-t)/b2m-km-k;/* 修正电压模*/for(i=1;i=m;i+)vi=vi-qqi;kk=kk+;goto a1;a3:for(i=1;i=n;i+)ei=vi*cos(wi);fi=vi*sin(wi);wwi=wi*180/3.1415926;/* 计算支路功率*/for(k=1;k0) ggk1=p1+(kb-1)*vi*vi*(-ygij); ttk1=q1-(kb-1)*vi*vi*(-ybij);ggk2=p2+(1-kb)*vj*vj*(-ygij)/kb; ttk2=q2-(1-kb)*vj*vj*(-ygij)/kb;elseggk1=p1;ggk2=p2;/* 有功 */ ttk1=q1+kb*vi*vi;/* 无功 */ ttk2=q2+kb*vj*vj;/* 计算支路功率损耗*/ for(k=1;k=ls;k+) pwk=ggk1+ggk2; qwk=ttk1+ttk2;tp=0;tq=0;for(k=1;k=ls;k+) tp=tp+pwk;tq=tq+qwk;printf( 迭代 %d 次 电压实部 e 电压虚部 f 电压模 V 电压角(度 )n,kk);for(i=1;i=n;i+) printf( 节点 %6d%14.7f%14.7f%14.7f%14.7fn,i,ei,fi,vi,wwi);printf(n);printf( 支路功率 : 有功 无功 n);for(k=1;k=ls;k+) i=zk1;j=zk2; printf ( 支路%d-%d%16.7f%18.7fn,i,j,ggk1,ttk1);printf ( 支路%d-%d%16.7f%18.7fn,j,i,ggk2,ttk2);printf( 支路损耗 : 有功 无功 n); for(k=1;k=ls;k+)i=zk1;j=zk2;printf ( 支路%d-%d%16.7f%18.7fn,i,j,pwk,qwk);printf( 总网耗 : 有功 无功 n); printf(%24.7f%18.7f,tp,tq);goto begin;N=11 LS=11 M=9支路数据潮流结果电阻电抗0.028100-0.1158000.0554501.0000000.150100-0.0096500.1109001.0000000.2079000.9570000.209100-0.0129000.225500-0.0134000.108200-0.0249000.1109001.0000000.027400-0.1127000.101300-0.207450注入无功变比或电纳节点I节点 J支路 1130.005670支路 23100.001476支路 31050.030200支路 4560.002951支路 5120.005670支路 6240.090000支路 7490.117500支路 8980.071400支路 9870.002950支路 10760.005500支路 11 6110.020400PQ节点注入有功12.3000001.60000020.8000000.5000003-2.60000