电力系统分析

《电力系统分析》由会员分享，可在线阅读，更多相关《电力系统分析（23页珍藏版）》请在金锄头文库上搜索。

1、2011-2012 学年度下学期电力系统分析课程设计电力系统短路故障的计算机算法程序设计姓名闫漫学号030940815班级0309408指导教师 钟建伟目录1 任务提出与方案论证11.1目的11.2课程设计说明11.3 MATLAB 语言介绍12总体设计32.1短路电流的计算原理32.1.1 节点导纳矩阵计算电压初始值的原理32.1.2节点阻抗矩阵计算短路电流的原理32.2 计算条件42.2.1.标么值的概念42.2.2 等值电路42.2.3 戴维南定理和诺顿定理的转换43详细设计53.1 程序流程图53.2 等值电路及参数53.2.1例 6-3 的等值电路53.2.2习题 6-10 的等值网

2、络及参数63.2.3习题 6-11 等值网络及参数73.3程序及说明 73.3.1.通用程序为：73.3.2例 6-3 的输入输出数据及具体程序83.3.3习题 6-10 的输入输出数据及具体程序103.3.4习题 6-11 输入输出数据及具体程序：134 总结15参考文献16信息工程学院课程设计任务书学生姓名闫漫学 号030940815成 绩设计题目电力系统短路故障的计算机算法程序设计设计 内 容用节点阻抗矩阵法计算短路电流，首先是建立等值网络，求出其节点导纳矩 阵，通过其逆的方法得到节点阻抗矩阵；其次将等值网路用诺顿定理进行变换, 得到电流源矩阵I,执行语句V=ZI就得到了各节点的初始电压

3、，然后利用电力 系统分析弟八早中的公式（6-10）和（6-11）分别计算短路电流和节点电压故障分 量,短路后的节点电压就等于短路前的初始电压和短路后节点电压故障分量的叠 加。求得了各节点电压即可求得各支路电流。设计要求给定任意等值网络能算出各节点的电压、短路节点的短路电流和各支路的短路电 流时 间 安 排3月5日-3月8日3月12日-3月14日参 考 资 料电力系统分析（上册）华中科技大学出版社何仰赞 MATLAB R2006基础篇 化学工业出版社 曹岩李明雨1 任务提出与方案论证1.1 目的根据所给的电力系统，编制短路电流计算程序，通过计算机进行调试，最后完成一 个切实可行的电力系统计算应用

4、程序。通过自己设计电力系统计算程序使同学对电力系 统分析有进一步理解，同时加强计算机实际应用能力的训练。1.2 课程设计说明本文根据电力系统三相对称短路的特点，建立了合理的三相短路的数学模型， 在此基础上， 形成电力系统短路电流实用计算方法。所谓短路，是指电力系统正常运行情况以外的相与相之间或相与地之间的“短 接”。在电力系统正常运行时，除中性点外，相与相或相与地之间是绝缘的。如果由 于某种原因使其绝缘破坏而构成了通路，我们就称电力系统是发生了短路故障。在 三相系统中，短路故障可分为两大类：即对称短路（三相短路）和不对称短路（两 相短路、两相接地短路、单相接地短路）。 电力系统短路电流计算是电

5、力系统运行分析，设计计算的重要环节，许多电业设计 单位和个人倾注极大精力从事这一工作。其中三相短路虽然发生的机会较少，但情 况严重，又是研究其它短路的基础。所以我们先研究最简单的三相对称短路时短路 电流的计算。1.3 MATLAB 语言介绍我选择的计算机语言是MATLAB。MATLAB是一门用途广泛的语言，它提供了专业水平的符号计算，文字处理，可视化建模仿真和实时控制等功能。MATLAB 的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十 分相似,故用 MATLAB 来解决问题要比用 C,FORTRAN 等语言方便的多。在新的版本中 也加入了对C,FORTRAN,c+ ,JAVA的

6、支持，可以直接调用，用户也可以将自己编写的 实用程序导入到MATLAB函数库中方便自己以后调用。MATLAB 开发环境是一套方便用户使用的 MATLAB 函数和文件工具集， 其中许多工具是图形化用户接口。它是一个集成的 用户工作空间，允许用户输 入输出数据，并提供了 M 文件的集成编译和调试环境，包括 MATLAB 桌面、 命令窗口、 M 文件编辑调试器、 MATLAB 工作空间和在线帮助文档。MATLAB 的程序极其简短且可移植性很好，基本上不做修改就可以在各种型号的 计算机和操作系统上运行。MATLAB既具有结构化的控制语句（如for循环，while 循环， break 语句和 if 语句

7、），又有面向对象编程的特性。但较之高级语言，由于 其程序不用编译等预处理，也不生成可执行文件，而是边解释边执行，所以速度较 慢。总之，MATLAB是一个高度集成的系统，集科学计算、图象处理、声音处理于一 体，具有极高的编程效率。近年来，MATLAB已经从最初的“矩阵实验室”，渗透到 科学与工程计算的多个领域，在自动控制、信号处理、神经网络、模糊逻辑、小波 分析等多个方向，都有着广泛的应用。2 总体设计2.1 短路电流的计算原理2.1.1 节点导纳矩阵计算电压初始值的原理在电力系统分析的第四章中学到了矩阵方程:YV=I,故由矩阵性质得：V=ZI2.1.2 节点阻抗矩阵计算短路电流的原理利用节点阻

