Modified reptile search algorithm for optimal integration of renewable energy sources in distribution networks

This paper introduces a Modified Reptile Search Algorithm (MRSA) designed to optimize the operation of distribution networks (DNs) considering the growing integration of renewable energy sources (RESs). The integration of RESs‐based Distributed Generation (DG) systems, such as wind turbines (WTs) and photovoltaics (PVs), presents a complex challenge due to its significant impact on DN operations and planning, particularly considering uncertainties related to solar irradiance, temperature, wind speed, consumption, and energy prices. The primary objective is cost reduction, encompassing electricity acquisition, PV and WTs unit costs, and annual energy losses. The proposed MRSA incorporates two strategies: the fitness‐distance balance method and Levy flight motion, enhancing its searching capabilities beyond standard Reptile Search Algorithm and mitigating local optima issues. The uncertainties in load demand, energy prices, and renewable energy generation are represented through probability density functions and simulated using Monte Carlo methods. Evaluation involves typical bentchmark functions and a real 112‐bus Algerian DN, comparing MRSA's efficacy with other optimization techniques. Results indicate that the proposed DN optimization program with WTs and PVs integration reduces annual costs by 21.43%, from 6.2715E + 06 to 4.9270E + 06 USD, reduce voltage deviations by 21.67%, from 77.1022 to 60.4007 USD, and enhance system stability by 2.59%, from 2.3699E + 03 to 2.4314E + 03 USD, compared with the base case.

