A New Fractional-Order Virtual Inertia Support Based on Battery Energy Storage for Enhancing Microgrid Frequency Stability

: Microgrids have a low inertia constant due to the high penetration of renewable energy sources and the limited penetration of conventional generation with rotating mass. This makes microgrids more susceptible to frequency stability challenges. Virtual inertia control (VIC) is one of the most effective approaches to improving microgrid frequency stability. Therefore, this study proposes a new model to precisely mimic inertia power based on an energy storage system (ESS) that supports low-inertia power systems. The developed VIC model considers the effect of both the DC-DC converter and the DC-AC inverter on the power of the ESS used. This allows for more precise and accurate modeling of the VIC compared to conventional models. Moreover, this study proposes a fractional-order derivative control for the proposed VIC model to provide greater flexibility in dealing with different perturbations that occur in the system. Furthermore, the effectiveness of the proposed fractional-order VIC (FOVIC) is verified through an islanded microgrid that includes heterogeneous sources: a small thermal power plant, wind and solar power plants, and ESSs. The simulation results performed using MATLAB software indicate that the proposed VIC scheme provides fast stabilization times and slight deviations in system frequency compared to the conventional VIC schemes. The proposed VIC outperforms the conventional load frequency control by about 80% and the conventional VIC model by about 45% in tackling load/RESs fluctuations and system uncertainty. Additionally, the studied microgrid with the proposed FOVIC scheme is noticeably more stable and responds faster than that designed with integer-order derivative control. Thus, the proposed FOVIC scheme gives better performance for frequency stability of low-inertia power systems compared to conventional VIC schemes used in the literature.


