High-Level Renewable Energy Integrated System Frequency Control with SMES-Based Optimized Fractional Order Controller

: The high-level penetration of renewable energy sources (RESs) is the main reason for shifting the conventional centralized power system control paradigm into distributed power system control. This massive integration of RESs faces two main problems: complex controller structure and reduced inertia. Since the system frequency stability is directly linked to the system’s total inertia, the renewable integrated system frequency control is badly affected. Thus, a fractional order controller (FOC)-based superconducting magnetic energy storage (SMES) is proposed in this work. The detailed modeling of SMES, FOC, wind, and solar systems, along with the power network, is introduced to facilitate analysis. The FOC-based SMES virtually augments the inertia to stabilize the system frequency in generation and load mismatches. Since the tuning of FOC and SMES controller parameters is challenging due to nonlinearities, the whale optimization algorithm (WOA) is used to optimize the parameters. The optimized FOC-based SMES is tested under ﬂuctuating wind and solar powers. The extensive simulations are carried out using MATLAB Simulink environment considering different scenarios, such as light and high load proﬁle variations, multiple load proﬁle variations, and reduced system inertia. It is observed that the proposed FOC-based SMES improves several performance indices, such as settling time, overshoot, undershoot compared to the conventional technique. parameters with WOA and facilitate several case studies. The system and the designed controllers are simulated in MATLAB Simulink environment considering light and heavy load disturbances, ﬂuctuating solar and wind powers, and reduced system inertia to test the effectiveness and robustness. The simulated results conﬁrm the promising performance in reducing system frequency deviations and improving frequency stability of the system. The proposed FOC-based controller is superior over conventional SMES controller employing feedback and proportional gains in reducing settling time, overshoot and undershoot as evident from the analysis. Moreover, the simulation outcomes demonstrate the potential beneﬁts of FOC-based energy storage in high-level renewable energy integration and endorse the green effort to improve the sustainability. Finally, the large-scale DFIG offshore wind farm detail model with FOC-based SMES virtual inertia controller can be studied as future work.


Introduction
Environmental concern and increasing energy demand lead to an extensive penetration of renewable energy sources (RESs) into distribution networks that force to reshape power systems where conventional fossil-fuel-based power plants are being replaced by RES rapidly.These RESs, being pollution-free, cheap, and conforming to the concept of sustainable development, gained considerable attention worldwide in the last two decades.Many countries have effectively integrated a large share of RES into the primary grid, and a high penetration level is targeted for the next three decades.For example, China targets 69%, whereas that for EU 74%, India 75%, and USA 63% from RES by the year 2050 [1].

Issues with Frequency Deviation
The most promising RESs are the solar and wind power generations are due to the available resources, lower cost of power generation, and maximum power point tracking capability over a wide range of wind and sunlight variation [2,3].The inherent volatility and uncertainty characteristics with these resources, these sources are generally connected via power electronics-based inverter/converters to the power network.Their interaction with the grid is substantially different from that of the conventional synchronous generator (SG)-based plants that use steam and hydro turbines.While intrinsic kinetic energy in the SG rotor inherently provides inertia to the system, this is not the case for the RES connected via power electronics converters.Consequently, during disturbances and supply/demand imbalances, the inertia, which slows down the natural reaction of the system and provides time for the responsible controllers to take actions, is significantly reduced because the rate of change of frequency (RoCoF) is much higher in systems with low inertia [4].This high level of RoCoF and large frequency deviation (∆f) can provoke the tripping of sensible loads, generating units and relays, thus affecting system frequency stability even at small load-generation mismatch [5,6].Although the variable wind turbines have inertia, they are effectively decoupled from the system due to the power electronics-based converter interface; thus, it cannot improve frequency response.
Similarly, solar PV plants also do not have any rotating parts, thus cannot provide any inertia to the power system.As a result, the high-level integration of RESs reduces the system's total inertia due to the replacement of conventional synchronous generators.Moreover, the reduction in reserve power due to the replacement of reserve generating units causes frequency deviation [7].With the RES level increase, the power oscillations and frequency stability issues increase under disturbances [8], which may set the critical limit for RES share [9].To promote ambitious RES penetration levels in the grid, some form of mitigation scheme is required, and this has received significant attention from the researchers.

