Residual Saturation under Pressure-Controlled Drainage

Technische Samenvatting: Residuele Verzadiging onder Drukgecontroleerde Drainage

Probleemstelling
Het artikel behandelt een fundamentele, nog onopgeloste vraag in de fysica van poreuze media: het gedrag van residuele verzadiging onder drukgecontroleerde drainage in de oneindige grootte-limiet. Terwijl standaard invasie-percolatie (IP)-modellen voorspellen dat de verzadiging van de invallende fase naar nul gaat wanneer de systeemgrootte naar oneindig gaat (waarbij slechts een eindig deel van de verdedigende vloeistof gevangen blijft door de fractale aard van de doorbraakcluster), blijft het onduidelijk of een drukgecontroleerd proces—dat de invallende fase toegang geeft tot alle poriën die verbonden zijn met de inlaat bij een bepaalde capillaire drempel—resulteert in een eindige residuele verzadiging van de verdedigende fase in de thermodynamische limiet. Dit onderscheid is cruciaal voor het begrijpen van fenomenen op veldschaal, zoals waterinjectie, waarbij de vraag rijst of de residuele olieverzadiging eindig blijft wanneer het geïndeerde gebied de afmetingen van het reservoir nadert.

Methodologie
De auteurs formuleren drukgecontroleerde drainage als Bond Percolation with Trapping (BPT) op een poriënetwerk-graaf.

• Graafrepresentatie: Het poreuze medium wordt gemodelleerd als een ongerichte graaf waarbij knopen porieregio's vertegenwoordigen en randen (edges) doorgangen (throats) zijn met toegewezen ingangsradii (capillaire drempels).
• Procesdefinitie: In tegenstelling tot standaard IP, dat zich door de meest lokaal toegankelijke weg door de doorgangen beweegt, opereert BPT onder een globale drukrestrictie. Naarmate de capillaire druk toeneemt, breidt de "actieve graaf" (bestaande uit doorgangen met ingangsdrempels onder de huidige druk) zich uit.
• Trapping-mechanisme: Invasie wordt gedefinieerd door connectiviteit met de inlaat op de actieve graaf. Trapping (vangst) vindt plaats wanneer een deel van de verdedigende fase de connectiviteit met de uitlaat verliest. Het algoritme werkt het systeem bij door: (1) de druk te verhogen om nieuwe doorgangen te activeren; (2) de geïndeerde regio te identificeren (inlaat-verbonden actieve poriën); (3) de verdedigende regio te identificeren (uitlaat-verbonden niet-geïndeerde poriën); en (4) alle resterende gedisconnecteerde poriën als "getrapt" te classificeren.
• Simulaties: De studie vergelijkt BPT met standaard Invasie-Percolatie (IP) op diverse 2D-roosters (vierkant, driehoekig, Voronoi) en 3D-roosters (diamant, simpel kubisch met variërende coördinatiegetallen). Simulaties maken gebruik van lognormaal verdeelde doorgangsradii en ensemble-gemiddelden om effecten van eindige grootte te analyseren.

Belangrijkste Bijdragen en Resultaten

1. Eindige Residuele Verzadiging in de Oneindige Limiet:
   De belangrijkste bevinding is dat drukgecontroleerde drainage leidt tot een ruimtevullende geïndeerde staat, waarbij de residuele verzadiging van de verdedigende fase een niet-nul, eindige limiet benadert ($0 < S_\infty < 1$) wanneer de systeemgrootte $L \to \infty$. Dit contrasteert scherp met standaard IP, waarbij de geïndeerde verzadiging bij doorbraak verdwijnt ($S_{IP} \to 0$) omdat de invallende cluster fractaal is en grote regio's van de verdedigende vloeistof omzeilt.

2. Finite-Size Scaling in Twee Dimensies:
   In 2D volgt de afwijking van de residuele verzadiging van de thermodynamische limiet een finite-size scaling-wet:
   $$S_r(L) \approx S_\infty - \left(\frac{L_0}{L}\right)^\delta$$
   waarbij de exponent $\delta \approx 0,25$. Deze exponent is universeel over verschillende rooster-typen (vierkant, driehoekig, Voronoi) en is onafhankelijk van microscopische details. De trage convergentie ($\Delta S_r \sim L^{-0,25}$) impliceert dat het bereiken van representatieve elementaire volumes in 2D extreem grote systeemgroottes vereist (bijvoorbeeld: het verdubbelen van de systeemgrootte vermindert de bias van de eindige grootte slechts met een factor $\approx 1/16$).

3. Dimensionaliteit en Coördinatiegetal-effecten in Drie Dimensies:
   In 3D vervallen de correcties voor eindige grootte aanzienlijk sneller, met een exponent $\delta \approx 0,75$. Bijgevolg is het representatieve volume dat nodig is om de oneindige grootte-limiet te benaderen in 3D veel kleiner dan in 2D. Bovendien, hoewel de asymptotische residuele verzadiging in 3D eindig blijft, hangt de waarde ervan af van het coördinatiegetal van het rooster; hogere connectiviteit biedt meer alternatieve invasiepaden, wat het gevangen deel vermindert, maar het niet volledig elimineert.

4. Onderscheid van Invasie-Percolatie:
   Het artikel verduidelijkt dat het verschil tussen BPT en IP niet louter een kwestie is van de vorm van de curve of burst-dynamiek, maar een fundamenteel verschil in asymptotische staten. IP levert een fractale invallende cluster bij doorbraak, terwijl BPT een ruimtevullende geïndeerde regio oplevert met een eindig getrapt deel dat in de oneindige limiet coëxisteert.

Betekenis
Het artikel legt een direct theoretisch verband tussen percolatietheorie en druk-verzadigingsrelaties door aan te tonen dat connectiviteit en trapping de bepalende factoren zijn voor residuele verzadiging onder drukgecontroleerde drainage. Door aan te tonen dat een eindig deel van de porieruimte gevangen blijft, zelfs in de oneindige systeemlimiet, breidt het werk het standaard invasie-percolatiebeeld uit voorbij de doorbraakstaat. Dit biedt een rigoureus kader voor het interpreteren van macroscopische constitutieve wetten en suggereert dat het "representatieve elementaire volume"-probleem zeer gevoelig is voor dimensionaliteit, waarbij het in 2D veel restrictiever is dan in 3D. De resultaten impliceren dat in drukgecontroleerde scenario's een eindig deel van de verdedigende vloeistof onweerlegbaar onherstelbaar is, ongeacht de systeemgrootte, een bevinding met directe implicaties voor de herverdeling van ondergrondse vloeistoffen en winningsprocessen.