Introduction
There is a worldwide interest in generating electricity from renewable energy sources (RESs).Various RESs are being integrated into power systems, such as hydro, solar, wind, geothermal, biomass, ocean energy, etc.Many countries aim to achieve more than 90% electricity generation from RESs by 2050 [1].Energy security concerns, environmental issues, and economic benefits drive this interest.However, the increase in intermittent RESs penetration results in tedious power system operation and control due to uncertainty and intermittency of RESs production [2].
Synchronous generators (SGs)-led conventional power systems have the inertia capability to maintain voltage and frequency deviations within standard ranges [3].Though the decreasing inertia makes future power systems driven by power converters more susceptible to system insecurity [4].Particularly for isolated small power systems and microgrids (MGs), the low-inertia feature adversely affects system stability when distributed generation (DG) penetration increases and loads suddenly change [5].For instance, the abrupt rise in loads or disconnection of generation units may exacerbate the transient response of MGs, such as the rate of change of frequency (RoCoF) and frequency deviation, leading to system instability [6], and protection devices such as under-frequency load shedding may intervene to prevent a blackout of the whole MG [7].Many studies investigated the frequency control of microgrids integrated with a high share of RESs [8,9].
Various secondary controllers were developed in previous studies to improve power systems' frequency stability in addition to conventional controllers.The proportional integral derivative derivative (PIDD) controller optimized by the fruit fly algorithm provided better performance than the integral (I), proportional integral (PI), proportional integral derivative (PID), and integral derivative derivative (IDD) controllers for a deregulated power system [10].The frequency stability of a two-area power system is improved using a new cascaded controller optimized using artificial rabbit optimizers [11].The authors of [12] enhanced islanded MG frequency stability using a prairie dog optimizationbased cascaded controller.Reference [13] used a new cascaded controller optimized by Barnacle Mating Optimizer for improving the frequency stability of two-area interconnected MGs.Another study used a fuzzy cascaded controller to improve frequency control of an islanded MG containing several types of generators and energy storage [14].
In an attempt to cope with the lack of inertia concern, virtual synchronous generators (VSGs) or virtual inertia concepts are proposed [15][16][17][18].With the aid of these techniques, the power electronic devices coupled with RESs or energy storage systems (ESSs) are controlled to emulate the real SGs swing equation and deliver the necessary inertia support to the grid.Consequently, the system frequency stability could be improved, and the frequency fluctuations and RoCoF were reduced following disturbances [19].
Several control techniques have been developed in the literature based on the virtual inertia concept.A conventional PI approach has been applied in [20] to enhance the dynamic security of isolated MGs.Furthermore, a PI-based virtual inertia control (VIC) of a wind turbine has been implemented [21].In [22], a fuzzy-based virtual inertia strategy is introduced to control frequency deviation in a hybrid power system.Methodologies considering robust control techniques, such as the coefficient diagram method and H∞ have been proposed in [23,24] for MGs considering high RESs penetration.Other methodologies have adopted model predictive control [25] and adaptive VIC [26,27] to improve system frequency stability.
To attain the best performance from the above-mentioned controllers, they usually contain parameters that need to be optimally fine-tuned.However, selecting the optimal parameters for various virtual inertia frequency controllers in MGs is challenging, as it either employs trial-and-error methods or depends on the designer's experience.To determine the optimal frequency controller parameters, many researchers have implemented a variety of optimization techniques, including particle swarm optimization [28], a modified gray wolf optimization algorithm [29], a chaotic crow search algorithm [30], a Newton-based eigenvalue optimization algorithm [31], a knee point-driven evolutionary algorithm [32], a jellyfish search optimizer [33], and other recently developed optimization algorithms that have been successfully implemented in frequency control applications.
It can be noticed from the literature that the adoption of virtual inertia increases the system order, which causes the output active power to oscillate and adversely affects frequency dynamics.In addition, in the scope of fine-tuning the prior controllers, it is not easy to realize a reasonable trade-off in performance between robustness and control.Hence, it is challenging to guarantee robust stability and efficiency with the aforementioned controlling approaches under a wide range of demand and generation perturbations.
In recent days, the fractional-order (FO) controller has received a lot of attention due to its greater flexibility for adjusting system dynamics, specifically for systems operating in uncertain environments [34][35][36].It has been adopted in a deregulated environment for automatic generation control (AGC) [37] and load frequency control (LFC) [38].The FO controller's superiority is due to its two additional tuning knobs, such as the FO of the differentiator (μ) and the integrator (λ).Because of their inherent flexibility, several researchers have recommended FO controllers over traditional PI and PID controllers for power system stability applications [39,40].On the other hand, rare work has considered the application of FO derivative in VIC (FOVIC) to improve the conventional VIC (i.e., which uses integer order derivative) response, as in [41,42].This new concept needs further study on different power systems with different resources and operating conditions.Interestingly, the FOVSG was created in [41,42] with the presumption that the DC link is ideal; nevertheless, the DC link's energy storage has to be further included in the VIC model due to its limitations.In addition, the authors of [41,42], did not study the influence of variation of the virtual inertia constant with different values of FO operators.
Inspired by the above considerations, this work proposes a FOVIC that is applied to an islanded MG, including a small thermal power plant, wind power plant, solar power plants, ESS, and MG loads.The contributions of this research are summarized as follows: The rest of this paper is structured as follows: Section 2 introduces system configuration.In Section 3, the design of the proposed FOVIC control system is developed.Section 4 verifies the effectiveness of the suggested model and FOVIC through simulation results.The discussion is presented in Section 5.In Section 6, this paper's main conclusions are summarized.

