Resistividad de Terreno con Método de Wenner y Estratificación Bicapa según IEEE 81 e IEEE 80

Fundamentos Electromagnéticos de la Conducción en Medios Geoeléctricos y Método de Wenner

La caracterización geoeléctrica del subsuelo constituye el pilar fundamental e ineludible para el diseño ingenieril de sistemas de puesta a tierra (SPT) en subestaciones de media, alta y extra alta tensión. A diferencia de los conductores metálicos puros donde el transporte de carga está regido exclusivamente por electrones libres bajo la ley de Ohm microscópica, el suelo es un medio heterogéneo, anisótropo y multifásico. La conducción eléctrica en el terreno se produce de manera predominante a través de un mecanismo electrolítico dentro del fluido poroso intersticial, gobernado por la salinidad, humedad volumétrica, temperatura, compacidad y composición mineralógica de la matriz sólida.

A nivel macroscópico, el comportamiento del campo electrostático cuasiestacionario en el semiespacio conductor infinito inferior ($z \ge 0$) se describe mediante las ecuaciones de Maxwell para medios conductores en régimen estacionario, donde la densidad de corriente $\mathbf{J}$ y el potencial escalar eléctrico $\Phi$ satisfacen:

Para un medio homogéneo e isótropo con conductividad escalar constante $\sigma = 1/\rho$, la ecuación diferencial gobernante se reduce a la ecuación de Laplace estándar $\nabla^2 \Phi = 0$, sujeta a las condiciones de frontera de Neumann en la interfaz suelo-aire ($z = 0$):

Derivación del Potencial para Fuentes Puntuales y Arreglo Tetrapolar Simétrico

Considerando una fuente puntual de corriente $I$ inyectada en la superficie de un semiespacio homogéneo ($z = 0$), la simetría esférica en el semiespacio impone que las líneas de corriente diverjan radialmente en un ángulo sólido de $2\pi$ estereorradianes. El potencial eléctrico $\Phi(r)$ a una distancia radial $r$ de la fuente se obtiene integrando la ley de Ohm en coordenadas esféricas:

En el método tetrapolar propuesto por Frank Wenner (1915), se disponen cuatro electrodos colineales equidistantes con un espaciamiento interelectródico $a$. Los dos electrodos exteriores ($C_1$ y $C_2$) actúan como inyectores de corriente ($+I$ y $-I$), mientras que los dos electrodos interiores ($P_1$ y $P_2$) actúan como sondas de medición del potencial diferencial $\Delta V = \Phi(P_1) - \Phi(P_2)$.

Aplicando el principio de superposición lineal para las fuentes de corriente situadas en $C_1$ ($x = 0$) y $C_2$ ($x = 3a$), los potenciales en los electrodos de potencial ubicados en $P_1$ ($x = a$) y $P_2$ ($x = 2a$) resultan:

La diferencia de potencial neta medida por el voltímetro de alta impedancia viene dada por:

Despejando el parámetro de resistividad, se define la resistividad aparente $\rho_a$ para una separación interelectródica $a$:

donde $R_m = \Delta V / I$ es la resistencia de transferencia mutua aparente leída por el telurómetro o teluróhmetro digital de baja frecuencia.

Corrección por Profundidad Finita de Enterramiento de Picas

En condiciones de campo reales, las picas de prueba se entierran hasta una profundidad finita $b$. Si la relación $b/a$ no es despreciable ($b > 0.1a$), la aproximación puntual de superficie introduce un error sistemático relevante. Considerando electrodos cilíndricos verticales de longitud $b$, la formulación matemática rigurosa corregida por el efecto de penetración de las picas viene expresada por:

Cuando $b \to 0$, el denominador de la fracción converge exactamente al valor de $2$, recuperando la ecuación clásica de Wenner $\rho_a = 2\pi a R_m$. En la práctica según la norma IEEE Std 81, se exige mantener $b \le 0.05a$ a fin de despreciar la corrección geométrica sin exceder un error analítico del $1\%$.

Fuentes de Incertidumbre Metrológica, Interferencias y Acoplamientos en Campo

La obtención de perfiles de resistividad aparente de alta fidelidad está sujeta a diversas perturbaciones electromagnéticas y no ideales del entorno que degradan severamente la función de transferencia del subsuelo. Los factores más críticos se resumen a continuación:

Fuente de Perturbación	Mecanismo Físico / Electromagnético	Efecto sobre $\rho_a$	Estrategia de Mitigación Normativa
Acoplamiento Inductivo entre Cables	Inducción mutua $M \cdot \frac{ d i}{ d t}$ entre los lazos de corriente $C_1-C_2$ y de potencial $P_1-P_2$.	Sobreestimación artificial de $\Delta V$, con desfase aparente a frecuencias superiores a 100 Hz.	Separación física transversal de cables de corriente y tensión; utilización de cables coaxiales o bifilares apantallados; inyección a frecuencias sub-armónicas (e.g., 55 Hz, 94 Hz, 105 Hz, 128 Hz).
Resistencias de Contacto Elevadas ($R_c$)	Impedancia de interfaz electrodo-suelo excesiva debida a grava seca, rocas o escasa humedad superficial.	Saturación de la etapa de entrada del voltímetro, atenuación drástica de $I$, ruido térmico Johnson-Nyquist y lecturas erráticas.	Humectación de picas con solución salina o bentonita; inserción de múltiples picas en paralelo en los terminales de corriente.
Corrientes Parásitas y Telúricas	Retornos de neutro desbalanceado, tracción eléctrica en CC/CA, gradientes galvánicos y actividad geomagnética.	Distorsión de la señal diferencial en $P_1-P_2$, deriva de cero en CC y fluctuaciones caóticas de baja frecuencia.	Filtrado síncrono por modulación de fase (detección homodina/Lock-in); inyección de corriente cuadrada con inversión de polaridad periódica.
Estructuras Metálicas Enterradas	Tuberías de gas/agua, mallas de tierra existentes y blindajes de cables que actúan como caminos de baja impedancia.	Subestimación crítica de la resistividad en espaciamientos grandes ($a > 10 m$), distorsionando la estratificación.	Trazado de perfiles ortogonales y diagonales; distanciamiento de infraestructuras metálicas enterradas mayor a $3a_{\max}$.

Teoría Matemática de la Estratificación Geoeléctrica Bicapa

En la inmensa mayoría de las formaciones geológicas de interés para la ingeniería eléctrica, el subsuelo no puede ser modelado con precisión mediante un semiespacio homogéneo. La sedimentación, el grado de compactación y el contenido freático generan gradientes verticales marcados. El modelo canónico fundamental es el modelo estratificado bicapa, compuesto por una primera capa superficial de espesor finito $h$ y resistividad intrínseca $\rho_1$, superpuesta sobre un semi-espacio infinito inferior de resistividad intrínseca $\rho_2$.

Resolución de la Ecuación de Laplace con Condiciones de Borde

Para una fuente puntual que inyecta una corriente $I$ en el origen $(0,0,0)$ sobre la superficie ($z=0$), el campo escalar de potencial en coordenadas cilíndricas $(\phi(r,z))$ satisface la ecuación de Laplace en cada medio homogéneo:

Aplicando el método de separación de variables mediante la transformada integral de Fourier-Bessel (orden cero), la solución general para el potencial en la primera capa ($0 \le z \le h$) y en la segunda capa ($z \ge h$) se formula como:

donde $J_0(\lambda r)$ es la función de Bessel de primera especie y orden cero, y los coeficientes $A(\lambda)$, $B(\lambda)$ y $C(\lambda)$ se determinan imponiendo las condiciones de frontera de continuidad del potencial y de la componente normal de la densidad de corriente:

1. Condición en la superficie aire-suelo ($z = 0$): $\left. \frac{\partial \Phi_1}{\partial z} \right|_{z=0} = 0 \implies A(\lambda) = B(\lambda)$
2. Continuidad del potencial en la interfaz ($z = h$): $\Phi_1(r, h) = \Phi_2(r, h)$
3. Continuidad de la densidad de corriente normal en la interfaz ($z = h$): $\frac{1}{\rho_1} \left. \frac{\partial \Phi_1}{\partial z} \right|_{z=h} = \frac{1}{\rho_2} \left. \frac{\partial \Phi_2}{\partial z} \right|_{z=h}$
4. Condición de regularidad en el infinito: $\lim_{z \to \infty} \Phi_2(r, z) = 0$

Definiendo el coeficiente de reflexión de resistividad $k$ como:

Al resolver el sistema algebraico lineal para $A(\lambda)$, se obtiene:

Expandiendo el denominador mediante su desarrollo en serie geométrica infinita $(1 - k e^{-2\lambda h})^{-1} = \sum_{n=1}^{\infty} k^n e^{-2n\lambda h}$, la función de potencial sobre la superficie libre ($z = 0$) adopta la forma:

