Simulation of Transients with a Modal-Domain

Based Transmission Line Model Considering Ground

as a Dispersive Medium
Alberto De Conti and Maique Paulo S. Emídio

Abstract-- In this paper, a modal-domain based transmission of grounding electrodes [7-11] and on the calculation of
line model available in popular electromagnetic transient switching and lightning transients in power systems [12-19].
programs is adapted to evaluate the effect of frequency-dependent In spite of this trend, the transmission line models available
ground conductivity and permittivity in the calculation of
in popular transient simulators rely on the use of Carson’s
transients on overhead power distribution lines. The calculation
of the line parameters considering ground as a dispersive medium integrals [20] or their approximation through the use of the
is performed in MATLAB making use of practical equations that complex ground return plane [21] for calculating the ground
are based on in situ measurements of ground conductivity and return impedance. Since both approaches are low-frequency
permittivity in a wide frequency range. The propagation function approximations in which σ>>ωε is assumed, the influence of
and characteristic impedance of the line are synthesized in the the ground permittivity in the ground-return impedance is not
frequency domain as the sum of rational functions using the taken into account properly. This feature, combined with the
vector fitting technique. The poles and residues of the synthesized
assumption of a constant value of ground conductivity,
functions are written in a .pch file that is read by the Alternative
Transients Program (ATP) as a JMarti model. Time domain suggests the possibility of errors in the simulation of cases
simulations are performed considering both switching and involving poor ground conductivities, high-frequency
lightning transients on single- and two-phase power distribution transients, or a combination of both.
lines. The results indicate that the consideration of the ground as In this paper, the transmission line model proposed by
a dispersive medium leads to a distortion of the calculated Marti [22] is adapted to include the effect of frequency
transient voltages that can be relevant in the case of very low dependent ground parameters in the time-domain simulation of
ground conductivity. However, such effect can be significantly
electromagnetic transients in power distribution lines. The
reduced in line topologies that include multiple branches and
grounding points. It is also shown that constant values of ground calculation of the line parameters is implemented in
conductivity and permittivity are able to lead to results MATLAB, where the poles and residues necessary to fit the
comparable to those obtained with frequency dependent characteristic impedance and propagation function of each
conductivity and permittivity provided a suitable value of ground mode are determined in the frequency domain via the vector
relative permittivity is selected. fitting technique [23]. The obtained poles and residues are
Keywords: Transmission line modeling, modal domain, then written in a .pch file that can be read by the Alternative
frequency-dependent ground parameters.
Transients Program (ATP) and solved as a JMarti model.
This paper is organized as follows. Section II discusses the
I. INTRODUCTION
calculation of the ground return impedance of transmission

T HERE has been an increasing interest in the simulation

of electromagnetic transients in transmission lines
considering the dispersive nature of the ground parameters.
lines considering different expressions under the assumption of
constant or frequency-dependent ground parameters. Section
III discusses the use of the J. Marti model available in ATP for
This is apparent not only from the number of publications calculating transients on power distribution lines considering
describing new methodologies to measure and model the ground as a dispersive medium. Results and analysis are
variation with frequency of the ground conductivity (σ) and presented in Section IV, followed by conclusions in Section V.
permittivity (ε) presented in the last fifteen years or so [1-6],
but also from the increasing amount of papers discussing the II. CALCULATION OF THE GROUND RETURN IMPEDANCE OF
prospective effect of such variation on the transient response TRANSMISSION LINES UNDER DIFFERENT ASSUMPTIONS

This work was supported by Brazilian agencies CAPES (Coordenação de A. Ground Return Impedance
Aperfeiçoamento de Pessoal de Nível Superior) and CNPq (National Council The self and mutual terms of the ground impedance matrix
for Scientific and Technological Development).
Alberto De Conti and Maique Paulo S. Emídio are with the Lightning
of overhead transmission lines can be calculated with the
Research Center, Federal University of Minas Gerais - UFMG, Belo equations proposed by Sunde [24]. Taking as reference the
Horizonte – MG, Brazil (conti@cpdee.ufmg.br; maiquepaulo@hotmail.com). two-conductor line illustrated in Fig. 1, Sunde’s equations read

Paper submitted to the International Conference on Power Systems

Transients (IPST2015) in Cavtat, Croatia June 15-18, 2015
 jωµ 0 ∞ e −2 hi λ can be shown that Carson’s integrals (6) and (7) are accurate
Z g′ ii =
π ∫ 0
λ2 + γ g2 + λ
dλ (1) up to few MHz, which is the frequency range of most power
system transients. However, for a relatively poor ground
− ( hi + h j ) λ conductivity of σ=0.001 S/m inaccuracies are observed in the
jωµ 0 ∞ e
Z ′g ij =
π ∫ 0
λ2 + γ g2 + λ
cos(rij λ )dλ (2) phase angle of Carson’s ground return impedance in
frequencies of few hundreds kHz [25]. This is because for such
where frequencies and such value of ground resistivity the
assumption σ>>ωε may no longer hold. In a worst case
scenario involving poorer ground conductivities, the validity
γg = jωµ 0 (σ + jωε r ε 0 ) (3)
of Carson’s integrals is even more limited. This is shown in
Fig. 4, which illustrates the percentage error in the magnitude
in which µ0=4π×10−7 H/m, ε0=8.85×10−12 F/m, ω is the angular and phase angle of the ground impedance calculated with (6)
frequency in rad/s, and σ and εr are the ground conductivity for σ=0.0001 S/m, taking as reference Sunde’s expression (1)
and relative permittivity, respectively. The following for hi=10 m, r=5 mm, and different values of ground relative
approximations to (1) and (2) in logarithmic form have been permittivity. It is seen Fig. 4(a) that the percentage error in the
proposed by Sunde [24] and Rachidi et al. [25], respectively magnitude and phase angle of the ground impedance can be
very significant for frequencies above 100 kHz.
jωµ 0 1 + γ g hi  A simplified way to obtain the ground return impedance
Z ii′ , g = ln   (4) consists in using the well-known approximate formulas
2π  γ g hi  proposed by Deri et al. [21], which leads to

