电力系统三相短路的编程计算

1、 电力系统分析课程专题报告 学 生 姓 名： 班 级： 学 号： 指 导 教 师： 所 在 单 位： 提 交 日 期： 评分 2电力系统三相短路的编程计算某某某（学院， ）摘要：在电力系统中，三相短路故障造成的危害是最大的，发生的几率也最高，故短路计算对电力系统的稳定运行具有十分重要的意义。作为电力系统三大计算之一，计算三相短路故障发生时的短路电流、各节点电压、各支路电流是短路计算的基本内容。在电力系统短路电流的工程计算中，由于快速继电保护的应用，最重要的是计算短路电流基频交流分量的初始值，即次暂态电流。在给定电源电势时，实际上就是稳态交流电路的求解。本文基于教材例3-2应用MATLAB编程计

2、算三相短路故障的电流电压情况，并于例3-3进行进一步的验证和完善。关键词：三相短路；MATLAB中图分类号：TM 713  文献标识码：A0 引言在电力系统的四种短路类型中，三相短路是其中最严重的，其短路电流可达数万安以至十几万安，随之产生的热效应和电动力效应将使电气设备遭受严重破坏。因此，计算短路电流的主要应用目的是电力系统设计中的电气设备选择，短路计算已成为电力系统运行分析、设计计算的重要环节。实际电力系统短路电流交流分量初始值的计算，小型系统可以手算，而对于结构复杂的大型系统，短路电流计算量较大，用计算机进行辅助计算成为大势所趋。1 解析法求三相短路电流1.1 参数说明(1)

3、为了元件参数标幺值计算方便，取基准容量为，可设任意值，但必须唯一值参与计算。(2) 取基准电压为平均额定电压,基于例3-2中系统的额定电压等级有、，平均额定电压分别为、，平均额定电压与线路额定电压相差5%，为简化计算，故取平均额定电压。(3) 为次暂态短路电流有效值，短路电流周期分量的初值等于时间时的有效值。满足产生最大短路电流的三个条件下的最大次暂态短路电流作为计算依据。(4) 为冲击电流，即为短路电流的最大瞬时值（满足产生最大短路电流的三个条件及）。一般取冲击电流。(5) 为短路电流冲击系数，主要取决于电路衰减时间常数和短路发生的时刻，。1.2 电抗(1) 发电机电抗 式（1-1）式中为发

4、电机电抗百分数，由发电机铭牌参数。(2) 负载电抗 式（1-2）式中U为负载所在网络的电压，为负载容量，为负载无功功率。(3) 变压器电抗 式（1-3）式中为变压器阻抗电压百分数。因变压器中,即可忽略。(4) 线路电抗 式（1-4）式中为线路单位长度电抗。输电线路的等值电路中有四个参数，一般电抗，故。不做特殊说明，电导、电纳一般不计，故只求电抗标幺值。(5) 电动机电抗（近似值） 式（1-5）1.3 短路电流标幺值(1) 次暂态短路电流标幺值（取） 式（1-6）(2) 冲击电流 式（1-7）(3) 短路容量 式（1-8）1.4 分析过程 (1) 等效电路 教材例3-2环形电力系统各元件的参数如

5、下：发电机：，；变压器：，；线路：三条线路完全相同，长，电抗。基于此环形电力系统的等效电路如图2所示。 图1 例3-2系统图 图2 例3-2等效电路由式（1-1）得发电机电抗标幺值由式（1-3）得变压器电抗标幺值由式（1-4）得线路电抗标幺值(2) 短路前运行状况的分析计算如果要计及负荷，则必须进行一次潮流计算，以确定短路点开路电压以及各待求量的正常值。采用近似计算，即忽略负荷，所有电动势、电压均为1，电流均为零。，(3) 计算故障分量故障分量网络，即将电源接地，在短路母线对地之间加一个负电压（-1），如图所示。在此电路中即可求得母线的短路电流（略去右上角的两撇）、发电机G1、G2和线路-、-

