33复杂电力网潮流计算的计算机解法

1、3.3复杂电力网潮流计算的计算机解法导纳矩阵的形成1 自导纳节点i的自导纳，亦称输入导纳，在数值上等于在节点i施加单位电压，其他节点全部接地时，经节点更具体地说,i注入网络的电流。扎就等于与节点 连接的所有支路导纳的和。主对角线兀素-| -2. 互导纳节点i、j间的互导纳，在数值上等于在节点i施加单位电压，其他节点全部接地时，经节点 j注入网络的电流。非对角线元素_. _ ,1,f- °更具体地说，.是连接节点j和节点i支路的导纳之和再加上负号而得。3. 导钠矩阵的特点：(1) 因为怜,导纳矩阵Y是对称矩阵；(2) 导纳矩阵是稀疏矩阵，每一非对角元素 兀.是节点i和j间支路导纳的负值

2、，当i和j间没有直接 相连的支路时，即为零，根据一般电力系统的特点，每一节点平均与3-5个相邻节点有直接联系，所以导纳矩阵是一高度稀疏的矩阵；(3 )导纳矩阵能从系统网络接线图直观地求岀。4. 节点导纳矩阵的修改(1 )从原有网络引出一支路，同时增加一节点，设i为原有网络结点，j为新增节点，新增支路ij的导纳为yj。如图3-17 (a)所示。因新增一节点，新的节点导纳阵需增加一阶。且新增对角元Y =yj，新增非对角元Y=Y= yj，同时对原阵中的对角元 Y进行修改，增加 AY, = yj °(2) 在原有网络节点i、j间增加一支路。如图3-17 (b)所示。设在节点i增加一条支路，由

3、于没有增加节点数，节点导纳矩阵Y阶次不变，节点的自导纳 Yi、Y和互导纳Yj分别变化量为斗y冲話(3-57)钙=卜'(* =_冶图3-17网络接线的变化图(a)网络引出一支路，(b)节点间增加一支路，(c)节点间切除一支路，(d)节点间导纳改变(3) 在原有网络节点i、j间切除一支路。如图3-17 (c)所示。设在节点i切除一条支路，由于没有增加节点数，节点导纳矩阵Y阶次不变，节点的自导纳 Yi、Y和互导纳Yj分别发生变化，其变化量为3-17 (d)所示。（4）原有网络节点i、j间的导纳改变为|v-r|o如图设节点i、j间的导纳改变为一丁，相当于在节点i、j间切除一条yj的支路，增加一

4、条(3-59 )（5）原有网络结点i、j间为变压器支路，其变比由 增一变比为K的变压器。K变为K',相当于切除一变比为K的变压器，新(3-60 )r = ?r十科当节点之间变压器等值电路如图（a）、（ b）时，该变压器变比的改变将要求节点ij有关元素作如下修改。骨片D1 L图3.18由导纳表示的变压器等值电路（a）导纳在低压侧，（b）网络等值电路，（c）导纳在高压侧，（d） 网络等值电路导纳阵的相应元素如下变化：(3-61 )A0吋乂二一当节点之间变压器等值电路如图（c）、（d）时，该变压器变比的改变将要求与节点ij有关元素作如下修帆丈一 g3.导纳矩阵的计算1）计算流程（1）导纳矩阵

5、的阶数等于电力系统网络中的节点数。（2）导纳矩阵各行非对角元素中非零元素的个数等于对应节点所连的不接地支路数。（3）导纳矩阵的对角元素，即各节点的自导纳等于相应节点所连支路的导纳之和。（4） 导纳矩阵非对角元素等于节点 与节点 之间的导纳的负数。例3-3已知5节点系统单线图如3-19所示，已知数据如表3-1、3-2、3-3所示，母线1与发电机相联，发电 机G惨数400MVA 15kV，选为平衡节点。发电机 G2参数800MVA 15kV，选为电压控制母线（PU节点），母线 3与发电机G2和负载相联，母线2、4、5为PQ节点。写岀各节点已知量和待求变量的关系，计算系统导纳矩 阵。母线类型U古度P

