Optimal Selection of Capacitors for Radial Distribution Systems to Improve Annual Cost Savings and Loss Minimization Using Optimization Techniques in Mat Lab Environment

M. Surendranatha Reddy*   M. Satish Kumar Reddy **  
*-** Assistant Professor, Department of Electrical and Electronics Engineering ,Vaagdevi Institute of Technology & Sciences,Andhra Pradesh, India.

Abstract

The problem of voltage deviation and power loss is mostly addressed in distribution system by growing domestic, industrial and commercial load day by day. To minimize these problems, effective planning of distribution system is required. This effective and reliable planning of distribution system is achieved by optimized placement of control device (capacitor) in distribution networks. The idea for optimal capacitor placement is determined with the location of the capacitor to be installed in the distribution network buses where power losses should be minimum and cost saving should be maximum. This paper presents a new and efficient approach for capacitor placement in radial distribution systems that determine the optimal locations and size of capacitor with an objective of improving the annual cost saving, the voltage profile and reduction of power loss. The solution methodology includes an algorithm that employs Plant Growth Simulation Algorithm (PGSA) and Genetic algorithm to estimate the optimal size of the capacitors at the optimal buses determined. The effectiveness of the proposed method is applied to 33 and 69 bus IEEE radial distribution systems. Result shows that after placement of capacitor in candidate buses the energy loss has been reduced.

Keywords :

Introduction

Dr N.Visali et.al, in their work “Loss Reduction in Radial Distribution System by Optimal Placement of Capacitor using Differential Evolution Method,” proposed optimization method (differential evolution) for capacitor placement in radial distribution system [1]. A.Kartikeya Sarma and K.Mahammad Rafi in their work “Optimal Selection of Capacitors for Radial Distribution Systems using Plant Growth Simulation Algorithm,” explained the importance of capacitor sizing and selection for radial distribution system to reduce the losses and cost savings [2]. The loss minimization in distribution systems has assumed greater significance recently since the trend towards distribution automation will require the most efficient operating scenario for economic viability variations. Studies have indicated that as much as 13% of total power generated is wasted in the form of losses in the distribution level. To reduce these losses, shunt capacitor banks are installed on distribution primary feeders. The advantages with the addition of shunt capacitors banks are to improve the power factor and feeder voltage profile. Power loss reduction increases the available capacity of the feeders. Therefore it is important to find out the optimal location and sizes of capacitors in the system to achieve the above mentioned objectives. Mr. Manish Gupta and Dr. Balwinder Singh Surjan in their work “Optimal sizing and Placement of Capacitors for Loss Minimization In 33-Bus Radial Distribution System Using Genetic Algorithm in MATLAB Environment,” proposed the method of genetic algorithm for capacitor placement and sizing in radial distribution system [3]. Since, the optimal capacitor placement is a complicated combinatorial optimization problem, many different optimization techniques and algorithms have been proposed in the past. Schmill [4] in his work “Optimum Size and Location of Shunt Capacitors on Distribution Feeders,” developed a basic theory of optimal capacitor placement. The author presented his well-known 2/3 rule for the placement of one capacitor assuming a uniform load and a uniform distribution feeder. H. Dura [5] in his work “Optimum Number Size of Shunt Capacitors in Radial Distribution Feeders,” considered the capacitor sizes as discrete variables and employed the dynamic programming to solve the problem. J.J. Grainger and S.H.Lee [6] in their work “Optimum Size and Location of Shunt Capacitors for Reduction of Losses on Distribution Feeders” developed a nonlinear programming based method in which capacitor location and capacity were expressed as continuous variables. J.J. Grainger and S. Civanlar [7] in their work “Volt/var control on Distribution systems with lateral branches using shunt capacitors as Voltage regulators” formulated the capacitor placement and voltage regulators problem and proposed decoupled solution methodology for general distribution system. M.E.Baran and F.F.Wu [8] in their work “Optimal Sizing of Capacitors Placed on a Radial Distribution System,” presented a method with mixed integer programming. Sundharajan and Pahwa [9] in their work “Optimal selection of capacitors for radial distribution systems using genetic algorithm,” proposed the genetic algorithm approach to determine the optimal placement of capacitors based on the mechanism of natural selection. In most of the methods mentioned above, the capacitors are often assumed as continuous variables. In distribution system to improve the voltage regulation it is required to minimize reactive power flow through the system. To overcome these difficulties the capacitor placement is essential in distribution system. The placement of capacitors in radial distribution systems, also provide power flow control, improving system stability, power factor correction, voltage profile management and losses minimization. In this paper, Capacitor Placement and Sizing is done by Loss Sensitivity Factors and Plant Growth Simulation Algorithm (PGSA) respectively. The loss sensitivity factor is able to predict which bus will have the biggest loss reduction when a capacitor is placed. Therefore, these sensitive buses can serve as candidate locations for the capacitor placement. PGSA is used for estimation of required level of shunt capacitive compensation to improve the voltage profile of the system. The proposed method is tested on 33 and 69 bus radial distribution systems and results are very promising. The advantages with the Plant Growth Simulation Algorithm (PGSA) is that it treats the objective function and constraints separately, which averts the trouble to determine the barrier factors and makes the increase/decrease of constraints convenient, and that it does not need any external parameters such as crossover rate, mutation rate, etc. It adopts a guiding search direction that changes dynamically as the change of the objective function. The remaining part of the paper is organized as follows: Section 1 gives the objective function; Section 2 Load Flows; Sections 3 explains system description, Sections 4 Solution Methodology gives brief description of the plant growth simulation algorithm and Genetic algorithm; Section 5 develops the test results and Section 6 gives conclusions.

