ESCUELA POLITÉCNICA NACIONAL

FACULTAD DE INGENIERÍA ELÉCTRICA Y

ELECTRÓNICA

MODELACIÓN, SIMULACIÓN Y MEDICIÓN DE LOS CAMPOS

ELECTROMAGNÉTICOS EN UNA LÍNEA DE EXTRA ALTO
VOLTAJE MEDIANTE EL MÉTODO DE ELEMENTOS FINITOS,
CASO DE APLICACIÓN A LÍNEA DE TRANSMISIÓN DE 500 kV
SAN RAFAEL EL INGA

TRABAJO DE TITULACIÓN PREVIO A LA OBTENCIÓN DEL TÍTULO DE

INGENIERA ELÉCTRICA E INGENIERO ELÉCTRICO

KRISSIA DALINA CABEZAS RUBIO

[email protected]

PEDRO DAVID JIMÉNEZ CORRALES

[email protected]

DIRECTOR: Ing. RAÚL ANTONIO CANELOS SALAZAR MSc.

[email protected]

Quito, diciembre de 2019

 AVAL

Certifico que el presente trabajo fue desarrollado por Krissia Dalina Cabezas Rubio y Pedro
David Jiménez Corrales, bajo mi supervisión.

Raúl Antonio Canelos Salazar

DIRECTOR DEL TRABAJO DE TITULACIÓN
 DECLARACIÓN DE AUTORÍA

Nosotros, Krissia Dalina Cabezas Rubio y Pedro David Jiménez Corrales, declaramos bajo
juramento que el trabajo aquí descrito es de nuestra autoría; que no ha sido previamente
presentada para ningún grado o calificación profesional; y, que hemos consultado las
referencias bibliográficas que se incluyen en este documento.

A través de la presente declaración cedemos nuestros derechos de propiedad intelectual

correspondientes a este trabajo, a la Escuela Politécnica Nacional, según lo establecido
por la Ley de Propiedad Intelectual, por su Reglamento y por la normatividad institucional
vigente.

KRISSIA DALINA CABEZAS RUBIO PEDRO DAVID JIMÉNEZ CORRALES

 DEDICATORIA

La culminación de este trabajo de titulación la dedico a las personas más

importantes en mi vida, quienes han sido mi motor, mi fuerza y mi inspiración diaria.
Los que con su amor, comprensión, paciencia, dedicación, apoyo incondicional y
motivación constante me han alentado a no rendirme jamás. Los amo

A mi familia,

Nelson, Conchita, Bibita, Pedro Isaac.

A mis hermanos,

Leito y Danielito

A mi colega y compañero de vida

Pedro David

A mi tío

Adrián

A mis abuelitos

Rosa y Hugo

Krissia Dalina
DEDICATORIA

A mis padres, Pedro y Nancy

A mi hermana Alejandra

A Krissia

Pedro David
 AGRADECIMIENTO

A mi Dios y Virgen de la Caridad porque sin ellos nada de esto hubiera sido posible ya que
me han ayudado a cumplir este maravilloso sueño, depositando en mí fuerza y sabiduría
para no rendirme ante los obstáculos que se me presentaban, acompañándome y
protegiéndome cada día en mi carrera universitaria.

A mi familia por ser un pilar fundamental en mi vida, mi motivación diaria que, pese a la
distancia, siempre han estado y estarán en mi mente y mi corazón. Gracias papi Nelson,
mami Conchi, ñaña Bibi, Pedro Isaac y Tío Adrián por estar siempre conmigo, entenderme,
confiar y creer en mí, brindarme su apoyo incondicional en cada una de las decisiones que
he tomado, sus consejos, su amor, su trabajo y su esfuerzo permanente para que yo pueda
llegar a cumplir esta meta tan anhelada que no es más que un peldaño de los muchos que
faltan. Gracias a ustedes soy lo que soy.

A mis hermanos Leito y Danielito, que con su inocencia y ocurrencias siempre me han
sacado una sonrisa y en los momentos que más estresada estaba han sido mi mejor
terapia. Gracias a ustedes mis chiquitos que me motivan a superarme cada día para que
vean en su hermana ese ejemplo a seguir y cuando crezcan se sientan orgullosos de mí.

A mi colega y compañero de vida Pedro David por la confianza que depositó en mí para la
realización de nuestro trabajo de titulación, gracias por ser muy bueno conmigo, por su
ayuda incondicional, su paciencia, su amor. Gracias por todas las aventuras y experiencias
que hemos vivido para llegar hasta donde estamos, el camino no fue fácil, pero juntos
apoyándonos mutuamente lo logramos. Sé que Dios nos tiene preparado grandes cosas y
también sé que las seguiremos experimentando y disfrutando juntos. Gracias a su familia
por acogerme en su hogar, por el cariño y todas las palabras y consejos de aliento que me
han brindado.

Al Ing. Juan Ramírez por brindarnos el apoyo necesario para la realización y revisión de
este proyecto, por su interés, dedicación y tiempo. Gracias por estar siempre presto a
atender y responder cualquier inquietud que se presentó.

Al Ing. Raúl Canelos por aceptar ser nuestro director, mi gratitud y respeto por su paciencia,
su ayuda y su profesionalismo.
A CELEC EP Transelectric por abrirnos las puertas y acogernos dentro de sus
instalaciones. Así como al Ing. Fredy Coyago, Ing. Cristhian Pazmiño y a todo el personal
de la Institución por compartirnos la información necesaria, conocimientos adicionales,
experiencias vividas y enseñanzas las cuales aportan y enriquecen nuestra vida
profesional.

Finalmente, agradezco a la Escuela Politécnica Nacional y a la Carrera de Ingeniería

Eléctrica por haber permitido formarme en ella y a su vez enseñarme que con dedicación,
esfuerzo y constancia todo es posible. A mis maestros que en cada clase aportaron con
conocimientos sólidos, útiles en la vida profesional, a mis amigos y colegas que durante
todos estos años conocí y fueron testigos de la entrega que hoy se ve reflejada con la
culminación de mi carrera universitaria.

Eternamente agradecida.

Krissia Dalina
 AGRADECIMIENTO

A mis padres, Pedro y Nancy, por su constante apoyo, sus consejos y regaños, su esfuerzo
y sacrificio por formar una linda familia. A mi padre, un ejemplo de vida y a quien admiro
mucho, sus frases de aliento y sus anécdotas que me ayudaron a crecer como persona,
simplemente gracias. A mi madre, quien supo brindarme su apoyo incondicional y siempre
estar ahí para cuando yo la necesitaba.

A mi hermana Alejandra, con quien compartí la mayor parte de mi vida, entre risas, llantos,
peleas, y que siempre estuvo ahí en las buenas y en las malas para aconsejarme.

A Krissia, quien se fue convirtiendo en un pilar fundamental de lo que soy ahora, por estar
a mi lado en todo momento y nunca dejarme a la deriva, por ir construyendo un sueño que
ahora se volvió realidad, por su amor incondicional.

A Transelectric por la apertura brindada. De manera especial a Fredy Coyago y Cristhian

Pazmiño quienes con su tiempo y conocimiento supieron orientarnos en el desarrollo de la
presente Tesis. De la misma manera a los compañeros de Gestión Social y Ambiental y
Diseño de Líneas de Transmisión por su amistad y constante apoyo.

A Juan Ramírez por su predisposición y ayuda para que la presente Tesis salga adelante.

A nuestro director Raúl Canelos, quien aceptó ser parte de este trabajo de titulación y por
guiarnos en la realización del mismo.

A la Escuela Politécnica Nacional, por habernos formado académicamente y realizarnos

íntegramente como profesionales.

Pedro David
 ÍNDICE DE CONTENIDO

AVAL ....................................................................................................................... I
DECLARACIÓN DE AUTORÍA............................................................................... II
DEDICATORIA ...................................................................................................... III
DEDICATORIA ...................................................................................................... IV
AGRADECIMIENTO ............................................................................................... V
AGRADECIMIENTO ............................................................................................. VII
ÍNDICE DE CONTENIDO .................................................................................... VIII
RESUMEN ........................................................................................................... XII
ABSTRACT ......................................................................................................... XIII
1. INTRODUCCIÓN ............................................................................................. 1
1.1 Objetivos ................................................................................................... 3
1.2 Alcance ..................................................................................................... 4
1.3 Marco Teórico ........................................................................................... 4
2. METODOLOGÍA ............................................................................................ 42
2.1 Elementos que conforman la línea de transmisión de 500 kV San Rafael - El
Inga ................................................................................................................... 42

2.2 Simulación de campos eléctricos y magnéticos en el software FEMM ....... 54

2.3 Procedimiento de medición de campos eléctricos y magnéticos en el campo
.......................................................................................................................... 57

2.4 Instrumentos de medición ........................................................................... 59

3. RESULTADOS Y DISCUSIÓN ...................................................................... 63
3.1 Medición de campos electromagnéticos bajo la línea de transmisión de 500
kV San Rafael El Inga .................................................................................... 63
3.2 Validación de las simulaciones realizadas .................................................. 71
 3.3 Casos de estudio ........................................................................................ 84

3.4 Pérdidas por efecto corona ....................................................................... 130

3.5 Configuración óptima en el software PLS-CADD ...................................... 133
4. CONCLUSIONES ........................................................................................ 137
5. RECOMENDACIONES ................................................................................ 139
6. REFERENCIAS BIBLIOGRÁFICAS ............................................................ 140
7. ANEXOS ...................................................................................................... 144
ANEXO I ............................................................................................................. 145
ANEXO II ............................................................................................................ 148
 RESUMEN

En el presente estudio se modelan y simulan los campos electromagnéticos en la línea de

transmisión de 500 kV San Rafael El Inga. Mediante el programa computacional FEMM
basado en el Método de Elementos Finitos se resuelven las ecuaciones diferenciales de
Maxwell que están implícitas en los fenómenos electromagnéticos. Así, los valores
obtenidos son contrastados con las mediciones realizadas en el campo.

Se realizaron diseños de diferentes estructuras a escala 1:1 en una herramienta de dibujo

tipo CAD para luego ser exportados al programa FEMM que resuelve el problema y del que
se obtienen los resultados relacionados a los campos electromagnéticos de la línea de
transmisión de extra alto voltaje, considerando estructuras ubicadas a diferentes altitudes
sobre el nivel del mar y observando cómo afectan las condiciones atmosféricas del medio
en la intensidad de dichos campos.

Se simulan diferentes alternativas de geometría de la línea de transmisión para obtener los

valores de campos tanto eléctricos como magnéticos con el fin de proponer una
configuración óptima o proponer modificaciones en la configuración actual que reduzcan
los valores de los campos electromagnéticos.

Por último, mediante el software de diseño de líneas de transmisión PLS-CADD y las

simulaciones realizadas en el programa FEMM se muestra el diseño óptimo resultante para
una línea de transmisión de extra alto voltaje.

PALABRAS CLAVE: método de elementos finitos, FEMM, PLS-CADD, campos

electromagnéticos, línea de transmisión, ecuaciones diferenciales de Maxwell, condiciones
atmosféricas.
 ABSTRACT

In the present study, electromagnetic fields are modeled and simulated through San Rafael
- El Inga 500 kV transmission line. Using software FEMM based on the Finite Element
Method the Maxwell differential equations that are implicit in electromagnetic phenomena
are solved. Thus, the values obtained are contrasted with the measurements made on the
site.

Designs of different structures are made at a 1:1 scale in a CAD type drawing tool to be
exported to the FEMM software that solves the problem and from which the results related
to the electromagnetic fields of the extra high voltage transmission line are obtained,
considering structures located at different altitudes above sea level and remarking how
atmospheric conditions affect the intensity of these fields.

Different alternatives of geometry of the transmission line are simulated in order to obtain
the values of both electric and magnetic fields in order to propose an optimal configuration
or suggest modifications in the current configuration that reduce the values of the
electromagnetic fields.

Finally, through the PLS-CADD transmission line design software and the simulations
carried out in FEMM program, the resulting optimal design for an extra-high voltage
transmission line is shown.

KEYWORDS: finite element method, FEMM, PLS-CADD, electromagnetic fields,

transmission line, Maxwell differential equations, atmospheric conditions.
 1. INTRODUCCIÓN
Las líneas de transmisión de alto voltaje transportan la energía eléctrica a grandes
distancias de una manera eficiente y segura desde los centros de generación hacia los
centros de consumo. En los últimos años con los avances tecnológicos, la incorporación
de grandes cargas al sistema y el cambio de la matriz productiva y energética, el Ecuador
ha tenido un crecimiento tendencial del consumo de energía eléctrica, siendo este un hito

Codo Sinclair que posee 8 unidades de generación y cuenta con una capacidad instalada
de 1500 MW. Para evacuar la generación de este proyecto se estableció como mejor
alternativa la implementación de un sistema de transmisión de 500 kV de extra alto voltaje,
que además de vincular Coca Codo Sinclair con Quito (subestación El Inga) permitirá
interconectarlo con el principal centro de carga del país, Guayaquil (subestación Chorrillos)
[1] [2].

El sistema de 500 kV está compuesto por dos líneas de transmisión de simple circuito, que
parten desde la Subestación del proyecto hidroeléctrico Coca Codo Sinclair que se
encuentra ubicada en el cantón El Chaco, provincia de Napo y finalizan en la Subestación
El Inga situado en el cantón Quito, provincia de Pichincha, con una longitud total de 126
km cada línea. Las rutas de las líneas atraviesan zonas de altitudes que varían entre 1200
metros sobre el nivel del mar (msnm) en el sector del proyecto de generación Coca Codo
Sinclair y 4100 msnm en el sector más alto conocido como el Paso de la Virgen localizado
a 10 km de la subestación El Inga. Cada una de las fases de la línea de transmisión está
constituida por un haz de 4 conductores de aleación de aluminio ACAR 1100 MCM 18/19
y dos cables de guarda, uno de acero y el otro OPGW. Los conductores se sostienen en
estructuras metálicas de celosía, auto-
horizontal, como se observa en la Figura 1.1. Las cadenas de aisladores en suspensión
[3].

Figura 1.1.
El transporte de la energía eléctrica ocasiona una alteración en el ambiente en la que se
originan ciertas perturbaciones debido a los campos eléctricos y magnéticos, esta
problemática actualmente causa preocupación sobre posibles afectaciones a la salud [4].
Por esta razón se han desarrollado varios estudios e investigaciones para analizar los
posibles efectos que tienen los campos electromagnéticos sobre las personas y el
ambiente. Actualmente no existen evidencias científicas de los posibles daños que los
campos electromagnéticos pueden ocasionar, sin embargo diferentes compañías han
preferido realizar técnicas de diseño que disminuyan la emisión de los campos
electromagnéticos en las líneas eléctricas [5]. En nuestro país, en la línea de extra alto
voltaje de 500 kV es evidente que la incidencia de los campos electromagnéticos será más
intensa, ya que, a mayor nivel de voltaje, mayor intensidad de campos electromagnéticos.
Estas consideraciones no son tomadas en cuenta en el diseño, ya que, al ser una línea de
transmisión nueva, se espera que los valores estén de acuerdo con los límites establecidos
en las normas nacionales e internacionales para campos eléctricos y magnéticos
provenientes de fuentes de baja frecuencia (60 Hz) [6] [7]. Esto con el fin de asegurar que
la exposición humana a los campos electromagnéticos no tenga efectos perjudiciales para
la salud.

La intensidad de los campos eléctricos y magnéticos se ve influenciada por las condiciones

atmosféricas (presión, temperatura, humedad, altura sobre el nivel del mar) y también por
las condiciones propias de la configuración (disposición de la geometría de los
conductores, distancia de separación entre subconductores, diámetro de cada conductor,
separación entre fases, etc.) [8].

En la actualidad, el desarrollo de programas que utilizan métodos finitos brinda la facilidad

de mostrar gráficamente las situaciones relacionadas con los problemas
electromagnéticos. Por ejemplo, el programa computacional FEMM (Finite Element Method
Magnetics) que en base a métodos finitos permite resolver las ecuaciones diferenciales de
Maxwell implícitas en temas de campos eléctricos y magnéticos. Es así como, para
diferentes configuraciones geométricas de las estructuras de líneas de transmisión, se
pueden obtener valores de campos electromagnéticos bastante aproximados a los reales.

En el Ecuador no se han realizado estudios de campos electromagnéticos en líneas de

extra alto voltaje, por lo que, al no contar con una referencia para la modelación de una
línea de transmisión considerando campos eléctricos y magnéticos, se desconoce si la
disposición geométrica utilizada en la actualidad es la adecuada en lo referente a la
intensidad de campos que éstas producen.
A través de las herramientas computacionales AutoCAD y FEMM, se modelan y se simulan
los campos electromagnéticos de la línea de transmisión de 500 kV San Rafael El Inga
comparando los resultados de las modelaciones y simulaciones mediante las mediciones
en campo utilizando equipos de medida proporcionados por CELEC EP TRANSELECTRIC.

Mediante las simulaciones para diferentes alternativas de geometría de torres de la línea

de transmisión se obtienen los valores de campos tanto eléctricos como magnéticos con el
fin de proponer una configuración óptima o modificaciones en la configuración actual que
reduzcan los valores de estos campos. Todo lo modelado se pone a prueba y se valida con
mediciones en campo, observando cómo afectan las diferentes condiciones atmosféricas
del medio en la intensidad de los campos electromagnéticos medidos y simulados. La
configuración óptima resultante se presenta en el software PLS-CADD como un aporte para
la implementación de futuros proyectos de líneas de transmisión de extra alto voltaje.

1.1 Objetivos
El objetivo general de este Estudio Técnico es:

Simular la distribución de los campos electromagnéticos a lo largo de la línea de

transmisión de 500 kV San Rafael - El Inga mediante la utilización de los métodos
de elementos finitos con el uso del software FEMM (Finite Element Method
Magnetics) y la contrastación con mediciones bajo la línea de transmisión.

Los objetivos específicos de este Estudio Técnico son:

Investigar fuentes bibliográficas referentes a la modelación, simulación y medición

de campos electromagnéticos en líneas de transmisión de extra alto voltaje.

Modelar las estructuras y los campos electromagnéticos de la línea de transmisión

de 500 kV San Rafael - El Inga, con los datos proporcionados por Transelectric con
el uso del software FEMM (Finite Element Method Magnetics).

Validar el modelo de simulación del software empleado, mediante las mediciones

de campos electromagnéticos realizadas a lo largo del tramo de la línea de
transmisión San Rafael El Inga.

Observar las variaciones en la intensidad de campos electromagnéticos que

ocasionan las diferentes distribuciones geométricas de la línea de transmisión y el
medio atmosférico que lo rodea mediante el modelo ya validado.
 Proponer una distribución óptima de la geometría de la línea de transmisión de 500
kV o posibles modificaciones en la configuración actual que reduzcan los valores
de campos electromagnéticos mediante el software PLS-CADD.

1.2 Alcance
En este proyecto se modelan y simulan los campos electromagnéticos de la línea de
transmisión de 500 kV San Rafael El Inga empleando el software FEMM (Finite Element
Method Magnetics) basado en el método de elementos finitos, usando dimensiones reales
de los planos de las líneas de transmisión. Estos valores obtenidos, se comparan y
contrastan con mediciones realizadas en el campo.

Se incluye una parte práctica en la cual se realizan mediciones puntuales de campos

electromagnéticos a lo largo de la línea de transmisión San Rafael El Inga, utilizando para
esto los instrumentos de medición que posee CELEC EP - TRANSELECTRIC, siendo estas
mediciones un respaldo para contrastar el modelo propuesto.

Una vez validado el modelo se simulan diferentes escenarios donde se varía la altitud,
condiciones atmosféricas y factores que afectan tanto a la permitividad eléctrica y
permeabilidad magnética del medio. De la misma manera se prueban diferentes
configuraciones en las torres de transmisión existentes en el Ecuador y otras existentes
para 500 kV y posibles modificaciones a la geometría actual de ser el caso.

Usando el software PLS-CADD que posee CELEC EP - TRANSELECTRIC se propone el

diseño de la distribución geométrica de la línea de transmisión óptima para este nivel de
voltaje, el cual puede reducir la intensidad de los campos eléctricos y magnéticos y ser un
aporte para futuros proyectos implementados en el Ecuador.

1.3 Marco Teórico

1.3.1 SISTEMA NACIONAL DE TRANSMISIÓN

El sector eléctrico ecuatoriano ha sido transformado y fortalecido en todas sus etapas,

constituyéndose así en uno de los grandes aciertos a través de la propuesta del
fortalecimiento de la matriz energética como eje de la productividad [9] enfocada en
aprovechar los recursos naturales como fuentes de generación, con miras al crecimiento
de este estratégico sector. La incorporación de nuevas centrales de generación como Coca
Codo Sinclair y Sopladora, viene acompañada de la construcción de líneas y subestaciones
de transmisión [1], las cuales permiten transportar la energía eléctrica a grandes distancias
de una manera eficiente y segura desde los centros de generación hacia los centros de
consumo de todo el país.

Actualmente, a noviembre de 2019, el sector eléctrico ecuatoriano cuenta con una

capacidad de generación instalada de 7153 MW y líneas de transmisión que operan en
niveles de voltaje de 500, 230 y 138 kV siendo 461 km de líneas de transmisión de 500 kV,
3188 km de líneas de 230 kV las mismas que se encuentran formando un anillo de 230 kV
entre las subestaciones Molino, Zhoray, Milagro, Dos Cerritos, Pascuales (Guayaquil),
Quevedo, Santo Domingo, Santa Rosa (Quito), Totoras (Ambato) y Riobamba, vinculando
de forma directa a los principales centros de generación con los grandes centros de
consumo del país. Finalmente, las líneas de transmisión de 138 kV suman 2155 km de
longitud y alimentan 15970 MVA de transformación [1] [10].

Además, existen líneas de interconexión con el sistema eléctrico colombiano a través de

un enlace de 230 kV desde el primero de marzo de 2003, mientras que las transacciones
de energía eléctrica con el sistema eléctrico peruano se realizan mediante transferencias
de bloques de carga de Ecuador a Perú, o viceversa, desde el año 2005, en el marco de la
integración regional de la Comunidad Andina de Naciones CAN [10].

El Sistema Nacional de Transmisión (SNT) cuenta con diferentes zonas operativas, como
se presenta en la Figura 1.2(a), así como en la Figura 1.2(b) se listan las centrales de
generación y subestaciones que conforman las diferentes zonas operativas [10].

(a)
 (b)
Figura 1.2. Zonas operativas del SNT (a) Distribución de las zonas operativas en
Ecuador; (b) Generadoras y Subestaciones que la conforman [12] [10]

1.3.2 SISTEMA DE TRANSMISIÓN DE 500 kV

El Sistema de Transmisión 500 kV, comprende el tendido de 461 kilómetros de líneas de

transmisión que enlazan a cuatro subestaciones, San Rafael, ubicada cerca de la central
Coca Codo Sinclair, El Inga en las cercanías de Quito, Tisaleo en el centro del país y
Chorrillos en las cercanías de Guayaquil, incorporando tecnología de punta y modernos
sistemas que son pioneros en el país. Este sistema se complementa con 287 kilómetros
de líneas de 230 mil voltios [10] [11]. La Figura 1.3 representa el sistema de transmisión de
500 kV y sus obras asociadas con la construcción del mismo.
 Figura 1.3. Sistema de 500 kV y obras asociadas [11]

A manera de resumen, se presenta la Tabla 1.1 en la que se incluye las características de

la línea de transmisión de 500 kV.

Tabla 1.1. Características técnicas de la línea de transmisión de 500kV del Ecuador

Línea de Longitud Descripción

Transmisión [km]
Coca Codo Sinclair- 2 líneas de simple circuito, haz de 4
San Rafael 7 conductores ACAR 1100 MCM por
fase.
2 líneas de simple circuito, haz de 4
San Rafael -El Inga 119 conductores ACAR 1100 MCM por
fase.
1 línea de simple circuito, haz de 3
El Inga- Tisaleo 135 conductores ACAR 1100 MCM por
fase.
1 línea de simple circuito, haz de 3
Tisaleo- Chorrillos 200 conductores ACAR 1100 MCM por
fase.

1.3.3 LÍNEA DE TRANSMISIÓN SAN RAFAEL EL INGA A 500 kV

1.3.3.1 Subestación San Rafael, 500/230 kV, 450 MVA

La Subestación San Rafael se encuentra ubicada en el cantón El Chaco, provincia de

Napo, a una altitud de 1800 msnm y a 7.8 km al noroccidente del proyecto de generación
Coca Codo Sinclair. Su capacidad de transformación es de 450 MVA.
1.3.3.2 Subestación El Inga, 500/230 kV, 1800 MVA

Es una de las obras más importantes que forma parte del Sistema de Transmisión de 500
kV. Fue construida en un terreno de 17 hectáreas, localizada en el sector de El Inga, al Sur
Oriente de Quito en la provincia de Pichincha.

La subestación el Inga cuenta con patios de 500, 230 y 138 kV ubicados a una altitud de
2831 msnm. El patio de 500 kV posee una capacidad de transformación de 1800 MVA, el
cual constituye el 15% de la capacidad total instalada en el Sistema Nacional de
Transmisión. La subestación El Inga, con sus tres patios, concentra una capacidad de
trasformación de 2100 MVA, característica que la convierte en la subestación más grande
del país [12].

1.3.3.3 Líneas de Transmisión Coca Codo Sinclair-San Rafael-El Inga a 500 kV

Este tramo se encuentra constituido por dos líneas de transmisión de 500 kV simple
circuito, Circuito A y Circuito B, las mismas que parten desde la Subestación del proyecto
hidroeléctrico Coca Codo Sinclair y llegan hacia la subestación San Rafael, posteriormente,
las líneas recorren por el cantón Quijos de la provincia de Napo, y finalizan en la
Subestación El Inga, como se presenta en la Figura 1.4. La longitud total aproximada es
de 126 km cada línea. Las rutas de las líneas atraviesan por zonas de altitud que varían
entre 1200 msnm en el sector del proyecto de generación Coca Codo Sinclair, 4100 msnm
en el sector más alto conocido como el Paso de la Virgen y, 2850 msnm a la llegada a la
Subestación El Inga [3].

El circuito A está constituido por 282 estructuras mientras que el circuito B está constituido
por 277 estructuras a lo largo del tramo [13].

Figura 1.4. Ruta de las L/T Coca Codo Sinclair- San Rafael El Inga 500 kV [13]
Cada fase de la línea de transmisión está constituida por un haz de 4 conductores de
aleación de aluminio ACAR 1100 MCM 18/19. Cada una de las líneas lleva dos cables de
guarda. Uno de ellos de tipo Alumoweld o de acero galvanizado extra-resistente y el otro
de fibra óptica (OPGW) de 48 hilos. El aislamiento de las líneas de transmisión está
diseñado para operar en las altitudes antes indicadas; los conductores se sostienen en
figuración
horizontal [3].

Las líneas de transmisión cuentan con tres transposiciones completas, con los respectivos
descargadores de línea (pararrayos) en cada sitio. Las cadenas de aisladores en
, dependiendo de las
pendientes laterales al eje de la línea, utilizando aisladores de porcelana [3].

