Medição de resistividade aparente do solo pelo método de Wenner de 4 eletrodos e estratificação em duas camadas conforme IEEE 81 e IEEE 80

Fundamentos Eletromagnéticos da Condução em Meios Geoelétricos e Método de Wenner

A caracterização geoelétrica do subsolo constitui o pilar fundamental e inquestionável para o projeto de engenharia de sistemas de aterramento em subestações de média, alta e extra-alta tensão. Ao contrário dos condutores metálicos puros, onde o transporte de carga é regido exclusivamente por elétrons livres sob a lei de Ohm microscópica, o solo é um meio heterogêneo, anisotrópico e multifásico. A condução elétrica no terreno ocorre de maneira predominante por meio de um mecanismo eletrolítico no fluido poroso intersticial, governado pela salinidade, umidade volumétrica, temperatura, compactação e composição mineralógica da matriz sólida.

Em nível macroscópico, o comportamento do campo eletrostático quase-estacionário no semiespaço condutor infinito inferior ($z \ge 0$) é descrito pelas equações de Maxwell para meios condutores em regime estacionário, onde a densidade de corrente $\mathbf{J}$ e o potencial escalar elétrico $\Phi$ satisfazem:

Para um meio homogêneo e isotrópico com condutividade escalar constante $\sigma = 1/\rho$, a equação diferencial governante reduz-se à equação de Laplace padrão $\nabla^2 \Phi = 0$, sujeita às condições de contorno de Neumann na interface solo-ar ($z = 0$):

Dedução do Potencial para Fontes Pontuais e Arranjo Tetrapolar Simétrico

Considerando uma fonte pontual de corrente $I$ injetada na superfície de um semiespaço homogêneo ($z = 0$), a simetria esférica no semiespaço impõe que as linhas de corrente divirjam radialmente em um ângulo sólido de $2\pi$ esterrorradianos. O potencial elétrico $\Phi(r)$ a uma distância radial $r$ da fonte é obtido integrando a lei de Ohm em coordenadas esféricas:

No método tetrapolar proposto por Frank Wenner (1915), dispõem-se quatro eletrodos colineares equidistantes com um espaçamento intereletródico $a$. Os dois eletrodos externos ($C_1$ e $C_2$) atuam como injetores de corrente ($+I$ e $-I$), enquanto os dois eletrodos internos ($P_1$ e $P_2$) atuam como sondas de medição do potencial diferencial $\Delta V = \Phi(P_1) - \Phi(P_2)$.

Aplicando o princípio da superposição linear para as fontes de corrente situadas em $C_1$ ($x = 0$) e $C_2$ ($x = 3a$), os potenciais nos eletrodos de potencial localizados em $P_1$ ($x = a$) e $P_2$ ($x = 2a$) resultam em:

A diferença de potencial líquida medida pelo voltímetro de alta impedância é dada por:

Isolando o parâmetro de resistividade, define-se a resistividade aparente $\rho_a$ para uma separação intereletródica $a$:

onde $R_m = \Delta V / I$ é a resistência de transferência mútua aparente lida pelo telurômetro digital de baixa frequência.

Correção por Profundidade Finita de Cravação das Hastes

Em condições de campo reais, as hastes de teste são cravadas até uma profundidade finita $b$. Se a relação $b/a$ não for desprezível ($b > 0.1a$), a aproximação pontual de superfície introduz um erro sistemático relevante. Considerando eletrodos cilíndricos verticais de comprimento $b$, a formulação matemática rigorosa corrigida pelo efeito de penetração das hastes é expressa por:

Quando $b \to 0$, o denominador da fração converge exatamente para o valor de $2$, recuperando a equação clássica de Wenner $\rho_a = 2\pi a R_m$. Na prática, segundo a norma IEEE Std 81, exige-se manter $b \le 0.05a$ a fim de desprezar a correção geométrica sem exceder um erro analítico de $1\%$.

Fontes de Incerteza Metrológica, Interferências e Acoplamentos no Campo

A obtenção de perfis de resistividade aparente de alta fidelidade está sujeita a diversas perturbações eletromagnéticas e não idealidades do entorno que degradam severamente a função de transferência do subsolo. Os fatores mais críticos estão resumidos a seguir:

Fonte de Perturbação	Mecanismo Físico / Eletromagnético	Efeito sobre $\rho_a$	Estratégia de Mitigação Normativa
Acoplamento Indutivo entre Cabos	Indução mútua $M \cdot \frac{ d i}{ d t}$ entre os laços de corrente $C_1-C_2$ e de potencial $P_1-P_2$.	Superestimativa artificial de $\Delta V$, com desfasamento aparente em frequências superiores a 100 Hz.	Separação física transversal de cabos de corrente e tensão; utilização de cabos coaxiais ou bifilares blindados; injeção em frequências sub-harmônicas (ex.: 55 Hz, 94 Hz, 105 Hz, 128 Hz).
Resistências de Contato Elevadas ($R_c$)	Impedância de interface eletrodo-solo excessiva devida a brita seca, rochas ou escassa umidade superficial.	Saturação do estágio de entrada do voltímetro, atenuação drástica de $I$, ruído térmico de Johnson-Nyquist e leituras erráticas.	Umidificação das hastes com solução salina ou bentonita; inserção de múltiplas hastes em paralelo nos terminais de corrente.
Correntes Parasitas e Telúricas	Retornos de neutro desbalanceado, tração elétrica em CC/CA, gradientes galvânicos e atividade geomagnética.	Distorção do sinal diferencial em $P_1-P_2$, deriva de zero em CC e flutuações caóticas de baixa frequência.	Filtragem síncrona por modulação de fase (detecção homódina/Lock-in); injeção de corrente quadrada com inversão periódica de polaridade.
Estruturas Metálicas Enterradas	Tubulações de gás/água, malhas de aterramento existentes e blindagens de cabos que atuam como caminhos de baixa impedância.	Subestimativa crítica da resistividade em grandes espaçamentos ($a > 10 m$), distorcendo a estratificação.	Traçado de perfis ortogonais e diagonais; distanciamento de infraestruturas metálicas enterradas superior a $3a_{\max}$.

Teoria Matemática da Estratificação Geoelétrica em Duas Camadas (Bicamada)

Na imensa maioria das formações geológicas de interesse para a engenharia elétrica, o subsolo não pode ser modelado com precisão por meio de um semiespaço homogêneo. A sedimentação, o grau de compactação e o nível freático geram gradientes verticais acentuados. O modelo canônico fundamental é o modelo estratificado em duas camadas (bicamada), composto por uma primeira camada superficial de espessura finita $h$ e resistividade intrínseca $\rho_1$, sobreposta a um semiespaço infinito inferior de resistividade intrínseca $\rho_2$.

Resolução da Equação de Laplace com Condições de Contorno

Para uma fonte pontual que injeta uma corrente $I$ na origem $(0,0,0)$ sobre a superfície ($z=0$), o campo escalar de potencial em coordenadas cilíndricas $(\phi(r,z))$ satisfaz a equação de Laplace em cada meio homogêneo:

Aplicando o método de separação de variáveis mediante a transformada integral de Fourier-Bessel (ordem zero), a solução geral para o potencial na primeira camada ($0 \le z \le h$) e na segunda camada ($z \ge h$) é formulada como:

onde $J_0(\lambda r)$ é a função de Bessel de primeira espécie e ordem zero, e os coeficientes $A(\lambda)$, $B(\lambda)$ e $C(\lambda)$ são determinados impondo as condições de contorno de continuidade do potencial e da componente normal da densidade de corrente:

1. Condição na superfície ar-solo ($z = 0$): $\left. \frac{\partial \Phi_1}{\partial z} \right|_{z=0} = 0 \implies A(\lambda) = B(\lambda)$
2. Continuidade do potencial na interface ($z = h$): $\Phi_1(r, h) = \Phi_2(r, h)$
3. Continuidade da densidade de corrente normal na interface ($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. Condição de regularidade no infinito: $\lim_{z \to \infty} \Phi_2(r, z) = 0$

Definindo o coeficiente de reflexão de resistividade $k$ como:

Ao resolver o sistema algébrico linear para $A(\lambda)$, obtém-se:

Expandindo o denominador por meio de seu desenvolvimento em série geométrica infinita $(1 - k e^{-2\lambda h})^{-1} = \sum_{n=1}^{\infty} k^n e^{-2n\lambda h}$, a função de potencial sobre a superfície livre ($z = 0$) assume a forma:

Fazendo uso da identidade integral de Lipschitz-Weber:

Chega-se à expressão analítica fechada do potencial eletrostático superficial devida ao método das imagens infinitas:

Equação Analítica da Resistividade Aparente de Wenner

Aplicando a fórmula geral do potencial de superfície ao arranjo de quatro pontos de Wenner, o potencial medido entre $P_1$ e $P_2$ perante as correntes $+I$ em $C_1$ e $-I$ em $C_2$ dá origem à formulação mestra da resistividade aparente para um solo bicamada:

Esta série alternada converge absolutamente para todo $|k| < 1$ e permite calcular teoricamente o perfil de $\rho_a(a)$ em função da separação intereletródica $a$. Apresenta duas assíntotas físicas fundamentais:

Algoritmos de Inversão Geoelétrica e Ajuste de Curvas Não Lineares

Na prática da engenharia forense e de projeto, o problema consiste em determinar o vetor de parâmetros desconhecidos do terreno $\mathbf{p} = [\rho_1, \rho_2, h]^T$ a partir de um conjunto discreto de $M$ medições experimentais $\{(\ln a_i, \ln \rho_{a,i}^{med})\}_{i=1}^M$. Este problema inverso é intrinsecamente não linear e malcondicionado.

Formulação do Problema Inverso por Mínimos Quadrados Não Lineares

Define-se a função objetivo do erro quadrático relativo normalizado $\chi^2(\mathbf{p})$ em escala logarítmica para balancear as sensibilidades sobre ordens de grandeza díspares:

A otimização é realizada por meio do Algoritmo de Levenberg-Marquardt (Gauss-Newton Amortecido), o qual interpola adaptativamente entre o método do gradiente descendente e o método de Gauss-Newton linearizado:

onde $\mathbf{J} \in \mathbb{R}^{M \times 3}$ é a matriz Jacobiana de sensibilidades parciais:

e $\mu_k \ge 0$ é o parâmetro de amortecimento (damping factor) ajustado em cada iteração em função da métrica de ganho de redução residual.

Matriz de Sensibilidade Analítica (Derivadas Parciais de Fréchet)

A diferença do cálculo numérico por diferenças finitas, que introduz erros de arredondamento, o cálculo analítico do gradiente acelera a convergência assintótica:

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

O critério de convergência estricto exige que o erro quadrático médio relativo residual (RMS) satisfaça:

Impacto da Estratificação em Duas Camadas no Projeto de Malhas de Aterramento (IEEE Std 80 / IEC 60479)

Omitir a estratificação do terreno e assumir um modelo de resistividade constante equivalente ponderado pela média aritmética gera erros catastróficos no dimensionamento de condutores de cobre, comprimento de hastes e avaliação das tensões de passo e toque em subestações.

Resistência de Aterramento ($R_g$) sob Modelo Bicamada

De acordo com as formulações estendidas de Sverak e Schwarz em meios bicamada, a resistência total de uma malha de aterramento quadrada ou retangular com eletrodos verticais depende criticamente da profundidade da primeira camada $h$ e do contraste de resistividades $\rho_2 / \rho_1$:

onde $\rho_{eq}$ é a resistividade aparente equivalente vista pela geometria da malha:

onde $Ltotal = L_c + L_v$ é o comprimento combinado de condutores horizontais ($L_c$) e hastes verticais ($L_v$), e $\psi(h)$ é a função de acoplamento de profundidade.

Potenciais de Passo e Toque: Comparação Fenomenológica de Cenários

Parâmetro de Projeto / Segurança	Cenário A: Terreno Favorável ($\rho_1 > \rho_2$, $k < 0$)	Cenário B: Terreno Crítico ($\rho_1 < \rho_2$, $k > 0$)	Implicação Crítica de Segurança
Distribuição da Corrente de Falta	A corrente injetada tende a drenar verticalmente em direção ao estrato inferior profundo de baixa resistividade.	A corrente é repelida pelo estrato inferior altamente resistivo (rocha), canalizando-se horizontalmente pela camada superficial.	No Cenário B, a densidade de corrente superficial $J_s$ aumenta de maneira exponencial.
Potencial de Malha ($V_m$) e de Toque ($Etouch$)	Gradientes de potencial superficial planos e atenuados. $Etouch$ mantém-se dentro dos limites admissíveis com condutores padrão.	Severa elevação do Gradiente de Potencial de Superficie (GPR / EPS) e cristas de tensão extremamente pronunciadas sobre a malha.	No Cenário B, supera-se frequentemente o limite de fibrilação ventricular segundo a IEEE 80 / IEC 60479.
Efetividade das Hastes de Aterramento	Altamente efetivas se penetrarem a camada profunda ($l_v > h$), reduzindo drasticamente $R_g$.	Pouco efetivas; as hastes quase não dissipam corrente adicional no estrato rochoso profundo de alta resistividade.	No Cenário B, exige-se densificar a grade horizontal e expandir a superfície perimetral.
Fator de Redução da Camada de Brita ($C_s$)	A espessura superficial de pedra britada ($h_s$) proporciona uma mitigação padrão ($C_s \approx 0.7 - 0.85$).	A interação da baixa $\rho_1$ com a rocha superficial e $\rho_2$ infinita altera a reflexão do leito de brita.	Requer recalcular $C_s$ por meio de séries infinitas de dupla reflexão (brita–camada 1–camada 2).

Cálculo Rigoroso do Fator de Reflexão da Camada Superficial ($C_s$)

O fator de correção $C_s$ para a camada de material de alta resistividade superficial (brita/pedra britada $\rho_s \approx 3000\,\Omega\cdot m$) de espessura $h_s$ instalada sobre um solo bicamada não pode ser computado com a clássica fórmula empírica de Sunde para solo homogêneo. A formulação generalizada mediante expansão em série de coeficientes combinados resulta em:

onde $k_s = (\rho_1 - \rho_s)/(\rho_1 + \rho_s)$ e $r_d = 0.08 m$ representa o raio equivalente do pé humano modelado como um disco circular metálico de raio padronizado.

Análise Forense de Falhas Associadas à Caracterização Geoelétrica Deficiente

Os erros metodológicos na campanha de prospecção de resistividade ou a simplificação para solo homogêneo quando subjaz um estrato altamente reflexivo ($k \to +1$) resultaram historicamente em falhas catastróficas em sistemas de potência. Apresentam-se a seguir mecanismos de falha e incidentes forenses documentados:

Disrupção Dielétrica em Transformadores de Medição e Serviços Auxiliares por Elevação Crítica de GPR

Em uma subestação de 230 kV construída sobre uma formação de arenito superficial ($\rho_1 = 180\,\Omega\cdot m$, $h = 1.8 m$) sobre um embasamento granítico maciço ($\rho_2 = 4200\,\Omega\cdot m$, $k = +0.918$), o projeto original assumiu um solo homogêneo médio de $350\,\Omega\cdot m$. Durante um curto-circuito monofásico à terra com corrente simétrica $I_f = 22 kA$, a resistência real da malha foi de $2.85\,\Omega$ (frente aos $0.68\,\Omega$ projetados sob a hipótese homogênea).

A densa corrente que não pôde penetrar no granito profundo retornou lateralmente pelas canaletas de cabos de controle, elevando o potencial local de terra dos armários externos acima da rigidez dielétrica dos enrolamentos secundários dos transformadores de corrente (TC) e transformadores de potencial (TP). A diferença de potencial de modo comum perfurou o isolamento galvânico de 10 kV dos relés de proteção IED microprocessados na sala de comando, queimando os módulos de entradas e saídas analógicas e inibindo a atuação da proteção de retardo temporizado (ANSI 50/51N), o que estendeu a extinção da falta de 80 ms para 1.2 segundo e causou a explosão destrutiva do transformador de potência principal.

Degradação Térmica e Fusão de Blindagens de Cabos Subterrâneos de Alta Tensão

Em linhas subterrâneas de cabos isolados em XLPE de 115 kV com blindagens metálicas aterradas em ambas as extremidades (solid bonding), a premissa errónea de uma baixa resistividade em estratos profundos distorce a impedância mútua de laço terra-blindagem $\underline{Z}_{m(c-g)}$ segundo as fórmulas de Carson-Clem:

Ao encontrar-se um estrato profundo de alta resistividade ($\rho_2 \gg \rho_1$), a profundidade de retorno equivalente $D_e$ aumentou em mais de 300%. Isso incrementou drasticamente a reatância indutiva do retorno pelo solo, forçando a corrente de desbalanceamento e de falta a circular quase na sua totalidade através das blindagens metálicas dos cabos de potência. A densidade de corrente térmica $Jscreen$ superou o limite admissível do cobre segundo a IEC 60287 e IEC 60949:

provocando a delaminação térmica da cobertura externa semicondutora, o amolecimento do polietileno reticulado (XLPE) por sobretemperatura superior a 250 °C e a consequente falha dielétrica fase-terra por perfuração de arborescências de água (water trees) aceleradas eletrotermicamente.

Fluxo Metodológico de Simulação e Projeto Avançado no Vexten Suite

O ecossistema de engenharia integral Vexten Suite incorpora um motor determinístico e multifísico de inversão geoelétrica acoplado a módulos de curto-circuito (IEC 60909 / IEEE 141), dimensionamento de condutores e cabos (IEC 60287 / NBR/NEC) e otimização de malhas de aterramento (IEEE 80 / IEC 60479). A seguir, detalha-se a integração algorítmica e o protocolo metodológico aplicado em projetos de infraestrutura crítica:

Exemplo Prático de Simulação e Inversão Numérica no Vexten Grounding Engine

Realiza-se uma campanha de prospecção geoelétrica tetrapolar pelo método de Wenner no local de instalação de uma nova subestação industrial de 138/13.8 kV. A aquisição de dados forneceu os seguintes valores de campo:

Espaçamento $a$ [m]	Resistência Medida $R_m$ [$\Omega$]	Resistividade Aparente Experimental $\rho_a^{med}$ [$\Omega\cdot m$]	Resistividade Aparente Calculada $\rho_a^{calc}$ [$\Omega\cdot m$]	Erro 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

Ao processar a matriz de dados no módulo de otimização inversa não linear do Vexten Grounding, o algoritmo de Levenberg-Marquardt converge em 7 iterações com um erro global RMS de 1.32%, determinando o seguinte vetor de estado estratificado:

Análise de Engenharia e Otimização do Projeto do Sistema de Aterramento (SAT)

Dado que o coeficiente de reflexão $k$ é fortemente negativo ($k = -0.6812$), o estrato profundo apresenta uma resistividade mais de cinco vezes inferior à da camada superficial ($\rho_2 \ll \rho_1$). A otimização computacional no Vexten Grounding explora este fenômeno físico por meio das seguintes diretrizes executivas:

1. Instalação Estratégica de Hastes Profundas: Especifica-se a cravação de hastes de aço cobreado (copperweld) com comprimento mínimo $l_v = 6.0 m$ nos quatro cantos e nós perimétricos críticos. Ao ultrapassar a profundidade da interface ($h = 4.35 m$), os eletrodos verticais penetram $1.65 m$ no estrato de baixa resistividade ($\rho_2 = 82.50\,\Omega\cdot m$). Isso reduz a resistência total do sistema $R_g$ em 42.7% em comparação com uma malha puramente horizontal de igual comprimento total de condutor.
2. Desgestionamento do Condutor Horizontal: Graças à alta admitância da camada inferior, o espaçamento entre condutores horizontais da malha pode ser ampliado de um reticulado denso de $5 m \times 5 m$ para $10 m \times 10 m$, alcançando uma economia significativa de cobre sem exceder as tensões de passo e de toque limite estipuladas pela norma IEEE Std 80.
3. Integração com o Módulo de Curto-Circuito do Vexten Suite (IEC 60909): A impedância real do terreno estratificado é exportada automaticamente para o motor de curto-circuito do Vexten Suite, permitindo recalcular de maneira acoplada os fatores de divisão de corrente por cabos para-raios e blindagens de fase ($S_f$), garantindo que a corrente que efetivamente retorna pelo solo ($I_g = S_f \cdot I_k''$) reflita com exatidão a física multifilar do sistema de transmissão e distribuição circundante.

A adoção deste fluxo metodológico integral garante a máxima rigorosidade analítica, segurança dielétrica do pessoal e dos ativos eletrotérmicos, bem como uma otimização econômica e operacional na construção de sistemas de aterramento para subestações e usinas de geração.