Description of case33bw

0001 function mpc = case33bw
0002 %CASE33BW  Power flow data for 33 bus distribution system from Baran & Wu
0003 %    Please see CASEFORMAT for details on the case file format.
0004 %
0005 %    Data from ...
0006 %       M. E. Baran and F. F. Wu, "Network reconfiguration in distribution
0007 %       systems for loss reduction and load balancing," in IEEE Transactions
0008 %       on Power Delivery, vol. 4, no. 2, pp. 1401-1407, Apr 1989.
0009 %       doi: 10.1109/61.25627
0010 %       URL: http://doi.org/10.1109/61.25627
0011 
0012 %% MATPOWER Case Format : Version 2
0013 mpc.version = '2';
0014 
0015 %%-----  Power Flow Data  -----%%
0016 %% system MVA base
0017 mpc.baseMVA = 100;
0018 
0019 %% bus data
0020 %    bus_i    type    Pd    Qd    Gs    Bs    area    Vm    Va    baseKV    zone    Vmax    Vmin
0021 mpc.bus = [  %% (Pd and Qd are specified in kW & kVAr here, converted to MW & MVAr below)
0022     1    3    0    0    0    0    1    1    0    12.66    1    1    1;
0023     2    1    100    60    0    0    1    1    0    12.66    1    1.1    0.9;
0024     3    1    90    40    0    0    1    1    0    12.66    1    1.1    0.9;
0025     4    1    120    80    0    0    1    1    0    12.66    1    1.1    0.9;
0026     5    1    60    30    0    0    1    1    0    12.66    1    1.1    0.9;
0027     6    1    60    20    0    0    1    1    0    12.66    1    1.1    0.9;
0028     7    1    200    100    0    0    1    1    0    12.66    1    1.1    0.9;
0029     8    1    200    100    0    0    1    1    0    12.66    1    1.1    0.9;
0030     9    1    60    20    0    0    1    1    0    12.66    1    1.1    0.9;
0031     10    1    60    20    0    0    1    1    0    12.66    1    1.1    0.9;
0032     11    1    45    30    0    0    1    1    0    12.66    1    1.1    0.9;
0033     12    1    60    35    0    0    1    1    0    12.66    1    1.1    0.9;
0034     13    1    60    35    0    0    1    1    0    12.66    1    1.1    0.9;
0035     14    1    120    80    0    0    1    1    0    12.66    1    1.1    0.9;
0036     15    1    60    10    0    0    1    1    0    12.66    1    1.1    0.9;
0037     16    1    60    20    0    0    1    1    0    12.66    1    1.1    0.9;
0038     17    1    60    20    0    0    1    1    0    12.66    1    1.1    0.9;
0039     18    1    90    40    0    0    1    1    0    12.66    1    1.1    0.9;
0040     19    1    90    40    0    0    1    1    0    12.66    1    1.1    0.9;
0041     20    1    90    40    0    0    1    1    0    12.66    1    1.1    0.9;
0042     21    1    90    40    0    0    1    1    0    12.66    1    1.1    0.9;
0043     22    1    90    40    0    0    1    1    0    12.66    1    1.1    0.9;
0044     23    1    90    50    0    0    1    1    0    12.66    1    1.1    0.9;
0045     24    1    420    200    0    0    1    1    0    12.66    1    1.1    0.9;
0046     25    1    420    200    0    0    1    1    0    12.66    1    1.1    0.9;
0047     26    1    60    25    0    0    1    1    0    12.66    1    1.1    0.9;
0048     27    1    60    25    0    0    1    1    0    12.66    1    1.1    0.9;
0049     28    1    60    20    0    0    1    1    0    12.66    1    1.1    0.9;
0050     29    1    120    70    0    0    1    1    0    12.66    1    1.1    0.9;
0051     30    1    200    600    0    0    1    1    0    12.66    1    1.1    0.9;
0052     31    1    150    70    0    0    1    1    0    12.66    1    1.1    0.9;
0053     32    1    210    100    0    0    1    1    0    12.66    1    1.1    0.9;
0054     33    1    60    40    0    0    1    1    0    12.66    1    1.1    0.9;
0055 ];
0056 
0057 %% generator data
0058 %    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
0059 mpc.gen = [
0060     1    0    0    10    -10    1    100    1    10    0    0    0    0    0    0    0    0    0    0    0    0;
0061 ];
0062 
0063 %% branch data
0064 %    fbus    tbus    r    x    b    rateA    rateB    rateC    ratio    angle    status    angmin    angmax
0065 mpc.branch = [  %% (r and x specified in ohms here, converted to p.u. below)
0066     1    2    0.0922    0.0470    0    0    0    0    0    0    1    -360    360;
0067     2    3    0.4930    0.2511    0    0    0    0    0    0    1    -360    360;
0068     3    4    0.3660    0.1864    0    0    0    0    0    0    1    -360    360;
0069     4    5    0.3811    0.1941    0    0    0    0    0    0    1    -360    360;
0070     5    6    0.8190    0.7070    0    0    0    0    0    0    1    -360    360;
0071     6    7    0.1872    0.6188    0    0    0    0    0    0    1    -360    360;
0072     7    8    0.7114    0.2351    0    0    0    0    0    0    1    -360    360;
0073     8    9    1.0300    0.7400    0    0    0    0    0    0    1    -360    360;
0074     9    10    1.0440    0.7400    0    0    0    0    0    0    1    -360    360;
0075     10    11    0.1966    0.0650    0    0    0    0    0    0    1    -360    360;
0076     11    12    0.3744    0.1238    0    0    0    0    0    0    1    -360    360;
0077     12    13    1.4680    1.1550    0    0    0    0    0    0    1    -360    360;
0078     13    14    0.5416    0.7129    0    0    0    0    0    0    1    -360    360;
0079     14    15    0.5910    0.5260    0    0    0    0    0    0    1    -360    360;
0080     15    16    0.7463    0.5450    0    0    0    0    0    0    1    -360    360;
0081     16    17    1.2890    1.7210    0    0    0    0    0    0    1    -360    360;
0082     17    18    0.7320    0.5740    0    0    0    0    0    0    1    -360    360;
0083     2    19    0.1640    0.1565    0    0    0    0    0    0    1    -360    360;
0084     19    20    1.5042    1.3554    0    0    0    0    0    0    1    -360    360;
0085     20    21    0.4095    0.4784    0    0    0    0    0    0    1    -360    360;
0086     21    22    0.7089    0.9373    0    0    0    0    0    0    1    -360    360;
0087     3    23    0.4512    0.3083    0    0    0    0    0    0    1    -360    360;
0088     23    24    0.8980    0.7091    0    0    0    0    0    0    1    -360    360;
0089     24    25    0.8960    0.7011    0    0    0    0    0    0    1    -360    360;
0090     6    26    0.2030    0.1034    0    0    0    0    0    0    1    -360    360;
0091     26    27    0.2842    0.1447    0    0    0    0    0    0    1    -360    360;
0092     27    28    1.0590    0.9337    0    0    0    0    0    0    1    -360    360;
0093     28    29    0.8042    0.7006    0    0    0    0    0    0    1    -360    360;
0094     29    30    0.5075    0.2585    0    0    0    0    0    0    1    -360    360;
0095     30    31    0.9744    0.9630    0    0    0    0    0    0    1    -360    360;
0096     31    32    0.3105    0.3619    0    0    0    0    0    0    1    -360    360;
0097     32    33    0.3410    0.5302    0    0    0    0    0    0    1    -360    360;
0098     21    8    2.0000    2.0000    0    0    0    0    0    0    0    -360    360;
0099     9    15    2.0000    2.0000    0    0    0    0    0    0    0    -360    360;
0100     12    22    2.0000    2.0000    0    0    0    0    0    0    0    -360    360;
0101     18    33    0.5000    0.5000    0    0    0    0    0    0    0    -360    360;
0102     25    29    0.5000    0.5000    0    0    0    0    0    0    0    -360    360;
0103 ];
0104 
0105 %%-----  OPF Data  -----%%
0106 %% generator cost data
0107 %    1    startup    shutdown    n    x1    y1    ...    xn    yn
0108 %    2    startup    shutdown    n    c(n-1)    ...    c0
0109 mpc.gencost = [
0110     2    0    0    3    0    20    0;
0111 ];
0112 
0113 
0114 %% convert branch impedances from Ohms to p.u.
0115 [PQ, PV, REF, NONE, BUS_I, BUS_TYPE, PD, QD, GS, BS, BUS_AREA, VM, ...
0116     VA, BASE_KV, ZONE, VMAX, VMIN, LAM_P, LAM_Q, MU_VMAX, MU_VMIN] = idx_bus;
0117 [F_BUS, T_BUS, BR_R, BR_X, BR_B, RATE_A, RATE_B, RATE_C, ...
0118     TAP, SHIFT, BR_STATUS, PF, QF, PT, QT, MU_SF, MU_ST, ...
0119     ANGMIN, ANGMAX, MU_ANGMIN, MU_ANGMAX] = idx_brch;
0120 Vbase = mpc.bus(1, BASE_KV) * 1e3;      %% in Volts
0121 Sbase = mpc.baseMVA * 1e6;              %% in VA
0122 mpc.branch(:, [BR_R BR_X]) = mpc.branch(:, [BR_R BR_X]) / (Vbase^2 / Sbase);
0123 
0124 %% convert loads from kW to MW
0125 mpc.bus(:, [PD, QD]) = mpc.bus(:, [PD, QD]) / 1e3;
case33bw

PURPOSE

SYNOPSIS

DESCRIPTION

CROSS-REFERENCE INFORMATION

SOURCE CODE
