Source Code for Module pypower.t.t_case9_opf

 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 9 bus, 3 generator case. 
18  """ 
19   
20  from numpy import array 
21   
22   
23 -def t_case9_opf(): 
24      """Power flow data for 9 bus, 3 generator case. 
25      Please see L{caseformat} for details on the case file format. 
26   
27      @return: Power flow data for 9 bus, 3 generator case, with OPF data. 
28      """ 
29      ##-----  Power Flow Data  -----## 
30      ## system MVA base 
31      baseMVA = 100.0 
32   
33      ## bus data 
34      # bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin 
35      bus = array([ 
36          [1,  3, 0,    0, 0, 0, 1, 1, 0, 345, 1, 1.1, 0.9], 
37          [2,  2, 0,    0, 0, 0, 1, 1, 0, 345, 1, 1.1, 0.9], 
38          [30, 2, 0,    0, 0, 0, 1, 1, 0, 345, 1, 1.1, 0.9], 
39          [4,  1, 0,    0, 0, 0, 1, 1, 0, 345, 1, 1.1, 0.9], 
40          [5,  1, 90,  30, 0, 0, 1, 1, 0, 345, 1, 1.1, 0.9], 
41          [6,  1, 0,    0, 0, 0, 1, 1, 0, 345, 1, 1.1, 0.9], 
42          [7,  1, 100, 35, 0, 0, 1, 1, 0, 345, 1, 1.1, 0.9], 
43          [8,  1, 0,    0, 0, 0, 1, 1, 0, 345, 1, 1.1, 0.9], 
44          [9,  1, 125, 50, 0, 0, 1, 1, 0, 345, 1, 1.1, 0.9] 
45      ]) 
46   
47      ## generator data 
48      # bus, Pg, Qg, Qmax, Qmin, Vg, mBase, status, Pmax, Pmin, Pc1, Pc2, 
49      # Qc1min, Qc1max, Qc2min, Qc2max, ramp_agc, ramp_10, ramp_30, ramp_q, apf 
50      gen = array([ 
51          [1,  0,   0, 300, -300, 1, 100, 1, 250, 90], 
52          [2,  163, 0, 300, -300, 1, 100, 1, 300, 10], 
53          [30, 85,  0, 300, -300, 1, 100, 1, 270, 10] 
54      ], float) 
55   
56      ## branch data 
57      # fbus, tbus, r, x, b, rateA, rateB, rateC, ratio, angle, status, angmin, angmax 
58      branch = array([ 
59          [1,  4, 0,      0.0576, 0,       0, 250, 250, 0, 0, 1], 
60          [4,  5, 0.017,  0.092,  0.158,   0, 250, 250, 0, 0, 1], 
61          [5,  6, 0.039,  0.17,   0.358, 150, 150, 150, 0, 0, 1], 
62          [30, 6, 0,      0.0586, 0,       0, 300, 300, 0, 0, 1], 
63          [6,  7, 0.0119, 0.1008, 0.209,  40, 150, 150, 0, 0, 1], 
64          [7,  8, 0.0085, 0.072,  0.149, 250, 250, 250, 0, 0, 1], 
65          [8,  2, 0,      0.0625, 0,     250, 250, 250, 0, 0, 1], 
66          [8,  9, 0.032,  0.161,  0.306, 250, 250, 250, 0, 0, 1], 
67          [9,  4, 0.01,   0.085,  0.176, 250, 250, 250, 0, 0, 1] 
68      ]) 
69   
70      ##-----  OPF Data  -----## 
71      ## area data 
72      # area refbus 
73      areas = array([ 
74          [1, 5] 
75      ]) 
76   
77      ## generator cost data 
78      # 1 startup shutdown n x1 y1 ... xn yn 
79      # 2 startup shutdown n c(n-1) ... c0 
80      gencost = array([ 
81          [1, 0, 0, 4,  0,        0,   100, 2500, 200, 5500, 250, 7250], 
82          [2, 0, 0, 2, 24.035, -403.5,   0,    0,   0,    0,   0,    0], 
83          [1, 0, 0, 3,  0,        0,   200, 3000, 300, 5000,   0,    0] 
84      ]) 
85   
86      return baseMVA, bus, gen, branch, areas, gencost 
87