Next Article in Journal
A Process Intensification Approach for CO2 Absorption Using Amino Acid Solutions and a Guanidine Compound
Next Article in Special Issue
Single-Machine Frequency Model and Parameter Identification for Inertial Constraints in Unit Commitment
Previous Article in Journal
A Review: Ion Transport of Two-Dimensional Materials in Novel Technologies from Macro to Nanoscopic Perspectives
Previous Article in Special Issue
Load Frequency Control of Multi-Region Interconnected Power Systems with Wind Power and Electric Vehicles Based on Sliding Mode Control
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Enhanced Stability Criteria of Network-Based Load Frequency Control of Power Systems with Time-Varying Delays

1
College of Automation Science and Technology, South China University of Technology, Guangzhou 510640, China
2
Guangzhou Institute of Standardization, Guangzhou 510050, China
3
School of Automation, Guangdong Polytechnic Normal University, Guangzhou 510665, China
4
School of Electrical Engineering, Yancheng Institute of Technology, Yancheng 224051, China
*
Author to whom correspondence should be addressed.
Energies 2021, 14(18), 5820; https://doi.org/10.3390/en14185820
Submission received: 18 July 2021 / Revised: 7 September 2021 / Accepted: 8 September 2021 / Published: 14 September 2021
(This article belongs to the Special Issue Advances in Power System Stability and Control)

Abstract

:
The stability problem for load frequency control (LFC) of power systems with two time-varying communication delays is studied in this paper. The one-area and two-area LFC systems are considered, respectively, which are modeled as corresponding linear systems with additive time-varying delays. An improved stability criterion is proposed via a modified Lyapunov-Krasovskii functional (LKF) approach. Firstly, an augmented LKF consisting of delay-dependent matrices and some single-integral items containing time-varying delay information in two different delay subintervals is constructed, which makes full use of the coupling information between the system states and time-varying delays. Secondly, the novel negative definite inequality equivalent transformation lemma is used to transform the nonlinear inequality to the linear matrix inequality (LMI) equivalently, which can be easily solved by the MATLAB LMI-Toolbox. Finally, some numerical examples are presented to show the improvement of the proposed approach.

1. Introduction

In order to maintain the power grid frequency (an important index of power quality) fixed or within a small allowable range, a load frequency control (LFC) strategy is a common technique equipped in the power systems [1,2,3,4,5,6,7]. With the development and expansion of the power grid, a dedicated independent communication network has been unable to meet the operation of the power grid, although the small transmission delay in the dedicated independent communication network can be ignored [8]. At present, the LFC scheme receives sensor signals and outputs control signals through an open communication network with a mass of data and extensive information exchange. However, random delays and data packets will be introduced into the LFC scheme through the open communication network, which are not negligible and important [9]. These network factors result in the LFC system performance degradation and even instability. Literature [10] pointed out that even if there is a time delay less than 100 ms in the process of information measurement and control output, the transient excitation controller of a generator cannot achieve the control goal. Thus, it is necessary to study the influence of time-varying delays on the performance of the LFC system in an open communication network.
The network controlled LFC system involves the transmitting data between controller and plant. Therefore, there are two main cases of time-varying delays: on the one hand, only the communication time-varying delay from the control center to the governor is considered [11,12,13,14,15,16,17,18], where delay-dependent stability analysis and controller design are investigated by using single delay to model all time delays arising in communication channels. In fact, the time delays arising in the feedback measurement channel and those in the forward control channel may have different properties. It has more useful guidelines with considering the different properties. Thus, on the other hand, not only the communication time-varying delays from the control center to the governor but also from the sensor to the LFC center are considered simultaneously [19,20,21], where general delay-dependent stability analysis is studied by using additive time-varying delays to model two different time delays arising in communication channels. In this way, the stability problems of the LFC system with two transmission delays can be investigated by using the method of stability analysis for general linear time-delayed systems.
The main analytical purpose is obtaining the stability condition and controller design method via Lyapunov stability theory application. The stability criterion based on Lyapunov stability theory is a sufficient condition, and inevitably has certain conservativeness. There are two main reasons for the conservativeness: the construction of LKF and techniques for estimating the upper bound of the derivative of LKF. Thus, there are many methods and techniques given to address these two aspects. For the construction of LKF, a LKF with the delay decomposition method [22,23,24], LKF with multiple integral items [25,26,27,28,29] and LKF with some augmented vectors [30,31] are proposed. On the other aspect, the Jensen inequality [32], B-L inequality [33] and relaxed integral inequality techniques [34,35,36,37,38] are used to estimate the upper bound of the derivative of LKF. In order to reduce the conservativeness of the LKF construction, a lot of coupling information between the system state variables and time delays is introduced into the LKF, which leads to some nonlinear terms in the final results. This makes the solution complex and even unsolvable. Recently, a novel negative definite inequality equivalent transformation lemma was proposed in [39], which improved the degree of freedom for solving the linear matrix inequality (LMI) in the main theorem without conservativeness. Thus, according to the development of stability methods for linear time-delayed systems, there is still space to further reduce the conservativeness of stability criteria for the LFC system.
Inspired by the above analysis, the contributions of this paper can be summarized as follows:
  • As mentioned above, only one-area LFC system with two different time-varying delays is considered. It is general and important to investigate the stability of two-area or multi-area LFC system with two or more time-varying delays. This paper studies one- and two-area LFC system with two time-varying delays.
  • The main improvements of the LKF are summarized as: (a) introducing four delay-dependent non-integral terms to the LKF, such as S i ( t ) , (i = 1, 2, 3, 4); (b) introducing some integral components to the single-integral terms under different time-varying delay subintervals, such as t h 1 t s x ( θ ) d θ , t h 1 t h 1 t x ( θ ) d θ , t h t s x ( θ ) d θ , t h t h t x ( θ ) d θ , and so on. These improvements make the LKF contain more information (the time-varying delays and the coupling information between the state variables and the time-varying delay) than the literature [11,17,18,21], which reduces the conservativeness caused by the LKF construction.
  • To overcome the nonlinear matrix inequality in the stability criterion, the novel negative definite inequality equivalent transformation lemma proposed in [39] is used to transform the nonlinear inequality to the LMI equivalently, which can be easily solved by the MATLAB LMI-Toolbox.
In this paper, the stability problem for LFC power systems with two time-varying communication delays is studied. Both one-area and two-area LFC systems are considered, respectively. An improved stability criterion is proposed via a modified LKF approach. Firstly, an augmented LKF consisting of delay-dependent matrices and some single-integral items containing time-varying delay information in two different delay subintervals is constructed. Secondly, a novel negative definite inequality equivalent transformation lemma is used to transform the nonlinear inequality to the LMI equivalently, which can be easily solved by the MATLAB LMI-Toolbox. Finally, some numerical examples are presented to show the improvement of the proposed approach. Moreover, the stability results can be applied to the LFC optimization to guarantee the stable operation of power system based on an open interconnect network control.
This paper is organized as follows. Section 2 gives the models of LFC schemes; Section 3 provides a stability assessment for the LFC system. Section 4 shows numerical examples. Conclusions are drawn in Section 5.
Notation 1.
Throughout this paper, the notations are standard. R n denotes the n-dimensional Euclidean space; R n × m is the set of all n × m real matrices; For P R n × n , P > 0 (respectively, P < 0 ) mean that P is a positive (respectively, negative) definite matrix. d i a g a 1 , a 2 , , a n denotes an n-order diagonal matrix with diagonal elements a 1 , a 2 , , a n . e i ( i = 1 , , m ) are block entry matrices. For example, e 2 = 0 I 0 0 m 2 . For a real matrix B and two real symmetric matrices A and C of appropriate dimensions, A B * C denotes a real symmetric matrix, where * denotes the entries implied by symmetry. S y m { A } = A + A T .

2. System Description and Problem Preliminaries

2.1. One-Area LFC System

