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

Resumen Técnico: Análisis de Estabilidad Lineal de Flujo Bifásico y Bicomponente en Medios Porosos

Planteamiento del Problema
Las inestabilidades de dedos viscosos durante el desplazamiento de fluidos en medios porosos comprometen significativamente la eficiencia de aplicaciones tales como la recuperación mejorada de petróleo, el secuestro de CO2 y la remediación de aguas subterráneas. Si bien existe una amplia literatura sobre el análisis de estabilidad lineal para desplazamientos totalmente inmiscibles y totalmente miscibles, el caso intermedio del flujo parcialmente miscible —donde ocurre una transferencia de masa limitada entre fases— permanece en gran medida inexplorado. Específicamente, existe una falta de marcos analíticos que den cuenta de la compleja interacción entre la gravedad, las fuerzas capilares, la dispersión mecánica y la transferencia de masa interfásica en un sistema de dos fases y dos componentes. Este estudio aborda la predicción de los modos de perturbación inestables en tales sistemas parcialmente miscibles, centrándose en el escenario donde un fluido gaseoso desplaza a un líquido.

Metodología
Los autores desarrollan un modelo matemático para el flujo bidimensional en un medio poroso isotrópico, asumiendo el equilibrio termodinámico para una mezcla binaria (componentes "a" y "b"). La metodología procede a través de los siguientes pasos:

1. Ecuaciones de Gobierno: El modelo deriva ecuaciones de balance de masa para un sistema de dos fases y dos componentes, incorporando la ley de Darcy, el flujo fraccional, la presión capilar y la dispersión mecánica. Un supuesto simplificador clave es que la densidad de la fase líquida es constante en la región de una sola fase, mientras que la viscosidad varía con la composición.
2. Formulación del Estado Base: El estado base se define como una solución de onda viajera unidimensional (frente de choque) derivada del enfoque de Buckley-Leverett, adaptado para la transferencia de masa. La configuración del choque se calcula utilizando el método de las características, teniendo en cuenta la discontinuidad en las funciones de flujo en la transición de la región de dos fases a la región de líquido puro.
3. Análisis de Estabilidad Lineal: Se superponen perturbaciones de pequeña amplitud de tipo ondulatorio sobre el estado base. Las ecuaciones de gobierno se linealizan para formar un problema de autovalores donde el autovalor representa la tasa de crecimiento de la perturbación ($\sigma$) y las autofunciones representan los perfiles de la perturbación.
4. Solución Numérica: Un desafío crítico identificado es la discontinuidad en las derivadas de las autofunciones en el punto de transición entre el flujo de dos fases y el de una sola fase. Para resolver esto, los autores derivan condiciones de salto específicas para las derivadas en esta interfaz. El problema de autovalores se resuelve numéricamente utilizando el método de Problema de Valor Inicial Emparejado (MIVP, por sus siglas en inglés). Esto implica integrar las ecuaciones diferenciales ordinarias linealizadas desde los límites de campo lejano hacia el punto de transición y emparejar las soluciones utilizando las condiciones de salto derivadas.

Contribuciones Clave

• Marco Novedoso: Este trabajo presenta, según los autores, el primer análisis de estabilidad lineal de un sistema genuino de dos fases y dos componentes que explica explícitamente la transferencia de masa interfásica.
• Condiciones de Salto: La derivación de las condiciones de salto para las derivadas de las autofunciones en la discontinuidad de la transición de fase permite la solución precisa del problema de autovalores en sistemas donde las funciones de flujo son discontinuas.
• Estudio Exhaustivo de Parámetros: El estudio investiga sistemáticamente los efectos de la composición inicial, las dispersividades (longitudinal y transversal), los niveles de miscibilidad (vía coeficientes de interacción binaria), las relaciones de viscosidad y la gravedad sobre la estabilidad del frente de desplazamiento.

Resultados

• Efecto Estabilizador de la Transferencia de Masa: El hallazgo principal es que la transferencia de masa tiene un efecto predominantemente estabilizador en el desplazamiento. Reduce la tasa de crecimiento máxima ($\sigma_{max}$) al disminuir la relación de viscosidad de equilibrio ($M_e$) y alterar las propiedades del choque (reduciendo la velocidad del choque $v_s$ y aumentando la saturación de gas $S_{gs}$). Esta estabilización es particularmente pronunciada para contrastes de viscosidad elevados.
• Interacciones de Gravedad: En desplazamientos ascendentes, la transferencia de masa amortigua la inestabilidad inducida por la gravedad. El contraste de densidad que impulsa la inestabilidad gravitacional también potencia la transferencia de masa, la cual contrarresta el efecto desestabilizador. En desplazamientos descendentes, tanto la gravedad como la transferencia de masa actúan como mecanismos estabilizadores.
• No Linealidad de la Dispersión: El estudio identifica una relación no lineal entre el coeficiente de dispersión longitudinal adimensional ($D^*_{\ell,xx}$) y la inestabilidad. Existe un valor crítico "más peligroso" de $D^*_{\ell,xx}$ donde tanto la tasa de crecimiento como el número de onda de corte ($n_{cut}$) se maximizan. Esto sugiere una interacción compleja donde niveles intermedios de dispersión mecánica, en relación con las fuerzas capilares, pueden de hecho aumentar la inestabilidad en comparación con regímenes de dispersión tanto muy baja como muy alta.
• Número de Onda de Corte: Aunque la transferencia de masa reduce significativamente la tasa de crecimiento, su efecto sobre el número de onda de corte es menos pronunciado y más complejo, resultando a veces en valores de $n_{cut}$ más altos para casos parcialmente miscibles en comparación con los inmiscibles bajo condiciones específicas.

Significancia y Reivindicaciones
El artículo afirma cerrar la brecha entre los casos límite bien estudiados de desplazamiento inmiscible y miscible. Al demostrar que la transferencia de masa actúa como un mecanismo estabilizador, los resultados sugieren que las condiciones parcialmente miscibles pueden ofrecer frentes de desplazamiento más robustos que los escenarios puramente inmiscibles, particularmente en entornos de alto contraste de viscosidad. La identificación del valor crítico del coeficiente de dispersión resalta una interacción compleja previamente no reportada entre las fuerzas capilares y la dispersión mecánica en el contexto de los dedos viscosos. Los autores señalan que, si bien su modelo simplifica la física de los sistemas multicomponentes, proporciona un marco necesario para comprender los mecanismos fundamentales de estabilidad en flujos parcialmente miscibles. El estudio reconoce limitaciones, como el supuesto de densidad de líquido constante en la región de una sola fase y curvas de permeabilidad relativa fijas, sugiriendo estas como áreas para futuras refinaciones.