Impact of Renewable Energy Sources into Multi Area Multi-Source Load Frequency Control of Interrelated Power System

Abstract

There is an increasing concentration in the influences of nonconventional power sources on power system process and management, as the application of these sources upsurges worldwide. Renewable energy technologies are one of the best technologies for generating electrical power with zero fuel cost, a clean environment, and are available almost throughout the year. Some of the widespread renewable energy sources are tidal energy, geothermal energy, wind energy, and solar energy. Among many renewable energy sources, wind and solar energy sources are more popular because they are easy to install and operate. Due to their high flexibility, wind and solar power generation units are easily integrated with conventional power generation systems. Traditional generating units primarily use synchronous generators that enable them to ensure the process during significant transient errors. If massive wind generation is faltered due to error, it may harm the power system’s operation and lead to the load frequency control issue. This work proposes binary moth flame optimizer (MFO) variants to mitigate the frequency constraint issue. Two different binary variants are implemented for improving the performance of MFO for discrete optimization problems. The proposed model was evaluated and compared with existing algorithms in terms of standard testing benchmarks and showed improved results in terms of average and standard deviation.

1. Introduction

Renewable energy sources (RESs) and especially wind technology are treated to be the most effective viable technology due to their environmental blessings, and their value of procedure and servicing have declined considerably within a few years. Hybrid plants secure the stability of supply commixture various sustainable energy resources like photovoltaic, wind generators, and even diesel generator sets (DGs) used for back-up purposes [1]. Therefore, grid connection of those together with traditional plants is adopted because of enhanced behavior with respect to effective load. It is ascertained that variation in frequency is induced because fluctuation in load is low, strengthening the insertion of inexhaustible resources. Load frequency control (LFC), in addition to the Proporional Integral Derivative (PID) controller [2,3] is suggested to overcome frequency inconsistency for a power system involving wind, hydro, and thermal units due to load and generating power variation induced due to the insertion of inexhaustible resources [4,5]. A system comprising thermal plants, hydro plants, and wind power plants will be designed with the help of MATLAB [6].

Our contributions in this paper are as follows: First, we recommend the two alternatives of binary moth flame optimizers to unravel the frequency restriction matter. We put into practice two diverse binary alternates for civilizing moth flame optimizer (MFO) behavior for distinct optimization tribulations [7]. In the primary variant, i.e., binary moth flame optimizer (BMFO1), coin flipping-based assortment probability of binary statistics is applicable. We applied the superior sigmoid transformation in the subsequent variant known as BMFO2. Along with Harris hawks optimizer (HHO) algorithms, these advanced algorithms are veteran and examined for a variety of unilateral, bilateral, and contract violation optimization problems. Secondly, Section 2 explores the impact of renewable energy sources on the load frequency control problem. Section 3 depicts the load frequency control issue’s mathematical behavior when integrated with a renewable energy source. Section 4 shows the transfer function model of the multiarea multisource power system integrated with a renewable energy source. Lastly, in Section 5, all these algorithms are estimated and evaluated in conditions of typical testing benchmarks in which the projected HHO model has superior consequences with regards to mean and standard deviation. Finally, Section 6 winds up the paper.

2. Impact of Renewable Energy Sources

Renewable energy sources (RESs) definitely disturb the vibrant performance of the power system in such a manner that may be diverse from predictable generating units. Traditional generating units primarily use synchronous generators that enable them to ensure the process during major transient errors. If, due to error, a massive amount of wind generation is faltered, then the adverse influence of that error on the power system’s operation with the LFC issue could also be expanded. High penetration of renewable energy in control arrangement may raise some reservations during the irregular procedure. It familiarizes numerous technical implications and exposes significant questions regarding what happens to LFC requirement [8,9] after the addition of various RESs to the present production system and whether the outdated mechanism methodologies are sufficient to operate in a fresh situation.

The influence on optimal flow of power, voltage, and management of frequency, power quality, and structure economics are increased due to the addition of RESs into control system grids. According to the behavior of RES power deviation, the impression on the LFC concern has involved rising research attention throughout the last era [10,11,12]. Substantial interrelated frequency fluctuations can cause over- or underfrequency transmitting and remove certain generations and loads. Under opposed situations, this may affect a dropping disappointment and miscarriage of the system. Figure 1 shows the block diagram of a wind energy generation system [13].

2.1. Advanced Optimization Methods

Optimization plays a major role in many fields of engineering. It is a path in which a suitable reaction to an exceptional matter is exposed via a search tool [14]. By means of upgradation in expertise, the novel creation of the matter fortitude optimization approach identified as metaheuristic has been used to deliberate mathematical humanity. Metaheuristic algorithms (MA) impersonate a systematic way to obtain the finest consequences for adilemma. MA act as a fantasy search for a good stipulation in an optimization concern [15]. In this paper, three advanced optimization techniques, BMFO1, BMFO2, and HHO, are considered and compared with conventional methods.

2.1.1. Binary Moth Flame Optimizer (BMFO1)

The basic MFO is a nature-inspired heuristic search technique that imitates moths’ navigational possessions about artificial lights. BMFO1 is a recently predictable metaheuristics search algorithm proposed by Mirjalili [16], which is re-energized by navigation actions of moth and their meeting close to the beam. It helps to recover the exploitation search of the moths and diminish the quantity of flames. Even if moths have a tough potential to sustain a protected loom with respect to the moon and grip a bearable assembly for traveling in a traditional scratch for broad remoteness. They are also caring in a serious/inoperative bowed pathway over a replicated source of illumination.

2.1.2. Modified SIGMOID Transformation (BMFO2)

The binary calibration of stable quest accommodation and spaces of investigating council, resolution to binary searching domicile could be obligatory for optimizing binary ecological issues such as LFC. A modified sigmoidal transfer function is assumed in the projected work, which has better presentation than any more substitute of it as shown in [17]. Basic moth flame optimizer applied with the modified sigmoidal transformation (BMFO2) is used to carry out the binary chart of actual moth value and flame location for fixing the LFC problem.

2.1.3. Harris Hawks Optimizer

HHO [17,18] is a gradient-free and population-centered algorithm containing unfair and investigative steps for wonder swoop, the fauna of prey assessment, and assorted ploy built on brutal speculation of Harris hawks.

3. LFC Model with Integration of Renewable Energy Source

When renewable energy sources are integrated into the power arrangement, a supplementary cause of the deviation is added to the arrangement’s adjustable behavior. To examine the deviations triggered by renewable energy source plants, the entire conclusion is significant, and every modification in renewable energy source output power does not requisite to be accorded with the conversion in other generating units running in the opposite way [19,20]. Sudden variations in load and renewable energy source output power might intensify each other either completely or partially. However, the sluggish renewable energy source power instability and entire average power fluctuation adversely subsidize power inequity and frequency distribution, which could be considered for the LFC arrangement. This energy variation must comprise a predictable LFC arrangement [21].

A general LFC model integrated with renewable energy source [22,23] is shown in Figure 2 in which the conforming blocks for governor dead-band, Generation Rate Constraints (GRC), and time delays are not involved. In this model, altered parametric standards are used for generator regulation and for turbine-governor to shelter the diversity of production categories in the control area. The expressed components and blocks are defined as follows: ΔP_m is mechanical power, Δf is frequency deviation, ΔP_L is load disturbance, ΔP_C is supplementary frequency control action, D_Sys is equivalent damping coefficient, H_Sys is equivalent inertia constant, βis frequency bias, α_i is participation factor, R_iis drooping characteristic, ΔP_P is primary frequency control action, ΔP_RESis renewable energy source power fluctuation, M_i(s) is a governor-turbine model, ACEis area control error [24], and finally, ΔP′_L and ΔP′_tie are amplified local load alteration and tie-line power vacillation signals, respectively.

In the revised LFC structure, the efficient ACE signal must signify the influences of renewable energy on planned stream over tie-line and local power variation through area frequency. The ACE signal is conventionally described as a linear amalgamation of tie-line power and frequency fluctuations as follows [25,26]:

$$\mathrm{ACE}=\mathrm{B}\Delta \mathrm{f}+\Delta {\mathrm{P}}_{\mathrm{tie}}$$

(1)

In the conventional power system, ΔP_tie(ΔP_tie−C) is a deviation between actual and scheduled energy flow over the tie-lines.

$$\Delta {\mathrm{P}}_{\mathrm{tie}-\mathrm{C}}={\displaystyle \sum }\left({\mathrm{P}}_{\mathrm{tie},\mathrm{actual}}-{\mathrm{P}}_{\mathrm{tie},\mathrm{Scheduled}}\right)$$

(2)

The difference among the updated LFC model in Figure 2 and the conventional one provides two new signals which represent the transient behavior of RESs on tie-line power and local load variations (ΔP′_L, ΔP′_tie−RES):

$$\Delta {\mathrm{P}}_{\mathrm{L}}^{\prime }\left(\mathrm{s}\right)=\Delta {\mathrm{P}}_{\mathrm{RES}}\left(\mathrm{s}\right)-\Delta {\mathrm{P}}_{\mathrm{L}}\left(\mathrm{s}\right)$$

(3)

$$\Delta {\mathrm{P}}_{\mathrm{tie}-\mathrm{RES}}= {\displaystyle \sum }\left(\Delta {\mathrm{P}}_{\mathrm{tie}-\mathrm{RES},\mathrm{actual}}-\Delta {\mathrm{P}}_{\mathrm{tie}-\mathrm{RES},\mathrm{Scheduled}}\right)$$

(4)

After adding a substantial amount of renewable power to the conventional power flow of tie-lines (ΔP_tie−C) in the power system, the updated renewable energy source power via tie-lines (ΔP_tie−RES) must be painstaking. Consequently, the restructured tie-line power fluctuation can be articulated as follows:

$$\begin{array}{c}\Delta {\mathrm{P}}_{\mathrm{tie}}^{\prime }= \Delta {\mathrm{P}}_{\mathrm{tie}-\mathrm{C}}+ \Delta {\mathrm{P}}_{\mathrm{tie}-\mathrm{RES}}\\ ={\displaystyle \sum }\left(\Delta {\mathrm{P}}_{\mathrm{tie}-\mathrm{C},\mathrm{actual}}-\Delta {\mathrm{P}}_{\mathrm{tie}-\mathrm{C},\mathrm{Scheduled}}\right)+{\displaystyle \sum }\left(\Delta {\mathrm{P}}_{\mathrm{tie}-\mathrm{RES},\mathrm{actual}}-\Delta {\mathrm{P}}_{\mathrm{tie}-\mathrm{RES},\mathrm{Scheduled}}\right)\end{array}$$

(5)

The entire renewable energy source power flow modification is generally smoother as compared to fluctuation in influences from personal renewable energy source elements. Using Equations (1) and (5), the restructured ACE signal can be accomplished by Equation (6):

$$\begin{array}{c}\mathrm{ACE}=\mathrm{B}\Delta \mathrm{f}+\Delta {\mathrm{P}}_{\mathrm{tie}}^{\prime }\\ =\mathrm{B}\Delta \mathrm{f}+ \left({\displaystyle \sum }\left(\Delta {\mathrm{P}}_{\mathrm{tie}-\mathrm{C},\mathrm{actual}}-\Delta {\mathrm{P}}_{\mathrm{tie}-\mathrm{C},\mathrm{Scheduled}}\right)+{\displaystyle \sum }\left(\Delta {\mathrm{P}}_{\mathrm{tie}-\mathrm{RES},\mathrm{actual}}-\Delta {\mathrm{P}}_{\mathrm{tie}-\mathrm{RES},\mathrm{Scheduled}}\right)\right)\end{array}$$

(6)

where P_tie-C,actual is actual conventional tie-line, P_{tie-C,scheduled} is scheduled conventional tie-line, P_{tie-RES,actual} is actual renewable energy source tie-line, and P_{tie-RES,scheduled} is scheduled renewable energy source tie-line powers.

To justify the investigation part of the expected algorithm, various conventional as well as advanced optimization techniques such as Particle Swarm Optimization (PSO), MFO, BMFO1, BMFO2, and HHO were compared in terms of best value, worst value, mean and standard deviation. The different gain values for various algorithms for multiarea modal were also compared.

Table 1. Performance of various algorithms in the proposed system.

Parameter	No. of Trials	Mean	Standard Deviation	Best	Worst
PSO	50	1.13807	1.00846	1.02774	1.05084
MFO	50	0.0489	0.00644	0.02835	0.05073
BMFO1	50	0.0284	0.00083	0.02745	0.03216
BMFO2	50	0.0272	0.00011	0.02764	0.02811
HHO	50	7.62 × 10−93	1.39 × 10−93	1.4 × 10−112	4.17 × 10−92

presents the performance of various algorithms in the proposed system and

Table 2. Improved gain values of Proportional Integral (PI) controller for the expected scheme with various algorithms.

Controller Type	Controller Parameters
	Area1		Area2
	K1p	K1Int	K2P	K2Int
PSO	0.53007	−7.12 × 10−6	0.28516	−0.84672
MFO	0.52333	−0.70889	0.31424	−0.92724
BMFO1	0.53126	−0.0034	0.28435	9.83 × 10−6
BMFO2	0.55025	−8.07 × 10−7	0.28545	−1.89 × 10−6
HHO	0.15025	−9.059 × 10−9	0.18293	−1.91 × 10−9

depicts the Proportional Controller (PI) controller’s gain values for various algorithms.

4. Transfer Function Model Multi Area Multi-Source Hydro-Thermal System with Wind Power Plant

The transfer function model of the hydro and thermal generating unit has been derived. The tie-line combines them in two methods, namely two-area, which is conventional, and then multiarea with renewable energy source, a new concept. On analyzing, the thermal system, when exposed to unit step load disturbance of 0.01 per unit in Area1 alone, the area frequency and the tie-line energy oscillates and settles with offset. The offset can be removed by including a secondary controller, which varies the governor’s power reference setting. An optimal secondary controller is to be developed for the effective operation of the hydro–thermal power system integrated with renewable energy sources like wind energy.

5. Results and Discussion

Sensitivity analysis with respect to frequency variations, fault in tie-line, actual power flow, and output reaction of dissimilar generators after instant load fluctuation in the planned structure in terms of dispersion of wind under various contracts, such as unilateral, bilateral, and contract-violation case have been publicized in Figure 2.

5.1. Unilateral Transaction

The simulations with penetration of wind have been achieved to check the expected reaction of power scheme with respect to area frequency, power flow among interrelated areas, and the reaction of generating units through an unexpected load modification state in terms of expected agreements of the decontrolled electricity market as shown in Figure 3, Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8.

Dynamic response of Area1 frequency with respect to time in seconds with various controllers like conventional, PSO, MFO, BMFO1, BMFO2, and HHO under unilateral contract are compared in Figure 3 and

Table 3. Graph analysis of Area1 frequency with various controllers under unilateral contract.

Controller	Delay Time	Rise Time	Peak Overshoot Time	Settling Time
Conventional	1.05	1.186	1.375	22
PSO	0.8833	1.05	1.186	21.430
MFO	0.8833	1.05	1.145	21.3
BMFO1	0.9667	1.05	1.145	21.02
BMFO2	0.8833	1.05	1.087	20.1
HHO	0.8833	0.9667	1.07	19.8

. The comparative outcomes obtained show improved results of the HHO controller with delay time 0.8633 s and settling time of 19.8 s.

Dynamic response of Area2 frequency with respect to time in seconds with various controllers like conventional, PSO, MFO, BMFO1, BMFO2, and HHO under unilateral contract are compared in Figure 4 and

Table 4. Graph analysis of Area2 frequency with various controller under unilateral contract.

Controller	Delay Time	Rise Time	Peak Overshoot Time	Settling Time
Conventional	1.474	1.654	1.868	23.39
PSO	1.336	1.414	1.594	21.17
MFO	1.336	1.474	1.594	21.01
BMFO1	−1.201	−1.414	−1.534	21.43
BMFO2	1.297	1.375	1.414	21.43
HHO	1.297	1.375	1.414	20.86

. The comparative outcomes obtained show improved results of the HHO controller with delay time 1.297 s and settling time of 20.86 s.

Dynamic response of deviation in actual tie-line power flow with respect to time in seconds with various controllers like conventional, PSO, MFO, BMFO1, BMFO2, and HHO under unilateral contract are compared in Figure 5 and

Table 5. Graph analysis of deviation in actual tie-line power flow with various controllers under unilateral contract.

Controller	Delay Time	Rise Time	Peak Overshoot Time	Settling Time
Conventional	1.935	3.478	3.865	23.39
PSO	0.664	1.195	1.325	21.17
MFO	0.535	0.975	1.068	21.01
BMFO1	0.435	−0.778	−0.865	21.43
BMFO2	0.535	0.957	1.063	21.43
HHO	0.162	0.292	0.324	20.86

. The comparative outcomes obtained show improved results of the HHO controller with delay time 0.162 s and settling time of 20.86 s.

Dynamic response of deviation in tie-line error with respect to time in seconds with various controllers like conventional, PSO, MFO, BMFO1, BMFO2, and HHO under unilateral contract are compared in Figure 6 and

Table 6. Graph analysis of deviation in tie-line error with various controller under unilateral contract.

Controller	Delay Time	Rise Time	Peak Overshoot Time	Settling Time
Conventional	2.935	5.280	5.865	22
PSO	2.118	3.814	4.235	21.430
MFO	1.955	3.498	3.886	21.3
BMFO1	1.115	2.002	2.224	21.02
BMFO2	1.988	3.555	3.950	20.1
HHO	1.157	2.085	2.315	19.8

. The comparative outcomes obtained show improved results of the HHO controller with delay time 1.157 s and settling time of 19.8 s [27].

Dynamic response of Generation Companies (GENCOs) generation of Area1 with respect to time in seconds with various controllers like conventional, PSO, MFO, BMFO1, BMFO2, and HHO under unilateral contract are compared in Figure 7 and

Table 7. Graph analysis of GENCOs generation response of Area1 with various controllers under unilateral contract.

Controller	Delay Time	Rise Time	Peak Overshoot Time	Settling Time
Conventional	3.665	6.595	7.325	23.39
PSO	3.598	6.478	7.195	21.17
MFO	3.508	6.315	7.015	21.01
BMFO1	3.335	6.002	6.668	21.43
BMFO2	3.503	6.305	7.005	21.43
HHO	3.335	6.002	6.668	20.86

. The comparative outcomes obtained show improved results of the HHO controller with delay time 3.335 s and settling time of 20.86 s.

Dynamic response of GENCOs generation of Area2 with respect to time in seconds with various controllers like conventional, PSO, MFO, BMFO1, BMFO2, and HHO under unilateral contract are compared in Figure 8 and

Table 8. Graph analysis of GENCOs generation response of Area2 with various controllers under unilateral contract.

Controller	Delay Time	Rise Time	Peak Overshoot Time	Settling Time
Conventional	0.888	1.598	1.775	22
PSO	0.725	1.305	1.450	21.430
MFO	0.734	1.333	1.468	21.3
BMFO1	0.755	1.355	1.505	21.02
BMFO2	0.723	1.300	1.445	20.1
HHO	0.695	1.251	1.390	19.8

. The comparative outcomes obtained show improved results of the HHO controller with delay time 0.695 s and settling time of 19.8 s.

5.2. Bilateral Based Transaction

The various dynamic responses with penetration of wind in mutual areas and the tie line are shown in Figure 9, Figure 10, Figure 11, Figure 12, Figure 13 and Figure 14.

Dynamic response of Area1 frequency with respect to time in seconds with various controllers like conventional, PSO, MFO, BMFO1, BMFO2, and HHO under bilateral contract are compared in Figure 9 and

Table 9. Graph analysis of Area1 frequency with various controllers under bilateral contract.

Controller	Delay Time	Rise Time	Peak Overshoot Time	Settling Time
Conventional	−2.414	−4.205	−4.668	23.39
PSO	−1.508	−2.715	−3.015	21.17
MFO	−1.503	−2.705	−3.005	21.01
BMFO1	−1.003	−1.805	−2.005	21.43
BMFO2	−0.978	−1.766	−1.955	21.43
HHO	−0.835	−1.499	−1.665	20.86

. The comparative outcomes obtained show improved results of the HHO controller with delay time −0.835 s and settling time of 20.86 s.

Dynamic response of Area2 frequency with respect to time in seconds with various controllers like conventional, PSO, MFO, BMFO1, BMFO2, and HHO under Bilateral Contract are compared in Figure 10 and

Table 10. Graph analysis of Area2 frequency with various controllers under bilateral contract.

Controller	Delay Time	Rise Time	Peak Overshoot Time	Settling Time
Conventional	0.0066	0.0118	0.0131	23.39
PSO	0.0031	0.0047	0.0052	21.17
MFO	0.0045	0.0074	0.0082	21.01
BMFO1	0.0026	0.0046	0.0051	21.43
BMFO2	0.0027	0.0048	0.0053	21.43
HHO	0.0023	0.0041	0.0045	20.86

. The comparative outcomes obtained show improved results of the HHO controller with delay time = 0.0023 s and settling time of 20.86 s.

Dynamic response of deviation in actual tie-line power flow with respect to time in seconds with various controllers like Conventional, PSO, MFO, BMFO1, BMFO2, and HHO under Bilateral Contract are compared in Figure 11 and

Table 11. Graph analysis of deviation in actual tie-line power flow with various controllers under bilateral contract.

Controller	Delay Time	Rise Time	Peak Overshoot Time	Settling Time
Conventional	0.0106	0.0190	0.0211	23.39
PSO	0.0043	0.0077	0.0085	21.17
MFO	0.0063	0.0114	0.0125	21.01
BMFO1	0.0033	0.0059	0.0065	21.43
BMFO2	0.0031	0.0055	0.0061	21.43
HHO	0.0029	0.0053	0.0058	20.86

. The comparative outcomes obtained show improved results of the HHO controller with delay time 0.0029 s and settling time of 20.86 s.

Dynamic response of Deviation in tie-line error with respect to time in seconds with various controllers like Conventional, PSO, MFO, BMFO1, BMFO2, and HHO under Bilateral Contract are compared in Figure 12 and

Table 12. Graph analysis of deviation in tie-line error with various controller under bilateral contract.

Controller	Delay Time	Rise Time	Peak Overshoot Time	Settling Time
Conventional	1.284	2.311	2.568	22
PSO	1.278	2.299	2.555	21.430
MFO	1.277	2.298	2.554	21.3
BMFO1	1.277	2.297	2.553	21.02
BMFO2	1.277	2.296	2.552	20.1
HHO	1.276	2.295	2.551	19.8

. The comparative outcomes obtained show improved results of the HHO controller with delay time 1.276 s and settling time of 19.8 s.

Dynamic response of GENCOs to DISCOs generation of Area-1 with respect to time in seconds with various controllers like Conventional, PSO, MFO, BMFO1, BMFO2, and HHO under Bilateral Contract are compared in Figure 13 and

Table 13. (a). Graphanalysis of GENCOs to Distribution Companies (DISCOs) generation response of Area1 with various controllers under bilateral contract. (b). Graph analysis of DISCOs to GENCOs generation response of Area1 with various controllers under bilateral contract.

(a)
Controller	Delay Time	Rise Time	Peak Overshoot Time	Settling Time
Conventional	3.130	5.630	6.255	23.39
PSO	3.116	5.608	6.231	21.17
MFO	3.056	5.514	6.112	21.01
BMFO1	3.056	5.514	6.112	21.43
BMFO2	3.044	5.478	6.088	21.43
HHO	3.006	5.411	6.012	20.86
(b)
Controller	Delay Time	Rise Time	Peak Overshoot Time	Settling Time
Conventional	4.508	8.114	9.015	22
PSO	4.335	7.799	8.665	21.430
MFO	4.335	7.799	8.665	21.3
BMFO1	4.258	7.663	8.514	21.02
BMFO2	4.258	7.663	8.514	20.1
HHO	4.224	7.600	8.445	19.8

a. The comparative outcomes obtained show improved results of the HHO controller with delay time 3.006 s and settling time of 20.86 s.

Dynamic response of Distribution Companies (DISCOs) to GENCOs generation of Area1 with respect to time in seconds with various controllers like conventional, PSO, MFO, BMFO1, BMFO2, and HHO under bilateral contract are compared in Figure 13 and Table 13b. The comparative outcomes obtained show improved results of the HHO controller with delay time 4.224 s and settling time of 19.8 s.

Dynamic response of GENCOs to DISCOs generation of Area2 with respect to time in seconds with various controllers like conventional, PSO, MFO, BMFO1, BMFO2, and HHO under bilateral contract are compared in Figure 14 and

Table 14. (a). Graph analysis of GENCOs to DISCOs generation response of Area2 with various controllers under bilateral contract. (b). Graph analysis of DISCOs to GENCOs generation response of Area2 with various controllers under bilateral contract.

(a)
Controller	Delay Time	Rise Time	Peak Overshoot Time	Settling Time
Conventional	3.114	5.591	6.212	23.39
PSO	2.943	5.297	5.885	21.17
MFO	2.998	5.397	5.996	21.01
BMFO1	2.857	5.143	5.714	21.43
BMFO2	2.663	4.792	5.324	21.43
HHO	2.506	4.511	5.012	20.86
(b)
Controller	Delay Time	Rise Time	Peak Overshoot Time	Settling Time
Conventional	6.434	11.579	12.865	22
PSO	6.006	10.811	12.012	21.430
MFO	6.066	10.901	12.112	21.3
BMFO1	6.066	10.901	12.112	21.02
BMFO2	5.885	10.595	11.770	20.1
HHO	5.564	10.001	11.112	19.8

a. The comparative outcomes obtained show improved results of the HHO controller with delay time 2.506 s and settling time of 20.86 s.

Dynamic response of DISCOs to GENCOs generation of Area2 with respect to time in seconds with various controllers like conventional, PSO, MFO, BMFO1, BMFO2, and HHO under bilateral contract are compared in Figure 14 and Table 14b. The comparative outcomes obtained show improved results of the HHO controller with delay time 5.564 s and settling time of 19.8 s.

5.3. Contract Violation Case

The various dynamic responses with penetration of wind energy in two areas and the tie-line are represented in Figure 15, Figure 16, Figure 17, Figure 18, Figure 19 and Figure 20.

Dynamic response of Area1 frequency with respect to time in seconds with various controllers like conventional, PSO, MFO, BMFO1, BMFO2, and HHO under contract violation case are compared in Figure 15 and

Table 15. Graph analysis of Area1 frequency with various controllers under contract violation case.

Controller	Delay Time	Rise Time	Peak Overshoot Time	Settling Time
Conventional	−2.414	−4.205	−4.668	23.39
PSO	−1.508	−2.715	−3.015	21.17
MFO	−1.503	−2.705	−3.005	21.01
BMFO1	−1.003	−1.805	−2.005	21.43
BMFO2	−0.978	−1.766	−1.955	21.43
HHO	−0.835	−1.499	−1.665	20.86

. The comparative outcomes obtained show improved results of the HHO controller with delay time −0.835 s and settling time of 20.86 s.

Dynamic response of Area2 frequency with respect to time in seconds with various controllers like conventional, PSO, MFO, BMFO1, BMFO2, and HHO under contract violation case are compared in Figure 16 and

Table 16. Graph analysis of Area2 frequency with various controller under contract violation case.

Controller	Delay Time	Rise Time	Peak Overshoot Time	Settling Time
Conventional	0.0066	0.0118	0.0131	23.39
PSO	0.0031	0.0047	0.0052	21.17
MFO	0.0045	0.0074	0.0082	21.01
BMFO1	0.0026	0.0046	0.0051	21.43
BMFO2	0.0027	0.0048	0.0053	21.43
HHO	0.0023	0.0041	0.0045	20.86

. The comparative outcomes obtained show improved results of the HHO controller with delay time 0.0023 s and settling time of 20.86 s.

Dynamic response of deviation in actual tie-line power flow with respect to time in seconds with various controllers like conventional, PSO, MFO, BMFO1, BMFO2, and HHO under contract violation case are compared in Figure 17 and

Table 17. Graph analysis of deviation in actual tie-line power flow with various controller under contract violation case.

Controller	Delay Time	Rise Time	Peak Overshoot Time	Settling Time
Conventional	0.0106	0.0190	0.0211	23.39
PSO	0.0043	0.0077	0.0085	21.17
MFO	0.0063	0.0114	0.0125	21.01
BMFO1	0.0033	0.0059	0.0065	21.43
BMFO2	0.0031	0.0055	0.0061	21.43
HHO	0.0029	0.0053	0.0058	20.86

. The comparative outcomes obtained show improved results of the HHO controller with delay time 0.0029 s and settling time of 20.86 s.

Dynamic response of deviation in tie-line error with respect to time in seconds with various controllers like conventional, PSO, MFO, BMFO1, BMFO2, and HHO under contract violation case are compared in Figure 18 and

Table 18. Graph analysis of deviation in tie-line error with various controllers under contract violation case.

Controller	Delay Time	Rise Time	Peak Overshoot Time	Settling Time
Conventional	1.284	2.311	2.568	22
PSO	1.278	2.299	2.555	21.430
MFO	1.277	2.298	2.554	21.3
BMFO1	1.277	2.297	2.553	21.02
BMFO2	1.277	2.296	2.552	20.1
HHO	1.276	2.295	2.551	19.8

. The comparative outcomes obtained show improved results of the HHO controller with delay time 1.276 s and settling time of 19.8 s.

Dynamic response of GENCOs to DISCOs generation of Area1 with respect to time in seconds with various controllers like conventional, PSO, MFO, BMFO1, BMFO2, and HHO under contract violation case are compared in Figure 19 and

Table 19. (a). Graph analysis of GENCOs to DISCOs generation response of Area1 with various controllers under contract violation case. (b). Graph analysis of DISCOs to GENCOs generation response of Area1 with various controllers under contract violation case.

(a)
Controller	Delay Time	Rise Time	Peak Overshoot Time	Settling Time
Conventional	3.130	5.630	6.255	23.39
PSO	3.116	5.608	6.231	21.17
MFO	3.056	5.514	6.112	21.01
BMFO1	3.056	5.514	6.112	21.43
BMFO2	3.044	5.478	6.088	21.43
HHO	3.006	5.411	6.012	20.86
(b)
Controller	Delay Time	Rise Time	Peak Overshoot Time	Settling Time
Conventional	4.508	8.114	9.015	22
PSO	4.335	7.799	8.665	21.430
MFO	4.335	7.799	8.665	21.3
BMFO1	4.258	7.663	8.514	21.02
BMFO2	4.258	7.663	8.514	20.1
HHO	4.224	7.600	8.445	19.8

a. The comparative outcomes obtained show improved results of the HHO controller with delay time 3.006 s and settling time of 20.86 s.

Dynamic response of DISCOs to GENCOs generation of Area1 with respect to time in seconds with various controllers like conventional, PSO, MFO, BMFO1, BMFO2, and HHO under contract violation case are compared in Figure 19 and Table 19b. The comparative outcomes obtained show improved results of the HHO controller with delay time 4.224 s and settling time of 19.8 s.

Dynamic response of GENCOs to DISCOs generation of Area1 with respect to time in seconds with various controllers like conventional, PSO, MFO, BMFO1, BMFO2, and HHO under contract violation case are compared in Figure 20 and

Table 20. (a). Graph analysis of GENCOs to DISCOs generation response of Area2 with various controllers under contract violation case. (b). Graph analysis of DISCOs to GENCOs generation response of Area2 with various controllers under contract violation case.

(a)
Controller	Delay Time	Rise Time	Peak Overshoot Time	Settling Time
Conventional	3.114	5.591	6.212	23.39
PSO	2.943	5.297	5.885	21.17
MFO	2.998	5.397	5.996	21.01
BMFO1	2.857	5.143	5.714	21.43
BMFO2	2.663	4.792	5.324	21.43
HHO	2.506	4.511	5.012	20.86
(b)
Controller	Delay Time	Rise Time	Peak Overshoot Time	Settling Time
Conventional	6.434	11.579	12.865	22
PSO	6.006	10.811	12.012	21.430
MFO	6.066	10.901	12.112	21.3
BMFO1	6.066	10.901	12.112	21.02
BMFO2	5.885	10.595	11.770	20.1
HHO	5.564	10.001	11.112	19.8

a. The comparative outcomes obtained show improved results of the HHO controller with delay time 2.506 s and settling time of 20.86 s.

Dynamic response of DISCOs to GENCOs generation of Area1 with respect to time in seconds with various controllers like conventional, PSO, MFO, BMFO1, BMFO2, and HHO under contract violation case are compared in Figure 20 and Table 20b. The comparative outcomes obtained show improved results of the HHO controller with delay time 5.564 s and settling time of 19.8 s.

The relative resultsillustrate that an HHO-based PI controller fabricated finer conclusions when compared with structures including conventional PSO, MFO, BMFO1, and BMFO2-centered PI controllers. The peaks’ overshoot/undershoot and settling period were compact and perturbations were quickly covered in the structure having a HHO-focused PI controller [28].

6. Conclusions

The power system’s consistent operation necessitates a constant balancing of source and load as per recognized operating criteria. MFO is a very promising and interesting algorithm due to its advantages like fast searching speed and simplicity, but hasdrawbacks like getting stuck in bad local optima because it focuses on exploitation rather than exploration. Therefore, it is of great importance to research and put forward advanced optimization algorithms like BMFO1, BMFO2, and HHO with better performance to supplement the algorithm. This paper discussed major issues regarding the addition of RESs into frequency regulation power structure, which is most noticeable recently. This work briefly studied the utmost significant problems with the current accomplishments reported in the literature with different contract cases of unilateral, bilateral, and contract violation. The analysis results showed the improved consequences compared to conventional regulators with the help of modern soft computing techniques like the Harris hawks optimizer. The revised LFC model was also presented, which maintains the system frequency without any steady-state error, unlike conventional PI and moth flame optimizer. It instantaneously responds to different load disturbances and makes the system stable within a short time.