Instabilidade por Ressonância Subsíncrona (SSR) em geradores com compensação em série

Fundamentos Electromecánicos y Origen de la Resonancia Subsíncrona

La inestabilidad por resonancia subsíncrona (SSR, por sus siglas en inglés) representa uno de los fenómenos más complejos y destructivos en los sistemas eléctricos de potencia de corriente alterna de muy alta tensión (HVAC). Ocurre típicamente en plantas de generación síncrona conectadas a redes de transmisión de ultra larga distancia que emplean compensación serie capacitiva. El propósito fundamental de la compensación serie es reducir la reactancia inductiva total de la línea de transmisión, mejorando así la capacidad de transferencia de potencia activa, optimizando el ángulo de fase y mejorando la estabilidad transitoria del sistema interconectado. Sin embargo, la introducción de bancos de capacitores en serie con las inductancias de la red crea circuitos LC resonantes serie cuyas frecuencias naturales se ubican por debajo de la frecuencia fundamental del sistema industrial (50 Hz o 60 Hz).

Desde la perspectiva de la teoría de circuitos distribuidos y concentrados, la frecuencia natural de oscilación eléctrica de un sistema de transmisión compensado en serie se define mediante la siguiente expresión fundamental:

fer = f_0 \sqrt{ \frac{X_C}{X_L} } = f_0 \sqrt{ kse }

Donde $f_0$ representa la frecuencia fundamental del sistema (50 Hz o 60 Hz), $X_C$ es la reactancia capacitiva total de los bancos serie, $X_L$ es la reactancia inductiva total equivalente del sistema vista desde los terminales del generador, y $kse$ es el grado de compensación serie expresado en por unidad o porcentaje. Dado que los grados de compensación en líneas largas varían comúnmente entre el 25% y el 70%, la frecuencia eléctrica resonante $fer$ resulta ser estrictamente menor que la frecuencia fundamental. En consecuencia, cualquier perturbación transitoria o armónica en la red excita esta frecuencia eléctrica natural, dando lugar a corrientes y tensiones subsíncronas.

Cuando estas corrientes subsíncronas circulan por los devanados estatóricos del generador síncrono, generan un campo magnético giratorio en el entrehierro de la máquina que gira a una velocidad angular menor que la velocidad síncrona del rotor. Este campo giratorio subsíncrono interactúa con el campo magnético de corriente continua del rotor, produciendo pares electromagnéticos oscilatorios en el eje del generador cuya frecuencia de pulsación $f_m$ viene dada por la diferencia entre la frecuencia fundamental y la frecuencia eléctrica resonante:

f_m = f_0 - fer = f_0 \left( 1 - \sqrt{kse} \right)

Si la frecuencia de este par oscilatorio coincide o se aproxima notablemente a una de las frecuencias naturales de vibración torsional del tren de ejes de la unidad generadora (el cual comprende la turbina de alta presión, media presión, baja presión y el rotor del generador y la excitatriz), se produce un acoplamiento mecánico-eléctrico severo. Este mecanismo se conoce formalmente como Interacción Torsional de Modos (Torsional Interaction - TI). La máquina actúa como un amplificador de energía mecánica a través de la realimentación negativa o positiva del amortiguamiento torsional, pudiendo llevar a la fatiga torsional extrema y a la fractura catastrófica del eje en cuestión de segundos.

Clasificación y Mecanismos Físicos de la SSR

El fenómeno global de la resonancia subsíncrona no es un evento monolítico, sino que se subdivide formalmente en tres mecanismos físicos distintos pero interrelacionados según el IEEE SSR Task Force:

Efecto de Inducción de Rotor (Rotor Induction Generator Effect - RIGE)

Este fenómeno es de naturaleza estrictamente eléctrica y puede ocurrir incluso con un rotor idealmente rígido (sin grados de libertad torsional). Cuando corrientes subsíncronas de frecuencia $fer$ fluyen por el estator, establecen un campo magnético que gira a una velocidad menor que la del rotor. Visto desde el sistema de referencia del rotor girando a la velocidad síncrona \omega_0, este campo aparece con una frecuencia de deslizamiento subsíncrona $fslip = f_0 - fer$ . Esta frecuencia relativa induce corrientes en los circuitos amortiguadores (devanados damper) y en el cuerpo masivo del rotor de acero macizo.

La impedancia equivalente vista desde el estator a la frecuencia subsíncrona incluye la resistencia efectiva del rotor reflejada. Debido al deslizamiento negativo para frecuencias menores que la fundamental ( $fer < f_0$ ), la resistencia aparente del rotor vista por el estator a esa frecuencia subsíncrona puede volverse negativa. Una resistencia total de estator y rotor negativa a la frecuencia $fer$ significa que el circuito absorbe energía de la red y la amplifica, o que el circuito amortigua negativamente las oscilaciones eléctricas, provocando autoexcitación eléctrica si la resistencia neta del circuito Rtotal(fer) < 0.

Interacción Torsional de Modos (Torsional Interaction - TI)

La Interacción Torsional ocurre cuando el sistema eléctrico y el tren de ejes mecánico se acoplan a través del campo magnético del entrehierro. Supóngase que el tren de ejes posee una frecuencia natural torsional $fmech$ . Si ocurre una oscilación torsional en el eje a dicha frecuencia, las variaciones angulares del rotor modulan las corrientes de estator, generando bandas laterales en el espectro eléctrico a las frecuencias f_0 \pm fmech. Si una de estas frecuencias laterales coincide con la frecuencia natural complementaria del circuito eléctrico serie ( $fer = f_0 - fmech$ ), se establece un ciclo de realimentación cerrado.

Matemáticamente, este acoplamiento se evalúa mediante el cálculo del coeficiente de amortiguamiento eléctrico torsional $D_e$ para cada modo torsional $i$ . El decremento logarítmico o amortiguamiento total $Dtotal$ de un modo torsional es la suma del amortiguamiento mecánico inherente $D_m$ y el amortiguamiento eléctrico inducido $D_e$ :

D_{total, i} = D_{m, i} + D_{e, i}

Si la resistencia de la red y las características de los circuitos LC producen que $D_{e, i}$ sea negativo y su magnitud supere al amortiguamiento mecánico positivo ( $|D_{e, i}| > D_{m, i}$ ), el amortiguamiento total del sistema se vuelve negativo ( $D_{total, i} < 0$ ). Esto conduce a una inestabilidad exponencial donde las oscilaciones torsionales crecen de manera descontrolada hasta la destrucción mecánica de los acoplamientos del eje.

Amplificación por Transitorios de Torque (Torque Amplification - TA)

A diferencia del RIGE y la TI que son fenómenos de inestabilidad dinámica y de pequeña señal, la Amplificación por Transitorios de Torque es un fenómeno transitorio de gran magnitud. Ocurre cuando se producen cortocircuitos severos en las proximidades de líneas compensadas con capacitores serie, seguidos por operaciones de apertura y recierre de interruptores de potencia. La brusca descarga y recarga de los bancos de capacitores serie interactúa con los armónicos transitorios y las corrientes de falla, inyectando un impulso de torque electromagnético severo en el eje del generador.

Si la frecuencia de este impulso coincide con los modos naturales del eje, el tren de ejes experimenta esfuerzos mecánicos que pueden superar ampliamente el par nominal de diseño, acumulando daño por fatiga plástica severa en un solo evento transitorio.

Modelado Matemático y Análisis de Pequeña Señal en Sistemas Multimodo

Para realizar un análisis riguroso de estabilidad por resonancia subsíncrona en entornos industriales y de transmisión avanzados, se emplea el modelado de espacio de estados linealizado alrededor de un punto de operación en régimen permanente. El sistema completo se descompone en tres subsistemas acoplados: la red eléctrica de transmisión con compensación serie, el generador síncrono con sus ecuaciones de Park (ejes d-q), y el modelo mecánico multimasas del tren de ejes.

El modelo mecánico del tren de ejes con $N$ masas rotatorias se representa mediante un sistema de ecuaciones diferenciales ordinarias de segundo orden:

J_i \frac{d^2 \theta_i}{dt^2} + Dmi \frac{d \theta_i}{dt} + \sum_{j=1}^{N} Kij (\theta_i - \theta_j) = Tmi - Tei, \quad i = 1, 2, \dots, N

