Electric Power Systems Research 80 (2010) 176183

Contents lists available at ScienceDirect

Electric Power Systems Research

journal homepage: www.elsevier.com/locate/epsr

Assessment of surge arrester failure rate and application studies in Hellenic high
voltage transmission lines
C.A. Christodoulou a , L. Ekonomou b, , G.P. Fotis a , I.F. Gonos a , I.A. Stathopulos a
a
b

National Technical University of Athens, School of Electrical and Computer Engineering, High Voltage Laboratory, 9 Iroon Politechniou St., Zografou Campus, 157 80 Athens, Greece
A.S.PE.T.E. - School of Pedagogical and Technological Education, Department of Electrical Engineering Educators, N. Heraklion, 141 21 Athens, Greece

a r t i c l e

i n f o

Article history:
Received 18 December 2007
Received in revised form 15 July 2009
Accepted 27 August 2009
Available online 25 September 2009
Keywords:
Surge arresters
Backashover
Ground ash density
Lightning performance
Lightning protection
Optical transient density
Overhead transmission lines
Shielding failures

a b s t r a c t
The use of transmission line surge arresters to improve the lightning performance of transmission lines
is becoming more common. Especially in areas with high soil resistivity and ground ash density, surge
arresters constitute the most effective protection mean. In this paper a methodology for assessing the
surge arrester failure rate based on the electrogeometrical model is presented. Critical currents that
exceed arresters rated energy stress were estimated by the use of a simulation tool. The methodology is
applied on operating Hellenic transmission lines of 150 kV. Several case studies are analyzed by installing
surge arresters on different intervals, in relation to the regions tower footing resistance and the ground
ash density. The obtained results are compared with real records of outage rate showing the effectiveness
of the surge arresters in the reduction of the recorded failure rate. The presented methodology can be
proved valuable to the studies of electric power systems designers intending in a more effective lightning
protection, reducing the operational costs and providing continuity of service.
2009 Elsevier B.V. All rights reserved.

1. Introduction
Lightning strikes are the main reason for outages in overhead
transmission lines. In an effort to maintain high power quality and
to avoid damages and disturbances, overhead ground wires and
surge arresters are used for the transmission lines protection.
Numerous different methodologies have been presented in the
technical literature in an effort to assess the lightning performance
and improve the lightning protection of transmission lines. These
methodologies are extended from the use of analogue computer
methods [1], Monte-Carlo simulation techniques [2] and travelling wave methods [3], to the use of electrogeometrical models
[4,5], simulation software [6,7] and articial neural networks [8,9].
Although there are many papers dealing with this subject, only
a few consider the presence of protection devices such as surge
arresters [1015]. The reason is obvious since the use of surge
arresters is not a common practice in all electric utilities. Even in
countries such as the USA and Japan where surge arresters have
installed since the 1980s [16,17], these are not utilized broadly by
all USA or Japanese electric utilities.
In the current work the contribution of surge arresters in the
lightning protection of transmission lines is studied. A methodology

Corresponding author. Tel.: +30 697 2702218; fax: +30 210 7723504.
E-mail address: [email protected] (L. Ekonomou).
0378-7796/$ see front matter 2009 Elsevier B.V. All rights reserved.
doi:10.1016/j.epsr.2009.08.018

is presented for the computation of the transmission lines lightning failures including the failure probability and the failure rate
of the surge arresters. Considering that a surge arrester is the last
protection mean in a transmission line, an arrester failure as well
as backashover and shielding failures were considered as transmission line faults. The proposed methodology, which is based on
the electrogeometrical model, uses also a simulation tool in order
to estimate the critical currents that exceed arresters rated energy
stress. The methodology is applied on six operating Hellenic transmission lines of 150 kV, of known outage rate, considering surge
arresters for every 1, 2 or 3 towers. It must be mentioned that none
of the Hellenic transmission lines, either the 150 kV or 400 kV, is
equipped with surge arresters. The obtained results show clearly
that the use of surge arresters reduces signicantly the lightning
faults and their use can improve the lightning performance of them.
2. Lightning parameters
The lightning parameters, which have been used in this study,
are based on the measurements performed by Berger in Monte
San Salvatore [18]. In order to simulate lightning strokes, lightning parameters such as the peak value of lightning current and
the steepness, are randomly selected from the statistical distributions, using the well known Monte-Carlo statistical technique [2]
and the methodology presented in detail in [19]. According to the
Lightning and Insulator Subcommittee of the T&D Committee [20],

C.A. Christodoulou et al. / Electric Power Systems Research 80 (2010) 176183

177

