PSOC Optimal Power
PSOC Optimal Power
& Control
Fall-2020
1
Optimal Power-Flow
2
Optimal Power-Flow
3
Optimal Power-Flow
4
Optimal Power-Flow
• The total losses of the power system are simply part of
the power flow calculation and are reflected in the
loading of the generation on the reference bus—thus,
• There is no need to specifically calculate the losses
since they are inherently a part of the power flow.
• More importantly, the ED can be constrained to meet
transmission system limits such as MW or MVA flow
limits on lines or transformers or voltage limits on buses.
• The result is the generation dispatch representing the
minimum $/h total generation cost and that also solves
the power flow at that optimum.
5
Optimal Power-Flow
6
Optimal Power-Flow
Does not Include e Cost
7
Optimal Power-Flow
8
Optimal Power-Flow
9
Optimal Power-Flow
WITHOUT INEQUALITY CONSTRAINTS
The primary objective of the optimal power flow
solution is to minimize the overall cost of generation.
This is represented by an objective function (or) cost
function as
Eq. 1
Eq. 2
Eq. 3
10
Optimal Power-Flow
WITHOUT INEQUALITY CONSTRAINTS
11
• The net real power injected into the system at
the ith bus
𝑃𝑖 = 𝑃𝐺𝑖 − 𝑃𝐷𝑖
𝑄𝑖 = 𝑄𝐺𝑖 − 𝑄𝐷𝑖
12
Optimal Power-Flow
WITHOUT INEQUALITY CONSTRAINTS
For load buses (or) P–Q buses, P and Q are specified and
hence Equations (2) and (3) form the equality constraints.
For P–V buses, P and |V| are specified as the function of
some vectors and are represented as
13
Optimal Power-Flow
WITHOUT INEQUALITY CONSTRAINTS
14
Optimal Power-Flow
WITHOUT INEQUALITY CONSTRAINTS
• Out of these independent variables certain variables are
chosen as control variables
• Which are to be varied to yield an optimal value of the
objective function.
• The remaining independent variables are called fixed or
disturbance or uncontrollable parameters
Let ‘u’ be the vector of control variables and ‘p’ the vector
of fixed or disturbance variables
15
• Hence, the vector of independent variables can
be represented as the combination of vector of
control variables ‘u’ and vector uncontrollable
variables ‘p’ and is expressed as
16
• The control parameters may be:
1. Voltage magnitude at P–V buses,
2. P at generator buses with controllable power
3. Slack bus voltage and regulating transformer
tap setting as additional control variables
4. In the case of buses with reactive-power control,
𝑄𝐺𝑖 is taken as a control variable.
Now, the optimal power flow problem can be stated as
17
Now, the optimal power flow problem can be stated as
20
21
22
Optimal Power Flow Solution by N–R method
23
Optimal Power Flow Solution by N–R method
+higher
Terms
24
Optimal Power Flow Solution by N–R method
Eq. 7
25
Optimal Power Flow Solution by N–R method
Eq.8
26
Optimal Power Flow Solution by N–R method
27
Optimal Power Flow Solution by N–R method
28
Optimal Power Flow Solution by N–R method
29
Optimal Power Flow Solution by N–R method
30
Algorithm for computational procedure
32
Algorithm for computational procedure
Eq. 9
33
Algorithm for computational procedure
36
Inequality constraints on control variables
otherwise set
37
Inequality constraints on control variables
38
Inequality constraints on control variables
39
OPTIMAL POWER FLOW USING THE
DC POWER FLOW
40
Power Flow Equations
N
PkI - å Vk Vi[Gki cos(q k - q i ) + B ki sin(q k - q i )] = 0
i=1
N
QkI - å Vk Vi[Gki sin(q k - q i ) - Bki cos(q k - q i )] = 0
i=1
N
PkI - å Vk Vi[Gki cos(q k - q i ) + B ki sin(q k - q i )] = 0
i=1
Assumptions
Neglect Resistance of the Branches
N
PkI - å Vk Vi[Gki cos(q k - q i ) + B ki sin(q k - q i )] = 0
i=1
N
PkI - å Vk Vi Bki sin(q k - q i ) = 0
i=1
Assumptions.
All Voltage Magnitudes = 1.0 p.u
N
PkI - å Vk Vi Bki sin(q k - q i ) = 0
i=1
N
PkI
- å Bki sin(q k - q i ) = 0
i=1
Assumptions
All angles are small N
PkI - å Bki sin(q k - q i ) = 0
i=1
N N
(q k - q i )
PkI - å Bki (q k - q i ) = 0 or PkI -å =0
i=1 i=1 x ki
Interpretation
N
(q k - q i )
PkI -å =0
i=1 x ki P1I P2I
P12
q1 q2
N
x12
PkI - å Pki = 0 x13 x 23
i=1 P31 P23
(q - q i ) q3
Pki = k
x ki P3I
(q1 - q 2 ) (q 2 - q 3 ) (q 3 - q1 )
P12 = ; P23 = ; P31 =
x12 x 23 x13
Why is it called dc power flow?
(q k - q i ) (Vk - Vi )
Pki = I ki =
x ki R ki
DC power flow
• Start with the objective function
89
Example
A B CA
1 2 10 $/MWh
PA
PAMAX=390MW
3
CB
450 MW
Economic dispatch:
PA = P = 390 MW
max 20$/MWh
A PB
PB = 60 MW PBMAX= 150 MW
90
Flows resulting from the economic dispatch
390MW 60MW
A B
Flowmax = 200MW
1 2
450 MW
91
Calculating the flows using superposition
Because we assume a linear model, superposition is
applicable
390 MW 60 MW
1 2
3
450 MW
390 MW 60 MW
1 2 1 2
3 3
390 MW
60 MW
92
Calculating the flows using superposition (1)
390 MW
FA FB
1 2
3
390 MW
FA + FB = 390 MW
FA = 260 MW
FB = 130 MW
93
Calculating the flows using superposition (2)
60 MW
FD FC
1 2
3
60 MW
FC + FD = 60 MW
FC = 40 MW
FD = 20 MW
94
Calculating the flows using superposition (2)
60 MW
FD FC
1 2
3
60 MW
FC + FD = 60 MW
FC = 40 MW
FD = 20 MW
95
Calculating the flows using superposition (3)
390 MW 60 MW
130 MW 20 MW
260 MW 40 MW
1 2 1 2
3 3
390 MW 60 MW
390 MW 60 MW
110 MW
280 MW 170 MW
1 2
Fmax = 260 MW
3
450 MW
96
Correcting unacceptable flows
390 MW 60 MW
280 MW
1 2
110 MW 40 MW
220 MW 80 MW
1 2 1 2
3 3
330 MW 120 MW
330 MW 120 MW
70 MW
260 MW 190 MW
1 2
Fmax = 260 MW
3
450 MW
98
Comments (1)
99
Comments (2)
100
Security Constrained OPF (SCOPF)
101
Example: base case solution of OPF
330MW 120MW
A 70 MW B
1 2
260 MW 190 MW
450 MW
102
Example: contingency case
330MW 120MW
A 0 MW B
1 2
330 MW 120 MW
450 MW
103
Formulation of the Security Constrained OPF
uk - u0 £ Du max