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

Dit artikel presenteert een nieuw open-source compositioneel stromingsmodel geïmplementeerd in het PorePy-framework dat niet-isotherme, multifasische stroming en halietprecipitatie in hoog-enthalpie gebarsten geothermische reservoirs simuleert, gebruikmakend van een robuuste primaire variabele formulering en een discrete breuk-matrix benadering om permeabiliteitsschade en operationele uitdagingen nauwkeurig te voorspellen.

De kern

Stel je een heet ondergronds reservoir voor als een gigantische, natuurlijke hogedrukpan gevuld met superheet zout water. Dit is de bron van geothermische energie. Om die energie eruit te krijgen, boren ingenieurs putten en pompen ze koud water naar binnen om het hete water weer omhoog te duwen. Echter, dit proces is lastig vanwege het "zout" in het water.

Denk aan het zout (haliet) als suiker in een kop warme thee. Als je de thee afkoelt of wat water laat verdampen, kan de suiker niet meer opgelost blijven en begint het weer terug te veranderen in vaste kristallen. In een geothermische put gebeurt dit wanneer het water afkoelt nabij een injectieput of wegkookt nabij een productieput. Het resultaat? Vaste zoutkristallen vormen zich en verstoppen de kleine gaatjes in de rotsen en de barsten (fracturen) waar het water doorheen stroomt. Het is als een verkeersopstopping veroorzaakt door suikerkristallen die een snelweg blokkeren.

Dit artikel presenteert een nieuwe computersimulatietool die ontworpen is om precies te voorspellen waar en hoe deze "suikeropstoppingen" zullen plaatsvinden in complexe, gebarsten ondergrondse gesteenten.

Hier is een overzicht van hoe de tool werkt en wat het heeft gevonden, met behulp van eenvoudige analogieën:

1. De Kaart: De Barsten Helder Zien

Traditionele kaarten van ondergrondse gesteenten maken alles vaak glad, waarbij ze het gesteente behandelen als een massieve spons. Maar in werkelijkheid stroomt het water vooral door een netwerk van barsten, zoals water door een gebarsten stoep stroomt in plaats van door het beton zelf.

• De Innovatie: Dit nieuwe model gebruikt een "Discrete Fracture-Matrix"-benadering. Stel je voor dat je de barsten als duidelijke, dunne lijnen op een kaart tekent, in plaats van ze simpelweg te vervagen in de achtergrond. Hierdoor kan de computer precies zien hoe de barsten met elkaar verbonden zijn (of niet) en hoe zout een specifieke barst kan verstoppen versus de omliggende rots.

2. De Motor: Een "Universele Afstandsbediening" voor Fysica

Het simuleren van kokend water, stoom en vast zout tegelijkertijd is extreem moeilijk voor computers. Meestal moet de computer constant van "modus" wisselen (bijv. "Oké, nu is het vloeistof; nu is het gas; nu is het vast"), wat kan leiden tot het crashen of vastlopen van de berekening.

• De Innovatie: De auteurs hebben een "verenigd" systeem gecreëerd. Denk aan een universele afstandsbediening die voor elk apparaat werkt zonder dat je van batterij of modus hoeft te wisselen. Het model gebruikt drie vaste "knoppen" (Druk, Warmte-energie en Zouthoeveelheid) die hetzelfde blijven of het water nu vloeibaar, stoom of vast zout is. Dit maakt de simulatie veel soepeler en stabieler, waardoor het de chaotische overgangen tussen toestanden kan afhandelen zonder vast te lopen.

3. De Snelheidstrik: Het "Spiekbriefje"

Het berekenen van de exacte fysica van zout water bij hoge temperaturen vereist meestal dat de computer keer op keer complexe wiskundige puzzels oplost, wat erg traag is.

• De Innovatie: Het team heeft een vooraf berekend "spiekbriefje" (zoeketabel) gemaakt. Voordat de simulatie begint, hebben ze alle mogelijke uitkomsten berekend voor hoe het zout zich onder verschillende omstandigheden gedraagt en deze opgeslagen. Tijdens de simulatie hoeft de computer in plaats van telkens de moeilijke wiskunde op te lossen, alleen maar het antwoord op het blaadje op te zoeken. Dit maakt de simulatie veel sneller terwijl de nauwkeurigheid behouden blijft.

4. Het Verstoppingseffect: De "Porositeit-Krimp"

Wanneer zoutkristallen ontstaan, nemen ze ruimte in beslag.

• De Innovatie: Het model laat de "buizen" (porositeit en permeabiliteit) automatisch krimpen naarmate het zout zich opbouwt. Het gebruikt een regel (Kozeny-Carman) die zegt: "Als zout 10% van het gat vult, wordt de buis aanzienlijk smaller." Hierdoor kan het model voorspellen hoe de doorstroming vertraagt of zelfs volledig stopt naarms de "suikeropstopping" erger wordt.

Wat de Simulatiesen toonden

Het team testte deze tool in twee hoofdscenario's:

Scenario A: De Gebroken Snelweg (Niet-verbonden Barsten)

