Analisi delle frequenze naturali e dei transitori elettromagnetici nelle reti di distribuzione

Fondamenti Fisico-Matematici dei Transitori Elettromagnetici e della Risposta in Frequenza

L'analisi rigorosa dei transitori elettromagnetici e delle frequenze naturali nelle reti di distribuzione di energia elettrica richiede una formulazione matematica avanzata che va oltre l'analisi di regime permanente fasoriale. Le moderne reti di distribuzione, caratterizzate da un'elevata penetrazione di cavi sotterranei a media tensione, trasformatori con nuclei ferromagnetici non lineari e carichi elettronici decentralizzati, agiscono come sistemi a parametri sia distribuiti che concentrati, soggetti a severe eccitazioni impulsive e armoniche.

Per modellare la propagazione delle onde elettromagnetiche lungo i conduttori della rete, si devono impiegare le equazioni differenziali parziali accoppiate dei telegrafisti (telegrapher's equations). Considerando una linea multiconductora simmetrica o asimmetrica di lunghezza infinitesimale $dx$ , il comportamento istantaneo delle tensioni v(x,t) e delle correnti i(x,t) è definito dal seguente sistema matriciale:

\frac{\partial}{\partial x} \mathbf{v}(x,t) + \mathbf{R} \mathbf{i}(x,t) + \mathbf{L} \frac{\partial}{\partial t} \mathbf{i}(x,t) = \mathbf{0}

\frac{\partial}{\partial x} \mathbf{i}(x,t) + \mathbf{G} \mathbf{v}(x,t) + \mathbf{C} \frac{\partial}{\partial t} \mathbf{v}(x,t) = \mathbf{0}

Dove $\mathbf{R}$ , $\mathbf{L}$ , $\mathbf{C}$ e $\mathbf{G}$ rappresentano rispettivamente le matrici per unità di lunghezza di resistenza, induttanza, capacità e conduttanza trasversale (perdita dielettrica). Trasformando queste equazioni nel dominio della frequenza mediante la trasformata di Laplace o di Fourier, si ottiene l'equazione di propagazione modale:

\frac{d^2}{dx^2} \mathbf{V}(x,s) = \mathbf{Z}(s) \mathbf{Y}(s) \mathbf{V}(x,s) = \mathbf{P}^2(s) \mathbf{V}(x,s)

Dove Z(s) = R + sL, Y(s) = G + sC, e P(s) è la matrice di propagazione. Gli autovalori della matrice prodotto Z(s)Y(s) determinano le costanti di propagazione modale \gamma_k(\omega) = \alpha_k(\omega) + j\beta_k(\omega), le quali governano l'attenuazione e la velocità di fase delle onde viaggianti. Le frequenze naturali del sistema (\omega_n) corrispondono ai poli della funzione di trasferimento dell'impedenza d'ingresso Zin(s) vista da qualsiasi nodo della rete di distribuzione, risolvendo la condizione caratteristica:

\det \left[ \mathbf{I} - \mathbf{M}_{end}(s) \right] = 0

Dove M_{end}(s) rappresenta la matrice di riflessione e trasmissione ai terminali della rete, accoppiando le impedenze dei trasformatori, dei carichi e dei dispositivi di protezione contro le sovratensioni (SPD).

Classificazione e Origine dei Transitori Elettromagnetici nelle Reti in Media e Bassa Tensione

I transitori elettromagnetici nelle reti di distribuzione traggono origine da perturbazioni operative e atmosferiche che iniettano spettri di frequenza a banda larga. La corretta categorizzazione di questi fenomeni è imperativa per la progettazione di schemi di isolamento coordinati ai sensi delle normative IEEE Std 1313.2 e IEC 60071-2.

Transitori di Origine Atmosferica (Fulmini)

Le scariche atmosferiche dirette o indirette sulle linee aeree di distribuzione inducono onde di sovratensione con fronti d'onda estremamente rapidi (dell'ordine da 0.5 \mu s a 5 \mu s) e code lunghe (fino a 50 \mu s). La corrente di fulmo viene tipicamente modellata mediante la doppia funzione esponenziale di Heidler:

i(t) = \frac{I_p}{\eta} \cdot \frac{(t/\tau_1)^n}{1 + (t/\tau_1)^n} \cdot e^{-t/\tau_2}

