Load Frequency Control Based on the Bees Algorithm for the Great Britain Power System

This paper focuses on using the Bees Algorithm (BA) to tune the parameters of the proposed Fuzzy Proportional–Integral–Derivative with Filtered derivative (Fuzzy PIDF), Fractional Order PID (FOPID) controller and classical PID controller developed to stabilize and balance the frequency in the Great Britain (GB) power system at rated value. These controllers are proposed to meet the requirements of the GB Security and Quality of Supply Standard (GB-SQSS), which requires frequency to be brought back to its nominal value after a disturbance within a specified time. This work is extended to employ the proposed fuzzy structure controller in a dual-area interconnected power system. In comparison with controllers tuned by Particle Swarm Optimization (PSO) and Teaching Learning-Based Optimization (TLBO) used for the same systems, simulation results show that the Fuzzy PIDF tuned by BA is able to significantly reduce the deviation in the frequency and tie-line power when a sudden disturbance is applied. Furthermore, the applied controllers tuned by BA including the Fuzzy PIDF prove their high robustness against a wide range of system parametric uncertainties and different load disturbances.


Introduction
A secure supply of power energy is the main concern in power system operation, i.e., reliable electricity with appropriate quality is supplied to the customers at all times. Thus, it is essential that energy production is continuously balanced with demand: this guarantees that the system's frequency is maintained within strict limits and centered at nominal value, through load frequency control (LFC) [1]. In the Great Britain power system, over 30% of the energy produced is based on a renewable energy source (RES), such as wind turbine generators (WTG) [2]. Over the last half-decade, capacity from the RES has increased by a factor of three, with a corresponding 33% drop in the contribution of fossil fuel, with power stations being closed due to a loss of economic viability or completion of their lifespan. RES is set to remain strategically significant to GB's power system, with predictions that it will contribute around 50% of energy production by 2025. The rationale for this growth is to achieve reductions in carbon emitted [3]. The increased use of energy from RESs decreases the system inertia in power systems, and this reduced inertia presents challenges for keeping frequency within a permissible limit. The reduced inertia can also cause major increases in the absolute value of the rate of change of frequency (RoCoF), with the risk that the system may be destabilized when subjected to suddenly occurring disturbances-for example, rising or dropping production or demand. This leads to the need for fast frequency response technologies that can address the challenges that lower inertia presents [4]. While providing a fast and stable response is required by the control system to the high RoCoF, it should be considered that system oscillation can result from Table 1. Frequency limits for the GB power system.

±0.2
System frequency in normal operational conditions and the acceptable frequency deviation following a generation loss or connecting demand to ±300 MW.

±0.5
The maximum deviation in frequency when generation units over 300 MW and of up to 1320 MW is lost. −0.8 The maximum deviation in frequency following a generation loss over 1320 MW and up to 1800 MW, requiring restoration of frequency to a minimum of 49.5 Hz in 60 s.
the control system to the high RoCoF, it should be considered that system oscilla result from very fast responses. In order to control the frequency in power syste ous control loops are involved: primary, secondary, tertiary and emergency in conditions. In the GB power system, the primary loop, which is also recognize namic power generation, reaches the maximum in ten seconds, where the second trol takes thirty seconds to reach its maximum operative capacity. Frequency services are classified into dynamic services, which respond automatically to quency alternation, and non-dynamic, which are triggered via load frequenc When the generation power and the demand is significantly imbalanced, the sys experience major frequency deviations; frequency fluctuation is required by GB remain within an acceptable range, as demonstrated in Table 1 [5]. Figure 1 sh actable frequency deviations in the normal operating conditions and when a ge loss of up to 1800MW (representing a substantial nuclear generator) or less sudd curs [1,6].

