短路电流计算的matlab算法之欧阳理创编

河南城建皇家学院

时间：2021.03.05

创作：欧阳理 电力系统短路计算的MATLAB算法

姓 名： 学 号： 专业班级： 指导老师：

所在院系：电气与信息工程学院 2014年11月 22 日

欧阳阳理创编 2021.03.04

欧阳阳理创编 2021.03.04

摘要

本设计是利用设计相关的MATLAB程序实现对不同类型电力系统短路故障进行计算机计算。

随着人类生产生活对电力日加依赖，电网正向着界的巨型电网发展，随着电网规模的扩大，人类对电力系统的稳定性要求也日益提高。电力系统的短路故障是电网故障中较为严重的故障，而且是发生几率最多的故障。当发生短路时，其短路电流可达数万安培，巨大的短路电流产生的热效应和电动力效应将使电气设备遭受严重破环，所以当发生短路时，开关电气设备必须经得起可能的最大短路电流而不致损坏，所以求取相关网络的短路电流对于电网的设计具有不可估量的巨大作用。

本文所设计的MATLAB算法是根据电力系统发生短路时的相关特点，建立不同故障所对应的短路数学模型，从而构造出网络的线性代数方程，最终得到电力系统短路电流的MATLAB汁算方法──节点阻抗矩阵的支路追加法，该办法适用于各种结构的电网，展现出了计算机计算的巨大优势。

目录

1 引言2 2 理论分析4

2.1计算条件4 2.2计算步骤4 3 仿真分析7

3.1程序主框图7 3.2程序代码7 3.3仿真12 4 结论16 5参考文献18

1 引言

电力系统短路电流计算是电力系统运行、分析的重要环节，是电力设计中最重要的计算之一。传统的短路电流计算是以手工计算为基

欧阳阳理创编 2021.03.04

欧阳阳理创编 2021.03.04

础进行的，先通过相关电路知识化简所求的电力网络,求出各电源点对短路点的转移阻抗，进而计算出电抗XJS，再查找运算曲线，最终求得短路电流的周期分量。传统的手工计算过程非常繁杂,工作量大, 容易出错。

随着电网规模的扩大传统的手工计算已经不能满足现代电力网络设计的要求。此时，采用计算机辅助计算显得势在必行。本文所设计的MATLAB算法是根据电力系统发生短路时的相关特点，建立不同故障所对应的短路数学模型，从而构造出网络的线性代数方程，最终得到电力系统短路电流的MATLAB汁算方法──节点阻抗矩阵的支路追加法。

欧阳阳理创编 2021.03.04

欧阳阳理创编 2021.03.04

2 理论分析

本文针对以下例题进行详细的计算机算法的分析和设计。 例：某三节点电力系统的简化等值电路如图2-1所示，阻抗参数标幺值编注在图上，发电机电压认为是1。计算：

1.节点3三相短路电流及各节点电压和各支路电流。

2.节点3发生单相短路接地、两相短路的瞬时，（1）节点1和2的电压；（2）线路1-2、1-3和2-3的电流；（3）发电机1、2 的端电压。

1j0.15①1j0.075j0.1②j0.1③j0.1

图2-1 某电力系统的简化等值电路

2.1计算条件

1) 假设系统有无限大的容量。短路后，系统母线电压能维持不变。

即计算阻抗比系统阻抗要大得多。

2) 在计算高压电器中的短路电流时，只需考虑发电机、变压器、电

抗器的电抗，而忽略其电阻；对于架空线和电缆，只有当其电阻大于电抗1/3时才需计入电阻，一般也只计电抗而忽略电阻。 3) 短路电流计算公式或计算图表，都以三相短路为计算条件。因为

单相短路或二相短路时的短路电流都小于三相短路电流。 2.2计算步骤

三相短路：

1.用节电阻抗矩阵计算短路电流

如果已经形成了故障分量网络的节点阻抗矩阵，则矩阵中的对角元素就是网络从f点看进去的等值阻抗，又称为f点的自阻抗。Zfi为

f点与i点的互阻抗均用大写Z表示。

由节点方程中的第f个方程：

ZIUff11ZffIfZfnIn。

欧阳阳理创编 2021.03.04

欧阳阳理创编 2021.03.04

Zff为其它节电电流为零时，节点f的电压和电流之比，即网络对

f点的等值阻抗。

根据故障分量网络，直接应用戴维南定理可求得直接短路电流（由故障点流出）为

IfUf0Zffzf （2-1）

为f点短路前的电压。 式中，zf为接地阻抗；Uf0如果短路点为直接短路，则zf=0，在实用计算中采用（2-2）式

