Instabilité par Résonance Sous-Synchrone (SSR) dans les systèmes de générateurs compensés en série

Fondaments Électromécaniques et Origine de la Résonance Subsynchrone

L'instabilité par résonance subsynchrone (SSR, pour Subsynchronous Resonance) représente l'un des phénomènes les plus complexes et destructifs dans les réseaux électriques de puissance à courant alternatif à très haute tension (THV). Elle se manifeste typiquement dans les centrales de production synchrone connectées à des réseaux de transport ultra-longue distance qui emploient une compensation série capacitive. L'objectif fondamental de la compensation série est de réduire la réactance inductive totale de la ligne de transport, améliorant ainsi la capacité de transit de puissance active, optimisant l'angle de phase et renforçant la stabilité transitoire du réseau interconnecté. Néanmoins, l'insertion de batteries de condensateurs en série avec les inductances du réseau génère des circuits LC résonants en série dont les fréquences propres se situent en deçà de la fréquence fondamentale du réseau industriel (50 Hz ou 60 Hz).

Du point de vue de la théorie des circuits distribués et localisés, la fréquence propre d'oscillation électrique d'un réseau de transport compensé en série se définit par l'expression fondamentale suivante :

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

Où $f_0$ représente la fréquence fondamentale du réseau (50 Hz ou 60 Hz), $X_C$ est la réactance capacitive totale des batteries en série, $X_L$ est la réactance inductive totale équivalente du système vue depuis les bornes du générateur, et $kse$ $ est le taux de compensation série exprimé par unité ou en pourcentage. Puisque les taux de compensation sur les longues lignes oscillent couramment entre 25 % et 70 %, la fréquence électrique résonnante $fer$ s'avère être strictement inférieure à la fréquence fondamentale. Par conséquent, toute perturbation transitoire ou harmonique sur le réseau excite cette fréquence électrique propre, générant des courants et des tensions subsynchrones.

Lorsque ces courants subsynchrones circulent dans les enroulements statoriques du générateur synchrone, ils créent un champ magnétique tournant dans l'entrefer de la machine qui tourne à une vitesse angulaire inférieure à la vitesse synchrone du rotor. Ce champ tournant subsynchrone interagit avec le champ magnétique à courant continu du rotor, produisant des couples électromagnétiques oscillatoires sur l'arbre du générateur dont la fréquence de pulsation $f_m$ est dictée par la différence entre la fréquence fondamentale et la fréquence électrique résonnante :

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

Si la fréquence de ce couple oscillatoire coïncide ou s'approche de manière critique de l'une des fréquences propres de vibration the torsion de la ligne d'arbres du groupe de production (comprenant les turbines à haute pression, moyenne pression, basse pression, le rotor du générateur et l'excitatrice), un couplage mécanico-électrique sévère se produit. Ce mécanisme est formellement désigné sous le terme d'Interaction Torsionnelle de Modes (Torsional Interaction - TI). La machine se comporte alors tel un amplificateur d'énergie mécanique par le biais de la contre-réaction ou de la réaction positive de l'amortissement torsionnel, pouvant conduire à une fatigue torsionnelle extrême et à la rupture catastrophique de l'arbre en l'espace de quelques secondes.

Classification et Mécanismes Physiques de la SSR

Le phénomène global de la résonance subsynchrone ne constitue point un événement monolithique, mais se subdivise formellement en trois mécanismes physiques distincts mais interdépendants selon l'IEEE SSR Task Force :

Effet d'Induction du Rotor (Rotor Induction Generator Effect - RIGE)

Ce phénomène revêt une nature strictement électrique et peut se manifester même en présence d'un rotor idéalement rigide (dépourvu de degrés de liberté torsionnels). Lorsque des courants subsynchrones de fréquence $fer$ circulent dans le stator, ils établissent un champ magnétique tournant à une vitesse inférieure à celle du rotor. Vu depuis le repère du rotor tournant à la vitesse synchrone \omega_0, ce champ apparaît à une fréquence de glissement subsynchrone $fslip = f_0 - fer$ . Cette fréquence relative induit des courants dans les circuits amortisseurs (enroulements damper) ainsi que dans la masse massive du rotor en acier forgé.

