Revisiting the hydromechanical formulation of a micromechanics-based phase-field model for poro-elastoplastic media

Tief unter der Erdoberfläche, in den porösen Gesteinen, die Öl, Gas und geothermische Wärme halten, wird ein ständiger, stiller Kampf zwischen festem Stein und den in seinen winzigen Poren gefangenen Fluiden ausgetragen. Wenn Ingenieure versuchen, diese Ressourcen zu extrahieren, pumpen sie oft Fluide in das Gestein, um es auseinanderzudrücken – ein Prozess, der als hydraulisches Aufbrechen (Hydraulic Fracturing) bekannt ist. Jahrzehntelang haben Wissenschaftler diesen Prozess modelliert, indem sie das Gestein als einen einfachen, elastischen Schwamm behandelten, der sich dehnt und wieder zurückfedert. Reale Gesteine sind jedoch komplexer; sie können sich auch dauerhaft verformen, zerbröckeln und entlang interner Risse gleiten, insbesondere wenn sie durch hohen Druck zusammengedrückt werden. Zu verstehen, wie genau diese Gesteine brechen, wenn sowohl das feste Skelett als auch das Innere des Fluids drücken und ziehen, ist entscheidend für das sichere Bohren von Brunnen, die Verwaltung der Abfalllagerung und sogar für die Vorhersage natürlicher geologischer Ereignisse wie dem Zerbrechen von Eisbergen oder dem Durchstoßen der Kruste durch Magma. Die Herausforderung besteht darin, eine mathematische Beschreibung zu schaffen, die die chaotische Realität des Gesteinsbruchs einfängt, ohne sich in unmöglichen Berechnungen zu verlieren.

Ein Forschungsteam hat die grundlegenden Gleichungen, die diesen Bruchprozess beschreiben, neu aufgerollt und dabei festgestellt, dass eine gängige Abkürzung, die in früheren Modellen verwendet wurde, zu einem signifikanten Fehler führt. In der Welt der Gesteinsmechanik verwenden Wissenschaftler einen „Phasenfeld“-Ansatz, bei dem ein Riss nicht als scharfe, unsichtbare Linie, sondern als eine verschwommene Zone behandelt wird, in der das Material allmählich schwächer wird. Dies ermöglicht es Computern, zu simulieren, wie ein Riss wächst und seine Form verändert. Um dies genau zu tun, muss das Modell zwei unterschiedliche Verhaltensweisen berücksichtigen: wenn das Gestein auseinandergezogen wird (Zugspannung) und wenn es zusammengedrückt wird (Druckspannung). Bei Zug entstehen winzige Risse; bei Druck schließen sich diese und die rauen Oberflächen im Inneren gleiten aneinander vorbei, was Reibung erzeugt. Die Forscher entdeckten, dass viele bestehende Modelle, die weit verbreitet in Industrie und Wissenschaft eingesetzt werden, die Variablen zur Beschreibung der Energie des Systems vertauschen. Insbesondere fanden sie heraus, dass die Behandlung des Fluiddrucks als direkter Input in der Energiegleichung – anstatt der Menge des Fluids in den Poren – zu einem mathematischen „Sprung“ oder einer Diskontinuität führt. Dieser Sprung bedeutet, dass das Modell vorhersagt, dass sich die Festigkeit des Gesteins abrupt und unrealistisch ändert, wenn es von Zug zu Druck wechselt, was zu ungenauen Vorhersagen darüber führt, wie und wo sich ein Bruch bilden wird.

Um dies zu beheben, entwickelten die Forscher eine neue, konsistente Formulierung, die die Beziehung zwischen der Deformation des Gesteins und dem Fluidgehalt sorgfältig verfolgt. Sie bauten ihr Modell auf einer mikroskopischen Sichtweise des Gesteins auf, indem sie es sich als eine Ansammlung winziger, pfenniggroßer Risse vorstellten, die mit Fluid gefüllt sind. Durch die Analyse, wie sich diese Mikro-Risse unter verschiedenen Spannungen öffnen und schließen, leiteten sie eine Reihe von Regeln ab, die sicherstellen, dass die Festigkeitsfläche des Gesteins – die mathematische Grenze, die definiert, wann es bricht – glatt und kontinuierlich bleibt, unabhängig davon, wie sich der Druck verändert. Sie bewiesen, dass das Modell durch die Verwendung einer spezifischen Art von Fließregel, die die permanente Deformation des Gesteins direkt mit den wirkenden Kräften verknüpft, und durch die korrekte Einbeziehung eines fehlenden Kopplungsterms in der Energiegleichung, den Übergang von Zug zu Druck genau beschreiben kann. Diese Korrektur ist nicht nur eine kleine Feinabstimmung; sie verändert grundlegend, wie das Modell die treibende Kraft berechnet, die einen Riss vorantreibt.

