SciELO - Scientific Electronic Library Online

 
vol.19 número2Módulo de Software en Lenguaje Python para Estudios de Flujo de Potencia InercialAnálisis de Indicadores de Calidad para la Actualización de la Regulación No. CONELEC - 003/08 Calidad del Transporte de Electricidad y del Servicio de Transmisión y Conexión en el Sistema Nacional Interconectado índice de autoresíndice de assuntospesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados

Journal

Artigo

Indicadores

Links relacionados

  • Não possue artigos similaresSimilares em SciELO

Compartilhar


Revista Técnica energía

versão On-line ISSN 2602-8492versão impressa ISSN 1390-5074

Revista Técnica energía vol.19 no.2 Quito Jan./Jun. 2023

https://doi.org/10.37116/revistaenergia.v19.n2.2023.554 

SISTEMAS ELÉCTRICOS DE POTENCIA

Análisis Estadístico de Sobretensiones por Maniobra de energización en Líneas de Transmisión de Extra Alta Tensión

Statistical Analysis of Overvoltages due to Energizing Switching in Extra High Voltage Transmission Lines

1Universidad Técnica de Cotopaxi, Latacunga, Ecuador E-mail: jhon.calle3254@utc.edu.ec ;

2 Universidad Técnica de Cotopaxi, Latacunga, Ecuador E-mail: wilian.guaman8956@utc.edu.ec


Resumen

En este trabajo se aplican las distribuciones Gaussiana y Weibull modificada para obtener las sobretensiones representativas de maniobras de energización en líneas de Extra Alto Voltaje y se comparan los resultados obtenidos en las sobretensiones representativas fase-tierra y fase-fase para tres escenarios de energización. Además, se expone la metodología para realizar un estudio estadístico de sobretensiones en el software ATPDraw, tomando como caso de estudio la interconexión Ecuador-Perú 500 kV. Los resultados demuestran que ambas distribuciones son aplicables para un análisis estadístico de sobretensiones si se consideran magnitudes de voltaje con igual probabilidad de aparecer en el sistema.

Palabras clave: Sobretensiones; Weibull; Gauss; Coordinación de aislamiento; ATPDraw.

Abstract:

In this paper, the Gaussian and modified Weibull distributions are applied to obtain the representative overvoltages of energizing switching operations in Extra High Voltage lines, and the results obtained in the representative phase-to-ground and phase-to-phase overvoltages for three energizing scenarios are compared. In addition, the methodology to perform a statistical study of overvoltages in the ATPDraw software is presented, using the Ecuador-Peru 500 kV interconnection as a case study. The results show that both distributions are applicable for statistical analysis of overvoltages if voltage magnitudes with equal probability of appearing in the system are considered.

Keywords: Overvoltages; Weibull; Gauss; Insulation Coordination; ATPDraw

INTRODUCCIÓN

La implementación de la interconexión Ecuador-Perú surge con el compromiso de la integración eléctrica regional y aprovecha la complementariedad hidrológica para realizar intercambios de energía entre ambos países. También se garantiza la integración de los futuros proyectos de generación que se desarrollan tanto en Ecuador como en Perú [1]. Esta interconexión se realizará a nivel de 500 kV desde la S/E Chorrillos hasta la S/E Pasaje en Ecuador y culminando en la S/E Piura en Perú.

La amplitud de las sobretensiones aumenta proporcionalmente en base al voltaje nominal de la red [2]. En sistemas de Extra Alta Tensión (EHV por sus siglas en inglés) las sobretensiones que podrían causar daños significativos, es el caso de la interconexión Ecuador-Perú 500 kV. Luego, la interconexión aumenta la posibilidad de que ocurran afectaciones en la infraestructura eléctrica que resulten de los sobrevoltajes en la red. Estas solicitaciones dieléctricas son alteraciones de potencial de un sistema con tendencia a superar su valor nominal, causadas principalmente por fallas a tierra, eventos de maniobra y descargas atmosféricas [2]. La exposición a esta elevación de tensión produce un desgaste en el aislamiento de los equipos que componen el sistema eléctrico de potencia (SEP).

Debido a la gran cantidad de parámetros que intervienen durante las sobretensiones de maniobra y su carácter aleatorio, es común utilizar métodos estadísticos para determinar la amplitud de las solicitaciones dieléctricas resultantes con alta probabilidad de aparecer en el sistema al momento de energizar la red, considerando como variable aleatoria continua a las sobretensiones que puedan ocurrir. La distribución normal o Gaussiana es la más utilizada en este tipo de análisis, tal y como se presenta en [3], [4] donde se aplica esta distribución para analizar sobretensiones por maniobra en EHV obteniendo resultados cercanos a los propuestos en las normativas IEEE 1313 e IEC 60071 [5], [6].