L'impédance équivalente vue depuis le stator à la fréquence subsynchrone intègre la résistance effective du rotor ramenée au stator. En raison du glissement négatif inhérent aux fréquences inférieures à la fondamentale ( $fer < f_0$ ), la résistance apparente du rotor perçue par le stator à cette fréquence subsynchrone est susceptible de devenir négative. Une résistance globale du stator et du rotor négative à la fréquence $fer$ signifie que le circuit absorbe de l'énergie du réseau pour l'amplifier, ou que le circuit amortit négativement les oscillations électriques, provoquant ainsi une auto-excitation électrique si la résistance nette du circuit Rtotal(fer) < 0.

Interaction Torsionnelle de Modes (Torsional Interaction - TI)

L'interaction torsionnelle survient lorsque le réseau électrique et la ligne d'arbres mécanique se couplent à travers le champ magnétique de l'entrefer. Supposons que la ligne d'arbres possède une fréquence propre torsionnelle $fmech$ . Si une oscillation torsionnelle se manifeste sur l'arbre à ladite fréquence, les variations angulaires du rotor modulent les courants statoriques, générant des bandes latérales dans le spectre électrique aux fréquences f_0 \pm fmech. Si l'une de ces fréquences latérales coïncide avec la fréquence propre complémentaire du circuit électrique série ( $fer = f_0 - fmech$ ), une boucle de rétroaction fermée s'établit.

Mathématiquement, ce couplage s'évalue par le calcul du coefficient d'amortissement électrique torsionnel $D_e$ pour chaque mode torsionnel $i$ . Le décrément logarithmique ou amortissement total $Dtotal$ d'un mode torsionnel correspond à la somme de l'amortissement mécanique intrinsèque $D_m$ et de l'amortissement électrique induit $D_e$ :

D_{total, i} = D_{m, i} + D_{e, i}

Si la résistance du réseau et les caractéristiques des circuits LC entraînent une valeur de $D_{e, i}$ négative dont l'amplitude excède l'amortissement mécanique positif ( $|D_{e, i}| > D_{m, i}$ ), l'amortissement global du système devient négatif ( $D_{total, i} < 0$ ). Cette situation mène à une instabilité exponentielle où les oscillations torsionnelles s'accroissent de façon incontrôlée jusqu'à la destruction mécanique des accouplements de l'arbre.

Amplification par Transitoires de Couple (Torque Amplification - TA)

À l'instar du RIGE et de la TI, qui relèvent de la dynamique de instabilités de petit signal, l'amplification par transitoires de couple constitue un phénomène transitoire de forte amplitude. Elle se produit lors de courts-circuits sévères à proximité de lignes compensées par des condensateurs de série, suivis d'opérations d'ouverture et de réenclenchement d'interrupteurs de puissance. La décharge et la recharge brusques des batteries de condensateurs série interagissent avec les harmoniques transitoires et les courants de défaut, injectant une impulsion de couple électromagnétique sévère dans l'arbre du générateur.

Si la fréquence de cette impulsion coïncide avec les modes propres de l'arbre, la ligne d'arbres subit des contraintes mécaniques susceptibles de dépasser largement le couple nominal de conception, accumulant ainsi des dommages par fatigue plastique sévère au cours d'un unique événement transitoire.

Modélisation Mathématique et Analyse de Petits Signaux dans les Systèmes Multimodes

Afin de mener une analyse rigoureuse de la stabilité face à la résonance subsynchrone au sein d'environnements industriels et de transport avancés, on recourt à la modélisation par espace d'états linéarisée autour d'un point de fonctionnement en régime permanent. Le système global se décompose en trois sous-systèmes couplés : le réseau de transport électrique à compensation série, le générateur synchrone régi par ses équations de Park (axes $d-q$ ), et le modèle mécanique multi-masses de la ligne d'arbres.

Le modèle mécanique de la ligne d'arbres dotée de $N$ masses tournantes s'exprime par un système d'équations différentielles ordinaires du second ordre :

J_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

Où $J_i$ représente le moment d'inertie de la masse $i$ , \theta_i l'angle de position du rotor, $Dmi$ le coefficient d'amortissement propre, $Kij$ la rigidité torsionnelle de l'arbre entre les masses $i$ et $j$ , $Tmi$ le couple mécanique appliqué par la turbine, et $Tei$ le couple électromagnétique induit sur la masse du générateur.

