-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathcase39mod3.m
executable file
·142 lines (136 loc) · 6.91 KB
/
case39mod3.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
function mpc = case39mod3
%% Modified 39 bus case in MATPOWER format
% W. A. Bukhsh, Feb 2013
%% References:
% [1] W. A. Bukhsh, Andreas Grothey, Ken McKinnon, Paul trodden, "Local Solutions of Optimal Power Flow Problem"
% submitted to IEEE Transactions on Power Systems, 2013
% [2] W. A. Bukhsh, Andreas Grothey, Ken McKinnon, Paul trodden, "Local Solutions of Optimal Power Flow Problem"
% Technical Report ERGO, 2011
%% 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 1 97.6 44.2 0 0 2 1.0393836 -13.536602 345 1 1.06 0.94;
2 1 0 0 0 0 2 1.0484941 -9.7852666 345 1 1.06 0.94;
3 1 322 2.4 0 0 2 1.0307077 -12.276384 345 1 1.06 0.94;
4 1 500 184 0 0 1 1.00446 -12.626734 345 1 1.06 0.94;
5 1 0 0 0 0 1 1.0060063 -11.192339 345 1 1.06 0.94;
6 1 0 0 0 0 1 1.0082256 -10.40833 345 1 1.06 0.94;
7 1 233.8 84 0 0 1 0.99839728 -12.755626 345 1 1.06 0.94;
8 1 522 176.6 0 0 1 0.99787232 -13.335844 345 1 1.06 0.94;
9 1 6.5 -66.6 0 0 1 1.038332 -14.178442 345 1 1.06 0.94;
10 1 0 0 0 0 1 1.0178431 -8.170875 345 1 1.06 0.94;
11 1 0 0 0 0 1 1.0133858 -8.9369663 345 1 1.06 0.94;
12 1 8.53 88 0 0 1 1.000815 -8.9988236 345 1 1.06 0.94;
13 1 0 0 0 0 1 1.014923 -8.9299272 345 1 1.06 0.94;
14 1 0 0 0 0 1 1.012319 -10.715295 345 1 1.06 0.94;
15 1 320 153 0 0 3 1.0161854 -11.345399 345 1 1.06 0.94;
16 1 329 32.3 0 0 3 1.0325203 -10.033348 345 1 1.06 0.94;
17 1 0 0 0 0 2 1.0342365 -11.116436 345 1 1.06 0.94;
18 1 158 30 0 0 2 1.0315726 -11.986168 345 1 1.06 0.94;
19 1 0 0 0 0 3 1.0501068 -5.4100729 345 1 1.06 0.94;
20 1 680 103 0 0 3 0.99101054 -6.8211783 345 1 1.06 0.94;
21 1 274 115 0 0 3 1.0323192 -7.6287461 345 1 1.06 0.94;
22 1 0 0 0 0 3 1.0501427 -3.1831199 345 1 1.06 0.94;
23 1 247.5 84.6 0 0 3 1.0451451 -3.3812763 345 1 1.06 0.94;
24 1 308.6 -92.2 0 0 3 1.038001 -9.9137585 345 1 1.06 0.94;
25 1 224 47.2 0 0 2 1.0576827 -8.3692354 345 1 1.06 0.94;
26 1 139 17 0 0 2 1.0525613 -9.4387696 345 1 1.06 0.94;
27 1 281 75.5 0 0 2 1.0383449 -11.362152 345 1 1.06 0.94;
28 1 206 27.6 0 0 3 1.0503737 -5.9283592 345 1 1.06 0.94;
29 1 283.5 26.9 0 0 3 1.0501149 -3.1698741 345 1 1.06 0.94;
30 2 0 0 0 0 2 1.0499 -7.3704746 345 1 1.06 0.94;
31 3 9.2 4.6 0 0 1 0.982 0 345 1 1.06 0.94;
32 2 0 0 0 0 1 0.9841 -0.1884374 345 1 1.06 0.94;
33 2 0 0 0 0 3 0.9972 -0.19317445 345 1 1.06 0.94;
34 2 0 0 0 0 3 1.0123 -1.631119 345 1 1.06 0.94;
35 2 0 0 0 0 3 1.0494 1.7765069 345 1 1.06 0.94;
36 2 0 0 0 0 3 1.0636 4.4684374 345 1 1.06 0.94;
37 2 0 0 0 0 2 1.0275 -1.5828988 345 1 1.06 0.94;
38 2 0 0 0 0 3 1.0265 3.8928177 345 1 1.06 0.94;
39 2 1104 250 0 0 1 1.03 -14.535256 345 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 = [
30 250 161.762 400*4 140*4 1.0499 100 1 1040*4 0 0 0 0 0 0 0 0 0 0 0 0;
31 677.871 221.574 300*4 -100*4 0.982 100 1 646*4 0 0 0 0 0 0 0 0 0 0 0 0;
32 650 206.965 300*4 150*4 0.9841 100 1 725*4 0 0 0 0 0 0 0 0 0 0 0 0;
33 632 108.293 250*4 0 0.9972 100 1 652*4 0 0 0 0 0 0 0 0 0 0 0 0;
34 508 166.688 167*4 0 1.0123 100 1 508*4 0 0 0 0 0 0 0 0 0 0 0 0;
35 650 210.661 300*4 -100*4 1.0494 100 1 687*4 0 0 0 0 0 0 0 0 0 0 0 0;
36 560 100.165 240*4 0 1.0636 100 1 580*4 0 0 0 0 0 0 0 0 0 0 0 0;
37 540 -1.36945 250*4 0 1.0275 100 1 564*4 0 0 0 0 0 0 0 0 0 0 0 0;
38 830 21.7327 300*4 -150*4 1.0265 100 1 865*4 0 0 0 0 0 0 0 0 0 0 0 0;
39 1000 78.4674 300*4 -100*4 1.03 100 1 1100*4 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.0035 0.0411 0.6987 600*4 600*4 600*4 0 0 1 -360 360;
1 39 0.001 0.025 0.75 1000*4 1000*4 1000*4 0 0 1 -360 360;
2 3 0.0013 0.0151 0.2572 500*4 500*4 500*4 0 0 1 -360 360;
2 25 0.007 0.0086 0.146 500*4 500*4 500*4 0 0 1 -360 360;
2 30 0 0.0181 0 900*4 900*4 2500*4 1.025 0 1 -360 360;
3 4 0.0013 0.0213 0.2214 500*4 500*4 500*4 0 0 1 -360 360;
3 18 0.0011 0.0133 0.2138 500*4 500*4 500*4 0 0 1 -360 360;
4 5 0.0008 0.0128 0.1342 600*4 600*4 600*4 0 0 1 -360 360;
4 14 0.0008 0.0129 0.1382 500*4 500*4 500*4 0 0 1 -360 360;
5 6 0.0002 0.0026 0.0434 1200*4 1200*4 1200*4 0 0 1 -360 360;
5 8 0.0008 0.0112 0.1476 900*4 900*4 900*4 0 0 1 -360 360;
6 7 0.0006 0.0092 0.113 900*4 900*4 900*4 0 0 1 -360 360;
6 11 0.0007 0.0082 0.1389 480*4 480*4 480*4 0 0 1 -360 360;
6 31 0 0.025 0 1800*4 1800*4 1800*4 1.07 0 1 -360 360;
7 8 0.0004 0.0046 0.078 900*4 900*4 900*4 0 0 1 -360 360;
8 9 0.0023 0.0363 0.3804 900*4 900*4 900*4 0 0 1 -360 360;
9 39 0.001 0.025 1.2 900*4 900*4 900*4 0 0 1 -360 360;
10 11 0.0004 0.0043 0.0729 600*4 600*4 600*4 0 0 1 -360 360;
10 13 0.0004 0.0043 0.0729 600*4 600*4 600*4 0 0 1 -360 360;
10 32 0 0.02 0 900*4 900*4 2500*4 1.07 0 1 -360 360;
12 11 0.0016 0.0435 0 500*4 500*4 500*4 1.006 0 1 -360 360;
12 13 0.0016 0.0435 0 500*4 500*4 500*4 1.006 0 1 -360 360;
13 14 0.0009 0.0101 0.1723 600*4 600*4 600*4 0 0 1 -360 360;
14 15 0.0018 0.0217 0.366 600*4 600*4 600*4 0 0 1 -360 360;
15 16 0.0009 0.0094 0.171 600*4 600*4 600*4 0 0 1 -360 360;
16 17 0.0007 0.0089 0.1342 600*4 600*4 600*4 0 0 1 -360 360;
16 19 0.0016 0.0195 0.304 600*4 600*4 2500*4 0 0 1 -360 360;
16 21 0.0008 0.0135 0.2548 600*4 600*4 600*4 0 0 1 -360 360;
16 24 0.0003 0.0059 0.068 600*4 600*4 600*4 0 0 1 -360 360;
17 18 0.0007 0.0082 0.1319 600*4 600*4 600*4 0 0 1 -360 360;
17 27 0.0013 0.0173 0.3216 600*4 600*4 600*4 0 0 1 -360 360;
19 20 0.0007 0.0138 0 900*4 900*4 2500*4 1.06 0 1 -360 360;
19 33 0.0007 0.0142 0 900*4 900*4 2500*4 1.07 0 1 -360 360;
20 34 0.0009 0.018 0 900*4 900*4 2500*4 1.009 0 1 -360 360;
21 22 0.0008 0.014 0.2565 900*4 900*4 900*4 0 0 1 -360 360;
22 23 0.0006 0.0096 0.1846 600*4 600*4 600*4 0 0 1 -360 360;
22 35 0 0.0143 0 900*4 900*4 2500*4 1.025 0 1 -360 360;
23 24 0.0022 0.035 0.361 600*4 600*4 600*4 0 0 1 -360 360;
23 36 0.0005 0.0272 0 900*4 900*4 2500*4 1 0 1 -360 360;
25 26 0.0032 0.0323 0.531 600*4 600*4 600*4 0 0 1 -360 360;
25 37 0.0006 0.0232 0 900*4 900*4 2500*4 1.025 0 1 -360 360;
26 27 0.0014 0.0147 0.2396 600*4 600*4 600*4 0 0 1 -360 360;
26 28 0.0043 0.0474 0.7802 600*4 600*4 600*4 0 0 1 -360 360;
26 29 0.0057 0.0625 1.029 600*4 600*4 600*4 0 0 1 -360 360;
28 29 0.0014 0.0151 0.249 600*4 600*4 600*4 0 0 1 -360 360;
29 38 0.0008 0.0156 0 1200*4 1200*4 2500*4 1.025 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.01*0 0.3 0.2*0;
2 0 0 3 0.01*0 0.3 0.2*0;
2 0 0 3 0.01*0 0.3 0.2*0;
2 0 0 3 0.01*0 0.3 0.2*0;
2 0 0 3 0.01*0 0.3 0.2*0;
2 0 0 3 0.01*0 0.3 0.2*0;
2 0 0 3 0.01*0 0.3 0.2*0;
2 0 0 3 0.01*0 0.3 0.2*0;
2 0 0 3 0.01*0 0.3 0.2*0;
2 0 0 3 0.01*0 0.3 0.2*0;
];