当前位置:文档之家› 牛顿-拉夫逊法-复杂电力系统-潮流计算-毕业论文

牛顿-拉夫逊法-复杂电力系统-潮流计算-毕业论文

基于极坐标的牛顿-拉夫逊法潮流计算

摘要

潮流计算是电力系统最基本的计算功能,其基本思想是根据电力网络上某些节点的已知量求解未知量,潮流计算在电力系统中有着独特的作用。它不仅能确保电力网络能够正常的运行工作、提供较高质量的电能,还能在以后的电力系统扩建中各种计算提供必要的依据。

计算潮流分布的方法很多,本设计主要用的是基于极坐标的牛顿-拉夫逊法。根据电力系统网络的基本知识,构建出能代表电力系统系统网络的数学模型,然后用牛顿—拉夫逊法反复计算出各个接点的待求量,直到各个节点的待求量满足电力系统的要求。我们可以画出计算框图,用MATLAB编写出程序,来代替传统的手算算法。复杂电力系统是一个包括大量母线、支路的庞大系统。对这样的系统进行潮流分析时,采用人工计算的方法已经不再适用。计算机计算已逐渐成为分析复杂系统潮流分布的主要方法。本设计中还用了一个五节点的电力系统网络来验证本设计在实际运行中的优越性。

关键词:牛顿-拉夫逊法,复杂电力系统,潮流计算

The method of Newton- Raphson based on polar

ABSTRACT

Power system load flow calculation is the most basic computing functions, the basic idea is based on some of the electricity network nodes to solve the unknown quantity of known volume,In power system, power flow, which can ensure that electrical net can work well and give the high quality power, but also later provide the necessary datas in the enlargement of the power system. has special function.

There are lots of methods about power flow. We mainly use the method of

Newton-Raphson based on polar in my design. According to the basic knowledge of the electrical network, we established the mathematics model which can presents the power system ,then computed again and again unknown members of the each bus with the method of Newton-RaphSon until the unknown numbers meet the demand of the power system. We can write down the block diagram and write the order with the Matlab in place of the traditional methods. Complex power system is a large system which involves lots of bus bars and branches. We also chose a five-bus power system for testing the advantages in the relity.

KEY WORDS: Newton-Raphson,power system,power flow

目录

前言 (1)

第一章电力系统潮流计算的基本知识 (2)

1.1潮流计算的定义及目的 (2)

1.2潮流计算方法的发展及前景................ 错误!未定义书签。第二章潮流计算的节点. (5)

2.1 节点的分类 (5)

2.2潮流问题变量的约束条件 (7)

第三章电力网络的数学模型 (8)

3.1 节点导纳矩阵的形成 (9)

3.2 节点导纳矩阵的修改 (9)

第四章潮流计算的原理 (12)

4.1 牛顿-拉夫逊法 (12)

第五章计算实例 (17)

5.1算例 (17)

5.2 节点导纳的形成 (17)

5.3 计算结果 (18)

结论 (20)

谢辞 (22)

参考文献 (23)

附录 (24)

计算程序 (25)

外文资料翻译 (41)

前言

潮流计算是电力系统中应用最广泛和最重要的一种电气计算。其任务是根据给定的网络结构及运行条件,求出整个网络的运行状态,其中个母线的电压、网络中的功率分布以及整个系统的功率损耗等。潮流计算可以分为简单网络的潮流计算和复杂系统的潮流计算。简单网络的潮流计算,比如:辐射型网络的潮流计算和闭式网络的潮流计算。它们是复杂电力系统潮流计算的基础。在复杂的电力系统潮流计算中需要对电力系统网络进行必要的计算,用来获得必要的数据。潮流计算在电力系统规划设计及运行方式分析的离线及在线计算中都发挥着重要的作用。在这个设计中,我们选折了MATLAB开发潮流计算程序,是因为潮流计算在数学上一般属于多元非线性代数方程组的求解,必须采用迭代计算其中涉及大量的向量和矩阵运算,使用传统的编程语言将十分麻烦。而MATLAB以复数矩阵为基本运算单元,且内置众多高精度、高可靠性矩阵、数组运算函数、、数值计算方法,可大大提高编程的效率。