La potencia generada en la central hidroeléctrica Coca Codo Sinclair es transportada por

las dos líneas de transmisión. El voltaje máximo de la línea en estado estable (fase-fase)
es de 525 kV, la corriente máxima por fase es de 1894 A y la capacidad de transporte
mínima de la configuración de un haz de 4 conductores es de 474 A cada fase [53].

1.3.4 CAMPOS ELÉCTRICOS Y MAGNÉTICOS

1.3.4.1 Campo Eléctrico

El concepto del campo eléctrico nació de la necesidad de explicar el fenómeno de las

fuerzas a distancia entre dos o más partículas por lo que se define al campo como la región
espacial donde se presenta una magnitud física y eléctrica. En consecuencia, el campo
eléctrico es la fuerza que se manifiesta por medio del rechazo o atracción de cuerpos
cargados; brinda una idea de la intensidad de fuerza eléctrica que sufriría una carga
ubicada en un sitio determinado [17].

Se puede decir también que el campo eléctrico es la interacción entre dos objetos cargados
eléctricamente, en otras palabras, es el efecto que produce una carga eléctrica en el
espacio sobre otra carga eléctrica ubicada en cualquier punto que a su vez experimenta
una fuerza, de la misma forma la segunda carga crea el mismo efecto que afecta a la
primera [17].

El origen de un campo eléctrico se da por la presencia de una carga eléctrica la cual

produce alteraciones en el espacio donde se encuentra ubicada. Esta carga eléctrica
modifica las propiedades físicas del espacio provocando así que el campo eléctrico se
produzca [47].
Matemáticamente el campo eléctrico se describe como una cantidad vectorial el cual está
presente en todo punto del espacio. Se lo puede definir como la fuerza que actúa sobre
una carga de prueba (positiva o negativa) por unidad de carga, ecuación 1.1 [17].

(1.1)

La unidad del campo eléctrico según el Sistema Internacional de Unidades está dada en
Newtons/Coulomb, o Voltio/metro, [17].

La fuerza sobre una carga positiva tiene la misma dirección del campo eléctrico mientras
que para una carga negativa la dirección es opuesta a la del campo [18].

1.3.4.1.1 Fuentes del campo eléctrico

El campo eléctrico se manifiesta debido a cargas puntuales aisladas, cargas no puntuales,

cargas distribuidas.

1.3.4.1.1.1 Campo eléctrico debido a carga puntual aislada

El campo eléctrico alrededor de una carga puntual aislada se calcula por medio de la
relación mostrada en la ecuación 1.2 [17]:

(1.2)

Donde es el vector unitario que indica la dirección de la recta que une cada carga q.

La Figura 1.5 presenta que la dirección del campo eléctrico en una carga puntual positiva
es hacia afuera mientras que la de la carga puntual negativa es hacia dentro.

Figura 1.5. Esquema que representa la dirección del campo debido a una carga puntual
positiva y una carga puntual negativa [47]
1.3.4.1.1.2 Campo eléctrico debido a carga no puntual

El campo eléctrico proveniente de una carga no puntual se calcula imaginando una

distribución de carga que se encuentra subdividida en muchos elementos de carga ,
(ecuación 1.3) considerando las diferentes distancias desde el punto del cual se desea
calcular el campo [17].

(1.3)

Donde:

: Distancia entre y la posición de interés.

: Vector unitario que indica la dirección de la recta que une cada carga

1.3.4.1.1.3 Campo eléctrico debido a un sistema de carga distribuidas

Al tener un sistema de cargas distribuidas en el espacio (Figura 1.6), el campo eléctrico se

calcula en base al principio de superposición1 [47].

(1.4)

Siendo así el campo eléctrico es el resultante de una suma vectorial de los campos
eléctricos de cada una de las cargas.

Figura 1.6. Campo eléctrico debido a un sistema de cargas puntuales [47]

1.3.4.1.2 Líneas del campo eléctrico

Las líneas del campo eléctrico no son objetos materiales, al contrario, son una
representación gráfica que describe cualitativamente al campo eléctrico. Relacionan al
campo eléctrico con una región en el espacio así [15]:

El vector del campo eléctrico es tangente a la línea del campo eléctrico en cada
punto. La dirección que tiene la línea se muestra por medio de la punta de la flecha,
es igual al vector del campo eléctrico. La dirección de la línea representa la fuerza
sobre la carga positiva que se encuentra en el campo.
El número de las líneas por unidad de área que atraviesan una superficie
perpendicular a dichas líneas es proporcional a la magnitud del campo eléctrico en
dicha región. Por lo que las líneas del campo estarán más juntas donde el campo
eléctrico es más intenso y estarán más separadas donde el campo eléctrico sea
más débil.

Figura 1.7. Líneas del campo eléctrico que atraviesan dos superficies [15]

Como se observa en la Figura 1.7 la densidad de las líneas de campo a través de la

superficie A es mayor a que la densidad de las líneas de campo de la superficie B debido
a que la magnitud del campo eléctrico en A es mayor que en B.

En caso de que las líneas en diferentes ubicaciones apunten a distintas direcciones se

tiene un campo no uniforme [15].

Las líneas de campo empiezan en una carga positiva y terminan en una carga negativa, en
caso de haber un exceso en cualquier carga, las líneas empezarán o terminarán en el
infinito. El número de líneas es proporcional a la magnitud del campo y a su vez dos líneas
de campo no se pueden cruzar entre ellas [15].

1.3.4.1.3 Flujo eléctrico

Es una propiedad de cualquier campo vectorial referida a una superficie hipotética la cual
puede ser abierta o cerrada [15].
El flujo eléctrico es una magnitud escalar definida por el producto escalar [15]:

(1.5)

Se utiliza el símbolo para representar una integral sobre una superficie cerrada.

Donde:

(phi mayúscula) : Flujo eléctrico.

: Campo eléctrico.

: Área superficial.

Como se indica en la ecuación 1.5, la ecuación del flujo eléctrico es una integral de
superficie cerrada; esta debe ser evaluada sobre la superficie en cuestión. La superficie
cerrada es aquella que divide el espacio en una región exterior y una región interior [15].

El valor del flujo eléctrico depende del campo como de la superficie. Es proporcional al
número de las líneas del campo eléctrico que penetran en una superficie [15]. Su unidad

de medida es newtons por metros al cuadrado entre coulomb

El flujo eléctrico de cualquier superficie cerrada depende de la carga que se encuentra al

interior de esta, en caso de que la superficie cerrada no rodee a ninguna carga el valor del
flujo es igual a cero [15].

Dado el caso de que la carga neta, que se encuentra al interior de la superficie, sea igual
a cero, el valor del flujo es igual a cero, ya que el flujo entrante anula al flujo saliente [15].

1.3.4.1.4 Ley de Gauss para el campo eléctrico

La Ley de Gauss es de mucha importancia para estudiar los campos eléctricos ya que
describe la relación que existe entre el flujo eléctrico neto a través de una superficie cerrada
y la carga encerrada por la superficie [15].

Al suponer una carga puntual positiva q que se encuentra ubicada en el centro de una
esfera de radio r, la magnitud del campo eléctrico en todos los puntos de una superficie es
igual a la mostrada en la ecuación 1.6 [15]:

(1.6)

Las líneas de campo se dirigen radialmente hacia afuera siendo así perpendiculares a la
superficie en todos sus puntos, como se presenta en la Figura 1.8. Por lo que, en cada
punto de la superficie el vector del campo eléctrico es paralelo al vector del área ,
donde representa un elemento de área local, el cual está rodeando al punto en la
superficie [15].

Figura 1.8. Superficie gaussiana esférica de radio r que rodea una carga puntual q [15]

Por lo tanto,

(1.7)

Y por la ecuación 1.5 se determina que la ecuación del flujo neto a través de la superficie
gaussiana es igual a la mostrada en la ecuación 1.8 [15]:

(1.8)

Donde E se encuentra fuera de la integral ya que es constante en la superficie.

Adicionalmente al considerar la superficie esférica se tiene que .

Así pues, el flujo neto a través de la superficie gaussiana es mostrada en la ecuación 1.9
[16]:

(1.9)

Donde , lo que permite escribir que:

(1.10)

La ecuación 1.10 establece que el flujo del campo eléctrico a través de una superficie
cerrada A es igual a la carga neta q contenida dentro de la superficie, dividida por la
constante de permitividad eléctrica [16].

La superficie cerrada que se emplea para determinar el flujo del campo eléctrico debido a
una cierta distribución de cargas se denomina superficie gaussiana.
Ahora bien, la ley de Gauss también puede ser representada de forma diferencial en la que
considera el valor de la densidad de carga contenida al interior de un volumen V
rodeado por una superficie cerrada [16].

Utilizando el primer teorema de la divergencia para transformar el flujo del campo eléctrico
en una integral de volumen, la ley de Gauss se escribe como la ecuación 1.11 [16]:

(1.11)

Como ambos lados de la igualdad tienen diferenciales volumétricas, se considera al

volumen arbitrario [16].

La forma diferencial de la ley de Gauss es [14] [16]:

(1.12)

Donde:

E : Campo eléctrico

: Operador matemático muy variable que se aplica a números normales,

corrientes, vectores, en este caso en operaciones vectoriales para el cálculo del campo
eléctrico [16].

En la ecuación 1.12 de la forma diferencial de la ley de Gauss

nacen y mueren las líneas de campo y cómo de intenso es el proceso de nacimiento o

muerte de líneas [14].

Al calcular la divergencia de E, pueden suceder tres cosas [14]:

1. Si , , significa que ninguna línea nace ni muere en ninguna parte, es

decir toda línea que entra en el entorno de un punto, sale otra vez de él, así mismo
toda línea que sale entró antes por ahí mismo, en otras palabras, al no haber carga
en el entorno del punto las líneas de campo que entran salen otra vez como si nada.
2. Si , si la divergencia es positiva, eso significa que en el entorno
más grande es el

intenso. En otras palabras, al ser divergencia positiva salen más líneas de las que
entran.
 3. Si , si la divergencia es negativa, eso significa que en el entorno

más intenso.

La letra griega rho , representa la densidad de carga eléctrica que a su vez es una medida
de cuanta carga eléctrica ya sea positiva o negativa se encuentra alrededor del punto. A
mayor más carga eléctrica se encuentra alrededor del punto [14].

La letra griega épsilon , es la constante eléctrica o permitividad eléctrica del vacío. Es

una constante física que según el Sistema Internacional de Unidades tiene el valor de

[14].

La constante eléctrica es la relación numérica entre la carga y la fuerza eléctrica, este valor
no varía jamás, es una constante universal [14].

En resumen, la ecuación 1.12 tiene solo dos factores variables que son la densidad de la
carga eléctrica y la divergencia del campo eléctrico.

Manipulando la ecuación 1.12 matemáticamente se obtienen las líneas del campo eléctrico
conociendo así donde nacen y mueren, observando cuando van cambiando en el espacio
y cuantas aparecen y desaparecen.

El campo decrece proporcionalmente al cuadrado de la distancia, el comportamiento de las

líneas de campo aplicando la ecuación en un protón y un electrón es mostrada en la Figura
1.9:

Figura 1.9. Líneas de un campo de un dipolo eléctrico [14]

La ecuación 1.12 es la ley de Gauss para el campo eléctrico, la misma que forma parte de
las ecuaciones de Maxwell.

Esta ecuación señala que las cargas eléctricas son los lugares donde nacen y mueren las
líneas del campo eléctrico, nacen en cargas positivas y mueren en cargas negativas.

1.3.4.2 Campo Magnético

El campo magnético no es más que la representación matemática que describe la cantidad

de fuerza ejercida en una carga en movimiento. Esta representación permite describir la
distribución de una fuerza magnética en el espacio dentro y fuera de algo magnético. Todas
las fuentes del campo magnético son bipolares por naturaleza, poseen un polo norte y un
polo sur (positivo y negativo) [19].

Al existir una fuente de energía, habrá a su alrededor un campo magnético, que describe
cómo actúan estas fuerzas magnéticas. El campo magnético es una magnitud vectorial que
posee módulo y dirección. Cada vector apunta la dirección en la que lo haría una brújula y
su magnitud depende de la fuerza magnética a la que se encuentra sometido [19].

La intensidad de campo magnético H, además se define como un vector de campo de la

densidad de flujo magnético B. Tanto H como B se relacionan mediante u, (ecuación 1.13),
que no es más que la permeabilidad magnética del aire [32].

(1.13)

1.3.4.2.1 Fuentes del campo magnético

Un campo magnético se puede generar por medio de una carga puntual o conjunto de
cargas en movimiento (corriente eléctrica) o por la cercanía de un imán [22].

1.3.4.2.1.1 Campo magnético proveniente de un imán

Un imán es un material capaz de producir un campo magnético exterior y atraer a

elementos o partículas ferromagnéticas. Estos campos son creados naturalmente por un
imán o por un metal magnetizado que son consecuencia del movimiento de los electrones
alrededor del núcleo atómico de dicho elemento y además sobre su propio eje, en un giro
conocido en la física como spin [19].

Como se observa en la Figura 1.10, el polo norte del imán tiende a alinearse con el campo
magnético, de modo que una aguja de brújula suspendida girará hasta alinearse con el
campo magnético [48].
 Figura 1.10. Campo magnético de un imán [48]

1.3.4.2.1.2 Campo magnético generado por una carga puntual en movimiento

Si una carga q se mueve con una velocidad , como se observa en la Figura 1.11, crea un
campo magnético en todo el espacio [21] [22].

Figura 1.11. Descripción de una carga en el espacio [21]

Este campo se calcula mediante la siguiente expresión de la ecuación 1.14 [21]:

(1.14)

Donde:

q : Carga creadora del campo, su unidad es Coulomb [C].

: Velocidad con la que se mueve la carga, su unidad es metro por segundo

r : Distancia desde el punto donde se encuentra la carga hasta el punto P donde se está
calculando el campo, su unidad es el metro [m]

: Vector unitario que va desde el punto donde se encuentra la carga hacia el punto donde
se calcula el campo.

: Permeabilidad magnética, la cual depende del medio en el que se encuentra la carga.

El módulo del vector de campo magnético se puede calcular por la expresión [22]:

(1.15)

Mientras que la dirección y sentido viene dado por la regla de la mano derecha [22]:

Figura 1.12. Ley de la mano derecha

En la Figura 1.12 se puede observar que cuando la carga q es positiva la orientación del
dedo pulgar es en el mismo sentido que el vector velocidad mientras que si la carga q es
negativa va en sentido contrario. Es decir, los dedos de la mano derecha apuntan en
dirección del campo cuando el pulgar apunte a la dirección de la corriente. El campo
magnético en la dirección del movimiento es nulo, ya que en este caso los vectores
son paralelos y su producto vectorial es cero [22].

1.3.4.2.1.3 Campo magnético proveniente por un conjunto de cargas en

movimiento (corriente eléctrica)

El movimiento de cargas eléctricas en un conductor genera un campo magnético alrededor

de este. Al incrementar la corriente que pasa por el conductor la magnitud del campo
magnético se incrementa y disminuye al aumentar la distancia con respecto al conductor.
Este campo, cuya unidad es el Tesla, puede ser representado y calculado mediante la Ley
de Biot-Savart relacionando el valor del campo magnético con las corrientes que lo
producen [21].

1.3.4.2.1.3.1 Ley de Biot-Savart

El físico danés Jean Biot dedujo una ecuación que permite calcular el campo magnético H
creado por un circuito recorrido por una corriente de intensidad I. Esta ley propone que el
diferencial de intensidad del campo magnético, , se produce por un elemento diferencial
de corriente . El campo varía inversamente proporcional con el cuadrado de la distancia,
es independiente del medio que le rodea y tiene la dirección dada por el producto cruz
de y [20]. Lo anterior se indica en la ecuación 1.16.

(1.16)

Figura 1.13. Filamento de un alambre infinito por el que circula una corriente [20]

En pocas palabras, esta ley establece que en cualquier punto P la magnitud de la intensidad
del campo magnético que produce el elemento diferencial es proporcional al producto de
la corriente, la magnitud del diferencial de longitud y el seno del ángulo formado entre el
filamento y la línea que lo conecta con el punto P en donde se busca el campo [16].
Gráficamente esta ley se encuentra indicada en la Figura 1.13.

Ahora bien, cada elemento que conforma el filamento de corriente contribuyen a

obteniéndose así una suma infinitesimal. Esta sumatoria conduce a la forma integral de la
ley de Biot-Savart, ecuación 1.17 [49]:

(1.17)

1.3.4.2.1.3.2 Ley Circuital de Ampère

La ley circuital de Ampère establece que la integral de línea de H sobre cualquier trayectoria
cerrada es exactamente igual a la corriente encerrada por dicha trayectoria. Se expresa en
la ecuación 1.18 [49]:

(1.18)

En 1826 Ampère publicó una ley matemática (ecuación 1.19) que explicó las experiencias
realizadas por Ørsted en las que postuló que las corrientes eléctricas son fuentes de
campos magnéticos [14]:

(1.19)
 es el rotacional del campo magnético que indica el modo en el que las líneas de
campo giran alrededor de cada punto del espacio, mientras que es la densidad de
corriente eléctrica. Si es muy grande en un punto dado, quiere decir que en ese punto se
encuentra concentrada una gran intensidad de corriente eléctrica; por otro lado, si
significa que allí no hay corriente alguna. Esta densidad de corriente indica la intensidad
de corriente que atraviesa cada metro cuadrado de superficie y su dirección estará dada
por la dirección de corriente que es la dirección que circulan los electrones. [15]

Como se observa en la ecuación 1.19, la dirección de la corriente no coincide con la del

campo magnético, sino con el eje de giro del rotacional.

1.3.4.2.1.4 Ley de Gauss para el campo magnético

(1.20)

La ecuación 1.20 describe el campo magnético a través de su divergencia. La divergencia

del campo magnético ( ) no hace más que indicar dónde nacen y mueren las líneas
de campo. En este caso, si la divergencia es positiva nacen más líneas de campo que las
que mueren mientras que si es negativa pasa lo contrario mueren más líneas de las que
nacen y si es cero no nacen ni mueren líneas de campo [14].

Sin embargo, como se presenta en la ecuación 1.20, esta ley expresa que no existen
cargas magnéticas (monopolos magnéticos) y que las líneas del campo magnético no
nacen ni mueren en ninguna parte de manera completa, es decir, las líneas del campo
magnético no tienen principio ni fin. Otras características que muestra esta ecuación es
que el número de líneas del campo que entran en una región siempre son iguales al número
que sale y que las líneas del campo son siempre cerradas como se observa en la Figura
1.14.

Figura 1.14. Líneas del campo magnético de un imán [14]

Ahora bien, para explicar de una manera clara y dinámica la divergencia de B, se considera
lo siguiente: se tiene un dipolo eléctrico y un dipolo magnético. Si se separara el dipolo
eléctrico se tendría un monopolo eléctrico, sin embargo, si se intentase separar el dipolo
magnético no se podría obtener un monopolo magnético. Por ejemplo, se tiene un imán
con polo norte y polo sur, si se corta ese imán de tal manera que se quiera separar estos
polos, se podría observar que cada mitad de imán quedaría con un polo norte y un polo
sur. Dicho de otras palabras, La existencia de una carga positiva no exige la de una carga
negativa, pero la existencia de un polo norte sí exige la de un polo sur, por lo tanto, la
divergencia es nula [15].

1.3.5 LAS ECUACIONES DE MAXWELL

Las ecuaciones de Maxwell son la descripción del campo electromagnético. El campo

eléctrico produce un campo magnético cambiante y viceversa el campo magnético produce
un campo eléctrico cambiante. Cuando los campos son variables con el tiempo el campo
magnético H no puede existir sin el campo eléctrico E ni E puede existir sin H [14].

Con estas ecuaciones se puede saber el comportamiento del campo electromagnético en

una región determinada.

Las dos primeras ecuaciones son similares cuando no tienen dependencia del tiempo.
Estas fueron descritas en los capítulos Ley de Gauss para el campo eléctrico (ecuación
1.12) y Ley de Gauss para el campo magnético (ecuación 1.20).

1.3.5.1 Ley de Faraday

Michael Faraday, al analizar las experiencias realizadas por Ørsted, manifestó su creencia
de que, si una corriente podía producir un campo magnético, entonces un campo
magnético debería ser capaz de producir una corriente [16].

Tras años de realizar experimentos similares, descubrió lo que a futuro le serviría a Maxwell
para plantear uno de sus postulados. Enrolló un alambre alrededor de un anillo de hierro y
lo conectó a una batería. En el otro lado del anillo enrolló un alambre conectado a un
galvanómetro como se aprecia en la Figura 1.15. La idea de Faraday fue que, si una
corriente eléctrica podía generar un campo magnético, tal vez un campo magnético podía
generar una corriente eléctrica.
 Figura 1.15. Experimento de Faraday [14]

Al cerrar el circuito, Faraday observó una ligera deflexión de la aguja del galvanómetro.
Situación similar ocurrió cuando se desconectaba la batería, pero ahora la deflexión era en
sentido contrario. Es decir, al momento de conectar y desconectar la batería aparecía una
corriente eléctrica en el segundo circuito. Pudo concluir que no era el hecho de estar
energizado o no el circuito para que aparezca una corriente, sino que era la transición del
encendido y apagado lo que causaba que la corriente aparezca, entonces, la corriente se
producía por la variación del campo magnético [14].

Se puede decir que la razón de cambio del flujo magnético con respecto al tiempo produce
una fuerza electromotriz (fem) en una trayectoria cerrada no necesariamente conductora.
Cuando un conductor se mueve a través de un campo magnético, cortando el flujo, se
induce un voltaje en el conductor. En forma similar, cuando el flujo cruza un conductor
estacionario, también se induce un voltaje. La ley de Faraday se acostumbra a expresar
como la ecuación 1.21 [49]:

(1.21)

Si la trayectoria es una superficie cerrada, la ecuación 1.21 se describiría como la ecuación

1.22 donde se interpreta el flujo que pasa a través de cualquiera de las N trayectorias
cerradas [49].

(1.22)

Ahora bien, expresando la ecuación 1.22 en su forma integral se tiene la ecuación 1.23
[49]:

(1.23)

Esta integral de línea da lugar a una diferencia de potencial con campos variantes en el
tiempo resultando una fuerza electromotriz o un voltaje. Así, se puede expresar esta
integral lineal en una integral de superficie de B, así [49]:

(1.24)
Si se considera una trayectoria estacionaria, el flujo magnético B es la única cantidad que
varía con el tiempo, pudiendo introducir la derivada parcial dentro de la integral [49]:

(1.25)

Aplicando el teorema de Stokes2 a la integral cerrada se tiene [49]:

(1.26)

Pudiendo expresarse estas integrales como diferenciales [49]:

(1.27)

Llegando así a una de las cuatro ecuaciones de Maxwell (ecuación 1.27)

La parte izquierda de la ecuación 1.27 es el rotacional del campo eléctrico, mientras que la
parte derecha es la razón de cambio del campo magnético con respecto al tiempo. Si esta
razón de cambio es cero, obviamente quiere decir que el campo magnético no cambia en
el tiempo, mientras que si es pequeño la variación es lenta y si es grande el cambio es muy

brusco. Esto es, es la rapidez de cambio en el campo magnético contrario a ese

cambio, representado por el signo negativo. De la misma manera, la ecuación 1.27 indica
que el campo magnético variable en el tiempo es capaz de producir un campo eléctrico de
la nada tal que su rotacional tenga sentido contrario al del cambio del campo magnético
[15].

Si la densidad de flujo magnético B no es una función del tiempo (ecuación 1.25) y

(ecuación 1.27) se reducen a la ecuación 1.28 [49]:

(1.28)

1.3.5.2 Ley Circuital de Ampère en función del tiempo

La ecuación 1.19 mostrada en el capítulo 1.3.4.2.1.3.2 fue descrita cuando se aplica a

campos magnéticos estables.

Si se aplica la divergencia de ambos lados de la ecuación 1.19, se tiene la ecuación 1.29

[49]:

(1.29)
Usando identidades del cálculo vectorial se tiene que la divergencia de un rotacional es
igual a cero [23] [49]:

(1.30)

Sin embargo, la ecuación 1.30 es inconsistente con la conservación de la carga3. Maxwell

propuso que la divergencia de la densidad de corriente es igual al negativo de la derivada
de la densidad de carga respecto del tiempo. Esto es, sólo podrá haber un flujo de corriente
si la cantidad de carga varía con el paso del tiempo, ya que esta disminuye o aumenta en
proporción a la carga que es usada para alimentar dicha corriente [24]. Esto se puede
observar en la ecuación 1.31 que no es más que la ecuación de continuidad [49]:

(1.31)

Maxwell corrigió la ley de Ampère al agregarle este último término a dicha ecuación. De la
primera ley de Maxwell y sustituyendo en la ecuación de continuidad, ecuación 1.31, obtuvo
[14]:

(1.32)

Para postular la ecuación 1.32, James Maxwell se percató de que, al igual que un campo
magnético variable produce un campo eléctrico de la nada, también sucede lo contrario: un
campo eléctrico variable produce un campo magnético [14].

Si , que es el desplazamiento eléctrico y es una generalización del campo eléctrico

en presencia de un dieléctrico, entonces la ecuación 1.32 queda [14]:

(1.33)

De manera resumida, en la Tabla 1.2. se presentan las ecuaciones de Maxwell en su forma

puntual:

Tabla 1.2. Ecuaciones de Maxwell

Nombre Forma diferencial
Ley de Gauss para el campo magnético
Ley de Gauss para el campo eléctrico
Ley de Faraday

Ley de Ampère - Maxwell

1.3.6 LÍNEAS DE TRANSMISIÓN Y CAMPOS ELECTROMAGNÉTICOS

El voltaje y corriente en los conductores de una línea de transmisión implica la presencia

de campos eléctricos y magnéticos alrededor del conductor a frecuencias bajas de 50 Hz
o 60 Hz (dependiendo de la región) [32].

El campo eléctrico puede variar dependiendo de las características de la línea, esto es,
número de conductores por fase, diámetro del conductor, número de circuitos, disposición
geométrica del conductor, entre otros. Por otro lado, la intensidad del campo magnético
dependiendo de la magnitud de corriente que pasa por los conductores [32].

1.3.6.1 Fuerza del campo eléctrico

El campo eléctrico es un vector de campo de la fuerza del campo eléctrico definido en el

espacio a lo largo de tres ejes perpendiculares. Para campos sinusoidales en estado
estable cada eje es un número complejo o fasor cuya magnitud es expresada por su valor
RMS y su dirección expresada en coordenadas rectangulares como se indica en la
ecuación 1.34, y viene dada por [30] [31]:

(1.34)

Una representación general alternativa de la ecuación 1.34 puede ser expresada como un
vector que se mueve en el espacio describiendo una elipse como en la Figura 1.16(a), cuyo
semieje principal representa la magnitud y dirección del valor máximo del campo eléctrico,
y el eje menor representa la magnitud y la dirección del campo un cuarto de ciclo después
cuando el campo es mínimo como en la Figura 1.16(b). El campo eléctrico en la dirección
perpendicular al plano de la elipse es cero [30].

