节点导纳矩阵及潮流计算

节点导纳矩阵及潮流计算

目录

摘要 ..................................................... 2

1任务及题目要求 ......................................... 2 2原理介绍 ............................................... 3

2. 1节点导纳矩阵 ...................................... 3 2. 2牛顿-拉夫逊法 ..................................... 4

2. 2. 1牛顿-拉夫逊法基本原理 ........................ 4 2. 2. 2牛顿一拉夫逊法潮流求解过程介绍 ............... 6

3分析计算 .............................................. 11 4结果分析 .............................................. 15 5总结 .................................................. 16

参考资料 ................................................ 17

1/171

节点导纳矩阵及潮流计算

节点导纳矩阵及潮流计算

摘要

电力网的运行状态可用节点方程或回路方程来描述。节点导纳矩阵是以系 统元件的等值导纳为基础所建立的、描述电力网络各节点电压和注入电流之间 关系的线性方程.潮流计算是电力系统分析中的一种最基本的计算，它的任务是 对给定的运行条件确定系统的运行状态，如各母线上的电压（幅值及相角）、网 络中的功率分布及功率损耗等。本文就节点导纳矩阵和潮流进行分析和计算。

1任务及题目要求

题L1初始条件：

其元件导纳参数为：y12= 0 o 5—j 3 , y：3=0. 8—j 4 , y13 = 0 • 75—j2. 5

任务及要求：1）根据给定的运行条件，确定图2所示电力系统潮流计算时 各节点的类型和待求量;

2）求节点导纳矩阵Y；

2/171

节点导纳矩阵及潮流计算

3) 给出潮流方程或功率方程的表达式；

4) 当用牛顿-拉夫逊法计算潮流时，给出修正方程和迭代收敛条件。

2原理介绍

2. 1节点导纳矩阵

节点导纳矩阵既可根据自导纳和互导纳的定义直接求取，也可根据电路知 识中找出改网络的关联矩阵，在节点电压方程的矩阵形式进行求解。本章节我们 主要讨论的是直接求解导纳矩阵。根据节点电压方程章节我们知道，在利用电 子数字讣算机讣算电力系统运行情况时，多采用IYV形式的节点方程式。其中阶 数等于电力网络的节点数。从而可以得到n个节点时的节点导纳矩阵方程组：

力\"+力\"+…+儿小二1「

Y2M + D2V2 +…二 12 … (2

>

KiV1 + K,2V2+-+LtV/i = L

一 1 )

山此可以得到n个节点导纳矩阵：

1 V ... 2

12

Y =

乙2… • • • 厶…

(2-2 )

它反映了网络的参数及接线情况，因此导纳矩阵可以看成是对电力网络电 气特性的一种数学抽象。由导纳短阵所联系的节点方程式是电力网络广泛应用 的一种数学模型.

通过上面的讨论，可以看岀节点导纳矩阵的有以下特点：

(I) 导纳矩阵的元素很容易根据网络接线图和支路参数直观地求得，形成 节点导纳

3/171

节点导纳矩阵及潮流计算

矩阵的程序比较简单。

(2) 导纳矩阵为对称矩阵。山网络的互易特性易知YH.

(3) 导纳矩阵是稀疏矩阵。它的对角线元素一般不为零，但在非对角线元素 中则存

在不少零元素。在电力系统的接线图中，一般每个节点与平均不超过 3、4个其他节点有直接的支路连接。因此，在导纳矩阵的非对角线元素中每行 仅有3〜4个非零元素，其余的都是零元素，而且网络的规模越大，这种现象越显 著。

节点导纳矩阵的形式可归纳如下：

(1) (2)

数.

导纳矩阵的阶数等于电力网络

导纳矩阵各行非对角元素中非零元素的个数等于对应节点所连得不接 地支路

(3) 导纳矩阵各对角元素，即节点的自导纳等于相应节点之间的支路导纳之 和。 (4) 导纳矩阵非对角元素，即节点之间的互导纳等于相应节点之间的支路导 纳的

负值.

2. 2牛顿一拉夫逊法

2. 2O 1牛顿-拉夫逊法基本原理

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

