Metodología para El Análisis de Transitorios Electromagnéticos de Maniobra en Sistemas de Alta t2

Metodología para el análisis de transitorios electromagnéticos de maniobra en sistemas de
alta tensión (230 kV/500 kV) con ayuda del EMTP/ATP Draw
Cristian C. ARIAS
Trabajo de grado para optar al título de Ingeniero Electricista
Director
Johan Andrés Sánchez Patiño
Especialista en Transmisión y Distribución de Energía Eléctrica
Magister en Ingeniería
Universidad Pontificia Bolivariana

Ingeniería Eléctrica
Medellín
2015
2
Arias (2015)
Contenido
INTRODUCCIÓN ...................................................................................................................................................... 12
1. CONCEPTOS BÁSICOS .......................................................................................................................... 13
1.1 Sobre el modelado de sistemas de potencia y los fenómenos transitorios ................................................ 13
1.2 El software ATP ......................................................................................................................................... 13
1.3 Aspectos importantes a considerar en la modelación de redes eléctricas ................................................ 14
1.4 Ondas viajeras ........................................................................................................................................... 16
1.5 Maniobras ejecutadas en Sistemas Eléctricos de Potencia. ...................................................................... 21
2. METODOLOGÍA PARA LA SIMULACIÓN DE TRANSITORIOS DE MANIOBRA CON AYUDA
DEL ATP/EMTP Y SU INTERFAZ GRÁFICA ATPDraw....................................................................................... 22
2.1 Información básica del sistema y las maniobras ....................................................................................... 22
2.2 Modelos implementados para la elaboración de la red ............................................................................ 26
2.3 Energización de líneas ............................................................................................................................... 39
2.4 Energización de transformador 220 kV/500 kV......................................................................................... 55
2.5 Recierre monofásico de líneas ................................................................................................................... 69
2.6 Análisis del arco secundario en línea a 500 kV ......................................................................................... 78
2.7 Verificación de TRV en interruptores de línea .......................................................................................... 89
CONLCLUSIONES ................................................................................................................................................... 99
ANEXOS .................................................................................................................................................................. 103
Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

3
Arias (2015)
Lista de Figuras
Figura 1. Representación de un pequeño tramo de una línea de transmisión ........................................................................................................................ 17
Figura 2. Reflexión y Refracción de la onda incidente (a) .................................................................................................................................................... 20
Figura 3. Reflexión y Refracción de la onda incidente (b) .................................................................................................................................................... 20
Figura 4. Reflexión y Refracción de la onda incidente (c) .................................................................................................................................................... 21
Figura 5. Reflexión y Refracción de la onda incidente (d) .................................................................................................................................................... 21
Figura 6. Red implementada en el ATP................................................................................................................................................................................. 23
Figura 7. Ajustes del ATP para las simulaciones .................................................................................................................................................................. 25
Figura 8. Localización del modelo de línea .......................................................................................................................................................................... 26
Figura 9. Ajustes del modelo de línea .................................................................................................................................................................................. 27
Figura 10. Representación gráfica del modelo de línea antes de ser configurado ................................................................................................................ 27
Figura 11. Representación gráfica del modelo final de la línea............................................................................................................................................ 27
Figura 12. Representación gráfica del modelo BCTRAN ..................................................................................................................................................... 29
Figura 13. Localización del modelo BCTRAN ..................................................................................................................................................................... 29
Figura 14. Ajustes y configuración del modelo BCTRAN .................................................................................................................................................... 30
Figura 15. Localización del inductor no lineal para representación de la saturación ........................................................................................................... 30
Figura 16. Representación gráfica del inductor no lineal tipo 98 ......................................................................................................................................... 31
Figura 17. Modelo final de transformador más saturación .................................................................................................................................................... 31
Figura 18. Ajustes del inductor no lineal tipo 98 .................................................................................................................................................................. 31
Figura 19. Característica de saturación ingresada en el inductor no lineal tipo 98 ............................................................................................................... 32
Figura 20. Representación de la resistencia no lineal tipo 92................................................................................................................................................ 32
Figura 21. Localización resistencia no lineal tipo 92 ........................................................................................................................................................... 33
Figura 22. Ajustes de la resistencia no lineal tipo 92 ........................................................................................................................................................... 33
Figura 23. Característica de maniobra del descargador ingresada en la resistencia no lineal tipo 92 ................................................................................... 34
Figura 24. Modelo final de descargador de sobretensiones .................................................................................................................................................. 34
Figura 25. Localización de la fuente Type 14....................................................................................................................................................................... 35
Figura 26. Ajustes de la fuente Type 14 ............................................................................................................................................................................... 35
Figura 27. Localización del Symmetric RL coupled line ...................................................................................................................................................... 36
Figura 28. Representación gráfica de la fuente Type 14....................................................................................................................................................... 36
Figura 29. Representación gráfica del Symmetric RL coupled line ...................................................................................................................................... 36

4
Arias (2015)
Figura 30. Ajustes del Symmetric RL coupled line................................................................................................................................................................ 37

Figura 31. Localización del modelo de compensación capacitiva........................................................................................................................................ 37
Figura 32. Representación de la compensación capacitiva luego de su configuración ......................................................................................................... 38
Figura 33. Localización del modelo de compensación inductiva ......................................................................................................................................... 38
Figura 34. Representación de la compensación inductiva luego de su configuración .......................................................................................................... 39
Figura 35. Localización de interruptores estadísticos en el menú rápido de elementos del ATP .......................................................................................... 39
Figura 36. Ajustes del ATP Settings para simulación con interruptores estadísticos (a) ....................................................................................................... 40
Figura 37. Ajustes del ATP Settings para simulación con interruptores estadísticos (b) ....................................................................................................... 41
Figura 38. Configuración del interruptor “Master” (Fase A) ............................................................................................................................................... 42
Figura 39. Configuración del interruptor “Slave” (Fase B) .................................................................................................................................................. 42
Figura 40. Configuración del interruptor “Slave” (Fase C) .................................................................................................................................................. 43
Figura 41. Ubicación del Splitter en el menú rápido de elementos del ATP ......................................................................................................................... 43
Figura 42. Splitter trifásico ................................................................................................................................................................................................... 44
Figura 43. Conexión del interruptor estadístico para efectos de la simulación ..................................................................................................................... 44
Figura 44. Modelo de línea más interruptor estadístico para simulación de energización de línea ....................................................................................... 44
Figura 45. Ejemplo de archivo .lis ........................................................................................................................................................................................ 45
Figura 46. Resultados de energización en el archivo .lis....................................................................................................................................................... 45
Figura 47. Correspondencia de resultados en el archivo .lis ................................................................................................................................................. 46
Figura 48. Configuración del resistor no lineal correspondiente al modelo de descargador de sobretensiones para registrar corriente ............................... 46
Figura 49. Modelo de línea más interruptor controlado por tiempo para simulación del caso crítico ................................................................................... 47
Figura 50. Selección del modelo de interruptor controlado por tiempo ................................................................................................................................ 47
Figura 51. Resultados de la energización de la línea 500 kV en la que se obtuvo la mayor sobretensión ............................................................................. 48
Figura 52. Resultados de la energización de la línea 220 kV en la que se obtuvo la mayor sobretensión ............................................................................. 48
Figura 53. Identificación de los polos del interruptor en el archivo .lis para asignación correcta de tiempos ....................................................................... 48
Figura 54. Configuración del interruptor con los tiempos críticos para energización de línea a 500 kV .............................................................................. 49
Figura 55. Configuración del interruptor con los tiempos críticos para energización de línea a 220 kV .............................................................................. 49
Figura 56. Tensión durante energización de línea a 500 kV ................................................................................................................................................. 50
Figura 57. Energía en descargadores de sobretensión de línea a 500 kV .............................................................................................................................. 50
Figura 58. Tensión durante energización de línea a 220 kV ................................................................................................................................................. 51
Figura 60. Modelo para energización con resistencias de preinserción ................................................................................................................................. 52

5
Arias (2015)
Figura 61. Configuración del interruptor 1 para energización con resistencias de preinserción ........................................................................................... 52
Figura 62. Configuración del interruptor 2 para energización con resistencias de preinserción ........................................................................................... 53
Figura 63. Selección del modelo de resistencia de preinserción ........................................................................................................................................... 53
Figura 64. Resistencia de preinserción .................................................................................................................................................................................. 53
Figura 65. Configuración del modelo de resistencia de preinserción para energización ....................................................................................................... 54
Figura 66. Tensión durante energización de línea a 500 kV con resistencias de preinserción .............................................................................................. 55
Figura 68. Modelo para Energización de transformador ....................................................................................................................................................... 56
Figura 69. Configuración de las fases del interruptor estadístico para monitorear corriente ................................................................................................. 57
Figura 70. Modelo transformador más interruptor controlado por tiempo para simulación de casos críticos ....................................................................... 58
Figura 71. Resultados de la energización del transformador desde el lado de alta en la que se obtuvo la mayor sobretensión ............................................ 58
Figura 72. Resultados de la energización del transformador desde el lado de baja en la que se obtuvo la mayor sobretensión............................................ 58
Figura 73. Configuración del interruptor con los tiempos críticos para energización del transformador desde el lado de alta ............................................. 59
Figura 74. Configuración del interruptor con los tiempos críticos para energización del transformador desde el lado de baja ............................................ 59
Figura 75. Resultados de la energización del transformador desde el lado de alta en la que se obtuvo la mayor corriente Inrush ....................................... 60
Figura 76. Resultados de la energización del transformador desde el lado de baja en la que se obtuvo la mayor corriente Inrush ...................................... 60
Figura 77. Configuración del interruptor con los tiempos críticos para energización del transformador desde el lado de alta ............................................. 60
Figura 78. Configuración del interruptor con los tiempos críticos para energización del transformador desde el lado de baja ............................................ 61
Figura 79. Tensión del transformador en el lado de alta al energizar desde el lado de baja .................................................................................................. 62
Figura 80. Energía en descargadores de sobretensión lado de alta ........................................................................................................................................ 62
Figura 81. Tensión en el lado de baja del transformador al energizar desde el lado de alta .................................................................................................. 63
Figura 82. Energia en descargadores de sobretensión lado de baja ....................................................................................................................................... 64
Figura 83. Corriente Inrush del transformador al energizar desde alta .................................................................................................................................. 64
Figura 84. Corriente Inrush del transformador al energizar desde baja ................................................................................................................................. 64
Figura 85. Tiempos óptimos para el mando sincronizado, lado 220 kV ............................................................................................................................... 65
Figura 86. Tiempos óptimos para el mando sincronizado, lado 500 kV ............................................................................................................................... 65
Figura 87. Configuración del interruptor para energización con mando sincronizado desde alta ......................................................................................... 66
Figura 88. Configuración del interruptor para energización con mando sincronizado desde baja ........................................................................................ 66
Figura 89. Tensión en el secundario del transformador al energizar desde el lado de alta .................................................................................................... 67
Figura 90. Corriente Inrush del transformador al energizar desde Alta ................................................................................................................................ 67
Figura 91. Tensión en el lado de alta del transformador al energizar desde el lado de baja .................................................................................................. 67

6
Arias (2015)
Figura 92. Corriente inrush del transformador al energizar desde baja ................................................................................................................................. 68
Figura 93. Energía en descargadores de sobretensión del lado de alta del transformador ..................................................................................................... 69
Figura 94. Energía en descargadores de sobretensión del lado de baja del transformador .................................................................................................... 69
Figura 95. Modelo para recierre monofásico de línea ........................................................................................................................................................... 70
Figura 96. Configuración del Interruptor 1 ........................................................................................................................................................................... 70
Figura 97. Configuración del Interruptor 2 ........................................................................................................................................................................... 71
Figura 98. Configuración de la Fase A encargada de la maniobra de recierre estadístico ..................................................................................................... 71
Figura 99. Selección de los interruptores monofásicos ......................................................................................................................................................... 72
Figura 100. Configuración de los interruptores monofásicos ................................................................................................................................................ 72
Figura 101. Modelo de línea e interruptor para recierre monofásico con tiempo crítico ....................................................................................................... 73
Figura 102. Resultados del recierre monofásico en la línea a 500 kV en la que se obtuvo la mayor sobretensión ............................................................... 74
Figura 103. Resultados del recierre monofásico en la línea a 220 kV en la que se obtuvo la mayor sobretensión ............................................................... 74
Figura 104. Configuración del interruptor monofásico controlado por tiempo para recierre monofásico con tiempo crítico en línea a 500 kV .................. 74
Figura 105. Configuración del interruptor monofásico controlado por tiempo para recierre monofásico con tiempo crítico en línea a 220 kV .................. 75
Figura 106. Tensión en las tres fases línea, 500 kV .............................................................................................................................................................. 76
Figura 107. Tensión en la fase en recierre, línea 500 kV ...................................................................................................................................................... 76
Figura 108. Energía en descargadores, línea 500 kV ............................................................................................................................................................ 76
Figura 109. Tensión en las tres fases, línea 220 kV .............................................................................................................................................................. 76
Figura 110. Tensión en la fase en recierre, línea 220 kV ...................................................................................................................................................... 77
Figura 111. Carga atrapada en la fase en recierre, línea 220 kV ........................................................................................................................................... 77
Figura 112. Energía en descargadores, línea 220 kV ............................................................................................................................................................ 77
Figura 113. Localización del MODEL en el menú rápido del ATP ...................................................................................................................................... 79
Figura 114. Representación gráfica del MODEL en el ATP ................................................................................................................................................. 79
Figura 115. Ventana de ajustes del MODEL......................................................................................................................................................................... 79
Figura 116. Editor de texto del MODEL ............................................................................................................................................................................... 80
Figura 117. Introducción del código en el editor de texto ..................................................................................................................................................... 81
Figura 118. Cuadro de diálogo emergente luego de salvar los cambios realizados en el editor de texto .............................................................................. 81
Figura 119. Ventana de edición de las variables del MODEL .............................................................................................................................................. 81
Figura 120. Representación gráfica final del MODEL.......................................................................................................................................................... 82
Figura 121. Localización de la resistencia TACS Type 91 en el menú rápido de elementos ................................................................................................ 82
Figura 122. Conexión del modelo de estudio de arco secundario a la línea de transmisión .................................................................................................. 83

7
Arias (2015)
Figura 123. Modelo para falla en extremo de la línea ........................................................................................................................................................... 84

Figura 124. Modelo de falla en la mitad de la línea .............................................................................................................................................................. 84
Figura 125. Ajustes del ATP para las simulaciones de arco secundario ................................................................................................................................ 85
Figura 126. Variables definidas para falla en extremo A ...................................................................................................................................................... 85
Figura 127. Variables definidas para falla en extremo B ...................................................................................................................................................... 86
Figura 128. Configuración interruptor asociado al extremo en falla ..................................................................................................................................... 86
Figura 129. Configuración del interruptor en el extremo remoto .......................................................................................................................................... 87
Figura 130. Configuración del interruptor auxiliar ............................................................................................................................................................... 87
Figura 131. Corriente de arco secundario. Falla en extremo A ............................................................................................................................................. 88
Figura 132. Corriente de arco secundario. Falla en extremo B ............................................................................................................................................. 88
Figura 133. Corriente de arco secundario. Falla en la mitad ................................................................................................................................................. 89
Figura 134. Localización de la resistencia que modela la falla trifásica y monofásica ......................................................................................................... 89
Figura 135. TRV Normalizado en interruptores .................................................................................................................................................................... 90
Figura 136. Ajustes de resistencia para falla monofásica sólida ........................................................................................................................................... 91
Figura 137. Ajustes de resistencia para falla trifásica sólida ................................................................................................................................................. 91
Figura 138. Ajustes de inerruptor para falla kilométrica ....................................................................................................................................................... 92
Figura 139. Ajustes de inerruptor para falla terminal ............................................................................................................................................................ 92
Figura 140. Modelo para verificación de TRV. Falla kilométrica ......................................................................................................................................... 93
Figura 141. Modelo para verificación de TRV. Falla terminal .............................................................................................................................................. 93
Figura 142. Ajustes del ATP Settings para obtener los valores de tensión punto a punto en el archivo .lis .......................................................................... 94
Figura 143. Información de la gráfica contenida en el archivo .lis........................................................................................................................................ 94
Figura 144. TRV interruptor línea a 220 kV, Fase B ............................................................................................................................................................. 95
Figura 145. TRV interruptor línea a 500 kV, Fase A ............................................................................................................................................................. 95
Figura 147. Punto de violación de la norma localizado en interruptor a 220 kV .................................................................................................................. 96
Figura 148. Punto de violación de la norma ampliado ......................................................................................................................................................... 96
Figura 150. Modelo interruptor más descargador para corrección del TRV .......................................................................................................................... 98
Figura 151. TRV interruptor línea a 500 kV con descargador de sobretensiones .................................................................................................................. 98
Figura 152. Modelo interruptor más capacitancias para corrección del TRV ........................................................................................................................ 98
Figura 153. TRV interruptor línea a 220 kV con capacitancias ............................................................................................................................................. 99

