Universidad Central “Marta Abreu” de

Las Villas
Facultad de Ingeniería Eléctrica
Departamento de Electroenergética

TRABAJO DE DIPLOMA

Modelación en ATPDraw del comportamiento

ante rayos de líneas de transmisión usando
las técnicas de Monte Carlo.

Autor: Adrian Alfonso Figueroa González

Tutor: Dr. Angel C. Valcárcel Rojas

Santa Clara

2015

"Año 57 de la Revolución”
 Universidad Central “Marta Abreu” de

Las Villas
Facultad de Ingeniería Eléctrica
Departamento de Electroenergética

TRABAJO DE DIPLOMA

Modelación en ATPDraw del comportamiento

ante rayos de líneas de transmisión usando
las técnicas de Monte Carlo.
Autor: Adrian Alfonso Figueroa González
[email protected]

Tutor: Dr. Angel C. Valcárcel Rojas

[email protected]

Santa Clara

2015

"Año 57 de la Revolución”
Hago constar que el presente trabajo de diploma fue realizado en la Universidad
Central “Marta Abreu” de Las Villas como parte de la culminación de estudios de la
especialidad de Ingeniería en Automática, autorizando a que el mismo sea
utilizado por la Institución, para los fines que estime conveniente, tanto de forma
parcial como total y que además no podrá ser presentado en eventos, ni
publicados sin autorización de la Universidad.

Firma del Autor

Los abajo firmantes certificamos que el presente trabajo ha sido realizado según
acuerdo de la dirección de nuestro centro y el mismo cumple con los requisitos
que debe tener un trabajo de esta envergadura referido a la temática señalada.

Firma del Autor Firma del Jefe de

Departamento donde se
defiende el trabajo

Firma del Responsable de

Información Científico-Técnica
 i

PENSAMIENTO.

La más larga caminata comienza con un paso.

Proverbio hindú.
 ii

DEDICATORIA.

A mis padres y a Roig, mis resultados siempre son suyos.

 iii

AGRADECIMIENTOS.

A mi mamá por su cariño y su confianza, a mi padre por su guía y su apoyo, a

Roig por estar siempre, a mi familia por creer en mí, a mi tutor por poner a mi
disposición sus conocimientos y a mis amigos.
 iv

TAREA TÉCNICA

Tarea Técnica:

1. Realizar una búsqueda bibliográfica sobre la modelación en ATPDraw del

comportamiento ante rayos de líneas de transmisión.

2. Realizar una búsqueda bibliográfica sobre las técnicas de Monte Carlo.

3. Familiarizarse con el trabajo con el ATPDraw.

4. Implementar en ATPDraw las técnicas de Monte Carlo para la determinación

de la razón de salida de líneas de transmisión ante descargas atmosféricas.

5. Modelar en ATPDraw tres líneas de transmisión.

6. Comparar los resultados obtenidos con las estadísticas reales de las líneas
de transmisión modeladas.

Firma del Autor Firma del Tutor

 v

RESUMEN.

El estudio del comportamiento de las líneas de transmisión frente al rayo ha sido

objeto de gran atención durante los últimos años, todavía existen muchos puntos
sobre los que no hay un conocimiento suficientemente preciso; entre los más
importantes se encuentran los principales parámetros que describen su
comportamiento lo que plantea la necesidad de identificar fallas en la transmisión,
nos planteamos como objetivo general modelar en ATPDraw, el comportamiento
de líneas de transmisión frente descargas atmosféricas empleando el método de
Monte Carlo comparando los resultados obtenidos con el modelo desarrollado con
las estadísticas de falla y con los resultados de un estudio anterior de las líneas de
transmisión modeladas. Se implementó el Método de Monte Carlo en la
modelación en ATPDraw para determinar el índice de salida ante descargas
atmosféricas en líneas de transmisión de las que se conocían las estadísticas de
fallas acumuladas durante 10 años, obteniendo un error relativo inferior al de otros
autores. Se sugiere el Método de Monte Carlo como herramienta estadística útil
para predecir el comportamiento de las líneas de transmisión, futuras o ya
existentes, ante descargas atmosféricas y para futuros trabajos tener en cuenta la
contaminación en la zona de estudio lo que posibilitaría una mayor exactitud de los
resultados.
 vi

TABLA DE CONTENIDOS

PENSAMIENTO. ...................................................................................................................i

AGRADECIMIENTOS. ...................................................................................................... iii

TAREA TÉCNICA ...............................................................................................................iv

RESUMEN. ........................................................................................................................... v

INTRODUCCIÓN. ................................................................................................................ 8

CAPÍTULO 1. REVISIÓN BIBLIOGRÁFICA. ............................................................. 11

CAPÍTULO 2. MODELACIÓN EN ATPDRAW. ......................................................... 13

2.1 Modelación de los elementos de las líneas .................................................... 13

2.1.1 Introducción ................................................................................................... 13

2.1.2 Modelación de la línea ................................................................................ 13

2.1.3 Modelación de la estructura ............................................................................... 16

2.1.4 Modelación de la cadena de aisladores .......................................................... 18

2.1.5 Modelación de la resistencia de puesta a tierra ............................................. 19

2.1.6 Modelación de las fuentes de tensión a frecuencia de potencia y de

corriente de impulso tipo rayo. ......................................................................................... 21

2.1.7 Modelos generales ................................................................................................... 22

2.1.8 Modelo insulexa........................................................................................................ 23

2.2 Determinación del punto de impacto de la descarga ........................................ 26

2.2.1 Modelo electrogeométrico ...................................................................................... 26

 vii

2.3 Zonas de estudio de CFE .......................................................................................... 29

CAPÍTULO 3. RESULTADOS Y DISCUSIÓN ........................................................... 34

CONCLUSIONES Y RECOMENDACIONES. ............................................................... 40

CONCLUSIONES. ......................................................................................................... 40

RECOMENDACIONES. ................................................................................................ 40

REFERENCIAS BIBLIOGRÁFICAS................................................................................ 41
 CAPÍTULO 1. REVISIÓN BIBLIOGRÁFICA
8

INTRODUCCIÓN.

El propósito fundamental de un sistema eléctrico de potencia es abastecer de

energía eléctrica a sus clientes que esperan un suministro ininterrumpido
independientemente de las fluctuaciones de la demanda, este propósito no es
siempre alcanzable por las fallas de origen aleatorio a las que está expuesto.

Existen muchas causas por las que se pueden originar sobretensiones en una red
de potencia; por esta razón su estudio se suele realizar atendiendo al origen, al
tipo de proceso transitorio y al tiempo de duración. La clasificación más aceptada
distingue tres tipos de sobretensiones:

• Sobretensiones temporales

Dentro de este grupo están aquellas sobretensiones de larga duración (varios

milisegundos), poco amortiguadas y de frecuencia igual o próxima a la frecuencia
de operación. Ejemplos de sobretensiones temporales son las debidas a un
cortocircuito entre una fase y tierra, o a un problema de ferrorresonancia.

• Sobretensiones de maniobra

Una sobretensión de este tipo es fuertemente amortiguada, de corta duración y

puede presentar un rango de frecuencias que varía entre los 2 y los 10 kHz. Su
origen puede estar en una maniobra de conexión o de desconexión, sin embargo
puede haber otras causas que den lugar a una sobretensión de este tipo; por
ejemplo, un cortocircuito puede provocar transitoriamente una sobretensión que se
clasificaría dentro de este grupo.

