Ferroresonancia en Reactores de Derivación por Conmutación Monofásica

Introducción y Marco Teórico de los Fenómenos No Lineales en Sistemas de Alta Tensión

La integración de sistemas de transmisión en corriente alterna de ultra alta tensión (UHVAC) y extra alta tensión (EHVAC) introduce complejidades operativas significativas asociadas a la gestión de la potencia reactiva capacitiva generada por las líneas de transmisión en vacío o con carga ligera. Para mitigar las sobretensiones a frecuencia industrial y controlar el perfil de tensión estacional, el uso de reactores de derivación (shunt reactors), ya sea de núcleo de hierro con entrehierros distribuidos o de núcleo de aire apantallados, es una práctica estándar indispensable. No obstante, la interacción entre la inductancia no lineal del núcleo magnético del reactor y la capacitancia de secuencia positiva y cero de la línea de transmisión crea un circuito LC resonante expuesto a fenómenos de ferroresonancia y saturación magnética transitoria severa durante maniobras operativas desbalanceadas, particularmente la conmutación monofásica.

Desde la perspectiva de la teoría de circuitos electromagnéticos avanzados, la ferroresonancia no es una condición de resonancia lineal clásica, sino un modo de oscilación no lineal caracterizado por la coexistencia de múltiples estados estables (multi-estabilidad) para los mismos parámetros de red. El sistema puede saltar abruptamente desde un estado de baja tensión y corriente cuasi-sinusoidal a un estado ferroresonante de alta tensión, caracterizado por sobretensiones temporales (TOV) sostenidas, distorsión armónica masiva, corrientes de excitación extremadamente elevadas y solicitaciones térmicas y dieléctricas extremas sobre los devanados del reactor y los equipos asociados de la subestación.

El análisis formal de este fenómeno requiere la formulación matemática rigurosa del flujo magnético en el núcleo del reactor acoplado con las ecuaciones diferenciales estocásticas y deterministas del sistema eléctrico trifásico con desequilibrios introducidos por la apertura o cierre secuencial de polos en interruptores automáticos de alta tensión (GIS o AIS).

Ecuaciones Gobernantes y Modelado No Lineal del Núcleo Magnético

Para modelar con precisión la saturación magnética y la ferroresonancia, es fundamental representar la característica de magnetización no lineal del núcleo de hierro del reactor de derivación. La relación entre el flujo magnético encadenado \lambda(t) y la corriente de excitación i(t) se expresa mediante funciones analíticas avanzadas que incorporan los efectos de histéresis y las corrientes parásitas (pérdidas en el núcleo). Una aproximación clásica utilizada en estudios de transitorios electromagnéticos (EMTP) se basa en la función de Fröhlich modificada o en expansiones en serie de potencias impares para evitar discontinuidades:

i(t) = a \lambda(t) + b \lambda(t)^{2n+1}

Donde $a$ representa la inversa de la inductancia diferencial en la zona lineal ( $L_0^{-1}$ ), $b$ es el coeficiente de saturación del núcleo, y $n$ es un entero positivo que determina la agudeza de la rodilla de saturación. La ecuación diferencial que rige el circuito serie equivalente compuesto por la resistencia de la línea $R$ , la inductancia parásita o la del propio reactor con su característica no lineal, y la capacitancia efectiva del sistema $C$ vista desde los terminales del reactor se formula como sigue:

\frac{d^2\lambda}{dt^2} + \frac{R}{L(\lambda)} \frac{d\lambda}{dt} + \frac{1}{C L(\lambda)} \lambda = \frac{1}{C} v_s(t)

Donde v_s(t) = V_m \sin(\omega t + \phi) representa la fuente de tensión sinusoidal del sistema, y la inductancia aparente o diferencial se define como L(\lambda) = \frac{d\lambda}{di}. Cuando el núcleo entra en saturación profunda, L(\lambda) disminuye drásticamente hacia valores cercanos a la inductancia de aire del devanado ( $Lair$ ), lo que altera de manera abrupta la frecuencia de resonancia natural del circuito:

\omega_0(\lambda) = \frac{1}{\sqrt{C L(\lambda)}}

Durante una maniobra de conmutación monofásica (por ejemplo, la apertura de una sola fase de un interruptor debido a la operación unipolares de recierre o una falta monofásica despejada), las fases sanas continúan acopladas capacitivamente a la fase abierta. El reactor conectado a la fase abierta queda flotando a través de una red capacitiva acoplada a las otras dos fases energizadas y a tierra. Esto establece un circuito sintonizado donde la capacitancia equivalente $Ceq$ y la inductancia no lineal L(\lambda) satisfacen las condiciones para la excitación de modos subarmónicos (\omega_0 = \omega/2, \omega/3, \dots), armónicos (3\omega, 5\omega, etc.) o cuasi-periódicos.

Análisis de Maniobras de Conmutación Monofásica y Origen del Desequilibrio