Donde $J_i$ es el momento de inercia de la masa $i$ , \theta_i es el ángulo angular del rotor, $Dmi$ es el coeficiente de amortiguamiento propio, $Kij$ es la rigidez torsional del eje entre las masas $i$ y $j$ , $Tmi$ es el par mecánico aplicado por la turbina, y $Tei$ es el par electromagnético inducido sobre la masa del generador.

El sistema eléctrico comprende las ecuaciones de transformaciones de Park para el estator y el rotor, acopladas con las ecuaciones diferenciales de los circuitos de transmisión R-L-C serie:

\frac{d i_d}{dt} = \frac{1}{L_d} \left( -R_a i_d + \omega \psi_q - v_d \right)

\frac{d i_q}{dt} = \frac{1}{L_q} \left( -R_a i_q - \omega \psi_d - v_q \right)

Las ecuaciones de la línea de transmisión compensada se expresan en coordenadas de ejes sincronizados con la red mediante ecuaciones diferenciales que incorporan la dinámica de carga y los capacitores serie:

\frac{d vcd}{dt} = \omega_0 i_d + \frac{1}{Cse} i_d

\frac{d vcq}{dt} = -\omega_0 i_q + \frac{1}{Cse} i_q

Al linealizar este conjunto de ecuaciones no lineales alrededor de un punto de operación estable ( $X_0, U_0$ ), se obtiene la matriz de estado del sistema cerrado:

\Delta \dot{X} = A \Delta X + B \Delta U

\Delta Y = C \Delta X + D \Delta U

El análisis de estabilidad de pequeña señal se efectúa calculando los valores propios (eigenvalues) \lambda_i = \sigma_i \pm j \omega_i de la matriz de estado $A$ . La parte real \sigma_i de cada valor propio determina la estabilidad del modo correspondiente:

Si \sigma_i < 0, el modo es asintóticamente estable (amortiguamiento positivo).
Si \sigma_i > 0, el modo es inestable (amortiguamiento negativo, crecimiento exponencial de la oscilación).
Si \sigma_i = 0, el modo se encuentra en el límite de estabilidad marginal.

Para los modos electromecánicos torsionales, la parte real del valor propio está directamente relacionada con el coeficiente de amortiguamiento total: \sigma_i = -\frac{D_{total, i}}{2 J_i}. Por lo tanto, un valor propio con parte real positiva indica directamente la ocurrencia de inestabilidad por Interacción Torsional (TI).

Matriz Comparativa de Normativa Internacional y Criterios de Evaluación

La evaluación y mitigación de la SSR están regidas por estrictos estándares internacionales emitidos por el IEEE y la IEC. La siguiente matriz comparativa detalla los parámetros críticos, límites normativos y consecuencias operativas asociados a la SSR en redes de alta tensión.

Norma / Estándar	Parámetro de Evaluación	Límite Crítico / Criterio de Aceptación	Condición de Falla Asociada	Consecuencia Operativa y Dieléctrica
IEEE Std 1133	Amortiguamiento Eléctrico Torsional ( $D_e$ )	Dtotal \ge 0 bajo cualquier condición de despacho ( $N-0, N-1$ )	Inestabilidad por Interacción Torsional (TI)	Fatiga torsional acumulada, deformación y fractura catastrófica de ejes.
IEEE Std 693	Respuesta Sísmica y Torsional Acoplada	Esfuerzos mecánicos máximos menores al 80% del límite elástico	Amplificación por transitorios de torque combinados	Daño estructural en soportes de turbogeneradores y cimentaciones de subestaciones.
IEC 60383 / IEC 60871	Sobretensiones en Bancos de Capacitores Serie	Vpeak \le 1.8 p.u. en régimen permanente subsíncrono	Resonancia subsíncrona prolongada y ferroresonancia	Perforación dieléctrica de los elementos capacitores, disparo de los VMS (Varistor Metal-Oxide).
IEEE Std 1547 / C50.13	Corrientes Subsíncronas en Devanados de Estator	Isub \le 0.05 p.u. continuo; Isub \le 0.15 p.u. por 10s	Efecto de inducción de rotor (RIGE) y calentamiento focalizado	Sobrecalentamiento crítico en cuñas del rotor, anillos de retención y núcleo estatórico.

Análisis Forense de Fallas en Equipos Críticos ante la SSR

