Ingeniería Eléctrica

Tensions de Pas et de Toucher dans les Postes selon IEEE 80

El diseño de mallas de puesta a tierra en subestaciones eléctricas no se limita a lograr un valor arbitrario de resistencia inferior a 1 ohmio. Bajo una falla a

Ing. Francisco Ramírez

Fondements Électrophysiologiques et Critères de Fibrillation Ventriculaire

La conception des systèmes de mise à la terre dans les sous-stations électriques à haute et très haute tension a pour objectif premier la préservation de la vie humaine et l'intégrité des actifs électromécaniques. Lors d'un événement de court-circuit à la terre, l'injection massive de courant dans le sol génère des gradients de potentiel de surface qui peuvent soumettre le corps d'un opérateur ou d'un passant à des différences de potentiel critiques.

La réponse physiologique de l'organisme face au passage d'un courant alternatif à fréquence industrielle (50/60 Hz) est dominée par l'induction d'une fibrillation ventriculaire, phénomène au cours duquel les myocytes ventriculaires perdent leur synchronisme contractile, annulant ainsi le débit cardiaque. La probabilité de fibrillation dépend de la magnitude du courant, de la durée de l'exposition, du trajet à travers le corps et de la phase du cycle cardiaque pendant laquelle se produit le choc électrique (notamment la phase vulnérable correspondant à l'onde T sur l'électrocardiogramme).

IB=ktsI_B = \frac{k}{\sqrt{t_s}}

Dans le cadre de la norme IEEE Std 80, Charles Dalziel a établi de manière empirique que l'énergie tolérable par le corps humain avant d'atteindre le seuil de 0,5 % de probabilité de fibrillation ventriculaire est régie par une constante d'énergie de choc SbS_b. Pour une population ayant un poids corporel d'environ 50 kg, la constante k50=0,116k50 = 0,116, tandis que pour 70 kg, k70=0,157k70 = 0,157. Par conséquent, le courant maximal admissible par le corps IBI_B pour une durée de choc tst_s (où 0,03 \le t_s \le 3,0 s ) est défini analytiquement par :

IB,50kg=0,116ts[A],IB,70kg=0,157ts[A]I_{B,50 kg } = \frac{0,116}{\sqrt{t_s}} \quad [ A ], \qquad I_{B,70 kg } = \frac{0,157}{\sqrt{t_s}} \quad [ A ]

D'autre part, la Commission Électrotechnique Internationale, dans son rapport technique IEC/TS 60479-1, adopte une approche probabiliste non linéaire basée sur les études de Biegelmeier. Au lieu de supposer une relation purement inversement proportionnelle à la racine carrée du temps pour toute la fenêtre temporelle, la norme IEC 60479-1 définit des courbes temps/courant (zones AC-1 à AC-4.3). L'impédance du corps humain ZBZ_B selon la norme IEC n'est pas une résistance fixe, mais un réseau résistif-capacitif dépendant de la tension de contact (UTU_T), du trajet du courant et de l'état d'humidité de la peau :

ZB(UT)=Rint+Rp1(UT)1+jωRp1Cp1+Rp2(UT)1+jωRp2Cp2Z_B(U_T) = Rint + \frac{Rp1(U_T)}{1 + j\omega Rp1 Cp1} + \frac{Rp2(U_T)}{1 + j\omega Rp2 Cp2}

Alors que la norme IEEE Std 80 simplifie de manière conservatrice la résistance du corps humain en une valeur constante purement résistive R_B = 1000\,\Omega (représentative du trajet main-pied ou main-main avec un contact ferme), les normes IEC 61936-1 et IEC 60479-1 appliquent des courbes dynamiques où l'impédance effective décroît drastiquement à mesure que la tension de contact dépasse 200 V, atteignant des valeurs asymptotiques proches de la résistance interne du corps (Rint \approx 500 à 750\,\Omega).

Cadre Mathématique et Formulation Rigoureuse selon IEEE Std 80

Circuits Équivalents de Thévenin et Limites de Tolérabilité