In this subsection, the model of one-area power system equipped with PI controllers and taking into account the time-varying communication delays is given. The basic diagram of the simplified LFC of one-area power system is shown in Figure 1, where e s d 1 and e s d 2 are time delays, respectively, arising during the measured signal Δ f transmitted from sensor to the LCF center and the control signal sent from the control center to the governor.
According to the LFC system shown in [19] and Figure 1, the common LFC scheme model of one-area can be expressed as follows:
x ¯ ˙ ( t ) = A ¯ x ¯ ( t ) + B ¯ Δ P c ( t ) , y ¯ ( t ) = C ¯ x ¯ ( t ) ,
where
x ¯ ( t ) = Δ f ( t ) Δ P m ( t ) Δ P v ( t ) , y ¯ ( t ) = A C E ( t ) , A ¯ = D M 1 M 0 0 1 T c h 1 T c h 1 R T g 0 1 T g , B ¯ = 0 0 1 T g , C ¯ = β 0 0 .
Here, Δ f ( t ) , Δ P m ( t ) and Δ P v ( t ) are the frequency deviation, the mechanical output change, and the valve position change, respectively; M and D are the moment of inertia of the generator and generator damping coefficient, respectively; T g and T c h are the time constant of the governor and the turbine, respectively; R is the speed droop; Δ P C ( t ) is the setpoint; and β is the frequency bias factor. The following PI controller is used as the LFC scheme:
u ( t ) = K P A C E ( t ) K I A C E ( t ) d t ,
where K P and K I are PI gains; and the A C E ( t ) is the area control error. Due to the existence of time-varying delays ( h 1 ( t ) and h 2 ( t ) ) in feedback and forward channels, respectively, the following is obtained
Δ P C ( t ) = u ( t h 2 ( t ) ) , A C E ( t ) = β Δ f ( t h 1 ( t ) ) .
By defining virtual state and measurement output vector as y ( t ) = col { A C E ( t ) , A C E ( t ) d t } and x ( t ) = col { Δ f ( t ) , Δ P m ( t ) , Δ P v ( t ) , A C E ( t ) d t } , the closed-loop LFC system can be expressed as the following linear system with two additive time-varying delays:
x ˙ ( t ) = A x ( t ) + A 1 x ( t h 1 ( t ) h 2 ( t ) )
and the system parameters are listed in the following form
A = D M 1 M 0 0 0 1 T c h 1 T c h 0 1 R T g 0 1 T g 0 β 0 0 0 , A 1 = 0 0 0 0 0 0 0 0 K P β T g 0 0 K I T g 0 0 0 0 .

2.2. Two-Area LFC System

In this subsection, models of two-area power system equipped with PI controllers and taking into account the time-varying communication delays are given. The basic diagram of the simplified LFC of two-area power system is shown in Figure 2, where Δ P 12 and T 12 are the tie-line power transfer and synchronising coefficient of the tie-line. According to the LFC system shown in [20] and Figure 2, the ACE can be expressed as follows:
A C E 1 ( t ) = Δ P 12 ( t ) + β 1 Δ f 1 ( t ) ,
A C E 2 ( t ) = Δ P 12 ( t ) + β 2 Δ f 2 ( t )
and the closed-loop LFC system with PI controllers and time-varying communication delays can be described as follows:
x ^ ˙ ( t ) = A ^ x ^ ( t ) + A ^ 1 x ^ ( t h 1 ( t ) h 2 ( t ) ) ,
where
x ^ ( t ) = col Δ f 1 ( t ) , Δ P m 1 ( t ) , Δ P v 1 ( t ) , A C E 1 ( t ) d t , Δ P 12 ( t ) , Δ f 2 ( t ) , Δ P m 2 ( t ) , Δ P v 2 ( t ) , A C E 2 ( t ) d t ,
A ^ = A 11 A 12 A 21 A 22 , A ^ 1 = A 111 A 112 A 121 A 122 ,
A 11 = D 1 M 1 1 M 1 0 0 1 M 1 0 1 T c h 1 1 T c h 1 0 0 1 R 1 T g 1 0 1 T g 1 0 0 β 1 0 0 0 1 2 π T 12 0 0 0 0 , A 12 = 0 0 0 0 0 0 0 0 0 0 0 0 2 π T 12 0 0 0 ,
A 21 = 0 0 0 0 1 M 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 , A 22 = D 2 M 2 1 M 2 0 0 0 1 T c h 2 1 T c h 2 0 1 R 2 T g 2 0 1 T g 2 0 β 2 0 0 0 , A 111 = 0 0 0 0 0 0 0 0 0 0 K P 1 β 1 T g 1 0 0 K I 1 T g 1 K P 1 T g 1 0 0 0 0 0 0 0 0 0 0 , A 112 = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 , A 121 = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 K P 2 T g 2 0 0 0 0 0 , A 122 = 0 0 0 0 0 0 0 0 K P 2 β 2 T g 2 0 0 K I 2 T g 2 0 0 0 0 .
As discussed in the Introduction section, the communication channel may encounter time delays, which are usually time-varying and bounded. Then, similar to the previous work, the delays are expressed by time-varying functions satisfying the following conditions:
0 h 1 ( t ) h 1 , 0 h 2 ( t ) h 2 , h ˙ 1 ( t ) μ 1 , h ˙ 2 ( t ) μ 2 , t 0 ,
where h 1 , h 2 , μ 1 and μ 2 are positive constants. Let h ( t ) = h 1 ( t ) + h 2 ( t ) and h = h 1 + h 2 .

3. Stability Assessment for LFC System

