Research Article  Open Access
Ombuki Mogaka, Roy Orenge, Julius Ndirangu, "Static Voltage Stability Assessment of the Kenyan Power Network", Journal of Electrical and Computer Engineering, vol. 2021, Article ID 5079607, 16 pages, 2021. https://doi.org/10.1155/2021/5079607
Static Voltage Stability Assessment of the Kenyan Power Network
Abstract
In recent years, the Kenyan Power Network has witnessed large growths in load demand. Although the increased load demand has somewhat been matched with an increase in transmission and generation capacity, the rate of expansion has not been matched with the rate of increase in load demand due to economic, environmental, and geographical constraints. This has led to the system being prone to instability since it is being operated under stressed conditions. In the recent past, several studies have been carried out on voltage stability analysis and improvement using various conventional methods. However, conventional methods have various limitations in their utilization for voltage stability analysis. One solution to overcome these limitations is to employ a combination of one or more methods so as to get more information and greater degree of accuracy in voltage stability studies. In this paper, a methodology is proposed involving the combination of QV modal analysis, sensitivity analysis (VQ) and powervoltage curves in assessing the static voltage stability analysis taking a case study of the Kenyan Power Network. VQ sensitivity analysis and QV modal analysis have been used to identify the load regions most susceptible to voltage instability and the corresponding weak buses in the network for various VQ responses. Reactive power loss sensitivities for branches in the network have been used to determine the critical (weak) lines in the network. Loading margins (LM) and voltage stability margins (VSM) have then been used to determine the proximity to voltage collapse of the voltage weak buses identified by QV modal analysis. The effect of tripping one the critical lines on the voltage weak buses is also investigated. The current high voltage power network under the average peak loading conditions during the year 2019 is considered for the study. The paper also reviews existing voltage stability analysis methods and their limitations.
1. Introduction
Due to economic and environmental constraints, bulk power systems are increasingly being operated close to their limits which have made them more prone to instability and power quality challenges. Voltage instability is one of the problems that has gained notoriety in recent years. It has been recognized as one of the major sources of power system insecurity and has been responsible for several partial or full power interruptions over the past decades. In fact, more than 50 cases of voltage collapse have been reported around the world in the past 50 years. Notable incidents of interruptions that arose from voltage instability include the ones of France (1978), Belgium (Aug 1982), Sweden (Dec 1983), Tokyo, Japan (July 1987), USA Tennessee (Aug 1987), Canada (March 1989), and the USA’s WSCC power grid blackouts (July 2 and Aug. 10 1996). These incidents led to power interruptions to millions of consumers in some of the countries and losses in the range of millions. Some of the incidents reported over the years are presented in Table 1 together with the time frames and total load interruptions reported [2–10].

Several definitions of voltage stability and related terms are provided in the literature [11–15]. In general, voltage stability is associated with the ability of a power system to maintain acceptable voltages at all the buses in a power system under normal operations conditions and after being subjected to a disturbance. Voltage instability on the other hand refers to the inability of a system to maintain acceptable voltages in some or all buses in the system under normal operating conditions and after being subject to disturbances, i.e., lowprobability single or multiple contingencies, which cause a gradual and uncontrollable decline in voltage. It can also result from a mismatch between generation and demand of reactive power, the case in which a power system fails to meet the reactive power requirements. If corrective measures are not initiated to check voltage instability, a significant decrease in system voltage occurs, due to heavy reactive power losses. Consequently, voltage collapse occurs. Voltage collapse is a process by which a sequence of events accompanying voltage instability leads to an unacceptably low voltage profile in a significant part of a system [12, 14]. It is characterised by heavy reactive power flows, inadequate reactive power support, and heavily loaded systems and results in partial or total system blackouts.
In Kenya, the government has been on a drive to connect more consumers to the national grid. The electricity access rate in Kenya has been on a steady rise over the last six years from 27% in 2013 to 77% in 2019. The number of connected customers had risen from 1.5 million in 2010 to 7 million as of 2019. This increase in demand has not been matched by similar completion of transmission line projects. Consequently, the system is being operated closer to its stability limit and is thus more prone to instability in case of disturbances. This has put the system controllers under more pressure as they have to operate the system within lower security margins and defensively operate the system during conditions of peak loads. Therefore, the voltage instability problem leads to degradation in a power systems’ performance and jeopardises the utility company’s objective of delivering uninterruptible and reliable power supply to its customers. In order to avoid such occurrences, it is important that voltage stability analysis of the power system be undertaken.
Voltage stability analysis can be divided into dynamic and static voltage stability analyses. Both methods involve the examination of how far a system is to voltage instability and the mechanism of voltage instability, i.e., how and why voltage instability arises, the key factors contributing to instability, identification of voltage weak areas the network under study, and identifying the most effective measures of counteracting instability and by extension improving voltage stability [15, 16].
Dynamic voltage stability analysis is based differential equations and deals with how bus voltages vary with changes in system operating parameters. It is useful in studying voltage collapse scenarios and understanding the chronology of events that lead to voltage collapse. Dynamic voltage stability methods include small signal stability analysis, time domain simulations, bifurcation analysis, and energy function method [16–22]. The features of the popular dynamic voltage stability methods are summarized in Table 2.

It is however noted that although voltage stability is a dynamic problem, much of the problem can be solved through static analysis since it yields most of the required information with regards to the voltage stability status of a system. Static analysis also takes less computational time and effort when compared to dynamic analysis and hence easier to carry out. Static voltage stability analysis involves computation of power flower equations. It is basically taking “snapshots” of probable operating conditions which can be used to validly point out the mechanism of voltage collapse for different operating conditions. Some of the existing methods used in static voltage stability analysis are PV and QV curves, sensitivity analysis, singular value decomposition, continuation power flow method (CPF), and modal analysis [12, 15, 16, 23]. Features of these methods are summarized in Table 3.

