Enhanced Particle Swarm Optimization-Based Feeder Reconfiguration Considering Uncertain Large Photovoltaic Powers and Demands

The Kyoto protocol recommended that industrialized countries limit their green gas emissions in 2012 to 5.2% below 1990 levels. Photovoltaic (PV) arrays provide clear and sustainable renewable energy to electric power systems. Solar PV arrays can be installed in distribution systems of rural and urban areas, as opposed to wind-turbine generators, which cause noise in surrounding environments. However, a large PV array (several MW) may incur several operation problems, for example, low power quality and reverse power. This work presents a novel method to reconfigure the distribution feeders in order to prevent the injection of reverse power into a substation connected to the transmission level. Moreover, a two-stage algorithm is developed, in which the uncertain bus loads and PV powers are clustered by fuzzy-c-means to gain representative scenarios; optimal reconfiguration is then achieved by a novel mean-variance-based particle swarm optimization. The system loss is minimized while the operational constraints, including reverse power and voltage variation, are satisfied due to the optimal feeder reconfiguration. Simulation results obtained from a 70-bus distribution system with 4 large PV arrays validate the proposed method.


Introduction
The Kyoto protocol mandated that industrialized countries should reduce their greenhouse gas emissions in 2012 to 5.2% below their 1990 levels [1].The Copenhagen Accord later called for the mitigation of global emissions in order to restrict rises in global temperatures to a maximum of 2 ∘ C [2].To attain a low carbon environment, utilizing renewable energy is essential to provide electric power in power systems.Of all renewable energies, wind-turbine generators (WTG), and solar photovoltaic (PV) can provide the cleanest energy to consumers.Although a single unit size of WTG is generally much larger than that of a PV array, the former requires additional effort in terms of installation and maintenance.WTG is also limited by its produced noise, possibly damaging the ecosystem by altering the migration patterns of birds or other animals.However, solar PV arrays can provide noise-free energy, which may be installed in rural fields and rooftops of houses in urban areas.Solar PV has been widely studied in a power system, including solar power prediction [3], capacity sizing [4], and maximum power point tracking [5,6].These works focus only on solar PV modules without addressing system operation problems, which the penetration of solar PV energies in a distribution system causes (e.g., reverse power and high voltages at downstream customers).
Designed as interconnected networks, electric distribution networks are operated in a radial structure.The feeder configuration problem attempts to find a radial operating structure that optimizes network performance while satisfying operation constraints.In practice, distribution system operators may tend to reduce the system real-power loss and enhance voltage quality under normal conditions [7][8][9][10][11][12].Jeon and Kim integrated tabu search with the simulated annealing algorithm to minimize system losses in [7].The tabu search 2 International Journal of Photoenergy attempts to determine a solution using the greatest-descentlike algorithm which could not ensure any convergence property [7].Lin et al. presented a refined genetic algorithm (RGA) to mitigate system losses [8].Traditional crossover and mutation strategies were adjusted through a competition mechanism in RGA.Morton and Mareels presented a bruteforce solution to determine the status of each switch in order to minimize system loss in a radial distribution system [9].The graph theory, which guarantees a global optimum yet requires an exhaustive search, was used in [9].Wu and Tsai presented an approach based on particle swarm optimization (PSO) with an integer code to determine the switch operation schemes, which minimizes the system loss via the feeder reconfiguration in a distribution system [10].Amanulla et al. presented a binary PSO-based algorithm to find the optimal statuses of the switches in order to maximize the reliability and minimize the real-power loss in a power system [11].Hong and Ho proposed a Prufer-number-based genetic algorithm to reconfigure a network to ensure that the voltage drops caused by faults in voltage-sensitive areas are feasible [12].Kavousi-Fard and Niknam proposed a selfadaptive modification method based on the clonal selection algorithm to reduce the cost of active power losses of the network and the customer interruption costs through the feeder reconfiguration [13].Kavousi-Fard and Akbari-Zadeh used the shuffled frog leaping algorithm to reconfigure the network by considering three reliability indices and total system losses [14].Kavousi-Fard et al. used a firefly algorithm to investigate the multiobjective probabilistic distribution feeder reconfiguration problem considering reliability [15].
The widely-installed renewable energy resources in the power system have issued new concerns regarding the feeder reconfiguration.Niknam et al. proposed a probabilistic method considering the uncertainty regarding the active and reactive load forecast errors as well as the wind-turbine output power variations concurrently to study the feeder reconfiguration.Self-adaptive modified teacher learning optimization algorithm was proposed to solve the multiobjective probabilistic problem [16].Niknam et al. assessed the distribution feeder reconfiguration strategy, solved by a selfadaptive bat algorithm, and examined the effect of renewable energy sources on the reliability of the power system [17].Niknam et al. proposed a two-stage approach: the roulette wheel mechanism in conjunction with the Weibull/Gaussian probability distribution functions of wind and demands were implemented to produce deterministic equivalents (scenarios); the feeder reconfiguration problem was then solved by the adaptive modified particle swarm optimization [18].Su et al. presented a binary PSO-based method aiming at enhancing feeder configuration and reactive power control to accommodate more distributed generation resources [19].Malekpour et al. presented an adaptive PSO algorithm to determine the optimal feeder reconfiguration while considering the uncertainty caused from the wind power generation and load demand [20].The point estimate method was used to calculate the stochastic power flow [20].Niknam et al. presented a honey bee mating optimization algorithm to investigate the distribution feeder reconfiguration problem while considering the effect of renewable energy sources on a system performance index [21].
Above literature survey reveals that system loss is a generally accepted performance index in the feeder reconfiguration problem.However, the above studies have the following limitations, which should be addressed in the modern smart grid: (1) renewable generation resources.Although they were not considered in conventional feeder reconfiguration problems [3][4][5][6][7][8][9][10][11][12][13][14][15], they should be considered in modern smart distribution systems; (2) uncertainty in renewable energies.Although uncertainty should be addressed in this area, owing to their intermittent characteristics, Kavousi-Fard and Niknam disregarded this essential factor [13,21].
Although the uncertainty factor in renewable energies was considered in [16,20], the point estimate method sampling two extra unsymmetrical operating points is considered to be less accurate than the Monte-Carlo simulation method [16,20]; (3) coordination between VAR controllers and switches.
Therefore, this work presents a novel method based on the mean-variance concept of particles (i.e., possible solutions) in PSO to determine the optimal status of each switch in the distribution system.Uncertainties in the power generation from solar PV arrays and demands are clustered to become scenarios with different weighting factors.Restated, uncertainties in the problem are transformed into scenarios.The inverters of solar PV arrays are then fully utilized to coordinate with the switches to regulate the voltage profile in order to minimize the system losses using the interior-point algorithm.More specifically, rather than stochastic power flow studies, the regulated voltages of solar PV arrays in each scenario are adjusted by the interior-point-based optimal power flow with fixed switches determined by the enhanced PSO algorithm.
The rest of this paper is organized as follows.Section 2 introduces the assumptions in the studied problem and presents the mathematical formulation of the problem.Section 3 then describes the proposed method based on fuzzy-c-means and enhanced PSO algorithm.Section 4 summarizes the simulation results.Conclusions are finally drawn in Section 5, along with recommendations for future research.

