Desplazamiento de Neutro y TOV a Gran Altitud en Minería: Análisis Térmico y Factor K según NFPA 70B

Contextualización Electrofísica de Operaciones Mineras a Gran Altitud

La operación de sistemas eléctricos de potencia en instalaciones mineras situadas a altitudes superiores a los 3000 metros sobre el nivel del mar (msnm) impone restricciones severas a la rigidez dieléctrica del medio aislante y a las capacidades de disipación térmica de los equipos. A altitudes elevadas, la densidad relativa del aire ( $\delta$ ) disminuye drásticamente de acuerdo con la ley barométrica fundamental:

\delta = \frac{P}{P_0} \cdot \frac{T_0}{T} = e^{-\left(\frac{g \cdot M \cdot h}{R \cdot T_0}\right)}

Donde $P$ y $T$ representan la presión y temperatura absoluta a la altitud $h$ (msnm), $P_0$ y $T_0$ son las condiciones estándar a nivel del mar (101.325 kPa y 288.15 K), $g$ es la aceleración de la gravedad, $M$ es la masa molar del aire seco y $R$ es la constante universal de los gases. Esta reducción de densidad altera la trayectoria media libre de los electrones libres, incrementando la probabilidad de ionización por colisión bajo la aplicación de un campo eléctrico, según la Ley de Paschen:

V_b = \frac{B \cdot (p \cdot d)}{\ln(A \cdot p \cdot d) - \ln\left[\ln\left(1 + \frac{1}{\gamma_{se}}\right)\right]}

Donde $V_b$ es la tensión de ruptura, $p$ la presión atmosférica, $d$ la distancia entre electrodos, $A$ y $B$ son constantes compuestas del gas y $\gamma_{se}$ es el segundo coeficiente de ionización de Townsend. Como resultado directo, las distancias de aislamiento en el aire (fases-tierra y fase-fase) así como las distancias de fuga (creepage distance) experimentan una degradación no lineal.

De acuerdo con la norma IEC 60071-2 e IEEE C37.100.1, el factor de corrección por altitud para el aislamiento externo ( $K_a$ ) se define según la formulación matemática:

K_a = e^{m \cdot \left(\frac{h - 1000}{8150}\right)}

Donde $m$ es un parámetro empírico que depende del tipo de tensión aplicada (frecuencia industrial, impulso de maniobra o impulso de rayo) y del valor de la tensión de soporte. Para instalaciones mineras situadas a 4500 msnm, el factor $K_a$ típicamente oscila entre 1.45 y 1.60, lo que exige sobredimensionar la tensión asignada soportable a la frecuencia industrial ( $V_w$ ) de interruptores, aisladores y celdas de media tensión según:

V_w \ge V_{req, sea\_level} \times K_a

Simultáneamente, la disipación térmica por convección natural se ve drásticamente atenuada debido a la masa volumétrica reducida del fluido refrigerante (aire). El factor de derating térmico ( $kalt$ ) aplicado a la corriente nominal de equipos como transformadores, barras conductoras e interruptores se aproxima mediante la norma IEEE C57.96:

kalt = 1 - \alpha \cdot \left( \frac{h - 1000}{100} \right)

Donde $\alpha \approx 0.005$ a $0.008$ según la geometría del elemento radiante y la relación entre pérdidas por radiación y convección. La interacción simultánea de la degradación dieléctrica y el derating térmico establece un escenario crítico para la aparición de sobretensiones temporales (TOV) y fenómenos de desplazamiento de neutro en redes con neutro aislado o puesto a tierra mediante alta impedancia (HRG).

Física del Desplazamiento de Neutro y Sobretensiones Temporales (TOV)

El desplazamiento de neutro es un fenómeno electrodinámico que ocurre en sistemas trifásicos no rígidamente aterrados cuando se rompe la simetría de las admitancias de fase a tierra o se presenta una falla monofásica a tierra (SLG - Single Line-to-Ground). Considérese un sistema de distribución minero en media tensión alimentando cargas no lineales donde el neutro del sistema $N$ está conectado a tierra a través de una admitancia $Y_n = 1/RNGR + j\omega C_n$ .

