简介几种潮流计算

简介几种潮流计算电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。

同时，为了实时监控电力系统的运行状态，也需要进行大量而快速的潮流计算。

因此，潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。

在系统规划设计和安排系统的运行方式时，采用离线潮流计算；在电力系统运行状态的实时监控中，则采用在线潮流计算，下面简单介绍三种潮流计算方法。

一、基于多口逆向矩阵的并行潮流计算方法多口逆向矩阵方法是求解线性方程组的普通并行方法，它只是修改了串行方法的几个部分，并且非常适用于从串行到并行的编程。

该方法已用于一些电力系统并行分析方法，比如说机电暂态稳定分析和小信号稳定性，并且并行效率高。

基于多口逆向矩阵方法，本文提出了一种并行牛顿潮流算法。

对一个划分几个网络的大型互联系统模型的仿真结果表明这种并行算法是正确的并且效率很高。

关键词:并行潮流计算，串行潮流计算，多口逆向矩阵方法，线性方程组，电力系统分析随着电力系统规模的扩大，尤其是区域互联网络，人们要求速度更快效率更高的功率计算，传统的串行计算越来越难满足要求，特别是对实时控制。

作为电力系统的基本计算，它的效率的提高会使其他为基础的计算速度都得到提高。

因为传统串行计算变的越来越难满足要求，并行计算成为提高潮流计算效率的需要。

潮流计算的主要步骤是求解稀疏线性方程组，因此对并行方法的研究主要集中在线性方程组的并行求解。

根据不同的实现方案，并行算法分为多因子方法、稀疏向量方法等等。

多口逆向矩阵方法在各种问题中是一种求解线性方程组的通用方法。

在这篇论文中，通过最常见的电力系统中的节点电压方程来说明这种方法。

多口逆向矩阵法不需要在矩阵中集中调整边界点，我们根据子网的密度把矩阵分裂并且把边界节点集中在顶部，整个网络的节点电压方程组如下：消去上矩阵中对应子网的部分，只保留边界部分。

经过网络分割，边界矩阵TT Y 注入电流向量T I 被分为主控制网和各个子网。

设定主控制网矩阵为TT Y ,子网i 的为TTiY 。

注入电流矩阵分割为子网i 为TiI ，即()310∑=+=ki TTiTT TT Y Y Y()41∑==ki TiT II定义：()5~~11iii Ti Ti Ti iT ii Ti TTi TTi I Y Y I I Y Y Y Y Y ---=-= TTi Y ~是子网i 的边界等效导纳矩阵，TiI ~是子网i 的边界等效注入电流矩阵。

()6~~110∑∑===⎪⎭⎫ ⎝⎛+ki TiT k i TTi TT I U Y Y去除与其他子网的接线，子网i 的网络方程如下：()7⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡''⎥⎦⎤⎢⎣⎡i Ti i Tiii iTTi TTi I I U U Y Y Y YTiU '是网络i 的边界电压向量，i U '是网络i 的内部节点电压向量，这是在去除与其他子网连接线情况下的解，TTiY ~和TiI ~的计算公式如下：()8~11--=-=TTiiT ii Ti TTi TTi Z Y Y Y Y Y()9~11TiTTi i ii Ti Ti Ti U Z I Y Y I I '=-=--TTiZ 是网络i 的导纳矩阵的边界节点子矩阵，TTi Z的维数是边界节点的总数，它通常较小。

当求解TTiZ 时，我们不需要对子网矩阵求逆，通过设置的右边网络方程的关联节点设为1，可以求解线性方程组得到TTiZ ，运用向量稀疏技术，计算量可以进一步减少。

TiU '可以通过求解方程组7得到，同样也可以运用稀疏技术减少计算量。

用6式解得边界节点电压矩阵T U ，把T U 带回到各个子网络，使用补偿注入电流方法计算各电网电压。

假设边界节点电压为TiU ，方程右侧的补偿电流可有下式计算得到：()10~~TiTi TTi i I U Y T -=∆网络的节点电压由下式计算得到：()11⎥⎦⎤⎢⎣⎡∆+=⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡i Ti Ti i Ti ii iT Ti TTi I I I U U Y Y Y Y网络矩阵方程和串行计算是一样，其他子网的影响由Ti I ∆表示，由此矩阵计算子网非常方便。

多口逆向矩阵方法可以用来求解大规模线性代数方程，它只是改变了串行步骤的一小部分。

这种方法只是增加了一点工作量，通信时间短，并行的效率高。

基于多口逆向矩阵方法的并行牛顿潮流计算当用多口逆向矩阵方法求解并行线性方程组时，对方程组有个基本的要求：各个子网的方程两侧之和再分解前后是一样的。

通过解列母线和线路，并使用串行计算步骤，导纳矩阵和注入电流向量满足（3）、（4）的要求，这样保证了节点电压方程组可以并行求解。

对于功率修正方程，左边和右边不相等，对应边界节点不满足方程（3）、（4），需要先进行处理。