第一章电力系统潮流计算的基本知识

1.1潮流计算的定义及目的

电力系统潮流计算分布计算,是指电力系统在某一稳定状态的正常运行方式下,电力网络各节点的电压和功率分布的计算。它的主要目的:

(1)检查电力系统各元件是否过负荷。

(2)检查电力系统各节点的电压是否满足电压质量的要求。

(3)根据对各种运行方式的潮流分布计算,可以帮助我们正确地选择系统的

接线方式,合理调整负荷,以保证电力系统安全、可靠地的运行,向用户供给高质量的电能。

(4)根据功率分布,可以选折电力系统的电气设备和导线截面积,可以为电

力系统继电保护整定计算提供必要的数据等。

(5)为电力系统的规划和扩建提供依据。

(6)为调整计算、经济运行计算、短路计算和稳定计算提供必要的数据。

在计算机技术还未发展以前,电力系统的潮流分布计算多采用“手工”近似计算,即按照电路的基本关系,用手工来推算各节点的功率和电压。随着电子计算机技术的进步,电力系统潮流分布的计算几乎已普遍采用计算机来进行,通过求解描述电力系统状态的数学模型,而得到较精确的解。

潮流计算是电力系统分析中的一种最基本的计算,它的任务是对给定的运行条件确定系统的运行状态,如各母线上的电压、网络中的功率分布以及功率损耗等。

1.2潮流计算方法的发展及前景

在用数字计算机求解电力系统潮流问题的开始阶段,人们普遍采用以节点导纳矩阵为基础的高斯-赛德尔迭代法。这个方法的原理比较简单,要求的数字计算机的内存量也比较小,适应当时的电子数字计算机制作水平和电力系统理论水平,于是电力系统计算人员转向以阻抗矩阵为主的逐次代入法。

20世纪60年代初,数字计算机已经发展到第二代,计算机的内存和计算速度发生了很大的飞跃,从而为高斯-赛德尔迭代法的采用创造了条件。阻抗矩阵是满矩阵,高斯-赛德尔迭代法要求计算机储存表征系统接线和参数的阻抗矩阵。这就需要较大的内存量。而且阻抗法每迭代一次都要求顺次取阻抗矩阵中的每一个元素进行计算,因此,每次迭代的计算量很大。

高斯-赛德尔迭代法改善了电力系统潮流计算问题的收敛性,解决了导纳法无法解决的一些系统的潮流计算,在当时获得了广泛的应用,曾为我国电力系统设计、运行和研究做出了很大的贡献。但是,高斯-赛德尔迭代法的主要缺点就是占用计算机的内存很大,每次迭代的计算量很大。当系统不断扩大时,这些缺点就更加突出。为了克服阻抗法在内存和速度方面的缺点,后来发展了以阻抗矩阵为基础的分块阻抗法。这个方法把一个大系统分割为几个小的地区系统,在计算机内只需存储各个地区系统的阻抗矩阵及它们之间的联络线的阻抗,这样不仅大幅度的节省了内存容量,同时也提高了节省速度。

克服高斯-赛德尔迭代法缺点的另一途径是采用牛顿-拉夫逊法。牛顿-拉夫逊法是数学中求解非线性方程式的典型方法,有较好的收敛性。解决电力系统潮流计算问题是以导纳矩阵为基础的,因此,只要在迭代过程中尽可能保持方程式系数矩阵的稀疏性,就可以大大提高牛顿潮流程序的计算效率。自从20世纪60年代中期采用了最佳顺序消去法以后,牛顿-拉夫逊法在收敛性、内存要求、计算速度方面都超过了阻抗法,成为直到目前仍被广泛采用的方法。

