Earth Fault Distance Computation Methods
Earth Fault Distance Computation Methods
Earth Fault Distance Computation Methods
par
tme
nto
fEl
ect
ric
alE
ngi
nee
ringa
ndAut
oma
ti
on
Mo
al
t
eart hfaul ts in me dium vo ltage( MV)
h
-
D
d Rafi Adz
Computation Me tho ds
D1
dist ribut io n fe e de rs.Singl eph aset oe art h
89
faul tis th emo stc o mmo n faul tt ypein po w er
/
20
dist ribut io n ne tw o rks.I nF inland and in
1
4
Basedo
nT ransients in
mo sto fc o unt rie s in Euro pet heMV
man
ne tw o rks areo pe rat edw ithne ut raliso l
ated
or re so nante art h e d syst e ms.T h ew ay the
PowerDistribution Syste
ms
Eart
ne ut ralis c o nne c te d is impo rtantsinc eit
w illaffe c tth et ypeo f pro t e ction syst em
hFaul
ne ede d and po w e r syst e mc o mpo ne nt
tD
se l
e c t
io n.I n iso lat e d ne ut ralne tw orks e arth
ist
faul tc urre nt s arere lative l y smal lc o mpare d
anc
w ithe arth e d syst e ms.T h einit ialt ransie nts
eCo
ofe art hfaul t s areimpo rt antfo r faul t Mo
hdRa
fiAd
zma
n
mput
locat io ne spe c ially in une art h e d and
co mpe nsat ed ne ut ralne tw o rk.T h e
atio
elec tro magne tict ransie ntpro videt h e
n Me
avail abl einfo rmat io n fo re st imat ing faul t
t
locat io n.T his disse rtat io n pro po se s fault
ho
ds Base
locat io n algo rit h ms base do ne art hfaul t
transie ntsignal s me asure d atse co ndary
sideo f MV and L V dist ribut io nt ransfo rme rs
donT
in orde rt ol oc at ea singl el inet oe art hfaul t
ransie
in une arth e d and c o mpe nsat ed ne ut ralMV
dist ribut io n ne tw o rks.
nts in P
owe
rDist
ribut
ion Syst
ems
I
S BN978-
95 2
-60-
5973-0(pri
nted) BUSI
N E
S S+
9HSTFMG*afjhda+
I
S BN978-
95 2
-60-
5974-7(pdf) ECONOM Y
I
S SN-
L1799-4934
I
S SN1799-
4934(pri
nted) ART+
I
S SN1799-
4942(pdf) DESI
G N+
ARCH I
TECT
U RE
A
Aa
lt
oU ni
versi
ty
a
Sc
hoolofElect
ric
alE ngi
ne e
ring SCI
ENCE+
lt
oU
Depa
rtmentofElect
rica
lE ngi
neeri
nga
ndAut
oma
ti
on T
EC HNOL
O GY
www.
aalt
o.f
i
ni
CROSSOV
ER
ve
DOCTORAL
r
DOCTORAL DI
S SE
R T
ATI
O NS
sit
DI
S SE
R T
ATI
O NS
y
Aal
toUni
vers
it
ypubl
icat
ionseri
es
DOC T
O RALDISSE
R TATIONS1 89
/20
14
Mo
hdRa
fiAd
zma
n
Ad octoraldi
ss ert
ati
onc o
mp let
e dforthed eg
reeo fDoct
oro f
Science( Tec
h no l
ogy)t
ob ed ef
e nded,wit
ht hepermi
ssi
o no ft
he
Aalt
oU nivers
ity,S choo
lofE l
ectric
alE ngi
neeri
ng,ata p
ub l
ic
examinationheldattheLectureh al
lS3o ft
h eschoolo
n8D ecember
2014at1 2no on.
Aa
lt
oU ni
versity
Sc
hoo lofE l
ec t
ric
alE ngineer
ing
Depa
rtme nto fElect
ricalEnginee
ringandAut
oma
ti
on
P
ow erS y
s t
ema ndH i
g hV o
lt
a geEngine
eri
ng
Super
vis
ingp ro
fessor
Pro
fes
sorM at
tiLeh
tonen
Prel
iminarye x
ami ne rs
Pro
fessorM ansour H.A bd
el-
R ah
man
Uni
versit
yofM anso ura
Egypt
Ass.P rof.Dr.I
ng .Pe
ter Brac
iní
k
UniversityofŽil
ina
Sl
ov ak i
a
Dr.SeppoHänni
nen
VTTT e
chni
calResearc
hCe
ntreo
fFi
nland
F
inland
Aal
toUni
vers
it
ypubl
icat
ionseri
es
DOC T
O RALDISSE
R TATIONS1 89
/20
14
©Mo
hdRaf
iBi
nAd
zman
I
S BN9 78-952-6
0-5973-0( pri
nted)
I
S BN9 78-952-6
0-5974 -
7( pdf)
I
S SN-L1 7
9 9-
4934
I
S SN1 799-4934(pri
nted)
I
S SN1 799-4942(pdf)
ht
tp:
//urn.
fi
/UR N:
ISBN :
97 8-
9 52-
60-59
74-
7
Uni
graf
iaOy
Hel
sink
i201
4
F
inl
and
A
bst
rac
t
A
alt
oUni
ver
sit
y,P
.O.Bo
x11
000
,FI
-00
076A
alow
t ww.
aal
to.
fi
Autho r
Mo hd Rafi Adzman
Nameo fthed oc t
o raldi
s sert
a t
ion
EarthF aultDistanc eCo mput atio
n Metho ds Base
donTransie
nts in Po
w e
rDist
ribut
ion
Syst
e ms
Publishe rScho olo f El
ectric
alEngineering
UnitD epartme nto f Elect
ricalEnginee
ring
Se
riesAal
toUnive
rsit
y publication se
riesDOCTORALD I
SSERT ATI
ONS 1
89/
201
4
F
iel
do fr
esea
rc hPo wer Syst
e m and HighVol
tageEngine
ering
Ma
nus
cri
pts
ubmi
ttd25 August20
e 14 Da
teo
fth
ede
fe e8 D
nc ec
ember 2014
P
ermi
ssi
ont
opub
li
shg
rant
ed(
dat
e)4 No
vembe
r 20
14 La
nguageEnglish
Mo
nog
rap
h Ar
ti
cled
iss
ert
ati
on(
summa
ry+ o
rig
ina
lar
ti
cle
s)
A b s tra c t
T h emo stc o mmo n faul tt ypein MV dist ribut io n ne tw ork is singl el inet oe art hfaul t.T h e
init ialt ransie nt so fe art hfaul ts areimpo rtante spe cial ly fo r une art h e d and c ompe nsat ed
ne ut ralne tw o rks.T h ee art hfaul tt ransie ntsignal sc onsisto f many diffe re ntfre que nc y
co mpo ne nt s, w h ic hre sul tfro mc h arging and disc h arging o ft h ene tw o rk c apac it anc es.T h e
transie ntc o mpo ne nt s pro videval uabl einfo rmat io n fo r faul tlo cat io n purpo se s.T h ec h arging
co mpo ne nth as h igh e r ampl itudeand l o w er fre que nc yt h an t hedisc h argec ompo ne ntand
h enc eis mo resuit abl et obeuse d fo r faul tl o c at io n purpo se s.
Int his t h e sis, w edisc uss al go rit h ms t ol o cat ean e art hfaul tin une art he do rac o mpe nsat e d
ne ut ralMV ne t w o rks using t h einfo rmat iono ft h eme asure dt ransie ntsignal .T h ene tw o rks
co nside re d areassume dt oberadial lyo pe rat ed and t h ey aremo de led using El e ct ro magne tic
T ransie ntP ro gram- Al t
e rnat iveT ransie ntP ro gram ( EMT P - AT P) .F ivet ype so f faul tl o c atio n
algo rith ms h avebe en de ve lo pe dw h ic harec al le d ge ne ralmo de l(GM)al go rith m, e xac tmo de l
(EM)al go rit hm, c o nt inuo us w ave le tt ransfo rm ( CWT )base d me t ho d, mul tiplere gre ssio n
anal ysis ( MRA)base d me th o d and art ific ialne uralne tw o rk ( ANN) .GM al go rith m is de ve lo pe d
base do n a simpl ifie d mo de lo f symme t ric alc o mpo ne nt sw h ileEM al go rith m is de ve l
o pe dw ith
exac t" pi" -mo de lo f symme t ric alc o mpo ne nt s.Bo thal go rit h ms ut ilizet h efre que nc yo fc h arging
transie ntt oe stimat et h efaul tdist anc e .CWTbase d al go rit h m re quire s bo thvo ltageand
curre nto ft ransie ntsignal st oe stimat et h efaul tpat hinduc tanc e.MRA and NN al go rit h ms
w e rede ve l o pe d using t h et ransie ntsignalme asure d fro mt h ese co ndary sideo ft h eMV/ L V
dist ribut io nt ransfo rme r.I n addit io n, an al go rit h mt ofind a c o rre c tpat ht ow ards t hepo sitio n
of faul tin ne tw o rk w h ichh as many branc h e s is pre se nt e d.
T h ere sul t s fro m int ensivesimul at io ns and e xpe rime nt s in ac tualdist ribut io n ne tw o rks are
alsopre se nt e d in t his t he sis.T h ere sul ts areanal yz ed using signalpro c essing t e c hnique s.T h e
algo rith ms appl yc o nt inuo us w ave le tt ransfo rm ( CWT )tol o cat et h edo minantc h arget ransie nt
fre que nc y and e xt rac tt h espe cificc o effic ie ntc o rre spo nding t ot h ec harget ransie ntfre que nc y.
Int h is th e sis, t hepro pe rtie so fH il be rtt ransfo rmat io n( H T )areuse dt oe stimat et h edamping
atte nuat io no ft het ransie ntsignal .
F inal lyt h epe rfo rmanc eo ft h epro po se d faul tl o cat io n al go rit hms is e val uate d and t h ere sul t s
arec ompare d.Base do nt h esimul at io n re sul t s, itis fo und t h att h epro po se d algo rit hms w o rk at
a re aso nabl el e ve lo f ac curac y.T here sul t s fro m re ale xpe rime ntdat a sh o wt hatbo thCWTand
GM al go rit h ms h avea c omparabl ere sul t .
KeyworsEart
d hfault, faul
tl ocation, dist
ribut io
n network, co ntinuous w avel
ettransfo
rm,
chargetransie nt, neuralne tw ork, multipl
eregre ssion analysis, ATP-
EMT P
I
S BN(pr
int
ed )978- 952- 60 -5973-0 ISBN( pdf)978- 952-60 -5974-7
I
S SNL1
- 799-4934 I
S S N( pr i
nted )1799-4934 I
S S N( pd)1
f 799-4942
L
ocat
iono fp
ubl
ish
erHe
lsinki Loca
ti
o nofp r
int
ingH el
sinki Year2014
Pa
ges234 n h
ur t
tp:/
/urn.
fi/U RN:
ISBN:978-
952-60
-5974-
7
Preface
First of all, thanks and praise is due to almighty Allah, The Most Gra-
cious, The Most Merciful, who guide me in every step I take.
This work would have been impossible without the help and guidance
of several people, whose contribution I would like to acknowledge. First
of all, I would like to express my deepest gratitude to my supervisor Pro-
fessor Matti Lehtonen, for accepting and giving me this wonderful re-
search project. His supervision and his guidance both helped me to chan-
nel and specify the discussed ideas and at the same time provided much
appreciated freedom and support to explore new ways and concepts. I
am also thankful to the pre-examiners Professor Dr. Mansour H. Abdel-
Rahman and Assoc. Professor Dr. Ing. Peter Braciník for reading the
thesis manuscript and I appreciate their valuable comments.
5
Preface
6
Contribution of the author
The author has developed the methodology presented and used in the
analyses and made all the analyses himself. The results and conclusions
are the sole work of the author.
7
Contribution of the author
8
Contents
Preface 5
Contents 9
List of Tables 13
List of Figures 17
List of Symbols 25
List of Abbreviations 31
1. Introduction 33
1.1 Background and Problem Statement . . . . . . . . . . . . . . 33
1.2 Research Objective . . . . . . . . . . . . . . . . . . . . . . . . 35
1.3 Research Contribution . . . . . . . . . . . . . . . . . . . . . . 36
1.4 Organization of Thesis . . . . . . . . . . . . . . . . . . . . . . 37
9
Contents
10
Contents
References 195
11
Contents
12
List of Tables
6.1 Test results of single line to earth fault for 90◦ fault incep-
tion angle using the WT algorithm. . . . . . . . . . . . . . . . 138
6.2 Test results of single line to earth fault for 45◦ fault incep-
tion angle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 139
6.3 Data generation for MRA and NN formulation. . . . . . . . . 150
6.4 Data generation for testing the MRA and NN algorithm. . . 150
13
List of Tables
14
List of Tables
15
List of Tables
16
List of Figures
17
List of Figures
18
List of Figures
19
List of Figures
20
List of Figures
21
List of Figures
7.1 Compensated neutral network test diagram for test 1-2 . . . 164
7.2 Compensated neutral network test diagram for test 3-4 . . . 165
7.3 Isolated neutral network test diagram for test 5-8 . . . . . . 166
7.4 Recorded signals of single line to earth fault at 0.46 km. . . 168
7.5 Recorded signals of single line to earth at 10.40 km. . . . . . 169
7.6 Recorded signals of single line to earth fault at 14.20 km.
Test number 5 and 6. . . . . . . . . . . . . . . . . . . . . . . . 170
7.7 Recorded signals of single line to earth fault at 14.20 km.
Test number 7 and 8. . . . . . . . . . . . . . . . . . . . . . . . 171
7.8 Example of output from pre-processing of test number 1
(0.76 km) transient signal using CWT analysis. . . . . . . . . 172
7.9 Example of output from fault distances algorithm based on
CWT for test number 1. . . . . . . . . . . . . . . . . . . . . . . 173
7.10 Results of fault distance estimation (lf ) using CWT algorithm.173
7.11 Lounais-Suomen Sahko OY network test diagram. . . . . . . 175
7.12 Recorded signals of single line to earth fault at 25.40 km.
Test number 1 and 2. . . . . . . . . . . . . . . . . . . . . . . . 176
7.13 Recorded signals of single line to earth fault at 25.40 km.
Test number 3 and 4. . . . . . . . . . . . . . . . . . . . . . . . 177
7.14 Recorded signals of single line to earth fault at 25.40 km
(test number 5) and 36.0 km (test number 6). . . . . . . . . . 178
7.15 Recorded signals of single line to earth fault at 36.0 km. . . 179
7.16 Example of CWT scalogram and energy percentage distri-
bution of test number 4 (25.4 km) transient signal. . . . . . . 180
7.17 Example of CWT extracted coefficients, FFT frequency spec-
trum of ectracted coefficent and, voltage and current linear
decaying of extracted CWT coefficient of test number 4 (25.4
km) transient signal. . . . . . . . . . . . . . . . . . . . . . . . 181
7.18 Fault distances calculated using CWT algorithm. . . . . . . . 182
7.19 Schematic diagram of Czech distribution network. . . . . . . 183
7.20 Photo taken during the experiment. . . . . . . . . . . . . . . 184
7.21 Recorded signals of experiments 1. . . . . . . . . . . . . . . . 186
7.22 Recorded signals of experiments 2. . . . . . . . . . . . . . . . 186
7.23 Recorded signals of experiments 3. . . . . . . . . . . . . . . . 187
22
List of Figures
23
List of Figures
24
List of Symbols
a Positive scale
am
0 Scaling (dilation) parameters
A Ampere
FFT amplitude of extracted CWT coefficient
Amplitude
b Time shifting factor
Time position
C Capacitance
Coefficient
Phase to phase capacitance
C(a, b) Wavelet coefficient
Ce Phase to earth capacitance
Ceq Equivalent capacitance
Cj,0 Zero sequence capacitance of sound line
Cj,p Positive sequence capacitance of sound line
Cl,n Capacitance per unit length of negative sequence system
Cl,0 Capacitance per unit length of zero sequence system
C l,0′ Zero sequence capacitance of faulty line located behind
of the fault point (F)
Cl,p Capacitance per unit length of positive sequence system
Cl,p′ Positive sequence capacitance of faulty line located behind
of the fault point (F)
Cp Positive sequence capacitances of the network
C0 Zero sequence capacitances of the network
CV (a, b) Wavelet coefficient of voltage transient signals
CI (a, b) Wavelet coefficient of current transient signals
CE Phase to earth capacitance
ESC Sum of percentage energy coefficient associated to same scale
25
List of Symbols
f Frequency
fc Charge transient frequency(Hz)
Calculated transient frequency(HZ)
Natural frequency(Hz)
fd Discharge transient frequency (Hz)
fdc Calculated damped transient frequency(Hz)
fdm Estimated damped transient frequency (HZ)
fm Estimated transient frequency(Hz)
fnc Calculated undamped transient frequency(Hz)
fnm Estimated undamped transient frequency (Hz)
fs Sampling frequency(Hz)
f (t) Original signal
F Fault point
Fc Pseudo-frequency
F0 Center frequency
g Output of filter
i Phase current
îc Initial amplitude of charge transient
îe Uncompensated steady state earth fault current
I Current amplitude coefficient
Ic Charging current
Iω Complex wavelet coefficient of current
k Integer variable parameters
Sample number
Node
kam
0 Translation (time shift) parameters
l Fault distance
lcl Length of connected line
lf Fault distance
L Inductance
Total loads
Line
Lc Inductance of suppression coil
Ldt Inductance of distribution transformer
Leq Equivalent inductance
Lj,0 Zero sequence inductance of sound line
Lj,p Positive sequence inductance of sound line
Ll Fault path inductance
26
List of Symbols
27
List of Symbols
28
List of Symbols
29
List of Symbols
30
List of Abbreviations
AC Alternating current
AM R Automatic meter reading
AN N Artificial neural network
AT P Alternative transient program
CW D Choi-Williams distribution
CW T Continuous wavelet transform
DC Direct current
DF T Discrete Fourier transform
DOM ’Domanin’ distribution substation
DT Distribution transformer
DW T Discrete wavelet transform
DW F Discrete wavelet filter
EM Exact model
EM T P Electromagnetic transient program
FD Frequency-dependent
FL Fault location
FT Fourier transform
GM General model
HT Hilbert transform
HV High voltage
Im Imaginary part of a complex number
IEDS Intelligent electronic devices
ln Natural logarithm
LV Low voltage
M LP Multilayer perceptron
M RA Multiple regression analysis
MV Medium voltage
N OV E ’Nove Mesto’ distribution substation
31
List of Abbreviations
32
1. Introduction
This thesis focuses on single phase to earth faults in medium voltage (MV)
distribution feeders. Single phase to earth fault is the most common fault
type in power distribution networks. In Nordic countries, about 50% to
80% of the faults are of this type [1, 2]. In Finland, about 70% of the per-
manent faults that occur in rural MV networks are caused by animals
and weather conditions, like thunder, snow, icing, storm, hard wind and
fallen trees [3]. One of the solution to reduce an earth fault in overhead
line feeders is by replacing the bare conductor to covered conductor type
of overhead line. However, the solution requires a huge amount of invest-
ment. Based on the statistics survey done by [4], a typical feeder construc-
tion in MV distribution network was a mixture of both covered conductors
and bare overhead lines on the same feeders. The share of covered conduc-
tors per feeder varies between 10% to 50% and all the permanent faults
at covered conductor feeders were earth faults [4].
In 2010, a fault statistic studied by utility company in Dresden Germany
shows that 75% of all faults are earth faults [5]. In Finland, over 80%
of the annual outage costs of customers are due to to faults in medium
voltage distribution networks [6]. In rural overhead networks, earth fault
currents are typically 10 - 40 A [7]. The small earth fault current is an
advantage to system operation and is unlikely to cause damage to the
power system components. However, the detection and localization of the
fault is difficult due to the small current that reduces the sensitivity of
the conventional relay which is normally based on currents and voltages
at power frequency.
Basically, electricity supply system differs from one country to another.
33
Introduction
34
Introduction
will provide the available information for estimating fault location. Even
though the transient may cause problems to the conventional protective
schemes it may work properly under transient protective scheme. The
transient disturbance recorder can be installed at the monitored network
to support the conventional signal measurement that is only based on the
root mean square (RMS) voltage and current of the power frequency com-
ponent. Nowadays, fault localization and detection based on fault tran-
sient signals in transmission and distribution network has been a subject
of keen interest to utility and researchers [32–35].
With the upgrading of the distribution network towards the concept of
“SmartGrids”, more measuring and recording equipment such as smart
meter, microprocessor digital relay, digital fault recorder and multi utility
controller will be installed in low and medium voltage distribution substa-
tions. In order to adapt and apply the concept of distribution automation
system in the network, those equipments are monitored and time syn-
chronized. Thus, the realization of this kind of system will open more op-
portunity to control, optimize and enhance the reliability of current sys-
tem. This will also affect the development and enhancement of current
protection system and fault location algorithms in power networks.
The term fault location can be represented as the determination of the
faulty feeder, branch or line section. It also can be used as a general
term when we discuss about fault distance computation or estimation.
This dissertation proposes fault location algorithms based on earth fault
transient signals measured at secondary side of MV and LV distribution
transformers in order to locate a single line to earth fault in unearthed
and compensated neutral MV distribution networks.
In order to reduce the outage costs and to satisfy the quality of supply
requirements, fault location algorithms have to locate fault as soon as
possible and with an efficient way. This can be achieved by utilizing the
information of initial transients of earthed faults. The idea of using tran-
sient signal as a basis of developing fault location algorithm is not a new
thing. However, the formulation used to calculate the exact fault distance
may be different. In addition, the implementation and the execution of
algorithm might also be dissimilar. The main objective of this research
is to developed accurate and efficient fault location algorithms for locat-
35
Introduction
The contribution of the research can be divided into three main cate-
gories. In the first category, there are fault location algorithms based
on measured transient signals using general model(GM) of symmetrical
components, exact model (EM) of symmetrical components and continu-
ous wavelet transform (CWT) based algorithm. The proposed algorithms
mainly use measured transient signal from MV side of the faulty feeder,
for example from MV main substation. In GM and EM algorithm, 4 gen-
eral models (GM) and 3 exact models (EM) of symmetrical components
has been presented and tested using simulated earth fault. Due to the
limitation of information of line parameters in real networks, only GM2
and CWT algorithm were tested using real recorded data. Based on the
results, GM2 and CWT algorithms managed to produce an acceptable ac-
curacy in estimating the fault distance in a real network.
In the second category, were the algorithms based on earth fault tran-
sient signals measured from the secondary side of the MV/LV distribution
transformer (DT). Fault location algorithm based on multiple regression
analysis (MRA) and artificial neural network (ANN) was developed in this
category. In addition, an algorithm which can be used as indication to se-
lect the correct path to the fault position when earth fault happens in a
network that have several branches was developed.
In the third category, a technique to measure the damping factor of
charge transient signal was proposed, based on continuous wavelet trans-
form (CWT) analysis and Hilbert transformation. CWT analysis is a pow-
erful tool to visualize the non-stationary signal such as earth fault tran-
sient. In this case, the calculated damping factor is used to estimate the
undamped charge transient frequency which in turn will be used to esti-
mate the fault distance in the GM algorithm.
36
Introduction
37
Introduction
38
2. Transients of Single Phase to Earth
Fault
This chapter studies the transient signals which are generated from sin-
gle line to earth fault. In [36], a technical analysis of transient phenomena
associated with earth fault in a three phase system is given. According
to [36], the transients are caused by redistribution of conductor voltage
throughout the whole system, with the fault location as the origin of the
change. When single line to earth fault occurs, two different processes are
superposed [37–39]. Both charge and discharge transient processes start
at the same time but with different duration. Figure 2.1 shows a typi-
cal single line to earth fault transient in the voltage and current of the
faulted phase. The depicted signals were recorded in a 20 kV overhead
line network.
In this chapter, we provide general knowledge about the transients which
arise during the single phase to earth fault. In the following section, the
process of charge and discharge transient will be explained using lumped
model of two feeders and the earth fault at phase 3 of faulty feeder as
shown in Figure 2.2 and Figure 2.3.
39
Transients of Single Phase to Earth Fault
4
x 10
4
3 Ua
Ub
2
Voltage [V]
Uc
1
0
−1
−2
−3
−4
1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2
Sample Number [k] x 10
4
600
Ia
Ib
400 Ic
Current [A]
200
−200
−400
1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2
4
Sample Number [k] x 10
Figure 2.1. An example of the earth fault transient signal in MV overhead line network.
LOAD B
C L2
HV/MV C L3
ZT1 L1
ZT2 L2
N Ce Ce Ce
ZT3 L3
C L2
C L3
Rf Ce Ce Ce
If
Ic3A
Ic3B
40
Transients of Single Phase to Earth Fault
Equation (2.1) has to be modified and the discharge transient can be cal-
culated as:
v
fd = (2.2)
4(lf + lcl )
where, v is the velocity of the traveling wave, lf is the fault distance and
lcl represent the length of the connected line. Generally, the discharge
transient components are short duration high frequency components. The
discharge transient components are highly damped and disappear much
quicker than the charge transient component. Based on the real tests
conducted in [40], the amplitude of discharge component was typically 5
to 10% of the amplitude of charge component and their frequencies var-
ied through the range of 500 to 2500 Hz for distribution network. Based
on the Equation (2.1) and Equation (2.2), the fault distance significantly
affects the frequency of the discharge transient. In HV transmission net-
works, the discharge transient frequency is about 300 to 5000 Hz [16].
Due to the ground fault, the voltage of the sound phases rapidly increases
and subsequently increase the charge stored in their capacitances. In
Figure 2.1, rapid increase of two sound phases right after the inception
of the earth fault can be seen clearly. Figure 2.3 show the illustration of
charging process of healthy phases in isolated neutral network. In the
process, all the capacitive charging currents (Ic ) of the healthy and faulty
feeder have to flow over the fault location back to the transformer in phase
3. Assuming an earth fault is located at high voltage or medium voltage
substation, according to [37, 38], based on the simplified lump circuit as
shown in Figure 2.4, the frequency of charge component in undamped
condition can be calculated as:
41
Transients of Single Phase to Earth Fault
LOAD B
L2
HV/MV C L3
ZT1 L1
ZT2 L2
N Ce Ce Ce
ZT3 L3
Ic1B
Faulty Feeder Ic2B L1
LOAD A
C L2
C L3
Rf Ce Ce Ce
If
Ic1A
Ic2A
Ic1B
Ic2B
ve
îc = = ve ωc Ceq (2.6)
ωc Leq
1.5LT
2(C+Ce)
Ceq ωc
îc = îe (2.7)
3Ce ωf
42
Transients of Single Phase to Earth Fault
4
x 10
1
Voltage[V]
0.5
−0.5
Fault Inception
−1
0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1
0.1
3250 Discharge transient component
0.09
1625 0.08
Charge transient component
Frequency[Hz]
1083.33 0.07
0.06
812.5
0.05
650 0.04
541.667 0.03
0.02
464.286
0.01
406.25
0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1
Time(s)
(a) Composite earth fault transient voltage (above) and its CWT Scalogram(below)
100
Current[A]
50
−50
Fault Inception
−100
0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1
0.1
Discharge transient component
3250 0.09
1625 0.08
Charge transient component
Frequency [Hz]
0.07
1083.33
0.06
812.5
0.05
650 0.04
0.03
541.667
0.02
464.286
0.01
406.25
0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1
Time(s)
(b) Composite earth fault transient current (above) and its CWT Scalogram (below)
Figure 2.5. An example of composite earth fault transient signal and its CWT analysis.
The red box referred most energy of charge transient component while the
black box is for discharge transient component.
In Figure 2.5, it shows that the charge transient component has lower
frequency and longer duration than the discharged transient component.
43
Transients of Single Phase to Earth Fault
di
Lp + Rp i = v̂sin(ωf t + φv ) (2.8)
dt
44
Transients of Single Phase to Earth Fault
LOAD B
C L2
HV/MV C L3
ZT1 L1
N ZT2 L2 Ce Ce Ce
ZT3 L3
Ic1B
Faulty Feeder Ic2B L1
LOAD A
Lc C L2
C L3
Rf If Ce Ce Ce
IL
Ic1A
Ic2A
Ic1B
Ic2B
t
v̂ (− )
i= [sin(ωf t + φv − φ) − e τ sin(φv − φ)] (2.9)
Z
where:
Z= (Rp )2 + (ωf Lp )2
ωf Lp
φ = arctan
Rp
Lp
τ=
Rp
45
Transients of Single Phase to Earth Fault
3Rf
E 3LC Ce
Figure 2.7. Equivalent circuit model of transient of the suppression coil [7].
S 2 + Sδ + ω02 (2.10)
where:
1
δ=
6Rf Ce
ω0 = 3Lc Ce
46
3. Fault Location Methods Based on
Earth Fault Transient Signals
Transient signals have been used and utilized for many purposed in power
systems. Protection system for example is one of the areas in power sys-
tems where researchers and engineers have studied the characteristics
of the transient signal for the purpose of fault location and detection.
The main idea is to improve the efficiency and reliability that can bring a
higher performance of the system itself. Actually, using earth fault tran-
sients for the purpose of fault location is not a new idea. However, the
idea and the technique proposed by the researchers for their algorithms
might be different by one and another. With new technologies which of-
fer better, faster, cheaper and more robust services than their outdated
technologies, studies on fault location becomes more challenging and in-
teresting. In this chapter, an overview of several existing ideas related
to single line to earth fault location algorithms using earth fault initial
transient signals is presented.
47
Fault Location Methods Based on Earth Fault Transient Signals
48
Fault Location Methods Based on Earth Fault Transient Signals
In case of a single phase to earth fault, the fault path inductance of the
faulty line is a composed of a series connection of zero, positive and nega-
tive sequence inductance as follows [40]:
1
Lf = (Ll,p + Ll,n + Ll,0 ) · l (3.3)
3
where, l is the fault distance and Ll,0 , Ll,p , Ll,n are the inductance per unit
length of zero-, positive- and negative-sequence systems.
A fault location algorithm based on wavelet transform was presented
by [43] to estimate the fault path inductance of the faulty line length us-
ing transient signal. According to [43], the fault path inductance can be
estimated with the following equation:
1 Uω (k∆t, f )
Lf = imag (3.4)
ω Iω (k∆t, f )
where Uω and Iω are the complex wavelet coefficients of voltage and cur-
rent and ω is the angular frequency. According to [44], the algorithm first
determines the maximum wavelet coefficient of the current including the
amplitude, frequency and location of the wavelet. Using this frequency
with different time translations, the equivalent fault inductance can be
calculated. The 2 ms inductance interval, corresponding to 10 subes-
timates, is then determined with the smallest standard deviation. The
mean value of the inductance, calculated in this interval, is finally used
to determine the fault distance.
Researcher in [16] deals with an earth fault in subtransmission net-
works. His work related to transient based ground fault location using
charge transient. He proposed of to use discrete wavelet filter (DWF), the
FIR filter with the Gaussian wavelet and amplitude correction in filtering
specific coefficient of transient component. Instead of using DWT in many
frequency scales, DWF only uses a single scale where the transient is of
interest. However, the method still required a pre-processing technique
to locate the specific transient frequency. This is done by using discrete
Fourier Transform (DFT). The equation of DWF is given as follows:
(n − b)/a)2
2 −( )+jω0 (n−b)/a
DW F [n] = √ x [n] · e 2δ 2 (3.5)
a 2πδ n
49
Fault Location Methods Based on Earth Fault Transient Signals
di
u(t) = Ri(t) + L (3.6)
dt
where u, i, R and L are the voltage, current, resistance and the inductance
of the faulty phase, respectively. According to [45], since the differentia-
tion is sensitive to higher frequency noise, the solution of above equation
is done by using the integration technique. With samples of voltage and
current signal, inductance L is calculated by using three equally spaced
pairs of voltage and current using trapezoidal rule as follows:
50
Fault Location Methods Based on Earth Fault Transient Signals
Then, the fault distance is solved as reactance of the faulty line length:
U 2 I1 − U 1 I2
L= (3.13)
I1 I2 (s2 − s1 )
For proper computation, least square fitting methods require signal fil-
tering and pre processing. Equation (3.13) is the basic equation based on
first order line model. For better result, the model can be replaced with
a higher order model which includes capacitances at the close end of the
line [47].
U (ω)
Z= (3.14)
I(ω)
3Im(Z(ω))
l= (3.15)
(Ll,0 (ω) + 2Ll,p (ω))ω
where, Z(ω), U (ω) and I(ω) are the frequency (ω) component of impedance,
measured voltage and current, Ll,0 and Ll,p are the inductance per unit
51
Fault Location Methods Based on Earth Fault Transient Signals
Researchers have long been studying traveling wave transients for fault
location and fault detection purposes [49–53]. A fundamental concept of
using traveling waves when earth fault happens in the network is de-
scribed in Figure 3.1.
Fault location
Amplitude
T1
T2
Time
Figure 3.1. Propagation of traveling wave in power network due to earth fault.
νtd
x= (3.16)
2
52
Fault Location Methods Based on Earth Fault Transient Signals
Over the past few decades, a number of fault location algorithms based
on artificial neural networks (ANN) have been developed, including event
signal characteristics neural network approach [63–65], wavelet analy-
sis and neural network approach [66–68], radial basis function neural
network approach [69–71], support vector machine and neural network
approach [72–74] and fuzzy and neural network approach [75–77].
A transient based ANN algorithm for locating single line to earth fault
has been developed in [78]. In this work, Multilayer Perceptron (MLP)
with one hidden layer was used as ANN-structure and was trained with
Backpropagation method. The ANN-structure has only one hidden layer
which is sufficient for representing any given input-output transformation
[78, 79].
Two algorithms has been developed and presented in [14, 78]. The first
algorithm uses either the phase voltage or the phase voltage and current
samples as input data and the second algorithm uses harmonic compo-
nents of the neutral voltage transient as input data. The transient re-
quired filtering and pre processing technique before it can be used as an
input data. Both algorithms use fault distance as an output value. The
first algorithm was compared with conventional algorithms which are dif-
ferential and wavelet algorithm. According to the author [78], the perfor-
mance of ANN was comparable to the conventional algorithms.
53
Fault Location Methods Based on Earth Fault Transient Signals
54
4. Determination of Transient
Components with Signal Processing
Methods
55
Determination of Transient Components with Signal Processing Methods
where g(t) is the output of the filter, f (t) the original signal and T is the
period of fundamental frequency. In case of the fundamental frequency is
50 Hz, T is 20 ms.
Discrete Fourier transform (DFT) is probably one of the best known meth-
ods to transform the signal from time domain into frequency domain. It
has been used for a decade to estimate the amplitude and phase of the fun-
damental frequency and associated harmonics of current and voltage [81].
The Fourier transform (FT) and discrete Fourier transform (DFT) is de-
fined as [82]:
∞
F T (f ) = x(t) · e−j2πf t dt (4.2)
−∞
N −1 2πkn
−j
DF T [k] = x[n] · e N (4.3)
n=0
56
Determination of Transient Components with Signal Processing Methods
1 n − kam
DW T (f, m, n) = m f (k) · ψ ∗ ( 0
) (4.4)
a0 k am0
57
Determination of Transient Components with Signal Processing Methods
In Equation (4.5), Ψa,b (t) are derived from a chosen function Ψ(t) called
"mother wavelet" given as [43]:
1 t−b
Ψa,b (t) = Ψ( )dt (4.6)
|a| a
The mother wavelet Ψ(t) must be short and oscillatory, and it must have
zero average and decay quickly at both ends [43]. In Equation (4.6), a and
b are the parameters for positive scale and time position (time shifting
factor). The Ψ in Equation (4.5) is known as the complex conjugate of Ψ
and the output of CWT would be the wavelet coefficient denoted as C(a, b).
In literature, several mother wavelets have been used to analyze transient
signals [8, 20, 22, 23, 43, 57, 123]. In [16, 43, 122, 124, 125], Morlet wavelet
has been applied to earth fault transient signal for the purpose of fault
detection and fault location. Therefore, on the basis of their investigation
results the Morlet wavelet has been chosen in this study.
The expression of Morlet wavelet is given in Equation (4.7) as follows:
2 /2)
Ψ(t) = e(−t ej2πF0 t (4.7)
58
Determination of Transient Components with Signal Processing Methods
F0
Fc = (4.8)
a∆
N
ESC (a) = SC(a, b) (4.10)
n=1
59
Determination of Transient Components with Signal Processing Methods
50
I1
Current[A] 0
−50
50 100 150 200 250 300
49 0.16
46
43 0.14
40
37 0.12
34
31 X= 41
Scales a
Y= 24 0.1
28
25 Level= 0.16376
22 0.08
19
16 0.06
13
10 0.04
7
4
0.02
1
50 100 150 200 250 300
Sample Number[k]
(a) Simulated transient current (above) and its CWT scalogram (below)
10
I
X: 24 EP
9
Y: 8.182
8
Energy percentage[%]
7
6
5
4
3
2
1
0
0 5 10 15 20 25 30 35 40 45 50
Scales[a]
Figure 4.1. Example of the ouput of pre-processing of transient signal using CWT.
60
Determination of Transient Components with Signal Processing Methods
where:
A(t) = u2 (t) + ũ2 (t) (4.14)
ũ(t)
ϕ(t) = tan−1 (4.15)
u(t)
dϕ(t)
ω(t) = (4.16)
dt
where, uv (t) and ui (t) are the extracted wavelet coefficients of voltage and
current signal respectively, while ũv and ũi are their Hilbert transforma-
tions.
61
Determination of Transient Components with Signal Processing Methods
(i) Assume that we have obtained the interesting part of the transient
signal which is extracted CWT coefficient of charge transient com-
ponent as shown in Figure 4.2(c)
Current[A]
20
0
−20 Fault inception
50 100 150 200 250 300 350 400 450 500
40
37
34 0.15
31
Scales a
28
25 Charge transient component
22 0.1
19
16
13
10 0.05
7
4
1
50 100 150 200 250 300 350 400 450 500
Sample Number[k]
(a)
10
I
X: 24 EP
9
Y: 8.182
8
Energy percentage[%]
7
6
5
4
3
2
1
0
0 5 10 15 20 25 30 35 40 45 50
Scales[a]
(b)
200
I
CWT
150
100
Current[A]
50
−50
−100
−150
−200
0 50 100 150 200 250 300 350 400 450 500
Sample Number[k]
(c)
Figure 4.2. Example of the output from pre-processing of current transient signals using
CWT. (a) Simulated transient current (above) and its CWT scalogram (below),
(b) scale-dependent energy percentage distribution and (c) extracted CWT
coefficient
62
Determination of Transient Components with Signal Processing Methods
(iii) The damping attenuation (αm ) is estimated from the plot of instan-
taneous amplitude of the signal (Ax(t)) versus time as given in Fig-
ure 4.3 (a). In Figure 4.3 (b), the attenuation (αm ) is the approxi-
mated slope of the logarithm of the amplitude peaks (InAx(t)). A
straight line can be fitted through a plot of the log of the peaks. Ba-
sically, the damping can be estimated slightly after the time of the
earth fault inception. For example, in Figure 4.3 (b), the damping
factor is estimated from the time of 0.03 seconds to 0.06 seconds. If
the damping is linear, a first order least square polynomial curve fit-
ting method will be sufficient to fit the curve. Otherwise, the damp-
ing can be estimated by calculating the slope of the amplitudes of
two successive peaks and then averaging it.
50
40
30
20
A x[t]
10
0
−10
−20
−30
−40
−50
0.03 0.035 0.04 0.045 0.05 0.055 0.06
Time[s]
(a) Instantaneous current amplitude signal decaying characteristic (red) and its envelope
(blue)
4
lnA x[t]
−1
−2
0.03 0.035 0.04 0.045 0.05 0.055 0.06
Time[s]
Figure 4.3. Example of instantaneous current amplitude, envelope and the linear decay-
ing of the extracted CWT coefficient of charge component.
63
Determination of Transient Components with Signal Processing Methods
where ωdm , fnm and αm are estimated damped charge transient frequency
(rad/s) , estimated undamped charged transient frequency (Hz) and atten-
uation of the extracted wavelet coefficients, respectively.
64
5. Fault Distance Estimation Using
Transient Signals
65
Fault Distance Estimation Using Transient Signals
The main principle of the general model (GM) algorithm is based on esti-
mating the fault distance using its correlation with the single frequency
of the charging process during the earth fault. As explained in Chap-
ter 2, the charging transient is usually of higher amplitude and lower fre-
quency than discharge transient. In the algorithm, in order to estimate
the earth fault distance, we use the estimated earth fault charge transient
frequency and compare it with the calculated earth fault charge transient
frequency. The earth fault charge transient frequency is estimated using
CWT and FFT analysis as explained in previous chapter, while the calcu-
lated charge transient frequency is formulated using developed simplified
lumped parameter model.
The idea of using simplified model for locating an earth fault in power
system have been proposed by several researcher [40,42,129,130]. Accord-
ing to [40], in case of single phase to earth fault, an interconnection of its
phase sequence networks can be simplified into a composite forms of se-
quence networks as shown in Figure 5.1 and its description is given in Ta-
ble 5.1. In case of a compensated neutral network, the compensation coil
can usually be omitted since at the transient frequencies its impedance is
high compared to that of the network earth capacitances.
The general models developed by [40] as given in Figure 5.1 were used
to formulate the equation of charge and discharge transient frequency. In
order to examine the best choice of model to be used for estimating the
earth fault distance, four general models have been examined. The mod-
els are called as general model 1 (GM1), general model 2 (GM2), general
model 3 (GM3) and general model 4 (GM4). The details of the models
is explained in the following section. The advantages of the algorithm is
that only one measurement is needed per primary substation transformer.
The developed algorithm is used in simple iteration to compare the esti-
mated and calculated earth fault charge transient frequencies. Figure 5.2
shows the basic idea of estimating earth fault distance based on general
model of symmetrical component algorithm, the details flowchart of the
algorithms is given in Figure 5.7 at the end of this section. The next is
subsection concerned with a formulated equation that was developed from
the general model of symmetrical components of earth fault transient.
66
Fault Distance Estimation Using Transient Signals
L2
2Ll,p+Ll,0
Llp,n Ll,0
R2
C0
C2
LT Cp
C0 R1
2LT (Cp/2) C1
L1
L2
LDT Llp,n
R2
RDT R1p,n
C2
LDT L3
RT R1
Cp RDT C1 R3
LL RL
LT L1
LL RL L4 R4
Figure 5.1. The general model for ground fault transients. Positive and negative se-
quence networks (a), zero sequence network for isolated networks (b), the
composite connection for isolated networks without (c) and with damping in-
cluded (d), model for load impedances (e), the equivalent circuit with loads
included for positive and negative sequence (f) and the combinations of all
sequence networks (g) [40].
Composite
Remarks
components
C1 0.5Cp
C2 Co
L1 2LT
L2 2Ll,p + Ll,0
L3 2LDT
L4 2LL
R1 2RT
R2 2Rl,p + Rl,0 + 3Rf
R3 2RDT
R4 2RL
Notes:
RT and LT = resistance and inductance of sub-
station transformer, Rl,p , Rl,0 and Ll,p , Ll,0 =
sequence resistances and inductances of faulty
line length, RDT and LDT = resistance and
inductance of distribution transformer, Cp and
C0 = sequences capacitances of the network, RL
and LL = resistance and inductance of compos-
ite loads and Rf = fault resistance.
67
Fault Distance Estimation Using Transient Signals
Figure 5.2. The basic idea of fault location algorithm using general model of symmetrical
components.
S 4 L1 L2 C1 C2 + S 3 (L1 R2 C1 C2 + L2 R1 C1 C2 )+
S 2 (L1 C1 + L2 C2 + L1 C2 + R1 R2 C1 C2 )+
S(R1 C1 + R2 C2 + R1 C2 ) + 1 = 0 (5.1)
68
Fault Distance Estimation Using Transient Signals
R2 L2 C2
L1
E i2 i1 C1
R1
Figure 5.3. The equivalent circuit for earth fault transient of general model 1 (GM1)
algorithm. [40].
By substituting the above Equations into Equation (5.2) and Equation (5.3),
S1,2 and S3,4 now can be written as:
2 − (α )2
S1,2 = −αc ± j ωnc c (5.6)
2
S3,4 = −αd ± j ωnd − (αd )2 (5.7)
where ωnc and ωnd are undamped natural frequency of charge and dis-
69
Fault Distance Estimation Using Transient Signals
charge transient. For damping factor ζ < 1, Equation (5.6) and Equa-
tion (5.7), can now be written as:
2 − (α )2
S1,2 = −ζc ωnc ± j ωnc c (5.8)
2
S3,4 = −ζd ωnd ± j ωnd − (αd )2 (5.9)
where αc = ζc ωnc and αd = ζd ωnd . And, the damping factors for charge
and discharge transient are given as:
1
ζd = (5.10)
ωd
1+
αd
1
ζc = (5.11)
ωc
1+
αc
−αd
ωnd = (5.12)
ζd
−αc
ωnc = (5.13)
ζc
ωnd
fnd = (5.14)
2π
ωnc
fnc = (5.15)
2π
70
Fault Distance Estimation Using Transient Signals
Let us consider the LC circuit as shown in Figure 5.4. The equivalent cir-
cuit is similar as in Figure 5.3 but without the effect of damping. From the
circuit, the characteristic equation can be obtained by solving the equa-
tions for voltages around the loops of the circuit. The derivation of the
characteristic equation is given in Appendix A Section A.2. The charac-
teristic equation of circuit in Figure 5.4 is given as follows:
S 4 L1 C1 L2 C2 + S 2 (L1 C1 + L2 C2 + L1 C2 ) + 1 = 0 (5.16)
L2 C2
E i2 i1 C1
L1
Figure 5.4. The equivalent circuit for earth fault transient of general model 2 (GM2)
algorithm. [40].
where ωnc and ωnd are undamped natural frequency of charge and dis-
charge transient. Equation (5.17) and Equation (5.18) can further be ex-
pressed as in Equation (5.15) and Equation (5.14), respectively.
71
Fault Distance Estimation Using Transient Signals
In this section, we consider the equivalent circuit with the effect of load
impedances on the earth fault transients as shown in Figure 5.5. The com-
posite load impedance is the combination of impedances of the low voltage
network and impedances of the loads. The combination is connected as a
series connection. In a low voltage network, the most significant compo-
nent is the reactance of the distribution transformer [40]. However, in this
model we consider the resistances of the distribution transformer as well.
We may also consider the impedances of the low voltage lines which can
be connected in series to the impedance of the distribution transformer. In
transient analysis, a single load, with reasonable accuracy, can be mod-
eled as resistance and inductance [40]. The load is modeled as a connec-
tion of these components. In isolated and compensated neutral networks,
there is no path for zero sequence current through the transformer, there-
fore the load impedance only affects the positive and negative sequence
networks. In Figure 5.1 (g), the impedance of the composite load is cou-
pled in parallel to the network positive sequence capacitances.
R2 L2 C2 R3 L3
L1
E i2 i1 C1 i3 L4 i4 R4
R1
Figure 5.5. The equivalent circuit for earth fault transient of general model 3 (GM3)
algorithm.
72
Fault Distance Estimation Using Transient Signals
S 6 (L3 L1 L2 C2 L4 C1 )+
S 5 (C1 L4 L3 L1 R2 C2 + C1 L3 R4 L1 L2 C2 + C1 L4 L3 R1 L2 C2
+R4 L4 L1 C1 L2 C2 + C1 L4 R3 L1 L2 C2 )+
4
S (C1 L4 R3 R1 L2 C2 + L4 L1 L2 C2 + C1 L4 R3 L1 R2 C2
+C1 L4 L3 R1 R2 C2 + R4 L4 R1 C1 L2 C2 + R4 L4 L1 C1 R2 C2
+C1 L3 R4 L1 R2 C2 + L4 L3 L2 C2 + C1 R3 R4 L1 L2 C2
+C1 L3 R4 R1 L2 C2 + L3 L4 L1 C2 + L1 L3 C1 L4 )+
S 3 (R4 L4 R1 C1 R2 C2 + C1 L3 R4 R1 R2 C2 + C1 R3 R4 R1 L2 C2
+C1 L4 R3 R1 R2 C2 + C1 R3 R4 L1 R2 C2 + L4 R3 L2 C2
+L4 R4 L2 C2 + L4 L1 R2 C2 + L4 L3 R2 C2 + R4 L1 L2 C2
+L3 R4 L2 C2 + L4 R1 L2 C2 + L1 L3 C1 R4 + L1 L4 C1 R4
+L1 R3 C1 L4 + L3 R4 L1 C2 + L4 R4 L1 C2 + R3 L4 L1 C2
+R1 L3 C1 L4 + L3 L4 R1 C2 )+
2
S (L4 R4 R2 C2 + L4 R3 R2 C2 + L4 R1 R2 C2
+R4 L1 R2 C2 + R3 R4 L2 C2 + L3 R4 R2 C2
+L1 L4 + R4 R1 L2 C2 + C1 R3 R4 R1 R2 C2 +
L 3 R 4 R 1 C 2 + R3 L 4 R 1 C 2 + R3 R 4 L 1 C 2 + L 1 R 3 C 1 R4
+L4 R4 R1 C2 + R1 L3 C1 R4 + R1 L4 C1 R4 + L3 L4
+R1 R3 C1 L4 ) + S(R3 R4 R2 C2 + R4 R1 R2 C2
+L4 R4 + L3 R4 + R3 L4 + L1 R4 + R3 R4 R1 C2 + R1 R3 C1 R4
+L4 R1 ) + R1 R4 + R3 R4 = 0 (5.19)
λ1 = −a1 (5.22)
λ2 = −a2 (5.23)
73
Fault Distance Estimation Using Transient Signals
L2 C2
E i2 L1 i1 C1 i3 L5
Figure 5.6. The equivalent circuit for earth fault transient of general model 4 (GM4)
algorithm.
S 4 C1 C2 L1 L2 L5 + S 2 (C1 L1 L5 + C2 L1 L2 + C2 L1 L5 + C2 L2 L5 )
+L1 + L5 = 0 (5.24)
By substituting all the required parameters as given in Table 5.1 and in-
cluding L5 into Equation (5.24) and simplifying the equation, the roots of
the equation have the expression similar as in Equation (5.17) and Equa-
tion (5.18).
74
Fault Distance Estimation Using Transient Signals
Figure 5.7. Flowchart for fault location algorithm using general model of symmetrical
components.
75
Fault Distance Estimation Using Transient Signals
In Figure 5.7, let us assume that an earth fault at one of the feeders has
been detected. First, the transient signal is recorded and sampled. Then,
signal processing is made to the recorded transient signal by using CWT
analysis as described in Chapter 4. From CWT, the wavelet coefficient
belonging to identified charge transient component is extracted. In the
next step, the extracted CWT coefficient is used to estimate the damping
factor (ζm ) and damped transient frequency (fdm ). Both estimated pa-
rameters are used to estimate the undamped transient frequency (fnm ).
Basically, the damped transient frequency (fdm ) can be used by general
model 1 (GM1) and general model 3 (GM3) algorithm provided that the
earth fault that happens has zero fault resistance. However, in a real
case, there is always some damping due to the fault resistance. There-
fore, in order to cancel the effect of fault resistance, we used the estimated
undamped transient frequency (fnm ) to be compared with the calculated
frequency (fnc ).
Assume that, general model 1 (GM1) algorithm is selected to estimate
the fault distance. The algorithm use simple iterative process to estimate
the fault distance. A low value of calculated fault distance (lf ) is used as
an initial condition of the iteration process. From the formulated equation
of general model 1 (GM1), we calculate the undamped natural frequencies
(fnc ) and compare them with the estimated undamped charge transient
frequency (fnm ) that was recorded in the previous step. If the estimated
frequency is lower than the calculated frequency the iteration process will
be repeated by calculating the next calculated frequency (fnc ) with the
increased fault distance by increment value (∆lf ). Finally, the estimated
fault distance (lf ) is successfully achieved if the estimated frequency (fnm )
is same or higher than the calculated frequency (fnc ). In order to reduce
the number of iterations required to reach the final solution, the fault
distance increment value (∆lf ) can be increased to accelerate the iteration
process.
76
Fault Distance Estimation Using Transient Signals
Early investigation done by [130, 132, 133] shows that the transient fre-
quency can be correlated to the fault distance by using the simple RLC-
circuit. However due to the reduction and simplification of the circuit the
accuracy of the methods is reduced. The accuracy can be improved if the
network can be modeled completely. Similar to the GM algorithm, the
main principle of exact model (EM) algorithm that we propose in this sec-
tion is based on estimating the fault distance using the single frequency of
the charging process during an earth fault. However, in exact model (EM)
algorithm, we take into account the effect of inductances of sound lines.
In addition, instead of lumped the line capacitances at the substation, we
model them as an exact π-model solution. In this solution, the lines are
modeled using exact π-model equivalent circuits as shown in Figure 5.8.
In Figure 5.8, the exact π-model based on sequence component equivalent
circuits represents the positive and zero sequence components of earth
fault transients where the distribution transformers, low voltage loads
and LV network have been neglected. In the model, we assume that the
negative sequence network and positive sequence network is identical.
(a)
Ll,0’ L1,0 Lj,0
F
j=1...n
(b)
Ll,p’ Ll,p Lj,p
F
j=1...n
Figure 5.8. Exact π-model equivalent circuit with damping included. Ll,0 , Ll,p are zero
sequence and positive sequence inductance of faulty line located in front of
the fault point (F ). Ll,0′ , Ll,p′ are zero sequence and positive sequence induc-
tance of faulty line located behind of the fault point (F ). Lj,0 , Lj,p are zero se-
quence and positive sequence inductance of sound line. LT is the inductance
of the transformer. n is the number of lines. Capacitances correspondingly.
77
Fault Distance Estimation Using Transient Signals
network in complex frequency domain. For single phase to earth fault, the
diagonal elements of the nodal transform impedance matrix in complex
frequency (s) is given as:
1 0 2 1
Zkk (s) = Zkk (s) + Zkk (s) (5.25)
3 3
1 2
Zkk (s) = Zkk (s) (5.26)
The fault resistance (Rf ) was not required to be added to the equation.
This is because, in this work we purposely formed the impedance ma-
trix for undamped case of ground fault transient. Then, in order to esti-
mate the fault distance, we substitute the estimated undamped complex
frequency (fnm ) of the charge transient for complex frequency (s) in Equa-
tion (5.25). Then, we solve the roots of a polynomial equation of Zkk (s) = 0.
In the computer model, the lowest positive value of the roots from the
simulation results gives the estimated fault distance. The explanation of
exact the π-model equivalent circuit that was used in exact model (EM)
algorithm is the main subject of the next subsection. The flowchart for
estimating the earth fault distance based on exact model (EM) algorithm
is given at the end of this section.
78
Fault Distance Estimation Using Transient Signals
the information of network components and real time data of their con-
nections, the equivalent circuit exact π-model can be used in computer
simulation model using nodal analysis. As an example, Equation (5.27)
shows the sequence admittance matrix in complex frequency (s) of Fig-
ure 5.9 that has k number of buses.
(a)
D*Ll,0’ F k*L1,0 Lj,0
N3 N2 N1 Nk
D*Cl,0’ k*Cl,0
2 2 Cj,0 Cj,0
D*Cl,0’ k*Cl,0 2 2
2 2
j=1...n
(b)
D*Ll,p’ F k*L1,p N1 Lj,p
N3 Nk
N2
D*Cl,p’ k*Cl,p LT
2 2 Cj,p Cj,p
D*Cl,p’ k*Cl,p 2 2
2 2
j=1...n
Figure 5.9. The equivalent circuit for earth fault transient of exact model 1 (EM1) algo-
rithm. Ll,0 , Ll,p are zero sequence and positive sequence inductance of faulty
line located in front of fault point (F ). Ll,0′ , Ll,p′ are zero sequence and pos-
itive sequence inductance of faulty line located behind of the fault point (F ).
Lj,0 , Lj,p are zero sequence and positive sequence inductance of sound line.
LT is the inductance of the transformer, k is a fault distance, T is the total
lenght of faulty feeder, D = T −k is the total length of faulty feeder (T ) minus
fault distance (k), N1 to Nk is the number of the buses and n is the number
of lines. Capacitances correspondingly.
⎡ ⎤
Y11 (s) Y12 (s) Y13 (s) · · · Y1k (s)
⎢ ⎥
⎢Y21 (s) Y22 (s) Y23 (s) · · · Y2k (s)⎥
⎢ ⎥
⎢ ⎥
Y012 (s) = ⎢Y31 (s) Y32 (s) Y33 (s) · · · Y3k (s)⎥ (5.27)
⎢ ⎥
⎢ .. .. .. .. .. ⎥
⎢ ⎥
⎢ . . . . . ⎥
⎣ ⎦
Yk1 (s) Yk2 (s) Yk3 (s) · · · Ykk (s)
79
Fault Distance Estimation Using Transient Signals
(a)
F k*L1,0 Lj,0
N2 Nk
N1
k*Cl,0 Cj,0
2 Cj,0
D*Cl,0
2
j=1...n
(b)
F k*L1,p Lj,p
N1 Nk
N2
k*Cl,p LT Cj,p
D*Cl,p 2 Cj,p
2
j=1...n
Figure 5.10. The equivalent circuit for earth fault transient of exact model 2 (EM2) al-
gorithm. Ll,0 , Ll,p are zero sequence and positive sequence inductance of
faulty line located in front of the fault point (F ). Lj,0 , Lj,p are zero sequence
and positive sequence inductance of sound line. LT is the inductance of
the transformer, k is a fault distance, T is the total lenght of faulty feeder,
D = T − k is the total length of faulty feeder (T ) minus fault distance (k), N1
to Nk is the number of the buses and n is the number of lines. Capacitances
correspondingly.
80
Fault Distance Estimation Using Transient Signals
We propose the exact model 3 (EM3) algorithm with the use of exact
π-model as shown in Figure 5.11. The model is constructed similar to
exact π-model in EM2 algorithm, except that the sequence capacitance
(Cl,0 ,Cl,p ) of its faulty feeder section is located at bus N1 as depicted in
Figure 5.11. Through the developed computer model, we simulate and
solve the impedance matrix in complex frequency (s) as given in Equa-
tion (5.25).
F (a)
k*L1,0 Lj,0
N1 Nk
j=1...n
F (b)
k*L1,p Lj,p
Nk
N1
LT
T*Cl,p Cj,p Cj,p
2 2
j=1...n
Figure 5.11. The equivalent circuit for earth fault transient of exact model 3 (EM3) al-
gorithm. Ll,0 , Ll,p are zero sequence and positive sequence inductance of
faulty line located in front of the fault point (F ). Lj,0 , Lj,p are zero sequence
and positive sequence inductance of sound line. LT is the inductance of
the transformer, k is a fault distance, T is the total lenght of faulty feeder,
D = T − k is the total length of faulty feeder (T ) minus fault distance (k), N1
to Nk is the number of the buses and n is the number of lines. Capacitances
correspondingly.
In Figure 5.12, the flowchart shows the steps used by exact model algo-
rithm for computing the undamped transient frequency (fnm ) is similar as
explained in previous section. Let us assume that Exact model 1 (EM1)
algorithm is selected to estimate the fault distance. Then, using a esti-
mated damped transient frequency (fdm ), we calculate the inductance per
unit length of line network using Frequency-Dependent(FD) Line Model
or known as JMarti FD model in EMTP environment. In EMTP envi-
ronment, we can calculate the sequence impedance of the line at specific
estimated charge transient frequency using line check calculation tools
provided that we have the details of line geometrical and electrical val-
ues. Next, the estimated undamped frequency (fnm ) and calculated line
81
Fault Distance Estimation Using Transient Signals
Figure 5.12. Flowchart for fault location algorithm using exact π-model of symmetrical
component.
82
Fault Distance Estimation Using Transient Signals
Lf
lf = 1 ′ ′ ′ (5.30)
3 (L0 + L1 + L2 )
′ ′ ′
where, L0 , L1 , L2 are the inductances per unit length of the zero, positive
and negative-sequence systems respectively. The final value of the fault
distance (lf ) was selected based on the smallest standard deviation of the
mean value of the fault distance calculated in Equation (5.30) for a set
of subsequent samples. The flow chart representing the CWT algorithm
developed in this work is shown in Figure 5.13.
83
Fault Distance Estimation Using Transient Signals
84
Fault Distance Estimation Using Transient Signals
85
Fault Distance Estimation Using Transient Signals
Figure 5.14. Flowchart for MRA and NN fault location algorithms based on recorded
transient signal of single MV/LV distribution transformer.
86
Fault Distance Estimation Using Transient Signals
When an earth fault happens in MV feeder, the transient signal can also
be detected from the secondary side of the MV/LV distribution trans-
former. The transient component can easily be distinguished from the
fundamental frequency component. Figure 5.15 shows an example of sim-
ulated transient signal recorded at the MV and LV side of the network.
In Figure 5.14, let us assume that an earth fault has been detected in
MV network. From the assumption that an earth fault has been detected,
the transient signal measured from MV/LV distribution transformers is
recorded and sampled. Then, the recorded signal will be converted from
phase to line voltage. Later, the signal is analyzed with CWT and the
CWT extraction features process is applied to the transient signal as has
been explained in Chapter 4, Section 4.4. Figure 5.16 shows an exam-
ple of extracted CWT coefficient of transient signal recorded from MV/LV
distribution transformer.
4
3 x 10
Va
Vb
2 Vc
1
Voltage [V]
−1
−2
−3
Fault Inception
−4
0 100 200 300 400 500 600 700 800 900 1000
Sample Number[k]
800
Vab
Fault Inception Vbc
600 Vca
400
Voltage [V]
200
−200
−400
−600
0 100 200 300 400 500 600 700 800 900 1000
Sample Number[k]
Figure 5.15. An example of simulated transient signal recorded at (a) the MV and (b) LV
side of the MV\LV distribution transformer
87
Fault Distance Estimation Using Transient Signals
It has been reported that the loads have a significant effect to the tran-
sient signal [40,139]. Therefore, due to this reason, the total load (MVA) of
the network prior an earth fault happens in MV network is also recorded.
Based on the extracted CWT coefficient characteristics, the peak value of
its FFT amplitude and its coefficient scale is recorded. Actually, the scale
of CWT coefficient is correlative to its frequency. Therefore, either one
of the parameters can be used. After computing and recording the two
characteristics of extracted CWT coefficients, with these two parameters
including the total load (MVA) of the network prior a to the earth fault
they are used as the inputs database of the MRA and ANN algorithms.
Multiple regression analysis (MRA), artificial neural network (ANN) and
correct path algorithms are the main subject of the following subsection.
500
V
250
Voltage[V]
−250
−500
100 150 200 250 300 350 400 450 500
Sample Number[k]
10
0
0 5 10 15 20 25 30 35 40
Scales[a]
100
FFT
90
80 X: 996.1
Y: 71.21
70
Amplitude[FFT]
60
50
40
30
20
10
0
0 500 1000 1500 2000 2500 3000
Frequency[Hz]
88
Fault Distance Estimation Using Transient Signals
89
Fault Distance Estimation Using Transient Signals
l = β0 + β1 A + β2 S + β3 L (5.33)
where ANN is an implicit non mathematical function that can only can be
designated as ’knowledge base’. Basically, using the same dependent and
independent variables as developed for MRA formula can also be used for
artificial neural network analysis. A proposed designed neural network
90
Fault Distance Estimation Using Transient Signals
is constructed with a two layers, feed forward network with tangent sig-
moid hidden neurons and linear output neuron. The developed neural
network training function updates weight and bias values according to
Levenberg-Marquardt optimization algorithm. ANN architecture of the
proposed ANN algorithm is depicted in Figure 5.17.
Output value
Output layer
(1 neuron)
Hidden layer
(20 neurons)
A S L Input Layer
(3 neurons)
Inputs
Figure 5.17. Architecture of neural network for estimate the fault distance.
The designed neural network has 3 inputs variables which contain valu-
able and relevant information for estimating fault distance. Each input
unit of neural network has a symbol associated with it where A is repre-
sented as the peak value of FFT amplitude of extracted CWT coefficient,
S is associated its coefficient scale and L is referred to total loads (MVA)
of the network prior to an earth fault. For fault distance estimation al-
gorithm in this work the ANN consists only one hidden layer. For this
kind of input-output task, one single hidden layer is assumed to be good
enough [79]. The outputs of hidden layer has sigmoid activation function
and are transferred to the output layer which is composed of only one
neuron which gives the estimated fault distance. It is important to train
and test the neural network with a large data set. The training process is
done in several times and each time the performance after each iteration
is saved where the performance is defined as the correlation coefficient
(r). Finally the saved trained network that has the highest value of cor-
relation coefficient (r) is used to for fault location algorithm. Since this
method is based on the ’knowledge base’ of the specific network, therefore
the NN have to retrain if there are any changes in terms of the network
topology.
The proposed algorithm have has tested using simulation model based
on a practical 20 kV unearthed network distribution network as shown
in Chapter 6 Section 6.4. The data set used in training and testing the
91
Fault Distance Estimation Using Transient Signals
neural network is given in Chapter 6, Table 6.3 and Table 6.4, respec-
tively. The result and performance of the algorithm will be discussed in
that chapter.
Figure 5.18. Flowchart for correct path location algorithm based on recorded transients
of MV/LV distribution transformers.
The line to line fault voltage transient signal measured from the sec-
92
Fault Distance Estimation Using Transient Signals
93
Fault Distance Estimation Using Transient Signals
94
6. Performance Evaluation of Proposed
Fault Location Algorithms Using
EMTP/ATP Simulation Model
In this chapter, the proposed algorithms are tested using transient signal
generated by a simulation model. The simulation of the voltage and cur-
rent transient is required and important in order to validate the proposed
fault location algorithms. In this research, we used Alternative Transient
Program (ATP/EMTP) as transient simulator. The program was chosen
because it is one of the most common simulation tools used by power en-
gineers and researcher for simulating transients [142]. The program soft-
ware contains extensive modeling capabilities for transmission lines, ca-
bles, breakers, loads, converters, protection devices, non-linear elements,
electromagnetic coupling, and major power electronics devices and equip-
ment. It also has an enhanced graphical user interface called ATPDraw
as a pre processor, which allows an easy entry of system topology and
data [143] . The program software will in future be referred to only as
EMTP.
In this chapter, four simulation models are presented. The first model
is according to general model of symmetrical component as we have dis-
cussed in Chapter 5, Section 5.1. In simulation model 2, we modeled the
unearthed neutral network with exact π-model as described in Chapter 5,
Section 5.2. Then in the following section, we investigate the earth fault
network using the line constants EMTP program and taking into account
their real geometrical and electrical values. In the last section of this
chapter, we tested the fault location algorithms with the modified system
from [144] and modeled the network using line constant EMTP program.
The earth faults are simulated with different scenarios such as fault dis-
tance, fault resistance, fault inception angle and also the effect of load.
Some results of this studies can be found in [137].
95
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
96
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
OH line
OH line
Figure 6.1. A Schematic(a) and ATPdraw simulated network model(b) for earth fault
simulation analysis in 20 kV unearthed medium votage network. The details
of the abbreviations in ATPdraw simulation model can be found in Table 5.1
n n
1 1
M AEGM = |GMi − li | = |∆GMi | (6.1)
n n
i=1 i=1
n n
1 1
M AE{f,ζ} = |∆{fm − fc , ζm − ζc }| = |∆{fi , ζi }| (6.2)
n n
i=1 i=1
97
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
where ∆|{fi , ζi }| is the absolute error of the parameters (f, ζ), fm is es-
timated transient frequency, fc is calculated transient frequency, ζm is
estimated damping factor and ζc is calculated damping factor and n is
numbers of observation of parameter error.
Figures 6.2 and 6.3 shows the example of earth fault transient voltage
and current for earth fault located at 2 km, 10 km and 18 km at faulty line
feeder with overhead line as shown in Figure 6.1. The transient signals
were captured at location "M" as shown in Figure 6.1. The earth fault was
simulated at time t=0 and the signal was captured before the time switch
is "on" at t=0.005 second.
8000
6000 V1
4000
Voltage [V]
2000 X: 100
Y: 0
0
−2000
−4000
−6000
−8000
0 100 200 300 400 500 600 700 800 900 1000
Sample Number [k]
8000
V1
6000
4000
Voltage [V]
2000 X: 100
Y: 0
0
−2000
−4000
−6000
−8000
0 100 200 300 400 500 600 700 800 900 1000
Sample Number [k]
8000
V1
6000
4000
Voltage [V]
2000 X: 100
Y: 0
0
−2000
−4000
−6000
−8000
0 100 200 300 400 500 600 700 800 900 1000
Sample Number [k]
Figure 6.2. Simulated transient voltage signal in the faulty phase when a single phase to
earth fault happens at a feeder in an example network with overhead lines.
The fault distances (lf ) are 2, 10 and 18 km from the substation. The fault
resistance is zero ohm.
98
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
100
I
80 1
60
40
Current [A]
20 X: 100
Y: −0
0
−20
−40
−60
−80
−100
0 100 200 300 400 500 600 700 800 900 1000
50
I
1
40
30
20
Current [A]
10 X: 100
Y: −0
0
−10
−20
−30
−40
0 100 200 300 400 500 600 700 800 900 1000
Sample Number [k]
40
I1
30
20
Current [A]
10
X: 100
Y: −0
0
−10
−20
−30
0 100 200 300 400 500 600 700 800 900 1000
Sample Number[k]
Figure 6.3. Simulated transient current signal in the faulty phase when a single phase to
earth fault happens at a feeder in an example network with overhead lines.
The fault distances (lf ) are 2, 10 and 18 km from the substation. The fault
resistance is zero ohm.
99
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
100
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
Voltage[V]
2000
0
−2000
Fault inception
50 100 150 200 250 300 350 400 450 500
40
37 0.12
34
31 0.1
Scales a
28
25 Charge transient component 0.08
22
19 0.06
16
13 Discharge transient component
10 0.04
7
4 0.02
1
50 100 150 200 250 300 350 400 450 500
Sample Number[k]
(a)
7
X: 23 VEP
X: 14
6 Y: 5.68 Y: 6.85
Energy percentage[%]
5
0
0 5 10 15 20 25 30 35 40 45 50
Scales[a]
(b)
4
x 10
1.5
V
CWT
1
0.5
Voltage[V]
−0.5
−1
−1.5
0 50 100 150 200 250 300 350 400 450 500
Saample Number[k]
(c)
4500
V
4000 X: 703.1
Y: 3632
3500
Magnitude[FFT]
3000
2500
2000
1500
1000
500
0
0 500 1000 1500 2000 2500 3000 3500 4000 4500 5000
Frequency (Hz)
(d)
10
V
HT
9.5
9
log(|V (t)|)
8.5
8
HT
7.5
6.5
5.5
0 50 100 150 200 250 300 350 400 450 500
Sample Number[k]
(e)
Figure 6.4. Example of the output from pre-processing of voltage transient signals using
CWT and Hilbert transformation. The fault distance is 10 km from the sub-
station. (a) Simulated transient voltage (above) and its CWT scalogram (be-
low), (b) scale-dependent energy percentage distribution, (c) extracted CWT
coefficient, (d) FFT frequency spectrum of extracted CWT coefficient, (e) lin-
ear decaying of extracted CWT coefficient.
101
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
Current[A]
20
0
−20 Fault inception
50 100 150 200 250 300 350 400 450 500
40
37
34 0.15
31
Scales a
28
25 Charge transient component
22 0.1
19
16
13
10 0.05
7
4
1
50 100 150 200 250 300 350 400 450 500
Sample Number[k]
(a)
9
I
EP
8 X: 24
Y: 8.182
Energy percentage[%]
0
0 5 10 15 20 25 30 35 40 45 50
Scales[a]
(b)
200
I
CWT
150
100
Current[A]
50
−50
−100
−150
−200
0 50 100 150 200 250 300 350 400 450 500
Sample Number[k]
(c)
80
I
70 X: 703.1
Y: 62.94
60
Magnitude[FFT]
50
40
30
20
10
0
0 500 1000 1500 2000 2500 3000 3500 4000 4500 5000
Frequency (Hz)
(d)
5.5
I
HT
5
4.5
log(|I (t)|)
3.5
HT
2.5
1.5
1
0 50 100 150 200 250 300 350 400 450 500
Sample Number[k]
(e)
Figure 6.5. Example of the output from pre-processing of current transient signals using
CWT and Hilbert transformation. The fault distance is 10 km from the sub-
station. (a) Simulated transient current (above) and its CWT scalogram (be-
low), (b) scale-dependent energy percentage distribution, (c) extracted CWT
coefficient, (d) FFT frequency spectrum of extracted CWT coefficient, (e) lin-
ear decaying of extracted CWT coefficient.
102
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
40 0.2
Current[A]
20
0.18
0
0.16
−20
0.14
−40
100 120 140 160 180 200 220 240 260 280 300
0.12
49
46
43
40 0.1
37
Scales a 34
31 0.08
28
25
22 0.06
19
16
13 0.04
10
7
4
1 0.02
100 120 140 160 180 200 220 240 260 280 300
Sample Number [k]
(a)
40 0.22
Current[A]
20 0.2
0
0.18
−20
0.16
−40
100 120 140 160 180 200 220 240 260 280 300
0.14
49
46
43 0.12
40
37
34 0.1
Scales a
31
28 0.08
25
22
19 0.06
16
13
10 0.04
7
4
1 0.02
100 120 140 160 180 200 220 240 260 280 300
Sample Number [k]
(b)
40 0.18
Current[A]
20
0.16
0
0.14
−20
100 120 140 160 180 200 220 240 260 280 300
0.12
49
46
43 0.1
40
37
34 0.08
Scales a
31
28
25
22 0.06
19
16
13 0.04
10
7
4
1 0.02
100 120 140 160 180 200 220 240 260 280 300
Sample Number [k]
(c)
Figure 6.6. Example of the output from pre-processing of current transient signals using
CWT with the effect of fault resistance. The fault distance is 10 km from the
substation. (a) Simulated transient current (above) and its CWT scalogram
(below) with Rf =25 Ω, (b) Rf =50 Ω and (c) Rf =100 Ω.
103
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
50 0.14
Current[A]
0 0.12
−50 0.1
100 120 140 160 180 200 220 240 260 280 300
49
46
43 0.08
40
37
34
Scales a
31
28 0.06
25
22
19
16 0.04
13
10
7
4
1 0.02
100 120 140 160 180 200 220 240 260 280 300
Sample Number [k]
(a)
50 0.12
Current[A]
0.11
0
0.1
0.09
−50
100 120 140 160 180 200 220 240 260 280 300 0.08
49
46
43 0.07
40
37
34 0.06
Scales a
31
28 Scale = 25
0.05
25
22
19 0.04
16
13
10 0.03
7
4
1 0.02
100 120 140 160 180 200 220 240 260 280 300
Sample Number [k]
(b)
30 0.11
Current[A]
20 0.1
10
0.09
0
0.08
−10
100 120 140 160 180 200 220 240 260 280 300 0.07
49
46
43 0.06
40
37
34 0.05
Scales a
31 Scale = 28
28 0.04
25
22
19 0.03
16
13
10 0.02
7
4
1 0.01
100 120 140 160 180 200 220 240 260 280 300
Sample Number [k]
(c)
Figure 6.7. Example of the output from pre-processing of current transient signals using
CWT with the effect of fault resistance. The fault distance is 10 km from the
substation. (a) Simulated transient current (above) and its CWT scalogram
(below) with Rf =200 Ω, (b) Rf =300 Ω and (c) Rf =400 Ω.
Figures 6.8 and 6.9 shows the scale-energy percentage distribution for
the simulated transient signal. Using the method for identification of
charge transient signal as explained in Chapter 4 Section 4.4, in this par-
ticular case the charge transient component can only be detected for earth
fault with fault resistance up to 200 Ω. However, if we observed the result
104
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
of CWT scalogram in Figure 6.7(b) and (c), we can make a guess estima-
tion what is the scale corresponding to the charge transient frequency by
inspecting the maxima of the energy. The maxima of the energy in Fig-
ure 6.7(b) is located between scale 22 and 28, while in Figure 6.7(c), it is
located between scale 22 and 37. The arrows in Figure 6.7(b) and (c) show
the area of the maxima of the energy. The purple line in Figure 6.7(b) and
(c), shows the guess estimation of the scales that correspond to the charge
transient frequency of the measured transient signal.
The calculated and estimated frequencies and damping factor of the
charge transient signal as a function of fault distance are shown in Fig-
ure 6.10, and the tabulated results are found in Appendix C Table C1.
Based on the result shown in Figure 6.10, the calculated and the esti-
mated values have a good agreement. The range of estimated damped
frequency is between 1132 Hz to 507 Hz. In Figure 6.10(a) and (b), the
graph shows that the damped and undamped charge transient frequency
is decreasing when the fault distance is increasing. The damping factor is
increased when the fault distance is increasing. The mean absolute error
(MAE) index of damped and undamped transient frequency is about 10
Hz, and 0.0006 for damping factor, as shown in Appendix C Table C1.
The results of fault distance estimation error using GM1 and GM2 algo-
rithms for simulated earth fault with different fault distances is given in
Figure 6.11, and the tabulated results are found in Appendix C Table C2.
As shown in Figure 6.11, basically the results from tested algorithms are
almost identical. In overall the distance calculation error is less than 1.5
km. The MAE index calculated for all algorithms is 0.36 km, as given in
Appendix C Table C2. The fault distance error is slightly higher when the
fault happens at the middle and at the end of the faulty feeder. From the
results, we found that the lowest fault distance estimation error is when
fault is simulated at 18 km from the substation. The reason why it has the
lowest fault distance estimation error is because based on the comparison
between the estimated (fdm ) and the calculated (fnm ) undamped charge
transient frequency, it has the lowest error (∆fn ) as given in Appendix C
Table C1.
105
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
4.5
4 X: 25
Y: 4.297
Energy percentage[%]
3.5
2.5
1.5
0.5
0
0 5 10 15 20 25 30 35 40 45 50
Scales[a]
(a)
4.5
4
X: 25
Y: 4.054
3.5
Energy percentage[%]
2.5
1.5
0.5
0
0 5 10 15 20 25 30 35 40 45 50
Scales[a]
(b)
4
3.5 X: 26
Y: 3.675
Energy Percentage[%]
2.5
1.5
0.5
0
0 5 10 15 20 25 30 35 40 45 50
Scales[a]
(c)
Figure 6.8. Example of scale-dependent energy percentage distribution. The fault dis-
tance is 10 km from the substation. (a) Scale-dependent energy percentage
distribution with Rf =25 Ω, (b) Rf =50 Ω and (c) Rf =100 Ω.
106
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
3.5
X: 31
3 Y: 3.253
Energy Percentage[%]
2.5
1.5
0.5
0
0 5 10 15 20 25 30 35 40 45 50
Scales[a]
(a)
3.5
3
Energy Percentage[%]
2.5
1.5
0.5
0
0 5 10 15 20 25 30 35 40 45 50
Scales[a]
(b)
4
3.5
Energy Percentage[%]
2.5
1.5
0.5
0
0 5 10 15 20 25 30 35 40 45 50
Scales[a]
(c)
Figure 6.9. Example of scale-dependent energy percentage distribution. The fault dis-
tance is 10 km from the substation. (a) Scale-dependent energy percentage
distribution with Rf =200 Ω, (b) Rf =300 Ω and (c) Rf =400 Ω.
107
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
1200
fdc
fdm
1100
Frequency [Hz]
1000
900
800
700
600
500
2 4 6 8 10 12 14 16 18 20
Fault Distance [km]
1200
fnc
fnm
1100
1000
Frequency [Hz]
900
800
700
600
500
2 4 6 8 10 12 14 16 18 20
Fault Distance [km]
0.05
ζc
ζm
0.04
0.03
0.02
0.01
0
2 4 6 8 10 12 14 16 18 20
Fault Distance [km]
Figure 6.10. The comparison of simulation test result to the theory of earth fault tran-
sient. fdc , fnc and ζc are calculated using the GM1 equations. fdm , fnm
and ζm are estimated transient parameters from the simulated earth fault
transient waveforms. fdm = estimated damped transient frequency (Hz), fdc
= calculated damped transient frequency (Hz), ζm = estimated damping fac-
tor, ζc = calculated damping factor, fn = undamped transient frequency (Hz),
fnm = estimated undamped transient frequency (Hz) and fnc = calculated
undamped transient frequency (Hz)
108
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
2
GM1
fdm
1.5 GM1
fnm
GM2
1
Error[km]
0.5
−0.5
−1
2 4 6 8 10 12 14 16 18 20
Fault Distance[km]
Figure 6.11. Comparison of the distance estimates obtained using general model (GM)
algorithms for the simulated single line to earth faults with varies fault dis-
tance. GM1f dm denotes a distance calculation error using general model 1
(GM1) with damped transient frequency (fdm ), GM1f nm denotes a distance
calculation error using general model 1 (GM1) with undamped transient
frequency (fnm ) and GM2 denotes a distance calculation error using gen-
eral model 2 (GM2) algorithm.
Figures 6.12, 6.13, 6.14, 6.15 and 6.16 show the results of simulated
single line to earth fault transients with the effect of fault resistance,
and the tabulated results are found in Appendix C Table C3. In gen-
eral, referring to general model (GM) equivalent circuit, the transient
frequencies become lower when the resistance (R2 ) is increased. In Fig-
ure 6.12, the graphs show that the calculated and estimated frequencies
have small reduction when fault resistance is increased. The reduction
can be seen more clearly in estimated frequency when high fault resis-
tance (Rf = 100Ω) is used. In Figure 6.13, the estimated and calculated
frequency is almost identical except when fault resistance is 100 Ω.
Based on the results in Figure 6.13 (a), (b) and (c), for fault close to
the substation and for the certain range of the fault resistances, the gap
between estimated and calculated frequencies can be seen more clearly.
Based on the MAE index shown in Appendix C Table C3, the mean abso-
lute error (MAE) index of damped and undamped transient frequency is
about 22 Hz. The value is slightly higher compared to the results without
the effect of fault resistance. The reason is due to the effect of fault resis-
tance to the transient signal. The fault resistance increases the damping
of the estimated signal. The MAE of damping factor is 0.01, this is about
17 times higher than MAE of the damping factor without the effect of fault
resistance. When damping is higher the transient almost disappears and
subsequently the identification of the charge transient component will be
difficult.
109
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
1200
fdc−0Ω
1100 fdc−25Ω
fdc−50Ω
fdc−100Ω
Frequency [Hz]
1000
900
800
700
600
500
2 4 6 8 10 12 14 16
Fault Distance [km]
1200
fdm−0Ω
1100 fdm−25Ω
fdm−50Ω
fdm−100Ω
Frequency [Hz]
1000
900
800
700
600
500
2 4 6 8 10 12 14 16
Fault Distance [km]
1200
fnc−0Ω
fnc−50Ω
1100 fnc−25Ω
fnc−100Ω
Frequency [Hz]
1000
900
800
700
600
500
2 4 6 8 10 12 14 16
Fault Distance [km]
1100 fnm−0Ω
fnm−25Ω
fnm−50Ω
1000 fnm−100Ω
Frequency [Hz]
900
800
700
600
2 4 6 8 10 12 14 16
Fault Distance [km]
Figure 6.12. Result of calculated and estimated damped (fdc , fdm ) and undamped
(fnc , fnm ) frequencies with the effect of fault resistance (Rf ). Fault resis-
tance (Rf )= 0, 25, 50, 100 Ω. For details of frequencies abbreviations see
Figure 6.10.
110
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
1200
fdc−0Ω
1100 fdm−0Ω
Frequency [Hz]
1000
900
800
700
600
500
2 4 6 8 10 12 14 16
Fault Distance [km]
1200
fdc−25Ω
1100 fdm−25Ω
Frequency [Hz]
1000
900
800
700
600
500
2 4 6 8 10 12 14 16
Fault Distance [km]
1200
fdc−50Ω
1100 fdm−50Ω
Frequency [Hz]
1000
900
800
700
600
500
2 4 6 8 10 12 14 16
Fault Distance [km]
1200
fdc−100Ω
1100 fdm−100Ω
Frequency [Hz]
1000
900
800
700
600
500
2 4 6 8 10 12 14 16
Fault Distance [km]
Figure 6.13. Comparison of result of calculated (fdc ) and estimated (fdm ) damped fre-
quencies with the effect of fault resistance (Rf ). Fault resistance (Rf )= 0,
25, 50, 100 Ω. For details of frequencies abbreviations see Figure 6.10.
111
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
fnc−0Ω
1100 fnm−0Ω
Frequency [Hz]
1000
900
800
700
600
2 4 6 8 10 12 14 16
Fault Distance [km]
fnc−25Ω
1100 fnm−25Ω
1000
Frequency [Hz]
900
800
700
600
2 4 6 8 10 12 14 16
Fault Distance [km]
fnc−50Ω
1100 fnm−50Ω
Frequency [Hz]
1000
900
800
700
600
2 4 6 8 10 12 14 16
Fault Distance [km]
fnc−100Ω
1100 fnm−100Ω
Frequency [Hz]
1000
900
800
700
600
2 4 6 8 10 12 14 16
Fault Distance [km]
Figure 6.14. Comparison of result of calculated (fnc ) and estimated (fnm ) undamped fre-
quencies with the effect of fault resistance (Rf ). Fault resistance (Rf )= 0,
25, 50, 100 Ω. For details of frequencies abbreviations see Figure 6.10.
112
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
0.3
ζc−0Ω
ζc−25Ω
0.25 ζc−50Ω
ζc−100Ω
0.2
0.15
0.1
0.05
0
2 4 6 8 10 12 14 16
Fault Distance [km]
0.3
ζm−0Ω
ζm−25Ω
0.25 ζm−50Ω
ζm−100Ω
0.2
0.15
0.1
0.05
0
2 4 6 8 10 12 14 16
Fault Distance [km]
Figure 6.15. Result of calculated (ζc ) and estimated (ζm ) damping factor with the effect
of fault resistance (Rf ). Fault resistance (Rf )= 0, 25, 50, 100 Ω. For details
of frequencies abbreviations see Figure 6.10.
113
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
0.06
ζc−0Ω
ζm−0Ω
0.05
0.04
0.03
0.02
0.01
0
2 4 6 8 10 12 14 16
Fault Distance
0.08
ζc−25Ω
ζm−25Ω
0.075
0.07
0.065
0.06
0.055
0.05
2 4 6 8 10 12 14 16
Fault Distance [km]
0.12
ζc−50Ω
ζm−50Ω
0.115
0.11
0.105
0.1
0.095
0.09
0.085
0.08
2 4 6 8 10 12 14 16
Fault Distance [km]
0.4
ζc−100Ω
ζm−100Ω
0.35
0.3
0.25
0.2
0.15
0.1
0.05
0
2 4 6 8 10 12 14 16
Fault Distance [km]
Figure 6.16. Comparison of result of calculated (ζc ) and estimated (ζm ) damping factor
with the effect of fault resistance (Rf ). Fault resistance (Rf )= 0, 25, 50, 100
Ω. For details of frequencies abbreviations see Figure 6.10.
114
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
−0.7
GM1fdm
0.7
GM1fnm
Error[km]
0.6
GM2
0.3
0.1
−0.1
−0.3
−0.5
−0.7
2 4 10 16
Fault Distance[km]
0.4
GM2
0.2
0
−0.2
−0.4
−0.6
−0.8
2 4 10 16
Fault Distance[km]
4
GM1
3.5 fdm
GM1fnm
3
GM2
Error[km]
2.5
2
1.5
1
0.5
0
2 4 10 16
Fault Distance[km]
4
GM1
3.5 fdm
GM1fnm
3
GM2
Error[km]
2.5
2
1.5
1
0.5
0
2 4 10 16
Fault Distance[km]
Figure 6.17. Distance estimates with effect of fault resistance for general model net-
works.
115
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
Figure 6.17, shows the results of the estimated of fault location using
GM1 and GM2 algorithms with the effect of fault resistance, Appendix
C Table C4 shows the tabulated simulation results. Referring to Fig-
ure 6.17, we can see that result of the distance estimation error is almost
identical for the tested algorithms. In overall, the calculated MAE index
of all algorithms is less than 0.6 km, as given in Appendix C Table C4.
The fault distance error is lower when the fault is simulated at the mid-
dle of the faulty feeder with the fault resistance between 0 to 50 Ω. As
expected, the fault estimation error is higher when the fault is located at
the end of the line. In overall, the fault distance error is less than 2 km.
Appendix C, Table C5 shows the results of simulated single line to earth
fault with the effect of fault inception angle. The fault inception angles
used are 0, 45 and 90 ◦ . In overall, we found that the fault inception angle
does not affect the transient frequencies. Unfortunately, we found that
the process to locate the transient component becomes much challenging
since the lower fault inception angle will reduce the peak amplitude of
the transient, subsequently effecting the process of identification of the
charge transient component as shown in Figure 6.18.
In Figure 6.19, the results show estimates of fault location using GM1
algorithm. The fault distance estimation is calculated using a estimated
damped and undamped transient frequency. A single line to earth fault
was simulated with zero fault resistance. The MAE index of the simula-
tion result was 0.34 km.
Figure 6.20 shows results of calculated damped transient frequencies
with the effect of load and fault resistance. The frequency slightly in-
creases when the load is increased and slightly decreases when the fault
resistance is increased. The effect of load on damped transient frequency
is higher when the fault located close to the substation. Figure 6.21,
shows the estimated damped transient frequencies with the effect of load
and fault resistance. The estimated damped frequency shows significant
decrease when the fault resistance is increased.
116
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
Current[A]
0
−5
50 100 150 200 250 300
30 0.2
Scales [a]
25
20 0.15
15
0.1
10
5 0.05
(a) Transient current (above) and its scalogram, fault inception angle 0◦
Current[A]
20
0
−20
−40
50 100 150 200 250 300
30
0.2
25
Scales [a]
20 0.15
15
0.1
10
5 0.05
50 100 150 200 250 300
Sample Number[k]
(b) Transient current (above) and its scalogram, fault inception angle 45 ◦
Figure 6.18. Comparison of simulated transient signal and its scalogram with the effect
of inception angle. The fault distance is 4 km and without fault resistance.
0.6
GM1fdm
0.4 GM1fnm
Error(km)
0.2
−0.2
−0.4
−0.6
490◦ 445◦ 40 ◦ 1690◦ 1645◦ 160◦
Fault Distance(km)Inception angle
Figure 6.19. Effect of fault inception angle on fault distance estimation. GM1f dm de-
notes a distance estimation error using GM1 with damped transient fre-
quency and GM1f nm denotes a distance estimation error using GM1 with
undamped transient frequency.
117
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
1100
fdc−0Ω−A
fdc−25Ω−A
1000 fdc−50Ω−A
Frequency [Hz]
fdc−100Ω−A
900
800
700
600
500
4 6 8 10 12 14 16
Fault Distance [km]
1100
fdc−0Ω−B
fdc−25Ω−B
1000 fdc−50Ω−B
fdc−100Ω−B
Frequency [Hz]
900
800
700
600
500
4 6 8 10 12 14 16
Fault Distance [km]
1100
fdc−0Ω−A
fdc−100Ω−A
1000 fdc−0Ω−B
fdc−100Ω−B
Frequency [Hz]
900
800
700
600
500
4 6 8 10 12 14 16
Fault Distance [km]
Figure 6.20. Comparison of calculated damped frequencies with the effect of load and
fault resistance
118
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
1000 fdm−0Ω−A
fdm−25Ω−A
fdm−50Ω−A
Frequency [Hz]
900 fdm−100Ω−A
800
700
600
4 6 8 10 12 14 16
Fault Distance [km]
1100
fdm−0Ω−B
fdm−25Ω−B
1000 fdm−50Ω−B
Frequency [Hz]
fdm−100Ω−B
900
800
700
600
500
4 6 8 10 12 14 16
Fault Distance [km]
1100
fdm−0Ω−A
fdm−100Ω−A
1000 fdm−0Ω−B
Frequency [Hz]
fdm−100Ω−B
900
800
700
600
500
4 6 8 10 12 14 16
Fault Distance [km]
Figure 6.21. Comparison of estimated damped frequencies with the effect of load and
fault resistance
119
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
Figure 6.22, shows the effect of load and fault resistance on calculated
undamped transient frequency. It can be seen that, with the effect of 5
MVA and 12.5 MVA loads, the calculated undamped transient frequency
is almost not affected by fault resistance. Based on the result, it is found
that the undamped frequency with a higher load (12.5 MVA) has higher
frequency than the undamped frequency with a lower load (5 MVA). The
results shown in Figure 6.23, tell us that the estimated undamped fre-
quency with high load and lower fault resistance is greater than the esti-
mated undamped frequency having low load and high fault resistance.
In Figure 6.24(d), it can be seen that, the estimated and calculated
damping factor proved to be identical with small error. When the fault
resistance in increased the damping factor is increased as well. The re-
sult shows that, the damping factor with load B is lower than the damping
factor with load A.
The bar charts shown in Figures 6.25 and 6.26, present a comparison
of earth fault distance estimation error by GM algorithms when calculat-
ing earth fault distance in the simulated network. Appendix C, Table C8
shows tabulated overall distance estimation results. A single line to earth
fault was simulated at 4 km, 10 km and 16 km from substation for loads
5 MVA and 12.5 MVA. Earth faults were simulated with fault resistance
values of 0 Ω , 25 Ω, 50 Ω and 100 Ω.
In overall, based on the results, we found that the fault distance esti-
mation error is less than 1 km when the fault resistance is below 50 Ω.
However, when fault resistance is 100 Ω, the estimation error is increased.
In Figures 6.25 and 6.26, when fault happens at 4 km and with the fault
resistance less than 25 Ω, GM3 algorithm calculates the fault distance
with the lower error compared with the other algorithms. Unfortunately,
when the fault distance is simulated at 10 km, GM3 algorithm gives the
highest fault estimation error among all algorithms. It can be seen as
well that GM3 algorithm has the highest estimation error when the earth
fault happens at fault distance 16 km and with fault resistance 100 Ω.
The results from GM1, GM2 and GM4 algorithm are almost identical for
all cases of simulated earth faults. In overall, the MAE error of estimated
fault distance using proposed GM algorithms is less than 0.6 km, when
fault resistance is 50 Ω or less.
120
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
1100
fnc−0Ω−A
fnc−25Ω−A
1000 fnc−50Ω−A
Frequency [Hz]
fnc−100Ω−A
900
800
700
600
500
4 6 8 10 12 14 16
Fault Distance [km]
1100
fnc−0Ω−B
fnc−25Ω−B
1000 fnc−50Ω−B
Frequency [Hz]
fnc−100Ω−B
900
800
700
600
500
4 6 8 10 12 14 16
Fault Distance [km]
1100
fnc−0Ω−A
fnc−100Ω−A
1000 fnc−0Ω−B
Frequency [Hz]
fnc−100Ω−B
900
800
700
600
500
4 6 8 10 12 14 16
Fault Distance [km]
Figure 6.22. Comparison of calculated undamped frequencies with the effect of load and
fault resistance
121
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
1100
fnm−0Ω−A
fnm−25Ω−A
1000 fnm−50Ω−A
Frequency [Hz]
fnm−100Ω−A
900
800
700
600
500
4 6 8 10 12 14 16
1100
fnm−0Ω−B
fnm−25Ω−B
1000 fnm−50Ω−B
fnm−100Ω−B
Frequency [Hz]
900
800
700
600
500
4 6 8 10 12 14 16
Fault Distance [km]
1100
fnm−0Ω−A
fnm−100Ω−A
1000 fnm−0Ω−B
fnm−100Ω−B
Frequency [Hz]
900
800
700
600
500
4 6 8 10 12 14 16
Fault Distance [km]
Figure 6.23. Comparison of estimated undamped frequencies with the effect of load and
fault resistance
122
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
ζc−0Ω−A
0.3 ζc−25Ω−A
ζc−50Ω−A
ζc−100Ω−A
0.25 ζc−0Ω−B
ζc−25Ω−B
0.2 ζc−50Ω−B
ζc−100Ω−B
0.15
0.1
0.05
4 6 8 10 12 14 16
Fault Distance[km]
ζm−0Ω−A
ζm−25Ω−A
ζm−50Ω−A
0.3 ζm−100Ω−A
ζm−0Ω−B
ζm−25Ω−B
0.25 ζm−50Ω−B
ζm−100Ω−B
0.2
0.15
0.1
0.05
4 6 8 10 12 14 16
Fault Distance [km]
ζc−0Ω−A
ζc−25Ω−A
ζc−50Ω−A
0.3 ζc−100Ω−A
ζm−0Ω−A
ζm−25Ω−A
0.25 ζm−50Ω−A
ζm−100Ω−A
0.2
0.15
0.1
0.05
4 6 8 10 12 14 16
Fault Distance [km]
ζc−0Ω−B
ζc−25Ω−B
ζc−50Ω−B
0.3 ζc−100Ω−B
ζm−0Ω−B
ζm−25Ω−B
0.25 ζm−50Ω−B
ζm−100Ω−B
0.2
0.15
0.1
0.05
4 6 8 10 12 14 16
Figure 6.24. Calculated and estimated damping factor (ζ) with the effect of load and fault
resistance.
123
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
2
GM1
GM2
1.5 GM3
GM4
1
Error[km]
0.5
−0.5
−1
4 10 16
Fault Distance[km]
2
GM1
GM2
1.5 GM3
GM4
1
Error[km]
0.5
−0.5
−1
4 10 16
Fault Distance[km]
2
GM1
GM2
1.5 GM3
GM4
1
Error[km]
0.5
−0.5
−1
4 10 16
Fault Distance[km]
2
GM1
GM2
1.5 GM3
GM4
1
Error[km]
0.5
−0.5
−1
4 10 16
Fault Distance[km]
Figure 6.25. Distance estimates with the effect of 5 MVA load and fault resistances for
general model networks.
124
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
2
GM1
GM2
1.5 GM3
GM4
1
Error[km]
0.5
−0.5
−1
4 10 16
Fault Distance[km]
2
GM1
GM2
1.5 GM3
GM4
1
Error[km]
0.5
−0.5
−1
4 10 16
Fault Distance[km]
2
GM1
GM2
1.5 GM3
GM4
1
Error[km]
0.5
−0.5
−1
4 10 16
Fault Distance[km]
2
GM1
GM2
1.5 GM3
GM4
1
Error[km]
0.5
−0.5
−1
4 10 16
Fault Distance[km]
Figure 6.26. Distance estimates with the effect of 12 MVA load and fault resistances for
general model networks.
125
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
In this section, we discuss the accuracy of the single line to earth fault lo-
cation algorithms that has been explained in Chapter 5, Section 5.2. The
fault location algorithms used in this section are the accurate model EM1
and for the comparison, the simple model GM2. The accuracy of the meth-
ods was tested using the simulation model and the network example with
overhead lines. The performance accuracy of EM1 and GM2 algorithm
was investigated in various fault conditions. For EM algorithm, the per-
formance accuracy was defined by the mean absolute error (MAEEM ) as
expressed as in Equation (6.3) below, while for GM algorithm, the perfor-
mance accuracy was defined as given in Equation (6.1).
n n
1 1
M AEEM = |EMi − li | = |∆EMi | (6.3)
n n
i=1 i=1
The single line diagram of the simulated unearthed medium voltage (MV)
network and its simulation model is shown in Figure 6.27. The simula-
tion model is based on the sequence component networks for the exact-π
solution for the earth fault transients. The system is 21 kV overhead
lines unearthed MV radial network and the line parameters are as fol-
lows: Rl,p =0.6 Ω/km, Rl,0 =1.3 Ω/km, Ll,p =1.0 mH/km, Ll,0 = 5.0 mH/km,
Cl,p =10.71 nF/km, Cl,0 =6.12 nF/km, RT =0.0216 Ω, LT =2.8 mH. We assume
the negative sequence parameters to be equal to the positive sequence
ones.
126
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
1 17km
2 15km
3 22km
110/21kV
110 kV
Supply 4 24km
25MVA
5 20km
M
Zf
E L1T
Cl|5 Cl|5 C1|j C1|j "1"
2 2 2 2
R1T
j=1...4
L2T
C2|5 C2|5 C2|j C2|j "2"
2 2 2 2
R2T
j=1...4
Figure 6.27. A schematic network diagram and its simulation model for a single phase
to earth fault in a distribution network. M denotes measurement point, F
refers to the fault location and Zf is the fault impedance. L0|5 , L1|5 , L2|5
are zero sequence, positive sequence and negative sequence inductance of
faulty line located in front of fault point (F ). L0|j , L1|j are zero sequence and
positive sequence inductance of sound line. LT and RT is the inductance and
resistance of the transformer, capacitances and resistances correspondingly.
The lines were modelled as a chain of exact π-model cells (1 cell per
10 km of the line). The network was 5 feeders and the earth fault was
simulated at feeder 5. The total length of the network is 98 km. All
faults with several different fault conditions were simulated with the soft-
ware package ATP (Alternative Transients Program), version of EMTP
127
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
program where the circuit was realized using ATPDraw. The sampling
frequency was 20 kHz. In the simulation analysis, the load was either
0.4 MW, 0.8 MW, 2 MW or zero.
Figure 6.28 shows the results of the simulation without the effect of load
and with the effect of fault distance, fault resistance and fault inception
angle. The tabulated results can be found in Appendix C Table C9, C10
and C11, respectively. In order to study the performance of the tested
algorithms, in Figure 6.28(a), an earth fault was simulated at 10 points
starting from 2 km from measuring point (M) with 2 kilometer steps. In
case of the effect of fault resistances as shown in Figure 6.28(b), 3 different
locations of earth fault was simulated which are 4 km, 10 km and 16 km
from the measuring point. Only two fault distances were simulated for
the case of fault inception angle variation, which are 4 km and 16 km, as
presented in Figure 6.28(c). Based on the MAE results given in Appendix
C Table C9, C10 and C11, it can be seen, that the EM1 algorithm provides
better results than the GM2 algorithm.
As shown in Figure 6.28(a), the estimation error is increased with the in-
crease of fault distance for both algorithms. GM2 algorithm recorded the
highest error, when earth fault was located at 20 km from the substation.
For the earth fault distance of 4km, 10km and 16km, the estimation error
is further increased when fault distance is increased with the increased
of fault resistance, as we seen in Figure 6.28(b).
In Figure 6.28(c), it seems that fault inception angle does not affect both
algorithms when fault happens at fault distance 4 km, however, when
fault is located at 16 km the error is increasing for EM1 algorithm but
decreasing for GM2 algorithm.
If we carefully observe the results in Figure 6.28, Table C10 and Ta-
ble C11 in Appendix C, in case of fault distance located at 16 km from
fault distance, for GM2 algorithm estimation result, indirectly the effect
of increased of fault resistance and fault inception angle at 0 degree is de-
creasing the estimated undamped charge transient frequency (fnm ) and
at the same time improving the algorithm accuracy.
128
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
2
EM1
GM2
1.5
0.5
Error(km)
0
−0.5
−1
−1.5
−2
2 4 6 8 10 12 14 16 18 20
Distance(km)
1
Error(km)
0.5
−0.5
−1
−1.5
−2
410 4 4 4100 10 10 10 10 1610 16 16 16
30 50 10 30 50 100 30 50 100
Distance(km)Rf (Ω)
2
EM1
GM2
1.5
1
Error(km)
0.5
−0.5
−1
−1.5
−2
490 445 40 1690 1645 160
Distance(km)Inception angle
◦
Figure 6.28. The effect of fault distance, fault resistance and inception angle on fault
distance calulation. EM1 denotes the distance calculation error using ex-
act model 1 and GM2 denotes the distance calculation error using general
model 2.
129
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
Figure 6.29 shows the results of fault location performance with the
effect of resistive load located at the end of each feeder. The tabulated
results can be found in Appendix C Table C12. In this part of simulation,
each feeder has same amount of load with three different cases which
were 400 kW, 800 kW and 2 MW.
Based on the MAE results given in Appendix C Table C12, it can be
seen, that the EM1 algorithm again provides better results than the GM2
algorithm. In overall, both algorithms have an increase of total MAE
when resistive load is added to the simulated network.
As shown in Figure 6.29, the fault distance error seems to be increased
when the load is increased. The more far the fault happens from the
substation the higher is the fault distance error. In Figure 6.29, in some
cases corresponding to GM2, we can see that the fault distance error is
slightly decreased when fault resistance is increased. In overall, the fault
distance estimation error is less than 1.5 km.
Figure 6.30 shows the results of fault fault location performance with
the effect of inductive load located at the end of each feeder. The tabulated
results can be found in Appendix C Table C13. In this part of simulation,
each feeder has same amount of load with three difference cases which
were 400 kVAr, 800 kVAr and 2 MVAr. Corresponding to the overall MAE
result shown in Appendix C Table C13, EM1 algorithm provides better
results than the GM2 algorithm.
It is found that the performance accuracy for EM1 algorithm is increased
with the effect of inductive load compared to than when it is been affected
by resistive load. However, in case of GM2 algorithm, the performance
accuracy is decreased with the effect of inductive load. In Figure 6.30, in
overall, the results show that the fault distance estimation error is less
than 1.5 km.
130
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
1.5
EM1
GM2
0.5
Error(km)
0
−0.5
−1
−1.5
230 250 2100 1030 1050 10100 1830 1850 18100
Distance(km)Rf (Ω)
(a) 400 kW
1.5
EM1
GM2
1
0.5
Error(km)
−0.5
−1
−1.5
230 250 2100 1030 1050 10100 1830 1850 18100
Distance(km)Rf (Ω)
(b) 800 kW
1.5
EM1
GM2
0.5
Error(km)
−0.5
−1
−1.5
230 250 2100 1030 1050 10100 1830 1850 18100
Distance(km)Rf (Ω)
(c) 2 MW
Figure 6.29. Fault location performance with the effect of resistive load located at the
end of each feeder. EM1 denotes a distance calculation error using exact
model 1 and GM2 denotes a distance calculation error using general model
2.
131
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
1.5
EM1
GM2
0.5
Error(km)
−0.5
−1
−1.5
230 250 2100 1030 1050 10100 1830 1850 18100
Distance(km)Rf (Ω)
0.5
Error(km)
−0.5
−1
−1.5
230 250 2100 1030 1050 10100 1830 1850 18100
Distance(km)Rf (Ω)
1.5
EM1
GM2
0.5
Error(km)
−0.5
−1
−1.5
230 250 2100 1030 1050 10100 1830 1850 18100
Distance(km)Rf (Ω)
(c) 2 MVAr
Figure 6.30. Fault location performance with the effect of inductive load located at the
end of each feeder. EM1 denotes a distance calculation error using 1 and
GM2 denotes a distance calculation error using general model 2.
132
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
Figure 6.31 shows the results of fault fault location performance with
the effect of load for 0.9 and 0.5 power factor. The tabulated results can be
found in Appendix C Table C14 and Table C15, respectively. In this part
of simulation, each feeder has same amount of load with two difference
cases which were 400 KVA and 800 KVA.
Based on the MAE result given in Appendix C Table C14 and Table C14,
it can be seen, that the EM1 algorithm provides again better results than
the GM2 algorithm. By comparison of MAE results in Appendix C Ta-
ble C14 and Table C14, we found that both algorithms have a decrease of
MAE when low power factor is used.
In this work, the most influential source of error was resistive load.
Theoretically, the load can be compensated by adding the corresponding
impedance to the algorithm models. However, the data of load devices and
their transient impedances is difficult to get. Another way is to use the
curve fitting technique to compensate the error with known load condi-
tions.
Another source of error is inaccuracy in estimating the undamped com-
plex frequency. This error is caused by the inaccuracy in estimation of the
damping attenuation through Hilbert transformation method and also by
an inaccurate detection of charge transient component using CWT. When
the magnitude of fault resistance or resistive loads is increased, the tran-
sients become more damped, with the result that the estimation of the
charge transient frequency becomes more difficult which affects the over-
all computation of fault distance. In following section, the investigation of
performance accuracy of proposed algorithms is performed using JMarti
(frequency-dependent) line model which is more practical compared to the
exact-π solution that we presented in this section.
Considering which one of the algorithms is easy to be implement in real
operation, the general model algorithm should be much easier to be imple-
ment compare to the exact model algorithm. It is because its required a
simple iteration process to estimate the fault distance. While exact model
algorithm required a heavy computation to solve a higher order equation,
which generally may not be easy.
133
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
1.5
EM1
GM2
1
0.5
Error(km)
0
−0.5
−1
−1.5
2 250 2100 1830 1850 18100
30
Distance(km)Rf (Ω)
1.5
EM1
GM2
1
0.5
Error(km)
−0.5
−1
−1.5
2 2 2100 1830 1850 18
30 50 100
Distance(km)Rf (Ω)
1.5
EM1
GM2
0.5
Error(km)
−0.5
−1
−1.5
2 2 2100 1830 1850 18100
30 50
Distance(km)Rf (Ω)
1.5
EM1
GM2
1
0.5
Error(km)
−0.5
−1
−1.5
230 250 2100 1830 1850 18
100
Distance(km)Rf (Ω)
Figure 6.31. Fault location performance with the effect of load with 0.9 and 0.5 PF lo-
cated at the end of each feeder. EM1 denotes the distance calculation error
using exact model 1 and GM2 denotes the distance calculation error using
general model 2.
134
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
n n
1 1
M AEW T = |W Ti − li | = |∆W Ti | (6.4)
n n
i=1 i=1
135
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
N23
5.6km
N1 N2
44km 6
N21 N22 0.1MVA
N3 N4 1MVA
3.54km 5.28km Zone C
32km
N19 N20
0.1MVA
3.89km 2.45km
N5 N6 1MVA
25MVA
36km 4
0.1MVA 0.1MVA
N18
N7 N8 1MVA
3.98km
28km
5
1MVA 0.1MVA
N9 N11 N12 N13
2.8km
N10
3.5km 5.4km 3.8km 2.1km
N14
3.9km
N15
2.7km
N16 N17
Zone B
2.6km
1 3
0.1MVA 0.1MVA 0.1MVA 0.1MVA 0.1MVA 0.1MVA
2 0.1MVA
8 0.1MVA 9
0.1MVA
N26 N27 N28
7
Legend: 3.98km 2.68km 8.58km
Fault position
11
Zone A
x Fault number 0.1MVA
N30 N31
10
NX Node number X 0.1MVA 2.76km 2.7km
Zone
0.1MVA 0.1MVA
Figure 6.32. A schematic diagram of 20 kV, 227.6 km, unearthed neutral MV network.
136
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
The results of single line to earth fault estimation with WT algorithm are
shown in Tables 6.1 and 6.2. In overall, MAE of WT algorithm is 1.09 km.
In Tables 6.1 and 6.2, it is shown that the fault location error is increas-
ing when the fault is located far from the measuring point and the fault
inception angle is small. Also the fault location error is higher when the
fault is located far from the measuring point and has a high fault resis-
tance. In overall, the fault distance estimation error is less than 5.2 km
and 6.52 km for simulation with test angle of 90◦ and 45◦ , respectively. It
is found that the test number with the longest fault distance (test num-
ber 5 (26.19 km)) calculated the highest fault location error. From our
experience in analyzing the transient signal, referring to Tables 6.1 and
6.2, when the magnitude of the fault resistance is increased the transient
signals become more damped and the fault location is more difficult.
Figures 6.34, 6.35 , 6.36 , 6.37 and 6.38, show results of calculated CWT
FFT amplitudes of transient at MV/LV substation for fault at inception
angles of 90◦ and 45◦ which are used in correct path identification algo-
rithm. Based on the correct path identification algorithm, the highest
recorded FFT amplitude will indicate which MV/LV substation should we
locate in order to determine the correct fault location if there are several
branches in the network. Also, based on the test results, another assump-
tion is that if the fault distance estimation shows a large error in its cal-
culation, the correct path algorithm can still give us a good guess about
the fault position. However, this assumption cannot be true for all fault
positions. For instance, in Table 6.1, the fault is located at test number 1
which is 1.4 km from measuring point and the LV substation detected is
at substation S23. However, the nearest LV substation to the fault point
is substation S10. Similarly, for fault located at test number 2, 3, 4, 7
and 10. Therefore, the assumption can be used only if the fault happens
exactly at the section which located at the end of the feeder branches such
as the faults located at test number 5, 6, 8, 9 and 11.
In this work, the zones were defined based on the location of the MV/LV
substation located in the same main branches. For examples, it can be
seen that MV/LV substations number S25, S26, S28, S29, S30 and S31
located in zone A, MV/LV substations number S14, S15, S16, S17 and S18
is located at zone B and MV/LV substations number S19, S20, S22 and
S23 is located at zone C.
137
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
Table 6.1. Test results of single line to earth fault for 90◦ fault inception angle using the
WT algorithm.
Test Test
Substation detected Fault location error (km) MAE
number Section Candidates
(distance) / Node Rf 0 Rf 25 Rf 50 Rf 0 Rf 25 Rf 50 (km)
√
1 (1.4km) N9-N10 N9-N10 S23 S23 S23 -0.07 -0.43 -1.75 0.75
√
N11-N13
2 (11.7km) N12 N24-N25 S20 S20 S20 0.53 0.52 -0.02 0.36
N24-N26
√
N15-N16
3 (22.8km) N15-N16 S18 S18 S18 -0.07 -0.38 -1.44 0.63
N19-N21
N13-N14
√
N13-N19
4 (17.5km) N13-N19 N27-N28 S20 S20 S20 0.58 -1.46 -0.84 0.96
N27-N29
N30-N31
N16-N17
√
N16-N18
5 (26.19km) N16-N18 S18 S18 S18 -2.09 -3.09 -5.11 3.43
N21-N23
N21-N22
N21-N22
√
N21-N23
6 (25.83km) N21-N23 S23 S23 S23 -0.87 -1.24 -1.64 1.25
N16-N17
N16-N18
N11-N12
√
7 (10.18km) N24 N11-N25 S29 S25 S25 0.75 0.05 0.91 0.57
√
N11-N26
N12-N13
√
8 (12.53km) N24-N25 N24-N25 S25 S25 S25 0.64 0.58 0.25 0.49
N24-N26
N14-N15
N19-N20
9 (19.72km) N27-N29 N19-N21 S29 S29 S29 -0.6 -1.44 -2.16 1.4
N27-N28
√
N27-N29
N12-N13
√
N24-N27
10 (14.16km) N26 S28 S28 S28 0.52 0.37 -0.78 0.52
N24-N30
N24-N25
√
N27-N28
N27-N29
11 (21.72km) N27-N28 N16-N17 S28 S28 S28 0.34 0.18 -3.11 1.21
N16-N18
N19-N21
138
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
Table 6.2. Test results of single line to earth fault for 45◦ fault inception angle
Test Test
Substation detected Fault location error (km) MAE
number Section Candidates
(distance) / Node Rf 0 Rf 25 Rf 50 Rf 0 Rf 25 Rf 50 (km)
√
1 (1.4km) N9-N10 N9-N10 S23 S23 S28 0.15 -1.27 -2.06 1.16
√
N11-N13
2 (11.7km) N12 N24-N25 S15 S15 S14 0.55 0.83 -0.23 0.54
N24-N26
√
N15-N16
3 (22.8km) N15-N16 S18 S18 S18 0.52 -0.08 -0.42 0.34
N19-N21
N13-N14
√
N13-N19
4 (17.5km) N13-N19 N27-N28 S20 S20 S20 -0.07 0.74 0.46 0.42
N27-N29
N30-N31
N16-N17
√
N16-N18
5 (26.19km) N16-N18 S18 S18 S18 -2 -2.26 -6.52 3.59
N21-N23
N21-N22
N21-N22
√
N21-N23
6 (25.83km) N21-N23 S23 S23 S23 -1.38 -1.9 0.03 1.1
N16-N17
N16-N18
N11-N12
7 (10.18km) N24 N11-N25 S31 S31 S31 0.96 0.45 -0.16 0.52
√
N11-N26
N12-N13
√
8 (12.53km) N24-N25 N24-N25 S25 S25 S25 0.68 0.92 0.11 0.57
N24-N26
N14-N15
N19-N20
9 (19.72km) N27-N29 N19-N21 S29 S29 S29 -0.67 -0.94 -3.84 1.82
N27-N28
√
N27-N29
N12-N13
√
N24-N27
10 (14.16km) N26 S29 S29 S29 -0.84 0.68 -2.32 1.28
N24-N30
N24-N25
√
N27-N28
N27-N29
11 (21.72km) N27-N28 N16-N17 S28 S28 S28 -0.19 -0.14 -2.44 0.92
N16-N18
N19-N21
139
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
200
Rf=0
R =25
f
150 R =50
f
Amplitude[FFT]
100
50
0
S10 S12 S13 S14 S15 S16 S17 S18 S19 S20 S22 S23 S25 S26 S28 S29 S30 S31
Subtation
80
60
40
20
0
S10 S12 S13 S14 S15 S16 S17 S18 S19 S20 S22 S23 S25 S26 S28 S29 S30 S31
Substation
Figure 6.34. CWT FFT amplitude of MV/LV substation for fault test number 1 with fault
inception angle of 90◦ and 45◦ . The arrow indicates the substation with the
highest amplitude.
From Figure 6.35, for example CWT FFT amplitude of MV/LV substa-
tions for test number 7, it can be seen that MV/LV substations number
S25, S26, S28, S29, S30 and S31 estimated almost the same CWT FFT
amplitude value, where substation S29 estimated the highest value. S25,
S26, S28, S29, S30 and S31 is actually placed inside zone A, as shown in
Figure 6.32. Therefore, if any one of these substations indicated the high-
est amplitude for test number 7 earth fault, we can make an assumption
that the correct path should go through that zone. Similar assumption
can be used for test number 8, 9, 10 and 11. The assumption also can
further be used for fault located at other zone, for instant test number 3,
referring to result in Figure 6.36, where it shows that CWT FFT ampli-
tude of MV/LV substations for S15, S16, S17 and S18 have close value and
fault is located at Zone B. Therefore, based from the earlier studies of the
simulation results, zoning the branches gives the first indication of the
correct path of fault passage in the network. However, this assumption
cannot be used for exact correct path of fault location since more branches
were located inside the zone. From the simulation results in Figures 6.34,
6.35, 6.36, 6.37 and 6.38, the correct path algorithm managed to identify
the correct LV substation to be used as a guidance of calculating the fault
distance from the measuring point towards the fault position. Figures
6.39, 6.40, 6.41, 6.42 and 6.43 show an example of schematic diagram of
correct path identification for test number 2, 3, 4, 6, and 8, respectively.
140
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
200
R =0
f
180 Rf=25
R =50
160 f
Amplitude [FFT]
140
120
100
80
60
40
20
0
S10 S12 S13 S14 S15 S16 S17 S18 S19 S20 S22 S23 S25 S26 S28 S29 S30 S31
Substation
Amplitude [FFT]
140
120
100
80
60
40
20
0
S10 S12 S13 S14 S15 S16 S17 S18 S19 S20 S22 S23 S25 S26 S28 S29 S30 S31
Substation
140
120
100
80
60
40
20
0
S10 S12 S13 S14 S15 S16 S17 S18 S19 S20 S22 S23 S25 S26 S28 S29 S30 S31
Substation
160
Amplitude [FFT]
140
120
100
80
60
40
20
0
S10 S12 S13 S14 S15 S16 S17 S18 S19 S20 S22 S23 S25 S26 S28 S29 S30 S31
Subtation
150
100
50
0
S10 S12 S13 S14 S15 S16 S17 S18 S19 S20 S22 S23 S25 S26 S28 S29 S30 S31
Substation
Figure 6.35. CWT FFT amplitude of MV/LV substation for fault test number 7, 2, 8, 10
and 4 with fault inception angle of 90◦ . The arrow indicates the substation
with the highest amplitude.
141
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
220
Rf=0
200 Rf=25
180 Rf=50
160
Amplitude [FFT]
140
120
100
80
60
40
20
0
S10 S12 S13 S14 S15 S16 S17 S18 S19 S20 S22 S23 S25 S26 S28 S29 S30 S31
Substation
Rf=0
250 R =25
f
Rf=50
Amplitude [FFT]
200
150
100
50
0
S10 S12 S13 S14 S15 S16 S17 S18 S19 S20 S22 S23 S25 S26 S28 S29 S30 S31
Substation
200
150
100
50
0
S10 S12 S13 S14 S15 S16 S17 S18 S19 S20 S22 S23 S25 S26 S28 S29 S30 S31
Substation
RF=0
RF=25
200 RF=50
Amplitude [FFT]
150
100
50
0
S10 S12 S13 S14 S15 S16 S17 S18 S19 S20 S22 S23 S25 S26 S28 S29 S30 S31
Substation
Rf=0
300
Rf=25
Rf=50
250
Amplitude [FFT]
200
150
100
50
0
S10 S12 S13 S14 S15 S16 S17 S18 S19 S20 S22 S23 S25 S26 S28 S29 S30 S31
Substation
Figure 6.36. CWT FFT amplitude of MV/LV substation for fault test number 9, 11, 3, 6
and 5 with fault inception angle of 90◦ . The arrow indicates the substation
with the highest amplitude.
142
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
200
Rf=0
180 Rf=25
R =50
160 f
Amplitude [FFT]
140
120
100
80
60
40
20
0
S10 S12 S13 S14 S15 S16 S17 S18 S19 S20 S22 S23 S25 S26 S28 S29 S30 S31
Substation
Amplitude [FFT]
140
120
100
80
60
40
20
0
S10 S12 S13 S14 S15 S16 S17 S18 S19 S20 S22 S23 S25 S26 S28 S29 S30 S31
Substation
140
120
100
80
60
40
20
0
S10 S12 S13 S14 S15 S16 S17 S18 S19 S20 S22 S23 S25 S26 S28 S29 S30 S31
Substation
140
120
100
80
60
40
20
0
S10 S12 S13 S14 S15 S16 S17 S18 S19 S20 S22 S23 S25 S26 S28 S29 S30 S31
Substation
160
140
120
100
80
60
40
20
0
S10 S12 S13 S14 S15 S16 S17 S18 S19 S20 S22 S23 S25 S26 S28 S29 S30 S31
Substation
Figure 6.37. CWT FFT amplitude of MV/LV substation for fault test number 7, 2, 8, 10
and 4 with fault inception angle of 45◦ . The arrow indicates the substation
with the highest amplitude.
143
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
220
Rf=0
200 Rf=25
180 Rf=50
160
Amplitude [FFT]
140
120
100
80
60
40
20
0
S10 S12 S13 S14 S15 S16 S17 S18 S19 S20 S22 S23 S25 S26 S28 S29 S30 S31
Substation
150
100
50
0
S10 S12 S13 S14 S15 S16 S17 S18 S19 S20 S22 S23 S25 S26 S28 S29 S30 S31
Substation
150
100
50
0
S10 S12 S13 S14 S15 S16 S17 S18 S19 S20 S22 S23 S25 S26 S28 S29 S30 S31
Substation
140
120
100
80
60
40
20
0
S10 S12 S13 S14 S15 S16 S17 S18 S19 S20 S22 S23 S25 S26 S28 S29 S30 S31
Substation
150
100
50
0
S10 S12 S13 S14 S15 S16 S17 S18 S19 S20 S22 S23 S25 S26 S28 S29 S30 S31
Substation
Figure 6.38. CWT FFT amplitude of MV/LV substation for fault test number 9, 11, 3, 6
and 5 with fault inception angle of 45◦ . The arrow indicates the substation
with the highest amplitude.
144
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
N1 N2
44km
N3 N4 1MVA N23
32km 5.6km
N9 N10 1MVA N11 N12 N13 N14 0.1MVA N15 0.1MVA N16 N17 0.1MVA
2.8km 3.5km 5.4km 3.8km 2.1km 3.9km 2.7km 2.6km
3.88km
N26 0.1MVA N27 N28 0.1MVA
3.98km 2.68km 8.58km
Legend:
N30 N31 0.1MVA
Fault position
0.1MVA 2.76km 2.7km
x Fault number
Correct path
NX Node number X 0.1MVA 0.1MVA
N1 N2
44km
N3 N4 1MVA N23
32km 5.6km
N9 N10 1MVA N11 N12 N13 N14 0.1MVA N15 0.1MVA N16 N17 0.1MVA
2.8km 3.5km 5.4km 3.8km 2.1km 3.9km 2.7km 2.6km
0.1MVA 3
0.1MVA N24 0.1MVA N25 0.1MVA N29 0.1MVA 0.1MVA
0.1MVA
3.88km 6.0km 5.48km
145
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
N1 N2
44km
N3 N4 1MVA N23
32km 5.6km
25MVA
N5 N6 1MVA N21 N22 0.1MVA
36km 3.54km 5.28km
4
N9 N10 1MVA N11 N12 N13 N14 0.1MVA N15 0.1MVA N16 N17 0.1MVA
2.8km 3.5km 5.4km 3.8km 2.1km 3.9km 2.7km 2.6km
Legend:
N30 N31 0.1MVA
Fault position
0.1MVA 2.76km 2.7km
x Fault number
Correct path
NX Node number x 0.1MVA 0.1MVA
N1 N2
44km
N3 N4 1MVA N23
32km 5.6km
6 0.1MVA
N5 N6 1MVA N21 N22
25MVA
36km 3.54km 5.28km
N9 N10 1MVA N11 N12 N13 N14 0.1MVA N15 0.1MVA N16 N17 0.1MVA
2.8km 3.5km 5.4km 3.8km 2.1km 3.9km 2.7km 2.6km
146
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
N1 N2
44km
N3 N4 1MVA N23
32km 5.6km
N9 N10 1MVA N11 N12 N13 N14 0.1MVA N15 0.1MVA N16 N17 0.1MVA
2.8km 3.5km 5.4km 3.8km 2.1km 3.9km 2.7km 2.6km
Legend:
N30 N31 0.1MVA
Fault position 0.1MVA 2.76km 2.7km
x Fault number
Correct path
NX Node number x 0.1MVA 0.1MVA
147
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
1 2
f1
DT25 Legend:
3 4 F1 a −Length a (km)
f2
A −Load A (kVA)
DT26 DTX −Distribution
5 6 F2 23 transformer(DT) X
f3
l
25MVA DT27 DT12
110kV
Supply 7 8 F3 L
f4
17 18 19 20 21 22 24 25 26 27 28 29
DT28 g h i j k m n o p q r
f
9 10 F4 DT6 DT7 DT8 DT9 DT10
DT11
DT13 DT14 DT15 DT16 DT17 DT18
f5
F G H I J K M N O P Q R
DT29
F5
11 12 13 14 15 16 s 30 31 32 33 34 35
a b c d e t u v w x
DT1 DT2 DT3 DT4 DT5 DT19 DT20 DT21 DT22 DT23 DT24
A B C D E S T U V W X
148
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
149
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
In this work, the transient signal is measured from the secondary side
of MV/LV substation number 24 (DT24), as shown in Figure 6.44. In order
to examine the performance of the algorithms, a large number of data sets
is required. The variation of the data sets used for formula construction
and testing in this work is given in Table 6.3 and Table 6.4, respectively .
In Table 6.3, the network model was simulated with 5 different load con-
ditions. For each load condition given in Table 6.3, the fault distance, fault
resistance and fault inception angle were varied. Similarly in Table 6.4,
the network model was simulated with load conditions as in Table 6.3 but
with different fault distance.
Table 6.3. Data generation for MRA and NN formulation. Fault resistances are 0, 25 and
50 Ω. Inception angles are 90, 60 and 30 ◦ for load 1 and load 2. 90 and 45 ◦
for load 3 , 4 and 5.
Load number 1 2 3 4 5
Number of variations
Fault distance 13 13 2 2 2
Fault resistance 3 3 3 3 3
Inception angle 3 3 2 2 2
Total 117 117 12 12 12
Table 6.4. Data generation for testing the MRA and NN algorithm. Fault resistances are
0, 25 and 50 Ω. Inception angles are 90 and 45 ◦ .
Load number 1 2 3 4 5
Number of variations
Fault distance 2 2 1 1 1
Fault resistance 2 2 2 2 2
Inception angle 3 3 3 3 3
Total 12 12 6 6 6
Data sets 42
Using the data in Table 6.3, the formula obtained for MRA algorithm,
which corresponds to Equation (5.33), is given as:
150
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
as:
l = AN N (A, S, L) (6.6)
n n
1 1
M AEM RA = |M RAi − li | = |∆M RAi | (6.7)
n n
i=1 i=1
n n
1 1
M AEN N = |N Ni − li | = |∆N Ni | (6.8)
n n
i=1 i=1
where ∆M RAi and ∆N Ni is the absolute error fault distance using MRA
and NN algorithm, M RAi and N Ni is estimated fault distance using MRA
and NN algorithm, li is exact fault distance and n is numbers of observa-
tion of parameter error.
151
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
200
180
160
Amplitude [FFT]
140
120
100
80
60
40
20 DT5 DT12 DT17 DT24
0
0 2 4 6 8 10 12 14 16 18 20
Fault Distance [km]
(a) Rf =0 Ω
200
180
160
Amplitude [FFT]
140
120
100
80
60
40
20 DT12 DT17 DT24 DT5
0
0 2 4 6 8 10 12 14 16 18 20
Fault Distance[km]
(b) Rf =25 Ω
200
180
160
Amplitude [FFT]
140
120
100
80
60
40
20 DT5 DT12 DT17 DT24
0
0 2 4 6 8 10 12 14 16 18
Fault Distance[km]
(c) Rf =50 Ω
152
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
200
180
160
Amplitude [FFT]
140
120
100
80
60
40
20
DT5 DT12 DT17 DT24
0
0 2 4 6 8 10 12 14 16 18 20
Fault Distance [km]
(a) Rf =0 Ω
200
180
160
Amplitude [FFT]
140
120
100
80
60
40
20 DT5 DT12 DT17 DT24
0
0 2 4 6 8 10 12 14 16 18 20
Fault Distance [km]
(b) Rf =25 Ω
200
180
160
Amplitude [FFT]
140
120
100
80
60
40
(c) Rf =50 Ω
153
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
200
180
(a) Rf =0 Ω
200
180
160
Amplitude [FFT]
140
120
100
80
60
40
20
DT5 DT12 DT17 DT24
0
0 2 4 6 8 10 12 14 16 18 20
Fault Distance[km]
(b) Rf =25 Ω
200
180
160
Amplitude [FFT]
140
120
100
80
60
40
(c) Rf =50 Ω
154
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
Figures 6.49 and 6.50 show the results of fault location performance by
using MRA and NN algorithms at fault inception angles of 90 and 45 de-
grees. The tabulated results can be found in Appendix C, Table C16 and
C17, respectively. In comparison of calculated MAE of both algorithms, it
shows that NN algorithm performed slightly better than MRA algorithm.
Both algorithms have fault location error less than 1 km. It can be seen
that, fault resistance and fault inception angle have a small effect on the
result of the algorithms. It seems that the result of fault distance estima-
tion from simulated earth fault at 45 degrees of fault inception angle is
slightly lower than in the case of 90 degrees of fault inception angle.
Figures 6.51 and 6.52 show the results of fault location performance
by using GM(GM1, GM2, GM4) algorithms at fault inception angles of
90 and 45 degrees. The tabulated results can be found in Appendix C
Table C18 and C19, respectively. The results show that GM2 algorithm
was a bit better than GM1 and GM4 algorithms for fault inception angle
of 90 degree, with calculated MAE 0.88 km. In case of fault inception
angle of 45 degree, GM1 performed slightly better than the other two GM
algorithms. It is found that the results of MAE for all GM algorithms
worked better at fault inception angle of 45 degrees. In case of simulated
fault with inception angle of 90 degrees, the distance estimation error is
highest at fault distances of 7.86 km and 13.56 km. However, in the case
of fault simulated with inception angle of 45 degrees, the error of fault
distance 15.59 km is the highest.
Comparison of distance estimation errors calculated using EM (EM1,
EM2, EM3) algorithms with the effect of fault inception angle and fault
resistance are shown in Figures 6.53 and 6.54. Calculated MAE of EM
algorithms can be found in Appendix C, Table C20 and C21. From the
results, with the effect of fault inception angle of 90 degrees, it shows
that the MAE of EM3 algorithm is the lowest compared with the other
EM algorithms. But EM2 algorithm shows a slightly better performance
when fault inception angle of 45 degrees is used. EM algorithms distance
estimated error is less than 2 km when simulated earth fault is not af-
fected by the fault resistance. In overall, EM algorithms fault distance
estimation error was less than 4 km.
155
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
1.5
MRA
NN
1
0.5
Error(km)
−0.5
−1
−1.5
1.84 2.34 4.45 7.86 8.07 13.56 15.59
Distance(km)
(a) Rf =0 Ω
1.5
MRA
NN
1
0.5
Error(km)
−0.5
−1
−1.5
1.84 2.34 4.45 7.86 8.07 13.56 15.59
Distance(km)
(b) Rf =25 Ω
1.5
MRA
NN
1
0.5
Error(km)
−0.5
−1
−1.5
1.84 2.34 4.45 7.86 8.07 13.56 15.59
Distance(km)
(c) Rf =50 Ω
Figure 6.49. Comparison of the distance estimate errors obtained from multiple regres-
sion (MRA) and neural network (NN) algorithm with the effect of fault resis-
tance (Rf ). Fault inception angle is 90◦ . MRA denotes a distance calculation
error using multiple regression and NN denotes a distance calculation error
using neural network algorithm.
156
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
1
MRA
NN
0.5
Error(km)
0
−0.5
−1
1.84 2.34 4.45 7.86 8.07 13.56 15.59
Distance(km)
(a) Rf =0 Ω
1
MRA
NN
0.5
Error(km)
−0.5
−1
1.84 2.34 4.45 7.86 8.07 13.56 15.59
Distance(km)
(b) Rf =25 Ω
1
MRA
NN
0.5
Error(km)
−0.5
−1
1.84 2.34 4.45 7.86 8.07 13.56 15.59
Distance(km)
(c) Rf =50 Ω
Figure 6.50. Comparison of the distance estimate errors obtained from multiple regres-
sion (MRA) and neural network (NN) algorithm with the effect of fault resis-
tance (Rf ). Fault inception angle is 45◦ . MRA denotes a distance calculation
error using multiple regression and NN denotes a distance calculation error
using neural network algorithm.
157
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
4
GM1
GM2
3 GM4
Error(km) 2
−1
−2
1.84 2.34 4.45 7.86 8.07 13.56 15.59
Distance(km)
(a) Rf =0 Ω
4
GM1
GM2
3 GM4
2
Error(km)
−1
−2
1.84 2.34 4.45 7.86 8.07 13.56 15.59
Distance(km)
(b) Rf =25 Ω
4
GM1
GM2
3 GM4
2
Error(km)
−1
−2
1.84 2.34 4.45 7.86 8.07 13.56 15.59
Distance(km)
(c) Rf =50 Ω
Figure 6.51. Comparison of the distance estimate errors obtained from general model
(GM1,GM2,GM4) algorithms with the effect of fault resistance (Rf ). Signal
measured from MV/LV substation and fault inception angle is 90◦ . GM1
denotes a distance calculation error using general model 1, GM2 denotes
a distance calculation error using general model 2 and GM4 denotes a dis-
tance calculation error using general model 4.
158
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
2
GM1
1.5 GM2
GM4
1
0.5
Error(km)
0
−0.5
−1
−1.5
−2
−2.5
−3
1.84 2.34 4.45 7.86 8.07 13.56 15.59
Distance(km)
(a) Rf =0 Ω
2
GM1
1.5 GM2
GM4
1
0.5
Error(km)
0
−0.5
−1
−1.5
−2
−2.5
−3
1.84 2.34 4.45 7.86 8.07 13.56 15.59
Distance(km)
(b) Rf =25 Ω
2
GM1
1.5 GM2
GM4
1
0.5
Error(km)
0
−0.5
−1
−1.5
−2
−2.5
−3
1.84 2.34 4.45 7.86 8.07 13.56 15.59
Distance(km)
(c) Rf =50 Ω
Figure 6.52. Comparison of the distance estimate errors obtained from general model
(GM1,GM2,GM4) algorithms with the effect of fault resistance (Rf ). Signal
measured from MV/LV substation and fault inception angle is 45◦ . GM1
denotes a distance calculation error using general model 1, GM2 denotes
a distance calculation error using general model 2 and GM4 denotes a dis-
tance calculation error using general model 4.
159
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
4
EM1
EM2
EM3
3
Error(km) 2
−1
−2
1.84 2.34 4.45 7.86 8.07 13.56 15.59
Distance(km)
(a) Rf =0 Ω
4
EM1
EM2
EM3
3
2
Error(km)
−1
−2
1.84 2.34 4.45 7.86 8.07 13.56 15.59
Distance(km)
(b) Rf =25 Ω
4
EM1
EM2
EM3
3
2
Error(km)
−1
−2
1.84 2.34 4.45 7.86 8.07 13.56 15.59
Distance(km)
(c) Rf =50 Ω
Figure 6.53. Comparison of the distance estimate errors obtained from exact model
(EM1-EM3) algorithms with the effect of fault resistance (Rf ). Signal mea-
sured from MV/LV substation and fault inception angle is 90◦ . EM1 denotes
a distance calculation error using exact model 1, EM2 denotes a distance
calculation error using exact model 2 and EM3 denotes a distance calcula-
tion error using exact model 3.
160
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
2
EM1
EM2
1.5 EM3
1
0.5
Error(km)
0
−0.5
−1
−1.5
−2
−2.5
−3
1.84 2.34 4.45 7.86 8.07 13.56 15.59
Distance(km)
(a) Rf =0 Ω
2
EM1
1.5 EM2
EM3
1
0.5
Error(km)
0
−0.5
−1
−1.5
−2
−2.5
−3
1.84 2.34 4.45 7.86 8.07 13.56 15.59
Distance(km)
(b) Rf =25 Ω
2
EM1
1.5 EM2
EM3
1
0.5
Error(km)
0
−0.5
−1
−1.5
−2
−2.5
−3
1.84 2.34 4.45 7.86 8.07 13.56 15.59
Distance(km)
(c) Rf =50 Ω
Figure 6.54. Comparison of the distance estimate errors obtained from exact model
(EM1-EM3) algorithms with the effect of fault resistance (Rf ). Signal mea-
sured from MV/LV substation and fault inception angle is 45◦ . EM1 denotes
a distance calculation error using exact model 1, EM2 denotes a distance
calculation error using exact model 2 and EM3 denotes a distance calcula-
tion error using exact model 3.
161
Performance Evaluation of Proposed Fault Location Algorithms Using EMTP/ATP Simulation Model
In overall, the proposed algorithms are able to perform well with the
use of transient signal measured from the secondary side of the MV/LV
distribution transformer.
Considering the fault distance estimation results in this simulation,
NN algorithm shows a better performance compared with the other al-
gorithms. The lowest MAE calculated was 0.34 km and the highest was
1.18 km. The proposed fault location algorithms proved much more accu-
rate in this simulation model compared with previous simulation model
(model 3). We noticed that the number of branches in this simulation net-
work model is less than in the previous simulation network model. Also,
the length of the branches in simulation network model 4 is short com-
pared with the simulation network model 3. In this section, we did not
investigate the performance of WT algorithm because the transient sig-
nals that we used were measured from the secondary side of the MV/LV
distribution substation transformer.
It is worth to mentioned here that, based on our experience when pro-
cessing the measured transient signals, we found that it is more difficult
to analyze those measured transient signals when the fault distance is
far from the measuring MV/LV substation (DT24). In addition, with high
resistance values and low fault inception angles, the task is more critical.
This is because the damping will be increased and the transient signal is
damped much faster. In the following chapter, the proposed algorithms
are investigated using transient signals measured in real networks.
162
7. Performance Evaluation of Proposed
Fault Location Algorithms Using
Transients Recorded in Real
Networks
Three real case studies are presented in this chapter. These consist of iso-
lated, compensated and partially compensated neutral MV distribution
networks. Two of the cases were conducted in Finnish distribution net-
works and the third is performed in a Czech distribution network. The
recorded data is used to test selected fault location algorithms that we
have proposed in Chapter 5. Due to the limitation of information of line
parameters in real networks, only GM2 and CWT algorithm were tested
using real recorded data.
In this case, an artificial single line to earth fault was carried out in a
21 kV Finnish overhead line network. The case consists of unearthed and
compensated neutral MV distribution networks in Vaasan Sähkö Oy. The
data was recorded with sampling frequency of 20 kHz. In this case study,
CWT and GM2 algorithms were used to estimate the fault distance. The
procedures used to estimate the fault distance based on the measured
transient signals are explained in Chapter 5 of this dissertation. The per-
formance accuracy of the selected algorithms was measured using MAE
as given in Equation (6.1) and Equation (6.4).
The parameters of the networks used in this analysis are given in Ta-
ble 7.1. The network diagrams of the field tests are shown in Figures 7.1,
7.2 and 7.3.
163
Performance Evaluation of Proposed Fault Location Algorithms Using Transients Recorded in Real Networks
Figure 7.1. Compensated neutral network test diagram for test 1-2, 0.76 km [40].
164
Performance Evaluation of Proposed Fault Location Algorithms Using Transients Recorded in Real Networks
Figure 7.2. Compensated neutral network test diagram for test 3-4, 10.40 km [40].
165
Performance Evaluation of Proposed Fault Location Algorithms Using Transients Recorded in Real Networks
Figure 7.3. Isolated neutral network test diagram for test 5-8, 14.20 km [40].
166
Performance Evaluation of Proposed Fault Location Algorithms Using Transients Recorded in Real Networks
The recorded fault voltages and currents are shown in Figures 7.4, 7.5,
7.6 and 7.7. As an example of pre processing used in this work, one single
phase to ground fault from test number 1 is presented in Figure 7.8 and
the output from fault distances algorithm based on CWT is depicted in
Figure 7.9. The results of estimated fault distance using CWT algorithm
for test number 1, 3, 5 and 7 are depicted in Figure 7.10.
The summary of the results are given in Table 7.2. In order to compare
the performance accuracy of both algorithms, we used similar information
and parameters determined from the measured earth fault transient sig-
nals for both algorithms. Based on the results, CWT algorithm produced
good results compared to GM2 algorithm. GM2 algorithm performed well
in compensated neutral network compared to its result in isolated neutral
network. The MAE calculated was 0.52 km. In case of GM2 algorithm,
the MAE calculated was 3.80 km. GM2 algorithm gave good results when
the fault distance is short and close to the measuring point and when the
effect of fault resistance was low. Unfortunately, when fault distance is
far from the measuring point, the fault distance estimation error is high.
The result is of no use if the fault distance and fault resistance is to high.
Table 7.2. Calculation results using CWT and GM2 algorithm.
Test
l Rf CWT ∆CWT GM2 ∆GM2
number
1 0.76 0 0.38 -0.38 0.10 -0.66
2 0.76 0 0.17 -0.59 0.10 -0.66
3 10.40 0 9.92 -0.48 10.00 -0.40
4 10.40 0 9.81 -0.59 9.90 -0.50
5 14.20 0 13.57 -0.63 16.70 2.50
6 14.20 0 14.77 0.57 21.20 7.00
7 14.20 50 14.12 -0.08 20.63 6.43
8 14.20 50 13.39 -0.81 26.46 12.26
MAE 0.52 3.80
167
Performance Evaluation of Proposed Fault Location Algorithms Using Transients Recorded in Real Networks
4
x 10
2
U
1.5
Voltage [V] 1
0.5
0
−0.5
−1
−1.5
−2
0 200 400 600 800 1000 1200 1400 1600 1800 2000
Sample Number [k]
1500
I
1000
Current [A]
500
−500
−1000
0 200 400 600 800 1000 1200 1400 1600 1800 2000
Sample Number [k]
0.5
0
−0.5
−1
−1.5
−2
0 200 400 600 800 1000 1200 1400 1600 1800 2000
Sample Number [k]
1500
I
1000
Current [A]
500
−500
0 200 400 600 800 1000 1200 1400 1600 1800 2000
Sample Number [k]
Figure 7.4. Recorded signals of single line to earth fault at 0.46 km.
168
Performance Evaluation of Proposed Fault Location Algorithms Using Transients Recorded in Real Networks
4
x 10
2
U
1.5
1
Voltage [V]
0.5
0
−0.5
−1
−1.5
−2
0 200 400 600 800 1000 1200 1400 1600 1800 2000
Sample Number [k]
300
I
200
Current [A]
100
−100
−200
−300
0 200 400 600 800 1000 1200 1400 1600 1800 2000
Sample Number [k]
0.5
0
−0.5
−1
−1.5
−2
0 200 400 600 800 1000 1200 1400 1600 1800 2000
Sample Number [k]
300
I
200
Current [A]
100
−100
−200
−300
0 200 400 600 800 1000 1200 1400 1600 1800 2000
Sample Number [k]
169
Performance Evaluation of Proposed Fault Location Algorithms Using Transients Recorded in Real Networks
4
x 10
2
V
1.5
Voltage [V] 1
0.5
0
−0.5
−1
−1.5
−2
0 200 400 600 800 1000 1200 1400 1600 1800 2000
Sample Number [k]
100
I
50
Current [A]
−50
−100
−150
0 200 400 600 800 1000 1200 1400 1600 1800 2000
Sample Number [k]
1
Voltage [V]
0.5
−0.5
−1
−1.5
−2
0 200 400 600 800 1000 1200 1400 1600 1800 2000
Sample Number [k]
150
I
100
Current [A]
50
−50
−100
0 200 400 600 800 1000 1200 1400 1600 1800 2000
Sample Number [k]
Figure 7.6. Recorded signals of single line to earth fault at 14.20 km. Test number 5 and
6.
170
Performance Evaluation of Proposed Fault Location Algorithms Using Transients Recorded in Real Networks
4
x 10
2
V
1.5
1
Voltage [V]
0.5
0
−0.5
−1
−1.5
−2
0 200 400 600 800 1000 1200 1400 1600 1800 2000
Sample Number [k]
80
I
60
40
Current [A]
20
0
−20
−40
−60
−80
0 200 400 600 800 1000 1200 1400 1600 1800 2000
Sample Number [k]
0.5
0
−0.5
−1
−1.5
−2
0 200 400 600 800 1000 1200 1400 1600 1800 2000
Sample Number [k]
80
I
60
40
Current [A]
20
−20
−40
−60
−80
0 200 400 600 800 1000 1200 1400 1600 1800 2000
Sample Number [k]
Figure 7.7. Recorded signals of single line to earth fault at 14.20 km. Test number 7 and
8.
171
Performance Evaluation of Proposed Fault Location Algorithms Using Transients Recorded in Real Networks
80 80
60
60 70
50
40 60
Scales[a]
20 40 50
X= 20
80 Y= 33 X: 33 40
Level= 0.11133 30 Y: 3.954
60 30
20
40 20
10
20 10
80 0.1
60
0.09
40
60 0.08
20
0.07
1
Scales[a]
0.06
40
80
0.05
X= 34 X: 34
Y= 34 Y: 3.918 60
0.04
Level= 0.10308
20
0.03 40
0.02 20
1 0.016
0 50 100 150 200 5 4 3 2 1 01
Sample Number[k] Energy percentage[%]
−1
−2
0 20 40 60 80 100 120 140 160 180 200
Sample Number [k]
4000
I
CWT
2000
Current[A]
−2000
−4000
−6000
0 20 40 60 80 100 120 140 160 180 200
Sample Number [k]
3000
Magnitude [FFT]
2500
2000
1500
X: 468.8
Y: 929.6
1000
500
0
0 500 1000 1500 2000 2500 3000 3500 4000
Frequency [Hz]
Figure 7.8. Example of output from pre-processing of test number 1 (0.76 km) transient
signal using CWT analysis.
172
Performance Evaluation of Proposed Fault Location Algorithms Using Transients Recorded in Real Networks
4
x 10
2
U [V]
0
v
−2
0 20 40 60 80 100 120 140 160 180 200
5000
U [A]
0
I
−5000
0 20 40 60 80 100 120 140 160 180 200
−3
x 10
4
L(H)
2
0
0 20 40 60 80 100 120 140 160 180 200
0.4
X: 86
xstd
0.2 Y: 0.06533
0
0 20 40 60 80 100 120 140 160 180 200
2
lf(km)
X: 86
1 Y: 0.3747
0
0 20 40 60 80 100 120 140 160 180 200
Sample Number [k]
Figure 7.9. Example of output from fault distances algorithm based on CWT for test
number 1.
100
80
l [km]
60
f
X: 51
40 Y: 30.34
20
0 20 40 60 80 100 120 140 160 180 200
Sample Number [k]
100
80
lf[km]
60
40 X: 14
Y: 25.78
20
0 20 40 60 80 100 120 140 160 180 200
Sample Number [k]
60
40
X: 13
Y: 13.57
lf[km]
20
−20
−40
0 20 40 60 80 100 120 140 160 180 200
Sample Number [k]
100
50
l [km]
0 X: 11
f
Y: 14.12
−50
−100
0 20 40 60 80 100 120 140 160 180 200
Sample Number [k]
Figure 7.10. Results of fault distance estimation (lf ) using CWT algorithm.
173
Performance Evaluation of Proposed Fault Location Algorithms Using Transients Recorded in Real Networks
In this section an artificial single line to earth fault was carried out in
a 22 kV Finnish overhead line network. The field test was recorded in
Lounais-Suomen Sähkö Oy. The network is partially compensated neu-
tral network. The data was recorded with sampling frequency of 20 kHz.
In this case study, CWT and GM2 algorithms were used to estimate the
fault distance. The procedure used to estimate the fault distance was
based on the measured transient signal as is explained in Chapter 5 of
this dissertation. The performance accuracy of the selected algorithms
were measured using MAE as given in Equation (6.1) and Equation (6.4).
The parameters of the network are given in Table 7.3. The network dia-
grams of the field test are shown in Figure 7.11.
The recorded fault voltages and currents are shown in Figures 7.12, 7.13,
7.14 and 7.15. Figure 7.16 and 7.17 shows an example of output from
pre-processing of test number 4 transient signal using CWT analysis and
Hilbert transformation. The results of estimated fault distance using
CWT algorithm for test number 1, 3, 6 and 8 is depicted in Figure 7.18.
The summary of the results are given in Table 7.4. By comparison of MAE
for both algorithms, CWT algorithm performed slightly better compared
to GM2 algorithm. There are 2 results showing unacceptable accuracy
performed by GM2 algorithm where MAE increased significantly. The
lowest fault estimation error was 0.1 km and the highest was 12.40 km.
In this case GM2 algorithm gave good results when the fault distance
is far from the measuring point. In comparison with the result of GM2
algorithm in previous case study for fault distance far from the measur-
ing point, it seems that the algorithm has performed well in compensated
neutral network compared to its result in isolated neutral network.
174
Performance Evaluation of Proposed Fault Location Algorithms Using Transients Recorded in Real Networks
4
x 10
4
U
3
2
Voltage [V]
1
0
−1
−2
−3
−4
0 200 400 600 800 1000 1200 1400 1600 1800 2000
Sample Number [k]
150
I
100
Current [A]
50
−50
−100
0 200 400 600 800 1000 1200 1400 1600 1800 2000
Sample Number [k]
4
x 10
4
U
3
2
Voltage [V]
1
0
−1
−2
−3
−4
0 200 400 600 800 1000 1200 1400 1600 1800 2000
Sample Number [k]
150
I
100
Current [A]
50
−50
−100
−150
0 500 1000 1500 2000 2500 3000 3500 4000
Sample Number [k]
Figure 7.12. Recorded signals of single line to earth fault at 25.40 km. Test number 1
and 2.
176
Performance Evaluation of Proposed Fault Location Algorithms Using Transients Recorded in Real Networks
4
x 10
4
U
3
2
Voltage [V]
1
0
−1
−2
−3
−4
0 500 1000 1500 2000
Sample Number [k]
150
I
100
50
Current [A]
0
−50
−100
−150
−200
−250
0 500 1000 1500 2000
Sample Number [k]
1
0
−1
−2
−3
−4
0 500 1000 1500 2000
Sample Number [k]
200
I
100
Current [A]
−100
−200
−300
0 500 1000 1500 2000
Sample Number [k]
Figure 7.13. Recorded signals of single line to earth fault at 25.40 km. Test number 3
and 4.
177
Performance Evaluation of Proposed Fault Location Algorithms Using Transients Recorded in Real Networks
4
x 10
4
U
3
2
Voltage [V]
1
0
−1
−2
−3
−4
0 500 1000 1500 2000
Sample Number [k]
300
I
250
200
Current [A]
150
100
50
0
−50
−100
−150
0 200 400 600 800 1000 1200 1400 1600 1800 2000
Sample Number [k]
1
0
−1
−2
−3
−4
0 200 400 600 800 1000 1200 1400 1600 1800 2000
Sample Number [k]
400
I
300
Current [A]
200
100
−100
−200
0 500 1000 1500 2000
Sample Number [k]
Figure 7.14. Recorded signals of single line to earth fault at 25.40 km (test number 5)
and 36.0 km (test number 6).
178
Performance Evaluation of Proposed Fault Location Algorithms Using Transients Recorded in Real Networks
4
x 10
4
U
3
2
Voltage [V]
1
0
−1
−2
−3
−4
0 500 1000 1500 2000 2500 3000 3500 4000
Sample Number [k]
300
I
200
Current [A]
100
−100
−200
−300
−400
0 500 1000 1500 2000
Sample Number [k]
1
0
−1
−2
−3
−4
0 500 1000 1500 2000 2500 3000
Sample Number [k]
400
I
300
Current [A]
200
100
−100
−200
0 500 1000 1500 2000 2500 3000
Sample Number [k]
Figure 7.15. Recorded signals of single line to earth fault at 36.0 km.
179
Performance Evaluation of Proposed Fault Location Algorithms Using Transients Recorded in Real Networks
4
x 10
2
Voltage[V}
V
0
−2
−4
20 40 60 80 100 120 140 160 180 200
115 0.03
109
103
97 0.025
91 X= 51
85 Y= 72
79 Level= 0.029502
Scales a
73 0.02
67
61
55
49 0.015
43
37
31
25 0.01
19
13
7
1 0.005
20 40 60 80 100 120 140 160 180 200
Sample Number[k]
2.5 X: 72
Y: 2.308
Energy percentage[%]
1.5
0.5
0
0 20 40 60 80 100 120
Scales[a]
200
I
Current[A]
100
0
−100
−200
−300
20 40 60 80 100 120 140 160 180 200
0.03
115
109
103
97 0.025
91
85
79
Scales[a]
73 0.02
67 X= 72
61 Y= 76
55 Level= 0.029556
49 0.015
43
37
31
25 0.01
19
13
7
1 0.005
20 40 60 80 100 120 140 160 180 200
Sample Number[k]
X: 76
Y: 2.1
2
Energy percentage[%]
1.5
0.5
0
0 20 40 60 80 100 120
Scales[a]
Figure 7.16. Example of CWT scalogram and energy percentage distribution of test num-
ber 4 (25.4 km) transient signal.
180
Performance Evaluation of Proposed Fault Location Algorithms Using Transients Recorded in Real Networks
5
x 10
1
V
CWT
0.5
Voltage[V]
0
−0.5
−1
0 20 40 60 80 100 120 140 160 180 200
Sample Number[k]
1500
I
CWT
1000
Current[A]
500
0
−500
−1000
−1500
0 20 40 60 80 100 120 140 160 180 200
Sample Number [k]
3.5
3
Magnitude[FFT]
2.5
1.5
0.5 X: 234.4
Y: 451.9
0
0 1000 2000 3000 4000 5000 6000
Frequency[Hz]
11.2
11.1
log(|V (t)|)
11
WT
10.9
10.8
10.7
10.6
10.5
10.4
100 110 120 130 140 150 160 170 180 190 200
Sample Number[k]
6.8
I
WT
6.7
6.6
6.5
log(|IWT(t)|
6.4
6.3
6.2
6.1
5.9
100 110 120 130 140 150 160 170 180 190 200
Sample Number[k]
Figure 7.17. Example of CWT extracted coefficients, FFT frequency spectrum of ec-
tracted coefficent and, voltage and current linear decaying of extracted CWT
coefficient of test number 4 (25.4 km) transient signal.
181
Performance Evaluation of Proposed Fault Location Algorithms Using Transients Recorded in Real Networks
100
80
lf[km]
60
X: 51
40 Y: 30.34
20
0 20 40 60 80 100 120 140 160 180 200
Sample Number [k]
100
80
lf[km]
60
40 X: 14
Y: 25.78
20
0 20 40 60 80 100 120 140 160 180 200
Sample Number [k]
100
80
lf[km]
60
X: 12
Y: 39.53
40
20
0 20 40 60 80 100 120 140 160 180 200
Sample Number [k]
100
80
lf[km]
60 X: 12
Y: 38.46
40
20
0 20 40 60 80 100 120 140 160 180 200
Sample Number [k]
Test
l Rf CWT ∆CWT GM2 ∆GM2
number
1 25.4 0 30.34 4.94 25.30 -0.10
2 25.4 0 29.82 4.42 25.30 -0.10
3 25.4 0 25.78 0.38 36.70 11.30
4 25.4 0 21.02 -4.38 25.30 -0.10
5 25.4 0 24.19 -1.21 37.80 12.40
6 36.0 0 39.53 3.53 36.70 0.70
7 36.0 0 38.70 2.70 36.10 0.10
8 36.0 0 38.46 2.46 35.50 -0.50
MAE 3.00 3.16
182
Performance Evaluation of Proposed Fault Location Algorithms Using Transients Recorded in Real Networks
In this section an artificial single line to earth fault was carried out in a
22 kV Czech compensated neutral network. The field test was recorded
by David Topolanek and his team from Brno University. The data was
used during his research visit in autumn 2012 at Aalto University. The
measured signal was recorded from MV substation and secondary side of
MV/LV distribution transformer. Due to the limitation of information, the
data was tested only with CWT algorithm described in Section 5.3.
OTHER FEEDERS
110/22kV .
.
.
.
.
.
.
.
.
.
Rp Lc
912m 695m
22/0.4kV 22/0.4kV 8.5m 22/0.4kV 15m 22/0.4kV
Legend:
Measuring point
Fault position
Domanin Olesna Nove Mesto Pohledec
183
Performance Evaluation of Proposed Fault Location Algorithms Using Transients Recorded in Real Networks
(a)
(b)
Figure 7.20. (a) Photo of one of the low voltage (LV) switch-boards and the digital oscil-
loscope used during the field test. (b) Connection of measuring probe at the
back of the LV distribution panel
184
Performance Evaluation of Proposed Fault Location Algorithms Using Transients Recorded in Real Networks
The recorded fault voltages and currents are shown in Figures 7.21, 7.22,
7.23, 7.24 and 7.25. The fault location accuracy is measured with mean
average error (MAE) in km as in Equation (6.4). Table 7.6 shows the
results of the proposed scheme tested on the real recorded data. In overall,
the MAE for fault distance is 0.65 km. The fault distance for test number
1 , 3, 4 and 5 calculated using CWT algorithm is depicted in Figure 7.26.
The summary of the results are given in Table 7.6.
In Table 7.6, the results show that by using the proposed algorithm we
managed to identify the correct LV substation to be used as indicator for
selecting a correct path towards the fault position. Based on the results,
the detected LV substation is located in front of the fault position. The dis-
tance between the detected LV substation and the fault location is about
118.5 meters, which is very close. The results from test number 1 to test
number 4 show that substation “Pohledec” record the second highest peak
amplitude. The substation is located 3.148 km behind the fault position.
In this experiment result, it seems that the first and second highest peak
amplitude of the LV substations gives us an indication that the fault is
located in the section between these two substations.
185
Performance Evaluation of Proposed Fault Location Algorithms Using Transients Recorded in Real Networks
4
x 10
4
Ua 400 Ia
3
Ub Ib
2 200 Ic
Current [A]
Voltage [V]
Uc
1 0
0
−200
−1
−2 −400
−3
−600
−4
1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2 1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2
Sample Number [k] x 10
4 Sample Number [k] x 10
4
Voltage [V]
Voltage [V]
200 200
100 100
0 0
−100 −100
−200 −200
−300 −300
−400 −400
−500 −500
1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2 1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2
Sample Number [k] x 10
4
Sample Number [k] x 10
4
200 200
100 100
0 0
−100 −100
−200 −200
−300 −300
−400 −400
−500 −500
1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2 1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2
4 4
Sample Number [k] x 10 Sample Number [k] x 10
4
x 10
5 400 Ia
4 Ua Ib
3 Ub 200 Ic
Current [A]
Uc
Voltage [V]
2
0
1
0 −200
−1
−400
−2
−3 −600
−4 1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2
1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2
4
Sample Number [k] x 10
4 Sample Number [k] x 10
Voltage [V]
200 200
100 100
0 0
−100 −100
−200 −200
−300 −300
−400 −400
−500 −500
1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2 1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2
4
Sample Number [k] x 10 Sample Number [k] x 10
4
Uc
Voltage [V]
100 200
0 100
0
−100
−100
−200
−200
−300
−300
−400
−400
−500 −500
1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2 1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2
4
Sample Number [k] x 10 Sample Number [k] 4
x 10
186
Performance Evaluation of Proposed Fault Location Algorithms Using Transients Recorded in Real Networks
4 400
x 10
5
Ia
4 Ua Ib
200
Ub Ic
3 Uc
Current [A]
Voltage [V]
2 0
1
−200
0
−1 −400
−2
−3 −600
−4 1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2
1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2 4
Sample Number [k] x 10
4 Sample Number [k] x 10
Voltage [V]
Voltage [V]
100 100
0 0
−100 −100
−200 −200
−300 −300
−400 −400
−500 −500
1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2 1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2
Sample Number [k] x 10
4 Sample Number [k] x 10
4
Voltage [V]
200
Voltage [V]
100 100
0 0
−100 −100
−200 −200
−300
−300
−400
−400
−500
−500 1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2
1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2 4
Sample Number [k] 4 Sample Number [k] x 10
x 10
4
x 10 600 Ia
4 Ua Ib
Ub 400 Ic
Current [A]
Voltage [V]
2 Uc
200
0
0
−2
−200
−4
−400
−6 1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2
1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2 x 10
4
Sample Number [k] 4 Sample Number [k]
x 10
Uc Uc
Voltage [V]
200 200
100 100
0 0
−100 −100
−200 −200
−300 −300
−400 −400
−500 −500
1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2 1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2
4 4
Sample Number [k] x 10 Sample Number [k] x 10
Voltage [V]
200 200
100 100
0 0
−100 −100
−200 −200
−300 −300
−400 −400
−500 −500
1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2 1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2
4
Sample Number [k] x 10 Sample Number [k] x 10
4
187
Performance Evaluation of Proposed Fault Location Algorithms Using Transients Recorded in Real Networks
4
x 10 600
4 Ia
3 Ua Ib
Ub 400 Ic
Current [A]
2
Voltage [V]
Uc
200
1
0 0
−1
−200
−2
−3 −400
−4 1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2
1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2 4
Sample Number [k] 4 Sample Number [k] x 10
x 10
200 Uc
Voltage [V]
200
100 100
0 0
−100 −100
−200 −200
−300 −300
−400 −400
−500 −500
1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2 1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2
4
Sample Number [k] x 10 4
Sample Number [k] x 10
200 200
Voltage [V]
100 100
0 0
−100 −100
−200 −200
−300 −300
−400 −400
−500 −500
1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2 1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2
4
Sample Number [k] x 10 Sample Number [k] x 10
4
60
lf[km]
40
X: 172
20 Y: 14.04
0
0 50 100 150 200 250
Sample Number[k]
40
lf[km]
X: 166
20 Y: 14.44
0
0 50 100 150 200 250
Sample Number[k]
40
lf[km]
X: 232
20 Y: 15.6
0
0 50 100 150 200 250
Sample Number[k]
50
40
l [km]
30 X: 164
20 Y: 14.52
f
10
0
0 50 100 150 200 250
Sample Number[k]
188
Performance Evaluation of Proposed Fault Location Algorithms Using Transients Recorded in Real Networks
Table 7.6. Test result of single line to earth fault using recorded signal.
Test
1 2 3 4 5
number
Fault
14.01 14.15 14.44 15.6 14.52
Location (km)
189
Performance Evaluation of Proposed Fault Location Algorithms Using Transients Recorded in Real Networks
190
8. Conclusions and Future Work
This thesis presented fault location algorithms for earth faults in un-
earthed and compensated neutral networks using measured transient sig-
nals. The initial transients of the earth faults are important especially in
unearthed and compensated neutral networks. In this work, five types
of fault location algorithms based on initial transients of the earth faults
have been developed. Two of the algorithms were developed using the
measured transient signals from secondary side of the MV/LV distribu-
tion transformers. Five types of fault location algorithms developed in
this thesis are known as general model (GM), exact model (EM), contin-
uous wavelet transform (CWT), multiple regression analysis (MRA) and
artificial neural network (ANN) algorithms.
Considering the results from the simulated network model 1, we can
conclude that the damped and undamped charge transient frequency and
damping factor from the simulation agreed with the result calculated
by theoretical formula of the general model algorithms, with acceptable
small error. Based on the fault distance result, it was shown that GM1
and GM2 algorithms have almost identical result. It shows that, in order
to estimate the fault distance, by using the undamped charge transient
frequency, GM2 algorithm managed to compensate the effect of fault re-
sistance as well as the unknown resistive parameters of the circuit.
We compared the performance of EM1 and GM2 algorithms using se-
quence component network in simulation model 2. In overall, we conclude
that EM1 algorithm has performed much better than GM2 algorithm. The
main reason why EM1 algorithm performs much better than GM2 is that
the EM1 circuit model has a similar structure as the simulation network,
except that we purposely neglected the resistive part in EM1 model. On
the other hand, GM2 circuit model is more robust and simple. Both al-
gorithms used same signal pre processing technique to estimate the un-
191
Conclusions and Future Work
192
Conclusions and Future Work
193
Conclusions and Future Work
194
References
[2] K. Winter, “Null point analysis - new method for detecting high resistance
earth faults,” Elektricitetens Rationella Användning (ERA), vol. 61, pp. 18–
20,23–24, 1988.
[5] G. Eberl and P. Schegner, “Earth fault location based on short time
low ohmic neutral earthing in cabled medium voltage networks of the
DREWAG - Stadtwerke Dresden GmbH,” in Seminar on Methods and
Techniques for Earth Fault Detection, Indication and Location. Espoo
Finland, Feb 2011, pp. 135–141.
[7] M. Lehtonen and T. Hakola, Neutral earthing and power system protection:
earthing solutions and protective relaying in medium voltage distribution
networks. ABB Transmit Oy, 1996.
[8] O. Chaari, M. Meunier, and F. Brouaye, “Wavelets: a new tool for the res-
onant grounded power distribution systems relaying,” IEEE Transactions
on Power Delivery, vol. 11, no. 3, pp. 1301–1308, July 1996.
[9] A. Nikander and P. Järventausta, “Safety aspects and novel technical solu-
tions for earth fault management in MV electricity distribution networks,”
in IET Conference Proceedings of 3rd IEE International Conference on Re-
liability of Transmission and Distribution Networks (RTDN 2005), vol. 4.
London UK, Jan 2005, pp. 207–211.
195
References
[17] A. Girgis, C. Fallon, and D. Lubkeman, “A fault location technique for rural
distribution feeders,” IEEE Transactions on Industry Applications, vol. 29,
no. 6, pp. 1170–1175, Nov 1993.
[20] P. Imriš and M. Lehtonen, “Transient based ground fault location using
wavelets,” in 4th IASTED International Conference on Power and Energy
Systems. Rhodes, Greece, 28-30 June 2004, pp. 507–511.
[22] F. H. Magnago and A. Abur, “A new fault location technique for radial
distribution systems based on high frequency signals,” in IEEE Power En-
gineering Society Summer Meeting, vol. 1. Edmonton, Alberta, Canada,
18-22 July 1999, pp. 426–431.
196
References
[30] Y.-Y. Hsu et al., “An expert system for locating distribution system faults,”
IEEE Transactions on Power Delivery, vol. 6, no. 1, pp. 366–372, Jan 1991.
[31] P. Järventausta, P. Verho, and J. Partanen, “Using fuzzy sets to model the
uncertainty in the fault location process of distribution networks,” IEEE
Transactions on Power Delivery, vol. 9, no. 2, pp. 954–960, Apr 1994.
[32] Y. D. Xue, B. Y. Xu, and Z. H. Wang, “The fault location technology us-
ing transient signals for single phase earth fault in non-solidly earthed
network,” in 22nd International Conference and Exhibition on Electricity
Distribution (CIRED 2013). Stockholm Sweden, June 2013, pp. 1–4.
[33] W.-R. Chen, Q.-Q. Qian, and X.-R. Wang, “Wavelet neural network based
transient fault signal detection and identification,” in Proceedings of 1997
International Conference on Information, Communications and Signal
Processing (ICICS), vol. 3. Singapore, 9-12 Sept 1997, pp. 1377–1381.
[34] X.-H. Zhang, S.-X. Fan, H. Duan, and B.-Y. Xu, “Analysis of transient dom-
inant frequency signal for single-phase earthed fault based on Prony algo-
rithm,” in China International Conference on Electricity Distribution (CI-
CED). Guangzhou China, 10-13 Dec 2008, pp. 1–6.
197
References
[35] S. Huang and X. G. Wang, “A fault location scheme based on spectrum char-
acteristic of fault-generated high-frequency transient signals,” in IEEE
Power Energy Society General Meeting (PES ’09). Calgary Canada, July
2009, pp. 1–5.
[37] G. Druml, A. Kugi, and O. Seifert, “A new directional transient relay for
high ohmic earth faults,” in CIRED 17th International Conference on Elec-
tricity Distribution. Barcelona Spain, 12-15 May 2003, pp. 1–6.
[41] P. Schegner and G. Eberl, “Overview of methods for earth fault distance
calculation,” in Seminar on Methods and Techniques for Earth Fault De-
tection, Indication and Location. Espoo Finland, Feb 2011, pp. 36–49.
[43] S. Hänninen, Single phase earth faults in high impedance grounded net-
works: characteristics, indication and location, ser. VTT Publications 453.
Espoo: Espoo : Technical Research Centre of Finland, 2001.
[45] A. G. Phadke and J. S. Thorp, Computer relaying for power systems, 2nd ed.
John Wiley & Sons Ltd, England, 2009.
198
References
[52] M. Aurangzeb, P. Crossley, and P. Gale, “Fault location using the high fre-
quency travelling waves measured at a single location on a transmission
line,” in Seventh International Conference on Developments in Power Sys-
tem Protection. Amsterdam Netherlands, Apr 2001, pp. 403–406.
[60] S. Hänninen and M. Lehtonen, Earth fault distance computation with fun-
damental frequency signals based on measurements in substation supply
bay, ser. VTT tiedotteita - Research Notes 2153. Espoo : Technical Re-
search Centre of Finland, 2002.
199
References
[63] C.-H. Kim et al., “A study on the fault indentification of underground cable
using neural networks,” in Proceedings of International Conference on En-
ergy Management and Power Delivery (EMPD ’95), vol. 2. Singapore, Nov
1995, pp. 571–576.
[65] Y. Assef, P. Bastard, and M. Meunier, “Artificial neural networks for single
phase fault detection in resonant grounded power distribution systems,”
in Proceedings of Transmission and Distribution Conference. Los Angeles
USA, Sep 1996, pp. 566–572.
[68] Z. Yujun and W. Xiaowei, “Research on wavelet neural network for fault lo-
cation in power distribution network,” in Electrical, Information Engineer-
ing and Mechatronics 2011, ser. Lecture Notes in Electrical Engineering,
X. Wang, F. Wang, and S. Zhong, Eds. Springer London, 2012, vol. 138,
pp. 1049–1057.
200
References
[73] R. Salat and S. Osowski, “Accurate fault location in the power transmission
line using support vector machine approach,” IEEE Transactions on Power
Systems, vol. 19, no. 2, pp. 979–986, May 2004.
[77] P. Dash, A. Pradhan, and G. Panda, “A novel fuzzy neural network based
distance relaying scheme,” IEEE Transactions on Power Delivery, vol. 15,
no. 3, pp. 902–907, July 2000.
[80] T. ElAli, Discrete Systems and Digital Signal Processing with MATLAB,
Second Edition, 2nd ed. Boca Raton : CRC Press, Dec 2011.
[83] W. Cochran et al., “What is the fast Fourier transform?” Proceedings of the
IEEE, vol. 55, no. 10, pp. 1664–1674, Oct 1967.
[84] G. Heydt, “A new method for the calculation of subtransmission and distri-
bution system transients based on the FFT,” IEEE Transactions on Power
Delivery, vol. 4, no. 3, pp. 1869–1875, July 1989.
201
References
[92] M. Gilany, E. El Din, M. Abdel Aziz, and D. Ibrahim, “An accurate scheme
for fault location in combined overhead line with underground power ca-
ble,” in IEEE Power Engineering Society General Meeting, vol. 3. San
Francisco USA, June 2005, pp. 2521–2527.
[98] L. Satish, “Short-time Fourier and wavelet transforms for fault detection
in power transformers during impulse tests,” IEE Proceedings - Science,
Measurement and Technology, vol. 145, no. 2, pp. 77–84, Mar 1998.
202
References
[99] B. Xu, H. Wang, L. Sun, and F. Yang, “Detection methods of broken rotor
bar fault in squirrel cage induction motor with light-load,” in International
Conference on Electrical Machines and Systems ( ICEMS 2008). Wuhan
China, Oct 2008, pp. 759–762.
[100] J.-Y. Lee, “Variable short-time Fourier transform for vibration signals with
transients,” Journal of Vibration and Control, vol. 1, pp. 1–15, July 2013.
[104] J. Rosero, L. Romeral, J. Ortega, and E. Rosero, “Short circuit fault de-
tection in PMSM by means of Empirical Mode Decomposition (EMD) and
Wigner Ville Distribution (WVD),” in IEEE Twenty-Third Annual Applied
Power Electronics Conference and Exposition (APEC 2008). Austin Texas
USA, Feb 2008, pp. 98–103.
[107] M. Prieto, A. Espinosa, J.-R. Ruiz, J. Urresty, and J. Ortega, “Feature ex-
traction of demagnetization faults in permanent-magnet synchronous mo-
tors based on box-counting fractal dimension,” IEEE Transactions on In-
dustrial Electronics, vol. 58, no. 5, pp. 1594–1605, May 2011.
203
References
[110] A. Quinquis and C. Dugal, “Using the wavelet transform for the detec-
tion of magnetic underwater transient signals,” in Proceedings Oceans En-
gineering for Today’s Technology and Tomorrow’s Preservation (OCEANS
’94), vol. 2. Brest, France, Sep 1994, pp. 538 – 543.
[114] S. Chen and H. Y. Zhu, “Wavelet transform for processing power quality
disturbances,” EURASIP Journal on Advances in Signal Processing, vol.
2007, pp. 1–20, Feb 2007.
[121] F. H. Magnago and A. Abur, “Fault location using wavelets,” IEEE Trans-
actions on Power Delivery, vol. 13, no. 4, pp. 1475–1480, Oct 1998.
204
References
[126] M. Misiti, Y. Misiti, G. Oppenheim, and J.-M. Poggi, Wavelet toolbox for use
with MATLAB, 1st ed. Natick, Mass. : Math Works, 1997.
[127] R. Bracewell, The Fourier Transform and Its Applications, 3rd ed., ser.
McGraw-Hill international editions. New Delhi : Tata McGraw-Hill, 2000.
[131] A. Agarwal and J. Lang, Foundations of Analog and Digital Electronic Cir-
cuits, 3rd ed., ser. The Morgan Kaufmann Series in Computer Architecture
and Design. Amsterdam ; Boston : Elsevier : Morgan Kaufman Publish-
ers, 2005.
[134] A. Chuang. (2008, October) Epri smart grid R&D overview. Pre-
sentation at HP sustainability innovation workshop. EPRI. [Online].
Available: http://www.hpl.hp.com/open_innovation/workshops/microgrid%
20and%20distributed%20power/HP_Chuang_102008_post.pdf
[136] D. Topolanek and P. Toman, “Possible ground fault location using data
recorded on the secondary side of distribution transformers,” in The 12th
international Scientific Conference on Electric Power Engineering (EPE12).
Kouty nad Desnou, Czech Republic, May 2011, pp. 209–212.
205
References
206
Appendices
A. Mathematical Derivations and
MATLAB Scripts
R2 L2 C2
L1
E i2 i1 C1
R1
Figure A1. The equivalent circuit for ground fault transient model 1
The equations for voltages around the loops of the circuit are in Laplace-
domain as follows
1
i1 (s)(R1 + SL1 + ) − i2 (s)(R1 + SL1 ) = 0 (1)
SC1
1 E
−i1 (s)(R1 + SL1 ) + i2 (s)(R1 + R2 + SL1 + SL2 + )= (2)
SC2 S
The current i2 (s) can be solved by rearranging i1 (s) from equation (1) and
substituting it into equation (2). After the substitutions and rearranging
the terms, equation (2) becomes:
where,
209
Mathematical Derivations and MATLAB Scripts
and
Q(s) = S 4 L2 C2 L1 C1 + S 3 (R2 C2 L1 C1 + L2 C2 R1 C1 )
+S 2 (R2 C2 R1 C1 + L1 C2 + L1 C1 + L2 C2 )
+S(R1 C2 + R2 C2 + R1 C1 ) + 1 (5)
L2 C2
E i2 i1 C1
L1
Figure A2. The equivalent circuit for ground fault transient model 2
The equations for voltages around the loops of the circuit are in Laplace-
domain as follows
1
i1 (s)(SL1 + ) − i2 (s)(SL1 ) = 0 (1)
SC1
1 E
−i1 (s)SL1 + i2 (s)(SL1 + SL2 + )= (2)
SC2 S
The current i2 (s) can be solved by rearranging i1 (s) from equation (1)
and substituting it into equation (2). After the substitutions and rear-
ranging the terms, equation (2) becomes,
where,
and
210
Mathematical Derivations and MATLAB Scripts
R2 L2 C2 R3 L3
L1
E i2 i1 C1 i3 L4 i4 R4
R1
Figure A3. The equivalent circuit for ground fault transient model 3
The equations for voltages around the loops of the circuit are in Laplace-
domain as follows
1 1
i1 (s)(R1 + SL1 + ) − i2 (s)(R1 + SL1 ) − i3 (s) =0 (1)
SC1 SC1
1 E
−i1 (s)(R1 + SL1 ) + i2 (s)(R2 + SL1 + SL2 + )= (2)
SC2 S
1 1
−i1 (s) + i3 (s)(R3 + SL3 + SL4 + ) − i4 (s)SL4 = 0 (3)
SC1 SC1
−i3 (s)SL4 + i4 (s)(R4 + SL4 ) = 0 (4)
From which the current i2 (s) can be solved in the following steps, First
Solving i4 (s) from equation (4) and substituting it into equation (3). Then,
solving i3 (s) from equation (3) and substituting it into equation (1). Next,
solving i1 (s) from equation (4) and substituting it into equation (2). After
the substitutions and rearranging the terms, equation (2) becomes:
where,
+ R1 L3 C1 L4 + L1 L4 C1 R4 ) + S 2 (L1 R3 C1 R4 + R1 L4 C1 R4 + L4 L1
+ R1 R3 C1 L4 + L3 L4 + R1 L3 C1 R4 ) + S(R3 L4 + R1 R3 C1 R4 + L1 R4
+ L4 R1 + L4 R4 + L3 R4 ) + R1 R4 + R3 R4 ) (6)
211
Mathematical Derivations and MATLAB Scripts
and,
+ C1 L4 L3 R1 L2 C2 + R4 L4 L1 C1 L2 C2 + C1 L4 R3 L1 L2 C2 )+
S 4 (C1 L4 R3 R1 L2 C2 + L4 L1 L2 C2 + C1 L4 R3 L1 R2 C2
+ C 1 L 4 L 3 R1 R 2 C 2 + R4 L 4 R 1 C 1 L 2 C 2 + R 4 L 4 L 1 C 1 R 2 C 2
+ C 1 L 3 R4 L 1 R 2 C 2 + L 4 L 3 L 2 C 2 + C 1 R3 R 4 L 1 L 2 C 2
+ C1 L3 R4 R1 L2 C2 + L3 L4 L1 C2 + L1 L3 C1 L4 )+
S 3 (R4 L4 R1 C1 R2 C2 + C1 L3 R4 R1 R2 C2 + C1 R3 R4 R1 L2 C2
+ C 1 L 4 R 3 R 1 R2 C 2 + C 1 R 3 R 4 L 1 R 2 C 2 + L 4 R 3 L 2 C 2
+ L 4 R 4 L 2 C 2 + L 4 L 1 R2 C 2 + L 4 L 3 R2 C 2 + R 4 L 1 L 2 C 2
+ L 3 R 4 L 2 C 2 + L 4 R 1 L 2 C 2 + L 1 L 3 C 1 R4 + L 1 L 4 C 1 R4
+ L 1 R3 C 1 L 4 + L 3 R4 L 1 C 2 + L 4 R4 L 1 C 2 + R 3 L 4 L 1 C 2
+ R1 L3 C1 L4 + L3 L4 R1 C2 ) + S 2 (L4 R4 R2 C2 + L4 R3 R2 C2
+ L 4 R1 R 2 C 2 + R4 L 1 R 2 C 2 + R3 R 4 L 2 C 2 + L 3 R 4 R2 C 2 + L 1 L 4
+ R4 R 1 L 2 C 2 + C 1 R 3 R 4 R 1 R2 C 2 + L 3 R 4 R 1 C 2 + R 3 L 4 R 1 C 2
+ R3 R 4 L 1 C 2 + L 1 R 3 C 1 R4 + L 4 R 4 R1 C 2 + R 1 L 3 C 1 R4
+ R1 L4 C1 R4 + L3 L4 + R1 R3 C1 L4 ) + S(R3 R4 R2 C2 + R4 R1 R2 C2
+ L 4 R 4 + L 3 R 4 + R3 L 4 + L 1 R4 + R 3 R 4 R1 C 2 + R 1 R3 C 1 R 4
+ L4 R1 ) + R1 R4 + R3 R4 (7)
212
Mathematical Derivations and MATLAB Scripts
L2 C2
E i2 L1 i1 C1 i3 L4
Figure A4. The equivalent circuit for ground fault transient model 4
The equations for voltages around the loops of the circuit are in Laplace-
domain as follows
1 1
i1 (s)(SL1 + ) − i2 (s)SL1 − i3 (s) =0 (1)
SC1 SC1
1 E
−i1 (s)SL1 + i2 (s)(SL1 + SL2 + )= (2)
SC2 S
1 1
−i1 (s) + i3 (s)(SL4 + )=0 (3)
SC1 SC1
From which the current i2 (s) can be solved in the following steps, First
Solving i3 (s) from equation (3) and substituting it into equation (1). Next,
solving i1 (s) from equation (1) and substituting it into equation (2). After
the substitutions and rearranging the terms, equation (2) becomes:
where,
and
Q(s) = S 4 L2 C2 L1 L4 C1 + S 2 (L1 C2 L4 + L2 C2 L4 + L2 C2 L1 + L1 L4 C1 )
+L1 + L4 (6)
213
Mathematical Derivations and MATLAB Scripts
214
B. ATPDraw OH Line Configurations
and Load Variation
1.3m 1.3m
1.1m 1.1m
1.1m
215
ATPDraw OH Line Configurations and Load Variation
216
C. Tables of Simulation Results
Table C1. Comparison of simulation test result to the theory of earth fault transient.
fdc , fnc and ζc are calculated using the GM1 equations. fdm , fnm and ζm
are estimated transient parameters from the simulated earth fault transient
waveforms.
l fd fn
ζ
(km) (Hz) (Hz)
Notes: l = Actual fault distance (km), fd = Damped transient frequency (Hz), fdm = Estimated
damped transient frequency (Hz), fdc = Calculated damped transient frequency (Hz), ∆ fd = Error
of damped transient frequency (Hz), ζ = Damping factor, ζm = Estimated damping factor, ζc = Cal-
culated damping factor, ∆ζ = Error of damping factor, fn = Undamped transient frequency (Hz),
fnm = Estimated undamped transient frequency (Hz), fnc = Calculated undamped transient fre-
quency (Hz), ∆ fn = Error of undamped transient frequency (Hz), MAE = Mean average error.
217
Tables of Simulation Results
Notes: l = Actual fault distance (km), GM1f dm = Estimated fault distance using
general model 1 and damped transient frequency (km), ∆GM1f dm = Error of esti-
mated fault distance using general model 1 and damped transient frequency (km),
GM1f nm = Estimated fault distance using General model 1 and undamped tran-
sient frequency (km), ∆GM1f nm = Error of estimated fault distance using general
model 1 and undamped transient frequency (km), GM2f nm = Estimated fault dis-
tance using General model 2 and undamped transient frequency (km), ∆GM2 = Er-
ror of estimated fault distance using general model 2 and undamped transient fre-
quency (km), MAE = Mean average error.
218
Tables of Simulation Results
Table C3. Comparison of simulation test result to the theory of earth fault transient with the
effect of fault resistance. fdc , fnc and ζc are calculated using GM1 equations. fdm
, fnm and ζm are estimated transient parameters from the simulated earth fault
transient waveforms.
l Rf fd fn
ζ
(km) (Ω) (Hz) (Hz)
Notes: l = Actual fault distance (km), fd = Damped transient frequency (Hz), fdm = Estimated damped
transient frequency (Hz), fdc = Calculated damped transient frequency (Hz), ∆ fd = Error of damped tran-
sient frequency (Hz), ζ = Damping factor, ζm = Estimated damping factor, ζc = Calculated damping factor,
∆ζ = Error of damping factor, fn = Undamped transient frequency (Hz), fnm = Estimated undamped tran-
sient frequency (Hz), fnc = Calculated undamped transient frequency (Hz), ∆ fn = Error of undamped
transient frequency (Hz), Rf = Fault resistance (Ω), MAE = Mean average error.
219
Tables of Simulation Results
Table C4. Comparison of the distance estimates obtained of the GM algorithm meth-
ods with effect of fault resistance.
Notes: l = Actual fault distance (km), GM1f dm = Estimated fault distance using general
model 1 and damped transient frequency (km), ∆GM1f dm = Error of estimated fault distance
using general model 1 and damped transient frequency (km), GM1f nm = Estimated fault dis-
tance using General model 1 and undamped transient frequency (km), ∆GM1f nm = Error of
estimated fault distance using general model 1 and undamped transient frequency (km), GM2
= Estimated fault distance using General model 2 and undamped transient frequency (km),
∆GM2 = Error of estimated fault distance using general model 2 and undamped transient fre-
quency (km), Rf = Fault resistance (Ω), MAE = Mean average error.
220
Tables of Simulation Results
Table C5. Comparison of simulation test result to the theory of ground fault transient
with effect of fault inception angle. fdc , fnc and ζc are calculated using the
GM1 equations. fdm , fnm and ζm are estimated transient parameters from
the simulated earth fault transient waveform.
l Rf fd fn
ζ
(km) (Ω) (Hz) (Hz)
Notes: l = Actual fault distance (km),fd = Damped transient frequency (Hz), fdm = Estimated
damped transient frequency (Hz), fdc = Calculated damped transient frequency (Hz), ∆ fd = Er-
ror of damped transient frequency (Hz), ζ = Damping factor, ζm = Estimated damping factor, ζc
= Calculated damping factor, ∆ζ = Error of damping factor, fn = Undamped transient frequency
(Hz), fnm = Estimated undamped transient frequency (Hz), fnc = Calculated undamped tran-
sient frequency (Hz), ∆ fn = Error of undamped transient frequency (Hz), Rf = Fault resistance
(Ω), MAE = Mean average error.
221
Tables of Simulation Results
Table C7. Comparison of simulation test result to the theory of ground fault transient with the effect
of composite loads. fdc , fnc and ζc are calculated using the general model 3 equations. fdm
, fnm and ζm are estimated transient parameters from the simulated earth fault transient
waveforms.
l Load Rf fd fn
ζ
(km) (MVA) (Ω) (Hz) (Hz)
Notes: l = Actual fault distance (km), fd = Damped transient frequency (Hz), fdm = Estimated damped transient
frequency (Hz), fdc = Calculated damped transient frequency (Hz), ∆ fd = Error of damped transient frequency (Hz),
ζ = Damping factor, ζm = Estimated damping factor, ζc = Calculated damping factor, ∆ζ = Error of damping factor,
fn = Undamped transient frequency (Hz), fnm = Estimated undamped transient frequency (Hz), fnc = Calculated
undamped transient frequency (Hz), ∆ fn = Error of undamped transient frequency (Hz), Rf = Fault resistance (Ω),
MAE = Mean average error.
222
Tables of Simulation Results
Table C8. Comparison of the distance estimates obtained of the GM algorithm methods with
effect of composite loads.
l (km) MVA Rf (Ω) GM1 ∆GM1 GM2 ∆GM2 GM3 ∆GM3 GM4 ∆GM4
Notes: l = Actual fault distance (km), GM1 = Estimated fault distance using General model 1 and
undamped transient frequency (km), ∆GM1 = Error of estimated fault distance using general model 1
and undamped transient frequency (km), GM2 = Estimated fault distance using General model 2 and un-
damped transient frequency (km), ∆GM2 = Error of estimated fault distance using general model 2 and
undamped transient frequency, GM3 = Estimated fault distance using General model 3 and undamped
transient frequency (km), ∆GM3 = Error of estimated fault distance using general model 3 and undamped
transient frequency (km), GM4 = Estimated fault distance using General model 4 and undamped tran-
sient frequency (km), ∆GM4 = Error of estimated fault distance using general model 4 and undamped
transient frequency (km), MAE = Mean average error.
223
Tables of Simulation Results
Notes:
l = Actual fault distance (km), fdm = Estimated damped transient fre-
quency (Hz), ζ = Damping factor, fnm = Estimated undamped transient
frequency (Hz), EM1 = Estimated fault distance using exact model 1 (km)
, ∆EM1 = Error of estimated fault distance using exact model 1 (km), GM2
= Estimated fault distance using General model 2, ∆GM2 = Error of es-
timated fault distance using general model 2 (km), MAE = Mean average
error.
224
Tables of Simulation Results
Table C11. Comparison of the distance estimates obtained between EM1 and
GM2 algorithm with the effect of fault inception angle. fdm , fnm
and ζ are estimated transient parameters from the simulated earth
fault transient waveform.
Notes: l = Actual fault distance (km), fdm = Estimated damped transient frequency (Hz),
ζ = Damping factor, fnm = Estimated undamped transient frequency (Hz), EM1 = Estimated
fault distance using exact model 1 (km) , ∆EM1 = Error of estimated fault distance using
exact model 1 (km), GM2 = Estimated fault distance using General model 2, ∆GM2 = Error
of estimated fault distance using general model 2 (km), MAE = Mean average error.
225
Tables of Simulation Results
Table C12. Comparison of the distance estimates obtained between EM1 and GM2 algorithm with
the effect of resistive loads. The loads are located at the end of each feeders. fdm ,
fnm and ζ are estimated transient parameters from the simulated earth fault transient
waveform.
Notes: l = Actual fault distance (km), fdm = Estimated damped transient frequency (Hz), ζ = Damping factor,
fnm = Estimated undamped transient frequency (Hz), EM1 = Estimated fault distance using exact model 1 (km)
, ∆EM1 = Error of estimated fault distance using exact model 1 (km), GM2 = Estimated fault distance using Gen-
eral model 2, ∆GM2 = Error of estimated fault distance using general model 2 (km), MAE = Mean average error.
226
Tables of Simulation Results
Table C13. Comparison of the distance estimates obtained between EM1 and GM2 algorithm
with the effect of inductive loads. The loads are located at the end of each feeders.
fdm , fnm and ζ are estimated transient parameters from the simulated earth fault
transient waveform.
Notes: l = Actual fault distance (km), fdm = Estimated damped transient frequency (Hz), ζ = Damping
factor, fnm = Estimated undamped transient frequency (Hz), EM1 = Estimated fault distance using exact
model 1 (km) , ∆EM1 = Error of estimated fault distance using exact model 1 (km), GM2 = Estimated fault
distance using General model 2, ∆GM2 = Error of estimated fault distance using general model 2 (km), MAE
= Mean average error.
227
Tables of Simulation Results
Table C14. Comparison of the distance estimates obtained between EM1 and GM2 algo-
rithm with the effect of loads of 0.9 power factor. The loads are located at the
end of each feeders. fdm , fnm and ζ are estimated transient parameters from
the simulated earth fault transient waveform.
Table C15. Comparison of the distance estimates obtained between EM1 and GM2 algorithm with the
effect of loads of 0.5 power factor. The loads are located at the end of each feeders. fdm
, fnm and ζ are estimated transient parameters from the simulated earth fault transient
waveform
Notes: l = Actual fault distance (km), fdm = Estimated damped transient frequency (Hz), ζ = Damping factor,
fnm = Estimated undamped transient frequency (Hz), EM1 = Estimated fault distance using exact model 1 (km) ,
∆EM1 = Error of estimated fault distance using exact model 1 (km), GM2 = Estimated fault distance using General
model 2, ∆GM2 = Error of estimated fault distance using general model 2 (km), MAE = Mean average error.
228
Tables of Simulation Results
229
Tables of Simulation Results
230
Tables of Simulation Results
Table C18. Comparison of the distance estimates obtained from general model
(GM1,GM2,GM4) algorithms with the effect of fault resistance
(Rf ). Signal measured from MV/LV side and fault inception an-
gle is 90◦ .
Notes: l = Actual fault distance (km), GM1 = Estimated fault distance using general
model 1 (km), ∆GM1 = Error of estimated fault distance using general model 1 (km),
GM2 = Estimated fault distance using general model 2 (km), ∆GM2 = Error of estimated
fault distance using general model 2 (km), GM4 = Estimated fault distance using general
model 4 (km), ∆GM4 = Error of estimated fault distance using general model 4 (km),
MAE = Mean average error (km).
231
Tables of Simulation Results
Table C19. Comparison of the distance estimates obtained from general model
(GM1, GM2,GM4) algorithms with the effect of fault resistance (Rf ).
Signal measured from MV/LV side and fault inception angle is 45◦ .
Notes: l = Actual fault distance (km), GM1 = Estimated fault distance using general model
1 (km), ∆GM1 = Error of estimated fault distance using general model 1 (km), GM2 = Esti-
mated fault distance using general model 2 (km), ∆GM2 = Error of estimated fault distance
using general model 2 (km), GM4 = Estimated fault distance using general model 4 (km),
∆GM4 = Error of estimated fault distance using general model 4 (km), MAE = Mean average
error (km).
232
Tables of Simulation Results
Table C20. Comparison of the distance estimates obtained from exact model
(EM1-GM3) algorithms with the effect of fault resistance (Rf ). Sig-
nal measured from MV/LV side and fault inception angle is 90◦ .
Notes: l = Actual fault distance (km), EM1 = Estimated fault distance using exact model 1
(km), ∆EM1 = Error of estimated fault distance using exact model 1 (km), EM2 = Estimated
fault distance using exact model 2 (km), ∆EM2 = Error of estimated fault distance using Ex-
act model 2 (km), EM3 = Estimated fault distance using exact model 3 (km), ∆GM3 = Error
of estimated fault distance using exact model 3 (km), MAE = Mean average error (km).
233
Tables of Simulation Results
Table C21. Comparison of the distance estimates obtained from exact model
(EM1-GM3) algorithms with the effect of fault resistance (Rf ). Sig-
nal measured from MV/LV side and fault inception angle is 45◦ .
Notes: l = Actual fault distance (km), EM1 = Estimated fault distance using exact model 1
(km), ∆EM1 = Error of estimated fault distance using exact model 1 (km), EM2 = Estimated
fault distance using exact model 2 (km), ∆EM2 = Error of estimated fault distance using Ex-
act model 2 (km), EM3 = Estimated fault distance using exact model 3 (km), ∆GM3 = Error
of estimated fault distance using exact model 3 (km), MAE = Mean average error (km).
234
De
par
tme
nto
fEl
ect
ric
alE
ngi
nee
ringa
ndAut
oma
ti
on
Mo
al
t
eart hfaul ts in me dium vo ltage( MV)
h
-
D
d Rafi Adz
Computation Me tho ds
D1
dist ribut io n fe e de rs.Singl eph aset oe art h
89
faul tis th emo stc o mmo n faul tt ypein po w er
/
20
dist ribut io n ne tw o rks.I nF inland and in
1
4
Basedo
nT ransients in
mo sto fc o unt rie s in Euro pet heMV
man
ne tw o rks areo pe rat edw ithne ut raliso l
ated
or re so nante art h e d syst e ms.T h ew ay the
PowerDistribution Syste
ms
Eart
ne ut ralis c o nne c te d is impo rtantsinc eit
w illaffe c tth et ypeo f pro t e ction syst em
hFaul
ne ede d and po w e r syst e mc o mpo ne nt
tD
se l
e c t
io n.I n iso lat e d ne ut ralne tw orks e arth
ist
faul tc urre nt s arere lative l y smal lc o mpare d
anc
w ithe arth e d syst e ms.T h einit ialt ransie nts
eCo
ofe art hfaul t s areimpo rt antfo r faul t Mo
hdRa
fiAd
zma
n
mput
locat io ne spe c ially in une art h e d and
co mpe nsat ed ne ut ralne tw o rk.T h e
atio
elec tro magne tict ransie ntpro videt h e
n Me
avail abl einfo rmat io n fo re st imat ing faul t
t
locat io n.T his disse rtat io n pro po se s fault
ho
ds Base
locat io n algo rit h ms base do ne art hfaul t
transie ntsignal s me asure d atse co ndary
sideo f MV and L V dist ribut io nt ransfo rme rs
donT
in orde rt ol oc at ea singl el inet oe art hfaul t
ransie
in une arth e d and c o mpe nsat ed ne ut ralMV
dist ribut io n ne tw o rks.
nts in P
owe
rDist
ribut
ion Syst
ems
I
S BN978-
95 2
-60-
5973-0(pri
nted) BUSI
N E
S S+
9HSTFMG*afjhda+
I
S BN978-
95 2
-60-
5974-7(pdf) ECONOM Y
I
S SN-
L1799-4934
I
S SN1799-
4934(pri
nted) ART+
I
S SN1799-
4942(pdf) DESI
G N+
ARCH I
TECT
U RE
A
Aa
lt
oU ni
versi
ty
a
Sc
hoolofElect
ric
alE ngi
ne e
ring SCI
ENCE+
lt
oU
Depa
rtmentofElect
rica
lE ngi
neeri
nga
ndAut
oma
ti
on T
EC HNOL
O GY
www.
aalt
o.f
i
ni
CROSSOV
ER
ve
DOCTORAL
r
DOCTORAL DI
S SE
R T
ATI
O NS
sit
DI
S SE
R T
ATI
O NS
y