# Optimal Operation of Distribution Networks
Considering Renewable Energy Sources Integration and
Demand Side Response

Ahmed Hachemi, Fares Sadaoui, Abdelhakim Saim, Mohamed Ebeed, Hossam Abbou, Salem Arif

► To cite this version:

Ahmed Hachemi, Fares Sadaoui, Abdelhakim Saim, Mohamed Ebeed, Hossam Abbou, et al.. Optimal Operation of Distribution Networks Considering Renewable Energy Sources Integration and Demand Side Response. Sustainability, 2023, 15 (24), pp.16707. 10.3390/su152416707. hal-04538463

HAL Id: hal-04538463

https://hal.science/hal-04538463

Submitted on 6 Jun 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Optimal Operation of Distribution Networks
Considering Renewable Energy Sources Integration and
Demand Side Response

Ahmed Hachemi, Fares Sadaoui, Abdelhakim Saim, Mohamed Ebeed, Hossam Abbou, Salem Arif

► To cite this version:

Ahmed Hachemi, Fares Sadaoui, Abdelhakim Saim, Mohamed Ebeed, Hossam Abbou, et al.. Optimal Operation of Distribution Networks Considering Renewable Energy Sources Integration and Demand Side Response. Sustainability, 2023, 15 (24), pp.16707. 10.3390/su152416707. hal-04538463

HAL Id: hal-04538463

https://hal.science/hal-04538463v1

Submitted on 6 Jun 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Optimal Operation of Distribution Networks Considering Renewable Energy Sources Integration and Demand Side Response

Ahmed T. Hachemi<sup>1,†</sup>, Fares Sadaoui<sup>1,‡</sup>, Abdelhakim Saim<sup>2,*</sup>†, Mohamed Ebeed<sup>3,†</sup>, Hossam E. A. Abbou<sup>4,‡</sup> and Salem Arif<sup>4</sup>

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Hachemi, A.T.; Sadaoui, F.; Saim, A.; Ebeed, M.; Abbou, H.E.A.; Arif, S. Optimal Operation of Distribution Networks Considering Renewable Energy Sources Integration and Demand Side Response. *Sustainability* **2023**, *15*, 16707. https://doi.org/10.3390/su152416707

Academic Editors:
Ron-Hendrik Hechelmann,
Henning Meschede, Florian Schlosser
and Alexander Schlüter

Received: 16 October 2023

Revised: 19 November 2023

Accepted: 7 December 2023

Published: 10 December 2023

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2023 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

**Abstract:** This paper demonstrates the effectiveness of Demand Side Response (DSR) with renewable integration by solving the stochastic optimal operation problem (OOP) in the IEEE 118-bus distribution system over 24 h. An Improved Walrus Optimization Algorithm (I-WaOA) is proposed to minimize costs, reduce voltage deviations, and enhance stability under uncertain loads, generation, and pricing. The proposed I-WaOA utilizes three strategies: the fitness-distance balance method, quasi-opposite-based learning, and Cauchy mutation. The I-WaOA optimally locates and sizes photovoltaic (PV) ratings and wind turbine (WT) capacities and determines the optimal power factor of WT with DSR. Using Monte Carlo simulations (MCS) and probability density functions (PDF), the uncertainties in renewable energy generation, load demand, and energy costs are represented. The results show that the proposed I-WaOA approach can significantly reduce costs, improve voltage stability, and mitigate voltage deviations. The total annual costs are reduced by 91%, from 3.8377 × 10⁷ USD to 3.4737 × 10⁶ USD. Voltage deviations are decreased by 63%, from 98.6633 per unit (p.u.) to 36.0990 p.u., and the system stability index is increased by 11%, from 2.444 × 10³ p.u. to 2.7245 × 10³ p.u., when contrasted with traditional methods.

**Keywords:** distribution network; Renewable Energy Sources; multi-objective function; optimal operation; Improved Walrus Optimization Algorithm; Demand Side Response

# 1. Introduction

## 1.1. Motivation

In recent years, the reliability of energy supplies to consumers has gained paramount significance. This is owing to its role as a pivotal determinant of both capacity quality and customer satisfaction within distribution networks. Enterprises wield a plethora of options to enhance distribution network reliability, including the utilization of Renewable Energy Sources (RESs). RESs, like wind and solar power, present promising low-carbon substitutes for conventional fossil fuels [1,2]. Nevertheless, the intermittent and uncertain characteristics of RESs present complex hurdles to grid design. To tackle these challenges, the integration of RESs and the implementation of Demand Side Management strategies come to the forefront. DSR, comprising compact power facilities located proximately to the point of consumption, holds the capacity to alleviate strain on centralized networks and bolster resilience through fault isolation for swift recovery from downtimes. Complementary to this, DSR techniques like real-time pricing (RTP) exhibit the potential to reduce energy consumption during peak demand intervals.

<footnote n="1">¹ Electrical Engineering Laboratory, University of Kasdi Merbah, Ouargla 30000, Algeria; hachemi.ahmedtidjani@univ-ouargla.dz (A.T.H.); sadaoui.fares@univ-ouargla.dz (F.S.)</footnote>

<footnote n="2">² IREENA Laboratory, Nantes University, 44035 Nantes, France</footnote>

<footnote n="3">³ Department of Electrical Engineering, Faculty of Engineering, Sohag University, Sohag 82524, Egypt; mebeed@eng.sohag.edu.eg</footnote>

<footnote n="4">⁴ LACoSERE Laboratory, University of Amar Telidji, Laghouat 03000, Algeria; h.abbou@lagh-univ.dz (H.E.A.A.); s.arif@lagh-univ.dz (S.A.)</footnote>

<footnote>* Correspondence: abdelhakim.saim@univ-nantes.fr</footnote>
Article

# Optimal Operation of Distribution Networks Considering Renewable Energy Sources Integration and Demand Side Response

Ahmed T. Hachemi<fnref n="1" />, Fares Sadaoui<fnref n="1" />, Abdelhakim Saim<fnref n="2" />,*, Mohamed Ebeed<fnref n="3" />, Hossam E. A. Abbou<fnref n="4" />
and Salem Arif<fnref n="4" />

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Hachemi, A.T.; Sadaoui, F.; Saim, A.; Ebeed, M.; Abbou, H.E.A.; Arif, S. Optimal Operation of Distribution Networks Considering Renewable Energy Sources Integration and Demand Side Response. *Sustainability* **2023**, *15*, 16707. https://doi.org/10.3390/su152416707

Academic Editors:
Ron-Hendrik Hechelmann,
Henning Meschede, Florian Schlosser
and Alexander Schlüter

Received: 16 October 2023

Revised: 19 November 2023

Accepted: 7 December 2023

Published: 10 December 2023

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2023 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

**Abstract:** This paper demonstrates the effectiveness of Demand Side Response (DSR) with renewable integration by solving the stochastic optimal operation problem (OOP) in the IEEE 118-bus distribution system over 24 h. An Improved Walrus Optimization Algorithm (I-WaOA) is proposed to minimize costs, reduce voltage deviations, and enhance stability under uncertain loads, generation, and pricing. The proposed I-WaOA utilizes three strategies: the fitness-distance balance method, quasi-opposite-based learning, and Cauchy mutation. The I-WaOA optimally locates and sizes photovoltaic (PV) ratings and wind turbine (WT) capacities and determines the optimal power factor of WT with DSR. Using Monte Carlo simulations (MCS) and probability density functions (PDF), the uncertainties in renewable energy generation, load demand, and energy costs are represented. The results show that the proposed I-WaOA approach can significantly reduce costs, improve voltage stability, and mitigate voltage deviations. The total annual costs are reduced by 91%, from 3.8377 × 10⁷ USD to 3.4737 × 10⁶ USD. Voltage deviations are decreased by 63%, from 98.6633 per unit (p.u.) to 36.0990 p.u., and the system stability index is increased by 11%, from 2.444 × 10³ p.u. to 2.7245 × 10³ p.u., when contrasted with traditional methods.

**Keywords:** distribution network; Renewable Energy Sources; multi-objective function; optimal operation; Improved Walrus Optimization Algorithm; Demand Side Response

# 1. Introduction

## 1.1. Motivation

In recent years, the reliability of energy supplies to consumers has gained paramount significance. This is owing to its role as a pivotal determinant of both capacity quality and customer satisfaction within distribution networks. Enterprises wield a plethora of options to enhance distribution network reliability, including the utilization of Renewable Energy Sources (RESs). RESs, like wind and solar power, present promising low-carbon substitutes for conventional fossil fuels [1,2]. Nevertheless, the intermittent and uncertain characteristics of RESs present complex hurdles to grid design. To tackle these challenges, the integration of RESs and the implementation of Demand Side Management strategies come to the forefront. DSR, comprising compact power facilities located proximately to the point of consumption, holds the capacity to alleviate strain on centralized networks and bolster resilience through fault isolation for swift recovery from downtimes. Complementary to this, DSR techniques like real-time pricing (RTP) exhibit the potential to reduce energy consumption during peak demand intervals.

<footnote n="1">¹ Electrical Engineering Laboratory, University of Kasdi Merbah, Ouargla 30000, Algeria; hachemi.ahmedtidjani@univ-ouargla.dz (A.T.H.); sadaoui.fares@univ-ouargla.dz (F.S.)</footnote>

<footnote n="2">² IREENA Laboratory, Nantes University, 44035 Nantes, France</footnote>

<footnote n="3">³ Department of Electrical Engineering, Faculty of Engineering, Sohag University, Sohag 82524, Egypt; mebeed@eng.sohag.edu.eg</footnote>

<footnote n="4">⁴ LACoSERE Laboratory, University of Amar Telidji, Laghouat 03000, Algeria; h.abbou@lagh-univ.dz (H.E.A.A.); s.arif@lagh-univ.dz (S.A.)</footnote>

<footnote>* Correspondence: abdelhakim.saim@univ-nantes.fr</footnote>
Collaboratively, the synergy of RESs and DSR stands poised to revamp the distribution system, allowing for a substantial integration of renewable energy. This transformation is indispensable for steering toward a sustainable energy future [3,4].

## 1.2. Related Work

Several studies have been carried out on the integration of renewable energies into electrical networks. For example, the authors in [5] present a novel and effective approach to the ideal siting and dimensioning of PV resources in distribution networks. In [6], a hybrid Harmony Search Algorithm approach is proposed to enhance bus voltage profiles and minimize power losses in the distribution network (DN). Integrating the Particle Artificial Bee Colony algorithm with Harmony Search addresses convergence issues. The hybrid algorithm determines optimal locations and sizes for distributed generators and capacitor banks (CB), proving superior loss reduction compared to current methods in the literature on test systems with 33 and 119 nodes. The authors in [7] present a method for optimal placement and sizing of PV and WT in distribution networks. Their approach, using a weighted aggregation Particle Swarm Optimization technique, addresses power loss minimization, voltage stability, and network security while considering the stochastic nature of renewable resources. Applied to an Indian rural distribution network, the method showcases effective enhancement of network operational benefits. In [8], a reinforcement learning-based method is proposed for optimal sizing of PV sources in unbalanced distribution networks, considering uncertainty and seasonal variations. The technique is validated on IEEE 13- and 37-bus feeders and is shown to be suitable for real-system implementation. The paper [9] addresses various technological challenges related to stability issues associated with high PV penetration into the power grid. These challenges include frequency regulation, active power reduction, reactive power injection, and energy storage. In [10], ramp-rate control algorithms have been proposed to minimize grid power fluctuations by using energy storage systems. The technical challenges arising from the extensive integration of PV systems into the electrical grid, like voltage limit violations and frequency disturbances, are explored in [11]. The optimal network reconfiguration is obtained using an enhanced whale optimization approach [12]. The mine blast method is used to find the best location and size of CB in DN [13] and the grey wolf optimization method [14]. The ant lion optimization algorithm in [15] is suggested for the best RES sizing and allocation in DN. The authors of the paper [16] have proposed a new method for optimal placement and sizing of RESs in DN, considering the uncertainties of RESs output powers and loading. In [17], the Salp Swarm Algorithm is applied to optimize the allocation of WT in DN to reduce overall energy and power loss. This approach supports the integration of commercially available renewable energy technologies, particularly wind turbines, and responds to global environmental concerns. In [18], the authors introduce a novel methodology, HHO-PSO, for optimizing the planning of RESs in active distribution systems. The results demonstrated that HHO-PSO provides superior solutions, maximizing techno-economic benefits. In [19], the authors propose a new algorithm for optimal allocation of distributed generation units in radial distribution networks. The proposed algorithm, called AEO-OBL, is a modified form of the optimization method modeled on artificial ecosystems. The authors in [20] propose an effective RESs allocation strategy for energy loss minimization in the presence of variable power demand. In [21], the study introduces a hybrid stochastic-IGDT strategy to improve the resilient functioning of intelligent DN with the incorporation of RESs and intelligent parking facilities. A comprehensive integration of storage approaches for the OOP of energy storage systems is proposed in DN [22]. In [23], a method is presented for concurrent positioning of distribution static compensators, RESs, and DN reconfiguration using stochastic multi-objective optimization. The problem is formulated as a mixed-integer, nonlinear programming problem. The authors in [24] propose a new approach for the optimal siting and sizing of wind turbines and devices for storing energy in distribution networks. The suggested methodology takes into account the impact of WTs on power quality, specifically flicker emission. In [4], the authors propose a technique for improving the reliability of DN using DSR programs and smart charging/discharging of plug-in electric vehicles (PEV). In [25], a DSR mechanism is introduced to enable energy
Collaboratively, the synergy of RESs and DSR stands poised to revamp the distribution system, allowing for a substantial integration of renewable energy. This transformation is indispensable for steering toward a sustainable energy future [3,4].

## 1.2. Related Work

Several studies have been carried out on the integration of renewable energies into electrical networks. For example, the authors in [5] present a novel and effective approach to the ideal siting and dimensioning of PV resources in distribution networks. In [6], a hybrid Harmony Search Algorithm approach is proposed to enhance bus voltage profiles and minimize power losses in the distribution network (DN). Integrating the Particle Artificial Bee Colony algorithm with Harmony Search addresses convergence issues. The hybrid algorithm determines optimal locations and sizes for distributed generators and capacitor banks (CB), proving superior loss reduction compared to current methods in the literature on test systems with 33 and 119 nodes. The authors in [7] present a method for optimal placement and sizing of PV and WT in distribution networks. Their approach, using a weighted aggregation Particle Swarm Optimization technique, addresses power loss minimization, voltage stability, and network security while considering the stochastic nature of renewable resources. Applied to an Indian rural distribution network, the method showcases effective enhancement of network operational benefits. In [8], a reinforcement learning-based method is proposed for optimal sizing of PV sources in unbalanced distribution networks, considering uncertainty and seasonal variations. The technique is validated on IEEE 13- and 37-bus feeders and is shown to be suitable for real-system implementation. The paper [9] addresses various technological challenges related to stability issues associated with high PV penetration into the power grid. These challenges include frequency regulation, active power reduction, reactive power injection, and energy storage. In [10], ramp-rate control algorithms have been proposed to minimize grid power fluctuations by using energy storage systems. The technical challenges arising from the extensive integration of PV systems into the electrical grid, like voltage limit violations and frequency disturbances, are explored in [11]. The optimal network reconfiguration is obtained using an enhanced whale optimization approach [12]. The mine blast method is used to find the best location and size of CB in DN [13] and the grey wolf optimization method [14]. The ant lion optimization algorithm in [15] is suggested for the best RES sizing and allocation in DN. The authors of the paper [16] have proposed a new method for optimal placement and sizing of RESs in DN, considering the uncertainties of RESs output powers and loading. In [17], the Salp Swarm Algorithm is applied to optimize the allocation of WT in DN to reduce overall energy and power loss. This approach supports the integration of commercially available renewable energy technologies, particularly wind turbines, and responds to global environmental concerns. In [18], the authors introduce a novel methodology, HHO-PSO, for optimizing the planning of RESs in active distribution systems. The results demonstrated that HHO-PSO provides superior solutions, maximizing techno-economic benefits. In [19], the authors propose a new algorithm for optimal allocation of distributed generation units in radial distribution networks. The proposed algorithm, called AEO-OBL, is a modified form of the optimization method modeled on artificial ecosystems. The authors in [20] propose an effective RESs allocation strategy for energy loss minimization in the presence of variable power demand. In [21], the study introduces a hybrid stochastic-IGDT strategy to improve the resilient functioning of intelligent DN with the incorporation of RESs and intelligent parking facilities. A comprehensive integration of storage approaches for the OOP of energy storage systems is proposed in DN [22]. In [23], a method is presented for concurrent positioning of distribution static compensators, RESs, and DN reconfiguration using stochastic multi-objective optimization. The problem is formulated as a mixed-integer, nonlinear programming problem. The authors in [24] propose a new approach for the optimal siting and sizing of wind turbines and devices for storing energy in distribution networks. The suggested methodology takes into account the impact of WTs on power quality, specifically flicker emission. In [4], the authors propose a technique for improving the reliability of DN using DSR programs and smart charging/discharging of plug-in electric vehicles (PEV). In [25], a DSR mechanism is introduced to enable energy
consumers to adjust their energy consumption patterns as proposed. The authors in [26] study the potential benefits of DSR resources for distribution networks. They consider two different types of customers: residential and commercial. They install data loggers on two real distribution feeders in Golpayegan city and collect load profiles. They then use these load profiles to develop two DSR programs: time-of-use and emergency DSR. They apply the DSR programs to the real network and simulate the system using the CYME software. In [27], the authors address the uncertainties introduced by the stochastic nature of RESs and DSR in distribution grids. The authors in [28] present a novel darts game theory-based optimization algorithm for microgrid formation, considering multiple objectives, topological and electrical constraints, and demand response. The modified IEEE 33-bus test system is used to validate the suggested methodology, demonstrating its superiority over existing methods and its potential to improve system resiliency under extreme events.

The addition of Cauchy mutations to optimization algorithms has been shown to improve their capability to address intricate global engineering optimization challenges. This enhancement leads to the discovery of better optimal solutions, as has been demonstrated by many research studies [29]. Furthermore, the refinement of optimization performance is achievable through the harmonization of fitness and distance considerations within the search process. This improvement has been substantiated through a diverse array of research inquiries [30].

In this context, quasi-opposite-based learning (QOBL) emerges as a widely adopted strategy to amplify population diversity within optimization algorithms. Proposed as an extension of the opposition-based learning approach, QOBL operates on a similar premise, involving the evaluation of agents at their opposite points within the parameter space. QOBL extends this principle by uniformly assessing points selected along a straight line connecting the search space's center with the opposing point of the agent's current location. The viability of QOBL holds promise for further enriching optimization algorithms' performance and diversification strategies, contributing to a holistic advancement in engineering optimization techniques [31]. Table 1 summarizes the work given and compares it to other relevant research on the optimal operation of DN with RESs.

<table><caption>**Table 1.** A comparative analysis of the presented work against other relevant studies on the optimal operation of DN.</caption><thead><tr><th rowspan="2">References</th><th>DG Type</th><th colspan="4">Uncertainty</th><th colspan="2">Improved Approach</th><th colspan="4">Objective Function</th><th colspan="3">DSR</th></tr><tr><th>Temperature</th><th>Irradiance</th><th>Price</th><th>Loading</th><th>Wind Speed</th><th>Cost of PV</th><th>Cost of WT</th><th>Cost of Energy loss</th><th>Cost of Grid</th><th>System Stability Index</th><th>Voltage Deviations</th><th></th><th></th></tr></thead><tbody><tr><td>[4]</td><td>PEV</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td></tr><tr><td>[5]</td><td>PV</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td></tr><tr><td>[6]</td><td>WT with CB</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td></tr><tr><td>[7]</td><td>WT with PV</td><td>×</td><td>✓</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>✓</td><td>×</td><td>×</td></tr><tr><td>[8]</td><td>PV</td><td>×</td><td>✓</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td></tr><tr><td>[17]</td><td>WT</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td></tr><tr><td>[18]</td><td>WT with PV</td><td>×</td><td>✓</td><td>×</td><td>×</td><td>×</td><td>×</td><td>✓</td><td>✓</td><td>✓</td><td>×</td><td>✓</td><td>✓</td></tr><tr><td>[32]</td><td>PV</td><td>×</td><td>✓</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>✓</td><td>×</td></tr><tr><td>[33]</td><td>WT with PV</td><td>×</td><td>✓</td><td>×</td><td>×</td><td>×</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>×</td></tr><tr><td>[27]</td><td>WT with PV</td><td>×</td><td>✓</td><td>×</td><td>×</td><td>×</td><td>×</td><td>✓</td><td>✓</td><td>×</td><td>×</td><td>×</td><td>✓</td></tr><tr><td>[28]</td><td>WT with PV</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>✓</td><td>✓</td><td>×</td><td>×</td><td>×</td><td>×</td><td>✓</td></tr><tr><td>This paper</td><td>WT with PV</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td></tr></tbody></table>
consumers to adjust their energy consumption patterns as proposed. The authors in [26] study the potential benefits of DSR resources for distribution networks. They consider two different types of customers: residential and commercial. They install data loggers on two real distribution feeders in Golpayegan city and collect load profiles. They then use these load profiles to develop two DSR programs: time-of-use and emergency DSR. They apply the DSR programs to the real network and simulate the system using the CYME software. In [27], the authors address the uncertainties introduced by the stochastic nature of RESs and DSR in distribution grids. The authors in [28] present a novel darts game theory-based optimization algorithm for microgrid formation, considering multiple objectives, topological and electrical constraints, and demand response. The modified IEEE 33-bus test system is used to validate the suggested methodology, demonstrating its superiority over existing methods and its potential to improve system resiliency under extreme events.

The addition of Cauchy mutations to optimization algorithms has been shown to improve their capability to address intricate global engineering optimization challenges. This enhancement leads to the discovery of better optimal solutions, as has been demonstrated by many research studies [29]. Furthermore, the refinement of optimization performance is achievable through the harmonization of fitness and distance considerations within the search process. This improvement has been substantiated through a diverse array of research inquiries [30].

In this context, quasi-opposite-based learning (QOBL) emerges as a widely adopted strategy to amplify population diversity within optimization algorithms. Proposed as an extension of the opposition-based learning approach, QOBL operates on a similar premise, involving the evaluation of agents at their opposite points within the parameter space. QOBL extends this principle by uniformly assessing points selected along a straight line connecting the search space's center with the opposing point of the agent's current location. The viability of QOBL holds promise for further enriching optimization algorithms' performance and diversification strategies, contributing to a holistic advancement in engineering optimization techniques [31]. Table 1 summarizes the work given and compares it to other relevant research on the optimal operation of DN with RESs.

<table><caption>**Table 1.** A comparative analysis of the presented work against other relevant studies on the optimal operation of DN.</caption><thead><tr><th rowspan="2">References</th><th>DG Type</th><th colspan="5">Uncertainty</th><th colspan="2">Improved Approach</th><th colspan="4">Objective Function</th><th colspan="3">DSR</th></tr><tr><th>Temperature</th><th>Irradiance</th><th>Price</th><th>Loading</th><th>Wind Speed</th><th>Cost of PV</th><th>Cost of WT</th><th>Cost of Energy loss</th><th>Cost of Grid</th><th>System Stability Index</th><th>Voltage Deviations</th></tr></thead><tbody><tr><td>[4]</td><td>PEV</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td></tr><tr><td>[5]</td><td>PV</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td></tr><tr><td>[6]</td><td>WT with CB</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td></tr><tr><td>[7]</td><td>WT with PV</td><td>×</td><td>✓</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>✓</td><td>×</td><td>×</td></tr><tr><td>[8]</td><td>PV</td><td>×</td><td>✓</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td></tr><tr><td>[17]</td><td>WT</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td></tr><tr><td>[18]</td><td>WT with PV</td><td>×</td><td>✓</td><td>×</td><td>×</td><td>×</td><td>×</td><td>✓</td><td>✓</td><td>✓</td><td>×</td><td>✓</td><td>✓</td></tr><tr><td>[32]</td><td>PV</td><td>×</td><td>✓</td><td>×</td><td>✓</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>✓</td><td>×</td></tr><tr><td>[33]</td><td>WT with PV</td><td>×</td><td>✓</td><td>×</td><td>✓</td><td>×</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>×</td></tr><tr><td>[27]</td><td>WT with PV</td><td>×</td><td>✓</td><td>×</td><td>×</td><td>✓</td><td>×</td><td>✓</td><td>×</td><td>×</td><td>×</td><td>×</td><td>✓</td></tr><tr><td>[28]</td><td>WT with PV</td><td>×</td><td>×</td><td>×</td><td>×</td><td>×</td><td>✓</td><td>✓</td><td>×</td><td>×</td><td>×</td><td>×</td><td>✓</td></tr><tr><td>This paper</td><td>WT with PV</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td></tr></tbody></table>
1.3. Contribution

As evident, the literature on addressing the OOP is extensive. However, to our knowledge, only a limited number of studies have tackled the stochastic OOP with optimal WT and solar PV system integration along with DSR while accounting for uncertainties in wind speed fluctuations, load, temperature shifts, variations in solar irradiation, and the volatility of energy purchasing prices. In summary, the research gaps identified in these studies are outlined below:

1. References [5–7] identified the optimal locations for PV and WT, but the shortcoming in these studies is the lack of consideration for uncertainty and DSR.

2. References [8,17,18,32,33] explored the issue of OOP, taking into consideration uncer- tainties in load, radiation, and wind speed. However, these studies omitted uncertain- ties related to temperatures and prices. Additionally, the Demand Side Response was not considered in these studies.

3. References [4,27,28] investigated the optimal operation problem (OOP) while relying exclusively on Demand Side Response to enhance the performance of the electrical dis- tribution network. However, a limitation of these studies is the lack of consideration for all uncertainties within the system.

In this study, we investigate the impact of integrating a system comprising photovoltaic
generators and wind turbines in conjunction with Demand Side Response.
We introduce a novel approach termed the Improved Walrus Optimization Algorithm,
aiming to address the optimization challenge within the IEEE 118 Distribution Network.
Our primary objective is twofold: to reduce costs and voltage deviations (VD) while enhanc-
ing the voltage stability index (VSI) for the duration of a day. The effectiveness of the model
is carefully evaluated, taking uncertainties into account in wind speed fluctuations, load,
temperature shifts, variations in solar irradiation, and the volatility of energy purchasing
prices. We compare our proposed method (I-WaOA) against established optimization tech-
niques, including Sand Cat swarm optimization (SCSO), Artificial Hummingbird Algorithm
(AHA), Dandelion Optimizer (DO), Harmony Search (HS), Artificial Rabbits Optimization
(ARO), Chernobyl Disaster Optimizer (CDO), Zebra Optimization Algorithm (ZOA), and
standard Walruses behavior (WaOA).

Figure 1 represents the graphical summary of this paper, incorporating the integration of RESs with DSR in line with the DN of impedance $R_{i,j} + jX_{i,j}$ connected between buses i and j.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Basic construction of the proposed RESs with a DSR-based system.</figcaption></figure>
### 1.3. Contribution

As evident, the literature on addressing the OOP is extensive. However, to our knowledge, only a limited number of studies have tackled the stochastic OOP with optimal WT and solar PV system integration along with DSR while accounting for uncertainties in wind speed fluctuations, load, temperature shifts, variations in solar irradiation, and the volatility of energy purchasing prices. In summary, the research gaps identified in these studies are outlined below:

1. References [5–7] identified the optimal locations for PV and WT, but the shortcoming in these studies is the lack of consideration for uncertainty and DSR.

2. References [8,17,18,32,33] explored the issue of OOP, taking into consideration uncertainties in load, radiation, and wind speed. However, these studies omitted uncertainties related to temperatures and prices. Additionally, the Demand Side Response was not considered in these studies.

3. References [4,27,28] investigated the optimal operation problem (OOP) while relying exclusively on Demand Side Response to enhance the performance of the electrical distribution network. However, a limitation of these studies is the lack of consideration for all uncertainties within the system.

In this study, we investigate the impact of integrating a system comprising photovoltaic generators and wind turbines in conjunction with Demand Side Response. We introduce a novel approach termed the Improved Walrus Optimization Algorithm, aiming to address the optimization challenge within the IEEE 118 Distribution Network. Our primary objective is twofold: to reduce costs and voltage deviations (VD) while enhancing the voltage stability index (VSI) for the duration of a day. The effectiveness of the model is carefully evaluated, taking uncertainties into account in wind speed fluctuations, load, temperature shifts, variations in solar irradiation, and the volatility of energy purchasing prices. We compare our proposed method (I-WaOA) against established optimization techniques, including Sand Cat swarm optimization (SCSO), Artificial Hummingbird Algorithm (AHA), Dandelion Optimizer (DO), Harmony Search (HS), Artificial Rabbits Optimization (ARO), Chernobyl Disaster Optimizer (CDO), Zebra Optimization Algorithm (ZOA), and standard Walruses behavior (WaOA).

Figure 1 represents the graphical summary of this paper, incorporating the integration of RESs with DSR in line with the DN of impedance $R_{i,j} + jX_{i,j}$ connected between buses i and j.

<figure><img src="image_4.png" /><figcaption>Figure 1. Basic construction of the proposed RESs with a DSR-based system.</figcaption></figure>
The key accomplishments of this research can be summarized as follows:

1. A novel Improved Walrus Optimization Algorithm (I-WaOA) is proposed to solve the optimal operation problem (OOP) of renewable energy resources (RESs) with Demand Side Response (DSR) in distribution networks (DN).

2. The uncertainties related to load demand, RESs generation, temperature, and energy pricing are represented through PDF and simulated using MCS.

3. The optimal size and location of wind turbines, solar photovoltaics, and DSR resources are determined simultaneously under uncertainties to minimize costs, reduce voltage deviations, and improve voltage stability in the IEEE 118-bus DN.

4. The proposed I-WaOA is demonstrated to outperform other optimization algorithms, like SCSO, AHA, DO, HS, CDO, ZOA, ARO, and standard WaOA, in handling uncertainties and solving the complex optimal operation problem.

### 1.4. Organization of Article

The paper presents the structure of an optimal planning study of a power system under uncertainties. Section 2 introduces the problem and its mathematical formulation. Section 3 discusses the modeling of systems, while Sections 4–6 present the uncertainty parameters, mathematical formulations of the WaOA, and I-WaOA problems, respectively. Section 7 is reserved for simulation results and discussions, and Section 8 concludes the paper with the implications of the results.

## 2. Formulating the Problem

This section presents the formulation of the problem, detailing both the optimal operation’s governing constraints and objective function.

### 2.1. Objective Function

Multiple objective functions are examined in this study, and they are as follows:

#### 2.1.1. The Minimization of Costs

The considered objective function includes the cost of WT ($C_{WT}$), the cost of PV units ($C_{PV}$), the annual price of energy loss ($C_{Loss}$), and the price of electricity procured from the grid ($C_{Grid}$). The total annual price can be written as:

$$C = C_{WT} + C_{PV} + C_{Loss} + C_{Grid} \quad (1)$$

In which,

$$C_{Grid} = 365 \times \sum_{h=1}^{24} P_{Grid(h)} \times U_{Grid(h)} \quad (2)$$

where $U_{Grid}$ is the price of buying energy from the network and $P_{Grid(h)}$ is the hourly power taken out of the grid.

$$C_{Loss} = 365 \times U_{Loss} \times \sum_{h=1}^{24} P_{T\_Loss(h)} \quad (3)$$

where $U_{Loss}$ represents the expenditure associated with energy loss and $P_{T\_Loss(h)}$ indicates the overall power losses during time $h$.

$$C_{PV} = C_{PV}^{inst} + C_{PV}^{O\&M} \quad (4)$$

where $C_{PV}^{O\&M}$ denotes the maintenance and operating expenses for the PV unit, while $C_{PV}^{inst}$ represents the installation cost of the PV system.

$$C_{PV}^{O\&M} = 365 \times U_{PV}^{O\&M} \times \sum_{i=1}^{Ns} \sum_{h=1}^{24} P_{PV(i,h)} \quad (5)$$
The key accomplishments of this research can be summarized as follows:

1. A novel Improved Walrus Optimization Algorithm (I-WaOA) is proposed to solve the optimal operation problem (OOP) of renewable energy resources (RESs) with Demand Side Response (DSR) in distribution networks (DN).

2. The uncertainties related to load demand, RESs generation, temperature, and energy pricing are represented through PDF and simulated using MCS.

3. The optimal size and location of wind turbines, solar photovoltaics, and DSR resources are determined simultaneously under uncertainties to minimize costs, reduce voltage deviations, and improve voltage stability in the IEEE 118-bus DN.

4. The proposed I-WaOA is demonstrated to outperform other optimization algorithms, like SCSO, AHA, DO, HS, CDO, ZOA, ARO, and standard WaOA, in handling uncertainties and solving the complex optimal operation problem.

### 1.4. Organization of Article

The paper presents the structure of an optimal planning study of a power system under uncertainties. Section 2 introduces the problem and its mathematical formulation. Section 3 discusses the modeling of systems, while Sections 4–6 present the uncertainty parameters, mathematical formulations of the WaOA, and I-WaOA problems, respectively. Section 7 is reserved for simulation results and discussions, and Section 8 concludes the paper with the implications of the results.

## 2. Formulating the Problem

This section presents the formulation of the problem, detailing both the optimal operation’s governing constraints and objective function.

### 2.1. Objective Function

Multiple objective functions are examined in this study, and they are as follows:

#### 2.1.1. The Minimization of Costs

The considered objective function includes the cost of WT ($C_{WT}$), the cost of PV units ($C_{PV}$), the annual price of energy loss ($C_{Loss}$), and the price of electricity procured from the grid ($C_{Grid}$). The total annual price can be written as:

$$C = C_{WT} + C_{PV} + C_{Loss} + C_{Grid} \quad (1)$$

In which,

$$C_{Grid} = 365 \times \sum_{h=1}^{24} P_{Grid(h)} \times U_{Grid(h)} \quad (2)$$

where $U_{Grid}$ is the price of buying energy from the network and $P_{Grid(h)}$ is the hourly power taken out of the grid.

$$C_{Loss} = 365 \times U_{Loss} \times \sum_{h=1}^{24} P_{T\_Loss(h)} \quad (3)$$

where $U_{Loss}$ represents the expenditure associated with energy loss and $P_{T\_Loss(h)}$ indicates the overall power losses during time $h$.

$$C_{PV} = C_{PV}^{inst} + C_{PV}^{O\&M} \quad (4)$$

where $C_{PV}^{O\&M}$ denotes the maintenance and operating expenses for the PV unit, while $C_{PV}^{inst}$ represents the installation cost of the PV system.

$$C_{PV}^{O\&M} = 365 \times U_{PV}^{O\&M} \times \sum_{i=1}^{N_s} \sum_{h=1}^{24} P_{PV(i,h)} \quad (5)$$
$$C_{PV}^{\text{inst.}} = CF \times U_{PV} \times P_{\text{rated\_PV}} \quad (6)$$

where $P_{\text{rated\_PV}}$ is the maximum amount of power that the PV system can generate, and CF is a factor that influences how quickly the cost of the PV system is paid back.

$$C_{WT} = C_{WT}^{\text{inst.}} + C_{WT}^{\text{O\&M}} \quad (7)$$

where $C_{WT}^{\text{O\&M}}$ is the wind’s maintenance and operating costs and $C_{WT}^{\text{inst.}}$ is the cost of installing the WT.

$$C_{WT}^{\text{O\&M}} = 365 \times U_{WT}^{\text{O\&M}} \times \sum_{i=1}^{Ns} \sum_{h=1}^{24} P_{WT(i,h)} \quad (8)$$

where the maintenance and operating expenses for WTs and PVs are represented by $U_{WT}^{\text{O\&M}}$ and $U_{PV}^{\text{O\&M}}$. $Ns$ is the number of RESs with the DSR system considered.

$$C_{WT}^{\text{inst.}} = CF \times U_{WT} \times P_{\text{rated\_WT}} \quad (9)$$

where $U_{WT}$ and $U_{PV}$ stand for the corresponding purchasing costs of WTs and PVs. The rated generated power of WT is denoted by $P_{\text{rated\_WT}}$. $P_{WT(i,h)}$ and $P_{PV(i,h)}$ denote the hourly produced power from WTs and PVs of the $i$th system.

$$CF = \frac{\beta \times (1 + \beta)^{NP}}{(1 + \beta)^{NP} - 1} \quad (10)$$

where $NP$ and $\beta$ are the system lifetime and interest rate of the PVs or WTs, respectively.

### 2.1.2. Increasing the Voltage Level

To ensure that the power grid operates efficiently and reliably, the voltage deviations should be maintained within an allowable range, ideally close to 1 p.u. The definition of the overall voltage deviation is [34,35]:

$$TVD = \sum_{h=1}^{24} \sum_{i=1}^{NB} |(V_{(i,h)} - 1)| \quad (11)$$

where $V_{(i,h)}$ represents the hourly voltage of the $i$th bus and $NB$ represents the number of buses in the network.

### 2.1.3. Enhanced System Stability

The index for voltage stability of the bus $j$ is as follows [36,37]:

$$VSI_j = |V_i|^4 - 4 (P_j X_{i,j} - Q_j R_{i,j})^2 - 4 (P_j R_{i,j} + Q_j X_{i,j})^2 |V_i|^2 \quad (12)$$

$$TVSI = \sum_{h=1}^{24} \sum_{j=2}^{NB} VSI_{(j,h)} \quad (13)$$

where $VSI_{(j,h)}$ is the voltage stability index for bus $j$. $P_j$ and $Q_j$ define the real and reactive power flow of the branch between $i$ and $j$ near at bus $j$, respectively. $X_{i,j}$ and $R_{i,j}$ represents the reactance and resistance of the transmission line.

The optimization takes into account the following three goal functions concurrently:

$$\min(F) = \min(\varepsilon_1 F_1 + \varepsilon_2 F_2 + \varepsilon_3 F_3) \quad (14)$$

$$F_1 = \frac{C_{\text{RERs/DSR}}}{C_{\text{Base}}} \quad (15)$$

$$F_2 = \frac{TVD_{\text{RERs/DSR}}}{TVD_{\text{Base}}} \quad (16)$$
$$C_{PV}^{\text{inst.}} = CF \times U_{PV} \times P_{\text{rated\_PV}} \qquad (6)$$

where $P_{\text{rated\_PV}}$ is the maximum amount of power that the PV system can generate, and CF is a factor that influences how quickly the cost of the PV system is paid back.

$$C_{WT} = C_{WT}^{\text{inst.}} + C_{WT}^{\text{O\&M}} \qquad (7)$$

where $C_{WT}^{\text{O\&M}}$ is the wind's maintenance and operating costs and $C_{WT}^{\text{inst.}}$ is the cost of installing the WT.

$$C_{WT}^{\text{O\&M}} = 365 \times U_{WT}^{\text{O\&M}} \times \sum_{i=1}^{Ns} \sum_{h=1}^{24} P_{WT(i,h)} \qquad (8)$$

where the maintenance and operating expenses for WTs and PVs are represented by $U_{WT}^{\text{O\&M}}$ and $U_{PV}^{\text{O\&M}}$. $Ns$ is the number of RESs with the DSR system considered.

$$C_{WT}^{\text{inst.}} = CF \times U_{WT} \times P_{\text{rated\_WT}} \qquad (9)$$

where $U_{WT}$ and $U_{PV}$ stand for the corresponding purchasing costs of WTs and PVs. The rated generated power of WT is denoted by $P_{\text{rated\_WT}}$. $P_{WT(i,h)}$ and $P_{PV(i,h)}$ denote the hourly produced power from WTs and PVs of the ith system.

$$CF = \frac{\beta \times (1 + \beta)^{NP}}{(1 + \beta)^{NP} - 1} \qquad (10)$$

where $NP$ and $\beta$ are the system lifetime and interest rate of the PVs or WTs, respectively.

### 2.1.2. Increasing the Voltage Level

To ensure that the power grid operates efficiently and reliably, the voltage deviations should be maintained within an allowable range, ideally close to 1 p.u. The definition of the overall voltage deviation is [34,35]:

$$TVD = \sum_{h=1}^{24} \sum_{i=1}^{NB} |(V_{(i,h)} - 1)| \qquad (11)$$

where $V_{(i,h)}$ represents the hourly voltage of the ith bus and $NB$ represents the number of buses in the network.

### 2.1.3. Enhanced System Stability

The index for voltage stability of the bus j is as follows [36,37]:

$$VSI_j = |V_i|^4 - 4 (P_j X_{i,j} - Q_j R_{i,j})^2 - 4 (P_j R_{i,j} + Q_j X_{i,j})^2 |V_i|^2 \qquad (12)$$

$$TVSI = \sum_{h=1}^{24} \sum_{j=2}^{NB} VSI_{(j,h)} \qquad (13)$$

where $VSI_{(j,h)}$ is the voltage stability index for bus j. $P_j$ and $Q_j$ define the real and reactive power flow of the branch between i and j near at bus j, respectively. $X_{i,j}$ and $R_{i,j}$ represents the reactance and resistance of the transmission line.

The optimization takes into account the following three goal functions concurrently:

$$\min(F) = \min(\varepsilon_1 F_1 + \varepsilon_2 F_2 + \varepsilon_3 F_3) \qquad (14)$$

$$F_1 = \frac{C_{RERs/DSR}}{C_{Base}} \qquad (15)$$

$$F_2 = \frac{TVD_{RERs/DSR}}{TVD_{Base}} \qquad (16)$$
$$F_3 = \frac{1}{TVSI} \qquad (17)$$

where Base and RERs/DSR are the base case and subscripts referring to RESs with DSR, respectively. The weighted values of $ε_1$, $ε_2$, and $ε_3$ were chosen to be 0.5, 0.25, and 0.25, respectively [15].

## 2.2. Inequality and Equality Constraints

### 2.2.1. Limitations of the Network (Inequality Constraints)

$$V_{min} \leq V_{(i,h)} \leq V_{max} \qquad (18)$$

$$\sum_{i=1}^{Ns} P_{PV\_rated,i} + \sum_{i=1}^{Ns} P_{WT\_rated,i} \leq \sum_{i=1}^{NB} P_{Load,i} \qquad (19)$$

$$PF_{min} \leq PF_{i} \leq PF_{max} \quad i = 1, 2, \dots, Ns \qquad (20)$$

$$I_i \leq I_{max,i} \quad i = 1, 2, \dots, NT \qquad (21)$$

where $V_{max}$ and $V_{min}$ are the upper and lower voltage limits, respectively. $P_{Load}$ represents the real load. $NT$ defines the number of lines. $I_{max,i}$ is the maximum allowable current limit of the line $i$. $PF_{min}$ and $PF_{max}$ are the minimum and maximum of the WT power factor, respectively. The power factor of the PV farms should be equal to one.

### 2.2.2. Equality Constraints

$$P_S(h) + \sum_{i=1}^{Ns} P_{PV,i}(h) + \sum_{i=1}^{Ns} P_{WT,i}(h) = \sum_{i=1}^{NT} P_{Loss,i}(h) + \sum_{i=1}^{NB} P_{Load,i}(h) \qquad (22)$$

$$Q_S(h) + \sum_{i=1}^{Ns} Q_{WT,i}(h) = \sum_{i=1}^{NT} Q_{Loss,i}(h) + \sum_{i=1}^{NB} Q_{Load,i}(h) \qquad (23)$$

where $Q_S(h)$ and $P_S(h)$ are the hourly reactive and real powers of the main network, respectively.

## 3. Modeling the Systems

### 3.1. PV System

The power output of a PV system can be determined using the following formula:

$$T_c(h) = T_a(h) + \frac{I(h)}{800} \cdot (NOCT - 20) \qquad (24)$$

$$P_{PV}(h) = A_{PV} \cdot \eta_{PV}(h) \cdot I(h) \qquad (25)$$

The total area used by the PV array, labeled as $A_{PV}$ in (m²), is multiplied by a number representing the efficiency of conversion of the PV panels, $η_{PV}$, to get the output produced by the PV, expressed as $P_{PV}$ (kW). The temperature of the cell is expressed in degrees Celsius as $T_c$. The temperature outside in degrees Celsius is represented by $T_a$. $I$ is the solar insolation in (kW/m²). The instantaneous efficiency of PV panels $η_{PV}$ is obtained using the following equation [38]:

$$\eta_{PV}(h) = \eta_r \cdot \eta_t \times \left[ 1 - \gamma \cdot (T_a(h) - T_r) - \gamma \cdot I(h) \cdot \left( \frac{NOCT-20}{800} \right) \cdot (1 - \eta_r \cdot \eta_t) \right] \qquad (26)$$

where $η_t$ is the effectiveness of the maximum power point tracking device, and $η_r$ is the PV panels' reference efficiency. It also takes into account the temperature coefficient of efficiency $γ$, which for silicon cells, normally ranges from 0.004 to 0.006 per (°C). Notable operating cell temperature NOCT in (°C), reference temperature $T_r$ in (°C) for PV cells, and ambient temperature $T_a$ in (°C) are also considered.
$$F_3 = \frac{1}{TVSI} \qquad (17)$$

where Base and RERs/DSR are the base case and subscripts referring to RESs with DSR, respectively. The weighted values of $ε_1$, $ε_2$, and $ε_3$ were chosen to be 0.5, 0.25, and 0.25, respectively [15].

## 2.2. Inequality and Equality Constraints

### 2.2.1. Limitations of the Network (Inequality Constraints)

$$V_{min} \leq V_{(i,h)} \leq V_{max} \qquad (18)$$

$$\sum_{i=1}^{Ns} P_{PV\_rated,i} + \sum_{i=1}^{Ns} P_{WT\_rated,i} \leq \sum_{i=1}^{NB} P_{Load,i} \qquad (19)$$

$$PF_{min} \leq PF_{i} \leq PF_{max} \quad i = 1, 2, \dots, Ns \qquad (20)$$

$$I_i \leq I_{max,i} \quad i = 1, 2, \dots, NT \qquad (21)$$

where $V_{max}$ and $V_{min}$ are the upper and lower voltage limits, respectively. $P_{Load}$ represents the real load. $NT$ defines the number of lines. $I_{max,i}$ is the maximum allowable current limit of the line $i$. $PF_{min}$ and $PF_{max}$ are the minimum and maximum of the WT power factor, respectively. The power factor of the PV farms should be equal to one.

### 2.2.2. Equality Constraints

$$P_S(h) + \sum_{i=1}^{Ns} P_{PV,i}(h) + \sum_{i=1}^{Ns} P_{WT,i}(h) = \sum_{i=1}^{NT} P_{Loss,i}(h) + \sum_{i=1}^{NB} P_{Load,i}(h) \qquad (22)$$

$$Q_S(h) + \sum_{i=1}^{Ns} Q_{WT,i}(h) = \sum_{i=1}^{NT} Q_{Loss,i}(h) + \sum_{i=1}^{NB} Q_{Load,i}(h) \qquad (23)$$

where $Q_S(h)$ and $P_S(h)$ are the hourly reactive and real powers of the main network, respectively.

## 3. Modeling the Systems

### 3.1. PV System

The power output of a PV system can be determined using the following formula:

$$T_c(h) = T_a(h) + \frac{I(h)}{800} \cdot (NOCT - 20) \qquad (24)$$

$$P_{PV}(h) = A_{PV} \cdot \eta_{PV}(h) \cdot I(h) \qquad (25)$$

The total area used by the PV array, labeled as $A_{PV}$ in (m²), is multiplied by a number representing the efficiency of conversion of the PV panels, $η_{PV}$, to get the output produced by the PV, expressed as $P_{PV}$ (kW). The temperature of the cell is expressed in degrees Celsius as $T_c$. The temperature outside in degrees Celsius is represented by $T_a$. $I$ is the solar insolation in (kW/m²). The instantaneous efficiency of PV panels $η_{PV}$ is obtained using the following equation [38]:

$$\eta_{PV}(h) = \eta_r \cdot \eta_t \times \left[1 - \gamma \cdot (T_a(h) - T_r) - \gamma \cdot I(h) \cdot \left(\frac{NOCT-20}{800}\right) \cdot (1 - \eta_r \cdot \eta_t)\right] \qquad (26)$$

where $η_t$ is the effectiveness of the maximum power point tracking device, and $η_r$ is the PV panels' reference efficiency. It also takes into account the temperature coefficient of efficiency $γ$, which for silicon cells, normally ranges from 0.004 to 0.006 per (°C). Notable operating cell temperature NOCT in (°C), reference temperature $T_r$ in (°C) for PV cells, and ambient temperature $T_a$ in (°C) are also considered.
## 3.2. WT System

The WT energy output ($P_{WT}$) can be determined using the following formula:

$$P_{WT}(W) = \begin{cases} 0 & \text{for } W < W_i \text{ and } W > W_o \\ P_{rated\_WT}\left(\frac{W-W_i}{W_r-W_i}\right) & \text{for } (W_i \le W \le W_r) \\ P_{rated\_WT} & \text{for } (W_r < W \le W_o) \end{cases} \quad (27)$$

The used wind turbine has a rated power of 250 kW, a rated velocity $W_r$ of 15 m/s, a cut-out speed $W_o$ of 25 m/s, and a cut-in speed $W_i$ of 2.5 m/s [39].

## 3.3. Demand Side Response

The power fluctuations in the distribution system can vary significantly on a daily basis, resulting in substantial differences between the minimum and maximum energy consumption. Demand Side Response is a viable approach to modifying customers' consumption patterns. In essence, DSR entails altering the power consumption of electric utility customers to align it more effectively with the system's supply. In other words, demand response programs focus on adjusting the demand for power rather than manipulating the supply [40,41].

In this research, the real-time pricing (RTP) program is employed as the demand-side response program due to its widespread and practical application. RTP is considered one of the most common and beneficial DSR programs, and its implementation in a real distribution system is relatively straightforward. As a result, the outcomes of this study closely align with the performance of an actual distribution system [42,43].

# 4. Representing the Uncertainties

The uncertainty in the parameters is depicted in the following manner:

## 4.1. A Probabilistic Representation of Solar Irradiance

The variations in sun irradiance have been modeled using the Beta PDF in the following way [44,45]:

$$f_b(G) = \begin{cases} \frac{\Gamma(\varphi+\tau)}{\Gamma(\varphi)\Gamma(\tau)} s^{(\varphi-1)} (1-G)^{(\tau-1)} & 0 \le G \le 1; \varphi, \tau \ge 0 \\ 0 & \text{otherwise} \end{cases} \quad (28)$$

where $\mu$ represents the mean value derived from past data and $\sigma$ is the standard deviation. The following equations can be used to determine $\tau$ and $\varphi$ [46,47]:

$$\tau = (1 - \mu) \times \left( \frac{\mu \times (1 + \mu)}{\sigma^2} - 1 \right) \quad (29)$$

$$\varphi = \frac{\mu \times \tau}{1 - \mu} \quad (30)$$

## 4.2. A Probabilistic Representation of Wind Speed

To represent the uncertainty of wind speed, the following is an explanation of how the Weibull PDF is applied [48,49]:

$$F(W) = \left(\frac{k}{c}\right) \left(\frac{W}{c}\right)^{k-1} \exp \left[ - \left(\frac{W}{c}\right)^k \right] \quad (31)$$

where W stands for the speed of the wind. The Weibull PDF scale and parameters for shape are represented by the c and k, respectively.
## 3.2. WT System

The WT energy output ($P_{WT}$) can be determined using the following formula:

$$ P_{WT}(W) = \begin{cases} 0 & \text{for } W < W_i \text{ and } W > W_o \\ P_{rated\_WT}\left(\frac{W-W_i}{W_r-W_i}\right) & \text{for } (W_i \le W \le W_r) \\ P_{rated\_WT} & \text{for } (W_r < W \le W_o) \end{cases} \quad (27) $$

The used wind turbine has a rated power of 250 kW, a rated velocity $W_r$ of 15 m/s, a cut-out speed $W_o$ of 25 m/s, and a cut-in speed $W_i$ of 2.5 m/s [39].

## 3.3. Demand Side Response

The power fluctuations in the distribution system can vary significantly on a daily basis, resulting in substantial differences between the minimum and maximum energy consumption. Demand Side Response is a viable approach to modifying customers' consumption patterns. In essence, DSR entails altering the power consumption of electric utility customers to align it more effectively with the system's supply. In other words, demand response programs focus on adjusting the demand for power rather than manipulating the supply [40,41].

In this research, the real-time pricing (RTP) program is employed as the demand-side response program due to its widespread and practical application. RTP is considered one of the most common and beneficial DSR programs, and its implementation in a real distribution system is relatively straightforward. As a result, the outcomes of this study closely align with the performance of an actual distribution system [42,43].

# 4. Representing the Uncertainties

The uncertainty in the parameters is depicted in the following manner:

## 4.1. A Probabilistic Representation of Solar Irradiance

The variations in sun irradiance have been modeled using the Beta PDF in the following way [44,45]:

$$ f_b(G) = \begin{cases} \frac{\Gamma(\varphi+\tau)}{\Gamma(\varphi)\Gamma(\tau)} s^{(\varphi-1)} (1-G)^{(\tau-1)} & 0 \le G \le 1; \varphi, \tau \ge 0 \\ 0 & \text{otherwise} \end{cases} \quad (28) $$

where $\mu$ represents the mean value derived from past data and $\sigma$ is the standard deviation. The following equations can be used to determine $\tau$ and $\varphi$ [46,47]:

$$ \tau = (1 - \mu) \times \left( \frac{\mu \times (1 + \mu)}{\sigma^2} - 1 \right) \quad (29) $$

$$ \varphi = \frac{\mu \times \tau}{1 - \mu} \quad (30) $$

## 4.2. A Probabilistic Representation of Wind Speed

To represent the uncertainty of wind speed, the following is an explanation of how the Weibull PDF is applied [48,49]:

$$ F(W) = \left(\frac{k}{c}\right) \left(\frac{W}{c}\right)^{k-1} \exp \left[ - \left(\frac{W}{c}\right)^k \right] \quad (31) $$

where W stands for the speed of the wind. The Weibull PDF scale and parameters for shape are represented by the c and k, respectively.
### 4.3. A Probabilistic Representation of Load Demand

The normal PDF can be used to model the uncertainty in load demand as follows. [47]:

$$f_n(L) = \frac{1}{\sigma_L \sqrt{2\pi}} \times \exp \left[ -\frac{(L - \mu_L)^2}{2\sigma_L^2} \right] \quad (32)$$

where the loading’s mean value is represented by $\mu_L$ and the standard deviation is represented by $\sigma_L$.

### 4.4. The Probabilistic Representation of Price

The price of electricity is a highly important random characteristic in a power system since it is an uncertain value derived from the grid. Equation (33) illustrates how the normal PDF, based on its mean value $\mu_P$ and standard deviation $\sigma_P$, can be used to simulate the PDF of the price of electricity [50,51]:

$$f(P) = \frac{1}{\sigma_P \sqrt{2\pi}} \exp \left[ -\frac{(P - \mu_P)^2}{2\sigma_P^2} \right] \quad (33)$$

### 4.5. The Probabilistic Representation of Temperature

Since the temperature outside is constantly changing, it is considered an erratic variable, and its uncertainty is expressed as follows using the normal probability distribution:

$$f(T) = \frac{1}{\sigma_T \sqrt{2\pi}} \exp \left[ -\frac{(T - \mu_T)^2}{2\sigma_T^2} \right] \quad (34)$$

where the mean temperature is represented by $\mu_T$ and the standard deviation by $\sigma_T$.

Monte Carlo simulations were used to generate 800 data points for solar irradiance, temperature, wind speed, load, and price, using the parameters of the probability density function (PDF) of each variable determined from the collected data [52]. However, to reduce the computational burden of the simulation, a data-based reduction method followed by conditional control was used to reduce the number of cases to a more manageable size (12 cases) [53,54]. Figure 2 shows generated scenarios of these parameters at 12:00. A.M.
### 4.3. A Probabilistic Representation of Load Demand

The normal PDF can be used to model the uncertainty in load demand as follows. [47]:

$$f_n(L) = \frac{1}{\sigma_L \sqrt{2\pi}} \times \exp \left[ -\frac{(L - \mu_L)^2}{2\sigma_L^2} \right] \quad (32)$$

where the loading’s mean value is represented by $\mu_L$ and the standard deviation is represented by $\sigma_L$.

### 4.4. The Probabilistic Representation of Price

The price of electricity is a highly important random characteristic in a power system since it is an uncertain value derived from the grid. Equation (33) illustrates how the normal PDF, based on its mean value $\mu_P$ and standard deviation $\sigma_P$, can be used to simulate the PDF of the price of electricity [50,51]:

$$f(P) = \frac{1}{\sigma_P \sqrt{2\pi}} \exp \left[ -\frac{(P - \mu_P)^2}{2\sigma_P^2} \right] \quad (33)$$

### 4.5. The Probabilistic Representation of Temperature

Since the temperature outside is constantly changing, it is considered an erratic variable, and its uncertainty is expressed as follows using the normal probability distribution:

$$f(T) = \frac{1}{\sigma_T \sqrt{2\pi}} \exp \left[ -\frac{(T - \mu_T)^2}{2\sigma_T^2} \right] \quad (34)$$

where the mean temperature is represented by $\mu_T$ and the standard deviation by $\sigma_T$.

Monte Carlo simulations were used to generate 800 data points for solar irradiance, temperature, wind speed, load, and price, using the parameters of the probability density function (PDF) of each variable determined from the collected data [52]. However, to reduce the computational burden of the simulation, a data-based reduction method followed by conditional control was used to reduce the number of cases to a more manageable size (12 cases) [53,54]. Figure 2 shows generated scenarios of these parameters at 12:00. A.M.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** The scenarios generated by MCS at 12:00 A.M. for (a) loading, (b) temperature, (c) solar irradiance, (d) price, and (e) wind speed.</figcaption></figure>

Figure 3 illustrates specific scenarios selected for Monte Carlo simulations of a power system over a 24 h period after data reduction. The figure provides a unique insight into the probabilistic analysis of the system, where the data have been refined to enhance simulation efficiency and present an accurate depiction of potential scenarios throughout the specified period.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** The scenarios generated by MCS at 12:00 A.M. for (a) loading, (b) temperature, (c) solar irradiance, (d) price, and (e) wind speed.</figcaption></figure>

Figure 3 illustrates specific scenarios selected for Monte Carlo simulations of a power system over a 24 h period after data reduction. The figure provides a unique insight into the probabilistic analysis of the system, where the data have been refined to enhance simulation efficiency and present an accurate depiction of potential scenarios throughout the specified period.
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Data-reduced scenarios for Monte Carlo simulations of a power system over 24 h for (a) loading, (b) temperature, (c) solar irradiance, (d) price, and (e) wind speed.</figcaption></figure>

**5. Walrus Optimization Algorithm (WaOA)**

WaOA is a bio-inspired algorithm that uses the behavior of walruses to search for food. By randomly initializing the population of walruses and then moving them around according to their fitness values, WaOA can effectively explore the search space and find better solutions [55].
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Data-reduced scenarios for Monte Carlo simulations of a power system over 24 h for (a) loading, (b) temperature, (c) solar irradiance, (d) price, and (e) wind speed.</figcaption></figure>

**5. Walrus Optimization Algorithm (WaOA)**

WaOA is a bio-inspired algorithm that uses the behavior of walruses to search for food. By randomly initializing the population of walruses and then moving them around according to their fitness values, WaOA can effectively explore the search space and find better solutions [55].
## 5.1. Feeding Strategy

The position of each walrus is updated by moving towards the strongest walrus, which is the best solution found so far. The mathematical representation of updating the walruses' positions is structured around their feeding behavior, guided by the dominant member of the group, as described by Equations (35) and (36).

$$x_{i,j}^{P_1} = x_{i,j} + rand_{i,j} \cdot (SW_j - I_{i,j} \cdot x_{i,j}) \quad (35)$$

$$X_i = \begin{cases} X_i^{P_1}, & F_i^{P_1} < F_i, \\ X_i, & \text{else,} \end{cases} \quad (36)$$

In the first phase of WaOA, new positions $x_{i,j}^{P_1}$ for each walrus $i^{th}$ are generated based on the feeding strategy. The $j^{th}$ dimension of the new position is $x_{i,j}^{P_1}$. After evaluating the function of objective $F_i^{P_1}$ at this new location, $x_{i,j}$ represents the location of particle $i$ in dimension $j$. Meanwhile, $X_i$ is a vector that encompasses the position across all dimensions. $rand_{i,j}$ introduces randomness by generating numbers between 0 and 1. The best solution so far $SW_j$ acts as the strongest walrus that guides movement. To boost exploration, $I_{i,j}$ are integers of 1 or 2 chosen randomly. $I_{i,j} = 2$ causes broader jumps versus 1 for normal steps. This expands the search space, helping to escape local optima and find new promising areas. In summary, the first phase updates positions via a randomized process biased towards the best solution so far. Occasional large steps enhance global exploration to avoid getting trapped.

## 5.2. Migration

The position of each walrus is updated by migrating toward the position of another randomly selected walrus. WaOA utilizes this migration procedure to steer the walruses throughout the exploration of the search space, aiding them in identifying favorable regions within it. This behavioral mechanism is formally expressed through Equations (37) and (38).

$$x_{i,j}^{P_2} = \begin{cases} x_{i,j} + rand_{i,j} \cdot (x_{k,j} - I_{i,j} \cdot x_{i,j}), & F_k < F_i, \\ x_{i,j} + rand_{i,j} \cdot (x_{i,j} - x_{k,j}), & \text{else,} \end{cases} \quad (37)$$

$$X_i = \begin{cases} X_i^{P_2}, & F_i^{P_2} < F_i, \\ X_i, & \text{else,} \end{cases} \quad (38)$$

Here, $x_{i,j}^{P_2}$ represents the newly generated position of the $i^{th}$ walrus derived from the second phase. $x_{i,j}^{P_2}$ signifies its position on the $j^{th}$ dimension. $F_i^{P_2}$ corresponds to its objective function value. $X_k$, where $k \in \{1, 2, ..., N\}$ and $k \neq i$, represents the position of the selected walrus to which the $i^{th}$ walrus is migrated. $x_{k,j}$ is the position of the $k^{th}$ walrus on the $j^{th}$ dimension, and $F_k$ indicates its objective function value.

## 5.3. Escaping Predators

WaOA uses the natural behavior of walruses to improve its ability to find food. By simulating the way walruses shift their positions in response to threats, WaOA can more effectively explore the search space and find better solutions. In the algorithm, a virtual neighborhood is established around each walrus, where new positions are randomly generated within this range using Equations (39) and (40). If the objective function improves with this novel location, it substitutes the old one, as outlined in Equation (41). This behavior simulation within WaOA enhances its efficacy in local search and adaptation to the problem's landscape.

$$x_{i,j}^{P_3} = x_{i,j} + (lb_{local,j}^{t} + (ub_{local,j}^{t} - rand \cdot lb_{local,j}^{t})) \quad (39)$$
## 5.1. Feeding Strategy

The position of each walrus is updated by moving towards the strongest walrus, which is the best solution found so far. The mathematical representation of updating the walruses' positions is structured around their feeding behavior, guided by the dominant member of the group, as described by Equations (35) and (36).

$$x_{i,j}^{P_1} = x_{i,j} + rand_{i,j} \cdot (SW_j - I_{i,j} \cdot x_{i,j}) \quad (35)$$

$$X_i = \begin{cases} X_i^{P_1}, & F_i^{P_1} < F_i, \\ X_i, & \text{else,} \end{cases} \quad (36)$$

In the first phase of WaOA, new positions $x_{i,j}^{P_1}$ for each walrus $i^{th}$ are generated based on the feeding strategy. The $j^{th}$ dimension of the new position is $x_{i,j}^{P_1}$. After evaluating the function of objective $F_i^{P_1}$ at this new location, $x_{i,j}$ represents the location of particle $i$ in dimension $j$. Meanwhile, $X_i$ is a vector that encompasses the position across all dimensions. $rand_{i,j}$ introduces randomness by generating numbers between 0 and 1. The best solution so far $SW_j$ acts as the strongest walrus that guides movement. To boost exploration, $I_{i,j}$ are integers of 1 or 2 chosen randomly. $I_{i,j} = 2$ causes broader jumps versus 1 for normal steps. This expands the search space, helping to escape local optima and find new promising areas. In summary, the first phase updates positions via a randomized process biased towards the best solution so far. Occasional large steps enhance global exploration to avoid getting trapped.

## 5.2. Migration

The position of each walrus is updated by migrating toward the position of another randomly selected walrus. WaOA utilizes this migration procedure to steer the walruses throughout the exploration of the search space, aiding them in identifying favorable regions within it. This behavioral mechanism is formally expressed through Equations (37) and (38).

$$x_{i,j}^{P_2} = \begin{cases} x_{i,j} + rand_{i,j} \cdot (x_{k,j} - I_{i,j} \cdot x_{i,j}), & F_k < F_i, \\ x_{i,j} + rand_{i,j} \cdot (x_{i,j} - x_{k,j}), & \text{else,} \end{cases} \quad (37)$$

$$X_i = \begin{cases} X_i^{P_2}, & F_i^{P_2} < F_i, \\ X_i, & \text{else,} \end{cases} \quad (38)$$

Here, $x_{i,j}^{P_2}$ represents the newly generated position of the $i^{th}$ walrus derived from the second phase. $x_{i,j}^{P_2}$ signifies its position on the $j^{th}$ dimension. $F_i^{P_2}$ corresponds to its objective function value. $X_k$, where $k \in \{1, 2, \dots, N\}$ and $k \neq i$, represents the position of the selected walrus to which the $i^{th}$ walrus is migrated. $x_{k,j}$ is the position of the $k^{th}$ walrus on the $j^{th}$ dimension, and $F_k$ indicates its objective function value.

## 5.3. Escaping Predators

WaOA uses the natural behavior of walruses to improve its ability to find food. By simulating the way walruses shift their positions in response to threats, WaOA can more effectively explore the search space and find better solutions. In the algorithm, a virtual neighborhood is established around each walrus, where new positions are randomly generated within this range using Equations (39) and (40). If the objective function improves with this novel location, it substitutes the old one, as outlined in Equation (41). This behavior simulation within WaOA enhances its efficacy in local search and adaptation to the problem's landscape.

$$x_{i,j}^{P_3} = x_{i,j} + (lb_{local,j}^{t} + (ub_{local,j}^{t} - rand \cdot lb_{local,j}^{t})) \quad (39)$$
$$ \text{Local bounds :} \quad \left\{ \begin{array}{l} lb_{\text{local},j}^{t} = \frac{lb_j}{t}, \\ ub_{\text{local},j}^{t} = \frac{ub_j}{t}, \end{array} \right. \qquad (40) $$

$$ X_i = \begin{cases} X_i^{P_3}, & F_i^{P_3} < F_i, \\ X_i, & \text{else,} \end{cases} \qquad (41) $$

In this context, $x_{i,j}^{P_3}$ refers to the newly created position of the $i^{th}$ walrus using the third phase. $x_{i,j}^{P_3}$ designates its location on the $j^{th}$ dimension. $F_i^{P_3}$ stands for its objective function value. $t$ represents the iteration cycle. $lb_j$ and $ub_j$ symbolize the lower and upper limits of the $j^{th}$ variable, respectively. Additionally, $lb_{\text{local},j}^{t}$ and $ub_{\text{local},j}^{t}$ denote the locally permitted lower and upper bounds for the $j^{th}$ variable, respectively, which simulate localized exploration within the vicinity of potential solutions.

# 6. Improved Walrus Optimization Algorithm (I-WaOA)

I-WaOA employs three strategies. The initial approach uses the Cauchy mutation, followed by the second method grounded in the fitness-distance balance (FDB) concept. The final strategy is the Quasi-Oppositional-Based Learning (QOBL) method.

## 6.1. The Cauchy Mutation (CM)

CM is a powerful technique that can be used to improve the efficiency of algorithm optimization. It has been shown to be effective in escaping from local optima and finding better solutions to problems [29,56,57]. In the proposed algorithm, CM is used to generate new solutions that are more likely to be better than the current solutions. The mathematical representation of CM is as follows:

$$ X_i^{t+1} = X_i^t + \alpha \otimes \text{Cauchy}(\mu, \sigma) \otimes (X_i^t - X_{best}) \qquad (42) $$

$$ F(x) = \frac{1}{\pi} \arctan \left( \frac{2(x - \mu)}{\sigma} \right) + \frac{1}{2} \qquad (43) $$

where $\otimes$ represent entrywise multiplication, $\alpha$ denotes to the step size, Cauchy $(\mu, \sigma)$ represent Cauchy distribution with mean $\mu$ and standard deviation $\sigma$, $F(x)$ is cumulative distribution function of the Cauchy distribution function, and $\sigma$ and $\mu$ are the scale parameters and the location, with values selected to be 4.5 and 0.8, respectively [58].

## 6.2. The Fitness-Distance Balance (FDB)

One way to choose the ideal optimization solution from the pool of available options is to use the FDB selection method. It works by measuring the value of the strength of each candidate solution and its distance from the optimal solution. The difference between two solutions is calculated by taking the absolute value. The FDB method then selects candidate solutions with their highest fitness values and shortest distances to the optimal solution. This helps the algorithm focus on candidate solutions that are good and search the search area around the optimal solution [59–62].

$$ DS_i = \sqrt{(x_i^1 - \text{Best } 1)^2 + (x_i^2 - \text{Best } 2)^2 + \dots + (x_i^d - \text{Best } d)^2} \qquad (44) $$

Then, build the fitness and distance value vectors in the following manner:

$$ DS = [DS_1, DS_2, \dots, DS_n] \qquad (45) $$

$$ F = [F_1, F_2, \dots, F_n] \qquad (46) $$
$$ \text{Local bounds :} \quad \left\{ \begin{array}{l} lb_{local,j}^{t} = \frac{lb_{i}}{t}, \\ ub_{local,j}^{t} = \frac{ub_{j}}{t}, \end{array} \right. \qquad (40) $$

$$ X_i = \begin{cases} X_i^{P_3}, & F_i^{P_3} < F_i, \\ X_i, & \text{else,} \end{cases} \qquad (41) $$

In this context, $x_{i,j}^{P_3}$ refers to the newly created position of the $i^{th}$ walrus using the third phase. $x_{i,j}^{P_3}$ designates its location on the $j^{th}$ dimension. $F_i^{P_3}$ stands for its objective function value. $t$ represents the iteration cycle. $lb_j$ and $ub_j$ symbolize the lower and upper limits of the $j^{th}$ variable, respectively. Additionally, $lb_{local,j}^{t}$ and $ub_{local,j}^{t}$ denote the locally permitted lower and upper bounds for the $j^{th}$ variable, respectively, which simulate localized exploration within the vicinity of potential solutions.

## 6. Improved Walrus Optimization Algorithm (I-WaOA)

I-WaOA employs three strategies. The initial approach uses the Cauchy mutation, followed by the second method grounded in the fitness-distance balance (FDB) concept. The final strategy is the Quasi-Oppositional-Based Learning (QOBL) method.

### 6.1. The Cauchy Mutation (CM)

CM is a powerful technique that can be used to improve the efficiency of algorithm optimization. It has been shown to be effective in escaping from local optima and finding better solutions to problems [29,56,57]. In the proposed algorithm, CM is used to generate new solutions that are more likely to be better than the current solutions. The mathematical representation of CM is as follows:

$$ X_i^{t+1} = X_i^t + \alpha \otimes \text{Cauchy}(\mu, \sigma) \otimes (X_i^t - X_{best}) \qquad (42) $$

$$ F(x) = \frac{1}{\pi} \arctan \left( \frac{2(x - \mu)}{\sigma} \right) + \frac{1}{2} \qquad (43) $$

where $\otimes$ represent entrywise multiplication, $\alpha$ denotes to the step size, Cauchy $(\mu, \sigma)$ represent Cauchy distribution with mean $\mu$ and standard deviation $\sigma$, $F(x)$ is cumulative distribution function of the Cauchy distribution function, and $\sigma$ and $\mu$ are the scale parameters and the location, with values selected to be 4.5 and 0.8, respectively [58].

### 6.2. The Fitness-Distance Balance (FDB)

One way to choose the ideal optimization solution from the pool of available options is to use the FDB selection method. It works by measuring the value of the strength of each candidate solution and its distance from the optimal solution. The difference between two solutions is calculated by taking the absolute value. The FDB method then selects candidate solutions with their highest fitness values and shortest distances to the optimal solution. This helps the algorithm focus on candidate solutions that are good and search the search area around the optimal solution [59–62].

$$ DS_i = \sqrt{(x_i^1 - \text{Best } 1)^2 + (x_i^2 - \text{Best } 2)^2 + \dots + (x_i^d - \text{Best } d)^2} \qquad (44) $$

Then, build the fitness and distance value vectors in the following manner:

$$ DS = [DS_1, DS_2, \dots, DS_n] \qquad (45) $$

$$ F = [F_1, F_2, \dots, F_n] \qquad (46) $$
After that, the fitness and distance values may be adjusted in the manner shown below:

$$normDS_i = \frac{DS_i - min(DS)}{max(DS) - min(DS)} \qquad (47)$$

$$normF_i = \frac{F_i - min(F)}{max(F) - min(F)} \qquad (48)$$

where the fitness vectors and distance's minimum and maximum are represented by the terms *min* and *max*, respectively. The FDB score is determined using the formula below:

$$FDBscore_i = \alpha * (1 - normF_i) + (1 - \alpha) * normDS_i \qquad (49)$$

In which

$$\alpha = 0.5 * \left( 1 + \frac{t}{T_{max}} \right) \qquad (50)$$

## 6.3. The Quasi-Opposite Based Learning (QOBL)

QOBL is a highly effective method designed to identify specific regions within a population that are likely to yield optimal solutions. By strategically employing the QOBL selection and promptly assigning it a collision count, the algorithm expedites the process, leading to improved efficiency in the context of meta-heuristic algorithms [63–68]. The following is an illustration of QOBL.

for i = 1 : No. Populations

for j = 1 : No. Dimensions

$$x_{i,j}^{o} = Lb_{j} + Ub_{j} - x_{i,j}$$

$$C_{i,j} = (Lb_{j} + Ub_{j}) / 2$$

if $x_{i,j} < C_{i,j}$

$$x_{i,j}^{qo} = C_{i,j} + (x_{i,j}^{o} - C_{i,j}) \times \text{rand}$$

else

$$x_{i,j}^{qo} = C_{i,j} + (C_{i,j} - x_{i,j}^{o}) \times \text{rand}$$

end

end

end

The Improved Walrus Optimization Algorithm represents an enhanced iteration of the WaOA algorithm, integrating three key strategies to elevate its performance. The utilization of Cauchy mutation, a more assertive operator than in the original WaOA, enables I-WaOA to break free from local optima, fostering exploration across new regions of the search space. The incorporation of Fitness-Distance Balance ensures a harmonized equilibrium between exploration and exploitation capabilities, preventing entrapment in local optima and facilitating the discovery of the global optimal solution. Additionally, the Quasi-Oppositional-Based Learning method empowers I-WaOA to navigate uncharted territories within the search space, contributing to its proficiency in locating the global optimal solution. The I-WaOA algorithm was developed to examine the weaknesses of the basic WaOA algorithm, such as its susceptibility to getting trapped in local optima and its potential for performance improvement.

The I-WaOA for the optimal operation solution is shown in Figure 4. Figure 5 shows the process for solving the stochastic OOP of a DN.
After that, the fitness and distance values may be adjusted in the manner shown below:

$$normDS_i = \frac{DS_i - min(DS)}{max(DS) - min(DS)} \qquad (47)$$

$$normF_i = \frac{F_i - min(F)}{max(F) - min(F)} \qquad (48)$$

where the fitness vectors and distance's minimum and maximum are represented by the terms *min* and *max*, respectively. The FDB score is determined using the formula below:

$$FDBscore_i = \alpha * (1 - normF_i) + (1 - \alpha) * normDS_i \qquad (49)$$

In which

$$\alpha = 0.5 * \left( 1 + \frac{t}{T_{max}} \right) \qquad (50)$$

## 6.3. The Quasi-Opposite Based Learning (QOBL)

QOBL is a highly effective method designed to identify specific regions within a population that are likely to yield optimal solutions. By strategically employing the QOBL selection and promptly assigning it a collision count, the algorithm expedites the process, leading to improved efficiency in the context of meta-heuristic algorithms [63–68]. The following is an illustration of QOBL.

for *i* = 1 : No. Populations

for *j* = 1 : No. Dimensions

$$x_{i,j}^{o} = Lb_{j} + Ub_{j} - x_{i,j}$$

$$C_{i,j} = (Lb_{j} + Ub_{j}) / 2$$

if $x_{i,j} < C_{i,j}$

$$x_{i,j}^{qo} = C_{i,j} + (x_{i,j}^{o} - C_{i,j}) \times \text{rand}$$

else

$$x_{i,j}^{qo} = C_{i,j} + (C_{i,j} - x_{i,j}^{o}) \times \text{rand}$$

end

end

end

The Improved Walrus Optimization Algorithm represents an enhanced iteration of the WaOA algorithm, integrating three key strategies to elevate its performance. The utilization of Cauchy mutation, a more assertive operator than in the original WaOA, enables I-WaOA to break free from local optima, fostering exploration across new regions of the search space. The incorporation of Fitness-Distance Balance ensures a harmonized equilibrium between exploration and exploitation capabilities, preventing entrapment in local optima and facilitating the discovery of the global optimal solution. Additionally, the Quasi-Oppositional-Based Learning method empowers I-WaOA to navigate uncharted territories within the search space, contributing to its proficiency in locating the global optimal solution. The I-WaOA algorithm was developed to examine the weaknesses of the basic WaOA algorithm, such as its susceptibility to getting trapped in local optima and its potential for performance improvement.

The I-WaOA for the optimal operation solution is shown in Figure 4. Figure 5 shows the process for solving the stochastic OOP of a DN.
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Flow diagram for the optimal operation solution for the suggested I-WaOA.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Steps for addressing the integration of RESs with DSR and system uncertainties for determining the optimal operation of a distribution grid.</figcaption></figure>

## 7. Results of Simulation

The performance of the proposed Improved Walrus Optimization Algorithm (I-WaOA) for solving the IEEE 118-bus system optimization problem was evaluated by comparing it against several established optimization algorithms. These included Chernobyl disaster optimizer (CDO) [69], Artificial hummingbird algorithm (AHA) [70], Dandelion Optimizer (DO) [71], Harmony Search (HS) [72], Artificial rabbits optimization (ARO) [73], Sand Cat swarm optimization (SCSO) [74], Zebra optimization algorithm (ZOA) [75], and standard walruses behavior (WaOA) [55]. All the algorithms were executed using MATLAB R2019b.
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Flow diagram for the optimal operation solution for the suggested I-WaOA.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Steps for addressing the integration of RESs with DSR and system uncertainties for determining the optimal operation of a distribution grid.</figcaption></figure>

## 7. Results of Simulation

The performance of the proposed Improved Walrus Optimization Algorithm (I-WaOA) for solving the IEEE 118-bus system optimization problem was evaluated by comparing it against several established optimization algorithms. These included Chernobyl disaster optimizer (CDO) [69], Artificial hummingbird algorithm (AHA) [70], Dandelion Optimizer (DO) [71], Harmony Search (HS) [72], Artificial rabbits optimization (ARO) [73], Sand Cat swarm optimization (SCSO) [74], Zebra optimization algorithm (ZOA) [75], and standard walruses behavior (WaOA) [55]. All the algorithms were executed using MATLAB R2019b.
on a computer with 6 GB of RAM and an Intel i7 CPU running at 2.5 GHz. The comparative analysis aimed to assess the performance of the proposed I-WaOA against current state-of-the-art optimization techniques. The different cases that were studied are presented here:

### 7.1. Testing the I-WaOA Technique on a Set of Commonly Used Test Functions

In this study, we employ the proposed I-WaOA method to evaluate 23 classic functions shown in Tables 2–4 [76–78]. To ensure consistency, the parameters for all situations are set based on Table 5, and the results are recorded after running the experiment 30 times.

<table><caption>**Table 2.** Unimodal functions.</caption><thead><tr><th>F<sub>min</sub></th><th>Range</th><th>Function</th></tr></thead><tbody><tr><td>0</td><td>[-100, 100]</td><td>f<sub>1</sub>(k) = &sum;<sub>j=1</sub><sup>n</sup> k<sub>j</sub><sup>2</sup></td></tr><tr><td>0</td><td>[-10, 10]</td><td>f<sub>2</sub>(k) = &sum;<sub>j=1</sub><sup>n</sup> |k<sub>j</sub>| + &prod;<sub>j=1</sub><sup>n</sup> |k<sub>j</sub>|</td></tr><tr><td>0</td><td>[-100, 100]</td><td>f<sub>3</sub>(k) = &sum;<sub>j=1</sub><sup>n</sup> (&sum;<sub>i=1</sub><sup>j</sup> k<sub>i</sub>)<sup>2</sup></td></tr><tr><td>0</td><td>[-100, 100]</td><td>f<sub>4</sub>(k) = max<sub>j</sub> |k<sub>j</sub>|, 1 &le; j &le; n</td></tr><tr><td>0</td><td>[-30, 30]</td><td>f<sub>5</sub>(k) = &sum;<sub>j=1</sub><sup>n-1</sup> [100(k<sub>j+1</sub> - k<sub>j</sub><sup>2</sup>) + (k<sub>j</sub> - 1)<sup>2</sup>]</td></tr><tr><td>0</td><td>[-100, 100]</td><td>f<sub>6</sub>(k) = &sum;<sub>j=1</sub><sup>n-1</sup> ([k<sub>j</sub> + 0.5]<sup>2</sup>)</td></tr><tr><td>0</td><td>[-1.28, 1.28]</td><td>f<sub>7</sub>(k) = &sum;<sub>j=1</sub><sup>n</sup> jk<sub>j</sub><sup>4</sup> + random(0, 1)</td></tr></tbody></table>

<table><caption>**Table 3.** Multimodal functions.</caption><thead><tr><th>F<sub>min</sub></th><th>Range</th><th>Function</th></tr></thead><tbody><tr><td>-12.56</td><td>[-500, 500]</td><td>f<sub>8</sub>(k) = &sum;<sub>j=1</sub><sup>n</sup> -k<sub>j</sub> sin(&radic;|k<sub>j</sub>|)</td></tr><tr><td>0</td><td>[-5.12, 5.12]</td><td>f<sub>9</sub>(k) = &sum;<sub>j=1</sub><sup>n</sup> [k<sub>j</sub><sup>2</sup> - 10 cos(2&pi;k<sub>j</sub>) + 10]</td></tr><tr><td>0</td><td>[-32, 32]</td><td>f<sub>10</sub>(k) = -20exp(-0.2&radic;<span style="text-decoration: overline;">1/n &sum;<sub>j=1</sub><sup>n</sup> k<sub>j</sub><sup>2</sup></span>) - exp(<span style="text-decoration: overline;">1/n &sum;<sub>j=1</sub><sup>n</sup> cos(2&pi;k<sub>j</sub>) + 20 + e</span>)</td></tr><tr><td>0</td><td>[-600, 600]</td><td>f<sub>11</sub>(k) = <sup>1</sup>&frasl;<sub>4000</sub> &sum;<sub>j=1</sub><sup>n</sup> k<sub>j</sub><sup>2</sup> - &prod;<sub>j=1</sub><sup>n</sup> cos(<sup>k<sub>j</sub></sup>&frasl;<sub>&radic;<span style="text-decoration: overline;">j</span></sub>) + 1</td></tr><tr><td>0</td><td>[-50, 50]</td><td>f<sub>12</sub>(k) = <sup>&pi;</sup>&frasl;<sub>n</sub> {10sin(&pi;z<sub>1</sub>) + &sum;<sub>j=1</sub><sup>n-1</sup> (z<sub>j</sub> - 1)<sup>2</sup> [1 + 10sin<sup>2</sup>(&pi;z<sub>j+1</sub>)] + (z<sub>n</sub> - 1)<sup>2</sup>}</td></tr><tr><td>0</td><td>[-50, 50]</td><td>+ &sum;<sub>j=1</sub><sup>n</sup> u(k<sub>j</sub>, 10, 100, 4)</td></tr><tr><td></td><td></td><td>f<sub>13</sub>(k) = 0.1 {sin<sup>2</sup>(3&pi;k<sub>1</sub>) + &sum;<sub>j=1</sub><sup>n</sup> (k<sub>j</sub> - 1)<sup>2</sup> [1 + sin<sup>2</sup>(3&pi;k<sub>j</sub> + 1)] + (k<sub>n</sub> - 1)<sup>2</sup> [1 + sin<sup>2</sup>(2&pi;k<sub>n</sub>)]}</td></tr><tr><td></td><td></td><td>+ &sum;<sub>j=1</sub><sup>n</sup> u(k<sub>j</sub>, 5, 100, 4)</td></tr><tr><td></td><td></td><td>u(k<sub>j</sub>, v, s, h) = <span style="font-size: 2em; vertical-align: middle;">{</span><table style="display: inline-table; vertical-align: middle; border-collapse: collapse; margin: 0 0.2em;"><tbody><tr><td style="border: none; padding: 0 0.2em;">s(k<sub>j</sub> - v)<sup>h</sup></td><td style="border: none; padding: 0 0.2em;">k<sub>j</sub> &gt; v</td></tr><tr><td style="border: none; padding: 0 0.2em;">0</td><td style="border: none; padding: 0 0.2em;">-v &lt; k<sub>j</sub> &lt; v</td></tr><tr><td style="border: none; padding: 0 0.2em;">s(-k<sub>j</sub> - v)<sup>h</sup></td><td style="border: none; padding: 0 0.2em;">k<sub>j</sub> &lt; -v</td></tr></tbody></table><span style="font-size: 2em; vertical-align: middle;">}</span></td></tr></tbody></table>

<table><caption>**Table 4.** Fixed-dimension multimodal benchmark functions.</caption><thead><tr><th>F<sub>min</sub></th><th>Range</th><th>Function</th></tr></thead><tbody><tr><td>1</td><td>[-65.536, 65.536]</td><td>f<sub>14</sub>(k) = (<sup>1</sup>&frasl;<sub>500</sub> + &sum;<sub>i=1</sub><sup>25</sup> <sup>1</sup>&frasl;<sub>i + &sum;<sub>j=1</sub><sup>2</sup> (k<sub>j</sub> - a<sub>ij</sub>)<sup>6</sup></sub>)<sup>-1</sup></td></tr><tr><td>0.00030</td><td>[-5, 5]</td><td>f<sub>15</sub>(k) = &sum;<sub>j=1</sub><sup>11</sup> [a<sub>j</sub> - <sup>k<sub>j</sub>(b<sub>j</sub><sup>2</sup> + b<sub>j</sub>k<sub>2</sub>)</sup>&frasl;<sub>b<sub>j</sub><sup>2</sup> + b<sub>j</sub>k<sub>3</sub> + k<sub>4</sub></sub>]<sup>2</sup></td></tr></tbody></table>
on a computer with 6 GB of RAM and an Intel i7 CPU running at 2.5 GHz. The comparative analysis aimed to assess the performance of the proposed I-WaOA against current state-of-the-art optimization techniques. The different cases that were studied are presented here:

### 7.1. Testing the I-WaOA Technique on a Set of Commonly Used Test Functions

In this study, we employ the proposed I-WaOA method to evaluate 23 classic functions shown in Tables 2–4 [76–78]. To ensure consistency, the parameters for all situations are set based on Table 5, and the results are recorded after running the experiment 30 times.

<table><caption>**Table 2.** Unimodal functions.</caption><thead><tr><th>F<sub>min</sub></th><th>Range</th><th>Function</th></tr></thead><tbody><tr><td>0</td><td>[-100, 100]</td><td>f<sub>1</sub>(k) = &sum;<sup>n</sup><sub>j=1</sub> k<sub>j</sub><sup>2</sup></td></tr><tr><td>0</td><td>[-10, 10]</td><td>f<sub>2</sub>(k) = &sum;<sup>n</sup><sub>j=1</sub> |k<sub>j</sub>| + &prod;<sup>n</sup><sub>j=1</sub> |k<sub>j</sub>|</td></tr><tr><td>0</td><td>[-100, 100]</td><td>f<sub>3</sub>(k) = &sum;<sup>n</sup><sub>j=1</sub> (&sum;<sup>j</sup><sub>i=1</sub> k<sub>i</sub>)<sup>2</sup></td></tr><tr><td>0</td><td>[-100, 100]</td><td>f<sub>4</sub>(k) = max<sub>j</sub> |k<sub>j</sub>|, 1 &le; j &le; n</td></tr><tr><td>0</td><td>[-30, 30]</td><td>f<sub>5</sub>(k) = &sum;<sup>n-1</sup><sub>j=1</sub> [100(k<sub>j+1</sub> - k<sub>j</sub><sup>2</sup>) + (k<sub>j</sub> - 1)<sup>2</sup>]</td></tr><tr><td>0</td><td>[-100, 100]</td><td>f<sub>6</sub>(k) = &sum;<sup>n-1</sup><sub>j=1</sub> ([k<sub>j</sub> + 0.5])<sup>2</sup></td></tr><tr><td>0</td><td>[-1.28, 1.28]</td><td>f<sub>7</sub>(k) = &sum;<sup>n</sup><sub>j=1</sub> jk<sub>j</sub><sup>4</sup> + random(0, 1)</td></tr></tbody></table>

<table><caption>**Table 3.** Multimodal functions.</caption><thead><tr><th>F<sub>min</sub></th><th>Range</th><th>Function</th></tr></thead><tbody><tr><td>-12.56</td><td>[-500, 500]</td><td>f<sub>8</sub>(k) = &sum;<sup>n</sup><sub>j=1</sub> -k<sub>j</sub> sin(&radic;|k<sub>j</sub>|)</td></tr><tr><td>0</td><td>[-5.12, 5.12]</td><td>f<sub>9</sub>(k) = &sum;<sup>n</sup><sub>j=1</sub> [k<sub>j</sub><sup>2</sup> - 10 cos(2&pi;k<sub>j</sub>) + 10]</td></tr><tr><td>0</td><td>[-32, 32]</td><td>f<sub>10</sub>(k) = -20exp(-0.2&radic;<span style="text-decoration: overline">1/n &sum;<sup>n</sup><sub>j=1</sub> k<sub>j</sub><sup>2</sup></span>) - exp(<span style="text-decoration: overline">1/n &sum;<sup>n</sup><sub>j=1</sub> cos(2&pi;k<sub>j</sub>) + 20 + e</span>)</td></tr><tr><td>0</td><td>[-600, 600]</td><td>f<sub>11</sub>(k) = <sup>1</sup>&frasl;<sub>4000</sub> &sum;<sup>n</sup><sub>j=1</sub> k<sub>j</sub><sup>2</sup> - &prod;<sup>n</sup><sub>j=1</sub> cos(<sup>k<sub>j</sub></sup>&frasl;<sub>&radic;<span style="text-decoration: overline">j</span></sub>) + 1</td></tr><tr><td>0</td><td>[-50, 50]</td><td>f<sub>12</sub>(k) = <sup>&pi;</sup>&frasl;<sub>n</sub> {10sin(&pi;z<sub>1</sub>) + &sum;<sup>n-1</sup><sub>j=1</sub> (z<sub>j</sub> - 1)<sup>2</sup> [1 + 10sin<sup>2</sup>(&pi;z<sub>j+1</sub>)] + (z<sub>n</sub> - 1)<sup>2</sup>}</td></tr><tr><td>0</td><td>[-50, 50]</td><td>+ &sum;<sup>n</sup><sub>j=1</sub> u(k<sub>j</sub>, 10, 100, 4)</td></tr><tr><td></td><td></td><td>f<sub>13</sub>(k) = 0.1 {sin<sup>2</sup>(3&pi;k<sub>1</sub>) + &sum;<sup>n</sup><sub>j=1</sub> (k<sub>j</sub> - 1)<sup>2</sup> [1 + sin<sup>2</sup>(3&pi;k<sub>j</sub> + 1)] + (k<sub>n</sub> - 1)<sup>2</sup> [1 + sin<sup>2</sup>(2&pi;k<sub>n</sub>)]}</td></tr><tr><td></td><td></td><td>+ &sum;<sup>n</sup><sub>j=1</sub> u(k<sub>j</sub>, 5, 100, 4)</td></tr><tr><td></td><td></td><td>u(k<sub>j</sub>, v, s, h) = <span style="font-size: 1.5em; vertical-align: middle;">{</span><table style="display: inline-table; vertical-align: middle; border-collapse: collapse;"><tbody><tr><td style="padding: 0 0.2em; text-align: center;">s(k<sub>j</sub> - v)<sup>h</sup></td><td style="padding: 0 0.2em; text-align: center;">k<sub>j</sub> &gt; v</td></tr><tr><td style="padding: 0 0.2em; text-align: center;">0</td><td style="padding: 0 0.2em; text-align: center;">-v &lt; k<sub>j</sub> &lt; v</td></tr><tr><td style="padding: 0 0.2em; text-align: center;">s(-k<sub>j</sub> - v)<sup>h</sup></td><td style="padding: 0 0.2em; text-align: center;">k<sub>j</sub> &lt; -v</td></tr></tbody></table><span style="font-size: 1.5em; vertical-align: middle;">}</span></td></tr></tbody></table>

<table><caption>**Table 4.** Fixed-dimension multimodal benchmark functions.</caption><thead><tr><th>F<sub>min</sub></th><th>Range</th><th>Function</th></tr></thead><tbody><tr><td>1</td><td>[-65.536, 65.536]</td><td>f<sub>14</sub>(k) = (<sup>1</sup>&frasl;<sub>500</sub> + &sum;<sup>25</sup><sub>i=1</sub> <sup>1</sup>&frasl;<sub>i + &sum;<sup>2</sup><sub>j=1</sub> (k<sub>j</sub> - a<sub>ij</sub>)<sup>6</sup></sub>)<sup>-1</sup></td></tr><tr><td>0.00030</td><td>[-5, 5]</td><td>f<sub>15</sub>(k) = &sum;<sup>11</sup><sub>j=1</sub> [a<sub>j</sub> - <sup>k<sub>j</sub>(b<sub>j</sub><sup>2</sup> + b<sub>j</sub>k<sub>2</sub>)</sup>&frasl;<sub>b<sub>j</sub><sup>2</sup> + b<sub>j</sub>k<sub>3</sub> + k<sub>4</sub></sub>]<sup>2</sup></td></tr></tbody></table>
<table><caption>**Table 4.** *Cont.*</caption><thead><tr><th>F<sub>min</sub></th><th>Range</th><th>Function</th></tr></thead><tbody><tr><td>-1.0316</td><td>[-5, 5]</td><td>f<sub>16</sub>(k) = 4k<sub>1</sub><sup>2</sup> - 2.1k<sub>1</sub><sup>4</sup> + <sup>1</sup>&frasl;<sub>3</sub>k<sub>1</sub><sup>6</sup> + k<sub>1</sub>k<sub>2</sub> - 4k<sub>2</sub><sup>2</sup> + 4k<sub>2</sub><sup>4</sup></td></tr><tr><td>0.398</td><td>[-5, 5]</td><td>f<sub>17</sub>(k) = (k<sub>2</sub> - <sup>5.1</sup>&frasl;<sub>4π<sup>2</sup></sub>k<sub>1</sub><sup>2</sup> + <sup>5</sup>&frasl;<sub>π</sub>k<sub>1</sub> - 6)<sup>2</sup> + 10(1 - <sup>1</sup>&frasl;<sub>8π</sub>)cos k<sub>1</sub> + 10</td></tr><tr><td>3</td><td>[-2, 2]</td><td>f<sub>18</sub>(k) = [1 + (k<sub>1</sub> + k<sub>2</sub> + 1)<sup>2</sup>(19 - 14k<sub>1</sub> + 3k<sub>1</sub><sup>2</sup> - 14k<sub>2</sub> + 6k<sub>1</sub>k<sub>2</sub> + 3k<sub>2</sub><sup>2</sup>)] * [30 + (2k<sub>1</sub> - 3k<sub>2</sub>)<sup>2</sup>(18 - 32k<sub>1</sub> + 12k<sub>1</sub><sup>2</sup> + 48k<sub>2</sub> - 36k<sub>1</sub>k<sub>2</sub> + 27k<sub>2</sub><sup>2</sup>)]</td></tr><tr><td>-3.86</td><td>[1, 3]</td><td>f<sub>19</sub>(k) = -Σ<sub>j=1</sub><sup>4</sup> c<sub>j</sub> exp { -Σ<sub>i=1</sub><sup>3</sup> a<sub>ji</sub>(k<sub>i</sub> - p<sub>ji</sub>)<sup>2</sup>}</td></tr><tr><td>-3.32</td><td>[0, 1]</td><td>f<sub>20</sub>(k) = -Σ<sub>j=1</sub><sup>4</sup> c<sub>j</sub> exp { -Σ<sub>i=1</sub><sup>6</sup> a<sub>ji</sub>(k<sub>i</sub> - p<sub>ji</sub>)<sup>2</sup>}</td></tr><tr><td>-10.153</td><td>[0, 10]</td><td>f<sub>21</sub>(k) = -Σ<sub>j=1</sub><sup>5</sup> [(k - a<sub>j</sub>)(k - a<sub>j</sub>)<sup>T</sup> + c<sub>j</sub>]<sup>-1</sup></td></tr><tr><td>-10.402</td><td>[0, 10]</td><td>f<sub>22</sub>(k) = -Σ<sub>j=1</sub><sup>7</sup> [(k - a<sub>j</sub>)(k - a<sub>j</sub>)<sup>T</sup> + c<sub>j</sub>]<sup>-1</sup></td></tr><tr><td>-10.536</td><td>[0, 10]</td><td>f<sub>23</sub>(k) = -Σ<sub>j=1</sub><sup>10</sup> [(k - a<sub>j</sub>)(k - a<sub>j</sub>)<sup>T</sup> + c<sub>j</sub>]<sup>-1</sup></td></tr></tbody></table>

<table><caption>**Table 5.** The parameters selected for the optimization algorithms.</caption><thead><tr><th>Algorithm</th><th>Parameter</th><th>Value</th></tr></thead><tbody><tr><td rowspan="4">SCSO [74]</td><td>Phases control range (R)</td><td>[-2rg, 2rg]</td></tr><tr><td>Sensitivity range (rg)</td><td>[2, 0]</td></tr><tr><td>Populations</td><td>30</td></tr><tr><td>Maximum iteration</td><td>300</td></tr><tr><td rowspan="3">AHA [70]</td><td>Migration coefficient</td><td>2n</td></tr><tr><td>Populations</td><td>30</td></tr><tr><td>Maximum iteration</td><td>300</td></tr><tr><td rowspan="4">DO [71]</td><td>α</td><td>[0, 1]</td></tr><tr><td>k</td><td>[0, 1]</td></tr><tr><td>Populations</td><td>30</td></tr><tr><td>Maximum iteration</td><td>300</td></tr><tr><td rowspan="3">CDO [69]</td><td>r</td><td>[0, 1]</td></tr><tr><td>Populations</td><td>30</td></tr><tr><td>Maximum iteration</td><td>300</td></tr><tr><td rowspan="4">HS [72]</td><td>HMCR</td><td>0.95</td></tr><tr><td>PAR</td><td>0.45</td></tr><tr><td>Populations</td><td>30</td></tr><tr><td>Maximum iteration</td><td>300</td></tr><tr><td rowspan="2">ZOA [75]</td><td>Populations</td><td>30</td></tr><tr><td>Maximum iteration</td><td>300</td></tr><tr><td rowspan="2">ARO [73]</td><td>Populations</td><td>30</td></tr><tr><td>Maximum iteration</td><td>300</td></tr><tr><td rowspan="2">WaOA [55]</td><td>Populations</td><td>30</td></tr><tr><td>Maximum iteration</td><td>300</td></tr><tr><td rowspan="2">I-WaOA</td><td>Populations</td><td>30</td></tr><tr><td>Maximum iteration</td><td>300</td></tr></tbody></table>

### 7.1.1. Statistical Results Analysis

This section presents a comparative analysis of the proposed I-WaOA with other optimization methods, including SCSO, AHA, DO, HS, CDO, ZOA, ARO, and the conventional WaOA. Statistical results, as shown in Table 6, encompass the mean, worst, SD, and best among all algorithm approaches.
<table><caption>**Table 4.** *Cont.*</caption><thead><tr><th>F<sub>min</sub></th><th>Range</th><th>Function</th></tr></thead><tbody><tr><td>-1.0316</td><td>[-5, 5]</td><td>f<sub>16</sub>(k) = 4k<sub>1</sub><sup>2</sup> - 2.1k<sub>1</sub><sup>4</sup> + <sup>1</sup>&frasl;<sub>3</sub>k<sub>1</sub><sup>6</sup> + k<sub>1</sub>k<sub>2</sub> - 4k<sub>2</sub><sup>2</sup> + 4k<sub>2</sub><sup>4</sup></td></tr><tr><td>0.398</td><td>[-5, 5]</td><td>f<sub>17</sub>(k) = (k<sub>2</sub> - <sup>5.1</sup>&frasl;<sub>4π<sup>2</sup></sub>k<sub>1</sub><sup>2</sup> + <sup>5</sup>&frasl;<sub>π</sub>k<sub>1</sub> - 6)<sup>2</sup> + 10(1 - <sup>1</sup>&frasl;<sub>8π</sub>)cos k<sub>1</sub> + 10</td></tr><tr><td>3</td><td>[-2, 2]</td><td>f<sub>18</sub>(k) = [1 + (k<sub>1</sub> + k<sub>2</sub> + 1)<sup>2</sup>(19 - 14k<sub>1</sub> + 3k<sub>1</sub><sup>2</sup> - 14k<sub>2</sub> + 6k<sub>1</sub>k<sub>2</sub> + 3k<sub>2</sub><sup>2</sup>)] * [30 + (2k<sub>1</sub> - 3k<sub>2</sub>)<sup>2</sup>(18 - 32k<sub>1</sub> + 12k<sub>1</sub><sup>2</sup> + 48k<sub>2</sub> - 36k<sub>1</sub>k<sub>2</sub> + 27k<sub>2</sub><sup>2</sup>)]</td></tr><tr><td>-3.86</td><td>[1, 3]</td><td>f<sub>19</sub>(k) = -Σ<sub>j=1</sub><sup>4</sup> c<sub>j</sub> exp { -Σ<sub>i=1</sub><sup>3</sup> a<sub>ji</sub>(k<sub>i</sub> - p<sub>ji</sub>)<sup>2</sup>}</td></tr><tr><td>-3.32</td><td>[0, 1]</td><td>f<sub>20</sub>(k) = -Σ<sub>j=1</sub><sup>4</sup> c<sub>j</sub> exp { -Σ<sub>i=1</sub><sup>6</sup> a<sub>ji</sub>(k<sub>i</sub> - p<sub>ji</sub>)<sup>2</sup>}</td></tr><tr><td>-10.153</td><td>[0, 10]</td><td>f<sub>21</sub>(k) = -Σ<sub>j=1</sub><sup>5</sup> [(k - a<sub>j</sub>)(k - a<sub>j</sub>)<sup>T</sup> + c<sub>j</sub>]<sup>-1</sup></td></tr><tr><td>-10.402</td><td>[0, 10]</td><td>f<sub>22</sub>(k) = -Σ<sub>j=1</sub><sup>7</sup> [(k - a<sub>j</sub>)(k - a<sub>j</sub>)<sup>T</sup> + c<sub>j</sub>]<sup>-1</sup></td></tr><tr><td>-10.536</td><td>[0, 10]</td><td>f<sub>23</sub>(k) = -Σ<sub>j=1</sub><sup>10</sup> [(k - a<sub>j</sub>)(k - a<sub>j</sub>)<sup>T</sup> + c<sub>j</sub>]<sup>-1</sup></td></tr></tbody></table>

<table><caption>**Table 5.** The parameters selected for the optimization algorithms.</caption><thead><tr><th>Algorithm</th><th>Parameter</th><th>Value</th></tr></thead><tbody><tr><td rowspan="3">SCSO [74]</td><td>Phases control range (R)</td><td>[-2rg, 2rg]</td></tr><tr><td>Sensitivity range (rg)</td><td>[2, 0]</td></tr><tr><td>Populations</td><td>30</td></tr><tr><td rowspan="2">AHA [70]</td><td>Maximum iteration</td><td>300</td></tr><tr><td>Migration coefficient</td><td>2n</td></tr><tr><td rowspan="3">DO [71]</td><td>Populations</td><td>30</td></tr><tr><td>Maximum iteration</td><td>300</td></tr><tr><td>α</td><td>[0, 1]</td></tr><tr><td rowspan="3">CDO [69]</td><td>k</td><td>[0, 1]</td></tr><tr><td>Populations</td><td>30</td></tr><tr><td>Maximum iteration</td><td>300</td></tr><tr><td rowspan="3">HS [72]</td><td>r</td><td>[0, 1]</td></tr><tr><td>Populations</td><td>30</td></tr><tr><td>Maximum iteration</td><td>300</td></tr><tr><td rowspan="2">ZOA [75]</td><td>HMCR</td><td>0.95</td></tr><tr><td>PAR</td><td>0.45</td></tr><tr><td rowspan="2">ARO [73]</td><td>Populations</td><td>30</td></tr><tr><td>Maximum iteration</td><td>300</td></tr><tr><td rowspan="2">WaOA [55]</td><td>Populations</td><td>30</td></tr><tr><td>Maximum iteration</td><td>300</td></tr><tr><td rowspan="2">I-WaOA</td><td>Populations</td><td>30</td></tr><tr><td>Maximum iteration</td><td>300</td></tr></tbody></table>

### 7.1.1. Statistical Results Analysis

This section presents a comparative analysis of the proposed I-WaOA with other optimization methods, including SCSO, AHA, DO, HS, CDO, ZOA, ARO, and the conventional WaOA. Statistical results, as shown in Table 6, encompass the mean, worst, SD, and best among all algorithm approaches.
<table><caption>Table 6. Statistical findings of various optimization algorithms on the standard functions.</caption><thead><tr><th>Fun.</th><th>I-WaOA</th><th>WaOA</th><th>SCSO</th><th>AHA</th><th>DO</th><th>HS</th><th>CDO</th><th>ZOA</th><th>ARO</th></tr></thead><tbody><tr><td>F1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>0</td><td>1.69 × 10<sup>-172</sup></td><td>6.03 × 10<sup>-74</sup></td><td>4.44 × 10<sup>-91</sup></td><td>0.000265</td><td>2519.226</td><td>5.30 × 10<sup>-85</sup></td><td>1.41 × 10<sup>-152</sup></td><td>7.05 × 10<sup>-37</sup></td></tr><tr><td>Average</td><td>0</td><td>2.66 × 10<sup>-168</sup></td><td>3.76 × 10<sup>-65</sup></td><td>2.46 × 10<sup>-82</sup></td><td>0.000975</td><td>3492.24</td><td>4.54 × 10<sup>-80</sup></td><td>1.53 × 10<sup>-148</sup></td><td>8.87 × 10<sup>-31</sup></td></tr><tr><td>Worst</td><td>0</td><td>2.46 × 10<sup>-167</sup></td><td>2.81 × 10<sup>-64</sup></td><td>1.83 × 10<sup>-81</sup></td><td>0.001716</td><td>4462.583</td><td>2.87 × 10<sup>-79</sup></td><td>5.97 × 10<sup>-148</sup></td><td>8.86 × 10<sup>-30</sup></td></tr><tr><td>SD</td><td>0</td><td>0</td><td>9.05 × 10<sup>-65</sup></td><td>5.88 × 10<sup>-82</sup></td><td>0.000597</td><td>566.3726</td><td>9.20 × 10<sup>-80</sup></td><td>2.39 × 10<sup>-148</sup></td><td>2.80 × 10<sup>-30</sup></td></tr><tr><td>F2</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>2.64 × 10<sup>-196</sup></td><td>4.70 × 10<sup>-90</sup></td><td>2.92 × 10<sup>-40</sup></td><td>8.08 × 10<sup>-50</sup></td><td>0.007229</td><td>10.28939</td><td>7.61 × 10<sup>-43</sup></td><td>3.37 × 10<sup>-81</sup></td><td>3.69 × 10<sup>-22</sup></td></tr><tr><td>Average</td><td>1.36 × 10<sup>-189</sup></td><td>8.26 × 10<sup>-86</sup></td><td>3.26 × 10<sup>-36</sup></td><td>1.18 × 10<sup>-42</sup></td><td>0.013546</td><td>13.19099</td><td>4.73 × 10<sup>-41</sup></td><td>6.18 × 10<sup>-78</sup></td><td>1.08 × 10<sup>-18</sup></td></tr><tr><td>Worst</td><td>2.70 × 10<sup>-188</sup></td><td>6.79 × 10<sup>-85</sup></td><td>2.40 × 10<sup>-35</sup></td><td>1.32 × 10<sup>-41</sup></td><td>0.020851</td><td>16.2306</td><td>3.79 × 10<sup>-40</sup></td><td>7.16 × 10<sup>-77</sup></td><td>9.32 × 10<sup>-18</sup></td></tr><tr><td>SD</td><td>0</td><td>1.55 × 10<sup>-85</sup></td><td>5.95 × 10<sup>-36</sup></td><td>3.63 × 10<sup>-42</sup></td><td>0.003821</td><td>1.883619</td><td>9.14 × 10<sup>-41</sup></td><td>1.70 × 10<sup>-77</sup></td><td>2.40 × 10<sup>-18</sup></td></tr><tr><td>F3</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>3.16 × 10<sup>-317</sup></td><td>1.01 × 10<sup>-133</sup></td><td>2.68 × 10<sup>-65</sup></td><td>1.72 × 10<sup>-92</sup></td><td>19.86247</td><td>21433.18</td><td>1.75 × 10<sup>-71</sup></td><td>4.66 × 10<sup>-104</sup></td><td>1.62 × 10<sup>-30</sup></td></tr><tr><td>Average</td><td>2.20 × 10<sup>-304</sup></td><td>6.13 × 10<sup>-123</sup></td><td>4.46 × 10<sup>-59</sup></td><td>2.10 × 10<sup>-73</sup></td><td>126.9791</td><td>35753.52</td><td>5.28 × 10<sup>-59</sup></td><td>1.97 × 10<sup>-88</sup></td><td>1.98 × 10<sup>-24</sup></td></tr><tr><td>Worst</td><td>2.00 × 10<sup>-303</sup></td><td>5.76 × 10<sup>-122</sup></td><td>4.39 × 10<sup>-58</sup></td><td>2.10 × 10<sup>-72</sup></td><td>283.6619</td><td>49787.79</td><td>5.28 × 10<sup>-58</sup></td><td>1.97 × 10<sup>-87</sup></td><td>1.83 × 10<sup>-23</sup></td></tr><tr><td>SD</td><td>0</td><td>1.81 × 10<sup>-122</sup></td><td>1.39 × 10<sup>-58</sup></td><td>6.63 × 10<sup>-73</sup></td><td>85.14181</td><td>8923.102</td><td>1.67 × 10<sup>-58</sup></td><td>6.22 × 10<sup>-88</sup></td><td>5.73 × 10<sup>-24</sup></td></tr><tr><td>F4</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>3.27 × 10<sup>-176</sup></td><td>1.51 × 10<sup>-82</sup></td><td>5.90 × 10<sup>-35</sup></td><td>1.07 × 10<sup>-43</sup></td><td>0.820293</td><td>34.00865</td><td>6.41 × 10<sup>-40</sup></td><td>6.48 × 10<sup>-71</sup></td><td>2.28 × 10<sup>-16</sup></td></tr><tr><td>Average</td><td>1.05 × 10<sup>-168</sup></td><td>1.40 × 10<sup>-79</sup></td><td>1.20 × 10<sup>-29</sup></td><td>4.83 × 10<sup>-39</sup></td><td>3.545085</td><td>41.78794</td><td>5.54 × 10<sup>-38</sup></td><td>1.14 × 10<sup>-68</sup></td><td>9.42 × 10<sup>-14</sup></td></tr><tr><td>Worst</td><td>9.88 × 10<sup>-168</sup></td><td>1.34 × 10<sup>-78</sup></td><td>1.19 × 10<sup>-28</sup></td><td>3.04 × 10<sup>-38</sup></td><td>6.640496</td><td>45.66198</td><td>2.52 × 10<sup>-37</sup></td><td>8.04 × 10<sup>-68</sup></td><td>7.15 × 10<sup>-13</sup></td></tr><tr><td>SD</td><td>0</td><td>4.22 × 10<sup>-79</sup></td><td>3.75 × 10<sup>-29</sup></td><td>9.54 × 10<sup>-39</sup></td><td>1.935587</td><td>3.361813</td><td>7.92 × 10<sup>-38</sup></td><td>2.58 × 10<sup>-68</sup></td><td>2.22 × 10<sup>-13</sup></td></tr><tr><td>F5</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>0</td><td>0</td><td>2.72 × 10<sup>1</sup></td><td>2.66 × 10<sup>1</sup></td><td>25.51494</td><td>1565934</td><td>2.78 × 10<sup>1</sup></td><td>2.83 × 10<sup>1</sup></td><td>1.76 × 10<sup>-1</sup></td></tr><tr><td>Average</td><td>0</td><td>0</td><td>2.84 × 10<sup>1</sup></td><td>2.74 × 10<sup>1</sup></td><td>29.92931</td><td>2104121</td><td>2.82 × 10<sup>1</sup></td><td>2.86 × 10<sup>1</sup></td><td>1.50</td></tr><tr><td>Worst</td><td>0</td><td>0</td><td>2.88 × 10<sup>1</sup></td><td>2.87 × 10<sup>1</sup></td><td>45.20891</td><td>3043646</td><td>2.87 × 10<sup>1</sup></td><td>2.89 × 10<sup>1</sup></td><td>6.61</td></tr><tr><td>SD</td><td>0</td><td>0</td><td>5.40 × 10<sup>-1</sup></td><td>6.14 × 10<sup>-1</sup></td><td>5.543385</td><td>420298.4</td><td>2.84 × 10<sup>-1</sup></td><td>1.82 × 10<sup>-1</sup></td><td>2.42</td></tr><tr><td>F6</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>0</td><td>0</td><td>1.49</td><td>3.43 × 10<sup>-2</sup></td><td>0.000109</td><td>2850.629</td><td>7.50</td><td>2.44</td><td>2.07 × 10<sup>-3</sup></td></tr><tr><td>Average</td><td>0</td><td>0</td><td>2.19</td><td>3.26 × 10<sup>-1</sup></td><td>0.000235</td><td>3806.061</td><td>7.50</td><td>2.99</td><td>1.85 × 10<sup>-2</sup></td></tr><tr><td>Worst</td><td>0</td><td>0</td><td>3.07</td><td>6.58 × 10<sup>-1</sup></td><td>0.000602</td><td>4724.147</td><td>7.50</td><td>3.72</td><td>4.01 × 10<sup>-2</sup></td></tr><tr><td>SD</td><td>0</td><td>0</td><td>6.00</td><td>2.37 × 10<sup>-1</sup></td><td>0.000138</td><td>607.6398</td><td>0</td><td>4.79 × 10<sup>-1</sup></td><td>1.27 × 10<sup>-2</sup></td></tr><tr><td>F7</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>9.92 × 10<sup>-7</sup></td><td>2.86 × 10<sup>-6</sup></td><td>8.83 × 10<sup>-7</sup></td><td>9.48 × 10<sup>-5</sup></td><td>0.018263</td><td>0.925058</td><td>9.82 × 10<sup>-6</sup></td><td>2.62 × 10<sup>-5</sup></td><td>2.06 × 10<sup>-4</sup></td></tr><tr><td>Average</td><td>5.88 × 10<sup>-5</sup></td><td>6.65 × 10<sup>-5</sup></td><td>5.40 × 10<sup>-4</sup></td><td>4.65 × 10<sup>-4</sup></td><td>0.04415</td><td>1.261838</td><td>2.01 × 10<sup>-4</sup></td><td>1.39 × 10<sup>-4</sup></td><td>9.75 × 10<sup>-4</sup></td></tr><tr><td>Worst</td><td>1.53 × 10<sup>-4</sup></td><td>1.86 × 10<sup>-4</sup></td><td>6.53 × 10<sup>-3</sup></td><td>1.16 × 10<sup>-3</sup></td><td>0.098175</td><td>1.897037</td><td>4.58 × 10<sup>-4</sup></td><td>3.49 × 10<sup>-4</sup></td><td>1.92 × 10<sup>-3</sup></td></tr><tr><td>SD</td><td>4.31 × 10<sup>-5</sup></td><td>5.29 × 10<sup>-5</sup></td><td>1.44 × 10<sup>-3</sup></td><td>3.00 × 10<sup>-4</sup></td><td>0.024073</td><td>0.268522</td><td>1.51 × 10<sup>-4</sup></td><td>9.20 × 10<sup>-5</sup></td><td>4.80 × 10<sup>-4</sup></td></tr><tr><td>F8</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>-1.26 × 10<sup>4</sup></td><td>-9.02 × 10<sup>3</sup></td><td>-7.41 × 10<sup>3</sup></td><td>-1.12 × 10<sup>4</sup></td><td>-8685.03</td><td>-11610.1</td><td>-4.72 × 10<sup>3</sup></td><td>-7.26 × 10<sup>3</sup></td><td>-1.05 × 10<sup>4</sup></td></tr><tr><td>Average</td><td>-1.26 × 10<sup>4</sup></td><td>-8.34 × 10<sup>3</sup></td><td>-6.43 × 10<sup>3</sup></td><td>-1.03 × 10<sup>4</sup></td><td>-7626.41</td><td>-11407.5</td><td>-3.71 × 10<sup>3</sup></td><td>-6.59 × 10<sup>3</sup></td><td>-9.49 × 10<sup>3</sup></td></tr><tr><td>Worst</td><td>-1.26 × 10<sup>4</sup></td><td>-7.52 × 10<sup>3</sup></td><td>-4.99 × 10<sup>3</sup></td><td>-8.99 × 10<sup>3</sup></td><td>-6129.31</td><td>-11106.4</td><td>-2.86 × 10<sup>3</sup></td><td>-5.85 × 10<sup>3</sup></td><td>-8.00 × 10<sup>3</sup></td></tr><tr><td>SD</td><td>7.52 × 10<sup>-10</sup></td><td>5.03 × 10<sup>2</sup></td><td>8.17 × 10<sup>2</sup></td><td>5.93 × 10<sup>2</sup></td><td>765.964</td><td>144.0269</td><td>5.46 × 10<sup>2</sup></td><td>5.75 × 10<sup>2</sup></td><td>6.43 × 10<sup>2</sup></td></tr><tr><td>F9</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>0</td><td>0</td><td>0</td><td>0</td><td>21.06603</td><td>50.57476</td><td>0</td><td>0</td><td>0</td></tr><tr><td>Average</td><td>0</td><td>0</td><td>0</td><td>0</td><td>40.21929</td><td>60.92757</td><td>9.60 × 10<sup>1</sup></td><td>0</td><td>0</td></tr><tr><td>Worst</td><td>0</td><td>0</td><td>0</td><td>0</td><td>59.40942</td><td>71.63848</td><td>2.58 × 10<sup>2</sup></td><td>0</td><td>0</td></tr><tr><td>SD</td><td>0</td><td>0</td><td>0</td><td>0</td><td>15.38431</td><td>6.732379</td><td>1.24 × 10<sup>2</sup></td><td>0</td><td>0</td></tr><tr><td>F10</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>8.88 × 10<sup>-16</sup></td><td>8.88 × 10<sup>-16</sup></td><td>8.88 × 10<sup>-16</sup></td><td>8.88 × 10<sup>-16</sup></td><td>0.003418</td><td>9.380315</td><td>4.44 × 10<sup>-15</sup></td><td>8.88 × 10<sup>-16</sup></td><td>8.88 × 10<sup>-16</sup></td></tr><tr><td>Average</td><td>8.88 × 10<sup>-16</sup></td><td>2.66 × 10<sup>-15</sup></td><td>8.88 × 10<sup>-16</sup></td><td>8.88 × 10<sup>-16</sup></td><td>0.006712</td><td>11.32595</td><td>4.44 × 10<sup>-15</sup></td><td>8.88 × 10<sup>-16</sup></td><td>8.88 × 10<sup>-16</sup></td></tr><tr><td>Worst</td><td>8.88 × 10<sup>-16</sup></td><td>4.44 × 10<sup>-15</sup></td><td>8.88 × 10<sup>-16</sup></td><td>8.88 × 10<sup>-16</sup></td><td>0.00951</td><td>12.76612</td><td>4.44 × 10<sup>-15</sup></td><td>8.88 × 10<sup>-16</sup></td><td>8.88 × 10<sup>-16</sup></td></tr><tr><td>SD</td><td>0</td><td>1.87 × 10<sup>-15</sup></td><td>0</td><td>0</td><td>0.00188</td><td>0.95227</td><td>0</td><td>0</td><td>0</td></tr><tr><td>F11</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.000995</td><td>27.86147</td><td>0</td><td>0</td><td>0</td></tr><tr><td>Average</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.016851</td><td>32.78121</td><td>1.61 × 10<sup>-3</sup></td><td>0</td><td>0</td></tr><tr><td>Worst</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.044359</td><td>47.68918</td><td>1.61 × 10<sup>-2</sup></td><td>0</td><td>0</td></tr><tr><td>SD</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.013777</td><td>6.341414</td><td>5.08 × 10<sup>-3</sup></td><td>0</td><td>0</td></tr><tr><td>F12</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>1.57 × 10<sup>-32</sup></td><td>1.57 × 10<sup>-32</sup></td><td>4.72 × 10<sup>-2</sup></td><td>2.13 × 10<sup>-3</sup></td><td>6.11 × 10<sup>-6</sup></td><td>176966.7</td><td>1.11</td><td>1.38 × 10<sup>-1</sup></td><td>2.74 × 10<sup>-4</sup></td></tr><tr><td>Average</td><td>1.57 × 10<sup>-32</sup></td><td>1.57 × 10<sup>-32</sup></td><td>1.27 × 10<sup>-1</sup></td><td>1.22 × 10<sup>-2</sup></td><td>0.324104</td><td>352773.6</td><td>1.44</td><td>2.38 × 10<sup>-1</sup></td><td>2.23 × 10<sup>-3</sup></td></tr><tr><td>Worst</td><td>1.57 × 10<sup>-32</sup></td><td>1.57 × 10<sup>-32</sup></td><td>3.54 × 10<sup>-1</sup></td><td>2.28 × 10<sup>-2</sup></td><td>2.799928</td><td>621526.7</td><td>1.67</td><td>3.64 × 10<sup>-1</sup></td><td>8.00 × 10<sup>-3</sup></td></tr><tr><td>SD</td><td>2.88 × 10<sup>-48</sup></td><td>2.88 × 10<sup>-48</sup></td><td>8.71 × 10<sup>-2</sup></td><td>6.39 × 10<sup>-3</sup></td><td>0.876357</td><td>142951.5</td><td>2.90 × 10<sup>-01</sup></td><td>7.64 × 10<sup>-2</sup></td><td>2.80 × 10<sup>-3</sup></td></tr><tr><td>F13</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>1.35 × 10<sup>-32</sup></td><td>1.35 × 10<sup>-32</sup></td><td>1.73</td><td>6.67 × 10<sup>-1</sup></td><td>4.38 × 10<sup>-5</sup></td><td>1360677</td><td>4.87 × 10<sup>-1</sup></td><td>1.83</td><td>9.98 × 10<sup>-4</sup></td></tr><tr><td>Average</td><td>1.35 × 10<sup>-32</sup></td><td>1.35 × 10<sup>-32</sup></td><td>2.29</td><td>2.20</td><td>0.006739</td><td>3046984</td><td>6.20 × 10<sup>-1</sup></td><td>2.25</td><td>2.60 × 10<sup>-2</sup></td></tr><tr><td>Worst</td><td>1.35 × 10<sup>-32</sup></td><td>1.35 × 10<sup>-32</sup></td><td>2.79</td><td>2.78</td><td>0.044155</td><td>6761426</td><td>7.69 × 10<sup>-1</sup></td><td>2.71</td><td>1.42 × 10<sup>-1</sup></td></tr><tr><td>SD</td><td>2.88 × 10<sup>-48</sup></td><td>2.88 × 10<sup>-48</sup></td><td>3.67 × 10<sup>-1</sup></td><td>5.95 × 10<sup>-1</sup></td><td>0.013908</td><td>1730831</td><td>1.01 × 10<sup>-1</sup></td><td>2.84 × 10<sup>-1</sup></td><td>4.17 × 10<sup>-2</sup></td></tr><tr><td>F14</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>9.98 × 10<sup>-1</sup></td><td>9.98 × 10<sup>-1</sup></td><td>9.98 × 10<sup>-1</sup></td><td>9.98 × 10<sup>-1</sup></td><td>9.98 × 10<sup>-1</sup></td><td>0.998004</td><td>2.02</td><td>9.98 × 10<sup>-1</sup></td><td>9.98 × 10<sup>-1</sup></td></tr><tr><td>Average</td><td>9.98 × 10<sup>-1</sup></td><td>9.98 × 10<sup>-1</sup></td><td>2.97</td><td>1.20</td><td>1.196809</td><td>0.998012</td><td>1.23 × 10<sup>1</sup></td><td>3.27</td><td>1.20</td></tr><tr><td>Worst</td><td>9.98 × 10<sup>-1</sup></td><td>9.98 × 10<sup>-1</sup></td><td>1.08 × 10<sup>-1</sup></td><td>2.98</td><td>1.992031</td><td>0.998071</td><td>1.83 × 10<sup>1</sup></td><td>6.90</td><td>2.98</td></tr><tr><td>SD</td><td>7.40 × 10<sup>-17</sup></td><td>0</td><td>2.91</td><td>6.27 × 10<sup>-1</sup></td><td>0.419119</td><td>2.11 × 10<sup>-5</sup></td><td>4.35</td><td>2.23</td><td>6.27 × 10<sup>-1</sup></td></tr></tbody></table>
<table><caption>Table 6. Statistical findings of various optimization algorithms on the standard functions.</caption><thead><tr><th>Fun.</th><th>I-WaOA</th><th>WaOA</th><th>SCSO</th><th>AHA</th><th>DO</th><th>HS</th><th>CDO</th><th>ZOA</th><th>ARO</th></tr></thead><tbody><tr><td>F1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>0</td><td>1.69 × 10<sup>-172</sup></td><td>6.03 × 10<sup>-74</sup></td><td>4.44 × 10<sup>-91</sup></td><td>0.000265</td><td>2519.226</td><td>5.30 × 10<sup>-85</sup></td><td>1.41 × 10<sup>-152</sup></td><td>7.05 × 10<sup>-37</sup></td></tr><tr><td>Average</td><td>0</td><td>2.66 × 10<sup>-168</sup></td><td>3.76 × 10<sup>-65</sup></td><td>2.46 × 10<sup>-82</sup></td><td>0.000975</td><td>3492.24</td><td>4.54 × 10<sup>-80</sup></td><td>1.53 × 10<sup>-148</sup></td><td>8.87 × 10<sup>-31</sup></td></tr><tr><td>Worst</td><td>0</td><td>2.46 × 10<sup>-167</sup></td><td>2.81 × 10<sup>-64</sup></td><td>1.83 × 10<sup>-81</sup></td><td>0.001716</td><td>4462.583</td><td>2.87 × 10<sup>-79</sup></td><td>5.97 × 10<sup>-148</sup></td><td>8.86 × 10<sup>-30</sup></td></tr><tr><td>SD</td><td>0</td><td>0</td><td>9.05 × 10<sup>-65</sup></td><td>5.88 × 10<sup>-82</sup></td><td>0.000597</td><td>566.3726</td><td>9.20 × 10<sup>-80</sup></td><td>2.39 × 10<sup>-148</sup></td><td>2.80 × 10<sup>-30</sup></td></tr><tr><td>F2</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>2.64 × 10<sup>-196</sup></td><td>4.70 × 10<sup>-90</sup></td><td>2.92 × 10<sup>-40</sup></td><td>8.08 × 10<sup>-50</sup></td><td>0.007229</td><td>10.28939</td><td>7.61 × 10<sup>-43</sup></td><td>3.37 × 10<sup>-81</sup></td><td>3.69 × 10<sup>-22</sup></td></tr><tr><td>Average</td><td>1.36 × 10<sup>-189</sup></td><td>8.26 × 10<sup>-86</sup></td><td>3.26 × 10<sup>-36</sup></td><td>1.18 × 10<sup>-42</sup></td><td>0.013546</td><td>13.19099</td><td>4.73 × 10<sup>-41</sup></td><td>6.18 × 10<sup>-78</sup></td><td>1.08 × 10<sup>-18</sup></td></tr><tr><td>Worst</td><td>2.70 × 10<sup>-188</sup></td><td>6.79 × 10<sup>-85</sup></td><td>2.40 × 10<sup>-35</sup></td><td>1.32 × 10<sup>-41</sup></td><td>0.020851</td><td>16.2306</td><td>3.79 × 10<sup>-40</sup></td><td>7.16 × 10<sup>-77</sup></td><td>9.32 × 10<sup>-18</sup></td></tr><tr><td>SD</td><td>0</td><td>1.55 × 10<sup>-85</sup></td><td>5.95 × 10<sup>-36</sup></td><td>3.63 × 10<sup>-42</sup></td><td>0.003821</td><td>1.883619</td><td>9.14 × 10<sup>-41</sup></td><td>1.70 × 10<sup>-77</sup></td><td>2.40 × 10<sup>-18</sup></td></tr><tr><td>F3</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>3.16 × 10<sup>-317</sup></td><td>1.01 × 10<sup>-133</sup></td><td>2.68 × 10<sup>-65</sup></td><td>1.72 × 10<sup>-92</sup></td><td>19.86247</td><td>21433.18</td><td>1.75 × 10<sup>-71</sup></td><td>4.66 × 10<sup>-104</sup></td><td>1.62 × 10<sup>-30</sup></td></tr><tr><td>Average</td><td>2.20 × 10<sup>-304</sup></td><td>6.13 × 10<sup>-123</sup></td><td>4.46 × 10<sup>-59</sup></td><td>2.10 × 10<sup>-73</sup></td><td>126.9791</td><td>35753.52</td><td>5.28 × 10<sup>-59</sup></td><td>1.97 × 10<sup>-88</sup></td><td>1.98 × 10<sup>-24</sup></td></tr><tr><td>Worst</td><td>2.00 × 10<sup>-303</sup></td><td>5.76 × 10<sup>-122</sup></td><td>4.39 × 10<sup>-58</sup></td><td>2.10 × 10<sup>-72</sup></td><td>283.6619</td><td>49787.79</td><td>5.28 × 10<sup>-58</sup></td><td>1.97 × 10<sup>-87</sup></td><td>1.83 × 10<sup>-23</sup></td></tr><tr><td>SD</td><td>0</td><td>1.81 × 10<sup>-122</sup></td><td>1.39 × 10<sup>-58</sup></td><td>6.63 × 10<sup>-73</sup></td><td>85.14181</td><td>8923.102</td><td>1.67 × 10<sup>-58</sup></td><td>6.22 × 10<sup>-88</sup></td><td>5.73 × 10<sup>-24</sup></td></tr><tr><td>F4</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>3.27 × 10<sup>-176</sup></td><td>1.51 × 10<sup>-82</sup></td><td>5.90 × 10<sup>-35</sup></td><td>1.07 × 10<sup>-43</sup></td><td>0.820293</td><td>34.00865</td><td>6.41 × 10<sup>-40</sup></td><td>6.48 × 10<sup>-71</sup></td><td>2.28 × 10<sup>-16</sup></td></tr><tr><td>Average</td><td>1.05 × 10<sup>-168</sup></td><td>1.40 × 10<sup>-79</sup></td><td>1.20 × 10<sup>-29</sup></td><td>4.83 × 10<sup>-39</sup></td><td>3.545085</td><td>41.78794</td><td>5.54 × 10<sup>-38</sup></td><td>1.14 × 10<sup>-68</sup></td><td>9.42 × 10<sup>-14</sup></td></tr><tr><td>Worst</td><td>9.88 × 10<sup>-168</sup></td><td>1.34 × 10<sup>-78</sup></td><td>1.19 × 10<sup>-28</sup></td><td>3.04 × 10<sup>-38</sup></td><td>6.640496</td><td>45.66198</td><td>2.52 × 10<sup>-37</sup></td><td>8.04 × 10<sup>-68</sup></td><td>7.15 × 10<sup>-13</sup></td></tr><tr><td>SD</td><td>0</td><td>4.22 × 10<sup>-79</sup></td><td>3.75 × 10<sup>-29</sup></td><td>9.54 × 10<sup>-39</sup></td><td>1.935587</td><td>3.361813</td><td>7.92 × 10<sup>-38</sup></td><td>2.58 × 10<sup>-68</sup></td><td>2.22 × 10<sup>-13</sup></td></tr><tr><td>F5</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>0</td><td>0</td><td>2.72 × 10<sup>1</sup></td><td>2.66 × 10<sup>1</sup></td><td>25.51494</td><td>1565934</td><td>2.78 × 10<sup>1</sup></td><td>2.83 × 10<sup>1</sup></td><td>1.76 × 10<sup>-1</sup></td></tr><tr><td>Average</td><td>0</td><td>0</td><td>2.84 × 10<sup>1</sup></td><td>2.74 × 10<sup>1</sup></td><td>29.92931</td><td>2104121</td><td>2.82 × 10<sup>1</sup></td><td>2.86 × 10<sup>1</sup></td><td>1.50</td></tr><tr><td>Worst</td><td>0</td><td>0</td><td>2.88 × 10<sup>1</sup></td><td>2.87 × 10<sup>1</sup></td><td>45.20891</td><td>3043646</td><td>2.87 × 10<sup>1</sup></td><td>2.89 × 10<sup>1</sup></td><td>6.61</td></tr><tr><td>SD</td><td>0</td><td>0</td><td>5.40 × 10<sup>-1</sup></td><td>6.14 × 10<sup>-1</sup></td><td>5.543385</td><td>420298.4</td><td>2.84 × 10<sup>-1</sup></td><td>1.82 × 10<sup>-1</sup></td><td>2.42</td></tr><tr><td>F6</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>0</td><td>0</td><td>1.49</td><td>3.43 × 10<sup>-2</sup></td><td>0.000109</td><td>2850.629</td><td>7.50</td><td>2.44</td><td>2.07 × 10<sup>-3</sup></td></tr><tr><td>Average</td><td>0</td><td>0</td><td>2.19</td><td>3.26 × 10<sup>-1</sup></td><td>0.000235</td><td>3806.061</td><td>7.50</td><td>2.99</td><td>1.85 × 10<sup>-2</sup></td></tr><tr><td>Worst</td><td>0</td><td>0</td><td>3.07</td><td>6.58 × 10<sup>-1</sup></td><td>0.000602</td><td>4724.147</td><td>7.50</td><td>3.72</td><td>4.01 × 10<sup>-2</sup></td></tr><tr><td>SD</td><td>0</td><td>0</td><td>6.00</td><td>2.37 × 10<sup>-1</sup></td><td>0.000138</td><td>607.6398</td><td>0</td><td>4.79 × 10<sup>-1</sup></td><td>1.27 × 10<sup>-2</sup></td></tr><tr><td>F7</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>9.92 × 10<sup>-7</sup></td><td>2.86 × 10<sup>-6</sup></td><td>8.83 × 10<sup>-7</sup></td><td>9.48 × 10<sup>-5</sup></td><td>0.018263</td><td>0.925058</td><td>9.82 × 10<sup>-6</sup></td><td>2.62 × 10<sup>-5</sup></td><td>2.06 × 10<sup>-4</sup></td></tr><tr><td>Average</td><td>5.88 × 10<sup>-5</sup></td><td>6.65 × 10<sup>-5</sup></td><td>5.40 × 10<sup>-4</sup></td><td>4.65 × 10<sup>-4</sup></td><td>0.04415</td><td>1.261838</td><td>2.01 × 10<sup>-4</sup></td><td>1.39 × 10<sup>-4</sup></td><td>9.75 × 10<sup>-4</sup></td></tr><tr><td>Worst</td><td>1.53 × 10<sup>-4</sup></td><td>1.86 × 10<sup>-4</sup></td><td>6.53 × 10<sup>-3</sup></td><td>1.16 × 10<sup>-3</sup></td><td>0.098175</td><td>1.897037</td><td>4.58 × 10<sup>-4</sup></td><td>3.49 × 10<sup>-4</sup></td><td>1.92 × 10<sup>-3</sup></td></tr><tr><td>SD</td><td>4.31 × 10<sup>-5</sup></td><td>5.29 × 10<sup>-5</sup></td><td>1.44 × 10<sup>-3</sup></td><td>3.00 × 10<sup>-4</sup></td><td>0.024073</td><td>0.268522</td><td>1.51 × 10<sup>-4</sup></td><td>9.20 × 10<sup>-5</sup></td><td>4.80 × 10<sup>-4</sup></td></tr><tr><td>F8</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>-1.26 × 10<sup>4</sup></td><td>-9.02 × 10<sup>3</sup></td><td>-7.41 × 10<sup>3</sup></td><td>-1.12 × 10<sup>4</sup></td><td>-8685.03</td><td>-11610.1</td><td>-4.72 × 10<sup>3</sup></td><td>-7.26 × 10<sup>3</sup></td><td>-1.05 × 10<sup>4</sup></td></tr><tr><td>Average</td><td>-1.26 × 10<sup>4</sup></td><td>-8.34 × 10<sup>3</sup></td><td>-6.43 × 10<sup>3</sup></td><td>-1.03 × 10<sup>4</sup></td><td>-7626.41</td><td>-11407.5</td><td>-3.71 × 10<sup>3</sup></td><td>-6.59 × 10<sup>3</sup></td><td>-9.49 × 10<sup>3</sup></td></tr><tr><td>Worst</td><td>-1.26 × 10<sup>4</sup></td><td>-7.52 × 10<sup>3</sup></td><td>-4.99 × 10<sup>3</sup></td><td>-8.99 × 10<sup>3</sup></td><td>-6129.31</td><td>-11106.4</td><td>-2.86 × 10<sup>3</sup></td><td>-5.85 × 10<sup>3</sup></td><td>-8.00 × 10<sup>3</sup></td></tr><tr><td>SD</td><td>7.52 × 10<sup>-10</sup></td><td>5.03 × 10<sup>2</sup></td><td>8.17 × 10<sup>2</sup></td><td>5.93 × 10<sup>2</sup></td><td>765.964</td><td>144.0269</td><td>5.46 × 10<sup>2</sup></td><td>5.75 × 10<sup>2</sup></td><td>6.43 × 10<sup>2</sup></td></tr><tr><td>F9</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>0</td><td>0</td><td>0</td><td>0</td><td>21.06603</td><td>50.57476</td><td>0</td><td>0</td><td>0</td></tr><tr><td>Average</td><td>0</td><td>0</td><td>0</td><td>0</td><td>40.21929</td><td>60.92757</td><td>9.60 × 10<sup>1</sup></td><td>0</td><td>0</td></tr><tr><td>Worst</td><td>0</td><td>0</td><td>0</td><td>0</td><td>59.40942</td><td>71.63848</td><td>2.58 × 10<sup>2</sup></td><td>0</td><td>0</td></tr><tr><td>SD</td><td>0</td><td>0</td><td>0</td><td>0</td><td>15.38431</td><td>6.732379</td><td>1.24 × 10<sup>2</sup></td><td>0</td><td>0</td></tr><tr><td>F10</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>8.88 × 10<sup>-16</sup></td><td>8.88 × 10<sup>-16</sup></td><td>8.88 × 10<sup>-16</sup></td><td>8.88 × 10<sup>-16</sup></td><td>0.003418</td><td>9.380315</td><td>4.44 × 10<sup>-15</sup></td><td>8.88 × 10<sup>-16</sup></td><td>8.88 × 10<sup>-16</sup></td></tr><tr><td>Average</td><td>8.88 × 10<sup>-16</sup></td><td>2.66 × 10<sup>-15</sup></td><td>8.88 × 10<sup>-16</sup></td><td>8.88 × 10<sup>-16</sup></td><td>0.006712</td><td>11.32595</td><td>4.44 × 10<sup>-15</sup></td><td>8.88 × 10<sup>-16</sup></td><td>8.88 × 10<sup>-16</sup></td></tr><tr><td>Worst</td><td>8.88 × 10<sup>-16</sup></td><td>4.44 × 10<sup>-15</sup></td><td>8.88 × 10<sup>-16</sup></td><td>8.88 × 10<sup>-16</sup></td><td>0.00951</td><td>12.76612</td><td>4.44 × 10<sup>-15</sup></td><td>8.88 × 10<sup>-16</sup></td><td>8.88 × 10<sup>-16</sup></td></tr><tr><td>SD</td><td>0</td><td>1.87 × 10<sup>-15</sup></td><td>0</td><td>0</td><td>0.00188</td><td>0.95227</td><td>0</td><td>0</td><td>0</td></tr><tr><td>F11</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.000995</td><td>27.86147</td><td>0</td><td>0</td><td>0</td></tr><tr><td>Average</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.016851</td><td>32.78121</td><td>1.61 × 10<sup>-3</sup></td><td>0</td><td>0</td></tr><tr><td>Worst</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.044359</td><td>47.68918</td><td>1.61 × 10<sup>-2</sup></td><td>0</td><td>0</td></tr><tr><td>SD</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.013777</td><td>6.341414</td><td>5.08 × 10<sup>-3</sup></td><td>0</td><td>0</td></tr><tr><td>F12</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>1.57 × 10<sup>-32</sup></td><td>1.57 × 10<sup>-32</sup></td><td>4.72 × 10<sup>-2</sup></td><td>2.13 × 10<sup>-3</sup></td><td>6.11 × 10<sup>-6</sup></td><td>176966.7</td><td>1.11</td><td>1.38 × 10<sup>-1</sup></td><td>2.74 × 10<sup>-4</sup></td></tr><tr><td>Average</td><td>1.57 × 10<sup>-32</sup></td><td>1.57 × 10<sup>-32</sup></td><td>1.27 × 10<sup>-1</sup></td><td>1.22 × 10<sup>-2</sup></td><td>0.324104</td><td>352773.6</td><td>1.44</td><td>2.38 × 10<sup>-1</sup></td><td>2.23 × 10<sup>-3</sup></td></tr><tr><td>Worst</td><td>1.57 × 10<sup>-32</sup></td><td>1.57 × 10<sup>-32</sup></td><td>3.54 × 10<sup>-1</sup></td><td>2.28 × 10<sup>-2</sup></td><td>2.799928</td><td>621526.7</td><td>1.67</td><td>3.64 × 10<sup>-1</sup></td><td>8.00 × 10<sup>-3</sup></td></tr><tr><td>SD</td><td>2.88 × 10<sup>-48</sup></td><td>2.88 × 10<sup>-48</sup></td><td>8.71 × 10<sup>-2</sup></td><td>6.39 × 10<sup>-3</sup></td><td>0.876357</td><td>142951.5</td><td>2.90 × 10<sup>-01</sup></td><td>7.64 × 10<sup>-2</sup></td><td>2.80 × 10<sup>-3</sup></td></tr><tr><td>F13</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>1.35 × 10<sup>-32</sup></td><td>1.35 × 10<sup>-32</sup></td><td>1.73</td><td>6.67 × 10<sup>-1</sup></td><td>4.38 × 10<sup>-5</sup></td><td>1360677</td><td>4.87 × 10<sup>-1</sup></td><td>1.83</td><td>9.98 × 10<sup>-4</sup></td></tr><tr><td>Average</td><td>1.35 × 10<sup>-32</sup></td><td>1.35 × 10<sup>-32</sup></td><td>2.29</td><td>2.20</td><td>0.006739</td><td>3046984</td><td>6.20 × 10<sup>-1</sup></td><td>2.25</td><td>2.60 × 10<sup>-2</sup></td></tr><tr><td>Worst</td><td>1.35 × 10<sup>-32</sup></td><td>1.35 × 10<sup>-32</sup></td><td>2.79</td><td>2.78</td><td>0.044155</td><td>6761426</td><td>7.69 × 10<sup>-1</sup></td><td>2.71</td><td>1.42 × 10<sup>-1</sup></td></tr><tr><td>SD</td><td>2.88 × 10<sup>-48</sup></td><td>2.88 × 10<sup>-48</sup></td><td>3.67 × 10<sup>-1</sup></td><td>5.95 × 10<sup>-1</sup></td><td>0.013908</td><td>1730831</td><td>1.01 × 10<sup>-1</sup></td><td>2.84 × 10<sup>-1</sup></td><td>4.17 × 10<sup>-2</sup></td></tr><tr><td>F14</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>9.98 × 10<sup>-1</sup></td><td>9.98 × 10<sup>-1</sup></td><td>9.98 × 10<sup>-1</sup></td><td>9.98 × 10<sup>-1</sup></td><td>9.98 × 10<sup>-1</sup></td><td>0.998004</td><td>2.02</td><td>9.98 × 10<sup>-1</sup></td><td>9.98 × 10<sup>-1</sup></td></tr><tr><td>Average</td><td>9.98 × 10<sup>-1</sup></td><td>9.98 × 10<sup>-1</sup></td><td>2.97</td><td>1.20</td><td>1.196809</td><td>0.998012</td><td>1.23 × 10<sup>1</sup></td><td>3.27</td><td>1.20</td></tr><tr><td>Worst</td><td>9.98 × 10<sup>-1</sup></td><td>9.98 × 10<sup>-1</sup></td><td>1.08 × 10<sup>-1</sup></td><td>2.98</td><td>1.992031</td><td>0.998071</td><td>1.83 × 10<sup>1</sup></td><td>6.90</td><td>2.98</td></tr><tr><td>SD</td><td>7.40 × 10<sup>-17</sup></td><td>0</td><td>2.91</td><td>6.27 × 10<sup>-1</sup></td><td>0.419119</td><td>2.11 × 10<sup>-5</sup></td><td>4.35</td><td>2.23</td><td>6.27 × 10<sup>-1</sup></td></tr></tbody></table>
<table><caption>Table 6. *Cont.*</caption><tbody><tr><td>F15</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>3.07 × 10<sup>-4</sup></td><td>3.07 × 10<sup>-4</sup></td><td>3.08 × 10<sup>-4</sup></td><td>3.07 × 10<sup>-4</sup></td><td>3.14 × 10<sup>-4</sup></td><td>0.000779</td><td>3.27 × 10<sup>-4</sup></td><td>3.08 × 10<sup>-4</sup></td><td>3.08 × 10<sup>-4</sup></td></tr><tr><td>Average</td><td>3.07 × 10<sup>-4</sup></td><td>3.07 × 10<sup>-4</sup></td><td>5.62 × 10<sup>-4</sup></td><td>3.08 × 10<sup>-4</sup></td><td>0.002615</td><td>0.001969</td><td>3.57 × 10<sup>-4</sup></td><td>2.44 × 10<sup>-3</sup></td><td>4.04 × 10<sup>-4</sup></td></tr><tr><td>Worst</td><td>3.07 × 10<sup>-4</sup></td><td>3.07 × 10<sup>-4</sup></td><td>1.38 × 10<sup>-3</sup></td><td>3.10 × 10<sup>-4</sup></td><td>0.020363</td><td>0.002977</td><td>4.08 × 10<sup>-4</sup></td><td>2.11 × 10<sup>-2</sup></td><td>7.7 × 10<sup>-4</sup></td></tr><tr><td>SD</td><td>1.45 × 10<sup>-15</sup></td><td>1.60 × 10<sup>-19</sup></td><td>4.04 × 10<sup>-4</sup></td><td>8.01 × 10<sup>-7</sup></td><td>0.006241</td><td>8.55 × 10<sup>-4</sup></td><td>2.88 × 10<sup>-5</sup></td><td>6.55 × 10<sup>-3</sup></td><td>1.56 × 10<sup>-4</sup></td></tr><tr><td>F16</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>-1.030</td><td>-1.030</td><td>-1.030</td><td>-1.030</td><td>-1.030</td><td>-1.03162</td><td>-1.030</td><td>-1.030</td><td>-1.030</td></tr><tr><td>Average</td><td>-1.030</td><td>-1.030</td><td>-1.030</td><td>-1.030</td><td>-1.03163</td><td>-1.03141</td><td>-1.030</td><td>-1.030</td><td>-1.030</td></tr><tr><td>Worst</td><td>-1.030</td><td>-1.030</td><td>-1.030</td><td>-1.030</td><td>-1.03163</td><td>-1.03106</td><td>-1</td><td>-1.030</td><td>-1.030</td></tr><tr><td>SD</td><td>0</td><td>7.40 × 10<sup>-17</sup></td><td>1.02 × 10<sup>-9</sup></td><td>7.35 × 10<sup>-15</sup></td><td>1.12 × 10<sup>-11</sup></td><td>2.02 × 10<sup>-04</sup></td><td>1.12 × 10<sup>-2</sup></td><td>9.95 × 10<sup>-10</sup></td><td>5.07 × 10<sup>-16</sup></td></tr><tr><td>F17</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td><td>0.397918</td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td></tr><tr><td>Average</td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td><td>0.397887</td><td>0.398031</td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td></tr><tr><td>Worst</td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td><td>0.397887</td><td>0.398279</td><td>4.00 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td></tr><tr><td>SD</td><td>0</td><td>0</td><td>6.43 × 10<sup>-8</sup></td><td>0</td><td>2.12 × 10<sup>-10</sup></td><td>1.05 × 10<sup>-4</sup></td><td>5.27 × 10<sup>-4</sup></td><td>1.62 × 10<sup>-8</sup></td><td>0</td></tr><tr><td>F18</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3.000039</td><td>3</td><td>3</td><td>3</td></tr><tr><td>Average</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3.00531</td><td>3.67 × 10<sup>1</sup></td><td>3</td><td>3</td></tr><tr><td>Worst</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3.017819</td><td>8.41 × 10<sup>1</sup></td><td>3</td><td>3</td></tr><tr><td>SD</td><td>4.44 × 10<sup>-16</sup></td><td>4.91 × 10<sup>-16</sup></td><td>3.81 × 10<sup>-5</sup></td><td>1.22 × 10<sup>-15</sup></td><td>8.17 × 10<sup>-8</sup></td><td>5.75 × 10<sup>-3</sup></td><td>3.54 × 10<sup>1</sup></td><td>1.72 × 10<sup>-5</sup></td><td>1.48 × 10<sup>-16</sup></td></tr><tr><td>F19</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>-3.86</td><td>-3.86</td><td>-3.86</td><td>-3.86</td><td>-3.86</td><td>-3.86278</td><td>-3.86</td><td>-3.86</td><td>-3.86</td></tr><tr><td>Average</td><td>-3.86</td><td>-3.86</td><td>-3.86</td><td>-3.86</td><td>-3.86278</td><td>-3.86274</td><td>-3.86</td><td>-3.86</td><td>-3.86</td></tr><tr><td>Worst</td><td>-3.86</td><td>-3.86</td><td>-3.85</td><td>-3.86</td><td>-3.86278</td><td>-3.86265</td><td>-3.86</td><td>-3.86</td><td>-3.86</td></tr><tr><td>SD</td><td>9.36 × 10<sup>-16</sup></td><td>8.63 × 10<sup>-16</sup></td><td>2.49 × 10<sup>-3</sup></td><td>7.40 × 10<sup>-16</sup></td><td>1.41 × 10<sup>-7</sup></td><td>4.33 × 10<sup>-5</sup></td><td>1.81 × 10<sup>-3</sup></td><td>1.32 × 10<sup>-3</sup></td><td>6.94 × 10<sup>-16</sup></td></tr><tr><td>F20</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>-3.32</td><td>-3.32</td><td>-3.32</td><td>-3.32</td><td>-3.32</td><td>-3.32197</td><td>-3.27</td><td>-3.32</td><td>-3.32</td></tr><tr><td>Average</td><td>-3.27</td><td>-3.27</td><td>-3.24</td><td>-3.30</td><td>-3.28631</td><td>-3.27421</td><td>-3.19</td><td>-3.32</td><td>-3.27</td></tr><tr><td>Worst</td><td>-3.20</td><td>-3.20</td><td>-3.02</td><td>-3.20</td><td>-3.20292</td><td>-3.20242</td><td>-3.08</td><td>-3.31</td><td>-3.20</td></tr><tr><td>SD</td><td>6.14 × 10<sup>-2</sup></td><td>6.14 × 10<sup>-2</sup></td><td>9.91 × 10<sup>-2</sup></td><td>5.01 × 10<sup>-2</sup></td><td>5.75 × 10<sup>-2</sup></td><td>6.15 × 10<sup>-2</sup></td><td>7.10 × 10<sup>-2</sup></td><td>4.96 × 10<sup>-3</sup></td><td>6.14 × 10<sup>-2</sup></td></tr><tr><td>F21</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>-1.02 × 10<sup>1</sup></td><td>-1.02 × 10<sup>1</sup></td><td>-1.02 × 10<sup>1</sup></td><td>-1.02 × 10<sup>1</sup></td><td>-1.02 × 10<sup>1</sup></td><td>-10.15</td><td>-6.7</td><td>-1.02 × 10<sup>1</sup></td><td>-1.02 × 10<sup>1</sup></td></tr><tr><td>Average</td><td>-1.02 × 10<sup>1</sup></td><td>-1.02 × 10<sup>1</sup></td><td>-5.59</td><td>-1.02 × 10<sup>1</sup></td><td>-4.91872</td><td>-4.90209</td><td>-4.99</td><td>-8.89</td><td>-9.40</td></tr><tr><td>Worst</td><td>-1.02 × 10<sup>1</sup></td><td>-1.02 × 10<sup>1</sup></td><td>-2.63</td><td>-1.01 × 10<sup>1</sup></td><td>-2.63047</td><td>-2.62769</td><td>-3.63</td><td>-2.63</td><td>-2.63</td></tr><tr><td>SD</td><td>1.03 × 10<sup>-15</sup></td><td>1.32 × 10<sup>-15</sup></td><td>2.60</td><td>6.52 × 10<sup>-3</sup></td><td>3.61</td><td>3.61</td><td>9.50 × 10<sup>-1</sup></td><td>2.72</td><td>2.38</td></tr><tr><td>F22</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>-1.04 × 10<sup>1</sup></td><td>-1.04 × 10<sup>1</sup></td><td>-1.04 × 10<sup>1</sup></td><td>-1.04 × 10<sup>1</sup></td><td>-1.04 × 10<sup>1</sup></td><td>-10.3979</td><td>-7.79</td><td>-1.04 × 10<sup>1</sup></td><td>-1.04 × 10<sup>1</sup></td></tr><tr><td>Average</td><td>-1.04 × 10<sup>1</sup></td><td>-9.34</td><td>-5.29</td><td>-9.87</td><td>-5.91376</td><td>-6.21638</td><td>-5.47</td><td>-9.87</td><td>-8.97</td></tr><tr><td>Worst</td><td>-1.04 × 10<sup>1</sup></td><td>-5.09</td><td>-1.84</td><td>-5.09</td><td>-2.75193</td><td>-2.7657</td><td>-2.26</td><td>-5.09</td><td>-2.77</td></tr><tr><td>SD</td><td>1.78 × 10<sup>-15</sup></td><td>2.24</td><td>2.07</td><td>1.68</td><td>3.87</td><td>3.58</td><td>1.89</td><td>1.68</td><td>3.03</td></tr><tr><td>F23</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>-1.050 × 10<sup>1</sup></td><td>-1.050 × 10<sup>1</sup></td><td>-1.050 × 10<sup>1</sup></td><td>-1.050 × 10<sup>1</sup></td><td>-1.050 × 10<sup>1</sup></td><td>-10.5328</td><td>-8.36</td><td>-1.050 × 10<sup>1</sup></td><td>-1.050 × 10<sup>1</sup></td></tr><tr><td>Average</td><td>-1.050 × 10<sup>1</sup></td><td>-9.45</td><td>-8.24</td><td>-1.050 × 10<sup>1</sup></td><td>-7.2443</td><td>-7.45963</td><td>-6.05</td><td>-8.91</td><td>-9.87</td></tr><tr><td>Worst</td><td>-1.050 × 10<sup>1</sup></td><td>-5.13</td><td>-3.84</td><td>-1.050 × 10<sup>-1</sup></td><td>-1.85948</td><td>-2.42579</td><td>-3.22</td><td>-5.13</td><td>-3.84</td></tr><tr><td>SD</td><td>1.78 × 10<sup>-15</sup></td><td>2.28</td><td>2.98</td><td>5.37 × 10<sup>-6</sup></td><td>3.64</td><td>3.95</td><td>1.65</td><td>2.61</td><td>2.12</td></tr></tbody></table>

Table 6 shows that for the majority of the functions under study, the suggested I-WaOA optimizer performs better than the others in terms of worst, best, and mean values. This means that the I-WaOA optimizer was able to find better solutions on average, as well as the best and worst solutions, than the other optimizers. All algorithms have a community size of 30 and a maximum iteration of 300.

### 7.1.2. Convergence Curve Analysis

The I-WaOA algorithm's convergence graphs and those of other previously published techniques, including AHA, SCSO, DO, HS, CDO, ZOA, ARO, and standard WaOA, can be observed in Figure 6.

The proposed I-WaOA algorithm was shown to be more accurate and faster than other algorithms based on convergence curves. The proposed modifications to the WaOA algorithm enhance both the exploration and exploitation phases of the algorithm, which allows it to converge to an optimal solution more quickly than the traditional WaOA algorithm.
<table><caption>Table 6. *Cont.*</caption><tbody><tr><td>F15</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>3.07 × 10<sup>-4</sup></td><td>3.07 × 10<sup>-4</sup></td><td>3.08 × 10<sup>-4</sup></td><td>3.07 × 10<sup>-4</sup></td><td>3.14 × 10<sup>-4</sup></td><td>0.000779</td><td>3.27 × 10<sup>-4</sup></td><td>3.08 × 10<sup>-4</sup></td><td>3.08 × 10<sup>-4</sup></td></tr><tr><td>Average</td><td>3.07 × 10<sup>-4</sup></td><td>3.07 × 10<sup>-4</sup></td><td>5.62 × 10<sup>-4</sup></td><td>3.08 × 10<sup>-4</sup></td><td>0.002615</td><td>0.001969</td><td>3.57 × 10<sup>-4</sup></td><td>2.44 × 10<sup>-3</sup></td><td>4.04 × 10<sup>-4</sup></td></tr><tr><td>Worst</td><td>3.07 × 10<sup>-4</sup></td><td>3.07 × 10<sup>-4</sup></td><td>1.38 × 10<sup>-3</sup></td><td>3.10 × 10<sup>-4</sup></td><td>0.020363</td><td>0.002977</td><td>4.08 × 10<sup>-4</sup></td><td>2.11 × 10<sup>-2</sup></td><td>7.7 × 10<sup>-4</sup></td></tr><tr><td>SD</td><td>1.45 × 10<sup>-15</sup></td><td>1.60 × 10<sup>-19</sup></td><td>4.04 × 10<sup>-4</sup></td><td>8.01 × 10<sup>-7</sup></td><td>0.006241</td><td>8.55 × 10<sup>-4</sup></td><td>2.88 × 10<sup>-5</sup></td><td>6.55 × 10<sup>-3</sup></td><td>1.56 × 10<sup>-4</sup></td></tr><tr><td>F16</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>-1.030</td><td>-1.030</td><td>-1.030</td><td>-1.030</td><td>-1.030</td><td>-1.03162</td><td>-1.030</td><td>-1.030</td><td>-1.030</td></tr><tr><td>Average</td><td>-1.030</td><td>-1.030</td><td>-1.030</td><td>-1.030</td><td>-1.03163</td><td>-1.03141</td><td>-1.030</td><td>-1.030</td><td>-1.030</td></tr><tr><td>Worst</td><td>-1.030</td><td>-1.030</td><td>-1.030</td><td>-1.030</td><td>-1.03163</td><td>-1.03106</td><td>-1</td><td>-1.030</td><td>-1.030</td></tr><tr><td>SD</td><td>0</td><td>7.40 × 10<sup>-17</sup></td><td>1.02 × 10<sup>-9</sup></td><td>7.35 × 10<sup>-15</sup></td><td>1.12 × 10<sup>-11</sup></td><td>2.02 × 10<sup>-04</sup></td><td>1.12 × 10<sup>-2</sup></td><td>9.95 × 10<sup>-10</sup></td><td>5.07 × 10<sup>-16</sup></td></tr><tr><td>F17</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td><td>0.397918</td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td></tr><tr><td>Average</td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td><td>0.397887</td><td>0.398031</td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td></tr><tr><td>Worst</td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td><td>0.397887</td><td>0.398279</td><td>4.00 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td><td>3.9888 × 10<sup>-1</sup></td></tr><tr><td>SD</td><td>0</td><td>0</td><td>6.43 × 10<sup>-8</sup></td><td>0</td><td>2.12 × 10<sup>-10</sup></td><td>1.05 × 10<sup>-4</sup></td><td>5.27 × 10<sup>-4</sup></td><td>1.62 × 10<sup>-8</sup></td><td>0</td></tr><tr><td>F18</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3.000039</td><td>3</td><td>3</td><td>3</td></tr><tr><td>Average</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3.00531</td><td>3.67 × 10<sup>1</sup></td><td>3</td><td>3</td></tr><tr><td>Worst</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3.017819</td><td>8.41 × 10<sup>1</sup></td><td>3</td><td>3</td></tr><tr><td>SD</td><td>4.44 × 10<sup>-16</sup></td><td>4.91 × 10<sup>-16</sup></td><td>3.81 × 10<sup>-5</sup></td><td>1.22 × 10<sup>-15</sup></td><td>8.17 × 10<sup>-8</sup></td><td>5.75 × 10<sup>-3</sup></td><td>3.54 × 10<sup>1</sup></td><td>1.72 × 10<sup>-5</sup></td><td>1.48 × 10<sup>-16</sup></td></tr><tr><td>F19</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>-3.86</td><td>-3.86</td><td>-3.86</td><td>-3.86</td><td>-3.86</td><td>-3.86278</td><td>-3.86</td><td>-3.86</td><td>-3.86</td></tr><tr><td>Average</td><td>-3.86</td><td>-3.86</td><td>-3.86</td><td>-3.86</td><td>-3.86278</td><td>-3.86274</td><td>-3.86</td><td>-3.86</td><td>-3.86</td></tr><tr><td>Worst</td><td>-3.86</td><td>-3.86</td><td>-3.85</td><td>-3.86</td><td>-3.86278</td><td>-3.86265</td><td>-3.86</td><td>-3.86</td><td>-3.86</td></tr><tr><td>SD</td><td>9.36 × 10<sup>-16</sup></td><td>8.63 × 10<sup>-16</sup></td><td>2.49 × 10<sup>-3</sup></td><td>7.40 × 10<sup>-16</sup></td><td>1.41 × 10<sup>-7</sup></td><td>4.33 × 10<sup>-5</sup></td><td>1.81 × 10<sup>-3</sup></td><td>1.32 × 10<sup>-3</sup></td><td>6.94 × 10<sup>-16</sup></td></tr><tr><td>F20</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>-3.32</td><td>-3.32</td><td>-3.32</td><td>-3.32</td><td>-3.32</td><td>-3.32197</td><td>-3.27</td><td>-3.32</td><td>-3.32</td></tr><tr><td>Average</td><td>-3.27</td><td>-3.27</td><td>-3.24</td><td>-3.30</td><td>-3.28631</td><td>-3.27421</td><td>-3.19</td><td>-3.32</td><td>-3.27</td></tr><tr><td>Worst</td><td>-3.20</td><td>-3.20</td><td>-3.02</td><td>-3.20</td><td>-3.20292</td><td>-3.20242</td><td>-3.08</td><td>-3.31</td><td>-3.20</td></tr><tr><td>SD</td><td>6.14 × 10<sup>-2</sup></td><td>6.14 × 10<sup>-2</sup></td><td>9.91 × 10<sup>-2</sup></td><td>5.01 × 10<sup>-2</sup></td><td>5.75 × 10<sup>-2</sup></td><td>6.15 × 10<sup>-2</sup></td><td>7.10 × 10<sup>-2</sup></td><td>4.96 × 10<sup>-3</sup></td><td>6.14 × 10<sup>-2</sup></td></tr><tr><td>F21</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>-1.02 × 10<sup>1</sup></td><td>-1.02 × 10<sup>1</sup></td><td>-1.02 × 10<sup>1</sup></td><td>-1.02 × 10<sup>1</sup></td><td>-1.02 × 10<sup>1</sup></td><td>-10.15</td><td>-6.7</td><td>-1.02 × 10<sup>1</sup></td><td>-1.02 × 10<sup>1</sup></td></tr><tr><td>Average</td><td>-1.02 × 10<sup>1</sup></td><td>-1.02 × 10<sup>1</sup></td><td>-5.59</td><td>-1.02 × 10<sup>1</sup></td><td>-4.91872</td><td>-4.90209</td><td>-4.99</td><td>-8.89</td><td>-9.40</td></tr><tr><td>Worst</td><td>-1.02 × 10<sup>1</sup></td><td>-1.02 × 10<sup>1</sup></td><td>-2.63</td><td>-1.01 × 10<sup>1</sup></td><td>-2.63047</td><td>-2.62769</td><td>-3.63</td><td>-2.63</td><td>-2.63</td></tr><tr><td>SD</td><td>1.03 × 10<sup>-15</sup></td><td>1.32 × 10<sup>-15</sup></td><td>2.60</td><td>6.52 × 10<sup>-3</sup></td><td>3.61</td><td>3.61</td><td>9.50 × 10<sup>-1</sup></td><td>2.72</td><td>2.38</td></tr><tr><td>F22</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>-1.04 × 10<sup>1</sup></td><td>-1.04 × 10<sup>1</sup></td><td>-1.04 × 10<sup>1</sup></td><td>-1.04 × 10<sup>1</sup></td><td>-1.04 × 10<sup>1</sup></td><td>-10.3979</td><td>-7.79</td><td>-1.04 × 10<sup>1</sup></td><td>-1.04 × 10<sup>1</sup></td></tr><tr><td>Average</td><td>-1.04 × 10<sup>1</sup></td><td>-9.34</td><td>-5.29</td><td>-9.87</td><td>-5.91376</td><td>-6.21638</td><td>-5.47</td><td>-9.87</td><td>-8.97</td></tr><tr><td>Worst</td><td>-1.04 × 10<sup>1</sup></td><td>-5.09</td><td>-1.84</td><td>-5.09</td><td>-2.75193</td><td>-2.7657</td><td>-2.26</td><td>-5.09</td><td>-2.77</td></tr><tr><td>SD</td><td>1.78 × 10<sup>-15</sup></td><td>2.24</td><td>2.07</td><td>1.68</td><td>3.87</td><td>3.58</td><td>1.89</td><td>1.68</td><td>3.03</td></tr><tr><td>F23</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Best</td><td>-1.050 × 10<sup>1</sup></td><td>-1.050 × 10<sup>1</sup></td><td>-1.050 × 10<sup>1</sup></td><td>-1.050 × 10<sup>1</sup></td><td>-1.050 × 10<sup>1</sup></td><td>-10.5328</td><td>-8.36</td><td>-1.050 × 10<sup>1</sup></td><td>-1.050 × 10<sup>1</sup></td></tr><tr><td>Average</td><td>-1.050 × 10<sup>1</sup></td><td>-9.45</td><td>-8.24</td><td>-1.050 × 10<sup>1</sup></td><td>-7.2443</td><td>-7.45963</td><td>-6.05</td><td>-8.91</td><td>-9.87</td></tr><tr><td>Worst</td><td>-1.050 × 10<sup>1</sup></td><td>-5.13</td><td>-3.84</td><td>-1.050 × 10<sup>-1</sup></td><td>-1.85948</td><td>-2.42579</td><td>-3.22</td><td>-5.13</td><td>-3.84</td></tr><tr><td>SD</td><td>1.78 × 10<sup>-15</sup></td><td>2.28</td><td>2.98</td><td>5.37 × 10<sup>-6</sup></td><td>3.64</td><td>3.95</td><td>1.65</td><td>2.61</td><td>2.12</td></tr></tbody></table>

Table 6 shows that for the majority of the functions under study, the suggested I-WaOA optimizer performs better than the others in terms of worst, best, and mean values. This means that the I-WaOA optimizer was able to find better solutions on average, as well as the best and worst solutions, than the other optimizers. All algorithms have a community size of 30 and a maximum iteration of 300.

### 7.1.2. Convergence Curve Analysis

The I-WaOA algorithm's convergence graphs and those of other previously published techniques, including AHA, SCSO, DO, HS, CDO, ZOA, ARO, and standard WaOA, can be observed in Figure 6.

The proposed I-WaOA algorithm was shown to be more accurate and faster than other algorithms based on convergence curves. The proposed modifications to the WaOA algorithm enhance both the exploration and exploitation phases of the algorithm, which allows it to converge to an optimal solution more quickly than the traditional WaOA algorithm.
<figure><img src="image_9.png" /><figcaption>Figure 6. *Cont.*</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>Figure 6. *Cont.*</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>Figure 6. Convergence analysis of various optimizers on test benchmark functions.</figcaption></figure>

### 7.1.3. Analysis Using Boxplots

Boxplot graphs effectively display data distribution patterns in four equal groups, enabling examination of data distribution features. Figure 7 shows the boxplots for the proposed I-WaOA method and other published optimization techniques. The boxplots for I-WaOA visibly display a smaller spread compared to the other optimization techniques.

<figure><img src="image_11.png" /><figcaption>Figure 7. Cont.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**Figure 6.** Convergence analysis of various optimizers on test benchmark functions.</figcaption></figure>

### 7.1.3. Analysis Using Boxplots

Boxplot graphs effectively display data distribution patterns in four equal groups, enabling examination of data distribution features. Figure 7 shows the boxplots for the proposed I-WaOA method and other published optimization techniques. The boxplots for I-WaOA visibly display a smaller spread compared to the other optimization techniques.

<figure><img src="image_11.png" /><figcaption>**Figure 7.** *Cont.*</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>Figure 7. *Cont.*</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>Figure 7. *Cont.*</figcaption></figure>
<figure><img src="image_13.png" /><figcaption>**Figure 7.** Boxplot analysis of test benchmark functions using various algorithms.</figcaption></figure>

## 7.2. Applying the Suggested Approach to Solve the IEEE 118 Distribution Network's Optimal Operating Problem

In this section, we implement the I-WaOA approach to address the OOP and ascertain the optimal placement of three systems composed of one PV unit rating and one WT in each system within the IEEE 118 distribution network. The configuration, depicted in Figure 8, encompasses 118 buses and 117 branches. Pertinent data are provided in [79]. Meteorological information, including wind speed and solar radiation, along with temperature data, is sourced from [7,80], respectively. To verify the suggested I-WaOA algorithm's effectiveness, we conduct a comparative analysis between its results and those obtained from the original method. For a meaningful and unbiased comparison, the population size and maximum iteration count of the suggested methods are fixed at 25 and 60, respectively. Table 7 lists the operational restrictions together with the RESs price variables.

<figure><img src="image_14.png" /><figcaption>**Figure 8.** Schematic diagram of the IEEE 118 distribution network.</figcaption></figure>
<figure><img src="image_13.png" /><figcaption>**Figure 7.** Boxplot analysis of test benchmark functions using various algorithms.</figcaption></figure>

## 7.2. Applying the Suggested Approach to Solve the IEEE 118 Distribution Network's Optimal Operating Problem

In this section, we implement the I-WaOA approach to address the OOP and ascertain the optimal placement of three systems composed of one PV unit rating and one WT in each system within the IEEE 118 distribution network. The configuration, depicted in Figure 8, encompasses 118 buses and 117 branches. Pertinent data are provided in [79]. Meteorological information, including wind speed and solar radiation, along with temperature data, is sourced from [7,80], respectively. To verify the suggested I-WaOA algorithm's effectiveness, we conduct a comparative analysis between its results and those obtained from the original method. For a meaningful and unbiased comparison, the population size and maximum iteration count of the suggested methods are fixed at 25 and 60, respectively. Table 7 lists the operational restrictions together with the RESs price variables.

<figure><img src="image_14.png" /><figcaption>**Figure 8.** Schematic diagram of the IEEE 118 distribution network.</figcaption></figure>
<table><caption>**Table 7.** The limitations and the coefficients of expenses.</caption><thead><tr><th>Parameter</th><th>Value</th></tr></thead><tbody><tr><td>WT cost [81]</td><td></td></tr><tr><td>The investment cost (U<sub>WT</sub>)</td><td>1400 USD/kW</td></tr><tr><td>The maintenance and operation costs (U<sub>WT</sub><sup>O&amp;M</sup>)</td><td>0.01 USD/kWh</td></tr><tr><td>The interest rate (β<sub>WT</sub>)</td><td>10%</td></tr><tr><td>The lifetime (NP<sub>WT</sub>)</td><td>20</td></tr><tr><td>PV cost [82]</td><td></td></tr><tr><td>The investment cost (U<sub>PV</sub>)</td><td>770 USD/kW</td></tr><tr><td>The maintenance and operation costs (U<sub>PV</sub><sup>O&amp;M</sup>)</td><td>0.01 USD/kWh</td></tr><tr><td>The interest rate (β<sub>PV</sub>)</td><td>10%</td></tr><tr><td>The lifetime (NP<sub>PV</sub>)</td><td>20</td></tr><tr><td>Cost coefficients [83]</td><td></td></tr><tr><td>The energy loss cost (U<sub>Loss</sub>)</td><td>0.06 USD/kWh</td></tr><tr><td>Constraints of grid and generators</td><td></td></tr><tr><td>Voltage boundaries [84]</td><td>0.9 p.u ≤ V ≤ 1.1 p.u</td></tr><tr><td>Area sizes</td><td>0 ≤ area ≤ 89057 m<sup>2</sup></td></tr><tr><td>WT sizes</td><td>0 ≤ WT ≤ 90 turbines</td></tr><tr><td>Power factor of the PV</td><td>1</td></tr><tr><td>Power factor of the WT</td><td>0.7 ≤ PF ≤ 1</td></tr></tbody></table>

The methodology presented in Section 4 took into account a variety of uncertainties, encompassing variations in electricity price, wind speed, temperature, solar irradiation, and load. Figure 9 shows the anticipated day-ahead load demand, market energy purchasing prices, wind speed, irradiance, and temperature patterns.

The results obtained from conducting case studies on the 118 bus system are presented in Table 8. This table provides information about the impact of integrating PVs and WTs both with and without Demand Side Response on the system. Overall, it offers insights into how PVs, WTs, and DSR can enhance the operation of distribution grids that utilize energy sources.

<table><caption>**Table 8.** The IEEE 118 distribution network's energy management results.</caption><thead><tr><th rowspan="2">Item</th><th rowspan="2">Without RESs and DSR</th><th colspan="2">With RESs Only</th><th colspan="2">With RESs and DSR</th></tr><tr><th>WAOA</th><th>I-WAOA</th><th>WAOA</th><th>I-WAOA</th></tr></thead><tbody><tr><td>Energy losses (kWh)</td><td>1.9613 × 10<sup>4</sup></td><td>1.8894 × 10<sup>4</sup></td><td>1.4501 × 10<sup>4</sup></td><td>7.8651 × 10<sup>3</sup></td><td>7.8020 × 10<sup>3</sup></td></tr><tr><td>Purchased power from the grid (kW)</td><td>4.5347 × 10<sup>5</sup></td><td>1.2795 × 10<sup>5</sup></td><td>1.3920 × 10<sup>5</sup></td><td>3.8089 × 10<sup>4</sup></td><td>5.8284 × 10<sup>4</sup></td></tr><tr><td rowspan="3">Optimal location of systems</td><td></td><td>2</td><td>104</td><td>38</td><td>66</td></tr><tr><td></td><td>67</td><td>68</td><td>36</td><td>31</td></tr><tr><td></td><td>7</td><td>35</td><td>69</td><td>26</td></tr><tr><td rowspan="3">Optimal area of the solar modules (m<sup>2</sup>)</td><td></td><td>3.2911 × 10<sup>3</sup></td><td>4.8552 × 10<sup>3</sup></td><td>5.4034 × 10<sup>4</sup></td><td>1.1217 × 10<sup>4</sup></td></tr><tr><td></td><td>7.0174 × 10<sup>3</sup></td><td>9.9422 × 10<sup>3</sup></td><td>7.4928 × 10<sup>3</sup></td><td>7.7015 × 10<sup>4</sup></td></tr><tr><td></td><td>4.2454 × 10<sup>3</sup></td><td>5.1878 × 10<sup>3</sup></td><td>1.5157 × 10<sup>4</sup></td><td>4242</td></tr><tr><td rowspan="3">Optimal-size WTs (kW)</td><td></td><td>1000</td><td>8500</td><td>2000</td><td>6250</td></tr><tr><td></td><td>9250</td><td>6500</td><td>5750</td><td>1500</td></tr><tr><td></td><td>10250</td><td>4250</td><td>2250</td><td>250</td></tr><tr><td rowspan="3">Optimal PF of WTs</td><td></td><td>0.8850</td><td>0.8378</td><td>0.7524</td><td>0.7566</td></tr><tr><td></td><td>0.8477</td><td>0.8428</td><td>0.7000</td><td>0.7351</td></tr><tr><td></td><td>0.8271</td><td>0.7156</td><td>0.7000</td><td>0.7113</td></tr><tr><td>Total annual energy loss cost (USD)</td><td>4.2952 × 10<sup>5</sup></td><td>4.1378 × 10<sup>5</sup></td><td>3.1758 × 10<sup>5</sup></td><td>1.7225 × 10<sup>5</sup></td><td>1.7086 × 10<sup>5</sup></td></tr><tr><td>Total annual purchased energy cost (USD)</td><td>3.8377 × 10<sup>7</sup></td><td>1.2646 × 10<sup>7</sup></td><td>1.3344 × 10<sup>7</sup></td><td>1.8089 × 10<sup>3</sup></td><td>4.6348 × 10<sup>4</sup></td></tr><tr><td>Total annual RESs cost (USD)</td><td></td><td>4.7564 × 10<sup>6</sup></td><td>4.5683 × 10<sup>6</sup></td><td>3.4428 × 10<sup>6</sup></td><td>3.2565 × 10<sup>6</sup></td></tr><tr><td>TVD(p.u)</td><td>98.6633</td><td>66.1004</td><td>49.9790</td><td>35.7749</td><td>36.0990</td></tr><tr><td>TVSI(p.u)</td><td>2.4440 × 10<sup>3</sup></td><td>2.5926 × 10<sup>3</sup></td><td>2.6965 × 10<sup>3</sup></td><td>2.7358 × 10<sup>3</sup></td><td>2.7245 × 10<sup>3</sup></td></tr><tr><td>Total annual cost (USD)</td><td>3.8806 × 10<sup>7</sup></td><td>1.7817 × 10<sup>7</sup></td><td>1.8230 × 10<sup>7</sup></td><td>3.6169 × 10<sup>6</sup></td><td>3.4737 × 10<sup>6</sup></td></tr><tr><td>Best MOF</td><td></td><td>0.3971</td><td>0.3616</td><td>0.1373</td><td>0.1363</td></tr></tbody></table>
<table><caption>**Table 7.** The limitations and the coefficients of expenses.</caption><thead><tr><th>Parameter</th><th>Value</th></tr></thead><tbody><tr><td>WT cost [81]</td><td></td></tr><tr><td>The investment cost (U<sub>WT</sub>)</td><td>1400 USD/kW</td></tr><tr><td>The maintenance and operation costs (U<sub>WT</sub><sup>O&amp;M</sup>)</td><td>0.01 USD/kWh</td></tr><tr><td>The interest rate (β<sub>WT</sub>)</td><td>10%</td></tr><tr><td>The lifetime (NP<sub>WT</sub>)</td><td>20</td></tr><tr><td>PV cost [82]</td><td></td></tr><tr><td>The investment cost (U<sub>PV</sub>)</td><td>770 USD/kW</td></tr><tr><td>The maintenance and operation costs (U<sub>PV</sub><sup>O&amp;M</sup>)</td><td>0.01 USD/kWh</td></tr><tr><td>The interest rate (β<sub>PV</sub>)</td><td>10%</td></tr><tr><td>The lifetime (NP<sub>PV</sub>)</td><td>20</td></tr><tr><td>Cost coefficients [83]</td><td></td></tr><tr><td>The energy loss cost (U<sub>Loss</sub>)</td><td>0.06 USD/kWh</td></tr><tr><td>Constraints of grid and generators</td><td></td></tr><tr><td>Voltage boundaries [84]</td><td>0.9 p.u ≤ V ≤ 1.1 p.u</td></tr><tr><td>Area sizes</td><td>0 ≤ area ≤ 89057 m<sup>2</sup></td></tr><tr><td>WT sizes</td><td>0 ≤ WT ≤ 90 turbines</td></tr><tr><td>Power factor of the PV</td><td>1</td></tr><tr><td>Power factor of the WT</td><td>0.7 ≤ PF ≤ 1</td></tr></tbody></table>

The methodology presented in Section 4 took into account a variety of uncertainties, encompassing variations in electricity price, wind speed, temperature, solar irradiation, and load. Figure 9 shows the anticipated day-ahead load demand, market energy purchasing prices, wind speed, irradiance, and temperature patterns.

The results obtained from conducting case studies on the 118 bus system are presented in Table 8. This table provides information about the impact of integrating PVs and WTs both with and without Demand Side Response on the system. Overall, it offers insights into how PVs, WTs, and DSR can enhance the operation of distribution grids that utilize energy sources.

<table><caption>**Table 8.** The IEEE 118 distribution network's energy management results.</caption><thead><tr><th rowspan="2">Item</th><th rowspan="2">Without RESs and DSR</th><th colspan="2">With RESs Only</th><th colspan="2">With RESs and DSR</th></tr><tr><th>WAOA</th><th>I-WAOA</th><th>WAOA</th><th>I-WAOA</th></tr></thead><tbody><tr><td>Energy losses (kWh)</td><td>1.9613 × 10<sup>4</sup></td><td>1.8894 × 10<sup>4</sup></td><td>1.4501 × 10<sup>4</sup></td><td>7.8651 × 10<sup>3</sup></td><td>7.8020 × 10<sup>3</sup></td></tr><tr><td>Purchased power from the grid (kW)</td><td>4.5347 × 10<sup>5</sup></td><td>1.2795 × 10<sup>5</sup></td><td>1.3920 × 10<sup>5</sup></td><td>3.8089 × 10<sup>4</sup></td><td>5.8284 × 10<sup>4</sup></td></tr><tr><td rowspan="3">Optimal location of systems</td><td></td><td>2</td><td>104</td><td>38</td><td>66</td></tr><tr><td></td><td>67</td><td>68</td><td>36</td><td>31</td></tr><tr><td></td><td>7</td><td>35</td><td>69</td><td>26</td></tr><tr><td rowspan="3">Optimal area of the solar modules (m<sup>2</sup>)</td><td></td><td>3.2911 × 10<sup>3</sup></td><td>4.8552 × 10<sup>3</sup></td><td>5.4034 × 10<sup>4</sup></td><td>1.1217 × 10<sup>4</sup></td></tr><tr><td></td><td>7.0174 × 10<sup>3</sup></td><td>9.9422 × 10<sup>3</sup></td><td>7.4928 × 10<sup>3</sup></td><td>7.7015 × 10<sup>4</sup></td></tr><tr><td></td><td>4.2454 × 10<sup>3</sup></td><td>5.1878 × 10<sup>3</sup></td><td>1.5157 × 10<sup>4</sup></td><td>4242</td></tr><tr><td rowspan="3">Optimal-size WTs (kW)</td><td></td><td>1000</td><td>8500</td><td>2000</td><td>6250</td></tr><tr><td></td><td>9250</td><td>6500</td><td>5750</td><td>1500</td></tr><tr><td></td><td>10250</td><td>4250</td><td>2250</td><td>250</td></tr><tr><td rowspan="3">Optimal PF of WTs</td><td></td><td>0.8850</td><td>0.8378</td><td>0.7524</td><td>0.7566</td></tr><tr><td></td><td>0.8477</td><td>0.8428</td><td>0.7000</td><td>0.7351</td></tr><tr><td></td><td>0.8271</td><td>0.7156</td><td>0.7000</td><td>0.7113</td></tr><tr><td>Total annual energy loss cost (USD)</td><td>4.2952 × 10<sup>5</sup></td><td>4.1378 × 10<sup>5</sup></td><td>3.1758 × 10<sup>5</sup></td><td>1.7225 × 10<sup>5</sup></td><td>1.7086 × 10<sup>5</sup></td></tr><tr><td>Total annual purchased energy cost (USD)</td><td>3.8377 × 10<sup>7</sup></td><td>1.2646 × 10<sup>7</sup></td><td>1.3344 × 10<sup>7</sup></td><td>1.8089 × 10<sup>3</sup></td><td>4.6348 × 10<sup>4</sup></td></tr><tr><td>Total annual RESs cost (USD)</td><td></td><td>4.7564 × 10<sup>6</sup></td><td>4.5683 × 10<sup>6</sup></td><td>3.4428 × 10<sup>6</sup></td><td>3.2565 × 10<sup>6</sup></td></tr><tr><td>TVD(p.u)</td><td>98.6633</td><td>66.1004</td><td>49.9790</td><td>35.7749</td><td>36.0990</td></tr><tr><td>TVSI(p.u)</td><td>2.4440 × 10<sup>3</sup></td><td>2.5926 × 10<sup>3</sup></td><td>2.6965 × 10<sup>3</sup></td><td>2.7358 × 10<sup>3</sup></td><td>2.7245 × 10<sup>3</sup></td></tr><tr><td>Total annual cost (USD)</td><td>3.8806 × 10<sup>7</sup></td><td>1.7817 × 10<sup>7</sup></td><td>1.8230 × 10<sup>7</sup></td><td>3.6169 × 10<sup>6</sup></td><td>3.4737 × 10<sup>6</sup></td></tr><tr><td>Best MOF</td><td></td><td>0.3971</td><td>0.3616</td><td>0.1373</td><td>0.1363</td></tr></tbody></table>
<figure><img src="image_15.png" /><figcaption>**Figure 9.** Displays of (a) the loading, (b) the solar irradiance, (c) temperature, (d) wind speed, and (e) market price for purchasing energy.</figcaption></figure>

The primary goal of this study was to optimize expenses while enhancing the performance of the system by incorporating renewable energy resources and Demand Side Response into an electrical DN. The findings indicated that the proposed I-WaOA algorithm proved effective in reducing costs, improving voltage stability, and minimizing voltage deviations. In the scenario where no RESs were integrated into the DN, the total annual cost amounted to 3.8806 × 10⁷ USD. However, by using the suggested I-WaOA algorithm and strategically allocating RES, this cost is significantly decreased to 1.8230 × 10⁷ USD, resulting in a reduction in voltage deviations from 98.6633 per unit (p.u.) to 49.9790 p.u. Moreover, there was an enhancement in Voltage
<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /></figure>

<figure><img src="image_18.png" /></figure>

<figure><img src="image_19.png" /><figcaption>**Figure 9.** Displays of (a) the loading, (b) the solar irradiance, (c) temperature, (d) wind speed, and (e) market price for purchasing energy.</figcaption></figure>

The primary goal of this study was to optimize expenses while enhancing the performance of the system by incorporating renewable energy resources and Demand Side Response into an electrical DN. The findings indicated that the proposed I-WaOA algorithm proved effective in reducing costs, improving voltage stability, and minimizing voltage deviations. In the scenario where no RESs were integrated into the DN, the total annual cost amounted to 3.8806 × 10⁷ USD. However, by using the suggested I-WaOA algorithm and strategically allocating RES, this cost is significantly decreased to 1.8230 × 10⁷ USD, resulting in a reduction in voltage deviations from 98.6633 per unit (p.u.) to 49.9790 p.u. Moreover, there was an enhancement in Voltage
Stability Index, increasing from 2.4440 × 10³ p.u. to 2.6965 × 10³ p.u. The best places for the systems in the DN are 35, 68, and 104, respectively. The WTs ratings are 4250 kW, 6500 kW, and 8500 kW, respectively; in contrast, the solar modules' comparable areas are 5.1878 × 10³ m², 9.9422 × 10³ m², and 4.8552 × 10³ m².

In a case involving the integration of RESs with DSR, there was a reduction in total annual costs to 3.4737 × 10⁶ USD along with a decrease in voltage deviations to 36.0990 p.u. while witnessing an increase in VSI to 2.7245 × 10³ p.u. The best places for the systems in the DN are 26, 31, and 66, respectively. The WTs ratings are 250 kW, 1500 kW, and 6250 kW, respectively; in contrast, the solar modules' comparable areas are 4242 m², 7.7015 × 10⁴ m², and 1.1217 × 10⁴ m².

The study of the most efficient integration of RESs with DSR in the IEEE 118-bus distribution network has resulted in the presentation of voltage profiles for all cases. Figures 10–12 provide a detailed look at the impact of integrating these systems on the overall voltage performance of the DN. The figures allow for a comparison of the voltage profiles before and after integration, providing insights into changes in voltage stability and overall system performance.

<figure><img src="image_16.png" /><figcaption>**Figure 10.** Voltage results for the base case.</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>**Figure 11.** Voltage results after integration of RESs only.</figcaption></figure>
Stability Index, increasing from 2.4440 × 10³ p.u. to 2.6965 × 10³ p.u. The best places for the systems in the DN are 35, 68, and 104, respectively. The WTs ratings are 4250 kW, 6500 kW, and 8500 kW, respectively; in contrast, the solar modules' comparable areas are 5.1878 × 10³ m², 9.9422 × 10³ m², and 4.8552 × 10³ m².

In a case involving the integration of RESs with DSR, there was a reduction in total annual costs to 3.4737 × 10⁶ USD along with a decrease in voltage deviations to 36.0990 p.u. while witnessing an increase in VSI to 2.7245 × 10³ p.u. The best places for the systems in the DN are 26, 31, and 66, respectively. The WTs ratings are 250 kW, 1500 kW, and 6250 kW, respectively; in contrast, the solar modules' comparable areas are 4242 m², 7.7015 × 10⁴ m², and 1.1217 × 10⁴ m².

The study of the most efficient integration of RESs with DSR in the IEEE 118-bus distribution network has resulted in the presentation of voltage profiles for all cases. Figures 10–12 provide a detailed look at the impact of integrating these systems on the overall voltage performance of the DN. The figures allow for a comparison of the voltage profiles before and after integration, providing insights into changes in voltage stability and overall system performance.

<figure><img src="image_20.png" /><figcaption>**Figure 10.** Voltage results for the base case.</figcaption></figure>

<figure><img src="image_21.png" /><figcaption>**Figure 11.** Voltage results after integration of RESs only.</figcaption></figure>
<figure><img src="image_18.png" /><figcaption>**Figure 12.** Voltage results after integration of RESs with DSR.</figcaption></figure>

Figures 10–12 demonstrate that integrating various technologies and system configurations can significantly impact voltage stability and overall distribution network performance. The observed improvements in voltage profiles provide valuable insights into the effectiveness of integrated approaches and can help identify the most suitable case for specific network conditions.

The power produced by the PV units and WT is depicted in Figure 13 in the case of only RESs and Figure 14 in the case of RESs with DSR.

<figure><img src="image_19.png" /><figcaption>**Figure 13.** The hourly power output of RESs for the case of only RESs: (a) output of PV units and (b) output of WT units.</figcaption></figure>

<figure><img src="image_20.png" /><figcaption>**Figure 14.** The hourly power output of RESs for the case of RESs with DSR: (a) output of PV units and (b) output of WT units.</figcaption></figure>
<figure><img src="image_22.png" /><figcaption>Figure 12. Voltage results after integration of RESs with DSR.</figcaption></figure>

Figures 10–12 demonstrate that integrating various technologies and system configurations can significantly impact voltage stability and overall distribution network performance. The observed improvements in voltage profiles provide valuable insights into the effectiveness of integrated approaches and can help identify the most suitable case for specific network conditions.

The power produced by the PV units and WT is depicted in Figure 13 in the case of only RESs and Figure 14 in the case of RESs with DSR.

<figure><img src="image_23.png" /><figcaption>Figure 13. The hourly power output of RESs for the case of only RESs: (a) output of PV units and (b) output of WT units.</figcaption></figure>

<figure><img src="image_24.png" /><figcaption>Figure 14. The hourly power output of RESs for the case of RESs with DSR: (a) output of PV units and (b) output of WT units.</figcaption></figure>
The variation in power output for both WT and PV units is visually demonstrated in Figures 13 and 14. These figures highlight how the power yields of these units experience continuous fluctuations due to variations in wind speed and irradiance.

Figure 15 visually demonstrates the convergence of the IEEE 118-bus system. This figure showcases the convergence curve for cases: (a) after integrating RESs only and (b) after integrating RESs with DSR.

<figure><img src="image_21.png" /><figcaption>**Figure 15.** Convergence curve of the IEEE 118-bus system benchmark cases: (a) after integration of RESs only and (b) after integration of RESs with DSR.</figcaption></figure>

The outcomes presented in Figure 15 demonstrate that the I-WaOA model exhibits a faster convergence rate and attains higher solution quality compared to the WaOA model. Specifically, for the scenarios involving RESs only and RESs with DSR, the I-WaOA algorithm converged to solutions of 0.3616 and 0.1363, respectively. Conversely, the same scenarios led the WaOA algorithm to converge at 0.3971 and 0.1373, respectively. In terms of solution quality, the I-WaOA algorithm yielded lower objective function values for all two scenarios. The objective function value serves as an indicator of how effectively the solution aligns with the optimization criteria, with a lower value signifying a superior solution. It should be highlighted here that the initial iteration process and the results of other algorithms are better. However, the result of the proposed algorithm is the best at the end of the iteration process compared to the other techniques. This is due to the Cauchy mutation, which increases the exploration process of the I-WaOA during the initial iteration process. Furthermore, these algorithms suffer from premature convergence.

## 8. Conclusions

This research offers an Improved Walrus Optimization Algorithm (I-WaOA) to solve the challenging Optimal Operation Problem (OOP) in distribution networks with high penetration of renewable generation. The uncertainties in load demand, wind and solar power generation, and energy prices pose significant challenges for optimal planning and operation. To tackle this, the I-WaOA incorporates the Quasi-Oppositional-Based Learning, the Cauchy mutation, and Fitness-Distance balance schemes to enhance optimization performance. The proposed approach is demonstrated on the IEEE 118-bus distribution network by optimizing three key objectives: minimizing the total system cost, reducing voltage deviations, and improving voltage stability under uncertainties. The optimal allocation and sizing of wind, solar PV, and DSR are simultaneously determined as part of the OOP solution. Monte Carlo simulations and probability density functions are used to model uncertainties in wind speed, load, solar irradiance, solar energy prices, and temperature. Results show that integrating optimally placed and sized renewable generation with DSR leads to significant improvements in all three objectives.
The variation in power output for both WT and PV units is visually demonstrated in Figures 13 and 14. These figures highlight how the power yields of these units experience continuous fluctuations due to variations in wind speed and irradiance.

Figure 15 visually demonstrates the convergence of the IEEE 118-bus system. This figure showcases the convergence curve for cases: (a) after integrating RESs only and (b) after integrating RESs with DSR.

<figure><img src="image_25.png" /><figcaption>**Figure 15.** Convergence curve of the IEEE 118-bus system benchmark cases: (a) after integration of RESs only and (b) after integration of RESs with DSR.</figcaption></figure>

The outcomes presented in Figure 15 demonstrate that the I-WaOA model exhibits a faster convergence rate and attains higher solution quality compared to the WaOA model. Specifically, for the scenarios involving RESs only and RESs with DSR, the I-WaOA algorithm converged to solutions of 0.3616 and 0.1363, respectively. Conversely, the same scenarios led the WaOA algorithm to converge at 0.3971 and 0.1373, respectively. In terms of solution quality, the I-WaOA algorithm yielded lower objective function values for all two scenarios. The objective function value serves as an indicator of how effectively the solution aligns with the optimization criteria, with a lower value signifying a superior solution. It should be highlighted here that the initial iteration process and the results of other algorithms are better. However, the result of the proposed algorithm is the best at the end of the iteration process compared to the other techniques. This is due to the Cauchy mutation, which increases the exploration process of the I-WaOA during the initial iteration process. Furthermore, these algorithms suffer from premature convergence.

## 8. Conclusions

This research offers an Improved Walrus Optimization Algorithm (I-WaOA) to solve the challenging Optimal Operation Problem (OOP) in distribution networks with high penetration of renewable generation. The uncertainties in load demand, wind and solar power generation, and energy prices pose significant challenges for optimal planning and operation. To tackle this, the I-WaOA incorporates the Quasi-Oppositional-Based Learning, the Cauchy mutation, and Fitness-Distance balance schemes to enhance optimization performance. The proposed approach is demonstrated on the IEEE 118-bus distribution network by optimizing three key objectives: minimizing the total system cost, reducing voltage deviations, and improving voltage stability under uncertainties. The optimal allocation and sizing of wind, solar PV, and DSR are simultaneously determined as part of the OOP solution. Monte Carlo simulations and probability density functions are used to model uncertainties in wind speed, load, solar irradiance, solar energy prices, and temperature. Results show that integrating optimally placed and sized renewable generation with DSR leads to significant improvements in all three objectives.
compared to the base case, despite uncertainties. The I-WaOA algorithm outperforms the classical Walrus Optimization Algorithm approach for the complex Optimal Operation Problem. This underscores the advantages of the proposed modifications in handling increased variability and optimizing for high renewable penetration. In conclusion, the paper provides an effective framework for stochastic optimal operation solutions, considering renewable uncertainties. The analysis against the original WaOA method highlights I-WaOA's superior performance, reducing total costs from 3.8377 × 10⁷ USD to 3.4737 × 10⁶ USD, decreasing TVD from 98.6633 p.u. to 36.0990 p.u., and increasing TVSI from 2.444 × 10³ p.u. to 2.7245 × 10³ p.u. compared to the base case. The Improved Walrus Optimization Algorithm can handle the growing complexity of modern power grids with large-scale renewable integration.

**Author Contributions:** Conceptualization, A.T.H.; Methodology, A.T.H., F.S., A.S., M.E., H.E.A.A. and S.A.; Software, A.T.H., M.E. and H.E.A.A.; Validation, A.T.H. and M.E.; Formal analysis, A.T.H.; Investigation, A.T.H., F.S. and A.S.; Writing—original draft, A.T.H.; Writing—review & editing, F.S., A.S., M.E. and S.A.; Visualization, A.T.H.; Supervision, F.S. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was partially supported by the National Research Agency (ANR), LEAP RE "MiDiNA—Microgrids Development in North Africa" project, grant number "ANR-23-LERE-0002-01".

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Data is contained within the article.

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

## References

1. Ammari, C.; Belatrache, D.; Touhami, B.; Makhloufi, S. Sizing, optimization, control and energy management of hybrid renewable energy system—A review. *Energy Built Environ.* **2022**, *3*, 399–411. [CrossRef]

2. Hachemi, A.T.; Sadaoui, F.; Arif, S.; Saim, A.; Ebeed, M.; Kamel, S.; Jurado, F.; Mohamed, E.A. Modified reptile search algorithm for optimal integration of renewable energy sources in distribution networks. *Energy Sci. Eng.* **2023**. [CrossRef]

3. Bakhshinejad, A.; Tavakoli, A.; Moghaddam, M.M. Modeling and simultaneous management of electric vehicle penetration and demand response to improve distribution network performance. *Electr. Eng.* **2021**, *103*, 325–340. [CrossRef]

4. Sadeghian, O.; Nazari-Heris, M.; Abapour, M.; Taheri, S.S.; Zare, K. Improving reliability of distribution networks using plug-in electric vehicles and demand response. *J. Mod. Power Syst. Clean Energy* **2019**, *7*, 1189–1199. [CrossRef]

5. Montoya, O.D.; Giral-Ramírez, D.A.; Hernández, J.C. Efficient integration of pv sources in distribution networks to reduce annual investment and operating costs using the modified arithmetic optimization algorithm. *Electronics* **2022**, *11*, 1680. [CrossRef]

6. Muthukumar, K.; Jayalalitha, S. Optimal placement and sizing of distributed generators and shunt capacitors for power loss minimization in radial distribution networks using hybrid heuristic search optimization technique. *Int. J. Electr. Power Energy Syst.* **2016**, *78*, 299–319. [CrossRef]

7. Kayal, P.; Chanda, C. Optimal mix of solar and wind distributed generations considering performance improvement of electrical distribution network. *Renew. Energy* **2015**, *75*, 173–186. [CrossRef]

8. KN, M.; EA, J. Optimal integration of distributed generation (DG) resources in unbalanced distribution system considering uncertainty modelling. *Int. Trans. Electr. Energy Syst.* **2017**, *27*, e2248.

9. Mansouri, N.; Lashab, A.; Guerrero, J.M.; Cherif, A. Photovoltaic power plants in electrical distribution networks: A review on their impact and solutions. *IET Renew. Power Gener.* **2020**, *14*, 2114–2125. [CrossRef]

10. Martins, J.; Spataru, S.; Sera, D.; Stroe, D.-I.; Lashab, A. Comparative study of ramp-rate control algorithms for PV with energy storage systems. *Energies* **2019**, *12*, 1342. [CrossRef]

11. Mansouri, N.; Lashab, A.; Sera, D.; Guerrero, J.M.; Cherif, A. Large photovoltaic power plants integration: A review of challenges and solutions. *Energies* **2019**, *12*, 3798. [CrossRef]

12. Abd Elazim, S.M.; Ali, E.S. Optimal network restructure via improved whale optimization approach. *Int. J. Commun. Syst.* **2021**, *34*, e4617. [CrossRef]

13. Abd Elazim, S.M.; Ali, E.S. Optimal locations and sizing of capacitors in radial distribution systems using mine blast algorithm. *Electr. Eng.* **2018**, *100*, e4617. [CrossRef]

14. Hachemi, A.; Sadaoui, F.; Arif, S. Optimal Location and Sizing of Capacitor Banks in Distribution Systems Using Grey Wolf Optimization Algorithm. In *Proceedings of the International Conference on Artificial Intelligence in Renewable Energetic Systems, Tamenrasset, Algeria, 20–22 November 2022*; Springer International Publishing: Cham, Switzerland, 2022; pp. 719–728.
compared to the base case, despite uncertainties. The I-WaOA algorithm outperforms the classical Walrus Optimization Algorithm approach for the complex Optimal Operation Problem. This underscores the advantages of the proposed modifications in handling increased variability and optimizing for high renewable penetration. In conclusion, the paper provides an effective framework for stochastic optimal operation solutions, considering renewable uncertainties. The analysis against the original WaOA method highlights I-WaOA's superior performance, reducing total costs from 3.8377 × 10⁷ USD to 3.4737 × 10⁶ USD, decreasing TVD from 98.6633 p.u. to 36.0990 p.u., and increasing TVSI from 2.444 × 10³ p.u. to 2.7245 × 10³ p.u. compared to the base case. The Improved Walrus Optimization Algorithm can handle the growing complexity of modern power grids with large-scale renewable integration.

**Author Contributions:** Conceptualization, A.T.H.; Methodology, A.T.H., F.S., A.S., M.E., H.E.A.A. and S.A.; Software, A.T.H., M.E. and H.E.A.A.; Validation, A.T.H. and M.E.; Formal analysis, A.T.H.; Investigation, A.T.H., F.S. and A.S.; Writing—original draft, A.T.H.; Writing—review & editing, F.S., A.S., M.E. and S.A.; Visualization, A.T.H.; Supervision, F.S. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was partially supported by the National Research Agency (ANR), LEAP RE "MiDiNA—Microgrids Development in North Africa" project, grant number "ANR-23-LERE-0002-01".

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Data is contained within the article.

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

## References

1. Ammari, C.; Belatrache, D.; Touhami, B.; Makhloufi, S. Sizing, optimization, control and energy management of hybrid renewable energy system—A review. *Energy Built Environ.* **2022**, *3*, 399–411. [CrossRef]

2. Hachemi, A.T.; Sadaoui, F.; Arif, S.; Saim, A.; Ebeed, M.; Kamel, S.; Jurado, F.; Mohamed, E.A. Modified reptile search algorithm for optimal integration of renewable energy sources in distribution networks. *Energy Sci. Eng.* **2023**. [CrossRef]

3. Bakhshinejad, A.; Tavakoli, A.; Moghaddam, M.M. Modeling and simultaneous management of electric vehicle penetration and demand response to improve distribution network performance. *Electr. Eng.* **2021**, *103*, 325–340. [CrossRef]

4. Sadeghian, O.; Nazari-Heris, M.; Abapour, M.; Taheri, S.S.; Zare, K. Improving reliability of distribution networks using plug-in electric vehicles and demand response. *J. Mod. Power Syst. Clean Energy* **2019**, *7*, 1189–1199. [CrossRef]

5. Montoya, O.D.; Giral-Ramírez, D.A.; Hernández, J.C. Efficient integration of pv sources in distribution networks to reduce annual investment and operating costs using the modified arithmetic optimization algorithm. *Electronics* **2022**, *11*, 1680. [CrossRef]

6. Muthukumar, K.; Jayalalitha, S. Optimal placement and sizing of distributed generators and shunt capacitors for power loss minimization in radial distribution networks using hybrid heuristic search optimization technique. *Int. J. Electr. Power Energy Syst.* **2016**, *78*, 299–319. [CrossRef]

7. Kayal, P.; Chanda, C. Optimal mix of solar and wind distributed generations considering performance improvement of electrical distribution network. *Renew. Energy* **2015**, *75*, 173–186. [CrossRef]

8. KN, M.; EA, J. Optimal integration of distributed generation (DG) resources in unbalanced distribution system considering uncertainty modelling. *Int. Trans. Electr. Energy Syst.* **2017**, *27*, e2248.

9. Mansouri, N.; Lashab, A.; Guerrero, J.M.; Cherif, A. Photovoltaic power plants in electrical distribution networks: A review on their impact and solutions. *IET Renew. Power Gener.* **2020**, *14*, 2114–2125. [CrossRef]

10. Martins, J.; Spataru, S.; Sera, D.; Stroe, D.-I.; Lashab, A. Comparative study of ramp-rate control algorithms for PV with energy storage systems. *Energies* **2019**, *12*, 1342. [CrossRef]

11. Mansouri, N.; Lashab, A.; Sera, D.; Guerrero, J.M.; Cherif, A. Large photovoltaic power plants integration: A review of challenges and solutions. *Energies* **2019**, *12*, 3798. [CrossRef]

12. Abd Elazim, S.M.; Ali, E.S. Optimal network restructure via improved whale optimization approach. *Int. J. Commun. Syst.* **2021**, *34*, e4617. [CrossRef]

13. Abd Elazim, S.M.; Ali, E.S. Optimal locations and sizing of capacitors in radial distribution systems using mine blast algorithm. *Electr. Eng.* **2018**, *100*, e4617. [CrossRef]

14. Hachemi, A.; Sadaoui, F.; Arif, S. Optimal Location and Sizing of Capacitor Banks in Distribution Systems Using Grey Wolf Optimization Algorithm. In *Proceedings of the International Conference on Artificial Intelligence in Renewable Energetic Systems, Tamenrasset, Algeria, 20–22 November 2022*; Springer International Publishing: Cham, Switzerland, 2022; pp. 719–728.
15. Ali, E.S.; Abd Elazim, S.M.; Abdelaziz, A.Y. Ant lion optimization algorithm for renewable distributed generations. *Energy* **2016**, *116*, 445–458. [CrossRef]

16. Ramadan, A.; Ebeed, M.; Kamel, S.; Ahmed, E.M.; Tostado-Véliz, M. Optimal allocation of renewable DGs using artificial hummingbird algorithm under uncertainty conditions. *Ain Shams Eng. J.* **2023**, *4*, 101872. [CrossRef]

17. Hasanov, M.; Boliev, A.; Suyarov, A.; Urinboy, J.; Jumanov, A. Optimal Integration of Wind Turbine Based Dg Units in Distribution System Considering Uncertainties. In *Khasanov, Mansur*; et al. "Rider Optimization Algorithm for Optimal DG Allocation in Radial Distribution Network.", *Proceedings of the 2020 2nd International Conference on Smart Power & Internet Energy Systems (SPIES)*, 15–18 September 2020, *Bangkok, Thailand*; IEEE: Piscataway, NJ, USA, 2020; pp. 157–159.

18. Elkadeem, M.R.; Elaziz, M.A.; Ullah, Z.; Wang, S.; Sharshir, S.W. Optimal planning of renewable energy-integrated distribution system considering uncertainties. *IEEE Access* **2019**, *7*, 164887–164907. [CrossRef]

19. Khasanov, M.; Kamel, S.; Houssein, E.H.; Rahmann, C.; Hashim, F.A. Optimal allocation strategy of photovoltaic-and wind turbine-based distributed generation units in radial distribution networks considering uncertainty. *Neural Comput. Appl.* **2023**, *35*, 2883–2908. [CrossRef]

20. Mahdavi, M.; Schmitt, K.; Bayne, S.; Chamana, M. An Efficient Model for Optimal Allocation of Renewable Energy Sources in Distribution Networks with Variable Loads. In *Proceedings of the 2023 IEEE Texas Power and Energy Conference (TPEC)*, College Station, TX, USA, 13–14 February 2023; IEEE: Piscataway, NJ, USA, 2023; pp. 1–6.

21. Khodadadi, A.; Abedinzadeh, T.; Alipour, H.; Pouladi, J. Optimal resilient operation of smart distribution network in the presence of renewable energy resources and intelligent parking lots under uncertainties. *Int. J. Electr. Power Energy Syst.* **2023**, *147*, 108814. [CrossRef]

22. Gangwar, T.; Padhy, N.P.; Jena, P. Storage allocation in active distribution networks considering life cycle and uncertainty. *IEEE Trans. Ind. Inform.* **2022**, *19*, 339–350. [CrossRef]

23. Safari, A.; Karimi, M.; Najmi, P.H.; Farrokhifar, M. Multi-objective model for simultaneous distribution networks reconfiguration and allocation of D-STATCOM under uncertainties of RESs. *Int. J. Ambient. Energy* **2022**, *43*, 2577–2586. [CrossRef]

24. Ghaffari, A.; Askarzadeh, A.; Fadaeinedjad, R. Optimal allocation of energy storage systems, wind turbines and photovoltaic systems in distribution network considering flicker mitigation. *Appl. Energy* **2022**, *319*, 119253. [CrossRef]

25. Jahani, M.T.G.; Nazarian, P.; Safari, A.; Haghifam, M. Multi-objective optimization model for optimal reconfiguration of distribution networks with demand response services. *Sustain. Cities Soc.* **2019**, *47*, 101514. [CrossRef]

26. Mohammad Hoseini Mirzaei, S.; Ganji, B.; Abbas Taher, S. Performance improvement of distribution networks using the demand response resources. *IET Gener. Transm. Distrib.* **2019**, *13*, 4171–4179. [CrossRef]

27. Osório, G.J.; Shafie-khah, M.; Lotfi, M.; Ferreira-Silva, B.J.; Catalão, J.P. Demand-side management of smart distribution grids incorporating renewable energy sources. *Energies* **2019**, *12*, 143. [CrossRef]

28. Osman, S.R.; Sedhom, B.E.; Kaddah, S.S. Optimal resilient microgrids formation based on darts game theory approach and emergency demand response program for cyber-physical distribution networks considering natural disasters. *Process Saf. Environ. Prot.* **2023**, *173*, 893–921. [CrossRef]

29. Li, K.; Li, S.; Huang, Z.; Zhang, M.; Xu, Z. Grey Wolf Optimization algorithm based on Cauchy-Gaussian mutation and improved search strategy. *Sci. Rep.* **2022**, *12*, 18961. [CrossRef]

30. Tabak, A.; Duman, S. Levy Flight and Fitness Distance Balance-Based Coyote Optimization Algorithm for Effective Automatic Generation Control of PV-Based Multi-Area Power Systems. *Arab. J. Sci. Eng.* **2022**, *47*, 14757–14788. [CrossRef]

31. Ghafari, R.; Mansouri, N. Improved Harris Hawks Optimizer with chaotic maps and opposition-based learning for task scheduling in cloud environment. *Clust. Comput.* **2023**, *1–49*. [CrossRef]

32. Hasanov, M.; Suyarov, A.; Urinboy, J.; Boliev, A. Optimal Integration of Photovoltaic Based DG Units in Distribution Network Considering Uncertainties. *Int. J. Acad. Appl. Res.* **2021**, *5*, 195–198.

33. Ebeed, M.; Ahmed, D.; Kamel, S.; Jurado, F.; Shaaban, M.F.; Ali, A.; Refai, A. Optimal energy planning of multi-microgrids at stochastic nature of load demand and renewable energy resources using a modified Capuchin Search Algorithm. *Neural Comput. Appl.* **2023**, *35*, 17645–17670. [CrossRef]

34. Asaad, A.; Ali, A.; Mahmoud, K.; Shaaban, M.F.; Lehtonen, M.; Kassem, A.M.; Ebeed, M. Multi-objective optimal planning of EV charging stations and renewable energy resources for smart microgrids. *Energy Sci. Eng.* **2023**, *11*, 1202–1218. [CrossRef]

35. Purlu, M.; Turkay, B.E. Optimal Allocation of Renewable Distributed Generations Using Heuristic Methods to Minimize Annual Energy Losses and Voltage Deviation Index. *IEEE Access* **2022**, *10*, 21455–21474. [CrossRef]

36. Amin, A.; Ebeed, M.; Nasrat, L.; Aly, M.; Ahmed, E.M.; Mohamed, E.A.; Alnuman, H.H.; El Hamed, A.M.A. Techno-Economic Evaluation of Optimal Integration of PV Based DG with DSTATCOM Functionality with Solar Irradiance and Loading Variations. *Mathematics* **2022**, *10*, 2543. [CrossRef]

37. Barnwal, A.K.; Yadav, L.K.; Verma, M.K. A Multi-Objective Approach for Voltage Stability Enhancement and Loss Reduction under PQV and P Buses through Reconfiguration and Distributed Generation Allocation. *IEEE Access* **2022**, *10*, 16609–16623. [CrossRef]

38. Diaf, S.; Diaf, D.; Belhamel, M.; Haddadi, M.; Louche, A. A methodology for optimal sizing of autonomous hybrid PV/wind system. *Energy Policy* **2007**, *35*, 5708–5718. [CrossRef]

39. Eltamaly, A.M.; Mohamed, M.A.; Alolah, A.I. A novel smart grid theory for optimal sizing of hybrid renewable energy systems. *Sol. Energy* **2016**, *124*, 26–38. [CrossRef]

40. Zhao, F.; Li, Z.; Wang, D.; Ma, T. Peer-to-peer energy sharing with demand-side management for fair revenue distribution and stable grid interaction in the photovoltaic community. *J. Clean. Prod.* **2023**, *383*, 135271. [CrossRef]