(a)
 (b)

Figura 1.16. (a) Representación del campo eléctrico describiendo una elipse, (b)
variación del campo eléctrico con el tiempo

1.3.6.2 Campo magnético

El campo magnético describe una elipse tal como se muestra en el campo eléctrico. La
resultante del campo magnético bajo una línea de transmisión viene dada por la ecuación
1.35 [30]:

(1.35)

En la cual son los valores RMS del campo magnético de las 3 componentes
perpendiculares de la elipse y son los valores RMS de los semiejes mayor y
menor de la elipse que representan los campos magnéticos máximos y mínimos,
respectivamente [30].

1.3.7 MÉTODOS DE RESOLUCIÓN DE ECUACIONES

Existen varios métodos los cuales permiten resolver las ecuaciones diferenciales de
Maxwell. Por un lado, se tienen métodos analíticos que restringen su aplicación en la
simulación de geometría de superficies complejas, mientras que, por otro lado, los más
utilizados, los métodos numéricos que, de una manera aproximada obtiene soluciones a
problemas utilizando netamente cálculos aritméticos, lógicos, de funciones y demás [25].

Existen diferentes técnicas para calcular los campos eléctricos y magnéticos en una línea
de transmisión, ya sean de forma analítica como numérica.

Métodos analíticos [25]:

Método de coeficientes de potencial de Maxwell
Método de Markt y Mengele
Método de las imágenes sucesivas
 Métodos numéricos [25]:
Método de simulación de cargas
Método de elementos de borde
Método de los elementos finitos

Para el presente caso de estudio se utiliza el método numérico de elementos finitos.

1.3.7.1 Método de elementos finitos

El presente método encuentra soluciones aproximadas a problemas de valor límite, es decir

soluciona problemas cuando se conocen los valores de sus condiciones de frontera. Es
comúnmente utilizado cuando se analizan problemas continuos (geometrías cerradas) [27].

Al utilizar este método, el área en análisis es discretizada en pequeñas regiones finitas

conocidas como elementos, generalmente en triángulos o cuadriláteros conectados entre
sí por puntos de unión llamados nodos. Las incógnitas del problema dejan de ser funciones
matemáticas y pasan a ser el valor de estas funciones en los nodos. Mediante la
interpolación de los valores conocidos en los nodos se puede llegar a conocer un valor
aproximado del modelo analizado. A mayor número de elementos finitos obtenidos mayor
exactitud, pero aumenta la complejidad en el planteamiento y resolución del problema [27].

Un elemento finito viene definido por sus nodos y por su contorno, formado por líneas que
los unen. Uno de los principales conflictos al aplicar este método es la determinación de
los coeficientes. Por este motivo, el problema se define de forma matricial, de tal manera
que al ingresar a un ordenador este arroje un conjunto de ecuaciones lineales o no, con el
fin de que pueda ser resuelto [27].

Considere un dominio con dos dieléctricos diferentes con voltajes aplicados a dos fronteras
según la Figura 1.17(a). Cuando el dominio es discretizado se crean tantos elementos con
sus respectivas fronteras como en la Figura 1.17(b).

(a) (b)

Figura 1.17. (a) Problema dieléctrico; (b) Problema discretizado [28]

El potencial de cada dominio varía linealmente con cada elemento de tal manera que [28]:

(1.36)

Si se considera un elemento finito del dominio discretizado en la Figura 1.17(a), se observa:

Figura 1.18. Elemento finito con sus coordenadas [28]

Este elemento debe satisfacer la ecuación de primer orden (ecuación 1.36), obteniéndose
así un sistema de 3 ecuaciones con 3 incógnitas:

(1.37)

El sistema de ecuaciones (ecuación 1.37) tendrá una solución no trivial si se consideran

condiciones de borde consistentes [28].

1.3.7.1.1 Condiciones de Borde

Una condición de borde, de frontera o de contorno son variables conocidas. Para el caso
de la solución de problemas electromagnéticos, las condiciones de Dirichlet.

1.3.7.1.1.1 Condiciones de borde de primer orden

Las condiciones de primer orden o también conocidas como condiciones de Dirichlet son
las que imponen el valor de la propia función (y no de sus derivadas) en ciertas regiones
(frontera) del dominio de estudio [29]. Cada nodo asociado al elemento tiene un potencial
conocido. En este caso, la variable o función que se va a determinar es constante.

1.3.7.1.1.2 Condiciones de borde de segundo orden

Las condiciones de segundo orden también conocidas como condiciones de Neuman

especifican el valor de la derivada normal en la frontera [29]. Aporta al conjunto de
derivadas parciales una solución de una derivada tomada sobre la frontera del elemento
finito.
1.3.7.1.1.3 Condiciones de borde de Robin

Son aquellas en las que se especifica una combinación del valor de la función y de su
derivada normal en la frontera [29].

1.3.8 SOFTWARE FEMM (FINITE ELEMENT METHOD MAGNETICS)

El software FEMM (Finite Element Method Magnetics) es un conjunto de programas que

permite resolver problemas electromagnéticos de baja frecuencia en dos dimensiones
planares y dominios asimétricos [31]. Se caracteriza por ser un software libre, rápido,
preciso y fácil de usar.

El programa computacional FEMM mediante el método de elementos finitos soluciona

problemas de campos magnéticos, electrostáticos, flujo de corriente y flujo de calor. La
Figura 1.19 presenta la estructura interna del programa FEMM.

1.3.8.1 Estructura interna del programa

Figura 1.19. Estructura interna del software FEMM

1.3.8.1.1 Interfaz

FEMM cuenta con una interfaz gráfica donde se puede interactuar con el programa por
medio de las herramientas de dibujo, selección, visualización, asignación de propiedades,
contornos, gráficos, etc. [31]

Dicha interfaz consta de dos bloques específicos: el preprocesador en el que se construye

una región de superficie cerrada donde se establece la geometría del problema a resolver
además de definir las propiedades de los materiales y las condiciones de contorno a utilizar,
y el post procesador que visualiza y analiza los resultados [31].

La interfaz se divide en 6 secciones principales mostradas en la Tabla 1.3.

Tabla 1.3. Interfaces de FEMM

Magnético
Preproceso Electrostático
Térmico
Magnético
Postproceso Electrostático
Térmico

1.3.8.1.2 Malla

El mallado consiste en segmentar la región de solución con una gran cantidad de

triángulos. Esto lo realiza el programa triangle.exe, el cual es independiente, pero a su vez
viene incluido en el software FEMM [31].

La malla se ejecuta al dar clic en el ícono Run mesh generator, Figura 1.20

Figura 1.20. Botón Run mesh generator

Las preferencias de la malla se pueden cambiar mediante dos parámetros:

Min Angle, donde se puede cambiar el ángulo mínimo de los triángulos que forman
la malla, por lo general, a mayor ángulo mayor es la calidad de la malla, sobre todo
en puntos críticos donde existen esquinas o irregularidades. Se puede poner
valores que se encuentren entre 1° y 33,8°.
Mesh size, al seleccionar un bloque y presionar la tecla espacio se accede al cuadro
de diálogo Properties for selected block en el que se puede modificar el tamaño de
la malla con el valor que se desee [33].

1.3.8.1.3 Cálculo

Los algoritmos de cálculo o solvers cambiarán dependiendo el tipo de problema que se

tenga que resolver. En la Tabla 1.4 se puede apreciar los tipos de problemas que el
software FEMM resuelve, así como sus respectivos solvers.
 Tabla 1.4. Algoritmos de Cálculo
Tipo de problema Solver
Magnético fkern.exe
Electrostático belasolv.exe
Térmico hsolv.exe
Flujo de corriente csolv.exe

Cada solver toma un conjunto de archivos de datos que describen un problema y resuelve
las ecuaciones de Maxwell para obtener valores del campo deseado en el dominio de la
solución [31].

El solver se ejecuta cuando se presiona el ícono Run análisis, Figura 1.21.

Figura 1.21. Botón Run analysis.

El tiempo que tarda en dar una solución depende netamente del problema a resolver y el
número de elementos finitos obtenidos al generar la malla. Los tiempos de solución varían
de menos de un segundo a varias horas esto se debe al tamaño y complejidad del
problema.

1.3.8.2 Etapas de programación de FEMM (Finite Element Method Magnetics)

Dentro del análisis a realizarse en el software FEMM se consideran tres bloques

específicos los cuales componen el proceso de resolución del problema [31]:

Preprocesador: en esta etapa se traza una región de superficie cerrada basada en

datos geométricos que pertenecen a los objetos a ser analizados, así mismo se
introducen las características propias del mismo y del ambiente que lo rodea.
Solver: A partir de la creación de nodos y elementos finitos, se empieza a resolver
los sistemas de ecuaciones.
Post procesador: Visualiza y analiza los resultados.

1.3.8.3 Tipos de problemas

El software FEMM resuelve algunos casos particulares de las ecuaciones de Maxwell,

ecuaciones que describen cómo va a ser y cómo va a comportarse el campo
electromagnético en una región determinada [14]. Las magnitudes que se utilizan en el
software FEMM son descritas en la Tabla 1.5.
 Tabla 1.5. Magnitudes en FEMM
Símbolo Magnitud Unidad SI
E Campo eléctrico N/C, V/m
B Inducción magnética T
H Intensidad del campo magnético A/m
J Densidad de corriente A/
Densidad de carga C/m, C/ , C/
Permitividad eléctrica F/m
Permeabilidad magnética H/m, N/
Conductividad media S/m
K Conductividad térmica W/(m*K)
T Temperatura K
q Calor generado por unidad de volumen W/

Para el caso que concierne al presente estudio, se analizarán problemas netamente

electrostáticos y magnéticos.

1.3.8.3.1 Problemas electrostáticos

Al analizar los problemas electrostáticos se considera el comportamiento de la intensidad

del campo eléctrico E, y la densidad del flujo eléctrico D (desplazamiento eléctrico).

Las ecuaciones de Maxwell relacionan a E y D, partiendo de la ley de Gauss que establece

que el flujo total que sale de una superficie cerrada es igual a la carga neta contenida dentro
de la superficie [20], se determinan las ecuaciones considerando campos estáticos,
ecuación 1.12 [31].

La densidad de flujo eléctrico D y el campo eléctrico E, se relaciona por medio de la

ecuación 1.38 [31]:

(1.38)

E y D tienen la misma forma solamente difieren por un factor que es una constante del
medio en el que se encuentran. El problema electrostático se resolverá siempre y cuando
la relación entre D y E sea lineal [31].

Para simplificar el cálculo de los campos, el programa emplea el potencial escalar eléctrico
V, ecuación 1.39, el cual se define por la relación existente con el campo eléctrico E [31].

(1.39)

Sustituyendo la ecuación 1.38 y ecuación 1.39 en ecuación 1.12 se tiene [31]:

(1.40)
Al ser el medio de la región de trabajo homogéneo, se retira la permitividad eléctrica de las
derivadas parciales involucradas en la divergencia, obteniéndose la ecuación 1.41 [31]:

(1.41)

Siendo la ecuación 1.41 la llamada Ecuación de Poisson4.

Cuando la región de trabajo contiene cargas en una distribución conocida , la ecuación

de Poisson puede usarse para determinar la función potencial [20]. Casi siempre la región
está libre de carga con una permitividad uniforme. Por lo que la ecuación de Poisson es:

(1.42)

Obteniéndose así la Ecuación de Laplace5.

El programa FEMM resuelve la ecuación 1.41, que no es más que el voltaje V sobre un
dominio con fuentes y condiciones de contorno definidos por el usuario.

1.3.8.3.2 Problemas magnetostáticos

Los problemas magnetostáticos son problemas con corrientes a frecuencia cero, donde las
corrientes inducidas se desprecian. Los campos no varían en el tiempo [31].

Las variables que se consideran para este tipo de problemas son la intensidad de campo
H y la densidad de flujo B las mismas que obedecen a la ecuación de Maxwell para campos
estáticos, mostrados en la ecuación 1.19 [31].

A su vez, al ser las líneas del flujo magnético curvas cerradas, no tienen ni principio ni fin
por lo que los campos B no tienen una fuente, ecuación 1.20 [31].

La intensidad del campo magnético H depende solo de las cargas en movimiento y es

independiente del medio, siendo así el campo de fuerza asociado con H es la densidad de
flujo magnético B, ecuación 1.13 [31].

Al tener una relación no lineal entre B-H, la permeabilidad está en función de B mostrada
en la ecuación 1.43. Para este efecto, FEMM puede resolver problemas no lineales [31].
 (1.43)

Para simplificar el cálculo de los campos, el programa emplea el potencial escalar

magnético A, el cual se define por la relación existente con densidad de flujo B, ecuación
1.44 [31].

(1.44)

Para una relación no lineal la ecuación 1.19 se reescribe para obtener la ecuación 1.45
[31]:

(1.45)

Para una relación lineal considerando que la ecuación 1.45 se reduce a la

ecuación 1.46 [31]:

(1.46)

FEMM trabaja con la ecuación 1.45 ya que el problema magnetostático si puede resolver
problemas no lineales que se presenten al relacionar B-H, a diferencia del problema
electrostático [31].

1.3.8.4 Características electromagnéticos usadas en el software FEMM

1.3.8.4.1 Permitividad eléctrica

La permitividad eléctrica de cualquier material describe cómo un campo eléctrico afecta y

es afectado por un medio. Se define como indica la ecuación 1.47 [49]:

(1.47)

Donde:

: permitividad en el vacío y es igual a 8.8541878176x10-12

: permitividad relativa del material

La permitividad relativa del material depende de diferentes factores como la temperatura,

la densidad, la humedad y composición [46].
1.3.8.4.2 Permeabilidad magnética

La permeabilidad magnética es un parámetro físico que describe cómo un material afecta

y es afectado por un campo magnético. Se define como indica la ecuación 1.48 [49]:

(1.48)

Donde:

-7

: permeabilidad relativa del material

1.3.8.4.3 Conductividad eléctrica

La conductividad eléctrica es la capacidad que tiene el material para permitir el flujo de

electrones. Su unidad es [siemens/metro]. Se define mediante la ecuación 1.49 [50]:

(1.49)

Donde:

La Tabla 1.6 presenta las diferentes permitividades eléctricas, permeabilidades magnéticas

y conductividades eléctricas de los materiales que componen la estructura de la línea de
transmisión.

Tabla 1.6. Características electromagnéticas [49] [51]

Permitividad Permeabilidad
Conductividad
Objeto Materiales eléctrica magnética
eléctrica [MS/m]
(relativa) (relativa)
Aislador Acero galvanizado 1000000 100 1,43
tipo Cemento 18 10 10-100
porcelana Porcelana 5,9-6,5 10 1,00E-19
Estructura
Acero galvanizado 1000000 100 1,43
Torre
Acero galvanizado 1000000 100 1,43
Aleación de aluminio 0 1-1,000032 40,8
Herrajes
Hierro fundido 1500-2000-
0 11,2
(comercial) 5000-200000
Aire Aire 1,0006 1,00000037 0
Moderadamente seco 3 12
Seco 15 12
Suelos <2E-7
Moderadamente
30 12
húmedo
1.3.9 FACTORES QUE INFLUYEN EN LA INTENSIDAD DE CAMPOS
ELECTROMAGNÉTICOS

La intensidad del campo eléctrico (gradiente potencial) alrededor del conductor puede
superar la rigidez dieléctrica del aire. Esto hace que el aire que lo rodea se ionice
volviéndose conductor (esto se puede representar como un radio ficticio del conductor,
Figura 1.22), provocándose así la conocida descarga corona [8] [52].

Figura 1.22. Conductor con efecto corona [52]

Existen diferentes factores que afectan el gradiente potencial y varían de la siguiente

manera [39]:

Humedad en el ambiente: a mayor humedad en el ambiente incrementa el gradiente

potencial y por ende el efecto corona.
Densidad relativa del aire: a menor densidad relativa del aire incrementa la
formación de efecto corona (campo eléctrico aumenta).
Presión atmosférica: si aumenta la presión atmosférica la densidad relativa del aire
también aumenta, disminuyendo el campo eléctrico.
Temperatura del ambiente: a mayor temperatura la densidad relativa del aire
disminuye, aumentando la formación del efecto corona.

Esto es, el campo eléctrico aumenta con cualquiera de los siguientes factores:

La humedad en el ambiente aumenta;

Temperatura del ambiente aumenta;
Densidad relativa del aire disminuye;
Presión atmosférica disminuye;

Los factores mostrados anteriormente además varían según la altitud sobre el nivel del
mar, esto es, a mayor altitud el gradiente potencial aumenta, y viceversa [42].
La forma de caracterizar los campos eléctricos y magnéticos se lo hace por la permitividad
y permeabilidad del aire, respectivamente. Así, en el presente estudio se analiza si
efectivamente, al cambiar la altitud sobre el nivel del mar varía el campo eléctrico y
magnético.

1.3.9.1 Gradiente Crítico Disruptivo

Este fenómeno juega un papel muy importante en el valor de la intensidad del campo
eléctrico en la superficie del conductor previo a la formación de efecto corona. Bajo
condiciones normales de temperatura y presión (25 °C y 760 mmHg, respectivamente) el
gradiente crítico disruptivo del aire es 29.8 kVpico/cm para causar la ionización del aire [42].

La fórmula de Peek6 es válida en condiciones normales de laboratorio. Si la densidad del

aire varía por efecto de temperatura o presión, esto afecta la rigidez dieléctrica del aire y
por consiguiente el gradiente crítico. La densidad del aire aumenta con la presión y
disminuye con la temperatura. Así mismo, la densidad del aire disminuye con la altitud. Es
decir, que a mayor altitud menor densidad del aire, mayor efecto corona, en otras palabras
[8]:

Altitud Densidad del aire

Altitud Efecto corona

La fórmula de Peek establece el voltaje crítico disruptivo, la misma que se describe en la

ecuación 1.50 [43]:

(1.50)

(1.51)

(1.52)

Donde:
r : radio del subconductor en cm
s : espaciamiento entre fases en cm
T : temperatura del ambiente en °C
P : presión atmosférica en cmHg
La ecuación 1.51 representa el gradiente crítico disruptivo del aire, mientras que la
ecuación 1.52 representa la densidad relativa del aire.

1.3.10 EFECTO CORONA Y PÉRDIDAS

El efecto corona es una manifestación de pérdidas en la línea de transmisión, por lo que

debe ser considerado en el diseño de una línea de transmisión de extra alto voltaje [8].

Este fenómeno se produce cuando el campo eléctrico en la superficie de los conductores

sobrepasa la rigidez dieléctrica del aire, produciendo emisión de energía acústica y energía
electromagnética en el rango de las radiofrecuencias, de forma que los conductores
pueden generar ruido e interferencias en la televisión y el radio. A una condición climática
lluviosa, el efecto corona es intenso, sin embargo, cuando la lluvia cesa y el conductor se
encuentre húmedo, el efecto corona es mucho más intenso [8].

La humedad relativa del aire, lluvia o niebla, hacen que las pérdidas por efecto corona
varíen de una manera proporcional, si aumentan estos factores, aumentan las pérdidas por
efecto corona [8].

Las pérdidas por efecto corona se expresan en kW/km-fase y resulta un poco difícil tener
un modelo analítico que permita tener valores exactos, ya que se involucran muchas
variables, sin embargo, para el cálculo de pérdidas por efecto corona se utilizan métodos
estadísticos y relaciones empíricas que se diferencian esencialmente por el tipo de
conductor de la línea y por las condiciones atmosféricas [8].

Para el cálculo de gradiente máximo superficial en líneas con conductores en haz se

considera un
se calcula el radio del círculo del conductor equivalente R [8]. Como se presenta en la
Figura 1.23.

(a) (b)

Figura 1.23. Radio del conductor único equivalente (a) 3 subconductores por haz, (b) 4
subconductores por haz
El radio del conductor único equivalente se calcula con la siguiente ecuación 1.53 [8]:

(1.53)

Donde:

: separación de los subconductores en cm

: radio del círculo de los subconductores en cm

: número de subconductores

El conductor único equivalente es un conductor ficticio por lo que tiene un radio ficticio
conocido como Radio Medio Geométrico. El RMG es calculado por la ecuación 1.54 [8]:

(1.54)

Donde:

: Radio Medio Geométrico (Radio ficticio) cm

: radio del subconductor en cm

La pérdida por efecto corona , en kW/km-fase, es calculada con la ecuación 1.55 [8]:

(1.55)

Donde:

: densidad relativa del aire

f : frecuencia del sistema, en Hz.

r : radio del subconductor, en cm.

DMG : distancia media geométrica entre fases, en cm.

V : voltaje efectivo fase-fase, en kV.

: Voltaje crítico disruptivo, en kV

El Voltaje de inicio de corona se calcula con la ecuación 1.56 [8]:

(1.56)

Donde:
 : factor de superficial del conductor

: factor ambiental

: es el radio del subconductor, en cm

: factor de disposición de haces.

El factor de disposición de haces se calcula con la ecuación 1.57 [8]:

(1.57)

Donde:

: radio del círculo de los subconductores en cm

N : número de subconductores por haz

La distancia media geométrica de un sistema trifásico se la calcula por medio de la

ecuación 1.58 [8]:

(1.58)

Donde son las distancias entre fases (1, 2, 3).

 2. METODOLOGÍA
En este capítulo se detallan las características de los elementos de las estructuras a
modelar de la línea de transmisión, se describe el procedimiento a seguir para realizar las
diferentes simulaciones de los campos eléctricos y magnéticos empleando el software
FEMM y a su vez se presenta el método aplicado que, según las normas se debe seguir
para medir los campos electromagnéticos a lo largo de la línea de transmisión. Las
mediciones obtenidas en campo serán comparadas con los valores obtenidos en la
simulación con el fin de validar el modelo y tener valores cercanos a los reales. En base al
modelo validado se observará si al cambiar diferentes parámetros como condiciones
atmosféricas, radio del conductor o geometría de los conductores, los valores de campos
eléctricos y magnéticos varían, obteniéndose así una configuración óptima que reduce
dichos valores.

2.1 Elementos que conforman la línea de transmisión de 500 kV

San Rafael - El Inga
La línea de transmisión está constituida básicamente por estructuras (torres), aisladores,
conductores eléctricos, cables de guarda, herrajes.

La línea de transmisión de 500 kV atraviesa diferentes altitudes sobre el nivel del mar
(Figura 2.1) y diversas condiciones atmosféricas.

(a)
 (b)

Figura 2.1. Perfiles de elevación de la ruta de la línea de transmisión de 500 kV San

Rafael El Inga. (a) Circuito A, (b) Circuito (B)

El tramo en estudio cuenta con tres diferentes tipos de estructuras como torre de
suspensión, retención y transposición.

2.1.1 ESTRUCTURAS (TORRES)

2.1.1.1 Torre de suspensión

Es una estructura en la que los conductores se encuentran suspendidos de las ménsulas7

a través de las cadenas de aisladores. Es diseñada para soportar el peso de los
conductores, cadenas de aisladores y herrajes además de la acción que ejerce el viento
sobre ellos. Las tensiones longitudinales son despreciables, son utilizadas en los tramos
rectos de la línea de transmisión. Se caracterizan por ser muy livianas, generalmente son
de celosía [36]. La Figura 2.2 presenta una estructura de suspensión del Circuito B ubicado
en el cantón Quito.
 Figura 2.2. Estructura de suspensión

2.1.1.2 Torre de retención

Soporta las mismas cargas de la torre de suspensión además que este tipo de estructura
también tiene cargas transversales que se producen por el ángulo entre los cables de
llegada y de salida.

Son usadas en los tramos que requieren un cambio de dirección en la línea de transmisión
como son vías, montañas o poblados. Se caracterizan por ser más robustas y tener una
apertura mayor entre las extensiones con el fin de soportar el momento del vuelco generado
por los cables [36]. La Figura 2.3 presenta una estructura de retención del Circuito B
ubicado en el cantón El Chaco.

Figura 2.3. Estructura de retención

Este tipo de torres se clasifican en dos clases:

Terminal: Se encuentra en el inicio o final de una línea, soporta la tensión que

ejercen los conductores ubicados de manera perpendicular a las ménsulas.
Angular: Se usa cuando hay cambio de dirección en la línea y se encarga de
soportar la tensión de los conductores producido por el cambio de dirección.

2.1.1.3 Torre de transposición

Es una estructura especial diseñada en aquellos puntos de la línea que por razones
técnicas de carácter eléctrico sea necesaria su instalación. Al viajar la energía eléctrica
distancias considerablemente largas se produce una autoinducción8 entre las fases, lo que
ocasiona un desequilibrio eléctrico en el circuito, por lo que es necesario cambiar la
disposición de los conductores. En este caso es necesario este tipo de estructura. La Figura
2.4 presenta una estructura de transposición del Circuito B ubicado en el cantón Quito.

Figura 2.4. Estructura de transposición

2.1.1.4 Estructuras seleccionadas para el caso de estudio

En base a la información de la tabla eléctrica AS BUILT proporcionada por CELEC EP-

TRANSELECTRIC, se escogen tres torres del circuito B considerando diferentes altitudes
sobre el nivel del mar, máxima, mínima y media. Determinándose para el estudio torres de
retención y una torre de transposición, las mismas que son de celosía de acero tipo cabeza
de gato.

Estas estructuras cuentan con dos numeraciones para ser distinguidas: número de
construcción y número AS BUILT. El número de construcción va incrementando
ascendentemente desde la central de generación Coca Codo Sinclair, GB1, hasta la
subestación El Inga, GB383N, mientras que la numeración AS BUILT inicia desde la
subestación San Rafael, 001, hasta la subestación El Inga, 258, de forma ascendente.

Cada torre cuenta con una nomenclatura normalizada propia que la caracteriza, por
ejemplo:

ST42(+12)-1

Donde:

El primer término indica que las estructuras que componen la línea de transmisión
es simple circuito.
El segundo término indica el tipo de la estructura, puede ser S para torres de
suspensión y T para torres de retención.
El tercer término indica el número de los conductores por haz, en este caso 4 indica
cuatro conductores por haz.
El cuarto término indica la zona en la que se encuentra la estructura. Existen dos
tipos de zonas, la zona 1 para altitudes entre 0-1000 msnm y zona 2 para altitudes
de 1000 msnm en adelante. La zona 2 a su vez está dividida entre zona 2.1 < 2000
msnm donde las estructuras se caracterizan por ser más livianas y la zona 2.2 <
4200 msnm donde las estructuras deben ser más robustas.
El quinto término indica la extensión del cuerpo de la torre.
El sexto término indica la extensión de la pata de la torre.

Las torres escogidas para analizar en este estudio, considerando las diferentes elevaciones
máxima, mínima y media, son las siguientes:
Altitud máxima: Torre de retención ST42(+12)-1, número de construcción GB353N,
número AS BUILT 229. Se encuentra ubicada a 4252.4 msnm en el cantón Quito,
provincia de Pichincha. Esta torre de retención se presenta en la Figura 2.5. La
ubicación en coordenadas UTM (Universal Transversal del Mercator) de esta
estructura es: Zona 17 M; Coordenada Este 808141.94 m E; Coordenada Norte
9962431.35 m S

