Mixing and spreading of gravity currents in heterogeneous porous media

Deze studie maakt gebruik van hoogwaardige numerieke simulaties om aan te tonen dat hoewel heterogeniteit in permeabiliteit de migratie van zwaartekrachtsstromen over het algemeen versnelt en de oplosbaarheid door dispersieve effecten verhoogt, het samenspel tussen de grootte van de instabiliteit, de correlatielengte en de dichtheidsstratificatie uiteindelijk de oplosbaarheidsefficiëntie bepaalt, waarbij homogene media vaak beter presteren dan heterogene, behalve bij lage Rayleigh-getallen.

De kern

Stel je voor dat je een dikke, zware siroop in een glas water giet. In een perfect, uniform glas zou de siroop langs de zijkant glijden in een glad, voorspelbaar vlak. Maar wat als het glas niet leeg was? Wat als het volgepropt zat met een spons die gaten had ter grootte van een speld en andere ter grootte van een marmeren balletje?

Dit is de wereld van zwaartekrachtstromen in porieuze media (zoals ondergrondse gesteenten of bodem) die Albert Jiménez-Ramos en Juan J. Hidalgo onderzochten in hun paper. Ze gebruikten krachtige computersimulaties om te observeren hoe vloeistoffen mengen en bewegen door deze "sponsachtige" ondergrondse lagen, met name gericht op hoe de ongelijkmatigheid van het gesteente het spel verandert.

Hier is het verhaal van hun bevindingen, opgesplitst in alledaagse concepten:

De Opzet: Twee Soorten Vloeistoffen

De onderzoekers bestudeerden twee scenario's, als twee verschillende soorten feesten:

1. Het Stabiele Feest: Stel je voor dat je olie in water giet. De olie is lichter en blijft bovenop, of als het zout water is, zakt het soepel. De vloeistoffen vechten niet tegen elkaar; ze mengen zich gewoon langzaam aan de grens.
2. Het Instabiele Feest: Dit is als het gieten van een vloeistof die zwaar is in het midden maar licht aan de randen. Het is een chaotische situatie waarbij de zware delen willen zakken en de lichte delen willen drijven, waardoor er "vingers" van vloeistof ontstaan die duiken of omhoog schieten, en alles gewelddadig mengen.

Het "Spons"-effect: Heterogeniteit

In de echte wereld is ondergronds gesteente niet uniform. Het is een heterogeen mengsel van harde, strakke plekken (lage permeabiliteit) en losse, open plekken (hoge permeabiliteit). De onderzoekers behandelden dit als een spons met willekeurige gaten.

Wat ze vonden:

• Het Barrière-effect: Wanneer de vloeistof een strakke, harde plek in het gesteente raakt, blijft het steken. Het is alsof je probeert door een menigte te rennen; als er een muur is (een zone met lage permeabiliteit), moet je eromheen. Deze "barrière" vertraagt het mengproces meestal, omdat de vloeistof er niet gemakkelijk doorheen kan komen.
• De Valstrik: Soms blijft de lichte vloeistof vastzitten in een strak zakje omringd door zware vloeistof. Het is als een bel die vastzit in een net. Uiteindelijk lost deze vastzittende bel snel op, waardoor er een kleine explosie van menging ontstaat.

De Grote Verrassing: Chaos versus Orde

De meest interessante ontdekking was hoe de "spons" (het gesteente) interacteerde met de "vingers" (de instabiele menging).

• In het Stabiele Geval: De ongelijkmatige gesteente fungeerde als een verspreider. Het spreidde de vloeistof uit, waardoor de mengzone breder en trager werd. Het was alsof je door een bos rent; je raakt verspreid en je komt niet erg ver, niet erg snel.
• In het Instabiele Geval: Je zou denken dat het gesteente de chaotische vingers zou vertragen, maar dat deed het niet. Het chaotische "vingeren" was zo sterk dat het de neiging van het gesteente om dingen te verspreiden overmeesterde. De vingers sloegen door de barrières van het gesteente heen.
  • Het Resultaat: De menging werd efficiënter in het instabiele geval dan in het stabiele. De vingers maakten het interface tussen vloeistoffen smaller en scherper, waardoor ze sneller in elkaar oplossen dan wanneer het gesteente perfect glad was.

De "Snelheid versus Menging"-afweging

Het paper benadrukt een trek- en duwgevecht tussen hoe snel de vloeistof beweegt en hoe goed het mengt:

• Hoge Snelheid (Hoog Rayleigh-getal): Wanneer de vloeistof zeer dicht is en snel beweegt, neigt het om in een strakke stroom te blijven. In een uniform gesteente mengt het goed. Maar in een hobbelig gesteente wint het "barrière-effect". De vloeistof wordt geblokkeerd, beweegt sneller langs de gemakkelijke paden, maar mengt minder in totaal.
• Lage Snelheid (Laag Rayleigh-getal): Wanneer de vloeistof langzaam beweegt, doet diffusie (de natuurlijke neiging om zich uit te spreiden) het werk. Hier helpt de ongelijkmatige gesteente eigenlijk. Het vroege chaos veroorzaakt door de hobbel van het gesteente zorgt ervoor dat de vloeistof beter mengt dan in een glad, uniform gesteente.