±0.2
System frequency in normal operational conditio the acceptable frequency deviation following a ge loss or connecting demand to ± 300 MW.
±0. 5 The maximum deviation in frequency when gen units over 300 MW and of up to 1320 MW is l −0.8 The maximum deviation in frequency following a tion loss over 1320 MW and up to 1800 MW, requi toration of frequency to a minimum of 49.5 Hz i Recently, a number of strategic approaches have been suggested for LFC i systems: this is to keep the frequency and tie-line power flow of the system with missible range even when the system experiences sudden perturbation. The G system in particular has received considerable attention in tackling the problem quency deviation, particularly with the increase in RES penetration and its impa total system inertia [7][8][9][10].
For LFC, recent studies have applied decentralized PID controllers based on optimization techniques. Genetic algorithm (GA) is used to tune the PID controll plied in interconnected power networks [11]. Meanwhile, [12] assesses Particle Sw timization (PSO) for PID applied for LFC. Work in [13] shows that the Teaching L Based Optimization (TLBO) algorithm could be used to tune I and PID controllers Recently, a number of strategic approaches have been suggested for LFC in power systems: this is to keep the frequency and tie-line power flow of the system within a permissible range even when the system experiences sudden perturbation. The GB power system in particular has received considerable attention in tackling the problem of frequency deviation, particularly with the increase in RES penetration and its impact on the total system inertia [7][8][9][10].
For LFC, recent studies have applied decentralized PID controllers based on different optimization techniques. Genetic algorithm (GA) is used to tune the PID controller as applied in interconnected power networks [11]. Meanwhile, [12] assesses Particle Swart Optimization (PSO) for PID applied for LFC. Work in [13] shows that the Teaching Learning-Based Optimization (TLBO) algorithm could be used to tune I and PID controllers applied to a single-and multi-area multi-source with or without an HVDC link. This algorithm performed more optimally than the Differential Evolution algorithm (DE)-based PID controller. A Fractional Order Proportional Integral (FOPID or PI λ D µ ) controller provides an extra two supplementary degrees of flexibility in comparison to the traditional PID controller, in addition to the main three parameters [14]. Researchers in [15] examined the potential capabilities of this controller for LFC in a two-area hydro-thermal power system. Additionally, the Big Bang Big Crunch algorithm was applied in [16] to optimally find the gains of an FOPID controller's parameters applied for a deregulated multi-area automatic generation control (AGC) scheme.
Researchers have also widely investigated fuzzy logic controllers to control the frequency in power systems. Researchers in [17] designed and implemented a fuzzy controller as a part of an AGC system for the National Control Centre of Eskom. In [18], the Imperialist Competitive Algorithm was applied in optimizing a fuzzy proportional-integral controller's output scaling factor, using the criterion of the Integral of Squared Error (ISE) to enhance the performance of LFC in a dual-area power system. Fuzzy fractional PI λ D with filtered derivative mode is investigated in [19], and the Cuckoo optimization algorithm is used to find the optimal values of the gains of the suggested controller. The authors of [20] proposed a novel approach to enhance the frequency performance of a hybrid dual-area power system via coordination between an optimized fuzzy fine-tuning scheme and Gate-Controlled Series Capacitors (GCSC). In [21], a design involving a robust fuzzy logic-based fine-tuning approach used to tune the parameters of the classical proportionalintegral controller is proposed for LFC in a multi-area power system. Further, a two-stage robust-intelligent controller design is proposed for LFC in a single-area power system: this design is based on the Kharitonov theorem and fuzzy logic to enhance the performance of a PI controller [22].
The research works discussed above aimed to determine optimized controllers' parameter values for speeding up power provision as well as reducing frequency deviations to enhance the stability of power systems. A state-of-the-art review of the recent techniques based on different theories proposed for LFC in power systems is provided in [23].
Accordingly, as demonstrated in the literature, in recent years, there has been considerable progress in intelligent algorithm-based controllers attempting to better control LFC systems, which have solved the problem to a great extent. However, to the best of the authors' knowledge, no attempt has been made to utilize the BA in designing the secondary frequency control of a power system. Since this algorithm has demonstrated a high level of superiority and effectiveness as an optimization tool in many different fields [24][25][26][27], this promising achievement motivated the authors to make use of this very powerful optimization technique to determine the optimum parameters' gains of the PID, FOPID and Fuzzy PIDF for LFC in power systems. The main investigations of the present work are:

•
To propose a metaheuristic algorithm, the Bees Algorithm (BA), inspired by the natural behavior of honeybees, for the LFC of the GB power system.

•
To optimize PID and FPID controllers' gains and study their dynamic performance for the GB power system.

•
To design and optimize the fuzzy logic controller (FLC) scale factor gains and study its dynamic performance for GB power system.

•
To compare the dynamic performance of BA-based PID, FPID and FLC controllers with those of PSO and TLBO-tuned different controllers for the same system.

•
To investigate the effects of parametric uncertainties of the system with different load disturbances when the proposed controllers are implemented for LFC.

•
This study is then extended to examine the validation of the proposed Fuzzy PIDF in a two-area interconnected power system; the robustness analysis of Fuzzy PIDF against parametric uncertainties in this system is also validated.
It may be worth mentioning that the reason behind choosing PSO and TLBO for comparison with BA is due to their wide use in the area of LFC and their superior performance as a tool of optimization. Therefore, if the proposed algorithm provides a frequency response similar to or better than these two algorithms, this will be another successful use of BA in engineering applications.
The Artificial Bee Colony (ABC) optimization algorithm has also different applications for LFC [28,29]. Despite some similarities in the mechanism of the proposed Bees Algorithm (BA) and the ABC algorithm, the proposed BA has outperformed the ABC algorithm as an optimization tool [30]. However, a new study to further investigate the differences and the superiority between these algorithms for LFC applications may be considered for future work in the authors' research.

The Simplified GB Power System Model
A simplified model of the GB power system shown in Figure 2 was developed using MATLAB Simulink; this model is utilized to analyze the power system frequency and then design an appropriate controller. The characteristics of the generators employed in the system are considered in this simplified model, as well as damping from the loads depending on frequency. Within this, synchronous coal-powered, gas-powered, hydropower and nuclear plants are responsive to any decline in frequency and increase their generated output power correspondingly. In the model, such synchronous generators are represented by first-order transfer function blocks used to model the governor and the turbine. The governor droop gain R represents the turbine velocity control; this gain is the combined value of all droops of generator speed governors in the system. T g is the typical time constant of the governor. Stable performance of the speed control is guaranteed by introducing transient droop compensation represented as lead-lag between the governor and turbine. The output mechanical power following the response of the governor which defines the turbine model is characterized in this simulation by the time constant T t . The system inertia in this simplified model is represented by the time constant H eq , which was considered to represent the current scenario of the GB power system with the high penetration of RERs. The damping obtained from frequency-dependent loads is represented by an equivalent gain value D. The effect of an electrical vehicle (EV) was also considered in this design, which was modeled as an aggregated value represented by a feedback gain in the primary loop, with an estimated aggregated load equal to 2.16 GW. The value of this load was considered to be =1.35 pu (EV load × f/(Network base)), where the network base value is equal to 79.2 GW [2]; the effect of EV gain on the primary loop response is provided in Appendix A. The secondary control applied in this model is the main study of this paper, which will be examined via different controllers tuned by proposed algorithms including the Bees Algorithm, which represents the main contribution of this work. The parameters applied in this simplified model are tabulated in Table 2 [2,6,31]. It should be noted that the model used in this study to examine the frequency response of the GB power system is simplified; accordingly, the small effect of nonlinearities such as governor dead band (GDB) and generation rate constrain (GRC) is neglected.

