Low frequency oscillations (LFOs) are load angle oscillations that have a frequency between 0.12.0 Hz. Power system stabilizers (PSSs) are very effective controllers in improvement of the damping of LFOs. PSSs are designed by linearized models of the power system. This paper presents a new model of the power system that has the advantages of the Single Machine Infinite Bus (SMIB) system and the multi machine power system. This model is named a single machine normalbus (SMNB). The equations that describe the proposed model have been linearized and a lead PSS has been designed. Then, particle swarm optimization technique (PSO) is employed to search for optimum PSS parameters. To analysis performance of PSS that has been designed based on the proposed model, a few tests have been implemented. The results show that designed PSS has an excellent capability in enhancing extremely the dynamic stability of power systems and also maintain coordination between PSSs.
1. Introduction
Since 1960, low frequency oscillations have been observed when large power systems are interconnected by relatively weak tie lines. These oscillations may sustain and grow to cause system separation if no adequate damping is available
[1]
.
Two types of oscillation phenomena can occur on the present power system. One is where the oscillation of one generator at a specific power plant has an influence on the system. This type of oscillation is called localmode oscillation and its behavior is mainly limited to the local area in the vicinity of the power plant. It has been known that the local oscillation is likely to occur when power is transmitted over longdistance transmission lines from a power plant at a remote location. This type of system can be accurately modeled using the SMIB system model
[2]
.
The other case has been known as interarea mode oscillation. This is the case where the lowfrequency oscillation is maintained between sets of generators in an interconnected power system. The simplest type of lowfrequency oscillation in the interarea mode is between two interconnected areas. The interarea mode oscillation has a long history. It has been observed in the tieline connecting the large Pacific Southwest and the Pacific Northwest in the United States. It has also been observed on the tieline connecting the northern Midwest and Canada
[3
,
4]
.
In order to avoid LFO, supplementary stabilizing signal have been proposed in the excitation systems through lead/lag power system stabilizers
[5]
or PI  PID PSSs
[6]
. The calculation of the parameters of these PSSs is based on the linearized model of the power system around a nominal operating point. PSSs are the most effective devices used to damp LFOs. For many years, conventional PSS (CPSS) have been widely used in the industry because of their simplicity
[7]
.
Modern control theory has been applied to the PSS design in recent years. These include optimal control, adaptive control; variable structure control and H
_{∞}
control
[8]
.
There have been proposed two power system linear models, which are the SMIB system and the multi machine power system. These models are used for designing PSSs. The SMIB system comprises a generator, a transmission line and an infinitebus which has been substituted instead of the rest of the power system and multi machine comprises all machines in the power system. The SMIB system is a simple model but in this model, interarea mode oscillations are not considered. Therefore, the designed PSS (especially CPSS) according to this model lead to interfere among PSSs, because the effects of the rest of the power system have been ignored. In practice, there are different kinds of power plants connected to the power system. Hence, LFOs have different frequency. For this reason, while PSSs is being designed, it is obligatory to consider the influences of the rest of the power system. Another model of the power system (the multi machine power system) comprises all machines that deigning PSS based on it is very precise and complete. Although this model is very complex. As a result, it is essential to have a model as simple as the SMIB system and as precise as the multi machine power system. For reaching this purpose, this paper proposed a new model of the power system that has mentioned characteristics. According to this model, a lead PSS has been designed where its parameters have been optimized by PSO algorithm. The parameters have been optimized for nominal load when the input power of a generator has been changed suddenly and the effects of the rest of the power system that have not been modeled by a SMIB system have been applied by two inputs. The results show that designed PSS has an excellent capability in improving extremely the dynamic stability of power systems and also maintain coordination between PSSs.
2. Power System Models for Studying Low Frequency Oscillations
This section presents a brief review of the different kinds of the power system models that has been recommended for studying low frequency oscillations.
 2.1 The single machine infinite bus system
