Ferroresonanz und magnetische Sättigung in Shunt-Drosseln bei einphasiger Schaltung

Einführung und Elektromagnetische Grundlagen von Querinduktivitäten (Drosselspulen)

Querinduktivitäten (Kompensationsdrosseln), die in Höchstspannungs- (HGÜ/UHV) und Höchstspannungs-Übertragungssystemen eingesetzt werden, spielen eine entscheidende Rolle bei der Kompensation der kapazitiven Blindleistung, die von langen Übertragungsleitungen unter Schwachlast- oder Leerlaufbedingungen erzeugt wird. Elektromechanisch bestehen diese Betriebsmittel aus ferromagnetischen Kernen mit verteilten Luftspalten, die so ausgelegt sind, dass sie unter Nennbedingungen und moderaten vorübergehenden Überspannungen eine lineare Magnetisierungskennlinie aufrechterhalten. Wenn jedoch einphasige Schaltvorgänge durchgeführt werden oder asymmetrische Öffnungen und transiente Fehler auftreten, erzeugt die Wechselwirkung zwischen der verteilten Kapazität der Leitung und der nichtlinearen Induktivität der Querinduktivität einen hochkomplexen Resonanzkreis.

Aus Sicht der Theorie der elektromagnetischen Felder und gekoppelter Stromkreise ist die Induktivität der Drosselspule kein zeitlich invariantes Parameter mehr, sondern wird zu einer Funktion, die vom verketteten magnetischen Fluss \lambda und dem Strom $i$ abhängt. Das Verhalten des Eisenkerns wird durch die folgende konstitutive Gleichung bestimmt:

