电力系统潮流计算.docx

1、电力系统潮流计算第四章 电力系统潮流分析与计算电力系统潮流计算是电力系统稳态运行分析与控制的基础，同时也是安全性分析、稳 定性分析电磁暂态分析的基础（稳定性分析和电磁暂态分析需要首先计算初始状态，而初 始状态需要进行潮流计算）。其根本任务是根据给定的运行参数，例如节点的注入功率，计 算电网各个节点的电压、相角以及各个支路的有功功率和无功功率的分布及损耗。潮流计算的本质是求解节点功率方程，系统的节点功率方程是节点电压方程乘以节点 电压构成的。要想计算各个支路的功率潮流，首先根据节点的注入功率计算节点电压，即 求解节点功率方程。节点功率方程是一组高维的非线性代数方程，需要借助数字迭代的计 算方法来

2、完成。简单辐射型网络和环形网络的潮流估算是以单支路的潮流计算为基础的。本章主要介绍电力系统的节点功率方程的形成，潮流计算的数值计算方法，包括高斯 迭代法、牛顿拉夫逊法以及 PQ解藕法等。介绍单电源辐射型网络和双端电源环形网络的 潮流估算方法。4-1潮流计算方程-节点功率方程1.支路潮流所谓潮流计算就是计算电力系统的功率在各个支路的分布、各个支路的功率损耗以及各个节点的电压和各个支路的电压损耗。 由于电力系统可以用等值电路来模拟， 从本质上说，电力系统的潮流计算首先是根据各个节点的注入功率求解电力系统各个节点的电压， 当各个节点的电压相量已知时，就很容易计算出各个支路的功率损耗和功率分布。假设支

3、路的两个节点分别为 k和丨，支路导纳为ykl，两个节点的电压已知，分别为V& 和V&，如图4-1所示。图4-1支路功率及其分布那么从节点k流向节点I的复功率为（变量上面的“一”表示复共扼）从节点I流向节点k的复功率为:功率损耗为:(4-3)Si 二 &i Sk =(V& -V&WM F) =-yki 因此，潮流计算的第一步是求解节点的电压和相位，根据电路理论，可以采用节点导 纳方程求解各个节点的电压。2.节点功率方程根据电路理论，要想求系统各个节点的电压，需要利用系统的节点导纳方程。I刃*t 1L ：f1测尸 V图4-2电网络示意图如图4-2所示的电网络，有 N个节点，假如已知各个节点的注入电

4、流源的电流，以及 各个支路的支路导纳，那么可以根据节点导纳方程求出电网各个节点的电压：YV =I s(4-4)其中YiiYl2LYiN -YY2221LY2NY =LLLLYniYn 2LYnn为电网络的节点导纳矩阵， Ykk( k=1,2,L N )为自导纳，是与 k节点所有连接支路导纳之和，Ykl( k =1 )为互导纳，等于负的连接 k和|节点的所有支路导纳之和。V=Vi,V2, L ,VnT为各个节点的电压相量，I s =lsi,ls2, L ,IsnT为注入到各个 节点的总电流。2.1节点复功率方程要想计算各个节点电压，除了需要知道系统参数及节点导纳矩阵以外，还需要知道节 点的注入电

5、流源的电流。然而电力系统中，节点的注入电流是不知道的，已知的是各个节 点的注入功率。这就需要将节点电压方程转化为节点功率方程。方程4-4中第k( k =1,2,L ,N )个节点的方程可以写作：N、YklV& 二 YdV& Yk2V&2 L YkkV&k L YrnV&n = l (4-5)&Sk在方程4-5两端乘以Vk，得到：(4-6)Vk YkiV&i Rk l&sk 三 Ssk = Psk 一 jQsk假如在电力系统中，各个节点的注入复功率都已知，那么就可以用方程 4-6组成的方程组求解各个节点的电压。然而实际情况并非如此，已知的条件是：有的节点的注入复功 率S是已知的，有的节点的电压幅

6、值和注入有功功率是已知的，有的节点的电压和相角是 已知的。根据这三种不同的情况，电力系统中各个节点分为三种类型： PQ节点、PV节点和V节点。所谓PQ节点，就是该节点的注入复功率 S是已知的，这样的节点一般为中间节点或 者是负荷节点。PV节点，指该节点已知的条件是注入节点的有功功率 P和该节点的电压幅值 V，这样的节点通常是发电机节点。V :.节点指的是该节点的电压幅值和相角是已知的，这样的节点通常是平衡节点，在 每个局部电网中只有一个这样的节点。当然，PQ节点和PV节点在一定条件下还可以互相转化，例如，当发电机节点无法维 持该节点电压时，发电机运行于功率极限时，发电机节点的有功和无功变成了已

