复杂网络N-R法潮流分析计算

编号

课程设计

（ 2012级本科）

系（部）院：物理与机电工程学院

专业：电气工程及其自动化

作者姓名：

指导教师：职称：副教授

完成日期： 2015 年 6 月 30 日

二○一五年六月

目录

1 牛顿－拉夫逊法概述 (1)

1.1牛顿-拉夫逊法基本原理 (1)

1.2牛顿--拉夫逊法潮流求解过程 (2)

2手算潮流计算 (6)

2.1用图中数据和等值网络形成节点导纳矩阵Y B (6)

2.2设各节点电压初始值为： (7)

2.3用公式 (7)

2.4求取雅可比矩阵 (7)

2.5求△修正量矩阵 (8)

2.6计算修正各节点电压 (8)

3计算机算法潮流计算 (9)

3.1牛顿—拉夫逊法的程序框图 (9)

3.2结果显示 (9)

总结 (19)

附件 (20)

参考文献 (25)

1 牛顿－拉夫逊法概述

电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷。各点电压是否满足要求，功率的分布和分配是否合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。潮流计算结果可用如电力系统稳态研究，安全估计或最优潮流等对潮流计算的模型和方法有直接影响。实际电力系统的潮流技术那主要采用牛顿-拉夫逊法。

1.1牛顿-拉夫逊法基本原理

牛顿--拉夫逊法(简称N —R 法)在数学上是求解非线性代数方程式的有效方法。其要点是把非线性方程式的求解过程变成反复地对相应的线性方程式进行求解的过程。即通常所称的逐次线性化过程。

对于非线性代数方程组：f (x )=0即

f i (x 1,x 2,…,x n )=0 (i =1,2,…,n ) (1-1-1)

在待求量x 的某一个初始估计值(0)

x 附近，将上式展开成泰勒级数并略去二阶及以上的高阶项，得

到如下的经线性化的方程组：

f(x (0))+f ′(x (0))∆x (0)=0 (1-1-2)

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

∆x (0)=−[f ′(x (0))]−1

f(x (0)) (1-1-3)

将(0)

x

和(0)

x

相加，得到变量的第一次改进值(1)x 。接着就从(1)

x 出发，重复上述计算过程。因此

从一定的初值(0)

x

出发，应用牛顿法求解的迭代格式为：

f ′(x (k ))∆x (k )=−f(x (k )) (1-1-4)

x (k+1)=x (k )+∆x (k ) (1-1-5)

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

有上式可见，牛顿法的核心便是反复形式并求解修正方程式。牛顿法当初始估计值(0)

x 和方程

的精确解足够接近时，收敛速度非常快，具有平方收敛特性。

牛顿潮流算法突出的优点是收敛速度快，若选择到一个较好的初值，算法将具有平方收敛特性，一般迭代4~5次便可以收敛到一个非常精确的解。而且其迭代次数与所计算网络的规模基本无关。牛顿法也具有良好的收敛可靠性，对于对以节点导纳矩阵为基础的高斯法呈病态的系统，牛顿法也能可靠收敛。牛顿法所需的内存量及每次迭代所需时间均较高斯法多。

牛顿法的可靠收敛取决于有一个良好的启动初值。如果初值选择不当，算法有可能根本不收

敛或收敛到一个无法运行的节点上。对于正常运行的系统，各节点电压一般均在额定值附近，偏移不会太大，并且各节点间的相位角差也不大，所以对各节点可以采用统一的电压初值(也称为平直电压)，如假定：

U i

(0)

=1 θi (0)

=0 或 e i (0)

=1 f i

(0)

=0 (i =q1,2,…,n,i ≠s ) (1-1-6)

这样一般能得到满意的结果。但若系统因无功紧张或其它原因导致电压质量很差或有重载线路而节点间角差很大时，仍用上述初始电压就有可能出现问题。解决这个问题的办法可以用高斯法迭代1~2次，以此迭代结果作为牛顿法的初值。也可以先用直流法潮流求解一次以求得一个较好的角度初值，然后转入牛顿法迭代。

1.2牛顿--拉夫逊法潮流求解过程

以下讨论的是用直角坐标形式的牛顿—拉夫逊法潮流的求解过程。当采用直角坐标时，潮流问题的待求量为各节点电压的实部和虚部两个分量e 1,f 1,e 2,f 2,…,e n ,f n 由于平衡节点的电压向量是给定的，因此待求两共需要2(n -1)个方程式。事实上，除了平衡节点的功率方程式在迭代过程中没有约束作用以外，其余每个节点都可以列出两个方程式。对PQ 节点来说P is 和Q is 是给定的，因而可以写出

{ ∆P i =P is −e i ∑(G ij e j −B ij f j )−f j ∑(G ij f j +B ij e j )=0

j∈i j∈i ∆Q i =Q is −f i ∑(G ij e j −B ij f j )+e j ∑(G ij f j +B ij e j )=0j∈i

j∈i

（1-2-1）

对PV 节点来说，给定量是P is 和Q is ，因此可以列出

{∆P i =P is −e i ∑(G ij e j −B ij f j )−f i ∑(G ij f j +B ij e j )=0j∈i

j∈i ∆V i 2=V is 2

−(e i 2+f i 2)=0

（1-2-2）

求解过程大致可以分为以下步骤： （1）形成节点导纳矩阵 （2）将各节点电压设初值U ，

（3）将节点初值代入相关求式，求出修正方程式的常数项向量 （4）将节点电压初值代入求式，求出雅可比矩阵元素 （5）求解修正方程，求修正向量 （6）求取节点电压的新值

（7）检查是否收敛，如不收敛，则以各节点电压的新值作为初值自第3步重新开始进行狭义次迭代，否则转入下一步

（8）计算支路功率分布，PV 节点无功功率和平衡节点注入功率。