Misura della resistività apparente del terreno con metodo Wenner a 4 elettrodi e stratificazione a due strati secondo IEEE 81 e IEEE 80

Fondamenti Elettromagnetici della Conduzione nei Mezzi Geoelettrici e Metodo di Wenner

La caratterizzazione geoelettrica del sottosuolo costituisce il pilastro fondamentale e imprescindibile per la progettazione ingegneristica dei sistemi di messa a terra (impianti di terra) in sottostazioni di media, alta ed altissima tensione. A differenza dei conduttori metallici puri, in cui il trasporto di carica è regolato esclusivamente dagli elettroni liberi secondo la legge di Ohm microscopica, il terreno è un mezzo eterogeneo, anisotropo e multifase. La conduzione elettrica nel sottosuolo si sviluppa prevalentemente attraverso un meccanismo elettrolitico all'interno del fluido poroso interstiziale, governato dalla salinità, dall'umidità volumetrica, dalla temperatura, dal grado di compattazione e dalla composizione mineralogica della matrice solida.

A livello macroscopico, il comportamento del campo elettrostatico quasi-stazionario nel semispazio conduttore infinito inferiore ($z \ge 0$) è descritto dalle equazioni di Maxwell per mezzi conduttori in regime stazionario, in cui la densità di corrente $\mathbf{J}$ e il potenziale scalare elettrico $\Phi$ soddisfano:

Per un mezzo omogeneo e isotropo con conducibilità scalare costante $\sigma = 1/\rho$, l'equazione differenziale fondamentale si riduce alla classica equazione di Laplace $\nabla^2 \Phi = 0$, soggetta alle condizioni al contorno di Neumann all'interfaccia terreno-aria ($z = 0$):

Derivazione del Potenziale per Sorgenti Puntiformi e Configurazione Tetrapolare Simmetrica

Considerando una sorgente puntiforme di corrente $I$ iniettata sulla superficie di un semispazio omogeneo ($z = 0$), la simmetria sferica nel semispazio impone che le linee di corrente divergano radialmente su un angolo solido di $2\pi$ steradianti. Il potenziale elettrico $\Phi(r)$ a una distanza radiale $r$ dalla sorgente si ottiene integrando la legge di Ohm in coordinate sferiche:

Nel metodo tetrapolare proposto da Frank Wenner (1915), vengono disposti quattro picchetti collineari ed equidistanti con interasse pari ad $a$. I due elettrodi esterni ($C_1$ e $C_2$) fungono da iniettori di corrente ($+I$ e $-I$), mentre i due elettrodi interni ($P_1$ e $P_2$) operano come sonde di misura per rilevare la differenza di potenziale $\Delta V = \Phi(P_1) - \Phi(P_2)$.

Applicando il principio di sovrapposizione degli effetti per le sorgenti di corrente posizionate in $C_1$ ($x = 0$) e $C_2$ ($x = 3a$), i potenziali in corrispondenza degli elettrodi voltmetrici ubicati in $P_1$ ($x = a$) e $P_2$ ($x = 2a$) risultano:

La differenza di potenziale netta misurata dal voltmetro ad alta impedenza è data da:

Esplicitando il parametro di resistività, si definisce la resistività apparente $\rho_a$ per una spaziatura interelettrodica $a$:

dove $R_m = \Delta V / I$ rappresenta la resistenza di trasferimento mutuo apparente rilevata dal tellurometro digitale a bassa frequenza.

Correzione per la Profondità Finita di Infissione dei Picchetti

Nelle reali condizioni di campo, i picchetti di prova vengono infissi fino a una profondità finita $b$. Qualora il rapporto $b/a$ non risulti trascurabile ($b > 0.1a$), l'approssimazione puntiforme di superficie introduce un errore sistematico non trascurabile. Considerando elettrodi cilindrici verticali di lunghezza $b$, la formulazione matematica rigorosa, corretta per l'effetto di penetrazione dei picchetti, è espressa come segue:

Per $b \to 0$, il denominatore della frazione converge esattamente al valore $2$, restituendo la classica equazione di Wenner $\rho_a = 2\pi a R_m$. Nella prassi operativa conforme alla norma IEEE Std 81, si raccomanda di mantenere $b \le 0.05a$ al fine di poter trascurare la correzione geometrica mantenendo l'errore analitico al di sotto dell'$1\%$.

Sorgenti di Incertezza Metrologica, Interferenze e Accoppiamenti in Campo

L'acquisizione di profili di resistività apparente ad alta fedeltà è soggetta a molteplici disturbi elettromagnetici e non-idealità ambientali che degradano severamente la funzione di trasferimento del sottosuolo. I fattori più critici sono riassunti di seguito:

Sorgente di Disturbo	Meccanismo Fisico / Elettromagnetico	Effetto su $\rho_a$	Strategia di Mitigazione Normativa
Accoppiamento Induttivo tra i Cavi	Induzione mutua $M \cdot \frac{ d i}{ d t}$ tra i circuiti di corrente $C_1-C_2$ e di potenziale $P_1-P_2$.	Sovrastima artificiale di $\Delta V$, con sfasamento apparente a frequenze superiori a 100 Hz.	Separazione fisica trasversale dei cavi di corrente e tensione; utilizzo di cavi coassiali o bipolari schermati; iniezione a frequenze sub-armoniche (es. 55 Hz, 94 Hz, 105 Hz, 128 Hz).
Elevate Resistenze di Contatto ($R_c$)	Impedenza di interfaccia elettrodo-terreno eccessiva causata da ghiaia asciutta, roccia o scarsa umidità superficiale.	Saturazione dello stadio d'ingresso del voltmetro, drastica attenuazione di $I$, rumore termico Johnson-Nyquist e letture erratiche.	Umettamento dei picchetti con soluzione salina o bentonite; inserimento di picchetti multipli in parallelo sui terminali amperometrici.
Correnti Parassite e Telluriche	Ritorno di neutro squilibrato, trazione elettrica in CC/CA, gradienti galvanici e attività geomagnetica.	Distorsione del segnale differenziale su $P_1-P_2$, deriva dello zero in CC e fluttuazioni caotiche a bassa frequenza.	Filtraggio sincrono con modulazione di fase (rilevamento omodina/Lock-in); iniezione di corrente a onda quadra con inversione periodica di polarità.
Strutture Metalliche Interrate	Condotte di gas/acqua, reti di terra esistenti e schermi di cavi che fungono da percorsi a bassa impedenza.	Sottostima critica della resistività per grandi spaziature ($a > 10 m$), con conseguente distorsione della stratificazione.	Esecuzione di stendimenti ortogonali e diagonali; distanziamento dalle infrastrutture metalliche interrate superiore a $3a_{\max}$.

Teoria Matematica della Stratificazione Geoelettrica a Due Strati (Bistrato)

Nella maggior parte delle formazioni geologiche di interesse per l'ingegneria elettrica, il sottosuolo non può essere modellato con precisione come un semispazio omogeneo. La sedimentazione, il grado di costipamento e la presenza di falde acquifere determinano marcati gradienti verticali. Il modello canonico fondamentale è il modello stratificato a due strati (bistrato), costituito da un primo strato superficiale di spessore finito $h$ e resistività intrinseca $\rho_1$, sovrapposto a un semispazio infinito inferiore di resistività intrinseca $\rho_2$.

Risoluzione dell'Equazione di Laplace con Condizioni al Contorno

Per una sorgente puntiforme che inietta una corrente $I$ nell'origine $(0,0,0)$ sulla superficie ($z=0$), il campo scalare di potenziale espresso in coordinate cilindriche $(\Phi(r,z))$ soddisfa l'equazione di Laplace all'interno di ciascun mezzo omogeneo:

Applicando il metodo di separazione delle variabili mediante la trasformata integrale di Fourier-Bessel (di ordine zero), la soluzione generale per il potenziale nel primo strato ($0 \le z \le h$) e nel secondo strato ($z \ge h$) si formula come segue:

dove $J_0(\lambda r)$ è la funzione di Bessel di prima specie e ordine zero, mentre i coefficienti $A(\lambda)$, $B(\lambda)$ e $C(\lambda)$ vengono determinati imponendo le condizioni al contorno di continuità del potenziale e della componente normale della densità di corrente:

1. Condizione alla superficie aria-terreno ($z = 0$): $\left. \frac{\partial \Phi_1}{\partial z} \right|_{z=0} = 0 \implies A(\lambda) = B(\lambda)$
2. Continuità del potenziale all'interfaccia ($z = h$): $\Phi_1(r, h) = \Phi_2(r, h)$
3. Continuità della densità di corrente normale all'interfaccia ($z = h$): $\frac{1}{\rho_1} \left. \frac{\partial \Phi_1}{\partial z} \right|_{z=h} = \frac{1}{\rho_2} \left. \frac{\partial \Phi_2}{\partial z} \right|_{z=h}$
4. Condizione di regolarità all'infinito: $\lim_{z \to \infty} \Phi_2(r, z) = 0$

Definendo il coefficiente di riflessione della resistività $k$ come:

Risolvendo il sistema algebrico lineare per $A(\lambda)$, si ottiene:

Sviluppando il denominatore in serie geometrica infinita $(1 - k e^{-2\lambda h})^{-1} = \sum_{n=1}^{\infty} k^n e^{-2n\lambda h}$, la funzione di potenziale sulla superficie libera ($z = 0$) assume la forma:

Facendo ricorso all'identità integrale di Lipschitz-Weber:

Si perviene all'espressione analitica in forma chiusa per il potenziale elettrostatico superficiale, riconducibile al metodo delle immagini speculari infinite:

Equazione Analitica della Resistività Apparente di Wenner

Applicando la formulazione generale del potenziale di superficie alla configurazione a quattro punti di Wenner, il potenziale differenziale misurato tra $P_1$ e $P_2$ a fronte delle correnti $+I$ in $C_1$ e $-I$ in $C_2$ fornisce l'equazione cardine per la resistività apparente di un terreno bistrato:

Tale serie alternata converge assolutamente per ogni $|k| < 1$ e consente di calcolare teoricamente il profilo di $\rho_a(a)$ in funzione della spaziatura interelettrodica $a$. Essa esibisce due asintoti fisici fondamentali:

Algoritmi di Inversione Geoelettrica e Fitting Non Lineare delle Curve

Nella pratica dell'ingegneria forense e di progettazione, il problema consiste nel determinare il vettore dei parametri incogniti del terreno $\mathbf{p} = [\rho_1, \rho_2, h]^T$ a partire da un insieme discreto di $M$ misurazioni sperimentali $\{(\ln a_i, \ln \rho_{a,i}^{mis})\}_{i=1}^M$. Questo problema inverso è intrinsecamente non lineare e mal condizionato.

Formulazione del Problema Inverso ai Minimi Quadrati Non Lineari

Si definisce la funzione obiettivo dell'errore quadratico relativo normalizzato $\chi^2(\mathbf{p})$ in scala logaritmica al fine di bilanciare le sensibilità su differenti ordini di grandezza:

L'ottimizzazione viene eseguita mediante l'Algoritmo di Levenberg-Marquardt (Damped Gauss-Newton), il quale interpola in modo adattivo tra il metodo del gradiente discendente e il metodo di Gauss-Newton linearizzato:

dove $\mathbf{J} \in \mathbb{R}^{M \times 3}$ rappresenta la matrice Jacobiana delle sensibilità parziali:

e $\mu_k \ge 0$ è il fattore di smorzamento (damping factor) ricalcolato a ogni iterazione in base al guadagno effettivo nella riduzione del residuo.

Matrice di Sensibilità Analitica (Derivate Parziali di Fréchet)

A differenza del calcolo numerico alle differenze finite, suscettibile a errori di troncamento e arrotondamento, la derivazione analitica esatta del gradiente accelera la convergenza asintotica:

dove $\frac{\partial k}{\partial \rho_1} = \frac{-2\rho_2}{(\rho_1 + \rho_2)^2}$ e $\frac{\partial k}{\partial \rho_2} = \frac{2\rho_1}{(\rho_1 + \rho_2)^2}$.

