Water Cycle Algorithm Optimized Type II Fuzzy Controller for Load Frequency Control of a Multi-Area, Multi-Fuel System with Communication Time Delays

: This paper puts forward the implementation of an intelligent type II fuzzy PID (T2-FPID) controller tweaked with a water cycle algorithm (WCA), subjected to an error multiplied with time area over integral (ITAE) objective index for regularizing the variations in frequency and interline power ﬂow of an interconnected power system during load disturbances. The WCA-based T2-FPID is tested on a multi-area (MA) system comprising thermal-hydro-nuclear (THN) (MATHN) plants in each area. The dynamical behavior of the system is analyzed upon penetrating area 1 with a step load perturbation (SLP) of 10%. However, power system practicality constraints, such as generation rate constraints (GRCs) and time delays in communication (CTDs), are examined. Afterward, a territorial control scheme of a superconducting magnetic energy storage system (SMES) and a uniﬁed power ﬂow controller (UPFC) is installed to further enhance the system performance. The dominancy of the presented WCA-tuned T2-FPID is revealed by testing it on a widely used dual-area hydro-thermal (DAHT) power system model named test system 1 in this paper. Analysis reveals the efﬁcacy of the presented controller with other approaches reported in the recent literature. Finally, secondary and territorial regulation schemes are subjected to sensitivity analysis to deliberate the robustness.


Introduction
In the present scenario, rapid industrialization demands more electrical energy consumption. Participation of electrical utilities in the existing power system has been rapidly increasing to meet this ever-changing power demand. With the penetration of these electrical utilities, an interconnected system becomes more complex. The electric utilities are sub-grouped as control areas, which connect the remaining areas via transmission lines named tie lines. As soon as a small load disturbance occurs on any of the control areas, there exists a considerable disturbance in the area frequency, making the entire interconnected system unhealthy. However, the stability and efficacy of the power system are more likely to be relied on damping out the variations in line flow and frequency that arise under disturbance conditions at the earliest. This task can be effectively addressed by a load frequency controller (LFC). An LFC can accomplish the mechanism of balancing the power through a secondary regulator by altering the set point valve of the generation unit, thereby regulating the system frequency. Thus, a supreme and robust secondary regulatory control technique is essential to neutralize the impact of load variations on system performance.
Tungadio et al. [1] reported various power system models that have been deliberated by researchers in LFC study. Elgerd et al. [2] initiated the visualization of isolated thermal plant behavior under load disturbances. Since then, more researchers are enlightened to contribute governing techniques in the LFC domain. Nanda et al. [3] investigated two-area and three-area thermal systems with unequal generation capacity. Researchers in [4] carried out the analysis of the performance assessment of a DAHT system on load disturbances. Ravi et al. [5] investigated the traditional dual-area system incorporated with traditional generating sources of thermal-hydro-gas units in each area. Authors in [6] analyzed DAHT system behavior without and with a generation limiter of a GRC in order to demonstrate its impact on system performance. Later, researchers concentrated on considering the integration of renewable plants with existing conventional units for the study of LFCs. Irrespective of the test system models, researchers are more concentrated on designing and optimizing secondary regulators (SRs) to improve power system performance [7].
Though the aforementioned optimizations are widely used by researchers, most of them get trapped into local minima and are unable to use proper search space [19]. This opens the space and scope for the use of novel meta-heuristic searching methods in the study of LFCs. So, in this paper, an effective approach of WCA is used to optimize T2-FPID regulator parameters to maintain the stability of interconnected systems. Though this WCA is used in LFC study, the combination of the WCA approach and the T2-FPID controller has not been studied so far and no literature is reported. This inspires the authors to implement this control mechanism along with territorial strategies of FACTS and energy storage devices (ESDs) for further enhancement. A complex interconnected system definitely needs an additional regulatory mechanism as SRs are unable to regulate the system deviations under sudden large disturbances. An additional strategy of a UPFC and an SMES is implemented in this work because of their significant features.
In view of the above, the contributions of this paper are: (a) The performance of the WCA-optimized T2-FPID controller is assessed for LFC study.
(b) The efficacy of the proposed WCA-tuned T2-FPID controller is revealed with other control approaches reported in the literature by implementing test system 1. (c) Impact of CTDs on the MATHN system is revealed and the importance of considering CTDs is discussed extensively. (d) A UPFC-SMES additional strategy is implemented to further improve system performance. (e) Sensitivity analysis is performed to showcase the presented coordinated regulation robustness.

Materials and Methods
Two power system models are used in this work to implement the presented T2-FPID controller based on the WCA approach. A DAHT model named test system 1 with unique generation capacities is analyzed by laying area 1 with an SLP of 10%. Another model is a dual-area THN system referred to as test system 2 in this paper. The constraints of practicality, such as GRCs and CTDs, are considered for test system 2 to exercise the research close to the nature of realistic practice. The test system models 1 and 2 are depicted in Figure 1 [4] and Figure 2 [21], respectively. The capacity of the DAHT model is 2000 MW, with areas of unique generation capacities, and the MATHN system has a 2000 MW power generation capacity for each area. Investigation of the LFC of a nuclear plant along with a conventional thermal-hydro system is not available in the literature so far and is one of the key issues that is addressed in this paper. The parameters of both the test systems are detailed in Appendix A. Normally in the LFC of a power system, the information from the sensing devices located at distant terminal points is transmitted to the command control room, and from the command room, the control signals are transmitted to the generation plant locations to alter the set point so that the stability of the system can be maintained. This transmission and receiving of signals among different devices in distant places can be possible only through the usage of communication channels. A communication channel always inherits the features of delay. With these delays, the performance of the power plants is adversely affected. Until now, considering CTDs for the power system LFC analysis is not done satisfactorily. In this work, CTDs are carried out to investigate their impact on the power system and are modeled as [24] given in Equation (1), and the parameter τ d indicates the time delay constant. Tie-Line  -

Type II Fuzzy Controller
In a type II fuzzy controller, the limitations of a type I controller, such as exhibiting the crisp value from the incurred two-dimensional membership functions (MFs), are eliminated. A T2-FPID controller is implemented for the LFC study as a secondary regulator in this work, and the scaling factors are also perceived. The MFs used in the type II fuzzy systems are three-dimensional and are created with an upper MF and lower MF combination like that in a type 1 fuzzy system [25]. In this work, the ACE and the change in the ACE are given as input variables to the T2-FPID controller, and the fuzzifier unit converts them into a 3D representation and then transfers them to the rule-base unit and the defuzzifier unit. The membership functions (MFs) perceived in this work for both error and change in error are five linguistic variables termed big positive (BP), small positive (SP), zero (Z), big negative (BN), and small negative (SN), as depicted in Figure 3. A Mamdani type of fuzzy engine is perceived, and the FLC output is calculated by using the defuzzification method of the center of gravity. The two-dimensional FLC rule base is noted in Table 1. The structure of the T2-FPID controller is shown in Figure 4. The parameters K 1 and K 2 are scaling factors, and the gains of the PID are indicated with K P , K I , and K D , and N represents the filter coefficient.

Type II Fuzzy Controller
In a type II fuzzy controller, the limitations of a type I controller, such as exhibiting the crisp value from the incurred two-dimensional membership functions (MFs), are eliminated. A T2-FPID controller is implemented for the LFC study as a secondary regulator in this work, and the scaling factors are also perceived. The MFs used in the type II fuzzy systems are three-dimensional and are created with an upper MF and lower MF combination like that in a type 1 fuzzy system [25]. In this work, the ACE and the change in the ACE are given as input variables to the T2-FPID controller, and the fuzzifier unit converts them into a 3D representation and then transfers them to the rule-base unit and the defuzzifier unit. The membership functions (MFs) perceived in this work for both error and change in error are five linguistic variables termed big positive (BP), small positive (SP), zero (Z), big negative (BN), and small negative (SN), as depicted in Figure 3. A Mamdani type of fuzzy engine is perceived, and the FLC output is calculated by using the defuzzification method of the center of gravity. The two-dimensional FLC rule base is noted in Table 1. The structure of the T2-FPID controller is shown in Figure 4. The parameters K1 and K2 are scaling factors, and the gains of the PID are indicated with KP, KI, and KD, and N represents the filter coefficient.
However, parameters such as scaling factors, filter, and PID gains are optimally rendered with the WCA mechanism subjected to the time domain objective index of the ITAE given in Equation (2), where t represents the simulation time; parameters Δ f 1 and Δf 2 indicate variations in the frequency of area 1 and area 2, respectively; and ΔP tie12 represents deviations in the tie-line power flow.   given in Equation (2), where t represents the simulation time; parameters Δ f 1 and Δf 2 indicate variations in the frequency of area 1 and area 2, respectively; and ΔP tie12 represents deviations in the tie-line power flow. Type-II Fuzzy ACE U(s)  However, parameters such as scaling factors, filter, and PID gains are optimally rendered with the WCA mechanism subjected to the time domain objective index of the ITAE given in Equation (2), where t represents the simulation time; parameters ∆ f 1 and ∆ f 2 indicate variations in the frequency of area 1 and area 2, respectively; and ∆P tie12 represents deviations in the tie-line power flow.

Territorial Strategy of the SMES-UPFC
In the course of greater load disturbances, the secondary regulator is not passable enough to handle the system deviations. So, a territorial control strategy is implemented in this work in order to sustain system stability even in worst-case scenarios.

UPFC
Through an intended path, real and reactive power can be diverted by the UPFC. Moreover, they are multivariable FACTS devices that can regulate system parameters and thereby enhance the interline flow of power. The combination of an SSSC and a STATCOM with a common DC link forms the device UPFC. Without additional sources, the UPFC can facilitate real power compensation and also have the capability to absorb/generate real power and hence enhance system stability. The UPFC structure implemented in this work as a damping controller is shown in Figure 5 [21].

Territorial Strategy of the SMES-UPFC
In the course of greater load disturbances, the secondary regulator is not passable enough to handle the system deviations. So, a territorial control strategy is implemented in this work in order to sustain system stability even in worst-case scenarios.

UPFC
Through an intended path, real and reactive power can be diverted by the UPFC. Moreover, they are multivariable FACTS devices that can regulate system parameters and thereby enhance the interline flow of power. The combination of an SSSC and a STAT-COM with a common DC link forms the device UPFC. Without additional sources, the UPFC can facilitate real power compensation and also have the capability to absorb/generate real power and hence enhance system stability. The UPFC structure implemented in this work as a damping controller is shown in Figure 5 [21].

SMES
An SMES play a handy role in providing real power during line outages at the load end. These devices store energy in the form of a DC magnetic field by passing current through the conductor wound around the magnet, which is maintained at cryogenic tem- Figure 5. UPFC damping controller [21].

SMES
An SMES play a handy role in providing real power during line outages at the load end. These devices store energy in the form of a DC magnetic field by passing current through the conductor wound around the magnet, which is maintained at cryogenic temperature to impart a lossless nature. SMES devices get charged up during off-peak demands and discharged at the time of power shortages. They have the feasibility to enhance the line capacity by softening inter-area deviations. An SMES acts as an excellent spinning reserve to meet the sudden demands on the system, thereby diminishing frequency fluctuations. An SMES has the capability of huge storage capacity, and its low charge/discharge time makes its application possible for LFC study. The mathematical model of the SMES [26] used in work is given in Equation (3).

Water Cycle Algorithm
Eskandar et al. [27] in the year 2012 proposed this WCA strategy inspired from the continuous water movement. Since then, this algorithm has been widely used by researchers to solve non-linear constrained optimization problems and also engineering optimizations. In this water cycle process, water collects at the banks of streams and rivers. Later, these rivers/streams continue to join the sea. Moreover, this strategy possesses strength to obtain minimum or maximum index values accurately within a few iterations. So, in this work, this WCA approach is followed to retrieve the parameters of the secondary regulator optimally so that variations in frequency are reduced. In this WCA, rivers or streams act as the initial population and the sea is the global best.
Every raindrop (RD) is generated as an array of Nvar acts as a solution for the problem initialized as Later, every RD cost is evaluated by using the ITAE equation, and the gains of the T2-FPID are initialized, subjected to the constraints modeled as follows:

Streams Flow into the River or Rivers Flow into the Sea
The water in the streams tries to flow toward the river, and in the same way, the water in the rivers flows toward the sea. The position of streams and rivers is updated using the following expressions: P River (t + 1) = P River (t) + rand() × k × (P Sea (t) − P River (t)) Here, k representing the constant usually takes a value between 1 and 2, rand() is a random number from 0 to 1, and t indicates iterations. In case the objective index developed by the stream is better than that of the river, then the positions of streams and rivers are altered. A similar strategy is applicable for rivers and the sea.

Evaporation and Rain
This loop has a vital role in the searching mechanism of WCA. This evaporation and rain loop is intended to avoid getting trapped into the local optima. This loop keeps this WCA approach in a hierarchical position compared to other techniques. The process of evaporation terminates if P Sea − P River < d max (13) After the process of evaporation ends, immediately the process of rain starts: The searching rate nearer to the sea is indicated with V. The global best solutions are displayed immediately when the main loop reaches the maximum iterations. The procedural flow of WCA is depicted in Figure 6. Every rain drop cost is calculated using Equation (2) Flow of rivers and streams to the sea Equation (11&12) Is cost value of sea (or river) < the river (or stream) ?
Flow of new raindrops to streams in different locatiobs-Equation (15) Is Equation (13)  To notice the execution of the presented regulator WCA-optimized T2-FPID, a conventional hydro-thermal dual-area system is considered by perceiving area 1 with a disturbance of 10% SLP. However, controlling strategies reported in the early and recent lit-  To notice the execution of the presented regulator WCA-optimized T2-FPID, a conventional hydro-thermal dual-area system is considered by perceiving area 1 with a disturbance of 10% SLP. However, controlling strategies reported in the early and recent literature, such as PI, a PID tuned with the respective GA [4], and hFA-PS [8] algorithms along with a fuzzy PID optimized with hPSO-PS [22], are implemented to demonstrate WCA-based T2-FPID control efficacy. The system nature under these controllers is assessed and compared in Figure 7, and respective responses' settling times along with regulator optimal parametric gains are noted in Table 2. Gathering from Figure 7 and interpreting the settling times noted in Table 2, it is confirmed that the presented WCA-based T2-FPID shows remarkable performance in damping out frequency and inter-area power flow deviations. Hence, the presented WCA-based T2-FPID approach is more dominating and outperforms the methodologies that are listed in the literature because of the inheritance of the evaporation and rain loop in the WCA mechanism to avoid solution divergence and because the usage of 3D MFs in the T2-FPID regulator enhances the degree of freedom. ulator optimal parametric gains are noted in Table 2. Gathering from Figure 7 and interpreting the settling times noted in Table 2, it is confirmed that the presented WCA-based T2-FPID shows remarkable performance in damping out frequency and inter-area power flow deviations. Hence, the presented WCA-based T2-FPID approach is more dominating and outperforms the methodologies that are listed in the literature because of the inheritance of the evaporation and rain loop in the WCA mechanism to avoid solution divergence and because the usage of 3D MFs in the T2-FPID regulator enhances the degree of freedom.

Case 2: Performance Analysis of Test System 2 without CTD Consideration
Later, the presented T2-FPID regulator is also implemented in the MATHN model, named test system 2 in this paper, depicted in Figure 2. The PID, fuzzy PID, and T2-FIPD controllers are used as secondary regulators optimized with the WCA mechanism one after the other, and the behavior of the system is assessed upon laying area 1 with a disruption of 10% SLP. A performance comparison of these controllers on the MATHN system without taking CTDs into consideration is shown in Figure 8. Visualizing Figure 8, the T2-FPID is more dominating than a conventional PID and fuzzy PID in dampening deviations and also providing the responses to reach a steady condition.

Case 2: Performance Analysis of Test System 2 without CTD Consideration
Later, the presented T2-FPID regulator is also implemented in the MATHN model, named test system 2 in this paper, depicted in Figure 2. The PID, fuzzy PID, and T2-FIPD controllers are used as secondary regulators optimized with the WCA mechanism one after the other, and the behavior of the system is assessed upon laying area 1 with a disruption of 10% SLP. A performance comparison of these controllers on the MATHN system without taking CTDs into consideration is shown in Figure 8. Visualizing Figure 8, the T2-FPID is more dominating than a conventional PID and fuzzy PID in dampening deviations and also providing the responses to reach a steady condition.

