第17章 电力系统潮流计算的计算机算法

1、本章提示19.1概述19.2潮流计算的基本方程19.3牛顿拉夫逊法潮流计算19.4PQ分解法潮流计算小结,第19章电力系统潮流计算的计算机算法,节点分类的概念；潮流计算的基本方程式；牛顿拉夫逊法潮流计算的计算机算法；PQ分解法潮流计算的计算机算法。,本章提示,类型：,19.1概述,19.2.1节点的分类19.2.2基本方程式,19.2潮流计算的基本方程,根据电力系统中各节点性质的不同，可把节点分成三种类型。,1.PQ节点事先给定的是节点功率（P、Q），待求的是节点电压向量（U、）。通常变电所母线都是PQ节点，当某些发电机的出力P、Q给定时，也可作为PQ节点。PQ节点上的发电机称之为PQ机（或P

2、Q给定型发电机）。在潮流计算中，系统大部分节点属于PQ节点。,19.2.1节点的分类,2.PU节点给出的参数是节点的有功功率P及电压幅值U，待求量为该节点的无功功率Q及电压向量的相角。通常选择有一定无功功率贮备的发电机母线或者有无功补偿设备的变电所母线作PU节点。PU节点上的发电机称之为PU机（或PU给定型发电机）。,3.平衡节点给定的运行参数是U和，而待求量是该节点的P、Q，因此又称为U节点。在潮流计算中，这类节点一般只设一个。,关于平衡节点的选择，一般选择系统中担任调频调压的某一发电厂（或发电机），有时也可能按其它原则选择。,任何复杂的电力系统都可以归结为以下元件（参数）组成：（1）发电机

3、（注入电流或功率）；（2）负荷（负的注入电流或功率）；（3）输电线支路（电阻、电抗）；（4）变压器支路（电阻、电抗、变比）；（5）母线上的对地支路（阻抗和导纳）；（6）线路上的对地支路（一般为线路充电电容导纳）。,19.2.2基本方程式,集中了以上各种类型元件的简单网络如图,19.2.2基本方程式,节点注入电流和节点电压构成以下线性方程组,其中,19.2.2基本方程式,节点功率与节点电流之间的关系为：,（19.3）,式中,PQ节点可以表示为,（19.4）,把这个关系式代入式（19.1）中，得,（19.5）,式（19.5）是一组共有n个非线性方程组成的复数方程式，如果把实部和虚部分开便得到2n个

4、实数方程，因此由该方程组可解出2n个运行参数。,19.2.2基本方程式,19.3.1牛顿拉夫逊法概要19.3.2牛顿拉夫逊法潮流计算19.3.3牛顿法的框图及求解过程19.3.4实例,19.3牛顿-拉夫逊法潮流计算,已知一个变量X的函数为：,（19.6）,解此方程式时，由适当的近似值X（0）出发，根据,（19.7）,反复进行计算，当X（n）满足适当的收敛判定条件时就是（19.6）式的根。这样的方法就是所谓的牛顿拉夫逊法。,19.3.1牛顿-拉夫逊法概要,几何意义:,图19.2函数曲线及切线示意图,19.3.1牛顿-拉夫逊法概要,用同样的方法考虑，给出对n个变量的n个方程式,（19.12）,对其

5、近似解的修正量可以解下面的方程式来确定,（19.13）,式（19.13）等号右边矩阵的等都是对于的值，这一矩阵称为雅可比（Jacobi）矩阵。,按上述得到修正量后，得到如下关系：,这比进一步接近于真值。这一步骤在收敛到希望的值以前重复进行。一般要反复计算到满足,时为止。为预先规定的小正数，此处是第n次迭代Xi的近似值。,19.3.1牛顿-拉夫逊法概要,1采用直角坐标结点电压和导纳可表示为：,将上述表示式代入的右端，展开并分出实部和虚部，便得：,（19.14）,19.3.2牛顿-拉夫逊法潮流计算,PQ节点的有功功率和无功功率是给定的，第i个节点的给定功率设为Pis和Qis。假定系统中的第1，2，

6、m号节点为PQ节点，对其中每一个节点可列方程,（i=1，2，m）（19.15）,19.3.2牛顿-拉夫逊法潮流计算,PU节点的有功功率和节点电压幅值是给定的，假定系统中的第m+1，m+2，n-1号节点为PU节点，则对其中每一节点可以列写方程：,（19.16）,第n号节点为平衡节点，其电压是给定的，故不参加迭代。,式（19.15）和（19.16）总共包含了2（n-1）个方程，待求的变量有也是2（n-1）个。,19.3.2牛顿-拉夫逊法潮流计算,方程式（19.15）和（19.16）具备方程组（19.12）的形式：,（19.17）,式中,19.3.2牛顿-拉夫逊法潮流计算,雅可比矩阵的各元素是对式（