在牛顿-拉夫逊法的基础上,根据电力系统的特点,抓住主要矛盾,对纯数学的牛顿法进行了改造,得到了P-Q分解法。P-Q分解法在计算速度方面有显著的提高,迅速得到了推广。

牛顿-拉夫逊法的特点是将非线性方程线性化。20世纪70年代后期,有人提出采用更精确的模型,即将泰勒级数的高阶项也包括进来,希望以此提高算法的性能,这便产生了保留非线性的潮流算法。另外,为了解决病态潮流计算,出现了将潮流计算表示为一个无约束非线性规划问题的模型,即非线性规划潮流算法。

近20多年来,潮流算法的研究仍然非常活跃,但是大多数研究都是围绕改进牛顿法和P-Q分解法进行的。此外,随着人工智能理论的发展,遗传算法、人工神经网络、模糊算法也逐渐被引入潮流计算。但是,到目前为止这些新的模型和算法还不能取代牛顿-拉夫逊法和P-Q分解法的地位。由于电力系统规模的不断扩大,对计算速度的要求不断提高,计算机的并行计算技术也将在潮流计算中得到广泛的应用,成为重要的研究领域

第2章 潮流计算的节点

2.1 节点的分类

节点电压方程是潮流计算的基础方程式。在电气网络理论中,一般是给出电压源或电流源,为求得网络内电流和电压的分布,只要直接求解网络方程就可以了。但是,在潮流计算中,在网络的运行状态求出以前,无论是电源的电势值,还是节点的注入的电流,都是无法准确给定的。

图2-2表示某个三节点的简单电力系统及其等值电路,其网络方程为

111213212223313233311232123312Y Y Y U U U I Y Y Y U U U I Y Y Y U U U I ???

?

??????

??

?=++???

=++???

=++??

1

2

3

312i i i i Y Y Y U U U I ?

???

=++(i =1,2,3) (2-1)

因为s ui ?*