8
Arias (2015)
Figura 154. Punto de violación de la norma ampliado. Caso con capacitores ...................................................................................................................... 99
Figura 155. Configuración geométrica líneas a 500 kV ..................................................................................................................................................... 103
Figura 156. Configuración geométrica líneas 220 kV ......................................................................................................................................................... 103
Figura 157. Característica de protección a maniobra ......................................................................................................................................................... 104
Figura 158. Característica de saturación transformador 500 kV/220 kV/33 kV ................................................................................................................. 105
Figura 159. Característica de saturación transformador 500 kV/220 kV/33 kV en el ATP ................................................................................................ 105

9
Arias (2015)
Glosario
Aislamiento externo: son las distancias en el aire y las superficies en contacto con aire del aislamiento sólido del equipo, que están
sujetas a los esfuerzos dieléctricos y a los efectos atmosféricos y otras condiciones externas.
Alta tensión: todo suministro con tensión igual o superior a 115 kV.
Arco primario: descarga en un gas caracterizada por una alta corriente y una baja tensión
Arco secundario: corriente residual que fluye a través de una falla monofásica durante el tiempo muerto del recierre.
Descargador de Sobretensión: dispositivo para la protección del sistema de potencia y sus componentes contra las sobretensiones, ya sea
producidas por descargas atmosféricas o por maniobras en el sistema durante fallas.
Energización/puesta en servicio: procedimiento que se realiza para la toma de tensión y la toma de carga de los equipos y sistemas de la
subestación y de los circuitos asociados, para disponer en operación comercial la instalación.
Interruptor: dispositivo de maniobra capaz de interrumpir, establecer y llevar las corrientes normales o asignadas del circuito y las
anormales o de cortocircuito, mediante la conexión o desconexión de circuitos.
Onda viajera: perturbación energética ocasionada por un transitorio que se propaga a lo largo de un sistema de potencia o equipo.
Sistema sólidamente aterrizado: un sistema donde todos los puntos de neutro están aterrizados directamente.
Sobretensión temporal: sobretensión a frecuencia industrial de duración relativamente larga.
Sobretensión transitoria: sobretensión de corta duración de unos pocos milisegundos o menos, oscilatoria o no oscilatoria, por lo general
altamente amortiguada.
Sobretensión: cualquier tensión entre un conductor de fase y tierra o entre conductores de fase cuyo valor pico exceda el correspondiente
valor pico de la tensión más alta del equipo.
Tensión máxima del sistema: es la máxima tensión de operación que se puede presentar durante operación normal en cualquier momento
y en cualquier punto del sistema.

10
Arias (2015)
Resumen
En este artículo se explica la metodología para la modelación y simulación de transitorios de maniobra en sistemas
de alta tensión (230 kV/500 kV) utilizando el software ATP y su interfaz gráfica ATPDraw. Se definen tiempos,
parámetros y factores que pueden tener influencia sobre maniobras como energización de líneas, energización de
transformadores, recierres monofásicos y TRV en interruptores. De igual forma, se presentan algunas de las medidas
implementadas para mitigar estos transitorios. Copyright © UPB 2015
Palabras clave: Alta Tensión, Transitorios de maniobra, modelación, simulación.

11
Arias (2015)
Abstract
This paper explains the methodology of how to model and simulate switching transients in high voltage systems
(230 kV/500 kV); using the ATP software and its graphic interface ATPDraw. This methodology aims to define
time, parameters and factors that could have an influence in the line energization, transformers energization,
single-phase reswitching and TRV interrupters. At the same time, it has been shown some measures that have been
implemented to mitigate these transients.
Keywords: High voltage, switching transients, modeling, simulation.

12
Arias (2015)
INTRODUCCIÓN De maniobra, generados por apertura y cierre de

interruptores en maniobras de conexión, reconexión,
Un fenómeno transitorio en un Sistema Eléctrico de Potencia
(SEP) puede ser resultado de un cambio inesperado o planeado despeje de fallas, etc.
con anticipación en la operación o configuración del mismo [1]. Estos eventos (descargas atmosféricas y maniobras de
Estos fenómenos someten los sistemas a esfuerzos que en muchas interruptores) generan un desplazamiento de las ondas
ocasiones pueden afectar la vida útil de los equipos e inclusive electromagnéticas a través de las líneas de transmisión aéreas a
dañar el aislamiento eléctrico, debido a la gran disipación una velocidad aproximada de 300 m/µs, que aunque se van
energética que se produce en estos eventos. atenuando a medida que aumenta su recorrido, la superposición de
Los transitorios en los SEP se pueden clasificar en dos ondas incidentes y ondas reflejadas provoca que en ciertos puntos
categorías [2]: se generen elevaciones temporales de tensión [2].
Transitorios con una interacción entre la energía En este orden de ideas, todo sistema eléctrico debe de ser
magnética almacenada en los inductores y la energía diseñado considerando los esfuerzos más severos que pueden
presentarse durante su operación, de forma que se garantice una
eléctrica almacenada en capacitores.
protección adecuada de los equipos y el sistema pueda continuar
Transitorios con una interacción entre la energía suministrando energía eléctrica a las cargas de forma constante
mecánica almacenada en las máquinas rotatorias y la y sin que se vea comprometida la calidad de la misma. Por lo
energía eléctrica almacenada en el sistema eléctrico. anterior, la simulación de fenómenos transitorios con ayuda de
herramientas especializadas es una necesidad que se convierte en
Los primeros corresponden a transitorios de naturaleza obligación para poder elaborar un diseño y dimensionamiento
completamente eléctrica y son denominados transitorios apropiado de los componentes de un sistema eléctrico, la
electromagnéticos. Los pertenecientes a la segunda categoría son coordinación de aislamiento, coordinación de protecciones, entre
denominados transitorios electromecánicos. otras, además que estas simulaciones permiten obtener los
Según su origen, los transitorios también se clasifican en [3]: tiempos en los cuales se deben ejecutar las maniobras de modo
que sean lo menos dañinas posibles.
Atmosféricos, generados por descargas.
En el presente documento se definen los principales transitorios
electromagnéticos de maniobra que pueden presentarse en los

13
Arias (2015)
sistemas de alta tensión, ejecutándolos luego en un equivalente de transitorio resulta ser más complejo, debido al papel que juegan
red llevado al ATP, para mostrar una metodología que permite los elementos ante variaciones de frecuencia [4], por los
simular estos fenómenos de forma adecuada. acoplamientos y las no linealidades [5].
De forma práctica, el comportamiento de un sistema de potencia
1. CONCEPTOS BÁSICOS
ante un evento de carácter transitorio, puede ser estudiado por
medio de modelos analógicos o programas digitales. Uno de los
1.1 Sobre el modelado de sistemas de potencia y los
métodos analógicos implementados es el TNA (Analizador de
fenómenos transitorios transitorios en redes), que se trata de una representación del
Como se mencionó en la Introducción, los fenómenos que puedan sistema en escala reducida. Los modelos reflejan la respuesta
alterar la operación normal de un sistema de potencia deben ser eléctrica del sistema real.
tenidos en cuenta en el momento de la elaboración de un diseño
Por otro lado, los métodos digitales, por tratarse de métodos
adecuado de la red. Para lograr este objetivo, existe la posibilidad
analíticos, están basados en las relaciones matemáticas que
de modelar, simular y analizar los sistemas y los fenómenos existen entre tensiones y corrientes, y tienen la capacidad de
transitorios asociados a estos. representar parámetros dependientes de la frecuencia,
Un modelo es una representación del circuito en forma de concentrados y distribuidos, y parámetros lineales y no lineales
resistencias, inductancias y capacitancias, que responden muy como descargadores de sobretensión y saturación de
bien a lo que se quiere representar. Los sistemas de potencia transformadores.
abarcan una gran cantidad de componentes. Un modelo apropiado
Cada método posee ventajas y desventajas, pero en general, los
puede ser construido conectando modelos de grandes elementos a resultados obtenidos sea de métodos analógicos o digitales, son lo
través de cables. suficientemente aproximados como para poder hacer un correcto
Para obtener un modelo que se comporte de forma correcta, es dimensionamiento de los equipos [5].
importante tener los valores precisos de resistencias, inductancias
y capacitancias. Algunos pueden ser obtenidos de las placas de 1.2 El software ATP
equipos, otros a partir de mediciones en campo.
El ATP (Alternative Transient Program) versión del EMTP
El grado de detalle del modelo dependerá del tipo de estudio que (ElectroMagnetic Transient Program), es la herramienta digital
se quiere realizar. Un modelo para un estudio de régimen

14
Arias (2015)
de uso universal para el modelado y simulación de transitorios El programa contiene una gran cantidad de elementos que
electromagnéticos y electromecánicos [6]. facilitan la modelación de la red, como por ejemplo
“El programa permite modelar redes monofásicas o polifásicas Elementos RLC
con cualquier topología, las cuales son traducidas sobre la forma Modelos para la representación de líneas aéreas y cables
de un conjunto de ecuaciones diferenciales y algebraicas que son Resistencias no lineales
resueltas en intervalos discretos de tiempo. La topología de la red
Interruptores estadísticos, sistemáticos y controlados por
es representada por la matriz de admitancia de barras y la
variables (tiempo, tensión)
formulación matemática es basada en las características (método
de Bergeron) para los elementos con parámetros distribuidos y en Fuentes de corriente y de tensión
la regla de integración trapezoidal para los elementos Máquinas sincrónicas
representados por parámetros concentrados.” [5]. Modelos de máquinas universales
Los fines de uso del ATP están orientados al diseño de sistemas de Sistemas de control en general
potencia y a la solución de problemas de operación. Casos típicos Algo para destacar del software es la implementación de TACS,
de estudio son: que permiten simular interacciones dinámicas entre la red
Transitorios de maniobra eléctrica y los sistemas de control de determinados equipos, y el
uso de MODELS, que es un lenguaje de programación que
Transitorios generados por descargas atmosféricas
permite al usuario desarrollar sus propios modelos según la
Coordinación de asilamiento necesidad. [5]
Transitorios mecánicos
HVDC 1.3 Aspectos importantes a considerar en la modelación de
Ferroresonancia redes eléctricas
Análisis de armónicos Con el avance de la ingeniería, la expansión y crecimiento de la
Arranque de motores demanda y la normatividad, los sistemas eléctricos se han ido
Sistemas de control en general desarrollando a tal punto que ya no es posible, o al menos se
Estudio de sistemas desbalanceados dificulta mucho y no es práctico, realizar análisis de forma
analítica, sino que es necesaria la implementación de software

15
Arias (2015)
especializado. La modelación de los sistemas eléctricos juega una modelado de sistemas eléctricos, cada una presenta ventajas y
parte fundamental en el diseño, planeación y operación de los desventajas respecto a otra y se comportan mejor en determinados
sistemas de potencia. Si se parte de modelos incorrectos y poco estudios. La elección del software será de acuerdo a las
precisos, los daños materiales podrían ser incalculables e inclusive necesidades.
muchas vidas podrían verse comprometidas. Aunque es imposible
Información necesaria para la simulación. Los modelos que tiene
diseñar una guía específica para la modelación y simulación de
por defecto cada software, están basados en la solución de algún
fenómenos eléctricos, ya que cada sistema tiene sus
modelo matemático que requiere ciertas variables de entrada,
particularidades, pueden establecerse pautas generales que
condiciones iniciales y/o parámetros geométricos y/o eléctricos,
garanticen la obtención de resultados congruentes con la realidad,
para generar ciertas salidas.
permitiendo hacer análisis correctos.
Para lograr que los resultados obtenidos sean correctos, deben
El proceso de la simulación. Gran parte del éxito de la
tenerse los valores precisos de los elementos del sistema. Si esta
simulación, se encuentra en la identificación exacta de lo que se
información no se conoce, debe ser solicitada.
requiere. Para llevar a cabo esa identificación deben seguirse estas
pautas: Validación del modelo. Es fundamental la verificación del modelo
que se ha construido. La verificación puede hacerse de varias
Delimitar la región de influencia o de interés dentro del
maneras. Una forma de verificar el modelo, es comprobando sus
sistema
resultados a medida que se va construyendo, lo que facilita la
Definir las variables objeto de estudio (tensiones, corrección de los errores localizados, ya que si la red es más
corrientes, cargabilidad, etc) compleja, la solución de errores puede tornarse muy complicada.
Debe hacerse uso del sentido común para evaluar ciertos
Identificar los parámetros que pueden tener influencia
resultados y seleccionar y configurar de forma correcta los
sobre esas variables (acoples, resistividad, velocidad del
elementos. El modelo siempre debe tratar de reproducir la
viento, etc)
realidad, aunque pueden presentarse discrepancias, deben
Tener completamente definidas las dudas que se quieren considerarse límites que no alteren de forma trascendental los
resolver con la simulación resultados. Siempre se debe ser conservativo sin sobrepasarse.
Otro punto crítico dentro del modelado de redes, es la selección Monitoreo de las variables. Deben tenerse claros las variables a
de la herramienta digital. Existen múltiples programas para el medir, los puntos de medida y el medidor seleccionado (corriente

16
Arias (2015)
pico, tensión pico, corriente rms, tensión rms, etc) para lograr una La expresión de valores en p.u resulta bastante útil,
medición correcta. La medición de variables está ligada al paso de además de que facilita la lectura de datos.
cálculo definido para la simulación y para el medidor (si aplica),
ya que si no se asigna un paso de muestreo correcto, puede 1.4 Ondas viajeras
perderse información importante en la medida. El uso de
Cuando se presentan perturbaciones energéticas como descargas
nomenclatura (es decir, nombrar los nodos en los que se hace la
atmosféricas o maniobras de conexión o desconexión de
medición), facilita la lectura de datos. Generar nombres
relacionados con lo que se ejecuta en la simulación para nombrar elementos en una línea de transmisión, se generan ondas de
los casos de estudio, resulta bastante útil al momento de la tensión y corriente que en las líneas aéreas se propagan a una
velocidad aproximada a la de la luz. Cuando estas ondas llegan al
identificación de la información.
extremo de la línea o en su recorrido encuentran discontinuidades
Configuración de la simulación. Uno de los aspectos más por efectos de apertura o desconexión de alguno de los elementos
importantes en la modelación de sistemas eléctricos, es el tiempo del sistema, se da lugar a la aparición de ondas incidentes,
de paso de cálculo. Este debe seleccionarse según el tipo de reflejadas y refractadas que pueden sobreponerse con las ondas
estudio a realizar. Por ejemplo, si es un estudio en estado estable iniciales y producir sobretensiones dañinas. Estas ondas, debido a
como un Flujo de Carga, el paso de cálculo no tiene que ser tan las pérdidas en la línea, van desapareciendo a medida que
pequeño como el de un estudio de fenómenos transitorios, que transcurre el tiempo [8] , [9].
involucra frecuencias mucho más altas. “Deben considerarse el
Una línea de transmisión en operación puede ser analizada
menor tiempo de viaje en elementos que representen el fenómeno
de la propagación y el menor tiempo requerido para la solución considerando un pequeño tramo de su longitud que tiene asociado
óptima de elementos no lineales que existan en la red”. [7] una impedancia z por unidad de longitud 𝑧 = 𝑟𝑥 + 𝑗𝜔𝐿, y una
admitancia y por unidad de longitud. En este sentido, rx y jωL
Tratamiento de los resultados representan la resistencia y la reactancia serie por unidad de
longitud de ese pequeño tramo. La impedancia para ese pequeño
Es adecuada la presentación de tablas y gráficas que
tramo de longitud dx es zdx y la admitancia es ydx [3].
contengan la información de la simulación.
Para obtener el modelo matemático de la onda viajera, se
La estadística es muy implementada para el análisis de
considera una línea de transmisión sin pérdidas (loss-free), es
volúmenes de datos, pero en algunos casos presenta
decir, sólo se considera z= jωL, debido a que el tratamiento del
restricciones.

17
Arias (2015)
modelo al incluir las pérdidas por resistencia se torna más

