Next Article in Journal
Power Output Enhancement of Straight-Bladed Vertical-Axis Wind Turbines with Surrounding Structures
Next Article in Special Issue
Real-Time Control of a Battery Energy Storage System Using a Reconfigurable Synchrophasor-Based Control System
Previous Article in Journal
Green Production of Biodiesel from High Acid Value Oil via Glycerol Esterification and Transesterification Catalyzed by Nano Hydrated Eggshell-Derived CaO
Previous Article in Special Issue
Long-Term Voltage Stability Bifurcation Analysis and Control Considering OLTC Adjustment and Photovoltaic Power Station
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

A Comprehensive Review on the Modelling and Significance of Stability Indices in Power System Instability Problems

by
Chandu Valuva
1,†,
Subramani Chinnamuthu
1,*,
Tahir Khurshaid
2,† and
Ki-Chai Kim
2,*
1
Department of Electrical and Electronics Engineering, SRM Institute of Science and Technology, Kattankulathur 603203, India
2
Department of Electrical Engineering, Yeungnam University, Gyeongsan 38541, Republic of Korea
*
Authors to whom correspondence should be addressed.
These authors equally contributed to this work.
Energies 2023, 16(18), 6718; https://doi.org/10.3390/en16186718
Submission received: 28 July 2023 / Revised: 6 September 2023 / Accepted: 14 September 2023 / Published: 20 September 2023
(This article belongs to the Special Issue Advanced Electric Power System 2023)

Abstract

:
Many technological advancements in the modern era have made actual use of electrical power and the constrained operating of power systems within stability limits. Some expeditious load variations and rising power demands initiate complications in voltage stability and can put stress on performance, leading to voltage instability. Voltage Stability Indices can be used to perform voltage stability assessment. This review evaluates various VSIs based on mathematical derivations, assumptions, critical values, and methodology. VSIs determine the maximum loadability, voltage collapse proximity, stability margin, weak areas, and contingency ranking. Stability indices can also specify the optimal placing and sizing of Distributed Generators. Thus, VSIs play a vital role in power system voltage stability. This review is a comprehensive survey of various indices and analyses their accuracy in determining the instability of power systems. Voltage stability is a crucial concern in operating a reliable power system, and the systematic evaluation of voltage stability is essential in a power system. This review considered and analyzed 34 indices from 138 articles from the literature for their significant performance in various power system stability problems. Of 33 indices, were 22 derived from transmission line parameters, referred to as line indices, and 12 from bus and line parameters, referred to as bus indices.

1. Introduction

Voltage instability is a crucial phenomenon that effects power systems because it is the main reason for blackouts and voltage collapse. Some unexpected consequences like an unexpected increase in load, line outage, and generator tripping may provoke an excessive load demand in the power system, leading to voltage instability. Consequently, a power system may be disrupted without appropriate action, resulting in cascading failures and large-scale blackouts. The notable blackouts around the world have primarily been caused by voltage instability. From [1], From 1965 to 2005, it can be seen that 22 significant backouts occurred; out of these 22 blackouts, 11 had voltage instability as a leading cause. Voltage instability is also known as voltage collapse; moreover, when instability is developing, the continuity of such events affects the system, resulting in low voltage and even more power outages [2].
According to an IEEE Power System Engineering Committee, “Voltage stability is the ability of a system to maintain voltage so that, when increasing load admittance, load power increases, and so that both power & voltage are manageable”. The consolidation definition of the IEEE and CIGRE is: “Voltage stability refers to the ability of a power system to maintain steady voltages at all buses in the system after being subjected to a disturbance from a given initial operating condition”.
Power system voltage instability may be due to the inadequacy of the source to deliver sufficient reactive power or to faulty power transmission lines not transmitting the desired reactive power. Voltage instability may affect the stability of a wide area or local area, resulting in the voltage collapse of a whole power system—mainly reliable active power delivered from the generators and reactive compensators like shunt capacitors. Voltage instability in a network is most probably the cause of unexpected load demand, transformer tripping, being over the limit of the generators’ reactive power, or defects in the on-load tap changing transformers. Consequently, the regular tracking and forecasting of instability are essential in the power system. In this respect, definite methods for suppressing voltage instability are:
  • Strengthen the weak bus while planning the power system, like accurately selecting distributed generators and balancing voltage.
  • Temporarily defending the supply by load shedding or improving the power factor by shunt capacitor switching.
  • Enhancing the Voltage Stability Margin (VSM) with FACTS devices.
  • Blocking on-load tap changers (OLTCs) to operate the transformer.
  • Build the generation station near to the load centre and reduce the transmission line length.
  • For executive controllers, make use of system voltage fluctuation controllers, Shunt Compensation, step-up transformer controllers, and Automatic Voltage Regulators (AVRs).
  • For Real-time Performance: Scheduling the generation according to load demand, evaluating voltage stability, and protecting the load shedding.
  • For safeguarding systems: Instant Load Tap Changer (LTC) control, Contingency of load demand, and High Voltage Direct Current (HVDC).
Indeed, over the past few decades, the reliability and stability of power systems have gone through considerable improvements. Stability analysis has evolved from simple static and transient stability analyses to more sophisticated and comprehensive methods.
Voltage Stability Indices (VSIs) generally measure voltage stability. The stability indices’ performance is essential in analyzing and illustrating the consequences of power system operation. Furthermore, the indices assist in predicting eventual improvement and determining long-term advancement. Especially in the last three decades, numerous techniques and methodologies have been developed for easy stability analysis. The performances of various indices for determining voltage stability have been explained in this review. For assessing a system’s stability, several indices consider the system impedance. However, few indices are unconstrained by impedance and performed using the current and voltage of the system. In real-time operation, it is impossible to determine the accurate impedance of a system because of troposphere conditions and inadequate knowledge of the system. Therefore, the accuracy of the impedance-operated indices is repeatedly accompanied by errors.
A taxonomy of comprehensive characteristics, variations in accuracy, and some parameters being inaccurate is presented below, as well as the terminology concerning VSIs, including mathematical derivation, instability conditions, and assumptions of specific parameters. This taxonomy can serve to make researchers aware of the advancements in the VSIs, and the respective operations including voltage stability enhancement, strengthening the weak bus, locating the reactive compensation devices, reducing the cost-effectiveness and location of FACTS devices, and counteracting the prevention of voltage collapse.
Determining the VSI values for a system is carried out in three modes of operation. One is the online mode: In this mode of operation, the engineer observes the voltage stability and operates accordingly to maintain the system’s stability. The statistics used for this in this paper were composed of mathematical models and power system measurements. The second mode of the stability indices is the offline mode: In this mode of operation, the designer gathers the simulated or historical data. The preserved data identify the voltage instability intensity and voltage collapse and do not contain real-time data on voltage instability. The real-time method is the third mode of stability indices, and is obtained using a Phasor Measurement Unit, Power Flow Analysis, and Artificial Intelligence.
This article aims to thoroughly review and scrutinize numerous Voltage Stability Indices. This review indicates the limitations, strengths, underlying principles, and applicability of these indices by sharply examining them. This article provides clarity to researchers, engineers, and power system operators in choosing suitable stability indices depending on the particular characteristics of a system and functional needs. Researchers will gain synthesized knowledge of Voltage Stability Indices, as this article presents the diverse range of stability indices with mathematical formulations, highlighting their peculiar characteristics and underlying principles. This review supports decision-making for system operators with the knowledge needed to estimate instabilities and ensure voltage stability, sustainable energy supply, and grid reliability.

2. Voltage Stability Indices

Voltage Stability Indices are values of an arithmetic sequence adopted for evaluating voltage stability. The indices specify the accessible and equitable regulation of the system’s stability so that voltage levels are within sustainable limits. Indices predominantly assist power engineers in determining voltage stability, predicting voltage collapse, and evaluating the significance of system characteristics. This review elaborates on the numerous methodologies for assessing voltage stability. Furthermore, categorizing the different VSIs in [3,4,5] as line VSIs, Bus VSIs, and overall VSIs. Jacobian matrix-based indices and PMU-based indices fall under overall VSIs.
The overall VSIs’ precision is superior compared to the other stability indices. The overall VSIs are challenging to determine and take more computation time. Line VSIs obtain easy and simple indices that are reasonable in determining all optimization issues through the proficiency of these indices and identifying the critical line and weak bus. These VSIs are derived from a mathematical equation by considering the standard two-bus concept. With these mathematical derivations, the system examines the stability of various stability indices. Similar aspects of voltage collapse points determine the VSIs’ formation. Figure 1 represents an illustration of the interconnected network with a single line in a standard two-bus model.
Jacobian matrix-based VSIs specify the scope of system stability, determine the voltage collapse point, and evaluate the VSM. However, they consume more time for computation, and in general, mathematical topology variations compel the transformation of the Jacobian matrix, causing the matrix to need to be again; therefore, they are inconvenient in real-time applications. Additionally, Jacobian matrix indices tend to increase the operational time of DG placing and sizing issues. However, the VSIs established by system variables claim lesser computational time and are acceptable for real-time operations. The detriment of these stability indices is that they are less accurate in evaluating the VSM. Consequently, they are not as accurate in determining the weak line and critical bus. In most applications, like placement of the DG and sizing issues, VSIs prevail for identifying the weak bus and critical line or stimulating the deterrents for voltage instability.
S r = P r + j Q r = V r I *
I = S r V r * = P r j Q r V r δ r
Vs, Vr: magnitude of bus voltage at sending and receiving bus; Ss, Ps, Qs: apparent, real, and reactive power at sending bus; Sr, Pr, Qr: apparent, real, and reactive power at receiving bus; δs, δr: voltage angle at sending and receiving bus; R, X, Z, θ = line resistance, line reactance, line impedance, and line impedance angle.
Considering the Kirchhoff Voltage Law (KVL), the current equation for two buses is,
I = V s δ s V r δ r R + j X
Substituting Equation (3) in Equation (2), we obtain:
V s δ s V r δ r R + j X = P r j Q r V r δ r
V s V r ( δ s δ r ) V r 2 = P r R j Q r R + j P r X + Q r X
Let δsδr = δ and simplify Equation (5) as:
V s V r δ V r 2 = P r R j Q r R + j P r X + Q r X
Transforming VsVr δ into a rectangular form, we obtain:
V s V r cos δ + j V s V r sin δ V r 2 = P r R j Q r R + j P r X + Q r X
Segmenting Equation (7) into real and imaginary parts, we obtain:
V s V r cos δ V r 2 = P r R + Q r X
V s V r sin δ = Q r R + P r X
Readjusting Equations (8) and (9), we obtain:
P r = Q r X + V s V r cos δ V r 2 R
Q r = P r X V s V r sin δ R
Substituting Equation (11) in Equation (10) and vice versa, we obtain:
P r = P r X V s V r sin δ R X + V s V r cos δ V r 2 R
Q r = Q r X + V s V r cos δ V r 2 R X V s V r sin δ R
Rearrange Equations (12) and (13) allows us to obtain:
V r 2 V s V r cos δ + X R sin δ + P r R + X 2 R = 0
V r 2 + V s V r R X sin δ cos δ + Q r R 2 X + X = 0
For a few VSIs, the line resistance (R) and reactance (X) increase with the impedance (Z) and impedance angle (θ); therefore, reducing the current equation, Equation (3) is redefined as:
I = V s δ s V r δ r Z θ
Substituting Equation (16) into Equation (1), the receiving end’s apparent power concerning the real and reactive power is:
P r + j Q r = V r V s δ s V r δ r Z θ *
Assuming δ s δ r = δ then Equation (17) is simplified as:
P r + j Q r = V s V r Z ( θ δ ) V r 2 Z θ
Equating the real and imaginary parts of Equation (18), we obtain:
P r = V s V r Z cos ( θ δ ) V r 2 Z cos θ
Q r = V s V r Z sin ( θ δ ) V r 2 Z sin θ
Simplifying Equations (19) and (20), we obtain:
V r 2 c o s θ V s V r cos ( θ δ ) + P r Z = 0
V r 2 s i n θ V s V r sin ( θ δ ) + Q r Z = 0
Most of the stability indices were derived mathematically; the quadratic equation’s voltage discriminant is greater than or equal to zero. The characteristics vary for the index, so different assumptions are considered for deriving the stability indices.

2.1. Line Voltage Stability Indices

2.1.1. Voltage Stability Load Index (VLSI)

The Voltage Stability Load Index (VLSI) derives from the mathematical derivation [6] of a standard two-bus system and Thevenin equivalent circuit. Figure 1 presents the voltage equation considered in this index, obtained from the network. Moreover, we determine the power flow by performing a power flow analysis directed toward speedy converging. Considering the square of the sending-end current, receiving-end active and reactive power with losses, we obtain:
I s 2 = P s 2 + Q s 2 V s 2
P r = P s P l o s s
Q r = Q s Q l o s s
P l o s s = P r 2 + Q r 2 V r 2 R s
Q l o s s = P r 2 + Q r 2 V r 2 X s
We then substitute Equations (24)–(27) in Equation (23) to simplify it, resulting in:
I 2 = P r + P r 2 + Q r 2 V r 2 R s 2 + Q r + P r 2 + Q r 2 V r 2 X s 2 V s 2
The analyzed equation is further reduced by equating it to zero and obtaining the sending end voltage equation:
V s 2 = V r 2 + 2 P r R s + Q r X s + P r 2 + Q r 2 V r 2 R s 2 + X s 2
V r 4 + V s 2 2 P r R s + Q r X s V s 2 + P r 2 + Q r 2 R s 2 + X s 2 = 0
If the obtained receiving end voltage quadratic equation contains real roots, then the discriminant is greater than or equal to zero; hence:
8 P r Q r R s X s 4 V r 2 P r R s + Q r X s + V s 4 4 P r 2 X s 2 + Q r 2 R s 2 0
4 V s 2 P r R s + Q r X s + P r X s Q r R s 2 V s 4 1
Accordingly, L represents the voltage stability load index,
L = 4 V s 2 P r R s + Q r X s + P r X s Q r R s 2 V s 4
since V s V r c o s δ s δ r V r 2 = P r R s + Q r X s   &   V s V r s i n δ s δ r = P r X s Q r R s
Therefore   L = 4 V s V r cos δ s δ r V r 2 cos δ s δ r 2 V s 2
L = 4 V s V r cos δ V r 2 cos δ 2 V s 2
Hence, the value L should be less than 1.0 for the system’s stability. Whenever the index value is more significant than 1.0, corresponding to Equation (32), the voltage emerges as imaginary and collapses. To defend the load instability, the load index value must be Li ≤ 0, and if Li is nearer to 1.0, then bus ‘i’ approaches the stability limit.

2.1.2. Line Stability Index (Lmn)

