INSTITUTO POLITECNICO NACIONAL

ESCUELA SUPERIOR DE INGENIERIA MECÁNICA Y ELÉCTRICA

UNIDAD PROFESIONAL “ADOLFO LÓPEZ MATEOS ZACATENCO
DEPARTAMENTO DE INGENIERIA ELÉCTRICA

ANÁLISIS COMPARATIVO DE INTERFERENCIA

ELECTROMAGNÉTICA EN LINEAS DE TRANSMISIÓN
CONVENCIONALES Y DE ALTA CARGABILIDAD

TESIS

QUE PARA OBTENER EL TITULO DE

INGENIERO ELECTRICISTA

PRESENTAN
C. DAVID ORTIZ CALLEJAS
C. ERICK BRAYN FERRER ABILES
C. JOSÉ IVÁN CRUZ VILLA

ASESOR TÉCNICO
M. EN C. CARLOS TEJADA MARTÍNEZ

ASESOR METODOLÓGICO
LIC. MARÍA DE LOURDES GARCÍA PIÑÓN

CDMX, 2024
ESPACIOS PARA ACTA DE APROBACIÓN
ESPACIO PARA CESIÓN DE DERECHOS (SI ES EL CASO) DEL TRABAJO
 Agradecimientos
En el trayecto académico que he recorrido en la Escuela Superior de Ingeniería
Mecánica y Eléctrica, Unidad Zacatenco, me siento dichoso y orgullosos por la
compañía de todos los involucrados en este proceso, momentos y emociones, por lo
que estoy agradecido
Primeramente, agradezco a mis padres el señor C. Diego Ferrer y la señora C. Tomasa
Abiles, así como a mis hermanos C. Alan y C. Cristian, por el apoyo que me han
brindado en diferentes aspectos de mi vida.
Al profesor M. en C. Carlos Tejada Martínez por su apoyo y gran esmero que dieron
como resultado el presente trabajo.
De igual manera, agradezco a la profesora Lic. María de Lourdes García Piñón por el
tiempo dedicado, la atención y apoyo en la asesoría del presente trabajo.
Así como a todos los profesores que aportaron conocimientos, atención y ayuda en
cada etapa de mi trayectoria escolar.
También a mis compañeros y amigos a los cuales estimo y agradezco por todas las
experiencias vividas.
Sin olvidar al Instituto Politécnico Nacional y a la Escuela Superior de Ingeniería
Mecánica y Eléctrica, Unidad Zacatenco, por haberme permitido crecer en mi
desarrollo profesional, brindándome sus instalaciones, así como el apoyo del personal
docente.
Gracias
Erick Brayn Ferrer Abiles
 Agradecimientos
Quiero expresar mi más sincero agradecimiento a todas las personas que han formado
parte de este proceso, por contribuir en mi desarrollo personal y profesional.
Agradezco a mis padres, el señor C. Gerardo Cruz Gómez y la señora C. Verónica
Villa Caporal, por su apoyo incondicional, amor y motivación constante, son las
personas a quienes les debo todo. Su sacrificio y dedicación han sido fundamentales
en el desarrollo de mi persona.
A mis hermanos el C. Gerardo y la C. Nadia, por su compañía y apoyo durante todo
este proceso, su presencia, consejos y guía siempre han sido una fuente de motivación
para alcanzar mis objetivos.
Al M. en C. Carlos Tejada Martínez, y a la Lic. María Lourdes García Piñón, por su
apoyo, colaboración y dedicación para la realización del presente trabajo.
A todos los profesores que tuve la fortuna de conocer en el Instituto politécnico
nacional, por sus valiosas lecciones, conocimiento y vocación docente. Siendo como
resultado en mi persona una buena formación académica.
Finalmente agradezco al Instituto Politécnico Nacional y a la Escuela Superior de
Ingeniería Mecánica y Eléctrica Unidad Zacatenco, por permitirme ser parte de esta
comunidad, que me brindo en forma grata y no tan grata memorables experiencias,
por brindarme educación de excelencia y por el privilegio de poder formarme
profesionalmente en la institución.

Gracias
José Iván Cruz Villa.
 Agradecimientos
A lo largo de mi formación académica en la Escuela Superior de Ingeniería Mecánica
y Eléctrica, Unidad Zacatenco, me considero muy afortunado por haber contado con
el apoyo de diversas personas, apoyo que, de alguna forma, se ve reflejado en este
trabajo. Quiero expresar mi más profundo agradecimiento a:
Mi familia, mis padres el señor C. Rufino Ortiz López y la señora C. María Victoria
Callejas Martínez, mis hermanas C. Karina Guadalupe y C. Alondra Jocelyn, quienes
con su apoyo incondicional y motivación constante me acompañaron en cada etapa
de este camino. A cada integrante de mi familia en general, quienes siempre me
brindaron su confianza y ánimo en los momentos difíciles.
Agradezco de manera especial al profesor M. en C. Carlos Tejada Martínez y a la
profesora Lic. María de Lourdes García Piñón, por su esencial asesoría y orientación,
así como a todos los profesores que, con su dedicación y profesionalismo,
contribuyeron al desarrollo de este trabajo.
No me olvido de mis compañeros de la escuela y amigos, quienes hicieron de esta
experiencia algo enriquecedor y llevadero, ofreciendo siempre una mano amiga.
También, expreso mi gratitud al Instituto Politécnico Nacional y a la Escuela Superior
de Ingeniería Mecánica y Eléctrica, Unidad Zacatenco, por su compromiso con la
excelencia académica y por darme la oportunidad de desarrollar mis habilidades
profesionales y académicas en un entorno de alto nivel.
Agradezco a cada uno de ustedes
David Ortiz Callejas
 ÍNDICE
Resumen .................................................................................................................................................... I
Introducción ............................................................................................................................................. II
Antecedentes .......................................................................................................................................... III
Planteamiento del problema................................................................................................................... IV
Objetivo general ....................................................................................................................................... V
Objetivos específicos ................................................................................................................................ V
Justificación ............................................................................................................................................. VI
Alcances y limitaciones........................................................................................................................... VII
Capítulo I Marco Teórico ......................................................................................................................... 8
Elementos de una línea de transmisión de alta tensión y carga habilidad. ............................................. 9
1.1 Líneas de transmisión ............................................................................................................... 9
1.2 Estructuras (torres) y tipos de estructuras ............................................................................... 9
1.2.1 Estructuras auto soportadas con celosía ........................................................................ 10
1.2.2 Estructuras auto soportadas tipo tubular ...................................................................... 10
1.2.3 Estructuras con retenida ................................................................................................ 11
1.3 Conductores para líneas de transmisión ................................................................................ 11
1.4 Conductores en haz ................................................................................................................ 14
1.5 Aisladores ............................................................................................................................... 15
1.6 Herrajes .................................................................................................................................. 18
1.7 Cargabilidad de una línea ....................................................................................................... 19
1.7.1 Líneas de transmisión de alta cargabilidad .................................................................... 20
1.7.2 Capacidad técnica del conductor ................................................................................... 23
1.7.3 Caída de tensión ............................................................................................................. 23
1.7.4 Margen de estabilidad .................................................................................................... 24
1.8 Parámetros de la línea ............................................................................................................ 25
1.8.1 Resistencia ...................................................................................................................... 26
1.8.2 Inductancia ..................................................................................................................... 27
1.8.3 Capacitancia ................................................................................................................... 28
1.8.4 Conductancia .................................................................................................................. 29
1.9 Derecho de vía........................................................................................................................ 30
Capítulo II Marco Metodológico............................................................................................................ 32
Métodos de cálculo de interferencia electromagnética en líneas de trasmisión de alta tensión.......... 33
 2.1 Efecto corona.......................................................................................................................... 33
2.1.1 Naturaleza y manifestación del efecto corona. .............................................................. 34
2.1.2 Descripción física del efecto corona ............................................................................... 34
2.2 Ecuaciones empíricas y semiempiricas .................................................................................. 35
2.2.1 Métodos empíricos......................................................................................................... 36
2.3 Interferencia Electromagnética. ............................................................................................. 36
2.4 Gradiente Crítico .................................................................................................................... 36
2.5 Cálculo del campo eléctrico en la superficie de los conductores. .......................................... 37
2.5.1 Método de los coeficientes de potenciales de Maxwell. ............................................... 38
2.5.2 Método de Mark y Mengele........................................................................................... 40
2.6 Ecuaciones para el Cálculo de la Radio Interferencia. ............................................................ 42
2.6.1 CIGRE .............................................................................................................................. 42
2.6.2 BPA.................................................................................................................................. 44
2.6.3 CREIPI ............................................................................................................................. 47
2.7 Límites máximos de RI en líneas de transmisión.................................................................... 49
Capítulo III Marco Técnico .................................................................................................................... 51
3.1 Simulación del campo eléctrico en COMSOL MULTIPHYSICS ................................................. 53
3.1.1 Caso 1: Línea de transmisión convencional con 4 subconductores por haz, operando a
500 kV. 55
3.1.2 Caso 2: Línea de transmisión de alta cargabilidad (HSIL) convencional con 6
subconductores por haz, operando a 500 kV. ................................................................................ 72
3.1.3 Caso 3: Línea de transmisión de alta cargabilidad (HSIL) no convencional con 4
subconductores por haz, operando a 500 kV. ................................................................................ 75
3.1.4 Caso 4: Línea de transmisión de alta cargabilidad (HSIL) no convencional con 6
subconductores por haz, operando a 500 kV. ................................................................................ 77
3.2 Cálculo de Radio Interferencia ............................................................................................... 79
3.2.1 CIGREE ............................................................................................................................ 79
3.2.2 BPA.................................................................................................................................. 81
3.2.3 CREIPI ............................................................................................................................. 83
Capítulo IV ............................................................................................................................................. 85
4.1 Análisis económico ................................................................................................................. 86
4.1.1 Costos unitarios .............................................................................................................. 88
4.2 Conclusiones........................................................................................................................... 90
4.3 Recomendaciones .................................................................................................................. 91
 4.4 Referencias. ............................................................................................................................ 92
4.5 Glosario .................................................................................................................................. 94
Índice de figuras ....................................................................................................................................... 1
Índice de tablas ........................................................................................................................................ 3
Apéndices ................................................................................................................................................. 4
 Resumen

Debido constante incremento de la demanda de energía eléctrica en México, se ha

impulsado la necesidad de aumentar la capacidad de los sistemas de transmisión. Esta
necesidad ha llevado a la utilización de tensiones de transmisión mayores, lo que
origina el fenómeno conocido como efecto corona en las líneas de transmisión (LT) de
corriente alterna (CA). Dicho fenómeno genera efectos significativos, tales como ruido
audible, pérdidas de potencia, distorsión de ondas viajantes e interferencia
electromagnética (RI). Las emisiones electromagnéticas resultantes del efecto corona
abarcan un amplio intervalo de frecuencias, lo que puede perturbar la operación
habitual de dispositivos electromagnéticos cercanos a las líneas.
Uno de los métodos para aumentar la capacidad de transmisión es la conversión de
LT de CA a líneas de alta cargabilidad (HSIL, por sus siglas en inglés). Esta técnica
optimiza la configuración geométrica de los conductores para disminuir la impedancia
característica de la línea, incrementando así su potencia natural. Entre las
configuraciones propuestas, se encuentran el aumento del número de conductores por
fase y la compactación de la línea, lo que implica reducir la distancia entre fases y su
altura. Sin embargo, aunque las líneas HSIL incrementan significativamente la
potencia transmitida en comparación con las líneas convencionales, también tienden
a elevar el campo eléctrico sobre la superficie de los conductores, lo que propicia la
aparición del efecto corona y, por ende, un aumento en los niveles de RI.
Este trabajo presenta un análisis comparativo de los perfiles de RI en cuatro casos de
estudio: uno correspondiente a una línea de transmisión convencional y tres a líneas
HSIL, incluyendo una convencional y dos no convencionales. Se seleccionan diversas
ecuaciones para el cálculo de RI, abarcando condiciones de buen ambiente, así como
ambientes húmedos o climas lluviosos. Estas ecuaciones se implementan en un
programa de cómputo (MATLAB®) utilizando datos de líneas de transmisión reales,
junto con el campo eléctrico sobre superficie de los conductores simulado en COMSOL
Multhipisics. Posteriormente, los resultados se comparan con los límites máximos
establecidos por norma, evidenciando que los tres casos de HSIL presentan niveles
de RI que superan lo recomendado.
Finalmente, se ajustan las alturas de las fases en los diseños HSIL con el objetivo de
reducir los valores máximos de RI, al mismo tiempo que se busca mantener la máxima
capacidad de transmisión posible. Sin embargo, este proceso de optimización no se
logra de manera satisfactoria en # de los casos analizados.

I
 Introducción

El incremento en la demanda de electricidad ha fomentado el desarrollo de sistemas

de transmisión más efectivos y el incremento en la capacidad. Este incremento ha
llevado al uso de niveles de tensión más altos en las líneas de transmisión (LT) de
corriente alterna (CA), originando fenómenos no deseados como el efecto corona. El
efecto corona genera interferencia electromagnética (RI), que influye en la operación
de dispositivos cercanos y comprometer la calidad del servicio eléctrico. Estas
emisiones abarcan un gran intervalo de frecuencias y se ven intensificadas en
situaciones climáticas desfavorables, como la humedad o la lluvia.
Para afrontar a las limitaciones de las LT convencionales, se han propuesto las líneas
de alta cargabilidad (HSIL), las cuales incrementan la capacidad de transmisión
mediante la compactación de la configuración geométrica de los conductores. Si bien
las líneas HSIL mejoran la capacidad de transmisión, también tienden a elevar el
campo eléctrico en la superficie, lo que resulta en mayores niveles de RI en
comparación con las líneas convencionales. Esta situación genera dudas sobre la
factibilidad de las líneas HSIL en términos de RI, especialmente bajo condiciones
ambientales desfavorables.
El presente trabajo tiene como objetivo llevar a cabo un análisis comparativo de la RI
en LT convencionales y en líneas HSIL. A través de la implementación de diversas
ecuaciones en MATLAB® y utilizando datos de líneas reales, junto con el campo
eléctrico sobre superficie de los conductores simulado en COMSOL Multhipisics, se
calculan los perfiles de RI en ambos tipos de líneas. Se analizan cuatro casos de
estudio, examinando si las líneas HSIL presentan niveles de RI que rebasen los límites
normativos establecidos y contrastando su comportamiento con el de las líneas
convencionales.

II
 Antecedentes
El famoso científico Claude H. Gary contribuyó de manera importante a la comprensión
de la interferencia electromagnética causada por corona en las líneas de transmisión.
En un artículo publicado en 1971, propuso un método para estimar los efectos de la
interferencia electromagnética.
En 1988, L. V. Timashova llevó a cabo mediciones de interferencia electromagnética
en transmisiones de 750 kV a una frecuencia de 500 kHz, evaluando diversas
condiciones climáticas. Sus hallazgos mostraron que en condiciones de lluvia y nieve
los niveles de interferencia electromagnética son más elevados.
En el año 2001, Allah y Alghamdi investigaron el efecto corona en una línea de CA de
circuito doble de 220 kV, la cual estaba siendo transformada en una línea de CD de
220 kV. Se evaluaron cuatro variantes para la conversión del circuito y así determinar
cuál era la más adecuada para minimizar la interferencia electromagnética y el ruido
audible.
Los conceptos sobre cargabilidad y comportamiento de líneas de transmisión, se
efectuaron con base a las tensiones de las líneas que entonces estaban en operación.
Entonces el estudio de cargabilidad de líneas de transmisión fue considerado por H.P.
St Clair donde, con base a dichos niveles de tensión y consideraciones prácticas sobre
las restricciones de límite térmico, caída de tensión y límite de estabilidad, así como la
evaluación de la cargabilidad para diferentes longitudes de línea. Asimismo, con base
en la experiencia sobre el comportamiento de las líneas se estableció un criterio
práctico que citaba una referencia de 1 pu. (Por unidad) de SIL (Carga de Impedancia
de Sobretensiòn) para una línea convencional de 330 kV, 300 millas, que opere a 60
Hertz; de tal manera que para esta carga la pérdida neta de potencia reactiva es igual
a cero [12].

III
 Planteamiento del problema

La interferencia electromagnética (RI) originada por el efecto corona en líneas de

transmisión (LT) de alta tensión (AT) implica una dificultad técnica. Este fenómeno,
presente en líneas de corriente alterna (CA), genera emisión de ruido electromagnético
que puede perturbar la operación habitual de dispositivos electrónicos cercanos e
impactar en la calidad del suministro eléctrico. Conforme aumenta la demanda de
energía, las LT convencionales enfrentan limitaciones en su capacidad para
transportar mayor energía, lo que ha llevado a la implementación de líneas de alta
cargabilidad (HSIL) como una opción. [3].
Las líneas HSIL, diseñadas para elevar la capacidad de transmisión a través de la
optimización geométrica de los conductores, suelen producir mayores niveles de RI
debido al incremento en el campo eléctrico superficial. Este aumento en las emisiones
electromagnéticas suscita cuestionamientos sobre el impacto comparativo entre las
líneas convencionales y las HSIL para procurar que los niveles de RI se mantengan
dentro de los límites establecidos por norma., especialmente en condiciones climáticas
adversas como la lluvia o alta humedad.

IV
 Objetivo general
Analizar los niveles de interferencia electromagnética en líneas de transmisión de alta
tensión convencionales y de alta cargabilidad para su comparación con los niveles
establecidos por norma

Objetivos específicos
• Describir los elementos que conforman las líneas de transmisión de alta
cargabilidad.
• Emplear métodos de cálculo de interferencia electromagnética en líneas de alta
tensión mediante los programas MATLAB Y COMSOL.
• Comparar los resultados de interferencia electromagnética en líneas de
transmisión de alta tensión convencionales de corriente alterna y de alta
cargabilidad.

V
 Justificación

Las líneas de transmisión de alta tensión de corriente alterna, híbridas y de alta

cargabilidad producen fenómenos electromagnéticos que pueden ser causantes sobre
la funcionalidad en las propias líneas de transmisión, afectando inclusive las
características de su capacidad, sin embargo, ante una creciente demanda de energía
eléctrica, el optar por un cambio de líneas de corriente alterna (CA) a corriente directa
(HVDC) resulta una opción limitante debido a los obstáculos y una inversión económica
mayor. El efecto corona es un fenómeno que se presenta en las líneas de transmisión
cuando el campo eléctrico en la superficie de los conductores sobrepasa la rigidez
dieléctrica del aire. Una de las consecuencias de este fenómeno es la interferencia
electromagnética, también llamada Radio Interferencia (RI). Existen criterios que
catalogan los niveles de RI, se llevará a cabo el análisis con los softwares MATLAB y
COMSOL que se explican brevemente; Software MATLAB, es esencial para realizar
análisis numéricos avanzadas y simulaciones complejas, ofreciendo herramientas
robustas para procesamiento de señales, análisis de datos y modelado matemático.
La flexibilidad y su basta biblioteca de funciones permite desarrollar modelos precisos
y algoritmos personalizados, asegurando la integridad y fiabilidad de las líneas de
transmisión. Por otro lado, COMSOL es otro software es esencial para simular la
interferencia electromagnética en líneas de transmisión de alta tensión debido a su
capacidad para integrar múltiples fenómenos físicos, el cual ofrece un análisis
detallado en dominios de frecuencia y tiempo, el cual permite la personalización
precisa de modelos complejos. Esto facilitará la identificación y la mitigación de efectos
adversos de la interferencia electromagnética (EMI), optimizando el diseño y
asegurando la fiabilidad de las líneas de transmisión.
Con el fin de estudiar una reducción en cuanto a estos niveles de radio interferencia
(RI) que puedan presentarse en diferentes casos de estudio. En cuanto a las líneas de
transmisión de alta tensión, se tiene un gran empleo de las llamadas líneas
convencionales, sin embargo, están latentes las líneas de transmisión de alta
cargabilidad, donde se busca la innovación y optimización en los valores SIL (Carga
de Impedancia de Sobretensiòn), es decir una compensación en los niveles de
potencia de las líneas de transmisión.

VI
 Alcances y limitaciones
Alcances
El presente trabajo analiza los niveles de interferencia electromagnética en líneas de
transmisión de alta tensión convencionales y de alta cargabilidad, comparando estos
niveles con los niveles establecidos por las normas vigentes. En primer lugar, se
describirán los componentes y características técnicas de las líneas de transmisión de
alta tensión, tanto convencionales como de alta cargabilidad. Este análisis se centrará
específicamente en cuatro tipos de líneas que son: una línea de transmisión de doble
circuito AC/AC de 400 kV, una línea de transmisión de circuito híbrido AC/DC de 400
kV/±500 kV, una línea de transmisión de doble circuito DC/DC de ±500 kV y una línea
de transmisión HSIL de 500 kV.

Se aplicarán diversos métodos de cálculo de interferencia electromagnética en líneas

de transmisión de alta tensión, se consideran cuatro estudios de caso detallados para
cada uno de los tipos de líneas de transmisión, donde se evaluarán los niveles de
interferencia electromagnética y comparándolos con los límites normativos
establecidos. Además, los resultados obtenidos a través de cálculos y simulaciones
serán validados mediante datos empíricos y experimentales disponibles, asegurando
así su consistencia y precisión.

Limitaciones
Una de las limitaciones es el estudio que se centrará en líneas de transmisión situadas
en México. La disponibilidad de datos específicos puede estar limitada por
restricciones de acceso, debido a que las líneas de transmisión de alta cargabilidad no
se encuentran en el país.

VII
Capítulo I
Marco
Teórico

8
Elementos de una línea de transmisión de alta tensión y carga habilidad.
En este capítulo se presenta una breve investigación de los elementos que
conforman las líneas de transmisión, como lo son las torres, conductores, aisladores
etc. Ya que estos repercuten en los niveles de radio interferencia que presenta el
propio sistema de transmisión. También se abordan las características y geometría
que pueden tomar los conductores en haz, así como los valores promedio máximo
y mínimo del gradiente eléctrico superficial.
1.1 Líneas de transmisión
Las líneas de transmisión son aquellos elementos que conforman al sistema
eléctrico de potencia, se encargan de transportar o transmitir grandes cantidades
de energía eléctrica a partir de estaciones de generación eléctrica hasta puntos o
estaciones receptoras para su pertinente distribución y utilización [2][3].
Las líneas de transmisión son identificadas en dos clases, las líneas aéreas y las
líneas subterráneas; en el caso de las líneas de transmisión aéreas, son aquellas
que se encuentran sujetadas por las torres de transmisión eléctrica, las cuales
soportan a los conductores principalmente, siendo el aire el aislamiento que hay
entre conductores, mientras que entre los conductores y la tierra se emplean
cadenas de aisladores, de tal modo de evitar el contacto entre conductores y las
estructuras de soporte [1][2].
Por otro lado, las líneas de transmisión subterráneas son aquellas se colocan en
zanjas en el terreno, estas líneas surgen por dificultades que en ocasiones se
presentan en las líneas aéreas, los conductores suelen ser cubiertos por hormigón
armado, por ende deben de considerarse todos los aspectos de la zona de
instalación de dichas líneas debido a factores como corrosión o incluso deterioro
mecánico por afecciones de la industria cercana, así como marcar y establecer
puntos de acceso al sistema de líneas de transmisión para ocasiones donde se
requiera accesar [6].
1.2 Estructuras (torres) y tipos de estructuras
Las estructuras en un sistema de líneas de transmisión son esenciales, ya que son
el sostén mecánico de cada elemento que constituyen a dichas líneas, comúnmente
metálicos; al ser fundamentales para la propia transmisión de energía, se asumen
como la mayor inversión económica para la implementación de las líneas de
transmisión [1].
Las proporciones para la instalación de las torres de transmisión dependen de
ciertos criterios a considerar, como lo son:

• Tipo de superficie: plano, irregular u ondulatorio, escarpado.

9
 • Distancia interpostal o claro horizontal: refiriéndose a la distancia que hay
entre dos puntos separados, en este caso, entre dos torres de transmisión
adyacentes.
• Finalidad de la torre: De acuerdo a la finalidad con la que se busca
implementar una torre, se pueden derivar de acuerdo a su tipo de
suspensión, tipo de tensión, tipo remates y de transposición.
Estableciendo los criterios anteriores, se pueden mencionar tres tipos de estructuras
para las líneas de transmisión [1].
1.2.1 Estructuras auto soportadas con celosía
Mayormente utilizadas en México para las líneas de transmisión debido a que no
son dependientes ni necesitan de demás elementos mecánicos para mantenerse
colocadas, es por ello que se deriva su nombre, siendo autosuficientes con sus
propias dimensiones, Figura 1 [1].

Figura 1. Estructura autosoportada [17].

1.2.2 Estructuras auto soportadas tipo tubular
De igual manera que las torres autosoportadas de celosía, estas no dependen
meramente de mecanismos externos para mantener su posición y brindar soporte,
están hechas con tubo de acero (Figura 2), siendo esta característica lo que define
su estructura más compacta, sin embargo, su costo es más elevado ya que pueden
ser diseñadas para soportar la tensión o suspensión, por lo que su implementación
provoca que solo sea en zonas específicas como lo son terrenos con dificultades
con la disponibilidad o por estética [1].

10
 Figura 2. Estructura autosoportada tipo tubular [17].
Las configuraciones más típicas de disposición de fases en el espacio son:

• Horizontal: Generalmente producen estructuras de menor peso.

• Vertical: Genera estructuras que demanda menor ancho de derecho de vía.
• Delta: Minimiza perdidas eléctricas y efectos de campo magnético para
estructuras de un solo circuito.
1.2.3 Estructuras con retenida
Este tipo de estructura funciona de acuerdo a una trabe horizontal que se encuentra
sostenida a uno o dos puntos, principalmente su estabilidad mecánica se ve
asegurada por tirantes (retenidas) adecuadas para sostenerse [1].
1.3 Conductores para líneas de transmisión
Los conductores eléctricos son básicamente el medio físico por donde se transporta
la energía eléctrica; de esta forma es posible transmitir la energía producida en
centrales generadoras a subestaciones para su distribución [2].
El material de construcción de un conductor eléctrico es cualquier sustancia que
permita de forma sencilla el flujo de electrones a través de su superficie, una vez le
se haya suministrado una diferencia de potencial en sus extremos. Los metales, son
excelentes conductores debido a su composición y alto nivel de conductividad, es
por ello su empleo en la fabricación de conductores para líneas de transmisión [2].
Los metales utilizados por excelencia en la fabricación de conductores eléctricos
son el aluminio y el cobre, no obstante, para conducir energía de forma general, en
la actualidad se pueden emplear materiales tales como aleaciones de cobre, plata,
hierro y acero [2].

11
El cobre y aluminio siendo uno de los materiales más utilizados en la fabricación de
conductores eléctricos, no son los mejores conductores, existen materiales con
mayor nivel de conductividad como lo es el caso de la plata, sin embargo, su alto
costo y proliferación en el hurtado de este metal hace inconveniente el uso del
mismo, ya que en la aplicación de las líneas de transmisión sería necesaria gran
cantidad de este material [2].
Por esta razón es que la selección del material conductor depende de diferentes
factores, no solo del aspecto económico o de las características eléctricas del
conductor. Si no también es importante considerar otros aspectos como lo son:
propiedades mecánicas, facilidad de hacer conexiones, su mantenimiento, la
cantidad de soportes necesarios, las limitaciones de espacio, resistencia a la
corrosión del material y otros [10].
En líneas de transmisión el material que se utiliza es el aluminio, ya que el costo de
este es menor en comparación con el del cobre, además de que es más ligero, lo
cual simplifica el trabajo de diseño de los soportes que sostienen los conductores,
esto refleja un costo menor ya que la mayor inversión que se realiza en los sistemas
de transmisión son precisamente en las torres de alta tensión. Debido a que el
aluminio presenta una baja rigidez mecánica, se recurre a la utilización de cable de
aluminio aleado o con refuerzos para su aplicación en líneas de transmisión.
Los conductores de aluminio más utilizados se mencionan en los siguientes incisos:
[2]

• AAC: All Aluminum Conductor (Conductor de Aluminio)

• AAAC: All Aluminum Alloy Conductor (Conductor de aluminio con aleación)
• ACSR: Aluminum Conductor Steel Reinforced (Conductor de Aluminio con
Refuerzo de Acero) 6
• ACAR: Aluminum Conductor Alloy Reinforced (Conductor de Aluminio con
Refuerzo de Aleación).
A continuación, se describen las características de los conductores antes
mencionados.
AAC – Conductor de aluminio: Cable constituido por alambres del mismo material.
Los manufacturados con aluminio relativamente puro son llamados All- Aluminum
Conductor (AAC). Posee enteramente las propiedades paralelas de los alambres
individuales. (ver Figura 3) ASTM B231, TS IEC 1089, DIN 48201, BS 215, UNE
21.018 [17]