Assumptions and Problem Formulation
2.1.Assumptions.The following assumptions are made for the problem to avoid high voltages at the downstream customers and prevent the reverse power injected into the transmission level at the swing buses (substations).
(a) The status of each tie switch is determined one day ahead by the preventive control using forecasted PV powers and demands with uncertainties.The voltages of inverter-based solar PV arrays are evaluated as supplemental results to support the preventive control.These inverter-controlled voltages should be recomputed in the real-time corrective control.For the preventive control, the scenarios are predetermined based on the empirical experience.No control action for the tie switch is needed once the scenarios occur.
That is, a predetermined network configuration is used for many scenarios.Thus, preventive control is highly applicable to uncertainty caused from the wind powers and demands.
(b) The inverter in either a rural field or an urban area at a bus can be considered as an aggregated one.Restated, hundreds of inverters of the solar PV arrays at a bus are operated with the same setting automatically or according to the commands of operators.Thus, a bus with large solar PV arrays is considered a voltagecontrolled bus rather than a constant PQ bus.
(c) The control center of a distribution system can communicate with the inverters through wire or wireless communication media in the demand management system (DMS) or SCADA.This assumption is applicable owing to the requirement of advanced ICT architecture in the smart distribution system.
(d) The DC voltage of a solar PV array is boosted to a higher DC voltage, which is generally supported by a capacitor.This higher DC voltage is then modulated to a nominal AC voltage via an inverter.To achieve the voltage regulation and prevent the injection of reverse power into the transmission level (swing buses), the DC bus between the boost converter and inverter must be supported by an energy storage; the inverter must be bidirectional as well.Alternatively, a capacitor is still used for the conventionally adopted inverter, and an additional inductor is used for an appended rectifier in parallel.When the PV arrays produce (absorb) the reactive power into (from) the power system, the corresponding bus operates with a lagging (leading) power factor.