𝜕𝑉 𝜕𝐼
complejo. = −𝐿 (3)
𝜕𝑥 𝜕𝑡
ΔL 𝜕𝐼 𝜕𝑉
= −𝐶 (4)
𝜕𝑥 𝜕𝑡
I x
V Δc Increasing Se usan derivadas parciales dado que la tensión 𝑉 y la corriente 𝐼,
son función ambas de la posición y del tiempo.
El signo menos en las ecuaciones surge de la convención usada.
La Figura 1 muestra el incremento de x hacia la derecha. Con la
Δx dirección de la corriente en la dirección indicada, V e I
disminuirán con el incremento de x.
Figura 1. Representación de un pequeño tramo de una línea de
transmisión Por diferenciación, I puede ser eliminada de la ecuación (3) y la
ecuación (4). Se deriva la ecuación (3) con respecto a x y la
De la Figura 1, el voltaje en el inductor está determinado por: ecuación (4) con respecto a t:
𝜕𝐼 𝜕2𝑉 𝜕2𝐼
−𝛥𝑉 = 𝐿 𝛥𝑥 (1) = −𝐿 (5)
𝜕𝑡 𝜕𝑥 2 𝜕𝑡 𝜕𝑥
𝜕𝑉 𝜕2𝐼 𝜕2𝑉
−𝛥𝐼 = 𝐶 𝛥𝑥 (2)
= −𝐶 (6)
𝜕𝑡 𝜕𝑥 2 𝜕𝑡 𝜕𝑥
El tramo de longitud 𝛥𝑥 puede expresarse como un tramo de 𝜕2 𝐼
longitud infinitesimal dx, por lo que la ecuación (1) y la Eliminando :
𝜕𝑡 𝜕𝑥
ecuación (2) quedan

18
Arias (2015)
𝜕𝑉 𝜕𝐼
𝜕2𝑉 𝜕2𝑉 = −𝐿
= 𝐿𝐶 (7) 𝜕𝑥 𝜕𝑡
𝜕𝑥 2 𝜕𝑡 2
= 𝐿ѵ[𝑓 ′1 (𝑥 − ѵ𝑡) − 𝑓 ′ 2 (𝑥 + ѵ𝑡)] (11)
Ahora, resolviendo para I en lugar de V, se obtiene una ecuación
idéntica para la corriente:
Integrando a ambos lados de la ecuación con respecto a x
𝜕2𝐼 𝜕2𝐼 𝑉(𝑥, 𝑡) = 𝐿ѵ[𝑓1 (𝑥 − ѵ𝑡) − 𝑓2 (𝑥 + ѵ𝑡)]
= 𝐿𝐶 (8)
𝜕𝑥 2 𝜕𝑡 2
𝑉(𝑥, 𝑡) = 𝑍0 𝑓1 (𝑥 − ѵ𝑡) − 𝑍0 (𝑥 + ѵ𝑡) (12)
La ecuación (7) y la ecuación (8) constituyen la ecuación de onda
viajera en una línea de transmisión. Considerando la ecuación De la ecuación (10) y la ecuación (12) se puede observar la
para la corriente, D’ Alembert mostró que esta cumple con la directa proporcionalidad que existe entre la corriente y el voltaje,
solución general: siendo ese factor de proporcionalidad la impedancia
característica, 𝑍0 .
−
𝑡 Como solución a la ecuación de onda, la única condición es que
𝐼 = 𝑓 [𝑥 + 1] (9) las funciones sean diferenciables.
(𝐿𝐶)2
Las funciones 𝑓1 (𝑥 − ѵ𝑡) y 𝑓2 (𝑥 + ѵ𝑡) son funciones de distancia
−
1 en tanto que x y ѵ𝑡 tengan dimensiones de longitud. En un
El término (LC) , es la velocidad ѵ de la onda, por esto, la
2
instante de tiempo t = t1, es posible obtener una ubicación espacial
ecuación (9) puede ser escrita: instantánea de la corriente o el voltaje en la línea. Una solución
igualmente válida para la ecuación de onda es
𝐼(𝑥, 𝑡) = 𝑓1 (𝑥 − ѵ𝑡) + 𝑓2 (𝑥 + ѵ𝑡) (10)
𝑥 𝑥
La ecuación (9) es solución de la ecuación (8). La solución para el 𝐼(𝑥, 𝑡) = 𝑓1 ( − 𝑡) + 𝑓2 ( + 𝑡) (13)
ѵ ѵ
voltaje puede ser obtenida de la ecuación (3):

19
Arias (2015)
𝑥 𝑥 𝑉1
𝑉(𝑥, 𝑡) = 𝑍0 𝑓1 ( − 𝑡) − 𝑍0 𝑓2 ( + 𝑡) (14) 𝐼1 = (15)
ѵ ѵ 𝑍𝐴
La ecuación (13) y la ecuación (14) permiten concentrarse en una Las ondas de tensión reflejadas y refractadas son 𝑉2 y 𝑉3
locación particular de la línea y observar lo que pasa como respectivamente. Ahora, las corrientes asociadas a esas ondas
función del tiempo. serán
Reflexión y refracción de ondas viajeras en una línea de
transmisión. La proporcionalidad existente entre las ondas de 𝑉2
𝐼2 = − (16)
tensión y corriente definidas por el parámetro 𝑍0 , sufre 𝑍𝐴
modificaciones cuando las ondas se encuentran con
discontinuidades en la línea. La onda reflejada de tensión y su 𝑉3
𝐼3 = (17)
compañera la corriente, viajan en sentido opuesto al de la onda 𝑍𝐵
incidente. La onda refractada penetra la discontinuidad. Las
amplitudes de las ondas reflejadas y refractadas, al igual que la Se Asignan los subíndices 1, 2 y 3 para las ondas incidente,
proporcionalidad existente entre la tensión y la corriente reflejada y refractada respectivamente, y por definición, la
incidentes, dependen de la impedancia característica del medio. dirección positiva de x es la dirección de la onda incidente. Note
Las ondas viajeras cumplen con las ecuaciones de continuidad en que I2, dado que está viajando en la dirección de menos x, tiene
una línea, es decir, se cumple la ley de conservación de energía signo opuesto a V2. Si el voltaje y la corriente son continuas en
[10]. una frontera, se deduce que
Considere la unión entre dos líneas de impedancias característica
diferentes, ZA y ZB, y suponga que ZA > ZB. Esta podría ser por 𝑉1 = 𝑉2 + 𝑉3 (18)
ejemplo la unión entre una línea aérea y un cable. Suponga que
una onda de tensión con una amplitud V1 viaja desde la fuente a 𝐼1 = 𝐼2 + 𝐼3 (19)
través de la línea aérea acercándose a la unión con el cable. La
corriente tendrá la misma forma y una amplitud Sustituyendo la ecuación (15), ecuación (16) y ecuación (17) en la
ecuación (19)

20
Arias (2015)
ZA≠ZB
𝑉1 𝑉2 𝑉3
− = (20)
𝑍𝐴 𝑍𝐴 𝑍𝐵
De la ecuación (18) y la ecuación (20) es posible escribir la onda

reflejada en términos de la onda incidente: V2
𝑍𝐵 − 𝑍𝐴
𝑉2 = ( )𝑉 (21)
𝑍𝐵 + 𝑍𝐴 1
I3
I1
ZB −ZA
El término ( ) es llamado coeficiente de reflexión y puede I2 V3
ZB +ZA V1
ser positivo o negativo dependiendo de los valores de ZA y ZB
Al sustituir la ecuación (21) en la ecuación (20) se obtiene la onda Figura 2. Reflexión y Refracción de la onda incidente (a)
refractada en términos de la onda incidente:
ZA≠ZB
2𝑍𝐵
𝑉3 = ( )𝑉 (22)
𝑍𝐵 + 𝑍𝐴 1
V1
I2 V3
2ZB
El término ( ) es llamado coeficiente de refracción y puede I1
I3
ZB +ZA
variar entre 0 y 2 dependiendo de los valores de ZA y ZB. [8]
V2
En la Figura 2, Figura 3, Figura 4 y Figura 5 se muestran varios
casos de reflexión y refracción de ondas incidentes de tensión y
corriente al encontrarse con una discontinuidad [10].
Figura 3. Reflexión y Refracción de la onda incidente (b)

21
Arias (2015)
ZA≠ZB 1.5 Maniobras ejecutadas en Sistemas Eléctricos de Potencia.

Energización de líneas de transmisión. Al energizar una línea de
transmisión en vacío, la sobretensión se origina por la
V1
V3
discrepancia de polos en el cierre del interruptor, lo que se traduce
V2
I1 en la no simultaneidad del cierre en sus contactos. Así, luego de
I3
haber cerrado la primera fase, se generan ondas de tensión en las
otras dos fases producto de su acoplamiento. Estas ondas se
I2
propagan a lo largo de la línea hasta alcanzar su otro extremo, en
donde al chocar con la impedancia de circuito abierto, se reflejan
para superponerse con las ondas que continúan propagándose,
produciéndose así la sobretensión [3], [11].
Recierre de líneas. El recierre de una línea de transmisión tiene
Figura 4. Reflexión y Refracción de la onda incidente (c) por objeto despejar fallas transitorias, y por lo tanto involucra los
procesos de apertura y cierre de los interruptores de la línea.
ZA≠ZB Considerando el caso que el interruptor tenga que desconectar una
línea en vacío, debido a la naturaleza capacitiva del circuito al
momento de interrumpir la corriente por su paso por cero, la
tensión estará pasando por su valor máximo, lo cual provoca que
I2 las tres fases queden con una tensión aproximadamente igual al
valor pico de la tensión fase. Por consiguiente, como producto de
esta maniobra, se genera una “carga atrapada”, la cual
permanecerá en la línea por varios periodos de tiempo, en el orden
I3
I1 de minutos, a no ser que sea drenada por medio de reactores en
V2 V3 los extremos de la línea.
V1
Energización de transformadores. Es un caso similar a la
energización de una línea de transmisión, sólo que el fenómeno se
Figura 5. Reflexión y Refracción de la onda incidente (d)

22
Arias (2015)
vuelve más complejo debido a las características no lineales de su 2. METODOLOGÍA PARA LA SIMULACIÓN DE
impedancia y la presencia de armónicos [11]. TRANSITORIOS DE MANIOBRA CON AYUDA DEL
ATP/EMTP Y SU INTERFAZ GRÁFICA ATPDraw
Inicio y despeje de falla. La ocurrencia de una falla fase-tierra en
algún punto del sistema, origina una sobretensión en las fases
2.1 Información básica del sistema y las maniobras
sanas (fases no falladas) que depende del grado de aterrizamiento
del sistema. En sistemas flotantes, las sobretensiones fase-tierra Para la representación y simulación de las maniobras objeto de
pueden superar el valor de la tensión fase-fase del sistema [12]. estudio se utiliza el programa digital ATP y su interfaz gráfica
ATPDraw en su versión 5.9p3 [13], [14]. La red modelada para
Rechazo de carga. Es una pérdida súbita de carga que puede
realizar los estudios se ilustra en la Figura 6.
generarse debido a la incorrecta actuación de un interruptor.
Debido a que el flujo de corriente por la línea se reduce, la caída El sistema modelado corresponde a una red 500 kV/220 kV que
de tensión debida a su impedancia también lo hace. Esto sumado contiene dos generadores y equivalentes de red en las
al efecto capacitivo de la línea, genera la sobretensión [12]. subestaciones frontera. Las líneas a 500 kV tienen una
compensación reactiva del 80%. Los descargadores de
Conmutación de corrientes inductivas y capacitivas. Dado que la
sobretensión son seleccionados de un catálogo y la característica
corriente y la tensión en cargas inductivas y capacitivas se
de saturación del transformador es extraída de los reportes de
encuentran desfasadas alrededor de 90°, cuando se interrumpe la
prueba de un transformador con condiciones similares al de la red.
corriente en su cruce por cero, la tensión estará en su valor
Las simulaciones son ejecutadas sobre una de las líneas a 500 kV,
máximo, lo cual ocasiona cambios en la tensión entre bornes del
otra a 220 kV y el único transformador 500 kV/220 kV existente
interruptor luego de haber sido interrumpida la corriente. La
en el sistema. En los Anexos se encuentran más detalladamente
tensión que aparece en los bornes del interruptor se denomina
las características de los elementos seleccionados para ejecutar
tensión de recuperación. Cuando la tensión de recuperación tiene
cada una de las simulaciones.
valores muy altos o crece de forma muy acelerada, se puede
ocasionar un reencendido del arco [11].

23
Arias (2015)
Eq. K Las maniobras generadoras de transitorios simuladas son:

Subestación K
500 kV
Energización de líneas
Energización de línea a 500 kV con resistencias de
LCC
LCC
Subestación A
preinserción
500 kV
Energización de transformador 220 kV/500 kV

Energización de transformador con mando sincronizado
Recierre monofásico de líneas
LCC
Eq. B
Subestación B
500 kV
Análisis del arco secundario en línea a 500 kV.
Verificación de TRV en interruptores de línea
Falla terminal
Falla kilométrica
Y
Transformador 220 kV/500 kV

seleccionado para maniobras
Los resultados de corriente y tensión para cada una de las
Y
BCT
simulaciones se expresan en p.u. En la Tabla 1 y Tabla 2 se

Subestación C Gen. E
muestran las corrientes y tensiones base para cada nivel de
220 kV
Subestación E
220 kV tensión.
Tabla 1. Corriente base de transformador para expresar valores
LCC
LCC
Eq. E
LCC
Subestación G
220 kV
Subestación I
220 kV
Eq. D
en p.u.
Subestación D
220 kV
Potencia nominal
Y
Tensión devanado Corriente

BCT
BCT
2 X 64 Mvar
Subestación H Subestación J nominal

18 kV 18 kV
(kV) (kA)
LCC
LCC
Gen. H Gen. J
(MVA)
Subestación F
220 kV
500 1,22
750
Eq. F 220 2,78
Figura 6. Red implementada en el ATP.

24
Arias (2015)
Tabla 2. Tensiones base del sistema para expresar valores en p.u. Los modelos y tipos utilizados en el ATP para representar cada
Tensión Tensión uno de los elementos del sistema se muestran en la Tabla 3 [14],
Tensión nominal [15], [16].
fase – tierra fase – fase
(kV)
(kV) (kV) Según el rango de frecuencias definido para este tipo de
500 408,20 707,10 transitorios [3], [2], debe configurarse el ATP de modo que el
220 179,60 311,10 paso de cálculo no sea superior a los 100 µs, ya que estos
fenómenos pueden alcanzar frecuencias de hasta 10 kHz. Sin
embargo, debido a que el sistema presenta algunas líneas con
Tabla 3. Modelos del ATP utilizados para representar los longitudes inferiores a los 2 km, es necesario ejecutar las
elementos de la red simulaciones con un paso mucho más pequeño para que no se
presenten problemas numéricos. Los ajustes realizados en el ATP
Elemento Modelo en ATP se muestran en la Figura 7.
Línea de Transmisión LCC Template, Bergeron Los elementos de la red seleccionados para la realización de las
maniobras están en la Tabla 4.
Source Type 14 más RL
Equivalente de red Coupled Line sequence
Descargador de sobretensiones
Resistencia no lineal Type 92
(DPS)
Transformador BCTRAN
Saturación Inductor no lineal Type 98
Interruptor Time 3-ph/statistic
Compensación Branch 3ph RLC-Y
Resistor Type 91. TACS
Resistencia de Arco Secundario
controlado por tiempo

25
Arias (2015)
Tabla 4. Elementos seleccionados para la ejecución de las

maniobras
Maniobra Elemento
-Línea D – F 220 kV (izquierda)

Energización -Línea A – B 500 kV
-Transformador 220 kV/500 kV

Recierre monofásico
-Línea A – B 500 kV
Análisis de arco secundario -Línea A – B 500 kV

TRV
-Línea A – B 500 kV
Figura 7. Ajustes del ATP para las simulaciones

26
Arias (2015)
2.2 Modelos implementados para la elaboración de la red

