Ferroresonancia y Saturación Magnética en Reactores de Derivación por Apertura Monofásica

Introducción y Fundamentos Electromagnéticos de los Reactores de Derivación

Los reactores de derivación conectados en sistemas de transmisión de extra alta tensión (EAT) y ultra alta tensión (UAT) desempeñan un papel crítico en la compensación de la potencia reactiva capacitiva generada por las líneas de transmisión largas bajo condiciones de baja carga o en vacío. Electromecánicamente, estos equipos están constituidos por núcleos ferromagnéticos con entrehierros distribuidos, diseñados para mantener una característica de magnetización lineal en condiciones nominales y de sobretensión temporal moderada. Sin embargo, cuando se ejecutan maniobras de conmutación monofásica, o cuando ocurren aperturas asimétricas y fallas transitorias, la interacción entre la capacitancia distribuida de la línea y la inductancia no lineal del reactor de derivación crea un circuito resonante altamente complejo.

Desde la perspectiva de la teoría de campos electromagnéticos y circuitos acoplados, la inductancia del reactor deja de ser un parámetro invariante en el tiempo para convertirse en una función dependiente del flujo magnético concatenado \lambda y de la corriente $i$ . El comportamiento del núcleo de hierro está gobernado por la ecuación constitutiva:

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

donde los coeficientes $c_n$ reflejan la no linealidad severa debida a la saturación magnética. Durante una apertura monofásica, los polos sanos de la línea de transmisión continúan induciendo tensiones electrostáticas y electromagnéticas sobre la fase desconectada a través de las capacidades mutuas interfasiales $Cm$ y las capacidades de tierra $Cs$ . Si la inductancia del reactor de derivación $L$ forma un lazo resonante con la capacitancia equivalente de la fase abierta $Ceq$ , las condiciones para la ferroresonancia serie o paralela quedan plenamente establecidas, desencadenando sobretensiones sostenidas de baja frecuencia, armónicos de orden sub-sincrónico y corrientes de excitación que superan en varios órdenes de magnitud la corriente nominal del equipo.

Ecuaciones Diferenciales y Modelado No Lineal del Núcleo Saturable

El modelado matemático riguroso de la ferroresonancia en reactores de derivación requiere la formulación de sistemas de ecuaciones diferenciales ordinarias y parciales que acoplen los fenómenos electrostáticos de la red con la histéresis y la saturación del núcleo magnético. Considerando un circuito equivalente monofásico simplificado que representa una fase desconectada de una línea de transmisión con un reactor de derivación conectado en los terminales, la dinámica del sistema se describe mediante la siguiente ecuación integro-diferencial:

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

Donde \lambda es el flujo concatenado en el reactor, $R$ es la resistencia equivalente de pérdidas en el cobre y en el núcleo, $L_0$ es la inductancia en el régimen lineal, \omega_0 es la frecuencia fundamental del sistema, y V_m \sin(\omega t + \theta) representa la tensión inducida desde las fases adyacentes energizadas a través del acoplamiento capacitivo. La función de no linealidad del núcleo magnético f(\lambda) se modela comúnmente mediante la aproximación de Fröhlich modificada o mediante curvas de saturación analíticas basadas en funciones arcotangente hiperbólica:

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

En este sistema, la presencia de la saturación provoca que la inductancia aparente del reactor varíe dinámicamente durante cada ciclo de la tensión aplicada. Cuando el núcleo entra en la región de saturación profunda, la inductancia diferencial del reactor Ldiff = d\lambda / di disminuye drásticamente, acercándose a los valores de la inductancia de dispersión del devanado. Esta reducción repentina desplaza la frecuencia de resonancia del circuito hacia valores cercanos o inferiores a la frecuencia fundamental, induciendo modos de oscilación subarmónicos (típicamente de orden 1/2, 1/3 o 1/5) y cuasi-periódicos.

El análisis de estabilidad de Lyapunov y los mapas de bifurcación demuestran que el sistema dinámico del reactor ferroresonante posee múltiples estados estables (multi-estabilidad). Dependiendo de las condiciones iniciales (como el ángulo de desconexión del interruptor automático y el valor instantáneo del flujo remanente \lambda_r en el núcleo), el sistema puede converger hacia un estado normal de baja tensión o saltar repentinamente a un estado ferroresonante caracterizado por corrientes distorsionadas y flujos magnéticos elevados que superan el punto de rodilla de la curva de saturación en más de un 200%.