• Sobretensiones de origen atmosférico

Son debidas a una descarga atmosférica, tienen una duración muy corta y una
amplitud que puede ser varias veces la tensión de pico nominal.

La importancia de las sobretensiones atmosféricas crece conforme disminuye la

tensión nominal de los componentes afectados por el rayo.
 CAPÍTULO 1. REVISIÓN BIBLIOGRÁFICA
9

Las descargas atmosféricas han sido uno de los principales problemas para el
diseño del aislamiento de sistemas de potencia y es todavía la causa principal de
salidas de servicio no programadas en líneas de transmisión.

Cuando las descargas atmosféricas impactan sobre líneas de transmisión

producen sobrevoltajes, los cuáles se calculan con mucha incertidumbre dada la
naturaleza aleatoria del rayo y el conocimiento poco preciso de sus parámetro Se
entiende por sobretensión una solicitación variable en el tiempo cuyo máximo valor
es superior al valor de pico de la tensión nominal del sistema en estudio.

El valor de las sobretensiones que se pueden producir en una red de transmisión,

originadas por un rayo es tan elevado frente al valor de la tensión nominal de la
red que tanto el nivel de aislamiento de los componentes como la selección y
coordinación de protecciones se realiza teniendo en cuenta el efecto de las
sobretensiones atmosféricas. En general, el nivel de las sobretensiones de origen
interno es inferior o muy inferior al nivel de aislamiento de los componentes de una
red de distribución, siendo el número de averías originadas por una sobretensión
de este tipo mucho más reducido que el originado por sobretensiones
atmosféricas directas o inducidas por el rayo. Un rayo puede originar una
sobretensión que termine provocando un fallo en la línea aérea afectada o que se
propague por la línea y pueda provocar una avería en algún otro equipo. En este
trabajo nos ocuparemos de la modelación de líneas de transmisión frente a
sobretensiones de origen atmosférico.

El cálculo de sobretensiones de origen atmosférico, y por tanto la evaluación del

comportamiento de líneas aéreas frente al rayo, puede realizarse con cierta
precisión mediante simulación digital y el empleo de modelos matemáticos
rigurosos para todos los componentes, teniendo como objetivo final determinar la
tasa de contorneos en por kilómetro de línea al año. La simulación de este
fenómeno puede ser fundamental no sólo en el análisis de líneas en explotación,
sino también en el diseño de nuevas líneas.
 CAPÍTULO 1. REVISIÓN BIBLIOGRÁFICA
10

En este trabajo se emplea el programa ATP (Alternative Transients Program), que

es un programa profesional para computador digital, utilizado internacionalmente
para la simulación de transitorios electromagnéticos, electromecánicos y funciones
de control en sistemas de potencia polifásicos, incluyendo los transitorios
asociados al efecto de las descargas atmosféricas en las líneas de transmisión
empleando el método de Monte Carlo, una técnica cuantitativa que hace uso de la
estadística y los ordenadores para imitar, mediante modelos matemáticos, el
comportamiento aleatorio de sistemas reales no dinámicos.

La transmisión de energía eléctrica debe llevarse a cabo de forma segura y al

menor costo posible, lo que plantea la necesidad de identificar fallas en la
transmisión, nos planteamos como objetivo general modelar en ATPDraw, el
comportamiento de líneas de transmisión frente descargas atmosféricas
empleando el método de Monte Carlo.

Para el logro del objetivo general se trazaron los siguientes objetivos específicos:

1. Implementar en ATPDraw las técnicas de Monte Carlo para la determinación de la

razón de salida de líneas de transmisión ante descargas atmosféricas.

2. Modelar en ATPDraw 16 líneas de transmisión de Comisión Federal de

Electricidad de México (CFM), de las cuales se dispone de las estadísticas de falla
acumuladas durante 10 años.

3. Comparar los resultados obtenidos con el modelo desarrollado con las estadísticas
de falla y con los resultados de un estudio anterior de las líneas de transmisión
modeladas.

Para alcanzar estos objetivos se realizó una revisión técnica y bibliográfica a

través de la búsqueda automatizada de información, la que permitió ampliar y
actualizar los conocimientos sobre la modelación en ATPDraw del comportamiento
ante rayos de líneas de transmisión y sobre las técnicas de Monte Carlo así como
familiarizarse con el trabajo con el ATPDraw.
 CAPÍTULO 1. REVISIÓN BIBLIOGRÁFICA
11

CAPÍTULO 1. REVISIÓN BIBLIOGRÁFICA.

Un sistema eléctrico de potencia está sometido a una serie de perturbaciones que

alteran su estado normal de operación. El paso de un estado a otro ocasionado
por una perturbación se hace en forma gradual, ya que las variables físicas como
voltajes y corrientes no pueden tener cambios bruscos debido a que las leyes de
la Física no lo permiten.

El estudio del estado transitorio debe hacerse en forma cuidadosa ya que las
variables físicas pueden llegar a tomar valores extremos y exigir a los equipos,
hasta ocasionar el deterioro de los mismos o la interrupción del suministro de
energía. El conocimiento de los estados transitorios permitirá tomar decisiones
para proteger adecuadamente los equipos, lo mismo que el aseguramiento de la
calidad de la potencia eléctrica, con unos estándares de calidad adecuados.

Los transitorios en un sistema de potencia en forma general son de tipo eléctrico,

mecánico y térmico. Los transitorios mecánicos y térmicos, a pesar de ser muy
importantes, son mucho más lentos que los de tipo eléctrico. Su frecuencia está
por debajo de la sincrónica (60 Hz).

Las líneas aéreas de transmisión de energía eléctrica se encuentran instaladas a

la intemperie, por lo que están expuestas a las descargas atmosféricas. Cuando
un rayo impacta una línea de transmisión genera un impulso de corriente, el cual a
su vez ocasiona una sobretensión transitoria. Esta sobretensión dependerá de la
magnitud del impulso de corriente del rayo, cuyo valor máximo esta usualmente en
el intervalo de 10kA a 100kA [1], pero también de la impedancia característica de
la línea y de las características de la torre de transmisión. Posterior al impacto de
la descarga sobre la línea. Se originan ondas viajeras hacia ambos extremos del
punto de impacto, las cuales serán distorsionadas y atenuadas conforme viajan
como consecuencia de las pérdidas propias de la línea y la conexión de
dispositivos de protección[2].Las descargas atmosféricas han sido uno de los
principales problemas para el diseño del aislamiento de sistemas de potencia y es
 CAPÍTULO 1. REVISIÓN BIBLIOGRÁFICA
12

todavía la causa principal de salidas de servicio no programadas en líneas de

transmisión [3]

En materia de fenómenos naturales y físicos; las descargas eléctricas

atmosféricas más comúnmente conocidas como rayos, son uno de los fenómenos
naturales con más impacto sobre la tierra, que constituyen la mayor de las
perturbaciones Electromagnéticas afectando el desempeño de los sistemas e
instalaciones eléctricas en todo el mundo, siendo la principal causa de salidas de
líneas de transmisión así como de daños en redes de distribución y equipos
electrónicos sensibles[4].

Un rayo produce sobre el aislamiento de una línea de transmisión un sobrevoltaje