8、抗矩阵计算短路电流如图所示，假定系统中的节点f经过过渡阻抗 发生短路。对于正常状态的网络而言，发生短路相当于在故障节点f增加了一个注 入电流。因此，网络中任一节点i的电压可表示为： f公式一：V=Y ZI Z I 公式（6-3）iij jif fj G由式可见，任一节点i的电压都由两项叠加而成。第一项是当时由网络内所有电源 在节点 i 产生的电压，也就是短路前瞬间正常运行状态下的节点电压，记为 V（0）。i 第二项是当网络中所有电流源都断开，电势源都短接时，仅仅由短路电流 I f 在节 点i产生的电压。这两个分量的叠加，就等于发生短路后节点i的实际电压，即 公式二：V i - V .（ ）-

9、z i f 公式（6-4）iiif f公式二也适用于故障节点f,于是有：V二（0）ZIZ ff I f 公式（ 6-5）Z是故障节点f的自阻抗，也称输入阻抗。ff方程式含有两个未知量 Vf,I f ，根据故障的边界条件：公式（6-6）V f f z I = ._VZ (0)由以上两个方程式解出： J = 公式(6-7)f Z + zff f2.2计算条件2.2.1. 标么值的概念计算时选定一个基准容量(Sjz)和基准电压(Ujz).将短路计算中各个参数都转化为 和该参数的基准量的比值(相对于基准量的比值)，称为标么值(这是短路电流计算最特别 的地方，目的是要简化计算).(1) 基准基准容量Sj

10、z =100 MVA基准电压 UJZ 规定为 8 级.230, 115, 37, 10.5, 6.3, 3.15 ,0.4, 0.23 KV 有了以上两项，各级电压的基准电流即可计算出，例：UJZ (KV)3710.56.30.4 因为 S=1.73*U*I 所以 IJZ (KA) 1.565.59.16144(2) 标么值计算容量标么值S* =S/SJZ例如：当10KV母线上短路容量为200 MVA时,其标么值容量 S* = 200/100=2.电压标么值U*= U/UJZ ;电流标么值I* =I/IJZ 有了标幺值的概念后，短路电流计算的过程就大大减少了，简化了运算。2.2.2 等值电路计

11、算前先要做出等值电路，并计算等值参数。在详细设计中会将例 6-3，习题 6-10及习题6-11的等值电路及参数列写出来。2.2.3 戴维南定理和诺顿定理的转换欲求初始节点电压，要将戴维南形式的等值电路化为诺顿形式的等值电路，以 便得到电流矩阵，从而利用公式V=ZI求得节点初始电压。诺顿定理得到一个理想电流源I与一个电导G并联联(习惯上与电流源对应使用的 都用电导，不用电阻，当然这只是叫法不同而已)戴维南定理得到一个理想电压源U与一个电阻R串联两者的转换关系为GR=1 (即诺顿定理得到的电阻与戴维南定理得到的电阻的阻 值相等)，且U=IR(即诺顿定理得到的电流源的电流、电阻与戴维南定理得到的电

12、压源电压形式上满足欧姆定律)3.1 程序流程图3详细设计3.2 等值电路及参数3.2.1 例6-3的等值电路输出结果z12=j0.105 z45=j0.184 z24=j0.08 z23=j0.065 z34=j0.05Y11=-j16.1905 Y12=j9.05238Y21=j905238 Y22=-j37.4084 Y23=j15.3846Y24=j12.500 Y32=J15.3846Y33=-j35.3846Y34=j20.000,Y42=j15.5000Y43=-j37.9348Y45=j5.4348Y54=j5.4348Y55=-j9.98023.2.2 习题 6-10 的等值网

13、络及参数选SB=120MV.A,VB=Vav,各元件标幺值参数如下：X10=0.105 X20=0.28 X24=0.21 X43=0.43554 X35=0.29036 X45=0.25406Y340=Y430=j0.01852Y350=Y530=j0.01234Y450=Y540=j0.01080故节点导纳为：Y11=-j13.8716Y13=Y31=j9.5238Y22=-j8.33343 Y24=Y42=j4.7619Y33=-j15.2329Y34=Y43=j2.2960Y55=-j7.3569Y35=Y53=j3.4440Y44=-j10.9646Y45=Y54=j3.93603.

14、2.3 习题 6-11 等值网络及参数节点导纳为：Y11=-j10 Y12=Y21=j5function dl63text63;Z=inv(Y)I=IV=Z*Ii=input(请输入节点i=，)Ii=1/Z(i,i)Y22=-j14.5 Y23=Y32=j10 Y33=-j17.66673.3 程序及说明3.3.1. 通用程序为：%定义一个函数名为dl63的函数%调用函数名为t ex t63的函数%求节点阻抗矩阵%求电流的转置矩阵% %求各节点的初始电压%输入短路节点i %求短路点电流%for k=1:mkVk=V(k,1)-Z(k,i)*Iiendfor p=1:mfor q=(1+p):mif(Y(p,q)=0)p,qVp=V(p,1)-Z(p,i)*Ii;Vq=V(q,1)-Z(q,i)*Ii;Ipq=-(Vp-Vq)*Y(p,q)%求各节点的电压%如果p, q节点间的导纳不为0 %输出P,q %分别求P,q点电压%求P,q两节点间的支路电流endendend3.3.2例6-3的输入输出数据及具体程序输入数据为：Y=-16.1905j9.5238j0009.5238j-37.4084j15.3846j12.5000j