牛拉法潮流计算程序附3机9节点结果对比.docx
- 文档编号:11210699
- 上传时间:2023-05-29
- 格式:DOCX
- 页数:13
- 大小:346.88KB
牛拉法潮流计算程序附3机9节点结果对比.docx
《牛拉法潮流计算程序附3机9节点结果对比.docx》由会员分享,可在线阅读,更多相关《牛拉法潮流计算程序附3机9节点结果对比.docx(13页珍藏版)》请在冰点文库上搜索。
牛拉法潮流计算程序附3机9节点结果对比
摘要
电力系统潮流计算是研究电力系统稳态运行的一种重要方法,它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态,包括各母线的电压、线路的功率分布以及功率损耗等等。
潮流计算主要用于电网规划和静态安全分析,它可为扩建电力网络,以达到规划周期内所需要的输电能力提供依据;也可以对预想事故进行模拟和分析,校核预想事故下的电力系统安全性。
本文简单介绍了牛顿-拉夫逊潮流计算的原理、模型与算法,然后用具体的实例,利用MATLAB对牛顿-拉夫逊法的算法进行了验证。
关键词:
电力系统潮流计算牛顿-拉夫逊法MATLAB
一、牛拉法的数学模型
对一个N节点的电力网路,列写节点电压方程,即
(1.1)
式中,
为节点注入电流列相量,Y为节点导纳矩阵,
为节点电压列相量。
由于异地测量的两个电流缺少时间同步信息,以注入功率替换注入电流作为已知量。
即
(1.2)
其中,
,带入上式,得到有功功率和无功功率方程
(1.3)
(1.4)
大部分情况下,已知PQ,求解Vθ。
考虑到电网的功率平衡,至少选择一台发电机来平衡全网有功功率,即至少有一个平衡节点,常选择调频或出线较多的发电机作为平衡节点。
具有无功补偿的母线能保持电压幅值恒定,这类节点可作为PV节点。
潮流计算中节点分类总结如下:
表1.1节点分类
PQ节点
负荷母线、变电站母线
绝大部分是PQ节点
已知PQ,求Vθ
PV节点
有无功储备的发电机母线
很小部分是PV节点
已知PV(V=1),求Qθ
平衡节点
选择一个发电机母线
已知Vθ(V=1,θ=0),求PQ
电压相量作为参考轴θ=0
已知电力系统有m个PQ节点,r个PV节点和1个平衡节点,则可以提取m+r个有功功率方程和m个无功功率方程,从而求解出m+r个θ和m个V,其余节点的有功和无功可通过式(1.3)、(1.4)求得,这样就完成了潮流计算。
二、潮流计算流程
用于潮流计算的有功和无功功率方程的阶数一般很高,常选择牛顿-拉夫逊迭代法进行求解高维的非线性方程组。
max
max
牛顿-拉夫逊法的迭代过程如下:
图2.1牛拉法迭代流程图
将功率方程改成
的形式
(2.1)
(2.2)
可以得到
(2.3)
为了便于计算,将上式改为
(2.4)
求取功率方程的雅克比矩阵
i≠j时,
(2.5)
(2.6)
(2.7)
(2.8)
i=j时,
(2.9)
(2.10)
(2.11)
(2.12)
牛拉法求解潮流的迭代过程如下
图2.2牛拉法潮流迭代流程图
三、MATLAB编程实例
实例采用王锡凡主编的《现代电力系统分析》p326-327的例题。
图3.1某小型电力网络图
例题中,1号发电机为平衡节点,2号和3号发电机为PV节点,其余为PQ节点。
为了便于编程,可将PV节点和平衡节点放在节点描述图的最后,故要重新对节点进行编号。
。
图3.2修改节点编号后的电力网络图
支路表如下:
%首节点末节点电阻电抗容纳之半变比
120.0100.0850.0881
130.0170.0920.0791
240.0320.1610.1531
360.0390.1700.1791
450.00850.0720.07451
560.01190.10080.10451
9100.057601
7400.062501
8600.058601
节点表如下:
%节点电压幅值电压相角节点有功节点无功节点类型/“1”为PQ节点,“2”为PV节点,“3”为平衡节点
11.00001
21.00-1.2500-0.51
31.00-0.9-0.31
41.00001
51.00-1-0.351
61.00001
71.02501.630002
81.02500.850002
91.0400003
程序清单:
取ε=10-5
a=load('zhilu.txt');
b=load('jiedian.txt');
N=size(a,1);%支路数
Nbus=size(b,1);%节点数
pq=0;
fork=1:
Nbus%PQ节点个数
ifb(k,6)==1
pq=pq+1;
end
end
Y=zeros(Nbus);%节点导纳矩阵
fork=1:
N
t1=a(k,1);t2=a(k,2);r=a(k,3);x=a(k,4);ban=a(k,5);K=a(k,6);
Y(t1,t1)=Y(t1,t1)+1/(r+j*x)+j*ban;
Y(t1,t2)=Y(t1,t2)-1/(K*(r+j*x));
Y(t2,t1)=Y(t2,t1)-1/(K*(r+j*x));
Y(t2,t2)=Y(t2,t2)+1/(K*K*(r+j*x))+j*ban;
end
G=real(Y);B=imag(Y);
precision=1;
t=0;%存储迭代次数
whileprecision>0.00001
P=zeros(Nbus,1);%存储所有节点的有功
Q=zeros(Nbus,1);%存储所有节点的无功
form=1:
Nbus%求pv、pq和平衡节点的有功
forn=1:
Nbus
P(m,1)=P(m,1)+b(m,2)*b(n,2)*(G(m,n)*cos(b(m,3)-b(n,3))+B(m,n)*sin(b(m,3)-b(n,3)));
end
end
form=1:
Nbus%求pq、pv和平衡节点的无功
forn=1:
Nbus
Q(m,1)=Q(m,1)+b(m,2)*b(n,2)*(G(m,n)*sin(b(m,3)-b(n,3))-B(m,n)*cos(b(m,3)-b(n,3)));
end
end
deltp=b(1:
Nbus-1,4)-P(1:
Nbus-1,1);%pq和pv节点有功差
deltq=b(1:
pq,5)-Q(1:
pq,1);%pq节点无功差
deltPQ=[deltp;deltq];
H=zeros(Nbus-1);
form=1:
Nbus-1%求H矩阵
forn=1:
Nbus-1
ifm~=n
H(m,n)=-b(m,2)*b(n,2)*(G(m,n)*sin(b(m,3)-b(n,3))-B(m,n)*cos(b(m,3)-b(n,3)));
else
H(m,m)=b(m,2)*b(m,2)*B(m,m)+Q(m,1);
end
end
end
N=zeros(Nbus-1,pq);
form=1:
Nbus-1%求N矩阵
forn=1:
pq
ifm~=n
N(m,n)=-b(m,2)*b(n,2)*(G(m,n)*cos(b(m,3)-b(n,3))+B(m,n)*sin(b(m,3)-b(n,3)));
else
N(m,m)=-b(m,2)*b(m,2)*G(m,m)-P(m,1);
end
end
end
J=zeros(pq,Nbus-1);
form=1:
pq%求J矩阵
forn=1:
Nbus-1
ifm~=n
J(m,n)=b(m,2)*b(n,2)*(G(m,n)*cos(b(m,3)-b(n,3))+B(m,n)*sin(b(m,3)-b(n,3)));
else
J(m,m)=b(m,2)*b(m,2)*G(m,m)-P(m,1);
end
end
end
L=zeros(pq,pq);
form=1:
pq%求L矩阵
forn=1:
pq
ifm~=n
L(m,n)=-b(m,2)*b(n,2)*(G(m,n)*sin(b(m,3)-b(n,3))-B(m,n)*cos(b(m,3)-b(n,3)));
else
L(m,m)=b(m,2)*b(m,2)*B(m,m)-Q(m,1);
end
end
end
Jacobi=[HN;JL];%雅克比矩阵
Correction=-Jacobi\deltPQ;%计算电压相角和幅值的修正量
form=1:
Nbus-1
b(m,3)=b(m,3)+Correction(m);
end
form=1:
pq
b(m,2)=b(m,2)+Correction(Nbus-1+m)*b(m,2);
end
precision=max(abs(deltPQ));
t=t+1
b
end
b(Nbus,4)=P(Nbus,1);
b(Nbus,5)=Q(Nbus,1);
form=pq+1:
Nbus-1
b(m,5)=Q(m,1);
end
b
运行结果:
四、结果分析与结论
运行结果与图3.2对应,例题给定的系统潮流与图3.1对应,对比之后发现结果完全一致,验证了牛顿-拉夫逊潮流计算的正确性。
程序迭代四次就完成,可见牛顿-拉夫逊算法的收敛性好,迭代次数少。
但每一次迭代都要更新雅克比矩阵,而且PQ方程联立求解,方程维数很高。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 牛拉法 潮流 计算 程序 节点 结果 对比