6、的故障电流、和、以及母线、电压的故障分量、。图3（b）（f）为化简网络的步骤。一般讲，网络化简总是从离短路点最远处开始逐步消去除短路点外的其他节点。图3 例3-2网络化简过程由图3（f）可得为了求得网络中各点电压和电流的分布，由短路点向网络中其他部分倒退回去计算，从图3（f）到图3（d）可求得已知各母线电压即可求得任意线中的电流为求得为小得多，实际上是故障分量中母线和之间的平衡电流。如果要计算短路之后的,不能假定正常时的为零，因为此时可能是同一数量级。(4) 计算各待求量的有名值。发电机送至短路点的电流为实际流过发电机的电流为2 节点导纳矩阵求解2.1导纳矩阵等值电路例3-2环形网络生成的节点

7、导纳矩阵等值电路如图4所示，节点数为3。图4 节点导纳矩阵等值电路2.2导纳计算公式用节点导纳矩阵表示的网络节点方程为 式（2-1）式（2-1）中节点导纳矩阵是节点阻抗矩阵的逆矩阵，其对角元素自导纳为式（2-2）即除i节点外其他节点均接“地”时，自节点i看进网络的等值导纳，显然等于与i节点相连的支路导纳之和。非对角元素互导纳为式（2-3）即除i节点外其他节点均接“地”时，j点注入电流与i点电压之比，显然等于i-j支路导纳的负值。 由于网络中任一节点一般只和相邻的节点有连接支路，所以有很多为零，即节点导纳矩阵是十分稀疏的。由上可见，节点导纳矩阵极易形成，网络结构变化时也易于修改。2.3导纳矩阵的

8、形成由此得到节点导纳矩阵3 计算机编程3.1MATLAB编程语言以前的潮流计算采用传统的FORTRAN过程性语言，具有不灵活、不易理解、难于扩展等缺点，不利于发展。使用C语言、BASIC等这些开发工具开发电力系统分析程序，要求开发者不但要有足够的对于电力系统分析的知识，还要求开发人员必须精通编程语言，才能够编制出合格的程序，这样就提高了电力系统分析程序的编制难度。同时由于忽视了软件工程的要求，使得程序虽然对于用户很友好，但却使后继的程序开发人员难于继续工作。随着科学技术的发展，电力系统变得越来越复杂，电气工程师掌握一种好的并能对电力系统进行计算和仿真的软件是学习和研究的必然要求。MATLAB是

9、一种交互式、面向对象的程序设计语言，广泛应用于工业界与学术界，主要用于矩阵运算，同时在数值分析、自动控制模拟、数字信号处理、动态分析、绘图等方面也具有强大的功能。MATLAB程序设计语言结构完整，且具有优良的移植性，它的基本数据元素是不需要定义的数组，可以高效率地解决工业计算问题，特别是关于矩阵和矢量的计算。MATLAB与C语言和FORTRAN语言相比更容易被掌握。通过MATLAB语言，可以用类似数学公式的方式来编写算法，大大降低了程序所需的难度并节省了时间，从而可把主要的精力集中在算法的构思而不是编程上，大大提高了电力系统编程计算和仿真模型的效率。3.2程序流程图5 用节点导纳矩阵计算短路电

10、流的程序流程图3.2程序脚本见附录。3.3程序结果及分析（1） 例3-2结果 由前述经MATLAB编程可得结果如图6所示。（2） 举例验证结果相关参数数据如下。发电机：，。变压器：，。线路：10km，。图7 系统图 图8 等效电路这里取，仍如前述取各电压等级平均额定电压，求得的各元件电抗标幺值为发电机变压器线路在例3-2程序基础上，将参数替换并稍作修改即可得系统中短路电流、各节点电压、各支路电流，结果如图9所示。图6 例3-2编程计算结果 图9 另一例编程计算结果4 结论随着电力系统日益扩大和复杂化，计算机在研究电力系统中的应用日益广泛，MATLAB更是本科阶段电力系统编程计算和模型仿真中的基

