See discussions, stats, and author profiles for this publication at: https://www.researchgate.

net/publication/328031789

Protection Coordination with Distributed Generation in Electrical System of

Iraqi Distribution Grid

Article · October 2016

CITATIONS READS

0 586

1 author:

Rashid H AlRubayi
University of Technology, Iraq
35 PUBLICATIONS   28 CITATIONS   

SEE PROFILE

Some of the authors of this publication are also working on these related projects:

Electrical Design of Multi Store Building View project

STABILITY ENHANCEMENT BASED ON DISTRIBUTED GENERATION UNITS IN AL-NAJAF AL-ASHRAF ELECTRICAL NETWORK View project

All content following this page was uploaded by Rashid H AlRubayi on 05 November 2018.

The user has requested enhancement of the downloaded file.

Eng. &Tech.Journal, Vol.34,Part (A), No.6,2016

 
Protection Coordination with Distributed Generation in Electrical
System of Iraqi Distribution Grid

Dr. Rashid H. Al-Rubayi

Electrical Engineering Department, University of Technology\Baghdad
Email:[email protected]
Ammar Abbas Majeed
Ministry of Electricity\Baghdad
Email:[email protected]

Received on:19/11/2015 & Accepted on:9/3/2016

ABSTRACT
The nature characteristic in conventional distribution networks is radial by single source
supplying a downstream network. The interest about the environmental impacts and
development in technologies have led to increase distributed generation (DG) interconnected in
distribution networks. Protective device coordination will be affected by adding DGs to the
existing network through participating to the change in direction of power flow and fault current
values and direction which cause loss in settings and mis-coordination for protective devices,
especially over current relays.
The effect of DG on coordination depend upon number, location and size of DG, so in
this work, the Particle Swarm Optimization (PSO) technicality utilization to locate optimal
location and sizeof DG to obtain minimum active power losses.The Time Current Characteristic
(TCC) curves represented which depended on the over current protection relays parameters to
find settings and limited any loss in it, in order to reset these relays to obtained the proper
operation without intersections in time of operation and satisfy optimal coordination between
primary and pack up over current protection relays.
In this work two soft wares are used, the first is Matlab R2014a for implementation of the PSO
algorithm while the second software is CYMDist program for load flow analysis, short circuit
current calculation and protection coordination device analysis. To verify the developed
algorithm parts from Iraqi distribution network (Baghdad Al-Rusafa 33KV distribution
networks). So, used two DG units with total capacity 50MW distributed in 33kV of South Al-
Rusafa distribution network which represented about 9.4% from the total load of this network
533.5MW, the total active power losses reduced from 11.597MW to 6.658MW with losses
reduction 6.96MW about 43% from total losses.  
Keywords: Distributed Generation (DG), Particle Swarm Optimization (PSO), Short Circuit
Level, Protection Coordination, and Time Current Characteristic (TCC).

INTRODUCTION

O
n the power system, the increasing demand has raised a difficult assignment to electrical
power engineers in preserving a reliable power system. With high loads network, when
the system losses and voltage drop increase in distribution system, the environmental
variation in design and operation of power systems has required the necessity to activate power
distribution network by penetration of an emerging technology in this new period and it
provides clean electric power by distributed generation (DG), which is also known in many
terms, “Dispersed generation”, “Embedded generation” and “Decentralized generation” [1].
Distributed generation refers any electricity generating technology units located in the effective
point of the electric distribution system near to the load center, in order to avoid the extend in
the distribution network to cover the increment load and future development for consumers
demand or reduce the necessary of energy transfers [2]. The optimal selecting size and location
1161
 
Eng. &Tech.Journal, Vol.34,Part (A), No.62016 Protection Coordination with Distributed Generation
in Electrical System of Iraqi Distribution Grid

of DGs in electrical distribution system is a complicated problem for optimization and if this
problem contain the multiple objective function, this problem become more complex, where
optimization is a mathematical tool which can be used to locate and size the DG units in the
system, so as to utilize these units optimally within certain limits and constraints [3].
One of the standard protections of distribution substation which detects fault circumstances on
feeders and commences a trip signal to separate the faulted equipment is the feeder protection.
The radial configuration has been completely used to facilitate the schemes of protective
relaying in the traditional feeder protection in distribution systems, especially in coordination
between downstream and upstream protective devices of resulting simple relaying. With the
DGs existence new sources of fault currents comes from DGs, and the advantages of protection
design obtained from radial currents flows are lost, sothe protection coordination between
upstream and downstream faced a big challenge[4].
Installation a DG into the distribution system makes a variation in the level of a system fault
current and cause many of problems in the protection of distribution system, like protection
blinding, short-circuit levels increase, false tripping in the protective devices and mal-
coordination between relays.To ensure selectivity, it is necessary to set a proper coordination
between relays and other protective equipment. Nevertheless, this coordination may be severely
obstructed if a DG is connected to a distribution network [5].

The Difficulties and Challenges with DG

i. Power loss: Depending on the DG location, the losses may increase or decrease, the
DG capacity and the proportional size of the load amount, the network topology and
other factors.DG unit's location is an important criterion that needs to be analyzed to
achieve best system reliability with losses reduction [6].
ii. Power quality :The quality of power are generally considered of the fluctuations of
transient voltage and the mains voltage harmonic distortion.The impact of increasing
fault level of the network during adding distributed generationusually leads to the
improvement of the power quality, [7].Important power quality problem is a flicker. It
refers to rapid voltage changes that can cause significant variations in lighting and
electronics to interrupt operations. To reducing flicker of the voltage during decrease
the power conditioning system output impedance in standalone mode [6].
iii. Reliability: The DG units can improve system reliability by ensuring the continuity of
power supply. In case of power failure occurring on the upstream main network, the
disconnecting a part of the system (containing DG) from the main network can be
performed by means of the islanding procedure. In this way, the DG units continuously
provide the power to customers of the islanding system. [8].
iv. Voltage Regulation: DG connection can lead to changes in the feeder's voltage profile
by changing the active and reactive power flow magnitude and direction. Also for
short-term abnormal voltage can occur when DG is switched off and on. The obtained
fluctuation result can have the disturbance impact on the distribution system. However,
voltage regulation can be effect positively or negatively when DG penetrated to the
distribution system depending on its characteristics and location [9].
v. Short Circuit Level: The short circuit levels of the distribution system affected by the
existence DG in this network. It establish high fault currents with respect to the normal
case when the network operate without DG.. The levels of short circuit can be changed
enough in the case of a few large units, or many small units, to cause mal coordination
between the protection devices such as relays or fuses [10].When a fault occurs, at fault
point the fault current is fed from both theutility power system and distributed
generation. The fault cannot be separated, and so continue If the total fault current
exceeds the capacity of the power supply circuit breaker [11].

