复杂电力系统潮流计算.doc

0华侨大学厦门工学院电力系统综合设计 课程设计报告题 目：复杂电力系统潮流计算专业、班级： 09 级电气(4)班 学生姓名： 学 号： 指导教师： 分 数 : 2012 年 5 月 26 日1目　录第一章 牛顿拉夫逊算法的基本资料.....................................21.1 牛顿拉夫逊算法定义..........................................21.2 牛顿拉夫逊算法法的发展与前景...............................第二章 电力网络的数学模型..........................................32.1 节点导纳矩阵的形成及修改....................................32.1.1 节点导纳矩阵的形成........................................32.1.2 节点导纳矩阵的修改........................................52.2 节点导纳矩阵元素的物理意义..................................7第三章 计算实例....................................................93.1 等值电路图.................................................113.2 节点导纳矩 阵...............................................113.3 设定所求变量的初值.........................................123.4 计算修正方程...............................................133.5 形成雅可比矩阵.............................................153.6 求解修正方程...............................................163.7 进行修正和迭代.............................................16 3.8 迭代精度的确认.............................................173.9 各节点电压计算功率分布.....................................17结　 论.............................................................19致 谢.............................................................19参考文献...........................................................202摘 要本次的课程设计主要针对复杂电力系统——用牛顿-拉夫逊法来进行潮流计算.牛顿-拉夫逊法对初值要求严格,迭代速度快的特点,利用电力网的结构特点,提出直角坐标和极坐标牛顿 -拉夫逊法潮流计算的三元素解法及相应的简化算法 ,并对其进行计算分析比较占用内存少,计算量小,且不影响其收敛性及准确性计算结果表明,综合算法在迭代次数和收敛速度上有优势。

关键词：牛顿-拉夫逊法 收敛迭代 潮流计算第一章 牛顿拉夫逊算法基础资料1、牛顿-拉夫逊法定义：牛顿迭代法（Newton's method）又称为牛顿-拉夫逊方法（Newton-Raphson method） ，它是牛顿在 17 世纪提出的一种在实数域和复数域上近似求解方程的方法多数方程不存在求根公式，因此求精确根非常困难，甚至不可能，从而寻找方程的近似根就显得特别重要方法使用函数 f(x)的泰勒级数的前面几项来寻找方程 f(x) = 0 的根牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程 f(x) = 0 的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根2、牛顿-拉夫逊法现状与前景： 利用电子计算机进行潮流计算从 20 世纪 50 年代中期就已经开始此后，潮流计算曾采用了各种不同的方法，这些方法的发展主要是围绕着对潮流计算的一些基本要求进行的对潮流计算的要求可以归纳为下面几点： （1）算法的可靠性或收敛性 （2）计算速度和内存占用量 （3）计算的方便性和灵活性 20 世纪 60 年代初，数字计算机已经发展到第二代，计算机的内存和计算速度发生了很大的飞跃，从而为阻抗法的采用创造了条件。

阻抗矩阵是满矩阵，阻抗法要求计算机储存表征系统接线和参数的阻抗矩阵这就 需要较大的内存量阻抗法改善了电力系统潮流计算问题的收敛性，解决了导纳法无法解决的一些系统的潮流计算，但是，阻抗法的主要缺点就是占用计算机的内存很大，3每迭代的计算量很大当系统不断扩大时，这些缺点就更加突出近 20 多年来，潮流算法的研究仍然非常活跃，但是大多数研究都是围绕改进牛顿法和 P-Q 分解法进行的此外，随着人工智能理论的发展，遗传算法、人工神经网络、模糊算法也逐渐被引入潮流计算但是，到目前为止这些新的模型和算法还不能取代牛顿法和 P-Q 分解法的地位由于电力系统规模的不断扩大，对计算速度的要求不断提高，计算机的并行计算技术也将在潮流计算中得到广泛的应用，成为重要的研究领域通过几十年的发展，潮流算法日趋成熟近几年，对潮流算法的研究仍然是如何改善传统的潮流算法，即高斯-塞德尔法、牛顿法和快速解耦法牛顿法，由于其在求解非线性潮流方程时采用的是逐次线性化的方法，为了进一步提高算法的收敛性和计算速度，人们考虑采用将泰勒级数的高阶项或非线性项也考虑进来，于是产生了二阶潮流算法后来又提出了根据直角坐标形式的潮流方程是一个二次代数方程的特点，提出了采用直角坐标的保留非线性快速潮流算法 【6】 。

