潮流计算的基本算法及使用方法

潮流计算的基本算法及

使用方法

Company number：【00WT-88YT-W8CCB-BUUT-202108】

潮流计算的基本算法及使用方法

一、 潮流计算的基本算法

1.

牛顿－拉夫逊法

1．1 概述

牛顿－拉夫逊法是目前求解非线性方程最好的一种方法。这种方法的特点就是把

对非线性方程的求解过程变成反复对相应的线性方程求解的过程，通常称为逐次线性化过程，就是牛顿－拉夫逊法的核心。

牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发，沿着该点的

一阶偏导数——雅可比矩阵，朝减小方程的残差的方向前进一步，在新的点上再计算残差和雅可矩阵继续前进，重复这一过程直到残差达到收敛标准，即得到了非线性方程组的解。因为越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性。而所谓“某一邻域”是指雅可比方向均指向解的范围，否则可能走向非线性函数的其它极值点，一般来说潮流由平电压即各母线电压(相角为0，幅值为1)启动即在此邻域内。

1．2 一般概念

对于非线性代数方程组

即 fix1,x2,,xn0 i1,2,n (1－1)

在待求量x的某一个初始计算值x0附件，将上式展开泰勒级数并略去二阶及以

上的高阶项，得到如下的线性化的方程组

fx0fx0x00 (1－2)

上式称之为牛顿法的修正方程式。由此可以求得第一次迭代的修正量

x0fx0fx (1－3)

10

将x0和x0相加，得到变量的第一次改进值x1。接着再从x1出发，重复上述

计算过程。因此从一定的初值x0出发，应用牛顿法求解的迭代格式为

fxkxkfxk (1－4)

xk1xkxk (1－5)

上两式中：fx是函数fx对于变量x的一阶偏导数矩阵，即雅可比矩阵J；k为迭代次数。

由式（1－4）和式子（1－5）可见，牛顿法的核心便是反复形成求解修正方程

式。牛顿法当初始估计值x0和方程的精确解足够接近时，收敛速度非常快，具有平方收敛特性。

1．3 潮流计算的修正方程

运用牛顿－拉夫逊法计算潮流分布时，首先要找出描述电力系统的非线性方程。

这里仍从节点电压方程入手，设电力系统导纳矩阵已知，则系统中某节点（i节点）电压方程为

从而得

SiUiYijUj

j1n•••n进而有 －6）

PijQiUiYijUj0

j1 （1

式（1－6）中，左边第一项为给定的节点注入功率，第二项为由节点电压求得的

节点注入功率。他们二者之差就是节点功率的不平衡量。现在有待解决的问题就是各节点功率的不平衡量都趋近于零时，各节点电压应具有的价值。

由此可见，如将式（1－6）作为牛顿－拉夫逊中的非线性函数FX0，其中节

点电压就相当于变量X。建立了这种对应关系，就可列出修正方程式，并迭代求解。但

由于节点电压可有两种表示方式——以直角做表或者极坐标表示，因而列出的迭代方程相应地也有两种，下面分别讨论。 1．3．1 直角坐标表示的修正方程

节点电压以直角坐标表示时，令Uieijfi、Ujejjfj，且将导纳矩阵中元素

••表示为YijGijjBij，则式（1－7）改变为

PijQieijfiGijjBijejjfj0

j1n （1－7）

再将实部和虚部分开，可得

PieiGijejBijfjfiGijfjBijej0j1  nQifiGijejBijfjeiGijfjBijej0j1n （1

－8）

这就是直角坐标下的功率方程。可见，一个节点列出了有功和无功两个方程。

对于PQ节点（i1，2,,m1），给定量为节点注入功率，记为Pi、Qi，则由

式（2－8）可得功率的不平衡量，作为非线性方程

PiPieiGijejBijfjfiGijfjBijejj1 nQiQifiGijejBijfjeiGijfjBijejj1n （1－9）

式中Pi、Qi——分别表示第i节点的有功功率的不平衡量和无功功率的不平衡量。

对于PV节点（im1,m2,,n），给定量为节点注入有功功率及电压数值，

记为Pi、Ui，因此，可以利用有功功率的不平衡量和电压的不平衡量表示出非线性方程，即有

PiPieiGijejBijfjfiGijfjBijejj1

Ui2Ui2ei2fi2n （1－10）

式中Ui为电压的不平衡量。

对于平衡节点（im），因为电压数值及相位角给定，所以USesjfS也确

