Sistemas de PotenciaSSRResonancia TorsionalEstabilidad de RedIngeniería Eléctrica

Subsynchrone Resonanz (SSR) Instabilität bei reihenkompensierten Generatorsystemen

Wir analysieren die Instabilität durch subsynchrone Resonanz (SSR) bei reihenkompensierten Generatoren, ihre mechanischen Auswirkungen und technische Minderungsmaßnahmen.

Ing. Francisco Ramírez

Elektromechanische Grundlagen und Ursprung der subsynchronen Resonanz

Die Instabilität durch subsynchrone Resonanz (SSR) stellt eines der komplexesten und destruktivsten Phänomene in Höchstspannungs-Wechselstrom-Übertragungssystemen (HVAC) dar. Sie tritt typischerweise in synchronen Erzeugungsanlagen auf, die an Weitverkehrs-Übertragungsnetze mit kapazitiver Reihenkompensation angeschlossen sind. Der grundlegende Zweck der Reihenkompensation besteht darin, die induktive Gesamtreattanz der Übertragungsleitung zu reduzieren, wodurch die Wirkübertragungsfähigkeit erhöht, der Phasenwinkel optimiert und die Transientenstabilität des vernetzten Systems verbessert wird. Die Einführung von Reihenkondensatorbänken in Serie mit den Netzinduktivitäten erzeugt jedoch elektrische Reihen-LC-Schwingkreise, deren Eigenfrequenzen unterhalb der industriellen Grundfrequenz des Systems (50 Hz oder 60 Hz) liegen.

Aus Sicht der Theorie verteilter und konzentrierter Schaltkreise wird die elektrische Eigenschwingungsfrequenz eines serienkompensierten Übertragungssystems durch folgenden Grundausdruck definiert:

fer=f0XCXL=f0ksefer = f_0 \sqrt{ \frac{X_C}{X_L} } = f_0 \sqrt{ kse }

Wobei f0f_0 die Grundfrequenz des Systems (50 Hz oder 60 Hz) darstellt, XCX_C die kapazitive Gesamtreattanz der Reihenbänke ist, XLX_L die äquivalente induktive Gesamtreattanz des Systems aus Sicht der Generatorklemmen ist und ksekse den Reihenkompensationsgrad, ausgedrückt in p.u. oder Prozent, bezeichnet. Da die Kompensationsgrade in langen Leitungen üblicherweise zwischen 25 % und 70 % variieren, ist die resonante elektrische Frequenz ferfer streng kleiner als die Grundfrequenz. Folglich regt jede transiente Störung oder Oberschwingung im Netz diese natürliche elektrische Frequenz an, was zu subsynchronen Strömen und Spannungen führt.

Wenn diese subsynchronen Ströme durch die Statorwicklungen des Synchrongenerators fließen, erzeugen sie ein rotierendes Magnetfeld im Luftspalt der Maschine, das mit einer Winkelgeschwindigkeit rotiert, die geringer ist als die Synchrondrehzahl des Rotors. Dieses subsynchrone Drehfeld interagiert mit dem Gleichstrom-Magnetfeld des Rotors und erzeugt oszillierende elektromagnetische Drehmomente auf der Generatorwelle, deren Pulsationsfrequenz fmf_m durch die Differenz zwischen der Grundfrequenz und der resonanten elektrischen Frequenz gegeben ist:

fm=f0fer=f0(1kse)f_m = f_0 - fer = f_0 \left( 1 - \sqrt{kse} \right)

Wenn die Frequenz dieses oszillierenden Drehmoments mit einer der Torsions-Eigenschwingungsfrequenzen des Wellenstrangs der Generatoreinheit übereinstimmt oder sich dieser stark annähert (welcher Hochdruck-, Mitteldruck- und Niederdruckturbine sowie den Generatorrotor und den Erreger umfasst), kommt es zu einer schweren mechanisch-elektrischen Kopplung. Dieser Mechanismus wird formal als Torsionsmoden-Wechselwirkung (Torsional Interaction - TI) bezeichnet. Die Maschine agiert als mechanischer Energieverstärker durch negative oder positive Rückkopplung der Torsionsdämpfung, was innerhalb von Sekunden zu extremer Torsionsermüdung und zum katastrophalen Bruch der Welle führen kann.