Figura 2.5. Vista longitudinal de la estructura ST42(+12)-1

Altitud media: Torre de transposición ST42(-3)+1, número de construcción GB378N,

número AS BUILT 254. Se encuentra ubicada a 2976.4 msnm en el cantón Quito,
provincia de Pichincha. Esta torre de transposición se presenta en la Figura 2.6. La
ubicación en coordenadas UTM (Universal Transversal del Mercator) de esta
estructura es: Zona 17 M; Coordenada Este 797874.10 m E; Coordenada Norte
9965482.34 m S
 Figura 2.6. Vista longitudinal de la estructura ST42(-3)+1

Altitud mínima: Torre de retención ST41(+9)-0, número de construcción GB103,

número AS BUILT 56. Se encuentra ubicada a 1345.64 msnm en el cantón El
Chaco, provincia de Napo. Esta torre de retención se presenta en la Figura 2.7. La
ubicación en coordenadas UTM (Universal Transversal del Mercator) de esta
estructura es: Zona 18 M; Coordenada Este 199369.40 m E; Coordenada Norte
9977423.37 m S

Figura 2.7. Vista longitudinal de la estructura ST41(+9)-0

 De lo mencionado anteriormente, se resume la Tabla 2.1:

Tabla 2.1. Lista AS BUILT de las estructuras en la línea de transmisión

Torre Altitud
Torre No. metros APA Cross H Vano Vano
No. AS Torre sobre (m) ARM (m) adelante atrás
construcción BUILT nivel del (m) (m) (m)
mar
msnm
GB353N 229 ST42(+12) -1 4252.40 37.30 37 0.30 308.88 690.84
GB378N 254 ST42(-3) +1 2976.40 24.30 24 0.30 628.94 651.04
GB103 56 ST41(+9) +0 1345.64 35.30 35 0.30 400.16 1078.7

Donde:

APA: Altura al punto de amarre del conductor.

CROSS ARM: Distancia del mojón a la cruceta.

H: Altura del mojón

Vano: Distancia horizontal entre las torres que soportan las líneas de transmisión.

*Nota: Esta relación se cumple para una torre de retención.

2.1.2 AISLADORES

En las líneas de transmisión los aisladores o cadenas de aisladores sirven de soporte a los
conductores eléctricos, asegurando el aislamiento eléctrico entre conductor y estructura.

La cadena de aisladores se caracteriza por una codificación de 5 partes, por ejemplo,

T6A4R

Donde [54]:

El primer término indica el tipo de cadena, el mismo que:

I- V- JI- JV-jumper
T- ES-suspensión para el cable de guarda, ET-retención
para el cable de guarda.

El segundo término indica la carga que soporta el aislador y es:

Para conductor: 1-160 kN; 2-210 kN; 3-2x160 kN; 4-2x210 kN; 6-2x320 kN; default
para 70 kN.
Para el cable de guarda: default.

El tercer término indica la altitud y grado de contaminación:

Para el conductor: A-Altitud<2200 msnm y grado de contaminación I.

B-Altitud<2200 msnm y grado de contaminación II.
C-Altitud>2200 msnm y grado de contaminación I.
Para el cable de guarda: default.
Para el cable OPGW: 1-OPGW-90; 2-OPGW-136.

El cuarto término indica el tipo de cable:

Para conductor: 4-4xACAR 1100 MCM; 3-3xACAR 1100 MCM.

Para el cable de guarda: 3-Alumweld 7 No.8; 4-Galvanized Steel 3-8.

El quinto término indica un uso especial:

Para conductor: S- a; L-
larga; R-cadena de retención opuesta, G-cadena de retención para pórtico, D-
cadena de retención para torres sin salida en el lado del pórtico.
Para el cable de guarda: 1-conjunto de suspensión de abrazadera simple, 2-
conjunto de suspensión de abrazadera doble.

2.1.2.1 Aisladores en las estructuras seleccionadas

Las estructuras de retención escogidas en el presente estudio cuentan con:

Cadenas de aisladores de retención (backward-forward) compuestas por 25, 27 y

47 aisladores tipo U300BP. Las características de esta cadena de aisladores se
presentan en la Tabla 2.2.
Cadenas de aisladores de suspensión compuestas por 29 y 31 aisladores tipo
U70BLP en I y V respectivamente. Las características de esta cadena de aisladores
se presentan en la Tabla 2.3.
Los conductores al estar sujetos por las cadenas de aisladores de retención se encuentran
a lados opuestos de la torre por lo que deben ser conectados entre sí por medio de un
jumper9. Las cadenas de aisladores de suspensión en este tipo de torres son las que
sujetan al jumper interconectado eléctricamente.

Tabla 2.2. Cadena de aisladores de retención backward/forward

Torre Torre Cadena de No. Aisladores de

No. No. Retención Cadena de Tipo de
construcción AS (backward-forward) Retención Aislador
BUILT Herraje (backward-forward)
GB353N 229 T6C4/T6C4 2x27/2x27 U300BP
L: CB, T6C4T/T6C4
GB378N 254 M: BA, T6C4T/T6C4T U300BP
R: AC, T6C4/T6C4T
GB103 56 T6A4/T6A4R 2x25/2x25 U300BP

Tabla 2.3. Cadena de aisladores de suspensión

No.
Torre Torre Cadena de Aisladores
No. No. Suspensión para de Cadena de Tipo de
construcción AS Jumper Suspensión Aislador
BUILT (Ext/Mid/Int) para Jumper
Herraje (Ext/Mid/Int)
GB353N 229 J1C4/JVC4L/J1C4 31/2x31/31 U70BLP
GB378N 254 J1C4/JVC4L/J1C4 31/2x31/31 U70BLP
GB103 56 J1A4/JVA4S/J1A4 29/2x29/29 U70BLP

Las especificaciones de los dos tipos de aisladores que forman las cadenas en las
estructuras seleccionadas se presentan en la Tabla 2.4.

Tabla 2.4. Tipos de aisladores usados en el estudio

Carga
Tipo de Espaciamiento Diámetro Distancia soportada
Aislador (H) mm (D) mm de Fuga 10 por el
según mm aislador,
IEC 305 kN
U70BLP 146 255 450 70
U300BP 195 330 505 300
En la Figura 2.8 y Figura 2.9 se puede apreciar las dimensiones de fábrica de los tipos de
aisladores que fueron caracterizados en la Tabla 2.4.

Figura 2.8. Aislador U70BLP

Figura 2.9. Aislador U300BP

Los aisladores en las líneas de transmisión aéreas pueden ser de vidrio o de porcelana,
siendo los de porcelana utilizados en el caso de estudio ya que presentan mejores
resultados de aislamiento al estar en la intemperie [38]. En la Figura 2.10 se observan los
materiales internos que tiene el aislador de porcelana.
 Figura 2.10. Materiales de un aislador de porcelana

2.1.3 HERRAJES

Son estructuras metálicas que acompañan la estructura de soporte de las líneas de

transmisión. Se utilizan los diferentes planos de herrajes preestablecidos para cada tipo de
estructura y aislador en cada caso de estudio respectivamente. Son de acero
galvanizado11.

2.1.4 CONDUCTORES ELÉCTRICOS

El conductor de fase que se utiliza en la línea de transmisión de 500 kV San Rafael - El

Inga es el conductor 1100 MCM 18/19 ACAR (Aluminum Conductor Alloy Reinforced) cuyas
especificaciones técnicas se presentan en la siguiente Tabla 2.5.

Tabla 2.5. Capacidad de transporte del conductor

Diámetro D Peso W Resistencia R Capacidad

Conductor (mm) (kg/m) mínima de
transporte (A)
ACAR 1100 MCM 18/19 30.65 1.529 0.0557 474

Además, se cuenta con 2 hilos de guarda

diámetro de alambres de acero y otro de fibra óptica denominado OPGW, el mismo que
contiene hilos de fibra óptica en su interior y recubiertos por cable de acero galvanizado.
Sus características se detallan en la Tabla 2.6.

Tabla 2.6. Especificaciones técnicas de los cables de guarda

Cable de guarda Diámetro D Peso W Resistencia R

(mm) (kg/m)
9.14 0.404 -
OPGW-48B4-90 13.4 0.509 0.0557

2.2 Simulación de campos eléctricos y magnéticos en el software

FEMM

Para la simulación de campos eléctricos y magnéticos en el software FEMM se considera

lo siguiente:

La simulación no es variable en el tiempo

El problema se lo realiza en dos dimensiones en coordenadas cartesianas.
Los conductores tienen longitud infinita
La tierra es plana y es un perfecto conductor.

Ahora bien, considerando las características detalladas anteriormente se dibujan las 3

torres seleccionadas con sus respectivas estructuras, aisladores, herrajes y conductores a
escala real en AutoCAD en formato .dxf, para luego ser importadas al preprocesador de
FEMM.

En la Figura 2.11 se observan los tipos de problemas que FEMM puede resolver; para el
presente estudio se escoge el problema electrostático y problema magnético.

Figura 2.11. Selección de problema en software FEMM

Una vez escogido el problema se procede a seleccionar la opción Import DXF, donde
aparece un cuadro de diálogo para introducir el valor de tolerancia, por lo general el valor
predeterminado es suficiente. Dependiendo del archivo .dxf FEMM toma cierto tiempo en
importarlo. En la Figura 2.12 se presenta el archivo CAD importado en FEMM.
 Figura 2.12. Estructura importada en el software FEMM

Al simular los campos eléctricos y magnéticos producidos en una línea de transmisión, es

importante conocer los valores de diferentes parámetros que pueden afectar la intensidad
de estos campos.

Se colocan los parámetros de los materiales que se van a utilizar. En el caso del problema
electrostático se requiere tener el valor de la permitividad mientras que para el problema
magnético se requiere tener el valor de la permeabilidad. Las Figuras 2.13 y Figura 2.14
muestran el proceso de asignación de materiales.

Figura 2.13. Selección de Materiales y Contornos

Figura 2.14 Lista de materiales a ser modificados

Adicionalmente se definen las condiciones de borde de los conductores según los voltajes
de fase más altos del sistema, para ello se considera el voltaje pico que corresponde al
voltaje máximo de operación de la línea por fase en estado estable, y la condición de borde
de voltaje igual a cero en el suelo.

Teniendo la información necesaria para la simulación se realiza la discretización del

dominio, donde se generan los triángulos de los elementos finitos en el software como se
presenta en la Figura 2.15.

Figura 2.15. Torre de transmisión mallada en FEMM

Finalmente, se obtienen los resultados proporcionados por el programa tanto del campo
eléctrico, campo magnético, así como voltaje, como se presenta en la Figura 2.16. Los
valores obtenidos pueden ser analizados mediante las gráficas que facilita el software
FEMM.

Figura 2.16. Resultados de la simulación de la torre de transmisión en FEMM

2.3 Procedimiento de medición de campos eléctricos y magnéticos
en el campo

Para poder medir tanto los campos eléctricos como magnéticos, es necesario seguir
procedimientos de tal manera que los valores medidos no sean erróneos. Hay normas
internacionales que citan métodos adecuados para medir estos campos como es la norma
IEEE Std. 644-1994 Standard Procedures for Measurement of Power Frequency Electric
and Magnetic Fields From AC Power Lines.

2.3.1 MEDICIÓN DE CAMPOS ELÉCTRICOS

Según la norma IEEE Std. 644-1994, la medición de intensidad de campos eléctricos se

debe realizar a una altura de 1 metro sobre el nivel del suelo considerando el mismo
procedimiento para la medición del campo magnético [30].

2.3.1.1 Perfil longitudinal

La medición debe realizarse a 1 metro sobre el nivel del suelo paralelo a la línea de
transmisión tomando al menos cinco mediciones equidistantes desde la mitad del tramo de
la línea de transmisión avanzando en ambas direcciones de tal manera que se cubra toda
la longitud del vano [6] [30], como se presenta en la Figura 2.17.

Figura 2.17. Medición lateral (transversal) y longitudinal. Vista superior del vano [30]

2.3.1.2 Perfil lateral

El perfil lateral debe medirse una vez obtenidos los puntos más altos de intensidad del
campo eléctrico medidos en el perfil longitudinal. Generalmente el campo eléctrico es
mayor en la mitad del vano. Se realiza la medición del campo eléctrico a 1 metro de altura
sobre el nivel del suelo en dirección perpendicular a la línea. Las mediciones deben ser
realizadas desde el centro de la línea del área de interés y recorrer al menos 30 metros
lateralmente desde la ubicación del último conductor tomando al menos cinco mediciones
equidistantes. Es recomendable que se realicen las mediciones lateralmente hasta que se
cubra la longitud de la franja de servidumbre12. La Figura 2.18 ilustra la medición del campo
eléctrico lateralmente desde la mitad del vano de una línea de transmisión [6] [30].

Figura 2.18. Medición perfil lateral [30]

2.3.2 MEDICIÓN DE CAMPOS MAGNÉTICOS

La medición de los campos magnéticos sigue el mismo procedimiento que los realizados
para medir la intensidad del campo eléctrico.

2.3.2.2 Perfil longitudinal

La medición deberá realizarse a 1 metro sobre el nivel del suelo paralelo a la línea de
transmisión tomando al menos cinco mediciones equidistantes desde la mitad del tramo de
la línea de transmisión avanzando en ambas direcciones de tal manera que se cubra toda
la longitud del vano [6] [30].

2.3.2.1 Perfil lateral

El perfil lateral debe medirse una vez obtenidos los puntos más altos de intensidad del
campo magnético medidos en el perfil longitudinal. Generalmente el campo magnético es
mayor en la mitad del vano. Se realiza la medición del campo magnético a 1 metro de altura
sobre el nivel del suelo en dirección perpendicular a la línea. Las mediciones deben ser
realizadas desde el centro de la línea del área de interés y recorrer al menos 30 metros
lateralmente desde la ubicación del último conductor tomando al menos cinco mediciones
equidistantes. Es recomendable que se realicen las mediciones lateralmente hasta que se
cubra la longitud de la franja de servidumbre. La Figura 2.18 ilustra la medición del campo
magnético lateralmente desde la mitad del vano de una línea de transmisión [6] [30].

Para el presente estudio se realizan las mediciones de campos eléctricos y magnéticos a

1 metro sobre el nivel del suelo de tal manera que se pueda validar la simulación. Así
mismo, las mediciones se las realiza vano adelante y vano atrás, es decir, de las tres
estructuras definidas para el estudio se hace el barrido longitudinal con dirección a la
siguiente torre y la que se encuentra posterior a la estructura bajo análisis, así como un
barrido transversal de tal manera que se cubre los 30 metros de franja de servidumbre.

2.4 Instrumentos de medición

Existen diferentes instrumentos que permiten medir campos electromagnéticos y campos

eléctricos y magnéticos por separado a diferentes frecuencias, en este caso, se consideran
instrumentos para medición de campos electromagnéticos en líneas de transmisión a 60
Hz.

Para el presente estudio se utilizan los instrumentos de medida de campos eléctricos y

magnéticos SPECTRAN NF-5035 y Survey Meter HI-3604 los cuales son medidores de
cuerpo libre (Free Body Meter13) y de uno, dos o tres ejes, según la opción seleccionada.
Las características técnicas de cada uno de los equipos se detallan en la Tabla 2.7 y Tabla
2.8, respectivamente, y los equipos se presentan en la Figura 2.19 y Figura 2.20.

Tabla 2.7. Especificaciones equipo de medición SPECTRAN NF-5035 [34]

SPECTRAN NF 5035
Specifications base unit
Frecuency Range (min) 1 Hz
Frecuency Range (max) 30 MHz
Electric Field [V/m] (min) (typical) 0.1 V/m
 Electric Field [V/m] (max) (typical) 5000 V/m
