MODELADO DE LÍNEAS MONOFÁSICAS APLICANDO EL

MÉTODO DE BERGERON
MODELING OF SINGLE-PHASE LINES APPLYING THE
BERGERON METHOD

Juan C. Escamilla, Filiberto Candia García, Ismael Albino, Ing. Miguel López Castillo
Colegio de Ingeniería Mecánica y Eléctrica, Facultad de Ingeniería, Benemérita Universidad Autónoma de Puebla
Facultad de Ingeniería, Ciudad Universitaria, Bvld. Valsequillo y Av. Sn. Claudio, Col. San Manuel, Puebla,
México Tel. (222) 2295500, ext. 7627
[email protected]

Resumen
En este artículo se presenta un modelo de parámetros de la línea de transmisión monofásica distribuidos
en el dominio de tiempo para el análisis de sobretensiones transitorias a lo largo de la línea de
transmisión. El método implementado es el modelo de Bergeron para el cual se considera la línea sin
pérdidas y con parámetros eléctricos constantes, por lo que el efecto de dispersión no está considerado
en este estudio. El modelo se utiliza en diferentes escenarios de simulación de línea abierta o cerrada a
diferentes instantes de tiempo.
Palabras clave: Línea de transmisión, sobretensiones transitorias, modelo de Bergeron, efecto de
dispersión.

Abstract
This article presents a model of monophasic transmission line parameters distributed in the time
domain for the analysis of transient overvoltages along the transmission line. The implemented method
is the Bergeron model for which the lossless line and constant electric parameters are considered, so the
dispersion effect is not considered in this study. The model is used for different simulation scenarios of
an open or closed line at different periods of time.
Key words: Transmission line, transient overvoltages, Bergeron model, dispersion effect.
 I. INTRODUCCIÓN están expuestas a fallas internas, sino son las que
están expuestas con más frecuencia a fallas
Para poder disponer de la electricidad, es
externas debido a la longitudinalidad que
fundamental contar con un sistema que permita su
presentan.
generación, transmisión y distribución. Dicho
sistema se conoce como Sistema Eléctrico de Determinar la magnitud de sobretensiones es de
Potencia (SEP), el cual está compuesto suma importancia en el análisis de las líneas de
principalmente por generadores, líneas de transmisión, debido a que de ello depende la
transmisión, subestaciones, transformadores y coordinación de aislamiento y el ajuste de
redes de distribución. Cada uno de los elementos protecciones. De esta manera, comprender el
enlistados anteriormente deben ser los adecuados tema representa una base importante en el estudio
para poder tener un servicio eficiente de energía; de transmisión de energía en el área de ingeniería
sin embargo, no suele ser suficiente contar con eléctrica. Enfocar un análisis sugiere conocer de
equipos bien diseñados, ya que pueden fallar en una manera más profunda las características y
cualquier momento, lo cual no necesariamente propiedades de un cuerpo u objeto, así que, para
significa que el equipo fallará por sí solo, sino una línea de transmisión, esto requiere de algo
que éste puede ser afectado por fallas internas más que la determinación de las sobretensiones.
(propias del sistema) o externas (descargas Por consiguiente, para llegar a una solución, se
atmosféricas) [1-6]. Este tipo de eventos generan deben considerar diferentes fenómenos que
transitorios electromagnéticos, los cuales se ocurren cuando la energía se propaga a lo largo
definen como un cambio repentino en las de las líneas de transmisión; es decir, es necesario
condiciones del sistema producido por la conocer todos aquellos fenómenos que permitan
interacción de los campos eléctricos y realizar un diseño oportuno y eficiente de la línea
magnéticos. El período de un transitorio suele ser para asegurar una transmisión de energía sin
muy corto, medido en milisegundos. No obstante, interrupciones por fallas en sus elementos.
estos períodos transitorios son extremadamente
Entre los efectos más comunes que producen
importantes, ya que es en esos momentos que los
alteraciones en la transmisión de la energía,
componentes del circuito están sometidos a
podemos enlistar los siguientes: a) efecto corona,
corrientes o tensiones excesivas [2-5]. Las líneas
de transmisión no sólo
b) efecto aislador, c) efecto piel, d) efecto de II. ALCULO DE
proximidad y e) efecto de pasividad. En la PARÁMETROS
actualidad, existen modelos que permiten analizar
ELÉCTRICOS PARA
las sobretensiones en las líneas de transmisión.
LÍNEAS MONOFÁSICAS
Estos modelos se pueden clasificar como modelos
de parámetros concentrados y parámetros El estudio y modelado de líneas de
distribuidos, donde este último es la transmisión requiere de aquellos parámetros
representación más detallada de las líneas de que definen la caída de tensión y la variación de
transmisión [2-12]. En este trabajo se presenta el corriente. Estos parámetros eléctricos son
método de Bergeron que será aplicado en una resistencia (𝑅) e inductancia (L), los cuales
línea monofásica para obtener las sobretensiones definen la caída de tensión formando la
transitorias. Este método es una técnica en el impedancia en serie (Z): 𝑍 =
dominio de tiempo que permite considerar la 𝑅 + 𝑗𝑤𝐿, donde w es la velocidad angular. Por
propagación de las ondas de voltaje y corriente, otro lado, la conductancia (𝐺) y la susceptancia
así que está clasificado como modelo de (𝐵) en derivación forman la admitancia
parámetros distribuidos. Tal método fue de los transversal de la línea definida de la siguiente
primeros modelos en ser implementados en forma: Y = 𝐺 + 𝑗𝐵. Debido a que se modelará
programas de simulación de tipo EMTP- la línea sin pérdidas, sólo se presenta el cálculo
ATP/Draw y, en la actualidad, se sigue de la inductancia y la capacitancia geométricas.
utilizando, principalmente cuando no se requiere Para el cálculo de la inductancia geométrica, se
mucha precisión en los resultados, ya que una de utiliza el método de las imágenes mostrado en la
las limitantes del modelo es que las pérdidas y el Fig. 1, asumiendo que el retorno por tierra es un
efecto piel no son considerados de forma conductor perfecto [2, 6, 7]. El concepto de
directa. Esta técnica fue inductancia se representa por la ecuación:
desarrollada originalmente por el analista de

