Bodenwiderstandsmessung nach dem Wenner-Vierpunktverfahren und Zweischicht-Bodenmodellierung gemäß IEEE 81 und IEEE 80

Elektromagnetische Grundlagen der Leitung in geoelektrischen Medien und die Wenner-Methode

Die geoelektrische Charakterisierung des Untergrunds bildet das fundamentale und unverzichtbare Fundament für die ingenieurmäßige Auslegung von Erdungsanlagen in Mittel-, Hoch- und Höchstspannungsschaltanlagen. Im Gegensatz zu reinen metallischen Leitern, in denen der Ladungstransport ausschließlich durch freie Elektronen gemäß dem mikroskopischen Ohmschen Gesetz geregelt wird, ist der Boden ein heterogenes, anisotropes und mehrphasiges Medium. Die elektrische Leitung im Erdreich erfolgt überwiegend über einen elektrolytischen Mechanismus innerhalb der interstitiellen Porenflüssigkeit, der durch den Salzgehalt, die volumetrische Feuchte, die Temperatur, die Kompaktizität und die mineralogische Zusammensetzung der festen Matrix bestimmt wird.

Auf makroskopischer Ebene wird das Verhalten des quasi-stationären elektrostatischen Feldes im unteren unendlichen leitfähigen Halbraum ( $z \ge 0$ ) durch die Maxwell-Gleichungen für leitfähige Medien im stationären Zustand beschrieben, wobei die Stromdichte $\mathbf{J}$ und das elektrische skalare Potenzial $\Phi$ folgender Beziehung genügen:

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

Für ein homogenes und isotropes Medium mit konstanter skalarer Leitfähigkeit $\sigma = 1/\rho$ reduziert sich die steuernde Differentialgleichung auf die Standard-Laplace-Gleichung $\nabla^2 \Phi = 0$ , unterworfen den Neumann-Randbedingungen an der Grenzfläche zwischen Boden und Luft ( $z = 0$ ):

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

Herleitung des Potenzials für Punktquellen und die symmetrische tetrapolare Anordnung

Betrachtet man eine Punktstromquelle $I$ , die an der Oberfläche eines homogenen Halbraums ( $z = 0$ ) eingespeist wird, erzwingt die Kugelsymmetrie im Halbraum, dass die Stromlinien radial in einem Raumwinkel von $2\pi$ Steradiant divergent verlaufen. Das elektrische Potenzial $\Phi(r)$ in einem radialen Abstand $r$ von der Quelle erhält man durch Integration des Ohmschen Gesetzes in Kugelkoordinaten:

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

Bei der von Frank Wenner (1915) vorgeschlagenen tetrapolaren Methode werden vier kollineare, äquidistante Elektroden mit einem Elektrodenabstand $a$ angeordnet. Die beiden äußeren Elektroden ( $C_1$ und $C_2$ ) dienen als Stromeinspeiser ( $+I$ und $-I$ ), während die beiden inneren Elektroden ( $P_1$ und $P_2$ ) als Messsonden für das Differenzpotenzial $\Delta V = \Phi(P_1) - \Phi(P_2)$ fungieren.

Unter Anwendung des linearen Superpositionsprinzips für die Stromquellen an den Positionen $C_1$ ( $x = 0$ ) und $C_2$ ( $x = 3a$ ) ergeben sich die Potenzialwerte an den Potenzialelektroden bei $P_1$ ( $x = a$ ) und $P_2$ ( $x = 2a$ ) wie folgt:

\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}

Die netto gemessene Potenzialdifferenz am hochohmigen Voltmeter ergibt sich zu:

\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}

Löst man nach dem Parameter des spezifischen Widerstands auf, wird der scheinbare spezifische Widerstand $\rho_a$ für einen Elektrodenabstand $a$ wie folgt definiert:

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

wobei $R_m = \Delta V / I$ den gegenseitigen scheinbaren Übertragungswiderstand darstellt, der vom digitalen Niederfrequenz-Erdungsmessgerät abgelesen wird.

Korrektur für die endliche Einbautiefe der Messpfähle

