Abstract

In this work, novel stability results for load frequency control (LFC) system considering time-varying delays, nonlinearly perturbed load, and time-varying disturbance of system parameters are proposed by using proportional-integral control strategy. Considering the nonlinearly exogenous load disturbance and system parameters disturbance, an improved stability criterion in the form of linear matrix inequalities (LMIs) is derived by novel simple Lyapunov–Krasovskii functionals (LKFs). The delay-dependent matrix in quadratic term, cross terms of variables, and quadratic terms multiplied by 1st, 2nd, and 3rd degrees of scalar functions are included in the new simple LKF. Taking the single-area and two-area LFC system installed with proportional-integral (PI) controller as example, our results surpass the previous maximum allowable size of time delay. Meanwhile, the relationship between time delay varying rate, load disturbance degree, gains of PI controller, and delay margin of the LFC system is researched separately. The results can provide guidance to tune the PI controller for achieving maximum delay margin, in which the LFC system can withstand without losing stability. At last, the simulation results verify the effectiveness and superiority of the proposed stability criterion.

1. Introduction

The usage of load frequency control (LFC) is widely in the power system, which restores the balance of the power system between load demand and generation supply [13]. The sudden change of load can bring threaten to the safe and economic operation of the power system. Thus, the analysis and research about LFC scheme of the power system are considered to be essential and important segment to guarantee the safe operation of the power system. The LFC scheme is also effectively applied in smart grid [4, 5].

A wide and dedicated open communication network is needed to transmit control signals, measurements between remote RTUs, and the control center in the conventional LFC system [3, 6]. The usage of the communication network leads to a series of inevitable problems, like time delays [711], packet losses, and so on. Hence, many researchers focus on the time delay phenomenon of the power control system and the impact of network delay on the communication-based power system [1216]. From the perspective of stability analysis, it is of significant to seek the maximum allowable network delay that the power system with LFC scheme can withstand without losing sable.

Generally speaking, there are mainly two approaches to seek the upper bound of the time delay system, which have different restriction and conservatism. One is the direct approach, which is frequency domain method, such as tracing eigenvalue [17, 18], or cluster treatment of characteristic roots [1921]. The advantage of these direct methods is that the accurate delay margin can be derived by calculating eigenvalues of the whole system. However, the disadvantage of direct methods is that they can only be applied to constant time delays situation. The other is time domain method, which is based on Lyapunov stability theory and linear matrix inequalities (LMIs) techniques [2226]. Though the time domain methods are more conservatism than the frequency domain methods, the time domain methods can be applied in both constant and time-varying delays situation. So, the time domain methods based on Lyapunov stability theory are abroad exploited.

The main effort to diminish the conservatism of time domain methods has focused on two aspects, methods of constructing L-K functional and analyzing techniques for bounding the derivatives of L-K functional with regard to time. The former includes delay-division functional, functional with matrices dependent on the time delays [27], functional with two-integral or triple-integral terms [28], and so on. The latter techniques include improved majorization technique, free weighting matrix method [29], integral inequality including Jesen inequality [23], Wirtinger inequality [30], auxiliary function-based integral inequality [31], and reciprocal convex technique [32]. More and more researchers focus on the delay-dependent criteria for stability of the power system because the delay-dependent criteria have less conservatism than delay-independent criteria [22, 33].

Until now, there are a few excellent research studies about stability analysis of the delayed LFC systems. Ramakrishnan and Ray [34] focused on the delay-dependent robust stability problem for the LFC system with multiple time-invariant delays in framework. The model of the LFC system with both sampling and transmission delay was proposed in [35]. The stability analysis problem was considered in [3], which focused on the LFC system with electric vehicles and time delays. The novel linear operator inequality approach was proposed, and a delay-dependent stability criterion with less conservative was obtained. Unfortunately, the load disturbance and random disturbance of system parameters were not considered either in the stability analysis or simulation in the above literatures. It is worth noting that the bounded nonlinear load perturbation model was built without considering random disturbance of system parameters in literatures [22, 36, 37]. The augmented L-K functional including single and double integral terms was developed to analyze stability of the LFC system in these literatures. The novel LKF with delay-dependent matrix in quadratic term, cross terms of variables, and quadratic terms multiple by a higher degrees of scalar function proposed by this paper contains more information about state variables than the routine L-K functional adopted in [22, 36, 37], and the results in this paper are less conservative than those of [22, 36, 37].

