Ferrorresonancia en Estaciones Megawatt Charging: Límites IEEE 519-2022 y Verificación I²t

¿Por qué la ferrorresonancia en transformadores MCS dispara la distorsión armónica y destruye cables por I²t? Desliza este dossier técnico para evaluar los lími

Ing. Francisco Ramírez

Fenomenología de la Ferrorresonancia en Sistemas de Media Tensión y Transformadores MCS

La ferrorresonancia es un fenómeno dinámico no lineal de resonancia que se presenta en sistemas de potencia eléctrica cuando un circuito que contiene una inductancia no lineal saturable (típicamente el núcleo ferromagnético de un transformador de subestación conectado por cable, o transformador MCS) y una capacitancia (asociada a cables de media tensión, líneas de transmisión o condensadores de acoplamiento) es excitado por una fuente de tensión alterna. A diferencia de la resonancia lineal, donde la reactancia inductiva es constante y la resonancia ocurre a una única frecuencia específica:
ω0=1LC\omega_0 = \frac{1}{\sqrt{L C}}
la ferrorresonancia se caracteriza por poseer múltiples estados estables de operación para un mismo conjunto de parámetros del sistema. Esto se debe a que la inductancia del transformador varía drásticamente en función del flujo magnético, transitando rápidamente de una inductancia lineal de alta impedancia en estado no saturado a una inductancia de muy baja impedancia en estado de saturación profunda. La relación no lineal entre el flujo magnético enlazado Ψ\Psi y la corriente de magnetización imi_m en el núcleo ferromagnético se puede modelar analíticamente mediante una función polinómica impar de orden superior:
im(Ψ)=aΨ+bΨ2n+1i_m(\Psi) = a \cdot \Psi + b \cdot \Psi^{2n+1}
Donde aa y bb son coeficientes característicos del material del núcleo y del diseño magnético, y nn es un número entero positivo (típicamente n2n \ge 2). Cuando el transformador es sometido a una sobretensión o a una maniobra de conmutación asíncrona, el flujo magnético supera el codo de saturación. En este punto, la inductancia incremental o tangente:
LΔ(Ψ)=dΨdimL_{\Delta}(\Psi) = \frac{d\Psi}{di_m}
se reduce en varios órdenes de magnitud. Esta reducción abrupta de la inductancia permite que la reactancia inductiva del transformador se iguale con la reactancia capacitiva del cable de alimentación a la frecuencia industrial o a sus armónicos, desencadenando el estado ferrorresonante. El comportamiento dinámico de este circuito no lineal se describe mediante la ecuación diferencial de Duffing forzada con amortiguamiento:
d2Ψdt2+γdΨdt+αΨ+βΨ3=Fcos(ωt+θ)\frac{d^2 \Psi}{dt^2} + \gamma \frac{d\Psi}{dt} + \alpha \Psi + \beta \Psi^3 = F \cos(\omega t + \theta)
Donde γ\gamma representa las pérdidas resistivas (amortiguamiento), α\alpha y β\beta definen la rigidez no lineal del circuito magnético, y FF es la magnitud de la fuerza electromotriz externa. Dependiendo de las condiciones iniciales de carga, el ángulo de fase de la maniobra θ\theta, y la magnitud de la tensión, el sistema puede converger a diferentes atractores en el espacio de fase, tales como:
  • Modo Fundamental: Las variables de estado oscilan con el mismo período T=2π/ωT = 2\pi/\omega de la fuente de excitación, presentando altos niveles de distorsión armónica.
  • Modo Subarmónico: Las oscilaciones se establecen a una fracción de la frecuencia fundamental (típicamente T/3T/3 o T/5T/5), lo que genera elevadas corrientes de magnetización de larga duración.
  • Modo de Cuasi-periodicidad: El sistema presenta oscilaciones cuyas frecuencias no son múltiplos enteros de la frecuencia de excitación, manifestándose como un espectro continuo con picos discretos.
  • Modo Caótico: Comportamiento aperiódico altamente sensible a las condiciones iniciales, donde la trayectoria en el espacio de fase describe un atractor extraño, provocando sobretensiones y sobrecorrientes erráticas e impredecibles.

Transición del Núcleo Ferromagnético al Estado de Saturación Profunda