Pour déterminer les tensions admissibles de pas (VstepVstep) et de contact (VtouchVtouch), on a recours à un circuit équivalent de Thévenin vu depuis les points de contact du corps avec la surface du sol et les structures mises à la terre. La résistance équivalente du sol vue par les pieds de l'individu est modélisée en considérant chaque pied comme un disque conducteur plan de rayon b=0,08mb = 0,08 m posé sur la surface d'un demi-espace homogène de résistivité apparente \rho_s.

La résistance propre d'un pied en contact avec le sol est donnée par la formulation de Maxwell :

Rfoot=ρs4bρs4(0,08)=3,125ρs3ρsRfoot = \frac{\rho_s}{4b} \approx \frac{\rho_s}{4(0,08)} = 3,125 \rho_s \approx 3\rho_s

Dans la condition de tension de contact, les deux pieds se trouvent en parallèle à une distance pratiquement identique du point de défaut, de sorte que la résistance équivalente est R_{eq,touch} = \frac{Rfoot}{2} = 1,5\rho_s. Pour la condition de tension de pas, les pieds se trouvent en série séparés par une distance normalisée de 1 mètre, ce qui donne R_{eq,step} = 2 Rfoot = 6\rho_s.

En incorporant une couche superficielle de matériau à haute résistivité (typiquement du gravier ou de la pierre concassée d'épaisseur hsh_s et de résistivité \rho_s) sur un sol naturel de résistivité \rho, il convient d'introduire un facteur de réduction ou de déclassement (derating) C_s(h_s, K), qui compense l'effet de réflexion provoqué par la discontinuité diélectrique :

Cs=1+16bπn=1Kn1+(2nhsb)210,09(1ρρs)2hs+0,09C_s = 1 + \frac{16 b}{\pi} \sum_{n=1}^{\infty} \frac{K^n}{\sqrt{1 + \left(2 n \frac{h_s}{b}\right)^2}} \approx 1 - \frac{0,09 \left(1 - \frac{\rho}{\rho_s}\right)}{2 h_s + 0,09}

Où le coefficient de réflexion KK entre le sol naturel et la couche de gravier superficielle s'exprime comme suit :

K=ρρsρ+ρsK = \frac{\rho - \rho_s}{\rho + \rho_s}

En intégrant ces déductions électrodynamiques, les potentiels maximaux tolérables calculés selon la norme IEEE Std 80 pour un opérateur de 50 kg et de 70 kg se traduisent par les équations fondamentales suivantes :

Vtouch,50=(1000+1,5Csρs)0,116ts,Vtouch,70=(1000+1,5Csρs)0,157tsV_{touch,50} = \left(1000 + 1,5 C_s \rho_s\right) \frac{0,116}{\sqrt{t_s}}, \qquad V_{touch,70} = \left(1000 + 1,5 C_s \rho_s\right) \frac{0,157}{\sqrt{t_s}}
Vstep,50=(1000+6,0Csρs)0,116ts,Vstep,70=(1000+6,0Csρs)0,157tsV_{step,50} = \left(1000 + 6,0 C_s \rho_s\right) \frac{0,116}{\sqrt{t_s}}, \qquad V_{step,70} = \left(1000 + 6,0 C_s \rho_s\right) \frac{0,157}{\sqrt{t_s}}

Élévation du Potentiel de Terre (GPR) et Paramètres de Conception de la Maille

Le courant de conception qui se dissipe effectivement à travers la grille vers le sol infini (IGI_G) ne correspond pas au courant de court-circuit symétrique triphasé ou monophasé total aux bornes (IfI_f), mais il est affecté par le facteur de division du courant de défaut (SfS_f) et le facteur de décrémentation (DfD_f) dû à l'asymétrie de la composante de courant continu :

IG=SfDfIfI_G = S_f \cdot D_f \cdot I_f

Le facteur de décrémentation DfD_f quantifie la teneur apériodique du courant de court-circuit en fonction du rapport réactance/résistance (X/RX/R) de l'impédance de Thévenin au point de défaut et de la durée d'élimination tft_f :

Df=1+Tatf(1e2tf/Ta),Ta=XωR=X2πfRD_f = \sqrt{1 + \frac{T_a}{t_f} \left(1 - e^{-2 t_f / T_a}\right)}, \qquad T_a = \frac{X}{\omega R} = \frac{X}{2\pi f R}

L'élévation maximale de potentiel du système par rapport à un point de référence distant est appelée élévation du potentiel de terre (Ground Potential Rise - GPRGPR) et se calcule au moyen de :

GPR=IGRgGPR = I_G \cdot R_g

La résistance du système de mise à la terre (RgR_g) pour des mailles complexes interconnectées avec des piquets ou des électrodes verticales profondes est déterminée avec une haute précision analytique grâce à la formule généralisée de Sverak :

Rg=ρ[1LT+120A(1+11+h20/A)]R_g = \rho \left[ \frac{1}{L_T} + \frac{1}{\sqrt{20 A}} \left( 1 + \frac{1}{1 + h \sqrt{20 / A}} \right) \right]

LTL_T est la longueur totale des conducteurs enterrés (grille horizontale plus électrodes verticales), AA est la surface occupée par la grille en mètres carrés, et hh est la profondeur d'enfouissement des conducteurs du réseau horizontal.

Tension de Maille et Tension de Pas Réelles dans la Sous-Station

La tension de maille réelle (EmE_m), qui représente la condition la plus défavorable de tension de contact au centre du quadrant extérieur le plus critique de la grille, est formulée rigoureusement par le produit de la résistivité apparente du sol (\rho), du facteur géométrique d'espacement KmK_m, du facteur d'irrégularité correctif KiK_i et de la densité linéaire de courant injecté :

Em=ρIGKmKiLME_m = \frac{\rho \cdot I_G \cdot K_m \cdot K_i}{L_M}

Les facteurs sans dimension KmK_m et KiK_i synthétisent la topologie du réseau et sont définis analytiquement comme suit :

Km=12π[ln(D216hd+(D+2h)28Ddh4d)+KiiKhln(8π(2n1))]K_m = \frac{1}{2\pi} \left[ \ln\left( \frac{D^2}{16 h d} + \frac{(D + 2h)^2}{8 D d} - \frac{h}{4d} \right) + \frac{Kii}{K_h} \ln\left( \frac{8}{\pi (2n - 1)} \right) \right]
Ki=0,644+0,148nK_i = 0,644 + 0,148 \cdot n

DD est l'espacement entre conducteurs parallèles, dd est le diamètre du conducteur de la grille, hh est la profondeur d'enfouissement, nn est le facteur géométrique équivalent dérivé des dimensions et de la forme périmétrique (n = n_a \cdot n_b \cdot n_c \cdot n_d), K_h = \sqrt{1 + h/h_0} (avec h0=1mh_0 = 1 m), et Kii=1Kii = 1 pour les grilles avec électrodes périmétriques ou Kii = 1/(2n)^{2/n} pour les grilles sans piquets.

La longueur effective pondérée du conducteur pour le calcul de la tension de maille (LML_M) incorpore l'apport des électrodes verticales de mise à la terre (LrL_r) :

LM=Lc+[1,55+1,22(LrLx2+Ly2)]LRL_M = L_c + \left[ 1,55 + 1,22 \left( \frac{L_r}{\sqrt{L_x^2 + L_y^2}} \right) \right] L_R

De manière analogue, la tension de pas calculée (EsE_s) au périmètre extérieur de la sous-station est évaluée à l'aide du coefficient géométrique de pas KsK_s :

Es=ρIGKsKiLSE_s = \frac{\rho \cdot I_G \cdot K_s \cdot K_i}{L_S}
Ks=1π[12h+1D+h+1D(10,5n2)]K_s = \frac{1}{\pi} \left[ \frac{1}{2h} + \frac{1}{D + h} + \frac{1}{D} \left( 1 - 0,5^{n-2} \right) \right]
LS=0,75Lc+0,85LRL_S = 0,75 \cdot L_c + 0,85 \cdot L_R

Comportement Transitoire et Impact de la Foudre selon la Norme IEC 62305

Impédance d'Impulsion vs Résistance à Fréquence Industrielle

La réponse d'une grille de mise à la terre face à l'injection de courants de type foudre régulés par la série de normes IEC 62305 diffère radicalement du comportement à 50/60 Hz. Les décharges atmosphériques directes ou induites présentent des fronts d'onde ultra-rapides normalisés (forme d'onde 10/350 \mu s pour le premier impact de foudre selon la norme IEC 62305-1, et 8/20 \mu s pour les impacts subséquents), avec des contenus fréquentiels s'étendant dans la plage des centaines de kilohertz aux mégahertz.