Dove $I_p$ è l'ampiezza della corrente di cresta, \tau_1 è la costante di tempo del fronte d'onda, \tau_2 è la costante di tempo della coda d'onda, $n$ è il fattore di ripidità (steepness), ed \eta è un fattore di correzione dell'ampiezza. Nel momento in cui queste onde penetrano da sezioni aeree a cavi sotterranei, l'impedenza caratteristica diminuisce drasticamente (tipicamente da Z_c \approx 400\,\Omega a Z_c \approx 30\,\Omega), provocando una rifrazione di tensione che può raddoppiare l'ampiezza dell'onda incidente nel punto di transizione se la lunghezza del cavo è inferiore alla lunghezza critica del fronte d'onda.

Manovre di Apertura e Chiusura degli Interruttori

L'azionamento di interruttori in vuoto (vacuum circuit breakers) o di sezionatori sotto carico genera transitori ad alta frequenza a causa del fenomeno della rignizione multipla dell'arco elettrico (multiple restrikes). Quando i contatti si separano, la corrente attraversa uno zero ad alta frequenza e la rigidità dielettrica del mezzo di interruzione si ripristina a una velocità superiore rispetto alla velocità di salita della tensione di recupero transitoria (TRV). Ciò provoca un troncamento abrupto della corrente (current chopping), la cui energia induttiva immagazzinata viene trasferita alla capacità parassita della rete secondo la relazione energetica:

\frac{1}{2} Leq ichop^2 = \frac{1}{2} Ceq vmax^2 \implies vmax = ichop \sqrt{\frac{Leq}{Ceq}}

Questa sovratensione oscillatoria ad alta frequenza (tipicamente compresa tra 1 kHz e 100 kHz) esercita una severa sollecitazione dielettrica sull'isolamento principale degli avvolgimenti dei trasformatori.

Risonanza e Ferroresonanza nei Sistemi di Distribuzione

La risonanza armonica e la ferroresonanza rappresentano due delle modalità di guasto più distruttive e complesse nelle reti di distribuzione interrate e miste. La loro analisi richiede la valutazione dell'interazione tra gli elementi lineari reattivi e le non linearità ferromagnetiche.

Risonanza Serie e Parallelo

La risonanza parallelo (antirisonanza) si verifica quando la reattanza induttiva equivalente della rete (fornita dai trasformatori di distribuzione e dalle reattanze limitatrici) eguaglia la reattanza capacitiva dei cavi sotterranei e dei banchi di condensatori per la rifasatura a una specifica frequenza armonica $h$ :

\omega_r = \frac{1}{\sqrt{Leq Ceq}} \implies h_r = \frac{f_r}{f_1} = \sqrt{\frac{X_c}{X_L}}

In questo stato, l'impedenza vista dalla fonte raggiunge un valore massimo teorico limitato esclusivamente dalla resistenza equivalente della rete. Qualsiasi iniezione di correnti armoniche (provenienti da carichi non lineari quali azionamenti a velocità variabile e raddrizzatori) vicina a $h_r$ genera massicce cadute di tensione armoniche, surriscaldamento per effetto Joule e grave deformazione dell'onda di tensione.

Ferroresonanza nei Trasformatori

La ferroresonanza è un fenomeno non lineare caratterizzato da sovratensioni sostenute a bassa frequenza (fondamentale, subarmonica o caotica), estrema distorsione armonica ed elevate correnti. Si verifica tipicamente quando un trasformatore con nucleo in ferro non saturabile viene energizzato a vuoto attraverso una capacità in serie (per esempio, la capacità di un cavo sotterraneo disconnesso a un'estremità o fasi aperte a causa dell'operazione monopolare dei fusibili).

Il comportamento non lineare del nucleo magnetico viene modellato mediante il flusso concatenato \lambda(i) o la relazione costitutiva di induttanza non lineare:

v(t) = R i(t) + \frac{d\lambda(i)}{dt}

La coesistenza di molteplici stati stabili (regime fondamentale e regimi ferroresonantici subarmonici) si spiega mediante la teoria delle biforcazioni nei sistemi dinamici non lineari. L'equazione differenziale che governa il circuito LC non lineare eccitato da una sorgente sinusoidale V_m \sin(\omega t) è:

L(i) \frac{d^2 i}{dt^2} + R \frac{d i}{dt} + \frac{1}{C} \ int i \, dt = V_m \omega \cos(\omega t)

Dove l'induttanza L(i) = \frac{d\lambda}{di} varia drasticamente quando si entra nella regione di saturazione magnetica, riducendo il valore di $L$ e spostando la frequenza di risonanza del circuito verso valori che coincidono con la frequenza di rete o con i suoi subarmonici.

Analisi Forense dei Guasti e Conseguenze Dielettriche

I transitori elettromagnetici ad alta frequenza e le risonanze prolungate deteriorano in modo cumulativo e catastrofico gli asset critici della rete di distribuzione. Di seguito viene dettagliata l'analisi forense dei meccanismi di guasto suddivisi per componente:

Trasformatori di Distribuzione e Potenza

Quando un fronte d'onda rapido (come quello generato da fulmini o manovre) penetra nell'avvolgimento in alta tensione di un trasformatore, la distribuzione iniziale della tensione lungo le spire non è lineare. A causa della presenza della capacità serie tra le spire ( $C_s$ ) e della capacità verso terra dell'avvolgimento ( $C_g$ ), il fattore di distribuzione \alpha = \sqrt{C_g / C_s} determina un gradiente di tensione estremamente elevato sui primi dischi o sulle prime spire dell'avvolgimento.

Ciò genera scariche parziali interne (Partial Discharges - PD) nella carta impregnata d'olio, perforazione dell'isolamento cellulosico e cortocircuiti tra le spire. Termicamente, le correnti armoniche indotte dalla risonanza parallelo provocano perdite supplementari negli avvolgimenti dovute all'effetto pellicolare (skin effect) e di prossimità, calcolate mediante il fattore di perdite addizionali $H$ :

PLL = PDC \left[ 1 + \sum_{h=2}^{n} \left( \frac{I_h}{I_1} \right)^2 H(h) \right]

Cavi Sotterranei in Media Tensione

I cavi isolati in Polietilene Reticolato (XLPE) sono altamente vulnerabili ai transitori di tensione prolungati e alla ferroresonanza. Le sovratensioni sostenute accelerano il processo di treeing elettrico (alberatura elettrica), un meccanismo di degradazione dell'isolamento solido in cui le microcavità sottoposte a campi elettrici intensi (> 10 kV/mm) generano scariche parziali continue che espandono canali dendritici fino alla rottura dielettrica totale (cedimento dell'isolamento primario).

Dispositivi di Manovra e Protezioni

Gli interruttori in vuoto e i richiusori subiscono una grave erosione dei contatti principali a causa dell'energia accumulata negli archi ad alta frequenza durante le operazioni di interruzione con carichi induttivi o capacitivi disadattati. Analogamente, gli scaricatori di sovratensione a ossido di zinco (MOV) possono incorrere in un guasto termico per valanga energetica se la densità di energia assorbita durante un transitorio supera la loro capacità di dissipazione nominale $Wth$ .

Parametri di Progettazione, Limiti Normativi e Strategie di Mitigazione

Per garantire l'affidabilità e il coordinamento dell'isolamento secondo le direttive internazionali, gli ingegneri devono applicare rigorosi criteri di progettazione preventiva e mitigare gli effetti delle frequenze naturali.

La seguente tabella consolida i parametri elettrici critici, i limiti normativi secondo gli standard IEEE e IEC, le condizioni di guasto critico e le conseguenze operative associate:

Parametro Elettrico / Fenomeno	Standard di Riferimento (IEEE / IEC)	Limite Normativo / Soglia Critica	Condizione di Guasto Critico	Conseguenze Operative e Dielettriche
Sovratensioni Temporanee (TOV)	IEEE Std 1313.2 / IEC 60071-2	VTOV \le 1.7 \times V_{ph-n} (per t \le 1\, s )	Ferroresonanza sostenuta o rigetto di carico in reti deboli	Degradazione termica e invecchiamento precoce dell'isolamento XLPE nei cavi.
Distorsione Armonica Totale di Tensione (THD-V)	IEEE Std 519-2022 / IEC 61000-2-2	THD-V \le 5.0\% (Sistemi in Media Tensione)	Risonanza parallelo vicina alla 5ª, 7ª o 11ª armonica	Surriscaldamento dei trasformatori, scatti errati delle protezioni e fatica dielettrica.
Tasso di Salita della Tensione (ripidità differenziata $dv/dt$ )	IEC 60034-18-41 / NEMA MG1	dv/dt \le 5\, kV /\mu s (ai terminali di motori/cavi)	Operazione di interruttori in vuoto con rignizioni multiple	Perforazione dell'isolamento delle spire e guasti inter-spira nei trasformatori e nelle macchine rotanti.
Capacità di Assorbimento Energetico negli SPD (MOV)	IEEE Std C62.11 / IEC 61643-11	Assorbimento nominale \ge 2.5\, kJ/kV di tensione assegnata	Impatti multipli di fulmini o fronti d'onda riflessi amplificati	Guasto termico catastrofico (esplosione) dello scaricatore e perdita di protezione.
Fattore di Derating per Armoniche nei Cavi	IEC 60287 / NEC 310.15	Fattore di riduzione $Kharm$ in base allo spettro di corrente	Correnti di neutro superiori a quelle di fase per componenti omopolari (3ª, 9ª)	Fusione dell'isolamento per eccessiva temperatura sul conduttore (> 90°C per XLPE).

