电力系统潮流计算课程设计论文Word格式.docx

电力系统潮流计算课程设计论文Word格式.docx

《电力系统潮流计算课程设计论文Word格式.docx》由会员分享，可在线阅读，更多相关《电力系统潮流计算课程设计论文Word格式.docx（30页珍藏版）》请在冰豆网上搜索。

建立数学模型，确定解算方法，制订计算流程，编制计算程序。

关键词

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

高斯消去法潮流计算Mtlab

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

在电力系统的正常运行中，随着用电负荷的变化和系统运行方式的改变，网络中的损耗也将发生变化。

要严格保证所有的用户在任何时刻都有额定的电压是不可能的，因此系统运行中个节点出现电压的偏移是不可避免的。

为了保证电力系统的稳定运行，要进行潮流调节。

随着电力系统及在线应用的发展，计算机网络已经形成，为电力系统的潮流计算提供了物质基础。

电力系统潮流计算是电力系统分析计算中最基本的内容，也是电力系统运行及设计中必不可少的工具。

根据系统给定的运行条件、网络接线及元件参数，通过潮流计算可以确定各母线电压的幅值及相角、各元件中流过的功率、整个系统的功率损耗等。

潮流计算是实现电力系统安全经济发供电的必要手段和重要工作环节，因此潮流计算在电力系统的规划设计、生产运行、调度管理及科学研究中都有着广泛的应用。

它的发展主要围绕这样几个方面：

计算方法的收敛性、可靠性；

计算速度的快速性；

对计算机存储容量的要求以及计算的方便、灵活等。

常规的电力系统潮流计算中一般具有三种类型的节点：

PQ、PV及平衡节点。

一个节点有四个变量，即注入有功功率、注入无功功率，电压大小及相角。

常规的潮流计算一般给定其中的二个变量：

PQ节点（注入有功功率及无功功率），PV节点（注入有功功率及电压的大小），平衡节点（电压的大小及相角）。

1、变量的分类：

负荷消耗的有功、无功功率——、、、

电源发出的有功、无功功率——、、、

母线或节点的电压大小和相位——、、、

在这十二个变量中，负荷消耗的有功和无功功率无法控制，因它们取决于用户，它们就称为不可控变量或是扰动变量。

电源发出的有功无功功率是可以控制的自变量，因此它们就称为控制变量。

母线或节点电压的大小和相位角——是受控制变量控制的因变量。

其中,、主要受、的控制,、主要受、的控制。

这四个变量就是简单系统的状态变量。

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

对控制变量

对没有电源的节点则为

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

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

2、节点的分类：

⑴第一类称PQ节点。

等值负荷功率、和等值电源功率、是给定的，从而注入功率、是给定的，待求的则是节点电压的大小和相位角。

属于这类节点的有按给定有功、无功率发电的发电厂母线和没有其他电源的变电所母线。

⑵第二类称PV节点。

等值负荷和等值电源的有功功率、是给定的，从而注入有功功率是给定的。

等值负荷的无功功率和节点电压的大小也是给定的。

待求的则是等值电源的无功功率，从而注入无功功率和节点电压的相位角。

有一定无功功率储备的发电厂和有一定无功功率电源的变电所母线都可以作为PV节点；

⑶第三类平衡节点。

潮流计算时一般只设一个平衡节点。

等值负荷功率、是给定的，节点电压的大小和相位也是给定的。

担负调整系统频率任务的发电厂母线往往被选作为平衡节点。

二．牛顿—拉夫逊法概要

1．首先对一般的牛顿—拉夫逊法作一简单的说明。

已知一个变量X函数为：

到此方程时，由适当的近似值出发，根据：

反复进行计算，当满足适当的收敛条件就是上面方程的根。

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

这一方法还可以做下面的解释，设第次迭代得到的解语真值之差，即的误差为时，则：

把在附近对用泰勒级数展开

上式省略去以后部分

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

比较两式，可以看出牛顿—拉夫逊法的休整量和的误差的一次项相等。

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

对其近似解得修正量可以通过解下边的方程来确定：

式中等号右边的矩阵都是对于的值。

这一矩阵称为雅可比（JACOBI）矩阵。

按上述得到的修正向量后，得到如下关系

这比更接近真实值。

这一步在收敛到希望的值以前重复进行，一般要反复计算满足

为预先规定的小正数，是第n次迭代的近似值

2．用牛顿法计算潮流时，有以下的步骤：

⑴输入线路，电气元件参数，形成节点导纳矩阵。

⑵给这各节点电压初始值。

⑶将以上电压初始值代入式（4—38a）~式（4—38c）或式（4—45c）、（4—45a），求出修正方程式中的不平衡量。

