Optimal Coordinated Design of Multiple Damping Controllers Based on PSS and UPFC Device to Improve Dynamic Stability in the Power System

Unified Power Flow Controller (UPFC) device is applied to control power flow in transmission lines. Supplementary damping controller can be installed on any control channel of the UPFC inputs to implement the task of Power Oscillation Damping (POD) controller. In this paper, we have presented the simultaneous coordinated design of themultiple damping controllers between Power System Stabilizer (PSS) and UPFC-based POD or between different multiple UPFC-based POD controllers without PSS in a singlemachine infinite-bus power system in order to identify the design that provided the most effective damping performance. The parameters of the damping controllers are optimized utilizing a Chaotic Particle Swarm Optimization (CPSO) algorithm based on eigenvalue objective function. The simulation results show that the coordinated design of the multiple damping controllers has high ability in damping oscillations compared to the individual damping controllers. Furthermore, the coordinated design of UPFC-based POD controllers demonstrates the superiority over the coordinated design of PSS and UPFC-based POD controllers for enhancing greatly the stability of the power system.


Introduction
When large power systems are interconnected through weak tie lines, Low Frequency Oscillations (LFO) in the range of 0.1-3 Hz are observed.These oscillations may sustain and grow to cause system separation if no adequate damping is available [1].
The installation of a Power System Stabilizer (PSS) appears as a simple and inexpensive technique for many years to produce an amount of damping torque through the injection of a supplementary stabilizing signal at a voltage reference input of an Automatic Voltage Regulator (AVR), which has increased the stability of the power system.However, their performance deteriorated when the system operating conditions varied widely [2][3][4].Hence, in order to operate power systems effectively without lessening the system safety and quality of supply, even in the case of emergency conditions such as loss of transmission lines and/or generating units, new control strategies need to be applied.Flexible AC Transmission System (FACTS) device is an effective technology that has appeared as an alternative way in controlling power flow and damping power system oscillations.Although the damping ratio of FACTS-based stabilizer controllers is not often their main function, but its ability to increase power system damping characteristics has been acknowledged [5,6].
Unified Power Flow Controller (UPFC) is one of the most important FACTS device families.Its primary function is to control and optimize the real and reactive power flow in a given line, voltage, and current at the UPFC bus [7,8].This is achieved by regulating the controllable parameters of the transmission system such as line impedance, phase angle, and voltage magnitude.Beside these primary functions, it can also provide appropriate damping effect to the tie interconnected modern power systems oscillation through its supplementary controller and increase system stability [9].
To improve general power system performance, possible interactions between PSSs and FACTS-damping controllers are considered, but uncoordinated local design of PSS and FACT-damping controller may cause destabilizing interaction on the damping of system oscillations [10].Therefore, coordinated design of PSSs and FACTS-based stabilizers is a necessity both to make use of the benefits of multiple stabilizers, which will enhance the system stability and prevent the drawbacks accompanied with their operation.
Several researches have been carried out for the coordination between PSSs and FACTS damping controllers.Some of these researches were depending on the complex nonlinear simulation while the others based on the linearized model for power system [11][12][13][14][15][16][17].One approach for achieving the required performance is to make the coordinated design of the controller as a constrained optimization problem.The constraints include limits of controller parameters.The most traditional optimization methods used some deterministic rules to move from one point to another in the decision hyperspace.Therefore, probability of getting stuck at a local optimum is the main problem with these methods [18].
Recently, Particle Swarm Optimization (PSO) technique has appeared as a useful tool for engineering global optimization.PSO is a population-based stochastic optimization method, which employs the swarm intelligence produced by the cooperation and rivalry between the particles in a swarm.PSO is unlike the other evolutionary algorithms having many advantages.The major advantages are as follows: it is fast and simple, does need to apply operators such as GA algorithm, and is easy to be implemented [19,20].The simple PSO greatly depends on its parameters, and it is not guaranteed to be global convergent.In order to improve the global searching ability and prevent a slide into the premature convergence to local minima, PSO and chaotic sequence techniques are combined to form a Chaotic Particle Swarm Optimization (CPSO) technique.
In this paper, we present the results of our comprehensive comparison and assessment of the damping function of multiple damping stabilizers using different coordinated designs in order to identify the design that provided the most effective damping performance.The two alternative designs we evaluated are listed below: (i) coordinated design between PSS and any one out of the four input control channels of the series and shunt structure of UPFC device because any control loop can superimpose a supplementary damping controller to implement the required damping, (ii) coordinated design between any two out of the four input control channels of the series and shunt structure of UPFC device as a multiple damping controller without using PSS.
The parameters of the damping controllers for individual and coordinated design are optimized utilizing CPSO technique based on eigenvalue objective function.The simulation results of the individual damping controllers show the best damping effects resulting from using the POD controllers   and   are higher than PSS response.In addition, the coordinated designs (  &   ), (  &   ), (  &   ), and (PSS &   ) are able to damp the LFO effectively and can accomplish this task faster and with lower overshoot in comparison to their individual control responses or the best individual controls   and   .

