Power Flow Analysis PDF
Power Flow Analysis PDF
Power Flow Analysis PDF
=
+ + =
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.