Durante la operación normal de un transformador MCS (Medium-Voltage Cable-Connected Substation), el núcleo opera en la región lineal de la curva de histéresis, donde la permeabilidad magnética relativa μr\mu_r es extremadamente alta (10,00010,000 a 100,000100,000). Bajo estas condiciones, la corriente de magnetización necesaria para mantener el flujo es despreciable (menor al 1% de la corriente nominal). Sin embargo, si se produce un evento de conmutación unipolar (por ejemplo, la apertura de una sola fase mediante un fusible o un interruptor monopolar), la corriente de retorno se ve obligada a circular a través de las capacitancias de acoplamiento de los cables de media tensión. Esto induce una sobretensión temporal que desplaza el punto de operación magnética más allá de la rodilla de saturación (B>1.7TB > 1.7 \, T para aceros de grano orientado). Al saturarse el núcleo, la permeabilidad μr\mu_r cae hacia la unidad (μr1\mu_r \approx 1), asemejándose al comportamiento del aire. La inductancia efectiva disminuye drásticamente, permitiendo que la capacitancia del cable cargue y descargue energía directamente sobre la bobina del transformador, lo que perpetúa la oscilación ferrorresonante sin necesidad de una fuente externa de gran potencia.

Análisis de Circuitos y Modelado Matemático de Ferrorresonancia Serie y Paralelo

La topología del circuito determina la naturaleza de la ferrorresonancia y los niveles de estrés dieléctrico y térmico resultantes en la instalación. Es crítico diferenciar rigurosamente entre los acoplamientos en serie y en paralelo.

Ferrorresonancia Serie en Transformadores MCS

La ferrorresonancia serie ocurre típicamente cuando una capacitancia se encuentra conectada en serie con la inductancia no lineal de magnetización del transformador. Este escenario es común en subestaciones MCS alimentadas por cables subterráneos blindados de gran longitud cuando se realiza una maniobra de energización o desenergización fase por fase, o cuando ocurre la fusión de un fusible en una de las fases del circuito primario. Considerando un transformador trifásico con conexión de primario en estrella flotante (neutro aislado) o en delta, alimentado a través de un cable de longitud LcL_c con capacitancia de secuencia cero C0C_0 y capacitancia de secuencia positiva C1C_1. Si la fase A se abre, la inductancia de magnetización de dicha fase queda conectada en serie con la combinación capacitiva equivalente del cable. El circuito equivalente de Thévenin visto desde los terminales de la inductancia saturable de la fase abierta presenta una capacitancia equivalente:
Ceq=3C0+CdCeq = 3 C_0 + C_d
Donde CdC_d representa la capacitancia de acoplamiento entre conductores del cable. Las ecuaciones de estado que gobiernan este sistema en el dominio del tiempo, considerando el acoplamiento magnético y despreciando las pérdidas para el análisis de peor caso, son:
dΨAdt=vC(t)vs(t)\frac{d\Psi_A}{dt} = v_C(t) - v_s(t)
dvCdt=1Ceqim(ΨA)\frac{dv_C}{dt} = -\frac{1}{Ceq} i_m(\Psi_A)
Sustituyendo la relación polinómica de la corriente de magnetización, obtenemos el sistema acoplado no lineal:
d2ΨAdt2+1Ceq(aΨA+bΨA2n+1)=dvs(t)dt\frac{d^2 \Psi_A}{dt^2} + \frac{1}{Ceq} \left( a \cdot \Psi_A + b \cdot \Psi_A^{2n+1} \right) = \frac{dv_s(t)}{dt}
Este sistema presenta una inestabilidad de bifurcación de tipo "silla-nodo" (saddle-node bifurcation). Cuando la tensión de la fuente supera un valor crítico VcritVcrit, el punto de operación salta bruscamente de la rama de baja tensión a la rama de alta tensión ferrorresonante, donde la tensión en los terminales del transformador y del cable puede alcanzar valores de hasta 3.0p.u.3.0 p.u. a 4.5p.u.4.5 p.u..

Ferrorresonancia Paralelo en Transformadores MCS

