复杂网络NR法潮流分析与计算的设计

word文档可自由复制编辑word文档可自由复制编辑word文档可自由复制编辑课程设计（论文）课程名称电力系统课程设计题目名称复杂网络N-R法潮流分析与计算的设计学生学部（系）电气工程系专业班级11电气工程及其自动化（）班学号学生姓名指导教师2014年6月12日课程设计（论文）任务书课程设计（论文）任务书一、课程设计（论文）的内容1.掌握潮流计算的基本原理2.根据电力系统网络推导电力网络数学模型，写出节点导纳矩阵；3.赋予各节点电压变量初值后，求解不平衡量；4.形成雅克比矩阵；5.求解修正量后，重新修改初值，从2开始重新循环计算；6.求解的电压变量达到要求的精度时，再计算各支路的功率分布、功率损耗和平衡节点功率；7.上机编程调试；8.计算分析给定系统潮流，并与手工计算结果进行比较分析；9.书写课程设计说明书（电子版），并打印纸质版和电子版同时上交。二、课程设计（论文）的要求与数据1.运用某种语言（C语言或Matlab语言等），给出程序源代码和程序说明；2.给定题目的输入输出文件；3.给出程序的计算过程；4.以某一系统为例，给出程序的手算过程（至少迭代2次）。三、课程设计（论文）应完成的工作1.按照格式规范，独立完成课程设计说明书的撰写；2.完成潮流计算程序的编写和调试；3.完成潮流计算的手算过程。题目名称复杂网络N-R法潮流分析与计算的设计学生学部（系）机械电气学部电气工程系专业班级11电气工程及其自动化（？）班姓名学号四、课程设计（论文）进程安排序号设计（论文）各阶段内容地点起止日期1资料收集，潮流计算的算法分析图书馆2014.6.3-6.42软件实现及相关程序的调试实验室机房6.5-6.83课程设计说明书撰写实验室机房6.94课程设计上交1-1166.10五、应收集的资料及主要参考文献[1]于永源，杨绮雯.电力系统分析[M].北京：中国电力出版社，2007[2]孟祥萍.电力系统分析[M].北京：高等教育出版社，2005[3]何仰赞.电力系统分析[M].武汉：华中理工大学出版社，2002[4]周开利，邓春晖.Matlab基础及其应用教程[M].北京：北京大学出版社，2007[5]谭浩强.C语言程序设计[M].北京：清华大学出版社，2005等等发出任务书日期：2014年6月3日指导教师签名：计划完成日期：2014年6月10日教学单位责任人签章：摘要随着电力事业的快速发展，电力电子新技术得到了广泛应用；出于技术、经济等方面的考虑，500kV及以上的超高压输电线路普遍不换位，再加上大量非线性元件的应用，电力系统的不对称问题日益严重。因此电力系统不对称故障分析与计算显得尤为重要。基于对称分量法的基本理论，对称分量法采取的具体方法之一是解析法，即把该网络分解为正，负，零序三个对称序网，这三组对称序分量可分别按对称的三相电路分解。计算机程序法。通过计算机形成三个序网的节点导纳矩阵，然后利用高斯消去法通过相应公式对他们进行数据运算，即可求得故障端点的等值阻抗。最后根据故障类型选取相关公式计算故障处各序电流，电压，进而合成三相电流电压。进行了参数不对称电网故障计算方法的研究。通过引计算机算法，系统介绍电网参数不对称的计算机算法方法。根据断相故障和短路故障的特点，通过在故障点引入计算机算法，，给出了各种断相故障和短路故障的仿真计算。此方法以将故障电网分为对称网络和不网络两部分，在程序法则下建立起不对称电网故障计算统一模型，根据线性电路的基本理论，并借助于相序参数变换技术完成故障计算。关键词:参数不对称;电网;故障计算

