基于Matlab电力系统潮流计算 - 图文

目 录

引言 ............................................................................................. 错误！未定义书签。 1 绪论 ........................................................................................................................... 1 1.1 电力系统叙述 ................................................................................................. 1 1.2潮流计算简介 .................................................................................................. 1 1.3 潮流计算的意义及其发展 ............................................................................. 2 2电力系统潮流计算基本原理 .................................................................................... 2 2.1电力网络的基本方程式 .................................................................................. 2 2.1.2 自导纳和互导纳的确定方法 ............................................................... 4 2.1.3 节点导纳矩阵的性质及意义 ............................................................... 5 2.1.4 非标准变比变压器等值电路 ............................................................... 6 2.2潮流计算的数学模型 ...................................................................................... 7 2.2.1 潮流计算的节点类型 ......................................................................... 7 2.2.2 潮流计算基本方程 ............................................................................. 7 2.3潮流计算方法 .................................................................................................. 9 2.3.1 牛顿——拉夫逊法 ............................................................................... 9 2.3.2 高斯——赛德尔法 ............................................................................... 9 2.3.3 PQ分解法 .......................................................................................... 10 3牛顿拉夫逊潮流计算理论分析 .............................................................................. 11 3.1概述 ................................................................................................................ 11 3.2牛顿法基本原理 ............................................................................................ 11 3.3牛顿法潮流计算方程 .................................................................................... 15

I

3.3.1节点功率方程 ...................................................................................... 15 3.3.2 修正方程 ........................................................................................... 16 4基于matlab的实例仿真 ........................................................................................ 19 4.1潮流计算原始资料参数 ................................................................................. 19 4.2参数计算及等值电路的绘制 ........................................................................ 20 4.2.1节点设置及分类 .................................................................................. 20 4.2.3支路参数计算并求解 .......................................................................... 21 4.3求解方法 ........................................................................................................ 22 4.4牛顿–拉夫逊潮流计算法的求解过程 ........................................................ 23 4.4.1牛顿–拉夫逊潮流计算法的计算框图 .............................................. 23 4.4.2牛顿法计算潮流的步骤如下 .............................................................. 24 4.4.3利用已知MATLAB程序求解，并修改相应程序变量 ........................ 25 4.4.4变电所负荷为题目所给数据进行求解 .............................................. 25 4.4.5修改程序 .............................................................................................. 26 4.5运行matlab程序输出结果 .......................................................................... 27 4.6 matlab仿真结果分析 .......................................................................................... 28 5总结 .......................................................................................................................... 30 致 谢 ........................................................................................................................... 31 参考文献 ..................................................................................................................... 31

II

基于Matlab电力系统潮流计算

摘要:电力系统潮流计算是研究稳态运情况一种重要分析计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态。在电力系统规划设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性，可靠性和经济性。MATLAB使用方便，有着其他高级语言无法比拟的强大的矩阵处理功能。这样使MATLAB成为电力系统潮流计算的首选计算机语言。牛顿-拉夫逊法是电力系统潮流计算的常用算法之一，它收敛性好，迭代次数少。介绍了电力系统潮流计算机辅分析的基本知识及潮流计算牛顿-拉普逊法,最后介绍了利用matlab进行实例仿真。

关键词：电力系统潮流计算；牛顿-拉普逊法；matlab

1 绪论

1.1 电力系统叙述

电力工业发展初期，电能是直接在用户附近的发电站(或称发电厂)中生产的，各发电站孤立运行。随着工农业生产和城市的发展，电能的需要量迅速增加，而热能资源(如煤田)和水能资源丰富的地区又往往远离用电比较集中的城市和工矿区，为了解决这个矛盾，就需要在动力资源丰富的地区建立大型发电站，然后将电能远距离输送给电力用户。同时，为了提高供电可靠性以及资源利用的综合经济性，又把许多分散的各种形式的发电站，通过送电线路和变电所联系起来。这种由发电机、升压和降压变电所，送电线路以及用电设备有机连接起来的整体，即称为电力系统。

电力系统加上发电机的原动机(如汽轮机、水轮机)，原动机的力能部分(如热力锅炉、水库、原子能电站的反应堆)、供热和用热设备，则称为动力系统。

现代电力系统提出了“灵活交流输电与新型直流输电”的概念。灵活交流输电技术是指运用固态电子器件与现代自动控制技术对交流电网的电压、相位角、阻抗、功率以及电路的通断进行实时闭环控制，从而提高高压输电线路的输送能力和电力系统的稳定水平。新型直流输电技术是指应用现电力电子技术的最新成果，改善和简化变流站的造价等。

1.2潮流计算简介

电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压，各元件中流过的功率，系统的功率损耗等等。在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性。可靠性和经济性。此外，电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。所以潮流计算是研究电力系统的一种很重要和基础的计算。

电力系统潮流计算也分为离线计算和在线计算两种，前者主要用于系统规划设计和安排系统的运行

1

方式，后者则用于正在运行系统的经常监视及实时控制。

利用电子数字计算机进行电力系统潮流计算从50年代中期就已经开始。在这20年内，潮流计算曾采用了各种不同的方法，这些方法的发展主要围绕着对潮流计算的一些基本要求进行的。对潮流计算的要求可以归纳为下面几点：

（1）计算方法的可靠性或收敛性； （2）对计算机内存量的要求； （3）计算速度；

（4）计算的方便性和灵活性。