12
 Figura 3. Conductor AAC [18].
AAAC – Conductor de aleación de aluminio. Los manufacturados con aleaciones de
aluminio son llamados All-Aluminum-Alloy Conductor (AAAC). Buena resistencia a
la corrosión. Tiene alta conductividad combinada con resistencia mecánica
moderada. (ver Figura 4) ASTM B399, ASTM B231, TS IEC 1089, DIN 48201, BS
215, UNE 21.018 [17]

Figura 4. Conductor AAAC [18].

ACSR - Conductor de aleación de aluminio reforzado con acero: Es el tipo más
común se compone por un núcleo de alambres de acero cubierto por una o varias
capas de alambres de aluminio. Posee las propiedades paralelas de los alambres
individuales. Buena resistencia a la corrosión. (ver Figura 5) ASTM B231, ASTM B
232, DIN 48201, BS 215, UNE 21.018. [17]

Figura 5. Conductor ACSR [18].

ACAR- Conductor de Aluminio con Refuerzo de Aleación: Es compuesto por la
combinación de alambres de aluminio relativamente puro y de aleaciones de
aluminio. Soporta alta resistencia mecánica, pero menor conductividad. (ver Figura
6) ASTM B230, ASTM B398 y ASTM B524 [17]

13
 Figura 6. Conductor ACAR [18].
En México, generalmente se utiliza conductor ACSR (Aluminum Conductor Steel
Reinforced) en las líneas de transmisión, este conductor posee en su interior alma
de acero y en su exterior se encuentran entrelazados los hilos conductores de
aluminio. El alma de acero proporciona soporte mecánico, por lo que este tipo de
conductor es ideal en longitudes largas de los claros entre torres. [1]
La selección del calibre del conductor dependerá de diversos factores, tales como
los que se muestran a continuación: [1]

• En primer lugar, debe tomarse en cuenta la capacidad de conducción de

corriente del conductor a una temperatura específica.
• La máxima caída de tensión que se puede presentar en la línea.
• Factores mecánicos como resistencia a la ruptura, coeficientes de dilatación
y peso.
• Las pérdidas que se presentan por efecto Joule, así como su límite
económico.

1.4 Conductores en haz

Anteriormente en los sistemas de transmisión de alta tensión, se consideró que el
campo eléctrico en los conductores era un factor importante en la generación del
efecto corona, además de que dicho campo decrece en proporción inversa al radio
del conducto, se optó por incrementar el diámetro de los conductores. No obstante,
esta alternativa no fue viable en términos económicos ya que la cantidad para la
construcción de los conductores era considerable. Una segunda opción fue utilizar
conductores huecos, pero las complicaciones en cuanto a la manufactura y el
mantenimiento de los mismos, hicieron de esta segunda alternativa un fracaso [10].
Hoy en día la solución que ha permitido transmitir tensiones cada vez más elevadas
fue la introducción de conductores en Haz. En 1930, Markt y Mengele obtuvieron la
patente para esa alternativa. Ellos argumentaron las ventajas de utilizar varios
conductores en cada una de las fases, tales como la disminución en los gradientes
de potencial que se presentan en la superficie de los conductores, así como
incrementar la capacidad de transmisión de la línea [10].

14
Actualmente se utilizan en la mayoría de los sistemas de transmisión, conductores
en Haz; de hecho, para niveles de tensión superiores a 420 kV no es factible utilizar
un solo conductor por fase [10].
En la Figura 7. se ilustra un haz de seis conductores. Las características principales
de un haz son el número de subconductores nc, su radio r, radio del haz R, radio
equivalente Req, los valores promedio, máximo y mínimo del gradiente eléctrico
superficial, Ea, Emáx, Emin respectivamente, y la separación entre subcontuctores
adyacentes que se calcula mediante la ecuación 1:
𝜋 (1)
𝑠 = 2𝑅𝑠𝑒𝑛 ( )
𝑛𝑐

Figura 7. Haz de seis conductores [3].

1.5 Aisladores
Los aisladores en líneas eléctricas de alta tensión son los elementos que cumplen
la función de sujetar mecánicamente a los conductores que forman parte de la línea,
manteniéndolos aislados de tierra y de otros conductores [2].
Teniendo en cuenta lo anterior los aisladores deben de mantener sus cualidades
mecánicas y eléctricas durante varios años después de haber sido instalados,
deben resistir ante las condiciones de todo tipo a los que serán sometidos.
La clasificación de los aisladores puede ser diversa dependiendo del punto de vista.
Es importante saber que existen en principio aisladores de corriente continua y de
corriente alterna.
Una forma de clasificarlos es según su uso, dentro de estos podemos encontrar:

• Aisladores de Intemperie.
• Aisladores de recinto cubierto.

15
 • Aisladores de suspensión.
• Aisladores de amarre.
• Aisladores de apoyo.
Otra forma de clasificarlos es de acuerdo con el material de construcción, por
ejemplo:

• Aisladores de vidrio,
• Porcelana.
• Plástico.
Al día de hoy los aisladores que han recibido la sanción favorable con base en la
experiencia son los del tipo cerámico y de vidrio (Figura 8). En la medida en la que
aumenta el nivel de tensión en las líneas de transmisión el material con el que se
diseñan y construyen los aisladores ha sido particularmente objeto de estudio, esto
se debe al gran esfuerzo eléctrico al que se ve sometido el aislador, por ello, en la
actualidad es necesario conocer el nivel de aislamiento de este elemento.

Figura 8. Aislador de Vidrio [2].

Por otro lado, el material que parece haber dado mejores resultados en su uso a la
intemperie es la porcelana. Por lo anterior se usa comúnmente la porcelana vidriada
la cual está constituida por una mescolanza de feldespato, cuarzo y caolín ya que
este material es de los mejores para satisfacer las condiciones requeridas de un
buen material aislante, sin embargo, este material posee algunas desventajas que
a continuación se presentan:

• El coeficiente de expansión térmica de la porcelana debe ser igual que el de

la capa de vidrio que la recubre, pues de lo contrario surgen tensiones
internas, que después de un cierto tiempo provocan daños mecánicos
manifestados en forma de grietas, acortando así la vida útil del elemento.
• Durante la segunda cocción de la porcelana el material experimenta una
contracción la cual afecta las dimensiones finales del elemento.

16
 • El material de porcelana es muy frágil durante descargas por arco, esto una
vez puesto en servicio el aislador.

Los materiales idóneos en la fabricación de aisladores en líneas de transmisión

aéreas son el vidrio y la porcelana (Figura 9). Por otro lado, los aisladores de plástico
son frecuentemente utilizados en instalaciones bajo techo, donde el grado de
contaminación y nivel de corrosión es menor.

Figura 9. Aislador de Porcelana [2].

En la Tabla 1 se muestran las características de los aisladores de vidrio y porcelana
[].
Tabla 1. Características de los aisladores de vidrio y porcelana [].

Vidrio Porcelana
Es sencilla la detección de defectos en Debido al material de porcelana no es
el elemento, debido a que el material es sencillo detectar defectos internos.
cristalino.
Se produce menor calentamiento ya Se produce un mayor calentamiento ya
que los rayos solares no son absorbidos que los rayos solares son absorbidos
por el aislador. por el material.
Fácil inspección desde el suelo o Es necesario realizar una inspección
helicóptero, se observa su estado sin minuciosa del elemento.
ningún error.
En situaciones de cambios de tensión, Se pueden producir perforaciones en el
no se producen perforaciones. material, debido a los esfuerzos
mecánicos.

17
En las líneas de transmisión eléctrica se manejan altas tensiones por lo que es
indispensable el diseño e instalación de cadenas de aisladores. Al aumentar los
niveles de tensión en la línea, también es necesario aumentar la altura de las torres,
esto con el fin de mantener el espacio requerido entre los conductores y tierra, de
la misma forma incrementa la cadena de los aisladores.
Como ya se mencionó, los aisladores al igual que las cadenas de aisladores tienen
como objetivo el soporte de los conductores eléctricos en la torre de transmisión, su
función es eludir la fuga de corriente y evitar en la medida de lo posible las pérdidas
de potencia eléctrica.
Estos elementos son sumamente relevantes ya que en ellos recae la
responsabilidad de soporte, y por ende la seguridad de las líneas cuando se
manejan niveles de tensión altos.
1.6 Herrajes
Los herrajes son componentes esenciales en la infraestructura de las lineas de
transmision, se emplean para la sujeción y conexión del cable al aislador y de este
a la estructura metálica.
Se dividen en 2 categorías:
1. Herrajes que conforman arreglos o conjuntos
Son grupos de piezas metálicas diseñadas para cumplir con una función específica,
algunos de los tipos mas comunes incluyen:

• Herrajes de fijación: Permiten la union de cadenas de aisladores a

estructuras. Ejemplos de esto son los grilletes y las articulaciones tipo
bisagra, que ofrecen estabilidad y movilidad controlada en las conexiones.
• Herrajes de unión: Unen los aisladores con los conductores. Entre estos se
encuentran los ojos, anillos y yugos, los cuales facilitan la conexión entre
diversos componentes, garantizando que la alineacion y sujecion sea la
correcta.
• Mordazas:Unen los conductores a las cadenas de aisladores, se distinguen
2 tipos principales, las mordazas de suspensión y las mordazas de amarre,
las cuales aseguran el conductor en su lugar, ya sea mediante corte o sin
corte del material conductor.
• Herrajes de protección: Mitigan los efectos del arco eléctrico y el efecto
corona, contribuyendo a la seguridad y al rendimiento eficiente de la línea de
transmisión.
• Herrajes para cables de guarda: Presentan similitudes con los que se utilizan
para los conductores de fase, la diferencia radica en su tamaño y materiales

18
 empleados parta proteger la línea de transmisión de eventos externos y
minimizar el impacto de descargas atmosféricas.

2. Herrajes independientes
Son elementos individuales que se sitúan a lo largo de las líneas de transmisión,
con el propósito de cumplir funciones específicas, como el soporte, la separación y
el aseguramiento de los conductores.
1.7 Cargabilidad de una línea
Es la cantidad máxima de potencia eléctrica que se puede trasladar en una línea de
transmisión, regularmente manteniendo las condiciones normales de operación o al
menos en condiciones aceptables de la misma línea. La cargabilidad es
dependiente del calibre y longitud de la línea, por lo que para comprender su
comportamiento se deben tener valores de tensión en ambos extremos de la línea,
siendo de 1 p.u [12][14].
La cargabilidad de una línea puede ser valorada de acuerdo a la relación entre la
potencia real y la potencia reactiva (P-Q), mediante la propia línea eléctrica, siendo
variante debido a la propia potencia de la línea, asimismo de la longitud que se
tenga. Propiamente, la cargabilidad tiene una unidad que, involucra la relación de
potencia (P-Q), el valor de la potencia real donde se requiere una potencia reactiva
de valor cero se le conoce como SIL (Potencia Natural de la Línea) [14].
Dentro de los conceptos involucrados en la cargabilidad se encuentran:

• Capacidad térmica del conductor.

• Caída de tensión.
• Margen de estabilidad.
El análisis de la cargabilidad puede ser planteado desde la operación del Sistema
Eléctrico de Potencia (SEP), tomándose en cuenta el comportamiento de las curvas
características dependientes de la cantidad de potencia sobre las líneas. Valiéndose
de dichas características para cuantificar los niveles de cargabilidad.
Por otra parte, un análisis más minucioso sobre la cargabilidad es para Sistemas
Eléctricos Longitudinales (SEL), ya que los niveles de cargabilidad son inferiores de
acuerdo a las curvas estándar de la relación de potencia, siendo exento de la
aplicación de sistemas de análisis robustos [12].
La cargabilidad puede estar determinada por la capacidad térmica de los
conductores o del equipo de terminales, debido a la caída de tensión o propiamente
por la separación angular entre fuentes de tensión, sin embargo, en sistemas
eléctricos de baja capacidad es inusual el encontrar una cargabilidad definida por la

19
capacidad térmica de los conductores, contrariamente a un sistema robusto, el cual
siendo aun de línea corta, se limita por dicho criterio [12].
1.7.1 Líneas de transmisión de alta cargabilidad
El término HSIL proviene del inglés que significa High Surge Impedance Loading y
se traduce como Líneas de Alta Cargabilidad. Estas líneas de transmisión optimizan
el SIL al ajustar la configuración geométrica de los conductores, lo que incrementa
la capacidad de transmisión [1].
El SIL es un factor clave en la evaluación de la carga máxima que puede transmitir
una línea de alta tensión.
El valor SIL tiene dependencia de la inductancia en serie y de la capacitancia en
derivación, así como de la tensión entre las fases, propiamente de la línea de
transmisión. Por lo que, el nivel de SIL es directamente proporcional a los valores
de tensión, por lo que a mayor tensión mayor es el SIL de la línea.
El SIL se calcula con la ecuación 2.
𝑘𝑉 2 𝐿−𝐿 (2)
𝑆𝐼𝐿 = [𝑀𝑊]
𝑍𝑐
(3)
𝐿
𝑍𝑐 = √ [Ω]
𝐶
Donde:
𝑆𝐼𝐿 = Potencia natural de la línea [MW].
kV 2 L−L = tension de operación de linea a linea [kV]
𝑍𝑐 = Impedancia característica de la línea [Ω].
𝐿 = Inductancia de la línea por unidad de longitud [H].
𝐶 = Capacitancia de la línea por unidad de longitud [F].

La ecuación 3 indica que se puede aumentar la potencia natural de la línea si se

disminuye la impedancia característica, esto es posible si se aumenta la
capacitancia (C) o disminuye la inductancia de la línea [1].
En la práctica es posible lograr esto de las siguientes formas:
• Añadiendo subconductores por fase
• Aumentando el tamaño del haz para cada fase, Tecnología de Haz Expandido
en inglés Expanded Bundle Technology (EXB).
• El uso de haces asimétricos
• Disminuyendo la distancia entre fases.

20
Ciertamente se puede establecer los siguientes puntos:
• Cuando los niveles de transferencia de potencia real sean menores al SIL, la
potencia capacitiva producida será menos que la potencia inductiva sumada,
por lo que bajo estos criterios la aportación de la línea es de potencia reactiva
capacitiva a la red eléctrica.
• Cuando los niveles de transferencia de potencia real sean igual al SIL, la
potencia capacitiva producida será igual que la potencia inductiva
consumida, por lo que, bajo estos criterios, la línea no inyecta, sin embargo,
tampoco absorbe potencia reactiva a la propia red eléctrica.
• Cuando los niveles de transferencia de potencia real son mayores al SIL, la
potencia reactiva generada en la línea eléctrica es mayor a la potencia
reactiva consumida, por lo que la demanda sería de potencia reactiva
inductiva, de tal manera que el suministro debe ser ajeno al Sistema Eléctrico
de Potencia [14].
En líneas de alta cargabilidad existen dos tipos de clasificación: [2].
• Líneas HSIL convencionales.
• Líneas HSIL no convencionales.
Líneas HSIL convencionales
Las líneas de transmisión HSIL convencionales se caracterizan por tener
conductores en configuración circular con una radio uniforme para todas las fases.
Aunque mantiene la misma configuración que las líneas de transmisión
tradicionales, su diseño permite una mayor capacidad de transmisión. En la Figura
10 se muestra un ejemplo de una línea HSIL convencional, donde los haces de
conductores son simétricos [2].

Figura 10. Línea de transmisión HSIL convencional [1].

21
Líneas HSIL no convencionales
Las líneas de transmisión HSIL no convencionales se diferencian por su disposición
asimétrica de los conductores, lo cual es evidente en las fases laterales de la línea
presentada en la Figura 11 [2].

Figura 11. Línea de transmisión HSIL no convencional [1].

Esta tecnología es innovadora y en la actualidad, solo hay seis líneas de transmisión
HSIL operativas en ciertos países. En la Tabla 2 se presenta el listado de las líneas
de transmisión HSIL instaladas actualmente, así como el país donde se encuentran
ubicadas [1].
Tabla 2. Líneas de transmisión HSIL existentes actualmente [1].

Ubicación Características de la línea

Línea de transmisión de 744 km de 500 kV con cuatro

Brasil.
subconductores por fase.

Línea de transmisión de 82.5 km de 500 kV con seis

China
subconductores por fase.

China Línea de 30 km de 220 kV con cuatro subconductores por fase.

Rusia Línea de 150 km de 330 kV con cuatro subconductores por fase.

Rusia Línea de 800 km de 500 kV con seis subconductores por fase.

Rusia Línea de 150 km de 220 kV con dos subconductores por fase.

22
1.7.2 Capacidad técnica del conductor
Refiriéndose a la limitación térmica, resultando del calor producido por las pérdidas
óhmicas, siendo dependiente de la máxima temperatura del conductor, de esta
manera, se ve afectado el valor de la flecha de las estructuras (torres) y propiamente
caídas de tensión. Superado dicho límite, puede verse afectada la distancia
permisible con respecto a la tierra, además de exceder el límite de elasticidad del
propio conductor, por lo que al enfriarse no tendería a regresar a su longitud original.
La temperatura del conductor depende de factores propios del conductor como lo
es la magnitud de la corriente, así como de factores como el viento, la temperatura
del ambiente y condiciones físicas sobre el propio conductor. Las duraciones de
estos límites pueden verse en corto y largo plazo. En líneas de alta tensión se tienen
consideraciones a efectos como el efecto corona y los efectos de campo para la
capacidad térmica, siendo de esta manera que el diseño se establezca con
capacidades térmicas elevadas [12][13].
1.7.3 Caída de tensión
La limitante de caída de tensión es muy importante, primordialmente en los sistemas
eléctricos longitudinales, ya que mientras las longitudes de una línea sean mayores,
mayor será la caída de tensión. De acuerdo al modelo de la Figura 12 la caída de
tensión entre el nodo de envió 𝑉𝐸 y el nodo de recepción 𝑉𝑅 se define mediante la
ecuación 4: [12][13][15]
𝑉𝐸 − 𝑉𝑅 (4)
𝐶𝑉 = ∗ 100
𝑉𝑅
Donde:
CV = caída de tensión [%]
VE = tensión en el nodo emisor [kV]
VR = tensión en el nodo receptor [kV]
Por otra parte, la tensión límite de la terminal receptora 𝑉𝑅𝐿 se determina por la
ecuación 5:
𝐶𝑉 (5)
𝑉𝑅𝐿 = 𝑉𝐸 [1 − ]
100
Donde:
VRL = tensión límite en el nodo receptor [kV]
Usualmente se le atribuye un 5% a la caída de tensión permisible, sin embargo,
puede ser catalogado como un valor estricto, sin embargo, al tener un amento por
la cargabilidad se puede considerar una caída de 7.5% [12] [16].

23
En la Figura 12 se puede observar el modelo de análisis de cargabilidad en líneas
de transmisión.

Figura 12. Modeló para análisis de cargabilidad en LT [12].

1.7.4 Margen de estabilidad
La limitante de estabilidad se refiere al margen que hay entre la potencia máxima
(𝑃𝑚á𝑥 ) y la potencia permisible de operación (𝑃𝑜𝑝 ). Por lo que el margen de
estabilidad (ME) se define de acuerdo con la ecuación 6:
𝑃𝑚á𝑥 − 𝑃𝑜𝑝 (6)
𝑀𝐸 = ∗ 100
𝑃𝑜𝑝
Donde:
ME = límite de estabilidad [%]
𝑃𝑚á𝑥 = potencia máxima [MW]
𝑃𝑜𝑝 = potencia permisible de operación [MW]

Dicho margen de estabilidad es seleccionado con el fin de tener una operación

estable, debido a las adversas eventualidades que puedan provocar cambios de
cargas transitorias y ciertamente, de régimen permanente en la línea. Tales cambios
pueden ser producidos por maniobras de apertura y cierre de transformadores o en
las mismas líneas, cambios del despacho de generación eléctrica o por
perturbaciones eléctricas como lo son: fallas o pérdidas en la generación de energía.
El margen es basado en los sistemas existentes, así como en los criterios de una
planificación de la confiabilidad de un sistema. Regularmente se utiliza un margen
de estabilidad del 30%-35%. Para un margen de estabilidad del 35% (0.35 pu) se
puede utilizar la ecuación 7:

24
 𝑃𝑜𝑝 𝑃𝑜𝑝 (7)
0.35 = 1 − ; = 0.65
𝑃𝑚á𝑥 𝑃𝑚á𝑥
Y considerando la simplificación de la línea sin pérdidas, se tiene que:
𝑃𝑜𝑝 = 𝑃𝑚á𝑥 sin(𝛿𝐸 )
(8)
𝛿𝐸 = 𝑠𝑖𝑛−1 (0.65) = 40°

1.8 Parámetros de la línea

Antes de llevar a cabo los cálculos de cualquier sistema eléctrico, es fundamental
contar con una serie de datos relevantes, como los componentes que lo integran,
su comportamiento, las interacciones entre ellos, los métodos de cálculo, las formas
de agruparlos, las técnicas numéricas.
Cualquier circuito eléctrico se compone de algunos de los siguientes elementos:
resistencia, inductancia, capacitancia y conductancia. Estos componentes forman
la base de todos los sistemas eléctricos actuales, desde circuitos simples hasta los
sistemas de potencia más complejos, como lo son las líneas de transmisión [1].
Los parámetros de las líneas de transmisión que se muestran en la Figura 13 son
los siguientes:
• Resistencia (R), representan perdidas a lo largo del conductor, expresada en
Ω/m.
• Inductancia (L), representa el flujo magnético debido a la corriente que fluye
a través de conductor, expresada en H/m.
• Capacitancia (C), representa la corriente del desplazamiento que fluye entre
dos conductores en el plano transversal, expresada en F/m.
• Conductancia (G=1/R), expresada en S/m, está en función de la
conductividad del dieléctrico de la línea (usualmente la conductancia (G)
puede despreciarse en líneas aéreas).

Figura 13. Parámetros de una línea de transmisión [].

25
Los parámetros eléctricos que afectan significativamente a lo largo de una línea
eléctrica son la resistencia y la inductancia. Aunque se combinan para formar la
impedancia (una vez que la inductancia se convierte en reactancia inductiva), su
comportamiento eléctrico es bastante diferente. Por un lado, el efecto inductivo de
las bobinas permite almacenar energía eléctrica en forma de campo magnético. Por
otro lado, la resistencia en un circuito eléctrico genera pérdidas debido al efecto
Joule
1.8.1 Resistencia
La resistencia es la medida de la oposición que presenta un material al flujo de
corriente eléctrica. Aunque su análisis datos de los primeros descubrimientos en el
ámbito eléctrico, no se conectó con otros parámetros eléctricos hasta que George
Simón Ohm desarrolló su ley fundamental. Esta ley, que es fundamental para la
electricidad, relaciona la resistencia con la tensión o diferencia de potencial y la
intensidad [1].
Conceptualmente, la resistencia de un conductor depende de sus dimensiones
físicas y de su resistividad, y se puede expresar como se muestra en la ecuación 9:
𝜌𝑙 (9)
𝑅𝑑𝑐 = [Ω]
𝑆
Donde:
𝑅𝑑𝑐 = 𝑅𝑒𝑠𝑖𝑠𝑡𝑒𝑛𝑐𝑖𝑎 𝑑𝑒𝑙 𝑐𝑜𝑛𝑑𝑢𝑐𝑡𝑜𝑟 𝑒𝑛 𝑐𝑜𝑟𝑟𝑖𝑒𝑛𝑡𝑒 𝑑𝑖𝑟𝑒𝑐𝑡𝑎, 𝑒𝑥𝑝𝑟𝑒𝑠𝑎𝑑𝑎 𝑒𝑛 Ω.
𝑙 = 𝐿𝑜𝑛𝑔𝑖𝑡𝑢𝑑 𝑑𝑒𝑙 𝑐𝑜𝑛𝑑𝑢𝑐𝑡𝑜𝑟, 𝑒𝑥𝑝𝑟𝑒𝑠𝑎𝑠 𝑒𝑛 𝑚.
𝑆 = 𝑆𝑒𝑐𝑐𝑖ó𝑛 𝑡𝑟𝑎𝑠𝑣𝑒𝑟𝑠𝑎𝑙 𝑑𝑒𝑙 𝑐𝑜𝑛𝑑𝑢𝑐𝑡𝑜𝑟, 𝑒𝑥𝑝𝑟𝑒𝑠𝑎𝑠 𝑒𝑛 𝑚𝑚2.
𝜌 = 𝑅𝑒𝑠𝑖𝑠𝑡𝑖𝑣𝑖𝑑𝑎𝑑 𝑑𝑒 𝑙𝑎 𝑚𝑎𝑡𝑒𝑟𝑖𝑎, 𝑒𝑥𝑝𝑟𝑒𝑠𝑎𝑑𝑎 𝑒𝑛 Ω · 𝑚𝑚2/𝑚.

En la Tabla 3 se muestra un listado de la resistividad de los materiales para

conductores.
Tabla 3. Resistividad del material de los conductores [11].

Fuerza a la
Resistividad (𝜌) Densidad (𝛿)
Material tracción (𝑋)
[Ω𝑚𝑚2 /𝑚] [kg/d𝑚3 ]
[kg/dc]
Aluminio 0.0260 2.7 14
Cobre 0.0176 8.9 28
Acero 0.0350 7.8 42

En los sistemas de potencia, que operan con altas tensiones e intensidades, se

requieren conductores de gran tamaño para transportar eficientemente esas

26
energías. Uno de los materiales más utilizados es el cobre, que, como otros metales,
cuenta con características de maleabilidad [1].
1.8.2 Inductancia
La inductancia es la propiedad de un componente del circuito que permite
almacenar energía en una bobina en forma de campo magnético. Una característica
clave de la inductancia es que solo se manifiesta en presencia de corriente alterna.
Aunque un elemento puede poseer inductancia debido a sus propiedades
geométricas y magnéticas, esto no tendrá un efecto en el circuito a menos que haya
un cambio en la corriente a lo largo del tiempo (corriente alterna, AC). El campo
magnético generado por un conductor depende de la corriente que fluye a través de
él en un momento específico como se muestra en la Figura 14[1].

Figura 14. Campo magnético y campo eléctrico asociados a un alinea de dos

conductores [].
Las líneas de flujo magnético se disponen en anillos cerrados alrededor de cada
conductor. En cambio, las líneas del campo eléctrico emergen de las cargas
positivas en un conductor y se dirigen hacia las cargas negativas en el otro.
Cualquier cambio en la corriente que fluye a través de los conductores genera una
variación en el número de líneas de flujo magnético que atraviesan el circuito.
Asimismo, cualquier cambio en el campo magnético induce una fuerza electromotriz
(fem) en el circuito, la cual es proporcional a la velocidad de cambio del flujo. La
inductancia es la propiedad de un circuito que relaciona la fem inducida con la
velocidad de cambio de la corriente [1].
La inductancia industrial de una línea se determina en Henrios (H), y se puede
expresar como lo indica la ecuación 10:
𝜇 𝐷𝑒 (10)
𝐿=[ + 2 ∙ ln ( )] 𝑙 ∙ 10−4 ∙ 𝑙𝑜𝑛𝑔 [𝐻]
2∙𝑛 𝑟𝑒
Donde:
𝑛 = Número de cables por fase
𝐷𝑒 = Distancia media geométrica entre fases (mm)

27
𝑟𝑒 = Radio equivalente (mm)
𝑙 = Longitud de la línea (km)
𝜇 = Permeabilidad

La inductancia en las líneas suele expresarse en (H/Km), por lo que la ecuación 11

se puede expresar de la siguiente manera:
𝜇 𝐷𝑒 (11)
𝐿𝐾 = [ + 2 ∙ ln ( )] ∙ 10−4 [𝐻/𝑘𝑚]
2∙𝑛 𝑟𝑒
Al convertir logaritmos neperianos a logaritmos decimales, se obtiene la ecuación
12:
𝜇 𝐷𝑒 (12)
𝐿𝐾 = [ + 4.6 ∙ 𝑙𝑜𝑔 ( )] ∙ 10−4 [𝐻/𝑘𝑚]
2∙𝑛 𝑟𝑒
La permeabilidad está influenciada por las propiedades del material y por las
condiciones eléctricas a las que se somete el conductor.
• Permeabilidad: 𝜇 = 𝜇0 ∙ 𝜇𝑟
• Permeabilidad absoluta: 𝜇 = 4 ∙ 𝜋 ∙ 10−7
Los conductores de cobre, aluminio y otras aleaciones tienen un valor de
permeabilidad de 𝜇 = 1. Por otro lado, los conductores de acero galvanizado tienen
una permeabilidad de 𝜇 = 200 [1].
Dado que se suelen emplear conductores de cobre, aluminio o sus aleaciones, y la
mayoría tiene un coeficiente de permeabilidad similar, podemos reemplazar este
valor en la ecuación 12.
1.8.3 Capacitancia
La capacidad de una línea de transmisión de energía eléctrica surge de la diferencia
de potencial entre los conductores que la componen. Esta diferencia provoca que
los conductores se carguen de manera similar a las placas de un condensador
cuando hay una diferencia de potencial entre ellas. La capacidad entre conductores
paralelos se define como la carga por unidad de diferencia de potencial, y es una
constante que depende del tamaño de los conductores y de la distancia que los
separan. El efecto de la capacidad suele ser pequeño y puede ser ignorado en
líneas eléctricas de menos de 80 km de longitud; Sin embargo, en líneas más largas,
se convierte en un parámetro importante a considerar [1].
En la Figura 15 las líneas de flujo eléctrico emergen de las cargas positivas y se
dirigen hacia las negativas. La densidad de carga que es perpendicular a la
superficie se representa como (B) y se calcula con la fórmula (ε·E), donde (ε)
corresponde a la permitividad del material que envuelve la superficie y (E) es la
intensidad de campo eléctrico [1].