7、知量，而 电压幅值则未知，此时，该节点由 PV节点转化为PQ节点。再比如某个负荷节点，运行要求电压不能越限，当该节点的电压幅值处于极限位置，或者电力系统调压要求该节点的电 压恒定，此时该负荷节点就由 PQ节点转化为PV节点。假如全系统有N个节点，其中有 M个PQ节点，N-M-1个PV节点，1个平衡节点， 每个节点有四个参数：电压幅值 V、相位角、：（用极坐标表示电压，如果用直角坐标表示电压相量则是e和f ）注入有功功率 Ps和无功功率Qs，任何一个节点的四个参数中总有 两个是已知的，因此 N个节点，有2N个未知变量，N个复数方程（即2N个实数方程， 实部和虚部各一个），通过解这个复数方程就可得

8、到另外 2N个参数。这就是潮流计算的本 质。但在实际求解过程中，由于我们求解的对象是电压，因此，实际上不需要 2N个功率方程，对于M个PQ节点，有2M个功率方程（M个实部有功功率方程， M个虚部无功功 率方程）；对于N-M-1个PV节点，由于电压有效值 V已知，因此只有 N-M-1个有功功率 方程；对于平衡节点，由于电压和相角已知，不需要功率方程。因此总计有 2M+N-M-仁N+M-1个功率方程。如果电压相量用极坐标表示，即 V& =VkE：k，则M个PQ节点有2M个未知数（M个电压有效值，M个电压相角），N-M-1个PV节点有N-M-1个未知数（电压有效值已知，未知数为电压相角），平衡节点没

9、有未知数，因此未知数的个数也是 N+M-1个，与方程数一致。如果复电压用直角坐标表示， V& =ek亠jfk，则有2（N-1）个未知数，还需要增加 N-M-1个电压方程，即 2 2。kk/feV+=2.2用直角坐标表示的电力系统节点功率方程对于PQ节点，已知的是注入节点的功率 p和Q,将Ykm =Gkm - jBkm和 V& =ek jfk带入节点功率方程的复数表示式中，可以得到有功功率和无功功率两个方程：N JN 1Psk 二 PGk - PLk 二 ek = (Gkmem - Bkm fm) fk * (Gkm fm BkmOm)(4-7)mA m 二- N _1 N _1Qsk 二 QG

10、k - QLk 二 fk (Gkmem - Bkm fm) - e (Gkm fm Bkmem)mA m 二上式中psk和QSk为注入到节点k的净功率，即注入和消耗的代数和。 PGk、QGk表示注入的功率，pLk和QLk为消耗的功率。对于PV节点，除了有功功率方程外，因为已知该节点的电压幅值，还有一个电压方程：2 2 2Vk ek f k （4-8）方程4-7可以抽象的表示为：APk（ei, fi,L 丄 f”）=0丿 （4-9）人Qk（e1,f1,L ,3 丄 fN）=0方程4-8可以抽象的表示为Vk（e1, f1,L ,eN，fN）=0 （4-10）因此，对于一个具有 N个节点的电力系统，

11、其中 M个PQ节点，N-M-1个PV节点, 1个平衡节点，有方程如下：也P（e1, f1,L ,eN，fN）=0.1-Q1 （e1, f1,L ,eN，f”）=0L L L L L L 2M个PQ节点的方程匚Pm G, 丄丄 f” v） = 0lQm G, f1,L , eN，fn）=0八Pm 1 (e1, f1,L ,eNv, fN4)=0八 Vm 1 (e1, f1,L , eN/, fN4)=0l Vn j (e1, f1,L , eN 4, fn j)N个节点，平衡节点的电压幅值和相角已知，即其横分量和纵分量已知，因此平衡节 点不参与计算。N-1个节点的电压的横分量和纵分量为未知数，共