Fig. 1
shows a power system that shows the SMIB system. The infinite bus depicts the thevenin equivalent of a large interconnected power system. The nonlinear equations that describe the generator and excitation system have been represented in following equations:
The above equations can be linearized for small oscillation around an operating condition
[9

11]
and be cast in the block diagram shown in
Fig. 2
. The linearized model of the power system as shown in
Fig. 1
is given as follows:
The Single machine infinite bus model (SMIB)
Block diagram of linearized model of the Single machine infinite bus model
 2.2 The Multimachine power system
Reference
[11]
used linearized model of the SMIB for every machine in the power system and has been proposed a new model that is called the multimachine power system. Block diagram of this model is shown in
Fig. 3
.
Block diagram of linearized model of the multimachine power model
3. Power System Stabilizers and Reasons for Loss of Their Coordination
Low frequency oscillations are one of the undesirable phenomena in the power systems. This phenomenon occurs when an electromechanical mode of generators is excited. In the power systems, switching, inserting loads and so on lead to excite this mode. Practically, prevention of an electromechanical mode of generators is impossible, but there is only one way to deal with it, which is increasing the low frequency oscillations damping. For achieving to this aim, It is sufficient to add D
_{e}
( D
_{e}
> 0 ) in parallel with D
_{m}
(
Fig. 4
). From a practical point of view, added this block means supplying a torque opposite direction of angular velocity variation that this act is done by PSSs. The basic function of PSS is added to damp the generator rotor oscillations by controlling its excitation by using auxiliary stabilizing signal(s).
Based on the Automatic Voltage Regulator (AVR) and using speed deviation, power deviation or frequency deviation as additional control signals, PSS is designed to introduce an additional torque coaxial with rotational speed deviation, so that it can be increased lowfrequency oscillation damping and enhance the dynamic stability of the power system
[12]
.
Fig. 5
shows the block diagram of the SMIB system along with classic PSS (CPSS). CPSS makes phase lead for counteracting phase lag of the electrical loop. This means that CPSS is designed so that in frequency of the electromechanical mode of the generator, its phase lead would be corresponded with the phase lag of the electrical loop. According to this design, Band width of CPSS is very narrow and only in frequency of the electromechanical mode of the generator, its phase lead be corresponded with the phase lag of the electrical loop. On the other hand, additional torque that supplying for more damping, only in that frequency would be the opposite direction of angular velocity variation. In addition, there are different kinds of power plants connected to the power system, such as fossil fuel, hydro and nuclear power plants and generators that have different characteristics. Therefore, generators have different electromechanical frequency mode relative to each other. If electromechanical mode of a generator (such as generator B
Fig. 6
) in the power system is excited, it will lead to produce oscillations that transmit to other generators (such the i
^{th}
generator
Fig. 6
) by transmission lines, but in this case produced phase lead of CPSS of the i
^{th}
generator is not corresponded with the phase lag of its electrical loop. Consequently, supplied torque for the increase damping will not be the opposite direction of angular velocity variation. Moreover if the frequency difference between the electromechanical mode of the i
^{th}
generator and oscillations that produced by generator B become too much, supplied torque will have the same direction with angular velocity variation. In this situation, CPSS lead to oscillation resonance. This condition is called interfering between PSSs. Hence, when PSSs are designed using of the SMIB system, we must expect this phenomenon. Because the SMIB model includes only a generator and other power systems are modeled with an infinite bus. Below are a few disadvantages and advantages of using the SMIB system in designing of PSSs:
Block diagram of linearized model of the SMIB system along with ideal PSS.
Block diagram of linearized model of the SMIB bus along with CPSS.