28
 Figura 15. Superficie gaussiana con cargas eléctricas en su interior [].
Existen otras fórmulas útiles para describir la capacitancia de un circuito, como se
puede apreciar en la ecuación 13. En particular, la capacitancia de una línea que
consta de dos conductores puede determinarse de manera indirecta si se conocen
la carga y la diferencia de potencial que se encuentra entre los conductores [1].
𝑄 (13)
𝐶= [𝐹]
𝑉
Donde:
𝐶 = Capacitancia (F)
𝑞 = carga sobre la línea (c)
𝑉 = es la diferencia de potencial entre los conductores (V)
1.8.4 Conductancia
La conductancia es otro parámetro clave en los circuitos eléctricos convencionales.
Se refiere a la facilidad con la que un material permite el flujo de corriente eléctrica,
lo que significa que es la inversa de la resistencia. Al igual que la capacitancia, es
un parámetro transversal, una diferencia de la resistencia o la inductancia. Su
unidad es el siemens (S) [1].
Una fórmula típica de representar la conductancia se presenta en la ecuación 14:
𝑃 (14)
𝑉 1 1 𝐼 𝑉 𝑃
𝑃 =𝑉∙𝐼; 𝑐𝑜𝑚𝑜 𝐼 = ; 𝑒𝑛𝑡𝑜𝑛𝑐𝑒𝑠 𝐺 = = = = = 2
𝑅 𝑅 𝑉 𝑉 𝑉 𝑉
𝐼
Donde:

𝐺 = Conductancia (S)
𝑃 = Potencia (W)
𝐼 = Corriente eléctrica (A)
𝑉 = Diferencia de potencial entre los conductores (V)

29
Al expresar la ecuación 14 en kilómetros de recorrido y en valores de fase
obtenemos la ecuación 15:
𝑃𝑓𝑎𝑠𝑒−𝑘𝑚 (𝑘𝑊) −3
𝑆 𝑈𝑙𝑖𝑛𝑒𝑎 (15)
𝐺𝑘 = ∙ 10 [ ∙ 𝑓𝑎𝑠𝑒] 𝐶𝑜𝑛 𝑉𝑓𝑎𝑠𝑒 =
𝑉 2𝑓𝑎𝑠𝑒 (𝑘𝑉)2 𝑘𝑚 √3
La conductancia considera las corrientes de fuga que ocurren tanto en los aisladores
que sostienen las líneas aéreas como en las pérdidas que se producen cuando los
electrones saltan a través del aire. Este parámetro está influenciado por diversos
factores, incluidos los climáticos y ambientales, que son difíciles de predecir y
pueden variar a lo largo de toda la línea.
Los valores de conductancia suelen ser relativamente bajos en comparación con los
efectos resistivos, inductivos o capacitivos mencionados anteriormente. Es
afortunado que la conductancia contribuya solo de manera limitada al total de los
efectos eléctricos en un circuito, ya que su cálculo preciso es extremadamente difícil
y, por esta razón, a menudo se desprecia. [1]
1.9 Derecho de vía
Es una franja de terreno que se ubica a lo largo de cada línea aérea, cuyo eje de
longitud coincide con el trazo topográfico de la línea, en caso de las líneas
construidas como tipo lindero, el ancho de derecho de vía se debe considerar a
partir del eje longitudinal del conductor. Su dimensión transversal varía de acuerdo
con el tipo de estructuras, con la magnitud y desplazamiento lateral de la flecha y
con la tensión de operación [19].
En caso de líneas aéreas multicircuito con diferentes niveles de tensión, el derecho
de vía debe ser definido por el circuito con la mayor tensión nominal de operación.
El derecho de Vía de una línea aérea se determina en función del tipo de estructura,
en el momento de que exista un cambio de tipo de estructura, el derecho de Vía se
ajusta a las nuevas condiciones. Para la determinación del derecho de Vía se debe
usar la torre de suspensión predominante sobre el mismo.
• Objetivo:
El objetivo del derecho de Vía es disponer del espacio físico bajo la línea de
transmisión, esto con la intención de que permita su adecuada operación con la
mayor confiabilidad y menores salidas, en beneficio del servicio eléctrico, de tal
manera que se faciliten las actividades de mantenimiento con el mínimo de
interferencias; además de salvaguardar la integridad de los residentes ubicados de
forma aledaña a la propia línea, evitando así la posibilidad de accidentes, debido a
los altos niveles de tensión o fenómenos de inducción.
El ancho del derecho de vía se compone del doble de la suma de las siguientes
distancias: la separación horizontal mínima de seguridad (distancia A), la proyección
horizontal de la flecha del conductor y la longitud de la cadena de aisladores de
suspensión, si corresponde, en función del ángulo de oscilación provocado por la

30
presión del viento (distancia B), y la distancia desde el eje de la estructura hasta el
conductor extremo en reposo (distancia C), tal como se aprecia en la Figura 16 y en
la ecuación 16.

Figura 16. Integración del derecho de Vía [20].

Ecuación para el cálculo del ancho de vía:
𝐴𝑛𝑐ℎ𝑜 𝑑𝑒 𝑉í𝑎 = 2(𝐴 + (𝐿𝑎 + 𝑓16℃ ) sin 𝛼 + 𝐶) (16)

Donde:
𝐴 = 𝑆𝑒𝑝𝑎𝑟𝑎𝑐𝑖ó𝑛 ℎ𝑜𝑟𝑖𝑧𝑜𝑛𝑡𝑎𝑙 𝑚í𝑛𝑖𝑚𝑎 𝑑𝑒 𝑠𝑒𝑔𝑢𝑟𝑖𝑑𝑎𝑑.
𝐵 = 𝑃𝑟𝑜𝑦𝑒𝑐𝑐𝑖ó𝑛 ℎ𝑜𝑟𝑖𝑧𝑜𝑛𝑡𝑎𝑙 𝑑𝑒 𝑙𝑎 𝑓𝑙𝑒𝑐ℎ𝑎 𝑚á𝑠 𝑐𝑎𝑑𝑒𝑛𝑎 𝑑𝑒 𝑎𝑖𝑠𝑙𝑎𝑑𝑜𝑟𝑒𝑠.
𝐶 = 𝐿𝑜𝑛𝑔𝑖𝑡𝑢𝑑 𝑜𝑠𝑐𝑖𝑙𝑎𝑛𝑡𝑒 𝑑𝑒 𝑙𝑎 𝑐𝑎𝑑𝑒𝑛𝑎 𝑑𝑒 𝑎𝑖𝑠𝑙𝑎𝑑𝑜𝑟𝑒𝑠 𝑒𝑛 𝑚𝑒𝑡𝑟𝑜𝑠.
𝐿𝑎 = 𝐹𝑙𝑒𝑐ℎ𝑎 𝑓𝑖𝑛𝑎𝑙 16 º𝐶.

Estos parámetros pueden variar según la tensión eléctrica nominal, el calibre del
conductor, la magnitud de la presión del viento, el tipo de estructura, la ubicación
geográfica y la altitud.

31
 Capítulo II
Marco
Metodológico

32
Métodos de cálculo de interferencia electromagnética en líneas de trasmisión
de alta tensión.
En las líneas de transmisión suele producirse interferencia electromagnética, debido
al fenómeno conocido como efecto corona.
El sistema de transmisión es diseñado con el objetivo de que el efecto corona sea
disminuido, debido a que éste, merma la capacidad de transmisión de energía.
Existen diversas condiciones que influyen en el efecto corona que se presenta en
las líneas de transmisión, como lo son; la temperatura, humedad, tamaño de los
conductores y geometría de su haz, tensión de red, entre otras.
En este capítulo se emplean los métodos de cálculo de interferencia
electromagnética, en su mayoría semiempiricos, en dichos métodos se consideran
diversos parámetros de la línea en cuestión, enfocados principalmente en las
dimensiones y cantidad de conductores, distancia entre ellos y gradiente máximo
de la superficie del conductor.
En México no existe una norma que regule los niveles de interferencia
electromagnética en las líneas de transmisión, por lo que en el presente proyecto
se tomó como referencia una norma canadiense EPRI (Electric Power Research
Institute, en español: Instituto de Investigación de Energía Eléctrica), que define los
límites máximos de RI en líneas de transmisión tomando como base la tensión
nominal de fase a fase.

2.1 Efecto corona.

Uno de los efectos presentes en las líneas de transmisión de alta tensión es el efecto
corona, debido a la ionización que se produce con el aire que rodea a los
conductores de estas, debido a que las secciones de los conductores son circulares
los cuales generan un haz luminoso similar a una corona, puede tomar un color
rojizo cuando los efectos sean tenues o un color azulado cuando los efectos sean
más drásticos sobre el conductor, tales colores dependen del nivel de tensión en la
línea de transmisión, mientras más altos sean dicho niveles más incrementa la
luminosidad de la corona y su calor. [1]

Por otra parte, sobre la superficie del conductor puede producir ozono, en caso de
presentarse niveles de humedad en el contexto, se genera ácido nitroso; por tales
motivos el efecto corona puede ser percibido debido a los chisporroteos y sonidos
silbantes que se producen. Este efecto puede ser presenciado tanto en corriente
alterna (C.A) como en corriente continua (C.D), así como en polaridad positiva y
polaridad negativa. En la Figura 17 se muestra el fenómeno [2].

33
 Figura 17. Conductor con efecto corona [3].

Durante el efecto corona se tienden a presentar pérdidas de potencia, consecuencia

de la adquisición de una diferencia de potencial elevado suficiente para dar hincapié
a lo que se le conoce como gradiente del campo eléctrico radial. [3]

Concretamente se pueden presentar los siguientes problemas en las líneas de

transmisión por efecto corona: [4].

• Pérdidas de energía.
• Deterioro del material.
• Producción de contaminación.
• Ruido audible.
• Interferencia electromagnética.
• Descargas luminosas.
• Generación de ozono.

2.1.1 Naturaleza y manifestación del efecto corona.

La intensidad del efecto corona es dependiente del campo eléctrico que se
encuentre en el área de los conductores, por ende, dicha intensidad es afectada de
igual manera por la tensión de la línea como se menciona anteriormente, así como
de del diámetro, disposición, estado físico de los conductores y las condiciones
climáticas del entorno. Los conductores de gran diámetro presentan un gradiente
eléctrico más bajo en la propia superficie de los conductores, por lo que el efecto
corona es mayor, contrariamente a los conductores de menor diámetro. En otro tipo
de afecciones se encuentran la presencia de insectos o incluso polvo (suciedad) ya
que provocan irregularidades en las condiciones de las líneas, asimismo con la
presencia de gotas de lluvia, nieve. humedad e incluso niebla [5].

2.1.2 Descripción física del efecto corona

Cuando una corriente alterna circula a través de dos conductores de una línea de
transmisión, y el espacio entre ellos es mayor que sus diámetros, el aire circundante
experimenta un estrés dieléctrico. A niveles bajos de tensión, no ocurre nada, ya

34
que el estrés es insuficiente para ionizar el aire. Sin embargo, al elevarse la tensión
por encima de 30 kV, conocido como voltaje crítico de ruptura, la intensidad del
campo aumenta, y el aire se somete a un estrés suficiente para ionizar el aire. Esto
provoca una descarga eléctrica alrededor de los conductores, debido a la ionización,
lo que también produce un brillo luminoso, ruido audible y en algunos casos emisión
de ozono, en la tabla 4 se muestran los factores que propician ele efecto corona [4].

Tabla 4. Factores que afectan el efecto corona [].

Factores Descripción

Durante una condición de tormentas o lluvias el

efecto corona ocurre a un voltaje mucho menor
en comparación con un buen tiempo. Debido a
Atmósfera
que durante una tormenta o lluvia la densidad
de iones alrededor de los conductores es
mucho mayor que en un buen tiempo.
La irregularidad de la superficie de los
conductores disminuye el voltaje de ruptura.
Por esta razón, con el bajo voltaje del
Tamaño del conductor aislamiento disruptivo se creará chispas y
corona. Es por eso que los conductores sólidos
se utilizan sobre todo en vez de conductores
trenzados para reducir la corona.
Mientras más grande sea la brecha entre los
Separación entre conductores
conductores menor será el efecto corona.
Cada voltaje de línea tiene un límite. Después
de que se alcance el límite disruptivo se
producirá y creará chispas y corona. Así, con la
Tensión de red
aplicación de una menor tensión de línea la
probabilidad de que el efecto corona ocurra
disminuye

2.2 Ecuaciones empíricas y semiempiricas

Es inviable determinar las características de rendimiento del RI únicamente a través
de métodos analíticos. A la par del avance de los métodos teóricos, se han realizado
numerosos estudios experimentales en los últimos cincuenta años, generando una
amplia base de datos sobre los niveles de RI de las líneas de transmisión para
clases distintas de voltaje y bajo diversas condiciones climáticas. Los datos
experimentales recopilados en estas investigaciones han servido para crear
métodos empíricos y semiempíricos [6].

35
Es fundamental diferenciar entre los métodos de cálculo empíricos y semiempíricos.
Los métodos empíricos, o comparativos, se basan en datos experimentales
recolectados de líneas de transmisión en funcionamiento y en líneas de prueba a
gran escala. Por otro lado, los métodos semiempíricos, a veces denominados
semianalíticos, integran la función de excitación de RI obtenida experimentalmente
con enfoques analíticos para el análisis de propagación de RI, permitiendo así
predecir el rendimiento de RI en nuevas configuraciones [6].

2.2.1 Métodos empíricos.

Las fórmulas empíricas o comparativas se derivan del análisis de regresión
estadística aplicado a datos experimentales de RI, considerando diversas variables
como el diámetro del conductor, el gradiente de su superficie y la distancia lateral
en la línea, entre otras. Contar con datos experimentales de buena calidad, libres
de interferencias ambientales, es un requisito crucial para que estas fórmulas
puedan prever con exactitud el nivel de RI en las nuevas configuraciones
propuestas.

Las fórmulas empíricas suelen elaborarse para líneas ubicadas en una determinada
región climática y bajo condiciones específicas, como clima favorable, adverso o
lluvioso, que son típicas de esa área. Así, las fórmulas basadas en datos de una
región particular pueden no proporcionar predicciones precisas de RI en otra zona,
donde la variabilidad anual de las condiciones climáticas es diferente [6].

2.3 Interferencia Electromagnética.

La interferencia electromagnética, conocida también como radio interferencia (RI),
se define como una reducción en la eficacia de un equipo, canal de transmisión o
sistema debido a una perturbación electromagnética.
La naturaleza pulsante de las corrientes generadas por el efecto corona crea
campos electromagnéticos que provocan alteraciones en un amplio espectro de
frecuencias por lo que interfiere con los sistemas de comunicación [4].

2.4 Gradiente Crítico

Cuando el campo eléctrico en la superficie de los conductores de una línea de
transmisión supera un umbral específico, se genera una descarga eléctrica parcial
en su entorno, conocida como descarga corona. Este umbral se conoce como
gradiente crítico y depende de factores como el diámetro del conductor, la rugosidad
de su superficie y la densidad del aire circundante. Además, el tipo de tensión
aplicada y la polaridad del conductor en corriente continua también influyen en este
valor. Gracias a estudios experimentales en laboratorio, se ha determinado el
gradiente crítico para conductores cilíndricos [7].

36
El investigador Whitehead realizó diversos estudios de laboratorio en este campo,
analizando el comportamiento del gradiente crítico en corriente continua para polos
positivos y negativos. Basándose en sus resultados, Whitehead derivó una fórmula
empírica (ecuación 17) para calcular el gradiente crítico en conductores cilíndricos
lisos de pequeño diámetro [7].

0.301
𝐸𝑐 = 30𝑚𝛿 [1 + ] [𝑘𝑉/𝑐𝑚] (17)
√𝛿𝑟

Donde:
𝐸𝑐 = Gradiente crítico [𝑘𝑉/𝑐𝑚].
𝑚 = Factor del estado superficial del conductor [kV/cm].
𝑟 = Radio del conductor [cm].
𝛿 = Densidad relativa del aire.

2.5 Cálculo del campo eléctrico en la superficie de los conductores.

Debido a la importancia que tiene el gradiente critico en el cálculo de la interferencia

electromagnética, resulta necesario con tal con un método confiable para
determinar este parámetro.

Existen numerosos métodos para evaluar el campo eléctrico en superficies de

conductores individuales y en configuraciones de haz. La exactitud de los gradientes
calculados no solo depende de la calidad del método utilizado, sino también de la
precisión de los datos de entrada. Entre los métodos disponibles, los que ofrecen
una mejor relación entre simplicidad y precisión son el método de Maxwell y el de
Mark y Mengele [8].

Para aplicar estos métodos, se consideran las siguientes suposiciones [8].

• Se asume que la superficie del suelo es plana, conductora, horizontal e

infinita.
• Los conductores se consideran cilíndricos, lisos, infinitamente largos y
paralelos entre sí y respecto a la tierra.
• Se establece que los conductores son superficies equipotenciales con un
potencial conocido, mientras que el plano de tierra tiene un potencial de cero.
• Se ignora la influencia de las estructuras que sostienen los conductores y de
cualquier otro objeto cercano.

37
2.5.1 Método de los coeficientes de potenciales de Maxwell.

Este método se puede utilizar para líneas que cuentan con solo un conductor por
fase en corriente alterna (CA) o por polo en corriente continua (CC). La técnica
implica sustituir la carga en la superficie de los conductores por una carga de
longitud infinita situada en el centro del conductor [9].

Para una línea con varios conductores por haz, el primer paso es calcular la carga
𝑄 de cada uno de ellos [10]:

[𝑄] = [𝑃]−1 ∙ [𝑉] [𝑘𝑉/𝑐𝑚] (18)

Donde:
[𝑉] = Vector de voltajes en los conductores reales con dimensión nx1 [V/m].
[𝑃] = Matriz de coeficientes de potencial o coeficientes de Maxwell con dimensión
nxn [m/F].
[𝑄] = Vector de las densidades de cargas lineales de los conductores reales con
dimensiones nx1 [C/m].

Los elementos de la matriz [𝑃]:

1 2𝐻𝑖 𝑚 2𝐻𝑖 𝑚
𝑃𝑖𝑖 = ln ( ) [ ] = 18 × 109 ln ( )[ ] (19)
2𝜋𝜀0 𝑅𝑖 𝐹 𝑅𝑖 𝐹

1 𝐷𝑖𝑗 𝑚 𝐷𝑖𝑗 𝑚
𝑃𝑖𝑗 = ln ( ) [ ] = 18 × 109 ln ( ) [ ] (20)
2𝜋𝜀0 𝑑𝑖𝑗 𝐹 𝑑𝑖𝑗 𝐹

Para determinar la distancia entre el conductor y la imagen del conductor, así como
las distancia entre subconductores se emplean las siguientes ecuaciones:

2 2
𝐷𝑖𝑗 = √(𝑥𝑗 − 𝑥𝑖 ) + (𝑦𝑗 − 𝑦𝑖 ) [𝑚] (21)

2 2
𝑑𝑖𝑗 = √(𝑥𝑗 − 𝑥𝑖 ) + (𝑦𝑗 − 𝑦𝑖 ) [𝑚] (22)

Donde:
𝑃𝑖𝑖 = Diagonal de la matriz de potenciales propios de cada haz de conductores.
𝑃𝑖𝑗 = Potenciales mutuos que rellenan la matriz.
1
𝜀0 = ∙ 10−9 = Permitividad del aire [F/m].
2𝜋𝜀0

38
𝑅𝑖 = Radio del conductor 𝑖 [m].
𝐻𝑖 = Altura del conductor 𝑖 al plano tierra [m].
𝐷𝑖𝑗 = Distancia entre el conductor 𝑖 y la imagen del conductor 𝑗 al plano tierra [m].
𝑑𝑖𝑗 = Distancia entre los conductores 𝑖 y 𝑗 [m].

En la Figura 18 se muestra un esquema del diagrama utilizado para el cálculo

de los coeficientes de Maxwell.

Figura 18. Diagrama para el cálculo de los coeficientes de Maxwell [10].

[Link] Radio Medio Geométrico.

El radio medio geométrico es un concepto matemático que resulta valioso para

calcular la inductancia. Se define como el radio de un conductor tubular con un
grosor de pared infinitesimal, considerando tanto el flujo interno como el flujo externo
a una distancia unitaria desde el centro del conductor [11].

El radio medio geométrico para una configuración específica se calcula como la

media geométrica del radio medio geométrico de cada subconductor de una fase,
teniendo en cuenta las distancias entre ese subconductor y los demás de la misma
fase [11].

Radio Medio Geométrico RMG [11].

• Para un circuito trifásico sencillo el RMG es el RMG del conductor.

• Para una línea trifásica de doble circuito, donde existen entre los conductores
de una misma fase distancias Daa’, Dbb’ y Dcc’ y es un circuito transpuesto.

39
El RMG de un conductor con n número de hilos puede calcularse por medio de la
Ecuación 23:

𝑛2
𝑅𝑀𝐺 = √𝑟′(𝐷𝑎𝑏 𝐷𝑎𝑐 … 𝐷𝑎𝑚 )𝑛 [m] (23)

Donde:
𝑅𝑀𝐺 = Radio medio geométrico del haz [m]
𝑟′ = Número de subcontuctores del haz
𝑟′ = Radio medio geométrico del conductor del conductor [m].
𝐷𝑎𝑏 = Distancia entre el subconductor de referencia y el resto de los subconductores
[m].

En la Figura 19 se puede apreciar un haz de 4 subconductores y el modelo del radio

medio geométrico para esa configuración:

Figura 19. Haz de 4 subconductores [11].

2.5.2 Método de Mark y Mengele.

Este método para calcular el campo eléctrico en la superficie de los conductores

considera por separado el efecto de la superficie del suelo y de los haces de
conductores. En el desarrollo del método, cada haz de conductores se reemplaza
por un conductor equivalente. Una vez calculada la carga de estos conductores
equivalentes, la carga de cada conductor del haz se determina dividiendo la carga
total del haz entre el número de conductores presentes, colocándose esta carga en
el centro de cada conductor. Esto implica que se supone que la carga en la
superficie de cada conductor del haz es la misma. Para que los valores del campo
eléctrico calculado sean válidos, es necesario que la distancia entre las fases y la
distancia entre las fases y el suelo sean considerablemente mayores que el radio
del haz [9].

40
 • Cálculo del radio equivalente:

El haz de conductores es reemplazado por un conductor equivalente de radio 𝑅𝑒 ,

que está dado por:

𝑁
𝑅𝑒 = √𝑁𝑟𝑅ℎ 𝑁−1 [𝑐𝑚] (24)
Donde:
𝑅𝑒 = Radio equivalente [cm].
𝑁 = Numero de conductores por haz.
𝑟 = Radio de cada conductor del haz o subconductor [cm].
𝑅ℎ = Radio del haz [cm].

𝑠
𝑅ℎ = [𝑐𝑚] (25)
2 𝑠𝑒𝑛(𝜋/𝑁)

Donde:
𝑠 = Distancia entre los conductores del haz [cm].

Después de substituir cada haz de conductores por su respectivo conductor

equivalente se procede a calcular los potenciales de Maxwell.

El campo eléctrico en medio de la superficie de cada uno de los conductores del

haz se concentra y se conoce como campo eléctrico medio.
El campo eléctrico medio se calcula con la siguiente ecuación:

𝑄
𝐸𝑚 = [𝑘𝑉/𝑐𝑚] (26)
2𝜋𝜀0 𝑁𝑟

Donde:
𝐸𝑚 = Campo eléctrico medio [kV/cm].
𝑁 = Numero de conductores por haz.
𝑟 = Radio de cada conductor del haz o subconductor [cm].
𝑄 = Potenciales de Maxwell [C/m].
1
𝜀0 = 2𝜋𝜀 ∙ 10−9 = Permitividad del aire [F/m].
0

El valor del campo eléctrico varía en torno del valor promedio, siendo su valor
Máximo y se calcula con la siguiente ecuación:

𝑟
𝐸𝑚𝑎𝑥 = 𝐸𝑚 [1 + (𝑁 − 1)] [𝑘𝑉/𝑐𝑚] (27)
𝑅ℎ

41
En la Figura 20 se puede apreciar el esquema del campo eléctrico en la superficie
de los conductores:

Figura 20. Esquema del campo eléctrico en la superficie de los conductores [9].

2.6 Ecuaciones para el Cálculo de la Radio Interferencia.

En este apartado se muestran algunas de las ecuaciones que han sido

desarrolladas por diversas organizaciones con el fin de determinar los niveles de RI
de los conductores individuales de las líneas de transmisión de corriente alterna.

2.6.1 CIGRE

Ecuación determinada por el CIGRE (Conseil International des Grands Réseaux

Électriques, en español: Consejo Internacional de Grandes Redes Eléctricas), sigue
las especificaciones: CP-9kHz-0.5Mhz) [12].

• Lluvia densa:

𝑫 𝒒
𝑹𝑰𝑳𝑳𝑫𝒇 = −𝟏𝟎 + 𝟑. 𝟓 𝒈 + 𝟔 𝒅 − 𝟑𝟑 𝐥𝐨𝐠 𝟏𝟎 ( ) + [𝒅𝑩(𝝁𝑽/𝒎)] (28)
𝟐𝟎 𝟑𝟎𝟎

Donde:
𝑅𝐼𝐿𝐿𝐷𝑓 = Radio interferencia de cada fase en condición de lluvia densa [𝑑𝐵(𝜇𝑉/𝑚)].
𝑔 = Gradiente de voltaje máximo promedio [kV/cm].
𝑑 = Diámetro de conductor [cm].
𝐷 = Distancia radial entre el conductor y el punto de medición [m].
𝑞 = Altitud [msnm].

42
Cálculo de distancias radiales:

Para el cálculo de las distancias de radiales de cada fase se emplea la ecuación 29.

𝐷𝑓 = √(𝑦𝑚 − 𝑦𝑖 )2 + (𝑥𝑚 − 𝑥𝑖 )2 [𝑚] (29)

Donde:
𝐷𝑓 = Distancia radial por fase [𝑚].
𝑦𝑚 = Altura de las fases [𝑚].
𝑥𝑖 = Distancias laterales de las fases fases [𝑚].

Con los resultados obtenidos de radio interferencia de cada fase, se realiza el

análisis de la magnitud de cada nivel de 𝑅𝐼𝐿𝐿𝐷 , ya que se debe de ordenar de manera
descendiente para poder obtener la RI resultante de acuerdo a las fórmulas 𝐶𝐼𝑆𝑃𝑅
30 y 31.

𝑆𝑖 𝑅𝐼1 − 𝑅𝐼2 ≥ 3 𝑑𝐵 𝐸𝑛𝑡𝑜𝑛𝑐𝑒𝑠 𝑅𝐼𝑇𝑂𝑇𝐴𝐿 = 𝑅𝐼1 (30)

O bien:

𝑆𝑖 𝑅𝐼1 − 𝑅𝐼2 ≤ 3 𝑑𝐵 𝐸𝑛𝑡𝑜𝑛𝑐𝑒𝑠 𝑅𝐼𝑇𝑂𝑇𝐴𝐿

𝑅𝐼1 + 𝑅𝐼2 (31)
= + 1.5
2

Donde:
𝑅𝐼1 = La magnitud mayor de 𝑅𝐼𝐿𝐿𝐷𝑓 de fase [𝑑𝐵(𝜇𝑉/𝑚)].
𝑅𝐼2 = La magnitud media de 𝑅𝐼𝐿𝐿𝐷𝑓 de fase [𝑑𝐵(𝜇𝑉/𝑚)].
𝑅𝐼3 = La magnitud menor de 𝑅𝐼𝐿𝐿𝐷𝑓 de fase [𝑑𝐵(𝜇𝑉/𝑚)].
𝑅𝐼𝑇𝑂𝑇𝐴𝐿 = 𝑅𝐼𝐿𝐿𝐷 = La magnitud total de radio interferencia del sistema en condición
de lluvia densa [𝑑𝐵(𝜇𝑉/𝑚)].