Unter realen Feldbedingungen werden die Testspieße bis zu einer endlichen Tiefe $b$ im Erdreich eingeschlagen. Wenn das Verhältnis $b/a$ nicht vernachlässigbar ist ( $b > 0,1a$ ), führt die Punktquellen-Näherung an der Oberfläche zu einem relevanten systematischen Fehler. Unter Berücksichtigung vertikaler zylindrischer Elektroden der Länge $b$ wird die mathematisch strenge Formulierung, korrigiert um den Penetrationseffekt der Spieße, wie folgt ausgedrückt:

\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}}}

Wenn $b \to 0$ , konvergiert der Nenner des Bruchs exakt gegen den Wert $2$ , wodurch die klassische Wenner-Gleichung $\rho_a = 2\pi a R_m$ wiederhergestellt wird. In der Praxis wird nach IEEE Std 81 gefordert, $b \le 0,05a$ einzuhalten, um die geometrische Korrektur ohne Überschreitung eines analytischen Fehlers von $1\%$ vernachlässigen zu können.

Quellen metrologischer Unsicherheit, Störungen und Feldkopplungen

Die Erfassung hochpräziser Profile des scheinbaren spezifischen Widerstands unterliegt diversen elektromagnetischen Störungen und nicht-idealen Umgebungsbedingungen, welche die Übertragungsfunktion des Untergrunds erheblich beeinträchtigen. Die kritischsten Faktoren sind nachfolgend zusammengefasst:

Störungsquelle	Physikalischer / Elektromagnetischer Mechanismus	Wirkung auf $\rho_a$	Normative Minderungsstrategie
Induktive Kopplung zwischen Kabeln	Gegeninduktion $M \cdot \frac{ d i}{ d t}$ zwischen den Stromschleifen $C_1-C_2$ und Potenzialschleifen $P_1-P_2$ .	Künstliche Überhöhung von $\Delta V$ mit scheinbarer Phasenverschiebung bei Frequenzen über 100 Hz.	Räumlich-transversale Trennung der Strom- und Spannungsleitung; Einsatz geschirmter Koaxial- oder Torsionskabel; Einspeisung mit subharmonischen Frequenzen (z. B. 55 Hz, 94 Hz, 105 Hz, 128 Hz).
Hohe Übergangswiderstände ( $R_c$ )	Übermäßige Impedanz an der Schnittstelle Elektrode-Boden aufgrund von trockenem Kies, Fels oder geringer Oberflächenfeuchte.	Sättigung der Eingangsstufe des Voltmeters, drastische Dämpfung von $I$ , thermisches Johnson-Nyquist-Rauschen und fehlerhafte Messwerte.	Befeuchtung der Messpfähle mit Salzlösung oder Bentonit; Parallelschaltung mehrerer Erder an den Stromklemmen.
Stör- und Tellurströme	Unsymmetrische Neutralleiter-Rückströme, Gleich-/Wechselstrom-Traktionssysteme, galvanische Gradienten und geomagnetische Aktivität.	Verzerrung des Differenzsignals an $P_1-P_2$ , Nullpunktdrift bei Gleichspannung und chaotische niederfrequente Schwankungen.	Phasensynchrone Filterung durch Homodyn- / Lock-in-Verfahren; Einspeisung von Rechteckstrom mit periodischer Umpolung.
Eingebettete Metallstrukturen	Unterirdische Gas-/Wasserleitungen, bestehende Erdungsnetze und Kabelschirme, die als niederohmige Pfade wirken.	Kritische Untersehätzung des spezifischen Widerstands bei großen Elektrodenabständen ( $a > 10\, m$ ), was die Schichtzuordnung verfälscht.	Aufnahme orthogonaler und diagonaler Messprofile; Einhaltung eines Mindestabstands von metallischen Untergrundinfrastrukturen von $> 3a_{\max}$ .

Mathematische Theorie der zweischichtigen geoelektrischen Stratifikation

