A Multi-Model Approach to Design a Robust SVC Damping Controller Using Convex Optimization Technique to Enhance the Damping of Inter-Area Oscillations Considering Time Delay

This paper introduces a multi-model approach to design a robust supplementary damping controller. The designed fixed-order supplementary damping controller adjusts the voltage reference set point of SVC. There are two main objectives of the controller design, damping low frequencies oscillations and enhancing power system stability. This method relies on shaping the closed-loop sensitivity functions in the Nyquist plot under the constraints of these functions. These constraints can be linearized by choosing a desired open-loop transfer function. The robust controller is designed to minimize the error between the open-loop of the original plant model and the desired transfer functions. These outcomes can be achieved by using convex optimization methods. Convexity of the problem formulation ensures global optimality. One of the advantages of the proposed approach is that the approach accounts for multi-model uncertainty. In contrast to the methods available in the literature, the proposed approach deals with full-order model (i.e., model reduction is not required) with lower controller order. The issue of time delay of feedback signals has been addressed in this paper for different values of time delay by applying a multi-model optimization technique. The proposed approach is compared to other existing techniques to design a robust controller which is based on H2 under pole placement. Both techniques are applied to the 68-bus system to evaluate and validate the robust controller performance under different load scenarios and different wind generations.


Introduction
The power system grid has increased rapidly, something which has added additional challenges to the process of reliable power transfer between interconnected systems in a large power network.The large-scale penetration of intermittent renewable energy increases uncertainty and variability to power system operation.For secure operation of power systems under variable conditions, a power system's damping controllers must be robust.Electromechanical oscillations in the range of 0.2 Hz to 1 Hz are categorized as inter-area modes.These modes primarily arise due to weak interconnections characterized by long transmission lines between different operating areas of an interconnected power system.One of the main challenges in the secure operation of an interconnected power system is the damping of these inter-area modes [1].System stability could be affected without adequate damping of these low frequency oscillations.
Event such as the 1996 western interconnect blackout is an example of one such event.
Recently, Flexible AC Transmission System (FACTS) devices have become widely used in power systems.The main purpose of using these devices is to increase the capability of transferred power between interconnected areas and enhance the voltage profile.A Static Var Compensator (SVC) is a shunt FACTS device which injects reactive power to maintain the voltage at a point of connection within a particular range to enhance system stability.Controlling SVCs helps to damp inter-area oscillations.A supplementary signal could be added to adjust the voltage reference set point of SVC to achieve the desired damping, [2] [3].The location of SVCs for damping inter-area oscillation is important and the SVCs are usually placed at either end of a tie-line, [4] [5].Depending on system configuration, multiple SVCs may be needed to improve overall system damping.
Recently, researchers have investigated H2, H ∞ optimizations [6] [7] [8] and µ-synthesis [9] in power systems to design a robust controller.These approaches are mainly used to improve power system stability and damp power system oscillations.Riccati's equation is used in these approaches to solve the H ∞ control design problem.Recently, the linear matrix inequalities (LMIs) approach has been used to solve this issue and has provided better performance.The new techniques are presented in [10] [11] and show how to design a robust controller for multi-model uncertainty using H 2 and H ∞ under pole placement; however, these techniques require reducing the order of the plant model.
The challenges of the existing approaches are as follows: 1) These approaches are based on reducing system order.The model reduction is the process of reducing the order of a given plant model in which the performance of the reduced model is identical to the original one [12] [13].As a result, there is a loss of information.The level of this information loss depends on the order to which the system is reduced and the method used to reduce the information.The proposed method does not require any model order reduction.
In addition, model order reduction is an O(n^3) operation and the computing model order reduction for large systems is computationally expensive.
2) The order of the controller based on existing approaches is usually high for large systems [4]- [10] since it is the sum of the orders of the reduced plant model plus the order of the weights filter as mentioned in [14].For example, in reference [15] the order of the controller is 10 and it is 7 in reference [10].Despite using an unreduced plant model, the proposed controller results in a lower order controller than existing methods.For the case study, the controller designed using the H 2 under pole placement technique resulted in a 9th order controller.
Using the proposed approach resulted in a 4th order controller with similar or better performance.
3) In most existing methods, the controller is designed based on one operating point [4]- [16], i.e., the robustness is guaranteed around this operating point, but is not guaranteed if the system operates far from this point.A power system is a non-stationary system where operating points change for every dispatch at the system operator level.Performance of the controller degrees depend on the deviation between the current operating point and the nominal operating point for which the controller was designed.
Considering these challenges, the contributions of this paper are outlined below:  The entire plant model is used and there is no need for model order reduction. The resulting controller order is less than the controller order used in other existing methods. Multi-model uncertainty is considered. Time delay of remote signals is accounted for. Convex formulation guarantees global optimal solutions while minimizing the error between open-loop and desired transfer functions.This paper is organized as follows:  Section 2 discusses controller design. Section 3 describes the test system with SVC and DFIG. Section 4 describes the step-by-step controller design procedures. Section 5 shows and discusses simulation results. Section 6 presents conclusions.