El trabajo de [7] sugiere utilizar la distribución de Weibull en lugar de la distribución Gaussiana, debido a que las expresiones de Weibull modificadas están definidas por los mismos parámetros que caracterizan a las expresiones gaussianas (media y desviación estándar). Además, en el Anexo C de la IEC 60071-2 [6], se recomienda utilizar la función de probabilidad modificada de Weibull tanto para el cálculo de sobretensiones como para determinar la descarga disruptiva de aislamiento. A pesar de esto, la mayoría de estudios se enfocan en utilizar la distribución Gaussiana para el cálculo de sobretensiones y en general para distintitas variables continuas. Sin embargo, se debe tener especial cuidado al suponer para una situación dada un modelo de probabilidad normal sin previa comprobación o comparación. Por ello, diversos estudios tales como [8], han comparado distintas distribuciones de probabilidad existentes para evaluar diversos fenómenos en su propio campo, por ejemplo, en la evaluación del potencial de energía eólica. llegando a la conclusión de que ambas distribuciones arrojan resultados aceptables.

En estudios de coordinación de aislamiento se requiere conocer la distribución estadística de sobretensiones y la confiabilidad de sus parámetros. En este trabajo se compara la aplicación de las distribuciones estadísticas Gaussiana y de Weibull modificada, a partir de las sobretensiones generadas en maniobras de energización mediante la herramienta de simulación ATPDraw, con la finalidad de determinar la magnitud de las sobretensiones fase-fase y fase-tierra con alta probabilidad de presentarse. Las principales contribuciones de este trabajo se detallan a continuación:

Se presenta un modelo la interconexión Ecuador-Perú 500 kV usando el software ATPDraw.

Se muestra el procedimiento para desarrollar la simulación de sobretensiones de maniobra por energización en el sistema de EHV de 500 kV.

Se comparan los resultados obtenidos a partir de las distribuciones estadísticas Gaussiana y de Weibull para sobretensiones de maniobra por energización.

El resto de este documento se organiza como sigue, en la Sección II, se revisan los fundamentos teóricos para el estudio estadístico de sobretensiones. Luego, en la Sección III se muestra la metodología utilizada para modelar el sistema y simular las sobretensiones de maniobra. En la Sección IV se presentan los resultados del análisis estadístico de las sobretensiones en la red. Finalmente, en la Sección V se exponen las conclusiones del trabajo.

MARCO TEÓRICO

En la coordinación de aislamiento es fundamental determinar la magnitud de las sobretensiones para seleccionar apropiadamente el aislamiento eléctrico que deberá estar en capacidad de soportar las solicitaciones dieléctricas durante la operación del sistema. A continuación, se revisan los fundamentos teóricos necesarios para estudiar las sobretensiones de la red.

Ondas Viajeras

Cuando se producen perturbaciones en las líneas de transmisión ya sea por fallas, maniobras o descargas atmosféricas, se generan ondas de corriente y voltaje, a estas ondas resultantes se las conoce como ondas viajeras, las cuales presentan velocidades aproximadas al desplazamiento de la luz con tiempos de propagación en el orden de los us (microsegundos). Se puede efectuar el cálculo del tiempo de propagación de onda en función de la longitud de la línea de transmisión evaluada mediante la Ecuación (1) [9].

Las ondas viajeras se propagan a lo largo de la línea hasta llegar a un punto de discontinuidad que está determinado por un cambio en la estructura geométrica o medio de la línea de trasmisión, este cambio produce una variación de parámetros de impedancia característica de la línea transmisión 𝑍 𝑐 , que determina el desplazamiento de la onda, esto se ilustra en la Fig. 1, donde se representan el comportamiento de las ondas viajeras al encontrarse con el punto de discontinuidad y las impedancias características ZA y ZB separadas por el punto de discontinuidad. Es en este punto de discontinuidad la onda viajera inicial o incidente Vi se divide en onda reflejada Vr y transmitida Vt, que obedecen a las expresiones descritas en la Ecuación (2) y Ecuación (3) respectivamente [10].

Entonces, la magnitud de la onda transmitida y reflejada depende del valor de ZA y ZB. Luego, si se evalúa el caso más severo de reflexión de onda, donde 𝑍 𝐵 →∞, se presenta un circuito abierto y consecuentemente la onda es reflejada en su totalidad, adicionándose a la onda incidente de tensión, lo que representa una significativa sobretensión en el extremo abierto de la línea [10].

Fuente: [11]

Figura 1: Ondas viajeras de corriente y voltaje 

Sobretensiones por Maniobra

