Pss Tlbo PDF
Pss Tlbo PDF
Pss Tlbo PDF
Abstract—This paper deals with an interesting application of Transmission devices (FACTS), which can enhance the power
recently evolved Teaching Learning based Optimization (TLBO) system stability and power transfer capability. They are
algorithm in designing coordinated Proportional-Integral (PI) economical, fast acting and can improve the efficiency and
controller based Power System Stabilizer (PSS) for single security of power system [4]. Thyristor Controlled Series
machine infinite bus power system equipped with Thyristor Compensator is one of the first generation FACTS devices. It
Controlled Series Compensator (TCSC). As the design is for is economical and effective means of enhancing dynamic
coordinated system, traditional TLBO results in suboptimal stability of power system by quick and flexible means of
solution and hence we propose a modified TLBO method based adjusting line reactance. It assures better control over power
on the concept of opposition based learning for designing
flow, improvement of transient stability limits and fault
coordinated controllers. Computer simulations of the proposed
approach on various loading conditions reveal the superiority of
current limitation [5-9]. For the small signal stability studies
modified TLBO in designing coordinated controller for of Single Machine Infinite Bus system (SMIB) linear model of
enhancing the dynamic stability of power system. Philip-Heffron is considered. To avoid the destabilizing
interactions the tuning of TCSC controller is coordinated with
Keywords: PSS, TCSC, TLOBA, small signal stability, SMIB. PSS. To further enhance the dynamic stability, PI controllers
are incorporated along with TCSC and PSS controllers.
I. INTRODUCTION As the coordinated controllers consist of more parameters
to be selected judiciously for better performance of the power
Owing to the growing complexity of modern day power
system, this calls for real parameter optimization in n-
systems, they are often interconnected with weak tie lines.
dimensional hyperspace. To carry out this optimization task
Fast acting, high gain Automatic Voltage Regulator’s (AVR)
we chose Teaching Learning Based Optimization (TLBO)
are being employed to the synchronous generators to maintain
algorithm, a newly evolved optimization algorithm. TLBO
the distantly located, inter connected power systems at
draws its inspiration from knowledge sharing phenomenon
constant operating voltage [1]. Though AVRs can enhance the
between students and teacher in a classroom. To further
overall transient stability, they are responsible for low
enhance the performance of TLBO method in designing the
frequency generator rotor angle oscillations (0.1-3 Hz). They
coordinated controllers we propose a new variant TLBO
may further grow in magnitude affecting the small signal
method based on the concept of Opposition. It is referred as
stability, which is the ability of the power system to remain in
TLOBA i.e., teaching learning opposition-based algorithm.
synchronism when the small disturbances due to variations in
The rest of paper is organized as follows. Section II deals
generation and loads occur [2].
with mathematical modeling of power system considered. In
In order to produce positive damping on these small
Section III a brief outline of problem is discussed. Section IV
frequency oscillations, Power System Stabilizers (PSS) are
summarizes the proposed approach followed by design
employed. The purpose of PSS is to introduce supplementary
perspectives in Section V. In section VI we elucidated the
signals (derived from speed deviation signal Δω) in the performance of modified TLBO over various loading
feedback loop of voltage regulator. Design of effective PSS is conditions and at the end we provide some conclusions and
very difficult when the frequency of oscillations begun to vary future scope in Section VII.
over a wide range. Also PSS causes variations in voltage
profiles and their operation is relatively slow [3].
The recent advancements in the high power semiconductor
technology lead to the development of Flexible AC
II. POWER SYSTEM MODELING B. PSS and Excitation system
The Single Machine connected to Infinite Bus through The conventional two-stage lead-lag Power System
transmission line with TCSC controller shown in Figure 1 is Stabilizer is considered in this study. IEEE Type-ST1A
being considered for small signal stability studies Excitation system is considered. The inputs to excitation
system are terminal voltage (VT), supplementary signal (Vs)
A. Generator Modeling from PSS and reference voltage (Vref). KA and TA are the gain
The 3rd order model consisting of the swing equation and and time constant of excitation system respectively.
the generator internal voltage equation describes the
generator. IEEE–ST1type excitation system is considered.
Figure 1: Single machine infinite bus power system with TCSC [6] Figure2: PSS and IEEE Type-ST1A Excitation system [6]
angle respectively.
ΔTD
=
TDω0Δω
≅
KPKTKDΔω (15) In the above equations, Δω(t, X) denotes the rotor speed
The
transfer
functions
of
the
PSS
and
the
TCSC
deviation for a set of controller parameters X. Here X
controller
are
(8)
and
(9)
respectively:
represents the parameters to be optimized. The optimization
is carried in two phases, initially the 10 parameters
⎛ sTWP ⎞⎛ 1 + sT1P ⎞⎛ 1 + sT3 P ⎞ corresponding to both TCSC and PSS controller are been
u PSS = K P ⎜⎜ ⎟⎟⎜⎜ ⎟⎟⎜⎜ ⎟⎟
(16)
tuned coordinately and in second phase by fixing the obtained
⎝ 1 + sTWP ⎠⎝ 1 + sT2 P ⎠⎝ 1 + sT4 P ⎠ parameters of TCSC and PSS controllers, the PI parameters
⎛ sTWT ⎞⎛ 1 + sT1T ⎞⎛ 1 + sT3T ⎞ Kp and Ki of both TCSC and PSS are tuned coordinately to
uTCSC = KT ⎜⎜ ⎟⎟⎜⎜ ⎟⎟⎜⎜ ⎟⎟ (17)
obtain optimum system response.
⎝ 1 + sTWT ⎠⎝ 1 + sT2T ⎠⎝ 1 + sT4T ⎠
In this structure, the washout time constants TWT and TWP
are usually pre-specified, TWT= TWP=5s. The controller gains
KT & KP and the time constants T1T , T2T, T3T, T4T, T1P , T2P,
T3P , T4P are to be determined.
Δω ⎛ sTw ⎞ ⎛ 1 + sT1 ⎞
Input K ⎜⎜ ⎟⎟ ⎜ ⎟
⎝ 1 + sTw ⎠ ⎝ 1 + sT 2 ⎠
u Ki ⎛ 1 + sT3 ⎞
Kp + ⎜ ⎟
Output s ⎝ 1 + sT 4 ⎠
two phases i.e., (i) Teacher Phase and (ii) Learner Phase X new
j,k,i
= X new
j,k,i (
+ rand i X new
j,k,i
− X new
j,Q,i )
i. Teacher Phase ELSE
In this his phase a teacher tries to ameliorate the mean X new
j,k,i
= X new
j,k,i (
+ rand i X new
j,Q,i
− X new
j,k,i )
result of class in the subject taught by him/her based on level End IF
of knowledge and skill he/she had in that particular subject. End FOR
For any ith iteration, let us consider there are m number of
Accept X new if it gives a better function value.
subjects (design variables), n number of learners (population j,k,i
size, k=1, 2,..,n) and Tj,i be the mean result of the learners in
B. Teaching Learning Opposition Based Optimization
jth subject (where j=1,2,…m). However the best overall result
i. Opposition-Based Learning
Xtotal-k-best,i (considering all the subjects together) in a class of
Most of the evolutionary optimization methods start with
learners can be considered as result of best leaner k-best and
some initial solutions and usually start with random guesses.
the best learner identified is replaced by the teacher. As the
The computational time depends upon the distance between
teacher Xtotal-k-best,i will try to move mean Ti towards its own
initial guess and optimal solution. Hence if the guess is not in
level, an adaptive heuristic is used to update the solution and
the vicinity of optimal solution computation time may
is done according to the difference between the existing mean
increase. The chance of improving our convergence can be
result of each subject and the corresponding result of the
done by starting with a fitter solution by simultaneously
teacher for each subject is given by.
checking the opposite solution [11]. If x is a obtained solution
Difference _ mean j,k,i = rand i ( X j,k−best,i − TFTj,i ) (21) of given function then the opposite solution x’ can be
where TF is termed as teaching factor, which decides calculated as follows
whether the value of mean is to be changed or not. The value x' = a + b − x (23)
of TF can be either 1 or 2, which is decided randomly with where x ∈ R within an interval of [a,b]
equal probability and randi is a random number in the range
[0, 1]. Xj,k-best,i is the result of the teacher in subject j. ii. Opposition-based Optimization
Difference_Meanj,k,i, defined in Eqn (21) is used in updating Let P = {x1 ,x 2 ,...,x D } be a point in D-dimensional space,
the existing solution according to the following expression.
where x1 ,x 2 ,...,x D ∈ R and x i ∈ [a i ,bi ] ∀i ∈ {1,2,...,D} Now
X new
j,k,l
= X j,k,l + Difference _ Mean j,k,i (22)
the opposite point P' = {x1' ,x '2 ,...,x 'D } is defined as
where X new
j,k,i
and X j,k,i are the new and existing values
corresponding to jth subject of kth learner of ith iteration. A x 'i = a i + bi − x i (24)
greedy mechanism is performed between X new and X j,k,i , the Now, with above definition of opposite point the
j,k,i
opposition based optimization can be formulated as follows.
learner with better function value is retained.
ii. Learner Phase
()
Assuming f ⋅ is fitness function via which candidate fitness
In the course of time a learner may interact randomly with is measured and according to the above given definitions of
other learners with the help of communications, discussions, P and P' , if f (P' ) ≥ f (P) then the point P can be replaced
etc. If a leaner interacts with other learner who has more
with P' ; hence, the point and its opposite point are evaluated
knowledge than him or her, he/she tries to learn new things simultaneously in order to go with the fitter one.
and tries to increase his/her knowledge. For a class of n
learners the learning phenomenon of this phase is expressed iii. Proposed Algorithm
with following pseudo code. Opposition scheme discussed above is applied two times
for the proposed TLOBO method at starting of teaching phase
and learning phase respectively. Once the algorithm has
started with random initial population simultaneously
opposite population are also calculated and then best n values
are picked up (based on the fitness value) and then passed in
to the teacher phase. Similarly before entering in to the
learning phase opposite population is evaluated and the best n
values are passed in to the learning phase and the rest is same
as that of TLBO. This is continued till the termination
0.16
criterion is reached. PSO - PSS
Instead of using predefined interval boundaries 0.14 PSO tuned TCSC-PSS
TLOBA tuned TCSC-PSS
[a i ,bi ] here we used the minimum and maximum values TLOBA tuned PI basedTCSC-PSS
0.12
( [a min ,b max ] ) of each dimension in current population to
as: 0.06
OPi, j = a min
j
+ b max
j
− Pi, j (25)
0.04
where Pi,j is the jth vector of the ith learner in the
0.02
population. OPi,j is the opposite position of Pi,j ; a min
j
and
0
b max
j
are the minimum and maximum values of the jth 0 1 2 3 4
time (sec)
5 6 7
T3T, T4T, T1P, T2P, T3P, T4P and PI controller parameters (Kp, Ki)
of both TCSC and PSS controllers. (the parameters with -2
0.1
responses obtained for nominal, heavy and light loaded Figure 8(a): Rotor Angle Deviation: Light loaded
systems are depicted in terms of speed deviation and rotor
angle deviations. Figure 8 shows the convergence
characteristics of TLOBA progressing towards optimum
values without and with PI controllers. Table 3 shows the
time domain indices values for different loading conditions in
terms of peak value and settling time.
x 10
-4
shoot and low settling time of speed deviation response.
10
PSO - PSS However the peak over shoot value for the rotor angle
PSO tuned TCSC-PSS deviation response with proposed controller is a bit high
TLOBA tuned TCSC-PSS when compared to case without coordinated tuned PI
TLOBA tuned PI based TCSC-PSS
controllers for TCSC and PSS controllers.
5
Similarly Figures 8(a) and 8 (b) shows the speed deviation
speed deviation
0.08
2.4
0.06
2.2
0.04
2
0.02
1.8
0
0 2 4 6 8 10
time (sec)
1.6
Figure 9(a): Rotor Angle Deviation: Nominal loaded 0 50 100 150 200
Number of Function Evaluations
x 10
-4
Figures 10(a): Convergence of TLOBA towards minimum: without
10
PI Controller
PSO-PSS -4
x 10
8 PSO tuned TCSC-PSS 1.6
TLOBA tuned TCSC-PSS NOMINAL SYSTEM
6 TLOBA tuned PI based TCSC-PSS LIGHT LOADED SYSTEM
1.4 HEAVY LOADED SYSTEM
Objective Function (In Log scale)
4
speed deviation
1.2
2
0 1
-2
0.8
-4
0.6
-6
0 1 2 3 4 5 6 7
time (sec)
0.4
Figure 9(b): Speed Deviation: Nominal loaded 0 50 100 150 200
Number of Function Evaluations
Figures 7(a) and 7(b) show the rotor angle deviation and Figure 10(b): Convergence of TLOBA towards minimum: without
the speed deviation responses respectively for the heavy PI Controller
loaded system. From Figures 7(a-b) and the time domain
indices recorded in Table 3 it is clear that proposed PI
controller based TCSC-PSS has produced less peak over
Table 1: Parametric Values Obtained for coordinated TCSC- PSS Using TLOBA and obj func. Minimization values
Table 2: Parametric Values Obtainedof coordinated PI controller TCSC- PSS Using TLOBA and obj func. Minimization values
Parameter Nominal System Heavy loaded System Light Loaded System
Table 3. Settling times and peak values for various loading conditions