Next Article in Journal
Optimising High-Rise Buildings for Self-Sufficiency in Energy Consumption and Food Production Using Artificial Intelligence: Case of Europoint Complex in Rotterdam
Next Article in Special Issue
New Hybrid Invasive Weed Optimization and Machine Learning Approach for Fault Detection
Previous Article in Journal
The Place of Energy Security in the National Security Framework: An Assessment Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Different Fuzzy Control Configurations Tuned by the Bees Algorithm for LFC of Two-Area Power System †

1
Wolfson Centre for Magnetics, School of Engineering, Cardiff University, Cardiff CF24 3AA, UK
2
High Value Manufacturing Group, School of Engineering, Cardiff University, Cardiff CF24 3AA, UK
*
Author to whom correspondence should be addressed.
This paper is an extended version of our paper published in IEEE 2021 56th International Universities Power Engineering Conference (UPEC), Middlesbrough, UK, 31 August–3 September 2021; pp. 1–6.
Energies 2022, 15(2), 657; https://doi.org/10.3390/en15020657
Submission received: 7 December 2021 / Revised: 27 December 2021 / Accepted: 14 January 2022 / Published: 17 January 2022

Abstract

:
This study develops and implements a design of the Fuzzy Proportional Integral Derivative with filtered derivative mode (Fuzzy PIDF) for Load Frequency Control (LFC) of a two-area interconnected power system. To attain the optimal values of the proposed structure’s parameters which guarantee the best possible performance, the Bees Algorithm (BA) and other optimisation tools are used to accomplish this task. A Step Load Perturbation (SLP) of 0.2 pu is applied in area one to examine the dynamic performance of the system with the proposed controller employed as the LFC system. The supremacy of Fuzzy PIDF is proven by comparing the results with those of previous studies for the same power system. As the designed controller is required to provide reliable performance, this study is further extended to propose three different fuzzy control configurations that offer higher reliability, namely Fuzzy Cascade PI − PD, Fuzzy PI plus Fuzzy PD, and Fuzzy (PI + PD), optimized by the BA for the LFC for the same dual-area power system. Moreover, an extensive examination of the robustness of these structures towards the parametric uncertainties of the investigated power system, considering thirteen cases, is carried out. The simulation results indicate that the contribution of the BA tuned the proposed fuzzy control structures in alleviating the overshoot, undershoot, and the settling time of the frequency in both areas and the tie-line power oscillations. Based on the obtained results, it is revealed that the lowest drop of the frequency in area one is −0.0414 Hz, which is achieved by the proposed Fuzzy PIDF tuned by the BA. It is also divulged that the proposed techniques, as was evidenced by their performance, offer a good transient response, a considerable capability for disturbance rejection, and an insensitivity towards the parametric uncertainty of the controlled system.

1. Introduction

Modern large power systems normally comprise multiple interconnected control areas that are based on diverse energy resources. It is extensively reported that any mismatch between the generated power and customer demand results in a deviation in the frequency and tie-line power interchange [1,2,3,4]. This deviation in some cases may lead to system performance degradation and even cause the power systems to collapse [5]. Load Frequency Control (LFC) is introduced in power systems as a key service that plays its role in maintaining the frequency and tie-line power within scheduled limits in normal conditions as well as in the event of any sudden disturbance. This results in the improvement of the stability of the power systems, forming a successful operation [6,7].
In classical LFC applications, the Proportional Integral (PI) and Proportional Integral Derivative (PID) controllers are commonly utilised. The classical PID with a filtered derivative mode optimized by an improved version of the Jaya algorithm, called the Self-Adaptive Multi-Population Elitist (SAMPE) Jaya optimizer, is proposed for the LFC in a two-area interconnected power system [8]. A novel Predictive Functional Modified PID (PFMPID), tuned by the Grasshopper Optimization Algorithm (GOA), is suggested to solve the LFC problem in a three-area power system under a restructured environment and including various generation units. Simulation results demonstrated that the proposed PFMPID outperformed the classical PID [9]. A cascade of PI and PD with a filtered derivative action (PI − PDn) tuned by an Enhanced Coyote Optimization Algorithm (ECOA), equipped in a dual-area power system that consists of a photovoltaic (PV) unit to overcome the problem of frequency and tie-line power deviation, is studied in [10]. The proposed PI − PDn has shown a great superiority over the other investigated controllers. However, due to the escalating complexities with high nonlinearities in current power systems, investigations into the possibility of applying other control techniques are necessitated [11]. Accordingly, several studies have proposed many strategies based on different control theories to address the challenges of the LFC in power systems. For example, Model Predictive Control (MPC) is employed for the LFC in a two-area power system that comprises a photovoltaic generation unit [12]. Sliding Mode Control (SMC) is used for the LFC in a simplified Great Britain (GB) power system [13]. In [14], a novel adaptive sliding mode control method is proposed for the LFC in a three-equal-area interconnected power system with non-reheat turbines; this design demonstrated a better performance in comparison with the classical adaptive sliding mode control scheme. A discrete LFC method has been suggested for power systems with a high penetration of wind power, based on a sampled-data control; this technique has been examined in a single-area power system, a dual-area interconnected power system, and a three-area restricted power system [15]. A Linear Matrix Inequality (LMI)-based LFC with a communication delay in a two-area electrical system is investigated in [16]. H control is used in an isolated, distributed generation system for load frequency control [17]. Moreover, due to its merits, Fuzzy Logic Control (FLC) has recently been widely addressed as a potential solution for LFC based on different structures. It is revealed that FLC can successfully handle the problem of load frequency control. However, there was no identified rule to be utilised in order to find the fuzzy parameters, i.e., the scaling factors of the inputs and outputs as well as the membership functions and the rule base [7]. Therefore, soft computing methods have emerged to deal with this issue. In [18], a fuzzy hierarchal scheme, tuned by Particle Swarm Optimization (PSO), is proposed for the LFC in the Great Britain (GB) power system. Differential Evolution (DE) was employed to optimize the parameters of a fuzzy PID with a derivative filter for the LFC in a multi-sourced, deregulated power system [19]. An optimized fuzzy self-tuning PID controller is proposed for the LFC in two- and three-area interconnected power systems [20]. To ameliorate the proposed controller, a Tribe-DE optimization algorithm was utilized to find the optimum values of the scaling factor and the membership function parameters of the fuzzy PID controllers. The fuzzy PID tuned by Teaching Learning Based Optimization (TLBO) for a dual-area power system is studied in [21]. A novel hybrid DE and Pattern Search (PS) has been used to tune the scaling factor gains of fuzzy PI/PID controllers employed for the LFC in a two-area power system [22]. The most recent controllers, based on the different strategies employed for LFCs in power systems, are concluded in [23].
In view of the above, this work proposes a design of a Fuzzy PID with a derivative filter action (Fuzzy PIDF) employed for the LFC in an unequal two-area interconnected thermal power system. Moreover, as the classical methods, including hit and trail, used to determine the optimum values of the scaling factor gains of fuzzy PIDF are time-consuming and may not provide desirable solutions, the parameters of the proposed controller are optimized by an algorithm called the Bees Algorithm (BA), which has demonstrated a successful implementation in diverse optimization problems [24,25,26,27]. TLBO and PSO are also used in this work for the same purpose. The aim of the work presented in this paper can be summarized as follows: (i) to propose a Fuzzy PIDF optimized by the BA and other two algorithms for load frequency control of a two-area power system and to investigate its dynamic performance; (ii) to assess the supremacy of the proposed technique by comparing the results with those of previously published works based on TLBO tuned Fuzzy PID [21] and Lozi map-based Chaotic Optimization Algorithm (LCOA) tuned PID [28]; and (iii) to investigate the robustness of the Fuzzy PIDF against a wide variation range in the parametric uncertainties of the investigated system.
Furthermore, from the comprehensive literature review, it is concluded that the proposed techniques based on different theories may provide the desired performance to overcome the problem of frequency deviation. However, the above-mentioned studies have not considered the reliability aspects in the design of the proposed schemes. This research gap has motivated the authors to suggest fuzzy control configurations for the LFC in power systems that offer different levels of reliability. Therefore, this study is then extended to propose three different fuzzy control structures, namely the Fuzzy Cascade PI − PD, the Fuzzy PI + PD, and the Fuzzy PI plus the Fuzzy PD. An extensive assessment of the robustness of these structures towards the parametric uncertainties of the testbed system, considering thirteen cases, is conducted.
The rest of this study is structured as follows. Section 2 introduces the testbed dual-area interconnected power system. Section 3 details the proposed Fuzzy PID with a filtered derivative action and the employed objective function. Section 4 gives a concise introduction about the suggested optimization technique—the Bees Algorithm. Section 5 details the main simulation results based on the proposed Fuzzy PIDF; this section also provides a comparative study between the Fuzzy PIDF tuned by the BA and the same controller tuned by PSO and TLBO, in addition to other previously published works. Then, Section 6 proposes three novel fuzzy configurations based on the suggested BA, while Section 7 investigates the robustness analyses of the proposed fuzzy structures towards the parametric uncertainties of the testbed system. Lastly, Section 8 encapsulates the main outcomes and proposes a clear path for future works based on this research.