Lemma 1
([33]). For a positive definite matrix R and differentiable function x in [ a , b ] R n , the following holds
a b x ˙ T ( s ) R x ˙ ( s ) d s 1 b a ω T R ω ,
where R = d i a g { R , 3 R , 5 R } , ω = c o l { ω 1 , ω 2 , ω 3 } with ω 1 = x ( b ) x ( a ) , ω 2 = x ( b ) + x ( a ) 2 b a a b x ( s ) d s , ω 3 = ω 1 6 b a a b x ( s ) d s + 12 ( b a ) 2 a b ( b s ) x ( s ) d s .
Lemma 2
([30]). For positive definite matrices R 1 , R 2 R n × n , vectors v 1 , v 2 R n and a scalar α [ 0 , 1 ] , if there exists symmetric matrices X 1 , X 2 R n × n and any matrices S 1 , S 2 R n × n such that
R 1 X 1 S 1 * R 1 0 , R 2 X 2 S 2 * R 2 0 ,
the following inequality holds
1 α v 1 T R 1 v 1 + 1 1 α v 2 T R 2 v 2 v 1 T [ R 1 + ( 1 α ) X 1 ] v 1 + v 2 T [ R 2 + α X 2 ] v 2 + 2 v 1 T [ α S 1 + ( 1 α ) S 2 ] v 2 .
Lemma 3
([39]). Let symmetric matrices A 0 , A 1 , A 2 R m × m and a vector ζ R m . Then, the following inequality
ζ T ( h t 2 A 2 + h t A 1 + A 0 ) ζ < 0
holds for all h t [ 0 , h ] if and only if there exists a positive definite matrix D R m × m and a skew-symmetric matrix G R k × k such that
A 0 1 2 A 1 * A 2 < C J T D G * D C J ,
where C = h 2 I 0 and J = h 2 I I .
Theorem 1.
Given scalars h 1 , h 2 , μ 1 and μ 2 , the system (7) satisfying condition (8) is asymptotically stable, if there exists positive definite matrices S i 2 , Q j R 3 n × 3 n , Q k R 8 n × 8 n , R i R n × n , D i , G i R 18 n × 18 n , symmetric matrices S i 1 , X i R 3 n × 3 n , any matrices Y i R 3 n × 3 n , ( i = 1 , 2 , 3 , 4 ; j = 1 , 2 ; k = 3 , 4 ), such that the following matrix inequalities hold for h ˙ 1 t μ j { μ 1 , μ 1 } , h ˙ t μ k { μ , μ } .
h 1 S j 1 + S j 2 > 0 , h S k 1 + S k 2 > 0 ,
R ¯ i X i Y i * R ¯ i > 0 ,
Ω 10 ( μ j ) 1 2 Ω 11 ( μ j ) * Ω 21 ( μ j ) C J T D j G j * D j C J < 0 ,
Ω 0 ( μ k ) 1 2 Ω 1 ( μ k ) * Ω 2 ( μ k ) C J T D k G k * D k C J < 0
with C = h 2 I 0 , J = h 2 I I , R ¯ i = diag { R i , 3 R i , 5 R i } ,
Ω 10 ( μ j ) = h ˙ 1 t Δ 11 T S 11 Δ 11 h ˙ 1 t Δ 21 T S 21 Δ 21 h 1 d Δ 53 T Q 1 Δ 53 + h 1 d Δ 63 T Q 2 Δ 63 + h 1 2 h 1 d e 7 T ( R 1 R 2 ) e 7 , Ω 11 ( μ j ) = Sym Δ 13 T S 11 Δ 12 Δ 23 T S 21 Δ 22 + Δ 92 T Q 1 Λ 1 + Λ 12 T Q 2 Λ 2 + h ˙ 1 t Δ 12 T S 11 Δ 12 h ˙ 1 t Δ 22 T S 21 Δ 22 + Δ 52 T Q 1 Δ 52 Δ 62 T Q 2 Δ 62 + Δ 72 T Q 3 Δ 72 Δ 82 T Q 4 Δ 82 h 1 d Δ 54 T Q 1 Δ 54 + h 1 d Δ 64 T Q 2 Δ 64 ,
Ω 21 ( μ j ) = Sym Δ 13 T ( S 11 Δ 11 + S 12 Δ 12 ) + h ˙ 1 t Δ 11 T S 11 Δ 12 + Δ 23 T [ ( h 1 S 21 + S 22 ) Δ 22 S 21 Δ 21 ] h ˙ 1 t Δ 22 T S 21 Δ 21 + Δ 51 Q 1 Δ 52 Δ 61 Q 1 Δ 62 h 1 d Δ 53 Q 1 Δ 54 + h 1 d Δ 63 Q 2 Δ 64 + Δ 91 T Q 1 Λ 1 + Λ 11 T Q 2 Λ 2 h 1 h 1 d e 7 T ( R 1 R 2 ) e 7 1 h 1 Γ 1 T X 1 Γ 1 Γ 2 T X 2 Γ 2 + Sym 1 h 1 Γ 1 T ( Y 1 Y 2 ) Γ 2 ,
Ω 0 ( μ k ) = Sym Δ 13 T S 12 Δ 11 + Δ 23 T ( h 1 S 21 + S 22 ) Δ 21 + Δ 33 T S 32 Δ 31 + Δ 43 T ( h S 41 + S 42 ) Δ 41 + Δ 90 T Q 1 Λ 1 + Λ 10 T Q 2 Λ 2 + Λ 20 T Q 3 Λ 3 + Λ 30 T Q 4 Λ 4 + h ˙ t Δ 31 T S 31 Δ 31 h ˙ t Δ 41 T S 41 Δ 41 + Δ 51 T Q 1 Δ 51 Δ 61 T Q 2 Δ 61 + Δ 71 T Q 3 Δ 71 Δ 81 T Q 4 Δ 81 h d Δ 73 T Q 3 Δ 73 + h d Δ 83 T Q 4 Δ 83 + e 6 T ( h 1 2 R 2 + h 2 R 4 ) e 6 + h 2 h d e 7 T ( R 3 R 4 ) e 7 + Γ 1 T R ¯ 1 Γ 1 + Γ 2 T ( R ¯ 2 + X 2 ) Γ 2 + Γ 3 T R ¯ 3 Γ 3 + Γ 4 T ( R ¯ 4 + X 4 ) Γ 4 + Sym { Γ 1 T Y 1 Γ 2 + Γ 3 T Y 3 Γ 4 } + Sym { Δ 0 T U Π 0 } ,
Ω 1 ( μ k ) = Sym Δ 33 T S 31 Δ 31 + Δ 33 T S 32 Δ 32 + h ˙ t Δ 31 T S 31 Δ 32 + Δ 43 T ( h S 41 + S 42 ) Δ 42 Δ 43 T S 41 Δ 41 h ˙ t Δ 42 T S 41 Δ 41 + Λ 22 T Q 3 Λ 3 + Λ 32 T Q 4 Λ 4 + Δ 72 T Q 3 Δ 72 Δ 82 T Q 4 Δ 82 h ˙ t Δ 74 T Q 3 Δ 74 + h d Δ 84 T Q 4 Δ 84 h h d e 8 T ( R 3 R 4 ) e 8 1 h [ Γ 3 T X 3 Γ 3 Γ 4 T X 4 Γ 4 ] + 1 h Sym { Γ 3 T ( Y 3 Y 4 ) Γ 4 } ,
Ω 2 ( μ k ) = Sym Δ 33 T S 31 Δ 32 Δ 43 T S 41 Δ 42 + Λ 22 T Q 3 Λ 3 + Λ 33 T Q 4 Λ 4 + h ˙ t Δ 32 T S 31 Δ 32 h ˙ t Δ 42 T S 41 Δ 42 + Δ 72 T Q 3 Δ 72 Δ 82 T Q 4 Δ 82 h d Δ 74 T Q 3 Δ 74 + h d Δ 84 T Q 4 Δ 84 ,
where the notations of other symbols and matrices can be found in Appendix A and Appendix B.
Proof. 
Refer to Appendix C. □
Remark 1.
The stability sufficient condition of the LFC system (4) is obtained in Theorem 1, which can guarantee the stability of the LFC system (4) in the range of the maximum allowable time delay. Compared with the literature [11,17,18,21], the Theorem 1 in this paper further reduces the conservativeness via the augmented LKF application. The main improvements of the LKF are summarized as: (a) introducing four delay-dependent non-integral terms to the LKF, such as S i ( t ) , ( i = 1 , 2 , 3 , 4 ); (b) introducing some integral components to the single-integral terms under different time-varying delay subintervals, such as t h 1 t s x ( θ ) d θ , t h 1 t h 1 t x ( θ ) d θ , t h t s x ( θ ) d θ , t h t h t x ( θ ) d θ , and so on. These improvements make the LKF contain more information of the time-varying delay and the coupling information between the state variables and the delay than the literature [11,17,18,21], which reduces the conservativeness caused by the LKF construction.
Remark 2.
It is well known that in order to improve the conservativeness of the stability criterion, the modified LKF must combine with some tight inequality technique [36]. Therefore, Lemmas 1 and 2 are used to estimate the derivative of the LKF, where some effective auxiliary functions and free-weight matrices are introduced into the upper bound on the derivative of the LKF. Finally, Theorem 1 based on the modified LKF and two proper inequality techniques is less conservative than those of the recently published literature [11,17,18,21].
Remark 3.
However, according to the inequality (A8) in the proof process, the final form of the upper bound on the derivative of the LKF is nonlinear due to h 2 ( t ) Ω 2 ( h ˙ t ) , which cannot be solved by MATLAB. The authors of [25] decomposed this nonlinear matrix inequality into three equivalent LMI conditions via the sufficient condition constraint lemma application. It reduced the degree of freedom for solving the matrix inequality in the main theorem, due to the introduction of two additional LMI constraints. Moreover, the sufficient condition constraint lemma is quite conservative [39]. To overcome the nonlinear matrix inequality in the stability criterion, the novel negative definite inequality equivalent transformation lemma (Lemma 3) proposed in [39] is used to transform the inequality (A8) to the LMI (12) equivalently, which can be easily solved by the MATLAB LMI-Toolbox.

4. Results and Discussions

In this section, the effectiveness of the stability criterion proposed in this paper is shown. For different K P and K I values, the maximum allowable time-delay upper bound values (MADUB) can be obtained by solving the LMIs in Theorem 1 via Matlab LMI-Toolbox. The LFC system parameters in Table 1 are given in [20]. One and two-area LFC systems will be discussed and comparatively analyzed in the following subsections. In the following subsections, ‘–’ in the tables indicates that the corresponding result is not given.

4.1. One-Area LFC System

4.1.1. Conservativeness Comparisons

In order to compare with the existing results, Table 2 and Table 3 give the MADUB values of the case of fixed K P and K I values, K P = 0.2 , K I = 0.1 , | h ˙ 1 ( t ) | 0.1 , | h ˙ 2 ( t ) | 0.8 . The corresponding results cannot be given in [11,17,18] since the time delays from sensor to the LFC center are ignored. Moreover, to show the PI controller gains influence on the MADUB of the LFC system, Table 4 obtains the MADUB values of the case of different K P , K I values, h 1 ( t ) = 0 and μ 2 = 0 ; Table 5 gives the MADUB values of the case of different K P , K I values, h 1 ( t ) = 0 and μ 2 = 0.9 . From these tables, it can observe the results of Theorem 1 are similar to those of [11], however, less conservative than those of [14,15,16,19,20,21]. Meanwhile, the MADUB increases with K P increment at a fixed K I value, and the MADUB decreases with K I increment at a fixed K P value.

4.1.2. Simulation Verification

Simulation studies are carried out under an increased step load of 0.1 pu occurring at 1 s, time-varying delays and the following assumptions. The simulation results are shown in Figure 3, Figure 4, Figure 5 and Figure 6, in which the LFC has achieved its objective, and the control system is stable. In Figure 4, the blue curve is close to the critical stability with K P = 0 , K I = 0.05 , h 1 ( t ) = 0 and h 2 ( t ) = 30.79 .
  • For Figure 3, fixed K P = 0.2 , K I = 0.1 :
    • h 1 ( t ) = 1 2 s i n 0.2 1 t + 1 2 , h 2 ( t ) = 6.155 2 s i n 1.6 6.155 t + 6.155 2 ;
    • h 1 ( t ) = 1.5 2 s i n 0.2 3 t + 1.5 2 , h 2 ( t ) = 6.003 2 s i n 1.6 6.003 t + 6.003 2 ;
    • h 1 ( t ) = 5.016 2 s i n 0.2 5.016 t + 5.016 2 , h 2 ( t ) = s i n 0.8 t + 1 ;
    • h 1 ( t ) = 3.217 2 s i n 0.2 3.217 t + 3.217 2 , h 2 ( t ) = 2 s i n 0.4 t + 2 ;
  • For Figure 4, different K P and fixed h 1 ( t ) = 0 , K I = 0.05 :
    • K P = 0 , h 2 ( t ) = 30.7 ;
    • K P = 0 , h 2 ( t ) = 30.79 ;
    • K P = 0.1 , h 2 ( t ) = 32.73 ;
  • For Figure 5, different K P and fixed h 1 ( t ) = 0 , K I = 0.6 :
    • K P = 0 , K I = 0.6 , h 2 ( t ) = 2.00 ;
    • K P = 0 , K I = 0.6 , h 2 ( t ) = 2.03 ;
    • K P = 0.1 , K I = 0.6 , h 2 ( t ) = 2.20 ;
  • For Figure 6, different K P and fixed h 1 ( t ) = 0 :
    • K P = 0 , K P = 0.2 , h 2 ( t ) = 7.12 2 s i n 1.8 7.12 t + 7.12 2 ;
    • K P = 0.1 , K P = 0.2 , h 2 ( t ) = 7.13 2 s i n 1.8 7.13 t + 7.13 2 ;
    • K P = 0 , K P = 1 , h 2 ( t ) = 0.79 2 s i n 1.8 0.79 t + 0.79 2 ;
    • K P = 0.1 , K P = 1 , h 2 ( t ) = 0.87 2 s i n 1.8 0.87 t + 0.87 2 .

