Variación Estacional de Puesta a Tierra: Sequía, Congelamiento y Tratamiento con Geles Conductivos (IEEE 80 / IEEE 142)

Termodinámica y Transporte Iónico del Suelo: Mecanismos Físicos de la Variación Estacional

La resistividad eléctrica del suelo ($\rho$) no es un parámetro estático de diseño; es una variable termodinámica y electroquímica altamente dependiente del estado higrométrico, la concentración de sales disueltas, la porosidad de la matriz geológica y la temperatura absoluta del medio. En sistemas de potencia, subestaciones eléctricas y plantas de generación, ignorar las oscilaciones estacionales introduce discrepancias críticas entre los cálculos nominales basados en normas como IEEE Std 80 o IEC 60364-5-54 y el comportamiento real del sistema de puesta a tierra (SPT) ante eventos de cortocircuito o descargas atmosféricas.

Conducción electrolítica, porosidad y saturación hídrica (Ecuación de Archie generalizada)

La conducción de corriente eléctrica en medios porosos no saturados ocurre predominantemente a través de la fase líquida contenida en el espacio poroso intersticial mediante transporte iónico migratorio. La matriz mineral de cuarzo, feldespatos o silicatos suele actuar como un aislante dieléctrico con resistividades intrínsecas superiores a $10^6\ \Omega\cdot m$. La conductividad global del terreno ($\sigma = 1/\rho$) se modela mediante una extensión de la Ley de Archie para medios no saturados y multicomponentes:

Donde:

• $\rho_w$: Resistividad del electrolito poral ($\Omega\cdot m$), gobernada por la actividad y movilidad de los iones ($Na ^+$, $Ca ^{2+}$, $Mg ^{2+}$, $Cl ^-$, $SO _4^{2-}$, $HCO _3^-$).
• $\phi$: Porosidad total del suelo (fracción volumétrica de vacíos, $0 < \phi < 1$).
• $S_w$: Grado de saturación de agua en los poros ($Vagua / V_{\text{vacíos}}$, $0 \le S_w \le 1$).
• $a$: Factor de tortuosidad lítica (típicamente entre 0.5 y 1.5).
• $m$: Factor de cementación de la matriz porosa ($1.3 \le m \le 2.5$).
• $n$: Exponente de saturación hídrica (generalmente $n \approx 2.0$).

Bajo condiciones de sequía severa, la succión mátrica ($\psi_m$) asciende exponencialmente según la curva de retención de agua en el suelo (SWRC, por sus siglas en inglés), expulsando el agua capilar libre. Conforme $S_w \to 0$, la continuidad percolativa de las trayectorias de conducción se rompe abruptamente (umbral de percolación), forzando a que la resistividad aparente del terreno ($\rho_{app}$) se incremente en hasta 2 a 3 órdenes de magnitud con respecto a su estado de saturación de diseño.

Termocriogenia del terreno: Transición de fase del agua y efecto del permafrost estacional

El descenso térmico por debajo del punto de congelación ($T < 0\ ^\circ C$) desencadena una transición de fase termodinámica de primer orden en el agua libre de los macroporos. El hielo posee una estructura cristalina hexagonal cerrada donde la movilidad protónica mediante el mecanismo de Grotthuss es extremadamente baja, exhibiendo una resistividad intrínseca en el orden de $10^5$ a $10^7\ \Omega\cdot m$.

Donde $\theta_u(T)$ representa el contenido volumétrico de agua líquida no congelada remanente en las microporosidades a la temperatura $T$, $\theta_0$ es el contenido de agua total previo al congelamiento, $T_f$ es la temperatura de congelación deprimida por la salinidad poral, y $\beta, \gamma$ son coeficientes empíricos texturales. La presencia de capas de hielo superficiales no solo anula el transporte por difusión iónica sino que induce el fenómeno de crio-succión, atrayendo humedad de horizontes inferiores hacia el frente de congelación activo, desecando y rigidizando dieléctricamente las capas subyacentes inmediatas al electrodo.

Dependencia térmica de la movilidad iónica y Ley de Arrhenius

Por encima del punto de congelación ($T > 0\ ^\circ C$), la resistividad del electrolito está controlada por la viscosidad dinámica ($\eta$) del solvente acuoso y la difusividad de los iones según la relación de Nernst-Einstein. La variación térmica de la resistividad puede expresarse formalmente mediante el modelo semi-linealizado o el formalismo de Arrhenius:

