Ferrorésonnance et saturation magnétique dans les réacteurs shunt lors de la commutation monophasée

Introduction et Fondements Électromagnétiques des Réactances de Compensation

Les réactances de compensation (reactors de dérivation) connectées aux systèmes de transmission à très haute tension (THT) et à ultra haute tension (UHT) jouent un rôle critique dans la compensation de la puissance réactive capacitive générée par les longues lignes de transmission en régime de faible charge ou à vide. Électromécaniquement, ces équipements sont constitués de noyaux ferromagnétiques à entrefairs distribués, conçus pour maintenir une caractéristique de magnétisation linéaire dans des conditions nominales et de surtension temporaire modérée. Cependant, lors de manœuvres de commutation monophasée, ou lors d'ouvertures asymétriques et de défauts transitoires, l'interaction entre la capacité distribuée de la ligne et l'inductance non linéaire de la réactance de compensation crée un circuit résonant hautement complexe.

Du point de vue de la théorie des champs électromagnétiques et des circuits couplés, l'inductance de la réactance cesse d'être un paramètre invariant dans le temps pour devenir une fonction dépendante du flux magnétique totalisé \lambda et du courant $i$ . Le comportement du noyau en fer est régi par l'équation constitutive :

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

où les coefficients $c_n$ reflètent la non-linéarité sévère due à la saturation magnétique. Pendant une ouverture monophasée, les pôles sains de la ligne de transmission continuent d'induire des tensions électrostatiques et électromagnétiques sur la phase déconnectée à travers les capacités mutuelles interfaciales $Cm$ et les capacités de terre $Cs$ . Si l'inductance de la réactance de compensation $L$ forme une boucle résonante avec la capacité équivalente de la phase ouverte $Ceq$ , les conditions de la féroresonance série ou parallèle sont pleinement établies, déclenchant des surtensions soutenues à basse fréquence, des harmoniques d'ordre sous-synchrone et des courants d'excitation qui dépassent de plusieurs ordres de grandeur le courant nominal de l'équipement.

Équations Différentielles et Modélisation Non Linéaire du Noyau Saturable

La modélisation mathématique rigoureuse de la féroresonance dans les réactances de compensation nécessite la formulation de systèmes d'équations différentielles ordinaires et aux dérivées partielles qui acouplent les phénomènes électrostatiques du réseau avec l'hystérésis et la saturation du noyau magnétique. En considérant un circuit équivalent monophasé simplifié représentant une phase déconnectée d'une ligne de transmission avec une réactance de compensation connectée aux bornes, la dynamique du système est décrite par l'équation intégro-différentielle suivante :

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

Où \lambda est le flux totalisé dans la réactance, $R$ est la résistance équivalente des pertes dans le cuivre et dans le noyau, $L_0$ est l'inductance en régime linéaire, \omega_0 est la fréquence fondamentale du système, et V_m \sin(\omega t + \theta) représente la tension induite par les phases adjacentes sous tension à travers le couplage capacitif. La fonction de non-linéarité du noyau magnétique f(\lambda) est couramment modélisée par l'approximation de Fröhlich modifiée ou par des courbes de saturation analytiques basées sur des fonctions arc-tangente hyperbolique :

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

Dans ce système, la présence de la saturation entraîne une variation dynamique de l'inductance apparente de la réactance au cours de chaque cycle de la tension appliquée. Lorsque le noyau entre dans la région de saturation profonde, l'inductance différentielle de la réactance Ldiff = d\lambda / di diminue drastiquement, s'approchant des valeurs de l'inductance de fuite de l'enroulement. Cette réduction soudaine déplace la fréquence de résonance du circuit vers des valeurs proches ou inférieures à la fréquence fondamentale, induisant des modes d'oscillation sous-harmoniques (typiquement d'ordre 1/2, 1/3 ou 1/5) et quasi-périodiques.

L'analyse de la stabilité de Lyapunov et les diagrammes de bifurcation démontrent que le système dynamique de la réactance féroresonante possède de multiples états stables (multistabilité). Selon les conditions initiales (telles que l'angle de coupure du disjoncteur et la valeur instantanée du flux rémanent \lambda_r dans le noyau), le système peut converger vers un état normal de basse tension ou basculer soudainement vers un état féroresonant caractérisé par des courants distordus et des flux magnétiques élevés qui dépassent le point de coude de la courbe de saturation de plus de 200%.

