Mathematical Modeling of Salt Precipitation and Multi-Phase Flow in High Enthalpy Fractured Geothermal Systems

Este artigo apresenta um novo modelo de fluxo composicional de código aberto implementado no framework PorePy que simula o fluxo multifásico não isotérmico e a precipitação de halita em reservatórios geotérmicos fraturados de alta entalpia, utilizando uma formulação de variável primária robusta e uma abordagem de matriz-fratura discreta para prever com precisão o dano de permeabilidade e os desafios operacionais.

O essencial

Imagine um reservatório subterrâneo de alta temperatura como uma panela de pressão gigante e natural, cheia de água salgada superquente. Esta é a fonte de energia geotérmica. Para extrair essa energia, engenheiros perfuram poços e bombeiam água fria para dentro para empurrar a água quente de volta para cima. No entanto, esse processo é complicado devido ao "sal" na água.

Pense no sal (halita) como o açúcar em uma xícara de chá quente. Se você esfriar o chá ou deixar parte da água evaporar, o açúcar não consegue permanecer dissolvido e começa a se transformar novamente em cristais sólidos. Em um poço geotérmico, isso acontece quando a água quente esfria perto de um poço de injeção ou evapora perto de um poço de produção. O resultado? Cristais de sal sólidos se formam e entopem os pequenos buracos nas rochas e as fendas (fraturas) que permitem o fluxo da água. É como um congestionamento causado por cristais de açúcar bloqueando uma rodovia.

Este artigo apresenta uma nova ferramenta de simulação computacional projetada para prever exatamente onde e como esses "engarrafamentos de açúcar" ocorrerão em rochas subterrâneas complexas e rachadas.

Aqui está uma análise de como a ferramenta funciona e o que ela descobriu, usando analogias simples:

1. O Mapa: Vendo as Rachaduras Claramente

Os mapas tradicionais de rochas subterrâneas costumam suavizar tudo, tratando a rocha como uma esponja sólida. Mas, na realidade, a água flui principalmente através de uma rede de rachaduras, como a água fluindo através de uma calçada rachada em vez de através do próprio concreto.

• A Inovação: Este novo modelo utiliza uma abordagem de "Matriz de Fratura Discreta". Imagine desenhar as rachaduras como linhas finas e distintas em um mapa, em vez de apenas misturá-las ao fundo. Isso permite que o computador veja exatamente como as rachaduras se conectam (ou não se conectam) e como o sal pode entupir uma rachadura específica versus a rocha ao redor.

2. O Motor: Um "Controle Remoto Universal" para a Física

Simular água fervente, vapor e sal sólido ao mesmo tempo é incrivelmente difícil para os computadores. Normalmente, o computador tem que mudar constantemente seu "modo" (ex: "Ok, agora é líquido; agora é gás; agora é sólido"), o que pode fazer o cálculo travar ou ficar preso.

• A Inovação: Os autores criaram um sistema "unificado". Pense nisso como um controle remoto universal que funciona para todos os dispositivos sem precisar trocar as pilhas ou os modos. O modelo usa três "ponteiros" fixos (Pressão, Energia Térmica e Quantidade de Sal) que permanecem os mesmos, quer a água seja líquida, vapor ou esteja se transformando em sal sólido. Isso torna a simulação muito mais suave e estável, permitindo que ela lide com a troca caótica entre estados sem quebrar.

3. O Truque de Velocidade: A "Folha de Cola"

Calcular a física exata da água salgada em altas temperaturas geralmente exige que o computador resolva quebra-cabeças matemáticos complexos repetidamente, o que é muito lento.

• A Inovação: A equipe criou uma "folha de cola" pré-calculada (tabela de consulta). Antes da simulação começar, eles calcularam todos os resultados possíveis sobre como o sal se comporta sob diferentes condições e os armazenaram. Durante a simulação, em vez de resolver a matemática difícil toda vez, o computador apenas consulta a resposta na folha. Isso faz com que a simulação rode muito mais rápido, mantendo a precisão.

4. O Efeito de Entupimento: O "Encolhimento do Poro"

Quando os cristais de sal se formam, eles ocupam espaço.

• A Inovação: O modelo reduz automaticamente o tamanho dos "canos" (porosidade e permeabilidade) conforme o sal se acumula. Ele utiliza uma regra (Kozeny-Carman) que diz: "Se o sal preencher 10% do buraco, o cano fica significativamente mais estreito". Isso permite que o modelo preveja como o fluxo diminuirá ou parará totalmente à medida que o "engarrafamento de açúcar" piora.

O Que as Simulações Mostraram

A equipe testou esta ferramenta em dois cenários principais:

Cenário A: A Rodovia Quebrada (Rachaduras Desconectadas)

• Configuração: Imagine um reservatório onde as rachaduras não se conectam umas às outras; a água tem que espremer-se através da rocha sólida entre elas.
• Resultado: Quando bombearam água fria, a água quente perto do poço de produção ferveu rapidamente. Isso fez com que o sal cristalizasse e entupisse a rocha logo ao redor do poço.
• A Reviravolta: Se eles bombeassem água mais rápido, o entupimento piorava muito, e a produção de energia caía significamente. O modelo mostrou que o "congestionamento" acontecia principalmente na rocha perto do poço, não apenas nas rachaduras.

Cenário B: A Rodovia Conectada (Rachaduras Conectadas)

• Configuração: Imagine um reservatório onde as rachaduras formam uma rodovia contínua de alta velocidade do poço de injeção ao poço de produção.
• Resultado: A água fria correu pelas rachaduras rapidamente. Como ela se movia muito rápido e permanecia fria, ela na verdade dissolveu o sal perto do poço de produção em vez de entupi-lo!
• A Reviravolta: A precipitação de sal moveu-se para um lugar diferente — exatamente na borda da zona de água fria — em vez de entupir o poço em si. Isso sugere que ter uma rede de rachaduras conectada pode, na verdade, proteger o poço contra o entupimento, embora mude o local onde o sal se acumula.

A Conclusão

Este artigo apresenta uma nova ferramenta de software de código aberto que ajuda engenheiros a entender a dança complexa entre calor, pressão e sal em poços geotérmicos. Ao mapear com precisão como as rachaduras se conectam e como o sal as entope, a ferramenta pode ajudar a prever:

1. Onde os poços podem ser bloqueados pelo sal.
2. Quanta energia pode ser extraída com segurança antes que os "canos" fiquem entupidos.
3. Se o layout das rachaduras subterrâneas ajudará ou prejudicará o processo de produção.

Os autores verificaram sua ferramenta em relação a um padrão estabelecido pela indústria e descobriram que ela coincidia perfeitamente, provando que é uma maneira confiável de simular esses ambientes de rocha rachada, salgada e de alta temperatura.

Resumo técnico

Resumo Técnico: Modelo de Fluxo Composicional Unificado para Simulação de Sistemas de Alta Entalpia Multifásicos e Fraturados

Definição do Problema
A simulação de reservatórios geotérmicos fraturados de alta entalpia apresenta desafios significativos devido ao complexo acoplamento de fluxo multifásico, multicomponente e não isotérmico, termodinâmica fortemente não linear e a influência dominante de redes de fraturas. Uma questão operacional crítica nestes sistemas é a incrustação mineral, especificamente a precipitação de halita (NaCl). Durante a produção, o declínio de pressão induz a ebulição, concentrando a salmoura residual e excedendo os limites de solubilidade, enquanto a injeção de água fria reduz a solubilidade local. Ambos os mecanismos levam à precipitação de sal nos espaços porosos e fraturas, prejudicando severamente a porosidade e a permeabilidade, reduzindo a produtividade dos poços e potencialmente causando o bloqueio do poço. Os modelos de contínuo existentes frequentemente falham em capturar a precipitação localizada e a evolução da permeabilidade na escala das fraturas, enquanto os simuladores composicionais tradicionais, que dependem de formulações de troca de variáveis (variable-switching), enfrentam dificuldades de estabilidade numérica próximo às fronteiras de fase.

Metodologia
Os autores apresentam um novo modelo de fluxo composicional unificado implementado dentro do framework de código aberto PorePy. A metodologia integra três componentes principais:

1. Formulação de Variável Persistente: Diferente das abordagens tradicionais que trocam as variáveis primárias com base na presença de fase, este modelo utiliza um conjunto fixo de variáveis primárias — pressão ($p$), entalpia específica ($h$) e fração de massa de sal global ($z$) — através de todas as configurações de fase (líquido, vapor e halita sólida). Esta abordagem lida naturalmente com transições de fase sem a necessidade de troca manual, aumentando a robustez numérica e a estabilidade.
2. Abordagem de Matriz-Fratura Discreta (DFM): O domínio físico é representado usando uma geometria de dimensão mista, onde as fraturas e suas interseções são tratadas como variedades de dimensão inferior (1D e 0D) inseridas na matriz porosa 2D. Esta representação explícita preserva a conectividade e a heterogeneidade das fraturas, permitendo a resolução de processos localizados, como a redução da permeabilidade induzida pela precipitação.
3. Fechamento Termodinâmico e Linearização: O modelo emprega as correlações abrangentes de Driesner e Christoph (2007) para a termodinâmica H2O–NaCl, cobrindo condições magmáticas. Para reduzir o custo computacional, o modelo utiliza uma estratégia de Linearização Baseada em Operador (OBL). As propriedades termodinâmicas e suas derivadas são pré-computadas em uma grade estruturada no espaço das variáveis primárias e recuperadas via interpolação multilinear durante a simulação, eliminando cálculos de separação de fase (flash) dispendiosos realizados em tempo real.
4. Evolução Dinâmica da Permeabilidade e da Abertura: O modelo incorpora relações do tipo Kozeny-Carman para atualizar dinamicamente a porosidade e a permeabilidade da matriz com base na saturação de halita. Para as fraturas, a abertura evolui de acordo com uma relação de potência dependente da saturação de halita, que por sua vez atualiza a permeabilidade da fratura via lei cúbica.

As equações governantes (balanço de massa e energia) são discretizadas usando um esquema de volume finito centrado na célula com Aproximação de Fluxo de Múltiplos Pontos (MPFA) para termos difusivos e upwinding para termos advectivos. O sistema não linear resultante é resolvido usando o método de Newton com uma busca de linha Armijo de retrocesso (backtracking) e passos de tempo adaptativos.

Principais Resultados
O modelo foi verificado e aplicado através de diversos experimentos numéricos:

• Verificação (Benchmark 1D): O modelo foi verificado contra o simulador de código fechado estabelecido CSMP++ utilizando um benchmark de dissolução de sal 1D. Os resultados mostraram forte concordância através de quatro regiões de fase distintas (vapor+halita, vapor+líquido, fase líquida única e líquido+halita) e suas transições, validando a formulação de variável persistente e a implementação termodinâmica.
• Rede de Fraturas Desconectada (2D): Em um reservatório com fraturas isoladas, o fluxo é dominado pela matriz. A simulação demonstrou que o declínio de pressão no poço de produção induz a ebulição, levando à precipitação significativa de halita na matriz circundante e nos segmentos de fratura adjacentes. Esta precipitação reduziu a permeabilidade da matriz em até 36% ($K/K_0 \approx 0,64$) e alterou as linhas de corrente do fluxo. O aumento da taxa de injeção intensificou a ebulição próxima ao poço, levando a uma maior saturação de halita ($s_{hal} \approx 0,46$) e uma redução de permeabilidade mais severa ($K/K_0 \approx 0,29$), impactando significativamente as taxas de produção.
• Rede de Fraturas Conectada (2D): Em um reservatório com uma cadeia de fraturas contínua conectando os poços de injeção e produção, o transporte tornou-se dominado pelas fraturas. A rede conectada facilitou a comunicação rápida de pressão e o transporte de fluido de baixa salinidade para o poço de produção. Consequentemente, embora a ebulição ocorresse próximo ao produtor, o fluido que chegava permanecia suficientemente diluído para evitar a saturação de halita local. Em vez de precipitação no produtor, observou-se dissolução, e a precipitação foi redistribuída para a matriz à frente da frente de injeção e ao longo da cadeia de fraturas devido aos efeitos de resfriamento. Esta configuração evitou a severa deterioração da permeabilidade vista no caso desconectado.

Significância e Alegações
O artigo afirma fornecer uma ferramenta numérica robusta e de código aberto para analisar o complexo transporte de calor e massa com efeitos de incrustação mineral em sistemas geotérmicos de alta entalpia fraturados. Sua principal significância reside em:

• Estabilidade Numérica: A formulação unificada de variável persistente evita as dificuldades de convergência associadas à troca de variáveis próximas às fronteiras de fase.
• Fidelidade Geométrica: A abordagem DFM resolve explicitamente a conectividade das fraturas, o que os autores demonstram ser um fator crítico que controla a distribuição espacial da precipitação e a resultante redução de transmissibilidade.
• Visão Operacional: Os resultados destacam que a conectividade das fraturas dita os desafios operacionais; redes desconectadas são propensas ao entupimento próximo ao poço e perda de produtividade, enquanto redes conectadas podem mitigar esses problemas ao manter a diluição do fluido no poço de produção.
• Reprodutibilidade: O trabalho fornece o modelo numérico completo e o fluxo de execução via um contêiner Docker, facilitando a reprodutibilidade e o desenvolvimento adicional dentro da comunidade geotérmica.

Os autores concluem que seu modelo prevê efetivamente os padrões de precipitação de halita e seu impacto no desempenho do reservatório, oferecendo uma ferramenta valiosa para compreender a interação entre fluxo, termodinâmica e incrustação mineral em reservatórios geotérmicos fraturados.