电力系统潮流计算问题在数学上是一组多元非线性方程式求解问题，其解法都离不开迭代。因此，对潮流计算方法，首先要求它能可靠地收敛，并给出正确答案。由于电力系统结构及参数的一些特点，并且随着电力系统不断扩大，潮流计算的方程式阶数也越来越高，对这样的方程式并不是任何数学方法都能保证给出正确答案的。这种情况成为促使电力系统计算人员不断寻求新的更可靠方法的重要因素。

1.3 潮流计算的意义及其发展

作为电力系统分析中最基本的计算，潮流计算能够在复杂电力系统故障和正常条件下对稳态运行状态进行分析计算。计算给定运行状态是其目标。潮流计算的基础是对电力系统的规划，对电力系统静、暂态稳定分析。潮流计算的结果可于稳态研究，安全估计或着最优潮流等。这些对其方法和模型有直接影响，牛-拉发是潮流计算主要方法。

在用数字计算机解统潮流问题的开始阶段，通常采取以节点导纳矩阵为基础的逐次代入法。此方法原理简单，占用计算机内存较小，适应50年代电子计算机制造水平和当时电力系统理论水平。它的缺点是收敛性比较差，系统规模扩大，迭代次数就会骤然上升，计算过程常常会出现不收敛现象。为了解决这个困扰，我们运用阻抗法。它改善了计算过程中收敛性问题，解决了导纳法无法求解的一些系统的潮流计算，在60年代应用广泛。占用计算机内存大是它的主要不足，随着迭代次数变大，系统扩大时，这些缺点就更加明显。60年代中期以阻抗矩阵为基础的分块阻抗法很好的克服了它的缺点。这中方法把大系统分割成几个小系统，这样大大减少了计算机内存，同时也提高了计算速度。

此外还有一种方法也可以很好的解决上面出现的各种问题。牛—拉发是导纳矩阵为基础，因此只要我们在计算过程中保持矩阵的稀疏性，就能提高这种方法的

2电力系统潮流计算基本原理

2.1电力网络的基本方程式

电力网络可以用结点方程式或回路方程式表示出来。在结点方程式中表示网络状态的变量是各节点的电压，在回路方程式中是各回路中的回路电流。

一般若给出网络的支路数b，结点数n，则回路方程式数m为m=b-n+1结点方程式数m?为m?=n-1因此，回路方程式数比结点方程式数多d=m-m?=b-2n+2

在一般电力系统中，各结点(母线)和大地间有发电机、负荷 、线路电容等对地支路，还有结点和结点之间也有输电线路和变压器之路，一般b>2n，用结点方程式表示比用回路方程式表示方程式数目要少。

2

而且如以下所示，用结点方程式表示容易建立直观的方程式，输电线的连接状态等变化时也很容易变更网络方程式。

把电力系统的发电机端子和负荷端子（同步调相机等的端子也作为发电机端来处理）抽出来，剩下的输电线路及其它输电系统概括为网络Ｎet表示 。在发电机结点和负荷结点上标出任意顺序的记号：1,2,…,I,…,n.在输电系统Net的内部不包含电源，并且各节点和大地间连接的线路对地电容、电力电容器等都作为负荷来处理。

?1,v?2,?,v?n，由各端子流向输电系统Net的电流相应为令端子1,2……,n的对地电压分别为v?,I?,?,I?，则此网络方程组可以表示为 I12n??YV?????I1111?Y12V2???Y1kVk??Y1nVn????YV????I2?Y21V1222???Y2kVk??Y2nVn? (2-1) ??????I??????n?Yn1V1?Yn2V2???YnkVk???YnnVn(2-1)式可以简单写成

???YV?Iiijj (I=1,2,…,n) (2-2)

j?1n或者写成

I＝YV (2-3)

其中

?????V?I1??Y111???????VI??2? V?2 Y??Y21I?????????????????????In???Vn???Yn1?Y12?Y1n??Y22?Y2n? (2-4)

???Yn2?Ynn?(2-4)的Y称为节点导纳矩阵。因输电系统Net只是由无源元件构成的，而导纳矩阵是对称矩阵，于是有以下关系

Yij?Yji （2-5）

电压V和电流Ｉ的关系用式(2-1)～(2-5) 表示时称为节点导纳方程式。这里电压Ｖ用电流Ｉ的方程式表示时，则(2-3)式化为

V＝ZI (2-6)

其中 Z?Y

(2-6)式称为结点阻抗方程式，当然，阻抗矩阵也是对称矩阵。

?13

2.1.2 自导纳和互导纳的确定方法

电力网络的节点电压方程：

IB?YBUB (2-7)

式(2-7)IB为节点注入电流列向量，注入电流有正有负，注入网络的电流为正，流出网络的电流为负。根据这一规定，电源节点的注入电流为正，负荷节点

为负。既无电源又无负荷的联络节点为零，带有地方负荷的电源节点为二者代数之和。

式(2-7)UB为节点电压列向量，由于节点电压是对称于参考节点而言的，因而需先选定参考节点。在电力系统中一般以地为参考节点。如整个网络无接地支路，则需要选定某一节点为参考。设网络中节点数为（不含参考节点），则IB，UB均为n*n列向量。YB为n*n阶节点导纳矩阵。

节电导纳矩阵的节点电压方程：

IB?YBUB，展开为：

???Y11?I1?????I2??Y21????Y31?I ?3????????I????n??Yn1Y12Y22Y32?Yn2??Y13?Y1n??U1????Y23?Y2n??U2???Y33?Y3n??U3??? (2-8)

?????????Yn3?Ynn???Un?YB是一个n*n阶节点导纳矩阵，其阶数就等于网络中除参考节点外的节点数。 节点导纳矩阵的对