2. Two-Area Power System—Model Understudy

The investigated system in this work is an unequal dual-area, non-reheat interconnected power system, as shown in Figure 1. This system is widely investigated in the literature on the design and analysis of the LFCs of interconnected power systems [21,28].
In this system, the generated mechanical power can be produced by gas, hydraulic, and thermal turbines. The governor senses the generator speed variation and varies the mechanical output power of the turbine by adjusting the position of the turbine input valve; this response is known as the primary frequency response. For stable operation, if the load increases, the governor will decrease the speed; this means moving the turbine input valve to a more open state (this also means increasing the torque). The coefficient of the speed regulation or the droop gain (R) represents the ratio of the frequency deviation to the generator output power variation. The relevant parameters of this model are given in Table 1.
In order to accomplish the task of the LFC, the term “Area Control Error (ACE)” is usually utilized. The ACE of each area is the input of the proposed fuzzy PIDF. Equations (1) and (2) express the ACEs of area one and area two, respectively. From the control design aspect, the desired operation of large power systems is to maintain the frequency and the tie-line power deviation, fixed on prespecified values, even in the case of load disturbance; this requires the term of ACE to be maintained at zero
ACE area   1 = P 12 +   B 1   F 1
ACE area   2 = P 21 +   B 2   F 2
where F 1 , and F 2 are the frequency deviation in areas one and two, respectively, P 12 and P 21 are the power flow deviation in areas one and two, and B 1 and B 2 are frequency biases.

3. Control Strategy and Objective Function

3.1. Fuzzy PIDF Controller

The structural design of the proposed controller equipped in area one is explained in Figure 2. As is evident from the figure, the controller has two inputs, ACE 1 and the derivative of ACE 1 , and one output. In this structure, K 1 and K 2 are the scaling factor gains of the input. While K P 1 ,   K I 1 , K D 1 , and K F 1 , which is the filter gain, are the scaling factors of the output. A similar controller is also employed in area two with K 3 and K 4 as the input scaling factors and K P 2 ,   K I 2 , K D 2 , and K F 2 for the output scaling gains. Accordingly, twelve parameters are to be optimized to obtain the desired dynamic response of the investigated system. Because of its simplicity and because less computation time is needed for this sort of membership, three triangular/two trapezoidal membership functions are utilized for the inputs and the output variables, as demonstrated in Figure 3, namely Negative Big (NB), Negative Small (NS), Zero (Z), Positive Small (PS), and Positive Big (PB). Therefore, as depicted in Table 2, twenty-five rule bases are required to generate the fuzzy output of the controller. As the controller’s performance depends on these rules, the tabulated rules are obtained by a comprehensive study of the investigated power system’s dynamic behaviour. Furthermore, the “Mamdani” interface method is utilized in this design for the fuzzification stage, while the “Centroid” tool is utilized for the defuzzification stage. Identical membership and rule base are utilised to design the Fuzzy PIDF employed in area two.

3.2. Cost Function

Prior to using the BA, TLBO, and PSO to achieve the best performance of the Fuzzy PIDF employed for the load frequency control in the two-area power system, a proper cost function should be selected. In this paper, the proposed Fuzzy PIDF is designed by minimizing the Integral Time Absolute Error (ITAE) cost function with the aid of the suggested algorithms. The used cost function is expressed in (3). The selection of this cost function is due to the fact that it was proven to reduce both the settling time and the overshoot/undershoot [7].
Objective   Function = ITAE = 0 t ( | F 1 | + | F 2 | + | P tie | ) . t .   dt  

4. 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 multiple fields have been successfully solved using the BA [24,25,26,27]. It may be worth mentioning that this is one of the first attempts to apply the BA in the field of LFC.
Similar to any optimization algorithm, a number of parameters are required to be set. Specifically, the parameters are the number of scout bees (n); number of sites selected for exploitation out of n visited sites (m); the number of best (elite) sites among the selected sites (e); the number of recruited bees in the best-selected e sites (nep); the number of bees sent to the other (m−e) selected sites (nsp); and the initial size of each patch (ngh) that includes the site and the neighbourhood area as well as the stopping criteria. The simplest pseudo-code for this algorithm is shown in Figure 4.
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 the sites visited by the n scout bees is performed. The m sites with the highest fitness in specified “chosen sites” in step 3 are selected for the neighbourhood or local search. The algorithm in steps 4 and 5 conducts searches in the neighbourhood of the selected sites, with more bees assigned to the best e sites. The selection of the best sites could be conducted directly, based on the fitness associated with them. Alternatively, using the fitness values, the probability of the sites being selected is determined. Searches in the neighbourhood 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 [29].

5. Results and Discussions

This study was carried out in MATLAB (2019a); the Bees Algorithm (BA), Teaching Learning Based Optimisation (TLBO), and Particle Swarm Optimisation (PSO) codes were programmed in (.m files); the examined dual-area power system was built in the MATLAB Simulink environment. The parameters of the BA and the PSO were set as shown in Table 3. Where C 1 and C 2 are the acceleration constants, w min and w max are the inertia weights, CR is the crossover rate, and No. Par is the number of particles. The TLBO has two parameters to be set, namely the population size and the number of iterations, which were set as 50 and 40, respectively.
A Step Load Perturbation (SLP) of 0.2 pu was applied in area one to study the dynamic performance of the system with the proposed Fuzzy PIDF. The optimum values of the Fuzzy PIDF parameters obtained by the BA, TLBO, and PSO are given in Table 4. The scaling factors of the proposed Fuzzy controller design were chosen in the limits of [0,1,2], and the filter coefficient K F was constrained in the range from 0 to 100.
Moreover, as is above mentioned, the results obtained from the proposed fuzzy structure are compared with those of the previously published studies for the same system investigated in [21] and [28]. The optimum gains of the PID tuned by the LCOA proposed in [28] and the TLBO optimized Fuzzy PID studied in [21] are given in Table 5.
The frequency deviation in both areas, F 1 and F 2 , following the implementation of 0.2 pu disturbance in area one, is shown in Figure 5 and Figure 6, respectively. The tie-line power deviation is given in Figure 7.
Figure 5, Figure 6 and Figure 7 summarize the main outcomes of the proposed Fuzzy PIDF, where it is obviously remarked that this controller offered the best response among the investigated methods. Furthermore, despite the clear similarity in the dynamic response obtained by the proposed fuzzy structure tuned by BA, TLBO, and PSO, it is observed that the BA optimized the proposed fuzzy controller and provided the best result in terms of the drop in the frequency represented by the peak undershoot occurring in area one after the 0.2 pu disturbance enforcement. However, the Fuzzy PIDF tuned by the TLBO and PSO offered the best drop in frequency in area two. The dynamic performance of the system based on the Fuzzy PIDF tuned by the suggested algorithms, the Fuzzy PID optimized by TLBO, and the PID controller tuned by LCOA, represented by undershoot, peak overshoot, and settling time in F 1 , F 2 , and P tie , is illustrated in Table 6; the value of the objective function based on each technique is also given.
Table 6 gives further verification of the superiority of the suggested controller over those presented in previous studies. The percentage of improvement in the undershoot ( U sh ), settling time ( T s ), and the ITAE for the Fuzzy configuration optimized by different algorithms and the Fuzzy PID proposed in [21], in comparison with the LCOA-based PID controller [28], is shown in Figure 8 (this figure is obtained by analyzing the characteristics provided in Table 6). From Figure 8, it is noted that with the proposed Fuzzy PIDF controller optimized by the suggested algorithms, the overall performance of the system has witnessed a remarkable improvement.
In order to examine the robustness of the Fuzzy PIDF towards the parametric uncertainties of the controlled system, several parameters of the investigated testbed system are simultaneously altered from their nominal values. The parameters Tg, Tt, and H in both areas are varied by +50%, while the parameters B and D are varied by −50%. A step load perturbation of 0.2 pu is suddenly applied (at time t = 0 s) in area one and the optimal gains of the Fuzzy PIDF obtained in the normal condition are not to be re-tuned to verify the robustness of the proposed controller. Figure 9, Figure 10 and Figure 11 and Table 7 show the dynamic performance of the two-area power system as it is exposed to a parametric deviation test with the recommended Fuzzy PIDF-based BA, TLBO, and PSO employed for LFC.
The results obtained from the robustness analysis demonstrate that the proposed fuzzy structure equipped in the testbed system for the LFC is robust towards the parametric variation of the controlled plant. It is also noticed that the same controller optimized by the TLBO has shown less robustness as compared with the same controller tuned by the BA and the PSO.
Moreover, for further robustness examination of the proposed fuzzy design at various load perturbations, a random disturbance is applied with a different magnitude in area one, as demonstrated in Figure 12. The dynamic responses of the testbed system when it is exposed to different load disturbances are shown in Figure 13, Figure 14 and Figure 15.
From Figure 13, Figure 14 and Figure 15, it is noticeable that the optimized fuzzy control structure continued to demonstrate its robustness even with random load disturbance applied in the system. Moreover, it is observed that the Fuzzy PIDF-based BA offers the best response as less peak undershoot and less oscillation are achieved in comparison with the same controller-based PSO and TLBO.