De "Anisotropie"-factor: Richting Maakt Uit

De onderzoekers keken ook naar de richting van de gaten in het gesteente.

• Horizontale Lagen (Zoals een Taart): Als het gesteente horizontale lagen heeft van harde en zachte plekken, werkt het als een reeks planken. De zakkende vingers raken een plank en stoppen. Dit zet de menging snel stil.
• Verticale Lagen (Zoals een Stapel Papier): Als de lagen verticaal zijn, kunnen de vingers er gemakkelijk langs glijden, maar de hele stroom beweegt langzamer omdat het de verticale muren moet navigeren.

De Conclusie

Het paper concludeert dat de efficiëntie van menging afhangt van een delicate balans:

1. De Grootte van het Chaos: Hoe groot zijn de "vingers" van menging?
2. De Grootte van de Hobbel van het Gesteente: Hoe groot zijn de gaten in de spons?

Als de vingers klein zijn en de hobbel van het gesteente groot (hoge snelheid, hoge variantie), werkt het gesteente als een barrière, vertraagt het de menging en laat het de vloeistof verder reizen.
Als de vingers groot zijn en de hobbel van het gesteente klein (lage snelheid), helpen de hobbel van het gesteente eigenlijk om de menging op gang te brengen, waardoor het efficiënter wordt dan een glad gesteente zou zijn.

Kortom: De "sponsachtige" ondergrond van de natuur vertraagt niet alleen dingen; het verandert de regels van het spel. Soms blokkeert het de stroom, en soms, als de vloeistof chaotisch genoeg is, helpt het de vloeistoffen sneller met elkaar te mengen dan ze dat zouden doen in een perfect, gladde wereld.

Technische samenvatting

1. Probleemstelling

De studie onderzoekt de complexe wisselwerking tussen zwaartekrachtsstromen (vloeistoffen aangedreven door dichtheidsverschillen) en heterogeniteit van poreuze media. Waar eerdere onderzoek de zwaartekrachtsstromen in homogene media of de effecten van heterogeniteit op eenvoudige menging heeft onderzocht, richt dit werk specifiek op het gecombineerde effect van:

• Permeabiliteitsheterogeniteit: Weergegeven door log-normale stochastische velden met variërende varianties en correlatielengtes.
• Dichtheidsstratificatie: Vergelijking van stabiele scenario's (lineair dichtheids-concentratie-relatie, bijvoorbeeld zoutwaterintrusie) met instabiele scenario's (niet-monotoon dichtheidsrelatie, bijvoorbeeld CO₂-injectie in pekel waar vingerinstabiliteiten optreden).
• Oplossing en Menging: Het begrijpen hoe heterogeniteit de oplossingsflux, migratiesnelheid en spreiding van de stroming verandert, wat cruciaal is voor toepassingen zoals aquifer-sanering en geologische CO₂-opslag.

2. Methodologie

De auteurs hebben hoogwaardige numerieke simulaties toegepast om de besturingsvergelijkingen voor tweedimensionale, mengbare stroming in poreuze media op te lossen.

• Besturingsvergelijkingen: Het systeem wordt gemodelleerd met de Darcy-vergelijking voor stroming en de advectie-diffusievergelijking voor concentratietransport. De vergelijkingen zijn genormaliseerd met het Rayleigh-getal ($Ra$).
• Permeabiliteitsmodellering: Heterogeniteit wordt gemodelleerd als een log-normaal stochastisch proces. Belangrijke variabele parameters zijn:
  • Variantie ($\sigma^2_{\log k}$): 1, 3 en 5.
  • Correlatielengtes ($\lambda_x, \lambda_z$): Gebruikt om isotrope en anisotrope velden te creëren (horizontale versus verticale stratificatie).
  • Anisotropieverhouding ($\gamma = \lambda_x / \lambda_z$): Geteste waarden van 0,5, 1, 2 en 10.
• Dichtheidswetten:
  • Stabiel geval: Lineaire relatie ($\rho = -c$).
  • Instabiel geval: Niet-monotoon relatie die CO₂-pekelsystemen nabootst, met een intermediair dichtheidsmaximum dat vingerinstabiliteiten triggert.
• Numerieke Opstelling:
  • Oplosser: SECUReFOAM (OpenFOAM-pakket) met de finite volume-methode.
  • Rooster: Hoogresolutie gestructureerd rooster ($16.000 \times 800$ cellen) om dunne vingers en concentratiegradiënten op te lossen.
  • Metrieken: De studie kwantificeerde drijvende massa, scalair dissipatiepercentage ($\chi$), oplossingsflux, interface-lengte en interface-breedte.

3. Belangrijkste Bijdragen

Het artikel biedt een uitgebreide analyse van hoe heterogeniteit de fysica van zwaartekrachtsstromen verandert, en daagt specifiek de aanname uit dat heterogeniteit altijd menging verbetert. Belangrijke bijdragen zijn:

• Tegenwerkende Mechanismen: Aantonen dat in instabiele (vinger) regime's het dispersieve effect van heterogeniteit wordt gecompenseerd door de vingerinstabiliteit, wat leidt tot andere uitkomsten dan in stabiele regime's.
• Anisotropie-effecten: Identificatie van een "barrièreeffect" in horizontaal gestratificeerde media (hoge $\gamma$) dat verticale convectie onderdrukt en oplossing reduceert, in contrast met het gedrag in verticaal gestratificeerde media.
• Variantiedrempels: Vaststellen dat de variantie van het permeabiliteitsveld alleen de oplossings-efficiëntie significant verbetert ten opzichte van homogene gevallen wanneer het Rayleigh-getal laag is.
• Blob-vangst: Gedetailleerde beschrijving van het mechanisme waarbij drijvende vloeistof wordt vastgehouden in zones met lage permeabiliteit, waardoor "blobs" ontstaan die snel oplossen en lokale opwaartse vingering creëren, zelfs bij stabiele stratificatie.

4. Belangrijkste Resultaten

A. Homogeen versus Heterogeen Medium

• Stabiel geval: Heterogeniteit werkt voornamelijk als een dispersief agent. Het verbreedt het menginterface, vermindert de concentratiegradiënt en vermindert hierdoor de oplossing. Het verhoogt echter de migratiesnelheid van de neus van de stroming.
• Instabiel geval: Heterogeniteit versnelt het begin van convectie. Hoewel dit aanvankelijk de oplossing verhoogt, verkleint het dispersieve effect van de heterogeniteit uiteindelijk het interface in vergelijking met het stabiele geval, maar creëert het een complexe interactie waarbij het homogene geval vaak beter presteert dan het heterogene geval in totale opgeloste massa voor hoge $Ra$.

B. Impact van Anisotropie ($\gamma$)

• Horizontale Stratificatie ($\gamma > 1$): Creëert een barrièreeffect. Horizontale lagen met lage permeabiliteit vangen drijvende vloeistof, onderdrukken vingerformatie en sluiten convectie lokaal af. Dit leidt tot een reductie in totale oplossing in vergelijking met homogene media, ondanks een initiële piek in oplossingsflux.
• Verticale Stratificatie ($\gamma < 1$): Werkt als een barrière voor horizontale migratie maar bevordert verticale vingerdaling. Dit resulteert in een intermittent oplossingsgedrag terwijl de stroming permeabiliteitsbarrières kruist, met een vroeger begin van het "shutdown"-regime (waarbij ophoping van dichte vloeistof aan de onderkant menging stopt).

C. Impact van Variantie ($\sigma^2_{\log k}$)

• Hoge Variantie: Verhoogt het dispersieve effect, verbreedt het interface en vermindert de concentratiegradiënt die oplossing aandrijft.
• Efficiëntie-trade-off: Voor hoge Rayleigh-getallen (sterke convectie) lossen heterogene gevallen over het algemeen minder massa op dan homogene gevallen omdat de barrière- en dispersieve effecten domineren.
• Lage Rayleigh-getallen: Heterogeniteit kan de oplossing verbeteren ten opzichte van het homogene geval vanwege het vroege begin van convectie en interface-elongatie, aangezien diffusie een dominantere rol speelt ten opzichte van de instabiliteitsgrootte.

D. Interface-dynamiek

• In instabiele gevallen is het interface over het algemeen smaller dan in stabiele gevallen vanwege de dominantie van vingeren boven diffusie.
• De aanwezigheid van vastgehouden "blobs" in gebieden met lage permeabiliteit veroorzaakt oscillerend gedrag in interface-lengte en oplossingsflux, met name in anisotrope media.

5. Betekenis en Implicaties

Dit onderzoek is cruciaal voor het optimaliseren van Geologische CO₂-opslag (GCS) en aquifer-sanering:

• CO₂-opslag: De bevindingen suggereren dat in scenario's met hoge permeabiliteit en hoge convectie (hoge $Ra$), natuurlijke heterogeniteit de oplosbare opslag van CO₂ daadwerkelijk kan belemmeren in vergelijking met geïdealiseerde homogene modellen. Ingenieurs moeten rekening houden met het "barrièreeffect" van horizontale stratificatie, wat verticale menging kan beperken.
• Sanering: Voor DNAPL-sanering helpt het begrijpen dat heterogeniteit vloeistoffen kan vasthouden in "blobs" die snel oplossen bij het voorspellen van de levensduur van saneringsinspanningen.
• Modelleringsrichtlijnen: De studie benadrukt dat eenvoudige upscaling van heterogene media onvoldoende is. De interactie tussen de grootte van convectieve instabiliteiten en de correlatielengte van het permeabiliteitsveld is de bepalende factor voor oplossings-efficiëntie. Optimale oplossing vereist een specifiek evenwicht tussen het Rayleigh-getal en de sterkte van de heterogeniteit.

Concluderend stelt het artikel vast dat heterogeniteit niet universeel menging verbetert; in plaats daarvan is de impact conditioneel afhankelijk van de stabiliteit van de dichtheidsstratificatie en de relatieve schalen van de stromingsinstabiliteiten versus de geologische heterogeniteit.