ferroresonanciatransformadores de medidaestabilidad de sistemassobretensiones transitoriasdiseño de subestaciones

Ferrorésonance dans les transformateurs de mesure : Analyse et atténuation

Analyse technique de la ferrorésonance dans les transformateurs de tension : causes physiques, modélisation de la saturation magnétique et stratégies d'atténuation.

Ing. Francisco Ramírez

Introduction et Fondements Physiques de la Ferrorésonnance

La ferrorésonnance représente l'un des phénomènes transitoires et permanents les plus complexes, destructifs et difficiles à prédire dans les réseaux électriques à moyenne et haute tension. Contrairement à la résonance linéaire classique, où un circuit RLC linéaire présente une impédance minimale ou maximale à une fréquence de sintonisation spécifique, la ferrorésonnance est un phénomène non linéaire caractérisé par la présence d'au moins un élément inductif avec un noyau de fer ferromagnétique fonctionnant dans des régions de saturation magnétique, couplé aux capacitances distribuées ou concentrées du système (typiquement les capacitances de câbles souterrains, de longues lignes aériennes, de condensateurs de couplage ou d'interrupteurs ouverts).

Le comportement non linéaire du noyau magnétique est modélisé mathématiquement par la relation constitutive entre le flux magnétique enlacé \lambda et le courant d'excitation ii, qui cesse d'être une fonction linéaire et est décrit par une courbe d'hystérésis et de saturation pouvant être approchée analytiquement par des fonctions potentielles ou polynomiales d'ordre élevé :

λ(i)=L0i+k=1naki2k+1\lambda(i) = L_0 i + \sum_{k=1}^{n} a_k i^{2k+1}

L0L_0 représente l'inductance dans la région linéaire (non saturée) et les coefficients aka_k déterminent la courbure et la sévérité de la non-linéarité dans la zone de saturation profonde. Lorsque le noyau du transformateur de mesure (qu'il s'agisse d'un transformateur de tension inductif - IT, ou d'un transformateur de potentiel capacitifs - TPC) entre en saturation, son inductance équivalente subit une réduction drastique, tombant à des valeurs proches de l'inductance de fuite de l'enroulement. Cet effondrement brusque de l'inductance transforme radicalement l'impédance du circuit résonnant, déplaçant la fréquence naturelle d'oscillation du système vers des valeurs subharmoniques (f0<fsysteˋmef_0 < f_{système}), harmoniques (f0>fsysteˋmef_0 > f_{système}) ou quasi-périodiques.

Dans les transformateurs de mesure, ce phénomène est typiquement déclenché dans des conditions de fonctionnement anormales mais hautement probables dans les réseaux électriques modernes : ouvertures monophasées de disjoncteurs dans des systèmes fonctionnant avec un neutre isolé ou compensé par une bobine de Petersen, défauts à la terre intermittents, enclenchement à vide de transformateurs à travers des lignes présentant une capacité significative, ou résonances soutenues après l'interruption de courants de court-circuit capacitifs.

Classification et Modes Dynamiques d'Oscillation dans les Transformateurs

La nature hautement non linéaire du circuit ferrorésonnant implique l'existence de multiples états stables possibles pour les mêmes paramètres de circuit (phénomène de ferro-multistabilité). La transition d'un état de fonctionnement normal vers un régime ferrorésonnant destructif dépend de manière critique des conditions initiales, de la magnitude de la tension d'excitation, de la capacitance série-parallèle et des pertes totales du système. Les modes d'oscillation sont formellement classés en quatre catégories fondamentales selon leur contenu spectral et leur signature transitoire :

Mode Fundamental (Résonance à la Fréquence du Réseau)

Caractérisé par des oscillations à la fréquence fondamentale du réseau (5050 ou 60Hz60 Hz). Le courant circulant dans l'enroulement du transformateur contient une composante fondamentale amplifiée massivement (jusqu'à 1010 à 2020 fois le courant d'excitation nominal) et une distorsion harmonique modérée. Les surtensions associées peuvent atteindre des valeurs de 2,02,0 à 3,53,5 par unité (p.u.), provoquant la rupture thermique progressive de l'isolement diélectrique en raison de l'augmentation exponentielle des pertes dans le noyau par hystérésis et par courants de Foucault parasitaires.

Mode Subharmonique (Oscillations à des Fractions de la Fréquence du Réseau)

Les oscillations se produisent à des fréquences sous-multiples exactes de la fréquence fondamentale, typiquement f0/2f_0 / 2, f0/3f_0 / 3 ou f0/nf_0 / n. Ce mode est extrêmement dangereux et déconcertant pour les opérateurs, car il génère des ondes de courant et de tension d'aspect chaotique sur l'oscillogramme. Les tensions induites peuvent dépasser 4,0p.u.4,0 p.u., et les courants alternatifs de basse fréquence forcent le noyau à entrer dans une saturation profonde bidirectionnelle à chaque cycle étendu, générant de graves contraintes mécaniques sur les enroulements en raison des forces électrodynamiques à basse fréquence.

Mode Harmonique (Oscillations à des Multiples Entiers de la Fréquence du Réseau)

Dominé par des fréquences supérieures à la fondamentale (2f0,3f0,5f02f_0, 3f_0, 5f_0). Il se manifeste fréquemment lors de processus d'enclenchement de transformateurs à vide connectés à des réseaux à forte capacité distribuée. Le spectre de fréquences montre des pics proéminents dans les harmoniques pairs et impairs, générant des surtensions transitoires de courte durée mais de très fort gradient dv/dt, ce qui sollicite sévèrement l'isolement entre spires des enroulements du transformateur de mesure.

Mode Quasi-Périodique (Chaotique)

Il représente un régime dynamique non linéaire où la tension et le courant contiennent de multiples fréquences incommensurables (non liées par des nombres entiers). Le système oscille de manière apériodique, sautant erratiquement entre différents niveaux d'énergie magnétique. Ce mode est visualisé dans le plan de phase par des attracteurs étranges et représente le scénario le plus sévère pour l'intégrité structurelle et thermique des transformateurs de potentiel, conduisant à une défaillance catastrophique en quelques minutes.

Analyse du Circuit Équivalent et Équations Différentielles Non Linéaires

Pour modéliser rigoureusement le phénomène, considérons un système monophasé équivalent composé d'une source de tension alternative idéale e(t) = E_m \sin(\omega t), d'une résistance série RR représentant les pertes dans les conducteurs et le réseau, d'une inductance linéaire de ligne LsL_s, d'une capacitance shunt CC (représentant la capacité de la ligne ou du câble par rapport à la terre), et d'un transformateur de tension inductif modélisé par une branche non linéaire reliant le flux \lambda(i) et une résistance de pertes dans le noyau RcR_c.

En appliquant les lois de Kirchhoff aux mailles et aux nœuds du circuit équivalent, on obtient le système d'équations intégro-différentielles non linéaires qui régit la dynamique du système :

LsdiLdt+RiL+vc(t)=Emsin(ωt)L_s \frac{di_L}{dt} + R i_L + v_c(t) = E_m \sin(\omega t)
Cdvcdt=iLiμ(vc)vcRcC \frac{dv_c}{dt} = i_L - i_\mu(v_c) - \frac{v_c}{R_c}

iLi_L est le courant traversant l'inductance série, vcv_c est la tension aux bornes du condensateur et du transformateur, et i_\mu(v_c) est la caractéristique de magnétisation non linéaire du noyau ferromagnétique. En exprimant la relation inverse du flux magnétique en fonction de la tension v_c = \frac{d\lambda}{dt}, le système peut être réduit à une équation différentielle ordinaire non linéaire du second ordre :

d2λdt2+(RLs+1RcC)dλdt+1LsCf(λ)=EmLssin(ωt)\frac{d^2\lambda}{dt^2} + \left(\frac{R}{L_s} + \frac{1}{R_c C}\right) \frac{d\lambda}{dt} + \frac{1}{L_s C} f(\lambda) = \frac{E_m}{L_s} \sin(\omega t)

Où la fonction non linéaire f(\lambda) représente le courant de magnétisation en fonction du flux concaténé. La solution analytique de cette équation est impossible par les méthodes classiques en raison du terme non linéaire, nécessitant l'utilisation de méthodes numériques avancées dans le domaine temporel (comme l'algorithme de Runge-Kutta d'ordre supérieur) ou de méthodes analytiques approchées telles que l'équilibrage harmonique et la théorie des perturbations.

L'analyse de stabilité de Lyapunov et l'étude des exposants de Lyapunov permettent de déterminer les conditions critiques de bifurcation (bifurcations de Hopf et selle-nœud) qui font basculer le système d'un état de fonctionnement stable vers la ferrorésonnance chaotique.

Paramétrage des Facteurs Critiques et Cadres Normatifs IEEE/IEC

Les organismes internationaux de normalisation (IEEE et IEC) établissent des directives strictes pour l'évaluation, la conception et l'essai des transformateurs de mesure face au risque de ferrorésonnance. Le tableau comparatif suivant détaille les paramètres électriques critiques, les seuils normatifs et les conséquences opérationnelles associées :

Paramètre Électrique / Condition Limite Normative (IEC 61869 / IEEE C57.13) Condition Critique de Défaillance Conséquence Opérationnelle et Diélectrique
Rapport Capacitance / Inductance (C/LC/L) C<CcritiqueC < Ccritique (Marge de sécurité > 300%) CC élevée par déconnexion de longues lignes ou de câbles souterrains Déplacement de la fréquence de résonance vers la bande subharmonique.
Facteur de Saturation du Noyau (ns=Bsat/Bopn_s = Bsat / Bop) n_s \geq 2,5 pour un fonctionnement nominal linéaire Fonctionnement avec V>1,5p.u.V > 1,5 p.u. ou rémanence magnétique élevée Effondrement de l'inductance de magnétisation et entrée en régime de saturation profonde.
Densité de Courant d'Excitation I_\mu \leq 2\% à 1,0 p.u. ; \leq 10\% à 1,2 p.u. \ Courants de magnétisation the distordus avec des pics supérieurs à 500\% Surcharge thermique sévère sur l'enroulement primaire et destruction de l'isolement solide.
Surtension Temporaire Transitoire (TOV) Supporter 1,9p.u.1,9 p.u. pendant 30secondes30 secondes (IEC 61869-3) Oscillations ferrorésonnantes soutenues supérieures à 3,0p.u.3,0 p.u. Décharges partielles incontrôlées, perforation du papier imprégné d'huile et explosion.
Résistance d'Amortissement Externe Dissipation thermique nominale calculée pour un régime de défaut continu Absence de résistance d'amortissement dans l'enroulement secondaire ouvert (triangle ouvert) Impossibilité de dissiper l'énergie réactive accumulée, perpétuant la ferrorésonnance.

Ingénierie Forensique des Pannes dans les Transformateurs de Mesure

L'apparition de la ferrorésonnance dans les sous-stations électriques laisse des empreintes physiques et métallurgiques indubitables qui permettent aux ingénieurs experts en analyse forensique de reconstituer la séquence des événements ayant conduit à la défaillance catastrophique de l'équipement. L'analyse forensique va de l'inspection visuelle externe jusqu'à la métallurgie analytique des enroulements.

Déformation Électrodynamique et Thermique des Enroulements

Dans des conditions ferrorésonnantes, les courants de court-circuit interne et d'excitation massive génèrent des forces de Lorentz extrêmes entre les spires concentriques du transformateur de tension. Les enroulements primaires, construits avec des conducteurs en cuivre de section réduite pour supporter des tensions élevées, subissent une compression radiale et des contraintes axiales qui dépassent la limite élastique du matériau, provoquant des courts-circuits entre spires contiguës. Sur le plan thermique, le flux magnétique de fuite et les courants de Foucault induits dans les écrans électrostatiques et les structures métalliques internes génèrent des points chauds (hot spots) avec des températures supérieures à 400°C400 °C, carbonisant l'isolement en papier cellulosique et décomposant l'huile isolante.

Analyse des Gaz Dissous (DGA) Post-Panne

Dans les transformateurs de potentiel immergés dans l'huile, l'analyse chromatographique des gaz dissous dans le fluide diélectrique révèle des schémas de défaillance caractéristiques. La ferrorésonnance prolongée génère une décomposition thermique sévère de l'huile minérale, caractérisée par des concentrations massives de :

  • Hydrogène (H2H _2) : Indicateur de décharges partielles et de contraintes électriques sévères.
  • Méthane (CH4CH _4) et Éthane (C2H6C _2 H _6) : Preuve d'une surchauffe locale modérée à sévère (300°C300 °C à 700°C700 °C).
  • Acétylène (C2H2C _2 H _2) : Marqueur indiscutable d'arcs électriques de forte puissance et de décharges disruptives entre spires ou vers la terre.

Comportement des Éléments de Protection Primaire

Les fusibles limiteurs de courant installés au primaire des transformateurs de tension jouent un rôle critique et souvent complexe lors de la ferrorésonnance. En raison du contenu harmonique élevé et de valeurs efficaces de courant qui n'atteignent pas le seuil de fusion rapide du fusible mais qui dépassent la capacité thermique de dissipation, les éléments fusibles subissent un vieillissement prématuré (« fatigue thermique »). Cela provoque la fusion asymétrique des phases, perpétuant la condition d'ouverture monophasée qui alimente et soutient le phénomène ferrorésonnant dans les phases saines restantes.

Stratégies de Conception Pratique et d'Atténuation Avancée

L'atténuation efficace de la ferrorésonnance nécessite une approche multidisciplinaire combinant le choix judicieux des paramètres de construction, la mise en œuvre de topologies de circuits d'amortissement et l'adoption de philosophies opérationnelles avancées pour le contrôle des disjoncteurs.

Calcul et Sélection des Résistances d'Amortissement (Amortisseurs Anti-Ferrorésonnance)

L'insertion d'une résistance d'amortissement dans le circuit secondaire connecté en triangle ouvert est la technique la plus économique et la plus fiable dans les transformateurs de tension inductifs. La puissance dissipée et la résistance optimale RdRd sont calculées à l'aide de la formulation analytique suivante, basée sur l'impédance caractéristique du système :

Rd=3Vsec2Scrit(1+ω2LeqCeq)Rd = \frac{3 Vsec^2}{Scrit} \left( 1 + \omega^2 Leq Ceq \right)

VsecVsec est la tension nominale secondaire, ScritScrit est la puissance apparente critique de dissipation requise pour absorber l'énergie réactive oscillante, et LeqLeq et CeqCeq représentent respectivement l'inductance équivalente de saturation et la capacitance shunt du système. Le dimensionnement thermique de cette résistance doit garantir qu'elle peut supporter le courant de séquence zéro résultant pendant un temps indéfini sans perdre ses propriétés ohmiques.

Critères de Conception des Transformateurs de Potentiel Capacitifs (TPC / CVT)

Dans les TPC, la présence du diviseur capacitif (C1,C2C_1, C_2) couplé au circuit de compensation électromagnétique inductive (LtL_t) crée un circuit intrinsèquement enclin à la ferrorésonnance. Pour neutraliser ce risque, un circuit de suppression de ferrorésonnance (circuit d'absorption ou ferro-détecteur) est obligatoirement incorporé, connecté au secondaire ou à un enroulement tertiaire auxiliaire du transformateur intermédiaire. Ce circuit est composé d'un élément non linéaire (varistance à oxyde de zinc - MOV ou circuit résonnant RLC) qui introduit une conductivité électrique élevée exclusivement lorsque des tensions contenant des composantes subharmoniques ou des surtensions supérieures à 1,2p.u.1,2 p.u. sont détectées, amortissant instantanément l'oscillation sans interférer avec la précision métrologique du TPC en fonctionnement normal.

Stratégies Opérationnelles de Manœuvre (Déclenchement Séquentiel et Commutation Contrôlée)

Du point de vue de l'automatisation des sous-stations, la prévention doit se concentrer sur l'élimination des causes profondes de la ferrorésonnance :

  • Utilisation de Disjoncteurs Tripolaires Synchronisés : L'ouverture et la fermeture coordonnées des pôles du disjoncteur à l'aide de contrôleurs de manœuvre synchronisée (Controlled Switching) éliminent la possibilité d'ouvertures monophasées prolongées dans les réseaux isolés.
  • Élimination des Réseaux à Vide Critiques : Éviter l'enclenchement de transformateurs de mesure à travers des tronçons excessivement longs de câbles souterrains sans charge.
  • Mise à la Terre Efficace du Neutre : Maintenir une impédance de mise à la terre suffisamment basse dans les transformateurs de puissance principaux pour éviter des déplacements sévères du point neutre qui sur-exciteraient les transformateurs de tension adjacents.

Application Pratique et Analyse Informatique via la Vexten Suite

Pour illustrer la rigueur analytique exigée dans la conception en ingénierie de pointe, un cas d'étude résolu par le moteur de calcul avancé de la plateforme Vexten Suite est présenté, intégrant les normes internationales pour le court-circuit, le dimensionnement des câbles et l'atténuation de la résonance harmonique.

Définition du Scénario de Réseau dans la Vexten Suite

Une sous-station haute à moyenne tension (138kV/13.8kV138 kV / 13.8 kV) alimentant un système dont la puissance de court-circuit triphasée symétrique est de Ssc=2500MVASsc = 2500 MVA est analysée. Le réseau moyenne tension dispose d'un système de longs câbles souterrains polymériques qui fournissent une capacitance totale à la terre par phase de C = 1,2\,\mu F . Un banc de transformateurs de tension inductifs avec une inductance de magnétisation linéaire L0=45HL_0 = 45 H et une inductance en saturation profonde Lsat=0,8HLsat = 0,8 H est installé.

Étape 1 : Validation du Court-Circuit et des Niveaux d'Excitation (IEC 60909 / IEEE 141)

En utilisant le module de court-circuit de la Vexten Suite, on calcule les composantes symétriques et asymétriques du courant de défaut ainsi que le facteur de choc thermique (ktk_t) afin de vérifier la résistance mécanique des transformateurs de mesure face aux événements transitoires :

Ik1=1,05Vn3R12+X12=1,05×1380003×0,122+2,452=33,85kAI_{k1''} = \frac{1,05 \, V_n}{\sqrt{3} \, \sqrt{R_1^2 + X_1^2}} = \frac{1,05 \times 138000}{\sqrt{3} \times \sqrt{0,12^2 + 2,45^2}} = 33,85 kA

Le facteur de puissance de court-circuit du réseau (\cos\varphi = R_1 / Z_1) détermine la décroissance de la composante unidirectionnelle de courant continu, laquelle interagit avec le noyau de fer et peut induire un état initial de rémanence magnétique critique évalué par Vexten à \psirem = 0,85\,\psisat.

Étape 2 : Dimensionnement des Câbles et Déclassement Harmonique (IEC 60287 / NEC 310)

Lors d'un événement ferrorésonnant, le spectre de courant est sévèrement pollué par des harmoniques de rang bas et moyen. Le module de dimensionnement des câbles de la Vexten Suite calcule le facteur de réduction dû au contenu harmonique (FharmFharm) applicable aux câbles de commande et de puissance associés aux transformateurs de mesure, garantissant que la température de fonctionnement de l'isolement XLPE ne dépasse pas les 90°C90 °C nominaux :

Fharm=11+h=2n(μh)2(h)2Fharm = \sqrt{ \frac{1}{1 + \sum_{h=2}^{n} (\mu_h)^2 (h)^2} }

Où \mu_h est le rapport harmonique du courant d'ordre hh par rapport à la fondamentale. Pour un spectre typique obtenu lors d'une simulation de ferrorésonnance où le troisième harmonique atteint 35\% et le cinquième 20\%, la Vexten Suite détermine un facteur de déclassement thermique strict de Fharm=0,74Fharm = 0,74, exigeant un surdimensionnement immédiat de la section transversale des conducteurs secondaires d'instrumentation pour éviter la dégradation thermique de l'isolement.

Étape 3 : Analyse de Résonance et Sintonisation de l'Amortissement

Le moteur d'analyse dynamique non linéaire de la Vexten Suite exécute l'intégration numérique dans le domaine temporel du système de circuit avec les paramètres saisis. Les résultats montrent que sans le réseau d'amortissement, la tension aux bornes du transformateur de tension atteint un pic soutenu de 3,42p.u.3,42 p.u. en mode subharmonique d'ordre 1/31/3, avec un courant d'excitation dépassant 45A45 A efficaces.

En appliquant les routines d'optimisation de la Vexten Suite, on calcule la résistance d'amortissement optimale pour le secondaire en triangle ouvert :

Rd,opt=125,4Ω(avecdissipationthermiquenominalePtherm=1,8kW)R_{d,opt} = 125,4\,\Omega \quad ( avec dissipation thermique nominale Ptherm = 1,8 kW )

L'incorporation virtuelle de cette valeur dans le modèle de simulation effondre instantanément l'oscillation ferrorésonnante en un intervalle inférieur à 2,5cycles2,5 cycles de la fréquence fondamentale, ramenant le système à son état de fonctionnement stable, avec des tensions comprises dans les marges métrologiques et protégeant de manière absolue l'intégrité diélectrique et thermique du transformateur de mesure.

Conclusions d'Ingénierie

La ferrorésonnance dans les transformateurs de mesure constitue un phénomène multiphysique critique qui défie les méthodes traditionnelles d'analyse linéaire. Sa compréhension exige une maîtrise approfondie de la théorie des circuits non linéaires, de la dynamique des réseaux électriques et du comportement thermoductif des matériaux ferromagnétiques. Par l'application rigoureuse des normes internationales IEC et IEEE, la conception préventive basée sur l'amortissement secondaire et l'utilisation d'outils de simulation avancés tels que la Vexten Suite, les ingénieurs électriciens peuvent garantir la fiabilité, la sécurité et la continuité d'exploitation des installations électriques à haute et moyenne tension face à ce régime transitoire complexe.