A New Dg Multiobjective Optimization Method Based on an Improved Evolutionary Algorithm

A distribution generation (DG) multiobjective optimization method based on an improved Pareto evolutionary algorithm is investigated in this paper. The improved Pareto evolutionary algorithm, which introduces a penalty factor in the objective function constraints, uses an adaptive crossover and a mutation operator in the evolutionary process and combines a simulated annealing iterative process. The proposed algorithm is utilized to the optimize DG injection models to maximize DG utilization while minimizing system loss and environmental pollution. A revised IEEE 33-bus system with multiple DG units was used to test the multiobjective optimization algorithm in a distribution power system. The proposed algorithm was implemented and compared with the strength Pareto evolutionary algorithm 2 (SPEA2), a particle swarm optimization (PSO) algorithm, and nondominated sorting genetic algorithm II (NGSA-II). The comparison of the results demonstrates the validity and practicality of utilizing DG units in terms of economic dispatch and optimal operation in a distribution power system.


Introduction
With the increasing demand for clean and renewable energy, the issue of distribution generation (DG) is drawing more attention worldwide.DG provides voltage support to largescale distribution power systems, which results in reliability improvements and reduction loss in the power grid.DG technology has become a hot research topic, given the increasing global concerns about environmental protection, energy conservation, and the increasing sophistication of wind power, photovoltaic power generation, and other renewable energy technologies.After DG is connected to a distribution network, the distribution network's structure, operation and control mode will tremendously change, and the distribution system automation and the demand-side management must consider the coordination between DG and distribution network control.Deciding the optimal DG output is a challenging research problem, especially considering the multiple optimal objectives associated with cases of multiple DG unit injections.
Traditionally, multi-objective DG optimization has been treated as a single-objective optimization problem using suitable weighting factors to form a weighted sum of single objectives.This approach has the disadvantage of finding only a single solution that does not express the tradeoffs with different weighting factors.Generating multiple solutions using this approach requires several runs with different factors, which leads to long running times [1,2].Recent studies have treated the objectives simultaneously and independently as a true multi-objective optimization problem.However, the optimization problem becomes more complicated due to such issues as continuity, local optima, linearization, and so forth.New optimization techniques, such as particle swarm optimization (PSO), different evolution (DE), and evolutionary programming (EP), have recently been introduced and applied in the field of power systems and with promising results [3][4][5][6][7][8][9][10][11][12].In a recent study [3], a differential evolution approach was proposed to solve an optimal power flow problem with multiple objectives.The active power dispatch and reactive power dispatch were considered.A nonlinear constrained multi-objective optimization problem was formulated.A general overview of evolutionary multi-objective optimization was provided in [4], and the most representative algorithms were discussed.In [5], Pareto-based multi-objective evolutionary algorithms were discussed and evaluated.A nondominated sorting genetic algorithm (NSGA), a niched Pareto genetic algorithm, and a strength Pareto evolutionary algorithm (SPEA) were developed and applied to an environmental/economic electric power dispatch problem.A multi-objective formulation for sitting and sizing of DG units was proposed in [6].The method involved searching for a compromise between the cost of network upgrades, cost of power losses, cost of energy not supplied, and cost of energy required by the served customers.A genetic algorithm was implemented to obtain a noninferior solution set.However, their method cannot guarantee a solution to be optimal solution.In [7], an improved swarm optimization (IPSO) method was presented to solve the multi-objective optimal power flow problem.The multiobjective optimal power flow considered the cost, loss, voltage stability and emission impacts as the objective functions.A fuzzy decision-based mechanism is used to select the best compromise Pareto set solution obtained by the proposed algorithm.In [8], a new penalty parameter-less constrainthandling scheme was employed to improve the performance of the evolutionary algorithm.The experiments in that paper revealed that PSO performs better in terms of solution quality and consistency, and DE performs better in terms of mean computation time.An improved Cai and Wang's method has been proposed to combine multi-objective optimization with differential evolution to address constrained optimization problems in [9].The method provided a novel infeasible solution replacement mechanism for differential evolution in theory.In [10], a robust DE algorithm was proposed for the control of selective harmonic distortion and total harmonic distortion.A fuzzy optimization technique and DE optimization method are described.
The literature includes several DG output studies that examined multiple objectives and applied evolutionary optimization techniques.From the perspective of mathematical optimization, DG unit injection is a complex multiobjective optimization problem that presents a challenge to the optimization analysis of a distribution power system.The objectives include optimal energy consumption, the minimum power consumer's electricity purchasing cost, and the minimum power loss based on the constraints of power grid security and DG power output.Multi-objective economic/emission dispatch algorithms were investigated in [11,12].In the optimization methods literature, the simulated annealing technique has been applied to optimize the proposed multi-objective model of DG planning [13].The multiobjective problems were solved by converting the original model into an equivalent model through calibration of the weighted factors method.In [14], a multi-objective Tabu search-(TS-) based method was utilized to optimize a DG allocation problem.In that paper, the TS-based approach was provided to find the optimal Pareto set.Fuzzy optimization was also used to solve the multi-objective optimization of DG allocation in [15].Voltage drop reduction, short circuit capacity augmentation, decrease operation cost, and system loss reduction were considered objectives for formulating fuzzy optimization.
In this paper, a DG multi-objective optimization method based on an improved evolutionary algorithm was investigated for a distribution power system.Adaptive crossover and a mutation operator were used in the evolutionary process, and simulated annealing was combined in the iterative process.A fuzzy clustering algorithm was applied to manage the size of the Pareto set.The rest of the paper is organized as follows.In Section 2, the formulation of the DG multiple-objective optimization for distribution management is presented.The Pareto-based algorithm and some basic concepts are introduced in Section 3. The improved Pareto evolution algorithm is described in Section 4. Section 5 provides the numerical results and comparison analysis with the proposed approach using the revised IEEE 33-bus system.The conclusion and future work are provided in Section 6.

