-
Notifications
You must be signed in to change notification settings - Fork 84
/
pglib_opf_case30_as.m
197 lines (191 loc) · 10.8 KB
/
pglib_opf_case30_as.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
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%% %%%%%
%%%% IEEE PES Power Grid Library - Optimal Power Flow - v23.07 %%%%%
%%%% (https://github.com/power-grid-lib/pglib-opf) %%%%%
%%%% Benchmark Group - Typical Operations %%%%%
%%%% 23 - July - 2023 %%%%%
%%%% %%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Power flow data for 30 bus, 6 generator case.
%
% 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.
%
% The generation MVA limit is used for Q upper bounds
% Shunts assumed to be in real value (not p.u. as specified in the paper)
%
% Copyright (c) 1974 by The Institute of Electrical and Electronics Engineers (IEEE)
% Licensed under the Creative Commons Attribution 4.0
% International license, http://creativecommons.org/licenses/by/4.0/
%
% Contact M.E. Brennan ([email protected]) for inquries on further reuse of
% this dataset.
%
function mpc = pglib_opf_case30_as
mpc.version = '2';
mpc.baseMVA = 100.0;
%% area data
% area refbus
mpc.areas = [
1 1;
];
%% bus data
% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin
mpc.bus = [
1 3 0.0 0.0 0.0 0.0 1 1.00000 0.00000 135.0 1 1.05000 0.95000;
2 2 21.7 12.7 0.0 0.0 1 1.02500 0.00000 135.0 1 1.10000 0.95000;
3 1 2.4 1.2 0.0 0.0 1 1.00000 0.00000 135.0 1 1.05000 0.95000;
4 1 7.6 1.6 0.0 0.0 1 1.00000 0.00000 135.0 1 1.05000 0.95000;
5 1 94.2 19.0 0.0 0.0 1 1.00000 0.00000 135.0 1 1.05000 0.95000;
6 1 0.0 0.0 0.0 0.0 1 1.00000 0.00000 135.0 1 1.05000 0.95000;
7 1 22.8 10.9 0.0 0.0 1 1.00000 0.00000 135.0 1 1.05000 0.95000;
8 1 30.0 30.0 0.0 0.0 1 1.00000 0.00000 135.0 1 1.05000 0.95000;
9 1 0.0 0.0 0.0 0.0 1 1.00000 0.00000 135.0 1 1.05000 0.95000;
10 1 5.8 2.0 0.0 5.26 1 1.00000 0.00000 135.0 1 1.05000 0.95000;
11 1 0.0 0.0 0.0 0.0 1 1.00000 0.00000 135.0 1 1.05000 0.95000;
12 1 11.2 7.5 0.0 0.0 1 1.00000 0.00000 135.0 1 1.05000 0.95000;
13 2 0.0 0.0 0.0 0.0 1 1.02500 0.00000 135.0 1 1.10000 0.95000;
14 1 6.2 1.6 0.0 0.0 1 1.00000 0.00000 135.0 1 1.05000 0.95000;
15 1 8.2 2.5 0.0 0.0 1 1.00000 0.00000 135.0 1 1.05000 0.95000;
16 1 3.5 1.8 0.0 0.0 1 1.00000 0.00000 135.0 1 1.05000 0.95000;
17 1 9.0 5.8 0.0 0.0 1 1.00000 0.00000 135.0 1 1.05000 0.95000;
18 1 3.2 0.9 0.0 0.0 1 1.00000 0.00000 135.0 1 1.05000 0.95000;
19 1 9.5 3.4 0.0 0.0 1 1.00000 0.00000 135.0 1 1.05000 0.95000;
20 1 2.2 0.7 0.0 0.0 1 1.00000 0.00000 135.0 1 1.05000 0.95000;
21 1 17.5 11.2 0.0 0.0 1 1.00000 0.00000 135.0 1 1.05000 0.95000;
22 2 0.0 0.0 0.0 0.0 1 1.02500 0.00000 135.0 1 1.10000 0.95000;
23 2 3.2 1.6 0.0 0.0 1 1.02500 0.00000 135.0 1 1.10000 0.95000;
24 1 8.7 6.7 0.0 25.0 1 1.00000 0.00000 135.0 1 1.05000 0.95000;
25 1 0.0 0.0 0.0 0.0 1 1.00000 0.00000 135.0 1 1.05000 0.95000;
26 1 3.5 2.3 0.0 0.0 1 1.00000 0.00000 135.0 1 1.05000 0.95000;
27 2 0.0 0.0 0.0 0.0 1 1.02500 0.00000 135.0 1 1.10000 0.95000;
28 1 0.0 0.0 0.0 0.0 1 1.00000 0.00000 135.0 1 1.05000 0.95000;
29 1 2.4 0.9 0.0 0.0 1 1.00000 0.00000 135.0 1 1.05000 0.95000;
30 1 10.6 1.9 0.0 0.0 1 1.00000 0.00000 135.0 1 1.05000 0.95000;
];
%% generator data
% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin
mpc.gen = [
1 125.0 115.0 250.0 -20.0 1.0 100.0 1 200.0 50.0;
2 50.0 40.0 100.0 -20.0 1.025 100.0 1 80.0 20.0;
5 32.5 32.5 80.0 -15.0 1.0 100.0 1 50.0 15.0;
8 22.5 22.5 60.0 -15.0 1.0 100.0 1 35.0 10.0;
11 20.0 20.0 50.0 -10.0 1.0 100.0 1 30.0 10.0;
13 26.0 22.5 60.0 -15.0 1.025 100.0 1 40.0 12.0;
];
%% generator cost data
% 2 startup shutdown n c(n-1) ... c0
mpc.gencost = [
2 0.0 0.0 3 0.003750 2.000000 0.000000;
2 0.0 0.0 3 0.017500 1.750000 0.000000;
2 0.0 0.0 3 0.062500 1.000000 0.000000;
2 0.0 0.0 3 0.008340 3.250000 0.000000;
2 0.0 0.0 3 0.025000 3.000000 0.000000;
2 0.0 0.0 3 0.025000 3.000000 0.000000;
];
%% branch data
% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax
mpc.branch = [
1 2 0.0192 0.0575 0.0264 130.0 130.0 130.0 0.0 0.0 1 -30.0 30.0;
1 3 0.0452 0.1852 0.0204 130.0 130.0 130.0 0.0 0.0 1 -30.0 30.0;
2 4 0.057 0.1737 0.0184 65.0 65.0 65.0 0.0 0.0 1 -30.0 30.0;
3 4 0.0132 0.0379 0.0042 130.0 130.0 130.0 0.0 0.0 1 -30.0 30.0;
2 5 0.0472 0.1983 0.0209 130.0 130.0 130.0 0.0 0.0 1 -30.0 30.0;
2 6 0.0581 0.1763 0.0187 65.0 65.0 65.0 0.0 0.0 1 -30.0 30.0;
4 6 0.0119 0.0414 0.0045 90.0 90.0 90.0 0.0 0.0 1 -30.0 30.0;
5 7 0.046 0.116 0.0102 70.0 70.0 70.0 0.0 0.0 1 -30.0 30.0;
6 7 0.0267 0.082 0.0085 130.0 130.0 130.0 0.0 0.0 1 -30.0 30.0;
6 8 0.012 0.042 0.0045 32.0 32.0 32.0 0.0 0.0 1 -30.0 30.0;
6 9 0.0 0.208 0.0 65.0 65.0 65.0 0.0 0.0 1 -30.0 30.0;
6 10 0.0 0.556 0.0 32.0 32.0 32.0 0.0 0.0 1 -30.0 30.0;
9 11 0.0 0.208 0.0 65.0 65.0 65.0 0.0 0.0 1 -30.0 30.0;
9 10 0.0 0.11 0.0 65.0 65.0 65.0 0.0 0.0 1 -30.0 30.0;
4 12 0.0 0.256 0.0 65.0 65.0 65.0 0.0 0.0 1 -30.0 30.0;
12 13 0.0 0.14 0.0 65.0 65.0 65.0 0.0 0.0 1 -30.0 30.0;
12 14 0.1231 0.2559 0.0 32.0 32.0 32.0 0.0 0.0 1 -30.0 30.0;
12 15 0.0662 0.1304 0.0 32.0 32.0 32.0 0.0 0.0 1 -30.0 30.0;
12 16 0.0945 0.1987 0.0 32.0 32.0 32.0 0.0 0.0 1 -30.0 30.0;
14 15 0.221 0.1997 0.0 16.0 16.0 16.0 0.0 0.0 1 -30.0 30.0;
16 17 0.0824 0.1932 0.0 16.0 16.0 16.0 0.0 0.0 1 -30.0 30.0;
15 18 0.107 0.2185 0.0 16.0 16.0 16.0 0.0 0.0 1 -30.0 30.0;
18 19 0.0639 0.1292 0.0 16.0 16.0 16.0 0.0 0.0 1 -30.0 30.0;
19 20 0.034 0.068 0.0 32.0 32.0 32.0 0.0 0.0 1 -30.0 30.0;
10 20 0.0936 0.209 0.0 32.0 32.0 32.0 0.0 0.0 1 -30.0 30.0;
10 17 0.0324 0.0845 0.0 32.0 32.0 32.0 0.0 0.0 1 -30.0 30.0;
10 21 0.0348 0.0749 0.0 32.0 32.0 32.0 0.0 0.0 1 -30.0 30.0;
10 22 0.0727 0.1499 0.0 32.0 32.0 32.0 0.0 0.0 1 -30.0 30.0;
21 22 0.0116 0.0236 0.0 32.0 32.0 32.0 0.0 0.0 1 -30.0 30.0;
15 23 0.1 0.202 0.0 16.0 16.0 16.0 0.0 0.0 1 -30.0 30.0;
22 24 0.115 0.179 0.0 16.0 16.0 16.0 0.0 0.0 1 -30.0 30.0;
23 24 0.132 0.27 0.0 16.0 16.0 16.0 0.0 0.0 1 -30.0 30.0;
24 25 0.1885 0.3292 0.0 16.0 16.0 16.0 0.0 0.0 1 -30.0 30.0;
25 26 0.2544 0.38 0.0 16.0 16.0 16.0 0.0 0.0 1 -30.0 30.0;
25 27 0.1093 0.2087 0.0 16.0 16.0 16.0 0.0 0.0 1 -30.0 30.0;
28 27 0.0 0.396 0.0 65.0 65.0 65.0 0.0 0.0 1 -30.0 30.0;
27 29 0.2198 0.4153 0.0 16.0 16.0 16.0 0.0 0.0 1 -30.0 30.0;
27 30 0.3202 0.6027 0.0 16.0 16.0 16.0 0.0 0.0 1 -30.0 30.0;
29 30 0.2399 0.4533 0.0 16.0 16.0 16.0 0.0 0.0 1 -30.0 30.0;
8 28 0.0636 0.2 0.0214 32.0 32.0 32.0 0.0 0.0 1 -30.0 30.0;
6 28 0.0169 0.0599 0.0065 32.0 32.0 32.0 0.0 0.0 1 -30.0 30.0;
];
% INFO : === Translation Options ===
% INFO : Phase Angle Bound: 30.0 (deg.)
% INFO : Setting Flat Start
% INFO :
% INFO : === Generator Bounds Update Notes ===
% INFO :
% INFO : === Base KV Replacement Notes ===
% INFO :
% INFO : === Transformer Setting Replacement Notes ===
% INFO :
% INFO : === Line Capacity Monotonicity Notes ===
% INFO :
% INFO : === Voltage Setpoint Replacement Notes ===
% INFO : Bus 1 : V=1.0, theta=0.0 -> V=1.0, theta=0.0
% INFO : Bus 2 : V=1.0, theta=0.0 -> V=1.025, theta=0.0
% INFO : Bus 3 : V=1.0, theta=0.0 -> V=1.0, theta=0.0
% INFO : Bus 4 : V=1.0, theta=0.0 -> V=1.0, theta=0.0
% INFO : Bus 5 : V=1.0, theta=0.0 -> V=1.0, theta=0.0
% INFO : Bus 6 : V=1.0, theta=0.0 -> V=1.0, theta=0.0
% INFO : Bus 7 : V=1.0, theta=0.0 -> V=1.0, theta=0.0
% INFO : Bus 8 : V=1.0, theta=0.0 -> V=1.0, theta=0.0
% INFO : Bus 9 : V=1.0, theta=0.0 -> V=1.0, theta=0.0
% INFO : Bus 10 : V=1.0, theta=0.0 -> V=1.0, theta=0.0
% INFO : Bus 11 : V=1.0, theta=0.0 -> V=1.0, theta=0.0
% INFO : Bus 12 : V=1.0, theta=0.0 -> V=1.0, theta=0.0
% INFO : Bus 13 : V=1.0, theta=0.0 -> V=1.025, theta=0.0
% INFO : Bus 14 : V=1.0, theta=0.0 -> V=1.0, theta=0.0
% INFO : Bus 15 : V=1.0, theta=0.0 -> V=1.0, theta=0.0
% INFO : Bus 16 : V=1.0, theta=0.0 -> V=1.0, theta=0.0
% INFO : Bus 17 : V=1.0, theta=0.0 -> V=1.0, theta=0.0
% INFO : Bus 18 : V=1.0, theta=0.0 -> V=1.0, theta=0.0
% INFO : Bus 19 : V=1.0, theta=0.0 -> V=1.0, theta=0.0
% INFO : Bus 20 : V=1.0, theta=0.0 -> V=1.0, theta=0.0
% INFO : Bus 21 : V=1.0, theta=0.0 -> V=1.0, theta=0.0
% INFO : Bus 22 : V=1.0, theta=0.0 -> V=1.025, theta=0.0
% INFO : Bus 23 : V=1.0, theta=0.0 -> V=1.025, theta=0.0
% INFO : Bus 24 : V=1.0, theta=0.0 -> V=1.0, theta=0.0
% INFO : Bus 25 : V=1.0, theta=0.0 -> V=1.0, theta=0.0
% INFO : Bus 26 : V=1.0, theta=0.0 -> V=1.0, theta=0.0
% INFO : Bus 27 : V=1.0, theta=0.0 -> V=1.025, theta=0.0
% INFO : Bus 28 : V=1.0, theta=0.0 -> V=1.0, theta=0.0
% INFO : Bus 29 : V=1.0, theta=0.0 -> V=1.0, theta=0.0
% INFO : Bus 30 : V=1.0, theta=0.0 -> V=1.0, theta=0.0
% INFO :
% INFO : === Generator Setpoint Replacement Notes ===
% INFO : Gen at bus 1 : Pg=23.54, Qg=0.0 -> Pg=125.0, Qg=115.0
% INFO : Gen at bus 1 : Vg=1.0 -> Vg=1.0
% INFO : Gen at bus 2 : Pg=60.97, Qg=0.0 -> Pg=50.0, Qg=40.0
% INFO : Gen at bus 2 : Vg=1.0 -> Vg=1.025
% INFO : Gen at bus 5 : Pg=21.59, Qg=0.0 -> Pg=32.5, Qg=32.5
% INFO : Gen at bus 5 : Vg=1.0 -> Vg=1.0
% INFO : Gen at bus 8 : Pg=26.91, Qg=0.0 -> Pg=22.5, Qg=22.5
% INFO : Gen at bus 8 : Vg=1.0 -> Vg=1.0
% INFO : Gen at bus 11 : Pg=19.2, Qg=0.0 -> Pg=20.0, Qg=20.0
% INFO : Gen at bus 11 : Vg=1.0 -> Vg=1.0
% INFO : Gen at bus 13 : Pg=37.0, Qg=0.0 -> Pg=26.0, Qg=22.5
% INFO : Gen at bus 13 : Vg=1.0 -> Vg=1.025
% INFO :
% INFO : === Writing Matpower Case File Notes ===