Classical Controllers
PID controllers are the most widely applied feedback controller for the process sector. Easy to understand and offering robustness, strong performance in spite of the differing and dynamic properties of the process plant and the low cost are the main merits of this controller. The PID controller essentially includes proportional, integral and derivative modes. While a proportional controller decreases the rise time, it cannot completely remove the steady-state error. Integral controls, on the other hand, eliminate steady-state error, yet they can worsen transient responses. Derivative controls improve the system stability, decrease the overshoot/undershoot and enhance transient responses. Current industrial applications most frequently rely upon the proportional-integral (PI) controller. Controls without derivative action are applied in the following circumstances: when it is not important for the system to respond rapidly; when significant noise and disturbance are experienced during the system operation. The overall stability of the system can be improved by adding the derivative mode as it enables an increase in the proportional gain and decrease in the integral, thus increasing the speed response of the controller. Accordingly, PID controllers are frequently applied in systems that require stability with a fast response. Considering these points, this study investigates the effectiveness of a PID controller for the LFC of the GB power system. Equation (1) illustrates the transfer function of this controller, where K P is the proportional gain, K I the integral gain, and K D is the derivative gain.
Generally, for LFC, conventional methods of control remain the main approach, but in more complex systems, such approaches can become inadequate. A recently introduced control approach used for LFC tasks is known as non-integer control, developed on the basis of fractional calculus, which is generalized from integer order calculus. Fractional calculus is a generalization of differentiation/integration to a non-integer order, and thus provides n degrees of additional freedom when designing controllers, which can make them more efficient, flexible and robust [14]. The FOPID controller's transfer function is shown in Equation (2), where λ is the order of integration, and μ is the order of the differentiator.
As highlighted in the literature, fractional-order control (FOC) has emerged to address the problem of LFC. However, this is the first attempt to tune the parameters of this controller using the Bees Algorithm.

Classical Controllers
PID controllers are the most widely applied feedback controller for the process sector. Easy to understand and offering robustness, strong performance in spite of the differing and dynamic properties of the process plant and the low cost are the main merits of this controller. The PID controller essentially includes proportional, integral and derivative modes. While a proportional controller decreases the rise time, it cannot completely remove the steady-state error. Integral controls, on the other hand, eliminate steady-state error, yet they can worsen transient responses. Derivative controls improve the system stability, decrease the overshoot/undershoot and enhance transient responses. Current industrial applications most frequently rely upon the proportional-integral (PI) controller. Controls without derivative action are applied in the following circumstances: when it is not important for the system to respond rapidly; when significant noise and disturbance are experienced during the system operation. The overall stability of the system can be improved by adding the derivative mode as it enables an increase in the proportional gain and decrease in the integral, thus increasing the speed response of the controller. Accordingly, PID controllers are frequently applied in systems that require stability with a fast response. Considering these points, this study investigates the effectiveness of a PID controller for the LFC of the GB power system. Equation (1) illustrates the transfer function of this controller, where K p is the proportional gain, K I the integral gain, and K D is the derivative gain.
Generally, for LFC, conventional methods of control remain the main approach, but in more complex systems, such approaches can become inadequate. A recently introduced control approach used for LFC tasks is known as non-integer control, developed on the basis of fractional calculus, which is generalized from integer order calculus. Fractional calculus is a generalization of differentiation/integration to a non-integer order, and thus provides n degrees of additional freedom when designing controllers, which can make them more efficient, flexible and robust [14]. The FOPID controller's transfer function is shown in Equation (2), where λ is the order of integration, and µ is the order of the differentiator.
As highlighted in the literature, fractional-order control (FOC) has emerged to address the problem of LFC. However, this is the first attempt to tune the parameters of this controller using the Bees Algorithm.

Fuzzy PID Logic Control
There has been a broad implementation of fuzzy PID controllers in different structures to solve the LFC problem and a significant enhancement in performance has been achieved, as reported in the literature. It is also proven that in order to gain further enhancement in the overall system performance and improve the stability, a simple filter for the derivative mode of the fuzzy PID controller can be applied [32,33]. Moreover, the performance of these controllers mainly relies on the selection of the scaling factor gains of the input and output of the controller, but it is difficult to find the optimum values of these gains using the trial-and-error technique. In view of the above, a fuzzy PID controller with derivative filter (Fuzzy PIDF), in which the scaling factor gains are tuned by the Bees Algorithm, is proposed in this section for LFC. The structural design of this controller is illustrated in Figure 3. As is clear from the figure, the controller has two inputs, ∆F and the derivative of ∆F and one output; in the case of employing this controller in a multi-area power system, the two inputs are the Area Control Error (ACE) and the derivative of (ACE).

