Análisis del fenómeno de voltaje transitorio de restablecimiento (TVR),
Omar Montufar Valer1 menor tiempo mediante simulación.
Las herramientas necesarias para esta vía
1. National University of San Antonio
de cálculo son las Leyes de Kirchoff y la
Abad of Cusco, Peru
Transformada de Laplace.
Abstract: Este trabajo presenta el 2.1.1. Leyes de Kirchoff
análisis de resultados de TRV, Para calcular corrientes y voltajes en una
simulación de falla kilométrica y red equivalente, se utilizan las leyes de
ferroresonancia por método Kirchoff, y con ellas se plantean
computacional usando simulador ecuaciones en las que intervienen
Digsilent. estas variables del circuito, generando
ecuaciones diferenciales.
En base a los resultados se verifica la
Ley 1: La suma de todas las corrientes
curva característica de soportabilidad
que entran a cualquier unión (o nodo)
para el sistema de prueba.
de un circuito es igual a la suma de las
corrientes que salen de dicha unión.
Key words: TRV, transitorio, análisis,
Ley 2: En una espira, la suma de todas
circuito, voltaje, método, Kirchoff,
las caídas de voltaje es igual al voltaje
Laplace.
total suministrado. Se sabe que los
voltajes o las diferencias de potencial a
1. Introduction través de los inductores (L), resistores
(R) o capacitores (C) están dados por las
Un transitorio es un fenómeno de alta siguientes fórmulas:
frecuencia y elevados niveles de corriente y 𝑉𝑅 = 𝑅𝑖 (𝑡)
tensión, son capaces de modificar la
topología de la red si no se controlan. Este 𝑖 (𝑡 )
𝑉𝐿 = 𝐿 ⅆ
fenómeno está presente en la red ⅆ𝑡
eléctrica por cada maniobra que se realice. 1
𝑉𝑐 = ∫ 𝑖 (𝑡) ⅆ𝑡
Las maniobras de apertura y cierre de 𝑐
interruptores de potencia por ejemplo
2.1.1. Transformada de Laplace
inyectan transitorios llamados transitorios
de voltaje de restablecimiento cuyo estudio Con esta herramienta podremos reducir las
y cálculo es bastante complejo debido a ecuaciones diferenciales a simples
que depende de múltiples factores. expresiones algebraicas de sencilla
resolución, donde las condiciones iniciales
2. Materials and Methods de éstas son previamente conocidas.
Definimos la transformada de Laplace de una
2.1 Método Analítico
función f(t) de la siguiente manera.
Este método requiere la resolución de
𝐹(𝑠) = 𝐿{𝑓 (𝑡)} = ∫ ⅇ −𝑠𝑡 𝑓(𝑡) ⅆ𝑡
procesos matemáticos complejos para una
situación en particular del sistema, esto
debido a las características y el número
elevado de elementos que los componen;
por lo que se utiliza programas
computacionales para el cálculo de
transitorios electromagnéticos que
proporcionan resultados satisfactorios y en
Donde “s” es una variable compleja y ⅇ −𝑠𝑡
es llamdo el núcleo de la transformación, Donde:
es decir la transformada de Laplace toma
una función f(t) en el dominio del tiempo
definida para 𝑡 ≥ 0 y le asocia una variable
compleja F(S) en el dominio de S llamada
dominio de frecuencia compleja.
La transformada de Laplace de f(t) existe
si la integral impropia converge.
El proceso inverso se llama INVERSA de
Laplace, que transforma una función de 𝐼𝑎 = 𝐼0 + 𝐼1 + 𝐼 2 𝑉𝐵 𝑉 = 𝑉𝐶 = 0
variable real f(t), donde 𝑡 ≥ 0. De acuerdo a la definición de componentes
Al resolver el problema en el dominio S, simétricas, los voltajes: V0, V1y V2, se
aplicaremos la transformada inversa de presentan a continuación en función del
Laplace para volver al dominio t, en este operador que es igual a:
caso la transformada de F(S) es la función 1
𝑉0 = (𝑣𝑎 + 𝑣𝑏 + 𝑣𝑐 )
3
𝑓(𝑡) = 𝐿−1 {𝐹 (𝑠)}
1
𝑉1 = (𝑉𝑎 + 𝑎 𝑉𝑏 + 𝑎2 𝑉𝑐 )
2.2. Componentes simétricas 3
Este método de cálculo reemplaza el circuito
con la asimetría de una falla por tres circuitos 1
balanceados para resolver el transitorio en 𝑉2 = (𝑣𝑎 + 𝑎2 𝑣𝑏 + 𝑎𝑉𝑐 )
3
componentes de secuencia positiva,
negativa y cero y que puede ser aplicado
en más de una forma. Reemplazando obtenemos la nueva
ecuación:
A continuación se presenta la aplicación de
este método para fallas trifásicas en un 𝑉0 = 𝑉1 + 𝑉2 + 𝑉3
sistema, el cual presenta a los voltajes 𝑣𝐴
entre los contactos del interruptor muy 𝐼0 + 𝐼1 + 𝐼2 =
𝑧
asimétricos cuando el primer polo del
interruptor es abierto y los otros aún
siguen cerrados como se muestra en la Las impedancias de secuencias designadas
Figura 2 en donde los voltajes y corrientes por Z0, Z1 y Z2 corresponde a las vistas en
del interruptor son designados como VA, el sistema equilibrado a través de los
VB, VC , IA, IB e IC en las fases contactos del interruptor en cualquiera de las
fases, por lo que en términos de las
respectivas y cuyoscomponentes
impedancias de secuencia reescribiendo las
simétricas son V0, V1 y V2, e I0, I1 y I2.
ecuaciones anteriores se tiene:
El polo abierto es simulado por la
introducción de una impedancia 𝑉0 = 𝐼0 𝑧0 𝑈 ′ 𝑣1 = 𝐼1 𝑍1 + 𝐸1 𝑈𝑉2 = 𝐼2 𝑍2
infinita Z entre los contactos del 𝑉𝑎 = 3𝑉𝑜 = 3𝑉1 = 3𝑉2
interruptor (VA=IAZ) cuya
representación y conexión en redes 3𝑉𝑜 = 𝐼𝑜 3𝑍𝑜 𝑈 3𝑉1 = 𝐼1 3𝑍1 +
de secuencia del interruptor se 3 𝐸1 𝑈 3 𝑉2 = 𝐼2 3𝑍2 + 3𝐸 1
representa en la Figura 1.
Donde la corriente en la fase abierta y
el voltaje que aparece a través del
circuito abierto en el primer polo del
interruptor corresponde al TRV, el cual
resulta de la suma de las respuestas
de las redes de secuencia y que se
interconectan de acuerdo para dar
valores verdaderos de respuesta del
sistema.
Sin embargo, no es recomendable para
fallas asimétricas ya que no está bien
adaptado para este tipo de fallas, debido a
Se ha simulado el sistema corrido flujo de potencia
que, visto desde los contactos ya no se
desbalanceado y condiciones iniciales produciéndose
busca un sistema equilibrado como nen una
falla trifásica y estos casos se trata de un corto circuito en la barra 03, cortocircuito trifásico,
resolver con cálculos de cortocircuito de en un tiempo de t=0 segundos, ocasionando un flujo
componentes simétricas. de corriente direccionada a tierra denominada
corriente de cortocircuito, ocasionando una apertura
de CB el cual se ha programado 100 milisegundos en
3.- SIMULACIONES caso ocurra un cortocircuito. Con el objetivo de aislar
3.1 SIMULACION Y ANALISIS FALLA EN la falla en la zona donde ha ocurrido.
TERMINALES Observamos el comportamiento de la barra 3.
Esta simulación se realiza en un SEP con
En la grafica observamos el comportamiento de la
las siguientes características
tensión en la barra 3, observamos un comportamiento
Red Slack: Smax = 8000 MVA; normal de la tensión en las tres fases, la cual esta
Smin = 2000 MVA representada por las graficas de color azul verde y
rojo y representa una onda senoidal constante. Al
T = 8 seg. momento de ocurrir el cortocircuito el nivel de
V1 = 1< 0º pu tensión en esta barra se hace cero (Vb=0). Este
Red PV : Smax = 8000 MVA; cambio de nivel de voltaje ocurre de manera
intempestiva, en la grafica observamos que al evaluar
Smin = 5000 MVA en el tiempo t=0, hay una variación de tensiones en
T = 8 seg. las tres fases en el momento de la apertura del
circuito
/V5/ = 1< 0º pu
PG5 = 100 MW La tensión de la fase A (color Azul) se encuentra en
𝛺
183.40 kv, la tensión de la fase B (color rojo) se
Linea 1-2: Z12 = 0.02+0.2j[𝑘 ] encuentra en 145.75 kv y la tensión de la fase C ( color
𝑚
𝛺
verde) se encuentra en 37.65 kv. Esto debido al
Z𝑜 = 0.3 + 1.1𝐽 𝑘𝑚 𝐿 = 100𝑘𝑚 desfasaje natural de 120 grados que tiene cada fase.
𝜇𝑠 𝜇𝑠
B12 = 2.8 𝐵0 = 1.8 𝑘𝑚
𝑘𝑚
𝐼 = 0.8𝑘𝐴 modelo II
𝑇𝑅 − 1: 𝑉𝑐𝑐 = 10%
𝑉 = 138⁄13.8𝑘𝑣 𝑐𝑜𝑛ⅇ𝑥𝑖𝑜𝑛 𝑌𝐷
𝑆 = 50𝑀𝑉𝐴
Carga 01: 𝑝 = 15𝑀𝑤
𝑓 𝑑 𝑝 = 0.96 𝑖𝑛ⅆ𝑢𝑐𝑡𝑖𝑣𝑜
Analizaremos ahora el comportamiento de las
corrientes, en la grafica N.º 4, se observan 3 niveles de
comportamiento, las graficas representan las
corrientes antes durante y después de ocurrida la falla
de cortocircuito.
Fundamentalmente el nivel de corriente a las tres
fases del sistema presenta un incremento siendo
Icc>>In, lo cual se observa claramente en la gráfica 4.
Analizando el comportamiento de las ondas antes del
cortocircuito, podemos observar que las corrientes de
las tres fases tienen forma senoidal dentro de sus
rangos de operación del sistema. Cuando en t=0
instante donde se presenta el cortocircuito, las
corrientes se incrementan de manera intempestiva,
donde sus amplitudes son muyo mayores a la
Observamos el comportamiento de la barra 4.
nominal. Ocasionando la presencia de componentes
El fenómeno que ocurre en esta barra es muy de corriente continua y sobre corrientes en el SEP, las
peligroso, la cual pasaremos a detallas a continuación: corrientes máximas de cortocircuito debido a la
Cuando el sistema esta se condiciones normales de
funcionamiento, el comportamiento de las ondas de
tensión en las tres fases es de naturaleza senoidal,
cuando ocurre un cortocircuito para un tiempo t=0
significa que la tensión en la barra numero 4 también
será cero la cual permanecerá en ese valor por un
periodo de tiempo donde dura el cortocircuito, para
t=0,1 segundo tiempo en el que demora la apertura
de interruptor Figura 2), las tensiones sufren un
incremento una diferencia de potencial en la barra 4,
este nivel de tensión se produce junto a un conjunto
de oscilaciones debido a las contaminaciones
componente continua empezarán a extinguirse
producidas por el evento de cortocircuito, estas
tensiones oscilan a velocidades angular más rápida
que la velocidad angular del sistema. Nuestro sistema
tiene una barra PV, y debido a la robustes de nuestro En la siguiente grafica podemos observar que en
sistema origina que estas ondas empiezan a interrupto deja pasar un poco de corriente de
amortiguar, en un punto en el cual la onda será cortocircuito antes de su apertura, esto se logra
nuevamente senoidal pura. El fenómeno de realizando una buena calibración y coordinación.
contaminación se ha presentado de manera Se observa en la simulación que el interruptor se
momentánea la cual empieza a amortiguarse. Figura 3 apertura en 100 milisegundos dejando pasar la
corriente hasta un instante en que son mínimas. Solo
en ese instante se debe apertura el interruptor CB
Analizando más a detalle el comportamiento de las
corrientes y el CB, en la siguiente grafica podemos
observar que al interruptor deja pasar un poco de
corriente de cortocircuito con el fin de evitar la
presencia de sobretensiones conocido como
CLOUSING, el interruptor debe buscar el punto de
apertura cuando las corrientes sean mínimas.
Si se apertura en un instante donde las corrientes no
son mínimas pueden presentarse eventos muy
peligrosos llamados arcos eléctricos.
Análisis de TRV, en la siguiente grafica observamos
que la sobre tensión no debe ser mayo a 2Vm , debido
a la robustes no se ha producido un TRV considerable,
si se observa que la tensión se incrementa,
semejándose a una onda triangular donde se puede
calcular el transitorio de recuperación de voltaje
RRRV, se observa que se ha producido oscilaciones
mucho más rápida s las oscilaciones de sistema
(Wo>>Wn), , estos dos parámetros darán origen a una
frecuencia natural la cual será mucho mayor a la
frecuencia de operación.