Las maniobras de apertura y cierre monopolar en sistemas EHVAC y UHVAC son operativas o de protección por naturaleza (ej. recierre automático unipolar tras fallas monofásicas temporales). En el intervalo de tiempo en el que una fase está abierta en el interruptor pero el reactor permanece conectado en el lado de la línea, la topología de la red cambia drásticamente. Las capacitancias mutuas entre fases ( $Cab, Cbc, Cca$ ) y las capacitancias a tierra ( $Can, Cbn, Ccn$ ) juegan un papel crítico.

Considerando una línea de transmisión transpuesta de longitud $l$ , la matriz de capacitancias por unidad de longitud se transforma, al abrirse la fase $A$ , en una red Thévenin equivalente vista desde los terminales del reactor de la fase $A$ . La tensión inducida electrostática y electromagnéticamente desde las fases sanas $B$ y $C$ energiza el circuito resonante. Si la reactancia inductiva del reactor a la frecuencia industrial se aproxima o cruza el valor de la reactancia capacitiva equivalente de la red desequilibrada, el sistema entra en una zona de inestabilidad severa.

Parámetro de Operación / Red	Condición Normal (Simétrica)	Condición de Conmutación Monofásica	Límite Normativo (IEC / IEEE)
Tensión de Secuencia Cero ( $V_0$ )	Prácticamente 0% a 1% de $V_1$	Puede alcanzar hasta 30% - 50% de $V_1$	IEEE Std 1036: Máx 5% bajo condiciones asimétricas
Corriente de Excitación del Reactor	Corriente nominal de diseño (I_n \le 1.0 p.u. )	Picos de hasta 8 - 15 p.u. por saturación profunda	IEC 60289: Tolerancia térmica de sobrecorriente transitoria
Distorsión Armónica Total de Tensión (THDv)	< 1.5\% en barras EHV	15\% - 40\% con predominio de armónicos 2°, 3° y 5°	IEEE Std 519: Límite estricto de 1.5% a 2.0% en EHV
Factor de Sobretensión Temporal (TOV)	$1.0 p.u.$	$1.8 p.u. - 2.6 p.u.$ sostenidos por varios ciclos	IEC 60071-2: Curvas de soportabilidad temporal a 10s - 30s

Análisis Forense de Fallas y Consecuencias en Equipos de Subestación

Cuando un reactor de derivación experimenta ferroresonancia sostenida inducida por conmutación monofásica, las consecuencias físicas y dieléctricas sobre los activos de la subestación son catastróficas si las protecciones no actúan con la velocidad adecuada. A continuación se desglosan los mecanismos de falla forense por cada tipo de equipo:

Fallas en Devanados de Reactores de Derivación

Las corrientes de saturación altamente distorsionadas generan fuerzas electrodinámicas radiales y axiales masivas sobre los conductores del devanado. Estas fuerzas superan los límites elásticos del cobre y de los soportes mecánicos, provocando el desplazamiento físico de las espiras, la deformación de los canales de refrigeración y la abrasión del aislamiento de papel celulósico. Asimismo, las pérdidas suplementarias por corrientes parásitas en los devanados y en los elementos estructurales metálicos del tanque debido al contenido armónico elevado (h \ge 3) producen puntos calientes locales (hotspots) con temperaturas que superan los 250 °C, degradando rápidamente el aceite aislante y generando gases combustibles característicos ( $CO, CO₂, CH_4, C_2H_2$ ).

Deterioro Dieléctrico en Cables de Alta Tensión y Pasatapas (Bushing)

Las sobretensiones temporales (TOV) ferroresonantes con componentes de alta frecuencia somete el aislamiento polimérico (XLPE) o impregnado en aceite de los cables de conexión y las pasatapas a esfuerzos dieléctricos repetitivos. Estos esfuerzos aceleran los fenómenos de descarga parcial (PD) interna, erosionando la estructura molecular del aislamiento hasta provocar la perforación dieléctrica en forma de árbol eléctrico (electrical treeing) y la consiguiente falla a tierra catastrófica con explosión del bushing.

Comportamiento Anómalo en Interruptores Automáticos (Circuit Breakers)

Durante la interrupción de corrientes con alto contenido armónico y componente de corriente continua (DC) decareciente generadas por el régimen ferroresonante, los interruptores SF6 enfrentan severas dificultades en el paso por cero de la corriente. Esto incrementa la probabilidad de reencendidos (reignitions) o reamoriamientos (restrikes), inyectando transitorios de alta frecuencia adicionales que agravan el colapso dieléctrico general de la red.

Estrategias de Diseño Práctico y Mitigación Avanzada

La mitigación efectiva de la inestabilidad por ferroresonancia y saturación magnética en reactores de derivación exige la implementación de criterios de diseño multidisciplinarios tanto en la fase de ingeniería de la red como en la especificación de los equipos. Las fórmulas y factores de diseño críticos se detallan a continuación:

Selección de la Característica de Linealidad del Núcleo

Para evitar la saturación prematura ante desequilibrios monofásicos, el diseño del reactor debe garantizar que la rodilla de la curva de saturación magnética se sitúe al menos entre el 140% y el 160% de la tensión nominal de operación continua:

V_{ knee } \ge 1.45 \cdot V_{n\_phase}

