O fraturamento hidráulico em formações rochosas pode ser capturado elegantemente com o método de campo de fase. Vamos explorra uma abordagem no COMSOL que combina a interface PDE, a lei de Darcy e a mecânica dos sólidos, utilizando o campo de fase como um "bloco de montar" matemático.
A variável de campo de fase η (eta) representa o estado do material: η = 0 indica rocha intacta, η = 1 fratura completa, e valores intermediários descrevem a transição. O funcional de energia é definido com termos de duplo poço. A equação governante do campo de fase é:
ε^2 * η_tt - G_f * (1 - η) + 2 * (1 - η) * H = 0
onde ε é o parâmetro de regularização que controla a espessura da zona de transição, G_f é a energia de fratura, e H é uma variável histórica que armazena a maior tensão principle, impedindo a autocura da fratura. Esse histórico funciona como uma "memória" do dano.
Na lei de Darcy, a permeabilidade deve variar com a fratura. Uma expressão modificada é:
k = k0 * (1 + 500 * η^3)
O termo cúbico garante que a permeabilidade aumente suavemente no início da fratura e se torne alta quando completamente aberta. Essa relação não linear melhora a correspondência com dados experimentais de pressão.
No módulo de mecânica dos sólidos, o campo de fase é incorporado à lei constitutiva. Em vez de um decaimento linear, utiliza-se um decaimento exponencial para o módulo de Young:
E = E0 * exp(-5 * η)
Isso proporciona uma perda de rigidez mais suave e numericamente estável.
O acoplamento entre as equações requer uma atualização defasada dos pasos de tempo. Um esquema iterativo simples é:
- Calcular a atualização do campo de fase usando o campo de deslocamento atual.
- Atualizar os parâmetros do material com o novo campo de fase.
- Resolver novamente as equações mecânicas.
Repetir três iterações é geralmente suficiente para a convergência.
Para visualização, recomenda-se exibir simultaneamente o mapa de tensões principais e uma superfície isovalor de η = 0,5 (semitransparente em vermelho) para representar a frente da fratura. Adicionar vetores de gradiente de pressão ajuda a entender o fluxo do fluido.
Dicas de depuração: Se houver problemas de penetração da malha, aumente o comprimento característico ε para o dobro do tamanho do elemento. Em caso de não convergência, reformule o termo de gradiente de pressão na lei de Darcy como forma fraca.