Las sobretensiones vienen definidas por su amplitud, tiempo de impulso y frecuencia de oscilaciones; y pueden ser de origen atmosférico, de maniobra o temporales como se identifica en la Figura 2. Las sobretensiones por maniobra (SOV por sus siglas en inglés) se usan para dimensionar el aislamiento en sistemas con voltaje nominal superior a 242 kV. Su origen es interno, de naturaleza oscilatoria y ocasionadas generalmente por la acción de apertura o cierre de los interruptores de potencia [12].

Fuente: [13]

Figura 2: Clasificación de Sobretensiones 

Los parámetros de estudio de sobretensiones dependen de la capacidad de cortocircuito en el lugar donde se ha conectado la línea y del instante en el que ocurre la unión en los contactos de los interruptores. Este instante es una variable aleatoria, por lo que el nivel de sobrevoltaje en el momento de encendido de la línea también es una variable aleatoria [3]. Las sobretensiones por maniobra se dan en líneas de transmisión, transformadores, bancos de reactores o condensadores, donde destacan los eventos de energización y reenganche de estos elementos [9].

Distribución Estadística de Sobretensiones

Para un correcto análisis estadístico se deben considerar a las sobretensiones como una variable aleatoria continua, ya que puede tomar un número infinito de valores de voltaje entre un intervalo definido [15]. Esto sumando a la gran cantidad de variables a analizar, resulta complejo caracterizar, por lo que, para tener una aproximación mayor al comportamiento real de las sobretensiones se suelen aplicar distribuciones estadísticas, también denominadas distribuciones de probabilidad [2]. Generalmente, un análisis estadístico de sobretensiones debe considerar tres aspectos:

Selección de una distribución estadística para todas las variables de carácter aleatorio.

Dependiendo de la sobretensión a calcular, se desarrolla un modelo matemático para el sistema de estudio.

Aplicación de un modelo de estudio estadístico, donde usualmente se usa el método de Monte Carlo que es un procedimiento iterativo de variables aleatorias.

La distribución estadística de sobretensiones originada por eventos de maniobra en líneas aéreas se obtiene a partir de un número determinado de maniobras generadas en simulación utilizando un interruptor estadístico en ATPDraw, estos resultados se ordenan para obtener un histograma similar al mostrado en la Figura 3.

Fuente: [2]

Figura 3: Histograma de sobretensiones 

Al reorganizarse el histograma de sobretensiones a partir de 𝑛 simulaciones, se obtiene una función de densidad de probabilidad [4], que si se considera una distribución Gaussiana, obedece a la siguiente expresión:

La Ecuación (4) representa la probabilidad de que se presente un determinado valor de tensión dentro de las muestras obtenidas, que resulta en la conocida campana de Gauss y sirve para calcular el riesgo de falla (𝑅) cuando se lo compara con los valores de la distribución acumulativa de aislamientos [2]. Sin embargo, la distribución estadística de sobretensiones también se puede expresar a través de una función complementaria denominada función de distribución acumulativa, presentada en la Ecuación (5), la cual representa la probabilidad de que la sobretensión resultante sea de una amplitud mayor a un valor de tensión dado.

No existe una expresión matemática analítica cerrada para la función de distribución acumulativa gaussiana evaluada con limites + ∞ 𝑦 − ∞, por lo que en la práctica se suelen acotar en intervalos limitados de − + 4 𝜎 a partir del valor medio o a su vez valorar la función de distribución acumulativa con limites desde − ∞ hasta el valor de sobretensión dado, este procedimiento lo realizan programas especializados en análisis de datos [2]. En la Figura 4 se muestra la representación gráfica de los valores de sobretensión por maniobra de las expresiones (4) y (5).

Figura 4: Caracterización de sobretensiones a través de sus funciones de distribución de probabilidad 

Debido a la simplicidad del procedimiento, con una precisión aceptable en los resultados, la utilización de la distribución normal o Gaussiana es empleada en todo tipo de aplicaciones que estén involucradas con variables continuas de naturaleza oscilatoria, en este caso esta variable es el valor de sobretensiones. Alternativamente, al considerar variables de este tipo también se puede emplear la función de densidad de probabilidad de Weibull [15], mostrada en la Ecuación (6):

Donde:

Cuando se considera valores de 𝜆=0 y 𝛼=1, la forma de la función de Weibull se asemeja a una campana, similar a la forma clásica de la distribución gaussiana, esto se da para valores positivos de 𝛽>2, como se muestra en la Fig. 5. Estos valores son recomendados por [2] debido a que las sobretensiones por lo general presentan un pico de resultados y luego decrecen, asemejándose a un campana estadística tradicional.

Figura 5: Distribución de Weibull, función de densidad de probabilidad 