PV and QV curves are among the most popular methods when it comes to static voltage stability analysis. PV curves are plots of load bus voltage magnitudes against load active power for a given power factor, whereas QV curves are plots of voltage against the reactive power injected into a bus. Static voltage stability analysis with these curves involves finding the extremes of these curves (nose points) by load flow calculations, where the distance between the current operating points and the extremes is the stability criterion. These curves only give the proximity to voltage collapse points but do not give information about the causes of voltage instability. Continuation power flow (CPF) method was proposed by Ajjarapu and Christy in [24, 25]. This method overcomes the challenges that arise with power flow convergence when using PV curves, as computation approaches the nose points or voltage collapse point. To resolve these power flow convergence issues, CPF utilizes the predictorcorrector process to find the next operating point. The major shortcomings of PV, QV curves, and CPF method are the required computational effort and time since they are only related to a specific bus.
The singular value decomposition method is used to determine voltage stability margin by observing how close the Jacobian matrix of the power flow equations is, to becoming singular. In the analysis, the load is increased and the power flow Jacobian matrix, J, computed [26]. The loadability limit is reached when the smallest singular value of determinant of J becomes 0. This method cannot find a specific cause for voltage instability but provides a relative proximity to the voltage stability limit. Even though singular value decomposition is based on linearized power flow equations, the Jacobian matrix of power flow equations has a high degree of nonlinearity after stable operating point to the bifurcation point; hence, it is very difficult to evaluate the relative degree of voltage stability.
Sensitivity analysis is based on the sensitivity matrix derived from power flow equations. The sensitivity parameters are determined by the relationship between state variables and control variables [17]. An increase in the sensitivity coefficients usually indicates a power network is moving towards instability. Sensitivity analysis is a simple and clear method that provides a good judge when assessing the static voltage stability status of a power system and identifying the voltage weak points of the network. However, in complex networks, the linear characteristic of the sensitivity index prevents it from clearly reflecting the critical state of a system.
Modal analysis was proposed by Gao et al. [27] to predict voltage collapse in complex power networks. It is based on a linear approximation of the steadystate system model. It involves computing the smallest eigenvalues and associated eigenvectors of the reduced Jacobian matrix obtained from the load flow solution. The eigenvalues are associated with a mode of voltage and load variation which can provide a relative measure of stability. The eigenvalues and eigenvectors are used to calculate bus, branch, and generator participation factors. Bus participation factors identify voltage weak buses and corresponding load regions most susceptible to instability in a power network; branch participation factors indicate branches that cause the various modes to be weak (they are either weak links or heavily loaded), and generator participation factors provide important information regarding the distribution of reactive reserves among all machines in a network so as to maintain an adequate stability margin. Modal analysis can either be approached through QV or PV modal analysis with most literature focusing on QV modal analysis since it produces more accurate results for participation factors. This is because Q is directly related to voltage magnitude [16, 27, 28].
In the recent past, several studies have been carried out on voltage stability analysis and improvement using the conventional methods presented in this paper. PV and QV curves via the voltage stability margins index, sensitivity analysis, and modal analysis are employed in [29–37] in the analysis of voltage stability. However, because of the drawbacks of some of these methods as discussed in Table 3, further studies have been carried out employing the combination of one or more methods so as to get more information and greater degree of accuracy in voltage stability studies. Mao et al. [38–40] use a combination of powervoltage curves and modal analysis in the voltage stability assessment of power systems with the later focusing on a microgrid system. In [41, 42], the voltage stability sensitivity index is used to identify the weakest bus in the IEEE 14 bus system which is further analysed using powervoltage curves to identify the voltage collapse point. Reactive power compensation is also explored on the voltage weak buses. Voltage stability assessment through the integration of modal analysis and sensitivity analysis techniques is further presented in [43–45]. Puppala and Chandrarao [46] investigate voltage stability analysis using CPF, modal analysis, linear regression, and artificial neural networks for different loading conditions. Similarly, in [47–54], the various methods and indices are used in assessing the voltage stability in power systems with renewable energy integration. Sizing and placement of FACTs devices are also investigated by the authors.
Several voltage stability indices (VSIs) have also been suggested over the past few decades, referred to either as bus or line have been proposed for both online and offline studies, although most literature focuses on standard IEEE bus test systems. They are used to determine the proximity of a given operating point to voltage collapse. Some of the VSIs include Voltage Collapse Point Indicator (VCPI) and Line Stability Index, L_{mn} [55], Voltage Stability Index, L [56], Fast Voltage Stability Index (FVSI)) [57], Voltage Reactive Power Index (VQI), and Line Stability Factor (LQP) in [58], and respectively [59]. The authors in [60–63] present a review of the voltage stability indices that have been proposed over the years.
In this paper, a methodology is proposed involving the combination of QV modal analysis, sensitivity analysis (VQ), and powervoltage curves used to assess the static voltage stability status of the Kenyan Power Network. Modal analysis is first used because it gives a system wide perspective on the voltage stability status of the system and identifies the causes of instability in the system, i.e., critical lines and weak buses. VQ sensitivity analysis is further employed to identify buses with large sensitivities in the system. The buses obtained by VQ sensitivity analysis are compared with those obtained through QV modal analysis. Powervoltage curves are then generated for the voltage weak buses identified by the two methods to get a measure of the proximity of these buses to voltage collapse using loading margins and voltage stability margins. The scripts used in generating the PV and QVs employ a predictorcorrector approach in generating the curves. A contingency study is also done to assess the static voltage stability status of the system by tripping one of the lines serving the main load region in the network.
This paper is organised as follows: Section 2 presents a brief literature review on the methods used, Section 3 presents an overview of the model of Kenyan Power Network under study and the methodology used, Section 4 presents the results obtained, the accompanying discussions, and an analysis of the results, and finally, Section 5 provides the conclusions of the study.
2. Literature Review
2.1. PV and QV Curves
PV and QV curves define the steadystate loadability limits which are related to voltage stability. They show the dependence of bus voltage magnitude on the active power load in a bus when increased from zero to its maximum possible value. Figure 1 shows a PV curve. The load is increased gradually and with each increment power flows are calculated until the nose of the PV curve is reached. At the nose point, further load increase provides no feasible operating voltage magnitude and it is the point where voltage collapse usually occurs. This limited load value is called the maximum loadability and the corresponding voltage is called the critical voltage [11, 12].
QV curves are plots of voltage against reactive power injected into a bus. They are generated by converting PQ buses into PV buses with open VAr limits, usually by applying fictitious VAr sources to the buses. QV curves show the reactive power injection or absorption for various scheduled voltages. Figure 2 shows a normalised QV curve. The yaxis shows the amount of reactive power that must be injected or absorbed for a bus to operate at a specific voltage. The operating point is the intersection of the curve with the xaxis. By themselves, QV curves can be used as an index for voltage instability; if the slope of the curve at the intersection is positive, then the system is stable since any additional reactive power will raise the voltage, and vice versa. The point where corresponds to the nose point of the curve. If the minimum point is above the horizontal axis, it indicates a reactive power deficiency and additional reactive power sources are needed to prevent a voltage collapse scenario [11, 12].
Since the stability criterion when using PV and QV curves in voltage stability analysis involves determining the margin between the current operating point and the nose points of the curves, analysis can be done using loading margins and voltage stability margins. For a PV curve, the active power margin is the “distance” of active power from the operating point to nose point of the curve, and for a QV curve, the reactive power margin is the “distance” of reactive power from the operating point to the nose point. These margins are computed using equations (1) and (2), respectively:where and are the values of load active and reactive powers at the nose point of the curves and and are the load values at the operating point.
Voltage stability margin (VSM) of a bus can be defined as the “distance” between the initial voltage operating point until the voltage critical/collapse point. It is obtained from both PV and QV curves [15]. Conventionally, VSM is calculated using the following equation:where is initial bus operating voltage and is the bus voltage at the voltage collapse point.
However, equation (3) fails to account for the load increments that occur during generation of PV and QV curves which can result in two dissimilar buses have same or close to same VSM values. To overcome this, some literatures [64] propose the determination of the voltage stability margin by measuring the hypotenuse distance of the triangle under the PV and QV curves as shown in Figures 3 and 4.
From Figures 3 and 4,where , , , and .
2.2. VQ Sensitivity Analysis
VQ Sensitivity analysis is based on the power flow equations and the Jacobian matrix. The power flow problem can be formulated as follows [12]:where and are the mismatch active and reactive powers, and are the incremental changes in the bus voltage angle and magnitude, respectively, and is the Jacobian matrix.
The diagonal element of the matrix represents the VQ sensitivity of the load bus when considering changes in active and reactive powers. System voltage stability is affected by both P and Q. However, at each operating point, if P is kept constant and voltage stability is evaluated by considering incremental relationship between Q and V. Based on these considerations and letting , equation (6) becomeswhere,
Equation (8) is the linearized power flow equation. is the reduced Jacobian matrix of the system. It is the matrix that directly relates bus reactive power injection and bus voltage magnitude. From equation (7),
The matrix is the reduced VQ Jacobian, and its diagonal element is the sensitivity at . The VQ sensitivity at a bus represents the slope of QV curve at a given operating point. A positive VQ sensitivity is indicative of stable operation; the smaller the sensitivity, the more stable the system. As the stability decreases, the magnitude of the sensitivity increases, becoming infinite at the stability limit. Conversely, a negative sensitivity is indicative of unstable operation. A small negative sensitivity represents a very unstable operation [12].
2.3. QV Modal Analysis
In QV modal analysis, the voltage stability characteristics of a system are identified by computing the eigenvalues and eigenvectors of the power flow reduced Jacobian matrix, [12, 16, 27]. It is represented by eigenvector matrices as shown in the following equation:where is the right eigenvector matrix of , is the left eigenvector matrix of , and is the diagonal eigenvalue matrix of
From equation (11), the reduced Jacobian matrix J_{R} can be expressed as
Equation (11) then becomesorwhere is the eigenvalue, is the column right eigenvector, and is the row left eigenvector of the reduced Jacobian matrix. Each eigenvalue, corresponding right and left eigenvectors, defines a mode of VQ response.
Similarly, from , . Therefore, equations (13) and (14) can be expressed asorwhere is the vector of modal voltage variations and is the vector of modal reactive power variations.
The corresponding modal voltage variation is given by
From equation (17), the stability of mode with respect to reactive power changes is defined by the modal eigenvalue. Large values of suggest small changes in modal voltage for reactive power changes. As the system is stressed, the value of becomes smaller and the modal voltage becomes weaker. If the magnitude of is equal to zero, the corresponding modal voltage collapses since it undergoes infinite changes for reactive power changes. A system is therefore defined as voltage stable if all the eigenvalues of the reduced Jacobian matrix are positive. The bifurcation or voltage stable limit is reached when at least one eigenvalue reaches zero, that is, one or more modal voltages collapses. If any of the eigenvalues is negative, the system is unstable. The magnitudes of eigenvalues provide a relative measure of the proximity of the system to instability [12].
There is no need to evaluate all the eigenvalues of the of a large power system, because once the minimum eigenvalue becomes zero, the system Jacobian matrix becomes singular and voltage instability occurs. On the other hand, evaluating only the minimum eigenvalue of is not sufficient because there is usually more than one mode associated with different parts of the system, and the mode associated with the minimum eigenvalue may not be the most troublesome mode as the system is stressed [12]. Therefore, evaluating 5–10 of the smallest eigenvalues usually gives sufficient information on the voltage stability status of the system.
Bus Participation Factors. In mode , the participation of bus is defined by
The bus participation factor determines the contribution of to the VQ sensitivity at bus . A high value of at bus for mode means this bus is close to voltage instability in this mode. Therefore, the bus participation factors of the critical modes can be used to determine weak buses and predict areas in the system most prone to instability.
Branch Participation Factors. The relative participation of branch in mode is given by the participation factors:
Branches with high participation factors indicate which branches consume the most reactive power in a network for a given incremental change in reactive load. They are either weak links or heavily loaded. Branches with high participation factors in the critical mode are identified as critical branches. Branch participation factors provide valuable information that can be used in contingency selections during contingency studies and how to enhance the transmission lines and redistribute power flows so as to alleviate loading on overloaded lines.
3. Case Study
3.1. Kenyan Power System
The Kenyan power system is rapidly developing, with loading conditions changing relatively fast. Currently, the system has an installed generation capacity of about 3 GW with a peak load demand slightly above 1.9 GW as of the last quarter of 2019. In this paper, a 330bus model of the system was developed and used as the study case. The model consists of 97 substations with 20 of them being generating substations. It has 185 transmission lines transmitting power at 66 kV, 132 kV and 220 kV with the aid 184 2winding transformers and 29 3winding transformers all modelled with load tap changers. The model has a total installed capacity of 2.9268 GW with an effective (grid connected) capacity of 2.14643 GW and a load demand of 1.58859 GW. Figure 5 shows a single line diagram of the model as seen in DigSILENT PowerFactory. A schematic for the same is provided alongside the supplementary data.
In carrying out this study, load flow analysis was performed while neglecting the effect of transformer load tap changers since this is a static voltage stability study. The aim of the load flow study was to determine the voltage profile of the test system. QV modal analysis was then carried out to identify the eigenvalues and eigenvectors of the system. As discussed in Section 2, the minimum eigenvalues which correspond to critical modes of the network which were then selected and evaluated to obtain weak buses in the system and also predict load regions most susceptible to instability. Reactive power loss sensitivity factors were then selected as a criterion for identifying critical lines in the network. The transmission lines with the highest reactive power loss sensitivity factors for each of the critical modes analysed were then determined. VQ sensitivity analysis was then carried out on the system to identify buses that undergo large voltage changes for reactive power imports in the system. A comparison was then made between the buses with the highest participation factors identified in each of the QModes and the sensitivities of these buses.
The load buses identified as having the highest participation factors were further investigated using voltage stability margins discussed in Section 2 to determine the proximity of the buses towards instability. VSM values of these buses were determined through generation of PV and QV curves for the respective load buses. The loading margins of these buses were also computed to determine how much active or reactive power could be absorbed from these buses before voltage instability occurs.
The effect on the static voltage stability of the system with the removal of a line was also investigated. The eigenvalues of the system after contingency were determined. The resultant effect on the voltage stability margins of the load buses with the highest participation factors was also determined.
3.2. Methodology
QV modal analysis was carried out using the load flow sensitivities tool in DigSILENT PowerFactory. By default, DigSILENT PowerFactory provides eigenvalues in the form ranking them in descending order of magnitude. Therefore, the eigenvalues with the largest magnitudes correspond to the critical modes of the system. In this study, the eigenvalues which correspond to the 5 critical modes were evaluated. Evaluation of these modes would suffice in getting a detailed “snapshot” of the voltage stability status of the system. Evaluation of the 5 critical modes involved determining the buses with the highest participation factors in each of the critical modes. These participation factors identify buses participating in making the identified modes weak. These are weak buses in the system. Since each mode can be associated with different parts of the system, load regions most susceptible to instability were also identified. In the identification of the critical lines in the network, reactive power loss sensitivities for the transmission lines were determined for each of the Qmodes. The sensitivities identified the lines that consumed more reactive power for reactive power imports in the system. Similarly, lines that supplied reactive power to the system for reactive power imports in each QMode were also determined.
VQ sensitivity analysis involved determining the diagonal elements of the inverse reduced Jacobian matrix of the system. The VQ sensitivities identify the buses that experience large voltage changes for reactive power variations in the system. These can be categorised as weak buses because a large absorption of reactive power from those buses results in large declines on the bus voltages.
The load buses (rated 3.3 kV to 33 kV) identified as having the highest participation factors in the critical QModes were further analysed using voltage stability margins obtained from the buses’ PV and QV curves. This was to determine the proximity to voltage collapse of these buses. It identified the buses closest to instability among the identified buses in the test network. The PV and QV curves of these buses were generated using DIgSILENT Powerfactory’s UP and UQ command scripts. The UP script generates a PV curve for a selected load bus by incrementally varying the selected load until the power flow does not converge. A QV curve on the other hand is generated by applying a fictitious VAr source at the selected bus, i.e., converting the bus to a PV bus with open VAr limits and then the UQ script is then executed. It automatically generates a QV curve by performing power flow for various set points until the load flow fails to converge. The loading margins of these buses were also determined.
A contingency study is also done to assess the voltage stability of the test system. The N1 contingency was selected as a criterion for the study, in which a line outage was simulated. A line was selected which is of crucial importance to the system but not directly associated with the load regions identified.
4. Results and Discussions
The system has 506 eigenvalues which correspond to different modes of VQ responses in the test system under this particular configuration of the system. The 5 critical modes were determined and analysed. The modes are presented as QModes 1–5 with Qmode 1 being the most critical. Negative eigenvalues were also obtained; they indicate the system is voltage unstable and in the event of a disturbance, voltage instability or even voltage collapse could occur without correct preventive measures. Table 4 presents the eigenvalues corresponding to the most critical modes and the modes with negative eigenvalues. A comparison of the buses with the highest participations in the critical Qmodes evaluated and VQ sensitivities of the buses was done as presented in Table 5.