Mathematical Model of Power System for SMIB Equipped with UPFC
Figure 1 shows a single-machine infinite bus (SMIB) equipped with a UPFC.The synchronous generator is supplying electric power to the infinite bus through a transmission line and UPFC.The UPFC consists of two three-phase GTO based voltage source converters (VSCs) coupled through a common dc terminal.VSC-1 is connected in shunt with the line through an excitation transformer (ET), and VSC-2 is inserted in series with the transmission line through a boosting transformer (BT).The dc voltage for both converters is provided by a common capacitor bank [21].The structure of the UPFC controller is shown in Figure 2, where  can be   ,   ,   , and   .In order to maintain the power flow between the series and shunt converters, a dc voltage regulator must be incorporated.The dc voltage is controlled through modulating the phase angle of the ET voltage   .Therefore, the   damping controller is to be considered as shown in Figure 3, where the dc voltage regulator is a PI-controller [23].

Nonlinear Dynamic Model of UPFC. Referring to
Figure 1, the nonlinear differential equations describe the dynamic behavior of the UPFC, the three-phase excitation current, three-phase boosting current, and dc link voltage which can be written as follows [24]: where V  ,   ,   , and   are the excitation transformer voltage, current, resistance, and inductance, respectively; V Bt ,   ,   , and   are the boosting voltage, current, resistance, and inductance respectively;  dc , V dc are the dc link capacitance and voltage respectively.By applying Park's transformation and ignoring the resistance and transient of the UPFC transformers, equations (1) become as follows [25]: where   is the ET reactance and   is the BT reactance.From Mathematical Problems in Engineering Figure 1, we can have where  is the armature current,   is the generator terminal voltage,   is the infinite-bus voltage,   and  V are the reactances of a transmission line.The terminal voltage and the armature current of the generator can be expressed in terms of the -axis and -axis components: where x  is the -axis transient reactance,   is the -axis reactance, and    is the transient generator internal voltage.Substituting (6  ) and (7  ) into (6) gives Equation ( 4) can be expressed in the - axis reference frame as: From ( 2), (8), and ( 9), we have obtained the current injection equations of UPFC: where where   is the -axis reactance.

Nonlinear Model of the Power System.
The nonlinear dynamic equations of the SMIB system shown in Figure 1 are described by [26]: where   is the synchronous speed,  and  are the angle and speed of the rotor, respectively,   and   are the input mechanical and output electrical powers of the generator, respectively,  and  are the machine inertia constant and damping coefficient, respectively,   is the generator field voltage, and    is the open-circuit field time constant.The excitation system is represented by a first-order model (IEEE type-ST1) [27]: where   and   are the gain and time constant of excitation system, respectively; V ref is the reference voltage.
The output power of the generator can be expressed in terms of the -axis and -axis components of the armature current and terminal voltage as