También, se puede emplear una distribución de Weibull modificada en función de los parámetros de 𝜇 y 𝜎. En la Ecuación (7) se muestra la función de distribución acumulativa, que parte de la función de densidad de probabilidad de Weibull, adaptada por [2] empleando el parámetro 𝑍.

Donde, 𝑍= 𝑉−𝜇 𝜎 , es el parámetro adaptado para la distribución de Weibull modificada, y 𝑉 es el valor de la tensión a ser evaluada en la función de distribución acumulativa de Weibull.

METODOLOGÍA

En este documento se toma como caso de estudio el Sistema de 500 kV de Ecuador. Para modelar los elementos del SEP en el software ATPDraw se tomaron los valores de referencia expuestos en [1] y el procedimiento de coordinación de aislamiento de la IEC 60071-2 [13]. Para la evaluación de los resultados de la distribución estadística de probabilidad, en [6] se proponen dos métodos:

Método Fase-Pico: En cada operación de maniobra se escoge el valor pico fase tierra de cada fase que compone la red, lo que contribuye con tres resultados de valores pico.

Método Caso-Pico: De cada operación de maniobra se escoge el valor pico fase-tierra más alto de las tres fases.

Una vez seleccionado el método, se debe tener en cuenta que, sin pararrayos la sobretensión viene dada por la amplitud, con una probabilidad 𝑓𝑠(𝑉)= 2%, de ser sobrepasada o del 98% de ocurrencia si se usa 𝐹𝑠 𝑉 [6]. A partir de estas consideraciones, se obtienen los valores de las sobretensiones 𝑉 𝑒2 (fase-tierra) y 𝑉 𝑝2 (fase-fase). También, se considera la desviación estándar 𝜎, el valor medio 𝜇 y el valor de sobretensiones de truncamiento, que es la más elevada que se puede presentar en el sistema tanto 𝑉 𝑒𝑒 (fase-tierra) como 𝑉 𝑝𝑡 (fase-fase) representados en las Ecuaciones (7) y (8) respectivamente [6].

Modelos de Líneas de Transmisión

Para determinar el modelo de línea más adecuado para la simulación, se consideran cuatro elementos:

exactitud del modelo.

tiempo de respuesta computacional.

la capacidad del modelo de simular líneas con alta simetría (transposiciones).

la complejidad del modelo para el estudio.

Además, se debe tener en cuenta que, para estudios de sobretensiones se consideran parámetros distribuidos dependientes de la frecuencia, debido a que, en las líneas de transmisión se presentan fenómenos como el efecto skin y el efecto de retorno a tierra a medida que aumenta la frecuencia, para así evaluar de manera más precisa la contribución natural de las perdidas en la línea de transmisión [16].

Existen diferentes modelos de líneas de transmisión utilizadas en estudios eléctricos como lo son: modelo Z, idempotente, directo Nguyen, modelo PI, Bergeron y J Martin, sin embargo, en el software ATP se tiene disponibles los tres últimos por lo que se hará especial énfasis en estos. El modelo PI es usado generalmente para simulaciones en estado estacionario de líneas cortas a la frecuencia de la red, por lo que su aplicación se limita a sobretensiones temporales. En cuanto al modelo Bergeron, no es aplicable en este análisis debido a que utiliza matrices de transformación cuyos parámetros contienen ciertas limitaciones para variar en función su frecuencia. Finalmente, el modelo J. Martin es el más utilizado para análisis de sobretensiones por maniobra, ya que ofrece una mayor exactitud en comparación con el modelo Bergeron, a pesar de que, se tiene una respuesta computacional más lenta, y es aplicable para los rangos de frecuencia establecidos en la Tabla 1. Cabe recalcar que, este modelo no es aplicable en bajas frecuencias y para secciones de línea muy cortas. Consecuentemente el modelo J. Martin es usado como modelo estándar de comparación para este tipo de análisis [17].

Tabla 1: Clasificaciones de modelos de líneas de transmisión [17] 

En la Figura 6 se muestra la comparación de la magnitud de voltaje durante una maniobra de energización para los tres modelos de línea, considerando una longitud de 100 km. Donde se puede apreciar que el número de oscilaciones en un ciclo es de 11 en el modelo Bergeron y J. Martin. Nótese que, la magnitud de tensión inferior en el modelo Bergeron. Mientras que, el modelo PI presenta un menor número de oscilaciones con respecto a los demás modelos, pero su magnitud de voltaje no varía considerablemente con respecto al modelo J. Martin.

Figura 6: Comparación de modelos de línea PI, J Martin y Bergeron 

