Source Code for Module pypower.t.t_auction_case

  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 testing auction code. 
 18  """ 
 19   
 20  from numpy import array 
 21   
 22   
 23 -def t_auction_case(): 
 24      """Power flow data for testing auction code. 
 25      """ 
 26      ##-----  Power Flow Data  -----## 
 27      ## system MVA base 
 28      baseMVA = 100.0 
 29   
 30      ## bus data 
 31      # bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin 
 32      bus = array([ 
 33                  [1, 3, 0, 0, 0, 0, 1, 1, 0, 135, 1, 1.05, 0.95], 
 34                  [2, 2, 5.04485072, 2.95251632, 0, 0, 1, 1, 0, 135, 1, 1.1, 0.95], 
 35                  [3, 1, 22.3182336, 11.1591168, 0, 0, 1, 1, 0, 135, 1, 1.05, 0.95], 
 36                  [4, 1, 8.8343008, 1.8598528, 0, 0, 1, 1, 0, 135, 1, 1.05, 0.95], 
 37                  [5, 1, 0, 0, 0, 0.19, 1, 1, 0, 135, 1, 1.05, 0.95], 
 38                  [6, 1, 0, 0, 0, 0, 1, 1, 0, 135, 1, 1.05, 0.95], 
 39                  [7, 2, 0, 0, 0, 0, 1, 1, 0, 135, 1, 1.05, 0.95], 
 40                  [8, 1, 13.948896, 13.948896, 0, 0, 1, 1, 0, 135, 1, 1.05, 0.95], 
 41                  [9, 1, 0, 0, 0, 0, 1, 1, 0, 135, 1, 1.05, 0.95], 
 42                  [10, 1, 6.7419664, 2.324816, 0, 0, 3, 1, 0, 135, 1, 1.05, 0.95], 
 43                  [11, 1, 0, 0, 0, 0, 1, 1, 0, 135, 1, 1.05, 0.95], 
 44                  [12, 1, 13.0189696, 8.71806, 0, 0, 2, 1, 0, 135, 1, 1.05, 0.95], 
 45                  [13, 2, 0, 0, 0, 0, 2, 1, 0, 135, 1, 1.1, 0.95], 
 46                  [14, 1, 7.2069296, 1.8598528, 0, 0, 2, 1, 0, 135, 1, 1.05, 0.95], 
 47                  [15, 2, 0, 0, 0, 0, 2, 1, 0, 135, 1, 1.05, 0.95], 
 48                  [16, 1, 4.068428, 2.0923344, 0, 0, 2, 1, 0, 135, 1, 1.05, 0.95], 
 49                  [17, 1, 10.461672, 6.7419664, 0, 0, 2, 1, 0, 135, 1, 1.05, 0.95], 
 50                  [18, 1, 3.7197056, 1.0461672, 0, 0, 2, 1, 0, 135, 1, 1.05, 0.95], 
 51                  [19, 1, 11.042876, 3.9521872, 0, 0, 2, 1, 0, 135, 1, 1.05, 0.95], 
 52                  [20, 1, 2.5572976, 0.8136856, 0, 0, 2, 1, 0, 135, 1, 1.05, 0.95], 
 53                  [21, 1, 3.39036344738, 2.1698326063232, 0, 0, 3, 1, 0, 135, 1, 1.05, 0.95], 
 54                  [220, 2, 0, 0, 0, 0, 3, 1, 0, 135, 1, 1.1, 0.95], 
 55                  [23, 2, 22.3182336, 11.1591168, 0, 0, 2, 1, 0, 135, 1, 1.1, 0.95], 
 56                  [24, 1, 10.1129496, 7.7881336, 0, 0.04, 3, 1, 0, 135, 1, 1.05, 0.95], 
 57                  [25, 1, 0, 0, 0, 0, 3, 1, 0, 135, 1, 1.05, 0.95], 
 58                  [26, 1, 4.068428, 2.6735384, 0, 0, 3, 1, 0, 135, 1, 1.05, 0.95], 
 59                  [27, 2, 0, 0, 0, 0, 3, 1, 0, 135, 1, 1.1, 0.95], 
 60                  [28, 1, 0, 0, 0, 0, 1, 1, 0, 135, 1, 1.05, 0.95], 
 61                  [29, 1, 2.7897792, 1.0461672, 0, 0, 3, 1, 0, 135, 1, 1.05, 0.95], 
 62                  [30, 2, 0, 0, 0, 0, 3, 1, 0, 135, 1, 1.05, 0.95] 
 63          ]) 
 64   
 65      ## generator data 
 66      # bus, Pg, Qg, Qmax, Qmin, Vg, mBase, status, Pmax, Pmin 
 67      gen = array([ 
 68                  [1, 23.54, 0, 60, -15, 1, 100, 1, 60, 35], 
 69                  [2, 60.97, 0, 60, -15, 1, 100, 1, 60, 12], 
 70                  [220, 21.59, 0, 60, -15, 1, 100, 1, 60, 12], 
 71                  [27, 26.91, 0, 60, -15, 1, 100, 1, 60, 12], 
 72                  [23, 19.2, 0, 60, -15, 1, 100, 1, 60, 12], 
 73                  [13, 37, 0, 60, -15, 1, 100, 1, 60, 12], 
 74                  [7, -30, -15, 0, -15, 1, 100, 1, 0, -30], 
 75                  [15, -30, -12, 0, -12, 1, 100, 1, 0, -30], 
 76                  [30, -30, -7.5, 0, -7.5, 1, 100, 1, 0, -30] 
 77          ]) 
 78   
 79      ## branch data 
 80      # fbus, tbus, r, x, b, rateA, rateB, rateC, ratio, angle, status, angmin, angmax 
 81      branch = array([ 
 82                  [1, 2, 0.02, 0.06, 0.03, 10000, 130, 130, 0, 0, 1], 
 83                  [1, 3, 0.05, 0.19, 0.02, 10000, 130, 130, 0, 0, 1], 
 84                  [2, 4, 0.06, 0.17, 0.02, 10000, 65, 65, 0, 0, 1], 
 85                  [3, 4, 0.01, 0.04, 0, 10000, 130, 130, 0, 0, 1], 
 86                  [2, 5, 0.05, 0.2, 0.02, 10000, 130, 130, 0, 0, 1], 
 87                  [2, 6, 0.06, 0.18, 0.02, 10000, 65, 65, 0, 0, 1], 
 88                  [4, 6, 0.01, 0.04, 0, 10000, 90, 90, 0, 0, 1], 
 89                  [5, 7, 0.05, 0.12, 0.01, 10000, 70, 70, 0, 0, 1], 
 90                  [6, 7, 0.03, 0.08, 0.01, 10000, 130, 130, 0, 0, 1], 
 91                  [6, 8, 0.01, 0.04, 0, 10000, 32, 32, 0, 0, 1], 
 92                  [6, 9, 0, 0.21, 0, 10000, 65, 65, 0, 0, 1], 
 93                  [6, 10, 0, 0.56, 0, 10000, 32, 32, 0, 0, 1], 
 94                  [9, 11, 0, 0.21, 0, 10000, 65, 65, 0, 0, 1], 
 95                  [9, 10, 0, 0.11, 0, 10000, 65, 65, 0, 0, 1], 
 96                  [4, 12, 0, 0.26, 0, 10, 65, 65, 0, 0, 1], 
 97                  [12, 13, 0, 0.14, 0, 10000, 65, 65, 0, 0, 1], 
 98                  [12, 14, 0.12, 0.26, 0, 10000, 32, 32, 0, 0, 1], 
 99                  [12, 15, 0.07, 0.13, 0, 10000, 32, 32, 0, 0, 1], 
100                  [12, 16, 0.09, 0.2, 0, 10000, 32, 32, 0, 0, 1], 
101                  [14, 15, 0.22, 0.2, 0, 10000, 16, 16, 0, 0, 1], 
102                  [16, 17, 0.08, 0.19, 0, 10000, 16, 16, 0, 0, 1], 
103                  [15, 18, 0.11, 0.22, 0, 10000, 16, 16, 0, 0, 1], 
104                  [18, 19, 0.06, 0.13, 0, 10000, 16, 16, 0, 0, 1], 
105                  [19, 20, 0.03, 0.07, 0, 10000, 32, 32, 0, 0, 1], 
106                  [10, 20, 0.09, 0.21, 0, 10000, 32, 32, 0, 0, 0], 
107                  [10, 17, 0.03, 0.08, 0, 10000, 32, 32, 0, 0, 0], 
108                  [10, 21, 0.03, 0.07, 0, 10000, 32, 32, 0, 0, 1], 
109                  [10, 220, 0.07, 0.15, 0, 10000, 32, 32, 0, 0, 1], 
110                  [21, 220, 0.01, 0.02, 0, 10000, 32, 32, 0, 0, 1], 
111                  [15, 23, 0.1, 0.2, 0, 10000, 16, 16, 0, 0, 1], 
112                  [220, 24, 0.12, 0.18, 0, 10000, 16, 16, 0, 0, 1], 
113                  [23, 24, 0.13, 0.27, 0, 10, 16, 16, 0, 0, 1], 
114                  [24, 25, 0.19, 0.33, 0, 10000, 16, 16, 0, 0, 1], 
115                  [25, 26, 0.25, 0.38, 0, 10000, 16, 16, 0, 0, 1], 
116                  [25, 27, 0.11, 0.21, 0, 10000, 16, 16, 0, 0, 1], 
117                  [28, 27, 0, 0.4, 0, 10000, 65, 65, 0, 0, 1], 
118                  [27, 29, 0.22, 0.42, 0, 10000, 16, 16, 0, 0, 1], 
119                  [27, 30, 0.32, 0.6, 0, 10000, 16, 16, 0, 0, 1], 
120                  [29, 30, 0.24, 0.45, 0, 10000, 16, 16, 0, 0, 1], 
121                  [8, 28, 0.06, 0.2, 0.02, 10000, 32, 32, 0, 0, 1], 
122                  [6, 28, 0.02, 0.06, 0.01, 10000, 32, 32, 0, 0, 1] 
123          ]) 
124   
125      ##-----  OPF Data  -----## 
126      ## area data 
127      # area refbus 
128      areas = array([ 
129                  [1,     8], 
130                  [2, 3], 
131                  [3,     26] 
132          ]) 
133   
134      ## generator cost data 
135      # 1 startup shutdown n x1 y1 ... xn yn 
136      # 2 startup shutdown n c(n-1) ... c0 
137      gencost = array([ 
138                  [1, 0, 0, 4, 0, 0, 12, 240, 36, 1200, 60, 2400], 
139                  [1, 0, 0, 4, 0, 0, 12, 240, 36, 1200, 60, 2400], 
140                  [1, 0, 0, 4, 0, 0, 12, 240, 36, 1200, 60, 2400], 
141                  [1, 0, 0, 4, 0, 0, 12, 240, 36, 1200, 60, 2400], 
142                  [1, 0, 0, 4, 0, 0, 12, 240, 36, 1200, 60, 2400], 
143                  [1, 0, 0, 4, 0, 0, 12, 240, 36, 1200, 60, 2400], 
144                  [1, 0, 0, 4, -30, -3000, -20, -2000, -10, -1000, 0, 0], 
145                  [1, 0, 0, 4, -30, -3000, -20, -2000, -10, -1000, 0, 0], 
146                  [1, 0, 0, 4, -30, -3000, -20, -2000, -10, -1000, 0, 0] 
147          ], float) 
148   
149      return baseMVA, bus, gen, branch, areas, gencost 
150