Optimal Placement and Sizing of DG in Radial Distribution Networks Using SFLA

Optimal placement and sizing of DG in distribution network is an optimization problem with continuous and discrete variables. Many researchers have used evolutionary methods for finding the optimal DG placement. This paper proposes a shuffled frog leaping algorithm (SFLA) for optimal placement and sizing of distributed generation (DG) in radial distribution system to minimize the total real power loss and to improve the voltage profile. The SFLA is a meta-heuristic search method inspired from the memetic evolution of a group of frogs when seeking for food. It consists of a frog leaping rule for local search and a memetic shuffling rule for global information exchange. The proposed SFL algorithm is used to determine optimal sizes and locations of multi-DGs. Test results indicate that SFLA method can obtain better results than the simple heuristic search method on the 33-bus radial distribution systems. Moreover, voltage profile improvement and branch current reduction are obtained.


Introduction
One of the most important motivation for the studies on integration of distributed resources to the grid is the exploitation of the renewable resources such as; hydro, wind, solar, geothermal, biomass and ocean energy, which are naturally scattered around the country and also smaller in size. Accordingly, these resources can only be tapped through integration to the distribution system by means of Distributed Generation. Distributed Generation (DG), which generally consists of various types of renewable resources, can be defined as electric power generation within distribution networks or on the customer side of the system [1].DG affects the flow of power and voltage conditions on the system equipment. These impacts may manifest themselves either positively or negatively depending on the distribution system operating conditions and the DG characteristics. Positive impacts are generally called 'system support benefits', and include voltage support and improved power quality; loss reduction; transmission and distribution capacity release; improved utility system reliability. On account of achieving above benefits, the DG must be reliable, dispatchable, of the proper size and at the proper locations [2], [3]. Distributed Generation (DG) is a small generator spotted throughout a power system network, providing the electricity locally to load customers. DG can be an alternative for industrial, commercial and residential applications. DG makes use of the latest modern technology which is efficient, reliable, and simple enough so that it can compete with traditional large generators in some areas [4], [5]. In fact, three types of DG are considered as follows: DG is capable of supplying only real power; DG is capable of supplying only reactive power; DG is capable of supplying real power but consuming proportionately reactive power. The methodology is proposed for optimal placement and sizing of only first type DG. Moreover, the heuristic search requires exhaustive search for all possible locations which may not be applicable to more than one DG. Therefore, in this paper, SFLA method is proposed to determine the optimal location and sizes of multi-DGs to minimize the total real power loss and improve the voltage profile of the distribution systems. The organization of this paper is as follows. Section 2 addresses the Load flow in radial distribution network. The problem formulation is in Section 3. The SFL algorithm is presented in Section 4. A SFLA computation procedure for the optimal placement and sizing of distributed generation problem is given in Section 5. Simulation result on the test systems are illustrated in Section 6.

Load Flow
On account of the some inherent features of distribution systems such as radial structure, unbalanced distributed loads, large number of nodes, a wide range of R/X ratios, the conventional techniques have been developed for transmis-sion systems generally fail on the determination of optimum size and location of distributed generations. In this study, the proposed methodology is based on the equivalent current injection that uses the Bus-Injection to Branch-Current (BIBC) and Branch-Current to Bus-Voltage (BCBV) matrices which were developed based on the topological structure of the distribution systems and is implemented for the load flow analysis of the distribution systems. The methodology proposed here requires only one base case load flow to determine the optimum size and location of DG. Detailed description of BIBC and BCBV matrix's building algorithm is omitted due to the lack of space.

Problem Formulation
The problem is to determine allocation and size of the DGs which minimizes the distribution power losses and to improve the voltage profile for a fixed number of DGs and specific total capacity of the DGs.

Objective Function
In this paper the objective function for the optimal placement and sizing of DG in distribution network problem is to minimize the real power losses and improve the voltage profile, which is calculated as follows: X=l 1 ,l 2 ,…,l bus ,x 1 ,x 2 ,…,x power limit (3) Where, R i and I i are resistance and actual current of the ith branch, respectively. N br is the number of the branches. V i is the magnitude of bus voltage i. V i,ref is the magnitude of voltage of slack bus .X is the vector of control variables. l is the number of DG location candidates. x is the number of capacity types of DGs. The objective function of the placement and sizing is to minimize the real power loss and improve the voltage profile. Mathematically, the objective function can be written as: W 1 and W 2 are weight factor.

Constraints
The constraints are listed as follows: a) Distribution line absolute power limits: Line are the absolute power and its corresponding maximum allowable value flowing over the distribution line between the nodes i and j, respectively.
Where N sc is total number of sections, P L is the real power loss in the system, P DGi is the real power generation at bus i, P Di is the power demand at bus i.

