Heuristic Optimization of Virtual Inertia Control in Grid-Connected Wind Energy Conversion Systems for Frequency Support in a Restructured Environment

: In the work reported in this paper, a novel application of the artificial bee colony algorithm is used to implement a virtual inertia control strategy for grid-connected wind energy conversion systems. The proposed control strategy introduces a new heuristic optimization technique that uses the artificial bee colony (ABC) algorithm to calculate the optimal gain value of an additional derivative control loop added to the control scheme of the machine side converter in a wind energy system to enable wind farms to participate in frequency control as specified by recent grid codes. This helps to minimize the frequency deviations, reduce active power deviation in the system, and increase the penetration level of wind energy in power systems. The study was performed in a restructured power system environment. The proposed control scheme and its robustness were evaluated using load–frequency analysis for three real-life transaction scenarios that can occur in an interconnected open-energy market and the validation was carried out using eigenvalue analysis. The results in this study show that the optimal gain of the proposed controller reduces the frequency deviations and improves stability and overall performance of the system.


Introduction
To meet increasing load demand and reduce environmental pollution caused by fossil-fuelled power plants, power generation from renewable sources has become a viable solution in many instances. Modern power systems consist of many control areas with different sources of generation [1]. Any variance between the generated power and required load demand results in the deviation of the system frequency from its nominal value. This also creates an inadvertent exchange of power between interconnected control areas. The intrinsic characteristic of conventional generators that make them applicable for frequency control is "inertia". The large moment of inertia present in synchronous generators allows conventional/traditional power systems to maintain stability during disturbance as a result of frequency regulation action. In the past, the impacts of renewable energy plants (REPs) in large power systems were considered negligible because of low-level penetration. With the increasing penetration of REPs, most countries have developed renewable energy grid codes that REPs seeking connection to the grid must comply with [2][3][4]. In a similar manner to codes for conventional plants, REPs grid codes specify that some ancillary services must be provided during transient conditions [5]. One of these ancillary services is frequency control. Depending on the amplitude and duration of the frequency deviation, frequency control can be classified into primary, secondary, tertiary, and emergency frequency control [6]. This study is concerned with the impacts of inertia on primary and secondary frequency control because the efficiency of these controls can mitigate against the need for a higher level of frequency control. Primary frequency control is able to attenuate small frequency deviation but is not able to restore the system to the nominal frequency; thus, secondary frequency control is activated by harnessing the available power reserve to restore the system to nominal operating frequency. Secondary frequency control is also known as automatic generation control (AGC) or load frequency control (LFC). This is now becoming a complex control problem in power system design, operation, and control because of the emerging complexities such as renewable energy penetration, deregulation, new standards, environmental constraints, and various uncertainties [7]. Early studies on LFC problems focused on a single-area power systems with conventional power plant generation [8][9][10]; LFC models for hydropower plants were introduced in [11][12][13]. By interconnecting power systems, reliability and resilience are improved because there is inter-area support during deviating conditions. High voltage AC (HVAC) tie lines are used for interconnection in traditional power systems; the effects of HVAC tie line models were investigated in [14,15]. Studies for a two-area power system considering conventional multi-source plants were presented in [16][17][18][19][20] and higher-order interconnected areas were studied in [21][22][23][24][25][26]. The application of HVAC lines for long-distance transmission is associated with several drawbacks such as effect of capacitance for underground transmission [27], line inductance [28], Ferranti effect [29], asynchronous interconnection [30], and reduce active power transmission-to-length ratio [31]. Modern power systems have adopted the application of high voltage DC (HVDC) transmission link for long distance transmission to provide solutions to inherent HVAC transmission problems. This leads to economic and technical advantages [30,32,33].
In the traditional power system environment, the vertically integrated utility (VIU) solely owns the generation companies (GENCOs), transmission companies (TRANSCOs), and distribution companies (DISCOs), which supply power to clients at regulated prices [34,35]. The physical boundaries of a VIU define the control areas which are interconnected via tie lines. In a restructured or deregulated environment, GENCOs may choose to participate in frequency regulation, and independent power producers (IPPs) or REPs seeking connection to the grid must meet the grid code requirements by providing frequency support to the system during disturbances. DISCOs have the freedom to individually contract power from GENCOs and IPPs in different areas. The transactions among the GENCOs, IPPs, and DISCOs are supervised by an independent system operator (ISO) or any other transparent organization [6]. Provision of ancillary services such as frequency control is optional and usually based on competitive pricing, thus it is regarded as an open energy market [36]. Frequency response models for power systems in a deregulated environment must consider the contracts between DISCOs and GENCOs; Sekhar et al. [37] investigated LFC models for hydro plants. Thermal plants were considered in [38] and multi-sources were studied in [39].
In modern power systems, the impacts of REPs can no longer be considered negligible because of the increased level of penetration. The increasing integration of renewable energy plants (almost zero-inertia) reduces the overall inertia of the system, which makes the system more susceptible to high transient frequency deviations, thus making LFC problems more complex [40]. The reduction in inertia causes some problems in the deregulated environment for the ISO, such as a reduction in frequency and system oscillation. Due to the rotational parts present in wind energy conversion systems, recent grid codes require that they participate in frequency regulation within a specified range [41]. Anaya-Lara et al. [42] investigated the capability of doubly-fed induction generators for providing short term frequency regulation through rotor flux manipulation. The kinetic energy in a WECS turbine and dc-link capacitor energy were extracted to provide extra inertia in [43]. A comparative evaluation of short term frequency response between fixed-speed and doubly-fed induction generators (FSIG and DFIG) was done in [44]. The results show that the DFIG can provide more kinetic energy than the FSIG. Several studies have proposed primary and secondary frequency control strategies for grid-connected WECS from the generator point of view; however, they failed to address the limitations of the power electronics devices that connect the generator to the grid [45]. For the power electronics-based devices, the emulation of inertia for frequency response can be achieved by the control of power that is proportional to the first-order derivative of the system frequency [46,47]. This control scheme is termed "virtual inertia control" because it imitates the inertia characteristics of a conventional synchronous generator. The concept of a virtual synchronous generator (VSG) or a virtual synchronous machine (VSM) was originally proposed by Project VSYNC [48]. The VSG consists of an energy storage device interfaced with a power electronics device (mostly VSCs) and a control algorithm to make it behave as a synchronous generator for intermediate energy requirements. This mimics the damping and inertia characteristics of a synchronous generator by injecting active power when the frequency deviation is negative and vice versa. In this work, a virtual inertia control strategy is proposed. A grid-connected PMSG is considered as an intermediate frequency controller for short-term inertia response, which is highlighted in the red box in Figure 1, given that the time constant of the power electronics devices is smaller when compared to conventional devices. However, for long-term frequency regulation, the use of ESS can be adopted. In frequency control studies, the control strategies which are applied in traditional power systems cannot be extended to modern power systems with renewable energy penetration and power system deregulation. This makes it necessary to search for new and efficient control strategies [49]. Early designs of load frequency controllers used classical PI, PID, robust and optimal control theory, a dual-mode PI control [50], and decentralized PI control [51]. These were proposed for LFC studies. The hydro-powered plant was investigated in [52]; the PI control strategy was used to design the LFC scheme. By analyzing the small-signal stability of a hybrid power system-diesel engine generator, wind turbine, PV, and ESSs in [53]-a PI-based scheme was proposed to mitigate frequency deviations during load and RES variations. In [54,55], optimal linear control theory was used to model linear regulators for AGC studies. To improve the robustness of the frequency control, the authors of [56,57] used linear matrix inequalities. The Ricatti equation was used in [58] and the authors of [39,[59][60][61] adopted optimal feedback theory in traditional and deregulated power system environments. The drawbacks of these approaches such as non-robustness to parametric variations, fault intolerance, and inaccurate estimation of state variables have led to the research in intelligent control approaches [62]. Recently, modern intelligent control methods such as artificial neural networks (ANNs), fuzzy logic, evolutionary algorithms (genetic algorithms (GAs) and differential evolution (DE) algorithms), and heuristic algorithms (ant colony optimization (ACO) algorithms, particle swarm algorithms (PSOs), and artificial bee colony (ABC) algorithms) have been utilized as efficient tools for solving complex and nonlinear engineering problems. In [63], a combination of fuzzy logic, genetic algorithm, and neural network was used to design a load frequency control for a single area power system. The flatness-based approach was applied to a three-area power system with a high penetration of wind energy to reduce frequency deviations [64]. With the penetration of electric vehicles, a PSO-based ANN technique [65] and model-free control, based on the sliding mode technique [66], were used for the secondary frequency control. An adaptive-VSG was used in [67] to improve the frequency response of a permanent magnet synchronous generator-based WECS. AGC of a DFIG-based WECS using the least square support vector machine approach was proposed in [68]. Donde et al. [35] used the gradient newton algorithm (GNA) for trajectory sensitivity analysis to determine the optimal parameters of the system for LFC studies. A differential evolution algorithm was applied to a fuzzy-PID controller to optimize the control parameter for ACE signals in a deregulated environment [69].
From the literature survey, the following omissions were observed: no study considered the combination of integration of renewable energy and deregulation; most reported research studies have focused on the demerits and limitations of increasing renewable energy penetration in terms of intermittent operation and variability [70]. However, since deregulation and grid codes require IPPs to participate in electricity pricing [6] and frequency support [71], respectively, this study investigated how renewable energy plants can perform frequency support operation in a deregulated environment given that modern wind energy conversion systems (WECS) are equipped with synchronous generators and power electronic converters that can produce virtual inertia control in a restructured power system environment.
Because the power system industry is conservative and reluctant to replace well-known classical controllers with modern control technologies, heuristic optimization techniques have been used to optimally tune the controller in order to adapt to the revolving power system environment. The artificial bee colony (ABC) algorithm is a relatively new member of the heuristic optimization techniques. It is based on foraging behavior, learning, memorizing, and information sharing characteristics of honeybees. Several research papers have already evaluated and established the performance of the ABC algorithm against other algorithms such as PSO, ACO, and other optimization techniques [72][73][74][75]. Based on these investigations, it was established that the ABC algorithm has less computational complexity, high solution accuracy, simplicity, and independent convergence. The authors of [76][77][78] used the ABC algorithm in the maximum power point tracking (MPPT) control of a PV farm. The proposed controller provides a better tracking reference than the PSO-based MPPT. For optimal allocation and sizing of DGs in a distributed system, the ABC algorithm was employed in [79]. Due to its advantages, the ABC algorithm has found application in this work. To the best of the authors' knowledge, the ABC optimization algorithm has not been applied to frequency control problems in a deregulated environment. The major contributions of this paper are: development of a small-signal frequency response model for HVAC/DC interconnected power systems with high integration of renewable energy in a restructured power system environment; development of a short-term virtual inertia control strategy for wind energy conversion systems in the restructured power system environment; and novel application of the ABC algorithm to determine the optimal gain value of the virtual inertia controller.
The rest of this paper is organized as follows. Section 2 presents a detailed explanation the small-signal model used for this study in a restructured environment. The proposed derivative control strategy for grid-connected WECS is presented in Section 3. The ABC optimization algorithm is explained in Section 4. In Section 5, the validation and discussion the results of the simulations are outlined.The conclusions are presented in Section 6.

Small-Signal Model of Interconnected Power System in a Restructured Environment
The modern power system is a multi-area, highly complex, nonlinear, and time-variant environment because there exists various dynamics with varying response times that need complex numerical computations to simultaneously handle various time steps. For the frequency response study, a linearized first-order model is sufficient to analyze the dynamic behavior of the system. The model used in this study is presented in Figure 2; it represents a two-area power system that is interconnected by HVAC and HVDC lines (modified Kundur model) [80,81]. Each area has two conventional powerplants and a wind farm; there are industrial and domestic loads in both areas. Each area has to operate at the same frequency within the area. In the frequency response, primary frequency control is provided locally by each generating plant (pink label), secondary frequency control is achieved by each area (blue label), and the wind farms are equipped with the proposed virtual inertia control for frequency support (red label). The power-frequency response of the system has to be demonstrated as being suitable for modeling the system in the frequency domain (Laplace domain). In this section, the small-signal model of the two-area power system in a deregulated environment is developed, as shown in Figure 3 [82]. There are four GENCOs, two IPPs, and four DISCOs in the model. The DISCOs in Areas 1 and 2 have an aggregated load of PL1 and PL2, respectively, and all GENCOs participate in secondary frequency control.
In the open energy market, DISCOs can purchase power from GENCOs and IPPs at competitive prices. DISCOs have the autonomy to choose which GENCOs to negotiate prices with; i.e., they may or may not have contracts with GENCOs in the same area. To allow for different types of contracts such as Poolco (Power pool company), bilateral, and combinations of both, the concept of the contract participation matrix (CPM) is introduced [83]. The CPM is a matrix with m rows and n columns, where m corresponds to the number of GENCOs and n corresponds to the number of DISCOs in the system. The ijth element of the CPM is a fraction of the total load contracted by the jth DISCO to the ith GENCO. The sum of all the elements in a column of the CPM must be equal to 1 as given by [34]: In this study, there are four GENCOs and four DISCOs that participate in the deregulation contract. The two wind farms participate in the frequency response control to reduce the power deviation caused by the load demand variations from the DISCOs. Hence, Since the load demand variations from the DISCOs must be reflected in the dynamics of the GENCO model, the GENCOs must follow the information signals sent by the DISCOs in contract with them. Figure 4 shows the typical schematic representation of a four GENCO-DISCO contract.  The change in frequency as result of mismatch between the generated power and load demanded in the interconnected system can be defined by: where H i and D i are the inertia and damping constant in area i, respectively; ∆Pmi is the change in mechanical power of GENCO i; ∆Pwi is the change in wind power output; ∆P Tie12 is the tie line power deviation between Areas 1 and 2; and ∆P Li is the lumped demand in area i.
To formulate the frequency response in a deregulated environment, the absent information from the traditional environment is introduced [84]. Firstly, the area control error (ACE) signals sent to the speed governor of the GENCOs are modified to include the scheduled demand from the DISCOs. For the model considered in this paper, the equations for the governor dynamics can be expressed as: and The scheduled steady-state power flow in the tie line can be defined as ∆P tiesch = (load demand variation of DISCOs in Area 1 to GENCOs in Area 2) − (load demand variation of DISCOs in Area 2 to GENCOs in Area 1). It can be mathematically expressed using The net tie line power flow can be expressed as: Where ∆P tieact is the total power flow of the ac tie line and the HVDC link which is Following on from this: The impact of adding an HVDC link in parallel with the HVAC line is shown in the frequency deviation plot in Figure 5. It can be observed that the parallel AC/DC system suppresses the oscillations associated with the AC system and settles more quickly with zero steady-state error.
To generate the ACE signal by the secondary controller for each area of the interconnected systems, the net tie line power deviation must be included, which can be by deduced using: In the steady-state, the power generated by the GENCOs must match the demand of DISCOs in contract with them. This can be expressed as: The load demand variation for the DISCOs in each area is given by:

Derivative Control Strategy in Grid-connected Wind Energy Systems for Frequency Support
The main configuration of a grid-connected WECS is depicted in Figure 6. It consists of a wind turbine directly coupled to the generator (PMSG); the generator is connected to the grid through the bidirectional converter topology and a transformer. Traditionally, WECS are not considered in frequency regulation and therefore do not to alter their power production during frequency deviations. More recently, due to the increasing penetration of the full power converter-based wind energy systems to the grid, most wind farms are expected to respond to frequency changes in the system by changing their active power set-points. For example, the South Africa renewable energy grid codes specify that grid-connected wind farms of category C (20 MVA and above) are required to regulate active power at 50 Hz and respond to frequency changes by regulating active power generation within its limits [85]. An active power reserve margin known as P delta , usually 3% above the available wind power, is used for primary control. For grid-connected wind energy systems to suppress the frequency perturbations during load demand variations, the machine side converters (MSCs) should be able to imitate the characteristics of the prime mover used in conventional synchronous generators [86]. To achieve this aim, an extra control loop is added to the MSC control scheme, as shown in Figure 7. The modified control scheme of the MSC includes the conventional active and reactive power control loops, and a virtual inertia control loop. The control scheme has d-q voltage controllers, in which the output powers (active and reactive) of the PMSG are controlled using the d and q components of the voltage references. The voltage references are obtained from the PI controllers used to regulate the d and q reference currents. The grid side converter (GSC) control scheme is usually used to regulate the dc-link voltage by considering power balance between the grid and the generator; the GSC control architecture is described in reference [87]. The additional control loop in the MSC control scheme forms the basis for virtual inertia power that can contribute to the increase of the overall inertia response of the system by injecting/absorbing active power during frequency dips/spikes. This work proposes a derivative control strategy to implement the control loop for virtual inertia emulation. The proposed control strategy is capable of reducing the transient frequency deviation by injecting short-term active power into the system. In this work, the wind farms in both areas are equipped with the proposed control scheme to mitigate frequency deviation in both areas. It is assumed that the wind farms are operating at a de-loading point, i.e. below the maximum power that can be extracted at the current wind speed.
To represent the proposed derivative control strategy in the frequency domain, the control law can be evaluated from [47]: where T wi is the time constant for imitating the dynamic response of the MSC in area i and K wi is the derivative virtual inertia control gain in area i. From Equation (20), it is clear that the gain of the controller will modify the active power set point of the MSC based on the rate of change of the frequency and this can be illustrated by the control diagram in Figure 8. To effectively support the system during disturbances, the control gain parameter must be tuned properly. As stated above, the ABC optimization algorithm is used to tune the control gain. To ascertain the robustness of this optimized control strategy, a comparative evaluation with the conventional tuning method (interior-point algorithm), as presented in [46], is designed.
Using the generalized state-space equation, Equations (3)- (20) can be written as: where ∆x is the state vector, A is the state or plant matrix, B is the input matrix, ∆u is the input vector, ∆y is the output vector, and C is the output matrix. The elements of the state, input, and output vectors are given by

Artificial Bee Colony Optimization Algorithm
The ABC optimization algorithm is a nature-inspired optimization approach proposed by Karaboga in 2005 [88]. It is a computational intelligence method that mathematically mimics the behavior of bees in search of food. The bee colony consists of three classes of bees. The first class are known as the employed bees (B em ); they randomly explore the search area for nectar positions (possible solutions). After finding a nectar position (NP), they memorize the details of this position (nectar amount) and share the information of the NP with the colony by dancing to the other bees in the hive. The duration of the dance denotes the quality of nectar (fitness value). The second class of bees are known as the onlooker bees (B on ); this set of bees watch the dance of the employed bees before choosing an NP. A rich NP attracts more onlooker bees than a poor NP. The third class of bees are called the scout bees; these are employed bees whose NP are abandoned because of poor quality after a certain number of trials. Onlooker bees can become employed bees if they discover a new NP. In this search process, exploration and exploitation are simultaneous processes [89].
In the ABC optimization algorithm, the number of employed bees is equal to the number of onlooker bees. The number of possible NPs is equal to the number of employed bees. The initialization of the scouting process is done by scout bees and each iteration involves three major steps: (1) search for various NPs and information on the quality of each NPs; (2) selection of NP by the onlooker bees through the information shared by the employed bees; and (3) employed bees with poor NPs become scout bees and are then sent to new NPs.
For the initialization stage, a random distribution of initial population of solutions x i (i = 1, 2, . . . , B em ) is generated, where i is the size of the population and B em is the number of employed bees. For each solution, dimension D n is the number of parameters to be optimized. After initialization, the population of the solutions is subjected to repeated cycles (C = 1, 2, . . . , MCN) of the search process for the three classes of bees, where MCN is the maximum number of cycles of the search process. For each cycle, the employed bees modify the NP based on the local information (visible content) and the amount of nectar. If the amount of nectar in the new position is greater than the previous one, the bee memorizes it and discards the previous solution; otherwise, it retains the previous position. When all the employed bees have finished the exploration process, they share the information with the onlooker bees in the hive. The onlooker bees select an NP by evaluating the nectar information shared. The probability of selecting an NP is related to the nectar amount. For this study, the roulette wheel selection method is adopted. Thus, the probability of selecting an NP can be evaluated using: where fitness i is the fitness value of the solution i. It can be deduced from Equation (18) that a rich NP will attract more onlooker bees than a poor NP. Before the onlooker bees select another NP, they compare the fitness value of position i to i + 1. This cycle is repeated until all onlooker bees are dispersed. If the fitness of the solution does not improve for a predetermined limit, the employed bees abandon this solution and become scout bees. When a new position is selected, the cycle starts again until the final conditions are met. To determine the surrounding NP relative to the current NP, the ABC algorithm uses: where k = i, j ∈ (1, 2, . . . , D n ) and k ∈ (1, 2, . . . , B em ). For each cycle, the scout bees produce a new solution, which is given by In summary, the ABC algorithm has three control parameters: the size of the colony, the limit value, and the maximum cycle number [90]. To apply the ABC optimization algorithm to the frequency response problem of the model described in Section 2, an objective function needs to be defined. The objective of this study is to minimize the frequency deviations in both areas. The integral of the squared error (ISE) function is used to define our objective function J, which is obtained from After the definition of the objective function, the following sequence must be followed to apply the ABC algorithm to the model. The flowchart for the algorithm is shown in Figure 9.
1. Generate initial population of solution x i (i = 1, 2, . . . , B em ) from Equation (19) 2. Evaluate the fitness of the population using fitness i = 1 1 + j 2 3. Using Equation (20), generate new NPs for the employed bees and evaluate the fitness 4. Apply the roulette wheel selection process to choose an NP 5. For the onlooker bees, calculate the probability P i for the solutions 6. If all onlooker bees are dispersed, go to Step 9; else, go to next step 7. Generate new NP for the onlooker bees and evaluate their fitness 8. Apply the roulette wheel selection process 9. If there is an abandoned solution for the scout bees, replace it with a new solution and evaluate its fitness 10. Memorize the best solution reached so far 11. Repeat Steps 1-10 for another cycle until C = MCN The initialization parameters of the ABC algorithm are provided in Table 1.

Results and Discussion
Rwo analyses were performed on the system model and the results are presented with discussion in this section. The first analysis is the eigenvalue analysis, which was used evaluate and validate the stability of the system model and proposed control strategy. The second analysis is the load-frequency study, which was used to analyze the frequency response of the system in order to demonstrate the robustness of the proposed strategy with respect to different realistic market transaction scenarios. The performance was compared without wind farm support and with the wind farm providing conventional frequency support (classical interior-point algorithm), as described in [46]. The simulations were done in the MATLAB/Simulink environment. The system was modeled in Simulink and the ABC algorithm was written as a script file in the MATLAB environment. The parameters used for the simulation are documented in Table A1 of Appendix A.

Eigenvalue Analysis
The state-space model described in Section 2 was used to analyze the stability of the study system by considering the location of the eigenvalues. To obtain the eigenvalues of the system, the characteristic equation is defined from the state matrix A as: where λ is the eigenvalues that can be real or complex and I is an identity matrix. The elements of state matrix A are given in Appendix B. The locations of λ in the s-plane determines the stability of the system. From Lyapunov's first method of stability [91], the following stability criteria apply: 1. If all of the eigenvalues have a negative real part, the system is asymptotically stable. 2. If one or more of the eigenvalues have a positive real part, the system is unstable.
The complex conjugates of eigenvalues can be defined by: where σ is the damping and ω is the frequency of oscillation. A positive σ corresponds to an increasing oscillation while a negative σ corresponds to a damped oscillation. The eigenvalues of the system used for this study are shown in Table 2. Table 3 shows the total amount of system oscillation damping. It was deduced that the addition of the wind farm adds two extra modes to the system. The total damping increased and a significant improvement was noticed using the proposed optimization method. In Figure 10, the positive effect of the proposed strategy is observed by the shifting of the eigenvalues to the left side of the s-plane.  From the presented results, it can be concluded that the participation of renewable energy sources in frequency control can improve the system performance through proper control. The improvement is largely reflected when the control parameters are obtained using the ABC algorithm through the oscillation damping.

Load-Frequency Analysis
To investigate the robustness of the proposed control strategy in this subsection, the response of the system in Figure 3 presents the simulation results of the proposed control strategy. For the initial execution of the algorithm, a wide search space is allotted to observe the convergence of the solution. This is reduced afterwards. In a restructured environment, the control system should be efficient for all possible power system transactions. Three transaction scenarios were considered in this study.

Scenario 1: Poolco Transaction
In a Poolco based contraction, the GENCOs participate in the frequency response control of their own control area only. In this scenario, the load demand variations of DISCO1 and DISCO2 are provided by GENCO1 and GENCO2, while the load demand variations of DISCO3 and DISCO4 are provided by GENCO3 and GENCO4. For a total load demand variation of 0.2 pu in each area, DISCO1 and DISCO2 demands are equal from GENCO1 and GENCO2, while DISCO3 and DISCO4 demands are equal from GENCO3 and GENCO4. The CPM for a Poolco based contract between the GENCOs and DISCOs used in this study is given by:   Figure 11 represents the plot of frequency deviations in both area. It can be observed that the amplitudes of the first undershoot and overshoot are highest when there is no frequency support provided by the wind farm. With the wind farm providing support using a conventional strategy, the undershoots and overshoots are reduced by 15% and 24%, respectively, for Area 1. Using the proposed control strategy, the deviations are quickly reduced to zero steady-state with minimal oscillations in both areas. In the steady-state, the power generated by the GENCOs must match the demand of the DISCOs in contract with them. From Equation (17), the power generated by each GENCO for the respective DISCOs in contract with it can be calculated from Figure 12 shows the various changes in active power generated by the GENCOs due to the load demand variations. The new steady-state value of each GENCO is the amount of contracted demand from the DISCOs in its area. It can be observed that the participation of wind farms in frequency control can reduce the amount of transient power needed by the GENCOs to reach the new steady-state values. By using the ABC algorithm for optimal control, the wind farm is capable of providing active power, thereby reducing the amplitude of the transient power provided by the GENCOs. It can be observed that overshoots were reduced by 16%, 14.3%, 16%, and 13.3% in the GENCOs when compared to the conventional control strategy. Since there is no exchange of power between the areas, the scheduled tie line power flow will be zero as defined by Equation (10), so that ∆P tiesch = 0 pu However, tie line perturbations will occur when there is load demand variation due to interconnection, as shown in Figure 13. It is generally observed that the proposed optimized controller helps the system to reach a faster settling time, with the system experiencing damping, which suppresses the oscillation, thereby providing a better dynamic response.

Scenario 2: Bilateral Transaction
In bilateral based transactions, DISCOs can be in contract with any GENCO within or outside their own areas; i.e., GENCOs in Area 1 can transfer power to DISCOs in Area 2 and vice versa. Therefore, inter-area exchange of power via the tie lines can be scheduled at a non-zero steady-state value. In this study, each DISCO demands 0.1 pu power from the GENCOs in both areas, as defined by the CPM below: The non-diagonal elements of the matrix correspond to the contract of a DISCO in one area with a GENCO in another area. Figures 14-16 show the response of the system as a result of bilateral transaction. Figure 14 shows the frequency deviation in Areas 1 and 2. It can be seen that the peak undershoot and overshoot are high when there is no support from the wind farms. When the wind farms participate in frequency control using the conventional control, the undershoot and overshoot are reduced by 15% and 18% in Area 1 and 13% and 50% in Area 2, respectively. A further reduction is noticed when the gain of the controller is optimized using the ABC algorithm. The oscillations are suppressed and the system frequency deviation settles faster to zero in the steady-state condition. From Equation (16), the change power generated by each GENCO in the bilateral transaction scenario can be calculated from The change in power generated by each GENCO is presented in Figure 15. It can be observed that, without the wind farm support, there is increased transient power generation by the GENCOs;, this deviation is reduced when there is support from the wind farm by providing active power for a short time. The effect of wind farm participation is largely noticed when the optimized derivative control is applied. This results in over 20% reduction in transient power from the GENCOs when the proposed strategy is implemented.
The scheduled tie line power flow is shown in Figure 16. This illustrates the exchange of power between both areas at steady-state. This can be calculated using

Scenario 3: Contract Violation
At times, DISCOs may demand more power than that specified in the contract. This extra variation is usually not contracted to any GENCO. When this happens, only GENCOs in the affected area can respond to the excess demand variations from the DISCO. The excess demand variation is always regarded as a local load demand variation. Modifying Scenario 2 (bilateral transaction) with an extra demand variation of 0.1 pu from DISCOs in Area 2 at 15 s, the CPM in Scenario 2 remains the same.
Beginning the analysis from 15 s in Figures 17-19, the effect of the un-contracted load demand variation in Area 2 can be observed. From the frequency deviation plots shown in Figure 17, it can be seen that the extra load demand variation in Area 2 causes frequency dips in Areas 1 and 2 due to interconnection. The oscillations associated with the system without wind farm support are damped when the wind farms participate in frequency control; however, the proposed optimized controller has a better damping response than the conventional control strategy.
In this transaction, the change in power generation due to un-contracted variation is the responsibility of the GENCOs in that area; that is, the GENCOs in other areas return to the initial steady-state values after participating in the frequency control. In Figure 18, the change in power generated by GENCO1 and GENCO2 remains the same; i.e., ∆P m1 = 0.14 pu and ∆P m2 = 0.08 pu. However, there are new steady-state values for GENCO3 and GENCO4. The changes after 15 s can be calculated using:

Scenario 4: Wind Power Fluctuation
Since wind power is largely dependent on wind speed, which is usually intermittent, it is expected that the contribution of wind farms in frequency response control will be limited. In this work, it is assumed that the de-loading margin is available in all range of wind speeds and the short-term support lasts for a fraction of second to few seconds. To test the robustness of the proposed control strategy, during wind power fluctuation, as shown in Figure 20, a change in load demand of 0.15 pu in Area 1 is simulated at 5 s. In Figure 21, the frequency deviations in both areas of the system is presented for cases where the wind farms support in frequency stability using the conventional tuning method (dashed blue) and ABC tuning method (red). It is observed that the proposed tuning method provides better damping characteristics and reduced frequency excursions. It is pertinent to point out that the control approach proposed in this work can be extended to energy storage systems to provide mid-and long-term frequency support to power systems.

Conclusions
In this study, the artificial bee colony algorithm was used in the implementation of a virtual inertia control strategy. This involves the formulation, development, and simulation of a system model that can be used in traditional power systems and can be extended to modern power systems with renewable energy penetration and power system deregulation. This optimization algorithm was used to investigate the positive impact that renewable energy has on the network with respect to frequency control. The heuristic method is used to optimally tune the gain of the virtual inertia control scheme that is proposed to reduce the frequency transient deviations in a restructured power system environment. The control strategy was tested in a two-area system having conventional generation, wind power plants, and parallel ac/dc links to model a simple modern power system. Three transaction scenarios (Poolco, bilateral, and contract violation) were implemented to show the robustness of the proposed control strategy. It is generally observed that peak overshoots and undershoots, settling time, and oscillations of the system were greatly reduced. In addition, the overall dynamic performance and stability of the system are improved. The proposed strategy can be applied to expanded power systems that include more interconnected areas with additional generating plants such as solar, geothermal, ocean, and tidal power plants.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The