Next Article in Journal
Towards Net Zero Carbon Economy: Improving the Sustainability of Existing Industrial Infrastructures in the UK
Previous Article in Journal
A Survey on Cybersecurity Challenges, Detection, and Mitigation Techniques for the Smart Grid
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Delay-Dependent DOFC for Damping Inter-Area Low-Frequency Oscillations in an Interconnected Power System Considering Packet Loss of Wide-Area Signals

1
School of Automation, Central South University, Changsha 410004, China
2
Hunan Xiangjiang Artificial Intelligence Academy, Changsha 410004, China
*
Author to whom correspondence should be addressed.
Energies 2021, 14(18), 5892; https://doi.org/10.3390/en14185892
Submission received: 26 June 2021 / Revised: 10 September 2021 / Accepted: 12 September 2021 / Published: 17 September 2021

Abstract

:
In this paper, the inter-area low-frequency oscillations are restrained in the interconnected power system by making use of the delay-dependent wide-area dynamic output feedback controller (DOFC). Modal analysis is adopted to obtain the modes of inter-area oscillation to be controlled and the Schur truncation model reduction technique is represented to reduce the order of the power system. The augmented closed-loop system model, where the transmission delay and packet loss of wide-area signals are considered, is established. The sufficient conditions of exponentially mean-square stable are obtained according to Lyapunov’s stability theory. Finally, case studies are carried out on a two-area four-machine power system, where our proposed controller, a conventional controller, and the wide-area damping controller in the existing references are installed, respectively. The simulation results under different external disturbances, packet loss rates, and delays are presented to show the effectiveness and advantages of our proposed controller.

1. Introduction

In recent years, with the increase of the scale and load ability of power systems, low-frequency oscillations can easily occur due to insufficient damping, which may threaten the stability and safety of the power grid and cause system islanding or even blackouts [1,2].
The conventional power system stabilizers (PSSs) have been used as the mainstream to damp the local low-frequency oscillations [3,4]. These PSSs can control the local modes well, but they are less effective at restraining inter-area modes because such modes are not observed or controlled directly from local signals of the generators [5].
In the modern power system, the wide-area measurement system (WAMS) makes the remote measurements available by using phasor measurement units (PMUs) [6,7], and with these measurements, some effective measures can be put forward to damp the inter-area modes of the power system [8,9]. Further, wide-area PSS (WAPSS) based on wide-area signals from WAMS can effectively damp these inter-area oscillation modes [10,11,12]. However, these signals are inherently subject to communication time delays which can adversely affect the performance of WAPSSs [13]. Due to the inherent reliability issues and unexpected factors related to the communication network, the wide-area signal transmission often suffers from fading and congestion, which will lead to delay and packet dropping [14].
Some methods have been used to deal with the problem of time delay. Pade approximation is applied to obtain the linear model with a certain order whose influence on the power system is approximately equivalent to the time delay [15]. It can be seen that the amount of calculation will be greatly increased because a higher order is adopted to reach better accuracy in the approximation. The method to compensate for time delay with the delay estimation is presented in [16,17,18], which, however, is sensitive to the system’s change. Some robust control methods are employed in the design of wide-area controllers to maintain the power system stability within an exact time delay [19,20,21,22]. The predictive control schemes are used for variable time-delays in [23]. It is worth noting that the packet loss problem may occur in remote signal transmission, which may disable the controllers and degrade the damping performance [24].
Some concerns and comments on packet loss can be seen in [25,26,27,28]. The influence of the packet loss on the interconnected power system is discussed in [25], but no strategy is presented to deal with this problem. The observer-based controllers in [26,27] are designed to restrain inter-area low-frequency oscillation in the case of packet loss, and the delay caused by the transmission of wide-area signals is ignored. Wang and al., in [28], propose a state feedback controller to damp the interconnected power system considering packet loss and time delay.
It is known that output feedback is a common feedback method which is easy to realize in technology, especially because dynamic output feedback has better properties and control effect. In addition, the PSSs controllers are used for feedback control. Therefore, wide-area DOFC has garnered a lot of interest. In [29], a wide-area output feedback damping controller is put forward and a delay-dependent dynamic output feedback control strategy is proposed in [30] to improve damping and enhance stability. However, the impact of packet loss has not been investigated. Although the DOFC is proposed to damp the inter-area oscillations in case of packet loss in [31], the way to obtain the parameters of the controller is not presented.
Inspired by the above investigation, the delay-dependent wide-area DOFC is put forward to restrain the inter-area oscillation with the consideration of the packet loss of wide-area signals in this paper. The rest of this paper is arranged as follows. The augmented power system is set up by considering time delay and packet loss in Section 2. Section 3 represents the proposed strategies that include the delay-dependent criterion of exponentially mean-square stable and the method to obtain the parameters of the DOFC. Case studies are shown in Section 4 to verify the effectiveness and advantage of the proposed method. Finally, the conclusions are drawn in Section 5.
The following contributions are made in this paper:
  • An augmented closed-loop system model considering time-delay and packet loss of remote signals is set up to deal with inter-area oscillations.
  • The sufficient conditions of exponentially mean-square stable for closed-loop power system with packet loss of remote signals are derived.
  • The method to obtain parameters of DOFC is given.
  • The relationship between delay margin and the packet loss rate is presented, and the effectiveness and advantage are verified by simulation results.
Notations: The notations used here are fairly standard except where otherwise stated. The superscripts “−1” and “T” stand for the inverse and transpose of a matrix, respectively; R n × m is the set of all n × m real matrices; R n denotes the n-dimensional Euclidean space; I is an identity matrix with appropriate dimension; d i a g   { } denotes a block-diagonal matrix;   P > 0 means that the matrix P is positive-definite; λ m a x ( A ) and λ m i n ( A )   mean the largest and smallest eigenvalue of the matrix A.   ·   is the spectral norm; p r {   ·   } means the probability and E {   ·   } denotes the expectation; the symbol “*” within a matrix represents a term that is induced by symmetry.

