Mesure de la résistivité apparente du sol par la méthode de Wenner à 4 piquets et modélisation bicouche selon IEEE 81 et IEEE 80

Fondements Électromagnétiques de la Conduction dans les Milieux Géoélectriques et Méthode de Wenner

La caractérisation géoélectrique du sous-sol constitue le piler fondamental et incontournable pour la conception ingénierique des systèmes de mise à la terre (SMT) dans les postes électriques de moyenne, haute et très haute tension. À la différence des conducteurs métalliques purs où le transport de charge est exclusivement régi par les électrons libres sous la loi d'Ohm microscopique, le sol est un milieu hétérogène, anisotrope et multiphasique. La conduction électrique dans le terrain s'effectue de manière prédominante à travers un mécanisme électrolytique au sein du fluide poreux interstitiel, gouverné par la salinité, l'humidité volumétrique, la température, la compacité et la composition minéralogique de la matrice solide.

À l'échelle macroscopique, le comportement du champ électrostatique quasi-stationnaire dans le demi-espace conducteur infini inférieur ( $z \ge 0$ ) est décrit par les équations de Maxwell pour les milieux conducteurs en régime stationnaire, où la densité de courant $\mathbf{J}$ et le potentiel scalaire électrique $\Phi$ satisfont :

\nabla \cdot \mathbf{J} = 0 \quad \implies \quad \nabla \cdot (\sigma(x,y,z) \nabla \Phi) = 0

Pour un milieu homogène et isotrope à conductivité scalaire constante $\sigma = 1/\rho$ , l'équation différentielle gouvernante se réduit à l'équation de Laplace standard $\nabla^2 \Phi = 0$ , soumise aux conditions aux limites de Neumann à l'interface sol-air ( $z = 0$ ) :

\left. \frac{\partial \Phi}{\partial z} \right|_{z=0} = 0

Dérivation du Potentiel pour Sources Ponctuelles et Dispositif Tétrapolaire Symétrique

En considérant une source de courant ponctuelle $I$ injectée à la surface d'un demi-espace homogène ( $z = 0$ ), la symétrie sphérique dans le demi-espace impose que les lignes de courant divergent radialement dans un angle solide de $2\pi$ stéradians. Le potentiel électrique $\Phi(r)$ à une distance radiale $r$ de la source s'obtient en intégrant la loi d'Ohm en coordonnées sphériques :

J(r) = \frac{I}{2\pi r^2} = -\sigma \frac{ d \Phi}{ d r} \implies \Phi(r) = \frac{\rho I}{2\pi r}

Dans la méthode tétrapolaire proposée par Frank Wenner (1915), quatre électrodes colinéaires équidistantes sont disposées avec un écartement inter-électrodes $a$ . Les deux électrodes extérieures ( $C_1$ et $C_2$ ) agissent comme injecteurs de courant ( $+I$ et $-I$ ), tandis que les deux électrodes intérieures ( $P_1$ et $P_2$ ) agissent comme sondes de mesure du potentiel différentiel $\Delta V = \Phi(P_1) - \Phi(P_2)$ .

En appliquant le principe de superposition linéaire pour les sources de courant situées en $C_1$ ( $x = 0$ ) et $C_2$ ( $x = 3a$ ), les potentiels aux électrodes de potentiel situées en $P_1$ ( $x = a$ ) et $P_2$ ( $x = 2a$ ) deviennent :

\Phi(P_1) = \frac{\rho I}{2\pi} \left[ \frac{1}{a} - \frac{1}{2a} \right] = \frac{\rho I}{4\pi a}

\Phi(P_2) = \frac{\rho I}{2\pi} \left[ \frac{1}{2a} - \frac{1}{a} \right] = -\frac{\rho I}{4\pi a}

La différence de potentiel nette mesurée par le voltmètre haute impédance est donnée par :

\Delta V = \Phi(P_1) - \Phi(P_2) = \frac{\rho I}{4\pi a} - \left( -\frac{\rho I}{4\pi a} \right) = \frac{\rho I}{2\pi a}

