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

Diese Arbeit präsentiert ein neues Open-Source-Modell für kompositionelle Strömungen, das im PorePy-Framework implementiert ist und nicht-isotherme, multiphasische Strömung sowie Halitpräzipitation in hochenthalpischen, geklüfteten geothermischen Reservoirs simuliert, wobei ein robustes Primärvariablen-Formulierungssystem und ein Diskretisierungsansatz für Kluft-Matrix-Systeme genutzt werden, um Permeabilitätsschäden und betriebliche Herausforderungen präzise vorherzusagen.

Das Wesentliche

Stellen Sie sich ein heißes unterirdisches Reservoir wie einen riesigen, natürlichen Schnellkochtopf vor, der mit superheißem Salzwasser gefüllt ist. Dies ist die Quelle der Geothermie. Um diese Energie zu gewinnen, bohren Ingenieure Brunnen und pumpen kaltes Wasser hinein, um das heiße Wasser wieder nach oben zu drücken. Dieser Prozess ist jedoch aufgrund des „Salzes“ im Wasser schwierig.

Betrachten Sie das Salz (Halit) wie Zucker in einer Tasse heißem Tee. Wenn Sie den Tee abkühlen lassen oder etwas Wasser verdampfen lassen, kann der Zucker nicht mehr gelöst bleiben und beginnt, wieder zu festen Kristallen zu werden. In einem geothermischen Bohrloch passiert dies, wenn das heiße Wasser in der Nähe eines Injektionsbrunnens abkühlt oder in der Nähe eines Produktionsbrunnens verdampft. Das Ergebnis? Feste Salzkristalle bilden sich und verstopfen die winzigen Löcher in den Gesteinen und den Rissen (Klüften), die das Wasser fließen lassen. Es ist wie ein Stau, der durch Zuckerkristalle verursacht wird, die eine Autobahn blockieren.

Dieses Paper stellt ein neues Computersimulationswerkzeug vor, das genau vorhersagen soll, wo und wie diese „Zuckerstaus“ in komplexen, rissigen Gesteinen auftreten werden.

Hier ist eine Aufschlüsselung, wie das Werkzeug funktioniert und was es herausgefunden hat, unter Verwendung einfacher Analogien:

1. Die Karte: Die Risse klar sehen

Traditionelle Karten von unterirdischen Gesteinen glätten oft alles aus und behandeln das Gestein wie einen festen Schwamm. In der Realität fließt das Wasser jedoch hauptsächlich durch ein Netzwerk von Rissen, ähnlich wie Wasser durch einen rissigen Gehweg fließt und nicht durch den Beton selbst.

• Die Innovation: Dieses neue Modell verwendet einen „Discrete Fracture-Matrix“-Ansatz (diskret-frakturiertes Matrix-Modell). Stellen Sie sich vor, man zeichnet die Risse als deutliche, dünne Linien auf eine Karte, anstatt sie einfach in den Hintergrund zu verschwimmen. Dies ermöglicht es dem Computer zu sehen, wie genau die Risse miteinander verbunden sind (oder nicht), und wie Salz einen spezifischen Riss im Vergleich zum umliegenden Gestein verstopfen könnte.

2. Der Motor: Eine „Universalfernbedienung“ für die Physik

Die gleichzeitige Simulation von kochendem Wasser, Dampf und festem Salz ist für Computer unglaublich schwierig. Normalerweise muss der Computer ständig seinen „Modus“ wechseln (z. B. „Okay, jetzt ist es flüssig; jetzt ist es Gas; jetzt ist es fest“), was dazu führen kann, dass die Berechnung abstürzt oder stecken bleibt.

• Die Innovation: Die Autoren haben ein „einheitliches“ System geschaffen. Denken Sie an eine Universalfernbedienung, die für jedes Gerät funktioniert, ohne dass man die Batterien oder Modi wechseln muss. Das Modell verwendet drei feste „Regler“ (Druck, thermische Energie und Salzmenge), die gleich bleiben, egal ob das Wasser flüssig, dampfförmig oder in festes Salz umgewandelt ist. Dies macht die Simulation viel reibungslender und stabiler, sodass sie den chaotischen Wechsel zwischen den Zuständen bewältigen kann, ohne abzustürzen.

3. Der Geschwindigkeitstrick: Das „Spickzettel“-Prinzip

