A Multi-Objective Hybrid Genetic Algorithm for Sizing and Siting of Renewable Distributed Generation

Renewable generation has been addressed in several aspects but it still represents a new paradigm for the expansion of the electricity supply. This paper aims to propose a new model for the sizing and siting problem of distributed generation (DG), based on renewable sources and considering three main aspects: technical, from the distribution utility viewpoint; economical, from the DG owner’s viewpoint, and environmental, from a sustainability perspective. A multi-objective Genetic Algorithm and the Maximin metric are implemented to obtain optimal Pareto sets; also, three decision criteria, considering the concept of preference, are applied to select a final solution from Pareto sets. Case-studies are carried out in medium voltage systems: the 69-bus distribution test system, known from literature, and a 918-bus Brazilian distribution system. Diversity of alternatives in the obtained Pareto sets testify algorithm effectiveness in searching for solutions to the distributed generation sizing and siting problem, in order to ensure power loss reductions, investment return, and environmental benefits. The proposed methodology contributes to the discussions and perspectives among electricity utilities, DG owners, society, and regulators regarding planning and decision


Introduction
Electrical power systems have been characterized by a centralized production, in which generation location is usually far from load centers, therefore requiring extensive and expensive transmission and distribution networks [1]. Currently, new opportunities provided by technological advances in reduced scale energy generation (on-site power generation) and regulatory changes have made possible the growth of the DG concept [2], leading to new challenges on power systems planning and operation due to the increased complexity of the systems resulting from DG insertions.
DG plants intrinsically bring some advantages to the electricity sector, such as load relief, reduction of technical losses, diversification of energy sources, stimulation of new technologies and a higher local equalization between energy offering and demand [3]. In addition, the DG concept represents new investment perspectives to the electricity market, foments the creation of new companies and originates new jobs and careers; furthermore, its growing represents a step to transform a passive customer into an active one.
Nonetheless, Distribution Companies (DisCo) have faced many challenges due to DG penetration, such as bidirectional power flows, rising of short-circuit levels, overvoltages, and overcurrents, islanding operation, and an increase in the operating scenarios caused by intermittence of some primary energy sources, e.g., sun, wind, and water [4]. As a result, the academic community has shown interest in these topics, dealing with them through a.
Three objective functions are considered regarding the proposed multi-objective model: energy losses as a technical aspect; DG plant installation costs as an economic aspect; and equivalent CO 2 emission as an environmental concern of society. To the best of the authors knowledge these objective functions have not been simultaneously considered in the way they are approached in this paper. b.
A new multi-objective solution methodology based on Genetic Algorithm combined with Maximin metric, named in this work as Multi-Objective Hybrid Genetic Algorithm (MOHGA), are implemented to obtain optimal Pareto sets of solution alternatives. c.
Max-Min approximation as non-preference criterion and Minimal power losses and maximum net present value as preference criteria are applied as decision making strategies to select final solutions from obtained optimal Pareto sets. To address the proposed modeling for DG sizing and siting problem and the adopted solution process, the following assumptions are made in this work: • Only renewable energy sources are considered (hydraulic, wind, and photovoltaic). • Network connection buses and maximum capacity (divided in modules) of generators are known data of the problem, in order to take into account proper conditions such as available source of energy, environmental restrictions, and occupation area. • Maximum capacity for DG plants is set in 5 MW, since it covers common regulations (specifically, in Brazil, this value corresponds to the maximum allowed for DG's).
This paper is organized as follows: the proposed multi-objective optimization model and solution method are addressed in Section 2; Section 3 presents the results and analysis performed in the two networks as case studies, whose main conclusions are drawn in Section 4.

Modeling and Solution Strategies
This section presents the modeling of the proposed hybrid genetic algorithm. The nomenclature can be consulted in Appendix A.
In general terms, a multi-objective optimization problem consists of solving two or more usually conflicting objective functions, in which the decision variables must satisfy a set of equality, inequality and border constraints, as summarized in (1): where n f -dimensional vector F is the set of objective-functions, n x -dimensional vector x corresponds to the decision variables of the optimization problem, g g and h h functions are the equality and inequality constraint functions, respectively; finally, x min and x max are related to lower and upper boundaries, respectively, of each decision variable.

Multi-Objective Optimization Model
The DG sizing and siting problem addressed in this work assumes that the locations of accessible primary energy sources and the maximum generation capacity available for DG projects are previously known. Each DG plant n for a total number of N plants can install NDG n modules for a total number of NDG max n modules, each one with capacity P DG n , such that the power output of each DG plant n is given by NDG n .P DG n . Therefore, decision variable vector x of the optimization problem is composed by network buses and number of modules of each DG plant n, with total number 2 · N (size and site) of variables.
The vector of multi-objective functions (F) of the proposed modeling is formed by three objective functions ( f 1 , f 2, f 3 ). The first objective-function ( f 1 ) represents an economical viewpoint and corresponds to minimize the total installation cost of DG plants (C instal ), obtained by Equation (2), where CO n instal is the investment cost of the n-th DG plant, NDG n is the quantity of generation modules of the n-th DG plant to be installed and N is the total number of DG plants to be connected in the network. In this paper, installation cost data are mean overnight cost: 1497 USD/kWe for onshore wind, 1210 USD/kWe for photovoltaic plants, and 1492 USD/kWe for hydro [25]. The second objective function ( f 2 ) corresponds to the minimization of daily active power losses (LAE daily in MWh), which is one of the objectives most commonly considered in DG sizing and siting planning [26]. The required data to calculate branch currents are obtained by a backward/forward sweep load flow known by Power Summation Method (PSM) and the daily energy losses are calculated by Equation (3), where r k−m and I k−m are the resistance and current of the k-m branch, respectively; N T is the number of time intervals considered in the daily analysis, Ω L is the set of distribution network branches, and T p is the duration of the time interval. In this case, time intervals of one hour are considered, aiming at a better discretization of daily consumption, although other time intervals can be considered in the modeling without significant changes in results.
The third objective function ( f 3 ) reflects environmental concerns and corresponds to the minimization of daily amount of equivalent CO 2 emission EqCO 2 E, given by Equation (4). CO 2 emission factors EqCO 2 EF are based on [27][28][29][30]: 21 gCO 2 eq/kWh for onshore wind, 17 gCO 2 eq/kWh for photovoltaic and 105 gCO 2 eq/kWh for hydro, while adopted capacity factors cf are [26]: 0.18 for photovoltaic, 0.34 for onshore wind and 0.47 for hydro. In this case, c f n is the capacity factor of the n-th DG plant and P DG n is the active power capacity for each module of the n-th DG plant.
It should be noted that a greater number of DG modules results in higher installation costs and CO 2 emissions; nonetheless, dispersed generation in distribution systems represents load reductions and, therefore, provides power loss reductions. Thus, minimization of installation costs and CO 2 emissions are conflicting with minimization of power losses.
Daily generation curves for each primary energy source (solar, hydro, and wind) and daily load curves for the three consumer classes (residential, commercial, and industrial) assumed in this work for power flow simulations are presented in Figure 1. These curves are based on Brazilian DG plants and load characteristics and can be changed according to local characteristics, although generation curves were based on the capacity factors. Overvoltage and overcurrent constraints are considered in the proposed model as they represent network limitations, which are generally set by local regulations. Concerning these constraints, the maximum voltage level ( ) allowed is 1.05 p.u. and overcurrent is identified using the network section ampacity (Ak-m). Another constraint considered in the model is related to the fast voltage variation due to renewable DG Equality constraints (g) are referred to the power balance in system nodes, as shown in [27,29,31], and defined in this work as PSM equations (x, p). Inequality constraints (h) usually represent inherent distribution systems characteristics or regulatory restrictions, such as voltage level constraints, cable capacity, and limits on long-time and short-time voltage variations. Due to the large number of constraints that can be modeled, three operational constraints were chosen: overvoltage, overcurrent, and penetration index.
Overvoltage and overcurrent constraints are considered in the proposed model as they represent network limitations, which are generally set by local regulations. Concerning these constraints, the maximum voltage level (V max ) allowed is 1.05 p.u. and overcurrent is identified using the network section ampacity (A k−m ). Another constraint considered in the model is related to the fast voltage variation due to renewable DG intermittence, such is calculated every hourly power flow simulation and then a ∆V max equal to 0.05 p.u. is defined to verify the compliance with referred constraint.
The penetration index (PI) represents a measure of impact (load relief) of DG plants on electrical systems. Various works found in the technical literature try to find or simply suggest a DG penetration value able to provide benefits to DisCo and DG owners without compromising the network operation. This index is calculated as shown in Equation (5) and a maximum penetration index is defined in each case study in this work.
Therefore, the complete mathematical model is presented in Equation (6).
For each candidate solution, two types of penalization for constraint violations are adopted. When the PI is violated, the number of generation modules NGD n of an n-th DG randomly chosen is decreased until the PI is attended. When any of the other constraints are violated, the objective function corresponding to the total installation cost is increased by an additional connection cost C con , aiming to simulate real situations where network reinforcements will be necessary. In this work, the connection cost is modeled as an exponential function multiplied by the installation cost, as given in Equation (7).

Multi-Objective Solution Methodology
The proposed multi-objective solution methodology (MOHGA) consists of two steps: a Genetic Algorithm routine to perform the optimization process in searching of feasible solutions (population in GA terminology) and the Maximin metric routine for dominance measure and classification of the solutions in population to form the Pareto-set.
For each i-th iteration of the optimization process, a feasible solution X i represented by variable vector x (called individual in the GA terminology) is evaluated and a vector of objective-functions F(X i ) is obtained. However, interpretations of the Maximin metric depends heavily on the values of the objective functions, which have different units and magnitudes. To mitigate this problem, a normalization of vector F(X i ) shown in Equation (8) is proposed, where F(X 0 ) corresponds to the vector of objective functions evaluated in the base-individual X 0 referring to the solution of mono-objective problem of the system loss minimization, so that no null value to any of the objectives can be obtained.
Maximin metric provides two useful informations [27,32]: (1) negative value indicates a non-dominated solution; (2) the more negative is the value, the less crowded is the region of Pareto frontier the solution lies. Thus, Maximin metric can be used as an additional fitness criterion to evaluate a population of solutions in order to provide observation of Pareto dominance along with the degree of clustering aiming to obtain diversity of solutions in Pareto sets. Consequently, the Pareto set can be updated at each iteration with alternatives from the population of solutions with the most negative (non-dominated) Maximin metric values. Equation (9) presents the Maximin metric of a normalized vector of objective functions Y i .
GA solution process starts with a population and applies genetic operators to search new feasible alternatives of solution in each iteration. The new population is inserted in the last Pareto set, the resulting set is evaluated by Maximin metric as ranking criterion and only solutions with negative Maximin metric are chosen to update the Pareto set. The process is repeated until the adopted maximum number of iterations is reached. Figure 2 exemplifies the adopted codification (chromosome) of alternative solutions (individuals) for four primary sources that can be inserted in two buses. In this example, the proposed solution suggests the installation of 4 DG units of type 1 at bus 13; for type 2, the installation of 5 DG units at bus 36 and 3 DG units at bus 58; for type 3, the installation of one DG unit at bus 3 and 2 DG units at bus 47; and, for type 4, the installation of 4 DG units at bus 62. The genetic operators applied in GA solution process are described as follow: • Elitism: best individuals from the current Pareto set are used as individuals in the next iteration. In this work, the number of best individuals to be considered in the new population is equal to 5% of the current population size. • Parental Selection: a tournament selection among three random individuals from the current iteration is applied. The individual with more negative Maximin metric wins the competition and becomes a father for the next generation. • Mutation: a fixed mutation tax of 5% is adopted. • Recombination: the process is always accomplished (recombination tax of 100%).
next iteration. In this work, the number of best individuals to be considered in the new population is equal to 5% of the current population size. • Parental Selection: a tournament selection among three random individuals from the current iteration is applied. The individual with more negative Maximin metric wins the competition and becomes a father for the next generation. • Mutation: a fixed mutation tax of 5% is adopted.

•
Recombination: the process is always accomplished (recombination tax of 100%). Three decision criteria considering the concept of preference were adopted to identify a final solution among Pareto frontier setpoints: minimum losses of active power point, maximum net present value (NPV) point, and Max-Min approximation point.
Minimum losses of active power criterion in this context represents the best solution for the DisCo viewpoint, since the utility usually seeks to improve the system efficiency.
Maximum net present value criterion represents the best solution for the DG owner viewpoint, since this one intends to sell the maximum energy with the lowest investment cost. Net present value for 10 years is calculated by Equation (10), where operation and maintenance costs are fixed as 14.5, USD/kWe.year for photovoltaic, 62.5 USD/kWe.year for onshore wind generation, and 44.76 USD/kWe.year for hydro generation [26].
In this case, the minimum attractiveness rate (MAR) is adopted by the DG decisionmaker.
Max-Min approximation (MMA) is a metric that calculates the gravity center of a set of points aiming the best compromise among all of the objective functions, being a benchmark for equilibrium among different criteria [27,29], since this metric is a nonpreference criterion for not considering the values of the decision-maker to choose a final solution. Given a Pareto set with l non-dominated solutions, the MMA metric for the i-th alternative within the final Pareto optimal set is calculated according to Equation (11), where index j is related to the j-th objective function.
Diversity of alternatives in Pareto sets can be evaluated according to the scattering of non-dominated solutions on the Pareto frontier. Spacing metric numerically evaluates how well the solutions in the Pareto frontier are distributed, so that value zero for this metric indicates that all the solutions are equidistantly spaced [33]. Given a Pareto set with l non-dominated solutions, the Spacing metric (S) of the set is determined by Equation (12). Three decision criteria considering the concept of preference were adopted to identify a final solution among Pareto frontier setpoints: minimum losses of active power point, maximum net present value (NPV) point, and Max-Min approximation point.
Minimum losses of active power criterion in this context represents the best solution for the DisCo viewpoint, since the utility usually seeks to improve the system efficiency.
Maximum net present value criterion represents the best solution for the DG owner viewpoint, since this one intends to sell the maximum energy with the lowest investment cost. Net present value for 10 years is calculated by Equation (10), where operation and maintenance costs are fixed as 14.5, USD/kWe.year for photovoltaic, 62.5 USD/kWe.year for onshore wind generation, and 44.76 USD/kWe.year for hydro generation [26].
In this case, the minimum attractiveness rate (MAR) is adopted by the DG decisionmaker.
Max-Min approximation (MMA) is a metric that calculates the gravity center of a set of points aiming the best compromise among all of the objective functions, being a benchmark for equilibrium among different criteria [27,29], since this metric is a nonpreference criterion for not considering the values of the decision-maker to choose a final solution. Given a Pareto set with l non-dominated solutions, the MMA metric for the i-th alternative within the final Pareto optimal set is calculated according to Equation (11), where index j is related to the j-th objective function.
Diversity of alternatives in Pareto sets can be evaluated according to the scattering of non-dominated solutions on the Pareto frontier. Spacing metric numerically evaluates how well the solutions in the Pareto frontier are distributed, so that value zero for this metric indicates that all the solutions are equidistantly spaced [33]. Given a Pareto set with l non-dominated solutions, the Spacing metric (S) of the set is determined by Equation (12).
where d i is the spacing factor of the i-th solution and d is the average value of all d i .

Results and Analysis
The modeling for the DG sizing and siting problem and the multi-objective solution methodology (MOHGA) proposed in this work are applied in two case-studies and main results are analyzed in graphs and tables. Case-studies are carried out in medium voltage systems: 69-bus test system known from literature (Case Study I) and a 918-bus Brazilian distribution system (Case Study II). Simulations were performed on a 2.21 GHz Intel Core i7-8750H CPU using MATLAB 2018b (9.5.0.944444). The elapsed time for Case Study I was 500 s and 4490 s for Case Study II, both with Popsize = 100 and 500 iterations.

Case Study I: 69-Bus Test System
First, the 69-bus and 12.66 kV distribution network [34] is simulated for a generalized analysis of the model and solution method. Table 1 presents the number of generation modules and power capacity of each DG plant, as well as the buses at which DG projects can be installed. Aiming to enlarge the search and solution spaces, a higher than usual 30~40% maximum penetration index PI max = 70% is considered in this case-study. The Pareto frontier obtained by the proposed multi-objective solution method based on Genetic Algorithm and Maximin metric (MOHGA) is presented in Figure 3, where y1, y2 and y3 are the normalized total installation cost of DG plants (C instal ), daily active power losses (LAE daily ), and daily amount of equivalent CO 2 emission (EqCO 2 E), respectively. Good scattering of Pareto frontier points (3.4%) certify the Maximin metric efficiency to aid multi-objective solution methods in the searching of diversity of solution alternatives. The scale of values for the axes in Figure 3 is related to the proposed normalization factors (value of objective functions referring to the solution obtained by miminization of daily active power losses considered as a mono-objective problem). Minimum power loss solution implies that normalized power losses of any alternative of solution achieve a range higher than 1.0, as observed in y2 axis. As minimization of daily active power losses seeks to maximum generation, the normalized values for the other two objective functions correspond to the maximum installation cost of DG plants and maximum daily amount of equivalent CO 2 emission, so that y1 and y3 axes attain values lower than 1.0.
Analysis of results identifies that the obtained Pareto-set presents 81% of individuals with photovoltaic in bus 13 with an average power of 272 kWe, 84% of individuals with wind source in bus 68 with an average power of 511 kWe, and 83% of its individuals with hydro source in bus 62 with an average power of 1.07 MWe. The Pareto-set also presents the set of buses 13 (PV)-68 (wind)-62 (hydro) in 60% of the individuals. Thus, analysis of Pareto-sets can also identify better locations for DG installations in distribution networks. Figure 4 presents the normalized objective functions plotted two by two aiming to verify conflict among them. It can be observed that the curvature in y1-y2 and y2-y3 planes indicates that, as mentioned, minimization of installation costs and minimization of CO 2 emissions are conflicting with minimization of network power losses. Figure 4 presents the normalized objective functions plotted two by two aiming to verify conflict among them. It can be observed that the curvature in y1-y2 and y2-y3 planes indicates that, as mentioned, minimization of installation costs and minimization of CO2 emissions are conflicting with minimization of network power losses.   Table 2 shows the phenotypes of solution alternatives chosen from Pareto-set, based on the three adopted decision criteria: minimum losses, maximum NPV and MMA point.
The solution based on the minimum losses of active power criterion attains comparatively larger quantities of DG generation modules, which reflects the DisCo viewpoint seeking to improve system efficiency by power loss reductions. Despite of higher installation costs of hydro and wind plants, the comparatively low capacity factor of photovoltaic plants may explain the absence of generation units from this primary energy source observed in the solution based on maximum net present value criterion, since the DG owner seeks to sell the highest amount of energy aiming a shorter payback time. Concerning the solution of the Max-Min approximation criterion, the comparatively  Table 2 shows the phenotypes of solution alternatives chosen from Pareto-set, based on the three adopted decision criteria: minimum losses, maximum NPV and MMA point. The solution based on the minimum losses of active power criterion attains comparatively larger quantities of DG generation modules, which reflects the DisCo viewpoint seeking to improve system efficiency by power loss reductions. Despite of higher installation costs of hydro and wind plants, the comparatively low capacity factor of photovoltaic plants may explain the absence of generation units from this primary energy source ob-served in the solution based on maximum net present value criterion, since the DG owner seeks to sell the highest amount of energy aiming a shorter payback time. Concerning the solution of the Max-Min approximation criterion, the comparatively lower number of generation modules reflects the needs for negotiations among conflicting goals to support a more efficient decision-making process.
The solution based on the criterion of minimum losses presents daily active power losses LAE daily = 1.93 MWh, which results a power loss reduction of 36.7% in comparison to the daily active power losses LAE daily = 3.05 MWh when no dispersed generation is allocated in the distribution network, although a installation cost of 4.74 × 10 6 USD will be necessary.
Regarding the maximum NPV criterion, a minimum attractiveness rate MAR = 15% and energy price 0.15 $/kWh are adopted in this work, so that the chosen solution presents a net benefit of 1.844 × 10 6 USD. If the DisCo will be the DG owner, LAE dayly = 1.95 MWh can be converted in USD using the same energy price and MAR as interest rate, which will result in savings of 7.4 × 10 6 USD to the DisCo. In this case, both minimum losses and maximum NPV can be adopted as preference criteria to the DisCo.
MMA decision criterion presents daily active power losses LAE daily = 2.3 MWh, which results in a power loss reduction of 24.6% (DisCo savings: 4.95 × 10 6 USD) compared to the case with no dispersed generation in the distribution network. MMA solution represents approximately 60% of installation cost and 50% of equivalent CO 2 emission compared to the minimum losses and maximum NPV criteria. Therefore, MMA can be adopted as a non-preference decision criterion aiming at common good among perspectives of concern.

Case Study II: 918-Bus Test System
The 918-bus/13.8 kV distribution network corresponds to a Brazilian urban feeder from Goiania, State of Goiás. The data for this system can be requested to the authors. Network buses on which DG projects can be inserted, as well as the number of generation modules and power capacity of each module are presented in Table 3. Only photovoltaic types were allowed since it is a downtown urban feeder. The maximum penetration index is fixed in PI max = 120%, so that power flow from the substation to the grid is possible. In PV-02 was given a higher efficiency and higher cost in comparison to previous cases, simulating a technological advance: c f = 20% and CO instal = 1500 USD/kWe, respectively.  Figure 5 shows the feeder picture (on Google Earth) and current during a typical day.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 11 of 16 Figure 5. Feeder picture (on Google Earth) and feeder current during a day for case II. Figure 6 depicts the Pareto frontier obtained by the proposed approach. In this case, y1, y2, and y3 are the normalized total installation cost, daily power losses and amount of equivalent CO2 emission, respectively. A well scattering of Pareto frontier points (1.3%) proves over again the efficiency of the Maximin metric to tackle multi-objective solution  Figure 6 depicts the Pareto frontier obtained by the proposed approach. In this case, y1, y2, and y3 are the normalized total installation cost, daily power losses and amount of equivalent CO 2 emission, respectively. A well scattering of Pareto frontier points (1.3%) proves over again the efficiency of the Maximin metric to tackle multi-objective solution methods in searching diversity of solution alternatives. The obtained Pareto-set presents two prevailing buses for installation of photovoltaic plants: for PV-01 type, 47% of the individuals propose bus 565 with 440 kWe of average power and 42% suggest bus 666 with 1440 kWe of average power; for PV-02 type, predominance is in bus 468 (95% of individuals) with an average power of 1265 kWe. Complementarily, the pair of buses 565-468 is in 47% of individuals and the pair 666-468 in 38%, totalizing 85% of solution alternatives from Pareto-set.
Normalized objective functions plotted two by two in Figure 7 and observation of curvature in y1-y2 and y2-y3 planes indicate over again that minimization of installation costs and of CO 2 emissions are conflicting with minimization of network power losses. Phenotypes of solution alternatives based on decision criteria are shown in Table 4, where MAR = 15% and energy price 0.15 $/kWh are adopted in maximum NPV criterion.
The comparatively larger quantities of generation modules observed in the solution based on minimum losses and maximum net present value criteria reflects an agreement Phenotypes of solution alternatives based on decision criteria are shown in Table 4, where MAR = 15% and energy price 0.15 $/kWh are adopted in maximum NPV criterion. The comparatively larger quantities of generation modules observed in the solution based on minimum losses and maximum net present value criteria reflects an agreement of goals to both DisCo and DG owner, since higher production of energy implies in lower losses to the DisCo seeking of efficiency while it results in higher sales and shorter payback time to the DG owner. The comparatively lower number of generation modules regarding the solution for the Max-Min approximation criterion reflects again the need for tradeoffs among conflicting goals in such a way as to ensure success in the decision-making process.
The solution based on both minimum losses and maximum NPV criteria presents daily active power losses LAE daily = 3.42 MWh at an installation cost of 5.57 × 10 6 USD and net benefit of 0.99 × 10 6 USD, resulting in a power loss reduction of 17.8% in comparison to the case with no dispersed generation (LAE daily = 4.16 MWh). When DG is owned by the DisCo, LAE dayly can also be converted in USD based on the same energy price and MAR as interest rate, which would result in savings of 4.88 × 10 6 USD to the DisCo.
MMA decision criterion results in daily active power losses (LAE daily ) of 3.72 MWh, corresponding with 10.6% of power loss reduction. In this case, the DisCo savings are 2.9 × 10 6 USD compared to the case with no dispersed generation. The solution obtained through MMA represents approximately 50% of equivalent CO 2 emission and 59% of installation cost compared to the minimum losses and maximum NPV criteria. It is important to mention that MMA criterion suggests installation of PV-02 type, since it is technologically superior, producing better results, even at a higher price. Thus, it can be concluded over again that MMA can be adopted as a non-preference decision criterion to satisfy varied perspectives of concern.

Conclusions
This paper presents a multi-objective modeling for the DG sizing and siting problem based on technical, economical and environmental aspects, where a multi-objective Genetic Algorithm and Maximin metric (MOHGA) were applied in the searching process to obtain Pareto sets with diversity of solutions. Locations (buses) for possible DG insertions and maximum capacity divided in generation modules were considered. Two case-studies were carried out to testify algorithm effectiveness and decision criteria considering the concept of preference were applied to choose a final solution from Pareto sets.
Management decisions typically involve considerations from a wide range of criteria and multi-objective modeling can be of interest to stakeholders, especially when tradeoffs among conflicting aspects are needed to support a more efficient decision-making.
Pareto frontier plots demonstrate that the applied methodology of solution process (MOHGA) obtained a satisfactory set of solution alternatives in terms of diversity, which can be understood as an acceptable sampling of the space of decision variables is achieved, thus providing an efficient option for solving multiobjective optimization problems.
The case-studies carried out demonstrate that several criteria (such as minimal losses solution, maximum NPV solution and MMA solution) allow decision makers (planners, DG owners, regulatory agencies, etc.) to develop strategies to meet their necessities (such as reduction of costs, CO 2 emissions and technical losses, etc.) at the same time, which is consistent with a more systematic strategy of generation expansion planning.
Intermittent sources with higher capacity factors also show a prevalence compared to the others, which can be noticed by the higher dimensions of wind and hydro DG's in case study I. This result points out the continuous need for a technological advance in DG technologies, especially photovoltaic, to increase their capacity factors before other energy sources. Finally, this work also illustrates the high complexity of renewables integration in distribution networks and emphasizes the necessity of research in this field.

A k−m
Cable ampacity of k-m branch c f n Capacity factor of n-th DG plant EqCO 2 EF n Equivalent CO 2 emission factor of n-th DG plant CO n instal