À ces fréquences, la réactance inductive distribuée longitudinale (LL') des conducteurs en cuivre ou en acier galvanisé devient le facteur dominant sur la conductance transversale du sol (GG'). Par conséquent, le courant n'a pas le temps physique de se distribuer dans toute la géométrie de la grille, se concentrant à proximité immédiate du point d'impact.

Zimp(t)=u(t)i(t)RgZimp(t) = \frac{u(t)}{i(t)} \neq R_g

La longueur effective du conducteur enterré lel_e, au-delà de laquelle toute extension supplémentaire de l'électrode ne réduit pas l'impédance d'impulsion apparente pour un temps de front T1T_1, s'exprime analytiquement selon la norme IEC 62305-3 comme suit :

le=ρT1μ01,6ρT1l_e = \sqrt{\frac{\rho \cdot T_1}{\mu_0}} \approx 1,6 \cdot \sqrt{\rho \cdot T_1}

Où \mu_0 = 4\pi \times 10^{-7} H/m est la perméabilité magnétique du vide et T1T_1 est le temps de front de l'impulsion en microsecondes (\mu s ).

Tensions de Pas et de Contact Induites par la Foudre

Sous l'action d'un courant impulsif de foudre i(t), le gradient de potentiel radial à la surface du sol génère des tensions transitoires extrêmes qui peuvent provoquer des décharges disruptives dans l'air le long des membres inférieurs du corps humain (arc superficiel ou flashover). La tension de contact impulsive transitoire (Utouch,impU_{touch,imp}) et la tension de pas transitoire (Ustep,impU_{step,imp}) selon la norme IEC 62305-3 Annexe E sont couplées par des composantes inductives et résistives pures :

