Development and application of equilibrium optimizer for optimal power flow calculation of power system

This paper proposes an enhanced version of Equilibrium Optimizer (EO) called (EEO) for solving global optimization and the optimal power flow (OPF) problems. The proposed EEO algorithm includes a new performance reinforcement strategy with the Lévy Flight mechanism. The algorithm addresses the shortcomings of the original Equilibrium Optimizer (EO) and aims to provide better solutions (than those provided by EO) to global optimization problems, especially OPF problems. The proposed EEO efficiency was confirmed by comparing its results on the ten functions of the CEC’20 test suite, to those of other algorithms, including high-performance algorithms, i.e., CMA-ES, IMODE, AGSK and LSHADE_cnEpSin. Moreover, the statistical significance of these results was validated by the Wilcoxon’s rank-sum test. After that, the proposed EEO was applied to solve the the OPF problem. The OPF is formulated as a nonlinear optimization problem with conflicting objectives and subjected to both equality and inequality constraints. The performance of this technique is deliberated and evaluated on the standard IEEE 30-bus test system for different objectives. The obtained results of the proposed EEO algorithm is compared to the original EO algorithm and those obtained using other techniques mentioned in the literature. These Simulation results revealed that the proposed algorithm provides better optimized solutions than 20 published methods and results as well as the original EO algorithm. The EEO superiority was demonstrated through six different cases, that involved the minimization of different objectives: fuel cost, fuel cost with valve-point loading effect, emission, total active power losses, voltage deviation, and voltage instability. Also, the comparison results indicate that EEO algorithm can provide a robust, high-quality feasible solutions for different OPF problems.


Introduction
Achieving an optimal power flow (OPF) represents an essential and complex non-linear optimization problem in power systems [1]. From an optimization standpoint, the OPF problem involves minimization of an objective, such as voltage instability, fuel cost, and total emission. For this minimization, researchers have applied algorithms in order to obtain an optimized adjustment of control variables that respect constraints, i.e., equality and inequality of operating conditions [2]. These variables include the real power generator, voltages of generation buses, tap ratios of Essam H. Houssein essam.halim@mu.edu.eg Extended author information available on the last page of the article. transformers, and reactive powers of shunt compensation capacitors [3].
Recently, several metaheuristic algorithms (MAs) [4] for solving the OPF problems have been proposed [5]. These include the adaptive multiple team perturbation-guiding Jaya (AMTPG-Jaya) algorithm [6], modified sine-cosine algorithm (MSCA) [7], modified grasshopper optimization algorithm (MGOA) [5], moth swarm algorithm (MSA) [8], electromagnetism mechanism algorithm (EM) [9], and colliding bodies optimization (CBO) [10]. Although these algorithms have solved the same kind of OPF problems, their objective functions, i.e., optimization goals, were different. Different goals lead to different optimized solutions and can therefore influence the optimization performance. Generally, optimization algorithm performance refers to the quality of the optimized solution and its computation time, i.e., time required for convergence of the algorithm. Although many MAs have provided satisfactory results, / Published online: 18 July 2022 Applied Intelligence (2023) 53:7232-7253 optimization problems have become increasingly challenging (as the number of optimized variables has increased), while satisfying many constraints and requirements. However, in spite of the advantages of algorithms, numerous existing meta-heuristic optimization algorithms do not constantly guarantee the globally optimum solution. In addition, owing to the variability of objectives Where, diverse functions are used for formulating the OPF problem, no algorithm can be considered the better-quality in solving all variants of OPF. Therefore, developing new meta-heuristic algorithms which can effectively handle diverse OPF formulations, is necessary [6]. Combining MAs -often referred to as hybridization-is effective in addressing the current optimization challenges [11]. Although hybridization improves optimization performance, it should be performed with adequate algorithms. Thus, selecting the algorithms is an important step. And common practice is to select them based on their standalone performance. Another way to improve algorithm performance is by adding optimization components to an original algorithm.
Therefore, to develop a more effective algorithm (than the methods typically employed) for solving OPF problems, we have studied recent algorithms and features. The Equilibrium Optimizer (EO) has attracted the attention of many researchers -in approximately one year, this method has been cited 500 times. In [12], the EO performed better than several other algorithms. This algorithm has been validated for over 58 benchmark functions, including composite functions and functions from the Congress On Evolutionary Computation 2017 (CEC '17), which are considered challenging optimization problems. To further validate its efficiency, the EO has also been applied to three canonical engineering problems, i.e., welded beam design, pressure vessel design, and tension/compression spring design. The EO optimized results yielded superior performance compared with those of seven other MAs. Although the EO has yielded promising results, the algorithm has some drawbacks. For example, depending on the optimization problem, slow convergence speed, convergence to a local minimum, performance dependence on algorithm parameters, and difficulty in achieving a balance between the exploration and exploitation phases, have been reported [13].
In this regard, the EO has gained much popularity in recent days in several fields of engineering and complex applications. Authors in [13] used EO in a similar context for solving combinatorial, global, engineering, and Multi-Objective problems. Authors in [14] applied Opposition Based Learning (OBL) at the initialization phase of EO for parameters identification of photovoltaic modules. In [15], authors combined the dimension learning hunting (DLH) with EO for multi-thresholding based COVID-19 CT images. Authors in [16], the support vector regression (SVR) method with equilibrium optimizer (EO) is combined for stock market prediction. In [17], authors developed a new variant of EO called general learning equilibrium optimizer (GLEO), they utilizes a general learning strategy to explore the promising regions, the GLEO is employed as a wrapper feature selection method, to select a subset of informative biological dataset's features. Authors in [18] proposed an enhanced EO version using ReliefF algorithm and the local search strategy, the introduced feature selection algorithm, is tested on 16 UCI datasets and 10 biological datasets. In [19], authors introduced an adaptive variant of EO called LWMEO for solving the engineering design problems, the Lévy flight random walk is utilized to enhance the traditional EO exploration, and spiral encirclement mechanism to enhance the exploitation process. In [20], an improved variant of the EO (IEO), to optimize the optimal power flow (OPF) problem, the IEO uses the chaotic equilibrium pool to improve the information sharing between individuals. In [21], authors suggested an algorithm that combines the a modified version of the EO and the extreme learning machine (ELM). To enhance the exploratory search, a gaussian mutation method is incorporated to the original EO.
Some studies have focused on overcoming these shortcomings. Recently, the Lévy Flights (LF) algorithmic feature has yielded excellent results in improving MA performance [22], and consequently has attracted attention from optimization algorithm developers. Indeed, LF has been integrated into the algorithms of some MAs, such as Grey wolf optimizer (GWO) [23], Particle swarm optimization [24,25], Evaporation rate water cycle algorithm [26], Whale Optimization Algorithm (WOA) [27], Chimp optimization algorithm [28], marine predators algorithm [29], and Lévy flight distribution [30]. The results have shown that (in general) LF improves standard MAs by strengthening the local search, escaping local minima, enhancing the convergence speed, or improving the exploration-exploitation balance. With these benefits in mind, incorporating LF features into an algorithm seems a promising avenue for addressing the aforementioned shortcomings. Therefore, the motivations of this paper are to: 1. Device a methodology for reinforcing the exploration and exploitation phases of MAs, thereby providing improved solutions to optimization problems; 2. Investigate the OPF and provide enhanced solutions; 3. Propose an optimization algorithm that exhibits better performance than recent and high-performance algorithms.
Thus, this paper proposes an enhanced EO (EEO), which includes an LF component and a new reinforcement strategy for solving OPF problems more efficiently than other methods. In particular, the proposed enhancing method consists of three components for improving the local and global searches. The aim is to reduce the potential weaknesses (such as premature convergence, unbalanced exploration and exploitation phases, and convergence to a local optimum) of the standard EO. Thus, the proposed algorithm aims to solve OPF problems more efficiently (than other methods) and to serve as a high-performance optimization tool [31]. To achieve this goal, the efficiency of the proposed EEO is evaluated on ten benchmark functions of the CEC'20 test suite. The proposed EEO is then used to solve the OPF problem of a standard IEEE 30-bus system. The fuel cost, fuel cost with value-point loading effect, emission, total active power losses, voltage stability enhancement, and voltage deviation are all individually optimized. Simulation results are compared with the results of the original EO, some of the most recent algorithms, and highperformance optimizers and winners of IEEE CEC competitions including; Moth-flame optimization algorithm (MFO), Sine Cosine Algorithm (SCA), Whale Optimization Algorithm (WOA), Grey wolf optimizer (GWO), Harris hawk optimization algorithm (HHO), Black Widow Optimization Algorithm (BWO), Evolution Strategy with Covariance Matrix Adaptation (CMA-ES), Ensemble Sinusoidal Differential Covariance Matrix Adaptation (LSHADE cnEpSin), Improved Multi-Operator Differential Evolution algorithm (IMODE), Adaptive Gaining Sharing Knowledge (AGSK) and the original EO.
In summary, the features that distinguish our work from previous studies are as follows: -We proposed an enhanced algorithm, i.e., the EEO, which includes a new exploitation-exploration method. -The EEO performance is analyzed over 10  The rest of this paper are organized as follows: Section 2 presents some preliminaries about EO and other used enhancement methods i.e., the original EO, the LF method. Section 3 presents the details of the proposed EEO algorithm, and its components. Section 4 introduces the formulation of OPF problem, i.e., the mathematical model and constraints. Section 5 presents the results obtained and analyses performed by the proposed EEO and competitive algorithms on the CEC'20 test suite and the OPF problems. Section 6 concludes the paper.

Preliminaries
The different components of the proposed algorithm, i.e., the original EO, the LF feature, are described in the following subsections.

Equilibrium optimizer (EO)
Inspired by physics observations, the authors of [12] have proposed an EO. Specifically, EO is based on the physics laws governing the balance of concentrations of nonreactive constituents in a controlled volume. An equation defines the conservation of mass that enters and leaves a specific volume and the system always tends to an equilibrium point. The EO algorithm is based on the ability to reach this point. Indeed, the algorithm tries to stabilize the concentration within the system. The three main mathematical steps of EO are: 1) Initialization, 2) Equilibrium pool and candidates, and 3) Concentration update (we describe these steps below and refer readers to [12] for additional details about the EO).
Initialization Similar to MAs based on population evolution, EO generates a population randomly. The population consists of particles and a uniform distribution is obtained. Particles are defined by concentration vectors. The initial population is generated from: Where, P initial i is the vector corresponding to the initial concentration of particle i, P max and P min are the upper and lower bounds, respectively, n is the number of particles in the population, and rand i generates a random value ∈ [0, 1].

Equilibrium pool and candidates
From an optimization standpoint, EO employs a pool of five particles to achieve the unknown state of equilibrium, which represents the optimal solution. The pool is composed of the four bestso-far particles for diversification purposes. The average of these four particles is employed for the exploitation process. The pool is defined as: − → P eq = − → P eq(1) , − → P eq(2) , − → P eq(3) , − → P eq(4) , − → P eq(avg) . (2) Concentration update At each iteration, the EO updates the particle population through the following equation: Where, F influences the exploration-exploitation balance and is defined as follows: Where, λ is a random value ∈ [0, 1], and and the value of t decreases with increasing iteration number iter, as follows, iter is the current iteration and Max iter is the maximum number of iterations. The constant a 2 controls the exploitation; as a 2 increases, the intensification process improves, but the exploration capability decreases. The vector t 0 is computed as follows: Where, the constant a 1 controls the diversification. The term sign(r − 0.5) designates the diversification and intensification directions. The exploration ability increases with increasing a 1 value, but the exploitation capability decreases. The vector R is referred to as the generation rate and is computed as follows: Where, − − → RCP is : Where, r 1 and r 2 are random values ∈ [0, 1] and RP is a variable that also influences the exploitation-exploration balance.

The proposed EEO algorithm
Initialization The proposed EEO integrates the abovementioned strategies. Particularly, the EEO initializes its population with LF distribution as follows: (9) Where, P initial i are the initial values of the i th particle, P min and P max are the lower and upper bounds, respectively. Levy (β) is LF random walk described. Previous work has confirmed that LF allows effective coverage of the search region, and hence candidate solutions will most likely converge to (near)-optimal solutions. Indeed, a previous study [32] has demonstrated that LF is an effective mechanism for escaping regions with local minima (even those with many deep local minima).
Reinforce exploration During the exploration stage, algorithm particles search the problem space broadly in order to identify promising areas. The exploration performed by the original EO algorithm involves only searches near the best particle, i.e., P eq , and can therefore be improved. Thus, this work proposes a new reinforcement exploration method that mutates the search particles by selecting two particles in the population as follows: 1. Based on the fitness value, the population is divided into two parts, i.e., the best and worst solutions; 2. The P r1 solution from the best fitness solutions and P r2 from the worst fitness solutions are selected via the tournament selection method, which is defined by: Where, x, y operators retain the stochastic nature and define the convergence direction of search particles, when generating P t solution; x and y are generated randomly through the following equations: Where, rand is a random number generated within average ∈ [0, 1].

Reinforce exploitation
During the exploitation process, the particles are crowded for searching around the identified promising areas from early exploration. As illustrated below, search-particle evolution is achieved by minimizing the distance from the best agent P eq . (3) of the original EO algorithm is slightly adapted in the proposed EEO and is given as follows: Where, vectors F and λ are defined by (4) and (3) respectively, and P eq(j) refers to the j-th portion of the best particle, P (i, j) refers the j-th portion of the i-th particle.

Balancing exploration and exploitation
During the first iterations, strengthening the exploration phase is essential for algorithm identification of promising areas in the problem space, while the subsequent iterations exploit the identified areas. Therefore, the following strategy is proposed for balancing the exploration and exploitation phases: Where, z is computed as follow: The value of z decreases with increasing iteration number iter; iter is the current iteration and Max iter is the maximum number of iterations. The variable a2 is a constant that controls the exploitation step. During the optimization process, the particles evolve by mutating some parts of each particle and keeping the best information parts; indeed, (14) retains the best parts of each particle P i,j .
Algorithm 1 presents the pseudo-code of the EEO algorithm, Where, iter refers to the current iteration number and Max iter refers to the total number of iterations.

Computation of complexity
The time complexity of the EEO relies mainly on the process of updating solutions' positions. Therefore, it can be formulated as follows: The Hence, the overall EEO time complexity in polynomial order.

General structure of OPF
We modeled the OPF problem as follows [33]: subject to: Where, f is the objective functions, g and h are sets of equality constraints and inequality constraints of the power system network, which are voltage and angle of load buses respectively. And x, which represents various state variables, is defined as follows: Where, P G 1 is the active power of the generator at the slack bus, V L i is the voltage magnitude of the i th load bus, Q G i is the reactive power output of the i th generator, and S l j is the line loading of the j th line. NL is the number of load buses. NG is the number of generators and ln is the number of transmission lines. Furthermore, u, the vector of control variables, which are included real and reactive power outputs from generators, bus voltages, series and/or shunt capacitors (reactors), tap-changer transformers setting defined as [33]: Where, P G i is the i th bus generator real power excluding the swing generator, V G i is the voltage magnitude of the i th generator, Q c d is the shunt compensation of the d th bus, T k is the k th branch transformer tap, NT is the number of regulating transformers, and NC is the number of shunt compensators. Any value within its range can be assumed as a control variable.

Objective functions for the OPF
The proposed EEO performance is evaluated over six case studies with various objective functions for a standard IEEE 30-bus system [33].

Case 1: fuel cost minimization
We can relate the total fuel cost ($/h) to the generated power (MW) as follows: Where, a i , b i , and c i aare the cost coefficients of the thermal generators P G i .

Case 2: fuel cost with value-point loading effect minimization
Minimization of the total fuel cost with value-point loading effect is achieved through the following relation: Where, d i and e i are the cost coefficients of the i th thermal generators. Table 1 lists the coefficient values presented in [34].

Case 3: emission minimization
The third objective function is used to minimize the emission produced by the thermal generation units: Where, α i , β i , γ i , ω i , μ are the emission coefficients of the thermal generators. The values presented in [34] are listed in Table 2.

Case 4: total real power losses minimization
The fourth objective minimizes the total real power loss and is expressed as follows: Where, G q (ij ) is the conductance transfer of branch q (ij ) and δ ij is the difference in voltage angles. V i is the voltage at bus i and V j is the voltage at bus j .

Case 5: voltage instability minimization
The power system stability refers to the ability of the system to maintain bus voltages within admissible limits. The voltage stability index (L-index) of each bus, which is an accepted metric for assessing this stability, is given as: LG and Y LL are submatrices of the admittance matrix at a specific bus. Thus, we can express the voltage stability as follows [35]: Where, L j is the L-index of the j th load bus.

Case 6: voltage deviation minimization
The last objective function is used to minimize the cumulative deviation of voltages obtained for the entire load bus:

Constraints
Usually, the OPF constraints are split into two categories: i) equality constraint and ii) inequality constraints [5]. The equality constraints are as follows: Where, θ ij is the difference in voltage angles. P Di and Q Di are the active load demand and the reactive load demand, respectively. G ij is the transfer conductance and B ij is the susceptance.
The inequality constraints, which are associated with five parts of the power system, are given as follows: -Generator constraints: -Transformer tap setting constraints -Shunt compensator constraints -Voltages at load bus constraints -Transmission line loading constraints We use the following penalty function to ensure feasible solutions, where all the constraints are respected: Where, K p , K q , K v , and K s are penalty factors. In this study, K p = K q = K v = 100, and K s = 100, 000. Figure 1 shows the flowchart of the OPF optimization performed by the proposed EEO.

Simulations results and discussion
Before applying the proposed algorithm to OPF problems, we assess the proposed EEO efficiency on the IEEE Congress on Evolutionary Computation 2020 (CEC'20) [36]. The simulation results of the proposed EEO are compared with those obtained from various classes of existing optimization methodsincluding; LSHADE cnEpSin, CMA-ES, IMODE, AGSK, MFO, SCA, WOA, GWO, HHO, BWO, and EO.

Parameter settings
To conduct a fair comparison, the EEO algorithm and the other competitors are investigated through 30 runs. The function evaluations (FEs) number is set to 150,000 for all considered problems. Table 3 shows the parameters' setting for each algorithm. The parameters of the proposed algorithm and competitive algorithms are tuned first. In order to get the suitable parameter values corresponding to the best performance when applying the test methods, the taguchi robust design parameter is used. The Taguchi [37] method utilizes the Orthogonal Array (OA) and the mean analysis to investigate the effects of the algorithm's parameters based on the statistical analysis of experiments. The OA is a fractional factorial matrix of numbers arranged so that each row represents the level of the factors in each run and each column represents a specific factor that can be changed from each run.

Statistical results and analysis
The algorithm results on the CEC'20 functions are compared. In particular, the efficiency of each algorithm is measured with the average of the best solutions obtained at each run and the corresponding standard deviation (STD). Table 4 presents the average and STD values of each algorithm for functions of 10-dimension, i.e., Dim = 10.
The best results are shown in boldface.

Convergence behavior analysis
The convergence of the algorithms is evaluated (see Fig. 2. As shown in the figure, the EEO algorithm converges to (near)-optimal solutions faster than most of the other algorithms, and is therefore a viable optimization technique for problems requiring fast computation, such as online optimization.

Wilcoxon rank test analysis
Wilcoxon's rank-sum test is performed to show the significance of the achieved results. Wilcoxon test demonstrates that the algorithm behavior is not random. Although, MAs are stochastic ones the probable performance is expected to be precise. For more details relating to Wilcoxon's test, interested reader can refer to [38]. Wilcoxon's rank sum test based on the average values is used to investigate the difference between EEO and the other optimization techniques. Table 5 compares the results obtained by EEO and the other algorithms. The "Better" column represents the sum of ranks on 10 problems when the EEO is better than the other methods. However, the "Worse" column denotes the sum of ranks on 10 functions when the EEO is worse than the other methods. The "p-value" is the significance test that decides whether the similarity hypothesis should be rejected. The significance test level should be less than 5%. Due to Friedman rank, EEO performs better than the MFO, SCA, WOA, GWO, HHO, BWO, and EO, in the same context, EEO compete with the CEC'20 competition winners including the IMODE and AGSK. In specific way,

Qualitative metrics analysis
The qualitative analysis of the proposed EEO algorithm are illustrated in Fig. 3. Notably, the agent's behaviors are displayed in Fig. 3, which include 2D views of the functions, search history, average fitness history, and convergence curves. The qualitative analysis depicts the exploration/exploitation balance of the optimization algorithm, through various metrics, especially the fluctuation of solutions diversity [39,40], over the course of the optimization process.
The following points are worthwhile from the qualitative analysis: -In terms of domain's topology -functions in 2D views: The first column of Fig. 3 Fig. 3 explains the average fitness history, such that the averages of fitness value as a function of the iteration number. This metric concentrates the light over the general behavior of the particles through the optimization process. Particularly, the result history curves are in a decreasing pattern, which is referring the improve of particles at each iteration. This stable improvement confirms a cooperative searching behavior between the EEO particles. -In terms of population diversity: The diversity plot curves are presented in the last column of Fig. 3, these curves depict the average distance between the population particles during the optimization process. The result diversity curves show that, at the iterations, the particles are most likely exploring the search space with a high diversity value. While the optimization progresses, the particles converge towards the best solution, in the exploitation phase, matched with a decreasing in the diversity value. The stable interchange in the particles diversity boost the exploration/exploitation balance strategy in the EEO algorithm.
In summary, from the results obtained, the following points can be observed: -The proposed EEO reached the optimal value for F1, F3, F4, F6 and F7 and near-optimal value for F5, F8, and F10. These results strongly suggest that the proposed EEO could perform well on other functions with similar characteristics. -The proposed EEO reached equivalent or better results than the other algorithms on most CEC'20 functions, as shown in Table 4. -The Wilcoxon's rank-sum test confirms that the EEO algorithm is statistically significant. -The convergence curves in Fig. 2 confirm that the proposed algorithm has better exploration and exploitation abilities than the other algorithms. For most functions, some of the other algorithms either get stuck in a local optimum or fail to converge to a lower value, indicating respectively poorer exploration ability and poorer exploitation capability than those of the proposed algorithm. The improved exploitation and exploration abilities result from addition of the LF and reinforcement exploration/exploitation strategies to the original EO algorithm.

Experimental series 2: applying EEO for solving OPF problems
The OPF problem-solving ability of the proposed EEO is evaluated. For this evaluation, we compared the proposed EEO to the original EO using the standard IEEE 30-bus system; Fig. 4 shows the single-line diagram of the system [41]. Its characteristics are presented in Tables 10, 11 and 12 in the Appendix. This system has 24 control variables which consist of the active power of PV buses, voltages magnitudes of generator buses, transformer ratio, and shunt reactive power compensating. Furthermore, the transformer tap and shunt reactive power compensating among the control variables both are discrete variables. We performed 20 independent runs for each objective function detailed in subsection 4.2; Table 6 lists the six case studies. We run the simulations with MATLAB R2016a on a computer that has a 2.4 GHz processor and 8 GB RAM. For both algorithms, the maximum number of iterations is 300 and the population size is 30. Table 7 presents the settings of all control and state (dependent) variables along with their allowable ranges for the best fitness value obtained for the objective function in a case study pertaining to a 30-bus test system using the proposed EEO and original EO algorithms. Active power of swing generator (PG1) and reactive power of all generators are states or dependent variables treated as constraints in the optimization. Values of these variables are listed here to show that the proposed EEO techniques duly comply with the limits of these constraining variables for the six cases. The statistical results of 20 runs for each study case Fig. 4 Single-line diagram of IEEE 30-bus system [42] performed (i.e. case 1 to case 6) using the EEO and EO algorithms are presented in Table 8. The columns indicate the best, mean, worst, and standard deviation (STD) values for the objective function in each case. The best solutions are shown in boldface.
For all objective functions, the proposed algorithm reached better solutions in regards to the average, best, worst, and STD values. Figure 5 presents the voltage profiles of the load buses obtained by EO and EEO for the six cases. Total active power loss (Ploss) 5 Voltage instability (L-index) 6 Voltage deviation (VD) We can see that EEO yields voltages within the lower and upper bounds. The safety margin (the largest absolute value of the voltage difference between the bounds and the value obtained by the algorithm) values achieved with EEO are higher than those achieved with EO. Thus, the optimized system by EEO will tolerate higher voltage perturbations than the system optimized by the original EO. Figures 6 to 8 compare the convergence curves of both algorithms; except for case 3, in the beginning, the original EO has lower objective function values than EEO, but after fewer than 50 iterations, the proposed EEO reaches lower values. Narrow data distributions are obtained for all cases (highest STD value obtained by EEO: 0.305491785). Figure 9 shows the 20 independent run distributions obtained by the proposed EEO for Case 1.
Fuel cost minimization is also performed. The EEO obtained a fuel cost of 800.415 $/h, which is lower than the value obtained by the original EO (800.433 $/h) and other techniques, as shown in Table 9. Figure 6a compares the convergence characteristics of EEO and EO; the EEO performs better than the EO. Table 7 presents the statistical analysis of the optimized solutions. In addition, Figs 6b and 7a show the load bus voltages of the IEEE 30-bus system for the cases considered.    This case aims to minimize the fuel cost associated with the value-point loading effect. Figure 6b shows the convergence of the studied algorithms; the EEO outperforms the original EO in terms of achieving the best solution. Table 7 shows that total fuel costs of 832.1817 $/h and 832.1969 $/h are achieved for EEO and EO, respectively, indicating the superiority of the proposed EEO to the original EO. Table 9 compares the best values of the costs obtained by the EEO with other counterparts.

Case 3: minimizing emission
This case is aimed at minimizing the total emission, thereby reducing pollution. The results obtained for this case are provided in Table 7, which reveals that EEO and EO achieve optimal total emissions of 0.2048212 t/h and 0.20482522 t/h, respectively. Figure 7a shows the convergence curves obtained by both algorithms. Furthermore, Table 9 compares the best result obtained by the proposed algorithm with the results available in the literature; the EEO provides one of the best results.

Case 4: minimizing total active power loss
In this case, the OPF solutions are optimized by considering the optimal total active power loss values. The convergence characteristics of the considered algorithms for this case are shown in Fig. 7b. From Table 9, the optimal power losses obtained by the EEO and EO are 3.088974 MW and 3.096271 MW, respectively. Compared with recent techniques (see Table 9), the EEO provides an approved solution and guarantees that all limits are always respected, unlike some of the published algorithms.

Case 5: voltage stability enhancement
The main objective of this case is to obtain the best values of voltage stability enhancement. The optimal settings of design variables for the finest stability enhancement by EEO and basic EO are listed in Table 7. As shown in the table, the L-index (the voltage stability indicator) obtained by EEO, 0.124236, is lower and therefore better than the value (0.124268) obtained by the original EO. A comparison of the EEO and EO convergence characteristics (Fig. 8a) reveals that the EEO outperforms the original EO in terms of convergence rate and optimized solution.

Case 6: minimizing the voltage deviation
The obtained results from the proposed algorithm and original algorithms for case 6 are listed in Table 7. Table 9 presents the results obtained by the proposed technique and other optimization techniques used for solving the same case. We see from the tables that, compared with the other algorithms, the proposed algorithm achieved a better solution. The convergence curves obtained from the proposed EEO and the original EO of case 6 are presented in Fig. 8b; the proposed method yields rapid convergence to the best solution.

Comparing EEO with published studies
We compared the EEO results to more than 20 published results as shown in Table 9. Compared with other methods, the EEO achieved lower function values in most cases, substantiating the EEO efficiency. Thus, at a large scale, substantial cost savings and emission reductions can potentially be achieved by EEO while improving the system stability.

Conclusion and future work
This paper proposed an enhanced version of the Equilibrium Optimizer (EO) called EEO, which relies mainly on reinforcing the algorithm exploration and exploitation process. The proposed EEO is applied for obtaining improved solutions to global optimization problems, and Optimal Power Flow (OPF) problems. During the process, we assess the performance of the proposed EEO with regard to ten functions of the CEC'20 test suite. The EEO achieved better or similar results than LSHADE cnEpSin, CMA-ES, IMODE, AGSK, MFO, SCA, WOA, GWO, HHO, BWO, and EO. Wilcoxon's rank-sum test confirms that the proposed EEO results are statistically significant. Moreover, we demonstrated the efficiency of the proposed EEO on OPF for the standard IEEE 30-bus system. We minimized different objectives, i.e., the fuel cost, fuel cost with value-point loading effect, total emission, active power loss, voltage deviation, and voltage instability. For most objectives, the EEO yielded better results than the original EO and the methods reported in 20 published studies. Substantial cost savings and emission reductions can potentially be achieved with EEO at a large scale while improving system stability. Thus, the proposed algorithm is a valuable optimization tool for engineers of power systems and a promising tool for solving more complex optimization problems than those associated with such systems. Indeed, as future work, the proposed EEO will be applied to more challenging problems (than the problem considered here), including multi-objective problems and problems such as prediction, image segmentation, Cloud Data Center [59,60], and prediction of cloud workloads [61][62][63]. Table 10 presents as in [5] the characteristics if the IEEE 30-bus system. Table 11 presents as in [64] the branch data used for the IEEE 30-bus system. Table 12 presents as in [64] the branch data used for the IEEE 30-bus system.

Conflict of Interests
The authors declare that there is no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.