Multi-Objective Optimization Approach for Placement of Multiple DGs for Voltage Sensitive Loads

This paper presents the optimal placement of multiple Dispersed Generators using multi-objective optimization. The optimization is carried out with objectives namely active power loss, reactive power loss, voltage deviation and overall economy. The multi-objective optimization and accounting conflicting objectives are realized through Particle Swarm Optimization with fuzzy decision approach to find the optimal sizes and sites of Dispersed Generators for voltage dependent residential, commercial and industrial loads. The clusters of buses are formulated from base case load flow to limit the search space for finding the placement of the Dispersed Generators. The effectiveness of the proposed approach is tested on a 69-bus radial distribution. It is found that the optimal placement of the Dispersed Generators improves the overall performance of the system and the optimal allocation is affected by the type of load.


Introduction
The distribution system is the final interconnection between the bulk power system and consumers.There is a large number of buses and lines with a high R/X ratio, which results into high power losses in distribution systems.These losses are even more significant at heavy loads.Moreover, the voltage of buses, farther from the substation, decreases.It is important to optimize the operation of distribution systems by reducing the losses and enhancing the voltages profile.To assure these improvements reconfiguration and capacitor placement have been employed earlier but in recent decade the placement of dispersed generation has attracted the researchers.Dispersed Generators (DGs) are small-scale generators that are located near to load centers.The size of DGs varies from less than a kW to tens of MW.The DGs have additional advantage of meeting load demand along with reduction in losses and improvement of voltage profile.To maximize these benefits, DG of optimal size should be placed at optimal site.Several optimization techniques have been utilized for DG placement by optimizing various objective functions, as single objective and multi-objective optimization problem.
The optimal placement of DG has been investigated through PSO [1,2] to minimize losses in distribution system.The effectiveness of three different types of DGs have been analyzed [1].The DG that is capable of supplying both active and reactive power is most effective for loss reduction.Loss minimization has been investigated by integrating multiple DGs [2].The realistic nature of loads has been accounted by time-varying characteristics.GA [3] has been applied to optimally allocate DGs to minimize losses.It has been concluded that the benefit of loss minimization by DG placement is limited up to certain locations only and after that trend reverses.The problem of DG placement to minimize losses has been solved by ABC algorithm [4] for two different load scenarios.Multiple DG placement has been obtained by Bat Algorithm [5] to achieve maximum reduction in annual energy losses in distribution system.
An analytical heuristic based on exact loss formula has been derived to attain the optimal size and location of DG [6] to reduce the power losses in distribution system.The location of DG, capable of supplying only active power, is found by sensitivity loss factors.The heuristic proposed in [6] is elaborated to consider different types of DG [7].Another classical technique utilizes loss formula based on bus-injection to branch-current (BIBC) and branch-current to bus-voltage (BCBV) matrices [8] to calculate the optimal size and site of DG.This method does not require calculation of impedance matrix.Another technique based on analytical expressions [9] obtain the optimal site and size of four different types of DG by minimizing the losses.A decision making approach [10] based on exhaustive load flow has been employed to allocate wind and photovoltaic (PV) based DGs in distribution system to reduce the power losses and improve voltage profile.The PV based DG is found more influential in voltage profile improvement in comparison to wind based DG.
The multi-objective optimization problem has been solved by two different approaches i.e., by representing multi-objective optimization as single objective optimization using weighed sum approach or by handling multiple objectives simultaneously.The common objectives for multi-objective DG placement are minimization of power loss, voltage deviation, customer interruptions, cost of purchased energy and cost of DG.
The weighted sum approach has been used in [11][12][13][14][15] to find the optimal size and optimal site of DG.The allocation has been investigated through GA [11] where the multi-objective problem is changed to single objective by assigning weights to four objectives namely real and reactive power loss index, voltage profile index and Mega Volt Ampere (MVA) capacity index.In addition to these four objectives short-circuit current is also considered as objective in [12].The optimization has been carried out using PSO after changing the multi-objective problem to single objective problem.Optimal site of DG for predefined size has been obtained by Simulated Annealing [13] for minimizing power loss, emissions and contingency.Optimal allocation of multiple DGs based on wind turbine and PV array has been done by PSO [14] for minimizing power loss and improving voltage stability.The multi-objective optimization problem is converted to single objective by assigning weights, where the weights were updated iteratively.Weighted sum approach has also been utilized in [15] for optimal sizing and siting of voltage controlled DG for balanced as well as unbalanced distribution networks for objectives as real and reactive power loss index, voltage profile index and reserve capacity of conductor index.The optimization has been attempted by supervised big-bang crunch method.The limitation of weighted sum approach is that the solution is dependent on weight assigned to each objective or the decision maker who assigns the weights to the objectives, influences the solution.Such solutions may be non-optimal solutions.
Multi-objective optimization problem has been solved by ε-constrained technique in [16][17][18] to suggest optimal allocation of DG in which one objective is selected as main objective and rest of objectives are treated as constraints.Minimization of cost of energy losses, voltage profile and power quality index are selected as objectives in [16] to attempt optimization by ε-constrained technique.The ε-constrained technique has been combined with GA [17] for optimal allocation of DG by minimizing four costs network upgrading, purchased energy, energy losses, energy not supplied.The combination of ε-constrained technique and GA [18] has been employed to allocate DG optimally considering voltage dependent loads.The minimization of interruptions to customers, cost of active power losses, cost of energy purchased, bus voltage deviation and maximization of security margin has been considered as multiple objectives.
Optimal location and size of DG has been obtained by Strength Pareto Evolutionary Algorithm and Non-dominated Sorting GA-II in [19] and Shuffled frog leaping algorithm in [20] to minimize real power loss and cost of DG.Optimization of weight factor is additional parameter to be optimized in [20] which balances the two main objectives.In addition to minimization of real power losses of system and cost of DG, voltage deviation minimization is also considered as objective in [21].Non-dominated Sorting GA-II was used to solve the problem of optimal allocation of specified number of DGs of multiple types.The pareto front is obtained for simultaneous optimization of three objectives and final solution is selected from non-dominated solutions by normalized membership functions.Optimal allocation of DG models has been done by minimizing voltage deviation index, real and reactive power loss index, line loading index and short-circuit index using fuzzy embedded GA [22].The multi-objective problem is simplified by assigning fuzzy membership function to each index.Goal optimization algorithm [23] has been used for optimal allocation of DG in distribution system by achieving the set goal for the indices used in [22].In [24] combination of GA and PSO has been utilized by where optimal site is obtained by GA and size of DG is obtained by PSO.Multiple DGs of optimal size has been placed to minimize losses, improve voltage profile and to optimize voltage stability index.Enhanced multi-objective PSO has been proposed to optimally allocate multiple DGs by minimizing power losses and maximize voltage stability index [25,26].Fuzzy decision-making has been employed to select single solution from Pareto front.Sequential quadratic programming has been employed in [27] to obtain set of optimal solutions to minimize total real power loss and cost of DG.In addition to loss and voltage deviation minimization, cost minimization is also considered as objective [28].Shuffled Bat algorithm has been employed in solving optimization.Pareto front has been obtained in both and compromised solution has been selected from set of solutions.The best solution has been selected by fuzzy decision-making in [27] and while it has been done by max-min technique in [28].
The GA and PSO are commonly used techniques for optimal placement of DG.The PSO has better diversity and exploration in comparison to GA and its convergence is fast due to momentum effect on particles [29].Therefore, the optimization in the presented work is carried out through PSO.
In this paper the four indices are selected for allocation of DG of optimal size at optimum location in radial distribution system (RDS).These indices are active power loss index (PLI), reactive power loss index (QLI), voltage deviation index (VDI) and overall economy index (OEI).The multi-objective optimization is implemented using fuzzy decision approach by assigning fuzzy memberships to these indices and optimizing the resultant membership.The optimal sizes and sites of DGs are obtained through PSO for voltage dependent residential, commercial and industrial loads.The optimization is carried out for placement of single as well as multiple DGs.The clustering of buses based on base case load flow is done to reduce the search space.

Objectives for DG Placement
The placement of DG in RDSs affects active power losses, reactive power losses and, voltage profile of buses.If DG of appropriate size is placed at optimum site then the losses will decrease and voltage profile of the buses gets improved.In the presented work four individual objectives are accounted simultaneously for the multi-objective formulation.Corresponding to these objectives, the values of indices before placement of DG is taken as unity and it is desired to reduce the value of these indices by placement of optimal sized DG at optimal site.These indices are formulated as:

Active Power Loss Index (PLI)
With the integration of DG in RDSs, active power losses are expected to reduce.The active power loss index (PLI) gives the effectiveness of DG integration on active power loss and is defined as: This objective attempts to minimize active power loss in distribution system by addition of DG.

Reactive Power Loss Index (QLI)
The reactive power loss also varies with integration of DG in RDSs.The effectiveness of DG integration on reactive power losses is given by QLI, which is defined as: This objective minimizes reactive power loss in distribution system by addition of DG.

Voltage Deviation Index (VDI)
The voltage deviation (VD) is a measure of the deviation of bus voltages from reference bus.The effectiveness of improvement in voltage profile of buses after optimal placement of DG is measured by voltage deviation index (VDI).It is defined as: where VD DG and VD are measure of deviation in the bus voltage with respect to reference voltage after and before integration of DG respectively.These are defined as: Economy is accompanying every activity in today's competitive market.The overall economy index (OEI) is minimized so as to make DG operation economical.The OEI is expressed as: CE DG and CE are stated as: Energies 2017, 10, 1733 5 of 17

Multi-Objective Planning for DG Placement
A general multi-objective problem can be formulated mathematically as: where X is solution vector, S is solution space, f i (X) is ith objective function, h(X) is equality constraint and C(X) is inequality constraint.
The objectives in the multi-objective optimization are usually conflicting.The objectives that are optimized simultaneously are expressed through indices namely PLI, QLI, VDI and OEI.Instead of using the true Pareto based optimization, the objectives are combined using the membership functions derived from the fuzzy decision approach.The fuzzy membership function values (µ) are computed for each objective and are defined as: , for y 1 < y j < y 2 ; (7) where y j is the value of performance index.
In this presented work the values of y 2 and y 1 are considered as 1.25 and 0 respectively.The multi-objective function (MOF) formulated by combining four objectives functions can be expressed as: where µ y is the fuzzy membership assigned to index y.The objective function represented by Equation ( 8) is subjected to inequality constraints of bus voltage limits and MVA capacity of lines.

Bus Voltage Constraints
The incorporation of DGs in RDS boosts the voltage level of all buses, which may result in over-voltage at some buses.Moreover, there are always limits on lower level of voltage so the bus voltage constraints are defined as: The limits for V min i and V max i are taken as 0.90 and 1.05 p.u. [30].

MVA Capacity Constraints
The integration of DG in RDS changes the nature of RDS from passive to active and hence the current levels change.This constraint is necessary to ensure thermal limits of lines.MVA capacity constraint is expressed as:

Limiting Search Space
The search space for optimal allocation of DG is limited by making the clusters of buses.Bus voltages obtained after the load flow for the base case are used to decide these clusters.These clusters represent the group of buses having significant low voltage.The potential buses for the placement of DGs are selected from the cluster of buses, and thus limiting the search space.Such identified cluster of buses although reduces the search space but always yields optimal solution and thereby does not compromise on the quality of the optimal solution.

Voltage Dependent Loads
The effect of node voltage on connected load is accounted while evaluating the optimal site and size of DG for single as well as multi-objective optimization.The practical loads such as industrial, commercial, residential are prevalent in the RDS.These practical loads affect the bus voltage or bus voltages affect the power drawn by these loads.These loads are modeled as: Here np and nq are the active and reactive power exponents.These exponents are considered as zero for constant load, which is conventionally, used in power flow studies.The value of these exponents for various types of loads is tabulated in Table 1 [18].

Algorithm for DGs Allocation
PSO is a population-based search and optimization technique in which individuals called particles change their position (state) with time.In PSO, particles move through the search space to keep itself in the best position pbest.The best position searched by the group particles gbest is assessable to each particle.The particles change their position with the help of their current positions, their pbest and gbest.The Figure 1 shows the flowchart for implementing the multi-objective DG allocation problem using PSO, where algorithm steps are described as follows: Step 1 Read the distribution system line and bus data and the type of load.
Step 2 Carryout the base case load flow using Bus injection to branch current (BIBC) and branch current to bus voltage (BCBV) matrices [31].Step 3 Calculate the indices i.e., PLI, QLI, VDI and OEI before placement of DG.
Step 4 Identify the cluster of buses for obtaining optimal site of DG.
Step 6 Initialize the particles positions that are random DGs at random buses in the search space and set the iteration counter k = 0. Step 6 For these positions, calculate load flow and check for constraints satisfaction.Compute the performance indices and corresponding membership functions using Equation (7).
Step 7 Obtain the objective function from Equation (8).
Step 8 For the positions, assign them as individual best positions pbest.Also identify the best of the pbest as global best gbest.
Step 9 Update the particles positions by accounting their adaptive weight and current positions [1], pbest and gbest and increase iterations as k = k + 1.
Step 10 Go to Step 6 until the convergence is achieved.

Results and Discussions
The methodology adopted for optimal allocation of DG in this presented work is tested on IEEE 69-bus RDS [32].This system has total load of 3.800 MW and 2.690 MVAr and it is resulting into 225.005kW active power losses and 102.173 kVAr reactive power losses.The algorithm is implemented under MATLAB (R2014b, MathWorks, Natick, MA, USA) in MAC environment.The presented algorithm has been executed while considering the swarm size as 50 for 50 iterations.For multi-objective optimization the optimal size and site of DG are found for the placement of single as well as multiple DGs.The DGs that are capable of supplying only active power, is considered for optimization.
The search space for obtaining optimal site of DG is obtained by forming the clusters of buses with significantly lower bus voltage from base-case load flow for all types of loads.The variation of

Results
The methodology adopted for optimal allocation of DG in this presented work is tested on IEEE 69-bus RDS [32].This system has total load of 3.800 MW and 2.690 MVAr and it is resulting into 225.005kW active power losses and 102.173 kVAr reactive power losses.The algorithm is implemented under MATLAB (R2014b, MathWorks, Natick, MA, USA) in MAC environment.The presented algorithm has been executed while considering the swarm size as 50 for 50 iterations.For multi-objective optimization the optimal size and site of DG are found for the placement of single as well as multiple DGs.The DGs that are capable of supplying only active power, is considered for optimization.
The search space for obtaining optimal site of DG is obtained by forming the clusters of buses with significantly lower bus voltage from base-case load flow for all types of loads.The variation of bus voltage for base case is shown in Figure 2. Two clusters of buses with minimum voltage are selected to obtain the optimal site of DG.First cluster of buses consists of buses from 57 to 65 and second cluster of buses ranges from bus 16 to 28.
bus voltage for base case is shown in Figure 2. Two clusters of buses with minimum voltage are selected to obtain the optimal site of DG.First cluster of buses consists of buses from 57 to 65 and second cluster of buses ranges from bus 16 to 28.

DG Placement Using Single Objective Optimization
Optimal size and site of DG are found when all four indices are considered independently for placement of single as well as multiple DGs.For placement of single DG, only first cluster of buses is considered while for placement of two DGs, one DG is selected from each of the respective cluster.The four cases that are investigated for DG placement through single objective optimization are: Minimizing PLI, Minimizing QLI, Minimizing VDI and Optimizing OEI.

Minimizing PLI
With PLI as objective, the optimal size and site of DG and corresponding indices after optimal allocation of DG for all types of loads are presented in Table 2.The value of PLI is 1.0 before placement of DG in RDS.Its value reduces significantly after integration of optimal sized DG at optimal site due to reduction in active power loss.For constant load the improvement is maximum, where its value reduces to 0.3696 and 0.3280 corresponding to placement of single and multiple DGs respectively.The values of all other indices also reduce after placement of DG.

DG Placement Using Single Objective Optimization
Optimal size and site of DG are found when all four indices are considered independently for placement of single as well as multiple DGs.For placement of single DG, only first cluster of buses is considered while for placement of two DGs, one DG is selected from each of the respective cluster.The four cases that are investigated for DG placement through single objective optimization are: Minimizing PLI, Minimizing QLI, Minimizing VDI and Optimizing OEI.

Minimizing PLI
With PLI as objective, the optimal size and site of DG and corresponding indices after optimal allocation of DG for all types of loads are presented in Table 2.The value of PLI is 1.0 before placement of DG in RDS.Its value reduces significantly after integration of optimal sized DG at optimal site due to reduction in active power loss.For constant load the improvement is maximum, where its value reduces to 0.3696 and 0.3280 corresponding to placement of single and multiple DGs respectively.The values of all other indices also reduce after placement of DG.  3 for placement of single as well as multiple DGs.In addition to reduction in QLI all other indices also reduces i.e., performance is improved.The reduction in QLI is almost same for both objectives PLI and QLI as both are dependent on current.The deviation in bus voltages in RDS is an important factor, which can be improved by integration of DG.The optimal size, site of DG and performance indices are summarized in Table 4 corresponding to minimization of VDI for placement of single as well as multiple DGs.The size of DG for VDI optimization is higher in comparison to the sizes of DGs obtained for optimizing PLI and QLI independently.The improvement in VDI results in improvement of voltage of all buses in RDS.The placement of multiple DGs is found more effective as the reduction in VDI is more as compared to single DG placement.The optimization of overall economy index OEI is very important objective as economy is major decision factor.For OEI formulation, the cost of power supplied by DG is taken as 46 US$/MWh and power tariff is taken as 44.5 US$/MWh [18].Corresponding to the optimization of OEI, the optimal placement and values of the indices are summarized in Table 5.The reduced value of OEI results in economic operation of DG.It is observed that placement of multiple DGs is more economical as compared to single DG.The obtained OEI for placement of two DG is always lower than its value for single DG.The reduction in OEI follows the pattern of PLI, as OEI is dependent on active power losses given by Equation (5).The size of DG is higher for constant load model among four types of loads for every objective and hence the reduction in all objectives is also maximum for constant load model.

DG Placement Using Multi-Objective Optimization
The placement of DG using multi-objective optimization is discussed for two cases.In first case PLI and VDI are considered simultaneously and in second case all the four objectives (PLI, QLI, VDI and OEI) are considered simultaneously.

Validation of Proposed Technique
The results obtained from proposed PSO based technique are compared with multi-objective GA for validation.The comparison of DG size, site and performance indices obtained by proposed technique with GA for single objective as well as multi-objective case is presented in Table 6.In addition to multi-objective case the comparison of two techniques is also done for single objective case.The constant load model is considered for comparison and is equally applicable to other load models.The optimal sizes of DG for the placement of single as well as multiple DGs obtained by both techniques are comparable.In implementation of GA true Pareto front is obtained while in PSO fuzzy memberships are assigned that yields in fast convergence.The true Pareto front for the placement of single DG as well as multiple DGs is shown in Figures 3 and 4

Case I Optimizing both PLI and VDI Simultaneously
The placement of single as well as multiple DGs has been investigated by considering PLI and VDI simultaneously by assigning fuzzy membership function to each performance index.For all types of voltage dependent loads optimal size, site and corresponding value of each index after placement of DGs are summarized in Table 7.These results are compared to the single objective optimization results presented in Tables 2  and 4, which summarize the results corresponding to minimization of PLI and minimization of VDI respectively.From Tables 2 and 4, it is clear that the value of PLI is minimum, while VDI is

Case I Optimizing both PLI and VDI Simultaneously
The placement of single as well as multiple DGs has been investigated by considering PLI and VDI simultaneously by assigning fuzzy membership function to each performance index.For all types of voltage dependent loads optimal size, site and corresponding value of each index after placement of DGs are summarized in Table 7.These results are compared to the single objective optimization results presented in Tables 2  and 4, which summarize the results corresponding to minimization of PLI and minimization of VDI respectively.From Tables 2 and 4, it is clear that the value of PLI is minimum, while VDI is

Case I Optimizing both PLI and VDI Simultaneously
The placement of single as well as multiple DGs has been investigated by considering PLI and VDI simultaneously by assigning fuzzy membership function to each performance index.For all types of voltage dependent loads optimal size, site and corresponding value of each index after placement of DGs are summarized in Table 7.These results are compared to the single objective optimization results presented in Tables 2  and 4, which summarize the results corresponding to minimization of PLI and minimization of VDI respectively.From Tables 2 and 4, it is clear that the value of PLI is minimum, while VDI is maximum when only PLI is the objective, whereas VDI is minimum and PLI is maximum when only VDI is the objective for all types of loads.
The results of multi-objective optimization, while compared with single objective optimization, PLI and VDI forms the non-dominated set of solutions.The variation of PLI with variation in VDI for placement of single as well as multiple DGs is shown in Figure 5a-d for different types of load models.From the Pareto front, a set of optimal solutions obtained by multi-objective optimization, the best compromised solution is obtained.The compromised solutions are marked as encircled solutions while; the extreme solutions are for individual objectives.As the optimal solution set corresponding to the multiple DG placement is more dominant i.e. near to the origin, which suggests that the placement of multiple DGs is more beneficial and multiple DG placement reduces values of indices as compared to single DG placement.
maximum when only PLI is the objective, whereas VDI is minimum and PLI is maximum when only VDI is the objective for all types of loads.
The results of multi-objective optimization, while compared with single objective optimization, PLI and VDI forms the non-dominated set of solutions.The variation of PLI with variation in VDI for placement of single as well as multiple DGs is shown in Figure 5a-d for different types of load models.From the Pareto front, a set of optimal solutions obtained by multi-objective optimization, the best compromised solution is obtained.The compromised solutions are marked as encircled solutions while; the extreme solutions are for individual objectives.As the optimal solution set corresponding to the multiple DG placement is more dominant i.e. near to the origin, which suggests that the placement of multiple DGs is more beneficial and multiple DG placement reduces values of indices as compared to single DG placement.

Case-II Optimizing All Indices Simultaneously
In second case of multi-objective optimization based placement of DG, all four indices PLI, QLI, VDI and OEI are combined by assigning fuzzy membership to each index.The optimal size and site of DG and all indices are tabulated in Table 8.It is observed from results that type of load is affecting DG size and site significantly.

Case-II Optimizing All Indices Simultaneously
In second case of multi-objective optimization based placement of DG, all four indices PLI, QLI, VDI and OEI are combined by assigning fuzzy membership to each index.The optimal size and site of DG and all indices are tabulated in Table 8.It is observed from results that type of load is affecting DG size and site significantly.
The placement of multiple DGs is more advantageous as compared to single DG placement and it results in reduced values of indices.The reduced values of PLI and QLI show reduction in active and reactive power losses respectively and reduced values of VDI show improvement in voltage profile of all buses.The reduction in OEI means that integration of DG in distribution system is economical as it reduces the losses and cost of DG is being overcome by reduction in losses.The placement of DGs of optimal size at optimal site improves the voltage of all buses of the RDS.The variation of voltage of each bus after and before placement of optimal sized DG at optimal site for multi-objective optimization after considering all indices is shown in Figures 6-9 for four different types of loads.The plot of voltage indicates that the improvement in bus voltage is better with placement of multiple DGs as compared to single DG.Also the voltage of all buses is within the limits of voltage constraints for all types of load after placement of single as well as multiple DG.
Energies 2017, 10, 1733 13 of 18 The placement of multiple DGs is more advantageous as compared to single DG placement and it results in reduced values of indices.The reduced values of PLI and QLI show reduction in active and reactive power losses respectively and reduced values of VDI show improvement in voltage profile of all buses.The reduction in OEI means that integration of DG in distribution system is economical as it reduces the losses and cost of DG is being overcome by reduction in losses.

Bus Voltage after DG Placement
The placement of DGs of optimal size at optimal site improves the voltage of all buses of the RDS.The variation of voltage of each bus after and before placement of optimal sized DG at optimal site for multi-objective optimization after considering all indices is shown in Figures 6-9 for four different types of loads.The plot of voltage indicates that the improvement in bus voltage is better with placement of multiple DGs as compared to single DG.Also the voltage of all buses is within the limits of voltage constraints for all types of load after placement of single as well as multiple DG.The placement of multiple DGs is more advantageous as compared to single DG placement and it results in reduced values of indices.The reduced values of PLI and QLI show reduction in active and reactive power losses respectively and reduced values of VDI show improvement in voltage profile of all buses.The reduction in OEI means that integration of DG in distribution system is economical as it reduces the losses and cost of DG is being overcome by reduction in losses.

Bus Voltage after DG Placement
The placement of DGs of optimal size at optimal site improves the voltage of all buses of the RDS.The variation of voltage of each bus after and before placement of optimal sized DG at optimal site for multi-objective optimization after considering all indices is shown in Figures 6-9 for four different types of loads.The plot of voltage indicates that the improvement in bus voltage is better with placement of multiple DGs as compared to single DG.Also the voltage of all buses is within the limits of voltage constraints for all types of load after placement of single as well as multiple DG.The effect of selection of DG on voltage profile of all buses is shown in Figure 10 for constant load model.Significant improvement in voltage of all buses is achieved after optimal integration of DG corresponding to all individual objectives and two cases of multi-objective optimization but it is maximum for placement of multiple DGs corresponding to objective VDI.The effect of selection of DG on voltage profile of all buses is shown in Figure 10 for constant load model.Significant improvement in voltage of all buses is achieved after optimal integration of DG corresponding to all individual objectives and two cases of multi-objective optimization but it is maximum for placement of multiple DGs corresponding to objective VDI.The effect of selection of DG on voltage profile of all buses is shown in Figure 10 for constant load model.Significant improvement in voltage of all buses is achieved after optimal integration of DG corresponding to all individual objectives and two cases of multi-objective optimization but it is maximum for placement of multiple DGs corresponding to objective VDI.The effect of selection of DG on voltage profile of all buses is shown in Figure 10 for constant load model.Significant improvement in voltage of all buses is achieved after optimal integration of DG corresponding to all individual objectives and two cases of multi-objective optimization but it is maximum for placement of multiple DGs corresponding to objective VDI.

Comparative Results
The proposed technique is validated by comparing with well-known multi-objective GA as presented in Table 6.Further the obtained results are compared with the available results in literature for placement of single as well as multiple DGs.The comparative analysis is presented in Table 9. 1.81 @ 61 62.95 GA [3] 1.87 @ 61 63.01 ABC [4] 1.90 @ 61 62.97 Analytical [6] 1.81 @ 61 62.86 Analytical [8] 1.81 @ 61 59.09 Analytical [9] 1

Discussion
The optimal placement of single as well as multiple DGs has been investigated for voltage sensitive loads through the optimization of objectives represented through indices PLI, QLI, VDI and OEI.The single and multi-objective optimization has been investigated for the placement of single DG and placement of multiple DGs.As summarized in Section 6.1 through Tables 2-5, the values of PLI, QLI, VDI and OEI through single objective optimization are higher for single DG placement then the respective values obtained with the placement of two DGs.This suggests that the placement of two DGs is yielding better optimization.The comparative results have been included in Table 9, where the results of proposed method have been compared with published results obtained through various heuristics PSO [1], GA [3], ABC [4], and Analytical methods [6,8,9] for minimizing losses.The proposed method is giving better performance.The multi-objective optimization, as summarized in Section 6.2 through Tables 7 and 8, yields the lower values of indices with two DG placement in comparison to single DG placement.The Pareto-front as shown in Figure 5 is closer to origin for two DG.This suggests that two DG placement is yielding better optimization compared to single DG placement.With multiple DGs of optimum size, the compensation is affected at multiple locations and thereby the load burden is not getting reflected up to substation.Such multiple DG placement thereby resulting into reduced losses, improved bus voltage profile.The ideal numbers of DG placement is an area of research and beyond the scope of presented manuscript.

Conclusions
In this paper a fuzzy decision based multi-objective optimization approach is presented to decide the optimal placement of single and multiple DGs.The optimization problem is solved by PSO to yield optimal solution.The search space for selection of optimal sites for DG placement is reduced by forming the cluster of buses.The proposed approach is implemented on 69-bus RDS for different voltage sensitive residential, commercial and industrial loads.It is observed that the load types significantly affects the sites and sizes of DGs.The significant improvement in the voltage profile is observed at all buses after integration of DG of optimal size at optimal site.Moreover, the placement of multiple DGs is found to be more effective in improving the system performance as compared to placement of single DG.Two DGs are sufficient to improve the overall performance for both single and multi-objective optimization for different loads.
Author Contributions: Navdeep Kaur has done this work and prepared the manuscript under supervision of Sanjay Kumar Jain.

Figure 1 .
Figure 1.Flowchart for the proposed DG allocation technique.

Figure 1 .
Figure 1.Flowchart for the proposed DG allocation technique.

Figure 2 .
Figure 2. Variation of base-case bus voltage for all loads. respectively.

Figure 3 .
Figure 3. Pareto front for the placement of single DG.

Figure 4 .
Figure 4. Pareto front for the placement of multiple DG.

Figure 3 .
Figure 3. Pareto front for the placement of single DG.

Figure 3 .
Figure 3. Pareto front for the placement of single DG.

Figure 4 .
Figure 4. Pareto front for the placement of multiple DG.

Figure 4 .
Figure 4. Pareto front for the placement of multiple DG.

Figure 5 .
Figure 5. Variation of PLI with variation in VDI for placement single and multiple DGs.(a) constant load; (b) commercial load; (c) residential load; (d) industrial load.

Figure 5 .
Figure 5. Variation of PLI with variation in VDI for placement single and multiple DGs.(a) constant load; (b) commercial load; (c) residential load; (d) industrial load.

Figure 6 .
Figure 6.Variation of bus voltage before and after placement of DG for Constant Load.

Figure 7 .Figure 6 .
Figure 7. Variation of bus voltage before and after placement of DG for Commercial Load.

Figure 6 .
Figure 6.Variation of bus voltage before and after placement of DG for Constant Load.

Figure 7 .Figure 7 .
Figure 7. Variation of bus voltage before and after placement of DG for Commercial Load.

Figure 8 .
Figure 8. Variation of bus voltage before and after placement of DG for Residential Load.

Figure 9 .
Figure 9. Variation of bus voltage before and after placement of DG for Industrial Load.

Figure 8 .
Figure 8. Variation of bus voltage before and after placement of DG for Residential Load.

Figure 8 .
Figure 8. Variation of bus voltage before and after placement of DG for Residential Load.

Figure 9 .
Figure 9. Variation of bus voltage before and after placement of DG for Industrial Load.

Figure 9 .
Figure 9. Variation of bus voltage before and after placement of DG for Industrial Load.

Figure 8 .
Figure 8. Variation of bus voltage before and after placement of DG for Residential Load.

Figure 9 .
Figure 9. Variation of bus voltage before and after placement of DG for Industrial Load.

Figure 10 .
Figure 10.Variation of bus voltage showing effect of objective for constant load.

Table 2 .
Optimal size and site of DG and Indices values for power loss index (PLI) as objective., the resulted optimal DG size, site and corresponding values of indices are presented in Table3for placement of single as well as multiple DGs.In addition to reduction in QLI all other indices also reduces i.e., performance is improved.The reduction in QLI is almost same for both objectives PLI and QLI as both are dependent on current.

Table 2 .
Optimal size and site of DG and Indices values for power loss index (PLI) as objective.With QLI minimization, the resulted optimal DG size, site and corresponding values of indices are presented in Table

Table 3 .
Optimal size and site of DG and Indices values for reactive power loss index (QLI) as objective.

Table 4 .
Optimal size and site of DG and Indices values for voltage deviation index (VDI) as objective.

Table 5 .
Optimal size and site of DG and Indices values for OEI as objective.

Table 6 .
Validation of proposed technique.

Table 7 .
Optimal size and site of DG for the placement of single DG and multiple DG and Indices values for multi-objective optimization (PLI and VDI).

Table 7 .
Optimal size and site of DG for the placement of single DG and multiple DG and Indices values for multi-objective optimization (PLI and VDI).

Table 7 .
Optimal size and site of DG for the placement of single DG and multiple DG and Indices values for multi-objective optimization (PLI and VDI).

Table 8 .
Optimal size and site of DG for single and multiple DG placement and Indices values for multi-objective optimization (all objectives).

Table 8 .
Optimal size and site of DG for single and multiple DG placement and Indices values for multi-objective optimization (all objectives).