前推回代法计算流程.docx

前推回代法计算流程.docx

《前推回代法计算流程.docx》由会员分享，可在线阅读，更多相关《前推回代法计算流程.docx（15页珍藏版）》请在冰豆网上搜索。

前推回代法计算流程

前推回代法计算流程

前推回代法计算流程

要看懂前推回代法计算程序，报告叙述计算原理及计算流程。

绘制计算流程框图。

确定前推回代支路次序（广度优先，或深度优先），编写前推回代计算输入文件。

进行潮流计算。

下列为节点配电网结构图及系统支路参数和系统负荷参数表。

图1-2节点配电网结构图

表1系统支路参数

支路

R（Ω）

X（Ω.）

B/2（TK）（S）

1~2

3.367

3.685

0.0

2~3

2.356

2.541

0.0

3~4

1.145

1.28

0.0

4~5

4.524

5.04

0.0

2~6

0.856

1.14

0.0

6~7

2.745

2.965

0.0

2~8

3.743

4.251

0.0

8~9

2.237

2.756

0.0

3~10

4.14

4.696

0.0

3~11

1.328

1.763

0.0

11~12

2.436

2.866

0.0

4~13

3.521

3.966

0.0

表2系统负荷参数

节点编号

节点类型

节点初始电压（kV）

Pi（MVA）

Qi（MVA）

1

根节点

10.4

0

0

2

中间节点

10.0

0.0342

0.0301

3

中间节点

10.0

0.0693

0.0642

4

中间节点

10.0

0.0845

0.0763

5

叶节点

10.0

0.0295

0.0261

6

中间节点

10.0

0.0474

0.0409

7

叶节点

10.0

0.1176

0.0957

8

中间节点

10.0

0.0946

0.0857

9

叶节点

10.0

0.0916

0.0859

10

叶节点

10.0

0.0271

0.0229

11

中间节点

10.0

0.0696

0.0643

12

叶节点

10.0

0.0676

0.0579

13

叶节点

10.0

0.0298

0.0242

主程序清单：

[PQ,FT,RX]=case114（）;%调用数据文件

NN=size（PQ,1）;%节点数

NB=size（FT,1）;%支路数数

V=PQ（:

1）;%V初始电压相量

maxd=1

k=1

whilemaxd>0.0001

PQ2=PQ;%每一次迭代各节点的注入有功和无功相同

PL=0.0;

fori=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;%计算支路首端有功/MWRX（i,1）~R

PQ1（i,2）=PQ2（kf,3）+RX（i,2）*x;%计算沿支路的无功损耗/MvarRX（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;

fori=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）;

fori=3:

1:

NN

ifabs（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]=case114（）

PQ=[

%节点电压有功无功

10.400

10.00.03420.0301

10.00.06930.0642

10.00.08450.0763

10.00.02950.0261

10.00.04740.0409

10.00.11760.0957

10.00.09460.0857

10.00.09160.0859

10.00.02710.0229

10.00.06960.0643

10.00.06760.0579

10.00.02980.0242

];

FT=[

%首端末端

54

134

43

103

1211

113

76

62

98

82

32

21

];

RX=[

%RX

4.5245.04

3.5213.966

1.1451.28

4.144.696

2.4362.866

1.3281.763

2.7452.965

0.8561.14

2.2372.756

3.7434.251

2.3562.541

3.3673.685

];

计算过程

maxd=

1

k=

1

maxd=

0.1780

k=

2

PQ1=

0.02960.0262

0.02990.0243

0.14430.1272

0.02720.0230

0.06780.0581

0.13780.1230

0.11820.0964

0.16600.1378

0.09200.0863

0.18900.1748

0.38470.3439

0.80990.7260

V=

Columns1through8

10.40009.88079.82209.96729.97349.97019.93909.8550

Columns9through13

9.95569.97809.96009.96689.9799

angle=

Columns1through8

00.30110.39860.42110.43870.34210.39160.3878

Columns9through13

0.42250.41730.44440.47470.4400

PL=

0.0471

maxd=

0.1787

k=

3

PQ1=

0.02960.0262

0.02990.0243

0.14430.1272

0.02720.0230

0.06780.0581

0.13780.1230

0.11820.0964

0.16600.1378

0.09200.0863

0.18900.1749

0.38490.3442

0.81120.7274

V=

Columns1through8