i(t) = f(\lambda(t) = \sum_{n=1, 3, 5}^{\infty} c_n \lambda(t)^n

wobei die Koeffizienten $c_n$ die starke Nichtlinearität aufgrund der magnetischen Sättigung widerspiegeln. Während einer einphasigen Öffnung induzieren die gesunden Pole der Übertragungsleitung weiterhin elektrostatische und elektromagnetische Spannungen über die gegenseitigen Zwischenphasenkapazitäten $Cm$ und die Erdkapazitäten $Cs$ auf die getrennte Phase. Wenn die Induktivität der Querinduktivität $L$ eine Resonanzschleife mit der äquivalenten Kapazität der offenen Phase $Ceq$ bildet, sind die Bedingungen für Serien- oder Parallelferroresonanz vollständig gegeben, was anhaltende niederfrequente Überspannungen, sub-synchron geordnete Oberschwingungen und Erregerströme auslöst, die den Nennstrom des Betriebsmittels um mehrere Größenordnungen übersteigen.

Differentialgleichungen und Nichtlineares Modellieren des Sättigbaren Kerns

Die rigorose mathematische Modellierung der Ferroresonanz in Querinduktivitäten erfordert die Formulierung von Systemen gewöhnlicher und partieller Differentialgleichungen, welche die elektrostatischen Phänomene des Netzes mit der Hysterese und der Sättigung des magnetischen Kerns koppeln. Unter Berücksichtigung eines vereinfachten einphasigen Ersatzschaltbildes, das eine getrennte Phase einer Übertragungsleitung mit einer an den Klemmen angeschlossenen Querinduktivität darstellt, wird die Dynamik des Systems durch die folgende Integrodifferentialgleichung beschrieben:

\frac{d^2 \lambda}{dt^2} + \frac{R}{L_0} \frac{d\lambda}{dt} + \omega_0^2 f(\lambda) = \omega_0^2 V_m \sin(\omega t + \theta)

Wobei \lambda der im Reaktor verkettete Fluss ist, $R$ der äquivalente Verlustwiderstand im Kupfer und im Kern, $L_0$ die Induktivität im linearen Bereich, \omega_0 die Fundamentalfrequenz des Systems und V_m \sin(\omega t + \theta) die über die kapazitive Kopplung von den benachbarten, unter Spannung stehenden Phasen induzierte Spannung darstellt. Die Nichtlinearitätsfunktion des magnetischen Kerns f(\lambda) wird üblicherweise durch die modifizierte Fröhlich-Approximation oder durch analytische Sättigungskurven auf Basis hyperbolischer Arkustangensfunktionen modelliert:

i(\lambda) = I_0 \left[ \frac{\lambda}{\lambda_s} + M \sinh\left( N \left(\frac{\lambda}{\lambda_s}\right)^k \right) \right]

In diesem System führt das Vorhandensein der Sättigung dazu, dass sich die scheinbare Induktivität der Drosselspule während jedes Zyklus der angelegten Spannung dynamisch verändert. Wenn der Kern in den Bereich der tiefen Sättigung eintritt, sinkt die differentielle Induktivität der Drosselspule Ldiff = d\lambda / di drastisch ab und nähert sich den Werten der Streuinduktivität der Wicklung. Diese plötzliche Reduzierung verschiebt die Resonanzfrequenz des Kreises auf Werte nahe oder unterhalb der Fundamentalfrequenz, wodurch subharmonische Schwingungsmoden (typischerweise der Ordnung 1/2, 1/3 oder 1/5) und quasi-periodische Schwingungen induziert werden.

Die Lyapunov-Stabilitätsanalyse und Bifurkationsdiagramme zeigen, dass das dynamische System der ferroresonanten Drosselspule mehrere stabile Zustände (Multistabilität) aufweist. Abhängig von den Anfangsbedingungen (wie dem Ausschaltwinkel des Leistungsschalters und dem Momentanwert des remanenten Flusses \lambda_r im Kern) kann das System zu einem normalen Niederspannungszustand konvergieren oder plötzlich in einen ferroresonanten Zustand springen, der durch verzerrte Ströme und hohe magnetische Flüsse gekennzeichnet ist, welche den Knickpunkt der Sättigungskurve um mehr als 200% überschreiten.

Topologie des Phänomens bei Einphasigen Schaltvorgängen

Einpolige Ein- und Ausschaltvorgänge in Hochspannungsschaltern sind Standardbetriebsmaßnahmen zur Beseitigung temporärer Fehler auf Übertragungsleitungen. Während des einphasigen Öffnungsintervalls bleibt die betroffene Phase jedoch von der Hauptstromquelle getrennt, ist jedoch weiterhin elektromagnetisch mit dem Rest des unter Spannung stehenden Systems gekoppelt. Die gegenseitigen Kapazitäten zwischen den Phasen $Cab, Cbc, Cca$ und die Erdkapazitäten $Can, Cbn, Ccn$ bilden ein komplexes kapazitives Netzwerk, das die offene Phase mit einer Spannung elektrostatischen Ursprungs speist.

Wenn eine Querinduktivität permanent an den Klemmen der Leitung angeschlossen ist, wird die Induktivität dieser Drossel parallel zur äquivalenten Null- und Gegen-Systemkapazität geschaltet, die von den offenen Klemmen der Phase aus gesehen wird. Das aus Sicht der Drosselklemmen gesehene Thevenin-Ersatzschaltbild besteht aus einer äquivalenten Spannungsquelle $Vth$ , die von den Spannungen der gesunden Phasen über die Gegenkopplungskondensatoren angetrieben wird, und einer äquivalenten Impedanz $Zth$ , die von der äquivalenten Kapazität $Ceq$ dominiert wird.

Schaltparameter	Normal- / Betriebsbedingung	Kritische Ferroresonanzbedingung	Elektromechanische Konsequenz
Schalterzustand	Symmetrisches dreiphasiges Schließen	Asymmetrisches einphasiges Öffnen	Schwere Spannungsasymmetrie und Mitsystemströme.
Äquivalente Kapazität ( $Ceq$ )	Geringe kapazitive Leitungsimpedanz	Optimale Kopplung mit $Lreactor$	Abstimmung des Resonanzkreises auf Subharmonische.
Remanenter Fluss (\lambda_r)	Gering oder symmetrisch verteilt	Hoher remanenter Fluss nach der Unterbrechung	Asymmetrie in der Hystereseschleife und unmittelbarer Eintritt in die Sättigung.
Systemverluste ( $Req$ )	Normale Lastdämpfung	Geringe Dämpfung im Leerlauf / bei Schwachlast	Unfähigkeit des Systems, die angesammelte Schwingungsenergie abzuführen.

Während des Öffnens eines Phasenpols kann der transiente Lichtbogen im Schalter hochfrequente Wiederzündungen (Re-strikes) erzeugen, die transiente Hochfrequenzkomponenten einspeisen, welche die Eigenmoden der Drossel-Leitungs-Kombination anregen. Wenn der remanente Fluss in der dem offenen Zweig entsprechenden Kernsäule hoch ist und in der Polarität mit der anfänglichen Halbperiode der induzierten Spannung übereinstimmt, erreicht der Kern im ersten Halbzyklus nach dem Erlöschen des Lichtbogens eine tiefe magnetische Sättigung, wodurch der Weg zur ferroresonanten Instabilität eingeleitet wird.

Forensische Fehleranalyse und Auswirkungen auf Unterstationsanlagen

Ferroresonanzereignisse und extreme magnetische Sättigung in Querinduktivitäten erzeugen schwere thermische, dielektrische und mechanische Belastungen, welche die Integrität der Unterstationsanlagen gefährden. Die forensische Analyse von Fehlern in Höchstspannungssystemen zeigt charakteristische Schadensmuster in den verschiedenen Komponenten des elektrischen Systems.

Mechanische Verformung und Elektrodynamische Kräfte in den Wicklungen

Die ferroresonanten Erregerströme bestehen massiv aus Oberschwingungen ungerader Ordnung (hauptsächlich der dritte und fünfte Oberschwingungsanteil) und subharmonischen Komponenten. Diese Ströme können das 3- bis 5-fache des konstruktiven Nennstroms der Drossel erreichen. Die elektromechanischen Kräfte auf die Wicklungsleiter sind proportional zum Quadrat des Momentanstroms (F \propto i(t)^2). Die zyklischen radialen und axialen Kräfte induzieren mechanische Ermüdung in den Isolierhalterungen, was zur Lockerung der Wicklungsdruckblöcke, zur Abriebbeschädigung der Kraftpapierisolierung und schließlich zu Windungsschlüssen führt.

Zusätzliche Erwärmung durch Wirbelströme und Kernverluste

Der Betrieb des magnetischen Kerns im Bereich der tiefen Sättigung verändert die Verteilung des magnetischen Flusses und zwingt den Fluss dazu, aus den kornorientierten Siliziumstahlblechen in die Haltestrukturen, magnetischen Abschirmungen und den Metalltank des Reaktors auszutreten. Die Wirbelstromverluste in den Abschirmplatten und Tankwänden steigen nichtlinear an:

Peddy = \frac{k_f \cdot d^2 \cdot f^2 \cdot Bmax^2}{\rho}

Wobei $d$ die Lamellendicke, $f$ die effektive Frequenz (einschließlich Oberschwingungen) und $Bmax$ die magnetische Spitzenflussdichte ist. Dieser massive Anstieg der Verluste erzeugt lokale Hotspots im Isolieröl, die die zulässigen thermischen Grenzwerte nach den Normen IEEE C57.109 und IEC 60076-6 überschreiten, was zum thermischen Abbau der Zelluloseisolation, zur Entstehung von im Öl gelösten Gasen (Wasserstoff, Methan, Acetylen) und zum Auslösen der Buchholzschutzeinrichtungen führt.

Dielektrische Überspannungen und Zusammenbruch der Isolation

Anhaltende ferroresonante Spannungen können Werte von bis zu 2,5 bis 3,5 p.u. erreichen, mit stark verzerrten Wellenformen, die abrupte Wellenfronten aufweisen. Diese Überspannungen wirken kontinuierlich auf das Hauptisolationssystem (Öl-Papier oder feste Trockenisolation). Die an den Drosselklemmen angeordneten Metalloxid-Ausklinkableiter (MOA) sind einer übermäßigen thermischen Energieabsorption ausgesetzt, wenn sie versuchen, diese niederfrequenten Überspannungen zu begrenzen, was zum thermischen Versagen und zur katastrophalen Explosion der Ableiter führen kann, wenn ihre Energieableitfähigkeit (Leitungsklasse) überschritten wird.

Fortgeschrittene Minderungswegsstrategien und Entwurfskriterien

Die wirksame Minderung von Ferroresonanz und magnetischer Sättigung in Querinduktivitäten erfordert die Implementierung von Gegenmaßnahmen auf Ebene des elektromagnetischen Designs, der Auswahl von Schaltgeräten und der Konfiguration fortschrittlicher Schutzsysteme.

Optimierung des Kerndesigns und der Luftspalte

Um den schnellen Eintritt in die Sättigung unter Bedingungen temporärer Überspannung und kapazitiver Kopplung zu verhindern, werden moderne Drosselspulen mit Panzerkernen oder Kernen mit mehreren nichtmagnetischen Luftspalten konstruiert, die gleichmäßig über die Höhe der Säulen verteilt sind. Dies linearisiert die Magnetisierungskennlinie bis zu 140% - 160% der Nennspannung. Die linearisierte Induktivität wird ausgedrückt durch:

L = \frac{\mu_0 N^2 A}{\sum lcore/\mu_r + \sum lgap} \approx \frac{\mu_0 N^2 A}{\sum lgap}

Indem die Gesamtlänge der Luftspalte \sum lgap zum dominanten Term im Nenner gemacht wird, wird der gesamte magnetische Reluktanzwiderstand unabhängig von der relativen Permeabilität des Stahls \mu_r, wodurch die Nichtlinearität unterdrückt und die Ferroresonanz verhindert wird.

Implementierung von Einschaltwiderständen und Schaltersynchronisation

Der Einsatz von Leistungsschaltern, die mit phasenwinkelgesteuerten Mechanismen (Point-on-Wave Switching) ausgestattet sind, ermöglicht die Eliminierung schwerer Ein- und Ausschalttransienten. Durch das Öffnen oder Schließen der Schalterpole bei Stromnulldurchgängen oder optimalen Spannungsmomenten wird die Erzeugung asymmetrischer remanenter Flüsse minimiert. Zusätzlich dämpft die Verwendung von Einschaltwiderständen in den Höchstspannungsschaltern transiente Schwingungen hoher und mittlerer Frequenz während der Schaltvorgänge.

Dämpfungsschemen und Passive Filter

In Unterstationen, in denen das Ferroresonanzrisiko aufgrund komplexer Netztopologien mit langen Leerlaufleitungen strukturell hoch ist, werden Dämpfungsnetzwerke implementiert, die auf Hochwiderstandskreisen basieren, welche in Reihe mit Hilfskondensatorbänken geschaltet sind, oder durch das Hinzufügen von Tertiärwicklungen in den Reaktoren, die in geschlossener Dreieckschaltung mit einem externen Dämpfungswiderstand verbunden sind. Dieser Widerstand disspiert die Schwingungsenergie der Subharmonischen und verhindert das Anwachsen der Amplitude der ferroresonanten Schwingung.

Praktische Anwendung und Analyse mit der Vexten Suite

Zur Veranschaulichung der Dimensionierung und technischen Validierung eines Blindleistungskompensationssystems, das transienten und harmonischen Regimen ausgesetzt ist, wird eine Fallstudie vorgestellt, die mit den analytischen Modulen der Vexten Suite unter Anwendung der internationalen Normen IEC 60909, IEEE 141 und IEC 60287 entwickelt wurde.

Spezifikationen des Leistungssystems

Systemnennspannung ( $U_n$ ): 400 kV (Drehstromnetz, 50 Hz).
Dreiphasige Kurzschlussleistung des Netzes ( $Sk3$ ): 25 GVA.
Dreiphasige Querinduktivität: 150 MVAR, 400/√3 kV pro Phase, Sternschaltung mit steif geerdetem Sternpunkt.
Äquivalente Nullsystemkapazität der Leitung ( $C_0$ ): 12,5 µF pro Phase.
Äquivalenter Wicklungswiderstand des Reaktors ( $Rcu$ ): 0,15 \Omega pro Phase bei 75 °C.

Kurzschlussberechnung und Transiente Parameter (IEC 60909 / IEEE 141)

Die von den Drosselklemmen aus gesehene äquivalente Netzimpedanz wird aus der Kurzschlussleistung berechnet:

Znet = \frac{U_n^2}{Sk3} = \frac{(400 \times 10^3)^2}{25 \times 10^9} = 6.4 \ \Omega

Die äquivalente Induktive Netzreaktanz ( $Xnet$ ) unter Annahme eines $X/R$ -Verhältnisses von 20:

Rnet = \frac{Znet}{\sqrt{1 + (X/R)^2}} = \frac{6.4}{\sqrt{1 + 400}} = 0.319 \ \Omega

Xnet = Rnet \times 20 = 0.319 \times 20 = 6.38 \ \Omega

Die Nenninduktivität der Querinduktivität pro Phase wird über die Nennreaktanz des Betriebsmittels bestimmt:

X_L = \frac{(U_n / \sqrt{3})^2}{Q_3 / 3} = \frac{(400 / \sqrt{3})^2}{50} = \frac{53333.33}{50} = 1066.67 \ \Omega

L = \frac{X_L}{2 \pi f} = \frac{1066.67}{2 \pi (50)} = 3.395 \ H

Bewertung des Resonanzfrequenzfaktors

Während der einphasigen Öffnung wird die Resonanzfrequenz des aus der Induktivität der Drossel $L$ und der äquivalenten Kapazität der offenen Phase Ceq \approx C_0 gebildeten Kreises berechnet durch:

fres = \frac{1}{2 \pi \sqrt{L \cdot Ceq}} = \frac{1}{2 \pi \sqrt{3.395 \times (12.5 \times 10^{-6})}}

fres = \frac{1}{2 \pi \sqrt{4.2437 \times 10^{-5}}} = \frac{1}{2 \pi (6.514 \times 10^{-3})} = 24.45 \ Hz

Der Wert von $fres = 24.45 Hz$ befindet sich genau im subsynchronen Bereich (ungefähr die Hälfte der Fundamentalfrequenz von 50 Hz), was die Anfälligkeit des Systems gegenüber der Anregung von Subharmonischen der Ordnung 1/2 bei asymmetrischen Schaltbedingungen bestätigt.

Oberschwingungsabsenkung und Leiterdimensionierung (IEC 60287 / NEC 310)

Als Folge der magnetischen Sättigung und der ferroresonanten Transienten enthalten die durch Verbindungskabel und Wicklungen fließenden Ströme einen signifikanten Oberschwingungsgehalt. Die Vexten Suite bewertet den Zusatzverlustfaktor und wendet den Oberschwingungs-Minderungsfaktor (Harmonic Derating Factor - HDF) gemäß der Norm IEC 60287 an:

I_{max\_derated} = Irated \cdot \sqrt{\frac{1 + \Delta_1}{1 + \sum_{h=2}^{n} M_h^2 \cdot (1 + Y_h)}}

Wobei $M_h$ die Oberschwingungsverzerrungsrate der Ordnung $h$ ist, $Y_h$ der Wirbelstrom-Verlustfaktor im Leiter für die Oberschwingung $h$ ist und \Delta_1 den Zuwachs der Hysterese- und Wirbelstromverluste im Kern darstellt. Für einen Gehalt an dritter Oberschwingung von 25% und fünfter Oberschwingung von 15%, abgeleitet von einem schweren transienten Ereignis, zeigt der vom Vexten-Suite-Werkzeug berechnete Minderungsfaktor an, dass die zulässige Strombelastbarkeit des Kabels um 18,5% reduziert werden muss, um den thermischen Abbau der polymeren Isolation (XLPE) zu verhindern.

Vexten Suite Modul	Referenznorm	Haupt-Eingabeparameter	Analytisches Ergebnis / Entwurfskriterium
Vexten Transient Studio	IEC 60076-6 / IEEE C57.109	Nichtlineare Induktivitätsmatrix und remanente Flüsse	Identifizierung von Bifurkationspunkten und ferroresonanten Überspannungen von 2,8 p.u.
Vexten Short-Circuit Pro	IEC 60909 / IEEE 141	Kurzschlussleistung ( $Sk3$ ) und Sequenzimpedanzen	Bestimmung elektrodynamischer Spitzenkräfte und asymmetrischer Ströme.
Vexten Thermal Cable Sizer	IEC 60287 / NEC 310	Stromoberschwingungsspektrum und Umgebungstemperatur	Berechnung des HDF-Faktors und thermische Überdimensionierung von Leistungskabeln.
Vexten Protection Coordination	IEC 60255 / IEEE C37.112	Erregungs- und temporäre Überstromkennlinien	Einstellung von zeitgesteuerten Überstromrelais und Restüberspannungsschutz.

Schlussfolgerungen und Normative Richtlinien

Die Instabilität durch Ferroresonanz und magnetische Sättigung in Querinduktivitäten während einphasiger Schaltvorgänge stellt eines der komplexesten und destruktivsten elektromagnetischen Phänomene in elektrischen Hochspannungs-Leistungssystemen dar. Die rigorose Analyse erfordert die Abkehr von traditionellen linearen Modellen und die Annahme nichtlinearer Simulationsansätze, welche die Hysterese, die Kernsättigung und die verteilte kapazitive elektromagnetische Kopplung der Übertragungsleitungen einbeziehen.

Die Implementierung von Kerndesigns mit verteilten Mehrfachluftspalten, der Einsatz phasenwinkelgesteuerter synchronisierter Schalter und die Validierung durch fortgeschrittene Engineering-Tools wie die Vexten Suite garantieren, dass Blindleistungskompensationsanlagen innerhalb der von den internationalen IEC- und IEEE-Normen festgelegten Sicherheitsmargen betrieben werden. Die strikte Einhaltung dieser Kriterien verhindert katastrophale Ausfälle, optimiert die Investitionen in Unterstationsanlagen und sichert die Zuverlässigkeit und Betriebsstabilität von elektrischen Großstrom-Übertragungsnetzen.