-
Notifications
You must be signed in to change notification settings - Fork 466
/
idx_gen.py
85 lines (74 loc) · 4.31 KB
/
idx_gen.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
# Copyright 1996-2015 PSERC. All rights reserved.
# Use of this source code is governed by a BSD-style
# license that can be found in the LICENSE file.
""" Defines constants for named column indices to gen matrix.
Some examples of usage, after defining the constants using the line above,
are::
Pg = gen[3, PG] # get the real power output of generator 4
gen[:, PMIN] = 0 # set to zero the minimum real power limit of all gens
The index, name and meaning of each column of the gen matrix is given
below:
columns 0-20 must be included in input matrix (in case file)
0. C{GEN_BUS} bus number
1. C{PG} real power output (MW)
2. C{QG} reactive power output (MVAr)
3. C{QMAX} maximum reactive power output (MVAr)
4. C{QMIN} minimum reactive power output (MVAr)
5. C{VG} voltage magnitude setpoint (p.u.)
6. C{MBASE} total MVA base of machine, defaults to baseMVA
7. C{GEN_STATUS} 1 - in service, 0 - out of service
8. C{PMAX} maximum real power output (MW)
9. C{PMIN} minimum real power output (MW)
10. C{PC1} lower real power output of PQ capability curve (MW)
11. C{PC2} upper real power output of PQ capability curve (MW)
12. C{QC1MIN} minimum reactive power output at Pc1 (MVAr)
13. C{QC1MAX} maximum reactive power output at Pc1 (MVAr)
14. C{QC2MIN} minimum reactive power output at Pc2 (MVAr)
15. C{QC2MAX} maximum reactive power output at Pc2 (MVAr)
16. C{RAMP_AGC} ramp rate for load following/AGC (MW/min)
17. C{RAMP_10} ramp rate for 10 minute reserves (MW)
18. C{RAMP_30} ramp rate for 30 minute reserves (MW)
19. C{RAMP_Q} ramp rate for reactive power (2 sec timescale) (MVAr/min)
20. C{APF} area participation factor
columns 21-24 are added to matrix after OPF solution
they are typically not present in the input matrix
(assume OPF objective function has units, u)
21. C{MU_PMAX} Kuhn-Tucker multiplier on upper Pg limit (u/MW)
22. C{MU_PMIN} Kuhn-Tucker multiplier on lower Pg limit (u/MW)
23. C{MU_QMAX} Kuhn-Tucker multiplier on upper Qg limit (u/MVAr)
24. C{MU_QMIN} Kuhn-Tucker multiplier on lower Qg limit (u/MVAr)
@author: Ray Zimmerman (PSERC Cornell)
@author: Richard Lincoln
"""
# define the indices
GEN_BUS = 0 # bus number
PG = 1 # Pg, real power output (MW)
QG = 2 # Qg, reactive power output (MVAr)
QMAX = 3 # Qmax, maximum reactive power output at Pmin (MVAr)
QMIN = 4 # Qmin, minimum reactive power output at Pmin (MVAr)
VG = 5 # Vg, voltage magnitude setpoint (p.u.)
MBASE = 6 # mBase, total MVA base of this machine, defaults to baseMVA
GEN_STATUS = 7 # status, 1 - machine in service, 0 - machine out of service
PMAX = 8 # Pmax, maximum real power output (MW)
PMIN = 9 # Pmin, minimum real power output (MW)
PC1 = 10 # Pc1, lower real power output of PQ capability curve (MW)
PC2 = 11 # Pc2, upper real power output of PQ capability curve (MW)
QC1MIN = 12 # Qc1min, minimum reactive power output at Pc1 (MVAr)
QC1MAX = 13 # Qc1max, maximum reactive power output at Pc1 (MVAr)
QC2MIN = 14 # Qc2min, minimum reactive power output at Pc2 (MVAr)
QC2MAX = 15 # Qc2max, maximum reactive power output at Pc2 (MVAr)
RAMP_AGC = 16 # ramp rate for load following/AGC (MW/min)
RAMP_10 = 17 # ramp rate for 10 minute reserves (MW)
RAMP_30 = 18 # ramp rate for 30 minute reserves (MW)
RAMP_Q = 19 # ramp rate for reactive power (2 sec timescale) (MVAr/min)
APF = 20 # area participation factor
# included in opf solution, not necessarily in input
# assume objective function has units, u
MU_PMAX = 21 # Kuhn-Tucker multiplier on upper Pg limit (u/MW)
MU_PMIN = 22 # Kuhn-Tucker multiplier on lower Pg limit (u/MW)
MU_QMAX = 23 # Kuhn-Tucker multiplier on upper Qg limit (u/MVAr)
MU_QMIN = 24 # Kuhn-Tucker multiplier on lower Qg limit (u/MVAr)
# Note: When a generator's PQ capability curve is not simply a box and the
# upper Qg limit is binding, the multiplier on this constraint is split into
# it's P and Q components and combined with the appropriate MU_Pxxx and
# MU_Qxxx values. Likewise for the lower Q limits.