短路电流计算的matlab算法

1、河南城建皇家学院电力系统短路计算的 MATLAB算法姓 名：学 号：专业班级：指导老师：所在院系： 电气与信息工程学院2014年11月22日摘要本设计是利用设计相关的MATLA程序实现对不同类型电力系统短路故障进 行计算机计算。随着人类生产生活对电力日加依赖， 电网正向着跨国界的巨型电网发展， 随 着电网规模的扩大， 人类对电力系统的稳定性要求也日益提高。 电力系统的短路 故障是电网故障中较为严重的故障， 而且是发生几率最多的故障。 当发生短路时， 其短路电流可达数万安培， 巨大的短路电流产生的热效应和电动力效应将使电气 设备遭受严重破环， 所以当发生短路时， 开关电气设备必须经得起可能的最大

2、短 路电流而不致损坏， 所以求取相关网络的短路电流对于电网的设计具有不可估量 的巨大作用。本文所设计的MATLA算法是根据电力系统发生短路时的相关特点，建立不同 故障所对应的短路数学模型， 从而构造出网络的线性代数方程， 最终得到电力系 统短路电流的MATLA汁算方法一一节点阻抗矩阵的支路追加法， 该办法适用于各 种结构的电网，展现出了计算机计算的巨大优势。目录1 引言 42 理论分析 52.1 计算条件 52.2 计算步骤 53 仿真分析 83.1 程序主框图 83.2 程序代码 93.3 仿真 154 结论 205 参考文献 211 引言 电力系统短路电流计算是电力系统运行、 分析的重要环

3、节， 是电力设计中最 重要的计算之一。 传统的短路电流计算是以手工计算为基础进行的， 先通过相关 电路知识化简所求的电力网络 , 求出各电源点对短路点的转移阻抗，进而计算出 电抗XJS,再查找运算曲线，最终求得短路电流的周期分量。传统的手工计算过 程非常繁杂 ,工作量大 , 容易出错。随着电网规模的扩大传统的手工计算已经不能满足现代电力网络设计的要求。此时，采用计算机辅助计算显得势在必行。本文所设计的 MATLA算法是根据 电力系统发生短路时的相关特点， 建立不同故障所对应的短路数学模型， 从而构 造出网络的线性代数方程，最终得到电力系统短路电流的MATLA汁算方法节 点阻抗矩阵的支路追加法。

4、j0.152理论分析本文针对以下例题进行详细的计算机算法的分析和设计。例：某三节点电力系统的简化等值电路如图2-1所示，阻抗参数标幺值编注在图上，发电机电压认为是1。计算：1. 节点3三相短路电流及各节点电压和各支路电流。2. 节点3发生单相短路接地、两相短路的瞬时，（1）节点1和2的电压；（2）线路1-2、1-3和2-3的电流；（3）发电机1、2的端电压。j0.075图2-1某电力系统的简化等值电路2.1计算条件1）假设系统有无限大的容量。短路后，系统母线电压能维持不变。即计算阻抗比系统阻抗要大得多。2）在计算高压电器中的短路电流时，只需考虑发电机、变压器、电抗器的电抗，而忽略其电阻；对于架

5、空线和电缆，只有当其电阻大于电抗1/3时才需计入电阻，一般也只计电抗而忽略电阻。3）短路电流计算公式或计算图表，都以三相短路为计算条件。因为单相短路或 二相短路时的短路电流都小于三相短路电流。2.2计算步骤三相短路：1. 用节电阻抗矩阵计算短路电流如果已经形成了故障分量网络的节点阻抗矩阵，则矩阵中的对角元素就是网 络从f点看进去的等值阻抗，又称为f点的自阻抗。乙为f点与i点的互阻抗均 用大写Z表示。由节点方程中的第f个方程：Uf Zf1I1 ZffI fzfnln。Zff为其它节电电流为零时，节点f的电压和电流之比，即网络对f点的等 值阻抗。根据故障分量网络，直接应用戴维南定理可求得直接短路电

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