角元素Yii (i=1，2，?n)成为自导纳。自导纳数Yii值上就等于在i节点施加的电压，其余节点接地时，经过i节点注入的电流，因此，它可以定义为：

?/U?(U??0,j?i) （2-9） Yii?Iiij节点i的自导纳Yii数值上就等于与节点直接连接的所有支路导纳的总和。

节点导纳矩阵的非对角元素Yij (j=1，2，?，n;i=1，2，?。，n;j=i)称互导纳，由此可得互导纳Yij数值上就等于在节点i施加单位电压，其他节点全部接地时，经节点j注入网络的电流，因此可定义为：

?/U?(U??0,j?i) (2-10) Yji?Ijiij节点j，i之间的互导纳Yij数值上就等于连接节点j，i支路到导纳的负值。显然，恒Yij等于Yji。互导纳的这些性质决定了节点导纳矩阵是一个对称稀疏矩阵。而且，由于每个节点所连接的支路数总有一个限度，随着网络中节点数的增加非零元素相对愈来愈少，节点导纳矩阵的稀疏度，即零元素数与总元素的比值就愈来愈高。

4

2.1.3 节点导纳矩阵的性质及意义

节点导纳矩阵的性质：

（1）YB为对称矩阵，Yij=Yji。如网络中含有源元件，如移相变压器，则对称性不再成立。

YB对无接地支路的节点，（2）其所在行列的元素之和均为零，即 ?Yi,j?0,?Yj,i?0。

j?1i?1nn对于有接地支路的节点，其所在行列的元素之和等于该点接地支路的导纳。利用这一性质，可以检验所形成节点导纳矩阵的正确性。

（3）YB具有强对角性：对角元素的值不小于同一行或同一列中任一元素。

（4）YB为稀疏矩阵，因节点i ，j 之间无支路直接相连时Yij=0，这种情况在实际电力系统中非常普遍。矩阵的稀疏性用稀疏度表示，其定义为矩阵中的零元素与全部元素之比，即 S?Z/n2， 式

中Z 为YB中的零元素。S 随节点数n 的增加而增加：n=50，S可达92%；n=100，S 可达90%；n=500，S可达99%，充分利用节点导纳矩阵的稀疏性可节省计算机内存，加快计算速度，这种技巧称为稀疏技术。

节点导纳矩阵的意义：

YB是n*n阶方阵，其对角元素 Yii(i=1，2，----n)称为自导纳，非对角元素Yij(i，j=1，2，?n，

i?j)称为互导纳。将节点电压方程IB?YBUB展开为：

???Y??I111??????I2???Y21???????????In????Yn1??/U?(U?可见Yii?Iiij?Y12?Y??Y22n2???U???Y1n1????2??YU2n????????? （2-11）

??????Ynn???Un??0,i,j?1,2,?,n,i?j) (2-12)

表明，自导纳Yii在数值上等于仅在节点i施加单位电压而其余节点电压均为零（即其余节点全部接地）时，经节点i注入网络的电流。其显然等于与节点i直接相连的所有支路的导纳之和。同时可见

?/U?(U??0,i,j?1,2,?n,j?i)。表明，互导纳在数值上等于仅在节点j施加单位电压而其Yij?Iiji余节点电压均为零时，经节点i注入网络的电流，其显然等于(?yij)即Yij=?yij。yij为支路的导纳，负号表示该电流流出网络。如节点ij之间无支路直接相连，则该电流为0，从而Yij=0。

注意字母几种不写法的不同意义：粗体黑字表示导纳矩阵，大写字母Yij代矩阵YB中的第i行第j列元素，即节点i和节点j之间的互导纳。小写字母i，j支路的导纳等于支路阻抗的倒数数，

5

yij?1/Zij。

2.1.4 非标准变比变压器等值电路

变压器型等值电路更便于计算机反复计算,更适宜于复杂网络的潮流计算.双绕组变压器可用阻抗与一个理想变压器串联的电路表示.理想变压器只是一个参数,那就是变比??U1/U2。现在变压器阻抗按实际变比归算到低压侧为例,推导出变压器型等值电路。

图2-2双绕组变压器原理图

图2-3变压器阻抗归算到低压侧等值模型

流入和流出理想变压器的功率相等

?I??? U11?U1I2/K

? I1?/K (2-13) ?I2式(2-13)中, ??U1/U2是理想变压器的变比,U1和 U2分别为变压器高,低绕组的实际电压.从图2-3直接可得：

?K?U??I?Z （2-14） U122T????UUYUYU12T1T2I1?2??? 从而可得: ? 2?ZT?ZT?????UUYTU121?? （2-15） I2????YTU2?ZTZT?式（2-14）中YT?1/ZT,又因节点电流方程应具有如下形式:

???I1?Y11U1+Y12U2

??YU??-I2211+Y22U2 （2-16）

将式（2-14）与（2-15）比较，得：Y11=YT/?，Y12=-YT/?；Y22=YT。 21=-YT/? ，Y6

2因此可得各支路导纳为：

????1??Y10?Y11?Y12?YT? （2-17） 2?????1Y20?Y22?Y21?YT???由此可得用导纳表示的变压器型等值电路：

Y12=-Y12?YT/?Y21=-Y21?YT/?

图2-4变压器型等值电路

2.2潮流计算的数学模型

2.2.1 潮流计算的节点类型

根据电力系统中各节点性质的不同，很自然地把节点分成三类： （1） PQ节点

对于这类节点，等值负荷功率PLi 、QLi和等值电源功率PGi、QGi是给定的，从而注入功率Pi、Qi

