View metadata, citation and similar papers at core.ac.

uk brought to you by CORE

provided by Biblioteca Digital do IPB

IEEE MELECON 2006, May 16-19, Benalmádena (Málaga), Spain

A Comparison of Voltage Stability Indices

Claudia Reis F.P. Maciel Barbosa
Polytechnic Institute of Braganca Faculty of Engineering of Porto University
Portugal Portugal
Email: claudia-reis@ipb.pt Email: fmb@fe.up.pt

Abstract— Voltage stability has become a very importante issue behaviour has been approached using static techniques, which
of power systems analysis. This paper discusses some important have been widely used on voltage stability analysis.
aspects related to voltage stability indices in electric power
An accurate knowledge of how close the actual system’s
systems. Some techniques previously studied in the literature
are analyzed and a comparison of the performance of several operating point is from the voltage stability limit is crucial
indices is presented. The effectiveness of the analyzed methods to operators. Therefore, to find a voltage stability index has
are demonstrated through numerical studies in IEEE 14 busbar become an important task for many voltage stability studies.
test system, using several different scenarios of load increase. These indices provide reliable information about proximity of
voltage instability in a power system. Usually, their values
changes between 0 (no load) and 1 (voltage collapse).
I. I NTRODUCTION
The voltage stability analysis, using different methods, will
Problems related to voltage stability in power systems are be highlighted in this paper and the results obtained from
one of the major concerns in power system planning and simulating on IEEE 14 busbar test system will be discussed.
operation.
Voltage stability is concerned with the ability of a power
system to maintain acceptable voltages at all nodes in the II. I NDICES F ORMULATION
system under normal condition and after being subject to a
disturbance [1]. The condition of voltage stability in a power system can be
A power system is said to have a situation of voltage known using voltage stability indices. This indices can either
instability when a disturbance causes a progressive and un- reveal the critical bus of a power system or the stability of
controllable decrease in voltage level. each line connected between two bus in an interconnected
During the last decades, the voltage stability problem has network or evaluate the voltage stability margins of a system.
been given more attention primarily due to a number of The indices used to examine the system stability are briefly
stability accidents that occurred in some countries. Some well- described in this section.
known examples of voltage stability incidents were reported
in France, Belgium, Sweden, Germany, Japan and in USA [2], 1) P-V and Q-V curves: The P-V curves are the most
[3]. used method of predicting voltage security. They are used to
As power systems become more complex and heavily determine the loading margin of a power system. The power
loaded, along with economical and environmental constraints, system load is gradually increased and, at each increment, is
voltage instability becomes an increasingly serious problem, necessary recompute power flows until the nose of the PV
leading systems to operate close to their limits. Voltage curve is reached. The margin between the voltage collapse
instability is essentially a local phenomenon, however its point and the current operating point is used as voltage stability
consequences may have widespread impact. criterion [5].
The study of voltage stability has been analyzed under dif- With Q-V curve is possible, for the operators, to know
ferent approaches that can be basically classified into dynamic which is the maximum reactive power that can be achieved or
and static analysis. added to the weakest bus before reaching minimum voltage
The static voltage stability methods depend mainly on the limit. The reactive power margin is the MVar distance from
steady state model in the analysis, such as power flow model the operating point to the bottom of the Q-V curve. The Q-V
or a linearized dynamic model described by the steady state curve can be used as an index for voltage instability. The point
operation. where dQ/dV is zero is the point of voltage stability limit [1].
The dynamic analysis implies the use of a model charac- 2) L index: Kessel et al. [6] developed a voltage stability
terized by nonlinear differential and algebric equations which index based on the solution of the power flow equations. The
include generators dynamics, tap changing transformers, etc, L index is a quantitative measure for the estimation of the
through transient stability simulations [4]. distance of the actual state of the system to the stability limit.
Although stability studies, in general, requires a dynamic The L index describes the stability of the complete system and
model of the power system, in this paper analysis of voltage is given by:

1-4244-0088-0/06/$20.00 ©2006 IEEE 1007

 Lines that presents values of Lmn close to 1, indicates that