4.2. Two-Area LFC System

4.2.1. Conservativeness Comparison

Table 6 and Table 7 give the MADUB values of the case of fixed K P and K I values, K P = 0.05 , K I = 0.4 and different μ 1 , μ 2 values. The corresponding results cannot be achieved in [11,17,18] since the time delays from sensor to the load frequency control center are ignored. Table 8 gives the MADUB values of the case of different K P and K I values, h 1 ( t ) = 0 and μ 2 = 0 ; and Table 9 give the MADUB values of the case of different K P and K I values, h 1 ( t ) = 0 and μ 2 = 0.5 . From these tables, it is obviously that the MADUB increases with K P increment at a fixed K I value, and the MADUB decreases with K I increment at a fixed K P value. In addition, the results of Theorem 1 are less conservative than those of [11,15,20].

4.2.2. Simulation Verification

Simulation studies are carried out under an increased step load of 0.1 pu occurring at 1 s, and time-varying delays and the following assumptions:
  • For Figure 7, fixed K P = 0.05 , K I = 0.4 :
    • h 1 ( t ) = 1 2 s i n 0.4 1 t + 1 2 , h 2 ( t ) = 2.018 2 s i n 0.4 2.018 t + 2.018 2 ;
    • h 1 ( t ) = 2 2 s i n 1 2 t + 2 2 , h 2 ( t ) = 0.908 2 s i n 1.6 0.908 t + 0.908 2 ;
    • h 1 ( t ) = 2.028 2 s i n 0.2 2.028 t + 2.028 2 , h 2 ( t ) = 1 2 s i n 0.4 1 t + 1 2 ;
    • h 1 ( t ) = 0.036 2 s i n 1 0.036 t + 0.036 2 , h 2 ( t ) = 3 2 s i n 1.6 3 t + 3 2 ;
  • For Figure 8, different K P and fixed h 1 ( t ) = 0 , K I = 0.05 :
    • K P = 0 , h 2 ( t ) = 30.89 ;
    • K P = 0 , h 2 ( t ) = 30.83 ;
    • K P = 0.1 , h 2 ( t ) = 31.67 ;
  • For Figure 9, different K P and fixed h 1 ( t ) = 0 , K I = 1 :
    • K P = 0 , h 2 ( t ) = 0.637 ;
    • K P = 0 , h 2 ( t ) = 0.62 ;
    • K P = 0.1 , h 2 ( t ) = 0.65 ;
  • For Figure 10, different K P , K I , time-varying h 2 ( t ) and fixed h 1 ( t ) = 0 :
    • K P = 0 , K I = 0.2 , h 2 ( t ) = 6.69 2 s i n 1 6.69 t + 6.69 2 ;
    • K P = 0.1 , K I = 0.2 , h 2 ( t ) = 7.13 2 s i n 1 7.13 t + 7.13 2 ;
    • K P = 0 , K I = 0.6 , h 2 ( t ) = 1.65 2 s i n 1 1.65 t + 1.65 2 ;
    • K P = 0.1 , K I = 0.6 , h 2 ( t ) = 1.78 2 s i n 1 1.78 t + 1.78 2 .
The simulation results are shown in Figure 7, Figure 8, Figure 9 and Figure 10, in which the LFC has achieve its objective, and the control system is stable. It can be seen from Figure 8 and Figure 9 that the red curves are the critical stability with K I = 0.05 , K P = 0 , h 1 ( t ) = 0 , h 2 ( t ) = 30.89 and K P = 0 , K I = 1 , h 1 ( t ) = 0 , h 2 ( t ) = 0.637 , respectively. Thus, the black curvess in Figure 8 and Figure 9 are close to the critical stability region, that is, Theorem 1 proposed in this paper is effective in estimating the upper bound of the maximum allowable time delay.

5. Conclusions

This paper mainly focuses on the stability analysis for load frequency control of power systems with time-varying delays. For the one-area and two-area LFC systems with two communication delays, stability criteria are obtained via Lyapunov stability theory application. Firstly, the one-area LFC system and two-area LFC system are described as linear systems with additive time-varying delays. Secondly, a modified LKF with some delay-dependent non-integral terms and augmented integral components in single integral term is constructed. Compared with the LKFs in some previous published literature, it contains more coupling information between time-varying delays and state variables, which reduces the conservativeness of the stability criterion. Thirdly, to overcome the nonlinear coupling in the stability criterion, the novel negative definite inequality equivalent transformation lemma (Lemma 3) is used to transform the nonlinear inequality to the LMI equivalently, which can be easily solved by the MATLAB LMI-Toolbox. Finally, the effectiveness of the proposed method is illustrated by comparisons and discussions in numerical examples. In addition, in order to approach the actual situation, the design and optimization of the controllers have always been a topic of concern, which inspires us to work around it in the future. The stability results can be applied to the LFC design and optimization to guarantee the stable operation of power system based on an open interconnect network control.
However, there are some limitations. Firstly, in order to reduce the conservativeness of the stability criterion, the dimensions of the decision variables for solving LMIs increase with the expanded dimensions of LKF. Secondly, in order to improve the degree of freedom of solving LMIs, the free weight matrices are introduced. Thirdly, some additional decision variables are introduced in LMIs when the nonlinear inequalities are transformed into LMIs by using Lemma 3, which increases the number of decision variables of the LMIs proposed in this paper. In conclusion, the improvement of stability results is in the cost of increasing computational complexity.
The derivation method of the stability criterion presented in this paper can be extended to multi-area LFC system. The relevant theories will be applied into practice, which is one of our further main topics.

Author Contributions

Conceptualization, W.F. and Y.X.; data curation, F.L.; methodology, W.F. and W.D.; software, X.Z.; writing—original draft, W.F.; writing—review and editing, F.L. and X.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported partly by the Natural Science Foundation of Guangdong Province under Grant 2018A0303130111, the Industry University Research Cooperation Project of Jiangsu Province under Grant BY2020649 and the Yellow Sea Rookie of Yancheng Institute of Technology.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

For the sake of simplicity on matrix representation, the notations of several symbols and matrices are defined as
h 1 t = h 1 ( t ) , h t = h ( t ) , h ¯ 1 t = h 1 h 1 ( t ) , h ¯ t = h h ( t ) , h 1 d = 1 h ˙ 1 ( t ) , h d = 1 h ˙ ( t ) , η 01 ( t ) = col x ( t ) , x ( t h 1 t ) , x ( t h 1 ) , η 02 ( t ) = col x ( t ) , x ( t h t ) , x ( t h ) , η 1 ( t ) = col x ( t ) , x ( t h 1 t ) , t h 1 t t x ( s ) d s , η 2 ( t ) = col x ( t h 1 t ) , x ( t h 1 ) , t h 1 t h 1 t x ( s ) d s , η 3 ( t ) = col x ( t ) , x ( t h t ) , t h t t x ( s ) d s , η 4 ( t ) = col x ( t h t ) , x ( t h ) , t h t h t x ( s ) d s , η 5 ( t , s ) = col x ˙ ( s ) , x ( s ) , η 01 ( t ) , s t x ( θ ) d θ , t h 1 t s x ( θ ) d θ , t h 1 t h 1 t x ( θ ) d θ , η 6 ( t , s ) = col x ˙ ( s ) , x ( s ) , η 01 ( t ) , s t h 1 t x ( θ ) d θ , t h 1 t t x ( θ ) d θ , t h 1 s x ( θ ) d θ , η 7 ( t , s ) = col x ˙ ( s ) , x ( s ) , η 02 ( t ) , s t h 1 t x ( θ ) d θ , t h 1 t t x ( θ ) d θ , t h 1 s x ( θ ) d θ , η 8 ( t , s ) = col x ˙ ( s ) , x ( s ) , η 02 ( t ) , s t h t x ( θ ) d θ , t h t t x ( θ ) d θ , t h s x ( θ ) d θ , ρ 1 ( t ) = t h 1 t h 1 t x ( s ) h ¯ 1 t , ρ 2 ( t ) = t h 1 t h 1 t ( t h 1 t s ) x ( s ) h ¯ 1 t 2 , ρ 3 ( t ) = t h t h t x ( s ) h ¯ t , ρ 4 ( t ) = t h t h t ( t h t s ) x ( s ) h ¯ t 2 , ρ 5 ( t ) = t h 1 t t x ( s ) h 1 t , ρ 6 ( t ) = t h 1 t t ( t s ) x ( s ) h 1 t 2 , ρ 7 ( t ) = t h t t x ( s ) h t , ρ 8 ( t ) = t h t t ( t s ) x ( s ) h t 2 ,
ξ ( t ) = col x ( t ) , x ( t h 1 t ) , x ( t h t ) , x ( t h 1 ) , x ( t h ) , x ˙ ( t ) , x ˙ ( t h 1 t ) , x ˙ ( t h t ) , x ˙ ( t h 1 ) , x ˙ ( t h ) , ρ 1 ( t ) , ρ 2 ( t ) , ρ 3 ( t ) , ρ 4 ( t ) , ρ 5 ( t ) , ρ 6 ( t ) , ρ 7 ( t ) , ρ 8 ( t ) .

