Simultaneous Distribution Network Reconﬁguration and Optimal Allocation of Renewable-Based Distributed Generators and Shunt Capacitors under Uncertain Conditions

: Smart grid technology has received ample attention in past years to develop the traditional power distribution network and to enable the integration of distributed generation units (DGs) to satisfy increasing demand loads and to improve network performance. In addition to DGs, integration of shunt capacitors (SCs) along with network reconﬁguration can also play an important role in improving network performance. Besides, network reconﬁguration can help to increase the distributed generation hosting capacity of the network. Some of the research in the literature have presented and discussed the problem of optimal integration of renewable DGs and SCs along with optimal network reconﬁguration, while the network load variability and/or the intermittent nature of renewable DGs are neglected. For the work presented in this paper, the SHADE optimization algorithm along with the SOE reconﬁguration method have been employed for solving the aforementioned optimization problem with consideration of uncertainty related to both the network load and the output power of the renewable DGs. Maximizing the hosting capacity (HC) of the DGs and reducing network power losses in addition to improving the voltage proﬁle have been considered as optimization objectives. Five different case studies have been conducted considering 33-bus and 59-bus distribution networks. The obtained results validate the effectiveness and the superiority of the employed techniques for maximizing the HC up to 17% and reducing power losses up to 95%. Besides, the results also depict the effect of SC integration and the consideration of uncertainties on achieving the optimization objectives with realistic modeling of the optimization problem.


