电力系统课程设计概要PPT学习教案

1、会计学1电力系统课程设计概要电力系统课程设计概要第1页/共42页第2页/共42页第3页/共42页第4页/共42页第5页/共42页第6页/共42页dXdXdX第7页/共42页line = 1 2 0.04 0.25 0.0 0.25 0; 1 3 0.10 0.35 0.0 0.0 0; 2 3 0.08 0.30 0.0 0.25 0; 5 3 0.00 0.03 0.0 0.0 1.05; 4 2 0.00 0.015 0.0 0.0 1.05;第8页/共42页第9页/共42页% (bus#)( volt )( ang )( p )( q )(bus type)(是否机组1是0否)(Xd)b

2、us= 1 1.00 0.0 -1.60 -0.80 1 0 0.000 ; 2 1.00 0.0 -2.00 -1.00 1 0 0.000 ; 3 1.00 0.0 -3.70 -1.30 1 0 0.000; 4 1.00 0.0 5.00 0.00 2 1 0.020; 5 1.00 0.0 0.00 0.00 3 1 0.020;% b#1 b#2 ( R )( X )( G )( B/2 ) ( tap )line = 1 2 0.04 0.25 0.0 0.25 0; 1 3 0.10 0.35 0.0 0.0 0; 2 3 0.08 0.30 0.0 0.25 0; 5 3 0

3、.00 0.03 0.0 0.0 1.05; 4 2 0.00 0.015 0.0 0.0 1.05;% 故障支路首端节点 故障支路末端短路 （距首端距离0-1）fault = 1 2 0.5; 第10页/共42页第11页/共42页第12页/共42页1 :1ijTymYTy2111Ty111TyijmY Y(I,I)=Y(I,I)+Yt+Ym;Y(I,I)=Y(I,I)+Yt+Ym; Y(J,J)=Y(J,J)+Yt; Y(J,J)=Y(J,J)+Yt; Y(I,J)=Y(I,J)-Yt; Y(I,J)=Y(I,J)-Yt; Y(J,I)=Y(I,J); Y(J,I)=Y(I,J); 第13页

4、/共42页n潮流计算的潮流计算的n支路参数：支路参数：nZt=R+jXZt=R+jX；nYt=1/Zt;Yt=1/Zt;nYm=G+jB/2;Ym=G+jB/2;n短路计算的短路计算的n近似支路参数近似支路参数1 1：nZt=R+jXZt=R+jX；nYt=imag(1/Zt);Yt=imag(1/Zt);nYm=1/(G+jB/2);Ym=1/(G+jB/2);n短路计算的短路计算的n近似支路参数近似支路参数2 2：nZt=jXZt=jX；nYt=1/Zt=1/jX;Yt=1/Zt=1/jX;nYm=jB/2;Ym=jB/2;第14页/共42页n短路计算的近似支路参数短路计算的近似支路参数1

5、 1：nZt=R+jXZt=R+jX；nYt=imag(1/Zt);Yt=imag(1/Zt);nYm=1/(G+jB/2);Ym=1/(G+jB/2);n短路计算的近似支路参数短路计算的近似支路参数2 2：nZt=jXZt=jX；nYt=1/Zt=1/jX;Yt=1/Zt=1/jX;nYm=jB/2;Ym=jB/2;第15页/共42页第16页/共42页v发电机作为有源支路通常表示为电势源Ei与阻抗zi （即jXd ）的串联支路（机端节点i与零电势点之间的“支路”），因此在短路电流计算时，使用等效电流源Ii与对地导纳yi来表示。/1/iididIEjXyjX第17页/共42页v等效变换后，机端

6、节点i增加了一条新的对地支路，因此需要对节点导纳矩阵对角线上第i个元素进行修正。1/ijijdfixyyjX第18页/共42页fixY = Y;for k=1:nb busType=bus(k,7); if (busType=1) fixY(bus(k,1),bus(k,1) = fixY(bus(k,1),bus(k,1) + 1/j/bus(k,8); endend第19页/共42页第20页/共42页ijyijiiyjjyjkykkkyijiiijyykkkjkyy接地第21页/共42页ijyijiiyjjy jkykkky接地第22页/共42页ijyijiiyjjy jkykkky接地9

7、99999999999jjyj第23页/共42页接地ijyijiiyjjyijiiyjjy0.5ijy0.5ijy第24页/共42页/0.5iiiiijyyy/0.5jjjjijyyy0ijjiyy第25页/共42页nf,mf= size(fault);for k=1:nf nodeI=fault(k,1); nodeJ=fault(k,2); dis=fault(k,3); if (nodeI=0) fixY(nodeJ,nodeJ) = 999999+j*999999; continue; end if (nodeJ=0) fixY(nodeI,nodeI) = 999999+j*9999

8、99; continue; end if (dis=0)&(nodeI*nodeJ=0) fixY(nodeI,nodeI) = 999999+j*999999; continue; end if (dis=1)&(nodeI*nodeJ=0) fixY(nodeJ,nodeJ) = 999999+j*999999; continue; end%接下页第26页/共42页%接上页if (dis=1)&(dis=0)&(nodeI*nodeJ=0) fixY(nodeI,nodeI) = fixY(nodeI,nodeI) - fixY(nodeI,nodeJ)/dis; fixY(nodeJ,n

9、odeJ) = fixY(nodeJ,nodeJ) - fixY(nodeI,nodeJ)/(1-dis); fixY(nodeI,nodeJ)=0; fixY(nodeJ,nodeI)=0; end end第27页/共42页第28页/共42页v节点注入电流 为N维列向量，其中N为电力系统节点数。将发电机从电势源等效变换为电流源后，可求出故障时发电机节点的注入电流，进而生成节点注入电流列向量/iidIEjXII第29页/共42页v在短路电流实用计算中，可以近似认为所有发电机的电势都为1，因此对公式作如下变换1/idIjX/iidIEjX第30页/共42页Iinj = zeros(nb,1);

10、for k=1:nb busType=bus(k,7); if(calcSettings(1)=1) v = 1; else v = bus(k,2); end if (busType=1) Iinj(bus(k,1),1) = Iinj(bus(k,1),1) + v/j/bus(k,8); endend第31页/共42页IYU第32页/共42页第33页/共42页v短路电流计算结果可以用一个短路电流计算结果可以用一个 的矩阵的矩阵表示，其中表示，其中nlnl为原始数据中的支路数量。为原始数据中的支路数量。矩阵的前两列分别为首末端节点编号，第三列为矩阵的前两列分别为首末端节点编号，第三列为短路

11、电流相量，第四列为短路电流有效值；短路电流相量，第四列为短路电流有效值；矩阵的前矩阵的前nlnl行为与原始数据中行为与原始数据中lineline矩阵对应的支矩阵对应的支路及其短路电流值；路及其短路电流值；矩阵的最后一行为故障点短路电流值。矩阵的最后一行为故障点短路电流值。14nl bI第34页/共42页ijyijiiyjjy jkykkky接地jijkjijijjkkjIIIyUUyUU第35页/共42页faultNode = 0;if(faultI=0) faultNode = faultJ;endif(faultJ=0) faultNode = faultI;endif(dis=1)&(f

12、aultI*faultJ=0) faultNode = faultJ;endif(dis=0)&(faultI*faultJ=0) faultNode = faultI;endif(faultNode=0) Ib(nl+1,1) = faultNode; Ib(nl+1,2) = faultNode; Ifault = 0; for k=1:nl I=line(k,1); J=line(k,2); Ib(k,1) = I; Ib(k,2) = J; if (I=faultNode) Ifault = Ifault - Ib(k,3); else if (J=faultNode) Ifault

13、= Ifault + Ib(k,3); end end end Ib(nl+1,3) = Ifault; Ib(nl+1,4) = abs(Ib(nl+1,3);end获取故障节点编号前两列都显示故障节点编号Ib(k,3)为已求得的非故障相关支路短路电流第36页/共42页接地ijiiyjjy0.5ijy0.5ijy0.50.5ijijijijyyIUU第37页/共42页if(dis=0)&(dis=1)&(faultI*faultJ=0) Ib(nl+1,1) = faultI; Ib(nl+1,2) = faultJ; Ib(nl+1,3) = U(faultI)*Y(faultI,faultJ)/dis + U(faultJ)*Y(faultI,faultJ)/(1-dis); Ib(nl+1,4) = abs(Ib(nl+1,3)