Siendo la ecuación no. 28 aplicable solo para condición de lluvia densa, para el
cálculo de RI en lluvia moderada y buen ambiente se deben aplicar las ecuaciones
32 y 33 respectivamente.

• Luvia moderada:

𝑅𝐼𝐿𝐿𝑀 = 𝑅𝐼𝐿𝐿𝐷 − 7 [𝑑𝐵(𝜇𝑉/𝑚)] (32)

43
• Buen ambiente:

𝑅𝐼𝐵𝐴 = 𝑅𝐼𝐿𝐿𝐷 − 24 [𝑑𝐵(𝜇𝑉/𝑚)] (33)

Donde:
𝑅𝐼𝐵𝐴 = Radio interferencia total en condición de buen ambiente [𝑑𝐵(𝜇𝑉/𝑚)].
𝑅𝐼𝐿𝐿𝑀 = Radio interferencia total en condición de lluvia moderada [𝑑𝐵(𝜇𝑉/𝑚)].
𝑅𝐼𝐿𝐿𝐷 = Radio interferencia total en condición de lluvia densa [𝑑𝐵(𝜇𝑉/𝑚)].

2.6.2 BPA

Ecuación determinada por la BPA (Bonneville Power Administration, en español:

Administración de Energía de Bonneville), sigue las especificaciones del CISPR:
(CP – 9 kHz – 0.5 MHz) [12].

• Buen ambiente:

𝒈 𝒅
𝑹𝑰𝑩𝑨𝒇 = 𝟒𝟔 + 𝟏𝟐𝟎 𝐥𝐨𝐠 ( ) + 𝟒𝟎 𝐥𝐨𝐠 ( )+
𝟏𝟕. 𝟓𝟔 𝟑. 𝟓𝟏
(34)
𝒒
𝟏𝟎 [𝟏 − (𝐥𝐨𝐠(𝟏𝟎 ∙ 𝒇))𝟐 ] + − 𝑪𝟏 + 𝑪𝟐
𝟑𝟎𝟎

Donde:
𝑅𝐼𝐵𝐴𝑓 = Radio interferencia de cada fase en condición de buen ambiente
[𝑑𝐵(𝜇𝑉/𝑚)].
𝑔 = Gradiente de voltaje máximo promedio en la superficie del conductor [kV/cm].
𝑑 = Diámetro de conductor [cm].
𝑓 = Frecuencia [MHz].
𝑞 = Altitud [msnm].
𝐶1 = Valor de la función de atenuación lateral Ca para un conductor de fase de
referencia [𝑚].
𝐶2 = Valor de la función de atenuación lateral Ca para el tipo de línea al que se
aplica la formula.

- Cálculo de la variable 𝐶𝑎 :

𝐶𝑎 = 10 ∙ log(𝐷𝑊 2 + 𝐸𝑆𝑈 2 + 𝐸𝐼𝑁𝐷2 )

(35)
Donde:
𝐶𝑎 = Función de atenuación lateral.
44
𝐷𝑊 = Componente de onda directa.
𝐸𝑆𝑈 = Componente de onda.
𝐸𝐼𝑁𝐷 = Componente decampo de inducción.

- Cálculo de la variable 𝐶1 :

Para el cálculo de la variable 𝐶1 se emplea la ecuación 35 referida con subíndice 1.

𝐶1 = 10 ∙ log(𝐷𝑊1 2 + 𝐸𝑆𝑈1 2 + 𝐸𝐼𝑁𝐷1 2 ) (36)

Se toman como referencia las siguientes ecuaciones para el cálculo de las variables
de la ecuación 36.

276.16
𝐷1 = 21.04 (2.7) 𝜌1 = (37)
𝜆2
2 + 0.3𝜌1
𝐷𝑊1 = 21.04 (2.7) 𝑓(𝜌1 ) = (38)
2 + 𝜌1 + 0.6𝜌1 2

70.55 31.11 ∙ 𝑓(𝜌1 )

𝐸𝐼𝑁𝐷1 = (2.7) 𝐸𝑆𝑈1 = (39)
𝑓2 𝑓
𝑣
𝜆= (2.7)
𝑓

Donde:
𝜆 = Longitud de onda [m].
𝑣 = Velocidad de la luz [m/s].
𝑓 = Frecuencia de operación [Hz].
𝜌1 = Resistividad del suelo con respecto a parámetros de referencia [Ωm].
𝑓(𝜌1 ) = Función evaluada en P1.

- Cálculo de la variable 𝐶2 :

Para el cálculo de la variable 𝐶2 se emplea la ecuación 35 referida con subíndice 2.

𝐶2 = 10 ∙ log(𝐷𝑊2 2 + 𝐸𝑆𝑈2 2 + 𝐸𝐼𝑁𝐷2 2 ) (40)

Nota: Se deben calcular las distancias radiales por fase “𝐷𝑓 ” con la ecuación 29 para
determinar D2 de cada fase.

De acuerdo con los valores de las distancias radiales calculadas, se calcula 𝐷𝑊2
conforme a las siguientes condiciones:

45
 12ℎ𝑎 ℎ𝑐 ℎ𝑐
Condición 1: Si 𝐷2 ≤ Entonces 𝐷𝑊2 = (41)
𝜆 𝐾𝐷2

12ℎ𝑎 ℎ𝑐 ℎ𝑐 12ℎ𝑎 ℎ𝑐
Condición 2: Si 𝐷2 > Entonces 𝐷𝑊2 = ∙ (42)
𝜆 𝐾 ∙ 𝐷2 𝜆𝐷2
Se aplican las siguientes ecuaciones para determinar el resto de las variables de la
ecuación 40 para cada fase.

2𝜋 2 + 0.3𝜌2
𝐾 = (2.7) 𝑓(𝜌2 ) = (43)
𝜆 2 + 𝜌2 + 0.6𝜌2 2

ℎ𝑐 ℎ𝑐 𝑓(𝜌2 )
𝐸𝐼𝑁𝐷2 = (2.7) 𝐸𝑆𝑈2 = (44)
𝐾 ∙ 𝐷22 𝐾 ∙ 𝐷2

52.5𝐷2
𝜌2 = (2.7)
𝛿𝜆2

Donde:
ℎ𝑐 = Altura del conductor [m].
ℎ𝑎 = Altura de la antena [m].
𝐷2 = Distancia radial entre los conductores y la antena [m].
𝑓 = frecuencia de operación [MHz].
𝛿 = Conductividad del suelo [mS/m].

Se sustituyen los valores de las variables correspondientes en la ecuación 40, para

determinar la RI por cada fase.

Con los resultados obtenidos de RI de cada fase, se realiza el análisis de la

magnitud de cada nivel de 𝑅𝐼𝐵𝐴 , ya que se debe de ordenar de manera descendiente
para poder obtener la RI resultante de acuerdo con las fórmulas 𝐶𝐼𝑆𝑃𝑅 45 y 46:

𝑆𝑖 𝑅𝐼1 − 𝑅𝐼2 ≥ 3 𝑑𝐵 𝐸𝑛𝑡𝑜𝑛𝑐𝑒𝑠 𝑅𝐼𝑇𝑂𝑇𝐴𝐿 = 𝑅𝐼1 (45)

O bien:
𝑅𝐼1 + 𝑅𝐼2
𝑆𝑖 𝑅𝐼1 − 𝑅𝐼2 ≤ 3 𝑑𝐵 𝐸𝑛𝑡𝑜𝑛𝑐𝑒𝑠 𝑅𝐼𝑇𝑂𝑇𝐴𝐿 = + 1.5 (46)
2

Donde:
𝑅𝐼1 = La magnitud mayor de 𝑅𝐼𝐵𝐴𝑓 de fase [𝑑𝐵(𝜇𝑉/𝑚)].
𝑅𝐼2 = La magnitud media de 𝑅𝐼𝐵𝐴𝑓 de fase [𝑑𝐵(𝜇𝑉/𝑚)].
𝑅𝐼3 = La magnitud menor de 𝑅𝐼𝐵𝐴𝑓 de fase [𝑑𝐵(𝜇𝑉/𝑚)].
𝑅𝐼𝑇𝑂𝑇𝐴𝐿 = 𝑅𝐼𝐵𝐴 La magnitud total de radio interferencia del sistema en condición de
buen ambiente [𝑑𝐵(𝜇𝑉/𝑚)].

46
Siendo la ecuación 34 aplicable solo para condición de buen ambiente, para el
cálculo de RI en lluvia moderada y lluvia densa se deben aplicar las ecuaciones 47
y 48 respectivamente:

• Luvia moderada:

𝑅𝐼𝐿𝐿𝑀 = 𝑅𝐼𝑇𝑂𝑇𝐴𝐿 + 17 [𝑑𝐵(𝜇𝑉/𝑚)] (47)

• Lluvia densa:

𝑅𝐼𝐿𝐿𝐷 = 𝑅𝐼𝑇𝑂𝑇𝐴𝐿 + 24 [𝑑𝐵(𝜇𝑉/𝑚)] (48)

Donde:
𝑅𝐼𝐵𝐴 = Radio interferencia total en condición de buen ambiente [𝑑𝐵(𝜇𝑉/𝑚)].
𝑅𝐼𝐿𝐿𝑀 = Radio interferencia total en condición de lluvia moderada [𝑑𝐵(𝜇𝑉/𝑚)].
𝑅𝐼𝐿𝐿𝐷 = Radio interferencia total en condición de lluvia densa [𝑑𝐵(𝜇𝑉/𝑚)].

2.6.3 CREIPI
Ecuación determinada por la CRIEPI (Central Research Institute of Electric Power
Industry, en español: Instituto Central de Investigaciones de la Industria Eléctrica),
sigue las especificaciones del CISPR: (CP – 9 kHz – 0.5 MHz) [12].

• Buen ambiente:

𝒉
𝑹𝑰𝑩𝑨𝒇 = 𝟒 + (𝟑. 𝟕 𝒈 − 𝟏𝟐. 𝟐) + 𝟒𝟎 𝒍𝒐𝒈𝟏𝟎 (𝒅) + 𝟐𝟎 𝐥𝐨𝐠𝟏𝟎 ( 𝟐 )
𝑫
(49)
𝟐 𝒒
−𝟏𝟐 (𝐥𝐨𝐠𝟏𝟎 (𝒇)) − 𝟏𝟕 𝐥𝐨𝐠𝟏𝟎 (𝒇) + [𝒅𝑩(𝝁𝑽/𝒎)]
𝟑𝟎𝟎

Donde:
𝑅𝐼𝐵𝐴𝑓 = Radio interferencia de cada fase en condición de buen ambiente
[𝑑𝐵(𝜇𝑉/𝑚)].
𝑔 = Gradiente de voltaje máximo promedio en la superficie del conductor [kV/cm].
𝑑 = Diámetro de conductor [cm].
ℎ = Altura de la fase [m].
𝐷 = Distancia radial por fase [𝑚].
𝑓 = Frecuencia [MHz].
𝑞 = Altitud [msnm].

Nota: Se deben calcular las distancias radiales por fase “𝐷𝑓 ” con la ecuación no. 29
para determinar D de cada fase.

47
Se sustituyen los valores de las variables correspondientes en la ecuación 49 para
calcular la radio interferencia en buen ambiente de cada fase 𝑅𝐼𝐵𝐴𝑓 .

Con los resultados obtenidos de RI de cada fase, se debe realizar el análisis de la

magnitud de cada nivel de 𝑅𝐼𝐵𝐴 , ya que se debe de ordenar de manera descendiente
para poder obtener la RI resultante de acuerdo con las fórmulas 𝐶𝐼𝑆𝑃𝑅 50 y 51:

𝑆𝑖 𝑅𝐼1 − 𝑅𝐼2 ≥ 3 𝑑𝐵 𝐸𝑛𝑡𝑜𝑛𝑐𝑒𝑠 𝑅𝐼𝑇𝑂𝑇𝐴𝐿 = 𝑅𝐼1 (50)

O bien:
𝑅𝐼1 + 𝑅𝐼2
𝑆𝑖 𝑅𝐼1 − 𝑅𝐼2 ≤ 3 𝑑𝐵 𝐸𝑛𝑡𝑜𝑛𝑐𝑒𝑠 𝑅𝐼𝑇𝑂𝑇𝐴𝐿 = + 1.5 (51)
2

Donde:
𝑅𝐼1 = La magnitud mayor de 𝑅𝐼𝐵𝐴𝑓 de fase [𝑑𝐵(𝜇𝑉/𝑚)].
𝑅𝐼2 = La magnitud media de 𝑅𝐼𝐵𝐴𝑓 de fase [𝑑𝐵(𝜇𝑉/𝑚)].
𝑅𝐼3 = La magnitud menor de 𝑅𝐼𝐵𝐴𝑓 de fase [𝑑𝐵(𝜇𝑉/𝑚)].
𝑅𝐼𝑇𝑂𝑇𝐴𝐿 = 𝑅𝐼𝐵𝐴 = La magnitud total de radio interferencia del sistema en condición de
buen ambiente [𝑑𝐵(𝜇𝑉/𝑚)].

Siendo la ecuación 49 aplicable solo para condición de buen ambiente, para el

cálculo de RI en lluvia moderada y lluvia densa se deben aplicar las ecuaciones 52
y 53 respectivamente:

• Luvia moderada:

𝑅𝐼𝐿𝐿𝑀 = 𝑅𝐼𝐵𝐴 + 17 [𝑑𝐵(𝜇𝑉/𝑚)] (52)

• Lluvia densa:

𝑅𝐼𝐿𝐿𝐷 = 𝑅𝐼𝐵𝐴 + 24 [𝑑𝐵(𝜇𝑉/𝑚)] (53)

Donde:
𝑅𝐼𝐵𝐴 = Radio interferencia total en condición de buen ambiente [𝑑𝐵(𝜇𝑉/𝑚)].
𝑅𝐼𝐿𝐿𝑀 = Radio interferencia total en condición de lluvia moderada [𝑑𝐵(𝜇𝑉/𝑚)].
𝑅𝐼𝐿𝐿𝐷 = Radio interferencia total en condición de lluvia densa [𝑑𝐵(𝜇𝑉/𝑚)].

48
2.7 Límites máximos de RI en líneas de transmisión.

Las compañías de servicios públicos canadienses, en colaboración con el

Departamento de Comunicaciones (DOC) de Canadá, establecieron límites para la
Interferencia de Radiofrecuencia (RFI) a partir de la investigación realizada en
HydroQuebec, según la Canadian Standards Association. Estos límites se detallan
en la Tabla 5.
Es importante destacar que los límites canadienses son ajenos a la fuente de la RFI,
lo que significa que son aplicables tanto a las descargas de corona como de brecha.
Canadá, al igual que Estados Unidos, enfrenta un problema similar en cuanto a la
principal causa de interferencia de radiofrecuencia (RFI). Las descargas de brecha
representan más del 90% de las reclamaciones relacionadas con RFI.
El Departamento de Comunicaciones, en colaboración con las compañías de
servicios públicos canadienses, invierte considerable tiempo y recursos en
identificar y resolver las fuentes de RFI. La tabla número 5 muestra los límites
máximos de RI establecidos por norma. [13].

Tabla 5. Límites máximos de condiciones climáticas favorables de la Asociación

Canadiense de Normas para campos de interferencia (0.15-30 MHz) [19].

Voltaje nominal de fase a fase Intensidad de campo de interferencia

[kV] [dBμ/m]
Debajo de 70 43
70 – 200 49
200 – 300 53
300 – 400 56
400 – 600 60
600 – 800 63

Se debe tener en cuenta que estos límites se basan en:

(1) Una frecuencia de 0.5 MHz.

(2) El medidor de tipo IEC/CISPR.
(3) El nivel máximo de RI que se produce durante el buen tiempo.
(4) Una distancia de 15 m de la fase exterior.

El medidor de tipo IEC/CISPR mide la corona alrededor de 2 dB menos que los

medidores de tipo ANSI más antiguos. Además, el nivel máximo de RI durante
condiciones climatológicas idóneas es aproximadamente 6 dB mayor que el nivel
medio de ruido de radio durante condiciones climatológicas favorables.

Desde una visión global, la IEC/CISPR ha tratado en varias ocasiones de establecer

límites de interferencia por radiofrecuencia (RI) para las líneas eléctricas aéreas. No
obstante, estos esfuerzos no han prosperado debido a que muchas naciones

49
consideran que la RI en este tipo de líneas es un asunto de interés nacional, más
que internacional. A pesar de ello, en 1986, la IEC/CISPR elaboró un "código de
prácticas para reducir la generación de ruido de radio" (IEC/CISPR 1986b). Aunque
algunos países cuentan con límites de RI, no es claro si estos límites están
actualizados. En un artículo publicado por la IEEE Power Engineering Society en
1980 (IEEE 1980), se informó que en países como Checoslovaquia (ahora dividida
en República Checa y República Eslovaca), la URSS (actualmente Rusia), Suiza y
Polonia existían regulaciones. Sin embargo, debido a la disolución de la URSS y
Checoslovaquia, no se tiene certeza sobre si estos límites siguen vigentes o si han
sido modificados o reemplazados [13].

La regulación suiza, promulgada en 1966, abarca una amplia variedad de

instalaciones y equipos eléctricos, incluidos los sistemas de alta tensión. Según esta
norma, el nivel de interferencia de radiofrecuencia (RI), medido a una frecuencia de
0,5 MHz ya una distancia de 20 metros de la fase externa en condiciones de tiempo
seco ya 10ºC, no debe exceder los 34 dB para líneas con menos de 100 kV, y los
46 dB para aquellas que abren por encima de dicho valor. De igual manera, la norma
polaca, establecida en 1969, cubre diversos tipos de equipos. Para las líneas de
transmisión, la intensidad del campo de interferencia, medida a 20 metros
lateralmente del conductor más externo durante su funcionamiento normal, no debe
superar los 750 µV/m (57,5 dB), siempre que la humedad relativa no sobrepase el
80 % , la temperatura sea superior a +5ºC, y la frecuencia esté dentro de los 500
kHz ±10 kHz [13].

En 1980, el IEEE (IEEE 1980) publicó una revisión exhaustiva que cubría todas las
consideraciones técnicas necesarias para definir los límites en líneas y estaciones
eléctricas aéreas. Una parte importante del contenido referente a la fijación de
límites en la banda de transmisión de amplitud modulada (AM) fue basada en los
estudios de Maruvada y Trinh, así como en trabajos previos del Subcomité de Ruido
de Radio del IEEE. [13].

50
Capítulo III
Marco
Técnico

51
En el presente capítulo se detalla el proceso de cálculo para determinar los perfiles
de interferencia electromagnética de cuatro líneas de transmisión. Inicialmente, se
calcula el campo eléctrico en la superficie de los conductores mediante simulaciones
utilizando el software COMSOL Multiphysics. Luego, se estiman los niveles de
interferencia electromagnética aplicando las ecuaciones () a (), implementadas en
un código desarrollado en MATLAB.
A continuación, se describen los cuatro casos de estudio analizados en este
capítulo:
Caso 1:
El primer caso se refiere a una línea de transmisión convencional de 500 kV en
configuración horizontal, con cuatro conductores por fase. Dicha línea de
transmisión se tomará como referencia para las configuraciones HSIL de los casos
2, 3 y 4.
Caso 2:
El segundo caso se refiere a una línea de transmisión HSIL convencional de 500 kV.
Partiendo de la configuración del caso 1, se realizó un estudio de optimización
previo para aumentar el SIL de la línea mediante la reconfiguración de las posiciones
de los conductores y aumentando el número de conductores de cuatro a seis por
fase, obteniéndose un haz simétrico, es decir, en forma circular y con la misma
distancia entre subconductores. Debido a lo anterior esta configuración se
denomina HSIL convencional.
Caso 3:
El tercer caso se refiere a una línea de transmisión HSIL no convencional de 500
kV. En este caso, se propone instalar 4 conductores por fase, dispuestos en un haz
asimétrico, es decir, no circular. Debido a lo anterior esta configuración se denomina
HSIL no convencional [18].
Caso 4:
El cuarto caso se refiere a una línea de transmisión HSIL no convencional de 500
kV. Finalmente, en este último caso se propone un haz de seis conductores, pero,
al igual que en el caso 3, los conductores de las fases laterales se colocan en una
configuración asimétrica.

52
3.1 Simulación del campo eléctrico en COMSOL MULTIPHYSICS
A continuación, se describen los pasos que se siguieron para simular el campo
eléctrico en los casos de estudio:
1. Inicializar COMSOL Multiphysics:
• Se requiere Iniciar el software COMSOL Multiphysics desde el acceso directo
en el escritorio o desde el menú de aplicaciones.

2. Configuración de un Nuevo Modelo:

• Seleccionar "Asistente de modeló" en la pantalla inicial, ya que es preferible
a comparación con “Modelo blanco” porque guía la configuración del modeló
paso a paso, asegurando la correcta selección de física, estudio, geometría,
materiales, condiciones de contorno y parámetros de solución, optimizando
así la precisión y eficiencia de la simulación.
• Elegir la dimensión del espacio de trabajo adecuado (3D, 2D asimétrico, 2D,
1D asimétrico,1D o 0D) conforme a las necesidades de la simulación, en
nuestro caso se opta por un sistema bidimensional (2D). Se lleva a cabo un
corte trasversal o axial de cada fase de la línea como se ilustra en la Figura
21.

Figura 21. Área de la sección transversal de las fases.

3. Selección del Módulo de Física (ver Figura 22):
• En el árbol de módulos de física, expandir el módulo "AC/DC".
• Expandir Campos eléctricos y corrientes
• Seleccionar "Electrostática (es)" para modelar el campo eléctrico estático,
esta interfaz además de calcular el campo eléctrico, calcula el campo de
desplazamiento eléctrico y el potencial en dieléctricos con una distribución
de carga prescrita. Utiliza la Ley de Gauss, resolviendo para el campo
eléctrico mediante el potencial escalar.
• Dar clic en "Agregar" para que esta aparezca en las interfaces de física
añadidas y finalmente seleccionar “Estudio”.
.

53
 Figura 22. Selección del Módulo de Física.
4. Definición del Estudio (ver Figura 23):
• Seleccionar el tipo de estudio adecuado, en este caso se opta por un estudio
“Estacionario”, ya que se busca analizar el comportamiento del campo
eléctrico en condiciones de estado estable, donde las variables del sistema
no cambian con el tiempo.
• Dar doble clic en “Estacionario” o un clic en “Hecho”.

Figura 23. Definición del estudio.

54
 • Tras seleccionar dicho estudio, como se presenta en las Figuras 24(a), 24(b),
24(c), 24(d), 24(e) y 24(f), se muestra la página inicial del programa.

Figura 24. (a) Barra de herramientas (b) Constructor de modeló (c) Barra de
ajustes (d) Barra de materiales (e) Gráficos (f) Barra de mensajes, estado de
cálculo, registro y tabla.
El análisis comienza con el caso 1, como se describe a continuación:
3.1.1 Caso 1: Línea de transmisión convencional con 4 subconductores por
haz, operando a 500 kV.

En la Tabla 6 se presentan los datos de los subconductores de la línea de

transmisión referida en este estudio como caso 1. La Figura 25 ilustra las distancias
y la distribución geométrica de la línea.
Tabla 6. Datos de la línea de transmisión convencional con 4 subconductores por
haz, operando a 500 Kv.

Caso 1 Número Radio [mm] Radio del Haz [m]

Subconductor 4 1.341 31.82

55
Figura 25. Distribución y dimensiones de la línea de transmisión convencional con
4 subconductores por haz, operando a 500 kV (Caso 1).
5. Posición de los subconductores
La posición de los subconductores en el software COMSOL Multiphysics se
determina especificando la altura, las distancias entre fases, el radio del haz y el
ángulo correspondiente en radianes, según se muestra en la en la Figura 26.

Figura 26. Círculo con su conversión de grados a radianes.

56
La ecuación 54 se emplea para determinar el desplazamiento angular de los
subconductores
360° (54)
𝐷𝑎𝑛𝑔 =
𝑁𝑠𝑢𝑏
Donde:
𝐷𝑎𝑛𝑔 = 𝐷𝑒𝑠𝑝𝑙𝑎𝑧𝑎𝑚𝑖𝑒𝑛𝑡𝑜 𝑎𝑛𝑔𝑢𝑙𝑎𝑟.

360° = 𝑅𝑎𝑑𝑖𝑜 𝑑𝑒𝑙 𝑐𝑖𝑟𝑐𝑢𝑙𝑜.

𝑁𝑠𝑢𝑏 = 𝑁ú𝑚𝑒𝑟𝑜 𝑑𝑒 𝑠𝑢𝑏𝑐𝑜𝑛𝑑𝑢𝑐𝑡𝑜𝑟𝑒𝑠.
Las ecuaciones 55, 56 y 57 se emplean para calcular la posición de los
subconductores en el eje X
𝑥𝑛 𝑓𝑎 = 𝑥1 + 𝑟𝑎𝑑𝑖𝑜ℎ𝑎𝑧 ∗ cos 𝜃 (55)
𝑥𝑛 𝑓𝑏 = 𝑥2 + 𝑟𝑎𝑑𝑖𝑜ℎ𝑎𝑧 ∗ cos 𝜃 (56)
𝑥𝑛 𝑓𝑐 = 𝑥3 + 𝑟𝑎𝑑𝑖𝑜ℎ𝑎𝑧 ∗ cos 𝜃 (57)

Donde:
𝑥𝑛 𝑓𝑎, 𝑥𝑛 𝑓𝑏, 𝑥𝑛 𝑓𝑐 = 𝑃𝑜𝑠𝑖𝑐𝑖𝑜𝑛 𝑒𝑛 𝑒𝑙 𝑒𝑗𝑒 𝑋 𝑑𝑒 𝑙𝑜𝑠 𝑠𝑢𝑏𝑐𝑜𝑛𝑑𝑢𝑐𝑡𝑜𝑟𝑒𝑠 𝑝𝑎𝑟𝑎 𝑙𝑎 𝑓𝑎𝑠𝑒 𝑎, 𝑏 𝑦 𝑐.
𝑥1, 𝑥2 𝑦 𝑥3 = 𝐷𝑖𝑠𝑡𝑎𝑛𝑐𝑖𝑎𝑠 𝑒𝑛𝑡𝑟𝑒 𝑓𝑎𝑠𝑒𝑠.
𝑟𝑎𝑑𝑖𝑜ℎ𝑎𝑧 = 𝑅𝑎𝑑𝑖𝑜 𝑑𝑒𝑙 ℎ𝑎𝑧 𝑑𝑒 𝑙𝑜𝑠 𝑠𝑢𝑏𝑐𝑜𝑛𝑑𝑢𝑐𝑡𝑜𝑟𝑒𝑠.

𝜃 = Á𝑛𝑔𝑢𝑙𝑜 𝑐𝑜𝑟𝑟𝑒𝑠𝑝𝑜𝑛𝑑𝑖𝑒𝑛𝑡𝑒 𝑒𝑛 𝑟𝑎𝑑𝑖𝑎𝑛𝑒𝑠

Las ecuaciones 58, 59 y 60 se emplean para calcular la posición de los
subconductores en el eje Y.
𝑦𝑛 𝑓𝑎 = ℎ𝑎 + 𝑟𝑎𝑑𝑖𝑜ℎ𝑎𝑧 ∗ sin 𝜃 (58)
𝑦𝑛 𝑓𝑏 = ℎ𝑏 + 𝑟𝑎𝑑𝑖𝑜ℎ𝑎𝑧 ∗ sin 𝜃 (59)
𝑦𝑛 𝑓𝑐 = ℎ𝑐 + 𝑟𝑎𝑑𝑖𝑜ℎ𝑎𝑧 ∗ sin 𝜃 (60)
Donde:
𝑦𝑛 𝑓𝑎, 𝑦𝑛 𝑓𝑏, 𝑦𝑛 𝑓𝑐 = 𝑃𝑜𝑠𝑖𝑐𝑖𝑜𝑛 𝑒𝑛 𝑒𝑙 𝑒𝑗𝑒 𝑌 𝑑𝑒 𝑙𝑜𝑠 𝑠𝑢𝑏𝑐𝑜𝑛𝑑𝑢𝑐𝑡𝑜𝑟𝑒𝑠 𝑝𝑎𝑟𝑎 𝑙𝑎 𝑓𝑎𝑠𝑒 𝑎, 𝑏 𝑦 𝑐.
ℎ𝑎, ℎ𝑏 𝑦 ℎ𝑐 = 𝐴𝑙𝑡𝑢𝑟𝑎 𝑑𝑒 𝑙𝑎𝑠 𝑓𝑎𝑠𝑒𝑠.
𝑟𝑎𝑑𝑖𝑜ℎ𝑎𝑧 = 𝑅𝑎𝑑𝑖𝑜 𝑑𝑒𝑙 ℎ𝑎𝑧 𝑑𝑒 𝑙𝑜𝑠 𝑠𝑢𝑏𝑐𝑜𝑛𝑑𝑢𝑐𝑡𝑜𝑟𝑒𝑠.