Haciendo uso de la identidad integral de Lipschitz-Weber:

Se arriba a la expresión analítica cerrada del potencial electrostático superficial debida al método de imágenes infinitas:

Ecuación Analítica de la Resistividad Aparente de Wenner

Aplicando la fórmula general del potencial de superficie al arreglo de cuatro puntos de Wenner, el potencial medido entre $P_1$ y $P_2$ ante las corrientes $+I$ en $C_1$ y $-I$ en $C_2$ da lugar a la formulación maestra de resistividad aparente para un suelo bicapa:

Esta serie alternada converge absolutamente para todo $|k| < 1$ y permite calcular teóricamente el perfil de $\rho_a(a)$ en función de la separación interelectródica $a$. Presenta dos asíntotas físicas fundamentales:

Algoritmos de Inversión Geoeléctrica y Ajuste de Curvas No Lineales

En la práctica de la ingeniería forense y de diseño, el problema consiste en determinar el vector de parámetros desconocidos del terreno $\mathbf{p} = [\rho_1, \rho_2, h]^T$ a partir de un conjunto discreto de $M$ mediciones experimentales $\{(\ln a_i, \ln \rho_{a,i}^{med})\}_{i=1}^M$. Este problema inverso es intrínsecamente no lineal y mal condicionado.

Formulación del Problema Inverso por Mínimos Cuadrados No Lineales

Se define la función objetivo de error cuadrático relativo normalizado $\chi^2(\mathbf{p})$ en escala logarítmica para balancear las sensibilidades sobre órdenes de magnitud dispares:

La optimización se lleva a cabo mediante el Algoritmo de Levenberg-Marquardt (Damped Gauss-Newton), el cual interpola adaptativamente entre el método de descenso de gradiente y el método de Gauss-Newton linealizado:

donde $\mathbf{J} \in \mathbb{R}^{M \times 3}$ es la matriz Jacobiana de sensibilidades parciales:

y $\mu_k \ge 0$ es el parámetro de amortiguamiento (damping factor) ajustado en cada iteración en función de la métrica de ganancia de reducción residual.

Matriz de Sensibilidad Analítica (Derivadas Parciales de Frechet)

A diferencia del cálculo numérico por diferencias finitas que introduce errores por redondeo, el cálculo analítico del gradiente acelera la convergencia asintótica:

donde $\frac{\partial k}{\partial \rho_1} = \frac{-2\rho_2}{(\rho_1 + \rho_2)^2}$ y $\frac{\partial k}{\partial \rho_2} = \frac{2\rho_1}{(\rho_1 + \rho_2)^2}$.

El criterio de convergencia estricto exige que el error cuadrático medio relativo residual (RMS) satisfaga:

Impacto de la Estratificación Bicapa en el Diseño de Mallas de Puesta a Tierra (IEEE Std 80 / IEC 60479)

Omitir la estratificación del terreno y asumir un modelo de resistividad constante equivalente promediado aritméticamente genera errores catastróficos en el dimensionamiento de conductores de cobre, longitud de picas y evaluación de tensiones de paso y contacto en subestaciones.

Resistencia de Puesta a Tierra ($R_g$) bajo Modelo Bicapa

De acuerdo con las formulaciones extendidas de Sverak y Schwarz en medios bicapa, la resistencia total de una malla de puesta a tierra cuadrada o rectangular con electrodos verticales depende críticamente de la profundidad de la primera capa $h$ y del contraste de resistividades $\rho_2 / \rho_1$:

donde $\rho_{eq}$ es la resistividad aparente equivalente vista por la geometría de la malla:

donde $Ltotal = L_c + L_v$ es la longitud combinada de conductores horizontales ($L_c$) y picas verticales ($L_v$), y $\psi(h)$ es la función de acoplamiento de profundidad.

Potenciales de Paso y Contacto: Comparativa Fenomenológica de Escenarios