Problem Formulation.
As described in Section 1, an operator has many performance indices available to determine the switch statuses.Assume that the operator attempts to minimize the real-power losses and satisfy all operational constraints.In particular, no reverse power is injected into the transmission levels.Further, assume that more than one substation (between the transmission and distribution systems) acts as swing buses in the studied problem.According to the above four assumptions, the studied problem can be formulated as follows: Min  ( Ṽ, θ, P, Q, ) P V      ≥  min , V = 1, 2, . . ., .
The symbols are defined as follows: ( Ṽ, θ, P, Q, ): stochastic real-power loss in the system ( Ṽ, θ, P, Q, ): vector of the stochastic real-power balance equations ( Ṽ, θ, P, Q, ): vector of the stochastic reactive power balance equations Ṽ: vector of stochastic voltage magnitudes at all buses θ: vector of stochastic phase angles at all buses P: vector of stochastic real-power injections at all buses Q: vector of stochastic reactive power injections at all buses : vector of statuses of all switches (1 or 0) Ṽ : stochastic voltage magnitude at bus ,  = 1, . . .,  Ĩℓ : stochastic line flow at line ℓ, ℓ = 1, . . .,  P  : total stochastic power injected into substation ,  = 1, 2, . . .,  P V : stochastic power factor at the Vth photovoltaic array, V = 1, 2, . . .,  : number of system buses : number of system lines : number of substations (swing buses) : number of voltage controlled buses with solar PV arrays.
The values of components in  must ensure that the structure of distribution system is radial, and no islands occur in the system.

Proposed Method
The mathematical formulation shown in ( 1)-( 7) incorporated with the radial and nonislanding constraints is a stochastic mixed binary programming problem.This problem has difficulty in treating stochastic equality constraints in ( 2) and ( 3) and stochastic inequality constraints in ( 4)- (7).Thus, by using the fuzzy-c-means (FCM), the proposed method identifies the best representative scenarios with no random variables.This work presents a novel mean-variancebased PSO algorithm to determine the statuses of switches for all scenarios.Once the status of each switch is given, the nonstochastic representative scenarios can be solved using traditional robust nonlinear programming referred to herein as interior-point algorithm.By using the "fmincon" developed in MATLAB software, this work implements the interior-point algorithm.A detailed description of "fmincon" can be found in [22], which will be neglected herein.

Fuzzy-c-Means to Identify Representative Scenarios. Both
PV power generations and demands are random variables with different stochastic distributions.Theoretically, the Monte-Carlo simulation can obtain an accurate solution in case the number of sampled PV power generations and demands from their stochastic distributions is sufficiently large.However, the Monte-Carlo simulation requires a very long CPU time.If these sampled PV power generations and demands can be clustered, then the center of each cluster can serve as a representative scenario.FCM attempts to identify the representative scenarios obtained from all random variables with their individual stochastic distributions.Restated, the scenarios with real values are generated according to the stochastic distributions of all random variables (PV power generations and demands).The data sets in these scenarios are clustered using FCM.Namely, the center vector in each cluster refers to the representative scenario, and the sum of membership values in each cluster indicates the corresponding weighting factor for this representative scenario.Bezdek et al. defined (, ) as an objective function in the FCM algorithm [23] as follows: where  signifies the clustering number, which is determined in Section 3.2, and is fixed for in FCM.Additionally,  denotes the number of data sets generated from the Monte-Carlo simulation herein;   refers to the vector of the center in the th cluster;   denotes the th data vector for clustering; the membership function value   serves as a weighting factor between   and   .After the minimum of (, ) is gained, the  data vectors are partitioned into  clusters.Let symbols ℎ and  represent the iterative index and convergence tolerance, respectively.Bezdek developed the following four algorithmic steps to attain the optimum.
Step 1. Initialize a matrix of membership functions as follows: where initial ℎ equals 0.

