Residual Saturation under Pressure-Controlled Drainage

Resumen Técnico: Saturación Residual bajo Drenaje Controlado por Presión

Planteamiento del Problema
El artículo aborda una cuestión fundamental aún no resuelta en la física de medios porosos: el comportamiento de la saturación residual bajo drenaje controlado por presión en el límite de tamaño infinito. Mientras que los modelos estándar de percolación por invasión (IP) predicen que la saturación de la fase invadida desaparece a medida que el tamaño del sistema tiende al infinito (dejando una fracción finita de fluido defensor atrapado solo debido a la naturaleza fractal del clúster de ruptura), no está claro si un proceso controlado por presión —que permite a la fase invadente acceder a todos los poros conectados a la entrada dada una presión capilar determinada— resulta en una saturación residual finita de la fase defensora en el límite termodinámico. Esta distinción es crítica para comprender fenómenos a escala de yacimiento, como la inundación de agua, donde surge la duda de si la saturación de petróleo residual permanece finita a medida que la región invadida se aproxima a las dimensiones del yacimiento.

Metodología
Los autores formulan el drenaje controlado por presión como Percolación de Enlaces con Atrapamiento (BPT) en un grafo de red de poros.

• Representación de Grafo: El medio poroso se modela como un grafo no dirigido donde los nodos representan regiones de poros y las aristas representan cuellos con radios de entrada asignados (umbrales capilares).
• Definición del Proceso: A diferencia de la IP estándar, que avanza cuello a cuello mediante la ruta más accesible localmente, la BPT opera bajo una restricción de presión global. A medida que la presión capilar aumenta, el "grafo activo" (compuesto por cuellos con umbrales de entrada por debajo de la presión actual) se expande.
• Mecanismo de Atrapamiento: La invasión se define por la conectividad con la entrada en el grafo activo. El atrapamiento ocurre cuando una porción de la fase defensora pierde la conectividad con el conjunto de salida. El algoritmo actualiza el sistema mediante: (1) el aumento de la presión para activar nuevos cuellos; (2) la identificación de la región invadida (poros activos conectados a la entrada); (3) la identificación de la región defensora (poros no invadidos conectados a la salida); y (4) la clasificación de todos los poros restantes desconectados como "atrapados".
• Simulaciones: El estudio compara la BPT frente a la Invasión por Percolación (IP) estándar en diversas redes 2D (cuadrada, triangular, Voronoi) y redes 3D (diamante, cúbica simple con diferentes números de coordinación). Las simulaciones utilizan radios de cuello distribuidos lognormalmente y promedios de ensamble para analizar los efectos de tamaño finito.

Contribuciones Clave y Resultados

1. Saturación Residual Finita en el Límite Infinito:
   La principal conclusión es que el drenaje controlado por presión conduce a un estado invadido de llenado de espacio donde la saturación residual de la fase defensora se aproxima a un límite finito y no nulo ($0 < S_\infty < 1$) a medida que $L \to \infty$. Esto contrasta marcadamente con la IP estándar, donde la saturación invadida desaparece en el momento de la ruptura ($S_{IP} \to 0$) porque el clúster de invasión es fractal y evade grandes regiones del fluido defensor.

2. Escalamiento de Tamaño Finito en Dos Dimensiones:
   En 2D, la desviación de la saturación residual de su límite termodinámico sigue una ley de escalamiento de tamaño finito:
   $$S_r(L) \approx S_\infty - \left(\frac{L_0}{L}\right)^\delta$$
   donde el exponente $\delta \approx 0.25$. Este exponente es universal entre diferentes tipos de redes (cuadrada, triangular, Voronoi) y es independiente de los detalles microscópicos. La convergencia lenta ($\Delta S_r \sim L^{-0.25}$) implica que alcanzar volúmenes elementales representativos en 2D requiere tamaños de sistema extremadamente grandes (por ejemplo, duplicar el tamaño del sistema reduce el sesgo de tamaño finito solo por un factor de $\approx 1/16$).

3. Efectos de Dimensionalidad y Número de Coordinación en Tres Dimensiones:
   En 3D, las correcciones de tamaño finito decaen significativamente más rápido, con un exponente $\delta \approx 0.75$. En consecuencia, el volumen representativo requerido para aproximar el límite de tamaño infinito es mucho menor en 3D que en 2D. Además, aunque la saturación residual asintótica sigue siendo finita en 3D, su valor depende del número de coordinación de la red; una mayor conectividad proporciona más rutas de invasión alternativas, reduciendo la fracción atrapada, aunque no la elimina por completo.

4. Distinción de la Invasión por Percolación:
   El artículo aclara que la diferencia entre la BPT y la IP no es simplemente una cuestión de la forma de la curva o la dinámica de ráfagas (bursts), sino una diferencia fundamental en los estados asintóticos. La IP produce un clúster de invasión fractal en el momento de la ruptura, mientras que la BPT produce una región invadida de llenado de espacio con una fracción atrapada finita que coexiste en el límite infinito.

Significancia
El artículo establece una conexión teórica directa entre la teoría de percolación y las relaciones presión-saturación al demostrar que la conectividad y el atrapamiento son los factores determinantes para la saturación residual bajo drenaje controlado por presión. Al mostrar que una fracción finita del espacio poroso permanece atrapada incluso en el límite de sistema infinito, el trabajo extiende la imagen estándar de la invasión por percolación más allá del estado de ruptura. Esto proporciona un marco riguroso para interpretar las leyes constitutivas macroscópicas y sugiere que el problema del "volumen elemental representativo" es altamente sensible a la dimensionalidad, siendo mucho más restrictivo en 2D que en 3D. Los resultados implican que, en escenarios controlados por presión, una fracción finita del fluido defensor es irrecuperable independientemente del tamaño del sistema, un hallazgo con implicaciones directas para la redistribución de fluidos en el subsuelo y los procesos de recuperación.