World Academy of Science, Engineering and Technology

International Journal of Electrical and Computer Engineering

Vol:7, No:10, 2013

Reliability Improvement with Optimal Placement of

Distributed Generation in Distribution System
N. Rugthaicharoencheep, T. Langtharthong

Distribution system reliability assessment can, in general,

Abstract—This paper presents the optimal placement and sizing be divided into the two basic tasks of assessing past
of distributed generation (DG) in a distribution system. The problem performance and predicting future performance. Predicting
is to reliability improvement of distribution system with distributed reliability performance is usually concerned with the supply
generations. The technique employed to solve the minimization
adequacy at the customer load points [5]. The conventional
problem is based on a developed Tabu search algorithm and
reliability worth analysis. The developed methodology is tested with approach to teaching distribution system reliability evaluation,
a distribution system of Roy Billinton Test System (RBTS) bus 2. It in either a university or industry based setting, is to use the basic
can be seen from the case study that distributed generation can reduce analytical equations to calculate load point failure rates, average
the customer interruption cost and therefore improve the reliability of outage durations and average annual outage times.
the system. It is expected that our proposed method will be utilized The technique employed to solve the minimization problem
Open Science Index, Electrical and Computer Engineering Vol:7, No:10, 2013 waset.org/Publication/17067

effectively for distribution system operator.

is based on a developed Tabu search algorithm and reliability
worth analysis. The Tabu algorithm systematically searches
Keywords—Distributed generation Optimization technique
Reliability improvement, Distribution system.
solutions expressed in forms of the location and size of DGs.
The solution obtained will then be passed to reliability worth
I. INTRODUCTION analysis to evaluate the quality of the solution. The process is
repeated until the best solution has been found. The developed
T HE distribution system is an important part that provides
the final link between the utility and customers. In
practice, most distribution systems have a single-circuit main
methodology is tested with a distribution system of Roy
Billinton Test System (RBTS) bus 2.
feeder and are configured radially. The radial distribution II. TABU SEARCH
system is widely used because of its simple design, generally
low cost and supportive protection scheme. This configuration Tabu search is a meta-heuristic that guides a local heuristic
suggests that all components between a load point and the search strategy to explore the solution space beyond local
supply point should completely operated and therefore poor optimality [6]. The basic idea behind the search is a move
reliability can be expected as the failure of any single from a current solution to its neighborhood by effectively
component causes the load points to be disconnected. An ideal utilizing a memory to provide an efficient search for
alternative on electric distributions to electric users is the optimality. The memory is called “Tabu list”, which stores
installation of a small sized generator or commonly known as attributes of solutions. In the search process, the solutions are
distributed generation (DG) [1]. in the Tabu list cannot be a candidate of the next iteration. As
DG is a small sized generator connected in parallel with the a result, it helps inhibit choosing the same solution many times
distribution system. DG is expected to play an increasing role and avoid being trapped into cycling of the solutions [7]. The
in emerging electric power systems. Studies have predicted quality of a move in solution space is assessed by aspiration
that DG will be a significant percentage of all new generation criteria that provide a mechanism for overriding the Tabu list.
going on line. There are several different types of resources Aspiration criteria are analogous to a fitness function of the
and technologies that can be used for DG such as wind, solar, genetic algorithm and the Bolzman function in the simulated
fuel cells, hydrogen, and biomass. DG can result in a network annealing as shown in Fig. 1.
operation and planning practices with economic implications. In the search process, a move to the best solution in the
The benefits of DG are classified into two groups: technical neighborhood, although its quality is worse than the current
and economics [2]. For example, loss, voltage profile, solution, is allowed. This strategy helps escape from local
reliability of supply, maintenance costs, and network optimal and explore wider in the search space.
connection reinforcement costs can be affected by the A Tabu list includes recently selected solutions that are
connection of DG to the distribution system [3], [4]. Electric forbidden to prevent cycling. If the move is present in the
utilities, therefore, can benefit from the installation of DG. Tabu list, it is accepted only if it has a better aspiration level
than the minimal level so far [8]. Fig. 2 shows the main
concept of a search direction in Tabu search.
N.Rugthaicharoencheep is with Department of Electrical Engineering, Faculty
of Engineering, Rajamangala University of Technology Phra Nakhon, Bangkok,
Thailand (e-mail: nattachote.r@rmutp.ac.th).
T. Langtharthong is with Department of Electrical Engineering, Faculty of
Engineering, Rajamangala University of Technology Phra Nakhon, Bangkok,
Thailand (e-mail: thong.l@rmutp.ac.th).

International Scholarly and Scientific Research & Innovation 7(10) 2013 1308 ISNI:0000000091950263
 World Academy of Science, Engineering and Technology