•定，不需要参加迭代求节点电压。

因此，对于n个节点的系统只能列出2n1个方程，其中有功功率方程n1个，无功功率方程m1个，电压方程nm个。将式（1－9）、式（1－10） 非线性方程联立，称为n个节点系统的非线性方程组，且按泰勒级数在fi0、ei0（i1,2,,n,im）展开，并略去高次项，得到以矩阵形式表示的修正方程如下

H11P1QJ111P2H21Q2J21PpHp1U2Rp1pPHnn12UnRn1N11L11N21L21Np1Sp1Nn1Sn1H12J12H22J22Hp2Rp2Hn2Rn2N12L12N22L22Np2Sp2Nn2Sn2H1pJ1pH2pJ2pHppRppHnpRnpN1pL1pN2pL2pNppSppNnpSnpH1nJ1nH2nJ2nHpnRpnHnnRnnN1nf1eL1n1N2nf2L2ne2 （1－11） NpnfpSpnepNnnfnSnnen上式中雅可比矩阵的各个元素则分别为 将（1－11）写成缩写形式

PHQJU2RNffLJe （1－12） eS对雅可比矩阵各元素可做如下讨论：

当ji时，对于特定的j，只有该特定点的fi和ei是变量，于是雅可比矩阵中各非对角元素表示为

当ji时，雅可比矩阵中各对角元素的表示式为

由上述表达式可知，直角坐标的雅可比矩阵有以下特点：

1) 雅可比矩阵是2n1阶方阵，由于HijHji、NijNji等等，所以它是一个不对称的方阵。

2) 雅可比矩阵中诸元素是节点电压的函数，在迭代过程中随电压的变化而不断地改变。

3) 雅可比矩阵的非对角元素与节点导纳矩阵YB中对应的非对角元素有关，当YB中的Yij为零时，雅可比矩阵中相应的Hij、Nij、Jij、Lij也都为零，因此，雅可比矩阵也是一个稀疏矩阵。

1．3．2 极坐标表示的修正方程

在牛顿－拉夫逊计算中，选择功率方程PijQiUiYijUj0作为非线性函数方

j1•n程，把式中电压向量表示为极坐标形式 则节点功率方程变为 将上式分解成实部和虚部

这就是功率方程的极坐标形式，由此可得到描述电力系统的非线性方程。

对于PQ节点，给定了

PiPiUiUjGijcosijBijsinijj1 nQiQiUiUjGijsinijBijcosijj1ni1、2、m1

（1－

13）

对于PV节点，给定了Pi、Ui，而Qi未知，式（1－13）中Qi将失去作用，于

是PV节点仅保留Pi方程，以求得电压的相位角。

（1－14）

对于平衡节点，同样因为Us、s已知，不参加迭代计算。

将式（1－13）、式（1－14）联立，且按泰勒级数展开，并略去高次项后，得出

矩阵形式的修正方程

H11P1QJ111P2H21Q2J21PpHp1PnHn1N11L11N21L21Np1Nn1H12J12H21J21Hp2Hn2N12L12N21L21Np2Nn2H1pJ1pH2pJ2pHppHnpH1nL1nN2nL2nHpnHnn1U1U12U2U

2pn （1－15）

雅可比矩阵终，对PV节点，仍可写出两个方程的形式，但其中的元素以零元素

U代替，从而显示了雅可比矩阵的高度稀疏性。式中电压幅值的修正量采用U的形

式，并没有什么特殊意义，仅是为了雅可比矩阵中各元素具有相似的表达式。

雅可比矩阵的各元素如下 将式（1－15）写成缩写形式

PHQJNU LU

（1－16）

以上得到了两种坐标系下的修正方程，这是牛顿－拉夫逊潮流计算中需要反复迭

代求解的基本方程式。 2.

快速分解法

2．1 概述

快速分解法的基本思想是：把节点功率表示为电压向量的极坐标方程式，抓主要矛盾，以有功功率误差作为修正电压向量角度的依据，以无功功率误差作为修正电压幅值的依据，把有功功率和无功功率的迭代分开来进行。快速分解法根据电力系统实际运行状态的物理特点，对牛顿-拉夫逊法潮流计算的数学模型进行合理的简化。

2．2 基本公式

在交流高压电网中，输电线路的电抗要比电阻大得多，系统中母线有功功率的变化主要受电压相位的影响，无功功率的变化主要受母线电压幅值变化的影响。在修正方程式的系数矩阵中，偏导数