Fuzzy PID Logic Control
There has been a broad implementation of fuzzy PID controllers in different structures to solve the LFC problem and a significant enhancement in performance has been achieved, as reported in the literature. It is also proven that in order to gain further enhancement in the overall system performance and improve the stability, a simple filter for the derivative mode of the fuzzy PID controller can be applied [32,33]. Moreover, the performance of these controllers mainly relies on the selection of the scaling factor gains of the input and output of the controller, but it is difficult to find the optimum values of these gains using the trial-and-error technique. In view of the above, a fuzzy PID controller with derivative filter (Fuzzy PIDF), in which the scaling factor gains are tuned by the Bees Algorithm, is proposed in this section for LFC. The structural design of this controller is illustrated in Figure 3. As is clear from the figure, the controller has two inputs, ∆F and the derivative of ∆F and one output; in the case of employing this controller in a multiarea power system, the two inputs are the Area Control Error (ACE) and the derivative of (ACE). The scaling factor gains of the input are (K 1 and K 2 ) and four scaling factors in the output, namely K P , K I , K D , and K F is the filter gain. Due to its simplicity and the lower computation time needed for this type of membership, five triangular membership functions are used for the inputs and the output variables shown in Figure 4, namely Negative Big, Negative Small, Zero, Positive Small and Positive Big. Thus, 25 rules are required to generate the fuzzy output of the controller. Table 3 depicts the rule base of the proposed controller. Since the performance of the controller depends on these rules, the tabulated rules are generated by a comprehensive study of the dynamic behaviors of the testbed power system. The "Mamdani" interface tool is used in this controller for the fuzzification stage, while the "Centroid" method is used in the defuzzification stage to convert the fuzzy output value of the controller to a real value.  The scaling factor gains of the input are (K 1 and K 2 ) and four scaling factors in the output, namely K p , K I , K D and K F is the filter gain. Due to its simplicity and the lower computation time needed for this type of membership, five triangular membership functions are used for the inputs and the output variables shown in Figure 4, namely Negative Big, Negative Small, Zero, Positive Small and Positive Big. Thus, 25 rules are required to generate the fuzzy output of the controller. Table 3 depicts the rule base of the proposed controller. Since the performance of the controller depends on these rules, the tabulated rules are generated by a comprehensive study of the dynamic behaviors of the testbed power system. The "Mamdani" interface tool is used in this controller for the fuzzification stage, while the "Centroid" method is used in the defuzzification stage to convert the fuzzy output value of the controller to a real value.

Fuzzy PID Logic Control
There has been a broad implementation of fuzzy PID controllers in different structures to solve the LFC problem and a significant enhancement in performance has been achieved, as reported in the literature. It is also proven that in order to gain further enhancement in the overall system performance and improve the stability, a simple filter for the derivative mode of the fuzzy PID controller can be applied [32,33]. Moreover, the performance of these controllers mainly relies on the selection of the scaling factor gains of the input and output of the controller, but it is difficult to find the optimum values of these gains using the trial-and-error technique. In view of the above, a fuzzy PID controller with derivative filter (Fuzzy PIDF), in which the scaling factor gains are tuned by the Bees Algorithm, is proposed in this section for LFC. The structural design of this controller is illustrated in Figure 3. As is clear from the figure, the controller has two inputs, ∆F and the derivative of ∆F and one output; in the case of employing this controller in a multiarea power system, the two inputs are the Area Control Error (ACE) and the derivative of (ACE). The scaling factor gains of the input are (K 1 and K 2 ) and four scaling factors in the output, namely K P , K I , K D , and K F is the filter gain. Due to its simplicity and the lower computation time needed for this type of membership, five triangular membership functions are used for the inputs and the output variables shown in Figure 4, namely Negative Big, Negative Small, Zero, Positive Small and Positive Big. Thus, 25 rules are required to generate the fuzzy output of the controller. Table 3 depicts the rule base of the proposed controller. Since the performance of the controller depends on these rules, the tabulated rules are generated by a comprehensive study of the dynamic behaviors of the testbed power system. The "Mamdani" interface tool is used in this controller for the fuzzification stage, while the "Centroid" method is used in the defuzzification stage to convert the fuzzy output value of the controller to a real value.

Objective Functions
In the design of modern controllers and for any controlled system, stable performance and fast response are required. However, in practice, both requirements are never achievable simultaneously. Therefore, a compromise between quick response and excellent stability is considered when designing a controller, which is achievable by adequately selecting an appropriate controller and designing it by minimizing a properly selected cost/objective function with the aid of an optimization algorithm. The objective function used to tune the controller mainly relies on a performance criterion that considers the overall closed-loop response of the system. Many objective functions have been proposed in the control design, in which four kinds are most often used for LFC. Because of their better performance compared to the other criteria [13], the Integral of Square Error (ISE) illustrated in Equation (3) and Integral Time Absolute Error (ITAE) expressed in Equation (4) are used in this paper. Therefore, PID, FOPID and Fuzzy PIDF controllers are designed for the LFC of the GB power system by minimizing the defined objective functions with the help of the Bees Algorithm and other two techniques.
It is proven that with ISE, large errors are more penalized than smaller ones. Thus, control systems designed to minimize ISE are more most likely to eliminate large errors quickly. However, they have to tolerate small ones that are continuous for a long period of time. ITAE calculates the integration of the absolute error multiplied by the time over the simulation period. This criterion is based on weighing errors that occur after a long time more largely than those that exist during the beginning of the response [34]. Control systems specified based on ITAE tuning tend to settle much more quickly than the ISE tuning methods.