Mitigation Schemes
To address the frequency excursion issues with the low inertia system, several techniques are presented in the literature, such as auxiliary load frequency control (LFC) technique, inertia emulation technique, deloading technique, droop technique, and energy storage-based techniques [10][11][12].Among these, LFC is most extensively used for several decades to meet two main objectives: maintaining the system frequency and tie-line power deviation within specified values [13].Because of fluctuating active power generation, integration of wind/PV system makes load frequency control more challenging.The incorporation of the energy storage system (ESS) in the LFC loop demonstrated a useful technique for enhancing frequency stability [13][14][15].With proper control, additional inertia can be imitated for frequency stability and resiliency [16,17].Among the various ESS studied, such as superconducting magnetic energy storage (SMES), supercapacitors, electric batteries, fuel cell, flywheel energy storage, SMES technology is recognized as the best candidate for inertia emulation [18][19][20].
Performance of emulated inertia and LFC system vastly depends on the control technique employed.Recently, various control techniques have been suggested for the LFC of the RES integrated power system.In [21], frequency control of the interconnected power system, including wind farms, is presented where the symbiotic organism search algorithm optimized PID parameters.However, in this study, no energy storage system is considered.In [22], a model predictive control (MPC)-based synthetic inertial control is proposed for a wind farm consisting of wind turbines (WTs) and a battery storage energy system (BESS).In [23], battery energy storage is used for frequency support of a doubly fed induction generator (DFIG)-based wind system.The battery is connected to the DC link of DFIG and controlled with the droop technique to reduce frequency deviation by scheduling active power exchange during system disturbances.In [24], supercapacitor energy storage system is proposed to emulate the dynamic inertia response of PV-based ac microgrid cluster.The detailed design of feedback and proportion gains; however, are not discussed in this work.
In [25], an auxiliary LFC technique is presented to control the frequency of the Egyptian grid considering high-level PV and wind integration employing a proportionalintegral-derivative (PID) controller.However, the LFC technique does not consider the detailed Egyptian grid model; also, it excludes tie-line power that needs further investigation.In [26], a virtual inertia support technique is presented for low inertia microgrid with particle swarm optimization (PSO)-based PI controller.In [14], SMES-based virtual inertia control is proposed for a microgrid system.The conventional derivative approach virtual inertia control loop is implemented.The detailed design of feedback and proportion gains; however, are not discussed in this work.In [27], a self-adaptive virtual inertia fuzzy controller is adopted for the high-level renewable integrated system.The proportional virtual gain is adapted by the fuzzy system, which uses the deviation of real power and frequency as inputs.In this scheme; however, the generalized energy storage is considered by a simple 1st order system.The sharing of active power from different energy storages are scheduled based on their capabilities in [28] for frequency control of renewable sources.In this capability coordinated frequency control (CCFC) approach, the total error signal is forwarded to the primary control loop of each unit based on their capabilities.The LFC for mass-less inertia PV systems is presented in [29] with PI controllers.The parameters are optimized with the PSO-WOA hybrid optimization technique in case of different step load changes.To stabilize a low inertia PV system, another virtual inertia synthetization of synchroconverter is reported in [30] with a machine learning technique.The optimized virtual inertia frequency control and protection schemes are developed in [31,32] for low frequency interconnected power systems.A combined SMES and thyristor controlled phase shifters (TCPS) [33] are applied in low inertia utility grid with an adaptive neuro-fuzzy inference system (ANFIS) controller.The detailed design of SMES negative feedback and proportional gains, however, is not considered.

Research Gaps and Contributions
In recent years, fractional calculus gained popularity and is recognized as a powerful tool in the control engineering field.It is proved that the fractional order approach can provide better overall closed-loop system performance than the integer order control [34].This is attributed to the possibility of a more robust tuning of fractional order controller (FOC) in the presence of additional real parameters [35,36].Analytic, numerical, rule-based, and self-tuning and auto-tuning methods are commonly used for tuning fractional order PI/PID controllers [37].For superior performance, different optimization methods have been proposed in the past like-genetic algorithm (GA) [38], adaptive genetic algorithm (AGA) [39], particle swarm optimization (PSO) [39], modified PSO (MPSO) [40], differential evolution (DE) [41], electromagnetic-like (EM) algorithm [41], non-dominated sorting genetic algorithm (NSGA) [42], artificial bee colony (ABC) algorithm [43], seeker optimization algorithm (SOA) [44], harmony search (HS) algorithm [45], Tabu Search (TS) algorithm [46].Each method employs the objective function of varying degrees of complexity and has a different convergence speed.
Based on the above critical review of several recent works, especially the works focusing on SMES or other energy storage applications in low inertia system, it is identified that the detailed design of SMES is missing with FOC to support virtual inertia for RESs.The integer order controller is not adequate to control the dynamic behavior of the power system integrating large-scale RESs.Thus, in this work, we propose an optimal design approach for the fractional order PI controller and SMES.The dynamic model of the system with SMES and fractional order PI controller is introduced, and parameters are tuned with the WOA algorithm based on frequency deviation cost function.The WOA has several advantages that make it attractive for its applicability to a wide range of optimization problems [47].Among these are (1) WOA is a gradient-free method, like many other population-based algorithms.This eliminates the calculation of the gradient and stepsize at each iteration of the optimization process; (2) the WOA has inherent adaptive exploration and exploitation mechanism which reduces the probability of getting trapped in local solution; (3) it is insensitive to the initial solution(s), which may significantly affect the convergence and performance of the traditional methods; (4) it is easy to implement and flexible.To the best of the authors' knowledge, the main contribution of the proposed optimized FOC-based SMES controller in RESs integrated system can be summarized as below.

•
The inertia of the system is supported virtually, which makes the system stable over a wide range of load-generation mismatch.

•
The system frequency deviation is greatly improved with the proposed approach.

•
The proposed optimized FOC-based SMES approach is robust against system parameter variations.

•
The overshoot, undershoot, and settling time of the response are improved compared to the conventional approach.

•
The proposed approach endorses the green effort to augment sustainability.
The rest of the paper is organized as follows.Section 2 describes the dynamic model of the PV/wind integrated system to facilitate analysis.The detailed modeling of SMES and the proposed controller design technique are presented in Section 3. The simulated results and discussions are presented in Section 4 with several scenarios.Finally, the outcomes and contributions of this work are summarized in Section 5.

System Configuration
A two area of wind/PV integrated hybrid power system, as shown in Figure 1, is considered for the frequency control study.Each area comprises the conventional thermal generating unit, solar PV, industrial load, and superconducting magnetic energy storage (SMES).The fractional order PI controller is designed for the SMES to virtually support inertia, thereby mitigating the high-frequency excursion problem caused by generation load mismatch.As mentioned earlier, this low system inertia is a consequence of high-level renewable energy penetration.A tie-line interconnects these two areas.Measured frequency deviation and tie-line power signals are accumulated in the control and monitoring center, and then it sends appropriate control signals to the controllable SMES of both areas.

System Dynamic Modeling
To analyze frequency control in the presence of RESs, adequate dynamic model of the studied system is constructed.For precise demonstration of the dynamic behavior of interconnected system, generally higher order nonlinear models for thermal generating units, wind system, solar PV, converters are considered.For large power system with power electronic converters, however, the simplified dynamic models are employed to study the frequency stability.The interested readers can find more details on such dynamic modeling in [14,48].For frequency stability analysis, the simplified dynamic model of the two area system can be developed as shown in Figure 2. Renewable energy integrated two area power system.

