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

Titel en Context

Titel: Warmteoverdrachtmodellering in Verbeterte Geothermische Systemen (EGS): Een Drie-Temperatuurbenadering voor Vaste Stof, Ingepompde en Aanwezige Vloeistoffen.
Auteurs: Yi-Yung Yang, Sanghyun Lee en Dmitri Kuzmin.
Doel: Het ontwikkelen van een geavanceerd numeriek model om warmteoverdracht in verbroken porieuze media (zoals die in EGS) nauwkeuriger te simuleren dan traditionele modellen.

---

1. Het Probleem

Verbeterte Geothermische Systemen (EGS) halen warmte uit diepe, laag-permeabele gesteentevormingen door koud werkende vloeistof in te pompen, deze door een verbroken porieuze medium te laten circuleren en het te produceren op verhoogde temperaturen.

• Beperkingen van bestaande modellen:
  • Lokale Thermische Evenwicht (LTE): De meeste simulators gaan ervan uit dat de vaste stof en de vloeistof dezelfde temperatuur hebben ($T_{vaste} = T_{vloeistof}$). Dit negeert de tijdsvertraging in warmteuitwisseling.
  • Traditionele LTNE (Local Thermal Non-Equilibrium): Zelfs geavanceerdere LTNE-modellen behandelen vaak de porievloeistof als één gemiddelde component. Hierdoor wordt de thermische evolutie van de ingepompte koude vloeistof niet expliciet opgelost, maar slechts indirect afgeleid via gemiddelden.
• Gevolg: Dit gemiddelde verbergt de gedetailleerde verwarmingsdynamiek van nieuwe vloeistofpakketten en beperkt het vermogen om multischaal-fenomenen zoals kanaalvorming in breuken, vingerinstabiliteiten en lokale warmteuitwisselingspunten op te lossen.

---

2. Methodologie

De auteurs ontwikkelen een nieuw drie-weg LTNE-framework (Three-way LTNE) dat onderscheid maakt tussen drie temperatuurvelden:

1. $T_i$: Temperatuur van de ingepompte vloeistof (injected fluid).
2. $T_r$: Temperatuur van de aanwezige vloeistof (resident fluid).
3. $T_s$: Temperatuur van het vaste gesteentematrix (solid matrix).

Kerncomponenten van het model:

• Concentratievariabele ($z$): Een variabele die de fractie van de ingepompte vloeistof in het porievolume trackt ($0 \le z \le 1$). Deze drijft de koppeling aan.
• Enthalpie-decompositie: De totale vloeistofenthalpie wordt opgesplitst: $C_f T_f = C_i T_i + C_r T_r$.
• Gekoppelde vergelijkingen: Het systeem omvat Darcy-stroming en aparte energie-vergelijkingen voor de drie fasen, inclusief termen voor warmteuitwisseling tussen:
  • Ingepompte vloeistof en vaste stof.
  • Aanwezige vloeistof en vaste stof.
  • Ingepompte en aanwezige vloeistof onderling (directe uitwisseling).
• Asymptotische analyse: Het artikel toont wiskundig aan dat dit model reduceert tot het klassieke LTE-model wanneer de warmteoverdrachtscoëfficiënten zeer groot zijn.

Numerieke Discretisatie:
Om het sterk gekoppelde, advectie-gedomineerde systeem op te lossen, gebruiken de auteurs een geavanceerde numerieke aanpak:

1. Verrijkte Galerkin (EG) Methode: Een eindige-elementenmethode die continue $Q_1$-elementen combineert met discontinu $Q_0$-verrijking. Dit garandeert lokale massabehoud (cruciaal voor transport) met een redelijk aantal vrijheidsgraden.
2. Flux-Corrected Transport (FCT): Een strategie om de transportvergelijkingen voor concentratie en temperatuur te stabiliseren. Dit dwingt het discrete maximumprincipe af en onderdrukt niet-fysische oscillaties (overshoots/undershoots) terwijl lokale behoudswetten worden gehandhaafd.
3. Tijdsintegratie:
   • IMPES-achtige splitting: Druk en snelheid worden impliciet opgelost, concentratie expliciet, en temperaturen impliciet.
   • SSP-RK2: Voor de concentratievergelijking wordt een Strong-Stability-Preserving Runge-Kutta methode van orde 2 gebruikt.
   • BDF2: Voor de temperatuurvergelijkingen wordt een impliciete tweede-orde methode gebruikt voor stabiliteit.

---

3. Belangrijkste Resultaten

De auteurs presenteren drie numerieke experimenten:

1. Convergentietest:

   • Het algoritme toont de verwachte optimale convergentie: tweede orde in de $L^2$-norm en eerste orde in de $H^1$-norm voor zowel concentratie als temperaturen.
   • De FCT-limiter blijft inactief bij gladde oplossingen, waardoor de hoge orde nauwkeurigheid behouden blijft.

2. Vergelijking LTE vs. LTNE:

   • Snelle warmteuitwisseling: Wanneer de warmteoverdrachtscoëfficiënt hoog is, gedragen LTE en LTNE zich vergelijkbaar (temperatuurverschillen zijn verwaarloosbaar).
   • Trage warmteuitwisseling: Bij lage coëfficiënten (en aanwezigheid van advectie) wijken de temperaturen van vaste stof en vloeistof sterk van elkaar af. Het klassieke LTE-model faalt hier om de werkelijke evenwichtstemperatuur of de temperatuurverschillen correct weer te geven.

3. Toepassing op Verbroken Media (Three-way LTNE):

   • In een scenario met 34 hoge-permeabiliteit breuken in een laag-permeabele matrix toont het model aan dat de ingepompte vloeistof ($T_i$) en de aanwezige vloeistof ($T_r$) verschillende thermische paden volgen.
   • Het model kan de verwarming van de ingepomste vloeistof langs stromingspaden expliciet volgen.
   • Er worden scherpe temperatuursprongen waargenomen bij concentratiefronten die door gemiddelde modellen worden gemist.
   • De enthalpie van de ingepomste vloeistof ($z \cdot T_i$) toont aan dat de warmteoverdracht sterk gekoppeld is aan de concentratiefronten.

---

4. Bijdrage en Significantie

• Nieuwe Modellering: Het introduceert een uniek "drie-temperatuur" framework dat de fysica van EGS beter weergeeft door de ingepomste vloeistof als een aparte thermische entiteit te behandelen.
• Verbeterde Voorspelling: Het stelt onderzoekers in staat om thermische doorbraak (thermal breakthrough) en de levensduur van reservoirs nauwkeuriger te voorspellen, omdat het de verwarmingsgeschiedenis van individuele vloeistofpakketten kan traceren.
• Robuuste Numeriek: De combinatie van EG, FCT en IMPES biedt een stabiele, lokaal behoudende oplossing voor complexe, advectie-gedomineerde problemen in de geothermie.
• Praktische Toepassing: Het model biedt een krachtig hulpmiddel voor het ontwerp en de optimalisatie van geothermische reservoirs, vooral in heterogene, verbroken formaties waar traditionele gemiddelde modellen tekortschieten.

Conclusie:
Dit werk biedt een significante stap voorwaarts in de modellering van geothermische systemen. Door de thermische scheiding tussen ingepomste en aanwezige vloeistof expliciet te modelleren, onthult het thermische structuren die anders onzichtbaar blijven, wat essentieel is voor het maximaliseren van de efficiëntie en veiligheid van EGS-projecten.