1 Objective Function

The aim of the present work is to find out annual cost savings & loss minimization using optimization techniques. Therefore the objective function consists of two terms: 1) capacitor cost and 2) loss minimization.

1.1 Capacitor cost

Capacitor cost is divided in to two terms: cost installation cost and variable cost which is proportional to the rating of capacitors. Capacitor cost ( c ) is

(1)

Where Ci is the installation cost of the capacitor
Cv is the rate of capacitor per KVar
Qk is the rating of capacitor in particular bus let kth bus

1.2 Loss minimization

If I is the current in particular section let I in time duration T, then energy loss in section is

(2)

ELi is the energy loss in section i

Ii is the current in section i

Ri is the resistance in section i

T is the total time duration

The energy loss in time (T) of a feeder with n sections can be calculated as

(3)

The energy loss cost (ELC) can be calculated by

(4)

Ce is the energy rate

Therefore from two terms the objective function can be written as

(5)

S is the cost function for minimization

By minimizing the cost function, the net saving due to the reduction of energy losses for a given period of time including the cost of capacitors is

(6)

2. Load Flows

Distribution power flow is used to determine the steady state of bus powers, bus voltages and losses of the system. In the earlier methods of load flow studies only PQ buses are considered and flat voltage profile of all buses is assumed at the starting of the load flow. In this paper a new approach is used for the load flow to determine the steady state voltages and powers of each bus without flat voltage profile. If the system consists of PV buses the same load flow is used by modifying certain equations. As we know the equation for complex power

(7)

According to ohms law

(8)
(9)

If lm is the branch element in the RDS then Eqn(9) is represented as follows

(10)

The term is taken as new variable and represented as β. So the equation for β is represented as

(11)

2.1 For Load buses

(12)
(13)

The size of the element incidence matrix is “Nb X Nb.” It is constructed by using the following procedure.

(14)
(15)
(16)

The bus currents are calculated by using Eqn (7)

As the transmission losses are ignored for the first iteration and the disparity between the specified powers and the calculated powers are the part of the transmission loss, .TLlm r is the transmission loss part for 'lmth' element for  'rth'iteration.

(17)
(18)
(19)

Where 'r' is the iteration count.

The criteria is continuing up to this load flow meeting of the convergence which is the tolerance of 0.0001.

2.2 Handling to the voltage controlled buses

When the power is supplied at multiple ends of the system, other supplying buses except slack bus are treated as voltage controlled buses. Equation (5) is modified for the mth voltage controlled bus

(20)
(21)

where Xlm is the real value of the βlm

Ylm is the imaginary value of the βlm

(22)
(23)

The trigonometric equation is to be solved to get the phase angle φ each PV bus m and the reactive power can be updated as

(24)

3. System Description

IEEE 33-bus system is considered to find the optimal location of capacitors PGSA & GA optimization methods. The system has 4 feeders which supply to 33 load centres. The system has 33 buses, 32 lines, 12.66KV as a base voltage and 100MVA as a base MVA as in Figure.1. IEEE 69-bus system is considered to find the optimal location of capacitors PGSA & GA optimization methods. The system has 69 buses, 68 lines. The system has 12.66KV as a base voltage and 100MVA as a base MVA as in Figure 2. V Solution Methodology

Figure 1. Single line diagram of 33 bus Radial Distribution System

Figure 2. Single Line Diagram of 69 bus Radial Distribution System

3.1 Plant Growth Simulation Algorithm

