Next Article in Journal
A CFD Study on the Effects of Injection Timing and Spray Inclusion Angle on Performance and Emission Characteristics of a DI Diesel Engine Operating in Diffusion-Controlled and PCCI Modes of Combustion
Previous Article in Journal
Anaerobic Digestion of Dairy Effluent in New Zealand, Time to Revisit the Idea?
Previous Article in Special Issue
Using Group Predictive Voltage and Frequency Regulators of Distributed Generation Plants in Cyber-Physical Power Supply Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Computation of Delay Margin Based on Grey Wolf Optimisation for a Load Frequency Control of Two-Area-Network Power Systems

1
Electrical & Electronic Engineering Programme Area, Universiti Teknologi Brunei, Bandar Seri Begawan BE1410, Brunei
2
DTU Engineering Technology, Technical University of Denmark, 2750 Ballerup, Denmark
3
School of Science and Engineering, Politeknik Brunei, Ministry of Education, Bandar Seri Begawan BA1311, Brunei
4
School of Computing and Informatics, Universiti Teknologi Brunei, Bandar Seri Begawan BE1410, Brunei
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Energies 2023, 16(6), 2860; https://doi.org/10.3390/en16062860
Submission received: 15 February 2023 / Revised: 13 March 2023 / Accepted: 16 March 2023 / Published: 20 March 2023
(This article belongs to the Special Issue Advanced Topics in Electrical Power Engineering)

Abstract

:
In classical power systems, frequency measurements are transferred via a specialised communication channel, resulting in time delay. The time delay plays a major role in a power system, which can reduce the dynamic performance of the load–frequency control (LFC) system and can destabilise the system. The research to date has tended to focus on developing a new algorithm to determine the delay margin (DM) rather than looking into a hybrid algorithm which includes a nature-inspired metaheuristic optimisation technique. This paper introduces a novel method for computing the DM based on grey wolf optimisation (GWO), specifically for the constant time delay. In the proposed method, GWO is employed to optimise the minimum error of the spectral radius and to determine the best design variable of the crossing frequency. With the help of the proposed method, the sweeping range is no longer required, which improves the accuracy of the result. To evaluate the proposed method, a two-area network power system is considered as a case study. Furthermore, the effect of the PI controller gains on the DM is taken into account. The proposed method efficacy is demonstrated by comparing it with the most recently published methods. The results demonstrate that the proposed method is remarkably better than the existing methods found in the literature, where the smallest percentage inaccuracy using the simulation-based DM based on GWO is found to be 0.000%.

1. Introduction

