Dear Prof.Ray,
The runpf and the runopf did not converge while handling a revised system
(i.e., IEEE 30 and IEEE 33). The revised system is attached. Could you give
me some suggestions when you feel convenient? Thanks.
The revising details: a standard IEEE 30 bus system [] is revised by
connecting the bus 2 to the bus 1 of an IEEE 33-bus radial system [] which
consists of 33 buses, 32 lines, a voltage of 12.66kV, load size of
3.715MW and 2.3MVar.
In the new combined system, the original bus-index of the IEEE 33-bus has
increased 30 (i.e., the original bus-index 1 corresponds the new index 31, and
so on). The resistance and reactance between bus 1 and bus 31 are 0.2030 (p.u.)
and 0.1034 (p.u.) respectively.
Regards,
Chixin
function mpc = case30 %CASE30 Power flow data for 30 bus, 6 generator case. % Please see CASEFORMAT for details on the case file format. % % Based on data from ... % Alsac, O. & Stott, B., "Optimal Load Flow with Steady State Security", % IEEE Transactions on Power Apparatus and Systems, Vol. PAS 93, No. 3, % 1974, pp. 745-751. % ... with branch parameters rounded to nearest 0.01, shunt values divided % by 100 and shunt on bus 10 moved to bus 5, load at bus 5 zeroed out. % Generator locations, costs and limits and bus areas were taken from ... % Ferrero, R.W., Shahidehpour, S.M., Ramesh, V.C., "Transaction analysis % in deregulated power systems using game theory", IEEE Transactions on % Power Systems, Vol. 12, No. 3, Aug 1997, pp. 1340-1347. % Generator Q limits were derived from Alsac & Stott, using their Pmax % capacities. V limits and line |S| limits taken from Alsac & Stott.
% MATPOWER
%% 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 0 0 0 0 1 1 0
135 1 1.05 0.95;
2 2 21.7 12.7 0 0 1 1 0
135 1 1.1 0.95;
3 1 2.4 1.2 0 0 1 1 0
135 1 1.05 0.95;
4 1 7.6 1.6 0 0 1 1 0
135 1 1.05 0.95;
5 1 0 0 0 0.19 1 1 0
135 1 1.05 0.95;
6 1 0 0 0 0 1 1 0
135 1 1.05 0.95;
7 1 22.8 10.9 0 0 1 1 0
135 1 1.05 0.95;
8 1 30 30 0 0 1 1 0
135 1 1.05 0.95;
9 1 0 0 0 0 1 1 0
135 1 1.05 0.95;
10 1 5.8 2 0 0 3 1 0
135 1 1.05 0.95;
11 1 0 0 0 0 1 1 0
135 1 1.05 0.95;
12 1 11.2 7.5 0 0 2 1 0
135 1 1.05 0.95;
13 2 0 0 0 0 2 1 0
135 1 1.1 0.95;
14 1 6.2 1.6 0 0 2 1 0
135 1 1.05 0.95;
15 1 8.2 2.5 0 0 2 1 0
135 1 1.05 0.95;
16 1 3.5 1.8 0 0 2 1 0
135 1 1.05 0.95;
17 1 9 5.8 0 0 2 1 0
135 1 1.05 0.95;
18 1 3.2 0.9 0 0 2 1 0
135 1 1.05 0.95;
19 1 9.5 3.4 0 0 2 1 0
135 1 1.05 0.95;
20 1 2.2 0.7 0 0 2 1 0
135 1 1.05 0.95;
21 1 17.5 11.2 0 0 3 1 0
135 1 1.05 0.95;
22 2 0 0 0 0 3 1 0
135 1 1.1 0.95;
23 2 3.2 1.6 0 0 2 1 0
135 1 1.1 0.95;
24 1 8.7 6.7 0 0.04 3 1 0
135 1 1.05 0.95;
25 1 0 0 0 0 3 1 0
135 1 1.05 0.95;
26 1 3.5 2.3 0 0 3 1 0
135 1 1.05 0.95;
27 2 0 0 0 0 3 1 0
135 1 1.1 0.95;
28 1 0 0 0 0 1 1 0
135 1 1.05 0.95;
29 1 2.4 0.9 0 0 3 1 0
135 1 1.05 0.95;
30 1 10.6 1.9 0 0 3 1 0
135 1 1.05 0.95;
%The following data is appended from case 33
31 1 0 0 0 0 1 1 0 12.66
1 1 1;
32 1 100 60 0 0 1 1 0
12.66 1 1.1 0.9;
33 1 90 40 0 0 1 1 0
12.66 1 1.1 0.9;
34 1 120 80 0 0 1 1 0
12.66 1 1.1 0.9;
35 1 60 30 0 0 1 1 0
12.66 1 1.1 0.9;
36 1 60 20 0 0 1 1 0
12.66 1 1.1 0.9;
37 1 200 100 0 0 1 1 0
12.66 1 1.1 0.9;
38 1 200 100 0 0 1 1 0
12.66 1 1.1 0.9;
39 1 60 20 0 0 1 1 0
12.66 1 1.1 0.9;
40 1 60 20 0 0 1 1 0
12.66 1 1.1 0.9;
41 1 45 30 0 0 1 1 0
12.66 1 1.1 0.9;
42 1 60 35 0 0 1 1 0
12.66 1 1.1 0.9;
43 1 60 35 0 0 1 1 0
12.66 1 1.1 0.9;
44 1 120 80 0 0 1 1 0
12.66 1 1.1 0.9;
45 1 60 10 0 0 1 1 0
12.66 1 1.1 0.9;
46 1 60 20 0 0 1 1 0
12.66 1 1.1 0.9;
47 1 60 20 0 0 1 1 0
12.66 1 1.1 0.9;
48 1 90 40 0 0 1 1 0
12.66 1 1.1 0.9;
49 1 90 40 0 0 1 1 0
12.66 1 1.1 0.9;
50 1 90 40 0 0 1 1 0
12.66 1 1.1 0.9;
51 1 90 40 0 0 1 1 0
12.66 1 1.1 0.9;
52 1 90 40 0 0 1 1 0
12.66 1 1.1 0.9;
53 1 90 50 0 0 1 1 0
12.66 1 1.1 0.9;
54 1 420 200 0 0 1 1 0
12.66 1 1.1 0.9;
55 1 420 200 0 0 1 1 0
12.66 1 1.1 0.9;
56 1 60 25 0 0 1 1 0
12.66 1 1.1 0.9;
57 1 60 25 0 0 1 1 0
12.66 1 1.1 0.9;
58 1 60 20 0 0 1 1 0
12.66 1 1.1 0.9;
59 1 120 70 0 0 1 1 0
12.66 1 1.1 0.9;
60 1 200 600 0 0 1 1 0
12.66 1 1.1 0.9;
61 1 150 70 0 0 1 1 0
12.66 1 1.1 0.9;
62 1 210 100 0 0 1 1 0
12.66 1 1.1 0.9;
63 1 60 40 0 0 1 1 0
12.66 1 1.1 0.9;
];
%% 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 23.54 0 150 -20 1 100 1 80
0 0 0 0 0 0 0 0 0 0
0 0;
2 60.97 0 60 -20 1 100 1 80
0 0 0 0 0 0 0 0 0 0
0 0;
22 21.59 0 62.5 -15 1 100 1 50
0 0 0 0 0 0 0 0 0 0
0 0;
27 26.91 0 48.7 -15 1 100 1 55
0 0 0 0 0 0 0 0 0 0
0 0;
23 19.2 0 40 -10 1 100 1 30
0 0 0 0 0 0 0 0 0 0
0 0;
13 37 0 44.7 -15 1 100 1 40
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.02 0.06 0.03 130 130 130 0
0 1 -360 360;
1 3 0.05 0.19 0.02 130 130 130 0
0 1 -360 360;
2 4 0.06 0.17 0.02 65 65 65 0
0 1 -360 360;
3 4 0.01 0.04 0 130 130 130 0
0 1 -360 360;
2 5 0.05 0.2 0.02 130 130 130 0
0 1 -360 360;
2 6 0.06 0.18 0.02 65 65 65 0
0 1 -360 360;
4 6 0.01 0.04 0 90 90 90 0
0 1 -360 360;
5 7 0.05 0.12 0.01 70 70 70 0
0 1 -360 360;
6 7 0.03 0.08 0.01 130 130 130 0
0 1 -360 360;
6 8 0.01 0.04 0 32 32 32 0
0 1 -360 360;
6 9 0 0.21 0 65 65 65 0
0 1 -360 360;
6 10 0 0.56 0 32 32 32 0
0 1 -360 360;
9 11 0 0.21 0 65 65 65 0
0 1 -360 360;
9 10 0 0.11 0 65 65 65 0
0 1 -360 360;
4 12 0 0.26 0 65 65 65 0
0 1 -360 360;
12 13 0 0.14 0 65 65 65 0
0 1 -360 360;
12 14 0.12 0.26 0 32 32 32 0
0 1 -360 360;
12 15 0.07 0.13 0 32 32 32 0
0 1 -360 360;
12 16 0.09 0.2 0 32 32 32 0
0 1 -360 360;
14 15 0.22 0.2 0 16 16 16 0
0 1 -360 360;
16 17 0.08 0.19 0 16 16 16 0
0 1 -360 360;
15 18 0.11 0.22 0 16 16 16 0
0 1 -360 360;
18 19 0.06 0.13 0 16 16 16 0
0 1 -360 360;
19 20 0.03 0.07 0 32 32 32 0
0 1 -360 360;
10 20 0.09 0.21 0 32 32 32 0
0 1 -360 360;
10 17 0.03 0.08 0 32 32 32 0
0 1 -360 360;
10 21 0.03 0.07 0 32 32 32 0
0 1 -360 360;
10 22 0.07 0.15 0 32 32 32 0
0 1 -360 360;
21 22 0.01 0.02 0 32 32 32 0
0 1 -360 360;
15 23 0.1 0.2 0 16 16 16 0
0 1 -360 360;
22 24 0.12 0.18 0 16 16 16 0
0 1 -360 360;
23 24 0.13 0.27 0 16 16 16 0
0 1 -360 360;
24 25 0.19 0.33 0 16 16 16 0
0 1 -360 360;
25 26 0.25 0.38 0 16 16 16 0
0 1 -360 360;
25 27 0.11 0.21 0 16 16 16 0
0 1 -360 360;
28 27 0 0.4 0 65 65 65 0
0 1 -360 360;
27 29 0.22 0.42 0 16 16 16 0
0 1 -360 360;
27 30 0.32 0.6 0 16 16 16 0
0 1 -360 360;
29 30 0.24 0.45 0 16 16 16 0
0 1 -360 360;
8 28 0.06 0.2 0.02 32 32 32 0
0 1 -360 360;
6 28 0.02 0.06 0.01 32 32 32 0
0 1 -360 360;
%The following data is appended from case 33(i.e. branch 1 to branch 32)
2 31 0.2030 0.1034 0 0 0 0 0 0
1 -360 360;
31 32 0.0922 0.0470 0 0 0 0 0 0
1 -360 360;
32 33 0.4930 0.2511 0 0 0 0 0
0 1 -360 360;
33 34 0.3660 0.1864 0 0 0 0 0
0 1 -360 360;
34 35 0.3811 0.1941 0 0 0 0 0
0 1 -360 360;
35 36 0.8190 0.7070 0 0 0 0 0
0 1 -360 360;
36 37 0.1872 0.6188 0 0 0 0 0
0 1 -360 360;
37 38 0.7114 0.2351 0 0 0 0 0
0 1 -360 360;
38 39 1.0300 0.7400 0 0 0 0 0
0 1 -360 360;
39 40 1.0440 0.7400 0 0 0 0 0
0 1 -360 360;
40 41 0.1966 0.0650 0 0 0 0 0
0 1 -360 360;
41 42 0.3744 0.1238 0 0 0 0 0
0 1 -360 360;
42 43 1.4680 1.1550 0 0 0 0 0
0 1 -360 360;
43 44 0.5416 0.7129 0 0 0 0 0
0 1 -360 360;
44 45 0.5910 0.5260 0 0 0 0 0
0 1 -360 360;
45 46 0.7463 0.5450 0 0 0 0 0
0 1 -360 360;
46 47 1.2890 1.7210 0 0 0 0 0
0 1 -360 360;
47 48 0.7320 0.5740 0 0 0 0 0
0 1 -360 360;
32 49 0.1640 0.1565 0 0 0 0 0
0 1 -360 360;
49 50 1.5042 1.3554 0 0 0 0 0
0 1 -360 360;
50 51 0.4095 0.4784 0 0 0 0 0
0 1 -360 360;
51 52 0.7089 0.9373 0 0 0 0 0
0 1 -360 360;
33 53 0.4512 0.3083 0 0 0 0 0
0 1 -360 360;
53 54 0.8980 0.7091 0 0 0 0 0
0 1 -360 360;
54 55 0.8960 0.7011 0 0 0 0 0
0 1 -360 360;
36 56 0.2030 0.1034 0 0 0 0 0
0 1 -360 360;
56 57 0.2842 0.1447 0 0 0 0 0
0 1 -360 360;
57 58 1.0590 0.9337 0 0 0 0 0
0 1 -360 360;
58 59 0.8042 0.7006 0 0 0 0 0
0 1 -360 360;
59 60 0.5075 0.2585 0 0 0 0 0
0 1 -360 360;
60 61 0.9744 0.9630 0 0 0 0 0
0 1 -360 360;
61 62 0.3105 0.3619 0 0 0 0 0
0 1 -360 360;
62 63 0.3410 0.5302 0 0 0 0 0
0 1 -360 360;
];
%%----- OPF Data -----%%
%% generator cost data
% 1 startup shutdown n x1 y1 ... xn
yn
% 2 startup shutdown n c(n-1) ... c0
mpc.gencost = [
2 0 0 3 0.02 2 0;
2 0 0 3 0.0175 1.75 0;
2 0 0 3 0.0625 1 0;
2 0 0 3 0.00834 3.25 0;
2 0 0 3 0.025 3 0;
2 0 0 3 0.025 3 0;
];
