PSA Chapter 2
PSA Chapter 2
By Dawit Adane
Load Flow Analysis
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.
Assumes balanced three phase system
Modeled as a single phase system (based on single line diagram)
A set of non-linear differential equations model both the Real (watts) and Reactive
(Vars) power flow
Matrices are developed for all impedances/admittances of transmission lines
interconnecting substations (buses)
Non-linear equations are solved through an iterative process, with an assumed initial
conditions
The behavior of the LF solution is often influenced by the bus chosen.
Importance of load-flow studies
• Great importance of power flow or load-flow studies is in the planning the future
expansion of power systems as well as in determining the best operation of
existing systems.
• Applications:
1. On-line analyses
• State estimation
The following questions are Criterion:
• Security
1. Will the iteration procedure converge to the
• Economic analyses
unique solution?
2. Off-line analyses
2. What is the convergence rate (how many
• Operation analyses
iterations are required)?
• Planning analyses
3. When using a digital computer, what are the
Network expansion planning
Power exchange planning computer storage and time requirements?
Security and adequacy analyses 4. Simplicity for programming and incorporating
-Faults different extra technique.
-Stability
Modeling of power system components
Loads
• can be classified into three categories of model;
i) constant power, 𝑃𝑝 = 𝐼2𝑅 ∝ 𝑉0 (𝑑𝑜𝑒𝑠 𝑛𝑜𝑡 𝑣𝑎𝑟𝑦 𝑤𝑖𝑡ℎ 𝑣𝑜𝑙𝑡𝑎𝑔𝑒)
ii) constant impedance, 𝑃𝑌 = 𝑉2/𝑅 ∝ 𝑉2 (𝑣𝑎𝑟𝑖𝑒𝑠 𝑤𝑖𝑡ℎ 𝑣𝑜𝑙𝑡𝑎𝑔𝑒 𝑠𝑞𝑢𝑎𝑟𝑒)
iii) constant current, 𝑃𝐼 = 𝐼𝑉 ∝ 𝑉1 (𝑣𝑎𝑟𝑖𝑒𝑠 𝑤𝑖𝑡ℎ 𝑣𝑜𝑙𝑡𝑎𝑔𝑒).
• To compute the AC power flow analysis under the normal steady-state values of the
bus voltages, the loads are always represented as constant power loads at a particular
time instant.
Transmission line
• The transmission lines are generally of medium length or of long length.
• A medium and long length lines are always represented by the nominal- model and the
equivalent- model respectively as shown in figure below.
Normal model of a medium line connected between buses ‘i’ and ‘j’
Equivalent model of a long transmission line connected between buses ‘i’ and ‘j’
Transformer
• The exciting current of the transformer is neglected.
• A two winding transformer, with tap ratio is 1:1, is represented by its per unit leakage
impedance as shown in figure below.
For a regulating transformer with transformation ratio 1:t or a:1, the equivalent circuit of
the transformer is shown in figures below. y is the per unit admittance of the transformer.
Where,
IBUS = [I1, I2, …, I5]T (5 × 1) is the vector of bus injection currents
VBUS = [V1, V2, …, V5]T (5 × 1) is the vector of bus voltages measured with
respect to the ground
YBUS (5 × 5) is the bus admittance matrix
• From the elements of the YBUS it can be observed that for i = 1, 2, …, 5;
Yii = sum total of all the admittances connected at bus ‘i’ called Self driving admittance
– diagonal elements
Yij = negative of the admittance connected between bus ‘i’ & ‘j’ (if these two buses are
physically connected with each other) called Mutual driving admittance – off-diagonal
element
Yij = 0; if there is no physical connection between buses ‘i’ and ‘j’
• Similarly, for a ‘n’ bus power system, where,
IBUS = [I1, I2, …, In]T (n × 1) is the vector of bus injection currents
VBUS = [V1, V2, …, Vn]T (n × 1) is the vector of bus voltages
YBUS (n × n) is the bus admittance matrix
Read about the Formation of YBUS matrix in the presence of mutually coupled elements
Basic Power Flow Equation
• For a ‘n’ bus system,
Or,
j=1 j=1
n n
Pi −jQi
• From the relation, Pi + jQ i = Vi Ii∗ we get, Ii = .
V∗i
Thus, 𝐧
𝟏 𝐏𝐢 − 𝐣𝐐𝐢
𝐕𝐢 = ∗ − 𝐘𝐢𝐤 𝐕𝐤
𝐘𝐢𝐢 𝐕𝐢
𝐤=𝟏, ≠𝐢
Here note that when computation is carried out for bus-i, updated values are already
available for buses 2,3, …, (i-1) in the current (k+1)st iteration. Hence these values are
used. For buses (i+1), ..., n, values from previous, kth iteration are used.
6. Continue iterations till
Where, ε is the tolerance limit. Generally, it is customary to use a value of 0.0001 pu.
7. Compute slack bus power after voltages have converged. [assuming bus 1 is slack bus].
9. The complex power loss in the line is given by Sik + Ski. The total loss in the system is
calculated by summing the loss over all the lines.
Ex
Ex
Case (b): Systems with PV buses also present:
For a system having multiple generators, the bus voltage initialization is carried out in a
two step procedure;
i) The load buses are initialized with flat start (i.e. Vj(0) = 1.000 for j = (m+ 1), (m+ 2),
…, n) and
ii) The magnitudes of the voltages of the PV buses are initialized with the corresponding
specified voltage magnitudes while initializing all these voltage angles to 00 (i.e. Vj(0) =
Vjsp00 for j = 2, 3, …, m, where Vjsp is the specified bus voltage magnitude of the jth
generator).
The reactive power supplied or absorbed by a generator (QG) is calculated by the load flow
procedure. However, any generator has a maximum and minimum limit on Q. If the
calculated QG from the generator is within these limits, then the generator excitation
system is able to maintain the terminal voltage at the specified value. On the other hand, if
the generator QG reaches or exceeds its limit on Q (either maximum or minimum), then
the generator excitation system would not be able to maintain the terminal voltage
magnitude at the specified value. In that case the generator bus would behave as a PQ bus
(P being already specified and QG is set at either maximum or minimum limiting value of
Q violated). Thus, both the magnitude and angle of the bus voltage are calculated in the
present iteration. This phenomenon (where the generator is behaving like a PQ bus) is
termed as ‘PV to PQ switching’.
Complete GSLF algorithm
Step 1: Initialize Vj(0) = Vjsp00 for j = 2, 3, …, m and Vj(0) = 1.000 for j = (m+ 1), (m+ 2), …,
n. Set iteration count k = 1.
Step 2: For PV buses, i = 2, 3, …, m, carry out the following operations.
a) Calculate,
∗ ∗
𝐤 𝐤−𝟏 𝐤 𝐤−𝟏
𝐐𝐢 = −𝐈𝐦 𝐕𝐢 σ𝐢−𝟏 𝐘 𝐕
𝐣=𝟏 𝐢𝐣 𝐣 + 𝐕𝐢 σ𝐧𝐣=𝐢 𝐘𝐢𝐣 𝐕𝐣 𝐤−𝟏 or
𝐧
𝐤 𝐤−𝟏 𝐤−𝟏 𝐤−𝟏 𝐤−𝟏
𝐐𝐢 = 𝐕𝐢 𝐕𝐣 𝐘𝐢𝐣 𝐬𝐢𝐧 𝛉𝐢 − 𝛉𝐣 − 𝛂𝐢𝐣
𝐣=𝟏
b) If, Qimin ≤ Qi(k) ≤ Qimax ; then assign |Vi(k)| = Visp and θi(k) = (Ai(k)). The quantity Ai(k) is given
by,
𝐤 𝐢−𝟏 𝐧 𝐤
𝐤 𝟏 𝐏𝐢 − 𝐣𝐐 𝐢 𝐤 𝐤−𝟏 𝐤 𝐕 𝐢
∠𝐕𝐢 = ∠ − 𝐘 𝐕
𝐢𝐣 𝐣 − 𝐘 𝐕
𝐢𝐣 𝐣 𝐎𝐫 𝐕 = 𝐕𝐢
𝐘𝐢𝐢 𝐕 𝐤−𝟏 ∗ 𝐢 𝐜𝐨𝐫𝐫𝐞𝐜𝐭𝐞𝐝
𝐕
𝐤
𝐢 𝐣=𝟏 𝐣=𝐢+𝟏 𝐢
c) If Qi(k) ≥ Qimax, then set Qi(k) = Qimax and calculate
𝒎𝒂𝒙 𝐢−𝟏 𝐧
𝐤 𝟏 𝐏𝐢 − 𝐣𝐐𝐢 𝐤 𝐤−𝟏
𝐕𝐢 = − 𝐘 𝐕
𝐢𝐣 𝐣 − 𝐘 𝐕
𝐢𝐣 𝐣
𝐘𝐢𝐢 𝐕 𝐤−𝟏 ∗
𝐢 𝐣=𝟏 𝐣=𝐢+𝟏
d) If Qi(k) ≤ Qimin, then set Qi(k) = Qimin and calculate
𝒎𝒊𝒏 𝐢−𝟏 𝐧
𝐤 𝟏 𝐏𝐢 − 𝐣𝐐𝐢 𝐤 𝐤−𝟏
𝐕𝐢 = − 𝐘 𝐕
𝐢𝐣 𝐣 − 𝐘 𝐕
𝐢𝐣 𝐣
𝐘𝐢𝐢 𝐕 𝐤−𝟏 ∗
𝐢 𝐣=𝟏 𝐣=𝐢+𝟏
Step 3: For PQ buses, i = (m+ 1), …, n, calculate
𝒌 𝐢−𝟏 𝐧
𝐤 𝟏 𝐏𝐢 − 𝐣𝐐𝐢 𝐤 𝐤−𝟏
𝐕𝐢 = − 𝐘 𝐕
𝐢𝐣 𝐣 − 𝐘 𝐕
𝐢𝐣 𝐣
𝐘𝐢𝐢 𝐕 𝐤−𝟏 ∗
𝐢 𝐣=𝟏 𝐣=𝐢+𝟏
Step 4: Compute ei(k) = |Vi(k) − Vi(k-1) | for all i = 2,3, …, n;
Step 5: Compute e(k) = max(e2(k), e3(k), …, en(k));
Step 6: If e(k) ≤ ε, stop and print the solution. Else set k = k + 1 and go to step 2.
Ex
Ex
Ex
Ex
Basic Newton - Raphson (NR) Techniques
Let there be ‘n’ equations in ‘n’ unknown variables x1, x2, …, xn as given below,
To solve the above equations, first we take an initial guess of the solution and let the
0 0 0 0 𝑇
vector of initial guesses be denoted as 𝑥 = 𝑥1 , 𝑥2 , … , 𝑥𝑛 . Subsequently, first
order Taylor’s series expansion (neglecting the higher order terms) is carried out for these
equation around the initial guess of solution as follow:
The above Equations can be written in matrix form as,
In the above equation, the matrix containing the partial derivative terms is known as the
Jacobin matrix (J). It is an n x n square matrix.
By rearranging and simplifying the above equation yields
This is the basic equation for solving the ‘n’ algebraic equations given at first place.
The steps of solution are as follow:
Step 1: Assume a vector of initial guess x(0) and set iteration counter k = 0.
Step 2: Compute f1(x(k)), f2(x(k)), …, fn(x(k)).
Step 3: Compute ∆m1, ∆m2, …, ∆mn.
Step 4: Compute error = max [|∆m1|, |∆m2|, …, |∆mn|]
Step 5: If error ≤ ε (pre - specified tolerance), then the final solution vector is x(k) and print
the results. Otherwise go to step 6.
Step 6: Form the Jacobin matrix analytically and evaluate it at x = x(k).
Step 7: Calculate the correction vector ∆𝑥 = ∆𝑥1 , ∆𝑥2 , … , ∆𝑥𝑛 𝑇 by using the above
equation.
Step 8: Update the solution vector x(k+1) = x(k) + ∆x and update k = k + 1 and go back to
step 2.
Newton Raphson Load Flow Technique (NRLF) in polar coordinates
For NRLF techniques, the starting equations are same as the basic load flow equations:
n 𝑛
Then
𝐽1 𝐽2 ∆𝜃
Where, the matrix 𝐽 = is known as the Jacobian matrix, the vector ∆𝑥 = is
𝐽3 𝐽4 ∆𝑉
known as the correction vector and
𝑃𝑠𝑝 − 𝑃𝑐𝑎𝑙
The vector ∆𝑀 = 𝑠𝑝 𝑐𝑎𝑙 is known as the mismatch vector.
𝑄 −𝑄
The elements of the Jacobian matrix can be calculated as follows:
𝜕𝑃
Matrix 𝐽1 = (in this case, i = 2, 3, …, n, & j = 2, 3, …, n)
𝜕𝜃
𝜕𝑃
Matrix 𝐽2 = (in this case, i = 2, 3, …, n, & j = m + 1, m + 2, …, n)
𝜕𝑉
𝜕𝑄
Matrix 𝐽3 = (in this case, i = m + 1, m + 2, …, n, & j = 2, 3, …, n)
𝜕𝜃
𝜕𝑄
Matrix 𝐽4 = (in this case, i = m + 1, m + 2, …, n, & j = m + 1, m + 2, …, n)
𝜕𝑉
Where, 𝐺𝑖𝑖 = 𝑌𝑖𝑖 cos 𝛼𝑖𝑖 and 𝐵𝑖𝑖 = 𝑌𝑖𝑖 sin 𝛼𝑖𝑖
Complete NRLF (polar) algorithm
Step 1: Initialize Vj(0) = Vjsp ˂ 0o for j = 2, 3, …, m and Vj(0) ≤ 1.0∟0o for j = (m+1),
(m+2), …, n. Let the vectors of the initial voltage magnitudes and angles be denoted as
V(0) and θ(0) respectively.
Step 2: Set iteration counter k = 1.
Step 3: For i = 2, 3, …, m, carry out the following operations.
a) Calculate,
𝑛
𝑘−1 𝑘−1 𝑘−1 𝑘−1
𝑄𝑖 = 𝑉𝑖 𝑉𝑗 𝑌𝑖𝑗 𝑠𝑖𝑛 𝜃𝑖 − 𝜃𝑗 − 𝛼𝑖𝑗
𝑗=1
b) If 𝑄𝑖 𝑚𝑖𝑛 ≤ 𝑄𝑖 ≤ 𝑄𝑖 𝑚𝑎𝑥 ; then assign 𝑉𝑖 𝑘 = 𝑉𝑖 𝑠𝑝 and the ith bus retained as PV
𝑘
Similarly,
Again,
The basic FDLF equation is:
Where,
Thus, there is a decoupling between ‘∆P - ∆θ’ and ‘∆Q - ∆V’ relations (i.e., ∆P depends
only on ∆θ and ∆Q depends only on ∆V).
The elements of J1 and J4 are:
Elements of J1:
Elements of J4:
We know ∆P = J1 ∆θ
Substituting J1 in to the above equation, results:
n
∆𝑷
= 𝑩′ ∆𝛉
𝑽
Matrix B’ is a constant matrix of (n−1) × (n−1). Its elements are the negative of the
imaginary part of the element (i, k) of the YBUS matrix where i & k = 2,3, …, n.
Similarly, 𝐧
No. GS NR FD
1 Provide large no. iteration to Provide less no. iteration to converge Provide even less no. iteration
converge to converge than NR
2 The computation time per iteration The computation time per iteration is The computation time per
is less more iteration is less
3 fewer memory requirement high memory requirement few memory requirement
For inverter
Therefore, for a two terminal HVDC link, the complete set of solution vector is;
Therefore, out of 9 unknown variables, any 4 variables need to be specified and thereafter,
the remaining 5 variables can be solved using 5 DC model equations. These 4 variables
can be specified using the control specification. There can be several combinations of
control specification and some of their combination are;