que depende principalmente de la magnitud de la descarga y de la resistencia de
puesta a tierra. El parámetro que clásicamente ha medido dicho comportamiento
es el número de salidas de la línea por cada 100 km de longitud y por cada año,
parámetro que constituye un criterio de diseño del aislamiento, o una medida del
desempeño para una línea que ya esté construida. El rayo es un fenómeno que
origina una sobretensión de origen externo o simplemente de origen atmosférico.

Las sobretensiones de origen atmosférico son independientes de los voltajes de

operación de los sistemas de potencia a diferencia de las que se presentan por
maniobra, que si dependen del nivel de tensión nominal de los sistemas. La
importancia relativa de las sobretensiones atmosféricas crece a medida que
disminuye el nivel de tensión de los sistemas. Los sistemas de transmisión
normalmente son diseñados con un apantallamiento con una probabilidad de falla
muy pequeña. Las descargas atmosféricas que caerían en forma directa serían de
muy poca probabilidad. El flameo inverso ocasionado por una descarga que cae
sobre la torre o sobre el cable de guarda es la causa más común de salida de la
línea. Se denomina flameo inverso porque este se presenta del neutro hacia la
fase, debido al aumento del voltaje en el mismo.

Aunque los efectos de las descargas atmosféricas pueden ser muy severos, la
probabilidad de ocurrencia es relativamente baja (un kilómetro cuadrado de la
 CAPÍTULO 1. REVISIÓN BIBLIOGRÁFICA
13

superficie de la tierra es impactado por una descarga atmosférica al menos una

vez al año) [5]

En cualquier instante dado, casi 1,800 tormentas eléctricas están en progreso

sobre la superficie de la Tierra. La intensidad media de la descarga de un rayo se
estima en 31,000 amperios, pero se han detectado rayos de hasta 250,000
amperios. Las características climáticas y montañosas de cada país determinan el
número y la intensidad de las tormentas que se producen, riesgo que varía dentro
de un mismo país.

No existe método alguno para evitar la formación de descargas atmosféricas.

Tampoco sería deseable, en vista que los rayos son responsables en gran parte
de la formación de vida en el universo. El propósito entonces es tratar de
protegerse contra las descargas atmosféricas, controlando el paso de la corriente
de las descargas eléctricas y así prevenir daños a la propiedad [6].

En condiciones normales, existe en la atmósfera un equilibrio entre las cargas

positivas y negativas, en el que la tierra está cargada más negativamente que el
aire y los elementos situados sobre el suelo. Pero al formarse las nubes de
tormenta se produce una polarización de las cargas: en la mayoría de los casos, la
parte baja de las nubes queda cargada negativamente induciendo una carga
positiva en la tierra y los elementos situados sobre ella, formándose en la
atmósfera un campo eléctrico que llega a alcanzar decenas de kilovoltios.

Esta carga positiva se manifiesta especialmente en objetos metálicos, elementos

terminados en punta y objetos con una buena conexión a tierra, incluyendo los
árboles. Cuando el campo eléctrico es suficientemente intenso, la nube comienza
a descargarse hacia tierra. El camino que forma esta descarga se denomina
trazador descendente y produce una variación muy brusca del campo eléctrico,
que afecta a las cargas positivas de los objetos situados sobre el suelo,
produciéndose el denominado efecto corona. Uno de estos objetos será el que
forme el trazador ascendente, que irá a encontrarse con el trazador descendente,
quedando así ya formado el camino de la descarga entre la nube y tierra. Éste
 CAPÍTULO 1. REVISIÓN BIBLIOGRÁFICA
14

será el objeto que recibirá el impacto del rayo. Toda la carga de la nube buscará el
camino más directo hacia tierra, camino que, si no está controlado, puede causar
graves daños.

Cuando las descargas atmosféricas impactan sobre líneas de transmisión

producen sobretensiones, las cuáles se calculan con mucha incertidumbre dada la
naturaleza aleatoria del rayo y el conocimiento poco preciso de sus parámetros.

Se entiende por sobretensión una solicitación variable en el tiempo cuyo máximo

valor es superior al valor de pico de la tensión nominal del sistema en estudio.
Existen muchas causas por las que se pueden originar sobretensiones en una red
de potencia; por esta razón su estudio se suele realizar atendiendo al origen, al
tipo de proceso transitorio y al tiempo de duración. La clasificación más aceptada
distingue tres tipos de sobretensiones:

a) Sobretensiones temporales: Dentro de este grupo están aquellas

sobretensiones de larga duración (varios milisegundos), poco amortiguadas y de
frecuencia igual o próxima a la frecuencia de operación. Ejemplos de
sobretensiones temporales son las debidas a un cortocircuito entre una fase y
tierra, o a un problema de ferrorresonancia.

b) Sobretensiones de maniobra: Una sobretensión de este tipo es fuertemente

amortiguada, de corta duración y puede presentar un rango de frecuencias que
varía entre los 2 y los 10 kHz. Su origen puede estar en una maniobra de conexión
o de desconexión, sin embargo puede haber otras causas que den lugar a una
sobretensión de este tipo; por ejemplo, un cortocircuito puede provocar
transitoriamente una sobretensión que se clasificaría dentro de este grupo.

c) Sobretensiones de origen atmosférico: Son debidas a una descarga

atmosférica, tienen una duración muy corta y una amplitud que puede ser varias
veces la tensión de pico nominal. No existe una frontera muy clara entre un tipo y
otro de sobretensión; por ejemplo, una sobretensión originada por un cortocircuito
es de tipo temporal, pero transitoriamente puede ser calificada como sobretensión
 CAPÍTULO 1. REVISIÓN BIBLIOGRÁFICA
15

de maniobra. Por otra parte, la severidad que puede alcanzar cualquier

sobretensión depende del tipo y, sobre todo, del nivel de tensión de la red; en
redes de tensión nominal inferior a 400 kV las sobretensiones debidas al rayo son
más peligrosas que las debidas a otra causa, mientras que por encima de 400 kV
las sobretensiones de maniobra suelen ser las más peligrosas.

La importancia de las sobretensiones atmosféricas crece conforme disminuye la

tensión nominal de los componentes afectados por el rayo. El valor de las
sobretensiones que se pueden producir en una red de distribución, originadas por
un rayo, es tan elevado frente al valor de la tensión nominal de la red que tanto el
nivel de aislamiento de los componentes como la selección y coordinación de
protecciones se realiza teniendo en cuenta el efecto de las sobretensiones
atmosféricas. En general, el nivel de las sobretensiones de origen interno es
inferior o muy inferior al nivel de aislamiento de los componentes de una red de
distribución, siendo el número de averías originadas por una sobretensión de este
tipo mucho más reducido que el originado por sobretensiones atmosféricas
directas o inducidas por el rayo. Un rayo puede originar una sobretensión que
termine provocando un fallo en la línea aérea dada la naturaleza del rayo, el
análisis debe ser estadístico porque el número de combinaciones de las variables
involucradas en el proceso de cálculo es mucho menor respecto a las que se
necesitarían mediante un cálculo sistemático[7, 8].

El análisis del comportamiento de una línea aérea frente al rayo está basado en
métodos estadísticos, que son los que se adaptan mejor a la naturaleza de este
fenómeno, y tiene como objetivo final determinar la tasa de contorneos en una
línea por kilómetro y por año. La simulación de este fenómeno puede ser
fundamental no sólo en el análisis de fallas actuales, sino también en el diseño de
líneas futuras.