对于非线性代数方程组：

/(x) = o (2—3)

即 …,心)=0 0 = 1,2,…加

在待求量X的某一个初始估计值附近，将上式展开成泰勒级数并略去二 阶及以上的高阶项，得到如下的经线性化的方程组：

4/171

节点导纳矩阵及潮流计算

/(x(o)) + /(x(o,)zkv(o)=0

(2-4)

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

Ar(o)=_[/(?o))r7(?o)) (2-5)

将心⑹和x⑹相加，得到变量的第一次改进值x⑴.接着就从2\"出发，重复 上述计算过程。因此从一定的初值x⑹出发，应用牛顿法求解的迭代格式为：

f0)3=-

(2—6)

严=卅>+2>

(2—7)

上两式中：f(x)是函数/(X)对于变量x的一阶偏导数矩阵,即雅可比矩阵 J；k为迭代次数.

由上式可见，牛顿法的核心便是反复形式并求解修正方程式.牛顿法当初始 估讣值兀⑹和方程的精确解足够接近时，收敛速度非常快，具有平方收敛特性。

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

5/171

节点导纳矩阵及潮流计算

牛顿法的可靠收敛取决于有一个良好的启动初值。如果初值选择不当，算法 有可能根本不收敛或收敛到一个无法运行的节点上。对于正常运行的系统，各节 点电压一般均在额定值附近，偏移不会太大，并且各节点间的相位角差也不大， 所以对各节点可以釆用统一的电压初值(也称为平直电压)，如假定：

t// = 1 (2—8)

0>

印。>=0 或 = 1 Z

，o>

=0 (山\"2,…,MHf)

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

2. 2. 2牛顿一拉夫逊法潮流求解过程介绍

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

丿 a

=

局-B,,/j)+£2(G/+B 局)=o

(2-9)

对PV节点来说，给定量是R和匕，因此可以列出 3 =代-注工(G/+B 局)= o]

r

化-&+/>o

J

求解过程大致可以分为以下步骤:

—10)

6/176

节点导纳矩阵及潮流计算

（1）形成节点导纳矩阵；

（2）将各节点电压设初值U

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

（7） 检查是否收敛，如不收敛，则以各节点电压的新值作为初值自第3步 重新

开始进行狭义次迭代，否则转入下一步；

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

以直角坐标系形式表示： 错误!未定义书签。迭代推算式

采用直角坐标时，节点电压相量及复数导纳可表示为：

S,

岭=Gij 十 jBij

（2

—11）

将以上二关系式代入上式中，展开并分开实部和虚部;假定系统中的第1,

2,…，m号为P-Q节点，第m+1 , m+2,…，n—1为P—V节点，根据节点性 质的不

同，得到如下迭代推算式：

（1）于PQ节点

»

勺-

+场勺）

7 = 1,2,…,加 （2

—12）

7/177

节点导纳矩阵及潮流计算

7=1 ）=• .

⑵对于PV节点

8/171

节点导纳矩阵及潮流计算

乂胡-e 龙(G応-B』)-恥) .

、、

J = 〃？ + 1昇7? + 2,・・・，

>1

八一1

J-I

13)

(3)对于平衡节点

平衡节点只设一个，电压为已知，不参见迭代，其电压为：

vn

(2—14)

错误!修正方程

两组迭代式中包括2(n-l)个方程。选定电压初值及变量修正量符号 之后代入，并将

其按泰勒级数展开，略去M禺二次方程及以后各项，得

到修正方程如下:

(2—15)

△Q

A/；

△u =

A/；. rn+l△e

其中，

AW =

△匕

△W =—丿△(/

LA/;,-.

9/171

节点导纳矩阵及潮流计算

㈣

药

可'

SAP

CAP

AT

dSP

VV 迪

■ ■ oAP

AT

號 6边 號 • ■

H-I AT

oL

號 % 號

叱 ■ ■ 6△

沁Q\\ ■ ■

■ ■ • ■ %

6

叽

• • •

% 苏•

君 號 ■ ■ 號 6町

■ ••

•

• •

■ ■ ■ • • ■

0

叽

6△&叽

心

2—1 6 )

