Residual Saturation under Pressure-Controlled Drainage

Technisches Resümee: Residuelle Sättigung unter druckgesteuerter Entleerung

Problemstellung
Die Arbeit befasst sich mit einer grundlegenden, bislang ungelösten Frage in der Physik poröser Medien: dem Verhalten der residuellen Sättigung unter druckgesteuerter Entleerung (drainage) im unendlich großen Grenzfall. Während Standardmodelle der Invasionsperkolation (Invasion Percolation, IP) vorhersagen, dass die Sättigung der eindringenden Phase gegen Unendlich verschwindet (wobei nur aufgrund der fraktalen Natur des Durchbruchsclusters ein endlicher Anteil der verteidigenden Flüssigkeit eingeschlossen bleibt), bleibt unklar, ob ein druckgesteuerter Prozess – der es der eindringenden Phase ermöglicht, auf alle Poren zuzugreifen, die bei einem gegebenen Kapillardruck mit dem Einlass verbunden sind – zu einer endlichen residuellen Sättigung der verteidigenden Phase im thermodynamischen Limes führt. Diese Unterscheidung ist entscheidend für das Verständnis großskaliger Phänomene wie der Wasserflutung (waterflooding), bei der die Frage aufkommt, ob die residuelle Öl-Sättigung endlich bleibt, wenn sich die eindringende Region den Dimensionen des Reservoirs annähert.

Methodik
Die Autoren formulieren die druckgesteuerte Entleerung als Bond-Perkolation mit Einschluss (Bond Percolation with Trapping, BPT) auf einem Porenetzwerk-Graphen.

• Graph-Repräsentation: Das poröse Medium wird als ungerichteter Graph modelliert, bei dem Knoten die Porenbereiche und Kanten die Engstellen (throats) mit zugewiesenen Eintrittsradien (Kapillarschwellenwerten) darstellen.
• Prozessdefinition: Im Gegensatz zur Standard-IP, die sich durch den lokal am leichtesten zugänglichen Pfad engstelle für Engstelle vorarbeitet, operiert BPT unter einer globalen Druckbeschränkung. Wenn der Kapillardruck steigt, dehnt sich der „aktive Graph“ (der aus Engstellen mit Eintrittsschwellenwerten unter dem aktuellen Druck besteht) aus.
• Einschlussmechanismus (Trapping): Die Invasion wird über die Konnektivität zum Einlass auf dem aktiven Graphen definiert. Ein Einschluss tritt ein, wenn ein Teil der verteidigenden Phase die Konnektivität zum Auslass verliert. Der Algorithmus aktualisiert das System durch: (1) Erhöhung des Drucks, um neue Engstellen zu aktivieren; (2) Identifizierung der eindringenden Region (einlass-verbundene aktive Poren); (3) Identifizierung der verteidigenden Region (auslass-verbundene nicht-invadierte Poren); und (4) Klassifizierung aller verbleibenden, nicht verbundenen Poren als „eingeschlossen“ (trapped).
• Simulationen: Die Studie vergleicht BPT mit der Standard-Invasionsperkolation (IP) auf verschiedenen 2D-Gittern (quadratisch, dreieckig, Voronoi) und 3D-Gittern (Diamant-Gitter, einfach kubisch mit variierenden Koordinationszahlen). Es werden Simulationen mit lognormal verteilten Engstellenradien und Ensemble-Mittelungen durchgeführt, um Finite-Size-Effekte zu analysieren.

Wesentliche Beiträge und Ergebnisse

1. Endliche residuelle Sättigung im unendlichen Limes:
   Die zentrale Erkenntnis ist, dass die druckgesteuerte Entleerung zu einem raumfüllenden eindringenden Zustand führt, bei dem die residuelle Sättigung der verteidigenden Phase gegen einen endlichen, von Null verschiedenen Grenzwert ($0 < S_\infty < 1$) konvergiert, wenn die Systemgröße $L \to \infty$ geht. Dies steht in starkem Kont contrast zur Standard-IP, bei der die eindringende Sättigung beim Durchbruch verschwindet ($S_{IP} \to 0$), da der eindringende Cluster fraktal ist und große Regionen der verteidigenden Flüssigkeit umgeht.

2. Finite-Size-Scaling in zwei Dimensionen:
   In 2D folgt die Abweichung der residuellen Sättigung von ihrem thermodynamischen Limes einem Finite-Size-Scaling-Gesetz:
   $$S_r(L) \approx S_\infty - \left(\frac{L_0}{L}\right)^\delta$$
   wobei der Exponent $\delta \approx 0,25$ ist. Dieser Exponent ist über verschiedene Gittertypen (quadratisch, dreieckig, Voronoi) hinweg universell und unabhängig von mikroskopischen Details. Die langsame Konvergenz ($\Delta S_r \sim L^{-0,25}$) impliziert, dass das Erreichen repräsentativer Elementarvolumina in 2D extrem große Systemgrößen erfordert (z. B. reduziert eine Verdoppelung der Systemgröße den Finite-Size-Bias nur um den Faktor $\approx 1/16$).

3. Dimensionalitäts- und Koordinationszahl-Effekte in drei Dimensionen:
   In 3D fallen die Finite-Size-Korrekturen signifikant schneller aus, mit einem Exponenten von $\delta \approx 0,75$. Folglich ist das repräsentative Volumen, das zur Annäherung an den unendlichen Limes erforderlich ist, in 3D wesentlich kleiner als in 2D. Darüber hinaus bleibt der asymptotische residuelle Sättigungswert in 3D zwar endlich, hängt jedoch von der Koordinationszahl des Gitters ab; eine höhere Konnektivität bietet mehr alternative Invasionspfade, was den eingeschlossenen Anteil reduziert, ihn jedoch nicht vollständig eliminiert.

4. Abgrenzung zur Invasionsperkolation:
   Die Arbeit stellt klar, dass der Unterschied zwischen BPT und IP nicht bloß eine Frage der Kurvenform oder der Burst-Dynamik ist, sondern einen fundamentalen Unterschied in den asymptotischen Zuständen darstellt. IP liefert beim Durchbruch einen fraktalen eindringenden Cluster, während BPT eine raumfüllende eindringende Region mit einem endlichen eingeschlossenen Anteil liefert, die im unendlichen Limes koexistieren.

Bedeutung
Die Arbeit stellt eine direkte theoretische Verbindung zwischen der Perkolationstheorie und den Druck-Sättigungs-Beziehungen her, indem sie zeigt, dass Konnektivität und Einschluss die bestimmenden Faktoren für die residuelle Sättigung unter druckgesteuerter Entleerung sind. Durch den Nachweis, dass selbst im unendlichen System ein endlicher Anteil des Porenraums eingeschlossen bleibt, erweitert die Arbeit das Standardbild der Invasionsperkolation über den Durchbruchzustand hinaus. Dies bietet einen rigorosen Rahmen für die Interpretation makroskopischer konstitutiver Gesetze und legt nahe, dass das Problem des „repräsentativen Elementarvolumens“ hochgradig sensitiv gegenüber der Dimensionalität ist und in 2D weita viel restriktiver ist als in 3D. Die Ergebnisse implizieren, dass in druckgesteuerten Szenarien ein endlicher Anteil der verteidigenden Flüssigkeit unabhängig von der Systemgröße ungewinnbar bleibt, was direkte Auswirkungen auf die Umverteilung von Fluiden im Untergrund und auf Förderprozesse hat.