电力系统的基本计算.pdf

《电力系统的基本计算.pdf》由会员分享，可在线阅读，更多相关《电力系统的基本计算.pdf（32页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第二篇 电力系统的基本计算 本篇进入电力系统分析课程的核心部分电力系统的基本计算。内容包括第三章电力系统的潮流计算、第四章电力系统的故障计算、第五章电力系统的稳定计算。常称电力系统的三大常规计算潮流、短路和稳定，尽管短路只是电力系统故障的一种。电力系统的基本计算既有自身独立的意义，又是电力系统设计、运行和研究的理论基础，其重要性自不待言。第三章 电力系统的潮流计算 电压（包括幅值U和相位）和功率（包括有功功率P和无功功率Q）是表征电力系统稳态运行的主要物理量。所谓电力系统的潮流计算就是采用一定的方法确定系统中各处的电压和功率分布（实为功率流，power flow，但电力界惯称潮流）。电力系统的

2、潮流计算和一般交流电路计算的根本差别在于：后者已知和待求的是电压与电流，而前者是电压与功率。正是这一差别决定了二者本质上的不同：描述交流电路特性的方程，如节点电压方程、回路电流方程，是线性方程，而描述电力系统稳态运行特性的潮流方程是非线性方程。以一条阻抗为 Z 的支路为例，描述其电路特性的方程IZU&=是线性方程，其中电压U&和电流I&之间的关系是线性关系；如果已知和待求的是电压与功率，因功率与电流之间的关系为=IUS&，则描述其特性的方程成为)/(*USZIZU=&，从而电压U&和功率 S&之间是非线性关系。由此使得求解方法有了根本不同：线性方程可直接采用消去法求解，而非线性方程只能采用迭代

3、法求解。这就是电力系统潮流问题的特点：已知和待求的是电压与功率，为非线性关系，需迭代求解。电力系统中进行电力系统潮流计算的目的在于：确定电力系统的运行方式；检查系统中的各元件是否过压或过载；为电力系统继电保护的整定提供依据；为电力系统的稳定计算提供初值；为电力系统规划和经济运行提供分析的基础。可见，电力系统的潮流计算是电力系统中一项最基本的计算，其既有一定的独立的实际意义，又是研究其它问题的基础。本章介绍利用计算机进行电力系统潮流计算的原理和方法。利用计算机解题，一般包括建立数学模型、确定解算方法、制订计算流程图和编程上机等步骤。本章主要介绍前两步：潮流计算的数学模型潮流方程和潮流方程的迭代求

4、解方法，至于编制程序上机计算可在上机实践环节中进行。另附带讨论潮流计算中的有关技术，包括潮流方程迭代求解时的初值设定、线性方程组的求解、稀疏技术简介和网络化简。应强调提出，用计算机求解潮流时，均采用标幺值。第一节第一节 潮流计算的数学模型潮流方程潮流计算的数学模型潮流方程 本节从电力网络的标幺值电压方程入手，介绍节点导纳矩阵和节点阻抗矩阵的性质、形成和修改方法，进而引出潮流计算的数学模型潮流方程。一、电力网络的节点电压方程一、电力网络的节点电压方程 电力网络是一种电路，因而求解电路的方法，如回路电流法、节点电压法和割集法等，原则上均可用于电力网络。但实际中割集法几乎不用于电力系统，回路电流法用

5、得也很少，广泛应用的是节点电压法，所以此处仅介绍节点电压法。节点电压方程形如 BBBUYI=（31）式中，下标B代表节点，是英文单词 bus 的第一个字母，因在电力系统中常以发电厂和变电所的母线（bus）作为节点；IB为节点注入电流列向量，注入电流有正有负，注入网络的电流为正，流出网络的电流为负。根据这一规定，电源节点的注入电流为正，负荷节点为负，既无电源又无负荷的联络节点为零，带有地方负荷的电源节点为二者之代数和；UB为节点电压列向量，由于节点电压是相对于参考节点而言的，因而需先选定参考节点。在电力系统中一般以地为参考节点，如整个网络无接地支路，则需选某一节点为参考。设网络中节点数为 n（不

6、含参考节点），则 IB、UB 均为1n 列向量；YB为 nn 阶节点导纳矩阵。节点电压方程可写成另一形式 BBBBBIZIYU=1 （32）式中，1=BBYZ称为节点阻抗矩阵。注意上式与回路电流方程 LLLIZU=是本质完全不同的两类方程，不应混淆。从而节点阻抗矩阵 YB和回路阻抗矩阵 ZL也有本质的不同。二、节点导纳矩阵二、节点导纳矩阵 YB 节点导纳矩阵 YB是nn 方阵，其对角元iiY（ni,1L=）称为自导纳，非对角元ijY（=i1，n，ji）称为互导纳。下面介绍节点导纳矩阵各元素的含义、性质及其形成和修改方法。1.节点导纳矩阵各元素的意义和性质 将节点电压方程 BBBUYI=展开为=

7、nnnnnnnYYYYYYYYYIIILMMLL&M&21222211121121 nUUU&M&21 （33）可见 0/=jUiiiiUIY&nji,1,L=；ji （34）表明，自导纳 Yii在数值上等于仅在节点i 施加单位电压而其余节点电压均为零（即其余节点全部接地）时，经节点i注入网络的电流。其显然等于与节点i直接相连的所有支路的导纳之和。注意不直接相连支路的导纳不应包括在内。同时可见 0/=iUjiijUIY&nji,1,L=；ij （35）表明，互导纳ijY在数值上等于仅在节点j施加单位电压而其余节点电压均为零（即接地）时，经节点 i 注入网络的电流。其显然等于（ijy），即ijY

8、=ijy。ijy为支路 ij的导纳，负号表示该电流流出网络。如节点ij之间无支路直接相连，则该电流为 0，从而ijY=0。注意字母 y 几种不同写法的不同意义：粗体黑字 Y 代表导纳矩阵，大写字母ijY代表矩阵 YB中的第i行第j列元素，即节点i和节点j之间的互导纳，小写字母ijy代表ij支路的导纳，等于支路阻抗的倒数，ijy=ijz/1。根据以上定义，可以容易地直接形成节点导纳矩阵，下面用一简单实例说明之。【例例 3-1】形成如图 31 所示系统的节点导纳矩阵。该系统是一个由三条输电线组成的环形网络，输电线用形等值电路表示。设三条线路参数的标幺值均相同：1.0jzL=，02.0jyL=。求系

9、统的节点导纳矩阵。解解：选地为参考节点。以节点 1 为例说明自导纳iiY 的形成。和节点 1 直接相连的支路有：支路 12 的阻 抗支路Lz，支路 13 的阻抗支路Lz以及和节点 1 直 接相连的两条并联导纳支路2Ly。将 Lz表为导纳，从而 98.1901.001.01.011.0111jjjjjY=+=图 31 例 31 网络图 以节点 1 与节点 2 之间的互导纳12Y为例说明互导纳的形成。1、2 节点间有直接支路，其导纳为Lz1，故 1011212jzyYL=照此办理，得到此系统的节点导纳矩阵=98.1910101098.1910101098.19jjjjjjjjjYB 节点导纳矩阵具

10、有如下性质：YB为对称阵，ijjiYY=。如网络中有含源元件，如移相变压器，则对称性不再成立。对无接地支路的节点，其所在行与列的元素之和均为零，即=jijiijYY0,0；对有接地支路的节点，其所在行与列的元素之和等于该点接地支路的导纳。利用这一性质，可以检验所形成节点导纳矩阵的正确性。YB具有强对角性：对角元的值不小于同一行或同一列中的任一元素。以上三点性质可在上例中得到验证。YB为稀疏阵，因节点i、j之间无支路直接相连时ijY=0，这种情况在实际电力系统中非常普遍。矩阵的稀疏性用稀疏度表示，其定义为矩阵中的零元素数与全部元素数之比，即 2nZS=（36）式中，Z为BY中的零元素数，S随节点

11、数 n 的增加而增加：n为 50 时，S可达 92%；n为100 时，S达 96%；n为 500 时，S达 99%。充分利用节点导纳矩阵的稀疏特性可节省计算机内存，加快计算速度，这种技巧称为稀疏技术。节点导纳矩阵 YB的形成除了上述的直接方法外，还可利用支路节点关联矩阵得到。此处不再介绍，有兴趣者可参阅有关文献。2、节点导纳矩阵的修改方法 电力系统的结线方式经常发生变化，从而节点导纳矩阵也随之发生变化。人们经过实践发现，系统结线情况变化后，不必完全重新形成节点导纳矩阵，而只需在原节点导纳矩阵的基础上稍作修改即可。这样可大大减少重复工作量。下面就介绍这类修改方法。电力系统结线方式情况的变化分为两

12、大类：一类是新增一条支路，如新建一条输电线给一个新负荷供电；另一类是原有网络某条支路的情况发生变化，如投入、退出一条支路或支路的参数发生变化。对应于这两类变化，节点导纳矩阵的修改方法亦分为两大类：第一类：从原网络节点i引出一新的支路，同时增加一个新的节点（编为第 n+1 节点）。此时节点导纳矩阵将增加一阶，从nn阶变为)1()1(+nn阶。由于节点n+1 只有这条支路和节点i相连，从而增加的对角元为 yYnn=+1,1（y为新增支路的导纳），增加的非对角元除yYYinni=+,11,外其余元素均为 0，同时原来的对角元iiY变为yYii+，即 i MLKMiiYii+yyyyYii00001M

13、MLLM nn )1()1(nn 上述修改方法由对角元和非对角元的定义直接得来。第二类：原网络i、j支路的参数发生变化。此时节点导纳矩阵 YB的阶数不变，仅需对其中的四个元素iiY、ijY、jiY 和jjY进行修正，即 ji jiYYYYjijjjiijiiMMLLLMLLLMM MMLLLMMLLLMMjjjiijiiYYYYji 式中：iiY=iiY+ijy，jjY=jjY+ijy，ijY=ijYijy，jiY=jiYijy，ijy为 ij支路导纳的变化量。当ijy=0ijy时表示该支路投入运行；当ijy=ijy0时表示该支路退出运行；当ijy=ijijyy时表示其支路参数发生变化，由原来

14、的ijy 变为ijy，如双回路输电线运行变单回线运行或反之。如某一变压器的变比由k改变为k，则由标准情况下变压器型等值电路（见图 32，其由图 232 得来，各支路参数表为导纳，TTZY1=，注意其不是变压器导纳支路的导纳 TTjBG），得到()()=TjiijTjjiiYkkyyYkkyy220 （37）图中节点i对应于变压器原方，j对应于副方。（三）节点阻抗矩阵BZ 节点阻抗矩阵BZ也是nn方阵，其对角元 iiZ（ni,1L=）称为节点自阻抗，非对角元ijZ（jinji=,1,L）称为节点i和节点j之间的互 阻抗。将BBBIZU=展开为 =nnnnnnnZZZZZZZZZUUULMMLL&

15、M&21222211121121 nIII&M&21 （38）可见 0=jIiiiiIUZ&jinji=,1,L （39）表明，自阻抗在数值上等于仅在节点i注入单位电流而其余节点均不注入电流（电源均开路）时节点i的电压。同时可见 图 32 变压器型等值 电路的导纳形式 0=iIjiiiIUZ&jinji=,1,L （310）表明，互阻抗在数值上等于仅在节点j注入单位电流而其余节点均不注入电流时节点i的电压。节点阻抗矩阵BZ在网络中无含源元件时也是对称阵，但与稀疏的节点导纳矩阵BY 不同，BZ是满阵。其原因在于当在某节点注入电流时网络的所有节点上均会感受到电压；或者从数学上看，稀疏矩阵的逆阵不再

16、是稀疏阵。与节点导纳矩阵不同的另一点是，iiY、ijY 均由具体支路的导纳组成，而iiZ、ijZ无具体支路阻抗相对应。形成节点阻抗矩阵的方法也有两类：一类是求逆法，即由定义：1=BBYZ，对已形成的节点导纳矩阵求取逆阵；另一类是根据自阻抗和互阻抗的定义直接一步步形成阻抗矩阵的方法，称为支路追加法。利用支路追加法还可方便地对节点阻抗矩阵进行修改。关于这方面的内容不再详叙，有兴趣者可参阅文献1。四、潮流方程四、潮流方程 前已提出，由于电力系统已知和待求的不是电流而是功率（原因是：对庞大的交流系统，电流相位的测定十分困难；而功率的测量十分方便，可由有功功率表和无功功率表得到），故将节点电压方程中的电

17、流代之以功率，*/USI=&，得到 YUUS=*/（311）即 =jjijiiUYUS&*（312）或 =jjijiiiUYUjQP&*ni,1L=（313）式中=njj1，下同。式（312）就是电力系统潮流计算的数学模型潮流方程。它具有如下特点：它是一组代数方程，因而表征的是电力系统的稳态运行特性；它是一组非线性方程，因而只能用迭代方法求其数值解；由于方程中的电压U&和导纳Y既可表为直角坐标，又可表为极坐标，因而潮流方程有多种表达形式极坐标形式、直角坐标形式和混合坐标形式。取iiiUU=&，ijijijyY=，得到潮流方程的极坐标形式 jjjijiiiiUYUjQP=（314）取 iiijf

18、eU+=&，ijijijjBGY+=，得到潮流方程的直角坐标形式()()()()+=+=jjijjijijjijjijiijjjijjijijijjijiieBfGefBeGfQeBfGffBeGeP （315）取iiiUU=&，ijijijjBGY+=得到潮流方程的混合坐标形式()()=+=jijijijijjiijijijijijjiiBGUUQBGUUPcossinsincos （316）式中 jiij=。不同坐标形式的潮流方程适用于不同的迭代解法：例如，利用牛顿拉夫逊迭代法求解时，以直角坐标和混合坐标形式的潮流方程为方便；而 P-Q 解耦法是在混合坐标形式的基础上发展而成，故当然采用混

19、合坐标形式。实际中真正极坐标形式的潮流方程用得很少，因而常将式（316）称为极坐标形式。它是一组n个复数方程，因而实数方程数为 2n，但方程中共含 4n个变量：iP、iQ、iU和i，ni,1L=，故必须预先指定 2n个变量才能求解。为将 2n个变量定为已知量，根据电力系统的实际情况，对每个节点指定两个变量，余下两个变量待求。通常将节点分为三种类型：PQ 节点、PV 节点和V节点。下面分别对这三类节点加以说明。PQ 节点：对这类节点指定P和Q，U和 待求。电力系统中绝大多数节点均属此类，如变电所母线节点，其无电源功率，负荷功率又已知，故该节点的节点注入功率DiDiijQPS=&已知；又如一些按指

20、定有功和无功功率发电的电厂，其GiP、GiQ指定，所带机端负荷也已知，从而节点的注入功率()()DiGiDiGiiQQjPPS+=&已知。PV节点：对这类节点指定P 和U，Q 和 待求。电力系统中此类节点属少数，个别小系统甚至没有。设置PV 节点是为了控制该点的电压为一定值从而保证系统的电压质量。为了控制电压必须要有一定的无功功率可供调节（其机理将在第七章阐述），故这类节点是有一定无功储备的发电厂和装有无功电源（电容器、调相机或静止无功补偿器）的变电所。这类节点也称为电压控制节点。V节点：、对这类节点指定U和，其有功功率P和无功功率Q由保证全系统功率平衡的条件确定，因而又称平衡节点。一般取其o

21、0=。电力系统潮流计算中必须有而且只有一个V节点，负责系统频率调整的主调频厂基本上起着平衡节点的作用。这三类节点各有不同的特点，在潮流方程的求解过程中有不同的处理方法，希引起注意。此外，由于电力系统在运行中必须满足一定的技术经济要求，如电压必须在允许范围内：miniUiU maxiU；各电源的功率必须在其所能发出的功率范围内：minGiP GiPmaxGiP，minGiQGiQ maxGiQ；某些节点电压间的相位差应在一定的范围内以满足系统运行稳定性的要求（其机理将在第五章阐述）：ijmaxij。这些便构成了潮流方程的约束条件。所以电力系统的潮流计算归结为求解一组非线性方程潮流方程，并满足一定

22、的约束条件。应指出，上述的节点分类方法并非唯一，还可有其它的分类方法。原则上讲，只要指定的变量总数为 2n且实际中有意义的方案均属可行，如P节点（只给定P）、U节点（均不给定）、PQV节点（均给定）等。第二节第二节 潮流方程的迭代求解 潮流方程是非线性方程，求解非线性方程的基本方法是迭代。有多种迭代算法：高斯迭代、牛顿拉夫逊迭代等。高斯（Gauss）迭代是最简单的一类迭代。在电力系统潮流计算中，牛顿拉夫逊迭代（Newton-Raphson 迭代，简记为 N-R 迭代）是占主导地位的有效方法，在其基础上结合电力系统的实际特点进行简化而成的 P-Q 解耦迭代（简记为 P-Q 迭代）得到了广泛应用。

23、下面分别予以介绍。一、潮流方程的一、潮流方程的 N-R 迭代迭代 先介绍一维情况下的 N-R 迭代，进而介绍n维情况下的 N-R 迭代，然后将其用于潮流方程的求解。1 一维情况下的 N-R 迭代 设非线性方程()0=xf，x为满足该方程的真解，其与所设初值)0(x 的差记为x，)0(xxx=。如x求出，则 xxx+=)0(。将0)()()0(=+=xxfxf在)0(x处展为台劳级数)()()()()0()0()0()0(xfxxfxfxxf+=+!22x+（317）如初值选择得当，x很小，则上式中的二次及以上高次项可略去，得到近似式 0)()()0()0(=+xxfxf （318）称为 N-R

24、 迭代的修正方程式，由其可得到修正量)0(x()()0)0(/)(xfxf=（319）注意到此时得到的)0(x并不是真正需要的x（因忽略了高次项），故()()1)0(0 xxx=+并不是真解x，只是向真解逼近了一步的改进值。以)1(x作为新的初值代入修正方程 ()()()()0)(=+kkkxxfxf L1,0=k （320）得到()1x，于是()()()112xxx+=，照此办理，当)(kx0 时便有()()0kxf，从而)(kx即为所求解。故 N-R 迭代的收敛判据为(k)x或()()kxf。N-R 迭代的核心是将非线性方程式的求解转换成相应线性修正方程式的多次求解。其迭代过程如图 33

25、所示。从直观上可见，N-R 迭代的收敛速度比 Gauss 迭代快，是一种滑梯式逼近过程。由于推导修正方程的前提是)0(x选择得当使x较小，故 N-R 迭代对初值要求较严，否则会不收敛。因而有时将 G-R 迭代和 N-R 迭代结合起来，利用前者得到较好的初值，然后转用 N-R 迭代加快收敛，能取得较好的效果。)2(x()1x()0 x()()()()xxfxf=00 图 33 NR 迭代【例 3-2】利用 N-R 迭代计算非线性方程0782=+xx的解。解：此时82)(=xxf，N-R 迭代公式为()()()()()()82/7821+=+kkkkkxxxxx L2,1,0=k 仍设初值为)0(

26、x=3，得到 )1(x=-1，)2(x=0.6，)3(x=0.9765，)4(x=0.9999084。四次迭代后5105.5)(=xf，收敛速度确实比较快。2n维情况下的 N-R 迭代 对n维非线性方程组()0=XF，其修正方程为()()XFXF+0=X （321）式中()=nnnnnnxfxfxfxfxfxfxfxfxfXF 212221212111LMMMLL()XJ=称为雅可比矩阵。迭代公式为 ()()()()()()+=+kkkkkkXXXXXJXF10 L2,1,0=k （322）收敛判据为()()kixfmax 。修正方程的求解常采用高斯消去法。3潮流方程的 N-R 迭代求解 前已

27、提及，利用 N-R 迭代求解潮流方程时，常采用直角坐标或混合坐标形式。此处仅介绍混合坐标形式（一般文献称其为极坐标形式），因下面将要介绍的另一种迭代方法 PQ 解耦迭代就是在其基础上结合电力系统的实际特点经过简化而得到。直角坐标形式潮流方程的求解原理和步骤与混合坐标形式基本相似，读者可自行推导之。将混合坐标形式的潮流方程（式（316）表为0)(=xf的形式，得到()()=+=jiiijijijijjiiiijiijijijijjiiiQQBGUUQQPPBGUUPP0cossin0sincos （323）上式的含义是：求一组节点电压iiU，使得由节点电压求得的功率iP、iQ与指定的节点注入功率

28、iP、iQ相等，或者说，使失配功率（mismatch power）P、Q满足给定的精度要求。此时三类节点的处理方法为：对V 节点，因其电压已给定，仍不参与迭代；对PQ节点，因其P和Q指定，U和待求，故既有有功失配功率，也有无功失配功率，即每一个PQ节点有两个迭代方程，并需设定电压的初值()0iU和()0i，mi,1L=；对PV节点，因其P和U指定，Q和待求，故仅有iP一个迭代方程，并需设定无功功率初值()0iQ 和电压相位初值()0i，1,1+=nmiL。每次迭代后，对 PV 节点，令()()kiikiUU=&，计算无功功 率()()ijijijijjjikiBGUUQcossin=，检 验

29、其 是 否 满 足 约 束 条 件 miniQiQmaxiQ：如不满足，“越限代限”，意思是“越过下限时代之以下限，越过上限时代之以上限”，此时PV 节点便转换成PQ 节点，转入下一次迭代。综上可知，利用 N-R 迭代求解混合坐标形式的潮流方程时共有（1n）个有功失配功率方程（1,1=niL）和m个无功失配功率方程（mi,1L=），方程总数为（1+mn），未知量有（1n）个电压相角i（1,1=niL）和m个电压幅值iU（mi,1L=），总数为（1+mn），方程数与未知量数相等，方程有定解。此时的迭代方程为+mmmnmmmnmnnnnnmnmnUQUQQQUQUQQQUPUPPPUPUPPPQQ

30、PPLLMMMMLLLLMMMMLLMMKKK111111111111111111111111111mnUUMMKKK111=0 ()()()()()()()()()()()()+=+kmkknkkmkknkkmkknkUUUUUUMMMMMMKKKKKKKK1111111111111 L2,1,0=k （324）简记为 0=+ULKNHQP ()()()()()()+=+kkkkkkUUU11 L2,1,0=k （325）收敛判据为 iiQP,max。式中的雅可比矩阵为=LKNHJ ，其中H为（1n）()1 n 矩阵，其各元素的表达式为 ()()=+=ijijijijjijiijiijiii

31、ijijijijjiiiiiBGUUPHQBUBGUUPHcossincossin2 （326）式中iQ的定义见式（323），为由节点电压求得的无功，此处应用它是因在计算失配功率iQ时已算出，直接引用可节省工作量；N为（1n）m矩阵，其各元素的表达式为 ()()+=+=ijijijijijiijiiijiiiijijijijjiiiiiiiBGUUPNUPGUBGUGUUPNsincos/sincos2 （327）式中iP的定义见式（323），为由节点电压求得的有功功率；K为m()1 n矩阵，其各元素为()()+=+=ijijijijjijiijiijiiiijijijijjiiiiiBGUU

32、QKPGUBGUUQKsincossincos2 （328）L为mm 阶 矩阵，各元素为()()=+=ijijijijijiijiiijiiiiiiijijijijjiiiiBGUUQLUQBUBUBGUUQLcossin/2cossin （329）观察上述各元素表达式发现：H、N、L、KH、N、L、K的非对角元 ijH、ijN、ijK、ijL的表达式均只有一项，且均含ijG、ijB。可以想见，如节点i和j间无支路直接连接，则 ijG、ijB为0，从而对应的ijH、ijN、ijK和ijL 均为 0。所以雅可比矩阵J是稀疏阵。利用稀疏技术可以节省计算机内存及提高计算速度。同时，J J具有强对角性

33、，但不对称。由于J J阵中的元素随 和U而改变，因而利用 N-R 迭代求解潮流方程时每次均需形成J J阵，每次要解修正方程，因而运算量大，但其收敛速度快，一般迭代 57 次便可得到满意的精度。且迭代次数不随节点数n明显增加。利用 N-R 迭代求解潮流的计算流程示于图 34。图 34 N-R 潮流迭代框图 【例例 33】利用 N-R 迭代求解【例 28】系统的潮流，设发电机 G1 的端电压为 1 p.u.，发出的有功、无功可调；发电机 G2 的端电压为 1 p.u ，按指定的有功 5.0=P p.u 发电，取 410=。节 点 编 号 输入原始数据 形成节点导纳矩阵 YB 设定初值，k=0 计算

34、失配功率 Pi、Qi iiQP,max。第六步，形成雅可比矩阵（阶数为 77）()=3882.280000.151743.90000.07523.70000.57523.20000.156418.303879.110000.00000.55587.85587.31743.93879.112422.200000.07523.25587.33110.60000.00000.00000.09889.40000.09889.40000.07523.70000.57523.20000.01632.290000.151743.90000.55587.35587.39889.40000.153768.313

35、879.117523.25587.33110.60000.01743.93879.115622.200MMMLLLLLLLLMMMMJ第七步，解修正方程：得到 ()o484819.701=,()o840433.502=；()o575785.503=,()o0981331.004=；()003449.001=U，()028523.002=U，()033880.003=U。从而 ()()()o484819.7010111=+=，()()()o840433.5020212=+=；()()()o575785.5030313=+=，()()()o0981331.0040414=+=；()()()0034

36、49.1010111=+=UUU，()()()028523.1020212=+=UUU，()()()03388.1030313=+=UUU。然后转入下一次迭代。经三次迭代后，iiQP,max104。迭代过程中失配功率的变化情况列于表 31，节点电压的变化情况列于表 32。表 31 迭代过程中失配功率变化情况 k 0 1 2 3 P1 P2 P3 P4 Q1 Q2 Q3-0.8055-0.18 0 0.5-3.3720 0.2475 0.3875 1.932210-2 4.004810-3-5.507610-3-1.340110-2-1.484810-2-3.857410-2-4.244010-

37、2 2.0010-4-4.7310-5-8.6610-5-8.3710-5-2.7510-4-4.0610-4-4.3410-4-7.710-7 9.3910-7-1.0110-6 ijsin，又由于高压网络中电阻R ijG，故ijijBcos ijGijsin，于是可将ijGijsin略去，并取ijcos1。又因式（326）中iiH 表达式的第一项 iiiBU2远大于第二项iQ，式（329）中iiL表达式的第一项iiiBU远大于第二项 iiUQ/，故均可将第二项略去（其物理解释是：iiiBU2代表除 i点外其余节点均接地时节点i的注入无功，而iQ代表所有节点均不接地时节点i的注入无功，显然前

38、者远大于后者）。于是，H和L各元素的表达式成为 ()=+=ijjiijijijijjiijiiiiiiiiiBUUBGUUHBUQBUHcossin22 （331）()=ijiijijijijiijiiiiiiiiiiBUBGULBUUQBULcossin/（332）从而有 =1,1212,1211,1111,212222221121,11112211121nnnnnnnnnnnBUBUUBUUBUUBUBUUBUUBUUBUHLMMMMLL =12100nUUUO 1,11,11,111nnnnBBBBLMML UBUUUUn=12100O （333）=mmmmmmmmmBUBUBUBUBU

39、BUBUBUBULLMMMMLL212222221211121111=mUUU0021OBUBBBBmmmm =LMML1111 （334）将其代入式（330），得到=+=+00UBUQUBUP （335）又因 U近似为一单位阵，故可取BUB，并各乘以1U和 1 U，从而上式成为=+=+0/0/UBUQBUP （336）此即 P-Q 解耦迭代的修正方程，注意上式中B和B 的元素均直接取原节点导纳矩阵相应元素的虚部，但阶数不同：前者为)1()1(nn，后者为mm。同理，U为1)1(n，U 为1m。P-Q 解耦的迭代公式为()()()()()()()()=+kkkkkkkkUQBUUUPB1111

40、 L2,1,0=k （337）【例例 34】利用 P-Q 解耦迭代求解【例 33】。解解：此时B和B 分别为 =2493.507757.289889.4150093.3101743.93879.114022.20B=7757.28150093.311743.93879.114022.20B 设 初 值 同 例3 3，算 得 的 失 配 功 率 仍 为 ()8055.001=P，()18.002=P，()003=P，()5.004=P，()3720.001=Q，()2475.002=Q，()3875.003=Q，从而()()8055.00101=UP，()()18.00202=UP，()()00

41、303=UP，()()5.00404=UP，()()3720.00101=UQ，()()2475.00202=UQ，()()3875.00303=UQ。代入式（345），有 ()()()=0307.10267.10105.13875.02475.03720.01111131211BUUU 继续迭代，所得结果示于表 35，迭代收敛后的计算同 N-R，不再重复。()()()()=oooooooo5691.18767.63933.74811.95.0018.08055.00000114131211B表 35 迭代过程中节点电压变化情况 k 1 2 3 4 1U 2U 3U 1 2 3 4 5 6 7

42、 8-9.4811-7.2731-7.4879-7.4639-7.4802-7.4743-7.4746-7.4748-7.3933-5.5498-5.9017-5.8538-5.8571-5.8535-5.8548-5.8548-6.8767-5.2948-5.6451-5.5844-5.5875-5.5854-5.5865-5.864-1.5691-0.0328-0.2233-2.2010-0.2043-0.2013-0.2022-0.20221.0105 0.9862 0.9905 0.9917 0.9917 0.9915 0.9916 0.9916 1.0267 1.0150 1.017

43、2 1.0175 1.0176 1.0175 1.0175 1.0175 1.0307 1.0213 1.0225 1.0228 1.0229 1.0229 1.0229 1.0229 由此例可见：由于上述两点改进：解耦降阶和用常数阵B、B 取代H、L，从而使计算大为简化。P-Q 解耦迭代的次数多于 N-R 法，但每次迭代费时少，约为原算法时间的 1/3，故总的速度快于 N-R 法。特别值得指出的是：虽然采用了一些简化假设，但丝毫不影响最终结果的精度（请比较表 33 和表 35 中节点电压的幅值和相角），因收敛判据和失配功率的计算公式与 N-R 时完全相同。这是一种新的思维方法目标控制的方法：

44、虽每一步不如常规时准确，但如果省时，可多走几步，仍能先期达到同样的目标。还需说明的是：由于其推导过程中 采用了一些简化假设，如实际系统中这些 假设不成立，会出现潮流求解不收敛的情 况，需采取一定的补救措施。例如，对配 电线路，其电阻电抗比 R/X较大，不符 合假设 R X，会导致潮流不收敛。此 时可采用下述方法：在支路中增加一个虚 拟节点k，使得每条支路的R/X均足够小，如图 36 所示。这样处理后可使潮流收敛。关于H、L矩阵的上述简化只是各种简化方法中的一种，形式上最简单，最易于理解，即B、B 均取节点导纳矩阵元素的虚部组成，但其收敛性能并不好（如上例中迭代了 8 次），因而可作一些改进。如

45、形成B时不计线路的充电电容及变压器非标准变比的影响，因其与系统中有功功率的分布基本无关；此外，对支路电阻可采用不同的处理方法：B和B 均计及电阻，取()22ijijijijxrxb+=，称为 B-B型；B不计电阻，取ijijxb1=，称为 X-B 型；B 不计电阻，称为 B-X 型；B和B 均不计电阻，称为 X-X 型。实践表明，X-B 和 B-X 型收敛性能较好。同时，B-X 型对R/X较大的网络有更好的收敛性。这方面的研究工作仍在进行中，有兴趣者可参阅有关文献。P-Q解耦迭代时可在 P迭代求出()1+k后立即代入Q迭代，以加快收敛。本节首先从节点电压方程入手导出了潮流计算的数学模型潮流方程

46、，然后介绍了解非线性潮流方程的两种迭代方法：N-R 迭代和 P-Q 解耦迭代。前一种方法是求解非线性方程的通用方法，后一种方法是利用电力系统的特点得出的简化快速方法，由于其简单快速，因而得到了广泛应用。图 36 支路 R/X 太大时的处理 回顾电力系统潮流计算的发展史，其经历了几个阶段：开始是手算，只能求解简单系统的潮流，而且还需采用一些简化假设。其后在 20 世纪 50 年代曾广泛采用计算台（包括直流计算台和交流计算台）解算潮流。这实际上是一种物理模拟方法。限于其规模，只能进行小型电力系统的潮流解算。自 50 年代开始，随着数字计算机的出现，电力工作者立即采用这一强有力工具，进行了计算机求解

47、潮流问题的研究。最初，并没有取得多少成功，因为它仅仅是手工劳动的替代而没有充分发挥计算机的潜力。采用的数学模型是回路电流方程，需要大量的数据准备工作。50 年代中期，美国学者 Ward 和 Hale 应用节点形式的数学模型取得了成功，基于节点导纳矩阵的迭代算法，即 G-S 迭代问世。到 60 年代，电力系统的规模越来越大，G-S 迭代的缺点日益突出；收敛性差，速度慢，于是出现了牛顿法和基于节点阻抗矩阵的算法。此后，潮流计算便沿着这两条路线向前发展：在 N-R 迭代的基础上于 70 年代提出了 P-Q 解耦迭代，后又在直角坐标 N-R 基础上提出了带二阶项的牛顿法；在节点阻抗矩阵法的基础上采用了

48、分块矩阵和分割系统的方法。直到今天，潮流计算仍在继续研究中，新的算法不断出现。其发展方向不外是改善收敛、提高速度和减少内存。同时，新的潮流问题，如动态潮流、谐波潮流、状态估计潮流、概率潮流、最优潮流等，也得到了研究。此处仅对基本的常规潮流，也称静态潮流，作了介绍。有了这个基础，其它潮流问题便可进一步深入。第三节第三节 潮流计算中的有关技术潮流计算中的有关技术 计算潮流时，除了上面介绍的基本原理和方法外，还要用到一些有关技术，以加快计算速度、节省计算机内存和提高求解的稳定性。这方面的技术具有重要的实用价值，内容十分丰富，而且在不断发展中。此处仅摘其要者介绍之，作为潮流问题的扩展和延伸。内容包括：

49、潮流方程迭代求解时的初值设定、线性方程组的求解、处理稀疏矩阵的稀疏技术和网络化简。下面分别予以介绍。一、潮流方程迭代求解时的初值设定一、潮流方程迭代求解时的初值设定 从计算机求解潮流的过程可看出，无论采用何种迭代方法，都需首先设定初值，然后方能进入迭代。在前述例题中，均设电压幅值为 1 p.u ，相角为o0，称为平起动。实际计算中，这样的初值设定一般可收敛，但不难想象，如果有更好的初值设定，一定可以减少迭代次数，加快计算进程，提高求解的稳定性。下面介绍这类方法中的一种。这种方法由两步组成：第一步设定电压相角的初值，第二步设定节点电压幅值的初值。第 一 步 中 采 用 如 下 两 点 假 设：节

50、 点 电 压 的 相 角 差ij很 小，从 而 可 取 ijijijsin,1cos；取所有节点电压的幅值均为 1。于是由混合坐标形式的潮流方程（316）)sincos(1ijijijijnjjiiBGUUP+=得到 1,1)(=+=niBGPjjiijijiK (338)这是一组线性方程，共（1n）个，可解出（1n）个未知量 121,nL。第二步中只需计算 PQ 节点电压幅值的初值，因 PV 节点和 V 节点的电压幅值已给定。对 PQ 节点，有 miUYUjQPSjjijiiii,1 L&=+=(339)即 =+=+jijijijiiiiiiiijUjBGUjjQPUjQP)sin(cos)