Optimal Sizing and Power System Control of Hybrid Solar PV-Biogas Generator with Energy Storage System Power Plant

: In this paper, the electrical parameters of a hybrid power system made of hybrid renewable energy sources (HRES) generation are primarily discussed. The main components of HRES with energy storage (ES) systems are the resources coordinated with multiple photovoltaic (PV) cell units, a biogas generator, and multiple ES systems, including superconducting magnetic energy storage (SMES) and pumped hydro energy storage (PHES). The performance characteristics of the HRES are determined by the constant power generation from various sources, as well as the shifting load perturbations. Constant power generation from a variety of sources, as well as shifting load perturbations, dictate the HRES’s performance characteristics. As a result of the ﬂuctuating load demand, there will be steady generation but also ﬂuctuating frequency and power. A suitable control strategy is therefore needed to overcome the frequency and power deviations under the aforementioned load demand and generation conditions. An integration in the environment of fractional order (FO) calculus for proportion-al-integral-derivative (PID) controllers and fuzzy controllers, referred to as FO-Fuzzy-PID controllers, tuned with the opposition-based whale optimization algorithm (OWOA), and compared with QOHSA, TBLOA, and PSO has been proposed to control the frequency deviation and power deviations in each power generation unites. The results of the frequency deviation obtained by using FO-fuzzy-PID controllers with OWOA tuned are 1.05%, 2.01%, and 2.73% lower than when QOHSA, TBLOA, and PSO have been used to tune, respectively. Through this analysis, the algorithm’s efﬁciency is determined. Sensitivity studies are also carried out to demonstrate the robustness of the technique under consideration in relation to changes in the sizes of the HRES and ES system parameters.


Introduction
Energy is one of the industries that is growing the fastest because it has to deal with a number of problems, such as traditional energy sources running out, global warming getting worse, and other problems.Because of these problems, it is more important than ever to use non-traditional or renewable energy sources, such as solar, wind, and biogas [1,2].There are a lot of these resources, but they are very random because they depend on the weather at any given time.There may be occasions when the sources' generation is not enough to satisfy the load demand at a particular moment in time.Batteries, ultra-capacitors, and flywheels might all be connected to hybrid power systems to correct this imbalance.These devices improve the quality of electricity while also lowering grid frequency variations [3,4].
Surplus power is available when the generation exceeds the load requirements, which most often occurs during peak weather conditions.At this time, storage devices temporarily store energy.This stored energy is used when the generated energy falls short of the rise in load demand.In the coming years, there will be a greater use of renewable energy sources, such as superconducting magnetic and pumped hydro energy storage systems, to reduce fuel prices, reduce CO 2 emissions, and improve network power quality [5][6][7].The load frequency control scale can be used to fix imbalances caused by changes in frequency caused by the fact that wind and solar energy are not always available [8,9].In a standby reserve, extra power is saved and can be used when the demand on the grid changes.The SMES, which are incorporated in the hybrid power system renewable energy generation analysis, have the same purpose of regulating the grid's variation in load as they are plugged in.Control strategies for coordination are required for proper performance.In [10], general classical controllers are utilized.The authors of [11] proposed the evolutionary optimization technique for a comparative analysis of various conventional controllers for an isolated wind-diesel hybrid power plant.In [12], the authors proposed employing linear matrix inequalities-based robust controller design and particle swarm optimization to manage frequency in hybrid power systems.The authors of [13] employed particle swarm optimization to control and simulate a wind-biomass isolated hybrid power system.For frequency regulation of hybrid power systems with electric vehicles, the authors in [14] developed the fractional-order adaptive fuzzy PID controller with modified salp swarm algorithm optimization.The fractional order fuzzy PID controller for frequency regulation of a solar-wind integrated power system with a hydrogen aqua equalizer-fuel cell unit was presented by the authors in [15].In [16], the authors suggested developing a fuzzy logic controller for load frequency management in an isolated hybrid power system.A type 2 fuzzy PID controller for frequency control in hybrid distributed power systems was presented by the authors in [17].A model of an island hybrid power system that uses artificial bee colony optimization for fuzzy logic-based load frequency management was described by the authors in [18].The authors in [19] suggested a novel scaling factorbased fuzzy logic controller for frequency regulation of an isolated hybrid power system.The authors of [10] provided integrated frequency and power regulation of a freestanding hybrid power system while considering a fuzzy classical controller with a scaling factor.The authors of [20] described frequency control on a remote island using parallel battery systems and the H control theory based on droop characteristics.The authors of [21] introduced the dependable H-infinity load frequency control in a hybrid distributed generating system.By incorporating SVC into power system networks, the writers in [22] explained how to employ the best load flow technique in MATLAB/Simulink settings to enhance voltage profiles.Reference [23] examined the best location for static VAr compensators in electric power systems.The authors discussed the integration of static variables to account for hybrid power systems and renewable-energy-producing systems in [24].
Fractional calculus [25,26] has been used in control systems and other domains in recent years.Modern contributions to this field include control systems that use computational intelligence [27] with arbitrary time delays, such as process control [28], reactive power control [29], a cascade fuzzy-fraction order (FO)-proportional-integral-derivative (PID) controller for automated generation control (AGC) [30], and using an FO-PID with a filter for AGC in linked thermal and thermal-hydro-gas power systems [31].In addition, ultra-capacitor effects on AGC in multi-area integrated thermal and thermal-hydro-gas power systems by utilizing a fuzzy-FO-integral-derivative controller and a PID with a filter [32], are recent.In integrated PV, thermal, and hydrothermal power plants, load frequency control (LFC) with fuzzy FO-integral-derivative with filter [33], fuzzy PID-filter-(1 + PI) controller in AGC PV, thermal, and hydrothermal power plants [34], and others.In this paper, commonly used implementations of traditional PID, fuzzy PID, FO-PID, and FO-fuzzy PID controllers are used to show the benefits of the FO-fuzzy PID controller system.
As already discussed, solar energy is mostly intermittent, this means that the frequency of both the proposed hybrid system and the country's grid is always changing.Installing a controller with a power loop, sends or receives signals to or from the grid to energy storage devices to store or release extra or not enough power, affects the power quality.The controller may even instruct the biogas generator to transfer unusually large bursts of electricity into the national grid in order to prevent the failure of the load to which it is connected.Even better, the suggested control architecture puts all of the hybrid power system's renewable energy production under a single set of rules.Centralization also removes the need for a separate controller and tuning for power storage devices.The study area is illustrated as shown in Figure 1 by using quantum geographic information system software (qgis software version 3.8).
power control [29], a cascade fuzzy-fraction order (FO)-proportional-integral-derivative (PID) controller for automated generation control (AGC) [30], and using an FO-PID with a filter for AGC in linked thermal and thermal-hydro-gas power systems [31].In addition, ultra-capacitor effects on AGC in multi-area integrated thermal and thermal-hydro-gas power systems by utilizing a fuzzy-FO-integral-derivative controller and a PID with a filter [32], are recent.In integrated PV, thermal, and hydrothermal power plants, load frequency control (LFC) with fuzzy FO-integral-derivative with filter [33], fuzzy PID-filter-(1 + PI) controller in AGC PV, thermal, and hydrothermal power plants [34], and others.In this paper, commonly used implementations of traditional PID, fuzzy PID, FO-PID, and FO-fuzzy PID controllers are used to show the benefits of the FO-fuzzy PID controller system.
As already discussed, solar energy is mostly intermittent, this means that the frequency of both the proposed hybrid system and the country's grid is always changing.Installing a controller with a power loop, sends or receives signals to or from the grid to energy storage devices to store or release extra or not enough power, affects the power quality.The controller may even instruct the biogas generator to transfer unusually large bursts of electricity into the national grid in order to prevent the failure of the load to which it is connected.Even better, the suggested control architecture puts all of the hybrid power system's renewable energy production under a single set of rules.Centralization also removes the need for a separate controller and tuning for power storage devices.The study area is illustrated as shown in Figure 1 by using quantum geographic information system software (qgis software version 3.8).The present study proposes an opposition-based improved whale optimization algorithm (WOA), also known as an opposition-based along whale optimization algorithm.In comparison to other meta-heuristic algorithms, WOA is a genuine parameter optimization approach that does not use derivatives and has fewer mathematical requirements for compression.The approach could also be used to figure out how big hybrid renewable-energy-generating systems should be devised, whether connected to the grid or not.This is just one of many engineering optimization problems that the approach could help with.
The WOA metaheuristic optimization method does a better job of quickly finding high-performance zones in single and multi-objective functions cases in hybrid renewa- The present study proposes an opposition-based improved whale optimization algorithm (WOA), also known as an opposition-based along whale optimization algorithm.In comparison to other meta-heuristic algorithms, WOA is a genuine parameter optimization approach that does not use derivatives and has fewer mathematical requirements for compression.The approach could also be used to figure out how big hybrid renewableenergy-generating systems should be devised, whether connected to the grid or not.This is just one of many engineering optimization problems that the approach could help with.
The WOA metaheuristic optimization method does a better job of quickly finding high-performance zones in single and multi-objective functions cases in hybrid renewable energy sources [35].Banerjee et al. [36] used a different method based on the HS algorithm to change the reactive power in a model of an autonomous hybrid power system.To improve performance, a new fractional order (FO) fuzzy control method is used, and its parameters are changed with the help of a particle swarm optimization (PSO) method and two chaotic maps [37].
The characteristics of LFC for renewable energy generation in hybrid power systems are thought to depend on both the optimization method and the configuration of the controller [38].This paper suggests the opposition-based whale optimization algorithm (OWOA) to optimize the gains of the controller.Therefore, this work proposes an OWOA method based on FO-fuzzy-PID control for frequency and power regulation of hybrid power systems for producing renewable energy.The results of the simulation show that, unlike PID, FO-PID, fuzzy-PID, and FO-fuzzy-PID controllers, the controller under consideration can remove frequency differences well, improve supplies and provide reliable electric power going forward.In light of the aforementioned, the following elements serve as the main sources of inspiration for this study:

•
Modeling the proposed HRES with an ESS system enables businesses and organizations such as hospitals, universities, and other healthcare facilities to deliver dependable power while meeting rising demand.

•
The applicability of fast-response energy storage devices such as superconducting magnetic energy storage (SMES) systems and long-sustaining energy storage systems such as pumped hydro energy storage (PHES) systems for supplying quality and reliable power to meet demand can be investigated.

•
The HRES with the ESS model's decision-making process is enhanced through a recently developed variation of the OWOA that is suitably tuned to manage both frequency and power variations in its major adjustable parameters.
The current study focuses on using a FO-fuzzy-PID controller with an oppositionbased whale optimization technique to reduce the frequency and power oscillations that a load perturbation causes in the HRES with ES system.Based on the previous interaction, the main contributions of this research are as follows:

•
The model of an HRES with ES system, which consists of two solar photovoltaic (PV) units, two biogas-generating units, two PHES units, and two SMES units, is set up in this work to demonstrate how it can be constructed.

•
The FO-fuzzy-PID controller's design parameters were optimized using the OWOA metaheuristic technique.

•
The proposed HRES with ES system frequency and power deviation were investigated by using PID, FO-PID, fuzzy-PID and FO-fuzzy-PID with different metaheuristic optimization techniques (such as OWOA, QOHSA, TLBOA and PSO) for tuned controller parameters.

•
Disturbances were investigated bay considering the variation of connected loads, HRES sources and both of two system unites.

•
The OWOA metaheuristic optimization method was adopted and used to fine-tune the different parameters that could be changed on the controllers that were being studied.

•
The OWOA method was used to tune the model's parameters that could be changed.
It was then compared to other optimization methods such as QOHSA, TLBOA, and PSO to see their effectiveness.
The completed work is organized as follows: Section 2 discusses the case study area, Section 3 presents how the proposed system will work, and Section 4 describes and sets up the proposed HRES with ES system.The FO-fuzzy-PID controller is introduced briefly in Section 5. Section 5 explains the problem formulation, objective function, and constraints for optimization.Section 6 presents the results and a discussion for structures in relation to variations in power output caused by load demand.The references come after Section 7, which is the conclusion of the paper.

Case Study Description
The developed HRES with ES system is used to power Debre Markos University, Debre Markos, Ethiopia (latitude 10.33 0 N and longitude 37.71 0 E).The HRES with ES system is installed to provide the required daily energy demand in the Debre Markos University.The site shown in Figure 1 has an inspiring renewable resource potential [39] that justifies its use for energizing the university, especially given its maximum interruption duration from the national gird.

Methodology
The proposed system is a hybrid system that is connected to the grid and can obtain power from solar PV, a biogas generator, a pumped hydro energy storage system, or a superconducting magnetic energy storage system.A first look was taken at the village's energy needs, the energy infrastructure that could be used, and the current and future renewable projects in the area.The architecture of the system is set based on the local economy and energy needs, and the best size for each power generation technology is figured out.However, the most important data, such as how much electricity was used, how much water was available for PHES, and how much biogas waste was made, came directly from the community.Solar radiation and temperature, on the other hand, are real-time meteorological data.

Connected Load Assessment
A case study is then conducted on one such remote hamlet in Debre Markos, Amhara, Ethiopia.The grid-connected hybrid renewable energy system (HRES) case study was developed to suit a community's electricity demands.Figure 2 depicts the load profile of the community on a typical worst-case-scenario day.The connected load profile is slightly increased from 1668.3000 kW to 1707.4031 kW (from 08:00-09:00 to 14:00-15:00 h), and the maximum load demand is connected during this period.Between 07:00 and 08:00, the electricity demand falls to 662.6801 kW.According to the graph below, the minimum load needs are connected to the national grid around 07:00-08:00 h in the morning.
Debre Markos, Ethiopia (latitude 10.33′0′' N and longitude 37.71′0′' E).The HRES with ES system is installed to provide the required daily energy demand in the Debre Markos University.The site shown in Figure 1 has an inspiring renewable resource potential [39] that justifies its use for energizing the university, especially given its maximum interruption duration from the national gird.

Methodology
The proposed system is a hybrid system that is connected to the grid and can obtain power from solar PV, a biogas generator, a pumped hydro energy storage system, or a superconducting magnetic energy storage system.A first look was taken at the village's energy needs, the energy infrastructure that could be used, and the current and future renewable projects in the area.The architecture of the system is set based on the local economy and energy needs, and the best size for each power generation technology is figured out.However, the most important data, such as how much electricity was used, how much water was available for PHES, and how much biogas waste was made, came directly from the community.Solar radiation and temperature, on the other hand, are real-time meteorological data.

Connected Load Assessment
A case study is then conducted on one such remote hamlet in Debre Markos, Amhara, Ethiopia.The grid-connected hybrid renewable energy system (HRES) case study was developed to suit a community's electricity demands.Figure 2 depicts the load profile of the community on a typical worst-case-scenario day.The connected load profile is slightly increased from 1668.3000 kW to 1707.4031 kW (from 08:00-09:00 to 14:00-15:00 h), and the maximum load demand is connected during this period.Between 07:00 and 08:00, the electricity demand falls to 662.6801 kW.According to the graph below, the minimum load needs are connected to the national grid around 07:00-08:00 h in the morning.

