Schritt- und Berührungsspannungen in Umspannwerken nach IEEE 80

Elektrophysiologische Grundlagen und Kriterien des Kammerflimmerns

Das primäre Ziel der Erdungssystemplanung in Hoch- und Höchstspannungsunterstationen ist der Schutz von Menschenleben sowie die Unversehrtheit elektromechanischer Betriebsmittel. Bei einem Erdschlussereignis erzeugt die massive Stromeinspeisung in das Erdreich oberflächliche Potenzialgradienten, die den Körper eines Bedieners oder Passanten kritischen Potenzialdifferenzen aussetzen können.

Die physiologische Reaktion des Organismus auf den Durchgang von Wechselstrom bei Industrie-frequenz (50/60 Hz) wird durch die Auslösung von Kammerflimmern dominiert – ein Phänomen, bei dem die Ventrikelmyozyten ihre Kontraktionssynchronie verlieren, wodurch das Herzminutenvolumen aufgehoben wird. Die Wahrscheinlichkeit von Kammerflimmern hängt von der Stromstärke, der Einwirkdauer, dem Strompfad durch den Körper sowie der Phase des Herzzyklus ab, in welcher der elektrische Schlag stattfindet (insbesondere der vulnerablen Phase, die der T-Welle im Elektrokardiogramm entspricht).

I_B = \frac{k}{\sqrt{t_s}}

Im Rahmen der Norm IEEE Std 80 stellte Charles Dalziel empirisch fest, dass die vom menschlichen Körper tolerierbare Energie vor Erreichen der Wahrscheinlichkeitsschwelle von 0,5 % für Kammerflimmern durch eine Stoßenergiekonstante $S_b$ geregelt wird. Für eine Population mit einem Körpergewicht von ca. 50 kg beträgt die Konstante $k50 = 0{,}116$ , während sie für 70 kg $k70 = 0{,}157$ beträgt. Daher wird der maximal zulässige Körperstrom $I_B$ für eine Einwirkdauer $t_s$ (wobei 0{,}03 \le t_s \le 3{,}0 s ) analytisch wie folgt definiert:

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 ]

Andererseits verfolgt die Internationale Elektrotechnische Kommission in ihrem technischen Bericht IEC/TS 60479-1 einen nichtlinearen, probabilistischen Ansatz, der auf den Studien von Biegelmeier basiert. Anstatt über das gesamte Zeitfenster eine rein umgekehrt proportionale Beziehung zur Quadratwurzel der Zeit anzunehmen, definiert IEC 60479-1 Zeit-Strom-Zonen (Zonen AC-1 bis AC-4.3). Die Impedanz des menschlichen Körpers $Z_B$ ist in der IEC kein fixer Widerstand, sondern ein ohmsch-kapazitives Netzwerk, das von der Berührungsspannung ( $U_T$ ), dem Strompfad und dem Feuchtigkeitszustand der Haut abhängt:

Z_B(U_T) = Rint + \frac{Rp1(U_T)}{1 + j\omega Rp1 Cp1} + \frac{Rp2(U_T)}{1 + j\omega Rp2 Cp2}

Während IEEE Std 80 den menschlichen Körper konservativ auf einen konstanten, reinen ohmschen Widerstand von R_B = 1000\,\Omega vereinfacht (reihm-Fuß- oder Hand-Hand-Pfad mit festem Kontakt), wenden die Normen IEC 61936-1 und IEC 60479-1 dynamische Kurven an, bei denen die effektive Impedanz drastisch abfällt, sobald die Berührungsspannung 200 V übersteigt, und asymptotische Werte nahe dem inneren Körperwiderstand erreicht (Rint \approx 500 bis 750\,\Omega).

Mathematischer Rahmen und strenge Formulierung nach IEEE Std 80

Thevenin-Äquivalentschaltungen und Toleranzgrenzen

Zur Bestimmung der zulässigen Schrittspannungen ( $Vstep$ ) und Berührungsspannungen ( $Vtouch$ ) wird auf eine Thevenin-Äquivalentschaltung zurückgegriffen, betrachtet von den Kontaktpunkten des Körpers mit der Bodenoberfläche und den geerdeten Strukturen. Der vom Erdreich aus gesehen äquivalente Fußwiderstand wird modelliert, indem jeder Fuß als flache, leitende Scheibe mit dem Radius $b = 0{,}08 m$ betrachtet wird, die auf der Oberfläche eines homogenen Halbraums mit der scheinbaren Leitfähigkeit bzw. dem spezifischen Widerstand \rho_s aufliegt.

