Tidal Supplementary Control Schemes-Based Load Frequency Regulation of a Fully Sustainable Marine Microgrid

The world is targeting fully renewable power generation by the middle of the century. Distributed generation is the way to increase the penetration level of renewable energies. This paper presents load frequency control of a hybrid tidal, wind, and wave microgrid to feed an isolated island. This research is a step towards 100% renewable energy communities in remote seas/oceans islands. The wave and tidal generation systems model are presented. The study presents load frequency control through three supplementary control strategies: conventional integrators, fractional order integrator, and non-linear fractional order integrator. All the controllers of the microgrid are designed by using a novel black widow optimization technique. The applied technique is compared to other existing state-of-the-art algorithms. The results show that the black widow non-linear fractional integrator has a better performance over other strategies. Coordination between the unloaded tidal system and blade pitch control of both wind and tidal systems are adopted in the microgrid to utilize the available reserve power for the frequency support. Simulation and optimization studies are performed using the MATLAB/SIMULINK 2017a software application.


Introduction
The world is targeting 100% renewable power generation by the middle of this century [1]. Distributed renewable energy generation worldwide is increasing due to its low carbon dioxide emission and cost. The integration of renewable energies to the grid or high penetration level of renewables standalone systems suffers in terms of stability, power quality, and reliability [2]. Load frequency control plays an important role in terms of improving the power quality of the power systems, which operates with high variability in both loads and generation due to the presence of renewables [3].
Tidal energy is considered one of the promising renewable energy technologies in the 100% renewable energies dream. In [4], the researchers presented a feasibility study of tidal energy applications. Tidal power plants were always associated with offshore wind turbines to cover consumption needs. The integration of wind and tidal power plants to a grid or a power system were welcomed despite the frequency stability problems. With the increase of the penetration level of renewables in standalone systems and interconnected power systems, renewable energies must participate in the load frequency control process. The contribution is made in tidal and wind generators 1 of 17 by deloading. Deloading can be defined as the operation of the wind or tidal system at a power below the maximum power to create a reserve. The reserve can be utilized to stabilize the power system for frequency regulation during active power mismatch between generation and demand [5]. However, using deloading only cannot drive to the observed reduction in frequency deviation, but involving inertia and damping controls can reduce the reduction [6].
There are a lot of technologies to convert tidal energy into electrical energy. The best known technologies are barrages, turbines, and fences. In this paper, turbines technology is assumed to be used to convert tidal energy into a mechanical one. Through doubly-fed induction generators, the mechanical energy is converted to an electric one. In [7], the researchers presented a complete study on the dynamic behavior of pitch and stall regulated tidal turbines. In [8], the effect of integration tidal power plant on a real grid is discussed. Maximum power point tracking under different tides speed is presented in [9]. In [10,11], hybrid offshore wind tidal generating system is presented. Based on the literature review, there is not much focus on load frequency control by using the deloading of tidal generators. In this work, the deloading is applied by using inertia, damping, and various supplementary control schemes.
Wave energy has good potential in the 100% renewable power generation goal. It can satisfy more than 10% of the total global demand for electricity [12]. The structure of the Archimedes wave swing is presented in [13]. In this paper, wave energy was simulated as an uncontrollable generation system.
Control of microgrids is an important topic addressed by many researchers and applied using different optimization techniques. In [14], a particle swarm algorithm (PSO) is presented. In [15], a genetic algorithm (GA) is presented. In [16], the teaching and learning algorithm (TLBO) is presented, and in [17], harmony search is presented. The black widow optimization algorithm is presented for the first time in [18] by Hayyolalam and Kazem as a novel meta-heuristic approach for solving engineering optimization problems.
Load frequency control was previously applied by using different control schemes: proportional integral (PI) control [19], proportional integral derivative (PID) control [20], H-infinity [21], fractional order proportional integral derivative (FOPID) control [22], and non-linear proportional integral derivative (NPID) control [23]. The researchers proved the simplicity of PID, while better performance can be achieved through FOPID and NPID. In this paper, the supplementary control is designed to be a non-linear fractional integrator. At the same time, other controllers will be set to be PID to ensure both simplicity and quick response to a frequency deviation.
In [24], the research presents load frequency control of standalone tidal and diesel microgrid. The research presents a contribution through supplementary control from tidal generation to demand change.
In [25], the research load frequency control for wind-diesel generation microgrid using the D-partition method (DPM). The research provided a single step /simplistic computing method for calculating the PI controller parameters of a dynamic system, such as the microgrid system comprising of the renewable energy sources without any further requirements of retuning.
In [26], a novel fractional-order model predictive control technique is presented to track the optimal frequency of a standalone microgrid through including fractional-order integral cost function into model predictive control (MPC) algorithm.
In [27], In this research, a simulation and control of tidal generation system has been presented. A fuzzy system has been used for the pitch controller to properly modify the gains of the PID at different tidal inputs.
The main contributions of this paper are: 1.
Simulation and control of a 100% renewable energy microgrid including tidal, wave and offshore wind hybrid generation system.

2.
The effect of different supplementary control schemes in terms of the integrator, fractional integrator, and non-linear fractional integrator on the dynamic performance of load frequency control (LFC) is examined. The paper is organized as follows. Section 2 presents the modelling of the microgrid, Section 3 illustrates the applied control schemes. Section 4 presents the controller's design process, Section 5 presents simulation results, and Section 6 is the main conclusion of this research.

Microgrid Modelling
The world is targeting 100% renewable energies by 2050. This paper presents a study on a standalone microgrid operated by fully sustainable marine generating systems, which are tidal, wind, and wave generators, as illustrated in Figure 1. The system's mathematical descriptions are presented in the following subsections.
Inventions 2020, 5, 53 3 of 18 3. Design of hybrid system controllers and tidal supplementary controller by using a novel black widow optimization technique and comparing it with other state-of-art optimization techniques.
The paper is organized as follows. Section 2 presents the modelling of the microgrid, Section 3 illustrates the applied control schemes. Section 4 presents the controller's design process, Section 5 presents simulation results, and Section 6 is the main conclusion of this research.

Microgrid Modelling
The world is targeting 100% renewable energies by 2050. This paper presents a study on a standalone microgrid operated by fully sustainable marine generating systems, which are tidal, wind, and wave generators, as illustrated in Figure 1. The system's mathematical descriptions are presented in the following subsections.

Modeling of the Tidal Generating System
Both wind generating systems and tidal generating systems have similar operation and control principles. The difference between the two systems is that the tidal speed and nominal turbine size are less than for the wind. The wind rated speed ranges from 12-15 m/sec and the tidal one from 2-3 m/sec. The mechanical power output ( ) can be modeled as illustrated in (1) [4].
where , , and are the seawater density, turbine blades area, and tidal speed flow, respectively; is the power coefficient in terms of tip speed ratio ( ) and blade pitch control angle ( ); while , , , , , , and are ( , ) equation parameters.
Such that and are the rotational speed of the tidal blades and radius of the blades, respectively. The tidal generating system has four modes of operation, as illustrated in Figure 2 and Table 1.