Linearized Model of a Power System Equipped with UPFC.
Linear dynamic model of the power system is obtained by linearizing the nonlinear equations ( 2)-( 14) around nominal operating point.The linearized model of the power system as shown in Figure 1 is given as follows: where where the linearization constants  1 - Referring to Figure 2, the UPFC dynamic models of Δ  , Δ  , and Δ  are described by the following state equations: In addition, from Figure 3, the UPFC dynamic model of a dc voltage regulator is described by the following state equations: where   and   are the proportional and integral gains of the dc voltage regulator controller, respectively; ( 1 - 4 ) and ( 1 - 4 ) are the gains and times constants of the main control loops for UPFC, respectively.Equation (15) describes the model of the machine.Equations ( 16) and ( 18) represent the models for the exciter and the UPFC dc voltage regulators link, respectively.Equations (19) and (20) represent the control action of the main control loops of the UPFC with damping controller, and (21) represents the control action of PI dc voltage regulator.In state-space representation, these equations can be arranged in compact form as where the state vector Δ and control vector Δ are The structure of the matrices  and  is Mathematical Problems in Engineering The modified Phillips-Heffron transfer function model including UPFC has 28 constants; on the other hand, the Phillips-Heffron model has only 6 constants as shown in Figure 4.This figure shows that the damping torque provided by the UPFC device consists of two parts.The first part is applied directly to the electromechanical oscillation loop of the generator.Its sensitivity is mainly measured by coefficients   ,   ,   ,   , and   , which is called the direct-damping torque.The second part consists of coefficients  V ,  V ,  V ,  V , and  V , which is applied through the field channel of the generator.Its sensitivity is related to the deviation of field voltage, which is referred to as the indirect-damping torque.

Eigenvalue Analysis without Stabilizer.
For nominal operating condition, the dynamic behavior of the system is recognized through the eigenvalues of the system matrix .By solving the system characteristic equation | − | = 0, the eigenvalues of the system are computed which are given below: It is clearly seen from eigenvalues of the matrix  that the system is unstable and needs a supplementary stabilizer for stability.

Structure of the Supplementary Damping Stabilizers PSS and POD Controller
In order to overcome the LFO problem, supplemental control action is applied to the generator excitation in the form of PSS or UPFC device as POD controller.The four main control parameters of the UPFC (  ,   ,   , and   ) can be modulated in order to produce the appropriate damping torque, which is illustrated in Figure 2.
The POD controller has a structure similar to that of the PSS. Figure 5 shows a sample block diagram of a POD controller.The controller contains three main blocks, that is, the gain block, the washout filter block, and two blocks (leadlag) phase compensators.It provides an electrical torque in phase with the speed deviation in order to improve damping of power system oscillations.The washout filter block acts as a high-pass filter to eliminate the dc offset of the POD output and prevent steady-state changes in the terminal voltage of the generator.From this perspective, the washout time   should have a value in the range of 1 to 20 seconds defined to the electromechanical oscillation modes [1].In this study, the time constants   ,  2 , and  4 were assigned specific values of 10 s, 0.1 s, and 0.1 s, respectively, while the parameters of the POD controller such as   ,  1 , and  3 need to be determined.The speed deviation Δ used as an input signal to the POD and Δ  is the controller output, where  = PSS,   ,   ,   , and   .With usage of the above dynamic damping controller, the number of matrix state variables increases from 10 to 13, by adding three state variables Δ Ẋ1 , Δ Ẋ2 , and Δ U .

Optimal Design of the Supplementary Damping Stabilizers (PSS or POD).
The main objective of optimization technique is to improve the dynamic stability of the power system against disturbances at different loading conditions.It can be achieved by suitable tuning of damping controller parameters.
The supplementary damping stabilizer (lead-lag type) can be described mathematically as where () is the transfer function of the supplementary damping stabilizer, () is the measurement signal, and () is the output signal of the supplementary damping stabilizer, which will provide additional damping by moving electromechanical mode to the left.Equation ( 26) can be expressed in state-space form as where Δ  is the controller state vector.Equation ( 22) describes a linear model of the power system extracted around a certain operating point.combining (22) with ( 27), we obtained a closed-loop system: where Δ ℓ is the state vector of the closed loop system,   is the th eigenvalue mode of the closed loop matrix  ℓ and   is the damping coefficient of the th eigenvalue.It is clear that the objective function  will identify the minimum value of the damping coefficient among all modes.The goal of optimization process is to maximize  in order to realize appropriate damping for all modes including electromechanical mode, by exploring the search space of admissible control parameters, which enhance the system damping characteristics, and maximum  is searched within the limited range of control parameters: Typical ranges of the optimized parameters are 0.01-100 for   and 0.001-1 for   .

