电力系统分析课程设计汇本_电力系统短路故障的计算机算法程序设计

电力系统分析课程设计电力系统短路故障的计算机算法程序设计佳琪学号2014409436班级20144094指导教师鲁明芳目录纳矩阵来计算短路电流直接利用导纳矩阵的元素来计算并不方便，实际的做法是利用的导纳矩阵来求出阻抗矩阵中的有关元素，然后仍然利用前述公式（

3.6）、（

3.7）、（

3.8）来进展计算一个网络的导纳矩阵与其阻抗矩阵存在互为逆矩阵的关系，所以一种方法是将导纳矩阵直接求逆，得出阻抗矩阵但当矩阵的阶数大时这样做计算量大，所以常用下面的方法当计算k点的短路电流时，从公式（

3.6）、（

3.7）、（

3.8）看到所需的阻抗矩阵元素是Z%Zk,…,Zkk,…,Z%而按定义当（二1时，4人“，也就是在k点注入单位电流，而其他节点注人电流都为零时，那么节点k的电压值就等于其自阻抗勿的其他各节点的电压值就等于各节点与节点k之间的互阻抗4七因此只要在计算机上进展下面的一次线性方程组的求解，就可以算出各节点的电压值

3.9解出的各点电压值就等于所需的阻抗值

3.10图

3.2给出了三相短路计算的原理框图例

3.1如图

3.3所示网络，母线

③发生三相直接短路，试作以下计算

（1）母线

③的故障电流

（2）故障后母线

1、2的电压各元件参数如下发电机Gl,100MV-A；G2,200MV•A额定电压均为io.5kV,次暂态电抗工1均为

0.2变压器Tl,100MV•A；T2,200MV•A变比均为

10.5kV/115kV,短路电压百分数均为10o电力线路三条电力线路（LI、L

2、L3）参数均为H5kV,60km,电抗/=

0.44Q/km,图

3.2三相短路计算的原理框图电容g=

0.008X10-6F/kmo负荷LD1,50MW,cos（p=

0.985；LD2,100MW,cos（p=l解元件参数的标幺值注于简化等值电路网络图

3.4中（元件参数的计算、网络的简化略），首先据图

16.4电路所示，形成节点导纳矩阵|LD24

③图

3.3例

3.1图图

3.4例

3.1简网络节点导纳矩阵为然后，对力求逆或解线性方程组，得节点阻抗矩阵为

③点短路电流为节点

①、

②的电压为例

3.2图

3.5所示电力系统，负荷全部略去，简化后的各电抗标幺值注于等值网络中试计算f点三相短路时的短路电流及网络中的电流分布

10.5kV TIi]5kV T

26.3kV CaL2|LD2LD1L3T3ld3GCH a系统图理

③域，j

0.2j

1.43b等值电路解题步骤1请输入短路点的数目NF2请输入节点数n3请输入支路数nl4输入各支路参数矩阵B矩阵B的每行是由以下参数构成的

①某支路的首端号P；

②末端号q，且Pq；

③支路的阻抗R+jX；

④支路的对地电纳；

⑤支路的变比K；

⑥折算到哪一侧的标志如果支路的首端P处于高压侧那么请输入“1〃，否那么请输入“0〃5输入由短路点号，短路点阻抗形成的矩阵D6请输入由各节点的初电压标幺值形成的列矩阵V07形成节点阻抗矩阵Z8求短路点电流⑼求网络中各节点的电压叵］10求网络中各支路电流解输入数据请输入短路点的数目NF=1请输入节点数n=4请输入支路数nl=5请输入由线路参数形成的矩阵B=[

010.2i010;0241010;

130.51i010;

230.59i010;

341.43i010]请输入由各节点的初电压标幺值形成的列矩阵vo=[l;l;l;l]请输入由短路号，短路点阻抗组成的矩阵口=[4,0]结果4点短路时的电流的标幺值If二0-

0.4890i各节点的电压标幺值u（节点号从小到大排）为

