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()
%% MATPOWER Case Format : Version 2
mpc.version = '2';
%%----- Power Flow Data -----%%
%% system MVA base
mpc.baseMVA = 100;
%% 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;
];
%% 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