Fig. 1. Electrogeometrical model: representation of ground wires and phase conductors.

90% of the lightning strokes are considered negative while 10% of

the lightning strokes are considered positive, something which has
also adopted for the calculations of this study.
3. Electrogeometrical model
The termination point of a lightning stroke to a transmission
line can be either a ground wire, phase conductor, tower or even
ground. The electrogeometrical model using the concept of striking
distance has got the ability to determine the termination point. In
general, the striking distance r in m is given by the equation:
r = AI b

(1)

where A and b are constants dependent on the termination point

and I is the prospective stroke current in kA.
Although there are several versions of electrogeometrical
model, where each one uses different values for the constants A
and b, all of them consider the following three concepts: (a) strokes
arrive vertically, (b) the lightning leader develops unaffected by
the existence of grounded objects until it arrives within striking
distance from the grounded object and (c) the striking distance is
related to the current of the return stroke. In this study the electrogeometrical model proposed by IEEE Working Group [5] was used
in the calculations.
The fractions of lightning strikes hA (IP ) and hB (IP ) that will terminate, respectively, on a phase conductor or on an overhead ground
wire, can be estimated for each value of the peak current IP using
Eqs. (2) and (3):
hA (Ip ) =

DC
D

(2)

hB (Ip ) =

D DC
D

(3)

where DC is the shielding failure exposure distance shown in Fig. 1

and D is the vertical distance between the intersection point C and
the intersection point A (Fig. 1).
4. Surge arresters
Surge arresters are designed to be insulators for nominal operating voltage, conducting at most a few milliamperes of current and
good conductors when the voltage of the line exceeds design specications to pass the energy of the lightning strike to the ground.

Several different types of arresters are available (e.g. gapped

silicon carbide, gapped or non-gapped metal-oxide) and all perform in a similar manner: they function as high impedances at
normal operating voltages and become low impedances during
surge conditions. Even though a great number of arresters which
are gapped arresters with resistors made of silicon carbide (SiC) are
still in use, the arresters installed today are almost all metal-oxide
(MO) arresters without gaps, something which means arresters
with resistors made of metal-oxide [21]. The distinctive feature of a
metal-oxide arrester is its extremely non-linear VI characteristic,
rendering unnecessary the disconnection of the resistors from the
line through serial spark gaps, as it is found in the arresters with
SiC resistors.
The most signicant technical characteristics of surge arresters
according to the IEC 60099-4 are [21,22]:
(a) Continuous operating voltage (Uc ). Designated rms value of
power frequency voltage that may be applied continuously
between the terminals of the arrester. MCOV of the arrester
must be higher than the maximum continuous operating voltage of the system.
(b) Rated voltage. Maximum permissible rms value of power frequency voltage between arrester terminals at which is designed
to operate correctly under temporary overvoltages.
(c) Discharge current. Impulse current which ows through the
arrester.
(d) Residual voltage (Ures ). Peak value of the voltage that appears
between arrester terminals when a discharge current is
injected.
(e) Rated discharge current. Peak value of lightning current
impulse, which is used to classify an arrester.
(f) Lightning impulse protective level. Voltage that drops across
the arrester when the rated discharge current ows through
the arrester.
(g) Energy absorption capability. Maximum level of energy injected
into the arrester at which it can still cool back down to its normal
operating temperature. Standards do not dene energy capability of an arrester. In IEC exists the term line discharge class,
but since this in not enough information, various manufacturers present thermal energy absorption capability in kJ/kV (Uc ),
dened as the maximum permissible energy that an arrester
may be subjected to two impulses according to IEC clause 8.5.5
[22], without damage and without loss of thermal stability
[23,21].

178

C.A. Christodoulou et al. / Electric Power Systems Research 80 (2010) 176183

5. Failure probability and failure rate of a surge arrester

The lightning energy E (in Joules) absorbed by a surge arrester
is computed by the relation:

E=

u(t) i(t)dt

(4)

to

where u(t) is the residual voltage of the arrester in kV and, i(t) is

the value of the discharge current through the arrester in kA.
When the arrester energy absorption exceeds its withstand
capability, the arrester is damaged (failure). Considering the energy
required to cause damage to an arrester (4), the failure probability
of an arrester is calculated by Eqs. (5)(7) [13,2426]:



PA =

f (Ip ) hA (IP )dIp

Tr



IA (Tt )

f (Ip ) hB (IP )dIp

Tr

(5)

g(Tt )dTt

(6)

PB =

g(Tt )dTt