This work proposes a novel stability criterion for the time-varying delay LFC power system considering nonlinearly perturbed and random disturbance of system parameters. Considering unknown exogenous load disturbance and random disturbance of system parameters as constrained time-varying nonlinear function related to current and delayed state vectors, the delayed LFC system with PI controller can be described as a time-varying delay system. Then, the novel stability criterion of the LFC power system is expressed in the form of LMIs by a simple LKF and integral inequalities. The main opinion of constructing the new LKF is the application of delay-dependent matrix, cross terms of variables, and quadratic terms multiplied by 1st, 2nd, and 3rd degrees of scalar function. Moreover, the case studies are carried out taking the single-area and two-area LFC system installed with PI controller as example with simulation examples to verify the effectiveness and superiority of the proposed stability criteria.

2. System Description and Problem Formulation

The notations of jth area for the LFC system are given in Table 1. The subscript j should be omitted for the single-area LFC system. For example, can be changed to in the single-area LFC system. It is worth noting that in the single-area LFC system.

Consider the uncertainties in the power parameters of the actual LFC system. The state-space model with disturbance parameters and load disturbance of the multiarea LFC system (N area) is presented as follows:where represents state vector. and represent the system matrices related to the current state and delayed state vectors. and represent perturbation parameters. is a time-varying matrix, which satisfies . , , and are the constant matrix. represents the network time-varying delay in th area. represents the known matrix related to the disturbance vector. The initial state can be described by a continuous function defined in .

2.1. LFC Modeling

In order to simplify the calculation, it supposes that , which means the multiple delays are supposed to be equal and shown as a unified time-varying delay . So, the state-space model of multiple-area can be presented as follows:where and . and represent the upper bound of time-varying delay and the delay derivative separately. The constraints of and can make time-varying delay systems well-posed because the fast time-varying delay will cause a series of problems, like causality, minimality, inconsistency, and so on.

Furthermore, considering and , the state-space model of multiple-area can be rewritten as follows:

For a typical two-area LFC system, the dynamic model is shown in Figure 1. is network time-varying delay notation. The state vectors, load disturbance vectors, and system matrices are presented as follows:

For the th area, , the notations are presented in Table 1.

The dynamic model of the single-area LFC system is shown in Figure 1 without dotted lines. The state vector, disturbance vector, and system matrices of the single-area LFC system are reduced as follows:

Remark 1. When there are subareas connected by tie-lines forming a multiarea power system, the dotted line will be copied by N times. The decentralized control strategy is applied in the multiarea power system, which means every control area is independent and has its own LFC center to maintain the balance of generation and load. As a result, the interactions between th area and other areas, , are treated as disturbances for th area. Then, define , where (4) is the expanded form of in two-area power system. And the system matrices of multiareas are listed as follows:

2.2. Disturbance Parameters and Load Disturbance Model of LFC System

This paper models the unknown disturbance of power system parameters as , which is nonlinear function related to current and delayed state vector. Owing to the restriction of , the following inequalities can be obtained:

Similarly, the load disturbance is considered as time-varying nonlinear function with current and delayed state vector as , which satisfy the following norm-bounded constraints:where and represent known scalars. and represent known constant matrices with appropriate dimensions. , , , and bound the magnitude of time-varying load disturbance of power system together.

3. Main Results

Before deriving the main results, Lemma 1 should be introduced. Lemma 1 bounds the integral of quadratic function multiplied by a 1st/2nd degree scalar function [26].

Lemma 1. For a positive definite matrix with appropriate dimensions and a vector , the following inequalities can be obtained for all scalar valued function :where and . And and represent appropriate dimensional matrices and represents arbitrary vector. Then, we can obtain the delay-dependent robust stability criterion of system (3).

Theorem 1. Given constants and , the considered system (3) is asymptotically stable if there exist positive matrices , , , and any matrices and satisfying the following LMIs:where represent and , and , and , and and , the four situations, respectively. Moreover, the details about other notations are listed in Appendix.

Proof. The Lyapunov functionals are adopted as follows:where .
Then, define as ; represents the th vector of , e.g., .
Computing the derivative of along the trajectories of system (3):The sum of all integral terms is expressed by :Owing to , the following inequalities can be obtained:Using Lemma 1 to estimate :The following inequality holds for any from (9) and (7), respectively:They can be expressed by the augmented state vector , respectively:By substituting (16) into (13) and adding the inequality (18), the following inequality can be obtained:where have been defined in Appendix, andFinally, considering and , inequality (9) is obtained by Schur complement, which ensures the following inequality:Hence, the system is asymptotically stable. The proof is finished.