也是给定的，待求的则是节点电压的大小Ui和相位角δi。这类节点称为PQ节点 （2） PU节点

对于这类节点，等值负荷和等值电源的有功功率PLi 、QGi，从而注入有功功率Pi是给定。等值负荷的无功功率QLi和节点电压Ui的大小事给定。待求的则是等值电源的无功功率QGi，从而注入无功功率Qi和节点电压相位角δi。这类节点称为PU节点 （3） 平衡节点

潮流计算时，一般都只设一个平衡节点。对这个节点，等值负荷功率PLS、QLS是给定的，节点电压大小和相位角Us、δs.也是给定的。待求的则是等值电源功率PGS、QGS，从而注入功率Ps、Qs。担负调整系统频率任务的电厂母线往往被选作平衡节点。 2.2.2 潮流计算基本方程

采用导纳矩阵时，节点注入电流和节点电压构成如式(2-7)所示线性方程组可展开如下形式：

???YV? Ii ij j j ?1n (i?1,2,?n ) (2-18)

由于实际电网中测量的节点注入量一般不是电流而是功率，因此必须将式中的注入电流用节点注入功率来表示。

7

节点功率与节点电流之间的关系为：

Si=Pi式中Pi?I??jQi?Uii (2-19)

?PGi?PLDi，Qi?QGi?QLDi

Pi?jQi

?Ui??S/U??因此用导纳矩阵时，PQ节点可以表示为Ii?ii把这个关系代入式中 ，得

nPi?jQi??YijUj(i?1,2,?n) （2-20）

?Uij?1式（2-20）就是电力系统潮流计算的数学模型-----潮流方程。它具有如下特点： 1：它是一组代数方程，因而表征的是电力系统的稳定运行特性。 2：它是一组非线性方程，因而只能用迭代方法求其数值解。

3：由于方程中的电压和导纳既可以表为直角坐标，又可表为极坐标，因而潮流方程有多种表达形式---极坐标形式，直角坐标形式和混合坐标形式。

?（1）取 Ui?Ui??i ，Yij?|yij|??ij，得到潮流方程的极坐标形式：

Pi?jQi?Ui??i?YijUj??i (2-21)

j?1n?(2) 取 Ui?ei?jfi， Yij?Gij?jBij，得到潮流方程的直角坐标形式：

?Pi?ei?(Gijej?Bijfj)?fi?(Gijfj?Bijej)?j?1j?1?? (2-22) nnQi?fi?(Gijej?Bijfj)?ei?(Gijfj?Bijej)??j?1j?1?nn?(3) 取 Ui?Ui??i Yij?Gij?jBij，得到潮流方程的混合坐标形式：

n??Pi?Ui?Uj(Gijcos?ij?Bijsin?ij)?j?1?? (2-23) n??U(Gsin??Bcos?)?Qi?Uijijijijij?j?1?不同坐标形式的潮流方程适用于不同的迭代解法。例如：利用牛顿---拉夫逊迭代法求解，以直角坐标和混合坐标形式的潮流方程为方便；而P-Q解耦法是在混合坐标形式的基础上发展而成，故当然采用混合坐标形式。

4: 它是一组n个复数方程，因而实数方程数为2n个但方程中共含4n个变量：P，Q，U和?，i=1，

8

这是变量?x1(0),?x2,?,?xn(0)(0)的线性方程组，称为牛顿法的修正方程，通过它可以解出

(0)(0)(0)?x1,?x2,?,?xn，并可以进一步求得

x1(1)?x1(0)??x1(0)?(1)(0)(0)?x2?x2??x2?? （3-14）

??(1)(0)(0)?xn?xn??xn?式中x1,x2,?,xn等到更接近真解的x1(2)(1)(1)(1)向真正解逼近了一步，如果再以它们作为初值重复解式（3-13）修正方程式，

,x2(2),?,xn(2)，如此迭代下去，并按式（3-14）进行修正，直到满足收敛要求

为止并停止迭代计算，这就构成了牛顿法的迭代过程。

一般第t次迭代式的修正方程为

??f1???x1(t)(t)(t)?f1(x1,x2,?,xn)??????f2(t)(t)(t)f(x,x,?,x)?21????x2n????1????(t)(t)(t)??fn(x1,x2,?,xn)?????fn???x1上式可以简写为

?f1?f??1??x2t?xnt?t(t)???x1??f2?f?(t)??2???x2? （3-15） ??x?x2nttt??????(t)???xn?????fn?fn???x?xnt?2tt?F(X(t))?J(t)?X(t) (3-16)

其中

(t)(t)?f1(x1(t),x2,?,xn)???(t)(t)(t)f(x,x,?,x)?，(t)212nF(X(t))?????J??(t)(t)(t)??fn(x1,x2,?,xn)????f1???x1???f2???x1??????fn???x1?f1?f??1??x2t?xnt?t??f2?f2???x?xnt? 2tt????fn?fn???x?xnt?2tt?其中的J(t)为第t次迭代时的雅克比矩阵；

同理可以得到第t次迭代时的修正量：

14

?X(t)??x1(t)??(t)??x2? （3-17） ???????(t)???xn??同样，也可以写出类似（3-14）的算式

X(t?1)?X(t)??X(t) （3-18）

这样反复交替的解式（3-16）及式（3-18）就可以使X(t?1)逐步趋近方程式的真正解。当满足人为收敛条件时，即

