Including Surge Arresters in the Lightning

Performance Analysis of 132kV Transmission

Line
Saeed Mohajeryami, Milad Doostan
University of North Carolina at Charlotte
Department of Electrical and Computer Engineering
Charlotte, NC
{smohajer,mdoostan}@uncc.com

Abstract- Line arresters are considered as an effective reasonable price and without tower structure
way to improve the lightning performance of reinforcement [4].
transmission lines, especially in parts of line that suffer Authors in [5] have investigated the effectiveness of
from high soil resistivity and lightning ground flash external and internal surge arresters on 500kV DC gas
density. This paper presents results of the application of insulated lines. In [6] also it has been shown that the
line surge arresters on the 132KV double circuit surge arresters are effective in controlling the
transmission line in EMTP-RV and all the practical switching overvoltages. The authors in [7] have
scenarios for installation of surge arresters. The study investigated the surge arrester application for 750kV
has shown that a significant level of improvement can be transmission lines and have obtained the right size. To
reached by installing arresters at all or only some of the assess the effectiveness of surge arresters, a
line phases. It can increase the strength of the line to probabilistic approach has been adopted by [8] to
withstand lightning currents up to -292kA. the estimate the failure risk of such surge arresters.
probability of having this lightning current, is practically This paper presents results of the application of
zero. transmission line surge arresters on the 132kV double
Index Terms—Surge Arrester, Backflashover,
circuit transmission line with using accurate models of
Lightning Overvoltage, Volt-time model, EMTP-RV line components in EMTP-RV environment.
At first, a brief description of test line has been
I. INTRODUCTION
presented and then, models used for modeling
lightning source, tower and transmission line, tower
Power equipment failure can incur huge financial footing resistance, insulator gap and surge arrester in
losses to the transmission companies and it has EMTP-RV have been introduced and their parameters
necessitated the need for paying attention to these have been calculated and finally simulation results and
widely exposed lines which are subject to vast variety their characterization and conclusion have been
of faults. One of the most frequent issues occurring at presented.
the transmission lines are outages caused by the II. BRIEF DESCRIPTION OF TEST LINE
lightning strikes [1]. These outages can be decreased
by the proper insulation design and the design of Arvandkenar-Abadan 132kV double circuit is a
varied devices hinges on lightning overvoltage in the transmission line located in Iran which has been used
network [2]. as a case study. The complete transmission line
To decrease the outage rate due to the lightning, information is tabulated in table I. Fig. 1 shows tower
many methods have been proposed such as tower structure, phase and shield wire configuration.
footing resistance reduction, increase of line insulation III. MODELING
level, installation of additional ground and guy wires, A. Lightning current source model
addition of under-built ground wires, etc. limited
In order to calculate lightning outage rate, the
effect or being expensive or difficult to use, are some
probability of exceeding stroke of current, I, is needed.
of the reasons that make some of them impractical [3].
This probability distribution function can
Thanks to the development of polymer-housed
approximately be computed from (2) [9].
metal-Oxide surge arresters, it became practical to use
arresters to control line lightning performance, with
 TABLE I
TRANSMISSION LINE DATA
General Data Shield Wire Data
Voltage Level 132 kV Core Type Hawk
Line Length 53 km Number 1 per tower
Location Mahshahr, Iran Diameter 8.04 mm
Number of Thunderstorms 15.2 day in year Cross Section Area 39.46mm2
DC Resistance 2.9 Ω/Km
Conductor Data Tower Data
Type Hawk/AL,26/3.87,477MCM Type HS2-10
Number of Bundles 2 Insulator Gap Length 1350mm
Bundle Spacing 45.7 cm Footing Resistance 20 ohm maximum
Diameter 21.8 mm Soil Resistivity 100-150 Ω.meter
Total Area 281.03 mm2 Lightning Discharge Path
400Ω
DC Resistance 0.1199 Ω/Km Impedance
GMR 0.8245cm Ruling Span 350m

1 B. Tower and Transmission line model