1162
 
Eng. &Tech.Journal, Vol.34,Part (A), No.62016 Protection Coordination with Distributed Generation
in Electrical System of Iraqi Distribution Grid

vi. Protection coordination: The purpose of protection is to isolate the faulty equipment
as quickly as possible, with a minimum of disruption to the supply of electricity, and
minimal impact system stability. The protection can be considered as a control system
as well. It monitors the operation components of the power system such as lines,
transformers and generators, and takes action if the component is in a potentially
dangerous state, This action is normally opening a circuit breaker, which eliminates
completely the supply from the component.[12].

Benefits of DG
The considerable advantages of the installation of DG into electric power networks are as
follows [13]:
a) Loss reduction in transmission and distribution lines, which minimizes costs related to loss.
b) Improved reliabilityof the utilitysystem as back-up or stand by power to customers.
c) Power quality improvementand voltage supportoffers customers a choice in meeting their
power needs.
d) Release Transmission and distribution lines capacity and reduces or Postponement the
needs of new infrastructure.
e) Loweringcoststhrough avoidinglong distanceof the high voltagetransmission lines.
f) Where renewable sources are used with some technologies such as solar or wind DG will be
friend of environment.
g) Can reduce the needfor the construction orupgradesoflargeinfrastructurebecause theDGcan
bebuilt on the siteof the load.
h) Spinning reserve support during generation outages.
One of the objectives for this paper is to obtain the optimal number, size and locations of DG
units to avoid the negative impacts of the ordinary installation of the DG .and then calculate the
settings for the relays protective device with and without DG to clarification the coordination
with new settings.

Mathematical Formulas
i. Power Flow Calculation: There are three ordinarily used iterative methods, precisely
Gauss-Seidel, Fast Decouple andNewton-Raphsonmethods of analysisto
thenonlinearalgebraic equation. These methods are usedto the powerflow problems
solution [14]. For considerable power system, the most active and practical method is
the Newton-Raphson found. [11].
ii. DG optimization using PSO Technique :TheDGplanning issuesconsists of twostages,to
find the optimallocationbusin the networkand theoptimalDGsize[15].The PSO-based
approach to find the optimal size and location of distributed generation unit to reduce
active power loss as in the below, Figure (1) represents flow chart of PSO.

1163
 
 Eng. &Tech.Journal, Vol.34,Part (A), No.62016 Protection Coordination with Distributed Generation
in Electrical System of Iraqi Distribution Grid

START
 
Enter system data  
 
Create initial population of particles (random size and location of DG)  
 
 
  Calculate real power loss of each particle
 
 
  Put the current
 
  Set the current overall best  
 
 
  Updating the velocity and position of particle

If k = max. Put new iteration

iteration k= k+1
No
Yes

  Print out the optimal solution

End

Figure (1): Flow chart of PSO

iii. Over Current Protection and Coordination [16],[17],[18]:Overcurrent relays are usually
supplied with time delay element and an instantaneous element in the same unit. An
earth-fault unitand three-phase overcurrent unit are the most modernistic microprocessor
protection has within the same case. Overcurrent protection relays setting include
choosing the time\current characteristic parameters of both instantaneous and the time
delay units. The operation achieved in two cases, first one to the phase overcurrent
protection relays and the earth fault protection relays secondly.
a) Instantaneous units: In distribution lines to set these units, used one of the two values
either between six and ten times the maximum circuit rating or 50 per cent of the maximum
short circuit current at the point connection of the CT supplying the relay.
.
(1)
b) IDMT unit: Inverse definite minimum time protection relays can be modify through
choosing two parameters, the time multiplier setting or time dial, and plug setting (tap
setting) or the pick-up.
 pick-up setting:Expressed to multiples of the pick-up relay and fault currents viewed
via use the relay. This value is ordinarily represented the Plug Setting Multiplier (PSM),
it is recognize as the plug setting or the ratio of fault current in secondary amps to pick-
up of the relay. A pick-up setting for phase relays is specified toward permitting a
margin higher than the nominal current in overload cases, as in the expressions below:

(2)

1164
 
Eng. &Tech.Journal, Vol.34,Part (A), No.62016 Protection Coordination with Distributed Generation
in Electrical System of Iraqi Distribution Grid

(3)
Where: PSM = plug setting multiplier, Inom = nominal circuit current rating, CTR = current transformer ratio
OLF = the overload factor which depends on the parameter being protected .It's rang between (50% and
200% in 25% steps). For generators, lines and transformers it recommendedordinarily in the range of 1.25 to
1.5.While in case of a possibility under emergency conditions to increment the load on feeders such as in
distribution systems, OLF it could be up to the order of 2.
In earth-fault overcurrent relays, under normal operating conditions the setting of pick-up is
specified taking into consideration the maximum disturbance that would present in the system.
Atypical disturbance permitted between (10% and 80% in 10% steps).
 Time dial setting:When the fault current arrives to a magnitude similar to or more than
the setting of relay current, before relay operates the time delay adjusted by time dial
setting .Also noted by the time multiplier setting (TMS) which its range between (0.05-
1 in 0.05 steps). The procedures for calculating (TMS), to find the convenient
coordination and protection to the system are explain in steps below:
Step1: Find the wanted operating time (t1) of the remote relay from the supply source through
utilizes the minimum TMS and dependence on a fault level. To coordinate with devices
installed downstream which it is necessary; this TMS may have to be higher.
Step2: Find the relay operating time related with the next substation breaker towards the source
t2a t1 td (4)
Where: t2a the back-up relay desired operating time, t1 is the operating time of the
downstream relay,
td is the discrimination (margin) time for breaker - breaker ( IDMT to IDMT relays).

td 0.25t1 0.25 (5)

