电力系统下课程设计_短路电流计算

1、下载可编辑电力系统分析 课程设计报告题目： 3G9bus 短路电流计算系别电气工程学院专业班级10 级电气四班学生姓名学号指导教师提交日期2012年12月10日.专业 .整理 .下载可编辑目录一、设计目的3二、短路电流计算的基本原理和方法32.1 电力系统节点方程的建立32.2 利用节点阻抗矩阵计算短路电流5三、3G9bus 短路电流在计算机的编程73.1、三机九节点系统73.3 输出并计算结果20四总结22.专业 .整理 .下载可编辑一、设计目的1.掌握电力系统短路计算的基本原理；2.掌握并能熟练运用一门计算机语言（MATLAB 语言或 FORTRAN 或 C 语言或 C+ 语言）；3.采用

2、计算机语言对短路计算进行计算机编程计算。二、短路电流计算的基本原理和方法2.1 电力系统节点方程的建立利用节点方程作故障计算，需要形成系统的节点导纳（或阻抗）矩阵。 一般短路电流计算以前要作电力系统的潮流计算，假定潮流计算的节点导纳矩阵已经形成，在此基础上通过追加支路的方式形成电力短路电流计算的节点导纳矩阵YN。1）对发电机节点在每一发电机节点增加接地有源支路Ei 与 Zi RijX i 串联求短路稳态解 : Ei EQiZiRij Xq i求短路起始次暂态电流解 : EiEiZiRi j Xi一般情况下发电机定子绕组电阻忽略掉，并将 Ei 与 ZiRijX i 的有源支路转化成电流源Yi G

3、iBi1I i Ei ZiRijX i与导纳并联的形式.专业 .整理 .下载可编辑2）负荷节点的处理负荷节点在短路计一算中一般作为节点的接地支路，并用恒定阻抗表示 ，其数值由短路前瞬间的负荷功率和节点实际电压算出，即首先根据给定的电力系统运行方式制订系统的等值电路 ，并进行各元件标么值参数的计算，然后利用变压器和线路的参数形成不含发电机和负荷的节点导纳矩阵YN 。2?VkZ LDkRLDk jX LDkYL D kG L D kSL D k?j B L D k2SLDkVk3）如果短路故障发生在输电线路内，在短路点要增加新节点 ，将一条输电线分成两段 ，并形成短路电流计算的节点导纳矩阵Y，最后

4、形成包括所有发电机支路和负荷支路的节点方程如下：YVI一般，节点导纳矩阵Y 阵与YN 阶次相同 。 其差别只在于YN 阵不含发电机和负荷的阻抗 ；节点注入电流向量I 中只有发电机端节点的电流不为零。有非零电流源注人的节点称为有源节点。系统中的同步调相机可按发电机处理。任进行起始次暂态电流计算时，大型同步电动机 、感应电动机以及以电动机为主要成分的综合负荷，特别是在短路点近处的这些负荷 ，必要时也可以用有源支路表示，并仿照发电机进行处理。必须指出 ，在计算机已普遍应用的情况下，如果有必要的话 ，只要能提供短路计算所需的准确的原始数据，对短路进行更精确的计算并不存在什么障碍。.专业 .整理 .下载

5、可编辑2.2 利用节点阻抗矩阵计算短路电流假定系统中的节点f 经过渡阻抗 zf 发生短路 。过渡阻抗渡阻抗zf，不参与形成网络的节点导纳 （或阻抗 ）矩阵 。图 6-3 中方框内的有源网络代表系统正常状态的单相等值网络。 现在我们保持故障处的边界条件不变，把网络的原有部分同故障支路分开（见图6-3 ）。容易看出 ，对于正常状态的网络而言，发生短路相当于在故障节点f 增加了一个注人电流一If （短路电流以流出故障点为正，节点电流则以注入为正）。因此，网络中任ViZij I jZifI f (i 1,2,., f ,., n)一节点 i的电压可表示为 ：j G（6-3 ）由式（6-3 ）可见，任一

6、节点 i的电压郁由两项叠加而成 。第一项表示当 If=0 时由网络内所有电源在节点 i 产生的电压 ，也就是短路前瞬间正常运行状态一F 的节点电压 ，这是节点电压的正常分量，记为 Vi(0)是电网的潮流解 。 第二项是当网络中所有电流源都断开，仅仅由短路电流 If 在节点 i产生的电压 ，这就是节点电压的故障分童 。 上述两个分量的叠加，就等于发生短路后节点，的实际电压 ，即Vi Vi(0)Zif I f( i 1,2,., f ,., n)V fV(f0 )Zi f I f （6-4 ）式中，Vi(0)Z ff是故障节点 f的自阻抗 ，也称是短路前故障点的正常电压 ；当 i=f 时，输入阻抗

