Two-branch retention behavior in unsaturated fractured rock driven by fracture-matrix flow partitioning

Die Studie zeigt durch numerische Simulationen und analytische Herleitungen, dass der ungesättigte Fluss in gebrochenem Gestein durch eine charakteristische Zweigstruktur im Retentionsverhalten gekennzeichnet ist, die einen kritischen Übergang zwischen matrix- und bruchdominierten Strömungsregimen aufweist.

Das Wesentliche

Wasser in zerklüftetem Gestein: Warum es zwei verschiedene Wege gibt

Stellen Sie sich vor, Sie gießen einen großen, trockenen Steinhaufen. Normalerweise denken wir, dass Wasser einfach durch den ganzen Haufen sickert, wie durch einen großen Schwamm. Aber in der Natur ist das Gestein oft nicht homogen. Es besteht aus zwei Teilen:

1. Der feste Stein (die Matrix): Ein dichtes, poröses Material, das wie ein feiner Schwamm wirkt.
2. Die Risse (die Frakturen): Spalten und Klüfte im Gestein, die wie kleine Röhren oder Kanäle wirken.

Die Forscher aus Schweden und der Schweiz haben herausgefunden, dass Wasser in diesem zerklüfteten Gestein nicht einfach nur „durchsickert". Stattdessen gibt es zwei völlig verschiedene Verhaltensweisen, je nachdem, wie viel Wasser bereits im Gestein ist. Man könnte es sich wie einen Zwei-Wege-Verkehr vorstellen.

Die zwei Phasen des Wasserverkehrs

Phase 1: Der langsame Schwamm (Wenig Wasser)
Wenn das Gestein noch sehr trocken ist, passiert Folgendes: Das Wasser wird von den winzigen Poren im festen Stein (der Matrix) wie durch einen starken Magneten angezogen. Das nennt man Kapillarkraft.

• Die Analogie: Stellen Sie sich vor, Sie gießen Wasser auf einen trockenen Schwamm. Das Wasser wird sofort vom Schwamm aufgesaugt und verteilt sich langsam und gleichmäßig. Die Risse (die Kanäle) bleiben fast leer, weil das Wasser lieber in den kleinen Poren des Steins bleibt.
• Das Ergebnis: Der Wasserfluss ist langsam und wird vom „Schwamm" (der Matrix) bestimmt.

Phase 2: Die Autobahn (Viel Wasser)
Sobald der Schwamm vollgesogen ist, passiert ein magischer Wechsel. Das Wasser kann nicht mehr so schnell in den Stein eindringen, weil dieser schon voll ist. Jetzt sucht sich das Wasser den Weg des geringsten Widerstands: die Risse.

• Die Analogie: Stellen Sie sich vor, der Schwamm ist jetzt voll. Wenn Sie weiter Wasser gießen, läuft es nicht mehr in den Schwamm, sondern fließt sofort über die Oberfläche oder durch große Röhren, die Sie in den Stein geschlagen haben. Plötzlich wird das Wasser schnell transportiert, wie auf einer Autobahn im Vergleich zu einem Feldweg.
• Das Ergebnis: Der Wasserfluss wird nun von den Rissen (den Frakturen) bestimmt und ist viel schneller.

Der „Kipppunkt"

Das Spannende an dieser Studie ist, dass die Forscher genau diesen Übergangspunkt gefunden haben. Sie nennen ihn die „kritische Sättigung".

• Es gibt einen genauen Moment, an dem das Wasser aufhört, sich wie ein langsamer Schwamm zu verhalten, und anfängt, wie eine schnelle Autobahn zu fließen.
• Vorher: Langsam, verteilt, vom Stein bestimmt.
• Danach: Schnell, konzentriert, von den Rissen bestimmt.

Warum ist das wichtig?

Früher waren Wissenschaftler sich uneinig. Manche dachten, das Wasser fließt immer langsam durch den Stein. Andere sahen, dass es in Rissen sehr schnell fließt. Diese Studie löst den Streit: Beide haben recht, aber zu unterschiedlichen Zeiten.

Das hängt davon ab, wie viel Wasser da ist und wie gut die Risse miteinander verbunden sind.

• Die Verbindung der Risse: Wenn die Risse wie ein gut vernetztes Straßennetz verbunden sind, funktioniert die „Autobahn" perfekt. Sind sie aber unterbrochen (wie eine Straße, die mitten im Nirgendwo endet), muss das Wasser wieder in den „Schwamm" (den Stein) ausweichen, was den Fluss wieder verlangsamt.