⑷将各节点电压的初值代入式（4-41a）、式（4-41b）或式（4-49a）~式（4-49d）,求修正方程式的系数矩阵——雅克比矩阵的各个元素。

⑸解修正方程式，求各节点电压的变化量，即修正量。

⑹计算各节点电压的新值，即修正后值

⑺运用各节点电压的新值自第三步开始进入下一次迭代。

⑻计算平衡节点功率和线路功率。

其中，平衡节点功率为

线路功率为

从而，线路上损耗的功率为

3．程序框图如下：

根据牛顿法，用MATLAB软件进行编程如下：

——————————————程——序————————————————

%本程序的功能是用牛顿——拉夫逊法进行潮流计算

functiony

clear;

clc;

n=input（'

请输入节点数：

n='

）;

nl=input（'

请输入支路数:

nl='

）

isb=input（'

请输入平衡母线节点号isb='

pr=input（'

请输入误差精度pr='

b=input（'

请输入由之路参数形成的矩阵B1='

V=input（'

请输入各节点参数形成的矩阵B2='

S=[0+0i;

-3.7-1.3i;

-2-1i;

-1.6-0.8i;

5];

%各节点的注入功率

%S=[0+0i;

-7-5i;

w1=zeros（2*n-2,1）;

P=real（S）;

Q=imag（S）;

e=zeros（1,n）;

f=zeros（1,n）;

E=zeros（1,n）;

Y=zeros（n）;

fori=1:

nl%导纳矩阵生成

p=b（i,1）;

q=b（i,2）;

Y（p,q）=Y（p,q）-1./（b（i,3）*b（i,5））;

Y（q,p）=Y（p,q）;

Y（q,q）=Y（q,q）+1./b（i,3）+b（i,4）./2;

Y（p,p）=Y（p,p）+1./（b（i,3）*b（i,5）^2）+b（i,4）./2;

end

disp（'

导纳矩阵Y：

'

disp（Y）

U=zeros（1,n）;

G=real（Y）;

B=imag（Y）;

n

e（i）=real（V（i,1））;

f（i）=imag（V（i,1））;

U（i）=V（i,3）;

B（i,i）=B（i,i）+V（i,3）;

T=0;

co=0;

d=0;

whileT==0

A=0;

co=co+1;

fori=2:

n%生成雅可比矩阵和功率修正量

forj=2:

x=0;

x1=0;

ifV（i,2）==2

forr=1:

x=x+（e（i）*（G（i,r）*e（r）-B（i,r）*f（r））+f（i）*（G（i,r）*f（r）+B（i,r）*e（r）））;

x1=x1+（f（i）*（G（i,r）*e（r）-B（i,r）*f（r））-e（i）*（G（i,r）*f（r）+B（i,r）*e（r）））;

end

w（2*i-1）=P（i）-x;

w（2*i）=Q（i）-x1;

elseifV（i,2）==3

w（2*i）=U（i）^2-（e（i）^2+f（i）^2）;

end

h=0;

h1=0;

ifi==j

ifr==i

continue

h=h+（G（i,r）*f（r）+B（i,r）*e（r））;

h1=h1+（G（i,r）*e（r）-B（i,r）*f（r））;

J（2*i-1,2*j-1）=2*G（i,i）*f（i）+h;

J（2*i-1,2*j）=2*G（i,i）*e（i）+h1;

J（2*i,2*j-1）=-2*B（i,i）*f（i）+h1;

J（2*i,2*j）=-2*B（i,i）*e（i）-h;

else

J（2*i-1,2*j-1）=-B（i,j）*e（i）+G（i,j）*f（i）;

J（2*i-1,2*j）=G（i,j）*e（i）+B（i,j）*f（i）;

J（2*i,2*j-1）=-G（i,j）*e（i）-B（i,j）*f（i）;

J（2*i,2*j）=-B（i,j）*e（i）+G（i,j）*f（i）;

J（2*i,2*j-1）=2*f（i）;

J（2*i,2*j）=2*e（i）;

J（2*i,2*j-1）=0;

J（2*i,2*j）=0;

%disp（J）

%disp（w）

fori=3:

2*n%高斯消去法求电压修正量

forj=3:

2*n

J1（i-2,j-2）=J（i,j）;

w1（i-2）=w（i）;

u=zeros（2*n-2,1）;

N=2*n-2;

fork=1:

N

m=0;

fori=k+1:

N

m=J1（i,k）./J1（k,k）;

w1（i）=w1（i）-m*w1（k）;

forj=k+1:

J1（i,j）=J1（i,j）-m*J1（k,j）;

u（N）=w1（N）./J1（N,N）;

fori=N-1:

-1:

1

c=0;

fork=i+1:

c=c+J1（i,k）*u（k）;

u（i）=（w1（