7、19.15）和（19.16）求偏导数,当时，对角元素是,（19.18）,19.3.2牛顿-拉夫逊法潮流计算,当时，矩阵中非对角元素是,（19.19）,19.3.2牛顿-拉夫逊法潮流计算,雅可比矩阵有以下特点：雅可比矩阵中的诸元素都是节点电压的函数，因此在迭代过程中，它们将随着各节点电压的变化而不断地改变；矩阵是不对称的；由式（19.19）可以看出，当导纳矩阵中的非对角元素Yij为零时，雅可比矩阵中相对应的元素也是零，即矩阵是非常稀疏的。,19.3.2牛顿-拉夫逊法潮流计算,2.采用极坐标结点电压和导纳可以表示为,（19.20）,将式（19.20）代入右端并将实部与虚部分开，得,（19.21）,

8、式中为i、j两结点电压相角差（）。,19.3.2牛顿-拉夫逊法潮流计算,雅可比矩阵中各元素可对式（19.21）取偏导数求得,19.3.2牛顿-拉夫逊法潮流计算,（19.22）,用分块矩阵的形式简化如下,修正方程为,（19.23）,19.3.3牛顿法的框图及求解过程,用牛顿法计算潮流时，有以下的步骤：给出各节点电压初始值将以上电压初始值代入式（19.15）和（19.16），求出修正方程式的常数项向量将电压初始值再代入式（19.18）和（19.19），求出修正方程式中系数矩阵（雅可比矩阵）的各元素解修正方程式（19.16），求出修正量修正各节点电压,19.3.3牛顿法的框图及求解过程,将再代入（1

9、9.15），（19.16）式求校验是否收敛，即如果收敛，迭代到此结束，进一步计算各线路潮流和平衡节点功率，并打印输出结果。如果不收敛，转回第（2）步进行下一次迭代计算，直到收敛为止。,19.3.3牛顿法的框图及求解过程,例19.1试用牛顿拉夫逊法计算图19.4所示电力系统的潮流分布。,图19.4例19.1网络图,19.3.4实例,解：1.需要输入的数据n节点数、nl支路数、isb平衡母线节点号（固定为1）、pr误差精度、如果要输入则输入eps即可。请输入由支路参数形成的矩阵B1矩阵B1的每行是由下列参数构成的：某支路的首端号P；某支路末端号Q，且PQ；支路的阻抗（R+jX）；支路的对地容抗；支

10、路的变比K；折算到哪一侧的标志（如果支路的首端P处于高压侧则请输入“1”，否则请输入“0”）。,19.3.4实例,请输入各节点参数形成的矩阵B2矩阵B2的每行是由下列参数构成的：节点所接发电机的功率SG；节点负荷的功率SL；节点电压的初始值。PU节点电压U的给定值。节点所接的无功补偿设备的容量。节点分类标号igl。1平衡节点igl=2PQ节点3PU节点（4）请输入由节点号及其容抗形成的矩阵X,19.3.4实例,2.先形成节点导纳矩阵。3.根据式（19.15）和（19.16）求出修正方程式的常数项向量4.据式（19.18）和（19.19）求出雅可比矩阵各元素值，即可得到第一次迭代时的修正方程式。

11、5.解方程式（19.16），求6.修正各节点电压，即得出第一次迭代后各节点的电压值。7.按以上计算步骤迭代下去，当收敛精度取（即PR=0.0001）时，需要进行四次迭代。求出了各节点电压后，即可求各支路的潮流分布。,19.3.4实例,针对例19.1输入数据如下：请输入节点数:n=5请输入支路数:nl=5请输入平衡母线节点号:isb=1请输入误差精度:pr=0.00001请输入由支路参数形成的矩阵:B1=120.03i01.050;230.08+0.3i0.5i10;240.1+0.35i010;340.04+0.25i0.5i10;350.015i01.051请输入各节点参数形成的矩阵:B2=

12、001.051.0501;03.7+1.3i1002;02+1i1002;01.6+0.8i1002;501.051.0503请输入由节点号及其对地阻抗形成的矩阵:X=10;20;30;40;50,19.3.4实例,结果如下：迭代次数5没有达到精度要求的个数78860各节点的实际电压标么值E为(节点号从小到大排列):1.05001.0335-0.0774i1.0260+0.3305i0.8592-0.0718i0.9746+0.3907i各节点的电压大小U为(节点号从小到大排列):1.05001.03641.07790.86221.0500各节点的电压相角为(节点号从小到大排列):0-4.28