IfUf0Zff1 （2-2） Zff因此，一旦形成了节点阻抗矩阵，任一点的短路电流即可方便地求出，即等于该点自阻抗的倒数。

节点导纳矩阵的特点是易于形成，当网络结构变化时也容易修改，而且矩阵本身是很稀疏的，但是应用它计算短路电流不如用节点阻抗矩阵那样直接。由于节点阻抗矩阵ZB是节点导纳矩阵YB的逆矩阵，可以先求YB再求ZB（等于YB1），或者ZB中的部分元素。

具体计算可以采用以下步骤：

应用YB计算短路点f的自阻抗和互阻抗Z1fZffZnf。 应用（2-1）式计算短路电流。 2．计算节点电压和支路电流

由故障分量网络可知，只有节点f有节点电流If，各节点电压的故障分量为

Z11U1=Zf1UfUZn1nZ1fZffZnfZfnZnn0If0Z1f=ZffZnf (2-3) If所以，各节点短路故障后的电压为

欧阳阳理创编 2021.03.04

欧阳阳理创编 2021.03.04

UΔUUΖΙU111f11010U0UΔUff0fUUUΔUnnn0n0ΖnfΙf

(2-4)

任一支路i-j的电流为

IijUUIzijj

(2-5)

式中，zij为i-j支路的阻抗。 不对称短路：

1.近似的实用计算中，对于短路故障可假设各节点短路前瞬间电压均为1。如果要求准确计算故障前的运行情况，则需要进行潮流计算。

，2.成正序、负序和零序节点导纳矩阵。发电机的正序电抗用xd。当不考虑负可计算故障后瞬时的量。发电机的负序电抗近似等于xd荷影响时，在正、序负序网络不接入负荷阻抗。因为负荷的中性点一般不接地，所以零序无通路。

3.形成三个序网的节点导纳矩阵后，可求得故障端点的等值阻抗。对于短路故障，只要令If1（其余节点电流均为零），分别应用三个序网的节点导纳矩阵求解一次即可得到三个序网和f点的有关阻抗。

4.根据不同的故障，分别利用表2-1列出的公式计算故障处各序电流、电压，进而合成得到三相电流、电压。

表2-1 三种不对称短路在短路点处的各序电流、电压计算公式 短路类型 短路点各序电流计算公式 IIf(1)f(2)I(0)Uf0Zff(1)Zff(2)Zff(0)短路点各序电压计算公式 Uf(1)UfUf(2)Uf(0)0单相短路  ZIf(1)ff(1)IZf(2)ZIf(0)ff(0)ZIf(1)ff(1)ZIf(2)ff(2) 两相短路 If(1)If(2)Uf0Zff(1)Zff(2) Uf(1)UfUf(2)0 ff(2)欧阳阳理创编 2021.03.04

欧阳阳理创编 2021.03.04

UfIf(1)0两相短路接地 Zff(1)Zff(2)Zff(0)(Zff(2)Zff(0))Zff(0)If(2)If(1)Zff(2)Zff(0)Zff(2)IIfff(1)Zff(2)Zff(0) 同单相接地 5.计算网络中任一点的电压，将用到以下相应的计算公式。

UUi(1)i0Ui(2)Ui(0)ZIf(1)if(1)If(2)Zif(2) ZIf(0)if(0)（2-6）

为短路前i点的电压。 式中，Ui06．对于短路故障，任一支路的各序电流均可用下式计算：

Iij(1)Iij(2)Iij(0)Ui(1)Uj(1)Zij(1)Ui(2)Uj(2)Zij(2) Ui(0)Uj(0)Zij(0)（2-7）

3 仿真分析

3.1程序主框图

3.2程序代码 三相短路： clear

形成节点阻抗矩阵 输入数据 图3.1程序框图 ZZ(1,2)=j*0.1; ZZ(1,3)=j*0.1; ZZ(2,3)=j*0.1;%节点i,j之间的阻抗（i输入总共节点数及故障节点号 YB=[-j*26.6266 j*10 j*10 j*10 -j*33.2933 j*10 j*10 j*10 -j*19.96];%输入节点导纳矩阵 计算短路电流If n=3;%输入网络的节点数 计算节点电压 欧阳阳理创编 2021.03.04 计算节点电流 欧阳阳理创编 2021.03.04

k=3;%确定短路点的节点号 for i=1:n if i==k

II(i)=1; else II(i)=0; end end