Appendix B

Notations of other symbols and matrices for the Theorem 1 are given as:
Δ 0 = col { e 1 , e 3 , e 6 } , Π 0 = A e 1 + A 1 e 3 e 6 , Δ 11 = col { e 1 , e 2 , e 0 } , Δ 12 = col { e 0 , e 0 , e 13 } , Δ 13 = col { e 6 , e 7 , e 1 h 1 d e 2 } , Δ 21 = col { e 2 , e 4 , h 1 e 11 } , Δ 22 = col { e 0 , e 0 , e 11 } , Δ 23 = col { h 1 d e 7 , e 8 , h 1 d e 2 e 4 } , Δ 31 = col { e 1 , e 3 , h 1 e 0 } , Δ 32 = col { e 0 , e 0 , e 17 } , Δ 33 = col { e 6 , h d e 3 , e 1 h d e 3 } , Δ 41 = col { e 3 , e 5 , h e 13 } , Δ 42 = col { e 0 , e 0 , e 13 } , Δ 43 = col { h d e 8 , e 10 , h d e 3 e 5 } , Δ 51 = col { e 6 , e 1 , e 1 , e 2 , e 4 , e 0 , e 0 , h 1 e 11 } , Δ 52 = col { e 0 , e 0 , e 0 , e 0 , e 0 , e 0 , e 15 , e 11 } , Δ 53 = col { e 7 , e 2 , e 1 , e 2 , e 4 , e 0 , e 0 , h 1 e 11 } , Δ 54 = col { e 0 , e 0 , e 0 , e 0 , e 0 , e 15 , e 0 , e 11 } , Δ 61 = col { e 9 , e 4 , e 1 , e 2 , e 4 , h 1 e 11 , e 0 , e 0 } , Δ 62 = col { e 0 , e 0 , e 0 , e 0 , e 0 , e 11 , e 15 , e 0 } , Δ 63 = col { e 7 , e 2 , e 1 , e 2 , e 4 , e 0 , e 0 , h 1 e 11 } , Δ 64 = col { e 0 , e 0 , e 0 , e 0 , e 0 , e 0 , e 15 , e 11 } , Δ 71 = col { e 6 , e 1 , e 1 , e 3 , e 5 , e 0 , e 0 , h e 13 } , Δ 72 = col { e 0 , e 0 , e 0 , e 0 , e 0 , e 0 , e 17 , e 13 } , Δ 73 = col { e 8 , e 3 , e 1 , e 3 , e 5 , e 0 , e 0 , h e 13 } , Δ 74 = col { e 0 , e 0 , e 0 , e 0 , e 0 , e 17 , e 0 , e 13 } , Δ 81 = col { e 10 , e 5 , e 1 , e 3 , e 5 , h e 13 , e 0 , e 0 } , Δ 82 = col { e 0 , e 0 , e 0 , e 0 , e 0 , e 13 , e 17 , e 0 } , Δ 83 = col { e 8 , e 3 , e 1 , e 3 , e 5 , e 0 , e 0 , h e 13 } , Δ 84 = col { e 0 , e 0 , e 0 , e 0 , e 0 , e 0 , e 17 , e 13 } , Δ 90 = col { e 1 e 2 , e 0 , e 0 , e 0 , e 0 , e 0 , e 0 , e 0 } , Δ 91 = col { e 0 , e 15 , e 1 , e 2 , e 4 , e 0 , e 0 , h 1 e 11 } , Δ 92 = col { e 0 , e 0 , e 0 , e 0 , e 0 , e 15 e 16 , e 16 , e 11 } , Λ 1 = col { e 0 , e 0 , e 6 , h 1 d e 7 , e 1 h 1 d e 7 , e 1 , h 1 d e 2 , h 1 d e 2 e 4 } , Λ 2 = col { e 0 , e 0 , e 6 , h 1 d e 7 , e 9 , h 1 d e 2 , e 1 h 1 d e 2 , e 4 } , Λ 3 = col { e 0 , e 0 , e 6 , h d e 8 , e 10 , e 1 , h d e 3 , h d e 3 e 5 } , Λ 4 = col { e 0 , e 0 , e 6 , h d e 8 , e 10 , h d e 3 , e 1 h d e 3 , e 5 } , Λ 10 = col { e 2 e 4 , h 1 e 1 , h 1 e 1 , h 1 e 2 , h 1 e 4 , h 1 2 ( e 11 e 12 ) , e 0 , h 1 2 e 12 } , Λ 11 = col { e 0 , e 11 , e 1 , e 2 , e 4 , 2 h 1 ( e 11 e 12 ) , h 1 e 15 , 2 h 1 e 12 } , Λ 12 = col { e 0 , e 0 , e 0 , e 0 , e 0 , e 11 e 12 , e 15 , e 12 } , Γ 1 = col { e 2 e 4 , e 2 + e 4 2 e 11 , e 2 e 4 + 6 e 11 + 12 e 12 } , Γ 2 = col { e 1 e 2 , e 1 + e 2 2 e 15 , e 1 e 2 + 6 e 15 + 12 e 16 } , Γ 3 = col { e 3 e 5 , e 3 + e 5 2 e 13 , e 3 e 5 + 6 e 13 + 12 e 14 } , Γ 4 = col { e 1 e 3 , e 1 + e 3 2 e 17 , e 1 e 3 + 6 e 17 + 12 e 18 } .

Appendix C