0.

91530.

73800.69930各支路短路电流的标幺值I（顺序同您输入B时一样）为:0+

4.5765i0+

0.1845i0-

0.4235i0-

0.0655i0-

0.4890i4附录程序清单

4.1形成节点导纳矩阵%本程序是形成节点导纳矩阵n=input请输入节点数n=，;nl=input请输入支路数nl-;B二input请输入由支路参数形成的矩阵B=;X=input请输入由节点号及其对地阻抗形成的矩阵X=;Y=zerosn;for i二l:nif Xi,2~=0;p=Xi,1;Yp,p=l./Xi,2;endendfor i=l:nlif Bi,6=0p=Bi,1;q=Bi,2;elsep=Bi,2;q=Bi,1;endY p,q=Yp,q-l./Bi,3*Bi,5;Y q,p=Y p,q;Y q,q=Yq,q+l./Bi,3*Bi,5厂2+Bi,

4./2;Y p,p=Yp,p+l./Bi,3+Bi,

4./2;enddispC导纳矩阵Y=;dispY

4.2形成节点阻抗矩阵为本程序是用支路追加法形成节点阻抗矩阵n=input请输入节点数n=，;nl=input C请输入支路数nl-;B=input请输入由支路参数形成的矩阵B=;m=0;Z=zerosn;for kl=l:nlp=Bkl,1;q=Bkl,2;if Bkl,6==01二L/Bkl,5;elsel=Bkl,5;endif p==0if qm%追加接地树枝Zq,q=Bkl,3;m=m+1;else%追加接地连支for k=l:m,Z k,m+l=-Z k,q;Z m+1,k=-Z q,k;endZ m+1,m+l=Z q,q+Bkl,3;for11=1:mfor k=l:mZ11,k=Z11,k-Z11,m+1*Z m+1,k,/Z m+1,m+1;endZll,m+l=0;endfor k=l:m+lZm+l,k=0endendelse if qm%追加不接地树枝for k=l:mZ k,q=Z k,p*1;Z q,k=Z p,k*1;endZq,q=l-2*Zp,p+l/2*Bkl,3;m=m+1;elsefor k=l:m%追加不接地连支Z k,m+1=1*Z k,p-Z k,q;Z m+1,k=1*Z p,k-z q,k;endZm+1,m+1=12*Zp,p+Zq,q-2*l*Zp,q+l2*Bkl,3;for11=1:mfor k=l:mZ11,k=Z11,k-Zll,m+1*Zm+1,k./Zm+1,m+1;endZll,m+l=0;endfor k=l:m+l Zm+1,k=0;endendendenddisp阻抗矩阵Z=;dispZ

11.

11.

22.

32.

32.

32.

32.

42.

42.

1.1关于电力系统短路故障的计算机算法程序设计目的电力系统正常运行的破坏多半是由于短路故障引起的，发生短路时，系统从一种状态剧变成另一种状态，并伴随复杂的暂态现象所谓短路故障，是指一切不正常的相与相之间或相与地发生通路的情况本文根据电力系统三相对称短路的特点，建立了合理的三相短路的数学模型，在此根底上，形成电力系统短路电流实用计算方法；节点阻抗矩阵的支路追加法编制了end endfork=l:NF/求各短路点的电流标幺值I D k,1,D k,1=V0D k,1,

1./Z D k,1,Dk,1+Dk,2;ft=num2strDk,1;tsl=点短路时’;ts2=C电流的标幺值If=;dn=strcatft,tsi,ts2;disp dn;dispIDk,l,Dk,1;for i=l:nVi,i=V0i,l-IDk,l,Dk,l*Zi,Dk,1;/求各节点的电压标end end对任意一个电力系统在任意点发生短路故障时三相短路电流及其分布的通用计算程序，该方法适用于各种复杂构造的电力系统从一个侧面展示了计算机应用于电力系统的广阔前景根据所给的电力系统，编制短路电流计算程序，通过计算机进展调试，最后完成一个切实可行的电力系统计算应用程序通过自己设计电力系统计算程序使同学对电力系统分析有进一步理解，同时加强计算机实际应用能力的训练电力系统的短路故障是严重的，而又是发生几率最多的故障，一般说来，最严重的短路是三相短路当发生短路时，其短路电流可达数万安以至十几万安，它们所产生的热效应和电动力效应将使电气设备遭受严重破环为此，当发生短路时，继电保护装置必须迅速切除故障线路，以防止故障局部继续遭受危害，并使非故障局部从不正常运行情况下解脱出来，这要求电气设备必须有足够的机械强度和热稳定度，开关电气设备必须具备足够的开断能力，即必须经得起可能最大短路的侵扰而不致损坏因此，电力系统短路电流计算是电力系统运行分析，设计计算的重要环节，许多电业设计单位和个人倾注极大精力从事这一工作的研究由于电力系统构造复杂，随着生产开展，技术进步系统日趋扩大和复杂化，短路电流计算工作量也随之增大，采用计算机辅助计算势在并行

1.2设计原理

1.

2.1计算机计算原理应用计算机进展电力系统计算，首先要掌握电力系统相应计算的数学模型；其次是运用合理的计算方法；第三那么是选择适宜的计算机语言编制计算程序建立电力系统计算的相关数学模型，就是建立用于描述电力系统相应计算的有关参数间的相互关系的数学方程式该数学模型的建立往往要突出问题的主要方，即考虑影响问题的主要因素，而忽略一些次要因素，使数学模型既能正确地反映实际问题，又使计算不过于复杂运用合理的计算方法，就是要求所选用的计算方法能快速准确地得出正确结果，同时还应要求在解算过程中占用存少，以利于提高计算机的解题规模选择适宜的语言编写程序，就是首先确定用什么计算机语言来编制程序；其次是做出计算的流程图；第三根据流程图用选择的语言编写计算程序然后上机调试，直到语法上无错误所编制的程序难免存在逻辑错误因此先用一个结果的系统作为例题进展计算用程序计算的结果和结果相比拟，如果结果相差甚远就要逐步分析程序的计算步骤，查出问题的出处；如果结果比拟接近，那么逐步分析误差来源；直到结果正确为止最后将所编制出的正确计算程序，用于电力系统的实际计算

1.

2.2电力系统短路计算计算机算法一般在电力系统短路计算中，多数情况下只要计算短路电流、电压的周期分量起始值因此，电力系统短路电流计算的数学模型，可归结为求解短路故障初始状态下的等值电路稳态解的问题对于三相对称短路，可建立一相等值电路的计算模型，对于不对称短路，那么可应用对称分量法建立系统的正、负、零序网络，从而建立故障计算的序网络模型计算机编程环境及编程语言的选择2本课程设计采用数学建模软件MATLAB2015a为主要设计工具MATLAB可以进展矩阵运算、绘制函数和数据、实现算法、创立用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域MATLAB的根本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C,FORTRAN等语言完成一样的事情简捷得多，并且MATLAB也吸收了像Maple等软件的优点，使MATLAB成为一个强大的数学软件在新的版本中也参加了对C,FORTRAN,C++,JAVA的支持

2.1优势特点1）高效的数值计算及符号计算功能，能使用户从繁杂的数学运算分析中解脱出来;2）具有完备的图形处理功能，实现计算结果和编程的可视化；3）友好的用户界面及接近数学表达式的自然化语言，使学者易于学习和掌握；4）功能丰富的应用工具箱（如信号处理工具箱、通信工具箱等），为用户提供了大量方便实用的处理工具

2.