Magnetic Field [Tesla] (min) (typical) 1 pT
Magnetic Field [Tesla] (max) (typical) 500 uT
Magnetic Field [Gauss] (min) (typical) 10 nG
Magnetic Field [Gauss] (max) (typical) 5G
Analog input [V] (min) typical 200 nV
Analog input [V] (max) typical 2V
RBW (resolution bandwidth) min 0.3 Hz
RBW (resolution bandwidth) max 1 MHz
Demodulator AM/FM
Units (additional unitis via PC software V, V/m, T, G, A/m
Detector RMS/MinMax
Internal Datalogger (size) 64K
FFT resolution (points) 1024
Lowest Sample Time 10 mS
Accuracy (typical) 3%

Figura 2.19. SPECTRAN NF-5035

Tabla 2.8. Especificaciones equipo de medición Survey Meter HI-3604 [35]

MODEL HI-3604 Survey Meter

Electric Specifications
Frequency Range
30-2000 Hz
(nominal)
Magnetic Field: +0,5, -2,0 dB (30-1000 Hz) -2,0,
Frequency Response
-6,0 dB (1000-2000 Hz)
(typical)
Electric Field: +0,5,-2,5 dB (39-2000 Hz)
Electric: 1 V/m-200 kV/m
Dynamic Range
Magnetic: 0.2 mG -20 G
Response True RMS
Logging On-Board, 112 Readings (max)
Electric fields: 1 V/m-199 kV/m
Sensitivity
Magnetic field: 0.1 mG-20 G
 Figura 2.20. Survey Meter HI-3604

Por otro lado, se utiliza el equipo Hybrid-Powered Environmental Quality Meter 850068 que
permite tomar medidas de temperatura y humedad relativa del aire. Las características de
este equipo se presentan en la Tabla 2.9, y el equipo se observa en la Figura 2.21.

Tabla 2.9. Especificaciones equipo Hybrid-Powered Environmental Quality Meter 850068

[41]

Medición Unidad Rango Resolución Precisión

Temperatura °C 0 a 50 °C 0.1 °C ±1.2 °C
ambiente °F 32 a 122 °F 0.1 °F ±2.5 °F

<70% RH
Humedad ±4% RH
%RH 10 a 95% RH 0.1% RH
relativa
±(4% reading ± 1.2% RH)

m/s 0.4 a 25.0 m/s 0.1 m/s

±3% F.S.
Velocidad del
viento
km/h 1.4 a 90 km/h 0.1 km/h >20 m/s
±4% F.S.
 Figura 2.21. Hybrid-Powered Environmental Quality Meter 850068

Para conocer la presión atmosférica a diferentes altitudes, se emplea la fórmula de Halley,

esta es la ecuación 2.1 [8]:

(2.1)

Donde:

: presión atmosférica en [mmHg]

: altura sobre el nivel del mar [m]

Una vez que se han simulado los campos eléctricos y magnéticos en el software FEMM,
se contrastan estos valores con las mediciones realizadas a lo largo de la línea de
transmisión en las estructuras escogidas de tal manera que se puedan ajustar los
parámetros del modelo a fin de que los valores simulados se acerquen a la realidad. Se
manejan varios escenarios en los cuales se varía la altitud con respecto al nivel del mar,
observando si este aspecto influye en la intensidad de los campos electromagnéticos
emitidos por una línea de transmisión.

Conocida la presión atmosférica del aire, empleando la ecuación 2.1, así como la
temperatura ambiente, se calcula el gradiente crítico disruptivo del aire con el fin de
establecer el voltaje de inicio de corona y el gradiente superficial del conductor, ya que un
menor valor del campo eléctrico en la superficie del conductor disminuye el efecto corona.
Tras determinar la configuración de línea de transmisión que mejor se ajuste en la
reducción de campos electromagnéticos, se calculan las pérdidas por efecto corona en
cada una de las estructuras analizadas y se determina si efectivamente, el nuevo diseño,
reduce las pérdidas por efecto corona comparado con la configuración actual.
 3. RESULTADOS Y DISCUSIÓN
En el presente capítulo se detallan los resultados obtenidos una vez medidos los campos
eléctricos y magnéticos en las diferentes estructuras seleccionadas que conforman la línea
de transmisión San Rafael El Inga. Con los datos obtenidos en la medición se validaron
las simulaciones realizadas en el software FEMM. Luego se realizaron simulaciones para
diferentes casos de estudio y diferentes escenarios. En base a esto, utilizando el software
PLS-CADD, se muestra la configuración óptima que reduce los valores de campos
electromagnéticos.

3.1 Medición de campos electromagnéticos bajo la línea de

transmisión de 500 kV San Rafael El Inga
La medición de los campos electromagnéticos se realizó de acuerdo con los
procedimientos presentados en las Normas Internacionales en la sección 2.3 con el fin de
obtener valores reales. Se mide el perfil longitudinal de la línea de transmisión tomando la
mayor cantidad de puntos equidistantes posibles, de esta manera se cubre tanto el vano
de adelante y el vano de atrás de la estructura escogida. Una vez obtenidos los valores
máximos del campo eléctrico y magnético longitudinalmente se mide el perfil lateral de la
línea de transmisión cubriendo el ancho de la franja de servidumbre que para el caso de la
línea de 500 mil voltios es 30 metros a cada lado. Así mismo para poder validar las
simulaciones realizadas se mide el perfil lateral tomando como referencia el centro de la
torre. La medición es realizada a 1 metro sobre el nivel del suelo.

Los instrumentos utilizados fueron: Spectran NF-5035, Survey Meter HI-3604 e Hybrid-
Powered Environmental Quality Meter 850068 descritos anteriormente, los cuales se
aprecian en la Figura 2.19, Figura 2.20 y Figura 2.21.

Para realizar las mediciones se consideró el procedimiento recomendado por la IEEE, ya

que al medir el campo eléctrico el equipo de medición en este caso, Survey Meter HI-3604
debe estar a 1 metro sobre el nivel del suelo y alejado a cualquier objeto o persona a una
distancia de al menos 2.5 metros, ya que la cercanía a estos va a afectar en los valores
obtenidos, causando errores en las mediciones. Mientras que, para medir el campo
magnético, no existe problema alguno ya que el valor del campo no es sensible a las
variaciones de la altura o presencia de objetos y personas [32].

Los valores medidos deben estar dentro de los límites establecidos en las normas
nacionales e internacionales para campos eléctricos y magnéticos provenientes de fuentes
de baja frecuencia (60 Hz) [6] [7]. En la Tabla 3.1 se presentan los valores límites de
exposición para campo eléctrico y magnético a un metro sobre el nivel del suelo.

Tabla 3.1. Niveles de referencia para la exposición a campos eléctricos y magnéticos de

60 Hz [6]
Intensidad del Intensidad del Densidad de Flujo
Tipo de campo Eléctrico campo Magnético Magnético (B)
exposición (E) (V/m) (H) (A/m) (Microteslas)
Público en General 4167 67 83
Personal
Ocupacionalmente 8333 333 417
expuesto

Las mediciones realizadas fueron hechas en tres escenarios, los cuales cuentan con
diferentes condiciones atmosféricas, altitud sobre el nivel del mar, vegetación, tipo de suelo
y planimetría. Por lo que se presentan 3 fichas técnicas con información relevante del medio
en el que se encuentran las estructuras analizadas y observaciones detalladas de
presencia de objetos, agua, vegetación y maleza que pudiesen afectar la medición de los
valores de campo electromagnético. Las mismas se presentan en la Tabla 3.2, Tabla 3.3,
y Tabla 3.4.

El tiempo utilizado para poder cubrir con las mediciones de los 3 casos de estudio
analizados fue de 5 días, en este caso son días escogidos aleatoriamente y a su vez con
condiciones climáticas diferentes dependiendo del día.

Tabla 3.2. Ficha técnica para medición de campos electromagnéticos SR-IN 2_254
Tabla 3.3. Ficha técnica para medición de campos electromagnéticos SR-IN 2_229
Tabla 3.4. Ficha técnica para medición de campos electromagnéticos SR-IN 2_56
Para realizar las respectivas simulaciones se consideró el voltaje y corriente en ese
instante. CELEC EP - TRANSELECTRIC se encarga de brindar el Reporte Post-Operativo
Diario del Sistema Nacional de Transmisión. Para el caso de estudio San Rafael - El Inga,
circuito B, se conocen los valores de voltaje de barra y corriente para el instante en el que
se midieron los valores de campos electromagnéticos. Este reporte se presenta en la Tabla
3.5.

Tabla 3.5. Resumen Post-Operativo Días de Medición

Resumen Post-Operativo Diario

Día de Voltaje Potencia
Estructura
medición (kV) (MVA)
27/8/2019 515.257 405.275
254
28/8/2019 513.291 360.906
229 30/8/2019 490.399 428.391
16/9/2019 500.000 414.472
56
17/9/2019 500.000 419.568

Para la simulación del campo eléctrico se requiere el voltaje fase-neutro por lo que debe
ser calculado con el voltaje de barra que se tiene del Resumen Post-Operativo como se
indica en la ecuación 3.1.

Ejemplo de cálculo para el voltaje fase-neutro de la estructura 254, considerando el

día 27/08/2019:

(3.1)

Donde:

Tabla 3.6. Voltaje fase-neutro para simulación del campo eléctrico

Voltaje
Estructura fase-neutro
(kV)
254 420.705
229 400.409
56 408.248

En la Tabla 3.6 se indica el valor calculado del voltaje fase - neutro para cada una de las
estructuras a analizar.

Para la simulación del campo magnético se requiere la densidad de corriente por lo que
debe ser calculada con la potencia que se tiene del Resumen Post-Operativo.

Ejemplo de cálculo para la densidad de corriente de la estructura 56, considerando

el día 17/09/2019:
Partiendo de la ecuación del flujo de potencia de la ecuación (3.2):

(3.2)

(3.3)

Donde:

La ecuación (3.3) es la corriente del conductor que circula en ese instante. Partiendo de
dicha ecuación se tiene:

Cálculo de la corriente máxima de todo el haz de conductores con la ecuación (3.4)

(3.4)

Cálculo de la corriente máxima de un subconductor del haz con la ecuación (3.5)

(3.5)

Tabla 3.7. Corriente del conductor, corriente máxima y corriente de un subconductor del
haz para simulación del campo magnético

Corriente
Corriente del Corriente máxima
máxima de un
Estructura haz de del haz de
subconductor
conductores (A) conductores
en el haz (A)
254 454.114 642.215 160.554
229 504.347 713.255 178.314
56 484.476 685.152 171.288
La Tabla 3.7 presenta los valores obtenidos del cálculo de la corriente del haz de
conductores, corriente máxima del haz de conductores y la corriente de un subconductor
del haz.

Para el cálculo de la densidad de carga, se necesita el valor de la corriente máxima de un

subconductor del haz y el área del conductor. Para el presente caso de estudio el conductor
es ACAR 1100 MCM 18/19, cuyo diámetro es de 30.65 mm y radio de 15.325 mm.

Se calcula el área del conductor con la ecuación (3.6)

(3.6)

La ecuación (3.7) es para calcular la densidad de corriente, parámetro requerido para la

simulación del campo magnético.

(3.7)

Donde:

Tabla 3.8. Densidad de la corriente para simulación del campo magnético

Densidad
Estructura de corriente
(MA/ )
254 0.217604768
229 0.241675682
56 0.232153423
La Tabla 3.8 presenta los valores obtenidos del cálculo de la densidad de corriente cuya
unidad es mega Amperio por metro cuadrado. Parámetro requerido para realizar la
simulación del campo magnético.

3.2 Validación de las simulaciones realizadas

Mediante el software FEMM se realizó la simulación del campo eléctrico y magnético de
las estructuras medidas. Se simuló la vista longitudinal y transversal de las torres. Esto
permite conocer los valores de permitividades y permeabilidades de los materiales de la
estructura, conductores, aire y suelo que mejor se ajustan a las diferentes condiciones
atmosféricas presentes en las torres analizadas.

3.2.1 VISTA LONGITUDINAL CAMPO ELÉCTRICO

3.2.1.1 Estructura 254

Figura 3.1. Distribución del campo eléctrico en la estructura 254 Vista longitudinal
 3.2.1.2 Estructura 229

|
Figura 3.2. Distribución del campo eléctrico en la estructura 229 Vista longitudinal

3.2.1.3 Estructura 56

Figura 3.3. Distribución del campo eléctrico en la estructura 56 Vista longitudinal

La Figura 3.1, Figura 3.2 y Figura 3.3 presentan la distribución del campo eléctrico desde
el centro de la torre y se observa que, entre más alta es la torre, el campo eléctrico a nivel
del suelo disminuye. Esto se puede comparar entre la Figura 3.1 y Figura 3.2, los cuales
las alturas del punto de amarre (APA) son de 37 metros y 24 metros respectivamente.

No existe presencia del campo eléctrico en el suelo, ya que como se observa en la Figura
3.4 las líneas de flujo del campo eléctrico no penetran la superficie del suelo.

3.2.2 VISTA LONGITUDINAL CAMPO MAGNÉTICO

3.2.2.1 Estructura 254

Figura 3.4. Distribución del campo magnético en la estructura 254 Vista longitudinal
3.2.2.2 Estructura 229

Figura 3.5. Distribución del campo magnético en la estructura 229 Vista longitudinal

3.2.2.3 Estructura 56

Figura 3.6. Distribución del campo magnético en la estructura 56 Vista longitudinal

La Figura 3.4, Figura 3.5 y Figura 3.6 presentan la distribución del campo magnético en las
estructuras seleccionadas. Se observa que existe una mayor presencia del campo
magnético cerca del conductor que en la propia superficie del suelo. Todo esto relacionado
con la cantidad de corriente que circula por el conductor.
A diferencia del campo eléctrico, las líneas de flujo penetran fácilmente la superficie del
suelo.

3.2.3 VISTA TRANSVERSAL CAMPO ELÉCTRICO

3.2.3.1 Estructura 254

Figura 3.7. Distribución del campo eléctrico en la estructura 254 Vista transversal

3.2.3.2 Estructura 229

Figura 3.8. Distribución del campo eléctrico en la estructura 229 Vista transversal
3.2.3.3 Estructura 56

Figura 3.9. Distribución del campo eléctrico en la estructura 56 Vista transversal

La Figura 3.7, Figura 3.8 y Figura 3.9 presentan la distribución del campo eléctrico en las
estructuras seleccionadas. Esto es, se realizan las mediciones del campo eléctrico vano
adelante (con dirección hacia El Inga) y vano atrás (con dirección hacia San Rafael). A
medida que hay un acercamiento mayor entre el conductor y el suelo, se registran valores
mayores del campo eléctrico.

En el caso del campo eléctrico, ciertos tipos de suelos facilitan la conducción eléctrica.
Diferentes tipos de sales, metales y minerales pueden ayudar a que el suelo sea un mejor
conductor. Los tipos de suelo presentes en las estructuras medidas son diferentes. Se
tienen suelos secos, moderadamente húmedos e incluso lugares en los cuales existe
presencia de vertientes de agua o ríos circundantes. Al medir la diferencia de potencial
existente entre los conductores y el suelo, se obtiene que, a menor altitud sobre el nivel del
mar el campo eléctrico es menor, mientras que, si se aumenta la altitud, el campo eléctrico
es mayor.
3.2.4 COMPARACIÓN DE RESULTADOS

A fin de poder validar la simulación, se redujo la diferencia entre los valores simulados y
medidos obteniendo un máximo del 7% [55]. Una vez que no se haya sobrepasado este
error, se pueden tomar los parámetros de la simulación como fijos, a fin de ser utilizados
en futuras simulaciones. La Tabla 3.9, Tabla 3.10 y Tabla 3.11 presentan los datos del
campo eléctrico y magnético medidos y obtenidos mediante la simulación.

Tabla 3.9. Resultados del campo eléctrico medidos y simulados de la estructura 254

CAMPO ELÉCTRICO
Perfil Longitudinal
Vano adelante

Ubicación del E [kV/m] a 1 metro de

Distancia recorrida del altura
punto de %Error
vano [m]
medición
Medido Simulado
0 Dentro de la torre 0.1731 0.180358 4.19
1 5 3.7146 3.84 3.29
2 10 3.58535 3.79 5.66
3 15 2.92425 2.85 2.42
4 20 2.72626 2.80 2.58
5 25 2.67547 2.81 5.10
6 30 2.13762 2.09 2.06
7 35 1.98801 2.12 6.53
8 40 1.95445 2.06 5.59
9 45 1.35278 1.44 6.63
10 50 1.3506 1.41 4.11
11 60 0.93152 0.92 1.53
12 70 0.55821 0.57 2.44
13 80 0.49477 0.51 3.52
14 90 0.28832 0.277848 3.63
15 100 0.35263 0.373109 5.81
16 110 0.21973 0.21 5.77
17 120 0.19725 0.20 0.28
18 130 0.15942 0.17 5.45
19 140 0.10136 0.10 5.38
Vano atrás

Ubicación del E [kV/m] a 1 metro de

Distancia recorrida del altura
punto de %Error
vano [m]
medición
Medido Simulado
0 Dentro de la torre 0.172214 0.180358 4.73
1 5 2.68541 2.85 6.08
 2 10 3.75285 4.01 6.96
3 15 3.86284 4.02 4.17
4 20 3.87641 4.10 5.65
5 25 3.77412 4.02 6.44
6 30 2.65432 2.71 2.17
7 35 1.31263 1.27445 2.91
8 40 1.35642 1.28 5.98
9 45 1.37136 1.28613 6.21
10 50 1.32782 1.27 4.32
11 60 1.0221 0.96 6.50
12 70 0.43623 0.41 6.35
13 80 0.38928 0.41 6.56
14 90 0.36425 0.39 6.38
15 100 0.25478 0.24 4.96
16 110 0.23621 0.25 3.88
17 120 0.14838 0.16 4.96
18 130 0.11025 0.11 2.85
19 140 0.10124 0.11 6.78
20 150 0.09615 0.10 3.90
DENTRO DE LA TORRE
Lado izquierdo Campo Eléctrico

Ubicación del E [kV/m] a 1 metro de

Distancia recorrida del altura
punto de %Error
vano [m]
medición
Medido Simulado
1 10 4.69 4.439 5.35
2 20 4.8 4.891 1.90
3 30 2.79 2.848 2.08
Lado derecho Campo Eléctrico

Ubicación del E [kV/m] a 1 metro de

Distancia recorrida del altura
punto de %Error
vano [m]
medición
Medido Simulado
1 10 3.812 4.069 6.74
2 20 2.321 2.251 3.02
3 30 1.71 1.605 6.14
Lado izquierdo Campo Magnético

Ubicación del B [uT] a 1 metro de

Distancia recorrida del altura
punto de %Error
vano [m]
medición
Medido Simulado
1 10 9.5725 9.91 3.50
2 20 7.65125 7.94 3.78
3 30 6.655 6.82 2.42
 Lado derecho Campo Magnético

Ubicación del B [uT] a 1 metro de

Distancia recorrida del altura
punto de %Error
vano [m]
medición
Medido Simulado
1 10 11.12 11.80 6.12
2 20 7.15625 7.38 3.07
3 30 5.3425 5.63 5.45

Tabla 3.10. Resultados del campo eléctrico medidos y simulados de la estructura 229

CAMPO ELÉCTRICO
Perfil Longitudinal
Vano adelante

Ubicación del E [kV/m] a 1 metro de

Distancia recorrida del altura
punto de %Error
vano [m]
medición
Medido Simulado
0 Dentro de la torre 0.059 0.06 5.16
1 5 0.082 0.08 2.19
2 10 1.36 1.34 1.28
3 15 2.65 2.68 1.22
4 20 2.29 2.40 4.78
8 40 0.458 0.43 6.04
9 50 0.246 0.23 5.29
10 55 0.133 0.12 6.26
11 60 0.121 0.13 4.97
12 80 0.045 0.05 1.47
14 100 0.02856 0.03 2.19
15 110 0.05625 0.06 4.79
16 125 0.0788 0.08 5.94
Vano atrás

Ubicación del E [kV/m] a 1 metro de

Distancia recorrida del altura
punto de %Error
vano [m]
medición
Medido Simulado
0 Dentro de la torre 0.059 0.06 5.16
1 5 0.066 0.06 5.99
2 10 1.829 1.95 6.56
3 15 1.819 1.87 2.67
4 20 1.31 1.24 4.96
5 25 1.097 1.10 0.44
6 30 0.98 0.94 4.18
8 40 0.41 0.39 6.02
 10 50 0.369 0.34 6.55
11 60 0.28 0.29 5.32
12 70 0.19 0.18 4.81
13 80 0.089 0.09 6.36
14 90 0.079 0.08 5.65
15 95 0.047 0.05 5.73
DENTRO DE LA TORRE
Lado izquierdo Campo Eléctrico

Ubicación del E [kV/m] a 1 metro de

Distancia recorrida del altura
punto de %Error
vano [m]
medición
Medido Simulado
1 10 1.373 1.45 5.28
2 20 1.597 1.71 6.98
3 30 1.179 1.25 6.28
Lado derecho Campo Eléctrico

Ubicación del E [kV/m] a 1 metro de

Distancia recorrida del altura
punto de %Error
vano [m]
medición
Medido Simulado
1 10 1.101 1.03 6.54
2 20 3.594 3.78 5.22
3 30 1.671 1.56 6.70
Lado izquierdo Campo Magnético

Ubicación del B [uT] a 1 metro de

Distancia recorrida del altura
punto de %Error
vano [m]
medición
Medido Simulado
1 10 9.7225 10.40 6.98
2 20 4.51375 4.25 5.91
3 30 1.865 1.97 5.43
Lado derecho Campo Magnético

Ubicación del B [uT] a 1 metro de

Distancia recorrida del altura
punto de %Error
vano [m]
medición
Medido Simulado
1 10 12.70625 13.518 6.39
2 20 9.498 10.155775 6.93
3 30 4.0165 4.295125 6.94
Tabla 3.11. Resultados del campo eléctrico medidos y simulados de la estructura 56

CAMPO ELÉCTRICO
Perfil Longitudinal
Vano adelante

Ubicación del E [kV/m] a 1 metro de

Distancia recorrida del altura
punto de %Error
vano [m]
medición
Medido Simulado
0 Dentro de la torre 0.00000288 0.000003 4.17
1 5 0.00000209 0.000002 4.31
2 10 0.57825 0.618476 6.96
3 15 1.26358 1.339790 6.03
4 88 0.30385 0.287382 5.42
5 206 0.0462 0.043063 6.79
6 226 0.03374 0.031430 6.85
7 246 0.025215 0.024079 4.51
8 266 0.019234 0.017996 6.43
9 276 0.019187 0.017944 6.48
10 281 0.01426 0.015220 6.73
11 291 0.014945 0.015887 6.30
12 301 0.012742 0.013375 4.97
13 311 0.010263 0.010958 6.77
14 316 0.009413 0.010036 6.62
Vano atrás

Ubicación del E [kV/m] a 1 metro de

Distancia recorrida del altura
punto de %Error
vano [m]
medición
Medido Simulado
0 Dentro de la torre 0.00000288 0.000003 4.17
1 5 0.00000213 0.000002 6.10
2 101 0.1508 0.155021 2.80
3 105 0.144843 0.154935 6.97
4 233 0.1458 0.136143 6.62
5 243 0.133674 0.140251 4.92
6 253 0.14315 0.152293 6.39
7 263 0.15632 0.165217 5.69
8 273 0.169753 0.179874 5.96
9 283 0.187287 0.196345 4.84
10 293 0.186689 0.199599 6.92
11 303 0.191634 0.199653 4.18
12 313 0.21952 0.230840 5.16
13 323 0.26952 0.282061 4.65
14 333 0.30874 0.329555 6.74
 15 343 0.40127 0.386085 3.78
16 353 0.41074 0.388887 5.32
17 358 0.41517 0.388880 6.33
18 403 0.73295 0.782323 6.74
19 410 0.829 0.880083 6.16
20 422 1.11 1.047060 5.67
21 442 0.96714 1.033730 6.89
22 834 0.01098 0.010220 6.92
DENTRO DE LA TORRE
Lado izquierdo Campo Eléctrico

Ubicación del E [kV/m] a 1 metro de

Distancia recorrida del altura
punto de %Error
vano [m]
medición
Medido Simulado
1 5 0.02852 0.0277689 2.63
2 10 0.21853 0.231538 5.95
3 15 0.31418 0.302889 3.59
4 20 0.66536 0.624359 6.16
5 25 1.48052 1.51343 2.22
6 30 1.35812 1.38086 1.67
Lado derecho Campo Eléctrico

Ubicación del E [kV/m] a 1 metro de

Distancia recorrida del altura
punto de %Error
vano [m]
medición
Medido Simulado
1 5 0.003679 0.00387277 5.27
2 10 1.45236 1.37116 5.59
3 15 1.18982 1.27271 6.97
4 20 0.57179 0.611131 6.88
Lado izquierdo Campo Magnético

Ubicación del B [uT] a 1 metro de

Distancia recorrida del altura
punto de %Error
vano [m]
medición
Medido Simulado
1 5 0.595 0.5657875 4.91
2 10 9.9975 10.6965125 6.99
3 15 7.7475 8.2554 6.56
4 20 7.6825 8.21855 6.98
5 25 8.98125 8.519 5.15
6 30 7.83125 8.151525 4.09
 Lado derecho Campo Magnético

Ubicación del B [uT] a 1 metro de

Distancia recorrida del altura
punto de %Error
vano [m]
medición
Medido Simulado
1 5 0.67625 0.6420575 5.06
2 10 9.57 10.153525 6.10
3 15 7.905 8.3068 5.08
4 20 6.36125 6.787325 6.70

En base a la ecuación 1.50, 1.51 y 1.51 descrita en el Capítulo 1 se calculó el voltaje crítico
disruptivo del aire, así como el voltaje de inicio de corona y de esta manera se observó la
relación que existe entre la altitud sobre el nivel del mar y la formación de efecto corona.

Ejemplo de cálculo para el voltaje crítico disruptivo en los conductores

adelante/atrás de la estructura 229:

P = 445.55 mmHg = 44.555 cmHg

T = 1 °C

Con la ecuación 1.52:

Con el valor de la densidad reemplazado en la ecuación 1.51, se obtiene el valor del

gradiente disruptivo del aire:

Con el gradiente crítico disruptivo reemplazado en la ecuación 1.50, se puede calcular el

voltaje crítico disruptivo:

Todos los valores obtenidos anteriormente son valores pico.

El factor m es un factor de corrección para irregularidades en la superficie del conductor y
condiciones ambientales que se obtiene del producto entre el coeficiente de estado
ambiental (mt) y el coeficiente de superficie (mc). La Tabla 3.12 detalla los valores que
puede tomar dichos coeficientes [8]:

Tabla 3.12. Factores de superficie m [8]

Factor ambiental mt Factor superficial mc

1 Ambiente seco 1 Conductor liso ideal
0.8 Ambiente húmedo 0.96 - 0.88 Conductor nuevo y limpio
0.25 - 0.16 Ambiente lluvioso 0.87 - 0.80 Conductor áspero envejecido
0.70 - 0.50 Conductor tratado
deficientemente

En este caso, para las condiciones registradas donde se encuentra la estructura 229, se
selecciona mt = 0.25, y mc = 0.92, obteniéndose así el voltaje crítico disruptivo:

Si el campo eléctrico (gradiente) en la superficie del conductor supera el voltaje crítico

disruptivo, se presenta el efecto corona.

La Tabla 3.13 presenta el voltaje crítico disruptivo y el gradiente crítico disruptivo del aire
para las estructuras seleccionadas que se encuentran a diferentes altitudes sobre el nivel
del mar.

Tabla 3.13. Gradiente crítico disruptivo del aire y el voltaje crítico disruptivo

3.3 Casos de estudio

En la presente sección se variaron los parámetros tanto de la permitividad y la

permeabilidad del aire y del suelo para observar la influencia de la altitud sobre el nivel del
mar sobre los campos eléctricos y magnéticos; de la misma manera se varió la
configuración geométrica de los conductores, distancia entre fases, número de
subconductores por haz, altura de los conductores desde el suelo con el fin de determinar
la configuración óptima la cual pueda reducir los valores de campos electromagnéticos.
3.3.1 VARIACIÓN DE PERMITIVIDADES DEL AIRE Y SUELO

3.3.1.1 Campo eléctrico al variar la permitividad del aire

Los campos eléctricos en las proximidades de una línea de transmisión se calculan

considerando que no existe carga libre en el espacio; de la misma manera se considera el
suelo como un perfecto conductor. A pesar de que, como se cita en [32], la permitividad es
prácticamente independiente de las condiciones climáticas del medio y es
aproximadamente igual a 1.0006, en el presente apartado se observará en qué porcentaje
el campo eléctrico varía si se toman diferentes escenarios como permitividades inferiores
y superiores a la unidad.

En este caso, se considera la permitividad del suelo igual a uno y se varía la permitividad
del aire.

Considérese una línea de transmisión de 525 kV, configuración horizontal, espaciamiento

de 14 metros entre fases, altura del suelo al conductor de 18.5 metros, 4 conductores por
haz, 3.065 centímetros de diámetro, 60 cm de espaciamiento entre subconductores.

Figura 3.10. Campo eléctrico con diferentes permitividades de aire; altura de los
conductores respecto al suelo = 18.5 m

En la Figura 3.10 se grafica el campo eléctrico con diferentes permitividades del aire. El
gráfico obtenido es el perfil transversal de la línea de transmisión, simulado a 1 metro sobre
el suelo. La simulación cubre los 60 metros de ancho de la franja de servidumbre más 20
metros adicionales a la misma con el fin de analizar el comportamiento del campo eléctrico
producido por la línea de transmisión. La Tabla 3.15 presenta el porcentaje de variación de
 valores del campo eléctrico con respecto a la variación de las permitividades, siendo =1
la constante física con la cual se comparan las variaciones. Se puede apreciar que,
mientras menor es la permitividad del aire, el campo eléctrico respecto a = 1 disminuye
con una variación aproximada al 8%, sin embargo, al aumentar la permitividad, el campo
eléctrico varía, pero esta variación no sobrepasa el 1%.

Estudios realizados demuestran que la permitividad del aire es inversamente proporcional

a la altitud sobre el nivel del mar, mientras que es directamente proporcional a la
temperatura y a la presión [46]. Cabe destacar que si bien es cierto la permitividad del aire
puede aumentar o disminuir respecto a la unidad, no lo hace en gran medida. La Tabla
3.14 presenta el cambio de valor del campo eléctrico al variar las permitividades de aire y
la Tabla 3.15 indica el valor de la variación en porcentaje del campo eléctrico al aumentar
o disminuir esta permitividad.

Tabla 3.14. Valor del campo eléctrico en [V/m] a diferentes permitividades del aire

Longitud*
[m] = 0.1 =1 = 10
0 1543.36325 1672.612 1683.09317
26.30872 12349.7489 12492.3832 12404.0556
52.61745 12739.582 12894.7635 12787.1086
78.92617 1800.77292 1961.33458 1981.79152
* La longitud representa el perfil transversal visto desde el centro de la torre hacia sus extremos, abarcando la franja de servidumbre (60 m)
más 20 m adicionales

Tabla 3.15. Porcentaje de variación del campo eléctrico respecto a cambio de permitividades
del aire

Longitud* [m] = 0.1 = 10

0 7.72736 0.62663
26.30872 1.14177 0.70705
52.61745 1.20345 0.83487
78.92617 8.18635 1.04301

Conservando los mismos datos del ejercicio planteado anteriormente, se aumenta la altura
de los conductores con respecto al suelo y se observa la variación de los valores del campo
eléctrico respecto a la permitividad del aire. La altura del conductor respecto al suelo es de
31.5 metros. La Figura 3.11 indica que a mayor permitividad del aire el campo eléctrico
aumenta; de la misma manera ocurre si la permitividad del aire disminuye, decreciendo el
valor del campo eléctrico.
 Figura 3.11. Campo eléctrico con diferentes permitividades de aire; altura de los
conductores respecto al suelo = 31.5 m

La Tabla 3.16 representa el valor del campo eléctrico en [V/m] con respecto al cambio en
la permitividad del aire y la Tabla 3.17 indica el valor de la variación en porcentaje del
campo eléctrico al aumentar o disminuir esta permitividad.

Tabla 3.16. Valor del campo eléctrico en [V/m] a diferentes permitividades del aire

Longitud*
[m] = 0.1 =1 = 10
0 671.57877 716.07396 797.12464
25.50336 4528.73528 4634.46497 4773.3846
50.33557 7666.04124 7770.89859 7840.91654
75.83893 4100.8049 4263.71472 4528.75657
100 628.99852 688.67415 811.33557
* La longitud representa el perfil transversal visto desde el centro de la torre hacia sus extremos, abarcando la franja de servidumbre (60 m)
más 40 m adicionales

Tabla 3.17. Porcentaje de variación del campo eléctrico respecto a cambio de

permitividades del aire

Longitud [m] = 0.1 = 10

0 6.21377 11.31876
25.50336 2.28138 2.99753
50.33557 1.34936 0.90103
75.83893 3.82084 6.21622
100.00000 8.66529 17.81124
En este caso se puede notar que, al igual que el caso anterior analizado (altura del suelo
a los conductores igual a 18.5 m), el valor del campo eléctrico a 1 metro sobre el suelo
aumenta hasta una variación de 17.8%, mientras que, si la permitividad disminuye, existe
una variación del 8.7% respecto a = 1.

3.3.1.2 Campo eléctrico al variar la permitividad del suelo

Dependiendo de la resistividad del suelo, existen diferentes permitividades relativas como

se presenta en la Tabla 3.18 [32]:

Tabla 3.18. Permitividades relativas del suelo

, Permitividad Relativa
Sue -m]
0.2 81
100 10
1000 5

Si se considera que, la Tierra es una fuente infinita de carga, el potencial de referencia de

este gran conductor es cero [44]. Por dicha razón, el campo eléctrico a nivel del suelo es
cero o tiende a ser cero. Cualquiera que sea la permitividad del suelo, esto no afectará al
campo eléctrico a lo largo de una línea de transmisión. Esto se puede observar en la Figura
3.12 en la que, a pesar de un aumento o disminución de la permitividad del suelo, el campo
eléctrico permanece siendo el mismo. Esto se da independientemente si la altura del
conductor al suelo aumenta, como se observa en la Figura 3.13 que, igual que la Figura
3.12, no hay cambio en el campo eléctrico a diferentes permitividades de suelo.

Figura 3.12. Campo eléctrico con diferentes permitividades de suelo, la altura del
conductor al suelo es 18.5 metros
 Figura 3.13. Campo eléctrico con diferentes permitividades de suelo, la altura del
conductor al suelo es 31.5 metros

3.3.2 VARIACIÓN DE PERMEABILIDADES DEL AIRE Y SUELO

3.3.2.1 Campo magnético al variar la permeabilidad del aire

Los campos magnéticos son producidos por la corriente que fluye por un conductor, en
otras palabras, este campo es generado por el movimiento de carga eléctrica en el
conductor. Al aumentar la carga, el campo magnético aumenta y disminuye con el aumento
de la altura del conductor sobre el nivel del suelo. El valor del campo magnético es
directamente proporcional a la corriente de la línea y no sufre perturbación alguna ya sea
con presencia de objetos y personas. Para está simulación se considera la corriente
máxima por fase de 1894 A.

La permeabilidad µ, que es una constante física para una temperatura y densidad de

material determinada, refleja la capacidad de un material o un medio para concentrar las
líneas de flujo magnético. Se denomina permeabilidad magnética a la capacidad de un
medio para atraer y que pasen por este los campos magnéticos.

La permeabilidad es baja en el vacío y es elevada en materiales ferromagnéticos como el

hierro.

En el presente apartado se observó el porcentaje de cambio en el campo magnético si se

toman diferentes escenarios como permeabilidades inferiores y superiores a la unidad.

En este caso, se considera la permeabilidad del suelo igual a la unidad y se varía la

permeabilidad del aire.
Considérese una línea de transmisión de 525 kV, configuración horizontal, espaciamiento
de 14 metros entre fases, altura del suelo al conductor de 18.5 metros, 4 conductores por
haz, 3.065 centímetros de diámetro, 60 cm de espaciamiento entre subconductores.

Figura 3.14. Campo magnético con diferentes permeabilidades del aire; altura de los
conductores respecto al suelo = 18.5 m

En la Figura 3.14 se grafica el campo magnético con diferentes permeabilidades del aire.
El gráfico obtenido es el perfil transversal de la línea de transmisión, simulado a 1 metro
sobre el suelo. La simulación cubre los 60 metros de ancho de la franja de servidumbre y
20 metros adicionales a la misma con el fin de analizar el comportamiento del campo
magnético producido por la línea de transmisión. Se toma de referencia la permeabilidad
del aire igual a 1 que en la Figura 3.14 es representada por la curva de color azul. Al
incrementar el valor de la permeabilidad del aire se tienen valores del campo magnético
bajos comparados con los valores de referencia como se puede apreciar de 0 a 14 metros,
mientras que los valores van incrementando al acercarse más a la fase central de la torre,
desde los 15 metros los valores del campo magnético incrementan y superan los valores
de referencia. El punto máximo es a 40 metros y tiene una variación del 53.5% en su valor.

La Tabla 3.19 presenta los valores del campo magnético obtenidos con diferentes
permeabilidades del aire, mientras que la Tabla 3.20 indica el porcentaje de variación de
valores del campo magnético con respecto a la variación de las permeabilidades, siendo
µ = 1 la constate física con la cual se compara las variaciones.
 Tabla 3.19. Valor del campo magnético en [A/m] a diferentes permeabilidades del aire

u=1 u = 10 u = 100
0 11.71507 6.09261 5.55757
20.9396 25.72179 31.64059 33.02097
40.80537 30.21098 46.38205 49.13412
60.67114 24.78296 29.01921 30.16288
80 11.75677 6.12967 5.58684

Tabla 3.20. Porcentaje de variación del campo magnético respecto a cambio de

permeabilidades del aire
Longitud
[m]
0.00000 47.99340 52.56050
20.93960 23.01084 28.37742
40.80537 53.52713 62.63664
60.67114 17.09340 21.70814
80.00000 47.86260 52.47983

Conservando los mismos datos del ejercicio planteado anteriormente, se aumenta la altura
de los conductores con respecto al suelo, y se observa la variación de valores del campo
magnético respecto a la permeabilidad del aire. La altura del conductor respecto al suelo
es de 31.5 metros.