Problem Formulation
2.1.Objective Functions.Three objectives are considered in the optimization model, which includes the fuel cost and the pollutant emission penalty, reducing consumer costs on electricity bills when DG units are injected into the distribution network and reducing transmission line losses.The first optimization objective is minimum energy consumption and a pollutant emission model, which is mainly based on government requirements.There will be more penalties if the system emits more pollutants and exhibits greater fuel consumption.The second objective is consumer related, where the consumer uses DG to maximize savings on their bills.The third objective is to lower system line losses, which is the demand objective of the power supply provider.The three objectives involve perspectives based on government requirements, consumer needs, and power supply enterprise needs, and the objectives can conflict.For example, when a consumer utilizes a micro gas turbine to maximize their savings on their energy bill, there is a subsequent increase in fuel cost and pollutant emission.In addition, the extra power from the micro gas turbine will increase or decrease the line losses, depending on the size and placement location of the micro gas turbine.

Fuel Cost and Pollutant Emission Minimization.
The first objective is to minimize the fuel cost and the pollutant emission penalty, which reflects the impact of energy utilization on the environment.It can be expressed as follows: where   is the energy consumption cost and   is the pollutant emission penalty.The fuel cost   normally can be further expressed as follows: Note that   ,   , and   are the quadratic cost coefficients of the th DG, and  DG is the number of distributed generators.
DG  is the real power output of the th generator. DG is the vector of real power outputs of generators and defined as follows: The pollutant emission quantity can be obtained based on DG output.Then, based on the penalty standard, the environmental penalty for pollutant emission is calculated as follows: where   is pollutant 's emission quantity and   is the penalty standard of pollutant .

Maximization of Cost Savings Using DG.
The second objective is to maximize the cost savings on electricity user bills when the DG is injected into the distribution network.The savings in electricity, which should have been purchased from the power supply enterprise, are the total power output of the DG units.Utilizing DG output and time-of-use (TOU) rate, consumer electricity purchasing costs could be reduced as follows: where   1 is peak price from  1 to  2 ,   2 is off-peak price, and  DG  is the DG total power output at moment .

Minimization of Line Losses.
The third objective is to minimize the system line losses after DG injection into the distribution network.This objective can be expressed as follows: where [] is the active power, [] is the reactive power at branch , [] is the voltage at branch  after DG injection, [] is the resistance of branch , and  is the number of branches in the distribution network.
In the previous three optimization models, the fuel cost and the pollutant emission penalty function  1 () and the system loss function  3 () should be minimized, whereas the cost-saving function  2 () should be maximized.

Constraints.
Three constraint conditions are considered in the optimization model, which includes constraints of power flow equations, nodal voltage, and DG capacity.

Equality Constraints. The constraint of power flow equations is described as follows:
where  DG  and  DG  are active and reactive generation outputs, whereas    and    are the active and reactive load at bus , respectively,   and   are the transfer conductance and susceptance between bus  and , respectively, and  is the number of buses.

Inequality Constraints
Generation limits: Load bus voltage constraints: Thermal limits: In the inequality constraints,  min DG  ,  max DG  ,  min DG  , and  max DG  are the lower/upper active and reactive generating unit limits of DG, respectively. max  is the apparent power thermal limit of the circuit between buses  and .
There is always a limit on penetration of DG for a distribution power system to ensure reliability.Different countries have different penetration factor values.The penetration factor indicates the aggregated DG rating on an electric power system (EPS) feeder, divided by the peak EPS feeder load.If we assume that the maximum DG penetration factor is 25%, then the maximum injected DG capacity should be limited to 25% of the maximum total load in the distribution network, which can be described as follows: where  DG  is the DG access capacity at node  and  max is the maximum load capacity of distribution network.

Overview Formulation.
Aggregating the objectives and constraints, the problem can be formulated as a nonlinear programming problem as follows: minimize [ 1 () ,  2 () , . . .,   ()] , subject to ℎ  () = 0,  = 1, . . ., ,   () ≤ 0,  = 1, . . ., , where  is the number of objectives and  is the vector of dependent variables consisting of slack bus power output and DG active power out  DG  , load bus voltage   , and generator reactive power outputs  DG  .Thus,  can be expresses as follows: where  is the number of inequality constraints,  is number of equation constraints,  is the number of objectives, and  is the number of DG units.

Basic Pareto-Based Evolutionary Algorithm.
The traditional Pareto-based evolutionary algorithm is shown in Figure 1.The detailed algorithm procedure is explained in [16].The main improvements on the Pareto-based algorithm can be generalized as follows.The penalty function is established to constrain the solution of the objective function.
The adaptive crossover and mutation are adopted in the evolution process, which improves the probability of global optimization.The simulated annealing algorithm is added to the iterative process, so that the algorithm is able to seek the optimal solution globally and rapidly converges to the optimal solution.

Proposed Improved Pareto Evolutionary Algorithm
where   is the population size and  is the number of control variables.Each variable  in a solution vector  in the generation  initialized within its boundaries   min and   max .

Fitness Evaluation.
The objective of each solution  1 (),  2 (), . . .,   () will be computed.The individual fitness values in both the population-based set  and nondominated archive set  will be calculated based on (17).The mismatch of each constraint value is multiplied by a large value and added to all objectives to remove infeasible solutions.The methodology is to evaluate the feasible solutions according to the value of objective function and remove the infeasible solutions according to the constraints.The individual's fitness () will be obtained from the sum of the primary fitness value () and the density () as follows: where () = ∑ ∈+,≻ (), () is the objective evaluation for the individual .(≻ indicates a dominated relation,   ≻   indicates   dominates   ,   is nondominated, and   is dominated).
() = 1/(   + 2),  = √  + .   represents the objective space distance between individual  in  and the th nearest neighbor individual in the .It is the Knearest neighbor (KNN) method, and the distance between the individual  in  and the other individuals in the  need to be computed, and then, the distance value can be sorted.

Adaptive Crossover and Mutation Probability.
The selection of crossover probability   and mutation probability   dominates the solution process.  and   determine the generation speed and the probability of new individuals, respectively.If   exceeds the threshold, the generation speed of the new population will be quicker, which means that there is a stronger capability to explore new space.If   is extremely small, the search process will be quite slow.If   is too large, the search process will be more random.The adaptive value of   and   is obtained from the following evaluated equations: where  avg is the average fitness value,  max is the biggest fitness value,   is the bigger fitness value of two genes in the crossover process,  is the mutating individual's fitness value, and the constants

Pareto Optimal Selection
Using Fuzzy Set Theory.In this paper, fuzzy set theory is used to select the optimal solution set among the obtained multiobjective solution sets.Fuzzy sets are sets whose elements have degrees of membership.Fuzzy set theory permits the gradual assessment of the membership of elements in a set.This membership is described with the aid of a membership function valued in the real unit interval [0, 1].First, define a linear membership function   as the weight of target  in a solution: where  max  is the maximum of th objective function,  min  is the minimum of th objective function, and   is the solution of th objective.The previous equation provides a measure of the degree of satisfaction for each objective function for a particular solution.
The dominant function   for each nondominated solution  in Pareto solution set is calculated as follows: where  is the number of the Pareto solution set and  0 is the number of the optimization objectives.
Because the value of   determines the capability of the solution, the solution with maximum   will be Pareto optimal.Moreover, the feasible priority sequence can be obtained by the value of   , in descending order.
The best Pareto optimal solution is the one achieving the maximum membership function   , as shown in (20).

Simulated Annealing in Population-Based Individual
Selection.Simulated annealing (SA) is a generic probabilistic metaheuristic for the global optimization problem of locating a good approximation to the global optimum of a given function in a large search space.It is often used when the search space is discrete.Here, SA is utilized in the individual selection.
Based on the individuals after selection, crossover and mutation steps, the simulated annealing operation is performed on the individuals of the population.The two genes in each individual will be selected and disturbed randomly.Then, the new individual will be evaluated to form new fitness values.If the fitness value of a new individual is larger than the old value, then the old individual will be replaced by the new individual.If the fitness value of the new individual is smaller than the old value, the new individual can also be accepted using the following probability: where  +1 and   are the fitness of the new individual and old individual, respectively, ( +1 ) is the acceptance probability at  +1 temperature, and  is the temperature descending coefficient.

Convergence Condition.
The iterative procedure can be terminated when any of the following conditions are met: (1) the true Pareto front is obtained, and (2) the iteration number of the algorithm reaches the predefined maximum number of iterations.However, the true Pareto front will not be known in advance in most practical multiobjective problems, so the convergence condition is to iterate to a predefined maximal iteration number.

Flowchart of Proposed Algorithm.
The flow chart of the proposed algorithm is illustrated in Figure 2. As shown in Figure 2, the steps of the proposed evolutionary algorithm are described as follows.
Step 1. Generate an initial set  randomly and an empty archive set  over the problem space; initialize the parameters of the population size , nondominated archive size , and maximum generation's number .
Step 2. Establish the penalty function to constrain each objective function, and then form new objective functions.
Step 3. Compute the fitness of individual in both the population-based set  and the nondominated archive set NDSet.The objective of each solution  1 (),  2 (), . . .,   () will be computed.
Step Step 6. Copy the superior dominated individual to  .
Step 7. Evaluate the convergence criteria.If the iteration number  ≥ , terminate the iteration to obtain the Pareto optimal solution and output the best solution; otherwise, set  =  + 1, and continue to Step 8.
Step 8. Perform adaptive crossover and mutation operation on the individuals of .
Step 9. Perform a simulated annealing operation, and then go to Step 3.

Experiments and Results
To demonstrate the effectiveness of the proposed method, the algorithm in Section 4 was implemented to obtain solutions for optimal active power dispatch of DG.The IEEE 33 bus distribution system was examined, and three objectives were considered in this study.These objectives were fuel cost/pollutant emission, transmission line loss, and cost savings on bills using DG.Photovoltaic (PV) panels, diesel turbine, and wind turbine distribution are injected into bus 7, bus 17, bus 21, and bus 32, respectively, as shown in Figure 3.In this paper,   and   are defined as follows:   1 = 0.4,   2 = 0.3, and   3 = 0.2 and   1 = 0.2,   2 = 0.1, and   3 = 0.05.The maximum iteration number is set to 200.The proposed algorithm was coded in C++ and run on an Intel i5-3210M 2.5 GHz notebook with 4 GB RAM.

Energy Utilization
Cost.Among the four DG units, only the two diesel turbine DG units have fuel cost.Because it would be difficult for market players to accept/implement  a central cost-based dispatch in the distribution system including DG units, the cost of fossil-fuel consumed by micro diesel turbine is calculated as follows: where   is fuel price at power unit  and (   ) is the required fuel quantity for power unit  at the moment .

Penalty on Pollutant Emission.
As global environmental pollution is growing, optimizing power generation and pollutant emission costs are two conflicting goals.These goals present a restrictive and coordinated relationship.Environmental cost mainly refers to the fines related to pollutant emission.Tables 1 and 2 show the pollutant emission data for various DG units and the standard electric power industry pollution fines, respectively.In reference report [19], there are similar pollutant cost coefficients for distributed generation.Based on the DG unit output  following multiple-objective optimization, the quantity of pollutant emission can be obtained.