Proof of Theorem 1. 
Construct an LKF candidate as
V ( t ) = V 1 ( t ) + V 2 ( t ) + V 3 ( t )
with
V 1 ( t ) = η 1 T ( t ) S 1 ( t ) η 1 ( t ) + η 2 T ( t ) S 2 ( t ) η 2 ( t ) + η 3 T ( t ) S 3 ( t ) η 3 ( t ) + η 4 T ( t ) S 4 ( t ) η 4 ( t ) , V 2 ( t ) = t h 1 t t η 5 T ( t , s ) Q 1 η 5 ( t , s ) d s + t h 1 t h 1 t η 6 T ( t , s ) Q 2 η 6 ( t , s ) d s + t h t t η 7 T ( t , s ) Q 3 η 7 ( t , s ) d s + t h t h t η 8 T ( t , s ) Q 4 η 8 ( t , s ) d s , V 3 ( t ) = h 1 t h 1 t h 1 t ( h 1 t t + s ) x ˙ T ( s ) R 1 x ˙ ( s ) d s + h 1 t h 1 t t ( h 1 t + s ) x ˙ T ( s ) R 2 x ˙ ( s ) d s + h t h t h t ( h t t + s ) x ˙ T ( s ) R 3 x ˙ ( s ) d s + h t h t t ( h t + s ) x ˙ T ( s ) R 4 x ˙ ( s ) d s ,
where S 1 ( t ) = h 1 t S 11 + S 12 , S 2 ( t ) = h ¯ 1 t S 21 + S 22 , S 3 ( t ) = h t S 31 + S 32 and S 4 ( t ) = h ¯ t S 41 + S 42 .
Calculating the derivative of V ( t ) , we can obtain the following formulas
V ˙ 1 ( t ) = 2 η ˙ 1 T ( t ) S 1 ( t ) η 1 ( t ) + η 1 T ( t ) S ˙ 1 ( t ) η 1 ( t ) + 2 η ˙ 2 T ( t ) S 2 ( t ) η 2 ( t ) + η 2 T ( t ) S ˙ 2 ( t ) η 2 ( t ) + 2 η ˙ 3 T ( t ) S 3 ( t ) η 3 ( t ) + η 3 T ( t ) S ˙ 3 ( t ) η 3 ( t ) + 2 η ˙ 4 T ( t ) S 4 ( t ) η 4 ( t ) + η 4 T ( t ) S ˙ 4 ( t ) η 4 ( t ) = 2 ξ T ( t ) Δ 13 T ( h 1 t S 11 + S 12 ) ( Δ 11 + h 1 t Δ 12 ) ξ ( t ) + ξ T ( t ) ( Δ 11 + h 1 t Δ 12 ) T h ˙ 1 t S 11 ( Δ 11 + h 1 t Δ 12 ) ξ ( t ) + 2 ξ T ( t ) Δ 23 T ( h ¯ 1 t S 21 + S 22 ) ( Δ 21 + h 1 t Δ 22 ) ξ ( t ) ξ T ( t ) ( Δ 21 + h 1 t Δ 22 ) T h ˙ 1 t S 21 ( Δ 21 + h 1 t Δ 22 ) ξ ( t ) + 2 ξ T ( t ) Δ 33 T ( h t S 31 + S 32 ) ( Δ 31 + h t Δ 32 ) ξ ( t ) + ξ T ( t ) ( Δ 31 + h t Δ 32 ) T h ˙ t S 31 ( Δ 31 + h t Δ 32 ) ξ ( t ) + 2 ξ T ( t ) Δ 43 T ( h ¯ t S 41 + S 42 ) ( Δ 41 + h t Δ 42 ) ξ ( t ) ξ T ( t ) ( Δ 41 + h t Δ 42 ) T h ˙ t S 41 ( Δ 41 + h t Δ 42 ) ξ ( t ) ,
V ˙ 2 ( t ) = η 5 T ( t , t ) Q 1 η 5 ( t , t ) η 6 T ( t , t h 1 ) Q 2 η 6 ( t , t h 1 ) + η 7 T ( t , t ) Q 3 η 7 ( t , t ) η 8 T ( t , t h ) Q 4 η 8 ( t , t h ) h 1 d η 5 T ( t , t h 1 t ) Q 1 η 5 ( t , t h 1 t ) + h 1 d η 6 T ( t , t h 1 t ) Q 2 η 6 ( t , t h 1 t ) h d η 7 T ( t , t h t ) Q 3 η 7 ( t , t h t ) + h d η 8 T ( t , t h t ) Q 4 η 8 ( t , t h t ) + 2 t h 1 t t η 5 T ( t , s ) Q 1 t η 5 ( t , s ) d s + 2 t h 1 t h 1 t η 6 T ( t , s ) Q 2 t η 6 ( t , s ) d s + 2 t h t t η 7 T ( t , s ) Q 3 t η 7 ( t , s ) d s + 2 t h t h t η 8 T ( t , s ) Q 4 t η 8 ( t , s ) d s = ξ T ( t ) ( Δ 51 + h 1 t Δ 52 ) T Q 1 ( Δ 51 + h 1 t Δ 52 ) ( Δ 61 + h 1 t Δ 62 ) T Q 2 ( Δ 61 + h 1 t Δ 62 ) ξ ( t ) + ξ T ( t ) ( Δ 71 + h t Δ 72 ) T Q 3 ( Δ 71 + h t Δ 72 ) ( Δ 81 + h t Δ 82 ) T Q 4 ( Δ 81 + h t Δ 82 ) ξ ( t ) h 1 d ξ T ( t ) ( Δ 53 + h 1 t Δ 54 ) T Q 1 ( Δ 53 + h 1 t Δ 54 ) ( Δ 63 + h 1 t Δ 64 ) T Q 2 ( Δ 63 + h 1 t Δ 64 ) ξ ( t ) h d ξ T ( t ) ( Δ 73 + h t Δ 74 ) T Q 3 ( Δ 73 + h 1 t Δ 74 ) ( Δ 83 + h t Δ 84 ) T Q 4 ( Δ 83 + h t Δ 84 ) ξ ( t ) + 2 ξ T ( t ) ( Δ 90 + h 1 t Δ 91 + h 1 t 2 Δ 92 ) T Q 1 Λ 1 + ( Λ 10 + h 1 t Λ 11 + h 1 t 2 Λ 12 ) T Q 2 Λ 2 ξ ( t ) + 2 ξ T ( t ) ( Λ 20 + h t Λ 21 + h t 2 Λ 22 ) T Q 3 Λ 3 + ( Λ 30 + h t Λ 31 + h t 2 Λ 32 ) T Q 4 Λ 4 ξ ( t ) ,
V ˙ 3 ( t ) = h 1 2 x ˙ T ( t ) R 2 x ˙ ( t ) + h 2 x ˙ T ( t ) R 4 x ˙ ( t ) + h 1 h 1 d h ¯ 1 t x ˙ T ( t h 1 t ) ( R 1 R 2 ) x ˙ ( t h 1 t ) + h h d h ¯ t x ˙ T ( t h t ) ( R 3 R 4 ) x ˙ ( t h t ) h 1 t h 1 t h 1 t x ˙ T ( s ) R 1 x ˙ ( s ) d s + t h 1 t t x ˙ T ( s ) R 2 x ˙ ( s ) d s h t h t h t x ˙ T ( s ) R 3 x ˙ ( s ) d s + t h t t x ˙ T ( s ) R 4 x ˙ ( s ) d s .
According to R i > 0 , ( i = 1 , 2 , 3 , 4 ), letting α = h 1 t h 1 and β = h t h , it follows from Lemmas 1 and 2 that
h 1 t h 1 t h 1 t x ˙ T ( s ) R 1 x ˙ ( s ) d s + t h 1 t t x ˙ T ( s ) R 2 x ˙ ( s ) d s 1 α ξ T ( t ) Γ 1 T R ¯ 1 Γ 1 ξ ( t ) 1 1 α ξ T ( t ) Γ 2 T R ¯ 2 Γ 2 ξ ( t ) ξ T ( t ) Γ 1 T ( R ¯ 1 + α X 1 ) Γ 1 + 2 Γ 1 T [ ( 1 α ) Y 1 + α Y 2 ] Γ 2 + Γ 2 T [ R ¯ 2 + ( 1 α ) X 2 ] Γ 2 ξ ( t ) ,
h t h t h t x ˙ T ( s ) R 3 x ˙ ( s ) d s + t h t t x ˙ T ( s ) R 4 x ˙ ( s ) d s 1 β ξ T ( t ) Γ 3 T R ¯ 3 Γ 3 ξ ( t ) 1 1 β ξ T ( t ) Γ 4 T R ¯ 4 Γ 4 ξ ( t ) ξ T ( t ) Γ 3 T ( R ¯ 3 + α X 3 ) Γ 3 + 2 Γ 3 T [ ( 1 β ) Y 3 + α Y 4 ] Γ 4 + Γ 4 T [ R ¯ 4 + ( 1 β ) X 4 ] Γ 4 ξ ( t ) .
For an appropriately matrix U R 3 n × n , we can get
0 = 2 x T ( t ) x T ( t h t ) x ˙ T ( t ) U A x ( t ) + A 1 x ( t h t ) x ˙ ( t ) = 2 ξ T ( t ) Δ 0 T U Π 0 ξ ( t ) .
Finally, from the above derivation (A2)–(A7), we have
V ˙ ( t ) ξ T ( t ) Ω 0 ( h ˙ t ) + h t Ω 1 ( h ˙ t ) + h t 2 Ω 2 ( h ˙ t ) ξ ( t ) + ξ T ( t ) Ω 10 ( h ˙ 1 t ) + h 1 t Ω 11 ( h ˙ 1 t ) + h 1 t 2 Ω 21 ( h ˙ 1 t ) ξ ( t ) .
According to Lemma 3, the matrix inequality (11) together with the matrix inequality (12) imply that V ˙ ( t ) < 0 . Therefore, by Lyapunov stability theorem, it can guarantee that the time-delayed system (4) is asymptotically stable. The proof is completed. □