The Proposed Algorithm
Over the last two decades, considerable interest has been shown by many researchers to propose algorithms inspired by the lifestyle of different animals and insects as well as the behavior of natural phenomena. Most of these algorithms have successfully demonstrated their potential in many applications to solve different optimization problems. This section precisely summarizes the main concept in steps that explain the Bees Algorithm (BA). Since its invention by Pham et al. in 2005 [24], several studies have been introduced based on this algorithm and different problems in many fields have been successfully solved using the BA [24][25][26][27]. The simplest pseudo-code for this algorithm is shown in Figure 5. As with any optimization algorithm, a number of parameters are required to be setspecifically, the number of scout bees (n), number of selected sites for search out of n visited patches (m), number of best (elite) sites among them selected sites (e), number of recruited bees in the best-selected e sites (nep), number of bees sent to the other (m-e) selected sites (nsp), the initial size of each patch (ngh), which includes the site and the neighborhood, area as well as the stopping criteria. The mechanism of this algorithm begins with placing the n scout bees randomly in the search space. In step 2, the evaluation of the fitness of sites visited by the n scout bees is done. The m sites with the highest fitness in specified "chosen sites" in step 3 are selected for neighborhood or local search. The algorithm in steps 4 and 5 conducts searches in the neighborhood of the selected sites, with more bees assigned to the best e sites. Selection of the best sites could be conducted directly based on the fitness associated with them. Alternatively, using the fitness values, the probability of sites being selected is determined. Searches in the neighborhood of the best e sites which represent the most promising solutions are made more prominent by recruiting more bees for them than for the other selected sites. Together with scouting, this differential recruitment is a key operation of the Bees Algorithm. For each patch, only the one bee that has found the site with the highest fitness (the "fittest" bee) is selected in step 5 to form part of the next bee population. In steps 6-8, the remaining bees in the population n are assigned randomly around the search space to scout for potential new solutions or to conduct the global search. These eight steps are repeated until a stopping criterion is met. The colony will have two parts to its new population at the end of each iteration: representatives from each selected patch and other scout bees assigned to conduct random searches [35].

Results and Discussion
This work was implemented in MATLAB (2019a), the BA, TLBO and PSO codes were programmed in (.m files), and the model of the GB power system was developed in the MATLAB Simulink environment. The parameters of BA and PSO were set as depicted in Tables 4 and 5, respectively. With TLBO, the population size was set to 50, and the maximum number of iterations was taken as 40 for all algorithms. As with any optimization algorithm, a number of parameters are required to be setspecifically, the number of scout bees (n), number of selected sites for search out of n visited patches (m), number of best (elite) sites among them selected sites (e), number of recruited bees in the best-selected e sites (nep), number of bees sent to the other (m-e) selected sites (nsp), the initial size of each patch (ngh), which includes the site and the neighborhood, area as well as the stopping criteria. The mechanism of this algorithm begins with placing the n scout bees randomly in the search space. In step 2, the evaluation of the fitness of sites visited by the n scout bees is done. The m sites with the highest fitness in specified "chosen sites" in step 3 are selected for neighborhood or local search. The algorithm in steps 4 and 5 conducts searches in the neighborhood of the selected sites, with more bees assigned to the best e sites. Selection of the best sites could be conducted directly based on the fitness associated with them. Alternatively, using the fitness values, the probability of sites being selected is determined. Searches in the neighborhood of the best e sites which represent the most promising solutions are made more prominent by recruiting more bees for them than for the other selected sites. Together with scouting, this differential recruitment is a key operation of the Bees Algorithm. For each patch, only the one bee that has found the site with the highest fitness (the "fittest" bee) is selected in step 5 to form part of the next bee population. In steps 6-8, the remaining bees in the population n are assigned randomly around the search space to scout for potential new solutions or to conduct the global search. These eight steps are repeated until a stopping criterion is met. The colony will have two parts to its new population at the end of each iteration: representatives from each selected patch and other scout bees assigned to conduct random searches [35].

Results and Discussion
This work was implemented in MATLAB (2019a), the BA, TLBO and PSO codes were programmed in (.m files), and the model of the GB power system was developed in the MATLAB Simulink environment. The parameters of BA and PSO were set as depicted in Tables 4 and 5, respectively. With TLBO, the population size was set to 50, and the maximum number of iterations was taken as 40 for all algorithms.  To study the dynamic performance of the GB power system, a step load perturbation of 0.03955 pu (at t = 5 s) represents a loss of generation unit equal to 1.32 GW (two of large generators, 660 MW) of the total generation power of the GB system used, which occurred in the GB system on 27th May 2008 [2].
ITAE and ISE are taken separately as objective functions to tune the parameters of the proposed controllers using the above-mentioned algorithms for LFC in the generalized GB power system model. Initially, the parameters of the PID controller are optimized; it is found that the PID tuned by the proposed algorithms performs satisfactorily to damp out the drop in the frequency. However, a reduced drop in frequency with a slow response is achieved when ISE is considered as an objective function, while with ITAE, the drop in frequency worsened with the fast response obtained, bringing the frequency back to the nominal value in a shorter period of time. Then, the gains of FOPID are tuned; in this regard, it is worth highlighting that, in general, the tuned FOPID tuned by BA designed via minimizing ITAE provides better results compared to the tuned PID. However, FOPID designed with ISE fails to bring the frequency back to its nominal value. Thereafter, Fuzzy PIDF parameters are tuned, where a significant improvement is achieved in comparison with PID and FOPID.

