1 Estabilidad de SEP
1 Estabilidad de SEP
1 Estabilidad de SEP
potencia
• Capacidad de
recuperar su
estado de
funcionamiento
después de haber
sido sometido a
una perturbación.
23/04/2016 cchilet@continental.edu.pe 2
Estabilidad en sistemas de potencia
• Estabilidad de ángulo: Es la
capacidad del sistema de
mantener a los generadores en
sincronismo después de una
perturbación.
• Estabilidad de tensiones
(colapso de tensión):
capacidad de mantener las
tensiones en los nudos
próximos al 1,0 p.u. tras una
perturbación.
23/04/2016 cchilet@continental.edu.pe 3
Tipos
La estabilidad de
ángulo se clasifica en:
– Estabilidad de régimen
permanente.
– Estabilidad transitoria.
– Estabilidad dinámica.
23/04/2016 cchilet@continental.edu.pe 4
Estabilidad de régimen permanente
• Estabilidad estática.
• Capacidad de recobrar el sincronismo tras
perturbaciones relativamente pequeñas y
graduales.
• Ejm. Potencia demandada por las cargas.
23/04/2016 cchilet@continental.edu.pe 5
Estabilidad transitoria (TS)
• Capacidad de recobrar el sincronismo tras
una gran perturbación.
• Las máquinas síncronas experimentan
desviaciones respecto a su valor de
sincronismo (60 Hz), cambiando sus ángulos
de potencia.
• Objetivo del Estudio TS: determinar si las
máquinas recuperan el sincronismo.
23/04/2016 cchilet@continental.edu.pe 6
Estudio de TS
23/04/2016 cchilet@continental.edu.pe 7
Estabilidad transitoria
• Viene determinada por la evolución del
ángulo de potencia.
• Condiciones:
– Duración de la primera oscilación=1s.
– Generadores operan a la Ns.
– U y I son AC sinusoidales.
– Potencia mecánica suministrada = constante.
– Iexc = constante.
– Carga constante
23/04/2016 cchilet@continental.edu.pe 8
Estabilidad dinámica
23/04/2016 cchilet@continental.edu.pe 9
Secuencia de eventos
• Fase 1: Prefallo
• Fase 2: Fallo y desconexión de la línea.
• Fase 3: Post fallo.
¿Luego de una secuencia de eventos dada, recupera
o nó el sistema la estabilidad?
23/04/2016 cchilet@continental.edu.pe 10
Equilibrio dinámico del grupo T-G
• Condición inicial:
δ= fijo
• En perturbación: se
presenta movimiento
relativo ⇒ δ cambia (rotor
se acelera o desacelera)
• Ecuación que describe
movimiento relativo
(ecuación de oscilación)
23/04/2016 cchilet@continental.edu.pe 11
Equilibrio dinámico del grupo T-G
Considerando:
• Que los rotores del grupo T-G forman un
sólido rígido
• Despreciable:
– Los pares debido al rozamiento y fricción de las
partes mecánicas y al devando amortiguador.
– Las pérdidas eléctricas.
23/04/2016 cchilet@continental.edu.pe 13
Introducción
23/04/2016 cchilet@continental.edu.pe 14
Diagrama unifilar Pre-falla en p.u.
• La tensión interna del generador durante el
régimen transitorio será igual a la que existía
antes de la perturbación y se denotará como
E´.
23/04/2016 cchilet@continental.edu.pe 15
Premisas
• Para facilitar el análisis se despreciarán:
– las pérdidas en el sistema de transporte y
generación,
– el efecto de la saturación del circuito magnético
del generador y el de los polos salientes.
– Durante el tiempo que dura la perturbación la
Iexc de la máquina es CTE, sin embargo, la NGEN
se desviará de su Ns y el ángulo δ variará
siguiendo el modelo definido por la ecuación de
oscilación del generador .
23/04/2016 cchilet@continental.edu.pe 16
Pre-falla
• Se calcula la Xeq del circuito en el régimen
previo a la perturbación (XT), como:
23/04/2016 cchilet@continental.edu.pe 17
Pre-falla
• El punto de
funcionamiento
de la máquina
corresponderá a
un ángulo δ
determinado, tal
que PI = Pm.
• Este ángulo es el
δI
23/04/2016 cchilet@continental.edu.pe 18
Ocurrencia de la falla
23/04/2016 cchilet@continental.edu.pe 19
En falla
• Como consecuencia de un cortocircuito en
el extremo de la LT2 las protecciones dan la
orden a los interruptores automáticos de
dejar abierta esta línea.
23/04/2016 cchilet@continental.edu.pe 20
En falla
• La inductancia total equivalente una vez
despejado el defecto tomará otro valor,
lógicamente mayor:
23/04/2016 cchilet@continental.edu.pe 21
En falla
• Por tanto, en Falla (II), la máxima potencia
eléctrica que podrá entregar el generador
que será inferior seguirá la expresión
siguiente:
23/04/2016 cchilet@continental.edu.pe 22
En falla
• Como la Pm de la
turbina = CTE, el punto
de funcionamiento de
la máquina pasará, una
vez transcurrido el
transitorio y siempre
que el sistema
permanezca estable, a
un ángulo mayor (δII).
23/04/2016 cchilet@continental.edu.pe 23
Ecuación de oscilación
23/04/2016 cchilet@continental.edu.pe 24
Ecuación de oscilación
• Durante cualquier fenómeno transitorio el par de
aceleración de la máquina en su conjunto (turbina
más generador) es igual al par motor mecánico de la
turbina menos el par electromagnético o de frenado
(supónganse medidos ambos pares en unidades de
N.m), es decir:
23/04/2016 cchilet@continental.edu.pe 25
Ecuación de oscilación
• En régimen permanente Mm = Me, es decir, el
par de aceleración será nulo, y el rotor gira a
una velocidad constante que es la velocidad
de sincronismo.
• Si el par de aceleración es positivo la
velocidad del rotor aumentará y si el par de
aceleración es negativo disminuirá.
23/04/2016 cchilet@continental.edu.pe 26
Ecuación de oscilación
• La variación de la velocidad de giro de la
máquina dependerá de su momento de
inercia, según la expresión siguiente:
donde:
J Momento de inercia combinado de la máquina
(generador y turbina) en kg.m2.
ωm velocidad angular mecánica del rotor, siendo su
valor de sincronismo ω0m.
23/04/2016 27
Ecuación de oscilación
• La ecuación (2.5) se puede transformar en
otra equivalente empleando la constante de
inercia normalizada de la máquina, H, que se
define como:
donde:
H constante de inercia normalizada de la máquina en segundos.
ω0m velocidad angular mecánica de sincronismo del rotor (rad/s)
SN potencia trifásica nominal del generador (VA).
23/04/2016 cchilet@continental.edu.pe 28
Constante de inercia normalizada (H)
• En la tabla siguiente se muestran valores típicos de la
constante de inercia normalizada en función del tipo
de máquina:
23/04/2016 cchilet@continental.edu.pe 29
Ecuación de oscilación
• Las expresiones 2.6 y 2.7 anteriores se pueden
combinar y transformar en las siguientes:
23/04/2016 cchilet@continental.edu.pe 30
Ecuación de oscilación
• Si se define momentos y velocidades
angulares en valores por unidad con las bases
siguientes:
23/04/2016 cchilet@continental.edu.pe 31
Ecuación de oscilación
• Se puede transformar la expresión anterior en la
siguiente:
donde:
• ω velocidad angular eléctrica del rotor en rad/s
• ω0 velocidad angular eléctrica inicial del rotor en
rad/s
• ω(pu) velocidad angular eléctrica del rotor en por
unidad
• p número de pares de polos.
23/04/2016 cchilet@continental.edu.pe 32
Ecuación de oscilación
• δ →la posición angular del rotor en rad.
eléctricos respecto de un eje de referencia y
δ0 → es su posición en un tiempo t = 0, se
tiene que:
23/04/2016 cchilet@continental.edu.pe 33
Ecuación de oscilación
23/04/2016 34
• De no existir amortiguamiento el sistema
seguiría oscilando continuamente.
• No obstante como para realizar el Estudio de
estabilidad transitoria es suficiente con
estudiar hasta la primera oscilación (t<1 s).
• Se puede prescindir del efecto del
amortiguamiento posterior.
23/04/2016 cchilet@continental.edu.pe 35
• Durante la perturbación, la oscilación del
ángulo δ, se superpone a la velocidad de
sincronismo ω0.
• No obstante, como la variación de la
velocidad (Δω = dδ/dt) es mucho más
pequeña que la velocidad inicial ω0, se puede
considerar que la velocidad permanece
prácticamente constante e igual a ω0.
23/04/2016 cchilet@continental.edu.pe 36
• En consecuencia, en valores por unidad, el
momento mecánico y el electromagnético
son prácticamente iguales a las potencias
mecánica y eléctrica respectivamente
(recuérdese que la potencia es igual al par
multiplicado por la velocidad angular).
23/04/2016 cchilet@continental.edu.pe 37
Ecuación de oscilación de potencia
• Con tal simplificación se tiene que la ecuación de la
oscilación se puede escribir finalmente como:
2 H d 2δm(t )
⋅ = Pm( PU ) (t ) − Pe( PU ) (t )
ωm 0 dt 2
2 H d δm(t )2
⋅ = Pm(t ) − Pe(t )
ωm 0 dt 2
23/04/2016 cchilet@continental.edu.pe 39
Criterio de igualdad de áreas
23/04/2016 cchilet@continental.edu.pe 40
Ecuación de oscilación
• La ecuación de oscilación de la máquina síncrona es:
Donde:
• H es la constante de inercia,
• ω0 la frecuencia de
sincronismo,
• Pm la potencia mecánica
Ecuación de oscilación
entrante y
despreciando torque amortiguante • Pe la potencia eléctrica
saliente.
Ecuación diferencial que describe el movimiento relativo del eje del rotor
respecto de los ejes del campo magnético resultante durante una
perturbación.
23/04/2016 cchilet@continental.edu.pe 41
Fundamento matemático
criterio de igualdad de áreas
23/04/2016 cchilet@continental.edu.pe 42
Criterio de igualdad de áreas
• En el instante en que el
• Después de la ángulo alcance su
perturbación el sistema máximo valor la
será estable si el pendiente de la curva
ángulo δ alcanza de variación del ángulo
rápidamente (en será nula, por tanto:
menos de 1 segundo)
un valor máximo (δm)
para disminuir a
continuación.
23/04/2016 cchilet@continental.edu.pe 43
Criterio de igualdad de áreas
• Por lo tanto la condición de estabilidad se
puede definir como:
23/04/2016 cchilet@continental.edu.pe 44
Criterio de igualdad de áreas
• Es decir que para que el
sistema sea estable el
área de la curva (Pm-Pe)
representada en función
del ángulo δ, debe ser
cero:
23/04/2016 cchilet@continental.edu.pe 45
Criterio de igualdad de áreas
• Si se separa el área positiva (A1) o de aceleración
(cuando Pm > Pe), del área negativa (A2) o de
deceleración (cuando Pm < Pe) se llega al criterio
de igualdad de áreas:
23/04/2016 cchilet@continental.edu.pe 46
Criterio de igualdad de áreas
23/04/2016 cchilet@continental.edu.pe 47
Criterio de igualdad de áreas
• El ángulo δ1 representa el ángulo del rotor en el
instante en que tras una perturbación se
reconfigura el sistema, por ejemplo se abren
interruptores para aislar la línea en defecto, con lo
que la potencia activa entregada por el generador
(Pe) que generalmente ha sido pequeña o incluso
nula durante la perturbación o cortocircuito, pasa
bruscamente a ser mayor que la mecánica y el
generador empieza a frenarse.
23/04/2016 cchilet@continental.edu.pe 48
Criterio de igualdad de áreas
• El mayor valor del ángulo δ1 que satisface la
condición de estabilidad se nomina ángulo crítico, y
el tiempo de actuación de las protecciones en caso
de cortocircuito o el tiempo de re conexión de una
línea en defecto que conduce a obtener dicho valor
de ángulo se denomina tiempo crítico.
• Si las protecciones tardan en actuar un tiempo
superior al crítico el área de aceleración es tan
grande que no es posible conseguir un área de
deceleración igual a la de aceleración y el sistema se
vuelve inestable.
23/04/2016 cchilet@continental.edu.pe 49
Ejemplo 1
• Un grupo generador de 60 Hz y H =5
segundos, se conecta a través de un
transformador de impedancia, Xt = j0,1 pu y
de una línea de impedancia XL = j0,3 a un
nudo de potencia infinita de tensión U=1p.u.
que consume una potencia activa de 0,8 pu
con un factor de potencia 0,7 inductivo.
• La reactancia transitoria del generador es X´d
= j 0,25 pu.
23/04/2016 cchilet@continental.edu.pe 50
Ejemplo 1
• En un instante determinado, t=0, se produce un
cortocircuito trifásico a tierra debido a la caída de
unas ramas de árbol sobre la línea. Las protecciones
actúan de forma instantánea disparando la línea, de
forma que el cortocircuito ha desaparecido cuando
transcurrido el tiempo necesario para el reenganche
o re conexión se vuelven a cerrar los interruptores
que ponen en servicio la línea. Calcular el ángulo
crítico del rotor y el tiempo crítico de re conexión
que garantizan la estabilidad transitoria del sistema.
23/04/2016 cchilet@continental.edu.pe 51
Ejemplo 1
23/04/2016 cchilet@continental.edu.pe 52
Ejemplo 1
23/04/2016 cchilet@continental.edu.pe 53
Solución 1
La potencia activa Teniendo en cuenta que:
entregada por el
generador al nudo de
potencia infinita en
función del ángulo del
rotor sigue la siguiente
expresión:
23/04/2016 cchilet@continental.edu.pe 54
Solución 1
Se debe determinar la fuerza electromotriz interna de la
máquina. Su valor permanecerá constante, ya que la excitación
permanece constante. Para ello se parte del dato de la potencia
activa y reactiva consumida por la carga en el nudo de potencia
infinita.
23/04/2016 cchilet@continental.edu.pe 55
Solución 1
Por otra parte la fuerza electromotriz interna (E) de la
máquina se calcula teniendo en cuenta el valor de la
reactancia serie total.
23/04/2016 cchilet@continental.edu.pe 56
Solución 1
A continuación se calcularán: el ángulo crítico (δc) del rotor y el
tiempo crítico de reconexión.
Como paso previo se calcula el ángulo inicial del rotor antes de la
falta:
23/04/2016 cchilet@continental.edu.pe 57
Solución 1
La potencia mecánica
permanece constante e igual a
la potencia mecánica inicial de
valor 0,8 pu.
Por tanto la igualdad de áreas
según la expresión.
23/04/2016 cchilet@continental.edu.pe 58
Solución 1
Igualando ambas áreas, integrando y sustituyendo se
obtiene que:
0,8 ( δcrítico - 0,328) = 2,487 [cos (δcrítico) -cos(π-0,328)] + 0,8
[δcrítico -(π-0,328)]
Si se resuelve la ecuación se obtiene que:
δcrítico = 1,718 radianes
Se puede observar como para valores del ángulo mayores de
(π−δ0) ya no es posible integrar más área de deceleración (A2),
ya que la potencia mecánica resulta mayor que la potencia
eléctrica activa entregada por el generador.
23/04/2016 cchilet@continental.edu.pe 59
Solución 1
Para calcular el tiempo crítico
hay que considerar que
cuando la línea está
desconectada, la potencia
eléctrica se reduce a cero y la
potencia mecánica permanece
constante (0,8 pu), por lo
tanto la ecuación de la Cuya solución se obtiene
simplemente integrando:
oscilación durante este tiempo
será:
23/04/2016 cchilet@continental.edu.pe 60
Solución 1
Como la derivada del ángulo Por otra parte para t=0, el
del rotor es nula para t=0 ángulo del rotor tomará el
(instante de inicio de la falta), valor δ0, por tanto:
se tiene que C=0. Integrando
de nuevo se obtiene que:
23/04/2016 cchilet@continental.edu.pe 61
Solución 1
Como caso particular de la Despejando se obtiene:
ecuación anterior, si se conoce
el ángulo del rotor crítico, se
puede calcular el tiempo
crítico de reconexión que
garantiza la estabalidad
transitoria del sistema, es
decir: Sustituyendo los valores
conocidos se tiene que:
23/04/2016 cchilet@continental.edu.pe 62
Solución 1
El tiempo crítico de reconexión engloba, contando
desde el instante t=0 en que se produce el
cortocircuito, tanto lo que tardan las protecciones en
abrir la línea en defecto, como el tiempo que
permanece abierta la línea hasta que al cerrar los
interruptores se restablece de nuevo el servicio.
El reenganche de la línea se produce con éxito gracias a
que en las líneas aéreas son muy frecuentes los
defectos temporales (Por ejem:caída de una rama de
árbol sobre la línea que se carboniza por la propia
circulación de corriente).
23/04/2016 cchilet@continental.edu.pe 63
Ejemplo 2
• El generador entrega una potencia activa de 1,0 pu
a un nudo de potencia infinita U0=1pu, con una
tensión en terminales del generador Ua =1pu.
23/04/2016 cchilet@continental.edu.pe 64
Ejemplo 2
23/04/2016 cchilet@continental.edu.pe 65
Ejemplo 2
• Calcular la potencia máxima que puede
suministrar el generador tanto durante la falta
como una vez despejada mediante la apertura
de la línea en defecto.
• Calcular también el ángulo crítico de despeje,
o ángulo máximo para el que los interruptores
situados en los extremos de la línea en
defecto pueden disparar con objeto de que se
mantenga el sincronismo.
23/04/2016 cchilet@continental.edu.pe 66
Factores que afectan a la estabilidad
23/04/2016 cchilet@continental.edu.pe 67
Observación
• Una forma de mejorar la estabilidad transitoria es
incrementando el límite de potencia a transmitir en
condiciones normales de funcionamientos,
(aumentar el límite de estabilidad estática.
• Entre los nudos 1 y 2 de un SEP, unidos mediante
una LT de capacidad despreciable e inductancia, X,
la Pmax que se puede transmitir es:
23/04/2016 cchilet@continental.edu.pe 68
Incremento de Pmax
23/04/2016 cchilet@continental.edu.pe 69
Factores que afectan a la estabilidad
transitoria
a) Actuación rápida de las protecciones y
disminución del tiempo de reenganche de la
línea en defecto.
b) Utilización de interruptores unipolares.
c) Utilización de sistemas de excitación o válvulas
de admisión rápidas.
d) La carga en el momento de producirse la
perturbación.
e) Las características del generador.
23/04/2016 cchilet@continental.edu.pe 70
Tiempo de protecciones
• En caso de un cortocircuito la estabilidad
transitoria mejorará cuanto más rápido
actúen las protecciones que desconectan la
línea en defecto. Hay que recordar que
mientras que la línea está en cortocircuito la
potencia eléctrica transmitida será nula o
muy pequeña , y que el rotor se acelerará.
23/04/2016 cchilet@continental.edu.pe 71
Tiempo de protecciones
23/04/2016 cchilet@continental.edu.pe 72
Tiempo de protecciones
23/04/2016 cchilet@continental.edu.pe 73
Tiempo de protecciones
23/04/2016 cchilet@continental.edu.pe 74
Tiempo de protecciones
23/04/2016 cchilet@continental.edu.pe 75
Tiempo de protecciones
Situaciones frecuentes de este tipo resultan:
• De la electrocución de pájaros que puentean
conductores de fases distintas o una fase con
elementos a tierra,
• De la caída de ramas de árbol sobre las líneas
construidas con conductores desnudos
• Del efecto de las sobretensiones
atmosféricas debidas a las caídas de rayos.
23/04/2016 cchilet@continental.edu.pe 76
Tiempo de protecciones
• En cualquier caso la aceleración del rotor
será tanto menor cuanto más rápidamente se
reenganche la línea afectada (mientras que la
LT está abierta la P transmitida es nula y por
tanto el rotor se acelera).
• En definitiva y suponiendo que el defecto se
ha extinguido, cuanto más corto sea el ciclo
de reenganche, mejores serán las
condiciones de estabilidad transitoria.
23/04/2016 cchilet@continental.edu.pe 77
Tiempo de protecciones
• Por esta razón y sobre todo en MT, las líneas
se explotan frecuentemente mediante un
sistema con dos o más ciclos de reenganche:
– 1ro. rápido (≤300 ms) que permite en muchos
casos reestablecer el servicio en la línea una vez
que el defecto se ha quemado, y
– 2do e incluso un 3ro más lento, (con un t variable
> 1s) que permite realizar un 2do o 3er intento de
cierre de los CB en caso de que el defecto hubiera
persistido tras el 1er ciclo de reenganche.
23/04/2016 cchilet@continental.edu.pe 78
Tiempo de protecciones
23/04/2016 cchilet@continental.edu.pe 79
Uso de interruptores unipolares
23/04/2016 cchilet@continental.edu.pe 80
Uso de sistemas de excitación o
válvulas de admisión rápidas
• Otra forma de mejorar la estabilidad
transitoria consiste en aumentar
rápidamente la Iexc (en <1 s), lo cual es
posible con empleo de modernos sistemas de
control.
• Por lo tanto en caso de una perturbación se
aumenta la tensión interna del generador (E´)
y también la potencia que se puede
transmitir.
23/04/2016 cchilet@continental.edu.pe 81
Uso de sistemas de excitación o
válvulas de admisión rápidas
23/04/2016 cchilet@continental.edu.pe 82
Reducción de la carga
23/04/2016 cchilet@continental.edu.pe 83
Características del generador