Analyse von Eigenfrequenzen und elektromagnetischen Transienten in Verteilungsnetzen

Physikalisch-mathematische Grundlagen der elektromagnetischen Ausgleichsvorgänge und des Frequenzgangs

Die rigorose Analyse elektromagnetischer Ausgleichsvorgänge und Eigenfrequenzen in elektrischen Verteilungsnetzen erfordert eine fortgeschrittene mathematische Formulierung, die über die klassische zeigerbasierte Stationärwertanalyse hinausgeht. Moderne Verteilungsnetze, die durch eine hohe Durchdringung von Mittelspannungs-Erdkabeln, Transformatoren mit nichtlinearen ferromagnetischen Kernen und dezentralen Leistungselektroniklasten gekennzeichnet sind, verhalten sich wie verteilte und konzentrierte Parametersysteme, die schweren Impuls- und Oberschwingungsanregungen ausgesetzt sind.

Zur Modellierung der Ausbreitung elektromagnetischer Wellen entlang der Netzkleiter müssen die gekoppelten partiellen Differentialgleichungen der Telegrafengleichungen (Telegrapher's equations) herangezogen werden. Unter Berücksichtigung einer symmetrischen oder asymmetrischen Mehrleitersystemleitung der infinitesimalen Länge $dx$ wird das Momentanverhalten von Spannungen v(x,t) und Strömen i(x,t) durch das folgende Matrizensystem definiert:

\frac{\partial}{\partial x} \mathbf{v}(x,t) + \mathbf{R} \mathbf{i}(x,t) + \mathbf{L} \frac{\partial}{\partial t} \mathbf{i}(x,t) = \mathbf{0}

\frac{\partial}{\partial x} \mathbf{i}(x,t) + \mathbf{G} \mathbf{v}(x,t) + \mathbf{C} \frac{\partial}{\partial t} \mathbf{v}(x,t) = \mathbf{0}

Wobei $\mathbf{R}$ , $\mathbf{L}$ , $\mathbf{C}$ und $\mathbf{G}$ die längenbezogenen Matrizen für Widerstand, Induktivität, Kapazität bzw. Querleitwert (dielektrische Ableitung) darstellen. Durch Transformation dieser Gleichungen in den Frequenzbereich mittels Laplace- oder Fourier-Transformation erhält man die modale Ausbreitungsgleichung:

\frac{d^2}{dx^2} \mathbf{V}(x,s) = \mathbf{Z}(s) \mathbf{Y}(s) \mathbf{V}(x,s) = \mathbf{P}^2(s) \mathbf{V}(x,s)

Wobei Z(s) = R + sL, Y(s) = G + sC und P(s) die Ausbreitungsmatrix ist. Die Eigenwerte der Produktmatrix Z(s)Y(s) bestimmen die modalen Ausbreitungskonstanten \gamma_k(\omega) = \alpha_k(\omega) + j\beta_k(\omega), welche die Dämpfung und Phasengeschwindigkeit der Wandelwellen steuern. Die Eigenfrequenzen des Systems (\omega_n) entsprechen den Polen der Eingangsimpedanz-Übertragungsfunktion Zin(s), gesehen von jedem beliebigen Knoten des Verteilungsnetzes aus, durch Lösung der charakteristischen Bedingung:

\det \left[ \mathbf{I} - \mathbf{M}_{end}(s) \right] = 0

Wobei M_{end}(s) die Reflexions- und Übertragungsmatrix an den Netzabschlüssen darstellt, die die Impedanzen von Transformatoren, Lasten und Überspannungsableitern (SPD) koppelt.

Klassifizierung und Ursprung elektromagnetischer Ausgleichsvorgänge in Mittel- und Niederspannungsnetzen

Elektromagnetische Ausgleichsvorgänge in Verteilungsnetzen entstehen durch Betriebs- und atmosphärische Störungen, die Breitband-Frequenzspektren einspeisen. Die korrekte Kategorisierung dieser Phänomene ist für den Entwurf koordinierter Isolationsschemata gemäß den Normen IEEE Std 1313.2 und IEC 60071-2 unerlässlich.

Ausgleichsvorgänge atmosphärischen Ursprungs (Blitzeinschläge)

Direkte oder indirekte atmosphärische Entladungen in Freileitungsnetze induzieren Überspannungswellen mit extrem schnellen Wellenfronten (in der Größenordnung von 0.5\,\mu s bis 5\,\mu s ) und langen Schwänzen (bis zu 50\,\mu s ). Der Blitzstrom wird typischerweise durch die doppelte Heidler-Exponentialfunktion modelliert:

i(t) = \frac{I_p}{\eta} \cdot \frac{(t/\tau_1)^n}{1 + (t/\tau_1)^n} \cdot e^{-t/\tau_2}

Wobei $I_p$ die Scheitelstromamplitude, \tau_1 die Zeitkonstante der Wellenfront, \tau_2 die Zeitkonstante des Wellenschwanzes, $n$ der Steilheitsfaktor (steepness factor) und \eta ein Amplitudenkorrekturfaktor ist. Beim Eindringen dieser Wellen von Freileitungsabschnitten in Erdkabel sinkt die Wellenimpedanz drastisch (typischerweise von Z_c \approx 400\,\Omega auf Z_c \approx 30\,\Omega), was eine Spannungsbrechung verursacht, die die Amplitude der einfallenden Welle am Übergabepunkt verdoppeln kann, wenn die Kabellegelänge kleiner als die kritische Frontwellenlänge ist.

Schalthandlungen beim Öffnen und Schließen von Leistungsschaltern

Der Betrieb von Vakuumleistungsschaltern (vacuum circuit breakers) oder Lasttrennschaltern erzeugt hochfrequente Ausgleichsvorgänge aufgrund des Phänomens mehrfacher Lichtbogenwiederzündungen (multiple restrikes). Wenn sich der Kontakt trennt, verläuft der Strom durch einen hochfrequenten Nulldurchgang, und die dielektrische Festigkeit des Löschmediums stellt sich schneller wieder her als die Anstiegsgeschwindigkeit der transienten Wiederkehrspannung (TRV). Dies führt zu einer abrupten Stromunterbrechung (current chopping), deren gespeicherte induktive Energie gemäß der Energiebeziehung auf die parasitäre Netzkapazität übertragen wird:

\frac{1}{2} Leq ichop^2 = \frac{1}{2} Ceq vmax^2 \implies vmax = ichop \sqrt{\frac{Leq}{Ceq}}

Diese oszillierende Hochfrequenzüberspannung (typischerweise zwischen 1\, kHz und 100\, kHz ) übt eine schwere dielektrische Beanspruchung auf die Hauptisolation von Transformatorwicklungen aus.

Resonanz und Ferroresonanz in Verteilungssystemen

Oberschwingungsresonanz und Ferroresonanz stellen zwei der destruktivsten und komplexesten Fehlerarten in unterirdischen und gemischten Verteilungsnetzen dar. Ihre Analyse erfordert die Auswertung des Zusammenspiels zwischen linearen reaktiven Elementen und ferromagnetischen Nichtlinearitäten.

Serien- und Parallelresonanz

Die Parallelresonanz (Antiresonanz) tritt auf, wenn die äquivalente induktive Reaktanz des Netzes (beigetragen durch Verteilungstransformatoren und Dämpfungsdrosseln) bei einer spezifischen Oberschwingungsfrequenz $h$ gleich der kapazitiven Reaktanz der Erdkabel und Blindleistungskompensations-Kondensatorbänke wird:

\omega_r = \frac{1}{\sqrt{Leq Ceq}} \implies h_r = \frac{f_r}{f_1} = \sqrt{\frac{X_c}{X_L}}

In diesem Zustand erreicht die von der Quelle aus gesehene Impedanz einen theoretischen Maximalwert, der nur durch den äquivalenten Netzwiderstand begrenzt wird. Jede Einspeisung von Oberschwingungsströmen (aus nichtlinearen Lasten wie Frequenzumrichtern und Gleichrichtern) nahe $h_r$ erzeugt massive Oberschwingungs-Spannungsabfälle, Überhitzung durch den Joule-Effekt und eine schwerwiegende Verformung der Spannungswelle.

Ferroresonanz in Transformatoren

Ferroresonanz ist ein nichtlineares Phänomen, das durch anhaltende niederfrequente Überspannungen (Grundschwingung, Subharmonische oder chaotisch), extreme Oberschwingungsverzerrungen und hohe Ströme gekennzeichnet ist. Sie tritt typischerweise auf, wenn ein Transformator mit ungesättigtem Eisenkern im Leerlauf über eine Seriencapazität erregt wird (z. B. die Kapazität eines an einem Ende getrennten Erdkabels oder offene Phasen durch einpoligen Sicherungsbetrieb).

Das nichtlineare Verhalten des magnetkrafterfüllten Kerns wird durch den verketteten Fluss \lambda(i) oder die konstitutive Beziehung der nichtlinearen Induktivität modelliert:

v(t) = R i(t) + \frac{d\lambda(i)}{dt}

Die Koexistenz mehrerer stabiler Zustände (Grundschwingungsregime und subharmonische Ferroresonanzregime) wird durch die Bifurkationstheorie in nichtlinearen dynamischen Systemen erklärt. Die Differentialgleichung, die den durch eine sinusförmige Quelle V_m \sin(\omega t) angeregten nichtlinearen LC-Kreis regelt, lautet:

L(i) \frac{d^2 i}{dt^2} + R \frac{d i}{dt} + \frac{1}{C} \ int i \, dt = V_m \omega \cos(\omega t)

Wobei die Induktivität L(i) = \frac{d\lambda}{di} beim Eintritt in den magnetischen Sättigungsbereich drastisch variiert, den Wert von $L$ verringert und die Resonanzfrequenz des Kreises zu Werten verschiebt, die mit der Netzhäufigkeit oder ihren Subharmonischen übereinstimmen.

Forensische Fehleranalyse und dielektrische Konsequenzen

Hochfrequente elektromagnetische Ausgleichsvorgänge und lang anhaltende Resonanzen verschlechtern die kritischen Betriebsmittel des Verteilungsnetzes kumulativ und katastrophal. Nachfolgend wird die forensische Analyse der Fehler- und Versagensmechanismen nach Komponenten aufgeschlüsselt:

Verteilungs- und Leistungstransformatoren

Wenn eine schnelle Wellenfront (wie sie durch Blitzeinschläge oder Schalthandlungen erzeugt wird) in die Hochspannungswicklung eines Transformators eindringt, ist die anfängliche Spannungsverteilung entlang der Windungen nicht linear. Aufgrund des Vorhandenseins der Serienkapazität zwischen den Windungen ( $C_s$ ) und der Erdkapazität der Wicklung ( $C_g$ ) bestimmt der Verteilungsfaktor \alpha = \sqrt{C_g / C_s} einen extrem hohen Spannungsgradienten in den ersten Scheiben oder Windungen der Wicklung.

Dies erzeugt interne Teilentladungen (Partial Discharges - PD) im ölgetränkten Papier, eine Durchschlagen der zellulosigen Isolation und Windungsschlüsse. Thermisch rufen die durch Parallelresonanz induzierten Oberschwingungsströme Zusatzverluste in den Wicklungen durch den Skin-Effekt (Stromverdrängung) und den Nähe-Effekt (proximity effect) hervor, berechnet über den zusätzlichen Verlustfaktor $H$ :

PLL = PDC \left[ 1 + \sum_{h=2}^{n} \left( \frac{I_h}{I_1} \right)^2 H(h) \right]

Mittelspannungs-Erdkabel

Mit vernetztem Polyethylen (VPE / XLPE) isolierte Kabel sind gegenüber anhaltenden Spannungstransitoren und Ferroresonanz äußerst anfällig. Anhaltende Überspannungen beschleunigen den Prozess des elektrischen Treeing (elektrische Baumstrukturen), einem Degradationsmechanismus der Feststoffisolation, bei dem Mikrokohlensäurehohlräume unter intensiven elektrischen Feldern (> 10\, kV/mm ) kontinuierliche Teilentladungen erzeugen, die dendritische Kanäle bis zum vollständigen dielektrischen Durchschlag (Zerstörung der Primärisolation) ausdehnen.

Schalt- und Schutzgeräte

Vakuumleistungsschalter und Wiedereinschalter leiden unter starker Erosion der Hauptkontakte aufgrund der in hochfrequenten Lichtbögen angesammelten Energie bei Unterbrechungsoperationen mit Fehlangepassten induktiven oder kapazitiven Lasten. Ebenso können Zinkoxid-Überspannungsableiter (MOV) bei Überschreitung ihrer Nennenergie-Dissipationskapazität $Wth$ durch einen energetischen Lawineneffekt thermisch versagen.

Designparameter, normative Grenzwerte und Minderungstrategien

Zur Gewährleistung der Zuverlässigkeit und Isolationskoordination gemäß internationalen Richtlinien müssen Ingenieure strenge präventive Designkriterien anwenden und die Auswirkungen der Eigenfrequenzen mindern.

Die folgende Tabelle fasst die kritischen elektrischen Parameter, normativen Grenzwerte nach IEEE- und IEC-Standards, kritische Ausfallbedingungen und die damit verbundenen operativen Konsequenzen zusammen:

Elektrischer Parameter / Phänomen	Referenzstandard (IEEE / IEC)	Normativer Grenzwert / Kritischer Schwellenwert	Kritische Ausfallbedingung	Operative und dielektrische Konsequenzen
Vorübergehende Überspannungen (TOV)	IEEE Std 1313.2 / IEC 60071-2	VTOV \le 1.7 \times V_{ph-n} (für t \le 1\, s )	Anhaltende Ferroresonanz oder Lastabwurf in schwachen Netzen	Thermische Degradation und vorzeitige Alterung der VPE-Kabelisolation.
Oberschwingungs-Spannungsgesamtfestigkeit (THD-V)	IEEE Std 519-2022 / IEC 61000-2-2	THD-V \le 5.0\% (Mittelspannungsnetze)	Parallelresonanz nahe der 5., 7. oder 11. Oberschwingung	Transformatorüberhitzung, Fehlauslösungen von Schutzeinrichtungen und dielektrische Ermüdung.
Spannungsanstiegsgeschwindigkeit ( $dv/dt$ )	IEC 60034-18-41 / NEMA MG1	dv/dt \le 5\, kV /\mu s (an Motoren-/Kabelklemmen)	Betrieb von Vakuumschaltern mit mehrfachen Wiederzündungen	Durchschlag der Windungsisolation und Zwischenwindungsfehler in Transformatoren und rotierenden Maschinen.
Energieabsorptionsvermögen von SPD (MOV)	IEEE Std C62.11 / IEC 61643-11	Nennabsorption \ge 2.5\, kJ/kV der Bemessungsspannung	Mehrfache Blitzeinschläge oder verstärkte reflektierte Wellenfronten	Katastrophaler thermischer Ausfall (Explosion) des Ableiters und Schutzverlust.
Oberschwingungs-Derating-Faktor für Kabel	IEC 60287 / NEC 310.15	Reduktionsfaktor $Kharm$ gemäß Stromspektrum	Nullleiterströme höher als Phasenströme durch Mit-/Nullsystemkomponenten (3., 9.)	Isolationsschmelze durch übermäßige Leitertemperatur ( $> 90 °C$ für VPE).

Minderungsstrategien und aktive/passive Kompensation

Zur Dämpfung resonanter Eigenfrequenzen und schwerer Transienten werden folgende ingenieurtechnische Gegenmaßnahmen eingesetzt:

Installation abgestimmter Oberschwingungsfilter (Passive LC-Filter): Sie werden parallel zum Netz geschaltet, um kritische Oberschwingungen (typischerweise die 5. und 7.) zu Entladewiderständen abzuleiten, wodurch die äquivalente Netzimpedanz modifiziert und Parallelresonanz vermieden wird. Ihre Abstimmung erfolgt knapp unterhalb der Ziel-Oberschwingungsfrequenz:
$fabgestimmt = \frac{1}{2\pi \sqrt{Lfilter Cfilter}} < foberschwingung$
Dämpfungsdrosseln und Voreinschaltwiderstände: In Hochspannungs- und kritischen Mittelspannungsschaltern dämpft der Einsatz von Voreinschaltwiderständen die transiente Wellenfront, wodurch der Überspannungsfaktor auf $vmax / Vph < 2.0$ reduziert wird.
Strategisch platzierte Überspannungsableiter auf Zinkoxidbasis (MOV): Sie werden an den Übergabepunkten zwischen Freileitungen und Erdkabeln sowie an den Transformatorbuchsen installiert, um transiente Spannungsspitzen zu kappen und die einfallende Energie zu begrenzen.

Praktische Anwendung und Analyse mittels Vexten Suite

Die fortgeschrittene Modellierung von Verteilungsnetzen zur Untersuchung elektromagnetischer Ausgleichsvorgänge und Eigenfrequenzen erfordert hochpräzise Rechenplattformen wie die Vexten Suite. Nachfolgend werden die im Softwarepaket implementierten algorithmischen Methoden und Engineering-Verfahren zur Erfüllung internationaler Standards dargelegt.

Kurzschluss- und Transienten-Berechnungsablauf (IEC 60909 / IEEE 141)

Die Transienten-Berechnungsengine der Vexten Suite löst das System der Differentialgleichungen im Zeitbereich unter Verwendung des Trapez-Integrationsverfahrens mit adaptivem Zeitschritt (\Delta t < 10\,\mu s ). Das integrierte analytische Verfahren umfasst folgende Schritte:

Topologische Erfassung und Modellierung: Das System importiert physikalische Kabelparameter (Geometrie, Halbleiterschichten, VPE-Isolation), Transformatoren (Kurzschlussimpedanzen, Leerlaufstrom, Fröhlich-Kennelly-Sättigungskurve) und Quellen dezentraler Erzeugung.
Eigenfrequenzanalyse (Impedanz-Scan): Eine Oberschwingungsstromquelle mit Einheitsamplitude und variabler Frequenz f \in [10\, Hz , 100\, kHz ] wird in kritische Netzwerkknoten eingespeist. Vexten berechnet die Netzknotenimpedanz Znn(j\omega) und generiert das Bode-Diagramm von Betrag und Phase, wobei Parallelresonanzspitzen und Serienresonanztäler automatisch identifiziert werden.
Schalt- und Blitzsimulation (EMTP-artiger Solver): Schaltvorgänge (Schalteröffnungen mit spezifischen Abschaltströmen) und genormte Blitzwellen (8/20\,\mu s und 1.2/50\,\mu s ) werden konfiguriert. Die Software bewertet die Wellenausbreitung unter Berücksichtigung der frequenzabhängigen modalen Dämpfung und des Skin-Effekts über die modifizierte Carson-Formulierung oder Bessel-Funktionen für zylindrische Leiter.
Thermische Derating-Bewertung von Kabeln (IEC 60287): Basierend auf dem Oberschwingungs-Stromspektrum aus der Oberschwingungs-Lastflussanalyse berechnet Vexten automatisch den Reduktionsfaktor für die Strombelastbarkeit ( $ampacity$ ) des Erdkabels und stellt sicher, dass die maximale Leitertemperatur unter verzerrten Bedingungen den normativen thermischen Grenzwert von $90 °C$ für VPE nicht überschreitet.

Somit ermöglicht die Vexten Suite dem Fachingenieur eine umfassende Diagnose, verhindert katastrophale Resonanzausfälle und garantiert eine optimale Isolationskoordination in komplexen Verteilungsnetzen.