Classical Controllers
The gains of the conventional PID and FOPID controllers obtained using BA, TLBO and PSO optimization algorithms using the suggested objective functions are depicted in Table 6. Results obtained based on the proposed algorithms "BA" tuning PID and FOPID are compared with those of TLBO and PSO; it is found to be a successful solution and provides excellent performance in many aspects, i.e., undershoot, overshoot and settling time, which are shown in Tables 7-10. The changes in the frequency of the GB power system when different conventional controllers tuned by different optimization algorithms are used to control the frequency deviation are shown in Figures 6-9. Results obtained based on the proposed algorithms "BA" tuning PID and FOPID are compared with those of TLBO and PSO; it is found to be a successful solution and provides excellent performance in many aspects, i.e., undershoot, overshoot and settling time, which are shown in Tables 7-10. The changes in the frequency of the GB power system when different conventional controllers tuned by different optimization algorithms are used to control the frequency deviation are shown in Figures 6-9.
From Figure 6 and Table 7, it is clear that when the PID is tuned using the proposed algorithms based on ISE, almost identical responses are obtained. On the other hand, when ITAE is used as an objective function to design the PID controller, as demonstrated in Figure 7 and Table 8, BA shows better performance in terms of undershoot and overshoot, with only a 0.1840 Hz drop in frequency.        As shown in Figure 8 and Table 9, FOPID designed by minimizing ISE is found to be less effective in eliminating the steady-state error, which made this technique a less preferable option for this system. In terms of undershoot and settling time T , Figure 9 and Table 10 prove that FOPID tuned by the proposed BA using ITAE as an objective function provides the best performance in comparison with the other classical controllers tuned by TLBO and PSO.

Fuzzy PIDF Controller
The optimal gains of the proposed Fuzzy PID with derivative filter obtained by the proposed BA, TLBO and PSO algorithms using the suggested objective functions are depicted in Table 11. Simulation results obtained with the BA are compared with those of TLBO and PSO and it is found to be an excellent tool and provides improved performance in many aspects.  From Figure 6 and Table 7, it is clear that when the PID is tuned using the proposed algorithms based on ISE, almost identical responses are obtained. On the other hand, when ITAE is used as an objective function to design the PID controller, as demonstrated in Figure 7 and Table 8, BA shows better performance in terms of undershoot and overshoot, with only a 0.1840 Hz drop in frequency.
As shown in Figure 8 and Table 9, FOPID designed by minimizing ISE is found to be less effective in eliminating the steady-state error, which made this technique a less preferable option for this system. In terms of undershoot and settling time T s , Figure 9 and Table 10 prove that FOPID tuned by the proposed BA using ITAE as an objective function provides the best performance in comparison with the other classical controllers tuned by TLBO and PSO.

Fuzzy PIDF Controller
The optimal gains of the proposed Fuzzy PID with derivative filter obtained by the proposed BA, TLBO and PSO algorithms using the suggested objective functions are depicted in Table 11. Simulation results obtained with the BA are compared with those of TLBO and PSO and it is found to be an excellent tool and provides improved performance in many aspects. Tables 12 and 13 and Figures 10 and 11 demonstrate the frequency response of the GB power system when the Fuzzy PIDF is optimized by different optimization algorithms implemented for LFC. It is observed that a significant improvement is achieved in comparison with classical controllers.  Tables 12 and 13 and Figures 10 and 11 demonstrate the frequency response of the GB power system when the Fuzzy PIDF is optimized by different optimization algorithms implemented for LFC. It is observed that a significant improvement is achieved in comparison with classical controllers.
Moreover, BA has proved to be a powerful technique to tune the Fuzzy PIDF as the results obtained from the proposed controller tuned by BA prove that the performance of the system experiences a clear improvement in terms of undershoot and settling time. Regarding the error and overshoot, almost similar results are obtained based on all algorithms.     Moreover, BA has proved to be a powerful technique to tune the Fuzzy PIDF as the results obtained from the proposed controller tuned by BA prove that the performance of the system experiences a clear improvement in terms of undershoot and settling time. Regarding the error and overshoot, almost similar results are obtained based on all algorithms.

Robustness Analysis against System Uncertainty
Parameters within the system, including the damping coefficient D, speed regulator R, system inertia coefficient H eq and turbine governor time constant T g , are subject to constant fluctuation, which can lead to a significant degradation in the performance of close-loop systems. There has been comparatively less focus in research on this issue within load frequency control; for example, the increase in the total system inertia will slow down the system response, while the frequency deviation decreases if the damping ratio increases, and if the governor time constant increases, the frequency deviation will increase. The impact of the variation in each parameter on the frequency response of the GB power system is provided in Appendix B. Therefore, investigations are carried out in order to study the consequences of parametric uncertainties in the system. For this, each parameter in the system is altered by ±50% from its nominal value. Two different scenarios of parameters' uncertainty, T g , D, R and H eq (listed in Table 14 and shown in Figure 12), are considered for the simplified GB power system model examination. In this sub-section, only controllers tuned by the proposed BA are examined. The optimal gains obtained during the normal condition will not be re-tuned when the model is subjected to variation in system parameters.

