An Upward-Continuation-Based Framework for Subsurface Structural Reconstruction from Gravity Anomalies

Resumo Técnico: Um Framework Baseado em Continuação para Cima para Reconstrução Estrutural de Subsuperfície a partir de Anomalias de Gravidade

Definição do Problema
A interpretação de dados de gravidade é fundamental para investigar a heterogeneidade da densidade de subsuperfície, mas é inerentemente não única; distribuições de densidade substancialmente diferentes podem gerar respostas anômalas semelhantes. Embora técnicas de inversão e modelagem direta existam, elas frequentemente dependem fortemente de restrições geológicas prévias, parametrização complexa ou otimização computacionalmente intensiva. Métodos baseados em transformações (ex: continuação para cima, sinal analítico, gradiente total normalizado) foram desenvolvidos para extrair informações estruturais sem inversão, mas a maioria foca no realce de anomalias, detecção de bordas ou separação de comprimento de onda, em vez de reconstruir a geometria espacial completa dos corpos causativos preservando a polaridade da anomalia. Métodos baseados em gradiente, como o Gradiente Total Normalizado (NFG), são eficazes para detecção de bordas e estimativa de profundidade, mas não representam diretamente a geometria completa ou a polaridade de densidade da fonte. Há uma necessidade de um framework independente de modelo que possa reconstruir a geometria estrutural de subsuperfície diretamente de anomalias de gravidade, fornecendo tanto a caracterização geométrica quanto restrições de profundidade sem exigir dados geológicos externos.

Metodologia
Os autores propõem o framework de Restrição de Decomposição por Continuação para Cima (UCRD), integrado ao método do Gradiente Total Normalizado (NFG), para abordar essas limitações. O fluxo de trabalho consiste em três estágios primários:

1. Continuação para Cima de Múltiplos Níveis e Decomposição Espectral:
   O campo de anomalia de gravidade de Bouguer é decomposto em níveis de comprimento de onda restrito (RWL) usando continuação para cima de múltiplos níveis. O processo começa com o campo observado ($U_0$) e aplica sucessivas continuações para cima ($U_n$) com alturas exponencialmente crescentes ($\Delta z_n$) para isolar fontes progressivamente mais profundas. O nível de comprimento de onda restrito para cada etapa é calculado como a diferença entre o campo na altura anterior e a altura atual ($RW L_n = U_{n-1} - U_n$). Isso separa o campo gravitacional em componentes associados a intervalos de profundidade específicos.

2. Transformação de Sinal Analítico e Normalização de Escala:
   Para aumentar a continuidade estrutural e preservar a interpretabilidade física, cada componente RWL passa por uma transformação de sinal analítico tridimensional. Um procedimento de normalização adaptável à escala é então aplicado para restaurar as unidades físicas (mGal) e preservar a polaridade da anomalia. Isso envolve o cálculo de um fator de escala específico para o nível baseado no alcance dinâmico do RWL e da amplitude do sinal analítico. O resultado é uma representação de Comprimento de Onda Restrito Escalado (SRWL). Os mapas SRWL capturam a magnitude e a direção dos contrastes laterais de densidade em cada nível de profundidade, fornecendo uma representação pseudo-3D da evolução da anomalia.

3. Integração com NFG para Restrições de Profundidade:
   Visto que a continuação para cima não fornece informações de profundidade física direta, a reconstrução geométrica do UCRD é complementada pelo método NFG. O NFG é aplicado usando uma formulação de série de Fourier com extensão de espelhamento de sinal para evitar descontinuidades de borda. Um procedimento adaptativo determina o número ideal de harmônicos ($N$) rastreando as respostas máximas de NFG através de intervalos harmônicos crescentes. Isso fornece restrições de profundidade quantitativas que são integradas com a geometria estrutural do UCRD.

Principais Contribuições

• Framework UCRD: Introdução de um método de decomposição que transforma anomalias de gravidade em imagens fisicamente interpretáveis da geometria estrutural de subsuperfície sem inversão.
• Representação SRWL: Desenvolvimento de um procedimento de escala que restaura as unidades físicas e a polaridade da anomalia aos resultados do sinal analítico, aumentando a visibilidade de características estruturais obscurecidas nos dados brutos.
• Fluxo de Trabalho Complementar: Uma abordagem unificada onde o UCRD delineia a geometria estrutural e a polaridade, enquanto o NFG fornece restrições de profundidade independentes, reduzindo a ambiguidade inerente às interpretações baseadas em gradiente.
• Independência de Modelo: A metodologia opera diretamente nos dados de gravidade, não exigindo informações geológicas prévias, dados sísmicos ou otimização computacionalmente intensiva.

Resultados
O framework foi validado usando três modelos sintéticos de complexidade crescente:

1. Corpos Isolados e Escalonados: O método delineou com sucesso seis corpos cúbicos e três estruturas escalonadas com profundidades e contrastes de densidade variados. O UCRD reconstruiu a distribuição espacial e a polaridade, enquanto o NFG identificou as profundidades correspondentes. Os resultados mostraram forte concordância com o modelo real, mesmo para corpos profundos parcialmente mascarados por anomalias rasas.
2. Fontes Sobrepostas: Em um modelo com prismas parcialmente sobrepostos, o UCRD separou com sucesso os componentes estruturais dominantes de fontes rasas e profundas, resolvendo as assinaturas de anomalia compostas. O NFG forneceu estimativas de profundidade estáveis para ambos os componentes, distinguindo geometrias sobrepostas que são difíceis de interpretar via análise de gravidade padrão.
3. Estruturas Lateralmente Extensas: Para um modelo composto por vinte e cinco blocos distribuídos em várias profundidades, o UCRD recuperou a geometria geral da estrutura de grande escala, revelando componentes tanto rasos quanto profundos. O NFG identificou níveis de profundidade distintos e limites estruturais dentro do sistema de fonte extensa.

Em todos os casos, as transformações SRWL produziram assinaturas de anomalia mais nítidas e localizadas em comparação com os componentes de comprimento de onda restrito padrão. O fluxo de trabalho combinado UCRD–NFG demonstrou a capacidade de recuperar padrões estruturais horizontais e verticais através de uma ampla gama de configurações de fonte.

Significância e Alegações
O artigo afirma que o framework UCRD–NFG proposto oferece uma ferramenta prática e computacionalmente eficiente para a interpretação de gravidade. Sua principal significância reside na capacidade de fornecer uma caracterização preliminar da geometria de subsuperfície e profundidade sem a necessidade de inversão ou restrições externas. Os autores posicionam o método não como um substituto para a modelagem direta ou inversão conjunta, mas como um passo preliminar valioso que pode gerar modelos iniciais, restrições geométricas e orientação de parametrização para investigações subsequentes mais intensivas. O framework é particularmente destacado para investigações de escala regional, onde a caracterização estrutural rápida é necessária. O estudo reconhece limitações inerentes à decomposição espectral, como a sobreposição parcial de comprimento de onda entre os níveis e efeitos de suavização da interpolação vertical, mas afirma que o método fornece uma abordagem robusta e intuitiva para o imageamento estrutural.