(t)(t)maxfi(x1(t),x2,?,xn??1或max?xi(t)??2 （3-19）

????迭代结束，式中?1,?2为预先给定的小正数。

3.3牛顿法潮流计算方程

3.3.1节点功率方程

电力系统的负荷习惯用功率表示，对于有n个节点的电力系统，系统中各节点注入电流与注入功率以标幺值表示的关系为

??S?IU?i?(Pi?jQi)*Ui* i=1，2，??，n （3-20）

式中?表示其共轭复数。将此关系式代入节点电压方程的通式，可得到以节点注入功率表示的节点电压方程:

?(P?jQ)ii

?Ui? (3-21) ??YijUjj?1n 上述的方程式，通常称为功率方程。根据方程中的节点电压向量表示的不同，可以得到不同形式的功率方程。

若节点电压向量以直角坐标表示，即以复数平面上实轴与虚轴上的投影表示可写成

Ui?ei?jfi （3-22）

其共轭值为

Ui?ei?jfi （3-23）

导纳表示为

Yij?Gij?jBij （3-24）

把这两关系式代回式（3-21）的功率方程中，展开后再将功率方程的实部和虚部分别写成有功、无

15

?功功率分离的节点方功率方程：

?Pi?ei?(Gijej?Bijfi)?fi?(Gijfj?Bijej)?j?1j?1?? （3-25） nnQi?fi?(Gijej?Bijfi)?ei?(Gijfj?Bijej)??j?1j?1?式中：i=1，2，??，n为各节点的编号。

nn??Uej?i 若节点电压以极坐标表示，则Uii或写成

??Ucos??jUsin? （3-26） Uiiiii 将其同导纳的复数表达式一起代入式（3-21）的功率方程，进整理可以得到

?Pi?Ui?Uj(Gijcos?ij?Bijsin?ij)?j?1?? （3-27） nQi?Ui?Uj(Gijsin?ij?Bijcos?ij)??j?1?式中：?ij??i??j——i与j节点电压的相角差。

由式（3-25）和（3-27）给出的功率方程表示方法避免了复数运算，因此，在潮流计算中普遍采用。 3.3.2 修正方程

采用牛顿法计算潮流时，需要对功率方程进行修改。下面将根据在不同坐标内的修改进行讨论： （1）在直角坐标系内时，由PQ节点功率方程（3-25）可知：节点i的注入功率是各点电压的函数，设节点的电压已知，代入式（3-25），可以求出节点i的有功及无功功率Pi,Qi，它们与给定的PQ 节点的注入功率Pis,Qis的差值应满足以下方程

n??Pi?Pis?Pi?Pis?ei?(Gijej?Bijfi)?fi?(Gijfj?Bijej)?0?j?1j?1?? （3-28） nn?Qi?Qis?Qi?Qis?fi?(Gijej?Bijfi)?ei?(Gijfj?Bijej)?0??j?1j?1?

对于PV 节点，已知节点的注入有功功率及节点电压大小，记作Pis,Uis，其节点的有功功率应满方程：

nn??Pi?Pis?Pi?Pis?ei?(Gijej?Bijfi)?fi?(Gijfj?Bijej)?0?j?1j?1? （3-29） ?2222?Ui?Uis?(ei?fi)?0?(i?m?1,m?2,.......,n?1)对于平衡节点，因为其电压给定，故不需要迭代求解。

通过以上分析可见，式（3-28）和式（3-29）共2（n-1）个方程，待求量e1,f1,e2,f2,?,en?1,fn?1共2（n-1）个。将上述2（n-1）个方程按泰勒级数展开，并略去修正量的高次方项后得到修正方程如下：

16

nn ?W??J?U （3-30）

22?W??P1?Q1??Pm?Qm?Pm?1?Um ?1??Pn?1?Un?1??T?U???e1?f1??em?fm?em?1?fm?1??en?1?fn?1?T??P??P??P??P??P??P??P???P?11111111????e??f1?em?fm?em?1?fm?1?en?1?fn?11?????Q1???Q1??Q1??Q1??Q1??Q1??Q1??Q1????e??f?e?f?e?f?e?f11mmm?1m?1n?1n?1?????????P??P??P??P??P??P??P??Pmmmmmmmm??????e1??f1?em?fm?em?1?fm?1?en?1?fn?1????Q??Q??Q??Q??Q??Q??Q??Q??mmmmmmmm????e?f1?em?fm?em?1?fm?1?en?1?fn?1?1?J????Pm?1??Pm?1??Pm?1??Pm?1??Pm?1??Pm?1???Pm?1??Pm?1?????e??f1?em?fm?em?1?fm?1?en?1?fn?11??22222222???Um??Um?1??Um?1??Um?1??Um?1??Um?1??Um?1??Um?1??1?????e?f?e?f?e?f?e?f11mmm?1m?1n?1n?1?????????P??P??P??P??P??P??P??P??n?1n?1n?1n?1n?1n?1n?1n?1????e?f1?em?fm?em?1?fm?1?en?1?fn?1?1??其中雅克22222222??Un?1??Un?11??Un?1??Un?1??Un?1??Un?1???Un?1??Un?1?????e??f1?em?fm?em?1?fm?1?en?1?fn?11??比矩阵的各元素可以对式（3-28）和式（3-29）求偏导数获得。

对于非对角元素（i?j）有

???Pi???Qi???(Gijei?Bijfi)??ej?fj????Pi??Qi???Bijei?Gijfi? （3-31） ?fj?ej????Ui2??Ui2???0?ej?fj??对于对角元素（i?j)有

17