6. Different Configurations of Fuzzy Control Tuned by BA

It is revealed that the membership function selection and the setting of the rule base are vital in designing a fuzzy controller. However, it is also a significant matter to investigate the impact of the configuration of the fuzzy controller. This is to explore how different structures of the scaling factor gains influence the performance of the controller. Based on this statement, this study is extended to propose three different structures of fuzzy control employed as LFC system for the two-area power system shown in Figure 1. The same membership function used to design the Fuzzy PIDF shown in Figure 3 is used with the proposed structures. Furthermore, the rule bases required to generate the fuzzy output of the controller are tabulated in Table 2.
The three novel proposed structures are shown in Figure 16, Figure 17 and Figure 18. Due to the superiority and robustness of the performance of cascading the fraction PI and the fractional PD proposed in [30,31], it is possible to use this idea of cascading to gain further improvement on the performance of the fuzzy controller by proposing a controller that benefits from the advantages of fuzzy control and the merits of the cascading PI and PD controllers. Accordingly, the configuration illustrated in Figure 16 is for the proposed Fuzzy Cascade PI − PD (Fuzzy C PI − PD) employed in area one. This controller has six scaling factor gains. Namely, K 1 and K 2 are the input gains of the fuzzy controller, K P 11 and K I 1 are the PI controller gains, and K P 12 and K D 1 are for the PD controller gains. The identical Fuzzy C PI − PD controller is employed in area two with the following scaling factor gains: K 3 , K 4 , K P 21 , K I 2 , K P 22 , and K D 2 .
Figure 17 demonstrates the structural diagram of the proposed Fuzzy PI plus Fuzzy PD (Fuzzy PI + Fuzzy PD) controller employed in area one for LFC purposes. As is obvious from the figure, two fuzzy controllers are equipped in each area. This hierarchal configuration should enhance the stability of the system and provide better reliability when any failure occurs in any part of this structure; the other part continues to provide its expected control action.
This controller has eight scaling factor gains. Namely, K 1 and K 2 are the input gains of the Fuzzy PI controller, K P 11 and K I 1 are the PI controller gains, K 3 and K 4 are the input gains of the Fuzzy PD controller, and K P 12 and K D 1 are for the PD controller gains. The identical Fuzzy PI + Fuzzy PD is employed in area two with the following parameters: K 5 , K 6 , K P 21 , K I 2 , K 7 , K 8 ,   K P 22 , and K D 2 .
Due to the number of fuzzy rules needed to implement the Fuzzy PI + Fuzzy PD structure in each area [(5 × 5) + (5 × 5) = 50], a longer execution time is required, which may result in slowing down the controller performance; it may be worthwhile to propose another structure that reduces the execution time and guarantees a satisfactory level of reliability. In order to accomplish this, the Fuzzy (PI + PD) shown in Figure 18 is suggested. This structure has less execution time [5 × 5 = 25] and still offers an acceptable range of reliability. This configuration has six scaling factor gains. Namely, K 1 and K 2 are the input gains of the fuzzy controller, K P 11 and K I 1 are the PI controller gains, and K P 12 and K D 1 are for the PD controller gains. A similar Fuzzy C PI − PD controller is employed in area two with the following scaling factor gains: K 3 , K 4 , K P 21 , K I 2 , K P 22 , and K D 2 .
In order to achieve the best possible performance of the proposed fuzzy control configurations, the Bees Algorithm (BA) is used to concurrently find the optimal values of the proposed controllers’ gains by minimising the ITAE of the frequency deviation in both areas and the tie-line power fluctuation. The optimum values of the suggested controllers’ gains are illustrated in Table 8 and Table 9. The scaling factors of the proposed fuzzy control configurations are chosen in the limits of [0,1,2].
To investigate the performance of the proposed fuzzy control structures, a load disturbance with a magnitude of 0.2 pu is applied in area one at time t = 0. The dynamic response of the system with the proposed controllers employed as the LFCs is demonstrated in Figure 19, Figure 20 and Figure 21. The frequency deviation in area one F 1 in (Hz) is given in Figure 19; the frequency deviation in area two F 2 in (Hz) is given in Figure 20; and the tie-line power deviation P tie in (pu) is illustrated in Figure 21. Furthermore, the characteristics of the dynamic response represented by the peak undershoot ( U sh ) , the peak overshoot ( O sh ) , the settling time ( T s ) and the values of the objective function are exemplified in Table 10.
The obtained simulation results show that the steady-state responses with the proposed controllers are similar because the frequency variation in both areas and the tie-line power deviation are ceased to zero. However, in terms of the transient response, the least drop in the frequency recorded in area one, as a consequence of the implication of the load disturbance, is −0.0346 Hz; this is achieved based on the proposed Fuzzy PI plus Fuzzy PD structure. Observably, this controller offered the slowest response for F 1 , with a 7.0632 s settling time. Furthermore, a negligible overshoot is observed in the dynamic response obtained based on the three suggested controllers. Regarding the drop in the frequency in area two, the proposed structures have provided satisfactory responses with a slight superiority of the Fuzzy C PI − PD over the other two structures. The tie-line power deviation of the system is illustrated in Figure 21, where the supremacy of the Fuzzy C PI − PD over the other controllers is observed. Moreover, the value of the ITAE is the smallest for the Fuzzy C PI − PD tuned by the BA as compared with the other structures tuned by the same algorithm.
Based on the simulation results provided in Figure 19, Figure 20 and Figure 21 and Table 10, it is evidenced that the proposed fuzzy configurations are serving as effective solutions for the issue of LFC as they provide different advantages, such as fast response with neglectable overshoot and zero steady-state error. Importantly, based on their structures, these controllers offer a wide range of reliability.

7. Robustness Investigation of Fuzzy Cascade PI − PD, Fuzzy PI + PD, and Fuzzy PI plus Fuzzy PD

