Artificial Electric Field Algorithm-Pattern Search for Many-Criteria Networks Reconfiguration Considering Power Quality and Energy Not Supplied

: Considering different objectives and using powerful optimization methods in the distribution networks reconﬁguration by accurately achieving the best network conﬁguration can further improve network performance. In this paper, reconﬁguration of radial distribution networks is performed to minimize the power loss, voltage sag, voltage unbalance, and energy not supplied (ENS) of customers using a new intelligent artiﬁcial electric ﬁeld algorithm-pattern search (AEFAPS) method based on the many-criteria optimization approach. The voltage sag and voltage unbalance are deﬁned as power quality indices and the ENS is the reliability index. In this study, the pattern search (PS) algorithm enhances the artiﬁcial electric ﬁeld algorithm’s (AEFA) ﬂexibility search both globally and locally. AEFAPS is applied to determine the decision variables as open switches of the networks considering the objective function and operational constraints. The proposed methodology based on AEFAPS is performed on an unbalanced 33-bus IEEE standard network and a real unbalanced 13-bus network. The reconﬁguration problem is implemented in single-criterion and many-criteria optimization approaches to evaluate the proposed methodology’s effectiveness using different algorithms. The single-criterion results demonstrated that some power quality indices might be out of range, while all indices are within the permitted range in the many-criteria optimization approach, proving the effectiveness of the proposed many-criteria reconﬁguration with logical compromise between different objectives. The results show that AEFAPS identiﬁed the network conﬁguration optimally and different objectives are improved considerably compared to the base network. The results conﬁrmed the superior capability of AEFAPS to obtain better objective values and lower values of losses, voltage sag, voltage unbalance, and ENS compared with conventional AEFA, particle swarm optimization (PSO), and grey wolf optimizer (GWO). Moreover, the better performance of AEFAPS is proved in solving the reconﬁguration problem compared with previous studies.


Introduction
Several studies have been performed in the past decade to improve the distribution network qualities through loss reduction, voltage deviations reduction, and voltage stability improvement in these networks [1]. Primary methods that suggest loss reduction in the distribution networks include distributed generations (DGs), the optimal use of capacitors, and network reconfiguration [2]. Network reconfiguration is one of the important benefits of network automation. Reconfiguration is a low-cost method for reducing the losses of distribution networks because this method has less need for additional equipment and high costs to reduce losses induced by the capacitor methods and use of DG resources [3]. Reconfiguration changes the distribution network topology by modifying the state of the applied for reconfiguration problem-solving. Many-criteria decision-making optimizationbased decomposition is a well-known and popular strategy for satisfactory performance and desirable efficiency of its implementation due to considering each objective as a vector. The meta-heuristic algorithm can solve the optimization problem by comprising all objectives into a participatory form.  [10] Not considering the power quality and reliability indices [4] Not evaluating the power quality and reliability indices [11] Without the power quality and reliability indices [12] Not studying the power quality and reliability indices [13] Not incorporating power quality and reliability indices [14] Without the power quality indices [15] Without the power quality and reliability indices [16] Not investigating the power quality and reliability indices [17] Without the power quality and reliability indices [18] Not considering the power quality and reliability indices [19] Without the power quality and reliability indices [20] Without studying the power quality and reliability indices [21] Without the power quality and reliability indices [22] Not considering the voltage sag and unbalanced network [23] Not considering the power quality and reliability indices [24] Without the power quality and reliability indices This paper Considering power quality, reliability, unbalanced network, and hybrid algorithm In this study, unbalanced networks reconfiguration is presented considering the manycriteria method-based decomposition with the aim of loss minimization, enhancement of the voltage sag and voltage unbalance (enhancing the power quality), and reducing the customer's energy not supplied (reliability enhancement) using the new AEFAPS metaheuristic intelligent algorithm. The decision-making variable of the problem is considered network open switches (tie-lines), which are determined by the AEFAPS method. Therefore, the best configuration of the network is found because of open and closed switches based on the proposed method that achieves the minimum loss and the highest reliability and power quality. This paper evaluates the results of loss, power quality indices, and reliability with and without many criteria for decision-making reconfiguration of unbalanced distribution networks using AEFAPS. Moreover, the superiority of AEFAPS is compared with wellknown PSO [29] and grey wolf optimizer (GWO) [30] in problem-solving. The highlights of this study are as follows: • Network automation based on optimal reconfiguration; • Network reconfiguration with many-criteria decision-making optimization; • Using artificial electric field algorithm-pattern search (AEFAPS) for problem-solving; • Using pattern search in the suggested algorithm to improve global and local searches; • Superiority of AEFAPS compared with conventional AEFA, GWO, and PSO methods.
The rest of the paper is structured as follows: Section 2 describes the formulation of the problem. Section 3 describes the load flow problem and the branch exchange method. Section 4 presents the optimization algorithm. Section 5 presents the results of the findings. Section 6 concludes with the major findings of the paper.