The plant growth simulation algorithm[10] is based on the plant growth process, where a plant grows a trunk from its root; some branches will grow from the nodes on the trunk; and then some new branches will grow from the nodes on the branches. Such process is repeated, until a plant is formed. Based on an analogy with the plant growth process, an algorithm can be specified where the system is to be optimized First one of the “grows” begins at the root of a plant and then the second “grows” begins at the branches continually until the optimal solution is found. By simulating the growth process of plant phototropism, a probability model is established. In the model, a function g(Y) is introduced for describing the environment of the node Y on a plant. The smaller the value of g(Y), the better is the environment of the node for growing a new branch. The outline of the model is as follows; A plant grows a trunk M, from its root B0. assuming there are k nodes BM1 , BM2 , Bm ,.....BMk that have been the better environment than the root on the trunk M, which means the function g(Y) of the nodes and satisfies g(BMi) < g(Bo) then morphactin concentrations CM1 ,CM2,….,CMk of nodes BM1 , BM2 , BM3 , …..,BMk are calculated using

(25)

The significance of the above equation (X) is that the morphactin concentration of a node is not only dependent on its environmental information but also depends on the environmental information of the other nodes in the plant, which really describes the relationship between the morphactin concentration and the environment.

A complete algorithm for the proposed method of the capacitor placement is given below[10]:

3.2 Genetic Algorithm

The genetic algorithm is a global search technique for solving optimization problems, which is based on the theory of natural selection, the process that drives biological evolution. Genetic Algorithm proved to be a very effective and efficient tool for operation and control of the power system. Among the various application schemes GA based control scheme has played a significant role in AGC.

3.2.1 Steps involved in determining the optimal parameters of the controller using genetic algorithm are given below.

Step 1. Start: Create random population of n Chromosomes

Step 2. Fitness: Evaluate fitness of each chromosome in the population

Step 3. New population

Selection: Based on fitness function

Step 4. Replace: old with new population and the new generation

Test: Test for problem criterion

Loop: Continue step II-V until criterion is satisfied.

3.2.2 Solution algorithm for capacitor placement [11]

The solution methodology is given by following steps

Step 1: Read system data's (Bus data's and Line data's).

Step 2: Calculate the Y-bus and perform the load flow analysis and find out the voltage magnitude and power flow in branches.

Step 3: Performs the TS Initialize Tabu list, Tabu size and initial solution, and find the optimal location of capacitor placements.

Step 4: Place the capacitor at appropriate location which was determined in the previous.

Step 5: Perform the load flow analysis and compare the results before and after placement of capacitors.

4. Test Results

The proposed method has been programmed using MATLAB. The effectiveness of the proposed method for loss reduction by capacitor placement is tested on 33 and 69 bus radial distribution systems. The results obtained in the proposed method are explained in the following sections.

4.1 33 Bus System

The first test case for the proposed method is a 33-bus radial distribution system single line diagram as shown in in Figure 1. The line and load data of the feeders are taken from the reference [12]. The base values of the system are taken as 12.66 kV and 100MVA. Figure 3 shows the flowchart for Capacitor Placement. Tables 1, 2 and Figures 4, 5 depict results of PGSA. The values in Table 5 show the voltage profile values before and after capacitor placement. The values in Table 6 explain the total power loss, total reactive power loss, cost of saving after capacitor placement and elapsed time values. In Figure 8 the Graph shows the comparison of voltage profile values before and after capacitor placement. The red colour curve of the graph shows the voltage profile before the capacitor placement. The blue colour curve of the graph shows the voltage profile after the capacitor placement. In Figure 9 the Graph shows the best cost reduction with respect to the number of iterations. Tables 5, 6 and Figures 8, 9 depict results of GA.

Figure 3. Genetic Flowchart of Solution Algorithm for Capacitor Placement

4.2 69 Bus System

The next test case for the proposed method is a 69-bus radial distribution system single line diagram is shown in Figure 2. The line data and feeder characteristics are taken from reference [13]. Tables 3, 4 and Figures 6, 7 depict results of PGSA. The values in Table 7 show the voltage profile values before and after capacitor placement. The values in Table 8 explain the total power loss, total reactive power loss, cost of saving after capacitor placement and elapsed time values. In Figure 10 the graph shows the comparison of voltage profile values before and after capacitor placement. The red colour curve of the graph shows the voltage profile before the capacitor placement. The blue colour curve of the graph shows the voltage profile after the capacitor placement. In Figure 11 the graph shows the best cost reduction with respect to number of iterations. Tables 7, 8 and Figures 10, 11 depict results of GA.

Figure 4. Voltage Profile of 33 Bus System Before and After Compensation

Figure 5. Best Cost Reduction of 33 Bus System

Table 1.Voltage Profile of 33 Bus Radial Distribution System

Table 2. Summary of results for 33 bus radial distribution system

Table 3. Voltage profile of 69 Bus Radial Distribution System

Table 4. Summary of results for 69 Bus Radial Distribution System

Table 5. Voltage profile of 33 Bus Radial Distribution

Table 6. Summary of results for 33 Bus Radial Distribution System

Figure 6. Voltage Profile of 69 Bus System Before and After Compensation

Figure 7. Best Cost Reduction of 69 Bus System

Figure 8. Voltage Profile of 33 Bus System Before and After Compensation