Z(:,k)=YB\\II';

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 iI(i,j)=(U(i)-U(j))/ZZ(i,j);%支路电流的实用计算 ij(1)=i;ij(2)=j; ij,Iij=I(i,j) end end End 不对称短路： clear

ZZ1(1,2)=j*0.1; ZZ1(1,3)=j*0.1;

ZZ1(2,3)=j*0.1;%节点m,n之间的正序阻抗（mZZ2(2,3)=j*0.1;%节点m,n之间的负序阻抗（m欧阳阳理创编 2021.03.04

欧阳阳理创编 2021.03.04

ZZ0(1,2)=j*0.2; ZZ0(1,3)=j*0.2;

ZZ0(2,3)=j*0.2;%节点m,n之间的零序阻抗（mj*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.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

欧阳阳理创编 2021.03.04

欧阳阳理创编 2021.03.04

Z1(:,k)=Y1\\I';Zk1=Z1(:,k);%正序网络中节点m的自阻抗和互阻抗 Z2(:,k)=Y2\\I';Zk2=Z2(:,k);%负序网络中节点m的自阻抗和互阻抗 Z0(:,k)=Y0\\I';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=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

I(p)=1; else

I(p)=0; end

欧阳阳理创编 2021.03.04

欧阳阳理创编 2021.03.04

end

Z1(:,k)=Y1\\I';Zk1=Z1(:,k);%正序网络中节点m的自阻抗和互阻抗 Z2(:,k)=Y2\\I';Zk2=Z2(:,k);%负序网络中节点m的自阻抗和互阻抗 Z0(:,k)=Y0\\I';Zk0=Z0(:,k);%零序网络中节点m的自阻抗和互阻抗 %第二部分：计算支路电流 for p=1:N1 U1(p)=1; end

u1=U1'+uu1(:,k);%计算所有节点正序电压 u2=uu2(:,k);%计算所有节点负序电压 u0=uu0(:,k);%计算所有节点零序电压 for m=1:N1 for n=1:N1 if mmn(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 end end

%第三部分：计算发电机的端电压 for p=1:N2 if p==k

欧阳阳理创编 2021.03.04

欧阳阳理创编 2021.03.04

II(p)=-Ik1; else II(p)=0; end end

vv1(:,k)=YY1\\II.'; vv2(:,k)=YY2\\II.'; 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 3.3仿真 三相短路

欧阳阳理创编 2021.03.04

欧阳阳理创编 2021.03.04

图3.2三相短路运行图

欧阳阳理创编 2021.03.04

欧阳阳理创编 2021.03.04

不对称短路

欧阳阳理创编 2021.03.04

欧阳阳理创编 2021.03.04

欧阳阳理创编 2021.03.04

欧阳阳理创编 2021.03.04

图3.3不对称短路4 结论

欧阳阳理创编 2021.03.04

欧阳阳理创编 2021.03.04

经过MATLAB计算机算法的计算，得到的短路电流参数与手算相同，证明了相关MATLAB程序的正确性。通过两种算法的比较，计算机算法与传统手算相比较的优势不言自明。

MATLAB基础及其应用是一门实践性很强的专业课，MATLAB在当今社会发展异常迅速，已经从最初的“矩阵实验室”，渗透到科学与工程计算的多个领域，在自动控制、信号处理、神经网络、模糊逻辑、小波分析等多个方向，都有着广泛的应用，因此学好MATLAB对我们非常重要。

通过这次运用MATLAB计算短路电流，我从中学到了许多知识，学会了怎样把课本理论知识运动到实际中去。从确定课题后开始着手准备，我查阅了很多资料。在做设计时，也复习了很多专业课的知识，发现了以前知识上存在的漏洞，这使得我的专业知识得到了巩固和提高。

经过此次设计，我深刻的感觉到了MATLAB功能的强大，尤其是它的计算能力。由于水平有限，在这次设计中遇到了不少困难，也走了不少弯路，但最终还是完成了此次的课程设计，这当然是与老师的教导是分不开的，借此机会对张老师表示深深的谢意！

欧阳阳理创编 2021.03.04

欧阳阳理创编 2021.03.04

5参考文献

[1] 李光琦.电力系统暂态分析.中国电力出版社，2007

[2] 周全仁.张清益.电网计算与程序设计.湖南科学技术出版社，1983 [3] 陈珩.电力系统稳态分析.中国电力出版社，2007

[4] 管爱红.MATLAB基础及其应用教程.电子工业出版社，2009 [5] 尚涛.MATLAB基础及其应用教程.电子工业出版社，2014

时间：2021.03.05 创作：欧阳理 欧阳阳理创编 2021.03.04