=,所以节点电流用功率和电压可以表示为

))((Gi

LDi Gi

LDi

j

i i

i i

Q

Q

P P S I U U *?

*

*

---==

(2-2)

将式(2-2)带入式(2-1)可得

))((1

2

3

312Gi

LDi Gi LDi

i i i j

i

Q

Q

P

P Y Y Y U U U U ?

??

*

---=

++

这是一组复数方程式,如果把实部和虚部分开,便得到6个实数方程。但是每个节点都有6个变量,即发电机发出的有功功率和无功功率、负荷需要的有功功率和无功功率,以及节点电压的幅值和相位(或对应与某一个参考直角坐标的实部和虚部)。对于n 个节点的网络,可以写2n 个方程,但是确有6n

个变量。因此,对于每个节点,必须给定这6个变量中的4个,使待求量的数目同方程的数目相等,才能对方程求解。

通常把负荷功率作已知量,并把节点功率i Gi LDi P P P =-和i

Gi

LDi

Q Q Q

=-引入网络方程。这样n 个节点的电力系统潮流方程的一般形式可以写为

1

n i

i

ij

j

j j Q P Y U U

*

=-?

=∑(i =1,2,….n )

1

n

i i

ij

j

j j i Q P U U U

*

=**

+=∑

(2-3)

将上述方程的实部和虚部分开,对每一个节点可得2个实数方程,但是变量仍还有4个,即P 、Q 、U 、δ。还要给定其中的2个,将剩下的2个作为待求变量,方程组才可以求解。根据电力系统的实际运行条件,按给定变量的不同,一般将节点分为以下三种类型。 1.PQ 节点

这类节点的有功功率P 和无功功率Q 是给定的。节点电压(U ,δ)是待求量。通常变电所都是在这一类型的节点,由于没有发电机设备,故发电机功率为零。若系统中某些发电厂送出的功率在一定时间内为固定时,则该发电厂母线可作为PQ 节点。可见电力系统的绝大多数节点属于这一类型。 2.PU 节点

这类节点的有功功率P 和电压幅值U 是给定的,节点的无功功率Q 和电压的相位δ是待求量。这类节点必须有足够的可调无功容量,用以维持给定的电压幅值,因而又称之为电压控制节点。一般是选折有一定无功储备的发电厂和具有可调无功电源设备的变电所作为PU 节点。在电力系统只能中,这一类的数目很少。 3.平衡节点

在潮流分布算出以前,网络中的功率损失是未知的,因此,网络中至少有一个节点的有功功律P 是不能给定的,这个节点承担了系统有功功率的平衡,故称之为平衡节点。另外,必须选定一个节点,指定其电压相位

为零,作为计算各节点电压相位的参考,这个节点称为基准节点。基准节点的电压幅值是给定的。(亦称为松弛节点、摇摆节点)。电力系统中平衡节点一般只有一个,它的电压幅值和相位是给定的,而其有功功率和无功功率是待求量。

一般选折主调频发电厂为平衡节点比较合适。但在进行潮流计算时也可以按照惯例的原则来选折,例如,为了提高导纳矩阵法潮流程序收敛性,也可以选折出线较多的发电厂母线做为平衡节点。

根据以上所述可以看到,尽管网络方程是线性方程但是由于在定解条件中不能给定节点电流,只能给出节点功率,这就使潮流方程变为非线性方程了。由于平衡节点的电压已给定,只需要计算其余(n -1)个节点的电压。所以方程式的数目实际上只有2(n -1)个。

2.2潮流问题变量的约束条件

通过求解方程得到了全部节点电压以后,就可以进一步计算各类节点的功率以及网络中功率的分布。这些计算结果代表了潮流方程在数学上的一组解答。但这组解答所反映的系统运行状态,在工程上是否具有实际意义还需要进行检验,因为电力系统运行必须满足一定技术上和经济上的要求。这些要求构成了潮流问题中某些变量的约束条件,通常的约束条件有:

1. 所有节点电压必须满足

min

max i i i U

U U ≤≤(i =1,2,……n)

这个条件是说各节点电压的幅值应限制在一定的范围之内。从保证电能质量和供电安全的要求来看,电力系统的所有电气设备都必须运行在额定电压附近。对于平衡节点的PU 节点,其电压幅值必须按上述条件给定。因此,这一约束条件主要是对PQ 节点而言。

2. 所有电源节点的有功功率和无功功率必须满足的条件

min

max

Gi Gi

Gi p

p

p

和min

max

Gi Gi

Gi Q

Q Q

≤≤

PQ 节点的有功功率和无功功率以及PU 节点的有功功率,在给定时就必须满足上式条件。因此对平衡节点的P 和Q 以及PU 节点Q 应按上述条件进行检验。

3. 某些节点之间电压的相位差应满足

max

i

j

i j δδ

δδ

-<-

为了保证系统运行的稳定性,要求某些输电线路两端电压相位差不超过一定的数值。这一约束的主要意义就在于此。如果计算出来的结果不满足这些约束条件,必须修改某些变量的给定值,甚至修改系统的运行方式。

第3章 电力网络的数学模型

电力网络的数学模型指的是将网络有关参数和变量及其相互关系归纳起来所组成的、可以反映网络性能的数学方程式组。也可以说是对电力系统的运行状态、变量和网络参数在电力系统潮流分布的计算中,广泛采用的是节点电压方程。

在电工原理课中,已讲过用节点导纳矩阵表示的节点电压方程为

n n n Y U I ??????

??=??????????

(3-1) 对于n 个节点的网络,它可以展开为

11111212122

221

2

3............2..................

.........

n n n n n n n U I Y Y Y Y Y

Y I U Y Y

Y I U ???????????

??

?????????????????????

??=?????????

??????????

?????????????????????

??

(3-2) 式(3-1)中的是节点注入电流的列向量。是节点电压的列向量。网络中有接地支路时,通常以大地为参考点,节点电压就是各节点的对地电压。

3.1 节点导纳矩阵的形成

节点导纳矩阵的对角线元素称为自导纳。由式(3-2)可见,自导纳ii Y 等于

在节点i 施加单位电压i U ?

,其它节点全部接地时,经过点i 向网络中注入的电

流,亦等于与节点相连支路的导纳之和。其表示式为

(1,0,)ii

ij

i j ij

j i

j i

i

i i j i I

y U U Y I

I U

?

?

∈∈======≠∑∑ (3-3)

节点导纳矩阵的非对角线元素ij Y (i =1、2、…n ,j=1、2、…n 但i ≠j)称为互导纳。由式(3-3)可见,互导纳ij Y 在数值上就等于节点i 施加单位电压,其它节点全部接地时,经节点j 注入网络的电流。亦等于节点i ,j 之间所连支路元件导纳的负值,其表示式为

(

1,0,)ij

i

j ji

i

j i j j I y I U

U Y U

?

?

===-==≠ (3-4)

依据互导纳的物理意义可知ij Y =ij

y -,即ij Y =ji Y ;特别地,当节点i 、j 之间无直接支路相连时,ij Y =ji Y =0。在复杂电力网中,这中情况较多,从而使矩阵中出现了大量的零元素、节点导纳矩阵称为稀疏矩阵。一般来说

ii

ij Y

Y >,即对角线元素的绝对值大于非对角线元素的绝对值,使节点导纳矩

阵称为具有对角线优势的矩阵。因此节点导纳矩阵是一个对称、稀疏且具有对角线优势的方阵。

3.2 节点导纳矩阵的修改

在电力系统中,接线方式或运行状态等均会发生变化,从而使网络接线改变。比如一台变压器支路的投入或切除,均会使与之相连的节点的自导纳或互导纳发生变化,而网络中其它部分结构并没有改变,因此不必从新形成节点导

纳矩阵,而只需对原有的矩阵作必要的修改就可以了。现在几种典型的接线变化说明具体的修改方法。

(1)从原有网络的节点i 引出一条导纳为ij Y 的支路(见图3-1,(a )),j 为新增加的节点,由于新增加了一个节点,所以节点导纳矩阵增加一阶,矩阵作如下修改:

1)原有节点i 的自导纳ii Y 的增量ii Y ?=ij

y

2)新增节点j 的自导纳ij Y =