of the ashover of the insulator strings, the vt curve is used recommended by IEEE [31]: VD = (400 + (710/t0.75 )) W, where VD is the
ashover voltage, t is the time to ashover and W is the insulators
string length. The tower footing resistance is modeled as a lumped
resistance. As far concerning the arresters it was used the element
provided by Simulink, due to the difculties and the limitations in
order to use one of the existing models [30,31].
As far concerning the estimation of the ground ash density, a
new method was used based on the optical transient density [27].
The ground ash density is estimated simple and with accuracy
using NASA satellites observations of average optical ash density for the geographical region of Hellas [33,34] and the nominal
relation of 4:1 between optical transient density and ground ash
density as it is proposed in [32]. It must be noticed the great advantage of this method, since NASAs databases contain not only annual
but also seasonal and monthly observation data.
6. Shielding and backashover failure rate

IB (Tt )

PT = PA + PB

(7)

AT = NL l PT

(8)

where PA is the probability that an arrester fails due to lightning

stroke on a phase conductor, PB is the probability that an arrester
fails due to lightning stroke on the overhead ground wire, IA (Tt ) is
the minimum stroke peak current in kA required to damage the
arrester, when lightning hits on a phase conductor, depending on
each time-to-half value, IB (Tt ) is the minimum stroke peak current
in kA required to damage the arrester, when lightning hits on the
overhead ground wire, depending on each time-to-half value, f(IP ) is
the probability density function of the lightning current peak value,
g(Tt ) is the probability density function of the time-to-half value of
the lightning current, Tr is the rise time of the incident waveform, PT
is the total failure probability of an arrester, AT is the arrester total
failure rate in failures per year per line, NL is the number of lightning
ashes to a line per 100 km per year, equal to: (Ng /10)(28 h0.6
+ g),
t
ht is the tower height in m, g is the horizontal spacing in m, between
the ground wires, Ng is the ground ash density in ashes per km2
per year and l is the line length in km.
A lightning strike to a phase conductor, to a tower or to an overhead ground wire can cause arrester failure. For direct strikes to
phase conductors, low tower footing resistance results in higher
energy through the arrester; for strikes on towers or on ground
wires, high resistances increase the arresters failure probability.
The probability of failure due to strike on phases is small, so
backashover is the dominant case for the determination of the
arresters failure rate. Thus, the arrester discharge current i(t) and
the energy E depend on the tower footing resistance value, with the
lightning peak current (that allows the arrester to reach its energy
absorption capability) to decrease as the tower footing resistance
increases.
The stroke peak current for each time-to-half value Tt (depending on the exceeded required energy for the arrester failure), can
be calculated analytically or using an appropriate simulation tool
[1113,2426]. In this work MATLAB & Simulink were used [27],
where the transmission lines were represented by a distributed
parameter line model, which is based on the Bergerons travelling
wave method [28] used by the EMTP [29] and represents wave
propagation phenomena and line end reections efciently. The
surge impedance values are 400  for phases and 700  for ground
wires. The towers are represented as distributed parameters lines
equal to 200 . The dielectric strength of the insulator was represented as a voltage-controlled switch. When the voltage U exceeds
the critical ashover voltage the switch closes. For the evaluation

Shielding failure rate NSF is associated to a required minimum

current Imin to cause a line insulation ashover [35,36]. NSF in failures per year per line is dened as follows:
NSF =

2Ng l
10

Imax

DC f (I)dI

(9)

Imin

where f(I) is the probability density function for the current, DC is

the shielding failure exposure distance shown in Fig. 1, Imax is the
maximum lightning current in kA that the ground wire will allow
to strike the phase conductor due to the placement of the ground
wire, Imin is the minimum current equal to 2Ua /Zsurge [35], Ua is
the insulation level of the transmission line in kV and Zsurge is the
conductor line surge impedance.
Backashover failure rate is estimated for transmission lines,
with or without ground wires, according to the method developed
by the members of High Voltage Laboratory of National Technical
University of Athens, and is presented in detail in [19,36,37,38,39].
NBF in failures per year per line is dened as follows:

(Ipeak )

NBF = NL l

max

(di/dt)max

P()dIpeak d(di/dt)
(Ipeak )

(10)

(di/dt)min

min

where P() is the probability distribution function of the random

variable , which is a function of the two random variables Ipeak
and di/dt as shown in the following relation:

Ipeak ,

di
dt

=R

Ipeak
2

0.85Ua + L

di
dt

(11)

with greater than zero when there is a backashover, R is the