In der überwiegenden Mehrheit der für die Elektrotechnik relevanten geologischen Formationen kann der Untergrund nicht präzise als homogener Halbraum modelliert werden. Sedimentation, Kompaktionsgrad und Grundwasserstand erzeugen ausgeprägte vertikale Gradienten. Das fundamentale kanonische Modell ist das zweischichtige Stratifikationsmodell, bestehend aus einer ersten Oberflächenschicht der endlichen Dicke $h$ und dem intrinsischen spezifischen Widerstand $\rho_1$ , die über einem unendlichen unteren Halbraum mit dem intrinsischen spezifischen Widerstand $\rho_2$ liegt.

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

Lösung der Laplace-Gleichung mit Randbedingungen

Für eine Punktquelle, die einen Strom $I$ am Ursprung $(0,0,0)$ an der Oberfläche ( $z=0$ ) einspeist, genügt das skalare Potenzialfeld in Zylinderkoordinaten ( $\Phi(r,z)$ ) der Laplace-Gleichung in jedem homogenen Medium:

\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)

Unter Anwendung der Methode der Trennung der Variablen mittels der Fouriertransformation nach Bessel-Funktionen (Ordnung Null) formuliert sich die allgemeine Lösung für das Potenzial in der ersten Schicht ( $0 \le z \le h$ ) und in der zweiten Schicht ( $z \ge h$ ) wie folgt:

\Phi_1(r, z) = \frac{\rho_1 I}{2\pi} \int_0^{\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} \int_0^{\infty} C(\lambda) e^{-\lambda z} J_0(\lambda r) \, d \lambda

wobei $J_0(\lambda r)$ die Bessel-Funktion erster Art und nullter Ordnung darstellt und die Koeffizienten $A(\lambda)$ , $B(\lambda)$ sowie $C(\lambda)$ durch Erzwingung der Randbedingungen bezüglich Kontinuität des Potenzials und der Normalkomponente der Stromdichte bestimmt werden:

Bedingung an der Oberfläche Luft-Boden ( $z = 0$ ): $\left. \frac{\partial \Phi_1}{\partial z} \right|_{z=0} = 0 \implies A(\lambda) = B(\lambda)$
Kontinuität des Potenzials an der Trennschicht ( $z = h$ ): $\Phi_1(r, h) = \Phi_2(r, h)$
Kontinuität der normalen Stromdichte an der Trennschicht ( $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}$
Regulativbedingung im Unendlichen: $\lim_{z \to \infty} \Phi_2(r, z) = 0$

Definiert man den Reflexionskoeffizienten des spezifischen Widerstands $k$ als:

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

Erhält man durch Lösung des linearen algebraischen Systems für $A(\lambda)$ :

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

Durch Entwicklung des Nenners über seine unendliche geometrische Reihe $(1 - k e^{-2\lambda h})^{-1} = \sum_{n=1}^{\infty} k^n e^{-2n\lambda h}$ nimmt die Potenzialfunktion an der freien Oberfläche ( $z = 0$ ) folgende Form an:

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

Unter Verwendung der Lipschitz-Weber-Integralidentität:

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

Gelangt man zum geschlossenen analytischen Ausdruck für das elektrostatische Oberflächenpotenzial infolge der Methode der unendlichen Spiegelbilder:

\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]

Analytische Gleichung des scheinbaren spezifischen Widerstands nach Wenner

Wendet man die allgemeine Oberflächenpotenzialformel auf die Vier-Punkt-Anordnung nach Wenner an, führt das zwischen $P_1$ und $P_2$ gemessene Potenzial bei Strömen von $+I$ an $C_1$ und $-I$ an $C_2$ zur Hauptformulierung des scheinbaren spezifischen Widerstands für einen zweischichtigen Boden:

\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]

Diese alternierende Reihe konvergiert absolut für alle $|k| < 1$ und erlaubt die theoretische Berechnung des Profils $\rho_a(a)$ in Abhängigkeit vom Elektrodenabstand $a$ . Sie weist zwei fundamentale physikalische Asymptoten auf:

\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 & ( Ansteigendes Profil: tiefere Schicht hochohmiger ) \\ k < 0 \implies \rho_2 < \rho_1 & ( Abfallendes Profil: tiefere Schicht niederohmiger ) \\ k = 0 \implies \rho_2 = \rho_1 & ( Homogener Boden: \rho_a(a) = \rho_1 = konstant ) \end{cases}