The real-time synchronisation of loads and demands is essential for the steady operation of a power system. The LFC technology is used to accomplish this objective. One of the most fundamental issues in managing power systems involves monitoring the frequency to detect any mismatch between the power generation and load. The three primary functions of the LFC system are to maintain a steady frequency, distribute the load among the generators, and manage the schedule for tie-line interchanges [1]. In an LFC system, a communication channel (i.e., between plant and controller) is typically utilised in the feedback loop to transmit and receive signals or information that are used to exercise control, which introduces temporal delays into the feedback path [2]. The delay is unavoidable, particularly if the power system adopts open communications [3]. Furthermore, this affects the performance of the LFC systems, which can cause instability [4,5]. Therefore, a less conservative and accurate delay-dependent stability analysis procedure must be used to compute a maximum DM for the LFC systems for various subsets of controller parameters so that it can be used as a practical guideline for fine-tuning controller parameters during implementation, even with only partial knowledge of network delays.
Several studies have recently been conducted to alleviate the issue of temporal delays in LFC systems. The current state-of-the-art in the literature to determine the DM of the LFC systems are namely the frequency-domain and the time-domain methods. In the frequency domain, the time delay of a system can be identified by determining the crossing frequency of a system, where the roots cross the imaginary axis. Meanwhile, in the time domain, the Lyapunov–Krasovskii functional (LKF) is employed to derive the stability condition in the form of a linear matrix inequality (LMIs).
The computational DM specifically for a constant time delay for LFC was developed by many researchers which will be discussed as follows. In [6], LKF is employed with truncated second-order Bessel–Legendre (BL) inequality to restrict the time derivatives of the states for the single-area and multi-area LFC systems with a PI controller. The stability margin is calculated using LMIs, where a binary search iteration algorithm is used to accurately determine the DM. In [7], the DM was examined by the generalised modified Mikhailov criterion in the time domain for a single-area LFC system with fractional order PI controller. The DM is obtained by solving the fractional polynomial equation and determining the intersection at the imaginary axis. In [8], the BL inequality and LKF are utilised to approximate the functional derivative and determine the largest possible lower bound for the LFC system that incorporates a PID controller. The time delay stability margin is solved using the MATLAB LMI toolbox. In [9], a novel augmented Lyapunov functional is constructed with the lower order and the derivative is bounded with the Wirtinger inequality for multi-area LFC schemes. In [10], LKF is employed with quadratic generalised free-weighting matrix inequality to tighten the boundary of the time derivatives for a single-area LFC system. In [11], LKF is employed with single- and double-integral terms for accurate integral inequalities and a new nonlinear optimisation technology was employed to solve the problem of nonlinear time-delay inequalities. In [12], LKF is employed with a triple integral term and using the method of integral inequality for one area and two area LFC. In [13], a modified LKF with delay-dependent matrices containing single-integral items and a unique negative definite inequality equivalent transformation lemma is utilised to translate the nonlinear inequality into LMIs for one and two-area LFC systems. In [14], a novel LKF is employed with quadratic terms multiplied by first, second, and third degrees of scalar functions for one and two-area LFC systems with the PI controller. In [15], a modified LKF is employed, involving the delay-dependent non-integral items and some augmented single-integral items for a single-area power system. In [16], the computation of DM is determined by converting the transcendental equation into nonlinear equations by introducing a new variable. It is reported that the results are more accurate, where the crossing frequency is determined by implementing the sweeping test. It has been shown in prior work that the LFC system has a time delay, and that computing it is crucial for practical purposes due to the majority of existing approaches having advanced complex mathematical analysis and involving numerous computations to determine the DM; however, engineers require a method that is quick and precise in practice [17].
A precise frequency domain technique for DM calculation is described in [18], where the stability condition is constructed and the characteristic equation is translated into a polynomial equation. Exact DM values are acquired by precisely determining the crossing frequencies. The sweeping test in [16,19] determines the precise values of the DM, however, the test is highly dependent on the sweeping range and the sweeping step size. It must be noted that employing the sweeping test with high accuracy (i.e., very small step size) can be time-consuming. Furthermore, the majority of researchers are focused on building a new algorithm for determining the exact DM through the iteration process. However, there is not much effort made to combine optimisation approaches with the existing DM algorithms. This study aimed to eliminate the iteration process and to determine the more accurate results of DM by introducing an optimisation technique. It is expected that by combining an optimisation technique with the existing DM algorithm, the DM results for a constant time delay problem will be more accurate and reliable. Furthermore, it has been proven in [20,21,22,23,24,25,26,27] that the combination of a metaheuristic optimisation technique and a conventional algorithm provides a significant improvement in many power system applications and other real-world problems.
In this paper, the method of improving the computation of the DM using GWO for a constant time delay is presented. The proposed method is developed to optimise the minimum error and thus, determine the best design variable (i.e., the crossing frequency). The crossing frequency obtained from the proposed method will lead to an accurate result of the DM. Furthermore, the contribution of this paper is that the proposed method can address the problem, as mentioned in [17], where the proposed method does not require any sweeping range, sweeping step size, or iteration method to determine the accurate crossing frequency. In addition, an objective function is introduced and formulated using the existing theorem to determine the optimal crossing frequency with minimal error. The proposed method is realised in a two-area network power system, and the DM values are compared with the methods in [18,28] for the different values of K P and K I . It is observed during the simulation that the DM results are more accurate and stable compared to [18,19,28,29]. The simulation was run for a prolonged period of time to avoid inaccuracy in the results. The remainder of the article is organised as follows: The dynamic model of a two-area LFC system with time delay is defined in Section 2, the overview of GWO is detailed in Section 3, the proposed computation of DM based on GWO is described in Section 4, the results are discussed in Section 5, and this article is concluded in Section 6.

2. Dynamic Model of Two-Area LFC System with Time Delay

Figure 1 depicts a standard model of a two-area LFC system. The network is made up of a governor, a turbine, a rotating mass, a PI controller, communication delay, and a load with the feedback of regulation constant and area control error (i.e., A C E 1 and A C E 2 ) for both areas (i.e., area 1 and area 2). The disturbance in the system happens when the load changes, as indicated by Δ P d 1 and Δ P d 2 . When the load on the system increases, the turbine speed decreases until the governor can adjust the steam input into the new load. The error signal decreases as the variation in speed value decreases, and the governor setting approaches the point required to keep the speed constant.
In terms of its state-space model, the model of Figure 1 with communication delay is defined as [28,30]:
x ˙ ( t ) = K x ( t ) + K d x ( t τ ) + F Δ P d
x ( t ) = [ Δ f 1 ( t ) Δ P m 1 ( t ) Δ P v 1 ( t ) A C E 1 d t Δ P 12 ( t ) Δ f 2 ( t ) Δ P m 2 ( t ) Δ P v 2 ( t ) A C E 2 d t ] T
K = D 1 M 1 1 M 1 0 0 1 M 1 0 0 0 0 0 1 T c h 1 1 T c h 1 0 0 0 0 0 0 1 R 1 T g 1 0 1 T g 1 0 0 0 0 0 0 β 1 0 0 0 1 0 0 0 0 2 π T 12 0 0 0 0 2 π T 12 0 0 0 0 0 0 0 1 M 2 D 2 M 2 1 M 2 0 0 0 0 0 0 0 0 1 T c h 2 1 T c h 2 0 0 0 0 0 0 1 R 2 T g 2 0 1 T g 2 0 0 0 0 0 1 β 2 0 0 0
K d = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 β 1 K P 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 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 K P 2 T g 2 β 2 K P 2 T g 2 0 0 K I 2 T g 2 0 0 0 0 0 0 0 0 0
F = 1 M 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 M 2 0 0 0 T
where Δ f 1 and Δ f 2 are the frequency deviations in area 1 and area 2, Δ P m 1 and Δ P m 2 are the mechanical output power deviation in area 1 and area 2, Δ P v 1 and Δ P v 2 are the valve position deviation in area 1 and area 2, Δ P 12 is the power transfer deviation in tie-line between the two areas, D 1 and D 2 are the generator damping coefficient in area 1 and area 2, M 1 and M 2 are the moment of inertia of generator in area 1 and area 2, T c h 1 and T c h 2 are the turbine time constant in area 1 and area 2, T g 1 and T g 2 are the governor time constant in area 1 and area 2, R 1 and R 2 are the speed drop in area 1 and area 2, β 1 and β 2 are the frequency bias factor in area 1 and area 2, T 12 is the synchronising coefficient of the tie-line between two areas, K P 1 and K P 2 are the proportional gain of the PI controller in area 1 and area 2, and K I 1 and K I 2 are the integral gain of the PI controller in area 1 and area 2. The area control error, ACE, is calculated for each control region as the sum or difference of the power transfer along the tie-line connecting the two areas and scaled by a frequency bias factor, which can be defined as [18,29]:
A C E 1 ( t ) = β 1 Δ f 1 ( t ) + Δ P 12 ( t )
A C E 2 ( t ) = β 2 Δ f 2 ( t ) Δ P 12 ( t )
The time delay τ is made up of the following components: the transducer delay, the delay in the communication connection, the delay in processing, the delay in multiplexing, and the delay in the analogue-to-digital conversion [31]. In dedicated communication channels, this delay is constant, but if open communication is utilised, it may fluctuate in length [29]. The two-area LFC system becomes a linear time delay system when Δ P d is set to zero, Equation (1) may then be stated as follows [16]:
x ˙ ( t ) = K x ( t ) + K d x ( t τ )
The maximum DM, τ d , can be determined by transforming Equation (8) into the s-domain using the Laplace transform and can be expressed as:
s I K K d e s τ = 0
The system is asymptotically stable for the given delay [16], if all of the roots of Equation (9) are situated on the left half plane. All of the roots of the free delay system are on the left half plane, and it is assumed the system is stable. In this paper, the frequency domain approach is employed to determine the maximum time delay using grey wolf optimisation, which will be discussed in the next section.