jωµ 0 [1 + 0.5γ g ( hi + h j )] + (0.5γ g rij ) 

2 2
Z ij′ , g = ln  jωµ 0  hi + p& 
2 
(5)
4π 2 Z g′ ii = ln   (8)
 [0.5γ g (hi + h j )] + (0.5γ g rij )  2π  hi 

In [25] it is shown that (4) and (5) reproduce (1) and (2) for  2 2 
0.01 ≥ σ ≥ 0.001 S/m with good accuracy. In Figs. 2 and 3, jωµ 0  (hi + h j + 2 p& ) + rij 
Z ′g ij = ln (9)
2π  
assuming two conductors with 5 mm radius, one located 10 m  (hi + h j ) 2 + rij2 
and the other 8.17 m above the ground, with rij=2 m (see Fig.
1), it is confirmed that (4) and (5) are also sufficiently accurate
If σ>>ωε, the complex penetration depth reads [21]
to reproduce (1) and (2) for σ=0.0001 S/m and different values
of εr, at least for heights and distances between conductors that
are typical of power distribution line configurations. 1
p& = (10)
If σ>>ωεrε0 is assumed in (3), then (1) and (2) reduce to jωµσ

jωµ 0 ∞ e −2 hi λ 3
Z g′ ii = ∫ dλ (6) Integral Expression (1)
Ω /m)

σ(ω), εr(ω)
2.5
π 0 2
λ + jωµ 0σ + λ Logarithmic Approx. (4) εr=10
Magnitude (Ω

2
1.5
− ( hi + h j ) λ
jωµ 0 ∞ e 1
Z ′g ij =
π ∫ 0
λ2 + jωµ 0σ + λ
cos(rij λ )dλ (7) 0.5 (a)
εr=40
εr=20
0
1.E+03 1.E+04 1.E+05 1.E+06 1.E+07 1.E+08

which are the equations proposed by Carson to represent the Frequency (Hz)

ground-return impedance of overhead lines [20]. They are

Phase Angle (Degrees)

80
implemented in popular transient simulators and for this reason Integral Expression (1)
60 Logarithmic Approx. (4)
are often used in the simulation of transients in power systems.
εr=10
40
εr=40
σ(ω), εr(ω)
dij 20 εr=20
j (b)
i 0
hi rij hj 1.E+03 1.E+04 1.E+05 1.E+06 1.E+07 1.E+08
Frequency (Hz)
Ground (σ, ε) Fig. 2 – Comparison between Sunde’s integral expression (1) and its
logarithmic approximation (4) for constant σ=0.0001 S/m and different
Fig. 1. Problem geometry. values of εr. Also included are curves obtained for frequency-dependent
values of ground resistivity and permittivity according to (12) and (13),
For a relatively high ground conductivity of σ=0.01 S/m, it obtained for ρ0=10000 Ω.m.
 3 simplified equations proposed by Deri et al. [21] to calculate
Integral Expression (2) σ(ω), εr(ω)
Ω /m)