Modeling of the Tidal Generating System
Both wind generating systems and tidal generating systems have similar operation and control principles. The difference between the two systems is that the tidal speed and nominal turbine size are less than for the wind. The wind rated speed ranges from 12-15 m/s and the tidal one from 2-3 m/s. The mechanical power output (P T ) can be modeled as illustrated in (1) [4].
where P, A, and V are the seawater density, turbine blades area, and tidal speed flow, respectively; C P is the power coefficient in terms of tip speed ratio (γ) and blade pitch control angle (β); while d 1 , d 2 , d 3 , d 4 , d 5 , d 6 , and d 7 are C P (γ, β) equation parameters.
Such that ω and R are the rotational speed of the tidal blades and radius of the blades, respectively. The tidal generating system has four modes of operation, as illustrated in Figure 2 and Table 1.
Turbine power variations are governed by the pitch control system comprising of a PID controller, which gets its input as error between measured turbine rotor speed and change in reference speed. The power varies with blade pitch control to what achieves the relationship between power output and rotor speed, as illustrated in Figure 2 [24].  Turbine power variations are governed by the pitch control system comprising of a PID controller, which gets its input as error between measured turbine rotor speed and change in reference speed. The power varies with blade pitch control to what achieves the relationship between power output and rotor speed, as illustrated in Figure 2 [24].
Participation for frequency support is a must for 100% renewable generating systems. To ensure that tidal generating systems participate in frequency control, it is necessary to operate them at a level well below the maximum power point. This phenomenon is called deloading such that the power output of the system is varied between deloading power and maximum power as shown in Figure 3. This could happen by varying the rotor speed from deloading speed ( ) to nominal or rated speed ( ). The maximum deloading percentage can be calculated based on the maximum allowed rotor speed for the generator, as illustrated in Figure 4.  Participation for frequency support is a must for 100% renewable generating systems. To ensure that tidal generating systems participate in frequency control, it is necessary to operate them at a level well below the maximum power point. This phenomenon is called deloading such that the power output of the system is varied between deloading power P del and maximum power P max as shown in Figure 3. This could happen by varying the rotor speed from deloading speed (ω del ) to nominal or rated speed (ω r ). The maximum deloading percentage x can be calculated based on the maximum allowed rotor speed for the generator, as illustrated in Figure 4.
Inventions 2020, 5, 53 4 of 18  Turbine power variations are governed by the pitch control system comprising of a PID controller, which gets its input as error between measured turbine rotor speed and change in reference speed. The power varies with blade pitch control to what achieves the relationship between power output and rotor speed, as illustrated in Figure 2 [24].
Participation for frequency support is a must for 100% renewable generating systems. To ensure that tidal generating systems participate in frequency control, it is necessary to operate them at a level well below the maximum power point. This phenomenon is called deloading such that the power output of the system is varied between deloading power and maximum power as shown in Figure 3. This could happen by varying the rotor speed from deloading speed ( ) to nominal or rated speed ( ). The maximum deloading percentage can be calculated based on the maximum allowed rotor speed for the generator, as illustrated in Figure 4.   Coordination between blade pitch control system and deloading system is required such that with the increase of the deloading percentage, the rotor speed may increase. At that instant, the blade pitch control system should adjust the blade pitch angle. The dynamic power reference ( ) at a specific rotor speed and reference speed ( ) can be defined as: where the measured mechanical toque is . To improve the frequency response, an additional signal ∆ will be added to the reference power output, as shown in Figure 5. The additional signal is coming from combined inertia and damping non-conventional machine equivalent controller. In other words, the frequency deviation and rate of change of frequency are represented in two signals, as illustrated in (7).
where and are the additional damping and the additional inertia, respectively. To improve the response of the tidal system to the frequency deviation, different control schemes are compared to be added in parallel to and .  Coordination between blade pitch control system and deloading system is required such that with the increase of the deloading percentage, the rotor speed may increase. At that instant, the blade pitch control system should adjust the blade pitch angle. The dynamic power reference (P r ) at a specific rotor speed and reference speed (ω re f ) can be defined as: where the measured mechanical toque is T m .
To improve the frequency response, an additional signal ∆P id will be added to the reference power output, as shown in Figure 5. The additional signal is coming from combined inertia and damping non-conventional machine equivalent controller. In other words, the frequency deviation and rate of change of frequency are represented in two signals, as illustrated in (7).
where D 1 and M are the additional damping and the additional inertia, respectively. To improve the response of the tidal system to the frequency deviation, different control schemes are compared to be added in parallel to D 1 and M.
To drive the speed to track the reference speed, a PID control is applied, as illustrated in (8) and shown in Figure 6.
where ∆ω e is the error in speed deviation while K ωTP , K ωTI , and K ωTD are the speed controller parameters. Thus, the change of the output power of the tidal generating system (∆P Tidal ) can be formulated as (9).
where ∂P ∂ω , ∂P ∂v , and ∂P ∂β are the change of tidal power with respect to the specific variation of turbine rotor speed, tidal speed, and blade pitch angle, respectively.
Tidal generating system parameters values are illustrated in Appendix A.
where and are the additional damping and the additional inertia, respectively. To improve the response of the tidal system to the frequency deviation, different control schemes are compared to be added in parallel to and .  To drive the speed to track the reference speed, a PID control is applied, as illustrated in (8) and shown in Figure 6.
where ∆ is the error in speed deviation while , , and are the speed controller parameters. Thus, the change of the output power of the tidal generating system (∆ ) can be formulated as (9).
where , , and are the change of tidal power with respect to the specific variation of turbine rotor speed, tidal speed, and blade pitch angle, respectively.