Die Berechnung der exakten Physik von Salzwasser bei hohen Temperaturen erfordert normalerweise, dass der Computer immer wieder komplexe mathematische Rätsel löst, was sehr langsam ist.

• Die Innovation: Das Team hat einen vorab berechneten „Spickzettel“ (Lookup-Tabelle) erstellt. Bevor die Simulation startet, haben sie alle möglichen Ergebnisse berechnet, wie sich das Salz unter verschiedenen Bedingungen verhält, und diese gespeichert. Während der Simulation muss der Computer nicht jedes Mal die schwere Mathematik lösen, sondern schaut einfach die Antwort auf dem Spickzettel nach. Dies macht die Simulation viel schneller, während sie gleichzeitig präzise bleibt.

4. Der Verstopfungseffekt: Das „Porengrößen-Schrumpfen“

Wenn sich Salzkristalle bilden, nehmen sie Platz ein.

• Die Innovation: Das Modell verkleinert die „Rohre“ (Porosität und Permeabilität) automatisch, wenn sich Salz ansammelt. Es verwendet eine Regel (Kozeny-Carman), die besagt: „Wenn Salz 10 % des Lochs ausfüllt, wird das Rohr signifikant schmaler.“ Dies ermöglicht es dem Modell vorherzusagen, wie der Fluss langsamer wird oder ganz zum Erliegen kommt, wenn der „Zuckerstau“ schlimmer wird.

Was die Simulationen zeigten

Das Team testete dieses Werkzeug in zwei Hauptszenarien:

Szenario A: Die kaputte Autobahn (Nicht verbundene Risse)

• Setup: Stellen Sie sich ein Reservoir vor, in dem die Risse nicht untereinander verbunden sind; das Wasser muss durch das feste Gestein zwischen ihnen quetschen.
• Ergebnis: Als sie kaltes Wasser einspeisten, kochte das heiße Wasser in der Nähe des Produktionsbrunnens schnell auf. Dies verursachte die Kristallisation von Salz und verstopfte das Gestein direkt um den Brunnen herum.
• Die Wendung: Wenn sie das Wasser schneller einspeisten, wurde die Verstopfung viel schlimmer und die Energieausbeute sank erheblich. Das Modell zeigte, dass der „Stau“ hauptsächlich im Gestein in der Nähe des Brunnens stattfand, nicht nur in den Rissen.

Szenario B: Die verbundene Autobahn (Verbundene Risse)

• Setup: Stellen Sie sich ein Reservoir vor, in dem die Risse eine kontinuierliche Hochgeschwindigkeitsautobahn vom Injektionsbrunnen zum Produktionsbrunnen bilden.
• Ergebnis: Das kalte Wasser strömte schnell durch die Risse. Da es sich so schnell bewegte und kühl blieb, löste es das Salz in der Nähe des Produktionsbrunnens tatsächlich auf, anstatt es zu verstopfen!
• Die Wendung: Die Salzpräzipitation verlagerte sich an einen anderen Ort – direkt an die Grenze der kalten Wasserzone, anstatt den Brunnen selbst zu verstopfen. Dies deutet darauf hin, dass ein vernetztes System von Rissen den Brunnen vor Verstopfungen schützen könnte, auch wenn es die Stelle verändert, an der sich das Salz ansammelt.

Das Faz-it

Dieses Paper stellt ein neues Open-Source-Softwaretool vor, das Ingenieuren hilft, den komplexen Tanz zwischen Hitze, Druck und Salz in geothermischen Brunnen zu verstehen. Durch die genaue Kartierung, wie Risse zusammenhängen und wie Salz sie verstopft, kann das Tool helfen vorherzusagen:

1. Wo Brunnen durch Salz blockiert werden könnten.
2. Wie viel Energie sicher extrahiert werden kann, bevor die „Rohre“ verstopfen.
3. Ob das Layout der unterirdischen Risse den Produktionsprozess unterstützen oder behindern wird.

Die Autoren haben ihr Werkzeug gegen einen etablierten Industriestandard überprüft und festgestellt, dass es perfekt übereinstimmt, was beweist, dass es ein zuverlässiger Weg ist, um diese hochtemperierten, salzhaltigen und rissigen Gesteinsumgebungen zu simulieren.

Technische Zusammenfassung

Technisches Resümee: Einheitliches kompositionelles Strömungsmodell zur Simulation mehrphasiger, frakturierter Hochenthalpie-Systeme