3. Overview of Grey Wolf Optimisation

A pack of grey wolves surrounding and assaulting prey for food is a simplistic depiction of the grey wolf optimisation algorithm, in which they have a disciplined social dominating hierarchy beginning with alpha, beta, delta, and omega [32]. Grey wolves’ hunting behaviour may be separated into three stages, namely the tracking, surrounding, and attacking prey. To begin, surrounding the prey may be defined as [32,33]:
D V = | C V . X V p ( t ) X V ( t ) |
X V ( t + 1 ) = X V p ( t ) A V . D V
where t denotes the iteration, X V p is the position vector of the prey, X V is the position vector of a grey wolf, A V and C V are the coefficient vectors, and D V is the distance between the position vectors of the prey and a grey wolf. These coefficients can be calculated by:
A V = 2 b . r 1 b
C V = 2 . r 2
where b decreases linearly over the length of iterations from 2 to 0 and r 1 , r 2 are random vectors in [0, 1]. Meanwhile, grey wolf hunting behaviour may be analytically described as [32,33]:
X V 1 = X V α A V 1 . ( D V α )
X V 2 = X V β A V 2 . ( D V β )
X V 3 = X V δ A V 3 . ( D V δ )
X V ( t + 1 ) = X V 1 + X V 2 + X V 3 3
where subscript α indicates the best candidate for searching the prey. Meanwhile, the subscripts β and δ represent the second and third best solution, respectively.

4. Delay Margin Computation Based on Grey Wolf Optimisation

As discussed in Reference [16], systems that deal with time delays can either be delay-dependent or delay-independent. A delay-independent system is said to be asymptotically stable only if the value of time delay is positive. However, in a delay-dependent system, the system is said to be asymptotically stable, marginally stable, and asymptotically unstable, if the value of τ < τ d , τ = τ d , and τ > τ d , respectively. The LFC system modelled in Equation (8) is said to be asymptotically stable independent of delay if it satisfies the condition stipulated in [34], as described in Equations (18) and (19):
d e t ( j ω I K K d e j ω τ ) 0 , ω > 0 , τ [ 0 , )
lim ω ρ ( ( j ω I K ) 1 K d ) = 0
If Equations (18) and (19) are satisfied, this means that the system does not intersect the imaginary axis. However, the system is said to be delay-dependent stable, if Equations (18) and (19) do not meet for some values of ω . Where for τ < τ d , all roots lie on the closed left half plane, for τ > τ d , some roots will lie on the right plane, and for τ = τ d , the roots will intersect with the imaginary axis. The following definition, as explained in Reference [16], uses the spectral radius to determine the crossover frequency, ω c , at which the roots cross the imaginary axis.
Definition 1
([34]). A pair of matrices’ combined spectral radius can be described as:
ρ ̲ ( K , K d ) : = m i n { | λ | | d e t ( K λ K d = 0 ) | }
where λ i ( K ) and λ i ( K , K d ) are the ith eigenvalue of the matrix K and the generalised eigenvalue of matrix pair K and K d , respectively.
To determine the maximum DM, we use the sweeping test [34] in the ω domain. The sweeping test is a very useful technique, especially in light of the improvements in computational power seen in modern computers. The sweeping test is superior since it requires less computation and produces accurate results. We apply the following theorem to determine the DM of the two-area LFC system.
Theorem 1
([34]). We can define a stable system (8) at τ d = 0, where K + K d is stable and r a n k ( K d ) = q.
τ ¯ i : = min 1 k n θ k i ω k i , i f λ i ( j ω k i I K , K d ) = e j θ k i f o r s o m e ω k i ( 0 , ) , θ k i [ 0 , 2 π ] , ρ ̲ ( j ω I K , K d ) > 1 ω ( 0 , )
Then,
τ ¯ d : = min 1 i q τ ¯ i
where the system in (8) is unstable at τ = τ d and stable for all τ∈ [ 0 , ) .
Proof 
([34,35,36,37]). The system (8) is considered stable for any time delay when Equation (22) is satisfied:
ρ ̲ ( j ω I K , K d ) = ρ ̲ ( j ω I K , K d e j ω τ ) > 1 , for ω > 0 , τ 0
According to the condition of Equation (22), the system is stable when τ = 0, det( K + K d ) ≠ 0, and ω = 0. However, it is unstable for a certain τ for τ d < . Then, the following assumption is examined:
det ( j ω I K K d e j ω τ ) 0 , ω [ 0 , )
This is true for ω ω k i , or else:
| λ i ( j ω I K , K d ) | 1 , i = 1 , 2 , 3 , . . . , m
For τ [ 0 , τ d ) , τ ω k i θ k i , Equation (25) is satisfied:
det ( j ω k i I K K d e j ω k i τ ) 0
Equation (26) is satisfied if τ = τ d , such that ( ω k i , θ k i ) is unique, which satisfies τ d = θ k i / ω k i , and consequently:
det ( j ω k i I K K d e j ω k i τ d ) = det ( j ω k i I K K d e j θ k i τ d ) = 0
Corollary 1
([34]). If and only if the system (8) is stable regardless of delay:
  • K is stable;
  • K + K d is stable, and;
  • ρ ̲ ( j ω I K , K d ) > 1, ω > 0.
