Analyse des fréquences naturelles et des transitoires électromagnétiques dans les réseaux de distribution

Fondements Physico-Mathématiques des Transitoires Électromagnétiques et de la Réponse Fréquentielle

L'analyse rigoureuse des transitoires électromagnétiques et des fréquences naturelles dans les réseaux de distribution d'énergie électrique exige une formulation mathématique avancée qui va bien au-delà de l'analyse en régime permanent phasoriel. Les réseaux de distribution modernes, caractérisés par une forte pénétration de câbles souterrains à moyenne tension, de transformateurs à noyaux ferromagnétiques non linéaires et de charges électroniques décentralisées, se comportent comme des systèmes à paramètres répartis et concentrés soumis à des excitations impulsives et harmoniques sévères.

Pour modéliser la propagation des ondes électromagnétiques le long des conducteurs du réseau, il convient d'utiliser les équations aux dérivées partielles couplées des télégraphistes. En considérant une ligne multiconductrice symétrique ou asymétrique de longueur infinitésimale $dx$ , le comportement instantané des tensions v(x,t) et des courants i(x,t) est défini par le système matriciel suivant :

\frac{\partial}{\partial x} \mathbf{v}(x,t) + \mathbf{R} \mathbf{i}(x,t) + \mathbf{L} \frac{\partial}{\partial t} \mathbf{i}(x,t) = \mathbf{0}

\frac{\partial}{\partial x} \mathbf{i}(x,t) + \mathbf{G} \mathbf{v}(x,t) + \mathbf{C} \frac{\partial}{\partial t} \mathbf{v}(x,t) = \mathbf{0}

Où $\mathbf{R}$ , $\mathbf{L}$ , $\mathbf{C}$ et $\mathbf{G}$ représentent respectivement les matrices par unité de longueur de résistance, d'inductance, de capacité et de conductance transversale (fuite diélectrique). En transformant ces équations dans le domaine fréquentiel à l'aide de la transformation de Laplace ou de Fourier, on obtient l'équation de propagation modale :

\frac{d^2}{dx^2} \mathbf{V}(x,s) = \mathbf{Z}(s) \mathbf{Y}(s) \mathbf{V}(x,s) = \mathbf{P}^2(s) \mathbf{V}(x,s)

Où Z(s) = R + sL, Y(s) = G + sC, et P(s) est la matrice de propagation. Les valeurs propres de la matrice produit Z(s)Y(s) déterminent les constantes de propagation modale \gamma_k(\omega) = \alpha_k(\omega) + j\beta_k(\omega), qui régissent l'atténuation et la vitesse de phase des ondes the voyagent. Les fréquences naturelles du système (\omega_n) correspondent aux pôles de la fonction de transfert de l'impédance d'entrée Zin(s) vue depuis n'importe quel nœud du réseau de distribution, en résolvant la condition caractéristique :

\det \left[ \mathbf{I} - \mathbf{M}_{end}(s) \right] = 0

Où M_{end}(s) représente la matrice de réflexion et de transmission aux bornes du réseau, couplant les impédances des transformateurs, des charges et des dispositifs de protection contre les surtensions (SPD).

Classification et Origine des Transitoires Électromagnétiques dans les Réseaux à Moyenne et Basse Tension

Les transitoires électromagnétiques dans les réseaux de distribution proviennent de perturbations opérationnelles et atmosphériques qui injectent des spectres de fréquences à large bande. La catégorisation correcte de ces phénomènes est impérative pour la conception de schémas d'isolation coordonnés selon les normes IEEE Std 1313.2 et IEC 60071-2.

Transitoires d'Origine Atmosphérique (Foudre)