Introduction
Electrical power has been provided efficiently via the electrical power systems for many long years, from large and centralized electricity generating stations to different consumers in homes and in industries.The traditional pattern of power generation and transmission has been changed by the constantly increasing global population and electrical demand loads which have led to saturation of the electric power network, negative environmental impacts from large centralized power generating plants such as global warming and intensive carbon emissions, changes in consumer behavior, and newly emerging techniques to get energy from renewable sources.The traditional electric power network has been developed to be a more intelligent, so-called smart grid technology (SGT) [1].
A smart grid is simply an intelligent electric power network with two ways of power flow and communications, which provide electric systems with better robustness, security, Energies 2022, 15, 2299 2 of 27 flexibility, reliability, and efficiency.The ability to incorporate one or more small power generators, called distributed generators (DGs), into a distribution network's feeders is considered one of the main benefits of smart grid technology.A distributed generator (DG) is simply known as an electric power source that is added to distribution network buses or feeders and it is allocated close to different consumers or points of consumption in order to improve network performance and the quality of the delivered power, enhance system reliability, reduce network power losses, and retard the need for new additional costs for the addition of new distribution feeders.There is only one unidirectional way for power flow from the power generating stations to distribution systems in conventional electric power systems.Recently, with the integration of DGs in distribution networks, generated power can flow from distribution feeders to substations, providing two ways of power flow (i.e., bidirectional power flow) [2].Distributed generators can be employed as prime, standby electric power sources, or even reactive power sources, incorporated in residential, industrial, or commercial areas.From a technological aspect, DGs might be renewable power sources, such as photovoltaic systems, wind turbines, and biomass, or non-renewable conventional sources, such as fuel cells, natural gas, and micro-turbines.Renewable energy sources (RESs) produce output power that is uncertain, very hard to predict with high accuracy, and intermittent as it depends on environmental factors, such as wind speed and solar irradiance.In the past few years, and it is predicted to increase in the near future, RESs are spreading throughout the world, and they have a major contribution to existing power systems.As a result, RESs, with their intermittent nature, are replacing conventional non-renewable sources, which leads to decreasing the power flow predictability in the network.The location and the size of DGs in a distribution network might lead to positive or negative impacts on network performance.Hence, they have to be incorporated into the network after comprehensive planning and study to determine their optimal locations and sizes in order to avert any negative impact, such as overvoltage, exceeding the acceptable loading limits of network feeders, harmonic overloading, etc.
In addition to DGs, installing shunt capacitors (SCs) in distribution systems can improve the power quality and the power factor, enhance the voltage profile, and reduce power losses by injecting reactive power.A shunt capacitor draws a leading current that modifies the characteristics of an inductive load by opposing some or all of the lagging components of the inductive load current at the installation point.Similar to DGs, adding shunt capacitors to distribution systems will be useful only when optimal locations and sizes are applied, and any inappropriate location or size may lead to negative impacts, such as increasing the power losses and/or decreasing the voltage value below the acceptable limit.
In the last few years, several different techniques have been presented in many studies for determining the optimal placement and sizing of DGs or SCs separately in distribution systems.Other researchers have proposed other techniques for simultaneous optimal integration of both DGs and SCs for compensating the real and the reactive power.In [3], different types of DGs, their various technologies, the economic and the technical advantages gained from incorporating DGs in electric power distribution systems and the imposed constraints were presented.The authors in [4,5] presented a comprehensive general view of distributed generation and the evolutions in DG technologies in addition to a review of the different optimization techniques employed for solving the problem of the optimal allocation and sizing of DGs.A novel multi-objective optimization strategy based on the PSO algorithm was presented in [6] for solving the problem of distribution system planning with optimal DG integration.In [7], the authors employed a simple heuristic optimization approach, which was based on loss sensitivity analysis for solving the problem of optimal DG allocation in distribution networks with the objective of minimizing network power losses.The problem of optimal allocation of DGs in distribution networks was handled in [8] using the GA optimization algorithm and considering economic and technical objectives.In [9], an advanced variant of the PSO algorithm was employed to solve the problem of optimal DG allocation, considering the objectives of maximizing annual savings and enhancing the bus voltage and feeder current profiles.The Kalman filter algorithm was employed in [10] for finding the optimal allocation of DGs in the distribution system.The authors in [11] employed a strategy, which was based on a voltage stability index and used the combined sensitivity factor analogy for solving the problem of the optimal allocation and sizing of DGs with multiple types in a 48-node distribution network, considering the objectives of reducing network power losses and enhancing the bus voltage profile.The genetic algorithm was been used in [12] for solving the optimal allocation problem of DGs in the distribution system for minimizing power losses.In [13], the problem of optimal DG sizing was handled and solved by achieving the objectives of minimizing the costs of DGs and maximizing network loadability.The authors in [14] employed the GA and the optimal power flow solution for solving the problem of simultaneous optimal allocation of both shunt capacitors and voltage regulators with achieving optimal investment costs.The fuzzy GA was employed in [15] for solving the problem of optimal allocation of shunt capacitors in distribution systems with the objective of power factor improvement.In [16], the mixedinteger nonlinear programming technique (MINLP) was employed to solve the problem of optimal placement and sizing of shunt capacitors in distribution systems.An optimization approach was presented and employed in [17] for solving the problem of optimal placement and sizing of shunt capacitors and active power conditioners in distribution systems for power quality improvement.To improve the performance of the power system, the grey wolf optimization technique was employed in [18] for determining the optimal sizing of reactive power.Optimal simultaneous placement of both DGs and shunt capacitors in the distribution system was presented in [19][20][21][22] for reducing or minimizing system power losses.In [23], a pareto multi-objective PSO technique was employed for determining the optimal placement and sizing of DGs and shunt capacitors in a distribution network with consideration of network load-related uncertainty, which was modeled using fuzzy logic.Enhanced meta-heuristic optimization methods were employed in [24] for solving the problem of optimal simultaneous allocation of DGs and shunt capacitors in distribution systems, considering demand load variability and with the objectives of minimizing annual power losses and bus voltage profile improvement.The authors in [25] employed the binary collective animal behavior optimization technique for optimal simultaneous placement and sizing of DGs and shunt capacitors in distribution systems for achieving the objectives of minimizing power losses and improving the voltage profile.In [26], the biogeography-based optimization technique was employed to solve the problem of DGs and shunt capacitors optimal allocation with consideration of the optimization objectives of minimizing the total harmonic distortion, minimizing power losses, and improving the voltage profile.The intersect mutation differential evolution (IMDE) optimization technique was presented in [27] for simultaneous optimal allocation and sizing of DGs and shunt capacitors, with the objective of minimizing the power losses and subject to satisfying the bus voltage and feeder current constraints.In [28], GA optimization strategy was employed for optimal simultaneous sizing of both DGs and shunt capacitors, with the objective of minimizing investment and maintenance costs.The authors in [29] presented and employed a sorting-based, non-dominant, multi-objective PSO optimization algorithm, with consideration of fuzzy decision principles for simultaneous optimal placement of renewable DGs and shunt capacitors.Different optimization techniques were employed in [30][31][32] for optimal allocation and sizing of DGs and shunt capacitors in distribution systems, with the consideration of minimizing power losses as the optimization objective.Optimal allocation of shunt capacitors in a distribution system integrated with wind energy-based DGs was performed in [33] by using a stochastic technique, which was based on a point estimation method.A hybrid optimization approach of weight improved PSO and the gravitational search algorithm (GSA) was employed in [34] for optimal DG and shunt capacitor allocation and sizing.In [35], simultaneous optimal placement of DGs and shunt capacitors was performed with consideration of the intermittency and uncertainty associated with DGs.
In addition to installing shunt capacitors (SCs) and distributed generation units (DGs), distribution network power loss reduction can also be achieved by reconfiguring the existing network topology.A distribution network includes two types of switches, namely sectionalizing and tie-switches.Although both types of switches are principally the same, tie switches are those switches in the network that remain open (normally opened) to preserve the radial structure of the network, while sectionalizing switches are those switches that remain closed (normally closed).To provide appropriate voltage control and protection coordination, radial topology is mostly recommended, and it is preferred for distribution systems.Each node in the radial distribution network has only one way to the substation node, and no loops exist in the network.Distribution network reconfiguration (DNR) can be simply defined as the process of changing the existing network topology or configuration via altering the opened/closed status of tie and sectionalizing switches, such that the radial structure of the network is preserved, all of the network loads are connected, all objectives are achieved, and all considered constraints are satisfied.The main and the conventional objective of network reconfiguration is power loss reduction and alleviating the overload of network feeders.However, distribution network reconfiguration has been employed in the literature for achieving other objectives, such as improving power quality and network loadability, enhancing the voltage profile, reducing operational costs and emissions, and providing more reliable and economic network operation [36][37][38][39][40][41][42][43][44][45][46][47].Moreover, DNR has been performed for unbalanced distribution networks [48,49], managing distribution system outages [50], reducing annual energy losses [51], estimation of active distribution system states [52], an active distribution network with electric vehicles [53], an active distribution network with DGs [54], etc.On the other hand, other research related to DNR has provided more efficient and effective computations via parallel computation on graphics processing units [55], an extended fast decoupled Newton-Raphson power flow methodology [56], an accelerating decomposition technique [57,58], and fast calculations for involved frequent power flow solutions [54].In addition, distribution network reconfiguration has also been employed in the literature as an approach for maximizing the distributed generation hosting capacity (i.e., the DG penetration that can be accommodated) of the network [59,60].
The limited electric power generation and variability associated with distribution system demand loads make the process of operating and controlling a distribution system a difficult and a complicated task, especially with high load densities.Minimized system power losses cannot be achieved considering the same network configuration for all levels of variable system loads, consequently, frequent reconfiguration of the network is required; and, for satisfying the required demand loads, DGs and SCs are incorporated in the network and this leads to a decrease in power losses, an enhanced voltage profile, alleviation of the feeders' overload, and improved reliability and efficiency of the power supply.For solving the problem of distribution network reconfiguration and DG and SC optimal integration, the variability associated with the network demand load in addition to the intermittent output power from renewable DGs must be considered.Allocation and sizing of integrated DGs and SCs in the distribution network, which is determined using voltage or loss indices, may not be the optimal solution for all network configurations; hence, distribution network reconfiguration and shunt capacitor and distributed generator integration need to be performed simultaneously to provide an optimal solution for minimizing power losses and enhancing network operation and the voltage profile.
Several different types of research in the literature have focused on the static singlehour DNR for providing an appropriate and fixed network configuration for a certain time duration-which can be an hour, a year, a season, or as required and applicable-without considering either the variations of DG unit output power with time or the variable network demand load where both of the generated power and demand loads are set to their peak or average values [61][62][63].Accordingly, the static single-hour DNR will not lead to an optimal solution for the DNR problem for distribution networks integrated with renewable DG units with intermittent output power, such as PV units or wind turbines.In light of the above, other research in the literature investigated the dynamic multi-hour DNR to account for the intermittent nature of renewable DGs and the variability associated with the network demand loads [64,65].Using dynamic DNR, the optimal configurations for multiple networks are provided for different hours or durations based on the variations of loads and the power generation of DGs with achieving the optimization objective, mainly the minimization of network power losses and satisfying the considered constraints for all provided configurations.In conclusion, the static DNR is more economically beneficial than dynamic reconfiguration, where a fixed configuration is provided and the cost of more switching actions is avoided.However, optimal and more accurate solutions are provided using the dynamic approach for the networks integrated with renewable DGs with intermittent output power.
Considering variable demand loads, distribution network reconfiguration problems were presented in [66,67].Seasonal distribution network reconfiguration was presented in [68], and the related switches to be opened or closed for changing the network configuration were determined.Although the variability of demand loads has been considered in [66][67][68] for distribution network reconfiguration, the output power of DG units has been considered a constant value disregarding the associated variability.In [69,70], optimal network configurations for the day ahead have been provided by applying multi-hour dynamic reconfiguration based on foreseen values for the power generation and demand loads.Multi-hour DNR has been performed in [71] with consideration of variable demand loads and renewable DGs with the objective of power loss minimization.Distribution network reconfiguration was performed in [72] to provide the optimal seasonal network configuration, considering variable demand load and renewable DG output power.In [73][74][75], the effect of considering the load variability and the uncertainty associated with renewable DG output power on the reconfiguration process of the distribution network was studied.Distribution network reconfiguration with the integration of distributed generators was performed in [76] as a mixed-integer linear programming problem.Optimal DG allocation and distribution network reconfiguration using a heuristic algorithm based on sensitivity indices for power loss minimization was presented in [77].The problem of distribution network reconfiguration along with distributed generation allocation was modeled in [78] as a multi-objective nonlinear mixed-integer programming problem, with a consideration of the demand load variability and the uncertainties associated with the output power of renewable DGs.Simultaneous optimal network reconfiguration and DG allocation in distribution networks was presented in [79] using a meta-heuristic harmony search algorithm (HAS), while the same problem was presented in [80], using an adaptive genetic algorithm and a graph theory.In [81], distribution network reconfiguration with distributed generation integration was presented for network reliability improvement and power loss minimization.A simulated annealing technique was used in [82] for distribution network feeder reconfiguration and capacitor setting to reduce power loss and to improve the voltage profile.In [83], a mixed-integer non-linear programming technique was used for optimal distribution network reconfiguration and shunt capacitors allocation with the objective of power loss minimization, while the same problem was presented in [84], using the ant colony search optimization technique.Simultaneous optimal network reconfiguration and shunt capacitor allocation in a distribution network was presented in the literature using a genetic algorithm [85], a modified particle swarm optimization algorithm [86], a hybrid optimization technique using a minimal nodal voltage approach and a genetic algorithm [87], the Chu-Beasley based genetic algorithm [88], a fuzzy-genetic algorithm [89], ordinal optimization [90], a minimum spanning tree algorithm [91], and a selective particle swarm optimization [92] and harmony search algorithm [93].
In view of the above and according to the authors' best knowledge, simultaneous distribution network reconfiguration with optimal allocation and sizing of renewable distributed generators and shunt capacitors for minimizing network power loss, maximizing the DG hosting capacity, and improving the voltage profile, as optimization objectives, with consideration of the network demand load variability and the uncertain parameters related to the considered renewable DGs, haven't been presented so far.Hence, for the work presented in this paper, the optimal allocation and sizing of renewable DGs (PV units and wind turbines) and shunt capacitors in the distribution network are implemented simultaneously with the network reconfiguration, such that network power losses are minimized, distributed generation hosting capacity of the network is maximized, and the voltage profile is enhanced, while keeping the considered constraints within their desired and acceptable limits.The 33-bus distribution network and 59-bus real distribution network in Cairo, Egypt, are considered for this study, taking into consideration the variability of demand loads in addition to the uncertainty associated with the output power from the considered renewable DGs in the optimization problem using the scenario-based technique.As indicated before, the output power from renewable DGs depends on environmental factors (such as wind speed and solar irradiance), which have an intermittent nature, and they are very difficult to predict.The uncertainty in the demand load is modeled using the normal probability density function (PDF), while the stochastic wind speed is modeled using Weibull PDF, and the stochastic solar irradiance is modeled using lognormal PDF.
For each of these uncertain parameters, one thousand scenarios are generated by running a Monte Carlo simulation, and then a scenario reduction technique called the backwards scenario reduction technique is applied to select a specific number from the generated scenarios.The optimization process shall be performed by using the success history-based adaptive differential evolution algorithm (SHADE) [94], which is an advanced variant of the differential evolution algorithm (DE).In the literature, there are three different categories for the reconfiguration techniques, namely, evolutionary algorithms (EAs) [61,62,95,96], mathematical programming (MP) [63,65,76,97,98], and heuristics-based techniques [99,100].For the work presented in this paper, the heuristics-based reconfiguration method presented in [101], called the switch opening and exchange method (SOE), will be adopted for distribution network reconfiguration.This method will be illustrated in detail within Section 4, using a three-feeder small example system.
The rest of the paper will be organized as follows: the adopted model for the variability and the uncertainty associated with demand load, wind speed, and solar irradiance is described in Section 2. The problem formulation with the definition of variables, optimization objectives, and considered constraints is given in Section 3. The employed reconfiguration method (SOE method) is illustrated in Section 4.An overview of the success history based adaptive differential evolution optimization algorithm (SHADE) is outlined in Section 5.The simulation results analysis and the discussion are given in Section 6. Final conclusions and recommendations for future work are listed in Section 7.