International Journal of Electrical and Computer Engineering
Vol:7, No:10, 2013

Average service unavailability index (ASUI)

ASUI = 1 − ASAI =
∑Ui Ni (5)
∑ Ni × 8760
Energy not supplied index (ENS)

ENS = ∑ La(i )U i (6)

Average energy not supplied index (AENS)

AENS =
∑ La(i)Ui
Fig. 1 Search direction of Tabu search (7)
∑ Ni
where
λi = failure rate of contingency i
Ni =
Open Science Index, Electrical and Computer Engineering Vol:7, No:10, 2013 waset.org/Publication/17067

total number of load points i

Ui = annual outage duration i
La (i ) = failure rate of contingency i
A basic approach to quantifying the worth of electric service
reliability is to estimate customer interruption costs due to
Fig. 2 Search direction of Tabu search electric power supply interruptions. One convenient way is an
interpretation of customer interruption costs in terms of
III. DISTRIBUTION SYSTEM RELIABILITY customer damage functions. The customer outage cost (
The basic distribution system reliability indices at a load ECOST ) is calculated from reliability indices of the load
point are average failure rate λ, average outage duration r, and point and customer damage function [9]-[10].
annual outage duration U. With these three basic load point
indices, the following system reliability indices can be IV. PROBLEM FORMULATION
calculated [8]. The objective is to minimize the customer outage cost that
can be written as follows:
Average interruption frequency index (SAIFI)
nh ni

λi N i Minimize ECOST = ∑ ∑ (La(i)Chi rh λh ) (8)

SAIFI = ∑ Ni
(1) h =1 i =1

where
System average interruption duration index (SAIDI) La (i ) = average load connected to load point i

C hi = outage cost ($/kW) of customer due to contingency h

Ui Ni
SAIDI = ∑ Ni
(2) λh = failure rate of contingency h
rh = average outage time of contingency h
nh = number of contingencies h
Customer average interruption duration index (CAIDI) ni = total number of load points i
Constraints:
UN
CAIDI = ∑ λiN i (3) 1) Power flow equations:
i i
NB

Average service availability index (ASAI)

Pk = ∑ YikViVk cos(θik + δk − δi ) (9)
i =1

ASAI =
∑ Ni × 8760 − ∑Ui Ni (4)
NB

∑ Ni × 8760 Qk = −∑ YikViVk sin(θik + δk − δi ) (10)

i =1

International Scholarly and Scientific Research & Innovation 7(10) 2013 1309 ISNI:0000000091950263
 World Academy of Science, Engineering and Technology
International Journal of Electrical and Computer Engineering
Vol:7, No:10, 2013

2) Voltage of each bus k must be within specified limits: iteration and iteration index m=1.
Vk min ≤ Vk ≤ Vk max (11) Step 2: Let the initial solution obtained in step 1 be the
current solution and the best solution: Sbest = S0, and
3) Current transfer capability of feeders l : Scurrent = S0.
Step 3: Perform a power flow analysis to determine whether the
Il ≤ Ilmax ; l ∈ {1, 2,..., Nl } (12) current solution satisfies the constraints defined in (9)
and (10). A penalty factor is applied for constraint
4) Maximum number of DGs (nDG ) to be installed: violation.
NB Step 4: Calculate ECOST using (8) with consideration of load
point restoration.
∑ e jk ≤ nDG j ∈ {1, 2,..., NC } (13)
Step 5: Calculate the aspiration level of Sbest : fbest = f(Sbest).
k =1
The aspiration level is the sum of ECOST and a
5) Maximum installed capacity of DGs: penalty function.
N B NC Step 6: Generate a set of solutions in the neighborhood of
∑ ∑ C j ⋅ e jk ≤G (14) Scurrent. This set of solutions is designated as Sneighbor.
k =1 j =1 Step7: Calculate the aspiration level for each member of
6) Decision variables for the installation of a DG: Sneighbor , and choose the one that has the highest
aspiration level, Sneighbor_best.
⎧⎪ 0 if the DG is not installed at bus k ⎫⎪⎪
Open Science Index, Electrical and Computer Engineering Vol:7, No:10, 2013 waset.org/Publication/17067

⎪⎪ Step 8: Check whether the attribute of the solution obtained