Klassifizierung und physikalische Mechanismen der SSR

Das globale Phänomen der subsynchronen Resonanz ist kein monolithisches Ereignis, sondern unterteilt sich gemäß der IEEE SSR Task Force formal in drei unterschiedliche, aber miteinander verknüpfte physikalische Mechanismen:

Rotorinduktionsgenerator-Effekt (Rotor Induction Generator Effect - RIGE)

Dieses Phänomen ist streng elektrischer Natur und kann selbst mit einem ideal starren Rotor (ohne Torsionsfreiheitsgrade) auftreten. Wenn subsynchrone Ströme der Frequenz ferfer durch den Stator fließen, erzeugen sie ein Magnetfeld, das mit einer geringeren Geschwindigkeit als der Rotor rotiert. Aus der Sicht des Referenzsystems des mit der Synchrondrehzahl \omega_0 rotierenden Rotors erscheint dieses Feld mit einer subsynchronen Schlupffrequenz fslip=f0ferfslip = f_0 - fer. Diese relative Frequenz induziert Ströme in den Dämpferkreisen (Dämpferwicklungen) und im massiven Körper des Vollstahlrotors.

Die vom Stator aus bei der subsynchronen Frequenz gesehene Äquivalenzimpedanz beinhaltet den reflektierten effektiven Rotorwiderstand. Aufgrund des negativen Schlupfs für Frequenzen unterhalb der Grundfrequenz (fer<f0fer < f_0) kann der vom Stator bei dieser subsynchronen Frequenz gesehene scheinbare Rotorwiderstand negativ werden. Ein negativer Gesamtwiderstand von Stator und Rotor bei der Frequenz ferfer bedeutet, dass der Stromkreis Energie aus dem Netz aufnimmt und verstärkt oder dass der Stromkreis die elektrischen Schwingungen negativ dämpft, was zu einer elektrischen Selbsterregung führt, wenn der Netzwiderstand Rtotal(fer) < 0 ist.

Torsionsmoden-Wechselwirkung (Torsional Interaction - TI)

Die Torsionswechselwirkung tritt auf, wenn das elektrische System und der mechanische Wellenstrang über das Magnetfeld des Luftspalts gekoppelt werden. Angenommen, der Wellenstrang besitzt eine mechanische Torsions-Eigenschwingungsfrequenz fmechfmech. Wenn eine Torsionsschwingung der Welle mit dieser Frequenz auftritt, modulieren die Winkeländerungen des Rotors die Statorströme, wodurch Seitenbänder im elektrischen Spektrum bei den Frequenzen f_0 \pm fmech erzeugt werden. Wenn eine dieser Seitenbandfrequenzen mit der komplementären Eigenfrequenz des elektrischen Reihenstromkreises übereinstimmt (fer=f0fmechfer = f_0 - fmech), wird ein geschlossener Rückkopplungskreis aufgebaut.

Mathematisch wird diese Kopplung durch die Berechnung des Torsions-Dämpfungskoeffizienten DeD_e für jeden Torsionsmodus ii bewertet. Das logarithmische Decrement oder die Gesamtdämpfung DtotalDtotal eines Torsionsmodus ist die Summe aus inhärenter mechanischer Dämpfung DmD_m und induzierter elektrischer Dämpfung DeD_e:

Dtotal,i=Dm,i+De,iD_{total, i} = D_{m, i} + D_{e, i}

Wenn der Netzwiderstand und die Eigenschaften der LC-Kreise dazu führen, dass De,iD_{e, i} negativ ist und dessen Betrag die positive mechanische Dämpfung übersteigt (De,i>Dm,i|D_{e, i}| > D_{m, i}), wird die Gesamtdämpfung des Systems negativ (Dtotal,i<0D_{total, i} < 0). Dies führt zu einer exponentiellen Instabilität, bei der die Torsionsschwingungen unkontrolliert anwachsen, bis es zur mechanischen Zerstörung der Wellenkupplungen kommt.

