电力系统潮流计算

第四章电力系统潮流的计算机算法第四章电力系统潮流的计算机算法第一节电力网络的数学模型第二节等值变压器模型及其应用第三节节点导纳矩阵的形成和修改第四节功率方程和变量及节点分类第五节高斯-赛德尔法潮流计算第六节牛顿-拉夫逊法潮流计算第七节P-Q分解法潮流计算第四章电力系统潮流的计算机算法第三章讨论简单电力网络的潮流分布计算，理解了与之相关的各种物理现象。对于复杂电力网络的潮流计算，一般必须借助电子计算机进行。运用电子计算机，一般要完成以下步骤：1、建立电力网络的数学模型2、确定解算方法3、制定计算流程和编制计算程序本章将着重讨论前两项，主要阐述在电力系统潮流的实际计算中常用的、基本的方法。第四章电力系统潮流的计算机算法第一节电力网络的数学模型电力网络的数学模型指的是将网络有关参数及其相互关系归纳起来，组成可以反映网络性能的数学方程式组。也就是对电力系统的运行状态、变量和网络参数之间相互关系的一种数学描述。有：节点电压方程回路电流方程割集电压方程等节点电压方程又分为以节点导纳矩阵表示的节点电压方程和以节点阻抗矩阵表示的节点电压方程。第四章电力系统潮流的计算机算法一、节点导纳矩阵的节点电压方程在电路理论中，已经讲过了节点导纳矩阵的节点电压方程对于n个节点的网络其展开为上式中，是节点注入电流的列向量。是节点电压的列向量。是一个n×n阶节点导纳矩阵。BIBUBY第四章电力系统潮流的计算机算法以网络节点导纳矩阵表示的节点电压方程在进行潮流计算时,可以减少计算机的内存,提高运算速度,因此是最为常用的。节点导纳矩阵各元素的物理意义：对角元：iiYijUiiiijUIY0|)/(非对角元：ijUijjijUIY0|)/(jiY第四章电力系统潮流的计算机算法例：如图所示网络以对角元和非对角元、为例。22Y12Y12y1231I2I10y20y30y23y31y13220y12y23y2I1I+-2U3I32Y第四章电力系统潮流的计算机算法二、节点阻抗矩阵的节点电压方程由的两边都左乘，可得，而，则节点电压方程为BBBUYI1BYBBBUIY1BBZY1BBBUIZ第四章电力系统潮流的计算机算法第二节等值变压器模型及其应用一、变压器为非标准变比时的修正无论采用有名制或标么制，凡涉及多电压级网络的计算，在精确计算时都必须将网络中所有参数和变量按实际变比归算到同一电压等级。实际上，在电力系统计算中总是有些变压器的实际变比不等于变压器两侧所选电压基准值之比，也就是不等于标准变比，而且变压器的变比在运行中是可以改变的。这将使每改变一次变比都要重新计算元件参数，很不方便。下面将介绍另一种可等值地体现变压器电压变换功能的模型。二、等值变压器模型第四章电力系统潮流的计算机算法1U2UZ12K:11I2ITZZ图一等值双绕组变压器1U2U1I2IKYT/KYKT)1(2)1(KKTY121I2ITKZ1KKZTKZKT121U2U12（a）非标准变比时的修正电路（c）以变压器阻抗表示（b）以变压器导纳表示第四章电力系统潮流的计算机算法对修正电路，其二端口网络的传输参数方程为2T211IKZUKU21IKI++--1U2U1I2I0Z1Z2Z对以阻抗表示的π型等值电路20201)1(IZUZZU210220121)1()]1(11[IZZUZZZZI比较，可得T0KZZT11ZKKZT221ZKKZ第四章电力系统潮流的计算机算法三、等值变压器模型的应用12341:K1K2:1T1T2如图所示网络，变压器T1、T2的电抗归算至1和4侧时，可用下图电路来模拟。1TZ1:1KlZ2lY2lY1:2K2TZ第四章电力系统潮流的计算机算法最终的等值电路如下图。lZ2lY2lY11TZK22TZK1111KZKT11211KZKT22221KZKT1222KZKT第四章电力系统潮流的计算机算法例等值电路如图，给出支路阻抗和对地导纳标么值、变压器变比标么值。求节点导纳矩阵。12345015.0j05.1:125.0j25.0j3.008.0j25.0j25.0j25.004.0j35.01.0j03.0j1:05.1T1T2第四章电力系统潮流的计算机算法第三节节点导纳矩阵的形成和修改一、节点导纳矩阵的形成节点导纳矩阵的计算归纳总结如下：1、节点导纳矩阵的阶数等于电力网络中除参考节点（一般为大地）以外的节点数。2、节点导纳矩阵是稀疏矩阵，其各行非对角非零元素的个数等于对应节点所连的不接地支路数。3、节点导纳矩阵的对角元素，即各节点的自导纳等于相应节点所连支路的导纳之和，即第四章电力系统潮流的计算机算法4、节点导纳矩阵的非对角元素等于节点和间支路导纳的负值，即5、节点导纳矩阵是对称方阵，因此一般只需要求取这个矩阵的上三角或下三角部分。ijYji节点导纳矩阵形成的基本方法：添加支路法。第四章电力系统潮流的计算机算法ijNijZijZijiijjNNNijZijZijZ二、节点导纳矩阵的修改在电力系统计算中，对于已知网络，其节点导纳矩阵已经形成。如果网络接线发生局部变化，此时不必重新计算节点导纳矩阵。仅仅需要在原节点导纳矩阵的基础上进行必要的局部修改就可以得到所求节点导纳矩阵。下面介绍几种情况。图三电力网络接线变更示意图(a)(b)(c)(d)第四章电力系统潮流的计算机算法（1）从原有网络中引出一条新的支路，图三（a）。同时增加一个新的节点。新增加节点的对角元素为：新增加非对角元素为：原有节点的自导纳增量为：第四章电力系统潮流的计算机算法（2）在原有节点和间增加一条支路，图三（b）。此情况下节点导纳矩阵的阶数不变。有关元素修改如下：ij第四章电力系统潮流的计算机算法（3）在原有节点间切除一条阻抗为的支路，见图三（c）这种情况下，相当于在节点和间增加阻抗为的支路，此时，节点导纳矩阵的阶数不变，其元素修正如下：jiijzijz第四章电力系统潮流的计算机算法（4）原有网络节点和之间支路阻抗由改变为，这种情况下，可以看作是在节点和间切除阻抗为的支路，并在节点和间增加阻抗为的支路，如图三（d）。此时，节点导纳矩阵的阶数不变，其元素修正如下：ijijzijzijijzijijz第四章电力系统潮流的计算机算法（5）原有网络节点和之间变压器的变比由变为时，相当于在原网络节点和之间切除一变比为的变压器支路，而又增加一个变比为的变压器支路。其元素修正如下：ijkkjikk第四章电力系统潮流的计算机算法第四节功率方程和变量及节点分类一、功率方程每个节点的注入功率方程式为：其中：对于N个节点的电力网络，可以列出2N个功率方程。每个节点具有四个变量，N个节点有4N个变量，但只有2N个关系方程式。第四章电力系统潮流的计算机算法二节点电力系统121U2U111LGIII222LGIII111LGSSS222LGSSS12y10y20y12111GGGjQPS222GGGjQPS1U2U111LLLjQPS222LLLjQPS第四章电力系统潮流的计算机算法节点电压方程：YUI，对上述二节点系统，即212221121121UUYYYYII展开112121111USUYUYI222221212USUYUYI即212111111UYUUYUS222212122UYUUYUS令111jeUU222jeUU，，并考虑到111jQPS222jQPS将有功功率，及无功功率，分列1P2P1Q2Q第四章电力系统潮流的计算机算法),,,(),(212111111UUPUPPPPLG),,,(),(212122222UUPUPPPPLG),,,(),(212111111UUQUQQQQLG),,,(),(212122222UUQUQQQQLG第四章电力系统潮流的计算机算法二、变量的分类1、负荷消耗的有功、无功功率（、）取决于用户，因而是无法控制的，故称为不可控变量或扰动变量。一般以列向量表示，即2、电源发出的有功、无功功率（、）是可以控制的变量，故称为控制变量，以列向量表示，即3、母线或节点电压和相位角（、），是受控制变量控制的因变量。其中主要受的控制，主要受的控制。故、称为系统的状态变量，以列向量表示，即LPLQGPGQduUUGQGPUxTLnLLnLQQPPd11TGnGGnGQQPPu11TnnUUx11第四章电力系统潮流的计算机算法三、节点的分类1、PQ节点：已知P、Q负荷、过渡节点，PQ给定的发电机节点，为大部分节点2、PV节点：已知P、V给定PV的发电机节点，具有可调电源的变电所，为少量节点3、平衡节点：也称为松弛节点，摇摆节点PQ节点12345PQ节点PV节点PQ节点平衡节点第四章电力系统潮流的计算机算法第五节高斯-塞德尔法潮流计算迭代法的基本思路考察下列形式的方程：这种方程是隐式的，因而不能直接得出它的根，但如果给出根的某个猜测值，代入上式的右端，即可求得：再进一步得到：第四章电力系统潮流的计算机算法如此反复迭代：确定数列{xk}有极限则称迭代过程收敛，极限值x*为方程的根。上述迭代法是一种逐次逼近迭代法，称为高斯迭代法。第四章电力系统潮流的计算机算法高斯-塞德尔迭代法在高斯法的每一次迭代过程中是用上一次迭代的全部分量来计算本次的所有分量，显然在计算第i个分量时，已经计算出来的最新分量并没有被利用，从直观上看，最新计算出来的分量可能比旧的分量要好些。因此，对这些最新计算出来的第k+1次近似分量加以利用，就是高斯-塞德尔迭代法。高斯-塞德尔迭代法计算潮流功率方程的特点：描述电力系统功率与电压关系的方程式是一组关于电压的非线性代数方程式，不能用解析法直接求解。第四章电力系统潮流的计算机算法高斯或高斯-塞德尔法的迭代方程。假设有n个节点的电力系统，没有PV节点，平衡节点编号为s。可列出节点电压方程nininnninniniiiininiIIIIUUUUYYYYYYYYYYYYYYYY21212121222212111211第四章电力系统潮流的计算机算法由iiiiIUjQPS得iiiiUjQPI对第i个节点iiininiiiiiUjQPUYUYUYUY2211即nijjiiijijiiiUjQPUYUY1第四章电力系统潮流的计算机算法可得以下迭代方程：对每一个PQ节点都可列出一个方程式，因而有n-1个方程式。在这些方程式中，注入功率Pi和Qi都是给定的，平衡节点电压也是已知的，因而只有n-1个节点的电压为未知量，从而有可能求得唯一解。第四章电力系统潮流的计算机算法设除平衡节点以外的其他节点的电压初值)0()0()0(jiieUUsini,2,1则高斯法潮流计算的迭代格式，即各节点电压的第k+1次计算值为nijjkjijkiiiiikiUYUjQPYU1)()()1(1sini,2,1在迭代过程中，平衡节点（编号为s）电压始终不变。第四章电力系统潮流的计算机算法高斯-塞德尔迭代法解潮流如下：如果系统存在PV节点，则对于PV节点（编号为p),迭代前应假设其注入无功功率和节点电压相位初值,)0(pQ)0(p完成第k次迭代后，做第k+1次迭代前，先求出节点p的注入无功功率PV节点的处理第四章电力系统潮流的计算机算法然后将其代入下式，求出