Strategie di Mitigazione e Compensazione Attiva/Passiva

Per mitigare le frequenze naturali di risonanza e i transitori severi, vengono implementate le seguenti contromisure ingegneristiche:

Installazione di Filtri Armonici Sintonizzati (Filtri Passivi LC): Vengono collegati in parallelo alla rete per deviare le armoniche critiche (tipicamente la 5ª e la 7ª) verso resistenze di scarica, modificando l'impedenza equivalente della rete ed evitando la risonanza parallelo. La loro sintonizzazione viene eseguita leggermente al di sotto della frequenza armonica target:
$fsintonizzata = \frac{1}{2\pi \sqrt{Lfilter Cfilter}} < farmonica$
Reattanze di Smorzamento e Resistenze di Pre-inserzione: Negli interruttori ad alta tensione e in media tensione critica, l'uso di resistenze di pre-inserzione smorza il fronte d'onda transitorio, riducendo il fattore di sovratensione a valori $vmax / Vph < 2.0$ .
Soppressori di Sovratensione a Ossido di Zinco (MOV) Posizionati Strategicamente: Installati nei punti di transizione tra linee aeree e cavi sotterranei, nonché sui morsetti dei trasformatori, per tagliare i picchi di tensione transitoria e limitare l'energia incidente.

Applicazione Pratica e Analisi mediante Vexten Suite

La modellazione avanzata delle reti di distribuzione per lo studio dei transitori elettromagnetici e delle frequenze naturali richiede piattaforme di calcolo computazionale ad alta precisione come Vexten Suite. Di seguito vengono dettagliate la metodologia algoritmica e le procedure ingegneristiche implementate nel software per ottemperare agli standard internazionali.

Flusso di Calcolo del Cortocircuito e dei Transitori (IEC 60909 / IEEE 141)

Il motore di calcolo transitorio di Vexten Suite risolve il sistema di equazioni differenziali nel dominio del tempo utilizzando il metodo di integrazione numerica trapezoidale con passo temporale adattativo (\Delta t < 10\,\mu s). La procedura analitica integrata comprende le seguenti fasi:

Acquisizione e Modellazione Topologica: Il sistema importa i parametri fisici dei cavi (geometria, strati semiconduttori, isolamento XLPE), dei trasformatori (impedenze di cortocircuito, corrente a vuoto, curva di saturazione di Fröhlich-Kennelly) e delle fonti di generazione distribuita.
Analisi della Frequenza Naturale (Scansione d'Impedenza - Frequency Scan): Viene iniettata una sorgente di corrente armonica ad ampiezza unitaria e frequenza variabile f \in [10\, Hz , 100\, kHz ] nei nodi critici della rete. Vexten calcola l'impedenza nodale Znn(j\omega) e genera il diagramma di Bode di magnitudine e fase, identificando automaticamente i picchi di risonanza parallelo e le valli di risonanza serie.
Simulazione di Manovre e Fulmini (Motore di tipo EMTP): Vengono configurati eventi di commutazione (apertura di interruttori con specifiche correnti di taglio) e onde di fulmine normalizzate (8/20 \mu s e 1.2/50 \mu s). Il software valuta la propagazione delle onde viaggianti considerando l'attenuazione modale dipendente dalla frequenza e l'effetto pellicolare (skin effect) mediante la formulazione di Carson modificata o le funzioni di Bessel per conduttori cilindrici.
Valutazione del Derating Termico dei Cavi (IEC 60287): Sulla base dello spettro della corrente armonica ottenuto dall'analisi del flusso di carico armonico, Vexten calcola automaticamente il fattore di riduzione della capacità di trasporto di corrente (ampacity) del cavo sotterraneo, assicurando che la temperatura massima del conduttore non superi il limite termico normativo di 90°C per l'XLPE in regime distorto.

In questo modo, Vexten Suite consente all'ingegnere specialista di effettuare una diagnostica completa, prevenendo guasti catastrofici dovuti a fenomeni di risonanza e garantendo il coordinamento ottimale dell'isolamento all'interno di reti di distribuzione complesse.