Les décharges atmosphériques directes ou indirectes sur les lignes aériennes de distribution induisent des ondes de surtension caractérisées par des fronts d'onde extrêmement rapides (de l'ordre de 0,5 \mu s à 5 \mu s) et de longues queues (jusqu'à 50 \mu s). Le courant de foudre est typiquement modélisé par la double fonction exponentielle de Heidler :

i(t) = \frac{I_p}{\eta} \cdot \frac{(t/\tau_1)^n}{1 + (t/\tau_1)^n} \cdot e^{-t/\tau_2}

Où $I_p$ est l'amplitude du courant de crête, \tau_1 est la constante de temps du front d'onde, \tau_2 est la constante de temps de la queue de l'onde, $n$ est le facteur de raideur (steepness), et \eta est un facteur de correction d'amplitude. Lorsque ces ondes se propagent depuis des sections aériennes vers des câbles souterrains, l'impédance caractéristique chute drastiquement (passant généralement de Z_c \approx 400\,\Omega à Z_c \approx 30\,\Omega), provoquant une réfraction de tension qui peut doubler l'amplitude de l'onde incidente au point de transition si la longueur du câble est inférieure à la longueur critique du front d'onde.

Manœuvres d'Ouverture et de Fermeture des Disjoncteurs

L'exploitation de disjoncteurs à vide (vacuum circuit breakers) ou de sectionneurs en charge génère des transitoires à haute fréquence en raison du phénomène de réamorçage multiple de l'arc électrique (multiple restrikes). Lorsque les contacts s'écartent, le courant passe par un zéro à haute fréquence et la rigidité diélectrique du milieu de coupure se rétablit à un rythme supérieur à la vitesse de montée de la tension de rétablissement transitoire (TRV). Cela provoque une interruption brusque du courant (current chopping), dont l'énergie inductive emmagasinée est transférée à la capacité parasite du réseau selon la relation énergétique :

\frac{1}{2} Leq ichop^2 = \frac{1}{2} Ceq vmax^2 \implies vmax = ichop \sqrt{\frac{Leq}{Ceq}}

Cette surtension oscillatoire à haute fréquence (généralement comprise entre 1 kHz et 100 kHz) soumet l'isolation principale des enroulements de transformateurs à une contrainte diélectrique sévère.

Résonance et Ferrosonance dans les Systèmes de Distribution

La résonance harmonique et la ferrosonance constituent deux des modes de défaillance les plus destructeurs et complexes dans les réseaux de distribution souterrains et mixtes. Leur analyse nécessite l'évaluation de l'interaction entre les éléments réactifs linéaires et les non-linéarités ferromagnétiques.

Résonance Série et Parallèle

La résonance parallèle (antirésonance) se produit lorsque la réactance inductive équivalente du réseau (fournie par les transformateurs de distribution et les réactances de limitation) égale la réactance capacitive des câbles souterrains et des batteries de condensateurs de correction du facteur de puissance à une fréquence harmonique spécifique $h$ :

\omega_r = \frac{1}{\sqrt{Leq Ceq}} \implies h_r = \frac{f_r}{f_1} = \sqrt{\frac{X_c}{X_L}}

Dans cet état, l'impédance vue de la source atteint une valeur maximale théorique limitée uniquement par la résistance équivalente du réseau. Toute injection de courants harmoniques (provenant de charges non linéaires telles que des variateurs de vitesse et des redresseurs) proche de $h_r$ entraîne des chutes de tension harmoniques massives, une surchauffe par effet Joule et une distorsion sévère de l'onde de tension.

Ferrosonance dans les Transformateurs

La ferrosonance est un phénomène non linéaire caractérisé par des surtensions soutenues à basse fréquence (fondamentale, sous-harmonique ou chaotique), une distorsion harmonique extrême et des courants élevés. Elle se produit typiquement lorsqu'un transformateur à noyau de fer insaturable est mis sous tension à vide à travers une capacité en série (par exemple, la capacité d'un câble souterrain déconnecté à une extrémité ou des phases ouvertes suite à un déclenchement monopolaire de fusibles).

Le comportement non linéaire du noyau magnétique est modélisé par le flux enlacé \lambda(i) ou la relation constitutive de l'inductance non linéaire :

v(t) = R i(t) + \frac{d\lambda(i)}{dt}

La coexistence de multiples états stables (régime fondamental et régimes ferrosonnants sous-harmoniques) s'explique par la théorie des bifurcations dans les systèmes dynamiques non linéaires. L'équation différentielle régissant le circuit LC non linéaire excité par une source sinusoïdale V_m \sin(\omega t) est :

L(i) \frac{d^2 i}{dt^2} + R \frac{d i}{dt} + \frac{1}{C} \int i \, dt = V_m \omega \cos(\omega t)

Où l'inductance L(i) = \frac{d\lambda}{di} varie radicalement lorsque le noyau entre dans sa région de saturation magnétique, ce qui réduit la valeur de $L$ et décale la fréquence de résonance du circuit vers des valeurs qui coïncident avec la fréquence du réseau ou ses sous-harmoniques.

Analyse Forensique des Défauts et Conséquences Diélectriques

Les transitoires électromagnétiques à haute fréquence et les résonances prolongées détériorent de manière cumulative et catastrophique les actifs critiques du réseau de distribution. L'analyse forensique des mécanismes de défaillance par composant est détaillée ci-dessous :

