Design of Robust Controller for LFC of Interconnected Power System Considering Communication Delays

The usage of open communication infrastructure for transmitting the control signals in the Load Frequency Control (LFC) scheme of power system introduces time delays. These time delays may degrade the dynamic performance of the power system. This paper proposes a robust method to design a controller for multi-area LFC schemes considering communication delays. In existing literature, the controller values of LFC are designed using time domain approach which is less accurate than the proposed method. In proposed method, the controller values are determined by moving the rightmost eigenvalues of the system to the left half plane in a quasi-continuous way for a preset upper bound of time delay. Then the robustness of the proposed controller is assessed by estimating the maximum tolerable value of time delay for maintaining system stability. Simulation studies are carried out for multi-area LFC scheme equipped with the proposed controller using Matlab/simulink. From the results, it has been concluded that the proposed controller guarantees the tolerance for all time delays smaller than the preset upper bound and provides a bigger delay margin than the existing controllers.


Introduction
For many years, the Load Frequency Control (LFC) plays a major role in power system operation and control.The main objective of LFC is to minimize the frequency variations when there is any change in load [1].In tra-ditional LFC framework, the control signals are transmitted through dedicated communication channels where the communication delays are very minimal.Nowadays, the open communication infrastructure is widely used for transmitting control signals.In such a case, time delays cannot be neglected.These time delays are certain to have an impact on the stability of the power system [2] [3].Therefore it is essential to consider the communication delays while analyzing the stability of power system.
At present, there is a rapid momentum in the advancement of research to deal LFC with communication delays.X. Yu et al. [4] proposed a Linear Matrix Inequalities (LMI) approach for LFC system with communication delays.Hassa Bevarani et al. [5] designed a robust decentralized PI controller based on H 2 /H ∞ control technique for three area interconnected power system with communication delay.L. Jiang et al. [6] examined the delay dependent stability of multi area LFC scheme with PI controllers using Lyapunov-theory based delay dependent criterion and LMI techniques.R. Dey et al. [7] investigated the delay dependent/independent design of H ∞ controller for LFC of two area system.Chuanke Zhang et al. [8] analyzed the delay dependent stability of LFC of multi area system by finding the delay margin using LMI technique and obtaining the relationship between delay margin and controller parameter.Chuanke Zhang et al. [9] designed a PID controller for delay dependent robust load frequency control.Sahin Sonmez et al. [10] discussed about the computation of time delay margin for single area LFC system using Routh array stability criterion.
J. Chen et al. [11] presented a method to compute the delay margin of a linear time delay system by determining eigenvalues and generalized eigenvalues of certain constant matrices.The method of tracing critical eigenvalue [12] and cluster treatment of characteristics roots [13]- [15] is a direct method of finding delay margin.Another method is the indirect method to determine delay margin based on Lyapunov stability theory and LMI [16]- [18].Wim Michiels et al. [19] suggested a new method for the determination of controller parameters in a broad class of linear control systems affected by time-delays.In this outlook, the research is focusing on designing a robust controller for LFC affected by communication delays gain prominence.
This paper proposes a new method to design a robust controller for the multi area LFC scheme considering communication delays in order to maintain the frequency and tie-line power between the areas.The controller is designed to guarantee the stability of power system for any delays smaller than the preset upper bound.The paper is organized as follows.In section 2, the multi area LFC structure is modelled considering time delays and it is represented in state space form.In section 3, the detailed description of a robust method to design the controller for multi area LFC affected by communication delays is presented.In section 4, the efficiency of the proposed controller is evaluated by computing the maximum tolerable value of time delay theoretically using Frequency Sweeping Test.In section 5, simulation is performed to prove the efficiency of the designed controller against delays for multi area LFC.

Dynamic Model of Multi Area LFC with Time Delay
This section illustrates the dynamic model of multi area LFC scheme with time delay.This is obtained by including an exponential term e −sτ in the secondary control loop of the conventional LFC model [1].Figure 1 shows the block diagram of i th control area of multi-area LFC scheme where i = 1, 2…N.The exponential term denotes the time delay.The turbine, governor and generator are modelled by a first order transfer function [1].
The notations used for i th control area are listed below.The two-area LFC scheme with time delay can be expressed in state space form as where x(t) is state vector and the state variables are The ACE signal is used as the input to load frequency controller, which is designed as The closed loop model of two area LFC system can be obtained by modifying Equation (1) using state output feedback method and expressed as where