El método más empleado actualmente en el análisis estadístico de sobretensiones

en sistemas de potencia es el método de Monte Carlo. Un cálculo de
sobretensiones mediante el método de Monte Carlo es un procedimiento numérico
iterativo que utiliza en cada nuevo cálculo un conjunto de valores distintos que
 CAPÍTULO 1. REVISIÓN BIBLIOGRÁFICA
16

varían según la distribución de probabilidad asociada a cada una de las variables

involucradas en el proceso transitorio. El uso del método de Monte Carlo como
herramienta de investigación, proviene del trabajo realizado en el desarrollo de la
bomba atómica durante la segunda guerra mundial en el Laboratorio Nacional de
Los ´Álamos en EE.UU.

En 1930 Enrico Fermi y Stanislaw Ulam desarrollaron las ideas básicas del
método.

A principios de 1947 John Von Neumann envió una carta a Richtmyer a Los
Álamos en la que expuso de modo exhaustivo tal vez el primer informe por escrito
del método de Monte Carlo [9, 10].

Una de las primeras aplicaciones de este método a un problema determinista fue

llevada a cabo en 1948 por Enrico Fermi, Ulam y Von Neumann cuando
consideraron los valores singulares de la ecuación de Schrödinger [9, 10].

Para determinar la razón de salida de una línea aérea de potencia por el impacto
de descargas atmosféricas, lo más recomendado en la literatura es el uso del
Modelo Electrogeométrico (MEG), tanto para impactos directos [11-14], como para
impactos indirectos.

El cálculo riguroso de sobretensiones de origen atmosférico es un aspecto

importante en el análisis y diseño de líneas aéreas de transmisión. Las diferentes
partes de una línea que están involucradas en las sobretensiones causadas por el
rayo deben ser representadas teniendo en cuenta el rango de frecuencia de los
transitorios causados por el impacto de una descarga atmosférica. Por otro lado,
los procedimientos a emplear en el cálculo de sobretensiones deben ser
desarrollados teniendo en cuenta la naturaleza aleatoria del rayo. Varias
herramientas de cálculo han sido empleadas hasta la fecha para estimar el
comportamiento de líneas aéreas de transmisión frente al rayo. Los
procedimientos más utilizados emplean una técnica basada en el dominio del
tiempo para la que se han de desarrollar y aplicar modelos adecuados de las
 CAPÍTULO 1. REVISIÓN BIBLIOGRÁFICA
17

distintas partes de una línea. Este artículo presenta un resumen del trabajo
realizado por el autor para desarrollar e implantar en una herramienta de
simulación tipo EMTP (ElectroMagnetic Transients Program) un procedimiento
basado en el método de Monte Carlo que permita analizar el comportamiento
frente al rayo de líneas aéreas de transmisión. Los resultados presentados en este
documento se basan en 18 líneas reales cuyo análisis permitirá ilustrar el alcance
del procedimiento y el tipo de estudios que se puede realizar con el procedimiento
desarrollado [11].EI programa ATP (Electromagnetic Transients Program), es un
programa para computador digital, utilizado para la simulación de transitorios
electromagnéticos, electromecánicos y funciones de control en sistemas de
potencia polifásicos[15, 16].

EI ATP fue desarrollado en su versión inicial por el Dr. Hermann Dommel a finales
de la década del 60 en Alemania (inicialmente con el nombre de EMTP), quien
posteriormente vendería los derechos a la Bonneville Power Administration (BPA)
de los Estados Unidos. A pesar de ser la BPA la que coordinara la distribución del
programa entre los diferentes grupos de usuarios, muchas universidades
contribuyeron al desarrollo de los diferentes modelos de este programa[17]. A
pesar de que el objetivo principal del programa es la obtención de la respuesta en
el tiempo de las diferentes variables electromecánicas del sistema bajo estudio,
también se puede obtener la respuesta en estado estacionario para un sistema
alimentado con corriente alterna.

El programa ATP en el ambiente de trabajo normal viene acompañado de otros

programas, como el ATPDRAW que permite realizar un diagrama de la red
eléctrica a ser analizada, para posteriormente interactuar con ATP [18, 19]. Los
modelos de líneas de la transmisión disponibles en el programa ATP son bastante
flexibles y los hay para cada necesidad de simulación de transitorios. Una línea en
forma general se puede simular mediante parámetros distribuidos o mediante
parámetros semidistribuidos.
 CAPÍTULO 1. REVISIÓN BIBLIOGRÁFICA
18

El programa permite la representación de resistencia e inductancias de tipo no

lineal para equipos como pararrayos, núcleos magnéticos con características de
saturación, resistencias no lineales, etc.

El programa ATP tiene una gama amplia de modelos de interruptores. Existen

modelos de interruptores controlado por tiempo, controlados por voltaje, de
corriente unidireccional (diodo), de corriente unidireccional controlado por señales
de compuerta, controladas por señal externa (por fuera de la red de potencia),
interruptores de tipo estadístico, interruptores de tipo sistemático.

Los interruptores estadísticos y los sistemáticos se usan para simular el

comportamiento aleatorio que tiene un interruptor en el tiempo de cierre ya que el
tiempo de cierre presenta una dispersión alrededor de un tiempo promedio
especificado por el fabricante.

Los interruptores controlados por tensión son utilizados para simular el

comportamiento de un gap o de un aislador, el cual presenta disrupción cuando
una tensión entre sus extremos es superada. En este tipo de interruptor se
controla el tiempo a partir del cual puede entrar en conducción y el tiempo mínimo
que debe permanecer cerrado.

El programa permite la representación de fuentes de excitación, de voltaje o de

corriente que están definidas analíticamente dentro del programa.

Estas fuentes por si solas tienen un comportamiento ideal, es decir que si se

requiere un modelo de la de fuente se debe agregar el correspondiente
equivalente de Thevenin o Norton según sea fuente de tensión o de corriente. Las
fuentes por definición se conectan entre la tierra y un nodo.

Cuando el rayo impacta directamente la línea eléctrica, ya sea en la torre, en el

vano entre torres o en los conductores de fase se habla de “impactos directos” .
Cuando el impacto se produce en la torre o en los conductores de blindaje en el
vano entre torres puede ocurrir lo que se conoce como una “descarga inversa”,
 CAPÍTULO 1. REVISIÓN BIBLIOGRÁFICA
19

provocando una circulación de la corriente de rayo que penetra en la estructura

hacia los conductores de fase.

Si el impacto se produce en los conductores de fase, ocurre una “falla de blindaje”,

la descarga que se puede desarrollar se conoce como “descarga directa”,
produciéndose ésta desde los conductores de fase hacia la estructura.

En el presente trabajo se propone un método de determinación de la razón de

salida de una línea de transmisión, aplicando las técnicas de Monte Carlo,
mediante el cual generan 1420 impactos de rayo, lo cual representa lo que
sucedería en una línea de 100 km de longitud, que corre por una zona con una
densidad de rayos a tierra de 4.375 rayos por km 2 al año, durante 10 años (la
misma cantidad de años de los que se tiene estadísticas de las líneas modeladas).
 13
CAPÍTULO 2. MODELACIÓN EN ATPDRAW

CAPÍTULO 2. MODELACIÓN EN ATPDRAW.

2.1 Modelación de los elementos de las líneas

