Heat Transfer Modeling in Enhanced Geothermal Energy: A Three-Temperature Approach for Solid, Injected, and Residing Fluids

Aquí presento un resumen técnico detallado del artículo en español, estructurado según los puntos solicitados:

Resumen Técnico: Modelado de Transferencia de Calor en Energía Geotérmica Mejorada (EGS)

Título del trabajo: Heat Transfer Modeling in Enhanced Geothermal Energy: A Three-Temperature Approach for Solid, Injected, and Residing Fluids (Modelado de transferencia de calor en energía geotérmica mejorada: Un enfoque de tres temperaturas para sólidos, fluidos inyectados y residentes).

1. Planteamiento del Problema

Los sistemas geotérmicos mejorados (EGS) extraen calor de formaciones rocosas profundas de baja permeabilidad inyectando un fluido de trabajo frío. La eficiencia y viabilidad a largo plazo dependen de procesos de flujo y transferencia de calor fuertemente acoplados en medios porosos fracturados, a menudo dominados por la advección.

El problema central identificado es la limitación de los modelos convencionales:

• Equilibrio Térmico Local (LTE): Asume una única temperatura efectiva para el sistema roca-fluido, ignorando las diferencias térmicas entre fases.
• No Equilibrio Térmico Local (LTNE) tradicional: Introduce temperaturas separadas para la fase sólida y la fase fluida, pero trata al fluido como un componente promedio único.

La limitación crítica: En ambos casos, la evolución térmica del fluido inyectado frío no se resuelve explícitamente. Se infiere indirectamente a través de promedios, lo que oculta dinámicas detalladas como el calentamiento de los paquetes de fluido recién inyectados, la canalización en fracturas y la inestabilidad de dedos térmicos. Esto dificulta la predicción precisa de la "ruptura térmica" (cuándo el fluido frío llega a los pozos de producción) y la vida útil del reservorio.

2. Metodología Propuesta

Los autores desarrollan un nuevo marco de modelado y una estrategia numérica robusta:

A. Modelo Matemático: LTNE de Tres Vías (Three-way LTNE)
Se propone un sistema que distingue explícitamente tres campos de temperatura:

1. $T_i$ (Fluido Inyectado): Temperatura del fluido frío que entra al sistema.
2. $T_r$ (Fluido Residente): Temperatura del fluido que ya ocupa los poros antes de la inyección.
3. $T_s$ (Matriz Sólida): Temperatura de la roca.

• Variable de Concentración ($z$): Se introduce una variable de concentración ($0 \le z \le 1$) que rastrea la fracción volumétrica del fluido inyectado. Esta variable acopla las ecuaciones de energía.
• Ecuaciones de Energía: Se derivan ecuaciones de conservación de energía separadas para $T_i$, $T_r$ y $T_s$. Estas incluyen:
  • Términos de transporte por advección (dependientes de la velocidad Darcy y la concentración).
  • Difusión térmica efectiva.
  • Intercambio de calor interfacial: Sólido-Fluido Inyectado, Sólido-Fluido Residente y, crucialmente, Fluido Inyectado-Fluido Residente.
• Consistencia Termodinámica: El modelo demuestra que, en el límite de intercambio de calor interfacial infinito, se reduce al modelo LTE clásico, y en el límite de dos temperaturas, al modelo LTNE estándar.

B. Discretización Numérica
Para resolver el sistema no lineal y dominado por la advección, se emplea una combinación de métodos avanzados:

• Método de Galerkin Enriquecido (EG): Se utiliza para la ecuación de Darcy, temperatura y concentración. Combina elementos continuos $Q1$ con una enriquecimiento discontinuo $Q0$.
  • Ventaja: Garantiza la conservación local de masa en cada elemento con un número moderado de grados de libertad, esencial para simulaciones a gran escala.
• Estrategia de Transporte Corregido por Flujo (FCT): Se implementa una estrategia FCT adaptada para la discretización EG de las ecuaciones de concentración y temperatura.
  • Objetivo: Imponer un principio de máximo discreto y suprimir oscilaciones no físicas (sobreshoots/undershoots) sin sacrificar la conservación local.