Robustness Analysis against System Uncertainty
Parameters within the system, including the damping coefficient D, speed regulator R, system inertia coefficient H eq and turbine governor time constant T g , are subject to constant fluctuation, which can lead to a significant degradation in the performance of close-loop systems. There has been comparatively less focus in research on this issue within load frequency control; for example, the increase in the total system inertia will slow down the system response, while the frequency deviation decreases if the damping ratio increases, and if the governor time constant increases, the frequency deviation will increase. The impact of the variation in each parameter on the frequency response of the GB power system is provided in Appendix B. Therefore, investigations are carried out in order to study the consequences of parametric uncertainties in the system. For this, each parameter in the system is altered by ±50% from its nominal value. Two different scenarios of parameters' uncertainty, T g , D, R and H eq (listed in Table 14 and shown in Figure 12), are considered for the simplified GB power system model examination. In this sub-section, only controllers tuned by the proposed BA are examined. The optimal gains obtained during the normal condition will not be re-tuned when the model is subjected to variation in system parameters.  Furthermore, as shown in Figure 12, and since, in the second scenario, the frequency response of the system is worse than the nominal case, the second scenario only is investigated, which also represents a possible decline in the total system inertia of the GB power system due to the increasing use of renewable energy resources. Figures 13 and 14 show the frequency response of the GB power system under parametric uncertainties when different controllers tuned by BA are employed as the LFC system. From Figures 13 and 14, it is noted that Fuzzy PIDF controllers provide high stability while classical controllers show less robustness against system uncertainty, with the worst drop in frequency recorded at −0.178 Hz when PID is applied for scenario 2. The frequency responses of the system with different BA-tuned controllers based on ISE and ITAE, respectively, are listed in Tables 15 and 16. for scenario 2.

Different Load Disturbances
To further investigate the robustness of Fuzzy PIDF, FOPID and PID tuned by BA, a loss of 1.8 GW (very large nuclear generator) in the generation unit representing around 0.053 pu is considered in this sub-section. The frequency response of the GB power system with the new load disturbance is shown in Figure 15 and the frequency response performances are tabulated in Table 17. Furthermore, in order to further examine the robustness of the proposed techniques, parameter uncertainties from scenario 2 are considered when a power generation of 0.053pu is lost and the frequency response of the system in this case is given in Figure 16; the frequency response performances are depicted in Table 18, from which it is obvious that the proposed controller "Fuzzy PIDF" tuned by BA is robust and performs satisfactorily even when a larger generator is lost with parameter uncertainties. Note that only controllers tuned by BA based on minimizing the ITAE objective function are considered in this part.

Different Load Disturbances
To further investigate the robustness of Fuzzy PIDF, FOPID and PID tuned by BA, a loss of 1.8 GW (very large nuclear generator) in the generation unit representing around 0.053 pu is considered in this sub-section. The frequency response of the GB power system with the new load disturbance is shown in Figure 15 and the frequency response performances are tabulated in Table 17. Furthermore, in order to further examine the robustness of the proposed techniques, parameter uncertainties from scenario 2 are considered when a power generation of 0.053pu is lost and the frequency response of the system in this case is given in Figure 16; the frequency response performances are depicted in Table 18, from which it is obvious that the proposed controller "Fuzzy PIDF" tuned by BA is robust and performs satisfactorily even when a larger generator is lost with parameter uncertainties. Note that only controllers tuned by BA based on minimizing the ITAE objective function are considered in this part.