La ferrorresonancia paralelo se manifiesta cuando la inductancia de magnetización no lineal se encuentra en paralelo con la capacitancia del sistema. Este fenómeno es característico de sistemas con neutro aislado o puestos a tierra a través de una alta impedancia (bobina de Petersen o resistencia de alto valor) durante la ocurrencia y posterior despeje de una falla monofásica a tierra. También se presenta comúnmente en transformadores de tensión (TTVV) conectados a barras acopladas mediante interruptores abiertos que poseen condensadores de gradiente de tensión en sus cámaras de extinción. En este caso, la capacitancia del cable a tierra CgC_g y la inductancia de magnetización LmL_m forman un circuito resonante paralelo excitado por una corriente de secuencia cero o por la acoplada capacitivamente a través de los contactos del interruptor. La ecuación diferencial que describe la tensión de nodo v(t)v(t) en el circuito paralelo es:
Cpdv(t)dt+1Rpv(t)+tim(Ψ)dτ=is(t)C_p \frac{dv(t)}{dt} + \frac{1}{R_p} v(t) + \int_{-\infty}^{t} i_m(\Psi) d\tau = i_s(t)
Donde CpC_p es la capacitancia equivalente en paralelo, RpR_p representa las pérdidas en el hierro del núcleo y las cargas conectadas al secundario del transformador, e is(t)i_s(t) es la corriente de inyección del sistema. Derivando respecto al tiempo y expresando en términos del flujo Ψ\Psi:
Cpd2Ψdt2+1RpdΨdt+(aΨ+bΨ2n+1)=ddt(is(t)dt)C_p \frac{d^2 \Psi}{dt^2} + \frac{1}{R_p} \frac{d\Psi}{dt} + \left( a \cdot \Psi + b \cdot \Psi^{2n+1} \right) = \frac{d}{dt} \left( \int i_s(t) dt \right)
En la ferrorresonancia paralelo, la corriente total absorbida por el conjunto puede ser relativamente baja, pero la tensión del neutro respecto a tierra (tensión de desplazamiento de neutro) se eleva a valores extremadamente altos y permanentes, destruyendo el aislamiento de los transformadores de medida y de los terminales del cable de MT.

Dinámica de Sobretensión y Degradación Térmica del Aislamiento: Límite Térmico I²t

El principal peligro operacional de la ferrorresonancia radica en la combinación sinérgica de sobretensiones transitorias y permanentes de alta frecuencia con un incremento masivo en el contenido armónico de la corriente. Esto somete tanto al aislamiento sólido del transformador como al aislamiento polimérico del cable (XLPE o EPR) a un estrés electrotérmico extremo.

Límite Térmico de Cables de Media Tensión (I²t)

Bajo condiciones de cortocircuito o transitorios severos de corta duración, el calentamiento de los conductores de los cables de media tensión se considera un proceso adiabático. Esto significa que no existe transferencia de calor significativa desde el conductor hacia los elementos externos (pantalla, cubierta, terreno) debido a la rapidez del evento (t<5st < 5 s). El límite térmico del conductor viene definido por la integral de Joule o energía específica admisible:
I2t=K2S2I^2 t = K^2 S^2
Donde:
  • II es la corriente eficaz (RMS) del transitorio o cortocircuito (AA).
  • tt es la duración del evento (ss).
  • SS es la sección transversal del conductor (mm2mm ^2).
  • KK es la constante de material térmico del conductor, calculada según la norma IEC 60287 e IEC 60949:
K=Qc(β+20)ρ20ln(θf+βθi+β)K = \sqrt{\frac{Q_c (\beta + 20)}{\rho_{20}} \ln\left( \frac{\theta_f + \beta}{\theta_i + \beta} \right)}
Donde QcQ_c es la capacidad calorífica volumétrica del material conductor (J/Km3J/K \cdot m ^3), β\beta es el recíproco del coeficiente de temperatura de la resistencia a 0C0 \, ^\circ C (234.5234.5 para cobre, 228228 para aluminio), ρ20\rho_{20} es la resistividad del conductor a 20C20 \, ^\circ C, θi\theta_i es la temperatura inicial de operación del conductor (típicamente 90C90 \, ^\circ C para XLPE bajo carga nominal), y θf\theta_f es la temperatura límite de cortocircuito admisible (250C250 \, ^\circ C para XLPE). Sin embargo, la ferrorresonancia es un fenómeno de larga duración (puede sostenerse durante minutos u horas si las protecciones no actúan). En este régimen no adiabático, la ecuación diferencial de balance térmico del cable debe incorporar la resistencia térmica del aislamiento T1T_1, de la cubierta T3T_3 y del medio circundante T4T_4:
Cthdθ(t)dt=I2Rca(θ)θ(t)θambRthCth \frac{d\theta(t)}{dt} = I^2 Rca(\theta) - \frac{\theta(t) - \theta_{amb}}{\sum Rth}
Donde CthCth es la capacidad térmica equivalente del cable por unidad de longitud, Rca(θ)Rca(\theta) es la resistencia en corriente alterna del conductor a la temperatura θ\theta (altamente influenciada por los efectos skin y de proximidad debidos a las altas frecuencias de la ferrorresonancia), y Rth\sum Rth es la resistencia térmica total del cable al ambiente. Durante la ferrorresonancia, el espectro de corriente presenta una distorsión armónica total de corriente (THDITHD _I) que supera con frecuencia el 100%, con componentes de alta frecuencia que incrementan la resistencia de CA debido al efecto skin:
Rca(f)=Rcc(1+xs(f)+xp(f))Rca(f) = Rcc \left( 1 + x_s(f) + x_p(f) \right)
Esto genera pérdidas por efecto Joule muy superiores a las nominales, acelerando la degradación térmica del aislamiento del cable según la ecuación de Arrhenius para el envejecimiento térmico:
Lvida=AeBθ(t)+273.15Lvida = A \cdot e^{\frac{B}{\theta(t) + 273.15}}
La exposición continua a temperaturas por encima del límite de diseño de emergencia (130C130 \, ^\circ C para XLPE) provoca la pérdida de las propiedades mecánicas del polímero, la formación de "arborescencias térmicas" (thermal treeing) y, finalmente, la perforación dieléctrica del cable.

Estrés Térmico y Dieléctrico en Devanados de Transformadores

En el transformador MCS, la saturación profunda del núcleo desvía el flujo magnético fuera del circuito magnético principal hacia las estructuras metálicas de soporte, los pernos de apriete y las paredes del tanque. Este flujo de dispersión induce corrientes de Foucault masivas en estas partes estructurales, provocando un calentamiento localizado extremo (hot-spots). Adicionalmente, las pérdidas por histéresis y corrientes de Foucault en el propio núcleo se incrementan exponencialmente con la inducción magnética BB y la frecuencia de las oscilaciones:
Pnucleo=ηfBmax1.6+σef2Bmax2d2Pnucleo = \eta \cdot f \cdot Bmax^{1.6} + \sigma_e \cdot f^2 \cdot Bmax^2 \cdot d^2
Donde η\eta es el coeficiente de Steinmetz, σe\sigma_e es la conductividad eléctrica del acero al silicio, y dd es el espesor de las laminaciones del núcleo. Bajo condiciones ferrorresonantes, la elevación de temperatura del punto más caliente del devanado (θH\theta_H) se calcula dinámicamente según la norma IEEE C57.91 mediante:
θH(t)=θA+ΔθTO(t)+ΔθH(t)\theta_H(t) = \theta_A + \Delta \theta_{\text{TO}}(t) + \Delta \theta_{\text{H}}(t)
Donde θA\theta_A es la temperatura ambiente, ΔθTO\Delta \theta_{\text{TO}} es la elevación de temperatura del aceite en la parte superior, y ΔθH\Delta \theta_{\text{H}} es la elevación de temperatura del punto más caliente sobre el aceite. La combinación de esta alta temperatura con el estrés por sobretensión de alta frecuencia destruye el papel aislante impregnado en aceite (pérdida del grado de polimerización), liberando gases inflamables e induciendo fallas entre espiras.

Tabla Comparativa de Parámetros Eléctricos, Límites Normativos y Consecuencias Dieléctricas