Resources Assessment
NASA's satellite meteorology and solar energy databases were used to research the temperature and amount of solar radiation in the chosen area.The monthly ambient temperature at the chosen site ranges from 16.24 to 21.46 °C.The monthly solar radiation ranges from 5.6011 to 7.8010 kWh/m 2 .Figure 3 shows the monthly solar radiation and temperature values of the selected site.Figure 4 shows how the amount of daily solar radiation changes each month and how clear the area is, which suggests that the chosen

Resources Assessment
NASA's satellite meteorology and solar energy databases were used to research the temperature and amount of solar radiation in the chosen area.The monthly ambient temperature at the chosen site ranges from 16.24 to 21.46 • C. The monthly solar radiation ranges from 5.6011 to 7.8010 kWh/m 2 .Figure 3 shows the monthly solar radiation and temperature values of the selected site.Figure 4 shows how the amount of daily solar radiation changes each month and how clear the area is, which suggests that the chosen site for the PV system has sufficient solar capacity.Many researchers investigated whether the selected site has sufficient solar radiation to generate solar power [39,40].
The site's potential for producing biogas was calculated using primary and secondary data received from various sources [41].Based on the data acquired, the biogas potential of the research area is estimated to be 7.9287 × 106 m 3 /year.Water resources are also available beyond the capacity of PHES energy storage technology.
site for the PV system has sufficient solar capacity.Many researchers investigated whether the selected site has sufficient solar radiation to generate solar power [39,40].The site's potential for producing biogas was calculated using primary and secondary data received from various sources [41].Based on the data acquired, the biogas potential of the research area is estimated to be 7.9287 × 106 m 3 /year.Water resources are also available beyond the capacity of PHES energy storage technology.

Proposed HRES with ES System Configuration and Description
The hybrid solar PV-biogas with SMES-PHES energy storage system that is connected to the national grid and, as shown in Figure 4, which consists of HRES for solar PV and biogas generators, HESS for SMES and PHES, a connection to the national grid, and AC loads connected in the system via their respective controlled power electronic converters, is being considered for quick response and long-term investigation.To obtain the most power out of solar PV systems, MPPT systems are used.When HRES's output is not sufficient to meet load demand, the primary storage system (PHES) adds power to the system.When HRES's output is too high, PHES takes power out of the system.While PHES has high reliability, it struggles with slow response times during transient conditions.SMES is therefore employed to make up for the PHES's slow dynamics.SMES has drawn energy from the system and supplied it.The hybrid system is connected to a variable load.A hybrid system's primary output is an electricity-generating energy source.If converters are part of the system model, a power converter system changes AC power to DC power and DC power to AC power.The converter converts the SMES's DC power output to the system's supplied AC power.The chopper's duty cycle and inductor current have also been limited.This plant model applies to both grid-connected and dis- site for the PV system has sufficient solar capacity.Many researchers investigated whether the selected site has sufficient solar radiation to generate solar power [39,40].The site's potential for producing biogas was calculated using primary and secondary data received from various sources [41].Based on the data acquired, the biogas potential of the research area is estimated to be 7.9287 × 106 m 3 /year.Water resources are also available beyond the capacity of PHES energy storage technology.

Proposed HRES with ES System Configuration and Description
The hybrid solar PV-biogas with SMES-PHES energy storage system that is connected to the national grid and, as shown in Figure 4, which consists of HRES for solar PV and biogas generators, HESS for SMES and PHES, a connection to the national grid, and AC loads connected in the system via their respective controlled power electronic converters, is being considered for quick response and long-term investigation.To obtain the most power out of solar PV systems, MPPT systems are used.When HRES's output is not sufficient to meet load demand, the primary storage system (PHES) adds power to the system.When HRES's output is too high, PHES takes power out of the system.While PHES has high reliability, it struggles with slow response times during transient conditions.SMES is therefore employed to make up for the PHES's slow dynamics.SMES has drawn energy from the system and supplied it.The hybrid system is connected to a variable load.A hybrid system's primary output is an electricity-generating energy source.If converters are part of the system model, a power converter system changes AC power to DC power and DC power to AC power.The converter converts the SMES's DC power output to the system's supplied AC power.The chopper's duty cycle and inductor current have also been limited.This plant model applies to both grid-connected and dis-

Proposed HRES with ES System Configuration and Description
The hybrid solar PV-biogas with SMES-PHES energy storage system that is connected to the national grid and, as shown in Figure 4, which consists of HRES for solar PV and biogas generators, HESS for SMES and PHES, a connection to the national grid, and AC loads connected in the system via their respective controlled power electronic converters, is being considered for quick response and long-term investigation.To obtain the most power out of solar PV systems, MPPT systems are used.When HRES's output is not sufficient to meet load demand, the primary storage system (PHES) adds power to the system.When HRES's output is too high, PHES takes power out of the system.While PHES has high reliability, it struggles with slow response times during transient conditions.SMES is therefore employed to make up for the PHES's slow dynamics.SMES has drawn energy from the system and supplied it.The hybrid system is connected to a variable load.A hybrid system's primary output is an electricity-generating energy source.If converters are part of the system model, a power converter system changes AC power to DC power and DC power to AC power.The converter converts the SMES's DC power output to the system's supplied AC power.The chopper's duty cycle and inductor current have also been limited.This plant model applies to both grid-connected and distributed power generation systems.If the power users are connected to the utility point, the hybrid system can bring in and send out power.When HRES systems have more power than they need and HESS is full, the hybrid system sends the extra power to the grid.The various technical details of the HRES component parts have been presented.

Mathematical Modeling of HRES with ES System
The HRES with ES system is illustrated schematically in Figure 4 using solar PV and biogas renewable energy sources with superconducting magnetic and pumped hydro energy systems.The mathematical modelling of energy sources integrated with energy storage system is characterized using transfer functions [42,43].Transfer functions of the first order can be used to represent the dynamics of the solar PV, biogas generator, and combinations of water pump and pump turbine analysis.The gain and time constants used to model the HRES with ES system's transfer function were adopted from different researchers [44,45].

Solar PV System Unit Modeling
The expression for the solar PV unit's output power (Figure 5) is formulated in Equation ( 1) [46].
where η PV stands for photovoltaic conversion efficiency, A stands for fixed panel area (m 2 ), ϕ stands for sunshine irradiation (kw/m 2 ), and T a stands for ambient temperature ( • C).
The η PV and ϕ primarily influential factors of the Solar PV output power.The inverter unit is the components of a PV system that generate power efficiently.
used to model the HRES with ES system's transfer function were adopted from different researchers [44,45].

Solar PV System Unit Modeling
The expression for the solar PV unit's output power (Figure 5) is formulated in Equation ( 1) [46].
where PV η stands for photovoltaic conversion efficiency, A stands for fixed panel area (m 2 ), ϕ stands for sunshine irradiation (kw/m 2 ), and Ta stands for ambient temperature (°C).The PV η and ϕ primarily influential factors of the Solar PV output power.The inverter unit is the components of a PV system that generate power efficiently.
Additionally, the PV subsystem can be depicted as Equation ( 5) with a first order transfer function and lag compensation.
( ) where TPV stands for the time constant carried by the PV system, and KPV stands for the gain constant.

Biogas Turbine Generator Unit Modeling
The primary source of biogas is locally produced, biodegradable human and animal waste, which can be creatively recycled into electrical production using a biogas turbine generator as shown in Figure 6.The linearized model of biogas turbine generator is determined via computing (3) while accounting for the actions of the biogas inlet valve, governor, combustor, and turbine [47].Additionally, the PV subsystem can be depicted as Equation ( 5) with a first order transfer function and lag compensation.
where T PV stands for the time constant carried by the PV system, and K PV stands for the gain constant.

Biogas Turbine Generator Unit Modeling
The primary source of biogas is locally produced, biodegradable human and animal waste, which can be creatively recycled into electrical production using a biogas turbine generator as shown in Figure 6.The linearized model of biogas turbine generator is determined via computing (3) while accounting for the actions of the biogas inlet valve, governor, combustor, and turbine [47].