1.1编程环境MATLAB由一系列工具组成这些工具方便用户使用MATLAB的函数和文件，其中许多工具采用的是图形用户界面包括MATLAB桌面和命令窗口、历史命令窗口、编辑器和调试器、路径搜索和用于用户浏览帮助、工作空间、文件的浏览器随着MATLAB的商业化以及软件本身的不断升级，MATLAB的用户界面也越来越精致，更加接近Windows的标准界面，人机交互性更强，操作更简单而且新版本的MATLAB提供了完整的联机查询、帮助系统，极大的方便了用户的使用简单的编程环境提供了比拟完备的调试系统，程序不必经过编译就可以直接运行，而且能够及时地报告出现的错误及进展出错原因分析

2.

1.2简单易用MATLAB是一个高级的矩阵/阵列语言，它包含控制语句、函数、数据构造、输入和输出和面向对象编程特点用户可以在命令窗口中将输入语句与执行命令同步，也可以先编写好一个较大的复杂的应用程序文件）后再一起运行新版本的MATLAB语言是基于最为流行的C++语言根底上的，因此语法特征与C++语言极为相似，而且更加简单,更加符合科技人员对数学表达式的书写格式使之更利于非计算机专业的科技人员使用而且这种语言可移植性好、可拓展性极强，这也是MATLAB能够深入到科学研究及工程计算各个领域的重要原因

2.

1.3强处理能力MATLAB是一个包含大量计算算法的集合其拥有600多个工程中要用到的数学运算函数，可以方便的实现用户所需的各种计算功能函数中所使用的算法都是科研和工程计算中的最新研究成果，而且经过了各种优化和容错处理在通常情况下，可以用它来代替底层编程语言，如C和C++在计算要求一样的情况下，使用MATLAB的编程工作量会大大减少MATLAB的这些函数集包括从最简单最根本的函数到诸如矩阵，特征向量、快速傅立叶变换的复杂函数函数所能解决的问题其大致包括矩阵运算和线性方程组的求解、微分方程及偏微分方程的组的求解、符号运算、傅立叶变换和数据的统计分析、工程中的优化问题、稀疏矩阵运算、复数的各种运算、三角函数和其他初等数学运算、多维数组操作以及建模动态仿真等

2.

1.4图形处理MATLAB自产生之日起就具有方便的数据可视化功能，以将向量和矩阵用图形表现出来，并且可以对图形进展标注和打印高层次的作图包括二维和三维的可视化、图象处理、动画和表达式作图可用于科学计算和工程绘图新版本的MATLAB对整个图形处理功能作了很大的改良和完善，使它不仅在一般数据可视化软件都具有的功能（例如二维曲线和三维曲面的绘制和处理等）方面更加完善，而且对于一些其他软件所没有的功能（例如图形的光照处理、色度处理以及四维数据的表现等），MATLAB同样表现了出色的处理能力同时对一些特殊的可视化要求，例如图形对话等，MATLAB也有相应的功能函数，保证了用户不同层次的要求另外新版本的MATLAB还着重在图形用户界面（GUI）的制作上作了很大的改善，对这方面有特殊要求的用户也可以得到满足

2.

1.5模块集和工具箱MATLAB对许多专门的领域都开发了功能强大的模块集和工具箱一般来说，它们都是由特定领域的专家开发的，用户可以直接使用工具箱学习、应用和评估不同的方法而不需要自己编写代码领域，诸如数据采集、数据库接口、概率统计、样条拟合、优化算法、偏微分方程求解、神经网络、小波分析、信号处理、图像处理、系统辨识、控制系统设计、LMI控制、鲁棒控制、模型预测、模糊逻辑、金融分析、地图工具、非线性控制设计、实时快速原型及半物理仿真、嵌入式系统开发、定点仿真、DSP与通讯、电力系统仿真等，都在工具箱〔Toolbox家族中有了自己的一席之地

2.