⎪⎪
e jk = ⎪⎨ 1 if the DG is installed at bus k ⎬ (15) in step 7 is in the Tabu list. If yes, go to step 9, or
⎪⎪ ⎪⎪ else Scurrent = Sneighbor_best and go to step 10.
⎪⎪ with the capacity at step j ⎪⎪
⎩ ⎭ Step 9: Accept Sneighbor_best if has a better aspiration level than
7) Only one DG can be installed at one position fbest and set Scurrent = Sneighbor_best , or else select a next-best
solution that is not in the Tabu list to become the
NC
current solution.
∑ e jk ≤ 1 k ∈ {1, 2,..., N B } (16)
Step 10: Update the Tabu list and set m = m+1.
j =1
Step 11: Repeat steps 6 to 10 until the specified maximum
where
Pk = iteration has been reached and report the best
power active power at bus k
solution.
Qk = power reactive power at bus k
Summing up, the performance of Tabu search depends on a
NB = number of buses proper choice of the neighbor of a solution, on the number of
Yik = element ( i , k ) in bus admittance matrix iterations for which a move is kept as Tabu, on the best
θik = angle of Yik combination of short- and long-term memory and on the best
δk = voltage angle at bus k balances of intensification and diversification mechanism [10].
Nl = number of feeders The solution algorithm for the problem is described step by
min = step as shown in Fig. 3.
V k minimum voltage at bus k
Vkmax = maximum voltage at bus k VI. CASE STUDY
Il = current flow in feeder l The developed Tabu search algorithm is tested with a
= distribution system of RBTS bus 2 [11] to minimize the customer
I lmax maximum current capability of feeder l
outage cost. There are 4 feeders and 22 load points. The peak
NB = number of buses loading level of bus 2 is 20 MW. The configuration of the system
NC = number of capacity steps of a DG is shown in Fig. 4. The maximum iteration for Tabu search is
Cj = capacity at step j of DG 100. The minimum and maximum voltages for each bus are 0.95 p.u.
G = maximum total installed capacity and 1.05 p.u., respectively. The sizes of DGs are 100 kW-1,500
NC = number of capacity steps of a DG kW. The failure of a transformer is recovered by repair. All
Cj = protective devices and DGs are assumed to be fully reliable.
capacity at step j of DG
Three cases are investigated in Table I. The results from the case
e jk = decision variable for installation of a DG at bus k study are shown in Tables II and III.
with the capacity at step j

V. SOLUTION ALGORITHM
The solution algorithm for the problem is described step by
step as follows:
Step 1: Randomly select a feasible solution from the search
space: S0 ∈ Ω. Set the size of a Tabu list, maximum

International Scholarly and Scientific Research & Innovation 7(10) 2013 1310 ISNI:0000000091950263
 World Academy of Science, Engineering and Technology
International Journal of Electrical and Computer Engineering
Vol:7, No:10, 2013

TABLE I
CASE STUDY FOR RELIABILITY ANALYSIS
Maximum number Total installed capacity,
Case
of DGs, nDG (unit) G (kW)
1 - -
2 1 ≤ 1000
3 3 ≤ 2000
4 3 ≤ 3000
5 4 ≤ 4000

TABLE II
OPTIMAL PLACEMENT AND SIZING OF DGS
Location of DG Capacity of DG installed Total capacity
Case
(bus) (kW) of DG (kW)
1 - - -
2 5 500 0500
3 3, 5, 11 100, 500, 500 1100
4 10, 11, 15 1000, 600, 1200 2800
5 5, 10, 12, 15 600, 1200, 200, 1200 3200
Open Science Index, Electrical and Computer Engineering Vol:7, No:10, 2013 waset.org/Publication/17067

All the cases have the same SAIFI because this index
depends only on the reliability of components (e.g., lines,
transformers) and is not affected distributed generations to be
installed. We can see that the overall reliability indices of
cases 2 to 5 in Table III are improved compared with that of
case 1 (base case). In cases 2, 3, 4 and 5, where the number of
DGs is limited at 1, 2, 3 and 4 unit respectively, see reductions
in the system ECOST.
TABLE III
RESULT OF CASE STUDY FOR RELIABILITY INDICES
Fig. 3 Flowchart for solution algorithm
Cases
Reliability indices
1 2 3 4 5
SAIFI
0.2482 0.2482 0.2482 0.2482 0.2482
(interruptions/customer yr)
SAIDI
3.7321 3.7290 3.7261 3.7253 3.7251
(hours/customer yr)
CAIDI
(hours/customer 15.036 15.024 15.012 15.009 15.008
interruption)
ASAI 0.9996 0.9996 0.9996 0.9996 0.9996
ASUI 0.0004 0.0004 0.0004 0.0004 0.0004
ENS (kWh/year) 40,775.30 40,509.70 40,256.00 39,455.40 39,443.60
AENS
21.37 21.23 21.10 20.68 20.67
(kWh/customer/year)
ECOST ($/year) 49,922.30 47,552.80 45,288.70 43,262.90 43,157.60
ECOST reduction (%) - 4.75 9.28 13.34 13.55