Aplicando la ley de corrientes de Kirchhoff en el nodo de neutro $N$ respecto al potencial de tierra de referencia $G$ , la tensión de desplazamiento de neutro ( $VNG$ ) ante un desequilibrio de admitancias fasoriales por fase ( $Y_A, Y_B, Y_C$ ) y tensiones generadas ( $E_A, E_B, E_C$ ) se expresa rigurosamente como:

VNG = \frac{E_A Y_A + E_B Y_B + E_C Y_C}{Y_A + Y_B + Y_C + Y_n}

Donde la admitancia de cada fase $i \in \{A, B, C\}$ incluye la capacitancia distribuida del sistema de cables subterráneos/apantallados ( $C_{0,i}$ ), las admitancias de aislamiento y los componentes de filtro pasivo conectadas a esa fase:

Y_i = G_i + j\omega C_{0,i} + \sumh \frac{1}{R_{h,i} + j\left(h\omega L_{h,i} - \frac{1}{h\omega C_{h,i}}\right)}

Bajo una condición de falla monofásica a tierra de alta impedancia o franca en la fase $A$ ( $R_f \to 0$ ), la admitancia $Y_A \to \infty$ . El límite de la ecuación de desplazamiento de neutro conduce a $VNG \to -E_A$ . Por lo tanto, el vector de tensión neutro-tierra se iguala a la magnitud de la tensión fase-neutro prefalla pero con signo opuesto. Las tensiones fasoriales de las fases no falladas ( $VBG$ y $VCG$ ) experimentan una elevación súbita alcanzando la tensión fase-fase prefalla ( $VLL$ ):

VBG = E_B - VNG = E_B + E_A = \sqrt{3} E_B e^{-j 30^\circ}

VCG = E_C - VNG = E_C + E_A = \sqrt{3} E_C e^{+j 30^\circ}

El factor de sobretensión temporal ( $kTOV$ ) para las fases no falladas en un sistema HRG se define mediante la razón entre la tensión pico fase-tierra durante la falla y la tensión pico nominal fase-neutro:

kTOV = \frac{|VBG|}{|E_B|} = \sqrt{3} \approx 1.732 p.u.

En el dominio de las componentes simétricas, la sobretensión temporal por desplazamiento de neutro se analiza evaluando las redes de secuencia positiva ( $Z_1$ ), secuencia negativa ( $Z_2$ ) y secuencia cero ( $Z_0$ ). La relación de tensiones durante una falla monofásica a tierra viene dada por:

VBG = a^2 E_A - \left( \frac{Z_2 + Z_0 + 3R_f}{Z_1 + Z_2 + Z_0 + 3R_f} a^2 + \frac{Z_2 - a Z_0}{Z_1 + Z_2 + Z_0 + 3R_f} a + \frac{Z_0 (1-a^2)}{Z_1 + Z_2 + Z_0 + 3R_f} \right) E_A

Cuando la relación $X_0 / X_1 \gg 3$ y $R_0 / X_1 \gg 1$ (característica de sistemas no aterradas o aterrados por alta resistencia), el sistema se clasifica como no efectivamente aterrado según el estándar IEEE 142. Bajo estas condiciones, el coeficiente $kTOV$ no se limita a $1.38$ p.u. (límite de redes efectivamente aterradas), sino que alcanza niveles superiores a $1.732$ p.u. estacionariamente, con transitorios de oscilación LC de alta frecuencia que pueden alcanzar picos de:

V_{peak, transitorio} = Vprefalla \cdot \left( 1 + \eta \cdot \sin(\omega_r t) e^{-\frac{t}{\tau}} \right) \approx 2.5 a 3.0 p.u.

La presencia de capacitancias distribuidas elevadas en alimentadores de minería subterránea o tajos abiertos (cables de arrastre SHD-GC de alta capacitancia lineal $C_0 \approx 0.3 - 0.6 \, \mu F/km$ ) combinada con la inductancia no lineal de transformadores de potencial o transformadores de potencia no cargados produce un acoplamiento ferroresonante. La ferroresonancia paralelo se desencadena cuando la reactancia capacitiva de secuencia cero ( $XC0 = 1 / (\omega C_0)$ ) iguala a la reactancia magnetizante no lineal del transformador saturado ( $X_m(i)$ ):

f_r = \frac{1}{2\pi \sqrt{L_m(i) \cdot C_0}}