ij

y

3)新增的非对角元素ij Y =ji Y =-

ij

y

;其它新增的非对角元均为零。(2)

在原有网络的节点i 与j 之间增加一条导纳为ij

y

的支路(见图3-1,(b )),则

与i 、j 有关的元素应作如下修改:

1)节点i 、j 的自导纳增量ii Y ?=jj Y ?=

ij

y

2)节点i 、j 的互导纳增量ij Y ?=ji Y ?=-ij

y

(3)在网络的原有节点i 、j 之间切除一条导纳为

ij

y

的支路,(见图3-1,(c ))),

其相当在i 、j 之间增加一条导纳为-ij

y

的支路,因此与i 、j 有关的元素应作

以下修改:

1)节点i 、j 的自导纳增量ii Y ?=jj Y ?=-

ij

y

2)节点i 、j 之间的互导纳增量ij Y ?=ji Y ?=ij

y

(4)原有网络节点i 、j 之间的导纳由

ij y 变成'

ij y 见图3-1,

(d )),相当于

在节点i 、j 之间切除一条导纳为

ij y 的支路,在增加一条导纳为'

ij y 的支路,

则与i 、j 有关的元素应作如下修改:

1)节点i 、j 的自导纳增量ii Y ?=jj Y ?='

ij y -ij y ;

2)节点i 、j 的互导纳增量ij Y ?=ji Y ?=ij y -'

ij y 。

y ij

ij

y′ij

