高斯_快速解耦潮流算法

第33卷第3期电网技术V ol. 33 No. 3 2009年2月Power System Technology Feb. 2009 文章编号：1000-3673（2009）03-0053-04 中图分类号：TM712 文献标志码：A 学科代码：470·4051

高斯–快速解耦潮流算法

彭谦1，胡国新2，张利3

（1．华北电力大学电气与电子工程学院，北京市昌平区 102206；

2．北京首钢设计院，北京市石景山区 100043；

3．北京信息科技大学自动化学院，北京市海淀区 100192）

A Load Flow Algorithm Based on Gauss Algorithm and Fast Decoupling Algorithm

PENG Qian1，HU Guo-xin2，ZHANG Li3

（1．School of Electrical and Electronic Engineering，North China Electric Power University，Changping District，

Beijing 102206，China；2．Beijing Shougang Design Institute，Shijingshan District，Beijing 100043，China；3．School of Automation，Beijing Information Science & Technology University，Haidian District，Beijing 100192，China）

ABSTRACT: In load flow calculation by traditional Gaussian algorithm, the convergence of the computation frequently becomes slower when the PV nodes are transformed into PQ nodes. For this reason, an approach using fast decoupling method to process PV nodes is proposed. In the proposed approach, traditional Gaussian algorithm is used to process PQ nodes and the PQ part of the network is eliminated by Gaussian elimination method, thus a network consisting of PV nodes and a balance bus is obtained, and then the obtained network is solved by fast decoupling method to implement fast convergence of the algorithm. The features of the proposed approach are as following: constant Jacobian matrix, less memory occupation, reliable convergence and high computing speed.

KEY WORDS：load flow；Gaussian algorithm；fast decoupling load flow；power system

摘要：高斯法潮流计算中，PV节点转化为PQ节点易造成计算收敛缓慢，对此文中提出了应用快速解耦法处理网络中PV节点的方法。该方法应用传统高斯法处理PQ节点，利用高斯消元法消去网络中的PQ部分，得到了一个由PV节点和平衡节点组成的网络，然后用快速解耦法求解这个网络，从而实现了算法的快速收敛。该方法有恒定的雅可比矩阵、且内存占用量小、收敛可靠、计算速度快。

关键词：潮流计算；高斯法；快速解耦法；电力系统

0 引言

电力系统中潮流计算最常用的方法有：高斯法、牛顿法和快速解耦法[1]。

高斯法根据迭代方程的不同，可分为以导纳矩阵为迭代矩阵的Y-bus法[2]和以阻抗矩阵为迭代矩阵的Z-bus法[3]。Y-bus法内存占用量少，但是收敛速度慢；Z-bus法内存占用量多，但收敛迅速。目前工程应用中广泛采用因子表技术进行潮流计算[4-5]。它使高斯法既有Y-bus法内存占用少的优点，又有Z-bus法收敛迅速的优点。但是对于含有PV 节点的网络，高斯法通常是把PV节点转化为PQ 节点进行处理，这使得迭代收敛缓慢甚至不收敛。因此，目前高斯法通常应用于PV节点较少的配电网络。

牛顿法是求解非线性代数方程的有效方法[6]。它把非线性方程的求解过程变成反复的相应线性方程求解[7]。牛顿法的突出优点是收敛速度快[8]，但是它必须反复形成修正方程[9]，因此对大规模电网进行潮流计算时，计算速度缓慢[10-11]。

由于交流高压电网电抗远远大于电阻，因此电力系统具有如下的物理特性，即有功功率的变化主要决定于电压相位角的变化，而无功功率的变化主要决定于电压模值的变化，因此可以得到解耦的方程组[12]。快速解耦法在进行输电网潮流计算时具有良好的收敛可靠性[13]，但是对于配电网络，由于系统网络不能满足电抗远远大于电阻的条件，会出现迭代次数大大增加甚至迭代不收敛的情况[14]，因此快速解耦法很少被用于配电网潮流计算。

前推回代法及其变形是配电网潮流计算的有效方法[15]，具有编程简单[16]、数值稳定性好[17]、效率高等优点。在辐射状配电网潮流计算中，如果配

54 彭谦等：高斯–快速解耦潮流算法 V ol. 33 No. 3

电馈电线负荷为恒功率负荷，前推回代法与牛顿法有相似的快速收敛性。但是前推回代法不是一种通用算法，对于含有多环网的电网，前推回代法不能保证迭代计算可靠收敛。