En este numeral se presentan los elementos del EMTP/ATPDraw
implementados para la modelación de la red sobre la cual se hacen
las simulaciones, con el objetivo de informar al lector el porqué
de los modelos seleccionados para el estudio de fenómenos
transitorios de frente lento, cuya principal fuente generadora son
las maniobras de apertura y cierre de interruptores. Algunos
elementos son definidos en este numeral, otros a medida que se
avanza en el documento, como por ejemplo los interruptores.
Modelo de línea. A pesar de que el ATP posee variedad de
modelos que permiten representar las líneas de transmisión aéreas,
no todos ellos son adecuados para el estudio de fenómenos
transitorios. En el numeral 1.4 se habla del fenómeno de onda
viajera en una línea de transmisión, y cómo la refracción de estas
ondas produce las sobretensiones. El modelo de la onda viajera es
uno de los más implementados para la modelación de transitorios
[17], [18]. Uno de los modelos contenidos en el ATP basado en
las ecuaciones para la onda viajera, es el modelo Bergeron, por lo
que resulta útil su implementación para estudio de maniobras en
sistemas de potencia. Se debe tener cuidado si se realizan estudios
donde se involucren frecuencias muy altas como por ejemplo
transitorios de frente rápido, ya que el modelo puede presentar
problemas.
Figura 8. Localización del modelo de línea

27
Arias (2015)
En la Figura 9 se ilustra la configuración para el modelo de línea.

Se pueden elegir opciones como considerar si es línea transpuesta,
efecto skin, etc. Debe conocerse la resistividad del terreno o de lo
contrario se deben usar valores típicos. El parámetro Freq init
corresponde a la frecuencia nominal de la red.
LCC
0
Figura 10. Representación gráfica del modelo de línea antes de
ser configurado
LCC
Figura 9. Ajustes del modelo de línea Figura 11. Representación gráfica del modelo final de la línea

28
Arias (2015)
Modelo de Transformador. El ATP tiene en su librería varios Se aclara que la obtención de la saturación del transformador
modelos de transformadores. Entre los más usados se encuentran: implementado en la red modelada, no se hizo por medio de la
el Saturable, el BCTRAN y el XFMR. El modelo seleccionado rutina SATURA, ya que desde el principio se contó con la
para representar los transformadores de la red es el BCTRAN. característica. La metodología planteada en este documento está
dirigida a personas que han tenido poca relación con el ATP, por
El modelo de transformador BCTRAN es una representación
lo que no se considera apropiada la explicación de la rutina, ya
lineal simplificada de transformadores monofásicos o trifásicos en
que su proceso es complejo y requiere de conocimientos más
forma de matrices de impedancias ó admitancias de rama,
avanzados. Si el lector desea profundizar en el manejo de la rutina
teniendo en cuenta los acoples mutuos.
SATURA, puede tener acceso a esta información consultando el
Para estudios de ferroresonancia, corrientes inrush y otros Theory Book del ATP.
comportamientos no lineales, es necesario considerar los efectos
Los valores de corriente y flujo en la configuración mostrada en la
de la saturación, una limitación que ha sido eliminada gracias al
Figura 18, corresponden a la corriente nominal y flujo nominal del
uso de ramas externas formadas por elementos no lineales
transformador, que se obtiene a partir de la expresión
conectadas al modelo BCTRAN. Dichos elementos pueden ser del
tipo 93 ó 98 (curva de saturación, implementada en este caso), ó
𝑉𝑝−𝐿𝑁
tipo 96 (curva de saturación que toma en cuenta el lazo de ∅𝑛 = (23)
hsitéresis). La rama que representa la saturación debe ser 𝜔
conectada en el devanado cuya corriente nominal fue la base para donde:
el cálculo de la característica. [15] ∅𝑛 : Flujo nominal
Cuando no se dispone de la característica de saturación del 𝑉𝑝−𝐿𝑁 : Tensión pico línea – neutro del lado del transformador en
transformador, debe hacerse uso de la rutina SATURA. El gran el cual se conecta la rama de saturación
problema que presenta la saturación, es que no es posible obtener 𝜔: Frecuencia nominal del sistema
la característica en la región saturada, es decir, después del codo El flujo obtenido a partir de la ecuación (23) debe coincidir con
de saturación, lo que introduce errores en la simulación de uno de los puntos de la característica introducida en el inductor, y
transitorios. En el numeral 2.4, se habla de cómo corregir este por ende, debe tener asociado una corriente, que es la corriente
problema. nominal.

29
Arias (2015)
BCT
Y Y
Figura 12. Representación gráfica del modelo BCTRAN
Para evitar que se presenten problemas numéricos en la

simulación debido a la delta del devanado terciario del
transformador, se conecta en este una capacitancia de 3 nF que
hace las veces de “tierra ficticia” y evita que la delta no quede
referenciada.
Figura 13. Localización del modelo BCTRAN

30
Arias (2015)
Figura 15. Localización del inductor no lineal para representación

Figura 14. Ajustes y configuración del modelo BCTRAN de la saturación

31
Arias (2015)
Figura 16. Representación gráfica del inductor no lineal tipo 98
BCT
Y Y
Figura 18. Ajustes del inductor no lineal tipo 98

Figura 17. Modelo final de transformador más saturación

32
Arias (2015)
Modelo de descargador de sobretensiones. Para el modelo del

descargador de sobretensiones, se utiliza el modelo de resistencia
no lineal Type 92. La ventaja de este modelo, es que a diferencia
de otros modelos no lineales, no requiere de una característica
exponencial exacta, sino que se comporta muy bien ante
características definidas por pequeños tramos de recta [13], que
precisamente son las suministradas en los catálogos de fabricantes
de estos dispositivos. La característica del descargador es extraída
de catálogos de fabricantes. Debe utilizarse la característica V-I
30/60 µs, ya que este frente corresponde a transitorios de
maniobra [19]. La configuración mostrada en la Figura 22 debe
dejarse como vienen por defecto.
PE
R(i)
Figura 19. Característica de saturación ingresada en el inductor Figura 20. Representación de la resistencia no lineal tipo 92
no lineal tipo 98

33
Arias (2015)
Figura 21. Localización resistencia no lineal tipo 92 Figura 22. Ajustes de la resistencia no lineal tipo 92

34
Arias (2015)
R(i)
R(i)
R(i)
Figura 24. Modelo final de descargador de sobretensiones
Equivalente de red. El equivalente de red o equivalente de
Thevenin es la manera de representar de forma reducida
determinada porción del sistema que puede que no sea de interés
en el modelo, pero que tiene influencia sobre la zona que pretende
ser estudiada. Los parámetros eléctricos que deben ser tenidos en
cuenta para el modelo del equivalente de red son: R0, X0, R1, X1 y
la tensión. Estos parámetros son los asociados al punto de
conexión [17] y son obtenidos de la reducción de red realizada
previamente en el software DIgSILENT.
Los componentes utilizados para modelar el equivalente de red
son la fuente Type 14 y el Symmetric RL coupled line
Figura 23. Característica de maniobra del descargador ingresada

en la resistencia no lineal tipo 92

35
Arias (2015)
Figura 25. Localización de la fuente Type 14 Figura 26. Ajustes de la fuente Type 14

36
Arias (2015)
Figura 28. Representación gráfica de la fuente Type 14
Figura 29. Representación gráfica del Symmetric RL coupled line
Figura 27. Localización del Symmetric RL coupled line

37
Arias (2015)
Otros elementos de la Red

Compensación capacitiva
Figura 30. Ajustes del Symmetric RL coupled line Figura 31. Localización del modelo de compensación capacitiva

38
Arias (2015)
Compensación de línea
Figura 32. Representación de la compensación capacitiva luego

de su configuración
Figura 33. Localización del modelo de compensación inductiva

39
Arias (2015)
Figura 34. Representación de la compensación inductiva luego de

su configuración
2.3 Energización de líneas

Selección, configuración y asignación de tiempo de los
interruptores. Para la ejecución de las maniobras de energización
es necesaria la implementación de interruptores monopolares
estadísticos. Para encontrar y seleccionar de forma rápida estos
dispositivos en el ATP, se despliega el menú ilustrado en la Figura
35 al dar clic derecho sobre la hoja de trabajo.
Figura 35. Localización de interruptores estadísticos en el menú
rápido de elementos del ATP

40
Arias (2015)
Los interruptores estadísticos se programan para efectuar 100

energizaciones, de modo que se pueda hacer un barrido amplio de
la onda de tensión [1]. Las simulaciones se programan con una
duración mínima de 100 ms para observar cómo se comporta el
sistema luego del transitorio. La configuración del menú
ATP Settings se muestra en la Figura 36 y Figura 37. Los
parámetros encerrados en los recuadros rojos, deben dejarse como
vienen por defecto. Para la implementación del interruptor
estadístico se define uno de los interruptores como “Master” y los
otros dos como “Slave”. La asignación de tiempos a cada uno de
los interruptores se hace en base a la norma IEC-60071-4,
referencia [17]: “The closing times of the three phases of the
breaker have a great influence on the generated overvoltage. Two
parameters are taken into account. First, the closing command
may happen at any time, which means anywhere on the power-
frequency sine. Then the three phases, which are supposed to
receive the command at the same time, respond to this command
with a random delay due to mechanical dispersions and
prestrikes. The strategy the most commonly adopted is as follows:
the closing command is supposed to be uniformly distributed on
the sine and real closing times are supposed to be around it, with
a Gaussian probability. If the transient behaviour of the system is
symmetric in terms of polarities, one can limit the excursion of the
switching order time to under one-third of a power voltage
period. Concerning the Gaussian distributions, the standard
deviation of the inter-phase time distribution is not always known.
A value of 0,8 ms to 2 ms may be used as a typical value. IEC Figura 36. Ajustes del ATP Settings para simulación con
62271-100 allows 5 ms.” interruptores estadísticos (a)

41
Arias (2015)
El parámetro Xopt en la Figura 37, debe tener el valor de la

frecuencia nominal del sistema para que los valores de las
inductancias sean tenidos tomados por el ATP en Ohmios y no en
Henrios.
Las variables de importancia en el interruptor estadístico son el
tiempo T, y la desviación estandar dev, que es la discrepancia
asociada al tiempo T. La variable T en el interruptor “Master”, en
este caso la Fase A, es el tiempo de orden de cierre del interruptor.
La variable T en los interruptores “slave”, en este caso las fases B
y C, corresponden a la dispersión de polos existente con respecto
al interruptor “Master”. La norma IEC-60071-4 sugiere tiempos
de orden de cierre aleatorios de la onda de tensión y considera
como típica una desviación estándar entre 0,80 ms y 2 ms. De
igual forma, la norma recomienda emplear una Distribución
Normal o Gaussiana.
La dispersión entre polos de un interruptor trifásico al momento
del cierre no puede conocerse, es un parámetro que se define entre
cliente y fabricante [20], pero puede estar alrededor de los 5 ms
máximo [21]. Para efectos de la simulación, se asume una
dispersión de polos de 5,55 ms y 2,77 ms para las fases B y C
respectivamente, con respecto a la Fase A. El tiempo de orden de
cierre del interruptor es 33,33 ms y la desviación estándar se
asume de 2 ms. En la Figura 38, Figura 39 y Figura 40 se muesta
la configuración de los interruptores estadísticos según los
criterios mencionados.
Figura 37. Ajustes del ATP Settings para simulación con
interruptores estadísticos (b)

42
Arias (2015)
Figura 38. Configuración del interruptor “Master” (Fase A)

Figura 39. Configuración del interruptor “Slave” (Fase B)

43
Arias (2015)
Con la configuración de tiempos del interruptor ya establecida, se

procede a realizar la conexión eléctrica de cada uno de los
interruptores para formar el bloque ilustrado en la Figura 43.
Luego de tener el interruptor estadístico completamente armado,
se conecta su borne de salida al modelo de línea de transmisión
para formar el modelo de línea más interruptor ilustrado en la
Figura 44. El elemento auxiliar implementado para realizar esta
conexión se denomina Splitter trifásico.
Figura 41. Ubicación del Splitter en el menú rápido de elementos

Figura 40. Configuración del interruptor “Slave” (Fase C) del ATP

44
Arias (2015)
0.0L 0.0
PE R(i)
PE R(i)
PE R(i)
S
STAT V
LCC
STAT
STAT
Figura 42. Splitter trifásico Figura 44. Modelo de línea más interruptor estadístico para
simulación de energización de línea
STAT Resultados esperados. Los valores de sobretensiones esperados al
ejecutar la simulación de energización de líneas, se presentan en
la Tabla 5. Estos valores son la base para determinar qué tan
severas o no son las sobretensiones obtenidas [22].
STAT
Tabla 5. Valores esperados por energización de línea
Maniobra Lo esperado
STAT
Sobretensiones
Fase – Tierra de 2 p.u a
Energización de líneas
3 p.u en el extremo en
Figura 43. Conexión del interruptor estadístico para efectos de la vacío.
simulación

45
Arias (2015)
Tratamiento de resultados. Luego de que el ATP finaliza la El archivo .lis, contiene los resultados de las variables medidas en
simulación, se genera un archivo plano con extensión .lis, que la simulación para cada una de las energizaciones realizadas por
contiene toda la información utilizada por el programa para la el interruptor. Estos resultados deben ser tabulados, llevados a
ejecución de las simulaciones y el informe de resultados producto valores en p.u y organizados de forma que se logre identificar la
de la energización. Para abrir el archivo basta con presionar la fase en la cual se presenta el valor más alto de tensión. Luego de
tecla F5 al finalizar el proceso de corrida. identificar la fase con la sobretensión más alta, debe hacerse un
tratamiento estadístico que permita obtener la media y la
desviación estándar (σ) de todas las sobretensiones registradas
asociadas a dicha fase. Con esta información se calcula el valor
estadístico, definido de forma que corresponda a las
sobretensiones con el 95% de probabilidad de ocurrencia [23].
𝑉𝑒𝑠𝑡 = 𝑉𝑚𝑒𝑑 + 1,64 ∗ 𝜎 (24)
donde:
𝑉𝑒𝑠𝑡 : Tensión estadística
𝑉𝑚𝑒𝑑 : Tensión media
𝜎: Desviación estándar
1,64: Coeficiente de confianza para una probabilidad de
ocurrencia del 95% [24]
Figura 45. Ejemplo de archivo .lis Figura 46. Resultados de energización en el archivo .lis

46
Arias (2015)
En la Figura 46 se ilustra la estructura en el que el archivo .lis significa que para la medición de la energía disipada por los
presenta los resultados de la simulación. Los recuadros rojos descargadores, es necesario registrar la corriente que circula por la
resaltan la información que debe identificarse. El primer recuadro resistencia no lineal que representa su modelo.
de arriba hacia abajo, encierra las variables que fueron
monitoreadas, y sus nombres corresponden a los nombres de los
nodos en los cuales se conectaron instrumentos de medida. Lo
delimitado por el óvalo dentro del recuadro, son las variables de
interés en este caso, ya que corresponden a las tensiones
Fase – Tierra. El segundo recuadro, resalta el número de la
energización, el tercero los tiempos en los que se realizó dicha
energización y el cuarto los resultados obtenidos. El .lis muestra
los resultados a modo de matriz, lo que implica que la lectura de
los datos debe hacerse como se indica en la Figura 47.
Figura 47. Correspondencia de resultados en el archivo .lis

La medición de la energía disipada por los descargadores de
sobretensión se hace solamente en la simulación en la que se
emplean los tiempos con los que se presenta la sobretensión más
severa, de tal manera que se obtenga la máxima energía disipada.
El ATP no muestra valores de energía en J como tal, sino un
equivalente de estos en A, es decir, 1 J equivale a 1 A, lo que Figura 48. Configuración del resistor no lineal correspondiente al
modelo de descargador de sobretensiones para registrar corriente

47
Arias (2015)
Simulación del caso crítico. Para llevar a cabo la simulación del

caso más crítico, se identifican los tiempos de cierre de cada fase
del interruptor con los cuales se obtuvo la sobretensión más alta
(siguiendo el procedimiento mencionado anteriormente) y se
realiza nuevamente la simulación pero empleando un interruptor
controlado por tiempo, programado con los tiempos ya
identificados.
0.0L 0.0
PE R(i)
PE R(i)
PE R(i)
S
V
LCC
Figura 50. Selección del modelo de interruptor controlado por

tiempo
Para los casos de energización de líneas simulados, se obtiene que
Figura 49. Modelo de línea más interruptor controlado por tiempo la mayor sobretensión en la línea a 500 kV se genera en la Fase C
para simulación del caso crítico con los tiempos de la energización número 30, y para la línea a
220 kV en la Fase A con los tiempos de la energización
número 41. En la Figura 51 y Figura 52 se muestra la localización
de estos tiempos en el archivo .lis.

48
Arias (2015)
número. En el caso de la Figura 52, los números son 31, 80 y 81