Utouch,imp(t)=Rimpi(t)+Lgriddi(t)dtU_{touch,imp}(t) = Rimp \cdot i(t) + Lgrid \cdot \frac{di(t)}{dt}
U_{step,imp}(r) = \frac{\rhoeff \cdot i(t)}{2\pi} \left( \frac{1}{r} - \frac{1}{r + \Delta r} \right)

Dans le cas d'impacts de foudre sur la structure de blindage (pointes de Franklin, câbles de garde ou portiques), les dérivées temporelles élevées (di/dt > 100 kA/ \mu s ) génèrent des potentiels inductifs parasites dans les descentes métalliques qui excédent la rigidité diélectrique de l'air, nécessitant des distances de séparation sécuritaires (ss) calculées conformément à la formulation :

skikckmls \ge k_i \cdot \frac{k_c}{k_m} \cdot l

kik_i est le coefficient qui dépend du niveau de protection contre la foudre (LPL I à IV), kck_c est le facteur de répartition géométrique du courant dans les descentes, kmk_m est le coefficient du matériau isolant (air ou solide), et ll est la distance linéaire le long de la descente jusqu'au point d'équipotentialisation le plus proche.

Matrice Comparative Normative : IEEE Std 80 vs IEC 60479-1 / IEC 61936-1 / IEC 62305-3