7、和互阻抗 乙彳Zf Znf。应用（2-1 ）式计算短路电流。2 计算节点电压和支路电流由故障分量网络可知，只有节点f有节点电流I f，各节点电压的故障分量 为U f :=Z fiZffZfnI f :0=Z ffI fUnZ n1ZnfZ nnZnfU 1Z11Z1 f0 乙f(2-3)所以，各节点短路故障后的电压为U i U io AU i U io Zi f 11Uf Uf|oAU f 0U nU 屮|AU nU 0Znf 1 f(2-4)任一支路i - j的电流为I.UiUj(2-5)Zij式中，Zj为i- j支路的阻抗。不对称短路：1. 近似的实用计算中，对于短路故障可假设各节点短路前

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

9、称短路在短路点处的各序电流、电压计算公式短路类型短路点各序电流计算公式短路点各序电压计算公 式单相短路I f (1)I f(2)I (0)Uf|0U f (1)U ”0I f (1) Zff (1)U f (2)I f (2) Zff (2)U f (0)I f (0) Zff (0)Z ff (1)Z ff (2)Z ff (0)两相短路IfIUqqU f(1) Uf|0I f (1)Zff (1)U f (2)I f (2)Zff (2)2)Zff ZffIUf|0|两相短路接地f(1)ZZff(1)IZ ff(2)Z ff(0)/(Zff(2)Zff (0)IZ ff (0)同单相接地1

10、 f(2)If(1) zz厶 ff (2)厶 ff (0)IZff Ifff(1) ZZff (2) Jff (0)5. 计算网络中任一点的电压，将用到以下相应的计算公式Ui(1)U10|I f (1)Zif (1)UiI f (2)Zif (2)Ui(0)I f (0)Zif (0)(2-6 )式中，5o为短路前i点的电压6对于短路故障，任一支路的各序电流均可用下式计算:3.1程序主框图I ij (1)I ij (2)I ij (0)Ui(1)Uj (1)Z ij (1)Ui(2)Uj (2 )Z ij (2)Ui(0)Uj (0 )Z ij (0)3仿真分析(2-7 )图3.1程序框图3.

11、2程序代码三相短路：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*10j*10j*10-j*33.2933 j*10j*10j*10-j*19.96;%俞入节点导纳矩阵n=3; %俞入网络的节点数k=3; %确定短路点的节点号for i=1:nif i=kII(i)=1;elseII(i)=0;endendZ(:,k)=YBII;Zk=Z(:,k)%节点m的自阻抗和互阻抗k,Ik=1/Z(k,k)for i=1:nU(i)=1-Z(i,k)*Ik;endUn=Ufor i=1

12、:nfor 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)endendEnd不对称短路：clearZZ1(1,2)=j*0.1; ZZ1(1,3)=j*0.1;ZZ1(2,3)=j*0.1;%节点m,n之间的正序阻抗(mn%节点m,n之间的负序阻抗(mnZZ0(1,2)=j*0.2; ZZ0(1,3)=j*0.2;ZZ0(2,3)=j*0.2;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.6

13、266 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 j*10 -j*59.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;%节点m,n之间的零序阻抗(mn%输入正序网络节点导纳矩阵%输入负序网络节点导纳矩阵%输入零序网络节点导纳矩阵YY2=YY1;%俞入包括发电机机端电压节点的正，负序网络节点导纳矩阵

14、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=kZZ2(1,2)=j*0.1; ZZ2(1,3)=j*0.1; ZZ2(2,3)=j*0.1;I(p)=1; elseI(p)=0;endendZ1(:,k)=Y1l;Zk仁Z1(:,k);%E序网络中节点m的自阻抗和互阻抗Z2(:,k)=Y2l;Zk2=Z2(:,k);%负序网络中节点m的自阻抗和互阻抗ZO(:,k)=Y

15、Ol;ZkO=ZO(:,k);%?序网络中节点m的自阻抗和互阻抗if fault=1 %根据故障类型选择不同的计算公式lk1=1/(Z1(k,k)+Z2(k,k)+ZO(k,k);lk2=lk1;lkO=lk1;elseif fault=2lk1=1/(Z1(k,k)+Z2(k,k); lk2=-lk1;lkO=O;elseif fault=3lk1=1/Z1(k,k);lk2=O;lkO=O;elseif fault=4lk1=1/(Z1(k,k)+Z2(k,k)*ZO(k,k)/(Z2(k,k)+ZO(k,k);lk2=-lk1*ZO(k,k)/(Z2(k,k)+ZO(k,k);lkO=-

16、lk1*Z2(k,k)/(Z2(k,k)+ZO(k,k);endendendendlk1%计算短路节点的正序电流for p=1:N1if p=kI(p)=1;elseI(p)=0;endendZ1(:,k)=Y1l;Zk仁Z1(:,k);%E序网络中节点m的自阻抗和互阻抗Z2(:,k)=Y2l;Zk2=Z2(:,k);%负序网络中节点m的自阻抗和互阻抗ZO(:,k)=YOl;ZkO=ZO(:,k);%?序网络中节点m的自阻抗和互阻抗%第二部分：计算支路电流for p=1:N1U1(p)=1;endu1=U1+uu1(:,k);%计算所有节点正序电压u2=uu2(:,k); %计算所有节点负序电

17、压uO=uuO(:,k); %计算所有节点?序电压for m=1:N1for n=1:N1if mnmn(1)=m;mn(2)=n;mnl1(m,n)=(u1(m)-u1(n)/ZZ1(m,n);%正序支路电流的实用计算l2(m,n)=(u2(m)-u2(n)/ZZ2(m,n);%负序支路电流的实用计算lO(m,n)=(uO(m)-uO(n)/ZZO(m,n);%?序支路电流的实用计算labc=T*l1(m,n) l2(m,n) lO(m,n).;labc%labc表示支路(m,n)的a,b,c三相电流abs(labc)endendend%第三部分：计算发电机的端电压for p=1:N2if

18、p=kII(p)=-Ik1;elseII(p)=0;endendvv1(:,k)=YY1II.;vv2(:,k)=YY2II.;for p=1:N2V1(p)=1;endv1=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:N2mVabc=T*(v1(m) v2(m) vO.*a1 a2 aO).;%考虑到变压器为 Y/ -11 接线VVabc=abs(Vabc)%VVab表示发电机机端a,b,c三相电压的有效值end3.3 仿真 三相短路min iels IIhQjutnd

19、 1（1 h NzYBltf ZksZ! I # 30 节苗沪白Hfi 魏如互蘭林 比 Ilr=I/2 g k）* ma_* 0.。” i0-1 73.I49KO. -4 0*19图3.2三相短路运行图w-忑 IL -q IP$4目(d 二2g -k s. 二啊I s- 用代rii z&pr oll30社冠护*#S2-03A3 出制吐.%FCS龙我 -2*3sfl3 因雷址H比护ff*t章二戈二t*二? z 口E 二IIA.L0冲*dMoA2 二 0弓上*二FU护丄上忙二*：9M 口3 EM 4*5 laz V-:E$J/【*MXOSA a _- s- * J HN s * If s H/ |

20、耳躍 h LIM-pu-QWAnUDrt fn 片吕 Ku 武更 xf #-. *r1w 半电MRgE嗟岂xtB筝Fs E4- FSSFCQ- tt- -9 el*R 乜 Q eft -Q7IT Qp 吗巧 ooa-l*lr s Gil;!- oo7r 0.r 沖Ou w0ronUIIacMAiJ?菲 Sh-TII 症总 Bp d g:H2 O-M*f m_Ha Hd置 guj哥 Hu .0 匸章二R )ffNM- *cvl* 髯 HHMUUMudtFyIEC巧-xflJ h;e rljmWOHN EN rplFH&NtaWWNIE*cftHW rdr-lE h-Nv啊- (u2 () -u2

21、 (n3 )/112 (, nl 站轉盍殆电庭的WRJ14K IQ11(斗M I2(vun:聽示妄略“皿的mg亡三相宙 廉rvdcn-dend苒稱三譚討i计xitvuniinrntffar p1rW2if Flsii j= n ；II lp)On4ervd,wl Giiii-TfYlMt. .watikli-YYSxIIu *fK pal iN2V| p)J ,endvlVl- +*vl /2-J*O S . ttC=QfairsJTiVaJbc= I 1 vi (k V3 &) ilO). (blI o2 b0八 * ;3l|*MftSSJt|T/ A* 11 W盂* 愛仏点打 电町机机曲叫S c三柑电住們眸0幢ndItaJl =0-3- | D&1UUaJbc Or,w4140. 9%14UUabc -o. -?Qr 3Z325mJ aib-c O *fl：5 1 3x0番0粛 二 D.丄 4 4 41.a- a 14441I abnr 壬0 - d 31fiTa aae - o. osiO0000 - Q04391L4C議： m oom脅6 0FOBI mib-d-a 仇 04001a aouit + _ oee