Por otro lado, como se señaló en la Sección 2.1, la geometría de la línea es un factor de relevancia para la propagación de ondas viajeras. En este estudio las dimensiones referenciales de la estructura de transmisión doble circuito para 500 kV se obtuvieron del trabajo de [18] y se muestran a detalle en la Figura 7. Además, de acuerdo a [1], [19] solamente el primer circuito será energizado, hasta evaluar el comportamiento de la red incluyendo la futura entrada al SNI del Sistema de transmisión Santiago.

Figura 7: Geometría de estructura doble circuito de 500 kV [18]. 

En la Tabla 2, se detallan los parámetros por fase de la línea de transmisión a ser utilizados en la Interconexión Ecuador-Perú 500kV, estos datos son necesarios para determinar valores de impedancia y admitancia en las líneas de transmisión.

Tabla 2: Parámetros de las líneas de Transmisión SNI 500 kV [1] 

Además, la línea presentará transposición de fases. La transposición en líneas con longitudes comprendidas entre 70 y 250 km se debe realizar en un ciclo completo que se compone de tres transposiciones simples en las longitudes 1/3, 3/6 y 5/6 de la línea, mientras que para longitudes mayores a 250 km la transposición se recomienda realizarla en dos ciclos completos que consiste en transposiciones simples en las longitudes 1/3, 3/12, 5/12, 7/12, 9/12, 11/12 [20].

Interruptores de Potencia

La configuración del interruptor de potencia es la pieza clave para un estudio estadístico de sobretensiones en el software ATPDraw, destacándose dos tipos de interruptores con tiempos de ejecución variables y útiles para análisis de sobretensiones [2]:

Interruptor Sistemático: el cierre es calculado de manera ordenada de un instante de tiempo mínimo a un tiempo máximo en incrementos de tiempo semejantes.

Interruptor Estadístico: el tiempo de cierre es calculado de manera aleatoria de acuerdo al tipo de distribución probabilística seleccionada.

Para la configuración del interruptor estadístico se deben considerar dos parámetros suministrados por el usuario, que son: tiempo medio de cierre o apertura (T) y la desviación estándar de tiempos de corte (σ). La selección del tipo de distribución probabilística (gaussiana o distribución uniforme) implica cambios en la dispersión de datos de simulación. Por ejemplo, para una distribución gaussiana se considera 𝑇+𝜎, y para una distribución uniforme 𝑇+ 1 2𝜎 3 , como se aprecia en la Figura 8. Así, se obtiene la variabilidad de datos que tendrá el programa al momento de ejecutar la simulación [15].

Si se escoge una distribución de tipo Gaussiana para evaluar las maniobras, la mayor cantidad de maniobras estará cerca del tiempo medio de cierre dependiendo del valor de desviación estándar. Mientras que, la distribución uniforme reparte la misma cantidad de maniobras a los extremos del tiempo medio en función de la desviación estándar que en este caso toma el nombre de rango de tiempo de amplitud de la distribución uniforme.

Fuente: [2]

Figura 8: a) Función de distribución gaussiana y b) Función de distribución uniforme 

Mientras que, para representación de interruptores se tiene en cuenta la relación que tienen los interruptores estadísticos entre sí, se pueden representar de tres maneras diferentes [21]:

Interruptores independientes: usados en simulaciones monofásicas ya, que el tiempo de cierre ( 𝑇 𝑐𝑖𝑒𝑟𝑟𝑒 ) es aplicable únicamente para el interruptor seleccionado.

Interruptores maestros y esclavos: se tiene un interruptor principal con configuración de 𝑇 𝑐𝑖𝑒𝑟𝑟𝑒 inicial, mientras que los interruptores esclavos dependen del tiempo de cierre del interruptor master y también de un tiempo de retardo, mostrados en la Ecuación 9. Se usa en unidades trifásicas, debido a que sus polos suelen estar unidos mecánicamente. Puede haber un retraso mínimo en el tiempo de cierre de los interruptores dependientes o esclavos.

Interruptores referencias (target): el tiempo medio de cierre de un interruptor se toma como referencia para otros interruptores relacionados con él.

En este estudio se ha seleccionado para cada fase un interruptor estadístico con distribución uniforme, formando así una sola unidad trifásica como se muestra en la Figura 9.

Figura 9: Interruptor Estadístico Trifásico en ATPDraw 

Una aproximación correcta de la distribución estadística de sobretensiones, resulta de someter al sistema a 200 simulaciones para obtener una precisión mayor en el cálculo de sobretensiones por maniobra como se propone en [2] y en [22].

Banco de Reactores Shunt y Reactor Neutro