Donde $\alpha_T$ es el coeficiente de temperatura de la resistividad del suelo (típicamente $\alpha_T \approx 0.02$ a $0.03\ K ^{-1}$ para electrólitos diluidos), $E_a$ es la energía de activación para el transporte iónico ($\approx 0.15 - 0.22\ eV$), y $k_B$ es la constante de Boltzmann. Una reducción de temperatura de $25\ ^\circ C$ a $1\ ^\circ C$ incrementa la resistividad del suelo en más del 70% exclusivamente por reducción de la agitación térmica y de la cinética de deriva electroforética.

Modelado Matemático de la Resistividad Equivalente y Parámetros de Red de Tierra

Las oscilaciones climatológicas modifican el suelo en profundidad conformando estratigrafías dinámicas. El modelado no puede considerarse homogéneo ($\rho = constante$), sino que requiere formularse como un problema de estratificación horizontal bicapa o multicapa dependiente del tiempo.

Estratificación estacional del suelo: Perfil bicapa dinámico variable en el tiempo

Durante la temporada de estiaje o congelamiento, la capa superficial de profundidad $h_1$ experimenta una alteración extrema en su resistividad ($\rho_1(t)$), mientras que los estratos profundos ($\rho_2$) permanecen relativamente protegidos de las influencias meteorológicas debido a la inercia térmica e hidrológica de la corteza. El coeficiente de reflexión electromagnético interfacial $K$ se define analíticamente:

El potencial electrostático $V(r)$ en la superficie de un suelo bicapa debido a una inyección puntual de corriente $I$ se obtiene resolviendo la ecuación de Laplace ($\nabla^2 V = 0$) en coordenadas cilíndricas mediante la integral de Hankel:

Cuando ocurre un congelamiento severo o sequía superficial extrema, $\rho_1 \gg \rho_2$, implicando que $K \to -1$. Por el contrario, en un lecho rocoso subyacente con superficie humedecida, $K \to +1$. En ambos casos límite, las series infinitas de coeficientes de reflexión condicionan la distribución de gradientes de potencial en la superficie y la dispersión tridimensional de corriente hacia tierra profunda.

Fórmulas analíticas de resistencia de puesta a tierra bajo condiciones límite (Schwarz, Sverak e IEEE 80)

Para una red de malla combinada con conductores horizontales y electrodos verticales embebidos en un suelo no homogéneo, el método analítico riguroso de C. W. Schwarz calcula la resistencia total de puesta a tierra ($R_g$) considerando el acoplamiento mutuo entre la retícula horizontal ($R_1$) y las picas verticales ($R_2$):

Donde los componentes se definen a partir de la longitud total del conductor de malla ($L_c$), el área ocupada ($A$), la profundidad de enterramiento ($h$), el radio del conductor ($r$), el número de picas verticales ($n_R$), y la longitud unitaria de pica ($L_r$):

Donde $k_1, k_2$ son coeficientes geométricos de forma de la malla, $b$ es el radio de la pica vertical y $\rho_{app}$ es la resistividad aparente ponderada que debe sintetizar la degradación de la capa superior. Si la capa superior de espesor $h_1 \ge h$ se congela o deseca, la contribución de la componente $R_1$ queda virtualmente aislada ($R_1 \to \infty$), forzando a que la totalidad de la corriente $I_G$ deba disiparse a través del subsistema de electrodos verticales $R_2$, sobredimensionando la densidad de corriente superficial $J_s$ y disparando la resistencia equivalente global.

Cálculo de potenciales de paso, contacto y GPR crítico durante eventos climáticos adversos

De acuerdo con la norma IEEE Std 80, la tensión de paso tolerable ($Estep$) y la tensión de contacto tolerable ($Etouch$) para un ser humano de 50 kg o 70 kg están gobernadas por el factor de reducción superficial $C_s$, el cual depende de la resistividad de la capa superficial del terreno ($\rho_s$) y del espesor de la grava/acondicionamiento ($h_s$):

Durante épocas de sequía, la deshidratación de la subcapa inmediata ($\rho \gg \rho_s$) degrada el valor del potencial de malla ($E_m$) y de paso ($E_s$). El incremento estacional de la resistencia global de la red $R_g(t)$ eleva el Potencial de Elevación de Tierra (GPR, Ground Potential Rise):