Geoelektrische Invertierungsalgorithmen und nichtlineare Kurvenanpassung

In der forensischen Engineering-Praxis und der Auslegungsberechnung besteht das Problem darin, den Vektor der unbekannten Bodenparameter $\mathbf{p} = [\rho_1, \rho_2, h]^T$ aus einem diskreten Satz von $M$ experimentellen Messungen $\{(\ln a_i, \ln \rho_{a,i}^{mess})\}_{i=1}^M$ zu bestimmen. Dieses inverse Problem ist von Natur aus nichtlinear und schlecht konditioniert.

Formulierung des inversen Problems mittels nichtlinearer kleinster Quadrate

Die Zielfunktion des normalisierten relativen quadratischen Fehlers $\chi^2(\mathbf{p})$ wird im logarithmischen Maßstab definiert, um die Sensitivitäten über unterschiedliche Größenordnungen hinweg auszugleichen:

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

Die Optimierung wird mit dem Levenberg-Marquardt-Algorithmus (gedämpfter Gauss-Newton) durchgeführt, der adaptiv zwischen der Methode des steilsten Abstiegs und der linearisierten Gauss-Newton-Methode interpoliert:

\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)

wobei $\mathbf{J} \in \mathbb{R}^{M \times 3}$ die Jacobi-Matrix der partiellen Sensitivitäten ist:

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]

und $\mu_k \ge 0$ der Dämpfungsfaktor ist, der in jeder Iteration auf der Grundlage der Kennzahl des gewollten Fehlerreduktionsgewinns angepasst wird.

Analytische Sensitivitätsmatrix (Partielle Fréchet-Ableitungen)

Im Gegensatz zur numerischen Berechnung über finite Differenzen, die Rundungsfehler einführt, beschleunigt die analytische Berechnung des Gradienten die asymptotische Konvergenz:

\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]

wobei $\frac{\partial k}{\partial \rho_1} = \frac{-2\rho_2}{(\rho_1 + \rho_2)^2}$ und $\frac{\partial k}{\partial \rho_2} = \frac{2\rho_1}{(\rho_1 + \rho_2)^2}$ gilt.

\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]

Das strikte Konvergenzkriterium verlangt, dass der mittlere relative quadratische Restfehler (RMS) folgende Bedingung erfüllt:

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

Einfluss der zweischichtigen Stratifikation auf das Design von Erdungsanlagen (IEEE Std 80 / IEC 60479)

Das Vernachlässigen der Bodenschichtung und das Annehmen eines homogenen Ersatzmodells mit arithmetisch gemitteltem spezifischem Widerstand führt zu katastrophalen Fehlern bei der Dimensionierung von Kupferleitern, der Länge von Tiefenerdern und der Bewertung von Schritt- und Berührungsspannungen in Umspannwerken.

Erdungswiderstand ( $R_g$ ) unter dem Zweischichtenmodell

Gemäß den erweiterten Formulierungen von Sverak und Schwarz für zweischichtige Medien hängt der Gesamtwiderstand einer quadratischen oder rechteckigen Erdungsanlage mit vertikalen Erdern kritisch von der Tiefe der ersten Schicht $h$ und dem Widerstandsverhältnis $\rho_2 / \rho_1$ ab:

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]

wobei $\rho_{eq}$ der äquivalente scheinbare spezifische Widerstand ist, den die Erdungsgeometrie sieht:

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

wobei $Lgesamt = L_c + L_v$ die Gesamtlänge der horizontalen Leiter ( $L_c$ ) und vertikalen Erder ( $L_v$ ) darstellt und $\psi(h)$ die Tiefenkopplungsfunktion ist.

Schritt- und Berührungsspannungen: Phänomenologischer Szenarienvergleich