2.5 Logarithmic Approx. (5) εr=10 the ground-return impedance. This procedure can be readily
Magnitude (Ω

2 used to include the effect of the ground permittivity in the

1.5 calculation of the series impedance of transmission lines, as
1
εr=20 done in, e.g., [12, 25, 29].
0.5 (a) εr=40 It is to be noted that the lack of explicit representation of
0
1.E+03 1.E+04 1.E+05 1.E+06 1.E+07 1.E+08 the ground permittivity in (6) and (7) can be related either to
Frequency (Hz) assuming εr=0 (if the current propagation in the wire is
neglected in the derivation of equations (6) and (7)) or to
Phase Angle (Degrees)

80
Integral Expression (2) assuming εr=1 (if a lossless propagation constant is considered
60 Logarithmic Approx. (5) in the derivation of the aforementioned equations) [26-28]. For
40 εr=10 example, it is known that Sunde’s equations (1) and (2) were
εr=40
σ(ω), εr(ω) derived neglecting propagation effects, which could be viewed
20
(b) εr=20 as a shortcoming. However, for including propagation effects
0 in (1) and (2) as well as in their approximate representations
1.E+03 1.E+04 1.E+05 1.E+06 1.E+07 1.E+08
(4) and (5), it suffices to use the product (εr – 1)ε0 instead of
Frequency (Hz) εrε0 in (3) and (11), as suggested in [26, 28]. Given that the
Fig. 3 – Same as Fig. 2, but for Sunde’s integral expression (2) and its differences observed in the calculated ground return
logarithmic approximation (5).
impedance with assuming either (εr – 1)ε0 or εrε0 in (3) are not
significant for the transient studies performed in this paper, all
250
Percentage Error (%)

(a) Magnitude εr=20 calculations presented here consider Sunde’s equations in their
200
εr=40 original form, which means to assume εrε0 in (3) and (11).
150
100 B. Ground Admittance
σ(ω), εr(ω) εr=10
50 Several expressions have been proposed to include the
0 ground admittance in the calculation of transmission line
1.E+03 1.E+04 1.E+05 1.E+06 1.E+07 parameters [26]. It can be shown that at high frequencies the
Frequency (Hz) ground admittance can be used to explain the transition from a
250 pure TEM propagation to a mixed TEM/TM/TE propagation
Percentage Error (%)

(b) Phase Angle

200 εr=20 [26]. However, different authors have come to the conclusion
150
εr=40 that for typical frequencies associated with power system
σ(ω), εr(ω) transients and realistic values of ground conductivity and
100
εr=10 permittivity the effect of the ground admittance can be
50
neglected without significant errors [13, 18-19]. For this
0 reason, only the ground impedance is assumed in this paper to
1.E+03 1.E+04 1.E+05 1.E+06 1.E+07
be affected by the non-perfectly conducting ground.
Frequency (Hz)
Fig. 4 – Percentage error in the (a) magnitude and (b) phase angle of the self C. Frequency-Dependent Ground Parameters
term of the ground impedance calculated with Carson’s expression (6), taking
The analysis in Section II-A assumes both σ and ε as
as reference the more general formula of Sunde (1) for σ=0.0001 S/m. Also
included is the error curve obtained considering frequency-dependent values constants. However, it is known that both parameters present a
of ground resistivity and permittivity according to (12) and (13), calculated wide variation with frequency [1-6]. Therefore, it is expected
for ρ0=10000 Ω.m. that an accurate simulation of electromagnetic transients in
transmission lines should consider the ground as a dispersive
In this particular case, the ground return impedance given medium, in which the conductivity and permittivity vary with
by (8) and (9) reproduces with good accuracy Carson’s frequency. In fact, some authors suggest that the variation of σ
formulas (6) and (7) [20]. In the more general case considering
and ε with frequency can affect switching and lightning
the explicit representation of the ground permittivity, the
transients in frequencies as low as tens of kHz, depending on
complex penetration depth is given by [29]
the soil characteristics [2].
In recent years, different procedures have been proposed
1 1 for measuring and modeling the variation of σ and ε with
p& = = (11)
γg jωµ 0 (σ + jωε r ε 0 ) frequency [1-6]. In this paper, the soil model of Visacro and
Alipio [4] is considered in the calculation of the ground-return
in which case (8) and (9) are equivalent to (4) and (5). In other impedance with (4) and (5). This soil model is based on
words, the errors observed in Fig. 4 with the direct application measurements performed at 31 different sites in Brazil, where
of Carson’s formula (6) become negligible if the complex low-frequency resistivity values ranging from 60 to 9100 Ω.m
penetration depth given by (11) is used instead of (10) in the were recorded. In all cases, a strongly frequency-dependent
behavior was observed for both σ and ε. dispersive medium in the time domain simulation of
By defining the relative resistivity as ρr(ω)=ρ(ω)/ρ0, where electromagnetic transients. For such, the vector fitting
ρ(ω) is the ground resistivity and ρ0 is the ground resistivity at technique [23] was used for synthesizing the characteristic
100 Hz, Visacro and Alipio [4] proposed the following impedance and the propagation function of the evaluated lines.
approximate expressions to represent the frequency-dependent For calculating the line parameters, the complex penetration
behavior of ρr(ω) and εr(ω) depth given by (11) was used in the expressions proposed by
Deri et al. [21], which is equivalent to considering Sunde’s
ρ r (ω ) = {1 + [1.2 ×10 −6 ρ 00.73 ] ⋅ [( f − 100) 0.65 ]}−1 (12) formulas (4) and (5). The variation of σ and ε with frequency
is assumed to be governed by the soil model described in
 7.6 ×103 f −0.4 + 1.3 f ≥ 10 kHz Section II-C.
ε r (ω ) =  (13) A dedicated set of poles was used to represent each
ε r (ω ) f =10 kHz
f < 10 kHz
transmission line mode. This was necessary because for a
poorly conductive ground the propagation function of the
Equations (12) and (13) are able to represent with very ground mode attenuates at a much lower frequency than the
good accuracy the behavior of the evaluated soil samples in aerial modes. The lossless time delay associated with each
the frequency range 100 Hz - 4 MHz [4]. Their use in the mode was used in the synthesis, which typically comprised a
simulation of the transient response of grounding electrodes frequency range from 0.1 Hz to 107-109 Hz, depending on the
also leads to a better agreement between the predictions of a considered line length. The real transformation matrix
rigorous electromagnetic model with measured data [7-9]. This necessary for the time domain simulations was calculated at
gives confidence about the accuracy and suitability of (12) and the upper frequency of the assumed frequency range.
(13) to evaluate the frequency dependence of σ and ε in the To include the effect of frequency dependent ground
simulation of transients in transmission lines. parameters in Marti’s model, real poles and residues calculated
Application examples of (12) and (13) are illustrated in the in MATLAB with the vector fitting technique were written as a
curves labeled as σ(ω) and εr(ω) shown in Figs. 2, 3, and 4, .pch file compatible with ATPDraw. Transient voltages and
where a ground conductivity σ0(ω)=1/ρ0 at low frequencies currents were then calculated at each time step by the ATP
was assumed for ρ0=10000 Ω.m. It is seen that the solver. A full code with a version of Marti’s model extended to
consideration of frequency-dependent ground parameters deal with complex poles was also written in MATLAB. The
affects considerably both the magnitude and phase angle of the implemented code was used to double check the results
ground-return impedance. obtained with ATP. In all cases discussed in Section IV, the
results obtained via the .pch files simulated in ATP were seen
III. TRANSMISSION LINE MODEL OF MARTI CONSIDERING to lead to results identical to those obtained by the MATLAB
DISPERSIVE GROUND PARAMETERS code using complex poles. For this reason, only the results
The transmission line model proposed by Marti [22] is obtained with the .pch files are presented.
possibly the most popular model for the digital simulation of
electromagnetic transients on overhead lines. It is a IV. RESULTS AND ANALYSIS
distributed-parameter model in which the variation of the line This section presents results of simulations of switching and
parameters with frequency is automatically considered in a lightning transients on overhead lines considering either
frequency range determined by the user. The solution of the frequency-dependent or constant ground parameters for
transmission line equations is performed in the modal domain, different values of εr. Conductor heights and distances typical
where a system of n coupled conductors is represented as n of power distribution lines are considered.
independent single-phase lines by means of a similarity
transformation. For the computation of the voltages and A. Single-phase distribution line
currents in time domain, a constant and real transformation Fig. 5 illustrates voltages calculated at the receiving end of
matrix calculated at a frequency determined by the user is a 10-m high, single-phase overhead line with radius of 5 mm
considered [22]. subjected to a lightning current impulse at the sending end.
In this paper, Marti’s model is used for evaluating the effect The injected current has a peak value of 12 kA, a front time of
of ground as a dispersive medium in the calculation of 0.3 µs (measured as the time from 0.3Ip to 0.9Ip, where Ip is the
electromagnetic transients on overhead transmission lines. current peak value), and a maximum steepness of 40 kA/µs
However, the J. Marti setup available in the LCC routine of [30, 31]. This current waveform is representative of
ATPDraw considers the expressions of Carson and Deri et al. subsequent strokes of negative downward lightning measured
[20, 21] to calculate the line parameters. As in their original at Mount San Salvatore, Switzerland [32]. Three different line
form such expressions do not explicitly consider the parameter lengths were considered, namely 600, 1800 and 3600 m. In all
ε, as discussed in Section II, and neglect the variation of σ cases, the line was grounded at both ends by means of a
with frequency, an alternative implementation of Marti’s matching resistance of 497.6 Ω. Three different conditions
model was necessary to evaluate the effect of ground as a were assumed to represent the ground parameters. One, whose
corresponding .pch file is listed in Appendix A, assumed the sending end of the other conductor was grounded and the
frequency-dependent parameters according to (12) and (13) for receiving ends of both conductors were left open. This line
σ0=0.0001 S/m. The other two cases assumed constant ground topology is typical of rural lines used in Brazil and serves to
parameters with σ0=0.0001 S/m calculated either with illustrate the influence of different assumptions regarding the
Carson’s expressions or with Sunde’s expressions for εr=40. ground parameters on switching overvoltages. In the
calculations, a 5-mm conductor radius was assumed together
3 with σ0=0.0001 S/m. The same conditions of the previous
Voltage (MV)

2.5
2
section were assumed for the calculation of the ground-return
1.5 impedance, namely frequency-dependent ground parameters
1
0.5
according to (12) and (13) for σ0=0.0001 S/m, or constant
(a)
0 ground parameters with σ0=0.0001 S/m assuming either
0 5 10 15 20 25 30
Carson’s expressions or Sunde’s expressions with εr=40. The
µs)
Time (µ
.pch file obtained for the frequency dependent case
3
(determined considering a maximum frequency of 107 Hz) is
Voltage (MV)

2.5
2 listed in Appendix B.
1.5
1
0.5 300
(b)
0
0 5 10 15 20 25 30 250
µs)
Time (µ µs)
Phase velocity (m/µ

3 200
Voltage (MV)

2.5
2 150
1.5 σ(ω), εr(ω)
1 σ0, εr=40, Sunde
0.5 σ0, Carson 100
(c)
0 σ(ω), εr(ω)
0 5 10 15 20 25 30 50 σ0, εr=40, Sunde
µs)
Time (µ σ0, Carson
Fig. 5 – Voltages at the receiving end of a single-phase line with length of 0
1 10 100 1000 10000
(a) 600 m, (b) 1800 m or (c) 3600 m assuming the injection of a lightning
Frequency (kHz)
current at the sending end and considering σ0=0.0001 S/m. Both line ends
were connected to a matching resistance of 497.6 Ω. Fig. 6 – Phase velocity associated with a single-phase overhead line with
height of 10 m and radius of 5 mm for σ0=0.0001 S/m and different ground
models.
Taking as reference the voltages calculated assuming
frequency-dependent ground parameters, the results shown in 1.6
Fig. 5 indicate peak values up to 15% higher if a constant
Voltage (V)

1.2
conductivity is assumed in Carson’s expressions. A noticeable 0.8 σ(ω), εr(ω)
difference is also observed in the propagation speed, which 0.4 σ0, εr=40, Sunde
seems to be slower if Carson’s expressions are considered. σ0, Carson (a)
0
This observation is consistent with Fig. 6, which shows that 0 10 20 30 40 50
the phase velocity associated with the use of Carson’s formula µs)
Time (µ

approaches the speed of light slower than the remaining curves 0.8
Voltage (V)

at high frequencies. Interestingly, a very good agreement is 0.4

observed between the voltage waveforms calculated assuming 0
ground as a dispersive medium and those obtained for constant -0.4
(b)
σ with εr=40. Although not shown, a similar agreement was -0.8
observed for ground conductivities above 0.0001 S/m. In any 0 10 20 30 40 50
µs)
Time (µ
case, for σ=0.001 S/m or higher, the influence of frequency-
Fig. 7 – Voltages at the receiving end of an 1800-m long two-phase line
dependent ground parameters and different values of ε was
assuming a 1 pu step voltage source with internal resistance of 150 Ω to
seen to be negligible on the calculated voltage waveforms. energize the topmost conductor while the sending end of the other conductor
was grounded and the receiving end of both conductors were left open: (a)
B. Two-phase distribution line
voltages at the topmost conductor; (b) voltages at the bottom conductor.
Fig. 7 illustrates voltages calculated at the receiving end of Ground conductivity σ0=0.0001 S/m.
an 1800-m long, two-phase distribution line consisting of two
vertically stacked conductors with heights of 10 m and 8.17 m, It is seen in Fig. 7 that the voltage waveforms calculated
respectively. In the calculations, a 1 pu voltage source with with Carson’s expressions present again the largest deviation
internal resistance of 150 Ω supplying a step waveform was from the waveforms calculated considering frequency-
connected to the sending end of the topmost conductor while dependent ground parameters according to (12) and (13). This
is apparent if the voltages induced on the grounded conductor
R2
are analyzed, in which differences of about 20% are observed B
V

in the induced peak values. Also, it is seen that assuming

PCH
V R2 V R2
σ0=0.0001 S/m and εr=40 leads again to voltage waveforms in A C

PCH

PCH

PCH
very good agreement with those calculated assuming
frequency-dependent ground parameters. Although not shown, R1 V

additional tests made for lateral distances up to 4 m between PCH PCH PCH PCH PCH PCH PCH D
both conductors (parameter rij in Fig. 1), which could be R2 R2 R2 R2 R2
considered representative of power distribution lines, were
seen to lead to similar conclusions. Again, as expected, the
differences between the calculated waveforms reduced Fig. 8 – Circuit implemented in ATPDraw to simulate a branched distribution
significantly with increasing ground conductivity. line. Each transmission line block is 150-m long, R1=150 Ω and R2=80 Ω.

C. Two-phase distribution line with branches and 1.6

grounding points

Voltage (V)
1.2
Most of the literature dealing with the influence of 0.8 σ(ω), εr(ω)
frequency-dependent ground parameters on switching and 0.4 σ0, εr=40, Sunde
σ0, Carson (a)
lightning transients on overhead lines disregards the presence 0
of line branches and multiple grounding points. Since this 0 10 20 30 40 50
µs)
Time (µ
condition is typical of power distribution lines, a final case is
presented here in which the branched distribution line 2
1.6
Voltage (V)

illustrated in Fig. 8 is subjected to the switching of a voltage

1.2
source with internal resistance R1=150 Ω supplying a step 0.8
waveform. The configuration of the power distribution line is 0.4
(b)
identical to the one considered in the previous section, 0
0 10 20 30 40 50
consisting of a two-phase line with vertically-stacked µs)
Time (µ
conductors. The neutral conductor is grounded periodically
2
with a grounding resistance R2=80 Ω. This value is adopted by 1.6
Voltage (V)

one of the major power utility companies in Brazil as the 1.2

maximum acceptable value of grounding resistance in their 0.8
distribution lines. Although the accurate modeling of 0.4
(c)
0
distribution transformers would require a detailed pole-residue 0 10 20 30 40 50
representation as well as the presence of surge arresters, the µs)
Time (µ
points A, B, C, and D in Fig. 8, which are left open, could be 2
interpreted as the primary of distribution transformers installed 1.6
Voltage (V)

to supply electrical energy to costumers in rural areas. Each 1.2

transmission line block in Fig. 8 has a length of 150 m and is 0.8
0.4
associated with a .pch file in ATP. As before, three different (d)
0
conditions were assumed in the calculation of the ground- 0 10 20 30 40 50
return impedance, namely frequency-dependent ground µs)
Time (µ

parameters according to (12) and (13) for σ0=0.0001 S/m, or Fig. 9 – Voltages at points (a) A, (b) B, (c) C, and (d) D of the power
constant ground parameters with σ0=0.0001 S/m assuming distribution line of Fig. 8.
either Carson’s expressions or Sunde’s expressions with εr=40.
V. CONCLUSIONS
The .pch file obtained for the frequency dependent case
(determined considering a maximum frequency of 108 Hz) is In this paper, a modal-domain based transmission line
listed in Appendix C. model was used to calculate lightning and switching transients
Fig. 9 shows the voltages calculated at points A, B, C, and considering frequency-dependent ground parameters. Tests
D for the case illustrated in Fig. 8. It is seen that the presence performed in single- and two-phase power distribution lines
of multiple branches and grounding points makes the results indicate that the consideration of ground as a dispersive
nearly independent on the ground-return impedance model, medium can be of some importance in the study of high-
which is an interesting result that might suggest the suitability frequency phenomena on overhead lines located above a
of simplified models for representing the ground parameters in poorly conducting ground (e.g., σ < 0.001 S/m).
certain types of analysis. Once again the results obtained By taking as reference a specific soil model that includes
assuming constant ground conductivity and εr=40 leads to the variation of the ground conductivity and permittivity with
voltage waveforms nearly coincident with the ones calculated frequency, it is shown that assuming a constant value for σ
assuming frequency-dependent ground parameters. together with a suitable value for ε is able to lead to voltage
waveforms in good agreement with those obtained with the C. PCH file of the 150-m long two-phase line of
frequency-dependent ground model. On the basis of the Section IV-C considering frequency-dependent ground
analysis presented in this paper, the use of εr=40 together with parameters
the low-frequency value of the ground conductivity is -1IN___AOUT__A 2. 0.00 -2 2
5 3.5562312867613576000E+02
recommended. In any case, it must be noted that additional 1.40366924312754850E+05
5.79144374266837620E+03
1.20280579296726370E+04
3.97251078211775480E+03
1.01033996388505820E+04

analyses are necessary to assess to what extent this assumption 2.03879833688935610E+04

6.71623531503889380E+00
9.31246882518187820E+01
2.38546185877424130E+00
2.50169970605745960E+01

5 5.0010361896812359000E-07
holds for different soil models. If that is the case, a simple 6.58093722985256200E+09
3.74773839204539220E+03
2.11127069118262040E+06
4.27736359182266880E+01
1.27876661948336260E+05

modification could be made in popular electromagnetic 6.69525582727721880E+09

2.35363219997971370E+06
2.09728303229752210E+08
9.14216626435442160E+04
2.60552706830312460E+07

-2IN___BOUT__B 2. 0.00 -2 2
transient simulators to accommodate the possibility of 10 6.2523360219713925000E+02
2.54382062355948330E+08 1.09120891224525930E+08 1.09851497441550920E+07
adjusting a suitable value for εr in the calculation of the 1.42639674133274950E+07
3.07688682057636800E+04
2.94692618246592670E+06
7.15373162912701990E+03
3.00291081286946840E+05
8.89594132964325040E+03
6.71087897067659650E+03
ground-return impedance of transmission lines. 9.06949432833862120E+06 1.78709259622575830E+06 5.87765236367034840E+05
2.70236454624685170E+05 5.00521462255167930E+04 5.95591089011073930E+03
Finally, the obtained results suggest that the presence of 6.12149900566291420E+02
4.21975457926969360E+00
4.74398752824588290E+01 1.20395534102758270E+01

5 4.9964038993728571000E-07
multiple branches and multiple grounding points is likely to -9.87211106052957920E+07 7.67339900090667960E+08 1.13095461526605150E+07
1.70281252045315550E+06 2.06398618397906680E+05
reduce the relative importance of the assumed ground- 3.20416100549760630E+09
1.08002289264442010E+07
1.23780467419306330E+09
1.55435466974295470E+06
9.55644856004598740E+07
0.74278792 0.62609170
impedance model in the simulation of transients in power 0.00000000 0.00000000
-0.66952476 0.77974458
distribution lines. In any case, a more definitive conclusion in -0.00000000 0.00000000

this direction requires further studies involving more

representative distribution line topologies, different soil VII. REFERENCES
models, and the investigation of other relevant transient [1] C. M. Portela, “Measurement and modeling of soil electromagnetic
phenomena. behavior,” in Proc. IEEE Int. Sym. Electromagnetic Compatibility,
Seattle, WA, pp. 1004–1009, 1999.
VI. APPENDIX [2] C. Portela, M. C. Tavares, and J. Pissolato, “Accurate Representation of
Soil Behaviour for Transient Studies”, IEE Proceedings on Generation,
A. PCH file of the 600-m long single-phase line of Transmission and Distribution, vol. 150, no. 6, pp. 734-744, Nov. 2003.
[3] C. Portela, J. B. Gertrudes, M. C. Tavares, and J. Pissolato, “Earth
Section IV-A considering frequency-dependent ground Conductivity and Permittivity Data Measurements – Influence in
parameters Transmission Line Transient Performance,” in Proceedings of
-1IN___AOUT__A 2. 0.00 -2 1 IPST’2005 – International Conference on Power Systems Transients,
12 4.9750968011831094000E+02
4.89040461178831810E+03 6.51545858335469800E+03 9.35961378202711420E+03 Montreal, Canada, Jun. 2005.
6.01139710016544090E+03 2.63495405613793850E+04 2.20748533525693900E+05
1.95671607001609800E+06 1.67217292504247610E+07 7.40079955983188750E+07 [4] S. Visacro and R. Alipio, “Frequency dependence of soil parameters:
1.11848745691922260E+08 -5.57306595871959460E+05 -1.33388694775600860E+05
3.23339052380764440E+00 7.88937973253849730E+00 2.53650703680675170E+01 experimental results, predicting formula and influence on the lightning
1.00735375837945360E+02 1.01584353861820910E+03 7.84889301241589780E+03
5.92334022100264560E+04 3.95056738199622370E+05 2.30536640092842420E+06 response of grounding electrodes,” IEEE Trans. Power Delivery, vol.
1.50274006472555360E+07 1.76533466748305010E+05 1.76533466748305010E+05
9 2.0007075820523406000E-06 27, no. 2, pp. 927–935, Apr. 2012.
-2.86006971760625270E+09 1.41665507916519280E+09 7.85117445818453280E+07
6.09349532161832790E+05 1.13390858030871430E+07 7.17538648221284970E+05 [5] R. Alipio and S. Visacro, “Modeling the frequency dependence of
1.42939449550848920E+06 5.73560263846185760E+05 2.55616523355094360E+04
2.04201565772296220E+10 4.72799999161449430E+09 5.64534931296067830E+08 electrical paramaters of soil,” IEEE Trans. Electromagnetic
1.60496735912643490E+08 7.88819246550548520E+07 1.09885607834759540E+07
1.04993484687795010E+07
1.00000000
2.10940636139622100E+06 3.23474779308710310E+05 Compatibility, vol. 56, no. 5, pp. 1163–1171, Oct. 2014.
0.00000000 [6] B. Zhu, W. Sima, T. Yuan, Q. Yang, P. Wu, “The influence of soil
permittivity’s frequency dependent characteristics on impulse transient
B. PCH file of the 1800-m long two-phase line of
resistance of grounding electrode,” in Proceedings of GROUND’2014
Section IV-B considering frequency-dependent ground International Conference on Grounding and Earthing, Manaus, Brazil,
parameters May 2014.
-1IN___AOUT__A 2. 0.00 -2 2
[7] S. Visacro, R. Alipio, M. H. Murta Vale, and C. Pereira, “The response
7 3.5525328338766809000E+02
5.27219799731382870E+05 5.96525869055192160E+04 3.51652379277314460E+03
of grounding electrodes to lightning currents: the effect of frequency-
1.25376214549482480E+04 8.03574766479353180E+03 4.84814218458455980E+03 dependent soil resistivity and permittivity,” IEEE Trans.
3.54456926582515870E+03
2.65592602606585250E+05 1.15751995068919150E+04 2.26931144121281530E+02
6.30093547280099240E+01 1.79877012447671450E+01 5.59894863868208330E+00 Electromagnetic Compatibility, vol. 53, no. 2, pp. 401–406, May 2011.
2.30847328455031290E+00
8 6.0016138667939422000E-06 [8] R. Alipio and S. Visacro, “Frequency dependence of soil parameters:
4.89823645263483940E+08 2.38948878808848090E+06 3.42913213363709450E+05
1.88987181593077660E+04 2.09539533277945470E+03 8.75141654430658780E+01 effect on the lightning response of grounding electrodes,” IEEE Trans.
9.36892396008977780E-01 1.19359700571623990E-02
5.43630779248099450E+08 4.72519781140828950E+07 1.13606180526303220E+07 Electromagnetic Compatibility, vol. 55, no. 1, pp. 132–139, Feb. 2013.
1.96248871920381440E+06 3.67746433290539780E+05 4.28077013056875150E+04
2.71205152918589370E+03 3.24632781453123090E+01 [9] R. Alipio and S. Visacro, “Impulse efficiency of grounding electrodes:
-2IN___BOUT__B 2. 0.00 -2 2
16 6.2404953291747472000E+02 effect of frequency-dependent soil parameters,” IEEE Trans. Power
5.34730260097013710E+03 6.07769973441367440E+03 7.05732368433331790E+03
4.20194652521304400E+03 1.28549711811596330E+04 7.25362251106686310E+04 Delivery, vol. 29, no. 2, pp. 716–723, Apr. 2014.
4.15570328443056900E+05 2.91022415981091610E+06 1.24412995455677660E+07
-1.62692694396426750E+07 7.06756743068040760E+07 1.59104622771032040E+08 [10] D. Cavka, N. Mora, and F. Rachidi, “A comparison of frequency-
1.60729530707831050E+07 2.33115750180559250E+08 -3.06987211672537410E+04
-1.54073656522314180E+04 dependent soil models: application to the analysis of grounding
3.98246324225728590E+00 8.23968024585554030E+00 2.10222198408786340E+01
6.79305585368173580E+01 3.64288866527566600E+02 1.94587866347706490E+03 systems,” IEEE Trans. Electromagnetic Compatibility, vol. 56, no. 1,
1.05518592569589930E+04 5.68208524525867110E+04 2.43360757661272710E+05
9.00948709131590210E+05
4.70549229097143280E+06
9.62284225664496540E+05
2.13077614834895580E+07
3.80751224521444600E+06
4.29050611052774470E+04
pp. 177–187, Feb. 2014.
4.29050611052774470E+04
10 5.9100550662727622000E-06
[11] J. Montaña, J. Herrera, J. Rios, and J. Silva, “High frequency behavior
5.42029485526780950E+02
1.99721686964451250E+05
1.49590966931294950E+04 -4.79351422571950590E+02
8.52060722736741300E+05 3.72060725713565840E+04
of grounding systems considering the frequency dependence of soil
-5.16093538375201050E+05 -2.62042446682419300E+06
-1.91223367449165350E+05
1.03994367943618110E+06 parameters,” in Proceedings of GROUND’2014 International
2.65076458732672040E+04 1.52927805642578750E+05 4.08697334562157570E+05
5.30208391726890000E+05 1.33822761907975820E+06 3.28894697890176720E+06 Conference on Grounding and Earthing, Manaus, Brazil, May 2014.
5.35479326034044850E+06 3.65129317464425040E+07 2.75746545196476800E+07
2.75746545196476800E+07 [12] A. C. S. Lima and C. Portela, “Inclusion of Frequency-Dependent Soil
0.74172450 0.62822302
0.00000000 0.00000000 Parameters in Transmission-Line Modeling”, IEEE Trans. Power Del.,
-0.67068771 0.77799217
-0.00000000 0.00000000 vol. 22, no. 1, pp. 492-499, Jan. 2007.
[13] J. B. Gertrudes, M. C. Tavares, and C. Portela, “Transient Performance [20] J. R. Carson, “Wave Propagation in Overhead Wires with Ground
Analysis on Overhead Transmission Line Considering the Frequency Return”, Bell Systems Technical Journal, vol. 5, pp. 539–556, 1926.
Dependent Soil Representation,” in Proceedings of IPST’2011 – [21] A. Deri, G. Tevan, A. Semlyen, and A. Castanheira, “The Complex
International Conference on Power Systems Transients, Delft, the Groud Return Plane – A Simplified Model for Homogeneous and Multi-
Netherlands, Jun. 14-17, 2011. Layer Earth Return”, IEEE Trans.Power. App. Sys., vol. PAS-100, no.
[14] F. H. Silveira, S. Visacro, R. Alipio, and A. De Conti, “Lightning- 8, pp. 3686-3095, Aug. 1981.
induced voltages over lossy ground: the effect of frequency dependence [22] J. R. Marti, “Accurate Modelling of Frequency-Dependent Transmission
of electrical parameters of soil,” IEEE Trans. Electromagnetic Lines in Eletromagnetic Transient Simulations”, IEEE Trans. on Power
Compatibility, vol. 56, no. 1, pp. 1129–1136, 2014. Apparatus and Systems, vol. PAS-101, no. 1, January 1982.
[15] M. Akbari, K. Sheshyekani, and M. R. Alemi, “The effect of frequency [23] Gustavsen, B., Semlyen, A. “Rational Approximation of Frequency
dependence of soil electrical parameters on the lightning performance of Domain Responses By Vector Fitting”, IEEE Trans. on Power Delivery,
grounding systems,” IEEE Trans. Electromagnetic Compatibility, vol. vol. 14, no. 3, July 1999.
55, no. 4, pp. 739–746, Aug. 2013. [24] E. D. Sunde, “Earth Conduction Effects in Transmission Systems”,
[16] M. Akbari, K. Sheshyekani, A. Pirayesh, F. Rachidi, M. Paolone, A. Dover Publications, New York, 1968.
Borghetti, and C. A. Nucci, “Evaluation of lightning electromagnetic [25] F. Rachidi, C. A. Nucci, and M. Ianoz, “Transient analysis of
fields and their induced voltages on overhead lines considering the multiconductor lines above a lossy ground,” IEEE Trans. Power Del.,
frequency dependence of soil electrical parameters,” IEEE Trans. vol. 17, no. 1, pp. 294-302, Jan. 1999.
Electromagnetic Compatibility, vol. 55, no. 5, pp. 1210–1219, Dec. [26] A. Ametani, Y. Miyamoto, Y. Baba, and N. Nagaoka, “Wave
2013. propagation on an overhead multiconductor in a high-frequency
[17] K. Sheshyekani, and M. Akbari, “Evaluation of lightning-induced region”, IEEE Trans. Electromagn. Compat., in press.
voltages on multi-conductor overhead lines located above a lossy [27] A. Ametani, “Stratified earth effects on wave propagation – frequency
dispersive soil,” IEEE Trans. Power Delivery, vol. 29, no. 2, pp. 683– dependent parameters,” IEEE Trans. Power App. Syst., vol. 93, no. 5,
690, Apr. 2014. pp. 1233-1239, Sep. 1974.
[18] R. A. R. Moura, M. A. O. Schroeder, P. H. L. Menezes, L. C. [28] W. H. Wise, “Propagation of high-frequency currents in ground return
Nascimento, A. T. Lobato, “Influence of the soil and frequency effects circuits,” in Proc. I. R. E., vol. 22, no. 4, pp. 522-527, 1934.
to evaluate atmospheric overvoltages in overhead transmission line – [29] F. Rachidi, C. A. Nucci, M. Ianoz, C. Mazzetti, “Influence of a Lossy
Ground on Lightning-Induced Voltages on Overhead Lines”, IEEE
part I: the influence of the soil in the transmission lines parameters,” in
Trans. Electromagn. Compat., Vol. 38, No. 3, pp.250-264, Aug. 1996.
Proceedings of XV International Conference on Atmospheric Electricity,
[30] C. A. Nucci, F. Rachidi, M. V. Ianoz, and C. Mazzetti, “Lightning-
Norman, USA, Jun. 15-20, 2014.
induced voltages on overhead lines”, IEEE Trans. Electromagn.
[19] R. A. R. Moura, M. A. O. Schroeder, P. H. L. Menezes, L. C.
Compat., Vol. 35, no. 1 , pp. 75-86, 1993.
Nascimento, A. T. Lobato, “Influence of the soil and frequency effects [31] A. De Conti and S. Visacro, “Analytical representation of single- and
to evaluate atmospheric overvoltages in overhead transmission line – double-peaked lightning current waveforms”, IEEE Trans. Electromagn.
part II: the influence of the soil in atmospheric overvoltages,” in Compat., Vol. 49, no. 2 , pp. 448-451, 2007.
Proceedings of XV International Conference on Atmospheric Electricity, [32] R. B. Anderson and A. J. Eriksson, “Lightning parameters for
Norman, USA, Jun. 15-20, 2014. engineering application,” Electra, vol. 69, pp. 65–102, 1980.