• Opzet: Stel je een reservoir voor waar de barsten niet met elkaar verbonden zijn; het water moet door het massieve gesteente tussen de barsten door persen.
• Resultaat: Wanneer ze koud water naar binnen pompten, kookte het hete water nabij de productieput zeer snel. Dit veroorzaakte dat zout kristalliseerde en de rots direct rond de put verstopte.
• De Wending: Als ze het water sneller naar binnen pompten, werd de verstopping veel erger en daalde de energieopbrengst aanzienlijk. Het model toonde aan dat de "verkeersopstopping" zich voornamelijk in de rots nabij de put bevond, en niet alleen in de barsten.

Scenario B: De Verbonden Snelweg (Verbonden Barsten)

• Opzet: Stel je een reservoir voor waar de barsten een continu, hoogwaardig snelwegnetwerk vormen van de injectieput naar de productieput.
• Resultaat: Het koude water stroomde snel door de barsten. Omdat het zo snel bewoog en koel bleef, loste het de zout nabij de productieput zelfs op in plaats van het te verstoppen!
• De Wending: De zoutprecipitatie (neerslag) verplaatste zich naar een andere plek — direct aan de rand van de koude waterzone — in plaats van de put zelf te verstoppen. Dit suggereert dat een netwerk van verbonden barsten de put juist kan beschermen tegen verstopping, ook al verandert het de plek waar het zout zich opbouwt.

De Kernboodschap

Dit artikel introduceert een nieuwe, open-source softwaretool die ingenieurs helpt de complexe dans tussen warmte, druk en zout in geothermische putten te begrijpen. Door nauwkeurig in kaart te brengen hoe barsten verbinding maken en hoe zout deze verstopt, kan de tool helpen voorspellen:

1. Waar putten mogelijk geblokkeerd raken door zout.
2. Hoeveel energie er veilig kan worden gewonnen voordat de "buizen" verstopt raken.
3. Of de lay-out van de ondergrondse barsten het productieproces zal helpen of juist zal hinderen.

De auteurs hebben hun tool geverifieerd tegenover een gevestigde industriestandaard en vonden dat deze perfect overeenkwam, wat bewijst dat het een betrouwbare manier is om deze hoogtemperatuur, zoute en gebarsten gesteentemilieus te simuleren.

Technische samenvatting

Technische Samenvatting: Een verenigd compositioneel stroommodel voor het simuleren van multifasische gefragmenteerde hoog-enthalpie systemen

Probleemstelling
Het simuleren van gefragmenteerde geothermische reservoirs met een hoge enthalpie vormt aanzienlijke uitdagingen vanwege de complexe koppeling van niet-isotherme, multifasische, multicomponentale stroming, sterk niet-lineaire thermodynamica en de dominante invloed van breuknetwerken. Een kritiek operationeel probleem in deze systemen is minerale schaling, specifnicamente haliet (NaCl) neerslag. Tijdens productie induceert drukverlaging koken, wat de resterende pekel concentreert en de oplosbaarheidslimieten overschrijdt, terwijl injectie van koud water de lokale oplosbaarheid verlaagt. Beide mechanismen leiden tot zoutneerslag in porieruimtes en breuken, wat de porositeit en permeabiliteit ernstig aantast, de putproductiviteit vermindert en potentieel verstoppingen in de put veroorzaakt. Bestaande continuümmodellen slagen er vaak niet in om gelokaliseerde neerslag en de evolutie van de permeabiliteit op breukniveau te vatten, terwijl traditionele compositionele simulators die vertrouwen op 'variable-switching' formuleringen worstelen met numerieke stabiliteit nabij fasegrenzen.

Methodologie
De auteurs presenteren een nieuw verenigd compositioneel stroommodel, geïmplementeerd binnen het open-source PorePy framework. De methodologie integreert drie kerncomponenten:

1. Persistent-Variable Formulering: In tegenstelling tot traditionele benaderingen die primaire variabelen wisselen op basis van de aanwezigheid van een fase, gebruikt dit model een vaste set primaire variabelen—druk ($p$), specifieke enthalpie ($h$) en de totale zoutmassa fractie ($z$)—over alle faseconfiguraties heen (vloeistof, damp en vaste haliet). Deze benadering handelt faseovergangen natuurlijk af zonder handmatige wisseling, wat de numerieke robuustheid en stabiliteit verbetert.
2. Discrete Fracture-Matrix (DFM) Benadering: Het fysieke domein wordt gerepresenteerd met behulp van een mixed-dimensionale geometrie waarbij breuken en hun snijpunten worden behandeld als lager-dimensionale manifolds (1D en 0D) ingebed in het 2D porieuze matrix. Deze expliciete representatie behoudt de connectiviteit en heterogeniteit van de breuken, waardoor processen zoals door neerslag veroorzaakte reductie van de permeabiliteit lokaal kunnen worden opgelost.
3. Thermodynamische Afsluiting en Linearisering: Het model maakt gebruik van de uitgebreide correlaties van Driesner en Christoph (2007) voor H2O–NaCl thermodynamica, die magmatische condities dekken. Om de computationele kosten te verlagen, gebruikt het model een Operator-Based Linearization (OBL) strategie. Thermodynamische eigenschappen en hun afgeleiden worden vooraf berekend op een gestructureerd rooster in de ruimte van de primaire variabelen en worden via multilineaire interpolatie opgehaald tijdens de simulatie, waardoor dure 'on-the-fly' fase-scheidingsberekeningen (flash) worden geëlimineerd.
4. Dynamische Permeabiliteit en Apertuur Evolutie: Het model incorporeert Kozeny-Carman-type relaties om de matrixporositeit en -permeabiliteit dynamisch bij te werken op basis van halietverzadiging. Voor breuken evolueert de apertuur volgens een machtswet die afhankelijk is van de halietverzadiging, wat op zijn beurt de breukpermeabiliteit bijwerkt via de wet van de kubus (cubic law).