References

  1. Kundur, P. Power System Stability and Control; McGraw-Hill: New York, NY, USA, 1994. [Google Scholar]
  2. Mello, R.; Apostolopoulou, D.; Alonso, E. Cost efficient distributed load frequency control in power systems. In Proceedings of the 21st IFAC World Congress, Berlin, Germany, 12–17 July 2020; Volume 53, pp. 8037–8042. [Google Scholar]
  3. Ali, D.; Mohammad, M.; Zeinolabedin, M.; Lieven, V. A novel technique for load frequency control of multi-area power systems. Energies 2020, 13, 2125. [Google Scholar]
  4. Giudice, D.; Brambilla, A.; Grillo, S.; Bizzarri, F. Effects of inertia, load damping and dead-bands on frequency histograms and frequency control of power systems. Int. J. Electr. Power Energy Syst. 2021, 129, 106842. [Google Scholar] [CrossRef]
  5. Ladygin, A.; Bogachenko, D.; Kholin, V. Efficient control of induction motor current in a frequency—Ccontrolled electric drive. Russ. Electr. Eng. 2020, 6, 362–367. [Google Scholar] [CrossRef]
  6. Baykov, D.; Gulyaev, I.; Inshakov, A.; Teplukhov, D. Simulation modeling of an induction motor drive controlled by an array frequency converter. Russ. Electr. Eng. 2019, 7, 485–490. [Google Scholar] [CrossRef]
  7. Aminov, R.; Garievskii, M. Effect of engagement in power and frequency control on the service life of steam–turbine power units. Power Technol. Eng. 2019, 53, 479–483. [Google Scholar] [CrossRef]
  8. Khalil, A.; Peng, A. An accurate method for delay margin computation for power system stability. Energies 2018, 11, 3466. [Google Scholar] [CrossRef] [Green Version]
  9. Jin, L.; Zhang, C.; He, Y.; Jiang, L.; Wu, M. Delay-dependent stability analysis of multi-area load frequency control with enhanced accuracy and computation efficiency. IEEE Trans. Power Syst. 2019, 34, 3687–3696. [Google Scholar] [CrossRef]
  10. Martin, G.; Begovic, M.; Taylor, D. Transient control of generator excitation in consideration of measurement and control delays. IEEE Trans. Power Deliv. 2002, 10, 135–141. [Google Scholar] [CrossRef]
  11. Manikandan, S.; Kokil, P. Stability analysis of load frequency control system with constant communication delays. IFAC PapersOnLine 2020, 53, 338–343. [Google Scholar] [CrossRef]
  12. Zhang, C.; Jiang, L.; Wu, Q.; He, Y.; Wu, M. Further results on delay-dependent stability of multi-area load frequency control. IEEE Trans. Power Syst. 2015, 28, 4465–4474. [Google Scholar] [CrossRef]
  13. Peng, C.; Zhang, J. Delay-distribution-dependent load frequency control of power systems with probabilistic interval delays. IEEE Trans. Power Syst. 2016, 31, 3309–3317. [Google Scholar] [CrossRef]
  14. Ramakrishnan, K.; Ray, G. Stability criteria for nonlinearly perturbed load frequency systems with time-delay. Electr. Power Compon. Syst. 2017, 45, 302–314. [Google Scholar] [CrossRef]
  15. Sönmez, A.; Nwankpa, C. An exact method for computing delay margin for stability of load frequency control systems with constant communication delays. IEEE Trans. Power Syst. 2016, 31, 370–377. [Google Scholar] [CrossRef]
  16. Yang, F.; He, J.; Wang, D. New stability criteria of delayed load frequency control systems via infinite-series-based inequality. IEEE Trans. Ind. Inform. 2018, 14, 231–240. [Google Scholar] [CrossRef]
  17. Chen, B.; Shangguan, X.; Jin, L.; Li, D. An improved stability criterion for load frequency control of power systems with time-varying delays. Energies 2020, 13, 2101. [Google Scholar] [CrossRef] [Green Version]
  18. Shen, H.; Jiao, S.; Park, J.; Sreeram, V. An improved result on H load frequency control for power systems with time delays. IEEE Syst. J. 2020, 15, 3238–3248. [Google Scholar] [CrossRef]
  19. Xu, H.; Zhang, C.; Jiang, L.; Smitha, J. Stability analysis of linear systems with two additive time-varying delays via delay-product-type Lyapunov functional. Appl. Math. Model. 2017, 45, 955–964. [Google Scholar] [CrossRef]
  20. Jiang, L.; Yao, W.; Wu, Q.; Wen, J.; Chen, S. Delay-dependent stability for load frequency control with constant and time-varying delays. IEEE Trans. Power Syst. 2012, 27, 932–941. [Google Scholar] [CrossRef]
  21. Shen, C.; Li, Y.; Zhu, X.; Duan, W. Improved stability criteria for linear systems with two additive time-varying delays via a novel Lyapunov functional. J. Comput. Appl. Math. 2019, 363, 312–324. [Google Scholar] [CrossRef]
  22. Hua, C.; Qiu, Y.; Wang, Y.; Guan, X. An augmented delays-dependent region partitioning approach for recurrent neural networks with multiple time-varying delays. Neurocomputing 2021, 423, 248–254. [Google Scholar] [CrossRef]
  23. Duan, W.; Li, Y.; Chen, J. Further stability analysis for time-delayed neural networks based on an augmented Lyapunov functional. IEEE Access 2019, 7, 104655–104666. [Google Scholar] [CrossRef]
  24. Duan, W.; Li, Y.; Chen, J. An enhanced stability criterion for linear time-delayed systems via new Lyapunov-Krasovskii functionals. Adv. Differ. Equ. 2020, 2020, 21. [Google Scholar] [CrossRef]
  25. Kim, J. Further improvement of Jensen inequality and application to stability of time-delayed systems. Automatica 2016, 64, 121–125. [Google Scholar] [CrossRef]
  26. Scopus, P.; Tian, Y.; Wang, Z. Stability analysis for delayed neural networks based on the augmented Lyapunov-Krasovskii functional with delay-product-type and multiple integral terms. Neurocomputing 2020, 410, 295–303. [Google Scholar]
  27. Gholami, Y. Existence and global asymptotic stability criteria for nonlinear neutral-type neural networks involving multiple time delays using a quadratic-integral Lyapunov functional. Adv. Differ. Equ. 2021, 112, 2021. [Google Scholar]
  28. Duan, W.; Fu, X.; Yang, X. Further results on the robust stability for neutral-type Lur’e system with mixed delays and sector-bounded nonlinearities. Int. J. Control Autom. Syst. 2016, 14, 560–568. [Google Scholar] [CrossRef]
  29. Duan, W.; Fu, X.; Liu, Z. Improved robust stability criteria for time-delay Lur’e system. Int. J. Control Autom. Syst. 2017, 19, 139–150. [Google Scholar] [CrossRef]
  30. Zhang, X.; Han, Q.; Seuret, A.; Gouaisbaut, F. An improved reciprocally convex inequality and an augmented Lyapunov-Krasovskii functional for stability of linear systems with time-varying delay. Automatica 2017, 84, 221–226. [Google Scholar] [CrossRef] [Green Version]
  31. Duan, W.; Li, Y.; Sun, Y.; Chen, J.; Yang, Y. Enhanced master-slave synchronization criteria for chaotic Lur’e systems based on time-delayed feedback control. Math. Comput. Simul. 2020, 177, 276–294. [Google Scholar] [CrossRef]
  32. Gu, K. An integral inequality in the stability problem of time-delay systems. In Proceedings of the Conference of the IEEE Industrial Electronics, Sydney, Australia, 12–15 December 2010. [Google Scholar]
  33. Seuret, A.; Gouaisbaut, F. Hierarchy of LMI conditions for the stability analysis of time-delay systems. Syst. Control Lett. 2015, 81, 1–7. [Google Scholar] [CrossRef] [Green Version]
  34. Feng, W.; Luo, F.; Duan, W.; Chen, J. An improved stability criterion for linear time-varying delay systems. Automatika 2020, 61, 229–237. [Google Scholar] [CrossRef]
  35. Kwon, N.; Lee, S. An affine integral inequality of an arbitrary degree for stability analysis of linear systems with time-varying delays. IEEE Access 2021, 9, 51958–51969. [Google Scholar] [CrossRef]
  36. Zhang, C.; He, Y.; Jiang, L.; Wu, M. Notes on stability of time-delay systems: Bounding inequalities and augmented Lyapunov-Krasovskii functionals. IEEE Trans. Autom. Control 2017, 62, 5331–5336. [Google Scholar] [CrossRef] [Green Version]
  37. Duan, W.; Du, B.; Li, Y.; Shen, C.; Zhu, X.; Li, X.; Chen, J. Improved sufficient LMI conditions for the robust stability of time-delayed neutral-type Lur’e systems. Int. J. Control Autom. Syst. 2018, 16, 2343–2353. [Google Scholar] [CrossRef]
  38. Duan, W.; Li, Y.; Chen, J. New results on stability analysis of uncertain neutral-type Lur’e systems derived from a modified Lyapunov-Krasovskii functional. Complexity 2019, 2019, 1706264. [Google Scholar] [CrossRef]
  39. de Oliveira, F.S.; Souza, F.O. Further refinements in stability conditions for time-varying delay systems. Appl. Math. Comput. 2020, 369, 124866. [Google Scholar] [CrossRef]