en ese orden respectivo. Estos números corresponden a la
identificación de cada uno de los polos o fases del interruptor. Es
muy importante considerar esta anotación, ya que si el interruptor
no se programa adecuadamente, el resultado será diferente al
esperado. Esta información también está contenida en el archivo
.lis.
Figura 51. Resultados de la energización de la línea 500 kV en la

que se obtuvo la mayor sobretensión
Figura 53. Identificación de los polos del interruptor en el archivo

.lis para asignación correcta de tiempos
En la Figura 53 se ilustra la identificación de polos del interruptor
correspondiente al caso de la Figura 52. Se observa que la Fase A
se identifica con el número 31, la Fase B con el número 80 y la
Fase C con el número 81.
Figura 52. Resultados de la energización de la línea 220 kV en la
que se obtuvo la mayor sobretensión En las opciones de configuración del interruptor controlado por
tiempo, los número 1, 2 y 3 corresponden a las fases A, B y C
Con los respectivos tiempos localizados en la Figura 51 y la respectivamente.
Figura 52, se realiza la programación de los interruptores
controlados por tiempo. Se debe notar que los tiempos en los que
se realiza la energización, se encuentran precedidos por un

49
Arias (2015)
Figura 54. Configuración del interruptor con los tiempos críticos Figura 55. Configuración del interruptor con los tiempos críticos
para energización de línea a 500 kV para energización de línea a 220 kV

50
Arias (2015)
Resultados de la simulación y análisis. En la Tabla 6 se presentan

los valores máximos registrados y los valores estadísticos de la
energización de líneas. En ambos casos de energización, la
disipación de energía en los descargadores de sobretensión es muy
baja, está entre 0% y 5% (Ver Anexos, energía base
descargadores de sobretensión). La mayor sobretensión
Fase – Tierra se presenta en la Fase C en 500 kV, con un valor de
2,16 p.u y un estadístico de 2,23 p.u. En 220 kV la mayor
sobretensión se presenta en la Fase A, con un valor 2,08 p.u y uno
estadístico de 2,30 p.u.
Las sobretensiones por energización se encuentran entre los
rangos definidos en la Tabla 5, lo que verifica la validez de los
modelos implementados y permite afirmar que las sobretensiones Figura 56. Tensión durante energización de línea a 500 kV
no son muy críticas.
Tabla 6. Resultados de energización de líneas
Tensión máxima
Tensión Energía en
Caso de (p.u)
estadística Descargadores
estudio Fase Fase Fase
(p.u) (J)
A B C
Energización
2,01 2,06 2,16 2,23 69807
línea 500 kV
Energización
2,08 2,06 2,07 2,30 5957,40
línea 220 kV Figura 57. Energía en descargadores de sobretensión de línea a
500 kV

51
Arias (2015)
Mitigación del transitorio. Después de tener los resultados más

críticos de energización de líneas, se procede a realizar la
energización con resistencias de preinserción para mitigar el
efecto de la maniobra. La energización con resistencias de
preinserción sólo se hace en la línea a 500 kV, ya que en ésta se
obtuvo el valor más alto de tensión (2,16 p.u en la Fase C). Se
emplea un valor de resistencia de preinserción y tiempo óptimo de
permanencia típicos para una línea a 500 kV. Estos valores son
320 Ω y 8 ms [22], [25]. En la Tabla 7 se encuentran los valores
esperados por la energización con resistencias de preinserción.
Tabla 7. Valores esperados por energización de línea con
resistencias de preinserción
Figura 58. Tensión durante energización de línea a 220 kV
Maniobra Lo esperado
Sobretensiones Fase – Tierra
Energización de líneas con
de 1,30 p.u a 2,20 p.u en el
extremo en vacío.
En la simulación, las tres fases del interruptor 1 cierran en los

tiempos con los que se obtuvo la sobretensión más alta, mostrados
en la Figura 51 y Figura 52. 8 ms después (tiempo de permanencia
de la resistencia), este se abre y al mismo tiempo el interruptor 2
cierra sus contactos, de tal forma que la resistencia de
preinserción queda desconectada del circuito y el flujo de energía
se da a través del interruptor 2.
Figura 59. Energía en descargadores de sobretensión de línea a
220 kV

52
Arias (2015)
R(i)
R(i)
R(i)
Int 2
Int 1 LCC V
Rpre
Figura 60. Modelo para energización con resistencias de

preinserción
Los tiempos de cierre de cada fase del interruptor 2 son
𝑡2𝑖 = 𝑡1𝑖 + 8 𝑚𝑠, (25)
donde:
t1i : Tiempo de cierre de la fase i del interruptor 1
t 2i : Tiempo de cierre de la fase i del interruptor 2
8 ms: Tiempo de permanencia de la resistencia de preinserción
Figura 61. Configuración del interruptor 1 para energización con


53
Arias (2015)
Figura 63. Selección del modelo de resistencia de preinserción
Figura 62. Configuración del interruptor 2 para energización con

resistencias de preinserción Figura 64. Resistencia de preinserción

54
Arias (2015)
En la Tabla 8 se presentan los resultados de energización de la

línea a 500 kV con resistencias de preinserción.
Tabla 8. Resultados de energización de línea con resistencias de
preinserción
Tensión máxima (p.u) Energía en

Caso de estudio Fase Fase Fase descargadores
A B C (J)
Energización
línea a 500 kV
1,30 1,26 1,15 0,60
con resistencias
de preinserción
Se observa una notoria disminución en la magnitud de las

sobretensiones con respecto al caso sin resistencias de
preinserción. La energía en descargadores es aproximadamente
0%, lo que indica que muy probablemente su disipación de
energía se debe a fugas y no al efecto de la maniobra. Las
tensiones se encuentran dentro del rango definido en la Tabla 7.
Se verifica la efectividad de la medida.
Figura 65. Configuración del modelo de resistencia de

preinserción para energización

55
Arias (2015)
600 2.4 Energización de transformador 220 kV/500 kV

[kV]
400 Se realizan energizaciones desde el lado de baja (220 kV) y el
lado de alta (500 kV) del transformador con el fin de determinar
200
la magnitud de la corriente Inrush y las sobretensiones en el
0
extremo que se encuentra en vacío.
Los parámetros que influencian la corriente Inrush son [26]:
-200
Valor instantáneo de la tensión aplicada al transformador

-400
en el instante de energización
-600
0.00 0.02 0.04 0.06 0.08 [s] 0.10 Magnitud y sentido del flujo residual en el núcleo
(file ENER_LINEA_500kV_DET_Rpre.pl4; x-var t) v:X0001A v:X0001B v:X0001C
magnético
Figura 66. Tensión durante energización de línea a 500 kV con
resistencias de preinserción Resistencia e inductancia equivalente serie del circuito
alimentador
0.7
[A] Resistencia e inductancia de dispersión del bobinado
0.6
primario del transformador
0.5
Características magnéticas y geométricas del núcleo
0.4
Valor de la resistencia de preinserción del disyuntor
0.3
Impedancia de carga conectada al secundario
0.2
Velocidad de cierre de los contactos del disyuntor
0.1
Existencia de bobinado terciario conectado en delta en

0.0
0.00 0.02 0.04 0.06 0.08 [s] 0.10 transformadores trifásicos.
(file ENER_LINEA_500kV_DET_Rpre.pl4; x-var t) c:X0001C- c:X0001B- c:X0001A-
Figura 67. Energía en descargadores de sobretensión de línea a Para obtener las corrientes de Inrush se implementa la
500 kV característica de saturación del transformador de forma externa.

56
Arias (2015)
No es necesaria la curva de histéresis, ya que las características Resultados esperados. Los valores de sobretensiones esperados al
fundamentales de la corriente de energización, son obtenidas de ejecutar la simulación de energización de transformadores se
modelos sin histéresis [27]. presentan en la Tabla 9. [22], [28], [29].
Los reportes de prueba de fabricantes de transformadores han
Tabla 9. Valores esperados por energización de transformador
mostrado que no es posible obtener la característica de saturación
en la región saturada, es decir, después del codo de saturación, ya Maniobra Lo esperado
que si las pruebas en vacío sobrepasan tensiones de 1,10 p.u, es
posible que se generen daños en el equipo [15]. Cuando no se Sobretensiones
dispone de la característica de saturación completa, su región Fase - Tierra de 2 p.u a
saturada puede ser modelada extrapolando los dos últimos puntos 3 p.u en el extremo en
Energización de
de la característica con una pendiente del 30% [17]. vacío y corrientes
Transformadores
Inrush entre 2 p.u y
4 p.u para trafos de
interruptores. Para la ejecución de las maniobras de energización
500 MVA y mayores.
de transformadores se implementan interruptores estadísticos. La
configuración de los interruptores es la descrita en el numeral 2.3.
Tratamiento de resultados. El tratamiento de resultados se hace de
la misma forma que se describe en el numeral 2.3. Hay que tener
en cuenta que en este caso ya no sólo van a identificarse, tabularse
y analizar estadísticamente los valores máximos de sobretensión
I en el extremo en vacío del transformador, sino que el fenómeno
V
S
STAT 0.0L 0.0

I Y
BCT
Y
de la corriente de Inrush también se incluye en el análisis, ya que
STAT
es una variable de muy alta importancia. El tratamiento de los
I
STAT resultados de la corriente de Inrush se hace exactamente igual que
como se hace con los resultados de las sobretensiones. Se aclara
PE
PE
PE
R(i)
R(i)
R(i)
que no necesariamente cuando se presenta la sobretensión más

alta se genera el valor de corriente de Inrush más grande, por lo
Figura 68. Modelo para Energización de transformador que hay que realizar análisis independientes.

57
Arias (2015)
La ecuación (24), también puede expresarse en términos de la

corriente, de tal forma que se calcule el valor estadístico, definido
de forma que corresponda a las corrientes de Inrush con el 95% de
probabilidad de ocurrencia
𝐼𝑒𝑠𝑡 = 𝐼𝑚𝑒𝑑 + 1,64 ∗ 𝜎, (26)

donde:
𝐼𝑒𝑠𝑡 : Corriente estadística
𝐼𝑚𝑒𝑑 : Corriente media
𝜎: Desviación estándar
1,64: coeficiente de confianza para una probabilidad de ocurrencia
del 95% [24]
Figura 69. Configuración de las fases del interruptor estadístico

para monitorear corriente
Simulación de los casos más críticos. Para llevar a cabo la
simulación de los casos más críticos en cuanto a tensión y
corriente de Inrush, se identifican los tiempos de cierre de cada
fase del interruptor con los cuales se obtuvo la sobretensión y la
corriente más alta, y se realizan las simulaciones empleando un

58
Arias (2015)
interruptor controlado por tiempo, programado con los tiempos

identificados.
V
S
0.0L 0.0
BCT
Y Y
Figura 71. Resultados de la energización del transformador desde

el lado de alta en la que se obtuvo la mayor sobretensión
PE
PE
PE
R(i)
R(i)
R(i)
Figura 70. Modelo transformador más interruptor controlado por
tiempo para simulación de casos críticos
Para los casos de energización simulados, se obtiene que la mayor
sobretensión en el extremo en vacío del transformador,
energizando desde el lado de alta tensión (500 kV), se genera en
la Fase B con los tiempos de la energización número 25. Para el
el lado de baja en la que se obtuvo la mayor sobretensión
extremo en vació del transformador energizando desde el lado de
baja tensión (220 kV), la sobretensión más alta se genera Con los respectivos tiempos localizados en la Figura 71 y
también la Fase B pero con los tiempos de la energización Figura 72, se realiza la programación de los interruptores
número 20. En la Figura 71 y Figura 72 se muestra la localización controlados por tiempo.
de estos tiempos en el archivo .lis.

59
Arias (2015)
Figura 73. Configuración del interruptor con los tiempos críticos Figura 74. Configuración del interruptor con los tiempos críticos
para energización del transformador desde el lado de alta para energización del transformador desde el lado de baja

60
Arias (2015)
Por otro lado, se obtiene que la mayor corriente de Inrush,

energizando desde el lado de alta tensión (500 kV), se genera en
la Fase B con los tiempos de la energización número 47, y
energizando desde el lado de baja tensión (220 kV), se genera en
la Fase C pero con los tiempos de la energización número 18. En
la Figura 75 y Figura 76 se muestra la localización de estos
tiempos en el archivo .lis.

el lado de alta en la que se obtuvo la mayor corriente Inrush

Figura 77. Configuración del interruptor con los tiempos críticos
el lado de baja en la que se obtuvo la mayor corriente Inrush
para energización del transformador desde el lado de alta

61
Arias (2015)
Resultados de la simulación y análisis. En la Tabla 10 y Tabla 11

se presentan las sobretensiones y corrientes obtenidas por la
energización del transformador desde el lado de alta y baja.
En los resultados expuestos se observa que las sobretensiones no
son considerablemente elevadas. La sobretensión más alta
registrada es de 1,89 p.u y la estadística es de 1,56 p.u, obtenida
en la Fase B al energizar desde el lado de baja tensión. La Energía
disipada por los descargadores de sobretensión se encuentra entre
el 0% y 5% de su capacidad, lo que indica que probablemente esta
disipación se debe a fugas.
El valor más alto de corriente Inrush se presenta en la Fase C al
energizar el transformador desde el lado de alta y corresponde al
valor de 2,87 p.u y un estadístico de 1,46 p.u. Este valor de
corriente Inrush, aunque no es crítico, puede generar problemas
con las protecciones del transformador al momento de ser
energizado, por lo que se debe ser cuidadoso con los resultados.
En cuanto a la energización de transformadores bajo carga, no se
realizan simulaciones en esta condición, ya que usualmente el
transformador opera en su región saturada durante maniobras de
energización con el secundario en vacío [22].
Figura 78. Configuración del interruptor con los tiempos críticos

para energización del transformador desde el lado de baja

62
Arias (2015)
Figura 79. Tensión del transformador en el lado de alta al Figura 80. Energía en descargadores de sobretensión lado de alta
energizar desde el lado de baja

63
Arias (2015)
Tabla 10. Resultados de tensión en la energización del transformador

Tensión máxima (p.u) Tensión Energía en
Caso de estudio Lado energizado
Fase A Fase B Fase C estadística (p.u) descargadores (J)
Energización de Alta (500 kV) 1,26 1,51 1,01 1,46 0,26
transformador
220/500 kV Baja (220 kV) 1,44 1,89 1,68 1,56 4,72
300
[kV]
200
100
-100
-200
0.03 0.04 0.05 0.06 0.07 0.08 0.09 [s] 0.10
(file ENER_TRF_H_DET_V.pl4; x-var t) v:SECA v:SECB v:SECC
Figura 81. Tensión en el lado de baja del transformador al energizar desde el lado de alta

64
Arias (2015)
0.30 3000
[A] [A]
2000
0.25
1000
0.20
0
0.15
-1000
0.10
-2000
0.05 -3000
0.00 -4000
0.00 0.02 0.04 0.06 0.08 [s] 0.10 0.00 0.02 0.04 0.06 0.08 [s] 0.10
(file ENER_TRF_H_DET_V.pl4; x-var t) c:SECA - c:SECB - c:SECC - (file ENER_TRF_H_DET_I.pl4; x-var t) c:SE-BA -PRIMA c:SE-BB -PRIMB c:SE-BC -PRIMC
Figura 82. Energia en descargadores de sobretensión lado de baja Figura 83. Corriente Inrush del transformador al energizar desde
alta
Tabla 11. Resultados de corriente en la energización del
3000
transformador [A]
2000
Corriente Inrush
Corriente 1000
Caso de Lado máxima (p.u)
estadística 0
estudio energizado Fase Fase Fase
(p.u) -1000
A B C
-2000
Energización Alta (500 kV) 2,10 0,58 2,87 1,46
-3000
de
transformador -4000
Baja (220 kV) 1,00 0,29 1,51 1,56
220/500 kV -5000
0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 [s] 0.10
(file ENER_TRF_L_DET_I.pl4; x-var t) c:X0022A-SECA c:X0022B-SECB c:X0022C-SECC
Figura 84. Corriente Inrush del transformador al energizar desde

baja

65
Arias (2015)
Mitigación del transitorio. La energización de equipos con mando

