To Whom It May Concern,
I built the network of RBTS and it is shown below.
Now, the OPF of the network is running correctly.
However, when I disconnected the line 70 which links Supply Point 1 (SP1) and
Supply Point 2 (SP2), the SP2 and SP3 cannot get the power supply.
The SP2 and SP3 suppose to get supply due to the line 71 which links SP2 and
SP3 at that condition.
Because the power can flow to SP3 via line 68 and 69 from SP1, then the SP2 can
get power via line 71 at that condition.
Do you know what happened to the network and how can I fix it ?
I list the case lines at the end of content body.
Best Regards,
Hui
[cid:a8f7b1a4-2893-4e1c-946d-68212f2fbd64]
function mpc = caseRBTS
%CASEcase33kv-11kv Power flow data of Distribution System for.RBTS BUS 4
%% MATPOWER Case Format : Version 4.1
mpc.version = '2.0';
%%----- 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 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
2 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94;
3 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
4 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94;
5 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
6 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94;
7 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
8 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94;
9 1 0.8137 0.2675 0 0 1 1 360 11 1 1.06 0.94;
10 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
11 1 0.6714 0.2207 0 0 1 1 360 11 1 1.06 0.94;
12 1 0.6714 0.2207 0 0 1 1 360 11 1 1.06 0.94;
13 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
14 1 1.63 0.5358 0 0 1 1 360 11 1 1.06 0.94;
15 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
16 1 2.445 0.8036 0 0 1 1 360 11 1 1.06 0.94;
17 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
18 1 1.63 0.5358 0 0 1 1 360 11 1 1.06 0.94;
19 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
20 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94;
21 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
22 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94;
23 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
24 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94;
25 1 0.8137 0.2675 0 0 1 1 360 11 1 1.06 0.94;
26 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
27 1 0.8137 0.2675 0 0 1 1 360 11 1 1.06 0.94;
28 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
29 1 0.6714 0.2207 0 0 1 1 360 11 1 1.06 0.94;
30 1 0.6714 0.2207 0 0 1 1 360 11 1 1.06 0.94;
31 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
32 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94;
33 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
34 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94;
35 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94;
36 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
37 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94;
38 1 0.8137 0.2675 0 0 1 1 360 11 1 1.06 0.94;
39 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
40 1 0.8137 0.2675 0 0 1 1 360 11 1 1.06 0.94;
41 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
42 1 0.6714 0.2207 0 0 1 1 360 11 1 1.06 0.94;
43 1 0.6714 0.2207 0 0 1 1 360 11 1 1.06 0.94;
44 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
45 1 1.63 0.5358 0 0 1 1 360 11 1 1.06 0.94;
46 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
47 1 1.63 0.5358 0 0 1 1 360 11 1 1.06 0.94;
48 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
49 1 1.63 0.5358 0 0 1 1 360 11 1 1.06 0.94;
50 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
51 1 1.63 0.5358 0 0 1 1 360 11 1 1.06 0.94;
52 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
53 1 1.63 0.5358 0 0 1 1 360 11 1 1.06 0.94;
54 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
55 1 2.445 0.8036 0 0 1 1 360 11 1 1.06 0.94;
56 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
57 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94;
58 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
59 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94;
60 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
61 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94;
62 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94;
63 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
64 1 0.8137 0.2675 0 0 1 1 360 11 1 1.06 0.94;
65 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
66 1 0.8137 0.2675 0 0 1 1 360 11 1 1.06 0.94;
67 1 0.6714 0.2207 0 0 1 1 360 11 1 1.06 0.94;
68 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
69 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
70 1 0 0 0 0 1 1 360 11 1 1.06 0.94;
71 3 0 0 0 0 1 1 360 33 1 1.03 0.97;
72 1 0 0 0 0 1 1 360 33 1 1.06 0.94;
73 1 0 0 0 0 1 1 360 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 = [
2 0.8869 0.2915 0.2915 -0.2915 1.00 100 1 0.8869 0 0 0 0 0
0 0 0 0 0 0 0;
4 0.8869 0.2915 0.2915 -0.2915 1.00 100 1 0.8869 0 0 0 0 0
0 0 0 0 0 0 0;
6 0.8869 0.2915 0.2915 -0.2915 1.00 100 1 0.8869 0 0 0 0 0
0 0 0 0 0 0 0;
8 0.8869 0.2915 0.2915 -0.2915 1.00 100 1 0.8869 0 0 0 0 0
0 0 0 0 0 0 0;
9 0.8137 0.2675 0.2675 -0.2675 1.00 100 1 0.8137 0 0 0 0 0
0 0 0 0 0 0 0;
11 0.6714 0.2207 0.2207 -0.2207 1.00 100 1 0.6714 0 0 0 0 0 0
0 0 0 0 0 0;
12 0.6714 0.2207 0.2207 -0.2207 1.00 100 1 0.6714 0 0 0 0 0 0
0 0 0 0 0 0;
14 1.63 0.5358 0.5358 -0.5358 1.00 100 1 1.63 0 0 0 0 0 0
0 0 0 0 0 0;
16 2.455 0.8036 0.8036 -0.8036 1.00 100 1 2.455 0 0 0 0 0 0
0 0 0 0 0 0;
18 1.63 0.5358 0.5358 -0.5358 1.00 100 1 1.63 0 0 0 0 0
0 0 0 0 0 0 0;
20 0.8869 0.2915 0.2915 -0.2915 1.00 100 1 0.8869 0 0 0 0 0 0
0 0 0 0 0 0;
22 0.8869 0.2915 0.2915 -0.2915 1.00 100 1 0.8869 0 0 0 0 0 0
0 0 0 0 0 0;
24 0.8869 0.2915 0.2915 -0.2915 1.00 100 1 0.8869 0 0 0 0 0 0
0 0 0 0 0 0;
25 0.8137 0.2675 0.2675 -0.2675 1.00 100 1 0.8137 0 0 0 0 0 0
0 0 0 0 0 0;
27 0.8137 0.2675 0.2675 -0.2675 1.00 100 1 0.8137 0 0 0 0 0 0
0 0 0 0 0 0;
29 0.6714 0.2207 0.2207 -0.2207 1.00 100 1 0.6714 0 0 0 0 0 0
0 0 0 0 0 0;
30 0.6714 0.2207 0.2207 -0.2207 1.00 100 1 0.6714 0 0 0 0 0 0
0 0 0 0 0 0;
32 0.8869 0.2915 0.2915 -0.2915 1.00 100 1 0.8869 0 0 0 0 0 0
0 0 0 0 0 0;
34 0.8869 0.2915 0.2915 -0.2915 1.00 100 1 0.8869 0 0 0 0 0 0
0 0 0 0 0 0;
35 0.8869 0.2915 0.2915 -0.2915 1.00 100 1 0.8869 0 0 0 0 0 0
0 0 0 0 0 0;
37 0.8869 0.2915 0.2915 -0.2915 1.00 100 1 0.8869 0 0 0 0 0 0
0 0 0 0 0 0;
38 0.8137 0.2675 0.2675 -0.2675 1.00 100 1 0.8137 0 0 0 0 0 0
0 0 0 0 0 0;
40 0.8137 0.2675 0.2675 -0.2675 1.00 100 1 0.8137 0 0 0 0 0 0
0 0 0 0 0 0;
42 0.6714 0.2207 0.2207 -0.2207 1.00 100 1 0.6714 0 0 0 0 0 0
0 0 0 0 0 0;
43 0.6714 0.2207 0.2207 -0.2207 1.00 100 1 0.6714 0 0 0 0 0 0
0 0 0 0 0 0;
45 1.63 0.5358 0.5358 -0.5358 1.00 100 1 1.63 0 0 0 0 0 0
0 0 0 0 0 0;
47 1.63 0.5358 0.5358 -0.5358 1.00 100 1 1.63 0 0 0 0 0 0
0 0 0 0 0 0;
49 1.63 0.5358 0.5358 -0.5358 1.00 100 1 1.63 0 0 0 0 0
0 0 0 0 0 0 0;
51 1.63 0.5358 0.5358 -0.5358 1.00 100 1 1.63 0 0 0 0 0
0 0 0 0 0 0 0;
53 1.63 0.5358 0.5358 -0.5358 1.00 100 1 1.63 0 0 0 0 0
0 0 0 0 0 0 0;
55 2.455 0.8036 0.8036 -0.8036 1.00 100 1 2.455 0 0 0 0 0 0
0 0 0 0 0 0;
57 0.8869 0.2915 0.2915 -0.2915 1.00 100 1 0.8869 0 0 0 0 0 0
0 0 0 0 0 0;
59 0.8869 0.2915 0.2915 -0.2915 1.00 100 1 0.8869 0 0 0 0 0 0
0 0 0 0 0 0;
61 0.8869 0.2915 0.2915 -0.2915 1.00 100 1 0.8869 0 0 0 0 0 0
0 0 0 0 0 0;
62 0.8869 0.2915 0.2915 -0.2915 1.00 100 1 0.8869 0 0 0 0 0 0
0 0 0 0 0 0;
64 0.8137 0.2675 0.2675 -0.2675 1.00 100 1 0.8137 0 0 0 0 0 0
0 0 0 0 0 0;
66 0.8137 0.2675 0.2675 -0.2675 1.00 100 1 0.8137 0 0 0 0 0 0
0 0 0 0 0 0;
67 0.6714 0.2207 0.2207 -0.2207 1.00 100 1 0.6714 0 0 0 0 0 0
0 0 0 0 0 0;
71 40 17 17 -17 1.02 100 1 41
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 = [
68 1 0.062 0.057 0 8.193 0 0 0 0 1 -360 360;
1 2 0.135 0.06 0 4.858 0 0 0 0 1 -360 360;
1 3 0.083 0.063 0 7.240 0 0 0 0 1 -360 360;
3 4 0.168 0.073 0 4.858 0 0 0 0 1 -360 360;
3 5 0.062 0.047 0 7.240 0 0 0 0 1 -360 360;
5 6 0.179 0.077 0 4.858 0 0 0 0 1 -360 360;
5 7 0.166 0.066 0 4.858 0 0 0 0 1 -360 360;
7 8 0.179 0.077 0 4.858 0 0 0 0 1 -360 360;
7 9 0.168 0.073 0 4.858 0 0 0 0 1 -360 360;
7 10 0.133 0.053 0 4.858 0 0 0 0 1 -360 360;
10 11 0.18 0.080 0 4.858 0 0 0 0 1 -360 360;
10 12 0.168 0.075 0 4.858 0 0 0 0 1 -360 360;
68 13 0.066 0.061 0 8.193 0 0 0 0 1 -360 360;
13 14 0.133 0.053 0 4.858 0 0 0 0 1 -360 360;
13 15 0.102 0.064 0 6.478 0 0 0 0 1 -360 360;
15 16 0.166 0.066 0 4.858 0 0 0 0 1 -360 360;
15 17 0.133 0.053 0 4.858 0 0 0 0 1 -360 360;
17 18 0.177 0.071 0 4.858 0 0 0 0 1 -360 360;
68 19 0.062 0.057 0 8.193 0 0 0 0 1 -360 360;
19 20 0.179 0.077 0 4.858 0 0 0 0 1 -360 360;
19 21 0.062 0.047 0 7.240 0 0 0 0 1 -360 360;
21 22 0.168 0.073 0 4.858 0 0 0 0 1 -360 360;
21 23 0.102 0.064 0 5.716 0 0 0 0 1 -360 360;
23 24 0.168 0.073 0 4.858 0 0 0 0 1 -360 360;
23 25 0.135 0.060 0 4.858 0 0 0 0 1 -360 360;
23 26 0.177 0.071 0 4.858 0 0 0 0 1 -360 360;
26 27 0.168 0.073 0 4.858 0 0 0 0 1 -360 360;
26 28 0.133 0.053 0 4.858 0 0 0 0 1 -360 360;
28 29 0.168 0.075 0 4.858 0 0 0 0 1 -360 360;
28 30 0.135 0.062 0 4.858 0 0 0 0 1 -360 360;
69 31 0.052 0.059 0 9.526 0 0 0 0 1 -360 360;
31 32 0.168 0.073 0 4.858 0 0 0 0 1 -360 360;
31 33 0.066 0.061 0 8.193 0 0 0 0 1 -360 360;
33 34 0.135 0.06 0 4.858 0 0 0 0 1 -360 360;
33 35 0.168 0.073 0 4.858 0 0 0 0 1 -360 360;
33 36 0.102 0.064 0 5.716 0 0 0 0 1 -360 360;
36 37 0.168 0.073 0 4.858 0 0 0 0 1 -360 360;
36 38 0.135 0.06 0 4.858 0 0 0 0 1 -360 360;
36 39 0.177 0.071 0 4.858 0 0 0 0 1 -360 360;
39 40 0.168 0.073 0 4.858 0 0 0 0 1 -360 360;
39 41 0.133 0.053 0 4.858 0 0 0 0 1 -360 360;
41 42 0.168 0.075 0 4.858 0 0 0 0 1 -360 360;
41 43 0.135 0.062 0 4.858 0 0 0 0 1 -360 360;
69 44 0.083 0.063 0 7.240 0 0 0 0 1 -360 360;
44 45 0.166 0.066 0 4.858 0 0 0 0 1 -360 360;
44 46 0.076 0.048 0 5.716 0 0 0 0 1 -360 360;
46 47 0.18 0.08 0 4.858 0 0 0 0 1 -360 360;
46 48 0.166 0.066 0 4.858 0 0 0 0 1 -360 360;
48 49 0.133 0.053 0 4.858 0 0 0 0 1 -360 360;
70 50 0.062 0.057 0 8.193 0 0 0 0 1 -360 360;
50 51 0.133 0.053 0 4.858 0 0 0 0 1 -360 360;
50 52 0.102 0.064 0 6.478 0 0 0 0 1 -360 360;
52 53 0.166 0.066 0 4.858 0 0 0 0 1 -360 360;
52 54 0.177 0.071 0 4.858 0 0 0 0 1 -360 360;
54 55 0.133 0.053 0 4.858 0 0 0 0 1 -360 360;
70 56 0.062 0.057 0 8.193 0 0 0 0 1 -360 360;
56 57 0.179 0.077 0 4.858 0 0 0 0 1 -360 360;
56 58 0.062 0.047 0 7.240 0 0 0 0 1 -360 360;
58 59 0.179 0.077 0 4.858 0 0 0 0 1 -360 360;
58 60 0.095 0.06 0 6.478 0 0 0 0 1 -360 360;
60 61 0.135 0.06 0 4.858 0 0 0 0 1 -360 360;
60 62 0.179 0.077 0 4.858 0 0 0 0 1 -360 360;
60 63 0.166 0.066 0 4.858 0 0 0 0 1 -360 360;
63 64 0.135 0.06 0 4.858 0 0 0 0 1 -360 360;
63 65 0.166 0.066 0 4.858 0 0 0 0 1 -360 360;
65 66 0.179 0.077 0 4.858 0 0 0 0 1 -360 360;
65 67 0.133 0.053 0 4.858 0 0 0 0 1 -360 360;
71 73 0.138 0.152 0 24.006 0 0 0 0 1 -360 360;
71 72 0.092 0.101 0 0.0001 0 0 0 0 1 -360 360;
72 73 0.092 0.101 0 24.006 0 0 0 0 1 -360 360;
71 68 0.054 1.087 0 11.5 23 0 0 0 1 -360 360;
71 68 0.054 1.087 0 11.5 23 0 0 0 1 -360 360;
72 69 0.05 1 0 15 0 0 0 0 1 -360 360;
72 69 0.05 1 0 15 0 0 0 0 1 -360 360;
73 70 0.05 1 0 15 0 0 0 0 1 -360 360;
73 70 0.05 1 0 15 0 0 0 0 1 -360 360;
];
%% 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 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 1 10 6000;
2 0 0 3 0 0 200;?
];