The delay-independent stability is represented by the three conditions in Corollary 1, where the first and second conditions state that the system is stable for τ = 0 and τ = , respectively. The third condition states that the system is stable for τ [ 0 , ) . Finally, we can compute the exact DM by using Theorem 1 to determine the system stability (i.e., delay-independent or delay-dependent stability). First, we can examine the condition in Equation (22) to see if the system (8) is delay independent. For some values of ω , if Equation (22) is not satisfied, the crossing frequency is determined using the proposed algorithm in Figure 2.
To begin, the system matrices for a two-area network are calculated using Equations (3) and (4). This is followed by the objective function to minimise the error of the spectral radius, ϵ , and can be expressed as:
f = m i n ( ϵ )
ϵ = | 1 ρ ̲ ( j ω I K , K d ) |
where integer 1 is the desired spectral radius at the crossing frequency, ω c . Furthermore, ω , is the only decision variable in this problem. The constraint can be expressed as:
0 ω ω m a x
where ω m a x is the maximum frequency desired by the user. Once the GWO algorithm found the optimal solution in the search space, the best design variable ω is then updated to determine the value of the crossing angle, θ , which can be expressed as:
λ i ( j ω I K , K d ) = e j θ
Finally, the desired DM is calculated as:
τ d = min 1 k n θ k i ω k i
In the event that the algorithm is unable to identify the optimum solution, the decision variable is randomly re-initialised in an effort to locate the optimal solution in the search space.

5. Case Study: Two-Area LFC System

