Revisiting the hydromechanical formulation of a micromechanics-based phase-field model for poro-elastoplastic media

Profundamente sob a superfície da Terra, nas rochas porosas que contêm petróleo, gás e calor geotérmico, uma batalha silenciosa é travada constantemente entre a pedra sólida e os fluidos presos em seus minúsculos poros. Quando engenheiros tentam extrair esses recursos, muitas vezes bombeiam fluidos para dentro da rocha para forçá-la a se separar, um processo conhecido como fraturamento hidráulico. Por décadas, cientistas modelaram esse processo tratando a rocha como uma esponja elástica simples que se estica e retorna à forma original. No entanto, as rochas reais são mais complexas; elas também podem se deformar permanentemente, esfarelar e deslizar ao longo de rachaduras internas, especialmente quando comprimidas por altas pressões. Compreender exatamente como essas rochas quebram quando tanto o esqueleto sólido quanto o fluido em seu interior estão empurrando e puxando é crucial para perfurar poços com segurança, gerenciar o armazenamento de resíduos e até mesmo prever eventos geológicos naturais, como o estalo de icebergs ou o magma empurrando através da crosta. O desafio reside em criar uma descrição matemática que capture a realidade desordenada da quebra da rocha sem se perder em cálculos impossíveis.

Uma equipe de pesquisadores revisitou as equações fundamentais usadas para descrever esse processo de quebra, descobrindo que um atalho comum usado em modelos anteriores leva a um erro significativo. No mundo da mecânica das rochas, os cientistas usam uma abordagem de "campo de fase" (phase-field), que trata uma rachadura não como uma linha nítida e invisível, mas como uma zona borrada onde o material está enfraquecendo gradualmente. Isso permite que computadores simulem como uma rachadura cresce e muda de forma. Para fazer isso com precisão, o modelo deve considerar dois comportamentos distintos: quando a rocha está sendo puxada (tensão) e quando está sendo espremida (compressão). Em tensão, pequenas rachaduras se abrem; em compressão, elas se fecham e as superfícies rugosas em seu interior deslizam umas contra as outras, criando fricção. Os pesquisadores descobriram que muitos modelos existentes, que são amplamente usados na indústria e na academia, confundem as variáveis usadas para descrever a energia do sistema. Especificamente, eles descobriram que tratar a pressão do fluido como uma entrada direta na equação de energia, em vez da quantidade de fluido nos poros, cria um "salto" ou descontinuidade matemática. Esse salto significa que o modelo prevê que a resistência da rocha muda abruptamente e de forma irrealista conforme ela alterna entre ser puxada e ser espremida, levando a previsões imprecisas de como e onde uma fratura se formará.

Para corrigir isso, os pesquisadores desenvolveram uma nova formulação consistente que rastreia cuidadosamente a relação entre a deformação da rocha e o conteúdo de fluido. Eles construíram seu modelo sobre uma visão microscópica da rocha, imaginando-a como uma coleção de minúsculas rachaduras em forma de moeda preenchidas com fluido. Ao analisar como essas micro-rachaduras se abrem e se fecham sob diferentes tensões, eles derivaram um conjunto de regras que garantem que a superfície de resistência da rocha — o limite matemático que define quando ela quebrará — permaneça suave e contínua, independentemente de como a pressão muda. Eles provaram que, ao usar um tipo específico de regra de fluxo que liga a deformação permanente da rocha diretamente às forças que atuam sobre ela, e ao incluir corretamente um termo de acoplamento ausente na equação de energia, o modelo pode descrever com precisão a transição da tensão para a compressão. Essa correção não é apenas um pequeno ajuste; ela altera fundamentalmente a forma como o modelo calcula a força motriz que impulsiona uma rachadura adiante.

A equipe testou seu novo modelo contra dois cenários muito diferentes para ver se ele se sustentava. Primeiro, simularam um benchmark clássico de fraturamento hidráulico, onde fluido é injetado em uma rocha para criar uma rachadura. Eles compararam seus resultados com uma solução analítica conhecida, que é uma resposta matemática precisa derivada para condições ideais. O novo modelo coincidiu com essa solução de forma muito mais próxima do que as formulações mistas antigas, prevendo corretamente a pressão necessária para iniciar a fratura e a largura da rachadura conforme ela crescia. Os modelos antigos tendiam a superestimar a pressão do fluido necessária, uma falha que poderia levar a operações de perfuração ineficientes ou inseguras. No segundo teste, simularam uma amostra de rocha sendo espremida de dois lados, uma condição que mimetiza o ambiente subterrâneo profundo, onde as rochas estão sob imensa pressão. Aqui, a rocha não apenas estalou; ela se deformou plasticamente, formando uma banda de cisalhamento onde o material deslizou sobre si mesmo antes de quebrar. O novo modelo reproduziu com sucesso esse comportamento complexo, capturando tanto as fraturas de cisalhamento causadas pelo esmagamento mecânico quanto as fraturas de tensão impulsionadas pela injeção de fluido.

Os resultados dessas simulações revelam que a deformação plástica da rocha atua como uma barreira significativa ao crescimento da fratura. Quando a rocha tem permissão para se deformar permanentemente, ela absorve mais energia, o que significa que a pressão do fluido dentro da rachadura deve ser maior para manter a fratura aberta e em crescimento. Essa resistência extra altera todo o perfil da fratura, tornando-a mais estreita e exigindo mais força para se propagar do que um modelo puramente elástico sugeriria. Além disso, nos testes de compressão, o modelo mostrou como a pressão de poro dentro da rocha muda conforme a fratura se forma. À medida que a rocha se deforma e cria espaço para o movimento do fluido, a pressão cai localmente, atraindo mais fluido para a zona danificada. Essa interação cria um ciclo de feedback onde o fluxo de fluido e a falha mecânica da rocha estão fortemente acoplados, uma dinâmica que o novo modelo captura com alta fidelidade.

Ao corrigir a base matemática desses modelos, os pesquisadores forneceram uma ferramenta mais confiável para engenheiros e cientistas que trabalham com meios porosos. A nova formulação garante que a transição entre diferentes tipos de falha da rocha seja suave e fisicamente realista, eliminando os saltos artificiais que prejudicavam as tentativas anteriores. Essa precisão é vital para aplicações que vão desde a otimização da produção de petróleo e gás até o design de armazenamento subterrâneo seguro para captura de carbono. O estudo demonstra que, mesmo em sistemas complexos onde sólidos e fluidos interagem, uma abordagem cuidadosa e consistente da física subjacente pode revelar a verdadeira natureza de como os materiais quebram. O modelo agora está pronto para ajudar a prever o comportamento da crosta terrestre com maior confiança, preenchendo a lacuna entre a mecânica teórica e a realidade plástica e desordenada do solo sob nossos pés.