Para compensar los efectos capacitivos que tiene la línea de transmisión, cuando esta se encuentre en vacío o con baja carga se instalan bancos de reactores en paralelo con la línea de transmisión para limitar las sobretensiones, aunque también se pueden conectar directamente en las barras para limitar sobretensiones por apertura y re cierre de líneas. A diferencia de los reactores de barra, los reactores de línea cuentan con un reactor neutro dimensionado en función de las características de la línea [20]. Los parámetros de compensación de línea, de barra y reactor neutro se presentan en la Tabla 3.

Tabla 3: Parámetros de compensadores [1], [20] 

El valor de reactor neutro se basa en la Ecuación (10), y cabe resaltar que la compensación debe ser realizada en un grado menor al grado critico de compensación presentado en la Ecuación (12), esto se realiza para evitar efectos de resonancia eléctrica [10].

Equivalente de Red y Configuraciones Adicionales

Resulta poco práctico modelar el SNI completo, y por esta razón para la mayoría de los cálculos se utilizan equivalentes de red, que consisten en la reducción de la red a partir de cálculos de cortocircuito en las barras seleccionadas [23]. Estos valores son clave para realizar simulaciones que proporcionen resultados coherentes y basados en los parámetros del sistema de estudio.

El equivalente de red en la barra Chorrillos se obtuvo a partir del trabajo publicado por [24], mientras que, el equivalente de red en la barra la Niña se obtuvo en base a la información publicada por el Comité de Operación Económica del Sistema Interconectado- Perú (COES) [25]. En la Tabla 4 se muestran los parámetros obtenidos para la red equivalente en la barra de Chorrillos en Ecuador y en la barra de La Niña en Perú.

Tabla 4: Equivalente de Red en la barra Chorrillos y la Niña 

Al realizar la implementación y configuración de los elementos del sistema de estudio mediante el software ATPDraw se obtiene el sistema mostrado en la Figura 10.

Figura 10: Interconexión Ecuador-Perú representado en software ATPDraw 

RESULTADOS

Para determinar las sobretensiones representativas por maniobra a través de un análisis estadístico, se someterá el sistema de estudio a tres escenarios de energización, en los cuales se verificarán las sobretensiones 𝑉 𝑒2 (fase-tierra) y 𝑉 𝑝2 (fase-fase):

Escenario 1: Energización de líneas de transmisión sin ningún control de voltaje

Escenario 2: Energización de líneas de transmisión con bancos de reactores conectados a la línea.

Escenario 3: Energización de líneas de transmisión con resistencias de pre-inserción de 400 ohm.

Luego, aplicando los factores de conversión (ver Sección 3) se determinan las sobretensiones de truncamiento 𝑉 𝑒𝑡 (fase-tierra) y 𝑉 𝑝𝑡 (fase-fase). En cada escenario se compararán los valores se tensión obtenidos mediante la distribución de Gauss con la distribución de Weibull modificada.

Escenario 1

En la Tabla 5 y Tabla 6 se muestran las sobretensiones fase-tierra y entre fases al energizar las líneas sin ningún tipo de control de voltaje, comparando los resultados entre la distribución Gaussiana y Weibull. La diferencia entre ambas distribuciones para el voltaje estadístico Ve2 de la L/T CHR-PSJ es de 2.8% y en la L/T PSJ-PRA es de 5.3% en sobretensiones fase-tierra.

Tabla 5: Resultados Escenario 1 de sobretensiones fase-tierra Ve2 y Vet 

Tabla 6: Resultados Escenario 1 de sobretensiones fase-fase Vp2 y Vpt 

Las Figura 11 y Figura 12 comparan las funciones de distribución Gaussiana y de Weibull para el Escenario 1 en la línea de transmisión CHR-PSJ y PSJ-PARA respectivamente.

Figura 11: Distribución Estadística Escenario 1 CHR-PSJ 

Figura 12: Distribución Estadística Escenario 1 PSJ-PRA 

Escenario 2

En la Tabla 7 y Tabla 8 se muestran los resultados de sobretensiones utilizando al banco de reactores con los valores establecidos en el PET, junto con los resultados de comparación de distribución estadística Gaussiana y Weibull. En este caso, la variación entre ambos métodos fue de 2% para la L/T CHR-PSJ y 5% en la L/T PSJ-PRA.

Tabla 7:/

Tabla 8: Resultados Escenario 2 de sobretensiones fase-fase Vp2 y Vpt

De manera análoga al Escenario 1, En la Figura 13 y la Figura 14 se comparan las funciones de distribución Gaussiana y de Weibull para el Escenario 2 de las dos líneas de transmisión en análisis.

Figura 13: Distribución Estadística Escenario 2 CHR-PSJ 

Figura 14: Distribución Estadística Escenario 2 PSJ-PRA 

Escenario 3