n??Pi????(Gijej?Bijfj)?Giiei?Biifi??eij?1?n???Pi???(Gijfj?Bijej)?Giifi?Biiei??fij?1?n???Qi??(Gijfj?Bijej)?Biiei?Giifi??eij?1?? （3-32） n??Qi???(Gijej?Bijfj)?Giiei?Biifi???fij?1???Ui2???2ei??ei?2???Ui??2fi??fi?由上述表达式可以看到，雅克比矩阵具有以下特点：

1） 各元素是各节点电压的函数，迭代过程中每迭代一次各节点电压都要变化，因而各元素每次也变化；

2） 雅克比矩阵不具有对称性；

??Ui2??Ui23） 互导纳Yij?0，与之对应的非对角元素亦为零，此外因非对角元素??0，故雅克

?ej?fj比矩阵是稀疏矩。

（1） 当在极坐标系内时，由功率方程（3-27）可知节点i的注入功率是各节点电压幅值和相

角的函数。代入式（3-27）可以求出节点i的有功功率和无功功率，它们与给定的PQ节点的注入功率Pis,Qis的差值满足下面方程：

??Pi?Pi?Ui?Uj(Gijcos?ij?Bijsin?ij)?0?j?1? ? （3-33） n?Qi?Qi?Ui?Uj(Gijsin?ij?Bijcos?ij)?0??j?1?式中：?ij??i??j——i与j节点电压的相角差。

在有n个节点的系统中，假定第1~m号节点为PQ节点，第m+1~n-1号节点为PV节点，第n号节点为平衡节点。Vn和?n是给定的，PV节点的电压幅值Vm?1~Vn?1也是给定的，因此，只剩下n-1个节点的电压相角?1,?2,?,?n?1和m个节点的电压幅值V1,V2,?Vm是未知量。由（3-33）可知一共包含了n-1+m方程式，正好同未知量的数目相等，而直角坐标形式的方程少了n-1-m个。

由方程（3-33）可以写出修正方程

n??P?????Q????K???式中

?N??????1? (3-34) L????VD2?V??18

??P1???Q1????1????P???Q??????22?;?Q???;????2???P????????????????????P?Q???n?1??m??n?1??? （3-35）

??V1??V1????V??V??2?2?????V?;V????D2??????????VVm??m???其中：H是(n?1)?(n?1)阶方阵，其元素为?ij?K是m?(n?1)阶矩阵，其元素为Kij???Pi??Pi；N是(n?1)?m阶矩阵，其元素为Nij?Vj；

?Vi??j??Qi??Qi；L是m?m阶矩阵，其元素为Lij?Vj。

?Vj??j 对式（3-33）求偏导数，可得雅克比矩阵元素的表达式如下： 非对角元素（i?j）

Hij??ViVj(Gijsin?ij?Bijcos?ij)??Nij??ViVj(Gijcos?ij?Bijsin?ij)?? （3-36）

Kij?ViVj(Gijcos?ij?Bijsin?ij)?Lij??ViVj(Gijsin?ij?Bijcos?ij)??对角元素（i?j）

Hij?Vi2Bii?Qi??2N??ViGii?Pi? ii? （3-37） 2Kii?ViGii?Pi??Lii?Vi2Bii?Qi?4基于matlab的实例仿真

4.1潮流计算原始资料参数

系统地理接线如图1所示

线长为：1－火厂：70km 2－火厂：120km4－火厂：80km 2－3：100km 3－水厂：60km 3－火厂：200km均为双回线，电压为220KV，单回线路参数为：r1=0.085 Ω/km s/km

19

x1=0.32Ω/km b1=3.5*10

-6

图1 系统地理接线图

四个变电所的负荷cosФ=0.85,年初负荷分配如下：①70MW ②90MW ③60MW ④80MW 10KV 35KV 10KV 35KV

变压器配置：P=60MW 时，配二台50MVA变压器；每台PK=254KW，uk%=14.5

P=70 MW 时，配二台63MVA变压器；每台PK=245KW，uk%=12.5 P=80 MW 时，配二台63MVA变压器；每台PK=245KW，uk%=12.5 P=90 MW 时，配二台75MVA变压器；每台PK=414KW，uk%=16.7

要求：①、变电所的负荷为年初数值，年末最大运行方式数值应比年初数值高10％年中最大运行方式数值应比年初数值低3％ 各时期最小运行方式数值应比当时的最大运行方式数值低30％。②、火电厂作为平衡节点，水电厂作为PV节点。③、各节点电压应在UN－1.05UN之间。④、水厂夏季满发P=120MW，－40

系统的网络接线图如下图所示：

图2 系统网络接线图

4.2参数计算及等值电路的绘制

4.2.1节点设置及分类

根据电力系统图将火电厂母线设为节点节点1，水电厂节点设为节点8，将变电所1、2、3、4的高压侧分别设为节点2、5、6、9、低压侧为3、4、7、10。并且，将节点1设为平衡节点，将节点8设为PV节点，其余节点设为PQ节点，即：

20

?1—平衡节点?节点分类标号igl??2、3、4、5、6、7、9、10—PQ节点

?8—PV节点?

平衡节点：1（电压初始设置为231V）

PV节点：2 （电压初始设置为231V，由于总负荷为200MW电厂1和电厂2装机容量比约为3:2，所

以设置电厂2的发电功率为120MW）

PQ节点：2、3、4、5、6、7、9、10（电压设置为220V）。 电压基准值为220V，视在功率基准值为100WVA。 4.2.3支路参数计算并求解

设定电压基准值为220V，视在功率基准值为100WVA，根据题目原始资料，计算发电厂、变压器及线路的参数有名值。

参数计算公式如下：

1.变压器阻抗：

2PkUN?3 RT??102SN2.变压器电抗：

2UK%UN XT?100SN3.线路电阻：