PHES System Unit Modeling.
In pumped-storage hydroelectric systems, the use of pump and turbine components is reversible, which can have many benefits.While energy is produced by using a reversible pump turbine to run the system as a turbine during peak hours, energy can also be stored by using the pump to transfer water to the upper reservoir.Instead of generating electricity, PHES units store the energy needed during times of peak demand when supply is low.They then provide this energy when demand goes up, maintaining the HRES system's frequency stable.When hydraulic reservoirs are used, energy is produced as a result of the hydraulic converter directly converting the potential energy into kinetic energy of the water in the channel.When the water flows, the following force is generated [48]:

PHES System Unit Modeling
In pumped-storage hydroelectric systems, the use of pump and turbine components is reversible, which can have many benefits.While energy is produced by using a reversible pump turbine to run the system as a turbine during peak hours, energy can also be stored by using the pump to transfer water to the upper reservoir.Instead of generating electricity, PHES units store the energy needed during times of peak demand when supply is low.They then provide this energy when demand goes up, maintaining the HRES system's frequency stable.When hydraulic reservoirs are used, energy is produced as a result of the hydraulic converter directly converting the potential energy into kinetic energy of the water in the channel.When the water flows, the following force is generated [48]: When utilizing the transit flow method with variable water availability, the gross energy output, expressed in kWh, at a nearby hydroelectric complex can be calculated as follows: a variable flow power expressed in kW, The actuator, governor, hydraulic turbine, and transient droop compensator make up the pumped-storage hydro system model.Figure 7 depicts the pump-hydro system's transfer function.

PHES System Unit Modeling.
In pumped-storage hydroelectric systems, the use of pump and turbine components is reversible, which can have many benefits.While energy is produced by using a reversible pump turbine to run the system as a turbine during peak hours, energy can also be stored by using the pump to transfer water to the upper reservoir.Instead of generating electricity, PHES units store the energy needed during times of peak demand when supply is low.They then provide this energy when demand goes up, maintaining the HRES system's frequency stable.When hydraulic reservoirs are used, energy is produced as a result of the hydraulic converter directly converting the potential energy into kinetic energy of the water in the channel.When the water flows, the following force is generated [48]: When utilizing the transit flow method with variable water availability, the gross energy output, expressed in kWh, at a nearby hydroelectric complex can be calculated as follows: ( ) (5) a variable flow power expressed in kW, ( ) The actuator, governor, hydraulic turbine, and transient droop compensator make up the pumped-storage hydro system model.Figure 7 depicts the pump-hydro system's transfer function.During the charging mode of operation, the transfer function of PHES will be [48]:

Angular velocity of the Generator
) whereas during the discharging mode of operation, the transfer function of PHES will be: where Tgh is the actuator time constant, Th1 and Th2 are the transient droop compensator system's parameters, Tw is the water's initial flow rate through the turbine, and -P-h is the rated power for water pumping unit.During the charging mode of operation, the transfer function of PHES will be [48]: whereas during the discharging mode of operation, the transfer function of PHES will be: where T gh is the actuator time constant, T h1 and T h2 are the transient droop compensator system's parameters, T w is the water's initial flow rate through the turbine, and -P-h is the rated power for water pumping unit.

Superconducting Magnetic Energy Storage System Modeling
The magnetic field created by the current flowing through a superconductor is used to store energy in the coil.This turns into a superconductor with no resistive losses when used at cryogenic temperatures between 20 K and 77 K.More than 98% of the time is spent efficiently.The superconducting coil, transformer, and power-training unit are crucial components of the energy storage system that make it up.Other advantages of SMES units include their enormous power, high energy density, quick response time, and low maintenance costs.The superconducting coil charges when the power system is operating normally because it is powered by the AC system.The current starts to conduct with fewer electric losses at cryogenic temperatures.SMES responds more slowly than the governor control (supplemental control) and vice versa.With the aid of a converter, the Sustainability 2023, 15, 5739 9 of 26 alpha communication angle is used to control the charge and discharge of the SMES.The SMES output energy can be calculated as in Equation ( 9) [49].
The SMES transfer function model is shown in Figure 8. Changes in the coil's current (∆I D ) are interpreted as negative feedback for SMES control, which increases the coil's ability to restore its current and enable quick responses to load disturbances in the future.The equations below in (10) and (11) are displayed for changes in DC voltage and current, respectively.
the alpha communication angle is used to control the charge and discharge of the SMES.
The SMES output energy can be calculated as in Equation ( 9) [49].
The SMES transfer function model is shown in Figure 8. Changes in the coil's current (∆ID) are interpreted as negative feedback for SMES control, which increases the coil's ability to restore its current and enable quick responses to load disturbances in the future.The equations below in (10) and (11) are displayed for changes in DC voltage and current, respectively.
( ) Finally, the SMES transfer function is expressed as in Equation ( 12) [21]: The following shows how the SMES unit's actual power output is represented in (13): where L stands for coil inductance (H), ∆ED for converter voltage deviation (kV), ∆ID for coil current variation (kA), KID for ∆ID feedback gain constant (kV/kA), and KF for gain constant.

Proposed HRES System Dynamics Modeling
The amount of power generated and consumed determines how the system's power and frequency are changed.Each unit's transfer functions displays the amount of power Finally, the SMES transfer function is expressed as in Equation ( 12) [21]: The following shows how the SMES unit's actual power output is represented in (13): where L stands for coil inductance (H), ∆E D for converter voltage deviation (kV), ∆I D for coil current variation (kA), K ID for ∆I D feedback gain constant (kV/kA), and K F for gain constant.

Proposed HRES System Dynamics Modeling
The amount of power generated and consumed determines how the system's power and frequency are changed.Each unit's transfer functions displays the amount of power generated by each generational subsystem.In Figure 9, the total output power (P Total ) of the HRES with ES system could be written as in (14).
Equation ( 15) is created when the dynamics of the HRES systems are represented in a configuration model [34].
In this work, the values of the inertia constant (M) and the damping constant (D) are used; M is 0.4 and D is 0.03.a configuration model [34].
( ) In this work, the values of the inertia constant (M) and the damping constant (D) are used; M is 0.4 and D is 0.03.