1.6程序接口新版本的MATLAB可以利用MATLAB编译器和C/C++数学库和图形库,将自己的MATLAB程序自动转换为独立于MATLAB运行的C和C++代码允许用户编写可以和MATLAB进展交互的C或C++语言程序另外，MATLAB网页效劳程序还容许在Web应用中使用自己的MATLAB数学和图形程序MATLAB的一个重要特色就是具有一套程序扩展系统和一组称之为工具箱的特殊应用子程序工具箱是MATLAB函数的子程序库，每一个工具箱都是为某一类学科专业和应用而定制的，主要包括信号处理、控制系统、神经网络、模糊逻辑、小波分析和系统仿真等方面的应用

2.

1.7应用软件开发在开发环境中，使用户更方便地控制多个文件和图形窗口；在编程方面支持了函数嵌套，有条件中断等；在图形化方面，有了更强大的图形标注和处理功能，包括对性对起连接注释等；在输入输出方面，可以直接向Excel和HDF5进展连接对称故障的计算机算法3对称故障是指三相短路故障这里首先通过一个简单的例子进展说明L13T1右图⑸所示的电力系统有三条母线，标为

①②③，其等值电路如图b所示，图中

41、「

2、L12代表和该母线连接的输电线的容抗总和，代表母ZL L23系统图a线

③的负荷阻抗，Zf代表故障阻抗，故障发生在母线

③图b也可简化为图c所示，是一个有三个节点的有源网络，□、

4、必是节点电压，

4、‘

2、L是节点注入电流，故障阻抗Zf是外加到节点

③的三相对称阻抗，当Zf=O时，代表三相直接短路右图c也可以看成是有三个端口的有源网络，每个节点与零电位点构成一个端口，按叠加定理可以通过阻抗型参数方程或导纳型参数方程来表示其电压和电流的关系以下分别进展叙述图

3.1电力系统接线网络及简化网络

3.1用阻抗矩阵计算三相短路电流图

16.2c网络的阻抗型参数方程如下:式中，9（）、%）、氏⑹为各节点的开路电压（也即注入电流/尸人工二时各节点的电压）,开路电压可由正常运行的潮流计算求得近似计算中，那么设各节点开路电压标幺值为l.Oo阻抗矩阵中的各元素Zn、Z.等，为各节点的自阻抗和各节点之间的互阻抗,据12自、互阻抗的物理意义可以确定其数值如下其余类推当在节点

③发生三相故障时，相当于在节点

③接上故障阻抗Zf,流过Zf的故障电流%其参考方向如图

16.2（c）所示，其他节点没有外接电路，所以其注入电流为零，因此节点

③故障时的边界条件是

3.2将公式⑶1）与⑶2）联立求解，得

3.3公式⑶3）就是计算故障电流的数学模型当Zf给定时，只要知道故障点的开路电压4（）和自阻抗Z33就可以算出人求出乙后代入公式⑶1）可求得各节点的电压为将上述关系推广到有n个节点的电力系统，那么其阻抗型参数方程为

3.4设在节点k发生三相故障，故障阻抗为Zf时，其边界条件

3.5联立求解式

3.4和

3.5,得

3.6故障电流求出后，代入式⑶4求出各节点电压

3.7各节点电压求得后，可按下式求各支路电流

3.8式中，Zu为连接节点i与节点j的支路阻抗在略去输电线电容电流的条件下，支路电流也就是输电线电流公式

3.

6、

3.

7、⑶8就是计算三相故障的根本数学模型，从式中看到当句给定后，只需知道节点的开路电压4⑹和阻抗矩阵中的元素Z,就可以求出需要的结果节点的开路电压可以由正常的潮流计算得出，阻抗矩阵中的所有元素可以用支路追加法求得当这些量都已求出并储存于计算机中，计算短路电流的工作就很简单要计算任一节点的短路电流和电压、电流分布时，只要按上述公式编好程序，取出有关的开路电压，有关的自阻抗、互阻抗进展计算便可通常利用支路追加法可直接形成节点阻抗矩阵

3.2用节点导纳矩阵计算三相短路电流导纳矩阵易于形成，而且是稀疏矩阵，所以占用计算机的存容量少，因此常利用导。