Análisis de Frecuencias Naturales y Transitorios Electromagnéticos en Redes de Distribución

Fundamentos Físico-Matemáticos de los Transitorios Electromagnéticos y la Respuesta en Frecuencia

El análisis riguroso de los transitorios electromagnéticos y las frecuencias naturales en las redes de distribución de energía eléctrica exige una formulación matemática avanzada que trasciende el análisis de régimen permanente fasorial. Las redes de distribución modernas, caracterizadas por una alta penetración de cables subterráneos de media tensión, transformadores con núcleos ferromagnéticos no lineales y cargas electrónicas descentralizadas, actúan como sistemas de parámetros distribuidos y concentrados sujetos a excitaciones impulsivas y armónicas severas.

Para modelar la propagación de ondas electromagnéticas a lo largo de los conductores de la red, se deben emplear las ecuaciones diferenciales parciales acopladas de los telegraphers (ecuaciones de los telegrafistas). Considerando una línea multiconductora simétrica o asimétrica de longitud infinitesimal $dx$ , el comportamiento instantáneo de las tensiones v(x,t) y corrientes i(x,t) se define mediante el siguiente sistema matricial:

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

Donde $\mathbf{R}$ , $\mathbf{L}$ , $\mathbf{C}$ y $\mathbf{G}$ representan las matrices por unidad de longitud de resistencia, inductancia, capacitancia y conductancia transversal (fuga dieléctrica), respectivamente. Al transformar estas ecuaciones al dominio de la frecuencia mediante la transformada de Laplace o Fourier, se obtiene la ecuación de propagación modal:

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

Donde Z(s) = R + sL, Y(s) = G + sC, y P(s) es la matriz de propagación. Las autovalores de la matriz producto Z(s)Y(s) determinan las constantes de propagación modal \gamma_k(\omega) = \alpha_k(\omega) + j\beta_k(\omega), las cuales gobiernan la atenuación y la velocidad de fase de las ondas viajeras. Las frecuencias naturales del sistema (\omega_n) corresponden a los polos de la función de transferencia impedancia de entrada Zin(s) vista desde cualquier nodo de la red de distribución, resolviendo la condición característica:

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

Donde M_{end}(s) representa la matriz de reflexión y transmisión en las terminales de la red, acoplando las impedancias de los transformadores, cargas y dispositivos de protección contra sobretensiones (SPD).

Clasificación y Origen de Transitorios Electromagnéticos en Redes de Media y Baja Tensión

Los transitorios electromagnéticos en las redes de distribución se originan a partir de perturbaciones operativas y atmosféricas que inyectan espectros de frecuencia de banda ancha. La correcta categorización de estos fenómenos es imperativa para el diseño de esquemas de aislamiento coordinados bajo las normativas IEEE Std 1313.2 e IEC 60071-2.

Transitorios de Origen Atmosférico (Rayos)

Las descargas atmosféricas directas o indirectas sobre las líneas aéreas de distribución inducen ondas de sobretensión con frentes de onda extremadamente rápidos (del orden de 0.5 \mu s a 5 \mu s) y colas largas (hasta 50 \mu s). La corriente de rayo se modela típicamente mediante la doble función exponencial de Heidler:

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

Donde $I_p$ es la amplitud de la corriente de cresta, \tau_1 es la constante de tiempo del frente de onda, \tau_2 es la constante de tiempo de la cola de la onda, $n$ es el factor de steepness, y \eta es un factor de corrección de amplitud. Al penetrar estas ondas desde tramos aéreos hacia cables subterráneos, la impedancia característica disminuye drásticamente (típicamente de Z_c \approx 400\,\Omega a Z_c \approx 30\,\Omega), provocando una refracción de tensión que puede duplicar la amplitud de la onda incidente en el punto de transición si la longitud del cable es menor que la longitud crítica del frente de onda.

Maniobras de Apertura y Cierre de Interruptores

La operación de interruptores en vacío (vacuum circuit breakers) o seccionadores bajo carga genera transitorios de alta frecuencia debido al fenómeno de reencendido múltiple del arco eléctrico (multiple restrikes). Cuando el contacto se separa, la corriente cruza por cero de alta frecuencia, y la rigidez dieléctrica del medio de interrupción se restablece a una tasa mayor que la velocidad de ascenso de la tensión de recuperación transitoria (TRV). Esto provoca un corte abrupto de corriente (current chopping), cuya energía inductiva almacenada se transfiere a la capacitancia parásita de la red según la relación energética:

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

Esta sobretensión oscilatoria de alta frecuencia (típicamente entre 1 kHz y 100 kHz) ejerce un esfuerzo dieléctrico severo sobre el aislamiento principal de los devanados de los transformadores.

Resonancia y Ferroresonancia en Sistemas de Distribución

La resonancia armónica y la ferroresonancia representan dos de los modos de fallo más destructivos y complejos en las redes de distribución soterradas y mixtas. Su análisis requiere evaluar la interacción entre los elementos lineales reactivos y las no linealidades ferromagnéticas.

Resonancia Serie y Paralelo