tower footing resistance in , L is the total equivalent inductance
of the system (tower and grounding systems inductance) in H,
calculated according to the simplied method presented in [5], di/dt
is a random variable denoting the lightning steepness in kA/s and
Ipeak is a random variable denoting the peak lightning current in kA.
7. Total lightning failure of transmission lines
The total lightning failure of a transmission line in failures per
year per line is the summation of any arrester failure, shielding
failure and failure due to backashover. Thus:
NTOTAL = AT + NSF + NBF

(12)

8. Application of the proposed methodology

The methodology presented in this paper has been applied and
tested on six 150 kV operating transmission lines of the Hellenic

C.A. Christodoulou et al. / Electric Power Systems Research 80 (2010) 176183

179

Table 1
Line characteristics of the analyzed transmission lines.
No.

Line

Voltage (kV)

Length (km)

No. of towers

Insulation level (kV)

Conductor dimensions (ACSR MCM)

No. of circuits

1
2
3
4
5
6

Ioannina-Kalpaki
Igoumenitsa-Sagiada
Kilkis-Serres
Arachthos-Igoumenitsa
Megalopoli-Sparti
Aktio-Argostoli

150
150
150
150
150
150

28
18
58
76
64
81

86
44
162
239
173
224

750
750
750
750
750
750

336.4
336.4
336.4
336.4
336.4
336.4

1
1
1
1
1
1

Table 2
Line design parameters of the analyzed transmission lines.
Line

Region

Towers

R (tower footing resistance in )

Ioannina-Kalpaki

I
II

140
4186

4.2
25.8

3.824
3.140

Igoumenitsa-Sagiada

I
II

129
3044

57.5
14.9

4.353
2.599

Kilkis-Serres

I
II
III

146
47106
107162

2.0
4.4
1.8

3.284
2.114
2.646

Arachthos-Igoumenitsa

I
II
III

180
81163
164239

5.2
13.0
45.4

3.540
4.188
2.950

Megalopoli-Sparti

I
II
III

145
4675
76173

5.1
39.7
11.2

3.700
3.430
3.381

Aktio-Argostoli

I
II
III

155
56137
138224

4.8
64.9
126.3

3.540
3.248
2.773

Table 3
Technical characteristics of the used surge arresters.
Continuous operating voltage (Uc )
Rated voltage
Rated discharge current

108 < Uc < 115 kV (rms)

144 kV (rms)
10 kA

Residual voltage

<330 kV (max) for 5 kA

<350 kV (max) for 10 kA
<390 kV (max) for 20 kA

Discharge energy class

Energy capability

3
8 kJ/kV

interconnected system. These lines, were carefully selected among

others, due to: (a) their high failure rates during lightning thunderstorms [39], (b) their consistent construction for at least 90% of
their length, (c) their sufcient length and their sufcient time in
service in order to present a reasonable exposure to lightning and
(d) the signicant different characteristics, such as the ground ash
density and the tower footing resistance, which exist through their
length, since they run at the same time through a plain region, a
coastline and/or a mountainous region.
Table 1 presents the Hellenic high voltage transmission lines
and their characteristics, which have been used in this study [39].
Table 2 presents the same lines divided into regions, showing
clearly the signicant different values of ground ash density [36]
and tower footing resistance that exist in each one of them [39].
It must be mentioned that none of the operating Hellenic transmission lines, either the 150 kV or 400 kV, is equipped with surge
arresters. Therefore the application of the proposed methodology
in these particular lines was a very challenging work.
The technical characteristics of the surge arresters used in the
analysis of the transmission lines are shown in Table 3. Using MATLAB & Simulink three different case studies were analyzed. In the
rst case, surge arresters were installed on all the towers of each
one region of the transmission lines and for each one of the three
phases. In the second and third case studies the surge arresters

Ng (average ground ash density 19952000)

interval for each one of the analyzed lines and for each one region
was set to every second and every third tower, respectively, for each
one of the three phases. In the boundary of two regions of a transmission line, where the tower footing resistance varies signicantly
(e.g. regions I and II of the transmission line Igoumenitsa-Sagiada)
or the ground ash density presents signicant differences (e.g.
regions II and III of the transmission line Arachthos-Igoumenitsa)
the surge arresters were installed on both towers, which were near
to the boundary, independently from the surge arrester interval.
It must be mentioned that in this study lightning ashes were
considered to strike only the towers in which surge arresters were
installed. Although this is not what really happens in the transmission lines (ashover can also occur on a tower without arrester
either due to on them or due to lightning strike on neighbouring
towers equipped with arresters), this is the worst case for the surge
arresters, which is the main issue of this study.
9. Results and discussion
In Table 4 are presented the eld observation data, i.e., actual
recorded transmission lines lightning failures, the estimated using
the proposed methodology transmission lines failure rate without
the use of arresters and the obtained failure rate of the examined
transmission lines for each one of the three analyzed cases, i.e.,
surge arresters are installed on every single, second or third tower.
In general, the use of surge arresters on transmission lines
certainly reduces the backashovers and the shielding failures.
However, in the total transmission line faults must also be included
the arrester damages, since arresters consist equipment of the line
and they need repair or replacement after their failure. Additionally,
a damaged arrester is not effective any more and can create new
problems if it continues to remain on the line. In the current work,
considering the arresters failures as line faults, the failure rates after
the arresters installation are compared with real failure records
and the estimated failures. The implementation of surge arresters

