Residual Saturation under Pressure-Controlled Drainage

Resumo Técnico: Saturação Residual sob Drenagem Controlada por Pressão

Definição do Problema
O artigo aborda uma questão fundamental e não resolvida na física de meios porosos: o comportamento da saturação residual sob drenagem controlada por pressão no limite de tamanho infinito. Embora os modelos padrão de percolação de invasão (IP) prevejam que a saturação da fase invadida desaparece conforme o tamanho do sistema se aproxima do infinito (deixando uma fração finita de fluido defensor retido apenas devido à natureza fractal do cluster de ruptura), permanece incerto se um processo controlado por pressão — que permite à fase invadida acessar todos os poros conectados à entrada em uma determinada limiar de capilaridade — resulta em uma saturação residual finita do fluido defensor no limite termodinâmico. Esta distinção é crítica para compreender fenômenos em escala de campo, como a lavagem de água (waterflooding), onde surge a questão se a saturação residual de óleo permanece finita à medida que a região invadida se aproxima das dimensões do reservatório.

Metodologia
Os autores formulam a drenagem controlada por pressão como Percolação de Ligação com Aprisionamento (BPT) em um grafo de rede de poros.

• Representação em Grafo: O meio poroso é modelado como um grafo não direcionado onde os nós representam regiões de poros e as arestas representam gargalos com raios de entrada atribuídos (limiares capilares).
• Definição do Processo: Diferente da IP padrão, que avança gargalo a gargalo via o caminho localmente mais acessível, a BPT opera sob uma restrição de pressão global. À medida que a pressão capilar aumenta, o "grafo ativo" (composto por gargalos com limiares de entrada abaixo da pressão atual) se expande.
• Mecanismo de Aprisionamento: A invasão é definida pela conectividade com a entrada no grafo ativo. O aprisionamento ocorre quando uma porção do fluido defensor perde a conectividade com o conjunto de saída. O algoritmo atualiza o sistema por: (1) aumentar a pressão para ativar novos gargalos; (2) identificar a região invadida (poros ativos conectados à entrada); (3) identificar a região defensora (poros conectados à saída que não foram invadidos); e (4) classificar todos os poros restantes desconectados como "aprisionados".
• Simulações: O estudo compara a BPT contra a Invasão de Percolação (IP) padrão em várias redes 2D (quadrada, triangular, Voronoi) e redes 3D (diamante, cúbica simples com diferentes números de coordenação). Simulações utilizam raios de gargalo distribuídos log-normalmente e média de conjunto para analisar efeitos de tamanho finito.

Principais Contribuições e Resultados

1. Saturação Residual Finita no Limite Infinito:
   A principal descoberta é que a drenagem controlada por pressão leva a um estado invadido de preenchimento de espaço onde a saturação residual do fluido defensor aproxima-se de um limite finito e não nulo ($0 < S_\infty < 1$) conforme o tamanho do sistema $L \to \infty$. Isso contrasta fortemente com a IP padrão, onde a saturação invadida desaparece no momento da ruptura ($S_{IP} \to 0$) porque o cluster de invasão é fractal e contorna grandes regiões do fluido defensor.

2. Escalonamento de Tamanho Finito em Dois Dímenisões:
   Em 2D, o desvio da saturação residual de seu limite termodinâmico segue uma lei de escalonamento de tamanho finito:
   $$S_r(L) \approx S_\infty - \left(\frac{L_0}{L}\right)^\delta$$
   onde o expoente $\delta \approx 0,25$. Este expoente é universal entre diferentes tipos de redes (quadrada, triangular, Voronoi) e é independente de detalhes microscópicos. A convergência lenta ($\Delta S_r \sim L^{-0,25}$) implica que alcançar volumes elementares representativos em 2D requer tamanhos de sistema extremamente grandes (por exemplo, dobrar o tamanho do sistema reduz o viés de tamanho finito por apenas um fator de $\approx 1/16$).

3. Efeitos de Dimensionalidade e Número de Coordenação em Três Dimensões:
   Em 3D, as correções de tamanho finito decaem significativamente mais rápido, com um expoente $\delta \approx 0,75$. Consequentemente, o volume representativo necessário para aproximar o limite de tamanho infinito é muito menor em 3D do que em 2D. Além disso, embora a saturação residual assintótica permaneça finita em 3D, seu valor depende do número de coordenação da rede; uma maior conectividade fornece mais caminhos alternativos de invasão, reduzindo a fração aprisionada, embora não a elimine totalmente.

4. Distinção da Percolação de Invasão:
   O artigo esclarece que a diferença entre BPT e IP não é meramente uma questão de forma da curva ou dinâmica de surtos (bursts), mas uma diferença fundamental nos estados assintóticos. A IP produz um cluster de invasão fractal no momento da ruptura, enquanto a BPT produz uma região invadida de preenchimento de espaço com uma fração aprisionada finita coexistindo no limite infinito.

Significância
O artigo estabelece uma conexão teórica direta entre a teoria de percolação e as relações pressão-saturação ao demonstrar que a conectividade e o aprisionamento são os fatores determinantes para a saturação residual sob drenagem controlada por pressão. Ao mostrar que uma fração finita do espaço poroso permanece aprisionada mesmo no limite de sistema infinito, o trabalho estende a imagem padrão da percolação de invasão além do estado de ruptura. Isso fornece um arcabouço rigoroso para interpretar leis constitutivas macroscópicas e sugere que o problema do "volume elementar representativo" é altamente sensível à dimensionalidade, sendo muito mais restritivo em 2D do que em 3D. Os resultados implicam que, em cenários controlados por pressão, uma fração finita do fluido defensor é irrecuperável, independentemente do tamanho do sistema, um achado com implicações diretas para a redistribuição de fluidos de subsuperfície e processos de recuperação.