11、础应用。基于教材例3-2，本文重在应用MATLAB对电力系统运行分析、设计计算的重要环节短路电流计算，根据其原理和方法进行编程计算，并在另一例题中得到验证和完善。在短路电流计算的基础上，扩充程序、建立设备数据库，构制对任意一个电力系统在任意点发生短路故障时的三相短路电流及其分布的通用程序，通过计算可以作为选择电气设备的依据，进行继电保护整定计算，将电力方面的专业课联结起来，这将是电力系统计算机应用的又一课题。参考文献1 方万良,李建华,王建学.电力系统暂态分析M.北京:中国电力出版社,2017.2 张志涌,杨祖樱.MATLAB教程M.北京:北京航空航天大学出版社.2015.3 秦媛媛,杨宛辉.

12、电力系统三相对称短路电流计算实用计算方法和程序J.维普,2015(12)：9-17.4 吴天明,赵新力.MATLAB电力系统设计与分析M.北京:国防工业出版社,2010.5 杨淑英.电力系统概论M.北京:中国电力出版社,2007.19附录1、例3-2源程序clear;clc;%三相短路计算%UB1=115;%基准值UB2=10.5;SB=60;IB1=SB/(sqrt(3)*UB1);XB1=UB12/SB;%电抗标幺值%变压器T1T3k=115/10.5;%变比Us=0.105;%短路电压标幺值UtN=115;%高压侧额定电压StN=120;%额定容量Xt=Us*SB/StN;%变压器电抗%

13、发电机SgN=100/0.86;%发电机额定容量UgN=10.5;%发电机额定电压xd=0.183;%同步电抗标幺值Xd=xd*SB/SgN;%折算后同步电抗值%线路xl=0.44;%单位长度电抗l=50;%线路长度Xl=xl*l/XB1;%线路阻抗%参数求取完毕%网络结构%1、首端节点；2、末端节点；3、支路电纳；4、首端节点对地支路电纳structure=1 2 1/Xl        2 3 1/Xl       1&#

14、160;3 1/Xl %1、节点；2、对地电纳      note=1 1/(Xt+Xd)   2 2/(Xt+Xd)   3 0;%由节点导纳阵生成阻抗阵%i=size(structure,1);Y=zeros(3,3);for x=1:1:3    for y=1:1:3        

15、;for m=1:1:i            if x=y  %对角线元素                if (structure(m,1)=x|structure(m,2)=y)       &

16、#160;            Y(x,y)=Y(x,y)+structure(m,3);                end            else  %非对角线元素

17、                if (structure(m,2)=y&structure(m,1)=x|structure(m,1)=y&structure(m,2)=x)                    

18、;Y(x,y)=-structure(m,3);                end            end        end    endendfor n=1:3 &

19、#160;  Y(n,n)=Y(n,n)+note(n,2);endZ=inv(Y);%阻抗阵%短路电流计算%节点3短路If=1/Z(3,3)*IB1;fprintf('3处短路电流为%f(kA)n',If);If13=Z(1,3)/Z(3,3)*If;fprintf('13短路电流为%f(kA)n',If13);If23=Z(2,3)/Z(3,3)*If;fprintf('23短路电流为%f(kA)n',If23);U1=If*Z(1,3)*XB1*sqrt(3);fprintf('1点电压为%f(kV)n'

20、,U1);U2=If*Z(2,3)*XB1*sqrt(3);fprintf('2点电压为%f(kV)n',U2);2、 另一例源程序clear;clc;%三相短路计算%UB1=230;%基准值UB2=10.5;SB=1000;IB1=SB/(sqrt(3)*UB1);XB1=UB12/SB;%电抗标幺值%变压器k=230/10.5;%变比Us=0.105;%短路电压标幺值UtN=230;%高压侧额定电压StN=370;%额定容量Xt=Us*SB/StN;%变压器电抗%发电机SgN=300/0.85;%发电机额定容量UgN=20;%发电机额定电压xd=0.168;%同步电抗标幺值Xd=xd*SB/SgN;%折算后同步电抗值%线路xl=0.321;%单位长度电抗l=10;%线路长度Xl=xl*l/XB1;%线路阻抗%参数求取完毕%网络结构%1、首端节点；2、末端节点；3、支路电纳；4、首端节点对地支路电纳structure=1 2 0 2 3 1/Xl 1 3 1/Xl ;%1、节点；2、对地电纳 note=1 1/(Xt