Donde $I_G$ es la corriente simétrica de falla a tierra de diseño, $D_f$ es el factor de decremento por asimetría de la componente de corriente directa (DC subtransitoria), y $S_f$ es el factor de división de corriente que drena por los cables de guarda. Un incremento estacional de $R_g(t)$ de un factor de $\times 4$ a $\times 10$ causa que el GPR exceda ampliamente las tensiones de soporte dieléctrico de los aislamientos secundarios, cables de control y cubiertas de telecomunicaciones acopladas al SPT.

Análisis Forense de Fallas Eléctricas Inducidas por Degradación Estacional del Terreno

El fallo de un sistema de puesta a tierra debido a variaciones climáticas no se manifiesta de forma lineal o predecible mediante mediciones estáticas tradicionales en primavera u otoño; su manifestación es destructiva y cataclísmica durante fallas asimétricas monofásicas ($1\Phi-G$) o descargas directas tipo rayo coincidentes con el pico de sequía estival o de congelamiento invernal.

Pérdida de referencia de neutro, sobretensiones dinámicas temporales (TOV) y falla de transformadores

En sistemas con transformadores de potencia conectados en estrella sólidamente aterrizada (Yg), la impedancia homopolar del sistema ($Z_0$) está intrínsecamente ligada a la impedancia del electrodo de puesta a tierra del neutro ($Z_g = R_g + j X_g$).

Bajo una falla a tierra en una de las fases (fase A) en un régimen de suelo degradado donde $R_g \gg X_1$, la relación de secuencias cambia drásticamente: $R_0 / X_1 > 2$ y $X_0 / X_1 > 3$. El sistema deja de operar como un sistema sólidamente puesto a tierra y transiciona dinámicamente a un régimen de neutro aislado o puesto a tierra por alta impedancia. La sobretensión temporal inducida en las fases sanas (fases B y C) se describe vectorialmente como:

Esta condición de sobretensión a frecuencia industrial sostenida ($TOV = \sqrt{3}\ p.u.$) excede la capacidad de tensión máxima continua de operación ($U_c$ / MCOV) de los pararrayos instalados y supera la rigidez dieléctrica del aislamiento menor de transformadores, devanados terciarios y terminales de cables MT/AT, provocando perforaciones dieléctricas destructivas por efecto corona y avalanchas de descargas parciales.

Inoperancia o recalcitrancia de descargadores de sobretensión (Surge Arresters / DPS)

Los descargadores de óxido de zinc ($ZnO$) dependen de una trayectoria de impedancia transitoria ultrabaja a tierra para evacuar la onda de corriente impulsiva ($8/20\ \mu s$ o $10/350\ \mu s$). La impedancia transitoria de impulso $Z_p(t)$ de un electrodo vertical bajo congelamiento severo o suelo seco no se comporta de forma resistiva pura; presenta una fuerte inductancia serie parásita ($Lp$) combinada con una resistencia impulsiva $R_i$ agravada por la falta de portadores iónicos libres:

Cuando el suelo superficial está congelado o seco, el campo eléctrico crítico de ionización del suelo ($E_0 \approx 300\ kV/m$) necesario para formar canales de micro-chispas de plasma que reduzcan la resistencia impulsiva no se alcanza localmente. Como consecuencia, el término $i(t) \cdot R_i(t)$ adquiere valores del orden de cientos de kilovoltios, transfiriendo la totalidad del frente de onda de la sobretensión hacia el aislamiento de barras de subestación, celdas blindadas en $SF _6$ o transformadores protegidos nominalmente.

Falla de coordinación de protecciones de falla a tierra (ANSI 50N/51N, 67N) por elevación de impedancia de retorno

La corriente de cortocircuito monofásica a tierra ($Ik1''$) calculada bajo la norma IEC 60909 se formula analíticamente a través de sus impedancias de secuencia directa ($Z_1$), inversa ($Z_2$) y homopolar ($Z_0$):