sincronizado es otra de las técnicas implementadas para la
mitigación de transitorios por efectos de energización. El mando
sincronizado consiste en la programación del tiempo de cierre de
los contactos del interruptor de tal forma que este se haga en
puntos específicos de la onda de tensión, logrando así que el
impacto sea menos severo y que el equipo conserve por mucho
más tiempo su vida útil [22], [5].
Se realiza la energización desde el lado de baja (220 kV) y el lado
de alta (500 kV) del transformador ajustando los tiempo de cierre
para mando sincronizado. Para las simulaciones de energización
con mando sincronizado, se emplean los tiempos recomendados Figura 85. Tiempos óptimos para el mando sincronizado, lado
por el COES (Comité de Operación Económica del Sistema 220 kV
Interconectado Nacional, Perú) para una conexión Yn:
Cierra primero la Fase A en el valor máximo de la onda de
tensión, luego cierran simultáneamente las fases B y C ¼ de ciclo
más tarde, el cual coincide con el cruce por cero de la Fase A, en
este instante las fases B y C tienen igual valor pero signo opuesto.
Estos tiempos óptimos de cierre para mando sincronizado desde
500 kV y 220 kV se ilustran en la Figura 85 y Figura 86. Ambas
figuras corresponden a la onda de tensión en el primario y
secundario del transformador en operación normal. Para la
simulación se emplean interruptores controlados por tiempo.
Figura 86. Tiempos óptimos para el mando sincronizado, lado

500 kV

66
Arias (2015)
Figura 87. Configuración del interruptor para energización con Figura 88. Configuración del interruptor para energización con
mando sincronizado desde alta mando sincronizado desde baja

67
Arias (2015)
En la Tabla 12 se presentan los resultados de la energización del 5.0
transformador con mando sincronizado. [A]

3.5
Se observa que al implementar la energización con mando

sincronizado, el valor de la corriente de Inrush se reduce a valores 2.0
menores a 1,00 p.u, lo que hace posible disminuir la probabilidad 0.5

de que se presenten problemas con las protecciones en el
momento de la energización del transformador. También se -1.0
observa que la onda de tensión mejora y no presenta distorsión

-2.5
como en el caso sin mando sincronizado (MS).
-4.0
0.00 0.02 0.04 0.06 0.08 [s] 0.10
(file ENER_TRF_H_MS.pl4; x-var t) c:SE-BA -PRIMA c:SE-BB -PRIMB c:SE-BC -PRIMC
300
Figura 90. Corriente Inrush del transformador al energizar desde
[kV]
200
Alta
100
-100
-200
-300
0.00 0.02 0.04 0.06 0.08 [s] 0.10
(file ENER_TRF_H_MS.pl4; x-var t) v:SECA v:SECB v:SECC
Figura 89. Tensión en el secundario del transformador al

energizar desde el lado de alta
Figura 91. Tensión en el lado de alta del transformador al
energizar desde el lado de baja

68
Arias (2015)
Tabla 12. Resultados de corriente y tensión en la energización del transformador con MS

Tensión máxima (p.u) Corriente Inrush máxima (p.u) Energía en
Lado
Caso de estudio descargadores
energizado Fase A Fase B Fase C Fase A Fase B Fase C
(J)
Energización de Alta (500 kV) 1,33 1,17 1,16 0,004 0,004 0,004 0,35
transformador
220/500 kV con
mando Baja (220 kV) 1,02 1,02 1,02 0,005 0,007 0,005 0,86
sincronizado
20
[A]
15
10
-5
-10
-15
-20
0.00 0.02 0.04 0.06 0.08 [s ] 0.10
(f ile ENER_TRF_L_MS.pl4; x-var t) c:X0022A-SECA c:X0022B-SECB c:X0022C-SECC
Figura 92. Corriente inrush del transformador al energizar desde baja

69
Arias (2015)
0.90 2.5 Recierre monofásico de líneas

[A]
0.75 Selección, configuración y asignación de tiempo de los
interruptores. Para la ejecución de las maniobras de recierre
0.60
monofásico se emplean interruptores estadísticos, pero en esta
0.45 ocasión, solo para la fase encargada de ejecutar el recierre.
0.30 Los parámetros a considerar en la programación de los tiempos de
los interruptores son la carga atrapada, que puede tener una
0.15
duración de hasta 3 minutos [30], y el tiempo de operación de las
0.00
protecciones [22], [3]. No se consideran fallas monofásicas, ya
0.00 0.02 0.04 0.06 0.08 0.10
(file ENER_TRF_L_MS.pl4; x-var t) c:PRIMA - c:PRIMB - c:PRIMC -
[s]
que estas drenan la carga atrapada y serían menos severas las
Figura 93. Energía en descargadores de sobretensión del lado de sobretensiones en esta condición [30]. Se selecciona la Fase A
alta del transformador para el estudio de la maniobra.
0.35
La secuencia para la realización de la maniobra es:
[A]
0.30 Se abren los polos de la Fase A en ambos extremos. El primer
interruptor abre en un tiempo de 30 ms; el segundo interruptor,
0.25
que está ubicado en el extremo remoto, abre 30 ms después del
0.20
primero (a los 63,33 ms) para dejar acumular carga atrapada en
0.15 esta fase. Luego del tiempo muerto, considerado de 500 ms
0.10
(tiempo típico en líneas a 500 kV) [31], se efectúa el recierre
monofásico estadístico (a los 563,33 ms).
0.05
Además del interruptor trifásico controlado por tiempo, se
0.00
0.00 0.02 0.04 0.06 0.08 [s] 0.10 emplean interruptores monofásicos igualmente controlados por
(file ENER_TRF_H_MS.pl4; x-var t) c:SECA - c:SECB - c:SECC -
tiempo, que ayudan a modelar la secuencia de los eventos. La
Figura 94. Energía en descargadores de sobretensión del lado de desviación estándar definida para la fase encargada del recierre
baja del transformador estadístico (Fase A) es 2 ms.

70
Arias (2015)
Interruptor 2
S
V
Int monofásico
LCC
Interruptor 1
STAT
Fase A
Encargada del recierre
monofásico estadístico
R(i)
R(i)
R(i)
Figura 95. Modelo para recierre monofásico de línea
En la Figura 95 se observa el modelo final de línea más
interruptor para el recierre monofásico. La resistencia conectada
entre bornes del interruptor 2 y el arreglo de interruptores
monofásicos, evita que se presenten problemas numéricos durante
la simulación. Su valor es muy pequeño, 0,001 Ω.
Figura 96. Configuración del Interruptor 1

71
Arias (2015)
Figura 97. Configuración del Interruptor 2 Figura 98. Configuración de la Fase A encargada de la maniobra
de recierre estadístico

72
Arias (2015)
Figura 99. Selección de los interruptores monofásicos Figura 100. Configuración de los interruptores monofásicos

73
Arias (2015)
Resultados esperados. En la Tabla 13 se presentan los valores

esperados por la maniobra de recierre.
Interruptor 2
Tabla 13. Valores esperados por maniobra de recierre monofásico
Maniobra Lo esperado S
V
Tensiones Fase – Tierra Int monofásico
de 1,80 p.u a 2,50 p.u LCC

Recierre monofásico
en el extremo de la fase
Interruptor 1
en falla Fase A
Encargada del recierre monofásico
Tratamiento de resultados. El tratamiento de resultados se hace de
la misma forma que se describe en el numeral 2.3. Se debe
R(i)
R(i)
R(i)
identificar la sobretensión más alta en el extremo de la fase que
ejecuta el recierre y calcular su valor estadístico a partir de la
ecuación (23).
Simulación del caso más crítico. Para llevar a cabo la simulación Figura 101. Modelo de línea e interruptor para recierre
del caso más crítico, se identifica el tiempo de cierre del monofásico con tiempo crítico
interruptor estadístico con el cual se obtuvo la sobretensión más
alta (siguiendo el procedimiento mencionado en el numeral 2.3) y Para los casos de recierre monofásico simulados, se obtiene que la
se realiza nuevamente la simulación pero empleando un mayor sobretensión en el extremo de la Fase A (fase del recierre)
interruptor monofásico controlado por tiempo en reemplazo de la en 500 kV, se genera con los tiempos de la energización
fase estadística, programado con el tiempo ya identificado. número 44. Para la línea a 220 kV, la sobretensión más alta se
produce con los tiempos de la energización número 6.
En la Figura 102 y Figura 103 se muestra la localización de estos
tiempos en el archivo .lis.

74
Arias (2015)
Figura 102. Resultados del recierre monofásico en la línea a

500 kV en la que se obtuvo la mayor sobretensión
Figura 103. Resultados del recierre monofásico en la línea a

220 kV en la que se obtuvo la mayor sobretensión
Con los respectivos tiempos localizados en la Figura 102 y Figura 104. Configuración del interruptor monofásico controlado
Figura 103, se realiza la programación del interruptor monofásico por tiempo para recierre monofásico con tiempo crítico en línea a
controlado por tiempo encargado de ejecutar el recierre. 500 kV

75
Arias (2015)
Resultados de la simulación y análisis. En la Tabla 14 se

presentan los resultados del recierre monofásico para ambas
líneas. Se obtienen valores cercanos a los límites establecidos en
la Tabla 13.
Tabla 14. Tensión en el extremo de la línea por recierre
monofásico
Tensión Tensión Energía en
Caso de estudio máxima estadística descargadores
(p.u) (p.u) (kJ)
Recierre
monofásico 2,40 2,39 4710
línea a 500 kV
Recierre
monofásico 1,26 2,37 73,77
línea a 220 kV
Figura 105. Configuración del interruptor monofásico controlado

por tiempo para recierre monofásico con tiempo crítico en línea a
220 kV

76
Arias (2015)
5
[MA]
0
0.0 0.2 0.4 0.6 0.8 [s] 1.0
(file RM_500kV_DET.pl4; x-var t) c:X0001C- c:X0001B- c:X0001A-
Figura 106. Tensión en las tres fases línea, 500 kV Figura 108. Energía en descargadores, línea 500 kV
1.00 400
[MV] [kV]
300
0.62
200
100
0.24
-0.14
-100
-200
-0.52
-300
-0.90 -400
0.0 0.2 0.4 0.6 0.8 [s] 1.0 0.0 0.1 0.2 0.3 0.4 0.5 0.6 [s] 0.7
(file RM_500kV_DET.pl4; x-var t) v:X0001A (file RM_220kV_DET.pl4; x-var t) v:X0004A v:X0004B v:X0004C
Figura 107. Tensión en la fase en recierre, línea 500 kV Figura 109. Tensión en las tres fases, línea 220 kV

77
Arias (2015)
400 80
[kV] [kA]
300 70
200 60
100 50
0 40
-100 30
-200 20
-300 10
-400 0
0.0 0.2 0.4 0.6 0.8 [s] 1.0 0.0 0.2 0.4 0.6 0.8 [s] 1.0
(file RM_220kV_DET.pl4; x-var t) v:X0004A (file RM_220kV_DET.pl4; x-var t) c:X0004C- c:X0004B- c:X0004A-
Figura 110. Tensión en la fase en recierre, línea 220 kV Figura 112. Energía en descargadores, línea 220 kV
Los resultados muestran que se presentan sobretensiones
considerables en la fase con falla, especialmente en la línea a
500 kV, donde se registra una tensión máxima de 2,4 p.u y una
estadística de 2,39 p.u. La disipación de energía en los
descargadores de sobretensión para el caso de la línea a 220 kV,
se encuentra alrededor del 2%, mientras que la energía disipada en
los descargadores de sobretensión de la línea a 500 kV se
encuentra alrededor del 80%, un valor elevado al que es necesario
prestarle atención, ya que está muy cerca del límite del equipo. En
base a esta información, se genera la necesidad de hacer un
análisis un poco más profundo en la línea a 500 kV para ver si es
posible o no implementar el recierre monofásico. En el numeral
2.6 se presenta el análisis de arco secundario para esta línea.
Figura 111. Carga atrapada en la fase en recierre, línea 220 kV

78
Arias (2015)
2.6 Análisis del arco secundario en línea a 500 kV Tabla 15. Comparación de algunos modelos para el estudio del
arco
El fenómeno de la corriente de arco secundario es una
consecuencia directa del recierre monofásico. Cuando se presenta Modelo Se basa en Concepto
un cortocircuito en una línea de transmisión, se origina la Ideal para simular
Propiedades
corriente de arco primario, luego de realizar el despeje de la falla, May’r pequeñas
esta sigue siendo alimentada por las corrientes producidas debido físicas del arco
corrientes
al acople electromagnético entre la fase en falla y las fases sanas. Propiedades Ideal para simular
Esta corriente es denominada corriente de arco secundario [3], Cassie
físicas del arco altas corrientes
[32]. Fácil de
Resistencia
Kizilcay implementar, es el
Algunos de los parámetros que influyen en la duración de la dinámica del arco
más usado
corriente de arco secundario son la corriente de arco primario,
Interruptores
acople entre las fases, tensión del sistema, condiciones Modificaciones al
KEMA especiales, difícil
atmosféricas, componente de cd y la geometría de las torres, modelo de May’r
especialmente en líneas a 500 kV y tensiones superiores [33], implementación
Combina las
[34].
características de
Propiedades de
Modelo para el estudio del arco secundario. Existen varios Habedank los modelos en los
Cassie y May’r
modelos desarrollados para el estudio de este fenómeno. Cada que se basa. Dificil
modelo emplea parámetros y supuestos diferentes. En la Tabla 15 implementación
se muestra la información básica de algunos de los modelos Modelo de May’r Fácil
desarrollados para estudio de corrientes de arco [34]. con Tau y P0 implementación.
Schavemaker como funciones Ideal para
de la potencia de representar la zona
entrada de corriente cero
En la Tabla 15 se observa que el modelo de Habedank es una
combinación de los modelos de Cassie y May’r, por lo que puede
representar mejor el fenómeno en zonas de corrientes altas como

79
Arias (2015)
en las zonas de corrientes cercanas a cero [34]. Por esta razón se

implementa este modelo para el estudio de arco secundario.
Implementación del modelo de Habedank en el lenguaje de
MODELS del ATP. El procedimiento para la creación e
implementación de MODELS, se describe a continuación.
Lo primero es buscar el elemento en el menú rápido de elementos
del ATP. Figura 114. Representación gráfica del MODEL en el ATP
Al dar doble clic sobre este elemento, aparece una ventana con los
ajustes del MODEL. Se selecciona la opción Edit.
Figura 113. Localización del MODEL en el menú rápido del ATP Figura 115. Ventana de ajustes del MODEL
Al seleccionar el elemento, aparece en la hoja de trabajo un Al sleccionar esta opción, aparece una nueva ventana como la que
gráfico como el mostrado en la Figura 114. se muestra en la Figura 116.

80
Arias (2015)
OUTPUT Rarc
DATA Uc {DFLT:9.6}, Po {DFLT:1}, TAUc
{DFLT:8.e-4},TAUm
{DFLT:8.e-4}
VAR g, gc, gm, G1, G2, Rarc
HISTORY gc {DFLT:1}
gm {DFLT:1}
g {DFLT:1}
G1 {DFLT:1}
G2 {DFLT:1}
Rarc {DFLT:1E-8}
INIT
Rarc:=1E-8
ENDINIT
EXEC
Figura 116. Editor de texto del MODEL G1:=(Iarc*Iarc)/(Uc*Uc*gc)
G2:=(Iarc*Iarc)/(Po)
Esta ventana es el editor de texto del MODEL, en la cual debe LAPLACE(gc/G1):=(1|s0)/(1|s0+TAUc|s1)
introducirse el siguiente código, que describe el modelo de LAPLACE(gm/G2):=(1|s0)/(1|s0+TAUm|s1)
Habedank para el estudio de arco secundario [34]. g:=(1/gc)+(1/gm)
MODEL Habedank IF (g<1e-10) THEN
Comment Rarc:=1E38
* Modelo de Habedank para representar el * ENDIF
* comportamiento del arco eléctrico primario y IF (g>1e-10) THEN
secundario * Rarc:=g
* Input: Tensión del sistema U y corriente de arco ENDIF
Iarc * ENDEXEC
* Output: comportamiento de Rarc, resistencia de ENDMODEL
arco *
endcomment
INPUT U, Iarc

81
Arias (2015)
A continuación se da clic en la opción Done para guardar los

cambios.
Figura 118. Cuadro de diálogo emergente luego de salvar los

cambios realizados en el editor de texto
Figura 117. Introducción del código en el editor de texto

Seguidamente, aparece un cuadro de diálogo que pregunta si se
está de acuerdo con las variables definidas en el código. Clic en
yes para finalizar el proceso.
Finalmente se abre una ventana que muestra y permite la
configuración de las variables que se definieron en el MODEL.
Clic en ok para terminar. Figura 119. Ventana de edición de las variables del MODEL