180

C.A. Christodoulou et al. / Electric Power Systems Research 80 (2010) 176183

Table 4
Field observation data versus obtained results for years 19952000.
Line

Region

Average eld observation

failure ratea

Failure rate without

arresters

Failure rate for

Case 1b

Failure rate for

Case 2c

Failure rate for

Case 3d

Ioannina-Kalpaki

I
II

0.5
1.0

0.462
0.974

0.102
0.449

0.255
0.776

0.374
0.790

Total

1.5

1.436

0.551

1.031

1.164

I
II

2.2
0.6

2.305
0.648

0.786
0.121

1.668
0.304

1.870
0.410

Total

2.8

2.953

0.907

1.972

2.280

I
II
III

0.9
1.9
0.7

0.955
1.884
0.724

0.063
0.147
0.049

0.218
0.513
0.255

0.356
0.705
0.436

Total

3.5

3.563

0.259

0.986

1.497

I
II
III

0.6
1.5
3.1

0.613
1.540
3.207

0.425
0.769
1.005

0.761
1.288
1.890

0.813
1.740
2.487

Total

5.2

5.360

2.199

3.939

5.040

I
II
III

0.5
0.9
2.3

0.461
0.836
2.199

0.277
0.439
0.753

0.424
0.723
1.332

0.400
0.898
1.815

Total

3.7

3.496

1.469

2.479

3.113

I
II
III

0.0
2.1
3.7

0.008
2.053
3.672

0.256
1.690
2.562

0.329
1.971
2.977

0.389
2.093
3.366

Total

5.8

5.733

4.508

5.277

5.848

Igoumenitsa-Sagiada

Kilkis-Serres

Arachthos-Igoumenitsa

Megalopoli-Sparti

Aktio-Argostoli

a
b
c
d

Average eld observation failure rate means recorded lightning failures (actual data).
Case 1: The surge arresters are installed on every single tower.
Case 2: The surge arresters are installed on every second tower.
Case 3: The surge arresters are installed on every third tower.

improves the lightning performance of the line, especially when

they are installed on every single or second tower. For arresters
installed every third tower, the reduction of the line failure rate is
not great enough and is possible for regions with high tower footing
resistance the total failure rate to increase (e.g. Aktio-Argostoli).
Fig. 2 shows the surge arresters failure probability in relation
to the tower footing resistances, for all the lines, which have the
same structural characteristics, for each one of the three analyzed
case studies, since the arresters failure probability depends only on
tower footing resistance for a given interval. As it was expected,
low resistance, corresponds to lower arresters failure probability,
with the failure rate to be depended on the ground ash density
and the line length. In any case the tower footing resistance must
be low enough, in order to minimize the surge arresters failure
probabilities and the transmission line lightning failures.

Fig. 3 shows the relationship between the installation interval

and arresters failure probabilities for the line Igoumenitsa-Sagiada.
It is observed that smaller interval decreases the arresters failure
rate. This is also obvious in Fig. 1, where for a given tower footing
resistance the surge arresters failure probability is higher as the
interval increases. Thus, for regions with high tower footing resistance (e.g. Aktio-Argostoli) arresters installation on every tower is
suggested.
Fig. 4 shows the percentage failure reduction for all the analyzed lines in each one of the three examined case studies, using
the equation:

F.R.R. (%) =

F.R.o F.R.i
100%
F.R.0

Fig. 2. The variation of surge arresters failure probability with tower footing resistance for each one of the three analyzed case studies.

(13)

C.A. Christodoulou et al. / Electric Power Systems Research 80 (2010) 176183

181

Fig. 3. The variation of surge arresters failure probability with the arresters interval for the transmission line Igoumenitsa-Sagiada.

Fig. 4. The obtained failure rates reduction for all the analyzed lines for each one of the three analyzed case studies.