Modeling of Load Demand Variability, Stochastic Wind Speed, and Solar Irradiance
This section discusses and explains modeling of the uncertainty associated with demand load, wind speed, and solar irradiance via scenario generation and reduction techniques.
The variability of load demand is modeled using the normal probability density function (PDF) with an associated mean (µ d ) and standard deviation (σ d ).The probability density of uncertain demand load (P d ) is expressed as follows: Based on this normal distribution of demand load, one thousand scenarios are generated using a Monte Carlo simulation.
The uncertainty associated with wind speed, which affects the output power from wind turbines, is modelled using the Weibull PDF, which is usually used for representing wind speed distribution.The probability density of uncertain wind speed (ν w ) is expressed as follows: where α is the scale parameter and β is the shape parameter for the Weibull PDF.Based on this distribution of wind speed, one thousand scenarios are generated using a Monte Carlo simulation.
The uncertainty associated with solar irradiance, which affects the output power from PV units, is modelled using the lognormal PDF, with mean (µ s ) and standard deviation (σ s ), which is usually used for representing solar irradiance distribution.The probability density of uncertain solar irradiance (G s ) is expressed as follows: Similar to demand load and wind speed, based on the lognormal distribution, one thousand scenarios are also generated using a Monte Carlo simulation.These generated scenarios for solar irradiance consider that sunlight is not available during the whole day and usually it is available only half of the day.Hence, the solar irradiance and, accordingly, the output power from PV units is zero during this period of the day.
The scenario generation technique followed in the work presented in this paper shall be the same technique presented in [102], where the generated one thousand scenarios, as discussed above, for demand load, wind speed and solar irradiance are combined, which results in one thousand scenarios with each scenario including three values (demand load, wind speed, and solar irradiance); hence, the i th scenario (S i ) can be represented as As it is impractical to handle so many scenarios, an appropriate scenario reduction technique is necessary.In this paper, the scenario reduction technique presented in [102]called the backward reduction technique-shall be applied in order to reduce the generated one thousand scenarios to thirty scenarios only, which is practically a good number of scenarios for optimization problems.In this technique, the procedure of scenario reduction depends on the distance between each pair of scenarios, where one scenario (the scenario with lower probability) is removed from the pair of scenarios having minimum distance and this procedure is repeated until the required number of scenarios is reached.The thirty scenarios resulting from the applied scenario reduction technique with their respective probabilities are listed below in Table 1.As mentioned before, each scenario represents values for demand load (as a percentage of total load), wind speed, and solar irradiance.Based on the values of wind speed and solar irradiance, wind and solar power (output power from wind turbines and PV units) are calculated as discussed in the following paragraphs.
For a wind turbine with rated power P WT rated , the output power generated as a function of wind speed (ν w ) can be calculated as follows: where v cut-in , v cut-out and v w,rated are the cut-in, cut-out, and rated wind speed of the turbine, respectively.
For a PV unit with rated power P PV rated , the output power generated as a function of solar irradiance (G s ) can be calculated as follows: where G std and R c are the standard environment solar irradiance and specified irradiance value, respectively.