QPQP和的数值相对于偏导数和是相当小的，作δVVδ为简化的第一步，可以将方程式（2－1）中的子块N和K略去不计，即认为它们的元素都等于零。这样，n1m阶的方程式便分解为一个n1阶和一个m阶的方程式，即将式（2－1）简化为式（2－2）和式（2－3）。

PHNQKLV1V (2－1) DPHδ (2－2)

1QLVDV (2－3)

上述的简化大大地节省了计算机的内存和解题时间，但是矩阵H和L的元素都是节点电压幅值和相角差的函数，其数值在迭代过程中是不断变化的。因此，快速分解法潮流计算的第二个简化，也是最关键的一步简化就在于把系数矩阵H和L简化成在迭代过程中不变的常数对称矩阵。在一般情况下，线路两端电压的相角差是不大的（通常不超过10~20）因此可以认为

cosij1 ， GijsinijBij （2－4）

此外，与系统各节点无功功率相适应的导纳BLDi必远小于该节点自导纳的虚部，即

BLDiQi2QVBii B 或 iiii2Vi考虑到上面的关系，矩阵H和L的元素的表达式便被简化为

HijViVjBij （i，j=1，2，…，n-1） （2－5） LijViVjBij (i，j=1，2，…，m) （2－6）

V1B11V1VBV2211HVn1Bn1,1V1V1B12V2V2B22V2Vn1Bn1,2V2V1B1,n1Vn1V2B2,n1Vn1 （2－7）

Vn1Bn1,n1Vn1V1B11V1VBVL2211VmBm1V1V1B12V2V2B22V2VmBm2V2V1B1mVmV2B2mVm （2－8）

VmBmmVm将式（2－7）和式（2－8）分别代入式（2－2）和（2－3），便得到：

11用VD1和VD2分别左乘以上两式便得简化了的修正方程式，可展开写成：

P1VB111BP221V2PBn1,1n1Vn1Q1VB111BQ221V2QBm1mVmB12B22Bn1,2B1,n1V1δ1VδB2,n122 （2－9）

Bn1,n1Vn1δn1B12B22Bm2B1mV1VB2m2 （2－10） BmmVm式（2－9）和式（2－10）就是快速分解法潮流计算的修正方程式，其中系数矩阵都是由节点导纳矩阵的虚部构成，只是阶次不同，矩阵B为n1阶，不含平衡节点对应的行和列，矩阵B为m阶，不含平衡节点和PV节点对应的行和列。

PiPisPiPisViVjGijcosδijBijsinδij （2－11）

j1nnQiQisQiQisViVj(GijsinδijBijcosδij) （2－12）

j1

修正方程式（2－9）和（2－10）与功率误差方程式（2－11）和（2－12）构成了快速分解法迭代的基本计算公式。 2．3 快速分解法的特点

快速分解法与牛顿法潮流计算的主要差别表现在它们的修正方程上。快速分解法通过对电力系统具体特点的分析，对牛顿法修正方程式的雅克比矩阵进行了有效的简化和改进，得到式（2－9）、式（2－10）所示的修正方程式。这两组方程式和牛顿法的修正方程相比主要有三个特点：

a）快速分解法的修正方程式用两个n阶线性方程组代替了一个2n阶线方程组。 b）快速分解法的修正方程式中系数矩阵的所有元素在迭代过程中维持常数不变。 c）快速分解法的修正方程式中系数矩阵是对称矩阵。

这些特点在提高计算速度和减少内存方面的作用是很明显的：首先，因为修正方程式的系数矩阵是导纳矩阵的虚部，因此在迭代过程中不必像牛顿法那样每次都要重新计算雅克比矩阵，这样不仅减少了运算量，而且也大大简化了程序；其次，由于系数矩阵在迭代过程中维持不变，因此在求解修正方程式时，不必每次都对系数矩阵进行消去运算，只需要在进入迭代过程以前，将系数矩阵用三角分解形成因子表，然后反复利用因子表对不同的常数项P/V或Q/V进行消去和回代运算，就可以迅速求得修正量，从而显着提高了迭代速度；第三，由于对称矩阵三角分解后，其上三角矩阵和下三角矩阵有非常简单的关系，所以在计算机中可以只存储上三角矩阵或下三角矩阵，从而也进一步节约了内存。

快速分解法所采用的一系列简化假定只影响了修正方程的结构，也就是说只影响了迭代过程，但未影响最终结果。因为快速分解法和牛顿法都采用同样的数学模型，最后

