Resumen Técnico: Un Marco de Fasores Dinámicos para el Análisis de Oscilaciones Subsíncronas en Sistemas Multimaquina con IBRs y Grandes Cargas

1. Planteamiento del Problema

La red eléctrica moderna enfrenta desafíos sin precedentes relacionados con la escalabilidad en el análisis de dominio del tiempo y de la frecuencia. Específicamente, la creciente integración de Recursos Basados en Inversores (IBRs) de Seguimiento de Red (Grid-Following, GFL) ha provocado un aumento de las Oscilaciones Subsíncronas (SSO). Simultáneamente, la rápida expansión de las cargas de Centros de Datos (DC) impulsados por Inteligencia Artificial (IA) introduce nuevas preocupaciones de fiabilidad, incluyendo rampas de potencia pronunciadas y ciclos de entrenamiento sincronizados que pueden inducir fluctuaciones de alta amplitud. Estas fluctuaciones amenazan la vida útil por fatiga de los ejes de las turbinas-generadores y afectan la respuesta de frecuencia primaria.

Si bien las herramientas de simulación de Transitorios Electromagnéticos (EMT) pueden modelar estos fenómenos, sufren de una baja escalabilidad debido a las altas demandas computacionales. Además, los modelos EMT son generalmente inadecuados para el desarrollo de los modelos de Tiempo Invariante Lineal (LTI) requeridos para el análisis de valores propios (eigenvalues), la identificación de la causa raíz y el diseño sistemático de controladores. Los marcos existentes de Fasores Dinámicos (DP) se han limitado mayoritariamente a sistemas de una sola máquina o un solo IBR, o no han abordado adecuadamente sistemas multimaquina con IBRs bajo condiciones desequilibradas.

2. Metodología

Los autores proponen un marco de Fasores Dinámicos (DP) generalizado diseñado para analizar las SSO en sistemas multimaquina que contienen IBRs de Seguimiento de Red (GFL) y de Formación de Red (GFM), generadores sincrónicos (SGs) de múltiples masas, y grandes cargas (específicamente centros de datos de IA) tanto en condiciones equilibradas como desequilibradas.

Estrategia de Modelado Central

El marco utiliza un sistema de coordenadas híbrido para aprovechar las fortalezas de diferentes marcos de referencia:

• IBRs (GFL y GFM): Modelados en el marco dq utilizando fasores dinámicos. Esto se alinea con las implementaciones estándar de control vectorial. Los modelos retienen coeficientes DP para $k=0$ (secuencia positiva) y $k=\pm2$ (secuencia negativa) para capturar la dinámica desequilibrada.
  • IBRs GFM: Descompuestos en mecanismos de droop, lazos de voltaje externos, lazos de voltaje/corriente internos y la dinámica del filtro de salida. Las estrategias de limitación de corriente (ángulo constante e $i_q$-prioridad) se modelan explícitamente.
  • IBRs GFL: Modelados con la dinámica del Lazo de Seguimiento de Fase (PLL), lazos de control de potencia activa/reactiva externos, lazos de control de corriente internos y la planta (filtro R-L en serie).
• Red y SGs: Modelados en el marco pnz (secuencias positiva, negativa y cero).
  • Líneas de Transmisión: Representadas mediante modelos de sección $\pi$ concentrada.
  • Generadores Sincrónicos: Se representan modelos detallados de turbina-generador de múltiples masas (incluyendo turbinas HP, IP, LPA, LPB, generador y excitatriz); la dinámica mecánica se captura usando el componente $k=0$, mientras que los transitorios eléctricos utilizan componentes $k=0, \pm2$.
  • Cargas: Las cargas de DC de IA se modelan como demandas de potencia constante con factor de potencia unitario (solo secuencia positiva).
• Acoplamiento: Las variables de interfaz (voltajes y corrientes terminales) se transforman entre el marco dq (IBRs) y el marco pnz (red/SGs) utilizando ecuaciones de transformación específicas derivadas en el artículo.

Propiedades Clave

• Invarianza Temporal: La formulación resulta en un sistema invariante en el tiempo, lo que permite la descomposición de valores propios.
• Linealizabilidad: El modelo puede linealizarse para identificar modos de SSO, realizar análisis de causa raíz y diseñar controladores de amortiguamiento.
• Escalabilidad: Al truncar la serie de Fourier a los coeficientes dominantes, el marco permite simulaciones en el dominio del tiempo más rápidas en comparación con EMT.

3. Contribuciones Clave

El artículo reclama las siguientes contribuciones específicas:

1. Marco Generalizado: Presentación de un marco de modelado basado en DP capaz de analizar fenómenos subsíncronos en sistemas con turbinas de múltiples masas, IBRs GFL y GFM, y la dinámica de líneas de transmisión bajo condiciones equilibradas y desequilibradas.
2. Benchmarking: Demostración de la aceleración del marco respecto a los modelos EMT en un sistema multimaquina.
3. Aplicación de Diseño de Control: Aprovechamiento de la invariancia temporal y la linealizabilidad del marco para diseñar un controlador descentralizado de bajo orden para amortiguar modos de SSO pobremente amortiguados en un sistema IEEE 68-bus modificado.
4. Análisis GFM vs. GFL: Demostración, mediante análisis de dominio del tiempo y de la frecuencia, de que reemplazar un IBR GFL por un IBR GFM elimina modos de SSO específicos pobremente amortiguados.
5. Estudio de Impacto de Carga de IA: Análisis del impacto de las cargas de DC de IA en la respuesta de frecuencia primaria y en las oscilaciones torsionales (esfuerzo cortante) de los ejes de turbina-generador de múltiples masas, proporcionando un método para evaluar los riesgos de vida útil por fatiga.

