Optimum Synthesis of a BOA Optimized Novel Dual-Stage PI − (1 + ID) Controller for Frequency Response of a Microgrid

: A renewable and distributed generation (DG)-enabled modern electriﬁed power network with / without energy storage (ES) helps the progress of microgrid development. Frequency regulation is a signiﬁcant scheme to improve the dynamic response quality of the microgrid under unknown disturbances. This paper established a maiden load frequency regulation of a wind-driven generator ( WG ), solar tower ( ST ), bio-diesel power generator ( BDPG ) and thermostatically controllable load (heat pump and refrigerator)-based, isolated, single-area microgrid system. Hence, intelligent control strategies are important for this issue. A newly developed butterﬂy algorithmic technique (BOA) is leveraged to tune the controllers’ parameters. However, to attain a proper balance between net power generation and load power, a dual stage proportional-integral- one plus integral-derivative PI − (1 + ID) controller is developed. Comparative system responses (in MATLAB / SIMULINK software) for di ﬀ erent scenarios under several controllers, such as a proportional-integral ( PI ), proportional-integral-derivative ( PID ) and PI − (1 + ID) controller tuned by particle swarm optimization (PSO), grasshopper algorithmic technique (GOA) and BOA, show the superiority of BOA in terms of minimizing the peak deviations and better frequency regulation of the system. Real recorded wind data are considered to authenticate the control approach. Validation, A.L.; Formal Analysis, A.L., S.M.S.H., D.C.D. and T.S.U..; Writing—Original Draft Preparation, A.L. and D.C.D.; Writing—Review and Editing, S.M.S.H., and T.S.U.; Visualization, A.L.; Supervision, D.C.D.; Funding Acquisition, T.S.U. All authors


Introduction
Electricity consumption is increasing in parallel with population and energy demand [1]. The increasing generation capacity with conventional energy sources has negative impacts on environment [2]. Research shows that the introduction of microgrids with renewable energy resources (RERs) is an environmentally friendly solution to this energy problem [3,4].
Microgrids can provide energy in a clean and optimal way when digital control technologies are coupled with sources such as wind and solar. However, the intermittent characteristics of RERs and low inertia of inverter-interfaced systems cause control and stability issues in microgrids [5]. Hence, combining diesel generators and RERs [6][7][8] is one possible solution to mitigate the detrimental effects (e.g., frequency fluctuation) of hybrid power systems. Additionally, a non-toxic bio-diesel power generator (BDPG) can be a more environment-conscious supplementary option for frequency regulation schemes. 2 of 12 To overcome the frequency fluctuation in a more reliable way, different storage devices (SDs) have been considered, such as battery (BSD), fuel cell (FCSD), ultra-capacitor (UCSD) and superconducting magnetic storage system (SMSD) [9]. There are maintenance and disposal concerns for BSD, while FCSD suffers from slow response and SMSU experiences leakage of expensive helium liquid [10]. In the following, a cost-effective, carbon-neutral-based, solid-oxide fuel cell (SOFC) can be utilized for the frequency regulation of isolated hybrid microgrids (IHµGS). In addition, different thermostatically controllable loads (TCLs) such as a heat pump (HP) and refrigerator (RFZ) are employed for smoothing the dynamic system responses. In the recent past, several studies have focused on the frequency regulation of IHµGS [11][12][13][14][15]. The study in [12] framed out a mathematical modeling of a system comprised of dish-stirling, solar thermal, diesel and SDs. The authors of [13] sketched a load frequency management for a hybrid power system that includes wind, solar PV and SDs. In fact, the application of a plug-in electric vehicle (PHEV)-battery (BSD) to contain frequency fluctuation is investigated in [14].
Several algorithmic techniques have been leveraged in order to optimally tune the controller parameters for IHµGSs, such as genetic algorithm (GA) [9], PSO [6], firefly technique (FA) [19], cuckoo search technique (CS) [14], mine blast technique (MBA) [20], grasshopper algorithmic technique (GOA) [21]. A comprehensive review of different algorithmic techniques for a load frequency controller is presented in [22]. In this regard, this work explores the application of a butterfly algorithmic tool (BOA) for designing the parameters of a frequency controller. This algorithm was recently developed and is considered in this paper due to its high convergence rate [23].
Therefore, the scope of this work and its contributions to the current body of knowledge can be summarized as follows: controllers under a BOA algorithmic tool; (d) Comparative system dynamic analysis of different algorithms (PSO, GOA and BOA), leveraging the acquired superior controller in (c); (e) Study system dynamics under real recorded wind data and other random disturbances.
The rest of the paper is organized as follows: Section 2 details the frequency response modeling steps. Section 3 gives an overview of the BOA technique and shows its adaptation for the purposes of this work. It also presents the proposed dual-stage controller. Simulation works and their analyses are presented in Section 4. The conclusions are given in Section 5.

