电力系统分析基础(第四章).pptx

1、电力工程系,Department of Electrical Engineering,电力系统分析基础 Power System Analysis Basis （四） 任 建 文 North China Electric Power University,TCSC、STATCOM、 UPFC、 FACTS,第四章 复杂电力系统潮流的计算机算法 复习上一章内容： 1. 潮流计算的目的及内容 2. 电压降落（损耗、偏移）、功率损耗的计算 3. 手算潮流的原理和方法,1)辐射型： 2)两端型： 3)环网型：,同一电压等级：已知末端电压或首末端电压 不同电压等级：归算电压或折算参数 计算自然功率（力矩

2、原理）、强迫功率 找功率分点、打开、按辐射型计算 单级：从电源点打开无强迫功率 多级：电磁环网归算法、等值法,4. 潮流调整：自然分布、串联电容、串联电抗、附加串联加压器,第四章 复杂电力系统潮流的计算机算法,本章主要内容：,1. 建立数学模型：节点电压方程、导纳矩阵的形成与修改 2. 功率方程、节点分类及约束条件 3. 迭代法计算潮流,功率方程的非线性性质 高斯塞德尔法,用于潮流计算速度慢、易于收敛,4. 牛顿拉夫逊法计算潮流,原理：局部线性化,直角座标法、极座标法、PQ分解法,用于潮流计算速度快、但注意初值选择,.,B,第一节 电力网的数学模型 一、节点电压方程,. IB = Y B U

3、二、导纳矩阵的形成 三、导纳矩阵的修改,YB节点导纳矩阵,&,&,&,第一节 电力网的数学模型,1、节点电压方程（示例）,U1,y12,U 2,. E 1, . U 1,. U, . U 2 3,. E 2,. I 1,y10,y13 y23 U 3 y30, y20,. I 2, 参考节点的选取接地点,Z12,. . . . . . I 1 = U 1 y 10 + ( U 1 U 2 )y 12 + ( U 1 U 3 )y 13,Z1 . E 1,Z13,Z3,Z23,Z2 . E 2,. . . . . . I 2 = U 2 y 20 + ( U 2 U 1 )y 21 + ( U

4、2 U 3 )y 23 . . . . . 0 = U 3 y 30 + ( U 3 U 1 )y 31 + ( U 3 U 2 )y 32,自,导,纳,互,导,纳,第一节 电力网的数学模型 . . . . . . . I 1 = (y10 + y 20 + y 30 ) U 1 y12 U 2 y13 U 3 = Y11 U 1 + Y12 U 2 + Y13 U 3 . . . . . . . I 2 = y 21 U 1 + (y 20 + y 21 + y 23 ) U 2 y 23 U 3 = Y21 U 1 + Y22 U 2 + Y23 U 3 . . . . . . 0 = y

5、 31 U 1 y 32 U 2 + (y 30 + y 31 + y 32 ) U 3 = Y31 U 1 + Y32 U 2 + Y33 U 3,Y 11 = y10 + y12 + y13 Y 22 = y20 + y21 + y23 Y 33 = y30 + y31 + y32,. I 1,y10,y12 y13 ,y23, y20,. I 2,Y 12 = Y 21 = y 12 Y 23 = Y 32 = y 23 Y 12 = Y 31 = y 13,y30 注：Y距阵的维数（n-1）, I. 1 Y 11, . , Y 21, ,M , .,Y 1 n U 1 , . , .

6、, ,节点i: 加单位电压 U = 1i,j,i, I Y ii = y i 0 + y,ij,ii,Y = . 节点 i 注入网络电流,j = 1, U,2、对称性,j i,i,i, I ,j,j,ij, U (U = 0, j i),自,导,纳,互,导,纳,1、阶数,节点i:,第一节 电力网的数学模型 2、导纳矩阵的形成, I 2 = K I n Y n 1,Y 12 Y 22 K Y n 2,K K K K,Y 2 n U 2 K M Y nn U n ,. . I = Y U 节点导纳距阵的特点：, . n 其余节点j: 全部接地 U = 0 ( U j = 0 , j i ) Yii