Para la evaluación del escenario 3 se utilizaron resistencias de pre-inserción que por lo general en los SEP pueden tomar valores de 300, 400 y 500 ohm, en este caso se ha seleccionado 400 ohm. De manera análoga a los dos casos anteriores, en la Tabla 9 y Tabla 10 se presentan variaciones en la distribución estadística de Weibull con respecto a Gauss en el orden de 2.5% y 3.1% para la L/T CHR-PSJ y L/T PSJ-PRA respectivamente.

Tabla 9: Resultados Escenario 3 de sobretensiones fase-tierra Ve2 y Vet

Tabla 10: Resultados Escenario 3 de sobretensiones fase-fase Vp2 y Vpt

Para el escenario 3, los resultados de la distribución Gaussiana y de Weibull están representados en la Figura 15 y Figura 16.

Figura 15: Distribución Estadística Escenario 3 CHR-PSJ 

Figura 16: Distribución Estadística Escenario 3 PSJ-PRA 

Finalmente, al obtener los valores de sobretensiones Ve2 y Vp2 aplicando las distribuciones estadísticas indicadas, y al ser evaluadas mediante su función completaría, se calcula la diferencia entre los resultados arrojados por la distribución Gaussiana y Weibull modificada mostrados en la Tabla 11. Donde, la variación viene dada en porcentaje y toma como referencia la distribución Gaussiana.

Tabla 11: Resultados de variación de sobretensiones

Al analizar la variación promedio de los voltajes Ve2 y Vp2 de los tres escenarios en las líneas de trasmisión CHR-PSJ y L/T PSJ-PRA, se evidencia que la diferencia es mayor para voltajes de fase-fase.

CONCLUSIONES

Se evidencia que la diferencia que existe entre ambas distribuciones es mínima. La dispersión de datos es menor a 3% en ambos casos, por lo tanto, ambas distribuciones son aplicables para un análisis estadístico de sobretensiones si se consideran magnitudes de voltaje con igual probabilidad de aparecer en el sistema, siendo recomendable por la IEC 60071-2 un 98% de probabilidad.

Las estrategias analizadas en este documento para la reducción de la amplitud de sobretensiones por energización contribuyen a amortiguar la elevación de tensión entre un 35.00 a 40.00 % con respecto al escenario 1. Siendo la mejor opción para controlar el pico de sobretensiones la utilización de resistencias de pre inserción, por otro lado, los bancos de reactores controlan de mejor manera las sobretensiones temporales después de la energización.

Las sobretensiones más elevadas se dieron en el escenario 1 de energización sin control de voltaje con una tensión de 2.688 p.u. con distribución Gaussiana y 2.635 p.u. con distribución Weibull. Estos resultados comprueban la necesidad de controles de voltaje antes de la energización y definen a la distribución de Weibull ligeramente más conservadora que la distribución Gaussiana.

REFERENCIAS

[1] CELEC TRANSELECTRIC, “Plan de Expansión de la Transmisión”, Quito, 2018 [ Links ]

[2] J. Martínez, Coordinación de Aislamiento en Redes Eléctricas de Alta Tensión, I. Madrid, MC Graw Hill, 2008. [ Links ]

[3] M. Benesz, W. Nowak, W. Szpyra, and R. Tarko, “Application of statistical methods in insulation coordination of Overhead Power Lines”, in Departmente of Electrical Enginneering and Power Engineering, 2017, pp. 2-5, doi: 10.1109/EPE.2017.7967316. [ Links ]

[4] A. H. Hamza, S. M. Ghania, A. M. Emam, and A. S. Shafy, “Statistical analysis of switching overvoltages and insulation coordination for a 500 kV transmission line,” in 2016 18th International Middle-East Power Systems Conference MEPCON 2016 - Proceedings ,2015, no. July, pp. 4-5, doi: 10.1109/MEPCON.2016.7836966. [ Links ]

[5] IEEE Std 1234, Guide for Improving the Lightning Performance of Transmission Lines. New York, 1997. [ Links ]

[6] International Standard IEC 60071-2, Insultation co-ordination part 2: Application guide, vol. I, no. Coordinacion de aislamiento Parte 2. 1997. [ Links ]

[7] A. R. Hileman, Insulation Coordination for Power Systems, vol. 19, no. 9. New York: Taylor and Francis, 1999. [ Links ]

[8] M. Sumair, T. Aized, S. A. R. Gardezi, M. M. A. Bhutta, S. M. S. Rehman, and S. U. ur Rehman, “Application of five continuous distributions and evaluation of wind potential at five stations using normal distribution,” Energy Explor. Exploit., vol. 39, no. 6, pp. 2214-2239, 2021, doi: 10.1177/0144598720939373. [ Links ]