The Line Stability Index is determined in [7] and derived from the voltage quadratic equations discriminant. Considering this index and Equation (20), we can derive Vr:
V r = V s sin θ δ ± V s sin θ δ 2 4 Z Q r s i n θ 0.5 2 s i n θ
Assuming Zsinθ = x,
V r = V s sin θ δ ± V s sin θ δ 2 4 x Q r 0.5 2 s i n θ
The above voltage equation should have real roots to obtain absolute values of Vr in terms of Qr. Therefore, the following condition is to manipulate the stability of the system.
V s s i n ( θ δ ) 2 4 x Q r 0 4 x Q r V s s i n ( θ δ 2 = L m n 1
Lmn represents the stability index of the line. Voltage collapse is predicted rigorously depending on the stability index’s value. If the obtained value of Lmn is lesser than the unity, then the system is in a stable condition; otherwise, if the value is more significant than the unity, then the system exhausts the stability; moreover, the voltage collapses.

2.1.3. Line Stability Factor (LQP)

The Line Stability Factor (LQP) developed in [8] makes use of a similar theory of the discriminant of the voltage quadratic equation, primarily from receiving end real power:
P r = P s R P s 2 Q s 2 V s 2
Reorganizing the power equation as R P s 2 V s 2 P s + P r + P r + R Q s 2 V s 2 = 0 allows us to obtain the roots of Ps:
1 4 R V s 2 P r + R Q s 2 V s 2 0
From Equation (40), it is clear that the system loses its stability when it violates this condition. Further, deriving the reactive power is similar to the derivation of active power:
1 4 X V s 2 Q r + X P s 2 V s 2 0
LQP = 4 X V s 2 Q r + X P s 2 V s 2
Since LQP < 1, the system is stable. Here, for deriving the index, it is assumed that the lines are lossless (R/X << 1) and neglect the shunt admittance.

2.1.4. Voltage Collapse Proximity Index (VCPI)

The VCPI [9] evaluates the line voltage stability depending on the theory of maximum power. The analysis assumes a frequent occurrence and varying of the absolute impedance value by keeping ϕ constant. With this inference, the performance may not be accurate, but it is easy to solve. During the operation, the power factor remained constant.
Along with the load demand increase, the current increases with a decrease in Zr, resulting in a receiving-end voltage drop:
I = V s Z s c o s θ + Z r c o s 2 + Z s s i n θ + Z r s i n 2
V r = Z r I = V s 1 + Z r / Z s 2 + 2 Z r / Z s c o s ( θ )
Hence, the receiving-end power is Pr = VrIcosϕ and Qr = VrIsinϕ
Simplifying Equation (44) with that for power, we obtain:
P r = V s 2 / Z s 1 + Z r / Z s 2 + 2 Z r / Z s c o s ( θ ) Z r Z s c o s
Q r = V s 2 / Z s 1 + Z r / Z s 2 + 2 Z r / Z s c o s ( θ ) Z r Z s s i n
Correspondingly, the power loss is P l = I 2 Z s c o s θ and Q l = I 2 Z s s i n θ .
Substituting the current equation, we then obtain:
P l = V s 2 / Z s 1 + Z r / Z s 2 + 2 Z r / Z s c o s ( θ ) c o s θ
Q l = V s 2 / Z s 1 + Z r / Z s 2 + 2 Z r / Z s c o s ( θ ) s i n θ
Considering the boundary condition P r / Z r = 0 allows us to calculate the maximum active power delivered at the receiving end. Hence, this condition results in the impedance Zr/Zs = 1. The maximum power transfer obtained by substituting this ratio in Equation (45) is:
P r ( m a x ) = V s 2 Z s c o s 4 cos 2 θ 2
Similarly, adopting this approach, the following power equation can be derived:
Maximum   transferable   reactive   power   Q r ( m a x ) = V s 2 Z s s i n 4 cos 2 θ 2
In   line   with   the   maximum   active   power   loss   P l ( m a x ) = V s 2 Z s c o s θ 4 cos 2 θ 2
In   line   with   the   maximum   reactive   loss   Q l ( m a x ) = V s 2 Z s s i n θ 4 cos 2 θ 2
Considering these maximum limits, Voltage Collapse Proximity Indicators can be designed as:
V C P I   ( 1 ) = P r P r ( m a x ) = R e a l   p o w e r   t r a n f e r r e d   t o   t h e   r e c e i v i n g   e n d M a x i m u m   r e a l   p o w e r   t h a t   c a n   b e   t r a n s f e r r e d
V C P I   ( 2 ) = Q r Q r ( m a x ) = R e a c t i v e   p o w e r   t r a n f e r r e d   t o   t h e   r e c e i v i n g   e n d M a x i m u m   r e a c t i v e   p o w e r   t h a t   c a n   b e   t r a n s f e r r e d
V C P I   ( 3 ) = P l P l ( m a x ) = R e a l   p o w e r i n   t h e   l i n e M a x i m u m   p o s s i b l e   r e a l   p o w e r   t h a t   l o s s
V C P I   ( 4 ) = Q r Q r ( m a x ) = R e a c t i v e   p o w e r   l o s s   i n   t h e   l i n e M a x i m u m   p o s s i b l e   r e a c t i v e   p o w e r   l o s s
The system may collapse if the values of the VCPIs are more significant than the unity. Therefore, VCPI (1) = VCPI (2) and VCPI (3) = VCPI (4), since only active or reactive terms were considered instead of all proximity indices.

2.1.5. Voltage Stability Index (Lp)

The line stability index proposed in [10] is a similar mathematical formulation to earlier line VSIs. Basically, from the active and reactive power equations, i.e., from Equations (19) and (20), δ can be derived as:
δ = cos 1 V s 2 P s R + Q s X V s V r
Reorganizing Equation (19) and deriving the voltage equation Vr, we obtain:
V r = V s cos ( θ δ ) ± V s cos ( θ δ ) 2 4 R P r 0.5 2 c o s θ
Hence, the Vr equation should possess real roots for obtaining the absolute values of Vr expressed in terms of Pr. Accordingly, it is then simplified to be enough to withstand the following conditions:
V s cos ( θ δ ) 2 4 R P r   0
4 R P r V s cos ( θ δ ) 2 = L p     1
LP is a line stability index that determines the condition of the transmission line and indicates instability limits. If the calculated value of Lp is more significant to the unity, the system approaches instability. The voltage collapses when the system exceeds the critical limits; moreover, the Vr becomes imaginary.

2.1.6. Fast Voltage Stability Index (FVSI)

The Fast Voltage Stability Index is derived in [11] primarily by considering the current through the line and then calculating the absolute roots of a receiving-end voltage (Equation (15)):
V r = R X s i n δ + c o s δ V s ± R X s i n δ + c o s δ V s 2 4 X + R 2 X Q r 2
For determining the absolute roots of Vr, the Vr discriminant is greater than or equal to zero.
4 Z 2 Q r X V s 2 R s i n δ + X c o s δ 2 1
The angular difference δ is minute; thus, δ 0 , Rsin δ   0, and Xcos δ  ≈ X. As a result, Equation (62) can be simplified as:
F V S I = 4 Z 2 Q r V s 2 X
The line is in-transit to the instability limit if the FVSI value is adjacent to 1. With further increases in a stability index value, that line may encounter an unexpected voltage drop accompanying system collapse.

2.1.7. Voltage Stability- Load Bus Index (VSLBI)

VSLBI evaluates the voltage stability by adopting the PMU [12]. Considering the maximum power condition, this stability index is subject to a voltage-drop ΔVr over the transmission impedance Zr that is equivalent to load bus voltage Vr:
ΔVr = Vr
Accordingly, to approach the voltage collapse exposure during the constant power loads, the Voltage Stability Load Bus Index (VLSBI) is represented as:
V S L B I = V r V
If the obtained value of the VSLBI is more significant than the unity, then the system is considered stable; if the VSLBI is less than 1, then the system is unstable, and the system may collapse.

2.1.8. Voltage Stability Margin Index (VSMI)

The VSMI, proposed in [13], depends on the correlation between the power transfer from the line and the angular difference between the sending- and receiving-end buses. To derive the fundamental equations of real and reactive power, the receiving end voltage is:
V r = V s   c o s ( r + δ ) c o s ( r )
P r = 1 2 V s 2 X s i n ( r + 2 δ ) c o s ( r ) t a n ( r )
Since t a n ( r ) = Q r P r , considering Equations (66) and (77), the power Pr maximum value can be calculated for whatever the value of Qr/Pr is. The respective angular differences of the sending- and receiving-end buses (δ), and the receiving-end voltage (Vr), are expressed as:
P r m a x = 1 2 V s 2 X 1 c o s ( r ) t a n ( r )
V r m a x = c o s π 2 + r c o s r V s
δ r m a x = π 2 r 2 for any ratio of Qr/Pr
With the correlation between the voltage stability and angular difference, the voltage stability margin that is determined depends on what degree of δ is adjacent to δrmax. Therefore, whatever the operating constraints, the receiving end index can be represented as:
V S M I = δ r m a x δ δ r m a x
The calculated value of the VSMI should be higher than zero to sustain the system’s stability and protection.

2.1.9. Voltage Collapse Proximity Index (VCPI_1)

The significant motive of VCPI_1 is that, during the voltage collapse, the voltage drop at the Thevenin impedance equals the load voltage. The proximity index derived in [14] is denoted as VCPI_1 to prevent confusion from the VCPI notation. An easy arithmetic prediction indicates that, during the critical situation, the generator phasor voltage (   V s   ¯ ) is double the load phasor voltage, V r ¯ . Hence, for the estimation of the possibility of voltage instability, VCPI_1 is represented as:
V C P I _ 1 = V r c o s δ 0.5 V s
If the VCPI is greater than or equal to zero, then the system is stable, otherwise it is unstable. In this index, the lines are complicated. However, considering that the sending-end bus connects to the ideal voltage source, the receiving end with the equivalent Thevenin impedance is neglected.

2.1.10. Critical Voltage (Vcr)

The Critical Voltage is an intelligible index that is derived from a single load, an infinite bus system utilizing the Eigen-value theorem, and load flow equations [15]. In this paper, it is used for adapting the power factor cosϕ to be constant. From the active and reactive power flow equations, the bus voltage angle δ is eliminated, and the resulting equation is:
V r 4 + 2 Z ( V r c o s θ + Q r s i n θ ) V s 2 V r 2 + P r 2 + Q r 2 Z 2 = 0
Corresponding to favorable circumstances, Equation (72) contains a specific solution; moreover, the correlation between the V s ¯ , V r ¯ ,   I ¯ vectors is developed as:
V c r = V + = V = V s 2 c o s θ
The load bus critical voltage is denoted by Vcr (P-V Curve nose point). V+ and V are the P-V curves’ upper and lower parts, respectively. Nearer to this, V+ and V are superposed values of Vcr. The indices concern the maximum power limit as a stability limiting point. Evaluating the simple system by the constant MVA load and cosϕ power factor, Equation (72) is:
S 2 + 2 S V 2 Z c o s ( θ ) + V 4 V s 2 V r 2 Z 2 = 0
Apparent   Power   S = P c o s
To obtain the maximum S, dS/dV is equated to zero. Accordingly:
V c r = V s 2 1 + c o s θ
After the phasor relationship, further reducing Equation (76) allows us to obtain:
V c r = E 2 c o s θ

2.1.11. Power Transfer Stability Index (PTSI)

The Power Stability Index proposed in [16] is derived with the help of a standard two-bus Thevenin-equivalent system; here, the slack bus and the load bus were connected on a single branch, as represented in Figure 2. The current delivery to load and load power is:
I ¯ = V ¯ T h e v Z ¯ T h e v + Z ¯ L
S ¯ L = Z ¯ L I ¯ I * ¯   = Z ¯ L I ¯ 2
Substituting Equation (78) in (79), we obtain:
S ¯ L = Z ¯ L V ¯ T h e v Z ¯ T h e v + Z ¯ L 2
Assuming that Z ¯ L = Z L φ and Z ¯ T h e v = Z T h e v θ , and alternating them into Equation (80), we obtain:
S ¯ L = Z L φ V T h e v Z T h e v θ + Z L φ 2
S L = V T h e v 2 Z L Z T h e v θ + Z L φ 2
S L = V T h e v 2 Z L Z T h e v 2 + Z L 2 + 2 Z T h e v Z L cos θ φ
To obtain the maximum value of the apparent power of the load, the SL is differentiated and equated to zero, S L / Z L = 0 .
S L m a x = V T h e v 2 2 Z T h e v 1 + 2 cos θ φ
The power margin at the load bus distance when it is approaching voltage collapse is SLmax – SL. If ZL = ZThev,, the power margin is equal to zero; moreover, it indicates no power transfer and may cause voltage collapse. Conversely, a voltage may collapse if the ratio of power is:
S L S L m a x = 1
Substituting Equations (83) and (84) in Equation (85), the power transfer stability index is:
P T S I = 2 S L Z   1 + c o s θ φ V s 2
The value of the PTSI is determined with Equation (86), considering the load’s power, impedance, Thevenin Voltage, Thevenin Impedances, and the Angular phase of the load. The obtained PTSI value lies between 0 and 1. Voltage collapse occurs when the index value reaches 1.

2.1.12. Voltage Stability Index (VSI_1)

The stability index predicts the voltage stability at the steady state operation; this index is proposed in [17]; however, for easy identification, the notation of the stability index is represented as VSI_1. Initially, this index calculates the maximum active, reactive, and apparent power transfer. To obtain the receiving-end voltage (Vr), considering the Equation (19) for active power (Pr) and Equation (20) for reactive power (Qr), we obtain:
V r = V s 2 2 Q X + P R ± A
Here, = V s 2 4 ( Q X + P R ) V s 2 ( P X Q R ) 2 , A ≥ 0.
Transmission lines containing a higher X/R ratio and ignoring the resistance value, resulting in Pmax, Qmax, and Smax, can be derived as:
P m a x = Q R X V S 2 R 2 X 2 + Z L V s V S   2 4 Q X 2 X 2
P m a x = V s 4 4 X 2 Q V s 2 X
Q m a x = P X R V S 2 X 2 R 2 + Z L V s V S 2 4 P R 2 R 2
Q m a x = V s 2 4 X P 2 X V s 2
S m a x = V s 2 Z L sin θ X + cos θ R 2 ( cos θ X s i n ( θ ) R ) 2
S m a x = 1 sin ( θ )   V s 2 2 c o s θ 2 X
The derived maximum powers identified that the maximum power transfer increases with an increase in the voltage magnitude (Vr) or a decrease in the impedance (ZL). Equivalently, Pmax decreases with an increase in Q, and Qmax decreases with an increase in P.
P m a r g i n = P m a x P , Q m a r g i n = Q m a x Q , S m a r g i n = S m a x S
From the maximum powers, the corresponding load margins can be obtained with Equation (94). Hence, the derived VSI from the obtained load margin is represented in Equation (95). If the obtained VSI is small, the load bus is adjacent to the stability margin and the lesser load margin that is available. The VSI equals zero if the load bus reaches its stability margin point.
V S I = m i n P m a r g i n P m a x , Q m a r g i n Q m a x , S m a r g i n S m a x

2.1.13. Novel Line Stability Index (NLSI)

The NLSI is determined with a similar power transmission approach [18], considering the quadratic Equations (8) and (9):
V r = V s c o s δ   ± V s 2 c o s 2 δ 4 ( P r R + Q r X ) 2
The discriminant should be greater than or equal to 0 to obtain the absolute value of Vr.
P R + Q X 0.25 V s 2 c o s 2 δ 1
As the angular difference (δ) between receiving and sending ends is usually very small, cosδ ≈ 1. Subsequently, the Novel Line Stability Index can be represented as:
N L S I = R s r P r + X s r Q r 0.25 V s 2
If the value of the NLSI of any line is adjacent to 1, that line is moving closer to the stability limit; consequently, the system may collapse. Hence, the NLSI allows us to find the stability limit, i.e., the index value is lesser than the unity.

2.1.14. Stability Index (SI)

For load flow analysis, the quadratic equation is preferable for calculating the sending-end voltage of a line and can be simplified in general form as:
V r 4 + 2 V r 2 P R + Q X V s 2 V r 2 + P 2 + Q 2 Z 2 = 0
Further, the receiving-end line active and reactive power is derived as:
P = c o s θ V r 2 ± c o s 2 θ V r 4 V r 4 Z 2 Q 2 2 V r 2 Q X + V s 2 V r 2 / Z
Q = s i n θ V r 2 ± s i n 2 θ V r 4 V r 4 Z 2 P 2 2 V r 2 P R + V s 2 V r 2 / Z
Therefore, to obtain the absolute value of the active and reactive power, the corresponding equation should satisfy the following:
c o s 2 θ V r 4 V r 4 Z 2 Q 2 2 V r 2 Q X + V s 2 V r 2 0
s i n 2 θ V r 4 V r 4 Z 2 P 2 2 V r 2 P R + V s 2 V r 2 0
Adding both equations, we obtain:
2 V s 2 V r 2 V r 4 2 V r 2 P R + Q X Z 2 P 2 + Q 2 0
The above equation illustrates that equation values decrease with increased power and impedance. Moreover, Equation (104) considers a stability index [19].
S I ( r ) = 2 V s 2 V r 2 V r 4 2 V r 2 P R + Q X Z 2 P 2 + Q 2
Therefore, the classic stability principle determines the stability index for individual lines of radial distribution networks. The minimum stability index value is the most vulnerable to collapse. If the index is equal to zero, then the system collapses. Here, also, neglecting the shunt admittance to determine the stability index will lead to collapse.

2.1.15. Voltage Stability Margin (VSM)

VSM technique [20] is when the admissible rise in load power occurs during the current operating mode or the critical mode. Considering Figure 1, the voltage (Vs) and load apparent power (SL) are derived as follows:
V r = V s Z r Z r 2 + Z 0 2 + Z L Z 0 c o s θ 0 θ r 1 / 2
S r = V r 2 Z r Z r 2 + Z 0 2 + Z L Z 0 c o s θ 0 θ r
Concerning the critical point, the load power attains the maximum value and can be written as:
S c r = V s 2 2 Z 0 1 + c o s θ 0 θ r
Zr is the available load impedance, and the variation between Zr and Z0 is assumed as an invulnerability limit. Accordingly, the VSM in terms of impedance is as follows:
V S M Z = Z r Z 0 Z 0
The allowable power increases until the stability limit is contingent to a higher degree on the angle of the load power factor with a similar VSMs. Therefore, it can be enhanced to signify the margin concerning the load apparent power, VSMs.
V S M s = S c r S L S c r
Hence, the margin of the apparent load power is derived from the VSM. Two assumptions for the stability margin are (i) that the line shunt admittance neglected and (ii) a constant power factor. The system is unstable when the apparent load power is more significant than the critical one. Even if the value of the VSM is less than zero, the system is unstable.

2.1.16. Voltage Reactive Power Index (VQI)

The line index proposed in [21] has a similar methodology to that of Lp. To insert the relation between the voltage (Vs) and reactive power (Qr) comprised from Equation (19) into Equation (18), we obtain:
V r V s Y s r . c o s θ δ V r 2 Y s r . c o s ( θ ) j Q r = V r V s Y s r θ δ V r 2 Y s r θ
V r 2 V r V s s i n θ δ s i n θ + Q r Y s r s i n θ = 0
As the value of δ is negligible and minimized to zero, then the entire term of (sin(θδ)/sin(θ)) is omitted:
V r 2 V r V s + Q r Y s r s i n θ = 0
Replacing the term Yrs sin(θ) with Brs,, the corresponding equation is:
V r 2 V r V s + Q r B s r = 0
Extracting the real roots of the quadratic equation in terms of Vr, we obtain:
V r = V r V r 2 4 Q r B s r 2
When the discriminant of Vr equates to zero, the respective equation contains one distinct real root, or two equal roots. Hence, the respective real roots are:
V r 2 4 Q r B s r 0 4 Q r B s r V r 2 1
The value of Vr ranges from 0 to 1, illustrating the limitation of real roots as the voltage stability limits. Moreover, real roots should be less than 1 and greater than 0, or the stability is affected. Thus, the derived equation evaluates the system stability and identifies the voltage collapse point as the voltage reactive power index.
V Q I = 4 Q r B s r V r 2 1
If the obtained value of VQI is nearer to 1, then the system is within stability limits. Otherwise, a system beyond the stability limits leads to instability. VQI also identifies the instability point of occurrence and voltage collapse.

2.1.17. Line Collapse Proximity Index (LCPI)

An improved LCPI was developed in [22] to overcome the limitations of specific indices [7,8,11] which ignored the line charging reactance while deriving their values. Since the line charging reactance is a critical factor supporting voltage stability, this proximity index was derived considering the definite transmission line model as well as the consequences of active and reactive power flows through lines. The definite transmission line model primarily represents an equivalent circuit of a two-port network considering the ABCD matrix. Thus, the proposed proximity index is derived from the ABCD matrix.
V s I s = A B C D V r I r
Hence, the transmission line parameters for a two-port circuit are A, B, C, and D.
A = 1 + Z × Y 2 , B = Z , C = Y × 1 + Z × Y 4 , D = A Receiving - end   current :   I r = P r j Q r / V r * = P r j Q r / V r δ r
Considering the receiving-end current, the sending-end voltage Vs is derived as:
V s δ s = A α * V r δ r + B β * I r 0 0
A and B are magnitudes, and their corresponding phase angles are α and β , respectively. Combining the current with the voltage equation, we obtain:
V s δ s = A α * V r δ r + B β * P r j Q r / V r δ r
V s V r δ = A α * V r 2 + B β * P r j Q r
Here, δ = δ s δ r . Ignoring the imaginary part and considering the real part of the above quadratic equation, we obtain:
V r 2 A c o s α V r V s c o s δ + P r B c o s β + Q r B s i n β = 0
Obtaining the roots of the quadratic equation allows us to obtain:
V r = V s c o s δ ± V s c o s δ 2 4 A c o s α P r B c o s β + Q r B s i n β 2 A c o s α
It should contain the actual non-zero values of the Voltage equation and the actual non-zero roots obtained from the voltage equation discriminant.
V s c o s δ 2 4 A c o s α P r B c o s β + Q r B s i n β > 0
Therefore, to maintain the system as stable and prevent voltage collapse, the following conditions should be satisfied:
4 A c o s α P r B c o s β + Q r B s i n β V s c o s δ 2 < 1
L C P I = 4 A c o s α P r B c o s β + Q r B s i n β V s c o s δ 2
Therefore, the index value should be less than 1 to maintain the system stability, i.e., LCPI < 1. At no load, the LCPI value is zero, and the LCPI value is adjacent to the unity during the system instability.

2.1.18. New Voltage Stability Index (NVSI)

The NLSI derived in [23] is a similar concept of the power transmission quadratic equation represented in Equation (8). Determining the NLSI is similar to the process for determining Lp. Neglecting δ, Equations (19) and (20) represented as:
V r 2 2 + 2 Q r X V s 2 V r 2 + X 2 P r 2 + Q r 2 = 0
The Voltage Equation Vr is a second-order equation. The limitation to possessing at least one solution is: 2 Q r X V s 2 2 4 X 2 P r 2 + Q r 2 0
2 X P r 2 + Q r 2 4 Q r X V s 2 1 N V S I = 2 X P r 2 + Q r 2 4 Q r X V s 2
Hence, to sustain stability, the value of the NLSI must be less than 1.00 in order toobtain the NLSI by neglecting the line resistance and line shunt admittances.

2.1.19. Integrated Transmission Line Transfer Index (ITLTI)

An integrated transmission line transfer index depends on the radial topology. The ITLTI’s purpose is identifying the power transfer; moreover, it identifies the weak line depending on the transmission line parameters, consequent to the unity, lagging, and leading power factor states. The transmission line in Figure 1 depicts a constant sending-end voltage V s δ transferring the power Ss. Hence, the functional ABCD parameters represent the receiving-end power Sr, Voltage Vr δ with the respective power factor angle θr. Further, the receiving-end complex power is derived as:
S r = V r I r = A V r 2 B β α + V s V r B β δ
Equivalently, the sending-end complex power is represented as:
S s = V s I s = A V s 2 B β α + V s V r B β + δ
Considering the radial transmission line [24], we obtain two-power circles with the same radius and two different centers. Receiving-end complex power circle polar coordinates are represented in Figure 3 with a variable power angle δ.
Considering the receiving-end circle and ΔOCR, we obtain:
C = δ = β α β δ = δ α O = θ r = 180 β α + θ r R = r = 180 δ θ r S r s i n δ = V s V r B s i n θ r = A V r 2 B s i n R
Therefore, the relation between the receiving-end and sending-end voltage can be derived by V r = V s A s i n r s i n θ r . Accordingly, an increased power angle at distinct power factors decreases the receiving-end voltage. The receiving-end power magnitude is:
S r = V s 2 A B s i n θ r + δ s i n δ s i n θ r 2
To deliver the maximum power, Sr is partially derived using δ’ and is equated to zero, which determines the critical value of the power angle.
δ r _ C R = 90 θ r 2 + α
Critical voltage value:
V r ( C R ) = V s 2 A 1 s i n θ r 2
Maximum Power Transfer:
S r ( m a x ) = V s 2 4 A B 1 s i n θ 2 2 2
The maximum limit of the receiving-end power for different power angles is as follows:
S r = S r ( m a x ) s i n θ r + δ s i n δ cos θ r 2 2
Moreover, the integrated transmission line transfer index (ITLTI) is expressed as:
S r _ i n d e x = s i n θ r + δ   s i n δ cos θ r 2 2
Hence, the ITLTI determines the system’s stability for all healthy operations. Determining the Sr_index with varying power angles, the maximum value should be 1; if it exceeds 1, the system is unstable. Therefore, the index value is maintained at less than 1. The loadability concerning the critical power angle increased with an increase in the receiving-end power factor, irrespective of leading or lagging operations. Separating the active and reactive power of Sr, the active power Pr is expressed as follows:
P r = A V r 2 B c o s ( β α ) + V s V r B c o s ( β δ )

2.1.20. Critical Boundary Index (CBI)

The Critical Boundary Index derived in [25] considers the active and reactive power variations and utilizes arithmetical methods to determine the critical boundaries. This index is preferable because of its high accuracy in predicting system performance and analyzing and adding the actual (Equation (8)) and imaginary (Equation (9)) parts of the power flow equation (Equation (7)).
P r R s r + Q r X s r + V r 2 2 + P r X s r Q r R s r 2 = V s 2 V r 2
V r 2 2 + 2 R s r P r + X s r Q r V s 2 2 V r 2 + R s r 2 + X s r 2 P r 2 + Q r 2 = 0
V r 2 = R s r P r + X s r Q r V s 2 2 ± A
Imagining that the intensity of the receiving-end and sending-end voltages is more significant than zero, and rearranging Equation (141), we obtain the following:
A = R s r P r + X s r Q r V s 2 2 2 R s r 2 + X s r 2 P r 2 + Q r 2
Hence, the above equation examines that the power transmission had limits for a specific time duration. Those limits are acknowledged as power and voltage stability limits. These limits are logical only when the magnitude of A is 0. Hence, determining the Qr when A = 0 is written as follows:
Q r = ± X s r 2 V s 4 + R s r 2 V s 4 4 X s r 2 R s r V s 2 P r 4 X s r 3 V s 2 P r 2 R s r 2 + 2 X s r R s r P r X s r V s 2 2 R s r 2
We plot the power stability limit P-Q characteristics with Qr. The P-V characteristic is the nose curve. The power system stability in Figure 4 includes the unstable and stable boundaries, starting from a current stable point K (Po, Qo) and reaching a critical boundary point C (X, Y). Considering the Lagrange Multipliers, Equation (141) simplifies as:
C X ,   Y = R s r X + X s r Y V s 2 2 2 R s r 2 + X s r 2 X 2 + Y 2 = 0
The function f (X, Y) is the length between the current stable point K (Po, Qo) and the adjacent point of voltage collapse C (X, Y). The displacement of f (X, Y) is shown as:
f 2 = X P 0 2 + Y Q 0 2
Here, the ΔPik, ΔQik, stability evaluation index, and displacement in the stability curve are represented in Figure 4. Further, as shown in Figure 1, the operating point between the sending-end and receiving-end transmission line can evaluated by:
P s r = X P 0   &   Q s r = Y Q 0 C B I s r = P s r 2 + Q s r 2
Due to its high-accuracy prediction, this index is preferable. If the index is close to zero, it is the worst-stability system.

2.1.21. Line Voltage Stability Index (LVSI)

A novel LVSI developed to overcome existing stability indices’ drawbacks, LVSI [26], is derived from the quadratic voltage equation. This index is also derived considering the transmission line ABCD parameters; moreover, the inclusion of line charging capacitance and resistance is ignored by existing stability indices. Consequently, this index evaluates the voltage sensitivity accurately, subject to all circumstances, and predicts the voltage collapse point. The receiving-end active power in terms of transmission line parameters is obtained using:
P r = V s V r c o s β δ B A V r 2 cos ( β α ) B
Rearranging the power equation according to a quadratic equation of voltage is obtained as:
V r 2 V s V r c o s ( β δ ) A c o s ( β α ) + P r B A c o s ( β α )
The derivation of the voltage concerning the active power to obtain the sensitivity is shown as:
d V r d P r = B 2 V r A cos β α V s cos β δ
Hence, the sensitivity should be negative to maintain the system’s stability.
B 2 V r A cos β α V s cos β δ < 0
V s cos β δ 2 V r A cos β α < 0
Further simplifying Equation (151), the stability index can be shown as:
L V S I = 2 V r A cos β α V s cos β δ > 1
It is highly stable if the index value is 2 at no load. If the index value is 1 at the maximum load, it is a collapse point. Therefore, the system to maintain the stability index value must be greater than the unity.

2.1.22. New Line Voltage Stability Index (BVSI)

The index BVSI derived in [27] depends on the power transmission approach. Considering Equation (14), the quadratic equation can be rearranged as follows:
R V r 2 V s V r R c o s δ + X s i n δ + P r R 2 + X 2 = 0
Deriving the quadratic equation roots, we obtain:
V r = V s R c o s δ + X s i n δ ± V s R c o s δ + X s i n δ 2 4 R P r R 2 + X 2 2 R
For determining the absolute and non-zero values of Vr, the discriminant is made greater than or equal to 0:
V s R c o s δ + X s i n δ 2 4 R P r R 2 + X 2   0
B V S I = 4 R P r Z 2 V s R c o s δ + X s i n δ 1
Therefore, the index value is less than the unity for the system to be at the stability limit, such that the index value lies between 0 and 1. If the BVSI value exceeds 1, the corresponding line encounters a voltage drop, and the system may collapse. The BVSI line ignores the shunt admittance and the effects of reactive power on stability.

2.2. Bus Voltage Stability Indices

The Bus Voltage Stability Indices evaluate the transmission bus’s voltage stability and do not derive the weak characteristics of the potential voltage issues across the system. Therefore, Bus VSIs are unable to be utilized for determining weak performances.

2.2.1. L-Index

The L-Index [28] is derived from the power flow and standard two-bus system and obtains a detailed analysis. The L index determines peculiar characteristics of power system susceptibility, measuring significant active power, restructuring weak areas or buses, predicting voltage collapse, and identifying instability. Considering the matrix Fij of an ith column and jth row, and elements obtained from the Y matrix, αL, and αG are the load bus and generator bus. The linear transmission system represents a Hybrid matrix (H).
V L V G = H . I L V G = Z L L F L G K G L Y G G . I L V G
Here, VL, IL: Voltage vector, Current vector at load nodes; VG, IG: Voltage vector, Current vector at generator nodes; ZLL, FLG, KGL, YGG: H submatrices. By partial inversion, the H is obtained from the Y-matrix. For any load node j, jαL, and the corresponding voltage equation Vj is:
V j = i α L Z j i . I i + i α G E j i V i
V j 2 + V _ o j V _ j * = S _ j + * Y _ j j +
Substituting the equivalent voltage Voj, admittance Yjj+, and power transferred Si+:
V _ o j = i α G F _ j i V _ i
Y _ j j + = 1 Z _ j j
Hence, Sj+ contains two parts:
S _ j + = S _ j + S _ j c o r r
S _ j c o r r = i α L i j   Z _ j i * Z _ j j * . S _ i V _ i . V _ j
The exceptional characteristic of Equation (159) is that its formation is similar to the admittance equation. Dissimilarities exist between the equivalent voltage Voj and power transferred Sj. The voltage Voj varies slightly, not remaining as a constant value, since the generator voltage remains mostly constant at varying loads. Adopting these limitations can evaluate and regulate local indicators, and Lj can allow us to design for every node j corresponding to the line.
L j = 1 + V _ o j V _ j = S _ j + Y _ j j + * V j 2 = . .
Under stable operations, the Lj ≤ l condition should never disobey any node j. Therefore, the global indicator L represents the stability of the entire subsystem, expressed as:
L = MAX j α L L j L = MAX j α L 1 i α G F _ j i V i V _ j
Here, a series of load nodes is αL, and a series of generator nodes is αG. If the L index value is 0, the system is stable; if the value is 1, the system is approaching instability.

2.2.2. Voltage Instability Proximity Index (VIPI)

The proximity index derived in [29] for estimating instability depends upon coordinating various load flow solutions and instability. The VIPI is a scalar index that analyses the stability margin, considering the angle separating the critical and specified-value vectors. The power flow equation in terms of rectangular coordinates is as follows:
Y s = Y x
Here, Ys is the specified value; x = e 1 ,   f 1 ,   e 2 ,   f 2 , , e n ,   f n T ; ei = Bus voltage real value; fi = Bus voltage imaginary value. Assume two types of voltage vectors (Figure 5): x is the operating solution, and x* is the fictitious solution; comply with the uniform specified-value vector. Employing critical vector (a) and deflection vector (b), x and x* can be represented as x = a + b, x* = ab. Further simplifying them, a = x + x * 2 ,   b = x x * 2 .
As this proximity index is the angle separating critical vector Y(a), and specified value vector Ys, the VIPI is as follows:
V I P I = θ = cos 1 Y s T Y ( a ) Y s . Y ( a )
Here, θ is the angle separating two vectors; critical vector Y(a) is in the space of node specification.

2.2.3. Voltage Collapse Proximity Index (VCPIBUS)

VCPIBUS, proposed in [30], is derived out of the standard power flow equation, and considering a network’s voltage phasor measurement and admittance matrix determines the VCPI at every bus. The complex power at bus k of an N bus network is represented as:
S k * = V k 2 V k c o s δ k j V k s i n δ k m = 1 m k N V m c o s δ m + j V m s i n δ m Y k k
The term V m in the above equation is expressed as:
V m = Y k m j = 1 j k N Y k j V m
Hence, the right-hand term in Equation (168) is complex in the form of a − jb. Let the two equations contain two unknown terms (Vk, δ).
f 1 V k , δ = V k 2 m = 1 m k N V m V k c o s δ
f 1 V k , δ = m = 1 m k N V m V k s i n δ
A partial derivative matrix is attained by simplifying two equations for evaluating the unknowns with the Newton-Raphson Methodology. During the voltage collapse, the matrix determinant is equal to 0, which results as:
V k c o s δ m = 1 m k N V m = 1 2
Further simplifying Equation (172), the voltage collapse prediction index (VCPI) determined at bus k is:
V C P I k t h b u s = 1 m = 1 m k N V m V k
Therefore, the above equation is the voltage collapse condition for bus k. The VCPI value fluctuates between 0 and 1, where 0 represents a stable system, and 1 represents voltage collapse.

2.2.4. S Difference Criterion (SDC)

The S Difference Criterion (SDC) [31] needs a protective relay to perform arithmetic derivation and utilizes the successive evaluation of apparent power (S) at the relay point. During the voltage instability environs, the inflation of sending-end power flow delivers the transmission losses. Consequently, the rise in the sending-end apparent power does not afford a rise in receiving-end power. It specifies that ΔS = 0 at the relay point during the voltage instability. Hence, the transfer of apparent power to the receiving end is:
S _ j ( k ) = V _ j ( k ) I _ j i ( k ) *
Considering the change in a time period Δt, the time interval changes to tk = tk + Δt. A respective rise in apparent power is as follows:
S _ j ( k + 1 ) = S _ j ( k ) + S _ j ( k + 1 ) = V j ( k ) + V _ j ( k + 1 ) . I _ j i ( k ) + I _ j i ( k + 1 ) * = S _ j ( k ) + V _ j ( k + 1 ) I _ j i ( k ) * + V _ j ( k ) I _ j i ( k + 1 ) * + V _ j ( k + 1 ) I _ j i ( k + 1 ) * 0
Here, the term V _ j ( k + 1 ) I _ j i ( k + 1 ) * is very small and can be ignored.
S _ j ( k + 1 ) = V _ j ( k + 1 ) I _ j i ( k ) * + V _ j ( k ) I _ j i ( k + 1 ) * = 0
1 + V _ j ( k + 1 ) I _ j i ( k ) * V _ j ( k ) I _ j i ( k + 1 ) * = 1 + a e j φ = 1 + a c o s φ + j s i n φ = 0
Thus, the SDC is a positive real value; then, considering the magnitude of Equation (177), we obtain:
S D C = 1 + a e j φ = 1 + V _ j ( k + 1 ) I _ j i ( k ) * V _ j ( k ) I _ j i ( k + 1 ) *
Intending to achieve voltage collapse, the value of ΔS = 0, and correspondingly, the SDC equates to 0. The relay operates if the SDC value is below a predetermined threshold.

2.2.5. Impedance Stability Index (ISI)

The ISI bus stability index depends on Tellegen’s rule, an index derived in [32] made simple in evaluating Thevenin’s variables. Let a bus r with a complex load (Sr) and impedance (Zr) be linked to a power system; the remaining system is expressed as Thevenin’s equivalent with the variables E and Zth. Let the current equation be:
Sr/Vr = Ir* = ((EVr)/Zth)*
Vr(EVr)* − SrZth* = 0
Two voltage solutions empower the phasor of Equation (180) with a power Sr; if the two solutions are equal, it allows maximum power transfer.
Vr = (EVr)*
Equation (181) is fundamental to a local voltage stability approach that conveys predetermined consequences through circuit theory.
Z _ r = Z _ t h
The system delivers the maximum power when the load-impedance magnitude equals Thevenin’s impedance magnitude. The adjacent Thevenin network N ^ possesses an identical analysis to the base network N. Hence, the derivation is similar to putting a hat (^) on every term, which signifies the adjoint network parameters. The divergence order of Tellegen’s theorem is equivalent to the following:
I ^ _ r * V _ r V ^ _ r I _ r = 0
Here, V _ r and I _ r * Express the advancement of the complex voltage and current at bus ‘r’ regarding the Vr and Ir* base values. The adjacent network load impedance Z ^ _ r is described as a ratio of phasor voltage V ^ _ r to phasor current I ^ _ r .
Z ^ r = Z _ ^ r = V ^ _ r / I ^ _ r
The load impedance specifies that it is the base case; moreover, it is determined at an instant earlier than that of the network that is regulated to power disturbances.
Z ^ r = V _ r / I _ r
For determining the Thevenin’s impedance, the current I ^ _ r * equation is substituted in Equation (183) from Equation (179).
E ^ _ V ^ _ r / Z _ ^ t h * V _ r V ^ _ r I ^ _ r * = 0
Z ^ _ t h * = E ^ _ V ^ _ r * / V ^ _ r I ^ _ r * V _ r
Similarly, from Equation (181), the condition for the maximum power transfer of adjacent network can be expressed as:
V ^ _ r = E ^ _ V _ ^ r *
Analyzing Equations (187) and (188), during the system collapse point, the equivalent Thevenin’s impedance is expressed as:
Z ^ _ t h * = V _ r / I r *  
Z ^ _ t h = V ^ _ r * / I _ r
However, the above two equations are conjugate with each other; for discarding the conjugate values, Thevenin’s impedance absolute value is as follows:
Z ^ t h = Z _ ^ t h = V _ r / I _ r
The contingency during the regular loading is:
Z ^ r Z ^ t h
During the commencement of voltage instability, the dissimilarities among these impedances converge to 0. Moreover, at that instant of collapse, the impedances are identical, i.e., Z ^ t h = Z ^ r .
V _ r / I _ r = V _ r / I _ r
This index precisely determines Thevenin’s impedance by evaluating the advancement in voltage and current following the base-case exposure to disturbances. Furthermore, the load impedance is evaluated through the ratio of the voltage to the current. Using Equations (192) and (193), a modest regulated impedance stability index (ISI) for identifying the voltage stability margin is:
I S I = Z ^ r Z ^ t h / Z ^ r
Further simplifying by substituting the impedances in terms of voltage and current, the stability index derives as:
I S I = 1 I _ r V _ r / V _ r I _ r

2.2.6. Voltage Stability Index (VSIBUS)

The VSI bus index relies on system elements to identify the distance between the operating and collapse points. This index [33] derives from power flow equations that are similar to the assumptions of the SDC. Let ‘i’ be one local bus of the system. The complex load power of bus i is:
S i = S i θ i = P i + j Q i = V i I i *
The complex bus voltage is Vi = Vi∠δi, and the complex load current is Ii = Ii(δi − θi). Employing Taylor’s theorem, the interrelation among the advancement difference between Vi and Ii as a consequence of advancement evolution in Si is expressed as:
S i = S i I i I i + S i V i V i + h i g h e r   o r d e r   t e r m s
Meanwhile, ignoring the higher-order terms,
S i = V i I i + I i V i
Usually, if the load increments, there is a corresponding increase in load current and decrease in load voltage. This means that, if the value of S i is positive, Current I i is positive and the V i Voltage is negative. Although the bus reaches a critical point or voltage collapse, adding additional load on a bus might not be able to raise the apparent power Si of the load due to rapid decrement in voltage in contrast to the current increment. Going through this case converges ΔSi to zero. Accordingly, the power S i through the voltage-stability limit is expressed as:
0   1 + I i V i V i I i
Considering that the power of Equation (199) is made to α (>1.0), an actual value, it determines almost the same level of linear characteristics. Therefore, excluding the general loss, the voltage stability index at the ‘i’ bus is as follows:
V S I i = 1 + I i V i V i I i α
Hence, at no load, the VSI varies from 1, and at voltage collapse, the value varies to 0.

2.2.7. ZL/ZS Ratio

The ZL/Zs ratio index proposed in [34] depends on a similar technique to the ISI. The correspondence index is the load bus impedance (ZL) ratio to Thevenin’s impedance (Zs). The voltage stability margin was generally evaluated from different scenarios. The evaluation of the stability margin depends on the determination of the voltage stability. For an equivalent circuit, the derivation of the apparent load power concerning load admittance is as follows:
d S d Y = V 2 1 Z S Y 2 1 + Z S Y 2 + 2 Z S Y c o s β
Whereas Y = 1/ZL, dividing both sides of the equation with V2:
d S * d Y * = 1 Z S Y 2 1 + Z S Y 2 + 2 Z S Y c o s β
While dS*/dY* = (Y/S) (dS/dY), thenthe Equation (202) together evaluates (ZL/Zs):
Z L Z S = M + 1 M c o s β + M c o s β 2 M 2 + 1 0.5
Here, M = (dS*/dY*), and the relevance between ZL/ZS and M relies upon the angle β. The estimation of β is a phase angular difference between the source and the load, i.e., β = φSφL. The definite value of φS operates on the existent form of the system. The parameter M is determined by evaluating the derivative of the apparent load power to admittance:
M = d S * d Y * = S 2   S 1 Y 2 +   Y 1 S 2 + S 1 Y 2   Y 1
Here, S and Y refer to the load power and admittance, respectively. Moreover, suffixes 1 and 2 refer to the starting of a change (time, t1) and completion of the change (time, t2), respectively. The variation between the two-time measurements needs to be around 500 ms. If the index value remains above 1, it is a reliable condition.

2.2.8. Equivalent Node Voltage Collapse Index (ENVCI)

The index ENVCI depends on the ESM (Equivalent System Model) and local voltage phasors [35]. The ENVCI has numerous advantages like accurate evaluation and modeling; moreover, this index is convenient for performing real-time and online operations. The outgoing power at node N of a single-line equivalent model (Figure 6) should perform the consequent equation of plain power flow.
P n + j Q n = V n . E k V n Z k n *
Let the phasor voltage equations for two nodes of a single-line equivalent network in terms of rectangular coordinates be E k = e k + j f k and V n = e n + j f n , then the corresponding line impedance is Zkn = Rkn + jXkn. Separating the real terms and imaginary terms of the power flow equation, we obtain:
P n R k n + Q n X k n = e n e k e n + f n f k f n P n X k n Q n R k n = f k e n e k f n
Hence, the power flow equation, Equation (206), is crucial for determining the phasor voltage at the receiving node, describing the system’s performance. In contrast, the phasor voltage of the sending node is known. Whenever the respective ESMs of all equivalent lines possess an arithmetical solution of the receiving node, signifying that they contain every node voltage of the system, at that instant, the system should possess a complete power flow solution; moreover, this should influence the voltage stability. The system’s stability relies on evaluating Equation (206) at every node. The evaluation of Equation (206) can be examined by the Jacobina matrix singularity:
J = e k 2 e n f k 2 f n f k e k
d e t J = 2 e k e n + f k f n e k 2 + f k 2 = 0
Therefore, Equation (208) specifies a novel stability index that is known as the Equivalent Node Voltage Collapse Index:
E N V C I = 2 e k e n + f k f n e k 2 + f k 2
The ENVCI expression in terms of polar coordinates is:
E N V C I = 2 E k V n c o s θ k n E n 2
Determining the ENVCI requires only the phasor voltage of two nodes. Every node contains ESM. The system converges to a collapse point if the ENVCI of a minimum of one node is approaching 0; the respective node is the weak node that induces system instability.

2.2.9. Power Stability Index (PSI)

The power stability index [36] was developed to identify the optimal position of the Distributed Generator (DG), including critical buses that are nearer voltage collapse. This index derives from a standard two-bus system, which is stable within the unity margin. Let the power flow in a standard two-bus system without DG be:
S L = P L + j Q L = V r I r *
V ¯ r = V ¯ r I r Z ¯
where
I r = P L j Q L V r *
Now, with DG, the current Ir concerning the active power and reactive power support are expressed as:
I r = P L P G j Q L Q G V r *
Combining Equations (212) and (214), and writing real and imaginary terms separately, we obtain:
P L P G = V r V s V r * c o s θ δ s + δ r V r 2 Z c o s ( θ )
  Q L Q G = V r V s V r * s i n θ δ s + δ r V r 2 Z s i n ( θ )
Rearranging Equation (215), we obtain:
V r 2 V r V s c o s θ δ c o s ( θ ) + Z ( P L P G ) c o s ( θ ) = 0
Hence, Equation (217) is in a quadratic form in terms of Vr. When stable, the node voltage must contain real roots, i.e., the discriminant is more significant than zero, and the index derives as:
P S I = 4 R i j P L P G V i c o s θ δ 2 1
If the PSI value is less than 1, the system is stable; if the value is nearer to zero, the system is more stable. Hence, the PSI determines the optimal location of the DG as the calculated value of the PSI for each line patterned from maximum to minimum value. The optimal location of the DG is at the line-end with a maximum value of the PSI.

2.2.10. Voltage Deviation Index (VDI)

The Voltage Deviation Index proposed in [37] determines the voltage deviation’s real value compared to 1 per unit. Assuming an N-bus system, the total VDI is the addition of all N buses’ voltage deviation index.
V D I j = 1 V j
V D I T = j = 1 N 1 V j

2.2.11. Simplified Voltage Stability Index (SVSI)

An advanced SVSI proposes [38] to evaluate the system’s stability. This SVSI depends on the Relative Electrical Distance (RED) concept, which considers identifying the convenient generator to a particular load bus; moreover, the relation of electrical parameters for enhancing the performances.
  • Relative electrical distance (RED)
The concerned system, the interrelation among the current (I), Generator bus voltage phasor (V), and the load buses (L) represented in the form of an admittance matrix can be shown as:
I G I L = Y G G Y G L Y L G Y L L V G V L
Rearranging Equation (221), we obtain:
V L I G = Z L L F G L K G L Y G G I L V G
FLG = −|YLL|−1|YLG| represents a composite matrix, which indicates the association between the voltages of the source and load buses. The REDs determined from the matrix FLG are shown as:
R L G = A a b s F L G = A a b s Y L L 1 Y L G
Here, the size of matrix A is (n-g) X g, the total number of buses is represented by n, and the total generator count is denoted by g. Matrix A is a particular matrix in that all elements are 1, i.e., a matrix of ones. The data obtained from matrix RLG are considered rather than path algorithms for determining the electrical distance from the generator bus to the load bus. Subsequently, after identifying the shortest distance of the generator to a particular load with the RLG matrix, the Thevenin impedance voltage drop ΔVi can be estimated with:
V i = b = 1 n j 1 V b V b + 1 V g V i
Vl and Vg are the phasor voltages at the load bus and closest generator. While an increase in demand for power reaches the maximum loading point in one particular bus or various buses in the system, some buses may experience high voltage drops. This process illustrates the potential of voltage collapse and comprises the formation of the SVSI in terms of a correction factor β.
β = 1 m a x V m V l 2
The derived correction factor is proportionate with higher variations between the voltage magnitude of bus m and bus l, which can be determined directly by measuring the PMU during the power system analysis, considering particular operating conditions. Thus, the SVSI is as follows:
S V S I i = V i β × V i
Suppose the obtained value of the proposed index is near the unity. In that case, the system considers the voltage instability on these terms: the Thevenin impedance voltage drop equals the voltage of the load bus.

2.2.12. P-Index

The P-index [39] evaluates the distance to collapse and measures the load shedding. The P-index is derived from a standard two-bus system, shown in Figure 1. The load at the receiving end bus is Pr + jQr, and Vr is the voltage magnitude. The equivalent load admittance is GrjBr,
G r = P r V r 2 , B r = Q r V r 2
Let ΔPr, ΔQr be an incremental increase in load without varying the power factor. The equivalent increment in admittance is ΔGr, ΔBr. The voltage magnitude changed by an amount of ΔVr, is typically negative; therefore, the updated bus voltage is Vr + ΔVr. The corresponding active power variation is as follows:
P r = V r + V r 2 G r + G r V r 2 G r = V r + V r 2 G r + 2 V r + V r G r V r
The substantial consequences of the two terms in Equation (228) are: the primary term is positive and describes the power attained because of the additional load ΔGr, while the secondary term is negative and is the power wasted on the original load Gr because of the voltage drop ΔV. Because of these opposite terms, the overall active power gained is balanced. At the stability limiting point, these two terms diminish, and the overall increment in power is zero.
The P-index bus stability index depends on the ratio of two terms in Equation (228), i.e., the ratio of power lost to power gained. When ΔV is negative for positive ΔGr, a minus sign is added to make the index positive.
P - index = 2 V r + V r G r V r + V r 2 . V r G r
In the limiting case, such as ΔGr, ΔV → 0,
P - index = 2 G r V r . d V r d G r
Although dVr/dGr is often not defined in terms of network terminology, it is simple to do so in terms of the power and voltage sensitivity of the system. Provided that dVr/dGr:
d V r d G r = d V r d P r . d P r d G r
Thus, differentiating the power, Pr = Vr2 Gr is as follows:
d P r = V r 2 d G r + 2 V r 2 G r d V r
Or:
d P r d G r = V r 2 + 2 V r G r d V r d G r
Substituting Equation (233) in Equation (231), we obtain:
d V r d G r = d V r d P r   V r 2 + 2 V r G r d V r d G r
Further simplifying Equation (234), it is reduced as:
d V r d G r = V r 2 d V r d P r 1 2 V r G r d V r d P r
Substituting Equation (235) in the P-index Equation (230), we obtain:
P - index = 2 V r G r d V r d P r 1 2 V r G r d V r d P r
Expressing Equation (236) in terms of active power, we obtain:
P - index = 2 P r V r ×   d V r d P r 1 2 P r V r ×   d V r d P r
Therefore, the P-index is expressed in terms of power and voltage sensitivities. When d V / d P r = , the theoretical value of the stability index is 1.0 at the stability limit. The stability limit varies from 0 (no load) to 1.0 at the collapse point of the system.

3. Summary

The Voltage Stability Indices that are preferred in power system analysis are derived in detail in this paper. Our review of stability indices examined the overall parameters of voltage stability analysis, including the voltage stability margin, stability indices, identifying voltage collapse proximity, and the maximum loadability. This comprehensive review features the significance of definite stability assessment and specifies the effectiveness and vulnerability of various bus and line Voltage Stability Indices.
The line stability indices, like the VSLI, Lmn, LQP, LP, FVSI, NLSI, VQI, and NVSI, and the bus stability indices like the L-Index, Voltage Instability, and VCPI, are derived from the power flow in a standard two-bus system. In contrast, the stability indices like the VCPI, PTSI, VSI_1, VSM, SDC, and PSI are determined by considering the maximum power transfer and maximum power losses. At the same time, stability indices like the SI, LCPI, LVSI, and BVSI derive from the quadratic voltage equation—the overall objectives of the listed stability indices are expressed in Figure 7.
The VSLI stability index depends on the load reactive power demand ratio to the available reactive power. It helps to determine the critical points and identifies the load bus nearer to the collapse point. The Lmn stability index is online monitoring that predicts the voltage collapse; the reliability test on the IEEE 24 bus system results in an accurate voltage collapse prediction. It also identifies the stressed lines and weak areas that are exposed to voltage collapse; moreover, it allows for flexible operation. The LQP is a static indicator that performs faster, evaluates stability, and determines the origin of voltage collapse.
The VCPI is a proximity index that determines the range of the collapse point; this index depends upon the systems’ generation and load characteristics and identifies exposure to voltage collapse. The VCPI is a method for the real-time prediction of voltage collapse. This proximity index illustrates the collapse position and is an alarming tool for preventing system collapse. LP performs voltage stability assessment in a radial distribution system. It predicts instabilities easily when changing the load. The FVSI can determine the collapse point, weakest bus, critical line, and maximum loadability. The VSLBI is used for local monitoring of threshold voltage collapse and emergency management during voltage-sensitive loads. The VSLBI operates within a small stability margin and maintains reactive reserve power. The VSMI presents data about weak locations, evaluates the voltage stability margin, and maintains stable voltage levels during load variations. Moreover, the VSMI depends on the relation between voltage stability and angular difference between the receiving- and sending-end buses. VCPI_1 is an online voltage stability evaluation index that depends on power transmission lines and identifies the weakest transmission lines by distinguishing the minor power outage buses, including the electrical distance.
The Critical Voltage (Vcr) specifies the minimum voltage where the system performs away from voltage collapse; the system is unstable if the Jacobian power flow matrix is singular. In Vcr, the Jacobian power matrix is a static voltage stability index specifying the range of the voltage stability limit. The PTSI is an essential technique that forecasts dynamic voltage collapse, calculates the effect of adding additional equipment like generators or transmission lines, and determines the weak points in the system. VSI_1 is a novel online stability index for predicting steady-state voltage stability. Coordinated with the network simplification method, it presents the stability margin of every load bus, sequentially determining the buses which are vulnerable to voltage collapse.
The NLSI performs by varying the active and reactive power and evaluates the voltage stability more accurately than other indices that vary only the reactive power. Moreover, the NLSI evaluates the voltage collapse point, specifies the order of transmission lines based on loadability, and considers the consistency of wind and solar power output and the response time of the system. The SI identifies the most vulnerable bus that is exposed to voltage collapse in the radial distribution system, is easily implemented for load flow analysis, and does not need extra parameters. The VSM is the variation between the minimal permissible voltage and absolute voltage at a particular bus; it determines the loss of voltage contribution during the voltage collapse. The VSM can efficiently prevent voltage collapse from voltage crisis, and can be measured by PV and QV curves. The VQI is an effective technique that determines the critical bus, detects instability in large-scale systems, and determines the distance of the collapse point.
The LCPI combines the influence of the proportional flow of active and reactive power with ABCD parameters, determines the consequences of different power system arrangements, and specifies the most sensitive lines. The NVSI evaluates the voltage stability effectively by varying the active and reactive power, determining the weakest bus and sensitive line, and utilizes the Reactive Power Sensitivity theory. The ITLTI was first derived for radial transmission networks and subsequently adapted for larger systems. The ITLTI contains extended integrated significance; moreover, it is probably a better resolution for online stability assessment and specifies the weakest line and the specific data about a critical voltage instability range. The CBI operates with the Lagrange Constant Computational Method (LCM); the main motive of this index is to determine the critical boundaries and voltage stability with lesser parameters. Varying the active and reactive power predicts voltage stability margin precisely with better accuracy. The LVSI determines the stability margin by considering the transmission line (ABCD) parameters and is expressed in terms of MVA. The LVSI can determine the voltage collapse point and critical lines, and involves only the information on the phasor bus voltage. The BVSI forecasts the voltage security level, considering various operating conditions and identifies the sensitive lines that are nearer to stability. This index can identify the capability of the bus to place a single DG and size effectively; moreover, it can accurately identify the critical lines and weakest bus.
The L-index in online testing depends on the eigenvalues of the Jacobian matrix of a network; these eigenvalues determine the gap between the voltage collapse point and the current operating point. Moreover, the L-index measures stability linearly, and non-linear behavior may not be measured thoroughly. The VIPI determines very sensitive buses that are exposed to voltage instability and finds the possible precautions for preventing voltage collapse; the VIPI depends on the relation between the instability and various load flow solutions. The VCPI is a real-time prediction that determines the system’s proximity to voltage collapse, which is evaluated by analyzing the present voltage and voltage levels at the point of collapse and utilizing the magnitude of voltage and angle; the admittance matrix forecasts the voltage collapse. The SDC’s utilization in the protection algorithm for voltage collapse depends on the local bus phasor of the voltage and current at every line’s relay point. The relay operation is local to voltage collapse, and the overall rise in apparent power is because of transmission losses to the supply.
The ISI calculates the system’s stability; this index illustrates the ratio of variation in system impedance to load impedance. The ISI is a local stability index that is derived using Tellegen’s theorem and adjoint network. The VSIBUS is a simple index which identifies the distance to the collapse point with the help of local data like the magnitude of voltage and load current, and thoroughly examines the influence of line tripping, reactive power limits, and non-linear variation in the load. The ZL/ZS measures the derivation of apparent power concerning admittance, quickly calculating both the power and admittance results in switching the impedance load tap-changing of a transformer. The SVSI utilizes the data of present operating conditions combined with phasor voltage measurement. The SVSI is easily accessible and less computationally derived.
P-index calculation is simple. It involves the system model, nodal power, and only the single datum of system states. The only computation performed is Jacobian matrix inversion, which can be attained smoothly. Therefore, the P-index is appropriate for online voltage stability analysis. Mainly, the load shedding is evaluated by the P-index. Instead of system-wide shedding, the method makes use of single-load shedding.
The Stability index’s concept, assumptions considered during derivation, mathematical equation, condition for index stability, objective, and some references are represented in Table 1 and Table 2 of the line Voltage Stability Indices and bus Voltage Stability Indices, respectively.
Figure 8 represents the number of articles published on the line and bus Voltage Stability Indices in the last decade. It shows an eventual rise in research in Voltage Stability Indices to enhance voltage stability and supply the desired power under various load conditions. Table 3 represents some distinct characteristics of Voltage Stability Indices that are the optimal location of the DG and sizing of the DG, as well as dependency on VSI.

4. Conclusions

Voltage stability indices are essential numerical parameters for measuring and evaluating stability points in the power system network. These indices’ significant roles help detect instability points early, preventing systems from voltage collapse and power outages, minimizing the actual power losses, and enhancing the system’s overall efficiency. The evaluation of these indices indicates the power system operators’ performances of a transmission system. Further, the stability of a system can improve through proper design, reactive power compensation, and optimal power flow. This article discusses the mathematical modeling of 34 line and bus stability indices. The significance of these indices for different power system problems with suitable objectives is also discussed with the help of suitable works from the literature, around 138. This article may be attractive to a researcher who has started researching power system stability analysis, the placement of reactive power compensators, and DG placement for stability enhancement. The main intention of this review is to summarize and highlight how crucial it is to choose suitable stability indices to preserve systems’ reliability and stability. Developing various stability indices for effective utilization in sustaining the instabilities across various operational scenarios has been undertaken.
The VSIs review various elements and perspectives and derive from the concepts of maximum power transfer through a line (VCPI, PTSI, VSI_1, VSMI), solution of voltage equation (VQIline, VSLI, L-index, SI, VCPIbus LCPI, FVSI, Lmn, Lp, LQP, NVSI, NLSI), P-V Curve (SDC, SI, VIPI), and maximum power transfer theorem (SDC, VSIbus, ISI, ZL/ZS, VCPI_1, VSLBI, SVSI). Most VSIs which correspond to the voltage equation concept should have a solution. A few VSIs neglect the sending-end bus Thevenin impedance (VCPI_1, VSLBI), and few VSIs consider bus Thevenin impedance to be extremely sensitive for a minute variation in parameters of two consecutive measurements (VSIbus, SDC, ZL/ZS ratio, and SDC).
In the procedure to enhance the efficiency and minimize the operational time of DG placement and sizing issues, the primary phase considers the simple line and bus VSIs, and the secondary phase considers the bus and overall VSIs operating with improved efficiency. This review article mainly focuses and presents information on the significant measurement of line and bus stability indices for maintaining the voltage and reactive power in a power system network.
Furthermore, this review can assist in choosing the desired stability index for a given scenario. The main findings of this paper can be highlighted as follows:
  • The power system stability analysis;
  • A comprehensive review of 34 Voltage Stability Indices derived mathematically;
  • Voltage Stability Indices evaluated the sizing and placement of distributed energy sources;
  • Various power system issues and their corresponding application of Voltage Stability Indices were presented;
  • The corresponding data: name, mathematical calculation, concept, assumptions, condition for stability, and objective for each VSI are listed;
  • This review article supports researchers, power system operators, and engineers regarding stability indices.

Funding

This work was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF), which is funded by the Ministry of Education under Grant NRF-2021R111A3053429.