Three case studies for a two-area LFC system are presented in this section. We utilise the parameters in References [18,28,29] to compare the outcomes of our suggested approach with those that have already been published. Table 1 tabulates the system DM for various PI controller gains, K P and K I , as well as the results of the methods in references [18,19,28,29]. It is noteworthy that the approach described in reference [18] offers the highest degree of accuracy-reported DMs. Moreover, the GWO algorithm has a maximum population size and maximum iterations of 50 and 200, respectively.
Table 1 shows that the proposed method yields remarkably similar results to the method described in Reference [18]. Furthermore, the proposed method provides fewer computations by reducing the mathematical complexity to obtain the system DMs. Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8 and Figure 9 compare the proposed method to the methods given in References [18,28]. For various values of proportional gain, K P , they show the DM, τ d , versus the integral gain, K I , curve. It is obvious that the proposed method precisely and accurately calculates the DM for the two-area LFC system for a variety of K I and K P values. Furthermore, the relative inaccuracy of the suggested approach and the one stated in [18] is rather minor. More importantly, the average relative percentage error is only 3.1611%, and thus, demonstrates that the DM values are nearly identical to those provided in [18].
Case study 1 with K P = 0.2 and K I = 0.4: The DM estimated using the method provided in [18] is 3.631 s, while it was 3.641 s with the proposed method. The proposed method, just like the method described in [18], gives precise DM values. Simulations with Matlab/Simulink are performed to validate the results. Figure 10 and Figure 11 showcase the frequency response of a two-area LFC system when subjected to different time delay values in area 1 and area 2, while utilising K P = 0.2 and K I = 0.4. At 10 s, the load changes, ( Δ P d 1 = Δ P d 2 ), by 0.1 p.u. Figure 10 and Figure 11 depict the frequency response with time delay, τ , of 3.5 s, 3.641 s, and 3.7 s, for area 1 and area 2, respectively. It is obvious that the system is stable with τ = 3.5 s and unstable with τ = 3.7 s. Meanwhile, the system is marginally stable with τ = 3.641 s.
The percentage inaccuracy while employing the simulation-based DM is 0.000%. Figure 12 depicts the spectral radius against the frequency for K P = 0.2 and K I = 0.4. It was found that the crossing frequency, ω c , was equal to 0.41914 rad/s with θ = 1.5261 rad. By solving Equation (31), we obtain the DM equal to 3.641 s.
Case study 2 with K P = 0.4 and K I = 0.6: The DM estimated using the method provided in [18] is 2.184 s, while it was 2.189 s with the proposed method. The frequency response of the two-area LFC system was illustrated in Figure 13 and Figure 14, taking into account the impact of various time delays on area 1 and area 2, while setting K P = 0.4 and K I = 0.6.
Figure 13 and Figure 14 depict the frequency response with time delay, τ , of 2.0 s, 2.189 s, and 2.2 s, for area 1 and area 2, respectively. It is obvious that the system is stable with τ = 2.0 s and unstable with τ = 2.2 s. Meanwhile, the system is marginally stable with τ = 2.189 s. The percentage inaccuracy while employing the simulation-based DM is 0.000%. Figure 15 depicts the spectral radius against the frequency for K P = 0.4 and K I = 0.6. It was found that the crossing frequency, ω c , was equal to 0.71530 rad/s with θ = 1.5658 rad. By solving Equation (31), we obtain the DM equal to 2.189 s.
Case study 3 with K P = 0.05 and K I = 0.05: The DM estimated using the method provided in reference [18] is 31.704 s, while it was 31.895 s with the proposed method. The frequency response of the two-area LFC system was illustrated in Figure 16 and Figure 17, taking into account the impact of various time delays on area 1 and area 2, while setting K P = 0.05 and K I = 0.05.
Figure 16 and Figure 17 depict the frequency response with time delay, τ , of 30.0 s, 31.800 s, and 33.0 s, for area 1 and area 2, respectively. It is obvious that the system is stable with τ = 30.0 s and unstable with τ = 33.0 s. Meanwhile, the system is marginally stable with τ = 31.800 s. The percentage inaccuracy while employing the simulation-based DM is 0.2987%. Figure 18 depicts the spectral radius against the frequency for K P = 0.05 and K I = 0.05. It was found that the crossing frequency, ω c , was equal to 0.05011 rad/s with θ = 1.5983 rad. By solving Equation (31), we obtain the DM equal to 31.895 s.
Four other published methods were compared with the proposed method. As indicated in [28,29], the DM method produces less traditional results than the LMI method, but the LMI method has the advantage of handling delays that fluctuate over time. The proposed method is easier to execute, yet the DM results are similar to those obtained using the frequency domain techniques, as described in [18,19,28,29].

6. Conclusions

In this paper, a metaheuristic algorithm, namely the GWO, has been successfully employed to determine the optimal DM for a constant time delay of a two-area network. The optimisation of the spectral radius subjected to the optimal frequency formulated with the GWO is regarded as the most effective of all the techniques listed as it provides the minimum error of the spectral radius and hence optimises the best design variable of the crossing frequency such that an accurate result of DM can be obtained. A case study employing a two-area LFC system suggested that the proposed method is significantly superior in comparison to the others, as the sweeping range and sweeping step size method in the frequency domain method are no longer needed to determine the crossing frequency. In comparison to the previous methods, the time delay simulation showed that this novel technique provides precise delay margins with the smallest percentage error of 0.000%. In future work, a detailed case study of time-varying delay for load–frequency control with the proposed method will be implemented and compared with the most recently published methods.

Author Contributions