Objective Function
In this section, a multi-objective optimization problem is formulated for maximizing the distributed generation probabilistic hosting capacity (PHC) of the distribution network, minimizing network power losses (i.e., maximizing the reduced power losses (RP loss )), and improving the voltage profile, while satisfying the operational constraints and limits for different scenarios of load demand, wind speed, and solar irradiance.The objective function for the optimization problem can be formulated as follows: Maximize RP loss (%) = P 0 loss − P loss P 0 loss × 100 (7) where HC s is the hosting capacity of DGs for the s th scenario, p s is the probability related to a scenario (s), N s is the total number of scenarios, P 0 loss is the total network power loss at normal loading conditions (Base Configuration), and P loss is the probabilistic total network power loss for all of the studied scenarios.The HC s and the P loss can be calculated as follows: Energies 2022, 15, 2299 9 of 27 where P WT n WT ,s is the output power of the wind turbine allocated at the bus or the node (n WT ) for scenario (s), P PV n PV ,s is the output power of the PV unit allocated at the bus or the node (n PV ) for scenario (s), P L n is the load power at the n th node, N nodes is the total number of network nodes, I b,s is the magnitude of the current flowing in the b th branch for scenario (s), r b is the b th branch resistance, and N line is the total number of network lines or branches.
For improving the voltage profile, minimizing the overall aggregated voltage deviation index (AVDI ov ) for all different considered scenarios is also set as an objective for the optimization problem.For each scenario, the aggregated voltage deviation index (AVDI s ) is calculated as the sum of voltage deviations of load buses from unity as follows: where V lb,s is the voltage at the load bus (lb) for scenario (s), and N lb is the total number of load buses.Hence, the objective for minimizing the overall aggregated voltage deviation index (AVDI ov ), for all considered scenarios, can be formulated as follows to ensure that voltages at the network load buses are close to unity:

Optimization Constraints
The operational and planning constraints applied for the considered optimization problem shall be as follows: HC s ≤ 100 ( 14) where P WT n WT ,rated and P PV n PV ,rated are the rated output power or capacity of the wind turbine and the PV unit allocated at nodes n WT and n PV , respectively.The P WT max and the P PV max are the maximum allowed rated power or capacity for the wind turbine and the PV unit, respectively.The I rated is the maximum allowed current of a network branch, |V n,s | is the magnitude of voltage at the n th node, V min and V max are the minimum and the maximum allowed node or bus voltage, respectively.

Switch Opening and Exchange Method
The switch opening and exchange method (SOE) employed in this paper for performing the distribution network reconfiguration is a novel heuristic method which was first proposed in [101] for efficient solving of large-scale multi-hour stochastic distribution network reconfiguration and it has proven its superiority to other heuristic algorithms, EAs, and MPs in terms of accuracy and/or solution speed.This method includes three stages and it adopts both of the sequential switch opening and branch exchanging techniques for solving the reconfiguration problem, where the first stage applies the sequential switch opening, while the second and the third apply the branch exchanging.These stages shall be illustrated in detail within the following paragraphs, using the three-feeder distribution system shown below in Figure 1 as an example system.and MPs in terms of accuracy and/or solution speed.This method includes three stages and it adopts both of the sequential switch opening and branch exchanging techniques for solving the reconfiguration problem, where the first stage applies the sequential switch opening, while the second and the third apply the branch exchanging.These stages shall be illustrated in detail within the following paragraphs, using the three-feeder distribution system shown below in Figure 1 as an example system.

Stage 1: Sequential Switch Opening
The status of all switches in this stage is initially set to be closed which forms loops in the distribution network.In each iteration, the status of only one switch (selected from the sectionalizing closed switches within the loops) is changed to be opened (i.e., changed to be a tie-switch) and the loop, where that switch is located, is opened.Hence, the number of iterations in this stage shall be equal to the number of loops in the network and the above procedure is repeated in each iteration till all loops are opened and the network

Stage 1: Sequential Switch Opening
The status of all switches in this stage is initially set to be closed which forms loops in the distribution network.In each iteration, the status of only one switch (selected from the sectionalizing closed switches within the loops) is changed to be opened (i.e., changed to be a tie-switch) and the loop, where that switch is located, is opened.Hence, the number of iterations in this stage shall be equal to the number of loops in the network and the above procedure is repeated in each iteration till all loops are opened and the network becomes radial.The criteria for selecting the sectionalizing switch to be opened or changed to a tie-switch in each iteration shall be as follows: (a) Based on the network configuration, close all of the network switches and determine the sectionalizing switches or branches within the formed loops.In the above example system, closing all switches forms three loops that include all of the system branches except for branch 20 (connecting nodes 9 and 12), which is not included within any of the formed loops.Hence, three iterations shall be required in thin stage for this system to select which sectionalizing switches to be changed to tie-switches.(b) For each switch within the loops: b.1.
Open the switch.b.2.
Check the network connectivity (i.e., ensure that there is a path for each node to the substation node).b.3.
If the network is not connected, assign a large value (e.g., 10 6 ) for the network power loss.b.4.
Else if the network is connected, provide the power loss value for the network after opening the switch.b.5.
Close the opened switch.b.6.
The above-described steps (from b.1 to b.5) are repeated for all other remaining sectionalizing switches within the loops.b. 7.
Permanently open the switch that achieves the minimum network power loss value among all other switches.This will consequently open the loop where this switch is located.b.8.
The above steps (from b.1 to b.7) are repeated for the sectionalizing switches within the remaining loops until the network becomes radial (i.e., all loops are opened).The obtained radial network configuration is called the "initial configuration".b.9.
Provide the power loss value associated with the obtained radial network configuration and record this value as "loss-0".