𝜃 = Á𝑛𝑔𝑢𝑙𝑜 𝑐𝑜𝑟𝑟𝑒𝑠𝑝𝑜𝑛𝑑𝑖𝑒𝑛𝑡𝑒 𝑒𝑛 𝑟𝑎𝑑𝑖𝑎𝑛𝑒𝑠

Ejemplo para la fase A del Caso 1:
Datos:
• ℎ𝑎 = 𝐴𝑙𝑡𝑢𝑟𝑎 = 28 [𝑚].

57
 • 𝑥1 = 𝑆𝑒𝑝𝑎𝑟𝑎𝑐𝑖ó𝑛 𝑒𝑛𝑡𝑟𝑒 𝑙𝑎𝑠 𝑓𝑎𝑠𝑒𝑠 𝐴 𝑦 𝐵 = −12.3 [𝑚]. 𝑐ó𝑚𝑜 𝑠𝑒 𝑖𝑙𝑢𝑠𝑡𝑟𝑎 𝑒𝑛 𝑙𝑎
𝑓𝑖𝑔𝑢𝑟𝑎 #.
• 𝑟𝑎𝑑𝑖𝑜ℎ𝑎𝑧 = 31.82 𝑐𝑚 = 0.3182 𝑚
• 𝜃=
𝜋 3𝜋 5𝜋 7𝜋
𝑐𝑜𝑛𝑑𝑢𝑐𝑡𝑜𝑟 1 (4 ) , 𝑐𝑜𝑛𝑑𝑢𝑐𝑡𝑜𝑟 2 ( 4 ) , 𝑐𝑜𝑛𝑑𝑢𝑐𝑡𝑜𝑟 3 ( 4 ) 𝑦 𝑐𝑜𝑛𝑑𝑢𝑐𝑡𝑜𝑟 4 ( 4 ),

Sustituyendo en la ecuación 54:

360° 360°
𝐷𝑎𝑛𝑔 = = = 90°
𝑁𝑠𝑢𝑏 4
Para este caso, los 90° se dividen entre 2 debido a la disposición de los
subconductores, como se muestra en la Figura 27.
90°
= 45°
2
Se suman 45° a cada 90°, convirtiendo los grados en radianes, como se muestra a
continuación:
𝜋
0° + 45° = 45° =
4
3𝜋
90° + 45° = 135° =
4
5𝜋
180° + 45° = 225° =
4
7𝜋
270° + 45° = 315° =
4

Figura 27. Posición de los cuatro subconductores por fase de la línea del caso 1.
58
Eje x: La posición de cada subconductor, se calcula sustituyendo los datos que
requiere la ecuación 55, de acuerdo a su fase correspondiente.
𝜋
𝑥1 𝑓𝑎 = −12.3 𝑚 + 0.3182 𝑚 ∗ cos ( ) = −12.075 𝑚
4
3∗𝜋
𝑥2 𝑓𝑎 = −12.3 𝑚 + 0.3182 𝑚 ∗ cos ( ) = −12.525 𝑚
4
5∗𝜋
𝑥3 𝑓𝑎 = −12.3 𝑚 + 0.3182 𝑚 ∗ cos ( ) = −12.525 𝑚
4
7∗𝜋
𝑥4 𝑓𝑎 = −12.3 𝑚 + 0.3182 𝑚 ∗ cos ( ) = −12.075 𝑚
4
Eje y: La posición de cada subconductor, se calcula sustituyendo los datos que
requiere la ecuación 58, de acuerdo a su fase correspondiente.
𝜋
𝑦1 𝑓𝑎 = 28 𝑚 + 0.3182 𝑚 ∗ cos ( ) = 28.225 𝑚
4
3∗𝜋
𝑦2 𝑓𝑎 = 28 𝑚 + 0.3182 𝑚 ∗ cos ( ) = 28.225 𝑚
4
5∗𝜋
𝑦3 𝑓𝑎 = 28 𝑚 + 0.3182 𝑚 ∗ cos ( ) = 27.775 𝑚
4
7∗𝜋
𝑦4 𝑓𝑎 = 28 𝑚 + 0.3182 𝑚 ∗ cos ( ) = 27.775 𝑚
4
Se aplica el mismo procedimiento para la fase B y C, así como para el caso 2, 3 y
4.
6. Parámetros 1
En la sección “Constructor de modelo”, como se muestra en la Figura 28(a), se
selecciona la opción “Parámetros 1”. En la Figura 28(b), bajo “Ajustes de
parámetros”, se configuran los parámetros de solución de la línea y se definen las
variables de salida pertinentes, según lo indica la Tabla 7.

59
 Figura 28. (a) Constructor de modelo “Parámetros 1”, (b) Ajustes de Parámetros
del caso 1.
Tabla 7. Parámetros del caso 1

7. Creación de la Geometría
Se selecciona la opción “Geometría 1” y al hacer clic derecho se despliega un menú
de herramientas de diseño geométrico para dibujar los conductores y otros

60
elementos del modeló, en cuanto a la forma se selecciona “Círculo”, ya que
representa un subconductor, como se muestra en la Figura 29(a) y de igual manera
se agregan 2 semicírculos, representando el domo de aire y el domo para el límite
de elementos infinitos.
Luego, se accede a “Ajustes de círculo”, como se ilustra en la Figura 29(b), donde
se ingresan los valores solicitados por el programa para definir las posiciones y
formas, los cuales ya fueron establecidos en la Tabla 7.
Los datos requeridos por el programa incluyen:
• Etiqueta: Nombre del subconductor o del domo.
• Tipo de objeto: Se selecciona “Sólido”.
• Tamaño y forma: Inserción del valor del radio y Ángulo de sector
(subconductores de 360°, domo de aire y domo para el límite de elementos
infinitos de 180°).
• Posición: Se establece la ubicación de cada subconductor en el plano (eje x,
y), donde el domo de aire y el domo para el límite de elementos infinitos se
colocan en x=0, y=0, con base en el centro.
Existen 2 métodos para ingresar la posición de cada subconductor:
1) En la Tabla 7 se nombran las variables para la posición en los ejes X e Y.
• x1fa: donde “x” indica el eje X, “1” al número del subconductor y “fa” la fase
A.
• y1fa: donde “y” indica el eje Y, “1” al número del subconductor y “fa” la fase
A.
Este mismo procedimiento se aplica para posicionar los demás subconductores.
2) Alternativamente, los valores de las posiciones en los ejes X e Y pueden ser
ingresados directamente si ya se conocen.
Para avanzar con la simulación de la línea de transmisión del Caso 1, se procede a
construir los siguientes elementos, aclarando que se insertan los valores que nos
solicita la ventanilla de “Ajustes de Círculo” de la Figura 29(b).
Nota: Los parametros que no se mencionen se quedan por defecto.
• Circulo 1
- Etiqueta: Domo de aire
- Radio: radioaire
- Ángulo de sector: 180°
La Figura 30(a), que se presenta mas adelante, ilustra las características
mencionadas anteriormente la cual corresponde al domo de aire.
• Circulo 2

61
 - Etiqueta: Domo de elementos infinitos
- Radio: radioeleminf
- Ángulo de sector: 180°
La Figura 30(b),ilustra las características mencionadas anteriormente la cual
corresponde al domo para el límite de elementos infinitos.
En la Figura 30(c), se ilustran los Circulos 3 a 14. Estos circulos representan los
subconductores de las Fases A, B y C, con las caracteristicas mencionadas
seguidamente en la Tabla 8.
En particular, los subconductores de la Fase A estan representados por los Circulos
3 a 6, los de la Fase B por los Círculos 7 a 10, y los de la Fase C por los Círculos
11 a 14.
Tabla 8. Características de los círculos que representan los subconductores de las
Fases A, B y C.

Circulo Etiqueta Radio Ángulo de x y

sector
3 Subconductor radioscond 360° x1fa y1fa
1-Fase A
4 Subconductor radioscond 360° x2fa y2fa
2-Fase A
5 Subconductor radioscond 360° x3fa y3fa
3-Fase A
6 Subconductor radioscond 360° x4fa y4fa
4-Fase A
7 Subconductor radioscond 360° x1fb y1fb
1-Fase B
8 Subconductor radioscond 360° x2fb y2fb
2-Fase B
9 Subconductor radioscond 360° x3fb y3fb
3-Fase B
10 Subconductor radioscond 360° x4fb y4fb
4-Fase B
11 Subconductor radioscond 360° x1fc y1fc
1-Fase C
12 Subconductor radioscond 360° x2fc y2fc
2-Fase C
13 Subconductor radioscond 360° x3fc y3fc
3-Fase C
14 Subconductor radioscond 360° x4fc y4fc
4-Fase C

62
Figura 29. (a)Constructor de modelo, (b) Ajustes del círculo para el domo de aire.

Figura 30. (a) Domo de Aire, (b) Domo para el límite de elementos infinitos, (c)
Subconductores de las Fases A, B y C.

63
El siguiente paso es convertir los conductores en huecos para optimizar el uso de
memoria en la computadora. Como se muestra en la Figura 29(a), se selecciona
"Geometría 1" y, al hacer clic derecho, se elige la opción "Booleanos y particiones".
Al seleccionar esta opción, aparece una ventana donde se elige "Diferencia 1". Para
añadir elementos, se utiliza la opción "Seleccionar casilla", ubicada en la parte
superior del gráfico, como se muestra en la Figura 31.
El procedimiento para llenar "Diferencia 1" es el siguiente:
• Objetos a añadir: Se añade el domo de aire (c1) y el domo para el límite de
elementos finitos (c2), los cuales se resaltan en color morado cuando son
seleccionados, como se muestra en la Figura 32.
• Objetos a sustraer: Se incorporan todos los subconductores del Caso 1: los
4 subconductores de la fase A (de c3 a c6), los 4 subconductores de la fase
B (de c7 a c10) y los 4 subconductores de la fase C (de c11 a c14), como se
ilustra en la Figura 33.

Figura 31. Seleccionar casilla.

Figura 32. Añadido del domo de aire (c2) y del domo para el límite de elementos
infinitos (c1).

64
 Figura 33. Sustracción para que los conductores de las fases A, B y C queden
huecos.
8. Definición del Dominio del Modelo:
Como se muestra en la Figura 29(a), se selecciona la sección "Definiciones"
haciendo clic derecho, como se muestra en la Figura 34(a), se elige la opción
"Dominio de elementos infinitos", lo que abre la configuración correspondiente.
Usando el seleccionador de casillas, se selecciona el domo para el límite de
elementos infinitos, que se resalta en color morado. En los ajustes, dentro del
apartado de geometría, se selecciona el tipo cilíndrico, como se ilustra en la Figura
34(b).

Figura 34. (a) Funciones de definiciones, (b) Ajustes del dominio de elemento
finito.
65
 9. Asignación de Materiales
Como se muestra en la Figura 29(a), se selecciona el apartado "Materiales"
haciendo clic derecho y se añade un material desde la biblioteca. Esto abre un
recuadro, como se muestra del lado izquierdo en la Figura 35. En la sección
"Básica", se hace clic para desplegar varios tipos de materiales, como aire y
aluminio. Se selecciona el material "Aire" y se asigna tanto al domo de aire (2) como
al domo para el límite de elementos infinitos (1).

Figura 35. Selección del material “Aire”.

10. Configuración de Condiciones de Contorno o fronteras
En la sección de electroestática, se definen las condiciones de contorno para
superficies y bordes debido a la resolución de ecuaciones diferenciales.
• Elementos a tierra: Como se muestra en la Figura 29(a) se selecciona la
opción “Electroestática (es)”, dando clic derecho se elige “Tierra”, como se
muestra en la Figura 36 se seleccionan las fronteras a tierra con el
seleccionador de casillas, las cuales se denominan con el 2 y 3.

Figura 36. Ajuste de las condiciones de contorno para la “Tierra”.

66
Potencial eléctrico: como se muestra en la Figura 29(a), se selecciona el apartado
de “Electroestática (es)”, dando clic derecho se elige “Potencial eléctrico”. Como se
muestra en la Figura 37, primero se seleccionan todos los subconductores de la
fase A utilizando el seleccionador de contorno, y se asigna el nombre de la variable
“vfasea” al potencial eléctrico. Luego, se repite el proceso para la fase B, asignando
“vfaseb”, y para la fase C, asignando “vfasec”.

Figura 37. Asignación del potencial eléctrico a cada subconductor de acuerdo a su

fase (A, B o C).
11. Generación de la malla
Para construir el mallado, como se muestra en la Figura 29(a), se selecciona el
apartado "Malla 1". En la opción de tipo de secuencia, se elige "Malla controlada por
el usuario". Dentro de esta sección, se selecciona la opción "Más fina" para crear
un mallado más detallado y garantizar mayor precisión en la simulación. Al hacer
clic derecho en "Malla 1", se selecciona "Triangular libre", lo que crea un nuevo
apartado llamado "Triangular libre 1". Al hacer clic derecho sobre este, se selecciona
la opción "Tamaño 1".
Configuración de tamaño 1:
a) Etiqueta: Se deja como predeterminado.
b) Selección de entidades geométricas:
• Nivel de entidad geométrico: Contorno.
• Selección: Manual. Usando el seleccionador de casillas, se seleccionan
todos los subconductores de cada fase.
c) Tamaño de elemento:
• Calibrar para: Física general (predeterminado).
• Seleccionar “Definido por el usuario”.

67
 d) Parámetros de tamaño de elemento
• Se selecciona la opción “Tamaño máximo de elemento” y se establece un
valor de 0.00005 m.
El mallado se refiere a la discretización del dominio de solución, donde el método
de elementos finitos divide el dominio en pequeños elementos, sobre los cuales se
definen ecuaciones. Las incógnitas son los potenciales en los vértices de los
triángulos que forman el mallado.
Al hacer clic en "Construir", se puede observar el mallado de la simulación en la
Figura 38. Además, en la Figura 39(a), se muestra un acercamiento del mallado en
la superficie de los subconductores de la fase A, y en la Figura 39(b), se presenta el
mallado de la superficie del subconductor 1 de la fase A.

Figura 38. Mallado

Figura 39. (a) Mallado en la superficie de los subconductores de la fase A, (b)

Mallado en la superficie del conductor 1 de la fase A

68
 12. Ejecución de la Simulación
Como se muestra en la Figura 40 se selecciona el apartado de “Estudio 1” y se hace
clic en “Calcular”. Esto inicia la simulación en el programa, que procede a calcular
el campo eléctrico con base en las condiciones y la geometría previamente
definidas.

Figura 40. Ejecución de la simulación.

13. Análisis de resultados
Al finalizar la simulación, primeramente, se selecciona el apartado de "Resultados
“en donde se desactiva la opción "Potencial eléctrico (es)", como se muestra en la
Figura 41.
A continuación, en el apartado "Campo eléctrico, norma (es)", se hace clic derecho
y se agrega la opción "Curva de nivel 1", la cual permite visualizar la distribución del
campo eléctrico, tal como se observa en la Figuras 42 y en la figura 43 para el haz
de la fase A.

Figura 41. Análisis de resultados.

69
Figura 42. Distribución del campo eléctrico en la superficie de los subconductores
del caso 1.

Figura 43. Distribución del campo eléctrico en la superficie de los subconductores

de la fase A del caso 1.
Posteriormente, como se muestra en la Figura 44(a), se selecciona la opción
"Definiciones" haciendo clic derecho, lo que abre una ventana donde se elige la
opción "Sondas". En la siguiente ventana, se selecciona "Sonda de contorno". En la
ventana de "Ajustes de sonda de contorno" como se muestra en la Figura 44(b), se
introducen los parámetros requeridos por el programa:
70
 • Etiqueta: EmaxFaseA.
• Nombre de variable: Se deja el valor predeterminado.
Tipo de sonda:
• Tipo: Máximo.
Selección de origen:
• Selección: Manual.
• Con el seleccionador de casillas, se seleccionan todos los subconductores
de la fase A.
Expresión:
• Expresión: [Link] (Presionando Ctrl+Espacio, en el apartado
"Componente 1", se selecciona "Electrostática", luego "Eléctrico", y
finalmente "[Link] - Campo eléctrico, norma - V/m").
• Unidad de tabla y gráfico: kV/cm.
Para las Fases A y B se sigue el mismo procedimiento anterior o de igual manera
se duplica la primera sonda de medición y se le cambia la etiqueta y la selección de
los subconductores con base a su fase correspondiente.

Figura 44. (a) Selección de la sonda de medición, (b) Ajuste de la sonda de

contorno con respecto a cada fase.
Finalmente, como se muestra en la parte superior izquierda de la Figura 44(b),
damos clic en "Actualizar resultados". Luego, en la parte inferior derecha de la
pantalla, aparecerá una ventana con el resultado obtenido, como se muestra en la
Figura 45.

71
 Figura 45. Resultado del Campo eléctrico para la Fase A.
Este procedimiento debe repetirse para cada sonda de medición, ya que cada una
corresponde a su respectiva fase.
En la Tabla 9 se presentan los resultados del campo eléctrico máximo obtenidos en
kV/cm para cada fase. Este valor es fundamental para calcular la interferencia
electromagnética generada por una línea de transmisión.
Tabla 9. Resultados del campo eléctrico máximo para el caso 1.

EmaxFaseA EmaxFaseB EmaxFaseC

14.143 kV/cm 15.364 kV/cm 14.140 kV/cm

3.1.2 Caso 2: Línea de transmisión de alta cargabilidad (HSIL) convencional

con 6 subconductores por haz, operando a 500 kV.
Este caso corresponde a una línea de transmisión de alta cargabilidad (HSIL) con
una configuración convencional de 6 subconductores por haz, operando a 500 kV,
como se ilustra en la Figura 46. En la Tabla 10 se detallan todos los parámetros
necesarios para su simulación empleando el software COMSOL Multiphysics,
siguiendo el mismo procedimiento descrito en el caso 1.

Figura 46. Distribución y dimensiones de la línea de transmisión de alta

cargabilidad (HSIL) convencional con 6 subconductores por haz, operando a 500
kV (Caso 2).
72
 Tabla 10. Parámetros del caso 2.

En la Figura 47 se muestra la distribución del campo eléctrico en la superficie de los

subconductores correspondientes al caso y en la Figura 48 para el haz de la fase A.

73
Figura 47. Distribución del campo eléctrico en la superficie de los subconductores
del caso 2.

Figura 48. Distribución del campo eléctrico en la superficie de los subconductores

de la fase A del caso 2.
En la Tabla 11 se presentan los resultados del campo eléctrico máximo obtenidos
en kV/cm para cada fase.
Tabla 11. Resultados del campo eléctrico máximo para el caso 2.

EmaxFaseA EmaxFaseB EmaxFaseC

20.290 kV/cm 20.527 kV/cm 20.293 kV/cm

74
3.1.3 Caso 3: Línea de transmisión de alta cargabilidad (HSIL) no
convencional con 4 subconductores por haz, operando a 500 kV.
Este caso corresponde a una línea de transmisión de alta cargabilidad (HSIL) con
una configuración no convencional de 4 subconductores por haz, operando a 500
kV. En las fases laterales de esta línea, los subconductores se encuentran
dispuestos de manera asimétrica, como se ilustra en la Figura 49. La Tabla 12
detalla los parámetros necesarios para su simulación empleando el software
COMSOL Multiphysics, siguiendo el mismo procedimiento descrito en el caso 1.

Figura 49. Distribución y dimensiones de la línea de transmisión de alta

cargabilidad (HSIL) convencional con 6 subconductores por haz, operando a 500
kV (Caso 2).
Tabla 12. Parámetros del caso 3.

75
En la Figura 50 se muestra la distribución del campo eléctrico en la superficie de los
subconductores correspondientes al caso y en la Figura 51 para el haz de la fase A.

Figura 50. Distribución del campo eléctrico en la superficie de los subconductores

del caso 3.

Figura 51. Distribución del campo eléctrico en la superficie de los subconductores

de la fase A del caso 3.
En la Tabla 13 se presentan los resultados del campo eléctrico máximo obtenidos
en kV/cm para cada fase.
Tabla 13. Resultados del campo eléctrico máximo para el caso 3

EmaxFaseA EmaxFaseB EmaxFaseC

17.249 kV/cm 17.571 kV/cm 17.248 kV/cm

76
3.1.4 Caso 4: Línea de transmisión de alta cargabilidad (HSIL) no
convencional con 6 subconductores por haz, operando a 500 kV.
Este caso corresponde a una línea de transmisión de alta cargabilidad (HSIL) con
una configuración no convencional de 6 subconductores por haz, operando a 500
kV. En las fases de esta línea, los subconductores se encuentran dispuestos de
manera asimétrica, como se ilustra en la Figura 52. La Tabla 14 detalla los
parámetros necesarios para su simulación empleando el software COMSOL
Multiphysics, siguiendo el mismo procedimiento descrito en el caso 1.

Figura 52. Distribución y dimensiones de la línea de transmisión de alta

cargabilidad (HSIL) no convencional con 6 subconductores por haz, operando a
500 kV (Caso 4)
Tabla 14. Parámetros del caso 4.

77
En la Figura 53 se muestra la distribución del campo eléctrico en la superficie de los
subconductores correspondientes al caso y en la Figura 54 para el haz de la fase A.

Figura 53. Distribución del campo eléctrico en la superficie de los subconductores

del caso 4.

Figura 54. Distribución del campo eléctrico en la superficie de los subconductores

de la fase A del caso 4.
En la Tabla 15 se presentan los resultados del campo eléctrico máximo obtenidos
en kV/cm para cada fase.
Tabla 15. Resultados del campo eléctrico máximo para el caso 4.

EmaxFaseA EmaxFaseB EmaxFaseC

17.592 kV/cm 17.674 kV/cm 17.591 kV/cm

78
3.2 Cálculo de Radio Interferencia
3.2.1 CIGREE
En la Tabla 16 se muestran los resultados de radio interferencia para las 4 líneas de
transmisión (Casos) y para cada tipo de ambiente, donde dichos perfiles tienen
valores máximos de RI de:
Tabla 16. Valores obtenidos de RI con la ecuación CIGREE para los 4 casos de
estudio y para cada tipo de ambiente analizado.

Perfiles de RI
Caso Buen ambiente Lluvia moderada Lluvia densa
dB dB dB
1 38.217652 55.217652 62.217652
2 53.484543 70.484543 77.484543
3 46.481588 63.481588 70.481588
4 43.970517 60.970517 67.970517

En las Figura 55, 56 y 57 se muestran los perfiles de radio interferencia de los 4

casos de estudio, se pueden apreciar las discrepancias de acuerdo a los niveles de
RI entre cada caso de estudio y para cada tipo de ambiente analizado.

Figura 55. Perfiles de radio interferencia de los 4 casos de estudio en Buen

ambiente (CIGREE).

79
 Figura 56. Perfiles de radio interferencia de los 4 casos de estudio en Lluvia
moderada (CIGREE).

Figura 57. Perfiles de radio interferencia de los 4 casos de estudio en Lluvia densa
(CIGREE).

80
3.2.2 BPA
En la Tabla 17 se muestran los resultados de radio interferencia para las 4 líneas de
transmisión (Casos) y para cada tipo de ambiente, donde dichos perfiles tienen
valores máximos de RI de:
Tabla 17. Valores obtenidos de RI con la ecuación BPA para los 4 casos de
estudio y para cada tipo de ambiente analizado.

Perfiles de RI
Caso Buen ambiente Lluvia moderada Lluvia densa
dB dB dB
1 47.938796 64.938796 71.938796
2 59.825357 76.825357 83.825357
3 55.490979 72.490979 79.490979
4 52.326711 69.326711 76.326711

En las Figura 58, 59 y 60 se muestran los perfiles de radio interferencia de los 4

casos de estudio, se pueden apreciar las discrepancias de acuerdo a los niveles de
RI entre cada caso de estudio y para cada tipo de ambiente analizado.

Figura 58. Perfiles de radio interferencia de los 4 casos de estudio en Buen

ambiente (BPA).

81
 Perfiles de radio interferencia de los 4 casos de estudio en Lluvia moderada
(BPA).

Figura 59. Perfiles de radio interferencia de los 4 casos de estudio en Lluvia densa
(BPA).

82
3.2.3 CREIPI
En la Tabla 18 se muestran los resultados de radio interferencia para las 4 líneas de
transmisión (Casos) y para cada tipo de ambiente, donde dichos perfiles tienen
valores máximos de RI de:
Tabla 18. Valores obtenidos de RI con la ecuación CREIPI para los 4 casos de
estudio y para cada tipo de ambiente analizado.

Perfiles de RI
Caso Buen ambiente Lluvia moderada Lluvia densa
dB dB dB
1 48.155709 65.155709 72.155709
2 63.843609 80.843609 87.843609
3 56.745161 73.745161 80.745161
4 53.674351 70.674351 77.674351

En las Figura 60, 61 y 62 se muestran los perfiles de radio interferencia de los 4

casos de estudio, se pueden apreciar las discrepancias de acuerdo a los niveles de
RI entre cada caso de estudio y para cada tipo de ambiente analizado.

Figura 60. Perfiles de radio interferencia de los 4 casos de estudio en Buen

ambiente (CREIPI).

83
 Figura 61. Perfiles de radio interferencia de los 4 casos de estudio en Lluvia
moderada (CREIPI).

Figura 62. Perfiles de radio interferencia de los 4 casos de estudio en Lluvia densa
(CREIPI).

84
Capítulo IV

85
4.1 Análisis económico
En esta parte se aborda el componente económico del proyecto. Dado que se trata
de una iniciativa conceptual, el análisis económico incluye los equipos, materiales,
insumos y la mano de obra necesaria requerida para la elaboración de este
proyecto.
En términos generales, los factores a tener en cuenta son: el costo del servicio
eléctrico, el servicio de internet, el equipo informático, la licencia de los softwares
MATLAB® y COMSOL Multiphysics®, así como el sueldo de un Ingeniero
Electricista.
• Gasto de energía eléctrica [1].
Para llevar a cabo este trabajo, se dedicaron alrededor de 3 horas diarias por
integrante, y el consumo de una computadora portátil es aproximadamente de 200
W. Esto resultaría en un consumo bimestral de 73 kWh.
Dado que no se supera el límite de 75 kWh, se aplicaría la tarifa básica de $0.902
MXN por cada kWh, lo que equivale a $65.76 MXN por bimestre, por lo que el total
de facturación de energía eléctrica en un lapso de 12 meses (tiempo en el que se
desarrolló el proyecto) es de $394.53 MXN.
• Servicio de internet.
El proyecto se llevó a cabo en aproximadamente 12 meses. El servicio de internet
contratado tiene un costo mensual de $350.00 MXN, por lo que el gasto total de
internet para el desarrollo de este proyecto es de $4,200.00 MXN.
• Equipo de cómputo:
Para llevar a cabo este proyecto, fue esencial contar con un equipo computo por
cada integrante, debido a que se emplearon métodos de cálculo y simulaciones en
los softwares MATLAB® y COMSOL Multiphysics® respectivamente, el monto total
de los 3 equipos es de $27,499.00 MXN. .
• Salario de un Ingeniero Electricista [2].
El salario mensual promedio de un Ingeniero Electricista en México es de
$11,000.00 MXN, por lo que el monto total durante los 12 meses es de $396,00.00
MXN [2].
• Licencia de MATLAB® [3].
Una parte significativa del desarrollo de este proyecto se llevó a cabo en MATLAB®,
utilizando la licencia proporcionada por el instituto. Sin embargo, si fuera necesario
adquirir una licencia, es importante tener en cuenta que la tienda de MathWorks
ofrece varias opciones (en algunas de las cuales la herramienta de optimización
“Optimization Toolbox” tiene un costo adicional, ya incluida en la lista de precios que
se presenta en la Tabla 19 a continuación):
86
 Tabla 19. Precios de las licencias de MATLAB® [3].

Licencia Tipo de licencia Precio en USD Precio en MXN

Licencia Anual 980.00 19,532.67
estándar/individual Perpetua 2,450.00 48,831.69
Licencia de uso Perpetua 149.00 2,969.76
Doméstico/individual
Licencia de uso Anual 275.00 5,481.11
académico/individual Perpetua 550.00 10,962.22
Licencia de Estándar 55.00 1,096.22
Estudiante (perpetua)
Suit (perpetua) 99.00 1,973.20

Tipo de cambio 1 USD = $ 19.9313 MXN (21/10/2024) [4].

