资源描述：
《电力系统短路电流计算例题与程序.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、20电力系统短路电流计算例题与程序编写佘名寰本文用导纳矩阵求逆计算节点阻抗矩阵，运用复合序网络图计算各节点对称故障和不对称故障时短路电流、节点电压和各支路故障电流。2.1用阻抗矩阵计算短路的基本公式：⑴节点三相对称短路，注入节点的短路电流Id=-Vd(0)/Zdd(2-1)式中Vd(0)故障点在短路发生前的电压，简化计算时Vd(0)=1Zdd故障点d的自阻抗负号表示电流从节点流出故障点短路电流在各节点所产生的电压分量V=ZI(2-2)式中Z节点阻抗矩阵I节点注入电流的列矩阵当只有一点故障时上述电压分量
2、为Vi(d)=ZdiId(i=1,2,3,………n)(2-3)式中Zdi故障点d与节点i的互阻抗短路故障后的节点电压Vi=Vi(0)+Vi(d)(2-4)式中VI(0)节点i故障发生前的电压短路故障时通过各支路的电流Iij=(Vi-VJ)/zij（2-5）式中zij联系节点i和节点j的支路阻抗⑵单相接地短路故障点的电流和电压：A相单相接地故障Ia0=Ia1=Ia2=-Va1(0)/(Zdd0+Zdd1+Zdd2)(2-6)Zdd0,Zdd1,Zdd2-------零序、正序、负序网络故障节点的自阻抗V
3、a0=Zdd0Ia0(2-7)Va1=Va1(0)+Zdd1Ia1(2-8)Va2=Zdd2Ia2(2-9)Ia=3Ia1(2-10)⑶两相接地短路：B．C相短路接地故障增广正序网的综合等值阻抗Z∑Z∑=Zdd0Zdd2/(Zdd0+Zdd2)(2-11)Ia1=-Va1(0)/(Zdd1+Z∑)(2-12)Ia0=-Ia1Zdd2/(Zdd0+Zdd2)(2-13)Ia2=-Ia1Zdd0/(Zdd0+Zdd2)(2-14)Ib=Ia0+a2Ia1+aIa2(2-15)a=(-1/2+j√3/2)2
4、0a2=(-1/2-j√3/2)⑷两相短路：B.C两相短路故障Ia1=Ia2=-Va1(0)/(Zdd1+Zdd2)(2-18)Ib=j√3Ia1(2-19)⑸支路i～j间的某一点d发生故障时，视d点为新的节点d点与节点k的互阻抗ZdkZdk=(1-L)ZIK+LZjk(2-20)d点的自阻抗ZddZdd=(1-L)2Zii+L2Zjj+2L(1-L)ZIJ+L(1-L)zij(2-21)式中L为端点i到故障点d的距离所占线路全长的百分数ZIK，Zjk分别为节点i和节点j与节点k的互阻抗Zii，,Zj
5、j为节点i和节点j的自阻抗ZIJ为节点i与节点j的互阻抗zij是节点i和节点j间的线路阻抗2.2短路电流计算时用导纳矩阵求逆计算节点阻抗矩阵参考文献①介绍了从网络的原始阻抗矩阵求节点导纳矩阵的方法和相关程序。与潮流计算时不同的是网络图参考节点为发电机中性点，计及发电机次暂态电抗；零序网络图参考节点为变压器中性点，计及变压器零序电抗，双回线有零序互感阻抗。包括中性点的节点导纳矩阵是不定导纳矩阵，在求逆矩阵前要去掉参考节点（中性点）转化为定导纳矩阵后再求阻抗矩阵。以下举例说明从导纳矩阵计算节点阻抗矩阵的方
6、法和程序。【例2.1】图2-1所示为一个3节点网络，两台发电机，四条线路。发电机一台中性点接地，另一台不接地。发电机次暂态电抗和线路阻抗原始数据在表2-1、表2-2中给出，线路3、4间有零序互感电抗。试计算该网络图的正序和零序节点阻抗矩阵。①①②③23451G1图2-1例2.1网络图G④6表2-1正序网络数据首端节点编号末端节点编号回路编号自感标么阻抗R(pu.)自感标么阻抗X(pu.)4160.00.20004132120.00.00.16000.0800202213333450.00.00.00.
7、06000.06000.1300表2-2零序网络数据首端节点编号末端节点编号回路编号自感标么阻抗R自感标么阻抗X互感标么阻抗RM互感标么阻抗XM4122132333123450.00.00.00.00.00.02000.14000.10000.12000.17000.00.00.05000.0500[例2.1]节点正序负序和零序阻抗矩阵计算源程序clearglobalNpNbNp=4;Nb=6;Np1=3;nstart=[1,2,2,1,4,4];nend=[2,3,3,3,1,3];mm=ffm(n
8、start,nend);%zb1zb0,theseriesimpedancesoftransmissionlinezb1i=[0.08,0.06,0.06,0.13,0.2,0.16];zb1=j*zb1i;yb1=zb1.^(-1);ys=diag(conj(yb1'));y=mm*ys*(mm)';Nb=5;nstart0=[1,2,2,1,4];nend0=[2,3,3,3,3];mm0=ffm(nstart0,nend0);zb0i=[0.14,