La siguiente tabla detalla los parámetros críticos de operación, los límites normativos internacionales y el impacto físico-químico sobre los activos de la subestación MCS cuando se superan dichos umbrales durante un evento ferrorresonante:
Parámetro / Escenario Límites Normativos (IEEE/IEC) Condición Crítica de Activación Consecuencias Dieléctricas y Operacionales
Sobretensión Temporal (TOV) IEC 60071-1: UmU_m (máxima tensión del sistema).
IEEE C62.11: Curva TOV de descargadores de sobretensión (típicamente <1.4p.u.< 1.4 p.u. para t>10st > 10 s).
Ferrorresonancia serie por apertura de fase o conmutación unipolar en cables de gran longitud (>1.5km> 1.5 km). Perforación del aislamiento sólido-líquido en transformadores; falla destructiva y explosión de descargadores de sobretensiones (MOVs) por sobrecarga térmica.
Límite Térmico de Cables (I2tI^2t) IEC 60287 / IEC 60949: Límite adiabático I2t=K2S2I^2t = K^2 S^2. θmax=250C\theta_{max} = 250 \, ^\circ C (cortocircuito) y 130C130 \, ^\circ C (sobrecarga de emergencia). Ferrorresonancia sostenida de modo fundamental o subarmónico con corrientes RMS superiores a 2.5p.u.2.5 p.u. Fusión de la pantalla metálica del cable, degradación irreversible del XLPE, deformación plástica del conductor y colapso del aislamiento de las terminaciones.
Factor de Pérdidas por Armónicos en Transformadores IEEE C57.110: Factor de reducción de capacidad (Derating Factor) para corrientes no sinusoidales:
FHL=h=1maxh2Ih2h=1maxIh2FHL = \frac{\sum_{h=1}^{max} h^2 I_h^2}{\sum_{h=1}^{max} I_h^2}.
Ferrorresonancia caótica o subarmónica con alto contenido de armónicos de orden impar (3º, 5º, 7º, 9º). Calentamiento severo de los devanados por pérdidas adicionales por corrientes parásitas (eddy currents); carbonización del papel aislante Kraft y aceleración del envejecimiento (pérdida de DP).
Tensión de Desplazamiento del Neutro IEEE C57.105: Conexiones de neutro en transformadores trifásicos. Tensión de neutro máxima recomendada <0.1p.u.< 0.1 p.u. en régimen continuo. Ferrorresonancia paralelo en sistemas con neutro aislado o puestos a tierra por alta impedancia durante fallas a tierra monofásicas. Saturación extrema de transformadores de tensión (TTVV) de barra, quema de fusibles primarios, fallas de fase a tierra secundarias por acoplamiento capacitivo.
Capacidad de Absorción de Energía (Descargadores) IEC 60099-4: Clasificación de energía de descarga de línea (Clases de transferencia de carga de 11 a 55). Energía máxima admisible en kJ/kVkJ/kV de MCOV. Inyección de energía continua por sobretensiones ferrorresonantes cuasi-periódicas que superan la tensión de cebado del pararrayos. Embalamiento térmico (thermal runaway) de las pastillas de ZnOZnO; cortocircuito interno franco en el descargador y falla catastrófica con proyección de fragmentos.

Análisis Forense de Fallas de Ingeniería

El análisis forense de fallas en subestaciones MCS que han experimentado eventos ferrorresonantes requiere una investigación multidisciplinaria que abarque metalurgia, química de aislamientos y modelado electromagnético transitorio. Los modos de falla se manifiestan de forma característica en los diferentes componentes del sistema.

Fallas en el Núcleo y Devanados del Transformador

Cuando un transformador MCS es sometido a ferrorresonancia, la firma forense interna es inconfundible y difiere significativamente de una falla por cortocircuito simétrico estándar:
  • Evidencia Física en el Núcleo: Se observan patrones de decoloración térmica localizados en las juntas de las esquinas del núcleo (juntas a inglete o mitred joints) y en los pernos de sujeción (core bolts). Esto es causado por el flujo magnético que sale del plano de las laminaciones de acero magnético, induciendo corrientes de Foucault circulares de gran magnitud directamente sobre la superficie de las chapas. El barniz aislante entre laminaciones se carboniza, lo que genera cortocircuitos entre chapas y un incremento permanente de las pérdidas en vacío.
  • Análisis de Gases Disueltos (DGA - IEC 60599): Las muestras de aceite dieléctrico extraídas inmediatamente después del evento muestran concentraciones críticas de gases clave. La ferrorresonancia se caracteriza por una firma de "Falla Térmica de Alta Temperatura" (T>700CT > 700 \, ^\circ C) combinada con descargas parciales de alta energía. Los gases predominantes son el Etileno (C2H4C _2 H _4) y el Acetileno (C2H2C _2 H _2), este último indicativo de arcos eléctricos locales entre espiras o entre el núcleo y la cuba debido a la pérdida de rigidez dieléctrica del aceite sobrecalentado. La relación de gases según el método de Rogers o el Triángulo de Duval sitúa inequívocamente el evento en la zona de falla térmica severa (T3T3) combinada con descargas de alta energía (D2D2).
  • Destrucción del Aislamiento de Devanados: La contracción térmica y la carbonización del papel aislante Kraft provocan una pérdida de la fuerza de apriete axial de los devanados. Bajo el estrés mecánico de las fuerzas de cortocircuito inducidas por las corrientes ferrorresonantes, las bobinas sufren deformaciones geométricas, pandeo radial (radial buckling) y, finalmente, un cortocircuito franco entre espiras adyacentes de la sección de entrada de la bobina de alta tensión.