Maungu, Ndhiwa, Sultan Hamud, Lokichar, and Garissa were identified to have the highest participation factors for QModes 1–5, respectively. The buses associated with QMode 1 and 3 referred to the same load region. Since QModes 1–3 had many buses participating in making them weak, only a sample of the buses was taken, i.e., those having the highest participation factors. In the case of QMode 4 and 5, all the buses are presented. The VQ sensitivities of the network identify Mtito Andei a 3.3 kV load bus as having the highest sensitivity of 0.02078821 p.u/MVar. The lowest VQ sensitivity determined was 0.0001433 p.u/MVar belongs to a 220 kV busbar located in one of the generating stations in the network. All voltagecontrolled buses had 0 sensitivities. Although, a direct comparison cannot be made between the VQ sensitivities and bus participation factors due the nonlinear nature of VQ relationships, the VQ sensitivities identified for the buses with the highest participation factors were the highest in the system. The supplementary data provided present the participation factors of all the buses associated with Qmodes 1–3.
The branches with the highest reactive power loss sensitivities in each of the QModes evaluated are identified are presented in Table 6. The positive reactive power loss sensitivities identify branches that consumed the most reactive power for reactive power imports in the system thus contributing to their respective modes being weak and the negative reactive power loss sensitivity identify branches that supplied reactive power in the network for their respective QModes Failures in some of these lines could lead to low voltage scenarios or voltage collapses in significant parts of the system. The buses identified and lines identified present candidate locations for placement of FACTs devices in the system, especially when improving the voltage stability of the system is the objective.