This approach presents multiple advantages, such as reducing greenhouse gas emissions, minimizing energy losses, improving voltage stability, and reducing the dependence on fossil fuels.
In this context, the optimization of RERs integration in DNs has gained significant attention in the literature, as it represents a crucial step towards achieving a sustainable and dependable energy system.For this aim, the integration of RERs-based DG units at different levels of the electrical networks has been recognized as a promising solution to address both the increasing demand for energy as well as to reduce the environmental impact.Additionally, the integration of such systems has shown multiple technoeconomic advantages reducing the total production cost of electricity, as well as enhancing system performance and voltage profiles. 1owever, the presence of DGs in DNs also introduces significant uncertainty, which greatly increases the intricacy of the optimal operation. 2The main sources of uncertainty in DNs are the load demand, the prices of electricity, and the output power of renewable DGs, which are subject to variations depending on weather conditions (solar radiation, temperature, and wind speed).As a result, determining the optimal placement of DGs in DNs is a challenging and strenuous task. 3The selection of appropriate combinations of PV and WTbased DGs has the potential to increase the efficiency and reliability of the DNs by addressing the issues caused by their variable nature. 4DGs can be strategically located and operated in the network to defer major system upgrades, improve voltage regulation, minimize distribution power losses, relieve heavily loaded feeders, and enhance equipment reliability. 5However, inappropriate DGs location and sizing can lead to various negative impacts on the DN, such as voltage instability, increased power losses, harmonic distortion, and even equipment damage, which can affect the reliability and quality of power supply, potentially leading to economic losses and customer dissatisfaction. 6Therefore, an optimal operation of the system is necessary to make full use of the benefits of DGs while mitigating their adverse effects.
A considerable amount of research has been conducted on optimal DG integration from various perspectives.In one study, Mansouri et al. 7 delved into numerous technological challenges that arise from the widespread penetration of PV into DNs, such as active power reduction, frequency regulation, energy storage, and reactive power injection.Another study 8 explored the technical obstacles that arise due to the extensive and intensive integration of PV into the DNs, such as frequency disturbances and voltage limit violations, as well as the resulting stability issues.The implications of integrating PV on a large scale into the DNs were also considered, along with potential solutions.The searchbased dragonfly algorithm was introduced as a proposed solution for the ideal allocation of DGs, taking into account uncertainties in load demand and DGs.In Prakash and Khatod, 9 various strategies for locating and sizing DGs in DNs were evaluated.In Zellagui et al., 10 two popular methods, the particle swarm optimization (PSO) and the firefly method, were proposed as means of achieving operational, financial, and environmental benefits in DNs.In a recent study, 11 a combined approach utilizing the PSO and hybrid enhanced gray wolf algorithms was proposed for the optimum location for DGs.This approach aimed to minimize active power losses, system costs, and emissions while improving the voltage stability index (VSI) and reducing the voltage deviation index (VDI).A bioinspired algorithm based on Monte Carlo simulation was proposed by Hemeida et al. 12 In Rathore and Patidar, 13 the PSO coefficient was employed to minimize total energy losses by optimizing the sizing and deployment of DGs.In El Sehiemy et al., 14 the Slap swarm optimizer was utilized to enhance the operational, financial, and environmental performance of power plants.The equilibrium optimizer method was introduced by Ahmed et al. 15 A highly effective algorithm was proposed for improving the size and placement of DGs within power networks and addressing the issues with microenergy management.In a study by Thokar et al., 16 a bilayer approach to the placement of energy storage systems and PV in DNs was proposed.In Biswal and Shankar, 17 the Strength Pareto Evolutionary Algorithm 2 was proposed as a solution for the problem of capacitor and DGs placement with load uncertainty.In Hadidian-Moghaddam et al., 18 an innovative ant lion optimizer approach was proposed with the goal of lowering energy costs, reducing losses and voltage deviation (VD), and improving reliability.In Ahmed et al., 19 the PSO approach was proposed, utilizing a probabilistic uncertainty modeling approach, for the siting and sizing of DGs.In Rao et al., 20 Monte Carlo simulation was used for the best placement of DGs within power systems.In this research, Ullah et al. 21reate an energy optimization framework for intelligent microgrids with the aim of reducing operational expenses, minimizing emissions, and enhancing availability.In Ali et al., 22 a demand side management strategy is proposed to optimize operational cost, pollution emission, and load coordination in smart grid using multiobjective wind-driven optimization and modeling consumer behavior the paper uses probability density function (pdf) to forecast wind speed for integrating wind power.Hafeez et al. 23 propose an energy management strategy using price-based demand response programs in an internet-of-things-enabled smart grid to schedule smart home appliances.In Hafeez et al., 24 a modular framework for efficient load scheduling in the smart grid, utilizing a hybrid algorithm and real-time pricing data to reduce costs and peak demands, benefiting both residents and power companies.
In Sakr et al., 25 a revised differential evolution algorithm was introduced for the best positioning of DGs.Additionally, the gorilla troops optimizer, which models gorilla social behavior and movement, has been applied to various engineering problems, including PV model extraction. 26he incorporation of the Levy flight distribution has been found to enhance the performance of optimization algorithms for engineering global optimization problems, leading to improved optimal solutions, as demonstrated in several studies. 27Additionally, balancing fitness and distance has been shown to enhance the algorithm's capabilities for searching, as reported in various studies. 28Overall, these findings highlight the potential of incorporating these techniques into optimization algorithms for engineering design problems.
Table 1 lists a comparison between the presented work and other related references for the Optimal Operation problem (OOP) of DNs with RERs.
This research paper proposes a modified version of the Modified Reptile Search Algorithm (MRSA) for the optimal operation of DNs.The proposed approach is compared against several established optimization techniques, including Sand Cat Swarm Optimization (SCSO), 29 PSO, 30 Dandelion Optimizer (DO), 31 Sine Cosine Algorithm (SCA), 32 improved harmony search (IHS), 33 and the conventional Reptile Search Algorithm (RSA). 28The study considers the integration of PV and WT generators into DNs with the primary objective of cost minimization and the summation of VD and VSI over a 24-h planning horizon.The model's performance is assessed in the presence of uncertainty in load, wind speed, solar irradiation, temperature, and energy prices for purchasing.Weather data, such as irradiance and temperature, are crucial input variables for simulating PV, whereas the output power of WT is influenced by several factors, such as turbine size and wind speed.Therefore, an accurate modeling approach is essential due to the sensitivity of both PV and WT.To this end, the study employs actual data on irradiance, temperature, and wind speed from the DN area, namely, from the real 112-bus Algerian DN, to guarantee the precision and practicality of the findings.
The novelty and significant contributions of the developed algorithm compared with existing works are articulated as follows: 1.This study proposes an MRSA specifically designed to solve the optimal planning problem involving RERs within DNs.

