EMTP RV Modelling

1
TRAINING DUBROVNIK, CROATIA - APRIL, 27 - 29 2009

SIMULATION & ANALYSIS OF POWER SYSTEM TRANSIENTS WITH
EMTP-RV

Modeling of Transmission Line and Substation for
Insulation Coordination Studies
Prof. Ivo Uglei

Faculty of Electrical Engineering and Computing
University of Zagreb, Croatia
2
OUTLINE OF PRESENTATION
INTRODUCTION
MODELING GUIDELINES
LIGHTNING MODEL
TOWER
INSULATOR
FOOTING RESISTANCE
LINE, CONDUCTORS AND EARTH WIRES
BOUNDARY CONDITIONS
SUBSTATION MODEL
SURGE ARRESTER
EXAMPLE
3
INTRODUCTION
Computer modeling of transmission lines and substation helps
engineers understand how protection systems behave during
disturbances and faults.

Any transient disturbance, such as lightning stroke terminating on a
phase conductor can be analyzed by use of traveling wave.

A lightning stroke to a conductor or the closing of a circuit breaker
produces traveling waves of voltage u(t) and current i(t) that are related
by a surge impedance Z equal to formula that travels along the
conductor at the speed of light c.
) (
) (
t i
t u
Z =
4
I, t
f
Back-flashover
I, t
f
Shielding failure
Induced overvoltage
I, t
f
I t
INTRODUCTION (Lightning overvoltages on HV transmission lines)
5
Definitions of insulation coordination:

Insulation coordination is the selection of the insulation strength.

Insulation coordination is the selection of the dielectric strength of
the equipment in relation to the voltages which can appear on the
system for which equipment is intended and taking into account
the service environment and the characteristics of the available
protective devices (*) .
(*) IEC 60071-1-1993-12: Insulation coordination Part 1: Definitions, principles and rules.
Line insulation coordination; transmission and distribution lines.

Substation insulation coordination; generation, transmission and
distribution substation.
INTRODUCTION
6
MODELLING GUIDELINES
There are various modeling strategies for lightning transient
studies have been presented elsewhere.

The summary of modeling guidelines that had been adapted:

IEC/TR 60071-4 Edition 1.0 (2004-06): Insulation co-ordination -
Part 4: Computational guide to insulation co-ordination and
modeling of electrical networks;

IEEE PES Task Force on Data for Modeling System Transients
of IEEE PES Working Group on Modeling and Analysis of System
Transients Using Digital Simulation: Parameter Determination for
Modeling System Transients, IEEE Transactions on Power
Delivery, Vol. 20, No. 3, July 2005.
7
MODELLING GUIDELINES
CIGRE, Working Group 01 of Study Committee 33: Guide to
Procedures for Estimating the Lightning Performance of
Transmission lines, Paris, October 1991.

IEEE Working Group 15.08.09: Modeling and Analysis of System
Transients Using Digital Programs, 1998.

IEEE Working Group: A Simplified Method for Estimating
Lightning Performance of Transmission Lines, IEEE Transactions
on Power Apparatus and System, Vol. 104, No. 4, April 1985.
8
LIGHTNING MODEL
Lightning stroke is represented as
a current source with magnitudes
between a few kA to over 200 kA.
Peak current magnitude and tail
time are important when
observing energy stresses of SA
(simplest representation is double
ramp).
Current wavefront is an important
parameter with regard to insulator
flashover.
CIGRE model describes well the
concave wavefront of a lightning
current.
t
90
t
30
t
h
I
30
I
50
I
90
I
100
I
p
t (s)
I (kA)
t
h t (s)
I
50
I
100
I
p
I (kA)
t
f
Double ramp shape
t
f
front time
t
h
time to half-value

9
A statistical approach considering the ground flash density at the
location is used for the determination of lightning parameters
such as:
crest value;
front time;
maximum current steepness;
duration.
The probability that a certain peak current will be equal or greater
than a current I can be determined by Andersons distribution:

6 . 2
31
1
1
|
\
|
+
=
I
P
Where:
P(I) = the probability that the peak current in any stroke will exceed I
I = the specified crest current of the stroke in kA.
LIGHTNING MODEL
10
Steepness can be determined as:
Coefficients o |
First stroke
S
30
3.2 0.25
S
m
3.9 0.55
Subsequent stroke
S
30
6.9 0.42
S
m
3.8 0.93

|
o I S =
LIGHTNING MODEL
11
LIGHTNING MODEL CIGRE model in EMTP RV
The model parameters are:

t
start
- start time, if t < t
start

the source is an open-
circuit;

I
max
- maximum current;

t
f
- from time;

S
m
- maximum steepness;

t
h
- time to half value;

t
stop
- stop time, if t > t
stop

the source is an open-
circuit. The stop time
must be greater than
the start time.
12
Extensive research has been performed on the response of
vertical towers to lightning strokes, and research is still
continuing.

The response of a tower is an electromagnetic problem, although
its study often relies on the circuit approach and models that are
simple to apply in transient simulations:
simple distributed line model,
multistory tower model.

Simple distributed line model provides a constant value of surge
impedance and the constant velocity of travel along the tower.

Different formulas are applied for various tower types.
TOWER
13
1
r
2
r
3
r
1
h
2
h
R
H
O half-angle of cone,
H tower height [m],
R tower base radius [m].
TOWER - Simple Distributed Line Model
The tower surge impedance depends on the direction of wave
propagation and the shape of a lightning current.
The average surge impedance recommended by IEEE and CIGRE:
|
|
\
|
|
|
\
|
|
\
|
=
|
|
\
|
|
\
|
O
=
H
R
Z
1
tan
2
1
cot ln 60
2
cot ln 60
H R
R
H
Z <<
(
\
|
= 1 ln 60
H
h r H r h r
R
) (
1 3 2 2 1
+ +
=
Radius R is calculated by dividing the tower into
upper and lower truncated cones:
An approximation of surge impedance equation is determined by
equivalently replacing the tower with a cylinder.
14
TOWER - Multistory Tower Model

Multistory tower model is developed
for representing towers of UHV
transmission lines (*). Its parameters
were revised according to the results
of experimental studies (**).

The model is composed of four
sections representing the tower
sections between cross-arms.

Each section consists of a lossless
line in series with a parallel R-L circuit,
included for attenuation of the
traveling waves.
1
l
2
l
3
l
4
l
1
R
1
L
1 t
Z
2 t
Z
4
R
4
L
2
R
3
R
f
R
2
L
3
L
1 t
Z
1 t
Z
(*) M. Ishii, T. Kawamura, T. Kouno, E. Ohsaki, K. Shiokawa, K. Murotani, and T. Higuchi, Multistory transmission tower model for
lightning surge analysis, IEEE Trans. Power Delivery, vol. 6, July 1991, pp. 13271335
(**) Yamada, T.; Mochizuki, A.; Sawada, J.; Zaima, E.; Kawamura, T.; Ametani, A.; Ishii, M.; Kato, S.; Experimental evaluation of a
UHV tower model for lightning surge analysis IEEE Trans. Power Delivery, Vol. 10, No. 1, Jan. 1995 pp 393 402
15
INSULATOR
The critical flashover voltage (CFO) is the impulse
voltage level at which the probability of flashover
of the insulator is 50%.

Flashover should not happen when the line
arrester is installed in parallel with the insulator
since the residual voltages developed across
surge arrester are much lower than the dielectric
strength of insulators, even for the highest stroke
currents.
Flashover voltage of line insulators should be randomly varied
according to the statistical distribution laws with the appropriate
standard deviation.
16
INSULATOR Flashover - Leader Propagation Model
v leader velocity (m/s)
U(t) - voltage across the gap (kV)
g - gap length (m)
l - leader length (m)
E
10
- critical leader inception gradient (kV/m)
k
l
- leader coefficient (m
2
V
-2
s
-1
)

The leader propagation stops if the gradient in the unbridged part
of the gap falls below E
10
.
(
= =
10
) (
) ( E
l g
t U
t U k
dt
dl
v
l
l
g
U(t)
Leader
The leader progression model is used to represent line
insulation flashovers:
17
INSULATOR Flashover - Volt-time Characteristic
The flashover voltage characteristic of the insulator string is
time dependent.

The experimental volt-time characteristic is only adequate for
relating the peak of the standard impulse voltage to the time of
flashover.

An open switch connected to insulator string terminals can control
the flashover voltage characteristic.
Time to breakdown (s)
18
INSULATOR Flashover - Area Criterion Model
The method allows the applied nonstandard waveform to be taken
into account.
It involves determining the instant of breakdown using a formula:
D dt V t V
k
t
T
gap
>
}
) ) ( (
0
0
(*) IEC 60071-4: Insulation co-ordination Part 4: Computational guide to insulation co-ordination and modeling of
electrical networks, 2004.
V
gap
(t) - voltage applied at the time t, to the terminals of the air gap,
V
0
- minimum voltage to be exceeded before any breakdown process
can start or continue,
T
0
- time from which V
gap
(t) > V
0
,
k, V
0
, D - constants corresponding to an air gap configuration and
overvoltage polarity (*).
Flashover occurs when the integral becomes greater or equal to D.
The parameters V
0
, k and D are determined by using the voltage-
time curve.
19
TOWER - Example

( ) H R
R
H
Z <<
)
`
\
|
= 1 ln 60
Tower surge impedances are calculated using equation:
Each tower is divided in four
parts. First part is from tower top
to upper arm, second one from
upper arm to middle arm, third
part from middle arm to lower arm
and the last part from lower arm to
ground. On this way it is possible
to calculate transient voltages of
tower arms.
20
TOWER

EMTP RV Modelling

Uploaded by

Copyright:

Available Formats

EMTP RV Modelling

Uploaded by

Document Information

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

EMTP RV Modelling

Uploaded by

Copyright:

Available Formats

How does computer modeling help engineers understand system behavior during disturbances?

What are the guidelines for modeling transmission lines and substations for insulation coordination studies?

1

TRAINING DUBROVNIK, CROATIA - APRIL, 27 - 29 2009

You might also like