82
Arias (2015)
El resultado gráfico final del MODEL se ilustra en la Figura 120.

Luego de tener listo el modelo en lenguaje de MODELS, se
procede a conectar el elemento a la red. La salida del modelo
(Iarc) controla una resistencia TACS Type 91 que varía con
respecto al tiempo de acuerdo al valor que le da el modelo (esta
resistencia no se configura, sus valores dependen del MODEL).
En la Figura 122 se muestra la conexión del MODEL al sistema.
Figura 120. Representación gráfica final del MODEL
Figura 121. Localización de la resistencia TACS Type 91 en el

menú rápido de elementos

83
Arias (2015)
evento), ubicando un medidor de tensión en ambos extremos de la

línea y haciendo el respectivo análisis de la onda para identificar
Interruptor 1 LCC Interruptor 2
los tiempos de cruce por cero. Para el caso en el que se aplica la
falla en la mitad de la línea, el tiempo de falla puede ser el medido
en cualquiera de los dos extremos de la línea
El tiempo de falla en las terminales del interruptor asociado a la
subestación A es 𝑡𝑓 = (0,005216) s y en las terminales del
interruptor asociado a la subestación B es 𝑡𝑓 = (0,005202) s. El
UI
MODEL
Interruptor auxiliar Habedank

polo del interruptor correspondiente a la fase en falla abre en
𝑡𝑐𝑏1 = (0,1 + 𝑡𝑓) s, luego el polo del interruptor en el extremo
U
RARC remoto abre en 𝑡𝑐𝑏2 = (0,033 + 𝑡𝑐𝑏1) s. Finalmente, el
T
interruptor auxiliar, que modela la extinción de la falla, abre en
𝑡𝑎 = (0,5 + 𝑡𝑐𝑏2) s. Los 0,1 s involucrados en el tiempo 𝑡𝑐𝑏1,
Figura 122. Conexión del modelo de estudio de arco secundario a corresponden al tiempo de detección de la falla.
la línea de transmisión Los criterios para evaluar la extinción del arco secundario son los
Configuración y asignación de tiempo de los interruptores. Para definidos por el CESI (Centro Electrotécnico Experimental
efectos de la simulación, se aplica una falla (representada por el Italiano): 40 A pico ó 28 A rms para la corriente de arco
MODEL y la resistencia no lineal) en la fase en la que se evalúa el secundario y 80 kV para el primer pico de la tensión de
recierre (Fase A). Se realizan 3 simulaciones: una en cada restablecimiento, con una pendiente para esta última de 8 kV/ms.
extremo de la línea y otra en la mitad. El instante en el que se Para este análisis es necesario modificar el paso de cálculo del
aplica la falla corresponde al momento en el que la onda de ATP, ya que con el paso definido para las demás simulaciones, el
tensión asociada a la Fase A cruza por cero, de este modo se modelo de arco secundario presenta problemas numéricos.
obtiene la peor condición, ya que en este instante la corriente de
falla presenta la máxima componente de cd [35]. En la Figura 125, Figura 126, Figura 127, Figura 128, Figura 129
y Figura 130 se muestran los ajustes realizados en ATP para cada
La medición de este tiempo se realiza ejecutando una simulación una de las simulaciones.
en la que el sistema opere normalmente (sin ningún tipo de

84
Arias (2015)
Hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh
hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh
h
Interruptor 1 LCC LCC Interruptor 2

Interruptor 1 LCC Interruptor 2
UI
UI
MODEL
MODEL
Interruptor auxiliar Habedank Interruptor auxiliar Habedank
U
U
RARC RARC
T
T
Figura 123. Modelo para falla en extremo de la línea Figura 124. Modelo de falla en la mitad de la línea

85
Arias (2015)
Figura 125. Ajustes del ATP para las simulaciones de arco Figura 126. Variables definidas para falla en extremo A
secundario

86
Arias (2015)
Figura 127. Variables definidas para falla en extremo B Figura 128. Configuración interruptor asociado al extremo en falla

87
Arias (2015)
Figura 129. Configuración del interruptor en el extremo remoto Figura 130. Configuración del interruptor auxiliar

88
Arias (2015)
Resultados de la simulación y análisis. Los resultados de la 0
corriente de arco secundario para cada una de las fallas definidas [A]
se ilustran en la Figura 131, Figura 132 y Figura 133. -100
En la Figura 131, Figura 132 y Figura 133 se observa que la

-200
corriente nunca cruza por cero en un periodo de tiempo de 900 ms
debido a su gran componente de directa, excepto en el caso de la
-300
Figura 133, donde la corriente llega a cero y se extingue en un
tiempo mayor a los 750 ms. Esto no es aceptable, ya que el tiempo
-400
muerto definido para las protecciones es de 500 ms. No es
necesario aplicar los criterios del CESI ya que estos resultados
-500
son suficientes para afirmar que las posibilidades de extinción de 0.00 0.15 0.30 0.45 0.60 0.75 [s] 0.90
(file AS_LINE500_A.pl4; x-var t) c:X0002A-XX0019
la corriente de arco secundario son nulas. Como consecuencia de
esto y lo obtenido en el numeral 2.5, se concluye que no es Figura 131. Corriente de arco secundario. Falla en extremo A
posible habilitar el recierre monofásico en esta línea con estos
0
tiempos y condiciones operativas. La problemática podría [A]
solucionarse aumentando el tiempo de operación de las
-100
protecciones y/o añadiendo reactores de neutro a las líneas [31].
-200
-300
-400
-500
0.00 0.15 0.30 0.45 0.60 0.75 [s] 0.90
(file AS_LINE500_B.pl4; x-var t) c:X0002A-XX0019
Figura 132. Corriente de arco secundario. Falla en extremo B

89
Arias (2015)
50 La fase en falla abre en un tiempo de 8 ms (aproximadamente

[A]
25 medio ciclo).
0 Para la falla terminal, se aplica una falla trifásica en los bornes del
-25 interruptor. Las 3 fases de este abren en un tiempo de 8 ms (medio
-50
ciclo).
-75 Los tiempos de apertura definidos son indiferentes, ya que el
-100 interruptor siempre abre en el tiempo de cruce por cero (más
-125
próximo por encima del tiempo de apertura definido) de la
corriente.
-150
0.00 0.15 0.30 0.45 0.60 0.75 [s] 0.90
(file AS_LINE500.pl4; x-var t) c:X0002A-XX0020 El modelo de la falla monofásica y trifásica se hace con una
Figura 133. Corriente de arco secundario. Falla en la mitad resistencia lineal conectada en Y.
2.7 Verificación de TRV en interruptores de línea

La comprobación del TRV se hace en base a la norma
IEC 62271-100. En la Figura 135 se presentan las curvas
normalizadas de TRV correspondientes a los niveles de tensión
nominal (Ur: rated voltage) de 245 kV y 550 kV, de acuerdo al
tipo de falla (terminal ó kilométrica). No se verifica el TRV en
oposición de fases ya que es una condición con muy baja
probabilidad de ocurrencia [36].
interruptores. La comprobación del TRV se realiza con la
implementación de interruptores controlados por tiempo.
Para simular la falla kilométrica, se aplica una falla monofásica Figura 134. Localización de la resistencia que modela la falla
sólida a tierra a una distancia de 1 km del interruptor analizado. trifásica y monofásica

90
Arias (2015)
TRV normalizado
900
876; 817
800
700
672; 629
600
Tensión [kV]
500
219; 438
400 392; 364
168; 337
300
98; 195 300; 280
200
75; 150
100
0
0 200 400 600 800 1000 1200 1400 1600 1800 2000
Tiempo [μs]
Falla terminal - Ur 245 kV Falla terminal - Ur 550 kV

Falla kilométrica - Ur 245 kV Falla kilométrica - Ur 550 kV
Figura 135. TRV Normalizado en interruptores

91
Arias (2015)
Figura 136. Ajustes de resistencia para falla monofásica sólida Figura 137. Ajustes de resistencia para falla trifásica sólida
Note que la configuración de la resistencia lineal tiene un valor En este caso la configuración de la resistencia lineal tiene todos
muy pequeño (0,001 Ω) y otros dos valores muy elevados sus valores muy pequeños (0,001 Ω), lo que corresponde al
(1E21 Ω). Esta configuración corresponde al modelo de falla modelo de falla trifásica sólida a tierra.
monofásica sólida a tierra

92
Arias (2015)
Figura 138. Ajustes de interruptor para falla kilométrica Figura 139. Ajustes de interruptor para falla terminal

93
Arias (2015)
El óvalo rojo en la Figura 138 y Figura 139 indica que el Tratamiento de resultados. La variable de interés en esta ocasión
interruptor debe registrar la tensión, ya que esta es la variable de es la tensión que presenta el interruptor en el momento de su
interés en este caso y es la forma más cómoda de medirla. apertura ante la falla. La tensión entre bornes del interruptor debe
ser comparada con las curvas establecidas por la norma, ilustradas
en la Figura 135. La mejor manera de realizar una comparación
entre los resultados y la norma con el fin de verificar el
Interruptor analizado
LCC LCC cumplimiento del TRV, se logra superponiendo ambas curvas.
U
Como no se tienen las herramientas para lograr generar en el ATP
las curvas definidas en la norma y al mismo tiempo superponerlas
1. km
con los resultados de la simulación, debe hacerse uso del archivo
.lis. El archivo .lis contiene, además de la información que se ha
mostrado en numerales anteriores, la información de la gráfica, es
Falla monofásica - tierra
decir, contiene los puntos correspondientes a las variables
medidas para cada instante de tiempo de la simulación, lo que
Figura 140. Modelo para verificación de TRV. Falla kilométrica permite al usuario obtener los resultados gráficos de la simulación
no solo por parte del ATP. Los datos del archivo .lis permiten ser
llevados a otro software, como por ejemplo Excel, de tal manera
que puedan ser graficados junto con las curvas definidas por la
Interruptor analizado norma y lograr la comparación deseada.
LCC
U Para poder obtener del archivo .lis los valores de tensión en
bornes del interruptor punto a punto para cada instante de tiempo,
deben realizarse algunos ajustes en el ATP Settings.
Falla trifásica - tierra
Figura 141. Modelo para verificación de TRV. Falla terminal

94
Arias (2015)
El parámetro encerrado en el óvalo rojo en la Figura 142, viene

por defecto configurado con un valor de 500. Esto hace que el
muestreo de la gráfica registrado en el archivo .lis no sea punto a
punto sino que aumente en pasos de 500 ms, lo que genera
pérdida de información importante de la gráfica. Para lograr
obtener cada uno de los puntos y evitar que se pierda información
gráfica, debe modificarse ese valor a 1.
Figura 142. Ajustes del ATP Settings para obtener los valores de Figura 143. Información de la gráfica contenida en el archivo .lis
tensión punto a punto en el archivo .lis

95
Arias (2015)
En la Figura 143 se observa la información que contiene el TRV falla terminal - 220 kV
archivo .lis referente a la gráfica. Se tiene el paso de cálculo, el 400000
tiempo y los valores de cada una de las variables que fueron 350000
monitoreadas en la simulación. Estos datos son los que deben ser 300000
Tensión [V]
llevados a Excel. 250000
200000
El análisis para la falla kilométrica, debido a que se lleva a cabo 150000
100000
aplicando una falla monofásica, sólo se hace en la fase del 50000
interruptor encargada de despejar la falla, mientras que el análisis 0
para falla terminal, que se lleva a cabo aplicando una falla 0,0144 0,0145 0,0146 0,0147 0,0148 0,0149 0,015
Tiempo [s]
trifásica, debe hacerse en cada una de las fases del interruptor. Por
simplicidad, se hace el análisis en una sola fase, la que presenta el Tensión Límite Norma
peor comportamiento.
Figura 144. TRV interruptor línea a 220 kV, Fase B
Resultados de la simulación y análisis
Falla terminal. En la Figura 144 y Figura 145 se presentan los TRV falla terminal - 500 kV
resultados para falla terminal. Para el interruptor de 220 kV, se 900000
800000
hace el análisis en la Fase B y para el de 500 kV se hace en la 700000
Fase A. Ambas fases fueron las que peor comportamiento
Tensión [V]
600000
mostraron entre las tres. 500000
400000
300000
200000
100000
0
0,009 0,0095 0,01 0,0105 0,011
Tiempo [s]
Tensión Límite Norma
Figura 145. TRV interruptor línea a 500 kV, Fase A

96
Arias (2015)
En la Figura 144 se observa que el TRV cumple para el interruptor

de 220 kV, ya que la tensión no sobrepasa la magnitud ni el
crecimiento en el tiempo permitido por la norma, mientras que en
la Figura 145, correspondiente al caso en 500 kV, hay una
violación en el límite de la tensión permitida.
Falla kilométrica. En la Figura 146, Figura 147, Figura 148 y
Figura 149 se ilustran los resultados para la falla kilométrica. Se
observa una violación en el límite de tiempo de crecimiento de la
tensión (RRRV) para el interruptor a 220 kV. El TRV para el
interruptor a 500 kV se encuentra cerca del límite en cuanto a
magnitud, pero no viola la norma.
Figura 147. Punto de violación de la norma localizado en
TRV Línea 220 kV interruptor a 220 kV
300000
280000
250000
200000
Tensión [V]
150.000
150000
100000
50000
0
0.0085 0.009 0.0095 0.01 0.0105 0.011
Tiempo [s]
Figura 146. TRV interruptor línea a 220 kV, Fase A

Figura 148. Punto de violación de la norma ampliado

97
Arias (2015)
TRV falla Kilométrica - 500 kV En muchas situaciones el RRRV puede ser mejorado al incluir en
700000
las simulaciones las capacitancias parásitas de los equipos. Esto
600000 ayuda a reducir e inclusive a eliminar la necesidad de colocar
500000 capacitores en paralelo con el interruptor [38].
Título del eje
400000
En la norma IEEE C37.011-2011, pueden encontrarse algunos
300000 valores de capacitancias parásitas típicas de equipos dependiendo
200000 de las características eléctricas del sistema.
100000
Las resistencias de apertura se implementan principalmente en
0 interruptores antiguos de aire comprimido, por lo que no son muy
0,0087 0,0092 0,0097 0,0102 0,0107 0,0112 0,0117
Título del eje utilizadas en la actualidad [39].
Tensión interruptor Límite Norma
Para efecto de las simulaciones, se implementan las medidas 3 y 4
Figura 149. TRV interruptor línea a 500 kV, Fase A en las situaciones en las que se presenten problemas con el TRV.
Mitigación del TRV. Para corregir el TRV en las ocasiones en las Es importante resaltar que si existe alguna metodología para el
que no se cumplen los criterios definidos por la norma en cuanto a dimensionamiento de los equipos que se van a utilizar para
magnitud y rata de crecimiento (RRRV), pueden aplicarse varios mitigar el efecto del TRV, no se aplica, ya que esto no hace parte
métodos [37]: de los objetivos del documento. Se emplean descargadores de
sobretensión y capacitancias con características que permiten
Cambiar el interruptor por uno de mayor capacidad ilustrar el efecto de estos dispositivos sobre el fenómeno.
Uso de resistencias de apertura Para mejorar la condición en 500 kV, se coloca un descargador de
Uso de capacitancias en paralelo con el interruptor sobretensión en paralelo con el interruptor. El descargador
implementado en este caso es seleccionado de un catálogo de
Descargadores de sobretensión en paralelo con los Siemens. El equipo se describe más detalladamente en los
interruptores [38]. Anexos.
La medida 3 es implementada para reducir el RRRV y las tres
restantes para disminuir la magnitud de la tensión.

98
Arias (2015)
Descargador de En la Figura 151 se ilustran los resultados para el caso con el

sobretensiones descargador de sobretensión en paralelo con el interruptor. Se
en paralelo con el observa como la magnitud de la tensión disminuye con respecto a
interruptor R(i)
lo observado en la Figura 145 y los valores de la curva no son
R(i)
sobrepasados. Se verifica la efectividad de esta medida
R(i)
Para mejorar la condición ilustrada en la Figura 148, se coloca
LCC
una capacitancia de 9 nF, obtenida por tanteo en las simulaciones,
U
en paralelo con el interruptor. En la Figura 153 y Figura 154 se
ilustran los resultados para el caso con la capacitancia en paralelo
Falla con el interruptor.
Figura 150. Modelo interruptor más descargador para corrección Capacitancia de 9 nF

en paralelo con el interruptor
del TRV
TRV falla terminal - 500 kV - Corregido

900000 LCC
800000 U
700000
Tensión [V]
600000
500000 Falla
400000
300000
200000
100000
0 Figura 152. Modelo interruptor más capacitancias para corrección
0.009 0.0095 0.01 0.0105 0.011
Tiempo [s]
del TRV
Tensión Límite Norma En la Figura 154 se observa que luego de agregar los capacitores
al interruptor, el RRRV disminuye, quedando por debajo de la
Figura 151. TRV interruptor línea a 500 kV con descargador de norma y solucionando el problema. Se verifica la efectividad de la
sobretensiones medida.

99
Arias (2015)
TRV Línea 220 kV corregido CONCLUSIONES

300000
El uso de herramientas digitales se ha convertido en una
250000
necesidad para llevar a cabo un correcto planeamiento, diseño,
200000 dimensionamiento y operación de sistemas eléctricos, ya que
Tensión [V]
150000
estos programas contienen elementos que permiten simular de
forma muy precisa y acertada cada una de las variables que
100000
influyen en el comportamiento de una red determinada.
50000
Antes de ejecutar cualquier tipo de simulación es necesario tener
0
0.0085 0.009 0.0095 0.01 0.0105 0.011 0.0115 0.012 claros y definidos los parámetros eléctricos y geométricos de cada
Tiempo [s]
uno de los elementos del sistema, la frecuencia de operación de la
red y de los fenómenos a estudiar, ya que a partir de éstos
Figura 153. TRV interruptor línea a 220 kV con capacitancias parámetros se establecen datos de entrada como tiempos de
configuración de interruptores y el paso de cálculo del software,
datos críticos de los que depende obtener resultados congruentes
y cercanos a lo real.
Conforme a los resultados obtenidos en cada una de las
simulaciones, se observa que uno de los parámetros que más
influye en la severidad del transitorio es el voltaje de la red, ya
que los resultados más críticos se obtuvieron para la tensión de
500 kV.
Los métodos de mitigación de fenómenos transitorios como el uso
de resistencias de preinserción, mando sincronizado,
descargadores de sobretensión y capacitancias en paralelo con los
interruptores, muestran ser muy eficaces en la reducción del
Figura 154. Punto de violación de la norma ampliado. Caso con impacto de las maniobras sobre el sistema, lo que sugiere su
capacitores implementación para conservar en mayor medida la vida útil de

100
Arias (2015)
los equipos. Cabe resaltar que algunos de los métodos como por REFERENCIAS
ejemplo el uso de capacitancias y descargadores de sobretensión
en paralelo con los interruptores para mitigar el efecto del TRV, [1] L. Martinez Torres, «ESTUDIOS ESTADÍSTICO DE SOBRETENSIONES
POR MANIOBRA DE INTERRUPTORES EN LINEAS DE
pueden requerir de cierta metodología para ser calculados y TRANSMISIÓN EMPLEANDO EL PSCAD/EMTP,» Mexico, D.F,
dimensionados correctamente. 2012.
[2] C. Fuerte Esquivel y J. Guardado Zavala, «Introducción al Análisis de
Una de las maniobras más importantes en los sistemas de potencia
Transitorios Electromagnéticos en Sistemas Eléctricos de Potencia,»
es el recierre monofásico, ya que el 90% o más de las fallas Ciencia Nicolaita, nº 6, pp. 141-148, 1994.
presentadas en los sistemas eléctricos son monofásicas a tierra, [3] J. Das, Transients in Electrical Systems. Analysis, Recognition, and
por lo que no es suficiente realizar solamente un estudio de Mitigation, Mc Graw Hill, 2010.
sobretensiones transitorias, sino que es necesario ejecutar análisis [4] L. F. Saá, «ANÁLISIS DE TRANSITORIOS ELECTROMAGNÉTICOS EN
como el de la corriente de arco secundario para determinar o LA ENERGIZACIÓN DE ELEMENTOS DE TRANSMISIÓN
verificar el tiempo muerto de las protecciones, esto sumado a los UTILIZANDO EL PROGRAMA DIGSILENT POWER FACTORY,»
Quito, 2010.
estudios de estabilidad transitoria, que permiten verificar si con
[5] J. L. Torres, «CONTROL DE TRANSITORIOS DE MANIOBRA CON
los tiempos de operación definidos el sistema no pierde su
INTERRUPTORES CON MANDO SINCRONIZADO - APLICACIÓN
estabilidad y si definitivamente puede o no habilitarse la A BANCOS DE CAPACITORES,» Quito, 2003.
maniobra. [6] R. Folkers, «DETERMINE CURRENT TRANSFORMER SUITABILITY
USING EMTP MODELS».
[7] J. F. Piñeros Saldarriaga, «EL PROCESO DE LA SIMULACIÓN DE
REDES ELÉCTRICAS,» Asosiación de Ingenieros Electricistas AIE
UdeA, pp. 45-52.
[8] A. Greenwood, Electrical Transients y Power Systems, John Wiley & Sons,
Inc., 1991.
[9] D. P. Landazuri y M. Rumipamba, «ESTUDIO Y PROPUESTA DE UN
TIPO DE INTERRUPTOR A SER IMPLEMENTADO EN LA LÍNEA
DE 500 kV MEDIANTE MODELACIÓN EN ATP,» Quito, 2011.
[10] I. A. Isaac Millan, Interviewee, Notas de clase Alta Tensión, UPB.
[Entrevista]. 2014.
[11] J. Navarro Matarrita, «Transitorios Electromagnéticos en los Sistemas de

101
Arias (2015)
Potencia,» 2012. SOBREVOLTAJES A FRECUENCIA INDUSTRIAL Y POR

[12] A. D´Ajuz, C. dos Santos Fonseca, F. M. Salgado Carvalho, J. Amon Filho, MANIOBRA EN LÍNEAS DE TRANSMISIÓN A NIVEL DE 500 kV,»
L. E. Nora Dias, M. P. Pereira, P. C. Vez Esmeraldo, R. Vaisman y S. de Quito, 2010.
Oliveira Frontin, TRANSITÓRIOS ELÉTRICOS COORDENAÇÃO de [26] M. Oliveira, A. Bretas, O. Perrone, J. Reversat, H. Muñoz y R. Orellana
ISOLAMENTO, Rio de Janeiro, 1987. Paucar, «Principales características de la corriente Inrush en
[13] A. Rifaldi, «Electronic Edition of the ATP Rulebook inPDF Format,» vol. 3, Transformadores».
nº 4, 2001. [27] V. Marines Castillo y G. Idarraga Ospina, «Análisis de los modelos de
[14] Electromagnetic Transients Program (EMTP) Theory Book. transformadores para la simulación de la protección diferencial,» vol.
XIII, nº 46, 2010.
[15] L. Dután, «Modelación de Transformadores para Estudios de Energización
usando ATP-EMTP». [28] IEB, «Coordinación de Aislamiento. Norma IEC 60071-2 vs Ingeniería por
desempeño,» 2009.
[16] L. Restrepo, G. Caicedo y F. Castro, «Modelos de línea de transmisión para
transitorios electromagnéticos en sistemas de potencia,» vol. 16, nº 1, [29] G. Ziegler, «Transformer Protection,» 2007-2008.
2008. [30] J. Vivas, V. Mendoza y F. Ferri, «ESTUDIO DETALLADO DE LA
[17] «IEC TR 60071-4. Insulation co-ordination - Part 4: Computational Guide to MANIOBRA DE DESPEJE TRIFÁSICA EN SISTEMAS DE
insulation co-ordination and modelling of electrical networks,» 2004. TRANSMISIÓN».
[18] N. Watson y J. Arrillaga, Power Systems Electromagnetic Transients [31] J. H. Vivas Nava, «CONSIDERACIONES ACERCA DE LOS ESQUEMAS
Simulation, London, 2007. DE RECONEXIÓN EN SISTEMAS DE TRANSMISIÓN,» 2002.
[19] M. Sánchez Tenorio, «Cálculos para la selección del pararrayos». [32] S.-H. Sohn, G.-J. Cho, J.-K. Park, Y.-S. Oh, C.-H. Kim, W.-J. Kim, H.-J. Oh,
J.-J. Yang, T. Senjyu y T. Funabashi, «Analysis of Secondary Arc
[20] «Zensol. Construcción Sólida, Respaldo Sólido,» [En línea]. Available: Extinction Effects according to the Application of Shunt Reactor and
http://www.zensol.com/es/pruebas-de-tiempo-analisis-de-la- High Speed Grounding Switches in Transmission Systems,» vol. 4, nº 4,
informacion#4221. [Último acceso: 3 Marzo 2015]. 2014.
[21] G. Susanibar Celedonio, «TRANSITORIOS ELECTROMAGNÉTICOS. [33] J. H. Vivas N, V. Mendoza M. y A. L. Villadar V., «ESTUDIO DE
Datos Necesarios para Realizar Simulaciones en ATP». SENSIBILIDAD DE FACTORES QUE INFLUYEN EN EL TIEMPO
[22] FURNAS/UFF, Transitórios Elétricos e Coordenação de Isolamento - DE EXTINCIÓN DEL ARCO SECUNDARIO».
Aplicação em Sistemas Elétricos de Alta Tensão, Rio de Janeiro, 1987. [34] D. Cubillos Rodríguez y L. J. Ibáñez Serrano, «DESCRIPCIÓN TEÓRICA Y
[23] R. M. Jimenez Velez, «ESTUDIO DE LA COORDINACIÓN DEL ANÁLISIS DE FALLAS DE ARCO MEDIANTE UN MODELO
AISLAMIENTO DE LA SUBESTACIÓN OBRAPIA - EMPRESA ADECUADO, CON EL PROPÓSITO DE DETERMINAR SU EFECTO
ELÉCTRICA REGIONAL DEL SUR S.A. - LOJA,» Cuenca, 2007. EN LÍNEAS DE TRANSMISIÓN,» Bogotá, 2008.
[24] Á. Juan, M. Sedano y A. Vila, «La Distribución Normal». [35] Schneider Electric, «Cuaderno Técnico No. 158. Cálculo de corrientes de
[25] K. F. Sailema Morales y X. A. Proaño Maldonado, «ANÁLISIS DE cortocircuito».

102
Arias (2015)
[36] «IEC 62271-100. High-voltage switchgear and controlgear - Part 100: AUTOR
Alternating-current circuit breakers,» 2008.
[37] W. V. S. Azevêdo, W. L. A. Neves, D. Fernandes Jr., B. A. Souza y R. M. Cristian Camilo ARIAS RODRÍGUEZ, nacido
Soares, «Device to Limit Transient Recovery Voltage». en Medellín, Colombia. Bachiller egresado del
[38] CIGRÉ Working Group A3.22, «Tenchnical Requirements for Substation Colegio Salesiano Santo Domingo
Equipment Exceeding 800 kV AC.,» 2008. Savio (2008). Diplomado en Gerencia de
[39] ABB, «Interruptores de Tanque Vivo. Guía para el comprador,» 2005.
Proyectos de la Universidad Pontificia
Bolivariana (2013). Actualmente, egresado no
[40] ATP Rulebook.
graduado de la Facultad de Ingeniería Eléctrica
de la Universidad Pontificia Bolivarana,
integrante del área de Estudios Eléctricos de la consultora HMV
Ingenieros ltda y docente en la Institución Educativa Politécnico
Interamericano del municipio de La Ceja.

103
Arias (2015)
ANEXOS
Tabla 16. Parámetros de las líneas objeto de las simulaciones
Línea Tensión R0 X0 B0 R1 X1 B1 Longitud Resistividad
(kV) (Ω/km) (Ω/km) (Ω/km) (Ω/km) (Ω/km) (ʊ/km) (km) (Ω*m)
A-B 500 2,84E-01 8,33E-01 3,10E-06 2,21E-02 3,16E-01 5,21E-06 150 400
D-F 220 2,78E-01 1,05 2,33E-06 5,46E-02 5,29E-01 3,15E-06 103,35 100
0 0
1 2 3
Figura 155. Configuración geométrica líneas a 500 kV
0 0
Figura 156. Configuración geométrica líneas 220 kV

104
Arias (2015)
Tabla 17. Características generales de los descargadores de ZnO

Descripción Unidad Tensión máxima del sistema (kV)
245 550 550
Tensión nominal (Ur) kV 198 444 396
Tensión continua de operación (Uc) kV 156 355 316
Corriente nominal de descarga (In) kA 20 20 20
Clase de descarga de línea Clase 5 5 5
Capacidad de disipación de energía kJ/kV 13 Ur 25 Ur 25 Ur
Capacidad de disipación de energía. MJ 2,57 11,10 9,90
Ubicación Línea 220 kV Línea 500 kV Interruptor 500 kV
Tabla 18. Características de protección a maniobra

Característica Tensión nominal Ur (kV)
Corriente (30/60 μs) 198 444 396
(A) (kVpico) (kVpico) (kVpico)
0.001 370000 814000 736000
500 380000 856000 754000
1000 390000 872000 772000
2000 409000 909000 810000
Característica de protección a maniobra

1000000
900000
800000
700000
Tensión [V]
600000
500000
400000
300000
200000
100000
0
0 0.0001 500 1000 2000
Corriente [I]
Ur=198 kV Ur=444 kV Ur=396 kV
Figura 157. Característica de protección a maniobra
Tabla 19. Parámetros básicos de la transformación

Parámetro Alta Media Baja
Capacidad (MVA) 750 750 250
Conexión YN YN D
Tensión, kV 500 220 33

105
Arias (2015)
Tabla 20. Parámetros eléctricos entre devanados

Base Impedancia Z
Devanados
(MVA) (%)
H-L 750 21,25
H-T 250 15,83
L-T 250 8,33
Característica de saturación
1800
1600
1400
Flujo [Wb-T]
1200
1000
800
600
400
200
0
-100 100 300 500 700 900 1100 1300 1500
Corriente [A]
saturación
Figura 158. Característica de saturación transformador 500 kV/220 kV/33 kV
Figura 159. Característica de saturación transformador 500 kV/220 kV/33 kV en el ATP
Tabla 21. Parámetros de los reactores de línea a 500 kV

Tensión nominal Potencia Tipo de
(kV) (Mvar) aislamiento
500/√3 156 Gradual

Metodología para El Análisis de Transitorios Electromagnéticos de Maniobra en Sistemas de Alta t2

Cargado por

Copyright:

Formatos disponibles

Metodología para El Análisis de Transitorios Electromagnéticos de Maniobra en Sistemas de Alta t2

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Metodología para El Análisis de Transitorios Electromagnéticos de Maniobra en Sistemas de Alta t2

Cargado por

Copyright:

Formatos disponibles

Metodología para el análisis de transitorios electromagnéticos de maniobra en sistemas de

Trabajo de grado para optar al título de Ingeniero Electricista

Universidad Pontificia Bolivariana

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

Figura 30. Ajustes del Symmetric RL coupled line................................................................................................................................................................ 37

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

Figura 123. Modelo para falla en extremo de la línea ........................................................................................................................................................... 84

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

INTRODUCCIÓN De maniobra, generados por apertura y cierre de

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

modelo al incluir las pérdidas por resistencia se torna más

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

De la ecuación (18) y la ecuación (20) es posible escribir la onda

Figura 3. Reflexión y Refracción de la onda incidente (b)

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

ZA≠ZB 1.5 Maniobras ejecutadas en Sistemas Eléctricos de Potencia.

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

Eq. K Las maniobras generadoras de transitorios simuladas son:

Energización de transformador 220 kV/500 kV

Transformador 220 kV/500 kV

simulaciones se expresan en p.u. En la Tabla 1 y Tabla 2 se

Tensión devanado Corriente

Subestación H Subestación J nominal

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

Tabla 4. Elementos seleccionados para la ejecución de las

-Línea D – F 220 kV (izquierda)

-Línea D – F 220 kV (izquierda)

Análisis de arco secundario -Línea A – B 500 kV

-Línea D – F 220 kV (izquierda)

Figura 7. Ajustes del ATP para las simulaciones

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

2.2 Modelos implementados para la elaboración de la red

Figura 8. Localización del modelo de línea

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

En la Figura 9 se ilustra la configuración para el modelo de línea.

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

Figura 12. Representación gráfica del modelo BCTRAN

Para evitar que se presenten problemas numéricos en la

Figura 13. Localización del modelo BCTRAN

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

Figura 15. Localización del inductor no lineal para representación

Facultad de Ingeniería Eléctrica y Electrónica Universidad Pontificia Bolivariana

Figura 16. Representación gráfica del inductor no lineal tipo 98