Proposed FO-Fuzzy-PID Controllers
For any traditional fuzzy-PID controller, the order is only taken into account as an integer value in the order of the rate of change of the error representation and the order representation of the integral at the output of the fuzzy logic controller.However, the FO-fuzzy-PID controller no longer uses integers for the rates of change of error at the input or the integral at the output.Instead, µ and λ, their equivalent tunable variables in FO, are used.A fuzzy logic controller's membership function structure has less of an impact on the effectiveness of the closed loop than Input-output variables do, which leads to the application of a scaling factor-based fuzzy logic controller.For this study, the OWOA optimization technique is used to tune the FO-fuzzy-PID controller parameters.
Das et al. introduced the FO fuzzy PID controller in [50], with [Ke, Kd] as its input scaling factors (SFs) and [KPI, KPD] as its output scaling factors (SFs).For applications in process control, it has been shown to produce positive results [32,[51][52][53].In the order of rate of change of error representation and order representation of the integral at the FLC output (i.e., µ = 1 and λ = 1 in Figure 10) for any classical fuzzy-PID controller of the FLC input level, the order is only taken into account as an integer value.

Proposed FO-Fuzzy-PID Controllers
For any traditional fuzzy-PID controller, the order is only taken into account as an integer value in the order of the rate of change of the error representation and the order representation of the integral at the output of the fuzzy logic controller.However, the FO-fuzzy-PID controller no longer uses integers for the rates of change of error at the input or the integral at the output.Instead, µ and λ, their equivalent tunable variables in FO, are used.A fuzzy logic controller's membership function structure has less of an impact on the effectiveness of the closed loop than Input-output variables do, which leads to the application of a scaling factor-based fuzzy logic controller.For this study, the OWOA optimization technique is used to tune the FO-fuzzy-PID controller parameters.
Das et al. introduced the FO fuzzy PID controller in [50], with [K e , K d ] as its input scaling factors (SFs) and [K PI , K PD ] as its output scaling factors (SFs).For applications in process control, it has been shown to produce positive results [32,[51][52][53].In the order of rate of change of error representation and order representation of the integral at the FLC output (i.e., µ = 1 and λ = 1 in Figure 10) for any classical fuzzy-PID controller of the FLC input level, the order is only taken into account as an integer value.In addition, the rule base and associated membership functions for the fuzzy controller are also shown in Table 1 and Figure 11.Fuzzy linguistic variables with the letters NL, NM, NS, ZR, PS, PM, and PL denote, respectively, negative large, negative medium, negative small, zero, positive small, positive medium, and positive large.The clear output of the fuzzy logic controller is found using the center-of-gravity method of defuzzification.The FO fuzzy controller SFs and integro-differential orders [Ke, Kd, KPI, KPD, µ, λ] are modified for a fixed rule base and membership function type using the opposition-based whale optimization algorithm.However, the FO-Fuzzy-PID controller's FLC no longer uses integers to determine the order of the "rates of change of error" at the input and the integral at the output.Instead, they are replaced by the corresponding tunable variables, namely µ and λ, of their FO counterparts.In addition, the rule base and associated membership functions for the fuzzy controller are also shown in Table 1 and Figure 11.Fuzzy linguistic variables with the letters NL, NM, NS, ZR, PS, PM, and PL denote, respectively, negative large, negative medium, negative small, zero, positive small, positive medium, and positive large.The clear output of the fuzzy logic controller is found using the center-of-gravity method of defuzzification.The FO fuzzy controller SFs and integro-differential orders [K e , K d , K PI , K PD , µ, λ] are modified for a fixed rule base and membership function type using the opposition-based whale optimization algorithm.However, the FO-Fuzzy-PID controller's FLC no longer uses integers to determine the order of the "rates of change of error" at the input and the integral at the output.Instead, they are replaced by the corresponding tunable variables, namely µ and λ, of their FO counterparts.put of the fuzzy logic controller is found using the center-of-gravity method of defuzzi-fication.The FO fuzzy controller SFs and integro-differential orders [Ke, Kd, KPI, KPD, µ, λ] are modified for a fixed rule base and membership function type using the opposition-based whale optimization algorithm.However, the FO-Fuzzy-PID controller's FLC no longer uses integers to determine the order of the "rates of change of error" at the input and the integral at the output.Instead, they are replaced by the corresponding tunable variables, namely µ and λ, of their FO counterparts.One of several continuous and discrete-time rational approximation methods for fractional order elements can be seen as the heart of the proposed fuzzy logic-based fractional order controller [54][55][56].The Oustaloup 5th order rational approximation is used in the current paper to continuously justify each guess value of the fractional order differ-integrals [µ, λ] within the optimization process.In essence, infinite-dimensional One of several continuous and discrete-time rational approximation methods for fractional order elements can be seen as the heart of the proposed fuzzy logic-based fractional order controller [54][55][56].The Oustaloup 5th order rational approximation is used in the current paper to continuously justify each guess value of the fractional order differ-integrals [µ, λ] within the optimization process.In essence, infinite-dimensional linear filters are what the FO differ-integrators are.However, for practical implementation, band-limited realizations of fractional order controllers are required.In this study, each fractional order element was interpreted using the Oustaloup recursive filter, which is illustrated in Equations ( 16) and ( 17).If we assume that the frequency range or expected fitting range of the controller is [ω b ,ω h ], then the higher order filter that most closely resembles the FO element S α can be expressed as follows [57]: The filter's poles, zeros, and gain can be assessed as follows: In Equations ( 16) and ( 17), α is the order of the differencing integration which ranges from 0 < α < 1 and |2N + 1| is the order of the realized analogue filter.This work considers a 5th-order Oustaloup approximation for all FO elements in the frequency range ω ∈ 10 −2 , 10 2 rad/s.
In the context of this work, it is assumed that an FLC's rule base and MFs are not both optimized for FO-fuzzy-PID controllers and conventional fuzzy-PID [58][59][60].To achieve the best enactment, the I/O SFs (K a , K b , K PI , and K PD ) are optimized in the fractional rate of error (i.e., µ) and FO integration (i.e., λ) of FLC output (FO-Fuzzy-PID).Table 1 provides the FLC's rule base, and Figure 11 presents fixed-shaped MFs.

Problem Formulations
The primary goal of the current research is to guarantee that the response profiles of the examined HRES and ES system models have a minimum frequency deviation and power deviation.

Objective Function
To set the controller parameters, the objective function is chosen.A decrease in overshoot, undershoot, and settling time is the system's desired dynamic response.The LFC control techniques that have been used state that the feedback system must meet the stability requirement and that the integral frequency error must be less than the specified low threshold value or close to zero.The main performance criteria for the present work is the integral of absolute error (IAE) considered.To design a controller, the objective function is primarily specified based on this parameter.Previous research has demonstrated that, of the three, the ISE objective function offers the best system performance.In this study, ISE was chosen as the objective function for the single area power system.It is written as [1,61,62]: where ∆ f represents a deviation in frequency.

Constraints
The optimization problem in the employed HRES with ES system model has constraints for optimizing a number of tunable variables.When minimizing the objective function for this study, the constraints of the objective function are given in Equation (19).
The minimum and maximum values of the parameters of the PID/FO-PID/fuzzy-PID/FO-fuzzy-PID controller in HRES with ES systems are denoted by the min and max, respectively.The optimized K P , K I , K D , K PI , and K PD parameters' boundaries are chosen from a range of 0 to 100, while the µ and λ are chosen from a range of 0 to 1.The algorithm used for this particular work is the opposition-based whale optimization algorithm (OWOA), which is introduced by Alamri et al. [63], who provided its pseudo-code.When the optimization process was complete, the program could output an ISE value.To evaluate the performance of different controllers, additional performance metrics are examined and compared in order to identify the most efficient controller [64].Finally, quasi-oppositional harmony search algorithm (QOHSA) [9,[65][66][67], teaching learning based optimization algorithm (TLBOA) [68] and particle swarm optimization (PSO) algorithm [13,33,36,69] are used for comparison with the OWOA results.

Result and Discussion
The investigation of the HRES with ES system model considers the outcomes of time-domain executions under various, distinct perturbation conditions.Under normal operating circumstances, it is assumed that the total load demand will be 1 p.u. Unexpected drops or increases in P L and PV are correctly used at specified time instants for the system being used to investigate the effects of a sudden change in demanded load or solar PV irradiation (ϕ) on the HRES with ES system model's closed-loop response.Figure 1 with a controlling action is referred to as FO-fuzzy-PID for simulation purposes in this work.The effectiveness of the opposition-based whale optimization algorithm in maximizing various controller (i.e., PID, FO-PID and fuzzy-PID) parameters is examined in this paper using different situations.The following scenarios are the input perturbation conditions:

Scenario 1: Disturbances caused by load fluctuation
In order to assess the effectiveness of the optimized values produced by the suggested opposition-based whale optimization algorithm, the studied power system model is operated under load variation ranging from 0.15 p.u. to 0.56 p.u., as shown in Figure 12.In all optimization techniques, the proposed FO-fuzzy-PID controller achieves minimum frequency and power deviations.Figure 15 shows how the frequency and power deviation parameters of the FO-fuzzy-PID controller are evaluated using OWOA, QOHSA, TLBOA, and PSO.As shown in Figure 16, OWOA achieves the best results compared to the other methods.Figures 13 and 14 show, separately, the patterns of ΔF (Hz) and ΔP (p.u.) with load demand for PID, FO-PID, fuzzy-PID, and FO-fuzzy-PID of the investigated HRES with ES system model by utilizing different metaheuristic optimization techniques such as OWOA, QOHSA, TLBOA and PSO.Figures 13 and 14 show that the proposed FO-fuzzy-PID model configuration outperforms the three other model configurations (i.e., PID, FO-PID and fuzzy-PID) in terms of performance criteria and has less maximum overshoot.The output power with the best system implementation is achieved with the suggested FO-fuzzy-PID controller using OWOA.As shown in Figures 13 and 14, the OWOA technique is better than other optimization methods in solving the optimization problem.In all optimization techniques, the proposed FO-fuzzy-PID controller achieves minimum frequency and power deviations.Figure 15 shows how the frequency and power deviation parameters of the FO-fuzzy-PID controller are evaluated using OWOA, QOHSA, TLBOA, and PSO.As shown in Figure 16, OWOA achieves the best results compared to the other methods.According to Figure 12, the OWOA results are far better than those of the other metaheuristic optimization techniques.As a result, OWOA will be used for additional research throughout the course of this work.As shown in Figures 16 and 17, power deviation is represented in the nonlinear response of SMES, biogas and PHES system utilizing OWOA.In all aspects, the FO-fuzzy-PID controller achieves the optimal solution for this study.According to Figure 12, the OWOA results are far better than those of the other metaheuristic optimization techniques.As a result, OWOA will be used for additional research throughout the course of this work.As shown in Figures 16 and 17, power deviation is represented in the nonlinear response of SMES, biogas and PHES system utilizing OWOA.In all aspects, the FO-fuzzy-PID controller achieves the optimal solution for this study.According to Figure 12, the OWOA results are far better than those of the other metaheuristic optimization techniques.As a result, OWOA will be used for additional research throughout the course of this work.As shown in Figures 16 and 17, power deviation is represented in the nonlinear response of SMES, biogas and PHES system utilizing OWOA.In all aspects, the FO-fuzzy-PID controller achieves the optimal solution for this study.It is important to keep in mind that there are specific HRES with ES system components where it can be determined that there is a significant gap between the output of the PID, FO-PID, fuzzy-PID and FO-fuzzy-PID controllers, with the FO-fuzzy-PID controller case showing the least amount of power differences.As a result, the analyzed HRES with ES system model's FO-fuzzy-PID is set up to be stronger than the other three model configurations.The numerical approximations of transient parameters for the well-known perturbation are listed in Table 2 alongside the other controllers.Furthermore, Table 3 shows the controller configurations that work well in this scenario for integrating different HRES and ES system configurations.

Scenario 2: Disturbance is limited to HRES with ES system
To evaluate the effectiveness of the optimized gains produced by the opposition-based whale optimization algorithm, the studied power system model was subjected It is important to keep in mind that there are specific HRES with ES system components where it can be determined that there is a significant gap between the output of the PID, FO-PID, fuzzy-PID and FO-fuzzy-PID controllers, with the FO-fuzzy-PID controller case showing the least amount of power differences.As a result, the analyzed HRES with ES system model's FO-fuzzy-PID is set up to be stronger than the other three model configurations.The numerical approximations of transient parameters for the well-known perturbation are listed in Table 2 alongside the other controllers.Furthermore, Table 3 shows the controller configurations that work well in this scenario for integrating different HRES and ES system configurations.

Scenario 2: Disturbance is limited to HRES with ES system
To evaluate the effectiveness of the optimized gains produced by the oppositionbased whale optimization algorithm, the studied power system model was subjected to a perturbation of solar PV generation from 0.12 to 0.55 p.u. without load perturbation.Figures 18 and 19 illustrate the characteristics of ∆F (Hz) and ∆P (p.u.) considering the perturbation of solar PV generation (shown in Figure 12) for the various configurations of HRES with ES system model (such as PID, FO-PID, fuzzy-PID, and FO-F-PID).The recommended FO-fuzzy-PID controller achieves the best HRES with ES system performance, as shown in Figures 18 and 19; in terms of steady-state error, overall overshoot (i.e., less maximum overshoot), and settling time, the proposed FO-fuzzy-PID configuration outperforms the PID, FO-PID and fuzzy-PID configurations.The power output of SMES, biogas and the PHES system is shown in Figure 20 and Figure 21, where it is possible to observe the significant differences in performance between the PID, FO-PID, fuzzy-PID and FO-fuzzy-PID controllers.The recommended FO-fuzzy-PID controller achieves the best HRES with ES system performance, as shown in Figures 18 and 19; in terms of steady-state error, overall overshoot (i.e., less maximum overshoot), and settling time, the proposed FO-fuzzy-PID configuration outperforms the PID, FO-PID and fuzzy-PID configurations.The power output of SMES, biogas and the PHES system is shown in Figure 20 and Figure 21, where it is possible to observe the significant differences in performance between the PID, FO-PID, fuzzy-PID and FO-fuzzy-PID controllers.The recommended FO-fuzzy-PID controller achieves the best HRES with ES system performance, as shown in Figures 18 and 19; in terms of steady-state error, overall overshoot (i.e., less maximum overshoot), and settling time, the proposed FO-fuzzy-PID configuration outperforms the PID, FO-PID and fuzzy-PID configurations.The power output of SMES, biogas and the PHES system is shown in Figures 20 and 21, where it is possible to observe the significant differences in performance between the PID, FO-PID, fuzzy-PID and FOfuzzy-PID controllers.
There are fewer power variations with the FO-fuzzy-PID controller.As a result, it is possible to conclude that FO-fuzzy-PID controllers with OWOA outperform PID, FO-PID and fuzzy-PID controllers with OWOA-tuned controllers.In Table 4, the numerical estimates of the transient controller parameters demonstrate how the FO-fuzzy-PID controller outperforms the fuzzy-PID, FO-PID and PID controllers in terms of performance.There are fewer power variations with the FO-fuzzy-PID controller.As a result, it is possible to conclude that FO-fuzzy-PID controllers with OWOA outperform PID, FO-PID and fuzzy-PID controllers with OWOA-tuned controllers.In Table 4, the numerical estimates of the transient controller parameters demonstrate how the FO-fuzzy-PID controller outperforms the fuzzy-PID, FO-PID and PID controllers in terms of performance.There are fewer power variations with the FO-fuzzy-PID controller.As a result, it is possible to conclude that FO-fuzzy-PID controllers with OWOA outperform PID, FO-PID and fuzzy-PID controllers with OWOA-tuned controllers.In Table 4, the numerical estimates of the transient controller parameters demonstrate how the FO-fuzzy-PID controller outperforms the fuzzy-PID, FO-PID and PID controllers in terms of performance.

Scenario 3: Disturbance from both the load and HRES with ES system
In order to evaluate the effectiveness of the optimized gains produced by the considered algorithm (i.e., OWOA), the undertaken HRES with ES system model was subjected to load perturbation and solar PV generation variation from 0.1 to 0.5 p.u. and 0.15 to 0.3 p.u., respectively, as shown in Figure 12.Figures 22 and 23 show the shapes for ∆F (Hz) and ∆P (p.u.) due to perturbation from the load and solar PV generation (as shown in Figure 12) for the different configurations of the controllers undertaken HRES with ES system, namely PID, FO-PID, fuzzy-PID and FO-fuzzy-PID.The suggested FO-fuzzy-PID controller enhances system performance, as shown in Figures 22 and 23.In terms of steady-state error, maximum overshoot (i.e., smaller maximum overshoot), and settling time, the suggested FO-fuzzy-PID configuration outperforms the PID, FO-PID and fuzzy-PID configurations.6.Table 7 shows the best value for the controller under consideration for the different ways HRES systems can be set up.    6.Table 7 shows the best value for the controller under consideration for the different ways HRES systems can be set up.

Scenario 4: Sensitivity analysis of the integrated power system
In this scenario, the proposed controller's sensitivity is analyzed against variations in the power system's block parameters of HRES with ES system from their nominal values, which is tabulated in Table 8.The suggested frequency and power deviation controller works by comparing the percentage of difference caused by changes to the parameters of the power system block (i.e., vectors M and D) with the value of the objective function.

HRES Units
The ∆F and ∆P profiles that are produced in relation to one of the agreed-upon constants are simultaneously subject to change, while the other constant is kept at its nominal value in Figure 26, and the performance measurements are noted.As a general rule in this section, when the values of the parameters of the power system block (i.e., vectors M and D) are changed significantly, there cannot be a change in significant values on frequency and power deviations underutilization of FO-fuzzy-PID controller tuned by OWOA.
The F Δ and P Δ profiles that are produced in relation to one of the agreed-upon constants are simultaneously subject to change, while the other constant is kept at its nominal value in Figure 26, and the performance measurements are noted.As a general rule in this section, when the values of the parameters of the power system block (i.e., vectors M and D) are changed significantly, there cannot be a change in significant values on frequency and power deviations underutilization of FO-fuzzy-PID controller tuned by OWOA.As a result, even with changing parameters, the suggested controller, optimized with the OWOA, convincingly demonstrates its robustness for controlling frequency and power in the HRES with ES system model.The results in Table 8 show that the FO-fuzzy-PID controller is robust against significant parameter fluctuations.
As shown in Table 9, the values of ISE are minimum, which means the proposed research findings have been achieved since the goal of this research was to minimize the numerical values of the objective function (ISE) under consideration of constraints.As a result, even with changing parameters, the suggested controller, optimized with the OWOA, convincingly demonstrates its robustness for controlling frequency and power in the HRES with ES system model.The results in Table 8 show that the FO-fuzzy-PID controller is robust against significant parameter fluctuations.
As shown in Table 9, the values of ISE are minimum, which means the proposed research findings have been achieved since the goal of this research was to minimize the numerical values of the objective function (ISE) under consideration of constraints.

Conclusions
This study proposes a HRES (i.e., solar PV and biogas generator) with an ES (superconducting magnetic and pumped hydro energy storage) system modelling and control system by using a recent controller as part of a centralized control scheme to minimize frequency and power oscillation in a hybrid power system.Classical controllers may not be as effective in addressing load frequency control issues, especially when nonlinearities and random load perturbations are present, along with electricity generated by solar PV production systems.A FO-fuzzy-PID controller is suggested for the current study for hybrid energy supply frequency control such as solar, biogas generators, and energy storage systems, such as superconducting magnetic energy storage or pumped hydro energy storage, in order to enhance the system's performance.In order to simplify the control system, save on maintenance costs, and reduce the number of control variables, a centralized control scheme has been implemented.The key conclusions of the study are listed as follows:

•
When the load and solar irradiation disrupt the system, the PID, FO-PID and fuzzy-PID controllers' configuration performs remarkably well to suppress the frequency and power oscillations of the system.

•
PID/FO-PID/fuzzy-PID/FO-fuzzy-PID controllers tuned with different metaheuristic optimization techniques (OWOA, QOHS, TLBOA, and PSO algorithms) have been proposed and implemented on HRES with ESS; as a result, OWOA is given priority and used throughout this work.
• Compared to the other system configurations that were taken into consideration, the FO-fuzzy-PID controller combination performs incredibly well and successfully reduces the oscillation of the system.

•
The FO-fuzzy-PID control technique outperforms the competition and is stable and robust in the face of random load perturbations, altered system parameters, and system nonlinearities.
When QOHSA, TBLOA, and PSO have been used to tune, respectively, the results of the frequency deviation obtained using FO-fuzzy-PID controllers with OWOA tuning are, respectively, 1.05%, 2.01%, and 2.73% lower.The proposed HRES with ES system model uses the OWOA technique, which has a highly effective computational performance, to enhance the FO-fuzzy-PID controller parameter for control of frequency and power.

Figure 1 .
Figure 1.Map of the study area.

Figure 1 .
Figure 1.Map of the study area.

Figure 2 .
Figure 2. Community's hourly load profile on a typical worst day.

Figure 2 .
Figure 2. Community's hourly load profile on a typical worst day.

Figure 3 .
Figure 3. Average daily temperature and solar radiation on the study area.

Figure 3 .
Figure 3. Average daily temperature and solar radiation on the study area.

Figure 3 .
Figure 3. Average daily temperature and solar radiation on the study area.

Figure 5 .
Figure 5.The Solar PV unit transfer function representation.

Figure 5 .
Figure 5.The Solar PV unit transfer function representation.

Figure 6 .
Figure 6.The biogas generator unit transfer function representation.

Figure 6 .
Figure 6.The biogas generator unit transfer function representation.

Figure 7 .
Figure 7.The pump-hydro unit transfer function representation.

Figure 7 .
Figure 7.The pump-hydro unit transfer function representation.

Figure 8 .
Figure 8.The SMES unit transfer function and subsystem block representation.

Figure 8 .
Figure 8.The SMES unit transfer function and subsystem block representation.

Figure 9 .
Figure 9. Schematic diagram of the HRES proposed system.

Figure 11 .
Figure 11.Schematic representation of membership function for error, FLC output, and fractional derivative of error.

• Scenario 1 :
Disturbances caused by load fluctuation.• Scenario 2: Disturbance is limited to HRES with ES system.• Scenario 3: Disturbance from both the load and HRES with ES system.• Scenario 4: Sensitivity analysis of the integrated power system.

Figures 13 and 14
Figures 13 and 14 show, separately, the patterns of ΔF (Hz) and ΔP (p.u.) with load demand for PID, FO-PID, fuzzy-PID, and FO-fuzzy-PID of the investigated HRES with ES system model by utilizing different metaheuristic optimization techniques such as OWOA, QOHSA, TLBOA and PSO.Figures 13 and 14 show that the proposed FO-fuzzy-PID model configuration outperforms the three other model configurations (i.e., PID, FO-PID and fuzzy-PID) in terms of performance criteria and has less maximum overshoot.The output power with the best system implementation is achieved with the suggested FO-fuzzy-PID controller using OWOA.As shown in Figures 13 and 14, the OWOA technique is better than other optimization methods in solving the optimization problem.

Figure 12 .
Figure 12.Graphical representation of connected load and PV generation step variation.

Figures 13 and 14
Figures 13 and 14 show, separately, the patterns of ∆F (Hz) and ∆P (p.u.) with load demand for PID, FO-PID, fuzzy-PID, and FO-fuzzy-PID of the investigated HRES with ES system model by utilizing different metaheuristic optimization techniques such as OWOA, QOHSA, TLBOA and PSO.Figures 13 and 14 show that the proposed FO-fuzzy-PID model configuration outperforms the three other model configurations (i.e., PID, FO-PID and fuzzy-PID) in terms of performance criteria and has less maximum overshoot.The output power with the best system implementation is achieved with the suggested FO-fuzzy-PID controller using OWOA.As shown in Figures 13 and 14, the OWOA technique is better than other optimization methods in solving the optimization problem.In all optimization techniques, the proposed FO-fuzzy-PID controller achieves minimum frequency and power deviations.Figure15shows how the frequency and power

Figure 12 .
Figure 12.Graphical representation of connected load and PV generation step variation.

Figure 13 .
Figure 13.Frequency deviation in the nonlinear response of HRES (a) with stepwise load changes (b) constant solar power production using OWOA, QOHSA, TLBOA and PSO tuned controllers.

Figure 14 .
Figure 14.Power deviation in the nonlinear response of HRES with (a) stepwise load changes (b) constant solar power production using OWOA, QOHSA, TLBOA and PSO tuned controllers.

Figure 16 .
Figure16.Power deviation in the nonlinear response of SMES system under constant solar output power and variable load conditions using OWOA-tuned controllers.

Figure 16 .
Figure 16.Power deviation in the nonlinear response of SMES system under constant solar output power and variable load conditions using OWOA-tuned controllers.

Figure 16 .
Figure 16.Power deviation in the nonlinear response of SMES system under constant solar output power and variable load conditions using OWOA-tuned controllers.

Figure 17 .
Figure 17.Power deviation in the nonlinear response of biogas and PHES system under constant solar output power and variable conditions using OWOA-tuned controllers.

Figure 17 .
Figure 17.Power deviation in the nonlinear response of biogas and PHES system under constant solar output power and variable conditions using OWOA-tuned controllers.
of solar PV generation from 0.12 to 0.55 p.u. without load perturbation.Figures 18 and 19 illustrate the characteristics of F Δ (Hz) and P Δ (p.u.) considering the perturbation of solar PV generation (shown in Figure 12) for the various configurations of HRES with ES system model (such as PID, FO-PID, fuzzy-PID, and FO-F-PID).

Figure 18 .
Figure 18.Frequency deviation in the nonlinear response of HRES with stepwise solar power production variations and constant connected load using OWOA-tuned controllers.

Figure 19 .
Figure 19.Power deviation in the nonlinear response of HRES with stepwise solar power production variations and constant connected load using OWOA-tuned controllers.

Figure 18 .
Figure 18.Frequency deviation in the nonlinear response of HRES with stepwise solar power production variations and constant connected load using OWOA-tuned controllers.

Figure 18 .
Figure 18.Frequency deviation in the nonlinear response of HRES with stepwise solar power production variations and constant connected load using OWOA-tuned controllers.

Figure 19 .
Figure 19.Power deviation in the nonlinear response of HRES with stepwise solar power production variations and constant connected load using OWOA-tuned controllers.

Figure 19 .
Figure 19.Power deviation in the nonlinear response of HRES with stepwise solar power production variations and constant connected load using OWOA-tuned controllers.

Figure 20 .
Figure 20.For Scenario 2, Power deviation in the nonlinear response of SMES system under constant load and variable solar output power conditions using OWOA-tuned controllers.

Figure 21 .
Figure 21.For Scenario 2, Power deviation in the nonlinear response of biogas and PHES system und constant load and variable solar output power conditions using OWOA-tuned controllers.

Figure 20 . 27 Figure 20 .
Figure 20.For Scenario 2, Power deviation in the nonlinear response of SMES system under constant load and variable solar output power conditions using OWOA-tuned controllers.

Figure 21 .
Figure 21.For Scenario 2, Power deviation in the nonlinear response of biogas and PHES system und constant load and variable solar output power conditions using OWOA-tuned controllers.

Figure 21 .
Figure 21.For Scenario 2, Power deviation in the nonlinear response of biogas and PHES system und constant load and variable solar output power conditions using OWOA-tuned controllers.

Figure 22 .
Figure 22.For Scenario 3, Frequency deviation in the nonlinear response of HRES with stepwise variations of solar power production and connected load using OWOA-tuned controllers.

Figure 22 .
Figure 22.For Scenario 3, Frequency deviation in the nonlinear response of HRES with stepwise variations of solar power production and connected load using OWOA-tuned controllers.Figures 24 and 25 depict the power output of the various HRES sources, demonstrating how the PID, FO-PID, and FO-fuzzy-PID controllers behave very differently, despite the FO-fuzzy-PID controller having less power fluctuation.In comparison with the OWOAtuned PID and FO-PID controllers, it is safe to say that the OWOA-tuned FO-fuzzy-PID controller provides the best system responses.The transient parameters for scenario 3 is

Figure 22 .
Figure 22.For Scenario 3, Frequency deviation in the nonlinear response of HRES with stepwise variations of solar power production and connected load using OWOA-tuned controllers.

Figure 23 .
Figure 23.For Scenario 3, Power deviation in the nonlinear response of HRES with stepwise variations of solar power production and connected load using OWOA-tuned controllers.

Figures 24
Figures 24 and 25 depict the power output of the various HRES sources, demonstrating how the PID, FO-PID, and FO-fuzzy-PID controllers behave very differently, despite the FO-fuzzy-PID controller having less power fluctuation.In comparison with the OWOA-tuned PID and FO-PID controllers, it is safe to say that the OWOA-tuned FO-fuzzy-PID controller provides the best system responses.The transient parameters for scenario 3 is contrasted with the parameters for other controllers in Table6.Table7shows the best value for the controller under consideration for the different ways HRES systems can be set up.

Figure 23 . 27 Figure 24 .
Figure 23.For Scenario 3, Power deviation in the nonlinear response of HRES with stepwise variations of solar power production and connected load using OWOA-tuned controllers.Sustainability 2023, 15, 5739 22 of 27

Figure 25 .
Figure 25.For Scenario 3, Power deviation in the nonlinear response of biogas and PHES system under constant load and variable solar output power conditions using OWOA-tuned controllers.

Figure 24 . 27 Figure 24 .
Figure 24.For Scenario 3, Power deviation in the nonlinear response of SMES system under constant load and variable solar output power conditions using OWOA-tuned controllers.

Figure 25 .
Figure 25.For Scenario 3, Power deviation in the nonlinear response of biogas and PHES system under constant load and variable solar output power conditions using OWOA-tuned controllers.

Figure 25 .
Figure 25.For Scenario 3, Power deviation in the nonlinear response of biogas and PHES system under constant load and variable solar output power conditions using OWOA-tuned controllers.

Figure 26 .
Figure 26.Sensitivity investigation of the HPS with the FO-fuzzy-PID controller for ± 10% variations in D&M.

Figure 26 .
Figure 26.Sensitivity investigation of the HPS with the FO-fuzzy-PID controller for ±10% variations in D&M.
Total = P PV + P B + P Dis PHES + P Dis SMES − P L + P Cha PHES + P Cha SMES

Table 1 .
Rule base for error, FLC output and fractional rate of error.

Table 1 .
Rule base for error, FLC output and fractional rate of error.

Table 1 .
Rule base for error, FLC output and fractional rate of error.
Figure 11.Schematic representation of membership function for error, FLC output, and fractional derivative of error.

Table 2 .
Transient response parameters of the HRES with controllers for scenario 1.

Table 3 .
Optimized values of the controller under consideration for various HRES with ES system configurations for scenario 1.

Table 2 .
Transient response parameters of the HRES with controllers for scenario 1.

Table 3 .
Optimized values of the controller under consideration for various HRES with ES system configurations for scenario 1.

Table 4 .
Transient response parameters of the HRES with controllers for scenario 2.

Table 5
shows the best values of the controllers for the different HRES with ES system configurations, as well as numerical estimates of the transient controller parameters for this scenario.

Table 5 .
Optimized values of the controller under consideration for various HRES with ES system configurations for Scenario 2.

Table 6 .
Transient response parameters of the HRES with controllers for scenario 3.

Table 6 .
Transient response parameters of the HRES with controllers for scenario 3.

Table 6 .
Transient response parameters of the HRES with controllers for scenario 3.

Table 7 .
Optimized values of the controller under consideration for various HRES with ES system configurations for Scenario 3.

Table 8 .
Nominal values of HRES components.

Table 9 .
Analysis of sensitivity to variations in system parameters.

Table 9 .
Analysis of sensitivity to variations in system parameters.