Bajo ferroresonancia, el neutro oscila caóticamente, generando sobretensiones extremas sostenidas ( $TOV > 3.0 p.u.$ ) con frecuencias subarmónicas ( $1/3 f_1$ , $1/2 f_1$ ) o armónicas que superan con frecuencia la tensión de soporte dieléctrico del aire corregida por altitud ( $V_b / K_a$ ).

Armónicos de Carga Minera y Factor K según NFPA 70B e IEEE C57.110

Las plantas de procesamiento minero operan con cargas no lineales de alta potencia, tales como variadores de frecuencia de media tensión (VFDs tipo VSI con topología neutral point clamped NPC de 3 niveles, 6/12/18 pulsos o AFE), accionamientos de molinos SAG y de bolas accionados por cicloconvertidores o motores síncronos alimentados por inversores LCI. Estas cargas inyectan corrientes armónicas de orden espectral:

h = p \cdot k \pm 1

Donde $p$ es el número de pulsos del convertidor y $k \in \{1, 2, 3, \dots\}$ . Para rectificadores de 6 pulsos, $h \in \{5, 7, 11, 13, 17, 19, \dots\}$ . En sistemas con desequilibrios o trifásicos de 4 hilos, emergen triplens armónicos ( $h = 3, 9, 15, 21, \dots$ ), los cuales poseen naturaleza de secuencia cero ( $Z_0$ ).

Las corrientes triplen armónicas no se cancelan en el nodo neutro de los transformadores en conexión estrella con neutro accesible; por el contrario, se suman aritméticamente en el conductor de neutro:

I_{N, rms} = \sqrt{3 \sum_{k=1}^{\infty} \left( I_{3(2k-1)} \right)^2 + I_{N, fundamental}^2}

La circulación de corrientes armónicas por los enrollamientos de los transformadores causa un calentamiento severo debido al incremento no lineal de las pérdidas por corrientes de Foucault en el cobre (eddy current losses, $PEC$ ) y pérdidas suplementarias en partes estructurales (stray losses, $POS$ ). El Factor K ( $K -Factor$ ), definido en IEEE C57.110 y adoptado en la norma NFPA 70B para la evaluación de integridad de aislamiento en mantenimiento preventivo/predictivo, cuantifica el efecto de calentamiento de corrientes armónicas en relación con la corriente fundamental:

K = \sum_{h=1}^{hmax} h^2 \left( \frac{I_h}{I_1} \right)^2 = \frac{\sum_{h=1}^{hmax} (h \cdot I_h)^2}{\sum_{h=1}^{hmax} I_h^2}

Donde $I_h$ es la corriente RMS de la armónica de orden $h$ , e $I_1$ es la corriente RMS fundamental. El factor de pérdidas por corrientes de Foucault ( $FHL$ ) utilizado para deratear la capacidad nominal del transformador se expresa como:

FHL = \frac{\sum_{h=1}^{hmax} h^2 \left( \frac{I_h}{I} \right)^2}{\sum_{h=1}^{hmax} \left( \frac{I_h}{I} \right)^2} = \frac{K}{1 + THD_I^2}

Donde $THD_I = \frac{\sqrt{\sum_{h=2}^{\infty} I_h^2}}{I_1}$ es la Distorsión Total de Armónicos de corriente. La capacidad de carga reducida del transformador en per-unit ( $Pmax(pu)$ ) bajo entorno armónico y de gran altitud se recalcula mediante la combinación de IEEE C57.110 y el factor de altitud:

Pmax(pu) = kalt \cdot \sqrt{ \frac{PLL(pu)}{1 + FHL \cdot P_{EC-R}(pu)} }

Donde $PLL(pu)$ son las pérdidas totales a plena carga calculadas en pu, e $P_{EC-R}(pu)$ son las pérdidas por corrientes de Foucault en el devanado calculadas a frecuencia fundamental en condiciones nominales.

Según el estándar NFPA 70B (Recommended Practice for Electrical Equipment Maintenance), la acumulación de distorsión armónica elevada combinada con elevaciones de temperatura aceleran el proceso de degradación del aislamiento térmico mediante la ecuación de Arrhenius para la vida útil esperada del aislamiento ( $\mathcal{L}$ ):

\mathcal{L} = \mathcal{L}_0 \cdot \exp\left( \frac{E_a}{R \cdot Thotspot} \right)

Donde $E_a$ es la energía de activación del dieléctrico, $R$ la constante de los gases y $Thotspot$ la temperatura del punto más caliente del devanado, la cual se eleva dramáticamente según:

Thotspot = T_{amb, alt} + \Delta T_{A-R} \cdot \left[ \frac{PDC + FHL \cdot P_{EC-R}}{PDC + P_{EC-R}} \right]^y + \Delta THVR

NFPA 70B exige para instalaciones industriales a gran altitud la medición obligatoria mediante termografía infrarroja ajustada por transmitancia atmosférica corregida por altitud, análisis de gases disueltos en aceite (DGA) con umbrales reducidos de formación de etileno y acetileno debido a descargas parciales inducidas por TOV, y la verificación periódica del Factor K del sistema de distribución.

Análisis Forense de Fallas Electromecánicas Integradas

La combinación del desplazamiento de neutro, las sobretensiones temporales (TOV), la atenuación dieléctrica por baja presión atmosférica y el esfuerzo térmico armónico desencadena mecanismos de falla multifásicos en el equipamiento de la mina. A continuación se realiza un análisis forense detallado por tipo de activo crítico:

Aislamiento en Cables de Potencia (XLPE / EPR)

En cables de alta y media tensión (p. ej. 13.8 kV o 34.5 kV aislados en XLPE), un desplazamiento de neutro sostiene una tensión fase-tierra de $\sqrt{3} VLN$ en los conductores de las fases sanas durante horas si el sistema HRG no elimina la falla rápidamente. La tensión de pico aplicada al dieléctrico interno es:

V_{peak, actual} = \sqrt{2} \cdot \left( \sqrt{3} VLN \right) \cdot \left( 1 + THD_V \right)

Esta sobretensión reduce drásticamente la Tensión de Incepción de Descargas Parciales (CIV - Partial Discharge Inception Voltage). La presencia de microcavidades (voids) dentro del aislante XLPE genera un gradiente de campo eléctrico local $Evoid = \epsilon_{r} \cdot E_0$ . Cuando $Evoid$ supera la tensión de ruptura del gas contenido dentro de la cavidad (cuya rigidez cae según la ley de Paschen si existe migración de presión o gradiente de altitud en terminaciones), se inician descargas parciales sostenidas. Estas descargas erosionan la matriz polimérica formando arborescencias eléctricas (electrical treeing), acelerando la ruptura dieléctrica catastrófica del cable.

Celdas de Media Tensión (Switchgear) y Aislamiento en Aire

El aire dentro del armario del switchgear actúa como medio aislante primario entre las barras fásicas y el bastidor aterrado. A 4200 msnm, la rigidez dieléctrica del aire se reduce en un factor aproximado de $\delta \approx 0.62$ . Durante un evento de TOV asociado a una falla a tierra en un sistema HRG con alta presencia de 3.º y 9.º armónicos:

La distancia de aislamiento (clearance) requerida para soportar impulsos de maniobra o transitorios de reencendido de arco debe incrementarse en un $45\% - 60\%$ según IEC 60664-1 Tabla A.2.
Si la celda fue diseñada para operación a nivel del mar ( $V_{\text{disrupción}} = 95 kV BIL$ ), a 4200 msnm la tensión de soporte BIL real desciende a $95 \times 0.62 = 58.9 kV BIL$ .
Un transitorio TOV con pico de 48 kV fase-tierra superpuesto con armónicos de alta frecuencia sobrepasa el límite dieléctrico del aire degradado, provocando un cebado de arco fase-tierra que evoluciona rápidamente a un cortocircuito trifásico destructivo por ionización masiva del volumen del armario (Arc Flash Event).

Transformadores de Potencia y Factor K

En transformadores, el aumento de la tensión de neutro eleva el punto de operación de la curva de magnetización del núcleo ( $B-H$ ). Si la TOV eleva la densidad de flujo $B$ por encima del punto de saturación ( $Bsat \approx 1.7 - 1.8 Tesla$ ):

B(t) = \frac{1}{N \cdot A} \int0^{t} VTOV(\tau) d\tau