Modeling of Wave Generating System
The wave generating system is one of the promising marine sustainable energy systems, which is not yet widely covered in the research. The wave system is assumed to be coupled with a permanent magnet synchronous generator. The following equations describe the dynamics of the system in terms of force ( ) and velocity ( ) [12,13].
where is the floater and translator displacement, is the total mass, is the damping constant of the generator, and is the damping constant of the wave swing. The system is simulated as a first-order generator, first-order converter, and first-order inverter, as illustrated in (12).
where is the transfer function of the system, and and are the gain and time constants of the wave generator. and are the time constants of converter and inverter, respectively.

Modeling of Wind Generating System
The offshore wind generating system applied in this study is assumed to participate in frequency

Modeling of Wave Generating System
The wave generating system is one of the promising marine sustainable energy systems, which is not yet widely covered in the research. The wave system is assumed to be coupled with a permanent magnet synchronous generator. The following equations describe the dynamics of the system in terms of force (F W ) and velocity (V W ) [12,13].
where x is the floater and translator displacement, m f t is the total mass, β G is the damping constant of the generator, and β W is the damping constant of the wave swing.
The system is simulated as a first-order generator, first-order converter, and first-order inverter, as illustrated in (12).
where G wave is the transfer function of the system, and K wave and T wave are the gain and time constants of the wave generator. T conv and T inv are the time constants of converter and inverter, respectively.

Modeling of Wind Generating System
The offshore wind generating system applied in this study is assumed to participate in frequency deviation by PID blade pitch control. The system is shown in Figure 7 and modeled in [22,25].

Modeling of Microgrid
The microgrid is now formed from three marine systems. The change of the total generated power (∆ ) can be formulated as (13).
The difference between the change in generation and the change in load is ∆ The transfer function of the power system ( ) in terms of changes in system frequency (∆ ) and ∆ is illustrated in (15). represents the frequency dependency of the load while is the microgrid moment of inertia.

Blade Pitch Controllers
Wind and tidal turbines have blade pitch controllers, which are selected to be proportional integral derivative (PID) controllers. PID controllers are the most widely used tool to minimize the error. The input of the PID controller is the error, and the output of the controller is the control action received by the actuators. Appropriate power injection is usually done through PID control, which usually leads to a minimum frequency deviation: where is the proportional gain, is the integral gain, and is the derivative gain.

Tidal Speed Controller
The tidal speed controller is selected to be PID controller, like the blade pitch controller.

Tidal Supplementary Control Schemes
Different tidal supplementary control schemes are applied to minimize the frequency deviation of the microgrid, which are as follows:

Modeling of Microgrid
The microgrid is now formed from three marine systems. The change of the total generated power (∆P G ) can be formulated as (13).
The difference between the change in generation and the change in load is ∆P Gd The transfer function of the power system (G PS ) in terms of changes in system frequency (∆ f ) and ∆P Gd is illustrated in (15). D represents the frequency dependency of the load while H is the microgrid moment of inertia.

Blade Pitch Controllers
Wind and tidal turbines have blade pitch controllers, which are selected to be proportional integral derivative (PID) controllers. PID controllers are the most widely used tool to minimize the error. The input of the PID controller is the error, and the output of the controller is the control action received by the actuators. Appropriate power injection is usually done through PID control, which usually leads to a minimum frequency deviation: where K P is the proportional gain, K I is the integral gain, and K D is the derivative gain.

Tidal Speed Controller
The tidal speed controller is selected to be PID controller, like the blade pitch controller.

Tidal Supplementary Control Schemes
Different tidal supplementary control schemes are applied to minimize the frequency deviation of the microgrid, which are as follows: where λ is the fractional-order operator of the integrator while G is the non-linearity gain, E is the error, and U NFI is the output of the non-linear fractional integrator. Figure 8 shows the structure of NFI. where is the fractional-order operator of the integrator while is the non-linearity gain, E is the error, and is the output of the non-linear fractional integrator. Figure 8 shows the structure of NFI.