The investigated two-area power system has several parameters that may vary due to different operating conditions. This variation influences the stability of the system. For example, the increase in the governor time constant (Tg) results in an increase in the frequency deviation, while decreasing the damping ratio (D) may lead to increasing the frequency deviation of the system, which may bring about a possibility of system instability. Moreover, increasing the inertia time constant (H) can slow down the system. Therefore, the LFC system should have the required control action towards the parametric uncertainties of the controlled system and provide an acceptable level of robustness.
Accordingly, to assess the robustness of the proposed fuzzy control configurations tuned by the BA equipped as an LFC system in the dual-area power system, thirteen scenarios are assumed for the parametric uncertainties of the testbed system, as given in Table 11. This assessment begins with individually varying each parameter in the system by (+ and −) 50% from their nominal values. As is understood, changing the parameters from their nominal values may have a positive impact on the overall system stability. Therefore, in order to make this analysis more credible, all the parameters are simultaneously varied from their nominal values. Accordingly, in case thirteen, we consider the negative impact of each parameter of uncertainty and change them simultaneously. This guarantees the assessment of the robustness in the worst scenario that the system may experience during the working time. An SLP of 0.2 pu is applied in area one to examine the effect of the system parametric uncertainties on the behavior of the Fuzzy C PI − PD, Fuzzy PI + PD, and Fuzzy PI plus Fuzzy PD; a comparative study is given based on the obtained results. Fuzzy control robustness analysis can also be carried out in different ways as investigated in [32,33,34].
From case 1 to case 12 in Table 11, only a single parameter is changed at a time by +50% and −50% from its nominal value. In case thirteen, the parameters Tg, Tt, H, and R in both areas are varied by 50%, while the parameters B and D are varied by −50%. The dynamic response of the system with the proposed controllers under parametric uncertainty conditions is demonstrated in Figure 22, Figure 23, Figure 24, Figure 25, Figure 26, Figure 27, Figure 28, Figure 29, Figure 30, Figure 31, Figure 32, Figure 33 and Figure 34. Moreover, the characteristics of the transient response are depicted in Table 12.
Figure 22 illustrates the dynamic behavior of the testbed system based on the proposed fuzzy control structures under parametric uncertainty case 1, where the inertia time constants in both areas are altered by +50% from their nominal values. It is noted that the increase in the time inertia has slowed the response of the system. For example, the settling time of the frequency in area one has increased from 2.1873 s, 7.0632 s, and 2.1384 s to 2.380 s, 7.768 s, and 2.305 s based on Fuzzy C PI − PD, Fuzzy PI plus Fuzzy PD, and Fuzzy PI + PD, respectively. The settling time of F 2 and P tie follow the same pattern, where a slight increase is observed. Moreover, it is concluded that the increase in the inertia time constant has led to a slight decrease in the drop of the frequency in both areas. Conversely, the decrease in the inertia time constant which is considered in case 2 brings about a further drop in the frequency and the tie-line power deviation in the system. Moreover, it leads to a slight improvement in terms of the settling time in F 1 , F 2 , and P tie . Figure 23 shows the dynamic performance of the system based on the proposed controllers under parametric uncertainty case 2.
The impact of uncertainty in the turbine time constant is investigated in case 3 and case 4. Figure 24 demonstrates the dynamic performance of the testbed system based on the proposed fuzzy configurations under parametric uncertainty case 3, where the turbine time constants in both areas are varied by +50% from their nominal values. From Figure 24 and Table 12, it is noticed that as a consequence of increasing the turbine time constants within the system, the drop in the frequency in area one ( F 1 ) has jumped from −0.0431 Hz, −0.0346 Hz, and −0.0792 Hz to −0.0585 Hz, −0.0486 Hz, and −0.1074 Hz. While the drop of the frequency in area two ( F 2 ) increased from −0.00099 Hz, −0.0024 Hz, and −0.0026 Hz to −0.0013 Hz, −0.0039 Hz, and −0.0042 Hz, based on Fuzzy C PI − PD, Fuzzy PI plus Fuzzy PD, and Fuzzy PI + PD, respectively. It is also obvious that due to the increase of the turbine time constant, the settling time of the F 1 , F 2 , and P tie has slightly decreased. In contrast, from Figure 25, the decline in the turbine time constant brings about a decrease in the frequency deviation and slightly slows the system.
Figure 26 indicates the dynamic response of the system under parametric uncertainty case 5. In this case, the frequency bias in both areas is altered by +50%. It is noticed that the increase in frequency bias has marginally improved the dynamic response in terms of the drop in the frequency, where the maximum undershoot of the frequency in area one ( F 1 ) has decreased from −0.0431 Hz, −0.0346 Hz, and −0.0792 Hz to −0.0377 Hz, −0.0307 Hz, and −0.0693 Hz, while the drop of the frequency in area two ( F 2 ) declined from −0.00099 Hz, −0.0024 Hz, and −0.0026 Hz to −0.00054 Hz, −0.0013 Hz, and −0.0014 Hz, based on Fuzzy C PI − PD, Fuzzy PI plus Fuzzy PD, and Fuzzy PI + PD, respectively. With regard to the influence of decreasing the frequency bias on the stability of the power systems which are considered in case 6 and illustrated in Figure 27, it is obvious that the decrease in the frequency bias has slightly worsened the dynamic response of the system in terms of the frequency variation, where it is noted that the maximum undershoot of the frequency in area one has increased from −0.0431 Hz, −0.0346 Hz, and −0.0792 Hz to −0.0577 Hz, −0.0546 Hz, and −0.1029 Hz, while the drop of the frequency in area two ( F 2 ) increased from −0.00099 Hz, −0.0024 Hz, and −0.0026 Hz to −0.0028 Hz, −0.0080 Hz, and −0.0074 Hz, based on Fuzzy C PI − PD, Fuzzy PI plus Fuzzy PD, and Fuzzy PI + PD, respectively.
The influence of uncertainty in the damping constant (D) is investigated in cases 7 and 8. Figure 28 and Figure 29 demonstrate the dynamic performance of the testbed system based on the proposed fuzzy structures under the parametric uncertainty conditions of cases 7 and 8, where the damping constants (D) in both areas are varied by +50% and −50% from their nominal values. Due to the change in this parameter, a negligible change in the dynamic performance of the system is observed based on these cases as compared with the results obtained using the nominal conditions.
The uncertainty in the governor time constant (Tg) has a similar impact of uncertainty as the turbine time constant on the stability of the system in terms of frequency variation and the speed of the response. Figure 30 reveals the dynamic performance of the dual-area power system when the proposed fuzzy structures are employed as LFCs in the system with the consideration of parametric uncertainty case 9, where the governor time constants in both areas are varied by +50% from their nominal values. As a result of uncertainty in the governor time constants within the system, the drop in the frequency in area one ( F 1 ) has incremented from −0.0431 Hz, −0.0346 Hz, and −0.0792 Hz to −0.0573 Hz, −0.0468 Hz, and −0.1038 Hz, while the drop in the frequency in area two ( F 2 ) increased from −0.00099 Hz, −0.0024 Hz, and −0.0026 Hz to −0.0012 Hz, −0.0035 Hz, and −0.0037 Hz, based on Fuzzy C PI − PD, Fuzzy PI + PD, and Fuzzy PI plus Fuzzy PD, respectively. Moreover, the settling time of F 1 decreased from 2.1873 s, 7.0632 s, and 2.1384 s to 1.724 s, 5.788 s, and 1.7013 s, based on Fuzzy C PI − PD, Fuzzy PI plus Fuzzy PD, and Fuzzy PI + PD, respectively. The dynamic response of the system under parametric uncertainty case 10 is illustrated in Figure 31. In this case of robustness analysis, the governor time constants in both areas are varied by −50% from their nominal values. The results obtained based on case 10 revealed that the decrease in the governor time constant results in a clear decrease in the frequency variation and tie-line power deviation.
Figure 32 and Figure 33 show the dynamic response of the testbed system for parametric uncertainties in case 11 and case 12, respectively. In case 11, the regulation constant in both areas is varied by +50% while it is altered by −50% in case 12. The obtained results based on cases 12 and 13 demonstrate that the uncertainty in the regulation constant has a small impact on the system stability when the proposed fuzzy controllers are equipped in the system for load frequency control.
The worst drop in frequency in both areas ( F 1 ) and ( F 2 ), as well as in the tie-line power deviation ( P tie ), is recorded based on the results obtained from case 13 of the robustness analysis, as shown in Figure 34, where the drop of the frequency in area one has increased from −0.0431 Hz, −0.0346 Hz, and −0.0792 Hz to −0.0791 Hz, −0.0958 Hz, and −0.1354 Hz. The drop of the frequency in area two increased from −0.00099 Hz, −0.0024 Hz, and −0.0026 Hz to −0.0045 Hz, −0.0197 Hz, and −0.0179 Hz, whilst the maximum overshoot of the tie-line power deviation increased from −0.00027 pu, −0.00064 pu, and −0.00072 pu to −0.00065 pu, −0.0019 pu, and −0.0020 pu, based on Fuzzy C PI − PD, Fuzzy PI plus Fuzzy PD, and Fuzzy PI + PD, respectively.
From Figure 22, Figure 23, Figure 24, Figure 25, Figure 26, Figure 27, Figure 28, Figure 29, Figure 30, Figure 31, Figure 32, Figure 33 and Figure 34 and Table 12, despite the wide range of parametric uncertainties of the investigated two-area system in the thirteen considered scenarios, the implementation of the three fuzzy control configurations tuned by the BA suggested in this study has shown an excellent level of robustness which preserved the stability of the system within acceptable limits. Furthermore, although there is the similarity of the performance of the proposed configurations, it is obvious that the proposed Fuzzy C PI − PD and Fuzzy PI plus Fuzzy PD have outperformed the Fuzzy PI + PD in all aspects.

8. Conclusions