Frequency Response Modeling of the Proposed Dual-Stage Controller
The hybrid system of the proposed work consists of wind generators (1.5 MW); solar-tower-based, solar-thermal power system (1 MW), BDPG (800 kW), SOFC (200 kW); thermostatically controllable HP; RFZ elements; and demanded loads (2.2 MW). The schematic layout and abbreviation of the relevant parameters are shown in Figure 1 and Table 1, respectively.

Solar Tower (ST)
The dual-axis (vertical and horizontal) heliostats-enabled central receiver system is ST, placed on a surface. Here, the reflected solar radiation is focused on the central receiver of ST with a higher concentration ratio (500-1000) and temperature (500-850 0 C) of fluid (steam or molten salt). The collected incident solar power (Pin) can be deliberated as . .
where Ah is the heliostats area, I incident solar radiation and ηh is the system constant. By solving the state equations, the linearized transfer function model of ST can be represented as [18] ( ) . .

Wind Generator (WG)
The kinetic energy of the wind converts into electrical energy through the wind generator (WG). As wind is a highly variable source, the power output through the WG depends on the instantaneous speed of the wind. Equation (1) formulates how wind energy is converted to the mechanical output power of WG.
where ρ, V WG , A bd , and C P are, in proper order, the air density, intermittent wind speed, blade-swept area and the extractable power co-efficient. Akkanayakanpatti station's recorded wind speed data are Energies 2020, 13, 3446 4 of 12 considered and modelled in the proposed work [24]. The rate of change in real recorded wind power (∆P WG ) and transfer function model of the WG are represented as shown in Equation (2) [25] (2)

Solar Tower (ST)
The dual-axis (vertical and horizontal) heliostats-enabled central receiver system is ST, placed on a surface. Here, the reflected solar radiation is focused on the central receiver of ST with a higher concentration ratio (500-1000) and temperature (500-850 • C) of fluid (steam or molten salt). The collected incident solar power (P in ) can be deliberated as where A h is the heliostats area, I incident solar radiation and η h is the system constant. By solving the state equations, the linearized transfer function model of ST can be represented as [18]

Biodiesel Power Generator (BDPG)
The combination valve regulator and combustion-engine-based biodiesel power generation (BDPG) was leveraged to offer support as a backup power generation. It has inherently biodegradable and non-toxic positive characteristics, which were the main reasons to incorporate it in the suggested work. Equation (5) details the transfer function of BDPG [21].

Solid-Oxide-Based Fuel Cell (SOFC)
Through the electrochemical reaction, the fuel cell produces dc power and, by using a DC-AC converter, this power is converted into AC. With a fast charging-discharging time and higher efficiency (~80%), the SOFC has gained much interest in recent years among all the categories of FCSDs. In view of the above, SOFC was selected as the storage device in the system. Its transfer function model is given in Equation (6) [18]

Thermostatically Controllable Loads (HP and RFZ)
In order to manage the energy consumption and to improve the system dynamic responses (by controlling operation cycles), two thermostatically controllable loads were considered, i.e., heat pump (HP) and RFZ. The transfer function models of HP [19] and RFZ [19] could be expressed as in Equations (7) and (8) G HP (s) = K HP sT HP + 1 G RFZ (s) = K RFZ sT RFZ + 1 Energies 2020, 13, 3446 5 of 12

IHµGS Dynamic Model
The instantaneous change in power (∆P G ) of the proposed IHµGS can be formulated as where The equivalent dynamic model of IHµGS could be illustrated as Refer to Table 1 for the nomenclature and abbreviations used for IHµGS system modelling.

Objective Function Formulation
The formulation of the objective function (J) has a great impact on system dynamics and the achieved results. Therefore, the proposed work considered the integral of square error (ISE) objective function. This could be formulated as Sub ject to : where i = 1, 2. The range of controller parameters is taken as (0-50).

Optimization Techniques
Three metaheuristic techniques were considered to optimally tune the controller parameters along with their comparative dynamic responses.

Particle Swarm Technique (PSO)
A swarm-based metaheuristic particle swarm technique (PSO) was developed by Eberhart and Kennedy in 1995 to solve the specified problem by improving the candidate solution with reference to the given measure quality [26]. The solution of PSO is termed as a particle. Every particle follows a track of coordinates in the problem space until the best solution is reached with respect to the suggested problem. The velocity of each particle is varied on the basis of best position (P best ) and I best location. The optimum values obtained by optimizer are termed as I best [26].

Grasshopper Algorithmic Technique (GOA)
A metaheuristic grasshopper algorithmic technique was proposed by Saremi et al. [27]. Its characteristics depend on the swarming and foraging characteristics of grasshopper, which could be modeled to form structural algorithmic techniques. The steps involved for initialization, exploitation and exploration are depicted in [27].