7、 。 V fzf I f（边界条件方程 ）带入 Vf V f(0)Zif I f 可得：I fV f(0)zfZ ffVi Vi(0)ZifVf(0)(i 1,2,., f ,., n)Z ffZ fkVpVqI pqZ pq.专业 .整理 .下载可编辑对于非变压器支路，令 k=1 。 从计算公式 （ 6-7 ）和（ 6-8 ）可以看到 ，式中所用到的阻抗矩阵元素都带有列标f。 这就是说 ，如果网络在正常状态下的节点电压为已知，为了进行短路计算 ，只须利用节点阻抗矩阵中与故障点f 对应的一列元素 。因此 ，尽管是采用了阻抗型的节点方程，但是并不需要作出全部阻抗矩阵。在短路的实际计一算中，一般只

8、需形成网络的节点导纳矩阵并根据具体要求 ，用第四章所讲的方法求出阻抗矩阵的某一列或某几列元素即可 。 在应用节点阻抗矩阵进行短路计算时，我们都将采用这种算法 。1） 解潮流计算 ， Vi(0)(i1,2,., f ,., n)2） 修正解潮流的 YN 形成 Y3） 指定短路点 f0YYYZ11121n1k0Y21Y22Y2nZ2 k1Yn1Yn 2YnnZnk04） 计算节点阻抗矩阵第 k 列0I fVf(0)Z ffzf5） 计算短路电流 ：6） 计算节点电压 ： ViVi(0)ZifI f (i1,2,., f ,., n)I pqkVpVq7） 计算支路电流 ：Z pq8） 输出计算结果

9、第 f 元素 ,( k1, 2,n).专业 .整理 .下载可编辑三、 3G9bus 短路电流在计算机的编程3.1 、三机九节点系统79283654f(3)1图 1 三机九节点系统表 1 九节点系统支路参数支路R（p.u.）X（ p.u.）B/2 （p.u.）1400.05761.02700.06251.03900.05861.0450.010.0850.088460.0170.0920.079570.0320.1610.153690.0390.170.179780.00850.0720.0745890.01190.10080.1045.专业 .整理 .下载可编辑表 2九节点系统发电机参数PGV