Problem Formulation
This study considers many-criteria decision-making optimization-based decomposition for unbalanced distribution networks reconfiguring using AEFAPS. The formulation of objective function and constraints are presented below.

Objective Function
The objective function based on many criteria is defined considering losses, voltage sag and unbalance, and ENS.

Power Losses
The power loss is calculated by determining the buses' voltage and also the line current by [11]. (1) where R k and I k indicate the ohmic resistance and current magnitude of the k-th line, respectively, and l refers to the line's number. Equation (1) is related to the active power losses of the distribution network based on the ohmic resistance and also current of the network lines.

Voltage Sag
This objective is achieved by evaluating the bus's remaining voltage in the voltage sag condition. The average voltage sag in all buses is applied to enhance this objective of the entire network [4,6,10].
where V j i is bus voltage magnitude i in the condition of the fault in bus j; i = 1, 2, . . . , n refers to the number of buses; j = 1, 2, . . . , m refers to the possible faults number, and V sag.av denotes the bus voltage average in the condition of the voltage sag.
The voltage sag model is implemented via fault performing in the network. The fault current for a fault at bus p is based on the admittance matrix with the fault (Y a,b,c F ), self-impedance of bus p (Z a,b,c pp ), and pre-fault voltage at bus p (E a,b,c p(o) ) and is determined as follows: where U refers to the unity matrix. Voltage drop is computed as follows: The rest of the bus voltage is computed as follows. The distribution network operation must minimize the voltage drop based on reconfiguration.
where E a,b,c Bus(o) is the bus voltages magnitude without fault as three-phase, which is achieved with power flow in system frequency.

Voltage Unbalance
The unbalanced network is due to load unbalanced. The network unbalance value is computed by calculating the voltage unbalance of the buses.
Voltage unbalance in three-phase networks is caused by network load unbalance, which is defined by the sum of voltage unbalances of all network buses as follows [31]: where V Pos,i indicates the bus i positive sequence and V Neg,i refers to the bus i negative sequence voltage magnitude.

Reliability
The objective function of the network reliability improvement or reduction of customers' energy not supplied (ENS) is defined considering the load point and its energy consumption. The electrical power of some load points may be cut off due to a line outage. Therefore, based on the probability of the network line outage, the amount of power lost by network customers is defined as ENS, which aims to reduce this index and improve network reliability. Based on reconfiguration, the ENS value can be reduced by changing the position of the switches and consequently by changing the path of the power flow in the case of a network line outage. The ENS should be minimized and its value is calculated by [32]: where nl is the load point, P i is an average of active power at load point i, λ i represents the line outage rate, U i is the unavailability of load point i, and r i indicates average outage time of the load point i.

Constraints
The many-criteria objective function should be optimized subject to the constraints below. These constraints are addressed under two categories: equality and inequality constraints [12].

•
Power balance • Voltage constraint The bus's voltage should be in the allowed range.
• Lines constraint The current of each line should not exceed a certain value as the maximum allowable current (I max pl ). I pl ≤ I max pl , p = a, b, c, l = 1, 2, 3, . . . , L • Voltage unbalance constraint The voltage unbalance should not exceed a certain value (V max unb ).
• Voltage sag constraint The voltage sag should not exceed a certain value (V min sag ).
V sag,av ≥ V min sag (13) • Network radiality In the reconfiguration process, the radial configuration of the network must be guaranteed and does not have a ring mode configuration. Matrix A has one row for each network line and one column for each network bus, with array a ij in row i and column j according to the following rule [16]: The matrix of bus incidence [16] is presented in the literature for a radiality check of the network during the reconfiguration process as follows: Radial and det(A) = 0 Not radial (14) where A is the bus incidence matrix.

Many-Criteria Decision-Making Optimization
In problem-solving based on the fuzzy method, each objective is presented as a fuzzy membership function. The fuzzy values given by the membership functions to each objective function are between 0 and 1 and depend on how much improvement has been achieved. In many-objective optimization, a compromise should be made between different objectives. Therefore, the fuzzy index (µ) is defined for each objective. In [28], the authors show how each of the fuzzy indices for the objectives is formulated which is shown in Figure 1. Since the four objectives, such as the loss, ENS, voltage sag, and unbalance, must be minimized, each of the corresponding membership functions must be defined. Moreover, for loss, ENS, voltage unbalance, and voltage sag, the membership function is formulated as follows: Since the four objectives, such as the loss, ENS, voltage sag, and unbalance, must be minimized, each of the corresponding membership functions must be defined. Moreover, for loss, ENS, voltage unbalance, and voltage sag, the membership function is formulated as follows: where x max and x min are the upper and lower values of each objective, respectively. x max is the objective value in the base network, while x min is extracted as single-objective optimization with reconfiguration; thus, x is placed among x max and x min . As expressed by Equations (15)- (18), the nearest objective value to x max has a lower µ. Suppose the minimum of all µ i is maximized. In that case, it means moving to the minimum value of all objectives and the minimization of different objectives is created simultaneously, which is formulated as follows [28]: where µ f 1 (x), µ f 2 (x), . . . , µ f m (x) are the vectors of individual objective functions based on fuzzy indices, x refers to the decision variables vector, and h(x) and g(x) are constraints of equality and inequality, respectively.