En isolant le paramètre de résistivité, on définit la résistivité apparente $\rho_a$ pour un écartement inter-électrodes $a$ :

\rho_a = 2\pi a \frac{\Delta V}{I} = 2\pi a R_m

où $R_m = \Delta V / I$ est la résistance de transfert mutuel apparente lue par le telluromètre ou tellurohmètre numérique basse fréquence.

Correction pour Profondeur Finie d'Enfouissement des Piquets

Dans les conditions réelles de terrain, les piquets d'essai sont enfoncés jusqu'à une profondeur finie $b$ . Si le ratio $b/a$ n'est pas négligeable ( $b > 0.1a$ ), l'approximation ponctuelle de surface introduit une erreur systématique significative. En considérant des électrodes cylindriques verticales de longueur $b$ , la formulation mathématique rigoureuse corrigée de l'effet de pénétration des piquets s'exprime par :

\rho_a = \frac{4\pi a R_m}{1 + \frac{2a}{\sqrt{a^2 + 4b^2}} - \frac{a}{\sqrt{a^2 + b^2}} - \frac{2a}{\sqrt{4a^2 + 4b^2}} + \frac{a}{\sqrt{4a^2 + b^2}}}

Lorsque $b \to 0$ , le dénominateur de la fraction meut exactement vers la valeur $2$ , retrouvant l'équation classique de Wenner $\rho_a = 2\pi a R_m$ . Dans la pratique selon la norme IEEE Std 81, il est exigé de maintenir $b \le 0.05a$ afin de négliger la correction géométrique sans dépasser une erreur analytique de $1\%$ .

Sources d'Incertitude Métrologique, Perturbations et Couplages sur le Terrain

L'obtention de profils de résistivité apparente à haute fidélité est sujette à diverses perturbations électromagnétiques et non-idéalités de l'environnement qui dégradent sévèrement la fonction de transfert du sous-sol. Les facteurs les plus critiques sont résumés ci-dessous :

Source de Perturbation	Mécanisme Physique / Électromagnétique	Effet sur $\rho_a$	Stratégie d'Atténuation Normative
Couplage Inductif entre Câbles	Induction mutuelle $M \cdot \frac{ d i}{ d t}$ entre la boucle de courant $C_1-C_2$ et la boucle de potentiel $P_1-P_2$ .	Surestimation artificielle de $\Delta V$ , avec déphasage apparent aux fréquences supérieures à 100 Hz.	Séparation physique transversale des câbles de courant et de tension ; utilisation de câbles coaxiaux ou bifilaires blindés ; injection à des fréquences sub-harmoniques (ex. 55 Hz, 94 Hz, 105 Hz, 128 Hz).
Résistances de Contact Élevées ( $R_c$ )	Impédance d'interface électrode-sol excessive due à des graviers secs, des roches ou une faible humidité superficielle.	Saturation de l'étage d'entrée du voltmètre, atténuation drastique de $I$ , bruit thermique de Johnson-Nyquist et lectures erratiques.	Humidification des piquets avec une solution saline ou de la bentonite ; insertion de multiples piquets en parallèle aux bornes de courant.
Courants Parasites et Telluriques	Retours de neutre déséquilibrés, traction électrique en CC/CA, gradients galvaniques et activité géomagnétique.	Distorsion du signal différentiel sur $P_1-P_2$ , dérive du zéro en CC et fluctuations chaotiques à basse fréquence.	Filtrage synchrone par modulation de phase (détection synchrone / Lock-in) ; injection de courant rectangulaire avec inversion périodique de polarité.
Structures Métalliques Enfouies	Conduites de gaz/eau, grilles de terre existantes et blindages de câbles agissant comme chemins de faible impédance.	Sous-estimation critique de la résistivité aux grands écartements ( $a > 10 m$ ), distordant la stratification.	Tracé de profils orthogonaux et diagonaux ; éloignement des infrastructures métalliques enfouies à une distance supérieure à $3a_{\max}$ .

Théorie Mathématique de la Stratification Géoélectrique Bi-couche

