Source Code for Module pypower.case14

  1  # Copyright (C) 1996-2011 Power System Engineering Research Center (PSERC) 
  2  # Copyright (C) 2011 Richard Lincoln 
  3  # 
  4  # PYPOWER is free software: you can redistribute it and/or modify 
  5  # it under the terms of the GNU General Public License as published 
  6  # by the Free Software Foundation, either version 3 of the License, 
  7  # or (at your option) any later version. 
  8  # 
  9  # PYPOWER is distributed in the hope that it will be useful, 
 10  # but WITHOUT ANY WARRANTY; without even the implied warranty of 
 11  # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 
 12  # GNU General Public License for more details. 
 13  # 
 14  # You should have received a copy of the GNU General Public License 
 15  # along with PYPOWER. If not, see <http://www.gnu.org/licenses/>. 
 16   
 17  """Power flow data for IEEE 14 bus test case. 
 18  """ 
 19   
 20  from numpy import array 
 21   
 22 -def case14(): 
 23      """Power flow data for IEEE 14 bus test case. 
 24      Please see L{caseformat} for details on the case file format. 
 25   
 26      This data was converted from IEEE Common Data Format 
 27      (ieee14cdf.txt) on 20-Sep-2004 by cdf2matp, rev. 1.11 
 28   
 29      Converted from IEEE CDF file from: 
 30      U{http://www.ee.washington.edu/research/pstca/} 
 31   
 32      08/19/93 UW ARCHIVE           100.0  1962 W IEEE 14 Bus Test Case 
 33   
 34      @return: Power flow data for IEEE 14 bus test case. 
 35      """ 
 36      ppc = {"version": '2'} 
 37   
 38      ##-----  Power Flow Data  -----## 
 39      ## system MVA base 
 40      ppc["baseMVA"] = 100.0 
 41   
 42      ## bus data 
 43      # bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin 
 44      ppc["bus"] = array([ 
 45          [1,  3,  0,    0,   0, 0,  1, 1.06,    0,    0, 1, 1.06, 0.94], 
 46          [2,  2, 21.7, 12.7, 0, 0,  1, 1.045,  -4.98, 0, 1, 1.06, 0.94], 
 47          [3,  2, 94.2, 19,   0, 0,  1, 1.01,  -12.72, 0, 1, 1.06, 0.94], 
 48          [4,  1, 47.8, -3.9, 0, 0,  1, 1.019, -10.33, 0, 1, 1.06, 0.94], 
 49          [5,  1,  7.6,  1.6, 0, 0,  1, 1.02,   -8.78, 0, 1, 1.06, 0.94], 
 50          [6,  2, 11.2,  7.5, 0, 0,  1, 1.07,  -14.22, 0, 1, 1.06, 0.94], 
 51          [7,  1,  0,    0,   0, 0,  1, 1.062, -13.37, 0, 1, 1.06, 0.94], 
 52          [8,  2,  0,    0,   0, 0,  1, 1.09,  -13.36, 0, 1, 1.06, 0.94], 
 53          [9,  1, 29.5, 16.6, 0, 19, 1, 1.056, -14.94, 0, 1, 1.06, 0.94], 
 54          [10, 1,  9,    5.8, 0, 0,  1, 1.051, -15.1,  0, 1, 1.06, 0.94], 
 55          [11, 1,  3.5,  1.8, 0, 0,  1, 1.057, -14.79, 0, 1, 1.06, 0.94], 
 56          [12, 1,  6.1,  1.6, 0, 0,  1, 1.055, -15.07, 0, 1, 1.06, 0.94], 
 57          [13, 1, 13.5,  5.8, 0, 0,  1, 1.05,  -15.16, 0, 1, 1.06, 0.94], 
 58          [14, 1, 14.9,  5,   0, 0,  1, 1.036, -16.04, 0, 1, 1.06, 0.94] 
 59      ]) 
 60   
 61      ## generator data 
 62      # bus, Pg, Qg, Qmax, Qmin, Vg, mBase, status, Pmax, Pmin, Pc1, Pc2, 
 63      # Qc1min, Qc1max, Qc2min, Qc2max, ramp_agc, ramp_10, ramp_30, ramp_q, apf 
 64      ppc["gen"] = array([ 
 65          [1, 232.4, -16.9, 10,   0, 1.06,  100, 1, 332.4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], 
 66          [2,  40,    42.4, 50, -40, 1.045, 100, 1, 140,   0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], 
 67          [3,   0,    23.4, 40,   0, 1.01,  100, 1, 100,   0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], 
 68          [6,   0,    12.2, 24,  -6, 1.07,  100, 1, 100,   0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], 
 69          [8,   0,    17.4, 24,  -6, 1.09,  100, 1, 100,   0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0] 
 70      ]) 
 71   
 72      ## branch data 
 73      # fbus, tbus, r, x, b, rateA, rateB, rateC, ratio, angle, status, angmin, angmax 
 74      ppc["branch"] = array([ 
 75          [1,   2, 0.01938, 0.05917, 0.0528, 9900, 0, 0, 0,     0, 1, -360, 360], 
 76          [1,   5, 0.05403, 0.22304, 0.0492, 9900, 0, 0, 0,     0, 1, -360, 360], 
 77          [2,   3, 0.04699, 0.19797, 0.0438, 9900, 0, 0, 0,     0, 1, -360, 360], 
 78          [2,   4, 0.05811, 0.17632, 0.034,  9900, 0, 0, 0,     0, 1, -360, 360], 
 79          [2,   5, 0.05695, 0.17388, 0.0346, 9900, 0, 0, 0,     0, 1, -360, 360], 
 80          [3,   4, 0.06701, 0.17103, 0.0128, 9900, 0, 0, 0,     0, 1, -360, 360], 
 81          [4,   5, 0.01335, 0.04211, 0,      9900, 0, 0, 0,     0, 1, -360, 360], 
 82          [4,   7, 0,       0.20912, 0,      9900, 0, 0, 0.978, 0, 1, -360, 360], 
 83          [4,   9, 0,       0.55618, 0,      9900, 0, 0, 0.969, 0, 1, -360, 360], 
 84          [5,   6, 0,       0.25202, 0,      9900, 0, 0, 0.932, 0, 1, -360, 360], 
 85          [6,  11, 0.09498, 0.1989,  0,      9900, 0, 0, 0,     0, 1, -360, 360], 
 86          [6,  12, 0.12291, 0.25581, 0,      9900, 0, 0, 0,     0, 1, -360, 360], 
 87          [6,  13, 0.06615, 0.13027, 0,      9900, 0, 0, 0,     0, 1, -360, 360], 
 88          [7,   8, 0,       0.17615, 0,      9900, 0, 0, 0,     0, 1, -360, 360], 
 89          [7,   9, 0,       0.11001, 0,      9900, 0, 0, 0,     0, 1, -360, 360], 
 90          [9,  10, 0.03181, 0.0845,  0,      9900, 0, 0, 0,     0, 1, -360, 360], 
 91          [9,  14, 0.12711, 0.27038, 0,      9900, 0, 0, 0,     0, 1, -360, 360], 
 92          [10, 11, 0.08205, 0.19207, 0,      9900, 0, 0, 0,     0, 1, -360, 360], 
 93          [12, 13, 0.22092, 0.19988, 0,      9900, 0, 0, 0,     0, 1, -360, 360], 
 94          [13, 14, 0.17093, 0.34802, 0,      9900, 0, 0, 0,     0, 1, -360, 360] 
 95      ]) 
 96   
 97      ##-----  OPF Data  -----## 
 98      ## generator cost data 
 99      # 1 startup shutdown n x1 y1 ... xn yn 
100      # 2 startup shutdown n c(n-1) ... c0 
101      ppc["gencost"] = array([ 
102          [2, 0, 0, 3, 0.0430293, 20, 0], 
103          [2, 0, 0, 3, 0.25,      20, 0], 
104          [2, 0, 0, 3, 0.01,      40, 0], 
105          [2, 0, 0, 3, 0.01,      40, 0], 
106          [2, 0, 0, 3, 0.01,      40, 0] 
107      ]) 
108   
109      return ppc 
110