. ! .
. F ji V i . those lines are closer to theirs instability points. To maintain
.
L = max {Lj } = max .1 − i∈αG . (1)
j∈αL j∈αL V . a secure condition, the Lmn index should be less than 1.
j
6) Line Stability Index FVSI: The line stability index FVSI
where αL is the set of consumer nodes and αG is the set of proposed by I.Musirin et al. [9] is based on a concept of power
generator nodes. flow through a single line. For a typical transmission line, the
Lj is a local indicator that determinates the busbars from stability index is calculated by:
where collapse may originate. The L index varies in a range
between 0 (no load) and 1 (voltage collapse). 4Z 2 Qj
F V SIij = (3)
3) V/V0 index: Assuming the bus voltage values (V) to Vi2 X
be known from load flow or state estimation studies, new where Z is the line impedance, X is the line reactance, Qj
bus voltages (V0) are obtained solving a load flow for the is the reactive power flow at the receiving end and Vi is the
system at an identical state but with all loads set to zero. sending end voltage.
The ratio V/Vo at each node yields a voltage stability map
The line that gives index value closest to 1 will be the most
of the system, allowing for immediate detection of weak and
critical line of the bus and may lead to the whole system
effective countermeasure spots [5].
instability. The calculated FVSI can also be used to determine
4) Modal analysis: Gao et al. [7] proposed a method that
the weakest bus on the system. The determination of the
computes the smallest eigenvalue and associated eigenvectors
weakest bus is based on the maximum load allowed on a load
of the reduced Jacobian matrix of the power system based on
bus. The most vulnerable bus in the system corresponds to the
the steady state system model. The eigenvalues are associated
bus with the smallest maximum permissible load.
with a mode of voltage and reactive power variation. If all
7) Line Stability Index LQP: The LQP index derived by
the eigenvalues are positive, the system is considered to be
A.Mohamed et al. [10] is obtained using the same concept
voltage stable. If one of the eigenvalues is negative, the system
as [8], [9], in which the discriminant of the power quadratic
is considered to be voltage unstable. A zero eigenvalue of the
equation is set to be greater or equal than zero. The LQP is
reduced Jacobian matrix means that the system is on the border
obtained as follows:
of voltage instability. The potential voltage collapse situation   
of a stable system can be predicted through the evaluation X X 2
LQP = 4 P + Qj (4)
of the minimum positive eigenvalues. The magnitude of each Vi2 Vi2 i
minimum eigenvalue provides a measure to know how close
where X is the line reactance, Qj is the reactive power flow
the system is to voltage collapse. By using the bus participation
at the receiving bus, Vi is the voltage on sending bus and Pi
factor, the weakest bus can be determined, which is the greatest
is the active power flow at the sending bus.
contributing factor for a system to reach voltage collapse
To maintain a secure condition, the value of LQP index
situation.
should be maintained less than 1.
5) Line Stability Index Lmn : M.Moghavemmi et al. [8]
8) Line Stability Indices VCPI: The VCPI indices proposed
derived a line stability index based on the power transmission
by M.Moghavvemi et al. [11] investigates the stability of
concept in a single line, in which discriminant of the voltage
each line of the system and they are based on the concept
quadratic equation is set to be greater or equal than zero to
of maximum power transfered through a line.
achieve stability. If the discriminant is small than zero, the
roots will be imaginary, which means that cause instability in PR
V CP I(1) = (5)
the system. Figure 1 ilustrates a single line of an intercon- PR(max)
nected network where the Lmn is derived from.
QR
V CP I(2) = (6)
i
Vi Vj
j QR(max)
Si,Pi,Qi
I
Sj, Pj, Qj where the values of PR e QR are obtained from conventional
power flows calculations, and PR(max) and QR(max) are the
R+jX maximum active and reactive power that can be transfered
through a line.
Fig. 1. Typical one-line diagram of transmission line.
The VCPI indices varies from 0 (no load condition) to 1
(voltage collapse).
The line stability index, for this model, can be defined as:
4XQj
Lmn = 2 (2)
[Vi sen (θ − δ)] III. T EST RESULTS AND DISCUSSION
where θ is the line impedance angle and δ is the angle The voltage stability analysis were performed on IEEE 14
difference between the supply voltage and the receiving end busbar test system. This system has 5 generator busbars, 9
voltage. load busbars and 20 interconnected branches.

1008
 The voltage stability margin can be easily calculated with all voltages in the busbars present values very near the ones