NOTA: Los precios de la Tabla 4.4 no incluyen impuestos.
Al comparar las licencias de MATLAB® de la Tabla 19, se puede observar que la
opción más económica es la licencia de estudiante estándar. Sin embargo, resulta
más conveniente adquirir la licencia suite, ya que incluye diversas herramientas
adicionales con un costo de $1,973.00 MXN.
• Licencia de COMSOL Multiphysics® [5].
Una parte fundamental en el desarrollo de este proyecto fue el empleo del software
COMSOL para la elaboración de simulaciones, la licencia fue proporcionada por el
asesor técnico.

Tabla 20. Precios de las licencias de COMSOL Multiphysics® [5].

Precio en Precio en
Licencia Tipo de licencia
USD MXN
Licencia Académica de Anual 998.00 19,891.44
COMSOL Multiphysics® Perpetua 1,995.00 39,762.94
Licencia de usuario único Anual 998.00 19,891.44
designada Perpetua 1,995.00 39,762.94
Anual 1,995.00 39,762.94
Licencia de red flotante
Perpetua 3,990.00 79,525.89
Anual 1,998.00 39,822.74
Licencia de COMSOL Server
Perpetua 3,995.00 79,625.54

Tipo de cambio 1 USD = $ 19.9313 MXN (21/10/2024) [4].

87
 NOTA: Los precios de la Tabla 4.4 no incluyen impuestos.
La licencia utilizada durante el desarrollo del proyecto fue la licencia académica de
COMSOL Multiphysics®, con un periodo de duración de un año por lo que el monto
de la licencia de acuerdo con la Tabla 20 es de $19,891.44 MXN.
A continuación, en la Tabla 21, se realiza un estimado del costo final del proyecto.

Tabla 21. Estimado del costo final del proyecto.

Insumos Costo unitario MXN Total, MXN

Gasto de energía eléctrica 65.76 (bimestral) 394.53
Internet 350.00 (mensual) 4,200.00
Equipo de cómputo 27,499.00 (pago único) 27,499.00
Licencia suite de estudiante de
1,973.20 (pago único) 1,973.20
MATLAB®
Salario de un Ingeniero
11,000.00 (mensual) 396,000.00
Electricista
Licencia Académica de COMSOL
19,891.44 (anual) 19,891.44
Multiphysics®
Total MXN 449,958.17

El costo total del proyecto ascendería a $449,892.42 MXN lo que lo convierte en

una opción viable. Los programas desarrollados pueden ser utilizados múltiples
veces, tanto para calcular la RI en líneas de transmisión convencionales como
líneas de alta cargabilidad. Esto podría permitir conocer los perfiles de RI en las
líneas ya existentes como en líneas próximas a implementar.
4.1.1 Costos unitarios

En la Taba 22 se pueden apreciar los costos unitarios del proyecto.

Se puede notar que los costos se dividen en servicios, salarios de mano de obra y
el costo de los equipos de cómputo utilizados, donde estos últimos se registran
como pago único.

88
 Tabla 22. Costos unitarios del proyecto.

Costo
Insumos Unidad Cantidad Total MXN
unitario MXN
Gasto de
energía Bimestre 6 65.76 394.53
eléctrica
Internet Mes 12 350.00 4,200.00
Equipo de
- 1 10,000.00 10,000.00
cómputo 1
Equipo de
- 1 8,999.00 8,999.00
cómputo 2
Equipo de
- 1 8,500.00 8,500.00
cómputo 3
Licencia suite
de estudiante - 1 1,973.20 1,973.20
de MATLAB®
Licencia
Académica de
Año 1 19,891.44 19,891.44
COMSOL
Multiphysics®
Salario de un
Ingeniero Mes 12 11,000.00 132,000.00
Electricista (1)
Salario de un
Ingeniero Mes 12 11,000.00 132,000.00
Electricista (2)
Salario de un
Ingeniero Mes 12 11,000.00 132,000.00
Electricista (3)
Total MXN 449,958.17

89
4.2 Conclusiones

90
4.3 Recomendaciones

91
4.4 Referencias.
[1] Enríquez Harper, G. (2008). Sistemas de Transmisión y Distribución de Potencia
Eléctrica. Ciudad de México, México: Limusa.

[2] Francisco M. Gonzalez-Longatt. (2007). Elementos de líneas de transmisión aéreas.

[3] Victor F. Nasimba Medina. (2020). Introducción a las líneas de transmisión.

[4] Efecto corona sobre lineas de transmision, CITE energia, Lima/Silicon Technology

[5] L. M. Checa, Líneas de Transporte de Energía, 3ra edición, Editorial Marcombo,

Barcelona España, 1988

[6] Redes subterránea de distribución. (1981).

[7] CIGRE (1974) "Interferences produced by corona effect electric systems (Description of
Phenomena and Practical Guide for Calculation).", CIGRE Working Grop 36.01. EMC
Aspects of Corona, Electric and Magnetic Fields.

[8] L. A. Siegert, Alta Tensión y Sistemas de Transmisión. Editorial Limusa, 1988.

[9] H. N. Carrazco Zuñiga (2023), ´´Efecto de las sobretensiones en el diseño de líneas

de transmisión de alta tensión mayores a 300 kV corriente alterna´´. Perú.

[10] C. Tejada, (2008), “Análisis de la Interferencia Electromagnética por Efecto

Corona en Líneas de Transmisión de Alta Tensión en CA”. Ciudad de México:
Instituto Politécnico Nacional.

[11] B. V. Mena Arias (2008) ´´Efecto corona en las líneas de transmisión de 500 kV´´. Quito
Escuela Politécnica Nacional.

[12] E. Montoya (2008). ´´Estudio estocástico de cargabilidad en líneas de transmisión´´,

Ciudad de México: Instituto Politécnico Nacional.

[13] M. Salazar (2013). “Alternativas de solución para el incremento de cargabilidad en

líneas de transmisión del sistema eléctrico nacional “. Ciudad de México: Universidad
Nacional Autónoma de México.

[14] ARTI SIL

92
[15] N. Roman, G. Arroyo, “Evaluación de cargabilidad en líneas de transmisión”, 4to.
Coloquio Académico de la Sección de Estudios de Postgrado e Investigación de la ESIME-
IPN, 25-27 Noviembre 1992, México, D.F.

[16] Miller, J. E. Timothy, “Reactive power control in electric systems”, New York, John Wiley,
1982.

[17] Stevenson, W. “Análisis de sistemas eléctricos de potencia”, 2da edición, Editorial

Libros McGRAW-HILL DE MEXICO S.A de C.V, 1985
[18] Henan Guowang Cable Group Co., LTD. Poder de yisainuo, Ciudad de Kaifeng,
provincia de Henan

[19] NORMA Oficial Mexicana NOM-001-SEDE-2012, Instalaciones Eléctricas (utilización)

[20] Derecho de Vía, Especificación CFE L 1000-10

[21] EPRI Project Manager (2005). "EPRI AC Transmission Line Reference Book - 200 kV
and Above, Third Edition", Palo Alto, California.

[22] P. Sarma Maruvada. (2000) Corona Performance of High-Voltage Transmission Lines,

Research Studies Press Ltd., England.

93
4.5 Glosario
Derecho de vía: es una franja de terreno que se ubica a lo largo de cada línea, cuyo
eje longitudinal coincide con el trazo topográfico de la línea. Su dimensión
transversal varía de acuerdo con el tipo de estructuras, con la magnitud y
desplazamiento lateral de la flecha y con la tensión de operación [16].
Sistema de transmisión: sistema de suministro eléctrico constituida por los
elementos necesarios para llevar hasta los puntos de consumo y a través de
grandes distancias, la energía eléctrica generada en las centrales eléctricas [2].
Línea de alta tensión: es una estructura metálica que permite la transmisión y
distribución de energía eléctrica a grandes distancias [3].
Líneas de alta cargabilidad (HSIL): líneas de transmisión convencionales y no
convencionales con configuraciones simétricas o asimétricas en sus conductores
con el fin de incrementar la cargabilidad de la línea [4].
Efecto corona: fenómeno que ocurre cuando en la superficie de los conductores
de una línea de transmisión se presenta una intensidad de campo eléctrico mayor
a la rigidez dieléctrica del aire, por lo que este se ioniza ocasionando descarga
luminosa, ruido audible, radio interferencia y distorsión en la propagación de ondas
[5].
SIL: la máxima capacidad de carga que soporta la línea sin experimentar
sobretensiones perjudiciales. Representa el punto máximo de transferencia de
potencia donde la impedancia característica de la línea se iguala a la carga
conectada en los extremos, es decir que el SIL representa la carga en MW de la
línea donde se obtiene el equilibrio de potencia reactiva natural [6].
Diámetro del conductor: el diámetro del conductor define el tamaño del conductor
utilizado en la línea de transmisión [6].
Conductores en haz: agrupaciones de dos o más subconductores que conforman
la configuración de cada una de las fases [7].
Permitividad: densidad de flujo eléctrico dividido por la intensidad de campo
eléctrico [8].
Intensidad de campo eléctrico: fuerza ejercida por un campo eléctrico sobre una
carga eléctrica puntual, dividida por el valor de la carga [8].
Capacitancia: carga dividida por la diferencia de potencial eléctrico.
Tensión eléctrica: la tensión entre dos puntos 1 y 2 es la integral de línea desde el
punto 1 hasta el punto 2 de la intensidad de campo eléctrico [8].
Corriente eléctrica: es la intensidad de una corriente constante que mantenida en
dos conductores paralelos rectilíneos de longitud infinita, cuya área de sección
circular es despreciable, colocados a un metro de distancia entre sí, en el vacío,

94
producirá entre estos conductores una fuerza igual a 2x10-7 newton por metro de
longitud [8].
Inducción magnética: elemento de corriente, es igual al producto vectorial de este
elemento y la densidad de flujo magnético [8].
Resistencia: la diferencia de potencial eléctrico dividido por la corriente, cuando no
haya fuerza electromotriz en el conductor [8].
Impedancia: la representación compleja de la diferencia de potencial, dividida por
la representación compleja de la corriente [8].
Reactancia: parte imaginaria de la impedancia [8].
Interferencia electromagnética: la interferencia electromagnética (EMI) es
una perturbación o ruido no deseado que afecta a los circuitos y rutas
eléctricas. También se le conoce como interferencia de radiofrecuencia (RFI). La
EMI puede provocar que los dispositivos electrónicos funcionen mal, presenten
fallas o dejen de funcionar por completo [9].
Herraje: contratuercas, pasacables (monitor) u otra parte de un sistema de
alambrado, destinado principalmente para desempeñar una función más mecánica,
que eléctrica [1].
Campo eléctrico: se define como la región del espacio donde una carga eléctrica
experimenta una fuerza debido a la presencia de otra carga [8].
Sistema eléctrico de potencia: el sistema eléctrico de potencia es un conjunto de
elementos que tiene como fin generar, transformar, transmitir, distribuir y consumir
la energía eléctrica de tal forma que se logre la mayor calidad al menor costo posible
[10].
Líneas de Distribución o Transmisión Aérea: es aquélla que está constituida por
conductores desnudos o aislados y circuitos, tendidos en espacios abiertos y que
están soportados por estructuras (postes o torres), con los accesorios necesarios
para la fijación, separación y aislamiento de los mismos conductores. [16].
Línea subterránea: línea que está constituida por uno o varios cables aislados que
forman parte de un circuito eléctrico o de comunicación, colocados bajo el nivel del
piso, ya sea directamente enterrados, en ductos o en cualquier otro tipo de
canalización aprobado [1].
Tensión nominal de un sistema: valor de la tensión eléctrica que se asigna a un
sistema, parte de un sistema, un equipo o cualquier otro elemento ya la cual se
refieren ciertas características de operación y/o comportamiento de éstos [11].

95
Frecuencia del sistema: número de ciclos o alternancias por segundo que efectúa
la tensión eléctrica de un sistema eléctrico; se expresa en hercios [11].
Tensiones eléctricas normalizadas: valores de tensión eléctrica preponderantes
para finos de lograr un desempeño óptimo de los sistemas eléctricos y de los
elementos que la componen, con el fin de satisfacer las necesidades de los usuarios
[11].
Aisladores: son los elementos que cumplen la función de sujetar mecánicamente
a los conductores que forman parte de la línea, manteniéndolos aislados de tierra y
de otros conductores [12].
Rigidez mecánica: la rigidez mecánica es la capacidad de un material o estructura
para resistir la deformación cuando se aplica una fuerza externa [13].
Rigidez dieléctrica: se entiende por rigidez dieléctrica o rigidez electrostática el
valor límite de la intensidad del campo eléctrico en el cual un material pierde su
propiedad aislante y pasa a ser conductor [14].
Mordazas: son herrajes de unión entre el conductor y las cadenas de aisladores,
se distinguen 2 tipos principales (suspensión ya amarre).
Cargabilidad de una línea: Es la cantidad máxima de potencia eléctrica que se
puede trasladar en una línea de transmisión, regularmente manteniendo las
condiciones normales de operación o al menos en condiciones aceptables de la
misma línea [15].
Capacidad térmica del conductor: limitación térmica, resultando del calor
producido por las pérdidas óhmicas, siendo dependiente de la máxima temperatura
del conductor.
Zona Urbana: son las localidades o áreas con 5 000 habitantes o más, o bien, las
cabeceras municipales independientemente del número de habitantes [16].
Zona Rural: son las localidades o áreas con menos de 5 000 habitantes [16].
Potencia real: representa la potencia eléctrica real utilizada por un dispositivo o
sistema de una red eléctrica para realizar un trabajo. En otras palabras, es la
potencia útil que acciona un motor, genera luz o hace funcionar aparatos eléctricos,
por ejemplo.
Potencia reactiva: representa la potencia eléctrica que fluye de un lado a otro entre
los conductores de fase y el conductor neutro de una red trifásica, pero que no
realiza ningún trabajo mecánico. Es generada principalmente por cargas inductivas
y capacitivas y es necesaria para el funcionamiento de estas cargas, pero no
contribuye a la realización de tareas.

96
Potencia aparente: es la potencia eléctrica total disponible en una red eléctrica,
independientemente de si se utiliza realmente para realizar tareas. La potencia
aparente se compone de potencia activa y potencia reactiva.
Factor de potencia: es la relación entre la potencia activa y la potencia aparente.
Su valor puede ir de 0 a 1, donde 1 es ideal e indica que toda la energía consumida
se transformó en trabajo [17].
Caída de tensión: la caída de tensión es la reducción de la tensión en un circuito
eléctrico entre la fuente y la carga [18].
Margen de estabilidad: los márgenes de estabilidad expresan cuánto pueden
cambiar los parámetros antes que el sistema asintóticamente estable se torne
críticamente estable [19].
Impedancia característica: se denomina impedancia característica de una línea
de transmisión a la relación existente entre la diferencia de potencial aplicada y la
corriente absorbida por la línea en el caso hipotético de que esta tenga una longitud
infinita [20].
Potencia eléctrica: es la proporción por unidad de tiempo, o ritmo, con la cual la
energía eléctrica es transferida por un circuito eléctrico, es decir, la cantidad de
energía eléctrica entregada o absorbida por un elemento en un momento
determinado [21].
Transferencia de potencia: la transferencia de potencia eléctrica es la cantidad de
energía eléctrica que se transfiere de una fuente generadora a un elemento
consumidor en un tiempo determinado [21].
Hilos de Guarda: son conductores colocados en la parte más alta de las torres de
transmisión, que se encuentran conectados a tierra para ofrecer un blindaje a los
conductores de fase, protegiéndolos contra las descargas atmosféricas [22].
Métodos semiempíricos: los métodos denominados semiempíricos son
caracterizados como aquellos que involucran supuestos, aproximaciones o
generalizaciones diseñadas para simplificar el cálculo o producir un resultado con
suficiente adecuación empírica [23].
Métodos empíricos: el procedimiento empírico es un modelo de investigación que
pretende obtener conocimiento a partir de la observación de la realidad. Por ende,
está basada en la experiencia [24].
CIGRE: consejo Internacional de Grandes Sistemas Eléctricos (CIGRE) [25].
CISPR: Comité Especial Internacional de Perturbaciones Radioeléctricas (CEIPR),
en idioma francés: Comité International Spécial des Perturbations Radioélectriques
(CISPR) [26].

97
IEE: instituto de Ingenieros Eléctricos y Electrónicos, una organización mundial que
tiene como uno de sus objetivos la definición de estándares que faciliten el
intercambio de tecnologías en diferentes áreas [27].
Caso de estudio: un caso de estudio es una metodología de investigación que
consiste en analizar un caso específico para comprenderlo en profundidad [28].
Método del elemento finito: es un método numérico general para la aproximación
de soluciones de ecuaciones diferenciales parciales muy complejas utilizado en
diversos problemas de ingeniería física [29].
Desplazamiento angular: el desplazamiento angular de un cuerpo es el ángulo con
el que se indica el giro de un elemento alrededor de un eje de rotación determinado,
respecto a una orientación inicial y en un sentido específico [30].
Dominio de un modelo: es un sistema de abstracciones que describe aspectos de
un dominio, es decir, una esfera de conocimiento, influencia o actividad [31].
Mallado: es un paso fundamental en la construcción de un modelo y se refiere a la
generación de una malla del modelo. La malla es básica para la discretización del
Método de los Elementos Finitos [32].
Gradiente eléctrico: el gradiente eléctrico es la diferencia de carga eléctrica que
se produce en dos puntos distintos. Se genera a partir de la separación de iones y
moléculas cargadas [33].
Banda de Radiofrecuencia: es un intervalo de frecuencias del espectro
radioeléctrico que se utiliza para comunicaciones por radio

98
 Índice de figuras
Figura 1. Estructura autosoportada [17]. ......................................................................................... 10
Figura 2. Estructura autosoportada tipo tubular [17]. .................................................................... 11
Figura 3. Conductor AAC [18]............................................................................................................ 13
Figura 4. Conductor AAAC [18]. ........................................................................................................ 13
Figura 5. Conductor ACSR [18]. ........................................................................................................ 13
Figura 6. Conductor ACAR [18]. ........................................................................................................ 14
Figura 7. Haz de seis conductores [3].............................................................................................. 15
Figura 8. Aislador de Vidrio [2]. ......................................................................................................... 16
Figura 9. Aislador de Porcelana [2]. ................................................................................................. 17
Figura 10. Línea de transmisión HSIL convencional [1]. ............................................................... 21
Figura 11. Línea de transmisión HSIL no convencional [1]........................................................... 22
Figura 12. Modeló para análisis de cargabilidad en LT [12]. ........................................................ 24
Figura 13. Parámetros de una línea de transmisión []................................................................... 25
Figura 14. Campo magnético y campo eléctrico asociados a un alinea de dos conductores [].
................................................................................................................................................................ 27
Figura 15. Superficie gaussiana con cargas eléctricas en su interior []...................................... 29
Figura 16. Integración del derecho de Vía [20]............................................................................... 31
Figura 17. Conductor con efecto corona [3]. ................................................................................... 34
Figura 18. Diagrama para el cálculo de los coeficientes de Maxwell [10]. ................................. 39
Figura 19. Haz de 4 subconductores [11]. ....................................................................................... 40
Figura 20. Esquema del campo eléctrico en la superficie de los conductores [9]. ................... 42
Figura 21. Área de la sección transversal de las fases. ................................................................ 53
Figura 22. Selección del Módulo de Física. .................................................................................... 54
Figura 23. Definición del estudio....................................................................................................... 54
Figura 24. (a) Barra de herramientas (b) Constructor de modeló (c) Barra de ajustes (d) Barra
de materiales (e) Gráficos (f) Barra de mensajes, estado de cálculo, registro y tabla. ............ 55
Figura 25. Distribución y dimensiones de la línea de transmisión convencional con 4
subconductores por haz, operando a 500 kV (Caso 1). ................................................................ 56
Figura 26. Círculo con su conversión de grados a radianes. ....................................................... 56
Figura 27. Posición de los cuatro subconductores por fase de la línea del caso 1. ................. 58
Figura 28. (a) Constructor de modelo “Parámetros 1”, (b) Ajustes de Parámetros del caso 1.
................................................................................................................................................................ 60
Figura 29. (a)Constructor de modelo, (b) Ajustes del círculo para el domo de aire. ................ 63
Figura 30. (a) Domo de Aire, (b) Domo para el límite de elementos infinitos, (c)
Subconductores de las Fases A, B y C............................................................................................ 63
Figura 31. Seleccionar casilla............................................................................................................ 64
Figura 32. Añadido del domo de aire (c2) y del domo para el límite de elementos infinitos
(c1). ........................................................................................................................................................ 64
Figura 33. Sustracción para que los conductores de las fases A, B y C queden huecos. ....... 65
Figura 34. (a) Funciones de definiciones, (b) Ajustes del dominio de elemento finito. ............ 65
Figura 35. Selección del material “Aire”. .......................................................................................... 66
Figura 36. Ajuste de las condiciones de contorno para la “Tierra”. ............................................. 66
Figura 37. Asignación del potencial eléctrico a cada subconductor de acuerdo a su fase (A, B
o C). ....................................................................................................................................................... 67

1
Figura 38. Mallado ............................................................................................................................... 68
Figura 39. (a) Mallado en la superficie de los subconductores de la fase A, (b) Mallado en la
superficie del conductor 1 de la fase A ............................................................................................ 68
Figura 40. Ejecución de la simulación.............................................................................................. 69
Figura 41. Análisis de resultados. ..................................................................................................... 69
Figura 42. Distribución del campo eléctrico en la superficie de los subconductores del caso 1.
................................................................................................................................................................ 70
Figura 43. Distribución del campo eléctrico en la superficie de los subconductores de la fase
A del caso 1. ......................................................................................................................................... 70
Figura 44. (a) Selección de la sonda de medición, (b) Ajuste de la sonda de contorno con
respecto a cada fase........................................................................................................................... 71
Figura 45. Resultado del Campo eléctrico para la Fase A. .......................................................... 72
Figura 46. Distribución y dimensiones de la línea de transmisión de alta cargabilidad (HSIL)
convencional con 6 subconductores por haz, operando a 500 kV (Caso 2).............................. 72
Figura 47. Distribución del campo eléctrico en la superficie de los subconductores del caso 2.
................................................................................................................................................................ 74
Figura 48. Distribución del campo eléctrico en la superficie de los subconductores de la fase
A del caso 2. ......................................................................................................................................... 74
Figura 49. Distribución y dimensiones de la línea de transmisión de alta cargabilidad (HSIL)
convencional con 6 subconductores por haz, operando a 500 kV (Caso 2).............................. 75
Figura 50. Distribución del campo eléctrico en la superficie de los subconductores del caso 3.
................................................................................................................................................................ 76
Figura 51. Distribución del campo eléctrico en la superficie de los subconductores de la fase
A del caso 3. ......................................................................................................................................... 76
Figura 52. Distribución y dimensiones de la línea de transmisión de alta cargabilidad (HSIL)
no convencional con 6 subconductores por haz, operando a 500 kV (Caso 4) ........................ 77
Figura 53. Distribución del campo eléctrico en la superficie de los subconductores del caso 4.
................................................................................................................................................................ 78
Figura 54. Distribución del campo eléctrico en la superficie de los subconductores de la fase
A del caso 4. ......................................................................................................................................... 78
Figura 55. Perfiles de radio interferencia de los 4 casos de estudio en Buen ambiente
(CIGREE).............................................................................................................................................. 79
Figura 56. Perfiles de radio interferencia de los 4 casos de estudio en Lluvia moderada
(CIGREE).............................................................................................................................................. 80
Figura 57. Perfiles de radio interferencia de los 4 casos de estudio en Lluvia densa
(CIGREE).............................................................................................................................................. 80
Figura 58. Perfiles de radio interferencia de los 4 casos de estudio en Buen ambiente (BPA).
................................................................................................................................................................ 81
Figura 59. Perfiles de radio interferencia de los 4 casos de estudio en Lluvia densa (BPA). . 82
Figura 60. Perfiles de radio interferencia de los 4 casos de estudio en Buen ambiente
(CREIPI)................................................................................................................................................ 83
Figura 61. Perfiles de radio interferencia de los 4 casos de estudio en Lluvia moderada
(CREIPI)................................................................................................................................................ 84
Figura 62. Perfiles de radio interferencia de los 4 casos de estudio en Lluvia densa (CREIPI).
................................................................................................................................................................ 84

2
 Índice de tablas
Tabla 1. Características de los aisladores de vidrio y porcelana []. ............................................. 17
Tabla 2. Líneas de transmisión HSIL existentes actualmente [1]. ............................................... 22
Tabla 3. Resistividad del material de los conductores [11]. .......................................................... 26
Tabla 4. Factores que afectan el efecto corona []. ......................................................................... 35
Tabla 5. Límites máximos de condiciones climáticas favorables de la Asociación Canadiense
de Normas para campos de interferencia (0.15-30 MHz) [19]. .................................................... 49
Tabla 6. Datos de la línea de transmisión convencional con 4 subconductores por haz,
operando a 500 Kv. ............................................................................................................................. 55
Tabla 7. Parámetros del caso 1 ......................................................................................................... 60
Tabla 8. Características de los círculos que representan los subconductores de las Fases A,
B y C. ..................................................................................................................................................... 62
Tabla 9. Resultados del campo eléctrico máximo para el caso 1. ............................................... 72
Tabla 10. Parámetros del caso 2. ..................................................................................................... 73
Tabla 11. Resultados del campo eléctrico máximo para el caso 2. ............................................. 74
Tabla 12. Parámetros del caso 3. ..................................................................................................... 75
Tabla 13. Resultados del campo eléctrico máximo para el caso 3 .............................................. 76
Tabla 14. Parámetros del caso 4. ..................................................................................................... 77
Tabla 15. Resultados del campo eléctrico máximo para el caso 4. ............................................. 78
Tabla 16. Valores obtenidos de RI con la ecuación CIGREE para los 4 casos de estudio y
para cada tipo de ambiente analizado. ............................................................................................ 79
Tabla 17. Valores obtenidos de RI con la ecuación BPA para los 4 casos de estudio y para
cada tipo de ambiente analizado. ..................................................................................................... 81
Tabla 18. Valores obtenidos de RI con la ecuación CREIPI para los 4 casos de estudio y para
cada tipo de ambiente analizado. ..................................................................................................... 83
Tabla 19. Precios de las licencias de MATLAB® [3]. ..................................................................... 87
Tabla 20. Precios de las licencias de COMSOL Multiphysics® [5]. ............................................. 87
Tabla 21. Estimado del costo final del proyecto. ............................................................................ 88
Tabla 22. Costos unitarios del proyecto. .......................................................................................... 89

3
Apéndices

4
Apéndice A: Cálculo del campo eléctrico por método de Mark y Mengele para los
casos 1 y 2.
Para el cálculo del campo eléctrico se considera el método de Mark y Mengele, que a
su vez se relaciona con los valores de carga, calculados con el método de los
potenciales de Maxwell.
1) Caso 1.
De la ecuación N, 𝑄 = 𝑃−1 ∗ 𝑉, se calculan primeramente los valores de potenciales
de Maxwell (𝑃), siendo la estructura de una matriz cuadrada de acuerdo al número de
haz de conductores, es decir una matriz 𝟑 ∗ 𝟑, de acuerdo a las ecuaciones N y N:
1 𝐻𝑖
𝑃𝑖𝑖 = ln ( )
2𝜋𝜀0 𝑅𝑖
1 𝐷𝑖𝑗
𝑃𝑖𝑗 = ln ( )
2𝜋𝜀0 𝑑𝑖𝑗

Siendo 𝑃𝑖𝑖 la diagonal de la matriz de potenciales propios de cada haz de conductores

y 𝑃𝑖𝑗 los valores de potenciales mutuos que rellenan la matriz.

Para este caso, los valores de alturas 𝐻 son de la misma magnitud para cada calculo,
la 𝑅𝑖 es el equivalente al valor obtenido de radio medio geométrico o en todo caso,
radio equivalente, siendo de acuerdo a las ecuaciones N y n:
𝑛2
𝑅𝑀𝐺 = √𝑟𝑅𝑀𝐺 (𝐷12 ∗ 𝐷12 … . 𝐷𝑚𝑛 )

𝑛
𝑅𝑒 = √𝑛𝑟𝑅ℎ𝑎𝑧 𝑛−1

Donde para determinar las distancias entre cada sub conductor de cada haz, se
emplea la ecuación n de distancias entre dos puntos:

𝐷 = √(𝑥2 − 𝑥1 )2 + (𝑦2 − 𝑦1 )2
De acuerdo a la configuración del haz:

• Mediante el método de RMG:

𝐷12 = √(−12.525 − (−12.075))2 + (28.225 − 28.225)2 = 0.45 𝑚

5
 𝐷13 = √(−12.525 − (−12.075))2 + (27.775 − 28.225)2 = 0.646 𝑚

𝐷14 = √(−12.075 − (−12.075))2 + (27.775 − 28.225)2 = 0.45 𝑚