System Dynamic Modeling
To analyze frequency control in the presence of RESs, adequate dynamic model of the studied system is constructed.For precise demonstration of the dynamic behavior of interconnected system, generally higher order nonlinear models for thermal generating units, wind system, solar PV, converters are considered.For large power system with power electronic converters, however, the simplified dynamic models are employed to study the frequency stability.The interested readers can find more details on such dynamic modeling in [14,48].For frequency stability analysis, the simplified dynamic model of the two area system can be developed as shown in Figure 2.
Power system of Area-1 Power system of Area-2

+ Industrial and Other Loads
Industrial and Other Loads - ( ) ( ) As the RESs generate energy from intermittent natural resources, the output power generated by them is variable.To reflect the realistic power output of fluctuating profile, white noise is added to the wind and PV system model, as shown in Figure 2.
For a wind turbine (WT) with a given physical dimension, the power output (PWT) mainly depends on wind velocity (Vw), which is variable.In this model, the wind-speed is multiplied by the random speed fluctuation (white noise) to estimate the WT output power fluctuations (ΔPWT).
Similarly, a PV system harvests energy from solar irradiation.For a given solar panel, the output power (Ppv) mainly depends on the solar irradiation (I) and atmospheric tem- As the RESs generate energy from intermittent natural resources, the output power generated by them is variable.To reflect the realistic power output of fluctuating profile, white noise is added to the wind and PV system model, as shown in Figure 2.
For a wind turbine (W T ) with a given physical dimension, the power output (P WT ) mainly depends on wind velocity (V w ), which is variable.In this model, the wind-speed is multiplied by the random speed fluctuation (white noise) to estimate the WT output power fluctuations (∆P WT ).
Similarly, a PV system harvests energy from solar irradiation.For a given solar panel, the output power (P pv ) mainly depends on the solar irradiation (I) and atmospheric temperature (T).Like WT, solar irradiation (I) is multiplied by the white noise to estimate PV output power fluctuations (∆P pv ).
Usually, renewable power plants deliver a considerable amount of power to the system without taking part in frequency control.Thus, solar PV, wind, and different loads (e.g., residential and industrial) are modeled as disturbances in the dynamic model.For further details on dynamic models of wind/PV systems, interested readers are directed to the literature [48].
In an interconnected power system, several physical constraints, such as generation rate constraints (GRC) and generator dead band (GDB), introduce nonlinearity and influence dynamic responses [49].It is shown that the dynamic responses of the system experience larger overshoot and longer settling times of frequency and tie-line power oscillation when these constraints are considered.Thus, controllers designed without considering these constraints may not perform well in practical application and even may lead to instability.Therefore, both constraints, GRC and GDB, are taken into consideration in these studies to reflect the practical implementation case.
As shown in Figure 2, the frequency deviation dynamics of the k-th area can be written as follows. where where H k , and D k are the inertia constant, and damping constant respectively in area-k.∆P TH,k , ∆P SA,k , and ∆P WF,k are the incremental power of thermal unit, solar farm, and wind farm, respectively, in area-k.T t,k , T g,k , T wind,k , and T pv,k are the time constant of turbine, governor, wind turbine, and solar PV, respectively, in area-k.∆P SMES,k is the variation in power output of the SMES.While conventional power generation block adjusts its power output depending on the change in frequency, the SMES block exchange (injects/absorbs) power with the grid that depends on system frequency deviation.

SMES Model with FOC
The superconducting magnetic energy storage (SMES) devices are capable of exchanging huge amount of power within very short duration.Its unique features, such as fast response, high efficiency, and long lifetime, have drawn the attention to the researchers in power system application and regarded as a very promising device for power system dynamic stability enhancement.
The SMES consists of two main parts: (i) the superconducting coil which is kept under extremely low temperature, and (ii) the power conversion system (PCS)-consisting of inverter/rectifier circuits [33].The energy exchange between the superconducting coil and the AC grid is facilitated by the three-phase transformers as shown in Figure 3. Two cascaded 12-pulse bridges filter the harmonic contents of the signals.During normal steady-state operating condition of the power system, the SMES coil is charged from the grid to its present value within very short period of time.In its charged state, the SMES coil starts conducting DC current with nearly zero losses as the coil temperature is maintained at extremely low value.The resulting DC voltage across the inductor is given by [14] where  0 is the maximum voltage of the bridge circuit,  is the thyristor firing angle,   is the superconducting coil current, and   is the commutating resistance or damping resistor.Thus, DC voltage appearing across the superconducting coil can be controlled with the variation of firing angle, .When the  is greater than 90°, the energy stored in the superconducting coil is released to the grid.Whereas, superconducting coil charges if the  is below 90°.In this way, inversion and rectification mode of converter discharges and charges, respectively, the superconducting coil almost instantly.It is also worth mentioning that the current through the coil has always positive direction, only altering voltage, either negative or positive, across the coil initiates the inversion or rectification mode.Now, during contingencies, as the power demand is initiated by the power system, the SMES discharges the stored energy through the PCS to the grid almost instantly, even before the governor and other supplementary control units start acting.After the contingencies, while the governors start working to support the power demand, the SMES again starts charging its preset value.
The detailed dynamic model of SMES for frequency stability studies along with the fractional order PI controller is shown in Figure 4.During excessive system loading, load surpasses the generation, the   becomes negative while the current   maintains the same direction.The incremental change in   can be written as where   is the SMES gain, E  is the output of fractional order PI controller,   is the negative feedback gain,   is the incremental change in superconducting coil current,   is the converter delay time.The incremental change in inductor current   is written as During normal steady-state operating condition of the power system, the SMES coil is charged from the grid to its present value within very short period of time.In its charged state, the SMES coil starts conducting DC current with nearly zero losses as the coil temperature is maintained at extremely low value.The resulting DC voltage across the inductor is given by [14] where V d0 is the maximum voltage of the bridge circuit, α is the thyristor firing angle, I D is the superconducting coil current, and R D is the commutating resistance or damping resistor.Thus, DC voltage appearing across the superconducting coil can be controlled with the variation of firing angle, α.When the α is greater than 90 • , the energy stored in the superconducting coil is released to the grid.Whereas, superconducting coil charges if the α is below 90 • .In this way, inversion and rectification mode of converter discharges and charges, respectively, the superconducting coil almost instantly.It is also worth mentioning that the current through the coil has always positive direction, only altering voltage, either negative or positive, across the coil initiates the inversion or rectification mode.Now, during contingencies, as the power demand is initiated by the power system, the SMES discharges the stored energy through the PCS to the grid almost instantly, even before the governor and other supplementary control units start acting.After the contingencies, while the governors start working to support the power demand, the SMES again starts charging its preset value.
The detailed dynamic model of SMES for frequency stability studies along with the fractional order PI controller is shown in Figure 4.During excessive system loading, load surpasses the generation, the E D becomes negative while the current I D maintains the same direction.The incremental change in E D can be written as where K SMES is the SMES gain, ∆E is the output of fractional order PI controller, K ID is the negative feedback gain, ∆I D is the incremental change in superconducting coil current, T DC is the converter delay time.The incremental change in inductor current I D is written as The active power deviation of SMES is given by Electronics 2021, 10, x FOR PEER REVIEW 9 of 22 The active power deviation of SMES is given by