Drehmomenttransienten-Verstärkung (Torque Amplification - TA)

Im Gegensatz zu RIGE und TI, bei denen es sich um dynamische Kleinzelldynamik- und Stabilitätsphänomene handelt, ist die Drehmomenttransienten-Verstärkung ein transientes Großsignalphänomen. Sie tritt auf, wenn es in der Nähe von mit Reihenkondensatoren kompensationierten Leitungen zu schweren Kurzschlüssen kommt, gefolgt von Öffnungs- und Wiedereinschaltvorgängen von Leistungsschaltern. Das abrupte Entladen und Wiederaufladen der Reihenkondensatorbänke interagiert mit transienten Oberschwingungen und Fehlerströmen und prägt der Generatorwelle einen starken elektromagnetischen Drehmomentimpuls auf.

Wenn die Frequenz dieses Impulses mit den Eigenmoden der Welle übereinstimmt, erfährt der Wellenstrang mechanische Beanspruchungen, die das Nenndrehmoment bei weitem überschreiten können, wodurch in einem einzigen transienten Ereignis schwerer plastischer Ermüdungsschaden akkumuliert wird.

Mathematische Modellierung und Kleinsignal-Analyse in Multimodesystemen

Zur Durchführung einer strengen Analyse der Stabilität durch subsynchrone Resonanz in fortgeschrittenen Industrie- und Übertragungsumgebungen wird die linearisierte Zustandsraummodellierung um einen Arbeitspunkt im stationären Betrieb verwendet. Das Gesamtsystem wird in drei gekoppelte Teilsysteme zerlegt: das elektrische Übertragungsnetz mit Reihenkompensation, den Synchrongenerator mit seinen Park-Gleichungen (d-q-Achsen) und das mechanische Mehr-Massen-Modell des Wellenstrangs.

Das mechanische Modell des Wellenstrangs mit NN rotierenden Massen wird durch ein System gewöhnlicher Differentialgleichungen zweiter Ordnung dargestellt:

Jid2θidt2+Dmidθidt+j=1NKij(θiθj)=TmiTei,i=1,2,,NJ_i \frac{d^2 \theta_i}{dt^2} + Dmi \frac{d \theta_i}{dt} + \sum_{j=1}^{N} Kij (\theta_i - \theta_j) = Tmi - Tei, \quad i = 1, 2, \dots, N

Wobei JiJ_i das Trägheitsmoment der Masse ii ist, \langle\theta_i der rotatorische Winkel der Masse, DmiDmi der Eigendämpfungskoeffizient, KijKij die Torsionssteifigkeit der Welle zwischen den Massen ii und jj, TmiTmi das von der Turbine aufgebrachte mechanische Drehmoment und TeiTei das auf die Generatormasse induzierte elektromagnetische Drehmoment ist.

Das elektrische System umfasst die Park-Transformationsgleichungen für Stator und Rotor, gekoppelt mit den Differentialgleichungen der R-L-C-Reihenübertragungsstromkreise:

diddt=1Ld(Raid+ωψqvd)\frac{d i_d}{dt} = \frac{1}{L_d} \left( -R_a i_d + \omega \psi_q - v_d \right)
diqdt=1Lq(Raiqωψdvq)\frac{d i_q}{dt} = \frac{1}{L_q} \left( -R_a i_q - \omega \psi_d - v_q \right)

Die Gleichungen der kompensierten Übertragungsleitung werden in netzsynchronisierten Achsenkoordinaten durch Differentialgleichungen ausgedrückt, die Lastdynamik und Reihenkondensatoren einbeziehen:

dvcddt=ω0id+1Cseid\frac{d vcd}{dt} = \omega_0 i_d + \frac{1}{Cse} i_d
dvcqdt=ω0iq+1Cseiq\frac{d vcq}{dt} = -\omega_0 i_q + \frac{1}{Cse} i_q

Durch Linearisierung dieses nichtlinearen Gleichungssystems um einen stabilen Arbeitspunkt (X0,U0X_0, U_0) erhält man die Zustandsmatrix des geschlossenen Systems:

ΔX˙=AΔX+BΔU\Delta \dot{X} = A \Delta X + B \Delta U
ΔY=CΔX+DΔU\Delta Y = C \Delta X + D \Delta U

Die Kleinsignal-Stabilitätsanalyse wird durch Berechnung der Eigenwerte (Eigenvalues) \lambda_i = \sigma_i \pm j \omega_i der Zustandsmatrix AA durchgeführt. Der Realteil \sigma_i jedes Eigenwerts bestimmt die Stabilität des entsprechenden Modus:

  • Wenn \sigma_i < 0, ist der Modus asymptotisch stabil (positive Dämpfung).
  • Wenn \sigma_i > 0, ist der Modus instabil (negative Dämpfung, exponentielles Anwachsen der Schwingung).
  • Wenn \sigma_i = 0, befindet sich der Modus an der Grenze der marginalen Stabilität.

Für die elektromechanischen Torsionsmoden steht der Realteil des Eigenwerts in direktem Zusammenhang mit der Gesamtdämpfung: \sigma_i = -\frac{D_{total, i}}{2 J_i}. Daher zeigt ein Eigenwert mit positivem Realteil direkt das Auftreten einer Instabilität durch Torsionswechselwirkung (TI) an.

Vergleichsmatrix internationaler Normen und Bewertungskriterien

Die Bewertung und Minderung der SSR unterliegen strengen internationalen Standards, die vom IEEE und der IEC herausgegeben werden. Die folgende Vergleichsmatrix detailliert kritische Parameter, normative Grenzen und operative Konsequenzen im Zusammenhang mit SSR in Hochspannungsnetzen.

Norm / Standard Bewertungsparameter Kritischer Grenzwert / Akzeptanzkriterium Zugehörige Fehlerbedingung Operative und dielektrische Konsequenz
IEEE Std 1133 Elektrische Torsionsdämpfung (DeD_e) Dtotal \ge 0 unter allen Lastflussbedingungen (N0,N1N-0, N-1) Instabilität durch Torsionswechselwirkung (TI) Akkumulierte Torsionsermüdung, Verformung und katastrophaler Wellenbruch.
IEEE Std 693 Gekoppelte seismische und Torsionsantwort Maximale mechanische Beanspruchung unter 80 % der Elastizitätsgrenze Verstärkung durch kombinierte Transienten-Drehmomente Strukturelle Schäden an Turbogeneratorträgern und Unterstationsfundamenten.
IEC 60383 / IEC 60871 Überspannungen in Reihenkondensatorbänken Vpeak \le 1.8 p.u. im stationären subsynchronen Betrieb Anhaltende subsynchrone Resonanz und Ferroresonanz Dielektrischer Durchschlag der Kondensatorelemente, Auslösung der VMS (Metalloxid-Varistoren).
IEEE Std 1547 / C50.13 Subsynchrone Ströme in Statorwicklungen Isub \le 0.05 p.u. dauerhaft; Isub \le 0.15 p.u. für 10 s Rotorinduktions-Effekt (RIGE) und lokalisierte Erwärmung Kritische Überhitzung von Rotorkeilen, Halteringen und Statorkern.

Forensische Fehleranalyse an kritischen Betriebsmitteln bei SSR

Das Vorhandensein von subsynchronen Strömen und Drehmomenten erzeugt schwere Beanspruchungen, die über den Synchrongenerator hinausgehen und die gesamte zugehörige Unterstation sowie die Transformierungs- und Schaltinfrastruktur tiefgreifend beeinträchtigen. Die forensische Fehleranalyse offenbart die Degradationsmechanismen in folgenden Komponenten:

Leistungstransformatoren und Autotransformatoren

Die der Grundfrequenz überlagerten subsynchronen Stromkomponenten (ferfer) erzeugen temporäre Asymmetrien im magnetischen Fluss des Transformatorenkerns. Da der magnetische Fluss das zeitliche Integral der Spannung ist, erhöht das Vorhandensein niederfrequenter magnetischer Flusskomponenten die Spitzenflussdichte \hat{B} und treibt den Kern tief in den Bereich der nichtlinearen magnetischen Sättigung.