Investigating the synergies between RERs (WTs and
PVs), thereby enhancing grid stability and reducing reliance on expensive grid energy purchases.
PV with energy storage This paper is organized as follows: In Section 2, the optimal planning problem formulation.Section 3 outlines the modeling approach for uncertainty parameters.Sections 4 and 5 explain the mathematical formulation of the RSA and the proposed MRSA, respectively.Section 6 provides the simulation results.Lastly, the conclusions of this paper are depicted in Section 7.

| PROBLEM FORMULATION
This section describes the problem formulation that includes the objective functions and the corresponding constraints of the optimal operation.

| Objective function
In this paper, three objective functions are considered which include the following:

| Cost minimization
The objective function taken into consideration comprises the cost of electricity acquired from the network (C Grid ), the cost of PV units (C PV ), the cost of WT (C Wind ), and the yearly cost of energy loss (C Loss ), and it may be expressed as follows 34 : In which, where U Grid represents the cost of purchasing electricity from the grid, and P Grid h ( ) refers to the power withdrawn from the grid. 35 where U Loss is the cost of energy, and P T_Loss h ( ) refers to the total power losses.
where C PV inst. is the cost of installing the PV, is the PV unit's operating and maintenance costs. 34 where CF is a factor affecting capital recovery, P rated_PV is the rated produced power of the PV. 34 where C Wind inst. is the cost of installing the WT, and & is the wind's operating and maintenance costs. 34 where where P rated_wind is the rated produced power of WT, U WT and U PV represent the WTs and PVs purchased cost ($/kW), and P Wind h ( ) and P PV h ( ) represent the yielded power from WTs and PVs systems.
where β and NP the interest rate and system lifetime of the winds turbine or PV unit. 36he generated power from the PV can be calculated as follows: pv pv pv (12)   The surface area utilized by the set of PV, denoted as A pv in (m 2 ), is multiplied by a constant value representing the conversion efficiency of the PV panels, for the purpose of calculating the power generated by the PV, denoted as P pv (kW), I solar insolation in (kW/m 2 ) and η pv represent the instantaneous efficiency of PV panels.The instantaneous efficiency of PV panels is obtained using the following equation 37 : ) .
pv r t a r r t (13)   The effectiveness of the maximum power point tracking equipment is denoted by η t , and η r represents the reference efficiency of the PV panels.The temperature coefficient of efficiency β, typically ranging from 0.004 to 0.006/°C for silicon cells, is also considered.The ambient temperature T a (°C), PV cell reference temperature T r (°C), and nominal operating cell temperature NOCT (°C) are also taken into account.
The WT's produced output power (P WT ) is computed as follows: the rated power of the used WT is 250 kW while its rated velocity W ( ) r is 15 m/s, the cut-in speed (W i ) is 2.5 m/s, and cut-out speed (W o ) the WT taken as 25 m/s. 38

| Enhancement of voltage level
To improve the network performance, the voltage deviations should be kept closed to 1 p.u. value which can be defined as 36,39 where NB represents the number of buses in the network, and V k represents the voltage of the kth bus.

| Improved system stability
The system stability index statement is as follows 40,41 : where VSI k is the voltage stability index, R km represent the resistance of the transmission lines while the X km is its reactance.P k and Q k define the real and reactive power at bus, respectively.The following three objective functions are taken into consideration simultaneously: where RERs and Base are subscripts refers to with RERs and the base case, respectively.ε 1 , ε 2 , and ε 3 are the weighted factors that were selected to be 0.5, 0.25, and 0.25, respectively. 42

| Limitations of the network
  where V min and V max are the lower and upper voltage limits.P Load and Q Load signified the real and reactive load, respectively; I max n , is the maximum allowable current limit of the line; NT defines the number of lines; PF max and PF min are the maximum and minimum of the WT power factor, respectively.
  where Q S and P S are the reactive and real powers of the main network.

| MODELING THE UNCERTAINTIES
In this section, the considered uncertain parameters are represented as follows:

| The probabilistic model of solar irradiance
The Beta pdf has been employed to model the intermittent nature of the solar irradiance, as follows 43,44 : where σ is the standard deviation while μ is the mean value which has been obtained from the historical data.β and α can be obtained from the following equations 45,46 :

| The probabilistic model of wind speed
For modeling the wind speed uncertainty, the Weibull pdf is utilized which can be described as follows 47,48 : (31)   where W represents the wind speed.c and k are the scale and shape parameters of the Weibull pdf.

| The probabilistic model of load demand
Normal pdf is employed to model the demanding uncertainty that can be mathematically described as follows 46 : where σ l and μ l are the standard deviation and the mean value of the loading, respectively.The probability of the loading is divided into subsegments according to the following equation:

| The probabilistic model of price
One of the most significant random characteristics in power systems is the price of electricity, which is an unreliable parameter that is acquired from the grid.To model the probability distribution of the electricity price, the Normal pdf can be used based on its mean value μ EP and standard deviation σ EP , as shown in Equation ( 33) 49 :

| The probabilistic model of temperature
As the surrounding temperature varies continuously.Thus, the temperature is considered as uncertain parameter and it is assumed that the uncertainty of the temperature is represented using the normal probability distribution for modeling the uncertainty of the temperature as follows: where σ T is the standard deviation while μ T is the mean value of temperature.

| REPTILE SEARCH ALGORITHM
The RSA draws inspiration from the hunting and encircling behaviors of crocodiles in the wild.These behaviors involve crocodiles working together to surround and capture their prey. 28SA uses mathematical modeling to mimic these behaviors and create an optimization process that is both gradient-free and population-based.This means that it can tackle optimization challenges of varying complexities, with or without specific constraints.

| Initialization phase
RSA begins the optimization by generating a set of possible solutions, represented as X in Equation ( 35), through a stochastic approach.The algorithm then identifies the best solution obtained and considers it as an approximation of the optimal solution in each iteration.
The set of candidate solutions, X , used in the RSA is generated randomly through the use of Equation ( 36), x i j , represents the value of the solution at the jth position within the ith candidate solution.The number of candidate solutions generated is denoted as N , while the dimension size of the problem is represented by n.

| Exploration phase
Two exploration search methods are utilized in the RSA.Each component is given a stochastic scaling factor to get more varied solutions and investigate different locations.
The position-updating equations used in the exploration phase are designed to imitate the surrounding behavior of crocodiles and are presented in Equation (37).It is noteworthy that the algorithm employs a simple rule to facilitate this behavior.
The values used in this calculation include Best t ( ) j , which denotes the position of the jth element in the optimal solution obtained until the current point.The integer rand is created at random and ranges from 0 to 1, T , the upper limit of iterations and t, the current iteration number.Equation ( 38) is used to calculate the hunting operator, denoted as (i,), for the jth location in the ith solution.Iterations are conducted with a fixed sensitive parameter, β, set to 0.1, to regulate the exploration accuracy during the encircling phase.A reduction function, denoted as R i j ( , ) , is employed to narrow down the search area and its calculation is determined by Equation (39).The random number r 1 , generated between [1 N ], is used to represent a stochastic location in the ith solution, represented by x r J ( , ) 1 .The value of N represents the number of candidate solutions.The Evolutionary Sense ES t ( ) probability ratio varies randomly between 2 and −2 during each iteration according to Equation (40).
This equation describes some variables used in RSA, including r 2 , which is a random number between [1 N ], and ϵ, which is a small value The equation also involves the use of r 3 , which represents a random value among −1 and 1. P i j ( , ) is another variable used in the algorithm, which represents the proportional difference between of the jth position of the best-obtained solution and the jth position of the current solution.This calculation is performed using Equation (41).
This excerpt explains the variables used in the hunting cooperation phase.The average position of the current solution is denoted by M x ( ) i , which is calculated using Equation (41).UB j ( ) and LB j ( ) represent the upper and lower boundaries of the jth position, respectively.The parameter α is fixed at 0.1 in this study.

| Exploitation phase
The RSA algorithm uses two main search strategies, namely, hunting cooperation and hunting coordination, to investigate the search space and locate an ideal solution.These strategies are modeled in Equation (43)  and are used for exploitation mechanisms.
The jth location in the optimal solution obtained until now is represented by the term Best t ( ) . The hunting operator applied to jth location in the ith solution, denoted as η i j ( , ) , is calculated using Equation (38).The disparity | 4653 percentage of the jth element in the current solution compared with the jth element in the best-obtained solution is denoted by P i j ( , ) , and is calculated using Equation (41).R i j ( , ) is a parameter utilized to shrink the search space, and its value is determined using Equation (39).

| MODIFIED REPTILE SEARCH ALGORITHM
The suggested MRSA is based on two improvement methods.The first modification aims to improve the exploration phase of the RSA using the fitness-distance balance (FDB) approach.1][52][53] The FDB selection method is based on score values of the solution candidates.Then, the locations of the crocodiles are modified based on the fitness values and the distance between the candidates and the best solution.The distance values are calculated as follows: Then, construct the vectors of the fitness and the distance values as follows: Then the distance and the fitness value can be normalized as follows: where min and max are the minimum and maximum in the distance and the fitness vectors.The following formula is used to calculate the FDB score: In which [56][57] F I G U R E 5 Schematic diagram of the 112-bus Algerian distribution network in the base case.


where r 3 and r 4 refer to a random value in the range (0, 1).x t ( ) represent the current location of the crocodiles and the best location, respectively.
represents an operator that measures the intensity of the Levy flight.L F is the Levy flight function which can be measured as follows: In which where u and v denote random values that can be obtained from the normal distribution.β refers to a constant that is 1.5.
F I G U R E 6 Expected load profile.
F I G U R E 7 Expected price profile.
Figure 1 shows the MRSA for the optimal operation solution.The methodology for solving the stochastic optimal operation of a distribution grid is depicted in Figure 2.

| SIMULATION RESULTS
The suggested MRSA is applied to solving the real 112-bus Algerian DN.Initially, The proposed approach's performance was assessed through a comparison with other established optimization methods, such as SCSO, 29 PSO, 30 DO, 31 SCA, 32 IHS, 33 and the conventional RSA. 28All the algorithms, including the proposed MRSA, were programmed in MATLAB software (MATLAB R2019b) and executed on a computer with an Intel i7, 2.5 GHz CPU, and 6 GB RAM.The studied cases are presented as follows:

| Validation of the MRSA technique on standard benchmark functions
Here, the proposed MRSA is utilized for 23 classic functions, including the multimodal, unimodal, and fixed-dimension multimodal functions which have been listed in Appendix A. 58,59 For all cases, the parameters are set according to Table 2 and the obtained results are represented over 30 run times.

| Statistical results analysis
This section depicts the performance of the proposed MRSA compared SCSO, 29 PSO, 30 DO, 31 SCA, 32 IHS, 33 and the conventional RSA. 28Table 3 shows the statistical results, including the mean, the worst, the average, the best, and the Wilcoxon p value between the proposed MRSA and the other algorithms.According to Table 3, the proposed MRSA optimizer is superior in terms of the mean, best, and worst values for the F1-F4, and F17-F23 while the obtained results for the reported algorithms are similar for F17-F23.In addition to that for F15 and F5, the RSA is better than MRSA in the value of the best score but the mean value of the MRSA is better.The p values of the Wilcoxon test compared with MRSA and the other optimizers are included in the 7th column of Table 3.When the p value is less than 5%, it is a significant difference between algorithms. 60On the opposite of that when is more than 5%, no significant difference between optimizers.In addition to that if the results of different algorithms are identical, the p value The lifetime of PV (NP PV ) 63 20 The lifetime of WT (NP WT ) 64 20 The energy loss cost K ( Loss ) 65 0.06 USD/kWh The system constraints

Voltage boundaries
  V 0.95 p. u.
1.05 p. u. | 4657 will be N/A.From the p value, it is clear that the proposed MRSA has significant difference compared with SCSO, PSO, DO, SCA, IHS, and the standard RSA for most of the studied benchmark functions, the population size and maximum iteration for all algorithms are 30 and 300, respectively.

| Convergence curves analysis
The convergence carves of the MRSA and other reported techniques including SCSO, 29 PSO, 30 DO, 31 SCA, 32 IHS, 33 and the conventional RSA 28 can be realized from Figure 3.According to the convergence carves, the proposed MRSA has the best convergence speed for the fixed-dimension functions, the unimodal functions, and F I G U R E 11 Voltage profile for the 112-bus Algerian distribution network without a hybrid system.
the multimodal functions.In addition to that, the suggested MRSA is superior and converged to the optimal solution faster than the conventional RSA due to the proposed modifications which can boost the exploration and exploitation phases of the proposed algorithm.

| Boxplot analysis
Boxplots are ideal for displaying data distributions in quartiles that can be used to show the characteristics of data distribution.Figure 4 depicts the Boxplots of the studied MRSA and the other reported optimization algorithms.It is obvious that the boxplots of the MRSA are narrower compared with the other optimizers.

| Solving the operation problem of Algerian DN using the proposed algorithm
In this section, the proposed MRSA is applied to solve the optimal operation and determine the best location for solar PV unit ratings and the WTs on 112-bus Algerian DN.The system, shown in Figure 5, consists of 112 buses and 111 branches with a cumulative load of 3367.60 kW and 3725.70 kVAR.The meteorological data for the wind speed, temperature, and solar radiation are collected over a period of 3 years from Ehsan and Yang. 61Also, the load demand data are collected for 3 years which given in Kaur et al. 62 To validate the effectiveness of the proposed MRSA algorithm, the obtained results have been with those by other meta-heuristic optimization algorithms, including RSA, SCSO, DE, PSO, SCA, and IHS.To ensure a fair and valid comparison, the maximum number of iterations and populations selected for the proposed algorithms were set at 60 and 25, respectively.The presence of uncertainties, including load variations, wind speed, energy price, temperature, and solar irradiance, were considered in accordance with the methodology outlined in Section 3. Figure 6 illustrates the expected day-ahead load demand, while Figure 7 displays the market price for purchasing energy.In addition, The main objective of this paper is to optimize the total annual cost while improving system performance.
Initially, the base case was considered without any RERs integrated into the DN.The results showed that the total annual purchased energy from the grid was 7.2657E + 04 kWh, with energy losses of 2.4901E + 03 kW.The total purchase energy cost was 6.2170E + 06 USD, and the energy loss cost was 5.4533E + 04 USD, resulting in a total annual cost of 6.2715E + 06 USD.The summation of the VD was 77.1022 p.u. and the VSI was 2.3699E + 03 p.u. Table 5 shows the simulation results which have been obtained by different algorithms.According to Table 5, with the application of the proposed MRSA with the ideal allocation of the RERs, the total annual cost has been from 6.2715E + 06 to 4.9270E + 06 USD and the voltage deviations have been reduced from 77.1022 to 60.4007 p.u., while VSI has been improved from 2.3699E + 03 to 2.4314E + 03 p.u.The system voltage profiles with and without RERs are provided in Figures 11 and 12. Judging from Figures 11 and 12, the voltage profile was enhanced considerably with the inclusion of the PV units and the WTs.Additionally, as seen in Figure 13, the power losses have been significantly decreased.The optimal locations of PV units in the DN are identified as 81, 82, and 102, and the optimal locations of WTs are identified as 5, 95, and 112, respectively as shown in Figure 14.The WTs ratings are 1000, 250, and 250 kW, respectively, while the  corresponding area of the solar modules are 4310, 4000, and 5000 m 2 .The statistical outcomes for the objective function achieved through various optimization techniques are tabulated in Table 5.As depicted in Table 6, the results obtained by the proposed MRSA algorithm outperform other optimization techniques in terms of the mean, best, and worst values.The power output of the PV units is illustrated in Figure 15, indicating that their yields fluctuate continuously with irradiance variations.Similarly, the generated powers of the WT vary with changes in wind speed, as demonstrated in Figure 16.obtained results were compared with other optimization algorithms, including the SCSO, DE, PSO, SCA, IHS, and the standard RSA.The numerical results of optimal inclusion of the RERs using the proposed MSA show that the total cost has been reduced from

F
I G U R E 14 Schematic diagram of the 112-bus Algerian distribution network after integration renewable energy resources.

7 |
CONCLUSIONSThis paper solved the OOP of a real 112-bus Algerian DN with ideal integration of the RERs, including the PV units and the WTs using a new MRSA.The proposed MRSA is based on the FDB method and the Levy flight motion selection strategies.The optimal operation was solved by considering the uncertainties of the price, the load, the wind speed, temperature, and the solar irradiance.The proposed MRSA has been applied and tested on standard benchmark functions and the F I G U R E 15 Hourly output power of PV units in the 112-bus Algerian DN.DN, distribution network; PV, photovoltaic; F I G U R E 16 The hourly output power of WT units for the 112-bus Algerian DN.DN, distribution network; WT, wind turbine.

1
Flowchart of the proposed MRSA for optimal operation.MRSA, Modified Reptile Search Algorithm.Steps procedure for solving the Operation of a Distribution Grid Considering RESs Integration.PV, photovoltaic; RESs, renewable energy sources; WT, wind turbine.Selected parameters of the optimizers.
| 4641 F I G U R E 2 T A B L E 2 Statistical results of different optimizers for the standard functions.
T A B L E 3 Constraints and the cost coefficients.
Energy management results of Algerian distribution network.
Note: Bold values indicate the best obtained solutions.Abbreviations: MRSA, Modified Reptile Search Algorithm; PF, power flow; PV, photovoltaic; RESs, Renewable Energy Sources; RSA, Reptile Search Algorithm; WT, wind turbine.
Voltage profile for the 112-bus Algerian distribution network with a hybrid system.
F I G U R E 13 Power losses without and with RERs for the 112-bus Algerian DN.DN, distribution network; RESs, Renewable Energy Sources.Figures 8-10 present the expected wind speed, irradiance, and temperature, respectively.Table4details the cost coefficients for RERs as well as the operational limitations.
Simulation results of objective function for applied algorithms.
T A B L E 6 RSA, in terms of performance and effectiveness.However, this study has certain limitations, such as the absence of considerations for energy storage systems and electric vehicles (EVs).Future work will expand to encompass energy management in distribution systems, including various energy storage systems, like, batteries, compressed air, and Superconducting Magnetic Energy Storage, as well as incorporating EV charging stations.How to cite this article: Hachemi A, Sadaoui F, Arif S, et al.Modified reptile search algorithm for optimal integration of renewable energy sources in distribution networks.Energy Sci Eng.2023;11:4635-4665.doi:10.1002/ese3.1605 T A B L E A1 Unimodal functions.