Remark 2. The condition of positive matrices , , and in Theorem 1 can guarantee the adopted Lyapunov functionals being positive. Thus, the positive definite of and negative definite of in time domain can guarantee the asymptotic stability of system (3) in .

Remark 3. Our work has investigated the stability condition of LFC system (3) under the limitation and . The limitation conditions are . That is to say, the stability condition should be satisfied under the four situations, , , , and . must be guaranteed under the four situations. Thus, inequality (9) under in Theorem 1 must be satisfied simultaneously for system asymptotic stability.

Remark 4. The main contribution of this section is to propose a novel Theorem 1 in the form of LMIs which can derive a less conservative result of maximum allowable network delay for LFC system (3). The LKF constructed by literature [36] is . Compared with that, the LKF proposed by this paper has four improvement points. (i) The augmented variable , rather than in [36], which can reduce the conservatism by increasing the link between variables. (ii) The time-delay dependent matrix is introduced in quadratic term, rather than constant matrix in [36], which can reduce the conservatism by introducing more messages about time-varying delay. (iii) The cross terms of variable in also can reduce the conservatism by increasing the link between variables. (iv) The 4th, 5th, and 6th terms of the novel LKF multiplied by 1st, 2nd, and 3rd degrees of a scalar function . The increasing of by one signifies the increasing the number of integral by one. Compared with the LKF of literatures [22, 36, 37], we increase the quadratic terms multiplied by 2nd and 3rd degrees of the scalar function , which may reduce the conservatism compared with literatures [22, 36, 37].

Remark 5. If there is a small disturbance or no disturbed load in the LFC system, the delay-dependent robust stability criterion can be developed from Theorem 1 by setting . If there is no or very small disturbance in LFC system parameters, set .

4. Case Studies

In this section, the typical single-area and two-area LFC systems are carried out to analyze the delay margins based on Theorem 1. The parameters of the real two-area LFC system are listed as Table 2. For verifying the less conservatism of our proposed method, the parameters of the two-area LFC system given in [37] are presented in Table 2 for comparison purposes. , , , , , and represent the inertia moments of generator, the generator damping coefficient, the frequency bias factor, time constant of the governor, the turbine time constant, and the speed drop separately. Delay margins with respect to different value of , , gains of PI controller are calculated for constant and time-varying time delays. Then, simulation studies based on MATLAB/Simulink are shown to investigate the effect of time delay on control performance and verify the improvement and effectiveness of the Theorem 1 proposed by this paper.

4.1. Theoretical Delay Margins
4.1.1. Single-Area LFC System

The parameters of Area 1 in Table 2 are adopted as single-area LFC system parameters.

Setting as very small, we investigate the effect of Integral (I) controller on delay margin under three situations of load disturbance (; ; and ) for both constant and time-varying delays. Suppose and in (9) and and in (8) as , where represents the dimension of the system state vector . It can be seen that in the single-area LFC system from (5). Theorem 1 is adopted to calculate the delay margins. The results are shown in Table 3 and Figures 24.

The following conclusions can be obtained from Table 3 and Figures 24:(i)The increasing of and represents that the load fluctuation becomes more pronounced. With the increase in and , the delay margin that PI controller can withstand without decrease in losing stability under the same value of and . The effect of load disturbance on delay-dependent stability of the LFC system can be obtained.(ii)It can be seen from Figures 24 that the increasing of causes the decreasing of delay margin with the same value of under the three situations of load disturbance (; ; ). The effect of time delay changing rate on delay-dependent stability of the LFC system can be obtained.(iii)Simultaneously, it can be seen from Figures 24 that maximum allowable network delay decreases with the increasing of with the same value of under the three situations of load fluctuation.

As for the nominal system condition with constant and time-varying delay , the maximum allowable network delays are listed in Table 4 compared with the other results from literatures [22, 36, 37]. It can be seen that the results calculated by Theorem 1 are less conservative than those of [22, 36, 37], owing to the progressiveness of Lyapunov functional and Lemma 1 adopted by this paper. Therefore, when compared with [22, 36, 37], the results proposed by this paper are more realistic.