错误!未定义书签。 雅可比矩阵各元素的算式

式(2-12)中，雅可比矩阵中的各元素可通过对式(2-8)和(2-9)进行 偏导而求得。当丿V,时，雅可比矩阵中非对角元素为 沁 P

= = _* = _(% +吋) 8△勺 沁打 6曲_沁、_R … = =BiQ — Gif f： 6側沁j dSU2 d^U1 _n

------- = -------- =u

6SQ

7 ;

(2-17)

当•时，雅可比矩阵中对角元素为:

10

节点导纳矩阵及潮流计算

dSP A

=一2?©勺一 沁P 二

)-一 BJ

―^-= 一乞(q/j + B冲 一 Gafi + Bj©

—= S(©./} + ◎勺)-G/ + B眉

Ge

i

川

= _X(GQ 勺一 ，jfj)+ ai + ,ifi YJ j J-l

B

Ge

B

(2-18)

由式(2-13)和(2-18)看出，雅可比矩阵的特点：

\\ o \\ac (0,1)矩阵中各元素是节点电压的函数，在迭代过程中，这些元 素随着节点电

压的变化而变化；

错误!未定义书签。导纳矩阵中的某些非对角元素为零时，雅可比矩阵中对 应的元素也是为零.若岭=0,则必有Jit = 0 ;

错误!未定义书签。雅可比矩阵不是对称矩阵；0 = §1,2,...,“;/工巧 雅可比矩阵各元素的表示如下:

一(％ +吋)

〈j 右)

~X (©勺一场Z/)一G”£ - BJi (j = i)

(2—19)

3/ — G/)

八如

(丿知)

一工 Gjfj + By」)+ B启-Giifi (j = i) g

(2—2 0 )

11/1710

节点导纳矩阵及潮流计算

工(Gijfj + d勺)+ B启 一 GJ； (j = i)

(2—21)

12/1710

节点导纳矩阵及潮流计算

(M)

一工 9泸 j _ By fj) + Gj© 十 B” 九(j = i)

沁从 =0 (j 知)

6勺 i -2eQ = i)

%+吋)

(2-22)

(2-23)

3分析计算

1.根据给定的运行条件，确定图中所示电力系统潮流计算时各节点的类型

和待求量

根据图中可以看出各节点的类型和待求量分别为:

节点1： PQ节点待求量：SO

节点2: PV节点 待求量：Q、e

节点3：平衡节点 待求量：

2o求节点导纳矩阵Y

诊=治+儿3 = 1・25-丿55

^2=>12 + >23 = 1-3-77 匕3 =

儿3 +儿3=155-丿65

K,=^2=-O.5 + J3 y；2=K3=-O.8 + y4 $=帛=_0.75 + J2.5

所以节点导纳矩阵为：

13/1711

节点导纳矩阵及潮流计算

-0.5 + j3

-0.75 + J2.5 -0.8 + j4 1.55-J6.5

Y = -0.5 + j3

-0.75 + J2.5

1.3-J7 -0.8 + J4

3o 潮流方程或功率方程的表达式 因为对n个节点的网络，电力系统的潮流方程一

般形式是:

R—jQ

*

V；

其中匕二P&i — PL护Qx二QG・・一 QL「，即PQ分别为节点的有功功率无 功功率。

所以代入得潮流方程：

-^―A =( 1 .25—j 5.5) • t/.ZJ. + (0. 5—j 3 )・1上久+(0。 75-j2o 5) 120°

0・5-仏

1\"

=(0。 4) • 1/0°

5 — J3)・ qzq+(l。3-J7)・ 1Z① + ( 0. 8-J

厶

竺一上 =(0. 75-j2o 5) • U.Zd +(0・ 8 — J4) • 1ZJ. +(1 o 55 — 1Z09 ・

j6・ 5) • IZ0°

4.用牛顿-拉夫逊法计算潮流时，给岀修正方程和迭代收敛条件 (1) 修正方程

计算1、2节点的不平衡量存

14/1712

节点导纳矩阵及潮流计算

节点3是平衡节点，其电压、是给定的，故不参加迭代。 根据给定的容许误差£=lCf,按收敛判据

行校验，以上节点1、2的不平衡量都未满足收敛条件，于是继续以下计算.

修正方程式为 ZVSLX (n=3)

㈣

■

场 场

■

沁P]

