电力系统分析潮流计算课程序设计及其MATLAB程序设计

电力系统分析潮流计算课程序

设计及其MATLAB程序设计

，自YN大多

广州学院

GuangzhouCollegeofSouthChinaUniversityofTechnology

电力系统分析潮流计算程序设计报告

题目：13节点配电网潮流计算

学院电气工程学院

专业班级_______________________________

学生姓名_______________________________

学号_______________________________

班内序号______________________________

指导教师____________房大中___________

提交日期2015年05月04日

目录

一、程序设计目的.............................2

二、程序设计要求.............................4

三、13节点配网潮流计算......................4

3.1主要流程........错误!未定义书签。

第一步的前推公式如下（式如-（1-5）：

....................错误!未定义书签。

第二步的回代公式如下（1-6）—（1-9）：

....................错误!未定义书签。

3.2配网前推后代潮流计算的原理......8

3.3配网前推后代潮流计算迭代过程………8

3.3计算原理.........................9

犯、计算框图流程............................11

五、确定前推回代支路次序….•…错误!未定义书签。

六、前推回代计算输.诩入文件12

主程序：....12

输入文件清单:15

计算结果：一17

据分析:18

七、配电网潮流计算的要点18

八、自我总结............18

0

九、参考文献..............................19

附录一MATLAB的简介......................20

1

华南理工大学广州学院电气工程学院课程设计报告

一、程序设计目的

开式网络潮流计算：配电网的结构特点呈辐射状，在正常运行时是开环的；配电网

的潮流计算采用的方法是前推回代法，本程序利用前推回代法的基本原理、收敛性。

(1)在电网规划阶段，通过潮流计算，合理规划电源容量及接入点，合理规划网架,

选择无功补偿方案，满足规划水平年的大、小方式下潮流交换控制、调峰、调相、调压

的要求。

(2)在编制年运行方式时，在预计负荷增长及新设备投运基础上，选择典型方式进

行潮流计算，发现电网中薄弱环节，供调度员日常调度控制参考，并对规划、基建部门

提出改进网架结构，加快基建进度的建议。

(3)正常检修及特殊运行方式下的潮流计算，用于日运行方式的编制，指导发电厂

开机方式，有功、无功调整方案及负荷调整方案，满足线路、变压器热稳定要求。

及电压质量要求。

