Re: PF and OPF did not converge

Carlos E Murillo-Sanchez Wed, 04 Oct 2017 15:49:24 -0700

Yes, both of them. Some reactances are larger than unity ... that's a lot for a branch. And for the OPF, you certainly need the slack generator with much larger PMAX. I tried a simple test reducing all branch parameters by a tenfold and the power flow converged.

Carlos.

Charalambos Ioannou wrote:

Dear Carlos,

So you are suggesting to check again the loads and also to reduce the reactances of the branches ?

Moreover to define a value for PMAX of the slack bus ?

Thanks

On 4 Oct 2017, at 18:12, Carlos E Murillo-Sanchez <[email protected]> wrote:

Dear Charalambos:

To begin with, you have a load of 894.31 MW (computed with sum(mpc.bus(:,PD)) and a total generation capacity of 30 (computed with sum(mpc.gen(:, PMAX)) . Your slack generator has PMAX specified as zero. And it seems to me that many of your branch reactances are way too high to allow a poser flow solution with a total load of 8.1 p.u. ,

Carlos.

Charalambos Ioannou wrote:

Hi Abhyankar,

I already did that you send before but i am facing the same problems.

Thanks

Regards

Charalampos

Get Outlook for iOS

From: [email protected] <[email protected]> on behalf of Abhyankar, Shrirang G. <[email protected]>
Sent: Tuesday, October 3, 2017 3:43:04 PM
To: MATPOWER discussion forum
Subject: Re: PF and OPF did not converge

http://www.pserc.cornell.edu/matpower/#pfconvergence

From: <[email protected]> on behalf of Charalambos Ioannou <[email protected]>
Reply-To: MATPOWER discussion forum <[email protected]>
Date: Tuesday, October 3, 2017 at 9:12 AM
To: "[email protected]" <[email protected]>
Subject: PF and OPF did not converge

Dear Ray,

I am trying to run PF and OPF for the following case but every time it doesnt converged the system. Please can you give me any advice??

The system consist 62 buses, 23 transformers, 8 photo voltaic plants and 1 Natural gus generator.

I consider all the buses as Type 1(PQ) and the natural gas generator as Type 2 (PV). The real power demand (Pd) for the photo voltaic plants i gave them a negative sign.

function mpc = IsleExample3()

TCL MERGE ERROR ( 10/03/2017 10:44:16 ): "invalid command name "MATPOWER"" OutmailID: 121902634, List: 'matpower-l', MemberID: 79270638 SCRIPT: "MATPOWER Case Format : Version 2

mpc.version = '2';

" ----- Power Flow Data -----TCL MERGE ERROR ( 10/03/2017 10:44:16 ): "invalid command name "
"" OutmailID: 121902634, List: 'matpower-l', MemberID: 79270638 SCRIPT: "

" system MVA base

mpc.baseMVA = 100;

TCL MERGE ERROR ( 10/03/2017 10:44:16 ): "invalid command name "bus"" OutmailID: 121902634, List: 'matpower-l', MemberID: 79270638 SCRIPT: "bus data

% bus_i type Pd Qd Gs Bs

area Vm

Va baseKV

zone Vmax

Vmin

mpc.bus = [

1 3 30.4 14.7 0 0 1 1 0 132 1 1.06 0.94;

2 1 30.4 14.7 0 0 1 1 0 132 1 1.06 0.94;

3 1 22.8 22.8 0 0 1 1 0 132 1 1.06 0.94;

4 1 22.8 22.8 0 0 1 1 0 132 1 1.06 0.94;

5 1 113 55 0 0 1 1 0 132 1 1.06 0.94;

6 1 113 55 0 0 1 1 0 132 1 1.06 0.94;

7 1 0 0 0 0 1 1 0 132 1 1.06 0.94;

8 1 113 55 0 0 1 1 0 132 1 1.06 0.94;

9 1 113 55 0 0 1 1 0 132 1 1.06 0.94;

10 1 113 55 0 0 1 1 0 33 1 1.06 0.94;

11 1 19 9 0 0 1 1 0 11 1 1.06 0.94;

12 1 113 55 0 0 1 1 0 33 1 1.06 0.94;

13 1 -6.89 0 0 0 1 1 0 33 1 1.06 0.94;

14 1 -5.00 0 0 0 1 1 0 33 1 1.06 0.94;

15 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

16 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

17 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

18 1 6 3 0 0 1 1 0 11 1 1.06 0.94;

19 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

20 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

21 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

22 1 -9.20 0 0 0 1 1 0 33 1 1.06 0.94;

23 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

24 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

25 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

26 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

27 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

28 1 17 8 0 0 1 1 0 11 1 1.06 0.94;

29 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

30 1 -5.00 0 0 0 1 1 0 33 1 1.06 0.94;

31 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

32 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

33 1 -5.00 0 0 0 1 1 0 33 1 1.06 0.94;

34 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

35 1 6 3 0 0 1 1 0 11 1 1.06 0.94;

36 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

37 1 -7.00 0 0 0 1 1 0 33 1 1.06 0.94;

38 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

39 1 6 3 0 0 1 1 0 11 1 1.06 0.94;

40 1 -1.5 0 0 0 1 1 0 33 1 1.06 0.94;

41 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

42 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

43 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

44 1 27 13 0 0 1 1 0 11 1 1.06 0.94;

45 1 27 13 0 0 1 1 0 11 1 1.06 0.94;

46 1 -5.00 0 0 0 1 1 0 33 1 1.06 0.94;

47 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

48 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

49 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

50 1 8 4 0 0 1 1 0 11 1 1.06 0.94;

51 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

52 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

53 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

54 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

55 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

56 1 12 6 0 0 1 1 0 11 1 1.06 0.94;

57 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

58 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

59 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

60 1 11 5 0 0 1 1 0 11 1 1.06 0.94;

61 1 0 0 0 0 1 1 0 33 1 1.06 0.94;

62 2 15.5 0 0 0 1 1 0 33 1 1.06 0.94;

];

" Generator Data

% bus Pg Qg Qmax

Qmin Vg

mBase status

Pmax Pmin

Pc1 Pc2

Qc1min Qc1max

Qc2min Qc2max

ramp_agc ramp_10

ramp_30 ramp_q

apf

mpc.gen = [

1 0 0 500 -500 1 100 1 0 0 0 0 0 0 0 0 0 0 0 0 0;

62 15.5 0 30 -30 1 100 1 30 0 0 0 0 0 0 0 0 0 0 0 0;

];

TCL MERGE ERROR ( 10/03/2017 10:44:16 ): "invalid command name "branch"" OutmailID: 121902634, List: 'matpower-l', MemberID: 79270638 SCRIPT: "branch data

% fbus tbus r x b rateA

rateB rateC

ratio angle

status angmin

angmax

mpc.branch = [

1 2 0.0012 0.0221 0.0000 99 0 0 0 0 1 -360 360;

1 3 0.0019 0.0232 0.0000 138 0 0 0 0 1 -360 360;

2 4 0.0019 0.0232 0.0000 138 0 0 0 0 1 -360 360;

5 2 0.0012 0.0221 0.0000 124 0 0 0 0 1 -360 360;

3 8 0.0146 0.0303 0.0830 124 0 0 0 0 1 -360 360;

4 9 0.0145 0.0303 0.0808 124 0 0 0 0 1 -360 360;

5 6 0.0039 0.0277 0.0000 45 0 0 0 0 1 -360 360;

5 10 0 0.2771 0 45 0 0 0.95 0 1 -360 360;

6 10 0 0.2251 0 45 0 0 0.95 0 1 -360 360;

6 7 0.0019 0.0101 0.0174 145 0 0 0 0 1 -360 360;

6 7 0.0026 0.0078 0.0184 145 0 0 0 0 1 -360 360;

7 8 0.0040 0.0281 0.0000 0 0 0 0 0 1 -360 360;

7 9 0.0040 0.0281 0.0000 0 0 0 0 0 1 -360 360;

8 12 0 0.2501 0 90 0 0 0.95 0 1 -360 360;

9 12 0 0.1986 0 90 0 0 0.95 0 1 -360 360;

12 13 0.0083 0.0073 0.0006 24 0 0 0 0 1 -360 360;

12 14 0.0262 0.0231 0.0020 27 0 0 0 0 1 -360 360;

12 15 0.0389 0.1302 0.0002 44 0 0 0 0 1 -360 360;

15 16 0.0008 0.0012 0.0002 35 0 0 0 0 1 -360 360;

16 18 0 1.0800 0 15 0 0 0.955 0 1 -360 360;

15 20 0.0315 0.1054 0.0001 44 0 0 0 0 1 -360 360;

20 21 0.0584 0.0667 0.0001 19 0 0 0 0 1 -360 360;

21 22 0.0238 0.0210 0.0018 27 0 0 0 0 1 -360 360;

23 21 0.0262 0.0299 0.0000 19 0 0 0 0 1 -360 360;

27 20 0.0093 0.0310 0.0000 44 0 0 0 0 1 -360 360;

27 28 0 0.9930 0 30 0 0 0.955 0 1 -360 360;

12 19 0.0398 0.1333 0.0002 44 0 0 0 0 1 -360 360;

19 17 0.0016 0.0031 0.0000 31 0 0 0 0 1 -360 360;

17 18 0 1.0730 0 15 0 0 0.955 0 1 -360 360;

19 24 0.0398 0.1333 0.0002 44 0 0 0 0 1 -360 360;

25 24 0.0027 0.0020 0.0001 12 0 0 0 0 1 -360 360;

26 24 0.0093 0.0310 0.0000 44 0 0 0 0 1 -360 360;

26 28 0 0.9930 0 30 0 0 0.955 0 1 -360 360;

10 11 0 1.0300 0 30 0 0 0.955 0 1 -360 360;

10 11 0 1.0270 0 30 0 0 0.955 0 1 -360 360;

10 12 0.0944 0.1853 0.0002 31 0 0 0 0 1 -360 360;

10 12 0.0912 0.1790 0.0002 31 0 0 0 0 1 -360 360;

10 29 0.0087 0.0037 0.0002 31 0 0 0 0 1 -360 360;

10 29 0.0087 0.0037 0.0002 31 0 0 0 0 1 -360 360;

29 30 0.1005 0.2384 0.0010 30 0 0 0 0 1 -360 360;

31 30 0.0672 0.1601 0.0002 30 0 0 0 0 1 -360 360;

31 32 0.0320 0.0628 0.0001 31 0 0 0 0 1 -360 360;

32 33 0.0011 0.0010 0.0001 27 0 0 0 0 1 -360 360;

34 32 0.0844 0.1555 0.0012 31 0 0 0 0 1 -360 360;

34 35 0 1.0100 0 15 0 0 0.955 0 1 -360 360;

31 38 0.0542 0.2397 0.0000 44 0 0 0 0 1 -360 360;

31 39 0 0.9910 0 15 0 0 1 0 1 -360 360;

38 39 0 0.9910 0 15 0 0 1 0 1 -360 360;

38 37 0.0843 0.1655 0.0002 31 0 0 0 0 1 -360 360;

37 36 0.0261 0.0379 0.0015 31 0 0 0 0 1 -360 360;

36 35 0 1.0140 0 15 0 0 0.955 0 1 -360 360;

38 40 0.0012 0.0010 0.0001 24 0 0 0 0 1 -360 360;

38 41 0.0268 0.3173 0.0004 21 0 0 0 0 1 -360 360;

12 41 0.0445 0.0730 0.0018 31 0 0 0 0 1 -360 360;

41 44 0 0.7210 0 30 0 0 0.945 0 1 -360 360;

12 42 0.0397 0.0683 0.0035 31 0 0 0 0 1 -360 360;

42 45 0 0.7760 0 30 0 0 0.955 0 1 -360 360;

12 47 0.0481 0.1612 0.0002 44 0 0 0 0 1 -360 360;

47 46 0.0199 0.0368 0.0003 31 0 0 0 0 1 -360 360;

43 47 0.0586 0.1024 0.0029 25 0 0 0 0 1 -360 360;

43 45 0 0.7210 0 30 0 0 0.955 0 1 -360 360;

47 48 0.0667 0.1289 0.0004 31 0 0 0 0 1 -360 360;

48 49 0.1870 0.2211 0.0003 20 0 0 0 0 1 -360 360;

49 50 0 0.4670 0 30 0 0 0.955 0 1 -360 360;

48 57 0.0368 0.0722 0.0001 31 0 0 0 0 1 -360 360;

57 56 0 1.0920 0 30 0 0 0.96 0 1 -360 360;

12 54 0.1168 0.2780 0.0003 30 0 0 0 0 1 -360 360;

61 54 0.0676 0.1313 0.0006 31 0 0 0 0 1 -360 360;

61 60 0 1.0510 0 15 0 0 0.95 0 1 -360 360;

54 55 0.0278 0.0662 0.0001 30 0 0 0 0 1 -360 360;

55 56 0 1.0920 0 30 0 0 0.96 0 1 -360 360;

12 53 0.0629 0.2108 0.0003 44 0 0 0 0 1 -360 360;

53 62 0.0011 0.0022 0.0002 44 0 0 0 0 1 -360 360;

52 53 0.0111 0.0372 0.0000 44 0 0 0 0 1 -360 360;

52 51 0.1980 0.2343 0.0003 20 0 0 0 0 1 -360 360;

51 50 0 0.4770 0 30 0 0 0.955 0 1 -360 360;

52 58 0.0789 0.1688 0.0007 31 0 0 0 0 1 -360 360;

59 58 0.0016 0.0030 0.0000 27 0 0 0 0 1 -360 360;

58 60 0 1.0510 0 15 0 0 0.95 0 1 -360 360;

];

" Generator Cost Data

% model Startup Shutdown

n x1 y1 y2

mpc.gencost = [

2 0 0 3 0 1 0;

2 0 0 3 0 1 0;

];

Thanks in Advance.

Kind Regards ,

Charalampos

Re: PF and OPF did not converge

Reply via email to