System Modelling and Configuration
This study uses an islanded MG system, i.e., a test system, to examine the effectiveness of the proposed VIC model and the impact of the FO control strategy on its performance.The studied MG includes a dispatchable distributed generator (i.e., small thermal power plant), non-dispatchable distributed generators (i.e., wind and solar power plants), ESS, and MG loads [20].Moreover, the studied system contains a 20 MW thermal power plant, a 4 MW solar power plant, an 8 MW wind power plant, a 1 MW ESS, and 15 MW local loads.The studied island, MG, has a 20 MW power base.A schematic diagram of the investigated system with the proposed control strategy is revealed in Figure 1. Figure 2 shows the frequency response model of the studied system with the proposed developed control strategy, and the studied system's parameters are given in Table 1.
The fifth-order linearized system for the islanded MG, which considers high penetration of RESs, can be effectively modeled using the state-space approach.The deviation of the examined system's frequency can be obtained by considering the governor action (i.e., the primary control loop) and the LFC (i.e., the secondary control loop), as follows: where,  Table 1.Values of the studied MG parameters [20].In this study, the power changes caused by the wind (∆P ), solar (∆P ), , and load (∆P ) are taken into account as disturbance signals.Using the state variables from (1) to ( 5), the linearized state-space model of the investigated MG could be easily developed as follows:

Parameter
where, where W is the input perturbation vector, U is the control input signal, X is the state vector, and Y is the control output signal.A is the state matrix of the studied system.B and E correspond to the control input signal and the disturbance inputs.The output measurement, or input to the load-frequency controller, is represented by C. D and Z are zero vectors with the same size as the input control signal and disturbance vector, respectively.As a result, the islanded MGs complete state-space representation, taking into account RESs, can be obtained in Equations ( 8) and (9).

Description of the Fractional-Order Calculus
Despite being proposed 300 years ago, the utilization of FO has only recently become prevalent due to its complexity [36].Generic integral and differential notations into any actual number are possible by means of the fractional operators in the controller.The FO differentiator is a mathematical operator that can be viewed as a more generalized form of integral and differential operators, as indicated in the following way: where q is the FO operator and ub and lb are the upper and lower bands to calculate operator D.
Two different theories can define the FO principle.The Riemann-Liouville (RL) method is the original one and is utilized to identify the order derivative of a function f(t) [41,42] as follows: where () =    , ℜ() > 0 is the function of Gamma, n ∈ N and the variable q is limited as n-1< q <n.The previous fractional derivative of RL in (11) can then be transformed using the well-known Laplace method as follows [43]: where s is the Laplace operator.
The second definition related to the essential notation of FO is Caputo's definition, where the time-domain representation for q order of a function f(t) can be stated as follows: Once more, we apply the Laplace transformation to Equation (13) to result in Equation (14), which has an initial condition and represents the integral order with a specific physical significance.
The implementation of FO operators in the time domain involves complicated mathematical calculations.The recursive estimation method is typically used to implement the FO definition [44].The Laplace transformation of the q th derivative is as follows: where N is an effective frequency range approximation order [ ,  ] that can be chosen as [−1000, 1000] rad/s.
The fractional calculus expands on the conventional PID controller by allowing more variables to be adjusted.Unlike the traditional PID controller, which only has three variables that can be fine-tuned, the FOPID controller provides five variables for tuning.Consequently, there are two additional variables of integral and differential FOs λ and μ, respectively, as shown in Figure 3.When compared with the conventional PID, these variables have the capability to improve the controller's stability, transient time, and steadystate error.Moreover, it offers the controller greater flexibility and enables them to manage disturbances in the system that occur across a broad range of The transfer function of the FOPID controller is fully represented below: where λ and μ lie in the range of 0 and 1, which can decrease the steady-state error and rise time in the system.

