⚡ ingeniería eléctrica

Modeling Capillary Heterogeneity in Geological CO2 Storage via Finite-volume Discretization with Interface Conditions

Este artículo presenta un esquema de volúmenes finitos mejorado llamado PPU-C, el cual incorpora condiciones de interfaz discretas para modelar con precisión la heterogeneidad capilar en simulaciones de almacenamiento geológico de CO2, logrando mejoras significativas en la precisión de la predicción del volumen de atrapamiento con menos del 25% de sobrecarga computacional en comparación con los métodos estándar.

Autores originales: Ammar Alali, Bradley Mallison, Omar Duran, Hamdi Tchelepi

Publicado 2026-10-22

📖 6 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Ammar Alali, Bradley Mallison, Omar Duran, Hamdi Tchelepi

Artículo original bajo licencia CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ Esta es una explicación generada por IA del artículo a continuación. No ha sido escrita ni avalada por los autores. Para mayor precisión técnica, consulte el artículo original. Leer descargo de responsabilidad completo

Muy bajo tierra, en las rocas porosas que yacen bajo nuestros pies, se está llevando a cabo una lucha silenciosa por el espacio. Este es el reino del almacenamiento geológico de carbono, una estrategia en la que los científicos bombean dióxido de carbono supercrítico hacia la tierra para mantenerlo fuera de la atmósfera. El objetivo es atrapar este gas de forma permanente, evitando que contribuya al cambio climático. Sin embargo, las formaciones rocosas utilizadas para este almacenamiento no son bloques uniformes de piedra; son mosaicos complejos de diferentes materiales, muy parecidos a un pastel de capas hecho de rebanadas alternas de bizcocho y arcilla densa. Cuando el gas inyectado se mueve a través de estas capas, encuentra cambios repentinos en las propiedades de la roca. Estos cambios crean barreras invisibles que pueden detener el gas en su camino o forzarlo a tomar rutas inesperadas. Comprender exactamente cómo se comporta el gas cuando golpea estos límites es crítico. Si no podemos predecir hacia dónde irá el gas, no podemos estar seguros de cuánto de él permanecerá atrapado de forma segura, o si podría filtrarse de nuevo hacia la superficie.

El desafío reside en la física de la interfaz entre estos diferentes tipos de roca. A medida que el gas empuja hacia arriba, impulsado por su propia flotabilidad, se encuentra con una capa de roca que es más difícil de penetrar. En el mundo real, el gas se detendrá y se acumulará debajo de esta barrera, o logrará atravesarla si tiene suficiente presión. Los modelos informáticos estándar utilizados para simular este proceso suelen tener dificultades para lograr esto con precisión. Para ver los detalles de cómo el gas se queda atrapado, estos modelos suelen requerir una cuadrícula increíblemente fina, dividiendo el espacio subterráneo en millones de celdas diminutas y microscópicas. Este nivel de detalle es computacionalmente costoso y a menudo impráctico para las vastas escalas de los sitios de almacenamiento del mundo real. Sin este detalle extremo, los modelos tienden a suavizar los bordes rugosos de las capas de roca, lo que los lleva a predecir que el gas invade nuevas áreas con demasiada facilidad y se queda atrapado demasiado poco.

En un estudio reciente, investigadores de la Universidad de Stanford desarrollaron una nueva forma de manejar estos límites complejos sin necesidad de dividir el problema en millones de piezas diminutas. Crearon un esquema numérico especializado, al que llaman PPU-C, diseñado para respetar las reglas físicas en la interfaz entre diferentes tipos de roca. En lugar de obligar al ordenador a observar cada poro microscópico, este método introduce un conjunto específico de reglas sobre lo que sucede exactamente donde se encuentran dos rocas diferentes. Trata la frontera como un lugar distinto donde la presión y el flujo del gas deben equilibrarse perfectamente, tal como ocurre en la naturaleza. Al resolver un problema pequeño y local en cada frontera, el modelo puede determinar con precisión si el gas queda atrapado o si está fluyendo a través de ella, incluso cuando las celdas de la cuadrícula utilizadas en la simulación son relativamente grandes.

Los investigadores probaron este nuevo enfoque contra los métodos estándar utilizando una serie de escenarios cada vez más complejos. Comenzaron con una columna vertical simple de dos capas de roca, donde se inyectó dióxido de carbono en la parte inferior y se permitió que ascendiera. En estas pruebas, el modelo estándar predijo que el gas entraría apresuradamente en la capa superior mucho antes de lo que debería, fallando al capturar la cantidad de gas que realmente era retenida. El nuevo método, sin embargo, identificó correctamente la barrera y mostró el gas acumulándose debajo de ella, coincidiendo con el comportamiento físico esperado incluso en una cuadrícula gruesa. Cuando los investigadores añadieron la complejidad de la histéresis —un fenómeno en el que la roca recuerda sus niveles de saturación pasados y se comporta de manera diferente dependiendo de si el gas se está moviendo hacia adentro o hacia afuera— el nuevo método continuó superando al anterior. Calculó con precisión el volumen de gas que quedaría atrapado permanentemente, una métrica crucial para determinar la seguridad y la capacidad de un sitio de almacenamiento.