P-V curves. This curves show the bus voltage level as the obtained with the system with all loads a zero (V0). The
loading factor λ increases. The loading factor is 1 at base situation B reveals that the voltage at the busbars have values
case and is gradually increased, in all busbars of the system, very far away of the ones obtained with the system with all
until maximum loading point is reached, maintaining constant loads a zero. So, in situation B, the system is near of the
the power factor. voltage instability. Figure 4 shows that the critical bus of this
system is bus 14, because it presents the smallest index value.
V (p.u.)

V/V0
Bus 1
100%
Bus 2
1 90%
Bus 3
Bus 4 80%
Bus 5
70%
Bus 6
60%
Bus 7
Bus 8 50%
0,8 Bus 9 40%
Bus 10
30%
Bus 11
Bus 12 20%
Bus 13 10%
Bus 14
0%
0,6 1 3 5 7 9 11 13
Bus
1 1,1 1,2 1,3 1,4 1,5 1,6 1,7 1,75 1,779
Loading factor Situation A Situation B

Fig. 2. P-V curves for IEEE 14 busbar test system.

Fig. 4. V/V0 index for IEEE 14 bus test system.
Analysing the Figure 2 we observe that, as the power
system load is gradually increased, the voltages at the busbars To illustrate the application of modal analysis, we choosed
decreases. As shown in Figure 2, the voltage stability margin three operating conditions on the P-V curves of Figure 2:
of this system is approximated 77, 9%. the base case, the case of loading factor 1,5 and the critical
operating case (loading factor 1,779). All the eigenvalues of
The Q-V curve calculates the voltage stability margin of a
reduced reactive Jacobian matrix are positive, indicating that
particular bus, i.e., this curves show how much reactive power
the system is voltage stable at all tested load conditions.
demand can be increased before the system suffers a voltage
The magnitudes of the eigenvalues decrease as the system
collapse.
approaches to instability.
MVAr At the critical operating point, the smallest eigenvalue is
250
1,4369. This value is considered the least stable mode for
200
the critical operating point and is used to determine the bus
participation factors. Figure 5 shows the bus participation
150 factors in the least stable mode for the critical operating point.
100
0,4
50
Bus participations

0,3
0
4 5 7 9 10 11 12 13 14
0,2
Bus
0,1
Fig. 3. Margin of reactive power in IEEE 14 bus test system.
0
In Figure 3, busbars 12 and 14 have the lowest margin of 1 2 3 4 5 6 7 8 9 10 11 12 13 14

reactive power, which indicates that they are the most critical Bus

busbars in the system.

Fig. 5. Bus participation factors in the least stable mode for critical operating
The V/Vo index is very simple to compute. In Figure 4, is case.
possible to observe the index value for two different loading
situations: for a condition near the base case (loading factor It is possible to know, with modal analysis, which bus of the
λ = 1, 2) - situation A - and for the critical operating case tested system contribute more to the voltage collapse. Figure
(loading factor λ = 1, 779) - situation B. 5 shows that the critical bus of this system is bus 14 because
As was expected, in situation A, the system is far away of busbars with high participation factors are considered the most
the voltage instability. This fact can be easily proved because vulnerables.

1009
 The local index Lj permits the determination of the weakest This table also illustrates that, for any state of loading,
bus in the system. Figure 6 shows the values of the local index VCPI(1) is always equal to VCPI(2).
Lj in the IEEE 14 busbar test system. Table 1 shows that the performance of the line indices
studied has high degree of accuracy, reliability and the results
Lj are very closed in agreement. Line stability indices can also
1
Bus 4
determine the weakest bus in the system and it is based on the
0,8 Bus 5
maximum permissible load. It is observed that buses 10, 11
and 14 indicates 0,948 p.u., 0,855 p.u. and 0,728 p.u. as the
Bus 7
0,6
maximum permissible reactive load, respectively. Since bus
Bus 9
14 has the smallest maximum loadability, it is considered the
Bus 10
0,4 most critical bus.
Bus 11

Bus 12 Table 1 - Line stability indices for reactive power load