Class of Controller
A linear parameterized controller can be represented as: [ ] and i ρ is the controller parameters, n is the number of controller parameters and ( ) i s ϕ is a basis function.One of the common basis functions is the Laguerre basis function given in (2) [17].
where 0 ζ > .It can be shown that for any finite order transfer function F(s), ar- bitrary Laguerre parameter 0 ζ > and an arbitrary constant Any finite order stable transfer function can be approximated using controller parameterization in (1) with a desired level of accuracy by changing the number of controller parameters n.A good approximation of F(s) can be achieved for a given controller order if the choice of ζ is proper.For more details of choosing the basic functions see [17].
To obtain a convex parameterization of a fixed order controller, the linearly parameterized controller in (1) is used.The reason for the use of ( 1) is that all the points of the open-loop transfer function ( ) , L jω ρ on the Nyquist plot can be written as a linear function of the controller parameters ρ as given in (4).

Multiplicative Uncertainty
Multiplicative uncertainty is represented in (5).Suppose that ( ) G jω is the normal plant frequency response and the actual plant that describes the normal plant with uncertainty is ( ) G jω as shown in Figure 1 and (5) [18].
where ( ) is an unknown stable transfer function with 1 ∞ ∆ < .

Robust Stability and Performance.
The closed loop system in Figure 1 can be represented by:  Let the open loop transfer function ( ) ( ) ( )  , and the sensitivity function ( ) ( )  be defined.It can be seen from ( 6) that ( ) T jω defines the relationship between the reference and the output signals.

( )
S jω defines the relationship between the reference and the error.These transfer functions define the main characteristic of the closed loop architecture.
The Nyquist diagram has been used to derive the criteria of robust performance as well as the robust stability.The point ( ) shown in Figure 2 is known as the critical point that is used to study the closed-loop system stability.The circle centered at the critical point ( ) ( )  the following is obtained: The normal performance condition of stable system can be given in the following standard form [18].
( ) ( ) To define the condition of the robust performance of the system given in Fig- ure 2, substitute (5) in to (8), as given in ( 9): ( ) Since from the above equation.This constraint is required for robust performance, and, by rearranging this constraint, the standard form of the robust performance is obtained as provided in ( 10):

The Proposed Approach
The constraints in (10) satisfy the robust stability as well as robust performance.These constraints are represented in the Nyquist plot.Robustness can be achieved by using a set of convex constraints on the frequency domain.The controller can be designed based on a convex optimization problem.The solution to this issue is to reduce the norm of the distance between the actual ( ) L jω ρ and desired ( ) L jω open-loop transfer function as shown in Figure 3 [17].
Multiplying (10) by ( ) The constraints in (11) are non-convex and ( ) L jω are used to linearize these constraints.By making the problem convex, global optimality is ensured.Line d is introduced as shown in Figure 3 which is tangent to the performance disc centered at ( ) and orthogonal to the line links at the center of the performance disc to L jω ρ has to be on the right side of line d for all frequencies as shown in Figure 3.
Line d is a straight line in the complex plane and can be represented by an infinite number of points.Each point in the complex plane has a real part x and imaginary part y.The equation of the straight line d is a function of ( ) L jω and 1 W .The line can be written at each point as:

Uncertainty circle
The critical point By substituting the ( ) ( ) cos α into the Equation ( 12), the follow- ing is obtained: The linear constraints of line d that exclude the performance disk are given in ( 14): The linear constraints in ( 14) can be written in a simpler way using the facts: The constraints in (14) become: where ( ) To satisfy the condition in (15) for a set of uncertainty models, the circle centered at ( ) L jω ρ should be approximated by a polygon with 2 v > vertices.
To satisfy the robust uncertainty in (10), all the vertices of the polygon located at -1 the uncertainty disk have to be on the right side of the line d.This condition can be represented by the linear constraints of ( 16): where It is observed that the number of the linear constraints is multiplied by v.
There is another way to satisfy the robust condition in (11) by increasing the radius of the circle ω ω ρ which leads to the following convex constraints: Considering these examinations, the quadratic optimization problem can be expressed as given in ( 19): { } For multi-model uncertainty case, the constraints in (18) can be repeated for all the plant models ( ) i G jω for 1, , i m =  .The constraints in (18) still can be used if the uncertainty weights filter 1 W , 2 W and the desired open-loop trans- fer function, di L , is different for each plant model since these constraints are convex with respect to ( ) i G jω for multi-model uncertainty.