d'P]

沁P] 沁P]

d^V；

■ 峯

5AV? 3AV；

■ ■ ■

以上雅可比矩阵J中的各元素值是通过求偏导数获得的，对PQ节点来说,

“是给定的，因而可以写出

对PV节点来说，给定量是只、和/訂因此可以列出

13/1713

节点导纳矩阵及潮流计算

当丿•幻时，雅可比矩阵中非对角元素为

B旳-GJ

当丿•“时，雅可比矩阵中对角元素为:

代入数值后的修正方程为

—1.25 -5L5 -^55 Q5 O

L25

Q5 3

3 \" -Q5

—2 —1 Q5 _o_

3 —1.3 O

—2

O _ 求解修正方程得

A/. -0.2547 ' -0.3611 0 _-0」015 （2）收敛条件

举=孝）+&。）=1-0254鼻 0T7453

16 / 1714

节点导纳矩阵及潮流计算

卅）=/+ 今严=O-Q36 14=-Q36 1 4° =閱+呈》=17=1

於）=理）+£» =O-Q 1 01S-Q1 01

—轮迭代结束，根据收敛条件收敛判据卜，若等

式成立，结果收敛，迭代结束，计算平衡节点的功率和线路潮流计算，否则继续 计算雅可比矩阵，解修正方程，直到满足收敛判据。

4结果分析

给定节点电压初值伞令*驾須二空去＞ 经过四次笔算

迭代过程后，得到节点电压和不平衡功率的变化情况分别于表4。1和表4。2所 示（取

£=1（/）:

迭代 计数k 节点电压 乙十+旗 认=色+必 1 0.7 4 5 3 —jOo 3611 1-jO. 1015 1 2 10 3 7 4 0 .413 1 —jO. 3 5 0 o 99 0 1-jO. 1 4 79 1 1.2973-j Oo 3 7 9 1 .0 0 8 3 - j 0. 018 5 1 Oo 8217-j Oo 3 4 0o 9986 —j 0。 1 088 0 表4o 1迭代过程中节点电压变化情况

17/1715

节点导纳矩阵及潮流计算

迭代 计数k 节点不平衡量 ■ —2 —Oo 14 8 2 0 1 -1 -Oo 9769 26 Oo 5 —0 o 07 03 -0.04 80 -0o 3610 0 -0. 0 1 2 -0 o 09 0 2 —0・6 0 7 1 —0 . 002 2 -0 o 017 1 3 -Oo 6 272 1 一4。325 4 -0.1816 -lo 2510 —0 1042 o 一 0 o 0049 表4. 2迭代过程中节点不平衡量变化情况

结果值与我的小组同学基本一样,也在预期之内.得到了基本一致的结果。 并且确定牛

顿法具有很好的二次收敛性，是求解多元非线性方程的正确算法。

5总结

这次的电力系统分析课程设讣让我对平时所学的专业知识有了更深刻更具体的 了解，明口了理论知识必须与实践相结合才能更好的发挥作用。在不停的翻书 上网查资料的过程中，我积累了大量的导纳矩阵和潮流计算以及电力系统的知识, 全面透彻的了解了相关知识的应用•使自己的知识更加牢固，并且有了更深的理 解。

通过这次毕业设计，我才明口学习是一个长期积累的过程，在以后的工作、 生活中都应该不断的学习，努力提高自己的知识和综合素质。总之，万事开头难, 知识必须通过应用才能实现其价值！有些东西以为学会了，但真正到用的时候才 发现是两回事，所以我认为只有到真正会用的时候才是真的学会了。

18/1716

节点导纳矩阵及潮流计算

参考资料

[1]何仰赞温増银.电力系统分析（第三版）。华中科技大学出版社 [2] 吴国炎。电力系统分析。浙江大学出版社

[3] 华智明岳湖山。电力系统稳态计算。重庆大学出版社

19/1717