From Table 5, the load buses (rated 3.3 kV to 33 kV) identified as having the highest participation factors in the critical QModes were further analysed using voltage stability margins obtained from the buses’ PV and QV curves. Figures 6 and 7 present the PV and QV curve of Maungu 3.3 kV load bus identified as having the highest participation factor in QMode 1. From the PV and QV curves, the critical voltages and loading margins were obtained which facilitated the determination of the voltage stability margins. The powervoltage curves for the other selected buses in the table were also plotted and the voltage stability margins determined. Table 7 presents the loading margins and voltage stability margins of these buses.

The results obtained in Table 7 identify Titanium, Chemosit2, Mtito Andei, Lokichar, and Mwingi as being the buses closest to instability in their respective QModes. Amongst the 5 buses, Chemosit2 was the buses closest to instability, having the lowest active and reactive power margins. The voltage of these load bus as identified by load flow analysis was 0.88131 pu, the lowest recorded in the system followed by Titanium identified in QMode 1 which had a voltage of 0.89531 pu which also happened to be the second bus closest to instability after Chemosit. These are buses in need of voltage support in the system and the weakest buses in the test network. Figure 8 presents a comparison between the bus voltages and voltage stability margins and of these weakest buses. The loading margins of these buses, also presented in Table 7, indicate the amount of active and reactive power that could be absorbed from these buses before voltage collapse occurs. It is worth noting that these margins exist in a simulation scenario, in a real situation, voltage instability occurs before the maximum load limit is reached.
The contingency study was also conducted to assess the voltage stability of the test system using N1 criterion. Nairobi North to Dandora line was chosen arbitrarily. This is a double circuit line that serves the major load region in the network. It evacuates power from one of the main generating stations in the country to the most loaded region. In recent years, it has been noted as being one of the most troublesome lines in the network, where an outage in one of the lines usually leads to a cascade of events that eventually lead to a total collapse of the network. One of the Nairobi North to Dandora lines was thus removed from the network and the resultant effect on the static voltage stability of the test network assessed. The effect on the voltage stability margins of the load buses is also noted. Table 8 presents the data obtained with the removal of the line from the network.