La presencia de corrientes y pares subsíncronos genera solicitaciones severas que trascienden el generador síncrono, afectando profundamente a toda la infraestructura de la subestación asociada y a los sistemas de transformación y maniobra. El análisis forense de fallas revela los mecanismos de degradación en los siguientes componentes:

Transformadores de Potencia y Autotransformadores

Las componentes de corriente subsíncrona ( $fer$ ) superpuestas a la frecuencia fundamental generan asimetrías temporales en el flujo magnético del núcleo del transformador. Dado que el flujo magnético es la integral temporal de la tensión, la presencia de componentes de baja frecuencia magnotectiva eleva la densidad de flujo pico \hat{B}, empujando al núcleo profundamente hacia la región de saturación magnética no lineal.

La saturación periódica del núcleo provoca:

Generación masiva de armónicos de orden superior (especialmente 2do, 3ro y 5to armónico) que distorsionan la onda de tensión.
Incremento drástico de las corrientes de excitación y de las pérdidas en el fierro por histéresis y corrientes parásitas (Foucault).
Vibraciones mecánicas severas en los tanques y devanados del transformador debido a la magnetostricción no lineal y a las fuerzas de Lorentz ponderadas por las corrientes asimétricas, lo que resulta en el aflojamiento de los aprietes mecánicos y la falla prematura del aislamiento sólido (papel impregnado en aceite).

Interruptores de Potencia y Sistemas de Maniobra

Los interruptores de potencia ubicados en las líneas compensadas con capacitores serie experimentan condiciones de restablecimiento de tensión extremadamente severas durante la ocurrencia de SSR. Las tensiones transitorias de recuperación (TRV) se ven moduladas por la frecuencia natural del circuito serie $fer$ , presentando tasas de crecimiento de tensión transitoria ( $du/dt$ ) muy elevadas y valores pico de tensión que pueden superar la capacidad nominal de interrupción del interruptor.

Adicionalmente, cuando un interruptor opera para despejar una falla en presencia de corrientes subsíncronas, la probabilidad de reencendido (re-strike) o reconducción del arco eléctrico aumenta considerablemente debido a que la corriente pasa por cero a velocidades anormales o contiene componentes bidireccionales superpuestas de baja frecuencia. Esto puede causar la destrucción catastrófica de las cámaras de extinción de arco de SF6 y la falla explosiva del polo del interruptor.

Cables de Alta Tensión y Sistemas de Medición

Los cables aislados de alta tensión interconectados a la subestación sufren un incremento en las pérdidas dieléctricas y en las pérdidas por corrientes inducidas en las pantallas metálicas (apantallamiento de cobre o plomo). Las corrientes subsíncronas inducen campos magnéticos alternos de baja frecuencia que atraviesan las pantallas, generando corrientes circulatorias elevadas que sobrecalientan el aislamiento polimérico (XLPE). Esto acelera la formación de arborizaciones eléctricas (water/electrical treeing) y reduce drásticamente la vida útil estimada del cable.

Por otro lado, los transformadores de corriente (CT) y de potencial (PT) sufren errores de relación y ángulo de fase significativos debido a la saturación magnética inducida por las componentes de corriente subsíncrona, lo que inhabilita o degrada severamente el desempeño de los relés de protección numéricos basados en fasores de frecuencia fundamental.

Estrategias Avanzadas de Mitigación y Diseño de Ingeniería

Para garantizar la viabilidad operativa de los sistemas de transmisión con compensación serie masiva, la ingeniería moderna implementa un conjunto de contramedidas tecnológicas avanzadas, divididas en mitigación a nivel de red, a nivel de equipos de generación y mediante sistemas electrónicos de potencia de alta velocidad.

Filtros de Bloqueo Subsíncrono (Subsynchronous Resonance Blocking Filters - SSBF)

Los filtros SSBF son circuitos sintonizados instalados en serie con los devanados del estator del generador o directamente en la línea de transmisión. Consisten en un arreglo de reactores y capacitores en paralelo sintonizados exactamente para presentar una impedancia extremadamente alta (un circuito abierto virtual) a la frecuencia eléctrica resonante crítica $fer$ .

La impedancia del filtro SSBF a la frecuencia de sintonía se expresa como:

ZSSBF(j\omegaer) = \frac{R + j\omegaer L}{1 - \omegaer^2 L C + j\omegaer R C} \approx \infty