Die periodische Kernsättigung verursacht:

  • Massive Erzeugung von Oberschwingungen höherer Ordnung (insbesondere 2., 3. und 5. Harmonische), welche die Spannungswelle verzerren.
  • Drastischen Anstieg der Erregerströme sowie der Eisenverluste durch Hystereseverluste und Wirbelströme (Foucaultsche Ströme).
  • Schwere mechanische Schwingungen in den Transformatortanks und -wicklungen infolge nichtlinearer Magnetostriktion und durch asymmetrische Ströme gewichteter Lorentzkräfte, was zur Lockerung mechanischer Verschraubungen und zum vorzeitigen Versagen der Feststoffisolation (ölgetränktes Papier) führt.

Leistungsschalter und Schaltanlagen

Leistungsschalter in mit Reihenkondensatoren kompensierten Leitungen erfahren während des Auftretens von SSR extrem schwere Spannungsverhalten. Die transienten Wiederkehrspannungen (TRV) werden durch die Eigenfrequenz des Reihenkreises ferfer moduliert und weisen sehr hohe Anstiegsgeschwindigkeiten der transienten Spannung (du/dtdu/dt) sowie Spannungsspitzenwerte auf, welche die nominale Ausschaltkapazität des Schalters überschreiten können.

Zusätzlich erhöht sich beim Schalten zur Fehlerklärung bei Vorhandensein subsynchroner Ströme die Wahrscheinlichkeit einer Rückzündung (Re-Strike) oder einer Wiederzündung des Lichtbogens erheblich, da der Strom den Nulldurchgang mit abnormalen Geschwindigkeiten durchläuft oder überlagerte bidirektionale niederfrequente Komponenten enthält. Dies kann zur katastrophalen Zerstörung der SF6-Löschkammern und zum explosionsartigen Versagen des Schalterpols führen.

Hochspannungskabel und Messsysteme

An die Unterstation angeschlossene Hochspannungskabel erleiden einen Anstieg der dielektrischen Verluste sowie der Verluste durch induzierte Ströme in den Metallschirmen (Kupfer- oder Bleischirmung). Die subsynchronen Ströme induzieren niederfrequente Wechselfelder, die die Schirmungen durchdringen und hohe Kreisströme erzeugen, welche die Polymerisolation (XLPE) überhitzen. Dies beschleunigt die Bildung von elektrischen Bäumen (Water/Electrical Treeing) und reduziert die geschätzte Lebensdauer des Kabels drastisch.

Andererseits leiden Stromwandler (CT) und Spannungswandler (PT) unter signifikanten Übersetzungs- und Phasenwinkel-Messfehlern infolge der durch subsynchrone Stromkomponenten induzierten Magnetsättigung, wodurch die Leistung von auf Grundfrequenz-Zeigern basierenden numerischen Schutzrelais beeinträchtigt oder unbrauchbar gemacht wird.

Fortgeschrittene Milderungsstrategien und Engineering-Design

Zur Gewährleistung der betrieblichen Machbarkeit von Übertragungssystemen mit massiver Reihenkompensation implementiert die moderne Ingenieurtechnik eine Reihe technologischer Gegenmaßnahmen, unterteilt in netzseitige Milderung, generatorseitige Milderung und den Einsatz von Leistungselektronik hoher Geschwindigkeit.

Subsynchrone Resonanz-Sperrfilter (Subsynchronous Resonance Blocking Filters - SSBF)

SSBF-Filter sind in Serie mit den Statorwicklungen des Generators oder direkt in der Übertragungsleitung installierte Abstimmkreise. Sie bestehen aus einer Anordnung von Parallelreaktoren und -kondensatoren, die so abgestimmt sind, dass sie der kritischen resonanten elektrischen Frequenz ferfer eine extrem hohe Impedanz (einen virtuellen Leerlauf) entgegensetzen.

Die Impedanz des SSBF-Filters bei Abstimmfrequenz wird wie folgt ausgedrückt:

ZSSBF(j\omegaer) = \frac{R + j\omegaer L}{1 - \omegaer^2 L C + j\omegaer R C} \approx \infty