为得到一种定雅可比迭代矩阵、适应各种网络类型、快速收敛的潮流计算方法，本文提出将高斯法与快速解耦法结合起来进行潮流计算。

1 高斯–快速解耦算法

1.1 PQ 节点的处理

对网络节点列写节点电压方程 = I YU (1) 设网络总节点数目为n ，展开式(1)得到 =1

=1,2,,n

i ij j

j I Y U i n =∑ " (2) 网络分析时，给定的运行变量一般不是注入电

流值，而是各节点的注入功率，对于节点i ，两者之间有如下关系

?j ??i i i

i i i

S P Q I U U ?== (3)

式中?S 、?U 分别为节点功率、节点电压相量的共轭。把式(3)带入式(2)得到

=1j =1,2,,?n

i i

ij j j i

P Q Y U i n U ?=∑ " (4)

潮流计算时，应用式(4)处理网络中的PQ 节点。

1.2 PV 节点的处理

网络中PV 节点给定的运行参数为节点有功功率和电压幅值，无功功率未知，无法得到节点功率

相量?S

，因此不能应用式(4)进行计算。 对于n 节点网络，设其中平衡节点的数目为1个，PQ 节点的数目为m 个，PV 节点的数目为n ?m ?1个，网络结构如图1所示。

图1 n 节点网络

Fig. 1 Sketch of n node network

图中节点1~m 为PQ 节点，节点m +1~n ?1为

PV 节点，节点n 为平衡节点。应用式(1)对图1网络列写方程得到

PQ PQ

11121321

2223PV PV 313233n n ????????????=????????????????????

Y I Y Y Y Y Y Y Y I Y Y Y Y I (5) 式中T PQ

1[]m U U = "U ，T PV

+1-1[]m n U U = "U 。 对式(5)进行初等行变换，应用第一行第一列元

素Y 11消去其余行第一列的元素，得到 PQ PQ 11121322

23PV PV 323300n n

????????????′′′=????????????′′′?????????? U I Y Y Y Y Y U I Y Y U I (6) 1111,,=2,3ij ij i 1j i j ?′=?Y Y Y Y Y 1PV PV 2111PQ ?′=? I I Y Y I 13111PQ n n ?′=? I I

Y Y I

PV 节点和平衡节点的注入电流变化为初始注入电流与PQ 节点的等效注入电流的代数和。等效系统网络如图2所示。

图2 等效网络 Fig. 2 Equivalent network

图2相对于图1消去了PQ 节点对应的部分，同时把PQ 部分电流等效注入到了PV 节点以及平衡节点中。图1所示网络转化为一个由PV 节点和平衡节点组成的等效网络。

系统中的PV 节点电压幅值已知，相角未知。系统运行时，一般相角的偏差量较小，所以整个系

统的收敛基本由PQ 节点决定，

PV 部分相角可以应用快速解耦法得到，PV 部分相角计算误差随着PQ 部分电压相量计算误差的减小而减小。

快速解耦法的有功功率迭代方程为

??′=P B θ (7)

式中1

=ij ij B X ′?，=ii

ij j i

B B ∈′′?∑，,=1,2,,i j n "。 对于图2所示网络应用式(7)进行计算，ij

B ′为22′Y 、23′Y 、32′Y 、33′Y 中非对角线元素，=ii ij j i

B B ∈′′∑ 。

第33卷 第3期 电 网 技 术 55

设节点k 为PV 节点，网络PQ 部分k 点等效

注入电流为k

I

′ ，k 点注入电流为?(j )k

k

k

P Q U ?，则 k 点的注入功率为

j ???k k k k k k P Q S U I U ??

?′=+??????

(8) 求解式(8)，则k 点初始有功功率变化为k P ′

?Re()k k k k P U I P ′′=+ (9)

2 高斯–快速解耦结合算法计算流程

高斯–快速解耦结合算法的流程如下：

1）输入原始数据，对网络重新编号，使得PQ 节点的编号范围为1~m ，PV 节点的编号范围为

m +1~n ?1，平衡节点编号为n ，修改支路首末节点编号。

2）计算生成导纳矩阵Y ，然后消元形成前m 行的因子表。

3）

用消元的导纳矩阵生成PV 部分的计算矩阵′B ，然后计算′B 的因子表。

4）应用式(3)计算所有节点的注入电流，然后应用导纳矩阵因子表前代m 行。

5）把前代的电流应用式(9)计算，得到PV 部分节点有功功率，然后应用消元的导纳矩阵和电压计算PV 节点的有功功率，两个计算结果相减作 差得到?P 。