La corriente de magnetización se convierte en impulsos picudos ricos en armónicos impares y triplens ( $Imag \propto B^{n}, n > 9$ ). Esto incrementa exponencialmente las pérdidas en el hierro ( $P_Fe \propto f \cdot B^2 + f^2 \cdot B^2$ ) e induce mayores pérdidas por corrientes de Foucault ( $PEC$ ). El calor no disipado debido a $kalt$ degrada térmicamente el papel aislante celulósico, liberando monóxido y dióxido de carbono ( $CO, CO₂$ ), y generando puntos calientes que inducen la formación de gas de furanos, reduciendo el Grado de Polimerización (DP) del papel de 1000 a valores críticos (< 200), lo que señala el fin de vida mecánico e inminente falla entre espiras del devanado.

Parámetro Electromecánico / Aislamiento	Condición Estándar (0 msnm, $THD_I < 5\%$ )	Condición Severa Minera (4500 msnm, $K \ge 13, TOV = 1.73 \, pu$ )	Consecuencia Térmica y Dieléctrica Impuesta
Rigidez Dieléctrica del Aire ( $E_b$ )	$3.0 kV/mm$	$\sim 1.85 kV/mm$ ( $\delta \approx 0.61$ )	Ruptura dieléctrica fase-tierra en celdas de media tensión y aisladores.
Tensión Fase-Tierra en Fases Sanas ( $VLG$ )	$1.0 p.u.$ ( $VLL / \sqrt{3}$ )	$1.732 p.u.$ sostenida (o $> 2.5 p.u.$ transitoria)	Envejecimiento acelerado del XLPE por descargas parciales (treeing).
Factor K de Pérdidas Armónicas	$K = 1.0$	$K = 13.0 - 20.0$	Pérdidas $PEC$ incrementadas en un $1300\% - 2000\%$ en devanados.
Temperatura de Punto Caliente ( $Thotspot$ )	$98^\circ C - 110^\circ C$	$> 145^\circ C$ (sin derating específico)	Destrucción de cadena polimérica del papel Kraft (DP < 200).
Tensión Incepción Descargas Parciales (CIV)	$100\%$ valor nominal	$\sim 60\% - 65\%$ del valor nominal	Generación sostenida de corona y degradación de aislamiento externo.
Capacidad de Disipación Radiadores ( $h_c$ )	$100\%$ eficiencia convectiva	$\sim 70\%$ eficiencia convectiva	Sobrecalentamiento generalizado del fluido aislante/refrigerante.

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

Para garantizar la continuidad operativa y la seguridad integral de la red eléctrica minera bajo condiciones combinadas de gran altitud, distorsión armónica y riesgo de TOV por desplazamiento de neutro, se deben implementar las siguientes metodologías de diseño e ingeniería de detalle:

Dimensionamiento de la Resistencia de Puesta a Tierra de Neutro (NGR) para Control de TOV y Altitud

En sistemas HRG a gran altitud, la corriente elegida para la resistencia de puesta a tierra ( $INGR$ ) debe ser estrictamente mayor que la corriente total de carga capacitiva de la red en secuencia cero ( $IC0$ ):

INGR \ge IC0 = 3 \cdot \omega \cdot C_{0, total} \cdot VLN

Donde $C_{0, total}$ es la capacitancia acumulada de todos los cables, bancos de condensadores y capacitancias parásitas de la red. Al garantizar $RNGR \le \frac{1}{3 \omega C_{0, total}}$ , se disipa la energía capacitiva almacenada durante la falla monofásica, amortiguando completamente los picos de sobretensión por reencendido de arco (arcing grounds) y limitando $kTOV \le 1.732 p.u.$ .

Adicionalmente, el elemento físico R (banco de resistencias de aleación de acero inoxidable o cromo-níquel) debe dimensionarse considerando el derating térmico por altitud para la disipación del calor generado durante una falla sostenida ( $tfalla = 10 s$ o continuo según NFPA 70B):

R_{NGR, derated} = R_{NGR, 20^\circ C } \left[ 1 + \alpha_{\text{T}} \cdot (\Delta Talt) \right]

P_{NGR, nominal} \ge \frac{VLN^2}{RNGR} \times \frac{1}{k_{alt, NGR}}

Selección y Coordinación de Aislamiento de Pararrayos (Surge Arresters)

La selección de pararrayos de óxido de zinc (ZnO) en instalaciones mineras a alta altitud sujetas a TOV requiere una evaluación rigurosa del margen entre la Tensión de Operación Continua Máxima (MCOV o $V_c$ ) y la curva de capacidad TOV del pararrayos ( $VTOV(t)$ ). La MCOV debe ser superior a la tensión máxima fase-tierra en presencia de desplazamiento de neutro:

VMCOV \ge V_{LL, max} \quad (Para redes HRG con eliminación de falla retardada)

El pararrayos debe ser capaz de absorber la energía de sobretensión sin sufrir embalamiento térmico (thermal runaway). La capacidad de absorción de energía específica ( $E_K$ en kJ/kV de MCOV) se descalifica por altitud debido a la reducida refrigeración del envolvente polimérico:

E_{K, min} = \frac{Etransitorio + ETOV}{k_{alt, arrester}}

Asimismo, la distancia de fuga externa del pararrayos se dimensiona aplicando un factor mínimo de $31 mm/kV$ de tensión máxima del sistema, corregido por el factor $K_a$ para ambientes con contaminación minera severa (polvo conductor metálico o salino).

Filtrado Armónico y Desintonización para Evitar Resonancia Paralelo

Para mitigar la inyección del espectro armónico y prevenir el sobrecalentamiento de transformadores por Factor K, se diseñan filtros pasivos desintonizados (detuned filters) o filtros tipo C (C-type damped filters). La frecuencia de sintonía ( $fst$ ) debe ubicarse por debajo de la quinta armónica ( $h=5$ ), típicamente en $h_r = 4.2 - 4.7$ , para evitar resonancia paralelo entre los bancos de condensadores y la inductancia de la red ( $Lsys + Ltrafo$ ):

f_r = f_1 \cdot \sqrt{ \frac{X_C}{X_L + Xsys} }

En presencia de variadores de velocidad de alta potencia, la implementación de Filtros Activos de Potencia (APF) de topología Multinivel conectados en paralelo inyecta en oposición de fase la corriente armónica exacta requerida por la carga ( $IAPF(t) = -Ih(t)$ ), garantizando que $THD_I < 5\%$ en el Punto de Acoplamiento Común (PCC) según IEEE 519, reduciendo efectivamente el Factor K visto por el transformador de alimentación a $K \approx 1.0$ .

Aplicación y Metodología de Cálculo con Vexten Suite

La plataforma avanzada de ingeniería eléctrica Vexten Suite integra módulos de simulación numérica y analítica diseñados para resolver simultáneamente las ecuaciones complejas de cortocircuito, derating de aislantes por altitud, flujo armónico y capacidad de conducción de cables en entornos críticos como la minería de gran altitud. A continuación se detalla la metodología de cálculo y el flujo de diseño aplicado mediante los módulos de Vexten Suite.

Módulo de Cortocircuito e Impedancias Vexten (IEC 60909 / IEEE 141)

El módulo Vexten Short-Circuit & Grounding Analyzer calcula la matriz de impedancias de secuencia ( $Z_1, Z_2, Z_0$ ) ajustada por la temperatura operativa de los conductores y la resistividad del terreno a alta altitud. Durante una falla monofásica a tierra en una red HRG minera, el software resuelve la corriente inicial de cortocircuito simétrica ( $Ik1''$ ) e ignora los efectos de la impedancia de carga convencional según IEC 60909-0:

Ik1'' = \frac{\sqrt{3} \cdot c \cdot V_n}{Z_1 + Z_2 + Z_0 + 3Z_f}

Dado que en una red HRG el término de impedancia de secuencia cero está dominado por la resistencia del neutro ( $Z_0 \approx 3 RNGR$ ), la corriente de falla a tierra se aproxima en Vexten por:

Ik1'' \approx \frac{c \cdot V_n}{\sqrt{3} RNGR}

Simultáneamente, el motor de Vexten calcula el perfil fasorial de tensiones en las barras sanas, computando directamente el factor $k_{TOV, bar}$ y evaluando la desviación del neutro $VNG$ para cada escenario de operabilidad de la mina.

Módulo de Sizing de Cables e Incremento por Armónicos Vexten (IEC 60287 / NEC 310)

El módulo Vexten Cable Sizing & Harmonic Derating Engine calcula la ampacidad efectiva ( $Iamp$ ) de cables unipolares o tripolares de media y baja tensión (p. ej. tipo SHD-GC o XLPE) bajo la acción simultánea de la altitud y la corriente de neutro armónica mediante la ecuación multidimensional de derating:

I_{amp, corregida} = Itabulada \cdot ktemp \cdot kagrup \cdot kalt \cdot kharm