Stage 2 & Stage 3: Branch Exchanging
In these stages, other different radial configurations are generated from the initial configuration obtained from the first stage searching for a better solution.This is achieved by opening a sectionalizing switch in the initial configuration (i.e., changing it to a tieswitch) and then applying the steps of the first stage (keeping the status of that switch to be opened) which leads to a different radial network configuration.This procedure is applied to the other sectionalizing switches in the initial configuration to generate other radial configurations.Hence, for all sectionalizing switches within the initial configuration, the following procedure shall be applied: (a) Open the sectionalizing switch and apply the steps of the first stage as described above, but keep the status of this switch at permanently opened, which leads to a different radial network configuration.(b) Provide and record the power loss value associated with the obtained radial configuration.(c) Close the opened switch.(d) The above steps (from a to c) are repeated for the other remaining sectionalizing switches in the initial configuration.In order to avoid the computational complexity, the following sets of sectionalizing switches (in the initial configuration) are excluded from the above procedure: Set-1.This set includes any sectionalizing switch in the initial configuration whose shortest path to the substation node includes n 1 or less sectionalizing switches.
If n 1 = 2, then for the above example system this set shall include the sectionalizing switches 1-4, 4-5, 4-6, 2-8, 8-10, 8-9, 3-13, 13-14, and 13-15.The sectionalizing switch (x-y) refers to the switch associated with the branch connecting the nodes x and y.Set-2.This set includes any sectionalizing switch, in the initial configuration, whose shortest path to an ending node includes n 2 or less sectionalizing switches and it is upstream of the same ending node.If node A locates in the shortest path of node B to the substation node, then it can be said that node A is upstream of node B and node B is downstream of node A. The node that has no downstream nodes is referred to as the "ending node".If n 2 = 1, then for the above example system, this set shall include the sectionalizing switches 4-5, 9-11, 8-10, 13-14, 6-7, and 15-16.In this work, the values assigned for the above-mentioned parameters (n 1 & n 2 ) shall be 3 & 2, respectively.Set-3.This set includes any sectionalizing switch, in the initial configuration, which is not located in a loop when all of the network switches are closed.For the above example system, this set shall only include the sectionalizing switch 9-12.
(e) For each radial configuration obtained from the above steps (from a to d), the following procedure shall be applied: e.1.Determine the set of sectionalizing switches as defined above in "Set-2".e.2.
Exclude the sectionalizing switches as defined above in "Set-1" and "Set-3" from the sectionalizing switches obtained from e.1., and apply the following for each of the remaining switches: e.3.
Execute the found available O1C1 actions, then provide the associated power loss value after executing each action, and record this value as a new element in a vector named "loss-1".e.5.Determine the O1C1 actions that lead to a decrease in the network power loss value.These O1C1 actions are referred to as "objective decreasing O1C1 actions".e.6.
Combine and execute any two independent O1C1 actions in the set of objective decreasing O1C1 actions, then provide the associated power loss value after executing the combined actions and record this value as a new element in a vector named "loss-2".Two O1C1 actions are called "independent actions" when none of the network feeder(s) involved in one O1C1 action are the same as any of the feeder(s) involved in the other action, the network feeder that includes a branch, directly connected to the substation node, in addition to the downstream branches.For the example system in Figure 1, there are 3 independent feeders which include the branches (1-4), (2-8), and (3-13) in addition to their downstream branches.e.7.
The above steps (from e.1 to e.6) are repeated for the other radial configurations.
(f) Find the minimum power loss value recorded in "loss-0", "loss-1", and "loss-2".The corresponding network configuration shall be the optimal solution for the reconfiguration problem.
The pseudo-code of SOE method, provided below in Algorithm 1, can help for a better understanding of the above-described steps:  12 Repeat above steps (from 2 to 12) till radial configuration is obtained "Initial Configuration".13 Calculate P loss for the obtained initial configuration and store the value as "loss-0" Branch Exchanging

14
For each sectionalizing switch in the initial configuration (excluding those included in set-1, set-2 and set-3) do 15 Permanently open the switch 16 Keeping that switch opened, apply above steps (from 1 to 13) 17 Calculate P loss for the obtained radial configuration and store the obtained value 18 Close the opened switch 19 end 20 For each radial configuration obtained after applying above steps (from 14 to 16) do

21
Determine "set-2" sectionalizing switches and exclude any switch related to "set-1" or "set-3" 22 For each switch obtained from previous step do 23 Determine the available open-1 close-1 (O1C1) actions 24 Execute each O1C1 action, calculate the related Ploss and store the obtained value in a vector "loss-1" 25 Determine the objective decreasing O1C1 actions 26 Combine and execute any two independent objective decreasing actions, calculate the related Ploss and store the obtained value in a vector "loss-2" 27 end 28 end

29
The configuration corresponding to minimum power loss value stored in "loss-0", "loss-1" and "loss-2" is con-sidered as the optimal network con-figuration.

Success History Based Adaptive Differential Evolution Algorithm (SHADE)
The SHADE algorithm is an advanced variant of the differential evolution optimization algorithm, which is employed for solving the optimization problem presented in this paper.Differential Evolution (DE) is a population-based stochastic optimization algorithm, which initializes with randomly generated individuals that evolve via probabilistic operators, such as mutation and recombination.The performance of the DE algorithm depends mainly on the population size N p , mutation and/or crossover strategies, the scale factor (F) associated with mutation, and the crossover rate (CR) associated with recombination.In order to do away with the tuning of parameters, adaptive mechanisms for online adjustment of the control parameters during the evolution have been investigated and a history based parameter adaptation scheme has been proposed (SHADE).The SHADE algorithm is discussed in detail within the following subsections, which starts with initialization, and then it is followed by an iterative procedure including mutation, crossover, and selection.

Initialization
The population of candidate solutions with randomly assigned values is created during the initialization.To ensure that the randomly generated values are within the defined upper and lower bounds, the j-th component of the i-th decision vector is generated as follows: x i,j = x min,j + rand i,j [0, 1] x max,j − x min,j where rand i,j [0, 1] is a uniform random number between 0 and 1 and the superscript (1)  represents initialization.For a population size (N p ) and dimensions of a candidate solution (d), i = 1, 2, 3, . . . . . . . . ., N p and j = 1, 2, 3, . . . . . . . . ., d.