It is very interesting to note that the constraint given in (14)

is binding for these four cases. The reason is that to minimize
the system ECOST, as many DGs as possible should be
installed. Distribution system of RBTS bus 2 with distributed
generation for case 5 shown in Fig. 5 It is observed that a DG,
if its size is large enough, tends to be installed at the end of a
feeder.
Fig. 4 RBTS bus 2 Radial Distribution system

International Scholarly and Scientific Research & Innovation 7(10) 2013 1311 ISNI:0000000091950263
 World Academy of Science, Engineering and Technology
International Journal of Electrical and Computer Engineering
Vol:7, No:10, 2013

[7] F. Glover, Tabu search-part I. ORSA J. Computing, vol. 1, no. 3. 1989.

[8] M. Hiroyuki, and O. Yoshihiro, “Parallel tabu search for capacitor
placement in radial distribution system,” in Proc. Power Engineering
Society Winter Meeting, vol. 4, pp. 2334-2339, 2000 .
[9] L. Goel, and R. Billinton, “Procedure for evaluating interrupted energy
assessment rates in an overall electric power system,” IEEE Trans.
Power Systems, vol. 6, no. 4, pp.1398-1403, 1991.
[10] L. Goel, and R. Billinton,. “Basic data and evaluation of distribution system
reliability worth,” in Proc. Computer, Power and Communication Systems in
a Rural Environment, pp. 271-277, 1991.
[11] Y. J. Jeon, and J. C. Kim, Application of simulated annealing and tabu
search for loss minimization in distribution systems, Electrical Power and
Energy Systems, vol. 26, no.1, pp. 9-18, 2004.

Nattachote Rugthaicharoencheep (M’10) received his Ph.D. in Electrical

Engineering from King Mongkut’s University of Technology North Bangkok
(KMUTNB), Thailand in 2010. He is currently a lecturer at the Department of
Electrical Engineering, Faculty of Engineering Rajamangala University of
Technology Phra Nakhon (RMUTP), Bangkok, Thailand. His research
interests include power system operation, optimization technique, and
distributed generation.
Open Science Index, Electrical and Computer Engineering Vol:7, No:10, 2013 waset.org/Publication/17067

Thong Lantharthong received his M.Eng in Electrical Engineering from

Rajamangala University of Technology Thanyaburi, Pathumthani, Thailand in
2010. He is currently a assistant professor at the Department of Electrical
Engineering, Faculty of Engineering Rajamangala University of Technology
Fig. 5 Distribution system of RBTS bus 2 with distributed generation Phra Nakhon (RMUTP), Bangkok, Thailand. His research interests include
for case 5 power system operation, optimization technique, and distributed generation.

VII. CONCLUSION
The search for the best compromise among the objectives is
achieved by Tabu search technique for optimal placement and
sizing of distributed generation in distribution systems.
Employing DG in a distribution system results in several benefits
such as increased overall system efficiency. The effectiveness of
the proposed method was demonstrated by a case study of a
distribution network of RBTS bus 2. It can be seen from the
case study that distributed generation can reduce the customer
interruption cost and therefore improve the reliability of the
system. It is expected that our proposed method will be
utilized effectively for distribution system operator.

ACKNOWLEDGMENT
The authors would like to express his gratitude to
Rajamangala University of Technology Phra Nakhon,
Thailand for support.

REFERENCES
[1] D. T. Wang, L. F. Ochoa, and G. P. Harrison, “DG impact on investment
deferral: network planning and security of supply,” IEEE Trans. Power
System, vol. 25, no. 2, pp. 1134-1141, May 2010.
[2] J. Zhang, H. Cheng, and C. Wang, “Technical and economic impacts of
active management on distribution network,” Electrical Power and
Energy Systems, vol. 31, pp. 130-138, 2009.
[3] J. Mutale, “Benefits of active management of distribution networks with
distributed generation,” in Proc. Power System Conf. and Exposition.
pp. 601 – 606, 2006.
[4] L. F. Ochoa, A. P. Feltrin, and G. P. Harrison, “Evaluating distributed time-
varying generation through a multiobjective index,” IEEE Trans. Power
Delivery, vol. 23, no. 2, pp. 1132-1138, April 2008.
[5] R. Billinton, and R. N. Allan, “Reliability evaluation of power systems,”
pitman advanced publishing program, 1984.
[6] D. Bernaand, and A. Cigdem, “Simulation optimization using tabu
search,” in Proc. Winter Simulation Conference, pp. 805-810. 2000.

International Scholarly and Scientific Research & Innovation 7(10) 2013 1312 ISNI:0000000091950263