C语言进行潮流计算

电力系统课程设计

C语言潮流计算

学院：电气工程

班级：电092班

学号：0912002020

学生姓名：闵凯

2013.3.7

电力系统的潮流计算是对电力系统分析的最基本步骤也是最重要的步骤，是指在一

定的系统结构和运行条件下，确定系统运行状态的计算，也即是对各母线（节点）电压，各元件（支路）传输电线或功率的计算。通过计算出的节点电压和功率分布用以检查系统各元件是否过负荷，各点电压是否合理，以及功率损耗等。

即使对于一个简单的电力系统，潮流计算也不是一件简单就可以完成的事，其运算量很大，因此如果对于一个大的、复杂的电网来说的话，由于其节点多，分支杂，其计算量可想而知，人工对其计算也更是难上加难了。特别是在现实生活中，遇到一个电力系统不会像我们期望的那样可以知道它的首端电压和首端功率或者是末端电压和末端功率，而是只知道它的首端电压和末端功率，更是使计算变的头疼万分。为了使计算变的简单，我们就可以利用计算机，用C 语言编程来实现牛顿-拉夫逊（Newton-Raphson ）迭代法，最终实现对电力系统潮流的计算。

一．

用牛顿-拉夫逊迭代法进行电力系统潮流计算的相关概念

1.节点导纳矩阵

如图所示的电力网络,将节点i 和j 的电压用•

U i 和•

.U j 表示,它们之间的支路导纳表示为y ij ,那么有基尔霍夫电流定律可知注入接点I 的电流•

.I i (设流入节点的电流为正)等于离开节点I 的电流之和,因此有

)(.

00•

•

≠=≠=•

•

-==∑∑j i n

i

j ij

n

i

j ij i U U I I y (1-1)

∴ •

≠=≠=•

•∑∑-=U y y U I n

i

j ij n

i

j ij i 00 (1-2)

如令

ii n

i

j ij

Y y

=∑≠=0

ij ij Y y =-则可将(1-2)改写为:

∑≠=•

•

=

n

i

j ij ij

i U Y

I 1 I=1,2,…,n. (1-3)

上式也可以写为: I =YU (1-4)

其中Y 为节点导纳矩阵，也称为稀疏的对称矩阵，它是n×n 阶方阵。对角元Y ii 称为自导纳，它等于与该节点I 直接相连的所有支路导纳总和；非对角元Y ij （i ≠j ）称为互导纳或转移导纳，它等于连结节点I ，j 支路导纳的负数，且有Y ij =Y ji ，当节点I ，j 之间没有支路直接相连时，Y ij =Y ji =0。

电力系统的分析计算中，往往要作不同的运行方式下的潮流计算，如果系统发生变化，如投切一条线路或一台变压器，由于改变了一条支路的状态或参数只影响该支路两端节点的自导纳和他们之间的互导纳，因而对每一种运行方式不必重新形成导纳矩阵，只需对原有导纳矩阵作相应的修改即可。

2．潮流计算的功率方程

在实际的电力系统中，已知的条件往往不是节点的注入电流而是负荷和发电机功率，且这些功率一般不随节点的电压变化而变化，而节点的电流则是随电压的变化而变化的，因此在已知节点导纳矩阵的情况下，必须用已知的节点功率来替代未知的节点注入电流，才能求

I i i δ∠

S ~Li jQ +

Li Gi i P P P -= （1-5）

Li Gi i Q Q Q -= （1-6） Θ ∑=•

•

=

n

j j ij

i U Y

I 1

(1-7)

节点注入电流用功率和电压表示为：

i

i

i i Li Gi Li Gi i i i U jQ P U Q Q j P P U S I ˆˆ)()(ˆˆ-=

---==•

（1-8） ∴

功率方程可以表示为：

∑=•=-n j j ij i i

i U Y U jQ P 1

ˆ （1-9）

3．节点分类

对于有n 个节点的电力网络，可以列出n 个功率方程，由图可知一个节点有四个变量：注入有功功率P i ，注入无功功率Q i ，节点电压幅值U i 和相角i δ。n 个节点有4n 个变量,但只有2n 个关系式,所以为了使潮流有确定解,必须给定其中2n 个变量。根据给定节点变量的不同，可以有以下三种类型的节点：

（1）

PQ 节点：给定注入功率P i ，Q i ，即已知P Gi ，P Li ，Q Gi ，Q Li ，待求U i ，δi 。例如：降压变电所母线（负荷节点），固定出力的发电厂母线。

（2） PV 节点：给定了注入有功功率P i （P Gi ，P Li ），U i 和Q Li ，待求Q Gi （Q i ），δi 。

例如：有一定无功电源的降压变电所母线，有一定储备的发电厂母线。

（3） 平衡节点：给定了U i ，δi 和P Li ，Q Li ，待求P Gi ，Q Gi ，即P i ，Q i ，用来平衡

全电网的功率，通常在一个独立的电力系统中只设一个平衡节点。

4．牛顿-拉夫逊迭代法

牛顿-拉夫逊迭代法将解非线性方程组的过程转化为反复求与之相对应的线性方程的求解过程。

对于一个n 维非线性方程组：i n i y x x x f =),...,,(21 n=1,2,3,…,n

假定其初值为x 1(0),x 2(0),…,x n (0),也即其近似解,它与真值之间的误差为

)

0()0(2)0(1,...,,n

x x x ∆∆∆也即各变量与真解之间的修正量。

将这n 个方程式都在初值的附近展开成Taylor 级数且忽略二次项及高次项,则可得修正方程

i n i y x x f =∆∂∂++∆∂∂+

⋯)

0(0

n n )0(1011(0)

n (0)2(0)1x f ...x f )x ,,x ,x (, I=1,2,…n. (1-10) 将修正方程写成矩阵形式:

(1-11)

其中令J= ⎥⎥⎥

⎥⎥⎥⎦

⎤

⎢⎢⎢⎢⎢⎢⎣⎡∂∂∂∂∂∂∂∂0n n 01n 0n 10

11x f ...x f .........x f ...x f ,称之为雅可比(Jacobi)方阵。

它的第I 行，第j 列交点的元素为第I 个函数),...,(21n i x x x f 对第j 个变量x j 的偏导数在点（x 1(0),x 2(0),…,x n (0)）的值，所以方程组是线性方程，可用于求出)

0()

0(2)

0(1,...,,n

x x x ∆∆∆，从

而得到新的近似解，

)

0()

0()

1(i i

i

x x x ∆+= （1-12）

于是得到一般迭代式：

⎥⎥

⎥

⎦

⎤

⎢⎢⎢⎣⎡∆∆⎥

⎥⎥⎥⎥⎥

⎦

⎤

⎢⎢

⎢⎢⎢⎢⎣⎡∂∂∂∂∂∂∂∂=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⋯-⋯-)0()0(10n n 01n 0n 10

11

(0)n (0)2(0)1(0)n (0)2(0)111...x f ...x f .........

x f ...x f )x ,,x ,x (:)x ,,x ,x (n n n x x f y f y