Mathematical Modeling of Salt Precipitation and Multi-Phase Flow in High Enthalpy Fractured Geothermal Systems

Este artículo presenta un nuevo modelo de flujo composicional de código abierto implementado en el marco de trabajo PorePy que simula el flujo multifásico no isotérmico y la precipitación de halita en reservorios geotérmicos fracturados de alta entalpía, utilizando una formulación de variable primaria robusta y un enfoque de matriz de fractura discreta para predecir con precisión el daño por permeabilidad y los desafíos operativos.

Lo esencial

Imagina un reservorio subterráneo de alta temperatura como una olla a presión gigante y natural llena de agua salada supercaliente. Esta es la fuente de la energía geotérmica. Para extraer esa energía, los ingenieros perforan pozos y bombean agua fría hacia adentro para empujar el agua caliente de regreso hacia arriba. Sin embargo, este proceso es complicado debido a la "sal" en el agua.

Piensa en la sal (halita) como el azúcar en una taza de té caliente. Si enfrías el té o dejas que parte del agua se evapore, el azúcar no puede permanecer disuelto y comienza a convertirse nuevamente en cristales sólidos. En un pozo geotérmico, esto sucede cuando el agua caliente se enfría cerca de un pozo de inyección o se evapora cerca de un pozo de producción. ¿El resultado? Se forman cristales de sal sólida que obstruyen los diminutos agujeros en las rocas y las grietas (fracturas) que permiten el flujo del agua. Es como un embotellamiento causado por cristales de azúcar bloqueando una autopista.

Este artículo presenta una nueva herramienta de simulación por computadora diseñada para predecir exactamente dónde y cómo ocurrirán estos "atascos de azúcar" en rocas subterráneas complejas y agrietadas.

Aquí hay un desglose de cómo funciona la herramienta y lo que encontró, utilizando analogías simples:

1. El Mapa: Viendo las Grietas con Claridad

Los mapas tradicionales de rocas subterráneas suelen suavizar todo, tratando la roca como una esponja sólida. Pero en la realidad, el agua fluye principalmente a través de una red de grietas, como el agua que fluye a través de una acera agrietada en lugar de a través del concreto mismo.

• La Innovación: Este nuevo modelo utiliza un enfoque de "Matriz de Fractura Discreta". Imagina dibujar las grietas como líneas finas y distintas en un mapa, en lugar de simplemente difuminarlas en el fondo. Esto permite que la computadora vea exactamente cómo se conectan las grizas (o no se conectan) y cómo la sal podría obstruir una grieta específica frente a la roca circundante.

2. El Motor: Un "Control Remoto Universal" para la Física

Simular agua hirviendo, vapor y sal sólida al mismo tiempo es increíblemente difícil para las computadoras. Usualmente, la computadora tiene que cambiar constantemente su "modo" (por ejemplo, "Bien, ahora es líquido; ahora es gas; ahora es sólido"), lo que puede causar que el cálculo falle o se detenga.

• La Innovación: Los autores crearon un sistema "unificado". Piensa en esto como un control remoto universal que funciona para cada dispositivo sin necesidad de cambiar baterías o modos. El modelo utiliza tres "perillas" fijas (Presión, Energía Térmica y Cantidad de Sal) que permanecen constantes ya sea que el agua sea líquida, vapor o se esté convirtiendo en sal sólida. Esto hace que la simulación sea mucho más fluida y estable, permitiéndole manejar el caótico cambio entre estados sin romperse.

3. El Truco de Velocidad: La "Hoja de Trucos"

Calcular la física exacta del agua salada a altas temperaturas suele requerir que la computadora resuelva acertijos matemáticos complejos una y otra vez, lo cual es muy lento.

• La Innovación: El equipo creó una "hoja de trucos" precalculada (tabla de búsqueda). Antes de que comience la simulación, calcularon todos los resultados posibles sobre cómo se comporta la sal bajo diferentes condiciones y los almacenaron. Durante la simulación, en lugar de resolver la matemática difícil cada vez, la computadora simplemente busca la respuesta en la hoja. Esto hace que la simulación corra mucho más rápido manteniendo la precisión.

4. El Efecto de Obstrucción: El "Encogimiento del Tamaño de los Poros"

Cuando se forman cristales de sal, estos ocupan espacio.

• La Innovación: El modelo reduce automáticamente los "tubos" (porosidad y permeabilidad) a medida que la sal se acumula. Utiliza una regla (Kozeny-Carman) que dice: "Si la sal llena el 10% del agujero, el tubo se vuelve significativamente más estrecho". Esto permite que el modelo prediga cómo el flujo se ralentizará o se detendrá por completo a medida que el "atasco de azúcar" empeore.

Lo que las Simulaciones Mostraron

