Transmission Line Modeling Techniques
Transmission Line Modeling Techniques
Line Modeling
Transmission
1. OBJECTIVE
This chapter aims to present the modeling of the lines of
transmission to be used in steady state studies (flow of car-
and short circuit in electrical systems.
The main aspect to be considered in the modeling of a line of
transmission is the difficulty of representing it by an equivalent circuit.
concentrated parameters. A transmission line consists of a
set of distributed parameters, inductance, resistance, capacitance and
dispersion conductance. According to circuit theory, for the lines of
transmissions to be represented by equivalent circuits with parameters
concentrated, they must have dimensions much smaller than ¼ of the
wavelength of the tensions and currents that circulate through it in steady state
permanent. In the particular case of power electrical systems in steady state
permanent sinusoidal at a frequency of 60 Hz, the value corresponding to ¼ of
wavelength is 1250 km, therefore transmission lines of some
Dozens of kilometers or more deserve special attention.
2. INTRODUCTION
Transmission lines are the 'arteries' of electrical systems, they
we transport all the electric energy produced at the generation plants to the
large consumption points. The transmission of energy in large blocks
they can be made through air lines, underground cables, and isolated lines
VI - 1
Chapter 7 - Transmission Line Modeling
P
Vnom@5,5. 0,62.L+ [1]
100
where:
L - length of the line in km
P - power to be transported in kW
VII - 2
Chapter 7 - Modeling of Transmission Lines
This value varies with temperature and with electrical frequency. The vari-
The action with temperature is mathematically expressed by the following equation:
where:
RT1electrical resistance of the conductor at temperature T1
RT2– electrical resistance of the conductor at temperature T2
atcoefficient of resistance increase with temperature
T2conductor temperature where the resistance is desired
T1– temperature of the conductor where the value of the resistance is known
tendency.
VII - 3
Chapter 7 - Transmission Line Modeling
the temperature of 20oC as a reference is worth 0.00385 (1/0C), for the cables
for aluminum, a value of 0.00403 can be adopted (1/0C).
The variation of resistance with frequency is related to the effect
skin or film. This effect causes a cylindrical conductor to be traversed
I slide longitudinally by an alternating current, the current density in the
its interior is smaller along its longitudinal axis and maximum along its su-
perfection.
The skin effect can be easily understood by imagining
the conductor composed of an infinite number of parallel longitudinal fibers
between itself and the longitudinal axis, each representing an infinite conductor
similar. Assuming two cross sections, at a certain distance from each other, the
voltage drop in any of the fibers must be the same. In alternating current-
nothing, in each wire there is not only an ohmic voltage drop, but also
an induced voltage or electromotive force by the alternating magnetic flux. The emf along with the
the surface of the conductor will be smaller than that induced in a fiber along with
a fiber closer to the conductor's axis, since the outer fiber is intertwined
for a magnetic flux smaller than that which encircles the most inter-
So that the voltage drops are equal in the smaller fibers.
inductive reactance that in those with higher inductive reactance, it is necessary that
if the currents in the first are greater than in the second, then the density-
Current density will be greater at the periphery of the conductors.
r=rcc+ra+rad [4]
where:
rcc(ohm/km) - resistance that the conductor presents to the flow
of direct current;
ra (ohms/km) - apparent resistance caused by the existence of
tendency of magnetic flows inside the con-
ducts;
rad (ohms/km) - additional apparent resistance.
VII - 4
Chapter 7 - Transmission Line Modeling
VII - 5
Chapter 7 - Modeling of Transmission Lines
VII - 6
Chapter 7 - Transmission Line Modeling
The conductors are cylindrical and their current densities are uniform.
The presence of the land does not affect the magnetic field, therefore it does not affect the
inductance formulas.
The presence of the land affects the electric field, therefore it affects the formulas of the
capacitance.
exterior of the conductor, in the space that surrounds it. The total magnetic flux
it will be equal to the sum of the internal and external flows.
VII - 7
Chapter 7 - Transmission Line Modeling
ò[Link] = i T
7
where :
H – magnetic field vector
ITtotal current involved in the closed loop
B=mR.m0.H [9]
where:
H - field intensity vector;
B - magnetic flux density vector
-7
0 magnetic permeability of vacuum (4π × 10 );
m-
mR - relative magnetic permeability of the medium
VII - 8
Chapter 7 - Transmission Line Modeling
f = ò[Link] 10
S
y=D y=D
m0.i m0.I D
ò
f = [Link]=
y=R
ò2.p.y .dy=
y=R 2.p
.Ln
R
[11]
By substituting equation [11] into equation [6], we can then obtain the in-
external inductance of this conductor LE
mR.m0 æD ö
THE
= .Ln ç ÷ 12
p èR ø
VII - 9
Chapter 7 - Transmission Line Modeling
x inside the conductor, the current inside is (x/R)2and using Am's law
father applied to this circle we obtain:
mR.m0.x.i
B = mR.m0.H= [13]
2.p.R2
The volume element shown in the figure has a volume of [Link] and the
the total energy stored in the magnetic field within the conductor is:
2
mR.m0 x =Ræ x.i ö mR.m0 2
wmf= .ò ç ÷ .[Link]= .i 14
2 x =0è 2.p.R2 ø 8.p
1
wmf=.LI.i2 [15]
2
m.m
R 0
LI = [16]
8.p
this is:
1 æ ö
mR.m0 mR.m0 æD ö mR.m0 é 4 æD öù mR.m0 ç D ÷
L=
T + .lnç ÷ = . êline +ln ç ÷ú = .lnç 1 ÷ [17]
4.p p èR ø p ë è R øû p ç - ÷
è R.e 4
ø
Analyzing equation [17], we find that by substituting the radius of the con-
by a ray equivalent this equation becomes identical to that of inductance
external to the conductor. This equivalent radius can be interpreted as sin-
of the radius of a fictional, theoretical conductor, which does not possess flux in-
VII - 10
Chapter 7 - Transmission Line Modeling
-1
r'=r.e40.7788 [18]
The term r' is referred to as the geometric mean radius (RMG) and is defined
as if being the radius of a fictitious solid cylindrical conductor that produces
a magnetic field equal to the magnetic field of the real conductor.
VII - 11
Chapter 7 - Transmission Line Modeling
q q
D= = [19]
s2px
where:
2px - area per linear meter of the conductor.
q - number of lines that cross this surface.
D q
E= = [20]
x 2xp x
x = x 0xr 21
where:
x is the electrical permittivity of air.
D2
q D
ò
v 12=Edx=
D1
ln 2
2pxD1
[22]
VII - 12
Chapter 7 - Transmission Line Modeling
D
q 1 q D
v ab = .ò .=dx .ln 23
2pxxR 2px R
And the capacitance of these two conductors is calculated from the charge
per unit of potential, that is:
q px
C ab= = 24
U ab D
ln
R
The dispersion conductance (g) is the line parameter that models the
losses due to dispersion. These losses include losses due to the
corona effect generated by the ionization of the air around the conductor and the losses
dielectrics in insulators. The losses due to the corona effect are distributed
along the line, and the dielectric losses in the insulators are concentrated in the
the same, however, as the distances between support structures are relatively
small is usually considered to be uniformly distributed as well.
From the theory of electrical materials, it is a known fact that no di-
Electric is perfect, that is, all of it has leakage currents. Thus, the insulators
the transmission lines, at the industrial frequency itself, are subject to
these leakage currents that cause energy losses. These losses vary
I am with the quality of the insulator material, with the surface conditions (
pollution, salinity), insulator geometry, frequency of the applied voltage,
weather conditions ...
The term corona is used to describe discharge phenomena.
electric events that occur internally or externally to equipment or with
electric components. The corona can be defined as a localized discharge,
resulting from the transient ionization of gases in an insulation system
when the applied voltage exceeds the critical value.
VII - 13
Chapter 7 - Transmission Line Modeling
There is a distinction between the term corona and partial discharges. The term
corona is associated with the visual phenomenon, partial discharges are phenomena
not visible, as it is internal to the material or component.
The corona effect can produce light, audible noise, and ozone, and can be
detected by our conventional sensory capabilities of vision, audio-
taste and smell.
The selection of drivers is one of the most important decisions to be made.
rem taken by the designer of the transmission lines. In medium voltage lines and
high voltages, the choice of conductor sections is usually related to
I need an economic equation involving losses due to Joule effect.
transport capacity and the necessary investments. In the lines in ten-
of EHV and UHV, the control of the manifestations of the Corona effect can be
the dominant element to guide this choice.
The corona effect appears on the surface of the conductors of an overhead line.
transmission area when the value of the potential gradient present there ex-
yield the value of the critical disruptive gradient of the air. In the non-uniform field in
around a conductor, the divergence of the field exerts additional influence, and
any contaminating particle, such as dust, for example, transforms
in a punctual discharge source.
Electrical discharges in gases are usually initiated by a field
electric that accelerates free electrons present there. When these electrons acquire
enough energy from the electric field, they can produce new electrons by
collision with other atoms. It is the process of ionization by impact. During the
acceleration in the electric field, each free electron collides with oxygen atoms
nitrogen and other gases present, losing part of it in this collision,
its kinetic energy. Occasionally an electron can hit an atom with
sufficient force, in order to excite it. Under these conditions, the affected atom reaches
to an initial state, releasing the excess energy in the form of heat, light,
acoustic energy and electromagnetic radiation. An electron can also
collide with a positive ion, converting it into a neutral atom.
VII - 14
Chapter 7 - Transmission Line Modeling
All the energy released must come from the electric field of the line,
so much, from the feeding system, for which it represents energy loss, due to
subsequent, loss. The losses are related to the geometry of the conduit.
res, operating tensions, potential gradients on the surfaces of the conductors
stores and with the local weather conditions.
Individual corona discharges cause pulses of voltage and color.
short-term rates that propagate along the lines, resulting in
electromagnetic fields in their vicinity. These discharges occur under
cover both semicircles of the applied voltage, however those that occur
during the positive semicircles is when noise is emitted that can interfere with
radio reception in the frequency bands of amplitude transmissions
sweets.
As described at the beginning of this chapter, due to the extent of the li-
transmission lines, model a line using its parameters concen-
trados can lead to significant errors in analyses involving the com-
her permanent regime behavior. From the theory of electric circuits, it is co-
It is known that in order to model a component by a circuit to
concentrated parameters need their dimensions to be very me-
no more than 1/4 of its wavelength, that is:
l c 300,000
= = =1250 Km
4 4. f 4.60
where:
c - speed of light propagation;
f - frequency in Hz
VII - 15
Chapter 7 - Transmission Line Modeling
i(x D
+ x, t)
i(x,t)
r. Dx l. Dx
v(x D
+ x, t)
v(x,t) [Link]
[Link]
Applying Kirchhoff's laws of voltage and current in the section of the line
of length Dx shown in the figure, it is shown that:
¶v
i(x,t) = i(x + Dx, t) + v(x + Dx, t) . g . Dx + c . Dx. [25]
¶t
¶i
v(x,t)=v(x+Dx,t)+i(x,t).[Link]+[Link]. [26]
¶t
v(x+Dx,t)-v(x,t) i
= -r.i-l. [27]
Dx ¶t
i(x+Dx,t)-i(x,t) ¶v
=-g.v-c. [28]
Dx ¶t
VII - 16
Chapter 7 - Transmission Line Modeling
¶i ¶v
= -g.v-c. [30]
¶x ¶t
¶ 2v ¶i ¶ 2i
= -r.-l. [31]
2
¶ x ¶x ¶x.¶t
¶ 2i ¶v ¶ 2v
= -g.-c. [32]
2
¶ x ¶x ¶x.¶t
¶ 2v ¶i ¶ 2I
= -r.-l. [33]
¶t.¶x ¶t .¶2 t
¶ 2i ¶v ¶ 2v
-g.-c. [34]
¶x.¶t ¶t ¶ 2t
Finally, substituting equations [29], [30], [33], and [34] into the equations
From [31] and [32], it is obtained:
¶ 2v é ¶v ù é ¶v ¶v2 ù
=-r.-gv-c
êë ú -l.-g.-c
ê ú 35
¶ x
2
¶t û ë ¶t ¶ 2t û
¶ 2i é ¶i ù é ¶i ¶i2 ù
= - -[Link]
êë l- ú - c.r.ê- -l 2 ú [36]
¶ x
2
¶t û ë ¶t ¶ tû
¶ 2v
(cr+lg ) ¶v +lc ¶ 2v
2
=r.g.v+ [37]
¶ 2x ¶t ¶ t
¶ 2I
=r.g.i+ (cr+lg ) i +lc
¶ 2i
38
¶ 2x ¶t ¶ 2t
VII - 17
Chapter 7 - Transmission Line Modeling
Applying the phasor method to the differential equations [37] and [38], obtain
mos:
¶ 2V
=z.y.V [39]
¶ 2x
¶ 2I
=z.y.I 40
¶ 2x
A gx B - gx
I= - e + e 42
Zc Zc
Z
Zc = [44]
y
VII - 18
Chapter 7 - Transmission Line Modeling
Aunt of phaseb.
Analyzing equations [41] and [42], it can be seen that the voltage phasor and
the current phasor consists both of the sum of its modulated sinusoids
measures that propagate along the axis "x" through the exponentials:
e ±g =e±ax e ± jbx .
Each sine wave "suffers":
x
a) a damping caused by the term eawhich depending on the signal of
exponent is positive or negative, this can increase the amplitude of the se-
nodes as the distance x varies or reduces the amplitudes of the se-
nodes, if it is positive or negative respectively.
b) At the same time, there is a progress of ebinxthe wave phase to which it is applied
I fell.
Thus we conclude that 'g' commands the way in which tensions and cor-
rents propagate along the line, hence the designation of pro- function
propagation constant or propagation constant.
From the known values of voltage and current at the beginning and
At the end of the transmission line, expressions for the constants A and can be obtained.
VII - 19
Chapter 7 - Transmission Line Modeling
(V2-ZcI 2 )
A= e - g1 [46]
2
æV-Z I c 2ö g1
B= ç 2 ÷e [47]
è 2 ø
Using equations [41] and [42] to obtain the values of the phasor ten-
are and the current phasor at x=0, substituting the values of A and B into the equations
from equations [46] and [47], we obtain:
y2 -y2 y1 -y1
æe+e ö æe-e ö
ç
V1=V2 ç ÷
÷ +ZcI 2 çç ÷
÷
2 2 [48]
è ø è ø
y2 -y2 y1 -y1
æe+e ö V2 æe-e ö
I 1I=2 çç ÷ +
÷ Z ç ç ÷
÷ [49]
è 2 ø c è 2 ø
e y1+e- y1 e y1-e- y1
how coshg1= e senhg1= finally we obtain a
2 2
general equation that allows obtaining the voltage and current at the beginning of a line
of transmission in steady state when the voltage and current at the end
the transmission line is known, then:
V=V
1 2 cosh(g.l0)+Zc I 2senh(g.l0 ) 50
V2
I1I2cosh(g.l0) + sinh(g.l0 ) [51]
Zc
8. SHORT LINES:
They are the transmission lines that we can entirely disregard.
mind the effects of capacitance and dispersion conductance. They are lines
whose performance in steady state can be determined by substituting
VII - 20
Chapter 7 - Transmission Line Modeling
in equations [50] and [51], the hyperbolic terms (hyperbolic cosine and sine
hyperbolic) by the first terms of the power series, without incurring errors
greater than 1%. So, doing: sin(h.g.l0)@g.l0ecosh(g.l0)@1, is located:
V=V
1 2+Zc I 2 g.l0 [52]
I1@I2 [53]
z
characteristic impedance and propagation constant, Zc.g= .zy z=, of
y
so that:
V=V
1 2+z.l0.I2=V2+ZLT.I2 [54]
based on the resistance of the cable and the reactance of the line in ohms per kilometer
R X
LT LT
V V
1 2
From a practical point of view, short lines are those that fit within the se-
following conditions:
They are lines up to 150kV, with maximum lengths of 80km
VII - 21
Chapter 7 - Transmission Line Modeling
They are lines with voltages greater than or equal to 150kV, but less than
9. AVERAGE LINES
They are transmission lines whose performance in steady state
can be determined by substituting into equations [50] and [51], the hyper- terms
hyperbolic (hyperbolic cosine and hyperbolic sine) by the first two terms of
power series, without incurring errors greater than 1%.
They are lines that can be characterized by:
a) have a length of 200km with nominal voltages
above 150kV and below 400kV.
b) They have a maximum length of 100km with nominal voltages
nice above 400kV.
The capacitance is included in the calculation of medium lines and the terms hi-
senh(g.l0)@(g.l0)+
(g.l0 )3 e cosh(g.l0)@1+
(g.l0)2
3! 2!
Substituting them into the general equations of performance analysis in
permanent regime of the transmission line (equations [41] and [42]) is found
se:
æ ZLTYLTö é ZLTYLT ù
V1 = V2ç .+ 1 ÷ + I 2.Z LT. ê1 + [55]
è 2 ø ë 6 úû
æ [Link]
ö é [Link]
ù
I=I1+2ç ÷ +V2.Y1+LT ê [56]
6 úû
1
è 2 ø ë
VII - 22
Chapter 7 - Modeling of Transmission Lines
Zp
Yp Yp
2 2
Figure 7 - Pi Circuit
æ Z Y ö æ Z Y ö
I1I= 12ç + ð ð ÷ + V2 Yð 1ç + ð ð ÷ [58]
è 2 ø è 4 ø
ZLT
YLT YLT
2 2
VII - 23
Chapter 7 - Transmission Line Modeling
ZLT/2 ZLT/2
YLT
In both equivalent circuits ZLTe YLT they are the concentration parameters
trados of the transmission line, so that, even in the medium lines the
the effect of the parameter distribution does not have a considerable influence on the analysis
it is necessary to use the complete equations [50] and [51], in the form ex-
exponential or in hyperbolic form.
For the modeling of these transmission lines in steady state
The equivalent pi and T circuits are usually used, whose parameters
are obtained by comparing the coefficients of the hyperbolic terms of the equations
and with the equations of the pi circuit. From this comparison, it can-
it is obtained that:
Zp=Zc senh(g.l0 )
VII - 24
Chapter 7 - Transmission Line Modeling
finally
senh(g.l0)
ZpZ=LT = Z [Link] [59]
g.l0
Zð Yð
cosh(g.l0)=1+
2
you
Yð cosh(g.l0) -1
=
2 Zð
likeZp=Zc senh(g.l0):
Yð cosh(g.l0)-1 1 æg.l0ö
= = .tanh ç ÷
2 Z [Link]( g.l )0 Z C è 2 ø
Yð YLT
= .FCA 60
2 2
VII - 25
Chapter 7 - Transmission Line Modeling
11. BIBLIOGRAPHY
Weedy, B.M. – Power Electrical Systems;
Stevenson, W. D. – Elements of Power System Analysis;
Fuchs, Rubens D. – Transmission of Electric Energy;
Elgerd, Olle - Introduction to Electric Power System Theory;
Beeman - Industrial Power Systems Handbook
Stagg - Applied Computing to Electric Power Systems;
Neuenswander - Modern Power System;
Celso Camargo – Transmission of Electric Energy
[9] Zaborsky, J and Rittenhouse, J. W. – Electric Power Transmission
[10] EHV Project, Transmission Line Reference Book
VII - 26