第二章 电力网络的数学模型 线性网络的常用解法有节点电压法和回路法，前者须列写节点电流平衡方程，后者则须列写回路方程本章重点介绍节点方程，以及节点导纳矩阵 【1】 2.1 节点导纳矩阵的形成及修改2.1.1 节点导纳矩阵的形成在图 2-1（a）的简单电力系统中，若略去变压器的励磁功率和线路电容，负荷用阻抗表示，便可以得到一个有 5 个节点（包括零电位点）和 7 条支路的等值网络，如图 2-1（b）所示将接于节点 1 和 4 的电势源和阻抗的串联组合变换成等值的电流源和导纳的并联组合，变得到图（c）的等值网络，其中和 分别称为节点 1 和 4 的注入电流源10IyE&404Iy&41243( a )Ė1Ė41234y1 0y1 2y2 4y2 0y2 3y3 4y4 0( b )İ1İ41234y1 2y2 4y2 3y3 4y4 0y2 0y ´1 0( c )图 2-1 电力系统及其网络以零电位点作为计算节点电压的参考点，根据基尔霍夫定律，可以写出 4个独立节点的电流平衡方程如下（ 2-1）10212103232442323444304()()()0()yUIyUyUy I&&&上述方程组经过整理可以写成（ 2-2）121324243440YUIYUI &&式中， ； ； ；1012Yy20321Yyy3234Yy； ； ； ；40243y142。

3一般的，对于有 个独立节点的网络，可以列写 个节点方程nn5（ 2-3）12112212nnnnYUYIUI&&LM也可以用矩阵写成（ 2-4）112212nnnnUIYI&LM或缩写为（2-5）YUI矩阵 称为节点导纳矩阵它的对角线元素 称为节点 的自导纳，其值YiYi等于接于节点 的所有支路导纳之和非对角线元素 称为节点 、 间的互导i ij j纳，它等于直接接于节点 、 间的支路导纳的负值若节点 、 间不存在直ij i接支路，则有 由此可知节点导纳矩阵是一个稀疏的对称矩阵0ijY2.1.2 节点导纳矩阵的修改在电力系统中，接线方式或运行状态等均会发生变化，从而使网络接线改变比如一台变压器支路的投入或切除，均会使与之相连的节点的自导纳或互导纳发生变化，而网络中其它部分结构并没有改变，因此不必从新形成节点导纳矩阵，而只需对原有的矩阵作必要的修改就可以了现在几种典型的接线变化说明具体的修改方法yij( a )iijyij( b )ij- yij( c )k * yT( 1 - K * ) yTk * ( k * - 1 ) yTijI 侧I I 侧( e )ij- yijy ´ij( d )6图 2-2 电力接线的改变（a）增加支路和节点；（b）增加支路；（c）切除支路；（d）改变支路参数；（e）改变变压器变比（1）从原有网络的节点 引出一条导纳为 的支路， 为新增加的节点，iijYj如图 2-2(a)所示。

由于新增加了一个节点，所以节点导纳矩阵增加一阶，矩阵作如下修改：1）原有节点 的自导纳 的增量 ＝ ；iiYiijY2）新增节点 的自导纳 ；jjijy3）新增的非对角元素 ；其它新增的非对角元均为零ijjiij（2）在原有网络的节点 与 j 之间增加一条导纳为 的支路，如图 2-2(b)所ijy示则与 、 有关的元素应作如下修改：ij1）节点 、 的自导纳增量； ij ijijY2）节点 、 的互导纳增量 ijjiijy（3）在网络的原有节点 、 之间切除一条导纳为 的支路，如图 2-ij2（c）所示，其相当在 、 之间增加一条导纳为 的支路，因此与 、 有关ij ijyij的元素应作以下修改：1）节点 、 的自导纳增量 ；ijijijY2）节点 、j 之间的互导纳增量 ；ijjiijy（4）原有网络节点 、 之间的导纳由 变成 ，相当于在节点 、 之ijijijij间切除一条导纳为 的支路，再增加一条导纳为 的支路，如图 2-2（d）所ijyijy示则与 、 有关的元素应作如下修改：ij1）节点 、 的自导纳增量 ；ijijijijY2）节点 、 的互导纳增量 。

ijjiijijy7（5）原有网络节点 、 之间变压器的变比由 变为 ，即相当于切除ij k一台变比为 的变压器，再投入一台变比为 的变压器，k ，如图 2-2（e）变压器Ⅱ型等值电路，图中 为与变压器())BUkCC Ty原边基准电压对应的变压器导纳标幺值，则与 、 有关的元素应作如下修改：ij1）节点 的自导纳增量 ；节点 的自导纳增量 ；i 0iYj2()jTYk2）节点 与 之间的互导纳增量 j ()ijjiTYky2.2 节点导纳矩阵元素的物理意义节点导纳矩阵的元素已在上一节作了说明，现在进一步讨论这些元素的物理意义如果令0kU&j(1,2,)jnjkL代入 2-3 的各式，可得ikiYI(,)i或（ 2-6）0,jiikkUkIY&当 时，公式 2-6 说明，当网络中除节点 以外所有节点都接地时，从ii节点 注入网络的电流同施加于节点 的电压之比，即等于节点 的自导纳 i iiY换句话说，自导纳 是节点 以外的所有节点都接地时节点 对地的总导纳显iYi然， 应等于与节点 相接的各支路导纳之和，即 iY（ 2-7）0iiijjy式中 ， 为节点 与零电位节点之间的支路导纳； 为节点 与节。