Butterfly Optimization Technique (BOA) and Proposed Dual-Stage Controller
Recently, since 2018, based on the food-probing approach and mating characteristics of butterflies, a natured-inspired algorithmic technique named BOA has been developed to solve several engineering problems. The unique food-probing strategy and mating characteristics of BOA are modeled in [23].
The main idea of BOA depends on three key parameters. These are sensor modality (M s ), impulsive intensity (I s ) and power component (γ). Moreover, the objective function of this technique depends on the distinction of I s and formation of fragrance (P), which could be formulated as where P i is the fragrance magnitude of ith butterfly. In the following, to investigate global search stage, a dominated fitted solution q* could be depicted as where m and q* are the solution vectors of ith butterfly and current best solution among all the solutions. The formulation of the social search could be illustrated as where, m g and m h are the gth and hth butterflies enabled in the search space [23]. The approximate flow diagram of BOA technique is framed out in Figure 2. All the parameters considered for algorithmic techniques are given in Appendix A.
Energies 2020, 13, x FOR PEER REVIEW 7 of 12 To reach above control target, a novel dual-stage proportional-integral-integral-derivative (PI-(1+ID)) controller is deployed, as framed in Figure 3. In the following, Δf is leveraged as an input signal, whereas C(s) is the output control signal of the controller. The control output signal and transfer function of the proposed controller are formulated as  To reach above control target, a novel dual-stage proportional-integral-integral-derivative (PI − (1 + ID)) controller is deployed, as framed in Figure 3. In the following, ∆f is leveraged as an input signal, whereas C(s) is the output control signal of the controller. The control output signal and transfer function of the proposed controller are formulated as To reach above control target, a novel dual-stage proportional-integral-integral-derivative (PI-(1+ID)) controller is deployed, as framed in Figure 3. In the following, Δf is leveraged as an input signal, whereas C(s) is the output control signal of the controller. The control output signal and transfer function of the proposed controller are formulated as ( ) . (

Frequency Response Studies and Analysis
In order to verify the proposed control strategy, two scenarios are simulated in a system with Core-i7-4770 CPU under MATLAB/SIMULINK (R2013a, MathWorks, Natick, USA) was environment. Three algorithmic techniques (PSO, GOA and BOA) have been considered. Furthermore, to validate the control strategy, real recorded wind speed data have been considered.

Scenario 1: Performance Analysis of All Controllers during Non-Accessibility of All RERs
In this scenario, assume that all the RERs are unavailable due to maintenance. Therefore, the extractable power forms WG (∆P WG ) and ST (∆P ST ) are zero (∆P WG = ∆P ST = 0%) during the entire period. A net constant critical load demand (∆P CL = 30%) is considered from t = 0 s onwards. The comparative performance of different controllers such as PI, PID and (PI − (1 + ID)) are displayed in Figure 4, where the tuned parameters are listed in Table 2. The system dynamics assessment of the abovementioned controllers under BOA and objective function (J ISE ) and figure of demerits (J FOD ) clearly depicts that the proposed (PI − (1 + ID)) controller is superior to the rest. To elaborate further, the performance indicators such as peak overshoot (+O P ), peak undershoot (-U P ) and settling time (T ST ) are tabulated in Table 2.   In this scenario, the proposed system is tested under real-recorded wind (obtained from National Institute of Wind Energy, India) [25], as displayed in Appendix B. The operating condition

Scenario 2: Performance analysis of Different Algorithms Under Concurrent Random Changes of WG (Utilization of Real-Recorded Data), ST and Critical Load Demand
In this scenario, the proposed system is tested under real-recorded wind (obtained from National Institute of Wind Energy, India) [25], as displayed in Appendix B. The operating condition is illustrated with 30% average power of ST and 50% critical load demand for the entire time duration. The wind speed and its corresponding output power is shown in Figure 5a. The results are depicted in Figure 5b-e, showing the comparative system dynamic responses of ∆f, ∆P BDPG , ∆P SOFC , ∆P HP , and ∆P RFZ . Figure 5b-e it clearly shows that the BOA-optimized PI − (1 + ID) controller performed better than the other suggested PSOs, GOA-tuned (PI − (1 + ID)) controller. The tuned values of the controller parameters are displayed in Table 3.

Conclusion
The present article develops a novel frequency regulation scheme for wind-solar-towerbiodiesel-based IHμGS. A novel dual-stage PI-(1+ID) controller is enabled to investigate the system dynamics under different scenarios. A recently developed BOA technique is utilized to optimally tune the proposed dual stage PI-(1+ID) controller gains and compare the system dynamics under real

Conclusions
The present article develops a novel frequency regulation scheme for wind-solar-tower-biodieselbased IHµGS. A novel dual-stage (PI − (1 + ID)) controller is enabled to investigate the system dynamics under different scenarios. A recently developed BOA technique is utilized to optimally tune the proposed dual stage (PI − (1 + ID)) controller gains and compare the system dynamics under real recorded wind data. The comparative system dynamic responses, as well as performance parameters such as peak deviation (+O P , −U P ) and settling time (T ST ), clearly indicate that the BOA-optimized (PI − (1 + ID)) controller performs better than other classical benchmark controllers. The simulation test results prove the effectiveness of the proposed control strategy. This control scheme could be further extended by integrating different RER technologies and storage devices, as well as electric vehicles, into the microgrid.