The effect of HVDC lines in power-grids via Kuramoto modelling

Resumen Técnico: El Efecto de las Líneas HVDC en las Redes Eléctricas mediante el Modelado de Kuramoto

Planteamiento del Problema
Las redes eléctricas modernas, tradicionalmente modeladas como grandes máquinas síncronas basadas en generadores rotativos, se están volviendo cada vez más vulnerables debido a la integración de fuentes de energía renovables y la reducción de la inercia rotacional. Esta reducción disminuye la capacidad del sistema para amortiguar las perturbaciones en el flujo de potencia y la frecuencia, lo que conduce a posibles inestabilidades, pérdida de sincronización y fallos en cascada (apagones). Si bien se sabe que las líneas de Corriente Continua de Alta Tensión (HVDC) transmiten la energía de manera eficiente a través de largas distancias y pueden segmentar las redes para limitar la propagación de cascadas, su impacto en la dinámica de sincronización y los posibles efectos adversos (como la paradoja de Braess) requieren una investigación más profunda que vaya más allá de las aproximaciones lineales. El artículo aborda la necesidad de comprender cómo el reemplazo de enlaces de CA específicos por conexiones HVDC afecta la estabilidad de la sincronización y el comportamiento de los fallos en cascada en un régimen no lineal y de gran escala.

Metodología
Los autores emplean un estudio numérico utilizando el modelo de Kuramoto de segundo orden adaptativo (KM) para simular la dinámica de la red de alta tensión (HV) europea (conjunto de datos EU2016, $N=13,478$ nodos). El modelo incorpora inercia ($I$), amortiguamiento ($\alpha$) e interacciones de fase no lineales, yendo más allá de las ecuaciones de oscilación linealizadas que se utilizan frecuentemente en ingeniería eléctrica.

Componentes metodológicos clave:

• Dinámica del Modelo: El sistema está gobernado por ecuaciones diferenciales de segundo orden donde los nodos representan osciladores con frecuencias propias extraídas de un conjunto gaussiano. El modelo incluye un mecanismo de fallo de línea basado en un umbral: si el flujo de potencia $|\sin(\theta_j - \theta_i)|$ supera un umbral $T=0.99$, la línea se elimina permanentemente, desencadenando una cascada.
• Implementación HVDC: El estudio simula el reemplazo de 36 cables submarinos que conectan la península escandinava con el resto de Europa (aprox. 0.2% del total de enlaces) por líneas HVDC. Se prueban dos enfoques principales de modelado:
  1. HVDC Estático: Los enlaces se eliminan del grafo de adyacencia y se aplican inyecciones/extracciones de potencia fijas en los extremos, desacoplando efectivamente las fases de las dos subredes.
  2. HVDC Adaptativo: Se añade un nuevo término a las ecuaciones de movimiento, $f(\dot{\theta}_j - \dot{\theta}_i) \cdot D_{ij}$, donde el flujo de potencia depende de la diferencia de frecuencia local. Se prueban varias funciones de activación ($f$), incluyendo lineales, $\tanh(x)$, $\text{softsign}(x)$ y $\text{sign}(x)$.
• Protocolo de Simulación: Las simulaciones se ejecutan en clústeres de GPU y CPU utilizando integradores de paso adaptativo (Bulirsch-Stoer y Runge-Kutta). El sistema se inicializa en un estado de sincronización de fase, se termaliza hacia un estado estacionario y luego se perturba mediante la eliminación de una línea aleatoria para observar la evolución de la cascada. Las métricas incluyen el parámetro de orden de Kuramoto ($R$), el parámetro de orden universal ($R_{uni}$), la varianza de la frecuencia ($\Omega$) y el tamaño promedio de la cascada ($\langle N_f \rangle$).

Contribuciones Clave