Description of the Proposed Fractional-Order Virtual Inertia Control
The FO controller-based PID has been widely applied for power system stability due to its wide stability and robustness.In contrast, most of the studies are dedicated to LFC [38,45], and AGC [40].In all prior work in the literature, the notion of FOPID is used as a controller to support the LFC.The synchronous generator, which is regarded as the principal dynamic source of the conventional power system, contributes inertia to the grid via its rotating mass in the traditional power system [46], as in (17).
where G denotes the whole number of generating units coupled to the grid, SB is the power system's rated capacity, and Hx and SB,x are the inertia value and rated power of the generation units, respectively.Moreover, the swing equation can be employed to explain the rotating dynamics of actual synchronous machines, as in (18).
where ΔPM is the change in mechanical power, ΔPL is the load power change, and Δf is the frequency deviation.
To enhance the stability of low system inertia in the presence of a significant amount of RESs, it is possible to recreate the synthetic damping attribute and synthetic inertia power in the power system.This action can effectively improve the frequency stability of the power system.Thus, the inverter-based ESS, which is used to inject/absorb active power into/from the power system, is regulated based on the VIC, as in (19).
where KH is the virtual inertia constant of the VIC.The typical VIC system on the basis of battery ESS (BESS) is shown in Figure 4, which simulates the real synchronous generator's inertia.When the RESs are penetrated, the VIC-based ESSs dynamic equation that reflects the power system's desired power can be stated as in (20).The amount of power that the inertia of the BESS can generate is constrained by two factors: the maximum capacity of the BESS during the processes of charging and discharging, and the state of charge (SOC) of the BESS.These limitations prevent the inertia output power from exceeding the rated capacity of the BESS and avoid the lifetime degradation of the battery.
where  refers to the inverter-based BESS time constant.The proposed VIC model separates the DC-DC converter stage and the DC-AC inverter stage, as shown in Figure 5.Each stage is represented by a first-order model with different time constants (TDC for the DC-DC converter and TBESS for the DC-AC inverter) and conversion gains (KDC for the DC-DC converter and unity for the DC-AC inverter).This allows for more precise and accurate modeling of the individual stages.Typically, when the inertia power is not being emulated, the base power of the BESS is directly proportional to the deviation in frequency.This relationship can be represented mathematically as: The complete transfer function of the proposed VIC can be expressed as follows: where KH is the virtual inertia constant.Moreover, this study intends to show how the order of the derivative term of the inertia term can influence the power system's maximum frequency deviation, as shown in Figure 6.After analyzing Figure 6, it can be concluded that the optimal order of the s operator is 0.4.This value is crucial as it ensures the minimum frequency deviation is achieved at various points of KH.In other words, the s operator plays a significant role in maintaining stability and consistency in the system's frequency response.As we already indicated, the FO controller typically served as the secondary controller in most prior work.The FOVIC has the significant advantage of minimizing the system's order, leading to a substantial decrease in frequency fluctuations and output power oscillations.The complete structure of the proposed FOVIC is shown in Figure 7.The proposed FOVIC will undergo a second modification compared to Figure 4, which involves replacing the entire order of the inertial emulation with a fractional order.To achieve this, the output power and its sign will be combined with the required power from the DC-DC stage.The resulting power will then be set as a command for the inverter to supply power to the power system.The complete transfer function of the proposed FOVIC can be expressed as in (23): where μ is the order of the inertia derivative.

Simulation Results
To simulate the researched MG, a MATLAB/Simulink program is employed.With the proposed developed VIC model, which is depicted in Figure 2, and the system's parameters listed in Table 1 [20], the simulation results have been extracted based on two case studies, as follows: A. Validate the proposed developed VIC model's superiority compared with the conventional known VIC model, which is presented in many studies, e.g., [2,8,9,12,14,16].
Impact of the FO derivative on the developed structure of the suggested VIC model Both case studies are investigated through different scenarios under various situations of RES, load variation, and system inertia changes, as follows:

A1. Evaluation of the developed VIC model's performance under sudden load changes and different values of the virtual inertia constant
In this case, the suggested developed VIC scheme is compared with both the known conventional VIC scheme and without VIC (i.e., using conventional frequency control only) in the studied islanded MG under different values of virtual inertia constant as well as a 10% step change in load that occurs at time = 0 s. Figure 8 and Table 2 exhibit the frequency deviation of the investigated system with diverse control techniques under the nominal system parameters.From this simulation result, the studied MG with the suggested developed VIC scheme is noticeably more stable and quicker than the conventional VIC scheme.Moreover, the studied system with conventional frequency control (i.e., without the VIC scheme) has high overshoots and slow responses compared to the case of the proposed developed VIC scheme.In contrast, although the frequency deviation of the system with VIC schemes (i.e., the proposed and conventional schemes) is more suppressed when the KH value is increased, the system with the known conventional VIC needs more time to achieve stability.Thus, by applying the proposed developedVIC scheme, the dynamic stability of the studied MG considering different values of KH has been preserved.