Figura 3.15. Campo magnético con diferentes permeabilidades del aire; altura de los
conductores respecto al suelo = 31.5 m
La Figura 3.15 es el perfil transversal de la línea de transmisión, simulado a 1 metro sobre
el nivel del suelo. La simulación cubre los 60 metros de ancho de la franja de servidumbre
y 40 metros adicionales a la misma con el fin de analizar el comportamiento del campo
magnético producido por la línea de transmisión.

En este caso al ser la distancia del conductor al suelo más alta comparado con el caso
anterior, el campo magnético es menor. Es decir, el campo magnético disminuye, al
incrementarse la altura del conductor sobre el nivel del suelo.

La Tabla 3.21 presenta los valores del campo magnético obtenidos con diferentes
permeabilidades del aire, mientras que la Tabla 3.22 indica el porcentaje de variación de
valores del campo magnético con respecto a la variación de las permeabilidades, siendo
µ = 1 la constate física con la cual se compara las variaciones.

Tabla 3.21. Valor del campo magnético en [A/m] a diferentes permeabilidades del aire
Longitud
[m] u=1 u = 10 u = 100
0 6.45323 3.20135 2.89762
25.50336 13.96742 16.71069 17.40175
50.33557 17.4885 27.00228 28.59222
75.83893 13.62331 15.76755 16.37757
100 6.45302 3.16793 2.85895

Tabla 3.22. Porcentaje de variación del campo magnético respecto a cambio de

permeabilidades del aire

Longitud
[m]

0.00000 50.39148 55.09810

25.50336 19.64048 24.58814
50.33557 54.40023 63.49161
75.83893 15.73946 20.21721
100.00000 50.90788 55.69598

3.3.2.2 Campo magnético al variar la permeabilidad del suelo

Ahora se considera la permeabilidad del aire igual a la unidad y se varía la permeabilidad

del suelo.
 Figura 3.16. Campo magnético con diferentes permeabilidades del suelo; altura de los
conductores respecto al suelo = 18.5 m

En la Figura 3.16 se grafica el campo magnético con diferentes permeabilidades del suelo.
El gráfico obtenido es el perfil transversal de la línea de transmisión, simulado a 1 metro
sobre el suelo. La simulación cubre los 60 metros de ancho de la franja de servidumbre y
20 metros adicionales a la misma con el fin de analizar el comportamiento del campo
magnético producido por la línea de transmisión. Se toma de referencia la permeabilidad
del suelo igual a 1 que en la Figura 3.16 es representada por la curva de color amarillo. Al
disminuir el valor de permeabilidad se observa que de 0 a 14 metros los valores del campo
magnético son menores a los valores de referencia, a los 15 metros los valores comienzan
a aumentar. El punto mínimo es a 40 metros y tiene una variación del 62,6% en su valor.
La Tabla 3.23 presenta los valores del campo magnético obtenidos con diferentes
permeabilidades del suelo, mientras que la Tabla 3.24 indica el porcentaje de variación de
valores del campo magnético con respecto a la variación de las permeabilidades, siendo
µ = 1 la constante física con la cual se compara las variaciones.

Tabla 3.23. Valor del campo magnético en [A/m] a diferentes permeabilidades del suelo

Longitud
[m] u = 0.01 u = 0.1 u=1
0 5.55755 6.0926 11.71507
20.9396 33.0209 31.64058 25.72179
40.80537 49.13428 46.38207 30.21098
60.67114 30.16284 29.01921 24.78296
80 5.58682 6.12967 11.75677
 Tabla 3.24. Porcentaje de variación del campo magnético respecto a cambio de
permeabilidades del suelo

Longitud [m]
0.00000 52.56068 47.99343
20.93960 28.37715 23.01080
40.80537 62.63716 53.52722
60.67114 21.70795 17.09337
80.00000 52.48002 47.86263

Conservando los mismos datos del ejercicio planteado anteriormente, se aumenta la altura
de los conductores con respecto al suelo, y se observa la variación de valores del campo
magnético respecto a la permeabilidad del suelo. La altura del conductor respecto al suelo
es de 31.5 metros.

Figura 3.17. Campo magnético con diferentes permeabilidades del suelo; altura de los
conductores respecto al suelo = 31.5 m

La Figura 3.17 es el perfil transversal de la línea de transmisión, simulado a 1 metro sobre

el suelo. La simulación cubre los 60 metros de ancho de la franja de servidumbre y 40
metros adicionales a la misma con el fin de analizar el comportamiento del campo
magnético producido por la línea de transmisión.

En este caso al ser la distancia del conductor al suelo más alta comparado con la distancia
del conductor al suelo del caso anterior, el campo magnético es menor. Es decir, el campo
magnético disminuye, al incrementarse la altura del conductor sobre el nivel del suelo.
La Tabla 3.25 presenta los valores del campo magnético obtenidos con diferentes
permeabilidades del aire, mientras que la Tabla 3.26 indica el porcentaje de variación de
valores del campo magnético con respecto a la variación de las permeabilidades, siendo
µ = 1 la constate física con la cual se compara las variaciones.

Tabla 3.25. Valor del campo magnético en [A/m] a diferentes permeabilidades del suelo
Longitud
[m] u = 0.01 u = 0.1 u=1
0 2.89762 3.20135 6.45323
25.50336 17.40176 16.71069 13.96742
50.33557 28.59231 27.00229 17.4885
75.83893 16.37756 15.76755 13.62331
100 2.85894 3.16793 6.45302

Tabla 3.26. Porcentaje de variación del campo magnético respecto a cambio de

permeabilidades del suelo

Longitud [m]

0.00000 55.09815 50.39148

25.50336 24.58824 19.64050
50.33557 63.49211 54.40032
75.83893 20.21715 15.73945
100.00000 55.69608 50.90789

3.3.3 EFECTO DEL CAMBIO DE PARÁMETROS EN LA LÍNEA DE

TRANSMISIÓN

Existen diferentes factores que pueden afectar el valor del campo eléctrico y magnético
bajo una línea de transmisión. Entre ellos está la altura desde el suelo al conductor, cambio
en el número de subconductores por haz, diferentes disposiciones geométricas de los
conductores, distancia entre fases.

3.3.3.1 Variación de la altura de los conductores sobre el nivel del suelo

En esta sección se trabaja con dos alturas diferentes desde el suelo hasta el conductor,
siendo estas 18.5 metros y 31.5 metros tanto para el campo eléctrico como campo
magnético. De la misma manera, se aumenta 1 metro a la distancia de 18.5 metros para
ver la variación existente en campo electromagnético si se incrementa la altura de los
conductores al suelo.
3.3.3.1.1 Efecto en el valor del campo eléctrico

Se considera una permitividad del aire igual a la unidad. En la Figura 3.18 se realizan dos
simulaciones del campo eléctrico a diferentes alturas de los conductores con respecto al
nivel del suelo. La curva de color azul es el resultado del campo eléctrico a una distancia
suelo conductor de 18.5 metros, mientras que la curva de color naranja representa el
campo eléctrico a una distancia de 31.5 metros. De las gráficas resultantes se observa que,
a mayor altura de los conductores, el campo eléctrico a un metro sobre el nivel del suelo
disminuye.

Figura 3.18. Efecto de la altura del conductor al suelo a un metro de altura del suelo con
respecto al campo eléctrico

Para una configuración horizontal, a mayor distancia entre los conductores y el suelo, el
campo eléctrico disminuye. En la gráfica se incrementó la altura de la línea en 13 metros.

En la Figura 3.19 se incrementa la altura de los conductores en 1 metro. Aquí se observa

que, el máximo valor del campo eléctrico a 18.5 metros sobre el suelo es de 14.808 [kV/m],
pero al aumentar en 1 metro la altura, es decir a 19.5 metros, el campo eléctrico disminuye
a 14.052 [kV/m]. Esto es, existe una variación de 757 [V/m] al aumentar la altura del
conductor al suelo en 1 metro.
 Figura 3.19. Efecto de la altura del conductor al suelo a un metro de altura del suelo,
variación de altura igual a 1 metro con respecto al campo eléctrico

Se puede encontrar el valor de cambio del campo eléctrico a medida que se aumenta o
disminuye la altura de la línea mediante la ecuación 3.8 [32]:

(3.8)

Donde:

: campo eléctrico máximo a una altura

: campo eléctrico máximo a una altura

: valor que depende de la configuración de la línea, para líneas simple circuito

con configuración horizontal; para líneas con configuración delta.

3.3.3.1.2 Efecto en el valor del campo magnético

A diferencia del campo eléctrico, el campo magnético no sufre perturbación alguna ya sea
en presencia de objetos o personas. Solamente objetos magnéticos como imanes pueden
cambiar la magnitud del campo magnético en un punto dado. El campo magnético depende
netamente de la corriente que circula por el conductor, por esta razón, para la presente
simulación se ha utilizado la corriente máxima por fase que tiene un valor de 1894 A. Los
mismos datos empleados para el campo eléctrico son utilizados en el presente ejemplo.
Línea de transmisión de 525 kV, configuración horizontal, espaciamiento de 14 metros
entre fases, altura del suelo al conductor de 18.5 metros, 4 subconductores por haz, 3.065
centímetros de diámetro, 60 centímetros de espaciamiento entre subconductores.
En la Figura 3.20 se observa la variación del campo magnético a dos diferentes alturas del
conductor, la primera representada por la línea de color azul indica una altura de la línea
de 18.5 metros, mientras que la de color naranja indica una altura de 31.5 metros. Si bien
es cierto se observa una diferencia entre ambas curvas, el valor del campo magnético en
números no es tan grande. Esto es ya que, el valor del campo magnético no es muy
sensible a las variaciones de altura de la línea.

Figura 3.20. Efecto de la altura del conductor al suelo a un metro de altura del suelo con
respecto al campo magnético

Lo dicho anteriormente se puede constatar con la Figura 3.21, la cual muestra que, el
máximo valor del campo magnético a 18.5 metros sobre el suelo es de 30.193 [A/m] (37.74
[uT]), pero al aumentar en 1 metro la altura, es decir a 19.5 metros, el campo magnético
disminuye a 28.695 [A/m] (35.85 [uT]). Esto es, existe una variación de 1.498 [A/m] (1.873
[uT]) al aumentar la altura del conductor al suelo en 1 metro.

Figura 3.21. Efecto de la altura del conductor al suelo a un metro de altura del suelo,
variación de altura igual a 1 metro con respecto al campo magnético
3.3.3.2 Efecto de la disposición geométrica de los conductores

En esta sección se observó cómo el diseño de la línea de transmisión de extra alto voltaje
puede reducir los impactos negativos de los campos eléctricos y magnéticos.

Una reducción significativa de campos electromagnéticos a lo largo de la línea de

transmisión se da seleccionando una configuración adecuada de la geometría en la
estructura de transmisión.

Se realizaron simulaciones con diferentes configuraciones como son horizontal, vertical,

delta, delta invertido y fase divida.

Considerando la línea de transmisión a simple circuito con la configuración horizontal

existente con el conductor ACAR 1100 MCM 18/19, voltaje máximo (fase-fase) de 525 kV,
corriente máxima por fase de 1894 A, 4 subconductores por haz, diámetro de cada
subconductor D = 3.065 cm, espacio entre subconductores 60 cm, espaciamiento entre
fases S = 14 m, altura del suelo al conductor H = 31.5 m.

En las siguientes simulaciones no se considera los cables de guarda, ya que estos no

tienen un efecto significativo en los campos electromagnéticos a 1 metro del suelo [32].

En la Figura 3.22, Figura 3.23, Figura 3.24, Figura 3.25 y Figura 3.26, se presentan los
esquemas con las diferentes configuraciones de conductores y sus respectivas
dimensiones.

Figura 3.22. Configuración horizontal, conductor ACAR 1100 MCM 18/19

 Figura 3.23. Configuración Vertical, conductor ACAR 1100 MCM 18/19, T = 0.5 S

Figura 3.24. Configuración Delta, conductor ACAR 1100 MCM 18/19, T = S

Figura 3.25. Configuración Delta Invertido, conductor ACAR 1100 MCM 18/19, T = S
Figura 3.26. Configuración Fase Dividida, conductor ACAR 1100 MCM 18/19, T = 0.5
S

Donde:

D : Diámetro del conductor.

H : Altura desde el suelo al conductor.

S : Espaciamiento horizontal entre fases.

T : Espaciamiento vertical entre fases.

En el diseño de las geometrías de los conductores se recomienda utilizar las diferentes

relaciones entre T y S, así [32]:

T=S

T = 0.5 S

T = 1.5 S

En base al conductor ACAR 1100 MCM utilizado en la configuración actual, se realiza una
comparación de los valores de campos eléctricos y magnéticos al aumentar y disminuir el
calibre del conductor a 1000 MCM y 1200 MCM respectivamente.
 ACAR 1000 MCM

Campo eléctrico

Figura 3.27. Campo Eléctrico con diferentes configuraciones de conductores utilizando

ACAR 1000 MCM

En la Figura 3.27 se presenta el comportamiento del campo eléctrico a 1 metro sobre el

nivel del suelo. Se puede observar en la gráfica que la curva de color marrón, la cual
representa a la configuración vertical, tiene los valores más bajos del campo eléctrico,
seguida por la curva de color amarrillo que representa a la configuración delta invertida
T=S. También se observa que la configuración horizontal es la que presenta los valores
más altos del campo eléctrico.

Tabla 3.27. Valores del campo eléctrico en [V/m] con las diferentes configuraciones de
conductores medidos a 1 metro desde el suelo

Delta Delta
Longitud Delta Delta Fase
Invertido Invertido Horizontal Vertical
[m] T=S T = 0.5 S Dividida
T=S T = 0.5 S
0.00000 678.97608 727.06526 499.34003 623.07646 618.90397 748.70304 341.65167
25.50336 4259.60077 4439.12054 3157.27507 3779.06930 4112.79745 4706.71874 2749.64771
50.33557 6511.63085 7038.18307 5701.61254 6524.64236 7349.57084 7800.98588 5543.63734
75.83893 3901.33256 4047.91805 2923.44328 3586.52636 3909.75121 4516.60745 2576.66907
100.00000 656.07591 684.52776 475.15018 631.93690 628.40108 778.62371 348.65600
En la Tabla 3.27, se indican los valores del campo eléctrico obtenidos de la simulación a 1
metro sobre el nivel del suelo de cada una de las configuraciones de conductores. Se
observa el punto medio a 50 metros donde el valor del campo eléctrico es el mayor,
teniendo como resultado en la configuración vertical un valor de 5543.64 [V/m] seguido por
la configuración delta invertido con 5701.61 [V/m], que difieren en 28.9% y 26.9%
respectivamente comparándolos con la configuración horizontal. La configuración
horizontal muestra 7800.98 [V/m] el cual es el valor más alto del campo eléctrico.

Campo magnético

Figura 3.28. Campo Magnético con diferentes configuraciones de conductores

utilizando ACAR 1000 MCM

En la Figura 3.28 se presenta el comportamiento del campo magnético a 1 metro sobre

el nivel del suelo. Se puede observar en la gráfica que la curva de color amarrillo
representa a la configuración delta invertido T=S que muestra los valores más bajos del
campo magnético, seguida por la curva de color verde que representa a la configuración
fase divida, la cual consiste en la división del conductor de fase en más conductores,
esta configuración es caracterizada por reducir significativamente el campo magnético.
También se observa que la configuración horizontal representada por la curva de color
celeste es la que presenta los valores más altos del campo magnético.
Tabla 3.28. Valores del campo magnético en [A/m] con las diferentes configuraciones
de conductores medidos a 1 metro desde el suelo

Delta Delta
Longitud Delta Delta Fase
Invertido Invertido Horizontal Vertical
[m] T=S T = 0.5 S Dividida
T=S T = 0.5 S
0.00000 5.40121 5.93638 3.97482 5.18040 4.64938 6.44734 4.62599
25.50336 11.34001 12.60581 8.58406 11.06701 10.13757 13.95388 10.42437
50.33557 13.45169 15.23549 11.26432 14.16484 13.06224 17.46563 14.15760
75.83893 11.09198 12.30461 8.34898 10.74439 9.86006 13.61733 10.11841
100.00000 5.40087 5.93612 3.97469 5.18108 4.66794 6.44614 4.62566

En la Tabla 3.28 se indican los valores del campo magnético obtenidos de la simulación
a 1 metro sobre el nivel del suelo de cada una de las configuraciones de conductores.
Se observa el punto medio a 50 metros donde el valor del campo magnético es el mayor,
teniendo como resultado en la configuración delta invertido T=S un valor de 11.264 [A/m]
seguido por la configuración fase dividida con 13.062 [A/m] los cuales difieren en 35.51%
y 25.21% respectivamente comparándolos con la configuración horizontal. La
configuración horizontal presenta 17.465 [A/m] el cual es el valor más alto del campo
magnético.

ACAR 1100 MCM

Campo eléctrico

Figura 3.29. Campo Eléctrico con diferentes configuraciones de conductores utilizando

ACAR 1100 MCM
En la Figura 3.29 se presenta el comportamiento del campo eléctrico a 1 metro sobre el
nivel del suelo, se puede observar en la gráfica que la curva de color marrón que
representa a la configuración vertical muestra los valores más bajos del campo eléctrico,
seguida por la curva de color amarrillo que representa a la configuración delta invertida
T=S. También se observa que la configuración horizontal es la que presenta los valores
más altos del campo eléctrico.

Tabla 3.29. Valores del campo eléctrico en [V/m] con las diferentes configuraciones de
conductores medidos a 1 metro desde el suelo

Delta Delta
Longitud Delta Delta Fase
Invertido Invertido Horizontal Vertical
[m] T=S T = 0.5 S Dividida
T=S T = 0.5 S
0.00000 717.04785 748.65930 517.50540 601.37497 620.09417 716.07396 423.10104
25.50336 4317.89088 4534.34844 3181.70788 3743.60846 4113.09957 4634.46497 2874.64916
50.33557 6651.36378 7088.90482 5678.59860 6487.53241 7278.35560 7770.89859 5565.94645
75.83893 4064.15824 4176.18637 3014.30235 3606.53432 3875.37828 4263.71472 2665.84815
100.00000 727.32593 717.50167 512.79929 620.66849 652.33804 688.67415 419.14359

En la Tabla 3.29, se indican los valores del campo eléctrico obtenidos de la simulación a 1
metro sobre el nivel del suelo de cada una de las configuraciones de conductores. Se
observa el punto medio a 50 metros donde el valor del campo eléctrico es el mayor,
teniendo como resultado en la configuración vertical un valor de 5565.94 [V/m] seguido por
la configuración delta invertido con 5678.598 [V/m] que difieren en 28.37% y 26.92%
respectivamente comparándolos con la configuración horizontal. La configuración
horizontal muestra 7770.898 [V/m] el cual es el valor más alto del campo eléctrico.

Campo magnético

Figura 3.30. Campo Magnético con diferentes configuraciones de conductores

utilizando ACAR 1100 MCM
En la Figura 3.30 se presenta el comportamiento del campo magnético a 1 metro sobre
el nivel del suelo. Se puede observar en la gráfica que la curva de color amarrillo
representa a la configuración delta invertido T=S que muestra los valores más bajos del
campo magnético, seguida por la curva de color verde que representa a la configuración
fase divida. También se observa que la configuración horizontal representada por la
curva de color celeste es la que presenta los valores más altos del campo magnético.

Tabla 3.30. Valores del campo magnético en [A/m] con las diferentes configuraciones
de conductores medidos a 1 metro desde el suelo

Delta Delta
Longitud Delta Delta Fase
Invertido Invertido Horizontal Vertical
[m] T=S T = 0.5 S Dividida
T=S T = 0.5 S

0.00000 5.39736 5.93257 3.98015 5.18716 4.57980 6.45323 4.63168

25.50336 11.34013 12.61127 8.59649 11.07592 10.07287 13.96742 10.43321
50.33557 13.44855 15.21186 11.22812 14.17298 13.05737 17.48850 14.16421
75.83893 11.08877 12.29029 8.36092 10.78755 9.74615 13.62331 10.09265
100.00000 5.39698 5.93125 3.98040 5.18751 4.58521 6.45302 4.63183

En la Tabla 3.30, se indican los valores del campo magnético obtenidos de la simulación
a 1 metro sobre el nivel del suelo de cada una de las configuraciones de conductores.
Se observa el punto medio a 50 metros donde el valor del campo magnético es el mayor,
teniendo como resultado en la configuración delta invertido T=S un valor de 11.228 [A/m]
seguido por la configuración fase dividida con 13.057 [A/m] que difieren en 35.8% y
25.34% respectivamente comparándolos con la configuración horizontal. La
configuración horizontal muestra 17.488 [A/m] el cual es el valor más alto del campo
magnético.

ACAR 1200 MCM

Campo eléctrico
Figura 3.31. Campo Eléctrico con diferentes configuraciones de conductores utilizando
ACAR 1200 MCM

En la Figura 3.31 se presenta el comportamiento del campo eléctrico a 1 metro sobre el

nivel del suelo. Se puede observar en la gráfica que la curva de color marrón representa
a la configuración vertical que muestra los valores más bajos del campo eléctrico,
seguida por la curva de color amarrillo que representa a la configuración delta invertida
T=S. También se observa que la configuración horizontal es la que presenta los valores
más altos del campo eléctrico.

Tabla 3.31. Valores del campo eléctrico en [V/m] con las diferentes configuraciones de
conductores medidos a 1 metro desde el suelo

Delta Delta
Longitud Delta Delta Fase
Invertido Invertido Horizontal Vertical
[m] T=S T = 0.5 S Dividida
T=S T = 0.5 S

0.00000 700.29333 729.58839 527.11833 650.31521 677.46097 729.53773 392.66066

25.50336 4261.45000 4448.51510 3252.19742 3914.29199 4215.19816 4651.26799 2845.78092
50.33557 6615.42371 7070.86419 5716.80295 6562.91555 7314.22116 7825.08194 5539.17573
75.83893 4021.93723 4257.36949 2979.49886 3648.72940 3972.20098 4455.87469 2590.68726
100.00000 721.59964 730.98097 515.52207 665.76587 681.70374 761.35213 400.46600

En la Tabla 3.31 se indican los valores del campo eléctrico obtenidos de la simulación a 1
metro sobre el nivel del suelo de cada una de las configuraciones de conductores. Se
observa el punto medio a 50 metros donde el valor del campo eléctrico es el mayor,
teniendo como resultado en la configuración vertical un valor de 5539.18 [V/m] seguido por
la configuración delta invertido con 5716.802 [V/m] los cuales difieren en 29.21% y 26.94%
respectivamente comparándolos con la configuración horizontal. La configuración
horizontal muestra 7825.082 [V/m] el cual es el valor más alto del campo eléctrico.

Campo magnético

Figura 3.32. Campo Magnético con diferentes configuraciones de conductores

utilizando ACAR 1200 MCM

En la Figura 3.32 se presenta el comportamiento del campo magnético a 1 metro sobre

el nivel del suelo. Se puede observar en la gráfica que la curva de color amarrillo
representa a la configuración delta invertido T=S que muestra los valores más bajos del
campo magnético, seguida por la curva de color verde que representa a la configuración
fase divida. También se observa que la configuración horizontal representada por la
curva de color celeste es la que presenta los valores más altos del campo magnético.

Tabla 3.32. Valores del campo magnético en [A/m] con las diferentes configuraciones
de conductores medidos a 1 metro desde el suelo

Delta Delta
Longitud Delta Delta Fase
Invertido Invertido Horizontal Vertical
[m] T=S T = 0.5 S Dividida
T=S T = 0.5 S

0.00000 5.35975 5.89055 3.95251 5.15161 4.48348 6.40742 4.59850

25.50336 11.37521 12.62284 8.60987 11.10649 9.81172 14.02805 10.45191
50.33557 13.52572 15.27474 11.49010 14.33208 12.73945 17.68919 14.36996
75.83893 11.09649 12.33950 8.35488 10.76532 9.50314 13.64795 10.01511
100.00000 5.36524 5.89692 3.95707 5.15706 4.43270 6.41398 4.60456
En la Tabla 3.32 se indican los valores del campo magnético obtenidos de la simulación
a 1 metro sobre el nivel del suelo de cada una de las configuraciones de conductores.
Se observa el punto medio a 50 metros donde el valor del campo magnético es el mayor,
teniendo como resultado en la configuración delta invertido T=S un valor de 11.49 [A/m]
seguido por la configuración fase dividida con 12.74 [A/m] los cuales difieren en 35.04%
y 27.98% respectivamente comparándolos con la configuración horizontal. La
configuración horizontal muestra 17.69 [A/m] el cual es el valor más alto del campo
magnético.

De la Figura 3.27, Figura 3.29 y Figura 3.31 se determina que la configuración óptima la
cual reduce considerablemente el valor de campos eléctricos es la configuración vertical
mientras que de la Figura 3.28, Figura 3.30 y Figura 3.32 se determina que la
configuración óptima para reducir el campo magnético es la configuración delta
invertido.

Tanto para el campo eléctrico como el campo magnético, la configuración horizontal es

la que produce los valores más altos.

Tabla 3.33. Valores del campo eléctrico en [V/m] con las diferentes configuraciones de
conductores ACAR 1000 MCM 18/19, ACAR 1100 MCM 18/19, ACAR 1200 MCM
18/19

Campo Delta Delta

Delta Delta Fase
Eléctrico Invertido Invertido Horizontal Vertical
T=S T = 0.5 S Dividida
[V/m] T=S T = 0.5 S
1000
MCM 6511.63085 7038.18307 5701.61254 6524.64236 7349.57084 7800.98588 5543.63734
18/19
1100
MCM 6651.36378 7088.90482 5678.59860 6487.53241 7278.35560 7770.89859 5565.94645
18/19
1200
MCM 6615.42371 7070.86419 5716.80295 6562.91555 7314.22116 7825.08194 5539.17573
18/19

En la Tabla 3.33 se indican los valores del campo eléctrico obtenidos en la simulación
en el punto medio a 50 metros donde el campo eléctrico es el mayor y a 1 metro sobre
el nivel del suelo de cada una de las configuraciones de conductores. Se comparan los
valores obtenidos con el calibre actual existente ACAR 1100 MCM 18/19 en la línea de
transmisión, aumentando el calibre a ACAR 1200 MCM 18/19 y disminuyendo el calibre
a ACAR 1000 MCM 18/19. Para la configuración delta se tiene que mientras menor sea
el calibre del conductor los valores del campo eléctrico son más bajos. Para la
configuración delta invertido, fase dividida y horizontal con el calibre ACAR 1100 MCM
18/19 los valores del campo eléctrico son más bajos. Para la configuración vertical se
tiene que mientras mayor sea el calibre del conductor el valor del campo eléctrico es
menor. Estas variaciones se dan por la posición en la que se encuentran los conductores
en las diferentes configuraciones.

Para el presente caso de estudio el que menor campo eléctrico produce comparado con
el calibre actual existente es el calibre 1200 MCM, configuración vertical, obteniéndose
una reducción del 28.7%.