Transformateurs de Distribution et de Puissance

Lorsqu'un front d'onde rapide (généralement généré par la foudre ou des manœuvres) pénètre dans l'enroulement à haute tension d'un transformateur, la distribution initiale de la tension le long des spires n'est pas linéaire. En raison de la présence de la capacité série entre spires ( $C_s$ ) et de la capacité à la terre de l'enroulement ( $C_g$ ), le facteur de distribution \alpha = \sqrt{C_g / C_s} induit un gradient de tension extrêmement élevé sur les premiers disques ou spires de l'enroulement.

Cela engendre des décharges partielles (Partial Discharges - PD) internes dans le papier imprégné d'huile, une perforation de l'isolation cellulosique et des courts-circuits entre spires. Du point de vue thermique, les courants harmoniques induits par résonance parallèle provoquent des pertes supplémentaires dans les enroulements par effet de peau (skin effect) et de proximité, calculées à l'aide du facteur de pertes additionnelles $H$ :

PLL = PDC \left[ 1 + \sum_{h=2}^{n} \left( \frac{I_h}{I_1} \right)^2 H(h) \right]

Cables Souterrains à Moyenne Tension

Les câbles isolés au Polyéthylène Réticulé (PRC / XLPE) se montrent hautement vulnérables aux transitoires de tension prolongés et à la ferrosonance. Les surtensions soutenues accélèrent le processus de treeing électrique (arborescence électrique), un mécanisme de dégradation de l'isolation solide où des microcavités soumises à des champs électriques intenses (> 10 kV/mm) génèrent des décharges partielles continues qui élargissent les canaux dendritiques jusqu'au claquage diélectrique total (rupture de l'isolation primaire).

Appareils de Manœuvre et Protections

Les disjoncteurs à vide et les réenclencheurs subissent une érosion sévère de leurs contacts principaux en raison de l'énergie accumulée dans les arcs à haute fréquence lors des opérations de coupure sur charges inductives ou capacitives mal adaptées. De même, les parafoudres à oxyde de zinc (MOV) peuvent faire l'objet d'une défaillance thermique par emballement énergétique si la densité d'énergie absorbée au cours d'un transitoire dépasse leur capacité de dissipation nominale $Wth$ .

Paramètres de Conception, Limites Normatives et Stratégies d'Atténuation

Afin de garantir la fiabilité et la coordination de l'isolement selon les directives internationales, les ingénieurs doivent appliquer des critères stricts de conception préventive et atténuer les effets des fréquences naturelles.

Le tableau suivant synthétise les paramètres électriques critiques, les limites normatives selon les normes IEEE et IEC, les conditions de défaillance critique ainsi que les conséquences opérationnelles associées :

Paramètre Électrique / Phénomène	Norme de Référence (IEEE / IEC)	Limite Normative / Seuil Critique	Condition de Défaillance Critique	Conséquences Opérationnelles et Diélectriques
Surtensions Temporaires (TOV)	IEEE Std 1313.2 / IEC 60071-2	VTOV \le 1,7 \times V_{ph-n} (pour t \le 1\, s )	Ferrosonnance soutenue ou rejet de charge dans des réseaux faibles	Dégradation thermique et vieillissement prématuré de l'isolation XLPE dans les câbles.
Taux de Distorsion Harmonique Global de Tension (THD-V)	IEEE Std 519-2022 / IEC 61000-2-2	THD-V \le 5,0\% (Réseaux à Moyenne Tension)	Résonance parallèle proche des harmoniques 5e, 7e ou 11e	Surchauffe des transformateurs, déclenchements intempestifs des protections et fatigue diélectrique.
Vitesse de Montée de la Tension (raideur différentielle $dv/dt$ )	IEC 60034-18-41 / NEMA MG1	dv/dt \le 5\, kV /\mu s (aux bornes des moteurs/cables)	Fonctionnement de disjoncteurs à vide avec réamorçages multiples	Perforation de l'isolation des spires et défauts inter-spires dans les transformateurs et machines tournantes.
Capacité d'Absorption d'Énergie des SPD (MOV)	IEEE Std C62.11 / IEC 61643-11	Absorption nominale \ge 2,5\, kJ/kV de tension assignée	Impacts multiples de foudre ou fronts d'onde réfléchis amplifiés	Panne thermique catastrophique (explosion) du parafoudre et perte de protection.
Facteur de Déclassement (Derating) Harmonique des Câbles	IEC 60287 / NEC 310.15	Facteur de réduction $Kharm$ selon le spectre de courant	Courants de neutre supérieurs au courant de phase par composantes homopolaires (3e, 9e)	Fusion de l'isolation par excès de température dans le conducteur (> 90°C pour le XLPE).

Stratégies d'Atténuation et de Compensation Active/Passive

Pour atténuer les fréquences naturelles de résonance et les transitoires sévères, les contromesures d'ingénierie suivantes sont déployées :

Installation de Filtres Harmoniques Sintonisés (Filtres Passifs LC) : Connectés en parallèle avec le réseau pour dévier les harmoniques critiques (typiquement le 5e et le 7e) vers des résistances de dissipation, modifiant ainsi l'impédance équivalente du réseau et évitant la résonance parallèle. Leur réglage (sintonisation) est effectué légèrement en dessous de la fréquence harmonique cible :
$fsintonise = \frac{1}{2\pi \sqrt{Lfilter Cfilter}} < fharmonique$
Réactances d'Amortissement et Résistances de Pré-insertion : Sur les disjoncteurs à haute tension et moyenne tension critique, l'utilisation de résistances de pré-insertion amortit le front d'onde transitoire, ramenant le facteur de surtension à des valeurs $vmax / Vph < 2,0$ .
Suppresseurs de Surtension à Oxyde de Zinc (MOV) Positionnés Stratégiquement : Installés aux points de transition entre lignes aériennes et câbles souterrains, ainsi qu'aux bornes des transformateurs, afin de tronquer les crêtes de tension transitoire et de limiter l'énergie incidente.

Application Pratique et Analyse via Vexten Suite

La modélisation avancée des réseaux de distribution pour l'étude des transitoires électromagnétiques et des fréquences naturelles nécessite des plateformes de calcul informatique de haute précision telles que Vexten Suite. La méthodologie algorithmique et les procédures d'ingénierie mises en œuvre dans le logiciel pour se conformer aux normes internationales sont détaillées ci-après.

Flux de Calcul de Court-Circuit et des Transitoires (IEC 60909 / IEEE 141)

Le moteur de calcul transitoire de Vexten Suite résout le système d'équations différentielles dans le domaine temporel en utilisant la méthode d'intégration numérique trapézoïdale avec un pas de temps adaptatif (\Delta t < 10\,\mu s). La procédure analytique intégrée comprend les étapes suivantes :

Acquisition et Modélisation Topologique : Le système importe les paramètres physiques des câbles (géométrie, couches semi-conductrices, isolation XLPE), des transformateurs (impédances de court-circuit, courant à vide, courbe de saturation de Fröhlich-Kennelly) et des sources de production décentralisée.
Analyse de Fréquence Naturelle (Balayage d'Impédance - Frequency Scan) : Une source de courant harmonique d'amplitude unitaire et de fréquence variable f \in [10\, Hz , 100\, kHz ] est injectée dans les nœuds critiques du réseau. Vexten calcule l'impédance nodale Znn(j\omega) et génère le diagramme de Bode en amplitude et en phase, identifiant automatiquement les pics de résonance parallèle et les creux de résonance série.
Simulation de Manœuvres et de la Foudre (Moteur de type EMTP) : Des événements de commutation sont configurés (ouverture de disjoncteurs avec des courants de coupure spécifiques) ainsi que des ondes de foudre normalisées (8/20 \mu s et 1,2/50 \mu s). Le logiciel évalue la propagation des ondes voyageuses en tenant compte de l'atténuation modale dépendante de la fréquence et de l'effet de peau (skin effect) par le biais de la formulation de Carson modifiée ou de fonctions de Bessel pour conducteurs cylindriques.
Évaluation du Déclassement Thermique des Câbles (IEC 60287) : Sur la base du spectre de courant harmonique obtenu à partir de l'analyse du flux de charge harmonique, Vexten calcule automatiquement le facteur de réduction de la capacité de transport de courant ( $ampacity$ ) du câble souterrain, garantissant ainsi que la température maximale du conducteur ne dépasse pas la limite thermique normative de 90°C pour le XLPE en régime distordu.

De cette manière, Vexten Suite permet à l'ingénieur spécialiste d'effectuer un diagnostic exhaustif, prévenant les défaillances catastrophiques dues à la résonance et assurant une coordination optimale de l'isolement au sein de réseaux de distribution complexes.