Dans la grande majorité des formations géologiques présentant un intérêt pour l'ingénierie électrique, le sous-sol ne peut pas être modélisé avec précision par un demi-espace homogène. La sédimentation, le degré de compactage et la teneur en eau phrreatique engendrent des gradients verticaux marqués. Le modèle canonique fondamental est le modèle stratifié bi-couche, composé d'une première couche superficielle d'épaisseur finie $h$ et de résistivité intrinsèque $\rho_1$ , superposée à un demi-espace infini inférieur de résistivité intrinsèque $\rho_2$ .

\rho(z) = \begin{cases} \rho_1, & 0 \le z \le h \\ \rho_2, & z > h \end{cases}

Résolution de l'Équation de Laplace avec Conditions aux Limites

Pour une source ponctuelle injectant un courant $I$ à l'origine $(0,0,0)$ à la surface ( $z=0$ ), le champ scalaire de potentiel en coordonnées cylindriques $(\phi(r,z))$ satisfait l'équation de Laplace dans chaque milieu homogène :

\frac{\partial^2 \Phi_i}{\partial r^2} + \frac{1}{r}\frac{\partial \Phi_i}{\partial r} + \frac{\partial^2 \Phi_i}{\partial z^2} = 0 \quad (i = 1, 2)

En appliquant la méthode de séparation des variables au moyen de la transformée intégrale de Fourier-Bessel (d'ordre zéro), la solution générale pour le potentiel dans la première couche ( $0 \le z \le h$ ) et dans la seconde couche ( $z \ge h$ ) se formule comme :

\Phi_1(r, z) = \frac{\rho_1 I}{2\pi} \int0^{\infty} \left[ e^{-\lambda z} + A(\lambda) e^{-\lambda z} + B(\lambda) e^{\lambda z} \right] J_0(\lambda r) \, d \lambda

\Phi_2(r, z) = \frac{\rho_2 I}{2\pi} \int0^{\infty} C(\lambda) e^{-\lambda z} J_0(\lambda r) \, d \lambda

où $J_0(\lambda r)$ est la fonction de Bessel de première espèce et d'ordre zéro, et les coefficients $A(\lambda)$ , $B(\lambda)$ et $C(\lambda)$ sont déterminés en imposant les conditions aux limites de continuité du potentiel et de la composante normale de la densité de courant :

Condition à la surface air-sol ( $z = 0$ ) : $\left. \frac{\partial \Phi_1}{\partial z} \right|_{z=0} = 0 \implies A(\lambda) = B(\lambda)$
Continuité du potentiel à l'interface ( $z = h$ ) : $\Phi_1(r, h) = \Phi_2(r, h)$
Continuité de la densité de courant normale à l'interface ( $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}$
Condition de régularité à l'infini : $\lim_{z \to \infty} \Phi_2(r, z) = 0$

En définissant le coefficient de réflexion de résistivité $k$ comme :

k = \frac{\rho_2 - \rho_1}{\rho_2 + \rho_1} \quad avec -1 \le k \le 1

En résolvant le système algébrique linéaire pour $A(\lambda)$ , on obtient :

A(\lambda) = B(\lambda) = \frac{k e^{-2\lambda h}}{1 - k e^{-2\lambda h}}

En développant le dénominateur en série géométrique infinie $(1 - k e^{-2\lambda h})^{-1} = \sum_{n=1}^{\infty} k^n e^{-2n\lambda h}$ , la fonction de potentiel sur la surface libre ( $z = 0$ ) adopte la forme :

\Phi_1(r, 0) = \frac{\rho_1 I}{2\pi} \left[ \frac{1}{r} + 2 \sum_{n=1}^{\infty} k^n \int0^{\infty} e^{-2n\lambda h} J_0(\lambda r) \, d \lambda \right]

En utilisant l'identité intégrale de Lipschitz-Weber :

\int0^{\infty} e^{-\alpha \lambda} J_0(\lambda r) \, d \lambda = \frac{1}{\sqrt{r^2 + \alpha^2}}

On parvient à l'expression analytique fermée du potentiel électrostatique superficiel issue de la méthode des images infinies :

\Phi(r) = \frac{\rho_1 I}{2\pi} \left[ \frac{1}{r} + 2 \sum_{n=1}^{\infty} \frac{k^n}{\sqrt{r^2 + (2nh)^2}} \right]

Équation Analytique de la Résistivité Apparente de Wenner

En appliquant la formule générale du potentiel de surface au dispositif à quatre points de Wenner, le potentiel mesuré entre $P_1$ et $P_2$ face aux courants $+I$ en $C_1$ et $-I$ en $C_2$ donne naissance à la formulation maîtresse de résistivité apparente pour un sol bi-couche :

\rho_a(a) = \rho_1 \left[ 1 + 4 \sum_{n=1}^{\infty} \frac{k^n}{\sqrt{1 + \left(\frac{2nh}{a}\right)^2}} - 2 \sum_{n=1}^{\infty} \frac{k^n}{\sqrt{4 + \left(\frac{2nh}{a}\right)^2}} \right]

Cette série alternée converge absolument pour tout $|k| < 1$ et permet de calculer théoriquement le profil de $\rho_a(a)$ en fonction de l'écartement inter-électrodes $a$ . Elle présente deux asymptotes physiques fondamentales :

\lim_{a \to 0} \rho_a(a) = \rho_1 \qquad \lim_{a \to \infty} \rho_a(a) = \rho_2

\begin{cases} k > 0 \implies \rho_2 > \rho_1 & ( Profil Ascendant : couche inférieure plus résistive ) \\ k < 0 \implies \rho_2 < \rho_1 & ( Profil Descendant : couche inférieure plus conductrice ) \\ k = 0 \implies \rho_2 = \rho_1 & ( Sol Homogène : \rho_a(a) = \rho_1 = constante ) \end{cases}

Algorithmes d'Inversion Géoélectrique et Ajustement de Courbes Non Linéaires

Dans la pratique de l'ingénierie forensique et de conception, le problème consiste à déterminer le vecteur de paramètres inconnus du terrain $\mathbf{p} = [\rho_1, \rho_2, h]^T$ à partir d'un ensemble discret de $M$ mesures expérimentales $\{(\ln a_i, \ln \rho_{a,i}^{mes})\}_{i=1}^M$ . Ce problème inverse est intrinsèquement non linéaire et mal conditionné.

Formulation du Problème Inverse par Moindres Carrés Non Linéaires

On définit la fonction objectif d'erreur quadratique relative normalisée $\chi^2(\mathbf{p})$ à l'échelle logarithmique afin d'équilibrer les sensibilités sur des ordres de grandeur dispars :

\chi^2(\mathbf{p}) = \sum_{i=1}^{M} \left[ \frac{\rho_{a,i}^{mes} - \rho_a(a_i, \mathbf{p})}{\rho_{a,i}^{mes}} \right]^2 = \sum_{i=1}^{M} r_i(\mathbf{p})^2 = \|\mathbf{r}(\mathbf{p})\|_2^2

L'optimisation est réalisée au moyen de l'Algorithme de Levenberg-Marquardt (Gauss-Newton amorti), lequel interpole de manière adaptative entre la méthode de la plus grande pente (descente de gradient) et la méthode de Gauss-Newton linéarisée :

\left[ \mathbf{J}^T(\mathbf{p}_k) \mathbf{J}(\mathbf{p}_k) + \mu_k \operatorname{diag}(\mathbf{J}^T(\mathbf{p}_k) \mathbf{J}(\mathbf{p}_k)) \right] \Delta \mathbf{p}_k = -\mathbf{J}^T(\mathbf{p}_k) \mathbf{r}(\mathbf{p}_k)

où $\mathbf{J} \in \mathbb{R}^{M \times 3}$ est la matrice jacobienne des sensibilités partielles :

Jij = \frac{\partial r_i(\mathbf{p})}{\partial p_j} \implies \left[ \frac{\partial \rho_a}{\partial \rho_1}, \frac{\partial \rho_a}{\partial \rho_2}, \frac{\partial \rho_a}{\partial h} \right]

et $\mu_k \ge 0$ est le facteur d'amortissement (damping factor) ajusté à chaque itération en fonction de la métrique de gain de réduction résiduelle.

Matrice de Sensibilité Analytique (Dérivées Partielles de Fréchet)

À la différence du calcul numérique par différences finies qui introduit des erreurs d'arrondi, le calcul analytique du gradient accélère la convergence asymptotique :

\frac{\partial \rho_a}{\partial \rho_1} = \frac{\rho_a}{\rho_1} + \rho_1 \frac{\partial k}{\partial \rho_1} \sum_{n=1}^{\infty} n k^{n-1} \left[ \frac{4}{\sqrt{1 + \left(\frac{2nh}{a}\right)^2}} - \frac{2}{\sqrt{4 + \left(\frac{2nh}{a}\right)^2}} \right]

où $\frac{\partial k}{\partial \rho_1} = \frac{-2\rho_2}{(\rho_1 + \rho_2)^2}$ et $\frac{\partial k}{\partial \rho_2} = \frac{2\rho_1}{(\rho_1 + \rho_2)^2}$ .

\frac{\partial \rho_a}{\partial h} = -4\rho_1 \sum_{n=1}^{\infty} k^n \left[ \frac{4n^2 h / a^2}{\left(1 + \left(\frac{2nh}{a}\right)^2\right)^{3/2}} - \frac{2n^2 h / a^2}{\left(4 + \left(\frac{2nh}{a}\right)^2\right)^{3/2}} \right]

Le critère de convergence strict exige que l'erreur quadratique moyenne relative résiduelle (RMS) satisfasse :

RMS _{ erreur } = \sqrt{\frac{1}{M} \sum_{i=1}^{M} \left( \frac{\rho_{a,i}^{mes} - \rho_a(a_i, \mathbf{p})}{\rho_{a,i}^{mes}} \right)^2} \times 100\% \le 5.0\%

Impact de la Stratification Bi-couche sur la Conception des Grilles de Mise à la Terre (IEEE Std 80 / IEC 60479)

Omettre la stratification du terrain et adopter un modèle de résistivité constante équivalente moyennée arithmétiquement engendre des erreurs catastrophiques dans le dimensionnement des conducteurs en cuivre, la longueur des piquets et l'évaluation des tensions de pas et de contact dans les postes électriques.

Résistance de Mise à la Terre ( $R_g$ ) sous Modèle Bi-couche

Selon les formulations étendues de Sverak et Schwarz dans les milieux bi-couches, la résistance totale d'une grille de mise à la terre carrée ou rectangulaire dotée d'électrodes verticales dépend de manière critique de la profondeur de la première couche $h$ et du contraste de résistivités $\rho_2 / \rho_1$ :

R_g(\rho_1, \rho_2, h) \approx \frac{\rho_{eq}}{\pi L} \left[ \ln\left(\frac{2L}{a'}\right) + k_1 \frac{L}{\sqrt{A}} - k_2 \right]

où $\rho_{eq}$ est la résistivité apparente équivalente vue par la géométrie de la grille :

\rho_{eq} = \rho_1 \rho_2 \frac{Ltotal}{\rho_1 (Ltotal - L_v) + \rho_2 L_v \cdot \psi(h)}

où $Ltotal = L_c + L_v$ est la longueur combinée des conducteurs horizontaux ( $L_c$ ) et des piquets verticaux ( $L_v$ ), et $\psi(h)$ est la fonction de couplage de profondeur.

Potentiels de Pas et de Contact : Comparaison Phénoménologique de Scénarios

Paramètre de Conception / Sécurité	Scénario A : Terrain Favorable ( $\rho_1 > \rho_2$ , $k < 0$ )	Scénario B : Terrain Critique ( $\rho_1 < \rho_2$ , $k > 0$ )	Implication Critique de Sécurité
Distribution du Courant de Défaut	Le courant injecté a tendance à s'écouler verticalement vers la couche inférieure profonde à faible résistivité.	Le courant est repoussé par la couche inférieure très résistive (roche), s'acheminant horizontalement dans la couche superficielle.	Dans le Scénario B, la densité de courant superficielle $J_s$ augmente de manière exponentielle.
Potentiel de Maille ( $V_m$ ) et de Contact ( $Etouch$ )	Gradients de potentiel superficiel plats et atténués. $Etouch$ reste dans les limites admissibles avec des conducteurs standards.	Élévation sévère du Gradient de Potentiel de Surface (GPR) et crêtes de tension extrêmement prononcées au-dessus de la grille.	Dans le Scénario B, la limite de fibrillation ventriculaire selon IEEE 80 / IEC 60479 est fréquemment dépassée.
Efficacité des Piquets de Mise à la Terre	Hautement efficaces s'ils pénètrent la couche profonde ( $l_v > h$ ), réduisant drastiquement $R_g$ .	Peu efficaces ; les piquets ne dissipent pratiquement aucun courant supplémentaire dans la couche rocheuse profonde à haute résistivité.	Dans le Scénario B, il est nécessaire de densifier le maillage horizontal et d'étendre la surface périmétrique.
Facteur de Réduction de la Couche de Gravier ( $C_s$ )	L'épaisseur superficielle de gravier ( $h_s$ ) offre une réduction standard ( $C_s \approx 0.7 - 0.85$ ).	L'interaction d'une faible résistivité $\rho_1$ avec la roche sous-jacente et un $\rho_2$ infini altère la réflexion du lit de gravier.	Nécessite de recalculer $C_s$ par des séries infinies de double réflexion (gravier-couche 1-couche 2).

Calcul Rigoureux du Facteur de Réduction de la Couche de Gravier ( $C_s$ )

Le facteur de correction $C_s$ pour la couche de matériau à haute résistivité superficielle (gravier/concassé $\rho_s \approx 3000\,\Omega\cdot m$ ) d'épaisseur $h_s$ installée sur un sol bi-couche ne peut pas être calculé avec la formule empirique classique de Sunde pour sol homogène. La formulation généralisée par expansion en série de coefficients combinés s'exprime sous la forme :

C_s = 1 + \frac{16}{\pi} \sum_{m=1}^{\infty} \frac{k_s^m}{2m-1} \arctan\left(\frac{2m-1}{2h_s/r_d}\right)

où $k_s = (\rho_1 - \rho_s)/(\rho_1 + \rho_s)$ et $r_d = 0.08 m$ représente le rayon équivalent du pied humain modélisé sous la forme d'un disque circulaire métallique de rayon normalisé.

Analyse Forensique des Défaillances Associées à une Caractérisation Géoélectrique Défaillante

Les erreurs méthodologiques commises lors de la campagne de prospection de résistivité ou la simplification à un sol homogène lorsqu'il existe une couche sous-jacente fortement réfléchissante ( $k \to +1$ ) ont historiquement conduit à des défaillances catastrophiques dans les réseaux électriques de puissance. On présente ci-après les mécanismes de défaillance et les incidents forensiques documentés :

Claquage Diélectrique dans les Transformateurs de Mesure et Services Auxiliaires par Élévation Critique du GPR

Dans un poste électrique de 230 kV construit sur une formation de grès superficiel ( $\rho_1 = 180\,\Omega\cdot m$ , $h = 1.8 m$ ) reposant sur un socle granitique massif ( $\rho_2 = 4200\,\Omega\cdot m$ , $k = +0.918$ ), la conception originale avait supposé un sol homogène moyen de $350\,\Omega\cdot m$ . Lors d'un court-circuit monophasé à la terre avec un courant symétrique $I_f = 22 kA$ , la résistance réelle de la grille s'est élevée à $2.85\,\Omega$ (contre $0.68\,\Omega$ projetés sous l'hypothèse du sol homogène).

GPR = I_g \cdot R_g = (0.75 \times 22000 A ) \times 2.85\,\Omega = 47.025 kV

Le courant dense qui n'a pas pu pénétrer le granite profond est revenu latéralement par les tranchées de câbles de contrôle, élevant le potentiel local de terre des armoires extérieures au-delà de la rigidité diélectrique des enroulements secondaires des transformateurs de courant (TC) et des transformateurs de tension (TT). La différence de potentiel en mode commun a perforé l'isolement galvanique de 10 kV des relais de protection numériques (IED) dans la salle de commande, détruisant les modules d'entrées/sorties analogiques et inhibant le déclenchement de secours temporisé (ANSI 50/51N). Cela a prolongé l'élimination du défaut de 80 ms à 1,2 seconde, provoquant l'explosion destructive du transformateur de puissance principal.

Dégradation Thermique et Fusion des Écrans de Câbles Souterrains de Haute Tension

Sur des liaisons de câbles isolés au XLPE de 115 kV comportant des écrans métalliques mis à la terre aux deux extrémités (solid bonding), la fausse hypothèse d'une faible résistivité dans les couches profondes distord l'impédance mutuelle de boucle terre-écran $\underline{Z}_{m(c-g)}$ selon les formules de Carson-Clem :

\underline{Z}_{m(c-g)} = \frac{\mu_0 \omega}{8} + j \frac{\mu_0 \omega}{2\pi} \ln\left(\frac{D_e}{Dij}\right) \quad avec D_e = 658.87 \sqrt{\frac{\rho_{eq}}{f}}

En présence d'une couche profonde de haute résistivité ( $\rho_2 \gg \rho_1$ ), la profondeur de retour équivalente $D_e$ s'est accrue de plus de 300 %. Cela a augmenté drastiquement la réactance inductive du retour par la terre, forçant le courant de déséquilibre et de défaut à circuler quasi intégralement à travers les écrans métalliques des câbles de puissance. La densité de courant thermique $Jscreen$ a dépassé la limite admissible du cuivre selon IEC 60287 et IEC 60949 :

Iadm^2 t > kmat^2 Sscreen^2 \ln\left( \frac{\theta_f + \beta}{\theta_i + \beta} \right)

provoquant le délaminage thermique de la gaine extérieure semiconductrice, le ramollissement du polyéthylène réticulé (XLPE) par une surtempérature supérieure à 250 °C et le claquage diélectrique phase-terre consécutif par perforation d'arborescences d'eau accélérées de façon électrothermique.

Flux Méthodologique de Simulation et Conception Avancée dans Vexten Suite

L'écosystème d'ingénierie intégrée Vexten Suite intègre un moteur déterministe et multiphysique d'inversion géoélectrique couplé à des modules de court-circuit (IEC 60909 / IEEE 141), de dimensionnement de conducteurs et câbles (IEC 60287 / NEC) et d'optimisation de réseaux de terre (IEEE 80 / IEC 60479). On détaille ci-après l'intégration algorithmique et le protocole méthodologique appliqué aux projets d'infrastructures critiques :

\begin{matrix} \boxed{\begin{array}{c} \textbf{1. Acquisition de Données de Terrain} \\ Mesure Wenner (a_i, R_{m,i}) \\ IEEE Std 81 \end{array}} & \longrightarrow & \boxed{\begin{array}{c} \textbf{2. Inversion Géoélectrique} \\ Levenberg-Marquardt dans Vexten \\ \mathbf{p} = [\rho_1, \rho_2, h]^T \end{array}} \\ \downarrow & & \downarrow \\ \boxed{\begin{array}{c} \textbf{4. Calcul des Tensions Admissibles} \\ Limites Etouch, Estep \\ IEEE 80 / IEC 60479 (50kg/70kg) \end{array}} & \longleftarrow & \boxed{\begin{array}{c} \textbf{3. Simulation de Défaut de Court-Circuit} \\ Injection I_k'' / I_g (IEC 60909) \\ Division de Courant S_f \end{array}} \\ \downarrow & & \downarrow \\ \boxed{\begin{array}{c} \textbf{5. Modélisation Électromagnétique (BEM/FEM)} \\ Grille SMT dans un Milieu Bi-couche \\ Calcul de R_g, GPR , Vmesh \end{array}} & \longrightarrow & \boxed{\begin{array}{c} \textbf{6. Optimisation et Atténuation} \\ Ajustement du Maillage, Piquets Profonds \\ Épaisseur de Gravier h_s et Écrans \end{array}} \end{matrix}

Exemple Pratique de Simulation et Inversion Numérique dans Vexten Grounding Engine

Une campagne de prospection géoélectrique tétrapolaire selon la méthode de Wenner est réalisée sur le site d'un nouveau poste électrique industriel 138/13.8 kV. L'acquisition des données a fourni les valeurs de terrain suivantes :

Écartement $a$ [m]	Résistance Mesurée $R_m$ [ $\Omega$ ]	Résistivité Apparente Expérimentale $\rho_a^{mes}$ [ $\Omega\cdot m$ ]	Résistivité Apparente Calculée $\rho_a^{calc}$ [ $\Omega\cdot m$ ]	Erreur Résiduelle Relative [\%]
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

En traitant la matrice de données dans le module d'optimisation inverse non linéaire de Vexten Grounding, l'algorithme de Levenberg-Marquardt converge en 7 itérations avec une erreur globale RMS de 1.32%, déterminant le vecteur d'état stratifié suivant :

\mathbf{p}^* = \begin{bmatrix} \rho_1 \\ \rho_2 \\ h \end{bmatrix} = \begin{bmatrix} 435.20\,\Omega\cdot m \\ 82.50\,\Omega\cdot m \\ 4.35 m \end{bmatrix} \implies k = \frac{82.50 - 435.20}{82.50 + 435.20} = -0.6812

Analyse d'Ingénierie et Optimisation de la Conception du SMT

Étant donné que le coefficient de réflexion $k$ est fortement négatif ( $k = -0.6812$ ), la couche profonde présente une résistivité plus de cinq fois inférieure à celle de la couche superficielle ( $\rho_2 \ll \rho_1$ ). L'optimisation numérique dans Vexten Grounding exploite ce phénomène physique selon les directives exécutives suivantes :

Installation Stratégique de Piquets Profonds : Il est prescrit d'enfoncer des piquets en acier cuivre (copperweld) d'une longueur minimale $l_v = 6.0 m$ aux quatre coins et aux nœuds périmétriques critiques. En dépassant la profondeur de l'interface ( $h = 4.35 m$ ), les électrodes verticales pénètrent de $1.65 m$ dans la couche à faible résistivité ( $\rho_2 = 82.50\,\Omega\cdot m$ ). Cela réduit la résistance totale du système $R_g$ de 42.7% par rapport à une grille purement horizontale possédant la même longueur totale de conducteur.
Allègement des Conducteurs Horizontaux : Grâce à la forte admittance de la couche inférieure, la maille entre conducteurs horizontaux peut être élargie d'un quadrillage dense de $5 m \times 5 m$ à $10 m \times 10 m$ , permettant d'économiser une quantité importante de cuivre sans dépasser les limites de tension de pas et de contact imposées par la norme IEEE Std 80.
Intégration avec le Module de Court-Circuit de Vexten Suite (IEC 60909) : L'impédance réelle du terrain stratifié est exportée automatiquement vers le moteur de court-circuit de Vexten Suite, permettant de recalculer de manière couplée les facteurs de répartition du courant dans les câbles de garde et écrans de phase ( $S_f$ ). Ceci garantit que le courant retournant effectivement par le sol ( $I_g = S_f \cdot I_k''$ ) reflète fidèlement la physique multifilaire du réseau de transport et de distribution environnant.

L'adoption de ce flux méthodologique global garantit la plus haute rigueur analytique, la sécurité diélectrique du personnel et des équipements électrothermiques, ainsi qu'une optimisation économique et opérationnelle lors de la construction des systèmes de mise à la terre pour postes électriques et centrales de production.