2.1.1 Introducción

Tomando como punto de partida algunos de los resultados obtenidos por John G.
Anderson [20, 21] y los obtenidos por A. Valcárcel [22] se desarrolló una modelación
en ATPDRAW [23] que permiten determinar la ocurrencia o no de descargas inversas
o directas en líneas de transmisión.

Con la modelación en ATPDRAW, se pueden aprovechar las prestaciones de este

programa y desarrollar los módulos necesarios para representar todos los
componentes involucrados.

A continuación se presenta un resumen de las características y objetivos de los

modelos desarrollados para llevar a cabo este tipo de cálculos.

En este caso se hace una representación de la línea, torres, cadenas de aisladores,

resistencia de puesta a tierra, fuentes de tensión a frecuencia de potencia y de
corriente de impulso.

2.1.2 Modelación de la línea

De acuerdo con las recomendaciones del Grupo de Trabajo 33.02 de CIGRE, en el

caso de sobretensiones debidas al rayo se trata de transitorios de frente rápido, por lo
que la representación de una línea aérea se debe hacer, de acuerdo con la Tabla 1
 14
CAPÍTULO 2. MATERIALES Y MÉTODOS
 15
CAPÍTULO 2. MATERIALES Y MÉTODOS

La representación debe extenderse a dos o tres vanos desde el punto de impacto de

la descarga atmosférica.
 Se debe incluir la representación de las estructuras y de sus impedancias de
puesta tierra.
 Es necesario desarrollar un modelo que evite reflexiones de ondas en
cualquiera de los extremos del tramo de línea representado; esto se puede
conseguir adaptando la línea en los dos extremos mediante la instalación de
una matriz de resistencias de valor igual a la impedancia característica, o
haciendo que las dos secciones de los extremos tengan una longitud tal que,
con el tiempo de simulación escogido, ninguna onda reflejada en los extremos
abiertos alcance la parte de la línea que es de interés.
 El modelo será, evidentemente, el de una línea no transpuesta pues se
analizará solo un tramo de la línea en estudio de entre 5 y 7 estructuras.
 Es aconsejable incluir el efecto de la tensión a frecuencia de operación.

Puesto que las tensiones más elevadas se producirán en el punto de impacto, sólo
será necesario considerar los apoyos cercanos a este punto. La Figura 2.1.2.1
muestra la representación escogida para la línea en ATPDRAW.

Fig. 2.1.2.1. Modelo que considera dos torres adyacentes más un tramo de línea
lo suficientemente grande como para que las reflexiones en los extremos no
afecten el análisis.

El modelo utilizado en ATPDRAW se muestra en la Figura 2.1.2.2.

 16
CAPÍTULO 2. MATERIALES Y MÉTODOS

Fig. 2.1.2.2. Modelo del ATPDRAW utilizado para representar los tramos de línea
entre estructuras.

2.1.3 Modelación de la estructura

Se puede suponer que la torre es una simple línea corta de transmisión con una
impedancia transitoria constante, aterrada mediante su resistencia de aterramiento a
la base. Esta impedancia transitoria se ha calculado y medido tanto en el campo como
en modelos de escala y que pueden ser idealizadas por un cono de altura h y un radio
de base r. Las torres convencionales de doble-circuito, poseen una impedancia
transitoria, ampliamente aceptada por la literatura [24, 25], que se expresa mediante:
 17
CAPÍTULO 2. MATERIALES Y MÉTODOS

 h2 
Z t  30 ln 21  2 
 r 
 2.1

Y para torres cilíndricas, la literatura [4][60][61] plantea que la expresión siguiente da

resultados satisfactorios:

h r
Z t  60 ln    90    60  2.2
r h

dónde h y r son la altura y el radio equivalente (significa la periferia dividida por 2 )

de la torre.

En el 2004, A. Soares Jr., S. Visacro [26], presentan un trabajo donde presentan una
nueva forma de análisis de impedancia característica para torres de transmisión.
Tomando la base cuadrada de la torre actual. Primero, consideran un círculo interior,
y en segundo, un círculo exterior a esta base cuadrada de la torre. A este método lo
llamaron HEM-Hybrid Electromagnetic Model, proponiendo nuevas fórmulas en esta
investigación. Posteriormente hacen una simulación con estas consideraciones y con
la ayuda del EMTP se lleva a cabo una comparación entre los método tradicionales:
el cilindro presentado por Wagner y Hileman [27] y la forma de cono por Sargent y
Darveniza [25] en los que encuentran desviaciones arriba de 22% y 13%
respectivamente.

En el desarrollo de este trabajo el autor adopta las formulaciones desarrolladas por

Sargen y Darveniza [25] por considerar que tienen una mayor validación práctica que
el resto de las formulaciones dada la extensa bibliografía en que son mencionados.

Con esta expresión de impedancia transitoria, se desarrolló un modelo de la torre en

ATPDRAW, véase la Figura 2.1.3.1.
 18
CAPÍTULO 2. MATERIALES Y MÉTODOS

Fig. 2.1.3.1. Representación de la estructura en el ATPDRAW.

2.1.4 Modelación de la cadena de aisladores

Los aisladores se representan mediante interruptores controlados como se puede ver

en la Figura 2.1.4.1.

Fig. 2.1.4.1. Representación de la cadena de aisladores en el ATPDRAW.

 19
CAPÍTULO 2. MATERIALES Y MÉTODOS

2.1.5 Modelación de la resistencia de puesta a tierra

La resistencia de puesta a tierra se puede modelar [28, 29] por una resistencia
constante e igual a la resistencia de estado estable [30] o por una resistencia variable
[31] con la corriente que penetra en el terreno, mostrado en la Figura 2.1.5.1.

Fig. 2.1.5.1. Representación de la resistencia de puesta a tierra en el ATPDRAW.

a) Resistencia constante e igual a la resistencia de estado estable, b)
Resistencia dinámica, variable con la corriente que penetra en el terreno.

La resistencia dinámica tiene un comportamiento que depende de la corriente que

penetra en el terreno, de la resistencia de estado estable y la resistividad del terreno.

R0
R 3.3
 
1  I 1 
 I 
 g 

E 0
Ig  3.4
2R02

Donde: R0 = Resistencia de estado estable (Ω)

E0 = Campo eléctrico de ionización del suelo (kV/cm)

ρ = Resistividad del terreno (Ω.cm)

I1 = Corriente que penetra en el terreno (A)

 20
CAPÍTULO 2. MATERIALES Y MÉTODOS

IG = Corriente necesaria para producir el gradiente eléctrico de potencial

E0 al que ocurre la avería a tierra (A)

En las Figuras 2.1.5.2 y 2.1.5.3 se muestra la variación de la corriente inyectada en el

terreno y de la resistencia de conexión a tierra en función de esta corriente.

Figura 2.1.5.2. Variación en el tiempo de la corriente inyectada al terreno.

Figura 2.1.5.3. Variación en el tiempo de la resistencia de conexión a tierra en

función de la corriente de la figura anterior.
 21
CAPÍTULO 2. MATERIALES Y MÉTODOS

2.1.6 Modelación de las fuentes de tensión a frecuencia de potencia y de

corriente de impulso tipo rayo.

Para la modelación de la fuente de alimentación se utilizó un modelo trifásico a la cual

se le genera aleatoriamente el ángulo al cual se produce el impacto, tal como muestra
la Figura 2.1.6.1.