A2. Evaluation of the developed VIC model's performance under several load step changes
In this case, a series of sudden load changes are applied in the investigated MG system with varied values of the virtual inertia.The step load change during the period of the simulation is revealed in Figure 9. Figure 10 demonstrates the frequency deviation of the studied MG with different control techniques under this disturbance.The results prove that the suggested developed VIC scheme gives the best performance for the studied islanded MG compared to the other control techniques, and the MG frequency variations are kept within the recommended levels during the load changes.Regarding the use of different values of the virtual inertia constant, the proposed developed VIC scheme showed fast stabilization time and slight deviations in system frequency compared to the conventional VIC scheme, which needs more time to stabilize.

A3. Evaluation of the developed VIC model's performance under step load change and RES presence
In this case, the MG with the suggested developed VIC scheme is studied under the following operational circumstances: a series of abrupt load changes shown in Figure 11 and power fluctuations of the RESs depicted in Figure 12. Figure 13 shows the frequency response of MG with the different control schemes, considering different values of virtual inertia.Figure 13 shows that the considered islanded MG with the suggested developed VIC scheme is more efficient and dependable.Additionally, in comparison with the conventional VIC scheme, the suggested scheme results in a better frequency response throughout all disturbances and operating conditions of the MG.On the other hand, conventional frequency control (i.e., without the VIC scheme) provided the poorest performance with high values of frequency deviation.Hence, with the proposed developed VIC scheme, the islanded MG system's overall performance is at its finest.

A4. Evaluation of the developed VIC model's performance considering communication delay.
Communication time delays are widely acknowledged in control systems for their potential to diminish system performance and induce instability.Consequently, these delays pose a significant challenge, emerging as a noteworthy uncertainty in analyzing the frequency stability of the power systems, given their ongoing expansion and increasing complexity.The effect of the communications delays with the proposed VIC strategy is tested at different time delays, as shown in Figures 14 and 15.This test considers both the communication time delay affecting the control input (i.e., before the LFC) and the delay in the control action (i.e., after the LFC).The effectiveness of the proposed VIC is evident, as it results in a lower frequency deviation compared to both the conventional VIC strategy and the LFC without VIC implementation.Additionally, with the proposed VIC, the system frequency achieves steady-state more rapidly than with the other two strategies.Notably, as the time delay increases from 0.1 s in Figure 14 to 0.2 s in Figure 15, both the conventional VIC and no VIC strategies exhibit heightened overshoot and settling time in frequency deviations.Consequently, applying the proposed VIC effectively mitigates the impact of communication delays without compromising the system's stability.

B1. Evaluation of the proposed FOVIC model's performance under sudden load change
In this case, the impact of the FO control technique has been considered in the suggested structure of the developed VIC scheme.Therefore, the MG with the suggested FO-VIC scheme is examined under different values of virtual inertia constant, different values of FO operators (μ) (i.e., 0.2, 0.4, 0.6, 0.8, and 1), and a 10% step change in load that occurs at time = 0 s. Figure 16 and Table 3 display the frequency deviation of the investigated system with FOVIC under the nominal system parameters.From Figure 16, it is clear that the higher the value of KH, the more the MG with the proposed FOVIC can dampen the frequency deviations, but it needs more time to achieve stability.On the other side, the role of the FO operator in the suggested FOVIC is shown to change the system's dynamic performance with respect to the overshoot, undershoot, and settling time.Therefore, the results indicate that the MG with the suggested FOVIC scheme at the fractional operator (μ) equals 0.4 is noticeably more stable and quicker than that designed with other fractional operator values.Moreover, it is obvious that this value is consistent with the range of values previously mentioned in the theoretical analysis section.