Where: Relay Error Factor = 0.25t1.
Downstream breaker interrupting Time = 100 msec.
Relay Overshoot Time = 50 msec.
Safety Margin = 100 msec.
Step3: The time dial setting calculation for back-up relay through the same current of fault as in
the previous 1 and 2 , and find t2a so as the pickup value.

TMS (6)
.
Step4: Find the actual operating time (t2b) of back-up relay, but now using the fault level and
pickup value where this relay is connected.
t2b Operating time at selected PSM and TMS 1.0 TMS (7)
Step5: Go on with the sequence, beginning from the step2.
Procedure explained above for earth and phase units can readily utilize when the operating
characteristics of relays are known by mathematical formulae in lieu of log-log paper curves.
Define the operating time mathematically by ANSI/IEEE and IEC Standards in the following
expression:
(8) Where: t = operating time of relay in seconds, Is = selected pick-up

current, I = secondary fault current level

in amps, k = time multiplier setting, or time dial, L = constant.

The parameters and determine the relay characteristics slope. These constants with L to
different standard of the overcurrent relays processed under IEC and ANSI-IEEE Standard with
typical characteristics for both types are given in Appendix (B).

1165
 
Eng. &Tech.Journal, Vol.34,Part (A), No.62016 Protection Coordination with Distributed Generation
in Electrical System of Iraqi Distribution Grid