Fallas en Cables de Media Tensión y Terminales

Los cables de media tensión que conectan la subestación actúan como el elemento capacitivo del circuito ferrorresonante. Sus modos de falla forenses incluyen:
  • Degradación de la Pantalla Metálica por Corrientes Homopolares: Durante la ferrorresonancia serie o paralelo con desplazamiento de neutro, circulan corrientes de alta frecuencia por la pantalla metálica (hilos de cobre o cinta de aluminio). Dado que estas pantallas están dimensionadas únicamente para conducir corrientes de falla a tierra durante fracciones de segundo (por ejemplo, 1s1 s o 3s3 s), el paso continuo de corrientes ferrorresonantes de varias decenas de amperios durante minutos provoca la fusión de los hilos de la pantalla, quemando la cubierta exterior de PVC o PE e interrumpiendo la continuidad de tierra del cable.
  • Perforación por Arborescencias Eléctricas (Electrical Treeing): El análisis microscópico de muestras de XLPE falladas revela la presencia de estructuras ramificadas carbonizadas que se originan en los puntos de concentración de estrés eléctrico (como protuberancias en la pantalla semiconductora o microcavidades en el aislamiento). Estas arborescencias se propagan aceleradamente bajo las sobretensiones de alta frecuencia características de la ferrorresonancia cuasi-periódica hasta perforar completamente el espesor del aislamiento.

Fallas en Aparamenta de Maniobra y Descargadores de Sobretensión

Los descargadores de sobretensión (pararrayos de óxido de zinc, ZnOZnO) son los elementos más vulnerables ante la ferrorresonancia:
  • Mecanismo de Falla del Descargador: Los pararrayos de ZnOZnO están diseñados para limitar sobretensiones transitorias de muy corta duración (microsegundos para descargas atmosféricas, milisegundos para sobretensiones de maniobra). No están diseñados para soportar sobretensiones temporales (TOV) continuas. Durante la ferrorresonancia, la tensión del sistema supera repetidamente la tensión de operación continua admisible (MCOVMCOV) del descargador, forzando a las pastillas de varistor a entrar en conducción profunda durante cada semiciclo. La energía disipada supera la capacidad de absorción térmica del varistor (medida en kJ/kVkJ/kV), provocando un embalamiento térmico: a medida que aumenta la temperatura del ZnOZnO, su resistencia disminuye, lo que incrementa la corriente de fuga y eleva aún más la temperatura, culminando en la destrucción explosiva del cuerpo del pararrayos.

Estrategias Prácticas de Diseño, Mitigación y Control

La mitigación de la ferrorresonancia debe abordarse desde la fase de diseño de la ingeniería de detalle de la subestación MCS, aplicando criterios de robustez topológica y sistemas de amortiguamiento activo o pasivo.

Eliminación de Dispositivos de Maniobra Monopolar

La medida de diseño más efectiva para prevenir la ferrorresonancia serie es la eliminación absoluta de fusibles unipolares y seccionadores de operación monofásica en el lado de media tensión del transformador. En su lugar, se debe prescribir obligatoriamente el uso de interruptores trifásicos de operación tripolar simultánea (gang-operated switchgear). Esto asegura que todas las fases se abran o cierren con una asincronía máxima entre polos inferior a 2ms2 ms, evitando la existencia de estados operativos donde una o dos fases queden abiertas mientras las restantes permanecen energizadas a través de la capacitancia de los cables.

Configuración del Neutro del Transformador