• If Power system be too small, an infinite bus cannot be substituted instead of the rest of the power system. Also for modeling the infinite bus instead of the rest of the power system (if be small) needs to have much time and many calculations.
Single line diagram of the multi machine power grid.
As we know, there are different frequency modes in the power system. If any of these is excited, it will produce low frequency oscillations. Hence, when PSSs are designed, it is essential to be considered all the frequency modes. Accordingly, the SMIB system is not a good model for designing PSSs. However, in the multi machine power system, the power system is modeled completely and has not disadvantages of infinite bus. Although using of this model is very complex and is not recommended. As mentioned above, two models that have been described are unsuitable models for designing of PSSs. Therefore, we need a new model for this work. So, a new model which has not disadvantages of the SMIB system and the multi machine power system has been proposed and simulated. This model is a single machine normal bus (SMNB) system.
4. The Single Machine Normal Bus System
To have a precise and complete model of the power system, it is essential that all the power systems are modeled (such as the multi machine power system) or the effects of the rest of the power systems that in the SMIB system are not modeled, via one or several inputs are modeled. Consider the multi machines power system, as shown in
Fig. 3
, the effects of different generators on the i
^{th}
generator are applied via
K
_{1ij}
to K
_{6ij}
(
i
≠
j
) , but as shown in
Fig. 6
all influences and counteractions transmit to the i
^{th}
generator via transmission line ij that has been attached to the power grid. It means that these effects can be considered with adding them to the bus j which is the front of the ith generator. Therefore, the effects of the rest of the power system can be modeled by measuring variables of the transmission ij line such as active currents and reactive currents or active powers and reactive powers or magnitude and angle (
θ_{j}
) of voltage (
V_{j}
) of the bus j. In the SMIB system, the magnitude and the angle of the infinite bus voltage is assumed constant. But from a practical point of view, the magnitude and the angle of the bus voltage which has been connected to the i
^{th}
generator by transmission lines is oscillating for the impacts of power system oscillations. Thus, if in the SMIB system changes of these two variables are considered yields the SMNB system so that any effects of generators are not neglected. Consider the i
^{th}
generator that shown in
Fig. 6
. Axis d of this model is 90 degree lag from its q axis. Relation between fi and its components in q and d axes are as following:
As drawn in
Fig. 7
,
f_{i}
is the indicator of the voltage, the current and the flow of the i
^{th}
generator. The power and the inductive voltage of the i
^{th}
generator and the voltage of the i
^{th}
bus are defined as follows:
Where
Also if Eqs. (1012) are linearized around a certain operating point result:
Therefore, if Eqs. (1012) are entered into Eqs. (14) and then obtained equations are linearized, block diagram of the SMNB is acquired (
Fig. 8
). Comparison between
Figs. 2
and
9
approve that Δ
V
_{j}
and Δ
θ
_{j}
are substituted instead of the rest of the power system which is not modeled in the SMIB system. It means that any kinds of oscillations with different frequencies can be applied with these inputs and then PSSs can be designed with considering the rest of the power system.
Phase diagram of the i^{th} generator
Block diagram of the SMNB system.

