Static Voltage Stability Analysis Using PSAT
Static Voltage Stability Analysis Using PSAT
ABSTRACT: Power System may lose stability and cause blackout when operating under loaded condition. In order to
protect the power system it becomes essential to predict and identify the voltage instability. This paper presents the
analysis of static voltage stability using continuation power flow method considering the power system contingencies
and its effect on Mega Watt Margin (MWM) and maximum loading point (MLP). The weak bus in the system is
identified by using contingencies ranking. On the basis of MWM percentage decrease, ranking of the system is done.
IEEE 6 and IEEE 30 are simulated using PSAT to verify the results.
KEYWORDS: Voltage collapse, voltage stability, megawatt margin, maximum loading point.
I. INTRODUCTION
The increasing demand of electricity results into overloaded condition of the system. The overloaded condition leads
the system to operate in unstable condition. In unstable condition voltages of the buses decreases and can cause voltage
collapse [1]. It is a phenomenon that causes severe changes in the power system hence voltage stability analysis
becomes necessary. Analysis of voltage stability is also essential from the power system protection point of view [2-3].
Continuation power flow (CPF) is the main analysis tool used for the IEEE 6 and IEEE 30 bus system. CPF technique
involves the use of predictor and corrector step for analysing static voltage stability [4]. In this paper, for the analysis of
static voltage stability maximum loading point (MLP) and megawatt margin (MWM) are used. Contingency analysis is
important factor while considering the power system protection [5]. Operation of system under no contingency gives
maximum loading point and maximum megawatt margin. Continuation power flow method is implemented on IEEE 6
and IEEE 30 bus systems for the computation of MLP and its corresponding decrease in MWM.
Contingency ranking is obtained by using percentage decrease in MWM [6]. On the basis of ranking system effect of
contingencies on the test system and proximity towards voltage collapse can be determined. As a result, weak bus in
the system can be determined and corrective action can be taken out before system collapse. Contingency identification
can be done by using channel component theorem [7]. Contingency ranking is done by using artificial neural method,
using line stability index method and performance index method [8-10]. Exact and precise method can be used to rank
the contingency of power system [11]. Analysis of contingency ranking gives the idea about corrective measure that
can be taken before system gets collapse [12]. FACTS devices can be used to enhance voltage stability during
contingency [13].
Static voltage stability includes slowly occurring changes in the system which results into the insufficient supply of
reactive power and the voltage drop in the buses. This can be determined from V – λ curve which gives proximity
towards voltage collapse. The V – λ curve shows that increase in loading factor (λ) reduces the voltage at receiving end
(V) resulting in disturbance of power handling capacity. The main reason behind the voltage drop is the presence of
contingencies in the system. Contingency may result into lack of reactive power in the system and smaller MWM
causing less voltage in the buses of the system. For the analysis of static voltage stability it becomes essential to rank
the contingency on the basis of MLP and percentage decrease in MWM.
Continuation power flow method is used to rank power system contingency. In the contingency analysis, maximum
loading point and corresponding decrease in percentage of MWM is computed. Contingencies ranking is done by
arranging the MLP in ascending order and its corresponding MWM decrease in decreasing order. Contingency having
lower MLP and higher percentage change in MWM is ranked higher than other contingencies. The values of MMWM
and MWM are calculated as:
MMWM = P imax-P base
MWM = P i+1max - P base
where, P max is maximum load active power corresponding with MLP and P base is base load active power.
The MWM decrease percent is calculated as:
MWM decrease percent = 100 × [1-( )]
The number of contingencies in the power system depends upon the number of elements subjected to failure. For the
event number of L with NCL: L=0,1,2,…,
!
NCL = ! ( )!
NC0 means no element is exposed to failure in the system called as zero level contingency. First level contingency
NC1means only one unique element is subjected to failure. Usually, first and second number of contingency is taken
into account. In this paper zero and first level contingency are considered.
NC0,1= 1+N.
The line and bus data for IEEE 6 and 30 bus systems are taken from case6ww and case 30 in Matpower 5.1. The P and
Q load powers are voltage independent and its changes can be assumed as:
PL= PL0(1+λ) &QL=QL0(1+λ)
Where, PL0 and QL0 are the active and reactive base loads and PL and QLare the active and reactive loads at bus L for
the operating point is given by λ.
A. IEEE 6 bus system
Continuation power flow is executed on IEEE 6 bus system with no contingency, maximum loading point λmaxis
obtained as 3.3208.The base load active power which is the summation of total active power at the buses during no
contingency is computed as 2.94p.u. The maximum load active power is the summation of total active power at the
buses during contingency is 6.9735p.u.
(01) SINGLE GENERATION UNIT OUTAGES
In single generation unit outage contingency the values obtained are tabulated in table 1. In generation unit outage the
bus number with lowest voltage magnitude is specified with loading factor, real and reactive power in the table 1.
Generation unit Bus no. with lowest Magnitude of lowest voltage in λmax Pload Qload
outage voltage magnitude MLP(p.u) (p.u) (p.u) (p.u)
The results obtained of MWM during zero level contingency are 4.0335. The values obtained of MWM and MWM %is
tabulated in the table 2 for generation unit outages.
The values obtained of percentage decrease in MWM are tabulated in table 3. In generation unit outage for bus 3
having λ=2.2179 and MWM % decrease= 57.42 is identified as critical contingency among the remaining contingency
of generation unit outage and it is tabulated in table 3
.
Table 03 CONTINGENCIES RANKING OF SINGLE GENERATION UNIT OUTAGES
Rank of the contingency Generation unit outage λmax(p.u) Percentage decrease in MWM
1 Bus 3 2.2179 57.42
2 Bus 2 2.2893 53.70
In generation unit outage connected tobus 3 ismarked as critical contingency in the system. Similarly other generator
outages are ranked on the basis of MWM percentage decrease, given in table 3.
The result obtained of MWM during zero level contingency is 4.0335. The values obtained of MWM and MWM% is
tabulated in the table 5.
In line unit outage, line 2 connected between bus 3 and 6 having λ= 1.7218 and MWM % decrease= 83.24 is identified
as critical contingency among the remaining contingencies of line unit outages and it is tabulated in table 6.
In line unit outage, line 2 connected to bus 3 and 6 ismarked as critical contingency in the system. Similarly other lines
are ranked on the basis of MWM percentage decrease, given in table 6.
The results obtained of MWM during zero level contingency are 7.7895. The values obtained of MWM and MWM %is
tabulated in the table 8 for generation unit outages.
The values obtained of percentage decrease in MWM are tabulated in table 9. In generation unit outage connected to
bus 27 having λ= 2.7428 and MWM % decrease= 67.39 is identified as critical contingency among the remaining
contingencies of generation unit outage.
In generation unit outage connected to bus 27 ismarked as critical contingency in the system. Similarly other generator
outages are ranked on the basis of MWM percentage decrease, given in table 3.
(02)LINE OUTAGES
In line unit outage values obtained are tabulated in table 10. The bus number with lowest voltage magnitude is
specified with loading factor in the table.
Outage line Bus no with lowest voltage Magnitude of lowest voltage in λmax
magnitude MLP (p.u) (p.u)
No contingency 8 0.9649 5.5171
Line 1 8 0.6196 4.8899
Line 2 8 0.5042 4.9981
Line 3 8 0.5023 4.8611
Line 4 8 0.5054 5.0676
Line 5 7 0.5129 4.5952
Line 6 8 0.5102 4.5042
Line 7 8 0.5208 4.3307
Line 8 7 0.4879 4.541
Line 9 7 0.5701 3.7153
Line 10 8 0.5069 2.0345
Line 11 8 0.5089 5.2937
Line 12 8 0.4877 5.4031
Line 13 8 0.5029 5.5172
Line 14 8 0.5073 5.2941
Line 15 8 0.5008 5.4244
Line 17 8 0.5003 5.4373
Line 18 8 0.5053 5.3496
Line 19 8 0.5166 5.3469
Line 20 8 0.5026 5.5146
Line 21 8 0.5042 5.4483
Line 22 18 0.4853 4.6828
Line 23 8 0.5012 5.3736
Line 24 19 0.4735 4.1476
Line 25 20 0.4886 3.5987
Line 26 17 0.4733 3.8712
Line 27 8 0.4980 5.2389
Line 28 8 0.5004 5.3691
Line 29 21 0.5015 4.5867
Line 30 8 0.5205 5.2848
Line 31 8 0.5013 5.4947
Line 32 8 0.5149 5.3756
Line 33 8 0.4967 5.495
Line 35 8 0.5116 5.3685
Line 36 8 0.4880 5.1482
Line 37 29 0.5334 3.2365
Line 38 30 0.6899 2.4635
Line 39 30 0.5563 4.2177
Line 40 8 0.4979 5.2915
Line 41 8 0.5035 5.498
The result obtained of MWM during zero level contingency is 7.7895. The values obtained of MWM and MWM% is
tabulated in table 11.
In line unit outage, line 10 connected between bus 6 and 8 having λ= 2.0345 and MWM % decrease = 84.59 is
identified as critical contingency among the remaining contingencies of line unitoutages and it is tabulated in table 12.
In line unit outage, line 10 connected to bus 6 and 8 ismarked as critical contingency in the system. Similarly other
lines are ranked on the basis of MWM percentage decrease, given in table 12.
IV.CONCLUSION
In this paper, we have analysed static voltage stability by using power system contingencies like generator outage and
line outage. Ranking of contingencies is done on the basis of percentage decrease in MWM and MLP by using
continuation power flow method in PSAT. The results obtained during contingencies shows that reduction in voltage
drop in the buses of the system, hence it gives reduce value of MWM and percentage decrease in MWM. Contingency
having higher percentage decrease in MWM is placed at the top. The weak bus in the system can be found out by this
method and corrective action can take place before system gets collapse.
REFERENCES
[01] PrabhaKundur, “Power System Stability and Control”, New York: McGraw-Hill, 1994.
[02] K. B Hemalatha, S.T Jayachristra, “Power System Analysis”, SCI TECH Publication, 2010.
[03] R. K. Rajput, “A Textbook of Power System Engineering”, Laxmi publication, 2012.
[04] MostafaAlinezhad and MehrdadAhmadiKamarposhti, “Static Voltage Stability Assessment Considering the Power System Contingencies using
Continuation Power Flow Method”, InternationalJournal of Electrical and Electronics Engineering, 2010.
[05] ArthitSode-Yome, NadarajahMithulananthan and Kwang Y. Lee, “A Maximum Loading Margin Method for Static Voltage Stability in Power
Systems”, IEEE Transactions on Power Systems, Vol. 21, Issue 2, May 2006.
[06] KusumVerma, K.R. Niazi, “Supervised Learning Approach to on Line Contingency Screening And Ranking In Power System,” Electrical
Power and Energy Systems, 2012.
[07] IrajRahimiPordanjani, YunfeiWang and WilsunXu, “Identification of Critical Components for Voltage Stability Assessment Using Channel
Components Transform” IEEE Transactions on Smart Grid, Vol. 4, Issue 2, pp: 1122-1132, June 2013.
[08] Sidhu T.S., Cui L., “Contingency screening for steady-state security analysis by using FFT and artificial neural networks”, IEEE Transaction on
Power System, Vol. 15, Issue 1, pp. 421–26, 2000.
[09] Subramani C., SubhransuSekhar Dash,Vivek Kumar and Harish Kiran, “Implementation of Line Stability Index for Contingency Analysis and
Screening in Power Systems”, Journal of Computer Science Science Publications, pp: 585-590, 2012.
[10] Ali Abdulwahhaband Abdulrazzaq, “Contingency ranking of power systems using a performance index”, International Research Journal of
Engineering and Technology, Vol 02, Issue 02, May 2015.
[11] Udaykumar and Udaykuma, “Contingency Ranking in Modern Power System by Exact and Precise Method” International Journal of Innovative
Research in Electrical, Electronics, Instrumentation and Control Engineering, Vol. 3, Issue 5, May 2015.
[12]MoussaReda Mansour, Luís Fernando Costa Alberto and Rodrigo Andrade Ramos, “Preventive Control Design for Voltage Stability Considering
Multiple Critical Contingencies” IEEE Transactions on Power systems, Vol. 31, Issue 2, pp: 1517-1525, March 2016.
[13]M Reddy Prasanna and C N Arpitha, “Voltage Stability Enhancement in Contingency Conditions Using Shunt Facts Devices” International
journal of electric and electronic engineering and telecommunication, Vol. 2, Issue 1, 2001.