ferroresonanciatransformadores de medidaestabilidad de sistemassobretensiones transitoriasdiseño de subestaciones

Ferrorresonancia en Transformadores de Medida: Análisis y Mitigación

Análisis técnico de la ferrorresonancia en transformadores de tensión: causas físicas, modelado de saturación magnética y estrategias de mitigación.

Ing. Francisco Ramírez

Introducción y Fundamentos Físicos de la Ferrorresonancia

La ferrorresonancia representa uno de los fenómenos transitorios y de estado estacionario más complejos, destructivos y difíciles de predecir en los sistemas de potencia eléctrica de media y alta tensión. A diferencia de la resonancia lineal clásica, donde un circuito RLC lineal exhibe una impedancia mínima o máxima a una frecuencia específica de sintonización, la ferrorresonancia es un fenómeno no lineal caracterizado por la presencia de al menos un elemento inductivo con núcleo de hierro ferromagnético que opera en regiones de saturación magnética, acoplado con capacitancias distribuidas o concentradas del sistema (típicamente capacitancias de cables subterráneos, líneas aéreas largas, condensadores de acoplamiento o interruptores abiertos).

El comportamiento no lineal del núcleo magnético se modela matemáticamente mediante la relación constitutiva entre el flujo magnético enlazado \lambda y la corriente de excitación ii, la cual deja de ser una función lineal y se describe mediante una curva de histéresis y saturación que puede aproximarse analíticamente mediante funciones potenciales o polinómicas de alto orden:

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

Donde L0L_0 representa la inductancia en la región lineal (no saturada) y los coeficientes aka_k determinan la curvatura y la severidad de la no linealidad en la zona de saturación profunda. Cuando el núcleo del transformador de medida (ya sea de tensión inductivo - IT, o transformador de potencial capacitivo - CVT) entra en saturación, su inductancia equivalente experimenta una reducción drástica, cayendo a valores cercanos a la inductancia de dispersión del devanado. Este colapso abrupto en la inductancia transforma radicalmente la impedancia del circuito resonante, desplazando la frecuencia natural de oscilación del sistema hacia valores subarmónicos (f0<fsistemaf_0 < fsistema), armónicos (f0>fsistemaf_0 > fsistema) o cuasi-periódicos.

En los transformadores de medida, este fenómeno se desencadena típicamente bajo condiciones de operación anormales pero altamente probables en redes eléctricas modernas: aperturas monofásicas de interruptores en sistemas operando con neutro aislado o compensado por bobina de Petersen, fallas de tierra intermitentes, energización en vacío de transformadores a través de líneas con capacitancia significativa, o resonancias sostenidas tras la interrupción de corrientes de cortocircuito capacitivas.

Clasificación y Modos Dinámicos de Oscilación en Transformadores

La naturaleza altamente no lineal del circuito ferrorresonante implica la existencia de múltiples estados estables posibles para los mismos parámetros circuitales (fenómeno de ferro-multiestabilidad). La transición de un estado operativo normal a un régimen ferrorresonante destructivo depende críticamente de las condiciones iniciales, la magnitud de la tensión de excitación, la capacitancia serie-paralelo y las pérdidas totales del sistema. Los modos de oscilación se clasifican formalmente en cuatro categorías fundamentales según su contenido espectral y firma transitoria:

Modo Fundamental (Resonancia a la Frecuencia del Sistema)

Caracterizado por oscilaciones a la frecuencia fundamental de la red (5050 o 60Hz60 Hz). La corriente que circula por el devanado del transformador contiene una componente fundamental amplificada masivamente (hasta 1010 a 2020 veces la corriente nominal de excitación) y una distorsión armónica moderada. Las sobretensiones asociadas pueden alcanzar valores de 2.02.0 a 3.53.5 por unidad (p.u.), provocando la ruptura térmica progresiva del aislamiento dieléctrico debido al incremento exponencial de las pérdidas en el núcleo por histéresis y corrientes parásitas de Foucault.

Modo Subarmónico (Oscilaciones a Fracciones de la Frecuencia del Sistema)