Problemstellung
Die Simulation von Hochenthalpie-Geothermie-Reservoirs mit Frakturen stellt eine erhebliche Herausforderung dar, da sie die komplexe Kopplung von nicht-isothermer, mehrphasiger, multikomponentiger Strömung, stark nichtlinearer Thermodynamik und dem dominierenden Einfluss von Kluftnetzwerken umfasst. Ein kritisches betriebliches Problem in diesen Systemen ist die Mineralienausfällung, insbesondere die Halit-Präzipitation (NaCl). Während der Produktion führt der Druckabfall zu Sieden, wodurch die verbleibende Sole konzentriert wird und die Löslichkeitsgrenzen überschritten werden. Gleichzeitig senkt die Kaltwasserinjektion die lokale Löslichkeit. Beide Mechanismen führen zur Salzabscheidung in Porenräumen und Klüften, was die Porosität und Permeabilität schwerwiegend beeinträchtigt, die Produktivität der Bohrung reduziert und potenziell zu Verstopfungen des Bohrlochs führt. Bestehende Kontinuumsmodelle versagen oft bei der Erfassung lokalisierter Ausfällungen und der frakturskaligen Permeabilitätsentwicklung, während traditionelle kompositionelle Simulatoren, die auf Variable-Switching-Formulierungen basieren, Schwierigkeiten mit der numerischen Stabilität nahe der Phasengrenzen haben.

Methodik
Die Autoren präsentieren ein neues, einheitliches kompositionelles Strömungsmodell, das innerhalb des Open-Source-Frameworks PorePy implementiert wurde. Die Methodik integriert drei Kernkomponenten:

1. Persistent-Variable-Formulierung: Im Gegensatz zu herkömmlichen Ansätzen, die die primären Variablen basierend auf der Anwesenheit von Phasen wechseln, verwendet dieses Modell einen festen Satz an Primärvariablen – Druck ($p$), spezifische Enthalpie ($h$) und den gesamten Salzmasseanteil ($z$) – über alle Phasenkonfigurationen (flüssig, dampfförmig und fester Halit) hinweg. Dieser Ansatz handhabt Phasenübergänge auf natürliche Weise ohne manuelles Umschalten, was die numerische Robustheit und Stabilität erhöht.
2. Discrete Fracture-Matrix (DFM)-Ansatz: Das physikalische Gebiet wird mittels einer gemischten dimensionalen Geometrie dargestellt, bei der Klüfte und deren Schnittpunkte als niederdimensionale Mannigfaltigkeiten (1D und 0D) in die 2D-poröse Matrix eingebettet behandelt werden. Diese explizite Darstellung bewahrt die Konnektivität und Heterogenität der Klüfte, was die Auflösung lokalisierter Prozesse wie der durch Präzipitation verursachten Permeabilitätsreduktion ermöglicht.
3. Thermodynamische Abschließung und Linearisierung: Das Modell nutzt die umfassenden Korrelationen von Driesner und Christoph (2007) für die H2O–NaCl-Thermodynamik, die magmatische Bedingungen abdecken. Zur Reduzierung des Rechenaufwands verwendet das Modell eine operatorbasierte Linearisierung (OBL). Thermodynamische Eigenschaften und deren Ableitungen werden auf einem strukturierten Gitter im Raum der Primärvariablen vorberechnet und während der Simulation mittels multilinearer Interpolation abgerufen, wodurch aufwendige Flash-Berechnungen zur Phasenseparation während der Laufzeit vermieden werden.
4. Dynamische Permeabilität und Apertur-Evolution: Das Modell beinhaltet Kozeny-Carman-Typ-Beziehungen, um die Matrixporosität und -permeabilität basierend auf der Halitsättigung dynamisch zu aktualisieren. Für Klüfte entwickelt sich die Apertur gemäß einer Potenzgesetz-Relation abhängig von der Halitsättigung, was wiederum die Kluftpermeabilität über das kubische Gesetz aktualisiert.

Die Erhaltungsgleichungen (Massen- und Energiebilanz) werden mittels eines zellzentrierten Finite-Volumen-Verfahrens diskretisiert, wobei die Multi-Point Flux Approximation (MPFA) für diffusive Terme und Upwinding für advektive Terme zum Einsatz kommen. Das resultierende nichtlineare System wird mittels der Newton-Methode mit einer Backtracking-Armijo-Line-Search und adaptiver Zeitschrittsteuerung gelöst.