Optimization Results
. Using the optimization model developed in Section 3, the optimized output of four DG units over 24 h and the power system losses after DG unit injection are shown in Figures 4 and 5, respectively.As shown in Figure 4, the four DG units have different active power outputs at different time periods in a day.The diesel power output will increase when the solar and wind power outputs are at a low level.When the PV output and wind power output increase to the peak, it will stop increasing and stay at the peak power output, and then, the diesel power output will gradually decrease.As shown in Figure 5, the line losses greatly decrease DG unit penetration into the distribution system.From the hours 8 to 17, the total output of the four DG units provides enough active power, which improves the voltage quality and reduces the line loss.
The forecasted and optimized solar power outputs based on the computed results are shown in Figure 6, and the forecasting and optimized wind power output are shown in Figure 7.The forecasted PV and wind generation values are based on historical distribution system data.Because of the cooperative optimization, the optimal real power values of PV and wind generation are smaller than the forecasted values in the peak time period.
Assuming that the coal consumption from the power plant is 0.35 kg/kWh and the highest coal price is 0.124 $/kg, the cost savings for coal consumption by using clean energy is illustrated in Figure 8, which shows that these increases in solar and wind power output results in greater coal consumption cost savings.A pollutant emission penalty reduction curve was obtained based on data from Tables 1 and 2, and the hourly penalty reduction for pollutant emission is shown in Figure 9.As there is no pollutant emission for solar and wind power generation, when the output of new energy power supply increases, the environment cost will decrease significantly.
Assuming that the time-of-use price is 0.095 $/kWh for peak time from 6:00 am to 22:00 pm and 0.054 $/kWh in other period, the cost saving for the electrical bills of users per hour is shown in Figure 10.Because the price is at a high level from 6:00 am to 18:00 pm, the bill saving increases with the increase of PV and wind power output.The results from the case study demonstrates that the system loss is greatly reduced by 65%, so that the users, in total, can save $1,671 per day on their electricity bills, and power plants can save   The forecast of wind generation The optimal computation of wind generation $870 and $9,906 on their coal costs and pollutant emission penalties per day, respectively.

Comparison of Different Algorithms.
The proposed algorithm was compared with the SPEA2 [17] and the particle swarm optimization method [20] and NSGA-II [21,22].The IEEE 33-bus system with four DG units was utilized as an example for this comparison.The load data at hour 11 is selected as the basic load data.The convergence condition was that the iteration number exceeded the preset maximum iteration number, which was set to 200.In PSO, the cognitive ratio and social ratio are all equal to 2.0.The number of swarm particles is 100.In NSGA-II, the crossover ratio is set to 0.8, and the mutation ratio is set to 0.2.The size of population is set to 100.Table 3 shows that the proposed algorithm performs better than the SPEA2 and the PSO algorithms with respect to calculating the multiple objective objectives in the same limited iterations, and the proposed algorithm has better convergence speed than SPEA2 and PSO because simulated annealing is added in addition to the adaptive crossover and mutation operations.Compared with NSGA-II, the proposed algorithm has approximate speed in searching the Pareto front.

Conclusion
This paper presented an improved Pareto-based evolutionary algorithm, which increases the global optimization ability with a simulated annealing iterative process and fuzzy set theory, to solve the multiobjective optimization problem for a distribution power system.The proposed algorithm was utilized to optimize a model of DG unit injection with objectives of maximizing the utilization of DG while minimizing the system loss and environmental pollution.The results indicate that the proposed optimization is applicable to practical multiobjective optimization problems that take

Figure 3 :
Figure 3: Schematic diagram of the IEEE 33-bus system with 4 DG units.

Figure 4 :Figure 5 :
Photovoltaic generation Wind power generationFigure 4: The optimized output of four DG units in 24 hours.
PV generationThe optimal computation of PV generation

Figure 6 :
Figure 6: The comparison of forecasted PV value and the computed optimal PV value.

Figure 7 :
Figure 7: The comparison of forecasted wind values and the computed optimal value.
[17,18]osed to solve the multiobjective optimization problem based on SPEA2[17,18].4.2.Initialization.In the improved SPEA2, an individual  at generation  is a multidimensional vector    = ( ,1 , . . .,  , ).The population is initialized by randomly generating individuals as 4.1.Overview.To solve the difficulties in traditional optimization techniques, a new evolutionary population-based searching technique 4. Duplicate the nondominated individuals in both the population and nondominated archive set to a new archive set  , if the size of   exceeds , then reduce   by means of the truncation operator; otherwise, fill   with dominated individuals in  and .Evaluate if the nondominated set   exceeds the predefined size .If the size of   is larger than , then truncate the nondominated individuals; otherwise, continue to Step 6.

Table 1 :
The pollutant emission rate of DG units (g/kWh).