Table 7. Voltage Profile of 69 Bus Radial Distribution System

Table 8. Summary of results for 69 Bus Radial Distribution System

Table 9. Summary of results for 33 Bus Radial Distribution System For PGSA and GA

Table 10. Summary of results for 69 Bus Radial Distribution System for PGSA and GA

Figure 9. Best cost reduction of 33 Bus System

Figure 10. Voltage Profile of 69 Bus System Before and After Compensation

Figure 11. Best Cost Reduction of 69 Bus System

Table 9 shows the values for real power loss, reactive power loss, time for iteration and cost saved after the capacitor placement for both the optimization methods GA and PGSA. For the present case study analysis GA shows the better result compared to PGSA for cost saving in 33 bus radial distribution system. The values in Table 10 show the values for real power loss, reactive power loss, time for iteration and cost saved after the capacitor placement for both of the optimization methods GA and PGSA. For the present case study analysis GA shows the better result when compared to PGSA for cost saving in 69 bus radial distribution system.

Conclusion

This work has been successfully applied in PGSA & GA optimization methods for determining the optimal location of capacitor placement in 33 & 69 bus radial distribution systems for minimum value of objective function. Here the power losses and cost savings are reduced by the placement of capacitors. By placing the capacitor the voltage level can always be in allowable range. So the power losses in the systems are reduced to a great extent by the optimal placement of capacitors in the distribution system.

References

[1]. Dr N Visali, M Satishkumar Reddy and M Surendranatha Reddy (2013). “Loss Reduction in Radial Distribution System by Optimal Placement of Capacitor Using Differential Evolution Method,” International Journal of Emerging Trends in Engineering Research, Vol 1, No 2, October , pp-53-57.
[2]. A.Kartikeya Sarma and K.Mahammad Rafi (2011). “Optimal Selection of Capacitors for Radial Distribution Systems Using Plant Growth Simulation Algorithm”, International Journal of Advanced Science and Technology, Vol. 30, May.
[3]. Mr. Manish Gupta and Dr. Balwinder Singh Surjan (2012). in their work “Optimal sizing and Placement of Capacitors for Loss Minimization In 33-Bus Radial Distribution System Using Genetic Algorithm in MATLAB Environment,” International Journal of Advanced Research in Computer Engineering & Technology, Volume 1, Issue 8, October.
[4]. J. V. Schmill, (1965). "Optimum Size and Location of Shunt Capacitors on Distribution Feeders," IEEE Transactions on Power Apparatus and Systems, Vol. 84, pp. 825-832, September.
[5]. H. Dura (1968). "Optimum Number Size of Shunt Capacitors in Radial Distribution Feeders: A Dynamic Programming Approach", IEEE Trans. Power Apparatus and Systems, Vol. 87, pp. 1769-1774, Sep.
[6]. J. J. Grainger and S. H. Lee, (1981). “Optimum Size and Location of Shunt Capacitors for Reduction of Losses on Distribution Feeders,” IEEE Trans. on Power Apparatus and Systems, Vol. 100, No. 3, pp.11051118, March.
[7]. J.J. Grainger and S. Civanlar, (1985). “Volt/var control on Distribution systems with lateral branches using shunt capacitors as Voltage regulators-part I, II and III,” IEEE Trans. Power Apparatus and systems, Vol. 104, No. 11, pp. 3278-3297, Nov.
[8]. M. E Baran and F. F. Wu, (1989). “Optimal Sizing of Capacitors Placed on a Radial Distribution System”, IEEE Trans. Power Delivery, Vol. No.1, pp. 1105-1117, Jan.
[9]. Sundharajan and A. Pahwa, (1994). “Optimal selection of capacitors for radial distribution systems using genetic algorithm,” IEEE Trans. Power Systems, Vol. 9, No.3, pp.1499-1507, Aug.
[10]. Chun Wang, H. Z. Chengand L. Z Yao, (2008). “Optimization of network reconfiguration in large distribution systems using plant growth simulation algorithm,” DRPT 2008 Conference, Nanjing, China, pp. 771-774, 6-9, April.
[11]. Mrs. Asha Gaikwad, Dr. Rakesh Ranjan, Dr. L. D. Arya (2010). ”Capacitor Placement for Loss Reduction of Reconfigured Radial Distribution Systems by Depth First Search Algorithm” International Journal of Advanced Engineering & Applications, Jan.
[12]. M.Chis, M. M. A. Salama and S. Jayaram, (1997). “Capacitor Placement in distribution system using heuristic search strategies,” IEE Proc-Gener, Transm, Distrib, Vol, 144, No.3, pp. 225-230, May.
[13]. Das et al, (1995). “Simple and efficient method for load flow solution of radial distribution network,” Electric Power and Energy Systems, Vol. 17, No.5, pp.335-346,