Optimization Process.
The problem of tuning the parameters for individual and coordinated design for multiple damping controllers, which would ensure maximum damping performance, was solved via a PSO optimization procedure that appeared to be a promising evolutionary technique for handling optimization problems.PSO is a populationbased, stochastic-optimization technique that was inspired by the social behavior of flocks of birds and schools of fish [28].
The advantages of PSO algorithm are that it is simple and easy to implement and it has a flexible and well-balanced mechanism to enhance the local and global exploration capabilities.Recently, it has acquired wide range of applications in solving optimization design problems featuring nonlinearity, nondifferentiability, and high dimensionality in many area search spaces [29,30].

Classical PSO Algorithm.
In the PSO, each possible solution is represented as a particle, and each set of particles comprises a population.Each particle keeps its position in hyperspace, which is related to the fittest solution it ever experiences in a special memory called best.In addition, the position related to the best value obtained so far by any particle in the population is called best.For each iteration of the PSO algorithm, the best and best values are updated   and each particle changes its velocity toward them randomly.This concept can be expressed as [18] V where V is the particle velocity,  is the particle position,  is the number of iterations,  is the inertia weight factor,  1 and  2 are the cognitive and asocial acceleration factors respectively,  is the number of particles, and  1 and  2 are the uniformly distributed random numbers in the range of 0 to 1. Figure 6 shows the flow chart of the PSO algorithm.

Gain Washout
Phase compensation 1 and 2

Chaotic Particle Swarm Optimization (CPSO).
The main disadvantage of the simple PSO algorithm is that the performance of it greatly depends on its parameters and it is not guaranteed to be global convergent.In order to improve the global searching ability and prevent a slide into  the premature convergence to local minima, PSO and chaotic sequence techniques are combined to form a chaotic particle swarm optimization (CPSO) technique, which practically combines the population-based evolutionary searching ability of PSO and chaotic searching behavior.The logistic equation employed for constructing hybrid PSO is described as [31]: where  is the control parameter with a real value between 0 to 4. Although (32) is deterministic, it exhibits chaotic dynamics when  = 4 and  0 ∉ {0, 0.25, 0.5, 0.75, 1}.It exhibits the sensitive dependence on initial conditions, which is the basic characteristic of chaos.The inertia weighting function in (30) is usually evaluated utilizing the following equation: where  max and  min are maximum and minimum values of , iter max is the maximum number of iterations, and iter is the current iteration number.The new weight parameter  new is defined by multiplying weight parameter  in (33) and logistic equation (32): To improve the global searching capability of PSO, we have to introduce a new velocity update equation as follows: We have observed that the proposed new weight decreases and oscillates simultaneously for total iteration, whereas the conventional weight decreases monotonously from  max to  min .The final choice of a parameter was considered to be the optimal choice: , iter max ,  1 ,  2 ,  min ,  max , , and  0 are chosen as 30, 100, 2, 2, 0.3, 0.9, 4, and 0.3, respectively.

Simulation and Comparison Results
In this section, the ability of UPFC in damping system oscillation and the dynamic interactions of UPFC-POD controllers are investigated intensively.The CPSO technique has been applied to design individual and coordinated damping controllers.To evaluate the performance of the proposed simultaneous coordinated design approaches, the responses with the proposed controllers were compared with the responses of the individual design of the PSS and UPFC-POD controllers for two schemes.The resultant optimal parameters of the individual controllers and coordinated designs are given in Tables 1, 2, and 3, respectively.4 and 5, respectively.The first and second rows represent the electromechanical mode and their damping ratio  EM using participation factor to identify the eigenvalue associated with electromechanical mode.It is clear that the proposed controllers stabilize the system effectively.The maximum overshoot as well as the settling time of the system responses has been compared to the  capability of the simulated coordinate approaches in oscillation damping.Figure 7 shows the speed deviation responses of the four UPFC-POD controllers and PSS individually.All these responses are shown in one figure for comparison purposes.It can be seen that by using the UPFC-POD controllers   and   provide a robust damping characteristic higher than that of PSS and enhance the stability.However, the UPFC-POD controllers   and   do not perform well and have poor capabilities in damping the system oscillation.Figures 8 and 9 show the improvement in overshoot and settling time of the system response while using the coordinated designs (PSS &   ) and (PSS &   ) over the individual control responses but less than the best individual controls   and   .Figure 10 shows that the coordinated design (PSS &   ) improves the system damping compared to their individual controls but no more than the best one   .Figure 11 shows the best response in Scheme 1 while using the coordinated design (PSS &   ) over their individual control  responses because it improved both damping parameters overshoot and settling time.