Paramètre / Critère de Conception IEEE Std 80 (2013) IEC 60479-1 / IEC 61936-1 (2021) IEC 62305-3 (Protection Foudre)
Modèle de Poids Corporel Discret : 50 kg (k=0,116k=0,116) ou 70 kg (k=0,157k=0,157) Statistique de population (Percentile 5 %, 50 %, 95 %) Non applicable (Risque évalué par choc impulsif)
Impédance du Corps (RBR_B / ZBZ_B) Constante : R_B = 1000\,\Omega (purement résistif) Non linéaire : Z_B = f(U_T), 500 à 1000\,\Omega R_B \approx 0\,\Omega (Arc électrique à haute fréquence)
Régime Temporel Critique Fréquence industrielle (50/60 Hz), 0,03 \le t_s \le 3,0 s Fréquence industrielle (50/60 Hz), 0,01 \le t_s \le 10,0 s Transitoire impulsif (10/350\,\mu s , 8/20\,\mu s )
Résistance de la Couche Superficielle Déclassement analytique formel C_s(h_s, K) Facteur de couche superficielle similaire ou empirique Isolation de contrôle : asphalte 5cm5 cm (\ge 100 k \Omega\cdot m )
Limite de Tension de Contact Vtouch = (1000 + 1,5 C_s \rho_s) \frac{k}{\sqrt{t_s}} Défini par la courbe limite U_{v,tol} = f(t_f, Z_B, Rfoot) Utouch \le 100 kV (avec équipotentialisation)
Facteur de Division de Défaut (SfS_f) Calcul explicite avec câbles de garde et neutres Facteur de réduction rr (fonction des écrans et gardes) Facteur de division kck_c basé sur la topologie des descentes
Comportement du Sol Stratification de résistivité multicouche à 50 Hz Modèles bi-couche et demi-espace homogène Effet d'ionisation du sol (E_0 \approx 300 kV/m )

Analyse Forensique des Défauts et Phénomènes de Transfert de Potentiel

Mécanismes de Propagation des Potentiels Transférés

L'un des modes de défaillance les plus dévastateurs dans l'ingénierie forensique des sous-stations correspond à l'exportation ou à l'importation par inadvertance de l'élévation du potentiel de terre (GPRGPR) hors des confins périmétriques de l'installation. Lorsqu'un court-circuit à la terre se produit à l'intérieur de la sous-station, toute la grille de terre s'élève à un potentiel V=GPRV = GPR par rapport au potentiel zéro de référence dans une terre éloignée.

Vtrans=GPRΔVdropIGRgVtrans = GPR - \Delta Vdrop \approx I_G \cdot R_g

Les principaux vecteurs physiques de couplage galvanique des potentiels transférés comprennent :

  • Conducteurs Neutres des Réseaux de Moyenne et Basse Tension : Si le neutre de l'enroulement secondaire du transformateur de services auxiliaires ou d'une ligne de distribution rurale est rigidement mis à la terre sur la grille interne et s'étend hors de la sous-station avec une mise à la terre multiple, le GPRGPR complet s'applique directement sur les isolateurs des réseaux des utilisateurs finaux, ce qui entraîne l'explosion de compteurs résidentiels, la perforation d'appareils électroménagers et un risque d'électrocution massive.
  • Écrans Métalliques des Câbles de Puissance et de Télécommunications : Les écrans des câbles souterrains MT connectés à la terre aux deux extrémités transfèrent le courant de retour du sol vers des terminaux distants, provoquant des courants circulaires de plusieurs centaines d'ampères qui fondent la gaine extérieure en polyéthylène (défaut thermique par I2tI^2 t et induisent des surtensions qui perforent les interfaces d'instrumentation numérique.
  • Conduites Hydrauliques ou de Gaz et Rails de Chemins de Fer : Des éléments métalliques continus enterrés sans joints d'isolation diélectrique agissent comme des électrodes d'extension infinie, exposant les opérateurs dans des stations de pompage ou des passages à niveau éloignés à l'intégralité du GPRGPR, où la tension de contact transférée V_{touch,trans} \approx GPR est dépourvue de l'effet amortisseur du gravier de la sous-station.

Comportement dans les Sous-Stations à Isolation Gazeuse (GIS)

Dans les sous-stations isolées auhexafluorure de soufre (SF6SF _6), la manœuvre de sectionneurs ou la survenue de décharges disruptives génère des transitoires très rapides (VFT, Very Fast Transients) avec des temps de montée de l'ordre de la nanoseconde (tr<10nst_r < 10 ns). Ces fronts d'onde ne parviennent pas à circuler vers la terre à travers les connexions conventionnelles de mise à la terre en raison de la réactance inductive parasite des barres de connexion (L \approx 1\,\mu H/m ).