7、0 . 3、稀疏性加单位电压 U = 1 Y = . Yij = Yji = yij 其余节点j: 全部接地 U = 0 由地流向节点j的电流 i j 稀疏性：当yij=0 时Yij=0,yij,i,j,i,j,第一节 电力网的数学模型 3、导纳矩阵的修改 1) 增加一节点 节点导纳矩阵增加一阶, Yii = yij,Yjj = yij,Yij = Yji = - yij 2) 增加一条支路,导纳矩阵的阶数不变 Yii = Yjj = yij Yij = Yji = - yij,yij,i,j,i,j,第一节 电力网的数学模型 3) 切除一条支路 相当于增加一导纳为（ -yij ）的支路,导纳

8、矩阵的阶数不变 Yii = Yjj = - yij Yij = Yji = yij 4) 修改一条支路的导纳值（ yij 改变为yij ） 导纳矩阵的阶数不变 Yii = Yjj = yij - yij Yij = Yji = yij - yij ,- yij -yij yij ,第一节 电力网的数学模型 5) 修改一条支路的变压器变比值（ k*改变为k* ） yT / k*,i yT(k*-1) / k*,j yT(1- k*) / k*2, Yii = 0 Yjj =(1/ k*2 - 1/ k*2 ) yT, Yij =,Yji =-(1/ k* - 1/ k*) yT,.,& = S

9、,I = Y,U I U ,&, ,.,n,= Yij U j,( Pi + jQi )*,j = 1,n,n,j =1 j =1,n,第二节 功率方程、节点分类及约束条件 . * 一、功率方程 * (i = 1,2,L n) (4 35) U i Pi = ei (Gij e j Bij f j ) + fi (G ij f j + Bij e j ) (4 36a) n n Qi = fi (Gij e j Bij f j ) ei (G ij f j + Bij e j ) j =1 j =1,Pi = U i U j (Gij cos ij + Bij sin ij ) j =1 n

10、Qi = U i U j (Gij sin ij Bij cos ij ) j =1,(4 36b),第二节 功率方程、节点分类及约束条件,二、节点分类, 一个电力系统有n个节点，每个节点可能有4个变量Pi,Qi ,ei,fi或Pi,Qi ,Ui, i,，则共有4n个变量，而上述功率方程只有,2n个，所以需要事先给定2n个变量的值。根据各个节点的已 知量的不同，将节点分成三类：PQ节点、PV 节点、平衡节 点。, 1、PQ节点(Load Buses), 已知Pi,Qi ，求,ei, fi（ Ui, i, ），负荷节点（或发固,定功率的发电机节点），数量最多。 2、PU节点(Voltage Co

11、ntrol Buses), 已知Pi, Ui ，求, Qi, i, ，对电压有严格要求的节点，,如电压中枢点。,第二节 功率方程、节点分类及约束条件,二、节点分类,3、平衡节点,（Slack Bus or Voltage Reference bus）,已知Ui ， i,，求, Pi, Qi, ，只设一个。,设置平衡节点的目的,在结果未出来之前，网损是未知的，至少需 要一个节点的功率不能给定，用来平衡全网 功率。,电压计算需要参考节点。,第二节 功率方程、节点分类及约束条件,三、约束条件, 实际电力系统运行要求：, 电能质量约束条件：Uimin Ui Uimax, 电压相角约束条件 |ij|=|