Results of the maximum allowable network delays with regard to gains of PI controller are listed in Table 5 and Figure 5 with constant delays and no load disturbance and in Table 6 and Figure 6 with time-varying delays and load disturbance , separately. Results demonstrate that smaller is obtained under the bigger , and a relatively larger is acquired for the smaller . is reduced with the increasing of when . For , increases at first and then decreases with the increasing of . This results implies that for a given , there is a optimal value of which can acquire the maximum , which can give guidance to tune PI controller to achieve maximum delay margins. Moreover, it can be found that of the time-varying delays is smaller and decreases faster than that of constant delays.

4.1.2. Two-Area LFC System

Maximum allowable network delay is calculated by Theorem 1 for the two-area LFC system listed in Table 2 by simplifying the multiple delay as a single delay like (3) and gains of PI controller for two areas as same (, ). In the two-area LFC system, both the values of , , , and are taken as , where .

In two-area LFC system case, for more clearly comparisons with literatures [22, 36, 37], the results are summarized in Table 7 with the condition of under constant and time-varying delay . Similarity, the listed results are less conservative than those of [22, 36, 37]. Also, the results of Table 7 show the effect of PI controller gains on the delay margins.

Also, the results of delay margin with respect to gains of PI controller are listed in Table 8 and Figure 7 with constant delay and no load fluctuation and in Table 9 and Figure 8 with time-varying delay and load fluctuation , respectively. The delay margin changing trend with regard to gains of PI controller of the two-area LFC system is consistent with that of the single-area LFC system. Hence, the details can be omitted here.

4.2. Simulation Verification

The models of the single-area and two-area LFC system are built in MATLAB/Simulink separately. The simulation results verify the above calculation results. The state vectors are selected as the observation since the load fluctuation directly affects the frequency of the LFC system.

4.2.1. Single-Area LFC System

It can be seen from Table 4 that the maximum allowable network delay calculated by Theorem 1 is under the condition of . The time-varying delay is supposed to be a sine function . Figure 9 shows that the state vector converges asymptotically to stable point for ; however, the delay margin calculated by [22] is 4.67, calculated by [36] is 5.94, and calculated by [37] is 6.59. Moreover, if exceeds beyond 6.64, the system becomes unstable as shown in Figure 10 for .

4.2.2. Two-Area LFC System

At first, when , the state-space model of two-area LFC system becomes . Figure 11 shows that the delay-free system is asymptotically stable. It can be seen from Table 9 that the maximum allowable delay calculated by Theorem 1 is under the condition of . Then, it can be seen from Figure 11 that the two-area LFC system is asymptotically stable under . Moreover, if exceeds beyond 6.78, the system becomes unstable as shown in Figure 12 for .

5. Conclusion

The stability criterion for the time-varying delay LFC system with nonlinearly load disturbance and random disturbance of system parameters has been investigated by adopting a novel simple LKF and integral inequality (Lemma 1). Considering the exogenous load disturbance and random disturbance of system parameters as constrained nonlinearly function, a time-varying delay system is modeled to describe the delayed LFC system. The novel simple LKF includes delay-dependent matrix in quadratic term, cross terms of variables, and quadratic terms multiplied by 1st, 2nd, and 3rd degrees of scalar function. Lemma 1 is adopted to bound the quadratic terms multiplied 1st and 2nd degrees of scalar function emerging in the derivative of LKF. The results of the novel simple LKF contributing to reduce conservatism and the relationships between the time delay varying rate, load disturbance degree, the gains of PI controller, and delay margins are derived by the case studied. Meanwhile, simulation verifies the effectiveness and superiority of our results. Nevertheless, the security threats in actual communication network such as denial of service (DOS) attacks and deception attacks are not considered in our work. In the future, how to analyze the stability of the LFC system under the simultaneous existence of various network attack, such as DOS and deception attacks, will be the focus of our following work.

Appendix

Data Availability

The raw/processed data required to reproduce these findings cannot be shared at this time as the data also form part of an ongoing study.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Authors’ Contributions

The manuscript was approved by all authors for publication.

Acknowledgments

This work was supported by the National Natural Science Foundation of China under Grant Nos. 61703060, 12061088, 61802036, and 61873305, the Sichuan Science and Technology Program under Grant No. 21YYJC0469, the project funded by China Postdoctoral Science Foundation under Grant Nos. 2020M683274 and 2021T140092, the Open Research Project of the State Key Laboratory of Industrial Control Technology, Zhejiang University, China, under Grant No. ICT2021B38, the Guangdong Basic and Applied Basic Research Foundation under Grant No. 2021A1515011692, and the Central Government Funds of Guiding Local Scientific and Technological Development for Sichuan Province of China under Grant No. 2021ZYD0012.