Fractional Order PI-Based SMES Controller Design
Fractional order PI-based SMES controller is proposed to improve frequency stability of low inertia power system considered in this study.The fractional order calculus involves generalized integration and differentiation of non-integer order, and these mathematical phenomena allow describing a real system more precisely than the classical integer order methods, especially when the dynamic system is of distributed parameter nature [50,51].As the fractional order controllers tuning is reaching to a matured state of practical use, their domination in the industry is increasing rapidly due to superior performance over conventional integer order controller.
The time domain fractional order PI controller can be represented as where () is the error signal,   is the proportional gain,   is the integral gain,  is the fractional order and a real number which lies between 0 to 2. The Laplace transformation gives the following transfer function for the fractional order PI controller.
( ) , (0, 2) The conventional integer order PI and fractional order PI can be understood by Figure 5 in  axis.The integer order controller is represented by two fixed points on the  axis.Whereas, the fractional order PI controller can be represented by infinite number of points between 0 and 2. Thus, it gives more degree of freedom and flexibility over the conventional integer order controller.

Fractional Order PI-Based SMES Controller Design
Fractional order PI-based SMES controller is proposed to improve frequency stability of low inertia power system considered in this study.The fractional order calculus involves generalized integration and differentiation of non-integer order, and these mathematical phenomena allow describing a real system more precisely than the classical integer order methods, especially when the dynamic system is of distributed parameter nature [50,51].As the fractional order controllers tuning is reaching to a matured state of practical use, their domination in the industry is increasing rapidly due to superior performance over conventional integer order controller.
The time domain fractional order PI controller can be represented as where e(t) is the error signal, K p is the proportional gain, K i is the integral gain, λ is the fractional order and a real number which lies between 0 to 2. The Laplace transformation gives the following transfer function for the fractional order PI controller.
The conventional integer order PI and fractional order PI can be understood by Figure 5 in λ axis.The integer order controller is represented by two fixed points on the λ axis.Whereas, the fractional order PI controller can be represented by infinite number of points between 0 and 2. Thus, it gives more degree of freedom and flexibility over the conventional integer order controller.
Electronics 2021, 10, x FOR PEER REVIEW 9 of 22 The active power deviation of SMES is given by

Fractional Order PI-Based SMES Controller Design
Fractional order PI-based SMES controller is proposed to improve frequency stability of low inertia power system considered in this study.The fractional order calculus involves generalized integration and differentiation of non-integer order, and these mathematical phenomena allow describing a real system more precisely than the classical integer order methods, especially when the dynamic system is of distributed parameter nature [50,51].As the fractional order controllers tuning is reaching to a matured state of practical use, their domination in the industry is increasing rapidly due to superior performance over conventional integer order controller.
The time domain fractional order PI controller can be represented as where () is the error signal,   is the proportional gain,   is the integral gain,  is the fractional order and a real number which lies between 0 to 2. The Laplace transformation gives the following transfer function for the fractional order PI controller.
( ) , (0, 2) The conventional integer order PI and fractional order PI can be understood by Figure 5 in  axis.The integer order controller is represented by two fixed points on the  axis.Whereas, the fractional order PI controller can be represented by infinite number of points between 0 and 2. Thus, it gives more degree of freedom and flexibility over the conventional integer order controller.As presented in Figure 4, the SMES virtual inertia based on fractional order PI is developed in this study to support the frequency of the low inertia interconnected system.WOA optimizes the feedback and proportional gains of SMES and the fractional order PI's proportional gain, integral gain, and fractional parameter.The following subsections describe the formulation of objective function and solution method employing WOA.

Description of the Cost Function
The appropriate cost function is vital in the application of nature inspired and heuristic optimization techniques in power system.In general, the cost function is defined to minimize the several performance indices like overshoot, undershoot, and settling time for various signals.In this work, several fractional order PI gains, fractional orders, SMES feedback gains, and proportional gains are designed based on the area frequency deviation and tie-line power fluctuation cost function.For the better comprehension of the optimization process, the following cost function, decision variables, and constraints are considered. Minimize: Decision Variables: Constraints: where subscripts 1 and 2 are to denote area-1 and area-2 for the interconnected power system.T is the simulation time, ∆ f is the frequency deviation, ∆P tie is the tie-line power deviation.K p is the fractional order PI proportional gain, K i is the FOC integral gain, K 1 & K 2 are the Area Error Controller (ACE) integral gains, K ID is the SMES negative feedback gain, and K SMES is the SMES proportional gain.Mainly, the upper and lower limits of the equation 14 are selected based on knowledge/experience of FOC and SMES applications in power system.The optimization algorithm is coded in MATLAB script (.m file) environment and linked with the MATLAB Simulink (.slx files) environment.

Controller Design with WOA
The minimization problem described by equation 12 is solved by whale optimization algorithm (WOA).Because of its inherent adaptive exploration and exploitation capability, WOA outperforms state-of-the-art metaheuristic algorithms, e.g., PSO, GA, GSA, ACO.Recently, WOA is being applied in many large-scale and real-world engineering optimization problems [47,52].
WOA is a metaheuristic algorithm that mimics the unique hunting strategy of humpback whales.The whale searches for the prey (e.g., krill, small fishes), and then attacks them by a technique called bubble-net hunting.Entire hunting mechanism consists of three processes, (i) encircling prey (ii) bubble-net feeding method, and (iii) search for prey.These processes are mathematically modeled in WOA as described below [53].

a. Encircling prey
Upon locating the prey, the humpback whale encircles them and advances their position toward the best search agent.During the course of iteration this behavior is expressed mathematically as where k is the current iteration, → X is the current position vector, → X * is the best position of the solution obtained so far.

→
A and → C are coefficient vectors which are updated as follows where → r is a random vector in [0, 1], and → a is linearly decreased from 2 to 0 over the course of iterations.
, where I max is the maximum number of iterations.

b. Bubble-net attacking mechanism
In this phase, the humpback whale starts to push bubbles and move toward the prey by shrinking circle and spiral-shaped path.The position of the humpback whale and the prey to mimic the helix-shaped motion is given by with where l is a random number in [−1,1], and b is a constant defining spiral-shaped motion.A probability p of 0.5 is assumed to distinguish between bubbles attacking with shrinking behavior and spiral motions as given by where p is a random number in [0, 1].

c. Search for prey
To increase the diversity of solutions and achieve a global solution, the search space is extended by using → A > 1 .Thus, in prey search phase, the humpback whales are forced to move away from a random whale position.The representative equation for the prey search is given by where → X rand is the position of a whale selected randomly from the current population.The overall flow chart for the WOA to design fractional order PI and SMES parameters is shown in Figure 6.

Simulation Results
The effectiveness and robustness of the proposed optimized FOC-based SMES to improve the frequency stability are presented in this section.The system modeled in Section 2 is considered for analytical analysis.The system parameters and optimized controller parameters listed in Tables 1 and 2, respectively, are used to conduct time domain simulations and facilitate analysis.The MATLAB Simulink is used to conduct analyses considering different scenarios such as light loading, heavy loading, solar and wind power variations, and reduced inertia.The frequency deviations in both areas are plotted to show the effectiveness of the proposed FOC.The system dynamic model and optimization code are interlinked to find the optimized parameters.The convergence of the optimization algorithm depicted in Figure 7 shows that most of the runs in simulation converge before 30 iterations.

Simulation Results
The effectiveness and robustness of the proposed optimized FOC-based SMES to improve the frequency stability are presented in this section.The system modeled in Section 2 is considered for analytical analysis.The system parameters and optimized controller parameters listed in Tables 1 and 2, respectively, are used to conduct time domain simulations and facilitate analysis.The MATLAB Simulink is used to conduct analyses considering different scenarios such as light loading, heavy loading, solar and wind power variations, and reduced inertia.The frequency deviations in both areas are plotted to show the effectiveness of the proposed FOC.The system dynamic model and optimization code are interlinked to find the optimized parameters.The convergence of the optimization algorithm depicted in Figure 7 shows that most of the runs in simulation converge before 30 iterations.In this scenario, the studied dynamic model of the system is tested under the load profile variations in area-1.The load profiles are step changed in area-1 to observe the frequency response improvement with the proposed FOC-based SMES.The frequency deviations for area-1 and area-2 are plotted in Figure 7a,b for a step load change of 0.1 p.u. at 25 s during time domain simulation.The positive effect of the proposed controller is clearly visualized through the reduction of frequency deviations.As visualized in Figure 8a, a step load change in area-1 causes a significant frequency deviation in area-1 without virtual inertia controller.The frequency deviation is around 1 Hz without any auxiliary controller.The conventional controller-based SMES improves the deviation to about 0.07 Hz.However, the proposed FOC-based SMES greatly improves the frequency deviation in area-1 which is around 0.01 Hz.The maximum overshoot without any auxiliary controller is 0.5 Hz whereas it is only 0.07 Hz with conventional SMES.The proposed optimized FOC-based SMES has almost zero overshoot.It is noticed that the proposed FOC-based SMES is capable of improving all the performance indices such as maximum overshoot, maximum undershoot and settling time.Likewise, the frequency deviation in area-2 is very high, around 0.04 Hz, without any inertial controller as depicted in Figure 8b.The conventional SMES controller improves frequency deviation to some extent.However, the proposed optimized FOC-based SMES reduces the frequency deviation to almost zero.A large step load change (0.3 p.u.) is also applied in area-1 as shown in Figure 9.It is observed in Figure 9a,b that the system cannot remain in synchronism without any auxiliary controller.The frequency deviation in both areas continues to increase and the system becomes unstable.A conventional SMES system improves the frequency response keeping the frequency deviation within acceptable limit.However, a great improvement in system frequency response is observed with optimized FOC-based SMES controller which is clearly visualized in Figure 9a,b.Since the system does not return to original stable point without any auxiliary controller, the settling time is infinite.However, the settling time is improved with the conventional SMES.The proposed optimized FOC-based controller is capable to greatly improve the settling time of the frequency responses in both areas.A large step load change (0.3 p.u.) is also applied in area-1 as shown in Figure 9.It is observed in Figure 9a,b that the system cannot remain in synchronism without any auxiliary controller.The frequency deviation in both areas continues to increase and the system becomes unstable.A conventional SMES system improves the frequency response keeping the frequency deviation within acceptable limit.However, a great improvement in system frequency response is observed with optimized FOC-based SMES controller which is clearly visualized in Figure 9a,b.Since the system does not return to original stable point without any auxiliary controller, the settling time is infinite.However, the settling time is improved with the conventional SMES.The proposed optimized FOC-based controller is capable to greatly improve the settling time of the frequency responses in both areas.The load disturbances, ranging from low to high, are also applied in area-2 with the system default inertia.As depicted in Figure 10a, the frequency deviation in area-1 is 0.65 Hz without any virtual inertia controller for a step load change of 0.2 p.u. in area-1.The