where F.R.R. is the failure rate reduction, F.R.o is the failure rate
without surge arresters, i is equal to 1, 2 or 3 denoting the three
examined case studies and F.R.i is the failure rate with installed
surge arresters for the examined case study.
In lines with low tower footing resistance (e.g. Kilkis-Serres)
the obtained failure rate after surge arresters installation on every
tower is approximately ten times lower than the failure rate
without arresters. In lines with high values of resistance (e.g. AktioArgostoli) the application of surge arresters reduces the failure rate,
but the reduction is not so impressive. Solution for these lines,
which have more failures due to lightning strikes, could be the use
of arresters with higher energy capability.
Based on the above analysis, it is clear that the application of
surge arresters in overhead high voltage transmission lines can
contribute to the improvement of their lightning performance and
to their failures reduction. Electrical engineer designers must pay
special attention to the surge arresters interval and to their energy
capability, which are related to the tower footing, resistance (dominant parameter) and the ground ash density. For lines with very
low tower footing resistance, surge arresters on every three towers
seems to be almost the same effective, as if they were installed on
every one or two towers, reducing a lot the number of failures. For
transmission lines with high soil resistivity, arresters with higher
withstand capability should be installed on every tower, in order to

achieve better results. In any case, the criteria are not only technical but economical too, since the installation and maintenance cost
and the expected benets have to be estimated for the optimum
and most economic transmission line design.
10. Conclusions
The paper describes in detail a methodology, which assesses
the lightning performance of high voltage transmission lines protected with surge arresters. The total transmission lines failure rate
is assessed including also in the lines faults the surge arresters
failures. The proposed methodology has been applied on six operating Hellenic transmission lines of 150 kV, of known outage rate,
which present signicant different line characteristics through
their length, i.e., ground ash density and tower footing resistance. Three different case studies were analyzed by installing surge
arresters at three different intervals, i.e., arresters were installed
on every tower, on every second tower and on every third tower.
The results have shown that the transmission lines with arresters
present a better lightning performance compared to the lines
without surge arresters. For surge arrester installed every third
tower, with the lines presenting high tower footing resistances,
the reduction of the total line failure rate (backashovers, shielding failures and arresters failures) is not signicant. Results have

182

C.A. Christodoulou et al. / Electric Power Systems Research 80 (2010) 176183

also shown that the arrester failure probability decreases, when

the tower footing resistance and the arresters interval decreases
and the failure rate is dependent on the ground ash density
and the transmission lines length. The proposed methodology can
be used by electric utilities as a useful tool for the design and
lightning protection of electric power systems and especially for
lines, which have signicant different characteristics through their
length, reducing the operational costs and providing continuity of
service.
Acknowledgements
The authors want to express their gratitude to the Hellenic Public Power Corporation S.A. for their kind supply of various technical
data and the Lightning and Atmospheric Electricity Research Team
at the GHCC for the supply of the OTD data.
References
[1] J.M. Clayton, F.S. Young, Estimating lightning performance of transmission
lines, IEEE Transactions on PAS 83 (1964) 11021110.
[2] C. Bouquegneau, M. Dubois, J. Trekat, Probabilistic analysis of lightning performance of high-voltage transmission lines, Electric Power Systems Research
102 (12) (1986) 518.
[3] L.V. Bewley, Travelling Waves on Transmission Systems, 2nd ed., John Wiley &
Sons, Inc., NY, 1951.
[4] F.A.M. Rizk, Modeling of transmission line exposure to direct lightning strokes,
IEEE Transactions on PWRD 5 (4) (1990) 19831997.
[5] IEEE Working Group on Lightning Performance of Transmission Lines, A simplied method for estimating lightning performance of transmission lines, IEEE
Transactions on PAS 104 (4) (1985) 919927.
[6] L.A. Kraft, Modelling lightning performance of transmission systems
using PSpice, IEEE Transactions on Power Systems 6 (2) (1991) 543
549.
[7] J.A. Martinez, F. Castro-Aranda, Lightning performance analysis of overhead
transmission lines using the EMTP, IEEE Transactions on PWRD 20 (3) (2005)
294300.
[8] L. Ekonomou, D.P. Iracleous, I.F. Gonos, I.A. Stathopulos, Lightning performance
identication of high voltage transmission lines using articial neural networks, International Journal of Engineering Intelligent Systems 13 (3) (2005)
219223.
[9] L. Ekonomou, P. Liatsis, I.F. Gonos, I.A. Stathopulos, Articial neural network
based software tool for calculating the lightning performance of high voltage transmission lines, IEE Proceedings: Science, Measurement and Technology
153 (5) (2006) 188193.
[10] L.C. Zanetta, C.E. de Morais Pereira, Application studies of line arresters in partially shielded 138-kV transmission lines, IEEE Transactions on Power Delivery
18 (1) (2003) 95100.
[11] M.S. Savic, Estimation of the surge arrester outage rate caused by lightning overvoltages, IEEE Transactions on Power Delivery 20 (1) (2005) 116
122.
[12] L.C. Zanetta, Evaluation of line surge arrester failure rate for multipulse
lightning stresses, IEEE Transactions on Power Delivery 18 (3) (2003) 796
801.
[13] J.A. Tarchini, W. Gimenez, Line surge arrester selection to improve lightning
performance of transmission lines, in: IEEE Bologna PowerTech Conference,
Bologna, Italy, 2003.
[14] J.A. Martinez, F. Castro-Aranda, Lightning ashover rate of an overhead transmission line protected by surge arresters, in: IEEE Power Engineering Society
General Meeting, 2007, pp. 16.
[15] K. Kannus, K. Lahti, Evaluation of the operational condition and reliability of
surge arresters used on medium voltage networks, IEEE Transactions on Power
Delivery 20 (2) (2005) 745750.
[16] C.H. Shih, R.M. Hayes, D.K. Nichols, R.E. Koch, J.A. Timoshenko, J.G. Anderson, Application of special arresters on 138 kV lines of Appalachian power
company, IEEE Trans Power Apparatus Systems PAS-104 (1985) 2857
2863.
[17] S. Furukawa, O. Usuda, T. Isokaki, T. Irie, Development and application of lightning arresters for transmission lines, IEEE Transactions on Power Delivery 4
(1989) 21212129.
[18] K. Berger, R.B. Anderson, H. Kroninger, Parameters of lightning ashes, Electra
41 (1975) 2337.
[19] I.F. Gonos, L. Ekonomou, F.V. Topalis, I.A. Stathopulos, Probability of backashover in transmission lines due to lightning strokes using Monte-Carlo
simulation, International Journal of Electrical Power and Energy Systems 25
(2) (2003) 107111.
[20] Lightning and Insulator Subcommittee of the T&D Committee, Parameters of
lightning strokes: a review, IEEE Transactions on Power Delivery 20 (1) (2005)
346358.