El equipo probó esta herramienta en dos escenarios principales:

Escenario A: La Carretera Rota (Grietas Desconectadas)

• Configuración: Imagina un reservorio donde las grietas no se conectan entre sí; el agua tiene que apretarse a través de la roca sólida entre ellas.
• Resultado: Cuando bombaron agua fría, el agua caliente cerca del pozo de producción hirvió rápidamente. Esto causó que la sal se cristalizara y obstruyera la roca justo alrededor del pozo.
• El Giro: Si bombeaban el agua más rápido, la obstrucción empeoraba considerablemente y la producción de energía caía significamente. El modelo mostró que el "embotellamiento" ocurrió principalmente en la roca cerca del pozo, no solo en las grietas.

Escenario B: La Carretera Conectada (Grietas Conectadas)

• Configuración: Imagina un reservorio donde las grietas forman una autopista continua de alta velocidad desde el pozo de inyección hasta el pozo de producción.
• Resultado: El agua fría pasó rápidamente a través de las grietas. Debido a que se movía tan rápido y se mantenía fresca, ¡en realidad disolvió la sal cerca del pozo de producción en lugar de obstruirlo!
• El Giro: La precipitación de sal se movió a un lugar diferente —justo en el borde de la zona de agua fría— en lugar de obstruir el pozo mismo. Esto sugiere que tener una red de grietas conectadas podría, de hecho, proteger al pozo de la obstrucción, aunque cambie el lugar donde se acumula la sal.

La Conclusión Final

Este artículo presenta una nueva herramienta de software de código abierto que ayuda a los ingenieros a comprender la compleja danza entre el calor, la presión y la sal en los pozos geotérmicos. Al mapear con precisión cómo se conectan las grietas y cómo la sal las obstruye, la herramienta puede ayudar a predecir:

1. Dónde podrían bloquearse los pozos por la sal.
2. Cuánta energía se puede extraer de forma segura antes de que los "tubos" se obstruyan.
3. Si la disposición de las grietas subterráneas ayudará o perjudicará el proceso de producción.

Los autores verificaron su herramienta contra un estándar establecido de la industria y encontraron que coincidía perfectamente, demostando que es una forma confiable de simular estos entornos de roca agrietada y salina de alta temperatura.

Resumen técnico

Resumen Técnico: Modelo de Flujo Composicional Unificado para la Simulación de Sistemas de Alta Entalpía Multifásicos y Fracturados

Planteamiento del Problema
La simulación de yacimientos geotérmicos fracturados de alta entalpía presenta desafíos significativos debido al complejo acoplamiento de flujo multifásico y multicomponente no isotérmico, la termodinámica fuertemente no lineal y la influencia dominante de las redes de fracturas. Un problema operativo crítico en estos sistemas es la incrustación mineral, específicamente la precipitación de halita (NaCl). Durante la producción, el descenso de presión induce la ebullición, concentrando la salmuera residual y excediendo los límites de solubilidad, mientras que la inyección de agua fría reduce la solubilidad local. Ambos mecanismos conducen a la precipitación de sal dentro de los poros y las fracturas, deteriorando severamente la porosidad y la permeabilidad, reduciendo la productividad del pozo y, potencialmente, causando el bloqueo del pozo. Los modelos de continuo existentes a menudo no logran capturar la precipitación localizada y la evolución de la permeabilidad a escala de fractura, mientras que los simuladores composicionales tradicionales que dependen de formulaciones de cambio de variables presentan dificultades de estabilidad numérica cerca de los límites de fase.

Metodología
Los autores presentan un nuevo modelo de flujo composicional unificado implementado dentro del marco de código abierto PorePy. La metodología integra tres componentes principales:

1. Formulación de Variables Persistentes: A diferencia de los enfoques tradicionales que cambian las variables primarias según la presencia de fases, este modelo utiliza un conjunto fijo de variables primarias —presión ($p$), entalpía específica ($h$) y fracción de masa de sal total ($z$)— a través de todas las configuraciones de fase (líquido, vapor y halita sólida). Este enfoque gestiona naturalmente las transiciones de fase sin necesidad de cambios manuales, mejorando la robustez y estabilidad numérica.
2. Enfoque de Matriz y Fractura Discreta (DFM): El dominio físico se representa mediante una geometría de dimensiones mixtas donde las fracturas y sus intersecciones se tratan como variedades de menor dimensión (1D y 0D) embebidas en la matriz porosa 2D. Esta representación explícita preserva la conectividad y heterogeneidad de las fracturas, permitiendo resolver procesos localizados como la reducción de permeabilidad inducida por la precipitación.
3. Cierre Termodinámico y Linealización: El modelo emplea las correlaciones exhaustivas de Driesner y Christoph (2007) para la termodinámica de H2O–NaCl, cubriendo condiciones magmáticas. Para reducir el costo computacional, el modelo utiliza una estrategia de Linealización Basada en Operadores (OBL). Las propiedades termodinámicas y sus derivadas se precalculan en una rejilla estructurada en el espacio de variables primarias y se recuperan mediante interpolación multilineal durante la simulación, eliminando los costosos cálculos de separación de fases (flash) en tiempo real.
4. Evolución Dinámica de la Permeabilidad y la Apertura: El modelo incorpora relaciones de tipo Kozeny-Carman para actualizar dinámicamente la porosidad y la permeabilidad de la matriz basándose en la saturación de halita. Para las fracturas, la apertura evoluciona según una relación de ley de potencia dependiente de la saturación de halita, lo que a su vez actualiza la permeabilidad de la fractura mediante la ley cúbica.

