Porous-Medium Scaling of CO$_2$ Plume Footprint Growth

Dit paper kwantificeert de groei van CO₂-pluimvoetafdrukken op veldschaal door analytische Barenblatt-oplossingen te vergelijken met seismische data van verschillende locaties, en biedt een fysiek transparant model voor de pluimdikte en -uitbreiding dat zowel injectie- als stilstandsscenario's beschrijft.

De kern

Stel je voor dat je een enorme, onzichtbare "luchtbel" van kooldioxide (CO2) in de grond injecteert, diep onder de zee of in een oude zoutmijn. Deze bel is niet statisch; hij beweegt, spreidt zich uit en verandert van vorm, net als een druppel inkt in een glas water. Maar in plaats van water hebben we hier te maken met poreus gesteente (zoals zandsteen) en zout water (brein).

Dit wetenschappelijke artikel, geschreven door Fernando Alonso-Marroquin en Christian Tantardini, probeert uit te leggen hoe snel en hoe groot deze CO2-bellen worden, en of we dat kunnen voorspellen met wiskunde.

Hier is de uitleg in simpele taal, met een paar creatieve vergelijkingen:

1. Het Probleem: Een onzichtbare olievlek

Wanneer bedrijven CO2 injecteren om het klimaat te helpen (door het onder de grond op te slaan), willen ze precies weten hoe ver die CO2-beweging gaat. Ze kijken naar foto's gemaakt door seismische scanners (een soort röntgenfoto's van de aarde). Maar deze foto's zijn vaak lastig te interpreteren. Hoe groot is de vlek precies? Hoe snel groeit hij?

De auteurs zeggen: "Laten we niet alleen kijken naar de foto's, maar laten we een wiskundig model gebruiken dat beschrijft hoe vloeistoffen zich door een spons bewegen."

2. De Wiskundige Sleutel: De "Poreuze Middel" Vergelijking

Stel je voor dat je een theedoek (de grond) hebt en je giet er een beetje siroop (de CO2) op. De siroop verspreidt zich niet als een perfecte cirkel die steeds sneller groeit. Omdat de doek vezels heeft, wordt de siroop langzamer en langzamer.

De auteurs gebruiken een wiskundige formule die ze de Barenblatt-oplossing noemen.

• De analogie: Denk aan een sneeuwpop die smelt. Als hij smelt, wordt hij niet alleen kleiner, hij verandert ook van vorm. De wiskunde in dit artikel beschrijft precies die vormverandering.
• Het mooie aan deze formule is dat hij voorspelt dat de rand van de CO2-bel een scherpe rand heeft. Het is niet zo dat de CO2 heel langzaam verdwijnt; er is een duidelijke lijn waar de CO2 stopt en het zoute water begint. Dit noemen ze "compact support" (een compacte steun).

3. Twee Verschillende Scenarios: De Kraan open of dicht?

Het artikel onderscheidt twee situaties, en dat is heel belangrijk:

Scenario A: De kraan staat open (Injectie loopt door)
Stel je voor dat je een tuinslang vasthoudt en water op de grond spuit.

• Direct onder de slang wordt de grond volledig verzadigd met water (de "kern").
• Verder weg wordt het water dunner.
• De les: Zolang je blijft spuiten, blijft er een dikke kern van CO2 onder de injectieput zitten. De buitenrand groeit, maar de binnenkant blijft "vol". De snelheid waarmee de buitenrand groeit, volgt een heel specifiek patroon (de wortel van de tijd).

Scenario B: De kraan gaat dicht (Injectie stopt)
Nu draai je de kraan dicht. Je hebt een bepaalde hoeveelheid water op de grond gespoten en nu moet het alleen maar verder verspreiden.

• De dikke kern onder de put begint te krimpen. Het water (of CO2) stroomt naar buiten om de druk te egaliseren.
• Op een gegeven moment is die dikke kern helemaal weg. De hele CO2-bel is nu een dunne, uitlopende vlek.
• De les: Zodra de kraan dicht is, verandert de groei-snelheid. De bel groeit dan trager dan tijdens het injecteren, en volgt de "Barenblatt-regel" voor langzame diffusie.

4. De Realiteit: Wat zeggen de foto's?

De auteurs hebben gekeken naar echte data van drie grote plekken waar CO2 wordt opgeslagen:

1. Sleipner (Noordzee, Nederland/Norwegen)
2. Aquistore (Canada)
3. Weyburn (Canada)

Ze hebben de foto's van deze plekken "uitgeprint" en de vlekken gemeten. Vervolgens hebben ze gekeken of de groeisnelheid van die vlekken overeenkwam met hun wiskundige voorspelling.

Het resultaat:
Het klopt! De CO2-bellen op deze plekken groeien precies zoals de wiskunde voorspelt voor een "poreus medium".

• De snelheid waarmee de vlek groeit, ligt ergens tussen de waarden die de formule voorspelt.
• Dit betekent dat hun simpele wiskundige model een heel goede manier is om te voorspellen hoe CO2 zich in de toekomst zal gedragen, zonder dat we duurdere en complexere computersimulaties hoeven te draaien voor elke kleine berekening.

5. Waarom is dit belangrijk?

Stel je voor dat je een verzekering wilt afsluiten voor een CO2-opslagproject. Je wilt weten: "Zal die CO2-bel over 50 jaar nog steeds veilig onder de grond zitten, of is hij ontsnapt?"

Dit artikel geeft ons een simpele, betrouwbare kompasnaald.

• Het laat zien dat we de groei van deze CO2-bellen kunnen begrijpen met een simpele wet: hoe dikker de grondlaag, hoe langzamer de verspreiding, en hoe de vorm verandert als we stoppen met injecteren.
• Het helpt ingenieurs om te zeggen: "Als we stoppen met injecteren, zal de kern van de CO2-bel langzaam verdwijnen en de rest van de vlek zal zich op een voorspelbare, veilige manier verspreiden."

Samenvattend:
De auteurs hebben een wiskundige "recept" gevonden voor hoe CO2 zich door de aarde verspreidt. Ze hebben bewezen dat dit recept werkt door het te vergelijken met echte foto's van ondergrondse CO2-bellen. Het is als het vinden van de perfecte formule om te voorspellen hoe een vlek inkt zich verspreidt in een spons, maar dan voor het redden van het klimaat.

Technische samenvatting

Titel: Schaling van de CO2-pluimvoetafdruk in poreuze media

Auteurs: Fernando Alonso-Marroquin en Christian Tantardini
Instituut: Center for Integrative Petroleum Research, King Fahd University of Petroleum and Minerals, Saoedi-Arabië

---

1. Probleemstelling

Bij geologische opslag van CO2 is het begrijpen van de migratie en uitbreiding van de CO2-pluim in ondergrondse porieuze formaties cruciaal voor monitoring en risicobeheer. Traditionele modellen voor grote schaal baseren zich vaak op Darcy-stroming met verticale middeling of scherpe-interface-aannames. Hoewel deze nuttig zijn, missen ze soms de nuance van niet-lineaire diffusie-effecten die voortvloeien uit de compressibiliteit van het gas en de interactie met het poreuze medium.

De kernvraag is hoe de "voetafdruk" (footprint) van een CO2-pluim in de tijd groeit en of deze groei kan worden beschreven door schaalwetten die voortkomen uit niet-lineaire diffusie-theorieën (specifiek de Porous Medium Equation of PME). Er is behoefte aan een fysiek transparant kader om velddata (zoals seismische monitoring) te vergelijken met analytische oplossingen, rekening houdend met zowel injectieperiodes als periodes waarin de injectie stopt (shut-in).

2. Methodologie

De auteurs hanteren een gelaagde aanpak die begint bij fundamentele transportvergelijkingen en eindigt bij een vergelijking met velddata:

• Van GBL naar q-PME: De studie start met de Global Buckley-Leverett (GBL) transportvergelijking voor meercomponenten- en meerfasenstroom. Door aannames te doen over verticale segregatie, advektie-dominantie en het verwaarlozen van adsorptie en diffusie, wordt dit gereduceerd tot een scalaire niet-lineaire diffusievergelijking.
• De q-Porous Medium Equation (q-PME): Er wordt een constitutieve sluiting geïntroduceerd waarbij de diffusiviteit afhangt van de dichtheid (of dikte) volgens een machtswet: $D(\rho) \propto \rho^{1-q}$. Dit leidt tot de q-PME:
  $$\frac{\partial \rho}{\partial t} = D_0 \nabla \cdot (\rho^{1-q} \nabla \rho)$$
  Hierbij bepaalt de exponent $q$ het diffusie-regime. Voor $0 < q < 1$ (trage diffusie) ontstaan oplossingen met een eindige voortplantingssnelheid en een compacte ondersteuning (de pluim heeft een scherpe rand).
• Analytische Oplossingen (Barenblatt): De auteurs gebruiken zelf-achtige (self-similar) Barenblatt-oplossingen voor de q-PME. Deze oplossingen beschrijven de vorm van de pluim en hoe de straal $R(t)$ schaalt met de tijd ($R(t) \propto t^\beta$).
• Gecompositeerd Profiel met een Kern: Een belangrijke innovatie is het introduceren van een model dat een "transiënte volledige-dikte kern" ($a(t)$) combineert met een Barenblatt-staart.
  • Tijdens injectie: Dicht bij de put kan de CO2 de volledige aquiferdikte $H$ innemen ($b=H$).
  • Na shut-in: Deze kern krimpt en verdwijnt uiteindelijk, waarna de pluim overgaat in een zuivere Barenblatt-regime.