Topologie du Phénomène lors des Manœuvres de Commutation Monophasée

Les manœuvres d'ouverture et de fermeture monopolaires dans les disjoncteurs à haute tension sont des procédures opérationnelles standard pour l'élimination des défauts transitoires sur les lignes de transmission. Cependant, pendant l'intervalle d'ouverture monophasée, la phase affectée reste isolée de la source d'alimentation principale, tout en restant couplée électromagnétiquement au reste du réseau sous tension. Les capacités mutuelles entre phases $Cab, Cbc, Cca$ et les capacités à la terre $Can, Cbn, Ccn$ forment un réseau capacitif complexe qui alimente la phase ouverte avec une tension d'origine électrostatique.

Lorsqu'une réactance de compensation est connectée de manière permanente aux bornes de la ligne, l'inductance de cette réactance est placée en parallèle avec la capacité équivalente de séquence zéro et de séquence inverse vue depuis les bornes ouvertes de la phase. Le circuit équivalent de Thévenin vu depuis les bornes de la réactance consiste en une source de tension équivalente $Vth$ alimentée par les tensions des phases saines à travers les condensateurs de couplage mutuel, et en une impédance équivalente $Zth$ dominée par la capacité équivalente $Ceq$ .

Paramètre de Manouvre	Condition Normale / Opérationnelle	Condition Critique de Féroresonance	Conséquence Électromécanique
État du Disjoncteur	Fermeture triphasée symétrique	Ouverture monophasée asymétrique	Déséquilibre sévère des tensions et des courants de séquence homopolaire.
Capacité Équivalente ( $Ceq$ )	Faible impédance capacitive de ligne	Couplage optimal avec $Lreactor$	Synthonisation du circuit résonant sur des fréquences sous-harmoniques.
Flux Rémanent (\lambda_r)	Faible ou distribué symétriquement	Flux rémanent élevé post-interruption	Asymétrie du cycle d'hystérésis et entrée immédiate en saturation.
Pertes du Système ( $Req$ )	Amortissement normal par la charge	Faible amortissement à vide / faible charge	Incapacité du système à dissiper l'énergie oscilatoire accumulée.

Pendant l'ouverture du pôle d'une phase, l'arc électrique transitoire dans le disjoncteur peut générer des réallumages à haute fréquence (ré-amorçages ou restrikes), injectant des composantes transitoires à haute fréquence qui excitent les modes naturels de l'ensemble réactance-ligne. Si le flux rémanent dans la colonne de la réactance correspondant à la phase ouverte est élevé et coïncide en polarité avec le demi-cycle initial de la tension induite, le noyau atteint la saturation magnétique profonde dès le premier demi-cycle suivant l'extinction de l'arc, amorçant la trajectoire vers l'instabilité féroresonante.

Analyse Forensique des Défauts et Conséquences sur les Actifs des Postes Électriques

Les événements de féroresonance et de saturation magnétique extrême dans les réactances de compensation génèrent des contraintes thermiques, diélectriques et mécaniques sévères qui compromettent l'intégrité des actifs du poste électrique. L'analyse forensique des avaries survenues dans les réseaux à THT révèle des motifs de dommages caractéristiques sur les différents composants du système électrique.

Déformation Mécanique et Contraintes Électrodynamiques dans les Enroulements

Les courants d'excitation féroresonante sont composés massivement d'harmoniques d'ordre impair (principalement le troisième et le cinquième harmonique) et de composantes sous-harmoniques. Ces courants peuvent atteindre des magnitudes de 3 à 5 fois le courant nominal de conception de la réactance. Les efforts électromécaniques sur les conducteurs de l'enroulement sont proportionnels au carré du courant instantané (F \propto i(t)^2). Les forces radiales et axiales cycliques induisent une fatigue mécanique dans les supports d'isolation, provoquant le desserrage des blocs de serrage de l'enroulement, l'abrasion de l'isolation en papier kraft et, éventuellement, des courts-circuits entre spires.

Échauffement Supplémentaire par Courants de Foucault et Pertes dans le Noyau

Le fonctionnement du noyau magnétique dans la région de saturation profonde altère la distribution du flux magnétique, forçant la sortie du flux hors des tôles d'acier au silicium à grains orientés vers les structures de serrage, les écrans magnétiques et la cuve métallique de la réactance. Les pertes par courants de Foucault dans les plaques de blindage et sur les parois de la cuve augmentent de manière non linéaire :

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

