潮流计算的快速分解法(可编辑修改word版)

潮流计算的快速分解法

摘要：本文采用快速分解法进行潮流计算，分析其基本理论，并使用 MATLAB 软件进行编程设计。最后运用实例进行验证。结果表明快速分解法具有较好的迭代速度。关键词：潮流计算快速分解法 MATLAB 编程，实例验证

1引言

潮流计算是电力系统分析最基本、最重要的计算，是电力系统运行、规划以及安全性、可靠性分析和优化的基础，也是各种电磁暂态和机电暂态分析的基础和出发点。潮流计算要求具有可靠的收敛性，占用内存少，计算速度快，调整和修改容易，使用灵活方便。各种算法的改进以及新算法的提出，很多都是为了使潮流计算能更好地满足计算要求。本文应用快速分解法进行潮流计算，并给出算例分析。

2潮流计算的快速分解法

研究表明，用牛顿-拉夫逊法计算潮流时，每次迭代都要重新形成雅可比矩阵，然后重新对它进行因子表分解并求解修正方程。为避免每次迭代重新形成雅可比矩阵及其因子表，人们研究用定雅可比矩阵取代随迭代过程不断变化的雅可比矩阵，这种方法叫定雅可比法。此外，人们还结合电力系统的物理特点，发展了各种版本的解耦潮流算法，20 世纪 70 年代初提出的快速分解法是这一阶段的主要研究成果。

关于快速分解潮流算法，有三项里程碑意义的研究成果。其一是 Stott 在1974年发现的 XB 型算法；其二是 Van Amerongen 在1989 年发现的 BX 型算法；

其三是Monticelli 等人在1990 年所作的关于快速分解潮流算法收敛机理的理论阐述。这些研究工作不仅是电力系统计算方面的典范，也揭示了这样一个事实：工程上有效的方法一定有其深刻的理论来支持。

2.1快速分解法的修正方程及迭代格式

将极坐标型定雅可比法的修正公式重写如下：

H L

- ⎡ B H - G N ⎤⎡V ∆⎤ = ⎡∆P V ⎤ ⎢G B ⎥⎢ ∆V ⎥ ⎢∆Q V ⎥

（2.1）

⎣ M L ⎦⎣ ⎦ ⎣

⎦ 经验表明，电力系统中有功功率主要受电压相角的影响，而无功功率主要受电压幅值的影响，同时由于高压电网大部分线路的电阻比电抗小，因此在牛顿- 拉夫逊迭代中可以忽略雅可比矩阵的非对角块，即将G N ， G M 设为零，从而实 现有功和无功潮流修正方程的解耦。Stott 通过大量的计算实践发现，为了获得最好的收敛性，还要对雅可比矩阵的对角块作特殊的常数化处理：对系数矩阵 B H

，忽略支路电阻和接地支路的影响，即用- 1 x 为支路电纳建立的节点电纳矩阵 B '

代替 B ；对系数矩阵 B ，用节点导纳矩阵中不包含 PV 节点的虚部 B '' 代替；

V ∆

前的电压幅值用标幺值 1 代替。于是可得简化的修正方程式如下：

- B '∆

= ∆P V

（2.2）

- B ''∆V = ∆Q V

（2.3）

在潮流计算中，上述两个修正方程式依次交替迭代，Stott 把在此基础发展起来的潮流算法称为快速分解法（fast decoupled load flow ）。假定当前点为 (

(k ) ,V (k ) ) ，则求解((k +1) ,V (k +1) ) 的连续迭代格式如下：

⎧∆V (k ) = - B ''-1∆Q ((k ) ,V (k ) ) V (k )

（2.4）

⎨ (k +1)

⎩V

= V (k ) + ∆V (k ) ⎧∆(k ) = - B '-1

∆P ((k )

,V (k +1) ) V (k +1)

⎨ (k +1) (k ) (k )

（2.5）

⎩ = + ∆

快速分解法公式的特点是：① P -

和Q - V 迭代分别交替进行；②功率偏

差计算时使用最近修正过得电压值，且有功无功偏差都用电压幅值去除；③ B ''

和B ' 的构成不同，B '应用-1x 建立，并忽略所有接地支路（对非标准变比变压

器支路，变比可取为 1），而B ''就是导纳矩阵的虚部，不包括PV 节点。在快速

分解法的实施中，这些技术细节缺一不可，否则程序的收敛性将受到影响。

1989 年，荷兰学者 Van Amerongen 通过大量仿真计算发现了另一版本的快速分解潮流算法，他把该算法称为 BX 型算法，而把 Stott 的算法称为 XB 型算法，用以区分二者。BX 型算法与 XB 型算法的主要不同在于雅可比矩阵对角块的形成

进行简化时，保留了支路电阻的上。BX 型算法的处理方式是：在对系数矩阵B

H

影响，但忽略了接地支路项。BX 型算法的迭代格式与 XB 型算法是相同的。计算

经验表明，BX 型和 XB 型两种快速分解潮流算法在大部分情况下性能接近，在某

些情况下 BX 型算法收敛性略好。

快速分解法只对雅可比矩阵作了简化，但节点功率偏差量的计算及收敛条

件仍是严格的，因此收敛后的潮流结果仍然是准确的。由于方程的维数减小了，

且B' 和B '' 是常数矩阵，只需在迭代计算之前形成一次，然后分解成因子表，并

一直在迭代过程中使用，所以计算效率大幅提高。快速分解法是一种定雅可比法，虽然只具有线性收敛速度，但由于其鲁棒性好，适应性强，在电力工业界被广泛采用，特别适合在线计算。

2.2快速分解法的理论基础

Stott 的快速分解法提出时并没有任何理论解释，它是计算实践的产物。多年来，人们普遍认为在满足r <

在许多实际应用中，当r >x 时，快速分解法也能很好收敛。因此，从理论上解

释快速分解法的收敛机理，便成为一个有趣的研究课题。20 世纪 80 年代末，Monticelli 等人的研究工作对这一问题做了比较完整的解释，在一定程度上阐明

了 XB 型和BX 型快速分解潮流算法的收敛机理。

Monticelli 等人的分析工作是以定雅可比牛顿-拉夫逊迭代方程为出发点的。具体过程如下：①通过高斯消去法，把牛顿-拉夫逊法的每一次迭代等价地

细分为三步计算；②对每一步计算作详细分析，证明了在连续的两次牛顿-拉夫

逊迭代中，上一次迭代的第三步和下一次迭代的第一步可以合并，从而导出等效的

两步式分解算法；③论证了该两步式分解算法的系数矩阵与快速分解法的系数