R?r?L

4.电抗：

X?x?L

5.线路电纳：

B=b·L

6.变电所负荷分别为：

变电所1 SL=47.53+j29.456 变电所2 SL=61.11+j37.873 变电所3 SL=40.74+j25.248 变电所4 SL=54.32+j33.665

将参数整理，双回路支路阻抗除以2，对地电纳乘以2，见下表：

表1 各支路等值参数

首端号

末端号 支路的阻抗（R+jx）/? 21

支路的对地电纳B/S 1 2 1 5 5 6 6 1 1 2 3 5 4 6 8 7 6 9 2.975+j11.2 1.494+j48.016 10.2+j38.4 1.718+j53.885 8.5+j32 2.55+j9.6 2.459+j70.18 17+j64 3.4+j12.8 1.494+j48.016 j0.00049 0 j0.00042 0 J0.00035 j0.00042 0 J0.0007 j0.00056 0 9 10 7.计算变压器分接头：0.95 -1.05

变压器有5个抽头，电压调节范围为UN?2*2.5%，UN对应的分接头开始时设变压器高压侧接主接头,降压变压器5个分接头时的非标准变比k*以备调压时选用。

4.3求解方法

画出系统的等效电路图，在计算出各元件参数的基础上，应用牛顿—拉夫逊Newton-Raphson法以及MATLAB软件进行计算对给定系统图进行潮流计算，经过调节均得到符合电压限制及功率限制的潮流分布。

采用标么值进行计算,根据题中最高电压等级为220KV,所以选择基准电压及功率分别为：

2UBS?100MVA,UB?220KV,ZB??484?SB

据此可以计算各部分标幺值。

绘制等效电路时，双回线阻抗减半，容抗加倍；并联变压器阻抗变为一半，无容抗。 绘制的等效电路如下：

22

③1.494+j48.016Ω④②j0.000245S⑤1.718+j52.885Ωj0.000175S10.2+j38.4Ω2.975+j11.2Ωj0.000245j0.00021S8.5+j32Ωj0.000175S⑥2.55+j9.6Ω⑧j0.00021SG①17+j64Ωj0.00035S3.4+j12.8ΩGj0.00028Sj0.00035Sj0.00021S2.459+j70.18Ω⑦j0.00021Sj0.00028S⑨1.494+j48.016Ω⑩图3 等效电路图

4.4牛顿–拉夫逊潮流计算法的求解过程

4.4.1牛顿–拉夫逊潮流计算法的计算框图 根据给出的程序画出的潮流计算框图如下图：

23

支路损耗对比表如下：

表5 支路损耗比较表

支路 损耗 1-2 2-3 1-5 5-4 5-6 6-8 6-7 1-6 1-9 9-10 总和 未调节 0.1499 - 25.4039i 0.0943 + 3.0321i 0.0845 - 21.7464i 0.1955 + 5.9146i 0.489 - 16.479i 0.6933 - 19.7206i 0.1153 + 3.2901i 0.2039 - 36.4503i 0.2267- 28.7566i 0.125 + 4.0183i 调节适当后 0.1499 - 25.4039i 0.0943 + 3.0321i 0.0832 - 21.7542i 0.1848 + 5.592i 0.4886 - 16.4832i 0.6935 - 19.7206i 0.1153 + 3.2899i 0.2038 - 36.4517i 0.2267 -28.7566i 0.125 + 4.0183i 过调节 0.1499 - 25.4039i 0.0943 - 3.0321i 0.0819 - 21.7616i 0.1745 + 5.2811i 0.4882 - 16.4872i 0.6936 - 19.7206i 0.1153 + 3.2897i 0.2037 - 36.453i 0.2267 - 28.7566i 0.125 + 4.0183i 2.3534 -132.96i 2.3775 -132.3i 2.3652 -132.64i 对比以上参数可知，分接头越低，对应节点电压越高，同时损耗越小 表1形成的图像如下

图5 变压器分接头未调节为1时仿真图形

从图5可以直观的看出4号节点电压不在1-1.05范围之内，所以节点4不满足题目要求。

29

表2形成的图像如下

图6 变压器分接头调节为0.975时仿真图形

从图中可以看出各节点电压都在UN－1.05UN之间，因此满足题目要求。 表3形成的图像如下

图7 变压器分接头调节为0.95时仿真图形

综上所述在分接头为0.95和1时均超出允许范围，在分接头为0.975时满足电压在1-1.05UN电压在允许的范围内，满足题目要求的条件。

另外,在调节电压时，如果单独调节分接头不能输出满足条件时，可以同时调节发电机的电压，以使输出参数满足要求的条件。

5总结

潮流计算是电力系统的最基本、最常用的计算。根据系统给定的运行条件、网络接线及元件参数，

30

通过潮流计算可以确定各母线的电压（相角及幅值）、各元件中流过的功率、整个系统的功率损耗等。潮流计算是实现电力系统安全经济发展的必要手段和重要环节。

牛顿—拉夫逊法是一种经典的数学方法，用来解决矩阵的计算问题，在数值分析课程中我们也接触过这种方法，这次在电力系统的潮流计算中得到了应用，更让我觉得知识灵活运用的重要性。

MATLAB是一套高性能的数学计算软件，它集数值分析、矩阵计算、信号处理和图形显示于一身，构成了一个方便的界面和友好的用户环境，其强大的扩展功能为各个领域的应用提供了方便。在潮流计算的大量数值计算过程中更显示出其优势，是我们应该掌握的一门基本技术。在这次课程设计中，我不仅学会了matlab的使用方法，还学会了很多的实用的matlab指令。