Conceptualisation, M.H.I., A.S.P. and A.K.; Formal analysis, M.H.I. and M.I.R.; Funding acquisition, A.S.P., M.N.D., K.H.L., S.Y., M.I.R. and T.W.A.; Investigation, M.H.I. and M.N.D.; Methodology, M.H.I. and A.S.P.; Supervision, A.S.P.; Validation, M.H.I., A.S.P., M.N.D., A.K. and M.I.R.; Visualisation, K.H.L. and S.Y.; Writing—original draft, M.H.I. and A.S.P.; Writing—review and editing, M.H.I., A.S.P., M.N.D. and M.I.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Saadat, H. Power System Analysis; McGraw-Hill Companies: New York, NY, USA, 1999. [Google Scholar]
  2. Xu, X.; Tomsovic, K. Application of Linear Matrix Inequalities for Load Frequency Control with Communication Delays. IEEE Trans. Power Syst. 2004, 19, 1508–1515. [Google Scholar]
  3. Mak, K.H.; Holland, B.L. Migrating electrical power network SCADA systems to TCP/IP and Ethernet networking. Power Eng. J. 2002, 16, 305–311. [Google Scholar]
  4. Gu, K.; Niculescu, S.I. Survey on Recent Results in the Stability and Control of Time-Delay Systems. J. Dyn. Sys. Meas. Control. 2003, 125, 158–165. [Google Scholar] [CrossRef]
  5. Yan, S.; Gu, Z.; Park, J.H.; Xie, X.; Dou, C. Probability-density-dependent load frequency control of power systems with random delays and cyber-attacks via Circuital implementation. IEEE Trans. Smart Grid 2022, 13, 4837–4847. [Google Scholar] [CrossRef]
  6. Yang, F.; He, J.; Pan, Q. Further Improvement on Delay-Dependent Load Frequency Control of Power Systems via Truncated B-L Inequality. IEEE Trans. Power Syst. 2018, 33, 5062–5071. [Google Scholar] [CrossRef]
  7. Çelik, V.; Özdemir, M.T.; Lee, K.Y. Effects of Fractional-Order Pi Controller on Delay Margin in Single-Area Delayed Load Frequency Control Systems. J. Mod. Power Syst. Clean Energy 2018, 7, 380–389. [Google Scholar] [CrossRef] [Green Version]
  8. Shen, L.; Xiao, H. Delay-Dependent Robust Stability Analysis of Power Systems with PID Controller. Chin. J. Electr. Eng. 2019, 5, 79–86. [Google Scholar] [CrossRef]
  9. Jin, L.; Zhang, C.-K.; 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. 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. 2020, 363, 312–324. [Google Scholar] [CrossRef]
  11. Chen, W.; Gao, F.; Liu, G. New Results on Delay-Dependent Stability for Nonlinear Systems with Two Additive Time-Varying Delays. Eur. J. Control. 2021, 58, 123–130. [Google Scholar] [CrossRef]
  12. Arunagirinathan, S.; Muthukumar, P. New Asymptotic Stability Criteria for Time-Delayed Dynamical Systems with Applications in Control Models. Results Control. Optim. 2021, 3, 100014. [Google Scholar] [CrossRef]
  13. 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. [Google Scholar] [CrossRef]
  14. Liu, X.; Shi, K.; Tang, L.; Tang, Y. Further Improvement on Stability Results for Delayed Load Frequency Control System via Novel Lkfs and Proportional-Integral Control Strategy. Discret. Dyn. Nat. Soc. 2022, 2022, 4798415. [Google Scholar] [CrossRef]
  15. Li, Y.; Qiu, T.; Yang, Y. Delay-Dependent Stability Criteria for Linear Systems with Two Additive Time-Varying Delays. Int. J. Control. Autom. Syst. 2022, 20, 392–402. [Google Scholar] [CrossRef]
  16. Khalil, A.; Swee Peng, A. A New Method for Computing the Delay Margin for the Stability of Load Frequency Control Systems. Energies 2018, 11, 3460. [Google Scholar] [CrossRef] [Green Version]
  17. Khalil, A.; Laila, D.S. An Accurate Method for Computing the Delay Margin in Load Frequency Control System with Gain and Phase Margins. Energies 2022, 15, 3434. [Google Scholar] [CrossRef]
  18. Sonmez, S.; Ayasun, S.; Nwankpa, C.O. 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]
  19. Thangaiah, J.M.; Parthasarathy, R. Delay-Dependent Stability Analysis of Power System Considering Communication Delays. Int. Trans. Electr. Energy Syst. 2016, 27, e2260. [Google Scholar] [CrossRef]
  20. Ibrahim, M.H.; Ang, S.P.; Dani, M.N.; Rahman, M.I.; Petra, R.; Sulthan, S.M. Optimizing Step-Size of Perturb & Observe and Incremental Conductance MPPT Techniques Using PSO for Grid-Tied PV System. IEEE Access 2023, 11, 13079–13090. [Google Scholar] [CrossRef]
  21. Shehata, A.A.; Tolba, M.A.; El-Rifaie, A.M.; Korovkin, N.V. Power System Operation Enhancement using a new hybrid methodology for optimal allocation of FACTS devices. Energy Rep. 2022, 8, 217–238. [Google Scholar] [CrossRef]
  22. Sheikhzadeh-Baboli, P.; Assili, M. A hybrid adaptive algorithm for power system frequency restoration based on proposed emergency demand side management. IETE J. Res. 2021, 1–13. [Google Scholar] [CrossRef]
  23. Nadimi-Shahraki, M.H.; Taghian, S.; Mirjalili, S.; Abualigah, L.; Abd Elaziz, M.; Oliva, D. EWOA-OPF: Effective Whale Optimization Algorithm to Solve Optimal Power Flow Problem. Electronics 2021, 10, 2975. [Google Scholar] [CrossRef]
  24. Nadimi-Shahraki, M.H.; Taghian, S.; Mirjalili, S.; Faris, H. MTDE: An effective multi-trial vector based differential evolution algorithm and its applications for engineering design problems. Appl. Soft Comput. 2020, 97, 106761. [Google Scholar] [CrossRef]
  25. Nadimi-Shahraki, M.H.; Taghian, S.; Mirjalili, S. An improved grey wolf optimizer for solving engineering problems. Expert Syst. Appl. 2021, 166, 113917. [Google Scholar] [CrossRef]
  26. Nadimi-Shahraki, M.H.; Moeini, E.; Taghian, S.; Mirjalili, S. DMFO-CD: A Discrete Moth-Flame Optimization Algorithm for Community Detection. Algorithms 2021, 14, 314. [Google Scholar] [CrossRef]
  27. Nadimi-Shahraki, M.H.; Taghian, S.; Mirjalili, S.; Ewees, A.A.; Abualigah, L.; Abd Elaziz, M. MTV-MFO: Multi-Trial Vector-Based Moth-Flame Optimization Algorithm. Symmetry 2021, 13, 2388. [Google Scholar] [CrossRef]
  28. Ramakrishnan, K. Delay-dependent stability criterion for delayed load frequency control systems. In Proceedings of the IEEE Annual India Conference (INDICON), Bangalore, India, 16–18 December 2016; pp. 1–6. [Google Scholar]
  29. 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. 2012, 27, 932–941. [Google Scholar] [CrossRef]
  30. Manikandan, S.; Kokil, P. Delay-Dependent Stability Analysis of Network-Based Load Frequency Control of One and Two Area Power System with Time-Varying Delays. Fluct. Noise Lett. 2019, 18, 1950007. [Google Scholar] [CrossRef]
  31. Naduvathuparambil, B.; Valenti, M.C.; Feliachi, A. Communication delays in wide area measurement systems. In Proceedings of the Thirty-Fourth Southeastern Symposium on System Theory, Huntsville, Alabama, 18–19 March 2002; pp. 118–122. [Google Scholar]
  32. Mirjalili, S.; Mirjalili, S.M.; Lewis, A. Grey Wolf Optimizer. Adv. Eng. Softw. 2014, 69, 46–61. [Google Scholar] [CrossRef] [Green Version]
  33. Ibrahim, M.H.; Ang, S.P.; Hamid, Z.; Mohammed, S.S.; Law, K.H. A Comparative Hybrid Optimisation Analysis of Load Frequency Control in a Single Area Power System Using Metaheuristic Algorithms and LQR. In Proceedings of the 2022 International Conference on Green Energy, Computing, and Sustainable Technology (GECOST 2022), Miri, Sarawak, Malaysia, 26–28 October 2022; pp. 227–232. [Google Scholar]
  34. Gu, K.; Kharitonov, V.L.; Chen, J. Stability of Time-Delay Systems; Springer: Berlin/Heidelberg, Germany, 2003. [Google Scholar]
  35. Chen, J. On Computing the Maximal Delay Intervals for Stability of Linear Delay Systems. IEEE Trans. Autom. Control 1995, 40, 1087–1093. [Google Scholar] [CrossRef]
  36. Chen, J.; Latchman, H.A. Frequency Sweeping Tests for Stability Independent of Delay. IEEE Trans. Autom. Control 1995, 40, 1640–1645. [Google Scholar] [CrossRef]
  37. Chen, J.; Gu, G.; Nett, C.N. A New Method for Computing Delay Margins for Stability of Linear Delay Systems. In Proceedings of the IEEE 33rd Conference on Decision and Control, Lake Buena Vista, FL, USA, 14–16 December 1994; pp. 433–437. [Google Scholar]