[21] V. Hinrichsen, Metal-Oxide Surge Arresters, Fundamentals, 1st ed., Siemens,

2001.
[22] IEC 60099-4, Surge ArrestersPart 4: Metal-Oxide Surge Arresters Without
Gaps for a.c. Systems, 2nd ed., 20042005.
[23] ABB, Surge Arresters-Buyers Guide, 5.1 ed., 2007.
[24] L. Montanes, M. Garcia-Gracia, M. Sanz, M.A. Garcia, An improvement
for the selection of surge arresters based on the evaluation of the failure probability, IEEE Transactions on Power Delivery 17 (1) (2002) 123
128.
[25] M. Garcia-Gracia, S. Baldovinos, M. Sanz, L. Montanes, Evaluation of the failure probability for gapless metal oxide arresters, in: IEEE Transmission and
Distribution Conference, 1999.
[26] K. Nakada, S. Yokoyama, T. Yokota, A. Asakawa, T. Kawabata, Analytical study
on prevention for distribution arrester outages caused by winter lightning, IEEE
Transactions on Power Delivery 13 (4) (1998) 13991404.
[27] A.H. Register, A Guide to MATLAB Object-Oriented Programming, Chapman &
Hall/CRC, 2007.
[28] L. Bergeron, Water Hammer in Hydraulics and Wave Surges in Electricity, Wiley,
New York, 1961.
[29] H. Dommel, Digital computer solution of electromagnetic transients in single
and multiple networks, IEEE Transactions on Power Apparatus and Systems
PAS-88 (4) (1969) 388389.
[30] IEEE Working Group 3.4.11, Modeling of metal oxide surge
arresters, IEEE Transactions on Power Delivery 7 (1) (1992) 302
309.
[31] P. Pinceti, M. Giannettoni, A simplied model for zinc oxide surge
arresters, IEEE Transactions on Power Delivery 14 (2) (1999) 393
398.
[32] W.A. Chisholm, Estimates of lightning ground ash density using optical transient density, in: IEEE PES Transmission and Distribution Conference and
Exposition, vol. 3, 2003, pp. 10681071.
[33] H.J. Christian, R.J. Blakeslee, D.J. Boccippio, W.L. Boeck, D.E. Buechler, K.T.
Driscoll, S.J. Goodman, J.M. Hall, W.J. Koshak, D.M. Mash, M.F. Stewart, Global
frequency and distribution of lightning as observed from space by the optical transient detector, Journal Geophysics Research 108 (D1) (2003) 4005
4019.
[34] NASA OTD Lightning, 2000. http://thunder.nsstc.nasa.gov.
[35] IEEE Working Group on Estimating the Lightning Performance of Transmission Lines, Estimating the lightning performance of transmission lines
IIupdates to analytical models, IEEE Transactions on PWRD 8 (3) (1993) 1254
1267.
[36] L. Ekonomou, I.F. Gonos, I.A. Stathopulos, Lightning performance assessment
of Hellenic high voltage transmission lines, Electrical Power Systems Research
78 (4) (2008) 703712.
[37] L. Ekonomou, I.F. Gonos, I.A. Stathopulos, E.N. Dialynas, A backashover model
for calculating the transmission lines lightning performance, vol. 1 (26), in:
ISH2001, Bangalore, India, 2001.
[38] L. Ekonomou, I.F. Gonos, I.A. Stathopulos, An improved backashover model for
estimating the lightning performance of transmission lines, in: UPEC, Thessaloniki, Greece, 2003, pp. 109112.
[39] PPC S.A., Transmission Lines Characteristics, Hellenic Public Power Corporation
S.A., Athens, 2006.