计算功率误差和判断收敛条件都是严格按照精确公式进行的，所以快速分解法和牛顿法一样都可以达到很高的精确度。

为了改善快速分解法的收敛特性，修正方程的系数矩阵B与B一般并不简单的是电力系统导纳矩阵的虚部，下面讨论一下B与B的构成。

B与B的阶数是不同的，B为 n1阶，B低于n1阶。因为式（2－10）不包含

于PV节点有关的项，所以，如果系统有r个PV节点，则B应为nr1阶。式（2－9）以有功功率误差为依据修正电压向量的角度，式（2－10）以无功功率误差依据修正电压幅值。为了加速收敛，使它们能够更有效地进行修正，可以考虑在B中尽量去掉那些与有功功率及电压向量角度无关或影响较小的因素，而在B中尽量去掉与无功功率及电压幅值影响较小的因素。所以，我们以电力系统导纳矩阵的虚部作为B和B时，可以在B去掉充电电容和变压器变比的影响，在B中去掉输电线路电阻对B的影响。

B和B的非对角元素和对角元素可分别按式（2－13）和（2－14）计算：

Bijxijrx2ij2ij Biijixijrx2ij2ij （）

B11bi0 （） Bxijjixij式（2－13）中rij和xij分别为支路ij的电阻和感抗，式（2－14）中bi0为节点i接地支路的电纳。

快速分解法改变了牛顿法迭代公式的结构，因此就改变了迭代过程的收敛性。牛顿法在迭代开始时收敛得较慢，当收敛到一定程度后，它的收敛速度非常之快，而快速分解法几乎是按同一速度收敛的，快速分解法每次迭代的计算量很小，因此快速分解法的计算速度比牛顿法有明显的提高。

二、 潮流计算的使用方法

1. 初始方式准备

对任何潮流模拟操作计算，总是在某一个初始的运行方式上进行。这种初始方式可

以是状态估计提供的实时运行方式，也可以是以往保存的历史运行方式。 2.

调度操作模拟

在准备好的初始潮流断面上，可以继续修改方式，模拟预想的潮流运行方式，再进行详细的潮流分析。模拟操作包括：

1) 开关刀闸变位模拟 2) 发电机功率调整 3) 负荷功率设置 4) 发电机分接头设置 5) 线路停运、投入 6) 变压器停运、投入 7) 母线停运、投入 8) 厂站停运、投入 3.

运行参数维护

潮流计算参数画面上可以设置算法、收敛判据、迭代次数、单/多平衡机等运行参数。

平衡发电机是电气岛内的电压相角参考点，当采用“单平衡机”模式时，电网的不平衡功率（包括发电、负荷和网损）都将由设定平衡机吸收。当采用“多平衡机”模式时，电网的不平衡功率将由多台发电机负责平衡，多台发电机之间的不平衡功率分配方式包括容量、系数和平均三种方式。选择容量时将根据发电机的可调容量分配，选择系数时

根据人工设置的系数按比例分配，选择平均时则平均分配不平衡功率。在分配过程中，确保发电机的出力在最大出力和最小出力范围内。

发电机参数中可以设置发电机的调节特性，包括节点类型（平衡节点、PQ 节点、PV节点等），对于PV节点可以设定控制机端电压还是高压侧母线电压以及控制的目标电压值，对于按指定系数参与有功调节的机组可以设置比例系数。 4.

计算结果分析

潮流计算结束后，计算结果分析包括： 1) 潮流计算状态 2) 电气岛、迭代信息 3) 潮流计算结果 4) 设备越限和重载监视 5) 运行信息 5.

误差统计

在潮流模拟计算完成后，如果现场很快发生了模拟的动作，可以从统计每个测点模拟计算值和实际量测值相比的误差，并统计出全网平均误差，统计方法遵循实用化考核细则。在表格中全部列出所有测点的SCADA量测值、潮流模拟计算值、考核基准值以及测点误差等内容。

如果只关心部分厂站的误差情况，局部误差统计中选择好需要关心的厂站，局部误差统计程序会过滤掉没有选中的厂站，只列出选中厂站的所有设备的误差统计情况。

如果想要保存当前误差统计的运行断面，可以保存误差断面，将当前的断面以及误差统计结果一起保存起来。误差统计历史记录包含所有保存的误差断面，内容包括统计时间、平均误差、执行用户、值班主机、断面名称以及操作信息等内容。