Entwurfs- / Sicherheitsparameter	Szenario A: Günstiger Boden ( $\rho_1 > \rho_2$ , $k < 0$ )	Szenario B: Kritischer Boden ( $\rho_1 < \rho_2$ , $k > 0$ )	Kritische Sicherheitsauswirkung
Fehlerstromverteilung	Der eingespeiste Strom fließt vorwiegend vertikal in die tiefe, niederohmige Schicht ab.	Der Strom wird von der hochohmigen tiefen Schicht (Fels) reflektiert und kanalisiert sich horizontal in der Oberflächenschicht.	In Szenario B steigt die Oberflächenstromdichte $J_s$ exponentiell an.
Maschenpotenzial ( $V_m$ ) und Berührungsspannung ( $Etouch$ )	Flache und gedämpfte Oberflächenpotenzialgradienten. $Etouch$ bleibt mit Standardleitern innerhalb zulässiger Grenzen.	Starke Erhöhung des Oberflächenpotenzialanstiegs (GPR) und extrem ausgeprägte Spannungspitzen über der Erdungsmasche.	In Szenario B wird die Herzkammerflimmergrenze nach IEEE 80 / IEC 60479 häufig überschritten.
Wirksamkeit von Tiefenerdern	Hocheffektiv, wenn sie die tiefe Schicht durchdringen ( $l_v > h$ ), was $R_g$ drastisch reduziert.	Geringe Wirksamkeit; die Erder können im hochohmigen Tiefengestein kaum zusätzlichen Strom ableiten.	In Szenario B muss das horizontale Maschennetz verdichtet und die Erdungsfläche erweitert werden.
Reduktionsfaktor der Kiesschicht ( $C_s$ )	Die Schotterdeckschicht ( $h_s$ ) bietet Standarddämpfung ( $C_s \approx 0,7 - 0,85$ ).	Die Interaktion von niedrigem $\rho_1$ mit dem Oberflächengestein und unendlichem $\rho_2$ verändert die Reflexion des Kiesbetts.	Erfordert die Neuberechnung von $C_s$ mittels unendlicher Doppelreflexionsreihen (Kies - Schicht 1 - Schicht 2).

Präzise Berechnung des Reduktionsfaktors der Oberflächenschicht ( $C_s$ )

Der Korrekturfaktor $C_s$ für die auf einem zweischichtigen Boden aufgebrachte Schicht aus hochohmigen Oberflächenmaterial (Kies/Schotter $\rho_s \approx 3000\,\Omega\cdot m$ ) mit der Dicke $h_s$ kann nicht mit der klassischen empirischen Formel von Sunde für homogenen Boden berechnet werden. Die verallgemeinerte Formulierung durch Reihenentwicklung kombinierter Koeffizienten lautet:

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)

wobei $k_s = (\rho_1 - \rho_s)/(\rho_1 + \rho_s)$ und $r_d = 0,08\, m$ den äquivalenten Radius des menschlichen Fußes darstellt, der als metallische Kreisscheibe modelliert wird.

Forensische Analyse von Ausfällen durch mangelhafte geoelektrische Charakterisierung

Methodische Fehler bei der Erfassung des spezifischen Widerstands oder die unzulässige Vereinfachung auf einen homogenen Boden bei Vorhandensein einer stark reflektierenden Tiefenschicht ( $k \to +1$ ) haben historisch zu katastrophalen Ausfällen in Energiesystemen geführt. Nachfolgend werden Schadensmechanismen und dokumentierte forensische Ereignisse dargestellt:

Dielektrischer Durchschlag in Messwandlern und Eigenbedarfsanlagen durch kritischen Erdpotenzialanstieg (GPR)

In einem 230-kV-Umspannwerk, das auf einer Sandsteinformation an der Oberfläche ( $\rho_1 = 180\,\Omega\cdot m$ , $h = 1,8\, m$ ) über massivem Granitsockel ( $\rho_2 = 4200\,\Omega\cdot m$ , $k = +0,918$ ) errichtet wurde, nahm der ursprüngliche Entwurf einen homogenen Boden von durchschnittlich $350\,\Omega\cdot m$ an. Bei einem einpoligen Erdschluss mit einem symmetrischen Kurzschlussstrom von $I_f = 22\, kA$ betrug der reale Widerstand der Erdungsanlage $2,85\,\Omega$ (verglichen mit den projizierten $0,68\,\Omega$ unter der Annahme eines homogenen Bodens).