Data Availability Statement

Not Applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Eremia, M.; Shahidehpour, M. Handbook of Electrical Power System Dynamics: Modeling, Stability, and Control; Wiley: Hoboken, NJ, USA, 2013. [Google Scholar]
  2. Kundur, P.; Paserba, J.; Ajjarapu, V.; Andersson, G.; Bose, A.; Canizares, C. Definition and classification of power system stability IEEE/CIGRE joint task force on stability terms and definitions. IEEE Trans. Power Syst. 2004, 19, 1387–1401. [Google Scholar]
  3. Modarresi, J.; Gholipour, E.; Khodabakhshian, A. A comprehensive review of the voltage stability indices. Renew. Sustain. Energy Rev. 2016, 63, 1–12. [Google Scholar] [CrossRef]
  4. Danish, M.S.S.; Senjyu, T.; Danish, S.M.S.; Sabory, N.R.; Narayanan, K.; Mandal, P. A recap of voltage stability indices in the past three decades. Energies 2019, 12, 1544. [Google Scholar] [CrossRef]
  5. Lim, Z.J.; Mustafa, M.W.; Muda, Z.B. Evaluation of the effectiveness of voltage stability indices on different loadings. In Proceedings of the International Power Engineering and Optimization Conference, Melaka, Malaysia, 6–7 June 2012. [Google Scholar]
  6. Rahman, T.K.A.; Jasmon, G.B. A new technique for voltage stability analysis in a power system and improved load flow algorithm for the distribution network. In Proceedings of the International Conference on Energy Management and Power Delivery EMPD, Singapore, 21–23 November 1995. [Google Scholar]
  7. Moghavvemi, M.; Omar, F.M. Technique for contingency monitoring and voltage collapse prediction. IEE Proc. Gener. Transm. Distrib. 1998, 145, 634–640. [Google Scholar] [CrossRef]
  8. Mohamed, A.; Jasmon, G.B.; Yusof, S. A static voltage collapse indicator using line stability factors. J. Ind. Technol. 1998, 7, 73–85. [Google Scholar]
  9. Moghavvemi, M.; Faruque, O. Real-time contingency evaluation and ranking technique. IEE Proc.—Gener. Transm. Distrib. 1998, 145, 517–524. [Google Scholar] [CrossRef]
  10. Moghavvemi, M.; Faruque, M.O. Technique for assessment of voltage stability in the ill-conditioned radial distribution network. IEEE Power Eng. Rev. 2001, 21, 58–60. [Google Scholar] [CrossRef]
  11. Musirin, I.; Khawa, T.; Rahman, A. Novel fast voltage stability index (FVSI) for voltage stability analysis in a power transmission system. In Proceedings of the Student Conference on Research and Development Proceedings, Shah Alam, Malaysia, 17 July 2002. [Google Scholar]
  12. Milosevic, B.; Begovic, M. Voltage-stability protection and control using a wide-area network of phasor measurements. IEEE Trans. Power Syst. 2003, 18, 121–127. [Google Scholar] [CrossRef]
  13. He, T.; Kolluri, S.; Mandal, S.; Galvan, F.; Rastgoufard, P. Identification of weak locations in bulk transmission systems using voltage stability margin index. In Proceedings of the IEEE Power Engineering Society General Meeting, Denver, CO, USA, 6–10 June 2004. [Google Scholar]
  14. Wang, L.; Liu, Y.; Luan, Z. Power transmission paths-based voltage stability assessment. In Proceedings of the IEEE/PES Transmission and Distribution Conference and Exhibition: Asia and Pacific, Dalian, China, 18 August 2005. [Google Scholar]
  15. Sun, H.; Zhou, X.; Li, R. Accuracy analysis of static voltage stability indices based on power flow model. In Proceedings of the 2005 IEEE/PES Transmission and Distribution Conference and Exposition: Asia and Pacific, Dalian, China, 18 August 2005. [Google Scholar]
  16. Nizam, M.; Mohamed, A.; Hussain, A. Performance evaluation of voltage stability indices for dynamic voltage collapse prediction. Appl. Sci. 2006, 6, 1104–1113. [Google Scholar] [CrossRef]
  17. Gong, Y.; Schulz, N.; Guzman, A. Synchro phasor-based real-time voltage stability index. In Proceedings of the PSCE Conference, Atlanta, GA, USA, 29 October–1 November 2006. [Google Scholar]
  18. Yazdanpanah-Goharrizi, A.; Asghari, R. A novel line stability index (NLSI) for voltage stability assessment of power systems. In Proceedings of the 7th WSEAS International Conference on Power Systems, Beijing, China, 15–17 September 2007. [Google Scholar]
  19. Eminoglu, U.; Hocaoglu, M.H. A voltage stability index for radial distribution networks. In Proceedings of the 42nd International Universities Power Engineering Conference, Brighton, UK, 4–6 September 2007. [Google Scholar]
  20. Guiping, D.; Yuanzhang, S.; Jian, X. A new index of voltage stability considering distribution network. In Proceedings of the Asia-Pacific Power and Energy Engineering Conference, Wuhan, China, 27–31 March 2009. [Google Scholar]
  21. Althowibi, F.A.; Mustafa, M.W. Line voltage stability calculations in power systems. In Proceedings of the IEEE International Conference on Power and Energy, Kuala Lumpur, Malaysia, 29 November–1 December 2010. [Google Scholar]
  22. Tiwari, R.; Niazi, K.R.; Gupta, V. Line collapse proximity index for prediction of voltage collapse in power systems. Int. J. Electr. Power Energy Syst. 2012, 41, 105–111. [Google Scholar] [CrossRef]
  23. Kanimozhi, R.; Selvi, K. A novel line stability index for voltage stability analysis and contingency ranking in power system using fuzzy based load flow. J. Electr. Eng. Technol. 2013, 8, 694–703. [Google Scholar] [CrossRef]
  24. Chuang, S.J.; Hong, C.M.; Chen, C.H. Improvement of integrated transmission line transfer index for power system voltage stability. Int. J. Electr. Power Energy Syst. 2016, 78, 830–836. [Google Scholar] [CrossRef]
  25. Furukakoi, M.; Adewuyi, O.B.; Danish, M.S.S.; Howlader, A.M.; Senjyu, T.; Funabashi, T. Critical Boundary Index (CBI) based on active and reactive power deviations. Int. J. Electr. Power Energy Syst. 2018, 100, 50–57. [Google Scholar] [CrossRef]
  26. Ratra, S.; Tiwari, R.; Niazi, K.R. Voltage stability assessment in power systems using line voltage stability index. Comput. Electr. Eng. 2018, 70, 199–211. [Google Scholar] [CrossRef]
  27. Ismail, B.; Wahab, N.I.A.; Othman, M.L.; Radzi, M.A.M. New Line Voltage Stability Index (BVSI) for Voltage Stability Assessment in Power System: The Comparative Studies. IEEE Access 2022, 10, 103906–103931. [Google Scholar] [CrossRef]
  28. Kessel, P.; Glavitsch, H. Estimating the voltage stability of a power system. IEEE Trans. Power Deliv. 1986, 1, 346–354. [Google Scholar] [CrossRef]
  29. Tamura, Y.; Sakamoto, K.; Tayama, Y. Voltage instability proximity index (VIPI) based on multiple load flow solutions in ill-conditioned power systems. In Proceedings of the 27th IEEE Conference on Decision and Control, IEEE, Austin, TX, USA, 7–9 December 1988. [Google Scholar]
  30. Balamourougan, V.; Sidhu, T.S.; Sachdev, M.S. Technique for online prediction of voltage collapse. IEEE Gener. Transm. Distrib. 2004, 151, 453–460. [Google Scholar] [CrossRef]
  31. Verbic, G.; Gubina, F. A new concept of voltage-collapse protection based on local phasors. IEEE Trans. Power Deliv. 2004, 19, 576–581. [Google Scholar] [CrossRef]
  32. Smon, I.; Verbic, G.; Gubina, F. Local voltage-stability index using Tellegen’s theorem. IEEE Trans. Power Syst. 2006, 21, 1267–1275. [Google Scholar] [CrossRef]
  33. Haque, M.H. Use of local information to determine the distance to voltage collapse. In Proceedings of the International Power Engineering Conference, Singapore, 3–6 December 2007. [Google Scholar]
  34. Wiszniewski, A. New criteria of voltage stability margin for the purpose of load shedding. IEEE Trans. Power Deliv. 2007, 22, 1367–1371. [Google Scholar] [CrossRef]
  35. Wang, Y.; Li, W.; Lu, J. A new node voltage stability index based on local voltage phasors. Electr. Power Syst. Res. 2009, 79, 265–271. [Google Scholar] [CrossRef]
  36. Aman, M.M.; Jasmon, G.B.; Mokhlis, H.; Bakar, A.H.A. Optimal placement and sizing of a DG based on a new power stability index and line losses. Int. J. Electr. Power Energy Syst. 2012, 43, 1296–1304. [Google Scholar] [CrossRef]
  37. Yang, C.F.; Lai, G.G.; Lee, C.H.; Su, C.T.; Chang, G.W. Optimal setting of reactive compensation devices with an improved voltage stability index for voltage stability enhancement. Int. J. Electr. Power Energy Syst. 2012, 37, 50–57. [Google Scholar] [CrossRef]
  38. Perez-Londõno, S.; Rodríguez, L.F.; Olivar, G. A Simplified Voltage Stability Index (SVSI). Int. J. Electr. Power Energy Syst. 2014, 63, 806–813. [Google Scholar] [CrossRef]
  39. Kamel, M.; Karrar, A.A.; Eltom, A.H. Development and Application of a New Voltage Stability Index for On-Line Monitoring and Shedding. IEEE Trans. Power Syst. 2018, 33, 1231–1241. [Google Scholar] [CrossRef]
  40. Makasa, K.J.; Venayagamoorthy, G.K. Online voltage stability load index estimation based on PMU measurements. In Proceedings of the 2011 IEEE Power and Energy Society General Meeting, Detroit, MI, USA, 24–28 July 2011. [Google Scholar]
  41. Alias, M.R.N.M.; Verayiah, R.; Yin, L.J. Investigation of Power Flow Tracing Strategy Using an Index in a Transmission System with Solar PV for Congestion Management. In Proceedings of the 11th International Conference on Power and Energy Systems (ICPES), Shanghai, China, 18–20 December 2021. [Google Scholar]
  42. Oukennou, A.; Sandali, A. Analysis and Comparison of Line Voltage Stability Indices. In Proceedings of the 6th International Renewable and Sustainable Energy Conference (IRSEC), Rabat, Morocco, 5–8 December 2018. [Google Scholar]
  43. Salah, M.; Wang, Y.; Chen, T. A novel collapse prediction index for voltage stability analysis and contingency ranking in power systems. Prot. Control. Mod. Power Syst. 2023, 8, 7. [Google Scholar]
  44. Shekhawat, N.; Gupta, A.K.; Sharma, A.K. Voltage Stability Assessment Using Line Stability Indices. In Proceedings of the 3rd International Conference and Workshops on Recent Advances and Innovations in Engineering (ICRAIE), Jaipur, India, 22–25 November 2018. [Google Scholar]
  45. Calma, E.R.B.; Pacis, M.C. Artificial Neural Network-based Voltage Stability Analysis of Power Transmission Networks with Distributed Generation Using Phasor Measurement Unit Synthetic Data. In Proceedings of the IEEE 12th Control and System Graduate Research Colloquium (ICSGRC), Shah Alam, Malaysia, 6 August 2021. [Google Scholar]
  46. Singh, P.; Parida, S.K.; Chauhan, B.; Choudhary, N. Online Voltage Stability Assessment Using Artificial Neural Network Considering Voltage Stability Indices. In Proceedings of the 21st National Power Systems Conference (NPSC), Gandhinagar, India, 17–19 December 2020. [Google Scholar]
  47. Chintakindi, R.; Mitra, A. Detecting Weak Areas in a Real-time Power System on a Static Load Model using Standard Line Voltage Stability Indices. In Proceedings of the IEEE 2nd International Conference on Smart Technologies for Power, Energy, and Control (STPEC), Chhattisgarh, India, 12 November 2021. [Google Scholar]
  48. Verayiah, R.; Abidin, I.Z. A Study on static voltage collapse proximity indicators. In Proceedings of the IEEE 2nd International Power and Energy Conference, Johor Bahru, Malaysia, 1–3 December 2008. [Google Scholar]
  49. Isaiah, A.; Yanxia, S. New performance indices for voltage stability analysis in a power system. Energies 2017, 10, 2042. [Google Scholar]
  50. Salama, M.M.; Saied, E.M.; Abou-Elsaad, M.M.; Ghariany, E.F. Estimating the voltage collapse proximity indicator using an artificial neural network. Energy Convers. Manag. 2001, 42, 69–79. [Google Scholar] [CrossRef]
  51. Chaithra, A.; Modi, S. Power System Vulnerability Assessment using Voltage Collapse Proximity Index. In Proceedings of the 2nd IEEE International Conference on Power Electronics, Intelligent Control and Energy Systems (ICPEICES), Delhi, India, 22–24 October 2018. [Google Scholar]
  52. Chaithra, A.; Modi, S. Evaluation of Power System Stability using Voltage Proximity Index. In Proceedings of the International Conference on Recent Innovations in Electrical, Electronics & Communication Engineering (ICRIEECE), Bhubaneswar, India, 27–28 July 2018. [Google Scholar]
  53. Khunkitti, S.; Premrudeepreechacharn, S. Voltage Stability Improvement Using Voltage Stability Index Optimization. In Proceedings of the 2020 International Conference on Power, Energy, and Innovations (ICPEI), Chiangmai, Thailand, 14–16 October 2020. [Google Scholar]
  54. Xi, Z.; Kong, W.C. Comparison of voltage stability indexes considering dynamic load. In Proceedings of the 2011 IEEE Electrical Power and Energy Conference, Winnipeg, MB, Canada, 3 October 2011. [Google Scholar]
  55. Soyemi, A.; Misra, S.; Oluranti, J.; Ahuja, R. Evaluation of Voltage Stability Indices. In Information Systems and Management Science; Springer: Cham, Switzerland, 2022; Volume 303. [Google Scholar]
  56. Hosseinzadeh, N.; Aziz, A.; Mahmud, A.; Gargoom, A.; Rabbani, M. Voltage Stability of Power Systems with Renewable-Energy Inverter-Based Generators: A Review. Electronics 2021, 10, 115. [Google Scholar] [CrossRef]
  57. Jirjees, M.A.; Al-Nimma, D.A.; Al-Hafidh, M.S.M. Voltage Stability Enhancement Based on Voltage Stability Indices Using FACTS Controllers. In Proceedings of the 2018 International Conference on Engineering Technology and their Applications (IICETA), Al-Najaf, Iraq, 8–9 May 2018. [Google Scholar]
  58. Pacis, M.C.; Antonio, I.J.M.; Bañaga, I.J.T. Under Voltage Load Shedding Algorithm using Fast Voltage Stability Index (FVSI) and Line Stability Index (LSI). In Proceedings of the IEEE 13th International Conference on Humanoid, Nanotechnology, Information Technology, Communication and Control, Environment, and Management, Manila, Philippines, 28–30 November 2021. [Google Scholar]
  59. Ali, M.B.; Kazmi, S.A.A.; Khan, Z.A.; Altamimi, A.; Alghassab, M.A.; Alojaiman, B. Voltage Profile Improvement by Integrating Renewable Resources with Utility Grid. Energies 2022, 15, 8561. [Google Scholar] [CrossRef]
  60. Meena, G.; Verma, K.; Mathur, A. FVSI-Based Meta-heuristic Algorithm for Optimal Load Shedding to Improve Voltage Stability. In Proceedings of the IEEE 10th Power India International Conference (PIICON), New Delhi, India, 25–27 November 2022. [Google Scholar]
  61. Ramírez-Perdomo, S.L.; Lozano, C.A. Evaluation of indices for voltage stability monitoring using PMU measurements. Ing. E Investig. 2014, 34, 44–49. [Google Scholar] [CrossRef]
  62. Arbelaez Ramirez Marcela, L. Monitoring Voltage Stability Margin across Several Transmission Lines Using Synchro Phasor Measurements. Ph.D. Thesis, Iowa State University, Ames, Iowa, 2016. [Google Scholar]
  63. Sajan, K.S.; Kumar, V.; Tyagi, B. Genetic algorithm-based support vector machine for online voltage stability monitoring. Int. J. Electr. Power Energy Syst. 2015, 73, 200–208. [Google Scholar] [CrossRef]
  64. Naganathan, G.S.; Babulal, C.K. Optimization of support vector machine parameters for voltage stability margin assessment in the deregulated power system. Soft Comput. 2019, 23, 10495–10507. [Google Scholar] [CrossRef]
  65. Reddy, M.V.; Murthy, V.S.K. Improvement of voltage stability based on static and dynamic criteria. In Proceedings of the 16th National Power Systems Conference, Hyderabad, India, 15–17 December 2010. [Google Scholar]
  66. Yunhwan, L.; Hwachang, S. Decentralized Load Shedding Method Based on Voltage Stability Margin Index Using Synchrophasor Measurement Technology. Electronics 2018, 7, 277. [Google Scholar]
  67. Dharmapala, K.D.; Rajapakse, A.; Narendra, K.; Zhang, Y. Machine Learning Based Real-Time Monitoring of Long-Term Voltage Stability Using Voltage Stability Indices. IEEE Access 2020, 8, 222544–222555. [Google Scholar] [CrossRef]
  68. Mesgarnejad, H.; Shahrtash, S.M. Power system voltage stability assessment employing phasor measurement units. In Proceedings of the 2nd International Conference on Computer and Automation Engineering (ICCAE), Singapore, 26–28 February 2010. [Google Scholar]
  69. Ankit, S.; Akash, S. Application of supervised learning for voltage stability assessment using GVSM. Electr. Electron. Technol. Open Access J. 2018, 2, 281–287. [Google Scholar]
  70. Hameedullah, Z.; Shah, D.M.S.; Tomonobu, S.; Mikaeel, A.; Abdul, N.; Mohebullah, W.; Mahdi, K.; Paras, M. A Contemporary Novel Classification of Voltage Stability Indices. Appl. Sci. 2020, 10, 1639. [Google Scholar]
  71. Ghassan, S.H.; Husham, I. Assessment of voltage stability based on power transfer stability index using com-putational intelligence models. Int. J. Electr. Comput. Eng. 2021, 11, 2790–2797. [Google Scholar]
  72. Izzri, N.; Mehdi, O.H.; Abdalla, A.N.; Jaber, A.S.; Shalash, N.A.; Noori, Y.L. Fast Prediction of Power Transfer Stability Index Based on Radial Basis Function Neural Network. Int. J. Phys. Sci. 2011, 6, 7978–7984. [Google Scholar]
  73. Nizam, M.; Mohamed, A.; Hussain, A. Dynamic Voltage Collapse Prediction on a Practical Power System Using Power Transfer Stability Index. In Proceedings of the 5th Student Conference on Research and Development, Selangor, Malaysia, 11–12 December 2007; pp. 1–6. [Google Scholar] [CrossRef]
  74. Woldu, T.A.; Ziegler, C.; Wolter, M. A New Method for Prediction of Static and Dynamic Voltage Collapse Using Node Parameters in Large Power Networks. In Proceedings of the IEEE PES Innovative Smart Grid Technologies Europe (ISGT-Europe), The Hague, The Netherlands, 26–28 October 2020. [Google Scholar]
  75. Jamroen, C.; Dechanupaprittha, S. Synchrophasor-based early-warning voltage instability index. In Proceedings of the IEEE Manchester PowerTech, Manchester, UK, 18–22 June 2017. [Google Scholar]
  76. Rohikaa, R.M.; Lakshmi, R.; Sunitha, R.; Ashok, S. Assessment of voltage stability in microgrid. In Proceedings of the International Conference on Electrical, Electronics, and Optimization Techniques (ICEEOT), Chennai, India, 3–5 March 2016. [Google Scholar]
  77. Oukennou, A.; Sandali, A. Assessment and analysis of Voltage Stability Indices in electrical networks using PSAT Software. In Proceedings of the Eighteenth International Middle East Power Systems Conference (MEPCON), Cairo, Egypt, 27–29 December 2016. [Google Scholar]
  78. Mahiraj, R.S.; Vadhera, S. Voltage Stability Assessment Techniques for Modern Power Systems. In Novel Advancements in Electrical Power Planning and Performance; IGI Global: Hershey, PA, USA, 2020; pp. 128–176. [Google Scholar]
  79. Gupta, A.R. Effect of optimal allocation of multiple DG and D-STATCOM in the radial distribution system for minimizing losses and THD. In Proceedings of the 7th International Symposium on Embedded Computing and System Design (ISED), Durgapur, India, 18–20 December 2017. [Google Scholar]
  80. Narayan, K.S.; Kumar, A. Probabilistic analysis of radial distribution network performance with varying wind speed levels. In Proceedings of the 6th IEEE Power India International Conference (PIICON), Delhi, India, 5–7 December 2014. [Google Scholar]
  81. Pattabhi, M.B.; Lakshmikantha, B.R.; Sundar, K.S. A Novel Method for Contingency Ranking Based on Voltage Stability Criteria in Radial Distribution Systems. Technol. Econ. Smart Grids Sustain. Energy 2022, 7, 9. [Google Scholar] [CrossRef]
  82. Sagwal, R.; Kumar, A. Voltage Stability Margin, based congestion management for a hybrid system, considering secure bilateral trading: MINLP approach. In Proceedings of the IEEE 6th International Conference on Power Systems (ICPS), New Delhi, India, 4–6 March 2016. [Google Scholar]
  83. Rao, A.R.N.; Vijaya, P.; Kowsalya, M. Voltage stability indices for stability assessment: A review. Int. J. Ambient Energy 2021, 42, 829–845. [Google Scholar]
  84. Hauqe, M.H. A linear static voltage stability margin for radial distribution systems. In Proceedings of the IEEE Power Engineering Society General Meeting, Montreal, QC, Canada, 18–22 June 2006. [Google Scholar]
  85. Ahmad, N.; Mohamad, F.; Sulaiman, M.; Fazliana, A.; Kadir, A.; Bin, R.O. Classifications of voltage stability margin (VSM) and load power margin (LPM) using probabilistic neural network (PNN). ARPN J. Eng. Appl. Sci. 2017, 12, 5591–5596. [Google Scholar]
  86. Merlyn, M.; Suresh, D.; Venkataramu, P.S.; Nagaraj, M.S. Effect of TCSC on Line Voltage Stability Indices under Single Line Outage Condition. Int. J. Innov. Res. Electr. Electron. Instrum. Control Eng. 2015, 3, 101–105. [Google Scholar]
  87. Singh, B.; Sharma, P.A.K. Voltage Stability Calculations in Power Transmission Lines: Indications and Allocations (IEEE 30 BUS SYSTEM). Int. J. Nov. Res. Electr. Mech. Eng. 2015, 2, 63–66. [Google Scholar]
  88. Mathew, M.; Gosh, S.; Suresh, D.; Ansari, A.A. An Assessment of Voltage Stability Based on Line Voltage Stability Indices and its Enhancement Using TCSC. IOSR J. Electr. Electron. Eng. 2015, 10, 81–88. [Google Scholar]
  89. Zhang, Y.; Huang, L.; Chen, X. An Enhanced Voltage-Reactive Power Index Considering Power Injection Limitations for Power System Stability Assessment. IEEE Trans. Power Syst. 2021, 36, 262–273. [Google Scholar]
  90. Kumar Kavuturu, K.V.; Sai Tejaswi, K.N.V.; Janamala, V. Performance and security enhancement using generalized optimal unified power flow controller under contingency conditions and renewable energy penetrations. J. Electr. Syst. Inf. Technol. 2022, 9, 18. [Google Scholar] [CrossRef]
  91. Alluri, A.; Lanka, R.S.; Rayapudi, S.R. Optimal Location for FACTS Devices: Hybrid Algorithm to Improve Power Flow and System Security. Cybern. Syst. 2022, 54, 938–960. [Google Scholar] [CrossRef]
  92. Ismail, B.; Wahab, N.I.A.; Othman, M.L.; Radzi, M.A.M.; Vijyakumar, K.N.; Naain, M.N.M. A Comprehensive Review of Optimal Location and Sizing of Reactive Power Compensation Using Hybrid-Based Approaches for Power Loss Reduction, Voltage Stability Improvement, Voltage Profile Enhancement, and Loadability Enhancement. IEEE Access 2020, 8, 222733–222765. [Google Scholar] [CrossRef]
  93. Prabhakar, P. Voltage Stability Indices: Formulation and classification perspectives. Int. J. Innov. Technol. Explor. Eng. 2019, 8, 5777–5785. [Google Scholar] [CrossRef]
  94. Salah, M.; Yifei, W.; Tiancong, C. Modern voltage stability index for prediction of voltage collapse and estimation of maximum load-ability for weak buses and critical lines identification. Int. J. Electr. Power Energy Syst. 2023, 145, 1–19. [Google Scholar]
  95. Gitanjali, S.; Kabir, C.; Priyanath, D. Enrichment of voltage stability in power system through novel generalized approximate reasoning based intelligent control with African buffalo optimization approach. Soft Comput. 2022, 27, 7473–7496. [Google Scholar]
  96. Kushwaha, P.K.; Bhattacharjee, C. Research on a selection of appropriate stability index under adverse system conditions for the assessment of voltage stability of an IEEE 14 bus power system. Eur. J. Electr. Eng. 2020, 22, 435–446. [Google Scholar] [CrossRef]
  97. Samy, A.; Venkadesan, A. Evaluation of Voltage Stability Margin in a Power System using an Artificial Neural Network. In Proceedings of the IEEE 10th Power India International Conference (PIICON), New Delhi, India, 25–27 November 2022. [Google Scholar]
  98. Ademola, A.S.; Sun, Y. Power system optimization approach to mitigate voltage instability issues: A review. Cogent Eng. 2023, 10, 2153416. [Google Scholar]
  99. Bode, O.A.; Lotfy, M.E.; Olabisi, B.A.; Rashid, H.O.; Senjyu, H.T.; Narayanan, K. Security-constrained optimal utility-scale solar PV investment planning for weak grids: Short reviews and techno-economic analysis. Appl. Energy 2019, 245, 16–30. [Google Scholar]
  100. Adewuyi, O.B.; Danish, M.S.S.; Howlader, A.M.; Senjyu, T.; Lotfy, M.E. Network Structure-Based Critical Bus Identification for Power System Considering Line Voltage Stability Margin. J. Power Energy Eng. 2018, 6, 97–111. [Google Scholar] [CrossRef]
  101. Akito, N.; Miyara, R.; Alharbi, T.; Prabaharan, N. Dynamic Voltage Stability Assessment in Remote Island Power System with Renewable Energy Resources and Virtual Synchronous Generator. Energies 2021, 14, 5851. [Google Scholar]
  102. Bode, O.A.; Shigenobu, R.; Ooya, K.; Senjyu, T.; Howlader, M.A. Static voltage stability improvement with battery energy storage considering optimal control of active and reactive power injection. Electr. Power Syst. Res. 2019, 172, 303–312. [Google Scholar]
  103. Ratra, S.; Singh, D.; Bansal, R.C.; Naidoo, R.M. Stochastic Estimation and Enhancement of Voltage Stability Margin considering Load and Wind Power Intermittencies. In Proceedings of the IEEE 6th International Conference on Computing, Communication, and Automation (ICCCA), Arad, Romania, 17–19 December 2021. [Google Scholar]
  104. Saifullah, M.K.; Kabir, M.M.; Mia, N.; Tanvir, A. Optimal Placement of FACTS Devices in the Power System Using ANN to Increase the System Loadability. Electronic Systems and Intelligent Computing. In Lecture Notes in Electrical Engineering; Mallick, P.K., Meher, P., Majumder, A., Das, S.K., Eds.; Springer: Singapore, 2015; Volume 686. [Google Scholar]
  105. Kushwaha, P.K.; Bhattacharjee, C. Role of Energy Storage System for improving System Stability considering Customer Behavior Uncertainties. In Proceedings of the International Conference on Contemporary Computing and Applications (IC3A), Lucknow, India, 5–7 February 2020. [Google Scholar]
  106. Adebayo, I.G.; Sun, Y. Voltage Stability Enhancement Capabilities of LTCT and STATCOM in a Power System. In Proceedings of the IEEE PES/IAS PowerAfrica, Cape Town, South Africa, 28–29 June 2018. [Google Scholar]
  107. Shuai, C.; Deyou, Y.; Weichun, G.; Chuang, L.; Guowei, C.; Lei, K. Global sensitivity analysis of voltage stability in the power system with correlated renewable energy. Electr. Power Syst. Res. 2021, 192, 106916. [Google Scholar] [CrossRef]
  108. Adebayo, I.G.; Jimoh, A.A.; Yusuff, A.A.; Sun, Y. An alternative method for the identification of critical nodes leading to voltage instability in a power system. Afr. J. Sci. Technol. Innov. Dev. 2018, 10, 323–333. [Google Scholar] [CrossRef]
  109. Lakshmi, C.; Manjula, S.S. Detection of Proximity to Voltage Collapse Indicator Using Artificial Neural Networks. Int. J. Sci. Eng. Res. 2019, 10, 59–63. [Google Scholar]
  110. Salama, H.S.; Vokony, I. Voltage stability indices–A comparison and a review. Comput. Electr. Eng. 2022, 98, 107743. [Google Scholar] [CrossRef]
  111. Nanba, M.; Huang, Y.; Kai, T.; Iwamoto, S. Studies on VIPI-based control methods for improving voltage stability. Int. J. Electr. Power Energy Syst. 1997, 20, 141–146. [Google Scholar] [CrossRef]
  112. Takehara, A.; Tanaka, Y.; Iwamoto, S. Voltage stability preventive and emergency-preventive control using VIPI sensitivity. Electr. Eng. Jpn. 2002, 143, 22–30. [Google Scholar] [CrossRef]
  113. Gubina, F.; Strmcnik, B. Voltage collapse proximity index determination using voltage phasors approach. IEEE Trans. Power Syst. 1995, 10, 788–794. [Google Scholar] [CrossRef]
  114. Adebayo, I.G.; Jimoh, A.A.; Yusuff, A.A. Prediction of voltage collapse through voltage collapse proximity index and inherent structural characteristics of a power system. In Proceedings of the IEEE PES Asia-Pacific Power and Energy Engineering Conference (APPEEC), Brisbane, QLD, Australia, 15–18 November 2015. [Google Scholar]
  115. Venu, Y.; Kanthala, S. Optimal placement and Sizing of Compensator for Voltage Stability improvement Using Voltage Collapse Proximity Indicator. In Proceedings of the 3rd International Conference on Electronics, Communication and Aerospace Technology (ICECA), Coimbatore, India, 12–14 June 2019. [Google Scholar]
  116. Adebayo, I.G.; Jimoh, A.A.; Yusuff, A.A. Analysis of Voltage Collapse in a Power System Using Voltage Stability Indices. J. Energy Technol. Policy 2018, 8, 40–45. [Google Scholar]
  117. Šmon, I.; Pantoš, M.; Gubina, F. An improved voltage-collapse protection algorithm based on local phasors. Electr. Power Syst. Res. 2008, 78, 434–440. [Google Scholar] [CrossRef]
  118. Maheshkumar, M. Online voltage stability analysis using synchrophasor technology. Int. J. Sci. Res. Dev. 2014, 2, 273–275. [Google Scholar]
  119. Amroune, M.; Bouktir, T.; Musirin, I. Power system voltage instability risk mitigation via emergency demand response-based whale optimization algorithm. Prot Control Mod Power Syst. 2019, 4, 1–14. [Google Scholar] [CrossRef]
  120. Gümüş, T.E.; Emiroglu, S.; Yalcin, M.A. Optimal DG allocation and sizing in distribution systems with Thevenin-based impedance stability index. Int. J. Electr. Power Energy Syst. 2023, 144, 108555. [Google Scholar] [CrossRef]
  121. Cuffe, P.; Milano, F. Validating Two Novel Equivalent Impedance Estimator. IEEE Trans. Power Syst. 2018, 33, 1151–1152. [Google Scholar] [CrossRef]
  122. Sadeghi, S.E.; Akbari Foroud, A. A new approach for static voltage stability assessment in distribution networks. Int. Trans. Electr. Energy Syst. 2019, 30, 1–21. [Google Scholar] [CrossRef]
  123. Stefan, P.; Herwig, R. Voltage stability monitoring methods for distribution grids using the Thevenin impedance. CIRED—Open Access Proc. J. 2017, 2017, 1535–1539. [Google Scholar]
  124. Zhou, Q.; Guan, Y.; Liu, Y.; Zhao, Z. A novel online local voltage stability monitoring approach based on clustering analysis. Int. J. Electr. Power Energy Syst. 2021, 132, 106807. [Google Scholar]
  125. Pan, J.; Dong, A.; Fan, J.; Li, Y. Online Static Voltage Stability Monitoring for Power Systems Using PMU Data. Math. Probl. Eng. 2020, 2020, 1–8. [Google Scholar] [CrossRef]
  126. Alzaareer, K.; Saad, M. Real-Time Voltage Stability Monitoring in Smart Distribution Grids. In Proceedings of the International Conference on Renewable Energy and Power Engineering (REPE), Toronto, ON, Canada, 24–26 November 2018. [Google Scholar]
  127. Lee, Y.; Han, S. Real-Time Voltage Stability Assessment Method for the Korean Power System Based on Estimation of Thévenin Equivalent Impedance. Appl. Sci. 2019, 9, 1671. [Google Scholar] [CrossRef]
  128. Wang, Y.; Wang, C.; Lin, F.; Li, W.; Wang, L.Y.; Zhao, J. A new transfer impedance-based system equivalent model for voltage stability analysis. Int. J. Electr. Power Energy Syst. 2014, 62, 38–44. [Google Scholar] [CrossRef]
  129. Akhil, P.; Atul, M.; Prabhakar, K.; Kothari, D.P. Assessment on various node voltage stability indices—A review. In Proceedings of the International Conference on Smart Technologies for Smart Nation (SmartTechCon), Bengaluru, India, 17–19 August 2017. [Google Scholar]
  130. Nagendra, P.; Datta, T.; Halder, S.; Paul, S. Power System Voltage Stability Assessment Using Network Equivalents-A Review. J. Appl. Sci. 2010, 10, 2147–2153. [Google Scholar] [CrossRef]
  131. Xinglong, W.; Zheng, X.; Zheren, Z. Power Stability Analysis and Evaluation Criteria of Dual-Infeed HVDC with LCC-HVDC and VSC-HVDC. Appl. Sci. 2021, 11, 5847. [Google Scholar]
  132. Xu, Z.; Xu, Z.; Xiao, L. Analysis and assessment standards of power stability of multi-send HVDC systems. J. Eng. 2018, 2019, 748–753. [Google Scholar] [CrossRef]
  133. Alam, M. Determination of Power System Contingency Ranking Using Novel Indices. Recent Advances in Power Systems. In Lecture Notes in Electrical Engineering; Gupta, O.H., Sood, V.K., Malik, O.P., Eds.; Springer: Singapore, 2022; Volume 812. [Google Scholar]
  134. Mehebub, A.; Shubhrajyoti, K.; Sankar, T.S.; Sumit, B. Determination of Critical Contingency Based on L-Index and Impact Assessment on Power System. In Smart Energy and Advancement in Power Technologies; Springer: Berlin, Germany, 2022; pp. 583–596. [Google Scholar]
  135. Khunkitti, S.; Boonluk, P.; Siritaratiwat, A. Optimal Location and Sizing of BESS for Performance Improvement of Distribution Systems with High DG Penetration. Int. Trans. Electr. Energy Syst. 2022, 2022, 6361243. [Google Scholar] [CrossRef]
  136. Rodríguez Lopez, D.K.; Perez-Londoño, S.; Rodriguez-Garcia, L. Optimal under voltage load shedding based on voltage stability index. Ing. E Investig. 2016, 36, 43–50. [Google Scholar]
  137. Shabir, M.; Nawaz, S.; Vijayvargiya, A. Voltage Stability Enhancement Using SVC in PSCAD Software. In Innovations in Electrical and Electronic Engineering. Lecture Notes in Electrical Engineering; Favorskaya, M.N., Mekhilef, S., Pandey, R.K., Singh, N., Eds.; Springer: Singapore, 2021; Volume 661. [Google Scholar]
  138. Akbarzadeh Aghdam, P.; Khoshkhoo, H. Voltage stability assessment algorithm to predict power system loadability margin. IET Gener. Transm. Distrib. 2020, 14, 1816–1828. [Google Scholar] [CrossRef]
