1 - 3a Load Flow Calculations - Theory PDF
1 - 3a Load Flow Calculations - Theory PDF
The load or PQ-bus: For these buses, the net real and reactive power is known. PQ-buses normally do not
have generators. However, if the reactive power of a generator reaches its limit, the corresponding bus is
treated as a PQ-bus. This is equivalent to adjusting the bus voltage until the generator reactive power falls
within the prescribed limits.
Distribution substations and feeders may be treated as generating buses in distribution networks.
The load flow equations are written in terms of the net power injection to each bus. With reference to
Figure 1.3.1, the net power injection into the kth bus is the combination of generated and load power. The
power flowing out of this bus must equal the net injected power. Therefore, the power balance equation at
the kth bus is written as follows in terms of the system voltage.
(1.3.1)
where
N is the number of network buses,
Pk is the net real power injected into the kth bus,
Qk is the net reactive power injected into the kth bus,
Yk,i is the total admittance between bus k and i: this total can be found from the bus admittance matrix,
Ybus, of the system,
Vi is the voltage of the ith bus.
The values of network admittances and voltages in Equation (1.3.1) are expressed in the pu system to
facilitate the calculations. To review the pu system definitions, see Section 1.1.
Equation (1.3.1) can be written separately for the real and reactive power.
Therefore,
P |V V Y | cos +
k
k
n
k , n|
k , n n k
n |
(1.3.2)
Q |V V Y | sin +
k
k
n
k , n|
k , n n k
n |
where k,n is the angle of the admittance, Yk,n, and j is the voltage phase angle at bus, j.
A real power equation is written for every PV- and PQ-bus and a reactive power equation is written for
every PQ-bus. Thus, for a power system with N buses of which L are PQ-buses there are (N-1) real
power equations (excluding the slack bus) and L reactive power equations (a total of N-1+L equations).
The unknowns are the magnitude and phase angle of the L PQ-bus voltages and the phase angle of the
(N-1-L) PV-bus voltages (a total of N-1+L unknowns).
The left-hand side of these equations are known and an iterative process is used for finding the
unknown voltages and phase angles such that Equations (1.3.2) are balanced.
Q fq ( , V)
(1.3.3)
where
P is the vector of the known net real power injections at PV- and PQ-buses,
Q is the vector of the known reactive power injections at PQ-buses,
V is the vector of the unknown bus voltage magnitudes,
is the vector of the unknown bus voltage phase angles, and
fp, fq are functions defined according to Equations (3.1.2).
Solution of the load flow problem requires finding the values of V and such that the right-hand side of
Equation (1.3.3) equals the known power injections at the network buses. For any estimation of V and ,
the difference between the known power injections, P and Q and the power injections calculated by
Equation (1.3.3) is called the power mismatch.
P fp ( , V)
S
Q fq ( , V)
(1.3.4)
P
DQ
The power mismatch is a measure of how close to the solution the estimations of V and are. A correction
to these estimations is obtained using the Newton-Raphson method, resulting in an iterative calculation
process.
j + 1
+ J S
(1.3.5)
j + 1
j + 1
+ J S
where the superscript, j, denotes variables calculated at the jth iteration step. J is the Jacobian matrix of
Equations (1.3.3):
d
p
p
d
dV
J
d
d
fq
fq
d
dV
(1.3.6)
At the beginning of a new iteration, (j+1), the power mismatch is calculated using Equation (1.3.4) with
the variables V and obtained from the previous iteration, j.
From the calculated mismatches, a new approximation to the system solution is found using Equation
(1.3.5).
The iteration process continues until the power mismatch at the jth step is smaller than a preset number .
To start the above iterative solution, an estimation of the unknown voltages and their phase angles is
required. This first solution approximation is called initial guess. Typically, the initial guess for the voltage
magnitudes is 1 pu and for their phase angles is 0 degrees (or radians).