Topología del Fenómeno durante Maniobras de Conmutación Monofásica

Las maniobras de apertura y cierre monopolar en interruptores de alta tensión son procedimientos operativos estándar para la eliminación de fallas transitorias en líneas de transmisión. Sin embargo, durante el intervalo de apertura monofásica, la fase afectada queda aislada de la fuente de alimentación principal, pero permanece acoplada electromagnéticamente al resto del sistema energizado. Las capacidades mutuas entre fases $Cab, Cbc, Cca$ y las capacidades hacia tierra $Can, Cbn, Ccn$ forman una red capacitiva compleja que alimenta la fase abierta con una tensión de origen electrostático.

Cuando un reactor de derivación está permanentemente conectado a los terminales de la línea, la inductancia de este reactor se coloca en paralelo con la capacitancia equivalente de secuencia cero y secuencia inversa vista desde los terminales abiertos de la fase. El circuito equivalente de Thévenin visto desde los terminales del reactor consta de una fuente de tensión equivalente $Vth$ impulsada por las tensiones de las fases sanas a través de los condensadores de acoplamiento mutuo, y una impedancia equivalente $Zth$ dominada por la capacitancia equivalente $Ceq$ .

Parámetro de Maniobra	Condición Normal / Operativa	Condición Crítica de Ferroresonancia	Consecuencia Electromecánica
Estado del Interruptor	Cierre trifásico simétrico	Apertura monofásica asimétrica	Desbalance severo de tensiones y corrientes de secuencia homopolar.
Capacitancia Equivalente ( $Ceq$ )	Baja impedancia capacitiva de línea	Acoplamiento óptimo con $Lreactor$	Sintonización del circuito resonante a frecuencias subarmónicas.
Flujo Remanente (\lambda_r)	Bajo o distribuido simétricamente	Alto flujo remanente post-interrupción	Asimetría en el ciclo de histéresis y entrada inmediata en saturación.
Pérdidas del Sistema ( $Req$ )	Amortiguamiento normal por carga	Bajo amortiguamiento en vacío/baja carga	Incapacidad del sistema para disipar la energía oscilatoria acumulada.

Durante la apertura del polo de una fase, el arco eléctrico transitorio en el interruptor puede generar reencendidos de alta frecuencia (re-strikes), inyectando componentes de alta frecuencia transitoria que excitan los modos naturales del conjunto reactor-línea. Si el flujo remanente en la columna del reactor correspondiente a la fase abierta es elevado y coincide en polaridad con el semi-ciclo inicial de la tensión inducida, el núcleo alcanza la saturación magnética profunda en el primer medio ciclo posterior a la extinción del arco, iniciándose la trayectoria hacia la inestabilidad ferroresonante.

Análisis Forense de Fallas y Consecuencias en Activos de Subestaciones

Los eventos de ferroresonancia y saturación magnética extrema en reactores de derivación generan solicitaciones térmicas, dieléctricas y mecánicas severas que comprometen la integridad de los activos de la subestación. El análisis forense de fallas ocurridas en sistemas de EAT revela patrones de daño característicos en los diferentes componentes del sistema eléctrico.

Deformación Mecánica y Esfuerzos Electrodinámicos en los Devanados

Las corrientes de excitación ferroresonante están compuestas masivamente por armónicos de orden impar (principalmente tercer y quinto armónico) y componentes subarmónicas. Estas corrientes pueden alcanzar magnitudes de 3 a 5 veces la corriente nominal de diseño del reactor. Los esfuerzos electromecánicos sobre los conductores del devanado son proporcionales al cuadrado de la corriente instantánea (F \propto i(t)^2). Las fuerzas radiales y axiales cíclicas inducen fatiga mecánica en los soportes de aislamiento, provocando el aflojamiento de los bloques de presión del devanado, la abrasión del aislamiento de papel kraft y, eventualmente, cortocircuitos entre espiras.

Calentamiento Suplementario por Corrientes de Foucault y Pérdidas en el Núcleo

La operación del núcleo magnético en la región de saturación profunda altera la distribución del flujo magnético, forzando la salida del flujo fuera de las chapas de acero al silicio de grano orientado hacia las estructuras de sujeción, pantallas magnéticas y el tanque metálico del reactor. Las pérdidas por corrientes parásitas (Foucault) en las placas de blindaje y en las paredes del tanque se incrementan de forma no lineal:

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