6）′B 的因子表前代n ?m ?1行，把第n ?m 行 对应的计算结果赋值为0，然后回代n ?m ?1行，对

PV 节点相角进行修正，然后计算PV 节点电压相量，把计算结果带入导纳矩阵因子表方程。

7）把导纳矩阵因子表回代m 行，计算结果替换PQ 部分对应元素。

8）判断是否收敛。如果收敛，转步骤9）；如果不收敛，转步骤4）

。 9）输出结果，程序结束。

3 算例分析

3.1 输电网计算结果比较

在处理器为1.7 GHz ，内存为512 MB 的个人机上应用VC 6.0分别编写牛顿法、快速解耦法和高斯–快速解耦结合算法程序，收敛条件为电压幅值偏差小于0.000 1。

对输电网络模型进行牛顿法、快速解耦法和高斯–快速解耦法迭代次数的比较。分别采用IEEE3、IEEE5、IEEE9、IEEE14、IEEE30、IEEE57、IEEE118

和IEEE162系统进行测试，迭代结果如表1所示。

表1 3种算法的迭代次数比较

Tab. 1 Iteration number comparison of

3 kinds of algorithm

迭代次数

测试系统牛顿法 快速解耦法

高斯–快速解耦法

IEEE 3 4 6 2 IEEE 5 9 7 15 IEEE 9 4 4 6 IEEE 14 4 5 5 IEEE 30 4 5 6 IEEE 57 4 5 7 IEEE 118

5

5

6

从表1可知，高斯–快速解耦结合算法对于输电网络具有良好的收敛性。 3.2 配电网计算结果比较

应用IEEE 43系统进行测试，牛顿法迭代8次收敛，计算耗时0.812 ms ；快速解耦法不收敛；高斯–快速解耦结合法迭代33次收敛，计算耗时0.656 ms 。

对于配电网络，快速解耦法不能保证可靠收敛，而高斯–快速解耦结合算法仍能保证收敛。虽然迭代次数相对于牛顿法有所增加，但是由于高 斯-快速解耦结合算法是定因子表计算，因此计算速度仍然比牛顿法快。

3.3 3种方法计算时间比较 比较3种方法的计算时间，计算结果如表2 所示。

表2 3种算法的迭代时间比较

Tab. 2 Iteration time comparison of 3 kinds of algorithm

迭代时间/ ms

测试系统牛顿法快速解耦法 高斯–快速解耦法

IEEE 14 0.609 0.219 0.219 IEEE 30 1.312 0.437 0.454 IEEE 118

5.078

1.016

1.265

由表2可知，高斯–快速解耦结合算法比快速

解耦法略慢，但比牛顿法快。随着网络规模的扩大，高斯–快速解耦结合算法的优势更为明显。

4 结论

本文提出了一种高斯–快速解耦相结合的潮流算法。应用高斯法计算网络中的PQ 节点，并用快速解耦法计算网络中的PV 节点。通过算例分析可知，该方法对各种网络均具有良好的收敛性，算法的通用性优于快速解耦法。虽然迭代计算次数较牛顿法有所增加，但由于该方法的雅可比矩阵为恒定的稀疏矩阵，因此计算速度快于牛顿法。高斯–快

56 彭谦等：高斯–快速解耦潮流算法V ol. 33 No. 3

速解耦相结合的潮流算法具有内存占用量小、收敛性强、计算速度快、通用性强的优点，在工程实际中可以广泛应用。

参考文献

[1] 陈珩．电力系统稳态分析[M]．北京：中国电力出版社，1995．

[2] 诸骏伟．电力系统分析[M]．北京：中国电力出版社，1995．

[3] 张伯明，陈寿孙．高等电力网络分析[M]．北京：清华大学出版社，

1996．

[4] 彭谦，姜彤，曲鹏．修正节点接入导纳潮流算法[J]．电网技术，

2007，31(10)：61-63．

Peng Qian，Jiang Tong，Qu Peng．Power flow algorithm by amending injected admittance[J]．Power System Technology，2007，31(10)：61-63(in Chinese)．

[5] 彭谦，张弘鲲，刘海燕，等．基于集抄系统的配电网潮流计算方

法[J]．电网技术，2007，31(7)：69-72．

Peng

Qian，Zhang Hongkun，Liu Haiyan，et al．Load flow algorithm for distribution network based on centralized meter reading system [J]．Power System Technology，2007，31(7)：69-72(in Chinese)．[6] Tinney W F，Hart C E．Power flow solutions by Newton's method