Virtual inertia constant=4
The system inertia (i.e., H) drops when the penetration level of RESs rises by replacing the synchronous generators, causing high-frequency deviations that may cause instability in the power system.In this case, the MG with the suggested FOVIC scheme is examined under different values of virtual inertia constant (KH), different values of FO operators (μ), and a 10% step change in load that occurs at time = 0 s.Moreover, the performance of the proposed FOVIC in the investigated system is evaluated by considering system inertia changes by ± 25% and ± 50%.show the frequency response of the MG with system uncertainties (i.e., increasing/decreasing the system inertia by ± 25% and ± 50%, respectively).According to the results, the suggested FOVIC scheme with different values of fractional operators significantly improves the MG frequency performance.It has lower system transients compared to the integer-order VIC.Furthermore, although the frequency variation of the system with the FOVIC scheme is more suppressed when the value of the virtual inertia is increased, the system may need more time to achieve stability.Hence, the proposed FOVIC scheme still maintains the stability of the MG frequency under uncertainties (i.e., system inertia changes).

B3. Evaluation of the proposed FOVIC model's performance under step load change and RESs
In this case, the MG with the suggested FOVIC scheme is investigated under different values of virtual inertia constant, different values of FO operators (μ), applying a series of abrupt load changes shown in Figure 21, and power fluctuations of the RESs shown in Figure 22. Figure 23 shows the frequency response of the MG with the suggested FOVIC scheme considering varied values of fractional operators.Out of the simulation results, the considered islanded MG with the proposed FOVIC scheme with different values of fractional operators is more efficient and dependable than the integer-order VIC.Additionally, compared to the proposed integer-order VIC scheme, the proposed FOVIC scheme results in a better frequency response throughout all circumstances of this case study.Although the frequency variation of the MG with the FOVIC scheme is more suppressed when the value of KH is increased, the system may need more time to achieve stability.Hence, the proposed FOVIC scheme still maintains the stability of the MG frequency under the high penetration of RESs.

Discussion
The proposed VIC model's developed structure can provide a virtual inertia facility by controlling the battery ESS (BESS) in a way that accurately mimics the behavior of a traditional synchronous generator during the emergency period.Moreover, to provide the proposed VIC model's developed structure with more flexibility in handling various microgrid disturbances, a fractional-order derivative control has been considered.In addition to enhancing frequency response, the proposed fractional-order virtual inertia control based on the BESS permits a greater level of renewable power penetration into the microgrid.The significance of the proposed control system has been supported by the findings obtained from the previous section.1.The proposed VIC model handles the sudden change in load demand better than the conventional VIC system by about 45% and the secondary frequency control (i.e., without VIC) system by about 80%.As a result, the proposed VIC model's developed structure has been chosen over alternative solutions for microgrid stabilization during sudden/series load changes.
2. The performance of the three control systems used-the proposed VIC model, the conventional VIC model, and the secondary load frequency control, i.e., without VIC-is evaluated under high RESs penetration as well as varying values of the system inertia constant.The system with secondary frequency control (i.e., without the VIC model) gives high oscillations that threaten the system's stability.Nonetheless, the proposed VIC and conventional VIC models are both able to stabilize the system frequency within an acceptable range.Among these two, the proposed VIC system exhibits the best frequency response under different values of the virtual inertia constant, with a maximum frequency deviation of ± 0.059 Hz (in comparison with ±0.11 Hz by the conventional VIC system).

