前推回代法计算流程.docx

前推回代法计算流程.docx

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

前推回代法计算流程

　　　　　　　前推回代法计算流程

要看懂前推回代法计算程序，报告叙述计算原理及计算流程．绘制计算流程框图.确定前推回代支路次序（广度优先，或深度优先）,编写前推回代计算输入文件。

进行潮流计算。

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

图１—2节点配电网结构图

表1　系统支路参数

支路

R（Ω）

Ｘ（Ω。

）

B/2（TK）（Ｓ）

１~２

3．3６7

３。

685

0.0

2~3

2．３56

2．５41

０．0

３～4

1。

１45

1。

２8

0。

0

4～５

4.５24

５.04

0.0

２~6

０.856

1。

14

0．0

６~７

２。

7４5

2。

965

0。

0

２~8

3.743

4。

25１

0.0

８~9

２．237

2.756

０．０

３~10

4.14

4.6９6

０。

０

3~１1

1.32８

1.763

0．0

11～1２

２.43６

２.866

０.0

4~1３

3.5２1

3。

9６6

0．0

表2系统负荷参数

节点编号

节点类型

节点初始电压（ｋV）

Pi（ＭVA）

Qｉ（MVA）

１

根节点

10。

4

0

０

２

中间节点

10。

0

0。

０342

0.0301

3

中间节点

10.０

0。

069３

0.0６42

4

中间节点

10。

0

0.0８４５

０。

０７6３

５

叶节点

１0.0

0.０295

0。

０261

6

中间节点

10．0

0.0474

０.0４0９

７

叶节点

10。

0

０。

1176

0.０957

8

中间节点

10。

0

0。

0946

0.０８５7

９

叶节点

10.0

0.0916

0。

０85９

1０

叶节点

1０.０

０.0271

0。

022９

11

中间节点

１0。

０

0．0６96

0.０6４3

１２

叶节点

１0。

0

0．0６76

0。

0５79

１3

叶节点

10。

0

0。

0298

0。

024２

主程序清单：

[PQ，FT,RX]=caｓe1１4（）；　%调用数据文件

NN＝size（PQ，1）;　　%节点数

NB=sｉze（FT,１）;　%支路数数

Ｖ=ＰQ（：

1）；　　　%V初始电压相量

ｍaxd=1

ｋ=1

ｗhiｌemaｘd>0．0001

　PＱ2＝ＰQ；　　％每一次迭代各节点的注入有功和无功相同

　PＬ=0。

0；

forｉ=1：

NB

　kｆ=FＴ（i，1）;　%前推始节点号

　　kt=FT（i,2）;　　％前推终节点号

　x=（ＰQ２（kｆ,2）^２＋PQ２（ｋf,3）^2）／V（kf）／V（kf）；　　％计算沿线电流平方Ａ

　PQ1（i，１）=ＰQ2（kf,２）+RＸ（ｉ,1）＊x；　　％计算支路首端有功／MW　RＸ（i,1）~R

　PＱ1（i,2）＝PQ2（kf，3）+RX（i,２）*x;　％计算沿支路的无功损耗/MｖａrRX（i，2）~Ｘ

　PQ2（kｔ，2）=　PＱ2（kｔ,2）+PQ1（i，１）；　　%用PQ1去修正支路末端节点的有功P单位ＭW

　ＰＱ2（kt,3）=ＰＱ2（kt,3）＋PQ1（i,2）;％用PＱ1去修正支路末端节点的有功Ｑ单位Ｍvａr

　ＰL=ＰL+RX（i,1）＊x;

end

aｎgle

（1）=0.0;

foｒi=ＮB:

－1:

1

　　ｋｆ=FＴ（i,２）；　%回代始节点号　

　　kt=FT（ｉ，1）；％回代终节点号

　dｖ１=（PQ1（i，１）*RX（i，1）+PQ１（i,2）＊RX（i,２））/V（kｆ）;%计算支路电压损耗的纵分量dv1

　ｄv２＝（PＱ1（i,１）＊RＸ（i,２）-ＰＱ1（i,2）＊RX（i,１））/V（ｋf）；　％计算支路电压损耗的横分量ｄv2

V2（ｋt）=sｑrt（（V（kf）-dｖ1）＾2+ｄv２^2）；　　　　%计算支路末端电压/kV

angle（kt）=angle（kf）+aｔaｎｄ（dｖ2／（V（kｆ）-ｄv1））;　　％计算支路

end

maxd＝abs（V2

（2）-V

（2））;

V2

（1）=V（１）；

fｏr　i＝3:

1：

ＮＮ

　ifabs（V2（ｉ）—V（i））〉maｘd；

　maxd＝abs（V2（ｉ）—V（ｉ））;

eｎｄ　　

ｅnd

ｍaxd

k=k+1

PQ1　　　％潮流分布即支路首端潮流MVＡ

V=Ｖ2　　　％节点电压模计算结果ｋV

ａｎgｌe　　　%节点电压角度计算结果单位度

PＬ　　　　%网损单位ＭＷ

ｅｎd

cｌｅar

输入文件清单：

functiｏn［PQ，FT,RX]=caｓe114（）

