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 % $Id: case39.m 1559 2010-03-10 18:08:32Z ray $ 0073 0074 %% MATPOWER Case Format : Version 2 0075 mpc.version = '2'; 0076 0077 %%----- Power Flow Data -----%% 0078 %% system MVA base 0079 mpc.baseMVA = 100; 0080 0081 %% bus data 0082 % bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin 0083 mpc.bus = [ 0084 1 1 97.6 44.2 0 0 2 1.0393836 -13.536602 345 1 1.06 0.94; 0085 2 1 0 0 0 0 2 1.0484941 -9.7852666 345 1 1.06 0.94; 0086 3 1 322 2.4 0 0 2 1.0307077 -12.276384 345 1 1.06 0.94; 0087 4 1 500 184 0 0 1 1.00446 -12.626734 345 1 1.06 0.94; 0088 5 1 0 0 0 0 1 1.0060063 -11.192339 345 1 1.06 0.94; 0089 6 1 0 0 0 0 1 1.0082256 -10.40833 345 1 1.06 0.94; 0090 7 1 233.8 84 0 0 1 0.99839728 -12.755626 345 1 1.06 0.94; 0091 8 1 522 176.6 0 0 1 0.99787232 -13.335844 345 1 1.06 0.94; 0092 9 1 6.5 -66.6 0 0 1 1.038332 -14.178442 345 1 1.06 0.94; 0093 10 1 0 0 0 0 1 1.0178431 -8.170875 345 1 1.06 0.94; 0094 11 1 0 0 0 0 1 1.0133858 -8.9369663 345 1 1.06 0.94; 0095 12 1 8.53 88 0 0 1 1.000815 -8.9988236 345 1 1.06 0.94; 0096 13 1 0 0 0 0 1 1.014923 -8.9299272 345 1 1.06 0.94; 0097 14 1 0 0 0 0 1 1.012319 -10.715295 345 1 1.06 0.94; 0098 15 1 320 153 0 0 3 1.0161854 -11.345399 345 1 1.06 0.94; 0099 16 1 329 32.3 0 0 3 1.0325203 -10.033348 345 1 1.06 0.94; 0100 17 1 0 0 0 0 2 1.0342365 -11.116436 345 1 1.06 0.94; 0101 18 1 158 30 0 0 2 1.0315726 -11.986168 345 1 1.06 0.94; 0102 19 1 0 0 0 0 3 1.0501068 -5.4100729 345 1 1.06 0.94; 0103 20 1 680 103 0 0 3 0.99101054 -6.8211783 345 1 1.06 0.94; 0104 21 1 274 115 0 0 3 1.0323192 -7.6287461 345 1 1.06 0.94; 0105 22 1 0 0 0 0 3 1.0501427 -3.1831199 345 1 1.06 0.94; 0106 23 1 247.5 84.6 0 0 3 1.0451451 -3.3812763 345 1 1.06 0.94; 0107 24 1 308.6 -92.2 0 0 3 1.038001 -9.9137585 345 1 1.06 0.94; 0108 25 1 224 47.2 0 0 2 1.0576827 -8.3692354 345 1 1.06 0.94; 0109 26 1 139 17 0 0 2 1.0525613 -9.4387696 345 1 1.06 0.94; 0110 27 1 281 75.5 0 0 2 1.0383449 -11.362152 345 1 1.06 0.94; 0111 28 1 206 27.6 0 0 3 1.0503737 -5.9283592 345 1 1.06 0.94; 0112 29 1 283.5 26.9 0 0 3 1.0501149 -3.1698741 345 1 1.06 0.94; 0113 30 2 0 0 0 0 2 1.0499 -7.3704746 345 1 1.06 0.94; 0114 31 3 9.2 4.6 0 0 1 0.982 0 345 1 1.06 0.94; 0115 32 2 0 0 0 0 1 0.9841 -0.1884374 345 1 1.06 0.94; 0116 33 2 0 0 0 0 3 0.9972 -0.19317445 345 1 1.06 0.94; 0117 34 2 0 0 0 0 3 1.0123 -1.631119 345 1 1.06 0.94; 0118 35 2 0 0 0 0 3 1.0494 1.7765069 345 1 1.06 0.94; 0119 36 2 0 0 0 0 3 1.0636 4.4684374 345 1 1.06 0.94; 0120 37 2 0 0 0 0 2 1.0275 -1.5828988 345 1 1.06 0.94; 0121 38 2 0 0 0 0 3 1.0265 3.8928177 345 1 1.06 0.94; 0122 39 2 1104 250 0 0 1 1.03 -14.535256 345 1 1.06 0.94; 0123 ]; 0124 0125 %% generator data 0126 % 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 0127 mpc.gen = [ 0128 30 250 161.762 400 140 1.0499 100 1 1040 0 0 0 0 0 0 0 0 0 0 0 0; 0129 31 677.871 221.574 300 -100 0.982 100 1 646 0 0 0 0 0 0 0 0 0 0 0 0; 0130 32 650 206.965 300 150 0.9841 100 1 725 0 0 0 0 0 0 0 0 0 0 0 0; 0131 33 632 108.293 250 0 0.9972 100 1 652 0 0 0 0 0 0 0 0 0 0 0 0; 0132 34 508 166.688 167 0 1.0123 100 1 508 0 0 0 0 0 0 0 0 0 0 0 0; 0133 35 650 210.661 300 -100 1.0494 100 1 687 0 0 0 0 0 0 0 0 0 0 0 0; 0134 36 560 100.165 240 0 1.0636 100 1 580 0 0 0 0 0 0 0 0 0 0 0 0; 0135 37 540 -1.36945 250 0 1.0275 100 1 564 0 0 0 0 0 0 0 0 0 0 0 0; 0136 38 830 21.7327 300 -150 1.0265 100 1 865 0 0 0 0 0 0 0 0 0 0 0 0; 0137 39 1000 78.4674 300 -100 1.03 100 1 1100 0 0 0 0 0 0 0 0 0 0 0 0; 0138 ]; 0139 0140 %% branch data 0141 % fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax 0142 mpc.branch = [ 0143 1 2 0.0035 0.0411 0.6987 600 600 600 0 0 1 -360 360; 0144 1 39 0.001 0.025 0.75 1000 1000 1000 0 0 1 -360 360; 0145 2 3 0.0013 0.0151 0.2572 500 500 500 0 0 1 -360 360; 0146 2 25 0.007 0.0086 0.146 500 500 500 0 0 1 -360 360; 0147 2 30 0 0.0181 0 900 900 2500 1.025 0 1 -360 360; 0148 3 4 0.0013 0.0213 0.2214 500 500 500 0 0 1 -360 360; 0149 3 18 0.0011 0.0133 0.2138 500 500 500 0 0 1 -360 360; 0150 4 5 0.0008 0.0128 0.1342 600 600 600 0 0 1 -360 360; 0151 4 14 0.0008 0.0129 0.1382 500 500 500 0 0 1 -360 360; 0152 5 6 0.0002 0.0026 0.0434 1200 1200 1200 0 0 1 -360 360; 0153 5 8 0.0008 0.0112 0.1476 900 900 900 0 0 1 -360 360; 0154 6 7 0.0006 0.0092 0.113 900 900 900 0 0 1 -360 360; 0155 6 11 0.0007 0.0082 0.1389 480 480 480 0 0 1 -360 360; 0156 6 31 0 0.025 0 1800 1800 1800 1.07 0 1 -360 360; 0157 7 8 0.0004 0.0046 0.078 900 900 900 0 0 1 -360 360; 0158 8 9 0.0023 0.0363 0.3804 900 900 900 0 0 1 -360 360; 0159 9 39 0.001 0.025 1.2 900 900 900 0 0 1 -360 360; 0160 10 11 0.0004 0.0043 0.0729 600 600 600 0 0 1 -360 360; 0161 10 13 0.0004 0.0043 0.0729 600 600 600 0 0 1 -360 360; 0162 10 32 0 0.02 0 900 900 2500 1.07 0 1 -360 360; 0163 12 11 0.0016 0.0435 0 500 500 500 1.006 0 1 -360 360; 0164 12 13 0.0016 0.0435 0 500 500 500 1.006 0 1 -360 360; 0165 13 14 0.0009 0.0101 0.1723 600 600 600 0 0 1 -360 360; 0166 14 15 0.0018 0.0217 0.366 600 600 600 0 0 1 -360 360; 0167 15 16 0.0009 0.0094 0.171 600 600 600 0 0 1 -360 360; 0168 16 17 0.0007 0.0089 0.1342 600 600 600 0 0 1 -360 360; 0169 16 19 0.0016 0.0195 0.304 600 600 2500 0 0 1 -360 360; 0170 16 21 0.0008 0.0135 0.2548 600 600 600 0 0 1 -360 360; 0171 16 24 0.0003 0.0059 0.068 600 600 600 0 0 1 -360 360; 0172 17 18 0.0007 0.0082 0.1319 600 600 600 0 0 1 -360 360; 0173 17 27 0.0013 0.0173 0.3216 600 600 600 0 0 1 -360 360; 0174 19 20 0.0007 0.0138 0 900 900 2500 1.06 0 1 -360 360; 0175 19 33 0.0007 0.0142 0 900 900 2500 1.07 0 1 -360 360; 0176 20 34 0.0009 0.018 0 900 900 2500 1.009 0 1 -360 360; 0177 21 22 0.0008 0.014 0.2565 900 900 900 0 0 1 -360 360; 0178 22 23 0.0006 0.0096 0.1846 600 600 600 0 0 1 -360 360; 0179 22 35 0 0.0143 0 900 900 2500 1.025 0 1 -360 360; 0180 23 24 0.0022 0.035 0.361 600 600 600 0 0 1 -360 360; 0181 23 36 0.0005 0.0272 0 900 900 2500 1 0 1 -360 360; 0182 25 26 0.0032 0.0323 0.531 600 600 600 0 0 1 -360 360; 0183 25 37 0.0006 0.0232 0 900 900 2500 1.025 0 1 -360 360; 0184 26 27 0.0014 0.0147 0.2396 600 600 600 0 0 1 -360 360; 0185 26 28 0.0043 0.0474 0.7802 600 600 600 0 0 1 -360 360; 0186 26 29 0.0057 0.0625 1.029 600 600 600 0 0 1 -360 360; 0187 28 29 0.0014 0.0151 0.249 600 600 600 0 0 1 -360 360; 0188 29 38 0.0008 0.0156 0 1200 1200 2500 1.025 0 1 -360 360; 0189 ]; 0190 0191 %%----- OPF Data -----%% 0192 %% generator cost data 0193 % 1 startup shutdown n x1 y1 ... xn yn 0194 % 2 startup shutdown n c(n-1) ... c0 0195 mpc.gencost = [ 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 2 0 0 3 0.01 0.3 0.2; 0206 ];