12、 2N-2个未知数。2M个PQ节点方程，2（N-M-1）个PV节点方程，共计2N-2个方程。解这个方程组，就可以得到电力系统 N个节点的电压相量，根据各个节点的电压相量和已知的注入功率，就可以计算出各个支路的潮流分布，及各个支路的功率损耗。2.3极坐标表示的节点功率方程对于PQ节点，已知的是注入节点的功率 p和Q,将Ykm二Gkm jBkm和V&二Vk/l：k带入节点功率方程的复数表示式中，可以得到实部和虚部两个方程：N - .-(4-12)Psk 二 PGk - PLk 二 Vk= Vm (Gkm COS km Bkm Sin 门 km)Qsk = QGk - Q Lk =Vkg Vm (G

13、kmSin km Bkm COSkm)m -1上式中，V代表电压幅值， 亦二、* - 5。对于PV节点，由于节点的电压幅值已知，因此只有有功功率方程而没有无功功率方 程。同样，方程4-12可以抽象的表示为：1-Pk (V1 L ,Vm , 1,L , ；!n)=0 (4-13a)Qk(V1 L ,Vm ,、1,L ,、n)=0 (4-13b)lQm (V1,L ,Vm , -1 ,L，:n 4) 0因此，对于一个具有 N个节点的电力系统，其中 M个PQ节点，N-M-1个PV节点， 1个平衡节点，有方程如下：P(V1,L ,Vm , 1,L , : n j0 N _M -1个 PV 节点方程Q(

14、V1,L ,Vm , 1,L =0L L L L L L 眩M个PQ节点方程：Pm(V1,L ,Vm,、1,LJ匚 Pm 1 (V1 ,L ,Vm , -1 ,L，: n 4) = 0L L L L L L (4-14) ：Pn-1(V1,L ,Vm ,、1,L ,、nG =01除了平衡节点外，N-1个节点中，有 W彳!个QT功点的率方幅值和相角都是未知数， N-M-1个PV节点的相角为未知数，因此共有 2M+N-M-仁N+M-1 个未知数，2M+N-M-仁N+M-1个方程。在方程4-14中，可以把N-1个有功功率方程放在一起， M个无功功率方程放在一起：P(V1,L ,Vm , 1,L ,、

15、N4)= 0L L L L L LPn 4(V1 ,L ,Vm , 1,L，:Nd) =0|lQ(V1,L ,Vm , 1,L ,、nv)=0 Qm (V1,L ,Vm ,L L L L L L :M 个无功功率方程 61,L ,6N_) = 0U(4-15) 解上述方程组，就可以得到电力系统中各个节点的电压幅值和相角，进而可以计算出 各个支路的潮流分布和损耗。3.小结潮流计算是计算电力网各个支路的功率潮流分布和功率损耗，同时也计算各个支路的 电压损耗。首先要求电力网各个节点的电压相量。根据电网络理论，节点电压通常采用节 点导纳方程来求解，即已知电网络的节点导纳矩阵和各个节点的注入电流源的电流

16、，求解 节点导纳方程。然而通常电力系统各个节点的注入电流是未知的，已知的是各个节点的注 入功率，因此需要将节点电压方程转化为节点功率方程。实际电力系统的节点注入功率并非都已知，有的已知注入有功功率 P和无功功率Q称为PQ节点；有的已知注入有功功率 P和节点电压有效值 V，称为PV节点；有的已知节点电压V和相角d,称为平衡节点或 V、节点。无论哪种类型节点，每一个节点均含有 4个参量P、Q、V、(或e、f)已知的是其中的两个，故而可以利用节点功率方程( 4-6)求解出另外两个参量。假设系统有 N个节点，必然有2N个未知数，同样有2N个节点功率方程(4-17中的实部和虚部各一个)。实际上，我们求解

17、的目标是电压，对于 PV节点和V、节点来说，前者电压有效值已知，后者电压相量已知，因此不存在 2N个未知数，当然也不需要 2N个方程。假设系统有 N个节点，M个PQ节点，1个平衡节点，对于直角坐标表示的节点电压来说，有 2(N-1)个未知数，2M+N-M-1个功率方程，只需要再补充 N-M-1个电压方程就可以了；对于极坐标 表示的电压来说，只有 N-1个未知数，M个V的未知数，因此只需要 N+M-1个功率方程 就足够了。无论怎样，潮流计算是解决这样的一组非线性代数方程组：F(X, C, U ) = 0 (4-16)其中，X代表系统状态，包括电压 V和相角J.； C代表参数，包括电导 G和电纳

18、B ； U表示系统激励，即注入的功率。求解这样的多维非线性代数方程组，需要利用计算机进行辅助迭代计算，即先给定一 个初值，然后不断迭代，逼近真实解。方法有：高斯 -赛德尔迭代法，牛顿-拉夫逊法和PQ解耦法。4-2高斯赛德尔叠代法1.基本原理为了方便理解这个n维方程组的叠代求解方法，先从一元非线性方程的求解开始。假 设有一维方程f(x) =0，高斯法的基本原理是，先将方程转化为：x=g(x) (4-17)那么给定一个初值x0，代入就可以得到一个新值 x二g(x0)，第k次叠代的值为:(4-18)(4-19)k 1 / 叫x g(x )一直叠代到误差满足要求为止，即N N -4x -X这个解方程的

19、方法称为高斯叠代法。这个叠代求解的过程可以这样来理解：x = g(x)的解可以认为是两个曲线 y = x和讨二g(x)的交点的横坐标x，首先给定一个初值x0， g(x0)与斜线y = x的交点的横坐标即为叠代后的新解 x1， g(x)与斜线y =x的交点的横坐标即为叠代后的新解 x2，如此围绕交点往复循环，不断地逼近方程的解，如图4-4所示。n为方程组为:(4-20)(4-21)高斯迭代法可以推广到 n维非线性代数方程组，假设 fl (xi,x2 L , xn) =0f2 (xi, x2 L , Xn ) = 0L L L L L Lfn(Xi,x2L ,Xn) =0首先将方程组4-20转化为

20、:Xi =g(Xi,x2,L ,xn)刈=g(Xi,x2,L ,Xn)L L L L L LM =g(Xi,X2 ,L ,Xn)给定一组初始值X0=，丄，Xn0T，带入上式，得到一组新值 x1 = g(X0)，不断叠代，循环往复，第 k次叠代为:(4-22)Xk 1二 g(Xk)其中第j个方程为k 1 / k k kx gj(xi ,X2 ,L ,Xn ) (4-23) 直到叠代前后的解的最大误差不超过允许的误差为止，即N 1 N| , 、maxxj -Xj (4-24)为了提高高斯叠代法的收敛速度，赛德尔提出将已经叠代出的新值代替旧值参与叠代 计算，如在第k次叠代中，第j个方程为k 1 /

21、k 1 1 k k、X = gj(Xi ,L ,x _ ,X ,L ,Xn ) (4-25) k 第1至j-1个元素已经叠代出k+1次的值，因此代替第 k次的值参与第j个元素的叠代, 就可以提高收敛速度。2.电力系统潮流计算的高斯-赛德尔迭代法电力系统潮流计算需要求解节点功率方程，其中第 m (m=1,2,N个节点功率方程为：N p 2 N q(4-26)(4-27)Vm YmlV &1 二 YmmVm fm、YmV & I = P Sm jQsm 14 l =1如上式变换为x=g(x)的形式，可以得到如下的方程：/& _ 1 ( PSm jQSm _、N YmlV&l)m Ymm Vm yl

22、 rm根据高斯-赛德尔迭代法，首先给定电压相量的初值，对于 PQ节点，不仅需要给定电压幅值的初值，还要给出相角的初值(设为零)。假如第m号节点为PQ节点，第k次叠代公式为(第 m个节点以前的节点第 k次叠 代已经完毕，因此用 k+1次的值取代k次的值，而在第 m个节点以后的节点尚未进行第 k次叠代)：对于PV节点，给定的初值的电压幅值为给定的电压，相角初值设为零。可是对于 PV节点来说，注入该节点的无功功率未知，因此第 k次叠代时，首先按照下式计算注入 PV节点(假设第m个节点是PV节点)的无功功率：k & kl k 1 & 国 N【k_1Qsm lmV m I Sm = lmV m ( Ym

23、lVl YmlVl (4-29)k ) l =1 lorn如果在叠代计算过程中，任意节点的电压和无功功率必须满足不等约束条件：如果在叠代过程中，PQ节点的电压幅值超出允许的范围，则该节点的电压幅值就固定为允许电压的上限(如果超出上限)或下限(如果越过下限) ，PQ节点就变为PV节点继续进行叠代。同样，对于 PV节点来说，如果在叠代过程中，无功功率 Q超出了允许的范围，则PV节点就变为PQ节点继续参与叠代。高斯-赛德尔叠代法的计算过程如下：(1)第一步：设置初始值，对于 PQ节点，由于其电压相量的幅值和相角都未知，因此初始的电压相量的幅值可以设定为各个点的额定电压，相角选择为零；对于 PV节点，

24、由于其电压相量的幅值已知，因此幅值用已知的设定电压，初始相角设定为零。(2) 第二步：对于 PQ节点，直接将设定的初始值代入，用 4-28求得下一次迭代的 电压值，然后判断是否电压越限，如果越限，则用其限值(越过上限用上限值，越过下限 则用下限值)，该节点在下一次迭代过程中转化为 PV节点；对于PV节点，则首先利用式 4-29求出注入的无功功率，然后校验无功功率是否越限，如果越限则采用上限值或者下限值，下一次迭代时该节点转化为 PQ节点，将求得的注入无功功率和已知的有功功率代入4-28求解下一次迭代的电压相量值。(3) 第三步：判断误差是否满足要求，用第 k次迭代的结果和k-1次迭代的结果进行

25、 比较，如果其最大的误差满足事先设定的误差要求，则输出计算结果，如果不满足要求，则返回第二步继续迭代。其计算流程图如图 4-5所示。没宜初值訂PQ召点？卜理一计算无功Q(k)*Q越限？fQ駆限值，菸为PQ节点计算电叭 gv.k)i取限值.转为PVP点亠T赵遍计算电图4-5高斯赛德尔迭代法求解电力系统潮流的计算流程图4-3 牛顿拉夫逊法因此，可以得到这个一元非线性方程的求解步骤为：首先给定解的初值 x0，然后根据公式4-32求出初始值的修正值.-；x0，由此可以得到该方程的新的解 X1 = X0 AX， 如此反复叠代，直到误差满足要求 Ax【N| g。迭代计算流程如图4-6所示。楠出结END图4

26、-6牛顿拉夫逊法计算流程其迭代求解过程的几何意义如图 4-7所示。图4-7牛顿拉夫逊法的几何解释可以把上述求解一元非线性代数方程的方法推广到 n维非线性代数方程(如 4-20)的求解。非线性代数方程组 4-20可以表示为矩阵形式：F(X) =0同样假定X0是该方程组的近似解，与真实解之间的误差为勒级数：(4-33)X，在X0处展开一阶泰F(X0 ：X) ： F (X0) J X = 0(4-34)其中:LcXCX21Cf2L依1CX2LLLCfnL01CX2:flXncf2-XnLCfnXn _x_x 0(4-35)X J F (X 0)(4-36)被称为雅克比矩阵。4-34称为修正方程，修正

27、方程可得到修正值 X :计算过程与一维方程的牛顿法求解类似，首先给定初值 X 0 = X10，L , Xn0 T，并计算出在初始值处的雅克比矩阵 J0，利用4-36式计算初始值的修正值X0二J ：F(X0)，根据这个差值可以得到修正后的解 X二X0 厶X0。如此循环往复，在第k次叠代时，计算雅克比矩阵 Jk，根据4-34计算修正值 X k = -J ：F(Xk)，得到第k+1次修正后的解：x k1 =xk 厶xk，重复上述过程，直到误差满足要求为 止。可见，牛顿拉夫逊法的关键在于求解雅克比矩阵 J,由于直角坐标表示和极坐标表示电压相量的节点功率方程有所不同，因此其雅克比矩阵也有很大的差异。(4

28、-37)2.直角坐标节点功率方程的牛顿-拉夫逊法仍然假设系统有 N个节点，其中 M个PQ节点，N-M-1个PV节点，1个平衡节点。 则M个PQ节点方程为(假设1号节点至M号节点为PQ节点)：NN-PSk - ek V (Gkl el - Bki fi ) - fk Xk -1,2,L ,M。(4-38) 72 2 2-Vk - Vk -e -fk 011 _11LmLMmLm |1Mm1Lmn1LMn1II-1佔LNrrtm1Hkfm |Hkltm1LNtn!II% Rh11SnLRmmSrr|m |IRr1m1Stpi!LFTn!nl fr1LLLL|LLLLL11H11LNmHUm |N1

29、m1H1m1LNt!1 111SnnLFR1mS1m |S1drLFR4n11 1(4-39)其中：Nkj=：Pk - -Gkjek - Bkj fk ( j - k):：Pk - -Gkkek - Bkk f k - (Gkl el - Bkl fl ) ,:ek i -1nA=Bkkek Gkk fk -、 (Gki fi 亠 Bki e) l丄n=Bkkek -Gkk fk +迟(Gkl fi + Bkie)i=1Lkj；:：Qkf=Bkj fk GkjekLkk Ql =Gkkek Bkk fkf kn .1(Gki eii d-Bki fi)Rkj-k)_ej基于直角坐标的牛顿-拉夫逊法求解潮流计算的步骤如下：(1) 第一步：设定初值，对于 PQ节点，其电压幅值的初值设定为