Load Frequency Control for Dual-Area Power System
This study is extended to a multi-area interconnected power system as demonstrated in Figure 17. This model is extensively investigated in the literature to study the dynamic behavior of different control concepts for LFC in power systems [36,37]. It consists of two unequal thermal interconnected areas; the parameters of this model are tabulated in Appendix C. Based on the function of the LFC loop in multi-area power systems, the main tasks of LFC are providing the required power from the generation units to meet the load variation and maintaining the interchanged power among interconnected control areas via tie-lines at pre-rated values. The specified objectives of LFC that contribute to improving the overall power system stability are to guarantee zero steady-state errors in frequency fluctuations of different control areas along with the tie-line power fluctuations. It is also responsible for reducing the overshoot and undershoot of the oscillation within a specified time; this depends on the size of the disturbance and the capacity of the power system. (6).
ACE area 2 ∆P 12 + B 1 ∆F 1 (6) where ∆F 1 and ∆F 2 are the frequency deviation in areas 1 and 2, respectively, ∆P 12 and ∆P 21 are the power flow deviation in areas 1 and 2, and B , B are frequency biases.  The term defined as Area Control Error (ACE) in each area is the input of the controller located in that area. In this model, the ACEs are represented in Equations (5) and (6).
where ∆F 1 and ∆F 2 are the frequency deviation in areas 1 and 2, respectively, ∆P 12 and ∆P 21 are the power flow deviation in areas 1 and 2, and B 1 , B 2 are frequency biases. The proposed Fuzzy PID with the filtered derivative mode is applied in this section as LFC in the dual-area interconnected power system. A step load disturbance of 0.2 pu is applied in area 1 to investigate the dynamic performance of the system with the proposed controller. The optimum values of the scaling factors of Fuzzy PID with Filtered derivative (Fuzzy PIDF) optimized by BA, PSO and TLBO by minimizing the ITAE objective function expressed in Equation (7) are depicted in Table 19. Due to the sudden subjection of the system to 0.2 pu disturbance in area 1, the frequency and tie-line power flow within the system experience a range of deviation from their nominal values. Frequency deviation in area 1 (∆F 1 ), frequency deviation in area 2 (∆F 2 ) and tie-line power flow deviation (∆P tie ) are illustrated in Figures 18-20, respectively. Table 20 provides the characteristics of the dynamic performance of the testbed systembased Fuzzy PIDF tuned by different algorithms, where undershoot (U sh ), overshoot (O sh ) and settling time (T s ) of frequency in both areas and tie-line power along with the values of the objective function are given. Due to the sudden subjection of the system to 0.2 pu disturbance in area 1, the fre quency and tie-line power flow within the system experience a range of deviation from their nominal values. Frequency deviation in area 1 (∆F 1 ), frequency deviation in area (∆F 2 ) and tie-line power flow deviation (∆P tie ) are illustrated in Figures 18-20, respec tively. Table 20 provides the characteristics of the dynamic performance of the testbed system-based Fuzzy PIDF tuned by different algorithms, where undershoot (U sh ), over shoot (O sh ) and settling time (T s ) of frequency in both areas and tie-line power alon with the values of the objective function are given.       Table 20, it is observed that the proposed Fuzzy PIDF controller provides an excellent dynamic response in different aspects. Moreover, despite the small difference in results obtained by applying the proposed controller based on the suggested optimization techniques, it is obvious that the Fuzzy PIDF tuned by BA offers the best response among the investigated techniques in terms of peak undershoot experienced in the area 1, with only −0.04 Hz recorded after applying the sudden disturbance, while the smallest drop in the frequency of area 2 was −0.0036 Hz based on the TLBO-and PSO-tuned Fuzzy PIDF.
To examine the robustness of the Fuzzy PIDF controller employed for LFC in the testbed dual-area power system, the parameters of the investigated system are simultaneously varied from their nominal values as shown in Table 21. The same disturbance of 0.2 pu is applied in area 1 and the optimum values of the Fuzzy PIDF parameters obtained in the normal condition are used to verify the robustness of the proposed controller.     The simulation results obtained from robustness analysis given in Figures 21-23 and Table 22 validate the robustness of the proposed Fuzzy PIDF-controlled LFC in the dualarea power system towards a wide range of parametric uncertainties. Furthermore, it is noticed that the TLBO-optimized Fuzzy PIDF is less robust against parametric uncertainties as compared to BA and PSO tuned with the same controller.  The simulation results obtained from robustness analysis given in Figures 21-23 and Table 22 validate the robustness of the proposed Fuzzy PIDF-controlled LFC in the dualarea power system towards a wide range of parametric uncertainties. Furthermore, it is noticed that the TLBO-optimized Fuzzy PIDF is less robust against parametric uncertainties as compared to BA and PSO tuned with the same controller.

Conclusions
In this paper, the Bees Algorithm is proposed for the first time in load frequency control (LFC). This algorithm was successfully used to find the optimal parameters of PID, FOPID and the proposed Fuzzy PID with filter for the derivative action (Fuzzy PIDF) in order to solve the problem of LFC for the GB power system. A step load disturbance was applied to study the frequency performance of the system, taking ISE and ITAE as objective functions to design the above-mentioned controllers tuned by BA. The parameters of PID and FOPID were first optimized, and the obtained results were compared with those from TLBO and PSO; it was evident from the results that the BA-tuned FOPID designed by minimizing ITAE offered the best performance among the investigated classical controllers. In terms of applying the proposed Fuzzy PIDF, an obvious improvement in the performance of the system was achieved, and the obtained results from this controller based on BA, TLBO and PSO were somewhat similar, with the lowest drop in frequency equal to −0.0028 Hz when BA was used. Furthermore, it was demonstrated that the Fuzzy PIDF tuned by BA is robust against system uncertainties and different load disturbances. Finally, the research was extended to examine the validation of the proposed Fuzzy PIDF as LFC system in a dual-area interconnected power system; the obtained results emphasize that the Fuzzy PIDF-based BA is an excellent and robust controller and can be successfully implemented for LFC in different power systems.

Appendix B
The effect of parametric uncertainties on the frequency response of the GB power system.

Appendix B
The effect of parametric uncertainties on the frequency response of the GB power system.

Appendix C
The parameters of the system as used in [6,12] are as follows (T = 2, F = 60 Hz, SLP = 0.2 pu):

Appendix C
The parameters of the system as used in [6,12] are as follows (T = 2, F = 60 Hz, SLP = 0.2 pu):

Appendix C
The parameters of the system as used in [6,12]   where T is the synchronization constant, F is the system frequency, SLP is the step load perturbation, Tg is the time constant of the governor (second), Tt is the time constant of the turbine (second), R is the regulation constant (MW/Hz), B is the frequency bias constant (Hz/MW), H is the inertia time constant (second), and D is the ratio of variation in load to variation in frequency.