Test System
The IEEE16 machines, 68-bus system is used in this paper.This test system is particularly suited for small signal stability studies.For instance, reference [14] uses the same test system for damping inter-area modes.There are five distinct areas in the test system with a total load of 18.23 GW.Areas NETS and NYPS are interconnected through two parallel tie-lines.Figure 4 shows the single line diagram of the test system.Parameters of the generators, exciters, governors, and transmission lines of the test system can be found in [14].
Power System Toolbox (PST) [19] is used to simulate the test system including the SVC and doubly-fed induction generator (DFIG).The controller was implemented in MATLAB based on the proposed approach and has been integrated in PST.To include renewable generation, a 500 MW wind farm is placed in area 2 at bus 39 as presented in Figure 4.The wind farm is installed to add more variability to the system due to the continuous change of the output power of the wind farm.The third order model of a DFIG is used.The dynamic model of the DFIG contains a set of differential algebraic equations that has been integrated in PST [20].A single model of DFIG is used to represent the wind farm.

Static Var Compensator
The block diagram of SVC is shown in Figure 5(a).The test system has an SVC installed at bus 50.Parameters of SVC are given in Table 1.The objective of designing the controller is to damp tie line oscillations by providing additional signal to the set point of the SVC.Control structure of the proposed approach is represented as shown in Figure 5(b).

Controller Design Procedure
In this section, the step by step procedure and rationale used in designing the controller is described in detail.

Selecting Inter-Area Modes
For the given test system, under nominal operating condition, two eigenvalue pairs have damping less than 5%.In fact, one of the eigenvalue pairs has damping very close to zero, hence, the system is close to instability point.Based on the eigenvalues for nominal operating point, the inter-area modes that need to be damped for the case study are listed in Table 2. Figure 6 shows the damping ratios for frequencies of interest corresponding to the nominal operating condition.

Selecting Input/Output Signal
Appropriate selection of the input signal for the designed controller is highly essential to guarantee that the inter-area eigenvalues are controllable and observable.The controllability metric is used to choose the best effective input signal to damp the inter-area modes.The controllability metric is defined as the amount of displacement that a pole would undergo due to small change in the feedback gain and this is given in (20) [21].Using the controllability metric as shown in Figure 7, the active power flow of the line 42 to 52 is found as the most controllable measurement to damp the inter-area modes.Therefore, the input signal that feeds the controller is provided from the tie-line (42 to 52), which connects areas 4 and 5.The controller output is used as additional control signal to the SVC.

Choice of Operating Points
A power system is a non-stationary system where a set of new dispatches are computed every 5 to 15 minutes.As a result, the total number of possible operating points are innumerable, hence, six different operating points which represent several stress levels of the system are used for controller design and validation.Stress levels of the system in this context are quantified using eigen-spectrum.Eigenvalues convey two very important attributes: oscillation frequencies and their corresponding damping ratio.Damping ratio illustrates how much energy is dissipated during each cycle for a given frequency.Six different operating points are created where the damping ratio of the eigenvalues that correspond to inter-area mode of the system are progressively made worse.The system has been extensively studied and these operating points listed in Table 3 are considered for this study as they greatly affect the inter-area modes.The generators G15 and G16 are adjusted to obtain different operating points.In addition, wind generation also varies between different operating points.All the values in Table 3 are in per-unit system.

Desired Open Loop Transfer Function (Ld)
Selecting d L is based on design specifications.d L normally has a high am- plitude in low frequencies for good tracking and the system follows the reference signal.At high frequencies, d L should have a small amplitude to provide ro- bustness and noise rejection characteristics.
( ) d L s could be chosen as c s ω where c ω is the desired closed-loop bandwidth [17] [22].Typically, the band- width is the range of frequencies for which the gain is significant.Generally, a high bandwidth will allow for a faster response.In the case study, the aim is to damp the inter-area modes in the range of frequencies (0.2 -1.0 Hz), so a bandwidth of more than 2π 2π 1.0 6.28 rad sec f = * = is needed.The desired bandwidth c ω should be more than 6.28 rad/sec.For the case study, there is a Energy and Power Engineering  strong resonance mode around 1 ω and 2 ω as shown in Figure 8 and these modes should be cancelled by the controller.So c ω is selected as