Christos A. Christodoulou was born on June 9, 1983 in

Lamia, Greece. He received his diploma in Electrical
and Computer Engineering from the National Technical
University of Athens in 2006. He is currently a Ph.D.
candidate in the N.T.U.A.s High Voltage Laboratory. His
research interests concern high voltages, lightning protection, distribution and transmission lines, measurements
and uncertainties, and electromagnetic compatibility. He
is member of the Technical Chamber of Greece and
the Hellenic Association of Electrical and Mechanical
Engineers.

Lambros Ekonomou was born on January 9, 1976 in Athens,

Greece. He received a Bachelor of Engineering (Hons)
in Electrical Engineering and Electronics in 1997 and a
Master of Science in Advanced Control in 1998 from University of Manchester Institute of Science and Technology
(U.M.I.S.T.) in United Kingdom. In 2006 he graduated with
a Ph.D. in High Voltage Engineering from the National
Technical University of Athens (N.T.U.A.) in Greece. He
has worked with various companies including Hellenic
Public Power Corporation S.A. and Hellenic Aerospace
Industry S.A., while in 2008 he was appointed Assistant Professor in A.S.PE.T.E. School of Pedagogical and
Technological Education. His research interests concern high voltage engineering, transmission and distribution lines, lightning performance and protection,
stability analysis, control design, reliability, electrical drives and articial neural
networks.

C.A. Christodoulou et al. / Electric Power Systems Research 80 (2010) 176183

Georgios P. Fotis was born on July 26, 1977 in Athens,

Greece. He received his diploma and his Ph.D. in Electrical Engineering from the National Technical University
of Athens in 2001 and 2006, respectively. He is currently a research associate in the N.T.U.A.s High Voltage
Laboratory. His research interests concern high voltages, electromagnetic compatibility and electrostatic
discharges. He is member of IEEE and the Technical Chamber of Greece. Dr. Fotis is the author of 25 papers in
scientic journals and conferences proceedings.
Ioannis F. Gonos was born on May 8, 1970 in Artemisio, Arcadia, Greece. He received his diploma in Electrical
Engineering and his Ph.D. from the National Technical University of Athens in 1993 and 2002, respectively. He was a
teaching assistant at the Hellenic Naval Academy and the
Technological Education Institute of Athens (19962001).
He is working at the High Voltage Laboratory of N.T.U.A.
(since 2001). His research interests concern grounding
systems, insulators, high voltages, measurements and
genetic algorithms. He is the author of more than 60
papers in scientic journals and conferences proceedings.

183

Ioannis A. Stathopulos was born in Volos, Greece in 1951.

He studied in the Faculty of Electrical and Mechanical Engineering of the National Technical University of
Athens (19691974). He carried out his doctor thesis
at the Technical University of Munich (19741978). He
become teaching assistant at the Technical University of
Munich (19741978), production engineer in the company VianoxFranke (19791980), teaching assistant at
the National Technical University of Athens (19791983)
and thereafter Lecturer (19831987), Assistant Professor
(19871991), Associate Professor (19911995) and Professor (since 1995) in the High Voltage Laboratory of the
N.T.U.A. He is the author of 8 books and more than 90 papers in scientic journals and conferences proceedings. He is lead assessor of the Hellenic Accreditation
Council.