4. Resultados y Validación

El marco fue validado y probado utilizando tres sistemas principales:

A. Validación contra Modelos EMT

• Modelo de Referencia Primero (FBM) de la IEEE para SSR: Los valores propios y las respuestas en el dominio del tiempo del modelo DP (tras una falla trifásica a tierra) coincidieron estrechamente con los resultados de EMT y los valores de la literatura.
• GFM IBR con Límites de Corriente: Validado contra un modelo EMTDC/PSCAD. El modelo DP replicó con precisión las respuestas a fallas de línea a línea (L-L) de auto-extinción bajo estrategias de limitación de corriente de ángulo constante e $i_q$-prioridad.
• Sistema IEEE de 4 Máquinas Modificado: Se simuló un sistema con dos IBRs GFL reemplazando generadores sincrónicos. El modelo DP (usando linealización) y el modelo EMT (usando análisis de Prony) mostraron una concordancia cercana en la frecuencia de los modos de SSO y en los coeficientes de amortiguamiento para diferentes anchos de banda de PLL.
• Aceleración: En el sistema IEEE de 4 máquinas, el modelo DP logró una aceleración promedio de 2x sobre el modelo EMT. En el contexto del sistema IEEE de 39 buses (referenciado de un trabajo previo), se observó una aceleración de 10x, aunque el presente artículo se centra en el sistema de 68 buses.

B. Amortiguamiento de SSO en el Sistema IEEE 68-Bus Modificado

Se propusieron dos soluciones para amortiguar un modo de SSO pobremente amortiguado (4.76 Hz, 0.2% de amortiguamiento) inducido por IBRs GFL:

1. Controlador Descentralizado: Un controlador de 5º orden diseñado mediante Optimización por Enjambre de Partículas (PSO). Utiliza la magnitud del voltaje terminal como retroalimentación para modular las referencias de corriente reactiva. Las simulaciones mostraron un amortiguamiento efectivo del modo de SSO tras diversas salidas de líneas de interconexión y fallas desequilibradas.
2. Reemplazo con GFM: Reemplazar un IBR GFL por un IBR GFM de la misma capacidad resultó en la desaparición completa del modo de SSO pobremente amortiguado, según confirmó el análisis de valores propios.

C. Impacto de las Cargas de DC de IA

• Respuesta de Frecuencia: Las simulaciones de rampas de carga de DC de IA (transicionando entre periodos de cómputo y de reposo) mostraron desviaciones de frecuencia significativas, activando el soporte de frecuencia ancilar pero no el deslastre de carga.
• Esfuerzo en el Eje: Al sintonizar las oscilaciones de la carga de DC de IA para que coincidan con los modos torsionales (ej. ~7.83 Hz), el estudio demostró que los esfuerzos en el eje podrían alcanzar los límites definidos por el diagrama de Goodman modificado aumentado, lo que podría conducir a la pérdida de vida útil por fatiga.
• Análisis de Sensibilidad: Simulaciones de Monte Carlo (100 ejecuciones) variando la amplitud, frecuencia y fase de las oscilaciones de la carga de IA revelaron que los esfuerzos en el eje son altamente sensibles a pequeñas desviaciones de frecuencia respecto al modo torsional.

D. Rendimiento Computacional

En el sistema IEEE de 68 buses modificado (que contiene 8 SGs de múltiples masas, 8 SGs de una sola masa, 1 GFL, 1 GFM y líneas de transmisión dinámicas):

• El marco DP simuló un escenario de oscilación forzada de 20 segundos proveniente de 7 DCs de IA en un promedio de 2.5 minutos.
• Esto establece la escalabilidad del marco para sistemas a gran escala donde las simulaciones EMT serían computacionalmente prohibitivas.

5. Significación y Reclamaciones

Los autores declaran explícitamente que no se reclama novedad respecto al diseño del control de amortiguamiento de SSO en sí mismo, ni al hecho general de que el control GFM puede remediar los problemas de SSO.

La principal significación de este trabajo radica en demostrar la capacidad del marco basado en DP para:

• Capturar estos fenómenos complejos (SSO multimaquina, fallas desequilibradas e impactos de carga de IA) por primera vez en un entorno unificado, escalable y linealizable.
• Permitir la transición de simulaciones EMT puramente en el dominio del tiempo hacia el análisis en el dominio de la frecuencia (descomposición de valores propios) para sistemas de gran escala con IBRs.
• Proporcionar una herramienta práctica para analizar la vida útil por fatiga de los ejes de turbina-generador bajo las condiciones de estrés específicas inducidas por las cargas modernas de centros de datos de IA.

El artículo concluye que el marco propuesto logra cerrar la brecha entre la precisión de EMT y la utilidad analítica de los modelos lineales, ofreciendo una solución escalable para el análisis y control de las futuras redes eléctricas.