Adicionalmente, el uso de reactores con núcleo de aire apantallados elimina por completo la no linealidad introducida por los núcleos ferromagnéticos, eliminando el riesgo de ferroresonancia. Sin embargo, su mayor tamaño, peso y campo magnético disperso requieren estudios detallados de compatibilidad electromagnética (EMC).

Implementación de Resistencias de Amortiguamiento y Descargadores de Sobretensión

Cuando se emplean reactores de núcleo de hierro, la instalación de descargadores de sobretensión de óxido metálico (MO Surge Arresters) de alta capacidad energética directamente en los terminales del reactor ayuda a limitar las crestas de tensión ferroresonante. La energía absorbida por el descargador durante un evento transitorio prolongado se calcula mediante la integral de Joule:

W_{ surge } = \int0^{t_{ clear }} v(t) i(t) dt

Asimismo, la incorporación de un esquema de puesta a tierra a través de una resistencia de amortiguamiento en el neutro del reactor (siempre que la configuración del sistema lo permita) disipa la energía oscilatoria de secuencia cero, rompiendo la condición de resonancia.

Aplicación Práctica y Análisis mediante Vexten Suite

En el entorno de ingeniería avanzada de Vexten Academy, la validación de los estudios de ferroresonancia y saturación se ejecuta utilizando módulos integrados que combinan los estándares internacionales IEC y IEEE. El flujo de trabajo computacional en la plataforma Vexten Suite abarca las siguientes etapas normativas:

Cálculo de Cortocircuito y Niveles de Contribución según IEC 60909 / IEEE 141

El primer paso consiste en determinar la potencia de cortocircuito de la red en el punto de interconexión del reactor para evaluar la relación de rigidez del sistema ( $SCR$ - Short Circuit Ratio):

SCR = \frac{S_{ sc\_3ph }}{Q_{ reactor }}

Si el $SCR$ es bajo ( $SCR < 3$ ), la red se considera débil, lo que incrementa exponencialmente la susceptibilidad a sobretensiones por ferroresonancia. El motor de cálculo de Vexten Suite procesa matrices de impedancia nodal $Zbus$ para extraer las componentes simétricas exactas durante la apertura monopolar del interruptor.

Derotación Armónica y Dimensionamiento de Cables según IEC 60287 / NEC 310

La presencia de corrientes armónicas elevadas inducidas por la saturación magnética del reactor altera las pérdidas por efecto Joule y por corrientes parásitas en los cables de potencia asociados a la subestación. El factor de derating por armónicos ( $FHD$ ) se calcula en Vexten Suite aplicando los coeficientes especificados en IEC 60287:

I_{ max\_cable } = I_{ base } \cdot \sqrt{\frac{1}{\sum_{h=1}^{k} I_h^2 \cdot (1 + Y_h)}}

Donde $I_h$ es la magnitud relativa del armónico $h$ y $Y_h$ representa el factor de pérdidas suplementarias dependiente de la frecuencia para los conductores y pantallas metálicas del cable.

Mitigación de Factor de Potencia y Resonancia Armónica

Para evitar la coincidencia entre la frecuencia de resonancia paralela del sistema de compensación reactiva y los armónicos característicos generados por la saturación, Vexten Suite realiza un barrido de frecuencia (Frequency Scan) calculando la impedancia de entrada Zin(j\omega) en el rango de $0 Hz$ a $3 kHz$ :

Zin(j\omega) = \left( Rsys(\omega) + j\omega Lsys(\omega) \right) \parallel \left( Rreact(\omega) + j\omega L_{)(\lambda)}(\omega) \right) \parallel \left( \frac{1}{j\omega Ceq} \end{right)

Los resultados de este barrido permiten a los ingenieros de diseño ajustar los parámetros de sintonización de los filtros de armónicos o modificar la inductancia nominal de los reactores de derivación para garantizar que los polos de impedancia queden estrictamente alejados de las frecuencias armónicas dominantes (2^\circ, 3^\circ, 5^\circ, 7^\circ), asegurando así la estabilidad operativa total del sistema de transmisión EHVAC/UHVAC bajo cualquier contingencia de conmutación monofásica.

Conclusiones y Recomendaciones de Ingeniería

La inestabilidad por ferroresonancia y la saturación magnética en reactores de derivación durante maniobras de conmutación monofásica representan un desafío crítico en la ingeniería de sistemas eléctricos de alta potencia. El análisis riguroso mediante modelos no lineales de inductancia dependiente del flujo magnético, combinado con la aplicación estricta de normativas internacionales (IEC 60289, IEC 60071, IEEE Std 1036 y IEEE Std 519), es indispensable para prevenir fallas catastróficas en los activos de la subestación. La utilización de herramientas de simulación avanzada como la plataforma Vexten Suite garantiza que las especificaciones de diseño, los factores de derating térmico y los estudios de barrido de frecuencia mitiguen eficazmente el riesgo de sobretensiones temporales y distorsión armónica extrema, consolidando la confiabilidad y seguridad operativa de las redes de transmisión del futuro.