Où $d$ est l'épaisseur des tôles, $f$ est la fréquence effective (y compris les harmoniques), et $Bmax$ est la densité de flux magnétique de crête. Cette augmentation massive des pertes génère des points chauds locaux (hot spots) dans l'huile isolante, dépassant les limites thermiques admissibles selon les normes IEEE C57.109 et IEC 60076-6, entraînant la dégradation thermique de l'isolation cellulosique, la génération de gaz dissous dans l'huile (hydrogène, méthane, acétylène) et le déclenchement des protections de pression sainte (relais Buchholz).

Surtensions Diélectriques et Effondrement de l'Isolation

Les tensions féroresonantes soutenues peuvent atteindre des valeurs allant jusqu'à 2,5 à 3,5 p.u. avec des formes d'onde fortement distordues présentant des fronts d'onde abrupts. Ces surtensions agissent de manière continue sur le système d'isolation principal (huile-papier ou isolation solide sèche). Les parafoudres à oxyde métallique (MOA) situés aux bornes de la réactance sont soumis à une absorption d'énergie thermique excessive en tentant de limiter ces surtensions à basse fréquence, ce qui peut conduire à la défaillance thermique et à l'explosion cataclysmique des parafoudres si leur capacité de dissipation d'énergie (classe de décharge de ligne) est dépassée.

Stratégies Avancées de Atténuation et Critères de Conception

L'atténuation efficace de la féroresonance et de la saturation magnétique dans les réactances de compensation nécessite la mise en œuvre de contre-mesures au niveau de la conception électromagnétique, de la sélection des appareils de manœuvre et de la configuration des systèmes de protection avancés.

Optimisation de la Conception du Noyau et des Entrefairs

Pour éviter l'entrée en saturation rapide dans des conditions de surtension temporaire et de couplage capacitif, les réactances modernes sont conçues avec des noyaux blindés ou à colonnes dotés de multiples entrefairs non magnétiques distribués uniformément sur la hauteur des colonnes. Cela linéarise la caractéristique de magnétisation jusqu'à 140% - 160% de la tension nominale. L'inductance linéarisée s'exprime par :

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

En faisant en sorte que la longueur totale des entrefairs \sum lgap soit le terme dominant au dénominateur, la réluctance magnétique totale devient indépendante de la perméabilité relative de l'acier \mu_r, supprimant la non-linéarité et prévenant la féroresonance.

Mise en Œuvre de Résistances de Pré-insertion et Synchronisation des Disjoncteurs

L'utilisation de disjoncteurs équipés de mécanismes de commande synchronisés par angle de phase (Point-on-Wave Switching) permet d'éliminer les transitoires sévères de mise sous tension et de mise hors tension. En ouvrant ou en fermant les pôles du disjoncteur lors des passages par zéro du courant ou aux instants optimaux de la tension, la génération de flux rémanents asymétriques est minimisée. De plus, l'utilisation de résistances de pré-insertion dans les disjoncteurs THT amortit les oscillations transitoires à haute et moyenne fréquence lors des manœuvres de commutation.

Schémas d'Amortissement et Filtres Passifs

Dans les postes où le risque féroresonant est structurellement élevé en raison de topologies de réseau complexes avec de longues lignes à vide, des réseaux d'amortissement basés sur des circuits à haute résistance connectés en série avec des batteries de condensateurs auxiliaires sont mis en œuvre, ou bien l'adjonction d'enroulements tertiaires dans les réactances connectés en triangle fermé avec une résistance d'amortissement externe. Cette résistance dissipe l'énergie oscillatoire des modes sous-harmoniques, empêchant l'accroissement de l'amplitude de l'oscillation féroresonante.

Application Pratique et Analyse avec la Suite Vexten

Pour illustrer le dimensionnement et la validation technique d'un système de compensation réactive soumis à des régimes transitoires et harmoniques, un cas d'étude est présenté, développé via les modules analytiques de la plateforme Vexten Suite, en appliquant les normes internationales IEC 60909, IEEE 141 et IEC 60287.

Spécifications du Système de Puissance

Tension nominale du système ( $U_n$ ) : 400 kV (Réseau triphasé, 50 Hz).
Puissance de court-circuit triphasée du réseau ( $Sk3$ ) : 25 GVA.
Réactance de compensation triphasée : 150 MVAR, 400/√3 kV par phase, connexion en étoile avec neutre rigidement mis à la terre.
Capacité équivalente de séquence homopolaire de la ligne ( $C_0$ ) : 12,5 µF par phase.
Résistance équivalente de l'enroulement de la réactance ( $Rcu$ ) : 0,15 \Omega par phase à 75 °C.