GPR = I_g \cdot R_g = (0,75 \times 22000\, A ) \times 2,85\,\Omega = 47,025\, kV

Der konzentrierte Strom, der nicht in den tiefen Granit eindringen konnte, floss seitlich über die Steuerkabelkanäle ab. Dies hob das lokale Erdpotenzial der Außenschränke über die dielektrische Festigkeit der Sekundärwicklungen der Stromwandler (ST) und Spannungswandler (VT) an. Die Gleichtakt-Potenzialdifferenz durchschlug die galvanische 10-kV-Isolierung der mikroprozessorbasierten IED-Schutzrelais im Schaltraum. Dies führte zum Ausbrennen der analogen Eingangs-/Ausgangsmodule und blockierte die zeitverzögerte Reserveauslösung (ANSI 50/51N), wodurch sich die Fehlerklärungszeit von 80 ms auf 1,2 Sekunden verlängerte und die destruktive Explosion des LeistungstrTransformators verursacht wurde.

Thermische Alterung und Schmelzen von Schirmen in Hochspannungskabeln

Bei VPE-isolierten 115-kV-Kabelstrecken mit beidseitig geerdeten Metallschirmen (Solid Bonding) verzerrt die fälschliche Annahme eines niedrigen spezifischen Widerstands in tiefen Schichten die Gegenseitigkeitsimpedanz der Schleife Erde-Schirm $\underline{Z}_{m(c-g)}$ gemäß den Formeln von 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 wobei D_e = 658,87 \sqrt{\frac{\rho_{eq}}{f}}

Beim Vorhandensein einer hochohmigen tiefen Schicht ( $\rho_2 \gg \rho_1$ ) vergrößerte sich die äquivalente Rückstromtiefe $D_e$ um mehr als 300 %. Dies erhöhte die induktive Reaktanz der Erdrückleitung drastisch und zwang den Unsymmterie- und Fehlerstrom fast vollständig über die metallischen Schirme der Leistungskabel zu fließen. Die thermische Stromdichte $JSchirm$ überschritt die zulässige Grenze für Kupfer nach IEC 60287 und IEC 60949:

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

Dies führte zur thermischen Delamination der äußeren leitfähigen Schicht, zur Erweichung des vernetzten Polyethylens (VPE) durch Übertemperaturen von mehr als 250 °C und zum anschließenden dielektrischen Phase-Erde-Durchschlag durch elektrothermisch beschleunigtes Water-Treeing.

Methodischer Simulations- und erweiterter Entwurfs-Workflow in Vexten Suite

Das durchgängige Engineering-Ökosystem Vexten Suite enthält eine deterministische, multiphysikalische Engine zur geoelektrischen Invertierung, gekoppelt mit Modulen zur Kurzschlussberechnung (IEC 60909 / IEEE 141), Leiter- und Kabeldimensionierung (IEC 60287 / NEC) und Erdungsanlagenoptimierung (IEEE 80 / IEC 60479). Nachfolgend ist die algorithmische Integration und das angewendete methodische Protokoll für kritische Infrastrukturprojekte dargestellt:

\begin{matrix} \boxed{\begin{array}{c} \textbf{1. Felddatenerfassung} \\ Wenner-Messung (a_i, R_{m,i}) \\ IEEE Std 81 \end{array}} & \longrightarrow & \boxed{\begin{array}{c} \textbf{2. Geoelektrische Invertierung} \\ Levenberg-Marquardt in Vexten \\ \mathbf{p} = [\rho_1, \rho_2, h]^T \end{array}} \\ \downarrow & & \downarrow \\ \boxed{\begin{array}{c} \textbf{4. Berechnung sicherer Spannungen} \\ Grenzwerte Etouch, Estep \\ IEEE 80 / IEC 60479 (50kg/70kg) \end{array}} & \longleftarrow & \boxed{\begin{array}{c} \textbf{3. Kurzschlussfehlersimulation} \\ Einspeisung I_k'' / I_g (IEC 60909) \\ Stromaufteilung S_f \end{array}} \\ \downarrow & & \downarrow \\ \boxed{\begin{array}{c} \textbf{5. Elektromagnetische Modellierung (BEM/FEM)} \\ Erdungsnetz im Zweischichtmedium \\ Berechnung von R_g, GPR , Vmesh \end{array}} & \longrightarrow & \boxed{\begin{array}{c} \textbf{6. Optimierung und Milderung} \\ Maschennetzanpassung, Tiefenerder \\ Kiesdicke h_s und Schirme \end{array}} \end{matrix}

