Linear Stability Analysis of Two-phase, Two-Component Flow in Porous Media

Resumo Técnico: Análise de Estabilidade Linear de Fluxo Bifásico e Bicomponente em Meios Porosos

Definição do Problema
As instabilidades de dedos viscosos (viscous fingering) durante o deslocamento de fluidos em meios porosos comprometem significativamente a eficiência de aplicações como a recuperação avançada de petróleo, sequestro de CO2 e remediação de águas subterrâneas. Embora exista uma vasta literatura sobre análise de estabilidade linear para deslocamentos totalmente imiscíveis e totalmente miscíveis, o caso intermediário de fluxo parcialmente miscível — onde ocorre transferência de massa limitada entre as fases — permanece amplamente inexplorado. Especificamente, há uma carência de estruturas analíticas que deem conta da complexa interação entre gravidade, forças capilares, dispersão mecânica e transferência de massa interfásica em um sistema de duas fases e dois componentes. Este estudo aborda a previsão de modos de perturbação instáveis em tais sistemas parcialmente miscíveis, focando no cenário em que um fluido gasoso desloca um líquido.

Metodologia
Os autores desenvolvem um modelo matemático para fluxo bidimensional em um meio poroso isotrópico, assumindo o equilíbrio termodinâmico para uma mistura binária (componentes "a" e "b"). A metodologia procede através das seguintes etapas:

1. Equações Governantes: O modelo deriva equações de balanço de massa para um sistema de duas fases e dois componentes, incorporando a lei de Darcy, fluxo fracionário, pressão capilar e dispersão mecânica. Uma premissa simplificadora fundamental é que a densidade da fase líquida é constante na região de fase única, enquanto a viscosidade varia com a composição.
2. Formulação do Estado Base: O estado base é definido como uma solução de onda viajante unidimensional (frente de choque), derivada da abordagem de Buckley-Leverett, adaptada para transferência de massa. A configuração do choque é calculada usando o método das características, considerando a descontinuidade nas funções de fluxo na transição da região bifásica para a região de líquido puro.
3. Análise de Estabilidade Linear: Pequenas perturbações de amplitude de onda são sobrepostas ao estado base. As equações governantes são linearizadas para formar um problema de autovalor, onde o autovalor representa a taxa de crescimento da perturbação ($\sigma$) e os autofunções representam os perfis de perturbação.
4. Solução Numérica: Um desafio crítico identificado é a descontinuidade nas derivadas das autofunções no ponto de transição entre o fluxo bifásico e o de fase única. Para resolver isso, os autores derivam condições de salto específicas para as derivadas no ponto de interface. O problema de autovalor é resolvido numericamente utilizando o método do Problema de Valor Inicial Casado (MIVP - Matched Initial Value Problem). Isso envolve a integração das equações diferenciais ordinárias linearizadas a partir das fronteiras de campo distante em direção ao ponto de transição e o ajuste das soluções utilizando as condições de salto derivadas.

Principais Contribuições

• Estrutura Inovadora: Este trabalho apresenta, tanto quanto aos autores sabem, a primeira análise de estabilidade linear de um sistema genuíno de duas fases e dois componentes que considera explicitamente a transferência de massa interfásica.
• Condições de Salto: A derivação das condições de salto para as derivadas das autofunções na descontinuidade da transição de fase permite a solução precisa do problema de autovalor em sistemas onde as funções de fluxo são descontínuas.
• Estudo Abrangente de Parâmetros: O estudo investiga sistematicamente os efeitos da composição inicial, dispersividades (longitudinal e transversal), níveis de miscibilidade (via coeficientes de interação binária), razões de viscosidade e gravidade sobre a estabilidade da frente de deslocamento.

Resultos

• Efeito Estabilizador da Transferência de Massa: O principal achado é que a transferência de massa possui um efeito predominantemente estabilizador no deslocamento. Ela reduz a taxa de crescimento máxima ($\sigma_{max}$) ao reduzir a razão de viscosidade de equilíbrio ($M_e$) e alterar as propriedades do choque (reduzindo a velocidade do choque $v_s$ e aumentando a saturação de gás $S_{gs}$). Esta estabilização é particularmente pronunciada para altos contrastes de viscosidade.
• Interações Gravitacionais: Em deslocamentos ascendentes, a transferência de massa amortece a instabilidade induzida pela gravidade. O contraste de densidade que impulsiona a instabilidade gravitacional também potencializa a transferência de massa, que atua neutralizando o efeito desestabilizador. Em deslocamentos descendentes, tanto a gravidade quanto a transferência de massa atuam como mecanismos estabilizadores.
• Não Linearidade da Dispersão: O estudo identifica uma relação não linear entre o coeficiente de dispersão longitudinal adimensional ($D^*_{\ell,xx}$) e a instabilidade. Existe um valor crítico "mais perigoso" de $D^*_{\ell,xx}$ onde tanto a taxa de crescimento quanto o número de onda de corte ($n_{cut}$) são maximizados. Isso sugere uma interação complexa onde níveis intermediários de dispersão mecânica, em relação às forças capilares, podem de fato aumentar a instabilidade em comparação com regimes de dispersão muito baixa ou muito alta.
• Número de Onda de Corte: Embora a transferência de massa reduza significamente a taxa de crescimento, seu efeito no número de onda de corte é menos pronunciado e mais complexo, resultando às vezes em valores de $n_{cut}$ mais elevados para casos parcialmente miscíveis em comparação com casos imiscíveis sob condições específicas.

Significância e Alegações
O artigo afirma preencher a lacuna entre os casos limites bem estudados de deslocamento imiscível e miscível. Ao demonstrar que a transferência de massa atua como um mecanismo estabilizador, os resultados sugerem que condições parcialmente miscíveis podem oferecer frentes de deslocamento mais robustas do que cenários puramente imiscíveis, particularmente em ambientes de alto contraste de viscosidade. A identificação do coeficiente de dispersão crítico destaca uma interação complexa anteriormente não relatada entre forças capilares e dispersão mecânica no contexto de dedos viscosos. Os autores observam que, embora seu modelo simplifique a física de sistemas multicomponentes, ele fornece uma estrutura necessária para compreender os mecanismos fundamentais de estabilidade em fluxos parcialmente miscíveis. O estudo reconhece limitações, como a suposição de densidade constante do líquido na região de fase única e curvas de permeabilidade relativa fixas, sugerindo estas como áreas para refinamento futuro.