• Veldvalidatie: De theoretische voorspellingen worden getoetst aan tijdreeksdata van drie grote opslaglocaties: Sleipner (Noordzee), Aquistore en Weyburn-Midale (Canada). De auteurs digitaliseren gepubliceerde seismische platen om een equivalente straal $R_{eq}(t)$ te extraheren en passen een machtwet ($R \propto t^\beta$) toe.

3. Belangrijkste Bijdragen

1. Analytisch Kader voor Pluimgroei: Het artikel biedt een gesloten analytische oplossing voor de groei van CO2-pluimen die de overgang beschrijft tussen injectie-gedreven groei en herverdeling na shut-in.
2. Transiënte Kernmodel: De ontwikkeling van een model dat een volledige-dikte kern ($a(t)$) en een Barenblatt-staart combineert. Dit verklaart waarom pluimen tijdens injectie anders gedragen dan na het stoppen van de injectie.
3. Koppeling van Dichtheid en Dikte: Het toont aan dat de verticale gemiddelde dikte van de CO2-laag dezelfde schalingswetten volgt als de dichtheid in de q-PME, mits een scherpe-interface-aanname wordt gebruikt.
4. Validatie met Velddata: Het is een van de eerste studies die directe kwantitatieve vergelijkingen maakt tussen de exponenten van de Barenblatt-oplossingen en daadwerkelijk gemeten pluimgroei uit seismische monitoring.

4. Resultaten

• Schaalwetten:
  • Voor constante injectie groeit de buitenrand van de pluim volgens een wortelwet ($R(t) \propto t^{1/2}$), zolang er een volledige-dikte kern aanwezig is.
  • Na shut-in krimpt de kern ($a(t)$) en verdwijnt deze op een kritieke tijd $t_c$. Daarna volgt de pluim de trage niet-lineaire diffusiewet: $R(t) \propto t^\beta$, waarbij $\beta = 1/[d(1-q) + 2]$. Voor een axiale symmetrie ($d=2$) en $0 < q < 1$ geldt $0.25 < \beta < 0.5$.
• Veldvergelijking:
  • De uit de velddata gefitte exponenten zijn:
    • Sleipner: $\beta \approx 0.52$ (ligt iets boven het theoretische bereik, mogelijk door migratie vanuit diepere lagen).
    • Aquistore: $\beta \approx 0.37$ (binnen het theoretische bereik).
    • Weyburn: $\beta \approx 0.46$ (binnen het theoretische bereik).
  • Deze waarden ondersteunen het idee dat de groei van de pluimvoetafdruk in het algemeen consistent is met trage niet-lineaire diffusie, hoewel lokale geologische factoren en injectiegeschiedenis de exacte exponent beïnvloeden.
• Kern-dynamiek: Het model voorspelt dat de volledige-dikte kern na shut-in lineair krimpt in het kwadraat van de straal en uiteindelijk verdwijnt, waarna de pluim de zuivere Barenblatt-vorm aanneemt.

5. Betekenis en Conclusie

Dit werk biedt een fysisch transparante basislijn voor het interpreteren van CO2-pluimgroei op veldschaal. Het bewijst dat complexe, niet-lineaire diffusieprocessen in porieuze media effectief kunnen worden gemodelleerd met de q-PME en Barenblatt-oplossingen.

• Praktische Toepassing: Het model stelt ingenieurs in staat om de interne structuur van een pluim (kern vs. staart) en de groeisnelheid te voorspellen op basis van de injectiegeschiedenis en de compressibiliteit van het gas.
• Toekomstperspectief: De auteurs suggereren dat dit kader kan worden uitgebreid met fractionele afgeleiden om niet-lokale effecten te modelleren, wat essentieel kan zijn voor het verklaren van afwijkende schaalgedrag in complexe geologische formaties.
• Validatie: De overeenkomst tussen de theoretische exponenten en de gemonitorde velddata (Sleipner, Aquistore, Weyburn) versterkt het vertrouwen in het gebruik van vereenvoudigde niet-lineaire diffusiemodellen voor het validatie van grootschalige CO2-opslagprojecten.

Kortom, het artikel sluit de kloof tussen fundamentele wiskundige theorie (niet-lineaire diffusie) en praktische monitoring van geologische CO2-opslag.