Governor
Turbine Generator the gain values of PI controller.The controller values are determined using Continuous Pole Placement method which is described in detail in the next section.

Description of the Algorithm
The idea behind the proposed stabilization method is to move the unstable eigenvalues to the left half plane in a quasi-continuous way by applying small changes to the controller gain, in the meanwhile monitoring the other eigenvalues with a large real part.The proposed stabilization method is based on the Theorem given in Appendix 1.
The algorithm for the proposed method is as follows: Step 1. Initialize the number of rightmost eigenvalues m = 1.
Step 2. Compute the rightmost eigenvalues for a particular preset upper bound of delay.
Step 3. Find the sensitivity of the m rightmost eigenvalues with respect to the changes in the controller gain K.
Step 4. Using the sensitivities computed in step 3, shift the m rightmost eigenvalues in the direction towards the left half plane by applying small changes to the controller gain K.
Step 5. Meanwhile monitor the uncontrolled eigenvalues.Stop when the stability is reached or go to step 2. The detailed explanation about the different steps involved in the algorithm is presented in the following sections.

Computation of the Rightmost Eigenvalues
In 1999, Engelborghs and Roose proposed a method which computes the rightmost eigenvalues of the characteristic equation.In this method, a discretization of the time integration operator of the linearized system is obtained.The eigenvalues of the linearized system are exponential transforms of the roots of the characteristic equation.Then, the selected eigenvalues of the resulting large matrix are computed.A step length heuristic is applied to ensure that all eigenvalues are approximated accurately by discretization and the accuracy is improved by employing Newton iteration on the characteristic equation taking the approximate eigenvalues as starting values.This method is implemented in the Matlab package DDE-BIFTOOL.This package is a collection of matlab routines used to find the right most eigenvalues of the system.

Sensitivity of Eigenvalues with Respect to the Controller Gain K
The characteristics equation of the two area LFC system Equation ( 5) can be written as ( ) ( ) where λ i is a solution of the characteristic equation and n(v i ) is a normalizing condition.Differentiating the Equation (7) and Equation ( 8) w.r.t. a component k j of K, ( ) From Equation ( 9),

Continuation of Eigenvalues as a Function of the Feedback Gain K
It is assumed that m eigenvalues ( ) , , , m λ λ λ  are to be controlled.The sensitivity of the eigenvalues with respect to changes in feedback gain is denoted as the sensitivity matrix , , , From Equation (10) the small changes in gain S is the Moore-Penrose inverse of sensitivity matrix (S m ).Thus with the availability of d m Λ and † m S , the small adjustments in feedback gain can be obtained one constraint that is imposed on the feedback gain is that its components must be real, but this can be obtained only by taking the components of d m Λ in complex conjugate pairs.For the new feedback gain m m K K + ∆ , the displacement of the controlled eigenvalues is not equal to d m Λ , because Equation ( 10) is based on linearization, and some correction is required.Since the eigenvalues and eigen functions are continuous with respect to parameter changes, the predictor for few Newton iterations on Equation ( 7) are required when m K ∆ is small.Since the research focuses on designing a controller in stability point of view, it is sufficient to control only the real parts of the eigenvalues i.e., . The real parts of the rightmost eigenvalues are adjusted in order to shift the rightmost eigenvalues to LHP.Therefore Equation ( 11) is modified as, where ( ) ∆ℜ Λ is the desired displacement of the real parts of the controlled eigenvalues.To illustrate the movement of the real parts of the rightmost eigenvalues of two area LFC scheme for the preset upper bound delay 4s, Figure 2 is shown.
Here the rightmost unstable eigenvalues whose real part is positive are shifted to negative real axis.The proposed algorithm converges to an optimum value at iteration on 142.At iteration 142 all the rightmost eigenvalues are moved to LHP.The final value of controller gain is K 1 = [−0.0233,−0.0146] T , K 2 = [0.0127,−0.0335] T .The spectrum of eigenvalues for final value of the controller gain is depicted in Figure 3.
The proposed controller guarantees the stability for the delays smaller than the preset upper bound.The robustness of the controller is validated by finding an index called delay margin.The delay margin of the system is defined as the maximum tolerable value of time delay after which the system goes unstable.The delay margin of the system is computed using Frequency Sweeping Test which is discussed in the next section.

