2.5-D Electrical Resistivity Forward Modelling with Undulating Topography using a Modified Half-Space Analytical Solution

Título: Modelagem Direta 2.5-D de Resistividade Elétrica com Topografia Ondulada Usando uma Solução Analítica Modificada de Meio Semi-Espaço

1. Problema e Contexto

A imagem de resistividade de corrente contínua (DC) é uma ferramenta geofísica essencial para investigações de subsuperfície em diversas áreas (hidrogeologia, geotecnia, arqueologia). No entanto, os levantamentos de campo frequentemente ocorrem em terrenos com topografia ondulada.

• O Desafio: A representação precisa das variações topográficas na modelagem direta é crucial para inversões confiáveis. Métodos numéricos existentes, como o Método dos Elementos Finitos (MEF), enfrentam dificuldades devido à singularidade do potencial elétrico nos eletrodos de fonte.
• Limitação dos Métodos Atuais: A técnica padrão de "remoção de singularidade" divide o potencial total em uma parte primária (analítica, para um meio homogêneo plano) e uma parte secundária (numérica). Para topografias planas, a solução analítica primária é exata. Para terrenos ondulados, métodos anteriores (como o de Penz et al., 2013) tentam usar a solução de meio plano, assumindo que a topografia local próxima à fonte é suave.
• A Falha Crítica: Em terrenos naturais com cantos agudos, curvaturas altas ou descontinuidades de inclinação (comuns em discretizações lineares de elementos finitos), a suposição de superfície plana falha. Isso gera um erro geométrico fundamental: o ângulo sólido subtendido pela topografia real na fonte difere do ângulo sólido de $2\pi$ assumido pela solução de meio plano. Isso resulta em erros significativos de modelagem, mesmo com refinamento de malha, e exige o alisamento artificial do terreno ou malhas extremamente finas, o que é computacionalmente custoso.

2. Metodologia

Os autores propõem uma nova estratégia de remoção de singularidade que deriva uma nova solução analítica primária baseada no ângulo sólido real subtendido pela topografia local na posição da fonte.

• Formulação Matemática:
  • Em vez de assumir um meio semi-espaço plano ($\beta = \pi$, ângulo sólido $S = 2\pi$), os autores derivam a solução para um potencial gerado por uma fonte pontual no vértice de uma cunha em V (um wedge).
  • A solução analítica para o potencial primário ($u_p$) é modificada para incluir o ângulo sólido $S$ da região exterior:
    $$u_p(r) = \frac{I}{\sigma_0 S r}$$
    onde $S = 4\pi - 2\beta$ (sendo $\beta$ o ângulo da cunha).
  • No domínio 2.5-D (após transformada de Fourier), a solução é expressa usando a função de Bessel modificada $K_0$, escalada pelo fator do ângulo sólido.
• Implementação no MEF:
  • O potencial primário modificado é usado para calcular o termo de fonte na equação do potencial secundário.
  • A condição de contorno de Neumann (fluxo zero) na superfície é satisfeita implicitamente pela nova solução analítica, que se adapta à geometria da cunha (ou à geometria linear aproximada da malha), eliminando a necessidade de suavizar a topografia artificialmente.
  • O método funciona tanto com elementos lineares (Q4) quanto de ordem superior, sem exigir refinamento excessivo da malha próximo à fonte.

3. Principais Contribuições

1. Derivação Analítica de Cunha V: Desenvolvimento de uma expressão analítica fechada para o potencial primário em uma geometria de cunha, que serve como base para lidar com descontinuidades de inclinação e curvaturas altas.
2. Correção Geométrica da Singularidade: Substituição da suposição de ângulo sólido constante ($2\pi$) pelo cálculo do ângulo sólido real subtendido pela malha de elementos finitos na fonte. Isso alinha a solução analítica com a geometria discretizada real.
3. Independência de Suavização: O método elimina a necessidade de modificar manualmente a topografia (alisamento) ou de usar malhas extremamente densas e elementos de alta ordem apenas para capturar a singularidade, mantendo a eficiência computacional.
4. Robustez com Elementos Lineares: Demonstra que elementos lineares (Q4), quando combinados com esta nova formulação, podem atingir precisão superior à de elementos de ordem superior com a formulação antiga, especialmente em terrenos complexos.

4. Resultados e Validação

Os autores validaram o método através de vários cenários sintéticos:

• Casos de Verificação (Plano e 2 Camadas): O método reproduziu soluções analíticas conhecidas com erros inferiores a 0,1%.
• Modelo de Trincheira em V (V-Trench):
  • Em uma trincheira com inclinação de 15°, o método tradicional (meio plano) apresentou erros de aproximadamente 16,5%, independentemente do refinamento da malha.
  • O método proposto manteve erros abaixo de 0,1% mesmo com malhas lineares grosseiras (elementos de 1 m).
  • A análise teórica mostrou que o erro do método antigo depende exclusivamente da razão entre o ângulo sólido real e $2\pi$, explicando a constância do erro.
• Modelo de Colina/Vale Senoidal:
  • Para topografias suaves mas com curvatura, o método tradicional falhou em malhas lineares (erros ~2%), enquanto o método proposto manteve erros < 0,1%.
  • Em pontos de transição abrupta (junção vale-plano), o método tradicional gerou erros de até 10%, enquanto o método proposto corrigiu isso calculando o ângulo sólido específico para aquela geometria.
• Eficiência: O método permite o uso de malhas mais grosseiras e elementos lineares, reduzindo drasticamente o custo computacional em comparação com métodos que exigem malhas refinadas ou elementos curvos de alta ordem para obter a mesma precisão.

5. Significado e Conclusão

Este trabalho oferece uma alternativa robusta e computacionalmente eficiente para a modelagem direta de resistividade elétrica em terrenos complexos. Ao integrar corretamente a geometria local (ângulo sólido) na solução analítica primária, o método resolve o problema fundamental de incompatibilidade geométrica que afeta as técnicas atuais.

• Impacto Prático: Facilita a implementação em frameworks de inversão 2.5-D em larga escala, onde a velocidade de cálculo é crítica.
• Conclusão: A abordagem baseada no ângulo sólido é superior para lidar com cantos agudos e descontinuidades, garantindo alta precisão (< 0,1%) sem os custos computacionais proibitivos de métodos anteriores. Os autores sugerem que esta abordagem pode ser estendida futuramente para modelagem 3-D completa.