Shuffled Frog Leaping Algrithm (SFLA)
The SFL algorithm, in essence, combines the benefits of the genetic-based MAs and the social behavior-based PSO algorithms. In the SFL, the population consists of a set of frogs (solutions) that is partitioned into subsets referred to as memeplexes. The different memeplexes are considered as different cultures of frogs, each performing a local search. Within each memeplex, the individual frogs hold ideas, that can be influenced by the ideas of other frogs, and evolve through a process of memetic evolution. After a defined number of memetic evolution steps, ideas are passed among memeplexes in a shuffling process. The local search and the shuffling processes continue until defined convergence criteria are satisfied [6].As described in the pseudocode of Appendix A, an initial population of P frogs is created randomly. For S-dimensional problems (S variables), a frog i is represented as X i = (x i1 , x i2 ,…, x iS ). Afterwards, the frogs are sorted in a descending order according to their fitness. Then, the entire population is divided into m memeplexes, each containing n frogs (P m×n). In this process, the first frog goes to the first memeplex, the second frog goes to the second memeplex, frog m goes to the mth memeplex, and frog m+1 goes back to the first memeplex, etc. Within each memeplex, the frogs with the best and the worst fitnesses are identified as X b and X w , respectively. Also, the frog with the global best fitness is identified as X g . Then, a process similar to PSO is applied to improve only the frog with the worst fitness (not all frogs) in each cycle. Accordingly, the position of the frog with the worst fitness is adjusted as follows: Change in frog position(D i ) = rand( )× ( X b X w ) (10) New position X w = current position Xw + Di (11) D max ≥Di≥-Dmax Where rand( ) is a random number between 0 and 1 and D max is the maximum allowed change in a frog's position. If this process produces a better solution, it replaces the worst frog, otherwise the calculations in (10) and (11) are repeated but with respect to the global best frog (X g replaces X b ). If no improvement becomes possible in this case, then a new solution is randomly generated to replace that frog. The calculations then continue for a specific number of iterations [6]. Accordingly, the main parameters of SFL are: number of frogs P, number of memeplexes, number of generation for each memeplex before shuffling, number of shuffling iterations, and maximum step size.

SFLA Procedure
The SFLA-based approach for solving the optimal placement and sizing of distributed generation problem to minimize the loss and improve the voltage profile takes the following steps: In SFLA, each possible solution X i = (x i1 , x i2 ,…, x iS ) that in this paper X i =l 1 ,l 2 ,…,l bus ,x 1 ,x 2 ,…,x power limit Where, l is the number of DG location candidates and x is the number of capacity types of DGs are is considered as a frog. The steps of the algorithm are as follows: Step 1: Create an initial population of P frogs generated randomly. SFLA_Population =[X 1 ,X 2 ,…,X p ] p×n Where, P=m×n, N is the number of DG, m is the number of memplexes and n is the number of frogs in memplex.
Step 2: Sort the population increasingly and divide the frogs into m memplexes each holding n frogs such that P=m×n. The division is done with the first frog going to the first memplex, second one going to the second memplex, the mth frog to the mth memplex and the m+lth frog back to the first memplex. Fig.1 illustrates this memeplex partitioning process. Step 3: Within each constructed memeplex, the frogs are infected by other frogs' ideas; hence they experience a memetic evolution. Memetic evolution improves the quality of the meme of an individual and enhances the individual frog's performance towards a goal. Below are details of memetic evolutions for each memeplex: Step 3-1: Set m 1 =0 where m 1 counts the number of me-meplexes and will be compared with the total number of memeplexes m. Set y 1 =0 where y 1 counts the number of evolutionary steps and will be compared with the maximum number of steps (y max ), to be completed within each memeplex.
Step 3-2: Set m 1 =m 1 +1 Step 3-3: Set y 1 =y 1 +1 Step 3-4: For each memplex, the frogs with the best fitness and worst fitness are identified as X w and X b respectively. Also the frog with the global best fitness X g is identified, and then the position of the worst frog X w for the memplex is adjusted such as (10) and (11). Fig.2 demonstrates the original frog leaping rule. If the evolutions produce a better frog (solution), it replaces the older frog, otherwise X b is replaced by X g in (10) and the process is repeated. If no improvement becomes possible in this case a random frog is generated which replaces the old frog.
Step 4: Check the convergence. If the convergence criteria are satisfied stop, otherwise consider the new population as the initial population and return to the step2. The best solution found in the search process is considered as the output results of the algorithm. The flowchart of the SFLA is illustrated in Fig. 3.

Simulation Results
The test system for the case study is radial distribution system with IEEE 33 buses as shown in Figure 4.The total loads for this test system are 3.72 MW and 2.3 MVR. The original total real power loss and reactive power loss in the system are 227.6934 KW and 150.1784 KVAR respectively. The substation voltage is 12.66 KV and the base of power is 10.00MVA.The current carrying capacity of branch No.1-9 is 400 A, and the other remaining branches including the tie lines are 200A. The minimum and maximum voltages are set at 0.95 and 1.05 p.u. respectively. The load data are given in Table A1 and branch data is in Table A2 [7,8]. For SFLA parameters, population size is 50. The maximum iteration for SFL algorithm is 5.The number of memplexes is 5. The number of frogs in memplex is 10.The number of iterations in each memeplexes is 5. The total number of algorithm iterations is 5.The maximum number of DG is 3. The maximum real power of DG is 1200KW. The improvement in the voltage profile after optimally placing the DGs is shown in Figure 5. Without DG, the bus no. 18 has the lowest voltage of 0.8889 p.u. and the bus voltage has improved to 0.9687 p.u. after installing DG. For the 33 bus systems, in table 1, the SFLA can obtain the same optimal size and location. For the 33bus system, first type DG can reduce the total real power loss by 48.09%. For two and three types DGs, they can further reduce the real power loss by 67.88% and 73.66%, respectively.

Conclusions
In this paper, a shuffled frog leaping algorithm (SFLA) for optimal placement and sizing of multi-DGs is efficiently minimizing the total real power loss and improve the voltage profile satisfying transmission line limits and constraints. The methodology is fast and accurate in determining the sizes and locations. DG regulating bus voltage will be considered in future research work.

Appendix
Pseudocode for a SFL procedure: Begin; Generate random population of P solutions (frogs); For each individual i∈P calculate fitness (i); Sort the population P in descending order of their fitness; Divide P into m memeplexes; For each memeplex determine the best and worst frogs; Improve the worst frog position using equations (10) or (11); Repeat for a specific number of iterations; End; Combine the evolved memeplexes; Sort the population P in descending order of their fitness; Check if termination = true; End;