Si la resistencia estacional del SPT de la subestación o de las estructuras de la línea se eleva drásticamente por congelamiento o sequedad ($R_g(t) \gg Z_1$), la magnitud de la corriente $Ik1''$ se reduce drásticamente por debajo del umbral de arranque (pickup) de los relés de sobrecorriente de tiempo inverso o definido (ANSI 51N/50N). La falla se convierte en una falla de alta impedancia (Hi-Z), permaneciendo sostenida en el sistema, lo que induce arcos eléctricos intermitentes, ignición de vegetación circundante y la elevación prolongada de potenciales de paso y contacto mortales en el perímetro exterior de las instalaciones.

Matriz Comparativa de Normativa y Límites Críticos de Operación

La siguiente matriz técnica condensa los límites operativos, implicaciones de seguridad y condiciones de falla según los principales estándares internacionales de ingeniería eléctrica frente al deterioro por condiciones ambientales extremas.

Parámetro / Norma Técnica	Límite Estándar Nominal (Condición Base)	Comportamiento Bajo Sequía Severa ($S_w < 0.1$)	Comportamiento Bajo Congelamiento ($T < -10\ ^\circ C$)	Modo de Falla Forense / Consecuencia Operacional
Resistencia de Puesta a Tierra ($R_g$) IEEE Std 81 / IEEE Std 142	Subestaciones AT: $R_g \le 1.0\ \Omega$ Distribución/MT: $R_g \le 5.0\ \Omega$	Incremento de $300\%$ a $1200\%$ ($3.0\ \Omega \le R_g \le 15.0\ \Omega$)	Incremento de $800\%$ a $3500\%$ ($10.0\ \Omega \le R_g \le 80.0\ \Omega$)	Violación del GPR; inoperancia de esquemas de disparo de neutro; flameo inverso (backflashover) en cadenas de aisladores.
Tensión de Contacto Tolerable ($Etouch$) IEEE Std 80 / IEC 61936-1	$Etouch \le \frac{1000 + 1.5 C_s \rho_s}{\sqrt{t_s}}$ Base: $\rho_s \approx 2500\ \Omega\cdot m$ (grava)	Aparente alza superficial en grava seca, pero colapso del perfil subsuperficial: $E_m$ supera el límite de seguridad.	La capa congelada satura el factor de reducción $C_s \to 0$. Caída drástica de la tensión admisible humana.	Fibrilación ventricular y muerte de operadores en maniobras de interruptores o seccionadores en patio AT.
Límite de Disparo ANSI 51N IEC 60255-151 / IEEE C37.112	Sensibilidad ajustada al $10\% - 20\%$ de la corriente nominal de carga o disparo rápido.	$Ik1''$ cae por debajo del umbral de arranque debido a la alta impedancia de lazo homopolar.	La corriente de falla desciende hasta confundirse con los desbalances normales de fase por carga asimétrica.	Falla no despejada; daño térmico irreversible en transformadores de potencia por TOV prolongada; colapso de cables.
Resistencia de Aislamiento y MCOV de DPS IEC 60099-4 / IEEE C62.11	$V_{ residual } \le BIL _{ equipo } / 1.4$ Tensión de corte continuo $U_c \ge 1.05 VLN$	TOV dinámica incrementa la tensión en bornes a $\sqrt{3} VLN$, sobrepasando el TOV capability curve.	La impedancia de dispersión $R_i$ se dispara; el pararrayos no absorbe energía, descargando sobre el devanado primario.	Explosión termomecánica de las envolventes poliméricas o de porcelana de los DPS de óxido de zinc ($ZnO$).
Gradiente de Tensión Perimetral IEC 62305 / NFPA 780	Separación equipotencial con $Vstep < 5\ kV$ en borde de malla a $1\ m$.	Expansión geométrica del cono de potencial hacia el exterior de la subestación.	La capa superficial congelada fuerza la salida de líneas de campo eléctrico concentradas lateralmente.	Descargas disruptivas hacia cercos perimetrales metálicos externos, electrocutando transeúntes o fauna silvestre.

Físico-Química de los Geles y Compuestos de Acondicionamiento Conductivo

Para mitigar la variabilidad climática estacional se utilizan compuestos de mejora electro-química que modifican localmente la interfaz electrodo-suelo, estabilizando el contenido de humedad, optimizando la concentración iónica y reduciendo la impedancia de contacto metal-terreno.

Clasificación de agentes: Bentonita sódica, geles a base de poliacrilamidas y matrices de carbono liofilizadas

