潮流计算的快速分解法(可编辑修改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 型算法的处理方式是：在对系数矩阵BH影响，但忽略了接地支路项。

BX 型算法的迭代格式与 XB 型算法是相同的。

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

快速分解法只对雅可比矩阵作了简化，但节点功率偏差量的计算及收敛条件仍是严格的，因此收敛后的潮流结果仍然是准确的。

由于方程的维数减小了，且B' 和B '' 是常数矩阵，只需在迭代计算之前形成一次，然后分解成因子表，并一直在迭代过程中使用，所以计算效率大幅提高。

快速分解法是一种定雅可比法，虽然只具有线性收敛速度，但由于其鲁棒性好，适应性强，在电力工业界被广泛采用，特别适合在线计算。

2.2快速分解法的理论基础Stott 的快速分解法提出时并没有任何理论解释，它是计算实践的产物。

多年来，人们普遍认为在满足r <<x 的系统中，快速分解法才能有较好的收敛性。

但在许多实际应用中，当r >x 时，快速分解法也能很好收敛。

因此，从理论上解释快速分解法的收敛机理，便成为一个有趣的研究课题。

20 世纪 80 年代末，Monticelli 等人的研究工作对这一问题做了比较完整的解释，在一定程度上阐明了 XB 型和BX 型快速分解潮流算法的收敛机理。

Monticelli 等人的分析工作是以定雅可比牛顿-拉夫逊迭代方程为出发点的。

具体过程如下：①通过高斯消去法，把牛顿-拉夫逊法的每一次迭代等价地细分为三步计算；②对每一步计算作详细分析，证明了在连续的两次牛顿-拉夫逊迭代中，上一次迭代的第三步和下一次迭代的第一步可以合并，从而导出等效的两步式分解算法；③论证了该两步式分解算法的系数矩阵与快速分解法的系数LMH ⎩ P P L 矩阵是一致的。

推导过程并未引用任何解耦的假设。

为以后书写方便，将式（2.1）中的∆P V 用∆P 代替， ∆Q V 用∆Q 代替， 而V ∆用∆代替，则给出的定雅可比法的修正公式改写如下- ⎡ H N ⎤⎡∆⎤ = ⎡∆P ⎤ ⎢M L ⎥⎢∆V ⎥ ⎢∆Q ⎥（2.6）式中⎣ ⎦⎣ ⎦ ⎣ ⎦ H = B H ≈∂∆P， N = -G ∂T N≈ ∂∆P∂V T M = G M≈ ∂∆Q ， L = B ∂T L≈ ∂∆Q ∂V T 整个推导分为三步。

1) 将原问题分解成 P ， Q 子问题首先，对式（2.6）用高斯消去法消去子块 N ，有- ⎡H - NL -1M 0⎤⎡∆⎤ = ⎡∆P - NL -1∆Q ⎤ ⎢ M L ⎥⎢∆V ⎥ ⎢ ∆Q ⎥（2.7）记⎣ ⎦⎣ ⎦ ⎣ ⎦~ -1 ~ -1 H = H - NL M ， ∆P = ∆P - NL ∆Q并定义∆V = -L -1∆Q ， ∆V= -L -1M ∆则式（2.7）的解可以表示为⎧∆= - ~-1∆P⎨∆V = ∆V + ∆VM上式中对∆~的计算可以采用较简单的方法。

在给定的电压幅值和相角初值附近，保持电压相角不变，考虑只有电压幅值的变化∆V L 时，有功功率的偏差量 为∆P (,V + ∆V L) ≈ ∆P (,V ) + ∂∆P∆V ∂V TL= ∆P (,V ) - NL -1∆Q = ∆~⎩H ∆P ( ,V （2.8）综合上述结果，如果当前的迭代点为((k ) ,V (k ) ) ，则第k 次迭代对式（2.6）的计算可以分解为以下三步。

①⎧∆V (k ) = -L -1∆Q ((k ) ,V (k ) ) ⎨ ~ L V (k +1) = V (k ) + ∆V (k )（2.9）②⎩ L⎧⎪∆(k ) = - ~ -1 (k )~(k +1)（2.10）③⎨⎪(k +1) = (k ) + ∆(k )⎧∆V (k ) = -L -1 M ∆(k )⎨ M ~ V (k +1) = V (k +1) + ∆V (k ) ⎩ M（2.11）2) 简化无功迭代步骤按①～③完成第k 次迭代后，下面再考察第k + 1次迭代的①，有⎧∆V (k +1) = -L -1∆Q ((k +1) ,V (k +1) ) ⎨ ~ LV (k +2) = V (k +1) + ∆V (k +1) ⎩ L（2.12）利用式（2.11），上式中的无功功率偏差为∆Q ((k +1) ,V (k +1) ) = ∆Q ((k +1)~(k +1)+ ∆V (k ) ) ,V M≈ ∆Q ((k +1)~(k +1) ) + ∂∆Q ∆V (k ),V ∂V TM = ∆Q ((k +1)~(k +1)) + L ∆V (k ) ,V M（2.13）代入式（2.12），经整理得∆V (k +1) + ∆V (k ) = -L -1∆Q ((k +1)~(k +1))L M,V（2.14）)V H ⎩H HH H 式（2.14）说明，如果将第k 次迭代的①计算出的 ~(k +1) 和②计算出的(k +1) ，用于计算第k + 1次迭代的无功偏差量，即式（2.14）中的∆Q ，则所求得的第k + 1次迭代的电压修正量将自动包含第k 次迭代的③的式（2.11）与第k + 1次迭代的①的式（2.12）合并，只需保留式（2.9）和式（2.10）。

因此，第k 次迭代对式（2.6）的计算可以用以下两步计算完成：⎧⎪∆V (k ) = -L -1∆Q ((k ) ,V (k ) )（2.15）⎨ ⎪⎩V (k +1)= V (k )+ ∆V (k ) ⎧⎪∆(k )= - ~ -1 ∆P ((k ),V (k +1) )（2.16）⎨⎪(k +1) = (k ) + ∆(k )在式（2.6）处已说明， ∆P 实际是∆P V ， ∆Q 实际是∆Q V ， ∆实际是 V ∆，式（2.15）和式（2.16）和快速分解法迭代格式相同。

显然，这种迭代算法是否与快速分解法等效，取决于系数矩阵 L 和 ~。

与 XB 型快速分解法的修正 方程相比，系数矩阵 L 是导纳矩阵的虚部，这与 B'' 相同，所以关键要看~是否与 B ' 有相同或相似的关系。

3) 简化有功迭代矩阵 ~~H 的定义，有~ -1 -1H = H - NL M = B H +G N B L G M（2.17）对于一般的电网， ~ 可能有较复杂的结构。