Load Profile Variation in Area-2
The load disturbances, ranging from low to high, are also applied in area-2 with the system default inertia.As depicted in Figure 10a, the frequency deviation in area-1 is 0.65 Hz without any virtual inertia controller for a step load change of 0.2 p.u. in area-1.The conventional SMES controller reduces the frequency deviation to 0.1 Hz, whereas the proposed optimized FOC controller-based SMES is capable of maintaining almost zero frequency deviation.Similarly, the frequency deviation in area-2 is 0.55 Hz without any auxiliary controller.The conventional SMES controller is capable of reducing frequency deviation by 81.8%.However, the proposed optimized FOC-based SMES controller reduces the frequency deviation by 99.21%.It is noteworthy that the settling time is slightly increased for conventional SMES controller while the frequency deviation is improved.The proposed optimized FOC-based SMES reduces the settling time.The system frequency oscillates over a wide range (greater than 0.5 Hz) without any inertia controller in both areas.Since these oscillations are beyond the acceptable limits, it mandates the system frequency protection relay to operate.The proposed optimized FOC-based control scheme is capable of avoiding under frequency relay operation, thus, it improves the system reliability.

Load Profile Variation in Area-2
The load disturbances, ranging from low to high, are also applied in area-2 with the system default inertia.As depicted in Figure 10a, the frequency deviation in area-1 is 0.65 Hz without any virtual inertia controller for a step load change of 0.2 p.u. in area-1.The conventional SMES controller reduces the frequency deviation to 0.1 Hz, whereas the proposed optimized FOC controller-based SMES is capable of maintaining almost zero frequency deviation.Similarly, the frequency deviation in area-2 is 0.55 Hz without any auxiliary controller.The conventional SMES controller is capable of reducing frequency deviation by 81.8%.However, the proposed optimized FOC-based SMES controller reduces the frequency deviation by 99.21%.It is noteworthy that the settling time is slightly increased for conventional SMES controller while the frequency deviation is improved.The proposed optimized FOC-based SMES significantly reduces the settling time.The system frequency oscillates over a wide range (greater than 0.5 Hz) without any inertia controller in both areas.Since these oscillations are beyond the acceptable limits, it mandates the system frequency protection relay to operate.The proposed optimized FOC-based control scheme is capable of avoiding under frequency relay operation, thus, it improves the system reliability.For the high step load changes in area-2 (0.3 p.u.), the frequency deviations of both areas are greater than 1 Hz which activates the under frequency relay operating set point of 59.5 Hz without any virtual inertia controller.However, the frequency deviation is well below the under frequency relay operating point with the conventional SMES controller as depicted in Figure 11a,b.In these figures, it is clearly visualized that the proposed controller is capable of keeping the frequency deviations in both areas to almost zero.Thus, the system stability and reliability are guaranteed with the proposed fractional order controller-based SMES.
areas are greater than 1 Hz which activates the under frequency relay operating set point of 59.5 Hz without any virtual inertia controller.However, the frequency deviation is well below the under frequency relay operating point with the conventional SMES controller as depicted in Figure 11a,b.In these figures, it is clearly visualized that the proposed controller is capable of keeping the frequency deviations in both areas to almost zero.Thus, the system stability and reliability are guaranteed with the proposed fractional order controller-based SMES.In this scenario, multiple load variations as shown in Figure 12 are applied in both areas of the studied system to investigate the system's ability to bring back the frequency deviation to zero before the next changes.A better performance of the proposed optimized FOC-based SMES is clearly visible from the system frequency response.As seen in Figure 13a,b, following the first step load change of 0.075 p.u. at 25 s in area-1, the proposed controller is faster to eliminate the frequency deviation before beginning the second step load change of 0.15 p.u. at 50 s compared to the conventional techniques.The frequency deviations in area-1 are very high at all points of step changes without virtual inertia controller.Although the conventional SMES controller improves the frequency response slightly, a notable improvement is achieved with the proposed technique.In this case, the proposed method also provides In this scenario, multiple load variations as shown in Figure 12 are applied in both areas of the studied system to investigate the system's ability to bring back the frequency deviation to zero before the next changes.
areas are greater than 1 Hz which activates the under frequency relay operating set point of 59.5 Hz without any virtual inertia controller.However, the frequency deviation is well below the under frequency relay operating point with the conventional SMES controller as depicted in Figure 11a,b.In these figures, it is clearly visualized that the proposed controller is capable of keeping the frequency deviations in both areas to almost zero.Thus, the system stability and reliability are guaranteed with the proposed fractional order controller-based SMES.In this scenario, multiple load variations as shown in Figure 12 are applied in both areas of the studied system to investigate the system's ability to bring back the frequency deviation to zero before the next changes.A better performance of the proposed optimized FOC-based SMES is clearly visible from the system frequency response.As seen in Figure 13a,b, following the first step load change of 0.075 p.u. at 25 s in area-1, the proposed controller is faster to eliminate the frequency deviation before beginning the second step load change of 0.15 p.u. at 50 s compared to the conventional techniques.The frequency deviations in area-1 are very high at all points of step changes without virtual inertia controller.Although the conventional SMES controller improves the frequency response slightly, a notable improvement is achieved with the proposed technique.In this case, the proposed method also provides A better performance of the proposed optimized FOC-based SMES is clearly visible from the system frequency response.As seen in Figure 13a,b, following the first step load change of 0.075 p.u. at 25 s in area-1, the proposed controller is faster to eliminate the frequency before beginning the second step load change of 0.15 p.u. at 50 s compared to the conventional techniques.The frequency deviations in area-1 are very high at all points of step changes without virtual inertia controller.Although the conventional SMES controller improves the frequency response slightly, a notable improvement is achieved with the proposed technique.In this case, the proposed method also provides much better performance in terms of overshoot, undershoot, and settling time.The frequency responses as visualized in Figure 13c,d show the better performance with the proposed control technique for the multiple load changes in area-2.much better performance in terms of overshoot, undershoot, and settling time.The frequency responses as visualized in Figure 13c,d show the better performance with the proposed control technique for the multiple load changes in area-2.The effectiveness of the proposed controller is also tested with fluctuating solar and wind powers in both areas.The intermittent solar and wind powers considered in this study are depicted in Figure 14a,b, respectively.The solar and wind powers have the mean value of 0.195 p.u. and 0.10 p.u., respectively.The solar power is integrated in area-1 at 50 s during 150 s simulation time which continues to inject fluctuating power during the entire simulation period.On the other hand, the intermittent wind generating unit is connected at 75 s which is kept connected throughout the entire simulation period.As shown in Figure 14c,d, the connection of varying solar and powers has detrimental effect on system frequency response without any auxiliary controller.The frequency of the system continues to vary during entire simulation period and does not settle to steady state value.The conventional SMES controller slightly improves the system frequency response.On the other hand, the proposed controller performance is superior, in terms of