Weighting Filters (W1 and W2)
Selection of 1 W and 2 W is essential for the controller design.In this work, 1 W is designed as a first-order low-pass filter to gain a good disturbance rejec- tion. 2 W is designed as a high-pass filter to guarantee robustness and minimize the controller effort in high frequencies.Frequency response of 1 W and 2 W is shown in Figure 9.

H2 Controller under Pole Placement
For comparison, a damping controller is designed using pole placement and H 2 optimization following two steps based on matrices [11].First, a state feedback controller is developed that uses the system states to generate a control signal.This is achieved by solving a set of Linear Matrix Inequalities (LMIs) that places the system poles into a cone area in the complex plane, while minimizing the amplitude of control signal represented by its H 2 norm.Thereafter, a state estimator is developed that constructs system states from the output.A similar set of LMIs is employed for this purpose.The controller can be obtained by a transfer function equivalence of the state-feedback controller and the state estimator combined.This approach is considered multi-model so the controller is designed based on different load conditions.However, this approach still suffers from drawbacks 1 and 2 listed in section 1.For the case study, the damping ratio is set to be 10% as the boundary of the pole placement region.Also, the weights filter is selected to be the same as the weight filter used in the proposed approach.
The same operating points listed in Table 3 are used to design the controller using this approach.Plant/system model needs to be reduced based on this approach in such a way that the response of the reduced system is similar to that of the original system in the frequency range of interest.The test system consists of 190 states including the DFIG and the SVC.For the frequency range of interest, the plant model can be reduced to at least 7th order.In addition, the total order of the controller based on [14] is equal to order of the reduced system plus the order of weighting filters.In this case, this equates to a controller order of 7 + 2 i.e., 9 states.Figure 10 shows the original and the reduced plant model and the figure shows that these models are identical in the frequency range of interest.No model order reduction is required for the proposed method.Using the proposed approach, a 4th order controller is designed which replicates the Figure 10.Frequency response of the original and the reduced system, OP1.frequencies of interest.This approach is applied to design a robust controller to compare it with the proposed method.Interested readers are referred to [11] for more details.

Results and Discussion
In this section, a two-part validation of the proposed approach is presented.
Comparisons are drawn for the proposed method for the numerical and time-based domains, with the base case only accounts for SVC and the methods presented in section 4. In the first part of validation, eigenvalue spectrum obtained using the different methods are compared.Specifically, comparisons for damping ratios are drawn for different modes of interest.
In the second part, time domain performance results are presented.Comparing controllers which have user-defined parameters is not straight forward.A given control methodology and design can result in a poorly performing controller based on bad parameter selection.One way to avoid this situation is to use a standard set of values for user defined parameters and use the same parameters in both approaches.The same parameters used in [11] [23], and [24] are used for H 2 under pole placement controller such as the damping ratio of the boundary of the pole placement region.The same controller operating points and weights are used in both approaches.Furthermore, the focus of the validation process is not to provide a quantitative comparison, but rather to show that the proposed controller gives comparable results to those of existing methods with the advantages listed in section 1.

Eigenvalue Analysis
Eigenvalue comparisons of the proposed controller with base case (i.e., with only the SVC for six different operating points) are given in Table 4. Substantial improvements in damping ratio are seen with the proposed controller.For instance, when mode 1 with operating point 4 is considered, the system has no controller and the damping ratio is negative.These circumstances lead to an unstable system.
With the addition of the proposed controller, the damping ratio is improved to (0.1814) from (−0.0008).A similar trend of improved damping ratio is seen across all six operating points.The modes of the test system under different load conditions are shown in Figure 11.

Robustness under Variability in Load Conditions and Wind
Generation To investigate the robustness of the proposed controller, a three-phase fault is placed at different areas with different operating points.Application of a fault in power systems results in differences between mechanical and electrical power which results in electromechanical oscillations.The tests used for validating A. Abdlrahem  Comparisons between the system with and without the proposed controller show that the maximum overshoot and damping are considerably better with the addition of the proposed controller under all three operating points tested for scenario one.Of particular note are the comparisons for operating point 4. Without the proposed controller, the system becomes marginally unstable as shown in Figure 12(c) and Figure 12(f).The addition of the proposed controller not only makes the system stable but also damps out oscillations fast.Both the H 2 under pole placement and the proposed controller have similar performance; however, in some cases the proposed controller has slightly better damping.
In scenario two, a 50 ms fault is applied at bus 49 in area 2. This results in significant drop in tie-line flow through lines 42 through 52 during the fault as can be seen in Figure 13(a)-(c).This scenario captures the performance of the proposed controller as the fault is applied relatively close to the SVC.Angular separation between areas 2 and 5 (i.e., between generators G16 and G10) is shown in Figure 13(d)-(f).