Optimization Problem Definition
The design of wind and tidal turbines with blade pitch controllers, tidal speed controllers, and integrator based tidal supplementary control schemes will be performed for test operation condition using the following algorithms: i. Black widow ii.
Quasi oppositional harmony search iii.
Teaching and learning-based optimization iv.
Particle swarm optimization v.
Genetic algorithm The best performance technique will be used for the rest of the subjected disturbances. The design is made in MATLAB to minimize frequency deviation. The optimization problem is described as follows: • Variables: PID control parameters, tidal additional damping, and inertia, in addition to tidal supplementary control schemes parameters. • Constraints: G and λ.
The multi-objective function is established such that weights of the ITAE and IAE are equal. The tidal supplementary control schemes parameters differ from one to another. In fractional integrators, there are two variables (KI and λ) and the constraint considered in its design is 0 ≤ λ ≤ 1 only. In nonlinear fractional integrators, there are three variables (KI, G, and λ) and the constraints considered are 0 ≤ λ ≤ 1 and 0 ≤ G ≤ 1 only.

Black Widow Optimization
Recently, due to the complexity of renewable energy-based power systems, the need for a viable meta-heuristic method has emerged. Nature-inspired optimization algorithms are used widely for

Optimization Problem Definition
The design of wind and tidal turbines with blade pitch controllers, tidal speed controllers, and integrator based tidal supplementary control schemes will be performed for test operation condition using the following algorithms: i.
Black widow ii. Quasi oppositional harmony search iii. Teaching and learning-based optimization iv. Particle swarm optimization v.
Genetic algorithm The best performance technique will be used for the rest of the subjected disturbances. The design is made in MATLAB to minimize frequency deviation. The optimization problem is described as follows: • Variables: PID control parameters, tidal additional damping, and inertia, in addition to tidal supplementary control schemes parameters. • Constraints: G and λ.
The multi-objective function is established such that weights of the ITAE and IAE are equal. The tidal supplementary control schemes parameters differ from one to another. In fractional integrators, there are two variables (K I and λ) and the constraint considered in its design is 0 ≤ λ ≤ 1 only. In non-linear fractional integrators, there are three variables (K I, G, and λ) and the constraints considered are 0 ≤ λ ≤ 1 and 0 ≤ G ≤ 1 only.

Black Widow Optimization
Recently, due to the complexity of renewable energy-based power systems, the need for a viable meta-heuristic method has emerged. Nature-inspired optimization algorithms are used widely for solving power system problems in an easy and flexible way. Their method is inspired by the unique mating behavior of the black widow spider. The proposed optimization technique includes the cannibalism stage, which omits inappropriate fitness from the selection circle, so that convergence comes earlier. Figure 9 illustrates the proposed algorithm steps, which are:

1.
Initial population: This is used in each optimization technique; it has other names, like chromosomes in the genetic algorithm and particle position in the particle swarm algorithm. In the black window, it has the name widow. To start the optimization, a candidate widow matrix with size N pop × N var is generated, where N var represents the solution of the problem array while N pop represents the number of populations.

2.
Procreate: In this step, an array called α is created such that the offspring is produced through (20) and (21): where x 1 and x 2 are parents, while y 1 and y 2 are offspring. The process repeated every Nvar 2 .

3.
Cannibalism Three parameters must be set for obtaining the desired results: procreating rate (pp), cannibalism rate, (CR) and mutation rate (pm). In this research, the parameters are selected as pp = 0.62, CR = 0.46, and pm = 0.4. Three parameters must be set for obtaining the desired results: procreating rate (pp), cannibalism rate, (CR) and mutation rate (pm). In this research, the parameters are selected as pp = 0.62, CR = 0.46, and pm = 0.4.

