导纳矩阵课程设计

南昌工程学院本科课程设计1第一章引言导纳矩阵计算是电力系统分析最基本的计算。除它自身的重要作用之外，在《电力系统分析综合程序》(PSASP)中，还是网损计算、静态安全分析、暂态稳定计算、小干扰静态稳定计算、短路计算、静态和动态等值计算的基础。对规划中的电力系统，通过计算可以检验所提出的电力系统规划方案能否满足各种运行方式的要求；对运行中的电力系统，通过计算可以预知各种负荷变化和网络结构的改变会不会危及系统的安全，系统中所有母线的电压是否在允许的范围以内，系统中各种元件(线路、变压器等)是否会出现过负荷，以及可能出现过负荷时应事先采取哪些预防措施等。1.1课程设计目的与要求1.1.1设计目的1.掌握电力系统潮流计算导纳矩阵的基本原理；2.掌握并能熟练运用一门计算机语言（MATLAB语言或C语言或C++语言）；3.采用计算机语言对潮流计算导纳矩阵进行计算机编程。1.1.2设计要求1.程序源代码；2.给定题目的输入，输出文件；3.程序说明；4.给定系统的程序计算过程；5.给定系统的手算过程（至少迭代2次）。1.2设计内容1.根据电力系统网络推导电力网络数学模型，写出节点导纳矩阵；2.赋予各节点电压变量（直角坐标系形式）初值后，求解不平衡量；3.形成雅可比矩阵；4.求解修正量后，重新修改初值，从2开始重新循环计算；5.求解的电压变量达到所要求的精度时，再计算各支路功率分布、功率损耗和平衡节点功率；6.上机编程调试；7.计算分析给定系统导纳矩阵分析并与手工计算结果做比较分析；8.书写课程设计说明书，准备答辩。第二章导纳矩阵计算概述2第二章导纳矩阵计算概述2.1导纳矩阵计算简介电力系统潮流导纳矩阵计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压。各元件中流过的功率，系统的功率损耗等等。在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量的分析比较供电方案或运行方式的合理性。可靠性和经济性。此外，电力系统的潮流计算也是计算机系统动态稳定和静态稳定的基础，所以潮流计算是研究电力系统的一种和重要和基础的计算。电力系统潮流计算也分为离线计算和在线计算两种，前者主要用于系统规划设计和安排系统的运行方式，后者则用于正在运行系统的经常监视及实时控制。利用电子数字计算机进行潮流计算从50年代中期就已经开始了。在这20年内，潮流计算曾采用了各种不同的方法，这些方法的发展主要围绕着对潮流计算的一些基本要求进行的，对潮流计算的要求可以归纳为以下几点：（1）计算方法的可靠性或收敛性；（2）对计算机内存量的要求；（3）计算速度；（4）计算的方便性和灵活性。2.2导纳矩阵计算的意义(1)在电网规划阶段,通过潮流计算,合理规划电源容量及接入点,合理规划网架,选择无功补偿方案,满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求。(2)在编制年运行方式时,在预计负荷增长及新设备投运基础上,选择典型方式进行潮流计算,发现电网中薄弱环节,供调度员日常调度控制参考,并对规划、基建部门提出改进网架结构,加快基建进度的建议。(3)正常检修及特殊运行方式下的潮流计算,用于日运行方式的编制,指导发电厂开机方式,有功、无功调整方案及负荷调整方案,满足线路、变压器热稳定要求及电压质量要求。(4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。南昌工程学院本科课程设计3第三章导纳矩阵设计题目3.1导纳矩阵计算题目表3-1在100MVA的基础上6条母线的系统图支路母线节点电阻电抗比率1120.0000.3001.0252140.0970.4073160.1230.5184250.2820.6405350.7231.0506430.0000.1331.1007460.0800.370图3-1系统图该系统图是一个6节点，7支路的电力网络，通过手工计算和计算机计算得到该系统图的导纳矩阵，计算结果进行比较，观察两次运算结果是否一样。55.0+j13.050.0+j5.030.0+j18.0P5=50.16125431.1:11.025:1第三章导纳矩阵设计题目43.2对课题的分析及求解思路此电力系统是一个6节点，7支路的电力网络。综合比较牛顿拉夫逊法（直角坐标、极坐标）、PQ分解法等多种求解方法的特点，最后确定采用牛顿拉夫逊法（极坐标）。因为此方法所需解的方程组最少。通过对本课程设计的理解分析，将整个系统的简化图分析如图3-2：图3-2系统简化图南昌工程学院本科课程设计5第四章导纳矩阵计算算法及手工计算4.1节点电压方程在电路中我们学过利用节点电压方程求解某几条支路的电流，现以下图4-1与图4-2为例推导节点电压方程组。图4-1节点电压法为例图4-2用电流源代替电压源为例图4-1表示了一个具有两个电源和一个等值负荷的系统。1e、2e为电源电势，1y、2y为电源的内部导纳，3y为负荷的等值导纳，4y、5y、6y为各支路的导纳。如果以地为电压参考点，设节点1、2、3的电压为，根据基尔霍夫电流KCL法对节点1、2、3列节点电流方程得式（4-1）：......1112113111......2162152322.....431532330yVyVVyVVyyVyVVyVVyyVVyVVyVee（4-1）上式中左端为节点1、2、3流出的电流，右端为注入个节点的电流。由上式可以得到一个等效的等值电路图4-1。图4-2中利用了电流源代替的电压源。在图4-2中可知的式（4-2）：..111..222IyIyee（4-2）第四章导纳矩阵计算算法及手工计算6为等值电流源向网络注入的电流。将与式（4-1）联立得式（4-3）：....1111221331....2112222332....3113223333YVYVYVIYVYVYVIYVYVYVI(4-3）上式中111462225633346,,YyyyYyyyYyyy称为节点1、2、3的自导纳，122161331423325,,YYyYYyYYy称为相应节点之间的互导纳。因此，在一般情况下，在电力网络中有n个节点，则可以按式（4-3）的形式列出n个节点方程式，也可用矩阵的形式表示IYV。其中.1.2....nIIII.1.2....nVVVV(4-4）分别为节点注入电流列向量及节点电压列向量；111212122212.....................nnnnnnYYYYYYYYYY(4-5）南昌工程学院本科课程设计7为节点导纳矩阵，其中对角元素iiY为节点i的自导纳，非对角线ijY为节点i与节点j之间的互导纳。4.2节点导纳矩阵节点导纳矩阵既可根据自导纳和互导纳的定义直接求取，也可根据电路知识中找出改网络的关联矩阵，在节点电压方程的矩阵形式进行求解。本章节我们主要讨论的是直接求解导纳矩阵。根据节点电压方程章节我们知道，在利用电子数字计算机计算电力系统运行情况是，多采用IYV形式的节点方程式。其中阶数等于电力网络的节点数。从而可以得到n个节点时的节点导纳矩阵方程组（4-6）如下：....11112211....211222232....11223.........................nnnnnnnnYVYVYVIYVYVYVIYVYVYVI(4-6）由此可以得到n个节点导纳矩阵：111212122212...............nnnnnnYYYYYYYYYY它反映了网络的参数及接线情况，因此导纳矩阵可以看成是对电力网络电气特性的一种数学抽象。由导纳短阵所联系的节点方程式是电力网络广泛应用的一种数学模型。通过上面的讨论，可以看出节点导纳矩阵的有以下特点：（1）导纳矩阵的元素很容易根据网络接线图和支路参数直观地求得，形成节点导纳矩阵的程序比较简单。（2）导纳矩阵为对称矩阵。由网络的互易特性易知ijjiYY。（3）导纳矩阵是稀疏矩阵。它的对角线元素一般不为零，但在非对角线元素中则存在不少零元素。在电力系统的接线图中，一般每个节点与平均不超过3~4个其他节点有直接的支路连接。因此，在导纳矩阵的非对角线元素中每行仅有3~4个非零元素，其余的都是零元素，而且网络的规模越大，这种现象越显著。导纳矩阵的对称性和稀疏性对于应用计算机求解电力系统问题有很大的影响。如果能充分地利用这两个第四章导纳矩阵计算算法及手工计算8特点，如在程序设计中储存导纳矩阵的对角元素和上三角元素（或下三角元素），排除零元素的储存和运算，就可以大大地节省储存单元和提高计算速度。节点导纳矩阵的形式可归纳如下：（1）导纳矩阵的阶数等于电力网络的节点数。（2）导纳矩阵各行非对角元素中非零元素的个数等于对应节点所连得不接地支路数。（3）导纳矩阵各对角元素，即节点的自导纳等于相应节点之间的支路导纳之和。（4）导纳矩阵非对角元素，即节点之间的互导纳等于相应节点之间的支路导纳的负值。而在电力系统中进行潮流计算时，往往要计算不同接线下的运行状况，例如，改变变压器主抽头时，潮流分布也随之变化，以及改变其他设备参数进行计算潮流分布，此时就需要导出变化时的导纳矩阵就需要对所设计的程序进行参数设定，而不需要重复上述步骤去导出所求的导纳矩阵。4.3变压器的∏型等值电路在电力系统潮流计算中，往往要计算节点导纳矩阵，而我们计算节点导纳矩阵采用节点电压法来实现，如在变压器构成的电力系统中，需要将变压器模型转变成变压器∏型等值电路，在利用电路知识列节点电压方程，从而导出所需的导纳矩阵。1U2U1:kTY/TYk2)1(kYkTkYkT)1(图4-3双绕组变压器的∏型等值电路（i，j为节点）而在电力系统潮流计算中一般采用标幺值进行计算，标幺值公式如下：有名值（任意单位）标幺值基准值（与有名值同单位）南昌工程学院本科课程设计9所以，变压器转变成∏型等值电路时，我们采用标幺值计算，使所求参数为变压器变比k的函数。而在一个已经归算好的电力系统网中，若改变变压器的分接头来进行调压，这时变压器的等值电路参数也会相应得改变，此时采用∏型等值电路进行折算就显得较为方便。下面是变压器的∏型等值电路分析过程：如不计如果采用标么值计算，元件参数都应归算到同一基准值时得标么值，才能在同一个等值电路上分析和计算。励磁支路的影响，双绕组变压器可用其阻抗与一个理想变压器串联的电路表示，如图所示。理想变压器只有一个参数，那就是变比12kUU。现以变压器阻抗按实际变比归算到低压侧的情况为例，推导出双绕组变压器的∏型等值电路。流入和流出理想变压器的功率相等：....1212/kUIUI(1U、2U分别为变压器高、低绕组的实际电压）（4-7）..12/kII（4-8）联立（4-7）、（4-8）两个公式解得：.....1212122TTTTYUYUUUIkZkZkk（4-9）.....11222TTTTYUUUIYUkZZk（4-10）根据《电路原理》节点1、2的节点电流方程具有如下形式：...1121112...2122122IYUYUIYUYU（4-11）将式（4-9）、（4-10）与式（4-11）比较得(4-12):211T12T21T12TYY/kYY/kYY/kYY（4-12）因此可以的得到各支路导纳为第四章导纳矩阵计算算法及手工计算101212T2121T2101112TTT2202221TTTYyY/kYyY/k1kyYyY/kY/kYkk1yYyYY/kYk（4-13）4.4导纳矩阵的手工计算在计