Figure 1. Two-area LFC scheme [28,30].
Figure 1. Two-area LFC scheme [28,30].
Energies 16 02860 g001
Figure 2. Flowchart of the proposed DM computation based on GWO.
Figure 2. Flowchart of the proposed DM computation based on GWO.
Energies 16 02860 g002
Figure 3. With K P = 0, the DM, τ d , is shown against the integral gain, K I , curve [18,28].
Figure 3. With K P = 0, the DM, τ d , is shown against the integral gain, K I , curve [18,28].
Energies 16 02860 g003
Figure 4. With K P = 0.05, the DM, τ d , is shown against the integral gain, K I , curve [18,28].
Figure 4. With K P = 0.05, the DM, τ d , is shown against the integral gain, K I , curve [18,28].
Energies 16 02860 g004
Figure 5. With K P = 0.1, the DM, τ d , is shown against the integral gain, K I , curve [18,28].
Figure 5. With K P = 0.1, the DM, τ d , is shown against the integral gain, K I , curve [18,28].
Energies 16 02860 g005
Figure 6. With K P = 0.2, the DM, τ d , is shown against the integral gain, K I , curve [18,28].
Figure 6. With K P = 0.2, the DM, τ d , is shown against the integral gain, K I , curve [18,28].
Energies 16 02860 g006
Figure 7. With K P = 0.4, the DM, τ d , is shown against the integral gain, K I , curve [18,28].
Figure 7. With K P = 0.4, the DM, τ d , is shown against the integral gain, K I , curve [18,28].
Energies 16 02860 g007
Figure 8. With K P = 0.6, the DM, τ d , is shown against the integral gain, K I , curve [18,28].
Figure 8. With K P = 0.6, the DM, τ d , is shown against the integral gain, K I , curve [18,28].
Energies 16 02860 g008
Figure 9. With K P = 1.0, the DM, τ d , is shown against the integral gain, K I , curve [18,28].
Figure 9. With K P = 1.0, the DM, τ d , is shown against the integral gain, K I , curve [18,28].
Energies 16 02860 g009
Figure 10. Frequency deviation for LFC−1: Δ f 1 was measured while using K P = 0.2 and K I = 0.4.
Figure 10. Frequency deviation for LFC−1: Δ f 1 was measured while using K P = 0.2 and K I = 0.4.
Energies 16 02860 g010
Figure 11. Frequency deviation for LFC−2: Δ f 2 , was measured while using K P = 0.2 and K I = 0.4.
Figure 11. Frequency deviation for LFC−2: Δ f 2 , was measured while using K P = 0.2 and K I = 0.4.
Energies 16 02860 g011
Figure 12. The spectral radius, ρ , against the frequency, ω , for K P = 0.2 and K I = 0.4.
Figure 12. The spectral radius, ρ , against the frequency, ω , for K P = 0.2 and K I = 0.4.
Energies 16 02860 g012
Figure 13. Frequency deviation for LFC−1: Δ f 1 was measured while using K P = 0.4 and K I = 0.6.
Figure 13. Frequency deviation for LFC−1: Δ f 1 was measured while using K P = 0.4 and K I = 0.6.
Energies 16 02860 g013
Figure 14. Frequency deviation for LFC−2: Δ f 2 was measured while using K P = 0.4 and K I = 0.6.
Figure 14. Frequency deviation for LFC−2: Δ f 2 was measured while using K P = 0.4 and K I = 0.6.
Energies 16 02860 g014
Figure 15. The spectral radius, ρ , against the frequency, ω , for K P = 0.4 and K I = 0.6.
Figure 15. The spectral radius, ρ , against the frequency, ω , for K P = 0.4 and K I = 0.6.
Energies 16 02860 g015
Figure 16. Frequency deviation for LFC−1: Δ f 1 was measured while using K P = 0.05 and K I = 0.05.
Figure 16. Frequency deviation for LFC−1: Δ f 1 was measured while using K P = 0.05 and K I = 0.05.
Energies 16 02860 g016
Figure 17. Frequency deviation for LFC−2: Δ f 2 was measured while using K P = 0.05 and K I = 0.05.
Figure 17. Frequency deviation for LFC−2: Δ f 2 was measured while using K P = 0.05 and K I = 0.05.
Energies 16 02860 g017
Figure 18. The spectral radius, ρ , against the frequency, ω , for K P = 0.05 and K I = 0.05.
Figure 18. The spectral radius, ρ , against the frequency, ω , for K P = 0.05 and K I = 0.05.
Energies 16 02860 g018
Table 1. The DM for various K P and K I values.
Table 1. The DM for various K P and K I values.
K P τ d , s K I
Method0.050.100.150.200.400.601.0
Proposed method30.93515.2009.9497.3223.2351.8510.586
[18]30.81215.0909.8427.2113.2251.8430.591
0.0[19]30.82715.178-7.2253.2751.930-
[28]30.75615.0729.8357.2103.2311.8490.586
[29]27.84813.6998.9746.6033.0021.7450.573
Proposed method31.89515.68010.2697.5613.3541.9300.631
[18]31.77215.57010.1627.4503.3451.9220.638
0.05[19]31.76315.587-7.5093.3992.008-
[28]31.70415.54710.1527.4483.3501.9280.631
[29]27.83014.0209.2056.7773.0951.8100.616
Proposed method32.77216.11710.5607.7803.4622.0000.669
[18]32.64716.00810.4537.6693.4531.9930.676
0.10[19]32.63216.021-7.7003.5072.079-
[28]31.08315.96810.4407.6643.4571.9980.669
[29]27.00113.6509.1666.8813.1741.8630.649
Proposed method34.24816.85211.0508.1463.6412.1130.716
[18]34.12216.74410.9438.0353.6312.1060.725
0.20[19]34.156316.768-8.0583.6942.1975-
[28]28.57915.10210.4957.9983.6342.1100.716
[29]25.09012.7028.5726.4973.2091.9310.692
Proposed method35.84517.64711.5748.5363.8122.1890.662
[18]35.72817.54211.4698.4243.8022.1840.684
0.40[19]35.722317.566-8.46733.8762.2997-
[28]22.84112.1968.6096.7813.7782.1840.662
[29]20.27810.3647.0145.3382.7351.7310.637
Proposed method34.91417.16211.2398.2753.5971.8740.454
[18]34.80917.06811.1368.1553.5881.8810.480
0.60[19]34.839317.103-8.21133.7102.1141-
[28]16.2548.8396.3875.1343.0891.8640.454
[29]14.2287.3324.9443.7681.9201.1980.443
Proposed method0.5550.5460.5370.5260.4820.4340.339
[18]0.5100.4980.4850.4720.4160.3570.243
1.0[19]-------
[28]0.4860.4740.4620.4500.3960.3390.229
[29]0.4650.4550.4440.4330.3840.3320.227
- Data are not available from the published paper.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ibrahim, M.H.; Peng, A.S.; Dani, M.N.; Khalil, A.; Law, K.H.; Yunus, S.; Rahman, M.I.; Au, T.W. A Novel Computation of Delay Margin Based on Grey Wolf Optimisation for a Load Frequency Control of Two-Area-Network Power Systems. Energies 2023, 16, 2860. https://doi.org/10.3390/en16062860

AMA Style

Ibrahim MH, Peng AS, Dani MN, Khalil A, Law KH, Yunus S, Rahman MI, Au TW. A Novel Computation of Delay Margin Based on Grey Wolf Optimisation for a Load Frequency Control of Two-Area-Network Power Systems. Energies. 2023; 16(6):2860. https://doi.org/10.3390/en16062860

Chicago/Turabian Style

Ibrahim, Mohammad Haziq, Ang Swee Peng, Muhammad Norfauzi Dani, Ashraf Khalil, Kah Haw Law, Sharina Yunus, Mohammad Ishlah Rahman, and Thien Wan Au. 2023. "A Novel Computation of Delay Margin Based on Grey Wolf Optimisation for a Load Frequency Control of Two-Area-Network Power Systems" Energies 16, no. 6: 2860. https://doi.org/10.3390/en16062860

APA Style

Ibrahim, M. H., Peng, A. S., Dani, M. N., Khalil, A., Law, K. H., Yunus, S., Rahman, M. I., & Au, T. W. (2023). A Novel Computation of Delay Margin Based on Grey Wolf Optimisation for a Load Frequency Control of Two-Area-Network Power Systems. Energies, 16(6), 2860. https://doi.org/10.3390/en16062860

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