Das Team testete sein neues Modell gegen zwei sehr unterschiedliche Szenarien, um zu sehen, ob es standhält. Zuerst simulierten sie ein klassisches Benchmark-Szenario des hydraulischen Aufbrechens, bei dem Fluid in ein Gestein injiziert wird, um einen Riss zu erzeugen. Sie verglichen ihre Ergebnisse mit einer bekannten analytischen Lösung, die eine präzise mathematische Antwort für ideale Bedingungen darstellt. Das neue Modell stimmte viel genauer mit dieser Lösung überein als die älteren, gemischten Formulierungen und sagte korrekt den Druck voraus, der benötigt wird, um den Bruch zu starten, sowie die Breite des Risses während seines Wachstums. Die älteren Modelle neigten dazu, den benötigten Fluiddruck zu überschätzen, ein Fehler, der zu ineffizienten oder unsicheren Bohrungen führen könnte. Im zweiten Test simulierten sie eine Gesteinsprobe, die von zwei Seiten zusammengedrückt wird – ein Zustand, der die Umgebung tief unter der Erde nachahmt, wo Gesteine unter immensem Druck stehen. Hier zerbrach das Gestein nicht einfach; es deformierte sich plastisch und bildete eine Scherbandstruktur, bei der das Material aneinander vorbeiglitt, bevor es brach. Das neue Modell reproduzierte dieses komplexe Verhalten erfolgreich und erfasste sowohl die durch mechanisches Zusammendrücken verursachten Scherbrüche als auch die durch Fluidinjektion getriebenen Zugbrüche.

Die Ergebnisse dieser Simulationen zeigen, dass die plastische Deformation des Gesteins als bedeutende Barriere für das Risswachstum wirkt. Wenn das Gestein in der Lage ist, sich permanent zu verformen, absorbiert es mehr Energie, was bedeutet, dass der Fluiddruck innerhalb des Risses höher sein muss, um den Bruch offen zu halten und wachsen zu lassen. Dieser zusätzliche Widerstand verändert das gesamte Profil des Bruchs und macht ihn schmaler und schwerer voranzutreiben, als es ein rein elastisches Modell vermuten ließe. Darüber hinaus zeigte das Modell in den Kompressionstests, wie sich der Porendruck innerhalb des Gesteins verändert, während der Bruch entsteht. Während sich das Gestein deformiert und Raum für die Bewegung des Fluids schafft, sinkt der Druck lokal, wodurch mehr Fluid in die beschädigte Zone gesogen wird. Diese Wechselwirkung erzeugt eine Rückkopplungsschleife, in der der Fluidfluss und das mechanische Versagen des Gesteins eng miteinander gekoppelt sind – eine Dynamik, die das neue Modell mit hoher Genauigkeit erfasst.

Durch die Korrektur des mathematischen Fundaments dieser Modelle haben die Forscher ein zuverlässigeres Werkzeug für Ingenieure und Wissenschaftler bereitgestellt, die mit porösen Medien arbeiten. Die neue Formulierung stellt sicher, dass der Übergang zwischen verschiedenen Arten des Gesteinsversagens glatt und physikalisch realistisch ist und eliminiert die künstlichen Sprünge, die bisherige Versuche plagten. Diese Genauigkeit ist entscheidend für Anwendungen, die von der Optimierung der Öl- und Gasproduktion bis hin zum Entwurf sicherer unterirdischer Speicher für die Kohlenstoffabscheidung reichen. Die Studie zeigt, dass selbst in komplexen Systemen, in denen Festkörper und Fluide interagieren, ein sorgfältiger, konsistenter Ansatz zur zugrunde liegenden Physik die wahre Natur des Materialbruchs offenbaren kann. Das Modell steht nun bereit, das Verhalten der Erdkruste mit größerem Vertrauen vorherzusagen und die Lücke zwischen theoretischer Mechanik und der chaotischen, plastischen Realität des Bodens unter unseren Füßen zu schließen.