Il criterio rigoroso di convergenza impone che lo scarto quadratico medio relativo (RMS) residuo soddisfi la condizione:

Impatto della Stratificazione Bistrato sulla Progettazione delle Reti di Terra (IEEE Std 80 / IEC 60479)

Omettere la reale stratificazione del terreno assumendo un modello a resistività costante equivalente calcolata come semplice media aritmetica conduce a errori sistematici gravissimi nel dimensionamento dei conduttori in rame, nella determinazione della lunghezza dei picchetti e nella verifica delle tensioni di passo e di contatto all'interno delle sottostazioni.

Resistenza di Terra ($R_g$) con Modello Bistrato

In accordo con le formulazioni estese di Sverak e Schwarz per mezzi a due strati, la resistenza globale di una maglia di terra quadrata o rettangolare integrata con dispersori verticali dipende in misura critica dallo spessore del primo strato $h$ e dal rapporto di contrasto di resistività $\rho_2 / \rho_1$:

dove $\rho_{eq}$ è la resistività apparente equivalente vista dalla specifica geometria dell'impianto di terra:

dove $L_{ totale } = L_c + L_v$ è lo sviluppo complessivo dei conduttori orizzontali ($L_c$) e dei picchetti verticali ($L_v$), mentre $\psi(h)$ rappresenta la funzione di accoppiamento di profondità.

Potenziali di Passo e Contatto: Confronto Fenomenologico degli Scenari

Parametro di Progetto / Sicurezza	Scenario A: Terreno Favorevole ($\rho_1 > \rho_2$, $k < 0$)	Scenario B: Terreno Critico ($\rho_1 < \rho_2$, $k > 0$)	Implicazione Critica di Sicurezza
Distribuzione della Corrente di Guasto	La corrente iniettata tende a defluire preferenzialmente in senso verticale verso lo strato inferiore profondo a bassa resistività.	La corrente viene riflessa dallo strato inferiore altamente resistivo (roccia), confinandosi e canalizzandosi orizzontalmente nello strato superficiale.	Nello Scenario B, la densità di corrente superficiale $J_s$ subisce un incremento esponenziale.
Potenziale di Rete ($V_m$) e di Contatto ($E_{ touch }$)	Gradienti di potenziale superficiale attenuati e uniformi. $E_{ touch }$ rimane entro i limiti ammissibili con orditure standard.	Elevazione severa del potenziale di terra (GPR/EPR) e picchi di gradiente di tensione estremamente ripidi sulla superficie della maglia.	Nello Scenario B, viene frequentemente oltrepassata la soglia di fibrillazione ventricolare secondo IEEE 80 / IEC 60479.
Efficacia dei Picchetti Verticali di Terra	Altamente efficaci se intersecano lo strato profondo ($l_v > h$), abbattendo sensibilmente la resistenza $R_g$.	Poco efficaci; i dispersori verticali dissipano una quota marginale di corrente nello strato roccioso profondo ad alta resistività.	Nello Scenario B, è imperativo infittire la maglia orizzontale ed estendere l'area perimetrale complessiva.
Fattore di Riduzione dello Strato di Ghiaia ($C_s$)	Lo strato superficiale di pietrisco ($h_s$) fornisce un'efficace mitigazione standard ($C_s \approx 0.7 - 0.85$).	L'interazione tra un $\rho_1$ moderato con il substrato roccioso superficiale ad elevato $\rho_2$ altera il meccanismo di riflessione dello strato di ghiaia.	Richiede il ricalcolo rigoroso di $C_s$ tramite lo sviluppo in serie a riflessione multipla (ghiaia-strato 1-strato 2).

Calcolo Rigoroso del Fattore di Riflessione dello Strato Superficiale ($C_s$)

Il fattore correttivo $C_s$ per lo strato di materiale ad alta resistività (ghiaia/pietrisco $\rho_s \approx 3000\,\Omega\cdot m$) di spessore $h_s$ installato al di sopra di un terreno bistrato non può essere calcolato con la classica formulazione empirica di Sunde valida per terreno omogeneo. La formulazione analitica generalizzata mediante sviluppo in serie dei coefficienti combinati risulta:

dove $k_s = (\rho_1 - \rho_s)/(\rho_1 + \rho_s)$ e $r_d = 0.08 m$ modella il raggio equivalente del piede umano come un disco circolare metallico a contatto con il suolo.

Analisi Forense dei Guasti Associati a una Caratterizzazione Geoelettrica Inadeguata

L'adozione di metodologie errate durante le campagne di prospezione o l'ingiustificata semplificazione a mezzo omogeneo in presenza di formazioni sottostanti fortemente riflettenti ($k \to +1$) sono state la causa primaria di guasti distruttivi in impianti di potenza. Di seguito si illustrano i meccanismi fisici di collasso e la casistica forense documentata:

Perforazione Dielettrica nei Trasformatori di Misura e Servizi Ausiliari dovuta a Sovraelevazione Critica dell'EPR/GPR

In una stazione elettrica a 230 kV realizzata su uno strato superficiale di arenaria ($\rho_1 = 180\,\Omega\cdot m$, $h = 1.8 m$) sovrastante un basamento granitico compatto ($\rho_2 = 4200\,\Omega\cdot m$, $k = +0.918$), il progetto esecutivo originario ipotizzava un terreno omogeneo equivalente di $350\,\Omega\cdot m$. Durante un cortocircuito monofase franco a terra con corrente simmetrica pari a $I_f = 22 kA$, la resistenza reale misurata dell'impianto di terra si è attestata a $2.85\,\Omega$ (contro il valore calcolato di progetto di appena $0.68\,\Omega$).

L'elevata densità di corrente, impossibilitata a penetrare il massiccio granitico profondo, è defluita trasversalmente lungo i cunicoli portacavi dei circuiti di controllo, provocando un'innalzamento del potenziale locale di terra dei quadri di campo di gran lunga superiore alla tenuta dielettrica dei secondari dei trasformatori di corrente (TA) e di tensione (TV). La sovratensione di modo comune ha perforato l'isolamento galvanico da 10 kV dei relè di protezione digitali (IED) nella sala comando, distruggendo i moduli di ingresso/uscita analogici e bloccando l'intervento di riserva temporizzato (ANSI 50/51N). Il tempo di estinzione del guasto si è conseguentemente esteso da 80 ms a 1.2 secondi, culminando nell'esplosione distruttiva del trasformatore di potenza principale.

Degrado Termico e Fusione degli Schermi dei Cavi Sotterranei in Alta Tensione

Nelle linee in cavo interrato a 115 kV con isolamento in XLPE e schermi metallici connessi a terra a entrambe le estremità (solid bonding), la sottostima della resistività degli strati profondi altera sensibilmente il calcolo dell'impedenza mutua dell'anello terra-schermo $\underline{Z}_{m(c-g)}$ secondo le formulazioni di Carson-Clem:

In presenza di uno strato profondo altamente resistivo ($\rho_2 \gg \rho_1$), la profondità di ritorno equivalente $D_e$ subisce una dilatazione superiore al 300%. Ciò determina un drastico incremento della reattanza induttiva del ritorno attraverso il terreno, costringendo la corrente di squilibrio e di guasto a richiudersi quasi integralmente attraverso gli schermi metallici dei cavi di potenza. La densità di corrente termica nello schermo $J_{ screen }$ ha superato la capacità termica limite del rame secondo le norme IEC 60287 e IEC 60949:

comportando la delaminazione termica della guaina semiconduttiva esterna, il rammollimento del polietilene reticolato (XLPE) per superamento della temperatura critica di 250 °C e il conseguente cedimento dielettrico fase-terra per perforazione causata dalla rapida proliferazione elettrotermica di water-treeing.

Flusso Metodologico di Simulazione e Progettazione Avanzata in Vexten Suite

L'ecosistema integrato di ingegneria Vexten Suite incorpora un solutore deterministico multifisico di inversione geoelettrica accoppiato a moduli per il calcolo delle correnti di cortocircuito (IEC 60909 / IEEE 141), dimensionamento di conduttori e cavi (IEC 60287 / NEC) e ottimizzazione avanzata delle reti di terra (IEEE 80 / IEC 60479). Di seguito viene schematizzata l'architettura algoritmica e il protocollo metodologico applicato nei progetti di infrastrutture critiche:

Esempio Pratico di Simulazione e Inversione Numerica nel Vexten Grounding Engine

È stata condotta una campagna di prospezione geoelettrica tetrapolare con metodo di Wenner presso l'area destinata a una nuova sottostazione industriale da 138/13.8 kV. L'acquisizione dei dati ha fornito i seguenti rilievi sperimentali:

Spaziatura $a$ [m]	Resistenza Misurata $R_m$ [$\Omega$]	Resistività Apparente Sperimentale $\rho_a^{mis}$ [$\Omega\cdot m$]	Resistività Apparente Calcolata $\rho_a^{calc}$ [$\Omega\cdot m$]	Errore Residuo Relativo [\%]
1.0	68.435	430.00	428.12	-0.44
2.0	32.786	412.00	418.54	+1.59
4.0	14.928	375.20	379.80	+1.23
8.0	5.431	273.00	269.45	-1.30
16.0	1.581	159.00	157.10	-1.19
32.0	0.507	102.00	103.85	+1.81

Elaborando la matrice dei dati nel modulo di ottimizzazione non lineare di Vexten Grounding, l'algoritmo di Levenberg-Marquardt converge in 7 iterazioni con un errore quadratico medio (RMS) globale dell'1.32%, determinando il seguente vettore dei parametri ottimali del terreno bistrato:

Analisi Ingegneristica e Ottimizzazione del Progetto dell'Impianto di Terra

Poiché il coefficiente di riflessione $k$ assume un valore fortemente negativo ($k = -0.6812$), lo strato inferiore evidenzia una resistività di oltre cinque volte inferiore rispetto a quella dello strato superficiale ($\rho_2 \ll \rho_1$). L'ottimizzazione computazionale implementata in Vexten Grounding sfrutta questa specifica risposta geoelettrica applicando i seguenti criteri progettuali:

1. Disposizione Strategica di Dispersori Profondi: Si stabilisce l'infissione di picchetti in acciaio ramato (copperweld) con lunghezza minima $l_v = 6.0 m$ ai quattro vertici e nei nodi perimetrali più sollecitati. Superando la profondità dello strato di interfaccia ($h = 4.35 m$), i dispersori verticali penetrano per $1.65 m$ all'interno dello strato altamente conduttivo ($\rho_2 = 82.50\,\Omega\cdot m$). Questa configurazione riduce la resistenza totale di terra $R_g$ del 42.7% rispetto a una rete puramente orizzontale di pari sviluppo lineare di corda di rame.
2. Diradamento dei Conduttori Orizzontali di Maglia: In virtù dell'elevata ammettenza dello strato profondo, l'interasse dei conduttori della griglia orizzontale può essere incrementato da una maglia fitta di $5 m \times 5 m$ a una maglia ottimizzata di $10 m \times 10 m$, garantendo un cospicuo risparmio di rame pur mantenendo i potenziali di passo e contatto ampiamente al di sotto delle soglie di sicurezza normate da IEEE Std 80.
3. Integrazione con il Modulo di Cortocircuito Vexten Suite (IEC 60909): L'impedenza complessa del terreno stratificato viene esportata in modo trasparente nel motore di calcolo dei cortocircuiti di Vexten Suite, consentendo il ricalcolo coerente dei fattori di ripartizione della corrente sui conduttori di guardia aerei e sugli schermi dei cavi ($S_f$). Questo assicura che l'effettiva frazione di corrente dispersa a terra ($I_g = S_f \cdot I_k''$) rispecchi fedelmente la risposta elettromagnetica multifase dell'intera rete di trasmissione e distribuzione interconnessa.

L'implementazione di questo flusso di calcolo integrato assicura la massima accuratezza analitica, la salvaguardia dielettrica del personale di esercizio e delle apparecchiature primarie e secondarie, nonché una marcata ottimizzazione tecnico-economica nella realizzazione degli impianti di terra per sottostazioni di trasformazione e centrali di generazione.