Figure 1. Standard two-bus model.
Figure 1. Standard two-bus model.
Energies 16 06718 g001
Figure 2. Simple two-bus Thevenin Equivalent Circuit.
Figure 2. Simple two-bus Thevenin Equivalent Circuit.
Energies 16 06718 g002
Figure 3. The receiving-end circle.
Figure 3. The receiving-end circle.
Energies 16 06718 g003
Figure 4. P-Q Characteristics.
Figure 4. P-Q Characteristics.
Energies 16 06718 g004
Figure 5. Concept of VIPI in the Node-Specification Space.
Figure 5. Concept of VIPI in the Node-Specification Space.
Energies 16 06718 g005
Figure 6. Equivalent System Model.
Figure 6. Equivalent System Model.
Energies 16 06718 g006
Figure 7. The objective of Voltage Stability Indices.
Figure 7. The objective of Voltage Stability Indices.
Energies 16 06718 g007
Figure 8. Number of voltage stability indices articles published.
Figure 8. Number of voltage stability indices articles published.
Energies 16 06718 g008
Table 1. Line Voltage Stability Indices.
Table 1. Line Voltage Stability Indices.
VSIConceptAssumptionEquationCondition for StabilityObjectiveReferences
VSLI
(1995)
Power flow in a transmission line of a two-bus system.Shunt admittance is ignored.VSLI = 4 V s V r cos δ V r 2 cos δ 2 V s 2 Stable: VLSI < 1
Unstable: VLSI > 1
Identify the critical buses near voltage collapse.[40,41,42]
Lmn
(1998)
Power flow in a transmission line of a two-bus system.The effect of active power is ignored, and the shunt admittance is approximately zero. L m n = 4 x Q r V s s i n ( θ δ 2 Stable: Lmn < 1
Unstable: Lmn > 1
Online monitoring that predicts the voltage collapse and identifies the stressed condition.[43,44,45]
LQP
(1998)
Power flow in a transmission line of a two-bus system.Y ≈ 0 &
R/Z << 1
LQP = 4 X V s 2 Q r + X P s 2 V s 2 Stable: LQP < 1
Unstable:
LQP > 1
Performs faster, evaluates the static voltage collapse.[46,47,48]
VCPI
(1998)
Maximum Power loss & Maximum Power transfer in a Transmission Line.Constant Power Factor Constant Shunt admittance is ignored. V C P I   ( 1 ) = P r P r ( m a x )
P r ( m a x ) = V s 2 Z s   c o s 4 cos 2 θ 2
V C P I   ( 2 ) = Q r Q r ( m a x )
Q r ( m a x ) = V s 2 Z s   s i n 4 cos 2 θ 2
V C P I   ( 3 ) = P l P l ( m a x )
P l ( m a x ) = V s 2 Z s   c o s θ 4 cos 2 θ 2
V C P I   ( 4 ) = Q r Q r ( m a x )
Q l ( m a x ) = V s 2 Z s   s i n θ 4 cos 2 θ 2
Stable: VCPI < 1
Unstable: VCPI > 1
Determines the range of the collapse point, which depends upon the systems’ generation load characteristics.[49,50,51,52]
Lp
(2001)
Power flow in a transmission line of a two-bus system.Shunt admittance is ignored, and the effect of reactive power. L p = 4 R P r V s cos ( θ δ ) 2 Stable: Lp < 1
UnstableLp > 1
Voltage stability assessment in a radial distribution. Simple index with higher accuracy.[53,54,55,56]
FVSI
(2002)
Power flow in a transmission line of a two-bus system.Shunt admittance is ignored and sinδ ≈ 0, cosδ ≈ 1, Rsinδ ≈ 0, Xcosδ ≈ X. F V S I = 4 Z 2 Q r V s 2 X Stable: FVSI < 1
Unstable:
FVSI > 1
Determines collapse point, weakest bus, critical line, and maximum loadability.[57,58,59,60]
VSLBI
(2003)
During the maximum power conditions, the voltage drop in impedance equals load bus voltage.Thevenin equivalent impedance connected to the sending end bus is ignored. V S L B I = V r V Stable: VSLBI > 1
Unstable:
VLSBI < 1
The Access to voltage collapse is local monitoring, watchful, and emergency control during voltage-sensitive load.[61,62]
VSMI
(2004)
VSMI considers the relation between the angular difference of voltage and maximum power transfer.The shunt admittance neglected V S M I = δ r m a x δ δ r m a x Stable: VSMI > 0
Unstable:
VLMI < 0
Evaluate stability margin and determine weak locations.[63,64,65,66]
VCPI_1
(2005)
The voltage drop across the Thevenin impedance is equal to the load; at the collapse point.Thevenin equivalent impedance connected to the sending end bus is ignored. V C P I _ 1 = V r c o s δ 0.5 V s VCPI_1 ≥ 0: Stable
VCPI_1 < 0: Unstable
Online evaluation index. Identify the weakest lines by distinguishing the minor power outage buses with distance.[67,68]
Vcr
(2006)
The load flow equations and Eigenvalue theorem.Constant Power factor and load. V c r = E 2 1 + c o s θ
V c r = E 2 c o s θ
The system is unstable if the Jacobian power matrix is singularSpecifies the minimum voltage where the system performs away from collapse. [69,70]
PTSI
(2006)
Maximum Power loss and maximum power transferable through a line are limited.Shunt admittance is ignored. P T S I = 2 S L Z 1 + c o s θ φ V s 2 Stable: PTSI < 1
Unstable:
PTSI > 1
Predicts the dynamic voltage collapse and calculates the effect of adding additional equipment. [71,72,73,74]
VSI_1
(2006)
Maximum Power loss and maximum power transferable through a line are limited.The resistance of line and shunt admittance is ignored. V S I = m i n P m a r g i n P m a x ,   Q m a r g i n Q m a x ,   S m a r g i n S m a x Stable: VSI_1 > 0
Unstable:
VSI_1 < 0
Predicts steady-state voltage stability, determines the stability margin of every load bus.[75,76]
NLSI
(2007)
Power flow in a transmission line of a two-bus system,
“Critical Clearing Time” (CCT).
The minimal angular difference between receiving and sending voltage and the shunt admittance is ignored. N L S I = R s r P r + X s r Q r 0.25 V s 2 Stable: NLSI < 1
Unstable: NLSI > 1
By varying active and reactive power, evaluating the collapse point, Rank transmission lines.[42,77,78]
SI
(2007)
Voltage Quadratic Equation.The shunt admittance is ignored. S I ( r ) = 2 V s 2 V r 2 V r 4 2 V r 2 P R + Q X Z 2 P 2 + Q 2 Stable: SI ≠ 0
Unstable: SI = 0
Predict the most vulnerable bus exposed to collapse in the radial distribution system.[79,80,81]
VSM
(2009)
Maximum Power loss and maximum power transfer in a transmission line.Constant power factor and shunt admittance is ignored. V S M Z = Z r Z 0 Z 0
S c r = V s 2 2 Z 0 1 + c o s θ 0 θ r
Stable: VSM > 0
Unstable:
VSM < 0
Determining the loss of voltage, VSM can efficiently prevent voltage collapse.[82,83,84,85]
VQI
(2010)
Power flow in a transmission line of a two-bus system.Zero angular difference and shut admittance are ignored. V Q I = 4 Q r B s r V r 2 1 Stable: VQI < 1
Unstable: VQI > 1
Determines the critical bus, detects the instability in large-scale systems, and the distance of the collapse point.[86,87,88,89]
LCPI
(2012)
Voltage Quadratic Equation.The transmission lines model is like π model. L C P I = 4 A c o s α P r B c o s β + Q r B s i n β V s c o s δ 2   Stable: LCPI < 1
Unstable: LCPI > 1
Combines the influence of the relative flow of reactive and active power flow with ABCD.[90,91,92,93]
NVSI
(2013)
Power flow in a transmission line. Reactive Power Sensitivity (RPS).The Shunt admittance and line resistance are ignored. N V S I = 2 X P r 2 + Q r 2 4 Q r X V s 2 Stable: NVSI < 1
Unstable: NVSI > 1
Evaluate the voltage stability effectively by varying the P & Q, determining the weak bus and sensitive line. [94,95,96,97]
ITLTI
(2016)
ABCD parameters, the power factor of the receiving end, and the power angle between the receiving and sending ends.Two power circles with two distinct centers but identical radius. S r _ i n d e x = s i n θ r + δ s i n δ cos θ r 2 2 Stable: ITLTI < 1
Unstable: ITLTI ≥ 1
For radial transmission networks and subsequently adapted for larger systems, the Weakest line. [4,98]
CBI
(2018)
Active and Reactive power changes.Negligible system impedance, linearized power flow model. C B I s r = P s r 2 + Q s r 2 Stable: CBI > 1
Unstable: CBI = 0
Operates with Lagrange Constant Computational Method (LCM); Determines the critical boundaries and voltage stability with lesser parameters.[99,100,101,102]
LVSI
(2018)
Voltage Quadratic Equation, ABCD Parameters.line’s resistance and charging capacitance are ignored. L V S I = 2 V r A cos β α V s cos β δ Stable: LVSI > 1
Unstable: LVSI < 1
Evaluate the stability margin considering the ABCD parameters expressed in MVA.[103,104,105]
BVSI
(2022)
Voltage Quadratic Equation.Line shunt admittance and the reactive power’s effects are ignored. B V S I = 4 R P r Z 2 V s R c o s δ + X s i n δ Stable: BVSI < 1
Unstable: BVSI > 1
Optimal location and sizing of distributed generations.[27]
Table 2. Bus Voltage Stability Indices.
Table 2. Bus Voltage Stability Indices.
VSIConceptAssumptionEquationCondition for StabilityObjectiveReferences
L–Index
(1986)
Power flow equation solution, Eigen Values.All generator voltages remain constant. L = MAX j α L 1 i α G F _ j i V i V _ j Stable: L-Index < 1
Unstable: L-Index > 1
Identify the critical points of the system.[106,107,108]
VIPI
(1989)
Power flow equation solution.steady-state condition, system impedance is negligible. V I P I = θ = cos 1 Y s T Y ( a ) Y s . Y ( a ) The operational solutions estimate critical points and a proximate fictitious solution.The potential voltage instability problems and the efficient control approach for avoiding instability.[109,110,111,112]
VCPI
(2004)
Power Flow Equation.It utilizes offline and online measurement. V C P I k t h b u s = 1 m = 1 m k N V m V k Stable: VCPI < 1
Unstable: VCPI > 1
The system’s proximity to voltage collapse.[113,114,115,116]
SDC
(2004)
Maximum Power Transfer Theorem.Minimal values, such as a change in voltage at receiving, are ignored. S D C = 1 + V _ j ( k + 1 ) I _ j i ( k ) * V _ j ( k ) I _ j i ( k + 1 ) * SDC > 0: Stable
SDC < 0: Unstable
Used to protect voltage collapse depends on local bus phasor current and voltage at every line’s relay point.[117,118,119]
ISI
(2006)
Maximum transferred when Thevenin impedance’s magnitude equals the load impedance’s amplitude.Constant system topology. I S I = Z ^ r Z ^ t h / Z ^ r
I S I = 1 I _ r V _ r / V _ r I _ r
ISI > 0: Stable
ISI < 0: Unstable
Calculates the stability of the system.[120,121,122,123]
VSIBUS
(2007)
The rise in sending apparent power no longer yields a rise in receiving line power. Minimal values, such as a change in voltage at receiving, are ignored. V S I i = 1 + I i V i V i I i α The VSIBUS value lies between unity and zeroIdentifying the distance to the collapse point with the help of local voltage.[124,125]
ZL/ZS Ratio
(2007)
Maximum transferred when the magnitude of Thevenin impedance is equal to the amplitude of the load impedance.730 ≤ ϕs ≤ 870 Z L Z S = M + 1 M c o s β + M c o s β 2 M 2 + 1 0.5
M = d S * d Y * = S 2 S 1 Y 2 + Y 1 S 2 + S 1 Y 2 Y 1
Stable :   Z L Z S Ratio > 1
Unstable :   Z L Z S Ratio < 1
ENCVI is accurate in design and measuring and simple in real-time implementation.[110,126,127]
ENVCI
(2009)
Equivalent system model (ESM) & Equivalent local network model (ELNM).Consider the effects of the local network and the system outside the local network. E N V C I = 2 e k e n + f k f n e k 2 + f k 2 Stable: ENVCI > 0
Unstable: ENVCI < 0
Optimal placement depends on ESM and considers only local voltage phasors.[4,128,129,130]
PSI
(2012)
Maximum power transfer theorem.DG depends on the most critical bus. P S I = 4 R i j P L P G V i c o s θ δ 2 Stable: PSI < 1
Unstable: PSI > 1
DG and sizing for distribution networks.[110,131,132]
VDI
(2012)
Negligible phase angle deviations.The real value of the deviation of the bus voltage. V D I T = j = 1 N 1 V j 0: Perfect voltage regulation
1: worst case voltage regulation.
Optimal placement of DG and sizing for distribution network.[133,134,135]
SVSI
(2014)
During the maximum Power, the load bus voltage equals the impedance voltage drop across the line.The voltage at the nearest generator to the load bus equals the Thevenin load voltage. S V S I i = V i β × V i Stable: SVSI < 1
Unstable: SVSI > 1
Uses the data of present operating conditions combined with phasor voltage measurement.[136,137,138]
P-Index
(2017)
The ratio of
power loss to
power gained.
For incremental increase, the power factor is unchanged. P - index = 2 P r V r d V r d P r 1 2 P r V r d V r d P r At stability limit when dVr/dPr = ∞, the value would be 1Dynamic voltage stability assessment and load shedding purposes.[39]
Table 3. Characteristics of Stability Indices.
Table 3. Characteristics of Stability Indices.
CharacteristicStability Indices
Optimal placement of DG & DG sizingLine VSIs, Bus VSIs except for SDC, VSIbus, ISI, and ZL/ZS ratio
Impedance dependent VSIVSLI, L VCPI_1, VSLBI, ISI, SDC, VSIbus, ZL/ZS ratio
Independent VSIVSMI, SI, LCPI, VCPIBUS, NLSI, VCPI, NVSI, SVSI, FVSI, Lmn, LQP, Lp, VIPI
Reduce Power LossesLine VSIs, Bus VSIs
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Valuva, C.; Chinnamuthu, S.; Khurshaid, T.; Kim, K.-C. A Comprehensive Review on the Modelling and Significance of Stability Indices in Power System Instability Problems. Energies 2023, 16, 6718. https://doi.org/10.3390/en16186718

AMA Style

Valuva C, Chinnamuthu S, Khurshaid T, Kim K-C. A Comprehensive Review on the Modelling and Significance of Stability Indices in Power System Instability Problems. Energies. 2023; 16(18):6718. https://doi.org/10.3390/en16186718

Chicago/Turabian Style

Valuva, Chandu, Subramani Chinnamuthu, Tahir Khurshaid, and Ki-Chai Kim. 2023. "A Comprehensive Review on the Modelling and Significance of Stability Indices in Power System Instability Problems" Energies 16, no. 18: 6718. https://doi.org/10.3390/en16186718

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop