
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

71672
10.1038/s41598-024-71672-0
Article
A new multi-objective-stochastic framework for reconfiguration and wind energy resource allocation in distribution network incorporating improved dandelion optimizer and uncertainty
Duan Fude 1
Basem Ali 2
Jasim Dheyaa J. 3
Belhaj Salem 4
Eslami Mahdiyeh mahdiyeh_eslami@yahoo.com

5
Khajehzadeh Mohammad mohammad.khajehzadeh@gmail.com

6
Palani Sivaprakasam shiva.palani08@gmail.com

7
1 School of Intelligent Transportation, Nanjing Vocational College of Information Technology, Nanjing, 210000 Jiangsu China
2 https://ror.org/03ase0085 0000 0004 7642 4328 Faculty of Engineering, Warith Al-Anbiyaa University, Karbala, 56001 Iraq
3 https://ror.org/021817660 grid.472286.d 0000 0004 0417 6775 Department of Petroleum Engineering, Al-Amarah University College, Maysan, Iraq
4 https://ror.org/03j9tzj20 grid.449533.c 0000 0004 1757 2152 Computer Science Department, College of Science, Northern Border University, 73222 Arar, Saudi Arabia
5 grid.466821.f 0000 0004 0494 0892 Department of Electrical Engineering, Kerman Branch, Islamic Azad University, Kerman, Iran
6 https://ror.org/02558wk32 grid.411465.3 0000 0004 0367 0851 Department of Civil Engineering, Anar Branch, Islamic Azad University, Anar, Iran
7 https://ror.org/02psd9228 grid.472240.7 0000 0004 5375 4279 College of Electrical and Mechanical Engineering, Addis Ababa Science and Technology University, 16417 Addis Ababa, Ethiopia
6 9 2024
6 9 2024
2024
14 208572 5 2024
29 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Improving the reliability and power quality of unbalanced distribution networks is crucial for ensuring consistent and reliable electricity supply. In this research, multi-objective optimization of unbalanced distribution networks reconfiguration integrated with wind turbine allocation (MORWTA) is implemented considering uncertainties of networks load, and also wind power incorporating a stochastic framework. The multi-objective function is defined by the minimization of power loss, voltage sag (VS), total harmonic distortion (THD), voltage unbalance (VU), energy not-supplied (ENS), system average interruption frequency index (SAIFI), system average interruption duration index (SAIDI), and momentary average interruption frequency (MAIFI). A new improved dandelion optimizer (IDO) with adaptive inertia weight is recommended to counteract premature convergence to identify decision variables, including the optimal network configuration through opened switches and the best location and size of wind turbines in the networks. The stochastic problem is modeled using the 2m + 1 point estimate method (PEM) combined with K-means clustering, taking into account the mentioned uncertainties. The proposed stochastic methodology is implemented on three modified 33-bus, and unbalanced 25-, and 37-bus distribution networks. The results demonstrated that the MORWTA enhanced all study objectives in comparison to the base networks. The results also demonstrated that the IDO had superior capability to solve the deterministic- and stochastic-MORWTA in comparison to the conventional DO, grey wolf optimizer (GWO), particle swarm optimization (PSO), and arithmetic optimization algorithm (AOA) in terms of achieving greater objective value. Moreover, the results demonstrated that when the stochastic-MORWTA model is considered, the power loss, VS, THD, VU, ENS, SAIFI, SAIDI, and MAIFI are increased by 18.35%, 9.07%, 10.43%, 12.46%, 11.90%, 9.28%, 12.16% and 14.36%, respectively for 25-bus network, and also these objectives are increased by 12.21%, 10.64%, 12.37%, 9.82%, 14.30%, 12.65%, 12.63% and 13.89%, respectively for 37-bus network compared to the deterministic-MORWTA model, which is related to the defined uncertainty patterns.

Subject terms

Electrical and electronic engineering
Energy storage
Fuel cells
Renewable energy
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Motivation and background

The issue that constantly occupies the focus of network administrators is the growing need to enhance network efficiency due to the proliferation of industrial nodes and distribution networks1,2. Electric energy distribution companies incur substantial expenses as a result of consequences such as rising system losses and unwelcome voltage reductions resulting from the growth of electric demands3,4. Reconfiguration of the distribution network and optimal resource allocation for distributed generation (DG) are approaches and tools that can be used to partially achieve the aforementioned objectives while spending the least amount of money5,6. Wind energy sources are one of the renewable sources of DG7. Compared to other renewable energy sources, wind utilization has significantly increased. As a source of electric power production, wind energy sources have been integrated into power systems so that they are implemented in the optimal location and with the largest possible size in the installation network8,9 to maximize their benefits. Obtaining the optimal network configuration while confronting the difficulty of network reconfiguration is necessary to modify the path of load flow and improve network efficiency10. Moreover, the optimal and concurrent use of both techniques in distribution networks can successfully increase the network's efficiency. This necessitates the use of a potent optimization technique with high search and estimation capabilities, as well as optimally describing the decision variables of each approach and establishing a comprehensive objective function while taking into consideration a variety of important objectives11,12. In addition, both the production of wind turbines and the demand requirements of the network are uncertain13. Due to the complexity of load flow computations caused by these distribution network uncertainties, it is necessary to use appropriate probabilistic instruments to calculate load flow and determine the optimal decision variables in light of the objective function established for network configuration and wind resource installation14,15.

Literature and research gaps

Numerous studies have addressed the problems of DG resource allocation and distribution network reconfiguration. Traditional methods, particularly heuristic methods, have demonstrated significant computational efficacy when addressing power optimization problems, despite their limited ability to reach global solutions. In recent years, meta-heuristic techniques have become increasingly popular for optimizing power systems. On the other hand, meta-heuristic algorithms can discover comprehensive solutions. Extremely complex real-world problems, which are typically addressed with numerical optimization techniques, can be addressed with meta-heuristic optimization techniques. In16, a genetic algorithm (GA) and PSO are used to optimize the network reconfiguration of a distribution network to enhance reliability and reduce losses. In17, discrete teaching–learning-based optimization (DTLBO) is employed to reconfigure the network to minimize the losses and enhance the voltage condition. To reduce losses and improve the network's voltage profile,18 constructs a multi-objective reconfiguration strategy for the distribution network. A method for reconfiguring a distribution network using a hybrid honey bee colony (HHBC) to reduce losses and unbalance is presented in19. In20, the reconfiguration problem is addressed using a discrete-learning optimization technique to minimize losses and voltage condition enhancement. An improved GA is used in21, to reconfigure the distribution network while minimizing losses. For network reconfiguration, a branch exchange heuristic and a method are employed in22, where they are combined with other methods such as evolutionary meta-heuristics and cluster evaluation. In23, a stochastic model based on Monte Carlo simulation (MCS) is applied to reconfigure the distribution network to reduce power losses and voltage sags as well as to improve reliability while taking into account load uncertainty. This is done by using an unscented transformation with an iterative mountaineering team-based optimization (IMTBO). A chaotic golden flower algorithm (CGFA) for the reconfiguration problem is presented in24. In25, a method for reconfiguring a distribution network is performed while also taking DG allocation into account using binary PSO (BPSO). Utilizing the MCS, the uncertainties of load fluctuation before network reconfiguration are also taken into account to reduce line loss costs and unsupplied energy. In26, the usage of a gray wolf optimization (GWO) is evaluated to minimize the system's power loss and generate a better voltage profile while meeting the operating limits specified by various operational situations. According to27, following an extreme event, the network is reconfigured using graph theory, and the least-cost option is found by using a minimum spreading tree and a radial topology that connects all grid users. In28, two forms of optimal reconfiguration are applied to reduce losses, improving the voltage profile by integrating the simulated annealing technique and modified PSO, and attempting to find an optimal solution while taking time-dependent alterations into account. The use of sectionalized switches, DG, and optimal reconfiguration of networks in tandem with improved PSO (IPSO) reduces switch costs; increases switch reliability, minimize losses, and improve voltage profiles, according to29. By using the multi-objective enhanced black widow optimization (MOIBWO), the integration of DG installation and network reconfiguration is described in30 to decrease losses, voltage variation, and emission of carbon. In31, an improved golden jackal optimization (IGJO) is applied for optimal planning of WTs and electric vehicle parking lots in a distribution network to achieve the lowest possible annual costs of losses, wind energy, battery degradation cost, and network voltage changes. An improved equilibrium optimization algorithm (IEOA) is used in32 to implement the sitting and sizing of wind turbines (WT) in electrical grid to reduce losses, and voltage condition and reliability improvement. The ideal positions are established by minimizing the placement index using the distribution load flow, whereas the optimum amount of WT is investigated using a beluga whale optimization algorithm (BWOA) in the optimization framework in33. In34, the utilization of PV and WT systems is investigated, with the identification of the best size and position of these resources to reduce losses and enhance voltage profile. A modified analytical energy approach is described in35 employing PSO and GA to locate energy resources in distribution networks while decreasing losses and yearly energy loss. In36, a multi-objective artificial electric field algorithm (MOAEFA) is used to determine the best solutions utilizing a fuzzy decision-making technique for the WT placement in the network to reduce losses and voltage variations. In37, a—chaotic sequence spotted hyena optimizer (SHO) is used to propose a deterministic and probabilistic allocation of WTs in the networks to minimize losses and enhance voltage characteristics and stability. A moth-flame optimization (MFO) technique is used in38 to perform an optimal multi-criteria reconfiguration of distribution networks with wind sources to minimize losses and enhancing the reliability. In39, distribution network WT allocation is assessed while attempting to minimize flicker output. In40, the network reconfiguration is studied to improve the power quality indices such as VS, VU and THD improvement. In41, the allocation of various energy resources is studied to determine the best functioning of a grid-connected microgrid while taking into account the energy resource uncertainties and applying the MCS. According to42, the best way to operate virtual power plants in a reconfigurable network is to use MINLP to solve a potential contingency problem in the system's lines. An artificial ecosystem-based optimization with power and energy loss minimization is provided in43 and improved utilizing opposition-based learning. The reconfiguration of a network that incorporates energy resources for minimizing the losses and enhance the voltage profile is studied in44. A northern goshawk optimization (NGO)-based allocation of WTs in distribution networks is described in45 to determine the WT's optimal size and position in order to minimize casualties.

The items that follow are a list of research voids in the literature:Several researchers utilized reconfiguration to take advantage of the network and improve its performance, according to a literature review. In some instances, they have utilized the optimal distribution of DG resources, particularly wind energy resources, to achieve this objective. In recent years, numerous studies have been published under the title "the method of simultaneous use of reconfiguration and resource allocation in distribution networks." However, simultaneous reconfiguration with wind energy sources to improve network performance has received less attention. In contrast, previous research for the hybrid-simultaneous reconfiguration-WT allocation approach did not adequately investigate the existence of uncertainty in the wind turbine power and network demand.

The remaining references, except Refs.23,26,32,41, did not examine how uncertainty affected the solution of the reconfiguration problem or the best site and size of WTs in the network. The majority of the uncertainty is associated with the use of the conventional Monte Carlo simulation (MCS). Due to the high scenario definition, this method takes a long time, has a high computing cost, and uses a wide number of samples from the probabilistic input space. When reducing the scenarios in the Monte Carlo approach, the clustering-based algorithm may result in inaccuracies when determining the output variables mean and standard deviation. While the number of clusters increases together with the elimination of these mistakes, increasing the computation time46,47.

The choice variables of the optimization problem have also been determined using meta-heuristic techniques, which work better than analytical and heuristic methods in terms of finding the global solution. However, it has been found that because of the complexity of the simultaneous hybrid reconfiguration and wind turbine allocation problem, it is preferable to enhance the efficiency of these algorithms by employing unique techniques to avoid premature convergence, which is presented in quantitative studies in the literature.

Contributions

The paper's contributions are summarized according to the identified research gaps in the literature:This study implements a coordinated solving process for multi-objective framework of distribution network reconfiguration and wind turbine allocation (MORWTA) in an unbalanced distribution network for minimizing the losses and increase quality of the power and reliability. Energy not-supplied (ENS), system average interruption frequency index (SAIFI), system average interruption duration index (SAIDI), and momentary average interruption frequency index (MAIFI) are defined as reliability indices, and voltage sag (VS), voltage unbalance (VU), and THD are power quality indices in the comprehensive multi-objective function with eight important objectives in the operation of the network.

The stochastic multi-objective solving process considering uncertainty is implemented on three modified 33-bus, and unbalanced 25-, and 37-bus distribution networks. With the use of PEM, this method accurately calculates the mean value and standard deviation while dependably assessing its statistical qualities, and it has lower processing costs than the Monte Carlo method48,49.

The problem is optimized and the best choice variables are discovered using an improved dandelion optimizer (IDO). In order to achieve better goals, the IDO's superiority over the traditional DO, GWO, PSO, and arithmetic optimization algorithm (AOA) is also evaluated. The traditional DO algorithm50is used in this research because of its straightforward computational structure and simplicity for comprehension, but also because its performance in the simultaneous-hybrid multi-objective problem has been enhanced in order to prevent getting stuck in Local the best through an adaptive inertia weight51. This is because, based on the no-free launch (NFL) theory, a meta-heuristic algorithm cannot be successful in solving optimization problems.

Structure of the paper

The rest of the study is organized as follows. The comprehensive objective function specified in this research is included in "Problem formulation"' section statement of the problem, together with any relevant restrictions. The stochastic model to model the uncertainties is offered in "Stochastic model" section. "Proposed meta-heuristic optimizer" section describes the suggested optimizer, how to improve it, and how to use it for the problem solving. In "Results and discussion" section, the simulation outcomes of the suggested approach using deterministic and stochastic methods are shown and contrasted. "model validation with previous studies" section presents a comparison with prior research, and "Conclusion" section concludes with a summary of the research's primary conclusions.

Problem Formulation

Multi-objective function

The multi-objective function of the coordinated solving process based on the MORWTA considering eight important objectives is defined as minimizing the objectives of power loss, VS, THD, VU, ENS, SAIFI, SAIDI, and MAIFI using the weighted coefficient technique as follows:1 F=ϖ1×PLossPLoss,max+ϖ2×VsagVsag,max+ϖ3×THDTHDmax+ϖ4×VUNVUNmax+ϖ5×ENSENSmax+ϖ6×SAIFISAIFImax+ϖ7×SAIDISAIDImax+ϖ8×MAIFIMAIFImax

where PLoss,max is maximum value of network losses, Vsag,max demonstrates the maximum of voltage sags, THDmax is the highest value of THD, ENSmax,SAIFImax, SAIDImax, and MAIFImax are maximum value of ENS, SAIFI, SAIDI, MAIFI indices. The values mentioned above are derived from the network's load flow in condition of base network before reconfiguration and WT allocation. Furthermore, it is essential that the absolute values of all weight coefficients (ϖ1,ϖ2,ϖ3,ϖ4,ϖ5,ϖ6,ϖ7andϖ8) collectively sum up to 1. In this study, the problem is simulated for a set of different weight coefficients. The authors have come to the conclusion that for a more accurate compromise and a better solution, the coefficients should be equally distributed among all objectives. Therefore, considering that the absolute value of the set of all coefficients should be equal to 1, the value of the weight coefficient of each target is considered equal to 0.125.

Power losses

Minimizing network line losses is a prevalent objective in network operation. The aim of decreasing active power losses is outlined as follows10,39:2 PLoss=Re∑m=1l∑j=acEjp-EjkCUR∗m

where let l represent the overall count of lines. Let Ejp denotes voltage at bus p, and phase j, and Ejk is voltage at bus k, and phase j, ddditionally, let CUR∗m represent the complex conjugate of the current flowing through the mth line between bus p and bus k. The user provides phases a, b, and c.

Ejpdenotethevoltageatbuspandphasej,andEjkrepresentthevoltageatbuskandphasej. Additionally, let CUR∗(m) represent the complex conjugate of the current flowing through the mth line between bus p and bus k. The user provides phases a, b, and c.

Voltage sag

Voltage sag (VS) is determined by quantifying the remaining voltage present in a bus system during a VS event. This research aims to enhance the VS of the entire network by raising the mean VS across all buses10,39,40,55.3 VSAG=1m∑j=1m1n∑i=1nEij

The objective is to determine the voltage at bus i during fault scenarios at bus j (Eij). where, i is the number of buses and the variable j shows the potential faults number. Additionally, VSAG denotes the average voltage of the bus under VS conditions. The analysis of VS is conducted through the execution of network fault analysis. The fault current resulting from the occurrence of a defect in bus p can be mathematically represented as stated in reference40.4 CURpFaulta,b,c=AdmFa,b,cEp(0)a,b,c(UnM+Impppa,b,cYFa,b,c)

The definition of voltage drop in buses is as follows:5 VBus(F)a,b,c=Impppa,b,cCURBus(F)a,b,c

The variable CURp(Fault) represents the p component of the bus injection current vector, ImpBUSa,b,c,AdmFa,b,c,andAdmFa,b,c are the impedance matrix of the three phases of the network, the admittance matrix considering the fault, the bus voltage wattage in all three phases of the grid, respectively. UnM is the unit matrix in the formulation.

The VS refers to the remaining voltage present in the bus system during a VS sag event. This study exclusively focuses on voltage dip occurrences just under fault circumstances. The computation of the residual voltage at the bus is determined using the following equation40:6 Esaga,b,c=EBus(0)a,b,c-EBus(F)a,b,c

Total harmonic distortion

The calculation of total harmonic distortion (THD) minimization involves determining the overall level of harmonic distortion present in the bus voltage, as described in reference10,40.7 \%VTHD,i=Ed,iErms,i×100,i=1,2,3,...,n

8 Vrms,i=E1,i2+Ed,i2

9 Ed,i=∑h=2mEh,i2

where VOLd,i is the distortion component of the bus voltage. In order to mitigate voltage harmonic distortion, efforts are made to limit the maximum THD value over all bus voltages.

Voltage unbalance

The presence of an unbalanced network can be attributed to load imbalances. The calculation of the network imbalance value in this study involves the measurement of voltage asymmetry at the specified buses10,40.10 VUNB=1n∑i=1n∑j=ac100E-,iE+,i

where E+, i represents the positive-sequence voltage of bus i, and E-,i is the negative-sequence voltage at the bus i.

Energy not-supplied

The interruption of grid lines can result in a partial disruption of network load, leading to a phenomenon known as energy not supplied (ENS) for network consumers. The dependability index of ENS is derived by taking into account the outage rate of the grid lines, the average duration of each outage, and the quantity of lost load, as outlined in references23,38:11 ENS=∑iNbrLi×ϑi×℧i

The rate of annual outage of the transmission line (Li) pertains to the mean duration of each outage experienced on a certain line (ϑi). The variable Nbr represents the overall number of lines, while ℧i denotes the demand at bus i.

SAIFI

The SAIFI9,14,16 represents the annual average frequency of power supply interruptions experienced by a client. The calculation is determined using the following mathematical expression.12 SAIFI=∑i∈ρτiβ

Let denote the count of customers impacted by each incidence as τi, the overall count of customers as β, and the set of all outage occurrences in the distribution network as ρ.

SAIDI

The SAIDI9,14 quantifies the mean duration of energy supply disconnection experienced by customers for a year. The calculation is determined using the subsequent expression:13 SAIDI=∑i∈ρ(riτi)β

where ri represents the restoration time for each incident i.

MAIFI

In the assessment of system reliability, the MAIFI9,14 takes into account the occurrence of temporary interruptions. The metric is denoted as interruptions per client per year. The MAIFI is a metric that quantifies the total number of short-term disruptions experienced by consumers, divided by the overall count of customers.14 MAIFI=∑i∈ρτpout,iτsh,iβ

where τpout,i denotes the quantity of short-term power outages, while τsh,i indicates the number of users affected by each individual short-term power outage within the given statistical timeframe.

Constraints

The optimizing the e multi-objective function has to be performed while considering the limitations outlined below. The limitations or restrictions of the problem are outlined in the following manner10,38–40:

Power balance

15 PPost+∑i=1WTPWT(i)=∑i=1brPloss(i)+∑q=1busPd(q)

16 QPost+∑i=1WTQWT(i)=∑i=1brQloss(i)+∑q=1busQd(q)

where PPost, PWT, Ploss and Pd indicate the active power supplied by the post, the active power provided to the network via the wind turbine, active losses of the line, and active power demand of the network demand, respectively. The variables QPost, QWT, Qloss and Qd denote the reactive power obtained from the post, the reactive power delivered into the grid by the wind turbine, the reactive power losses in the transmission line, and the reactive power requirement, respectively. The variable NWT denotes the number of hybrid systems, which, in the context of this study, is equal to 1.

Bus voltage

It is imperative that the voltage levels of the buses are maintained within a predetermined threshold.17 Emin≤Epi≤Emax

Lines current

It is imperative that the current flowing to the network does not beyond the maximum permissible current of the grid lines.18 CURpl≤CURplmax

Voltage unbalance

The electrical potential difference the set value for bus imbalance should not be exceeded.19 E+,iE-,i≤Vunbmax,i=1,2,3,⋯,n

Voltage sag

The extent of VS in the network should not exceed the allowable threshold as defined by the distribution network standards.20 Vsag,av≥Vsagmin

Radiality

To guarantee the radial topology of the distribution network, no loops should be formed, and also every bus of the system should be connected to a substation point which is described fully in Ref.21.

WT power

21 PWTMin≤PWT≤PWTMax

where PWGMin and PWGMax represents the lower and upper limit for wind turbine power generation.

For each of the constraints (15) to (21), a penalty, i.e. a very large number, is embedded in the objective function. During the simulation, if any of the constraints violates the allowed value or the logic governing it, a penalty is added to the objective function and the solution becomes impractical, and infeasible. In this way, unacceptable solutions are removed.

Power flow analysis

Forward–backward power flow is employed for research on distribution networks. Every line is represented using series impedance. A forward–backward sweeping method (FBSM) that iterates performs the power flow in order to obtain convergence and determine line and bus voltages10,12,40. The FBSM consists of two fundamental phases that are simple to comprehend and apply across all radial networks. These two stages function as a ring. In the FBSM, the bus voltage is set to 1 p.u, and the current for line, and also voltage of each buses are simply determined. The voltage is adjusted through the feeder's final position. As previously indicated, the voltage has been set to 1 p.u to acquire the bus flow. This current generates a fresh voltage at the end of the feeder. The processes for power flow are outlined below: The voltage three-phase across all phases is regarded equivalent to the slack bus voltage. Each phase’s load is adjusted by a positive sequence and a 120-degree phase variance. In this stage, the forward–backward power flow is completed. Begin at the ending bus and work your way to the sending bus, updating the upstream current of lines to ensure that the upstream current of line corresponds to the all downstream lines current the fact that start at the ending bus of upstream branch, and this stage gets repeated until the current of branch is identified. Next, the forward load flow is carried out. The calculation moves from the sending bus to the ending bus of the network to determine the voltage at the ending bus of each line, denoted as Vdowna,b,c=Vupa,b,c-Zup,downa,b,cIup,downa,b,c. Here, Vdowna,b,c represents the voltage of each phase at the load side (referred to as the lower bus), Vupa,b,c refers to the phase voltage of the sending bus in the line between the up and down buses, Zup,downa,b,c denotes the branch impedance, and Iup,downa,b,c signifies the current in the line between the up and down buses.

The forward load flow process continues until the voltages of all buses are determined. The forward–backward load flow steps are repeated until the system reaches convergence. This method calculates the voltages at all buses at the main frequency. The harmonic flow problem is addressed using direct harmonic analysis. Harmonic-generating loads are represented by their Norton equivalents, while the lines and linear loads are modeled by their respective admittances at harmonic frequencies.

To maintain the radial structure of the network, the branch exchange method is applied. In this method, when a line is closed, another line must be opened to prevent the formation of loops and ensure the network remains radial. The advantage of this approach is that it automatically ensures the network's radiality without the need for additional computations or algorithms to detect the network's configuration.

Stochastic model

Various schemes can be implemented within the PEM approach48, each characterized by a different number of deterministic methodologies to solve. Specifically, the schemes discussed in this paper require 2m + 1 deterministic runs, where m represents the number of random input variables. Additionally, 2m or 4m + 1 runs can be used. The more deterministic methodologies run, the more accurate the estimation of the statistical moments of the output random variables. However, an increase in the number of runs necessitates more information about the input random variables. For instance, in the 4m + 1 scheme, it is necessary to know up to the eighth standard central moment of the input random variables, resulting in a higher computational burden. Nevertheless, the PEM approach offers a significant reduction in computational efforts compared to the MCS method. Therefore, this paper uses the 2m + 1 PEM method to model uncertainties with fewer scenario definitions, striking a balance between the accuracy of estimating statistical moments of the output random variables and computational cost49.

The stochastic solving process for multi-objective optimization of network reconfiguration and WT sitting and sizing in a distribution network is conducted considering uncertainties in WT generation and network load demand. This study utilizes a 2m + 1 PEM in conjunction with the IDO method to address the challenges of distribution network layout and wind turbine allocation within an uncertain environment. The primary objective is to determine the outputrandomvariableORV uncertainty level by utilizing the inputrandomvariableIRV central moments at k specific locations for each variable. These points are commonly referred to as concentrations and are assessed using various mathematical functions to establish a relationship between IRVs and ORVs. This approach offers the benefit of precisely calculating the mean and std through the use of the PEM while also reliably assessing statistical features. The 2m + 1 approach often employs three concentration sites for each random variable, with one of these sites positioned at its mean48,49.

Adhering to the principles of the PEM methodology, the ORV as Y is evaluated 3 times for each IRV, while the mean value is used for the remaining input random variables. Yi,m denotes the array of ORVs linked to the mth point of concentration related to the ith IRV, its computation involves the following procedure22 Yi,m=F(μ1,⋯,pi,m,⋯,μn)

where μ1,⋯μn correspond to the means of the IRVs, excluding the ith input variable. The concentration point corresponding to the mth instance of the ith IRV is denoted by the pair (pi,m,wi,m), with pi,m indicating the location and wi,m representing the concentration point weight.

The calculation of position pi,m can be expressed in relation to the mean and std of the ith IRV, denoted as μi and σi, in the following manner48,49:23 pi,m=μi+ξi,mσi

where ξi,m represents the location of the mth quantity of the ith IRV.

The determination of the position (ξi,m) and weights (wi,m) for the corresponding IRV is presented by48,4924 ξi,m=0.5λi,3+(-1)3-m×λi,4-0.75λi,32m=1,2,ξi,3=0;

25 λi,3=Epi-μi3σi-3

26 ωi,m=-13-mξi,mξi,1-ξi,2-1,m=1,2,

27 ωi,3=n-1-λi,4-λi,32-1

In this instance, λi,jcorrespondstothejthcentralmomentofstandardrelatedtotheithIRV. Subsequently, a deterministic power flow computation is executed using the mean of point related to the concentration as {μ1,⋯,μn}, resulting in the determination of Yμ. Consequently, the qth moment of ORVs, denoted as μyq can be acquired using the following equation48,49:28 μyq=∑i=1n∑m=13ωi,mYi,mq+ωμYμq

where ωμ=1-∑i=1n∑m=13ωi,m denotes the weighted coefficients for average point of concentration.

The uncertain input data set, consisting of network load and WT power, is grouped using the K-means clustering approach. A set of k starting agents is selected at random to serve as input sample points. Every individual data point is allocated to the cluster that corresponds to the nearest centroid. Subsequently, the means of each cluster are computed based on the individuals inside each cluster, and the agent associated with each cluster is subsequently updated. The latter two steps must be iteratively performed until the alterations in the clusters are below a specified threshold. Once the agents {a1,⋯,ak} have been identified, the approach based on clustering in the MCS is applied to each cluster agent, resulting in the acquisition of the output random variable Y′. Hence, the utilization of cluster agents is limited to serving as a condensed representation of the complete data points in the input space for MCS computations.

The outcomes derived from the clustering technique are subsequently integrated with the average and std acquired from the PEM:29 Yj=μy+Yj′-μy′σy′σyj=1,…,k

By utilizing the final findings Y, along with the average and std of the ORV, it becomes possible to compute the statistical characteristics of the ORV including distribution functions, probabilistically. The standard location in 2m + 1 schemes is independent of the IRVs total number. In this framework, it is evident that setting λi,3=0 restricts the variable to be centered on its mean, randomly. Consequently, all potential estimations of 2m + 1 schemes exhibit three deterministic points, accompanied by the distinctive mean value of the corresponding IRV. Flowchart of the 2m + 1 PEM method is depicted in Fig. 1.Fig. 1 Flowchart of 2m + 1 PEM method.

Proposed meta-heuristic optimizer

This section provides a description of the framework of the suggested hybrid optimization approach for problem resolution, followed by an explanation of its implementation in problem solving.

Overview of DO

Dandelions undergo population evolution through a three-stage process wherein they transfer their seeds to subsequent generations. The primary source of inspiration for the concept of Dandelion optimizer (DO) is derived from the aforementioned three stages50.

Initialization

The DO follows the principles of iterative optimization and population evolution through population initialization. Within The recommended DO takes each dandelion seed into account a potential candidate solution, and the entire population is represented by5030 Population=x11⋯x1Dim⋮⋱⋮xpop1⋯xpopDim

where "pop" represents the population size, and "Dim" signifies the variable dimension. Each potential solution is randomly generated within the constraints of the upper limit (UB) and lower limit (LB) specific to the problem at hand. The expression for the ith individual, denoted as Xi, is formulated as follows50:31 Xi=rand×UB-LB+LB

where i is an integer ranging from 1 to pop, and rand represents a randomly generated number within the interval of 0 to 1. The LB and UB are defined as follows:32 LB=lb1,…,lbDimUB=ub1,…,ubDim

During the initialization process, the DO selects the individual corresponding thighest fitness quantity in of the primary elite candidate. This candidate is perceived as the maximum favorable dandelion position grows to thrive. Using the lower limit as an illustration, the representation of the primary elite candidate, denoted as Xelite, is50:33 fbest=minfXiXelite=Xfindfbest==fXi

where find() indicates two indices with same quantities.

Stage of rising

During the initial phase of development, dandelion seeds must attain a specific altitude in order to disperse from their parental plant. Dandelion seeds ascend to varying altitudes due to factors such as wind velocity and atmospheric moisture content. In this location, the weather can be categorized into two distinct cases50.

Case 1: On a clear day, the wind velocities can be represented by a lognormal distribution lnY ∼ N(μ,σ2)50. Within this distribution, values are more densely spread randomly during the Y-axis, leading to an increased likelihood of dandelion grows being carried to distant areas. As a result, in this situation, the DO approach places a greater emphasis on exploration. In the search domain, the wind disperses dandelion seeds in a random manner, causing them to land in various locations. The elevation to which a dandelion seed rises is influenced by the wind's velocity. Stronger winds propel dandelions higher, causing their seeds to scatter over greater distances. Due to the variable nature of wind speed, the vortex formations above dandelion seeds undergo continuous adjustments, inducing a spiral ascent. This phenomenon can be mathematically expressed as follows48:34 Xt+1=Xt+α* vx*vy∗lnY∗Xs-Xt

where Xt is position of the dandelion kernel during iteration t, whereas Xs is position randomly chosen from the search space throughout iteration t. Equation (35) furnishes the formula for generating this randomly selected position.35 Xs=rand1,Dim∗UB-LB+LB

ln Y represents a lognormal distribution with parameters μ=0 and σ2=1, and its mathematical formulation is as follows50:36 lnY=1y2πexp-12σ2lny2y≥00y<0

where, y stands for an independent variable conforming to the normal distribution N (0,1). The parameter α is adaptively utilized to modify the search step length, and its mathematical representation is50:37 α=rand()∗1T2t2-2Tt+1

where variable α represents a random perturbation ranging from 0 to 1. The aforementioned oscillations prompt the algorithm to prioritize global search during the initial phase and transition to local search during the subsequent phase. The variables vx and vy are used to clear the lift component coefficients associated with the separated eddy action on a dandelion. The magnitude of the force applied on the variable dimension is calculated by5038 r=1eθ

39 vx=r∗cosθ

40 vy=r∗sinθ

Let θ is a random variable that follows a uniform distribution in the interval [− π,π].

Case 2: During inclement weather conditions, such as a wet day, the ability of dandelion seeds to effectively disperse through the wind is hindered due to the presence of air resistance, humidity, and various other influencing factors. In this particular scenario, the utilization of dandelion seeds occurs within the immediate vicinity, and the associated mathematical representation is5041 Xt+1=Xt∗k

where the parameter k is employed to control the extent of the local search domain for a dandelion, and Eq. (42) is utilized to compute the domain50.42 q=1T2-2T+1t2-2T2-2T+1t+1+1T2-2T+1

43 k=1-rand()∗q

The variable k demonstrates a pattern of "downwards convex" oscillation, which facilitates the algorithm's ability to exploit local conditions by employing a large stride during the initial phase and a smaller stride during the subsequent phase. At the end of each cycle, the parameter k demonstrates a gradual approach towards the value of 1.

The formulation of dandelion grows during the rising step can be presented as:44 Xt+1=Xt+α∗vx∗vy∗lnY∗Xs-Xtrandn<1.5Xt∗kelse

where randn () represents a number following the standard normal distribution, random.

Figure 2, illustrates the flight patterns of dandelion seeds in varying meteorological circumstances. The figure provides an approximation of the regeneration locations for dandelion seeds. In clear weather conditions, dandelion seeds undergo updates that prioritize the exploration process by incorporating randomly picked location information. The eddy located above the seed exerts an influence on the vector in motion by performing a multiplication operation on its x and y components, thereby adjusting the trajectory of the spiral movement pattern of the dandelion. In the second scenario, the dispersal of dandelion seeds occurs in many directions within the immediate vicinity of the community. The utilization of standard distribution of arbitrary numbers is employed to dynamically regulate the processes of utilization and discovery. To improve the method's focus on search globally, a cut-off point of 1.5 is established. In the initial step of iterative optimization, the chosen setting facilitates the extensive traversal of dandelion seeds throughout the whole search space, with the aim of obtaining accurate directional guidance for subsequent stages.Fig. 2 The schematic diagram depicting the dandelion seeds rising stage50.

Descending stage

During this stage, the proposed DO continues to prioritize exploration. Following their ascent to a specific height, dandelion seeds experience a gradual descent. The algorithm employs Brownian motion to simulate the movement trajectory of these seeds. This approach enables individuals to discover a broader range of search spaces during iterative updates, given that Brownian motion adheres to a standard distribution with each alteration. In order to assess the dandelion descent stability, the mean location data obtained during the rising phase is employed. This strategy contributes to steering the entire population toward favorable areas for further development. The corresponding mathematical representation is as follows50:45 Xt+1=Xt-α∗βt∗Xmean_t-α∗βt*Xt

where βt represents Brownian movement and is a number drawn from the standard distribution, randomly. Xmean-t represents the average position of the population at the ith iteration, and the formulation is given by50:46 Xmean_t=1pop∑i=1popXi

The repeated updating of individuals relies on the average positional data of the population. This information plays a key role in determining the evolution direction of individuals.

Landing stage

During this phase, the DO shifts its focus towards exploitation. Building on the insights from the preceding stages, the dandelion seed undergoes a random selection process to determine its landing site. As iterations unfold, the algorithm aims to gradually converge towards the optimum solution, globally. Consequently, the identified optimum solution represents the estimated position where dandelion spores are likely to be found to thrive. To facilitate precise convergence to the global optimum, search agents utilize the valuable information from the current elite to exploit their local neighborhoods. As the population evolves, the ultimate revelation of the optimum solution globally becomes achievable. This behavior is mathematically captured by5047 Xmean_t=Xelite+levyλ∗α∗(Xelite-Xt∗δ)

In this equation, Xelite denotes the optimum location of the dandelion seed in the ith iteration. The function Levy(λ) encapsulates the Levy flight and can be determined by5048 levyλ=s×ω×σt1β

where β is a randomly generated number within the range [0, 2] (β = 1.5 in this study). The constant value of s is set to 0.01. Moreover, w and t represent random numbers chosen from the interval [0, 1]. The mathematical formulation for α is provided by5049 σ=Γ(1+β)×sin(πβ2)Γ(1+β2)×β×2(β-12)

where β remains constant at 1.5. Furthermore, δ is determined as a linearly increasing function within the range [0, 2], and its computation is defined as follows:50 δ=2tT

In order to effectively approach the solution globally, the utilization of a linear rising function is employed on individuals as a means to prevent excessive utilization. The Levy flight coefficient is employed at this stage to simulate the magnitude of individual movement steps.

Execution of the DO

This subsection presents a thorough explanation of the precise sequence of steps involved in the execution of the DO algorithm. In the case of DO, vectors clusters produced randomly within searching for space to begin optimization. Subsequently, dandelion grows undergo three distinct stages optimization by iteration cycles: rising, descending, and landing. It is important to highlight the protocols for generating new populations prior to the subsequent iteration. To elucidate this approach, consider the scenario where each dandelion grows is ranked in increasing order based on their fitness values. The individual corresponding to minimum fitness assumes the role of the elite individual for the subsequent production of the individuals. Consequently, the population is sorted, serving as the foundational population for the ensuing iteration. This arranging technique proves advantageous for preserving valuable information and prevents the algorithm from inadvertently discarding accurate optimization directions due to the iterative progression. The algorithm concludes its optimization process by reaching the predefined maximum iteration count.

Flowchart of the DO is depicted in Fig. 3.Fig. 3 Flowchart of the DO.

The pseudo-code of the DO is presented in Algorithm 1.

Overview of improved DO (IDO)

The adaptive inertia weight (∅)51, adapted automatically based on the present physiological condition of the agent, is added to the local search and exploitation phase to prevent the algorithm from getting trapped in the local optimum. In this study, ∅ is used in the exploitation phase to increase diversity and avoid entering the local optimal. The position of the clouded leopard in the second phase, i.e., the exploitation phase, is updated using ∅ as follows:51 Xmean_t=∅*Xelite+levyλ∗α∗Xelite-Xt∗δ

The converging modification is governed by the operative ([0, 1]), which is automatically adjustable based on the present operator location. While the efficiency of a factor exceeds the mean quantity of the entire population, the weight operator ∅ of is assigned to its utmost quantity so that the population can quickly converges on the optimal location. Alternatively, as the current agent's fitness is below average, ∅ is designated a lesser quantity. The less the weight quantity, the weaker the impact of revising the dandelions' status, deteriorating the performance will be. The approach ∅ enhances the algorithm's localized search capabilities and prevents abrupt convergence. The accompanying equation explains depending on the improved Versoria mapping function.52 ∅mt=1-1φ.ait-122+2ifait≤0.51φ.ait-122+2otherwise

where ai(t),ai(t) ∈ [0,1] indicates the dandelion's classification among all other dandelions within the lowest fitness and the mean fitness, as specified below:53 favet=∑m=1Npfi(t)Np

54 fmint=minf1t,f2t,⋯,fNpt

55 ait=fit-fmintfavet-fmint

where fi(t) is the fitness of the ith dandelion at iteration t, fave(t) and fmin(t) represent the mean and smallest fitness values of all Np size populations, respectively.

The IDO steps to solve the problem

The steps of the IDO implementation for addressing the deterministic MORWTA problem solving. Then, to solve the stochastic problem, the steps described in "Stochastic model" section have been solved considering the deterministic model. Flowchart of the IDO implementation to solve the deterministic approach is depicted in Fig. 4.Fig. 4 Flowchart of the IDO implementation to solve the deterministic approach.

Step 1: Apply network data, which includes load demand data, network lines data, tie-lines, and WT maximum capacity.

Step 2: Identify the upper bound and lower boundary of each tie-line, and additionally the lower and upper size, power factor, and installation location of the WT (decision variable set) for each Dandelion of the algorithm's population. To solve simultaneous reconfiguration and WT allocation, the first component of the decision variable vector is the number of open distribution network branches, while the second, third, and fourth components are the WT installation location (bus), WT size, and power factor (p.f) of the WT, respectively.

Consequently, the solution vector (X) for the simultaneous-hybrid MORWTA problem corresponds to follows:56 X=DVWT¯DVTLs¯

57 DVWT=ILWTSWTp.fWT

58 DVTLs=TL1TL2TL3TL4TL5

where DVWT¯ and DVTLs¯ are decision vectors of WT and network tie-lines, respectively. ILWT, SWT and p.fWT are site, capacity and power factor of the WT, and TL1, TL2, TL3, TL4, and TL5 are the five network tie-lines.

Step 3: On the basis of the initial population of the Dandelion, the backward-forward power flow meeting the network radiality is utilized for calculating the objective function (Eq. (1)) of each Dandelion satisfying the constraints of Eqs. (15)–(21). The Dandelion with the best fitness is saved to the best member of the algorithm according to the fitness values.

Step 4: Update the population based on the DO in the fourth step. The solution vector has been chosen at random for each member of the updated population, and the fitness function (Eq. (1)) is computed for each set of variables by satisfying the operational and radial constraint (Eqs. (15)–(21)). The finest population member is the Dandelion with the highest value among the updated Dandelions. If its fitness function is lesser, the current best member of the population is substituted with the best acquired in step 3.

Step 5: Utilizing adaptive inertia weight, update the algorithm's population.

Step 6: The solution vector is selected randomly for each updated population according to the adaptive inertia weight (Eqs. (51)–(55)), and the fitness function (Eq. (1)) is calculated for each variable set satisfying the operational and radial constraints (Eqs. (15)–(21)). The best population member is the Dandelion with the highest value among the updated Dandelions. If its fitness function is lesser, the current best member of the population is replaced with the best gained in step 4.

Step 7: Convergence criteria. By executing the utmost number of iterations, repetition and achieving the best objective function value convergence criteria of the IDO are examined. If yes, go to step 4; otherwise, go to step 8.

Step 8: Terminate the algorithm and save the results.

Results and discussion

The results of many-objective reconfiguration and wind turbine allocation (MORWTA) for three modified 33-bus, unbalanced 25-bus, and 37-bus distribution networks under two different scenarios as Scenario#1) deterministic-MORWTA (without uncertainty), and Scenario#2) stochastic-MORWTA (with uncertainty) are presented. In this research, the line outage rate is considered to be 0.046 per km/yr54. It is also incorporated that all other components are 100% availability. In the event of a branch outage, the time within the outage and repair is considered to be 8 h, and interruption duration for downstream customers is 8 h54. Moreover, the MORWTA based on the IDO is compared with conventional DO, GWO, PSO, and AOA to demonstrate its superiority. Each algorithm's general and regulatory parameters are regarded as being identical to their respective reference papers. The number of population, the maximum number of iterations, and the independent runs number have been set to 50, 100, and 30 (according to the trial and error and user experience) for a fair comparison of the effectiveness of various optimization methods. In this study, besides the rearrangement of the distribution network, the allocation of a 1 MW turbine is also considered. In the stochastic framework, considering the uncertainty, the network loading probability distribution function (PDF) as well as the wind turbine power is considered as Fig. 5, based on which the output random variables are considered for the 2m + 1 PEM method.Fig. 5 PDF of (a) distribution network loading, and (b) wind power considering 1000 scenarios.

The control parameters of each algorithm are considered as in their reference article for fair competition according to Table 1. Table 1 The control parameters of different algorithms.

Algorithm	Parameter	Value	
Dandelion optimizer (DO)50	K	[0 1]	
α	[0 1]	
Particle swarm optimization (PSO)16,19	W1, W2	0.8, 0.4	
C1, C2	2, 2	
Grey wolf optimizer (GWO)52	a→	[2, 0]	
Arithmetic optimization Algorithm (AOA)53	α	5	
μ	0.5	

The simulation results are presented in two below scenarios:

Scenario#1: Deterministic-MORWTA without uncertainty.

Scenario#2: Stochastic-MORWTA with uncertainty.

Table 1 The control parameters of different algorithms.

Results of deterministic, and stochastic-MORWTA for modified 33-bus network

The schematic of the 33-bus network is depicted in Fig. 6. There are 33 buses, and 37 lines with five tie-lines within buses 8–21, 9–15, 12–22, 18–33, and 25–29. The 33-bus network has 3715 kW of active load and 2300 kW of the reactive load. In this simulation, the base voltage and the base apparent power are considered to be 12.66 kV and 100 MVA, respectively. The load data of this network is modified and an unbalanced network is extracted. The load and lines data of this modified network are derived from Ref.55.Fig. 6 Schematic of the 33-bus distribution network.

Deterministic results for modified 33-bus network (Scenario#1)

In the base network, the objective values for power loss, VS, total harmonic distortion, VU, ENS, SAIFI, SAIDI, and MAIFI are calculated 207.44 kW, 0.96 p.u, 0.1835%, 4.92%, 6.44 MWh, 1.36 failure/customer/yr, 1.21 h/customer/yr and 4.64 mom/customer/yr, respectively.

In this section, the simulation outcomes of the first deterministic-MORWTA scenario utilizing the IDO are presented and compared to the conventional DO, GWO, PSO, and AOA methods. Figure 7, depicts the convergence process derived from various optimization techniques. The IDO algorithm with the lowest convergence tolerance has obtained the greatest solution value with the least objective function value in comparison to other algorithms. Also, Fig. 7 demonstrated that improving the conventional DO has increased the system's effectiveness in locating the optimal solution with a quicker convergence rate.Fig. 7 Convergence process to solve scenario#1 for modified 33-bus distribution network.

Table 2, presents the numerical results of the first scenario based on deterministic-MORWTA derived from various algorithms. IDO is preferable to other algorithms in solving the deterministic-MORWTA problem because the outcomes of the objective function and the eight problem objectives achieved through IDO are smaller than those derived by other algorithms. The IDO has put in 894 kW of WT generation with a p.f of 0.86 on bus 7 and determined the network configuration of 33-bus by considering tie lines as lines 10, 20, 28, 30, 35. As a consequence, the value of each of lines 10, 20, 28, 30 and 35. The values of power loss, VS, THD, VU, ENS, SAIFI, SAIDI, and MAIFI are obtained 112.79 kW, 0.42 p.u, 0.1725%, 2.07%, 5.13 MWh, 1.16 failure/customer/yr, 1.03 h/customer/yr, and 2.38 mom/customer/yr, respectively. In addition, the results of the statistical analysis presented in Table 3 demonstrate that the IDO-solved deterministic-MORWTA scenario has yielded more effective results than those of other algorithms by attaining superior indicators of statistical effectiveness. Table 2 Numerical results of scenario#1 optimization for modified 33-bus distribution network.

Method/Item	Base Net	IDO	DO	GWO	PSO	AOA	
Opened switches	33, 34, 35, 36, 37	10, 20, 28, 30, 35	7, 12, 21, 28, 32	7, 9, 12, 27, 32	7, 10, 12, 28, 32	7, 12, 21, 28, 36	
PV location (@bus)/Size (kW)/p.f	–	@7/894/0.86	@29/938/0.87	@25/965/0.83	@8/947/0.84	@29/971/0.85	
Power loss (kW)	207.44	112.79	120.22	115.63	114.11	117.06	
Voltage sag (p.u)	0.96	0.42	0.50	0.47	0.46	0.51	
THD (%)	0.1835	0.1725	0.1793	0.1791	0.1731	0.1734	
Voltage unbalance (%)	4.92	2.07	2.26	2.28	2.14	2.17	
Energy not-supplied (MWh)	6.44	5.13	5.29	5.26	5.24	5.33	
SAIFI (failure/customer/yr)	1.36	1.18	1.22	1.23	1.21	1.31	
SAIDI (hour/customer/yr)	1.21	1.03	1.14	1.12	1.06	1.09	
MAIFI (mom/customer/yr)	4.64	2.38	2.47	2.44	2.45	2.56	
Objective function	1	0.6715	0.7163	0.7074	0.6894	0.7157	

Table 3 Statistical analysis of different algorithms performance to solve Scenario#1 for modified 33-bus network.

Method	Best	Mean	Worst	std	
IDO	0.6715	0.6738	0.6752	0.0331	
DO	0.7163	0.7185	0.7204	0.0402	
GWO	0.7074	0.7083	0.7097	0.0355	
PSO	0.6894	0.6909	0.6914	0.0367	
AOA	0.7157	0.7175	0.7196	0.0384	

Stochastic results for modified 33-bus network (Scenario#2)

This section depicts the power demand of all loads as a normal covariable with a mean of the nominal value and a standard deviation of 20% of the mean30. In order to simulate the uncertainty of wind power, a Weibull PDF is applied to the speed change curve's parameters. In the stochastic approach, the sample points are grouped into a clusters prior substantial number to applying the K-mean clustering to find the optimal agent for each cluster. Then, based on the gathered agents, transmission computations are performed. Then, according to the achieved agents, calculations of power transmission are performed. Using the IDO, the simulation results of the next scenario based on stochastic-MORWTA are presented, with WT generation and load uncertainties, and the performance of this optimization technique is compared to the characteristics of conventional DO, GWO, PSO, and AOA with the application of stochastic simulation based on 2m + 1 PEM and established clustering. In the network's fundamental state (without reconfiguration and WT allocation), the values of each objective, including losses, VS, THD, VU, ENS, SAIFI, SAIDI, and MAIFI, are calculated to be 218.43 kW, 1.02 p.u, 0.2001%, 5.40%, 6.63 MWh, 1.42 failure/customer/year, 1.26 h/customer/yr, and 5.03 mom/customer/yr, respectively. Figure 8 illustrates the convergence process of various stochastic-MORWTA problem-solving strategies. According to the convergence curves, the proposed IDO method has achieved a lower objective function value than other algorithms. In addition, it has a faster convergence rate and a lower tolerance for convergence than other problem-solving algorithms. As a result of solving the stochastic-MORWTA problem, the WT size has increased.Fig. 8 Convergence process to solve scenario#2 for modified 33-bus distribution network.

Table 4 illustrates the numerical results of the second stochastic-MORWTA scenario via the IDO and the conventional DO, GWO, PSO, and AOA. Assessing the results demonstrated that, in comparison to other algorithms, the value of each objective is decreased more than their base values. The IDO has installed 926 kW of WT generation with a p.f of 0.88 on bus 7 and determined the configuration of the 33-bus network by considering the opened lines as lines 7, 11, 28, 34, and 36. As a result, the values of power loss, VS, THD, VU, ENS, SAIFI, SAIDI, and MAIFI are determined to be 137.94 kW, 0.47 p.u, 0.1964%, 2.39%, 5.68 MWh, 1.30 failure/customer/yr, 1.15 h/customer/yr, and 2.75 mom/customer/yr, respectively. Furthermore, the statistical analysis depicted in Table 5 validates the superiority of the IDO in resolving the stochastic-MORWTA problem by achieving lower values for the Best, Mean, Worst, and std criteria (See Fig. 9). Also, the results of Table 5, showed the IDO has obtained lower convergence iteration (CI) similar to PSO, and convergence time (CT) to obtain the best solution. Table 4 Numerical results of scenario#2 optimization for modified 33-bus distribution network.

Method/Item	Base Net	IDO	DO	GWO	PSO	AOA	
Opened switches	33, 34, 35, 36, 37	7, 11, 28, 34, 36	10, 17, 28, 33, 34	7, 10, 32, 34, 37	5, 6, 10, 14, 15	10, 17, 28, 33, 34	
PV location (@bus)/Size (kW)/p.f	–	@7/926/0.88	@6/957/0.85	@7/936/0.86	@29/960/0.85	@25/981/0.84	
Power loss (kW)	218.43	137.94	147.46	140.78	139.44	144.64	
Voltage sag (p.u)	1.02	0.49	0.57	0.53	0.52	0.54	
THD (%)	0.2001	0.1964	0.1977	0.1974	0.1974	0.1975	
Voltage unbalance (%)	5.40	2.39	2.54	2.42	2.44	2.52	
Energy Not-supplied (MWh)	6.63	5.68	5.82	5.73	5.72	5.81	
SAIFI (failure/customer/yr)	1.42	1.30	1.38	1.35	1.35	1.38	
SAIDI (hour/customer/yr)	1.26	1.15	1.25	1.20	1.18	1.23	
MAIFI (mom/customer/yr)	5.03	2.75	2.86	2.81	2.79	2.83	
Objective function	1	0.7213	0.7645	0.7421	0.7382	0.7550	

Table 5 Statistical analysis of different algorithms performance to solve Scenario#2 for modified 33-bus network.

Method	Best	Mean	Worst	std	CI	CT (s)	
IDO	0.7213	0.7235	0.7254	0.0291	40	137	
DO	0.7645	0.7675	0.7689	0.0417	53	164	
GWO	0.7421	0.7449	0.7472	0.0335	48	153	
PSO	0.7382	0.7405	0.7426	0.0304	39	139	
AOA	0.7550	0.7572	0.7593	0.0368	42	161	
*CI, Convergence iteration; and CT, convergence time.

Fig. 9 Statistical analysis comparison of different algorithms to solve Scenario#2 for modified 33-bus network.

The overall objective function according to Eq. (1) has eight objectives, which are solved as a multi-objective structure based on the method of weight coefficients. If we solve the problem for each objective as a single objective and a total multi-objective problem, then we will face nine objective functions. Although the general objective of this study is to solve the multi-objective function (Eq. (1)). In Table 6, the results of the nonparametric test analysis of different algorithms for the 33-bus network are presented as Friedman's test rank, which shows the superiority of the proposed optimization method. Table 6 Friedman’s test rank summary of statistical analysis results for different algorithms for modified 33-bus network.

Objective function/algorithm	IDO	DO	GWO	PSO	AOA	
PLoss	1	5	3	2	4	
Vsag	1	5	2	3	4	
THD	1	5	2	3	4	
VUN	1	4	3	2	5	
ENS	1	5	3	2	4	
SAIFI	1	5	3	2	4	
SAIDI	1	4	3	2	5	
MAIFI	1	5	3	2	4	
F	1	5	3	2	4	
MFr	1	4.7778	2.7778	2.2222	4.2222	
Final rank	1	5	3	2	4	

Scenarios comparison for modified 33-bus network

In this part, the effectiveness of the MORWTA problem according to the IDO is contrasted with two deterministic and stochastic methods, illustrated in Table 7. Due to the WT generation and load uncertainties, the network configuration has been updated in addition to the change in installed WT generation in bus 7. The results demonstrated considering the WT generation and load uncertainties in the MORWTA problem in the form of stochastic-MORWTA, the value of each of the goals of power loss, VS, THD, VU, ENS, SAIFI, SAIDI, and MAIFI is increased 22.30%, 11.90%, 13.86%, 15.46%, 10.72%, 10.17%, 11.65% and 15.55%, respectively in contrast to the deterministic scenario without uncertainty (Deterministic-MORWTA). Table 7 Results comparison of Scenarios# 1, and 2 for modified 33-bus distribution network.

Method/Item	Deterministic-MORWTA	Stochastic-MORWTA	Change (%)	
Opened switches	10, 20, 28, 30, 35	7, 11, 28, 34, 36	–	
PV location (@bus)/Size (kW)/p.f	@7/894/0.86	@7/926/0.88	–	
Power loss (kW)	112.79	137.94	22.30 ↑	
Voltage sag (p.u)	0.42	0.49	16.67 ↑	
THD (%)	0.1725	0.1964	13.86 ↑	
Voltage unbalance (%)	2.07	2.39	15.46 ↑	
Energy not-supplied (MWh)	5.13	5.68	10.72 ↑	
SAIFI (failure/customer/yr)	1.18	1.30	10.17 ↑	
SAIDI (hour/customer/yr)	1.03	1.15	11.65 ↑	
MAIFI (mom/customer/yr)	2.38	2.75	15.55 ↑	
Objective function	0.6715	0.7213	7.42 ↑	
*↑ refers to the increasing.

Figure 10, depicts a visual comparison of each objective, including power loss, VS, THD, voltage imbalance, ENS, SAIFI, SAIDI, and MAIFI, for both deterministic- and stochastic-MORWTA scenarios. Based on Fig. 10, in the deterministic and stochastic scenarios, the power loss is 112.79 kW and 137.94 kW, the voltage sag is 0.42 p.u and 0.49 p.u, the THD is 1.725% and 1.647%, and the VU is 2.07% and 2.39%, respectively. ENS has increased from 5.13 MWh in the deterministic scenario to 5.68 MWh in the stochastic scenario. Moreover, the SAIFI, SAIDI, and MAIFI indices have increased from 1.18 failures per customer per year, 1.03 h per customer per year, and 2.38 under the deterministic scenario to 1.30 failures per customer per year, 1.15 h per customer per year, and 2.75 under the stochastic scenario. Consequently, under the conditions of contemplating the uncertainty of the network demand as well as the wind turbine's electricity production, the value of each objective has been raised, which is related to the uncertainty patterns that have been defined. As a result, the uncertainties considered based on the power changes of the load and the power fluctuations of the wind turbine have an effect on the system optimal configuration, and with the change of the size and generation factor of the WT, it led to significant changes in each of the investigated objectives. Consequently is crucial to account for uncertainties in the operation of distribution networks based on the reconfiguration and optimal sitting and sizing of WT so that network operators can know the exact values of the loss measures, power quality, and reliability and make the right decision. Otherwise, if the uncertainties of demand and power of renewable resources are not factored into the operation strategy, it is not practicable to achieve reliable network utilization and superior power quality.Fig. 10 Results comparison (a) Power loss, (b) VS and THD, (c) VU and ENS, (d) SAIFI, SAIDI, MAIFI for modified 33-bus network.

The variations in VS, VU, and THD of the buses are depicted in Figs. 11, 12 and 13. Clearly, the value of each of these goals has increased in the stochastic model considering uncertainty compared to the deterministic model. As depicted in Fig. 11, the value of VS in the stochastic-MORWTA model has increased in the majority of buses compared to the deterministic-MORWTA model (0.42–0.49 p.u). According to Fig. 12, the stochastic-MORWTA model has a higher value of voltage unbalance in all network circuits than the deterministic-MORWTA model (from 2.07 to 2.39%). In addition, Fig. 13 demonstrates that the value of THD in all network buses for the stochastic-MORWTA model has increased relative to the deterministic-MORWTA model (from 17.25 to 19.64%).Fig. 11 Voltage sag changes comparison using the IDO for modified 33-bus distribution network.

Fig. 12 VU changes comparison using the IDO for modified 33-bus distribution network.

Fig. 13 THD changes comparison using the IDO for modified 33-bus distribution network.

Deterministic, and stochastic results for unbalanced 25-bus network (Scenarios#1, and 2)

Figure 14 illustrates the schematic of the 25-bus network, consisting of 25 buses and 27 lines, with three tie-lines connecting buses 5–22, 8–15, and 12–25 (lines 25, 26, and 27). The network has an active load of 3239 kW and a reactive load of 2393 kW. In this simulation, the base voltage and base apparent power are set at 4.16 kV and 30 MVA, respectively. The data for loads and lines in this network are sourced from references56,57, while reliability data are taken from reference55. The following section presents the deterministic and stochastic results for the unbalanced 25-bus network.Fig. 14 Schematic of the unbalanced 25-bus distribution network.

Deterministic results for unbalanced 25-bus network (Scenario#1)

In the base unbalanced 25-bus network, lines 25, 26, and 27 are open. The outcomes for Scenario #1 of the 25-bus network are as follows: power loss is 150.12 kW (adapted from Ref.40,58), voltage stability (VS) is 0.8059 p.u., total harmonic distortion (THD) is 0.4835%, voltage unbalance (VU) is 0.4300%, energy not supplied (ENS) is 3.82 MWh, system average interruption frequency index (SAIFI) is 1.09 failures/customer/year, system average interruption duration index (SAIDI) is 0.88 h/customer/year, and momentary average interruption frequency index (MAIFI) is 3.27 moments/customer/year. In Scenario #1 (deterministic approach), solved using the IDO, switches 15, 17, and 22 are opened, and a wind turbine with a capacity of 611 kW and a power factor of 0.89 is installed at bus 7. As a result, the values for each objective are power loss of 75.42 kW, VS of 0.5546 p.u., THD of 0.3450%, VU of 0.2791%, ENS of 2.94 MWh, SAIFI of 0.97 failures/customer/year, SAIDI of 0.74 h/customer/year, and MAIFI of 2.02 moments/customer/year, demonstrating improvements in various objectives via the IDO.

Stochastic results for unbalanced 25-bus network (Scenario#2)

Figure 15 shows the convergence process for solving Scenario #2 of the unbalanced 25-bus distribution network using the proposed IDO, conventional DO, PSO, GWO, and AOA algorithms. The convergence curves indicate that the proposed IDO attained a lower objective function value (best solution) than the other algorithms and did so with a faster convergence rate.Fig. 15 Convergence process to solve scenario#2 for unbalanced 25-bus distribution network.

The numerical outcomes of Scenario #2 using the IDO, along with the conventional DO, GWO, PSO, and AOA algorithms are presented in Table 8. The evaluation of the results shows that, compared to other algorithms, each objective's value has decreased more than their base values. The IDO installed 673 kW of wind turbine generation with a power factor of 0.88 on bus 7 and determined the configuration of the 33-bus network by opening lines 13, 17, and 22. Consequently, the values for power loss, VS, VS, voltage unbalance (VU), ENS, SAIFI, SAIDI, and MAIFI are determined to be 89.26 kW, 0.6049 p.u., 0.3810%, 0.2801%, 3.29 MWh, 1.06 failures/customer/year, 0.83 h/customer/year, and 2.31 mom/customer/year, respectively. Additionally, the statistical analysis shown in Table 9 confirms the IDO's superiority in solving the stochastic-MORWTA problem by achieving lower values for the Best, Mean, Worst, and standard deviation criteria, as depicted in Fig. 16. The results in Table 9 also indicate that the IDO achieved lower CI and CT to find the best solution. Table 8 Numerical results of scenario#2 optimization for 25-bus distribution network.

Method/item	Base Net	IDO	DO	GWO	PSO	AOA	
Opened switches	25, 26, 27	13, 17, 22	14, 9, 22	13, 9, 22	6, 13, 17	14, 17, 20	
PV location (@bus)/Size (kW)/p.f	–/–	@7/673/0.88	@14/703/0.84	@7/685/0.85	@18/711/0.86	@14/692/0.85	
Power loss (kW)	157.85	89.26	95.03	92.88	92.13	93.47	
Voltage sag (p.u)	0.8440	0.6049	0.6064	0.6058	0.6054	0.6061	
THD (%)	0.5136	0.3810	0.3819	0.3813	0.3813	0.3817	
Voltage unbalance (%)	0.4501	0.2801	0.2810	0.2802	0.2803	0.2806	
Energy Not-supplied (MWh)	4.06	3.29	3.35	3.33	3.30	3.36	
SAIFI (failure/customer/yr)	1.13	1.06	1.14	1.11	1.09	1.14	
SAIDI (hour/customer/yr)	0.91	0.83	0.92	0.87	0.85	0.90	
MAIFI (mom/customer/yr)	3.44	2.31	2.45	2.39	2.36	2.44	
Objective function	1	0.7473	0.7807	0.7656	0.7580	0.7765	

Table 9 Statistical analysis of different algorithms performance to solve Scenario#2 for 25-bus network.

Method	Best	Mean	Worst	std	CI	CT (s)	
IDO	0.7473	0.7482	0.7488	0.0201	12	118	
DO	0.7807	0.7822	0.7860	0.0340	26	147	
GWO	0.7656	0.7665	0.7674	0.0284	21	132	
PSO	0.7580	0.7590	0.7597	0.0226	22	139	
AOA	0.7765	0.7773	0.7795	0.0319	41	155	
*CI, convergence iteration; and CT, convergence time.

Fig. 16 Statistical analysis comparison of different algorithms to solve Scenario#2 for 25-bus network.

In Table 10, the nonparametric test analysis results for different algorithms to solve the different objectives in the unbalanced 25-bus network are given as Friedman's test rank, which clears the better performance of the proposed IDO method. Table 10 Friedman’s test rank summary of statistical analysis results for different algorithms for 25-bus network.

Objective function/algorithm	IDO	DO	GWO	PSO	AOA	
PLoss	1	5	2	3	4	
Vsag	1	5	3	2	4	
THD	1	5	2	3	4	
VUN	1	4	3	2	5	
ENS	1	5	3	2	4	
SAIFI	1	4	3	2	5	
SAIDI	1	4	2	3	5	
MAIFI	1	5	3	2	4	
F	1	5	3	2	4	
MFr	1	4.6667	2.667	2.3333	4.3333	
Final rank	1	5	3	2	4	

Scenarios comparison for unbalanced 25-bus network

This section compares the effectiveness of the MORWTA problem using the IDO with two deterministic and stochastic methods, as shown in Table 11. Due to the uncertainties in wind turbine (WT) generation and load, the network configuration and the installed WT generation at bus 7 have been adjusted. The results indicate that when accounting for these uncertainties in the MORWTA problem, represented as stochastic-MORWTA, the objectives of power loss, VS, THD, VU, ENS, SAIFI, SAIDI, and MAIFI increased by 18.35%, 9.07%, 10.43%, 12.46%, 11.90%, 9.28%, 12.16%, and 14.36%, respectively, compared to the deterministic scenario without uncertainty, as illustrated in Fig. 17. Table 11 Results comparison of the deterministic- and stochastic- stochastic-MORWTA for unbalanced 25-bus network.

Method/item	Deterministic-MORWTA	Stochastic-MORWTA	Change (%)	
Opened switches	15, 17, 22	13, 17, 22	–	
WT location (@bus)/size (kW)/p.f	@7/611/0.89	@7/673/0.88	–	
Power loss (kW)	75.42	89.26	18.35 ↑	
Voltage sag (p.u)	0.5546	0.6049	9.07 ↑	
THD (%)	0.3450	0.3810	10.43 ↑	
Voltage Unbalance (%)	0.2791	0.2801	12.46 ↑	
Energy Not-supplied (MWh)	2.94	3.29	11.90 ↑	
SAIFI (failure/customer/yr)	0.97	1.06	9.28 ↑	
SAIDI (hour/customer/yr)	0.74	0.83	12.16 ↑	
MAIFI (mom/customer/yr)	2.02	2.31	14.36 ↑	
Objective Function	0.7089	0.7835	10.52 ↑	
*↑ refers to the increasing.

Fig. 17 % of increasing value of each objective in stochastic approach compared with the deterministic approach for unbalanced 25-bus network.

Figure 18 provides a visual comparison of each objective such as power loss, VS, THD, voltage imbalance, ENS, SAIFI, SAIDI, and MAIFI for both deterministic and stochastic MORWTA scenarios in the unbalanced 25-bus network. The comparison shows that under the conditions of considering the uncertainties in network demand and wind turbine electricity generation, the value of each objective has increased. This increase is attributed to the defined uncertainty patterns. Additionally, the percentage increase in each objective's value under the stochastic approach, compared to the deterministic approach, is presented based on the results in Table 11 for the unbalanced 25-bus network. This comparison highlights the impact of uncertainty on the network's performance, demonstrating the importance of accounting for stochastic variables in network optimization.Fig. 18 Results comparison (a) Power loss, (b) VS, THD, VU, ENS, SAIFI, SAIDI, and MAIFI for unbalanced 25-bus network.

Deterministic, and stochastic results for unbalanced 37-bus network (Scenarios#1, and 2)

The layout of the 37-bus network is shown in Fig. 19. This network is unbalanced and consists of 37 buses along with three tie-lines connecting buses 735–741, 701–722, and 727–732 (specifically lines 36, 37, and 38). The network has a total active load of 2472 kW and a reactive load of 1207 kVAr. For this simulation, the reference voltage and apparent power are set at 4.8 kV and 30 MVA, respectively. Data on the loads and lines for this network are sourced from Ref.59, while the reliability information is obtained from Ref.55.Fig. 19 Schematic of the unbalanced 37-bus distribution network.

Deterministic results for unbalanced 37-bus network (Scenario#1)

In the initial 37-bus network configuration, the objective values for power loss, VS, THD, VU, ENS, and reliability indices such as SAIFI, SAIDI, and MAIFI are calculated at 81.56 kW (based on Ref.59), 0.6352 p.u, 0.2075%, 0.13%, 2.186 MWh, 0.73 failures/customer/year, 0.61 h/customer/year, and 1.39 moments/customer/year, respectively. In Scenario #1, which employs a deterministic approach, the IDO is used, leading to the disconnection of three lines between buses 704–720, 708–732, and 734–737. Additionally, a wind turbine with a capacity of 548 kW and a power factor of 0.85 is installed at bus 722. As a result, the updated values for power loss, VS, THD, VU, ENS, SAIFI, SAIDI, and MAIFI are 34.56 kW, 0.4142 p.u, 0.1310%, 0.09%, 1.478 MWh, 0.63 failures/customer/year, 0.49 h/customer/year, and 1.20 moments/customer/year, respectively.

Stochastic results for unbalanced 37-bus network (Scenario#2)

The convergence process for solving Scenario #2 of the unbalanced 37-bus distribution network using the proposed IDO method, along with conventional DO, PSO, GWO, and AOA algorithms is presented in Fig. 20. The convergence curves indicate that the proposed IDO method not only reaches a lower objective function value (best solution) compared to the other algorithms but also does so at a faster rate.Fig. 20 Convergence process to solve scenario#2 for unbalanced 37-bus distribution network.

The numerical results for Scenario #2 of the 37-bus network are as follows: power loss is obtained 83.17 kW, voltage stability (VS) is 0.6755 p.u, THD is 0.2183%, VU is 0.139%, ENS is 2.257 MWh, SAIFI is 0.76 failures/customer/year, SAIDI is 0.64 h/customer/year, and MAIFI is 1.48 moments/customer/year.

For Scenario #2, which employs a stochastic approach, the IDO method is applied, and according to Table 12, three lines between buses 704–720, 708–732, and 710–734 are disconnected, and a wind turbine with a 579 kW capacity and a power factor of 0.83 is installed at bus 706. With these changes, the updated values are: power loss is 38.78 kW, VS is 0.4583 p.u, THD is 0.1472%, VU is 0.101%, ENS is 1.689 MWh, SAIFI is 0.71 failures/customer/year, SAIDI is 0.55 h/customer/year, and MAIFI is 1.36 moments/customer/year. Additionally, statistical analysis in Table 13 demonstrates the effectiveness of the IDO method in addressing the stochastic-MORWTA problem. The IDO method achieves lower values for Best, Mean, Worst, and standard deviation criteria compared to other algorithms, as shown in Fig. 21. Table 13 also indicates that the IDO method requires fewer CI and less CT to reach the optimal solution. Table 12 Numerical results of scenario#2 optimization for unbalanced 37-bus distribution network.

Method/Item	Base Net	IDO	DO	GWO	PSO	AOA	
Opened switches	735–741, 701–722, 727–732	704–720, 708–732, 710–734	704–713, 708–709, 710–734	704–720, 708–732, 737–738	707–720, 708–732, 710–734	704–720, 709–730, 734–737	
WT location (@bus)/Size (kW)/p.f	–	@706/579/0.83	@707/566/0.84	@737/559/0.84	@704/583/0.82	@720/573/0.82	
Power loss (kW)	83.17	38.78	39.71	39.54	39.15	39.58	
Voltage sag (p.u)	0.6755	0.4583	0.4596	0.4591	0.4587	0.4594	
THD (%)	0.2183	0.1472	0.1478	0.1475	0.1473	0.1478	
Voltage Unbalance (%)	0.139	0.101	0.101	0.101	0.101	0.101	
Energy Not-supplied (MWh)	2.257	1.689	1.697	1.693	1.691	1.695	
SAIFI (failure/customer/yr)	0.76	0.71	0.75	0.74	0.72	0.75	
SAIDI (hour/customer/yr)	0.64	0.55	0.62	0.58	0.57	0.60	
MAIFI (mom/customer/yr)	1.48	1.36	1.42	1.41	1.39	1.42	
Objective function	1	0.7505	0.7787	0.7675	0.7599	0.7745	

Table 13 Statistical analysis of different algorithms performance to solve Scenario#2 for unbalanced 37-bus network.

Method	Best	Mean	Worst	std	CI	CT (s)	
IDO	0.7505	0.7514	0.7523	0.0274	10	126	
DO	0.7787	0.7804	0.7825	0.0375	25	140	
GWO	0.7675	0.7686	0.7698	0.0313	25	138	
PSO	0.7599	0.7612	0.7621	0.0286	22	133	
AOA	0.7745	0.7759	0.7778	0.0370	20	151	
* CI: Convergence Iteration, and CT: Convergence Time.

Fig. 21 Statistical analysis comparison of different algorithms to solve Scenario#2 for 37-bus network.

In Table 14, the Friedman's test rank results as a nonparametric test analysis are presented for different algorithms for unbalanced 37-bus network and the obtained results are proved the superior capability of the IDO to solve different objectives compared with the other algorithms. Table 14 Friedman’s test rank summary of statistical analysis results for different algorithms for unbalanced 37-bus network.

Objective function/algorithm	IDO	DO	GWO	PSO	AOA	
PLoss	1	5	3	2	4	
Vsag	2	4	3	1	5	
THD	1	5	3	2	4	
VUN	1	5	3	2	4	
ENS	1	5	3	2	4	
SAIFI	1	5	2	3	4	
SAIDI	1	4	3	2	5	
MAIFI	1	5	3	2	4	
F	1	5	3	2	4	
MFr	1.2222	4.7778	2.8889	2.0000	4.2222	
Final rank	1	5	3	2	4	

Scenarios comparison for unbalanced 37-bus network

In this section, the performance of the IDO in addressing the MORWTA problem is compared with two other approaches, deterministic and stochastic, as shown in Table 15. The network configuration has been adjusted to account for uncertainties in WT generation and load, including changes in the WT generation at bus 7. The results reveal that incorporating WT generation and load uncertainties into the MORWTA problem, referred to as stochastic-MORWTA, leads to increases in the objectives of power loss, VS, THD, VU, ENS, SAIFI, SAIDI, and MAIFI by 12.21%, 10.64%, 12.37%, 9.82%, 14.30%, 12.65%, 12.63%, and 13.89%, respectively, compared to the deterministic scenario without uncertainties, as depicted in Fig. 22. Table 15 Results comparison of the deterministic- and stochastic- stochastic-MORWTA for unbalanced 37-bus network.

Method/Item	Deterministic-MORWTA	Stochastic-MORWTA	Change (%)	
Opened Switches	704–720, 708–732, 734–737	704–720, 708–732, 710–734	–	
WT Location (@bus)/Size (kW)/p.f	@722/548/0.85	@706/579/0.83	–	
Power Loss (kW)	34.56	38.78	12.21 ↑	
Voltage Sag (p.u)	0.4142	0.4583	10.64 ↑	
THD (%)	0.1310	0.1472	12.37 ↑	
Voltage Unbalance (%)	0.09	0.101	9.82 ↑	
Energy Not-supplied (MWh)	1.478	1.689	14.30 ↑	
SAIFI (failure/customer/yr)	0.63	0.71	12.65 ↑	
SAIDI (hour/customer/yr)	0.49	0.55	12.63 ↑	
MAIFI

(mom/customer/yr)

	1.20	1.36	13.89 ↑	
Objective function	0.7007	0.7505	7.11 ↑	
*↑ refers to the increasing.

Fig. 22 % of increasing value of each objectives in stochastic approach compared with the deterministic approach for unbalanced 37-bus network.

A visual comparison of each objective—such as power loss, VS, THD, voltage imbalance, ENS, SAIFI, SAIDI, and MAIFI—for both deterministic and stochastic-MORWTA scenarios in the unbalanced 37-bus network is depicted in Fig. 23. The results show that when accounting for uncertainties in network demand and wind turbine electricity production, the values for each objective increase according to the defined uncertainty patterns. Without incorporating these uncertainties into the operational strategy, achieving reliable network performance and high power quality is impractical. Additionally, the percentage increase of the stochastic approach compared to the deterministic method for the unbalanced 37-bus network is detailed in Table 15.Fig. 23 Results comparison (a) Power loss, (b) VS, THD, VU, ENS, SAIFI, SAIDI, and MAIFI for unbalanced 37-bus network.

Model validation with previous studies

Comparison for the modified 33-bus network

The numerical results of the deterministic-MORWTA method are compared with those from previous studies15,40,60,61. Study60 introduces the reconfiguration of unbalanced distribution networks using a fuzzy multi-criteria approach and the new improved coronavirus herd immunity optimizer algorithm (ICHIOA). The focus is on reducing losses, VU, VS, and energy that is not provided for customers. Voltage imbalances and VS are considered power quality metrics, and ENS represents energy not delivered. In61, the network reconfiguration is optimized for minimizing losses, VS, VU, and ENS on an unbalanced 33-bus network using an artificial electric field algorithm-pattern search (AEFAPS). Additionally, studies15,40 implement the 33-bus network reconfiguration to minimize power losses and improve power quality indices. In Table 16, the performance of IDO is compared with ICHIOA and AEFAPS. The results show that IDO achieves lower losses, VS, VU, and ENS, demonstrating the superior performance of the proposed method. Table 16 Comparison with previous studies for the modified 33-bus network.

Method/item	IDO Based-MORWTA	15	DTLBO20	40	ICHIOA60	AEFAPS61	
Opened switches	10, 20, 28, 30, 35	7, 9, 14, 32, 37	7, 9, 14, 28, 32	7, 9, 14, 32, 37	9, 13, 16, 28, 33	7, 11, 13, 28, 32	
WT location (@bus)/Size (kW)/p.f	@7/894/0.86	–	@3/50/0.80, @6/100/0.90, @24/200/0.90, @29/100/0.90	–	–	–	
Power loss (kW)	112.79	139.50	115.70	144.48	167.371	146.96	
Voltage sag (p.u)	0.42	–	–	0.91	0.665	0.653	
THD (%)	0.1725	–	–	0.1753	–	–	
Voltage Unbalance (%)	2.07	–	–	4.79	3.379	2.42	
Energy Not-supplied (MWh)	5.13	–	–	–	6.421	6.69	
SAIFI (failure/customer/yr)	1.18	–	–	–	–	–	
SAIDI (hour/customer/yr)	1.03	–	–	–	–	–	
MAIFI (mom/customer/yr)	2.38	–	–	–	–	–	

Comparison for the unbalanced 25-bus network

A comparison of the deterministic and stochastic numerical results obtained in this study is presented in Table 7 for the unbalanced 25-bus distribution network with those from Refs.40,56,62. In56, the reconfiguration of the unbalanced 25-bus network focuses on minimizing power losses, while in62, both DG allocation and network reconfiguration are carried out simultaneously. According to Table 17, the proposed method in our study achieves lower power losses in both deterministic and stochastic approaches compared to Refs.56,62. Table 17 Comparison with previous studies for the unbalanced 25-bus network.

Method/item	IDO based-MORWTA (Deterministic)	IDO based-MORWTA (Stochastic)	40	56	62	
Opened switches	15, 17, 22	13, 17, 22	17, 22, 27	15, 17, 22	15, 17, 22	
WT location (@bus)/size (kW)/p.f	@7/611/0.89	@7/673/0.88	@13/215/0.85	–	@12/100/0.9

@22/100/0.9

	
Power loss (kW)	75.42	89.26	85.525	136.13	113.72	

Comparison for the unbalanced 37-bus network

The numerical results of the deterministic and stochastic MORWTA methods are compared in Table 18 with the results of Ref.19, which utilized a hybrid honey bee colony algorithm (HHBCA) to allocate capacitors in the unbalanced 37-bus distribution network. The objective in Ref.19 was to minimize power system losses and unbalances while maintaining voltage levels and THD at the buses. As shown in Table 18, the proposed methodology in this study achieves lower values of power loss, THD, and voltage unbalance. Table 18 Comparison with previous studies for the unbalanced 37-bus network.

Method/Item	IDO Based-MORWTA (Deterministic)	IDO Based-MORWTA (Stochastic)	19	
Opened switches	704–720, 708–732, 734–737	704–720, 708–732, 710–734	–	
WT location (@bus)/size (kW)/p.f	@722/548/0.85	@706/579/0.83	900kVAr Capacitor	
Power loss (kW)	34.56	38.78	53.42	
THD (%)	0.1310	0.1472	0.3550	
Voltage unbalance (%)	0.09	0.101	0.107	

Conclusion

In this study, a solving approach based on the MORWTA is proposed for minimizing the losses, VS, THD, VU, ENS, SAIFI, SAIDI, and MAIFI in three modified 33-bus, and unbalanced 25-, and 37-bus distribution networks by employing the IDO. For stochastic problem-solving and modeling, the 2m + 1 PEM integrated with clustering was employed. Two deterministic and stochastic models without and with uncertainties, respectively, were applied to the problem. The findings of the investigation are as follows:The results cleared that the MORWTA improved all objectives by simultaneously reconfiguring the distribution networks and allocating wind turbines in comparison to the network's initial state. Consequently, the simultaneous approach of determining the network configuration optimally and the optimal injection of wind power into the network has enhanced the efficiency of the network and met all of the objectives.

For the modified 33-bus network, in the deterministic model utilizing the IDO, the losses, VS, THD, VU, ENS, SAIFI, and SAIDI were decreased by 45.63%, 56.25%, 5.99%, 57.93%, 20.34%, 13.24%, 14.88%, and 48.71%, respectively, compared to the base network.

The statistical analysis demonstrates that the IDO-solved deterministic-MORWTA had produced more effective results than the conventional DO, GWO, PSO, and AOA techniques by achieving superior statistical efficacy criteria for all networks. The comparison of deterministic- and stochastic MORWTA results confirms the superiority of the IDO by achieving a higher objective value and lower values for the Best, Mean, Worst, and standard criteria.

For the modified 33-bus network, comparison of results revealed that the power loss, VS, THD, VU, ENS, SAIFI, SAIDI, and MAIFI were increased by 20.52%, 11.90%, 13.86%, 15.46%, 10.72%, 10.17%, 11.65% and 15.55%, respectively, when the stochastic model is used in place of the deterministic model without uncertainty consideration.

For the unbalanced 25-bus network, the results demonstrated considering the WT generation and load uncertainties in the MORWTA problem in the form of stochastic-MORWTA, the value of each of the goals of power loss, VS, THD, VU, ENS, SAIFI, SAIDI, and MAIFI is increased 18.35%, 9.07%, 10.43%, 12.46%, 11.90%, 9.28%, 12.16% and 14.36%, respectively in contrast to the deterministic scenario without uncertainty.

For the unbalanced 37-bus network, the results demonstrated considering the WT generation and load uncertainties in the MORWTA problem in the form of stochastic-MORWTA, the value of each of the goals of power loss, VS, THD, VU, ENS, SAIFI, SAIDI, and MAIFI is increased 12.21%, 10.64%, 12.37%, 9.82%, 14.30%, 12.65%, 12.63% and 13.89%, respectively in contrast to the deterministic scenario without uncertainty.

The results of the stochastic model demonstrated that the 2m + 1 PEM with clustering was an efficient method for modeling uncertainties and resolving the stochastic MORWTA problem. In addition, the results cleared that the uncertainties altered both the network configuration and the wind turbine capacity injected into the network, which had a significant impact on the value of each objective.

In order for network operators to know the exact values of power losses, power quality, and reliability indices and make the correct decision, it is crucial to account for uncertainties in operation of the network via the MORWTA. In the absence of this, it is impossible to provide the network with reliable performance and optimal power quality in the presence of uncertainties.

The robust unbalanced network reconfiguration with hybrid photovoltaic/wind energy system allocation utilizing the information gap decision theory are recommended for future research.

List of symbols

ALO Ant lion optimizer

AOA Arithmetic optimization algorithm

BWOA Beluga whale optimization algorithm

CGFA Chaotic golden flower algorithm

DG Distributed generation

DO Dandelion optimization

DTLBO Discrete teaching learning-based optimization

DVTLs¯ Decision vectors of network load

DVWT¯ Decision vectors of wind turbine

ENS Energy not-supplied

GA Genetic algorithm

GWO Gray wolf optimization

IGJO Improved golden jackal optimization

IMTBO Iterative mountaineering team-based optimization

IRV Input random variable

MAIFI Momentary average interruption frequency index

MCS Monte Carlo Simulation

MFO Moth-flame optimization

IBWO Improved black widow optimization

MORWTA Multi-objective reconfiguration and wind turbine Allocation

NGO Northern Goshawk Optimization

NFL No-free launch

ORV Output random variable

PEM Point estimate method

p.fWT Wind turbine power factor

Ploss Active power loss

PPost Active post power

PSO Particle swarm optimization

PWT Active wind power

PWGMax Upper limit for wind active generation

PWGMin Lower limit for wind active generation

Qloss Rective power loss

QPost Rective post power

QWT Rective wind power

SAIDI System average interruption duration index

SAIFI System average interruption frequency index

SHO Spotted hyena optimizer

THD Total harmonic disorder

VOLd,i Distortion component of the bus voltage

VU Voltage unbalance

VS Voltage sag

WT Wind turbine

Xelite Optimum location of the dandelion seed in the ith iteration

Xmean-t Average position of the population at ith iteration

∅ Adaptive inertia weight

β Overall number of customers

Acknowledgements

The authors extend their appreciation to the Deanship of Scientific Research at Northern Border Univer-sity, Arar, KSA for funding this research work through the project number “NBU-FFR-2024-2200-02“.

Author contributions

F Duan: investigation, visualization, writing—original draft preparation, A. Basem: reviewing and editing original draft, formal analysis, D. J. Jasim: visualization, writing—original draft preparation, project administration, S. Belhaj: reviewing and editing original draft, formal analysis, M. Eslami: supervision, software, data curation, resources, validation, formal analysis, software, M Khajehzadeh: conceptualization, reviewing and editing original draft, formal analysis, S. Palani: resources, investigation, visualization, writing—original draft preparation, project administration, funding acquisition.

Data availability

The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Ali ZM Diaaeldin IM El-Rafei A Hasanien HM Aleem SHA Abdelaziz AY A novel distributed generation planning algorithm via graphically-based network reconfiguration and soft open points placement using Archimedes optimization algorithm Ain Shams Eng. J. 2021 12 2 1923 1941 10.1016/j.asej.2020.12.006
Ali, Z. M. et al. A novel distributed generation planning algorithm via graphically-based network reconfiguration and soft open points placement using Archimedes optimization algorithm. Ain Shams Eng. J. 12(2), 1923–1941 (2021).10.1016/j.asej.2020.12.006
2. Baker, P. Challenges facing distribution system operators in a decarbonised power system. Regulatory Assistance Project (2020).
3. Caballero-Peña J Cadena-Zarate C Parrado-Duque A Osma-Pinto G Distributed energy resources on distribution networks: A systematic review of modelling, simulation, metrics, and impacts Int. J. Electr. Power Energy Syst. 2022 138 107900 10.1016/j.ijepes.2021.107900
Caballero-Peña, J., Cadena-Zarate, C., Parrado-Duque, A. & Osma-Pinto, G. Distributed energy resources on distribution networks: A systematic review of modelling, simulation, metrics, and impacts. Int. J. Electr. Power Energy Syst. 138, 107900 (2022).10.1016/j.ijepes.2021.107900
4. Mrzljak V Poljak I Mrakovčić T Energy and exergy analysis of the turbo-generators and steam turbine for the main feed water pump drive on LNG carrier Energy Convers. Manage. 2017 140 307 323 10.1016/j.enconman.2017.03.007
Mrzljak, V., Poljak, I. & Mrakovčić, T. Energy and exergy analysis of the turbo-generators and steam turbine for the main feed water pump drive on LNG carrier. Energy Convers. Manage. 140, 307–323 (2017).10.1016/j.enconman.2017.03.007
5. Davoudkhani IF Zishan F Mansouri S Abdollahpour F Grisales-Noreña LF Montoya OD Allocation of renewable energy resources in distribution systems while considering the uncertainty of wind and solar resources via the multi-objective Salp Swarm algorithm Energies 2023 16 1 474 10.3390/en16010474
Davoudkhani, I. F. et al. Allocation of renewable energy resources in distribution systems while considering the uncertainty of wind and solar resources via the multi-objective Salp Swarm algorithm. Energies 16(1), 474 (2023).10.3390/en16010474
6. Nezhad EH Ebrahimi R Ghanbari M Fuzzy Multi-objective allocation of photovoltaic energy resources in unbalanced network using improved manta ray foraging optimization algorithm Expert Syst. Appl. 2023 234 121048 10.1016/j.eswa.2023.121048
Nezhad, E. H., Ebrahimi, R. & Ghanbari, M. Fuzzy Multi-objective allocation of photovoltaic energy resources in unbalanced network using improved manta ray foraging optimization algorithm. Expert Syst. Appl. 234, 121048 (2023).10.1016/j.eswa.2023.121048
7. Ali ZM Diaaeldin IM Abdel Aleem SHE El-Rafei A Abdelaziz AY Jurado F Scenario-based network reconfiguration and renewable energy resources integration in large-scale distribution systems considering parameters uncertainty Mathematics 2020 9 1 26 10.3390/math9010026
Ali, Z. M. et al. Scenario-based network reconfiguration and renewable energy resources integration in large-scale distribution systems considering parameters uncertainty. Mathematics 9(1), 26 (2020).10.3390/math9010026
8. Fathi R Tousi B Galvani S Allocation of renewable resources with radial distribution network reconfiguration using improved salp swarm algorithm Appl. Soft Comput. 2023 132 109828 10.1016/j.asoc.2022.109828
Fathi, R., Tousi, B. & Galvani, S. Allocation of renewable resources with radial distribution network reconfiguration using improved salp swarm algorithm. Appl. Soft Comput. 132, 109828 (2023).10.1016/j.asoc.2022.109828
9. Ghadi-Sahebi M Ebrahimi R Parvin-darabad V Optimal probabilistic operation management of smart parking lot and renewable sources in microgrid to reduce cost and improve system reliability considering demand response program Int. Trans. Electr. Energy Syst. 2021 31 12 13108 10.1002/2050-7038.13108
Ghadi-Sahebi, M., Ebrahimi, R. & Parvin-darabad, V. Optimal probabilistic operation management of smart parking lot and renewable sources in microgrid to reduce cost and improve system reliability considering demand response program. Int. Trans. Electr. Energy Syst. 31(12), 13108 (2021).10.1002/2050-7038.13108
10. Moghaddam MJH Kalam A Shi J Nowdeh SA Gandoman FH Ahmadi A A new model for reconfiguration and distributed generation allocation in distribution network considering power quality indices and network losses IEEE Syst. J. 2020 14 3 3530 3538 10.1109/JSYST.2019.2963036
Moghaddam, M. J. H. et al. A new model for reconfiguration and distributed generation allocation in distribution network considering power quality indices and network losses. IEEE Syst. J. 14(3), 3530–3538 (2020).10.1109/JSYST.2019.2963036
11. Huy THB Enhancing distribution system performance via distributed generation placement and reconfiguration based on improved symbiotic organisms search J. Control Sci. Eng. 2023 2023 1 10 10.1155/2023/6081991
Huy, T. H. B. Enhancing distribution system performance via distributed generation placement and reconfiguration based on improved symbiotic organisms search. J. Control Sci. Eng. 2023, 1–10 (2023).10.1155/2023/6081991
12. Vai V Suk S Lorm R Chhlonh C Eng S Bun L Optimal reconfiguration in distribution systems with distributed generations based on modified sequential switch opening and exchange Appl. Sci. 2021 11 5 2146 10.3390/app11052146
Vai, V. et al. Optimal reconfiguration in distribution systems with distributed generations based on modified sequential switch opening and exchange. Appl. Sci. 11(5), 2146 (2021).10.3390/app11052146
13. Nowdeh SA Naderipour A Davoudkhani IF Guerrero JM Stochastic optimization–based economic design for a hybrid sustainable system of wind turbine, combined heat, and power generation, and electric and thermal storages considering uncertainty: A case study of Espoo, Finland Renew. Sustain. Energy Rev. 2023 183 113440 10.1016/j.rser.2023.113440
Nowdeh, S. A., Naderipour, A., Davoudkhani, I. F. & Guerrero, J. M. Stochastic optimization–based economic design for a hybrid sustainable system of wind turbine, combined heat, and power generation, and electric and thermal storages considering uncertainty: A case study of Espoo, Finland. Renew. Sustain. Energy Rev. 183, 113440 (2023).10.1016/j.rser.2023.113440
14. Taheri S Akbari A Ghahremani B Razban A Reliability-based energy scheduling of active buildings subject to renewable energy and demand uncertainty Therm. Sci. Eng. Progr. 2022 28 101149 10.1016/j.tsep.2021.101149
Taheri, S., Akbari, A., Ghahremani, B. & Razban, A. Reliability-based energy scheduling of active buildings subject to renewable energy and demand uncertainty. Therm. Sci. Eng. Progr. 28, 101149 (2022).10.1016/j.tsep.2021.101149
15. Mahdavi M Alhelou HH Hesamzadeh MR An efficient stochastic reconfiguration model for distribution systems with uncertain loads IEEE Access 2022 10 10640 10652 10.1109/ACCESS.2022.3144665
Mahdavi, M., Alhelou, H. H. & Hesamzadeh, M. R. An efficient stochastic reconfiguration model for distribution systems with uncertain loads. IEEE Access 10, 10640–10652 (2022).10.1109/ACCESS.2022.3144665
16. Kahouli O Alsaif H Bouteraa Y Ben Ali N Chaabene M Power system reconfiguration in distribution network for improving reliability using genetic algorithm and particle swarm optimization Appl. Sci. 2021 11 7 3092 10.3390/app11073092
Kahouli, O., Alsaif, H., Bouteraa, Y., Ben Ali, N. & Chaabene, M. Power system reconfiguration in distribution network for improving reliability using genetic algorithm and particle swarm optimization. Appl. Sci. 11(7), 3092 (2021).10.3390/app11073092
17. Nguyen TT Truong AV Phung TA A novel method based on adaptive cuckoo search for optimal network reconfiguration and distributed generation allocation in distribution network Int. J. Electri. Power Energy Syst. 2016 78 801 815 10.1016/j.ijepes.2015.12.030
Nguyen, T. T., Truong, A. V. & Phung, T. A. A novel method based on adaptive cuckoo search for optimal network reconfiguration and distributed generation allocation in distribution network. Int. J. Electri. Power Energy Syst. 78, 801–815 (2016).10.1016/j.ijepes.2015.12.030
18. Nguyen TT Truong AV Distribution network reconfiguration for power loss minimization and voltage profile improvement using cuckoo search Int. J. Electri. Power Energy Syst. 2015 68 233 242 10.1016/j.ijepes.2014.12.075
Nguyen, T. T. & Truong, A. V. Distribution network reconfiguration for power loss minimization and voltage profile improvement using cuckoo search. Int. J. Electri. Power Energy Syst. 68, 233–242 (2015).10.1016/j.ijepes.2014.12.075
19. Taher SA Bagherpour R A new approach for optimal capacitor placement and sizing in unbalanced distorted distribution systems using hybrid honey bee colony algorithm Int. J. Electri. Power Energy Syst. 2013 49 430 448 10.1016/j.ijepes.2013.02.003
Taher, S. A. & Bagherpour, R. A new approach for optimal capacitor placement and sizing in unbalanced distorted distribution systems using hybrid honey bee colony algorithm. Int. J. Electri. Power Energy Syst. 49, 430–448 (2013).10.1016/j.ijepes.2013.02.003
20. Lotfipour A Afrakhte H A discrete teaching–learning-based optimization algorithm to solve distribution system reconfiguration in presence of distributed generation Int. J. Electri. Power Energy Syst. 2016 82 264 273 10.1016/j.ijepes.2016.03.009
Lotfipour, A. & Afrakhte, H. A discrete teaching–learning-based optimization algorithm to solve distribution system reconfiguration in presence of distributed generation. Int. J. Electri. Power Energy Syst. 82, 264–273 (2016).10.1016/j.ijepes.2016.03.009
21. Swarnkar A Gupta N Niazi KR A novel codification for meta-heuristic techniques used in distribution network reconfiguration Electric Power Syst. Res. 2011 81 7 1619 1626 10.1016/j.epsr.2011.03.020
Swarnkar, A., Gupta, N. & Niazi, K. R. A novel codification for meta-heuristic techniques used in distribution network reconfiguration. Electric Power Syst. Res. 81(7), 1619–1626 (2011).10.1016/j.epsr.2011.03.020
22. Pereira EC Barbosa CH Vasconcelos JA Distribution network reconfiguration using iterative branch exchange and clustering technique Energies 2023 16 5 2395 10.3390/en16052395
Pereira, E. C., Barbosa, C. H. & Vasconcelos, J. A. Distribution network reconfiguration using iterative branch exchange and clustering technique. Energies 16(5), 2395 (2023).10.3390/en16052395
23. Zhu M Arabi Nowdeh S Daskalopulu A An improved human-inspired algorithm for distribution network stochastic reconfiguration using a multi-objective intelligent framework and unscented transformation Mathematics 2023 11 17 3658 10.3390/math11173658
Zhu, M., Arabi Nowdeh, S. & Daskalopulu, A. An improved human-inspired algorithm for distribution network stochastic reconfiguration using a multi-objective intelligent framework and unscented transformation. Mathematics 11(17), 3658 (2023).10.3390/math11173658
24. Swaminathan D Rajagopalan A Montoya OD Arul S Grisales-Noreña LF Distribution network reconfiguration based on hybrid golden flower algorithm for smart cities evolution Energies 2023 16 5 2454 10.3390/en16052454
Swaminathan, D., Rajagopalan, A., Montoya, O. D., Arul, S. & Grisales-Noreña, L. F. Distribution network reconfiguration based on hybrid golden flower algorithm for smart cities evolution. Energies 16(5), 2454 (2023).10.3390/en16052454
25. Liu KY Sheng W Liu Y Meng X A network reconfiguration method considering data uncertainties in smart distribution networks Energies 2017 10 5 618 10.3390/en10050618
Liu, K. Y., Sheng, W., Liu, Y. & Meng, X. A network reconfiguration method considering data uncertainties in smart distribution networks. Energies 10(5), 618 (2017).10.3390/en10050618
26. Alqahtani M Marimuthu P Moorthy V Pangedaiah B Reddy CR Kiran Kumar M Khalid M Investigation and minimization of power loss in radial distribution network using gray wolf optimization Energies 2023 16 12 4571 10.3390/en16124571
Alqahtani, M. et al. Investigation and minimization of power loss in radial distribution network using gray wolf optimization. Energies 16(12), 4571 (2023).10.3390/en16124571
27. Quintana E Inga E Optimal reconfiguration of electrical distribution system using heuristic methods with geopositioning constraints Energies 2022 15 15 5317 10.3390/en15155317
Quintana, E. & Inga, E. Optimal reconfiguration of electrical distribution system using heuristic methods with geopositioning constraints. Energies 15(15), 5317 (2022).10.3390/en15155317
28. Sellami R Khenissi I Guesmi T Alshammari BM Alqunun K Alshammari AS Tlijani K Neji R Optimal reconfiguration of distribution network considering stochastic wind energy and load variation using hybrid SAMPSO optimization method Sustainability 2022 14 18 11208 10.3390/su141811208
Sellami, R. et al. Optimal reconfiguration of distribution network considering stochastic wind energy and load variation using hybrid SAMPSO optimization method. Sustainability 14(18), 11208. 10.3390/su141811208 (2022).10.3390/su141811208
29. Alanazi A Alanazi TI Multi-objective framework for optimal placement of distributed generations and switches in reconfigurable distribution networks: An improved particle swarm optimization approach Sustainability 2023 15 11 9034 10.3390/su15119034
Alanazi, A. & Alanazi, T. I. Multi-objective framework for optimal placement of distributed generations and switches in reconfigurable distribution networks: An improved particle swarm optimization approach. Sustainability 15(11), 9034 (2023).10.3390/su15119034
30. Yan X Zhang Q Research on combination of distributed generation placement and dynamic distribution network reconfiguration based on MIBWOA Sustainability 2023 15 12 9580 10.3390/su15129580
Yan, X. & Zhang, Q. Research on combination of distributed generation placement and dynamic distribution network reconfiguration based on MIBWOA. Sustainability 15(12), 9580 (2023).10.3390/su15129580
31. Yang J Xiong J Chen YL Yee PL Ku CS Babanezhad M Improved golden jackal optimization for optimal allocation and scheduling of wind turbine and electric vehicles parking lots in electrical distribution network using Rosenbrock’s direct rotation strategy Mathematics 2023 11 6 1415 10.3390/math11061415
Yang, J. et al. Improved golden jackal optimization for optimal allocation and scheduling of wind turbine and electric vehicles parking lots in electrical distribution network using Rosenbrock’s direct rotation strategy. Mathematics 11(6), 1415 (2023).10.3390/math11061415
32. Alanazi A Alanazi M Nowdeh SA Abdelaziz AY Abu-Siada A Stochastic-metaheuristic model for multi-criteria allocation of wind energy resources in distribution network using improved equilibrium optimization algorithm Electronics 2022 11 20 3285 10.3390/electronics11203285
Alanazi, A., Alanazi, M., Nowdeh, S. A., Abdelaziz, A. Y. & Abu-Siada, A. Stochastic-metaheuristic model for multi-criteria allocation of wind energy resources in distribution network using improved equilibrium optimization algorithm. Electronics 11(20), 3285 (2022).10.3390/electronics11203285
33. Rehman N Mufti MUD Gupta N Metaheuristic method for a wind-integrated distribution network to support voltage stabilisation employing electric vehicle loads Appl. Sci. 2023 13 4 2254 10.3390/app13042254
Rehman, N., Mufti, M. U. D. & Gupta, N. Metaheuristic method for a wind-integrated distribution network to support voltage stabilisation employing electric vehicle loads. Appl. Sci. 13(4), 2254 (2023).10.3390/app13042254
34. Ali ES Abd Elazim SM Hakmi SH Mosaad MI Optimal allocation and size of renewable energy sources as distributed generations using shark optimization algorithm in radial distribution systems Energies 2023 16 10 3983 10.3390/en16103983
Ali, E. S., Abd Elazim, S. M., Hakmi, S. H. & Mosaad, M. I. Optimal allocation and size of renewable energy sources as distributed generations using shark optimization algorithm in radial distribution systems. Energies 16(10), 3983 (2023).10.3390/en16103983
35. Aref M Oboskalov V El-Shahat A Abdelaziz AY Modified analytical technique for multi-objective optimal placement of high-level renewable energy penetration connected to egyptian power system Mathematics 2023 11 4 958 10.3390/math11040958
Aref, M., Oboskalov, V., El-Shahat, A. & Abdelaziz, A. Y. Modified analytical technique for multi-objective optimal placement of high-level renewable energy penetration connected to egyptian power system. Mathematics 11(4), 958 (2023).10.3390/math11040958
36. Naderipour A Abdul-Malek Z Mustafa MWB Guerrero JM A multi-objective artificial electric field optimization algorithm for allocation of wind turbines in distribution systems Appl. Soft Comput. 2021 105 107278 10.1016/j.asoc.2021.107278
Naderipour, A., Abdul-Malek, Z., Mustafa, M. W. B. & Guerrero, J. M. A multi-objective artificial electric field optimization algorithm for allocation of wind turbines in distribution systems. Appl. Soft Comput. 105, 107278 (2021).10.1016/j.asoc.2021.107278
37. Naderipour A Nowdeh SA Saftjani PB Abdul-Malek Z Mustafa MWB Kamyab H Davoudkhani IF Deterministic and probabilistic multi-objective placement and sizing of wind renewable energy sources using improved spotted hyena optimizer J. Clean. Product. 2021 286 124941 10.1016/j.jclepro.2020.124941
Naderipour, A. et al. Deterministic and probabilistic multi-objective placement and sizing of wind renewable energy sources using improved spotted hyena optimizer. J. Clean. Product. 286, 124941 (2021).10.1016/j.jclepro.2020.124941
38. Jafar-Nowdeh A Babanezhad M Arabi-Nowdeh S Naderipour A Kamyab H Abdul-Malek Z Ramachandaramurthy VK Meta-heuristic matrix moth–flame algorithm for optimal reconfiguration of distribution networks and placement of solar and wind renewable sources considering reliability Environ. Technol. Innov. 2020 20 101118 10.1016/j.eti.2020.101118
Jafar-Nowdeh, A. et al. Meta-heuristic matrix moth–flame algorithm for optimal reconfiguration of distribution networks and placement of solar and wind renewable sources considering reliability. Environ. Technol. Innov. 20, 101118 (2020).10.1016/j.eti.2020.101118
39. Ghaffari A Askarzadeh A Fadaeinedjad R A power quality-based framework for optimal siting and sizing of wind turbines by search space reduction via power losses and flicker emission sensitivity indices IET Gener. Trans. Distrib. 2023 17 17 3801 3813 10.1049/gtd2.12932
Ghaffari, A., Askarzadeh, A. & Fadaeinedjad, R. A power quality-based framework for optimal siting and sizing of wind turbines by search space reduction via power losses and flicker emission sensitivity indices. IET Gener. Trans. Distrib. 17(17), 3801–3813 (2023).10.1049/gtd2.12932
40. Ch Y Goswami SK Chatterjee D Effect of network reconfiguration on power quality of distribution system Int. J. Electri. Power Energy Syst. 2016 83 87 95 10.1016/j.ijepes.2016.03.043
Ch, Y., Goswami, S. K. & Chatterjee, D. Effect of network reconfiguration on power quality of distribution system. Int. J. Electri. Power Energy Syst. 83, 87–95 (2016).10.1016/j.ijepes.2016.03.043
41. Dashtaki AA Hakimi SM Hasankhani A Derakhshani G Abdi B Optimal management algorithm of microgrid connected to the distribution network considering renewable energy system uncertainties Int. J. Electri. Power Energy Syst. 2023 145 108633 10.1016/j.ijepes.2022.108633
Dashtaki, A. A., Hakimi, S. M., Hasankhani, A., Derakhshani, G. & Abdi, B. Optimal management algorithm of microgrid connected to the distribution network considering renewable energy system uncertainties. Int. J. Electri. Power Energy Syst. 145, 108633 (2023).10.1016/j.ijepes.2022.108633
42. Aghdam FH Javadi MS Catalão JP Optimal stochastic operation of technical virtual power plants in reconfigurable distribution networks considering contingencies Int. J. Electri. Power Energy Syst. 2023 147 108799 10.1016/j.ijepes.2022.108799
Aghdam, F. H., Javadi, M. S. & Catalão, J. P. Optimal stochastic operation of technical virtual power plants in reconfigurable distribution networks considering contingencies. Int. J. Electri. Power Energy Syst. 147, 108799 (2023).10.1016/j.ijepes.2022.108799
43. Khasanov M Kamel S Halim Houssein E Rahmann C Hashim FA Optimal allocation strategy of photovoltaic-and wind turbine-based distributed generation units in radial distribution networks considering uncertainty Neural Comput. Appl. 2023 35 3 2883 2908 10.1007/s00521-022-07715-2
Khasanov, M., Kamel, S., Halim Houssein, E., 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. 35(3), 2883–2908 (2023).10.1007/s00521-022-07715-2
44. Wu R Liu S Multi-objective optimization for distribution network reconfiguration with reactive power optimization of new energy and EVs IEEE Access 2023 11 10664 10674 10.1109/ACCESS.2023.3241228
Wu, R. & Liu, S. Multi-objective optimization for distribution network reconfiguration with reactive power optimization of new energy and EVs. IEEE Access 11, 10664–10674 (2023).10.1109/ACCESS.2023.3241228
45. Dinh, B. H., Nguyen, T. T., & Nguyen, T. T. Optimal placement of photovoltaic systems and wind turbines in distribution systems by using northern goshawk optimization algorithm. In Advances in Artificial Intelligence for Renewable Energy Systems and Energy Autonomy 221–241 (Springer, 2023).
46. Aien M Hajebrahimi A Fotuhi-Firuzabad M A comprehensive review on uncertainty modeling techniques in power system studies Renew. Sustain. Energy Rev. 2016 57 1077 1089 10.1016/j.rser.2015.12.070
Aien, M., Hajebrahimi, A. & Fotuhi-Firuzabad, M. A comprehensive review on uncertainty modeling techniques in power system studies. Renew. Sustain. Energy Rev. 57, 1077–1089 (2016).10.1016/j.rser.2015.12.070
47. Abud TP Augusto AA Fortes MZ Maciel RS Borba BS State of the art monte Carlo method applied to power system analysis with distributed generation Energies 2022 16 1 394 10.3390/en16010394
Abud, T. P., Augusto, A. A., Fortes, M. Z., Maciel, R. S. & Borba, B. S. State of the art monte Carlo method applied to power system analysis with distributed generation. Energies 16(1), 394 (2022).10.3390/en16010394
48. Maulik A Probabilistic power management of a grid-connected microgrid considering electric vehicles, demand response, smart transformers, and soft open points Sustain. Energy Grids Netw. 2022 30 100636 10.1016/j.segan.2022.100636
Maulik, A. Probabilistic power management of a grid-connected microgrid considering electric vehicles, demand response, smart transformers, and soft open points. Sustain. Energy Grids Netw. 30, 100636 (2022).10.1016/j.segan.2022.100636
49. Che Y Wang X Lv X Hu Y Probabilistic load flow using improved three point estimate method Int. J. Electri. Power Energy Syst. 2020 117 105618 10.1016/j.ijepes.2019.105618
Che, Y., Wang, X., Lv, X. & Hu, Y. Probabilistic load flow using improved three point estimate method. Int. J. Electri. Power Energy Syst. 117, 105618 (2020).10.1016/j.ijepes.2019.105618
50. Zhao S Zhang T Ma S Chen M Dandelion Optimizer: A nature-inspired metaheuristic algorithm for engineering applications Eng. Appl. Artif. Intell. 2022 114 105075 10.1016/j.engappai.2022.105075
Zhao, S., Zhang, T., Ma, S. & Chen, M. Dandelion Optimizer: A nature-inspired metaheuristic algorithm for engineering applications. Eng. Appl. Artif. Intell. 114, 105075 (2022).10.1016/j.engappai.2022.105075
51. Cao D Xu Y Yang Z Dong H Li X An enhanced whale optimization algorithm with improved dynamic opposite learning and adaptive inertia weight strategy Complex Intell. Syst. 2023 9 1 767 795 10.1007/s40747-022-00827-1
Cao, D., Xu, Y., Yang, Z., Dong, H. & Li, X. An enhanced whale optimization algorithm with improved dynamic opposite learning and adaptive inertia weight strategy. Complex Intell. Syst. 9(1), 767–795 (2023).10.1007/s40747-022-00827-1
52. Mirjalili S Mirjalili SM Lewis A Grey wolf optimizer Adv. Eng. Softw. 2014 69 46 61 10.1016/j.advengsoft.2013.12.007
Mirjalili, S., Mirjalili, S. M. & Lewis, A. Grey wolf optimizer. Adv. Eng. Softw. 69, 46–61 (2014).10.1016/j.advengsoft.2013.12.007
53. Abualigah L Diabat A Mirjalili S Abd Elaziz M Gandomi AH The arithmetic optimization algorithm Comput. Methods Appl. Mech. Eng. 2021 376 113609 10.1016/j.cma.2020.113609
Abualigah, L., Diabat, A., Mirjalili, S., Abd Elaziz, M. & Gandomi, A. H. The arithmetic optimization algorithm. Comput. Methods Appl. Mech. Eng. 376, 113609 (2021).10.1016/j.cma.2020.113609
54. Haghifam MR Falaghi H Malik OP Risk-based distributed generation placement IET Gener. Transm. Distrib. 2008 2 2 252 260 10.1049/iet-gtd:20070046
Haghifam, M. R., Falaghi, H. & Malik, O. P. Risk-based distributed generation placement. IET Gener. Transm. Distrib. 2(2), 252–260 (2008).10.1049/iet-gtd:20070046
55. Swarnkar A Gupta N Niazi KR Adapted ant colony optimization for efficient reconfiguration of balanced and unbalanced distribution systems for loss minimization Swarm Evolut. Comput. 2011 1 3 129 137 10.1016/j.swevo.2011.05.004
Swarnkar, A., Gupta, N. & Niazi, K. R. Adapted ant colony optimization for efficient reconfiguration of balanced and unbalanced distribution systems for loss minimization. Swarm Evolut. Comput. 1(3), 129–137 (2011).10.1016/j.swevo.2011.05.004
56. Subrahmanyam JBV Radhakrishna C A simple method for feeder reconfiguration of balanced and unbalanced distribution systems for loss minimization Electric Power Compon. Syst. 2009 38 1 72 84 10.1080/15325000903273270
Subrahmanyam, J. B. V. & Radhakrishna, C. A simple method for feeder reconfiguration of balanced and unbalanced distribution systems for loss minimization. Electric Power Compon. Syst. 38(1), 72–84 (2009).10.1080/15325000903273270
57. Samal, P., Mohanty, S., & Ganguly, S. Planning of distributed generation and capacitor in an unbalanced radial distribution system using cuckoo search algorithm. In 2016 19th International Conference on Electrical Machines and Systems (ICEMS) 1–5 (IEEE, 2016).
58. Gerez C Coelho Marques Costa E Sguarezi Filho AJ Distribution network reconfiguration considering voltage and current unbalance indexes and variable demand solved through a selective bio-inspired metaheuristic Energies 2022 15 5 1686 10.3390/en15051686
Gerez, C., Coelho Marques Costa, E. & Sguarezi Filho, A. J. Distribution network reconfiguration considering voltage and current unbalance indexes and variable demand solved through a selective bio-inspired metaheuristic. Energies 15(5), 1686 (2022).10.3390/en15051686
59. Montoya OD Arias-Londoño A Grisales-Noreña LF Barrios JÁ Chamorro HR Optimal demand reconfiguration in three-phase distribution grids using an MI-convex model Symmetry 2021 13 7 1124 10.3390/sym13071124
Montoya, O. D., Arias-Londoño, A., Grisales-Noreña, L. F., Barrios, J. Á. & Chamorro, H. R. Optimal demand reconfiguration in three-phase distribution grids using an MI-convex model. Symmetry 13(7), 1124 (2021).10.3390/sym13071124
60. Naderipour A Abdullah A Marzbali MH Nowdeh SA An improved corona-virus herd immunity optimizer algorithm for network reconfiguration based on fuzzy multi-criteria approach Expert Syst. Appl. 2022 187 115914 10.1016/j.eswa.2021.115914 34566274
Naderipour, A., Abdullah, A., Marzbali, M. H. & Nowdeh, S. A. An improved corona-virus herd immunity optimizer algorithm for network reconfiguration based on fuzzy multi-criteria approach. Expert Syst. Appl. 187, 115914 (2022).34566274 10.1016/j.eswa.2021.115914
61. Alanazi A Alanazi M Artificial electric field algorithm-pattern search for many-criteria networks reconfiguration considering power quality and energy not supplied Energies 2022 15 14 5269 10.3390/en15145269
Alanazi, A. & Alanazi, M. Artificial electric field algorithm-pattern search for many-criteria networks reconfiguration considering power quality and energy not supplied. Energies 15(14), 5269 (2022).10.3390/en15145269
62. Zidan, A., & El-Saadany, E. F. Service restoration in balanced and unbalanced distribution systems with high DG penetration. In 2011 IEEE power and energy society general meeting 1–8 (2011, IEEE).
