13节点配电网潮流计算—上机
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
配网前推后代潮流计算,
对给定的13节点配电网,熟悉MATLAB语言前腿后代潮流计算程序。
13节点配电网结构如图1所示,支路参数见表1;节点类型及迭代初始电压及各节点所接负荷见表2。
要看懂前推回代法计算程序,报告叙述计算原理及计算流程。绘制计算流程框图。确定前推回代支路次序(广度优先,或深度优先),编写前推回代计算输入文件。进行潮流计算。
整理专利计算结果。
内容包括节点电压;线路潮流分布,网损
总结配电网潮流计算的要点
图1 13节点配电网结构图
表1 系统支路参数
表2 系统负荷参数
主程序清单:
[PQ,FT,RX]=case113(); %调用数据文件
NN=size(PQ,1); %节点数
NB=size(FT,1); %支路数数
V=PQ(:,1); %V初始电压相量
maxd=1
k=1
while maxd>0.0001
PQ2=PQ; %每一次迭代各节点的注入有功和无功相同
PL=0.0;
for i=1:NB
kf=FT(i,1); %前推始节点号
kt=FT(i,2); %前推终节点号
x=(PQ2(kf,2)^2+PQ2(kf,3)^2)/V(kf)/V(kf); %计算沿线电流平方A
PQ1(i,1)=PQ2(kf,2)+RX(i,1)*x; %计算支路首端有功/MW RX(i,1)~R
PQ1(i,2)=PQ2(kf,3)+RX(i,2)*x; %计算沿支路的无功损耗/Mvar RX(i,2)~X PQ2(kt,2)= PQ2(kt,2)+PQ1(i,1); %用PQ1去修正支路末端节点的有功P 单位MW
PQ2(kt,3)= PQ2(kt,3)+PQ1(i,2); %用PQ1去修正支路末端节点的有功Q 单位Mvar
PL=PL+RX(i,1)*x;
end
angle(1)=0.0;
for i=NB:-1:1
kf=FT(i,2); %回代始节点号
kt=FT(i,1); %回代终节点号
dv1=(PQ1(i,1)*RX(i,1)+PQ1(i,2)*RX(i,2))/V(kf); %计算支路电压损耗的纵分量dv1 dv2=(PQ1(i,1)*RX(i,2)-PQ1(i,2)*RX(i,1))/V(kf); %计算支路电压损耗的横分量dv2 V2(kt)=sqrt((V(kf)-dv1)^2+dv2^2); %计算支路末端电压/kV
angle(kt)=angle(kf)+atand(dv2/(V(kf)-dv1)); %计算支路
end
maxd=abs(V2(2)-V(2));
V2(1)=V(1);
for i=3:1:NN
if abs(V2(i)-V(i))>maxd;
maxd=abs(V2(i)-V(i));
end
end
maxd
k=k+1
PQ1 %潮流分布即支路首端潮流MVA V=V2 %节点电压模计算结果kV
angle %节点电压角度计算结果单位度
PL %网损单位MW
end
clear
输入文件清单:
function [PQ,FT,RX]=case113()
PQ=[
%节点电压有功无功
10.5 0 0
10. 0.6 0.45
10 0.4 0.3
10 0.4 0.28
10 0.6 0.4
10 0.4 0.3
10 0.5 0.35
10 0.5 0.4
];
FT=[
%首端末端
4 3
6 5
7 5
3 2
5 2
8 2
2 1
];
RX=[
% R X
0.6 0.35
1.0 0.55
0.65 0.35
0.62 0.5
0.72 0.75
0.90 0.5
0.54 0.65
];
计算结果清单:
k =7
PQ1 =
0.4014 0.2808
0.4025 0.3014
0.5025 0.3513
0.8074 0.5857
1.5294 1.0782
0.5036 0.4020
3.5356 2.6303
V =10.5000 10.1557 10.0776 10.0439 9.9677 9.9107 9.9226 10.0913
angle = 0 0.4716 0.4944 0.4785 0.6815 0.6351 0.6511 0.4101
PL =0.1356MW