• The advantages of this model are:
5. The Calculation of the Constants of the SMNB System
In order to obtain the constants of the SMNB system Eqs. (1012) should be interred in Eqs. (14) and then are linearized around a certain operating point. Hence, linearization of these equations around a certain operating point results in constants
K
_{1}
,
,
K
_{2}
,
,
K
_{3}
,
,
K
_{4}
and
K
'
_{4}
as follows:
To separate the real and the image parts of Eq. (12) for Vi yields:
If Eq. (13) and Eq. (14) are substituted into Eq. (26) and Eq. (27), yields:
where
and
.
Relation between
V_{i}
with
V_{qi}
and
V_{di}
and also Δ
V_{i}
with Δ
V_{qi}
and Δ
V_{di}
is as below:
To linearize Eq. (28) and Eq. (29) around a certain operating point and are substituted into Eq. (31) constants
K
_{5}
,
,
K
_{6}
and
are defined as follows:
6. Problem Formulation
After proposing a new model of the power system for studying low frequency oscillations, it is time that a PSS is designed. Thus a synchronous machine with an IEEE type 2 static AVR system (AVRST2) (
Fig. 9
and
Table 1
) connected to a normal bus through a transmission line has been selected to illustrate the deviation of simplified linear models of the power system for dynamic stability analysis. The excitation system that is used for designing of PSS in linearized model of power system (SMNB), a second separated mode IEEE type 2 static (ST2) excitation system
[13
14]
(
Table 2
and
Fig. 8
) is used. As shown in
Fig. 10
, the model consists of a generator supplying bulk power to a normal bus through a transmission line, with a PSS. The linearized model of the SMNB system that has been linearized around the certain operating point is given as follows:
The parameters of the IEEE type 2 static AVR system (ST2) (Fig. 9)
The parameters of the IEEE type 2 static AVR system (ST2) (Fig. 9)
The IEEE Type 2 static AVRompoundsource rectifier exciter [14]
The Single machine normalbus system
System parameters
Various operating points
The constants K1 to K6 and
to
represent the system parameters at a certain operating condition. Analytical expressions for K1 to K6 as a function of loading (P, Q) are derived in
[15
16]
. Typical data for such a system are as
Table 2
. According to Eqs. (3639), nominal operating point (
Tables 3
and
4
), 100% changes of input power for 18 cycles, ΔVj = 0.05sin (20t) and Δ
θ
j = 0.02sin (20t) a lead PSS as the following form is designed, which stabilizes the system while minimizing Eq. (41).
The objective function that has been used in this paper is demonstrated in Eq. (41). This robust controller can be obtained by solving Eq. (41). Where J
_{1}
= d
ω
, J
_{2}
= dvi and
γ
= t
_{sim}
which t
_{sim}
is simulation time, d
ω
is the angular velocity deviation,
α
and
β
are weighting parameters that in this case
α
is 10 and
β
is 1 and dvi is the terminal voltage deviation of the i
^{th}
generator. Therefore, the design problem can be formulated as the following optimization problem (
Fig. 11
):
PSS controller
Minimize fitness subject to
7. PSO Algorithm
The particle swarm optimization (PSO) algorithm was first proposed by Kennedy and Eberhart
[17]
. Where is a novel evolutionary algorithm paradigm which imitates the movement of birds flocking or fish schooling looking for food. Each particle has a position and a velocity, representing the solution to the optimization problem and the search direction in the search space. The particle adjusts the velocity and position according to the best experiences which are called the pbest found by itself and gbest found by all its neighbors. In PSO algorithms each particle moves with an adaptable velocity within the regions of decision space and retains a memory of the best position it ever encountered. The best position ever attained by each particle of the swarm is communicated to all other particles. The updating equations of the velocity and position are given as follows
[17]
.
Where v is the velocity and x is the position of each particle. c
_{1}
and c
_{2}
are positive constants referred to as acceleration constants and must be c
_{1}
+ c
_{2}
≤ 4 , usually c
_{1}
= c
_{2}
= 2. r
_{1}
and r
_{2}
are random numbers between 0 and 1, w is the inertia weight, p refers to the best position found by the particle and pg refers to the best position found by its neighbors. The flowchart of PSO algorithm is shown in
Fig. 12
.
Flowchart of the PSO algorithm [17]
8. Simulation Results
Typical range of the optimized parameters is [0500] for K, Z and P. The convergence of PSO algorithm is in
Fig. 13
. The optimized parameters are found via PSO algorithm are K = 500, Z = 2.3302 and P = 81.7212. Optimized parameters have been obtained when the input power of the ith generator and the magnitude and angle of the normal bus has been changed 100% for 18 cycles (one cycle in frequency of twenty) instantaneously, ΔVj = 0.04sin (20t) and Δ
θ
j = 0.015sin (20t) respectively and the normal operating point are Pe = 0.9 (pu) and Qe = 0.256 (pu) (
Table 3
and
4
) (the inputs ( ΔVj and Δ
θ
j ) for all tests have been applied for one cycle in their frequency). The dynamic response of the rotor angle deviation, angular velocity deviation and terminal voltage deviation are shown in
Fig. 14
. It is clearly seen that with inclusion of PSS, electromechanical damping characteristics of the system is improved. The parameters of PSS have been obtained using PSO based on the SMNB system. It means that this PSS must have coordination with other PSSs. On the other hand, if any low frequency oscillations reach to the i
^{th}
generator, PSS have to can damp them. So, it is necessary for carrying out some tests on the system to assure that designed PSS has good performance. Hence, three different kinds of tests have been implemented. First test is robustness of designed PSS opposite low frequency oscillations with different frequency. Therefore, different sine oscillations have been applied by ΔVj and Δ
θ
j in second operating point (
Table 4
 tests 12 and 13). The dynamic responses (the rotor angle deviation, angular velocity deviation and terminal voltage deviation) according to these inputs are shown in
