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

Sintesi Tecnica: L'effetto delle linee HVDC nelle reti elettriche tramite modellazione di Kuramoto

Problematica
Le moderne reti elettriche, tradizionalmente modellate come grandi macchine sincrone basate su generatori rotanti, stanno diventando sempre più vulnerabili a causa dell'integrazione di fonti di energia rinnovabile e della riduzione dell'inerzia rotazionale. Questa riduzione diminuisce la capacità del sistema di smorzare le perturbazioni nel flusso di potenza e nella frequenza, portando a potenziali instabilità, perdita di sincronizzazione e guasti a cascata (blackout). Sebbene le linee ad alta tensione in corrente continua (HVDC) siano note per trasmettere energia efficientemente su lunghe distanze e possano segmentare le reti per limitare la propagazione dei guasti a cascata, il loro impatto sulla dinamica di sincronizzazione e i potenziali effetti avversi (come il paradosso di Braess) richiedono un'indagine più approfondita che vada oltre le approssimazioni lineari. Il documento affronta la necessità di comprendere come la sostituzione di specifici collegamenti AC con connessioni HVDC influenzi la stabilità di sincronizzazione e il comportamento dei guasti a cascata in un regime non lineare e su larga scala.

Metodologia
Gli autori utilizzano uno studio numerico utilizzando il modello di Kuramoto del secondo ordine adattivo (KM) per simulare la dinamica della rete elettrica europea ad alta tensione (dataset EU2016, $N=13.478$ nodi). Il modello incorpora l'inerzia ($I$), lo smorzamento ($\alpha$) e interazioni di fase non lineari, andando oltre le equazioni di oscillazione linearizzate spesso utilizzate nell'ingegneria elettrica.

Componenti metodologiche chiave:

• Dinamica del modello: Il sistema è governato da equazioni differenziali del secondo ordine dove i nodi rappresentano oscillatori con frequenze proprie estratte da un insieme gaussiano. Il modello include un meccanismo di guasto delle linee basato su una soglia: se il flusso di potenza $|\sin(\theta_j - \theta_i)|$ supera una soglia $T=0,99$, la linea viene rimossa permanentemente, innescando una cascata.
• Implementazione HVDC: Lo studio simula la sostituzione di 36 cavi sottomarini che collegano la penisola scandinava al resto d'Europa (circa lo 0,2% dei collegamenti totali) con linee HVDC. Vengono testati due approcci di modellazione primari:
  1. HVDC Statica: I collegamenti vengono rimossi dal grafo di adiacenza e vengono applicate iniezioni/prelievi di potenza fissi ai terminali, deccoppiando efficacemente le fasi delle due sottoreti.
  2. HVDC Adattiva: Viene aggiunto un nuovo termine alle equazioni del moto, $f(\dot{\theta}_j - \dot{\theta}_i) \cdot D_{ij}$, dove il flusso di potenza dipende dalla differenza di frequenza locale. Diverse funzioni di attivazione ($f$) sono testate, tra cui lineare, $\tanh(x)$, $\text{softsign}(x)$ e $\text{sign}(x)$.
• Protocollo di simulazione: Le simulazioni sono eseguite su cluster GPU e CPU utilizzando integratori a passo adattivo (Bulirsch-Stoer e Runge-Kutta). Il sistema viene inizializzato in uno stato di fase sincronizzata, termalizzato verso uno stato stazionario, e poi perturbato dalla rimozione di una linea casuale per osservare l'evoluzione della cascata. Le metriche includono il parametro d'ordine di Kuramoto ($R$), il parametro d'ordine universale ($R_{uni}$), la varianza di frequenza ($\Omega$) e la dimensione media della cascata ($\langle N_f \rangle$).

Contributi Chiave

