matlab实验 电力系统短路分析(10页).doc

《matlab实验 电力系统短路分析(10页).doc》由会员分享，可在线阅读，更多相关《matlab实验 电力系统短路分析(10页).doc（10页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、-matlab实验 电力系统短路分析-第 - 17 - 页实验二 短路电流计算程序的实现一、三相短路电流计算程序计算短路电流周期分量，如（）时，实际上就是求解交流电路的稳态电流，其数学模型也就是网络的线性代数方程，一般选用节点电压方程。方程的系数矩阵是对称的。在短路电流计算中变化的量往往是方程的常数项，需要多次求解线性方程组。1等值网络图2-1给出了不计负荷情况下计算短路电流的等值网络。在图2-1（a）中G代表发电机端电压节点，发电机等值电势和电抗分别为和，D表示负荷节点，点为直接短路点。应用叠加原理如图2-1所示。正常运行方式为空载运行，网络中各点电压均为1；在故障分量网络中。只需作故障分量

2、的计算。由图2-1的故障分量网络可见，这个网络与潮流计算的网络的差别在于发电机节点上多接了对地电抗。当然如果短路计算中可以忽略线路电阻和电纳，而且不计变压器的实际变比，则短路计算网络较潮流计算网络简化，而且网络本身是纯感性的。图2-1 在不计负荷情况下计算短路电流I的等值电路2 用节电阻抗矩阵计算短路电流如果已经形成了故障分量网络的节点阻抗矩阵，则矩阵中的对角元素就是网络从点看进去的等值阻抗，又称为点的自阻抗。为点与点的互阻抗，均用大写表示。由节点方程中的第个方程：。为其它节电电流为零时，节点的电压和电流之比，即网络对点的等值阻抗。 根据故障分量网络，直接应用戴维南定理可求得直接短路电流（由故

3、障点流出）为 （2-1）式中，为接地阻抗；为点短路前的电压。如果短路点为直接短路，则=0，在实用计算中采用（2-2）式 （2-2）因此，一旦形成了节点阻抗矩阵，任一点的短路电流即可方便地求出，即等于该点自阻抗（该点对角元素）的倒数。节点导纳矩阵的特点是易于形成，当网络结构变化时也容易修改，而且矩阵本身是很稀疏的，但是应用它计算短路电流不如用节点阻抗矩阵那样直接。由于节点阻抗矩阵是节点导纳矩阵的逆矩阵，可以先求再求（等于），或者中的部分元素。具体计算可以采用以下步骤：（1） 应用计算短路点的自阻抗和互阻抗。（2） 应用（2-1）式计算短路电流。3计算节点电压和支路电流由故障分量网络可知，只有节点

4、有节点电流，各节点电压的故障分量为= (2-3)所以,，各节点短路故障后的电压为 (2-4)任一支路-的电流为 (2-5)式中，为-支路的阻抗。这种计算方法实际上就是利用节点导纳矩阵一次求得与故障点有关的一列节点阻抗矩阵元素，应用节点导纳矩阵计算短路电流的原理框图如图2-2所示。图2-2 应用节点导纳矩阵计算短路电流的原理框图例2-1 某三节点简单电力系统的简化等值电路如图2-3所示，阻抗参数标幺值已经标在图上，发电机电压近似认为是1。应用计算机算法计算节点3三相短路电流及各节点电压和各支路电流。图2-3 某电力系统的简化等值电路解：不计负荷的影响。下面给出计算步骤和源程序。计算步骤：（1）形

5、成节点导纳矩阵；（2）因为，所以。取，即节点3注入单位电流，求得电压向量，即节点3的自阻抗和互阻抗。（3）应用（2-2）式计算短电流；（4）应用（2-4）式计算各点电压；（5）应用（2-5）式计算线路故障电流。源程序：clearZZ(1,2)=j*0.1; ZZ(1,3)=j*0.1; ZZ(2,3)=j*0.1;%节点i,j之间的阻抗（ij）YB=-j*26.6266 j*10 j*10 j*10 -j*33.2933 j*10 j*10 j*10 -j*19.96;%输入节点导纳矩阵 n=3;%输入网络的节点数 k=3;%确定短路点的节点号 for i=1:n if i=k II(i)=1