对于pq 类型的边界点，节点电压是给定的，满足：()12,1)(22,1)(11∑∑====nk k ii nk k ii F F F F()13,1)(∑==nk k ijij J J它可不需要修改应用于基于多口逆向矩阵法的并行计算。

)(1k i F 是子网k 的有功功率偏差，)(2k i F 是子网k 的无功功率偏差，)(k ij J 是子网k 的雅克比矩阵的元素。

对于pv 类型的边界节点，有功功率偏差i F 1和pq 节点的一样不需要修正，i F 2是电压修正不满足∑==nk k ii FF ,1)(22，对应于i F 2的雅克比元素不满足（13）。

在相同电压下，每个子网的pv 类型的电压修正和对应的雅克比元素和串行计算时一样。

因此为了满足（12）、（13）需要把电压修正为原来的1/n 。

对于θV 节点，i F 1和i F 2是电压的实部和虚部，可以使用与pq 节点相似的方法进行处理。

使用以上方案形成各子网的修正方程，可以使用多口逆向矩阵方法进行并行潮流计算，主要步骤如下:（1）设置节点初始电压和类型，迭代此次n=1，迭代收敛误差ε，最大迭代次数N 。

（2）根据节点类型用串行计算步骤同时形成各子网的功率修正方程。

（3）对于pv 和θV 节点，修正雅克比矩阵元素并将修正方程右侧乘以1/n 。

（4）使用多口逆向矩阵方法进行并行计算，得到电压修正值。

（5）找出各子网绝对值最大的电压修正值max V ∆。

（6）在各子网中如果max V ∆小于ε，则该子网迭代收敛，转到（7），如果n N ，则迭代不收敛，转到（8），如果n ≤N ，转到（2）。

（7）输出当前结果。

（8）结束计算。

二、孤网运行的潮流计算方法电力系统实际运行中，某些原因会导致局部电网解列，形成孤网运行的情况，为保证孤网的安全稳定运行，需要知道孤网的潮流分布。

然而，目前常规的潮流模型无法反映孤网运行的重要特征，得出的潮流分布与实际运行情况不符。

为此，文章建立了孤网运行潮流计算的数学模型。

新模型不需设置平衡节点，负荷变化由装设有调速系统和励磁系统的发电机组按照开环设置的参数共同承担。

通过模拟中国电力科学研究院 22 节点测试系统解列为 2 个孤网运行时的潮流分布，验证了新模型的有效性和实用性。

常规潮流计算的基本假定是：系统中所有控制已按所需的情况完成，即发电机组的调速系统、励磁系统的一、二次调节已按所需的输出条件完成了其设置与调节。

因此，发电机成为了 PV 或 PQ 节点，并将一台发电机设置为给定电压、注入功率松弛的平衡节点形式，以满足网络方程中的有功、无功损耗特性。

这样，就可方便的讨论电力系统规划设计中的问题。

孤岛电网潮流计算数学模型的建立及求解数学模型的建立。

电力系统的运行状态，可用如下的微分代数方程组来描述x = f ( x , y , u) (1) g ( x , y ) = 0 (2)式中：变量 x 为节点电压，y 为节点功率，u 为同步发电机组调速系统和励磁系统的空载频率 0f 、空载电压 0VG 。

式(1)为发电机组调节的动态方程，表征 了同步发电机组的动态过程；式(2)为网络方程，表征了网络的运行约束。

在完成一次调节后，系统将稳定于稳态运行点，这时描述发电机组动态性能的微分方程式(1)退化为代数方程式(3)，此时，发电机组动态方程达到平衡点，即发电机组的有功频率、无功电压静态特性曲线与负荷功率特性曲线相交。

x = f ( x , y , u) = 0 (3)如果设全系统共有 g 台同步发电机组，并设其有功频率静态特性的斜率为 KGi 。

第 i 台同步发电机自动调速系统的有功频率静态特性可表示为()()40i Gi Gi f f K P -'-=即()5/0GiGi i K P f f -='式中：i f 0、Gi P 分别表示第 i 台发电机的空载频率和输出的有功功率。

因稳态运行时，各台机组的频率应相同，所以有.......3,2,1//00=-=-i K P f K P f GgGg g Gi Gi i 式中：g f 0、Gg P 和 Gg K 分别表示第g 台发电机组的空载频率、有功输出和功频静态特性的斜率。

发电机无功功率的励磁调节可表示为gi QGi i i G K I V V +=0。

在全系统频率 f ′ 相同的稳态运行点，同步发电机自动调速系统需满足的方程为1g 2,10//00-==+--、、i K P K P f f Gg Gg Gi Gi g i励磁系统调节需满足的方程为g i V Q K V i Gi gi i 、、,2,10/V G0i ==--节点有功平衡需满足的方程为n ,2,10、、==--i P P P i Di Gi节点无功平衡需满足的方程为n ,2,10Q Q Q 、、==--i i Di Gi式中：Di P 、Di Q 、i P 、i Q 分别表示节点 i 的有功负荷、无功负荷、有功注入功率、无功注入功率，其具体列式如下：()[]s i i i pi i i pi pi Di P V V C V V B A P 0200//++=()[]si i i qi i i qi qi DiQ V V C V V B AQ 0200//++=()ij ij ij ij j i i B G V V P θθsin cos +=∑()ij ij ij ij j i i B G V V Q θθcos sin -=∑式中：下标“0”表示初始运行条件时相关变量的值；系数 Pi A 、pi B 、pi C 和 qi A 、qi B 、qi C 表示负荷模型静态电压特性的参数；ij G 和 ij B 分别为节点导纳矩阵元素的实部和虚部；ij θ为节点 i 和节点 j 的电压相角差。

求解方法建立的新模型是非线性代数方程组的形式，目前，牛顿法是求解非线性代数方程组的有效算法，并且牛顿法在电力系统潮流计算中得到了广泛的应 用,故而新模型可以采用牛顿法求解。