Optimal Number of Clusters.
Yang and Wu proposed a partition separation (PS) validity index, in which a normalized partition coefficient is aggregated with an exponential separation measure for each cluster [24].
The PS validity index for cluster  is defined as follows: where and  denotes the mean of all   s.
The PS validity index for  clusters is defined as follows: A large PS() refers to a situation in which cluster  is compact inside and far away from other (−1) clusters.Consequently, the optimal number  * (i.e.,  defined in Section 3.1) of clusters can be obtained as follows:  the behavior of schooling fish or flocking birds [25,26].PSO comprises a population of particles (individuals) that represent possible solutions moving in a -dimensional search space.The iterative temporary solution and updated value of a particle  are called the position and velocity, respectively.Let vectors   and Δ  be the -dimensional position and velocity of particle ,  = 1, 2, . . .,  (i.e., number of population size), respectively.The proposed method is based on the mean of the particles and standard deviation (squared root of variance) of the distances between any two particles in the th iterations.The updated value of a particle  at iteration  is expressed as follows: where the inertia weight   is decreased linearly from 0.
Since the last term on the right hand side of ( 15) is appended, compared to the original PSO,   = 4 −   −   .Also,    is a -dimensional vector of the mean values of all particles.The last term introduced in ( 15) is used to coordinate the global and local searches.
Let   be the standard deviation of all distances among particles in the th iteration.If all distances follow a Gaussian distribution, then    ± 3  covers 99.7% of all particles.Let   = 3 − 2  .Because   in (15) decreases to zero, the values of   are increased to 3.This work recommends that the particles outside the distance of    ±     in iteration  should be discarded.These discarded particles are substituted by   best .Restated, if    is inside the    ±     , then If    is outside    ±     , then In this work,  +1  denotes the vector of statuses of all switches.All elements of  +1  are bounded to 0 or 1 through the sigmoid function.Because the number of switches with "closed" statuses is fixed (say ) for a distribution system to ensure a radial structure, the largest  elements in  are rounded to 1 and the remaining elements are reset to 0.   After the status of each switch is determined using the mean-variance-based PSO, the topology of a distribution system is fixed.Additionally, the uncertainties in wind powers and demands are modeled by many scenarios with real numbers obtained by FCM.Thus, each scenario becomes a traditional nonlinear "optimal power flow" (OPF) problem, which adjusts the inverter-controlled voltages to minimize the weighting real-power losses and meets the operational constraints, including the reverse power.The weighting factor of the objective function for each scenario  is evaluated by ∑    where the th data set belongs to the th cluster.In this work, the optimal power flow problem is solved using the interior-point algorithm developed in MATLAB software.