图(3-1)

第4章 潮流计算的原理

4.1 牛顿-拉夫逊法

设有单变量非线性方程

()0f x = (4-1)

求解此方程时。先给出解的近似值

(0)

x 它与真解的误差为(0)

x

?,则

(0)

(0

x x

x =

+?将满足方程,即

(0)

(0)

()0f x

x +?= (4-2)

将(3-8)式左边的函数在(0)

x

附近展成泰勒级数,于是便得

2

'

''

()

(0)

(0)

(0)

(0)(0)

(0)

(0)

(0)

(0)

()()()()

......()

....

2!

!

()

()

n

n f f n x x f

f

f

x x x x x x x +?=+

?+

++

+?? (4-3)式中,

'

(0)

()f x

,……

()

(0)()n f

x 分别为函数()f x 在(0)

x

处的一阶导数,….,n 阶导

数。

如果差值(0)

x

?很小,(3-9)式右端(0)

x

?的二次及以上阶次的各项均

可略去。于是,(3-9)便简化为

'

(0)(0)(0)

(0)(0)

()()()f f f

x x x x x +?=+

?=0 (4-4)

这是对于变量的修正量(0)

x

?的现行方程式,亦称修正方程式。解此方程

可得修正量 (0)

(0)

'

(0)

()

()

f x x f x

?=-

(4-5)

用所求的(0)

x

?去修正近似解,变得

(0)

(1)

(0)(0)

(0)

'

(0)

()

()

f x x

x x

x f x

=+?=-

(4-6)

由于(3-10)是略去高次项的简化式,因此所解出的修正量(0)

x ?也只是

近似值。修正后的近似解

(1)

x

同真解仍然有误差。但是,这样的迭代计算

可以反复进行下去,迭代计算的通式是

()

(1)

()

'

()

()

()

k k k k f x x

x f x

+=-

(4-7)

迭代过程的收敛判据为

()

1()k f x ε< (4-8) 或 ()2k x

ε?< (4-9)

式中

1

ε

2

ε

为预先给定的小正数。

这种解法的几何意义可以从图3-1得到说明。函数y =f(x)为图中的曲线。f(x)=0的解相当于曲线与x 轴的交点。如果第k 次迭代中得到()

k x

则过()

()

()

,()k k k f y x

x ??=????

点作一切线,此切线同x 轴的交点便确定了下一个近似值(1)

k x

+。由此可见,牛顿-拉夫逊法实质上就是切线法,是一种逐步

线性化的方法。

应用牛顿法求解多变量非线性方程组(3-1)时,假定已给出各变量的初值

1

(0)x ,2

(0)x

….

(0)n

x

,令

1

(0)x ?,2(0)x ?,….. (0)n

x

?分别为各

变量的修正量,使其满足方程(3-1)即

11

12221

1221122(0)(0)(0)(0)(0)(0)(,,....,)0

(0)(0)(0)(0)(0)(0)(,,....,)0......(0)(0)(0)(0)(0)(0)(,,....,)0n n n n n

n n f x x x x x x f x x x x x x f x x x x x x ?+?+?+?=??+?+?+?=??

?

?

+?+?+?=?? (4-10) 将上式中的n 个多元函数在初始值附近分别展成泰勒级数,并略去含有

1

(0)x ?,2

(0)x ?,……,(0)n

x ?二次及以上阶次的各项,便得

11100011212121110002121212101211(0)(0)(0)(0)(0)(0)(,,...,)...0(0)(0)(0)(0)(0)(0)(,,...,)...0......(0)(0)(0)(0)(,,...,)|||||||n n n n n n n n f f f f x x x x x x x x x f f f f x x x x x x x x x f f x x x x x ???+?+?++?=??????+?+?++?=?????+?+?110022(0)(0) 0

||n

n f f x x x x ????

???

?

?

?

???++?=????

(4-11) 方程式(3-17)也可以写成矩阵形式

11100

0121122

2221200

1

2

1200

01

2

...(0)(0)(0)(,,...,)(0)(0)(0)(,,...,)..................

...

(0)(0)(0)(,,...,)...|||||||||n

n n n

n n n n n n f f

f x x x f x x x f f f f x x x x x

x f

x x x f f f x x

x ???

?

??????????

?????????=-????????

??

?

?

??

?????1

2

(0)(0)...(0)n x x x ???????????

?

?????

??

????

?????????

?????

(4-12) 方程式(3-18)是对于修正量1

(0)x ?,2

(0)x ?,……, (0)n

x ? 的线性方程组,称

为牛顿法的修正方程式.利用高斯消去法或三角分解法可以解出修正量

1

(0)x ?,2

(0)x ?,……, (0)n

x ?。然后对初始近似值进行修正

(1)(0)(0)i

i

i

x

x x

=+? (i=1,2,….,n) (4-13)

如此反复迭代,在进行k +1次迭代时,从求解修正方程式

111121122

222121

2

121

2

...()()()(,,...,)()()()(,,...,)..................

...

()()()(,,...,)...|||||||||k k

k n

n n k k

k

n

n n n n n k k

k n k k k k k k k k k f f

f x x x f x x x f f f f x x x x x

x f

x x x f f f x x

x ???

?

??????????

????????

?=-????????

??

?

?

??

?????1

2

()()...()n k k k x x x ???????????

?

?????

??

????

?????????

?????

(4-14) 得到修正量1

()k x ?,2

()k x ?,()n

k x ?,并对各变量进行修正

(1)()()i

i

i

k k k x

x x

+=+? (i=1,2,…,n) (4-15)

式(3-20)和(3-21)也可以缩写为 ()()

()

()k k k F J X X

=-? (4-16) 和

(1)

()

()

k k k X

X

X

+=

+? (4-17)

式中的X 和X ?分别是由n 个变量和修正量组成的n 维列向量;F(X)是由n 个多元函数组成的n 维列项量;J 是n 阶方阵,称为雅可比矩阵,它的第i 、j 个元素i ij i

f J x

=

??是第n 个函数

1

2

(,,...,,)n

i

f x x x 对第j 个变量j

x

的偏导数;

上角标(k)表示J 阵的每一个元素都在点,

,

,

()()()(...,)12i

k k k n f x x x 处取值。

迭代过程一直到满足收敛判据 {}1

1

2

()()()max

(,,...,)i

n

k k k f x x x ε< (4-18)

或 {

}

2

()max i

k x

ε

?<

(4-19)

为止。

1

ε

2

ε

为预先给定的小正数。

将牛顿-拉夫逊法用于潮流计算,要求将潮流方程写成形如方程式(3-1)的形式。由于节点电压可以采用不同的坐标系表示,牛顿-拉夫逊法潮流计算也将相应的采用不同的计算公式。

图(4-1)牛顿-拉夫逊方法的几何意义

第5章 计算实例

5.1算例

图1为一五结点系统,各支路参数均为标么值。假定结点1、2、3为PQ 节点,结点4为PV 节点、结点5为平衡结点,试分别用直角坐标和极坐标牛顿-拉夫逊法计算其潮流。取收敛判据为|?P i |≤10-5和|?Q i (?V i 2)|≤10-5。给定:

S 1=-1.6-j0.8 S 2=-2.0-j1.0 S 3=-3.7-j1.3 P 4=5.0 |V 1(0)|=|V 2(0)|=|V 3(0)|=|V 4(0)|=1.00 |V 4|=|V 5|=1.05

00.1e e e (0)3(0)2(0)1=== 05.1e e 5(0)4==

图1 - 5节点系统

5.2 节点导纳的形成

根据节点导纳矩阵的定义,可求的节点导纳矩阵各元素,即

5(0)4(0)3(0)2(0)1=====f f f f f

相关主题
相关文档 最新文档