电力系统 潮流计算仿真报告

Beijing Jiaotong University

电力系统潮流计算

仿真报告

姓名：TYP

班级：电气0906

学号：09291183

指导老师：吴俊勇

完成日期：2012.6.24

一、实验内容

电力系统潮流计算是研究电力系统稳态运行情况的一种基本电气计算。它的任务是根据给定的运行条件和网路结构确定整个系统的运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等。电力系统潮流计算的结果是电力系统稳定计算和故障分析的基础。

对于简单系统，可以将其分为开式网络和闭式网络手工计算。对于复杂电力系统，根据定解条件，应用牛顿—拉夫逊法进行计算，在手工计算中，由于涉及大量变量、微分方程、矩阵计算，求解很烦琐，而且容易出错，计算不同系统时需要重新计算。故而我们可以借助计算机来进行潮流计算，方便快捷且准确率高。

二、计算机潮流计算方法

我们常用牛顿—拉夫逊法来进行潮流计算。

牛顿—拉夫逊法(简称牛顿法)在数学上是求解非线性代数方程

式的有效方法，其要点是把非线性方程式的求解过程变成反复地对相应的线性方程式进行求解的过程，即通常所称的逐次线性化过程。

1、基本原理

从几何意义上，牛顿—拉夫逊法实质上就是切线法，是一种逐步线性化的方法。

2、牛顿—拉夫逊法潮流求解过程

以下讨论的是用直角坐标形式的牛顿—拉夫逊法潮流的求解过程。当采用直角坐标时，潮流问题的待求量为各节点电压的实部和虚部两个分量，由于平衡节点的电压向量是给定的，因此待求量共2(n-1)需要2(n-1)个方程式。事实上，除了平衡节点的功率方程式在迭代过程中没有约束作用以外，其余每个节点都可以列出两个方程式。

求解过程大致可以分为以下步骤：

(1)形成节点导纳矩阵；

(2)将各节点电压设初值；

(3)将节点初值代入相关求式，求出修正方程式的常数项向量；

(4)将节点电压初值代入求式，求出雅可比矩阵元素；

(5)求解修正方程，求修正向量；

(6)求取节点电压的新值；

(7)检查是否收敛，如不收敛，则以各节点电压的新值作为初值自第3步重新开始进行狭义次迭代，否则转入下一步；

(8)计算支路功率分布，PV节点无功功率和平衡节点功率。

其用程序仿真的过程可以由以下流程图来简单表示出来：

三、Matlab程序

1、原代码

%潮流计算

fprintf('开始潮流计算\n');

fprintf('请输入待求网络的相应参数\n');

%参数输入部分

n=input('网络中的节点数:n=');

L=input('网络中的支路数:L=');

ss=input('平衡节点ss=');

pr=input('误差精度:pr=');

X1=input('支路参数:X1=');

X2=input('节点参数:X2=');

X=input('节点号和对地参数:X=');

fprintf('参数输入部分结束\n\n');

Y=zeros(n);

%置迭代次数

mm=1;

%创建节点导纳矩阵

for i=1:L

if X1(i,6)==0 %不含变压器的支路

p=X1(i,1);

q=X1(i,2);

Y(p,q)=Y(p,q)-1/X1(i,3);

Y(q,p)=Y(p,q);

Y(p,p)=Y(p,p)+1/X1(i,3)+0.5*X1(i,4); Y(q,q)=Y(q,q)+1/X1(i,3)+0.5*X1(i,4);

else%含有变压器的支路

p=X1(i,1);

q=X1(i,2);

Y(p,q)=Y(p,q)-1/(X1(i,3)*X1(i,5)); Y(q,p)=Y(p,q);

Y(p,p)=Y(p,p)+1/X1(i,3);

Y(q,q)=Y(q,q)+1/(X1(i,5)^2*X1(i,3));

end

end

Y;

OrgS=zeros(2*n-2,1);

DetaS=zeros(2*n-2,1); %将OrgS、DetaS初始化%创建OrgS，用于存储初始功率参数

j=0;

for i=1:n %对PQ节点的处理

if i~=ss&&X2(i,6)==2

h=h+1;

for j=1:n

OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(X2(i,3))*(real(Y(i,j))*real(X2(j,3)) -imag(Y(i,j))*imag(X2(j,3)))+imag(X2(i,3))*(real(Y(i,j))*imag(X2(j,3) )+imag(Y(i,j))*real(X2(j,3)));

OrgS(2*h,1)=OrgS(2*h,1)+imag(X2(i,3))*(real(Y(i,j))*real(X2(j,3))-ima g(Y(i,j))*imag(X2(j,3)))-real(X2(i,3))*(real(Y(i,j))*imag(X2(j,3))+im ag(Y(i,j))*real(X2(j,3)));

end

end

end

for i=1:n %对PV节点的处理，注意这时不可再将h初始化为0 if i~=ss&&X2(i,6)==0

h=h+1;

for j=1:n

OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(X2(i,3))*(real(Y(i,j))*real(X2(j,3)) -imag(Y(i,j))*imag(X2(j,3)))+imag(X2(i,3))*(real(Y(i,j))*imag(X2(j,3) )+imag(Y(i,j))*real(X2(j,3)));

OrgS(2*h,1)=OrgS(2*h,1)+imag(X2(i,3))*(real(Y(i,j))*real(X2(j,3))-ima g(Y(i,j))*imag(X2(j,3)))-real(X2(i,3))*(real(Y(i,j))*imag(X2(j,3))+im ag(Y(i,j))*real(X2(j,3)));

end

end

end

OrgS;

%创建PVU 用于存储PV节点的初始电压

PVU=zeros(n-h-1,1);

t=0;

for i=1:n

if X2(i,6)==0

t=t+1;

PVU(t,1)=X2(i,3);

end

end

PVU;

%创建DetaS，用于存储有功功率、无功功率和电压幅值的不平衡量