Time Delay
The major problem with remote signal use is the time delay and the range of the time delay which varies depending on various factors such as the distance of the  remote signals.These signals can be delayed up to 100 ms.Therefore, it is very important to account for an uncertain time delay to ensure the robustness under various conditions.In this paper, a multi-model optimization method is used to include the effect of time delay.In the previous section, no time delay is considered since the main concern was to show that the method works for different scenarios and is comparable with the existing approach.The time delay in this section has been approximated by the second order Pade approximation [23].To design a robust controller based on multi-model optimization approach for the uncertainty in time delay, the time delay incorporates the worst case (OP3).Operating points 1 and 2 are also chosen to design the controller.The new controller is designed based on three operating points using the same procedure in section 4. Figure 14 shows the black diagram for incorporating the time delay with the feedback signal.hand, the second controller can maintain system stability and is able to damp these oscillations as seen in Figure 16(b) and Figure 16(d).

Conclusion
A multi-model approach is used in this paper to design a robust supplementary damping controller.The designed fixed-order supplementary damping controller provides a supplementary signal to the voltage reference set point of SVC.
The main objectives achieved in this paper are damping low frequencies oscillations and enhancing power system stability.The controller design relies on shaping the closed-loop sensitivity functions in the Nyquist plot under the constraints of these functions.scenarios are run where disturbances are applied to different areas of the test system under different operating conditions.In all tested cases, the proposed controller significantly improved the system's dynamic response and compared favorably with an existing control technique H 2 under pole placement.Improved controller performance with a lower order controller and without the need for model order reduction are the primary advantages of the proposed method that has been validated using both numerical and time domain analysis.The issue of delaying the feedback signal has been addressed using multi-model optimization.

)
Energy and Power Engineeringwhere, G is a scalar function used in the case of a single model and in the case of a multi-model controller, a set of the multi-model uncertainty can be defined as represents the open-loop transfer function for the i-th model.

6 )
Energy and Power Engineering

Figure 1 .
Figure 1.Black diagram representing an uncertain feedback system.

1 W 1 W
jω is known as the performance disc.The uncertainty disc is represented by the circle with radius stability is achieved if, and only if, the uncertainty disc centered at the original open-loop transfer function with radius jω on the Nyquist plot.
. A sufficient condition for constraints in(11) is that the circle centered at the actual open-loop transfer function

Figure 4 .
Figure 4. Single line diagram of 68 bus test system.
Energy and Power Engineering

Figure 5 .
Figure 5.The block diagram of (a) SVC and (b) control representation.

Figure 6 .
Figure 6.Damping ratios and frequencies of eigenvalues for OP1, normal operating point.

Figure 7 .
Figure 7. Controllability indices of controllable eigenvalues based on selecting the line 42 to 52.
(rad/s) Original plant model Reduced plant model Energy and Power Engineering

Figure 11 .
Figure 11.Modes of the test system under three different operation points.(a) Modes of the test system, OP1; (b) modes of the test system, OP3; (c) modes of the test system, OP4.

Figure 15 shows
Figure15shows the dynamic response of the test system with the designed controller for different values of time delay.It can be seen that, the controller is able to damp the power system oscillations under variety of operating points and time delay values.A comparison between the controller designed in section 4 (without incorporating the time delay) and the new controller designed based on incorporating the time delay is shown in Figure16.Both controllers behave almost the same when the feedback signal is delayed by 200 ms.However, the first controller is not able to damp the inter-area oscillations and the feedback signal is delayed by 300 ms as shown in Figure16(b) and Figure16(d).On the other

Figure 14 .Figure 15 .
Figure 14.Black diagram of output signal time delay.
300ms time delay Energy and Power Engineering

Figure16.
Figure16.Dynamic response of the test system with the two controllers under different time delay.(a) Tie-line power, fault at 41, OP 3200 ms; (b) tie-line power, fault at 41, OP 3300 ms; (c) tie-line power, fault at 41, OP 6200 ms; (d) tie-line power, fault at 41, OP 3300 ms.

Table 2 .
Eigenvalues, damping ratios, and the frequencies of the inter-area modes of the test system.

Table 3 .
Different operating points for 68 bus system.
*It is not used in the control design, it is used to validate the controller.

Table 4 .
Damping and frequencies of the inter-area modes under different load conditions of the 68-bus system.
The IEEE 68-bus system with wind farm is used to demonstrate the controller performance.Test scenarios are designed to emulate real life scenarios seen at system operator level; specifically, uncertainties in operating conditions and changes to system topology are considered.Several test