sistemas hidráulicos L. Bergeron en 1949 y ha 

L (1)
sido adaptada en el área de ingeniería eléctrica I

para el análisis y los estudios de los transitorios

electromagnéticos en líneas de transmisión [2, 7].
 rcond < 2h, por lo que la inductancia geométrica
+
(Lg) se aproxima como:

0  2h 
(4)
L  2 ln r
g  

 cond 

Para el cálculo de la capacitancia también se hace

uso del método de las imágenes. El concepto de
capacitancia se representa por la ecuación:
-
qc
C (5)
vc
Figura 1. Método de las imágenes para la línea monofásica

La relación del campo eléctrico E y el voltaje se

puede obtener a partir de:
donde Ψ representa el flujo magnético producido a

por la corriente que circula en el conductor e I (6)

Vab   E
b
representa la corriente en el conductor. Para
dl
obtener un valor adecuado y preciso, se deben donde qc representa la carga del conductor y vc
considerar los flujos internos y externos del representa el voltaje. A partir del método de las
conductor bajo estudio. A partir del método de las imágenes de la Fig. 1 se obtiene:
imágenes se obtiene la siguiente ecuación para el q  2h  r  (7)
vc  ln c cond

flujo magnético: 20  rcond 

1
0 I 0 I donde ε0 es la constante de permitividad en el

2hr2
dr  
h
2hr1
2r dr
(2)
 r1 2r
2
vacío. Para una línea de transmisión, se tiene que

De la suma de Ψ1 y Ψ2 se obtiene: rcond < 2h. Combinando (7) y (5), Cg se define por:

 I  2h  r 
(3) 20 (8)
Cg   2h 
 0
ln  cond

2  rcond  ln  
r
 cond 

donde h representa la altura del conductor, rcond es el vacío. Para una línea de transmisión se tiene
el radio del conductor y μ0 es la permeabilidad en que
III. DESCRIPCIÓN DEL
MÉTODO DE BERGERON
Las ecuaciones que describen el comportamiento
de las sobretensiones que se propagan a lo largo
de la línea de transmisión despreciando las
pérdidas se definen por:
i (x, t) v (x, t) (9)
 c
x t Figura 3. Atenuación de la onda viajera
v (x, t) i (x, t)
 L (10) La solución de las ecuaciones diferenciales
x t