Par conséquent, on observe le phénomène d'élévation transitoire du potentiel de l'enveloppe (TEV, Transient Enclosure Voltage), où la coque extérieure en aluminium de la GIS subit des oscillations de potentiel atteignant plusieurs dizaines de kilovolts par rapport aux structures métalliques de support et aux sols environnants, générant des micro-arcs secondaires qui détruisent les câbles de commande du système de contrôle et de protection numérique (SAS) et provoquent des décharges par contact pour les inspecteurs de site.

Méthodologie de Stratification du Sol et d'Optimisation de la Grille

Inversion des Courbes de Sondage Électrique Vertical (SEV)

La précision des calculs de VmV_m, VsV_s et RgR_g est strictement liée à la formulation du modèle géoélectrique du sous-sol. La prospection par la méthode tétra-électrodique de Wenner consiste à injecter un courant continu ou à très basse fréquence par l'intermédiaire d'électrodes extérieures (C1,C2C_1, C_2) et à mesurer la chute de potentiel au niveau des électrodes intérieures (P1,P2P_1, P_2), disposées symétriquement avec un espacement aa :

ρa(a)=2πaΔVI\rho_a(a) = 2\pi a \frac{\Delta V}{I}

Étant donné que le sol est rarement homogène, la résistivité apparente \rho_a varie en fonction de la profondeur de pénétration des lignes de courant, laquelle est proportionnelle à l'espacement aa. Pour modéliser une structure bi-couche (couche supérieure de résistivité \rho_1 et d'épaisseur hh, reposant sur un demi-espace profond de résistivité \rho_2), on résout analytiquement l'équation intégrale de Stefanescu au moyen de la transformée de Hankel d'ordre zéro :

ρa(a)=ρ1[1+4n=1(Kn1+(2nha)2Kn4+(2nha)2)]\rho_a(a) = \rho_1 \left[ 1 + 4 \sum_{n=1}^{\infty} \left( \frac{K^n}{\sqrt{1 + \left(\frac{2 n h}{a}\right)^2}} - \frac{K^n}{\sqrt{4 + \left(\frac{2 n h}{a}\right)^2}} \right) \right]

Où le facteur de réflexion géoélectrique est :

K=ρ2ρ1ρ2+ρ1K = \frac{\rho_2 - \rho_1}{\rho_2 + \rho_1}

L'algorithme d'optimisation numérique (moindres carrés amortis non linéaires de type Levenberg-Marquardt) cherche à minimiser la fonction objectif d'erreur quadratique moyenne RMS :

\varepsilonRMS = \sqrt{\frac{1}{N} \sum_{i=1}^{N} \left( \frac{\rho_{a,meas}(a_i) - \rho_{a,calc}(a_i)}{\rho_{a,meas}(a_i)} \right)^2} \times 100\%

Topologies Avancées de Atténuation des Gradients de Surface

Lorsque la résistivité superficielle est défavorable (\rho_1 \gg \rho_2) ou que le courant de défaut IGI_G est massif, la simple augmentation de conducteurs uniformes dans le réseau présente des rendements décroissants. Les techniques avancées d'optimisation électromagnétique comprennent :

  • Grille à Densification Graduelle (Compression Périmétrique) : Réduction exponentielle de l'espacement entre conducteurs à mesure que l'on s'approche de la limite extérieure de la sous-station. Cela permet d'uniformiser la densité de courant expulsée par mètre linéaire, en aplatissant les pics périphériques de potentiel qui déclenchent la tension de contact extérieure et de pas.
  • Électrodes Profondes Verticales en Périmétrie et aux Sommets : L'insertion de piquets profonds (15 à 30 mètres) forés dans les angles de la grille permet de drainer préférentiellement le courant vers les strates inférieures les plus conductrices (\rho_2), en faisant chuter le GPRGPR global et en réduisant substantiellement le gradient de surface sur les bords critiques.
  • Couche Délictrique à Haute Performance : Application de couches de gravillon lavé de granit (\rho_s \ge 3000\,\Omega\cdot m ) ou d'un tapis d'asphalte (\rho_s \ge 10000\,\Omega\cdot m ) d'épaisseur hs=0,10aˋ0,15mh_s = 0,10 à 0,15 m. Il convient de vérifier strictement que le facteur de réduction CsC_s conserve son efficacité tout au long de la durée de vie de l'installation, en atténuant la contamination par les sédiments fins et la végétation au moyen de géotextiles imperméables.