Le sous-système électrique englobe les équations de transformation de Park pour le stator et le rotor, couplées aux équations différentielles des circuits de transport R-L-C série :

\frac{d i_d}{dt} = \frac{1}{L_d} \left( -R_a i_d + \omega \psi_q - v_d \right)

\frac{d i_q}{dt} = \frac{1}{L_q} \left( -R_a i_q - \omega \psi_d - v_q \right)

Les équations de la ligne de transport compensée sont formulées dans des coordonnées d'axes synchronisés avec le réseau par le biais d'équations différentielles intégrant la dynamique de charge et les condensateurs de série :

\frac{d vcd}{dt} = \omega_0 i_d + \frac{1}{Cse} i_d

\frac{d vcq}{dt} = -\omega_0 i_q + \frac{1}{Cse} i_q

La linéarisation de cet ensemble d'équations non linéaires au voisinage d'un point de fonctionnement stable ( $X_0, U_0$ ) conduit à la matrice d'état du système bouclé :

\Delta \dot{X} = A \Delta X + B \Delta U

\Delta Y = C \Delta X + D \Delta U

L'analyse de la stabilité de petit signal s'effectue par le calcul des valeurs propres (eigenvalues) \lambda_i = \sigma_i \pm j \omega_i de la matrice d'état $A$ . La partie réelle \sigma_i de chaque valeur propre conditionne la stabilité du mode associé :

Si \sigma_i < 0, le mode est asymptotiquement stable (amortissement positif).
Si \sigma_i > 0, le mode est instable (amortissement négatif, croissance exponentielle de l'oscillation).
Si \sigma_i = 0, le mode se trouve à la limite de la stabilité marginale.

Concernant les modes électromécaniques torsionnels, la partie réelle de la valeur propre est directement corrélée au coefficient d'amortissement total : \sigma_i = -\frac{D_{total, i}}{2 J_i}. Par conséquent, une valeur propre dotée d'une partie réelle positive indique de manière univoque l'occurrence d'une instabilité par Interaction Torsionnelle (TI).

Matrice Comparative des Normes Internationales et des Critères d'Évaluation

L'évaluation et la mitigation de la SSR obéissent à des normes internationales strictes édictées par l'IEEE et la CEI. La matrice comparative ci-dessous détaille les paramètres critiques, les limites normatives et les conséquences opérationnelles inhérentes à la SSR dans les réseaux à haute tension.

Norme / Standard	Paramètre d'Évaluation	Limite Critique / Critère d'Acceptation	Condition de Défaut Associée	Conséquence Opérationnelle et Diélectrique
IEEE Std 1133	Amortissement Électrique Torsionnel ( $D_e$ )	Dtotal \ge 0 sous toute condition de transit ( $N-0, N-1$ )	Instabilité par Interaction Torsionnelle (TI)	Fatiga torsionnelle accumulée, déformation et rupture catastrophique des arbres.
IEEE Std 693	Réponse Sismique et Torsionnelle Couplée	Contraintes mécaniques maximales inférieures à 80 % de la limite élastique	Amplification par transitoires de couple combinés	Dommages structurels sur les supports de turbogénérateurs et les fondations de sous-stations.
CEI 60383 / CEI 60871	Surtensions sur les Batteries de Condensateurs Série	Vpeak \le 1.8 p.u. en régime permanent subsynchrone	Résonance subsynchrone prolongée et féro-résonance	Perforation diélectrique des éléments capacitifs, déclenchement des VMS (Varistor Metal-Oxide).
IEEE Std 1547 / C50.13	Courants Subsynchrones dans les Enroulements Statoriques	Isub \le 0.05 p.u. en continu ; Isub \le 0.15 p.u. pendant 10 s	Effet d'induction du rotor (RIGE) et échauffement localisé	Surchauffe critique dans les cales du rotor, les anneaux de retenue et le noyau statorique.

Analyse Forensique des Défauts sur les Équipements Critiques face à la SSR

La présence de courants et de couples subsynchrones engendre des sollicitations sévères qui dépassent le cadre du générateur synchrone pour affecter profondément l'ensemble de l'infrastructure de la sous-station associée ainsi que les systèmes de transformation et de commutation. L'analyse forensique des avaries met en lumière les mécanismes de dégradation sur les composants suivants :

Transformateurs de Puissance et Autotransformateurs

Les composantes de courant subsynchrone ( $fer$ ) superposées à la fréquence fondamentale induisent des asymétries temporelles dans le flux magnétique du noyau du transformateur. Le flux magnétique représentant l'intégrale temporelle de la tension, la présence de composantes de basse fréquence magnétomotrice élève la densité de flux de crête \hat{B}, poussant le noyau profondément au cœur de la zone de saturation magnétique non linéaire.

La saturation périodique du noyau provoque :

Une génération massive d'harmoniques d'ordre supérieur (notamment les 2e, 3e et 5e harmoniques) qui déforment l'onde de tension.
Une augmentation drastique des courants d'excitation ainsi que des pertes dans le fer par hystérésis et par courants de Foucault.
Des vibrations mécaniques sévères dans les cuves et les enroulements du transformateur en raison de la magnétostriction non linéaire et des forces de Lorentz pondérées par les courants asymétriques, ce qui entraîne le desserrage des serrages mécaniques et la défaillance prématurée de l'isolation solide (papier imprégné d'huile).

