2012 IEEE Students Conference on Electrical, Electronics and Computer Science

Optimal Coordination Of Directional Overcurrent

Relays: A Genetic Algorithm Approach
Dharmendra Kumar Singh

Dr. S. Gupta

M.Tech.(Computer Technology),Student
Department Of Electrical Engineering
National Institute Of Technology
Raipur, India
e-mail: [email protected]

Associate Professor, MIEEE, MIET

Department Of Electrical Engineering
National Institute Of Technology
Raipur, India
e-mail: [email protected]

AbstractIn this paper, the optimization of coordination of

directional overcurrent relays in an interconnected power system
is presented. The objective of protective relay coordination is to
achieve selectivity without sacrificing sensitivity and quick fault
clearance time. The calculation of the time dial setting (TDS) and
the pickup current (Ip) setting of the relays is the core of the
coordination study. The inequality constraints guarantee the
coordination margin for each primary/backup relay pair having a
fault very close to the primary relay. Using this formulation, the
size of the optimization problem is significantly reduced. Genetic
algorithm is the algorithm being applied to minimize the operating
times of the relays. In this paper both linear and non linear
equations are framed for the test bus system used and optimized
using the genetic algorithm.

suitable settings such that their fundamental protective function

is met under the requirements of sensitivity, selectivity,
reliability and speed.

Justification for Using Genetic Algorithms

This paper describes a systematic overcurrent (OC)
protection grading method based on genetic algorithm (GA).
Genetic algorithms are computerized search and optimization
algorithms based on the mechanics of natural genetics and
natural selection. Professor Holland of University of Michigan
envisaged the concept of these algorithms in the mid 60s.
Thereafter a number of students and other researchers have
contributed to the development of this field. Genetic algorithm
Keywords-Directional overcurrent relay; genetic algorithm; is good at taking larger, potentially huge, search space and
navigating them looking for optimal combinations of things and
operating time of relays; pickup current; time dial setting.
solutions which we might not find in a life time.
I.
INTRODUCTION
C. Solution Approach
Abnormal conditions (faults, overload, overvoltage, etc.) can
For linear programming, the values of time dial setting
frequently occur in an interconnected power system. Due to
(TDS)
have been found for given values of pickup currents (Ip)
this, interruption of the supply and damage of equipments
[3]
and
for non-linear programming, the values of pickup
connected to the power system may occur. The faulted
currents
(Ip)
have been found for given values of TDS [4].
components must be readily identified and isolated in order to
guarantee the energy supply to the largest number of consumers Using GA subject to the constraints and hence the operating
possible and to maintain the system stability. Therefore a times of the relays is minimized.
reliable protective system is required. In-case the main
protective system fails; a backup protective system should exist
II.
OPTIMAL COORDINATION PROBLEM
to ensure reliability. This backup protection should act as a
Directional overcurrent (DOCR) coordination problem is a
backup either in the same station or in the neighboring lines parametric optimization problem, where different constraints
with time delay according to the selectivity requirement.
have to be considered in solving the objective function [1].
Directional over current relay (DOC) are commonly used for
Here the objective function to be minimized is the sum of the
power system protection. For coordination of these relays
operating times of the relays connected to the system, subject to
optimization methods are used. In these methods, at first, the
coordination constraints for each main and backup relay are the following constraints.
determined [1]. The other optimization methods that have been
A. Relay Characterstics
used for the above problems are simplex method [2].
A typical inverse time directional over current relay consists
of two elements, an instantaneous unit and a time overcurrent
A. Problems Associated With Above Techniques
A protective relay should trip for a fault in its zone and unit. The overcurrent unit has two values to be set: Pick up
should not, for a fault outside its zone, except to backup a failed current value (Ip) and the time dial setting (TDS). The pickup
relay or circuit breaker. The problem of coordinating protective current value (Ip) is the minimum current value for which the
relays in electric power systems consists of selecting their relay operates [5]. The time dial setting defines the operation

978-1-4673-1515-9/12/$31.00 2012 IEEE

B.

To ensure reliability of the protective system, the back-up

time (T) of the device for each current value, and is normally
given as a curve T vs M, where M is the ratio of relay fault scheme shouldnt come into action unless the primary (main)
fails to take the appropriate action. Only when CTI is exceeded,
current I, to the pickup current value i.e.
backup relay should come into action.
M = I/IP

(1) This case is expressed as:

In general, overcurrent relays respond to a characteristic

function of the type:

Tbackup Tprimary CTI

Where, Tbackup is the operating time of the backup relay

T = f(TDS, Ip, I)

(2)

After considering all these criteria, this problem can be

formulated mathematically as:

This function can be approximated as:

T=

K1* TDS
I
K2
CT Ratio* Ip
+K3

Tprimary is the operating time of the primary relay

(3)

min

n
i=1

Ti

Where K1, K2 and K3 are constants that depend upon the Where, n represents the number of relays.
specific device being simulated.
III. GENETIC ALGORITHM
The following two cases have been considered in this paper
Genetic algorithm (GA) is a search technique used
for obtaining the objective function and minimizing it.
in computing, to
find
exact
or approximate solutions
Case I
to optimization and search problems. Genetic algorithms
are categorized as global search heuristics [6]-[7]. Genetic
The linear equation is formulated in terms of TDS by taking
algorithms are a particular class of evolutionary algorithms (EA)
Ip to be constant at a particular value in (3)
that use techniques inspired by evolutionary biology such
Case II
as inheritance, mutation, selection, and crossover [8]-[9].
The non-linear equation is formulated in terms of Ip by taking
A. Terminologies Related To Genetic Algorithm
TDS to be constant at a particular value in (3).
x Fitness functionA fitness function is a particular type
B. Relay Settings
of objective function that prescribes the optimality of a
solution (that is, a chromosome) in a genetic algorithm
The calculation of the two settings, TDS and Ip, is the
so that particular chromosome may be ranked against
essence of the directional overcurrent relay coordination study.
all the other chromosomes.
In general, directional overcurrent relays allow for continuous
x
ChromosomeIn genetic algorithms, a chromosome is
time dial settings but discrete pickup current settings.
a set of parameters which define a proposed solution to
Therefore this constraint can be formulated as:
the problem that the genetic algorithm is trying to
solve. The chromosome is often represented as a simple
string, although a wide variety of other data structures
TDSi min TDSi TDSi max
are also used.
x SelectionDuring each successive generation, a
proportion of the existing population is selected to
IPimin IPi IPimax
breed a new generation. Individual solutions are
selected through a fitness-based process, where fitter
C. Coordination Problem
solutions are typically more likely to be selected.
x
ReproductionThe next step is to generate a second
In any power system, a primary protection has its own
generation population of solutions from those selected
backup one for guaranteeing a dependable power system. The
through crossover and mutation.
two protective systems (primary and back-up) should be
For each new solution to be produced, a pair of
coordinated together. Coordination time interval (CTI) is the
"parent" solutions is selected for breeding from the pool
criteria to be considered for coordination. Its a predefined
selected previously. By producing a "child" solution
coordination time interval and it depends on the type of relays.
using the above methods of crossover and mutation, a
For electromagnetic relays, CTI is of the order of 0.3 to 0.4 s,
new solution is created which typically shares many of
while for a microprocessor based relay, it is of the order of 0.1
the characteristics of its "parents". New parents are
to 0.2 s.
selected for each new child, and the process continues

SCEECS 2012

until a new population of solutions of appropriate size

is generated.
CrossoverIn genetic algorithms, crossover is genetic
operator used to vary the programming of a
chromosome or chromosomes from one generation to
the next. It is analogous to reproduction and biological
crossover, upon which genetic algorithms are based.
MutationIn genetic algorithms of computing,
mutation is a genetic operator used to maintain genetic
diversity from one generation of a population of
algorithm chromosomes to the next. It is analogous to
biological mutation.
The purpose of mutation in GAs is preserving and
introducing diversity. Mutation should allow the
algorithm to avoid local minima by preventing the
population of chromosomes from becoming too similar
to each other, thus slowing or even stopping evolution.

B. Outline Of The Genetic Algorithm

It is very clear from the flow chart given below in Fig. 1

Step 5: Minimization of objective function is carried out and

optimum values of TDS are determined using linear
programming technique in MATLAB.
Step 6: The values of TDS and minimization of objective
function is further optimized using genetic algorithm (GA)
using MATLAB toolbox.
IV.

SIMULATION AND RESULTS

The performance of the proposed method is evaluated by, a

6-bus overcurrent relay network shown in Fig. 2. The aim is to
find an optimal setting for network relays in order to minimize
the final operating time. The coordination time interval (CTI) is
assumed to be 0.2 seconds. The ratios of the current
transformers (CTs) are indicated in Table I. In order to solve
the relay coordination problem first, it needs to determine the
primary and backup relay pairs. Table II shows the
primary/backup relays pairs and faults currents.
Fig. 3. shows the convergence of the proposed method. The
results of the proposed method are shown in Table III and
Table IV for linear and nonlinear function respectively.

Create a Population of
Chromosomes
Determine The Fitness of
Each Evaluation
> n Gen.

yes

Display
Results

no
Select Next Generation
Perform Reproduction Using
crossover

Perform Mutation
Figure 2. Six bus test system
Figure 1. Flow chart
TABLE I.

V.

PROPOSED ALGORITHM

Relay No.

CT Ratio

CT RATIO
Relay No.

CT Ratio

Step 1: Generate P/B relay pairs. The knowledge of

primary/backup relay pairs is
essential in the formulation of
the coordination constraints.

1200/5

800/5

800/5

800/5

800/5

10

600/5

Step 2: Load flow analysis is done using Newton-Raphson

method to determine line currents. This analysis can be done
using any simulation softwares and here we have used ETAP
simulation software.

800/5

11

800/5

800/5

12

800/5

1200/5

13

1200/5

800/5

14

800/5