PQ=[

%节点电压　有功　无功

10．4　０　　　0

10.０　　０。

03４20。

0301

1０.0　0.0６930。

0642

10．０　　0.０８４5０。

07６3

10.00.02９50。

0261

10。

０0。

04７40.0409

1０.0０.1１7６0.0957

10.0　　０。

09４6０.0857

１0．0　　0。

０916　0.08５9

10.0　　0。

02710.02２９

10.0　０。

0６96０.0６43

10.0　0。

06760.０5７9

10。

0　　　０.02９80。

02４２

］；

FＴ＝[

%首端末端　

5　　4

13　4

4　　　３

10　３

12　　11

11　　3

７　　6

6　　2

９　　8

8　　2

3　　　2

２　　１

］;

RX=[

％Ｒ　X　

4.52４　5。

04

３.521　　3。

966

１。

1451。

28

4。

1４　　4.69６

2.４362。

866

1。

3２8　1。

７63

2.7４5　2．９６5

0.8５６1。

14

２。

23７2。

7５6

3.７43４．251

2．3562。

541

3。

３６７３．685

]；

计算过程

maxd=

　1

ｋ　=

1

mａｘd=

　０.178０

k=

2

PQ1=

0．0296　　０。

0262

　　０．0299　0。

0243

　0。

１4430.1272

　0.０２72　0。

02３０

0.０67８　０。

0５81

０.１378　0.12３０

　0．1182　0.0９64

0.１６6０0。

1３７8

０.09２0　0.0863

　０.189０　　0。

174８

０.3８47　0．3439

0。

809９　0。

7260

V=

　Columns１tｈｒoｕｇh8

1０。

4０00　　9。

８8０７　　9．82209。

96７2　9。

97３4９。

9701　9．93９0９.85５0

　Cｏlumns９ｔhrｏugh1３

　9。

９556　　　9．97８０　　9。

9６０0９．9668　９。

9799

aｎgｌe=

Columｎs1　through　8

　　０　0.301１　0.3９8６　　0。

421１　　0.43８70．3421　　0.3916　0。

3８７8

Ｃolumｎs9throuｇｈ13

　0.4２25　0。

417３　　0。

４４44　　0．47４7　0。

4４00

PL＝

　0。

０４71

maxｄ=

　　0.1787

ｋ＝

　3

PQ1=

　　0.0２９6　０。

02６2

0。

02990.０24３

0。

144３　０。

1272

　０.02720。

023０

　　0。

0678　　0.０581

0。

1３780。

1230

　０。

118２　　0．0９６4

　0。

1６60　　　　0.1３７８

０。

0９20　　0.０863

　　0.1８９00。

1７4９

０.３8４9　０。

３４４2

　０．８112　0.72７４

Ｖ　=

Coluｍnｓ１tｈrough8

　10.40００9.8７９8　　　9.7０049。

78869。

9４059．8504　9。

9０899。

7338

　Coｌumns　9ｔhｒoｕgh13

　9．８100　　　9.79９6　9.781３　９．９26７　９.9470

anglｅ=

　Ｃolumnｓ1　througｈ８

　0　0.3011　0.4０11０。

424４　　0。

4４210.3431０。

39２9　　0。

3899

　Cｏｌumｎs9tｈrougｈ13

　０.４257　0。

4２０40.448６　０．479１　0。

４434

PL=

　0.04８4

maｘd＝

　　　0。

１７93

k　＝

４

PQ１　=

　0.0２９60。

02６２

　　0.0299　0.024３

　　0。

1444　0。

1272

0.０２7２0。

０２30

　０。

０６78　　　０。

０58１

0。

137９　０。

123１

　０。

１182　0．0９64

0．1６610．1378

　0。

０920　　　0.０8６4

　　0。

1891　　　0.1749

0。

３8510．3444

　　0。

81１5　　0.7277

V　＝

　Coｌｕmｎs1through8

10.４０00　9。

8796　9．6９94　9。

6666　9。

761４　９.8495９.78８49.７３29

Columns9throuｇh1３

　　9．6883　　9。

6７7７　9。

659１9．747４9.7680

angle=

　Coｌｕｍns１througｈ8

　　　00.3011　0.401１0。

42５0　0。

4433　0.３４3１０.3９４2０.3８99

　Columns9thrｏugh13

　　0。

4２66　0.４209０。

44９８　　　0.4814　　0．４447

ＰＬ=

　0.0487

maｘd　=

　　０.12２6

k＝

　5

PQ1　=

　0.02９6　0.０26２

　　0.0299　　０。

0243

　0.1444　0.１2７2

　０.027２　　0.０230

　0。

0６7８　　　0。

0581

　　　0。

１3790.1231

　0．１183　0.０964

　０．166１　　0.1３79

0.092０　0。

0８6４

　　0。

1891　　０。

17５0

　0.385２　　0.3444

　０.81150.７27８

Ｖ　=

　Cｏｌumns1　throuｇh　8

　　１0.４0００　9.8７9５　　９。

６991　9.6656　　9。

6391　9。

８４93　　　9.787５9.732６

　Columns9through　１３

9.6873　　9.6767　９。

658１9.62４8　　９.6４57

angle　=

Cｏlumns　1　ｔhrouｇｈ８

　　0　　　0.30110．4０11０.42５0　０。

443８　０．3４310．3942　０。

3899

　Coluｍns９tｈｒough13

　０.42６６　　0.42０9０。

4498　　　0.48２2　0．4452

PL　=

　0.０487

maｘd=

　　0。

0０1０

k=

　6

PQ1=

　　０.０29６　　0。

0262

　　　0。

0299　0。

02４3

0。

1４４4　0．12７3

　0.02７2　0．０２30

　0。

06780。

0581

　0。

13７90。

1231

　０.1１８３０。

0964

　　　0.１6６１　0.1379

0。

092０　0.0864

0.１８91　　0.1750

　　0。

38520.344５

　0．８116　　0。

７278

V=

　Cｏlumns1throｕgh　8

　　１0。

4０00　　　9．8７9５　9。

699１9.6653　　　9.6３81　　９.8４92　　９．７8739.７３26

　Ｃoｌumns９through1３

　９.６８709。

6７64９。

6579　9。

6238　　9。

６4４7

anglｅ=

Columｎｓ1through8

　0　　0。

30１1　0。

4０110.42500。

443８0.343１０.39４2　0.38９9

　Cｏluｍns9thｒougｈ　1３

0.４2６６　0。

42０90。

449８　0。

48２3　0。

４４５2

PL=

0.０4８８

ｍaｘd=

2。

6０21e—0０４

k＝

　７

PＱ1=

　0.029６0。

０２６２

　0。

02９９0.０２43

　０．1444　0.１2７3

0．０27２　0。

0２３０

　　0。

０６78　　0.058１

0。

13７9　0．1231

　０。

11８3０．０964

　　0。

1６６1０。

1３79

　0.09２０　　0.0８64

0.189１0.1750

　　0。

38５2　0.3445

　０.８11６０。

7278

V＝

　Columｎs1ｔhrｏｕgｈ　8

10.４000　９。

8７959。

699１9。

6６５29。

63７8　９。

84９29。

78729。

7326

　Columns　９tｈrough1３

　　９。

６８70　９.6764　9。

657８　9．６235　9。

6４45

angle＝

Columns1thｒouｇh８

　　０　　　0。

30１1０.４0110．425０　0.4438　0。

34３1　0.３942　0.389９

Coluｍnｓ9tｈｒoｕｇh13

０．4266　０.4２0９　0。

449８　0．４8230。

445２

ＰＬ=

　　０。

04８8

mａｘd=

　６。

1０４6e—00５

k　＝

　8

PQ１　＝

０．０296　０。

0２62

0。

０２99　0.0２43

0。

1444　0。

1273

0.０2７２0．0230

0。

067８　0.0581

　0.１379　　0.12３１

　　　0．118３0。

０９6４

　　0．１６61　0.1379

0。

０９２0　0．086４

　　0。

1８9１０.1750

　0。

3852　　0.３44５

　0。

8116　0。

7278

V＝

　Cｏluｍns１　thｒｏｕｇh8

10.4000　９。

87９５9。

699１　　9。

6６52９.6377　　　９。

84９2　　9。

7８7２　　9．7３2６

Ｃoｌuｍns9thｒough１3

9．6870　　9.6７６4　　　9。

65７89.6２3５9。

6444

ａngle=

Columｎs　１throuｇｈ8

　00.30110。

4011　0．４２5０　0.４43８0。

34３1　　０。

3942　0.38９9

　Ｃolumns9through13

　　0。

4266　　0.４209　0。

4498　0。

4823　　0.4452

ＰＬ＝

　０。

04８８

计算结果清单:

mａxd　=1

k=１

maxｄ=０。

178０

ｋ=　2

mａxd=　０.1７87

k=3

ｍaｘd=0.1７９3

k＝４

maxｄ＝０。

12２6

k=　5

maｘd=0。

０010

k＝6

ｍaxd＝2．60２1e—0４

k　=７

mａxd=6．1０４6e—05

k＝８

ＰQ1=

　0．0２96　　０。

０262

0.029９0。

0243

　0。

144４　　　０．12７3

0.02720。

0230

0。

0678　０．058１

　0。

1３７９　　0。

12３1

　　０。

1１830。

0964

　0．1661　0。

1３7９

　　0.09200。

0８６４

0。

1８９1　０。

１75０

　0.38５20.34４5

　0。

8１1６０。

7278

Ｖ=

10。

4000９.８79５9。

69９1　　9。

６652　　　9.637７　9。

8492９。

7８72　9。

７３26　　　　9。

68７09。

6７64　9。

６57８　9.６２35　　　9。

６４44

anｇｌe＝

　　0　　０.３0110。

4011　0.4２500.44380。

3431　0。

39４2　0。

３8９９　　0.４266　０。

4２090。

44９8　0.４82３　　０。

4452

ＰL=　0.０488

　　参考文献

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

ﻩ［２］　李维波．《ＭAＴLAB在电气工程中应用》.中国电力出版社。

2007