Communication time delay:
With the proposed VIC model, the best frequency nadir (about 0.06 Hz) was recorded.In comparison, the conventional VIC model exhibits a frequency drop of about 0.12 Hz, whereas the load frequency control displays about 0.33 Hz.Because of this, the proposed VIC model works best in situations where there is a varying random time delay within the range of (0, 0.2), as in the case of the practical system.
4. The performance of the proposed VIC based on fractional-order derivative control has been evaluated under load/RESs fluctuations, system uncertainties (i.e., increasing/decreasing the system inertia by ± 25% and ± 50%), and different values in the virtual inertia constant and FO operators.The results show that the microgrid with the proposed FOVIC scheme at the fractional operator equal to 0.4 is significantly more stable and faster than that designed with other fractional operator values.Furthermore, it is clear that this value falls within the range of values that were previously mentioned in the section on theoretical analysis.

Conclusions
The decreasing inertia makes modern power systems more prone to system instability.The isolated microgrids (MGs) are more susceptible to instability due to the high penetration of non-dispatchable renewable energy sources, low inertia, and small number of dispatchable generation units.Recent studies proved the efficacy of virtual inertia control (VIC) in enhancing the frequency stability of low-inertia power systems.Therefore, this study proposes a new VIC model to accurately mimic inertia power based on an energy storage system (ESS).The new VIC model considers the effect of both the DC-DC converter and the DC-AC inverter on the power of the ESS used.The obtained results prove that the proposed new VIC scheme provides fast stabilization times and slight deviations in system frequency compared to the conventional VIC schemes for all MG operating conditions and disturbances.Moreover, this study proposes a fractional-order VIC to provide greater flexibility in dealing with MG disturbances.The obtained results prove that the proposed VIC outperforms the conventional LFC by about 80% and the conventional VIC model by about 45% in tackling loads/RESs fluctuations and system uncertainty.Hence, the proposed FOVIC scheme is noticeably more stable and faster in response compared to that designed with integer-order VIC for all MG operating conditions and disturbances.

Figure 1 .
Figure 1.A schematic diagram of the studied islanded MG with the proposed FOVIC strategy.

Figure 2 .
Figure 2. A LFC model of the studied islanded MG with the proposed developed FOVIC.

Figure 3 .
Figure 3.The basic structure of the FOPID controller as an example of applying FO calculus.

Figure 4 .
Figure 4.The dynamic structure of conventional VIC-based ESS.

Figure 5 .
Figure 5.The proposed VIC-based ESS dynamic structure.

Figure 6 .
Figure 6.Effect of the virtual inertia constant and FO operator on the system frequency deviation.

Figure 7 .
Figure 7.The proposed FOVIC to support the dynamic stability of MG.

Figure 8 .
Figure 8.The frequency response of the MG with different values of KH for case A1.

Figure 10 .
Figure 10.The frequency response of the MG with different values of KH for case A2.

Figure 12 .
Figure 12.Renewables power generation profile for case A3.

Figure 13 .
Figure 13.The frequency response of the MG with RES generation and different values of KH for case A3.

Figure 14 .
Figure 14.Effect of the 0.1s communications delay before and after LFC on the MG frequency response for case A4.

Figure 15 .
Figure 15.Effect of the 0.2s communications delay before and after LFC on the MG frequency response for case A4.

Figure 16 .
Figure 16.The frequency response of the MG with different values of KH and fractional operators for case B1.

Figure 17 .
Figure 17.The frequency response of the MG with different values of KH and fractional operators for case B2 with +25% system inertia (H) variation.

Figure 18 .
Figure 18.The frequency response of the MG with different values of KH and fractional operators for case B2 with -25% system inertia (H) variation.

Figure 19 .
Figure 19.The frequency response of the MG with different values of KH and fractional operators for case B2 with +50% system inertia (H) variation.

Figure 20 .
Figure 20.The frequency response of the MG with different values of KH and fractional operators for case B2 with -50% system inertia (H) variation.

Figure 23 .
Figure 23.The frequency response of the MG with RES generation, different values of KH, and fractional operators for case B3.

Table 2 .
Values of performance indices for case A1.

Table 3 .
Values of performance indices for case B1.