Power Flow and Branch Exchange Method
Forward-backwards load flow is used to understudy the distribution networks. Each line is modeled by series impedance (resistance and inductance). A forward-backward sweeping iterative algorithm produces the load flow to achieve convergence and find the current of lines and bus voltages [12]. There are two basic steps in the forward-backward sweeping iterative algorithm; they are very easy to learn. These two steps are used for all radial networks and act as a ring. In the forward sweeper, first, the bus voltage is equal to 1 p.u and the current for each bus is easily obtained. In the forward sweeper, the current of each bus is obtained. The amount of voltage is updated from the end of the feeder. As stated, the voltage is set to 1 p.u and the bus flow is obtained. From this current, a new voltage is again obtained from the end of the feeder. The steps of load flow are as follows: The three-phase voltage for all phases is considered equal to the bus voltage of the slack (post). A load of each phase is determined as a positive sequence and 120 degrees of phase difference. In this step, the backward load flow is achieved. Start from the ending bus and proceed to the sending bus and update the upstream lines' current as follows: the current of the upstream branch equals the current of all the downstream lines that begin from the ending bus of the upstream branch (this step is repeated until the determination of the branches' current). The forward load flow is then achieved. Move from the sending to the ending bus of the network to calculate the ending bus voltage of each line as follows [13]:  branches' current). The forward load flow is then achieved. Move from the sending to the ending bus of the network to calculate the ending bus voltage of each line as follows [13]: where , , is the voltage of each phase of the bus (load side), which is referred to the lower bus ( Figure 2 To evaluate the radiality of the network, the branch exchange method has been applied. In order to maintain radiality, by closing one line another line should be opened so that the ring does not form and the network is still radial. This process is presented in Figure 2a, where t represents a tie-line. The load that connects to the bus n feeds through p-q-r lines. The power supply path to the bus n can be modified by closing the tie-line t and opening the line r to the p-s-t ( Figure 2b). The advantage of this method is that the radiality of the network is automatically met, and it requires no computing and other algorithms to detect the radial mode of the network [13].

Proposed Optimizer
The fuzzy many-criteria reconfiguration is presented for minimizing the losses, power quality, and reliability enhancement using the meta-heuristic AEFAPS algorithm. One of the optimization variables is the situation of distribution network opened lines determined using AEFAPS. In this section, the AEFAPS method is formulated. To evaluate the radiality of the network, the branch exchange method has been applied. In order to maintain radiality, by closing one line another line should be opened so that the ring does not form and the network is still radial. This process is presented in Figure 2a, where t represents a tie-line. The load that connects to the bus n feeds through p-q-r lines. The power supply path to the bus n can be modified by closing the tie-line t and opening the line r to the p-s-t ( Figure 2b). The advantage of this method is that the radiality of the network is automatically met, and it requires no computing and other algorithms to detect the radial mode of the network [13].

Proposed Optimizer
The fuzzy many-criteria reconfiguration is presented for minimizing the losses, power quality, and reliability enhancement using the meta-heuristic AEFAPS algorithm. One of the optimization variables is the situation of distribution network opened lines determined using AEFAPS. In this section, the AEFAPS method is formulated.

