A Nonmatching-Grid Discontinuous Galerkin Method for Nonlinear Reactive Transport in Hybrid-Dimensional Fractured Media

No mundo oculto sob nossos pés, a água e os produtos químicos dissolvidos movem-se através de um labirinto complexo de rocha. Esta rocha raramente é uniforme; é frequentemente um bloco sólido, chamado de matriz, atravessado por finas fendas conhecidas como fraturas. Estas fraturas são as rodovias do subsolo. Elas podem atuar como canais super-rápidos que aceleram fluidos através da rocha, ou como barreiras que bloqueiam o movimento inteiramente, dependendo de sua natureza. Como estas rachaduras são frequentemente microscópicas em largura, mas podem estender-se por quilômetros em comprimento, modelar como as substâncias viajam através delas é um desafio massivo para os cientistas. Se tentarem construir um modelo computacional que capture cada detalhe minúsculo da abertura da rachadura, o cálculo torna-se tão pesado que nem mesmo os supercomputadores mais poderosos conseguem resolvê-lo. Para contornar isso, os pesquisadores usam um atalho inteligente chamado modelo híbrido-dimensional. Em vez de tratar a rachadura como uma fina fatia tridimensional, eles a achatam em uma linha bidimensional ou uma linha unidimensional que reside dentro da rocha. Isso permite que simulem o fluxo sem a necessidade de resolver a espessura minúscula da rachadura em si.

A dificuldade surge quando os cientistas tentam conectar a matemática da rocha sólida com a matemática destas rachaduras achatadas. No mundo real, a rocha e a rachadura estão em constante conversa, trocando fluidos e produtos químicos através de sua fronteira compartilhada. Para simular isso, os computadores geralmente dividem a rocha e a rachadura em pequenos pedaços, como um mosaico, para calcular o fluxo. Durante anos, a abordagem padrão exigia que os pedaços da rocha e os pedaços da rachadura se alinhassem perfeitamente ao longo de sua borda compartilhada, como duas peças de um quebra-cabeça se encaixando. Se a rocha fosse dividida em grandes quadrados e a rachadura fosse dividida em segmentos minúsculos, o computador tinha que forçá-los a coincidir, muitas vezes criando uma camada complexa e artificial de dados para fazer a ponte. Este requisito tornava os modelos rígidos e difíceis de usar quando diferentes partes de um sistema precisavam ser refinadas em diferentes escalas.

Um novo estudo de Alvaro H. Salas e Carlos D. Acosta, da Universidad Nacional de Colombia, aborda essa rigidez ao fazer uma pergunta simples, mas profunda: e se a rocha e a rachadura não precisarem se alinhar de forma alguma? Eles desenvolveram um método onde a rocha e a fratura podem ser divididas em pedaços de tamanhos e formas completamente diferentes, gerados independentemente um do outro, e ainda assim trocar informações perfeitamente. Imagine tentar despejar água de um balde grande e grosseiro em uma peneira fina e detalhada. Normalmente, você precisaria de um funil especial para fazer a água ser transferida suavemente sem derramar. Os pesquisadores encontraram uma maneira de despejar a água diretamente, calculando a quantidade exata que passa por cada ponto onde a borda do balde toca a peneira, independentemente de quão desajustados estejam os furos.

A equipe criou uma estrutura matemática que trata a troca de produtos químicos entre a rocha e a fratura como uma interação física direta. Em vez de forçar as duas grades a coincidir ou usar um intermediário para traduzir os dados, eles construíram um sistema que calcula a transferência nas pequenas seções sobrepostas onde a borda da rocha e a superfície da fratura realmente se encontram. Esta abordagem é chamada de método de grade não correspondente (nonmatching-grid method). Em suas simulações, eles testaram este sistema com reações químicas complexas, onde as substâncias mudam conforme se movem, um cenário notoriamente difícil de modelar. Eles provaram matematicamente que esta conexão direta é estável e precisa, mesmo quando a rocha é dividida em blocos muito grandes e a fratura é dividida em segmentos muito finos, ou vice-versa. A estabilidade de seu método não depende da razão entre os tamanhos desses blocos, o que significa que o computador não trava ou produz erros selvagens apenas porque uma grade é muito mais fina que a outra.