Disjoncteurs de Puissance et Systèmes de Manœuvre

Les disjoncteurs de puissance situés sur les lignes compensées par des condensateurs en série subissent des conditions de rétablissement de tension extrêmement sévères lors de l'apparition d'une SSR. Les tensions transitoires de rétablissement (TTR) se trouvent modulées par la fréquence propre du circuit série $fer$ , affichant des vitesses de montée de tension transitoire ( $du/dt$ ) très élevées ainsi que des valeurs de crête de tension susceptibles de surpasser le pouvoir de coupure assigné du disjoncteur.

En outre, lorsqu'un disjoncteur intervient pour éliminer un défaut en présence de courants subsynchrones, la probabilité d'un réamorçage (re-strike) ou d'une reconduction de l'arc électrique s'accroît considérablement du fait que le passage par zéro du courant s'effectue à des vitesses anormales ou contient des composantes bidirectionnelles superposées à basse fréquence. Cette situation peut causer la destruction catastrophique des chambres d'extinction d'arc à SF6 et la défaillance explosive du pôle du disjoncteur.

Câbles à Haute Tension et Systèmes de Mesure

Les câbles isolés à haute tension raccordés à la sous-station subissent un accroissement des pertes diélectriques ainsi que des pertes par courants induits dans les écrans métalliques (blindage en cuivre ou en plomb). Les courants subsynchrones induisent des champs magnétiques alternatifs à basse fréquence qui traversent les écrans, générant des courants de circulation élevés qui surchauffent l'isolation polymère (XLPE). Ce phénomène accélère l'apparition de structures de type arborescence électrique ou aquatique (water/electrical treeing) et réduit drastiquement la durée de vie prévisionnelle du câble.

Par ailleurs, les transformateurs de courant (TC) et de tension (TP) essuient des erreurs de rapport et de déphasage significatives imputables à la saturation magnétique induite par les composantes de courant subsynchrone, ce qui paralyse ou dégrade sévèrement les performances des relais de protection numériques fondés sur les phaseurs de fréquence fondamentale.

Stratégies Avancées de Mitigation et de Conception en Ingénierie

Afin de garantir la viabilité opérationnelle des réseaux de transport à compensation série massive, l'ingénierie moderne met en œuvre un ensemble de contre-mesures technologiques de pointe, ventilées entre la mitigation au niveau du réseau, au niveau des équipements de production et par le truchement de systèmes électroniques de puissance à grande vitesse.

Filtres de Blocage Subsynchrone (Subsynchronous Resonance Blocking Filters - SSBF)

Les filtres SSBF constituent des circuits accordés installés en série avec les enroulements du stator du générateur ou directement insérés sur la ligne de transport. Ils consistent en un agencement de réacteurs et de condensateurs montés en parallèle et dimensionnés pour présenter une impédance extrêmement élevée (un circuit ouvert virtuel) à la fréquence électrique résonnante critique $fer$ .

L'impédance du filtre SSBF à la fréquence de résonance s'exprime par la relation :

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