Simulation Results
In this part, the dynamic performance of the microgrid in terms of load frequency control is presented under different operating conditions. It is assumed that the tidal system is operating at tide speed 2.4 m/sec having a deloading effect with a 3• blade pitch control angle. The microgrid is simulated and optimization processes are applied on MATLAB/SIMULINK 2017a works in a Core i5, 2.50 GHz Samsung laptop with 6 GB of RAM. The parameters of all controllers in the wind and tidal systems are optimized using black widow optimization and its code is established using MATLAB 2017a. In order to establish the supremacy of black widow for the present work, system performance for 1% step load increase using integrator tidal supplementary control scheme is compared to other state-of-the-art optimization methods. Figures 10-12 show the frequency deviation, change in tidal power, and change in wind power. Table 2 shows a comparison between the performance of each optimization technique on ITAE, IAE, (∆ )2, transient response of ∆f in terms of undershoot (Ush), overshoot (Osh), settling time (ts) in addition to peak time (tp), and the number of iterations performed using each optimization technique. The results show that the black widow algorithm has the best performance over other algorithms. Therefore, the tests applied on the system to compare between different supplementary controllers will be applied using the black widow algorithm only. The dynamic study of the microgrid in terms of load frequency control, under the action of black widow tuned control schemes installed in the studied system, was subjected to the following tests:

Simulation Results
In this part, the dynamic performance of the microgrid in terms of load frequency control is presented under different operating conditions. It is assumed that the tidal system is operating at tide speed 2.4 m/sec having a deloading effect with a 3• blade pitch control angle. The microgrid is simulated and optimization processes are applied on MATLAB/SIMULINK 2017a works in a Core i5, 2.50 GHz Samsung laptop with 6 GB of RAM. The parameters of all controllers in the wind and tidal systems are optimized using black widow optimization and its code is established using MATLAB 2017a. In order to establish the supremacy of black widow for the present work, system performance for 1% step load increase using integrator tidal supplementary control scheme is compared to other state-of-the-art optimization methods. Figures 10-12 show the frequency deviation, change in tidal power, and change in wind power. Table 2 shows a comparison between the performance of each optimization technique on ITAE, IAE, (∆ f ) 2 , transient response of ∆f in terms of undershoot (U sh ), overshoot (O sh ), settling time (t s ) in addition to peak time (t p ), and the number of iterations performed using each optimization technique. The results show that the black widow algorithm has the best performance over other algorithms. Therefore, the tests applied on the system to compare between different supplementary controllers will be applied using the black widow algorithm only. The dynamic study of the microgrid in terms of load frequency control, under the action of black widow tuned control schemes installed in the studied system, was subjected to the following tests: Test 1: A step increase in the demand; Test 2: Real demand variation at a certain day; Test 3: Sinewave variation of the wave generation system;

System Performance under Test 1
In this test, the studied microgrid is subjected to a step increase in demand by 10% 1 second after starting the simulation. Four control strategies have been carried out on the system:  Figures 13-15 show the frequency deviation, change in tidal power, and change in wind power when the system is subjected to Test 1. The results show that the non-linear fractional integrator achieved a better performance than the fractional integrator, followed by the conventional integrator. To avoid mechanical oscillations, the optimal parameters of the supplementary control, added damping (D 1 ) and added inertia (M) in addition to the rest of the controllers of the system, are tuned at the same time to get the optimal parameters for the whole system considering the mechanical oscillations. If the integrator is not inserted in the system, there is a steady-state error and a very long settling time. System controllers and different supplementary control schemes parameters are illustrated in Table 3. Table 4 illustrates the performance of the microgrid at different frequency dependencies of the load (D); the results show that as D increases, the performance of the microgrids gets better in terms of ITAE, IAE, and settling time: in other words, as the loads used are more frequency dependent, the microgrid performance improves.  In this test, the studied microgrid is subjected to a step increase in demand by 10% 1 second after starting the simulation. Four control strategies have been carried out on the system: The results show that the non-linear fractional integrator achieved a better performance than the fractional integrator, followed by the conventional integrator. To avoid mechanical oscillations, the optimal parameters of the supplementary control, added damping ( ) and added inertia ( ) in addition to the rest of the controllers of the system, are tuned at the same time to get the optimal parameters for the whole system considering the mechanical oscillations. If the integrator is not inserted in the system, there is a steady-state error and a very long settling time. System controllers and different supplementary control schemes parameters are illustrated in Table 3. Table 4 illustrates the performance of the microgrid at different frequency dependencies of the load ( ); the results show that as increases, the performance of the microgrids gets better in terms of ITAE, IAE, and settling time: in other words, as the loads used are more frequency dependent, the microgrid performance improves. Figure 13. Change in frequency when the system was subjected to test 1. Figure 13. Change in frequency when the system was subjected to test 1.