Tabla 3.34. Valores del campo magnético en [A/m] con las diferentes configuraciones
de conductores ACAR 1000 MCM 18/19, ACAR 1100 MCM 18/19, ACAR 1200 MCM
18/19

Campo Delta Delta

Delta Delta Fase
Magnético Invertido Invertido Horizontal Vertical
T=S T = 0.5 S Dividida
[A/m] T=S T = 0.5 S
1000 MCM
13.45169 15.23549 11.26432 14.16484 13.06224 17.46563 14.15760
18/19
1100 MCM
13.44855 15.21186 11.22812 14.17298 13.05737 17.48850 14.16421
18/19
1200 MCM
13.52572 15.27474 11.49010 14.33208 12.73945 17.68919 14.36996
18/19

En la Tabla 3.34 se indica los valores del campo magnético obtenidos de la simulación
en el punto medio a 50 metros donde el campo magnético es el mayor y a 1 metro sobre
el nivel del suelo de cada una de las configuraciones de conductores. Se comparan los
valores obtenidos con el calibre actual existente ACAR 1100 MCM 18/19 en la línea de
transmisión, aumentando el calibre a ACAR 1200 MCM 18/19 y disminuyendo el calibre
a ACAR 1000 MCM 18/19. Para la configuración delta y delta invertido se tiene que con
el calibre ACAR 1100 MCM 18/19 el valor del campo magnético es más bajo, mientras
que para la configuración fase dividida se tiene el menor valor del campo magnético
utilizando un calibre mayor. Para la configuración horizontal y vertical se tiene que
mientras menor sea el calibre del conductor el valor del campo magnético es menor.

Para el presente caso de estudio el que menor campo magnético produce comparado
con el calibre actual existente es el calibre 1100 MCM, configuración delta invertido T=S,
obteniéndose una reducción del 35.8%.

En resumen, la geometría óptima que podría reducir los valores de campos

electromagnéticos es la configuración vertical. A pesar de que para campo magnético la
mejor geometría es delta invertida, este campo depende exclusivamente de la corriente
que circula por el conductor, por lo que, en este sentido, las simulaciones del campo
eléctrico son las que deciden la geometría adecuada. Por otro lado, se puede pensar
que, al disminuir el valor del campo eléctrico cerca del conductor, el gradiente potencial
en la superficie del conductor reduciría, haciendo que el efecto corona disminuya y por
ende reduzcan las pérdidas asociadas con la formación de este fenómeno.

En configuración vertical no existen estructuras para simple circuito con la disposición S=0,
por lo que se realiza con la disposición S=T, como se aprecia en la Figura 3.33.

Figura 3.33. Configuración vertical con S = T

De las simulaciones realizadas con disposición S=0 y S=T se tiene como resultado que los
valores obtenidos son los mismos tanto del campo eléctrico como del campo magnético.
Como se observa en la Figura 3.34 y Figura 3.35, respectivamente.

Figura 3.34. Campo eléctrico en [V/m] con la configuración vertical con S=0 y S=T a 1
metro sobre el nivel del suelo.
Figura 3.35. Campo magnético en [A/m] con la configuración vertical con S=0 y S=T a 1
metro sobre el nivel del suelo.

Tabla 3.35. Valores obtenidos de la simulación del campo eléctrico en [V/m]

configuración vertical

Longitud Vertical Vertical

[m] S=0 S=T
0.00000 423.10104 423.10104
25.50336 2874.64916 2874.64916
50.33557 5565.94645 5565.94645
75.83893 2665.84815 2665.84815
100.00000 419.14359 419.14359

Tabla 3.36. Valores obtenidos de la simulación del campo magnético en [A/m]

configuración vertical

Longitud Vertical Vertical

[m] S=0 S=T
0.00000 4.63168 4.63168
25.50336 10.43321 10.43321
50.33557 14.16421 14.16421
75.16779 10.23082 10.23082
100.00000 4.63183 4.63183
En la Tabla 3.35 y Tabla 3.36 se presentan los valores de campos eléctricos y magnéticos
obtenidos a ciertas distancias en el perfil longitudinal de la torre de transmisión.

Para este caso, y como se aprecia en la Tabla 3.35 y Tabla 3.36, la disposición de los
conductores en la configuración vertical no representa ningún cambio en los valores de
campo electromagnético.

3.3.3.3 Separación entre fases

El espaciamiento entre fases juega un papel muy importante a lo que se refiere al campo
electromagnético a nivel del suelo. Líneas más compactas producen menores campos
eléctricos a nivel del suelo [32]. Se realiza una simulación de campos eléctricos y
magnéticos en las diferentes configuraciones de la línea de transmisión a una distancia
entre fases de 14 metros y se aumenta en un metro la distancia entre fases para observar
si existe cambio en el campo electromagnético a 1 metro sobre el nivel del suelo.

Tabla 3.37. Campo eléctrico en [V/m] a 1 metro sobre el nivel del suelo a diferentes
configuraciones; separación entre fases
Campo eléctrico con separación entre fases S = 14 m

Delta Delta
Longitud Delta Delta Fase
Invertido Invertido Horizontal Vertical
[m] T=S T=S/2 Dividida
T=S T=S/2

0 717.047852 748.659301 517.505397 601.374968 620.09417 716.073958 621.736141

25.50 4317.89088 4534.34844 3181.70788 3743.60846 4113.09957 4634.46497 3918.33993
50.34 6651.36378 7088.90482 5678.5986 6487.53241 7278.3556 7770.89859 6025.01671
75.17 4173.94169 4292.5604 3100.95767 3688.97685 3996.44857 4394.42004 2838.14889
100 727.325933 717.501674 512.79929 620.668492 652.338043 688.674153 426.803081
Campo eléctrico con separación entre fases S = 15 m
0 756.272597 769.436128 483.268931 672.67578 681.481371 760.670965 618.096031
25.50 4381.48951 4565.06422 3086.2138 3941.29461 4250.64001 4822.51531 3824.70012
50.34 6562.03105 6996.70595 5627.87924 6568.62294 7337.31614 7788.67675 6040.20004
75.17 4281.74788 4464.22991 3004.49998 3669.53248 4110.48329 4628.79866 2986.6049
100 760.216184 792.020308 493.530806 612.111019 659.659284 729.963638 473.964041

Tabla 3.38. Variación del campo en [V/m] eléctrico a 1 metro sobre el nivel del suelo a
diferentes configuraciones; separación entre fases

Delta Delta
Longitud Delta Delta Fase
Invertido Invertido Horizontal Vertical
[m] T=S T=S/2 Dividida
T=S T=S/2

0 -39.2247 -20.77683 34.236466 -71.300812 -61.3872 -44.597007 3.64011

25.50 -63.5986 -30.71578 95.49408 -197.68615 -137.5404 -188.05034 93.63981
50.34 89.33273 92.19887 50.71936 -81.09053 -58.96054 -17.77816 -15.18333
75.17 -107.806 -171.6695 96.45769 19.44437 -114.0347 -234.37862 -148.4560
100 -32.8903 -74.51863 19.268484 8.557473 -7.321241 -41.289485 -47.16096
En la Tabla 3.37 se puede observar el valor del campo eléctrico a 1 metro sobre el nivel
del suelo cuando se tienen una separación entre fases de 14 y 15 metros. En la Tabla 3.38
se puede observar la variación del campo eléctrico si se aumenta en un metro la separación
entre fases. Cuando se tiene una configuración en delta, exceptuando el punto central
(50.34 m) el campo eléctrico aumenta. Es decir, si se tiene una configuración en delta por
cada aumento de distancia entre las fases, el campo eléctrico aumenta. En el caso de delta
invertido, para cualquier distancia transversal, si se aumenta en 1 metro la separación entre
fases, el campo eléctrico a nivel del suelo disminuye. Para los casos de fase dividida y
configuración horizontal, si la distancia entre fases aumenta, el campo eléctrico aumenta.
Se puede decir que, no para todas las configuraciones el aumento o disminución en la
distancia entre las fases conlleva a la reducción del campo eléctrico a nivel del suelo.
Dependiendo del tipo de configuración, como se describió anteriormente, el campo
eléctrico aumenta o disminuye.

Tabla 3.39. Campo magnético en [A/m] a 1 metro sobre el nivel del suelo a diferentes
configuraciones; separación entre fases
Campo magnético con separación entre fases S = 14 m

Longitud Delta Delta Delta Inv. Delta Inv. Fase

Horizontal Vertical
[m] T=S T=S/2 T=S T=S/2 Dividida

0 5.39736 5.93257 3.98015 5.18716 4.57980 6.45323 5.11284

25.50 11.34013 12.61127 8.59649 11.07592 10.07287 13.96742 11.03794
50.34 13.44855 15.21186 11.22812 14.17298 13.05737 17.48850 13.20087
75.17 11.21119 12.43572 8.47544 10.92880 9.90519 13.79590 9.27678
100 5.39698 5.93125 3.98040 5.18751 4.58521 6.45302 4.21398
Campo magnético con separación entre fases S = 15 m
0 5.42879 5.96282 3.95923 5.18609 4.58598 6.48388 5.06871
25.50 11.25286 12.50536 8.50783 10.98014 10.00172 13.85944 10.92613
50.34 13.09011 14.84792 11.13918 13.99114 12.80304 17.13514 13.11340
75.17 11.13016 12.36333 8.37415 10.83960 9.86556 13.70640 9.31815
100 5.42839 5.96331 3.95900 5.18578 4.59099 6.48430 4.26030

Tabla 3.40. Variación del campo magnético en [A/m] a 1 metro sobre el nivel del suelo a
diferentes configuraciones; separación entre fases

Longitud Delta Delta Delta Inv. Delta Inv. Fase

Horizontal Vertical
[m] T=S T=S/2 T=S T=S/2 Dividida

0 -0.03143 -0.03025 0.02092 0.00107 -0.00618 -0.03065 0.04412

25.50 0.08727 0.10591 0.08866 0.09578 0.07114 0.10798 0.11181
50.34 0.35844 0.36393 0.08894 0.18185 0.25433 0.35335 0.08748
75.17 0.08103 0.07239 0.10129 0.08919 0.03963 0.08950 -0.04136
100 -0.03142 -0.03206 0.02140 0.00172 -0.00578 -0.03127 -0.04632
La Tabla 3.39 representa el campo magnético a 1 metro sobre el nivel del suelo cuando se
tiene una separación entre fases de 14 y 15 metros. Se puede observar que no existe
mayor cambio si se aumenta la distancia entre fases ya que, para el campo magnético, el
valor de este depende de la corriente que circula por el conductor. En la Tabla 3.40 se
observa la variación del campo magnético. La configuración delta invertida es la que mejor
respuesta a la separación entre fases presenta, esto es si se aumenta en 1 metro la
distancia entre fases el campo magnético disminuye.

Los resultados mostrados en la Tabla 3.39 y Tabla 3.40 muestran una relación directa entre
la separación de fases, es decir que, para algunas configuraciones, si se aumenta la
distancia entre fases el campo eléctrico y/o magnético también lo hace, y viceversa.

3.3.3.4 Cambio número de subconductores por Haz

Para el presente estudio se han considerado diferentes calibres de conductor según el

número de subconductores en el haz a fin de poder observar cómo afecta este parámetro
en los valores de campos electromagnéticos.

Para determinar el número de subconductores por haz se consideran las condiciones más
desfavorables de operación de la línea de transmisión:

Voltaje máximo de la línea en estado estable (fase - fase): 525 kV

Corriente máxima por fase: 1894 A

Para calcular la corriente mínima de cada subconductor del haz se divide la corriente
máxima por fase entre el número de subconductores, según se indica la ecuación 3.9:

(3.9)

Donde:

: Corriente mínima por subconductor

: Corriente máxima por fase

: número de subconductores por fase

Para la simulación del campo magnético se requiere ingresar el valor de la densidad de

corriente. Se realiza el cálculo de la densidad de corriente para cada uno de los
conductores utilizados empleando la ecuación 3.9:
Ejemplo de cálculo:

Cálculo del área de la sección transversal del conductor a utilizar en este caso ACAR
1100 MCM 18/19 a través de la ecuación 3.10.

(3.10)

Cálculo de la densidad de corriente se realiza utilizando la ecuación 3.11

(3.11)

Para la selección del calibre del conductor se considera que la capacidad de corriente del
subconductor sea mayor a la capacidad mínima de corriente que puede transportar el
subconductor.

La Tabla 3.41 muestra los calibres de conductor seleccionados para la simulación y análisis
de los campos electromagnéticos bajo la línea de transmisión.

Tabla 3.41. Calibre de conductor para diferentes números de subconductores en el haz

No. Subconductores en el haz Conductor ACAR

3 1700 MCM 42/19
4 1100 MCM 18/19
5 800 MCM 18/19

La Tabla 3.42 muestra los parámetros que se obtienen considerando diferentes calibres
para el haz de subconductores basados en los cálculos descritos anteriormente.
 Tabla 3.42. Parámetros calculados para los calibres de los subconductores en el haz

Densidad
Calibre del Diámetro Radio Área Corriente de
Conductores [m] [m] [m2] [A] corriente
[MA/m2]
800 MCM 0.02815 0.014075 0.000622367 378.800 0.60864392
1000 MCM 0.02923 0.014615 0.000671039 473.500 0.70562258
1100 MCM 0.03065 0.015325 0.000737821 473.500 0.64175483
1200 MCM 0.03202 0.01601 0.000805253 473.500 0.58801370
1700 MCM 0.03816 0.01908 0.001143685 631.333 0.55201656

3.3.3.4.1 Campo eléctrico, configuración horizontal, diferente número de

subconductores en el haz

En la Figura 3.36 se observa el campo eléctrico obtenido bajo el conductor a 1 metro de

altura sobre el suelo para diferentes subconductores del haz. Se puede apreciar que, la
curva que representa el menor campo eléctrico es la que posee 3 subconductores por haz
los mismos que tienen un calibre de 1700 MCM. Mientras que para n = 4 y n = 5, el campo
eléctrico es superior.

Figura 3.36. Campo eléctrico con conductor ACAR 1100 MCM, n = 4; ACAR 800 MCM,
n = 5; ACAR 1700 MCM, n = 3

En efecto, conforme el diámetro del conductor disminuye, el campo eléctrico a nivel del
suelo aumenta. Si se aumenta el número de subconductores, el campo eléctrico a nivel del
suelo aumenta. En la Tabla 3.43 se puede observar que el incremento del campo eléctrico
es proporcional al número de subconductores.
Tabla 3.43. Valores del campo eléctrico en [V/m] para n = 3, n = 4, n = 5 subconductores
por haz

Longitud n=3 n=4 n=5

[m] 1700 MCM 1100 MCM 800 MCM
0 696.37661 716.07396 758.94513
25.50336 4437.58852 4634.46497 4878.99648
50.33557 7397.32597 7770.89859 8119.45915
75.16779 4153.28192 4394.42004 4514.84775
100.00000 646.44359 688.67415 670.77404

3.3.3.4.2 Campo magnético, configuración horizontal, diferente número de

subconductores en el haz

En la Figura 3.37 se aprecia el cambio del valor del campo magnético respecto al número
de subconductores del haz. El campo magnético al depender únicamente de la corriente
que circula por el conductor no tiene una dependencia directa del calibre del conductor. En
la Figura 3.37 se indica que las curvas de color azul y rojo no presentan mayores cambios,
mientras que la curva de color amarillo que es la que posee 5 subconductores por haz
presenta la mayor brecha en comparación con n = 3 y n = 4. Esto quiere decir que, el
arreglo de 5 subconductores por haz es el que menor campo magnético presenta a un
metro sobre el nivel del suelo.

Figura 3.37. Campo magnético con conductor ACAR 1100 MCM, n = 4; ACAR 800 MCM,
n = 5; ACAR 1700 MCM, n = 3
La Tabla 3.44 muestra que no existe un mayor cambio en los valores del campo magnético
al cambiar el número de subconductores por haz, por lo que, este factor no decide al
momento de seleccionar el calibre adecuado para una línea de transmisión.

Tabla 3.44. Valores del campo magnético en [A/m] para n = 3, n = 4, n = 5

subconductores por haz
n=3 n=4 n=5
Longitud [m]
1700 MCM 1100 MCM 800 MCM
0 7.86938 7.90872 6.80197
25.50336 17.02900 17.11769 14.90329
50.33557 21.32948 21.43292 18.79797
75.16779 16.77567 16.90748 14.68478
100 7.86965 7.90847 6.81020

El número de subconductores del haz que mejor se ajusta en la reducción de los campos
electromagnéticos, para el presente estudio es el haz conformado por 3 subconductores
los mismos que tienen un calibre de 1700 MCM.

Por lo tanto, la configuración óptima a modelar tendrá las siguientes características:

3 subconductores por haz

Configuración vertical S = T
Calibre del conductor: 1700 MCM 42/19

3.3.3.5 Variación del número de subconductores en diferentes fases

Con los casos de simulación presentados anteriormente y a su vez conociendo que

mientras menor es el número de sub conductores en el haz, el calibre es mayor y se
reducen los valores de campos electromagnéticos. Se tienen los siguientes casos de
estudio donde se considera:

Caso 1: 3 subconductores por haz, calibre de 1700 MCM para la Fase A y Fase C; 4
sub conductores por haz, calibre de 1100 MCM en la Fase B.
Caso 2: 3 subconductores por haz, calibre de 1700 MCM para la Fase B; 4 sub
conductores por haz, calibre de 1100 MCM en la Fase A y Fase C.
3.4.3.5.1 Variación del número y calibre de subconductores en la Fase A y
Fase C

Figura 3.38. Variación de la disposición de sub conductores de la Fase A y Fase C

En la Figura 3.38 se tiene la disposición geométrica para cada una de las tres fases. Para
la Fase A y Fase C se usan 3 subconductores por haz con un calibre de 1700 MCM, y así
mismo para la Fase B se conservan los 4 subconductores por haz con un calibre de 1100
MCM.

Campo Magnético

Figura 3.39. Campo Magnético con diferentes configuraciones de conductores

utilizando la disposición de 4 subconductores por haz para la Fase B y 3
subconductores por haz para la Fase A y Fase C.

En la Figura 3.39 se tiene las curvas resultantes de los valores del campo magnético
con diferentes configuraciones geométricas teniendo como resultado que la
configuración óptima para este caso es delta invertido T=S la cual contiene los valores
del campo magnético más bajos.
Tabla 3.45. Valores del campo magnético en [A/m] para las diferentes configuraciones
de conductores medidos a 1 metro desde el suelo
4 subconductores por haz Fase A, Fase B y Fase C (1100 MCM)
Delta Delta
Longitud Delta Delta
Invertido Invertido Horizontal Vertical
[m] T=S T = 0.5 S
T=S T = 0.5 S
0 5.39736 5.93257 3.98015 5.18716 6.45323 4.63168
25.50336 11.34013 12.61127 8.59649 11.07592 13.96742 10.43321
50.33557 13.44855 15.21186 11.22812 14.17298 17.4885 14.16421
75.83893 11.08877 12.29029 8.36092 10.78755 13.62331 10.09265
100 5.39698 5.93125 3.9804 5.18751 6.45302 4.63183
Caso 1: 3 subconductores por haz Fase A (1700 MCM), Fase C (1700 MCM) y 4
subconductores por haz Fase B (1100 MCM)
0 6.41533 7.31501 3.96987 6.3709 6.39899 1.75223
25.50336 10.50024 11.70507 8.68327 10.49809 13.47628 10.5334
50.33557 13.22571 14.68732 11.22854 13.85578 17.109 13.26099
75.83893 10.79378 11.88441 7.75162 9.6814 13.32661 9.88522
100 5.52272 6.05073 4.16132 6.01259 6.5202 4.71712

En la Tabla 3.45 se encuentran los valores obtenidos de la simulación a 1 metro de altura

sobre el nivel del suelo, tanto para la disposición de las tres fases con 4 subconductores,
así como cambiando la Fase A y Fase C por 3 subconductores por haz. Se tiene como
resultado que la configuración delta invertido T=S tiene los valores más bajos en
comparación al resto de configuraciones.

Tabla 3.46. Variación del campo magnético en [A/m] comparando el caso de las tres
fases con 4 sub conductores por haz con el caso de la Fase A y Fase C con 3 sub
conductores por haz.
Delta Delta
Longitud Delta Delta
Invertido Invertido Horizontal Vertical
[m] T=S T = 0.5 S
T=S T = 0.5 S
0 -1.01797 -1.38244 0.01028 -1.18374 0.05424 2.87945
25.50336 0.83989 0.9062 -0.08678 0.57783 0.49114 -0.10019
50.33557 0.22284 0.52454 -0.00042 0.3172 0.3795 0.90322
75.83893 0.29499 0.40588 0.6093 1.10615 0.2967 0.20743
100 -0.12574 -0.11948 -0.18092 -0.82508 -0.06718 -0.08529

La Tabla 3.46 presenta la variación existente en los valores del campo magnético en [A/m]
obtenidos de la simulación con 3 subconductores por haz en la Fase A y Fase C
comparándolos con la simulación de las tres fases con 4 subconductores por haz. Se
evidencia que para la configuración delta invertido T=S la variación no supera los 0.7 [A/m].
A 25, 50 y 100 metros se tiene una variación negativa lo que significa que el caso con 3
subconductores por haz en la Fase A y Fase C muestra valores de campos magnéticos
más bajos comparados con el caso de 4 subconductores por haz en las tres fases.

En el resto de las configuraciones la variación en cada caso es mínima, no supera el 1.5

[A/m].
Campo Eléctrico

Figura 3.40. Campo Eléctrico con diferentes configuraciones de conductores utilizando

la disposición de 4 subconductores por haz para la Fase B y 3 subconductores por haz
para la Fase A y Fase C.

En la Figura 3.40 se tiene las curvas resultantes de los valores del campo eléctrico con
diferentes configuraciones geométricas. Se tiene como resultado que la configuración
óptima para este caso es vertical la cual contiene los valores del campo eléctrico más
bajos.

Tabla 3.47. Valores del campo eléctrico en [V/m] con las diferentes configuraciones de
conductores medidos a 1 metro desde el suelo
4 subconductores por haz Fase A, Fase B y Fase C (1100 MCM)
Delta Delta
Longitud Delta Delta
Invertido Invertido Horizontal Vertical
[m] T=S T = 0.5 S
T=S T = 0.5 S
0 717.04785 748.6593 517.5054 601.37497 716.07396 423.10104
25.50336 4317.89088 4534.34844 3181.70788 3743.60846 4634.46497 2874.64916
50.33557 6651.36378 7088.90482 5678.5986 6487.53241 7770.89859 5565.94645
75.83893 4064.15824 4176.18637 3014.30235 3606.53432 4263.71472 2665.84815
100 727.32593 717.50167 512.79929 620.66849 688.67415 419.14359
Caso 1: 3 subconductores por haz Fase A (1700 MCM), Fase C (1700 MCM) y 4
subconductores por haz Fase B (1100 MCM)
0 1025.8177 1029.84011 773.54939 916.4074 1092.149 465.11538
25.50336 4076.9084 4404.45936 2960.40594 3427.71885 4572.07324 2756.57829
50.33557 6386.23313 6897.49611 5452.97502 6354.74434 7394.52353 5322.17356
75.83893 4352.43103 4273.42679 3021.52769 3661.20056 5211.84974 2829.39468
100 617.53676 666.85175 693.61113 833.11521 1106.35512 381.02711

En la Tabla 3.47 se encuentran los valores obtenidos de la simulación a 1 metro de altura

sobre el nivel del suelo, tanto para la disposición de las tres fases con 4 subconductores,
así como cambiando la Fase A y Fase C por 3 subconductores por haz. Se tiene como
resultado que la configuración vertical tiene los valores más bajos seguida por delta
invertido T=S, en comparación al resto de configuraciones.

Tabla 3.48. Variación del campo eléctrico en [V/m] comparando el caso de las tres fases
con 4 subconductores por haz con el caso de la Fase A y Fase C con 3 subconductores
por haz

Delta Delta
Longitud Delta Delta
Invertido Invertido Horizontal Vertical
[m] T=S T = 0.5 S
T=S T = 0.5 S
0 -308.76985 -281.18081 -256.04399 -315.03243 -376.07504 -42.01434
25.50336 240.98248 129.88908 221.30194 315.88961 62.39173 118.07087
50.33557 265.13065 191.40871 225.62358 132.78807 376.37506 243.77289
75.83893 -288.27279 -97.24042 -7.22534 -54.66624 -948.13502 -163.54653
100 109.78917 50.64992 -180.81184 -212.44672 -417.68097 38.11648

La Tabla 3.48 presenta la variación existente en los valores del campo eléctrico en [V/m]
obtenidos de la simulación con 3 subconductores por haz en la Fase A y Fase C
comparándolos con la simulación de las tres fases con 4 subconductores por haz. Se
evidencia que para la configuración vertical la variación no supera los 250 [V/m]. A 0 y 75
metros se tiene una variación negativa lo que significa que el caso con 3 subconductores
por haz en la Fase A y Fase C muestra valores de campos magnéticos más bajos
comparados con el caso de 4 subconductores por haz en las tres fases.

En el resto de las configuraciones la variación en cada caso no supera los 1000 [V/m].

3.3.3.5.2 Variación del número y calibre de subconductores en la Fase B

Figura 3.41. Variación de la disposición de subconductores de la Fase B

En la Figura 3.41 se tiene la disposición geométrica para cada una de las tres fases,
teniendo que para la Fase B se usan 3 subconductores por haz con un calibre de 1700
MCM, y así mismo para la Fase A y Fase C se conserva los 4 subconductores por haz con
un calibre de 1100 MCM.

Campo Magnético
 Figura 3.42. Campo Magnético con diferentes configuraciones de conductores
utilizando la disposición de 4 subconductores por haz para la Fase A y Fase C y 3
subconductores por haz para la Fase B

En la Figura 3.42 se tienen las curvas resultantes de los valores del campo magnético
con diferentes configuraciones geométricas. Se tiene como resultado que la
configuración óptima para este caso es delta invertido T=S la cual contiene los valores
del campo magnético más bajos.

Tabla 3.49. Valores del campo magnético en [A/m] con las diferentes configuraciones
de conductores medidos a 1 metro desde el suelo

4 subconductores por haz Fase A, Fase B y Fase C (1100 MCM)

Delta Delta
Longitud Delta Delta
Invertido Invertido Horizontal Vertical
[m] T=S T = 0.5 S
T=S T = 0.5 S
0 5.39736 5.93257 3.98015 5.18716 6.45323 4.63168
25.50336 11.34013 12.61127 8.59649 11.07592 13.96742 10.43321
50.33557 13.44855 15.21186 11.22812 14.17298 17.4885 14.16421
75.83893 11.08877 12.29029 8.36092 10.78755 13.62331 10.09265
100 5.39698 5.93125 3.9804 5.18751 6.45302 4.63183
4 subconductores por haz Fase A (1100 MCM), Fase C (1100 MCM) y 3
subconductores por haz Fase B (1700 MCM)
0 5.18468 6.73263 4.91793 6.05544 7.42683 5.08903
25.50336 10.93883 8.8621 8.28486 10.45656 12.9421 9.03153
50.33557 13.39116 13.31067 11.53893 13.46855 17.17257 12.73967
75.83893 10.86213 11.89067 8.381 9.86962 13.29604 8.92791
100 6.28358 5.87211 4.36548 5.38411 6.46704 4.46379