Calcul de Court-circuit et Paramètres Transitoires (IEC 60909 / IEEE 141)

L'impédance équivalente du réseau vue depuis les bornes de la réactance est calculée à partir de la puissance de court-circuit :

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

La réactance inductive équivalente du réseau ( $Xnet$ ) en supposant un rapport $X/R = 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

L'inductance nominale de la réactance de compensation par phase est déterminée par la réactance nominale de l'équipement :

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

Évaluation du Facteur de Fréquence de Résonance

Pendant l'ouverture monophasée, la fréquence de résonance du circuit formé par l'inductance de la réactance $L$ et la capacité équivalente de la phase ouverte Ceq \approx C_0 est calculée par :

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

La valeur de fres = 24.45 \ Hz se situe exactement dans la région sous-synchrone (environ la moitié de la fréquence fondamentale de 50 Hz), ce qui confirme la vulnérabilité du système face à l'excitation de sous-harmoniques d'ordre 1/2 lors de conditions de commutation asymétriques.

Déclassement Harmonique et Dimensionnement des Conducteurs (IEC 60287 / NEC 310)

En raison de la saturation magnétique et des régimes transitoires féroresonants, les courants circulant dans les câbles de raccordement et les enroulements contiennent un contenu harmonique significatif. La Suite Vexten évalue le facteur de pertes supplémentaires et applique le coefficient de réduction pour harmoniques (Harmonic Derating Factor - HDF) conformément à la norme IEC 60287 :

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

Où $M_h$ est le taux de distorsion harmonique d'ordre $h$ , $Y_h$ est le facteur de pertes par courants de Foucault dans le conducteur pour l'harmonique $h$ , et \Delta_1 représente l'accroissement des pertes par hystérésis et Eddy dans le noyau. Pour un contenu de troisième harmonique de 25% et de cinquième harmonique de 15% issus d'un événement transitoire sévère, le facteur de réduction calculé par l'outil Vexten Suite indique que la capacité de transport de courant admissible du câble doit être réduite de 18,5% pour éviter la dégradation thermique de l'isolation polymérique (XLPE).

Module Vexten Suite	Norme de Référence	Paramètre d'Entrée Principal	Résultat Analytique / Critère de Conception
Vexten Transient Studio	IEC 60076-6 / IEEE C57.109	Matrice d'inductances non linéaires et flux rémanents	Identification des points de bifurcation et surtensions féroresonantes de 2,8 p.u.
Vexten Short-Circuit Pro	IEC 60909 / IEEE 141	Puissance de court-circuit ( $Sk3$ ) et impédances de séquence	Détermination des efforts électrodynamiques de crête et courants asymétriques.
Vexten Thermal Cable Sizer	IEC 60287 / NEC 310	Spectre harmonique de courant et température ambiante	Calcul du facteur HDF et surdimensionnement thermique des câbles de puissance.
Vexten Protection Coordination	IEC 60255 / IEEE C37.112	Caractéristiques d'excitation et surintensité temporaire	Réglage des relais de surintensité temporisée et protection contre les surtensions résiduelles.

Conclusions et Directrices Normatives

L'instabilité par féroresonance et saturation magnétique dans les réactances de compensation lors de manœuvres de commutation monophasée représente l'un des phénomènes électromagnétiques les plus complexes et destructifs dans les réseaux électriques de puissance à haute tension. L'analyse rigoureuse exige d'abandonner les modèles linéaires traditionnels et d'adopter des approches de simulation non linéaire qui incorporent l'hystérésis, la saturation du noyau et le couplage électromagnétique capacitif distribué des lignes de transmission.

La mise en œuvre de conceptions de noyaux à entrefairs multiples distribués, l'utilisation de disjoncteurs synchronisés avec contrôle d'angle de phase, et la validation par des outils d'ingénierie avancés tels que la Suite Vexten garantissent que les installations de compensation réactive opèrent dans les marges de sécurité établies par les normes internationales IEC et IEEE. Le respect strict de ces critères prévient les pannes catastrophiques, optimise l'investissement dans les actifs des postes électriques et assure la fiabilité et la stabilité opérationnelle des réseaux de transport d'énergie électrique à grande échelle.