2. Control Problem Formulation

2.1. Modelling of Power System

The following state-space model can be obtained by linearizing power system:
{ x ˙ 0 ( t ) = A 0 x 0 ( t ) + B 0 u 0 ( t ) y 0 ( t ) = C 0 x 0 ( t )
where x 0 ( t ) R n × 1 , u 0 ( t ) R m × 1 , and   y 0 ( t ) R p × 1   are the states, input, and output vectors of the system, respectively. A 0 R n × n , B 0 R n × m , and C 0 R p × n are the state, input, and output matrices with appropriate dimensions, respectively.
The order of the linearized model is comparatively high in a large-scale power system, which will make the controller design difficult or even infeasible. Hence, it is not necessary to use a full model for analyzing the system and designing a controller. The low-frequency oscillation occurs over the frequency range from 0.1 Hz to 2.5 Hz. A reduced-order model can be employed to guarantee the required frequency range. Thus, the Schur model-reduction strategy has been applied to reduce the order of the power system and the details are also seen in [32]. The reduced-order model of power system is given as follows:
{ x ˙ ( t ) = A x ( t ) + B u ( t ) y ( t ) = C x ( t )
where u ( t ) , y ( t ) ,   and x ( t ) are the input, output, and state vectors, A, B, and C are the system, control, and output matrices of the reduced-order power system, respectively.

2.2. The Augmented Closed-Loop Power System

Figure 1 represents the scheme of closed-loop power system, including the power system, our proposed wide-area controller, and delay time during the transmission of wide-area signals. In addition, the packet loss of remote signals will occur during the transmission. If the power system experiences packet loss, the previous data will be used. Thus, we have:
y 2 ( k ) = { y 1 ( k 1 d τ s )      experience packet loss y 1 ( k d τ s ) does not experience packet loss
where, d   is the transmission delay and τ s is the sampling time.
The following events can describe the packet loss of wide-area signals:
γ i ( i = 1 , 2 , , l ) = {   1               successful   transmission   of   i t h   remote   signal 0               unsuccessful   transmission   of   i t h   remote   signal
There is a certain probability of packet loss, which may not happen. The probability of packet loss can be described as follows:
{   p r {   γ i = 0 } = ρ i   p r {   γ i = 1 } = 1 ρ i
where, ρ i is called packet loss rate.
Furthermore, it will be yielded that:
y 2 ( k ) = ( I γ ) y 1 ( k 1 d τ s ) + γ y 1 ( k d τ s )
where γ = d i a g ( γ 1 ,   γ 2 , , γ l ) .
Obviously, y 1 ( k ) is the discrete value of y ( t ) and u ω ( t ) is the continuous value of y 2 ( k ) . Their relationship can be given as:
{ u ( t ) = y ω ( t ) u ω ( t ) = ( I γ ) y ( t τ s d ) + γ y ( t d )
The controller in type of dynamic output feedback control is given as follows:
{   x ˙ ω ( t ) = A ω x ω ( t ) + B ω u ω ( t ) y ω ( t ) = C ω x ω ( t ) + D ω u ω ( t )
where, u ω ( t ) ,     y ω ( t ) and x ω ( t ) are the input, output and state vectors of the controller. A ω ,   B ω , C ω ,     and   D ω are the parameters of the proposed controller to be designed.
Hence, the augmented closed-loop power system can be represented as:
z ˙ ( t ) = A 1 z ( t ) + A 2 z ( t τ 1 ) + A 3 z ( t τ 2 )
where,
A 1 = [   A 0   B C ω A ω   ] ,   A 2 = [   B D ω γ C B ω γ C 0 0   ] ,   A 3 = [   B D ω ( I γ ) C B ω ( I γ ) C 0 0   ] , z ( t ) = [   x T ( t )   x ω T ( t )   ] T   ,   τ 1 = d ,     τ 2 = d + τ s ,   γ = d i a g ( γ 1 , γ 2 , , γ l ) .
Definition 1. 
For constants a > 0 and b > 0 , system (8) is exponentially mean-square sta ble (EMS) if the following condition holds:
E { ψ 2 } a e b t sup τ 2 s , v 0 E { ψ ( s ) 2 + ψ ˙ ( v ) 2 } ,   t 0
Definition 2. 
Define the infinitesimal generator L of V ( x t ) as follows:
L V ( x t ) = lim Δ 0 + 1 Δ E { ( V ( x t + Δ ) | x t ) V ( x t ) }

3. Main Results

Taking an infinitesimal operator L from (9), the event i can be represented by the packet loss rate in LMI (8). Thus, the sufficient conditions of the exponentially mean-square stable for the closed-loop power system (7) are discussed below.
Theorem 1.
For given positive number d , matrices   A ω ,   B ω ,   C ω ,   D ω , if there are positive definite matrices P ,       R j ( j = 1 ,   2 ,   3 ,   4 ) ,     Q j ( j = 1 , 2 )   and any matrices N ,     M ,     ρ n with appropriate dimensions, such that the following LMIs are satisfied, the closed-loop power system (7) is said to be exponentially mean-square stable:
[   Ξ 11 Ξ 12 Ξ 22   ] < 0
where,   Ξ 11 = [   Q 1 + Q 2 0 Q 1 0 0 Q 2 P 0 0 Δ ] + F W f + W f T F , Ξ 12 = [   τ 1 N    ( τ 2 τ 1 ) M τ 2 U   ] , Ξ 22 = d i a g [ τ 1 R 1 ( τ 1 τ 2 ) R 2 τ 2 R 3   ]
Δ = τ 1 R 1 + ( τ 2 τ 1 ) R 2 + τ 2 R 3 , F = [   N M U V   ] , ρ n = d i a g ( ρ 1 , ρ 2 , , ρ l )
W f = [ I 0 I A 1 I I 0 A ^ 2 0 I I A ^ 3 0 0 0 I ] , A ^ 2 = [ B D ω ρ n C 0 B ω ρ n C 0 ] , A ^ 3 = [ B D ω ( I ρ n ) C 0 B ω ( I ρ n ) C 0 ] .
Proof. 
Choose the following Lyapunov–Krasovskii function:
V ( z ( t ) ) = i = 1 3 V i ( z ( t ) )
where,
V 1 ( z ( t ) ) = z T P z ( t ) ,   V 2 ( z ( t ) ) = t τ 1 t z T ( s ) Q 1 z ( s ) d s + t τ 2 t z T ( s ) Q 2 z ( s ) d s ,                                                                                                                                         V 3 ( z ( t ) )                                               = t τ 1 t s t z ˙ T ( v ) R 1 z ˙ ( v ) d v d s + t τ 2 t τ 1 s t z ˙ T ( v ) R 2 z ˙ ( v ) d v d s                                                                                                                                                                               + t τ 2 t s t z ˙ T ( v ) R 3 z ˙ ( v ) d v d s .
The infinitesimal generator L of (9) is obtained as:
{ L V 1 ( z ( t ) ) = 2 z T P z ˙ ( t )    L V 2 ( z ( t ) ) = z T ( t ) Q 1 z ( t ) z T ( t τ 1 ) Q 1 z ( t τ 1 )   + z T ( t ) Q 2 z ( t )   z T ( t τ 2 ) Q 1 z ( t τ 2 ) L V 3 ( z ( t ) ) = z ˙ T Δ z ˙ ( t ) t τ 1 t z ˙ T ( s ) R 1 z ˙ ( s ) d s   t τ 2 t τ 1 z ˙ T ( s ) R 2 z ˙ ( s ) d s t τ 2 t z ˙ T ( s ) R 3 z ˙ ( s ) d s
For any suitable dimensioned matrices N, M, and U, the following equations will satisfy by virtue of the Newton–Leibniz formula:
{ 2 δ T ( t ) N [ z ( t ) z ( t τ 1 ) t τ 1 t z ˙ ( s ) d s ] = 0 2 δ T ( t ) M [ z ( t τ 1 ) z ( t τ 2 ) t τ 2 t τ 1 z ˙ ( s ) d s ] = 0 2 δ T ( t ) U [ z ( t ) z ( t τ 2 ) t τ 2 t z ˙ ( s ) d s ] = 0 2 δ T ( t ) V [ A 1 z ( t ) + A 2 z ( t τ 1 ) + A 3 z ( t τ 2 ) z ˙ ( t ) ] = 0
where δ ( t ) = [ z T ( t ) , z T ( t τ 1 ) , z T ( t τ 2 ) , z ˙ ( t ) ] T .
Then, we will obtain:
{ 2 δ T ( t ) N t τ 1 t z ˙ ( s ) d s τ 1 δ T ( t ) N R 1 1 N T δ ( t ) + t τ 1 t z ˙ T ( s ) R 1 z ˙ ( s ) d s   2 δ T ( t ) M t τ 2 t τ 1 z ˙ ( s ) d s ( τ 2 τ 1 ) δ T ( t ) M R 2 1 M T δ ( t ) + t τ 2 t τ 1 z ˙ T ( s ) R 2 z ˙ ( s ) d s 2 δ T ( t ) V t τ 2 t z ˙ ( s ) d s τ 1 δ T ( t ) V R 3 1 V T δ ( t ) + t τ 2 t z ˙ T ( s ) R 3 z ˙ ( s ) d s
By substituting (11) and (12) into (10), we can get the following inequality:
L V ( z ( t ) ) δ T ( t ) [ Ξ 11 + Ξ 12 T Ξ 22 1 Ξ 12 ] δ ( t )
By applying the Schur complement theorem to (13), we can obtain:
L V ( z ( t ) ) λ δ T ( t ) δ ( t )
where λ = λ m i n ( ϕ ) , ϕ = [ Ξ 11 Ξ 12 Ξ 22 ] .
Taking expectation on both sides of (14), it will be yielded that:
E { L V ( x ( t ) ) } λ E { z ( t ) 2 + z ˙ ( t ) 2 }
Now, let us define:
W ( z ( t ) ) = e ε t V ( z ( t ) )
where ε > 0 .
Its infinitesimal operator L will be given by:
L W ( z ( t ) ) = e ε t ( ε V ( z ( t ) ) + L V ( z ( t ) ) )
Integrating and taking expectation on both sides of (15), the following can be obtained:
  E { W ( z ( t ) ) W ( z ( ( 0 ) ) } = 0 t E { e ε t ( ε V ( z ( s ) ) + L V ( z ( s ) ) ) d s } 0 t E { e ε t   [ ( ε λ m a x ( P ) + ε τ 1 λ m a x ( Q 1 ) + ε τ 2 λ m a x ( Q 2 ) λ ) z ( s ) 2 + ( α ε λ ) z ˙ ( s ) 2 ] d s }
where, α = τ 1 2 λ m a x R 1 + ( τ 2 τ 1 ) 2 λ m a x R 2 + τ 2 2 λ m a x R 3 .
Select suitable positive ε to make the following inequality hold.
ε λ m a x ( P ) + ε τ 1 λ m a x ( Q 1 ) + ε τ 2 λ m a x ( Q 2 ) λ < 0 ,   α ε λ < 0
which implies
E { W ( z ( t ) ) } E { W ( z ( 0 ) ) }
Then, we can have
e ε t E { V ( z ( t ) ) } ( λ m a x ( P ) + τ 1 λ m a x ( Q 1 ) + τ 2 λ m a x ( Q 2 ) ) E { z ( 0 ) 2 } + α E { z ˙ ( 0 ) 2 }
From (19), we can obtain
E { V ( x ( t ) ) } β s u p τ 2 s , v 0 E { ψ ( s ) 2 + ψ ˙ ( v ) 2 } e ε t , t 0
where β = max { λ m a x ( P ) + τ 1 λ m a x ( Q 1 ) + τ 2 λ m a x ( Q 2 ) ,   α } .
It is observed that
V ( z ( t ) ) λ m i n ( P ) z T ( t ) z ( t )
Thus, (20) can be transformed into:
E { x T ( t ) x ( t ) } β ¯ e ε t s u p τ 2 s , v 0 E { ψ ( s ) 2 + ψ ˙ ( v ) 2 } , t 0
where β ¯ = β / λ m i n ( P ) .
It is concluded from (22) that system (7) is exponentially mean-square stable according to definition 1. This completes the proof. Next, we will discuss how to obtain the parameters of WADC.
Theorem 2. 
For preset scalars d ,   θ 1 , , θ 8 ,   ε 1 , , ε 8 , ρ , the closed-loop system (7) is exponentially mean-square stable, if there are positive definite matrices P ^ , Q ^ j ( j = 1 , 2 ) ,   R ^ j ( j = 1 , 2 , 3 ) , any suitably dimensioned matrices N ^ , M ^ , U ^ , A ^ , B ^ , C ^ , D ^ , and revertible matrix X 1 ,   X 2 , such that the following LMIs hold:
[ Ξ ^ 11 Ξ ^ 12 Ξ ^ 22 ] < 0
where,
  Ξ   ^ 11 = [ Q ^ 1 + Q ^ 2 0 Q ^ 1 0 0 Q ^ 2 P ^ 0 0 Δ ^ ] + F ^ W ^ f + W ^ f T F ^   ,     Ξ   ^ 12 =        [   τ 1 N ^ ( τ 2 τ 1 ) M ^ τ 2 U ^   ] ,                           Ξ   ^ 22 = d i a g [   τ 1 R ^ 1 ( τ 2 τ 1 ) R ^ 2 τ 2 R ^ 3   ] ,   Δ ^ = τ 1 R ^ 1 + ( τ 2 τ 1 ) R ^ 2 + τ 2 R ^ 3 ,                                    T = d i a g [ Y , Y , Y , Y ] , Y = d i a g [ x 1 , x 2 ] , F ^ = [   N ^ M ^ U ^ V ^ ] ,            S = d i a g [ Y , Y , Y ] , P ^ = Y T P Y , Q ^ 1 = Y T Q 1 Y ,   Q ^ 2 = Y T Q 2 Y ,   R ^ 1 = Y T R 1 Y ,     R ^ 2 = Y T R 2 Y ,   R ^ 3 = Y T R 3 Y ,   N ^ = T T N Y ,   M ^ = T T M Y ,   U ^ = T T U Y V = [ θ 1 I ε 1 I θ 2 I ε 2 I θ 3 I ε 3 I θ 4 I ε 4 I θ 5 I ε 5 I θ 6 I ε 6 I θ 7 I ε 7 I θ 8 I ε 8 I ] T , W ^ f = [ I 0 0 0 I 0 A X 1 0         0 I 0 0 0 I B C ^ ω A ^ ω     I 0 I 0 0 0 ( 1 ρ ) B D ^ ω ( 1 ρ ) B ^ ω 0 I 0 I 0 0 0 0 0 0 I 0 I 0 ρ B D ^ ω ρ B ^ ω 0 0 0 I 0 I 0 0 0 0 0 0 0 0 X 1 0 0 0 0 0 0 0 0 X 2 ]      
Furthermore, the parameters of controllers are given by:
A ω = A ^ ω X 2 1 ,   B ω = B ^ ω X 1 1 C 1 ,   C ω = C ^ ω X 2 1 ,   D ω = D ^ ω X 1 1 C 1
Proof. 
Pre-and post-multiply (11) with d i a g { T T , S T } and d i a g { T , S } , respectively, which yields:
[ T T Ξ 11 T T T Ξ 12 S S T Ξ 12 T T S T Ξ 22 S ] < 0
It is observed that:
T T F W f T = [   T T N Y         T T M Y         T T U Y         T T V ] × [ I 0 0 0 I 0 A X 1 0         0 I 0 0 0 I B C ω X 2 A ω X 2     I 0 I 0 0 0 ( 1 ρ ) B D ω C X 1 ( 1 ρ ) B ω C X 1 0 I 0 I 0 0 0 0 0 0 I 0 I 0 ρ B D ω C X 1 ρ B ω C X 1 0 0 0 I 0 I 0 0 0 0 0 0 0 0 X 1 0 0 0 0 0 0 0 0 X 2 ]
Let V be [ θ 1 X 1 1 ε 1 X 1 1 θ 2 X 2 1 ε 2 X 2 1 θ 3 X 1 1 ε 3 X 1 1   θ 4 X 2 1 ε 4 X 2 1 θ 5 X 1 1 ε 5 X 1 1 θ 6 X 2 1 ε 6 X 2 1   θ 7 X 1 1 ε 7 X 1 1 θ 8 X 2 1 ε 8 X 2 1 ] T , where θ 1 , ,   θ 8 and ε 1 , ,   ε 8 are given constants.
Defining. 
A ^ ω = A ω X 2 , B ^ ω = B ω C X 1 , C ^ ω = C ω X 2 , D ^ ω = D ω C X 1 ,   Q ^ j ( j = 1 , 2 ) ,   R ^ j ( j = 1 , 2 , 3 ) , P   ^ as mentioned above and substituting (25) into (24), it can be easily known that (24) is equal to (23), which completes the proof.
Remark 1. 
The system (7) is proven to be exponentially mean-square stable with mutually independent ρ 1 , ρ 2 , , ρ l in Theorem 1. Meanwhile, in Theorem 2, for simplicity, we assume ρ 1 = ρ 2 = = ρ l = ρ and ρ is the packet loss rate of the remote signal.
Remark 2. 
In Theorem 2, when the parameters θ 1 , , θ 8 ,     ε 1 , , ε 8   a n d   ρ are preset, the conditions in Theorem 2 are strict LMIs and will be solved via the LMI Toolbox in MATLAB.
Remark 3. 
If the reduced-order power system is of n   order, the augmented closed-loop system will be of 2 n order. Thus, P ^ ,     Q ^ 1 ,   Q ^ 2 ,   R ^ 1 ,   R ^ 2 ,   R ^ 3 n × n , X 1 ,   X 2 n × n and N ^ , M ^ , U ^ 8 n × 2 n , the dimensions of A ^ , B ^ , C ^ , D ^ are n × n , n × n , 1 × n , 1 × n , respectively.

4. Simulation in Four Machine Two-Area Power System

As shown in Figure 2, the test system consists of four machines in two areas, the details of which are found in [11]. Table 1 shows the results of the modal analysis. There are two local modes, which can be improved by installing local PSSs in generator 1 and 3 with detailed parameters are shown in [11]. Though the local mode is suppressed by local PSSs, the system still has an inter-area mode with a damping ratio is 0.0823 and a frequency of 0.6567 Hz. Therefore, G1 is installed by our proposed controller in the form of DOFC to restrain the inter-area mode. According to the residue approach in [18], the difference between the rotor speed of G1 and G3 is selected as the input signal of our proposed controller.
The open-loop system, where the local controller is not concluded, is of 76 orders and the reduced-order system is obtained by the Schur balanced truncation reduction technique. The frequency responses of a system with 76 orders, 9 orders and 8, orders are presented in Figure 3, and they are very close, over 0.1–2.5 Hz. Hence, the 8th order system model is used for the simulation.
For the following given numbers, θ 1 = θ 3 = θ 5 = θ 7 = ε 2 = ε 4 = ε 6 = ε 8 = 1 , θ 2 = θ 8 = ε 3 = ε 5 = 1 ,   θ 4 = θ 6 = ε 1 = ε 7 = 0 , the parameters of the controller are obtained according to Theorem 2 by virtue of the LMI toolbox in MATLAB.
A ω =                                            [ 6.770 0.954 0.193 0.257 1.106 0.249 0.242 0.302 2.260 0.051 0.814 0.528 0.896 0.040 0.060 0.235 43.921 11.434 5.0883 4.3086 9.6593 1.4041 1.3725 3.565 28.574 8.650 1.409 1.671 7.314 0.911 1.134 2.875 6.107 4.579 0.401 0.046 5.358 0.550 0.483 0.518 3.108 1.026 0.197 1.512 2.388 3.295 1.183 0.066 35.261 3.044 0.444 1.703 10.664 1.844 4.565 3.754 0.163 0.208 0.034 0.102 0.174 0.112 0.025 0.662 ] , B ω = [   5.451 7.547 5.601 4.316 2.588 5.653 3.191 10.945   ] T , C ω = [   1.077 0.348 5.304 4.733 1.435 0.182 6.297 0.103   ]   , D ω = [   0.059   ] .
The relationship between packet loss rate and delay bound is shown in Table 2. It is clear that the increase in the packet loss rate ρ will cause a decrease in delay bound d ^ , which means that the stability of the power system is degraded.
To further clarify this, the simulation is carried out under different packet loss rates where G1 is equipped with our proposed controller. The difference between angular velocity and rotor angle of G1 and G3 is selected as performance index. Three-phase short-circuit occurs at t = 1 s and is removed at t = 1.2 s . For the same delay d = 100 m s , the response caves of δ 13 and ω 13 under different packet loss rates are shown in Figure 4a and for the   same   packet   loss   rate   ρ = 0.1 , the response caves of δ 13 and ω 13 under different delays are Figure 4b.
It can be seen from Figure 4a that the system equipped with our proposed controller has adequate damping to suppress inter-area low-frequency oscillations. Even if the packet loss rate reaches 0.9, the low-frequency oscillation can be still suppressed. Moreover, there are only slow changes in the amplitude and the regulating time of the system response as the packet loss rate increases.
From Figure 4b, we can know that the effect of damping the low-frequency oscillation will become worse as the delay increases, which also shows the importance of considering delay during the design of controller. In addition, when the delay time increases to a certain value, the system will become instable.
To illustrate that our controller is more effective in the case of packet loss, the comparison simulations among our controller, a conventional controller, and WADC proposed in [20] are carried out. Case studies are carried out under the following three scenarios. Scenario1: three-phase short-circuit occurs at t = 1 s and it is removed at t = 1.2 s . Scenario 2: at   t = 1 s , the terminal voltage of G1 increases by 5%. Scenario 3: the power flow on tie-lines increases from 413 MW to 460 MW. For delay time d = 0.1 s , the system responses under different packet loss rates are shown in the figures below.
When the packet loss of wide-area signals does not occur, that is to   say , ρ = 0 , the response curves are shown in Figure 5, Figure 6 and Figure 7, corresponding to scenarios 1, 2, 3, and some key indexes are presented in Table 3. The system with our controller, WADC in [20], and the conventional controller can all ensure the stability of the closed-loop power system under different external disturbances, but our dynamic performances are slightly better than those of WADC in [20] and the conventional controller—at least not poorer than theirs.
It can be easily seen from Figure 8, Figure 9, Figure 10 and Table 4 that except for the slightly larger overshoot of ω 13 in Figure 9, our proposed controller can stabilize the system in a shorter time and have smaller overshoot under different external disturbances when packet loss occurs.
Next, we will present other simulation results with three controllers at different rates of packet loss.
It is obviously seen from Figure 11, Figure 12, Figure 13 that though all the controllers can stabilize the system when external disturbances occur, our proposed controller has a better dynamic performance and can stabilize the system in a shorter time. The higher the packet loss rate, the more obvious the advantage. Due to the limitation of space, we will not present other results in this paper. In all, our proposed controller, compared with WADC in [20] and conventional controller, can provide more adequate damping and have better dynamic performance.

5. Conclusions

The proposed wide-area DOFC can effectively damp inter-area oscillations for the power system with packet loss and time delay. The delay-dependent stability criterion exponentially mean-square stable and the techniques to obtain the delay margin and DOFC parameters with different packets loss rates are presented. The simulation results in a four-machine and two-area power system show that under different external disturbances, compared with the WADC in [18] and the conventional controller, though all the controllers can stabilize the system, our proposed controller has better performance. The higher the packet loss rate, the more obvious the advantages.

Author Contributions

Conceptualization, M.S. and Y.G. and S.S.; methodology, M.S. and Y.G.; software, M.S.; validation, M.S., Y.G. and S.S.; formal analysis, M.S. and Y.G.; investigation, Y.G.; resources, S.S.; data curation, S.S.; writing—original draft preparation, M.S. and Y.G.; writing—review and editing, M.S. and Y.G.; visualization, M.S.; supervision, supervision; project administration, M.S.; funding acquisition, M.S. and Y.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Postgraduate Innovative Project of Central South University, grant number 2020XQLH02 and The APC was funded by Postgraduate Innovative Project of Central South University.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article and supplementary material. The data presented in this study are available in Reference [20].

Conflicts of Interest

We promise that our research does not involve conflict of interest.

References

  1. Kundur, P. Power System Stability and Control; McGraw Hill Education Private Limited: New Delhi, India, 2010. [Google Scholar]
  2. Prasertwong, K.; Mithulananthan, N.; Thakur, D. Understanding low frequency oscillation in power systems. Int. J. Electr. Eng. Educ. 2010, 47, 248–262. [Google Scholar] [CrossRef]
  3. Surinkaew, T.; Ngamroo, I. Hierarchical Co-Ordinated Wide Area and Local Controls of DFIG Wind Turbine and PSS for Robust Power Oscillation Damping. IEEE Trans. Sustain. Energy 2016, 7, 943–955. [Google Scholar] [CrossRef]
  4. Zhang, J.; Chung, C.Y.; Han, Y. A Novel Modal Decomposition Control and Its Application to PSS Design for Damping Interarea Oscillations in Power Systems. IEEE Trans. Power Syst. A Publ. Power Eng. Soc. 2012, 27, 2015–2025. [Google Scholar] [CrossRef]
  5. Ray, S.; Venayagamoorthy, G.K. Wide area signal based optimal neuro-controller for a UPFC. IEEE Trans. Power Deliv. 2008, 23, 1597–1605. [Google Scholar] [CrossRef]
  6. Kim, D.-I.; White, A.; Shin, Y.-J. PMU-Based Event Localization Technique for Wide-Area Power System. IEEE Trans. Power Syst. 2018, 33, 5875–5883. [Google Scholar] [CrossRef]
  7. Fesharaki, F.H.; Hooshmand, R.A.; Khodabakhshian, A. Simultaneous Optimal Design of Measurement and Communication Infrastructures in Hierar-chical Structured WAMS. IEEE Trans. Smart Grid 2014, 5, 312–319. [Google Scholar] [CrossRef]
  8. Agnihotri, P.; Kulkarni, A.M.; Gole, A.M.; Archer, B.A.; Weekes, T. A Robust Wide-Area Measurement-Based Damping Controller for Networks with Embedded Multiterminal and Multiinfeed HVDC Links. IEEE Trans. Power Syst. 2017, 32, 3884–3892. [Google Scholar] [CrossRef]
  9. Jiang, L.; Yao, W.; Wu, Q.H.; Wen, J.Y.; Cheng, S.J. Delay-Dependent Stability for Load Frequency Control With Constant and Time-Varying Delays. IEEE Trans. Power Syst. 2011, 27, 932–941. [Google Scholar] [CrossRef]
  10. Kamwa, I.; Samantaray, S.R.; Joos, G. Compliance Analysis of PMU Algorithms and Devices for Wide-Area Stabilizing Control of Large Power Systems. IEEE Trans. Power Syst. 2012, 28, 1766–1778. [Google Scholar] [CrossRef]
  11. Yao, W.; Jiang, L.; Wu, Q.H.; Wen, J.Y.; Cheng, S.J. Delay-Dependent Stability Analysis of the Power System With a Wide-Area Damping Controller Embedded. IEEE Trans. Power Syst. 2010, 26, 233–240. [Google Scholar] [CrossRef]
  12. Yao, W.; Jiang, L.; Wen, J.; Wu, Q.H.; Cheng, S. Wide-Area Damping Controller of FACTS Devices for Inter-Area Oscillations Considering Communication Time Delays. IEEE Trans. Power Syst. 2013, 29, 318–329. [Google Scholar] [CrossRef] [Green Version]
  13. Prakash, T.; Singh, V.P.; Mohanty, S.R. A synchrophasor measurement based wide-area power system stabilizer design for inter-area oscillation damping considering variable time-delays. Int. J. Electr. Power Energy Syst. 2019, 105, 131–141. [Google Scholar] [CrossRef]
  14. Zhu, K.; Rahimi, S.; Nordström, L.; Zhang, B. Design phasor data concentrator as adaptive delay buffer for wide-area damping control. Electr. Power Syst. Res. 2015, 127, 22–31. [Google Scholar] [CrossRef]
  15. Niu, X.; Ye, H.; Liu, Y.; Liu, X. Padé approximation based method for computation of eigenvalues for time delay power system. In Proceedings of the 48th International Universities' Power Engineering Conference (UPEC), Dublin, Ireland, 2–5 September 2013; pp. 1–4. [Google Scholar] [CrossRef]
  16. Majumder, R.; Chaudhuri, B.; Pal, B.; Zhong, Q.-C. A unified Smith predictor approach for power system damping control design using remote signals. IEEE Trans. Control. Syst. Technol. 2005, 13, 1063–1068. [Google Scholar] [CrossRef] [Green Version]
  17. Yu, S.S.; Chau, T.K.; Fernando, T.L.; Iu, H.H.-C. An Enhanced Adaptive Phasor Power Oscillation Damping Approach with Latency Compensation for Modern Power Systems. IEEE Trans. Power Syst. 2017, 33, 4285–4296. [Google Scholar] [CrossRef]
  18. Meng, L.; Yong, C.; Chen, A. Wide-Area Dynamic Damping Controller Based on Robust H∞ Control for Wide-Area Power Systems With Random Delay and Packet Dropout. IEEE Trans. Power Syst. 2018, 33, 4026–4037. [Google Scholar]
  19. Zhang, C.K.; Jiang, L.; Wu, Q.H.; He, Y.; Wu, M. Delay-Dependent Robust Load Frequency Control for Time Delay Power Systems. IEEE Trans. Power Syst. 2013, 28, 2192–2201. [Google Scholar] [CrossRef]
  20. Sun, M.; Nian, X.; Dai, L.; Guo, H. The design of delay-dependent wide-area DOFC with prescribed degree of stability alpha for damping inter-area low-frequency oscillations in power system. ISA Trans. 2017, 68, 82. [Google Scholar] [CrossRef]
  21. Alghamdi, S.; Markovic, U.; Stanojev, O.; Schiffer, J.; Hug, G.; Aristidou, P. Wide-area oscillation damping in low-inertia grids under time-varying communication delays. Electr. Power Syst. Res. 2020, 189, 106629. [Google Scholar] [CrossRef]
  22. Darabian, M.; Bagheri, A. Design of adaptive wide-area damping controller based on delay scheduling for improving small-signal oscilla-tions. Int. J. Electr. Power Energy Syst. 2021, 133, 107224. [Google Scholar] [CrossRef]
  23. Molina-Cabrera, A.; Ríos, M.A.; Besanger, Y.; Hadjsaid, N. A latencies tolerant model predictive control approach to damp Inter-area oscillations in delayed power systems. Int. J. Electr. Power Energy Syst. 2018, 98, 199–208. [Google Scholar] [CrossRef]
  24. Padhy, B.P.; Srivastava, S.C.; Verma, N.K. A Wide-Area Damping Controller Considering Network Input and Output Delays and Packet Drop. IEEE Trans. Power Syst. 2016, 32, 166–176. [Google Scholar] [CrossRef]
  25. Mekki, K.; Snyder, A.F.; HadjSaid, N.; Feuillet, R.; Georges, D.; Margotin, T. Damping controller input-signal loss effects on the wide-area stability of an interconnected power system. In Proceedings of the 2000 Power Engineering Society Summer Meeting (Cat. No.00CH37134), Seattle, WA, USA, 16–20 July 2000; pp. 1015–1019. [Google Scholar]
  26. Sevilla, F.R.S.; Jaimoukha, I.M.; Chaudhuri, B.; Korba, P. A Semidefinite Relaxation Procedure for Fault-Tolerant Observer Design. IEEE Trans. Autom. Control. 2015, 60, 3332–3337. [Google Scholar] [CrossRef] [Green Version]
  27. Segundo Sevilla, F.R.; Jaimoukha, I.; Chaudhuri, B.; Korba, P. Fault-tolerant control design to enhance damping of inter-area oscillations in power grids. Int. J. Robust Nonlinear Control. 2014, 24, 1304–1316. [Google Scholar] [CrossRef] [Green Version]
  28. Wang, S.; Meng, X.; Chen, T. Wide-Area Control of Power Systems through Delayed Network Communication. IEEE Trans. Control. Syst. Technol. 2011, 20, 495–503. [Google Scholar] [CrossRef]
  29. Li, Y.; Zhou, Y.; Liu, F.; Cao, Y.; Rehtanz, C. Design and Implementation of Delay-Dependent Wide-Area Damping Control for Stability Enhancement of Power Systems. IEEE Trans. Smart Grid 2016, 8, 1831–1842. [Google Scholar] [CrossRef]
  30. Li, Y.; Liu, F.; Rehtanz, C.; Luo, L.; Cao, Y. Dynamic output-feedback wide area damping control of HVDC transmission considering signal time-varying delay for stability enhancement of interconnected power systems. Renew. Sustain. Energy Rev. 2012, 16, 5747–5759. [Google Scholar] [CrossRef]
  31. Zhang, Y.; Dong, Y.; Hu, S. WAMS-based dynamic output-feedback control of power systems with network-induced delay and random data missing. In Proceedings of the 2015 34th Chinese Control Conference (CCC), Hangzhou, China, 28–30 July 2015. [Google Scholar]
  32. Chaudhuri, B.; Majumder, R.; Pal, B. Wide-Area Measurement-Based Stabilizing Control of Power System Considering Signal Transmission Delay. IEEE Trans. Power Syst. 2004, 19, 1971–1979. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The structure of the augmented closed-loop power system.
Figure 1. The structure of the augmented closed-loop power system.
Energies 14 05892 g001
Figure 2. Power system with four machines in two areas.
Figure 2. Power system with four machines in two areas.
Energies 14 05892 g002
Figure 3. The response of frequency for system with 76–order, 9–order and 8-order.
Figure 3. The response of frequency for system with 76–order, 9–order and 8-order.
Energies 14 05892 g003
Figure 4. (a). The trajectories under different packet   loss   rate   ρ , (b). The trajectories under different time delay with   ρ = 0.1 .
Figure 4. (a). The trajectories under different packet   loss   rate   ρ , (b). The trajectories under different time delay with   ρ = 0.1 .
Energies 14 05892 g004
Figure 5. Response of the system in Scenario 1 when ρ = 0 .
Figure 5. Response of the system in Scenario 1 when ρ = 0 .
Energies 14 05892 g005
Figure 6. Response of the system in Scenario 2 when ρ = 0 .
Figure 6. Response of the system in Scenario 2 when ρ = 0 .
Energies 14 05892 g006
Figure 7. Response of the system in Scenario 3 when ρ = 0 .
Figure 7. Response of the system in Scenario 3 when ρ = 0 .
Energies 14 05892 g007
Figure 8. Response of the system in Scenario 1 when ρ = 0.3 .
Figure 8. Response of the system in Scenario 1 when ρ = 0.3 .
Energies 14 05892 g008
Figure 9. Response of the system in Scenario 2 when ρ = 0.3 .
Figure 9. Response of the system in Scenario 2 when ρ = 0.3 .
Energies 14 05892 g009
Figure 10. Response of the system in Scenario 3 when ρ = 0.3 .
Figure 10. Response of the system in Scenario 3 when ρ = 0.3 .
Energies 14 05892 g010
Figure 11. Response of the system in Scenario 1 when ρ = 0.1 .
Figure 11. Response of the system in Scenario 1 when ρ = 0.1 .
Energies 14 05892 g011
Figure 12. Response of the system in Scenario 1 when ρ = 0.2 .
Figure 12. Response of the system in Scenario 1 when ρ = 0.2 .
Energies 14 05892 g012
Figure 13. Response of the system in Scenario 1 when ρ = 1 .
Figure 13. Response of the system in Scenario 1 when ρ = 1 .
Energies 14 05892 g013
Table 1. The results of modal analysis in a power system without a controller.
Table 1. The results of modal analysis in a power system without a controller.
Mode No.Mode TypeFrequency (Hz)Damping Ratio
1Inter-area0.424−0.0249
2Local1.1620.0952
3Local1.1140.0995
Table 2. The relationship between the rate of packet loss and delay margin (ms).
Table 2. The relationship between the rate of packet loss and delay margin (ms).
ρ 00.010.050.10.150.20.25
d ^ 431390362283264247209
Table 3. The key index of the response curve when there is no packet loss.
Table 3. The key index of the response curve when there is no packet loss.
Figure No.ControllerOvershot (dec\pu)
Figure 5Ours 33.464 \ 9.279 × 10 4
WADC in [20] 34.517 \ 9.859 × 10 4
Conventional 33.375 \ 9.267 × 10 4
Figure 6Ours 25.156 \ 2.205 × 10 4
WADC in [20] 29.696 \ 7.538 × 10 4
Conventional 26.079 \ 8.759 × 10 4
Figure 7Ours 24.900 \ 2.147 × 10 4
WADC in [20] 24.900 \ 4.099 × 10 4
Conventional 24.900 \ 2.169 × 10 4
Table 4. The key index of the response curve in case of packet loss.
Table 4. The key index of the response curve in case of packet loss.
Figure No.ControllerOvershot (dec\pu)
Figure 8Ours 33.437 \ 9.276 × 10 4
WADC in [20] 34.550 \ 9.912 × 10 4
Conventional 33.357 \ 9.263 × 10 4
Figure 9Ours 25.891 \ 1.100 × 10 3
WADC in [20] 28.843 \ 6.625 × 10 4
Conventional 26.400 \ 8.973 × 10 4
Figure 10Ours 24.900 \ 2.136 × 10 4
WADC in [20] 24.900 \ 3.880 × 10 4
Conventional 24.900 \ 2.459 × 10 4
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sun, M.; Guo, Y.; Song, S. The Delay-Dependent DOFC for Damping Inter-Area Low-Frequency Oscillations in an Interconnected Power System Considering Packet Loss of Wide-Area Signals. Energies 2021, 14, 5892. https://doi.org/10.3390/en14185892

AMA Style

Sun M, Guo Y, Song S. The Delay-Dependent DOFC for Damping Inter-Area Low-Frequency Oscillations in an Interconnected Power System Considering Packet Loss of Wide-Area Signals. Energies. 2021; 14(18):5892. https://doi.org/10.3390/en14185892

Chicago/Turabian Style

Sun, Miaoping, Yaosheng Guo, and Seyha Song. 2021. "The Delay-Dependent DOFC for Damping Inter-Area Low-Frequency Oscillations in an Interconnected Power System Considering Packet Loss of Wide-Area Signals" Energies 14, no. 18: 5892. https://doi.org/10.3390/en14185892

APA Style

Sun, M., Guo, Y., & Song, S. (2021). The Delay-Dependent DOFC for Damping Inter-Area Low-Frequency Oscillations in an Interconnected Power System Considering Packet Loss of Wide-Area Signals. Energies, 14(18), 5892. https://doi.org/10.3390/en14185892

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop