Optimal distributed generation location and sizing for loss minimization and voltage profile optimization using ant colony algorithm

A system of power generation whereby the generating equipment is located close to the point of usage, thereby reducing losses and operation cost is called distributed generation (DG). However, it is imperative that DGs are sited such that the quality of power delivered is optimized and the total real power loss within the system minimized. This paper proposes an approach for optimum sizing and siting of DGs sizing in a power distribution system using Ant Colony Optimization (ACO) algorithm. To validate the algorithm the IEEE 30 bus standard test system was employed. A 92% decrease in real power loss within the system relative to the value before the connection of DGs was observed, while the minimum bus voltage increased from 0.656 per unit to 0.965 per unit. The results obtained from ACO are further verified by creating an ETAP model of the IEEE 30 bus system and simulating the impact of DG on the system. A significant reduction in total real power losses within the system and improvement in voltage profile was observed when the DGs are placed at the ACO derived sites relative to at other locations. Therefore, Ant Colony Algorithm can be used in deriving the optimum sites and sizes of DGs in a power distribution system.


Introduction
Distributed Generators (DGs) are power generation systems whose output are not connected to a central grid structure for transmission over long distances to the point of usage, rather they are located at the point of usage, with limited or no requirement for transportation [1,2]. Due to their modularity and portability and little or no requirement for the transportation of the energy resources (fuel) from point of production to the point of usage, renewable energy technologies such as solar photovoltaic and solar thermal systems which harnesses solar energy from the sun and transform it into a suitable form of energy such as heat and electricity are prime candidates for distributed energy generation [3]. Other forms of renewable energy technologies applied in distributed generation systems are small hydro power systems, wind turbines, biogas and biothermal systems and so on. Also, due to the harmful effect of fossil fuel usage on the environment such as the greenhouse effect and global warming, there has been a significant increase in the usage of renewable energy technologies.
Thus, due to the variable nature of solar and other renewable energy resources, as their penetration rate into the grid increases, they pose some technical challenges on the system such as feeder overloading, harmonic pollution and so on. Also, high capital cost, low reliability and low efficiency limits the widespread usage of these technologies [1]. Besides, voltage flickers and power fluctuations due to variations in solar irradiation are some of the adverse effects of high penetrated PV in power systems [4]. In their research, Miller and Ye [5] evaluated the impact of using DG as a significant portion of the total energy resources by creating a model of Jeju Island, Korea and made use of three distributed generation systems to assess the potential negative power quality effects of high penetration of DG. Also, Ajan and Nirmala [6] modelled a stand-alone and grid-tied PV system using Simulink which was used to evaluate the relationship between PV system parameters and the output. The various Voltages, Current and Power levels are recorded. The various advantages of optimally allocated DGs in a distribution network such as improvement in overall system voltage, reduction in total system power loss, improvement in system power quality and frequency regulation are highlighted in [7,8].
However, to derive maximum benefits from DG, the generators must be appropriately sized and located within the power distribution network [7,9,10]. To achieve this objective, various researchers have proposed different approaches. These methods are classified based on the driving algorithms into metaheuristic, analytical or numerical based methods. Analytical techniques for solving the DG allocation problem are presented in [11][12][13]. Several authors have also published research studies on optimal allocation of DG and reactive power sources such as capacitors in a distribution system using numerical based techniques [14]. One of such study carried out by Fawzi et al. [15] made use of Dynamic programming. Resener et al. [16] applied Mixed-integer linear programing (MILP). However, these analytical and numerical based methods are often computationally intensive as all possible combinations of DG sites are required to be evaluated to derive the optimum solution. This often limits the number of variables considered in most researches applying these methods. Also, due to the non-linear nature of the problem, the linear programming methods often perform poorly in finding the optimal solutions [17].
Metaheuristic population-based algorithms such as genetic algorithm, firefly, grey wolf, bee colony, particle swarm and ant colony algorithms have an advantage over the numerical and analytical based methods because it doesn't have to cover the entire possible search space in deriving the optimal solution. This helps to minimize the computational resources required and allows for consideration of a large number of variables and buses. Other advantages of these type of algorithms in solving optimization problems include their simplicity and flexibility. They however only derive a near-optimal solution rather than an exact one which is often within acceptable limits for most problems [18]. A metaheuristic cost-benefit approach for the siting and sizing of DGs to ensure the peak demand forecast is optimally met was presented by El-khattam et al. [19]. Kuri et al. [20] applied genetic algorithm for optimal positioning and sizing of DGs in a distribution network, they took into consideration the number of planning years in setting the optimization objectives and constraints. Golshan and Arefifar [21] used an algorithm based on Tabu Search (TS) in solving a comprehensive planning objective for determining the size, locations and mode of operation of Reactive Power Sources (RPSs) and Distributed Generation Resources (DGRs) in a power system. Gandomkar et al. [22], presented a hybrid of simulated annealing and genetic algorithm for solving the problem of DG siting in power distribution networks. Kalkhambkar et al. [23][24][25] had also implemented the grey wolf algorithm to minimize energy loss in energy storage and renewable distributed generations. Celli et al. [26] proposed a tool to plan and operate DGs for deferment of upgrades in transmission and distribution systems and so on. In the proposed algorithm, multi-objective programming and genetic algorithm (GA) were utilized in a finding the optimum location and number of distributed energy resources to be placed in a network. Similarly, a method for reduction in real power loss, improvement in voltage profile and release of substation in a power system based on analysis of voltage sensitivity index was presented by Gopiya et al. [27].
In simplifying the problem of DG location and sizing, most researchers in the literatures often either fix the number of DGs or the sizes. The proposed method in this paper eliminates this restriction by allowing the algorithm to heuristically search through the entire search space without restricting the possible number of DGs or the capacity of the DGs in deriving the optimum combination of the number, the site and size of DGs to minimize the real and reactive power loss in the system. Due to the relative novelty of ant colony optimization (ACO) algorithm compared to other population-based optimization algorithm and also because ACO was initially designed for discrete problems, not many researches on the suitability of ACO for DG allocation problem has been done. However, ACO has been chosen over other metaheuristic based approach for solving the DG allocation problem for its ability to derive the global optimum solution with minimal number of iterations and its better comparison with other similar algorithms [28,29]. This study helps to further explore the application of the ACO algorithm in discrete and continuous problem domain. The ACO algorithm will be applied to optimally site and size DGs in a power distribution system such that the total real power loss within the system is minimized while the system voltage profile is significantly enhanced. In order to demonstrate the suitability of the proposed algorithm, it was applied on the IEEE 30-Bus Distribution System.
The remainder of the paper is organized as follows: Sect. 2 provides the mathematical formulation of the problem of real power loss and voltage profile deviation of a power distribution system. The proposed Ant Colony Optimization algorithm for solving the DG optimization problem is presented in Sect. 3, 4. Section 5 presents the model of the case study power distribution system. Analysis and simulation results are presented in Sect. 6. Finally, the conclusion is drawn in Sect. 7.

Problem formulation
The total real power loss within the distribution system and the deviation of the bus voltages from the nominal voltage magnitude, is given by equations below [30]: where ℜ(•) denotes the real component of a complex number,y ij is the admittance of the ij th branch, N is the number of buses, V i and V j denotes the bus voltages, ij The objective function to be minimized to derive the best DG allocation for minimum power losses and enhanced voltage profile using ACO is formulated as: The inequality constraints of the minimization algorithm are outlined below: (1) where: nDG is the number of DG sources, P DG i is the real power output of the ith DG source while Q DG i is its reactive power output, S ij is the apparent power flow in the ijth branch and S ij max is the maximum apparent power flow in the ijth branch.

Optimization algorithm
In this study, the Ant Colony Optimization (ACO) is proposed for solving the DG allocation problem. ACO is an algorithm that is inspired by the biological behaviour of ants, which enable them to discover the shortest route between their nest and food. Ants do not communicate directly with one another, but they exchange information using what is known as pheromones [31]. When an ant traverses a path in search for food, it lays pheromones on the path, as more ants pass through the shorter path between the nest and the food per unit time relative to the longer path, the quantity of pheromone laid on the shorter path becomes higher than that on the longer path, thereby increasing the probability that an ant will choose the shorter path. The likelihood that an ant will choose a given path, therefore, depends on the quantity of pheromone laid on the path relative to other paths which in turn depends on the number of ants that had previously chosen the same path [28,32]. The ACO algorithm procedure comprises majorly of three steps: In the first step, the candidate variable values are randomly generated and the pheromone trail is initialized. Secondly, each ant uses the probabilistic state transition rule to find solution to the problem. Thirdly, the pheromone value at each edge is updated, first by evaporating a proportion of the pheromone at all the edges, thereafter the amount of pheromone on paths whose solutions are of high fitness values is reinforced.

Implementation of the algorithm
The Ant Colony Optimization (ACO) algorithm applied for this study was as presented in [33]. The main steps of the applied ACO algorithm are outlined below:

Generation of the candidate variable values
The solution variables are the real and reactive power generation capacities of the DGs and the bus number at which the generators will be connected. The number of DG is initialized to a certain number. Initially, the candidate variable values for the generators' real and reactive (8) min i ≤ i ≤ max i power capacity, which are real continuous variables, are randomly sampled subject to the equality and inequality constraints as: where (m + ϑ) is the number of initialized randomly generated solution for variable i, rand j i is a uniformly distributed random number between 0 and 1, i = 1, 2, …, n, while j = 1, 2, …, (m + ϑ). The initial buses for siting the DGs are uniformly randomly chosen from the total number of buses. The randomly generated initial solutions are evaluated and ranked according to their fitness values. The solution with the highest fitness value is captured as the initial global best solution.
Four sources are utilized in generating the variable values from which the ants choose during each iteration, they are: the best solutions derived from all previous iterations, which are also referred to as the global best solution; solutions generated from a random exploration process across the entire solution space; solutions derived from a dynamic exploitation of the solution space around the global best solution; and finally from the values selected by the ants in the last iteration of the algorithm. In each iteration of the ACO algorithm, the m number ants create m solutions of the problems which serves as candidate values for the next iteration of the algorithm. The best solution from the previous iteration is thereafter compared to the current global best solution in terms of their fitness value, and the global best solution is updated accordingly.

Dynamic exploitation process
This involves the searching of the solution space around x 0 = (x 0 1 , x 0 2 , ..., x 0 n ) , the global best solution, in the interval .., n . The value of the variables in the global best solution is either increased, decreased or left unchanged according to the following equations: where i ∈ (0, 1] and q ∈ [0, 1) are uniform random values, i = 1, 2, ..., n. Then the resulting solution x = (x 1 ,x 2 , ...,x n ) . If this process generates a solution with a better fitness value than the global best solution, the recorded global best solution is replaced with the derived solution. The process is repeated for ϑ number of times.

Ants' solution construction
In this step, m ants construct a solution from the generated candidate variable values. An ant k selects a variable with an index l (k) i for the ith variable from the candidate variable values according to (11).
where q is a uniform random number based on which an ant either selects a variable with the highest pheromone value among the ones generated in the previous iteration or randomly chooses an index L (k) l ∈ {0, 1, ..., m + g i } . The probability distribution for the random variable choice is given in (12).
The solution derived by each ant from the step above is denoted by

Pheromone update
The solutions generated by all the ants in step 3 are evaluated at the end of each iteration of the algorithm and are sorted in descending order of their fitness value. Thereafter, the amount of pheromones on a predetermined number of variable values are evaporated according to (13).
where is a real number between 0 and 1 which defines the applied rate of pheromone evaporation per iteration of the algorithm, T min is a predefined constant representing the minimum pheromone value. Also, to ensure convergence of the algorithm towards the global minima, the pheromone values of the best solutions' variable values are reinforced as where is the defined rate of pheromone reinforcement and it is any number between 0 and 1. The implementation flowchart for the ACO algorithm is shown in Fig. 1.

Test system
The proposed ACO algorithm was implemented in MATLAB R2016a and executed on an Intel® Pentium® N3540 2.16 GHz personal computer, while the power flow algorithm for evaluating the fitness value of each solution generated by the ants was implemented using Matpower® library in MATLAB. It was tested on an IEEE 30-bus system. The total loads of the 30-bus system are 189.2 MW and 107.2 MVar, with its bus and branch parameters given in [34]. A schematic diagram of the test feeder is shown in Fig. 2.
The following assumptions were made in formulating the proposed method: • Each bus within the system can take a maximum of one DG unit; • The total DG penetration within the distribution system will not exceed 100% of the total load; • The size of a single DG unit connected to a bus in the algorithm is allowed to can range from zero (no DG connected to the bus) to the total system load (100% DG penetration with only 1 DG source); • All the DGs have unlimited capability to supply real and reactive power according to the demand; • All the buses, except bus 1 which is the reference bus, has the capacity for DG connection; • Continuous values are assumed for the size of the DGs.

Validation
In order to validate the impact of results of the ACO algorithm, a model of the IEEE 30 bus system was developed using the Electrical Transient Analyzer Programme (ETAP) and load flow analysis of the system under the ACO derived optimal DG number, sites and sizes were carried out. The system was initially connected to the grid at bus 1 and the voltage profile at each of the 30 buses as well as the total system losses were recorded. Thereafter, 6 number generators were connected to the system, at first, the generators were connected at bus numbers 7, 8, 12, 19, 21 and 30 derived as the optimum sites for the DGs from the ACO analysis, the observed losses and system voltage profile are compared to the values obtained under grid power supply. To further revalidate the selection of the optimum site by the ACO algorithm, three additional case studies with the DGs sited at randomly selected buses with in the modelled system were evaluated. The IDs of the buses at which the DGs were sited under the various case studies are outlined in Table 1. The system voltage profile