P1 = (2) Double circuit transmission line is modeled with
1
1+ ( ) 2.6 frequency-dependent model in EMTP-RV. Phase
31 coupling is considered in this model. Multistory model
[11] is used for tower model. It is multistoried as
In this paper, double exponential formula which is shown in Fig. 2. Tower surge impedance is calculated
provided in (3) has been employed to model the from CIGRE recommended equation (6) and then
current source [10]. because of similar structure of HS2-10” vertical tower
t t
and [11], the same ratio has been chosen for dividing
− − surge impedance between upper and lower halves.
=i (t ) 1.04 I m (e T1
−e T2
) (3) R
T1 = 1.36543TR (4) Z surge = 60 ln[cot(0.5 tan −1 ( ))] (6)
h
Ts Where
T2 = (5) r1h2 + r2 h + r3 h1
2.282835 R= , h= h1 + h2 (7)
Where 𝑇𝑇𝑠𝑠 is “rise time” and 𝑇𝑇𝑅𝑅 is “tail time”. h
Indeed, as the positive polarity lightning stroke is Height and radius has been shown in Fig. 3.
Calculated tower model parameters and their formulas
nearly 5% of all strokes, it is neglected in this work.
have been listed in Table II.
Moreover, the impedance of the lightning discharge
channel is assumed to be 400Ω.

Fig. 2 Multistory tower model Fig. 3 Tower schematic

C. Tower footing resistance model

Several methods have been proposed to calculate
the tower footing resistance while there are fast
transient surges which can ionize the soil and change
the resistance. However, one of the most simple and
practical approaches which takes the surges
Fig. 1 HS2-10” type tower structure considerations into account is provided by CIGRE [9].
 TABLE II. Fig. 4 shows the V-t relationship of a sample
MULTISTORY TOWER MODEL PARAMETERS
𝑍𝑍𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠 = 179.2 [Ω] insulator.
𝑍𝑍𝑡𝑡1 = 𝑍𝑍𝑡𝑡2 = 𝑍𝑍𝑡𝑡3 = 1.19𝑍𝑍 = 213.2 [Ω]
𝑍𝑍𝑡𝑡4 = 0.81𝑍𝑍 = 145.15 [Ω]
𝑣𝑣𝑡𝑡1 = 𝑣𝑣𝑡𝑡2 = 𝑣𝑣𝑡𝑡3 = 𝑣𝑣𝑡𝑡4 = 𝑣𝑣 = 300 [m/𝜇𝜇𝜇𝜇]
𝛾𝛾 = 0.8944
𝜏𝜏 = 2𝐻𝐻⁄𝑣𝑣=0.28 [𝜇𝜇𝜇𝜇]
𝑟𝑟1 = −(2 × 𝑍𝑍𝑡𝑡1 × ln 𝛾𝛾)/(ℎ1 + ℎ2 + ℎ3 ) = 2.5 [Ω/m]
𝑟𝑟2 = −(2 × 𝑍𝑍𝑡𝑡4 × ln 𝛾𝛾)/ℎ4 = 1.25 [Ω/m]
𝑅𝑅1 = 𝑟𝑟1 × ℎ1 = 14 [Ω]
𝑅𝑅2 = 𝑟𝑟1 × ℎ2 = 16.8 [Ω]
𝑅𝑅3 = 𝑟𝑟1 × ℎ3 = 16.8 [Ω]
𝑅𝑅4 = 𝑟𝑟2 × ℎ4 = 32.4 [Ω]
𝐿𝐿1 = 𝑅𝑅1 × 𝜏𝜏 = 3.92 [𝜇𝜇𝜇𝜇]
𝐿𝐿2 = 𝑅𝑅2 × 𝜏𝜏 = 4.07 [𝜇𝜇𝜇𝜇]
𝐿𝐿3 = 𝑅𝑅3 × 𝜏𝜏 = 4.07 [𝜇𝜇𝜇𝜇]
𝐿𝐿4 = 𝑅𝑅4 × 𝜏𝜏 = 9.1 [𝜇𝜇𝜇𝜇] Fig. 4 Volt-time model
E. surge arrester model
According to this method the tower footing Metal Oxide surge arresters protect power
resistance could be calculated from (8). equipment insulation from system overvoltage. These
R0 arresters have high resistance at the normal system
RF = (8)
I operation but when facing system overvoltage, they
1+
Ig will show low resistance. There are several models
proposed in literature to describe this nonlinear
Where: 𝑅𝑅0: Tower footing resistance at low current
behavior but in this study IEEE frequency dependent
and low frequency (ohm), 𝐼𝐼𝑔𝑔 : Limiting current model [13] has been chosen because of its enough
initiating soil ionization (kA) which can be calculated accuracy for lightning studies.
by (9):
1 E0 ⋅ ρ
=
Ig × (9)
2π R0 2
Where 𝜌𝜌 is soil resistivity (ohm-meter), and 𝐸𝐸0 is soil
ionization gradient (about 400 kV/m).
Fig. 5 IEEE frequency dependent model
D. Backflashover model
The inductance L0 in the model represents the
In volt-time model [9], it is assumed that inductance associated with magnetic fields near the
Backflashover occurs when tower voltage is higher arrester. The resistor R0 is used to stabilize the
than lightning impulse withstand voltage or Basic numerical integration in the software. The capacitance
insulation level (BIL) of the insulator strings. The C represents the terminal-to-terminal capacitance of
lightning impulse withstand voltage of the insulator the arrester.
string, is not a unique number. The insulator string The inductance L1 and the resistance R1 make a
may withstand a high magnitude impulse voltage with filter between the two nonlinear resistances.
short duration while it fails to withstand a lower For slow-front surges, this R-L filter has very little
magnitude impulse voltage with longer duration. A impedance and the two non-linear sections of the
simplified expression of withstand voltage capability model are in parallel. For fast-front surges the
for an insulator string can be calculated as in (10): impedance of the R-L filter becomes more significant.
𝐾𝐾
𝑉𝑉𝑓𝑓0 = 𝐾𝐾1 + 2� 0.75 (10) This results in more current in the non-linear section
𝑡𝑡
Where Vfo = flashover voltage (kV), named A0 than in the section named A1.
K1 = 400*L d 
L1 = 15   µ H (11)
K2 = 710*L n
L = insulator length (m) d 
=R1 65   Ω (12)
t = elapsed time after lightning stroke, μs. n
 Fig. 7 Simulated line schematic