Was bedeutet das für uns?

Dieses Verständnis ist wie ein neuer Bauplan für die Natur. Es hilft uns besser zu verstehen:

• Grundwasser: Wie schnell verschmutztes Wasser in den Boden gelangt.
• Kernkraftwerk-Abfälle: Wie sicher wir Abfälle in tiefen Gesteinsschichten lagern können, ohne dass Wasser sie schnell wegschwemmt.
• Wasserversorgung: Wie wir in trockenen Gebieten Wasser effizienter finden und nutzen können.

Zusammenfassend: Wasser in zerklüftetem Gestein ist kein einfaches „Sickern". Es ist ein dynamischer Prozess, der zwischen einem langsamen „Schwamm-Modus" und einem schnellen „Riss-Modus" wechselt. Sobald man weiß, wo dieser Wechsel stattfindet, kann man viel besser vorhersagen, was mit dem Wasser unter der Erde passiert.

Technische Zusammenfassung

Titel: Zweig-Verhalten der Retention in ungesättigtem, zerklüftetem Gestein, angetrieben durch die Strömungsaufteilung zwischen Klüften und Matrix

Autoren: M.R. Andiva, C. Jiang, M. Ziegler, Q. Lei (Uppsala Universität & swisstopo)

---

1. Problemstellung und Hintergrund

Die Vorhersage von ungesättigten Strömungen in zerklüftetem Gestein ist eine langjährige Herausforderung in der Hydrologie und Geophysik, mit weitreichenden Implikationen für das Grundwassermanagement und die Endlagerung radioaktiver Abfälle.

• Das Paradoxon: Unter gesättigten Bedingungen dominieren Klüfte die Strömung aufgrund ihrer hohen Permeabilität. Unter ungesättigten Bedingungen hingegen wird theoretisch erwartet, dass die Kapillarkräfte das Wasser in die feinen Poren der Gesteinsmatrix saugen (Imbibition), was die Strömung durch Klüfte unterdrückt und zu einer diffusen Infiltration führt.
• Der Konflikt: Feld- und Laborstudien (z. B. am Yucca Mountain) zeigen jedoch, dass Wasser auch unter ungesättigten Bedingungen schnell durch bevorzugte Pfade in Klüften wandern kann, oft mit minimaler Wechselwirkung mit der Matrix.
• Die Frage: Wie lässt sich der scheinbare Widerspruch zwischen kapillar-dominiertem Matrix-Fluss und klüft-dominiertem Transport in einem konsistenten physikalischen Rahmen erklären?

2. Methodik

Die Autoren kombinieren dreidimensionale direkte numerische Simulationen (DNS) mit analytischen Herleitungen, um das Verhalten zu entschlüsseln.

• Numerische Simulation:

  • Es wurden synthetische, dreidimensionale Kluftnetzwerke in einem kubischen Bereich generiert.
  • Die Kluftradien folgen einer Potenzgesetz-Verteilung ($f(r) \propto r^{-a}$).
  • Untersucht wurden verschiedene Kombinationen aus Potenzgesetz-Exponenten ($a$) und Kluftintensitäten ($\gamma$), insgesamt 25 parametrische Szenarien mit je 10 Realisierungen.
  • Die Strömung wird durch die Richards-Gleichung beschrieben. Für Matrix und Kluft werden separate Stoffgesetze (Brooks-Corey-Modell) verwendet, die die Beziehung zwischen Kapillardruck ($p_c$), effektiver Sättigung ($S_e$) und relativer Permeabilität ($k_r$) definieren.
  • Stationäre Simulationen wurden durchgeführt, wobei der Druckkopf am Einlass systematisch variiert wurde, um den gesamten Bereich der effektiven Sättigung abzudecken.

• Analytische Herleitung:

  • Basierend auf den numerischen Ergebnissen wurde eine verallgemeinerte Retentionsformel entwickelt.
  • Es wurde eine kritische Sättigung ($S_c$) und ein kritischer Druckkopf ($h_c$) analytisch hergeleitet, die den Übergang zwischen den beiden Strömungsregimen markieren.
  • Die Herleitung nutzt das Kubische Gesetz für die Kluftpermeabilität und berücksichtigt die geometrische Aufteilung zwischen Kluft und Matrix.

3. Schlüsselbeiträge und Ergebnisse

A. Emergenz eines Zweig-Verhaltens (Two-Branch Behavior)

Die numerischen Simulationen zeigen eindeutig, dass sowohl die relative Permeabilität als auch die Kapillardruckkurven eine robuste Zweig-Struktur aufweisen:

1. Matrix-dominiertes Regime (Niedrige Sättigung): Bei niedrigen effektiven Sättigungen ($S_e < S_c$) wird der Fluss primär durch die Matrix getragen. Die Kapillarkräfte ziehen das Wasser in die Matrix, und der Kluftfluss ist gering.
2. Kluft-dominiertes Regime (Hohe Sättigung): Bei höheren Sättigungen ($S_e > S_c$) wechselt das System abrupt in ein Regime, das von der Kluftströmung dominiert wird.
3. Kritischer Punkt ($S_c$): Der Übergang erfolgt bei einer kritischen Sättigung $S_c$ (im Beispiel ca. 0,37). An diesem Punkt übersteigt der Fluss durch die Klüfte den Fluss durch die Matrix.

B. Verallgemeinerte analytische Formulierung

Die Autoren leiten eine neue, verallgemeinerte Brooks-Corey-Formel ab, die dieses Zweig-Verhalten mathematisch beschreibt:

• Die Formel unterteilt die Sättigungsbereiche in zwei Teile, getrennt durch $S_c$.
• Für $S_e \leq S_c$ werden die Parameter der Matrix verwendet.
• Für $S_e > S_c$ werden die Parameter der Klüfte verwendet.
• Diese analytische Kurve stimmt hervorragend mit den numerischen Simulationsergebnissen überein und erklärt den physikalischen Mechanismus des Übergangs.

C. Physikalische Mechanismen des Übergangs

• Der Übergang wird durch den Wettbewerbsmechanismus zwischen Kapillarkräften und viskosen Kräften gesteuert.
• Die Matrix erreicht ihre volle Sättigung und maximale relative Permeabilität ($k_r \approx 1$) früher als die Klüfte.
• Sobald die Matrix gesättigt ist, nimmt der Kluftfluss weiter zu, da die Klüfte aufgrund ihrer größeren Porenweite (geringere Kapillareintrittsspannung) noch nicht ihre maximale Kapazität erreicht haben.
• Der kritische Druckkopf $h_c$ wird analytisch als Funktion der Kluftrippen (Apertur), der Kluftintensität und der Permeabilitätskontraste hergeleitet.

D. Rolle der Netzwerkkonnektivität

• Die Analyse zeigt, dass das Zweig-Verhalten auch in Netzwerken mit unterschiedlicher Konnektivität auftritt.
• Die analytische Lösung für den kritischen Druckkopf gilt streng genommen für vollständig vernetzte Systeme (oberhalb der Perkolationsschwelle $\chi_c \approx 0,7 - 2,8$).
• Unterhalb dieser Schwelle (dispergierte Netzwerke) weichen die numerischen Ergebnisse systematisch ab, da die Klüfte keine durchgehenden Pfade bilden und der Fluss stark von der Matrix-Kluft-Austauschströmung abhängt.

4. Bedeutung und Schlussfolgerungen

• Auflösung des Paradoxons: Die Studie löst den langjährigen Konflikt zwischen der Theorie des matrix-dominierten Flusses und der Beobachtung von klüft-dominiertem Transport. Sie zeigt, dass beide Regime existieren, aber durch den Sättigungszustand gesteuert werden.
• Upscaling: Die Ergebnisse bieten einen mechanistischen Rahmen für das "Upscaling" (Hochskalierung) von ungesättigten Strömungen in zerklüftetem Gestein. Anstatt ein homogenes äquivalentes poröses Medium anzunehmen, kann nun ein Modell verwendet werden, das den Übergang zwischen den beiden Regimen explizit abbildet.
• Allgemeingültigkeit: Das Konzept des Zweig-Verhaltens wurde auch mit dem van-Genuchten-Modell getestet und bleibt bestehen, was die Robustheit des Ansatzes unterstreicht.
• Praktische Implikationen: Das Verständnis der kritischen Sättigung $S_c$ ist entscheidend für die Modellierung von Grundwasserneubildung, der Migration von Schadstoffen oder radioaktiven Abfällen in ungesättigten Zonen, da sich die Transportgeschwindigkeit und -pfade drastisch ändern, sobald dieser Schwellenwert überschritten wird.

Zusammenfassend liefert diese Arbeit einen physikalisch fundierten Rahmen, der die komplexe Interaktion zwischen Kluft und Matrix in ungesättigten Systemen durch die Identifizierung eines kritischen Übergangspunkts und die Entwicklung einer entsprechenden analytischen Retentionsformel neu definiert.