Parámetro de Diseño / Seguridad	Escenario A: Terreno Favorable ($\rho_1 > \rho_2$, $k < 0$)	Escenario B: Terreno Crítico ($\rho_1 < \rho_2$, $k > 0$)	Implicación Crítica de Seguridad
Distribución de Corriente de Falla	La corriente inyectada tiende a drenar verticalmente hacia el estrato inferior profundo de baja resistividad.	La corriente es repelida por el estrato inferior altamente resistivo (roca), canalizándose horizontalmente por la capa superficial.	En el Escenario B, la densidad de corriente superficial $J_s$ aumenta de manera exponencial.
Potencial de Malla ($V_m$) y de Contacto ($Etouch$)	Gradientes de potencial superficial planos y atenuados. $Etouch$ se mantiene dentro de los límites admisibles con conductores estándar.	Severa elevación del Gradiente de Potencial de Superficie (GPR) y crestas de tensión extremadamente pronunciadas sobre la malla.	En el Escenario B, se supera frecuentemente el límite de fibrilación ventricular según IEEE 80 / IEC 60479.
Efectividad de las Picas de Puesta a Tierra	Altamente efectivas si penetran la capa profunda ($l_v > h$), reduciendo drásticamente $R_g$.	Poco efectivas; las picas apenas disipan corriente adicional en el estrato de roca profundo de alta resistividad.	En el Escenario B, se requiere densificar la cuadrícula horizontal y expandir la superficie perimetral.
Factor de Reducción de la Capa de Grava ($C_s$)	El espesor superficial de piedra partida ($h_s$) proporciona una mitigación estándar ($C_s \approx 0.7 - 0.85$).	La interacción de $\rho_1$ baja con la roca superficial y $\rho_2$ infinita altera la reflexión del lecho de grava.	Requiere recalcular $C_s$ mediante series infinitas de doble reflexión (grava-capa 1-capa 2).

Cálculo Riguroso del Factor de Reflexión de Capa Superficial ($C_s$)

El factor de corrección $C_s$ para la capa de material de alta resistividad superficial (grava/chancado $\rho_s \approx 3000\,\Omega\cdot m$) de espesor $h_s$ instalada sobre un suelo bicapa no puede computarse con la clásica fórmula empírica de Sunde para suelo homogéneo. La formulación generalizada mediante expansión en serie de coeficientes combinados resulta:

donde $k_s = (\rho_1 - \rho_s)/(\rho_1 + \rho_s)$ y $r_d = 0.08 m$ representa el radio equivalente del pie humano modelado como un disco circular metálico de radio normalizado.

Análisis Forense de Fallas Asociadas a Caracterización Geoeléctrica Deficiente

Los errores metodológicos en la campaña de prospección de resistividad o la simplificación a suelo homogéneo cuando subyace un estrato altamente reflexivo ($k \to +1$) han derivado históricamente en fallas catastróficas en sistemas de potencia. A continuación se presentan mecanismos de falla e incidentes forenses documentados:

Disrupción Dieléctrica en Transformadores de Medición y Servicios Auxiliares por Elevación Crítica de GPR

En una subestación de 230 kV construida sobre una formación de arenisca superficial ($\rho_1 = 180\,\Omega\cdot m$, $h = 1.8 m$) sobre un basamento granítico masivo ($\rho_2 = 4200\,\Omega\cdot m$, $k = +0.918$), el diseño original asumió un suelo homogéneo promedio de $350\,\Omega\cdot m$. Durante un cortocircuito monofásico a tierra con corriente simétrica $I_f = 22 kA$, la resistencia real de la malla fue de $2.85\,\Omega$ (frente a los $0.68\,\Omega$ proyectados bajo el supuesto homogéneo).

La densa corriente que no pudo penetrar el granito profundo retornó lateralmente por las trincheras de cables de control, elevando el potencial local de tierra de los gabinetes exteriores por encima de la rigidez dieléctrica de los devanados secundarios de los transformadores de corriente (TC) y transformadores de potencial (TP). La diferencia de potencial de modo común perforó el aislamiento galvánico de 10 kV de los relés de protección IED basados en microprocesador en la sala de control, quemando los módulos de entradas y salidas analógicas e inhibiendo el disparo de respaldo temporizado (ANSI 50/51N), lo que extendió el despeje de la falla de 80 ms a 1.2 segundos y causó la explosión destructiva del transformador de potencia principal.

Degradación Térmica y Fusión de Pantallas de Cables Subterráneos de Alta Tensión

En tendidos de cables aislados XLPE de 115 kV con pantallas metálicas puestas a tierra en ambos extremos (solid bonding), la asunción errónea de una baja resistividad en estratos profundos distorsiona la impedancia mutua de bucle tierra-pantalla $\underline{Z}_{m(c-g)}$ según las fórmulas de Carson-Clem:

Al encontrarse un estrato profundo de alta resistividad ($\rho_2 \gg \rho_1$), la profundidad de retorno equivalente $D_e$ se incrementó en más del 300%. Esto incrementó drásticamente la reactancia inductiva del retorno por tierra, forzando a que la corriente de desbalance y de falla circulara casi en su totalidad a través de las pantallas metálicas de los cables de potencia. La densidad de corriente térmica $Jscreen$ superó el límite admisible del cobre según IEC 60287 e IEC 60949:

provocando la delaminación térmica de la cubierta exterior semiconductora, el reblandecimiento del polietileno reticulado (XLPE) por sobretemperatura superior a 250 °C y la consecuente falla dieléctrica fase-tierra por perforación de árboles de agua acelerados electrotérmicamente.

Flujo Metodológico de Simulación y Diseño Avanzado en Vexten Suite

El ecosistema de ingeniería integral Vexten Suite incorpora un motor determinístico y multifísico de inversión geoeléctrica acoplado a módulos de cortocircuito (IEC 60909 / IEEE 141), dimensionamiento de conductores y cables (IEC 60287 / NEC) y optimización de mallas de puesta a tierra (IEEE 80 / IEC 60479). A continuación se detalla la integración algorítmica y el protocolo metodológico aplicado en proyectos de infraestructura crítica:

Ejemplo Práctico de Simulación e Inversión Numérica en Vexten Grounding Engine

Se realiza una campaña de prospección geoeléctrica tetrapolar con método de Wenner en el emplazamiento de una nueva subestación industrial de 138/13.8 kV. La adquisición de datos arrojó los siguientes valores de campo:

Espaciamiento $a$ [m]	Resistencia Medida $R_m$ [$\Omega$]	Resistividad Aparente Experimental $\rho_a^{med}$ [$\Omega\cdot m$]	Resistividad Aparente Calculada $\rho_a^{calc}$ [$\Omega\cdot m$]	Error Residual Relativo [\%]
1.0	68.435	430.00	428.12	-0.44
2.0	32.786	412.00	418.54	+1.59
4.0	14.928	375.20	379.80	+1.23
8.0	5.431	273.00	269.45	-1.30
16.0	1.581	159.00	157.10	-1.19
32.0	0.507	102.00	103.85	+1.81

Al procesar la matriz de datos en el módulo de optimización inversa no lineal de Vexten Grounding, el algoritmo de Levenberg-Marquardt converge en 7 iteraciones con un error global RMS del 1.32%, determinando el siguiente vector de estado estratificado:

Análisis de Ingeniería y Optimización del Diseño del SPT

Dado que el coeficiente de reflexión $k$ es fuertemente negativo ($k = -0.6812$), el estrato profundo presenta una resistividad más de cinco veces inferior a la de la capa superficial ($\rho_2 \ll \rho_1$). La optimización computacional en Vexten Grounding explota este fenómeno físico mediante las siguientes directrices ejecutivas:

1. Instalación Estratégica de Picas Profundas: Se especifica la hincadura de picas de acero cobreado (copperweld) con una longitud mínima $l_v = 6.0 m$ en las cuatro esquinas y nodos perimetrales críticos. Al sobrepasar la profundidad de la interfaz ($h = 4.35 m$), los electrodos verticales penetran $1.65 m$ en el estrato de baja resistividad ($\rho_2 = 82.50\,\Omega\cdot m$). Esto reduce la resistencia total del sistema $R_g$ en un 42.7% en comparación con una malla puramente horizontal de igual longitud total de conductor.
2. Descongestionamiento de Conductor Horizontal: Gracias a la alta admitancia de la capa inferior, la separación entre conductores horizontales de la retícula se puede incrementar de un reticulado denso de $5 m \times 5 m$ a $10 m \times 10 m$, logrando un ahorro significativo de cobre sin exceder las tensiones de paso y de contacto límite estipuladas por la norma IEEE Std 80.
3. Integración con el Módulo de Cortocircuito de Vexten Suite (IEC 60909): La impedancia real del terreno estratificado se exporta automáticamente hacia el motor de cortocircuito de Vexten Suite, permitiendo recalcular de manera acoplada los factores de división de corriente por cables de guarda y blindajes de fase ($S_f$), asegurando que la corriente que efectivamente retorna por el terreno ($I_g = S_f \cdot I_k''$) refleje con exactitud la física multifilar del sistema de transmisión y distribución circundante.

La adopción de este flujo metodológico integral garantiza la máxima rigurosidad analítica, seguridad dieléctrica del personal y de los activos electrotérmicos, así como una optimización económica y operativa en la construcción de sistemas de puesta a tierra para subestaciones y plantas de generación.