10、GX d发电机编号节点类型（p.u.）（ p.u.）（ p.u.）E （p.u.）1V1.040.31.1372PV1.631.0250.31.2113PV0.851.0250.31.043表 3九节点系统负荷参数节点编号节点类型PiQi（ p.u.）（ p.u.）4PQ005PQ1.250.56PQ0.90.37PQ008PQ10.359PQ00主程序3.2 程序设计主函数Sbase_MVA=100.fid=fopen('Nodedata.txt');.专业 .整理 .下载可编辑N=textscan(fid, '%s %u %d %f %f %f %f %f %f&#

11、39;)fclose(fid);busnumber=size(N1,1)for i=1:busnumberBus(i).name=N1(i);Bus(i).type=N2(i);Bus(i).no=i;Bus(i).Base_KV=N3(i);Bus(i).PG=N4(i);Bus(i).QG=N5(i);Bus(i).PL=N6(i);Bus(i).QL=N7(i);Bus(i).pb=N8(i);Bus(i).V=1.0;Bus(i).angle=0;endfid=fopen('Aclinedata.txt');A=textscan(fid, '%s %s %f %

12、f %f %f')fclose(fid);aclinenumber=size(A1,1)for i=1:aclinenumberAcline(i).fbname=A1(i);.专业 .整理 .下载可编辑Acline(i).tbname=A2(i);Acline(i).Base_KV=A3(i);Acline(i).R=A4(i);Acline(i).X=A5(i);Acline(i).hB=A6(i);for k=1:busnumberif strcmp(Acline(i).fbname, Bus(k).name)Acline(i).fbno=Bus(k).no;endif strcm

13、p(Acline(i).tbname, Bus(k).name)Acline(i).tbno=Bus(k).no;endendendfid=fopen('Transdata.txt');T=textscan(fid, '%s %f %f %s %f %f %f %f')fclose(fid);tansnumber=size(T1,1)for i=1:tansnumberTrans(i).fbname=T1(i);Trans(i).fbBase_KV=T2(i);Trans(i).fbrated_KV=T3(i);.专业 .整理 .下载可编辑Trans(i).tb

14、name=T4(i);Trans(i).tbBase_KV=T5(i);Trans(i).tbrated_KV=T6(i);Trans(i).R=T7(i);Trans(i).X=T8(i);for k=1:busnumberif strcmp(Trans(i).fbname, Bus(k).name)Trans(i).fbno=Bus(k).no;endif strcmp(Trans(i).tbname, Bus(k).name)Trans(i).tbno=Bus(k).no;endendTrans(i).k=Trans(i).tbrated_KV*Trans(i).fbBase_KV/Tr

15、ans(i).fbrated_KV/Trans(i).tbBase_KV;tempx=Trans(i).fbrated_KV2/Trans(i).fbBase_KV2;Trans(i).X=tempx*Trans(i).X;Trans(i).R=tempx*Trans(i).R;end%N=0%Trans(1)%Trans(2).专业 .整理 .下载可编辑% for Y=G+jB matrix G,B,B2=FormYmatrix(Bus,busnumber,Acline,aclinenumber,Trans,tansnumber); %B:=B'B2:=B"dlmwrite

16、('Gmatrix.txt', G, 'delimiter', 't','precision', 6); dlmwrite('Bmatrix.txt', B, 'delimiter', 't','precision', 6);G B B2 pause JP,JQ=FormJPQmatrix(Bus,B,B2,busnumber);JP iJP=-inv(JP)JQ iJQ=-inv(JQ) pause %maxiteration=0 for i=1:busnumber

17、NodeV(i)=Bus(i).V; Nodea(i)=Bus(i).angle; VX(i)=Bus(i).V*cos(Bus(i).angle); VY(i)=Bus(i).V*sin(Bus(i).angle); dQGQL(i)=Bus(i).QG-Bus(i).QL;.专业 .整理 .下载可编辑dPGPL(i)=Bus(i).PG-Bus(i).PL;endNodeV=NodeV'Nodea=Nodea'%VX=VX'%VY=VY'dQGQL=dQGQL'dPGPL=dPGPL'pause%for nointer=1:10maxdP=1

18、.;maxdQ=1.;epsilon=0.000001;noiteration=0;while (maxdP>epsilon)&(maxdP>epsilon)deltaP,deltaQ,maxdP,maxdQ=FormdPQvector(Bus,NodeV,Nodea,dQGQL,dPGPL,B,G,busnumber);deltaP;deltaQ;maxdP;maxdQ;da=iJP*deltaP;.专业 .整理 .下载可编辑dV=iJQ*deltaQ;Nodea=Nodea+da;NodeV=NodeV+dV;noiteration=noiteration+1;if n

19、oiteration>20breakendendfor i=1:busnumberBus(i).V=NodeV(i);NodeV(i)=NodeV(i)*Bus(i).Base_KV;Bus(i).angle=Nodea(i);Nodea(i)=Nodea(i)*180/pi;endnoiterationNodea=Nodea'NodeV=NodeV'Clear子函数%生成 G、B 矩阵functionG,B,X=FormYmatrix(Bus,busnumber,Acline,aclinenumber,Trans,tansnumber).专业 .整理 .下载可编辑Y=z

20、eros(busnumber);X=zeros(busnumber);for i=1:busnumberY(i,i)=Y(i,i)+Bus(i).pb*j;endfor i=1:aclinenumberf=Acline(i).fbno;t=Acline(i).tbno;Y(f,f)=Y(f,f)+Acline(i).hB*j+1/(Acline(i).R+Acline(i).X*j);Y(t,t)=Y(t,t)+Acline(i).hB*j+1/(Acline(i).R+Acline(i).X*j);Y(f,t)=Y(f,t)-1/(Acline(i).R+Acline(i).X*j);Y(t

21、,f)=Y(t,f)-1/(Acline(i).R+Acline(i).X*j);X(f,f)=X(f,f)-1/Acline(i).X;X(t,t)=X(t,t)-1/Acline(i).X;X(f,t)=1/Acline(i).X;X(t,f)=1/Acline(i).X;endfor i=1:tansnumberf=Trans(i).fbno;t=Trans(i).tbno;Y(f,f)=Y(f,f)+1/(Trans(i).R+Trans(i).X*j);Y(t,t)=Y(t,t)+1/(Trans(i).R+Trans(i).X*j)/Trans(i).k2;.专业 .整理 .下载可

22、编辑Y(f,t)=Y(f,t)-1/(Trans(i).R+Trans(i).X*j)/Trans(i).k;Y(t,f)=Y(t,f)-1/(Trans(i).R+Trans(i).X*j)/Trans(i).k;X(f,f)=X(f,f)-1/Trans(i).X;X(t,t)=X(t,t)-1/Trans(i).X;X(f,t)=1/Trans(i).X;X(t,f)=1/Trans(i).X;endG=real(Y);B=imag(Y);end%生成 JP、JQ 矩阵function JP ,JQ=FormJPQmatrix(Bus,B,B2,busnumber)JP=B;JQ=B2;

23、for i=1:busnumberif Bus(i).type=1for k=1:busnumberJQ(i,k)=0.;JQ(k,i)=0.;JP(i,k)=0.;JP(k,i)=0.;.专业 .整理 .下载可编辑endJQ(i,i)=1.;JP(i,i)=1.;endif Bus(i).type=3for k=1:busnumberJQ(i,k)=0.;JQ(k,i)=0.;endJQ(i,i)=1.;endendend%计算偏节点 PQ 差量functiondeltaP,deltaQ,maxdP,maxdQ=FormdPQvector(Bus,NodeV,Nodea,dQGQL,dPGP

24、L,B,G,busnumber)deltaQ=dQGQL;deltaP=dPGPL;maxdP=0.;maxdQ=0.;.专业 .整理 .下载可编辑for i=1:busnumberif Bus(i).type=1deltaQ(i)=0.;deltaP(i)=0.;endif Bus(i).type=3deltaQ(i)=0.;%y1=0;%y2=0;y3=0;for k=1:busnumberif (B(i,k)=0|G(i,k)=0)%y1=y1+(G(i,k)*VX(k)-B(i,k)*VY(k);%y2=y2+(G(i,k)*VY(k)+B(i,k)*VX(k);y3=y3+NodeV

25、(k)*(G(i,k)*cos(Nodea(i)-Nodea(k)+B(i,k)*sin(Nodea(i)-Nodea(k);endenddeltaP(i)=deltaP(i)-y3*NodeV(i);%deltaP2(i)=(deltaP2(i)-(y1*VX(i)+y2*VY(i)/Bus(i).V;endif Bus(i).type=2.专业 .整理 .下载可编辑%y1=0;%y2=0;y3=0;y4=0;for k=1:busnumberif (B(i,k)=0|G(i,k)=0)%y1=y1+(G(i,k)*VX(k)-B(i,k)*VY(k);%y2=y2+(G(i,k)*VY(k

26、)+B(i,k)*VX(k);y3=y3+NodeV(k)*(G(i,k)*cos(Nodea(i)-Nodea(k)+B(i,k)*sin(Nodea(i)-Nodea(k);y4=y4+NodeV(k)*(G(i,k)*sin(Nodea(i)-Nodea(k)-B(i,k)*cos(Nodea(i)-Nodea(k);endenddeltaP(i)=deltaP(i)-y3*NodeV(i);%deltaP2(i)=(deltaP2(i)-(y1*VX(i)+y2*VY(i)/Bus(i).V;deltaQ(i)=deltaQ(i)-y4*NodeV(i);%deltaQ2(i)=(de

27、ltaQ2(i)-(y1*VY(i)-y2*VX(i)/Bus(i).V;endif maxdP<abs(deltaP(i);maxdP=abs(deltaP(i);end.专业 .整理 .下载可编辑if maxdQ<abs(deltaQ(i);maxdQ=abs(deltaQ(i);enddeltaP(i)=deltaP(i)/NodeV(i);deltaQ(i)=deltaQ(i)/NodeV(i);endend3.3 输出并计算结果由以上程序可以得出 ：（ 1）进行系统正常运行状态的潮流计算，求得 Vi (0) ：己知公式 Vi (0) 1- （ Zif /(Z ff +Z f)）因为金属性短路时 Zf 0 公式为 Vi (0) 1- （Zif /Z ff ）再根据公式可得 ，V1&(0)=0.2774 ； V2&(0) 0.6770 ； V3&(0) 0.6770 ； V4&(0) 0；V5&(0)0.1870 ； V6&(0) 0