Praktisches Beispiel für Simulation und numerische Invertierung in der Vexten Grounding Engine

Am Standort eines neuen industriellen Umspannwerks (138/13,8 kV) wird eine geoelektrische Erkundung nach der Wenner-Methode durchgeführt. Die Datenerfassung ergab folgende Feldwerte:

Elektrodenabstand $a$ [m]	Gemessener Widerstand $R_m$ [ $\Omega$ ]	Experimenteller scheinbarer spezifischer Widerstand $\rho_a^{mess}$ [ $\Omega\cdot m$ ]	Berechneter scheinbarer spezifischer Widerstand $\rho_a^{ber}$ [ $\Omega\cdot m$ ]	Relativer Restfehler [\%]
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

Bei der Verarbeitung der Datenmatrix im Modul für nichtlineare inverse Optimierung von Vexten Grounding konvergiert der Levenberg-Marquardt-Algorithmus in 7 Iterationen mit einem globalen RMS-Fehler von 1,32 % und ermittelt folgenden geschichteten Zustandsvektor:

\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

Ingenieurtechnische Analyse und Optimierung der Erdungsanlagen-Auslegung

Da der Reflexionskoeffizient $k$ stark negativ ist ( $k = -0,6812$ ), weist die tiefe Schicht einen spezifischen Widerstand auf, der mehr als fünfmal niedriger ist als der der Oberflächenschicht ( $\rho_2 \ll \rho_1$ ). Die rechnerische Optimierung in Vexten Grounding nutzt dieses physikalische Phänomen durch folgende Richtlinien aus:

Strategische Installation von Tiefenerdern: Es wird das Einbringen von kupferummantelten Stahlerdern (Copperweld) mit einer Mindestlänge von $l_v = 6,0\, m$ an den vier Ecken und kritischen Umfangsknoten vorgegeben. Durch das Überschreiten der Grenzschichttiefe ( $h = 4,35\, m$ ) dringen die vertikalen Elektroden $1,65\, m$ in die niederohmige Schicht ein ( $\rho_2 = 82,50\,\Omega\cdot m$ ). Dies verringert den Gesamtwiderstand des Systems $R_g$ um 42,7 % im Vergleich zu einer rein horizontalen Erdungsmasche mit gleicher Gesamtleiterlänge.
Optimierte Verminderung horizontaler Erdungsleiter: Dank der hohen Admittanz der unteren Schicht kann der Abstand zwischen den horizontalen Maschenleitern von einem dichten Netz von $5\, m \times 5\, m$ auf $10\, m \times 10\, m$ vergrößert werden. Dies führt zu erheblichen Kupfereinsparungen, ohne dass die nach IEEE Std 80 festgelegten Grenzwerte für Schritt- und Berührungsspannungen überschritten werden.
Integration mit dem Kurzschlussmodul der Vexten Suite (IEC 60909): Die reale Impedanz des geschichteten Bodens wird automatisch an die Kurzschlussberechnung der Vexten Suite übergeben. Dies ermöglicht die gekoppelte Neuberechnung der Reduktionsfaktoren für die Stromaufteilung über Erdseile und Phasenschirme ( $S_f$ ) und stellt sicher, dass der tatsächlich über die Erde zurückfließende Strom ( $I_g = S_f \cdot I_k''$ ) die Physik des umliegenden Übertragungs- und Verteilungsnetzes exakt widerspiegelt.

Die Umsetzung dieses ganzheitlichen methodischen Workflows garantiert höchste analytische Strenge, dielektrische Sicherheit für Personal und Anlagen sowie eine wirtschaftliche und betriebliche Optimierung beim Bau von Erdungsanlagen für Umspannwerke und Kraftwerke.