在这次毕业设计中我还将平时学的调压方法运用在电力系统的调节，真正做到了学以致用。 实际的情况远比我们计算的情况复杂，这让我深刻了解了潮流计算的重要性。精准的潮流计算不仅可以节省能源并充分发挥能源的作用，还可以使电网处于稳定且平衡的状态，这对电网的安全运行起到关键性的作用。所以我认为学好电力系统的关键在于学好潮流计算。

致 谢

毕业设计是培养学生综合运用所学知识，发现，提出，分析和解决实际问题,锻炼实践能力的重要环节,是对学生实际工作能力的具体训练和考察过程。回顾此次毕业设计，从选题到定稿，从理论到实践，在近几个月的日子里，虽然工作辛苦，但是我掌握了很多新的东西，这都是在四年的理论学习中没有接触过的，不仅巩固了以前所学过的知识，而且拓展了本专业的知识面。通过这次毕业设计我懂得了理论与实际相结合的重要性，只有把所学的理论知识与实践相结合起来，从实验中得出结论，验证理论，应用理论，才能提高自己的实际动手能力和独立思考的能力。在设计的过程中不断遇到各种问题，正是这些问题使我发现了自己的不足之处，解决每一个问题的过程使我收获颇丰，也督促我还要继续深入学习。

这次毕业设计的顺利完成离不开刘玉娟老师的悉心指导，刘老师在设计方向以及很多细节方面都给予了我很大帮助，刘老师对待工作一丝不苟的严谨工作作风深深影响了我，特在此表示感谢！同时，对给过我帮助的所有同学和各位指导老师再次表示忠心的感谢！ 参考文献

[1]于永源, 杨绮雯，电力系统分析（第三版）［M]. 中国电力出版社，2007

[2] 张勇军， 任 震等， 基于灾变遗传算法的无功规划优化， 〔 J〕 电力系统自动化， 2002 年第 23 期。

[3]王守相，刘玉田 电力系统潮流计算研究现状--《山东电力技术》1996年05期

[4]何仰赞，温增银. 电力系统分析（上册）第三版［M]. 湖北：华中科技大学出版社，2002 [5] 刘同娟. MATLAB在电路分析中的应用.电气电子教学学报.2002 [6] 西安交通大学等.电力系统计算[M].北京：水利电力出版社，1993.12 [7] 李光琦.电力系统暂态分析[M].北京： 水利电力出版社，2002.5

[8]何仰赞，温增银. 电力系统分析（下册）第三版［M]. 湖北：华中科技大学出版社，2002

31

[9]韦化，李滨，杭乃善，等．大规模水一火电力系统最优潮流的现代内点算法实现[J]．中国电机工程学报，2003．23(6)：13一l8．

[10] 刘卫国.MATLAB 程序设计与应用 (第二版）[M].北京:高等教育出版社,2006. [11]王锡凡.现代电力系统分析[M].北京：北京科学出版社,2003.

[12]电力系统分析课程设计综合实验，祝淑萍 等，中国电力出版社，2001 [13]电力系统稳态分析，陈珩，中国电力出版社，2003

[14]《MATLAB实用指南》，苏金明，王永利，电子工业出版社，2005 [15] 孟祥萍，电力系统分析．高等教育出版社.2004

[16] 李维波，MATLAB在电气工程中应用，中国电力出版社，2007 [17] 邱晓燕，刘天琪，电力系统分析的计算机算法，中国电力出版社，2009

[18] 吴天明，赵新力，刘建存，MATLAB电力系统设计与分析，国防工业出版社，2007 [19] 张伯明,陈寿孙，高等电力网络分析，北京:清华大学出版社,1996 [20] 薛定宇，控制系统计算机辅助设计.北京:清华大学出版社，1996

[21] 刘军，MATLAB在电力系统分析中的应用，电力系统及自动化学报，2000 [22]刘卫国，陈昭平, 张颖，MATLAB程序设计与应用，北京:高等教育出版社，2002.6 [23]Tsai-Hsiang Chen, Mo-Shing Chen, Kab-ju Hwang, etc. Distribution system power flow analysis:arigid approach. IEEE Trans. On Power Delivery. 1991,6(30):1146-1152 [24]ZU Weiyu.The Proposed CLP Method for Allocation Real Power Losses of Multiple Flows.IEEE Power Engineering Review,2001.

[25] Wu Z Q, Chen G Z. MVA Power Flow and Loss Analysis for Electricity Market. IEE Proc on Ge ne r,Transm&Distrib, 2001, 148(2):153-158

32

Based on Matlab power system load flow calculation

Lu Jie

Department of Binjiang Auto Control, Nanjing University of Information Science & Technology, Nanjing,

210044

ABSTRACT

Power Flow Calculation of Power System is an important analysis and calculation of power system steady-state operation, which according to the given operating conditions and system wiring to determine the various parts of the power system running state. In the study of power system design and the current operation mode are required Power Flow Calculation to quantitatively analyzed and compared to the program or run mode power supply reasonable, reliability and economy or not.MATLAB is easy to use, the powerful matrix processing is the other high-level language can not be compared with. This allows MATLAB to become the preferred computer language of power flow calculation.Newton Raphson power flow calculation is one of the most commonly used algorithms, which has good convergence and fewer iterations .This article describes the power flow computer assisted analysis of the basic knowledge and power flow Newton - Raphson method, Finally it describe how to use M language to simulate .

Keywords: power flow calculation; Newton - Raphson method; matlab

33

本文来源：https://www.bwwdw.com/article/dqi7.html