Mutation
For a given generation (t), mutation operator generates a mutant vector v (t) i corresponding to each population member x (t) i as follows: where r i 1 and r i 2 are mutually exclusive integers randomly chosen from the population range, x (t) pbest is randomly selected from best individuals of current generation, and i is a positive scale factor at t th generation.
After performing the mutation, if an element v (t) i,j violates the search boundaries [x min,j , x max,j ], it is corrected as follows:

Parameter Adaptation
For a generation (t), each included individual has a scale factor F (t) i and a crossover rate CR (t) i which are used to provide a new population of solutions.These parameters are updated as follows: r i , one of the means of successful crossover rates from previous generations stored in the memory M, and 0.1 variance.All µF and µCR are set to 0.5 during initialization and at the end of each generation, the memory is updated at positions determined by indices k.

Crossover
For a generation (t), the decision variables of mutant vector v (t) i are combined with those of target vector x (t) i and provide offspring vector u where j rand is any randomly selected number from {1, 2, . . .., d}, with d as the problem size.
The pseudo-code of SHADE algorithm can be organized as provided below in Algorithm 2: F i,G = randc i (M F,r i , 0.1); 12 p i,G = rand p min , 0.2 ; 13 Generate offspring vector u i,G 14 end

Results Analysis, Comparison and Discussion
Simulation results of the optimization problem presented and discussed in this paper, with detailed analysis and comparison against other previously reported results in the literature, are provided and discussed in this section for the considered 33-bus and 59-bus distribution networks whose basic configurations and data are as shown below in Figure 2 and Table 2, respectively.The SHADE optimization algorithm is employed for solving the optimization problem of optimal allocation and sizing of renewable DGs (PV units and wind turbines (WTs)) and shunt capacitors (SCs) in the considered networks simultaneously with the network reconfiguration using the SOE method, with consideration of the demand load variability and uncertainty associated with the renewable DG output power for addressing the optimization objectives and satisfying the operational and planning constraints.The parameters related to the renewable DGs, shunt capacitors, and SHADE algorithm considered in this study are listed below in Tables 3 and 4. The following case studies are performed and the related simulation results are obtained using MATLAB R2020a on a laptop computer with the features as Core i7 CPU, 2.60 GHz, and 16 GB of RAM.
Energies 2022, 15, x FOR PEER REVIEW 16 of 29 units and wind turbines (WTs)) and shunt capacitors (SCs) in the considered networks simultaneously with the network reconfiguration using the SOE method, with consideration of the demand load variability and uncertainty associated with the renewable DG output power for addressing the optimization objectives and satisfying the operational and planning constraints.The parameters related to the renewable DGs, shunt capacitors, and SHADE algorithm considered in this study are listed below in Tables 3 and 4. The following case studies are performed and the related simulation results are obtained using MATLAB R2020a on a laptop computer with the features as Core i7 CPU, 2.60 GHz, and 16 GB of RAM.
(a) (b)   (a) Case 1: In this case study, PV units and WTs, as the considered renewable DGs, are allocated at certain predetermined buses of the 59-bus distribution network as listed below in Table 5.The optimal sizing of these DGs is determined, using the SHADE algorithm, in addition to the optimal network configuration, considering different scenarios of load demand, wind speed and solar irradiance for maximizing the probabilistic hosting capacity (PHC) of the network, minimizing the network power losses (i.e., maximizing the reduced power losses), and improving the voltage profile.This case study is conducted in order to provide a feasible comparison with the previously presented study in [103] using multi-objective optimization techniques (NSGA-II, MOPSO, MOMVO and MOFPA).The obtained simulation results of this case study regarding the PHC, RP loss and voltage profile, in comparison to those presented in [103], are listed below in Table 6.Optimal sizing of the allocated PV units and WTs, in MW, is provided below in Table 7.In addition, the tie-switches selected by the SOE method for the optimal network reconfiguration for various scenarios are supplied below in Table 8.Finally, the voltage profile improvement after the optimal DG integration and network reconfiguration for various scenarios is shown below in Figure 3.(b) Case 2: In this case study, optimal allocation and sizing of two shunt capacitors in the considered 59-bus distribution network are added to the optimization problem of the previous case study (Case 1) to step on the effect of integrating SCs on system performance.The obtained simulation results of this case study regarding the PHC, RP loss and voltage profile are listed below in Table 9. Optimal locations and sizes of SCs, in MVAR, in the distribution network in addition to optimal sizing of the allocated PV units and WTs, in MW, are provided below in Table 10.Besides, the tie-switches selected by the SOE method for the optimal network reconfiguration, for various scenarios, are supplied below in Table 11.Finally, the voltage profile improvement after the optimal DG and SC integration and network reconfiguration for various scenarios is shown below in Figure 4. (c) Case 3: In this case study, the optimal allocation and sizing of three DGs and three SCs in the considered 33-bus distribution network are determined, using the SHADE algorithm simultaneously with the optimal network configuration, using the SOE method, without considering the demand load variability or the DG output power uncertainty for minimizing network power loss and improving the voltage profile.The simulation results obtained from this case study are compared with those provided by the previously presented studies in [104,105], using the DE and the BPSO optimization algorithms, and listed below in Table 12.Voltage profile enhancement after the optimal DG and SC integration in addition to optimal network reconfiguration is shown below in Figure 5.(d) Case 4: In this case study, the variability or uncertainty associated with the network demand load is considered based on different loading scenarios, as indicated in Table 1, in solving the optimization problem of the previous case study (Case 3) in order to investigate the effect of considering realistic variable loads on the optimization results.
The obtained simulation results of this case study are listed below in Table 13.Voltage profile improvement after the optimal DG and SC integration in addition to optimal network reconfiguration is as depicted below in Figure 6.(e) Case 5: In this case study, the variability or uncertainty associated with both the network demand load and DG output power is considered in solving the optimization problem of Case 3. One wind turbine and two PV units are considered as the DG units in this case study with different scenarios for wind speed and solar irradiance as listed in Table 1.This case study is performed to investigate the effect of realistic uncertain demand load and DG output power on optimization results.The obtained simulation results of this case study are provided below in Table 13.Voltage profile enhancement is shown below in Figure 7 after the optimal DG and SC integration along with optimal network reconfiguration.For the problem of case 1, as shown above in Tables 6 and 7, it has been found that using the SHADE optimization algorithm and the SOE reconfiguration technique for the different considered scenarios, the probabilistic hosting capacity () of the network reached 16.851%, the reduced network power losses reached 93.99%, while the voltage deviation index ( ) decreased to 0.0814.From Table 7, it can be noted that the optimal sizes of WTs allocated at buses 24 and 31 (0.58 MW and 0.91 MW), obtained using the SHADE algorithm, are less than those obtained using the NSGA-II algorithm (5.5 MW and 1.7 MW) and the MOPSO algorithm (1 MW for both WTs).Besides, the optimal sizes of PV units allocated at buses 7 and 50 (3.9MW and 1.7 MW), obtained by the SHADE algorithm, are less than those obtained by the MOPSO algorithm (4.4 MW and 12.2 MW), while the obtained optimal size of the PV unit allocated at bus 22 (1.8 MW), using the SHADE algorithm, is less than that obtained using both the NSGA-II algorithm (2.3 MW) and the MOPSO algorithm (2.9 MW).Despite that, the optimal sizes of PV units and WTs at certain buses of the network, obtained by the SHADE algorithm, are less than those obtained using other techniques as indicated above; the employed SHADE optimization algorithm, in association with the SOE reconfiguration method, achieves PHC of 16.85%, which is almost touching that provided by the NSGA-II algorithm (18.08%) and it exceeds that achieved using MOPSO, MOMVO, and MOFPA algorithms.It is also notable, from Table 6, that the SHADE and the SOE methods achieve the maximum reduced power losses (93.995%) which surpasses the values provided using NSGA-II, MOPSO, MOMVO, and MOFPA algorithms by 10.7%, 12.2%, 21.1%, and 12.8%, respectively.In addition, the NSGA-II algorithm achieves the maximum value for the minimum bus voltage (0.9944 p.u.), then the SHADE algorithm (0.9938 p.u.), then the MOPSO algorithm (0.9931 p.u.), while the MOPSO algorithm achieves the maximum value for the maximum bus voltage (1.0029 p.u.), then the SHADE algorithm (1.0013 p.u.), then the NSGA-II algorithm (1.0004 p.u.).However, the SHADE algorithm, along with the SOE method, achieves the best (lowest) value for the overall aggregated voltage deviation index (0.0814).In view of those aforementioned results, it can be stated that the SHADE optimization algorithm and the SOE reconfiguration method provide the best performance and solutions, compared to other techniques, which indicates their effectiveness and their superiority for solving the considered optimization problem.The effect of adding optimal allocation and sizing of SCs to the optimization problem of case 1 has been investigated in the second case study (case 2) and the results, listed in Tables 9 and 10, indicate that optimal integration of SCs and DGs with the optimal sizes and locations listed in  For the problem of case 1, as shown above in Tables 6 and 7, it has been found that using the SHADE optimization algorithm and the SOE reconfiguration technique for the different considered scenarios, the probabilistic hosting capacity (PHC) of the network reached 16.851%, the reduced network power losses reached 93.99%, while the voltage deviation index (AVDI ov ) decreased to 0.0814.From Table 7, it can be noted that the optimal sizes of WTs allocated at buses 24 and 31 (0.58 MW and 0.91 MW), obtained using the SHADE algorithm, are less than those obtained using the NSGA-II algorithm (5.5 MW and 1.7 MW) and the MOPSO algorithm (1 MW for both WTs).Besides, the optimal sizes of PV units allocated at buses 7 and 50 (3.9MW and 1.7 MW), obtained by the SHADE algorithm, are less than those obtained by the MOPSO algorithm (4.4 MW and 12.2 MW), while the obtained optimal size of the PV unit allocated at bus 22 (1.8 MW), using the SHADE algorithm, is less than that obtained using both the NSGA-II algorithm (2.3 MW) and the MOPSO algorithm (2.9 MW).Despite that, the optimal sizes of PV units and WTs at certain buses of the network, obtained by the SHADE algorithm, are less than those obtained using other techniques as indicated above; the employed SHADE optimization algorithm, in association with the SOE reconfiguration method, achieves PHC of 16.85%, which is almost touching that provided by the NSGA-II algorithm (18.08%) and it exceeds that achieved using MOPSO, MOMVO, and MOFPA algorithms.It is also notable, from Table 6, that the SHADE and the SOE methods achieve the maximum reduced power losses (93.995%) which surpasses the values provided using NSGA-II, MOPSO, MOMVO, and MOFPA algorithms by 10.7%, 12.2%, 21.1%, and 12.8%, respectively.In addition, the NSGA-II algorithm achieves the maximum value for the minimum bus voltage (0.9944 p.u.), then the SHADE algorithm (0.9938 p.u.), then the MOPSO algorithm (0.9931 p.u.), while the MOPSO algorithm achieves the maximum value for the maximum bus voltage (1.0029 p.u.), then the SHADE algorithm (1.0013 p.u.), then the NSGA-II algorithm (1.0004 p.u.).However, the SHADE algorithm, along with the SOE method, achieves the best (lowest) value for the overall aggregated voltage deviation index (0.0814).In view of those aforementioned results, it can be stated that the SHADE optimization algorithm and the SOE reconfiguration method provide the best performance and solutions, compared to other techniques, which indicates their effectiveness and their superiority for solving the considered optimization problem.The effect of adding optimal allocation and sizing of SCs to the optimization problem of case 1 has been investigated in the second case study (case 2) and the results, listed in Tables 9 and 10, indicate that optimal integration of SCs and DGs with the optimal sizes and locations listed in Table 10, simultaneously with the optimal network reconfiguration, lead to increasing the network PHC to be 17.1% instead of 16.85% provided in the first case study (without SCs), increasing the reduced network power losses (RP loss ) to 95.3% instead of 93.9% and decreasing the voltage deviation index (AVDI ov ) to 0.0724 instead of 0.0814.It can be noted that the optimal sizes of WTs allocated at buses 13 and 55 (0.95 MW and 2.08 MW), in case 2, are less than those obtained in case 1 (3.3 MW and 3 MW).Besides, the optimal sizes of PV units allocated at buses 2, 7, and 43 (0.29 MW, 1.3 MW and 3.3 MW), in case 2, are less than those obtained in case 1 (3.1 MW, 3.9 MW and 4.1 MW).However, the network PHC obtained in case 2 (17.089%) is larger than that obtained in case 1 (16.851%).It is also notable that despite the minimum bus voltage obtained from case 1 being slightly less than that obtained in case 2, the aggregated voltage deviation index (AVDI ov ) obtained in case 2 is less than that obtained in case 1.Hence, it can be concluded that optimal integration of SCs in addition to DGs and network reconfiguration has a positive effect on improving system performance.For the problem of case 3, as listed in Table 12, it has been found that using the SHADE algorithm along with the SOE method, minimum network power losses were achieved (12.7 kW) compared to other techniques, namely DE (15.63 kW) and BPSO (15.47 kW).In other words, the SHADE algorithm in association with the SOE method provide the maximum reduced power loss (93.7%) compared to DE (92.3%) and BPSO (92.4%).It can also be noted that the SHADE algorithm provides the maximum value for the minimum bus voltage (0.9936 p.u.) compared to DE (0.9891 p.u.) and BPSO (0.9887 p.u.).Those demonstrated results in Table 12 indicate the efficiency and the effectiveness of the employed SHADE algorithm and the SOE method in addition to their superiority over other used techniques for solving the considered optimization problem.The obtained simulation results of cases 4 and 5, as indicated above in Table 13, show that considering the variability and uncertainty associated with the network demand load and DG output power in solving the optimization problem of case 3 has a major effect on the optimization results which notably differ from those obtained in case 3 (without considering uncertainties).In case 4, the obtained network power loss (probabilistic power losses) decreased to 2.75 kW (98.6% power loss reduction) with the optimally integrated DGs and SCs of total capacity 1.04 MW and 0.64 MVAR; respectively, while in case 5 the obtained network power loss has decreased to 7.30 kW (96.4% power loss reduction) with optimally integrated DGs and SCs of total capacity 1.50 MW and 2.12 MVAR, respectively, instead of 12.7 kW obtained power loss (93.7% power loss reduction) with optimally integrated DGs (2.89 MW total allocated capacity) and SCs (1.69 MVAR total allocated capacity) in case 3, where the uncertainties are neglected.The aforementioned results of cases 4 and 5, in comparison to those of case 3, highlight the importance and effect of considering the uncertainties related to network load and DG output power for realistic modeling of the considered optimization problem and obtaining optimal solutions for the different network operational scenarios [106].