Existen diferencias estructurales y electrodinámicas críticas entre los diversos agentes acondicionadores empleados en ingeniería:

• Bentonita Sódica (Montmorillonita Hydratada): Arcilla mineral filosilicatada compuesta por capas octaédricas de alúmina intercaladas entre dos capas tetraédricas de sílice ($2:1$). Posee la capacidad de hincharse hasta 12 a 15 veces su volumen seco original en presencia de agua libre. Su conductividad es predominantemente iónica por intercambio de cationes $Na ^+$. Sin embargo, bajo condiciones de sequía térmica persistente, sufre contracción volumétrica irreversible (agrietamiento reticular), perdiendo el contacto físico-eléctrico directo con el electrodo metálico y aumentando drásticamente la resistencia de contacto.
• Geles de Reticulación Polimérica (Poliacrilamidas e Hidrogeles Iónicos): Polímeros superabsorbentes formados por cadenas hidrofílicas tridimensionales entrecruzadas covalentemente. Retienen soluciones electrolíticas concentradas de sales higroscópicas ($MgSO _4$, $CaCl _2$, silicatos de potasio) mediante enlaces de hidrógeno. No experimentan agrietamiento mecánico volumétrico durante la desorción hídrica y retienen la humedad a tensiones mátricas de hasta $1.5\ MPa$, impidiendo la congelación a temperaturas inferiores a $-18\ ^\circ C$ debido al descenso crioscópico del soluto atrapado en la matriz.
• Cementos Conductivos y Matrices Carbonosas Base Grafito: Compuestos de curado hidráulico a base de cemento Portland dopados con micro-partículas de grafito cristalino desulfurado ($C \ge 99\%$). La conducción es predominantemente electrónica (túnel cuántico / teoría de percolación cuántica) y no iónica, con resistividades intrínsecas secas de $0.02$ a $0.1\ \Omega\cdot m$. No dependen del agua poral para conducir, haciéndolos impermeables a la sequía y al congelamiento estructural, además de proporcionar una barrera física contra la corrosión galvánica.

Higroscopia, capacidad de intercambio catiónico (CIC) y cinéticas de desorción hídrica

La Capacidad de Intercambio Catiónico (CIC) cuantifica la densidad de cargas negativas en la superficie de la matriz coloidal capaces de retener iones intercambiables en equilibrio:

Mientras que un suelo arenoso estándar posee un $CIC \le 5\ meq /100\ g$, un gel polimérico dopado o una bentonita tratada supera los $80 - 120\ meq /100\ g$. La cinética de evaporación y retención de humedad a través del tiempo $S_w(t)$ bajo estrés térmico se modela a través de la ecuación de transporte difusivo no lineal de Richards:

Los geles electroquímicos de última generación modifican localmente la conductividad hidráulica no saturada $K(\theta)$, actuando como una válvula osmótica inversa que previene la migración ascendente de agua por termomigración hacia la superficie evaporativa durante períodos de sequía extrema.

Corrosión galvánica inducida, pH y pasivación de electrodos de cobre y acero cobreado

Un acondicionador químico no debe generar la degradación acelerada del electrodo metálico. El potencial electroquímico de corrosión ($E_{ corr }$) y la densidad de corriente anódica de pasivación ($J_{ corr }$) se calculan a través de las curvas de polarización de Tafel:

Para asegurar una vida útil de diseño superior a 30 años conforme a la norma IEC 62561-7 (Requisitos para compuestos de mejora de tierra):

• El pH del gel curado debe situarse estrictamente en la región de inmunidad / pasivación termodinámica de Pourbaix para el cobre: $7.0 \le pH \le 10.5$.
• La concentración lixiviable de iones sulfato ($SO _4^{2-}$) y cloruro ($Cl ^-$) debe ser menor al $0.05\%$ en peso seco, para evitar el picado por picaduras (pitting corrosion) del recubrimiento de cobre en varillas de acero cobreado (copper-clad steel, grosor mínimo nominal $\ge 254\ \mu m$).
• La resistividad propia del compuesto compactado debe ser menor a $0.2\ \Omega\cdot m$ ($20\ \Omega\cdot cm$) en saturación y menor a $1.0\ \Omega\cdot m$ en condiciones deshidratadas de ensayo de laboratorio acelerado a $105\ ^\circ C$.

