Curvelet-Regularized SPDE Inversion on Piecewise-Planar Fractures with Trace-Graph Coupling

Resumo Técnico: Inversão de EDP com Regularização por Curvelets em Fraturas Planares por Partes com Acoplamento de Grafo-Traço

Enunciado do Problema
Métodos geoestatísticos clássicos, como krigagem e modelos de Campo Aleatório Gaussiano (GRF), tipicamente assumem que a correlação é uma função da distância euclidiana. Essa premissa falha em meios fraturados, onde o transporte e o fluxo se propagam preferencialmente ao longo de redes de fraturas conectadas, em vez de através da rocha intacta. Consequentemente, pontos próximos no espaço euclidiano podem ter um acoplamento fraco, enquanto pontos distantes conectados por um corredor de fratura podem estar fortemente acoplados. Além disso, campos controlados por fraturas frequentemente exibem características anisotrópicas, curvilíneas e multiescala (ex: cristas, frentes, canais) que são mal representadas apenas por priors de suavidade isotrópica. O artigo aborda o desafio de reconstruir um campo escalar latente (ex: log-transmitividade, teores geoquímicos) em superfícies de fratura, dados conjuntos de medições pontuais esparsas e polilinhas de traços de fraturas 3D.

Metodologia
O framework proposto opera em duas etapas principais: inferência geométrica e reconstrução esparsa-para-densa.

1. Inferência Geométrica (Traços 3D para Mapas de Planos):

   • Ajuste de Planos: Polilinhas de traços 3D são ajustadas a planos usando PCA/SVD. Traços degenerados (quase colineares) são sinalizados e excluídos ou atribuídos post-hoc.
   • Agrupamento (Clustering): Hipóteses de planos são agrupadas usando uma métrica combinada de desvio de vetor normal e erro de offset (estilo DBSCAN) para identificar planos de fratura dominantes.
   • Coordenadas Locais: Para cada plano dominante, um sistema de base ortonormal local $(u, v)$ é estabelecido. Amostras 3D são projetadas nesses mapas locais e um grid 2D é construído para discretizar o campo latente $z_p(u, v)$.

2. Modelos de Reconstrução:

   • Linha de Base (Gaussiano/EDP): Um estimador MAP Gaussiano padrão é formulado usando uma Matriz de Precisão Esparsa derivada de uma Equação Diferencial Parcial Estocástica (EDP), onde o operador de precisão é definido como $Q = (\kappa^2 I + L_{uv})^2$, sendo $L_{uv}$ o Laplaciano 2D discreto. Isso fornece uma linha de base suave e isotrópica.
   • Regularização por Curvelets: Para capturar estruturas direcionais alinhadas às fraturas, o modelo adiciona uma penalidade $\ell_1$ sobre os coeficientes de curvelets discretos não decimados (UDCT). A função objetivo minimiza o erro de dados, o prior quadrático da EDP e o termo de esparsidade de curvelets:
     $$\min_z \frac{1}{2}\|W^{1/2}(Hz - y)\|_2^2 + \frac{1}{2}z^\top Q z + \lambda \|Cz\|_1$$
     onde $C$ é a transformada de curvelet.
   • Acoplamento Grafo-Traço (Extensão Versão 2): Para impor a distância "ao longo da fratura", um Laplaciano de Grafo $L_G$ é construído na rede de traços. Um termo de acoplamento $z^\top P^\top L_G P z$ é adicionado, onde $P$ mapeia valores de grid para nós do grafo. Isso penaliza as diferenças entre vértices conectados ao longo da rede de fraturas, independentemente da distância euclidiana.

3. Otimização:

   • O termo de curvelet não suave é resolvido usando o Método Alternado de Multiplicadores Dirigidos (ADMM).
   • O algoritmo divide o problema em uma resolução linear esparsa (bloco quadrático da EDP) e um passo de contração (shrinkage) por coeficiente (bloco de curvelet).
   • A convergência é garantida sob suposições padrão de convexidade, com critérios de parada baseados em normas de resíduos primal e dual.

Principais Contribuições

• Framework para Meios Fraturados: Um pipeline reproduzível que infere suportes de fraturas planas por partes a partir de traços 3D e reconstrói campos de propriedades a partir de dados esparsos.
• Prior de Esparsidade Direcional: A integração formal de transformadas de curvelet em modelos geoestatísticos baseados em EDP para visar especificamente singularidades curvilíneas e características anisotrópicas inerentes às fraturas.
• Regularização Impulsionada pela Conectividade: A definição de uma métrica de distância "ao longo da fratura" via Laplacianos de grafo, permitindo que o modelo imponha suavidade ao longo das redes de fraturas em vez de no espaço euclidiano.
• Arquitetura Modular: O artigo distingue entre uma "Versão 1" (inferência por plano com priors Gaussianos e de curvelet) e uma extensão "Versão 2" (acoplamento de grafo explícito e restrições de interseção), fornecendo um roteiro claro para aumentar a complexidade física.
• Benchmarks Reproduzíveis: Os autores fornecem benchmarks sintéticos (S1: anisotropia curvilínea; S2: conectividade impulsionada por grafo) e uma série completa de estudos de ablação para isolar as contribuições da suavidade euclidiana, esparsidade direcional e priors de conectividade.

Resultos
Experimentos numéricos em dados sintéticos demonstram:

• Esparsidade Direcional (S1): O modelo regularizado por curvelets supera significativamente a krigagem ordinária e a linha de base pura de EDP na reconstrução de anomalias do tipo crista alinhadas aos traços de fratura, alcançando menor RMSE e maior $R^2$ em amostras de teste.
• Conectividade (S2): Em cenários onde a verdade de campo é impulsionada pela conectividade do grafo em vez da distância euclidiana, a adição do prior de grafo (Versão 2) melhora a precisão da reconstrução. O modelo completo (EDP + Curvelets + Grafo) produz o melhor desempenho, capturando efetivamente a estrutura de correlação não euclidiana.
• Eficiência: O solver ADMM converge eficientemente, aproveitando álgebra linear esparsa para o bloco quadrático e transformadas rápidas ($O(M \log M)$) para o bloco de curvelet.

Significância e Escopo
O artigo posiciona-se como um framework de métodos, e não como um solver físico completo. Sua principal significância reside em fornecer uma camada de interpolação/inversão estatisticamente controlada que pode ser acoplada a modelos de fluxo ou transporte subsequentes.

• Alegações Modestas: Os autores declaram explicitamente que não reivindicam "completude física". O framework é projetado para reconstruir campos condicionados à geometria de fraturas, não para resolver as EDPs governantes do fluxo de fluidos diretamente.
• Fundação para Trabalhos Futuros: A extensão "Versão 2" e o roteiro (Seção 14) delineiam a futura integração com física não linear e inversões restritas por EDP.
• Validação: O artigo enfatiza um protocolo de validação rigoroso, incluindo estudos de sensibilidade sobre hiperparâmetros ($\kappa, \lambda, \rho$), resolução de grid e densidade de amostragem, estabeleando uma base para geoestatística consciente de redes em meios fraturados.