12、 i - j | ijmax,稳定运行的一个重要条件。, 有功、无功约束条件 Pimin Pi Pimax,Qimin Qi Qimax,&,& & &,i i,n,n,j =1 j =1,直角坐,n,极坐标形式：,第三节 高斯塞德尔迭代法潮流计算 一、功率方程的非线性, * I = YU and I = S U ,(P + jQ)* * Ui,n . = Yij U j j=1,(i =1,2,Ln),(4 35),Pi = ei (Gij e j Bij f j ) + fi (G ij f j + Bij e j ) (4 36a ) n n 标形式： Qi = fi (Gij e j

13、Bij f j ) ei (G ij f j + Bij e j ) j =1 j =1 Pi = U i U j (Gij cos ij + Bij sin ij ) j =1 n (4 36b) Qi = U i U j (Gij sin ij Bij cos ij ) j =1 非线性方程组，不能用常规代数求解方程方法求解 。,第三节 高斯塞德尔迭代法潮流计算,两种常见的求解非线性方程的方法：, 高斯-塞德尔迭代法, 牛顿-拉夫逊迭代法,= 0 = 0.6667,(1),2,x,( 2 ),(,= ,x2 解：（1）将方程组 ( k + 1 ) 2 ( k ) ( k ),x1 = 0.

14、5817,x23) = 0.8167, 2,1,2,1,3,= 0.3333x ( ) = 0 + 3,x = 4815,( (,第三节 高斯塞德尔迭代法潮流计算,二、高斯-塞德尔迭代法原理及求解步骤,例6-11 已知方程组 3x1 + 2 x1 x2 1 = 0 2 3 3x2 x1 x2 + 2 = 0,用高斯-塞德尔求解（0.01）。,( 3) x 1 x 1 x 2 +,3,( 改写成迭代公式：, x ( k + 1 ) = 1 x ( k ) x ( k ) 2,直到|x(k+1)-x(k)| 3 3,（2）设初值 x1 0 ) = x20 ) = 0；代入上述迭代公式,第三节 高斯

15、塞德尔迭代法潮流计算 二、高斯-塞德尔迭代法原理及求解步骤, 设有非线性方程组 的一般形式： 将其改写成下述便 于迭代的形式：,f1 ( x1 , x2 , x3 ,L xn ) = 0 f2 ( x1 , x2 , x3 ,L xn ) = 0 LLLLLLLL fn ( x1 , x2 , x3 ,L xn ) = 0 x1 = g1 ( x1 , x2 , x3 ,L xn ) x2 = g2 ( x1 , x2 , x3 ,L xn ) LLLLLLLL xn = gn ( x1 , x2 , x3 ,L xn ),(0),(0),第三节 高斯塞德尔迭代法潮流计算 二、高斯-塞德尔迭代

16、法原理及求解步骤 假设变量（x1, x2, .,xn）的一组初值,（,x1(0) , x2 ,L, xn,）, 将初值代入迭代格式（6-18），完成第一次迭代 将第一次迭代的结果作为初值，代入迭代公式，进行第二 次迭代 检查是否满足收敛条件： | xi( k +1) xi( k ) |max ,(1),(0),(0),(0),(0),(1),(1) (0) (0),(0),(0) ,(1),(1) (1) (1),(1) (0),xn = gn(x , x2 , x3 Lxn1, xn ) ,1 1,1,第三节 高斯塞德尔迭代法潮流计算 二、高斯-塞德尔迭代法原理及求解步骤, 迭代公式：,x

17、= g1(x , x2 , x3 Lxn ) x2 = g2(x , x2 , x3 Lxn ) LLLLLLLLLLLLLLL xi = gi (x1 , x2 , x3 Lxi1, xi ,Lxn ) LLLLLLLLLLLLLLL (1) (1) (1) (1) (1) (0) ,(,(,(, 更,( k + 1) ( k + 1),( k ),( k ),( k ),( k ) ,( k ),( k +1) ( k +1) ( k +1) ( k +1) ( k +1),( k ),x,(k+1) (k+1) (k+1),(k),(k),i , ,ii i,第三节 高斯塞德尔迭代法潮流

18、计算 二、高斯-塞德尔迭代法原理及求解步骤 x1( k +1) = g1 ( x1( k ) , x2k ) , x3k ) L xnk ) ) 一 x2 = g 2 ( x1 , x2 , x3 L xn ) 般 LLLLLLLLLLLLLLL 的 ( k +1) ( k +1) ( k +1) ( k +1) ( k +1) 形 xi = gi ( x1 , x2 , x3 L xi 1 , xi ,L xn ) 式： LLLLLLLLLLLLLLL xn = g n ( x1 , x2 , x3 L xn 1 , xn ) , 简化形式：,(k+1) i,=g (x ,x2 ,Lx1 ,

19、x ,Lxn ),=,1 i n,第三节 高斯塞德尔迭代法潮流计算 二、高斯-塞德尔迭代法原理及求解步骤 迭代收敛条件： | xi( k +1) xi( k ) | max .(i = 1,2,L, n) 同一道题可能存在多种迭代格式，有的迭代格式收敛， 有的迭代式不收敛。下面讨论收敛条件：,当迭代格式为 定理,n xi( k +1) = bij x (j k ) + g i j =1,i = 1, 2,L , n,n 如果 L max | b ij | 1 n j = 1 则迭代格式 xi = bij xi + gi j=1 对任意给定的初值都收敛。,i = 1,2,L,n,.,n,= Yi

20、j U j,(Pi + jQ i )*,j=1, j =n ,1 P jQ,.,.,Yii U,i,第三节 高斯塞德尔迭代法潮流计算 三、高斯-塞德尔迭代法潮流计算 1. 方程表示： 用高斯-塞德尔法计算电力系统潮流首先要将功率方程 改写成能收敛的迭代形式 Q : 设系统有n个节点，其中 m个是PQ节点，n-(m+1) 个是PV节点，一个平衡节点，且假设节点1为平衡节点,（电压参考节点） 功率方程改写成： * U i U i = i * i Yij U j j =1 j i ,(i = 1,2,Ln), ,* ( k +1),* ( k ),* ( k ),* ( k ),*,1 P 2 jQ

21、 2 ,Y22 ,(,* ( k +1) * ( k +1) ,* ( k ),* ( k ),*,1 P3 jQ3,Y33 ,* ( k ),1 Pn jQn,*, ,Ynn ,(,第三节 高斯塞德尔迭代法潮流计算,三、高斯-塞德尔迭代法潮流计算,或更具体的形式为：,U 2 = * Y21 U1 Y23 U 3 Y24 U 4 L Y2 n U n , U 2k ) ,U 3 = * Y31 U1 Y 32 U 2 Y24 U 4 L Y2 n U n , U 3 k ) , LLLLLLLLLLLLLLLLLLLLLLLLLLL ,LLLLLLLLLLLLLLLLLLLLLLLLLLL ,

22、* ( k +1) * ( k +1) * ( k +1) ,U n = * Yn1 U1 Y n 2 U 2 Y n 3 U 3 Y24 U 4 L ,( k ), U n , ,第三节 高斯塞德尔迭代法潮流计算 三、高斯-塞德尔迭代法潮流计算 2. 求解的步骤： 上述迭代公式假设n 个节点全部为PQ节点。 式中等号右边采用第k次迭代结果，当ji时，采用第k次迭代结果。 用G-S迭代法求解的步骤： 第一步：形成节点导纳距阵； 第二步：设除平衡节点外的其它节点的初值，一般都 * ( 0 ) 设 U i = 10 (i = 2,3,L n) ； 第三步：迭代求解，判断收敛与否？若满足收敛条件，,

23、则迭代停止,* * ( k +1) * ( k ) U i max =| U i U i | ,第三节 高斯塞德尔迭代法潮流计算 三、高斯-塞德尔迭代法潮流计算 3. PV节点的处理： 由于该类节点的V已知，Q未知，故在给定初值时，对 该类节点增加初值 Q i( 0 ) = 0 . 5 Pi ； 增加计算无功的迭代公式：,Qi( k ), . ( k ) i 1 * * ( k +1) n * * ( k ) = ImUi ( Yij U j + Yij U j ) j=1 j=i ,(6 25), 对于PV节点的计算步骤： 除了完成（6-24）的迭代计算外，还要执行（6-25） 的迭代计算

24、对（6-25）得到的结果要进行下列三种情况的校核：,( k ),( k ),第三节 高斯塞德尔迭代法潮流计算,三、高斯-塞德尔迭代法潮流计算, (a) Qi( k ) Qi min 这种情况由于计算得到的结果比允许的最小,值还小，所以不允许以计算得到的结果再代入进行迭代，以,Qi min 作为PV节点的无功功率，此时，PV节点就转化为PQ节点 (b) Qi Qi max 这种情况由于计算得到的结果比允许的最大,值还大，所以不允许以计算得到的结果再代入进行迭代，而 是以 Q i max 作为PV节点的无功功率，此时，PV节点就转化 为PQ节点。, (c) Qi min Qi Qi max 因求出

25、的无功功率满足要求，所以迭,代得到的结果继续代入公式（6-25）进行计算,( k +1),),( k +1),第三节 高斯塞德尔迭代法潮流计算 三、高斯-塞德尔迭代法潮流计算 对于PV节点，由于它的U值是给定的，每次用公式（6-24）,得到的结果,* ( k +1) U i, (U,( k +1) i, ,( k +1) i,中的 U i,一般不,等于给定的值，这种情况要用给定的U代替计算得到的幅,值，用 U i i,组成新的电压初值。如果通过迭代,得到的与限值比较已经越限，则转化为PQ节点后，就不必 做电压幅值的更换了。,n,*,*,*,*, ,*,*,*,*,*,*,ij,第三节 高斯塞德

26、尔迭代法潮流计算 三、高斯-塞德尔迭代法潮流计算 4. 潮流计算： 平衡节点的功率：,支路功率：, * * * S1 = U1 Y1 j U j = P1 + jQi j=1, * * S ij = U i Iij = U i Ui yi0 + (Ui U j ) yij = P + jQij S iji = U j I ji = U j U j y j 0 + (U j Ui ) y ji = Pji + jQji 支路功率损耗： S ij = S ij + S ji = Pij + jQij,(6 27) (6 29),第三节 高斯塞德尔迭代法潮流计算,四、高斯-塞德尔迭代法潮流计算流程图

27、,见书上P155,1,1.17-j4.71,3, ,第三节 高斯塞德尔迭代法潮流计算,四、例题：用G-S计算潮流分布 平衡节点 U1=1.00,y13,PU节点 P3=0.4, U3=1.1 ,y12,5.88-j23.5 y30 j0.33,PQ节点 S2=-0.8-j0.6 2 解：网络的节点导纳距阵为： Y11 Y12 Y13 7.05 j 28.21 5.88 + j 23.5 1.17 + j 4.71 YB = Y21 Y22 Y23 = 5.88 + j 23.5 5.88 j 23.5 0 Y31 Y32 Y33 1.17 + j 4.71 0 1.17 j 4.38 ,* (

28、 0 ),U 2,=,=,* ( 0 ),Y33 ,=,=,第三节 高斯塞德尔迭代法潮流计算,设, ( 0 ) U 2, ( 0 ) = 1 .0 0 , U 3,= 1 .1 0 ,Q,( 0 ) 3,= 0 .20,，代入式（6-24）求, (1) U 2, (1) U 2,1 P2 jQ 2 Y 22 , ( 0 ) Y 21 U 1 Y 23 U 3 ,1 0 .8 + j 0 .6 5 .88 j 23 .5 1 .0 0 , ( 5 .88 + j 23 .5 ) 1 .0 0 0 (1 .1 0 ) ,= 0 .9680 j 0 .0260 = 0 .9683 1 .539 ,

29、(1 ) U 3,1 P3 jQ 3( 0 ) U 3, (1 ) Y31 U 1 Y32 U 2 ,1 0 .4 j 0 .2 1 .17 j 4 .38 1 .1 0 , ( 1 .17 + j 4 .71) 1 .0 0 0 ( 0 .9683 1 .539 ) ,= 1 .1298 + j 0 .0484 = 1 .1310 2 .451 ,* (1),*,Q31),(,= Im = 0.0685,* (1 ),Y 22 ,=,=,第三节 高斯塞德尔迭代法潮流计算,修正U3为, (1) U 3 = U 3 3(1) = 1.12.451,，再用式（6-25）计算：, (1) * (1)

30、 = ImU 3 ( Y 33 U 3 + Y 31 U 1 + Y 32 U 2 ) 1.1 2.451 (1.17 j4.38) 1.12.451 + 1.1 2.451 (1.17 + j4.71) 1.00 + 1.1 2.451 0 0.9683 1.539 然后开始第二次迭代：, ( 2 ) U 2,1 P2 jQ 2 U 2, (1 ) Y 21 U 1 Y 23 U 3 ,1 0 . 8 + j 0 . 6 5 . 88 j 23 . 5 0 . 9683 1 . 539 , ( 5 . 88 + j 23 . 5 ) 1 . 0 0 0 (1 . 1 2 . 451 ) ,=

31、 0 . 9662 j 0 . 0260 = 0 . 9665 1 . 541 ,Y33,=,=,* (2),*,(2),= Im = 0.0596,再,计,算, ) , ),第三节 高斯塞德尔迭代法潮流计算, ( 2 ) U 3,1 P3 jQ 3(1 ) * (1 ) U 3, ( 2 ) Y 31 U 1 Y32 U 2 ,1 0 .4 j 0 .0685 1 .17 j 4 . 38 1 . 1 2 . 451 , ( 1 . 17 + j 4 . 71 ) 1 . 0 0 0 ( 0 .9665 1 .541 ) ,= 1 .1011 + j 0 .0566 = 1 .1026 2