[J]．IEEE Trans-PAS，1967，86(11)：1449-1460．

[7] 周勇．电力系统三相不对称潮流计算[J]．电网技术，1996，20(1)：

24-29．

Zhou Yong．Three phase load flow calculation for unsymmetrical power system[J]．Power System Technology， 1996，20(1)：24-29(in Chinese)．

[8] 朱辉，郭志忠．电流注入模型的电力系统潮流计算[J]．电网技术，

1998，22(11)：39-41．

Zhu Hui，Guo Zhizhong．Power flow algorithm based on current influx model[J]．Power System Technology，1998， 22(11)： 39-41(in Chinese)．

[9] Lrving M R，Sterling M J H．Efficient Newton-Raphson algorithm for

load-flow calculation in transmission and distribution networks [J]．IEE Proc-Generation，Transmission and Distribution，1987，134(5)：325-327．

[10] 孙德胜．电力系统三相潮流的研究[D]．哈尔滨：哈尔滨工业大学

硕士学位论文，1999．[11] 吴政球，荆勇．考虑负序零序非线性求解的三相潮流计算方法

[J]．中国电机工程学报，2001，22(4)：77-81．

Wu Zhengqiu，Jing Yong．Three phase power flow solution based on node injection currents[J]．Proceedings of the CSEE， 2001， 22(4)：77-81(in Chinese)．

[12] Stott B．Review of load-flow calculation methods [J]．Proceedings of

IEEE，1974，62(7)：916-929．

[13] Stott B，Alsac O．Fast decoupled load-flow[J]．IEEE Trans- PAS，

1974，(93)：859-869

[14] 褚文捷，张小平，陈珩．不对称三相P-Q分解法潮流的并行计算

[J]．中国电机工程学报，1995，15(3)：185-192．

Chu Wenjie，Zhang Xiaoping，Chen Heng．Unsymmetrical three- phase P-Q decoupled load flow study[J]．Proceedings of the CSEE，1995，15(3)：185-192(in Chinese)．

[15] 王守相，王成山．基于区间算法的配电网三相潮流计算模型[J]．中

国电机工程学报，2002，22(2)：52-58．

Wang Shouxiang，Wang Chengshan．Distribution three-phase power flow models based on interval algorithm[J]．Proceedings of the CSEE，2002，22(2)：52-58(in Chinese)．

[16] 王成山，王守相．基于区间算法的配电网三相潮流计算及算例分

析[J]．中国电机工程学报，2002，22(3)：58-62．

Wang Chengshan，Wang Shouxiang．Distribution three-phase power flow based on interval algorithm and test results [J]．Proceedings of the CSEE，2002，22(3)：58-62(in Chinese)．

[17] 车仁飞，李仁俊．一种少环配电网三相潮流计算新方法[J]．中国

电机工程学报，2003，23(1)：74-79．

Che Renfei，Li Renjun．A new three-phase power flow method for weakly meshed distribution systems[J]．Proceedings of the CSEE，2003，23(1)：74-79(in Chinese)．

收稿日期：2008-05-25。

作者简介：

彭谦(1980—)，男，博士研究生，研究方向为

配电网自动化，E-mail：pq.ncepu@https://www.360docs.net/doc/1d14356146.html,；

胡国新(1965—)，男，教授级高级工程师，从

事电力系统电气设计方面的工作。

（编辑王晔）彭谦

我国首条750 kV同塔双回输电线路全线贯通我国首条750kV同塔双回架设的输电线路——兰州东—平凉—乾县线路于1月9日正式全线贯通，这标志着我国已通过自主创新掌握了具有国际先进水平的高电压等级同塔双回关键技术。这项技术可以减少输电线路走廊征地30%，降低工程造价20%，在大幅度提高输电线路自然输送功率、节约土地资源、降低工程造价等方面具有重大意义，将为我国750kV主干电网架的快速发展提供宝贵的经验。

750kV兰州东—平凉—乾县输电线路全长442km，途经甘肃、宁夏、陕西3省(自治区)14个县市，是我国西北地区750kV主电网架的枢纽工程，也是西北水、火、风电打捆实施“西电东送”战略的主通道，对于加强陕西、甘肃电网之间的联络，推进西北能源资源在全国范围内的优化配置具有重大意义。该线路是同塔双回架设技术在我国750kV电压等级输变电项目中的第一次工程实践，从关键技术研究到工程建设，完全由我国电网科研人员自主创新完成。