6、; else II(i)=0; end endZ(:,k)=YBII; Zk=Z(:,k)%节点m的自阻抗和互阻抗 k,Ik=1/Z(k,k) for i=1:n U(i)=1-Z(i,k)*Ik; end Un=U for i=1:n for j=1:n if ij I(i,j)=(U(i)-U(j)/ZZ(i,j);%支路电流的实用计算 ij(1)=i;ij(2)=j; ij,Iij=I(i,j) end end end二、不对称短路故障的计算程序（一） 不对称短路故障的计算步骤1 近似的实用计算中，对于短路故障可假设各节点短路前瞬间电压均为1。如果要求准确计算故障前的运行情况，则需要进行

7、潮流计算。2 成正序、负序和零序节点导纳矩阵。发电机的正序电抗用，可计算故障后瞬时的量。发电机的负序电抗近似等于。当不考虑负荷影响时，在正、序负序网络不接入负荷阻抗。因为负荷的中性点一般不接地，所以零序无通路。3 形成三个序网的节点导纳矩阵后，可求得故障端点的等值阻抗。对于短路故障，只要令（其余节点电流均为零），分别应用三个序网的节点导纳矩阵求解一次即可得到三个序网和点的有关阻抗。4 根据不同的故障，分别利用表2-1列出的公式计算故障处各序电流、电压，进而合成得到三相电流、电压。表2-1 三种不对称短路在短路点处的各序电流、电压计算公式短路类型短路点各序电流计算公式短路点各序电压计算公式单相短

8、路两相短路两相短路接地同单相接地5 计算网络中任一点的电压，将用到以下相应的计算公式。 （2-6）式中，为短路前点的电压。6对于短路故障，任一支路的各序电流均可用下式计算： （2-7）将各序分量合成相量的问题，涉及到计算点和故障点之间的变压器的连接方式，在例7-10的源程序中有所反映。（二）计算原理框图应用对称分量法计算不对称短路故障的计算步骤是很简明的。图2-4给出计算短路故障的计算程序原理框图。图2-4 不对称短路计算程序框图（三）例题例2-2 编写程序计算例2-1中节点3发生单相短路接地、两相短路的瞬时，（1）节点1和2的电压；（2）线路1-2、1-3和2-3的电流；（3）发电机1、2

9、的端电压。解：瞬时的含义：发电机的正序电抗用，负序电抗近似等于。计算结果如下表2-2 节点电压节点1电压节点2电压节点3电压单相接地=0两相短路表2-3 支路电流支路12支路13支路23单相接地=0.6613两相短路=0=0=0表2-4 发电机的机端电压发电机1#发电机2#单相接地两相短路源程序：clearZZ1(1,2)=j*0.1; ZZ1(1,3)=j*0.1; ZZ1(2,3)=j*0.1;%节点m,n之间的正序阻抗（mn）ZZ2(1,2)=j*0.1; ZZ2(1,3)=j*0.1; ZZ2(2,3)=j*0.1;%节点m,n之间的负序阻抗（mn）ZZ0(1,2)=j*0.2; ZZ

10、0(1,3)=j*0.2; ZZ0(2,3)=j*0.2;%节点m,n之间的零序阻抗（mn）Y1=-j*26.6266 j*10 j*10 j*10 -j*33.2933 j*10 j*10 j*10 -j*19.96;%输入正序网络节点导纳矩阵Y2=-j*26.6266 j*10 j*10 j*10 -j*33.2933 j*10 j*10 j*10 -j*19.96;%输入负序网络节点导纳矩阵Y0=-j*30 j*5 j*5 j*5 -j*50 j*5 j*5 j*5 -j*10;%输入零序网络节点导纳矩阵 YY1=-j*39.96 j*10 j*10 j*20 0 j*10 -j*59.

11、96 j*10 0 j*40 j*10 j*10 -j*19.96 0 0 j*20 0 0 -j*30 0 0 j*40 0 0 -j*60; YY2=YY1;%输入包括发电机机端电压节点的正，负序网络节点导纳矩阵 N1=3;%输入网络的节点数 N2=5;%输入包括所有发电机节点的网络的节点数 k=3;%输入短路点的节点号 fault=1;%输入短路类型f(3)=3;f(1)=1;f(2)=2;f(1,1)=4 %第一部分：计算所有节点的a,b,c三相电压 for p=1:N1 if p=k I(p)=1; else I(p)=0; end end Z1(:,k)=Y1I;Zk1=Z1(:,

12、k);%正序网络中节点m的自阻抗和互阻抗 Z2(:,k)=Y2I;Zk2=Z2(:,k);%负序网络中节点m的自阻抗和互阻抗 Z0(:,k)=Y0I;Zk0=Z0(:,k);%零序网络中节点m的自阻抗和互阻抗 if fault=1%根据故障类型选择不同的计算公式 Ik1=1/(Z1(k,k)+Z2(k,k)+Z0(k,k); Ik2=Ik1;Ik0=Ik1; else if fault=2 Ik1=1/(Z1(k,k)+Z2(k,k); Ik2=-Ik1;Ik0=0; else if fault=3 Ik1=1/Z1(k,k);Ik2=0;Ik0=0; else if fault=4 Ik1=

13、1/(Z1(k,k)+Z2(k,k)*Z0(k,k)/(Z2(k,k)+Z0(k,k); Ik2=-Ik1*Z0(k,k)/(Z2(k,k)+Z0(k,k); Ik0=-Ik1*Z2(k,k)/(Z2(k,k)+Z0(k,k); end end end end Ik1 %计算短路节点的正序电流 for p=1:N1 if p=k I1(p)=-Ik1; I2(p)=-Ik2; I0(p)=-Ik0; else I1(p)=0; I2(p)=0; I0(p)=0; end end uu1(:,k)=Y1I1.; uu2(:,k)=Y2I2.; uu0(:,k)=Y0I0.; for p=1:N1

14、 U1(p)=1; end u1=U1+uu1(:,k);%计算所有节点正序电压 u2=uu2(:,k); %计算所有节点负序电压 u0=uu0(:,k); %计算所有节点零序电压 a=-0.5+j*sqrt(3)/2; T=1 1 1 a2 a 1 a a2 1; for p=1:N1 U=u1(p) u2(p) u0(p); p Uabc=T*U.; %T为对称分量法的合成矩阵 UUabc=abs(Uabc)%UUabc表示i节点的a,b,c三相电压有效值 end%第二部分：计算支路电流 for p=1:N1 U1(p)=1; end u1=U1+uu1(:,k);%计算所有节点正序电压

15、u2=uu2(:,k);%计算所有节点负序电压 u0=uu0(:,k);%计算所有节点零序电压 for m=1:N1 for n=1:N1 if mn mn(1)=m;mn(2)=n; mn I1(m,n)=(u1(m)-u1(n)/ZZ1(m,n);%正序支路电流的实用计算 I2(m,n)=(u2(m)-u2(n)/ZZ2(m,n);%负序支路电流的实用计算 I0(m,n)=(u0(m)-u0(n)/ZZ0(m,n);%零序支路电流的实用计算 Iabc=T*I1(m,n) I2(m,n) I0(m,n).; Iabc%Iabc表示支路(m,n)的a,b,c三相电流 abs(Iabc) end

16、 end end %第三部分：计算发电机的端电压 for p=1:N2 if p=k II(p)=-Ik1; else II(p)=0; end end vv1(:,k)=YY1II.; vv2(:,k)=YY2II.; for p=1:N2 V1(p)=1; end v1=V1+vv1(:,k);v2=vv2(:,k); v0=0; a1=sqrt(3)/2+j*0.5;a2=sqrt(3)/2-j*0.5;a0=0; for m=N1+1:N2 m Vabc=T*(v1(m) v2(m) v0.*a1 a2 a0).;%考虑到变压器为Y/-11接线 VVabc=abs(Vabc)%VVabc表示发电机机端a,b,c三相电压的有效值 end