ACO results
The voltage magnitude across all the buses in the distribution system before and after ACO optimized DGs siting are presented in Fig. 3. It shows a significant improvement in the system voltage profile. In the base case, with the system only fed from the grid via Bus 1, there is a decline in the per unit (pu) voltage magnitude at the buses as we progress along the system from bus 1 to bus 30. However, with the DGs sited at the optimum buses as derived using ACO, there is a significant improvement in the voltage profile of the system with all the buses approaching the nominal voltage. Before DG siting, the least voltage magnitude across all the buses was obtained at bus 26 with a per unit voltage of 0.656, however, after DG siting, this value increased by 47% to 0.965 per unit.
Similarly, as illustrated in Fig. 4 and Fig. 5, in the base case when the only power source of the network was situated at bus 1, we observed significant real and reactive power losses from bus 1 to 10, however, with the optimally sized and sited DGs, the losses reduced drastically across all the branches, with values approaching zero as close as possible.

Comparison of ACO with ETAP
The analysis of the distribution system voltage profile as derived from ETAP under the various distributed generation arrangement are illustrated in Fig. 6. The ACO results   Fig. 3 with the obtained results from ETAP shown in Fig. 6 for optimized site plot are in good correlation with an average difference of 0.52%. It is also shown in Fig. 6 that there is a significant improvement in the voltage profile when the system is fed by DGs compared to when it is only fed from the grid. Also, while the voltage at some of the buses under the other three case studies are better than when the DGs are sited at the ACO optimized sites, the best overall voltage profile is obtained at the ACO optimized sites. Table 2 presents the summary of results of the ACO optimized sites with the DG sites of randomly selected buses simulated in ETAP. It revealed that the total real and reactive power losses within the system reduced by 92% and 97%, respectively when the distribution system is fed by six DGs sited at the ACO optimized locations relative to the values under grid power supply. Also, the highest reduction in losses was obtained with the DGs sited at the ACO derived sites as against when they are sited at other locations within the distribution system, as obtained in Case 1, 2 and 3. This is also illustrated in Fig. 7 which shows a significant reduction in the total real and reactive power losses across the system when DGs are sited at various locations within the system under the four case studies, however, the least losses are obtained when the DGs are sited at the ACO derived optimum buses.
The results showed that the ACO based approach presented in this paper was able to achieve a 92% reduction in real power loss and a 97% reduction in reactive power loss within the distribution system. Compared to other reviewed literatures on the subject, this method resulted in the highest reduction in real and reactive power losses. A 67.01% reduction in losses was derived using Differential evolution technique [35]. Using a hybrid of GA and an analytical solution, a reduction in system losses of 83.74% was derived in [36]. Using closed-form analytical expressions, a maximum reduction in real power losses of 95% and 94% reduction in the reactive power losses, which is comparable with the result obtained in this study with the ACO method presented achieving this with a much smaller search space [12]. The superior performance of the method could be attributed to the fact that the algorithm is not excessively constrained with the optimum combination of DG power factors, the number, the capacity and the location of DGs heuristically derived by the algorithm coupled with the superior capacity of ACO to converge as the global minima compared with other metaheuristic algorithms.

Conclusion
The proper placement and sizing of DGs in a power distribution system plays a significant role in reducing the total real and reactive losses as well as improves the voltage profile within the system. In this paper, an optimization method for determining the optimal sites and capacities of distributed generation resources in a power distribution network using Ant Colony Optimization was presented. The problem consists of two types of variable, the DG placement problem which is discrete in nature and generators' sizing which is of a continuous variable type. The results obtained showed that the ACO based approach resulted in a significant reduction in the total real and reactive power losses as well as improvement in the voltage profile of the power distribution network as a result of the DGs. The impact of the optimally sized and sited DGs was also significantly higher than the figures achieved in most past research work on the same subject.
In the future, this research work will be extended by considering and modelling the different types of renewable energy sources in formulating the problem rather than using a generic DG source with the capacity to generate any amount of real power and supply or sink any amount of reactive power as presented in this work. Furthermore, more work still needs to be done to solve the planning and dispatch problem associated with a high penetration of renewable energy DGs in a power distribution system due to the variable nature of renewable energy sources.

Compliance with ethical standards
Conflict of interest The authors declare that they have no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.