El estudio no se detuvo en columnas simples. El equipo expandió sus pruebas a espacios tridimensionales y luego a una sección bidimensional de una formación geológica real basada en un estándar de la industria. En estos entornos más grandes y realistas, el modelo estándar mostró nuevamente una tendencia a permitir que el gas se extendiera demasiado lejos y demasiado rápido, subestimando la cantidad de gas que sería almacenado de forma segura. El nuevo método mantuvo el gas en su lugar apropiado, mostrando una ruta de migración más realista y un volumen significemente mayor de carbono atrapado. Quizás lo más importante para la aplicación práctica, los investigadores descubrieron que este aumento en la precisión no conllevó un costo prohibitivo. Aunque el nuevo método requería un poco más de potencia de cálculo, el incremento fue modesto, nunca superando el veinticinco por ciento en sus pruebas. En algunos casos, debido a que el nuevo método era más estable y requería menos intentos para resolver las ecuaciones, de hecho funcionaba más rápido que el enfoque estándar.

Los hallazgos sugieren que, al prestar más atención a las reglas específicas que gobiernan los límites entre diferentes rocas, podemos simular el comportamiento de los sitios de almacenamiento de carbono con mucha mayor confianza. Esto significa que los ingenieros y científicos pueden diseñar proyectos de almacenamiento que sean tanto más seguros como más eficientes, sabiendo exactamente cuánto gas será retenido por la geología natural. La nueva técnica ofrece una forma de modelar estos entornos subterráneos complejos sin la necesidad de cuadrículas imposiblemente finas, haciendo que las simulaciones de alta fidelidad sean accesibles para proyectos a gran escala. Mientras el mundo busca expandir el almacenamiento geológico de carbono para cumplir con los objetivos climáticos, contar con herramientas que puedan predecir con precisión el destino del gas inyectado es esencial. Este trabajo proporciona un camino práctico hacia adelante, asegurando que la promesa del almacenamiento subterráneo esté respaldada por una comprensión clara y fiable de la física en juego.

Resumen Técnico: Modelado de la Heterogeneidad Capilar en el Almacenamiento Geológico de CO2 mediante la Discretización de Volumen Finito con Condiciones de Interfaz

Planteamiento del Problema
La heterogeneidad capilar es un factor crítico para determinar las rutas de migración de la pluma de CO2 y el volumen de CO2 atrapado de forma segura en formaciones geológicas. Sin embargo, los esquemas de simulación numérica estándar, como el método de Upwinding de Potencial de Fase (PPU, por sus siglas en inglés), tienen dificultades para resolver con precisión el atrapamiento impulsado por capilaridad en las interfaces de materiales sin un refinamiento excesivo de la malla. En medios porosos heterogéneos, la presión capilar y la permeabilidad relativa son espacialmente discontinuas debido a cambios abruptos en las propiedades de la roca (por ejemplo, el tamaño de los poros y la mojabilidad). Las aproximaciones de flujo estándar a menudo fallan al intentar satisfacer la continuidad de la presión capilar discreta en estas interfaces en mallas gruesas prácticas, lo que genera errores significativos en la predicción de los volúmenes de atrapamiento por residualidad y heterogeneidad capilar. Aunque se han utilizado condiciones de interfaz discretas en otros contextos (por ejemplo, atrapamiento de DNAPL y sistemas matriz-fractura), aún no son comunes en los simuladores de volumen finito totalmente implícitos para el Almacenamiento Geológico de Carbono (GCS) a gran escala, como lo evidencia su ausencia en los estudios de referencia recientes de SPE-11.

Metodología
Los autores presentan un nuevo esquema de volumen finito, denominado PPU-C, diseñado para imponer el equilibrio capilar local en los límites de tipos de roca dentro del marco de simulación de código abierto GEOS. La metodología se basa en el método PPU estándar, pero introduce condiciones de interfaz discretas para manejar interfaces heterogéneas.

Componentes metodológicos clave:

Variable Primaria de Interfaz: En lugar de utilizar saturaciones como la variable primaria (lo que a menudo requiere un cambio de variable complejo entre estados de una sola fase y de dos fases), el PPU-C utiliza una presión capilar de interfaz ( $P_{c,ij}$ ) como la incógnita primaria en cada interfaz de material.
Solucionador No Lineal Local: En cada interfaz heterogénea, se emplea un solucionador no lineal local para imponer dos condiciones:
1. Continuidad de Flujo: Conservación de la tasa de flujo de masa para cada fase a través de la interfaz.
2. Continuidad de Presión Capilar Extendida: Imposición de la continuidad de las presiones de fase (y, por ende, de la presión capilar) solo cuando ambas fases están fluyendo. Si una fase está atrapada, la presión de la fase atrapada puede ser discontinua. Esto se gestiona mediante una formulación de "presión capilar extendida" que transiciona suavemente entre los estados de atrapamiento y de flujo.
Integración de Histéresis: El esquema se extiende para incorporar la histéresis de la presión capilar y la permeabilidad relativa utilizando el modelo de Killough. Esto permite el modelado simultáneo del atrapamiento residual (impulsado por el cierre de poros o snap-off y la historia de reversión de flujo) y el atrapamiento por heterogeneidad capilar. El solucionador local evita el cambio de variables al resolver directamente la presión de la interfaz, incluso cuando existen curvas de escaneo.
Implementación: El método se implementa dentro del marco de Newton-Krylov totalmente implícito de GEOS. El problema escalar no lineal local en las interfaces se resuelve mediante un método de Newton salvaguardado con bisección y amortiguamiento, el cual es computacionalmente económico y paralelizable.

Contribuciones Clave

Mejora de la Precisión en Mallas Gruesas: El estudio demuestra que el PPU-C mejora significativamente la precisión de los volúmenes predichos de atrapamiento capilar y residual en mallas prácticas en comparación con el PPU estándar, sin requerir el refinamiento excesivo de la malla que típicamente se necesita para resolver la física de la interfaz.
Manejo Eficiente de la Histéresis: Los autores extienden las formulaciones previas mejorando la eficiencia del solucionador local e integrando con éxito modelos de histéresis sin necesidad de cambio de variables, lo que simplifica la implementación y la robustez del esquema.
Escalabilidad y Desempeño: El esquema se implementa en un marco de computación de alto rendimiento (HPC) llamado GEOS y se valida en configuraciones multidimensionales, incluyendo una sección transversal 2D a escala de campo basada en el benchmark SPE-11B.

Resultados
El método fue evaluado mediante cuatro ejemplos numéricos de complejidad creciente:

Segregación por Gravedad 1D: El PPU-C convergió a la solución analítica de estado estacionario mucho más rápido que el PPU. En mallas gruesas (por ejemplo, 10 celdas), el PPU estándar no logró predecir ningún atrapamiento, mientras que el PPU-C capturó el volumen atrapado con precisión.
Segregación por Gravedad 3D: En entornos 3D, el PPU-C mantuvo su ventaja de precisión. Notablemente, en simulaciones 3D de alta resolución donde las soluciones de PPU y PPU-C coinciden, el PPU-C demostró una convergencia no lineal superior, resultando en aceleraciones de hasta 2.8× en comparación con el PPU estándar debido a la reducción de cortes en los pasos de tiempo e iteraciones descartadas.
Dominio 2D con Múltiples Barreras: En un dominio con barreras capilares aleatorias, el PPU-C produjo distribuciones de saturación finales diferentes al PPU, mostrando un mayor atrapamiento en las capas inferiores y un retraso en la invasión de las capas superiores.
Benchmark SPE-11B (Escala de Campo): En una sección transversal 2D del benchmark SPE-11B (8.4 km × 1.2 km), el PPU-C predijo rutas de migración de la pluma distintas, con el PPU invadiendo las capas superiores más rápido que el PPU-C.

Sobrecarga Computacional
Se determinó que el costo computacional de imponer las condiciones de interfaz fue moderado. En todos los casos de prueba, la sobrecarga del PPU-C en comparación con el PPU no excedió el 25%. En varios casos 3D, la mejora en la convergencia no lineal del PPU-C resultó en una aceleración neta, haciendo que el método fuera más eficiente que el enfoque estándar en ciertos escenarios de alta resolución.

Significancia y Reivindicaciones
El artículo concluye que el uso de condiciones de interfaz discretas es una técnica numérica práctica y eficiente para modelar con precisión la heterogeneidad capilar en simulaciones de almacenamiento de CO2 a gran escala. Los autores afirman que:

Los modelos PPU estándar subestiman consistentemente las cantidades de CO2 atrapadas de forma segura en mallas prácticas, a menos que se utilice un refinamiento excesivo.
El PPU-C reduce los errores de heterogeneidad capilar en un orden de magnitud en comparación con el PPU estándar en las mismas mallas.
El método maneja eficazmente la combinación de la física de la heterogeneidad capilar y la histéresis, lo cual es crítico para una evaluación precisa del GCS.
La modesta sobrecarga computacional (≤25%) y el potencial de aceleración en escenarios 3D complejos eliminan la barrera para la adopción práctica de las condiciones de interfaz en las simulaciones industriales de GCS.

Los autores señalan que, aunque los resultados actuales se centran en flujos inmiscibles, la formulación se extiende a simulaciones composicionales, y el trabajo futuro abordará el atrapamiento por disolución.

Más como este