Frequency Response Analysis for Solar and Wind Power Variations
The effectiveness of the proposed controller is also tested with fluctuating solar and wind powers in both areas.The intermittent solar and wind powers considered in this study are depicted in Figure 14a,b, respectively.The solar and wind powers have the mean value of 0.195 p.u. and 0.10 p.u., respectively.The solar power is integrated in area-1 at 50 s during 150 s simulation time which continues to inject fluctuating power during the entire simulation period.On the other hand, the intermittent wind generating unit is connected at 75 s which is kept connected throughout the entire simulation period.As shown in Figure 14c,d, the connection of varying solar and wind powers has detrimental effect on system frequency response without any auxiliary controller.The frequency of the system continues to vary during entire simulation period and does not settle to steady state value.The conventional SMES controller slightly improves the system frequency response.On the other hand, the proposed controller performance is superior, in terms of settling time, overshoot, and undershoot, over the conventional SMES controller.The improvement of several performance indices are listed in Table 3 to demonstrate the superiority of the proposed controller.
settling time, overshoot, and undershoot, over the conventional SMES controller.The improvement of several performance indices are listed in Table 3 to demonstrate the superiority of the proposed controller.In this case, the robustness of the proposed controller is verified with the system inertia variations.The inertia in both areas is reduced by 50% and a step load change of 0.2 p.u. is applied in area-1and area-2 at 50 s.The frequency response for step load change in area-1 is depicted in Figure 15a,b.As depicted in Figure 8a,b, the system is capable of maintaining stable operation with a step load change of 0.2 p.u. in case of default inertia (100%).However, Figure 15a,b show that the frequency deviations in both areas gradually increase leading to instability.The conventional SMES controller reduces the frequency deviations and stabilizes the system.However, the proposed control method augments the system stability greatly by reducing frequency deviations to almost zero even with the 50% system inertia.Thus, the proposed controller is more robust compared to the conventional technique.The robustness of the proposed optimized FOC-based SMES is also depicted in Figure 15c,d for a 0.2 p.u. step load change in area-2.The frequency deviations in both areas are unbounded without any auxiliary controller.Although the conventional  In this case, the robustness of the proposed controller is verified with the system inertia variations.The inertia in both areas is reduced by 50% and a step load change of 0.2 p.u. is applied in area-1and area-2 at 50 s.The frequency response for step load change in area-1 is depicted in Figure 15a,b.As depicted in Figure 8a,b, the system is capable of maintaining stable operation with a step load change of 0.2 p.u. in case of default inertia (100%).However, Figure 15a,b show that the frequency deviations in both areas gradually increase leading to instability.The conventional SMES controller reduces the frequency deviations and stabilizes the system.However, the proposed control method augments the system stability greatly by reducing frequency deviations to almost zero even with the 50% system inertia.Thus, the proposed controller is more robust compared to the conventional technique.The robustness of the proposed optimized FOC-based SMES is also depicted in Figure 15c,d for a 0.2 p.u. step load change in area-2.The frequency deviations in both areas are unbounded without any auxiliary controller.Although the conventional SMES makes the system stable, the proposed FOC-based SMES is superior to stabilize the system after step load disturbance.

Discussions
In this subsection, the overall performance of the proposed controller is discussed and compared to the existing methods.Moreover, the several limitations of the proposed techniques are mentioned along with the possible solution techniques.The addition of this virtual inertia with the FOC-based SMES makes the system stable over a wide range of load-generation mismatch.Since the FOC is superior over the conventional controller, the proposed technique performs better in reducing system frequency deviation.Simulation results show that the proposed controller is capable of improving the settling time, maximum overshoot, and maximum undershoot of frequency deviations.
It is noteworthy that oscillations are beyond the acceptable limits in some cases without auxiliary devices and with conventional SMES.Thus, it mandates the system frequency protection relay to operate.However, the proposed technique is capable of avoiding under frequency relay operation, which improves the system reliability.The system shows robust performance with reduced inertia as discussed in Section 4.1.5.The main disadvantage of the proposed technique is that it does not consider the optimization of the number of SMES required to stabilize the system.Thus, further study is needed to analyze the cost-benefit of the proposed FOC-based SMES technique as well as reduce the number of SMES for large power system using any suitable modern optimization technique.

Discussion
In this subsection, the overall performance of the proposed controller is discussed and compared to the existing methods.Moreover, the several limitations of the proposed techniques are mentioned along with the possible solution techniques.The addition of this virtual inertia with the FOC-based SMES makes the system stable over a wide range of load-generation mismatch.Since the FOC is superior over the conventional controller, the proposed technique performs better in reducing system frequency deviation.Simulation results show that the proposed controller is capable of improving the settling time, maximum overshoot, and maximum undershoot of frequency deviations.
It is noteworthy that oscillations are beyond the acceptable limits in some cases without auxiliary devices and with conventional SMES.Thus, it mandates the system frequency protection relay to operate.However, the proposed technique is capable of avoiding under frequency relay operation, which improves the system reliability.The system shows robust performance with reduced inertia as discussed in Section 4.1.5.The main disadvantage of the proposed technique is that it does not consider the optimization of the number of SMES required to stabilize the system.Thus, further study is needed to analyze the cost-benefit of the proposed FOC-based SMES technique as well as reduce the number of SMES for large power system using any suitable modern optimization technique.

Conclusions
This work proposes a load frequency control technique with FOC-based SMES for highly renewable energy integrated low inertia power system.A combined dynamic model of FOC, SMES, and two area power system is developed to design optimized parameters with WOA and facilitate several case studies.The system and the designed controllers are simulated in MATLAB Simulink environment considering light and heavy load disturbances, fluctuating solar and wind powers, and reduced system inertia to test the effectiveness and robustness.The simulated results confirm the promising performance in reducing system frequency deviations and improving frequency stability of the system.The proposed FOC-based controller is superior over conventional SMES controller employing feedback and proportional gains in reducing settling time, overshoot and undershoot as evident from the analysis.Moreover, the simulation outcomes demonstrate the potential benefits of FOC-based energy storage in high-level renewable energy integration and endorse green effort to improve the sustainability.Finally, the large-scale DFIG offshore wind farm detail model with FOC-based SMES virtual inertia controller can be studied as future work.

Figure 1 .
Figure 1.Renewable energy integrated two area power system.Electronics 2021, 10, x FOR PEER REVIEW 6 of 22

Figure 2 .
Figure 2. The dynamic model of two area system with proposed controller.

Figure 2 .
Figure 2. The dynamic model of two area system with proposed controller.

Figure 4 .
Figure 4. Dynamic model of SMES along with FOC.

2 (Figure 5 .
Figure 5.The integer order and fractional order PI.

Figure 4 .
Figure 4. Dynamic model of SMES along with FOC.

Figure 4 .
Figure 4. Dynamic model of SMES along with FOC.

2 (Figure 5 .
Figure 5.The integer order and fractional order PI.Figure 5.The integer order and fractional order PI.

Figure 5 .
Figure 5.The integer order and fractional order PI.Figure 5.The integer order and fractional order PI.

Figure 6 .
Figure 6.The WOA flow chart for optimizing control parameter.

Figure 6 .
Figure 6.The WOA flow chart for optimizing control parameter.

Figure 7 .
Figure 7. Convergence of WOA algorithm for FOC-based SMES design.

Electronics 2021 ,
10, x FOR PEER REVIEW 14 of 22virtual inertia controller.The frequency deviation is around 1 Hz without any auxiliary controller.The conventional controller-based SMES improves the deviation to about 0.07 Hz.However, the proposed FOC-based SMES greatly improves the frequency deviation in area-1 which is around 0.01 Hz.The maximum overshoot without any auxiliary controller is 0.5 Hz whereas it is only 0.07 Hz with conventional SMES.The proposed optimized FOC-based SMES has almost zero overshoot.It is noticed that the proposed FOCbased SMES is capable of improving all the performance indices such as maximum overshoot, maximum undershoot and settling time.Likewise, the frequency deviation in area-2 is very high, around 0.04 Hz, without any inertial controller as depicted in Figure8b.The conventional SMES controller improves frequency deviation to some extent.However, the proposed optimized FOC-based SMES reduces the frequency deviation to almost zero.

Figure 13 .
Figure 13.Frequency response for both areas under multiple load profile variations.4.1.4.Frequency Response Analysis for Solar and Wind Power Variations

Figure 13 .
Figure 13.Frequency response for both areas under multiple load profile variations.

Figure 14 .
Figure 14.Frequency response for both areas under multiple load profile variations.

Figure 14 .
Figure 14.Frequency response for both areas under multiple load profile variations.

Electronics 2021 ,
10, x FOR PEER REVIEW 19 of 22SMES makes the system stable, the proposed FOC-based SMES is superior to stabilize the system after step load disturbance.

Figure 15 .
Figure 15.Frequency response for both areas for step load change with 50% reduced inertia.

Figure 15 .
Figure 15.Frequency response for both areas for step load change with 50% reduced inertia.

Table 3 .
Performance indices improvement with FOC-based SMES.

Table 3 .
Performance indices improvement with FOC-based SMES.