La conexión del neutro del devanado primario del transformador juega un papel decisivo en la formación del circuito ferrorresonante:
  • Neutro Sólidamente a Tierra: Conectar el neutro del primario sólidamente a tierra elimina la posibilidad de ferrorresonancia serie al proporcionar un camino de retorno de baja impedancia para las corrientes de secuencia cero, impidiendo que la capacitancia de los cables quede conectada en serie con la inductancia de magnetización. Sin embargo, esto incrementa la magnitud de las corrientes de cortocircuito monofásicas a tierra.
  • Uso de Transformadores de Conexión Estrella-Estrella con Terciario en Delta: El devanado terciario en delta estabiliza el neutro y proporciona un camino de baja impedancia para las corrientes de tercer armónico, reduciendo significativamente la impedancia de magnetización equivalente de secuencia cero y dificultando la entrada en saturación del núcleo.

Sistemas de Amortiguamiento Pasivo (Damping Resistors)

Para instalaciones donde no es posible modificar la configuración de puesta a tierra o la aparamenta de maniobra, se deben diseñar e implementar circuitos de amortiguamiento pasivo. Para transformadores de tensión (TTVV) conectados a barras con riesgo de ferrorresonancia paralelo, se utiliza una resistencia de amortiguamiento (RdR_d) conectada en los terminales del devanado secundario en delta abierta (open-delta secondary):
Rd3ωC0R_d \le \frac{3}{\omega C_0}
Donde C0C_0 es la capacitancia equivalente de secuencia cero del sistema de barras y cables visto desde el primario del transformador, referida al secundario mediante la relación de transformación at2a_t^2:
C0,sec=C0at2C_{0,sec} = C_0 \cdot a_t^2
La resistencia disipa la energía acumulada en el circuito resonante, forzando a los polos del sistema en el espacio de fase a colapsar hacia el atractor de operación normal de baja tensión. El valor de potencia nominal de la resistencia de amortiguamiento debe calcularse para soportar la corriente de circulación continua inducida por el desequilibrio normal del sistema y la corriente transitoria durante una falla monofásica a tierra sin sufrir sobrecalentamiento:
PRd(3Vsec,nom)2Rdtfalla3600PRd \ge \frac{\left( 3 V_{sec,nom} \right)^2}{R_d} \cdot \frac{tfalla}{3600}

Implementación de Filtros Sintonizados y Supresores Activos

En subestaciones MCS de gran envergadura con presencia de armónicos significativos, se pueden instalar dispositivos de mitigación activos, tales como el "Supresor Electrónico de Ferrorresonancia" (EFS). Estos equipos monitorean continuamente la tensión de neutro y, al detectar la firma espectral característica de una oscilación ferrorresonante (especialmente subarmónicos de 20Hz20 \, Hz o 16.67Hz16.67 \, Hz en sistemas de 50Hz50 \, Hz), insertan de forma instantánea mediante tiristores de potencia una resistencia de amortiguamiento de muy bajo valor durante unos pocos ciclos, extinguiendo el fenómeno sin introducir pérdidas permanentes en el sistema.

Aplicación Práctica y Simulación con Vexten Suite

El análisis de susceptibilidad a la ferrorresonancia y la verificación del límite térmico I2tI^2t se realizan de manera integrada utilizando la plataforma avanzada de ingeniería eléctrica Vexten Suite. El flujo de trabajo metodológico se divide en tres módulos de simulación interconectados. ``` +-----------------------------------------------------------------+ | VEXTEN SUITE | | | | +---------------------+ +-------------------------------+ | | | MÓDULO 1: | | MÓDULO 2: | | | | Vexten Short- | | Vexten Cable Sizing | | | | Circuit | | & Ampacity | | | | - IEC 60909 | | - IEC 60287 | | | | - S_sc, X_s/R_s | | - I^2t adiabatic limit | | | +----------+----------+ +---------------+---------------+ | | | | | | +----------------+---------------+ | | | | | v | | +-----------------------+ | | | MÓDULO 3: | | | | Vexten Power Factor | | | | & Resonance | | | | - Harmonic scan | | | | - De-tuned filters | | | +-----------------------+ | +-----------------------------------------------------------------+ ```

Módulo 1: Vexten Short-Circuit (IEC 60909 / IEEE 141)