Conclusions and Future Works
In the work presented in this paper, the SHADE optimization algorithm has been employed for optimal integration of renewable DGs (WTs and PV units) and shunt capacitors in the considered distribution networks, simultaneously with the optimal network reconfiguration using the SOE method and with consideration of demand load and DG output power uncertainties.Optimization objectives have included maximizing the probabilistic hosting capacity of the network, minimizing the network power losses, and enhancing the voltage profile.Five different case studies were conducted, considering the real 59-bus distribution network of Cairo and the 33-bus test distribution network.The first case study has discussed optimal DG integration simultaneously with optimal network reconfiguration for the 59-bus distribution network and considering load and DG output uncertainties.The employed SHADE algorithm and the SOE method have provided maximum hosting capacity (16.851%) and reduced power losses (93.99%) as well as minimum voltage deviation index (0.0814) compared to techniques presented in other previous literature [103]-as indicated in Table 6, which validates their effectiveness and superiority.The second case study has investigated the effect of adding SC integration to the optimization problem of the first case study, which is the main contribution of this paper.As shown in Table 9, SC integration has led to improving (increasing) hosting capacity (17.1%) and reduced power loss (95.3%) in addition to further improvement of the voltage profile with an improved (reduced) voltage deviation index (0.0724).The third case study presented optimal integration of DGs and SCs in the 33-bus distribution network, simultaneously with the network reconfiguration and without considering the uncertainties related to demand load and DG output.The employed SHADE algorithm and the SOE method, as listed in Table 12, provided the best results, compared to other techniques used in previously presented research [104,105] regarding power loss reduction (93.7% compared to 92.3% and 92.4%) and voltage profile improvement (0.994 min.bus voltage compared to 0.989 and 0.988).The fourth and fifth case studies investigated the effect of considering uncertainties related to network load and DG output power in solving the optimization problem of case 3.These case studies also reflect the main contribution of the work presented in this paper as the optimal integration of DGs and SCs, simultaneously with network reconfiguration and with consideration of demand load and DG output uncertainty, hasn't been presented in the available literature.The obtained results, as depicted in Table 13, have indicated the major effect of considering such uncertainties on the optimization results to provide optimal solution for the different network operational scenarios with realistic modeling of the considered optimization problem.There are some aspects that haven't been considered in work presented within this research and can be addressed in future works.These aspects include integration of a suitable type energy storage system (ESS) into the considered distribution network to support the renewable DGs in case of low or zero output power, in addition to investigating the effect on system performance.Besides, system reliability assessment and considering other indices for evaluating distribution system performance.Future works could also include considering other distribution networks with larger scale and employing other optimization techniques with a comparative study between different considered techniques.