Scheme 2.
In this scheme, the power system is considered to possess a UPFC device without PSS.Figures 12-17 show the effect of applying the individual controllers and coordinated designs between different UPFC-POD controllers in a SMIB of speed deviation of Δ with 10% step change in mechanical input power.The eigenvalues of the system with the proposed controllers in scheme 2 are given in Table 6.It is clear that the proposed controllers stabilize the system effectively.
From Figures 12,13, and 14, it is clear that the system damping resulting from using coordinated designs (  &   ), (  &   ), and (  &   ) are much higher than those resulting from using the best individual controls   and   .Figures 15 and 16 show that the coordinated designs (  &   ) and (  &   ) improve the system damping compared to their individual controls but not more than the best one   .Figure 17 shows the improvement in overshoot  and settling time of the system response while using the coordinated design (  &   ) over the individual control responses but less than the best individual controls   and   .

Conclusions
This paper is concerned with the damping of LFO via PSS and UPFC-POD controllers applied independently and also through the simultaneous coordinated designs in a SMIB power system.To improve the global searching ability and prevent a slide into the premature convergence to local minima, PSO and chaos theory are combined to form a CPSO.For the proposed controller design problem, a CPSO algorithm was used as the optimization technique to search for the optimal damping controller parameters in both the individual and the coordinated designs.The simulation results of the individual damping controllers showed the best damping effects resulting from using   and   as POD controllers are much higher than those resulting from using PSS,   , and   responses.In addition, the coordinated designs (PSS &   ) in scheme

Figure 3 :
Figure 3: UFPC with dc voltage regulator and damping controller.

Figure 4 :
Figure 4: Modified Phillips-Heffron transfer function model of SMIB system with UPFC.

Figure 5 :
Figure 5: Structure of the supplementary damping stabilizers (PSS or POD).

Figure 6 :
Figure 6: PSO algorithm for the tuning parameters of an individual and coordinated design.

4. 1 .
Scheme 1. show the effect of applying the individual controllers and coordinated designs for PSS and different UPFC-POD controllers in a SMIB of speed deviation of Δ with 10% step change in mechanical input power.The system eigenvalues with the proposed individual stabilizers and coordinated designs for nominal operating condition are given in Tables

Figure 8 :
Figure 8: Speed variation responses for individual damping controllers PSS,   and coordinated design PSS &   .

Figure 9 :Figure 10 :
Figure 9: Speed variation responses for individual damping controllers PSS,   and coordinated design PSS &   .

Figure 11 :
Figure 11: Speed variation responses for individual damping controllers PSS,   and coordinated design PSS &   .

Figure 12 :
Figure 12: Speed variation responses for individual damping controllers   ,   and coordinated design   &  E .

Figure 13 :Figure 14 :
Figure 13: Speed variation responses for individual damping controllers   ,   and coordinated design   &   .

Figure 15 :
Figure 15: Speed variation responses for individual damping controllers   ,   and coordinated design   &   .

Figure 16 :
Figure 16: Speed variation responses for individual damping controllers   ,   and coordinated design   &   .

Table 1 :
The optimal parameters of the individual controllers.

Table 2 :
The optimal parameters of the coordinated designs between PSS and UPFC-POD controllers.

Table 3 :
The optimal parameters of the coordinated designs between different UPFC-POD controllers.

Table 4 :
System eigenvalues of the individual controllers.
1; (  &   ), (  &   ) and (  &   ) in scheme 2 provide superior performance in comparison with the best individual controls   and   while the coordinated designs (PSS &   ) in scheme 1; (  &   ), and (  &   ) in scheme 2 improve the system damping compared to the individual controls except   , which is