El primer paso consiste en determinar la rigidez dieléctrica y la impedancia equivalente de la red de alimentación en el punto de conexión común (PCC) de la subestación MCS. Utilizando el motor de cálculo de cortocircuito de Vexten Short-Circuit bajo la norma IEC 60909, se modela la red aguas arriba mediante su potencia de cortocircuito trifásica (SscSsc) y la relación X/RX/R:
Zk=cUn3IscZ_k = \frac{c \cdot U_n}{\sqrt{3} Isc}
Donde cc es el factor de tensión según IEC 60909. Una potencia de cortocircuito baja (red débil) implica una impedancia de fuente ZkZ_k elevada. Esta condición incrementa drásticamente la susceptibilidad a la ferrorresonancia serie, ya que la impedancia de la fuente no es capaz de estabilizar la tensión ante variaciones rápidas de la corriente de magnetización del transformador. El módulo calcula la matriz de impedancias de transferencia del sistema, la cual es exportada directamente al solver de transitorios electromagnéticos de Vexten para definir las condiciones de frontera de la simulación dinámica.

Módulo 2: Vexten Cable Sizing & Ampacity (IEC 60287 / NEC 310)

Para evaluar el impacto térmico de las corrientes ferrorresonantes sobre los cables de MT de la subestación, se utiliza el módulo Vexten Cable Sizing & Ampacity. Este módulo implementa las ecuaciones de transferencia de calor en estado estacionario y transitorio de la norma IEC 60287. El usuario introduce la geometría del cable (por ejemplo, cable tripolar XLPE de 15kV15 \, kV, sección de 150mm2150 \, mm ^2 de cobre, con pantalla de hilos de cobre de 16mm216 \, mm ^2), la disposición de instalación (enterrado directamente, en ducto o al aire) y los parámetros del suelo (resistividad térmica del terreno g=1.0Km/Wg = 1.0 \, K \cdot m/W, temperatura del suelo θs=25C\theta_s = 25 \, ^\circ C). Bajo condiciones de ferrorresonancia detectadas en la simulación dinámica, se ingresa el espectro de corriente armónica calculado. El módulo aplica el factor de reducción de capacidad por armónicos (Harmonic Derating Factor) ajustando la resistencia de CA para cada orden armónico hh:
Rca,h=Rcc(1+ysh2+yph2)R_{ca,h} = Rcc \left( 1 + y_s \cdot h^2 + y_p \cdot h^2 \right)
El software evalúa la integral de Joule real del evento ferrorresonante prolongado:
0tduri2(t)dt\int0^{tdur} i^2(t) dt
y la compara con el límite térmico adiabático del cable K2S2K^2 S^2 y con el límite no adiabático de deformación térmica de la cubierta. Si la integral excede el límite permitido, Vexten Suite emite una alerta de violación de diseño y calcula automáticamente el incremento requerido en la sección del conductor o la necesidad de reespecificar la pantalla metálica para evitar la falla del cable.

Módulo 3: Vexten Power Factor & Resonance Mitigation

Este módulo se utiliza para diseñar los sistemas de mitigación de resonancia armónica y ferrorresonancia. Permite realizar un "Harmonic Scan" (barrido de frecuencia de la impedancia del sistema) desde la frecuencia fundamental hasta el armónico 50. Si el barrido revela un punto de resonancia paralelo (pico de alta impedancia) cerca de las frecuencias armónicas típicas de la ferrorresonancia (como el 3º o 5º armónico), el módulo asiste al ingeniero en el diseño de un banco de condensadores con reactancias de desintonización (de-tuned reactor banks). El factor de desintonización pp se calcula mediante:
p=XLXC=(f1fr)2p = \frac{X_L}{X_C} = \left( \frac{f_1}{f_r} \right)^2
Típicamente se selecciona un factor p=7%p = 7\% (frecuencia de sintonía fr=189Hzf_r = 189 \, Hz para sistemas de 50Hz50 \, Hz) o p=5.67%p = 5.67\%. Esto desplaza la frecuencia de resonancia del sistema por debajo del tercer armónico, impidiendo de forma efectiva que las corrientes de magnetización armónicas generadas por la saturación transitoria del transformador MCS encuentren un camino de alta impedancia que pueda iniciar o sostener un fenómeno de ferrorresonancia paralelo. Adicionalmente, el software simula la inserción de las resistencias de amortiguamiento calculadas en el devanado en delta abierta del transformador de tensión, graficando el plano de fase del flujo frente a la corriente de magnetización y demostrando visualmente la convergencia del sistema hacia el atractor estable normal (extinción completa de la ferrorresonancia en menos de 150ms150 \, ms).