潮流计算报告教材

一、系统结构图：

二、网络参数：

网络

类型

支路

编号

电压

等级

（kV）

装机

容量

(MW)

导线的技术参数

1

r（Ω）

1

x（Ω）1

b

（）

环网1-2

220

100

13.6 125.5 67.85 1-3 8.321 130.5 52.24 3-5 10.2 128.8 74.99 2-3 8.5 105.4 28.36 1-4 7.579 129.6 51.45 4-5 13.84 125.31 2.78

辐射网1-2

————

2 4

——2-3 3 6

3-4 4 8

2-5 1 2

5-6 4 4

节点类型节点编

号

发电功

率（MW）

负荷视在功率

1 0 未知(平衡节点)

环网2 100 0（PV节点）

3 0 15+9.4i

4 0 27+6i

5 0 35.5+25.5i

辐射网1 0 未知(平衡节点)

2 0 4+2i

3 0 6+3.2i

4 0 3+1.44i

5 0 4+3.2i

6 0 2+1.1i

四、matlab程序：

clear;

Un=input('请输入Un:'); %输入所需的额定电压PQ=[

%节点电压有功无功

Un 0 0

Un 4 2

Un 6 3.2

Un 3 1.44

Un 4 3.2

Un 2 1.1

];

FT=[

%首端末端

4 3

3 2

6 5

5 2

2 1

];

RX=[

% R X

4 8

3 6

4 4

1 2

2 4

];

NN=size(PQ,1); %节点数

NB=size(FT,1); %支路数数

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

maxd=1

k=1

while maxd>0.0001

k=k+1;

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

losss(i,1)=RX(i,1)*x; %计算线路有功损耗 /MW

losss(i,2)=RX(i,2)*x; %计算线路无功损耗/MW

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; %计算沿支路的无功/MW 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

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

fullloss(1,1)=0;%计算线路总损耗

fullloss(1,2)=0;

finalPQ=max(PQ1);

for i=1:NB

fullloss(1,1)=fullloss(1,1)+losss(i,1);

fullloss(1,2)=fullloss(1,2)+losss(i,2);

end

disp('辐射网迭代次数：')

k

disp('辐射网系统电压差精度：')

maxd

disp('辐射网系统末端节点有功和无功：')

finalPQ %潮流分布即支路首端潮流MVA

disp('辐射网系统总功率损耗：')

fullloss %线路总损耗MVA

disp('辐射网系统各支路功率损耗：')

losss %各支路损耗MVA

disp('辐射网系统各节点电压幅值：')

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

disp('辐射网系统各节点电压相角：')

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

end

clc

disp('辐射网迭代次数：')

k

disp('辐射网系统电压差精度：')

maxd

disp('辐射网系统末端节点有功和无功/MVA：')

FinPQ=finalPQ(1,1)+finalPQ(1,2)*j %潮流分布即支路首端潮流MVA

disp('辐射网系统总功率损耗/MVA：')

Fulloss=fullloss (1,1)+fullloss(1,2)*j %线路总损耗MVA

disp('辐射网系统各支路功率损耗/MVA：')

for(a=1:5)

LOSS=losss (a,1)+losss(a,2)*j %各支路损耗MVA