This study proposes a design and implementation of a Fuzzy PID with a filtered derivative action (Fuzzy PIDF) for the LFC to enhance the stability of a dual-area interconnected power system. The TLBO, PSO, and BA optimization tools were utilized to obtain the optimum values of the proposed controller parameters by reducing the ITAE objective function. A disturbance of 0.2 pu was subjected in area one to study the dynamic performance of the testbed system. The results obtained from the suggested Fuzzy PIDF controller employed for the LFC in the investigated system have been compared with those of previously published studies based on the classical PID and another design of Fuzzy PID. The obtained results revealed that the Fuzzy PIDF provides better dynamic performance as it gives the best objective function values and less undershoot for frequency and tie-line power in comparison with other controllers proposed in previous studies. For example, based on the Fuzzy PIDF tuned by the BA results, as compared with the results based on the classical PID tuned by the LCOA reported in previous study, the peak undershoot and the settling time of the frequency deviation in area one were improved by 90.345% and 40.698%, respectively, while the same characteristics of the frequency deviation in area two were improved by 94.277% and 8.403%, respectively. Furthermore, notwithstanding a wide range of variations in the power system parameters and implementing random load disturbance, it is proven that the Fuzzy PIDF is robust and has successfully kept the system stable. It is also concluded that the BA, PSO, and TLBO have demonstrated to be effective techniques for soft computing (TLBO to a lesser extent as the LFC system with Fuzzy PIDF based on the TLBO is less robust against the system parametric uncertainties in comparison with the Fuzzy PIDF tuned by the BA and PSO).
This study was further extended to propose three other fuzzy configurations for the LFC. Namely, Fuzzy Cascade PI − PD (Fuzzy C PI − PD), Fuzzy PI plus Fuzzy PD (Fuzzy PI + Fuzzy PD), and Fuzzy (PI + PD). These configurations have shown several strengths in their performance. For example, in addition to offering a robust control action with a quick response, they guarantee a higher range of reliability as compared with other structures. The Bees Algorithm was employed to find the optimum values of the scaling factor gains of the suggested configurations. An extensive examination of the impact of the parametric uncertainties of the testbed system on the performance of the proposed fuzzy control structures was conducted considering different scenarios that the system may experience in real-time operation. The obtained results based on these three structures showed that the lowest drop of the frequency in area one was −0.0431 Hz, which was achieved by the proposed Fuzzy PI + Fuzzy PD, while the lowest drop of frequency in area two was −0.00099 Hz, which was obtained by employing Fuzzy C PI − PD. The simulation results revealed that the proposed fuzzy controllers showed a high level of robustness towards the parametric uncertainties of the two-area power system (Fuzzy (PI + PD) to a lesser extent).
This research may be further extended in future work in two directions:
  • To assess the validity of the proposed fuzzy structures as LFCs in power systems that comprise Renewable Energy Resources (RESs) and to assess the impact of some nonlinear aspects within the system, such as the Governor Dead Band (GDB) and the Generation Rate Constraint (GRC) on the behavior of the suggested controllers.
  • Due to the supremacy of the Fractional Order PID controller over the traditional PID, it may be worth investigating the possible improvement on the performance of the proposed structures if the fractional-order scheme is used instead of the traditional PI − PD. Moreover, an extra improvement might be achieved if another optimization technique is used to tune the parameters of the proposed fuzzy configurations.

Author Contributions

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

Funding

This paper is part of the research of the corresponding author, M.S., who is sponsored by the Ministry of Higher Education and Scientific Research in Libya.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank Cardiff University/School of Engineering for accepting to pay the APC towards publishing this paper.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
Fuzzy PIDFFuzzy Proportional Integral Derivative with Filtered derivative mode
LFCLoad Frequency Control
BAThe Bees Algorithm
SLPStep Load Perturbation
SAMPESelf-Adaptive Multi-Population Elitist
PFMPIDPredictive Functional Modified PID
GOAGrasshopper Optimization Algorithm
PI − PDnCascade of Proportional Integral and Proportional Derivative with filter
ECOAEnhanced Coyote Optimization Algorithm
PVPhotovoltaic
MPCModel Predictive Control
SMCSliding Mode Control
LMILinear Matrix Inequalities
FLCFuzzy Logic Control
PSOParticle Swarm Optimization
DEDifferential Evolution
TLBOTeaching Learning Based Optimization
PSPattern Search
LCOALozi map- based Chaotic Optimization Algorithm
ACEArea Control Error
NBNegative Big
NSNegative Small
ZZero
PSPositive Small
PBPositive Big
ITAEIntegral Time Absolute Error
Fuzzy C PI − PDFuzzy Cascade PI − PD
Fuzzy PI + Fuzzy PDFuzzy PI plus Fuzzy PD
RESsRenewable Energy Resources
GDBGovernor Dead Band
GRCGeneration Rate Constraint

References

  1. Shouran, M.; Anayi, F.; Packianather, M. Fuzzy PID with Filtered Derivative Mode Based Load Frequency Control of Two-Area Power System. In Proceedings of the 2021 56th International Universities Power Engineering Conference (UPEC), Virtuell, UK, 31 August–3 September 2021; pp. 1–6. [Google Scholar] [CrossRef]
  2. Barisal, A. Comparative performance analysis of teaching learning based optimization for automatic load frequency control of multi-source power systems. Int. J. Electr. Power Energy Syst. 2015, 66, 67–77. [Google Scholar] [CrossRef]
  3. Sekhar, G.C.; Sahu, R.K.; Baliarsingh, A.; Panda, S. Load frequency control of power system under deregulated environment using optimal firefly algorithm. Int. J. Electr. Power Energy Syst. 2016, 74, 195–211. [Google Scholar] [CrossRef]
  4. Kumar, N.; Tyagi, B.; Kumar, V. Application of Fractional Order PID Controller for AGC Under Deregulated Environment. Int. J. Autom. Comput. 2018, 15, 84–93. [Google Scholar] [CrossRef]
  5. Gheisarnejad, M.; Khooban, M.H. Design an optimal fuzzy fractional proportional integral derivative controller with derivative filter for load frequency control in power systems. Trans. Inst. Meas. Control 2019, 41, 2563–2581. [Google Scholar] [CrossRef]
  6. Alhelou, H.H.; Hamedani-Golshan, M.E.; Zamani, R.; Heydarian-Forushani, E.; Siano, P. Challenges and opportunities of load frequency control in conventional, modern and future smart power systems: A comprehensive review. Energies 2018, 11, 2497. [Google Scholar] [CrossRef] [Green Version]
  7. Kumar, B.; Kumar, T.; Ranjan, J.; Panda, S.; Kumar, S. A novel hybrid LUS–TLBO optimized fuzzy-PID controller for load frequency control of multi-source power system. Int. J. Electr. Power Energy Syst. 2016, 74, 58–69. [Google Scholar]
  8. El Ela, A.A.A.; El-Sehiemy, R.A.; Shaheen, A.M.; Diab, A.E.G. Optimal Design of PID Controller Based Sampe-Jaya Algorithm for Load Frequency Control of Linear and Nonlinear Multi-Area Thermal Power Systems. Int. J. Eng. Res. Afr. 2020, 50, 79–93. [Google Scholar] [CrossRef]
  9. Nosratabadi, S.M.; Bornapour, M.; Gharaei, M.A. Grasshopper optimization algorithm for optimal load frequency control considering Predictive Functional Modified PID controller in restructured multi-resource multi-area power system with Redox Flow Battery units. Control Eng. Pr. 2019, 89, 204–227. [Google Scholar] [CrossRef]
  10. El-Ela, A.A.A.; El-Sehiemy, R.A.; Shaheen, A.M.; Diab, A.E.-G. Enhanced coyote optimizer-based cascaded load frequency controllers in multi-area power systems with renewable. Neural Comput. Appl. 2021, 33, 8459–8477. [Google Scholar] [CrossRef]
  11. Arya, Y. Automatic generation control of two-area electrical power systems via optimal fuzzy classical controller. J. Frankl. Inst. 2018, 355, 2662–2688. [Google Scholar] [CrossRef]
  12. Zeng, G.Q.; Xie, X.Q.; Chen, M.R. An adaptive model predictive load frequency control method for multi-area inter-connected power systems with photovoltaic generations. Energies 2017, 10, 1840. [Google Scholar] [CrossRef] [Green Version]
  13. Shouran, M.; Anayi, F.; Packianather, M. Design of sliding mode control optimised by the Bees algorithm for LFC in the Great Britain power system. Mater. Today Proc. 2021, in press. [CrossRef]
  14. Guo, J. Application of a novel adaptive sliding mode control method to the load frequency control. Eur. J. Control 2020, 57, 172–178. [Google Scholar] [CrossRef]
  15. Shangguan, X.; He, Y.; Zhang, C.; Jiang, L.; Spencer, J.W.; Wu, M. Sampled-data based discrete and fast load frequency control for power systems with wind power. Appl. Energy 2019, 259, 114202. [Google Scholar] [CrossRef]
  16. Yu, 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] [CrossRef]
  17. Singh, V.P.; Mohanty, S.R.; Kishor, N.; Ray, P.K. Robust H-infinity load frequency control in hybrid distributed generation system. Int. J. Electr. Power Energy Syst. 2013, 46, 294–305. [Google Scholar] [CrossRef]
  18. Obaid, Z.A.; Cipcigan, L.; Muhssin, M.T. Fuzzy hierarchal approach-based optimal frequency control in the Great Britain power system. Electr. Power Syst. Res. 2016, 141, 529–537. [Google Scholar] [CrossRef]
  19. Sahu, R.K.; Sekhar, G.T.C.; Panda, S. DE optimized fuzzy PID controller with derivative filter for LFC of multi source power system in deregulated environme nt. Ain Shams Eng. J. 2015, 6, 511–530. [Google Scholar] [CrossRef] [Green Version]
  20. Jalali, N.; Razmi, H.; Doagou-Mojarrad, H. Optimized fuzzy self-tuning PID controller design based on Tribe-DE optimization algorithm and rule weight adjustment method for load frequency control of interconnected multi-area power systems. Appl. Soft Comput. 2020, 93, 106424. [Google Scholar] [CrossRef]
  21. Sahu, B.K.; Pati, S.; Mohanty, P.K.; Panda, S. Teaching-learning based optimization algorithm based fuzzy-PID controller for automatic generation control of multi-area power system. Appl. Soft Comput. J. 2015, 27, 240–249. [Google Scholar] [CrossRef]
  22. Sahu, R.K.; Panda, S.; Yegireddy, N.K. A novel hybrid DEPS optimized fuzzy PI/PID controller for load frequency control of multi-area interconnected power systems. J. Process Control 2014, 24, 1596–1608. [Google Scholar] [CrossRef]
  23. Shouran, M.; Anayi, F.; Packianather, M. A State-of-the-Art Review on LFC Strategies in Conventional and Modern Power Systems. In Proceedings of the 2021 International Conference on Advance Computing and Innovative Technologies in Engineering (ICACITE), Greater Noida, India, 4–5 March 2021; pp. 268–277. [Google Scholar]
  24. Pham, D.; Ghanbarzadeh, A.; Koç, E.; Otri, S.; Rahim, S.; Zaidi, M. The Bees Algorithm—A Novel Tool for Complex Optimisation Problems. In Intelligent Production Machines and Systems; Elsevier: Amsterdam, The Netherlands, 2006; pp. 454–459. [Google Scholar] [CrossRef]
  25. Pham, D.T.; Castellani, M. The bees algorithm: Modelling foraging behaviour to solve continuous optimization problems. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2009, 223, 2919–2938. [Google Scholar] [CrossRef]
  26. Castellani, M.; Otri, S.; Pham, D.T. Printed circuit board assembly time minimisation using a novel Bees Algorithm. Comput. Ind. Eng. 2019, 133, 186–194. [Google Scholar] [CrossRef]
  27. Pham, D.T.; Koç, E. Design of a two-dimensional recursive filter using the bees algorithm. Int. J. Autom. Comput. 2010, 7, 399–402. [Google Scholar] [CrossRef]
  28. Farahani, M.; Ganjefar, S.; Alizadeh, M. PID controller adjustment using chaotic optimisation algorithm for multi-area load frequency control. IET Control. Theory Appl. 2012, 6, 1984–1992. [Google Scholar] [CrossRef]
  29. Shouran, M.; Anayi, F.; Packianather, M. The Bees Algorithm Tuned Sliding Mode Control for Load Frequency Control in Two-Area Power System. Energies 2021, 14, 5701. [Google Scholar] [CrossRef]
  30. Shouran, M.; Alsseid, A.M. Cascade of Fractional Order PID based PSO Algorithm for LFC in Two-Area Power System. In Proceedings of the 2021 3rd International Conference on Electronics Representation and Algorithm (ICERA), Yogyakarta, Indonesia, 29–30 January 2021; pp. 1–6. [Google Scholar] [CrossRef]
  31. Çelik, E. Design of new fractional order PI–fractional order PD cascade controller through dragonfly search algorithm for advanced load frequency control of power systems. Soft Comput. 2020, 25, 1193–1217. [Google Scholar] [CrossRef]
  32. Hosseinzadeh, M.; Sadati, N.; Zamani, I. H∞ disturbance attenuation of fuzzy large-scale systems. In Proceedings of the 2011 IEEE International Conference on Fuzzy Systems (FUZZ-IEEE 2011), Taipei, Taiwan, 27–30 June 2011; pp. 2364–2368. [Google Scholar] [CrossRef]
  33. Zhu, J.; Nguang, S.K. Fuzzy Model Predictive Control with Enhanced Robustness for Nonlinear System via a Discrete Disturbance Observer. IEEE Access 2020, 8, 220631–220645. [Google Scholar] [CrossRef]
  34. Petrovic, D.; Kalata, M.; Luo, J. A fuzzy scenario-based optimisation of supply network cost, robustness and shortages. Comput. Ind. Eng. 2021, 160, 107555. [Google Scholar] [CrossRef]
