电力系统潮流计算课程设计论文

课程设计论文

基于MATLAB的电力系统潮流计算

学院:电气工程学院

专业:电气工程及自动化

班级:电自0710班

学号:**********

**: ***

内容摘要

潮流计算是电力系统最基本最常用的计算。根据系统给定的运行条件，网络接线及元件参数，通过潮流计算可以确定各母线的电压（幅值和相角），各支路流过的功率，整个系统的功率损耗。潮流计算是实现电力系统安全经济发供电的必要手段和重要工作环节。因此，潮流计算在电力系统的规划计算，生产运行，调度管理及科学计算中都有着广泛的应用。

潮流计算在数学上是多元非线性方程组的求解问题，牛顿—拉夫逊Newton-Raphson法是数学上解非线性方程组的有效方法，有较好的收敛性。运用电子计算机计算一般要完成以下几个步骤：建立数学模型，确定解算方法，制订计算流程，编制计算程序。

关键词

牛顿-拉夫逊法（Newton-Raphson）变压器及非标准变比无功调节

高斯消去法潮流计算Mtlab

一 .电力系统潮流计算的概述

在电力系统的正常运行中，随着用电负荷的变化和系统运行方式的改变，网络中的损耗也将发生变化。要严格保证所有的用户在任何时刻都有额定的电压是不可能的，因此系统运行中个节点出现电压的偏移是不可避免的。为了保证电力系统的稳定运行，要进行潮流调节。

随着电力系统及在线应用的发展，计算机网络已经形成，为电力系统的潮流计算提供了物质基础。电力系统潮流计算是电力系统分析计算中最基本的内容，也是电力系统运行及设计中必不可少的工具。根据系统给定的运行条件、网络接线及元件参数，通过潮流计算可以确定各母线电压的幅值及相角、各元件中流过的功率、整个系统的功率损耗等。潮流计算是实现电力系统安全经济发供电的必要手段和重要工作环节，因此潮流计算在电力系统的规划设计、生产运行、调度管理及科学研究中都有着广泛的应用。它的发展主要围绕这样几个方面：计算方法的收敛性、可靠性；计算速度的快速性；对计算机存储容量的要求以及计算的方便、灵活等。

常规的电力系统潮流计算中一般具有三种类型的节点：PQ 、PV 及平衡节点。一个节点有四个变量，即注入有功功率、注入无功功率，电压大小及相角。常规的潮流计算一般给定其中的二个变量：PQ 节点（注入有功功率及无功功率），PV 节点（注入有功功率及电压的大小），平衡节点（电压的大小及相角）。 1、变量的分类：

负荷消耗的有功、无功功率——1L P 、1L Q 、2L P 、2L Q 电源发出的有功、无功功率——1G P 、1G Q 、2G P 、2G Q 母线或节点的电压大小和相位——1U 、2U 、1δ、2δ

在这十二个变量中，负荷消耗的有功和无功功率无法控制，因它们取决于用户，它们就称为不可控变量或是扰动变量。电源发出的有功无功功率是可以控制的自变量，因此它们就称为控制变量。母线或节点电压的大小和相位角——是受控制变量控制的因变量。其中, 1U 、2U 主要受1G Q 、2G Q 的控制, 1δ、2δ主要受

1G P 、2G P 的控制。这四个变量就是简单系统的状态变量。

为了保证系统的正常运行必须满足以下的约束条件：

对控制变量

max min max min ;Gi Gi Gi Gi Gi Gi Q Q Q P P P <<<<

对没有电源的节点则为

0;0==Gi Gi Q P

对状态变量i U 的约束条件则是

max min i i i U U U <<

对某些状态变量i δ还有如下的约束条件 max

j i j i δδδδ-<-

2、节点的分类：

⑴ 第一类称PQ 节点。等值负荷功率Li P 、Li Q 和等值电源功率Gi P 、Gi Q 是给定的，从而注入功率i P 、i Q 是给定的，待求的则是节点电压的大小i U 和相位角i δ。属于这类节点的有按给定有功、无功率发电的发电厂母线和没有其他电源的变电所母线。

⑵ 第二类称PV 节点。等值负荷和等值电源的有功功率Li P 、Gi P 是给定的，从而注入有功功率i P 是给定的。等值负荷的无功功率Li Q 和节点电压的大小i U 也是给定的。待求的则是等值电源的无功功率Gi Q ，从而注入无功功率i Q 和节点电压的相位角i δ。有一定无功功率储备的发电厂和有一定无功功率电源的变电所母线都可以作为PV 节点；

⑶ 第三类平衡节点。潮流计算时一般只设一个平衡节点。等值负荷功率Ls P 、Ls Q 是给定的，节点电压的大小和相位也是给定的。担负调整系统频率任务的发电厂母线往往被选作为平衡节点。

二．牛顿—拉夫逊法概要

1．首先对一般的牛顿—拉夫逊法作一简单的说明。 已知一个变量X 函数为：0)(=X f

到此方程时，由适当的近似值)

0(X

出发，根据：

,......)2,1()

()()()()

()

1(='-=+n X f X f X

X

n n n n 反复进行计算，当)

(n X 满足适当的收敛条件就是上面方程的根。这样的方法

就是所谓的牛顿—拉夫逊法。

这一方法还可以做下面的解释，设第n 次迭代得到的解语真值之差，即)

(n X 的

误差为ε时，则：

0)()(=+εn X f

把)()

(ε+n X

f 在)

(n X

附近对ε用泰勒级数展开

......)(!

2)()()()(2

)

()

()

(=+''+

'+=+n n n n X f X

f X f X f εεε

上式省略去2ε以后部分

0)()()()(≈'+n n X f X f ε

)

(n X 的误差可以近似由上式计算出来。

)

()

()()(n n X f X f '-

≈ε 比较两式，可以看出牛顿—拉夫逊法的休整量和)

(n X 的误差的一次项相等。

用同样的方法考虑，给出n 个变量的n 个方程：

⎪⎪⎩⎪⎪⎨

⎧===0

),,,(0),,,(0),,,(21212211n n n

n X X X f X X X f X X X f 对其近似解1X '得修正量1

X '∆可以通过解下边的方程来确定： ⎥⎥⎥⎥⎦⎤

⎢⎢⎢

⎢⎣⎡∆∆∆⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤

⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂-=⎥⎥⎥⎥⎥⎥⎦

⎤⎢⎢⎢⎢⎢⎢⎣⎡'''''''''n n n n n n n n n n n X X X x f x f x f x f x f x f x f x f x f X X X f X X X f X X X f

212

122212

121

1

1

21212211),,,(),,,(),,,( 式中等号右边的矩阵

n

n

x f ∂∂都是对于n X X X ''',,,2

1 的值。这一矩阵称为雅可比（JACOBI ）矩阵。按上述得到的修正向量n X X X '∆'∆'∆,,,2

1 后，得到如下关系 n n n

X X X ∆+'='' 这比n X X X ''',,,21 更接近真实值。这一步在收敛到希望的值以前重复进行，

一般要反复计算满足