Lattice-Boltzmann para Meios Porosos: 100M+ Horas em GPU

Por que os espaços minúsculos nas rochas importam

Quando armazenamos energia no subsolo, capturamos dióxido de carbono ou operamos células a combustível, a forma como dois fluidos distintos serpenteiam por espaços minúsculos na rocha pode determinar o sucesso da tecnologia. No entanto, observar esses fluidos se moverem dentro de rochas reais é extremamente lento e caro. Este artigo apresenta um enorme conjunto de dados abertos de simulações computacionais que fazem o segundo melhor: recriam como água e fluidos semelhantes a óleo se deslocam um pelo outro em amostras de rocha realistas, usando mais de 100 milhões de horas de processamento em unidades gráficas. O resultado é um recurso compartilhado que qualquer pesquisador pode usar para testar ideias sobre escoamento subterrâneo sem precisar de um acelerador de partículas ou de seu próprio supercomputador.

Explorando o movimento dos fluidos em rochas digitais

Os autores concentram-se em “escoamento bifásico”, em que dois fluidos que não se misturam, como água e óleo, se movem juntos por um labirinto de poros na rocha. Uma grandeza-chave para engenheiros é a permeabilidade relativa, que indica com que facilidade cada fluido se desloca quando o outro está presente. Normalmente, medições completas exigem semanas de trabalho laboratorial cuidadoso para cada amostra de rocha e conjunto de condições. Em vez disso, a equipe usou um pacote de simulação especializado, o LBPM, para calcular o escoamento diretamente em imagens 3D de rochas reais. Essas rochas digitais vieram de exames de microtomografia por raios X de arenito e vidro sinterizado, capturando formas e tamanhos realistas dos poros até algumas micrômetros.

Um experimento virtual massivo

Executar simulações realistas nessa escala fina continua sendo muito dispendioso. A equipe aproveitou computadores de alto desempenho com milhares de GPUs para varrer condições que seriam impraticáveis em laboratório. Eles variaram o quanto a superfície da rocha prefere um fluido em relação ao outro (seu comportamento de molhabilidade) e com que intensidade os fluidos são empurrados (capturado por um parâmetro chamado número capilar). Para quatro materiais porosos diferentes, realizaram protocolos de escoamento tanto em “estado estacionário” quanto em “estado não estacionário” que imitam experimentos padrão de enchimento de núcleo usados na indústria de petróleo e gás. No total, produziriam 50 curvas completas de permeabilidade relativa e mais de 25.000 configurações fluidas distintas.

Vendo formas, não apenas médias

Além das taxas médias de escoamento, as simulações acompanham as formas detalhadas e as conexões dos fluidos ao longo do tempo — informações quase impossíveis de obter em escala a partir de experimentos. O LBPM separa caminhos de fluxo conectados de manchas isoladas, às vezes chamadas de gânglios, e mede seu volume, área de superfície, curvatura, conectividade, pressão e movimento. Essas grandezas são registradas a cada passo da simulação em tabelas de texto simples, de modo que os usuários podem reconstruir como bolsões aprisionados se formam, destacam-se, reconectam-se ou drenam lentamente. Ao comparar diferentes tipos de rocha e padrões de molhabilidade, o conjunto de dados revela como mudanças sutis na preferência de superfície podem deslocar onde o fluido fica preso e quão facilmente ele se move, ajudando a explicar tendências observadas em medições laboratoriais.

Lições sobre rocha molhada e fluido aprisionado

Usando o conjunto de dados, os autores validam que suas simulações se comportam de forma sensata verificando padrões conhecidos. Por exemplo, à medida que a rocha se torna mais oleofílica, a quantidade de óleo remanescente após uma lavagem com água tende a diminuir, coincidindo com experimentos anteriores. Em casos mais hidrofílicos com garganta estreita mal resolvida, as simulações mostram rompimento precoce do fluido não-molhante e curvas de escoamento incomumente planas, ilustrando como a resolução da imagem pode levar a interpretações erradas. Em outras rochas com poros melhor resolvidos, a mudança nas curvas de escoamento com variação da molhabilidade alinha-se com estudos laboratoriais anteriores. Medidas de conectividade confirmam que, sob certas condições intermediárias, ambos os fluidos permanecem altamente interconectados, um estado associado a formas de interface complexas e aglomerados aprisionados de longa duração.

Uma base compartilhada para trabalhos futuros

Em termos simples, este artigo entrega um mapa detalhado de como dois fluidos podem compartilhar os mesmos espaços minúsculos dentro de uma rocha, sob uma grande variedade de condições, tudo codificado em dados reutilizáveis. As simulações são cruzadas com experimentos avançados de imagem por raios X e organizadas de modo que outros possam facilmente calcular novas curvas resumo, treinar modelos de aprendizado de máquina ou testar novas teorias sobre como detalhes em escala de poro se somam ao escoamento em grande escala. Para quem trabalha com armazenamento subterrâneo de carbono, armazenamento de hidrogênio, remediação de águas subterrâneas ou componentes de células a combustível, este conjunto de dados aberto oferece um atalho poderoso: em vez de começar do zero, é possível construir diretamente sobre mais de 100 milhões de horas de GPU de experimentos digitais cuidadosamente curados.

Citação: Armstrong, R.T., Tavakkoli, O., Da Wang, Y. et al. Lattice-Boltzmann for Porous Media: 100M⁺ GPU Hours. Sci Data 13, 697 (2026). https://doi.org/10.1038/s41597-026-06823-1

Palavras-chave: meios porosos, escoamento bifásico, rocha digital, lattice boltzmann, permeabilidade relativa