From Table 8, it is observed that load regions 1 and 3 associated with QModes 1and 3 are the most affected with the outage of one of the Nairobi North to Dandora lines. Maungu and Sultan Hamud identified in QModes 1 and 3, respectively, have the greatest percentage decrease in active powervoltage stability margins. In terms of reactive power, Jomvu and Maungu in QMode 1 and Mtito Andei in QMode 3 have the greatest decrease in voltage stability margins. A recalculation of the system eigenvalues identifies the magnitude of the smallest eigenvalue, as 0.12195, an increase from 0.11958. This coupled with the decreased voltage stability margins in the studied load buses indicates the system moves even closer to voltage instability. This is attributed to the fact that, with removal of the line, the effective power transferred decreases in the network.
5. Conclusion
In this paper, a methodology is proposed for carrying out static voltage stability assessment in bulk power systems. The methodology is employed in the assessment of the static voltage stability status of the Kenyan Power Network based on a model of the system as of 2019 has been carried out. QV Modal analysis, a practical based approach towards static voltage stability assessment of bulk power systems, has been effectively employed to determine load regions most susceptible to instability in the network. Maungu, Ndhiwa, Sultan Hamud, Lokichar, and Garissa have been identified as load buses contributing the most towards voltage instability in the various load regions identified. The bus VQ sensitivities of the system have also been determined to identify buses that have the largest reactive power changes for reactive power imports in the network. QV modal analysis and VQ sensitivity analysis complement each other as the buses having the highest participation factors in their respective QModes also have the highest VQ sensitivities in the network. The critical lines in the network, based on the reactive power loss sensitivities for the QModes analysed, have also be identified. These are lines that consume or supply the most reactive power for reactive power imports in the system.
Further, the load buses with the highest participation factors in the various load regions have been analysed using voltage stability margins with an aim of determining the proximity of these buses towards voltage instability. Titanium, Chemosit2, Mtito Andei, Lokichar, and Mwingi are identified as the buses closest to instability in the network, since they have the lowest voltage stability margins, with Chemosit2 having the lowest VSM values among the 5 buses. These are buses that require voltage support in the network. This study has provided a useful approach towards static voltage stability assessment of bulk power systems, especially in identifying regions most susceptible instability, the contributing factors, and in the identification of buses that are in dire need of voltage support in the system. The proposed methodology can be particularly useful especially when it comes to carrying out studies aimed at maintaining the voltage security of a system.
Data Availability
Research data used to support the findings of this study are available from the corresponding authors upon request.
Conflicts of Interest
The authors declare they have no conflicts of interest.
Supplementary Materials
Tables SISIII present the participation factors in QModes 1–3. (Supplementary Materials)
References
 M. Glavic and D. Novosel, Voltage Stability Monitoring, Instability Detection and Control, Quanta Technology, Raleigh, NC, USA, 2011.
 V. Ajjarapu and B. Lee, “Bibliography on voltage stability,” IEEE Transactions on Power Systems, vol. 13, no. 1, pp. 115–125, 1998. View at: Publisher Site  Google Scholar
 North American Electric Reliability Council, System Disturbances: Review of Selected 1996 Electric System Disturbances in North America, North American Electric Reliability Council, Princeton, NJ, USA, 1996.
 J. M. Ordacgi Filho, Brazilian Blackout 2009: Blackout Watch Protection, Automation and Control World, 2010, http://www.pacw.org/fileadmin/doc/MarchIssue2010/Brazilian_Blackout_march_2010.pdf.
 Y. Hain and I. Schweitzer, “Analysis of the power blackout on june 8, 1995 in the Israel electric corporation,” IEEE Transactions on Power Systems, vol. 4, no. 12, pp. 1752–1758, 1997. View at: Publisher Site  Google Scholar
 P. Pourbiek, P. Kundur, and C. W. Taylor, “The anatomy of power grid blaclout,” IEEE Power and Energy Magazine, vol. 5, no. 4, pp. 22–29, 2006. View at: Google Scholar
 L. D. Vargas, R. D. Miranda, and V. H. Quintana, “Voltage collapse scenario in the Chilean interconnected system,” IEEE Transactions on Power Systems, vol. 4, no. 14, pp. 1415–1421, 1999. View at: Publisher Site  Google Scholar
 C. D. Vournas, V. C. Nikolaidis, and A. A. Tassoulis, “Postmortem analysis and data validation in the wake of the 2004 athens blackout,” IEEE Transactions on Power Systems, vol. 3, no. 21, pp. 1331–1339, 2006. View at: Publisher Site  Google Scholar
 D. M. Kosterev, C. W. Taylor, and W. A. Mittelsadt, “Model validation for the august 10, 1996 WSCC system outage,” IEEE Transactions on Power Systems, vol. 14, no. 3, pp. 967–979, 1996. View at: Google Scholar
 J. G. Ndirangu, J. N. Nderu, A. M. Muhia, and C. M. Maina, “Power quality challenges and mitigation meaures in grid integration of wind energy conversion system,” in Proceedings of the IEEE International Energy Conference (ENERGYCON), Limassol, Cyprus, June 2018. View at: Google Scholar
 T. V. Custem and C. Vournas, Voltage Stability Analysis in Electrical Power Systems, Springer Science, New York, NY, USA, 1998.
 P. Kundur, Power System Stability and Control, McGrawHill, New York, NY, USA, 1994.
 IEEE Committee Report, Voltage Stability of Power Systems: Concepts, Analytical Tools and Industrial Experience, IEEE/PES, Piscataway, NJ, USA, 1990.
 IEEE/CIGRE Joint Task Force on Stability Terms and Definitions, “Definition and classification of power system stability,” IEEE Transactions on Power Systems, vol. 19, no. 2, pp. 1387–1401, 2004. View at: Google Scholar
 C. W. Taylor, Power System Voltage Stability, McGrawHill, New York, NY, USA, 1993.
 B. Gao, G. K. Morison, and P. Kundur, “Volatge stability analysis using static and dynamic approaches,” IEEE Transactions on Power Systems, vol. 8, no. 3, pp. 1159–1171, 1993. View at: Google Scholar
 M. Youjie, L. Shaofeng, Z. Xuesong, and G. Zhiqiang, “Review analysis of voltage stability in power system,” in Proceedings of the IEEE International Conference on Mechatronics and Automations, Takamatsu, Japan, August 2017. View at: Google Scholar
 X. Y. Zhao, X. B. Zhang, and X. L. Su, “Voltage stability and bifurcation theory in power systems,” Transactions of China Electrotechnical Society, vol. 23, pp. 87–95, 2008. View at: Google Scholar
 A. S. Mendoza, F. Esquivel, and C. R. Becerril, “A numerical study on the effect of degenerate hopf bifurcations on the voltage stability of power systems,” Electrical Power Systems Research, vol. 101, no. 8, pp. 102–109, 2013. View at: Google Scholar
 Y. J. Ma, X. S. Li, and X. S. Zhou, “The comments on dynamic bifurcation of voltage stability in power systems,” ICIE, vol. 355, IEEE, 2010. View at: Google Scholar
 T. J. Overbye and C. L. Demarco, “Improved technique for power system voltage stability assesment using energy methods,” IEEE Transactions on Power Systems, vol. 9, no. 1, pp. 1446–1452, 1991. View at: Publisher Site  Google Scholar
 T. J. Overbye, I. Dobson, and C. L. DeMarco, “QV curve interpretations of energy measures for voltage security,” IEEE Transactions on Power Systems, vol. 9, no. 1, pp. 331–340, 1994. View at: Publisher Site  Google Scholar
 C. A. Canizares, I. Dobson, T. Van Cutsem et al., 2003, Voltage stability assesment: concepts, practices and tools,” IEEEPES power systems stability subcommittee special publication.
 V. Ajjarapu and C. Christy, “The continuation power flow: a tool for steady state voltage stability analysis,” IEEE Transactions on Power Systems, vol. 8, no. 1, pp. 415–423, 1992. View at: Google Scholar
 V. Ajjarapu, Computational Techniques for Voltage Stability Assessment and Control, Springer Science, New York, NY, USA, 2006.
 P.A. Lof, G. Andersson, and D. J. Hill, “Voltage stability indices for stressed power systems,” IEEE Transactions on Power Systems, vol. 8, no. 1, pp. 326–335, 1993. View at: Publisher Site  Google Scholar
 B. Gao, G. K. Morison, and P. Kundur, “Voltage stability evaluation using modal analysis,” IEEE Transactions on Power Systems, vol. 7, no. 4, pp. 1529–1542, 1992. View at: Publisher Site  Google Scholar
 P. Khampariya and B. Telang, “Voltage stability analysis using modal analysis,” International Journal of Scientific Research, Engineering and Technology (IJSRET), vol. 4, no. 4, pp. 408–411, 2015. View at: Google Scholar
 M. Sulaiman, A. F. Bujal, and N. R. Bujah, “Voltage instability analysis on PV and QV curves for radialtype and meshtype electrical power networks,” International Review of Electrical Engineering (IREE), vol. 10, no. 1, pp. 109–115, 2015. View at: Publisher Site  Google Scholar
 R. Bujal, E. A. Hassan, and M. Sulaiman, “Analysis of voltage stability in power systems,” in Proceedings of the International Conference of Engineering Technology and Technoprenuership, Kuala Lumpur, Malaysia, August 2014. View at: Google Scholar
 A. Parma and G. Radman, “A new approach for estimating voltage collapse point based on quadratic approximation of PV curves,” Electric Power Systems Research, vol. 79, pp. 653–659, 2009. View at: Publisher Site  Google Scholar
 F. O. Enemuoh, J. C. Onuegbu, and E. A. Anazia, “Modal based analysis and evaluation of voltage stability of bulk power systems,” International Journal of Engineering Research and Development, vol. 6, no. 12, pp. 71–79, 2013. View at: Google Scholar
 B. Maya and R. Sunitha, “Voltage stability assessment of large power systems using modal analysis,” in Proceedings of the International Conference on Emerging Trends and Electrical Computer Technology, Nagercoil, Tamilnadu, India, March 2011. View at: Google Scholar
 C. Mitali and S. Patra, “Voltage stability analysis using conventional methods,” in Proceedings of the International Conference on Signal Processing, Communication, Power and Embedded Systems (SCOPES), Odisha, India, October 2016. View at: Google Scholar
 R. Toma and M. Gavrilas, “Volatge stability assessment based on load flow sensitivity method,” in Proceedings of the International Conference and Exposition on Electrical and Power Engineering, Romania, Balkans, October 2018. View at: Google Scholar
 M. M. A. Seedahmed, S. A. K. Mokred, and G. Kamara, “Voltage stability estimation for complex power system based on modal analysis technique,” in Proceedings of the 6th International Conference on Signal Processing and Intergrated Networks, Noida, India, March 2019. View at: Google Scholar
 J. Patino, C. A. Ramires, and J. Espinoza, “Modal analysis for power system benchmark with topographical changes,” Communications in Computer and Information Science, vol. 1052, 2019. View at: Google Scholar
 W.X. Mao, Q. Yang, M.Y. Yang, S. Zhang, F. Gao, and Y.T. Song, “Static voltage stability analysis of power system using the combination of PV curve and modal analysis,” in Proceedings of the International Power, Electronics and Materials Engineering Conference, Dalian, China, May 2015. View at: Google Scholar
 Y. A. Jabri, N. Houseinzadeh, R. A. Abri, and A. A. Hinai, “Voltage stability assessement of a microgrid,” in Proceedings of the IEEE 8th GCC Conference and Exhibition, Muscat, Oman, February 2015. View at: Google Scholar
 B. Ombuki and R. Orenge, “Steadystate voltage stability assessment of the Kenyan high voltage transmission network,” in Proceedings of the 2019 IEEE AFRICON, Accra, Ghana, September 2019. View at: Google Scholar
 B. Sourav, G. Chandramouli, and C. Surmana, “Voltage stability margin and maximum loading point of a multibus system before and after compensation,” International Journal of Engineering Research and Development, vol. 5, no. 7, pp. 30–35, 2013. View at: Google Scholar
 F. Mahmood, S. Ahmad, G. Mukita, S. Razwan, M. Maruf, and F. Abatsh, “Weakest location exploration in IEEE14 bus system for voltage stability improvement using STATCOM, synchronous condenser and static capcitor,” in Proceedings of the IEEE Conference on Electrical, Computer and Communication, Cox’s Bazar, Bangladesh, February 2017. View at: Google Scholar
 M. K. Jalboub, H. S. Rajamani, R. A. AbdiAlhameed, and A. M. Ihbal, “Weakest bus identification based on modal analysis and singular value decomposition techniques,” in Proceedings of the 1st International Conference on Energy, Power and Control (EPCIQ), Basrah, Iraq, December 2010. View at: Google Scholar
 M. A. Djari, L. Benasal, and W. Rahnouni, “Voltage stability assessment using VQ sensitivity and modal analyses methods,” in Proceedings of the 5th International Conference on Electrical Engineering, Boumerdes, Algeria, October 2017. View at: Google Scholar
 J. G. Mercado, V. Reyes, H. T. Hernandz, and J. Z. Ayala, “Voltage stability assessment by modal analysis and load flow linear sensitivity techniques,” in Proceedings of the IEEE International Autumn Meeting on Power Electronics and Computing, Ixtapa, Mexico, November 2018. View at: Google Scholar
 V. Puppala and P. Chandrarao, “Analysis of continuous power flow method, model analysis, linear regression and ANN for voltage stability assessment for different loading conditions,” Procedia Computer Science, vol. 47, pp. 168–178, 2015. View at: Publisher Site  Google Scholar
 A. B. Almeida, E. V. De Lorenci, R. C. Leme, A. C. Z. Souza, B. I. L. Lopes, and K. Lo, “Probabilistic voltage stability assessment considering renewable sources with the help of the PV and QV curves,” IET Renewable Power Generation, vol. 7, no. 5, pp. 521–530, 2013. View at: Publisher Site  Google Scholar
 V. N. Sewdien, R. Preece, J. R. Torres, and M. Meijder, “Evaluation of PV and QV based voltage stability analyses in converter dominated power systems,” in Proceedings of the IEEE PES AsiaPacific Power Energy Engineering, Kota Kinabalu, Malaysia, October 2018. View at: Google Scholar
 N. Manjul and M. S. Rawat, “PV/QV curve based optimal placement of statc VAR system in power network using DigSILENT PowerFactory,” in Proceedings of the IEEE 8th Power India International Conference, Kurukshetra, India, December 2018. View at: Google Scholar
 X. Liang, X. Yan, and N. Khan, “Measurement based characteristic curves at point of interconnection of wind farms,” in Proceedings of the IEEE Electrical Power and Energy Conference, Toronto, Canada, October 2018. View at: Google Scholar
 R. Kumal, A. Mittal, N. Sharma, A. Kumal, and I. Duggal, “PV and QV curve analysis using series and shunt compensation,” in Proceedings of the IEEE 9th Power India International Conference, Sonepat, India, March 2020. View at: Google Scholar
 I. C. Guradin, Z. Muslimin, A. M. Ilyas, and A. Siswan, “Comparison of voltage stability index before and after wind turbine penetrated to sulsetrebar interconnection power system using modal analysis method,” Material Science and Engineering Conference Series, vol. 875, 2020. View at: Google Scholar
 Y. Burak, “Investigation with modal analysis of effects of high PV penetration on power system voltage stability,” International Journal of Energy and Smart Grid, vol. 2, no. 1, pp. 17–26, 2017. View at: Google Scholar
 A. Phadke, S. K. Banisal, and K. Rheman, “A comparison of voltage stability indices for placing shunt FACTs controllers,” in Proceedings of the 1st International Conference on Emerging Trends in Engineering and Technology, Thrissur, Kerala, India, January 2018. View at: Google Scholar
 M. Moghavemmi and F. Omar, “Technique for contigency monitoring and voltage collapse prediction,” IEE Proceedings on Generation, Transmission and Distribution, vol. 145, pp. 634–640, 1998. View at: Google Scholar
 P. Kessel and H. Glavitsch, “Estimating the voltage stability of a power system,” IEEE Transactions on Power Delivery, vol. 1, no. 3, pp. 346–354, 1986. View at: Publisher Site  Google Scholar
 I. Musirin and A. Rahman, “Estimating maximum loadability for weak bus identification using FVSI,” IEEE Power Engineering Review, vol. 22, no. 11, pp. 50–52, 2002. View at: Google Scholar
 M. W. Mustafa and F. A. Althowibi, “Voltage stability calculations in power transmission lines: indication and allocations,” in Proceedings of the IEEE International Conference on Power and Energy, Kuala Lumpur, Malaysia, December 2010. View at: Google Scholar
 A. Mohammed, G. Jasmon, and S. Yusoff, “A static voltage collapse indicator using line stability factors,” Journal of Industrial Technology, vol. 7, pp. 73–85, 1989. View at: Google Scholar
 C. Reis and M. Barbosa, “A comparison of voltage stability indices,” in Proceedings of the IEEE Mediteranean Electrotechnical Conference, Malaga, Spain, May 2006. View at: Google Scholar
 E. Gholipour, J. Modarresi, and A. Khodabakhshian, “A comprehensive review of voltage stability indices,” Renewable Energy and Sustainable Energy Review, vol. 63, pp. 1–12, 2016. View at: Google Scholar
 S. S. Danish, T. Senjyu, S. S. Danish, R. N. Sabory, K. Narayan, and P. Mandal, “A recap on volatge stability indices in the past three decades,” Energies, vol. 12, no. 8, 2019. View at: Google Scholar
 H. Zaheb, M. S. S. Danish, S. Tomonobu et al., “A contemporary novel classification of voltage stability indices,” Applied Sciences, vol. 10, no. 5, p. 1639, 2020. View at: Google Scholar
 R. Omar, M. Sulaiman, A. Kadir, and A. F. Mohamad, “Determining voltage stability margin values by measuring hypotenuse under PV and QV curves,” International Journal of Electrical Engineering and Applied Sciences, vol. 1, no. 1, 2018. View at: Google Scholar
Copyright
Copyright © 2021 Ombuki Mogaka et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.