Donde $d$ es el espesor de las laminaciones, $f$ es la frecuencia efectiva (incluyendo armónicos), y $Bmax$ es la densidad de flujo magnético pico. Este incremento masivo de pérdidas genera puntos calientes locales (hot spots) en el aceite aislante, superando los límites térmicos permisibles según las normas IEEE C57.109 e IEC 60076-6, provocando la degradación térmica del aislamiento celulósico, la generación de gases disueltos en el aceite (hidrógeno, metano, acetileno) y el disparo de las protecciones de presión súbita (Buchholz).

Sobretensiones Dieléctricas y Colapso del Aislamiento

Las tensiones ferroresonantes sostenidas pueden alcanzar valores de hasta 2.5 a 3.5 p.u. con formas de onda altamente distorsionadas que presentan frentes de onda abruptos. Estas sobretensiones actúan de manera continua sobre el sistema de aislamiento principal (aceite-papel o aislamiento sólido en seco). Los descargadores de sobretensión de óxido de metal (MOA) ubicados en los terminales del reactor se ven sometidos a una absorción de energía térmica excesiva al intentar limitar estas sobretensiones de baja frecuencia, lo que puede llevar a la falla térmica y explosión catastrófica de los pararrayos si su capacidad de disipación de energía (classe de descarga de línea) es rebasada.

Estrategias Avanzadas de Mitigación y Criterios de Diseño

La mitigación efectiva de la ferroresonancia y la saturación magnética en reactores de derivación requiere la implementación de contramedidas a nivel de diseño electromagnético, selección de equipos de maniobra y configuración de sistemas de protección avanzada.

Optimización del Diseño del Núcleo y Entrehierros

Para evitar la entrada en saturación rápida bajo condiciones de sobretensión temporal y acoplamiento capacitivo, los reactores modernos se diseñan con núcleos acorazados o de columnas con múltiples entrehierros no magnéticos distribuidos uniformemente a lo largo de la altura de las columnas. Esto linealiza la característica de magnetización hasta el 140% - 160% de la tensión nominal. La inductancia linealizada se expresa mediante:

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

Al hacer que la longitud total de los entrehierros \sum lgap sea el término dominante en el denominador, la reluctancia magnética total se vuelve independiente de la permeabilidad relativa del acero \mu_r, suprimiendo la no linealidad y previniendo la ferroresonancia.

Implementación de Resistores de Preinserción y Sincronización de Interruptores

La utilización de interruptores automáticos equipados con mecanismos de control sintonizados por ángulo de fase (Point-on-Wave Switching) permite eliminar los transitorios severos de conexión y desconexión. Al abrir o cerrar los polos del interruptor en los cruces por cero de la corriente o en los instantes óptimos de tensión, se minimiza la generación de flujos remanentes asimétricos. Adicionalmente, el uso de resistores de preinserción en los interruptores de EAT amortigua las oscilaciones transitorias de alta y media frecuencia durante las maniobras de conmutación.

Esquemas de Amortiguamiento y Filtros Pasivos

En subestaciones donde el riesgo ferroresonante es estructuralmente elevado debido a topologías de red complejas con líneas largas en vacío, se implementan redes de amortiguamiento basadas en circuitos de alta resistencia conectados en serie con bancos de condensadores auxiliares o la adición de devanados terciarios en los reactores conectados en delta cerrada con una resistencia de amortiguamiento externa. Esta resistencia disipa la energía oscilatoria de los modos subarmónicos, impidiendo el crecimiento de la amplitud de la oscilación ferroresonante.

Aplicación Práctica y Análisis con la Vexten Suite

Para ilustrar el dimensionamiento y la validación técnica de un sistema de compensación reactiva sujeto a regímenes transitorios y armónicos, se presenta un caso de estudio desarrollado mediante los módulos analíticos de la plataforma Vexten Suite, aplicando los estándares internacionales IEC 60909, IEEE 141 e IEC 60287.

Especificaciones del Sistema de Potencia

Tensión nominal del sistema ( $U_n$ ): 400 kV (Red trifásica, 50 Hz).
Potencia de cortocircuito trifásica de la red ( $Sk3$ ): 25 GVA.
Reactor de derivación trifásico: 150 MVAR, 400/√3 kV por fase, conexión en estrella con neutro rígidamente aterrado.
Capacitancia equivalente de secuencia homopolar de la línea ( $C_0$ ): 12.5 µF por fase.
Resistencia equivalente del devanado del reactor ( $Rcu$ ): 0.15 \Omega por fase a 75 °C.

