Source Code for Module pypower.case6ww

 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 6 bus, 3 gen case from Wood & Wollenberg. 
18  """ 
19   
20  from numpy import array 
21   
22 -def case6ww(): 
23      """Power flow data for 6 bus, 3 gen case from Wood & Wollenberg. 
24      Please see L{caseformat} for details on the case file format. 
25   
26      This is the 6 bus example from pp. 104, 112, 119, 123-124, 549 of 
27      I{"Power Generation, Operation, and Control, 2nd Edition"}, 
28      by Allen. J. Wood and Bruce F. Wollenberg, John Wiley & Sons, NY, Jan 1996. 
29   
30      @return: Power flow data for 6 bus, 3 gen case from Wood & Wollenberg. 
31      """ 
32      ppc = {"version": '2'} 
33   
34      ##-----  Power Flow Data  -----## 
35      ## system MVA base 
36      ppc["baseMVA"] = 100.0 
37   
38      ## bus data 
39      # bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin 
40      ppc["bus"] = array([ 
41          [1, 3,  0,  0, 0, 0, 1, 1.05, 0, 230, 1, 1.05, 1.05], 
42          [2, 2,  0,  0, 0, 0, 1, 1.05, 0, 230, 1, 1.05, 1.05], 
43          [3, 2,  0,  0, 0, 0, 1, 1.07, 0, 230, 1, 1.07, 1.07], 
44          [4, 1, 70, 70, 0, 0, 1, 1,    0, 230, 1, 1.05, 0.95], 
45          [5, 1, 70, 70, 0, 0, 1, 1,    0, 230, 1, 1.05, 0.95], 
46          [6, 1, 70, 70, 0, 0, 1, 1,    0, 230, 1, 1.05, 0.95] 
47      ]) 
48   
49      ## generator data 
50      # bus, Pg, Qg, Qmax, Qmin, Vg, mBase, status, Pmax, Pmin, Pc1, Pc2, 
51      # Qc1min, Qc1max, Qc2min, Qc2max, ramp_agc, ramp_10, ramp_30, ramp_q, apf 
52      ppc["gen"] = array([ 
53          [1, 0,  0, 100, -100, 1.05, 100, 1, 200, 50,   0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], 
54          [2, 50, 0, 100, -100, 1.05, 100, 1, 150, 37.5, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], 
55          [3, 60, 0, 100, -100, 1.07, 100, 1, 180, 45,   0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0] 
56      ]) 
57   
58      ## branch data 
59      # fbus, tbus, r, x, b, rateA, rateB, rateC, ratio, angle, status, angmin, angmax 
60      ppc["branch"] = array([ 
61          [1, 2, 0.1,  0.2,  0.04, 40, 40, 40, 0, 0, 1, -360, 360], 
62          [1, 4, 0.05, 0.2,  0.04, 60, 60, 60, 0, 0, 1, -360, 360], 
63          [1, 5, 0.08, 0.3,  0.06, 40, 40, 40, 0, 0, 1, -360, 360], 
64          [2, 3, 0.05, 0.25, 0.06, 40, 40, 40, 0, 0, 1, -360, 360], 
65          [2, 4, 0.05, 0.1,  0.02, 60, 60, 60, 0, 0, 1, -360, 360], 
66          [2, 5, 0.1,  0.3,  0.04, 30, 30, 30, 0, 0, 1, -360, 360], 
67          [2, 6, 0.07, 0.2,  0.05, 90, 90, 90, 0, 0, 1, -360, 360], 
68          [3, 5, 0.12, 0.26, 0.05, 70, 70, 70, 0, 0, 1, -360, 360], 
69          [3, 6, 0.02, 0.1,  0.02, 80, 80, 80, 0, 0, 1, -360, 360], 
70          [4, 5, 0.2,  0.4,  0.08, 20, 20, 20, 0, 0, 1, -360, 360], 
71          [5, 6, 0.1,  0.3,  0.06, 40, 40, 40, 0, 0, 1, -360, 360] 
72      ]) 
73   
74      ##-----  OPF Data  -----## 
75      ## generator cost data 
76      # 1 startup shutdown n x1 y1 ... xn yn 
77      # 2 startup shutdown n c(n-1) ... c0 
78      ppc["gencost"] = array([ 
79          [2, 0, 0, 3, 0.00533, 11.669, 213.1], 
80          [2, 0, 0, 3, 0.00889, 10.333, 200], 
81          [2, 0, 0, 3, 0.00741, 10.833, 240] 
82      ]) 
83   
84      return ppc 
85