De onderliggende vergelijkingen (massa- en energiebalans) worden gediscretiseerd met behulp van een cel-gecentreerd eindige volume schema met Multi-Point Flux Approximation (MPFA) voor diffusieve termen en upwinding voor advectieve termen. Het resulterende niet-lineaire systeem wordt opgelost met de Newton-methode met een terugkerende Armijo line search en adaptieve tijdstappen.

Belangrijkste Resultaten
Het model is geverifieerd en toegepast via verschillende numerieke experimenten:

• Verificatie (1D Benchmark): Het model is geverifieerd tegenover de gevestigde closed-source simulator CSMP++ met behulp van een 1D zout-oplossingsbenchmark. De resultaten vertoonden sterke overeenkomst over vier verschillende fasegebieden (damp+haliet, damp+vloeistof, enkelvoudige fase vloeistof, en vloeistof+haliet) en hun overgangen, wat de persistent-variable formulering en de thermodynamische implementatie valideert.
• Niet-verbonden Breuknetwerk (2D): In een reservoir met geïsoleerde breuken is de stroming matrix-gedomineerd. De simulatie toonde aan dat drukverlaging bij de productiewet koken induceert, wat leidt tot significante halietneerslag in de omliggende matrix en aangrenzende breuksegmenten. Deze neerslag verminderde de matrixpermeabiliteit met maximaal 36% ($K/K_0 \approx 0,64$) en veranderde de stroomlijnen. Het verhogen van de injectiegraad intensiveerde het koken nabij de put, wat leidde tot een hogere halietverzadiging ($s_{hal} \approx 0,46$) en een ernstiger reductie van de permeabiliteit ($K/K_0 \approx 0,29$), wat de productiecijfers aanzienlijk beïnvloedde.
• Verbonden Breuknetwerk (2D): In een reservoir met een continue breukketen die injectie- en productiewells verbindt, werd het transport breuk-gedomineerd. Het verbonden netwerk faciliteerde snelle drukcommunicatie en het transport van laag-saliniteitsvloeistof naar de productiewel. Hierdoor vond weliswaar koken plaats nabij de producent, maar bleef de arriverende vloeistof voldoende verdund om lokale halietverzadiging te voorkomen. In plaats van neerslag bij de producent, werd dissolutie waargenomen, en de neerslag werd herverdeeld naar de matrix vóór de injectiefront en langs de breukketen door koelingseffecten. Deze configuratie vermeed de ernstige beperking van de permeabiliteit die in het niet-verbonden geval werd gezien.

Betekenis en Claims
Het artikel claimt een robuuste, open-source numerieke tool te bieden voor het analyseren van complexe warmte- en massatransport met minerale schalingseffecten in gefragmenteerde hoog-enthalpie geothermische systemen. De primaire betekenis ligt in:

• Numerieke Stabiliteit: De verenigde persistent-variable formulering voorkomt de convergentieproblemen die gepaard gaan met 'variable-switching' nabij fasegrenzen.
• Geometrische Getrouwheid: De DFM-benadering lost de connectiviteit van breuken expliciet op, wat de auteurs hebben aangetoond een cruciale factor is die de ruimtelijke distributie van neerslag en de resulterende transmissiviteitsreductie bepaalt.
• Operationeel Inzicht: De resultaten benadrukken dat breukconnectiviteit de operationele uitdagingen bepaalt; niet-verbonden netwerken zijn gevoelig voor verstoppingen nabij de put en verlies van productiviteit, terwijl verbonden netwerken deze problemen kunnen mitigeren door de vloeistofverdunning bij de productiewel te handhaven.
• Reproduceerbaarheid: Het werk levert het volledige numerieke model en de uitvoering workflow via een Docker-container, wat reproduceerbaarheid en verdere ontwikkeling binnen de geothermische gemeenschap faciliteert.

De auteurs concluderen dat hun model effectief patronen van halietneerslag en hun impact op de reservoirprestaties kan voorspellen, en daarmee een waardevolle tool biedt voor het begrijpen van de interactie tussen stroming, thermodynamica en minerale schaling in gefragmenteerde geothermische reservoirs.