Der Eigenwiderstand eines Fußes im Kontakt mit dem Boden ergibt sich aus der Maxwellschen Formulierung:

Rfoot = \frac{\rho_s}{4b} \approx \frac{\rho_s}{4(0{,}08)} = 3{,}125 \rho_s \approx 3\rho_s

Unter der Bedingung der Berührungsspannung befinden sich beide Füße parallel in nahezu identischem Abstand zum Fehlerort, weshalb der äquivalente Widerstand R_{eq,touch} = \frac{Rfoot}{2} = 1{,}5\rho_s beträgt. Für die Bedingung der Schrittspannung befinden sich die Füße in Reihe, getrennt durch einen genormten Abstand von 1 Meter, was zu R_{eq,step} = 2 Rfoot = 6\rho_s führt.

Wird eine Deckschicht aus hochohmigem Material (typischerweise Schotter oder Splitt mit der Dicke $h_s$ und dem spezifischen Widerstand \rho_s) über einem nativen Boden mit dem spezifischen Widerstand \rho aufgebracht, muss ein Reduktions- oder Derating-Faktor C_s(h_s, K) eingeführt werden, welcher den Reflexionseffekt durch die dielektrische Diskontinuität ausgleicht:

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

Wobei der Reflexionskoeffizient $K$ zwischen dem nativen Boden und der oberen Schotterschicht ausgedrückt wird als:

K = \frac{\rho - \rho_s}{\rho + \rho_s}

Durch Integration dieser elektrodynamischen Herleitungen ergeben sich die nach IEEE Std 80 berechneten maximal zulässigen Potenziale für einen Operateur von 50 kg und 70 kg in den folgenden Grundgleichungen:

V_{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}}

V_{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}}

Erdererhöhung (GPR) und Maschen-Designparameter

Der Bemessungsstrom, der effektiv über das Erdnetz in das unendliche Erdreich abgeleitet wird ( $I_G$ ), entspricht nicht dem gesamten dreiphasigen oder einphasigen symmetrischen bzw. unsymmetrischen Kurzschlussstrom an den Klemmen ( $I_f$ ), sondern wird durch den Fehlerstrom-Teilungsfaktor ( $S_f$ ) und den Dekrementfaktor ( $D_f$ ) aufgrund der Asymmetrie der Gleichstromkomponente beeinflusst:

I_G = S_f \cdot D_f \cdot I_f

Der Dekrementfaktor $D_f$ quantifiziert den aperiodischen Anteil des Kurzschlussstroms in Abhängigkeit vom Reaktanz-Widerstands-Verhältnis ( $X/R$ ) der Thevenin-Impedanz am Fehlerort und der Abschaltzeit $t_f$ :

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

Die maximale Potenzialanhebung des Systems in Bezug auf einen entfernten Bezugspunkt wird als Ground Potential Rise ( $GPR$ ) bezeichnet und wie folgt berechnet:

GPR = I_G \cdot R_g

Der Ausbreitungswiderstand des Erdungssystems ( $R_g$ ) für komplexe, miteinander verbundene Netze mit Erdstäben oder tiefen vertikalen Elektroden wird mit hoher analytischer Genauigkeit über die verallgemeinerte Formel nach Sverak bestimmt:

R_g = \rho \left[ \frac{1}{L_T} + \frac{1}{\sqrt{20 A}} \left( 1 + \frac{1}{1 + h \sqrt{20 / A}} \right) \right]

Wobei $L_T$ die Gesamtlänge der vergrabenen Leiter (horizontaler Erdungskranz plus vertikale Elektroden) darstellt, $A$ die von der Masche eingenommene Fläche in Quadratmetern und $h$ die Verlegetiefe der horizontalen Gitterleiter.

Reale Maschen- und Schrittspannung in der Unterstation

Die reale Maschenspannung ( $E_m$ ), welche die ungünstigste Bedingung der Berührungsspannung in der Mitte des kritischsten äußeren Quadranten der Masche darstellt, wird rigoros formuliert durch das Produkt aus dem scheinbaren spezifischen Bodenwiderstand (\rho), dem geometrischen Abstandsfaktor $K_m$ , dem korrigierenden Unregelmäßigkeitsfaktor $K_i$ und der linearen Dichte des eingespeisten Stroms:

E_m = \frac{\rho \cdot I_G \cdot K_m \cdot K_i}{L_M}

Die dimensionslosen Faktoren $K_m$ und $K_i$ fassen die Topologie des Gitters zusammen und sind analytisch definiert als:

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]

K_i = 0{,}644 + 0{,}148 \cdot n

Wobei $D$ der Abstand zwischen parallelen Leitern, $d$ der Durchmesser des Maschenleiters, $h$ die Verlegetiefe, $n$ der geometrische Äquivalenzfaktor abgeleitet aus den Abmessungen und der Umfangsform (n = n_a \cdot n_b \cdot n_c \cdot n_d), K_h = \sqrt{1 + h/h_0} (mit $h_0 = 1 m$ ) und $Kii = 1$ für Maschen mit Umrandungselektroden oder Kii = 1/(2n)^{2/n} für Maschen ohne Erdstäbe ist.

Die gewichtete effektive Leiterlänge für die Berechnung der Maschenspannung ( $L_M$ ) berücksichtigt den Beitrag der vertikalen Erdungselektroden ( $L_r$ ):

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

Analog wird die berechnete Schrittspannung ( $E_s$ ) am Außenumfang der Unterstation unter Verwendung des geometrischen Schrittkoeffizienten $K_s$ bewertet:

E_s = \frac{\rho \cdot I_G \cdot K_s \cdot K_i}{L_S}

K_s = \frac{1}{\pi} \left[ \frac{1}{2h} + \frac{1}{D + h} + \frac{1}{D} \left( 1 - 0{,}5^{n-2} \right) \right]

L_S = 0{,}75 \cdot L_c + 0{,}85 \cdot L_R

Transientes Verhalten und Blitzeinwirkung nach IEC 62305

Impuls-Impedanz vs. Netzfrequenz-Widerstand

Das Ansprechverhalten eines Erdungssystems auf die Injektion von Blitzströmen, die durch die Normenreihe IEC 62305 geregelt werden, unterscheidet sich grundlegend vom Verhalten bei 50/60 Hz. Direkte oder induzierte atmosphärische Entladungen weisen ultra-schnelle, genormte Wellenfronten auf (Wellenform 10/350 \mu s für den ersten Blitzschlag nach IEC 62305-1 und 8/20 \mu s für nachfolgende Einschläge) mit Frequenzanteilen, die sich über den Bereich von Hunderten von Kilohertz bis Megahertz erstrecken.

Bei diesen Frequenzen wird die verteilte longitudinale Induktivität ( $L'$ ) der Kupfer- oder feuerverzinkten Stahlleiter zum dominierenden Faktor gegenüber der Querleitfähigkeit des Bodens ( $G'$ ). Infolgedessen hat der Strom physikalisch keine Zeit, sich über die gesamte Maschengeometrie zu verteilen, sondern konzentriert sich auf die unmittelbare Umgebung des Einschlagspunkts.

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

Die effektive Länge des vergrabenen Leiters $l_e$ , jenseits derer jede weitere Verlängerung der Elektrode die scheinbare Impuls-Impedanz für eine Frontzeit $T_1$ nicht mehr reduziert, wird nach IEC 62305-3 analytisch ausgedrückt als:

l_e = \sqrt{\frac{\rho \cdot T_1}{\mu_0}} \approx 1{,}6 \cdot \sqrt{\rho \cdot T_1}

Wobei \mu_0 = 4\pi \times 10^{-7} H/m die magnetische Feldkonstante des Vakuum und $T_1$ die Frontzeit des Impulses in Mikrosekunden (\mu s ) ist.

Durch Blitzschlag induzierte Schritt- und Berührungsspannungen

Unter der Einwirkung eines Impulsblitzstroms i(t) erzeugt der radiale Potenzialgradient an der Bodenoberfläche extreme transiente Spannungen, die Durchschläge in der Luft entlang der unteren Extremitäten des menschlichen Körpers verursachen können (Oberflächenüberschlag oder Flashover). Die transiente Impuls-Berührungsspannung ( $U_{touch,imp}$ ) und die transiente Schrittspannung ( $U_{step,imp}$ ) gemäß IEC 62305-3 Anhang E koppeln sich über rein induktive und ohmsche Komponenten ein:

U_{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)

Bei Blitzschlägen in die Abschirmungsstruktur (Fangstangen, Erdseile oder Portale) erzeugen die hohen zeitlichen Ableitungen (di/dt > 100 kA/ \mu s ) parasitäre induktive Potenziale in den Ableitungen, die die dielektrische Festigkeit der Luft überschreiten, was sichere Trennungsabstände ( $s$ ) erfordert, die nach folgender Formel berechnet werden:

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

Wobei $k_i$ der Koeffizient in Abhängigkeit vom Blitzschutzniveau (LPL I bis IV) ist, $k_c$ der geometrische Aufteilungsfaktor des Stroms in den Ableitungen, $k_m$ der Koeffizient des Isoliermaterials (Luft oder fest) und $l$ der lineare Abstand entlang der Ableitung bis zum nächstgelegenen Potenzialausgleichspunkt.

Normative Vergleichsmatrix: IEEE Std 80 vs. IEC 60479-1 / IEC 61936-1 / IEC 62305-3

Parameter / Designkriterium	IEEE Std 80 (2013)	IEC 60479-1 / IEC 61936-1 (2021)	IEC 62305-3 (Blitzschutz)
Körpergewichtsmodell	Diskret: 50 kg ( $k=0{,}116$ ) oder 70 kg ( $k=0{,}157$ )	Statistische Population (Perzentil 5 %, 50 %, 95 %)	Nicht anwendbar (Risikobewertung durch Impulsstoß)
Körperimpedanz ( $R_B$ / $Z_B$ )	Konstant: R_B = 1000\,\Omega (rein ohmsch)	Nichtlinear: Z_B = f(U_T), 500 bis 1000\,\Omega	R_B \approx 0\,\Omega (Lichtbogen bei Hochfrequenz)
Kritischer Zeitbereich	Industriefrequenz (50/60 Hz), 0{,}03 \le t_s \le 3{,}0 s	Industriefrequenz (50/60 Hz), 0{,}01 \le t_s \le 10{,}0 s	Transstanter Impuls (10/350\,\mu s , 8/20\,\mu s )
Oberflächenschichtwiderstand	Formales analytisches Derating C_s(h_s, K)	Ähnlicher oder empirischer Oberflächenschichtfaktor	Isolierende Steuerung: Asphalt $5 cm$ (\ge 100 k \Omega\cdot m )
Berührungsspannungsgrenze	Vtouch = (1000 + 1{,}5 C_s \rho_s) \frac{k}{\sqrt{t_s}}	Definiert durch Grenzkurve U_{v,tol} = f(t_f, Z_B, Rfoot)	Utouch \le 100 kV (mit Potenzialausgleich)
Fehlerteilungsfaktor ( $S_f$ )	Explizite Berechnung mit Erdseilen und Leitern	Reduktionsfaktor $r$ (Funktion von Schirmen und Erdseilen)	Teilungsfaktor $k_c$ basierend auf Ableitungstopologie
Bodenverhalten	Mehrschicht-Widerstandsverlauf bei 50 Hz	Zweischicht- und homogener Halbraummodelle	Bodenionisationseffekt (E_0 \approx 300 kV/m )

Forensische Fehleranalyse und Potenzialverschleppungsphänomene

Ausbreitungsmechanismen verschleppter Potenziale

Einer der verheerendsten Ausfallmodi in der forensischen Unterstationsingenieurskunst entspricht dem unabsichtlichen Export oder Import der Erdererhöhung ( $GPR$ ) über die Umzäunungsgrenzen der Anlage hinaus. Wenn es innerhalb der Unterstation zu einem Erdschluss kommt, hebt sich das gesamte Erdnetz auf ein Potenzial $V = GPR$ an, bezogen auf das Null-Referenzpotenzial in der fernen Erde.

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

Die Hauptvektoren für die galvanische Kopplung verschleppter Potenziale umfassen:

Neutralleiter von Mittel- und Niederspannungsnetzen: Wenn der Neutralleiter der Sekundärwicklung des Eigenbedarfs transformators oder einer ländlichen Verteilungsleitung starr mit dem internen Erdnetz verbunden ist und mit mehrfacher Erdung aus der Unterstation herausgeführt wird, wirkt sich das vollständige $GPR$ direkt auf die Isolatoren der Endnutzernetze aus, was zu Explosionen von Hauszählern, Durchschlägen in Haushaltsgeräten und der Gefahr von Massenelektrokutionen führt.
Metallische Schirme von Leistungskabeln und Telekommunikationsleitungen: An beiden Enden geerdete Schirme von Mittelspannungs-Erdkabeln übertragen den Erdfehlerstrom zu fernen Endverschlüssen, was Kreisströme von Hundertern von Ampere verursacht, die den äußeren Polyethylenmantel thermisch zerstören (thermischer Fehler durch $I^2 t$ ) und Überspannungen induzieren, welche die digitalen Schnittstellen der Leittechnik zerstören.
Wasser- oder Gasrohrleitungen sowie Eisenbahnschienen: Kontinuierliche, im Erdreich vergrabene metallische Elemente ohne dielektrische Trennstellen wirken als Elektroden unendlicher Ausdehnung und setzen Bediener in Pumpstationen oder fernen Bahnübergängen dem vollen $GPR$ aus, wobei die verschleppte Berührungsspannung V_{touch,trans} \approx GPR den Dämpfungseffekt des Schotterbettes der Unterstation vermissen lässt.

Verhalten in gasisolierten Schaltanlagen (GIS)

In mit Schwefelhexafluorid isolierten Schaltanlagen ( $SF _6$ ) erzeugt die Betätigung von Trennern oder das Auftreten von Überschlägen sehr schnelle transiente Überspannungen (VFT, Very Fast Transients) mit Anstiegszeiten im Nanosekundenbereich ( $t_r < 10 ns$ ). Diese Wellenfronten vermögen aufgrund der parasitären induktiven Reaktanz der Verbindungsschienen (L \approx 1\,\mu H/m ) nicht über herkömmliche Erdungsverbindungen zur Erde abzufließen.

Folglich entsteht das Phänomen der transienten Gehäusepotenzialanhebung (TEV, Transient Enclosure Voltage), bei dem das äußere Aluminiumgehäuse der GIS Potenzialschwingungen von bis zu mehreren Dutzend Kilovolt in Bezug auf die metallischen Tragstrukturen und umliegenden Böden erfährt, was sekundäre Mikrolichtbögen erzeugt, die Steuerkabel des digitalen Leitsystems (SAS) zerstören und Entladungen durch Berührung für Anlagenprüfer verursachen.

Methodik der Bodenstratifizierung und Maschenoptimierung

Inversion von Vertikalsondierungskurven (VES)

Die Genauigkeit der Berechnungen von $V_m$ , $V_s$ und $R_g$ ist streng an die Formulierung des geoelektrischen Untergrundmodells gebunden. Die Prospektion mittels der Wenner-Vierpunktmethode besteht darin, Gleichstrom oder Niederfrequenzstrom über die äußeren Elektroden ( $C_1, C_2$ ) einzuspeisen und den Spannungsabfall an den inneren Elektroden ( $P_1, P_2$ ), die symmetrisch mit einem Abstand $a$ angeordnet sind, zu messen:

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

Da der Boden selten homogen ist, variiert der scheinbare spezifische Widerstand \rho_a als Funktion der Eindringtiefe der Stromlinien, die proportional zum Abstand $a$ ist. Um eine Zweischichtstruktur (obere Schicht mit spezifischem Widerstand \rho_1 und Dicke $h$ über einem tiefen Halbraum mit spezifischem Widerstand \rho_2) zu modellieren, wird die Stefanescu-Integralgleichung mittels der Hankel-Transformation nullter Ordnung analytisch gelöst:

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

Wobei der geoelektrische Reflexionsfaktor wie folgt lautet:

K = \frac{\rho_2 - \rho_1}{\rho_2 + \rho_1}

Der numerische Optimierungsalgorithmus (gedämpfte nichtlineare kleinste Quadrate nach Levenberg-Marquardt) zielt darauf ab, die quadratische RMS-Fehlerzielfunktion zu minimieren:

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

Erweiterte Topologien zur Minderung oberflächlicher Gradienten

Wenn der spezifische Oberflächenwiderstand ungünstig ist (\rho_1 \gg \rho_2) oder der Fehlerstrom $I_G$ massiv ausfällt, zeigt die einfache Erhöhung uniformer Leiter im Gitter abnehmende Erträge. Zu den fortgeschrittenen Techniken der elektromagnetischen Optimierung gehören:

Abgestufte Verdichtungmasche (Perimeterverdichtung): Exponentielle Verringerung des Abstands zwischen den Leitern bei Annäherung an die Außenbegrenzung der Unterstation. Dies vergleichmäßigt die pro Längeneinheit ausgestoßene Stromdichte und glättet die peripheren Potenzialspitzen, welche die äußere Berührungs- und Schrittspannung auslösen.
Tiefe vertikale Erder an Perimeter und Ecken: Das Einbringen tiefer Erdstäbe (15 bis 30 Meter), die an den Ecken der Masche abgeteuft werden, ermöglicht die bevorzugte Ableitung des Stroms in die tieferen, leitfähigeren Schichten (\rho_2), wodurch das Gesamte $GPR$ abgesenkt und der Oberflächengradient an den kritischen Rändern substanziell reduziert wird.
Hochleistungs-Dielektrikumsschicht: Aufbringen von Schichten aus gewaschenem Granitsplitt (\rho_s \ge 3000\,\Omega\cdot m ) oder Asphaltdecken (\rho_s \ge 10000\,\Omega\cdot m ) mit einer Dicke von $h_s = 0{,}10 bis 0{,}15 m$ . Es ist strikt zu verifizieren, dass der Reduktionsfaktor $C_s$ seine Wirksamkeit über die gesamte Lebensdauer der Anlage beibehält, wobei Schlamm- und Unkrautverschmutzungen durch dichte Geotextilien vermieden werden.

Implementierung und computergestützte Modellierung in der Vexten Suite

Integration von Kurzschlüssen und parametrische Auswertung

Im ingenieurtechnischen Entwurfsworkflow innerhalb der integrierten Plattform Vexten Suite wird die Dimensionierung der Erdung nativ mit den Berechnungsmotoren für dreiphasige und unsymmetrische Kurzschlüsse nach den Normen IEC 60909 und IEEE 141 (Red Book) synchronisiert. Die exakte Bestimmung der Nullsystem-Stromkomponente ( $3I_0$ ) und der Thevenin-Äquivalenzimpedanz am Knoten der Unterstation speist direkt den Erdungsoptimierungsmotor.

Die folgende Matrizenformulierung wird in der Vexten Suite zur Lösung der induktiven und galvanischen Kopplung mehrerer Freileitungen mit Erdseil (OPGW/Stahl), die an die Unterstation angeschlossen sind, implementiert, was die Bestimmung des realen Teihungsfaktors $S_f$ mit absoluter Strenge ermöglicht:

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

D-Mapping von Oberflächengradienten und normativer Nachweis

Der elektromagnetische Motor der Vexten Suite unterteilt das Erdnetz und die vertikalen Elektroden in diskrete zylindrische Segmente und wendet die Momentenmethode (MoM) an, um die longitudinale Ladungsdichte $q_j$ und den radialen Leckstrom in jedem Segment zu berechnen. Die Potenzialverteilung an jeder beliebigen Oberflächenkoordinate (x, y, 0) wird mittels numerischer Integration analytisch ausgewertet:

\Phi(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'

Sobald das dreidimensionale Potenzialfeld berechnet ist, verarbeitet die Vexten Suite automatisch die Anlagentopologie durch Überlagerung von Schaltgeräten, Leistungstransformatoren, Umzäunungen und Betriebswegen. Das System generiert räumliche Vektoren von 1 Meter in alle radialen Richtungen ausgehend von den berührbaren massiven Leitern, um Punkt für Punkt abzugleichen:

\max\left( \Phistructure - \Phi(xfoot, yfoot, 0) \right) \le V_{touch,adm}(t_s)

\max\left( \Phi(x+1, y, 0) - \Phi(x, y, 0) \right) \le V_{step,adm}(t_s)

Jede Überschreitung der normativen Toleranzgrenzen nach IEEE Std 80 oder IEC 61936-1 wird in der grafischen Umgebung mittels hochauflösender Farbgradienten (Heatmaps) hervorgehoben, was es dem Ingenieur ermöglicht, interaktiv Hilfserdstäbe, periphere Potenzialausgleiche oder Modifikationen der Schotterschichtdicke vorzunehmen, bis die vollständige Konformität der Anlage gegenüber den strengsten Anforderungen an die Personensicherheit und operative Robustheit zertifiziert ist.