Wesentliche Ergebnisse
Das Modell wurde durch mehrere numerische Experimente verifiziert und angewendet:

• Verifizierung (1D-Benchmark): Das Modell wurde gegen den etablierten proprietären Simulator CSMP++ anhand eines 1D-Salzdissolutions-Benchmarks verifiziert. Die Ergebnisse zeigten eine starke Übereinstimmung über vier verschiedene Phasenregionen (Dampf+Halit, Dampf+Flüssigkeit, einphasige Flüssigkeit und Flüssigkeit+Halit) sowie deren Übergänge, was die Persistent-Variable-Formulierung und die thermodynamische Implementierung validiert.
• Diskonnektiertes Kluftnetzwerk (2D): In einem Reservoir mit isolierten Klüften ist die Strömung matrixdominiert. Die Simulation zeigte, dass der Druckabfall am Produktionsbrunnen das Sieden induziert, was zu einer signifikanten Halit-Präzipitation in der umgebenden Matrix und angrenzenden Kluftsegmenten führt. Diese Präzipitation reduzierte die Matrixpermeabilität um bis zu 36 % ($K/K_0 \approx 0,64$) und veränderte die Strömungsstrukturen. Eine Erhöhung der Injektionsrate intensivierte das Sieden nahe des Bohrlochs, was zu einer höheren Halitsättigung ($s_{hal} \approx 0,46$) und einer schwerwiegenderen Permeabilitätsreduktion ($K/K_0 \approx 0,29$) führte, was die Produktionsraten erheblich beeinflusste.
• Verbundenes Kluftnetzwerk (2D): In einem Reservoir mit einer kontinuierlichen Kluftkette, die die Injektions- und Produktionsbrunnen verbindet, wurde der Transport kluftdominiert. Das verbundene Netzwerk ermöglichte eine schnelle Druckkommunikation und den Transport von fluiden mit geringem Salzgehalt zum Produktionsbrunnen. Infolgedessen trat zwar nahe dem Produzenten Sieden auf, aber das ankommende Fluid war ausreichend verdünnt, um eine lokale Halitsättigung zu verhindern. Anstatt einer Präzipitation am Produzenten wurde eine Dissolution beobachtet, und die Präzipitation wurde aufgrund von Abkühlungseffekten in die Matrix vor der Injektionsfront und entlang der Kluftkette umverteilt. Diese Konfiguration vermeidete die schwere Beeinträchtigung der Permeabilität, wie sie im diskonnektierten Fall auftrat.

Bedeutung und Ansprüche
Das Paper beansprucht, ein robustes, Open-Source-Numerikwerkzeug zur Analyse komplexer Wärme- und Stofftransportprozesse mit Mineralienausfällung in hochenthalpischen, frakturierten Geothermie-Systemen bereitzustellen. Die primäre Bedeutung liegt in:

• Numerischer Stabilität: Die einheitliche Persistent-Variable-Formulierung vermeidet die Konvergenzprobleme, die mit dem Wechsel der Variablen nahe der Phasengrenzen einhergehen.
• Geometrischer Treue: Der DFM-Ansatz bildet die Konnektivität der Klüfte explizit ab, was die Autoren demonstrieren, dass ein entscheidender Faktor für die räumliche Verteilung der Ausfällung und die daraus resultierende Transmissibilitätsreduktion ist.
• Operativer Erkenntnisgewinn: Die Ergebnisse verdeutlichen, dass die Kluftkonnektivität die operativen Herausforderungen bestimmt; diskonnektierte Netzwerke sind anfällig für Verstopfungen im Nahbereich des Bohrlochs und Produktivitätsverluste, während verbundene Netzwerke diese Probleme mildern können, indem sie die Fluidverdünnung am Produktionsbrunnen aufrechterhalten.
• Reproduzierbarkeit: Die Arbeit stellt das vollständige numerische Modell und den Ausführungs-Workflow über einen Docker-Container bereit, was die Reproduzierbarkeit und die Weiterentwicklung innerhalb der Geothermie-Gemeinschaft erleichtert.

Die Autoren kommen zu dem Schluss, dass ihr Modell die Halit-Präzipitationsmuster und deren Auswirkungen auf die Reservoirleistung effektiv vorhersagen kann und somit ein wertvolles Werkzeug für das Verständnis des Zusammenspiels von Strömung, Thermodynamik und Mineralienausfällung in frakturierten Geothermie-Reservoirs darstellt.