System Performance under Test 2
In this test, the studied microgrid is subjected to real-time load variations (shown in Figure 16). The same control strategies of test 1 are applied in this test. The same parameters as test 1 are kept for all controllers and integrators. The competitive profiles of frequency deviation, in addition to the change in tidal generation and change in wind generation are presented in Figure 17, Figure 18, and Figure 19, respectively. The results show that the non-linear fractional integrator has a better performance than the other control schemes.

System Performance under Test 2
In this test, the studied microgrid is subjected to real-time load variations (shown in Figure 16). The same control strategies of test 1 are applied in this test. The same parameters as test 1 are kept for all controllers and integrators. The competitive profiles of frequency deviation, in addition to the change in tidal generation and change in wind generation are presented in Figure 17, Figure 18, and Figure 19, respectively. The results show that the non-linear fractional integrator has a better performance than the other control schemes.

System Performance under Test 3
In this test, the wave system has sinusoidal power variation as expressed in the following equation [12,13].
The same parameters as test 1 are kept for all controllers and integrators. The competitive profiles of frequency deviation, in addition to change in tidal generation and change in wind generation all measured per unit (p.u.), are presented in Figure 20, Figure 21, and Figure 22, respectively. The results show that the non-linear fractional integrator has a better performance than the other control schemes.

System Performance under Test 3
In this test, the wave system has sinusoidal power variation as expressed in the following equation [12,13].
The same parameters as test 1 are kept for all controllers and integrators. The competitive profiles of frequency deviation, in addition to change in tidal generation and change in wind generation all measured per unit (p.u.), are presented in Figure 20, Figure 21, and Figure 22, respectively. The results show that the non-linear fractional integrator has a better performance than the other control schemes.

System Performance under Test 3
In this test, the wave system has sinusoidal power variation as expressed in the following equation [12,13].
∆P wave = 0.002 sin(4.36t) + 0.005 sin(5.3t) − 0.01 sin(6t) The same parameters as test 1 are kept for all controllers and integrators. The competitive profiles of frequency deviation, in addition to change in tidal generation and change in wind generation all measured per unit (p.u.), are presented in Figure 20, Figure 21, and Figure 22, respectively. The results show that the non-linear fractional integrator has a better performance than the other control schemes.

Conclusions
This paper has presented the load frequency control of a 100% renewable energy marine microgrid in terms of wind, tidal, and waves generators. The results show that using a tidal supplementary controller in the presence of an integrator drives the microgrid to zero steady-state frequency deviation in different operating conditions. The results show that the contribution of the

Conclusions
This paper has presented the load frequency control of a 100% renewable energy marine microgrid in terms of wind, tidal, and waves generators. The results show that using a tidal supplementary controller in the presence of an integrator drives the microgrid to zero steady-state frequency deviation in different operating conditions. The results show that the contribution of the

Conclusions
This paper has presented the load frequency control of a 100% renewable energy marine microgrid in terms of wind, tidal, and waves generators. The results show that using a tidal supplementary controller in the presence of an integrator drives the microgrid to zero steady-state frequency deviation in different operating conditions. The results show that the contribution of the

Conclusions
This paper has presented the load frequency control of a 100% renewable energy marine microgrid in terms of wind, tidal, and waves generators. The results show that using a tidal supplementary controller in the presence of an integrator drives the microgrid to zero steady-state frequency deviation in different operating conditions. The results show that the contribution of the tidal supplementary controller to the load or generation variation is more effective than the wind blade pitch controller in different operating conditions. The results also proved that the proposed non-linear fractional order integrator (NFI) based supplementary control achieves better performance than fractional order integrator (FI) and conventional integrator (I) schemes in different operating conditions. The results also proved that the FI control scheme drives the system to a better performance than the I control scheme. The results also proved that the controller design using a black widow optimization algorithm drives the system to a better performance than other existing state-of-the-art algorithms, in terms of ITAE, IAE, number of iterations, and change in frequency transient response. The results also show that using an NFI control scheme will lead the system to a better performance when subjected to sinusoidal wave power generation than will FI and I control schemes. The paper presents a technique to solve the frequency deviation in sea/ocean isolated microgrid. To practically apply this technique, integration between tidal, wind and wave technologies and software and microcontroller facilities, are required.