En la Tabla 3.49 se encuentran los valores obtenidos de la simulación a 1 metro de altura

sobre el nivel del suelo, tanto para la disposición de las tres fases con 4 subconductores,
así como cambiando la Fase B por 3 subconductores por haz. Se tiene como resultado que
la configuración delta invertido T=S tiene los valores más bajos en comparación al resto de
configuraciones.

Tabla 3.50. Variación del campo magnético en [A/m] comparando el caso de las tres
fases con 4 subconductores por haz con el caso de la Fase B con 3 subconductores por
haz

Delta Delta
Longitud Delta Delta
Invertido Invertido Horizontal Vertical
[m] T=S T = 0.5 S
T=S T = 0.5 S
0 0.21268 -0.80006 -0.93778 -0.86828 -0.9736 -0.45735
25.50336 0.4013 3.74917 0.31163 0.61936 1.02532 1.40168
50.33557 0.05739 1.90119 -0.31081 0.70443 0.31593 1.42454
75.83893 0.22664 0.39962 -0.02008 0.91793 0.32727 1.16474
100 -0.8866 0.05914 -0.38508 -0.1966 -0.01402 0.16804

La Tabla 3.50 presenta la variación existente en los valores del campo magnético en [A/m]
obtenidos de la simulación con 3 subconductores por haz en la Fase B comparándolos con
la simulación de las tres fases con 4 subconductores por haz. Se evidencia que para la
configuración delta invertido T=S la variación no supera el 1 [A/m]. A 0, 50, 75 y 100 metros
se tiene una variación negativa lo que significa que el caso con 3 subconductores por haz
en la Fase B muestra valores de campos magnéticos más bajos comparados con el caso
de 4 subconductores por haz en las tres fases.

En el resto de las configuraciones la variación en cada caso es mínima, no supera los 5

[A/m].

Campo Eléctrico

Figura 3.43. Campo Eléctrico con diferentes configuraciones de conductores utilizando

la disposición de 4 subconductores por haz para la Fase A y Fase C y 3
subconductores por haz para la Fase B
En la Figura 3.43 se tiene las curvas resultantes de los valores del campo eléctrico con
diferentes configuraciones geométricas. Se tiene como resultado que la configuración
óptima para este caso es vertical la cual contiene los valores del campo eléctrico más
bajos.

Tabla 3.51. Valores del campo eléctrico en [V/m] con las diferentes configuraciones de
conductores medidos a 1 metro desde el suelo

4 subconductores por haz Fase A, Fase B y Fase C (1100 MCM)

Delta Delta
Longitud Delta Delta
Invertido Invertido Horizontal Vertical
[m] T=S T = 0.5 S
T=S T = 0.5 S
0 717.04785 748.6593 517.5054 601.37497 716.07396 423.10104
25.50336 4317.89088 4534.34844 3181.70788 3743.60846 4634.46497 2874.64916
50.33557 6651.36378 7088.90482 5678.5986 6487.53241 7770.89859 5565.94645
75.83893 4064.15824 4176.18637 3014.30235 3606.53432 4263.71472 2665.84815
100 727.32593 717.50167 512.79929 620.66849 688.67415 419.14359
4 subconductores por haz Fase A (1100 MCM), Fase C (1100 MCM) y 3 subconductores
por haz Fase B (1700 MCM)
0 997.38903 1073.49747 792.73893 884.95137 1042.40214 627.2808
25.50336 4240.80066 4369.25986 2955.3149 3989.75248 4835.71535 3145.68256
50.33557 6710.12246 6987.74603 5796.19826 6443.44026 7704.92752 5475.35604
75.83893 4138.61328 4237.76118 3491.38421 4294.38042 4334.36789 2430.94142
100 1081.94141 668.65337 489.53489 887.16118 1090.9144 541.80744

En la Tabla 3.51 se encuentran los valores obtenidos de la simulación a 1 metro de altura

sobre el nivel del suelo, tanto para la disposición de las tres fases con 4 subconductores,
así como cambiando la Fase B por 3 subconductores por haz. Se tiene como resultado que
la configuración vertical tiene los valores más bajos seguida por delta invertido T=S, en
comparación al resto de configuraciones.

Tabla 3.52. Variación del campo eléctrico en [V/m] comparando el caso de las tres fases
con 4 subconductores por haz con el caso de la Fase B con 3 subconductores por haz

Delta Delta
Longitud Delta Delta
Invertido Invertido Horizontal Vertical
[m] T=S T = 0.5 S
T=S T = 0.5 S
0 -280.34118 -324.83817 -275.23353 -283.5764 -326.32818 -204.17976
25.50336 77.09022 165.08858 226.39298 -246.14402 -201.25038 -271.0334
50.33557 -58.75868 101.15879 -117.59966 44.09215 65.97107 90.59041
75.83893 -74.45504 -61.57481 -477.08186 -687.8461 -70.65317 234.90673
100 -354.61548 48.8483 23.2644 -266.49269 -402.24025 -122.66385

La Tabla 3.52 presenta la variación existente en los valores del campo eléctrico en [V/m]
obtenidos de la simulación con 3 subconductores por haz en la Fase B comparándolos con
la simulación de las tres fases con 4 subconductores por haz. Se evidencia que para la
configuración vertical la variación no supera los 280 [V/m]. A 0 ,25 y 75 metros se tiene una
variación negativa lo que significa que el caso con 3 subconductores por haz en la Fase B
muestra valores de campos magnéticos más bajos comparados con el caso de 4
subconductores por haz en las tres fases.

En el resto de las configuraciones la variación en cada caso no supera los 700 [V/m].

Con el análisis realizado anteriormente, se concluye que tanto para la disposición de 3

subconductores en el haz para la Fase A y Fase C y 4 subconductores en el haz para la
Fase B, como para la disposición de 4 subconductores en el haz para la Fase A y Fase C
y 3 subconductores en el haz para la Fase B, no se observa un gran variación en los valores
obtenidos de campo electromagnético por lo que se usa la disposición de 3 subconductores
en las tres fases, donde se observa que presenta una disminución relevante en sus
resultados de valores de campo electromagnético.

Las tablas de la sección 3.3 presentan valores de campos eléctricos que exceden, en
algunos casos, los límites máximos establecidos por la norma IEEE Std. 644-1994,
mientras que, para todos los casos, el campo magnético no supera los límites máximos de
exposición.

3.3.3.6 Variación del número de subconductores para la configuración actual

existente con calibre ACAR 1100 MCM

En esta sección se realiza la simulación de los campos eléctricos y magnéticos con 3 y 5

subconductores por haz con calibre 1100 MCM, configuración horizontal, con el fin de
comparar la variación de los valores de campos eléctricos y magnéticos con la disposición
geométrica actual existente en las líneas de transmisión de 500 kV, que tiene 4
subconductores por haz, calibre 1100 MCM, configuración horizontal.

Los parámetros tanto de voltaje como de densidad de corriente utilizados para la simulación
del campo eléctrico y del campo magnético se presentan en la Tabla 3.53.

Tabla 3.53. Parámetros utilizados en la simulación de los campos eléctricos y magnéticos

Campo Campo
Eléctrico Magnético
Número de
Corriente
subconductores Voltaje RMS Densidad de
máxima por
en el haz fase-neutro corriente
subconductor
[kV] [MA/m2]
[A]
3 428.660 631.33 0.85533
4 428.660 473.50 0.64175
5 428.660 378.80 0.51340
Campo eléctrico

Figura 3.44. Campo eléctrico para una configuración horizontal, con 3, 4 y 5

subconductores por haz, ACAR 1100 MCM

En la Figura 3.44 se presenta el campo eléctrico a 1 metro sobre el nivel del suelo para
una configuración horizontal, tanto para 3, 4 y 5 subconductores por haz, con el mismo
calibre en cada subconductor. Se observa que, a menor número de subconductores el
campo eléctrico disminuye. Estos valores se encuentran en la Tabla 3.54.

Tabla 3.54. Valores del campo eléctrico en [V/m] para diferente número de
subconductores por haz, configuración horizontal, calibre 1100 MCM

Longitud
n=3 n=4 n=5
[m]
0.00000 1776.73322 1887.60879 2075.60576
20.40268 5330.24332 5649.25685 6062.59531
40.26846 7325.37734 7771.43144 8243.27406
60.13423 5131.87857 5444.70485 5762.29103
80.00000 1692.38413 1773.68899 1876.59569
 Tabla 3.55. Variación de los valores del campo eléctrico en [V/m] para diferente número
de subconductores por haz, configuración horizontal, calibre 1100 MCM

Longitud
[m]

0.00000 5.87386 -9.95953

20.40268 5.64700 -7.31669
40.26846 5.73966 -6.07150
60.13423 5.74551 -5.83294
80.00000 4.58394 -5.80185

En comparación con 4 subconductores por haz, calibre 1100 MCM, si el número de

subconductores disminuye el campo eléctrico a 1 metro sobre el nivel del suelo también
disminuye. La variación se presenta hasta en un 5.8%. De la misma manera, si se aumenta
el número de subconductores el campo eléctrico también lo hace hasta alcanzar una
variación aproximada del 10%. El porcentaje de variación para n = 3 y n = 5 comparado
con n = 4 se presenta en la Tabla 3.55.

Campo magnético

Figura 3.45. Campo magnético para una configuración horizontal, con 3, 4 y 5

subconductores por haz, ACAR 1100 MCM
En la Figura 3.45 se presenta el campo magnético a 1 metro sobre el nivel del suelo para
una configuración horizontal, tanto para 3, 4 y 5 subconductores por haz, con el mismo
calibre en cada subconductor. Se indica que las tres curvas no presentan mayores
cambios, sin embargo, para n = 3 subconductores el campo magnético es menor, mientras
que para n = 5 el campo magnético es mayor. Estos valores se encuentran en la Tabla
3.56.

Tabla 3.56. Valores del campo magnético en [A/m] para diferente número de
subconductores por haz, configuración horizontal, calibre 1100 MCM

Longitud
n=3 n=4 n=5
[m]
0.00000 11.550227 11.613264 11.666999
20.40268 18.497283 18.552227 18.715415
40.26846 21.294388 21.433928 21.509613
60.13423 18.355177 18.459608 18.613738
80.00000 11.560907 11.620960 11.677404

Tabla 3.57. Variación de los valores del campo magnético en [A/m] para diferente número
de subconductores por haz, configuración horizontal, calibre 1100 MCM

Longitud
[m]

0.00000 0.54281 -0.46270

20.40268 0.29616 -0.87962
40.26846 0.65102 -0.35311
60.13423 0.56572 -0.83496
80.00000 0.51676 -0.48571

En la Tabla 3.57 se puede observar que la variación del campo magnético al aumentar o
disminuir el número de subconductores por haz es despreciable puesto que no llega a
superar el 1% de variación.

3.4 Pérdidas por efecto corona

En la sección 3.3.3.6 se realizó un análisis de los valores simulados tanto de campo

eléctrico como de campo magnético para una línea de transmisión de 500 kV con
configuración horizontal con 3, 4 y 5 subconductores por haz respectivamente. Ahora se
calculan las pérdidas por efecto corona para cada caso analizado y se comparan dichos
valores con la configuración actual existente a fin de determinar si aumentan o disminuyen
dichas pérdidas variando el número de subconductores en el haz.

Tabla 3.58. Cálculos de pérdidas por efecto corona para n = 3, 4, 5; ACAR 1100 MCM

Factor Voltaje
Número de Radio círculo Voltaje Pérdida
de de inicio
sub de sub RMG DMG efectivo corona
disposi- de
conductores conductores [cm] [cm] fase-fase kW/km-
ción de corona
por haz [cm] [kV] fase
haces B [kV]

3 34.641 17.670 1763.889 371.231 0.363 104.653 1053.182

4 42.426 26.157 1763.889 371.231 0.277 125.357 1404.243
5 51.039 34.930 1763.889 371.231 0.224 144.405 1755.304

En la Tabla 3.58 se presentan las pérdidas por efecto corona conservando el mismo calibre
para cada subconductor del haz, variando el número de subconductores a 3, 4 y 5. Los
parámetros utilizados para el cálculo de pérdidas por el efecto corona son en base a las
condiciones atmosféricas que se presentan a mayor altitud sobre el nivel del mar, en este
caso estas condiciones se dan donde se encuentra ubicada la estructura 229 analizada
anteriormente.

Los valores resultantes de pérdidas por efecto corona varían en 351.061 kW/km-fase
respectivamente obteniéndose que con 3 subconductores las pérdidas por efecto corona
son menores y con 5 subconductores por haz estas pérdidas son mayores. El valor de
pérdidas por efecto corona aumenta conforme incrementa el número de subconductores.

Una vez conocida la mejor configuración de disposición geométrica de línea de transmisión

de extra alto voltaje que reduce los valores de campos eléctricos y magnéticos, se analiza
si efectivamente con la configuración propuesta se reduce el efecto corona. Es decir, se
presentan los valores de pérdidas por efecto corona con la configuración actual (disposición
horizontal, 4 subconductores por haz ACAR 1100 MCM) y con la configuración propuesta
(disposición vertical, 3 subconductores por haz ACAR 1700 MCM), y se observa cuál de
estas configuraciones tiene menores pérdidas por efecto corona.

Aplicando las ecuaciones descritas en el subcapítulo 1.3.10 (Efecto corona y pérdidas), se

tiene lo que se describe a continuación:
 Tabla 3.59. Cálculos de los parámetros del conductor de la línea de transmisión

Número de Distancias (m) Diámetro Radio Radio círculo de

DMG RMG
subconductores Conductor Conductor subconductores
Entre Entre Entre (cm) (cm)
por haz n (cm) (cm) (cm)
AyB ByC AyC

3 15.652 15.652 14 3.816 1.908 34.641 1508.074 19.009

4 14 14 28 3.065 1.5325 42.426 1763.889 26.157

En la Tabla 3.59 se indican las distancias de separación entre las fases tanto para 4
subconductores y 3 subconductores por haz, el diámetro y radio de cada uno de los
conductores utilizados, la distancia media geométrica DMG y el radio medio geométrico
RMG.

Tabla 3.60. Cálculo del Voltaje de inicio de corona, en kV

Número de Factor de Factor Factor Voltaje de

subconductores Estructura disposición ambiental superficial inicio de
por haz n de haces mt mc corona, en kV
4 254 0.2771 1 0.92 365.892
4 229 0.2771 0.25 0.92 75.174
4 56 0.2771 0.8 0.92 404.916
3 254 0.3701 1 0.92 322.950
3 229 0.3701 0.25 0.92 66.352
3 56 0.3701 0.8 0.92 357.394

En la Tabla 3.60 se presenta el valor del Voltaje de inicio de corona calculado para cada
una de las estructuras analizadas, considerando un haz de 4 subconductores y un haz de
3 subconductores.

Tabla 3.60. Cálculo de Pérdidas por Efecto Corona

Voltaje Voltaje de Pérdida

Número de Densidad Radio del
efectivo inicio de Corona
subconductores Estructura relativa subconductor
fase- corona kW/km-
por haz n del aire (cm)
fase (kV) (kV) fase
4 254 0.707 1.5325 371.231 365.892
4 229 0.637 1.5325 371.231 75.174
4 56 0.839 1.5325 371.231 404.916
3 254 0.707 1.908 371.231 322.950
3 229 0.637 1.908 371.231 66.352
3 56 0.839 1.908 371.231 357.394
En la Tabla 3.60 se presentan los valores de las pérdidas por efecto corona para cada una
de las torres analizadas, se realiza una comparación tanto para 4 como para 3
subconductores por haz. Obteniendo como resultado que con 3 subconductores por haz
se tiene valores de pérdida por efecto corona menores comparándolos con los valores con
4 subconductores por haz.

Los conductores de gran diámetro tienen gradientes más bajos del campo eléctrico en la
superficie del conductor por lo que el efecto corona es más bajo en relación con
conductores con menor diámetro. Al considerar las estructuras ubicadas en diferentes
altitudes sobre el nivel del mar se puede observar cómo varían estas pérdidas. Existe una
relación directamente proporcional con la altura, es decir a mayor altura sobre el nivel del
mar, mayor es el valor de pérdidas por efecto corona. El punto más crítico de estudio es la
estructura 229, que se encuentra a 4252.4 msnm el cual presenta valores de pérdida por
efecto corona mayores en relación con los demás, sin embargo, al reducir el número de
subconductores en el haz y aumentar el calibre de estos, se logra reducir dicho valor en
133.759 kW/km-fase.

3.5 Configuración óptima en el software PLS-CADD

En esta sección se presenta la configuración óptima del diseño de línea de transmisión de

500 kV para diseños futuros. Para esto se utiliza PLS-CADD el cual es un programa muy
potente y completo, utilizado para el diseño estructural y geométrico de líneas aéreas.
Integra uniformemente, en un solo medio, todos los aspectos de diseño de líneas de
transmisión, incluyendo modelado y representación del terreno, selección de rutas,
ubicación de estructuras a un mínimo costo manual o automático, flechas y tensiones,
espacios disponibles y esfuerzos, cálculos del campo eléctrico y magnético, generación de
lista de materiales, planos de planta y perfil y mucho más [45].

Para una línea de transmisión de extra alto voltaje, se considera el diseño óptimo con
configuración vertical, simple circuito, 3 subconductores por haz, ACAR 1700 MCM 42/19,
como se presenta en la Figura 3.46.
 Figura 3.46. Configuración vertical, simple circuito

Con el uso de PLS-CADD se simula el diseño resultante de la línea de trasmisión. En base

al archivo existente de la configuración actual se modifica la disposición geométrica de las
fases, la tensión de los conductores, así como la altura de las torres. De la misma manera
se cambian los subconductores por haz de ACAR 1100 MCM 18/19 a ACAR 1700 MCM
42/19. El modelo propuesto debe cumplir la distancia de seguridad del suelo al conductor
de al menos 11.5 metros.

Figura 3.47. Estructura GB103, disposición vertical

En la Figura 3.47 se observa el perfil de la estructura GB103, con una altura de la torre de
42 metros. Esta torre se encuentra a 1345.60 msnm. La estructura cuenta con un cable de
guarda de color verde y un cable óptico OPGW de color azul. De color rojo se muestran
cada una de las fases. De la misma manera se observa que se respeta la distancia de
seguridad mínima establecida desde el suelo hacia el conductor.

Figura 3.48. Estructura GB353N, disposición vertical

En la Figura 3.48 se observa el perfil de la estructura GB353N, con una altura de la torre
de 52 metros. Esta torre se encuentra a 4252.40 msnm. La estructura cuenta con un cable
de guarda de color verde y un cable óptico OPGW de color azul. De color rojo se muestran
cada una de las fases. De la misma manera se observa que se respeta la distancia de
seguridad mínima establecida desde el suelo hacia el conductor.
 Figura 3.49. Estructura GB378N, disposición vertical

En la Figura 3.49 se observa el perfil de la estructura GB378N, con una altura de la torre
de 40 metros. Esta torre se encuentra a 2976.40 msnm. La estructura cuenta con un cable
de guarda de color verde y un cable óptico OPGW de color azul. De color rojo se muestran
cada una de las fases. De la misma manera se observa que se respeta la distancia de
seguridad mínima establecida desde el suelo hacia el conductor.

En el Anexo I se muestra el tendido de los conductores considerando cierto número de

estructuras adyacentes a las torres estudiadas.
 4. CONCLUSIONES
Si bien es cierto, las simulaciones permiten obtener valores aproximados a la realidad, la
medición en el campo ayudó a contrastar los valores simulados y medidos. En este caso,
se utilizó una tolerancia de hasta el 7% ya que en las simulaciones el parámetro ajustable
es la permitividad y permeabilidad del aire y suelo que son valores aproximados a la unidad
a condiciones estándar de presión y temperatura. Por otro lado, al medir en el lugar donde
se encuentran las líneas de transmisión se tienen diferentes temperaturas, presión del aire,
densidad del aire, que tienen influencia directa en los campos electromagnéticos.

Factores como la altitud sobre el nivel del mar, la configuración geométrica de los
conductores, la distancia entre fases, el número de subconductores por haz, la altura desde
el nivel del suelo al conductor produce una variación en los valores del campo eléctrico y
magnético bajo una línea de transmisión. Un adecuado diseño involucra tomar en cuenta
estos aspectos de tal manera que los valores de los campos eléctricos y magnéticos se
encuentren dentro de los límites de exposición establecidos en las normas relacionadas.

En las estructuras 56 y 229 los valores de campos eléctricos y magnéticos a 1 metro sobre
el nivel del suelo se encuentran dentro de los valores establecidos en la norma IEEE Std.
644-1994. Sin embargo, para la estructura 254 el valor medido es de 4800 V/m el cual
excede el valor de campo eléctrico dado por la norma que es de 4167 V/m. En este caso,
dicho valor está fuera de los límites puesto que la altura del punto de amarre APA de los
conductores es de 24 metros sobre el nivel del suelo, mientras que el APA de las
estructuras 56 y 229 son de 35 y 37 metros, respectivamente. Por otro lado, el campo
magnético en las tres estructuras analizadas se encuentra dentro del límite establecido por
la norma que es de 67 A/m.

Al medir la diferencia de potencial existente entre los conductores y el suelo, se obtiene

que, en la estructura 56 que se encuentra a menor altitud sobre el nivel del mar, el campo
eléctrico es menor con un valor de 218.53 V/m, mientras que, la estructura 229 que se
encuentra a mayor altitud tiene un valor de campo eléctrico de 1373 V/m, comparados a
10 metros del centro de la torre, transversalmente, lado izquierdo. Esto viene relacionado
con la permitividad del aire que es inversamente proporcional a la altitud sobre el nivel del
mar, mientras que es directamente proporcional a la temperatura y a la presión La
densidad del aire aumenta con la presión y disminuye con la temperatura. Así mismo, la
densidad del aire disminuye con la altitud. La densidad relativa del aire que rodea la
estructura 56 es de 0.839 mientras que la densidad relativa del aire que rodea la estructura
229 es de 0.637.
Las diferentes configuraciones simuladas de la línea de transmisión demuestran que, para
una reducción en el valor del campo eléctrico y magnético, la disposición geométrica más
adecuada es el arreglo de fases de manera vertical. El campo magnético al depender
principalmente de la corriente que circula por el conductor no es un factor que decide al
momento de seleccionar una configuración adecuada; caso contrario ocurre con el campo
eléctrico que al depender del voltaje de la línea de transmisión puede reducirse según la
disposición geométrica de fases utilizada.

Al variar 4 subconductores en el haz con calibre ACAR 1000 MCM, 1100 MCM, 1200 MCM
con la configuración horizontal se tiene que, para el calibre ACAR 1100 MCM, los valores
del campo eléctrico son más bajos; mientras que para el campo magnético el calibre 1000
MCM es el que menor valor del campo magnético presenta.

Cambiando el número de subconductores del haz en una configuración horizontal con

calibre 1100 MCM se tiene que a menor número de subconductores el campo eléctrico es
menor. En este caso, para n = 3 el campo eléctrico disminuye en 5.8%, mientras que para
n = 5 el campo eléctrico aumenta en 10%, comparados con n = 4. Lo mismo ocurre para el
campo magnético, con 3 subconductores en el haz el campo magnético disminuye,
mientras que para n = 5 el campo magnético aumenta. Sin embargo, las variaciones
respecto a n = 4 no sobrepasan el 1%.

Las líneas más compactas presentan un menor campo eléctrico a nivel del suelo. En la
configuración horizontal si la distancia entre fases aumenta, el campo eléctrico aumenta.
Por otro lado, si la distancia entre las fases aumenta, el campo magnético disminuye bajo
las fases mientras que aumenta en los extremos de la franja de servidumbre.

Al variar el número de subconductores en el haz considerando diferentes calibres de

conductores según la corriente máxima que pasa por cada uno de ellos, se tiene que
mientras menor sea el número de subconductores en el haz el campo eléctrico es menor.
Por otro lado, a mayor número de subconductores en el haz el campo magnético es menor.
Para este caso, las simulaciones fueron realizadas para 3, 4 y 5 subconductores en el haz
con calibres ACAR 800 MCM, 1100 MCM y 1700 MCM, respectivamente.

El efecto corona provoca pérdidas en las líneas de transmisión, por lo que al momento de
diseñarlas se deben considerar todos los parámetros necesarios para que el efecto corona
sea mínimo, la configuración propuesta de 3 subconductores por haz con calibre ACAR
1700 MCM, al disminuir los valores de campos electromagnéticos, disminuye también las
pérdidas por efecto corona en la línea, reduciendo en 133.759 kW/km-fase comparada con
la configuración actual.
 5. RECOMENDACIONES
Para futuros diseños de líneas de transmisión de extra alto voltaje se sugiere implementar
una configuración vertical con 3 subconductores por haz, calibre ACAR 1700 MCM 42/19
que permite obtener los valores de campo electromagnético más bajos.

Para realizar una correcta medición se deben seguir los procedimientos establecidos en la
norma IEEE Std. 644-1994. En este caso, para el campo eléctrico el equipo empleado para
la medición debe estar a un metro sobre el nivel del suelo y al menos 2.5 metros distantes
de la persona que realiza la medición, de esta manera el valor del campo eléctrico no se
ve perturbado por la presencia de alguna carga. Por otro lado, el campo magnético al
depender netamente de la corriente que circula por el conductor no se ve afectada por la
presencia o cercanía del suelo, personas u objeto; no obstante, ningún material
ferromagnético debe estar cerca al instrumento de medida porque podría causar errores
en la medición.

La operación de una línea de transmisión de alto voltaje provoca interferencias

electromagnéticas. Esto es, existe una inducción de voltaje y corriente en superficies
conductoras. A consecuencia de esto, puede existir corrosión y deterioro en tuberías
metálicas que comparten el mismo derecho de vía de la línea de transmisión. Por esta
razón, para futuros diseños se recomienda realizar un estudio previo para que voltajes y
corrientes inducidas estén dentro de los límites establecidos en normas nacionales e
internacionales y de no ser el caso reducirlos para mitigar la afectación que ocasionan.

Como trabajo complementario al presente se sugiere realizar un análisis técnico

económico de la implementación de la configuración óptima propuesta, ya que al tener un
calibre de subconductor mayor implica más peso del conductor lo que se traduce,
probablemente, en estructuras más robustas y altas.
6. REFERENCIAS BIBLIOGRÁFICAS
 7. ANEXOS

ANEXO I. PLS-CADD

ANEXO II. Registro Fotográfico Medición In situ

ANEXO III. Los modelos simulados en FEMM se encuentran en la siguiente dirección web:
https://github.com/KDPJPedro/FEMM.git
 ANEXO I

GB103 (As Built: 56)

GB353N (As Built: 229)
GB378N (As Built: 254)
 ANEXO II

Estructura 229.
Estructura 254.
Estructura 56.
ORDEN DE EMPASTADO