1. Analisi Non Perturbativa: Lo studio fornisce una simulazione dinamica del modello di Kuramoto del secondo ordine su una rete realistica su larga scala, catturando gli effetti non lineari e le grandi differenze di fase che i modelli lineari trascurano.
2. Contesto della Dimensione Spettrale: I risultati sono interpretati attraverso la lente della dimensione spettrale ($d_s \approx 2,34$ per la rete EU). Gli autori collegano i miglioramenti della sincronizzazione al comportamento di scaling delle dimensioni finite dei modelli di Kuramoto del secondo ordine su grafi con $d_s < 4$, dove la forza di accoppiamento critica $K_c$ diverge con la dimensione del sistema.
3. Modellazione Comparativa HVDC: Il documento confronta sistematicamente la segmentazione HVDC statica con varie strategie di controllo adattivo (basate sulle differenze di frequenza), offrendo una prospettiva da fisico sui complementi alla stabilità di rete delle analisi di carico dell'ingegneria elettrica esistenti.

Risultati

• Sincronizzazione e Stabilità:
  • HVDC Statica: Sostituire i collegamenti AC con connessioni DC statiche migliora la sincronizzazione di fase (maggiore $R$) e riduce le dimensioni delle cascate per forze di accoppiamento intermedie. Tuttavia, ciò avviene a scapito di una maggiore dispersione di frequenza ($\Omega$) poiché le linee DC non trasmettono informazioni di fase o frequenza, creando efficacemente due sottoreti indipendenti con maggiori varianze di frequenza individuali.
  • HVDC Adattiva: I metodi adattivi migliorano generalmente la sincronizzazione a regime stazionario. La funzione di attivazione $\text{sign}(x)$ (funzione a gradino) fornisce la stabilizzazione più robusta in tutti i valori di accoppiamento, eliminando il paradosso di Braess della frequenza osservato nei modelli statici ad alto accoppiamento. Tuttavia, questo metodo introduce discontinuità che aumentano drasticamente i tempi di rilassamento (di ordini di grandezza) e causano una convergenza lenta.
  • Attivazioni Lineari e Smooth: Le funzioni lineare, $\tanh$ e $\text{softsign}$ migliorano la sincronizzazione per valori di accoppiamento inferiori ($K < 6000$), ma non riescono a mitigare le dimensioni delle cascate ad accoppiamenti più elevati, mostrando un comportamento di crossover. Restringere l'intervallo di attivazione lineare spinge il punto di saturazione più in basso, riducendo le dimensioni delle cascate.
• Comportamento delle Cascate: Sebbene l'HVDC riduca generalmente le dimensioni delle cascate rispetto al modello AC completo, l'efficacia dipende dalla forza di accoppiamento e dalla specifica strategia di controllo HVDC. I metodi adattivi mostrano potenziale ma soffrono di lunghi tempi di rilassamento.
• Effetti Braess: Lo studio nota che variare la potenza totale trasmessa ad accoppiamenti globali grandi e piccoli può portare a svantaggi nella dispersione di frequenza e negli effetti Braess, particolarmente quando le fluttuazioni sono piccole, causando un "congelamento" nella dinamica.

Significatività e Rivendicazioni
Gli autori affermano che il loro lavoro fornisce un complemento più profondo, basato sulla fisica, alle analisi esistenti dell'ingegneria elettrica delle reti elettriche. Risolvendo il modello di Kuramoto del secondo ordine adattivo senza approssimazioni lineari, dimostrano che:

• La dimensione spettrale della rete ($d_s < 4$) influenza fondamentalmente le transizioni di sincronizzazione, suggerendo che segmentare la rete (tramite HVDC) in sottografi più piccoli può migliorare la stabilità di sincronizzazione per un dato accoppiamento.
• Le linee HVDC adattive, in particolare quelle basate sul controllo della differenza di frequenza, sono più efficienti nello stato stazionario rispetto alle sostituzioni statiche, ma comportano la penalità di tempi di rilassamento molto lunghi.
• Il compromesso tra sincronizzazione di fase e dispersione di frequenza è un fattore critico; mentre l'HVDC migliora l'allineamento di fase, può degradare la coerenza di frequenza.

Il documento conclude che, sebbene la segmentazione tramite HVDC sia una strategia percorribile per mitigare i guasti a cascata e migliorare la sincronizzazione in regimi specifici, la scelta della strategia di controllo (statica vs adattiva) comporta significativi compromessi tra stabilità, tempo di rilassamento e coerenza di frequenza. Gli autori suggeriscono con modestia che lavori futuri dovrebbero investigare gli effetti dei rumori stocastici additivi e confrontare le distribuzioni di frequenza con dati reali.