32、.940 ,( 2) 再修正U3为：U 3,= U 3 3( 2) = 1.12.940, (2) * (2) Q3 = ImU3 (Y33 U3 +Y31U1 +Y32 U 2 ) 1.1 2.940 (1.17 j4.38 1.12.940 +1.1 2.940 (1.17+ j4.71 1.00 +1.1 2.940 00.96651.541 ,因此，第二次迭代结束时节点2的电压为, U 2 = 0 . 9662 j0 . 0260 = 0 . 9665 1 . 541 ,节点3的电压相位角为3=2.940,与之对应的节点3的无功功率为Q3=0.0596.,第四节 牛顿拉夫逊法潮流计算

33、一、N-R原理 1. 非线性方程的求解： f(x)=0 设：x(0)为的初始近似解， x(0)为与真实解的偏差 则：x= x(0) x(0) f(x(0) x(0)=0,按Taylors展开,f(x(0) x(0)=,f(x(0)- f(x(0) x(0)+.+(-1)n fn(x(0) ( x(0)n/n!+.=0,第四节 牛顿拉夫逊法潮流计算 由于 x(0）较小，故忽略高次项后：,f(x(0),x(0)= f(x(0)- f(x(0) x(0) =0,x(0) = f(x(0)/ f(x(0),x(1) = x(0) -,x(0) = x(0) - f(x(0)/ f(x(0),k次迭代时

34、修正方程为： f(x(k)- f(x(k) x(k) =0 x(k) = f(x(k)/ f(x(k) x(k+1) = x(k) - f(x(k)/ f(x(k),第四节 牛顿拉夫逊法潮流计算 结束迭代的条件（收敛）：|f(x(k)|1 或 | x(k)|2,物理意义,初值不当不收敛,第四节 牛顿拉夫逊法潮流计算 2. 非线性方程组的求解：,推广于（6-16）表示的 多变量非线性方程组,f1 ( x1 , x 2 , x 3 ,L x n ) = 0 f 2 ( x 1 , x 2 , x 3 , L x n ) = 0 LLLLLLLL f n ( x1 , x 2 , x 3 ,L x

35、n ) = 0 ,(4 30),(4-31),| 0,第四节 牛顿拉夫逊法潮流计算,(4-32),fi,式中 x j 为函数 fi(x1,x2,.xn) 对自变量 xj 的偏导在初,始值处的值,用,矩,阵,表,示,第四节 牛顿拉夫逊法潮流计算,(4-33),得到新的近似解：,更,一,般,的,表,示,第四节 牛顿拉夫逊法潮流计算 (4-33a),第k+1次迭代后的解为：,(4-33b), f i,=,jij, x j,第四节 牛顿拉夫逊法潮流计算,式(4-33a)可简写为：,F(x(k)=J(k),x(k),(4-34),J(k) 为n n阶雅可比矩阵,其元素 ( k ) | x ( k ) 为

36、函数 fi(x1,x2,.xn) 对自变量 xj 的偏导在点（ x(k) ）的值,式(4-33b)可简写为：,x(k+1)= x(k) x(k),(4-34b),第k次迭代后用下面的公式检查是否收敛 (4-35a) (4-35b),第四节 牛顿拉夫逊法潮流计算 二、潮流计算时的修正方程式（直角座标法） 1. 直角座标法： PQ,节 点 PU 节 点,(4-38a) (4-38b), u 2 R S ,第四节 牛顿拉夫逊法潮流计算,n1 直角坐标的缩写形式： m1 nm, P H N f Q = J L e ,Hij,=, Pi f j,J,ij,=, Qi f j,Rij,=,2 Ui f j

37、,Nij,=, Pi ej,Lij,=, Qi ej,Sij,=,2 Ui ej,第四节 牛顿拉夫逊法潮流计算,直角座标法矩阵表示,(4-37),第四节 牛顿拉夫逊法潮流计算,雅可比矩阵元素值,对角元素(i=j),非对角元素(ij),(4-41a) (4-41b),第四节 牛顿拉夫逊法潮流计算 矩阵的特点及计算步骤 雅可比矩阵的特点： 各元素是各节点电压的函数 不是对称矩阵 Yij=0， Hij= Nij= Jij= Lij= 0，另Rij= Sij= 0，故稀疏 牛顿拉夫逊法潮流计算的基本步骤： 1. 输入原始数据和信息：y、C、Pis、Qis、Uis、约束条件 2. 形成节点导纳矩阵YB=

38、CTyC 3. 设置各节点电压初值ei(0), fi(0),4. 将初始值代入(4-38)求不平衡量 Pi(0),Qi(0),Ui2(0),fi(k),s sj j s s,第四节 牛顿拉夫逊法潮流计算 5. 计算雅可比矩阵各元素（Hij、Lij、Nij、Jij、Rij、Sij）,6. 解修正方程(4-37) ，求,ei(k),fi(k),7. 求节点电压新值ei(k+1) =ei(k) -,ei(k), fi(k+1) = fi(k) -,8. 判断是否收敛：Max|,fi(k) |, Max|,ei(k) |,9. 重复迭代第4、5、6、7步，直到满足第8步的条件 10. 求平衡节点的功率

39、和PV节点的Qi及各支路的功率 . . n * * Ss = U Y U = P + jQ n j=1 Qi = fi(Gijej Bijfj) ei(Gijfj + Bijej) j=1,潮流计算时 的修正方程 式（框图）,见书上P165,i i, , , ,Pi,ij,i,U,ij,j,ij,i,U,ij,j,i,第四节 牛顿拉夫逊法潮流计算 2. 极座标法：,n P = P Ui U j (Gij cosij + Bij sinij ) = 0 j = 1 n Qi = Qi Ui Uj (Gij sinij Bij cosij ) = 0 j=1,(4 45a) (4 45b),n1

40、P H m1 Q = J,N L U / U ,H = j,P Q N = j U J = j,Q L = j U,第四节 牛顿拉夫逊法潮流计算,极座标法矩阵表示,(4-44),展,开,式,计,及,第四节 牛顿拉夫逊法潮流计算,极座标法系数推导,(4-45a),(4-45b),(4-48),第四节 牛顿拉夫逊法潮流计算,极座标法系数推导,当ij ，对特定的j，只有特定节点的j，从而ij= i- j 是变量,(4-49a),对特定的j，只有该特定节点的Uj是变量,(4-49b),第四节 牛顿拉夫逊法潮流计算,极座标法系数推导,当i=j ，由于i是变量，从而所有ij= i- j 都是变量，可得,(

41、4-49c),相似地，由于Ui是变量，可得,(4-49d),Qi(0),i(k),第四节 牛顿拉夫逊法潮流计算 极座标法潮流计算的基本步骤： 1. 输入原始数据和信息：y、C、Pis、Qis、Uis、约束条件 2. 形成节点导纳矩阵YB=CTyC 3. 设置各节点电压初值Ui(0), i(0) 4. 将初始值代入(4-45) 求不平衡量 Pi(0), 5. 计算雅可比矩阵各元素（Hij、Lij、Nij、Jij）,6. 解修正方程(4-44) ，求,Ui(k),i(k),7. 求节点电压新值Ui(k+1) = Ui(k) -,Ui(k), i(k+1) = i(k) -,2,第四节 牛顿拉夫逊法

42、潮流计算 极座标法潮流计算的基本步骤：,8. 判断是否收敛：Max|,Ui(k) |, Max|,i(k) |,9. 重复迭代第4、5、6、7步，直到满足第8步的条件 10. 求平衡节点的功率和PV节点的Qi及各支路的功率 . . n * * S s = U s Ysj U j = Ps + jQ s j=1 n Q i = U i U j (G ij sin ij B ij cos ij ) j= 1 . . * * * * S ij = U i ( U i U j ) y ij + U i y i 0,第五节 P-Q分解法潮流计算,一、P-Q分解法原理,所谓P-Q分解法就是利用牛顿-拉夫逊

43、法修正方程的极标形 式，考虑了电力系统的一些特性（如网络参数XijRij，,BijGij，ij0。P ，Q U），得出的一种简化形式。,图 形 解 释,第五节 P-Q分解法潮流计算,二、P-Q分解法的修正方程式,重写极座标方程,(4-53),第五节 P-Q分解法潮流计算,简写为 进一步 计及cosij1, Gij sinij Bij,(4-54) (4-55),第五节 P-Q分解法潮流计算,(4-49a),(4-49b),(4-56a),(4-49c),(4-49d) (4-43b),(4-56b),第五节 P-Q分解法潮流计算,(4-57),第五节 P-Q分解法潮流计算,(4-58a),(6

44、-75),(4-58b),P1/U1,B11,B12,B1n,U11,P2/U2,B21,B22,B2n,U22,(4-59a),Pn/Un,Bn1 Bn2,Bnn,Unn,Q1/U1,B11,B12,B1m, U1,Q2/U2,B21,B22,B2m, U2,(4-59b),Qm/Um,Bm1 Bm2,Bmm, Um,简写为：,P/U=BU Q/U=B U,(4-60a) (4-60b),P-Q分解法的修正方程式的特点：,以一个(n-1)阶和一个(m-1)阶系数矩阵B、B 替代原有的(n+m-2)阶系数矩阵J，提高了计算 速度，降低了对存储容量的要求。,以迭代过程中不变的系数矩阵B、B替代变

45、化 的系数矩阵J，显著地提高了计算速度。,以对称的系数矩阵B、B替代不对称的系数矩 阵J，使求逆等运算量和所需的存储容量大为 减少。,牛顿拉夫逊法 和PQ分解法的 特性：,PQ分解法,牛顿拉夫逊法,三、P-Q分解法的潮流计算的基本步骤,1) 形成系数矩阵B、B ，并求其逆矩阵。,2) 设各节点电压的初值I(0)(i=1,2,n,is)。UI(0)(i=1,2,m,is) 3) 按式（445a）计算有功不平衡量PI(0)(i=1,2,n,is)。,4) 解修正方程式，求各节点电压相位的变量 I(0)(i=1,2,n,is) 5) 求各节点电压相位的新值I(1) = I(0) + I(0)(i=1,2,