Mise en Œuvre et Modélisation Informatique dans Vexten Suite

Intégration des Courts-Circuits et Évaluation Paramétrique

Dans le flux de travail de conception en ingénierie au sein de la plateforme intégrée Vexten Suite, le dimensionnement des mises à la terre est synchronisé de manière native avec les modules de calcul de court-circuit triphasé et asymétrique selon les normes IEC 60909 et IEEE 141 (Red Book). La détermination exacte de la composante de courant de séquence zéro (3I03I_0) et de l'impédance équivalente de Thévenin au nœud de la sous-station alimente directement le moteur d'optimisation des terres.

La formulation matricielle suivante est implémentée dans Vexten Suite pour la résolution du couplage inductif et galvanique de multiples lignes aériennes équipées d'un câble de garde (OPGW/acier) raccordées à la sous-station, permettant de déterminer avec une rigueur absolue le facteur de division réel SfS_f :

[Ig1Ig2Ign]=[Z11Z12Z1nZ21Z22Z2nZn1Zn2Znn]1[VGPREind,1VGPREind,2VGPREind,n]\begin{bmatrix} \mathbf{I}_{g1} \\ \mathbf{I}_{g2} \\ \vdots \\ \mathbf{I}_{gn} \end{bmatrix} = \begin{bmatrix} Z11 & Z12 & \cdots & Z1n \\ Z21 & Z22 & \cdots & Z2n \\ \vdots & \vdots & \ddots & \vdots \\ Zn1 & Zn2 & \cdots & Znn \end{bmatrix}^{-1} \begin{bmatrix} \mathbf{V}_{GPR} - \mathbf{E}_{ind,1} \\ \mathbf{V}_{GPR} - \mathbf{E}_{ind,2} \\ \vdots \\ \mathbf{V}_{GPR} - \mathbf{E}_{ind,n} \end{bmatrix}

Cartographie 3D des Gradients de Surface et Vérification Normative

Le moteur électromagnétique de Vexten Suite subdivise la grille de terre et les électrodes verticales en segments cylindriques discrets, en appliquant la méthode des moments (MoM) pour calculer la densité de charge longitudinale qjq_j et le courant de fuite radial sur chaque segment. La distribution de potentiel en n'importe quelle coordonnée de surface (x, y, 0) est évaluée analytiquement par intégration numérique :

Phi(x,y,0)=14πσj=1MLjIj()(xx)2+(yy)2+(z)2dPhi(x, y, 0) = \frac{1}{4\pi \sigma} \sum_{j=1}^{M} \int_{L_j} \frac{I_j(\ell')}{\sqrt{(x - x')^2 + (y - y')^2 + (z')^2}} \, d\ell'

Une fois la carte tridimensionnelle des potentiels calculée, Vexten Suite traite automatiquement la topologie de la centrale en superposant les appareils de manœuvre, les transformateurs de puissance, les clôtures périmétriques et les chemins de circulation. Le système génère des vecteurs spatiaux de 1 mètre dans toutes les directions radiales à partir des masses conductrices accessibles afin de confronter point par point :

\max\left( \Phistructure - \Phi(xpied, ypied, 0) \right) \le V_{touch,adm}(t_s)
max(Φ(x+1,y,0)Φ(x,y,0))Vstep,adm(ts)\max\left( \Phi(x+1, y, 0) - \Phi(x, y, 0) \right) \le V_{step,adm}(t_s)

Toute transgression des marges normatives des normes IEEE Std 80 ou IEC 61936-1 est mise en évidence dans l'environnement graphique au moyen de gradients de couleur à haute résolution (cartes thermiques ou heatmaps), permettant à l'ingénieur d'appliquer de manière interactive des piquets auxiliaires, des égalisateurs de gradient périmétriques ou des modifications de l'épaisseur de la couche de gravier jusqu'à certifier la conformité totale de l'installation face aux exigences maximales de sécurité humaine et de robustesse opérationnelle.