Fig. 8 (upper left) insulator induced voltage without arrester, Fig. 9 (upper middle) insulator induced voltage with arrester, Fig. 10 (upper right) installation
scenarios, Fig. 11 (middle left) Scenarios with one arrester, Fig. 12 (middle middle) Scenarios with two arresters, Fig. 13 (middle right) Scenarios with three
arresters, Fig. 14 (lower left) Scenarios with four arresters, Fig. 15 (lower middle) Scenarios with five arresters, Fig. 16 (lower right) Lightning current
flashover in all scenarios

d 
L0 = 0.2   µ H (13)
n
d 
=R0 100   Ω (14)
n
n
C = 100   pF (15)
d 
Where d is the estimated height of the arrester in
meters (from catalog data) and n is the number of
parallel columns of metal oxide in the arrester
The non-linear V-I characteristics A0 and A1 can be
estimated from per unitized curves given in Figure 6. Fig. 6 V- I Relationships for Non-Linear Resistor Part of Model
For A0 in Fig. 6, The study has shown that, as expected, a significant
level of improvement can be reached by installing
𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷ℎ𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎, 𝐾𝐾𝐾𝐾 = [𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅 𝐼𝐼𝐼𝐼 𝑖𝑖𝑖𝑖 𝑝𝑝. 𝑢𝑢. 𝑓𝑓𝑓𝑓𝑓𝑓 𝐴𝐴0 (𝑖𝑖)]
∗ [𝑉𝑉10 ⁄1.6] (16) arresters at all or only some of the line phases. The
Likewise for A1, the improvement of lightning performance can be very
𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷ℎ𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎, 𝐾𝐾𝐾𝐾 = [𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅 𝐼𝐼𝐼𝐼 𝑖𝑖𝑖𝑖 𝑝𝑝. 𝑢𝑢. 𝑓𝑓𝑓𝑓𝑓𝑓 𝐴𝐴1 (𝑖𝑖)] significant when arresters are installed at two phases,
∗ [𝑉𝑉10 ⁄1.6] (17) but even with the installation of a single arrester per
Where V10 is discharge voltage for 10kA current with tower at the upper phase, a significant reduction of the
8/20 µs waveshape per kV. Table III shows line total flashover rate can be achieved.
arrester specification used for this simulation.
TABLE III REFERENCES
SIEMENS 3EL5 108-0LK23 SPECIFICATION [12] [1] S.T. Mobarakei, T. Sami, B. Porkar, "Back Flashover
Height(d) 1.53 m phenomenon analysis in power transmission substation for
Parallel column(n) 1 insulation coordination," 11th International Conference on
Maximum value of Residual voltage at Environment and Electrical Engineering, pp.170-174, May 2012
308kV [2] M. Babaei, M. Babaei, M. Niasati, "Parametric analysis of
discharge current 8/20 µs and 10 kA(V10)
overvoltages caused by back-flashover in “Siah-bishe” 400kV GIS
Maximum value of Residual voltage at
259kV substation," 3rd International Conference on Electric Power and
discharge current 30/60 µs and 2kA
Energy Conversion Systems, pp.1-6, Oct. 2013
Maximum Continuous Operation Voltage 86 kV [3] K. Munukutla et al, “A Practical Evaluation of Surge Arrester
Placement for Transmission Line Lightning Protection”, IEEE
IV. SIMULATION AND RESULTS Trans. on Power Delivery, Vol. 25, No. 3, July 2010
[4] J.A. Martinez, F. Castro-Aranda,” Modeling Overhead
Simulation has been done with accurate modeling Transmission Lines for Line Arrester Studies”, IEEE Power
of components in EMTP-RV software. Fig. 7 shows Engineering Society General Meeting, pp.1125-1130 June 2004
the schematic of test line. Fig. 8 shows the induced [5] F. Goll, R. Witzmann, "Lightning Protection of 500-kV DC Gas-
Insulated Lines (GIL) With Integrated Surge Arresters," IEEE
voltage on insulator without surge arrester for standard Trans. on Power Delivery, vol.30, no.3, pp.1602-1610, June 2015
lightning current of -100 kA that as it can be observed, [6] P. Bunov, et al "Transmission line arresters application for
flashover occurs across insulator. Fig. 9 shows the control of switching overvoltages on 500-kV transmission line,"
same situation for insulator but with surge arrester IEEE T&D Conference and Exposition, pp.1-5, April 2014
[7] J. Wang, et al "The research and design of parameters for 750kV
installed at the upper phase. For installation of surge transmission line metal-oxide surge arresters with series gap,"
arresters a comprehensive plan must be studied to International Conference on High Voltage Engineering and
consider all scenarios. Fig. 10 illustrates all the Application (ICHVE), pp.1-4, Sept. 2014
different practical scenarios for installation of [8] R. Shariatinasab, F. Ajri, H. Daman-Khorshid, "Probabilistic
evaluation of failure risk of transmission line surge arresters caused
insulators. Scenarios with the same number of surge by lightning flash," IET Generation, Transmission & Distribution,
arrester per tower have been compared in Figs. 11-15. vol.8, no.2, pp.193-202, February 2014
Finally, in Fig. 16 lightning current causing flashover [9] CIGRE WG 33-01: “Guide to Procedures for Estimating the
in all scenarios, have been shown graphically. Lightning Performance of Transmission Lines”, Technical
Brochure, October 1991.
V. CONCLUSION [10] B.Vahidi, M.Ghorat, E.Goudarzi,”Overvoltage calculation on
According to the simulation results, the standard Bam substation by monte carlo method with accurate substation
components modeling”, IEEE conference on power technology,
lightning current of -77kA can cause flashover across pp.1807-1811, July 2007
insulators of the line without surge arrester. This [11] M. Ishii, T. Kawamura, T. Kouno, E. Ohsaki, K. Muromi,
current results 3.7 line outage per 100km per year, but T.Higuchi, “Multistory Transmission Tower Model for Lightning
with installing surge arrester, this line can withstand surge Analysis”, IEEE Trans. on Power Delivery, Vol. 6, No. 3,
pp.1327-1335, July 1991
lightning currents up to -292kA. The probability of [12] SIEMENS surge arrester 3EL, Germany, 2010.
having this lightning current, is practically zero. [13] IEEEWorking Group 3.4.11, “Modeling of metal oxide surge
As shown in Fig. 16, scenarios 1, 4 and 10 are the arrester,” IEEE Trans. on Power Del., vol. 7, no. 1, pp. 302–309,
most efficient among other installation scenarios; as it Jan. 1992.

is shown in Figs. 11-16, inefficient strategies for

installing surge arresters can waste a lot of money even
with more arresters per tower.