Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Power Flow Analysis PDF

Download as pdf or txt
Download as pdf or txt
You are on page 1of 33

Power Flow Analysis

Well known as : Load Flow


The Power Flow Problem
Power flow analysis is fundamental to the study of power
systems.
In fact, power flow forms the core of power system analysis.
power flow study plays a key role in the planning of additions
or expansions to transmission and generation facilities.
A power flow solution is often the starting point for many other
types of power system analyses.
In addition, power flow analysis is at the heart of contingency
analysis and the implementation of real-time monitoring
systems.
Problem Statement
For a given power network, with known complex
power loads and some set of specifications or
restrictions on power generations and voltages, solve
for any unknown bus voltages and unspecified
generation and finally for the complex power flow in the
network components.
Network Structure
Power Flow Study Steps
1. Determine element values for passive network
components.
2. Determine locations and values of all complex power
loads.
3. Determine generation specifications and constraints.
4. Develop a mathematical model describing power flow
in the network.
5. Solve for the voltage profile of the network.
6. Solve for the power flows and losses in the network.
7. Check for constraint violations.
Formulation of the Bus Admittance Matrix
The first step in developing the mathematical model describing
the power flow in the network is the formulation of the bus
admittance matrix.
The bus admittance matrix is an n*n matrix (where n is the
number of buses in the system) constructed from the
admittances of the equivalent circuit elements of the segments
making up the power system.
Most system segments are represented by a combination of
shunt elements (connected between a bus and the reference
node) and series elements (connected between two system
buses).
Bus Admittance Matrix
Formulation of the bus admittance matrix follows two
simple rules:
1. The admittance of elements connected between node k
and reference is added to the (k, k) entry of the
admittance matrix.
2. The admittance of elements connected between nodes
j and k is added to the (j, j) and (k, k) entries of the
admittance matrix.
The negative of the admittance is added to the (j, k) and
(k, j) entries of the admittance matrix.
Bus Admittance Matrix
Bus Admittance Matrix
Node-Voltage Equations
Applying KCL at each node yields: Defining the Ys as
The Y-Bus
The current equations reduced to
In a compact form
Where,
The General Form of the Load-Flow
Equations
In Practice, bus powers S
i
is specified rather than the bus
currents I
i
.
As a result, we have
) ( | | ) (
1 1
*
*
i n in n i in
N
n
N
n
n in
i
i i i i
V V Y V Y V I V jQ P + = = =

= =
ij ij ij ij ij ij ij ij ij
jB G Y j Y Y Y + = + = = sin | | cos | | | |
i
i
*
V
S
I =
Load-Flow Equations
These are the static power flow equations. Each equation is
complex, and therefore we have 2n real equations. The nodal
admittance matrix current equation can be written in the power
form:
) ( | | ) (
1 1
*
*
i n in n i in
N
n
N
n
n in
i
i i i i
V V Y V Y V I V jQ P + = = =

= =
Let,
ij ij ij ij ij ij ij ij ij
jB G Y j Y Y Y + = + = = sin | | cos | | | |
Load-Flow Equations
Finally,

=
+ + =
N
i n
n
i n in in n i ii i i
Y V V G V P
1
2
) cos( | | | |
o This is known as NR (Newton Raphson) formulation

=
+ =
N
i n
n
i n in in n i ii i i
Y V V B V Q
1
2
) sin( | | | |
Gauss Power Flow
*
* * *
i
1 1
* * * *
1 1
*
*
1 1,
*
*
1,
We first need to put the equation in the correct form
S
S
S
S 1
i i
i
i
n n
i i i ik k i ik k
k k
n n
i i i ik k ik k
k k
n n
i
ik k ii i ik k
k k k i
n
i
i ik k
ii k k i
V I V Y V V Y V
V I V Y V V Y V
Y V Y V Y V
V
V Y V
Y
V
= =
= =
= =
=
| |
= = =
|
\ .
= = =
= = +
|
=


|
|
|
\ .
Difficulties
Unless the generation equals the load at every bus, the
complex power outputs of the generators cannot be arbitrarily
selected.
In fact, the complex power output of at least one of the
generators must be calculated last, since it must take up the
unknown slack due to the uncalculated network losses.
Further, losses cannot be calculated until the voltages are
known.
Also, it is not possible to solve these equations for the
absolute phase angles of the phasor voltages. This simply
means that the problem can only be solved to some arbitrary
phase angle reference.
Difficulties
For a 4- bus system, suppose that S
G4
is arbitrarily allowed
to float or swing (in order to take up the necessary slack
caused by the losses) and that S
G1,
S
G2,
S
G3
are specified.
Remedies
Now, with the loads known, the equations are seen as
four simultaneous nonlinear equations with complex
coefficients in five unknowns. (V
1
, V
2
, V
3
, V
4
and S
G4
).
Designating bus 4 as the slack bus and specifying the
voltage V
4
reduces the problem to four equations in four
unknowns.
Remedies
The slack bus is chosen as the phase reference for all
phasor calculations, its magnitude is constrained, and the
complex power generation at this bus is free to take up
the slack necessary in order to account for the system
real and reactive power losses.
Systems of nonlinear equations, cannot (except in rare
cases) be solved by closed-form techniques.
Load Flow Solution
There are four quantities of interest associated with each bus:
1. Real Power, P
2. Reactive Power, Q
3. Voltage Magnitude, V
4. Voltage Angle,
At every bus of the system, two of these four quantities will be
specified and the remaining two will be unknowns.
Each of the system buses may be classified in accordance with
which of the two quantities are specified
Bus Classifications
Slack Bus The slack bus for the system is a single bus for which the
voltage magnitude and angle are specified.
The real and reactive power are unknowns.
The bus selected as the slack bus must have a source of both real
and reactive power, since the injected power at this bus must swing
to take up the slack in the solution.
The best choice for the slack bus (since, in most power systems,
many buses have real and reactive power sources) requires
experience with the particular system under study.
The behavior of the solution is often influenced by the bus chosen.
Bus Classifications
Load Bus (P-Q Bus) : A load bus is defined as any bus of the system
for which the real and reactive power are specified.
Load buses may contain generators with specified real and reactive
power outputs;
however, it is often convenient to designate any bus with specified
injected complex power as a load bus.
Voltage Controlled Bus (P-V Bus) : Any bus for which the voltage
magnitude and the injected real power are specified is classified as a
voltage controlled (or P-V) bus.
The injected reactive power is a variable (with specified upper and
lower bounds) in the power flow analysis.
(A P-V bus must have a variable source of reactive power such as a
generator.)
Solution Methods
The solution of the simultaneous nonlinear power flow
equations requires the use of iterative techniques for
even the simplest power systems.
There are many methods for solving nonlinear equations,
such as:
- Gauss Seidel.
- Newton Raphson.
- Fast Decoupled.
Guess Solution
It is important to have a good approximation to the load-
flow solution, which is then used as a starting estimate
(or initial guess) in the iterative procedure.
A fairly simple process can be used to evaluate a good
approximation to the unknown voltages and phase
angles.
The process is implemented in two stages: the first
calculates the approximate angles, and the second
calculates the approximate voltage magnitudes.
Gauss Iteration Method
Gauss Iteration Example
( 1) ( )
(0)
( ) ( )
Example: Solve - 1 0
1
Let k = 0 and arbitrarily guess x 1 and solve
0 1 5 2.61185
1 2 6 2.61612
2 2.41421 7 2.61744
3 2.55538 8 2.61785
4 2.59805 9 2.61798
v v
v v
x x
x x
k x k x
+
=
= +
=
Stopping Criteria
=
Gauss Power Flow
*
* * *
i
1 1
* * * *
1 1
*
*
1 1,
*
*
1,
We first need to put the equation in the correct form
S
S
S
S 1
i i
i
i
n n
i i i ik k i ik k
k k
n n
i i i ik k ik k
k k
n n
i
ik k ii i ik k
k k k i
n
i
i ik k
ii k k i
V I V Y V V Y V
V I V Y V V Y V
Y V Y V Y V
V
V Y V
Y
V
= =
= =
= =
=
| |
= = =
|
\ .
= = =
= = +
|
=


|
|
|
\ .
Example
A 100 MW, 50 Mvar load is connected to a generator through a line
with z = 0.02 + j0.06 p.u. and line charging of 0.05 p.u on each end
(100 MVA base). Also, there is a 0.25 p.u. capacitance at bus 2. If the
generator voltage is 1.0 p.u., what is V
2
?
100+j50
0.25 p.u.
Z = 0.02 + j0.06
V= 1 0
V
2
Y-Bus
2
2 bus
bus
22
The unknown is the complex load voltage, V .
To determine V we need to know the .
1
5 15
0.02 0.06
5 14.95 5 15
Hence
5 15 5 14.70
( Note - 15 0.05 0.25)
j
j
j j
j j
B j j j
=
+
+
(
=
(
+

= + +
Y
Y
Solution
*
2
2
*
22 1,
2
2
*
2
(0)
2
( ) ( )
2 2
1 S
1 -1 0.5
( 5 15)(1.0 0)
5 14.70
Guess 1.0 0 (this is known as a flat start)
0 1.000 0.000 3 0.9622 0.0556
1 0.9671 0.0568 4 0.9622 0.0556
2 0
n
ik k
k k i
v v
V Y V
Y
V
j
V j
j
V
V
v V v V
j j
j j
=
| |
=
|
\ .
| |
+
= +
|

\ .
=
+

.9624 0.0553 j
Solution (cont.)
Gauss-Seidel Iteration
( 1) ( ) ( ) ( )
2 1
2 2 3
( 1) ( 1) ( ) ( )
2 1 3 2 3
( 1) ( 1) ( 1) ( ) ( )
2 1 4 2 3 4
( 1) ( 1) ( 1) ( 1) ( )
2 1 2 3 4
Immediately use the new voltage estimates:
( , , , , )
( , , , , )
( , , , , )
( , , , ,
v v v v
n
v v v v
n
v v v v v
n
v v v v v
n n
V h V V V V
V h V V V V
V h V V V V V
V h V V V V V
+
+ +
+ + +
+ + + +
=
=
=
=

)
The Gauss-Seidel works better than the Gauss, and
is actually easier to implement. It is used instead
of Gauss.

You might also like