Algorithm 2 : 1 G
SHADE Optimization Algorithm Initialization Phase = 0; 2 Initialize population P 0 3 Set all values in M CR and M F to 0.5; 4 Archive A = ϕ; 5 Index counter k = 1; Main Loop 6 While the termination criteria are not met do 7 S CR = ϕ, S F = ϕ; 8 for i = 1 to N do 9 r i = Select from [1, H] randomly; 10 CR i,G = randn i (M CR,r i , 0.1); 11

Figure 3 .
Figure 3. Voltage profile of the 59-bus distribution network for various scenarios (Case 1).

Figure 3 .
Figure 3. Voltage profile of the 59-bus distribution network for various scenarios (Case 1).

Figure 4 .
Figure 4. Voltage profile of the 59-bus distribution network for various scenarios (Case 2).

Figure 4 .
Figure 4. Voltage profile of the 59-bus distribution network for various scenarios (Case 2).

Figure 7 .
Figure 7. Voltage profile of the 33-bus distribution network for various scenarios (Case 5).

Figure 7 .
Figure 7. Voltage profile of the 33-bus distribution network for various scenarios (Case 5).

Table 1 .
Resulting scenarios from the applied scenario reduction technique (30 Scenarios).

Table 2 .
Considered distribution networks basic data.

Table 3 .
Considered renewable DGs and shunt capacitors input parameters.

Table 5 .
PV units and WTs buses of the 59-bus distribution network (Case 1).

Table 6 .
Simulation results of the first case study (Case 1).

Table 7 .
Optimal sizing of PV units and WTs in MW (Case 1).

Table 8 .
Selected tie-switches for optimal network configuration for various scenarios (Case 1).

Table 8 .
Selected tie-switches for optimal network configuration for various scenarios (Case 1).

Table 9 .
Simulation results of the second case study (Case2).

Table 10 .
Optimal allocation and sizing of PV units, WTs and SCs in MW / MVAR (Case 2).

Table 9 .
Simulation results of the second case study (

Table 10 .
Optimal allocation and sizing of PV units, WTs and SCs in MW / MVAR (Case 2).

Table 11 .
Selected tie-switches for optimal network configuration for various scenarios (Case 2).

Table 11 .
Selected tie-switches for optimal network configuration for various scenarios (Case 2).

Table 12 .
Simulation results of the third case study (Case 3).

Table 12 .
Simulation results of the third case study (

Table 13 .
Simulation results of the fourth and fifth case studies (Cases 4 & 5).

Table 13 .
Simulation results of the fourth and fifth case studies (Cases 4 & 5).

Table 13 .
Simulation results of the fourth and fifth case studies (Cases 4 & 5).