La resonancia paralela (antiresonancia) ocurre cuando la reactancia inductiva equivalente de la red (aportada por transformadores de distribución y reactancias limitadoras) se iguala a la reactancia capacitiva de los cables subterráneos y bancos de condensadores de corrección de factor de potencia a una frecuencia armónica específica $h$ :

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

En este estado, la impedancia vista desde la fuente alcanza un valor máximo teórico limitado únicamente por la resistencia equivalente de la red. Cualquier inyección de corrientes armónicas (provenientes de cargas no lineales como variadores de velocidad y rectificadores) cercana a $h_r$ genera caídas de tensión armónicas masivas, sobrecalentamiento por efecto Joule y deformación severa de la onda de tensión.

Ferroresonancia en Transformadores

La ferroresonancia es un fenómeno no lineal caracterizado por sobretensiones sostenidas de baja frecuencia (fundamental, subarmónica o caótica), distorsión armónica extrema y corrientes elevadas. Se produce típicamente cuando un transformador con núcleo de hierro ferroso insaturable se energiza en vacío a través de una capacitancia en serie (por ejemplo, la capacitancia de un cable subterráneo desconectado en un extremo o fases abiertas por operación monopolar de fusibles).

El comportamiento no lineal del núcleo magnético se modela mediante el flujo enlazado \lambda(i) o la relación constitutiva de inductancia no lineal:

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

La coexistencia de múltiples estados estables (régimen fundamental y regímenes ferroresonantes subarmónicos) se explica mediante la teoría de bifurcaciones en sistemas dinámicos no lineales. La ecuación diferencial que gobierna el circuito LC no lineal excitado por una fuente sinusoidal V_m \sin(\omega t) es:

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)

Donde la inductancia L(i) = \frac{d\lambda}{di} varía drásticamente al entrar la región de saturación magnética, reduciendo el valor de $L$ y desplazando la frecuencia de resonancia del circuito hacia valores que coinciden con la frecuencia de la red o sus subarmónicos.

Análisis Forense de Fallas y Consecuencias Dieléctricas

Los transitorios electromagnéticos de alta frecuencia y las resonancias prolongadas deterioran de manera acumulativa y catastrófica los activos críticos de la red de distribución. A continuación, se detalla el análisis forense de los mecanismos de fallo por componente:

Transformadores de Distribución y Potencia

Cuando un frente de onda rápido (como el generado por rayos o maniobras) penetra en el devanado de alta tensión de un transformador, la distribución inicial de la tensión a lo largo de las espiras no es lineal. Debido a la presencia de la capacitancia serie entre espiras ( $C_s$ ) y la capacitancia a tierra del devanado ( $C_g$ ), el factor de distribución \alpha = \sqrt{C_g / C_s} determina un gradiente de tensión extremadamente elevado en los primeros discos o espiras del devanado.

Esto genera descargas parciales (Partial Discharges - PD) internas en el papel impregnado en aceite, perforación del aislamiento celulósico, y cortocircuitos entre espiras. Térmicamente, las corrientes armónicas inducidas por resonancia paralela provocan pérdidas suplementarias en los devanados por efecto skin (pelicular) y de proximidad, calculadas mediante el factor de pérdidas adicionales $H$ :

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

Cables Subterráneos de Media Tensión

Los cables aislados con Polietileno Reticulado (XLPE) son altamente vulnerables a los transitorios de tensión prolongados y a la ferroresonancia. Las sobretensiones sostenidas aceleran el proceso de treeing eléctrico (arborización eléctrica), un mecanismo de degradación del aislamiento sólido donde microcavidades bajo campos eléctricos intensos (> 10 kV/mm) generan descargas parciales continuas que expanden canales dendríticos hasta la rotura dieléctrica total (ruptura del aislamiento primario).

Dispositivos de Maniobra y Protecciones

Los interruptores en vacío y reconectadores sufren erosión severa de los contactos principales debido a la energía acumulada en los arcos de alta frecuencia durante las operaciones de corte con cargas inductivas o capacitivas desadaptadas. Asimismo, los apartarrayos de óxido de zinc (MOV) pueden experimentar falla térmica por avalancha energética si la densidad de energía absorbida durante un transitorio supera su capacidad de disipación nominal $Wth$ .

Parámetros de Diseño, Límites Normativos y Estrategias de Mitigación

Para garantizar la confiabilidad y la coordinación de aislamiento bajo las directrices internacionales, los ingenieros deben aplicar estrictos criterios de diseño preventivo y mitigar los efectos de las frecuencias naturales.

La siguiente tabla consolida los parámetros eléctricos críticos, los límites normativos bajo estándares IEEE e IEC, las condiciones de fallo crítico y las consecuencias operativas asociadas:

Parámetro Eléctrico / Fenómeno	Estándar de Referencia (IEEE / IEC)	Límite Normativo / Umbral Crítico	Condición de Fallo Crítico	Consecuencias Operativas y Dieléctricas
SobRETENSIONES Temporales (TOV)	IEEE Std 1313.2 / IEC 60071-2	VTOV \le 1.7 \times V_{ph-n} (para t \le 1\, s )	Ferroresonancia sostenida o rechazo de carga en redes débiles	Degradación térmica y envejecimiento prematuro del aislamiento XLPE en cables.
Distorsión Armónica Total de Tensión (THD-V)	IEEE Std 519-2022 / IEC 61000-2-2	THD-V \le 5.0\% (Sistemas de Media Tensión)	Resonancia paralela cercana a los armónicos 5°, 7° u 11°	Sobrecalentamiento de transformadores, disparos en falso de protecciones y fatiga dieléctrica.
Tasa de Subida de Tensión (differentiated steepness $dv/dt$ )	IEC 60034-18-41 / NEMA MG1	dv/dt \le 5\, kV /\mu s (en terminales de motores/cables)	Operación de interruptores en vacío con múltiples reencendidos	Perforación del aislamiento de espiras y fallas inter-espiras en transformadores y máquinas rotativas.
Capacidad de Absorción de Energía en SPD (MOV)	IEEE Std C62.11 / IEC 61643-11	Absorción nominal \ge 2.5\, kJ/kV de tensión asignada	Impactos múltiples de rayos o frentes de onda reflejos amplificados	Falla térmica catastrófica (explosión) del apartarrayos y pérdida de protección.
Factor de Derating por Armónicos en Cables	IEC 60287 / NEC 310.15	Factor de reducción $Kharm$ según espectro de corriente	Corrientes de neutro superiores a la fase por componentes homopolares (3°, 9°)	Fusión del aislamiento por exceso de temperatura en el conductor (> 90°C para XLPE).

Estrategias de Mitigación y Compensación Activa/Pasiva

Para mitigar las frecuencias naturales resonantes y los transitorios severos, se despliegan las siguientes contramedidas de ingeniería:

Instalación de Filtros Armónicos Sintonizados (Filtros Pasivos LC): Se conectan en paralelo con la red para desviar los armónicos críticos (típicamente 5° y 7°) hacia resistencias de descarga, modificando la impedancia equivalente de la red y evitando la resonancia paralela. Su sintonización se realiza ligeramente por debajo de la frecuencia armónica objetivo:
$fsintonizado = \frac{1}{2\pi \sqrt{Lfilter Cfilter}} < farmonico$
Reactancias de Amortiguamiento y Resistencias de Pre-inserción: En interruptores de alta tensión y media tensión crítica, el uso de resistencias de pre-inserción amortigua el frente de onda transitorio, reduciendo el factor de sobretensión a valores $vmax / Vph < 2.0$ .
Suporesores de Sobretensión Basados en Óxido de Zinc (MOV) Ubicados Estratégicamente: Instalados en los puntos de transición entre líneas aéreas y cables subterráneos, así como en los bornes de los transformadores, para recortar los picos de tensión transitoria y limitar la energía incidente.

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

El modelado avanzado de redes de distribución para el estudio de transitorios electromagnéticos y frecuencias naturales requiere plataformas de cálculo computacional de alta precisión como Vexten Suite. A continuación se detalla la metodología algorítmica y los procedimientos de ingeniería implementados en el software para cumplir con los estándares internacionales.

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

El motor de cálculo transitorio de Vexten Suite resuelve el sistema de ecuaciones diferenciales en el dominio del tiempo utilizando el método de integración numérica trapezoidal con paso de tiempo adaptativo (\Delta t < 10\,\mu s). El procedimiento analítico integrado comprende los siguientes pasos:

Adquisición y Modelado Topológico: El sistema importa los parámetros físicos de cables (geometría, capas semiconductoras, aislamiento XLPE), transformadores (impedancias de cortocircuito, corriente de vacío, curva de saturación de Fröhlich-Kennelly) y fuentes de generación distribuida.
Análisis de Frecuencia Natural (掃描 de Impedancia - Frequency Scan): Se inyecta una fuente de corriente armónica de amplitud unitaria y frecuencia variable f \in [10\, Hz , 100\, kHz ] en los nodos críticos de la red. Vexten calcula la impedancia nodal Znn(j\omega) y genera el diagrama Bode de magnitud y fase, identificando automáticamente los picos de resonancia paralela y los valles de resonancia serie.
Simulación de Maniobras y Rayos (EMTP-like Engine): Se configuran eventos de conmutación (apertura de interruptores con corrientes de corte específicas) y ondas de rayo normalizadas (8/20 \mu s y 1.2/50 \mu s). El software evalúa la propagación de ondas viajeras considerando la atenuación modal dependiente de la frecuencia y el efecto pelicular (skin effect) mediante la formulación de Carson modificada o funciones de Bessel para conductores cilíndricos.
Evaluación de Derating Térmico de Cables (IEC 60287): Con base en el espectro de corriente armónica obtenido del análisis de flujo de carga armónico, Vexten calcula automáticamente el factor de reducción de capacidad de conducción de corriente ( $ampacity$ ) del cable subterráneo, asegurando que la temperatura máxima del conductor no supere el límite térmico normativo de 90°C para XLPE bajo régimen distorsionado.

De este modo, Vexten Suite permite al ingeniero especialista realizar un diagnóstico integral, previniendo fallas catastróficas por resonancia y garantizando la coordinación óptima del aislamiento en redes de distribución complejas.