10.40009.87989.70049.78869.94059.85049.90899.7338

Columns9through13

9.81009.79969.78139.92679.9470

angle=

Columns1through8

00.30110.40110.42440.44210.34310.39290.3899

Columns9through13

0.42570.42040.44860.47910.4434

PL=

0.0484

maxd=

0.1793

k=

4

PQ1=

0.02960.0262

0.02990.0243

0.14440.1272

0.02720.0230

0.06780.0581

0.13790.1231

0.11820.0964

0.16610.1378

0.09200.0864

0.18910.1749

0.38510.3444

0.81150.7277

V=

Columns1through8

10.40009.87969.69949.66669.76149.84959.78849.7329

Columns9through13

9.68839.67779.65919.74749.7680

angle=

Columns1through8

00.30110.40110.42500.44330.34310.39420.3899

Columns9through13

0.42660.42090.44980.48140.4447

PL=

0.0487

maxd=

0.1226

k=

5

PQ1=

0.02960.0262

0.02990.0243

0.14440.1272

0.02720.0230

0.06780.0581

0.13790.1231

0.11830.0964

0.16610.1379

0.09200.0864

0.18910.1750

0.38520.3444

0.81150.7278

V=

Columns1through8

10.40009.87959.69919.66569.63919.84939.78759.7326

Columns9through13

9.68739.67679.65819.62489.6457

angle=

Columns1through8

00.30110.40110.42500.44380.34310.39420.3899

Columns9through13

0.42660.42090.44980.48220.4452

PL=

0.0487

maxd=

0.0010

k=

6

PQ1=

0.02960.0262

0.02990.0243

0.14440.1273

0.02720.0230

0.06780.0581

0.13790.1231

0.11830.0964

0.16610.1379

0.09200.0864

0.18910.1750

0.38520.3445

0.81160.7278

V=

Columns1through8

10.40009.87959.69919.66539.63819.84929.78739.7326

Columns9through13

9.68709.67649.65799.62389.6447

angle=

Columns1through8

00.30110.40110.42500.44380.34310.39420.3899

Columns9through13

0.42660.42090.44980.48230.4452

PL=

0.0488

maxd=

2.6021e-004

k=

7

PQ1=

0.02960.0262

0.02990.0243

0.14440.1273

0.02720.0230

0.06780.0581

0.13790.1231

0.11830.0964

0.16610.1379

0.09200.0864

0.18910.1750

0.38520.3445

0.81160.7278

V=

Columns1through8

10.40009.87959.69919.66529.63789.84929.78729.7326

Columns9through13

9.68709.67649.65789.62359.6445

angle=

Columns1through8

00.30110.40110.42500.44380.34310.39420.3899

Columns9through13

0.42660.42090.44980.48230.4452

PL=

0.0488

maxd=

6.1046e-005

k=

8

PQ1=

0.02960.0262

0.02990.0243

0.14440.1273

0.02720.0230

0.06780.0581

0.13790.1231

0.11830.0964

0.16610.1379

0.09200.0864

0.18910.1750

0.38520.3445

0.81160.7278

V=

Columns1through8

10.40009.87959.69919.66529.63779.84929.78729.7326

Columns9through13

9.68709.67649.65789.62359.6444

angle=

Columns1through8

00.30110.40110.42500.44380.34310.39420.3899

Columns9through13

0.42660.42090.44980.48230.4452

PL=

0.0488

计算结果清单：

maxd=1

k=1

maxd=0.1780

k=2

maxd=0.1787

k=3

maxd=0.1793

k=4

maxd=0.1226

k=5

maxd=0.0010

k=6

maxd=2.6021e-04

k=7

maxd=6.1046e-05

k=8

PQ1=

0.02960.0262

0.02990.0243

0.14440.1273

0.02720.0230

0.06780.0581

0.13790.1231

0.11830.0964

0.16610.1379

0.09200.0864

0.18910.1750

0.38520.3445

0.81160.7278

V=

10.40009.87959.69919.66529.63779.84929.78729.73269.68709.67649.65789.62359.6444

angle=

00.30110.40110.42500.44380.34310.39420.38990.42660.42090.44980.48230.4452

PL=0.0488

参考文献

[1]何仰赞温增银《电力系统分析》．华中科技大学出版社.

[2]李维波.《MATLAB在电气工程中应用》．中国电力出版社.2007