La fuente de impulso se representa mediante una fuente controlada por el model

insulexa, a la cual se le genera la onda de corriente aleatoriamente tanto en valor
máximo como en tiempo de frente, además el model determina si el impacto se
produce en el tope de la torre o en el vano entre torres y cierra los interruptores
Close1 (Tope) o Close2 (Vano), ver la Figura 2.1.6.2.

Fig. 2.1.6.1. Fuente de tensión a frecuencia de potencia.

 22
CAPÍTULO 2. MATERIALES Y MÉTODOS

Fig. 2.1.6.2. Fuente de corriente de impulso.

2.1.7 Modelos generales

Para la representación del impacto en el tope o el vano de la torre se utiliza el modelo

de la Figura 2.1.7.1
 23
CAPÍTULO 2. MATERIALES Y MÉTODOS

Fig. 2.1.7.1. Modelo para el impacto en el tope y el vano de la torre.

2.1.8 Modelo insulexa

El método de Monte Carlos se implementa en el modelo insulexa de la figura 2.1.8.1

Los datos de entrada del modelo son:

Td: días tormenta al año

tower_h: altura de los hilos de guarda m

shw_sep: separación entre los hilos de guarda

clearance: distancia de la cadena de aisladores

 24
CAPÍTULO 2. MATERIALES Y MÉTODOS

sigma_str:

imean mean:

sigma_gap:

front_st: datos de la onda del rayo

tstart:

tail:

ksm:

Lv: distancia entre las torres

Entradas:

vph_a: tensión en el conductor de la fase A

vph_b: tensión en el conductor de la fase B

vph_c: tensión en el conductor de la fase C

arm_a:

arm_b tensión en la torre

arm_c

knt: número de descargas atmosféricas

tmax: tiempo de muestreo

vsh_h: tensión en el tope

 25
CAPÍTULO 2. MATERIALES Y MÉTODOS

Salidas:

bfla:

bflb: determina cual de las fases falla

bflc:

ltn: onda de la corriente del rayo

pot_a:

pot_b: tensión en las fases

pot_c:

resis: resistencia de puesta a tierra de la estructura

istrok: valor máximo de onda del rayo

Fig. 2.1.8.1. Modelo insulexa

 26
CAPÍTULO 2. MATERIALES Y MÉTODOS

2.2 Determinación del punto de impacto de la descarga

2.2.1 Modelo electrogeométrico

Basado en el modelo electrogeométrico (MEG), el ancho del área de atracción de una

línea aérea viene dado como se muestra en la Figura 2.2.1.1 por [32]:

A  b 2x (m) 2.2

o lo que es lo mismo:


A  b  2 rs2  rsg  h 
2
(m) 2.3

siendo:

rs  8  I s0.65 (m) 2.4

rsg  k sg  rs (m) 2.5

rs = Distancia de impacto (m)

rsg = Distancia de impacto al plano de tierra (m)

I s = Corriente de retorno (kA)

k sg = Entre 0.6 y 1 según el criterio de diferentes autores dependiendo del nivel de

tensión y la altura de la línea.

h= Altura de los cables protectores de la línea en caso de que esté desprotegida.

b = Separación entre los cables protectores para una línea protegida con dos cables o
la separación entre las fases extremas de una línea desprotegida e igual a cero para
una línea protegida con un solo cable protector.
 27
CAPÍTULO 2. MATERIALES Y MÉTODOS

La razón de salida total de la línea es la suma de las salidas por impactos en la torre,
por impactos en el vano entre torres y por fallas de blindaje.

Según lo planteado anteriormente, el área total de impacto en una línea está dado
por:

A
Atotal  L (km2) 2.6
1000

Donde: L = longitud de la línea (km)

rs
rsg-h rsg

x
h

Figura 2.2.1.1. Área de atracción a ambos lados de la línea.

El área de impacto en las torres será, según la Figura 2.2.1.2:

Atorres  At Nt (km2) 2.7

Donde: At = área de atracción de una torre (km2)

Nt = número de torres de la línea

Ax
At  (km2) 2.8
2

x  rs2  rsg  h 
2
(km) 2.9
 28
CAPÍTULO 2. MATERIALES Y MÉTODOS

Figura 2.2.1.2. Áreas de impacto en una línea.

El área de impacto en el vano entre torres se determina como la diferencia entre el

área total (Atotal) y el área de las torres (Atorres):

Avano  Atotal  Atorres (km2) 2.10

Una situación especial se presenta para valores de corriente de rayo donde las áreas
de atracción de las torres se solapan, la Figura 2.2.1.3 muestra esta situación; en este
caso, al área de las torres, calculada anteriormente, se le resta las áreas de
solapamiento:

Atorres2  Atorres1  ( N t  1) Asolapamiento (km2) 2.11

Donde: Atorres1 = área de las torres sin solapamiento (km2)

Atorres2 = área de las torres con solapamiento (km2)

Además:

  
   L2  
  A 2 1  v 

       4 rs2 
1
Asolapamiento   2 Lv A2  4 
L2v  
 4 L2 2 L2
v   ArcTan   
 2   v  A 4     2.12
8
 rs    rs2   Lv 
  
  
  

Donde: LV = longitud del vano entre torres (km)

 29
CAPÍTULO 2. MATERIALES Y MÉTODOS

Figura 2.2.1.3. Solapamiento del área de atracción de las torres

2.3 Zonas de estudio de CFE

Para el análisis del comportamiento ante rayos se seleccionaron dos zonas del
territorio mexicano con igual densidad de rayos a tierra figuras 2.3.1.

Fig. 2.3.1. Mapa de densidad de rayos a tierra.

En cada zona se seleccionaron 8 líneas. Las líneas seleccionadas y sus

características se muestran en la Tabla 2.3.1.
 30
CAPÍTULO 2. MATERIALES Y MÉTODOS

En ambas zonas se seleccionaron líneas con dos tipos de estructuras la 1A11 y la

1YR21, cuyas características se muestran en las Figuras 2.3.3 y 2.3.4
respectivamente.

Las estadísticas de falla de cada una de las líneas, según CFE [33], se muestran en
las Tablas 2.3.2 y 2.3.3.

Fig. 2.3.3. Estructura tipo 1A11. Fig. 2.3.4. Estructura tipo 1YR21.
 CAPÍTULO 2. MATERIALES Y MÉTODOS
31

Tabla 2.3.1. Datos característicos de las líneas seleccionadas.

 CAPÍTULO 2. MATERIALES Y MÉTODOS
32

Tabla 2.3.2. Estadísticas de falla por descargas atmosféricas de las líneas de la zona Baja California.

AÑO Salidas Tipo

LÍNEA IFSLT
1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 totales Torre

INSURGENTES-LORETO 3 2 4 2 1 2 2 3 2 3 24 2.40 1A11

MEXICALI II-RUIZ
1 3 2 1 0 1 1 1 1 1 12 1.20 1A11
CORTINES

PUNTA PRIETA II-V.

5 2 7 2 2 2 7 3 4 2 36 3.60 1YR21
CONSTITUCIÓN

PUNTA PRIETA II-V.

4 4 5 3 2 3 3 4 3 4 35 3.50 1YR21
CONSTITUCIÓN

PUNTA PRIETA II-EL

1 1 2 1 2 1 1 2 1 1 13 1.30 1YR21
TRIUNFO

SANTIAGO-SAN JOSE 2 1 0 1 0 1 2 2 1 1 11 1.10 1YR21

TRIUNFO-SANTIAGO 2 1 2 0 2 0 1 2 1 1 12 1.20 1YR21

V. CONSTITUCIÓN-
0 1 1 2 2 1 2 1 1 1 12 1.20 1YR21
SANTO DOMINGO
 CAPÍTULO 2. MATERIALES Y MÉTODOS
33

Tabla 2.3.3. Estadísticas de falla por descargas atmosféricas de las líneas de la zona Peninsular.

AÑO SAL. Tipo

LÍNEA IFSLT
1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 TOT Torre

LERMA-CHAMPOTON 0 1 0 1 1 2 1 1 3 1 11 1.1 1YR21

LERMA-MAXCANU 1 3 4 3 2 2 1 3 4 3 26 2.6 1YR21

LERMA-MÉRIDA II 0 3 5 3 4 3 4 3 4 5 32 3.2 1YR21

LERMA-TICUL II 3 3 6 2 6 3 5 4 5 2 39 3.9 1YR21

TICUL II-KAMBUL 1 5 6 4 4 5 5 2 4 3 39 3.9 1A11

TIZIMÍN-CANCÚN 0 1 2 3 3 3 4 4 5 6 31 3.1 1YR21

VALLADOLID-NIZUC 2 3 3 4 3 1 5 4 3 6 34 3.4 1A11

VALLADOLID-PLAYA
DEL CARMEN 3 7 3 3 7 2 3 5 2 6 41 4.1 1A11
 CAPÍTULO 3. RESULTADOS Y DISCUSIÓN
34

CAPÍTULO 3. RESULTADOS Y DISCUSIÓN

A continuación se muestra un resumen a partir de los resultados obtenidos en el

ATPDraw.

El índice de salida real de las líneas, el hallado por otros autores [34] y el que se
calcula a través de la implementación del método de Monte Carlo así como el error
relativo se muestran en las tablas 3.1 y 3.2, como se puede apreciar existe
coincidencia con las estadísticas y los errores relativos son inferiores a los que se
obtienen en [34]

Tabla 3.1 Baja California.

IFSL de Error de Tipo de

LÍNEA IFSLT(MC) IFSLT Error(MC)
Alfredo Alfredo Torre

INSURGENTES-LORETO 2.60 2.40 7.69 2.7244 4.41 1A11

MEXICALI II-RUIZ CORTINES 1.40 1.20 14.28 1.4976 19.46 1A11

PUNTA PRIETA II-V.

3.79 3.60 5.01 3.2572 10.76 1YR21
CONSTITUCIÓN

PUNTA PRIETA II-V.

3.86 3.50 9.32 3.3156 5.74 1YR21
CONSTITUCIÓN

PUNTA PRIETA II-EL

1.10 1.30 18.18 0.9499 38.29 1YR21
TRIUNFO

SANTIAGO-SAN JOSE 0.83 1.10 32.53 0.7162 54.92 1YR21

TRIUNFO-SANTIAGO 1.17 1.20 2.56 1.0083 19.01 1YR21

V. CONSTITUCIÓN-SANTO
0.91 1.20 31.86 0.7893 53.84 1YR21
DOMINGO
 CAPÍTULO 3. RESULTADOS Y DISCUSIÓN
35

Tabla 3.2 Zona Peninsular.

IFSLT de Error de Tipo de

LÍNEA IFSLT(MC) IFSLT Error(MC)
Alfredo Alfredo Torre

LERMA-CHAMPOTON 0.97 1.10 13.4 1.8343 39.89 1YR21

LERMA-MAXCANU 2.23 2.60 16.59 4.2147 38.24 1YR21

LERMA-MÉRIDA II 2.96 3.20 8.1 5.5979 28.41 1YR21

LERMA-TICUL II 2.36 3.90 65.25 4.4720 12.75 1YR21

TICUL II-KAMBUL 2.60 3.90 50 4.0631 3.94 1A11

TIZIMÍN-CANCÚN 2.50 3.10 24 4.7293 34.32 1YR21

VALLADOLID-NIZUC 3.85 3.40 11.68 6.0027 43.35 1A11

VALLADOLID-PLAYA
3.90 4.10 5.12 6.0758 32.51 1A11
DEL CARMEN

Los datos generales del comportamiento de las líneas de transmisión ante rayos
se muestran en la tabla 3.3

Tabla3.3 Comportamiento ante rayos de líneas de transmisión usando el

método de Monte Carlo

Fallas por Impactos Fallas por Fallas Fallas Fallas

Impactos Tipo de
impacto en en la impacto en la en la en la
en el vano estructura
el vano torre en la torre fase A fase B fase C

903 6 474 29 35 0 0 1A11

747 1 422 22 16 7 0 TYR21

 CAPÍTULO 3. RESULTADOS Y DISCUSIÓN
36

Las tablas 3.4 y 3.5 muestran las magnitudes de corriente y resistencia a tierra
para ambos tipos de estructura 1A11 y TYR21 respectivamente.

Tabla3.4 1A11

Corriente del Fase Resistencia de la

KNT Impacto
rayo (kA) fallada torre(Ω)
73 104 A 14 T
113 86 A 18 T
178 159 A 11 T
258 191 A 11 T
289 96 A 11 V
310 139 A 13 V
342 209 A 13 V
363 241 A 17 T
420 91 A 14 T
425 115 A 19 T
426 173 A 13 T
477 167 A 15 T
510 123 A 18 T
525 185 A 18 T
536 97 A 12 T
540 149 A 15 V
798 116 A 16 T
806 95 A 17 T
819 111 A 15 V
858 146 A 15 T
880 113 A 17 T
902 125 A 10 T
909 156 A 11 T
1000 106 A 16 T
1044 182 A 14 T
1047 147 A 13 T
1067 194 A 16 T
1099 154 A 19 V
1180 115 A 18 T
1221 105 A 18 T
1245 132 A 14 T
1280 77 A 19 T
1284 96 A 13 T
1293 187 A 14 T
1304 128 A 17 T
 CAPÍTULO 3. RESULTADOS Y DISCUSIÓN
37

Tabla 3.5 TYR21

Corriente
Fase Resistencia de
KNT del rayo Impacto
fallada la torre(Ω)
(kA)
69 116 A 18 T
178 159 B 11 T
258 191 A 11 T
342 209 A 13 V
363 241 A 17 T
425 115 A 19 T
426 173 A 13 T
477 167 A 15 T
510 123 B 18 T
525 185 B 18 T
540 149 A 15 T
798 116 B 16 T
858 146 A 15 T
880 113 A 17 T
1000 106 A 16 T
1044 182 A 14 T
1047 147 A 13 T
1067 194 A 16 T
1180 115 A 18 T
1221 105 B 18 T
1245 132 B 14 T
1293 187 A 14 T
1304 128 B 17 T

Las figuras 3.1 y 3.3 muestran la distribución de probabilidad de corrientes de rayo

que hacen fallar el aislamiento de las líneas, para las estructuras 1A11 y TYR21
respectivamente.

Las figuras 3.2 y 3.4 muestran la distribución de probabilidad de la resistencia a

tierra de las estructuras bajo estudio 1A11 y TYR21 respectivamente.
 CAPÍTULO 3. RESULTADOS Y DISCUSIÓN