Case 3: Performance Analysis of Test System 2 with CTD Consideration
In this case, the MATHN model is used with practicality constraints of CTDs and the behavior of the system is monitored under the supervision of a PID, fuzzy PID, and T2-FPID based on the WCA mechanism one at a time. Figure 9 depicts the comparative analysis of these controllers on the MATHN system with CTDs and reveals that even when the system is considered with CTDs, the presented T2-FPID outperforms other controllers. Optimal parameters of controllers and numerical results of the system responses are noted in Tables 3 and 4.     To showcase the impingement of CTDs on MATHN system performance, MATHN system responses under the governing of the WCA-optimized T2-FPID regulator are compared in Figure 10, without and with adopting CTDs. From the above cases, as the T2-FPID controller is proven to be the most dominating, only the system responses with this controller have been considered in this case. From Figure 10, CTDs have a considerable effect on system stability. Reponses are much deviated upon considering CTDs with the system compared to the case of not taking them into account. With the delay in communication between sensors situated at distant locations and power generation units, altering of the turbine set point valve position will get delayed. Due to this, there will be widening of the real power mismatch, leading to more variations in system dynamics. However, CTDs are needed to be taken into account while assessing the performance of interconnected models. While doing so, the designed regulator can handle the system deviations more effectively compared to the case of a controller designed for a system not taking CTDs into account.

Case 5: Performance Assessment of the UPFC-SMES on Test System 2
Later, the mechanism of a territorial control scheme is implemented for test system 2 under governing of the WCA-tuned T2-FPID regulator to improve the system dynamical nature. First, only SMES devices are merged in each area and the system performance is assessed. Later, a UPFC device is placed in the inter-area along with SMES device positions, which are kept unchanged, and the responses are shown in Figure 11. As shown in Figure 11, there is an enhancement in the performance of the MATHN model with the UPFC-SMES strategy. This is possible only because of the nature of the SMES in meeting a sudden load demand and the UPFC in regulating the inter-area power flow.

Case 6: Sensitivity Analysis
Further, the MATHN model is subjected to different load variations and the system parameter of the tie-line coefficient in the range of ±50% from the nominal value to validate the robustness of secondary and territorial coordinated control mechanisms. The responses for testing the robustness are shown in Figures 12 and 13. From this case, it is finalized Energies 2021, 14, 5387 13 of 19 that the system dynamical behavior is not much affected, even if the system is subjected to parameter deviations and load variations. under governing of the WCA-tuned T2-FPID regulator to improve the system dynamical nature. First, only SMES devices are merged in each area and the system performance is assessed. Later, a UPFC device is placed in the inter-area along with SMES device positions, which are kept unchanged, and the responses are shown in Figure 11. As shown in Figure 11, there is an enhancement in the performance of the MATHN model with the UPFC-SMES strategy. This is possible only because of the nature of the SMES in meeting a sudden load demand and the UPFC in regulating the inter-area power flow.

Case 6: Sensitivity Analysis
Further, the MATHN model is subjected to different load variations and the system parameter of the tie-line coefficient in the range of ±50% from the nominal value to vali- Further, the MATHN model is subjected to different load variations and the system parameter of the tie-line coefficient in the range of ±50% from the nominal value to validate the robustness of secondary and territorial coordinated control mechanisms. The responses for testing the robustness are shown in Figures 12 and 13. From this case, it is finalized that the system dynamical behavior is not much affected, even if the system is subjected to parameter deviations and load variations.

Conclusions
This paper addressed the potentiality of a WCA-optimized T2-FPID controller in mitigating the frequency and interline flow of power fluctuations in the MATHN power system. However, the efficacy of the T2-FPID is revealed by comparing its performance with a fuzzy PID and a conventional PID. Further, the performance of the MATHN power system is scoped by considering it without and with CTDs. The presented T2-FPID regulator enhances the performance index value by 35.6% with a fuzzy PID and 42.2% with a traditional PID in the case of not considering CTDs and by 23.3% with a fuzzy PID and 31.4% with a traditional PID in the case of considering CTDs. The impingement of CTDs on the MATHN system is showcased and justified. Subsequently, implementing the territory control of the UPFC-SMES scheme under the governing of a WCA-tuned type II fuzzy regulator shows remarkable enhancement in mitigation of variations. Later, a robustness test is performed and it reveals the robustness of the presented secondary and territory regulation schemes. Moreover, the presented WCA-based type II fuzzy regulator is applied on the widely accepted test system of DAHT and the presented approach outperforms other controllers such as the PI and PID tuned with GA and hFA-PS techniques, along with the fuzzy PID tuned with hPSO-PS mentioned in the recent literature.