Step 3: Short circuit analysis is done using the same simulation

software to find fault currents. Step 4: Pickup current is
calculated based on the load current. Here it will be set at 1.5
times the maximum load current, for phase protection.

SCEECS 2012

TABLE II. PRIMARY/ BACKUP RELAY PAIRS AND FAULT CURRENTS

Backup
Relay

Fault Current
(KA)

Primary
Relay

Fault Current
(KA)

TABLE IV.

RESULT USING GA (FOR NONLINEAR FUNCTION)

Pickup
Current

Obtained
Values (A)

Pickup
Current

Obtained
Values (A)

4.5890

6.0720

Ip1

0.6700

Ip8

0.5000

4.5890

14

5.4570

Ip2

0.5000

Ip9

0.5000

0.8680

2.3510

Ip3

0.5000

Ip10

0.5000

0.8680

14

5.4570

Ip4

0.5000

Ip11

0.5000

1.3650

4.8030

Ip5

0.5000

Ip12

0.8660

0.5528

30.5470

Ip6

0.5000

Ip13

0.5000

0.9770

13

17.8210

Ip7

0.5350

Ip14

0.5680

0.8610

12

17.7050

1.4830

2.3510

1.4830

6.0720

0.7670

18.3380

0.6390

10

4.0770

10

0.9455

11

30.9390

11

1.0740

2.8380

11

1.0740

4.4960

12

3.4220

5.1860

12

3.4220

4.4960

13

0.6010

18.1720

14

1.7640

5.1860

14

1.7640

2.8380

Fitness Function Value= 10.67856 sec

VI.

CONCLUSION

In this paper, an optimization methodology is presented to

solve the problem of coordinating directional overcurrent relays
in an interconnected power system. The operating time of the
relays was determined using genetic algorithm for 6 bus system
for both linear and non linear objective functions. For linear
function value of time dial setting (TDS) and for nonlinear
function value of pick up current (Ip) are found. This method
increases coordination and the operation speed of relays.
Finding the absolute optimal point, the ability to be applied on
large networks, the ability to consider both linear and nonlinear characteristics of relays are some of advantages of the
proposed method.
REFERENCES

Figure 3. Convergence of the proposed method

TABLE III.
Time Dial
Setting

RESULT USING GA (FOR LINEAR FUNCTION)

Obtained
Values (A)

Time Dial
Setting

Obtained
Values (A)

TDS1

0.55647

TDS8

0.50000

TDS2

0.50010

TDS9

0.50003

TDS3

0.50000

TDS10

0.50000

TDS4

0.50004

TDS11

0.50000

TDS5

0.50002

TDS12

0.59127

TDS6

0.50000

TDS13

0.50000

TDS7

0.51584

TDS14

0.52268

[1] L. Perez, A. J. Urdaneta, "Optimal Computation of

Directional
Overcurrent Relays Considering Definite Time Backup Relaying ", IEEE
Transaction on Power Delivery, Vol. 4, No.4, October 1999, PP. 385-388.
[2] L. Perez, A. J. Urdaneta, "Optimal Computation of Distance Relays
Second Zone Timing in a Mixed Protection Scheme with Directional
Relays", IEEE Transaction on Power Delivery, Vol. 16, No.3, July 2001,
PP. 1276-1284.
[3] B. Chattopadhyay, M. S. Sachdev and T.S. Sidhu, "An On-Line Relay
Coordination Algorithm for Adaptive Protection Using Linear
Programming Technique," IEEE Trans. Power Del, Vol 11, No. 1, PP
165-173, Jan. 1996.
[4] AI. Urdaneta, R. Nadira and L. G. Perez Jimenez, "Optimal Coordination
of Directional Overcurrent Relays in Interconnected Power Systems,"
IEEE Trans. Power Del, Vol 3, No. 3, PP 903-91 I, July.1988.
[5] Farzad Razavi, Hossein Askarian Abyaneh - A new comprehensive
Genetic Algorithm method for Optimal overcurrent relays coordination
International Journal on Electric Power Systems Research, 78(2008) Pg:
713-720.
[6] David E. Goldberg - Genetic Algorithms in Search, optimization and
Machine Learning - Addison-Wesley, Reading MA, 1989.
[7] David E. Goldberg - Genetic Algorithm Approach: Why, How, and What
Next? -Proceeding of conference on Adaptive and Learning System:
Theory and Applications, May 1985, New Heaven, CT, USA.
[8] Cheng-Hung Lee and Chao-Rong Chen - Using Genetic Algorithm for
overcurrent Relay Coordination in Industrial Power System - The 14th
National Conference on Intelligent Systems Applications to Power
Systems, ISAP 2007 National Taipei University of Technology, Taiwan.
[9] C. W. So, K. K. Li, K. T. Lai and K. Y. Fung, "Application of Genetic
Algorithm to Overcurrent Relay Grading Coordination," APSCO. 97, Vol
1, PP. 283-287 ,Nov. 1997.

Fitness Function Value=10.729477sec

SCEECS 2012