13、1917.8535-4.778521.8433各节点的功率S为(节点号从小到大排列):2.5794+2.2994i-3.7000-1.3000i-2.0000-1.0000i-1.6000-0.8000i5.0000+1.8131i,19.3.4实例,各条支路的首端功率Si为(顺序同您输入B1时一样):2.5794+2.2994i-1.2774+0.2032i0.1568+0.4713i1.5845+0.6726i5.0000+1.8131i各条支路的末端功率Sj为(顺序同您输入B1时一样):-2.5794-1.9745i1.4155-0.2443i-0.1338-0.3909i-1.4662

14、-0.4091i-5.0000-1.4282i各条支路的功率损耗DS为(顺序同您输入B1时一样):-0.0000+0.3249i0.1381-0.0412i0.0230+0.0804i0.1184+0.2635i0.0000+0.3849i,19.3.4实例,每次迭代后各节点的电压值如图,19.3.4实例,19.4.1PQ分解法的基本方程式19.4.2计算步骤和程序框图19.4.3实例,19.4PQ分解法潮流计算,（19.24）,（19.25）,快速分解法（又称PQ分解法）是从简化牛顿法极坐标形式计算潮流程序的基础上提出来的。节点功率方程在用极坐标形式表示时，它的修正方程式可简化为：,在一般情

15、况下，线路两端电压的相角是不大的（不超过1020）。因此，,19.4.1PQ分解法的基本方程式,与系统各节点无功功率相应的导纳B远远小于该节点自导纳的虚部，即,式（19.24）的系数矩阵中的各元素可表示为：,（i,j=1,2,，n-1）（19.27）,（i,j=1,2,，m）(19.28),19.4.1PQ分解法的基本方程式,而系数矩阵H和L则可以分别写成：,将（19.29）和（19.30）式代入（19.24）中，得到,（19.31）,（19.32）,这就是简化了的修正方程式，它们也可展开写成：,（19.33）,（19.34）,用和分别左乘以上两式便得：,19.4.1PQ分解法的基本方程式,在

16、这两个修正方程式中系数矩阵元素就是系统导纳矩阵的虚部，因而系数矩阵是对称矩阵，且在迭代过程中保持不变。这就大大减少了计算工作量。用极坐标表示的节点功率增量为：,（19.35）,式（19.33）、（19.34）和（19.35）构成了PQ分解法迭代过程的基本方程式。,19.4.1PQ分解法的基本方程式,（1）给定各节点电压的初始值（2）代入式（19.35）计算各节点有功功率，并求出（3）解修正方程式（19.33），得出各节点电压相角修正量（4）修正各节点电压的相角，（5）代入式（19.35）计算各节点无功功率误差，并求出（6）解修正方程式（19.34），得出各节点电压幅值的修正量；（7）修正各节点

17、电压的幅值，（8）返回（2）进行迭代，直到各节点功率误差及都满足收敛条件,19.4.2计算步骤和程序框图,例19.2用PQ分解法计算图19.4所示网络的潮流分布。输入数据为：请输入节点数:n=5请输入支路数:nl=5请输入平衡母线节点号:isb=1请输入误差精度:pr=0.00001请输入由支路参数形成的矩阵:B1=120.03i01.050;230.08+0.3i0.5i10;240.1+0.35i010;340.04+0.25i0.5i10;350.015i01.051请输入由节点参数形成的矩阵:B2=001.051.0501;03.7+1.3i1.05002;02+1i1.05002;0

18、1.6+0.8i1.05002;501.051.0503请输入由节点号及其对地阻抗形成的矩阵:X=10;20;30;40;50请输入PQ节点数na=3结果如下：迭代次数10,19.4.3实例,每次没有达到精度要求的有功功率个数为4444444430每次没有达到精度要求的无功功率个数为3333333310各节点的电压标么值E为(节点号从小到大排列):1.05001.0335-0.0774i1.0260+0.3305i0.8592-0.0718i0.9746+0.3907i各节点的电压大小U为(节点号从小到大排列):1.05001.03641.07790.86221.0500各节点的电压相角为(节点号从小到大排列):0-4.281917.8535-4.778521.8433各节点的功率S为(节点号从小到大排列):2.5794+2.2994i-3.7000-1.3000i-2.0000-1.0000i-1.6000-0.8000i5.0000+1.8131i各条支路的首端功率S