Cálculo de Cortocircuito y Parámetros Transitorios (IEC 60909 / IEEE 141)

La impedancia equivalente de la red vista desde los terminales del reactor se calcula a partir de la potencia de cortocircuito:

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

La reactancia inductiva equivalente de la red ( $Xnet$ ) asumiendo una relación $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

La inductancia nominal del reactor de derivación por fase se determina mediante la reactancia nominal del equipo:

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

Evaluación del Factor de Frecuencia de Resonancia

Durante la apertura monofásica, la frecuencia de resonancia del circuito formado por la inductancia del reactor $L$ y la capacitancia equivalente de la fase abierta Ceq \approx C_0 se calcula mediante:

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

El valor de $fres = 24.45 Hz$ se encuentra exactamente en la región sub-sincrónica (aproximadamente la mitad de la frecuencia fundamental de 50 Hz), lo que confirma la vulnerabilidad del sistema ante la excitación de subarmónicos de orden 1/2 cuando se producen condiciones asimétricas de conmutación.

Derotación Armónica y Dimensionamiento de Conductores (IEC 60287 / NEC 310)

Como consecuencia de la saturación magnética y los regímenes transitorios ferroresonantes, las corrientes que circulan por los cables de conexión y devanados contienen un contenido armónico significativo. La Vexten Suite evalúa el factor de pérdidas suplementarias y aplica el coeficiente de reducción por armónicos (Harmonic Derating Factor - HDF) conforme a la norma IEC 60287:

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

Donde $M_h$ es la tasa de distorsión armónica de orden $h$ , $Y_h$ es el factor de pérdidas por corrientes de Foucault en el conductor para el armónico $h$ , y \Delta_1 representa el incremento de pérdidas por histéresis y Eddy en el núcleo. Para un contenido de tercer armónico del 25% y quinto armónico del 15% derivados de un evento transitorio severo, el factor de reducción calculado por la herramienta Vexten Suite indica que la capacidad de conducción de corriente admisible del cable debe reducirse en un 18.5% para evitar la degradación térmica del aislamiento polimérico (XLPE).

Módulo Vexten Suite	Norma de Referencia	Parámetro de Entrada Principal	Resultado Analítico / Criterio de Diseño
Vexten Transient Studio	IEC 60076-6 / IEEE C57.109	Matriz de inductancias no lineales y flujos remanentes	Identificación de puntos de bifurcación y sobretensiones ferroresonantes de 2.8 p.u.
Vexten Short-Circuit Pro	IEC 60909 / IEEE 141	Potencia de cortocircuito ( $Sk3$ ) e impedancias de secuencia	Determinación de esfuerzos electrodinámicos pico y corrientes asimétricas.
Vexten Thermal Cable Sizer	IEC 60287 / NEC 310	Espectro armónico de corriente y temperatura ambiente	Cálculo del factor HDF y sobredimensionamiento térmico de cables de potencia.
Vexten Protection Coordination	IEC 60255 / IEEE C37.112	Características de excitación y sobrecorriente temporal	Ajuste de relés de sobrecorriente temporizada y protección contra sobretensión residual.

Conclusiones y Directrices Normativas

La inestabilidad por ferroresonancia y saturación magnética en reactores de derivación durante maniobras de conmutación monofásica representa uno de los fenómenos electromagnéticos más complejos y destructivos en los sistemas eléctricos de potencia de alta tensión. El análisis riguroso exige abandonar los modelos lineales tradicionales y adoptar enfoques de simulación no lineal que incorporen la histéresis, la saturación del núcleo y el acoplamiento electromagnético capacitivo distribuido de las líneas de transmisión.

La implementación de diseños de núcleo con entrehierros múltiples distribuidos, el uso de interruptores sincronizados con control de ángulo de fase, y la validación mediante herramientas de ingeniería avanzada como la Vexten Suite garantizan que las instalaciones de compensación reactiva operen dentro de los márgenes de seguridad establecidos por las normativas internacionales IEC e IEEE. El cumplimiento estricto de estos criterios previene fallas catastróficas, optimiza la inversión en activos de subestaciones y asegura la confiabilidad y estabilidad operativa de las redes de transmisión de energía eléctrica a gran escala.