En este caso, al ser una configuración simétrica se comprende que los valores
obtenidos son los mismos para cada conductor tomado como referencia, además, el
valor de 𝑟𝑅𝑀𝐺 es el radio geométrico medio del conductor, tomado del anexo N, siendo
para un calibre de 715 500 CM:
42
𝑅𝑀𝐺 = √(0.01082)4 (0.45 ∗ 0.646 ∗ 0.45)4 = 0.1939 𝑚 = 19.39 𝑐𝑚

• Mediante el método de 𝑹𝒆 :
4
𝑅𝑒 = √4(1.3147)(31.82)4−1 = 20.288 𝑐𝑚
Con lo anterior, se calculan los coeficientes de potenciales de Maxwell propios, siendo
3 coeficientes, uno para cada haz:
1 28
𝑃11 = ln ( ) = 8.9383𝑒 10
2𝜋𝜀0 0.1939
1 28
𝑃22 = ln ( ) = 8.9383𝑒 10
2𝜋𝜀0 0.1939
1 28
𝑃33 = ln ( ) = 8.9383𝑒 10
2𝜋𝜀0 0.1939
Ahora, para los coeficientes de potenciales de Maxwell mutuas, se deben de calcular
de acuerdo al método de imágenes, siendo para el caso el siguiente diagrama:

6
Donde para las distancias del haz de conductores real (𝐴, 𝐵, 𝐶) se calcula con respecto
a los haces de conductores imaginarios (𝐴′ , 𝐵 ′ , 𝐶 ′ ) aplicando teorema de Pitágoras:

𝐷𝐴𝐵′ = √(12.3)2 + (56)2 = 57.3348 𝑚

𝐷𝐴𝑐 ′ = √(24.6)2 + (56)2 = 61.165 𝑚

𝑑𝐴𝐵 = 12.3 𝑚
𝑑𝐴𝐶 = 24.6 𝑚
1 57.3348 1 61.165
𝑃12 = ln ( ) = 2.7669𝑒 10 𝑃13 = ln ( ) = 1.6372𝑒 10
2𝜋𝜀0 12.3 2𝜋𝜀0 24.6
1 57.3348 1 57.3348
𝑃21 = ln ( ) = 2.7669𝑒 10 𝑃23 = ln ( ) = 2.7669𝑒 10
2𝜋𝜀0 12.3 2𝜋𝜀0 12.3
1 61.165 1 57.3348
𝑃31 = ln ( ) = 1.6372𝑒 10 𝑃32 = ln ( ) = 2.7669𝑒 10
2𝜋𝜀0 24.6 2𝜋𝜀0 12.3

Siendo la matriz de potenciales de Maxwell [𝑃] la siguiente:

8.9383𝑒 10 2.7669𝑒 10 1.6372𝑒 10
[ 2.7669𝑒 10 8.9383𝑒 10 2.7669𝑒 10 ]
1.6372𝑒 10 2.7669𝑒 10 8.9383𝑒 10

Invirtiendo la matriz [𝑃]:

1.2490𝑒 10 −3.4928𝑒 10 −1.2065𝑒 10
[ −3.4928𝑒 10 1.2490𝑒 10 −3.4928𝑒 10 ]
−1.2065𝑒 10 −3.4928𝑒 10 1.2490𝑒 10
Para el vector de tensiones [𝑉]
2𝜋 2𝜋
500𝑒 3 (cos 3 + 𝑖 sin 3 )
√3
3 0 0
500𝑒 (cos 3 + 𝑖 sin 3) −1.4434𝑒 5 + 2.5𝑒 5
=[ 2.8868𝑒 5 ]
√3 5 5
−1.4434𝑒 − 2.5𝑒
4𝜋 4𝜋
500𝑒 3 (cos 3 + 𝑖 sin 3 )
[ √3 ]
Calculando las cargas:
𝑄 = [𝑃]−1 ∗ [𝑉]

7
 1.2490𝑒 −10 −3.4928𝑒 −10 −1.2065𝑒 −10 −1.4434𝑒 5 + 2.5𝑒 5
𝑄 = [ −3.4928𝑒 −10 1.2490𝑒 −10 −3.4928𝑒 −10 ] ∗ [ 2.8868𝑒 5 ]
−1.2065𝑒 −10 −3.4928𝑒 −10 1.2490𝑒 −10 5
−1.4434𝑒 − 2.5𝑒 5

−2.1848𝑒 − 06 + 2.9521𝑒 − 6𝑖
𝑄=[ 4.0537𝑒 − 6 ]
−2.1848𝑒 − 6 − 2.9521𝑒 − 6𝑖

Con el valor de las cargas se procede a calcular el campo eléctrico medio:

𝑄1 −2.1848𝑒 − 06 + 2.9521𝑒 − 6𝑖
𝐸𝑚𝑒𝑑 𝐴 = = = −7.467𝑒 3 + 10.09𝑒 3
2𝜋 ∗ 𝜀0 ∗ 𝑛 ∗ 𝑟 2𝜋 ∗ 𝜀0 ∗ 4 ∗ 1.3147
𝑄1 4.0537𝑒 − 6
𝐸𝑚𝑒𝑑 𝐵 = = = 13.855𝑒 3
2𝜋 ∗ 𝜀0 ∗ 𝑛 ∗ 𝑟 2𝜋 ∗ 𝜀0 ∗ 4 ∗ 1.3147
𝑄1 −2.1848𝑒 − 06 − 2.9521𝑒 − 6𝑖
𝐸𝑚𝑒𝑑 𝐶 = = = −7.467𝑒 3 − 10.09𝑒 3
2𝜋 ∗ 𝜀0 ∗ 𝑛 ∗ 𝑟 2𝜋 ∗ 𝜀0 ∗ 4 ∗ 1.3147
Determinando sus valores absolutos:
𝑘𝑉
𝐸𝑚𝑒𝑑 𝐴 = 12.553
𝑐𝑚
𝑘𝑉
𝐸𝑚𝑒𝑑 𝐵 = 13.855
𝑐𝑚
𝑘𝑉
𝐸𝑚𝑒𝑑 𝐶 = 12.553
𝑐𝑚
Con los valores de campo eléctrico medio, se calcula el campo eléctrico máximo en
cada haz de conductores:
𝑟 1.3147 𝑘𝑉
𝐸max 𝐴 = 𝐸𝑚𝑒𝑑 𝐴 ∗ 1 + ( ) (𝑛 − 1) = 12.553𝑒 3 ∗ 1 + ( ) (4 − 1) = 14.108
𝑅ℎ𝑎𝑧 31.82 𝑐𝑚
𝑟 1.3147 𝑘𝑉
𝐸max 𝐵 = 𝐸𝑚𝑒𝑑 𝐵 ∗ 1 + ( ) (𝑛 − 1) = 13.855𝑒 3 ∗ 1 + ( ) (4 − 1) = 15.572
𝑅ℎ𝑎𝑧 31.82 𝑐𝑚
𝑟 1.3147 𝑘𝑉
𝐸max 𝐶 = 𝐸𝑚𝑒𝑑 𝐶 ∗ 1 + ( ) (𝑛 − 1) = 12.553𝑒 3 ∗ 1 + ( ) (4 − 1) = 14.108
𝑅ℎ𝑎𝑧 31.82 𝑐𝑚

8
 kV kV
Emed [ ] Emax [ ]
Fase cm cm

A 12.553 14.108

B 13.855 15.572

C 12.553 14.108

2) Caso 2.
De la ecuación N, 𝑄 = 𝑃−1 ∗ 𝑉, se calculan primeramente los valores de potenciales
de Maxwell (𝑃), siendo la estructura de una matriz cuadrada de acuerdo al número de
haz de conductores, es decir una matriz 3 ∗ 3, de acuerdo a las ecuaciones N y N:
1 𝐻𝑖
𝑃𝑖𝑖 = ln ( )
2𝜋𝜀0 𝑅𝑖
1 𝐷𝑖𝑗
𝑃𝑖𝑗 = ln ( )
2𝜋𝜀0 𝑑𝑖𝑗

Siendo 𝑃𝑖𝑖 la diagonal de la matriz de potenciales propios de cada haz de conductores

y 𝑃𝑖𝑗 los valores de potenciales mutuos que rellenan la matriz.

Para este caso, los valores de alturas 𝐻 son de la misma magnitud para cada calculo,
la 𝑅𝑖 es el equivalente al valor obtenido de radio medio geométrico o en todo caso,
radio equivalente, siendo de acuerdo a las ecuaciones N y n:
𝑛2
𝑅𝑀𝐺 = √𝑟𝑅𝑀𝐺 (𝐷12 ∗ 𝐷12 … . 𝐷𝑚𝑛 )

𝑛
𝑅𝑒 = √𝑛𝑟𝑅ℎ𝑎𝑧 𝑛−1

Donde para determinar las distancias entre cada sub conductor de cada haz, se
emplea la ecuación n de distancias entre dos puntos:

𝐷 = √(𝑥2 − 𝑥1 )2 + (𝑦2 − 𝑦1 )2
De acuerdo a la configuración del haz:

9
 • Mediante el método de RMG:

𝐷12 = √(−3.725 − (−2.975))2 + (32.65 − 32.65)2 = 0.75 𝑚

𝐷13 = √(−4.1 − (−2.975))2 + (33.3 − 32.65)2 = 1.29 𝑚

𝐷14 = √(−3.725 − (−2.975))2 + (33.95 − 32.65)2 = 1.5 𝑚

𝐷15 = √(−2.975 − (−2.975))2 + (33.95 − 32.65)2 = 1.3 𝑚

𝐷16 = √(−2.6 − (−2.975))2 + (33.3 − 32.65)2 = 0.75 𝑚

En este caso, al ser una configuración simétrica se comprende que los valores
obtenidos son los mismos para cada conductor tomado como referencia, además, el
valor de 𝑟𝑅𝑀𝐺 es el radio geométrico medio del conductor, tomado del anexo N, siendo
para un calibre de 477 000 CM:

62
𝑅𝑀𝐺 = √(0.00883)6 (0.75 ∗ 1.29 ∗ 1.5 ∗ 1.3 ∗ 0.75)6 = 0.4817 𝑚 = 48.17 𝑐𝑚

• Mediante el método de 𝑹𝒆 :
6
𝑅𝑒 = √6(1.088)(75)6−1 = 49.92 𝑐𝑚
Con lo anterior, se calculan los coeficientes de potenciales de Maxwell propios, siendo
3 coeficientes, uno para cada haz:

1 33.3
𝑃11 = ln ( ) = 7.6142𝑒 10
2𝜋𝜀0 0.4817
1 27.5
𝑃22 = ln ( ) = 7.2702𝑒 10
2𝜋𝜀0 0.4817
1 33.3
𝑃33 = ln ( ) = 7.6142𝑒 10
2𝜋𝜀0 0.4817

Ahora, para los coeficientes de potenciales de Maxwell mutuas, se deben de calcular

de acuerdo al método de imágenes, siendo para el caso el siguiente diagrama:

10
Donde para las distancias del haz de conductores real (𝐴, 𝐵, 𝐶) se calcula con respecto
a los haces de conductores imaginarios (𝐴′ , 𝐵 ′ , 𝐶 ′ ) aplicando teorema de Pitágoras:

𝐷𝐴𝐵′ = √(3.35)2 + (60.8)2 = 60.8922 𝑚

𝐷𝐴𝑐 ′ = √(6.7)2 + (66.6)2 = 66.9361 𝑚

𝑑𝐴𝐵 = 6.7 𝑚
𝑑𝐴𝐶 = 6.7 𝑚
1 60.8922 1 66.9361
𝑃12 = ln ( ) = 3.967𝑒 10 𝑃13 = ln ( ) = 4.1371𝑒 10
2𝜋𝜀0 6.7 2𝜋𝜀0 6.7
1 60.8922 1 60.8922
𝑃21 = ln ( ) = 3.967𝑒 10 𝑃23 = ln ( ) = 3.967𝑒 10
2𝜋𝜀0 6.7 2𝜋𝜀0 6.7
1 66.9361 1 60.8922
𝑃31 = ln ( ) = 4.1371𝑒 10 𝑃32 = ln ( ) = 3.967𝑒 10
2𝜋𝜀0 6.7 2𝜋𝜀0 6.7

Siendo la matriz de potenciales de Maxwell [𝑃] la siguiente:

7.6142𝑒 10 ′ 3.967𝑒 10 4.1371𝑒 10

[ 3.967𝑒 10 2.1778𝑒 − 11 3.967𝑒 10 ]
4.1371𝑒 10 3.967𝑒 10 7.6142𝑒 10

11
Invirtiendo la matriz [𝑃]:
2.1116𝑒 − 11 −7.3517𝑒 − 12 −7.6432𝑒 − 12
[ −7.3517𝑒 − 12 2.1778𝑒 − 11 −7.3517𝑒 − 12 ]
−7.6432𝑒 − 12 −7.3517𝑒 − 12 2.1116𝑒 − 11

Para el vector de tensiones [𝑉]

2𝜋 2𝜋
500𝑒 3 (cos 3 + 𝑖 sin 3 )
√3
0 0
500𝑒 3 (cos 3 + 𝑖 sin 3) −1.4434𝑒 5 + 2.5𝑒 5
=[ 2.8868𝑒 5 ]
√3 5 5
−1.4434𝑒 − 2.5𝑒
3 4𝜋 4𝜋
500𝑒 (cos 3 + 𝑖 sin 3 )
[ √3 ]
Calculando las cargas:
𝑄 = [𝑃]−1 ∗ [𝑉]

2.1116𝑒 − 11 −7.3517𝑒 − 12 −7.6432𝑒 − 12 −1.4434𝑒 5 + 2.5𝑒 5

𝑄 = [ −7.3517𝑒 − 12 2.1778𝑒 − 11 −7.3517𝑒 − 12 ] ∗ [ 2.8868𝑒 5 ]
−7.6432𝑒 − 12 −7.3517𝑒 − 12 2.1116𝑒 − 11 5 5
−1.4434𝑒 − 2.5𝑒

−4.0670𝑒 − 6 + 7.1899𝑒 − 6𝑖
𝑄 = [8.4091𝑒 − 6 + 1.0588𝑒 − 21𝑖 ]
−4.0670𝑒 − 6 − 7.1899𝑒 − 6𝑖

Con el valor de las cargas se procede a calcular el campo eléctrico medio:

𝑄1 −4.0670𝑒 − 6 + 7.1899𝑒 − 6𝑖
𝐸𝑚𝑒𝑑 𝐴 = = = −11.198𝑒 3 + 19.797𝑒 3
2𝜋 ∗ 𝜀0 ∗ 𝑛 ∗ 𝑟 2𝜋 ∗ 𝜀0 ∗ 6 ∗ 1.088
𝑄2 8.4091𝑒 − 6 + 1.0588𝑒 − 21𝑖
𝐸𝑚𝑒𝑑 𝐵 = = = 23.154𝑒 3 + 2.9154𝑒 −12
2𝜋 ∗ 𝜀0 ∗ 𝑛 ∗ 𝑟 2𝜋 ∗ 𝜀0 ∗ 6 ∗ 1.088
𝑄3 −4.0670𝑒 − 6 − 7.1899𝑒 − 6𝑖
𝐸𝑚𝑒𝑑 𝐶 = = = −11.198𝑒 3 − 19.797𝑒 3
2𝜋 ∗ 𝜀0 ∗ 𝑛 ∗ 𝑟 2𝜋 ∗ 𝜀0 ∗ 6 ∗ 1.088

Determinando sus valores absolutos:

12
 𝑘𝑉
𝐸𝑚𝑒𝑑 𝐴 = 22.745
𝑐𝑚
𝑘𝑉
𝐸𝑚𝑒𝑑 𝐵 = 23.154
𝑐𝑚
𝑘𝑉
𝐸𝑚𝑒𝑑 𝐶 = 22.745
𝑐𝑚
Con los valores de campo eléctrico medio, se calcula el campo eléctrico máximo en
cada haz de conductores:
𝑟 1.088 𝑘𝑉
𝐸max 𝐴 = 𝐸𝑚𝑒𝑑 𝐴 ∗ 1 + ( ) (𝑛 − 1) = 22.745𝑒 3 ∗ 1 + ( ) (6 − 1) = 24.394
𝑅ℎ𝑎𝑧 75 𝑐𝑚
𝑟 1.088 𝑘𝑉
𝐸max 𝐵 = 𝐸𝑚𝑒𝑑 𝐵 ∗ 1 + ( ) (𝑛 − 1) = 23.154𝑒 3 ∗ 1 + ( ) (6 − 1) = 24.833
𝑅ℎ𝑎𝑧 75 𝑐𝑚
𝑟 1.088 𝑘𝑉
𝐸max 𝐶 = 𝐸𝑚𝑒𝑑 𝐶 ∗ 1 + ( ) (𝑛 − 1) = 22.745𝑒 3 ∗ 1 + ( ) (6 − 1) = 24.394
𝑅ℎ𝑎𝑧 75 𝑐𝑚
kV kV
Emed [ ] Emax [ ]
Fase cm cm

A 22.745 24.394

B 23.154 24.833

C 22.745 24.394

13
Apéndice B: Cálculo de Radio Interferencia (RI) para el caso 1.
A. CIGREE
Para el cálculo de RI, se emplea la formula siguiente:
𝐷 𝑞
𝑪𝑰𝑮𝑹𝑬𝑬: 𝑅𝐼𝐿𝐿𝐷 = −10 + 3.5 𝑔 + 6 𝑑 − 33 log10 ( )+
20 300
Siendo la ecuación anterior para lluvias densas, por lo que para regular los resultados
en condiciones de ambiente adecuado y lluvias moderadas se considera lo siguiente:
𝐿𝑙𝑢𝑣𝑖𝑎 𝑚𝑜𝑑𝑒𝑟𝑎𝑑𝑎: 𝑅𝐼𝐿𝐿𝑀 = 𝑅𝐼𝐿𝐿𝐷 − 7
𝐵𝑢𝑒𝑛 𝑎𝑚𝑏𝑖𝑒𝑛𝑡𝑒: 𝑅𝐼𝐵𝐴 = 𝑅𝐼𝐿𝐿𝐷 − 24
Siendo:

𝑔𝐴 = 14.13 𝑘𝑉 𝑟𝑚𝑠 ⁄𝑐𝑚 𝑞 = 2000 𝑚𝑠𝑛𝑚

𝑔𝐵 = 15.36 𝑘𝑉 𝑟𝑚𝑠 ⁄𝑐𝑚 𝐷 =?

𝑔𝐶 = 14.13 𝑘𝑉 𝑟𝑚𝑠 ⁄𝑐𝑚

𝑑 = 2 ∗ 1.341 𝑐𝑚 = 2.682 𝑐𝑚 (𝑑𝑖𝑎𝑚𝑒𝑡𝑟𝑜 𝑑𝑒𝑙 𝑐𝑜𝑛𝑑𝑢𝑐𝑡𝑜𝑟)
Calculando las distancias radiales:

𝐷 = √(𝑦𝑚 − 𝑦𝑖 )2 + (𝑥𝑚 − 𝑥𝑖 )2
𝑦𝑚 = 𝐴𝑙𝑡𝑢𝑟𝑎 𝑑𝑒 𝑙𝑎𝑠 𝑓𝑎𝑠𝑒𝑠
𝑥𝑖 = 𝐷𝑖𝑠𝑡𝑎𝑛𝑐𝑖𝑎𝑠 𝑙𝑎𝑡𝑒𝑟𝑎𝑙𝑒𝑠 𝑑𝑒 𝑙𝑎𝑠 𝑓𝑎𝑠𝑒𝑠
Siendo:
𝑦𝐴 = 28 𝑚 𝑥𝑚 = −50: 1: 50 𝑚
𝑦𝐵 = 28 𝑚 𝑥𝐴 = −12.3 𝑚
𝑦𝐵 = 28 𝑚 𝑥𝐵 = 0 𝑚
𝑦𝑖 = 1 𝑚 𝑥𝐴 = 12.3 𝑚

Para cada fase y con el primero valor del vector 𝑥𝑚 :

𝐷𝐴 = √(28 − 1)2 + (0 − (−12.3))2 = 29.6696 𝑚

𝐷𝐵 = √(28 − 1)2 + (0 − 0)2 = 27 𝑚

𝐷𝐶 = √(28 − 1)2 + (0 − (12.3))2 = 29.6696 𝑚

14
Calculando la 𝑅𝐼 por fase:
29.6696 2000
𝑅𝐼𝐿𝐿𝐷−𝐹𝐴𝑆𝐸 𝐴 = −10 + 3.5(14.13) + 6(2.682) − 33 log10 ( )+
20 300
= 𝟓𝟔. 𝟓𝟔𝟏𝟑 𝒅𝑩
27 2000
𝑅𝐼𝐿𝐿𝐷−𝐹𝐴𝑆𝐸 𝐵 = −10 + 3.5(15.36) + 6(2.682) − 33 log10 ( ) +
20 300
= 𝟔𝟐. 𝟐𝟏𝟕𝟔 𝒅𝑩
29.6696 2000
𝑅𝐼𝐿𝐿𝐷−𝐹𝐴𝑆𝐸 𝐶 = −10 + 3.5(14.13) + 6(2.682) − 33 log10 ( )+
20 300
= 𝟓𝟔. 𝟓𝟔𝟏𝟑 𝒅𝑩

Con los resultados anteriores se realiza el análisis de la magnitud de cada resultado,

ya que se debe de ordenar de manera descendiente para poder obtener la RI
resultante de acuerdo a la formula 𝐶𝐼𝑆𝑃𝑅:
𝑆𝑖 𝑅𝐼1 − 𝑅𝐼2 ≥ 3 𝑑𝐵 𝑜 𝑆𝑖 𝑅𝐼1 − 𝑅𝐼2 ≤ 3 𝑑𝐵
𝑅𝐼1 + 𝑅𝐼2
𝑅𝐼𝑇𝑂𝑇𝐴𝐿 = 𝑅𝐼1 𝑅𝐼𝑇𝑂𝑇𝐴𝐿 = + 1.5
2
Por lo que lo que ordenando los resultados:
𝑅𝐼1 = 62.2176 𝑑𝐵
𝑅𝐼2 = 56.5613 𝑑𝐵
𝑅𝐼3 = 56.5613 𝑑𝐵

Llevando el análisis:
𝑆𝑖 62.2176 − 56.5613 ≥ 3 𝑑𝐵 𝑜 𝑆𝑖 57.9126 − 56.5613 ≤ 3 𝑑𝐵
5.6563 ≥ 3 𝑑𝐵 5.6563 ≤ 3 𝑑𝐵

𝑅𝐼𝑇𝑂𝑇𝐴𝐿 = 𝑅𝐼1 𝑅𝐼𝑇𝑂𝑇𝐴𝐿 = (𝑁𝑂 𝐴𝑃𝐿𝐼𝐶𝐴)

𝑅𝐼𝑇𝑂𝑇𝐴𝐿 = 𝟔𝟐. 𝟐𝟏𝟕𝟔 𝒅𝑩
De lo anterior se concluye:
𝑅𝐼𝐿𝐿𝐷 = 𝟔𝟐. 𝟐𝟏𝟕𝟔 𝒅𝑩
𝑅𝐼𝐿𝐿𝑀 = 52.3572 − 7 = 𝟓𝟓. 𝟐𝟏𝟕𝟔 𝒅𝑩
𝑅𝐼𝐵𝐴 = 52.3572 − 24 = 𝟑𝟖. 𝟐𝟏𝟕𝟔 𝒅𝑩

15
Asimismo, el cálculo se realiza para cada punto del vector 𝑥𝑚 .

B. BPA
Para el cálculo de RI, se emplea la formula siguiente:
𝑔 𝑑
𝑩𝑷𝑨: 𝑅𝐼𝐵𝐴 = 46 + 120 log10 ( ) + 40 log10 ( )+
17.56 3.51
𝑞
10(1 − (log10 (𝑓))2 ) + − 𝐶𝟏 + 𝐶2
300
Siendo la ecuación anterior para ambiente adecuado, por lo que para regular los
resultados en condiciones de lluvias moderadas y lluvias densas se considera lo
siguiente:
𝐿𝑙𝑢𝑣𝑖𝑎 𝑚𝑜𝑑𝑒𝑟𝑎𝑑𝑎: 𝑅𝐼𝐿𝐿𝑀 = 𝑅𝐼𝐵𝐴 + 17
𝐿𝑙𝑢𝑣𝑖𝑎 𝑑𝑒𝑛𝑠𝑎: 𝑅𝐼𝐿𝐿𝐷 = 𝑅𝐼𝐵𝐴 + 24
Siendo:

𝑔𝐴 = 14.13 𝑘𝑉 𝑟𝑚𝑠 ⁄𝑐𝑚 𝑞 = 2000 𝑚𝑠𝑛𝑚

𝑔𝐵 = 15.36 𝑘𝑉 𝑟𝑚𝑠 ⁄𝑐𝑚

𝑔𝐶 = 14.13 𝑘𝑉 𝑟𝑚𝑠 ⁄𝑐𝑚

𝑑 = 2 ∗ 1.341 𝑐𝑚 = 2.682 𝑐𝑚 (𝑑𝑖𝑎𝑚𝑒𝑡𝑟𝑜 𝑑𝑒𝑙 𝑐𝑜𝑛𝑑𝑢𝑐𝑡𝑜𝑟)
𝑓 = 0.5 𝑀ℎ𝑧

𝐶𝟏 = 𝑅𝑒𝑓𝑒𝑟𝑖𝑑𝑜 𝑎 𝑙𝑜𝑠 𝑠𝑖𝑔𝑢𝑖𝑒𝑛𝑡𝑒𝑠 𝑝𝑎𝑟á𝑚𝑒𝑡𝑟𝑜𝑠:

31.1 31.1
𝐷1 = 21.04 𝑚 𝐷𝑊1 = = = 62.2
𝑓 0.5
70.55 70.55
𝐸𝐼𝑁𝐷1 = = = 282.2
𝑓2 0.52
𝑣 300000000 276.16 276.16
𝜆= = = 600 𝑚 𝜌1 = = = 7.6711𝑒 −4
𝑓 500000 𝜆2 6002
2 + 0.3 𝜌1 2 + 0.3 (7.6711𝑒 −4 )
𝑓(𝜌1 ) = = = 0.9997
2 + 𝜌1 + 0.6 𝜌1 2 2 + 7.6711𝑒 −4 + 0.6 (7.6711𝑒 −4 )2

31.1 𝑓(𝜌1 ) 31.1( 0.9997 )

𝐸𝑆𝑈1 = = = 62.1833
𝑓 0.5
16
 𝐶1 = 10 log10 ((62.2)2 + (62.1833)2 + (282.2)2 ) = 49.4137

𝐶2
𝐷2 = ?
Calculando las distancias radiales:

𝐷2 = √(𝑦𝑚 − 𝑦𝑖 )2 + (𝑥𝑚 − 𝑥𝑖 )2
𝑦𝑚 = 𝐴𝑙𝑡𝑢𝑟𝑎 𝑑𝑒 𝑙𝑎𝑠 𝑓𝑎𝑠𝑒𝑠
𝑥𝑖 = 𝐷𝑖𝑠𝑡𝑎𝑛𝑐𝑖𝑎𝑠 𝑙𝑎𝑡𝑒𝑟𝑎𝑙𝑒𝑠 𝑑𝑒 𝑙𝑎𝑠 𝑓𝑎𝑠𝑒𝑠
Siendo:
𝑦𝐴 = ℎ𝑐1 = 28 𝑚 𝑥𝑚 = −50: 1: 50 𝑚
𝑦𝐵 = ℎ𝑐2 = 28 𝑚 𝑥𝐴 = −12.3 𝑚
𝑦𝐵 = ℎ𝑐3 = 28 𝑚 𝑥𝐵 = 0 𝑚
𝑦𝑖 = ℎ𝑎 = 1 𝑚 𝑥𝐴 = 12.3 𝑚
Para cada fase y con el primero valor del vector 𝑥𝑚 :

𝐷2𝐴 = √(28 − 1)2 + (0 − (−12.3))2 = 29.6696 𝑚

𝐷2𝐵 = √(28 − 1)2 + (0 − 0)2 = 27 𝑚

𝐷2𝐶 = √(28 − 1)2 + (0 − (12.3))2 = 29.6696 𝑚