0,2
Load (p.u.) Line Lmn FVSI LQP VCPI(1) VCPI(2)
Bus 13
Q10=0,948 9-10 0,5238 0,5405 0,4741 0,495 0,495
0 Bus 14 10-11 0,5213 0,5346 0,4557 0,4947 0,4947
1 1,1 1,2 1,3 1,4 1,5 1,6 1,7 1,779
Loading factor
Q11=0,855 6-11 0,7723 0,8272 0,681 0,7311 0,7311
10-11 0,6265 0,6673 0,5692 0,5969 0,5969
Q14=0,728 9-14 0,8475 0,8994 0,756 0,8184 0,8184
Fig. 6. Evaluation of Lj index versus load variation. 13-14 0,801 0,8493 0,7 0,7668 0,7668

Figure 6 shows that bus 14 exhibits the highest Lj index,

which indicates that it is the most vulnerable bus on the IV. C ONCLUSION
system.
This paper presents a comparative study and analysis of the
In Figure 7, L index and the voltage at bus 14 (critical bus)
performance of some online static voltage collapse indices.
are plotted as a function of loading factor.
The application of those indices on IEEE 14 bus testing system
gave accurate results. The shown simulations indicate that
the bus 14 of IEEE 14 busbar test system is considered the
weakiest bus in the system. Line indices provide an accurate
information with regard to the stability condition of the lines.
The research shows an agreement between the different line
stability indices.
R EFERENCES
[1] P.Kundur “Power System Stability and Control” McGraw-Hill, New
York, 1994
[2] “Voltage Stability of Power Systems: Concepts, Analytical Tools and
Industry Experience”, IEEE Committee Vol.IEEE/PES 93TH0358-2-
Fig. 7. Stability indicator L and its relation to the critical voltage. PWR 1990.
[3] K. Takahashi and Y. Nomura “The Power System Failure on July 23rd
In the critical operating point L=0,958, so the voltage 1987 in Tokyo” CIGRE SC-37 Meeting 37.87(JP)07(E) 1987.
[4] J.C. Chow, R. Fischl and H. Yan “On the Evaluation of Voltage
stability of this system is guaranteed. The stability limit is Collapse Criteria” IEEE Trans., PWRS-5, pp. 612-620, May 1990
reached for L=1. [5] Editor/Coordinator: Claudio Canizares “Voltage Stability Assessment:
Concepts, Practices and Tools” IEEE/PES Power System Stability
Subcommittee Special Publication, August 2002
To investigate the effectiveness of the line stability indices, [6] P.Kessel, H.Glavitsch “Estimating the Voltage Stability of a Power
the reactive load were gradually increased, only in one bus at System” IEEE, Transactions on Power Delivery, Vol.PWRD-1, N3, July
a time, from the base case until its maximum allowable load, 1986
[7] B.Gao, G.K.Morison, P.Kundur “Voltage Stability Evaluation Using
keeping the load at the other busbars fixed at base load. The Modal Analysis” IEEE, Transactions on Power Systems, Vol.7, N4,
test was carried out for several buses but only the cases of November 1992
buses 10, 11 and 14 are presented. Table 1 shows the stressed [8] M.Moghavvemi, F.M.Omar “Technique for Contingency Monitoring
and Voltage Collapse Prediction” IEEE Proceeding on Generation,
conditions of the lines for the maximum loadability of the Transmission and Distribution, Vol. 145, N6, pp. 634-640 November 1998
buses. The values of the line stability indices are maximum [9] I.Musirin, T.K.A.Rahman “Novel Fast Voltage Stability Index (FVSI)
at the maximum loadability of each individual load bus. The for Voltage Stability Analysis in Power Transmission System” 2002
Student Conference on Research and Development Proceedings, Shah
line that presents the largest index with respect to a bus is Alam, Malasia, July 2002
considered the most critical line of that bus. From Table 1 it [10] A.Mohamed, G.B.Jasmon, S.Yusoff “A Static Voltage Collapse Indi-
is observed that the line that connect bus 9 to bus 10 is the cator using Line Stability Factors” Journal of Industrial Technology,
Vol.7, N1, pp. 73-85, 1989
most critical line referred to bus 10, the line 6-11 is the most [11] M.Moghavvemi, O.Faruque “Real-Time Contingency Evaluation and
critical line with respect to bus 11 and line 9-14 is the most Ranking Technique” IEEE Proceeding on Generation, Transmission and
critical line of bus 14. Distribution, Vol. 145, N5, September 1998

1010