6、GQGPLQLPGmaxPGmin1平衡1.00002负荷008.02.83电压常量1.055.20.80.44.0-2.84负荷00005负荷0000容量基准值，Sb=100MVA母线1、3电压基准值=15kV,母线2、4、5电压基准值UB=345kV表3.2例5线路输入数据（线路参数均为标幺值）母线-母线R/(V长度M最大值MVA2-40.00900.10001.7220012.02-50.00450.05000.8810012.04-50.002250.02500.445012.0表3.3例5变压器输入数据（变压器参数均为标幺值）母线-母线RX变比容量MVA最大值MVA抽头最大值设置1-

7、50.001500.0215/345kV4006003-40.000750.01345/15kV8001000解：输入数据和待求变量列于表 3.4。对于母线1,选为平衡节点，Pi和Q是待求变量。对于母线3，电压 受控母线（PU节点），Q3和待求变量。母线2、4和 5，与负荷相联（PQ节点），U2、U4、U和、是待求变表3.4例3-3母线输入数据和待求变量母线输入数据待求变量1U = 1.0,务=0P1, Q 12P2 = P G2- P L2 = -8Q = Q G2- Q L2 = -2.83U3 = 1.05P3= PG3- P L3 = 4.4Q,晒4P4 =0, Q 4 =0U4, %

8、5P5 =0, Q 5 =0U5,邑计算导纳矩阵：导纳矩阵Y的元素可由自导纳和互导纳的定义计算得到，以母线2为例写岀互导纳与自导纳的计算式，由于母线1和3不是直接连接到母线2，所以Y21 = Y 23 = 0-1 -1抵 0.89276+ J9.919M-9.9597295.143"R；4+M20.009+>0.1-1 -1得以s 二二二-1.78552 十J19.S3932 = M 十丿0.0045 七0,05y T I _+ 严1 TO fl QQ (0.89276- J9.91964) + (l 78552-J19.83932) +j + J-2 22.67828j；28

9、.4590=2B.5fr47-84.624"其中，连接到母线2的每条线路的并联导纳的一半包含在丫22中（另一半置于这些线路的另一端）同理可计算出导纳矩阵其他元素。S岭牯21：!£耳'3.73-J49.72000-3.73 + J49.72-02,68 - >28.4(50-0 89+J9.92-1.79 + J19.84=007.46-j99.44-7 46 + J99.44007屈十用92-7.46 + J99.4411.92-J147.96-3,57 + j39.CS-3.73+J49.72-1.79+J19.840-3 57 + j39.689.09-j

10、lOS.58_高斯-赛德尔法高斯-塞德尔法潮流计算（1 ）功率方程的特点描述电力系统功率与电压关系的方程式是一组关于电压的非线性代数方程式，不能用解析法直接求解。（2）迭代计算式如式（3-65 ）中的閃以乜替代（i=1,2,3 ）,就可用以解非线性节点电压方程(3-66)显然，式（3-66 ）中的.就对应于式（3-67 ）中的二.，叼就对应于巧就对应于呂一，內就对应于厂；或。但需指岀，按式（3-67 ）进行迭代时，除平衡节点外，其他节点的电压都将变化，而这一情况不符合PV节点电压大小不变的约定。因此，每次迭代求得这些节点的电压后，应对他们的大小按给定值修正，并据此调整这些节点注入的无功功率。这

11、是潮流计算运用高斯-塞德尔法时的特殊之处。（3）高斯-塞德尔潮流计算算法假设有n个节点的电力系统，没有 PU节点，若平衡节点编号为1,功率方程可写成下列复数方程式：(3-68 )对每一个PQ节点都可列出一个方程式, 定的，平衡节点电压也是已知的，因而只有高斯-塞德尔迭代法解潮流公式如下：因而有n-1个方程式。在这些方程式中，注入功率P和Q都是给 n-1个节点的电压为未知量，从而可以求得唯一解。上式可展开为:炉）=渤駕浮-（纬心4匕肿亠斗纭呼）、 诸Z厝-肉玩亠岭卫严亠+岭肿），y33 L 5一亏孕-区禹7必严）十巴严诃叽丄式中U是平衡节点的电压，k为迭代次数，上式是按高斯-赛的法解方程式组的标

12、准是书写的，对于 PQ节点，由于其功率是给定的， 故只要写出节点电压初值刖0）,即可利用（3-69 ）式迭代计算各节点节点电压。式中等号右侧的U采用经k次迭代值，等号右侧的 Uj，当jvi时，采用经（k+1）次迭代后的值，当j>i时, 采用经k次迭代后的值。迭代过程可进行多次，当某次迭代的解与前一次迭代后的解相差小于事先给定的允 许误差£时，即Jb"左附| m （j 2乡，璋），迭代终止，这就是迭代收敛的条件。一般系统内存在 PU节点，这种PU节点注入的无功功率受电源供应无功功率的限制。假设节点p为PU节点，设定的节点电压为 二、1，因其无功功率是未知量，只能在迭代开

13、始时给定初值必须在逐次迭代的过程中计算得岀。假定高斯-塞德尔迭代法已完成第 k次迭代，接着要做第 k+1次迭代前,先按下式求岀节点p的注入无功功率：101 （呼土与百严J1(3-70 )然后将其代入下式，求出节点 p的电压：-i w(3-71 )在迭代过程中，按上式求得的节点p的电压大小不一定等于设定的节点电压 中，应以设定的.对电压进行修正，但其相角仍保持上式所求得的值，使得心f V'，所有在下一次的迭代如果系统中有多个PU节点，可按上述相同计算方法处理。在迭代过程中往往求得 PU节点的无功功率会出现越限，即按式（3-70 ）求得的十-，不能满足约束条件 “一 '廿时，考虑到

14、实际工程中对节点电压的限制不如对节点功率的限制严格，这时可用纸或纸代入式（3-71）计算碍"】），此时不再需要修正电压的数值。换言之，这时只能满足约束条件Q眄达5©，而不能满足约束条件 =定值。事实上，此时该节点已由 PU节点转化为PQ节点。Lj-1i和节点j，在节点i测量支路电流规定由节点i流向节点j时为正求取线路潮流，线路连接节点其值为I1f （ 3 72 ）同理在节点j测量支路电流Iji规定由节点j流向节点i时为正。其值为:I厂一心+5三打十刀叼（3-73）复功率Sj表示又节点i流向节点j，Si表示由节点j流向节点i。其值为：去=巧7=6【5 % +丽磅剧=W 皿+q

15、 （U；-町）（3-74）-（3-75）以及各线路的功率损耗可由下式算岀：图3-20计算线路潮流的线路模型（4）高斯-塞德尔迭代法计算潮流的步骤：1）设定各节点电压的初值，并给定迭代误差判据；2）对每一个PQ节点，以前一次迭代的节点电压值代入功率迭代方程式求岀新值；3） 对于PV节点，求出其无功功率，并判断是否越限，如越限则将PV节点转化为PQ节点；4） 判别各节点电压前后二次迭代值相量差的模是否小于给定误差，如不小于，则回到第2步，继续进行 计算，否则转到第5步；5）根据功率方程（3-5）求出平衡节点注入功率；6）求支路功率分布和支路功率损耗。需注意：按高斯-塞德尔法进行迭代时，除平衡节点外

16、，其它节点的电压都将变化，这一情况不符合PU节点电压大小不变的约定。 因此，每次迭代求得这些节点的电压后，应对PU节点电压的大小按给定值进行修正，并据此调整这些节点注入的无功功率，如上面算法中所述。这是潮流计算中，运用高斯-塞德尔法时的特殊之处。图3-21高斯赛德法潮流计算流程例3-4利用高斯-赛德尔法计算例3-3系统潮流分布情况。解：对于例3-3所示的电力系统，用高斯-赛德尔法计算时，首先需要对各节点赋初值，之后除去平衡节 点外，按从小到大编号的节点进行迭代计算，迭代收敛后计算平衡节点的功率及网络损耗。（1）赋初值对 PQ节点赋电压初值：U=1.0/0° U=1.0/0°

17、 U5=1.0Z0°对PU节点赋电压（相位）初值： U=1.05Z0°（2）迭代求解PQ节点电压、PU节点电压相角和无功功率二 09215- J0.2737 = 0.9613Z-16.5WQ叽血耳口切二坯摯士写巧）>-1= m1.05x04-0+(7.46-j99.44/x(1.05ZO+(-?.46+J99.44)(1.0ZOj=5.220 a Qg - Qn = 36（射+珀普+环呼弋普卜44 少-fo+ 0+一246 十 JS9 轴)x 1. QZ0°+ 0 7.4-jS9.44 L L05Z0° L二 1.0374+/0.0393 = 1

18、0381Z2J7130由于求得的不等于给定的U3,将.修正为1 05Z2.17l3o-l JJ492+用茁98斷兀炉十车即十窿訓>70上工-卩/七砂.或 bOL 先 1址-1血乂打+(-丁心亠丹 9.+0j.l.(K4.m3%C357rjffibl-1. UJ544 JU.UU= J. UL554ZJC.4174°y ro_ylK53 占諜Tf"373十拠卫>d 0Z叫(T下+j«,3fl)c(IJSlNTH5«CP+a+FmJ.OKP-jO Q碑T QD& 1Z-1713P求得各节点电压新值后， 再按式计算I，开始第二次迭代，各节点

19、电压（标幺值）迭代结果示于表3.5,迭代过程中PU节点无功功率（标幺值）示于表3.6。手算时迭代误差可适当设置，应用计算机程序求解时，误差一般设为10-5表3.5例3-4迭代过程中各节点电压（标幺值）迭代次数kUU3U4U50I.OOOO+jO.OOOO1.0500+j0.00001.0000+j0.00001.0000+j0.00001.0000 ZO.OO °1.0500 ZO.OO °1.0000 ZO.OO °1.0000 ZO.OO °10.9215-j0.27371.0492+j0.03981.0354+j0.00751.0049-j0.04

20、780.9613 / -16.54 °1.0500 Z2.17 °1.0354 Z0.42 °1.0061 Z -2.72 °20.8577-j0.26711.0489+j0.04901.0323+j0.00140.9921-j0.04890.8983 / -17.29 °1.0500 Z2.68 °1.0323 Z0.08 °0.9933 Z -2.82 °30.8302- j0.28341.0491+j0.04291.0272-j0.00420.9853-j0.05400.8773 / -18.84 °

21、;1.0500 Z2.34 °1.0272 Z -0.23 °0.9868 Z -3.14 °40.8116- j0.28851.0493+j0.03711.0243-j0.00980.9808-j0.05700.8614 / -19.57 °1.0500 Z2.03 °1.0243 Z -0.55 °0.9825 Z -3.33 °经过49次迭代，最大误差精度为 8.9861e-6<1.0e-5490.7708-j0.31781.0499-j0.01091.0181-j0.05040.9712-j0.07720.83

22、38 / -22.40 °1.0500 Z - 0.59 °1.01 93Z -2.83 °0.9743 Z -4.55 °表3.6例3-4迭代过程中PU节点无功功率（标幺值）迭代次数k012349Q3.60001.32511.59032.12522.9747（3）求平衡节点1的功率 迭代收敛后，就可计算平衡节点的功率,5 = 2；(；;VJ)* = l Lr(r?x(3 73->4972)x1+ 04-0+0+(-3 73 J4S 72)x(0 9712-;U.U/72)J* = 35473+/I 1427(4) 求各条支路的功率及损耗利用式(3

23、-74 )、( 3-75 )及(3-76 )计算支路功率及支路损耗。以支路1-5为例求解。F -十儿 1441OLOOL5+JO(E)吊-马遲 _ 甲 珀十砖彳-1x0 < h：l - (19712- 7(110772)1 筑-393D5->D 3066电呱弘/!国&姑-硏血4耳妬-辽)(0.?712-0772X002-jQC7Z3T-lx笛厂& + 骂厂 3 国 5S+yl.144U-39aS-;ttQ02534 川昭t其余支路求解过程略，结果见表3.7表3.7 3-4各支路功率及损耗（标幺值）支路i-j耳%1-53.9458+j1.1441-3.9205-j0.

24、80650.0253+j0.33762-4-2.9185-j1.39103.0369+j1.21530.1184-j0.17572-5-5.0814-j1.40905.2564+j2.63010.1750+j1.22113-44.4008+j2.9747-4.3816-j2.71880.0192+j0.25594-51.3449+j1.5034-1.3344-j1.82510.0104-j0.3217网络总损耗为:i= (0.0253-Fj037d)+(Q11M-；11757)4-(101750+；12211)4-® 0192 + Q255ST)+(001M-；0 3217)= 0.

25、3433+71.3172以及这一网络的输电效率为：哗=>4Q0% = 96.219%至此本例潮流计算全部解算完毕。牛顿-拉夫逊法1.牛顿-拉夫逊算法原理牛顿-拉夫逊（Neton-Raphson ）法是求解非线性代数方程有效的迭代计算方法。在每一次的迭代过程中,非线性问题通过线性化逐步近似。以一个变量为x的非线性函数求解过程加以说明。设一维非线性方程-（ 3-77 ）求解x，设真值为x*。首先在x*附近选一初值x（0），则误差为沖二少-科式（3-77 ）写为齊-卅）=6（3-78）将上式展开成Talor级数,并丄皿牛打一列禺十卢理)冷如果初值x(0)接近真值，则误差足够小，可略去上式中的高

26、阶项(3-79)可得(3-80)广严)将x(0)代入上式，求得误差修正量，即可得到所求解。(3-81)如此继续下去，则可得到充分逼近解:3-82)(3-83)Ax3-84)理论上收敛条件(3-85)图3.22( a)中示出牛顿-拉夫逊法的解算过程，可见的接近于精确解，否则迭代过程可能不收敛，如图3-22x(v+1)更接近于真值。运用这种方法时，初值要选取(b)所示。（a）初始值选取合适收敛（b）初始值选取不合适不收敛对于n维非线性方程组(3-86 )Axp =-X* J =1, 2r环申云町轧M詡一两u （3-87） 石岸匸闵时却酬屮&竹=0I工岸）一阿阈却心怛屮出紳展成Talor级数

27、，并略去二阶以上项B'却打璽P利° T就0 .=0（3-88）粼絆津出卜疋）十亟,盹虜）十十並o丸0 =0>K川利岀却卜莹9禺関）十十凱0兀» -=0整理成为如下的矩阵方程:d(徵申用)-国罢鱼 矗% 1 4«(3-89 )式（3-89 ）等号右边矩阵中的是分别对于xi,x 2,.,xn求导的值，这一矩阵称为雅克比（Jacobi ）矩阵。上式简记为:己（0）=卅）Ax可解出4辭7罰求解形式如下：(3-90 )收敛条件（3-91）(3-92 )2. 直角坐标系下的牛顿-拉夫逊算法-赛德法时相同，不同的只是迭代过程。迭代过程中，两种方法应用的基本方程都是

28、¥U-，运用高斯-赛德法时，将运用牛顿-拉夫逊法计算潮流时，节点导纳矩阵的形成、平衡节点和线路功率的计算与高斯其展开为电压方程，运用牛顿-拉夫逊法时，将其展开为功率方程(3-93)“也-u, T啊L J司式中，第一部分为给定的节点注入功率，第二部分为由节点电压求得的节点注入功率，二者之差就是节点功率的误差，当节点功率误差趋近于零时，各节点电压即为所求方程的解。在采用直角坐标系下，节点电压和导纳可表示成式:（ 3-94）卩二务+止将上式代入式（3-93），展开取出实部和虚部，得到式:L土 吗勺-3血 4 VP 严 1J-1Ql灶g昭）-罚丈（G+ S）>1J-3(3-95)根据节

29、点分类，若第i个节点为PQ节点，给定功率设为 刃和C功率误差方程可列为:伞胡-乌-昭理炉利-电门-倉炉訂广¥山02J"也二2 -Q二3 煌®厂仆吐阳十巧引二o问>1(3-96)若第i个节点为PU节点，：-和i给定，功率和电压的误差方程可列为:臂乂厂£ = & -色£（&护厂塢吾）-血G羔+B洱）=0J-1>-1(3-97)式中，节点电压大小（模数）的误差表示为给定的节点电压的平方与求得的节点电压平方之差。上述功率和电压误差方程即为牛顿-拉夫逊潮流计算所要求解的非线性方程组。非线性方程组的待求量为各节点的电压的实部 e

30、i和虚部fi。对于n个节点的系统而言，第s=n号节点为平衡节点，除去平衡节点电压为已知外，式（3-96 ）和（3-97 ）共包含2 （n1）个方程，待求变量也有2 （ n- 1 ）个。除去平衡节点外 的所有节点有功功率不平衡量的表达式有（n-1 ）个，即i=1，2，n，;除去平衡节点外的所有节点无功功率不平衡量也jg的表达式有（m-1 ）个，即i=1，2,m ；所有PU节点电压大小不平衡量的表达式有（n-1 ） - （m-1） =n-m个，即i=m+1，m+2n，;平衡节点功率和电压方程不包括在这组方程之内，其电压向量是给定的故不需要求取，当上述各点电压迭代收敛后再求取平衡节点的注入功率。把各

31、节点的电压变量用初始值与修正量的形式表示为将此关系代入到式（3-96 ）、（ 3-97 ）中，在附近的 、無范围内将其展开为泰勒级数并略去咼阶项，可得3V-W （3-98 ）式中：価心屿出辄叫叽 4也讯丿沁Q 龛1»恥Q¥n网% 恥Q% *AA%*«鳩+i4页也iB!兀丄*44!咤L1辱ds、眾%號+1亀一 i前一 1酬a.迢恥a.J-龛1%臂i沁恥氏a.恥沁环1.恥凡+】沁+1龛%、%(rj恥沁沁呱恥璟1龛臥驾！鸥！ 41各1彰.MJai叱i啤1鬼%Bg-1驰几处g沁*、-弼%沁其中，雅克比矩阵J中的各元素可以通过对式（3-96 ）和（3-97 ）求偏导数得到

32、当j=i时，对角线元素为：(3-99 )rtAF “= -£ （爲勺-G禹-% j-iSAP. «= - £趕）+ B 耗- 】卬!/-1Qhf）, K=工灿心+零/+隅。屈J如 A1学二 一£（% - 57料 + gdst = 闫遊=一叮曲当 时，非对角线元素为：(3-100)(3-101 )由以上表达式可得岀雅克比矩阵的特点：（1）矩阵中的元素是节点电压的函数，在迭代过程中将随着节点电压的变化而改变。（2）矩阵是不对称的。（3 ）当导纳矩阵中的非对角线元素.为零时，雅克比矩阵中相应的元素也为零。矩阵是稀疏的，可以应用稀疏矩阵的求解技巧。3. 采用极坐

33、标下牛顿-拉夫逊法潮流算法 以极坐标表示时，节点电压和导纳可表示为,q 二q評 “g皆ysm）< ! L 1功率误差方程可表示为(3-102)>1对一个具有n个独立节点，其中有（n-m-1 ）个PU节点的网络，式（3-102 ）组成的方程组共有（n-1） +m个方程式。采用极坐标时，方程组个数较采用直角坐标表示时少了（n-m-1）个。因为PU节点，采用极坐标时，待求的只有电压的相位角和注入的无功功率，而采用直角坐标时，待求量为电压的实数部分、虚数部分 和注入的无功功率，因此采用极坐标可使未知变量少了（n-m-1）个，方程数也少了（n-m-1）个，这样建立修正方程式的矩阵形式为(3-

34、103 )式中AP =鼻珂-QQ =AQafAA =3“5AU =| 是（n- 1）x（n -1） 阶方阵，S是（n- 1） Xm阶矩阵， 划.是 mX （n -1） I阶矩阵，是mXm阶矩阵。各矩阵种元素分别为(3-104 )丹，坐"，竺LM =直叫'JJ 1在式（3-103 ）中，电压幅值的修正量采用厂 的形式，是为了使雅克比矩阵中各元素有比较相似的pc j-尽=討=陋%沖殆-坷小毎）(3-105a)形式。矩阵中各元素可对式（3-83）（ 3-92 ）取偏导数求得，雅克比矩阵中各元素具有比较整齐的形式。岛二等6二一 56 ©血% -塢阿岛）3（3-105b ）q

35、（C?护in吗-cos吗）4 2耳鸟Jrl片】吗=等6 = 7曲2泅十血Slii如3AP«恥二春E TY埠+吗迪令）-血愆(3-105C )恋萨警二叩佩cos爲+ %创陽）(3-105d )式中门为i、j两节点电压相角之差3 以直角坐标系为例说明牛顿-拉夫逊法的程序计算步骤及流程图(1) 形成网络导纳矩阵 丫，设定各节点电压的初值 总卩)y(0)；(2) 将以上电压初始值代入式(3-96 )和(3-97)，求取修正方程式中的误差函数值(3) 将电压初始值再代入式(3-99)，求取雅可比矩阵中的各个元素；(4) 解修正方程式，求出节点电压修正量 厂二 ,<；(5) 修正各节点电压

36、，萃.严A捫(6) 将新值总/再代入式(3-96 )和(3-97 )，计算新的各节点功率及电压误差函数值 '十.-) -rr- "；(7) 检查计算是否收敛，当电压趋于真实值时，其功率误差趋于零。收敛判据为品聘 瓦忖或応翩 砰9詞,其中£为是先给定的小数。(8) 若收敛，迭代到此结束，计算平衡节点功率、 各支路潮流、损耗及输电效率，并输岀结果；若不收敛,则转回第 步，以 代替严进行下一次迭代，直至收敛。平衡节点功率计算如式(3-95 )，计算线路功率则调用式(3-74 )和(3-75 )，计算线路损耗调用式(3-76 )。牛顿拉夫逊法潮流计算程序流 程如图3-23所

37、示。利用极坐标系计算潮流的过程类似。1檜出菇察)图3-23牛顿拉夫逊法程序流程图运用牛顿拉夫逊法计算潮流时，由于初值要选取的比较接近于精确解，否则可能是迭代过程不收敛，实 际计算程序中，往往采用高斯赛德尔法与牛顿拉夫逊算法配合使用的方案，即在前几次迭代时采用高斯赛德 尔法，得到牛顿拉夫逊算法的初值，之后再利用后一种方法求解。若计算过程中每次迭代计算得到的 x变化不大，也可以经多次迭代后才重新计算一次雅克比矩阵各元素。 因此牛顿-拉夫逊法获得了广泛的应用。例3-5利用牛顿-拉夫逊法直角坐标方式计算例 3-3所示网络潮流分布情况。解：确定例3-3系统雅可比矩阵的维数。系统有n= 5条母线(节点)，

38、采用直角坐标方法求解时组成 2(n-1) =8个方程，J(i)维数为8X8。按题意要求，该系统中，节点 1为平衡节点，保持 U=1+jO为定值，2, 4，5为 PQ节点，3为 PU 节点，U3=1.O5+jO。(1)赋初值由已知可知平衡节点：岂(,£对PQ PU节点赋电压初值：評屮).屮=1.炉=0,<-L05f 炉-0时=& - $ = *曲刨7肿-曲Z +*)= -S.D-l OxD + (2.783xl.00)+ 0+(0 3520x1 十。 =-3.055站=2 - G化Q广炉工縄鉀-亏多“第乞色忙加J-JH-2.3-(0 1.0x0 +()tOt(0 i-95

39、l?7x1.0)+(0i-)J - 1.5砂=.炉詔厂甥£ (陽畀p#卜炉去务护q陽叨)j-i问-44-(1.05ko+0+(7.4580x1 05-C)4(7.£530x1 .0-0)0+0)-4.D0S5A?fJ2 =苗-u严-耳？-他曲+护"卜1防- 1+ n) = n时=廿胪=耳厂半壬妙膚'珂学)-炉$(知+時忙x 05-0)4-(：1.929xl.0-0)+(-3.5711 xl.D-Q)J + O = 0.3729- a-曲-a -用T %娜)卜屮丈趴外窗j-JH= 0-0-：+ (0-1479529x1 0) +(0 + 39.5768x1

40、0) = d 05255碑七-叩=E -*1戈cv畀迢护xrn爼炉+册)J-lJ-l=U - 1.0 (-3.7230X 1 J - 0)+1 - 1.7S55x 1-U - J +0+(-3.5711x1.0- 0) + 9.085dxl D - 0)+ 0)= D 缈迢-Q犁应G再叽坯炉"理虫惩甲谒创巧= O-(0-l.0x(Oi-4?.7ZO3KL.O|+(OlW93xL0)+0+(0+39.ff7®d xl. 0) + (0-108 5732x1 0)J = QCd(3 )计算雅可比矩阵以节点2 ( PQ有功、无功功率和节点 3 ( PU电压幅值分别对各节点电压实部