Figure 1. The basic diagram of the simplified LFC of one-area power system.
Figure 1. The basic diagram of the simplified LFC of one-area power system.
Energies 14 05820 g001
Figure 2. The basic diagram of the simplified LFC of two-area power system.
Figure 2. The basic diagram of the simplified LFC of two-area power system.
Energies 14 05820 g002
Figure 3. Frequency deviation and control error responses of one-area deregulated LFC scheme with K P = 0.2 , K I = 0.1 and different time-varying delays.
Figure 3. Frequency deviation and control error responses of one-area deregulated LFC scheme with K P = 0.2 , K I = 0.1 and different time-varying delays.
Energies 14 05820 g003
Figure 4. Frequency deviation and control error responses of one-area deregulated LFC scheme with K I = 0.05 , h 1 ( t ) = 0 and different K P , h 2 ( t ) .
Figure 4. Frequency deviation and control error responses of one-area deregulated LFC scheme with K I = 0.05 , h 1 ( t ) = 0 and different K P , h 2 ( t ) .
Energies 14 05820 g004
Figure 5. Frequency deviation and control error responses of one-area deregulated LFC scheme with K I = 0.6 , h 1 ( t ) = 0 and different K P , h 2 ( t ) .
Figure 5. Frequency deviation and control error responses of one-area deregulated LFC scheme with K I = 0.6 , h 1 ( t ) = 0 and different K P , h 2 ( t ) .
Energies 14 05820 g005
Figure 6. Frequency deviation and control error responses of one-area deregulated LFC scheme with h 1 ( t ) = 0 and different K P , K I and time-varying h 2 ( t ) .
Figure 6. Frequency deviation and control error responses of one-area deregulated LFC scheme with h 1 ( t ) = 0 and different K P , K I and time-varying h 2 ( t ) .
Energies 14 05820 g006
Figure 7. Frequency deviation and control error responses of two-area deregulated LFC scheme with K P = 0.05 , K I = 0.4 and time-varying delays.
Figure 7. Frequency deviation and control error responses of two-area deregulated LFC scheme with K P = 0.05 , K I = 0.4 and time-varying delays.
Energies 14 05820 g007
Figure 8. Frequency deviation and control error responses of two-area deregulated LFC scheme with K I = 0.05 , h 1 ( t ) = 0 and different K P , h 2 ( t ) .
Figure 8. Frequency deviation and control error responses of two-area deregulated LFC scheme with K I = 0.05 , h 1 ( t ) = 0 and different K P , h 2 ( t ) .
Energies 14 05820 g008
Figure 9. Frequency deviation and control error responses of two-area deregulated LFC scheme with K I = 1 , h 1 ( t ) = 0 and different K P , h 2 ( t ) .
Figure 9. Frequency deviation and control error responses of two-area deregulated LFC scheme with K I = 1 , h 1 ( t ) = 0 and different K P , h 2 ( t ) .
Energies 14 05820 g009
Figure 10. Frequency deviation and control error responses of two-area deregulated LFC scheme with h 1 ( t ) = 0 and different K P , K I , h 2 ( t ) .
Figure 10. Frequency deviation and control error responses of two-area deregulated LFC scheme with h 1 ( t ) = 0 and different K P , K I , h 2 ( t ) .
Energies 14 05820 g010
Table 1. LFC systems parameters with T 12 = 0.1986 pu/rad.
Table 1. LFC systems parameters with T 12 = 0.1986 pu/rad.
T ch ( s ) β R T g ( S ) D M ( s )
Area 10.3210.050.1110
Area 20.421.50.050.171.512
Table 2. MAUBs h 2 for given h 1 , K P and K I under one-area LFC system.
Table 2. MAUBs h 2 for given h 1 , K P and K I under one-area LFC system.
Methods h 1 1.01.21.5
[20]4.8034.6034.303
[19]5.8825.6825.383
[21]5.9875.8785.559
Th. 16.1556.0466.003
Table 3. MAUBs h 1 for given h 2 , K P and K I under one-area LFC system.
Table 3. MAUBs h 1 for given h 2 , K P and K I under one-area LFC system.
Methods h 2 2.03.04.0
[20]3.8032.8031.803
[19]4.8923.8862.885
[21]4.9973.9793.001
Th. 15.0164.0313.217
Table 4. MAUBs h 2 for h 1 ( t ) = 0 and μ 2 = 0 under one-area LFC system.
Table 4. MAUBs h 2 for h 1 ( t ) = 0 and μ 2 = 0 under one-area LFC system.
μ 2 = 0
K P Methods K I 0.050.10.150.20.40.61
0[20]27.9313.789.066.693.121.910.89
[14]27.9213.776.693.121.910.88
[15]30.9215.209.967.343.382.040.92
[16]30.9115.207.343.392.050.93
[11]30.7915.149.927.313.372.310.92
Th. 130.7915.179.957.333.382.030.92
0.05[20]27.8714.069.286.873.221.970.93
[15]31.8815.6810.287.583.502.120.97
[11]31.7415.6210.247.553.492.120.97
Th. 131.7715.6710.277.573.502.120.97
0.1[14]27.0313.696.943.292.020.96
[20]27.0313.686.943.292.020.96
[16]31.6116.027.793.612.191.01
Th. 132.7316.0810.317.813.622.191.04
Table 5. MAUBs h 2 for h 1 ( t ) = 0 and μ 2 = 0.9 under one-area LFC system.
Table 5. MAUBs h 2 for h 1 ( t ) = 0 and μ 2 = 0.9 under one-area LFC system.
μ 2 = 0.9
K P Methods K I 0.050.10.20.40.61
0[14]20.459.9634.591.811.010.48
[20]26.3712.966.252.851.680.74
[16]27.2613.396.432.911.710.75
[17]27.5013.736.613.021.800.78
Th. 128.0314.247.123.201.910.79
0.1[14]17.399.164.671.851.050.48
[20]20.2511.075.932.871.750.74
[16]22.0012.326.593.111.840.75
[17]29.5114.527.023.231.940.86
Th. 130.0714.617.133.241.950.87
Table 6. MAUBs h 2 for given h 1 , K P and K I under two-area LFC system.
Table 6. MAUBs h 2 for given h 1 , K P and K I under two-area LFC system.
μ 2 0.20.50.8
Methods μ 1 h 1 11.5211.5211.52
Th. 10.12.1221.7741.3232.0511.6371.2732.0151.5821.098
0.22.0181.5591.0642.0031.5071.0181.9971.4991.001
0.51.9921.5291.0301.9731.4810.9951.8951.3980.908
Table 7. MAUBs h 1 for given h 2 , K P and K I under two-area LFC system.
Table 7. MAUBs h 1 for given h 2 , K P and K I under two-area LFC system.
μ 2 0.20.50.8
Methods μ 1 h 2 123123123
Th. 10.12.0281.1830.2511.9391.0780.1891.9251.0510.081
0.22.0161.0270.0791.9211.0170.0771.9101.0090.054
0.52.0091.0120.0531.9071.0070.0371.8971.0070.036
Table 8. MAUBs h 2 for h 1 ( t ) = 0 and μ 2 = 0 under two-area LFC system.
Table 8. MAUBs h 2 for h 1 ( t ) = 0 and μ 2 = 0 under two-area LFC system.
μ 2 = 0
K P Methods K I 0.050.10.150.20.40.61
0[20]27.8513.708.976.603.001.750.57
[15]30.8115.099.847.213.231.840.59
[11]30.69515.0429.8167.1963.2251.8460.585
Th. 130.8315.119.877.243.231.850.62
0.05[20]27.8314.029.216.783.101.810.62
[15]31.2715.5710.167.453.351.920.64
[11]31.64315.52610.1327.4333.3341.9240.638
Th. 131.6715.5510.157.443.341.930.65
0.1[14]13.656.883.171.86
[20]13.656.883.171.86
[16]15.977.673.472.03
[17]16.017.683.472.03
Th. 131.6716.0510.167.713.482.030.65
Table 9. MAUBs h 2 for h 1 ( t ) = 0 and μ 2 = 0.5 under two-area LFC system.
Table 9. MAUBs h 2 for h 1 ( t ) = 0 and μ 2 = 0.5 under two-area LFC system.
μ 2 = 0.5
K P Methods K I 0.20.40.6
0[20]5.552.361.18
[14]6.142.681.40
[16]6.412.811.54
[17]6.662.951.65
Th. 16.692.971.65
0.1[20]5.352.551.30
[14]6.342.831.51
[16]6.752.841.53
[17]7.103.161.78
Th. 17.133.191.78
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Feng, W.; Xie, Y.; Luo, F.; Zhang, X.; Duan, W. Enhanced Stability Criteria of Network-Based Load Frequency Control of Power Systems with Time-Varying Delays. Energies 2021, 14, 5820. https://doi.org/10.3390/en14185820

AMA Style

Feng W, Xie Y, Luo F, Zhang X, Duan W. Enhanced Stability Criteria of Network-Based Load Frequency Control of Power Systems with Time-Varying Delays. Energies. 2021; 14(18):5820. https://doi.org/10.3390/en14185820

Chicago/Turabian Style

Feng, Wenxi, Yanshan Xie, Fei Luo, Xianyong Zhang, and Wenyong Duan. 2021. "Enhanced Stability Criteria of Network-Based Load Frequency Control of Power Systems with Time-Varying Delays" Energies 14, no. 18: 5820. https://doi.org/10.3390/en14185820

APA Style

Feng, W., Xie, Y., Luo, F., Zhang, X., & Duan, W. (2021). Enhanced Stability Criteria of Network-Based Load Frequency Control of Power Systems with Time-Varying Delays. Energies, 14(18), 5820. https://doi.org/10.3390/en14185820

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