Ingeniería Eléctrica

Ferrorisonanza e saturazione magnetica nei reattori di shunt durante la commutazione monofase

La apertura monofásica de líneas de transmisión compensadas con reactores de derivación es uno de los escenarios operativos más desafiantes en subestaciones de

Ing. Francisco Ramírez

Introduzione e Fondamenti Elettromagnetici dei Reattori di Derivazione

I reattori di derivazione connessi nei sistemi di trasmissione in alta tensione (HV), extra alta tensione (EHV) e altissima tensione (UHV) svolgono un ruolo critico nella compensazione della potenza reattiva capacitiva generata dalle linee di trasmissione lunghe in condizioni di basso carico o a vuoto. Elettromeccanicamente, questi apparecchi sono costituiti da nuclei ferromagnetici con traferri distribuiti, progettati per mantenere una caratteristica di magnetizzazione lineare in condizioni nominali e di sovratensione temporanea moderata. Tuttavia, quando vengono eseguite manovre di commutazione monofase, o quando si verificano aperture asimmetriche e guasti transitori, l'interazione tra la capacità distribuita della linea e l'induttanza non lineare del reattore di derivazione crea un circuito risonante altamente complesso.

Dalla prospettiva della teoria dei campi elettromagnetici e dei circuiti accoppiati, l'induttanza del reattore cessa di essere un parametro invariante nel tempo per diventare una funzione dipendente dal flusso magnetico concatenato \lambda e dalla corrente ii. Il comportamento del nucleo di ferro è governato dall'equazione costitutiva:

i(t)=f(λ(t))=n=1,3,5cnλ(t)ni(t) = f(\lambda(t) = \sum_{n=1, 3, 5}^{\infty} c_n \lambda(t)^n

dove i coefficienti cnc_n riflettono la severa non linearità dovuta alla saturazione magnetica. Durante un'apertura monofase, i poli sani della linea di trasmissione continuano a indurre tensioni elettrostatiche ed elettromagnetiche sulla fase disconnessa attraverso le capacità mutue interfase CmCm e le capacità verso terra CsCs. Se l'induttanza del reattore di derivazione LL forma un anello risonante con la capacità equivalente della fase aperta CeqCeq, le condizioni per la ferroresonanza serie o parallelo vengono pienamente stabilite, scatenando sovratensioni sostenute a bassa frequenza, armoniche di ordine sub-sincrono e correnti di eccitazione che superano di diversi ordini di grandezza la corrente nominale dell'apparecchio.

Equazioni Differenziali e Modellazione Non Lineare del Nucleo Saturabile

La modellazione matematica rigorosa della ferroresonanza nei reattori di derivazione richiede la formulazione di sistemi di equazioni differenziali ordinarie e parziali che accoppino i fenomeni elettrostatici della rete con l'isteresi e la saturazione del nucleo magnetico. Considerando un circuito equivalente monofase semplificato che rappresenta una fase disconnessa di una linea di trasmissione con un reattore di derivazione connesso ai terminali, la dinamica del sistema è descritta dalla seguente equazione integro-differenziale:

d2λdt2+RL0dλdt+ω02f(λ)=ω02Vmsin(ωt+θ)\frac{d^2 \lambda}{dt^2} + \frac{R}{L_0} \frac{d\lambda}{dt} + \omega_0^2 f(\lambda) = \omega_0^2 V_m \sin(\omega t + \theta)

Dove \lambda è il flusso concatenato nel reattore, RR è la resistenza equivalente delle perdite nel rame e nel nucleo, L0L_0 è l'induttanza nel regime lineare, \omega_0 è la frequenza fondamentale del sistema, e V_m \sin(\omega t + \theta) rappresenta la tensione indotta dalle fasi adiacenti energizzate attraverso l'accoppiamento capacitivo. La funzione di non linearità del nucleo magnetico f(\lambda) è comunemente modellata mediante l'approssimazione di Fröhlich modificata o tramite curve di saturazione analitiche basate su funzioni arcotangente iperbolica:

i(λ)=I0[λλs+Msinh(N(λλs)k)]i(\lambda) = I_0 \left[ \frac{\lambda}{\lambda_s} + M \sinh\left( N \left(\frac{\lambda}{\lambda_s}\right)^k \right) \right]

In questo sistema, la presenza della saturazione fa sì che l'induttanza apparente del reattore varí dinamicamente durante ogni ciclo della tensione applicata. Quando il nucleo entra nella regione di saturazione profonda, l'induttanza differenziale del reattore Ldiff = d\lambda / di diminuisce drasticamente, avvicinandosi ai valori dell'induttanza di dispersione dell'avvolgimento. Questa improvvisa riduzione sposta la frequenza di risonanza del circuito verso valori vicini o inferiori alla frequenza fondamentale, inducendo modi di oscillazione subarmonici (tipicamente di ordine 1/2, 1/3 o 1/5) e quasi-periodici.

L'analisi di stabilità di Lyapunov e i diagrammi di biforcazione dimostrano che il sistema dinamico del reattore ferroresonante possiede molteplici stati stabili (multi-stabilità). A seconda delle condizioni iniziali (come l'angolo di disconnessione dell'interruttore automatico e il valore istantaneo del flusso rimanente \lambda_r nel nucleo), il sistema può convergere verso uno stato normale a bassa tensione o saltare improvvisamente a uno stato ferroresonante caratterizzato da correnti distorte e flussi magnetici elevati che superano il punto di ginocchio della curva di saturazione di oltre il 200%.

Topologia del Fenomeno durante le Manovre di Commutazione Monofase

Le manovre di apertura e chiusura unipolare negli interruttori ad alta tensione sono procedure operative standard per l'eliminazione di guasti transitori nelle linee di trasmissione. Tuttavia, durante l'intervallo di apertura monofase, la fase interessata rimane isolata dalla fonte di alimentazione principale, ma permane accoppiata elettromagneticamente al resto del sistema energizzato. Le capacità mutue tra le fasi Cab,Cbc,CcaCab, Cbc, Cca e le capacità verso terra Can,Cbn,CcnCan, Cbn, Ccn formano una rete capacitiva complessa che alimenta la fase aperta con una tensione di origine elettrostatica.

Quando un reattore di derivazione è permanentemente connesso ai terminali della linea, l'induttanza di questo reattore viene a trovarsi in parallelo con la capacità equivalente di sequenza zero e di sequenza inversa vista dai terminali aperti della fase. Il circuito equivalente di Thévenin visto dai terminali del reattore consiste in una fonte di tensione equivalente VthVth pilotata dalle tensioni delle fasi sane attraverso i condensatori di accoppiamento mutuo, e un'impedenza equivalente ZthZth dominata dalla capacità equivalente CeqCeq.

Parametro di Manovra Condizione Normale / Operativa Condizione Critica di Ferroresonanza Conseguenza Elettromeccanica
Stato dell'Interruttore Chiusura trifase simmetrica Apertura monofase asimmetrica Severo squilibrio di tensioni e correnti di sequenza omopolare.
Capacità Equivalente (CeqCeq) Bassa impedenza capacitiva di linea Accoppiamento ottimale con LreactorLreactor Sintonizzazione del circuito risonante su frequenze subarmoniche.
Flusso Rimanente (\lambda_r) Basso o distribuito simmetricamente Elevato flusso rimanente post-interruzione Asimmetria nel ciclo di isteresi e immediata entrata in saturazione.
Perdite di Sistema (ReqReq) Normale smorzamento per carico Basso smorzamento a vuoto / carico leggero Incapacità del sistema di dissipare l'energia oscillatoria accumulata.

Durante l'apertura del polo di una fase, l'arco elettrico transitorio nell'interruttore può generare riaccensioni ad alta frequenza (re-strikes), iniettando componenti transitorie ad alta frequenza che eccitano i modi naturali dell'insieme reattore-linea. Se il flusso rimanente nella colonna del reattore corrispondente alla fase aperta è elevato e coincide in polarità con il semiciclo iniziale della tensione indotta, il nucleo raggiunge la saturazione magnetica profonda nel primo mezzo ciclo successivo all'estinzione dell'arco, avviando la traiettoria verso l'instabilità ferroresonante.

Analisi Forense dei Guasti e Conseguenze sugli Asset delle Sottostazioni

Gli eventi di ferroresonanza e di estrema saturazione magnetica nei reattori di derivazione generano severe sollecitazioni termiche, dielettriche e meccaniche che compromettono l'integrità degli asset della sottostazione. L'analisi forense dei guasti verificatisi nei sistemi EHV rivela tipici pattern di danno nei diversi componenti del sistema elettrico.

Deformazione Meccanica e Sforzi Elettrodinamici negli Avvolgimenti

Le correnti di eccitazione ferroresonante sono composte massicciamente da armoniche di ordine dispari (principalmente terza e quinta armonica) e componenti subarmoniche. Queste correnti possono raggiungere magnitudini da 3 a 5 volte la corrente nominale di progetto del reattore. Gli sforzi elettromeccanici sui conduttori dell'avvolgimento sono proporzionali al quadrato della corrente istantanea (F \propto i(t)^2). Le forze radiali e assiali cicliche inducono fatica meccanica sui supporti di isolamento, provocando l'allentamento dei blocchi di pressione dell'avvolgimento, l'abrasione dell'isolamento in carta kraft e, infine, cortocircuiti tra le spire.

Riscaldamento Supplementare per Correnti Parassite e Perdite nel Nucleo

Il funzionamento del nucleo magnetico nella regione di saturazione profonda altera la distribuzione del flusso magnetico, costringendo il flusso a fuoriuscire dalle lamiere di acciaio al silicio a grano orientato per dirigersi verso le strutture di serraggio, gli schermi magnetici e il serbatoio metallico del reattore. Le perdite per correnti parassite (Foucault) sulle piastre di schermatura e sulle pareti del serbatoio aumentano in modo non lineare:

Peddy=kfd2f2Bmax2ρPeddy = \frac{k_f \cdot d^2 \cdot f^2 \cdot Bmax^2}{\rho}

Dove dd è lo spessore delle laminazioni, ff è la frequenza effettiva (inclusi gli armonici), e BmaxBmax è la densità di flusso magnetico di picco. Questo incremento massiccio delle perdite genera punti caldi locali (hot spots) nell'olio isolante, superando i limiti termici ammissibili secondo le norme IEEE C57.109 e IEC 60076-6, provocando la degradazione termica dell'isolamento cellulosico, la generazione di gas disciolti nell'olio (idrogeno, metano, acetilene) e l'intervento delle protezioni a pressione improvvisa (Buchholz).

Sovratensioni Dielettriche e Collasso dell'Isolamento

Le tensioni ferroresonanti sostenute possono raggiungere valori compresi tra 2,5 e 3,5 p.u. con forme d'onda altamente distorte che presentano fronti d'onda ripidi. Queste sovratensioni agiscono continuativamente sul sistema di isolamento principale (olio-carta o isolamento solido a secco). I))&&(I scaricatori di sovratensione a ossido metallico (MOA) situati ai terminali del reattore sono sottoposti a un assorbimento di energia termica eccessivo nel tentativo di limitare tali sovratensioni a bassa frequenza, il che può portare al cedimento termico e all'esplosione catastrofica degli scaricatori qualora la loro capacità di dissipazione di energia (classe di scarica della linea) venga superata.

Strategie Avanzate di Mitigazione e Criteri di Progettazione

La mitigazione efficace della ferroresonanza e della saturazione magnetica nei reattori di derivazione richiede l'implementazione di contromisure a livello di progetto elettromagnetico, selezione delle apparecchiature di manovra e configurazione dei sistemi di protezione avanzata.

Ottimizzazione del Progetto del Nucleo e dei Traferri

Per evitare la rapida entrata in saturazione in condizioni di sovratensione temporanea e accoppiamento capacitivo, i reattori moderni sono progettati con nuclei blindati o a colonne dotati di molteplici traferri non magnetici distribuiti uniformemente lungo l'altezza delle colonne. Ciò linearizza la caratteristica di magnetizzazione fino al 140% - 160% della tensione nominale. L'induttanza linearizzata è espressa da:

L=μ0N2Alcore/μr+lgapμ0N2AlgapL = \frac{\mu_0 N^2 A}{\sum lcore/\mu_r + \sum lgap} \approx \frac{\mu_0 N^2 A}{\sum lgap}

Facendo in modo che la lunghezza totale dei traferri \sum lgap sia il termine dominante al denominatore, la riluttanza magnetica totale diventa indipendente dalla permeabilità relativa dell'acciaio \mu_r, sopprimendo la non linearità e prevenendo la ferroresonanza.

Implementazione di Resistenze di Pre-inserzione e Sincronizzazione degli Interruttori

L'utilizzo di interruttori automatici equipaggiati con meccanismi di controllo sintonizzati sull'angolo di fase (Point-on-Wave Switching) consente di eliminare i transitori severi di inserzione e disinserzione. Aprendo o chiudendo i poli dell'interruttore in corrispondenza degli zeri di corrente o negli istanti ottimali di tensione, viene minimizzata la generazione di flussi rimanenti asimmetrici. Inoltre, l'uso di resistenze di pre-inserzione negli interruttori EHV smorza le oscillazioni transitorie ad alta e media frequenza durante le manovre di commutazione.

Schemi di Smorzamento e Filtri Passivi

Nelle sottostazioni in cui il rischio ferroresonante è strutturalmente elevato a causa di topologie di rete complesse con linee lunghe a vuoto, vengono implementate reti di smorzamento basate su circuiti ad alta resistenza collegati in serie con banchi di condensatori ausiliari, oppure mediante l'aggiunta di avvolgimenti terziari nei reattori collegati in delta chiusa con una resistenza di smorzamento esterna. Questa resistenza dissipa l'energia oscillatoria dei modi subarmonici, impedendo la crescita dell'ampiezza dell'oscillazione ferroresonante.

Applicazione Pratica e Analisi con la Vexten Suite

Per illustrare il dimensionamento e la validazione tecnica di un sistema di compensazione reattiva soggetto a regimi transitori e armonici, viene presentato un caso studio sviluppato tramite i moduli analitici della piattaforma Vexten Suite, applicando gli standard internazionali IEC 60909, IEEE 141 e IEC 60287.

Specifiche del Sistema di Potenza

  • Tensione nominale del sistema (UnU_n): 400 kV (Rete trifase, 50 Hz).
  • Potenza di cortocircuito trifase della rete (Sk3Sk3): 25 GVA.
  • Reattore di derivazione trifase: 150 MVAR, 400/√3 kV per fase, collegamento a stella con neutro rigidamente collegato a terra.
  • Capacitanza equivalente di sequenza omopolare della linea (C0C_0): 12,5 µF per fase.
  • Resistenza equivalente dell'avvolgimento del reattore (RcuRcu): 0,15 \Omega per fase a 75 °C.

Calcolo di Cortocircuito e Parametri Transitori (IEC 60909 / IEEE 141)

L'impedenza equivalente della rete vista dai terminali del reattore viene calcolata a partire dalla potenza di cortocircuito:

Znet=Un2Sk3=(400×103)225×109=6,4 ΩZnet = \frac{U_n^2}{Sk3} = \frac{(400 \times 10^3)^2}{25 \times 10^9} = 6,4 \ \Omega

La reattanza induttiva equivalente della rete (XnetXnet) assumendo un rapporto X/R=20X/R = 20:

Rnet=Znet1+(X/R)2=6,41+400=0,319 ΩRnet = \frac{Znet}{\sqrt{1 + (X/R)^2}} = \frac{6,4}{\sqrt{1 + 400}} = 0,319 \ \Omega
Xnet=Rnet×20=0,319×20=6,38 ΩXnet = Rnet \times 20 = 0,319 \times 20 = 6,38 \ \Omega

L'induttanza nominale del reattore di derivazione per fase è determinata tramite la reattanza nominale dell'apparecchio:

XL=(Un/3)2Q3/3=(400/3)250=53333,3350=1066,67 ΩX_L = \frac{(U_n / \sqrt{3})^2}{Q_3 / 3} = \frac{(400 / \sqrt{3})^2}{50} = \frac{53333,33}{50} = 1066,67 \ \Omega
L=XL2πf=1066,672π(50)=3,395 HL = \frac{X_L}{2 \pi f} = \frac{1066,67}{2 \pi (50)} = 3,395 \ H

Valutazione del Fattore di Frequenza di Risonanza

Durante l'apertura monofase, la frequenza di risonanza del circuito formato dall'induttanza del reattore LL e dalla capacità equivalente della fase aperta Ceq \approx C_0 viene calcolata tramite:

fres=12πLCeq=12π3,395×(12,5×106)fres = \frac{1}{2 \pi \sqrt{L \cdot Ceq}} = \frac{1}{2 \pi \sqrt{3,395 \times (12,5 \times 10^{-6})}}
fres=12π4,2437×105=12π(6,514×103)=24,45 Hzfres = \frac{1}{2 \pi \sqrt{4,2437 \times 10^{-5}}} = \frac{1}{2 \pi (6,514 \times 10^{-3})} = 24,45 \ Hz

Il valore di fres = 24,45 \ Hz si trova esattamente nella regione sub-sincrona (circa la metà della frequenza fondamentale di 50 Hz), confermando la vulnerabilità del sistema nei confronti dell'eccitazione di subarmoniche di ordine 1/2 quando si verificano condizioni di commutazione asimmetrica.

Derating Armonico e Dimensionamento dei Conduttori (IEC 60287 / NEC 310)

Come conseguenza della saturazione magnetica e dei regimi transitori ferroresonanti, le correnti che circolano nei cavi di connessione e negli avvolgimenti contengono un significativo contenuto armonico. La Vexten Suite valuta il fattore di perdite supplementari e applica il coefficiente di riduzione per armoniche (Harmonic Derating Factor - HDF) in conformità alla norma IEC 60287:

Imax_derated=Irated1+Δ11+h=2nMh2(1+Yh)I_{max\_derated} = Irated \cdot \sqrt{\frac{1 + \Delta_1}{1 + \sum_{h=2}^{n} M_h^2 \cdot (1 + Y_h)}}

Dove MhM_h è il tasso di distorsione armonica di ordine hh, YhY_h è il fattore di perdite per correnti parassite (Foucault) nel conduttore per l'armonica hh, e \Delta_1 rappresenta l'incremento delle perdite per isteresi ed Eddy nel nucleo. Per un contenuto di terza armonica del 25% e di quinta armonica del 15% derivanti da un severo evento transitorio, il fattore di riduzione calcolato dallo strumento Vexten Suite indica che la capacità di trasporto di corrente ammissibile del cavo deve essere ridotta del 18,5% per evitare il degrado termico dell'isolamento polimerico (XLPE).

Modulo Vexten Suite Norma di Riferimento Parametro d'Ingresso Principale Risultato Analitico / Criterio di Progetto
Vexten Transient Studio IEC 60076-6 / IEEE C57.109 Matrice di induttanze non lineari e flussi rimanenti Identificazione di punti di biforcazione e sovratensioni ferroresonanti di 2,8 p.u.
Vexten Short-Circuit Pro IEC 60909 / IEEE 141 Potenza di cortocircuito (Sk3Sk3) e impedenze di sequenza Determinazione degli sforzi elettrodinamici di picco e delle correnti asimmetriche.
Vexten Thermal Cable Sizer IEC 60287 / NEC 310 Spettro armonico di corrente e temperatura ambiente Calcolo del fattore HDF e sovradimensionamento termico dei cavi di potenza.
Vexten Protection Coordination IEC 60255 / IEEE C37.112 Caratteristiche di eccitazione e sovracorrente temporanea Taratura dei relè di sovracorrente a tempo e protezione da sovratensione residua.

Conclusioni e Linee Guida Normative

L'instabilità dovuta a ferroresonanza e saturazione magnetica nei reattori di derivazione durante le manovre di commutazione monofase rappresenta uno dei fenomeni elettromagnetici più complessi e distruttivi nei sistemi elettrici di potenza in alta tensione. L'analisi rigorosa richiede di abbandonare i modelli lineari tradizionali e di adottare approcci di simulazione non lineare che incorporino l'isteresi, la saturazione del nucleo e l'accoppiamento elettromagnetico capacitivo distribuito delle linee di trasmissione.

L'implementazione di progetti di nucleo con traferri multipli distribuiti, l'uso di interruttori sincronizzati con controllo dell'angolo di fase e la validazione mediante strumenti di ingegneria avanzata come la Vexten Suite garantiscono che gli impianti di compensazione reattiva operino entro i margini di sicurezza stabiliti dalle normative internazionali IEC ed IEEE. Il rispetto rigoroso di questi criteri previene guasti catastrofici, ottimizza gli investimenti negli asset delle sottostazioni e assicura l'affidabilità e la stabilità operativa delle reti di trasmissione di energia elettrica su larga scala.