Figs. 15
and
16
. As shown these figures, the designed PSS gives very good damping over a range of frequencies. This means that PSS can damp any kinds of low frequency oscillations with different frequencies. Second test is robustness of designed PSS opposite disturbances and harmonics. Thus an input comprises several harmonics with different magnitudes and frequencies have been applied by ΔVj and Δ
θ
j (
Table 4
 test 2). The dynamic responses according to this input are shown in
Fig. 17
. This figure shows that designed PSS has robustness opposite harmonics for damping low frequency oscillations. At last, third test is robustness of designed PSS opposite operating point changes. For this reason, two operating points (Heavy and light) have been used (
Table 4
 test 3). The dynamic responses according to these operating points are shown in
Figs. 18
and
19
. These figures show that designed PSS has good performance opposite changes of operating points and can damp oscillations rapidly. According to the tests, it is clearly observed that designed PSSs based on the SMNB system, have a good performance to damp LFOs and there is not any interfering between PSSs.
Convergence of PSO algorithm
System dynamic response for test 11. (a) rotor angle deviation (pu), (b) angular velocity deviation (pu), (c) terminal voltage deviation (pu).
The values of Constants K1to K6and K′1to K′6according to operating points
The values of Constants K_{1} to K_{6} and K′_{1} to K′_{6} according to operating points
System dynamic response for test 12. (a) rotor angle deviation (pu), (b) angular velocity deviation (pu), (c) terminal voltage deviation (pu).
Applied tests
System dynamic response for test 13. (a) rotor angle deviation (pu), (b) angular velocity deviation (pu), (c) terminal voltage deviation (pu).
9. Conclusion
This paper investigates the stability enhancement problem of the power system in order to maintaining coordination between PSSs via proposing a new model of the power system. A lead PSS has been used for analyzing this model. PSO was employed to search for optimal PSS parameters. This new model is called the single machine normal bus system. Furthermore this model improves performance of the SMIB system and the multi machine power system by regarding the rest of the power system and making simplicity respectively. The effectiveness of. the proposed model and lead PSS have been illustrated under different tests. The results show that designed PSS based on the proposed model achieve good robust performance, provide superior coordination between PSSs in comparison with SMIB model and enhance greatly the dynamic stability of the power systems.
System dynamic response for test 2. (a) rotor angle deviation (pu), (b) angular velocity deviation (pu), (c) terminal voltage deviation (pu).
System dynamic response for test 3 and first operating point (heavy load). (a) rotor angle deviation (pu), (b) angular velocity deviation (pu), (c) terminal voltage deviation (pu).
System dynamic response for test 3 and third operating point (light load). (a) rotor angle deviation (pu), (b) angular velocity deviation (pu), (c) terminal voltage deviation (pu).
BIO
Alireza Fereidouni was born in Behshahr, Iran, in 1987. He received his B.Sc. degree in Electrical Engineering from Babol University of Technology, Babol, Iran, in 2009. He is currently a M.Sc. student in the Department of Electrical Engineering at Amikabir University of Technology, Tehran, Iran. His research interests are power system stability and protection, power electronics, and high voltage.
Behrooz Vahidi was born in Abadan, Iran in 1953. He received the B.S. in electrical engineering from Sharif University of Technology, Tehran, Iran in 1980 and M.S. degree in electrical engineering from Amirkabir University of Technology, Tehran, Iran in 1989. He also received his Ph.D. in electrical engineering from UMIST, Manchester, UK in 1997. From 1980 to 1986 he worked in the field of high voltage in industry as chief engineer. From 1989 to present he has been with the department of electrical engineering of Amirkabir University of Technology where he is now a professor. Prof. Vahidi is selected by the ministry of higher education of Iran and by IAEEE (Iranian Association of Electrical and Electronics Engineers) as the distinguished researcher of Iran. Prof. Vahidi is senior member of IEEE and member of CIGRE WG C4.26. His main fields of research are high voltage, electrical insulation, power system transient, lightning protection and pulse power technology. He has authored and coauthored more than 250 papers and five books on high voltage engineering and power system.
AbdelMagid Y. L.
,
Abido M.A.
2004
“Robust Coordinated Design of Excitation and TCSCBased Stabilizers Using Genetic Algorithms”
Electric Power Systems Research
69
129 
141
DOI : 10.1016/j.epsr.2003.06.009
F. P. Demello
,
C. Concordia
1969
“Concepts of synchronous machine stability as affected by excitation control”
IEEE Trans. Power Apparatus Syst.
88
316 
29
Schleif F. P.
,
White J. H.
1966
“Damping for the NorthwestSouthwest tie line oscillations an analog Study”
IEEE Trans. Power Apparatus Syst.
85
1239 
47
1995
IEEE Working Group on System Oscillations, Inter Area Oscillation in Power Systems, IEEE Special Publications
Larsen E. V.
,
Swann D. A.
1981
“Applying power system stabilizers”, part I and II
IEEE Trans. PAS.
100
3017 
3046
Phiri A.
,
Folly K. A.
2008
“Application of Breeder GA to Power System Controller Design”
the IEEE Swarm Intelligence Symposium
Saint Louis, USA
Sheetekelai S.
,
Folly K.
,
Malik O.
2009
“Design and Implementation of Power System Stabilizers Based on Evolutionary Algorithms”
The IEEE conference
Nairobi
Anderson P. M.
,
Fouad A. A.
1977
“Power System Control and Stability”
Ames, Iowa state University Press
The IEEE conference
Bergen A. R.
1986
“Power system Analysis”
Prentice Hall
Englwood Cliffs
Yu Y. N.
1983
“Electric Power System Dynamic”
Academic Press
New York
Soliman H. M.
,
Saker M. M. F.
1988
“WideRange Power System Pole Placer”
IEE Proc., Pt. C
195 
201
Diao Y.
,
Wasynczuk O.
,
Krause P. C.
1986
“Solution of State Equations in Terms of Separated Modes With Applications to Synchronous Machines”
IEEE Trans. EC. Ec1
IEEE Recommended Practice for Excitation System Models for Power System Stability Studies, IEEE Std. 421.51992
Soliman H. M.
,
Shafei El A. L.
,
Shaltout A. A.
,
Morsi M. F.
2000
“Robust Power System Stabilizer”
IEE Proc. Pt. C
147
285 
291
DOI : 10.1049/ipcom:20000661
Miyuan S.
,
Juan W.
,
Danni P.
2007
“Particle Swarm and Ant Colony Algorithms Hybridized for Multimode Resourceconstrained Project Scheduling Problem with Minimum Time Lag”
IEEE Conference
147
5898 
902
Zhao L.
,
Yang Y.
2009
“PSObased single multiplicative neuron model for time series prediction”
International journal of Expert Systems with Applications
2805 
2812