Table (1) Electrical distribution network size for part of Baghdad Al-Rusafa 33KV Network
Items substations Substation Transformers North Net. Center Net. South Net.
Total No. Total No. S.S. T.R S.S. T.R S.S. T.R
1 132/33/11KV 6 17 -- -- -- -- 6 17
2 132/33KV 15 22 7 11 7 8 1 3
3 132/33KV mobile 1 1 -- -- -- -- 1 1
4 33/11kv 57 118 18 35 16 34 23 49
5 33/11kv mobile 2 2 -- -- -- -- 2 2
Simulation Results and Discussion for case study
The work is dedicated to the simulation results obtained from the different analysis and
calculations for part from Baghdad electric distribution network. Mainly the works are divided
in to stages as below:
a. Solving the optimization problem of DG by used Particle Swarm Optimization
(PSO) methodswhich has been discussed previously through use (MATLAB
R2014a). So the optimal number, sizes and locations of DG will be accomplish in
this stage.
b. Find short circuit calculations and protection coordination device analysis by
CYMEDist program package with and without DG.
i. Baghdad Al-Rusafa Distribution Network: Part of Baghdad Al-Rusafa distribution
network, the region restricted between army channel and Tigris river and the
extended from north to the south of the city with its outskirt.Table (1) illustrate the
size of choice part from of Baghdad Al-Rusafa distribution network.
ii. South of Al-Rusafa33KV Distribution Network: The total number of feeder in this
part of network is (51) feeders with two types transmission lines, underground and
overhead with substations and transformers illustratein the bus and line data of the
distribution network are given in Appendix(A)and diagrammatically taken from
Iraqi MOE. This part show the highest total loads for south of Rusafa distribution
network (533.5MW) compares with two others parts (north and center).The PSO
algorithm which was used for the optimal size and locations of DG units near the
loads in high voltage side at the 33/11KV transformers with objective function of
minimum active power losses, the results to these optimizations in each case as
shown in Table (2).The power factor is kept constant (cos ∅ 0.85 .

Table (2) Optimal Number, Locations and Sizes of DGwith Total Power Losses for South of Baghdad –Al-
Rusafa 33KV Network
Case Numbers DGs Size MW Locations (Bus No.) Losses MW Losses Reduction %
1 Without DGs ----- ------ ١١٫٥٩٧ ------
2 1 30 61 ٨٫٣١١ 28
3 1 50 59 9.35 19
4 2 30 62 6.927 40.3
20 61
5 2 25 61 6.658 43
25 62
6 3 15 62 6.910 40.4
15 60
15 61
7 3 20 68 6.947 40.1
15 61
15 62

1166
 
Eng. &
&Tech.Journall, Vol.34,Part (A),
(A No.62016 Protection
n Coordinatio on with Distriibuted Generration
in Electtrical System of Iraqi Distrribution Grid
d

Fromm the results the optimall locations aand size app pears in casee (5) when uuse two DG G with
25MW W for each one
o of them ini location 661(Dawaa sub bstation) and
d 62 (Zaafranniya substatiion) to
obtainned the miniimum active power lossees to this parrt of network k. The first loocation 61 lo
ocated
at thee end of longg 33KV feedder which hass three distriibution transfformers 33/111 KV at thee same
feedeer with high total loads and
a high dropp voltage. Also
A the second location 662 located at a high
loadss distributionn transformer 33/11KV V with high h drop voltaage too.Figuure (2) com mpares
betweeen these casses and show
ws the reductiion of powerr losses for all
a cases.

TOTTAL POWER
R LOSSES
LOSSES(MW)
POWERL 

15
10
5 1
11.597 8.311
8 9.35 6.927
6 6.658
6 6.91 6.947
6
0
witthoutDG with 1 DG w ith 1DG with 2DG with 2DG wiith 3DG with 3DG
(3
30MW) (550MW) (30
0+20)MW (2*
*25)MW (3**15) MW (20+
+2*15)MW

Figu
ure (2) : Tottal power lossses for South of Al-Russafa 33 Netw
work

Two DG units wiith total sizess 50 MW (255 MW to eacch unit) whicch represent 9.4% with reespect
to thee total load and 5MW losses reducttion about 43% from thee total activee power losses of
southh distributionn network.
iiii. Shorrt Circuit Current Calcullation forEm
mbedded Two o Distributedd
Geneeration(25M MW) Near Al--Dawaa and d Al-Zaafraniya 33/11KV V Substationss

In thiis case, the effect in the short

s circuit ccurrent was determined
d through
t embbedded two DGsD in
size oof 25MW foor each one which repressents the optimal numbeer, size and location thatt PSO
limiteed, one of thhem near Al-D Dawaa substtation and th he other one near
n Al-Zaaffrainya substtation.
Figurre (3) showss the 33kv feeders
fe and ssubstations and
a the resullts of short ccircuit curreent for
differrent type of faults
f after DGs
D embeddeed. All resullts are summ marized in Tabble (3).

Figurre(3) 33KV Feeders (Qu

usaiba_Zaaffraniya and
d Qusaiba_D
Dawaa)withT
Two DGs

1167
 
 Eng. &Tech.Journal, Vol.34,Part (A), No.62016 Protection Coordination with Distributed Generation
in Electrical System of Iraqi Distribution Grid

Table(3) Short Circuit Currents Magnitude for Embedded two DGs

Fault QusaibaS.S.Amp. Dawaa S.S. Amp. Zaafra S.S. Amp. 25MW DG 25MWDG
Without With Without With Without With Dawaa Amp. Zaafra.
Type
DG DG DG DG DG DG Amp.
L LL 8241 10725 2431 4054 4214 6055 4029 5954
LG 215 604 191 533 202 566 532 565
LL 7137 9315 2105 3534 3650 5268 3513 5181
LLG 7137 9318 2124 3568 3661 5293 3546 5205

iv. Setting and Protection Coordination

Through using protective device coordination analysis in CYMEDist program to represent
Time Current Characteristic (TCC) curves by branch device coordination feature, for south
of Al-Rusafa distribution network with this simulation steps for relay coordination:
a. The system is simulated without DG to get the fault Currents.
b. The relays are set for the fault currents and coordinated without DG.
c. DG is connected to the system to get the new fault currents.
d. The operating time is checked after DG connection.
e. The relays are coordinated with new fault currents and operating time.
f. The new coordination margin is calculated.
g. Loss of coordination is checked.
h. New settings are selected for the relays to accommodate the DG

Case (1): Without DG

The current network depend on one circuit breaker in 33KV side of Qusaiba substation contain
over current relays with specifications in Appendix( C ) which are used practically in
distribution network.
a. Qusaiba-Al-Dawaa 33kv feeder which protected by circuit breaker in Qusaiba 33kv bus
bar with its over current relay. Tables(4), (5), shows the parameters for (I.D.M.T.) relay
and operating time for 3- phase fault and earth fault (E/F) respectively. To represent the
points of TCC curves assuming load current increase 110% from pickup current at
beginning point and the fault current in Qusaiba, Sayloo and Dawaa as the other points
in curves.
Table(4) Parameters and Operating Timefor Over Current Relay (3-phase fault & instantaneous) in
C.B. of Qusaiba S.S. Without DG
TMS = 0.3 PS = 6.33 Ipu = 760A C.T.ratio = 600/5
Fault locations L-L-L Fault current(A) Operating Time(second)
110% Ipu 836 22.012
Dawaa B.B 2431 1.787
Sayloo B.B 2964 1.523
Qusaiba B.B 8241 0.860
Instantaneous : Ipu = 4200 Delay Time = 0.06 sec.

Table(5) Parameters and Operating Timefor Over Current Relay (E/F & instantaneous) in C.B. of
Qusaiba S.S. Without DG
TMS = 0.1 PS = 0.2 Ipu = 24A C.T.ratio = 600/5
Fault locations L-G Fault current(A) Operating Time(second)
110% Ipu 26.4 7.337
Dawaa B.B 191 0.331
Sayloo B.B 197 0.326
Qusaiba B.B 215 0.312
Instantaneous : Ipu = 120A Delay Time = 0.06 sec.

1168
 
Eng. &
&Tech.Journall, Vol.34,Part (A),
(A No.62016 Protection
n Coordinatio on with Distriibuted Generration
in Electtrical System of Iraqi Distrribution Grid
d

Figurre (4) repressent(TCC) cuurves in Logg-Log paper illustrate op perating timee for over current
relays in differennt values andd locations (points in Tables(4),
T (5)) above), thrree phase faault in
phasee protection and earth faaults(E/F) ffor ground protection,
p allso shows thhe informatio on for
relays used.
Asweell Figure(5) shows the instantaneou
i us feature thaat makes the relays operrating time will
w be
set ass(600 ms) whhen fault excceed (4200AA) for three phase
p fault an
nd(120A) forr earth fault,, these
valuees represent 50%
5 maximu um faults cuurrently approoximately, all values willl be instantaaneous
operaating.

Figu
ure (4) :Time Current Cha
aracteristic (T
TCC) Curvess for Qusaiba
a S.S. O/C (3--phase fault&
&E/F)
Rellays Without DG

1‐Relay: 3 ph /F   with Inst.Q
Qusaiba  
2‐Relay: E/F with Inst. Qusa aiba 
Figure (55) :Time Currrent Charactteristic (TCC) Curves Qussaiba S.S. O/C
C(3-phase
fault&E
E/F)Relays W
With instantaaneous (Witho
out DG)

b
b. Qusaiba-A Al-Zaafraniyya 33kv feedder which protected throu ugh circuit bbreaker in Qu
usaiba
33kv buss bar with itss over currennt relay. Tabbles(6), (7), shows the pparameters foor this
relays andd operating time,
t Assumming load currrent increasee 110% from m pickup currrent at
beginningg point and the
t fault currrent in Qusaiba, Zaafraniya and ten ttimes from pickup
p
current (110*Ipu) as th
he other poinnts in curves.

1169
 
Eng. &
&Tech.Journall, Vol.34,Part (A),
(A No.62016 Protection
n Coordinatio on with Distriibuted Generration
in Electtrical System of Iraqi Distrribution Grid
d

Table (6) Parameters and Operatin ng Timefor Over

O Current Relay (3-pphase fault &
insta
antaneous) in C.B. of Qusaaiba S.S. Withhout DG
T
TMS = 0.3 PS =5.5 Ipu = 6660A C.T.ratio
C = 60
00/5
Fault locations L-L-L FFault current((A) Operatinng Time(secon
nd)
110% % Ipu 726 22.012
Zaafraniyya B.B 4214 1.112
10* Ipu B.B 6600 0.891
Qusaiba B.B 8241 0.811
Insstantaneous : Ipu
I = 4200 Deelay Time = 0.06 sec.

Table(7) Parrameters and Operating T Timefor Overr Current Rellay (E/F & innstantaneous)) in
C.B. of Q
Qusaiba S.S. Without
W DG
T
TMS = 0.1 PS = 0.2 Ipu = 24A C.T.ratio = 600
0/5
Fault locations L-G Fau ult current(A
A) Operatinng Time(secon
nd)
110% % Ipu 26.4 7.337
Zaafraniyya B.B 202 0.322
10*Ipu B.B 209 0.317
Qusaiba B.B 215 0.312
Insstantaneous : Ipu
I = 120A Deelay Time = 0.06
0 sec.

Figurres (6), (7) presented

p (TC CC) curves ffor over curreent relay withh different m
magnitudes of
nce also illusttrated.
three phase faultss and earth faaults, operatinng time, relaay performan

1‐Relayy: 3 ph /F   Qusa
aiba 
2‐Relay: E/FQusaiba 

Figgure (6):Time Current Characteristic (TCC

C) Curves for Qusaiba S.S.O
O/C (3-phase faault &E/F) Rellays
Without DG G 

1170
 
Eng. &
&Tech.Journall, Vol.34,Part (A),
(A No.62016 Protection
n Coordinatio on with Distriibuted Generration
in Electtrical System of Iraqi Distrribution Grid
d

1‐Relay: 3 ph /Fw
with Inst.Qusaiba
a
2‐Relay: E/F with
h Inst.Qusaiba 

Figgure (7) :Time Current

C Chara
acteristic (TCCC) Curves Qussaiba S.S. O/C(3-phase fault&
&E/F)Relays With
W
instanttaneous (Witho
out DG)

Case (2) Embedd ded Two Disstributed Geeneration (2 25MW) Near Al-Dawaaa and Al-
Zaafrfraniya 33/111 KV Substtations:
In caase of penetraate the optim mal number, location and d size DGs to
t the 33 KVV network. Tables
T
(8), ((9), (10), (111), (12), (13)) illustrate thhe setting paarameters and
d the operatiing time resu
ults to
the oover current protection reelays (3-phaase fault and nstantaneouss) in Qusaiba, Al-
d E/F with in
Dawaaa and DG station
s circuiit breakers reespectively when
w two DG with (25M MW) for eacch one
embeedded near thhe load buses of Al-Daw waa and Al- Zaafraniya
Z 33/11
3 KV sub
ubstations wh hich is
represent the optimal location ns to south off Al-Rusafa 33 3 KV netwo ork.

Tablle(8) Parameteers and Operatting Timefor O

Over Current Relay
R (3-phase fault& instanttaneous) in C.B.
C of
Qusaibba S.S. With Two
T DG
TMMS = 0.55 PS = 4 Ipu = 480A C.T.ratio
o = 600/5
Fault locaations L-L-L Faault current(A)) Operatinng Time(second
d)
1100% Ipu 528 40.356
Dawaa B.B. 4054 1.767
Qusaiba B.B. 10725 1.247
Instaantaneous : Ipu = 5400A Delay Timme = 0.05 sec.

Tablle(9) Parameters and Operatting Timefor O

Over Current Relay
R (E/F & in
nstantaneous) in C.B. of Qu
usaiba
S..S. With Two DG
D
TMMS = 0.2 P = 0.4
PS Ipu = 48A C.T.ratio = 600/5
Fault locaations L-G Faultt current(A) Operatinng Time(second
d)
110% Ipu
I ٥٢٫٨ ١٤٫٦٧٥
Dawaa B.B. ٥٣٣ 0.568
Qusaiba B.B. ٦٠٤ 0.539
Instaantaneous : Ipu = 300A Delay Tim
me = 0.05 sec.

Tablle(10) Parametters and Opera

ating Time forr Over Currentt Relay (3-phase fault & instaantaneous) in C.B. of
Daw waa S.S With Two
T DG
T
TMS = 0.3 PS = 2.5 Ipu =300A
= C.T.ratio
C = 6000/5
Fault locations L-L-L Faault current(AA) Operatingg Time (secon
nd)
110% Ipu ٣٣٠ 222.012
Dawaa B.B B. 4054 0.786
Qusaiba B.B B. 10725 0.680
Instaantaneous : Ipuu = 2040A Delay Time = 0.005 sec.

1171
 
Eng. &Tech.Journal, Vol.34,Part (A), No.62016 Protection Coordination with Distributed Generation
in Electrical System of Iraqi Distribution Grid

Table(11) Parameters and Operating Timefor Over Current Relay(E/F & instantaneous) in C.B. of Dawaa S.S
With Two DG
TMS = 0.1 PS = 0.3 Ipu = 36A C.T.ratio = 600/5
Fault locations L-G Fault current(A) Operating Time(second)
110% Ipu 39.6 7.337
Dawaa B.B. ٥٣٣ 0.253
Qusaiba B.B. ٦٠٤ 0.241
Instantaneous : Ipu = 240A Delay Time = 0.05 sec.

Table(12) Parameters and Operating Timefor Over Current Relay (3-phase fault& instantaneous) in C.B. of DG
Station Near Dawaa S.S.
TMS = 0.4 PS = 4 Ipu =480A C.T.ratio = 600/5
Fault locations L-L-L Fault current(A) Operating Time(second)
110% Ipu 528 29.350
DG station ٤٠٢٩ ١٫٢٨٩
Dawaa B.B. 4054 ١٫٢٨٥
Instantaneous : Ipu = 2400A Delay Time = 0.06 sec.

NOTE:
 Appendix(D) shows Figures (10), (11), (12), (13)
 All protected schemes data are practically used in Rusafa distribution network
according to Iraqi MOE.

Table(13) Parameters and Operating Time for Over Current Relay (E/F & instantaneous) in C.B. of DG
Station Near Dawaa S.S.
TMS = 0.2 PS = 0.3 Ipu = 36A C.T.ratio = 600/5
Fault locations L-G Fault current(A) Operating Time (second)
110% Ipu 39.6 14.675
DG station ٥٣٢ 0.506
Dawaa B.B. ٥٣٣ 0.506
Instantaneous : Ipu = 264A Delay Time = 0.06 sec.

Figures (8),and (9) represent the (TCC) curves which is shows the coordination in operating
timebetween the over current relays against the faults. This coordination satisfies DG continued
to supply the load when the utility network is separated during faults. The short circuit current
magnitude with two DGs increases more than its magnitude in cases without DG, so it makes
(TMS) will be increase. The pickup current will be reducing from utility source because the
DGs supplied part of loads. These two changes effect on the setting and coordination. Tables
(14), (15), (16), (17), (18), (19) shows the setting parameters and the operating timeresults to the
over current protection relays (3-phase fault and E/F with instantaneous) in Qusaiba, Al-
Zaafraniya and DG station circuit breakers respectively when two DG (25MW) embedded near
the load buses of Al-Dawaa and Al- Zaafraniya 33/11 KV substations which is represent the
optimal locations to south of Al-Rusafa 33 KV network.

1172
 
Eng. &
&Tech.Journall, Vol.34,Part (A),
(A No.62016 Protection
n Coordinatio on with Distriibuted Generration
in Electtrical System of Iraqi Distrribution Grid
d

1‐Rellay : E/F DG statio
on 3‐Relay: E/F A Al‐Dawaa5‐ Relay  : E/F Qusaiba
2‐ Reelay: 3 ph /FDG stta on                 4‐‐Relay: 3 ph /F  A
Al‐Dawaa 6‐ Relay : 3 ph /F  Qu
R usaiba
F
Figure (8) :Tim
me Current Cha
aracteristic (TC
CC) Curves Qusaiba,
Q Dawaa
a and DG statiion O/C (3-pha
ase
fault&E/F F) Relays With
h Two DG

1173
 
Eng. &
&Tech.Journall, Vol.34,Part (A),
(A No.62016 Protection
n Coordinatio on with Distriibuted Generration
in Electtrical System of Iraqi Distrribution Grid
d

1‐Reelay : E/F with insst. DG station                           3‐Rellay: E/F with inst. Al‐Dawaa5‐ Rela
ay : E/F with inst.. Qusaiba 
2‐ R
Relay; 3 ph /Fwith h inst. DG station   4‐Relay: 33 ph /F  with inst. Al‐Dawaa 6  ‐‐ Relay : 3 ph /F  w ba
with inst. Qusaib
 
Figgure (9) :Time Current Charracteristic (TC C) Curves for Qusaiba, Daw
waa and DG staation O/C(3-ph
hase
fault&
&E/F)Relays W With instantaneeous (With Tw
wo DG)

Taable(14) Parammeters and Opeerating Timefoor Over Curreent Relay (3-phhase fault& insstantaneous) in
n C.B. of
Qussaiba S.S. Withh Two DG
TMMS = 0.6 P =2
PS Ipu = 180A C.T.ratio = 600/5
Fault locations L-L-L FFault current((A) Opeerating Time(ssecond)
110%% Ipu 198 44.025
Zaaf. B.B. 6055 1.380
Qusaibaa B.B. 10725
5 1.380
Instaantaneous : Ipu = 5400A Delay Timme = 0.05 sec.

Tablee(15) Parameteers and Opera

ating Timefor O
Over Current Relay
R (E/F & instantaneous)
i ) in C.B. of Qusaiba
S..S. With Two DG
D
TMMS = 0.2 P = 0.4
PS Ipu = 48A C.T.ratio = 600/5
Fault locatiions L-G Fault cuurrent(A) Operatingg Time(second))
110% Ipuu 52,8 1 4,675
Zaaf. B.B
B. 566 00.554
Qusaiba B.B.
B 604 00.539
Instaantaneous : Ipu = 300A Delay Timme = 0.05 sec.
Tablee(16) Parameteers and Operating Timefor O
Over Current Relay
R (3-phasee fault& instanntaneous) in C.B.
C of Zaaf.
S.S With Tw
wo DG
TMMS = 0.3 P =2
PS Ipu =240A C.T.ratio = 600/5
Fault locations L-L-L F
Fault current((A) Opeerating Time(ssecond)
110% % Ipu 264 22.012
Zaaf. B.B. 6055 0.680
Qusaibaa B.B. 107255 0.680
Instaantaneous : Ipu = 3000A Delay Tim me = 0.05 sec.

1174
 
Eng. &Tech.Journal, Vol.34,Part (A), No.62016 Protection Coordination with Distributed Generation
in Electrical System of Iraqi Distribution Grid

Table(17) Parameters and Operating Timefor Over Current Relay (E/F & instantaneous) in C.B. of Zaaf. S.S With
Two DG
TMS = 0.1 PS = 0.3 Ipu = 36A C.T.ratio = 600/5
Fault locations L-G Fault current(A) Operating Time(second)
110% Ipu 39.6 7.337
Zaaf. B.B. 566 0.247
Qusaiba B.B. 604 0.241
Instantaneous : Ipu = 240A Delay Time = 0.05 sec.
Table(18) Parameters and Operating Timefor Over Current Relay (3-phase fault& instantaneous) in C.B. of DG
Station Near Zaaf. S.S.
TMS = 0.4 PS = 4 Ipu =480A C.T.ratio = 600/5
Fault locations L-L-L Fault current(A) Operating Time(second)
110% Ipu 528 29.350
DG station 5954 1,084
Dawaa B.B. 6055 1,077
Instantaneous : Ipu = 3000A Delay Time = 0.06 sec.

Table(19) Parameters and Operating Timefor Over Current Relay (E/F & instantaneous) in C.B. of DG Station
Near Zaaf. S.S.
TMS = 0.2 PS = 0.3 Ipu = 36A C.T.ratio = 600/5
Fault locations L-G Fault current(A) Operating Time(second)
110% Ipu 39.6 14.675
DG station 565 0.495
Dawaa B.B. 566 0.494
Instantaneous : Ipu = 240A Delay Time = 0.06 sec.

Figures (10), (11) represented the (TCC) curves, the increases in short circuit current after two
DG are embedded in 33KV network, for Qusaiba over current relay(I.D.M.T.) in 3-phase faults,
lead to (TMS) increase to (0.6) compares to(0.3) in case without DGand for (E/F) increase to
(0.2) compares to(0.1)for without DG. But the pickup current will be decrease to(180A) with
two DGs compares to (660A) without DG, because DGs supplied a large part of loads and the
load current will be little from utility source lead to mal-coordination between primary and
backup protection.
Figure (12) above represented the (TCC) curves with coordination curves through
assuming current multiplier is (3.67) to re-coordinate the over current protection relay of
Qusaiba circuit breaker for 33KV feeder supplied Al-Zaafraniya loads, the pickup current will
increased to be the same as in case without DG (660A) instead of (180A) with two DG. But, we
should pay attention to the pickup current of instantaneous time of operation situation. It must
be set as (1473A) instead of (5400A) taking into consideration the current multiplier above to
stay on the same setting.
Figure (13) represented the (TCC) curves after re-coordination, it is clear the over current
protection relays (3-phase faultand E/F) for Zaafraniya 33KV circuit breaker will be operate fist
when three phase faults and earth faults occur in the zone between Qusaiba and Zaafraniya and
then Qusaiba over current relay will be operate to separate the feeder and clear the fault

CONCLUSIONS
The conclusions obtained from this work can be summarized as follows:
1. In radial distribution system, utilization many DGs with optimal number, location and size
are better than use one DG alone.
2. To choose the optimal locations, there are geographical and environmental constraints depend
on the nature of electric distribution network and its location with respect to the city center,
these constraints effect on the availability of appropriate area to construct DG power plants.

1175
 
Eng. &Tech.Journal, Vol.34,Part (A), No.62016 Protection Coordination with Distributed Generation
in Electrical System of Iraqi Distribution Grid

3. DGs penetrated will effect on over current protection relays (phase and ground) settings
during disturbance, like sudden increase in loads or during three phase and earth fault issues
which leads to increase in short circuit current and then mal-coordination between primary
and backup over current relays.
4. Maximum load current from utility supply reduced when DGs embedded to the network to
the fact that a great part of loads will be powered from DG, the pickup current will effect in
the utility side which cause intersection in time of operation, so it is necessary to reset the
affected over current relays, the optimal setting in some situations is take the pickup current
in case without DG.
5. In general the Time Setting Multiplier (TMS) must be increase after DG penetrated.
6. Because DG station almost be close to load bus substation, the parameter of overcurrent
relays will be convergent.
7. Concentricity on activate instantaneous feature in over current protection relays to avoid the
risks from high current faults.

REFERENCE
[1]Deepak Pandey, Jitendra Singh Bhadoriya, “Optimal Placement & Sizing of
DistributedGeneration (DG) to Minimize Active Power Loss Using Particle Swarm
Optimization (PSO) ”, International Journal of Scientific & Technology Research Vol.
3,Issue 7, PP. 246-254, July, 2014.
[2] O. Amanifar, M.E. HamedaniGolshan, “Optimal Distributed Generation Placement
andSizing for Loss and THD Reduction and voltage Profile Improvement in Distribution ”,
Journal on Technical and Physical Problems of Engineering (IJTPE), Published by
Internationa Organization on TPE (IOTPE) ,Issue 7,Vol.3,No.2,PP.4753,June,2011.
[3] ArashAfrazet al., “Optimal Sizing and Sitting in Radial Standard System using PSO”,
American Journal of Scientific Research, Issue 67, pp. 50-58, 2012.
[4] Helen Cheung et al.,“Investigations of Impacts of Distributed Generations on Feeder
Protections”, IEEE Indicon, 2009.
[5] HadiZayandehroodi et al., “A Comprehensive Review of Protection Coordination Methods
in Power Distribution Systems in the Presence of DG”, Paper, the university of Islamic
Azad, Iran, and, the university of kebangsaan, Malaysia (UKM) , 2011.
[6] Vikas Singh Bhadoria, et al., “A Review on Distributed Generation Definitions and DG
Impacts on Distribution System”,Department of Electrical Engineering ,Gautam Buddha
University and National Institute of Technology, India, Paper, 2011.
[7] J.A. Pec¸as Lopes, et al., “Integrating Distributed Generation into Electric Power Systems: a
Review of Drivers, Challenges and Opportunities”, Electrical Power and Energy Systems
(ELSEVIER), pp.1189-1203, 2007.
[8] AngeloL'Abbate, et al., “Distributed Power Generation in Europe: Technical Issues for
Further Integration”,scientific and technical Report for European Commission
JointResearch Centre Institute for Energy, Netherlands, 2008.
[9] Carmen L.T. Borges, Djalma M. Falcao, “Optimal Distributed Generation Allocation for
Reliability, Losses, and Voltage Improvement”, Electrical Power and Energy
Systems(ELSEVIER), pp. 413-420, 2006.
[10] AngelFernándezSarabia, “Impact of Dstributed Generation on Distribution System”, M.Sc.
Thesis, Department of Energy Technology, Aalborg University, June, 2011.
[11] Zaid Tariq Hasan, “Stability Enhancement Using Distributed Generation in Al-Najaf Al-
Ashraf Electrical Network”, M. Sc. Thesis, Electrical Engineering department, University
of Technology, Baghdad, June, 2014.
[12] Stephen J. Conner, BEng, “Distributed Dispatching for Embedded Generation”, Ph.D.
Thesis, Department of Mechanical Engineering, University of Strathclyde, Glasgow, U.K.,
August, 2003.

1176
 
 Eng. &Tech.Journal, Vol.34,Part (A), No.62016 Protection Coordination with Distributed Generation
in Electrical System of Iraqi Distribution Grid

[13] Nitin Singh, “Optimal Sizing and Placement of DG in Radial Distribution Network using
Sensitivity Based Methods”, M. Sc. Thesis, Electrical and Instrumentation Engineering
Department, Thapar University, India, July, 2014.
[14] ManasRanjanNayak, “Optimal Feeder Reconfiguration of Distribution System with
Distributed GenerationUnits using HC-ACO”, International Journal on Electrical
Engineering and Informatics ‐ Vol. 6, No. 1, pp. 107-128, March, 2014.
[15] Hani M. Baniodeh, et al., “The Pest Location And Size Of Distributed Generation (DG)
inA Radial Distribution Systems”,International Journal of Advanced Research in
Electrical,Electronics and Instrumentation Engineering, Vol. 2, Issue 10, October, 2011
[16] Juan M. Gers, Edward J. Holmes, “Protection of Electricity Distribution Networks ”, 2n
edition, The Institution of Electrical Engineers, London, United Kingdom, 2004.
[17] C. R. Bayliss, B. J. Hardy, “Transmission and Distribution Electrical Engineering”, Third
edition,Elsevier Ltd., Book, 2007.
[18] H. C. Mehta, “Power Engineering Training Course on Fault Current Calculations, Relay
Setting and Relay Co-ordination”Power-Linker Training Center, July, 2004.

1177
 
Eng. &Tech.Journal, Vol.34,Part (A), No.62016 Protection Coordination with Distributed Generation
in Electrical System of Iraqi Distribution Grid

APPENDIX ( A )
Data Base for South of Al-Rusafa 33KV Distribution NetworkTaken From Ministry of Electricity
(MOE) in Iraq
From To R X PL QL From To R X PL QL
Bus Bus P.U. P.U MW Mvar Bus Bus P.U. P.U MW Mvar
1 2 0.00 0.18 0 0 8 37 0.06 0.05 11.75 7.282
1 3 0.00 0.18 0 0 8 38 0.034 0.0643 9.75 6.043
1 4 0.00 0.18 0 0 8 39 0.03 0.025 9.5833 5.9392
1 5 0.00 0.18 0 0 9 40 0.06 0.05 5 3.1
1 6 0.00 0.18 0 0 9 41 0.03 0.025 6.75 4.183
1 7 0.00 0.18 0 0 9 42 0.034 0.0643 2.6 1.611
1 8 0.00 0.18 0 0 9 43 0.035 0.03 1 0.62
1 9 0.00 0.18 0 0 10 44 0.035 0.03 7 4.34
1 10 0.00 0.18 0 0 10 45 0.036 0.069 19.67 12.19
1 11 0.00 0.18 0 0 10 46 0.035 0.03 1 0.62
1 12 0.00 0.18 0 0 11 47 0.017 0.032 9 5.58
1 13 0.00 0.18 0 0 12 48 0.0705 0.01355 6.667 4.138
1 14 0.00 0.18 0 0 12 49 0.027 0.0514 15.333 9.5
1 15 0.00 0.18 0 0 12 50 0.058 0.049 6.517 4.04
1 16 0.00 0.18 0 0 13 51 0.017 0.032 16.28 10.09
1 17 0.00 0.18 0 0 13 52 0.0096 0.0184 10.667 6.61
1 18 0.00 0.18 0 0 14 53 0.035 0.03 9.833 6.094
1 19 0.00 0.18 0 0 14 54 0.0955 0.01842 15.667 9.71
1 20 0.00 0.18 0 0 14 55 0.176 0.3384 18.75 11.62
2 21 0.0174 0.033 23.75 14.72 14 56 0.111 0.2134 10.83 6.712
3 22 0.0183 0.035 10.58 6.557 56 57 0.066 0.126 0.833 0.516
3 23 0.0082 0.0156 21.58 13.374 15 58 0.116 0.099 9.65 5.98
3 24 0.0174 0.033 16.83 10.43 15 59 0.045 0.087 13.083 8.11
3 25 0.0082 0.0156 19.25 11.93 59 60 0.1825 0.2443 5.167 3.2
4 26 0.0275 0.0523 6.083 3.77 60 61 0.079 0.0115 14.67 9.092
5 27 0.0232 0.02 10.533 6.53 15 62 0.098 0.1884 27.75 17.2
5 28 0.0327 0.063 6.583 4.08 15 63 0.116 0.099 0.583 0.361
5 29 0.039 0.074 10.5 6.5 16 64 0.035 0.03 8 5
5 30 0.0087 0.0165 9.333 5.78 17 65 0.29 0.247 11.333 7.024
5 31 0.0232 0.02 4.583 2.84 17 66 0.174 0.148 4.6667 2.892
6 32 0.0246 0.047 6.75 4.18 18 67 0.13 0.25 13.583 8.42
6 33 0.06 0.05 8.667 5.37 19 68 0.039 0.075 9.333 5.78
6 34 0.0232 0.02 11.62 7.21 19 69 0.078 0.15 5.6667 3.512
7 35 0.006 0.011 6.167 3.822 20 70 0.185 0.35 18.467 11.444
7 36 0.0145 0.0276 12.333 7.64 20 71 0.185 0.35 11.92 7.387

1178
 
Eng. &Tech.Journal, Vol.34,Part (A), No.62016 Protection Coordination with Distributed Generation
in Electrical System of Iraqi Distribution Grid

South of Al-Rusafa 33kV distribution network

APPENDIX ( B )
ANSI/IEEE and IEC constants for standard overcurrent relays
Curve description Standard α β L
Moderately inverse IEEE 0.02 0.0515 0.114
Very inverse IEEE 2.0 19.61 0.491
Extremely inverse IEEE 2.0 28.2 0.1217
Inverse CO8 2.0 5.95 0.18
Short-time inverse CO2 0.02 0.0239 0.0169
Standard inverse IEC 0.02 0.14 0
Very inverse IEC 1.0 13.5 0
Extremely inverse IEC 2.0 80.0 0
Long-time inverse UK 1.0 120 0

1179
 
 Eng. &
&Tech.Journall, Vol.34,Part (A),
(A No.62016 Protection
n Coordinatio
in Electtrical System

APPPENDIX ( C )
Relays Tyypes and Specifications
For Q
Qusaiba: Phase Fault (50 , 51) SIEMENS 7S Sj600 IEC I, E/F (50N , 51
For D
Dawaa, Zaafranniya and DG sttations :Phaase Fault (50 , 51)
5 ALSTOM
AL
LSTOM P120 I

A
APPENDIX( D)
D

1‐Relaay : E/F DG station     3‐Relay: E/F  Al‐Z Zaafraniya 5‐ Relay :: 3 ph /F   Qusaiba
2‐ Rellay: 3 ph /FDG sta on                  4‐Relay: 3 ph /F  Al‐Zaafraaniya6‐ Relay : E/F Q

ure (10) :Tim

Figu haracteristic (TCC) Curv
me Current Ch ves Qusaiba, Z
phase
p fault& E/F) Relays Two
T With DG

1‐Relay : E/F with inst. D
DG station                 3‐Relay: E/F w
with inst. Al‐Zaafran
2‐ Reelay: 3 ph /Fwith insst. DG sta on                   4‐Relay: 3 ph //F  with inst. Al‐Zaaf

Figurre (11) :Time Current

C Characteristic (TCC
C) Curves for Qusaiba,
Q Zaafr
fault&
&E/F)Relays W
With instantaneeous (With Tw

1180

View publication stats