• Esquema de Integración Temporal:
  • Se utiliza un esquema de tipo IMPES (Presión Implícita, Concentración Explícita, Calor Implícito) para desacoplar el sistema no lineal.
  • Presión/Velocidad: Implícito.
  • Concentración: Explícito con un esquema Runge-Kutta de segundo orden que preserva la estabilidad fuerte (SSP-RK2).
  • Temperaturas: Implícito con el esquema BDF2 (diferenciación hacia atrás de segundo orden) para manejar la rigidez de los términos difusivos.

3. Contribuciones Clave

1. Resolución Explícita del Fluido Inyectado: Por primera vez en este contexto, el modelo permite rastrear cómo los paquetes de fluido inyectado se calientan a medida que viajan, diferenciándolos del fluido residente, algo imposible con modelos de temperatura promedio.
2. Marco LTNE de Tres Temperaturas: Formalización matemática rigurosa que incluye el intercambio de calor directo entre el fluido inyectado y el residente, escalado por la concentración y la porosidad.
3. Algoritmo Numérico Estable y Conservativo: Desarrollo de un esquema EG-FCT que logra alta precisión, conservación local y estabilidad en regímenes dominados por la advección, evitando oscilaciones espurias en frentes de concentración abruptos.
4. Validación de Límites: Demostración teórica y numérica de cómo el modelo se reduce correctamente a LTE y LTNE de dos temperaturas bajo condiciones específicas.

4. Resultados Numéricos

Los autores presentan tres experimentos numéricos utilizando la biblioteca deal.II:

• Prueba de Convergencia: Se verificó la precisión del acoplamiento entre la ecuación de concentración y el modelo LTNE de tres vías. Los resultados mostraron tasas de convergencia óptimas (cercanas a segundo orden en norma $L^2$ y primer orden en $H^1$), confirmando la validez del esquema EG-FCT.
• Comparación LTE vs. LTNE:
  • En regímenes de intercambio de calor rápido, el modelo LTNE converge rápidamente al LTE, validando la teoría asintótica.
  • En regímenes de intercambio lento (típicos de EGS), se observaron diferencias significativas entre $T_f$ y $T_s$. El modelo LTE falló en predecir correctamente la temperatura de equilibrio mixto, especialmente cuando hay advección, demostrando la necesidad del enfoque LTNE.
• Caso de Estudio EGS con Fracturas:
  • Se simuló un medio heterogéneo con 34 canales de fractura de alta permeabilidad.
  • El modelo de tres vías reveló estructuras térmicas ocultas: mostró cómo el fluido inyectado ($T_i$) se calienta rápidamente al pasar por las fracturas, mientras que el fluido residente ($T_r$) se enfría gradualmente.
  • Se observó un "salto" de temperatura en los frentes de concentración que los modelos promediados no capturan.
  • La simulación permitió visualizar la evolución de la entalpía y la ruptura térmica con mayor detalle espacial.

5. Significado e Impacto

Este trabajo representa un avance significativo en la simulación de reservorios geotérmicos:

• Precisión Predictiva: Proporciona una herramienta más informativa para el diseño y optimización de EGS, permitiendo predecir con mayor exactitud cuándo el fluido frío alcanzará los pozos de producción (ruptura térmica), lo cual es crítico para la gestión de la vida útil del reservorio.
• Comprensión de Mecanismos: Permite entender fenómenos multiescala como la canalización en fracturas y la inestabilidad de dedos térmicos, que son invisibles para los modelos de temperatura única.
• Viabilidad Numérica: Demuestra que es posible resolver sistemas complejos de tres temperaturas con conservación local y estabilidad, haciendo viable su aplicación en escenarios reales de ingeniería geotérmica.

En conclusión, el enfoque propuesto supera las limitaciones de los modelos tradicionales al descomponer la temperatura del fluido, ofreciendo una visión más realista y detallada de la dinámica térmica en sistemas geotérmicos mejorados.