Simulation Results
A case study involving a 70-bus distribution system is used as an example [27], as shown in Figure 1, to show the simulation Four large solar PV arrays are installed at buses 13, 26, 44, and 61.Characteristics of the solar PV arrays are the Weibull distributions; that is, pdf (Pg) = (/)(Pg/) − 1 exp(−(Pg/)), where Pg denotes the PV generation.Table 1 lists the installed capacities of all PV arrays and their parameters of the Weibull distributions.
In this system, the loads at 17 buses are modeled by the Gaussian distributions.They are buses 15-17, 19, 31, 33-35, 50-54, and 66-69.The mean and standard deviation of each stochastic load are 0.8 p.u. and 5% (based on its original MVA), respectively.The loads at other buses are fixed.
Three hundred sets of PV generations and loads are randomly produced using the parameters of the Weibull and Gaussian distributions, respectively; that is,  = 300, in this work.Next, by using FCM, these 300 sets of data are clustered together by considering different .Theoretically, the maximum number of clusters is √ 300 (≅17, see ( 14)).However, the computational burden in a PC is considered here; let the maximum number of clusters be 10.The total CPU time ( = 2, 3, . . ., 10) to cluster these sampled PV generations and loads is 1.38 s, which was estimated by a PC with Intel Core i5 2.50 GHz and 8 GB RAM.Table 2 lists the PS validity indices for  = 2, 3, . . ., 10.This table reveals that the optimal clustering number is 8 (i.e.,  * = 8 in ( 14)) because its corresponding PS validity value (3.039) is the largest.
The 8 representative scenarios obtained by FCM were examined first by running OPF.The voltage profile and reactive powers were adjusted using 4 inverter-based solar PV at buses 13, 26, 44, and 61.Table 3 lists the 17 load buses with their representative real-and reactive-power loads (p.u.) in the 8 representative scenarios.The 4 solar PV arrays generate real power only initially with the unity power factor.This table also shows the weighting factor obtained by FCM for each scenario.
First, the interior-point-based OPF was used to regulate the voltage profile and minimize the system losses individually for these 8 scenarios.Table 4 summarizes the optimal results before the feeder reconfiguration.This table reveals that a reverse power −95.092 kW was injected into bus 1 in scenario 7.There is also reverse reactive power (−361.796kVAr) at bus 26.The other scenarios had no reverse power at buses 1 and 70.
By using the proposed method, the statuses of 79 switches were determined using the mean-variance-based PSO; the voltage profile was also regulated using the interior-pointbased OPF.The population size was 10, and the inertia was decreased linearly from 0.5 to 0.3 in 100 iterations in the proposed PSO algorithm.Table 5 lists the optimal solutions obtained by the proposed method.The injection powers −95.092 (bus 1) and 2,800.633kW (bus 70) before feeder reconfiguration were changed to 702.954 (bus 1) and 1,893.876kW (bus 70) after feeder reconfiguration in scenario 7, respectively.The system losses were also reduced from 157.874 to 92.875 kW.The proposed method required 11,835.41CPU seconds to attain the optimal solution.Figure 2 shows the optimal feeder reconfiguration using the proposed method.
Table 6 shows the optimal solutions gained by using chaos-based PSO and GA.The results gained by chaos-based PSO and GA are the same.The numerical parameters used in the chaos-based PSO are the same as those in the proposed method.The crossover and mutation rates in GA are  The aggregating loss attained by the proposed method is 92.875 kW, which is smaller than 93.738 kW obtained by both chaos-based PSO and GA. Figure 3 illustrates the optimal feeder reconfiguration obtained by the chaos-based PSO and GA algorithms.

Conclusions
This work presents a novel method that incorporates FCM, partition separation validity index, mean-variance-based PSO, and interior-point-based optimal power flow to conduct feeder reconfiguration in order to prevent the revere power caused by uncertain solar PV arrays from injecting into the transmission level.The novelties are as follows: (a) the uncertainty in the PV powers and demands is modeled by multiple scenarios gained by FCM in order to avoid probabilistic calculations; (b) the novel mean-variance-based PSO implements an adaptive elite strategy that replaces the particles far away from   best with   best .The number of replacement decreases according to the standard deviation of distances among particles.Simulation results demonstrate that the proposed method is faster than the chaos-based PSO method.Furthermore, the proposed method can solve the problem with the reverse power and high voltage at the customers' side while minimizing the system loss.

Figure 1 :
Figure 1: One-line diagram for the 70-bus distribution system.

Figure 3 :
Figure 3: Optimal feeder reconfiguration obtained by chaos-based PSO and GA.

Table 1 :
Installed capacities of PV arrays and parameters of the Weibull distributions.

Table 2 :
PS validity indices for different numbers of clusters.
) 3.3.Mean-Variance-Based PSO Algorithm to Determine .As an evolutionary optimization method, PSO minimizes an objective function (performance index) by imitating International Journal of Photoenergy 5 5 to 0.3 because a large   is used for global searches, whereas local searches require a small   .The symbols  1 ,  2 , and  3 are random numbers within [0, 1].Additionally,   best and [25,26]e the best position of a particle in the iteration  and the best known position up to present, respectively.The learning factors  and  should meet  +  ≤ 4[25,26]and are defined as

Table 4 :
Optimal results before feeder reconfiguration.

Table 5 :
Optimal results after the feeder reconfiguration using the proposed method.

Table 6 :
Optimal results after feeder reconfiguration using chaos-based PSO and GA.obtained by the proposed method.Buses 1 and 70 are two substations (the swing buses).The solid lines denote the normally closed sectionalized switches, and the dot lines are normally open tie switches.Totally, the statuses of 79 switches must be determined in this work by the proposed meanvariance-based PSO. results