38

300

250

200

150

100

50

0
1 3 5 7 9 11 13 15 17 19 21 23 25 27 29 31 33 35

Figura 3.1 Distribución de probabilidad de magnitud de corriente del rayo de

la estructura 1A11

20
18
16
14
12
10
8
6
4
2
0
1 3 5 7 9 11 13 15 17 19 21 23 25 27 29 31 33 35

Figura 3.2 Distribución de probabilidad de resistencia a tierra de la

estructura 1A11
 CAPÍTULO 3. RESULTADOS Y DISCUSIÓN
39

250

200

150

100

50

0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17

Figura 3.3 Distribución de probabilidad de magnitud de corriente del rayo de

la estructura TYR21

20

18

16

14

12

10

0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20

Figura 3.4 Distribución de probabilidad de resistencia a tierra de la

estructura TYR21
 CONCLUSIONES Y RECOMENDACIONES
40

CONCLUSIONES Y RECOMENDACIONES.

CONCLUSIONES.

Se implementó el Método de Monte Carlo en la modelación en ATPDraw para

determinar el índice de salida ante descargas atmosféricas en 16 líneas de
transmisión de CFM, de las que se conocían las estadísticas de fallas acumuladas
durante 10 años, existiendo una gran coincidencia con las estadísticas reales y
con un error relativo inferior al de otros autores.

RECOMENDACIONES.

Perfeccionar el Método de Monte Carlo implementado para tener en cuenta la

contaminación en la zona de estudio lo que posibilitaría una mayor exactitud de los
resultados.
 REFERENCIAS BIBLIOGRAFICAS
41

REFERENCIAS BIBLIOGRÁFICAS

[1] I. S. 4. (1978 Modificación de la Norma IEEE Técnicas de Pruebas de alta tensión.

[2] T. I. U. García Lagos. FA, "Estudio estadístico de sobretensiones transitorias debidas a
descargas atmosféricas en líneas de transmisión. ," 2012.
[3] Y. F. Ruales L. Cálculo del Comportamiento de una Línea de Transmisión frente al Flameo
Inverso basado en el Método Montecarlo. .
[4] M. J. Grisales Garcia V. (2013, Diseño de un sistema de protección contra
descargaseléctricas atmosféricas: caso bloque e de la universidad tecnológica de Pereira.
Universidad Tecnológica de Pereira. Facultad de Tecnología. .
[5] Available: http://vmwl1.iie.org.mx/sitioIIE/sitio/indice.php
[6] M. G. L. Angeles Sánchez JM. (2010, Implementación del método de volumen de colección
para la protección contra descargas eléctricas atmosféricas en la TDGL Poza Rica
[7] G. M. F., "Evaluación estadística del comportamiento de líneas aéreas de distribución
frente a sobretensiones de origen externo. ," Tesis doctoral.
, Barcelona. , Barcelona. , Barcelona. , 2001.
[8] A. S. F., "Fenómeno del rayo. Gestión de riesgos. ."
[9] Y. G. FG. (2010, Cálculo del comportamiento de una línea de transmisión frente al flameo
inverso basado en el método de Monte Carlo. .
[10] G. P. G. Fuentes Rodríguez M, "Método de Monte Carlo. ," 2010.
[11] C.-A. F. Martínez-Velasco JA, "Implantación en EMTP de un procedimiento basado en el
método de Monte Carlo para analizar el comportamiento de líneas de transmisión frente
al rayo Revista chilena de ingeniería.," vol. vol. 16 Nº 1, pp. 169-180 2008.
[12] R. B. a. A. J. E. Anderson, "Lightning parameters for engineering applications," 1980.
[13] F. Elizarraraz R., VillalobosR., A., "Ensayo alternativo de tensión de perforación para
aisladores de vidrio y porcelana," 2002.
[14] D. Qi-Bin Zhou, Y., "Using EMTP for evaluation of surge current distribution in metallic
gridlike structures," 2005.
[15] C. C.L. (1995). Teoría y práctica con el ATP. .
[16] R. C. L. Vázquez O. K, "Cálculo de sobrevoltajes originados por descargas atmosféricas en
lineas de transmisión utilizando el ATP (Alternative transcients program) aplicado a la línea
de transmisión Santa Rosa- Pomasqui a 220kV. ," 2005.
[17] A. Br. Bande, "Estudio de las sobretensiones por descargas atmosféricas directas en las
líneas aéreas de distribución de 13.8kV de os distritos Múcura y san Tomé de PDVSA.
Maracay. ," 2008.
[18] H. H. Prikler L. ( 1998 ATPDraw for Windows 3.1x/95/NT version 1.0.
[19] G. C. H. A. Villa R, "Simulación con ATP de pararrayos para las líneas a 115 y 400 kVde
Edelca. Rev Ciencia e Ingeniería. ," No 1., vol. Vol 22, 2001.
[20] A. J. G., "Transmission Line Reference Book," pp. Pag 545–595, 1987.
[21] J. G. a. T. A. S. Anderson, "Algorithms for calculation of lightning induced voltages on
distribution lines”, presented at the 1992 summer power meeting for publication in the
IEEE Trans. on Power Delivery," 1992.
[22] M. A. Sargent, "The Frequency Distribution of Current Magnitudes of Lightning Strokes to
Tall Structures”, IEEE Transactions on Power Apparatus and Systems
Volume PAS-91," 1972.
 REFERENCIAS BIBLIOGRAFICAS
42

[23] T. M. Mekhaldi A, "Algorithm for pollution flashover calculation," 1996.

[24] B. K. Gautam, Ito, M., Marungsri, B., Matsuoka, R., Ito, S., Arakawa, K, "Contamination
flashover, performance of hydrophobic polymer insulators with different core diameters,"
2004.
[25] J. Mrazek, "Physics of natural and artificial lightning discharges”, ," 2002.
[26] K. F. Portella, Piazza F., Inone, P. C., Ribeiro S., Cabussú, M. S, Pedreira C. D., Da Silva Ch.
C, "Atmosferic pollution (coastal and industrial) effects on electric energy distribution
insulators in Salvador, Brazil, metropolitan region," 2008.
[27] S. C. L.A., "Alta tensión y sistemas de transmisión," 1989.
[28] C. Fierro, J. L., Tapia, O. O y Arbiol, A. M. A., "Comportamiento de la tensión de flameo en
60 Hz de diferentes tipos de aisladores a una altitud de 2600 m," 2001.
[29] IEEE, "Guide for the application of insulation coordination," 1999.
[30] J. D. Ramírez, Martínez, M. L., "Influencia de las sobretensiones por descargas
atmosféricas y comportamiento de los sistemas de puesta a tierra en líneas aéreas de
transmisión y distribución," 2004.
[31] T. E. McDermott, T.A. Short, J.G. Anderson, "Lightning Protection of Distribution Lines”,
IEEE Trans. on Power Delivery," 1994.
[32] CFE-LAPEM, "Curso de coordinación de aislamiento," 2001.
[33] C.F.E., "Estadísticas de falla de líneas de transmisión," 1993-2002.
[34] A. R. RAMÍREZ, "ANÁLISIS DEL COMPORTAMIENTO DE LÍNEAS DE TRANSMISIÓN DE 115 KV
ANTE DESCARGAS ATMOSFÉRICAS," 2008.