Frequency Sweeping Test (FST)
Frequency sweeping Test [2] is applied to analyze the asymptotic stability of the system and determine the value of delay margin.The stability of the system exists only for a subset of nonnegative delays.The necessary and sufficient condition for delay-dependent stability of the system is based on the following theorem.
Let d τ be the delay margin.It is assumed that the system is stable at τ = 0.The rank of matrix A d is k.The delay margin can be defined as That is, the system remains stable for all [ ] 0, d τ τ and the system becomes unstable at τ = τ d.The explana- tion for the above theorem is presented below  The generalized eigenvalues of the matrix pencil ( )  The generalized eigenvalues of the matrix pencil ( ) = .There exists a pair ( i q ω , i q θ ) at which the absolute value of each eigenvalue variation reaches one such that delay margin

Algorithm for Delay Margin Calculation
The algorithm of Frequency sweeping Test for computing the delay margin is given below: 1) Obtain the maximum of real parts of all eigenvalues of the matrix A + A d .If it is less than zero proceed.
2) Find the rank of matrix A d .Consider rank(A d ) = k.then the number of crossover points from right half plane to left half plane is k. 3) Choose the frequency range and step size of frequency range.where 1 i k ≤ ≤ .

Case Study
Simulation studies have been carried out for multi-area LFC equipped with PI controller assuming the load change of 0.1 p.u in Area 1.The system parameters of each control area are listed in Table 1.
The upper bounds of time delay in each area are preset as same value while designing the controllers.First the controller values are determined by using continuous pole placement method and then for the designed controller values, delay margin is theoretically calculated using FST.The theoretical results are presented in Table 2.
The theoretical value of delay margin is calculated as 34.612 s.That is, the controller designed using the proposed method for the preset upper bound of delay 8 s not only retains stability for time delays up to 8s, it can also ensure stability till 34.612 s (delay margin).To validate the theoretical results, simulation is performed using MATLAB/SIMULINK by increasing the delay step by step from zero until the LFC system becomes unstable.The simulation results are shown in Figure 4 from which it is observed that the delay margin of the system is 34.6 s which is very closer to the theoretical delay margin (34.612).
Similarly, the controller values of three-area LFC are determined using the proposed method for preset upper bound of time delay 10 s.The results are compared with the controller gain reported in [9] and shown in Table 3.
The delay margin of the three-area LFC system with PI controller (designed using proposed method) is theoretically computed as 21.115 s whereas the method reported in (9) can maintain the stability only up to 14 s.This   clearly shows that the proposed controller is highly robust compared to the existing method.
To validate the theoretical delay margin values, simulation is performed for three area LFC scheme keeping the preset upper bound of time delay as 10 s and the results are shown in Figure 5.The solid line shows the response of system for preset upper bound of delay 10 s and dashed line shows the response of the system for delay margin.The simulation results reveal that, the stability of the system is guaranteed for all time delays smaller than the delay margin.From the figure, it can be concluded that the delay margin of the system is 21.1s which is nearer to the theoretical delay margin 21.115 s.

Conclusion
In this paper, a robust controller based on continuous pole placement method is designed for multi area LFC scheme affected by communication delays.The controller values are determined by shifting the rightmost eigenvalues to left plane in a quasi continuous way for any particular preset upper bound of time delay.The proposed controller is highly robust in sustaining the stability of the system even for delays greater than the preset upper bound of time delay.Case studies have been carried out for two area and three area LFC schemes.The efficiency of the controller is validated by finding the value of time delay margin theoretically using Frequency Sweeping test.The theoretical value of delay margin is verified using simulation studies.Simulation result shows that the proposed controller gives larger stability margin.

Figure 1 .
Figure 1.Dynamic model of i th control area of multi area LFC scheme.
obtained.e j is j th unity vector. of the eigenvalues with respect to changes in feedback gain k.

Figure 2 .
Figure 2. Illustration of the movement of the real parts of the rightmost eigenvalues for each iteration.
the system goes unstable.

4 ). 5 ) 6 )
For different frequencies ω , find the absolute values of all the eigenvalues of the matrix pencil Determine the angle and frequency at which the absolute value of each eigenvalue variation reaches one.Otherwise go to step 3 and vary the frequency range.Calculate the delay margin

Figure 5 .Table 3 .
Figure 5. Frequency deviation of three area LFC for preset time delay 10 s.(a) Area 1, (b) Area 2, (c) Area 3.Table3.Controller parameters of three area LFC system for preset upper bound of time delay 10 s.

Table 2 .
Controller parameters of two area LFC system or preset upper bound of time delay 8 s.