[9] E. Chaj, “Análisis De Sobretensiones Debido a Transitorios Por Maniobras En Sistemas Eléctricos De Potencia Mayores De 300 Kv,” Coast. Estuar. Process. , pp. 1-360, 2009. [ Links ]

[10] L. Jaimis and INEL, “Ondas viajeras en Lineas de Transmision,” Ondas Viajeras, 2020. https://inelinc.com/ondas-viajeras-en-lineas-de-transmision-avhqn. [ Links ]

[11] K. Sailema, X. Proaño , and L. Ruales , “Analisis de Sobrevoltajes a Frecuencia industrial y por Maniobra en Lineas de Transmision a Nivel de 500 kV,” Esc. Politec. Nac., vol. 0, no. 0, pp. 1-10, 2010 ;Available: http://bibdigital.epn.edu.ec/handle/15000/2479.Links ]

[12] Power Engineering Society IEEE 1313, “Guide for the Application of Insulation Coordination,” New York, 1999. [ Links ]

[13] International Standard IEC 60071-1, Insulation coordination. Part 1: Definitions, principles and rules. New York, 2006. [ Links ]

[14] N. Morales, “Coordinación de aislamiento en Subestaciones a Nivel de 500 kV.,” Escuela Politecnica Nacional, 2008. [ Links ]

[15] G. Canavos, Applied Probability and Stadistical Methods, VII. Mexico: Camara Nacional de La industria Editorial, 2005. [ Links ]

[16] R. Collaguazo, “Modelación de Transitorios Electromagnéticos en Sistemas de Transmisión por medio del Modelo J. Martí,” Univ. Politécnica Sales., vol. 0, no. 0, p. 20, 2017. [ Links ]

[17] L. H-Restrepo Gladys Caicedo Delgado Ferley Castro-Aranda, “Transmission line models for electromagnetic transients in power systems,” Rev. Energía y Comput., vol. 16, no. 1, pp. 14-26, 2016. [ Links ]

[18] J. F.Rios, “Coordinacion de aislamiento para la linea a 500kV , doble circuito , de la interconexion Ecuador-Peru por sobrevoltajes de maniobra.”, Escuela Politecnica Nacional, Quito, 2019. [ Links ]

[19] G. Uquillas, “Plan de Desarrollo Electrico y Oportunidades de Inversion,” 2020, p. 34, [Online]. Available: https://issuu.com/artesjournal/docs/presentacion_ing._gonzalo_uquillas._plan_de_desarr.Links ]

[20] COES SINAC and CELEC EP, “Servicio de Consultoria para el Desarrollo del Anteproyecto Interconexion en 500 kV Ecuador-Perú,” 2014. [ Links ]

[21] Y. Suarez, “Simulacion de Tensiones Transitorias de Recuperacion Mediante ATPDraw.” Universidad Tecnologica Nacional, Santiago de Chile, p. 38, 2015. [ Links ]

[22] COES SINAC, “Indice general enlace 500 kV la Niña - Piura, subestaciones, líneas y ampliaciones asociadas,” Lima, 2016. [ Links ]

[23] G. Mejia, ATPDraw Guia Basica de uso. Antioquia, slideshare, 2016. [ Links ]

[24] D. Reinoso, “Desarrollo de Casos de Estudio en Base al Sistema Nacional Interconectado de Ecuador para el Análisis de Sistemas Eléctricos de Potencia.,” Universidad Técnica de Cotopaxi, 2022. [ Links ]

[25] COES SINAC, “Comité de Operación Económica del Sistema Interconectado,” Diagrama Unifilar SEIN, 2021. https://www.coes.org.pe/Portal/Operacion/CaractSEIN/DiagramaUnifilar. [ Links ]

Recibido: 06 de Noviembre de 2022; Aprobado: 13 de Enero de 2023

John Henry Calle. - Nació en Ambato, Ecuador en 1998. Recibió su título bachiller en físico-matemático de la U.E “Bolívar” en 2017; de Ingeniero Eléctrico en Sistemas Eléctricos de Potencia de la Universidad Técnica de Cotopaxi en 2022. Sus campos de investigación están relacionados con los Sistemas Eléctricos de Potencia y Diseño en Alto Voltaje.

Wilian Guamán Cuenca. - Nació en Riobamba en 1989. Recibió su título de Ingeniero en Electromecánica por la Universidad de las Fuerzas Armadas ESPE en 2013 y de máster en Ingeniería de la Energía por la Universidad Politécnica de Madrid en 2017. Su campo de investigación se encuentra relacionado con la Planificación de los Sistemas de Eléctricos de Potencia.

Creative Commons License Este es un artículo publicado en acceso abierto bajo una licencia Creative Commons