Overview of AEFA
The AEFA is modeled with Coulomb's law for electrostatic force. This law describes the electrostatic reactions between the electrical charges. In the conventional AEFA, the charged particles are selected as agents, and each agent's resistance is evaluated based on its charges. The conventional AEFA algorithm is modeled based on electrostatic attraction force. The first law of Coulomb says that the particles of the same name have repellent properties. Otherwise, the particles pull each other, with the second law of Coulomb also stating that there is an attractive force among opposing charges and repulsive force between same-name charges, which is directly related to the multiplication of the charges and conversely related to their distances. Suppose the position of the i-th particle considering d-dimension searching as ( , where x i d refers to the i-th particle position. The best position with the best fitness obtained from particles is determined. The best fitness position achieved by any particle i is defined by [26]: Optimal fitness is defined as P best = X best . The force of the charge i inserted by the charge j is calculated by [26]: where Q i (t) and Q j (t) refer to the charged particles of i and j, K(t) indicates the constant of Coulomb, ε is a small constant, and R ij (t) is defined as the Euclidean distance among two charged particles of i and j and is calculated by [26]: The K(t) is based on the iteration number and max iter action (max iter), which is defined by [26]: where α refers to the parameter and K 0 is the initial value. The inserted electric force on particle i via the other particles is defined in d search space at time t as follows [26]: where rand() indicates a uniform random number in [0, 1], N clears the number of particles, and F i refers to the force inserted into charged particle i. The particle (i) electrical force in d-th dimension search space is defined by: So, applying Newton's Second Law of Motion, the particle i acceleration is defined by [26]: M i (t) clears the particle mass of i considering the iteration t. The velocity of the charged particle and its position are updated by [27]: For the minimization or maximization problem, the fitness should have a downward or upward trend, respectively [26].
where fit i refers to the fitness of particle i. The values of best and worst fitness for fitness maximization are formulated by [27]: For the problem with the minimization approach, the values of best and worst fitness are presented by [26]:

Overview of AEFAPS
Direct search techniques are referred to the intelligent algorithms that are applied for solving constrained problems. One of these methods is the pattern search (PS) algorithm, a desirable algorithm for optimization problem solving [27]. The algorithm begins by creating the points around the initial starting point that can be provided by the user or derived from the calculation results of the previous steps. All processes for creating this set of points are called meshes, being mesh points performed by adding the initial point size to pattern vectors. Therefore, the algorithm steps are presented below: Step 1: PS begins with the initial point of X 0 given by the user.
Step Step 3: The initial point X 0 is added to the unit vectors, with the fitness computed for the obtained new points. The formation of the mesh is shown in Figure 3.

Overview of AEFAPS
Direct search techniques are referred to the intelligent algorithms that are applied for solving constrained problems. One of these methods is the pattern search (PS) algorithm, a desirable algorithm for optimization problem solving [27]. The algorithm begins by creating the points around the initial starting point that can be provided by the user or derived from the calculation results of the previous steps. All processes for creating this set of points are called meshes, being mesh points performed by adding the initial point size to pattern vectors. Therefore, the algorithm steps are presented below: Step 1: PS begins with the initial point of given by the user. Step 3: The initial point is added to the unit vectors, with the fitness computed for the obtained new points. The formation of the mesh is shown in Figure 3.  Step 4: At this step, it is determined whether the objective function value for one of these points is lower than the fitness value for the initial point of the X 0 . If such a point exists, it is said that a successful selection occurs, and the algorithm places this point as X 1 . Therefore, the mesh changes as follows: Step 5: If at each step the fitness value for that step is larger than the value of the objective functions generated by each of the generated X (mesh points), then there will be an unsuccessful selection. For this unsuccessful selection, the algorithm multiplies the size of the mesh by a decreasing factor (typically 0.5).
The algorithm continues to find an optimal solution to minimize the objective function and stops when one of the following factors occurs: (a) The size of the mesh is lowering its defined value. (b) The iteration number reaches the predetermined value. (c) The distance between points determined in a desirable selection and the point determined in the following desirable selection is lower than the allowable value. (d) Changing the fitness from a desirable selection to another successful selection is lower than the allowable value.
In AEFAPS, at first, the optimal variables determined by the conventional AEFA are applied as initial values for the PS method. The determination of the globally optimal values is continued using the PS method so that the best decision variables and the best objective function of the problem can be obtained. The flowchart of the AEFAPS method is shown in Figure 4.

The AEFAPS Method Implementation
The AEFAPS method is applied to solving the network reconfiguration. The steps involved in the AEFAPS implementation in solving the reconfiguration are as follows: Step 1: The number of optimization variables is equal to the number of tie-lines. The 13-and 33-bus networks have five tie-lines. Therefore, in opening the tie-lines, another line must be closed so that no loop is formed in the network and the radial state of the network is maintained.
Step 2: In the reconfiguration process, with the opening of the tie-lines, the network constraints, as well as the network radiality constraints, must be evaluated. The tie-line numbers must not be duplicated, otherwise, the selected variable set must be deleted. To evaluate the radiality state of the network, the bus incidence matrix (A) is used, in which the row number indicates the lines number and the column number indicates the number of buses. If line i is not connected to bus j, the matrix element is considered 0; if line i exits bus j, the matrix element is considered 1; and if line i enters bus j, the matrix element is considered +1. Finally, the network is radial if its determinants are not zero. The set of variables with matrix 0 is removed in the reconfiguration process.
Step 3: The objective function is computed for a set of decision variables (selected tie-lines) using the conventional AEFA by considering the network constraints by executing the load flow. The best network configuration is determined by the best set of variables with the best value of the objective function (minimum value).
Step 4: In this step, by updating the population of the conventional AEFA, steps 1 to 3 are performed for a set of new variables. Suppose the objective function value is less than step 3. In that case, the decision variables corresponding to the new configuration will be replaced.
Step 5: The best solution includes the best optimization variables set obtained by the conventional AEFA applied as initial values for the PS method. So, the best decision variables corresponding to the lowest objective function are determined using the PS method. If the value of the objective function is lower than the obtained value at step 4, the decision variables set corresponding to the new configuration will be replaced.
Step 6: Go to step 7 if the minimum objective function value with the algorithm's maximum iteration is achieved. Otherwise, go to step 4.
Step 7: Stop the algorithm and save the optimal variables (optimal network configuration).

Results and Discussion
The results of the 13-and 33-bus distribution networks reconfiguration are performed to minimize power losses, improving power quality and reliability using the combined AEFAPS method. The unbalanced 33-bus IEEE network has five lines between buses 8-21, 9-15, 12-22, 18-33, and 25-29 on which the proposed methodology is implemented. This network has an active load of 3715 kW, a reactive load of 2300 kVAr, and a voltage of 12.66 kV. The unbalanced 33-bus network information is taken from [33] and network reliability information, including outage rate and line repair time, is taken from [34]. The reconfiguration is also implemented on a real unbalanced 13-bus network based on the proposed methodology. This network has a 10,536 kW active load and a 5992 kVAr reactive load with a voltage of 20 kV. The 13-bus network has five tie-lines between buses 13-4, 6-9, 13-12, 12-10, and 9-10. The 13-bus network information is taken from [35]. Moreover, the reconfiguration results in different loadings, including nominal load and 62.5% and 125% of nominal load peak, and the capability of the AEFAPS method is evaluated. In addition, the superiority of AEFAPS is evaluated with conventional AEFA, well-known PSO, and GWO algorithms. The population number of algorithms, the maximum number of iterations based on user experience, repeated execution of the optimization program, and evaluation of results are equal to 50 and 100, respectively.

Results of Unbalanced 33-Bus Network
After performing the load flow for the base model of the 33-bus network without the reconfiguration, the results are given in Table 2. The values of power losses, minimum voltage, ENS, and voltage unbalance are calculated as 207.44 kW, 0.9566 p.u, 7.43 MWh and 4.92%, respectively. The problem of reconfiguration is performed as a single and fuzzy many-criteria, with the method results given. The optimization method's population is set to 50, and the maximum iteration is fixed at 100. The algorithms' parameters are found based on the trial-and-error method. The convergence process of different algorithms is depicted in Figure 5. This graph shows that AEFAPS obtained better fuzzy fitness than the conventional AEFA, PSO, and GWO methods. It is clear that conventional AEFA has many convergence tolerances and that the combination with the PS method improves its convergence process to achieve a better objective function value.

Results of Unbalanced 33-Bus Network
After performing the load flow for the base model of the 33-bus network without the reconfiguration, the results are given in Table 2. The values of power losses, minimum voltage, ENS, and voltage unbalance are calculated as 207.44 kW, 0.9566 p.u, 7.43 MWh and 4.92%, respectively. The problem of reconfiguration is performed as a single and fuzzy many-criteria, with the method results given. The optimization method's population is set to 50, and the maximum iteration is fixed at 100. The algorithms' parameters are found based on the trial-and-error method. The convergence process of different algorithms is depicted in Figure 5. This graph shows that AEFAPS obtained better fuzzy fitness than the conventional AEFA, PSO, and GWO methods. It is clear that conventional AEFA has many convergence tolerances and that the combination with the PS method improves its convergence process to achieve a better objective function value. The simulation results of the 33-bus network at nominal load are presented in Table  3, including optimal network configuration based on opened lines. The values of each objective include the power loss, voltage sag, voltage unbalance, and ENS considering single-criterion and many-criteria optimization using the AEFAPS algorithm. The best network combination is obtained in many-criteria optimization by opening 7, 11, 13, 28, and 32 switches based on the proposed AEFAPS. According to Table 3, the losses, voltage sag, voltage unbalance, and ENS is obtained at 146.96 kW, 0.653 p.u, 2.42%, and 6.69 MWh/year, respectively, using the AEFAPS based on many-criteria reconfiguration. The obtained results showed that the AEFAPS is superior compared to the other algorithm to obtain the lower values of different objectives. Moreover, the results demonstrated that single-criterion results could not generate a logical compromise between all the objectives. The simulation results of the 33-bus network at nominal load are presented in Table 3, including optimal network configuration based on opened lines. The values of each objective include the power loss, voltage sag, voltage unbalance, and ENS considering single-criterion and many-criteria optimization using the AEFAPS algorithm. The best network combination is obtained in many-criteria optimization by opening 7, 11, 13, 28, and 32 switches based on the proposed AEFAPS. According to Table 3, the losses, voltage sag, voltage unbalance, and ENS is obtained at 146.96 kW, 0.653 p.u, 2.42%, and 6.69 MWh/year, respectively, using the AEFAPS based on many-criteria reconfiguration. The obtained results showed that the AEFAPS is superior compared to the other algorithm to obtain the lower values of different objectives. Moreover, the results demonstrated that single-criterion results could not generate a logical compromise between all the objectives. To investigate the ability of the AEFAPS method, its performance in solving reconfiguration problems has been compared with the conventional AEFA, PSO, and GWO methods. According to the results, it is obvious that the value of each objective in optimizing it as a single criterion is less than the other objectives. For example, the minimum amount of power loss is in a single-criterion reconfiguration to minimize power loss. This is true for the other objectives (see Figures 6-9). The results support that in a single-criterion reconfiguration to minimize the voltage unbalance. The network power loss is close to the base value (it does not change much) and even increases the voltage sag, plus weakens the ENS (increased from 7.432 MWh to 7.810 MWh). Therefore, the results clarify that single-criterion reconfiguration cannot establish a logical and acceptable compromise between all the objectives. However, the results of many-criteria reconfiguration demonstrate that all the objectives of the problem have been established and improved by a compromise between them. In the AEFAPS method, opened network lines, including switches 7, 11, 13, 28, and 32, are determined as the optimal configuration of the 33-bus network. In the many-criteria reconfiguration, the power loss value decreased from 207.44 kW to 146.96 kW. The sag voltage value decreased from 0.9566 p.u in the base network to 0.6530 p.u, the voltage unbalance declined from 4.922% to 2.42%, and the ENS value decreased from 7.432 MWh to 6.69 MWh. Therefore, the reconfiguration based on many-criteria optimization using the AEFAPS method has improved all the objectives. Furthermore, comparing the results obtained from AEFAPS in solving the problem of single-criterion and many-criteria reconfiguration with conventional AEFA, PSO, and GWO methods shows the superiority of AEFAPS in achieving lower values of power losses, voltage sag, voltage unbalance, and ENS.

Before Reconfiguration
After Reconfiguration by HCSAPS Figure 9. The ENS in single-criterion and many-criteria optimizations using AEFAPS. Figures 10-12 show the variations of the sag voltage, voltage unbalance, and voltage profile of the 33-bus network before and after many-criteria reconfiguration based on the AEFAPS method, respectively. After many-criteria reconfiguration, the improvement of each mentioned objective confirmed the optimal performance of the proposed AEFAPSbased methodology.    . The ENS in single-criterion and many-criteria optimizations using AEFAPS. Figures 10-12 show the variations of the sag voltage, voltage unbalance, and voltage profile of the 33-bus network before and after many-criteria reconfiguration based on the AEFAPS method, respectively. After many-criteria reconfiguration, the improvement of each mentioned objective confirmed the optimal performance of the proposed AEFAPSbased methodology.

Effect of Load Increasing
The results of single-criterion and many-criteria reconfiguration of the unbalanced 33-bus network in terms of load increase equal to 125% of the network peak load based on the AEFAPS method are given in Table 4. The convergence process obtained from the proposed method, conventional AEFA, PSO, and GWO are shown in Figure 13. The superiority of the AEFAPS method in solving the many-criteria reconfiguration problem is confirmed with lower convergence tolerance and achieving less objective function value than the other methods.

Effect of Load Increasing
The results of single-criterion and many-criteria reconfiguration of the unbalanced 33-bus network in terms of load increase equal to 125% of the network peak load based on the AEFAPS method are given in Table 4. The convergence process obtained from the proposed method, conventional AEFA, PSO, and GWO are shown in Figure 13. The superiority of the AEFAPS method in solving the many-criteria reconfiguration problem is confirmed with lower convergence tolerance and achieving less objective function value than the other methods.  Table 4 shows that in the single-criterion reconfiguration to minimize voltage unbalance, the power loss, voltage sag, and ENS have been weakened and increased compared with the base state poor performance of the single-criterion approach. In contrast, the results obtained from many-criteria reconfiguration indicate the creation of a desirable compromise between different objectives and thus the improvement of all objectives based on the AEFAPS method compared to the base network. In the AEFAPS method, opened network lines including switches 7, 10, 28, 32, and 34 are determined optimal network configuration. In many-criteria reconfiguration, the power loss, voltage sag, voltage unbalance, and the ENS values are obtained at 232.98 kW, 0.833 p.u, 2.932%, and 8.074 MWh, respectively. In addition, according to Table 4, the AEFAPS method has performed better in single-criterion and many-criteria reconfiguration in terms of lower values of power losses, voltage sag, voltage unbalance, and ENS values compared to the conventional AEFA, PSO, and GWO methods. Figure 12. Voltage profile of a 33-bus network using AEFAPS before and after fuzzy many-criteria reconfiguration (nominal load).

Effect of Load Increasing
The results of single-criterion and many-criteria reconfiguration of the unbalanced 33-bus network in terms of load increase equal to 125% of the network peak load based on the AEFAPS method are given in Table 4. The convergence process obtained from the proposed method, conventional AEFA, PSO, and GWO are shown in Figure 13. The superiority of the AEFAPS method in solving the many-criteria reconfiguration problem is confirmed with lower convergence tolerance and achieving less objective function value than the other methods.

Effect of Load Decreasing
The results of single-criterion and many-criteria reconfiguration of the 33-bus distribution network in terms of network load reduction (62.5% of network nominal load) are presented in Table 5. The convergence process of the optimization methods is shown in Figure 14. As can be observed, the AEFAPS method achieved a lower fuzzy fitness compared to the other methods with less convergence tolerance. Therefore, it succeeded in enhancing the probability of achieving the optimal response. In this case, in the same way as the nominal load conditions and the increase in network load, it is observed that the loss values, power quality indices, and ENS values obtained in many-criteria optimization are the closest values to the best values in single-criterion optimization obtained by AEFAPS. The AEFAPS results were compared with the conventional AEFA, PSO, and GWO methods' performances, which confirmed the superiority of AEFAPS given an improvement in the objectives regarding loss reduction, power quality, and reliability enhancement.

Comparison of the Effect of Load Variations
The results of many-criteria reconfiguration in different loading conditions as losses, power quality indices, and ENS and improvement percentage of each objective are presented in Table 6 compared with the base state before reconfiguration. The results showed that as demand increased, net loss also increased. In addition, the value of the voltage sag and unbalanced voltage also increased with the increasing demand, and vice versa. Furthermore, with the increase in network load, the reliability objective weakened. In the nominal load network, the objectives of losses, power quality indices, and ENS are improved by 29.15%, 31.73%, 50.83%, and 9.98%, respectively. The greatest improvement is for unbalanced voltage and the least improvement is for ENS. Table 6 shows that when increasing and decreasing the network load, similar to the nominal load condition, the most improvement is for unbalanced voltage and the least is for ENS. In Figure 15, the improvement percentage of each objective is presented in different loadings. It is clear that increasing the network load demand increased loss and weakened power quality and reliability. Figure 16 shows that, in fuzzy many-criteria reconfiguration, the highest improvement percentage is obtained on voltage unbalance and the least improvement percentage is on network reliability. In Figure 15, the improvement percentage of each objective is presented in different loadings. It is clear that increasing the network load demand increased loss and weakened power quality and reliability. Figure 16 shows that, in fuzzy many-criteria reconfiguration, the highest improvement percentage is obtained on voltage unbalance and the least improvement percentage is on network reliability. The network voltage profile for various loadings using the AEFAPS method for an unbalanced 33-bus network is shown in Figure 16. The graph shows that with demand increasing, voltage deviations increased and the network voltage profile weakened, and vice versa.   The network voltage profile for various loadings using the AEFAPS method for an unbalanced 33-bus network is shown in Figure 16. The graph shows that with demand increasing, voltage deviations increased and the network voltage profile weakened, and vice versa. The network voltage profile for various loadings using the AEFAPS method for an unbalanced 33-bus network is shown in Figure 16. The graph shows that with demand increasing, voltage deviations increased and the network voltage profile weakened, and vice versa.

Results of Unbalanced 13-Bus Network
The many-criteria reconfiguration is performed to minimize the losses, voltage sag, voltage unbalance, and ENS on a real 13-bus network. The convergence process of AEFAPS in various loadings of the many-criteria reconfiguration is shown in Figure 17. Moreover, the reconfiguration results are presented in different loadings based on the AEFAPS method in Tables 7-9.

Results of Unbalanced 13-Bus Network
The many-criteria reconfiguration is performed to minimize the losses, voltage sag, voltage unbalance, and ENS on a real 13-bus network. The convergence process of AEFAPS in various loadings of the many-criteria reconfiguration is shown in Figure 17. Moreover, the reconfiguration results are presented in different loadings based on the AEFAPS method in Tables 7-9. The obtained results showed that according to Tables 7-9, the objectives of losses, voltage sag and unbalance, and ENS of the 13-bus network are enhanced after reconfiguration. The value of the net losses in the base case in 100% of the network load is 175.58 kW, which after reconfiguration reached 127.5 kW. The switches 6,8,9,10, and 14 are determined as open switches by AEFAPS. Moreover, the voltage sag, voltage unbalance, and ENS values before reconfiguration are 0.980 p.u, 0.81%, and 4.348 MWh, respectively, which decreased to 0.31 p.u, 0.552%, and 1.860 kWh, respectively, after many-criteria reconfiguration. The power quality indices variations without and with many-criteria reconfiguration are displayed in Figures 18 and 19, respectively. It is clear that the voltage sag and voltage unbalance values are reduced after reconfiguration in different network buses due to the positive effect of AEFAPS on the power quality of the 13-bus network. The obtained results showed that according to Tables 7-9, the objectives of losses, voltage sag and unbalance, and ENS of the 13-bus network are enhanced after reconfiguration. The value of the net losses in the base case in 100% of the network load is 175.58 kW, which after reconfiguration reached 127.5 kW. The switches 6,8,9,10, and 14 are determined as open switches by AEFAPS. Moreover, the voltage sag, voltage unbalance, and ENS values before reconfiguration are 0.980 p.u, 0.81%, and 4.348 MWh, respectively, which decreased to 0.31 p.u, 0.552%, and 1.860 kWh, respectively, after many-criteria reconfiguration. The power quality indices variations without and with many-criteria reconfiguration are displayed in Figures 18 and 19, respectively. It is clear that the voltage sag and voltage unbalance values are reduced after reconfiguration in different network buses due to the positive effect of AEFAPS on the power quality of the 13-bus network.  Figure 18. Voltage sag before and after many-criteria reconfiguration of an unbalanced 13-bus network using AEFAPS. Figure 19. Voltage unbalances before and after many-criteria reconfiguration of an unbalanced 13bus network using AEFAPS.

Comparison the Results with Previous Studies
The superior ability of AEFAPS in many-criteria reconfiguration of unbalanced 13-and 33-bus networks to minimize power losses and improve power quality and reliability indices    Voltage unbalances before and after many-criteria reconfiguration of an unbalanced 13bus network using AEFAPS.

Comparison the Results with Previous Studies
The superior ability of AEFAPS in many-criteria reconfiguration of unbalanced 13-and 33-bus networks to minimize power losses and improve power quality and reliability indices After Reconfiguration Before Reconfiguration Figure 19. Voltage unbalances before and after many-criteria reconfiguration of an unbalanced 13-bus network using AEFAPS.

Comparison the Results with Previous Studies
The superior ability of AEFAPS in many-criteria reconfiguration of unbalanced 13and 33-bus networks to minimize power losses and improve power quality and reliability indices is confirmed considering different loadings compared to conventional AEFA, PSO, and GWO methods as given in Table 10. The performance of AEFAPS is analyzed and compared with previous studies. In [33], the optimal reconfiguration of the unbalanced 33-bus network is performed to minimize power losses based on adaptive ant colony optimization (AACO). Moreover, in [13], reconfiguration of the 33-bus network has been developed to improve power quality using the MINLP approach. A comparison of the proposed method's results in terms of power losses with AACO and MINLP methods has shown that AEFAPS has fewer losses than the other methods, indicating its better performance. In [34], the balanced 33-bus network reconfiguration to minimize power losses and ENS is presented using ant colony optimization (MACO). According to Table 11, the ENS value is obtained using AEFAPS for the unbalanced 33-bus network (unbalanced) equal to 6.48 MWh and the balanced 33-bus network (balanced) equal to 6.21 MWh. A comparison of Table 11 shows that AEFAPS performs better than the MACO method in [34] in achieving a lower ENS value for the balanced 33-bus network. Moreover, the ENS results of balanced and unbalanced 33-bus networks indicate lower reliability (higher ENS value) in the unbalanced distribution network. In addition, according to Table 12, it can be seen that the value of power losses using the AEFAPS method for the unbalanced 33-bus network (unbalanced) was 142.85 kW, while for the balanced 33-bus network (balanced) it was 139.55 kW. A comparison of the results according to Table 11 shows that AEFAPS has better capability than the RGA [18], AACO [35], MINLP [36], and CSA [37] methods in achieving lower power losses for the balanced 33-bus network. The results also show that network unbalances increase network losses compared to those with a balanced load. The results comparison of the proposed method with previous studies has confirmed the superiority of AEFAPS in solving the reconfiguration problem with desirable capability.

Conclusions
In this study, many-criteria reconfiguration of unbalanced networks was performed to minimize the losses and enhance power quality indices and ENS minimization using the AEFAPS method considering different loadings. A hybrid meta-heuristic optimization algorithm named AEFAPS was presented to determine the open switches of the network and the optimal configuration of the distribution network. The proposed methodology was implemented on a 33-bus IEEE unbalanced network and a real unbalanced 13-bus network. The reconfiguration results, including losses, voltage sag, and unbalance, and the ENS are presented without and with the reconfiguration based on single-criterion and many-criteria reconfiguration approaches. The results highlighted the better performance of the proposed AEFAPS in achieving a lower (better) value of different objectives than the conventional AEFA, PSO, and GWO methods based on the many-criteria reconfiguration. AEFAPS obtained optimal network configuration with open switches of 7, 11, 13, 28, and 32. In this condition, the losses, voltage sag, voltage unbalance, and ENS were determined at 146.96 kW, 0.653 p.u, 2.42%, and 6.69 MWh/year, respectively. Furthermore, the results showed that the combination of the pattern search with the AEFA achieved a better fuzzy fitness than the traditional AEFA in solving the reconfiguration problem in single-criterion and many-criteria reconfiguration approaches while achieving better values of different objectives. The results confirmed that the many-criteria decision-making optimization-based reconfiguration, in contrast to the single-criterion reconfiguration, provided satisfactory results by compromising different parts of the objective function, and each objective was within its allowable range. The results also showed that the network demand increasing caused losses, increasing and weakening of the voltage sag and unbalance, and also lower reliability of the network, and vice versa. The superior capability of AEFAPS was proved compared with the AEFA, PSO, and GWO methods with lower objective function value, lower convergence tolerance, and more convergence accuracy when compared with previous studies. A reconfiguration of the distribution network integrated with the planning of renewable energy resources to improve the reliability and power quality indices is suggested for future work. Data Availability Statement: Data sharing is not applicable to this article as no new data were created or analyzed in this study.

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

P loss
Power loss R k Ohmic resistance of the k-th line