目录 前言.........................................................11.电力系统短路故障的基本知识...................................21.1短路故障的概述...........................................21.2标幺制...................................................2对称分量法在不对称短路计算中的应用...........................2不对称三相量的分解.......................................3对称分量法在不对称短路计算中的应用.......................4简单不对称短路的分析与计算...................................4单相（a相）接地短路......................................7两相（b，c相）短路......................................7两相（b相和c相）短路接地...............................7简单不对称短路的分析与计算计算机计算程序法...................8简单故障的计算程序原理...................................9网络节点方程的形成......................................10电力系统不对称短路计算实例..................................11单相接地短路和两相短路不对称故障分析与计算..............11两种计算方法的对比......................................18结语.........................................................19参考资料.....................................................19附录：不对称短路电流计算程序.................................23前言电力系统的安全、稳定、经济运行无疑是历代电力工作者所致力追求的，但是从电力系统建立之初至今电力系统就一直伴随着故障的发生而且电力系统的故障类型多样。所以电力系统故障分析计算方法年来一直是学术研究的热点，其为预防及消除电力系统的故障准备，必要的理论知识。引起电网不对称的原因，超高压架空输电线路不换位.变压器结构不对称.交直流变换器的存在，系统负荷不平衡，系统中存在非线性元件。建立用于故障分析的网络的数学模型是计算机编程计算的基础．选取合适的数学模型直接关系到故障分析计算的优劣。当前一般采用节点导纳矩阵或节点阻抗矩阵作为原网数学模型．节点导纳矩阵描述的是网络的短路参数，只包含网络的局部信息，因而具有较好的稀疏性。因此，节点导纳矩阵主要运用于解题规模为主要矛盾的情况,节点阻抗矩阵要运用于进行大量故障计算的情况。随着计算机技术的快速发展，特别是计算机内存容量的扩大和计算速度的提高，节点阻抗矩阵的解题规模和解题速度都有了很大的提高，一般可以满足大多数电网故障分析计算的需要。电网中的故障可分为两大类：简单故障和复杂故障。复杂故障一般可以由两种或两种以上的简单故障组合而成，简单故障又分为对称型故障和不对称型故障；而不对称故障又可以分为短路故障(横向故障)和断路故障(纵向故障)。故障计算方法总体上可以划分为两种：一种解析法，另一种计算机程序法。1.电力系统短路故障的基本知识1.1短路故障的概述在电力系统运行过程中，时常发生故障，其中大多数是短路故障。所谓短路：是指电力系统正常运行情况以外的相与相之间或相与地（或中性线）之间的连接。除中性点外，相与相或相与地之间都是绝缘的。电力系统短路可分为三相短路，单相接地短路。两相短路和两相接地短路等。三相短路的三相回路依旧是对称的，故称为不对称短路。其他的几种短路的三相回路均不对称，故称为不对称短路。电力系统运行经念表明，单相短路占大多数，上述短路均是指在同一地点短路，实际上也可能在不同地点同时发生短路，例如两相在不同地点接地短路。依照短路发生的地点和持续时间不同，它的后果可能使用户的供电情况部分地或全部地发生故障。当在有由多发电厂组成的电力系统发生端来了时，其后果更为严重，由于短路造成电网电压的大幅度下降，可能导致并行运行的发电机失去同步，或者导致电网枢纽点电压崩溃，所有这些可能引起电力系统瓦解而造成大面积的停电事故，这是最危险的后果。产生短路的原因很多主要有如下几个方面（1）原件损坏，例如绝缘材料的自然老化，设计，安装及维护不良所带来的设备缺陷发展成短路，(2)气象条件恶化，例如雷电造成的闪络放电或避雷针动作，架空线路由于大风或导线覆冰引起电杆倒塌等。(3)违规操作，例如运行人员带负荷拉刀闸(4)其他，例如挖沟损伤电缆。在电力系统和电气设备的设计和运行中，短路计算是解决一系列技术问题所不可缺少的基本计算，比如在选择发电厂和电力系统的主接线时为啦比较不同方案接图，进行电力系统暂态稳定计算，研究短路对用户的影响。合理配置各种继电保护和自动装置并正确整定其参数都必须进行短路的计算和分析.1.2标幺制在短路计算中，各电气量如电流.电压.电阻.电抗.功率等数值可以用名值表示，也可以用标幺值表示，为了计算方便，通常在1KV以下的低压系统用有名值，而高压等级中有多个电压等级，存在电抗换算问题，宜采用标幺值。1.21标幺制概念所谓标幺制，就是把各个物理量均用标幺值来表示的一种相对单位制。某一物理量的标幺值A,等于它的实际值A与所选定的基准值A的比值，即：dAA（1-1）Ad在进行标幺值计算时，首先选定基准值。基准值原则上可以任意选定，但因物理量之间有内在的必然联系，所以并非所有的基准值都可以任意选取。在短路计算中经常用到的四个物理量容量S，电压U，电流I，和电抗X。通常选定基准容量S和基准电压Ud d则基准电流和基准电抗X分别为：dSId（1-2）d 3Ud U U2Xdd（1-3）d 3I Sd d为了计算方便，常取基准容量S=100MV.A，基准电压用各级线路的平均额定电压，即dUU.。d av所谓线路平均额定电压，是指线路始端最大额定电压与线路末端最小额定电压的平均值，如表1-1所示。表1-1线路的额定电压与平均额定电压额定电压U/kv0.220.3836103560110220330平U/kv0.230.43.156.310.53763115230345在，电抗器等所给出的标幺值，都是以其本身额定值为基准的标幺值或百分值。由于各电气设备的额定值往往不尽相同，基准值不同的标幺值是不能直接进行运算的，因此，必须把不同基准值的标幺值换算成统一基准值的标幺值。1.2.2电力系统各原件电抗标幺值的计算发电机电抗标幺值：X%SXGB（1-4） G 100S式中X%--------发电机电抗百分数，有发电名牌参数的X100%X G d G;S----已设定的基准容量（基准功率），MV.A；S-----发电机额定容量。B N负载电抗标幺值：U2XQ（1-5） LS LLU2----原件所在网络的电压标幺值；S----负载容量标幺值；Q----负载无功率标幺 L L值变压器电抗标幺值： U% SXKB（1-6） T 100S在变压器中电抗XR，即R忽要是指电抗，由变压器电抗有名值推出抗标值为： S U% XBU2K （1-7） TU2 NT100式中U%-----变压器阻抗百分；S-----基准容量，MV.A;SU------变压器铭牌K B NT,NT参数给定额定容量,MV.A;U-----基准电压U取平均电压U,KV。 B B av4线路电抗标幺值：SXBxL（1--8）U20B式中：x--线路单位长度电抗；L---线路单位长度，KM。S--基准容量.MV.AU-线路额定均电压。基准电压UU,kv输电线路的等值路有四个参数，一般抗R,故R0。由，故电导.电纳一般不计，故而只求电抗标幺值。W W W2.对称分量法在不对称短路计算中的应用2.1不对称三相量的分解人们在长期的实践中发现，在三相电路中，任意一组不对称的三相相量（电压或电流），可以分解为三组三相对称的相量分量，式（2-1）。在线性电路中，可以用叠加原理对这三组对称分量按照三相电路去解，然后将其结果叠加起来。就是不对称三相电路的解答，这个方法就叫做对称分量法。设F，F，F为三相系统中任意一组不对称的三相量，可以分解为三组对称的三 a b c序分量如下：FFFFFFFF（2-1）FFFFa(1) a(2) a(b(1) b(2) b(0)c(1) c(2) c(三相序分量如图2-1所示：)1(a)1(cI)1(a)1(cIIb(0)I Ic(0)b(1)II c(2) b(2) 正序分量负序分量零序分量图2-1三序分量正序分量：F.F.F三相的正序分量大小相等，彼此相位相差120，与系统正a(1)b(1)c(1)常对称运行对称运行方式下的相序相同，达到最大值abc，在电机内部产生正转磁场，这就是正序分量。此正序分量为一平衡的三相系统，因此有：FFF0. a(1) b(1) c(1)负序分量：F.F.F三相的负序分量大小相等，彼此相位相差120，与系统正a(2)b(2)c(2)常对称运行对称运行方式下的相序相反，达到最大值acb，在电机内部产生反转磁场，这就是正序分量。此正序分量为一平衡的三相系统，因此有：FFF0. a(2) b(2) c(2)零序分量：F.F.F三相的零序分量大小相等，相位相同，三相的零序分量同a(0)b(0)c(0)时达到最大值，在电机内部产生漏磁，其合成磁场为零，这就是零序分量。如果以A相为基准相，各序分量有如下关系：正序分量：FFa2FFa2FaF a(1) b(1) a(1) c(1) b(1) a(1)负序分量：FFaFFaFa2F零序分量：FFFFF a(2) b(2) a(2) c(2) b(2) a(2) a(0) b(0) a(0) c(0) a(0)331133其中：aej120j；a2ej240j2 2 21aa20aa21；a31于是有：FFFFa a(1) a(2) a(0)Fa2FaFF（2-2） b a(1) a(2) a(0)FaFa2FFc a(1) a(2) a(0)F111Faa2a1Faa(2)(1)（2-3）FbFcaa21Fa(0)则有：Ib(1)a2Ia(1)Ic(1)aIa(1)Ib(2)aIa(2)Ic(2)a2Ia(2)Ib(0)Ic(0)Ia(0)其逆关系式为：F1aa(1)1Fa(2)31a2F11a(0)a2FaaFb（2-4）1Fc这样根据式（2-3）可以把三组三相对称向量合成三个不对称向量，而根据式（2-4）可以把三个不对称向量分解成三组对称量。2.2对称分量法在不对称短路计算中的应用电力系统的正常运行一般是对称的，它的三相电路的参数相同，各相的电流，电压对称，这就是说只有正序分量存在。当电力系统的某一点发生不对称故障时，三相电路的对称条件受到破坏，三相对称电路就成为不对称的了。此时，可用对称分量法，将实际的故障系统变成三个互相独立的序分量系统，而每个序分量系统本身又是三相对称的，从而就可以用进行电路计算了。图2-2简单系统单相接地故障图如图2-2所示的简单系统发生单相接地短路故障。应用对称分量法，可绘出三序网图（三序等值电路图），如图2-3所示为最简化的三序网图，三序网的参数可分为正序，负序，零序。图中Z,Z ,Z 分为正序阻抗，负序阻抗，零序阻抗。 1 2 0图2-3简化三序网图列出电压方程：EIZUa1a11a1IZUa2（2-5）a22IZUa0a00由此可见，应用对称分量法进行不对称故障计算时，其关键问题是先求出各序网络的等效电抗（即要求出系统中各主要原件（发电机，变压器，线路等）的各序电抗值），然后根据短路的类型，边界条件，把正，负，零序网连接成串，并联的形式，从而可求出电流，电压的各序分量，再应用对称分量法进而可求出各相电流和电压等.3简单不对称短路的分析与计算（解析法）电力系统发生不对称故障时，无论是单相接地短路，还是两相短路，两相短路接地，只是短路点的电压，电流出现不对称，利用对称分量法将不对称的电流电压分解为三组对称的序分量，由于每一序系统中三相对称，则在选好一相为基准后，每一序只需要计算一相即可，用对称分量法计算电力系统的不对称故障。其大概步骤如下：计算电力系统各个原件的序阻抗；制定电力系统的各序网络；由各序网络和故障列出对应方程；从联立方程组解出故障点电流和电压的各序分量，将相对应的各序分量相加，以求得故障点的各相电流和电压；计算各序电流和各序电压在网络中的分布，进而求得各指定支路的各相电流和指定节点的各相电压。【4】3.1单相（a相）接地短路单相接地短路时的系统接线图如图3-1所示，故障处的三相的边界条件：U0Ib0Ic0a图3-1图3-1单相接地短路时的等值接线图UUU0 a1 a2 a0a2IaII0 a1 a2 a0aIa2II0（3-1） a1 a2 a0即有：1IIII（3-2）a1a2a03a根据单相接地短路时的边界条件（式3-1，式3-2）连接复合网，如图3-2所示。由复合网图可以写出各序分量：图3-2图3-2单相接地短路时的复合序网120120aaaEIIIZZZ111aaUEIZ222aaUIZ000aaUIZ（3-3）于是可以用对称分量法得到短路点的各相电流电压：IIII3Ia1 a2 a0 a1II0cUUUU0 a a1 a2 a0Uba2Ua1aUa2Ua0Ia1a2aZ2a21Z0UaUa2UUIaa2Z2a1Z0（3-4） c a1 a2 a0 a1短路点电流，电压的相量图如图3-3所示。这里按纯电感性电路画的，电流滞后电压90，若不是纯电感电路，则电流与电压角度由Z Z 的阻抗角决定，一般小于90。在 2 0相量图中，将每相的序分量相加，得各相电流电压的大小和相位。图3-3单相接地时短路处的电流，电压的相量图3.2两相（b,c相）短路两相短路的系统接线如图3-4所示，在k点发生b，c两相短路。短路点的边界条件：图3-4两相短路等值接线图I0IIU0UU0 a b c bc b c用对称分量表示为：III0a2II(aIa2II)（3-5） a1 a2 a0 a1 a2 a0a2UaUUaUa2UU a1 a2 a0 a1 a2 a0于是有：III0UU(3-6)图3-5两相短路的复合序网图从复合序网中可以直接求出电流，电压的各序分量：1221a图3-5两相短路的复合序网图从复合序网中可以直接求出电流，电压的各序分量：1221aaEIIZZ121122aaaaUUEIZIZ（3-7）III0a1 a2IIa2IaI(a2a)Ij3Ic a1 a2 a1a1UUUU2U2IZ a a1 a2 a0 a1 a121UUaUUa2UUU（3-8） b c a2 a0 a1 a1 2a短路点电压和电流的相量图。如图3-6所示，这里任然是纯电感电路。电流滞后电压90图3-6两相接地短路点电流和电压相量图3.3两相（b相和c相）短路接地两相短路接地时系统接线图如图3-7所示，短路点的边界条件为：图3-7两相短路接地的等值接线图I0UU0 a b c用对称分量表示为：III0 a1 a2 a01UUUU（3-9） a1 a2 a03a由（3-9）可以画出两相短路接地的复合序网图是三个序网并联，如图3-8所示，根据复合序网可求出电流，电压各序分量：图3-8两相短路接地的复合序网图1图3-8两相短路接地的复合序网图123//EZZZa1ZII0（3-10）a2a1ZZ 2 0ZII2a0a1ZZ 2 0 Z ZUUUI 20EIZ a1 a2 a0 a1Z Z a112 0用对称分量法合成各相电流电压为：IIII0a a1 a2 a0 Z aZ Ia2IaIIIa22 0b a1 a2 a0 a1 Z Z 2 0 Z a2Z IaIa2IIIa2 0（3-11）c a1 a2 a0 a1Z Z 2 0UUUU3Ua a1 a2 a0 a1UU0b c短路点流入地中的电流为：ZIII3I 2 （3-12）b b c a0Z Z 2 0短路点电压，电流的相量图，如图3-9所示。这里任然是纯电感电路。电流滞后电压90图3-9两相接地短路处的电压电流相量图4简单不对称短路的分析与计算计算机计算程序法4.1简单故障的计算程序原理框图﹝1﹞如果要求准确计算故障前的运行情况，则需要进行潮流计算，在近似使用计算中，对短路故障可假设节点U.均为1IGI﹝2﹞这里采用形成节点导纳矩阵的方法。发电机的正序电抗用Xn,可计算故障瞬时量。发电机的负序的电抗近似等于Xn.当计算中不计负何影响时，正可通过其额定功率和电压计算。负何阻抗很难却确定，一般取X=0.35﹝以负何额定功率为基准。负何的中性点一般不接地，零序无通路﹝3﹞形成三序网的节点导纳矩阵后，利用公式﹝3-19﹞和高斯消去法可求得U～U,即为ZZ.对于短路故障，只 1 n 1f～nf需令I=1（其余节电电流均为零），分别应用高斯消去法求解一次所得电压，即为三f序网和f点有关节点阻抗（4）根据不同短路故障，可分别应用表3.1中相应公式计算故障处各序电流，电压，进而合成得三相电流，电压。（5）计算网络中任一点电压时，负序和零电压只须计算由故障点电流引起的电压。对于正序则还需要加上正常运行时的电压表4-1不对出故障处的各序电流电压计算公式。下图4-1即为不对称短路计算原理图：输入数据形成网络节点导纳矩阵 选择短路点 短路点注入单位电流时应用高斯消去法解各点电压即Z....Z 1f nf应用式表（应用式表（4-1）中各种故障对应电流公式计算FI应用节点电压公式计算任一节点电压应用节点电压公式计算任一节点电压结束图4-1不对称短路计算程序框图4.2网络节点方程的形成计算机编程计算中，考虑了对地电容和变压器的实际变比〔7〕导纳计算公式如下：对角线元素：IY（jI,U=1,U=0,ij）（4-1）YiI iiU i ii i ji非对角线元素：YIjy（U=1,U=0,ij）（4-2）ijU ji i ji式中i，j------第i,j节点II，U，U----第i,j节点电流电压变压器支路的导纳 i，，j i j-变化YYT（高压侧）（4-3）iiK2yYT（4-4）ijK Y ..Y4）网络节点方程的形成矩阵表达式.11. 1.n . .. Yn1..Ynn网络节点方程：IYU（4-5）IY11..Y1nU .. .U II. ..U（4-6） IYn1 YnnU5电力系统不对称短路计算实例5.1单相接地短路和两相短路不对称故障分析与计算已知量如图5--1中所示数据，现求取节点K1单相接地短路时和节点K2两相短路时的短路电流，冲击电流以及短路容量。图5-1图5-1系统电路图（一）：解析法分析与计算单相接地短路和两相短路时情况下：短路电流，冲击电流取基准容量：BS100MV.A;基准电压：avBUU（kv）图5-2等值电路图avB等值电路图如图5-2：各原件电抗标幺值计算：取：S100MV.A,UUG1电机电抗标幺值由式（1-4）得：100X0.2640.850.204 11 110正序电抗标幺值X等于负序电抗标幺值X等于零序电抗X即： 11 12 10X=X=X=X G1 111210发电机G2的正序电抗标幺值由式（1-4）得到100X10.130250.800.4165正序电抗标幺值X等于负序电抗标幺值X等于负序电抗X即：X=X=X==0.416.对下标的说各原件下标括的数字表示为第几个G2 515250原件，括号内的数字1，2，0分别为正序，负序，零序。例如，XX12第一个原件（负序）U2负载。负载电抗标幺值X由式（1-5）得XQX;负载（负何）正序 7 LS2 L 7电抗标幺值X6;负序电抗标幺值X 0.35;零X 1.2. 71 72 7010.5100变压器。由式（1-6）得：T变压器正序电抗标幺值X:X 0.656; 1 2 2110016正序电抗标幺值X等于负序电抗标幺值X 等于零序电抗标幺值X 。即 21 22 2010.5100XXX。T变压器正序电抗标幺值：XX0.333;2122202414110031.5即XXX（每台）。T3变压器正序电抗标幺值X: 414240 910.5100X1.05;XXX.9110010919290（4）线路。由式（1-8）得：L1线路的正序电抗标幺值X3:10.5X 0.31000.25;XX;X 3X0.885；L2线路的正 311152 3132303110.5序电抗标幺值X6:X 0.3841000.290;XX ; 61115 616210.5X 3X0.87L3线路的正序电抗标幺值X8:X 0.384500.145;6061 811152XX;X 3X0.435.818280812不对称短路电流的计算：各序等值电路图，正序等值电路图如图5-3所示图5-3图5-3正序等值电路图 X XX X //X X //XX）//XX）=0.22131112131415161718191短路点K2正序标幺值为： X XX //XXX //XX）//XX）=0.33921112131415161718191（1）负序等值电路图：负载的负序电抗标幺值X0.35；其他原件的负序电抗标幺72值等于正序标幺值，如图5-4所示：图5-4图5-4负序等值电路图2827262524232221//////XXXXXXXX X X 32 92短路点K2负序标幺值为： X XX //X X X //X X //X X =0.31322 122232425262728292零序等值电路图：负载的零序电阻抗标幺值X 1.2；线路的零序电抗标幺值等于分70别为：X 3X0.885;X 3X0.87;X 3X0.435.其他原件的零序3031 6061 8081阻抗标幺值等于正序标幺值，如图5-5所示图5-5图5-5零序等值电路图//////XXXXXXXX0.128 30 20304060708090短路点K1零序标幺值为: X X //X X //X X //X X 0.40020 203040607080902K1点单相接地短路K111 1 1 1.934 I 3a1X X X 0.2210.1680.128 31 32 301I131.9345.803IISB5.8031002.913KAIImK 3a 3a1 K K 3U 3115Bi2.55I2.552.9137.4282KA；SSI5.803100508.3MVAsh K K BKK2点两相短路K2：2 2 1 1 1.534I 2a1X X 0.3390.313 21 22IIIm2I21.7321.5342.657 K b c 2a1 S 100KA；i2.55I2.551.3343.402KAIIB2.6571.334K K 3U 3115 sh KBSSI2.657100265.7MVA二计算机程序法析与计算单相接地短路和两相短路时情况下：短路电流，冲击电流以及短路容量：110%U/U1.由式K15%U1/UB1,,K11KT1,K22KT2,K33KT3,得到：11 2 B2 110%110/115 110%110/115 110/115K15%10/10.51.05；K2215%10/10.51.05；K336.3/6.30.96112.正序导纳标幺值Y矩阵Y矩阵的形成：正序电抗标幺值和对地电容标幺值的等值电路图如图5-6所示：图5-6正序电抗标幺值和对地电容标幺值的等值电路图图5-6正序电抗标幺值和对地电容标幺值的等值电路图对地导纳式如下（5-1）：lbSUYBB021 1 11152 YY Y 2.921061000.019 20.0 20"..022100 1 11152 YY Y 2.971061000.020 21.0 21"..022100 1 11152 YY Y 2.97106500.010 22.0 22"..022100其中，Y是图5-6节点1—3得对地导纳之半，Y，Y分别为节点3—5，3—6的 22.0 21..0 22..0y导纳之半【9】.Y矩阵的元素为：Yyyj6.402;Y12j1.45; 11 10 12 12 k11y yYYYYY0;Y12j1.45;Yy21yj4.754 13 14 15 16 17 21 k 22 22.0k2 23 11 22Yyj3.39；YYYY0;Y0;yyyy0.049； 23 23 24 25 26 27 31 30 20.0 21.0 22.0yYyj3.39;Yyy34yyj19.116; 32 23 33 30 32y2 35 3622yY34j5.7;Yyj3.44;8Yyj6.897；Y0；YY0; 34 k 35 35 36 36 37 41 4222YYj3.448;Y0;Yyyyj3.6;YY0;YY0; 53 35 54 55 7 53 21.0 56 57 61 62YYj5.7;Yyyj10.8;YYY0;YY0;43 34 44 34 5 45 46 47 51 52y y图5-7负序阻抗标幺值和对地电纳标幺值的等值电路图同理的得：如图5-7负序阻抗标幺值和对地电纳标幺值的等值电路图，图5-7负序阻抗标幺值和对地电纳标幺值的等值电路图同理的得：如图5-7负序阻抗标幺值和对地电纳标幺值的等值电路图，63 36 64 65 66y2 63 22.0 67k 33 33YYYYY0;YYj0.99;Yj0.95. 71 72 73 74 75 76 67 77 j6.424 j1.45 0 0 0 0 0 j1.45 j4.754 j3.39 0 0 0 0 0 j3.39j19.116 j5.7 j3.448j6.897 0 Y10 0 j5.7 j10.8 0 0 0（5-1） 0 0 j3.448 0 j3.6 0 0 0 0 j6.897 0 0 j7.9 j0.99 0 0 0 0 0 j0.99j0.953负序阻抗Y矩阵及Y负序导纳矩阵如下式如下（5-2）：2 j6.424 j1.45 0 0 0 0 0 j1.45 j4.754 j3.39 0 0 0 0 0 j3.39 j19.116 j5.7 j3.448 j6.897 0 Y20 0 j5.7 j10.8 0 0 0（5-2） 0 0 j3.448 0 j6.2854 0 0 0 0 j6.897 0 0 j7.9 0 0 0 0 0 j0.99j0.99j0.953.零序列Y矩阵：同理零序列电抗标幺值的等值电路图如图5-8所示：图5-8所示零序列电抗标幺值的等值电路图图5-8所示零序列电抗标幺值的等值电路图及Y负序导纳矩阵如下如下式（5-3）：00000000 0 0 0j2.513j1.130 0 j1.13 j100j1.149 j2.30 Y0 0 0 0 0 0 0（5-3）0 0 0 j1.1490j1.983 0 0 0 0 j2.30 0 j3.330 0 0 0 0 0 0 04求不对称短路电流及相关参数：由式3-19及3-21节相关知识，运用高斯消去程序【10】解得：K1点单相接地短路时的转移阻抗为如下（5-4）：ZZ132311jj00..23320526ZZ132322jj00..01366381ZZ132300j00.039Zj0.3046Zj0.2130Z0j0.127331 332 33Z431j0.1607;Z432j0.1124;Z4300 ;（5-4）Z1j0.2917Z532j0.1169Z530j0.05053Z631j0.3058Z632j0.2139Z630j0.060Z731j0.3187Z732j0.2229Z7300K2点两相短路时的转移阻如下式（5-5）：ZZ121jj00..40450913ZZ221212jj00..07922495ZZ122200j0.0380221Zj0.2332Z322j0.1631Z320j0.039j0.1231;Z422j0.0861;Z4200 ;（5-5）Zj0.2234Z522j0.0895Z520j0.022521Z621j0.2342Z622j0.1638Z620j0.027Z721j0.2441Z722j0.1707Z7200 由式表（3-1）相关公式得：K1点单相接地短路K11 1 1 1 j1.956I 33X X X j0.3046j0.2130j0.1275331333330 S 100 KA I5.868;II；II B5.8682.943 3a K 3a K K3U 3115Bi2.55I2.552.9437.5047KA; sh KSSI5.868100586.8MV.A; K B KK2点两相短路K2:2 2 1 1 j1.529；I2I2;I 22X X j0.4045j0.2495 2a1 22221222IIIm2I21.7321.5292.648； K b c 2a1 S 100 KA; II B2.6481.329 K K3U 3115Bi2.55I2.551.3293.389KA; M KSSI2.648100264.8MVA;程序流图输入数据形成网络节点导纳矩阵 选择短路点短路点注入单位电流时应用高斯消去法应用式表（4.1）中各种故解各点电压即Z....Z障对应电流公式计算I. 1f nf F应用公式应用公式1122fffffZIUU222ffffZIU000ffffZIU1100fffffZIUU222ffffZIU计算任意点节点电压应用公式，计算任一点节点电流应用公式，计算任一点节点电流021fffIIIffffffffZZZZU3021)0(fffffffffZZZUII21)0(21应用公式1-9，1-13计算次暂态电流和容量结束结束5-10不对称短路计算程序框图5.2两种计算方法的对比下面将解析法与计算机程序法进行简单比较。将例一和例二结果如表5-1下：表5-1解析法与计算机程序法数据 比较项目故障处短路电流 短路冲击电流（KA） 短路容量（MV.A）短路类型 1 2 1 2 1 2单相接地短路5.8035.8687.42827.5047508.3568.3两相短路2.6572.6483.4023.389265.7264.81--代表解析法结果；2--代表计算机程序法结果。由例一和例二计算过程以及上面的表格数据及比较可得：解析法误差大，计算过程数据大多需要进行四舍五入。计算机程序法中，利用Y矩阵计算时考虑对地电容，变压器实际变比误差小，通过变换式（3-19）中的电流的取值可以将网络各个节点转移阻抗求出，使分析其他点的短路故障提供了更容易更直观的参数值，因此计算机程序法有通用性强的特点。结束语电力系统发生不对称短路后，由于短路点对地故障支路的不对称，使得整个网络电流电压三相不对称：。本论文解决不对称短路的问题核心是对称分量法。根据对称分量法采取的具体方法之一是解析法，即把该网络分解为正，负，零序三个对称序网，这三组对称序分量可分别按对称的三相电路分解。然后将其结果叠加起来。word文档可自由复制编辑word文档可自由复制编辑word文档可自由复制编辑求解不对称短路。首先应该计算各原件的序参数和画出等值电路。然后制定各序网络。根据不同的故障类型，确定出以相分量表示的边界条件，进而列出以序分量表示的边界条件，按边界条件将三个序网联合成复合网，由复合网求出故障处各序电流和电压，进而合成三相电流电压。方法二是计算机程序法。通过计算机形成三个序网的节点导纳矩阵，然后利用高斯消去法通过相应公式对他们进行数据运算，即可求得故障端点的等值阻抗。最后根据故障类型选取相关公式计算故障处各序电流电压，进而合成三相电流电压。通过上面两种方法放入计算进行对比可以得出，计算机程序法比较解析法具有计算过程及结果更精确，通用性更强的优点。参考资料孙丽华．电力工程基础．北京：机械工业出版社，2006李光琦．电力系统暂态分析(第二版)．北京:中国电力出版社,1995于永源，杨绮雯．电力系统分析(第二版)．北京:中国电力出版社,2004刘万顺.电力系统故障分析．北京：中国电力出版社，1999[5]刘从爱．电力工程基础．济南：山东科学技术出版社，1997《电气工程师手册》第二版编辑委员会．电气工程师手册．北京：机械工业出版社，2000西安交通大学，电力工业部西北电力设计院，电力工业部西北勘测设计院．短路电流实用计算方法.北京：电力工业出版社，1982[7]何仰赞温增银，电力系统分析(第三版).武汉华中科技大学出版社，2001刘万顺．电力系统故障分析．北京：中国电力出版社，1998张俊才．电力系统故障的分析计算一对称分量法与相分量法．江西工业大学．邱关源．电路．北京：高等教育出版社，1999罗建军杨琦.MATLAB教程北京：电子工业出版社2005潘小辉陈强．MATLAB5.1全攻略宝典.北京：中国水利水电出版社，2000BillintonRoy.ALLANRONALD(RONALDNORMAN).RELIABILITYEVALUATIONOFPOWERSYSTEMS;NEWYORE;PLENUMPRESS;1996IJ.NAGRATH.D.P.K.OTHARI.MODERN.POWERSYSTEMANALYSIS.NEWDELHI:TATAMCGRAW-HILLCO.TD.1980附录：不对称短路电流计算程序#include“stdio.h”#defineSB100#defineUav115Main(){intK,j,n,d,SFloatx1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,Y[7][7],U[7],I[7],UK,SN,P,Q,X,b,L,Y1,Y2,Y3,K1,K2,K3,IK,IM,SD,PN;K1=K2=1.05;K3=0.96;Clear();电抗.对地导纳Printf(“pleaseinputthepatamentofgenerater:Xd,PN,cosa\n”);Scanf(“%f,%f,%f”,&Xd,&PN,&X);X1=Xd*Sn/(PN/X);Scanf(“%f,%f,%f”,&Xd,&PN,&X);X10=Xd*Sn/(PN/X);Printf(“X1=%f,X10=%f\n”,X0,X10);Printf(“pleaseinputtheputamenteroftramsfer:UKandSN\n1:”);Scanf(“%f,%f,”,&UK,&SN);X2=UK*Sn/SN;Printf(“\n2:\n”);Scanf(“%f,%f,”,&UK,&SN);X3=UK*Sn/SN;Printf(“\n3:\n”);Scanf(“%f,%f,”,&UK,&SN);X4=UK*Sn/SN;Printf(“X2=%f,X3=%f,X4=%f\n”,X2,X3,X4);Printf(“pleaseinputtheputamenterofload:P,Q\n”);Scanf(“%f,%f,”,&P,&Q);X8=1/((Q/100)*Q/100+(P/100)*(P/100))*Q/100);X11=0.35;X12=1.2;Printf(“pleaseinputtheMofload:P,cosa\n”);Scanf(“%f,%f,”,&P,&X);X9=1/6.5*100/(P/X);Printf(“X8=%f,X9=%f\n”,X8,X9);Printf(“pleaseinputtheputamenterofline:\n”);Printf(“forexample:X=0.39,b=02.92e-6,i=100”);Printf(“\n1:”);Scanf(“%f,%f,”,&X,b,&&L);X5=X*L*Sn/(Uav*Uav);Y1=(double)1\2*b*L*(Uav*Uav)/100;Printf(“\n2:”);Scanf(“%f,%f,”,&X,b,&&L);X6=X*L*Sn/(Uav*Uav);Y2=(double)1\2*b*L*(Uav*Uav)/100;Printf(“\n3:”);Scanf(“%f,%f,”,&X,b,&&L);X7=X*L*Sn/(Uav*Uav);Y3=(double)1\2*b*L*(Uav*Uav)/100;Printf(“X5=%f,X6=%f,X7=%f\n”,X5,X6,X7);Printf(“Y1=%f,Y2=%f,Y3=%f\n”,Y1,Y2,Y3);Printf(“Y1=Y1,Y2=Y2,Y3=Y3\n”,Y1,Y2,Y3);Printf(“pleaseinputtheK(7)toseeYarray:\n”);Scanf(“&d”,&K);