Donde el factor de degradación por corrientes armónicas en el neutro ( $kharm$ ) se determina internamente en Vexten evaluando el calentamiento del conductor neutro y de los conductores de fase de acuerdo con la norma NEC 310.15(E) e IEC 60287-1-1. Si el porcentaje de armónicos de tercer orden ( $h=3$ ) supera el $33\%$ , la corriente en el neutro supera la corriente de fase, y Vexten ajusta automáticamente el tamaño del conductor de neutro a $150\% - 200\%$ de la sección de la fase, recalculando el gradiente térmico capa por capa:

\Delta \theta = (Pcond + Pdiel) \cdot T_1 + \left[ Pcond(1 + y_s + y_p) + Pdiel \right] \cdot n \cdot (T_2 + T_3 + T_4)

Donde $T_1, T_2, T_3, T_4$ representan las resistencias térmicas de las capas del cable (aislante, cubierta, lecho y terreno/aire) modificadas en el código de Vexten según la densidad volumétrica del aire a la altitud ingresada por el usuario.

Módulo de Factor K y Análisis de Resonancia Vexten (IEEE C57.110 / IEEE 519)

El módulo Vexten Harmonic & Resonance Analyzer realiza el escaneo de impedancia espectral del sistema ( $Z(f)$ ) desde $50 Hz$ hasta $2500 Hz$ en pasos de $1 Hz$ . El algoritmo detecta los puntos de inflexión donde la derivada de la impedancia cambia de signo ( $\frac{d|Z|}{df} = 0$ ), identificando frecuencias de resonancia paralelo ( $f_p$ ) y serie ( $f_s$ ).

A partir del espectro armónico de corriente introducido o exportado desde medidores de calidad de energía (archivos PQDIF/COMTRADE), Vexten calcula automáticamente:

El Factor K del sistema según IEEE C57.110.
El factor de pérdidas por corrientes de Foucault $FHL$ .
La capacidad de carga derateada del transformador en kVA ( $kVA_{derated} = kVA_{nominal} \cdot Pmax(pu)$ ).
El perfil de elevación de temperatura del punto más caliente ( $Thotspot$ ) e índice de pérdida de vida útil según el modelo Arrhenius de la NFPA 70B.

Algoritmo de Cálculo Integrado en Vexten Suite

El siguiente procedimiento describe el flujo de cálculo computacional ejecutado por Vexten Suite para resolver la interacción entre el aislamiento por altitud, el desplazamiento de neutro y la distorsión por Factor K:

\begin{matrix} \text{\textbf{Paso 1: Lectura de Datos}} & \longrightarrow & h (msnm) , Tamb, Espectro Armónico \{h, I_h\}, Topología de Neutro (HRG/Solido) \\ \Downarrow & & \\ \text{\textbf{Paso 2: Corrección Dieléctrica}} & \longrightarrow & \delta = e^{-\frac{g M h}{R T_0}} \quad \implies \quad K_a = e^{m \frac{h-1000}{8150}} \quad \implies \quad V_{BIL, req} = V_{BIL, std} \cdot K_a \\ \Downarrow & & \\ \text{\textbf{Paso 3: Análisis Desplazamiento}} & \longrightarrow & VNG = \frac{\sum E_i Y_i}{\sum Y_i + Y_n} \quad \implies \quad kTOV = \frac{|VBG|}{|E_B|} \quad \implies \quad Check MCOV Arresters \\ \Downarrow & & \\ \text{\textbf{Paso 4: Evaluación Factor K}} & \longrightarrow & K = \frac{\sum (h I_h)^2}{\sum I_h^2} \quad \implies \quad FHL = \frac{K}{1+THD_I^2} \quad \implies \quad Pmax(pu) = kalt \sqrt{\frac{PLL}{1+FHL P_{EC-R}}} \\ \Downarrow & & \\ \text{\textbf{Paso 5: Optimización de Diseño}} & \longrightarrow & Sizing NGR (INGR \ge 3 IC0) \quad \& \quad Diseño Filtro Desintonizado (fst < 5th) \end{matrix}

Mediante esta arquitectura computacional, Vexten Suite provee una herramienta integral de alta precisión que elimina los riesgos de fallas por sobretensiones temporales, arcos eléctricos dieléctricos y colapsos térmicos de aislamiento en sistemas de potencia mineros de alta complejidad a gran altitud.