Figure 1. Transfer function model of the testbed system.
Figure 1. Transfer function model of the testbed system.
Energies 15 00657 g001
Figure 2. Structural diagram of Fuzzy PIDF controller of area 1.
Figure 2. Structural diagram of Fuzzy PIDF controller of area 1.
Energies 15 00657 g002
Figure 3. Membership functions of the two inputs and output.
Figure 3. Membership functions of the two inputs and output.
Energies 15 00657 g003
Figure 4. The Bees Algorithm flowchart.
Figure 4. The Bees Algorithm flowchart.
Energies 15 00657 g004
Figure 5. Frequency variation in area one ( F 1 in Hz).
Figure 5. Frequency variation in area one ( F 1 in Hz).
Energies 15 00657 g005
Figure 6. Frequency variation in area two ( F 2 in Hz).
Figure 6. Frequency variation in area two ( F 2 in Hz).
Energies 15 00657 g006
Figure 7. Tie-line power variation ( P tie in pu).
Figure 7. Tie-line power variation ( P tie in pu).
Energies 15 00657 g007
Figure 8. Percentage of improvement with different techniques.
Figure 8. Percentage of improvement with different techniques.
Energies 15 00657 g008
Figure 9. Frequency deviation in area one ( F 1 in Hz) under parametric uncertainties of the testbed system.
Figure 9. Frequency deviation in area one ( F 1 in Hz) under parametric uncertainties of the testbed system.
Energies 15 00657 g009
Figure 10. Frequency deviation in area two ( F 2 in Hz) under parametric uncertainties of the testbed system.
Figure 10. Frequency deviation in area two ( F 2 in Hz) under parametric uncertainties of the testbed system.
Energies 15 00657 g010
Figure 11. Tie-line power deviation ( P tie in pu) under parametric uncertainties of the testbed system.
Figure 11. Tie-line power deviation ( P tie in pu) under parametric uncertainties of the testbed system.
Energies 15 00657 g011
Figure 12. Random load profile.
Figure 12. Random load profile.
Energies 15 00657 g012
Figure 13. Frequency deviation in area one: (A) based on BA tuning and (B) based on TLBO and PSO tuning.
Figure 13. Frequency deviation in area one: (A) based on BA tuning and (B) based on TLBO and PSO tuning.
Energies 15 00657 g013
Figure 14. Frequency deviation in area two: (A) based on BA tuning and (B) based on TLBO and PSO tuning.
Figure 14. Frequency deviation in area two: (A) based on BA tuning and (B) based on TLBO and PSO tuning.
Energies 15 00657 g014
Figure 15. Tie-line power deviation ( P tie in pu): (A) based on BA tuning and (B) based on TLBO and PSO tuning.
Figure 15. Tie-line power deviation ( P tie in pu): (A) based on BA tuning and (B) based on TLBO and PSO tuning.
Energies 15 00657 g015
Figure 16. Block diagram of Fuzzy Cascade PI − PD controller configuration equipped in area one.
Figure 16. Block diagram of Fuzzy Cascade PI − PD controller configuration equipped in area one.
Energies 15 00657 g016
Figure 17. Block diagram of Fuzzy PI plus Fuzzy PD controller configuration equipped in area one.
Figure 17. Block diagram of Fuzzy PI plus Fuzzy PD controller configuration equipped in area one.
Energies 15 00657 g017
Figure 18. Block diagram of Fuzzy (PI + PD) controller configuration equipped in area one.
Figure 18. Block diagram of Fuzzy (PI + PD) controller configuration equipped in area one.
Energies 15 00657 g018
Figure 19. Frequency deviation in area one ( F 1 in Hz).
Figure 19. Frequency deviation in area one ( F 1 in Hz).
Energies 15 00657 g019
Figure 20. Frequency deviation in area two ( F 2 in Hz).
Figure 20. Frequency deviation in area two ( F 2 in Hz).
Energies 15 00657 g020
Figure 21. Tie-line power deviation ( P tie in pu).
Figure 21. Tie-line power deviation ( P tie in pu).
Energies 15 00657 g021
Figure 22. Dynamic response of the testbed power system based on different fuzzy controllers under parametric uncertainty conditions, case 1. (A) Frequency variation in area 1; (B) frequency variation in area 2; (C) tie-line power variation.
Figure 22. Dynamic response of the testbed power system based on different fuzzy controllers under parametric uncertainty conditions, case 1. (A) Frequency variation in area 1; (B) frequency variation in area 2; (C) tie-line power variation.
Energies 15 00657 g022
Figure 23. Dynamic response of the testbed power system based on different fuzzy controllers under parametric uncertainty conditions, case 2. (A) Frequency variation in area 1; (B) frequency variation in area 2; (C) tie-line power variation.
Figure 23. Dynamic response of the testbed power system based on different fuzzy controllers under parametric uncertainty conditions, case 2. (A) Frequency variation in area 1; (B) frequency variation in area 2; (C) tie-line power variation.
Energies 15 00657 g023
Figure 24. Dynamic response of the testbed power system based on different fuzzy controllers under parametric uncertainty conditions, case 3. (A) Frequency variation in area 1; (B) frequency variation in area 2; (C) tie-line power variation.
Figure 24. Dynamic response of the testbed power system based on different fuzzy controllers under parametric uncertainty conditions, case 3. (A) Frequency variation in area 1; (B) frequency variation in area 2; (C) tie-line power variation.
Energies 15 00657 g024
Figure 25. Dynamic response of the testbed power system based on different fuzzy controllers under parametric uncertainty conditions, case 4. (A) Frequency variation in area 1; (B) frequency variation in area 2; (C) tie-line power variation.
Figure 25. Dynamic response of the testbed power system based on different fuzzy controllers under parametric uncertainty conditions, case 4. (A) Frequency variation in area 1; (B) frequency variation in area 2; (C) tie-line power variation.
Energies 15 00657 g025
Figure 26. Dynamic response of the testbed power system based on different fuzzy controllers under parametric uncertainty conditions, case 5. (A) Frequency variation in area 1; (B) frequency variation in area 2; (C) tie-line power variation.
Figure 26. Dynamic response of the testbed power system based on different fuzzy controllers under parametric uncertainty conditions, case 5. (A) Frequency variation in area 1; (B) frequency variation in area 2; (C) tie-line power variation.
Energies 15 00657 g026
Figure 27. Dynamic response of the testbed power system based on different fuzzy controllers under parametric uncertainty conditions, case 6. (A) Frequency variation in area 1; (B) frequency variation in area 2; (C) tie-line power variation.
Figure 27. Dynamic response of the testbed power system based on different fuzzy controllers under parametric uncertainty conditions, case 6. (A) Frequency variation in area 1; (B) frequency variation in area 2; (C) tie-line power variation.
Energies 15 00657 g027
Figure 28. Dynamic response of the testbed power system based on different fuzzy controllers under parametric uncertainty conditions, case 7. (A) Frequency variation in area 1; (B) frequency variation in area 2; (C) tie-line power variation.
Figure 28. Dynamic response of the testbed power system based on different fuzzy controllers under parametric uncertainty conditions, case 7. (A) Frequency variation in area 1; (B) frequency variation in area 2; (C) tie-line power variation.
Energies 15 00657 g028
Figure 29. Dynamic response of the testbed power system based on different fuzzy controllers under parametric uncertainty conditions, case 8. (A) Frequency variation in area 1; (B) frequency variation in area 2; (C) tie-line power variation.
Figure 29. Dynamic response of the testbed power system based on different fuzzy controllers under parametric uncertainty conditions, case 8. (A) Frequency variation in area 1; (B) frequency variation in area 2; (C) tie-line power variation.
Energies 15 00657 g029
Figure 30. Dynamic response of the testbed power system based on different fuzzy controllers under parametric uncertainty conditions, case 9. (A) Frequency variation in area 1; (B) frequency variation in area 2; (C) tie-line power variation.
Figure 30. Dynamic response of the testbed power system based on different fuzzy controllers under parametric uncertainty conditions, case 9. (A) Frequency variation in area 1; (B) frequency variation in area 2; (C) tie-line power variation.
Energies 15 00657 g030
Figure 31. Dynamic response of the testbed power system based on different fuzzy controllers under parametric uncertainty conditions, case 10. (A) Frequency variation in area 1; (B) frequency variation in area 2; (C) tie-line power variation.
Figure 31. Dynamic response of the testbed power system based on different fuzzy controllers under parametric uncertainty conditions, case 10. (A) Frequency variation in area 1; (B) frequency variation in area 2; (C) tie-line power variation.
Energies 15 00657 g031
Figure 32. Dynamic response of the testbed power system based on different fuzzy controllers under parametric uncertainty conditions, case 11. (A) Frequency variation in area 1; (B) frequency variation in area 2; (C) tie-line power variation.
Figure 32. Dynamic response of the testbed power system based on different fuzzy controllers under parametric uncertainty conditions, case 11. (A) Frequency variation in area 1; (B) frequency variation in area 2; (C) tie-line power variation.
Energies 15 00657 g032
Figure 33. Dynamic response of the testbed power system based on different fuzzy controllers under parametric uncertainty conditions, case 12. (A) Frequency variation in area 1; (B) frequency variation in area 2; (C) tie-line power variation.
Figure 33. Dynamic response of the testbed power system based on different fuzzy controllers under parametric uncertainty conditions, case 12. (A) Frequency variation in area 1; (B) frequency variation in area 2; (C) tie-line power variation.
Energies 15 00657 g033
Figure 34. Dynamic response of the testbed power system based on different fuzzy controllers under parametric uncertainty conditions, case 13. (A) Frequency variation in area 1; (B) frequency variation in area 2; (C) tie-line power variation.
Figure 34. Dynamic response of the testbed power system based on different fuzzy controllers under parametric uncertainty conditions, case 13. (A) Frequency variation in area 1; (B) frequency variation in area 2; (C) tie-line power variation.
Energies 15 00657 g034
Table 1. The parameters of the testbed system.
Table 1. The parameters of the testbed system.
ParametersDefinitionValues in Area 1Values in Area 2
RRegulation constant0.05 MW/Hz0.0625 MW/Hz
BFrequency bias20.6 Hz/MW16.9 Hz/MW
DThe ratio of change in load to change in frequency0.60.9
HSystem inertia time constant54
T g Governor time constant0.2 s0.3 s
T t Turbine time constant0.5 s0.6 s
TSynchronization coefficient2
FFrequency of the system60 Hz
SLPStep Load Perturbation0.2 pu
Table 2. Fuzzy rule base of the proposed controller.
Table 2. Fuzzy rule base of the proposed controller.
ACE   1 ACE   1 ˙
NBNSZPSPB
NBNBNBNBNSZ
NSNBNBNSZPS
ZNBNSZPSPB
PSNSZPSPBPB
PBZPSPBPBPB
Table 3. The BA and PSO parameters.
Table 3. The BA and PSO parameters.
ControllerParameters
BANmenepnspngh
301261170.011
PSONo. Par w min w max C 1 C 2 CR
300.40.9220.65
Table 4. Gains of Fuzzy PIDF controllers tuned by BA, TLBO, and PSO.
Table 4. Gains of Fuzzy PIDF controllers tuned by BA, TLBO, and PSO.
ControllerController Gains of Area 1Controller Gains of Area 2
K 1 K 2 K P 1 K I 1 K D 1 K F 1 K 3 K 4 K P 2 K I 2 K D 2 K F 2
Fuzzy PIDF-BA0.403222298.48410.26481.00810.91331.97301.988993.8922
Fuzzy PIDF-TLBO 0.0351.99921.99861.998681.999599.06061.96020.037070.44351.30030.01999.7446
Fuzzy PIDF-PSO0.0222221002220.0151.403511.21
Table 5. Gains of PID controllers tuned by LCOA and Fuzzy PID controllers tuned by TLBO.
Table 5. Gains of PID controllers tuned by LCOA and Fuzzy PID controllers tuned by TLBO.
ControllerController Gains of Area 1Controller Gains of Area 2
PID-LCOA [28] K P 1 K I 1 K D 1 K P 2 K I 2 K D 2
0.939 0.7998 0.5636 0.5208 0.4775 0.7088
Fuzzy PID-TLBO [21] K 1 K 2 K 3 K 4 K 5 K 6 K 7 K 8
1.98571.99681.68701.98761.34691.55120.80980.5043
Table 6. Characteristics of the testbed system with several controllers.
Table 6. Characteristics of the testbed system with several controllers.
ControllerFrequency in Area 1Frequency in Area 2Tie-Line Power DeviationITAE
U sh   in   Hz O sh   in   Hz T s   in   s U sh   in   Hz O sh   in   Hz T s   in   s U sh   in   pu O sh   in   pu T s   in   s
BA-Fuzzy PIDF0.04140.00416.94010.0038019.29910.0010019.3600.0361
TLBO-Fuzzy PIDF0.08680.00405.75440.0036019.32730.00099018.8930.0304
PSO-Fuzzy PIDF0.08900.00405.71750.0036019.10200.0010019.1540.0330
TLBO-Fuzzy [21]0.18850.00364.99360.019023.51880.0042023.9370.3264
LCOA-PID [28]0.42880.015511.7030.0664021.06980.0134021.9780.7842
Table 7. Frequency response performances with different controllers for parametric uncertainties analysis.
Table 7. Frequency response performances with different controllers for parametric uncertainties analysis.
ControllerFrequency in Area 1Frequency in Area 2Tie-Line Power DeviationITAE
U sh   in   Hz O sh   in   Hz T s   in   s U sh   in   Hz O sh   in   Hz T s   in   s U sh   in   pu O sh   in   pu T s   in   s
BA-Fuzzy PIDF0.11400.01315.98580.020309.37810.0026010.4530.03094
TLBO-Fuzzy PIDF0.14580.01115.43780.02780.0018314.8180.00260.0000659.37690.0511
PSO-Fuzzy PIDF0.14650.01155.44680.0175010.2690.0024010.4210.02535
Table 8. The optimum values of the proposed Fuzzy C PI − PD and Fuzzy (PI + PD) controllers obtained by the BA.
Table 8. The optimum values of the proposed Fuzzy C PI − PD and Fuzzy (PI + PD) controllers obtained by the BA.
ControllerController Gains of Area 1Controller Gains of Area 2
K 1 K 2 K P 11 K I 1 K P 12 K D 1 K 3 K 4 K P 21 K I 2 K P 22 K D 2
Fuzzy C PI − PD0.0833222221.53070.00121.64720.0120.99021.1494
Fuzzy (PI + PD) 0.05942221.5673220.00020.0011.07160.41880.092
Table 9. The optimum Fuzzy PI + Fuzzy PD gains optimized by BA.
Table 9. The optimum Fuzzy PI + Fuzzy PD gains optimized by BA.
ControllerParameters
Fuzzy PI + Fuzzy PDController gains of area 1
K 1 K 2 K 3 K 4 K P 11 K I 1 K P 12 K D 1
0.002020.598122222
Controller gains of area 2
K 5 K 6 K 7 K 8 K P 21 K I 2 K P 22 K D 2
20.0010.00040.00150.34290.75111.19972
Table 10. Frequency response performances with different fuzzy control structures tuned by BA.
Table 10. Frequency response performances with different fuzzy control structures tuned by BA.
ControllerFrequency in Area 1Frequency in Area 2Tie-Line Power DeviationITAE
U sh   in   Hz O sh   in   Hz T s   in   s U sh   in   Hz O sh   in   Hz T s   in   s U sh   in   pu O sh in   pu T s   in   s
Fuzzy C PI − PD0.04310.000382.18730.00099021.17030.00027021.75220.01351
Fuzzy PI + Fuzzy PD0.03460.000897.06320.0024020.50110.00064020.81530.02576
Fuzzy (PI + PD)0.07920.001202.13840.0026021.06810.00072021.08850.03104
Table 11. Investigated scenarios of system parametric variations.
Table 11. Investigated scenarios of system parametric variations.
CaseParametersNominal ValuesVariation RangeNew Values
NumberArea 1Area 2Area 1Area 2
Case 1H54+50%7.56
Case 2H54−50%2.52
Case 3Tt0.50.6+50%0.750.9
Case 4Tt0.50.6−50%0.250.3
Case 5B20.616.9+50%30.925.35
Case 6B20.616.9−50%10.38.45
Case 7D0.60.9+50%0.91.35
Case 8D0.60.9−50%0.30.45
Case 9Tg0.20.3+50%0.30.45
Case 10Tg0.20.3−50%0.10.15
Case 11R0.050.0625+50%0.0750.0937
Case 12R0.050.0625−50%0.0250.0312
Case 13B20.616.9−50%10.38.45
H54+50%7.56
R0.050.0625+50%0.0750.0937
D0.60.9−50%0.30.45
Tt0.50.6+50%0.750.9
Tg0.20.3+50%0.30.45
Table 12. Dynamic response of the system under different parametric uncertainties scenarios with different controllers.
Table 12. Dynamic response of the system under different parametric uncertainties scenarios with different controllers.
Case NoControllerFrequency in Area 1Frequency in Area 2Tie-Line Power Deviation
U sh   in   Hz O sh   in   Hz T s   in   s U sh   in   Hz O sh   in   Hz T s   in   s U sh   in   pu O sh   in   pu T s   in   s
Case 1Fuzzy C PI − PD−0.03390.000372.380−0.00095021.4962−0.00026021.697
Fuzzy PI + Fuzzy PD−0.02970.000877.768−0.0023020.7037−0.00061020.895
Fuzzy (PI + PD)−0.06150.00112.305−0.0024021.107−0.00067021.117
Case 2Fuzzy C PI − PD−0.06950.000441.713−0.0012020.714−0.00032021.688
Fuzzy PI + Fuzzy PD−0.05790.00122.455−0.0032019.835−0.00080020.300
Fuzzy (PI + PD)−0.12680.00101.765−0.0030020.947−0.00085020.972
Case 3Fuzzy C PI − PD−0.05850.00051.5501−0.0013020.157−0.00036021.1190
Fuzzy PI + Fuzzy PD−0.04860.00236.0375−0.0039018.896−0.00090019.8623
Fuzzy (PI + PD)−0.10740.00171.5144−0.0042019.826−0.0010020.4457
Case 4Fuzzy C PI − PD−0.02600.000322.955−0.00082022.077−0.00023022.223
Fuzzy PI + Fuzzy PD−0.02200.000759.0899−0.0018021.720−0.00049021.757
Fuzzy (PI + PD)−0.04830.000922.851−0.0019021.611−0.00054021.633
Case 5Fuzzy C PI − PD−0.03770.000252.0857−0.00054031.659−0.00022025.188
Fuzzy PI + Fuzzy PD−0.03070.000561.6032−0.0013027.7717−0.00051028.1434
Fuzzy (PI + PD)−0.06930.000772.1015−0.0014028.6358−0.00058028.660
Case 6Fuzzy C PI − PD−0.05770.000742.1637−0.0028019.4344−0.00040026.209
Fuzzy PI + Fuzzy PD−0.05460.00186.3768−0.0080011.6084−0.00099012.1848
Fuzzy (PI + PD)−0.10290.00214.8285−0.0074011.8962−0.0010012.1313
Case 7Fuzzy C PI − PD−0.04300.000382.1892−0.00099020.927−0.00027022.192
Fuzzy PI + Fuzzy PD−0.03460.000987.0758−0.0024020.531−0.00064020.837
Fuzzy (PI + PD)−0.07900.00122.141−0.0026021.084−0.00072021.104
Case 8Fuzzy C PI − PD−0.04310.000382.1854−0.0010021.418−0.00027021.325
Fuzzy PI + Fuzzy PD−0.03470.000897.0507−0.0024020.470−0.00064020.792
Fuzzy (PI + PD)−0.07940.00122.1359−0.0026020.052−0.00072021.072
Case 9Fuzzy C PI − PD−0.05730.000471.724−0.0012020.633−0.00034021.414
Fuzzy PI + Fuzzy PD−0.04680.00135.788−0.0035019.235−0.00079020.257
Fuzzy (PI + PD)−0.10380.00321.7013−0.0037020.298−0.00094020.738
Case 10Fuzzy C PI − PD−0.02670.000342.7138−0.00086021.6312−0.00024021.995
Fuzzy PI + Fuzzy PD −0.02240.00808.9968−0.0019021.285−0.00053021.362
Fuzzy (PI + PD)−0.05070.000982.5984−0.0021021.339−0.00057021.378
Case 11Fuzzy C PI − PD−0.04310.000382.1736−0.001024.738−0.00027017.791
Fuzzy PI + Fuzzy PD−0.03460.00896.953−0.0025020.2308−0.00064020.496
Fuzzy (PI + PD)−0.07920.00122.1284−0.0026020.848−0.00072020.863
Case 12Fuzzy C PI − PD−0.04300.000412.212−0.00095014.328−0.00027032.2229
Fuzzy PI + Fuzzy PD−0.03460.000897.3773−0.0023021.268−0.00063021.7101
Fuzzy (PI + PD)−0.07900.00112.1689−0.0024021.693−0.00071021.7307
Case 13Fuzzy C PI − PD−0.07910.010271.435−0.00450.000049.9188−0.00065014.922
Fuzzy PI + Fuzzy PD−0.09580.01035.1487−0.0197010.577−0.001900.00000710.393
Fuzzy (PI + PD)−0.13540.0204.9656−0.0179010.1364−0.00200.00004710.681
Values that represent the best characteristics are indicated in bold.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shouran, M.; Anayi, F.; Packianather, M.; Habil, M. Different Fuzzy Control Configurations Tuned by the Bees Algorithm for LFC of Two-Area Power System. Energies 2022, 15, 657. https://doi.org/10.3390/en15020657

AMA Style

Shouran M, Anayi F, Packianather M, Habil M. Different Fuzzy Control Configurations Tuned by the Bees Algorithm for LFC of Two-Area Power System. Energies. 2022; 15(2):657. https://doi.org/10.3390/en15020657

Chicago/Turabian Style

Shouran, Mokhtar, Fatih Anayi, Michael Packianather, and Monier Habil. 2022. "Different Fuzzy Control Configurations Tuned by the Bees Algorithm for LFC of Two-Area Power System" Energies 15, no. 2: 657. https://doi.org/10.3390/en15020657

APA Style

Shouran, M., Anayi, F., Packianather, M., & Habil, M. (2022). Different Fuzzy Control Configurations Tuned by the Bees Algorithm for LFC of Two-Area Power System. Energies, 15(2), 657. https://doi.org/10.3390/en15020657

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