CASE39 Power flow data for 39 bus New England system.
Please see CASEFORMAT for details on the case file format.
Data taken from [1] with the following modifications/additions:
- renumbered gen buses consecutively (as in [2] and [4])
- added Pmin = 0 for all gens
- added Qmin, Qmax for gens at 31 & 39 (copied from gen at 35)
- added Vg based on V in bus data (missing for bus 39)
- added Vg, Pg, Pd, Qd at bus 39 from [2] (same in [4])
- added Pmax at bus 39: Pmax = Pg + 100
- added line flow limits and area data from [4]
- added voltage limits, Vmax = 1.06, Vmin = 0.94
- added identical quadratic generator costs
- increased Pmax for gen at bus 34 from 308 to 508
(assumed typo in [1], makes initial solved case feasible)
- re-solved power flow
Notes:
- Bus 39, its generator and 2 connecting lines were added
(by authors of [1]) to represent the interconnection with
the rest of the eastern interconnect, and did not include
Vg, Pg, Qg, Pd, Qd, Pmin, Pmax, Qmin or Qmax.
- As the swing bus, bus 31 did not include and Q limits.
- The voltages, etc in [1] appear to be quite close to the
power flow solution of the case before adding bus 39 with
it's generator and connecting branches, though the solution
is not exact.
- Explicit voltage setpoints for gen buses are not given, so
they are taken from the bus data, however this results in two
binding Q limits at buses 34 & 37, so the corresponding
voltages have probably deviated from their original setpoints.
- The generator locations and types are as follows:
1 30 hydro
2 31 nuke01
3 32 nuke02
4 33 fossil02
5 34 fossil01
6 35 nuke03
7 36 fossil04
8 37 nuke04
9 38 nuke05
10 39 interconnection to rest of US/Canada
This is a solved power flow case, but it includes the following
violations:
- Pmax violated at bus 31: Pg = 677.87, Pmax = 646
- Qmin violated at bus 37: Qg = -1.37, Qmin = 0
References:
[1] G. W. Bills, et.al., "On-Line Stability Analysis Study"
RP90-1 Report for the Edison Electric Institute, October 12, 1970,
pp. 1-20 - 1-35.
prepared by E. M. Gulachenski - New England Electric System
J. M. Undrill - General Electric Co.
"generally representative of the New England 345 KV system, but is
not an exact or complete model of any past, present or projected
configuration of the actual New England 345 KV system.
[2] M. A. Pai, Energy Function Analysis for Power System Stability,
Kluwer Academic Publishers, Boston, 1989.
(references [3] as source of data)
[3] Athay, T.; Podmore, R.; Virmani, S., "A Practical Method for the
Direct Analysis of Transient Stability," IEEE Transactions on Power
Apparatus and Systems , vol.PAS-98, no.2, pp.573-584, March 1979.
URL: http://ieeexplore.ieee.org/stamp/stamp.jsp?arnumber=4113518&isnumber=4113486
(references [1] as source of data)
[4] Data included with TC Calculator at http://www.pserc.cornell.edu/tcc/
for 39-bus system.0001 function mpc = case39 0002 %CASE39 Power flow data for 39 bus New England system. 0003 % Please see CASEFORMAT for details on the case file format. 0004 % 0005 % Data taken from [1] with the following modifications/additions: 0006 % 0007 % - renumbered gen buses consecutively (as in [2] and [4]) 0008 % - added Pmin = 0 for all gens 0009 % - added Qmin, Qmax for gens at 31 & 39 (copied from gen at 35) 0010 % - added Vg based on V in bus data (missing for bus 39) 0011 % - added Vg, Pg, Pd, Qd at bus 39 from [2] (same in [4]) 0012 % - added Pmax at bus 39: Pmax = Pg + 100 0013 % - added line flow limits and area data from [4] 0014 % - added voltage limits, Vmax = 1.06, Vmin = 0.94 0015 % - added identical quadratic generator costs 0016 % - increased Pmax for gen at bus 34 from 308 to 508 0017 % (assumed typo in [1], makes initial solved case feasible) 0018 % - re-solved power flow 0019 % 0020 % Notes: 0021 % - Bus 39, its generator and 2 connecting lines were added 0022 % (by authors of [1]) to represent the interconnection with 0023 % the rest of the eastern interconnect, and did not include 0024 % Vg, Pg, Qg, Pd, Qd, Pmin, Pmax, Qmin or Qmax. 0025 % - As the swing bus, bus 31 did not include and Q limits. 0026 % - The voltages, etc in [1] appear to be quite close to the 0027 % power flow solution of the case before adding bus 39 with 0028 % it's generator and connecting branches, though the solution 0029 % is not exact. 0030 % - Explicit voltage setpoints for gen buses are not given, so 0031 % they are taken from the bus data, however this results in two 0032 % binding Q limits at buses 34 & 37, so the corresponding 0033 % voltages have probably deviated from their original setpoints. 0034 % - The generator locations and types are as follows: 0035 % 1 30 hydro 0036 % 2 31 nuke01 0037 % 3 32 nuke02 0038 % 4 33 fossil02 0039 % 5 34 fossil01 0040 % 6 35 nuke03 0041 % 7 36 fossil04 0042 % 8 37 nuke04 0043 % 9 38 nuke05 0044 % 10 39 interconnection to rest of US/Canada 0045 % 0046 % This is a solved power flow case, but it includes the following 0047 % violations: 0048 % - Pmax violated at bus 31: Pg = 677.87, Pmax = 646 0049 % - Qmin violated at bus 37: Qg = -1.37, Qmin = 0 0050 % 0051 % References: 0052 % [1] G. W. Bills, et.al., "On-Line Stability Analysis Study" 0053 % RP90-1 Report for the Edison Electric Institute, October 12, 1970, 0054 % pp. 1-20 - 1-35. 0055 % prepared by E. M. Gulachenski - New England Electric System 0056 % J. M. Undrill - General Electric Co. 0057 % "generally representative of the New England 345 KV system, but is 0058 % not an exact or complete model of any past, present or projected 0059 % configuration of the actual New England 345 KV system. 0060 % [2] M. A. Pai, Energy Function Analysis for Power System Stability, 0061 % Kluwer Academic Publishers, Boston, 1989. 0062 % (references [3] as source of data) 0063 % [3] Athay, T.; Podmore, R.; Virmani, S., "A Practical Method for the 0064 % Direct Analysis of Transient Stability," IEEE Transactions on Power 0065 % Apparatus and Systems , vol.PAS-98, no.2, pp.573-584, March 1979. 0066 % URL: http://ieeexplore.ieee.org/stamp/stamp.jsp?arnumber=4113518&isnumber=4113486 0067 % (references [1] as source of data) 0068 % [4] Data included with TC Calculator at http://www.pserc.cornell.edu/tcc/ 0069 % for 39-bus system. 0070 0071 % MATPOWER 0072 0073 %% MATPOWER Case Format : Version 2 0074 mpc.version = '2'; 0075 0076 %%----- Power Flow Data -----%% 0077 %% system MVA base 0078 mpc.baseMVA = 100; 0079 0080 %% bus data 0081 % bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin 0082 mpc.bus = [ 0083 1 1 97.6 44.2 0 0 2 1.0393836 -13.536602 345 1 1.06 0.94; 0084 2 1 0 0 0 0 2 1.0484941 -9.7852666 345 1 1.06 0.94; 0085 3 1 322 2.4 0 0 2 1.0307077 -12.276384 345 1 1.06 0.94; 0086 4 1 500 184 0 0 1 1.00446 -12.626734 345 1 1.06 0.94; 0087 5 1 0 0 0 0 1 1.0060063 -11.192339 345 1 1.06 0.94; 0088 6 1 0 0 0 0 1 1.0082256 -10.40833 345 1 1.06 0.94; 0089 7 1 233.8 84 0 0 1 0.99839728 -12.755626 345 1 1.06 0.94; 0090 8 1 522 176.6 0 0 1 0.99787232 -13.335844 345 1 1.06 0.94; 0091 9 1 6.5 -66.6 0 0 1 1.038332 -14.178442 345 1 1.06 0.94; 0092 10 1 0 0 0 0 1 1.0178431 -8.170875 345 1 1.06 0.94; 0093 11 1 0 0 0 0 1 1.0133858 -8.9369663 345 1 1.06 0.94; 0094 12 1 8.53 88 0 0 1 1.000815 -8.9988236 345 1 1.06 0.94; 0095 13 1 0 0 0 0 1 1.014923 -8.9299272 345 1 1.06 0.94; 0096 14 1 0 0 0 0 1 1.012319 -10.715295 345 1 1.06 0.94; 0097 15 1 320 153 0 0 3 1.0161854 -11.345399 345 1 1.06 0.94; 0098 16 1 329 32.3 0 0 3 1.0325203 -10.033348 345 1 1.06 0.94; 0099 17 1 0 0 0 0 2 1.0342365 -11.116436 345 1 1.06 0.94; 0100 18 1 158 30 0 0 2 1.0315726 -11.986168 345 1 1.06 0.94; 0101 19 1 0 0 0 0 3 1.0501068 -5.4100729 345 1 1.06 0.94; 0102 20 1 680 103 0 0 3 0.99101054 -6.8211783 345 1 1.06 0.94; 0103 21 1 274 115 0 0 3 1.0323192 -7.6287461 345 1 1.06 0.94; 0104 22 1 0 0 0 0 3 1.0501427 -3.1831199 345 1 1.06 0.94; 0105 23 1 247.5 84.6 0 0 3 1.0451451 -3.3812763 345 1 1.06 0.94; 0106 24 1 308.6 -92.2 0 0 3 1.038001 -9.9137585 345 1 1.06 0.94; 0107 25 1 224 47.2 0 0 2 1.0576827 -8.3692354 345 1 1.06 0.94; 0108 26 1 139 17 0 0 2 1.0525613 -9.4387696 345 1 1.06 0.94; 0109 27 1 281 75.5 0 0 2 1.0383449 -11.362152 345 1 1.06 0.94; 0110 28 1 206 27.6 0 0 3 1.0503737 -5.9283592 345 1 1.06 0.94; 0111 29 1 283.5 26.9 0 0 3 1.0501149 -3.1698741 345 1 1.06 0.94; 0112 30 2 0 0 0 0 2 1.0499 -7.3704746 345 1 1.06 0.94; 0113 31 3 9.2 4.6 0 0 1 0.982 0 345 1 1.06 0.94; 0114 32 2 0 0 0 0 1 0.9841 -0.1884374 345 1 1.06 0.94; 0115 33 2 0 0 0 0 3 0.9972 -0.19317445 345 1 1.06 0.94; 0116 34 2 0 0 0 0 3 1.0123 -1.631119 345 1 1.06 0.94; 0117 35 2 0 0 0 0 3 1.0494 1.7765069 345 1 1.06 0.94; 0118 36 2 0 0 0 0 3 1.0636 4.4684374 345 1 1.06 0.94; 0119 37 2 0 0 0 0 2 1.0275 -1.5828988 345 1 1.06 0.94; 0120 38 2 0 0 0 0 3 1.0265 3.8928177 345 1 1.06 0.94; 0121 39 2 1104 250 0 0 1 1.03 -14.535256 345 1 1.06 0.94; 0122 ]; 0123 0124 %% generator data 0125 % 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 0126 mpc.gen = [ 0127 30 250 161.762 400 140 1.0499 100 1 1040 0 0 0 0 0 0 0 0 0 0 0 0; 0128 31 677.871 221.574 300 -100 0.982 100 1 646 0 0 0 0 0 0 0 0 0 0 0 0; 0129 32 650 206.965 300 150 0.9841 100 1 725 0 0 0 0 0 0 0 0 0 0 0 0; 0130 33 632 108.293 250 0 0.9972 100 1 652 0 0 0 0 0 0 0 0 0 0 0 0; 0131 34 508 166.688 167 0 1.0123 100 1 508 0 0 0 0 0 0 0 0 0 0 0 0; 0132 35 650 210.661 300 -100 1.0494 100 1 687 0 0 0 0 0 0 0 0 0 0 0 0; 0133 36 560 100.165 240 0 1.0636 100 1 580 0 0 0 0 0 0 0 0 0 0 0 0; 0134 37 540 -1.36945 250 0 1.0275 100 1 564 0 0 0 0 0 0 0 0 0 0 0 0; 0135 38 830 21.7327 300 -150 1.0265 100 1 865 0 0 0 0 0 0 0 0 0 0 0 0; 0136 39 1000 78.4674 300 -100 1.03 100 1 1100 0 0 0 0 0 0 0 0 0 0 0 0; 0137 ]; 0138 0139 %% branch data 0140 % fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax 0141 mpc.branch = [ 0142 1 2 0.0035 0.0411 0.6987 600 600 600 0 0 1 -360 360; 0143 1 39 0.001 0.025 0.75 1000 1000 1000 0 0 1 -360 360; 0144 2 3 0.0013 0.0151 0.2572 500 500 500 0 0 1 -360 360; 0145 2 25 0.007 0.0086 0.146 500 500 500 0 0 1 -360 360; 0146 2 30 0 0.0181 0 900 900 2500 1.025 0 1 -360 360; 0147 3 4 0.0013 0.0213 0.2214 500 500 500 0 0 1 -360 360; 0148 3 18 0.0011 0.0133 0.2138 500 500 500 0 0 1 -360 360; 0149 4 5 0.0008 0.0128 0.1342 600 600 600 0 0 1 -360 360; 0150 4 14 0.0008 0.0129 0.1382 500 500 500 0 0 1 -360 360; 0151 5 6 0.0002 0.0026 0.0434 1200 1200 1200 0 0 1 -360 360; 0152 5 8 0.0008 0.0112 0.1476 900 900 900 0 0 1 -360 360; 0153 6 7 0.0006 0.0092 0.113 900 900 900 0 0 1 -360 360; 0154 6 11 0.0007 0.0082 0.1389 480 480 480 0 0 1 -360 360; 0155 6 31 0 0.025 0 1800 1800 1800 1.07 0 1 -360 360; 0156 7 8 0.0004 0.0046 0.078 900 900 900 0 0 1 -360 360; 0157 8 9 0.0023 0.0363 0.3804 900 900 900 0 0 1 -360 360; 0158 9 39 0.001 0.025 1.2 900 900 900 0 0 1 -360 360; 0159 10 11 0.0004 0.0043 0.0729 600 600 600 0 0 1 -360 360; 0160 10 13 0.0004 0.0043 0.0729 600 600 600 0 0 1 -360 360; 0161 10 32 0 0.02 0 900 900 2500 1.07 0 1 -360 360; 0162 12 11 0.0016 0.0435 0 500 500 500 1.006 0 1 -360 360; 0163 12 13 0.0016 0.0435 0 500 500 500 1.006 0 1 -360 360; 0164 13 14 0.0009 0.0101 0.1723 600 600 600 0 0 1 -360 360; 0165 14 15 0.0018 0.0217 0.366 600 600 600 0 0 1 -360 360; 0166 15 16 0.0009 0.0094 0.171 600 600 600 0 0 1 -360 360; 0167 16 17 0.0007 0.0089 0.1342 600 600 600 0 0 1 -360 360; 0168 16 19 0.0016 0.0195 0.304 600 600 2500 0 0 1 -360 360; 0169 16 21 0.0008 0.0135 0.2548 600 600 600 0 0 1 -360 360; 0170 16 24 0.0003 0.0059 0.068 600 600 600 0 0 1 -360 360; 0171 17 18 0.0007 0.0082 0.1319 600 600 600 0 0 1 -360 360; 0172 17 27 0.0013 0.0173 0.3216 600 600 600 0 0 1 -360 360; 0173 19 20 0.0007 0.0138 0 900 900 2500 1.06 0 1 -360 360; 0174 19 33 0.0007 0.0142 0 900 900 2500 1.07 0 1 -360 360; 0175 20 34 0.0009 0.018 0 900 900 2500 1.009 0 1 -360 360; 0176 21 22 0.0008 0.014 0.2565 900 900 900 0 0 1 -360 360; 0177 22 23 0.0006 0.0096 0.1846 600 600 600 0 0 1 -360 360; 0178 22 35 0 0.0143 0 900 900 2500 1.025 0 1 -360 360; 0179 23 24 0.0022 0.035 0.361 600 600 600 0 0 1 -360 360; 0180 23 36 0.0005 0.0272 0 900 900 2500 1 0 1 -360 360; 0181 25 26 0.0032 0.0323 0.531 600 600 600 0 0 1 -360 360; 0182 25 37 0.0006 0.0232 0 900 900 2500 1.025 0 1 -360 360; 0183 26 27 0.0014 0.0147 0.2396 600 600 600 0 0 1 -360 360; 0184 26 28 0.0043 0.0474 0.7802 600 600 600 0 0 1 -360 360; 0185 26 29 0.0057 0.0625 1.029 600 600 600 0 0 1 -360 360; 0186 28 29 0.0014 0.0151 0.249 600 600 600 0 0 1 -360 360; 0187 29 38 0.0008 0.0156 0 1200 1200 2500 1.025 0 1 -360 360; 0188 ]; 0189 0190 %%----- OPF Data -----%% 0191 %% generator cost data 0192 % 1 startup shutdown n x1 y1 ... xn yn 0193 % 2 startup shutdown n c(n-1) ... c0 0194 mpc.gencost = [ 0195 2 0 0 3 0.01 0.3 0.2; 0196 2 0 0 3 0.01 0.3 0.2; 0197 2 0 0 3 0.01 0.3 0.2; 0198 2 0 0 3 0.01 0.3 0.2; 0199 2 0 0 3 0.01 0.3 0.2; 0200 2 0 0 3 0.01 0.3 0.2; 0201 2 0 0 3 0.01 0.3 0.2; 0202 2 0 0 3 0.01 0.3 0.2; 0203 2 0 0 3 0.01 0.3 0.2; 0204 2 0 0 3 0.01 0.3 0.2; 0205 ];