(4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方

案。

1o

O--------OO--------O--------O

76289

O------O------O------O

1()31112

Q------O

50

图113节点配电网结构图

表1系统支路参数

2

华南理工大学广州学院电气工程学院课程设计报告

支路R(。)X(Q)B/2(TK)(S)

1-23.3673.6850.0

2~32.3562.5410.0

3-41.1451.280.0

4~54.5245.040.0

2~60.8561.140.0

6~72.7452.9650.0

2〜83.7434.2510.0

8-92.2372.7560.0

3-104.144.6960.0

3〜111.3281.7630.0

11-122.4362.8660.0

4~133.5213.9660.0

表2系统负荷参数

节点编号节点类型节点初始电(kV)Pi(MVA)Qi(MVA)

1根节点10.400

2中间节点10.00.03420.0301

3中间节点10.00.06930.0642

4中间节点10.00.08450.0763

5叶节点10.00.02950.0261

6中间节点10.00.04740.0409

7叶节点10.00.11760.0957

8中间节点10.00.09460.0857

9叶节点10.00.09160.0859

10叶节点10.00.02710.0229

11中间节点10.00.06960.0643

12叶节点10.00.06760.0579

13叶节点10.00.02980.0242

3

华南理工大学广州学院电气工程学院课程设计报告

二、程序设计要求

1.看懂前推回代法计算程序；

2.报告叙述计算原理及计算流程；

3.绘制计算流程框图；

4.确定前推回代支路次序(广度优先，或深度优先)并编写前推回代计算

输入文件，然后进行潮流计算；

5.整理专利计算结果；

6.总结配电网潮流计算的要点；

三、设计内容

1.根据电力系统网络推导电力网络数学模型，写出节点导纳矩阵；

2.赋予各节点电压变量(直角坐标系形式)初值后，求解不平衡量；

3.形成雅可比矩阵；

4.求解修正量后，重新修改初值，从2开始重新循环计算；

5.求解的电压变量达到所要求的精度时，再计算各支路功率分布、功率损耗和节点

6.上机编程调试；连调；

7.计算分析给定系统潮流分析并与手工计算结果作比较分析。

8.准备计算机演示答辩，书写该课程设计说明书(必须计算机打印)

四、13节点配网潮流计算

4.1牛拉法的原理及其基本方程

牛顿迭代法(Newton'smethod)又称为牛顿-拉夫逊(拉弗森)方法(Newton-Raphson

method),它是牛顿在17世纪提出的一种在实数域和复数域上近似求解方程的方法。

设r是f(x)=0的根，选取x0作为r初始近似值，过点(x0,f(xO))做曲线y=

f(x)的切线L,L的方程为y=f(xO)fr(xO)(x-xO),求出L与x轴交点的横坐标xl=

xO-f(xO)/ff(xO),称xl为r的一次近似值。过点(xl,f(xl))做曲线y=f(x)的切

线，并求该切线与x轴的横坐标x2=xl-f(xl)/f,(xl),称x2为r的二次近似值。重

复以上过程，得r的近似值序列，其中x(n+1)=x(n)-f(x(n))/f(x(n)),称为r的n+1

次近似值，上式称为牛顿迭代公式。解非线性方程f(x)=O的牛顿法是把非线性方程线

性化的一种近似方法。把f(x)在xO点附近展开成泰勒级数f(x)=f(x0)+(x-

4

华南理工大学广州学院电气工程学院课程设计报告

xO)f(x0)+(x-x0)^2*f,,(x0)/2!+…取其线性部分，作为非线性方程f(x)=0的近

似方程，即泰勒展开的前两项，则有f(xO)+f'(xO)(x-xO)=f(x)=O设方(xO)WO则其

解为xl=xO—f(xO)/f'(xO)这样，得到牛顿法的一个迭代序列：x(n+l)=x(n)-

f(x(n))/f'(x(n))。

4.2PQ分解法的原理及其基本方程

PQ分解法的基本思想是根据电力系统实际运行特点：通常网络上的电抗远大于电

阻值，则系统母线电压幅值的微小变化对母线有功功率的改变△尸影响很小。同样,

母线电压相角的少许改变A。,也不会引起母线有功功率的明显改变，因此，节

点功率方程在用极坐标形式表示时，它的修正方程式可简化为

一△OR。L^U/U\(2.1)

这就是把2(n-1)阶的线性方程组变成了两个n-L阶的线性方程组，将P和Q分

开来进行迭代计算，因而大大地减少了计算工作量。但是H、L在迭代过程中仍然在不

断的变化而且又都是不对称矩阵。对牛顿法的进一步简化，即把式(2-1)中的系数矩

阵简化为在迭代过程中不变的对称矩阵。

在一般情况下，线路两端的电压相角“是不大的(不超过10°〜20°),因此，

cos%«1

Gysin«与J

可以认为(2-2)

此外，与系统各节点无功功率相应的导纳远小于该节点自导纳的虚部，即

Q

BLDi=方"线

因而2<<U科(2-3)

考虑到以上关系，式(2-0的系数矩阵中的个元素可以表示为

居

(i,j=l,2,…，n-1)

5

华南理工大学广州学院电气工程学院课程设计报告

4=0"乌

(i,j=l,2,…，m)

而系数矩阵H和L则可以分别写成:

3"…。％_亿一|

U282MU2B22U2...3色，一〃,1

UnIu.U；,2U”——“T

n-\A/i—IJti口

%…4,n=lU、

B?2-^2,n=l

J■

♦♦•

B“T2…BnUz

=UDPUD\

q4M

u、BM...u^inum

=uju、U2B22U2...U2B2mUtn

••••••••••••

SAMUmBm2U2...U,〃B.JJ〃

U、42

3

B223

=Un2B-UD1

(2-5)

%=G尹丹

Ui=U/&=Uj(cosa+Jsin0中,得到

将式（2-4）和（2-5）代人式

叫］陷

[AQ]=-[UD2][^][AU]

华南理工大学广州学院电气工程学院课程设计报告

用和分别左乘以上两式，便得

[〃俏阳=-网网必句

(2-6)

(2-7)

这就是简化了的修正方程式,它们也可以展开写成

空

a%…"I■UQ4

Afi%

“…U2\02

——

u2••••••••••••

尻一\…纥.If

一%」

(2-8)

1

512

U\pH…B,1,/n一

必B?2…B>

一B”2,7wu?

U?

屹,纥1,2•••BUm

1nijn_

(2-9)

在这两个修正方程式中系数矩阵元素就是系统导纳矩阵的虚部，因而系数矩阵是对

称矩阵，且在迭代过程中保持不变。这就大大减少了计算工作量。

用极坐标表示的节点功率增量为

△片=匕-uEUj®cos%+8sin%)=0

△2=Qk-4XUcos0.-Bsin%)=0

J=i

(2-10)

式(2-8)、(2-9)和(2-10)构成了PQ分解法迭代过程的基本方程式。

7

华南理工大学广州学院电气工程学院课程设计报告

4.3配网前推后代潮流计算的原理

前推回代法在配电网潮流计算中简单实用，所有的数据都是以矢量形式存储，因此

节省了大量的计算机内存，对于任何种类的配电网只要有合理的R/X值，此方法均可

保证收敛。算法的稳定性也是评价配电网潮流算法的重要指标。一般情况下，算法的收

敛阶数越高，算法的稳定性越差，前推回代法的收敛阶数为一阶，因此它也具有较好的

稳定性。比较而言，前推回代法充分利用了网络呈辐射状的结构特点，数据处理简单，

计算效率高,具有较好的收敛性,被公认是求解辐射状配电网潮流问题的最佳算法之一。

4.4配网前推后代潮流计算迭代过程

配网潮流前推回代潮流算法第n+1步的迭代公式如下：节点i的前推计算公式为

由图1可知节点i+1的注入有功功率和无功功率分别为：

NBNB-1

P(i+1)=ZPL(j)十ZLP(j)

Q(i+l)=ZQL(j)+gLQ(j)

其中：i=L2,000NB-1,NB为节点数，PL⑴为第i节点所带负荷有功功率；QL(i)

为第i节点所带负荷无功功率；LP⑴为第i条线段上的有功功率损耗；LQ(i)为第i条线

段上的无功功率损耗。

收敛条件

以前后两次迭代的电压偏差作为迭代收敛条件，(8)式表明节点电压幅值最大偏差

小于设定阈值，即认为迭代收敛，则迭代过程结束。

max{|Vk(i)-Vk-l(i)|}<e.(8)

在开始计算的时候先把各节点上的LP⑴、LQ⑴分别设为零，利用式⑴计算出

P(i+1)、Q(i+1),分别代入式⑹计算各节点电压幅值，作为初始化的节点电压值.然后把

此电压值代入公式⑵来修正线路的损耗LP⑴、LQ(i),用经过修正的线路损耗求新的节

8

华南理工大学广州学院电气工程学院课程设计报告

点注入功率，再根据节点注入功率求新的节点电压，根据前后两次迭代的电压偏差是否

小于设定阈值判断是否收敛，如果满足，则输出结果，如果不满足，必须重复上述过程

直到满足收敛条件为止。

4.5计算原理

第一步，从与叶节点联系的支路开始，该支路的末端功率即等于叶节点功率，利用这

个功率和对应的节点电压计算支路功率损耗，求得支路的首端功率。当以某节点为始节

点的各支路都计算完毕后，便想象将这些支路都拆去，是的该节点成为新的叶节点，其

节点功率等于原有的负荷功率与该节点为始节点的各支路首端功率之和。于是计算便可

延续下去，直到全部支路计算完毕。

第二步，利用第一步所得的支路首端功率和本步骤刚算出的本支路始节点的电压

（对电源点为已知电压），从电源点开始逐条支路进行计算，求得各支路终节点的电压。

在迭代计算开始之前，要先处理好支路的计算顺序问题。介绍两种确定之路计算顺

序的方法。

第一种方法是，按与叶节点联接的之路顺序，并将已排序的之路拆除，在此过程中

将不断出现新的叶节点，而与其联接的支路又加入排序行列。这样就可以全部排列好从

叶节点向电源点计算功率损耗的支路顺序。其逆序就是进行电压计算的支路顺序。

第二种是逐条追加支路的方法。首先从根节点（电源点）开始接出第一条支路，引出一

个新节点，以后每次追加的支路都必须从已出现的节点接出遵循这个原则逐条追加支

路，直到全部支路追加完毕。所得到的支路追加顺序即是进行电压计算的支路顺序，其

逆序便是功率损耗计算的支路顺序。

9

华南理工大学广州学院电气工程学院课程设计报告

10

华南理工大学广州学院电气工程学院课程设计报告

五、计算框图流程

11

华南理工大学广州学院电气工程学院课程设计报告

六、前推回代计算输入文件

主程序：

[PQ,FT9RX]=casell5();%调用数据文件

%节点数

NN=size(PQ9l)；

NB=size(FT,l)；%支路数数

V=PQ(:,1)；%V初始电压相量

maxd=l

whilemaxd>0.0001

PQ2=PQ;%每一次迭代各节点的注入有

功和无功相同

PL=O.O;

fori=l:NB

kf=FT(i,l);%前推始节点号

kt=FT(i,2);%前推终节点号

x=(PQ2(kf,2)A2+PQ2(kf,3)A2)/(V(kf)A2);%计

算沿线电流平方A

PQl(i,l)=RX(i,l)*x+PQ2(kf,2);%计

算支路首端有功/MWRX(i,l)-R

PQl(i,2)=RX(i,2)*x+PQ2(kf,3);%计算沿支路

的无功损耗/MvarRX(i,2)~X

12

华南理工大学广州学院电气工程学院课程设计报告

PQ2(kt,2)=PQ2(kt,2)+PQl(i,l);%用PQ1去

修正支路末端节点的有功P单位MW

PQ2(kt,3)=PQ2(kt,3)+PQl(i,2);%用PQ1去

修正支路末端节点的有功Q单位Mvar

PL=PL+RX(i,l)*x;%累积计算

整个网络的功率损耗

end

angle(l)=O.O;

fori=NB:-l:l

kf=FT(i,2);%回代始节点号

kt=FT(i,l);%回代终节点号

dvl=(PQl(i,l)*RX(i,l)+PQl(i,2)*RX(i,2))/V(kf);

%计算支路电压损耗的纵分量dvl

dv2=(PQl(i1)*RX(i,2)・PQl(i,2)*RX(i1))/V(kf)；

%计算支路电压损耗的横分量dv2

V2(kt)=sqrt((V(kf)-dvl)A2+(dv2)A2);

%计算支路末端电压/kV

angle(kt)=

atand(dv2/(V(kf)-dvl))+angle(kf);%计算支

华南理工大学广州学院电气工程学院课程设计报告

路电压相角(单位度)

end

maxd=abs(V2(2)-V(2));

V2(1)=V(1);

fori=3:l:NN

ifabs(V2(i)-V(i))>maxd;

maxd=abs(V2(i)-V(i));

end

end

maxd

k=k+l

PQl%潮流分布即支路首端潮流

MVA

V=V2%节点电压模计算结果kV

angle%节点电压角度计算结果单位

度

PL%网损单位MW

end

clear

14

华南理工大学广州学院电气工程学院课程设计报告

输入文件清单：

function[PQ9FT,RX]=casell5()

PQ=[

%节点电压有功无功

10.400

10.00.03420.0301

10.00.06930.0642

10.00.08450.0763

10.00.02950.0261

10.00.04740.0409

10.00.11760.0957

10.00.09460.0857

10.00.09160.0859

10.00.02710.0229

10.00.06960.0643

10.00.06760.0579

10.00.02980.0242

];

FT=[

%首端末端

54

134

15

华南理工大学广州学院电气工程学院课程设计报告

1211

43

103

113

76

98

32

62

82

21

];

RX=[

%RX

4.5245.04

3.5213.966

2.4362.866

1.1451.28

4.144.696

1.3281.763

2.7452.965

2.2372.756

2.3562.541

16

华南理工大学广州学院电气工程学院课程设计报告

0.8561.14

3.7434.251

3.3673.685

];

计算结果：

k=8

PQ1=

0.02960.0262

0.02990.0243

0.06780.0581

0.14440.1273

0.02720.0230

0.13790.1231

0.11830.0964

0.09200.0864

0.38520.3445

0.16610.1379

0.18910.1750

0.81160.7278

V=10.40009.87959.69919.6652

17

华南理工大学广州学院电气工程学院课程设计报告

9.63779.84929.78729.7326

9.68709.67649.65789.62359.6444

angle00.30110.40110.4250

0.44380.34310.39420.38990.4266

0.42090.44980.48230.4452

PL=0.0488

数据分析：

经过8次迭代，各节点的电压误差均在O.OOOlkv以内，计算到此结束。

七、配电网潮流计算的要点

1.需要注意系统中所有母线的电压是否在允许的范围以内。

2.系统中各种元件（线路、变压器等）是否会出现过负荷。

3.在进行电压和功率分布计算以前，先要对网络等值电路作些简化处理。

4.确立前推回代的计算次序。

5.在开始编写程序之前，要了解matlab仿真软件，包括其程序的编写，程序的运

行，我们都要清楚的知道

6.为了提高计算精度，重复计算，给出一个容许误差作为计算的收敛的判据。

八、自我总结

此次课程设计首先让我明白了要使电力系统运行的稳定，必须经过精密的设计和计

算。在进行课题设计的过程中，加深了我对潮流计算的认识，尤其是对牛顿拉夫逊潮流

计算的求解思路有了比较透彻的理解。

在此次程序设计中，我在MATLAB编程，13节点配电网潮流计算，WORD文档的编辑

18

华南理工大学广州学院电气工程学院课程设计报告

方面均有所提高，之前学过的潮流计算仅限于书本，通过完成这个作业使我对这块知识

了解得更加全面，但也暴露出了一些问题：理论知识还是不足，基本知识不充足，做起

这个大作业难度还是很大的，对MATLAB的性能和特点还不能有一个全面的把握，相信

通过以后的学习能弥补这些不足，从而达到一个新的层次。

在过程设计中我发现了一些要点，要求得支路的首端功率，从与叶节点联接的支路

开始，该支路的末端功率即等于叶节点功率，利用这个功率和对应的节点电压计算支路

功率损耗，才能求首段功率。还有在迭代计算开始之前，先要确定支路的计算顺序问题,

即确定好是深度优先还是广度优先。要求得各支路终结点的电压的话利用上述所得的支

路首端功率和本步骤刚算出的本支路始节点的电压，从电源点开始逐条支路进行计算。