Uma das descobertas mais significativas de seu trabalho é que este método direto é superior às abordagens indiretas mais antigas que tentam projetar os dados de uma grade sobre a outra. Quando a troca de produtos químicos é forte, os métodos de projeção antigos perdem informações. Eles essencialmente suavizam os detalhes da fratura, descartando as variações de pequena escala que são cruciais para uma previsão precisa. Os pesquisadores demonstraram que seu método direto retém todas as informações necessárias. Em testes onde a transferência de produtos químicos era intensa, o método antigo produziu erros vinte e um por cento maiores do que o novo método. Além disso, quando refinaram a grade da fratura para ser ainda mais detalhada enquanto mantinham a grade da rocha fixa, o método antigo parou de melhorar, estagnando-se. O novo método, no entanto, continuou a tornar-se mais preciso, provando que captura a verdadeira física da interação sem limitações artificiais.

Os pesquisadores também mostraram que este método funciona eficientemente em computadores. Eles escreveram um programa de software para testar sua teoria, executando simulações em um computador desktop padrão. Eles compararam o tempo necessário para configurar os cálculos para grades correspondentes versus grades não correspondentes. Os resultados mostraram que o trabalho extra necessário para lidar com as peças desajustadas foi insignificante. O tempo para montar as equações foi quase idêntico, o que significa que os cientistas agora podem usar grades muito mais flexíveis e realistas sem pagar um custo em tempo de computação. Esta flexibilidade é vital para aplicações do mundo real, onde diferentes partes de uma formação geológica podem ser modeladas por diferentes equipes usando softwares diferentes, ou onde a rede de fraturas é tão complexa que uma única grade perfeitamente alinhada é impossível de criar.

O estudo também explorou como o método se comporta sob condições extremas. Eles testaram cenários onde as reações químicas eram muito fortes e a transferência de fluidos era rápida. Nesses casos difíceis, o método permaneceu robusto, convergindo para uma solução sem falhar. Eles até testaram a geometria das fraturas, simulando rachaduras que dobram ou se interceptam umas com as outras. Embora os testes principais tenham usado uma rachadura reta, a lógica subjacente de seu método sustentou-se quando aplicada a essas formas mais complexas, sugerindo que a abordagem pode ser estendida para as redes emaranhadas encontradas na natureza. Os pesquisadores foram cuidadosos ao notar que seu trabalho atual foca em tipos específicos de rochas e geometrias de fraturas, e que estender o método para redes tridimensionais completas, curvas e que se interceptam é uma tarefa para o futuro.

O que torna este trabalho particularmente importante é que ele remove uma barreira de longa data na simulação de meios fraturados. Por décadas, a necessidade de grades correspondentes limitou como os cientistas podiam modelar esses sistemas, muitas vezes forçando-os a comprometer o detalhamento da fratura ou da rocha para fazer a matemática funcionar. Ao provar que os dois podem ser acoplados diretamente e de forma conservativa sem uma correspondência perfeita, os pesquisadores abriram as portas para simulações mais realistas e flexíveis. Isso significa que, no futuro, modelos de contaminação de águas subterrâneas, armazenamento de carbono ou extração de energia geotérmica podem ser construídos com maior liberdade, permitindo que a rocha e as rachaduras sejam representadas exatamente como são, em vez de como devem ser para se ajustarem a uma restrição matemática rígida. O método fornece uma base nova e estável para entender como fluidos e produtos químicos se movem através da terra fraturada, garantindo que a conexão entre a rocha sólida e as rachaduras seja tão fluida no computador quanto é no solo.