Home > matpower4.1 > case39.m

case39

PURPOSE ^

CASE39 Power flow data for 39 bus New England system.

SYNOPSIS ^

function mpc = case39

DESCRIPTION ^

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.

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

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,v 1.14 2010/03/10 18:08:13 ray Exp $
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 ];

Generated on Mon 26-Jan-2015 15:00:13 by m2html © 2005