Source Code for Module pypower.case57

  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 57 bus test case. 
 18  """ 
 19   
 20  from numpy import array 
 21   
 22 -def case57(): 
 23      """Power flow data for IEEE 57 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      (ieee57cdf.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      Manually modified C{Qmax}, C{Qmin} on generator 1 to 200, -140, 
 33      respectively. 
 34   
 35      08/25/93 UW ARCHIVE           100.0  1961 W IEEE 57 Bus Test Case 
 36   
 37      @return: Power flow data for IEEE 57 bus test case. 
 38      """ 
 39      ppc = {"version": '2'} 
 40   
 41      ##-----  Power Flow Data  -----## 
 42      ## system MVA base 
 43      ppc["baseMVA"] = 100.0 
 44   
 45      ## bus data 
 46      # bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin 
 47      ppc["bus"] = array([ 
 48          [1, 3, 55, 17, 0, 0, 1, 1.04, 0, 0, 1, 1.06, 0.94], 
 49          [2, 2, 3, 88, 0, 0, 1, 1.01, -1.18, 0, 1, 1.06, 0.94], 
 50          [3, 2, 41, 21, 0, 0, 1, 0.985, -5.97, 0, 1, 1.06, 0.94], 
 51          [4, 1, 0, 0, 0, 0, 1, 0.981, -7.32, 0, 1, 1.06, 0.94], 
 52          [5, 1, 13, 4, 0, 0, 1, 0.976, -8.52, 0, 1, 1.06, 0.94], 
 53          [6, 2, 75, 2, 0, 0, 1, 0.98, -8.65, 0, 1, 1.06, 0.94], 
 54          [7, 1, 0, 0, 0, 0, 1, 0.984, -7.58, 0, 1, 1.06, 0.94], 
 55          [8, 2, 150, 22, 0, 0, 1, 1.005, -4.45, 0, 1, 1.06, 0.94], 
 56          [9, 2, 121, 26, 0, 0, 1, 0.98, -9.56, 0, 1, 1.06, 0.94], 
 57          [10, 1, 5, 2, 0, 0, 1, 0.986, -11.43, 0, 1, 1.06, 0.94], 
 58          [11, 1, 0, 0, 0, 0, 1, 0.974, -10.17, 0, 1, 1.06, 0.94], 
 59          [12, 2, 377, 24, 0, 0, 1, 1.015, -10.46, 0, 1, 1.06, 0.94], 
 60          [13, 1, 18, 2.3, 0, 0, 1, 0.979, -9.79, 0, 1, 1.06, 0.94], 
 61          [14, 1, 10.5, 5.3, 0, 0, 1, 0.97, -9.33, 0, 1, 1.06, 0.94], 
 62          [15, 1, 22, 5, 0, 0, 1, 0.988, -7.18, 0, 1, 1.06, 0.94], 
 63          [16, 1, 43, 3, 0, 0, 1, 1.013, -8.85, 0, 1, 1.06, 0.94], 
 64          [17, 1, 42, 8, 0, 0, 1, 1.017, -5.39, 0, 1, 1.06, 0.94], 
 65          [18, 1, 27.2, 9.8, 0, 10, 1, 1.001, -11.71, 0, 1, 1.06, 0.94], 
 66          [19, 1, 3.3, 0.6, 0, 0, 1, 0.97, -13.2, 0, 1, 1.06, 0.94], 
 67          [20, 1, 2.3, 1, 0, 0, 1, 0.964, -13.41, 0, 1, 1.06, 0.94], 
 68          [21, 1, 0, 0, 0, 0, 1, 1.008, -12.89, 0, 1, 1.06, 0.94], 
 69          [22, 1, 0, 0, 0, 0, 1, 1.01, -12.84, 0, 1, 1.06, 0.94], 
 70          [23, 1, 6.3, 2.1, 0, 0, 1, 1.008, -12.91, 0, 1, 1.06, 0.94], 
 71          [24, 1, 0, 0, 0, 0, 1, 0.999, -13.25, 0, 1, 1.06, 0.94], 
 72          [25, 1, 6.3, 3.2, 0, 5.9, 1, 0.982, -18.13, 0, 1, 1.06, 0.94], 
 73          [26, 1, 0, 0, 0, 0, 1, 0.959, -12.95, 0, 1, 1.06, 0.94], 
 74          [27, 1, 9.3, 0.5, 0, 0, 1, 0.982, -11.48, 0, 1, 1.06, 0.94], 
 75          [28, 1, 4.6, 2.3, 0, 0, 1, 0.997, -10.45, 0, 1, 1.06, 0.94], 
 76          [29, 1, 17, 2.6, 0, 0, 1, 1.01, -9.75, 0, 1, 1.06, 0.94], 
 77          [30, 1, 3.6, 1.8, 0, 0, 1, 0.962, -18.68, 0, 1, 1.06, 0.94], 
 78          [31, 1, 5.8, 2.9, 0, 0, 1, 0.936, -19.34, 0, 1, 1.06, 0.94], 
 79          [32, 1, 1.6, 0.8, 0, 0, 1, 0.949, -18.46, 0, 1, 1.06, 0.94], 
 80          [33, 1, 3.8, 1.9, 0, 0, 1, 0.947, -18.5, 0, 1, 1.06, 0.94], 
 81          [34, 1, 0, 0, 0, 0, 1, 0.959, -14.1, 0, 1, 1.06, 0.94], 
 82          [35, 1, 6, 3, 0, 0, 1, 0.966, -13.86, 0, 1, 1.06, 0.94], 
 83          [36, 1, 0, 0, 0, 0, 1, 0.976, -13.59, 0, 1, 1.06, 0.94], 
 84          [37, 1, 0, 0, 0, 0, 1, 0.985, -13.41, 0, 1, 1.06, 0.94], 
 85          [38, 1, 14, 7, 0, 0, 1, 1.013, -12.71, 0, 1, 1.06, 0.94], 
 86          [39, 1, 0, 0, 0, 0, 1, 0.983, -13.46, 0, 1, 1.06, 0.94], 
 87          [40, 1, 0, 0, 0, 0, 1, 0.973, -13.62, 0, 1, 1.06, 0.94], 
 88          [41, 1, 6.3, 3, 0, 0, 1, 0.996, -14.05, 0, 1, 1.06, 0.94], 
 89          [42, 1, 7.1, 4.4, 0, 0, 1, 0.966, -15.5, 0, 1, 1.06, 0.94], 
 90          [43, 1, 2, 1, 0, 0, 1, 1.01, -11.33, 0, 1, 1.06, 0.94], 
 91          [44, 1, 12, 1.8, 0, 0, 1, 1.017, -11.86, 0, 1, 1.06, 0.94], 
 92          [45, 1, 0, 0, 0, 0, 1, 1.036, -9.25, 0, 1, 1.06, 0.94], 
 93          [46, 1, 0, 0, 0, 0, 1, 1.05, -11.89, 0, 1, 1.06, 0.94], 
 94          [47, 1, 29.7, 11.6, 0, 0, 1, 1.033, -12.49, 0, 1, 1.06, 0.94], 
 95          [48, 1, 0, 0, 0, 0, 1, 1.027, -12.59, 0, 1, 1.06, 0.94], 
 96          [49, 1, 18, 8.5, 0, 0, 1, 1.036, -12.92, 0, 1, 1.06, 0.94], 
 97          [50, 1, 21, 10.5, 0, 0, 1, 1.023, -13.39, 0, 1, 1.06, 0.94], 
 98          [51, 1, 18, 5.3, 0, 0, 1, 1.052, -12.52, 0, 1, 1.06, 0.94], 
 99          [52, 1, 4.9, 2.2, 0, 0, 1, 0.98, -11.47, 0, 1, 1.06, 0.94], 
100          [53, 1, 20, 10, 0, 6.3, 1, 0.971, -12.23, 0, 1, 1.06, 0.94], 
101          [54, 1, 4.1, 1.4, 0, 0, 1, 0.996, -11.69, 0, 1, 1.06, 0.94], 
102          [55, 1, 6.8, 3.4, 0, 0, 1, 1.031, -10.78, 0, 1, 1.06, 0.94], 
103          [56, 1, 7.6, 2.2, 0, 0, 1, 0.968, -16.04, 0, 1, 1.06, 0.94], 
104          [57, 1, 6.7, 2, 0, 0, 1, 0.965, -16.56, 0, 1, 1.06, 0.94] 
105      ]) 
106   
107      ## generator data 
108      # bus, Pg, Qg, Qmax, Qmin, Vg, mBase, status, Pmax, Pmin, Pc1, Pc2, 
109      # Qc1min, Qc1max, Qc2min, Qc2max, ramp_agc, ramp_10, ramp_30, ramp_q, apf 
110      ppc["gen"] = array([ 
111          [1, 128.9, -16.1, 200, -140, 1.04, 100, 1, 575.88, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], 
112          [2, 0, -0.8, 50, -17, 1.01, 100, 1, 100, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], 
113          [3, 40, -1, 60, -10, 0.985, 100, 1, 140, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], 
114          [6, 0, 0.8, 25, -8, 0.98, 100, 1, 100, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], 
115          [8, 450, 62.1, 200, -140, 1.005, 100, 1, 550, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], 
116          [9, 0, 2.2, 9, -3, 0.98, 100, 1, 100, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], 
117          [12, 310, 128.5, 155, -150, 1.015, 100, 1, 410, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0] 
118      ]) 
119   
120      ## branch data 
121      # fbus, tbus, r, x, b, rateA, rateB, rateC, ratio, angle, status, angmin, angmax 
122      ppc["branch"] = array([ 
123          [1, 2, 0.0083, 0.028, 0.129, 9900, 0, 0, 0, 0, 1, -360, 360], 
124          [2, 3, 0.0298, 0.085, 0.0818, 9900, 0, 0, 0, 0, 1, -360, 360], 
125          [3, 4, 0.0112, 0.0366, 0.038, 9900, 0, 0, 0, 0, 1, -360, 360], 
126          [4, 5, 0.0625, 0.132, 0.0258, 9900, 0, 0, 0, 0, 1, -360, 360], 
127          [4, 6, 0.043, 0.148, 0.0348, 9900, 0, 0, 0, 0, 1, -360, 360], 
128          [6, 7, 0.02, 0.102, 0.0276, 9900, 0, 0, 0, 0, 1, -360, 360], 
129          [6, 8, 0.0339, 0.173, 0.047, 9900, 0, 0, 0, 0, 1, -360, 360], 
130          [8, 9, 0.0099, 0.0505, 0.0548, 9900, 0, 0, 0, 0, 1, -360, 360], 
131          [9, 10, 0.0369, 0.1679, 0.044, 9900, 0, 0, 0, 0, 1, -360, 360], 
132          [9, 11, 0.0258, 0.0848, 0.0218, 9900, 0, 0, 0, 0, 1, -360, 360], 
133          [9, 12, 0.0648, 0.295, 0.0772, 9900, 0, 0, 0, 0, 1, -360, 360], 
134          [9, 13, 0.0481, 0.158, 0.0406, 9900, 0, 0, 0, 0, 1, -360, 360], 
135          [13, 14, 0.0132, 0.0434, 0.011, 9900, 0, 0, 0, 0, 1, -360, 360], 
136          [13, 15, 0.0269, 0.0869, 0.023, 9900, 0, 0, 0, 0, 1, -360, 360], 
137          [1, 15, 0.0178, 0.091, 0.0988, 9900, 0, 0, 0, 0, 1, -360, 360], 
138          [1, 16, 0.0454, 0.206, 0.0546, 9900, 0, 0, 0, 0, 1, -360, 360], 
139          [1, 17, 0.0238, 0.108, 0.0286, 9900, 0, 0, 0, 0, 1, -360, 360], 
140          [3, 15, 0.0162, 0.053, 0.0544, 9900, 0, 0, 0, 0, 1, -360, 360], 
141          [4, 18, 0, 0.555, 0, 9900, 0, 0, 0.97, 0, 1, -360, 360], 
142          [4, 18, 0, 0.43, 0, 9900, 0, 0, 0.978, 0, 1, -360, 360], 
143          [5, 6, 0.0302, 0.0641, 0.0124, 9900, 0, 0, 0, 0, 1, -360, 360], 
144          [7, 8, 0.0139, 0.0712, 0.0194, 9900, 0, 0, 0, 0, 1, -360, 360], 
145          [10, 12, 0.0277, 0.1262, 0.0328, 9900, 0, 0, 0, 0, 1, -360, 360], 
146          [11, 13, 0.0223, 0.0732, 0.0188, 9900, 0, 0, 0, 0, 1, -360, 360], 
147          [12, 13, 0.0178, 0.058, 0.0604, 9900, 0, 0, 0, 0, 1, -360, 360], 
148          [12, 16, 0.018, 0.0813, 0.0216, 9900, 0, 0, 0, 0, 1, -360, 360], 
149          [12, 17, 0.0397, 0.179, 0.0476, 9900, 0, 0, 0, 0, 1, -360, 360], 
150          [14, 15, 0.0171, 0.0547, 0.0148, 9900, 0, 0, 0, 0, 1, -360, 360], 
151          [18, 19, 0.461, 0.685, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
152          [19, 20, 0.283, 0.434, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
153          [21, 20, 0, 0.7767, 0, 9900, 0, 0, 1.043, 0, 1, -360, 360], 
154          [21, 22, 0.0736, 0.117, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
155          [22, 23, 0.0099, 0.0152, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
156          [23, 24, 0.166, 0.256, 0.0084, 9900, 0, 0, 0, 0, 1, -360, 360], 
157          [24, 25, 0, 1.182, 0, 9900, 0, 0, 1, 0, 1, -360, 360], 
158          [24, 25, 0, 1.23, 0, 9900, 0, 0, 1, 0, 1, -360, 360], 
159          [24, 26, 0, 0.0473, 0, 9900, 0, 0, 1.043, 0, 1, -360, 360], 
160          [26, 27, 0.165, 0.254, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
161          [27, 28, 0.0618, 0.0954, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
162          [28, 29, 0.0418, 0.0587, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
163          [7, 29, 0, 0.0648, 0, 9900, 0, 0, 0.967, 0, 1, -360, 360], 
164          [25, 30, 0.135, 0.202, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
165          [30, 31, 0.326, 0.497, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
166          [31, 32, 0.507, 0.755, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
167          [32, 33, 0.0392, 0.036, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
168          [34, 32, 0, 0.953, 0, 9900, 0, 0, 0.975, 0, 1, -360, 360], 
169          [34, 35, 0.052, 0.078, 0.0032, 9900, 0, 0, 0, 0, 1, -360, 360], 
170          [35, 36, 0.043, 0.0537, 0.0016, 9900, 0, 0, 0, 0, 1, -360, 360], 
171          [36, 37, 0.029, 0.0366, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
172          [37, 38, 0.0651, 0.1009, 0.002, 9900, 0, 0, 0, 0, 1, -360, 360], 
173          [37, 39, 0.0239, 0.0379, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
174          [36, 40, 0.03, 0.0466, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
175          [22, 38, 0.0192, 0.0295, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
176          [11, 41, 0, 0.749, 0, 9900, 0, 0, 0.955, 0, 1, -360, 360], 
177          [41, 42, 0.207, 0.352, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
178          [41, 43, 0, 0.412, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
179          [38, 44, 0.0289, 0.0585, 0.002, 9900, 0, 0, 0, 0, 1, -360, 360], 
180          [15, 45, 0, 0.1042, 0, 9900, 0, 0, 0.955, 0, 1, -360, 360], 
181          [14, 46, 0, 0.0735, 0, 9900, 0, 0, 0.9, 0, 1, -360, 360], 
182          [46, 47, 0.023, 0.068, 0.0032, 9900, 0, 0, 0, 0, 1, -360, 360], 
183          [47, 48, 0.0182, 0.0233, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
184          [48, 49, 0.0834, 0.129, 0.0048, 9900, 0, 0, 0, 0, 1, -360, 360], 
185          [49, 50, 0.0801, 0.128, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
186          [50, 51, 0.1386, 0.22, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
187          [10, 51, 0, 0.0712, 0, 9900, 0, 0, 0.93, 0, 1, -360, 360], 
188          [13, 49, 0, 0.191, 0, 9900, 0, 0, 0.895, 0, 1, -360, 360], 
189          [29, 52, 0.1442, 0.187, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
190          [52, 53, 0.0762, 0.0984, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
191          [53, 54, 0.1878, 0.232, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
192          [54, 55, 0.1732, 0.2265, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
193          [11, 43, 0, 0.153, 0, 9900, 0, 0, 0.958, 0, 1, -360, 360], 
194          [44, 45, 0.0624, 0.1242, 0.004, 9900, 0, 0, 0, 0, 1, -360, 360], 
195          [40, 56, 0, 1.195, 0, 9900, 0, 0, 0.958, 0, 1, -360, 360], 
196          [56, 41, 0.553, 0.549, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
197          [56, 42, 0.2125, 0.354, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
198          [39, 57, 0, 1.355, 0, 9900, 0, 0, 0.98, 0, 1, -360, 360], 
199          [57, 56, 0.174, 0.26, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
200          [38, 49, 0.115, 0.177, 0.003, 9900, 0, 0, 0, 0, 1, -360, 360], 
201          [38, 48, 0.0312, 0.0482, 0, 9900, 0, 0, 0, 0, 1, -360, 360], 
202          [9, 55, 0, 0.1205, 0, 9900, 0, 0, 0.94, 0, 1, -360, 360] 
203      ]) 
204   
205      ##-----  OPF Data  -----## 
206      ## generator cost data 
207      # 1 startup shutdown n x1 y1 ... xn yn 
208      # 2 startup shutdown n c(n-1) ... c0 
209      ppc["gencost"] = array([ 
210          [2, 0, 0, 3, 0.0775795, 20, 0], 
211          [2, 0, 0, 3, 0.01, 40, 0], 
212          [2, 0, 0, 3, 0.25, 20, 0], 
213          [2, 0, 0, 3, 0.01, 40, 0], 
214          [2, 0, 0, 3, 0.0222222, 20, 0], 
215          [2, 0, 0, 3, 0.01, 40, 0], 
216          [2, 0, 0, 3, 0.0322581, 20, 0] 
217      ]) 
218   
219      return ppc 
220