Las oscilaciones ocurren a frecuencias submúltiplas exactas de la frecuencia fundamental, típicamente f0/2f_0 / 2, f0/3f_0 / 3 o f0/nf_0 / n. Este modo es extremadamente peligroso y desconcertante para los operadores, ya que genera ondas de corriente y tensión de aspecto caótico en el oscilograma. Las tensiones inducidas pueden superar los 4.0p.u.4.0 p.u., y las corrientes alternadas de baja frecuencia fuerzan al núcleo a entrar en saturación profunda bidireccional en cada ciclo extendido, generando esfuerzos mecánicos severos en los bobinajes debido a las fuerzas electrodinámicas de baja frecuencia.

Modo Armónico (Oscilaciones a Múltiplos Enteros de la Frecuencia del Sistema)

Dominado por frecuencias superiores a la fundamental (2f0,3f0,5f02f_0, 3f_0, 5f_0). Se manifiesta frecuentemente durante procesos de energización de transformadores en vacío conectados a redes con alta capacitancia distribuida. El espectro de frecuencia muestra picos prominentes en armónicos pares e impares, generando sobretensiones transitorias de corta duración pero de altísimo gradiente dv/dt, lo que estresa severamente el aislamiento entre espiras de los devanados del transformador de medida.

Modo Cuasi-Periódico (Caótico)

Representa un régimen dinámico no lineal donde el voltaje y la corriente contienen múltiples frecuencias incommensurables (no relacionadas por números enteros). El sistema oscila de manera aperiódica, saltando erráticamente entre diferentes niveles de energía magnética. Este modo se visualiza en el plano de fase mediante atractores extraños y representa el escenario más severo para la integridad estructural y térmica de los transformadores de potencial, conduciendo a la falla catastrófica en cuestión de minutos.

Análisis del Circuito Equivalente y Ecuaciones Diferenciales No Lineales

Para modelar rigurosamente el fenómeno, considérese un sistema monofásico equivalente compuesto por una fuente de tensión alterna ideal e(t) = E_m \sin(\omega t), una resistencia serie RR que representa las pérdidas en los conductores y la red, una inductancia lineal de línea LsL_s, una capacitancia shunt CC (representando la capacitancia de la línea o cable hacia tierra), y un transformador de tensión inductivo modelado por una rama no lineal que relaciona el flujo \lambda(i) y una resistencia de pérdidas en el núcleo RcR_c.

Aplicando las leyes de Kirchhoff para las mallas y nodos del circuito equivalente, se obtiene el sistema de ecuaciones integro-diferenciales no lineales que rige la dinámica del sistema:

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}

Donde iLi_L es la corriente por la inductancia serie, vcv_c es la tensión en los bornes del condensador y del transformador, e i_\mu(v_c) es la característica de magnetización no lineal del núcleo ferromagnético. Expresando la relación inversa del flujo magnético en función de la tensión v_c = \frac{d\lambda}{dt}, se puede reducir el sistema a una ecuación diferencial ordinaria de segundo orden no lineal:

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)

Donde la función no lineal f(\lambda) representa la corriente de magnetización en función del flujo concatenado. La solución analítica de esta ecuación es imposible mediante métodos clásicos debido al término no lineal, requiriendo el uso de métodos numéricos avanzados en el domínio del tiempo (como el algoritmo de Runge-Kutta de orden superior) o métodos analíticos aproximados como el balance armónico y la teoría de perturbaciones.

El análisis de estabilidad de Lyapunov y el estudio de los exponentes de Lyapunov permiten determinar las condiciones críticas de bifurcación (bifurcaciones de Hopf y de silla-nodo) que precipitan al sistema desde un estado operativo estable hacia la ferrorresonancia caótica.

Parametrización de Factores Críticos y Marcos Normativos IEEE/IEC

Los organismos internacionales de normalización (IEEE e IEC) establecen directrices estrictas para la evaluación, diseño y prueba de transformadores de medida frente al riesgo ferrorresonante. La siguiente tabla comparativa detalla los parámetros eléctricos críticos, los umbrales normativos y las consecuencias operativas asociadas:

Parámetro Eléctrico / Condición Límite Normativo (IEC 61869 / IEEE C57.13) Condición Crítica de Falla Consecuencia Operativa y Dieléctrica
Relación Capacitancia / Inductancia (C/LC/L) C<CcrıˊticaC < C_{crítica} (Margen de seguridad > 300%) CC elevada por desconexión de líneas largas o cables soterrados Desplazamiento de la frecuencia de resonancia hacia la banda subarmónica.
Factor de Saturación del Núcleo (ns=Bsat/Bopn_s = Bsat / Bop) n_s \geq 2.5 para operación nominal lineal Operación con V>1.5p.u.V > 1.5 p.u. o remanencia magnética elevada Colapso de la inductancia magnetizante y entrada en régimen de saturación profunda.
Densidad de Corriente de Excitación I_\mu \leq 2\% a 1.0 p.u. ; \leq 10\% a 1.2 p.u. \ Corrientes de magnetización distorsionadas con picos superiores al 500\% Sobrecarga térmica severa en el devanado primario y destrucción del aislamiento sólido.
Sobretensión Temporal Transitoria (TOV) Soportar 1.9p.u.1.9 p.u. durante 30segundos30 segundos (IEC 61869-3) Oscilaciones ferrorresonantes sostenidas superiores a 3.0p.u.3.0 p.u. Descargas parciales descontroladas, perforación del papel impregnado en aceite y explosión.
Resistencia de Amortiguamiento Externa Disipación térmica nominal calculada para régimen continuo de falla Ausencia de resistencia de amortiguamiento en el devanado secundario abierto (broken delta) Imposibilidad de disipar la energía reactiva acumulada, perpetuando la ferrorresonancia.

Ingeniería Forense de Fallas en Transformadores de Medida

La ocurrencia de ferrorresonancia en subestaciones eléctricas deja huellas físicas y metalúrgicas inconfundibles que permiten a los ingenieros forenses reconstruir la secuencia de eventos que condujeron a la falla catastrófica del equipo. El análisis forense abarca desde la inspección visual externa hasta la metalurgia analítica de los devanados.

Deformación Electrodinámica y Térmica de los Devanaos

Bajo condiciones ferrorresonantes, las corrientes de cortocircuito interno y de excitación masiva generan fuerzas de Lorentz extremas entre las espiras concéntricas del transformador de tensión. Los devanados primarios, construidos con conductores de cobre de sección reducida para soportar altas tensiones, sufren compresión radial y esfuerzos axiales que superan el límite elástico del material, provocando cortocircuitos entre espiras contiguas. Térmicamente, el flujo magnético disperso y las corrientes parásitas inducidas en las pantallas electrostáticas y estructuras metálicas internas generan puntos calientes (hot spots) con temperaturas superiores a 400°C400 °C, carbonizando el aislamiento de papel celulósico y descomponiendo el aceite aislante.

Análisis de Gases Dissueltos (DGA) Post-Falla

En transformadores de potencial sumergidos en aceite, el análisis cromatográfico de los gases disueltos en el fluido dieléctrico revela patrones característicos de falla. La ferrorresonancia prolongada genera una descomposición térmica severa del aceite mineral, caracterizada por concentraciones masivas de:

  • Hidrógeno (H2H _2): Indicador de descargas parciales y estrés eléctrico severo.
  • Metano (CH4CH _4) y Etano (C2H6C _2 H _6): Evidencia de sobrecalentamiento local moderado a severo (300°C300 °C a 700°C700 °C).
  • Acetileno (C2H2C _2 H _2): Marcador indiscutible de arcos eléctricos de alta potencia y descargas disruptivas entre espiras o hacia tierra.

Comportamiento de los Elementos de Protección Primaria

Los fusibles limitadores de corriente instalados en el primario de los transformadores de tensión desempeñan un papel crítico y a menudo complejo durante la ferrorresonancia. Debido al contenido armónico elevado y a los valores eficaces de corriente que no alcanzan el umbral de fusión rápida del fusible pero que sobrepasan la capacidad térmica de disipación, los elementos fusibles sufren un envejecimiento prematuro ("fatiga térmica"). Esto provoca la fusión asimétrica de las fases, perpetuando la condición de apertura monofásica que alimenta y sostiene el fenómeno ferrorresonante en las fases sanas restantes.

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

La mitigación efectiva de la ferrorresonancia requiere un enfoque multidisciplinario que combina la selección juiciosa de parámetros constructivos, la implementación de topologías de circuitos de amortiguamiento y la adopción de filosofías operativas avanzadas en el control de interruptores.

Cálculo y Selección de Resistencias de Amortiguamiento (Anti-Ferroresonance Dampers)

La inserción de una resistencia de amortiguamiento en el circuito secundario conectado en triángulo abierto ("broken delta") es la técnica más económica y confiable en transformadores de tensión inductivos. La potencia disipada y la resistencia óptima RdRd se calculan mediante la siguiente formulación analítica basada en la impedancia característica del sistema:

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

Donde VsecVsec es la tensión nominal secundaria, ScritScrit es la potencia aparente crítica de disipación requerida para absorber la energía reactiva oscilante, y LeqLeq y CeqCeq representan la inductancia equivalente de saturación y la capacitancia shunt del sistema, respectivamente. El dimensionamiento térmico de esta resistencia debe garantizar que pueda soportar la corriente de secuencia cero resultante durante un tiempo indefinido sin perder sus propiedades óhmicas.

Criterios de Diseño en Transformadores de Potencial Capacitivos (CVT)

En los CVTs, la presencia del divisor capacitivo (C1,C2C_1, C_2) acoplado al circuito electromagnético de compensación inductiva (LtL_t) crea un circuito intrínsecamente propenso a la ferrorresonancia. Para neutralizar este riesgo, se incorpora obligatoriamente un circuito de supresión de ferrorresonancia (suction circuit o ferro-detector) conectado en el secundario o en un devanado terciario auxiliar del transformador intermedio. Este circuito consta de un elemento no lineal (varistor de óxido de zinc - MOV o un circuito sintonizado RLC) que introduce una alta conductancia eléctrica exclusivamente cuando se detectan tensiones con componentes subarmónicas o sobretensiones superiores a 1.2p.u.1.2 p.u., amortiguando instantáneamente la oscilación sin interferir con la precisión metrológica del CVT bajo operación normal.

Estrategias Operativas de Maniobra (Sequential Tripping & Controlled Switching)

Desde la perspectiva de la automatización de subestaciones, la prevención debe enfocarse en evitar las causas raíces de la ferrorresonancia:

  • Uso de Interruptores Tripolares Sincronizados: La apertura y cierre coordinado de los polos del interruptor mediante controladores de maniobra sincronizada (Controlled Switching) elimina la posibilidad de aperturas monofásicas prolongadas en sistemas aislados.
  • Eliminación de Redes en Vacío Críticas: Evitar la energización de transformadores de medida a través de tramos excesivamente largos de cables subterráneos sin carga.
  • Puesta a Tierra Efectiva del Neutro: Mantener una impedancia de puesta a tierra suficientemente baja en los transformadores de potencia principales para evitar desplazamientos severos del punto neutro que sobreexciten los transformadores de tensión adyacentes.

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

Para ilustrar el rigor analítico exigido en el diseño de ingeniería de alta gama, se presenta un caso de estudio resuelto mediante el motor de cálculo avanzado de la plataforma Vexten Suite, integrando normativas internacionales para cortocircuito, dimensionamiento de cables y mitigación de resonancia armónica.

Definición del Escenario de Red en Vexten Suite

Se analiza una subestación de alta a media tensión (138kV/13.8kV138 kV / 13.8 kV) alimentando un sistema con una potencia de cortocircuito trifásica simétrica de Ssc=2500MVASsc = 2500 MVA. La red de media tensión cuenta con un sistema de cables subterráneos poliméricos de gran longitud que aportan una capacitancia total a tierra por fase de C = 1.2\,\mu F . Se encuentra instalado un banco de transformadores de tensión inductivos con inductancia magnetizante lineal L0=45HL_0 = 45 H y una inductancia en saturación profunda Lsat=0.8HLsat = 0.8 H.

Paso 1: Validación de Cortocircuito y Niveles de Excitación (IEC 60909 / IEEE 141)

Utilizando el módulo de cortocircuito de Vexten Suite, se calculan las componentes simétricas y asimétricas de la corriente de falla y el factor de choque térmico (ktk_t) para verificar la resistencia mecánica de los transformadores de medida ante eventos transitorios:

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

El factor de potencia de cortocircuito de la red (\cos\varphi = R_1 / Z_1) determina el decaimiento de la componente unidireccional de corriente continua, la cual interactúa con el núcleo de hierro y puede inducir un estado inicial de remanencia magnética crítica evaluado por Vexten como \psirem = 0.85\,\psisat.

Paso 2: Dimensionamiento de Cables y Derating Harmónico (IEC 60287 / NEC 310)

Durante un evento ferrorresonante, el espectro de corriente se encuentra severamente contaminado por armónicos de orden bajo y medio. El módulo de dimensionamiento de cables de Vexten Suite calcula el factor de reducción por contenido armónico (FharmFharm) aplicable a los cables de control y potencia asociados a los transformadores de medida, garantizando que la temperatura de operación del aislamiento XLPE no supere los 90°C90 °C nominales:

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

Donde \mu_h es la relación armónica de la corriente de orden hh respecto a la fundamental. Para un espectro típico obtenido en simulación ferrorresonante donde el tercer armónico alcanza el 35\% y el quinto el 20\%, Vexten Suite determina un factor de derating térmico estricto de Fharm=0.74Fharm = 0.74, exigiendo un sobredimensionamiento inmediato de la sección transversal de los conductores secundarios de instrumentación para evitar la degradación térmica del aislamiento.

Paso 3: Análisis de Resonancia y Sintonización de Amortiguamiento

El motor de análisis dinámico no lineal de Vexten Suite ejecuta la integración numérica en el dominio del tiempo del sistema circuital con los parámetros ingresados. Los resultados muestran que sin la red de amortiguamiento, la tensión en los bornes del transformador de tensión alcanza un pico sostenido de 3.42p.u.3.42 p.u. en modo subarmónico de orden 1/31/3, con una corriente de excitación que supera los 45A45 A eficaces.

Aplicando las rutinas de optimización de Vexten Suite, se calcula la resistencia de amortiguamiento óptima para el secundario en triángulo abierto:

Rd,opt=125.4Ω(condisipacioˊnteˊrmicanominalPtherm=1.8kW)R_{d,opt} = 125.4\,\Omega \quad ( con disipación térmica nominal Ptherm = 1.8 kW )

La incorporación virtual de este valor en el modelo de simulación colapsa instantáneamente la oscilación ferrorresonante en un intervalo menor a 2.5ciclos2.5 ciclos de la frecuencia fundamental, devolviendo el sistema a su estado operativo estable, con tensiones dentro de los márgenes metrológicos y protegiendo de manera absoluta la integridad dieléctrica y térmica del transformador de medida.

Conclusiones de Ingeniería

La ferrorresonancia en transformadores de medida constituye un fenómeno multifísica crítico que desafía los métodos tradicionales de análisis lineal. Su comprensión exige un dominio profundo de la teoría de circuitos no lineales, la dinámica de sistemas de potencia y el comportamiento termoductivo de los materiales ferromagnéticos. Mediante la aplicación rigurosa de las normativas internacionales IEC e IEEE, el diseño preventivo basado en amortiguamiento secundario y el uso de herramientas de simulación avanzada como Vexten Suite, los ingenieros electricistas pueden garantizar la confiabilidad, seguridad y continuidad operativa de las instalaciones eléctricas de alta y media tensión frente a este complejo régimen transitorio.