En la Fig. 2 se muestra la propagación de las parciales (9) y (10) son (11) y (12),
ondas viajeras a lo largo de la línea de respectivamente [2, 6, 7],
transmisión. El
rebote de las ondas que se muestra en el extremo 
 x
v(x, t)  t     x  (11)
  v
t
final de la línea ocasiona las sobretensiones v  v
v
transitorias. Para determinar y cuantificar tanto su 
i (x, t)    x   x  (12)

Yc v t  
t    v
vv
forma de onda como sus magnitudes, se resuelven     

las ecuaciones (9) y (10) en el dominio de tiempo.

donde 𝑢 es la velocidad de propagación de las
La atenuación de la forma de onda dependerá de
ondas viajeras a través de la línea, t es el tiempo,
dos parámetros eléctricos, las pérdidas de la línea
Yc es la admitancia característica y x es el eje de
de trasmisión y la carga conectada. La Fig. (3)
propagación [2, 6, 7]. A partir de la combinación
muestra la atenuación de la forma de onda
de (11) y (12) se obtiene una igualación que
durante su propagación.
permite relacionar la onda incidente con la onda
reflejada.
 x  x
2v  t   2v  t 
  



(13)
 v  v

v(x, t)  Z c i (x, t)  v(x, t)  Z c i (x, t) (14)

Figura 2. Representación del rebote de una onda viajera donde Zc es la impedancia característica.
Aplicando condiciones de frontera a las
ecuaciones (13) y (14), según la Fig. 4, se obtiene
lo siguiente:
 Figura 5. Circuito dual de Norton para una línea de
transmisión
Figura 4. Condiciones de frontera

v (t  )  Z
i (t  )  v (t)  Z (t) (15)
i
0 c0 L c L
La Fig. 6 muestra la representación de una línea

v (t  )  Z i
(t  )  v0 (t)  Zc i0 (t) (16) radial excitada con una fuente de tensión y con
L c L

una carga conectada en el extremo receptor.

Despejando las corrientes correspondientes a
ambos extremos de la línea de (15) y (16), se
obtiene:

i  (t)  Y v (t)  I 0
(17)
L 0 L Hist

i  (t)  Y v (t)  I L
(18)
0 0 0 Hist
Figura 6. Circuito equivalente de una línea de transmisión
utilizado por el Modelo de Bergeron
donde:

V  (t  )
0
I Hist  0
 i 0 (t  ) (19) Por medio de la Fig. 6 se obtiene el siguiente
Zc
sistema de ecuaciones que relaciona los voltajes
V (t  )

en ambos extremos de la línea y los términos de
L
I Hist  L
 i L (t  ) (20)
Zc
historia correspondientes.

Las ecuaciones (19) y (20) son términos de

historia V0   Yc 0 I IL  (21)

Yf I
  inv  f Hist
0 
que representan voltajes y corrientes desfasados VL   0 YL  Yc  Hsit 
un
tiempo de viaje. Con las ecuaciones (17) y (18) se El método de Begeron es un método gráfico que
forma el circuito dual de Norton para una línea de realiza el cálculo de voltajes y corrientes en cierto
transmisión como el que se muestra en la Fig. 5. instante de tiempo en función de variables
desfasadas un tiempo de viaje con respecto al representación gráfica se muestra en la Fig. 7.
instante de tiempo de referencia. Esta

tener una carga con impedancia resistiva muy

grande. En las Figuras 9 y 10 se observa que la
tensión se duplica y se mantiene prácticamente
constante. Estas oscilaciones no se atenúan
debido a que el método no incluye las pérdidas en
su modelado. Al conectar una carga en el extremo
receptor de la línea, las oscilaciones y
sobretensiones se atenúan debido a que la
resistencia las amortigua. Esto se muestra en las
Figuras. 11 y 12. Por otro lado, cuando se conecta
una carga con un valor muy cercano a la
impedancia característica de la línea,
Figura 7. Posición de las corrientes de historia y ondas de
energía durante el tiempo de observación prácticamente no se presenta sobretensión
transitoria debido a que las reflexiones en el
IV. APLICACIÓN DEL extremo receptor se anulan. Este comportamiento
MÉTODO DE BERGERON se muestra en las Figuras 13 y 14.

El algoritmo fue implementado en el lenguaje de

programación Matlab. Los casos de simulación
que se presentan son los siguientes: a) simulación de línea abierta. La programación
energización de una línea en vacío, b) es equivalente a
energización de una línea con carga y c) línea con
descarga atmosférica. La secuencia del algoritmo
se muestra en la Fig. 8. Los datos de la línea que
se simularon son de fuente función escalón
unitario y una función cosenoidal, donde la
resistencia de la fuente es de 0.001 Ω. La
resistencia de la carga se varió durante las
simulaciones, el radio del conductor fue de 0.02
m, la altura fue de 28 m y la longitud fue de 100
km. Para el primer caso, se realiza la
 Inicio

Figura 8. Diagrama de flujo del método de Bergeron

Ingreso de datos de simulación para la línea de transmisión

Cálculo de parámetros eléctricos de la línea

Cálculo de corrientes de historia

Cálculo de tensiones transitorias

Gráfica de tensión transitoria

Fin
 2
2.5 V - Fuente
V - Fuente
V - ZL = 10 KOhms
V - ZL = 1 GOhm 1.5
2
1

1.5

Voltaje (p.u.)
0.5
Voltaje (p.u.)

1 0

-0.5
0.5

-1
0
-1.5
0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04
-0.5 Tiempo (s)
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
Tiempo (s)
Figura 12. Sobretensión por energización de línea con
Figura 9. Sobretensión transitoria por energización de línea
impedancia de carga de 10 𝐾Ω
abierta con fuente función escalón unitario
1.5
V - Fuente
V - ZL = 1000 Ohms
2 V - ZL = 800 Ohms V - ZL = 600 Ohms
V - Fuente
1.8 V - ZL = 10 KOhms
1
1.6

Voltaje (p.u.)
1.4

1.2
Voltaje (p.u.)

1 0.5

0.8

0.6

0.4 0

0.2
0 1 2 3 4 5 6
-3
0 Tiempo (s)
x 10
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
Tiempo (s)
Figura 13. Sobretensión transitoria por energización de
Figura 10. Sobretensión por energización de línea abierta línea con impedancia de carga de 1000, 800 y 600 Ω

3
V - Fuente
V - ZL = 1 GOhm
2

1 1.5
Voltaje (p.u.)

0 1

-1 0.5
Voltaje (p.u.)

-2 0

-3 -0.5
0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04
Tiempo (s)
-1

Figura 11. Sobretensión transitoria por energización de

-1.5
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
línea con impedancia de carga de 10 𝐾Ω Tiempo (s)

Figura 14. Sobretensión transitoria por energización de

línea con impedancia de carga de 1000, 800 y 600 Ω

En los siguientes tres ejemplos de aplicación se

V - Fuente
V - ZL = 1000 Ohms V - ZL = 800 Ohms V - ZL = 600 Ohms

realizan simulaciones para el caso donde el

interruptor se cierra en un instante de tiempo
diferente de cero para la onda cosenoidal. En el
primer ejemplo, el cierre se realiza en un tiempo t

= 16 ms. En la Fig. 15 se muestra la sobretensión que se presenta ante la energización de la línea y

se observa que la señal no se atenúa debido a que 2 V - Cierre en t = 0.0123 s
V - Inicio de la línea.
1.5 V - Cierre en t = 0.016 s
la línea se encuentra en el vacío (sin carga). 1

0.5

Voltaje (p.u.)
0

-0.5

-1

-1.5
2
V - Fuente -2
V - ZL = 10 KOhms
1.5 0 0.005 0.01 0.015 0.02 0.025 0.03 0.035
Tiempo (s)
1

Figura 17. Comparación de sobretensión transitoria por

Voltaje (p.u.)

0.5

0 energización de línea abierta

-0.5

-1
Una descarga atmosférica tiene la peculiaridad de
-1.5
0 0.005 0.01
0.015
Tiempo (s)
0.02 0.025 0.03 0.035
ser irregular; es decir, en la realidad no tiene una
Figura 15. Sobretensión transitoria por energización de forma definida debido a que la tensión tiene
línea con impedancia de carga de 10 kΩ
múltiples variaciones durante el tiempo que se
En la Fig. 16 se muestra la actuación de la señal presenta, pero para su estudio teórico se dispone
al ir variando la carga. Este efecto se presenta de una forma de tensión normalizada, la cual es
debido a que los valores se aproximan al de la un impulso de tensión con un tiempo de subida
impedancia característica de la línea. Otra forma hasta el valor cresta de 1.2 𝜇𝑠 y un tiempo de cola
de evitar la sobretensión transitoria es realizar el de 50
cierre lo más cercano al cruce de la señal 𝜇𝑠, como se muestra en la Fig. 18 [7]. Las
cosenoidal con cero. La Fig. 17 muestra la sobretensiones que se presentan en el extremo
diferencia entre energizar un valor pico de la receptor de la línea de transmisión debido a una
señal con respecto a un valor cercano al cruce con descarga atmosférica se presentan en la Fig. 19.
cero, incluso estando la línea en el vacío. La
atenuación de la señal es mayormente
pronunciada.
Voltaje (KV)

V - Fuente
V - ZL = 1000 Ohms
V - ZL = 800 Ohms V - ZL = 600 Ohms

0.5
Voltaje (p.u.)

0
Tiempo (μs)

-0.5

-1
Figura 18. Forma de onda normalizada
0 0.005 0.01 0.015 0.02 0.025 0.03 0.035
Tiempo (s)

Figura 16. Sobretensión transitoria por energización de

línea con variación de impedancia de carga
 2

1.5 V - Decarga Atm. V - ZL = 10 KOhms presentan a través de su técnica gráfica, la cual ha

1

0.5
sido de las primeras en ser implementadas en
Voltaje (p.u.)

0
programas de simulación de tiempo EMTP y
-0.5

-1 sigue siendo utilizada en la actualidad. Se modeló

una línea de transmisión sin incluir las pérdidas y
-1.5

-2

la dependencia frecuencial en los parámetros

0 1 2 3 4 5 6
-3
Tiempo (s)
x 10

Figura 19. Sobretensión transitoria por descarga eléctricos. Sin embargo, la solución obtenida en
atmosférica en línea con carga de 10 𝐾Ω
los diferentes casos de simulación es consistente
V. CONCLUSIÓN y presenta el comportamiento esperado,
demostrando que aunque no se consideren dichos
En este trabajo se presentó el método de Bergeron
efectos, la técnica de Bergeron es confiable y de
aplicado al análisis de sobretensiones transitorias
fácil implementación para líneas aéreas,
en sistemas de transmisión monofásicos. El
obteniendo así resultados confiables.
método permite describir las sobretensiones que
se

BIBLIOGRAFÍA

[1] Harper, G. E. (s. f.). Elementos de Diseño de Subestaciones Eléctricas. México: Limusa.
[2] Paul, C. R. (1994). Analysis of Multiconductor Transmission Lines. New York: John Wiley & Sons
Inc.
[3] L. H. Delgado, G. C. y Aranda, F. C. (2008). Modelos de línea de transmisión para transitorios
electromagnéticos en sistemas de potencia. Energía y Computación.
[4] Harper, G. E. (2014). Sistemas de Transmisión y Distribución de Potencia Eléctrica. México: Limusa.
[5] Sluis, L. V. (2001). Transients in Power Systems. New York: John Wiley & Sons Inc.
[6] Stevenson, W. D. (1996). Análisis de Sistemas de Potencia. México: McGraw Hill.
[7] A. Greenwood. Electrical Transients in Power Systems (Second Edition). John Wiley & Sons Inc.
[8] L. A. Siegert C. (1989). Alta Tensión y Sistemas de Transmisión. México D. F.: Limusa.
[9] Electromagnetic Transient Program (EMTP-Theory Book). Portland, USA. July 1995.
[10] J. Martí. Accurate Modelling of Frequency-Dependent Transmission Lines in Electromagnetic
Transient Simulations. IEEE Trans. Power Apparatus and Systems, Vol. PAS-101, No. 1, January
1982.
[11] B. Gustavsen and A. Semlyen, Combined Phase Domain and Modal Domain Calculation of
Transmission Line Transients Based on Vector Fitting. IEEE Trans. Power Delivery, Vol. 13, No. 2,
April 1998.
[12] A. Morched, B. Gustavsen and M. Tartibi, A Universal Model for Accurate Calculation of
Electromagnetic Transients on Overhead Lines and Underground Cables. IEEE Trans. Power Delivery,
Vol. 14, No. 3, July 1999.