41、、虚部求导为例，其他节点的求解过程略。讐卜-号 - D-G於-BM = -2.6783 韻+也亠 “297590憎|闵1 -G蛉-Bf2 08阀第1=十2 = °讐|= 为=0.8928= -G»E 十為勺=9.9197d讐卜-G遏-鸟八7855j 一內石 + 19 *393陋Q龛2一士勺卜甩卄-27 1590Q J-1沁2-=-£（令再-夠/） +十% = 2.G7830 J"SB卄签=0o黔卜2 m-叱=0lo= 9.9197-1 =-1 = -0.8928 筋I 叫Io= 19.8393观| 二-洱=-1.7855% I 施L哲纠=-23 =-2.

42、1000殉L10asaQITaQ讯目XI £-25783-29.7590-27J5902 £7 83II0J3239.91979.3157-OS92S1.705519 8393弘2莎 會 彰| 01訳丸歸10筋1 禺0臥阳k钳|J娥0您10°赵L1竺10 睥Qo珂酬Q1昭L須0L割I00.39289.91970弘Q-a2O39-210007.4520994406(.9.9197-1S92®-39+4067.83101C4.4127(J9? 4406-7 45S好一二QDk-Q3231姐*)0=-於戸0.0037车0>)0.0423也-03810 0

43、1160-0.073019,3393-17855求解各节点新值，_ 1.1二一.人.-I,二._：一-：1.5430-141993 5711IHo-：54.Q1Q912 294S39.6706-3.571 110IQ1.735519.83933.57H29 67S69.0856-107.918219.8353-173550039,刃跖 25711-109 23E29 Q85&迭代次数k节点电压U2U3U4U510.9429-j0.32311.0500+j0.00371.0423-j0.03811.0116-j0.0730收敛判断，不满足收敛要求，进入下一次迭代。每次迭代所得结果示于表3.8-3.9表3.8例3-5迭代过程中节点电压变化情况20.8075-j0.31801.0500-j0.00801.0229-j0.04820.9797-j0.07643 二0.7734-j0.31781.0499-j0.01081.0184-j0.05020.9718-j0.07724 二0.7708-j0.31781.0499-j0.01091.0181-j0.05040.9712-j0.077350.7708-j0.31781.0499-j0.01091.0181-j0.