1. Análisis No Perturbativo: El estudio proporciona una simulación dinámica del modelo de Kuramoto de segundo orden en una red realista y de gran escala, capturando efectos no lineales y grandes diferencias de fase que los modelos lineales omiten.
2. Contexto de la Dimensión Espectral: Los resultados se interpretan a través de la lente de la dimensión espectral ($d_s \approx 2.34$ para la red de la UE). Los autores relacionan las mejoras de sincronización con el comportamiento de escalamiento de tamaño finito de los modelos de Kuramoto de segundo orden en grafos con $d_s < 4$, donde la fuerza de acoplamiento crítica $K_c$ diverge con el tamaño del sistema.
3. Modelado Comparativo de HVDC: El artículo compara sistemáticamente la segmentación estática de HVDC frente a diversas estrategias de control adaptativo (basadas en diferencias de frecuencia), ofreciendo una perspectiva de la física sobre los complementos de estabilidad de la red respecto a los estudios de flujo de carga de la ingeniería eléctrica existentes.

Resultados

• Sincronización y Estabilidad:
  • HVDC Estático: Reemplazar los enlaces de CA por conexiones de CC estáticas mejora la sincronización de fase (mayor $R$) y reduce los tamaños de cascada para fuerzas de acoplamiento intermedias. Sin embargo, esto conlleva un aumento en la dispersión de la frecuencia ($\Omega$) porque las líneas de CC no transmiten información de fase o frecuencia, creando efectivamente dos subredes independientes con mayores varianzas de frecuencia individuales.
  • HVDC Adaptativo: Los métodos adaptativos generalmente mejoran la sincronización en estado estacionario. La función de activación $\text{sign}(x)$ (función escalón) proporciona la estabilización más robusta en todos los valores de acoplamiento, eliminando la paradoja de Braess de frecuencia observada en los modelos estáticos a altos niveles de acoplamiento. No obstante, este método introduce discontinuidades que aumentan drásticamente los tiempos de relajación (por órdenes de magnitud) y causan una convergencia lenta.
  • Activaciones Lineales y Suaves: Las funciones lineales, $\tanh$ y $\text{softsign}$ mejoran la sincronización para valores de acoplamiento bajos ($K < 6000$), pero no logran mitigar los tamaños de cascada en acoplamientos más altos, mostrando un comportamiento de cruce. Estrechar el intervalo de activación lineal desplaza el punto de saturación hacia valores más bajos, reduciendo los tamaños de cascada.
• Comportamiento de Cascada: Aunque el HVDC generalmente reduce los tamaños de cascada en comparación con el modelo completo de CA, la efectividad depende de la fuerza de acoplamiento y de la estrategia específica de control de HVDC. Los métodos adaptativos resultan prometedores pero sufren de largos tiempos de relajación.
• Efectos de Braess: El estudio señala que variar la potencia total transmitida en acoplamientos globales grandes y pequeños puede derivar en desventajas en la dispersión de la frecuencia y efectos de Braess, particularmente cuando las fluctuaciones son pequeñas, causando un "congelamiento" en la dinámica.

Significancia y Afirmaciones
Los autores afirman que su trabajo proporciona un complemento más profundo, basado en la física, a los análisis de ingeniería eléctrica existentes para las redes eléctricas. Al resolver el modelo de Kuramoto de segundo orden adaptativo sin aproximaciones lineales, demuestran que:

• La dimensión espectral de la red ($d_s < 4$) influye fundamentalmente en las transiciones de sincronización, sugiriendo que segmentar la red (vía HVDC) en subgrafos más pequeños puede mejorar la estabilidad de la sincronización para una fuerza de acoplamiento dada.
• Las líneas HVDC adaptativas, particularmente aquellas que utilizan el control basado en la diferencia de frecuencia, son más eficientes en el estado estacionario que los reemplazos estáticos, pero conllevan la penalización de tiempos de relajación muy largos.
• El compromiso entre la sincronización de fase y la dispersión de la frecuencia es un factor crítico; si bien el HVDC mejora el alineamiento de fase, puede degradar la coherencia de la frecuencia.

El artículo concluye que, aunque la segmentación mediante HVDC es una estrategia viable para mitigar los fallos en cascada y mejorar la sincronización en regímenes específicos, la elección de la estrategia de control (estática vs. adaptativa) implica compromisos significativos entre estabilidad, tiempo de relajación y coherencia de frecuencia. Los autores sugieren modestamente que trabajos futuros deberían investigar los efectos de ruidos estocásticos aditivos y comparar las distribuciones de frecuencia con datos del mundo real.