De acuerdo a los valores de distancias radiales calculadas, se calcula 𝐷𝑊2 de acuerdo
a sus condiciones:
12 ℎ𝑎 ℎ𝑐 ℎ𝑐
𝐶𝑜𝑛𝑑𝑖𝑐𝑖ó𝑛 1 → 𝑆𝑖 𝐷2 ≤ ∴ 𝐷𝑊2 =
𝜆 𝐾 𝐷2
12 ℎ𝑎 ℎ𝑐 ℎ𝑐 12 ℎ𝑎 ℎ𝑐
𝐶𝑜𝑛𝑑𝑖𝑐𝑖ó𝑛 2 → 𝑆𝑖 𝐷2 > ∴ 𝐷𝑊2 =
𝜆 𝐾 𝐷2 𝜆 𝐷2
Para la fase 𝐴:
12 ℎ𝑎 ℎ𝑐 12(1)(28)
= = 0.56
𝜆 600
29.6696 ≤ 0.56 (𝑁𝑂 𝐴𝑃𝐿𝐼𝐶𝐴)
29.6696 > 0.56

17
Por lo que de acuerdo a la segunda condición:
Para 𝐷𝑊2 :
2𝜋
𝐾= = 0.0104
600
28 12(1)(28)
𝐷𝑊2 = = 1.7127
0.0104(29.6696) 600 (29.6696)
Para 𝐸𝑆𝑈2 :
52.5 𝐷2 52.5 (29.6696)
𝜌2 = 2
= = 1.0817𝑒 −3
𝛿𝜆 4 (6002 )
2 + 0.3 𝜌2 2 + 0.3 (1.0817𝑒 −3 )
𝑓(𝜌1 ) = = = 0.9996
2 + 𝜌2 + 0.6 𝜌2 2 2 + 1.0817𝑒 −3 + 0.6 (1.0817𝑒 −3 )2
ℎ𝑐 𝑓(𝜌2 ) 28( 0.9996 )
𝐸𝑆𝑈1 = = = 90.7066
𝐾 𝐷2 0.0104(29.6696)
Para 𝐸𝐼𝑁𝐷2 :
ℎ𝑐 28
𝐸𝐼𝑁𝐷2 = = = 294.0816
(𝐾 𝐷2)2 (0.0104 ∗ 29.6696)2
𝐶2 = 10 log10 ((1.7127)2 + (90.7066)2 + (294.0816)2 ) = 49.7641
Para la fase 𝐵:
12 ℎ𝑎 ℎ𝑐 12(1)(28)
= = 0.56
𝜆 600
29.6696 ≤ 0.56 (𝑁𝑂 𝐴𝑃𝐿𝐼𝐶𝐴)
29.6696 > 0.56
Por lo que de acuerdo a la segunda condición:
Para 𝐷𝑊2 :
2𝜋
𝐾= = 0.0104
600
28 12(1)(28)
𝐷𝑊2 = = 2.0681
0.0104(27) 600 (27)
Para 𝐸𝑆𝑈2 :
52.5 𝐷2 52.5 (27)
𝜌2 = 2
= = 9.8437𝑒 −4
𝛿𝜆 4 (6002 )

18
 2 + 0.3 𝜌2 2 + 0.3 (9.8437𝑒 −4 )
𝑓(𝜌1 ) = = = 0.9996
2 + 𝜌2 + 0.6 𝜌2 2 2 + 9.8437𝑒 −4 + 0.6 (9.8437𝑒 −4 )2
ℎ𝑐 𝑓(𝜌2 ) 28( 0.9996 )
𝐸𝑆𝑈1 = = = 99.6752
𝐾 𝐷2 0.0104(27)
Para 𝐸𝐼𝑁𝐷2 :
ℎ𝑐 28
𝐸𝐼𝑁𝐷2 = = = 3.0584
(𝐾 𝐷2)2 (0.0104 ∗ 27)2
𝐶2 = 10 log10 ((2.0681)2 + (99.6752)2 + (355.1107)2 ) = 51.3367
Para la fase 𝐶:
12 ℎ𝑎 ℎ𝑐 12(1)(28)
= = 0.56
𝜆 600
29.6696 ≤ 0.56 (𝑁𝑂 𝐴𝑃𝐿𝐼𝐶𝐴)
29.6696 > 0.56
Por lo que de acuerdo a la segunda condición:
Para 𝐷𝑊2 :
2𝜋
𝐾= = 0.0104
600
28 12(1)(28)
𝐷𝑊2 = = 1.7127
0.0104(29.6696) 600 (29.6696)
Para 𝐸𝑆𝑈2 :
52.5 𝐷2 52.5 (29.6696)
𝜌2 = 2
= = 1.0817𝑒 −3
𝛿𝜆 4 (6002 )
2 + 0.3 𝜌2 2 + 0.3 (1.0817𝑒 −3 )
𝑓(𝜌1 ) = = = 0.9996
2 + 𝜌2 + 0.6 𝜌2 2 2 + 1.0817𝑒 −3 + 0.6 (1.0817𝑒 −3 )2
ℎ𝑐 𝑓(𝜌2 ) 28( 0.9996 )
𝐸𝑆𝑈1 = = = 90.7066
𝐾 𝐷2 0.0104(29.6696)
Para 𝐸𝐼𝑁𝐷2 :
ℎ𝑐 28
𝐸𝐼𝑁𝐷2 = = = 3.0584
(𝐾 𝐷2)2 (0.0104 ∗ 29.6696)2
𝐶2 = 10 log10 ((1.7127)2 + (90.7066)2 + (294.0816)2 ) = 49.7641
Calculando la 𝑅𝐼 por fase:

19
 14.13 2.682
𝑅𝐼𝐵𝐴−𝐹𝐴𝑆𝐸 𝐴 = 46 + 120 log10 ( ) + 40 log10 ( )+
17.56 3.51
2000
10(1 − (log10 (10 ∗ 0.5))2 ) + − 49.4137 + 49.7641 = 𝟒𝟐. 𝟏𝟑𝟏𝟔 𝒅𝑩
300
15.36 2.682
𝑅𝐼𝐵𝐴−𝐹𝐴𝑆𝐸 𝐵 = 46 + 120 log10 ( ) + 40 log10 ( )+
17.56 3.51
2000
10(1 − (log10 (10 ∗ 0.5))2 ) + − 49.4137 + 51.3367 = 𝟒𝟔. 𝟒𝟖𝟏𝟓 𝒅𝑩
300
14.13 2.682
𝑅𝐼𝐵𝐴−𝐹𝐴𝑆𝐸 𝐶 = 46 + 120 log10 ( ) + 40 log10 ( )+
17.56 3.51
2000
10(1 − (log10 (10 ∗ 0.5))2 ) + − 49.4137 + 49.7641 = 𝟒𝟐. 𝟏𝟑𝟏𝟔 𝒅𝑩
300
Con los resultados anteriores se realiza el análisis de la magnitud de cada resultado,
ya que se debe de ordenar de manera descendiente para poder obtener la RI
resultante de acuerdo a la formula 𝐶𝐼𝑆𝑃𝑅:
𝑆𝑖 𝑅𝐼1 − 𝑅𝐼2 ≥ 3 𝑑𝐵 𝑜 𝑆𝑖 𝑅𝐼1 − 𝑅𝐼2 ≤ 3 𝑑𝐵
𝑅𝐼1 + 𝑅𝐼2
𝑅𝐼𝑇𝑂𝑇𝐴𝐿 = 𝑅𝐼1 𝑅𝐼𝑇𝑂𝑇𝐴𝐿 = + 1.5
2
Por lo que lo que ordenando los resultados:
𝑅𝐼1 = 46.4815 𝑑𝐵
𝑅𝐼2 = 42.1316 𝑑𝐵
𝑅𝐼3 = 42.1316 𝑑𝐵
Llevando el análisis:
𝑆𝑖 46.4815 − 46.111 ≥ 3 𝑑𝐵 𝑜 𝑆𝑖46.4815 − 46.111 ≤ 3 𝑑𝐵
4.3499 ≥ 3 𝑑𝐵 4.3499 ≤ 3 𝑑𝐵

𝑅𝐼𝑇𝑂𝑇𝐴𝐿 = 𝑅𝐼1 𝑅𝐼𝑇𝑂𝑇𝐴𝐿 = (𝑁𝑂 𝐴𝑃𝐿𝐼𝐶𝐴)

𝑅𝐼𝑇𝑂𝑇𝐴𝐿 = 𝟒𝟔. 𝟒𝟖𝟏𝟓
De lo anterior se concluye:
𝑅𝐼𝐵𝐴 = 𝟒𝟔. 𝟒𝟖𝟏𝟓 𝒅𝑩
𝑅𝐼𝐿𝐿𝑀 = 46.4815 + 17 = 𝟔𝟑. 𝟒𝟖𝟏𝟓 𝒅𝑩
𝑅𝐼𝐿𝐿𝐷 = 46.4815 + 24 = 𝟕𝟎. 𝟒𝟖𝟏𝟓 𝒅𝑩

20
Asimismo, el cálculo se realiza para cada punto del vector 𝑥𝑚 .

C. CREIPI
Para el cálculo de RI, se emplea la formula siguiente:
ℎ
𝑪𝑹𝑬𝑰𝑷𝑰: 𝑅𝐼𝐵𝐴 = 4 + (3.7 𝑔 − 12.2) + 40 log10 (𝑑) + 20 log10 ( 2 )
𝐷
2
𝑞
−12 (log10 (𝑓)) − 17 log10 (𝑓) +
300
Siendo la ecuación anterior para ambiente adecuado, por lo que para regular los
resultados en condiciones de lluvias moderadas y lluvias densas se considera lo
siguiente:
𝐿𝑙𝑢𝑣𝑖𝑎 𝑚𝑜𝑑𝑒𝑟𝑎𝑑𝑎: 𝑅𝐼𝐿𝐿𝑀 = 𝑅𝐼𝐵𝐴 + 17
𝐿𝑙𝑢𝑣𝑖𝑎 𝑑𝑒𝑛𝑠𝑎: 𝑅𝐼𝐿𝐿𝐷 = 𝑅𝐼𝐵𝐴 + 24
Siendo:

𝑔𝐴 = 14.13 𝑘𝑉 𝑟𝑚𝑠 ⁄𝑐𝑚 𝑞 = 2000 𝑚𝑠𝑛𝑚

𝑔𝐵 = 15.36 𝑘𝑉 𝑟𝑚𝑠 ⁄𝑐𝑚

𝑔𝐶 = 14.13 𝑘𝑉 𝑟𝑚𝑠 ⁄𝑐𝑚

𝑑 = 2 ∗ 1.341 𝑐𝑚 = 2.682 𝑐𝑚 (𝑑𝑖𝑎𝑚𝑒𝑡𝑟𝑜 𝑑𝑒𝑙 𝑐𝑜𝑛𝑑𝑢𝑐𝑡𝑜𝑟)
𝑓 = 0.5 𝑀ℎ𝑧

Calculando las distancias radiales:

𝐷 = √(𝑦𝑚 − 𝑦𝑖 )2 + (𝑥𝑚 − 𝑥𝑖 )2
𝑦𝑚 = 𝐴𝑙𝑡𝑢𝑟𝑎 𝑑𝑒 𝑙𝑎𝑠 𝑓𝑎𝑠𝑒𝑠
𝑥𝑖 = 𝐷𝑖𝑠𝑡𝑎𝑛𝑐𝑖𝑎𝑠 𝑙𝑎𝑡𝑒𝑟𝑎𝑙𝑒𝑠 𝑑𝑒 𝑙𝑎𝑠 𝑓𝑎𝑠𝑒𝑠
Siendo:
𝑦𝐴 = ℎ𝑐1 = 28 𝑚 𝑥𝑚 = −50: 1: 50 𝑚
𝑦𝐵 = ℎ𝑐2 = 28 𝑚 𝑥𝐴 = −12.3 𝑚
𝑦𝐵 = ℎ𝑐3 = 28 𝑚 𝑥𝐵 = 0 𝑚
𝑦𝑖 = ℎ𝑎 = 1 𝑚 𝑥𝐴 = 12.3 𝑚
Para cada fase y con el primero valor del vector 𝑥𝑚 :

21
 𝐷𝐴 = √(28 − 1)2 + (0 − (−12.3))2 = 29.6696 𝑚

𝐷𝐵 = √(28 − 1)2 + (0 − 0)2 = 27 𝑚

𝐷𝐶 = √(28 − 1)2 + (0 − (12.3))2 = 29.6696 𝑚

Calculando la 𝑅𝐼 por fase:
28
𝑅𝐼𝐵𝐴−𝐹𝐴𝑆𝐸 𝐴 = 4 + (3.7 (14.13) − 12.2) + 40 log10 (2.682) + 20 log10 ( )
(29.6696)2
2000
− 12 (log10 (0.5))2 − 17 log10 (0.5) + = 𝟒𝟏. 𝟗𝟔𝟔𝟕 𝒅𝒃
300
28
𝑅𝐼𝐵𝐴−𝐹𝐴𝑆𝐸 𝐵 = 4 + (3.7 (15.36 ) − 12.2) + 40 log10 (2.682) + 20 log10 ( )
(27)2
2000
− 12 (log10 (0.5))2 − 17 log10 (0.5) + = 𝟒𝟖. 𝟏𝟓𝟓𝟕 𝒅𝒃
300
28
𝑅𝐼𝐵𝐴−𝐹𝐴𝑆𝐸 𝐶 = 4 + (3.7 (14.13) − 12.2) + 40 log10 (2.682) + 20 log10 ( )
(29.6696)2
2000
− 12 (log10 (0.5))2 − 17 log10 (0.5) + = 𝟒𝟏. 𝟗𝟔𝟔𝟕 𝒅𝒃
300
Con los resultados anteriores se realiza el análisis de la magnitud de cada resultado,
ya que se debe de ordenar de manera descendiente para poder obtener la RI
resultante de acuerdo a la formula 𝐶𝐼𝑆𝑃𝑅:
𝑆𝑖 𝑅𝐼1 − 𝑅𝐼2 ≥ 3 𝑑𝐵 𝑜 𝑆𝑖 𝑅𝐼1 − 𝑅𝐼2 ≤ 3 𝑑𝐵
𝑅𝐼1 + 𝑅𝐼2
𝑅𝐼𝑇𝑂𝑇𝐴𝐿 = 𝑅𝐼1 𝑅𝐼𝑇𝑂𝑇𝐴𝐿 = + 1.5
2
Por lo que lo que ordenando los resultados:
𝑅𝐼1 = 48.1557 𝑑𝐵
𝑅𝐼2 = 41.9667 𝑑𝐵
𝑅𝐼3 = 41.9667𝑑𝐵
Llevando el análisis:
𝑆𝑖 48.1557 − 41.9667 ≥ 3 𝑑𝐵 𝑜 𝑆𝑖 48.1557 − 41.9667 ≤ 3 𝑑𝐵
6.1890 ≥ 3 𝑑𝐵 6.1890 ≤ 3 𝑑𝐵

𝑅𝐼𝑇𝑂𝑇𝐴𝐿 = 𝑅𝐼1 𝑅𝐼𝑇𝑂𝑇𝐴𝐿 = (𝑁𝑂 𝐴𝑃𝐿𝐼𝐶𝐴)

𝑅𝐼𝑇𝑂𝑇𝐴𝐿 = 𝟒𝟖. 𝟏𝟓𝟓𝟕
De lo anterior se concluye:

22
𝑅𝐼𝐵𝐴 = 𝟒𝟖. 𝟏𝟓𝟓𝟕 𝒅𝑩
𝑅𝐼𝐿𝐿𝑀 = 46.4815 + 17 = 𝟔𝟓. 𝟏𝟓𝟓𝟕 𝒅𝑩
𝑅𝐼𝐿𝐿𝐷 = 46.4815 + 24 = 𝟕𝟐. 𝟏𝟓𝟓𝟕 𝒅𝑩

Asimismo, el cálculo se realiza para cada punto del vector 𝑥𝑚 .

23
Apéndice C: Programa para el cálculo de perfiles de radio interferencia en líneas
de transmisión convencionales y de alta cargabilidad (HSIL).

Datos caso 1
%HSIL_simetrica_4_cond

V=500;
Fases=3;
Dc=2*1.341;
Rsub=0.3182;
f=0.5;
v=3e8;
x0=50;

for k=1:Fases

if k==1
g(k)=14.13;
end

if k==2
g(k)=15.36;
end

if k==3
g(k)=14.13;
end
end
xlat=12.3;

for k=1:Fases

h(k)=28;

if k==1
x(k)=-xlat;
end
if k==2
x(k)=0;

24
 end
if k==3
x(k)=xlat;
end

end

Datos caso 2
%HSIL_simetrica_6_cond

V=500;
Fases=3;
Dc=2*1.083;
Rsub=0.75;
xlat=3.35;
f=0.5;
v=3e8;
x0=50;

for k=1:Fases

if k==1
x(k)=-xlat;
h(k)=27.5+5.8023;
g(k)=20.28;
end
if k==2
x(k)=0;
h(k)=27.5;
g(k)=20.53;
end
if k==3
x(k)=xlat;
h(k)=27.5+5.8023;
g(k)=20.28;
end

end

25
Datos caso 3
%HSIL_asimetrica_6_cond

V=500;
Fases=3;
Dc=2*1.341;
Rsub=0.75;

xlat=9.2;
f=0.5;
v=3e8;
x0=50;

for k=1:Fases

h(k)=28;
if k==1
x(k)=-xlat;
g(k)=17.21;
end
if k==2
x(k)=0;
g(k)=17.60;
end
if k==3
x(k)=xlat;
g(k)=17.21;
end

end

Datos caso 4
%HSIL_asimetrica_7_cond

V=500;
Fases=3;
Dc=2*1.083;
Rsub=0.75;

26
xlat=7.2;
f=0.5;
v=3e8;
x0=50;

for k=1:Fases

h(k)=28;
if k==1
x(k)=-xlat;
g(k)=17.53;
end
if k==2
x(k)=0;
g(k)=17.53;
end
if k==3
x(k)=xlat;
g(k)=17.53;
end

end

Programa de la radio interferencia en líneas de

transmisión (CIGREE)
clc,clear all

%HSIL_CASO1_4_cond
%HSIL_CASO2_6_cond
%HSIL_CASO3_6_cond
HSIL_CASO4_7_cond

x1=-[Link];
y1=1;
q=2000/300;
% calculo de las distancias radiales
for k1=1 : length(x1)

27
 for m=1:Fases

D(m)= sqrt ((h(m)-y1)^2 + (x1(k1) - x(m))^2);

if k1==51
fprintf('\nLa distancia radial: %f',D(m));
end
%Calculo de radio interferencia
%fprintf("\n\nLa radio interferencia entre fases es la siguiente\n")
RI(m)=-10+(3.5*g(m))+(6*Dc)-(33*log10(D(m)/20))+q;

end

B = sort(RI,'descend');
%fprintf('\nLa radio interferencia de la fase A es RI_1 = %f',B(1));
%fprintf('\nLa radio interferencia de la fase B es RI_2 = %f',B(2));
%fprintf('\nLa radio interferencia de la fase C es RI_3 = %f',B(3));

if B(1)-B(2)>3
RA_T= B(1);
%RA_T
else
RA_T=((B(1)+B(2))/2)+1.5;
end
R_total_LLD(k1)=RA_T;

end
%mean(RA_T);
R_BA=R_total_LLD-24;
R_LLM=R_total_LLD-7;
plot( x1, R_total_LLD , x1, R_BA , x1, R_LLM,'LineWidth',1.4), grid;
yline(60,'--k');
xline(-25,'-.g');
xline(25,'-.g');
legend('Lluvia densa','Ambiente favorable','Lluvia moderada','Limite
regulatorio','Derecho de vía','Location','south');
hold off
title('Radio interferencia');
xlabel('Distancia Lateral (m)');
ylabel('Radio interferencia (dB (µV/m))');
fontsize(8, 'points');
xticks(-[Link]);
%ylim([10 65]);
ylim([10 80]); %Para caso 2
yticks([Link]);

28
fprintf('\nLa radio interferencia maxima en buen ambiente: %f',R_BA(51));
fprintf('\nLa radio interferencia maxima en lluvia moderada: %f',R_LLM(51));
fprintf('\nLa radio interferencia maxima en lluvia densa:
%f',R_total_LLD(51));

Programa de la radio interferencia en líneas de

transmisión (BPA)
clc,clear all

%HSIL_CASO1_4_cond
%HSIL_CASO2_6_cond
%HSIL_CASO3_6_cond
HSIL_CASO4_7_cond

x1=-[Link];
y1=1;
q=2000/300;
lamb=v/(f*1e6);
K=(2*pi)/lamb;
% Parámetros de referencia para C1
% Ca=10*log10((DW)^2+(ESU)^2+(EIND)^2);

DW1=31.1/f;
D1=21.04;
EIND1=70.55/((f)^2);
p1=(52.5*D1)/(4*lamb^2);
fp1=(2+(0.3*p1))/(2+p1+(0.6*(p1)^2));
ESU1=(31.1*fp1)/f;
C1=10*log10((DW1)^2+(ESU1)^2+(EIND1)^2);

for k1=1 : length(x1)

for m=1:Fases

% Para el cálculo de C2:

% calculo de las distancias radiales para C2
D2= sqrt ((h(m)-y1)^2 + (x1(k1) - x(m))^2);

% Para DW2
if D2<=((12*h(m)*y1)/lamb)

29
 DW2=h(m)/(K*D2);
else
DW2=(h(m)/(K*D2))*(((12*h(m)*y1))/(lamb*D2));
end

% Para ESU2
p2=(52.5*D2)/(4*lamb^2);
fp2=(2+(0.3*p2))/(2+p2+(0.6*(p2)^2));
ESU2=(fp2*h(m))/(K*D2);

% Para EIND2
EIND2=h(m)/(K*D2)^2;

% C2:
C2=10*log10((DW2)^2+(ESU2)^2+(EIND2)^2);
%if k1==51
%fprintf('\nLa distancia radial: %f',D2);
%fprintf('\nLa DW2: %f',DW2);
%fprintf('\nLa ESU2: %f',ESU2);
%fprintf('\nLa EIND2: %f',EIND2);
%fprintf('\nLa C2: %f',C2);
%end

%Calculo de radio interferencia

%fprintf("\n\nLa radio interferencia entre fases es la siguiente\n")
RI(m)=46+(120*log10(g(m)/17.56))+(40*log10(Dc/3.51))+(10*(1-
(log10(10*f))^2))+q-C1+C2;

end

B = sort(RI,'descend');
%if k1==51
%fprintf('\nLa radio interferencia de la fase A es RI_1 = %f',B(1));
%fprintf('\nLa radio interferencia de la fase B es RI_2 = %f',B(2));
%fprintf('\nLa radio interferencia de la fase C es RI_3 = %f',B(3));
%end

if B(1)-B(2)>3
RA_T= B(1);
%RA_T
else
RA_T=((B(1)+B(2))/2)+1.5;

30
end
R_total_BA(k1)=RA_T;

end
%mean(RA_T);
R_LLM=R_total_BA+17;
R_LLD=R_total_BA+24;
plot( x1, R_LLD , x1, R_total_BA , x1, R_LLM,'LineWidth',1.4), grid;
yline(60,'--k');
xline(-25,'-.g');
xline(25,'-.g');
legend('Lluvia densa','Ambiente favorable','Lluvia moderada','Limite
regulatorio','Derecho de vía','Location','south');
hold
title('Radio interferencia');
xlabel('Distancia Lateral (m)');
ylabel('Radio interferencia (dB (µV/m))');
fontsize(8, 'points');
xticks(-[Link]);
%ylim([10 65]);
ylim([10 85]); %Para caso 2
yticks([Link]);
fprintf('\nLa radio interferencia maxima en buen ambiente:
%f',R_total_BA(51));
fprintf('\nLa radio interferencia maxima en lluvia moderada: %f',R_LLM(51));
fprintf('\nLa radio interferencia maxima en lluvia densa: %f',R_LLD(51));

Programa de la radio interferencia en líneas de

transmisión (CREIPI)
clc,clear all

%HSIL_CASO1_4_cond
%HSIL_CASO2_6_cond
%HSIL_CASO3_6_cond
HSIL_CASO4_7_cond

x1=-[Link];
y1=1;
q=2000/300;

for k1=1 : length(x1)

31
 for m=1:Fases

% calculo de las distancias radiales para

D(m)= sqrt ((h(m)-y1)^2 + (x1(k1) - x(m))^2);
if k1==51
fprintf('\nLa distancia radial: %f',D(m));
end
%Calculo de radio interferencia
RI(m)=4+((3.7*g(m))-12.2)+(40*log10(Dc))+(20*log10(h(m)/(D(m)^2)))-
(12*(log10(f))^2)-(17*log10(f))+q;

end

B = sort(RI,'descend');
if k1==51
fprintf('\nLa radio interferencia de la fase A es RI_1 = %f',B(1));
fprintf('\nLa radio interferencia de la fase B es RI_2 = %f',B(2));
fprintf('\nLa radio interferencia de la fase C es RI_3 = %f',B(3));
end
if B(1)-B(2)>3
RA_T= B(1);
%RA_T
else
RA_T=((B(1)+B(2))/2)+1.5;
end
R_total_BA(k1)=RA_T;

end
%mean(RA_T);
R_LLM=R_total_BA+17;
R_LLD=R_total_BA+24;
plot( x1, R_LLD , x1, R_total_BA , x1, R_LLM,'LineWidth',1.4), grid;
yline(60,'--k');
xline(-25,'-.g');
xline(25,'-.g');
legend('Lluvia densa','Ambiente favorable','Lluvia moderada','Limite
regulatorio','Derecho de vía','Location','south');
hold
title('Radio interferencia');
xlabel('Distancia Lateral (m)');
ylabel('Radio interferencia (dB (µV/m))');
fontsize(8, 'points');
xticks(-[Link]);

32
%ylim([10 65]);
ylim([10 90]); %Para caso 2
yticks([Link]);
%fprintf('\nLa distancia radial 0',D(51));
fprintf('\nLa radio interferencia maxima en buen ambiente:
%f',R_total_BA(51));
fprintf('\nLa radio interferencia maxima en lluvia moderada: %f',R_LLM(51));
fprintf('\nLa radio interferencia maxima en lluvia densa: %f',R_LLD(51));

33
Apéndice D: Simulación en ATP Draw para la obtención de los parámetros de las
torres de transmisión para el cálculo del SIL.
- Caso 1.

Figura 63. Ícono de la torre LCC

template

Figura 64. Modelado de la torre - caso 1

34
 Figura 65. Parámetros de la estructura de la torre caso 1

Figura 66. Parámetros obtenidos caso 1

De acuerdo a la ecuación 1.2, y con los parámetros obtenidos en la simulación de ATP Draw, se obtiene
el siguiente valor de impedancia característica:

0.309
𝑍𝐶 = √ = 248.6109 Ω
4.9994 ∗ 10−6

Con el valor de la impedancia característica, sustituimos de acuerdo a la ecuación 1.3 para obtener el
valor SIL de la línea de transmisión, siendo:

35
 (500 𝑘𝑉)2
𝑆𝐼𝐿 = = 1005.5871 𝑀𝑊
248.6109 Ω
- Caso 2.

Figura 67. Modelado de la torre - caso 2

36
 Figura 68. Parámetros de la estructura de la torre caso 2

Figura 69. Parámetros obtenidos caso 2

De acuerdo a la ecuación 1.2, y con los parámetros obtenidos en la simulación de ATP Draw, se obtiene
el siguiente valor de impedancia característica:

0.184
𝑍𝐶 = √ = 149.3311 Ω
8.2512 ∗ 10−6

37
Con el valor de la impedancia característica, sustituimos de acuerdo a la ecuación 1.3 para obtener el
valor SIL de la línea de transmisión, siendo:

(500 𝑘𝑉)2
𝑆𝐼𝐿 = = 1674.1318 𝑀𝑊
149.3311 Ω

- Caso 3.

Figura 70. Modelado de la torre - caso 3

38

Figura 71. Parámetros de la estructura de la torre caso 3

De acuerdo a la ecuación 1.2, y con los parámetros obtenidos en la simulación de ATP Draw, se obtiene
el siguiente valor de impedancia característica:

0.2561
𝑍𝐶 = √ = 207.45 Ω
5.9509 ∗ 10−6

Con el valor de la impedancia característica, sustituimos de acuerdo a la ecuación 1.3 para obtener el
valor SIL de la línea de transmisión, siendo:

(500 𝑘𝑉)2
𝑆𝐼𝐿 = = 1205.109 𝑀𝑊
217.5782 Ω

- Caso 4.

39
 Figura [Link]ámetros de la estructura de la torre caso 4

Figura 75. Parámetros obtenidos caso 4

De acuerdo a la ecuación 1.2, y con los parámetros obtenidos en la simulación de ATP Draw, se obtiene
el siguiente valor de impedancia característica:

40
 0.2058
𝑍𝐶 = √ = 164.5189 Ω
7.6035 ∗ 10−6

Con el valor de la impedancia característica, sustituimos de acuerdo a la ecuación 1.3 para obtener el
valor SIL de la línea de transmisión, siendo:

(500 𝑘𝑉)2
𝑆𝐼𝐿 = = 1519.5818 𝑀𝑊
186.1861 Ω

41