Metodología de Ingeniería, Diseño y Simulación Avanzada con Vexten Suite

La integración del análisis numérico avanzado de tierras con la simulación electromecánica de potencia es indispensable para garantizar el cumplimiento normativo en todo el ciclo estacional.

Inversión de sondeos Wenner/Schlumberger y parametrización de capas en Vexten Grounding Engine

El primer paso metodológico reside en la ejecución rigurosa de ensayos de sondeo eléctrico vertical (SEV) mediante la configuración tetraelectródica de Wenner según IEEE Std 81, con espaciamientos crecientes $a = [1, 2, 4, 8, 16, 32]\ m$. La resistencia medida $R(a)$ arroja una resistividad aparente:

En el módulo de inversión geofísica no lineal de Vexten Grounding Engine, se implementa el algoritmo de optimización Gauss-Newton con regularización de Tikhonov para resolver el problema inverso y obtener la estratigrafía multicapa de diseño:

Donde $\mathbf{p} = [\rho_1, \rho_2, \dots, \rho_n, h_1, H₂, \dots, h_{n-1}]^T$ representa el vector de parámetros geoeléctricos, $\mathbf{L}$ es la matriz de suavizado diferencial de primer orden, y $\lambda$ es el multiplicador de Lagrange del parámetro de regularización.

Con Vexten Suite, el ingeniero parametriza dos escenarios envolventes (worst-case envelopes):

• Escenario Estival / Sequía Crítica: Modificación de la capa 1 ($h_1 = 1.5\ m$) con un factor multiplicador hidrológico $\kappa_{ dry } \in [3.0, 8.0]$.
• Escenario Invernal / Congelamiento Extremo: Modificación de la capa superficial ($h_1 = 0.8\ m$) con resistividad criogénica $\rho_1 \ge 10\,000\ \Omega\cdot m$.

Co-simulación de cortocircuito asimétrico IEC 60909 / IEEE 141 acoplado a la impedancia estacional de malla

A diferencia de los flujos de trabajo desacoplados tradicionales, el solucionador matricial de cortocircuito de Vexten Short-Circuit IEC 60909 ejecuta una co-simulación acoplada iterativa con el cálculo de malla:

El valor de $R_g(t)$ calculado por el módulo electrostático retroalimenta dinámicamente la matriz nodal de admitancias de secuencia cero. De este modo, el software recalcula con absoluta precisión el factor de división de corriente real ($S_f$), integrando la contribución real de los hilos de guarda, pantallas de cables aislados dimensionados según IEC 60287 / NEC 310, y las corrientes de retorno por el terreno profundo.

Optimización topológica de electrodos profundos y zanjas mejoradas con gel conductivo

Cuando la simulación en terreno natural arroja una violación de los potenciales máximos admisibles de contacto ($E_m > E_{ touch }$), Vexten Grounding Engine permite formular e iterar la estrategia de mitigación mediante acondicionamiento electroquímico perimetral y pozos profundos encapsulados.

El radio equivalente hidrodinámico-eléctrico ($req$) de una pica vertical de radio metálico $r_0$ alojada en una perforación de radio $r_b$ e inyectada con gel conductivo de resistividad $\rho_{ gel }$ dentro de un suelo de resistividad local $\rho_{ suelo }$ se modela cuantitativamente como:

Dado que $\rho_{ gel } \ll \rho_{ suelo }$ (por ejemplo, $0.1\ \Omega\cdot m$ frente a $500\ \Omega\cdot m$), el término exponencial tiende a $e^0 = 1$, provocando que $req \to r_b$. De forma analítica, el pozo de acondicionamiento transforma el radio de inyección de corriente desde el radio del metal ($r_0 \approx 9.5\ mm$) al radio total de la excavación ($r_b \approx 75 - 150\ mm$), reduciendo la resistencia intrínseca de contacto del electrodo en más de un $40\%$ a $60\%$ e independizándolo radicalmente de la desecación superficial.

Al forzar mecánicamente la disipación hacia la capa inferior protegida ($\rho_2$) a través del electrodo encapsulado en matriz coloidal, se asegura la estabilidad del sistema frente a fluctuaciones climáticas, garantizando una operación confiable, la correcta coordinación de las protecciones de sobrecorriente y la seguridad total del personal de campo durante toda la vida útil de la subestación.