Las ecuaciones gobernantes (balance de masa y energía) se discretizan utilizando un esquema de volúmenes finitos de centro de celda con Aproximación de Flujo de Múltiples Puntos (MPFA) para términos difusivos y upwinding para términos advectivos. El sistema no lineal resultante se resuelve mediante el método de Newton con una búsqueda de línea Armijo de retroceso y pasos de tiempo adaptativos.

Resultos Clave
El modelo fue verificado y aplicado a través de varios experimentos numéricos:

• Verificación (Benchmark 1D): El modelo fue verificado contra el simulador de código cerrado CSMP++ utilizando un benchmark de disolución de sal en 1D. Los resultados mostraron una fuerte concordancia a través de cuatro regiones de fase distintas (vapor+halita, vapor+líquido, fase líquida simple y líquido+halita) y sus transiciones, validando la formulación de variables persistentes y la implementación termodinámica.
• Red de Fracturas Desconectada (2D): En un yacimiento con fracturas aisladas, el flujo está dominado por la matriz. La simulación demostró que el descenso de presión en el pozo de producción induce la ebullición, provocando una precipitación significativa de halita en la matriz circundante y en los segmentos de fractura adyacentes. Esta precipitación redujo la permeabilidad de la matriz hasta en un 36% ($K/K_0 \approx 0.64$) y alteró las líneas de corriente de flujo. El aumento de la tasa de inyección intensificó la ebullición cerca del pozo, lo que llevó a una mayor saturación de halita ($s_{hal} \approx 0.46$) y una reducción de permeabilidad más severa ($K/K_0 \approx 0.29$), impactando significativamente las tasas de producción.
• Red de Fracturas Conectada (2D): En un yacimiento con una cadena de fracturas continua que conecta los pozos de inyección y producción, el transporte pasó a estar dominado por las fracturas. La red conectada facilitó la comunicación rápida de presión y el transporte de fluido de baja salinidad hacia el pozo de producción. En consecuencia, aunque ocurrió ebullición cerca del productor, el fluido que llegaba permaneció suficientemente diluido para evitar la saturación de halita local. En lugar de la precipitación en el productor, se observó disolución, y la precipitación se redistribió hacia la matriz delante del frente de inyección y a lo largo de la cadena de fracturas debido a los efectos de enfriamiento. Esta configuración evitó la severa degradación de la permeabilidad observada en el caso desconectado.

Significancia y Reivindicaciones
El artículo afirma proporcionar una herramienta numérica robusta y de código abierto para analizar el complejo transporte de calor y masa con efectos de incrustación mineral en sistemas geotérmicos de fractura de alta entalpía. Su principal significancia radica en:

• Estabilidad Numérica: La formulación unificada de variables persistentes evita las dificultades de convergencia asociadas con el cambio de variables cerca de los límites de fase.
• Fidelidad Geométrica: El enfoque DFM resuelve explícitamente la conectividad de las fracturas, lo cual los autores demuestran es un factor crítico que controla la distribución espacial de la precipitación y la consecuente reducción de transmisibilidad.
• Perspectiva Operativa: Los resultados resaltan que la conectividad de las fracturas dicta los desafíos operativos; las redes desconectadas son propensas al colmatado cerca del pozo y a la pérdida de productividad, mientras que las redes conectadas pueden mitigar estos problemas manteniendo la dilución del fluido en el pozo de producción.
• Reproducibilidad: El trabajo proporciona el modelo numérico completo y el flujo de ejecución mediante un contenedor Docker, facilitando la reproducibilidad y el desarrollo posterior dentro de la comunidad geotérmica.

Los autores concluyen que su modelo predice eficazmente los patrones de precipitación de halita y su impacto en el rendimiento del yacimiento, ofreciendo una herramienta valiosa para comprender la interacción entre el flujo, la termodinámica y la incrustación mineral en yacimientos geotérmicos fracturados.