Les 3: Reservoirmodellering en Petrofysica

De Essentie van Reservoirmodellering

Reservoirmodellering is het proces van het construeren van een digitaal, driedimensionaal model van een ondergronds koolwaterstofreservoir. Dit model dient als een geïntegreerde representatie van alle beschikbare geologische, geofysische, petrofysische en engineering data. Het is geen doel op zich, maar een fundamenteel instrument dat ons in staat stelt om het reservoir te begrijpen, de hoeveelheid olie en gas te schatten (reserves), en de meest efficiënte strategie voor winning te plannen. We onderscheiden twee hoofdcategorieën modellen: statische en dynamische modellen.

• Statische (of geologische) modellen: Deze representeren het reservoir in zijn oorspronkelijke, ongestoorde staat, vóór de aanvang van productie. Ze beschrijven de geometrie (structuur, stratigrafie) en de verdeling van gesteente-eigenschappen (facies, porositeit, permeabiliteit, waterverzadiging). Dit model vormt de fundering voor alle verdere analyses.
• Dynamische (of simulatie) modellen: Deze bouwen voort op het statische model door er vloeistofeigenschappen en de fysica van vloeistofstroming aan toe te voegen. Met een dynamisch model kunnen we de beweging van olie, gas en water door het reservoir simuleren als gevolg van productie en injectie. Dit is essentieel voor het voorspellen van productieprofielen, het optimaliseren van putlocaties en het evalueren van verschillende winningsstrategieën (bv. waterinjectie, gasinjectie).

De Workflow van Statische Reservoirmodellering

Het bouwen van een robuust statisch reservoirmodel is een stapsgewijs, iteratief proces dat een multidisciplinaire aanpak vereist. De typische workflow omvat de volgende hoofdfasen:

1. Structurele Modellering

De eerste stap is het definiëren van het 'container' van het reservoir. Dit omvat de externe geometrie en de interne architectuur, voornamelijk bepaald door tektonische processen.

• Horizonten: De boven- en ondergrenzen van het reservoir, evenals belangrijke interne lagen, worden geïnterpreteerd op basis van seismische data. Deze reflectoren, die significante veranderingen in akoestische impedantie vertegenwoordigen, worden omgezet in dieptekaarten met behulp van een snelheidsmodel.
• Breuken (Faults): Breuken zijn cruciaal omdat ze het reservoir kunnen opdelen in afzonderlijke compartimenten, wat de vloeistofstroming aanzienlijk kan beïnvloeden. De interpretatie van breuken uit seismische data is vaak een uitdaging. Eenmaal gemodelleerd, moet de eigenschap van de breuk worden geanalyseerd: is deze 'sealing' (on-doorlatend) of 'leaking' (doorlatend)? Dit wordt geëvalueerd met technieken zoals 'Shale Gouge Ratio' (SGR) analyse, die de hoeveelheid klei in de breukzone kwantificeert.

Het resultaat van deze fase is een structureel raamwerk dat de 3D-ruimte definieert waarbinnen het reservoir zich bevindt.

2. Stratigrafische Modellering en Gridding

Nadat het structurele raamwerk is vastgesteld, wordt de interne gelaagdheid gedefinieerd. Dit volgt de principes van de sequentie-stratigrafie, waarbij het model wordt opgedeeld in genetisch gerelateerde eenheden. Vervolgens wordt dit volume onderverdeeld in een raster van cellen, een proces dat 'gridding' wordt genoemd.

• Grid Type: De keuze van het grid is belangrijk. Eenvoudige 'Cartesian grids' zijn rekenkundig efficiënt, maar kunnen complexe geometrieën, zoals hellende lagen en breuken, niet nauwkeurig weergeven. Daarom wordt in de praktijk vaak gebruik gemaakt van 'Corner-Point Grids', waarbij de cellen vervormd kunnen worden om de geologische structuren nauwkeurig te volgen.
• Grid Resolutie: De grootte van de cellen (bv. 50m x 50m lateraal, 1m verticaal) is een compromis. Een fijner grid kan geologische heterogeniteit beter vastleggen, maar leidt tot een zeer groot aantal cellen, wat de rekentijd voor simulaties drastisch verhoogt. Een te grof grid kan belangrijke details missen, zoals kleine kanalen of barrières voor stroming.

3. Faciesmodellering

Facies zijn gesteentelichamen met specifieke kenmerken (lithologie, textuur, sedimentaire structuren) die een bepaalde afzettingsomgeving weerspiegelen. Aangezien petrofysische eigenschappen sterk worden bepaald door de facies, is het modelleren van hun ruimtelijke verdeling van essentieel belang. Omdat we slechts op enkele locaties (de putten) directe informatie hebben, moeten we de verdeling tussen de putten in schatten. Hiervoor worden geostatistische methoden gebruikt.

• Pixel-gebaseerde methoden (bv. Sequential Indicator Simulation - SIS): Deze methode vult het grid cel voor cel in. Voor elke lege cel wordt de waarschijnlijkheid van elke facies berekend op basis van de data in nabijgelegen cellen en een statistisch model van ruimtelijke continuïteit (het variogram). Vervolgens wordt willekeurig een facies gekozen op basis van deze waarschijnlijkheidsverdeling. Dit proces wordt herhaald totdat alle cellen zijn gevuld. SIS is robuust en eert de putdata, maar kan soms geologisch onrealistische 'peper en zout'-patronen produceren.
• Object-gebaseerde methoden: Deze methoden proberen direct geologische objecten, zoals rivierkanalen of delta-lobben, in het model te plaatsen. De geometrie, oriëntatie en afmetingen van deze objecten worden gedefinieerd op basis van geologische analogieën en seismische data. Deze aanpak levert vaak geologisch meer realistische modellen op, vooral voor afzettingssystemen met duidelijke geometrische vormen.

4. Petrofysische Modellering

Nadat elke cel in het grid een facies heeft gekregen, moeten de continue petrofysische eigenschappen—porositeit, permeabiliteit en waterverzadiging—worden toegewezen.

• Modelleren van Porositeit: Porositeit is vaak relatief goed verdeeld en kan goed gemodelleerd worden met methoden als 'Sequential Gaussian Simulation' (SGS). Deze geostatistische techniek werkt vergelijkbaar met SIS maar dan voor continue variabelen. De modellering wordt vaak per facies uitgevoerd, omdat elke facies zijn eigen karakteristieke porositeitsverdeling heeft.
• Modelleren van Permeabiliteit: Permeabiliteit is veel moeilijker te modelleren. Het kan over korte afstanden ordes van grootte variëren en is sterk afhankelijk van de textuur en structuur van het gesteente. Een directe geostatistische modellering is mogelijk, maar vaak wordt permeabiliteit gemodelleerd door een relatie met porositeit te gebruiken (co-kriging of co-simulatie). De porositeit-permeabiliteit relatie is echter zelden eenduidig en vaak de grootste bron van onzekerheid.
• Modelleren van Waterverzadiging: Initiële waterverzadiging wordt doorgaans niet direct geostatistisch gemodelleerd. In plaats daarvan wordt het functioneel berekend op basis van de hoogte boven het vrije waterniveau, de porositeit en de capillaire druk-eigenschappen van het gesteente.

5. Upscaling en Onzekerheidsanalyse

Het uiteindelijke statische model, met zijn miljoenen of zelfs miljarden cellen, is te gedetailleerd voor dynamische simulatie. Het moet worden 'opgeschaald' (upscaled) naar een grover simulatiegrid. Tijdens dit proces worden de eigenschappen van de fijne cellen gemiddeld om een effectieve eigenschap voor de grovere simulatiecel te verkrijgen. Het correct middelen van permeabiliteit is hierbij de grootste uitdaging en vereist speciale technieken die de stroming correct representeren.

Een cruciaal aspect van reservoirmodellering is het erkennen en kwantificeren van onzekerheid. Er is onzekerheid in elke stap: de diepte van de horizonten, de positie van de breuken, de faciesverdeling, de porositeit-permeabiliteit relatie, etc. In plaats van één 'beste' model te bouwen (een deterministische aanpak), worden vaak honderden of duizenden mogelijke modellen (realisaties) gegenereerd, die allemaal consistent zijn met de beschikbare data. Door simulaties op al deze modellen uit te voeren, krijgen we een probabilistische voorspelling van de reserves en productie (bv. P10, P50, P90 schattingen), wat een veel robuustere basis biedt voor investeringsbeslissingen.

Introductie tot Petrofysica

Petrofysica is de wetenschap die de fysische en chemische eigenschappen van gesteenten en hun interacties met vloeistoffen bestudeert. Voor de petroleumgeologie is het de discipline die de opslagcapaciteit (porositeit) en de doorstromingscapaciteit (permeabiliteit) van een reservoirgesteente kwantificeert. Een gesteente kan pas een commercieel reservoir zijn als het zowel voldoende porositeit heeft om een significante hoeveelheid koolwaterstoffen op te slaan, als voldoende permeabiliteit om deze koolwaterstoffen met een economisch rendabele snelheid naar een productieput te laten stromen. We zullen de fundamentele eigenschappen en hun interacties in detail bespreken.

1. Porositeit (Φ)

Porositeit is de meest fundamentele reservoireigenschap. Het is de fractie van het totale gesteentevolume dat wordt ingenomen door poriënruimte, en wordt meestal uitgedrukt als een percentage of een fractie. De wiskundige definitie is:

Φ = (Volume van poriën / Totaal bulkvolume) = (V_bulk - V_grain) / V_bulk

Porositeit kan worden geclassificeerd op basis van haar oorsprong:

• Primaire Porositeit: Dit is de porositeit die ontstaat tijdens het afzetten van het sediment. Bij klastische gesteenten (zandstenen) is dit voornamelijk de ruimte tussen de korrels (intergranulaire porositeit). De kwaliteit ervan wordt bepaald door de sortering, rondheid en pakking van de korrels. Goed gesorteerde zanden hebben een hogere porositeit dan slecht gesorteerde zanden.
• Secundaire Porositeit: Dit is porositeit die ontstaat na de afzetting door geologische processen (diagenese). Voorbeelden zijn:
  • Oplossingsporositeit (dissolution): Mineralen (zoals veldspaten of carbonaatcement) worden opgelost door stromend formatiewater, waardoor nieuwe poriënruimte ontstaat. Dit kan de reservoirkwaliteit aanzienlijk verbeteren.
  • Breukporositeit (fracture): Tektonische krachten kunnen het gesteente doen breken, waardoor een netwerk van breuken ontstaat. Hoewel het volume van de breuken zelf vaak klein is, kunnen ze de permeabiliteit van een anderszins 'dicht' gesteente (zoals een kalksteen of schalie) drastisch verhogen.

Een cruciale distinctie wordt gemaakt tussen totale porositeit (alle poriënruimte) en effectieve porositeit (de onderling verbonden poriënruimte die bijdraagt aan vloeistofstroming). Geïsoleerde poriën (bv. binnenin een korrel of omgeven door cement) dragen bij aan de totale porositeit maar niet aan de effectieve porositeit. Voor reservoir-engineering is alleen de effectieve porositeit van belang.

2. Permeabiliteit (k)

Permeabiliteit is een maat voor het vermogen van een gesteente om vloeistoffen door te laten. Het is een eigenschap van het poreuze medium zelf, onafhankelijk van de vloeistof. De formele definitie komt voort uit de empirische wet van Henry Darcy (1856):

Q = - (k * A / μ) * (dP/dL)

Waarbij:

• Q = Stroomsnelheid (debiet) (m³/s)
• k = Permeabiliteit (m²)
• A = Dwarsdoorsnede-oppervlak (m²)
• μ = Dynamische viscositeit van de vloeistof (Pa·s)
• dP/dL = Drukgradiënt (Pa/m)
Open visualization in new tab

De eenheid voor permeabiliteit is de Darcy (D). Eén Darcy is de permeabiliteit die het mogelijk maakt dat een vloeistof met een viscositeit van 1 centiPoise (cP) met een snelheid van 1 cm/s stroomt door een oppervlak van 1 cm² onder een drukgradiënt van 1 atm/cm. In de praktijk zijn de meeste reservoirgesteenten veel minder permeabel, dus wordt de milliDarcy (mD) (1 D = 1000 mD) veel gebruikt.

De relatie tussen porositeit en permeabiliteit is complex. Een hoge porositeit garandeert geen hoge permeabiliteit. Een puimsteen heeft bijvoorbeeld een zeer hoge porositeit, maar de poriën zijn niet verbonden, dus de permeabiliteit is bijna nul. Permeabiliteit wordt niet alleen bepaald door het volume van de poriën, maar vooral door de grootte en de verbondenheid van de 'poriënkeeltjes' (pore throats). Kleimineralen in de poriën kunnen de permeabiliteit drastisch verlagen, zelfs bij een redelijke porositeit.

Wanneer er meerdere, niet-mengbare vloeistoffen (bv. olie en water) in de poriën aanwezig zijn, wordt het concept complexer:

• Absolute Permeabiliteit (k): De permeabiliteit gemeten wanneer het gesteente 100% verzadigd is met één enkele vloeistof.
• Effectieve Permeabiliteit (k_o, k_w, k_g): De permeabiliteit voor een specifieke vloeistof (olie, water of gas) wanneer er andere vloeistoffen aanwezig zijn. De effectieve permeabiliteit is altijd lager dan de absolute permeabiliteit.
• Relatieve Permeabiliteit (k_ro, k_rw, k_rg): De verhouding van de effectieve permeabiliteit tot de absolute permeabiliteit. Dit is een dimensieloze waarde tussen 0 en 1 die een functie is van de vloeistofverzadiging.
Open visualization in new tab

3. Vloeistofverzadiging (S)

Verzadiging (saturatie) is de fractie van de poriënruimte die wordt ingenomen door een specifieke vloeistof. De som van alle vloeistofverzadigingen is altijd gelijk aan 1 (of 100%).

S_w + S_o + S_g = 1

Waarbij S_w, S_o en S_g respectievelijk de verzadiging van water, olie en gas zijn.

Twee belangrijke concepten zijn:

• Irreduceerbare Waterverzadiging (S_wi): Zelfs in een olie- of gaszone is er altijd een bepaalde hoeveelheid water aanwezig die niet door de koolwaterstoffen kan worden verdrongen. Dit water wordt door capillaire krachten vastgehouden in de kleinste poriën en als een dun laagje op het oppervlak van de korrels (in water-natte systemen). Deze S_wi is niet produceerbaar.
• Residuele Olie/Gas Verzadiging (S_or / S_gr): Tijdens de productie, bijvoorbeeld door een water-sweep, zal niet alle olie worden verdrongen. Een deel blijft achter in de poriën, eveneens door capillaire krachten. Dit is de residuele verzadiging en bepaalt de uiteindelijke recovery factor.

4. Gesteente-Vloeistof Interacties

De verdeling en stroming van vloeistoffen in het reservoir worden niet alleen bepaald door de poriëngeometrie, maar ook door de interactie tussen de vloeistoffen en het gesteente-oppervlak.

• Wettability (Bevochtigbaarheid): Dit beschrijft de neiging van een vloeistof om zich te hechten aan een vast oppervlak in de aanwezigheid van andere, niet-mengbare vloeistoffen. Een reservoir kan water-nat (water-wet) zijn, waarbij water het korreloppervlak bedekt, of olie-nat (oil-wet), waarbij olie het oppervlak bedekt. Wettability heeft een diepgaande invloed op de capillaire druk, de relatieve permeabiliteit en de verdeling van S_wi en S_or. De meeste zandsteenreservoirs zijn oorspronkelijk water-nat, maar kunnen olie-nat worden door adsorptie van zware componenten uit de ruwe olie over geologische tijd.
• Capillaire Druk (Pc): Dit is het drukverschil over het grensvlak tussen twee niet-mengbare vloeistoffen in een porie. Het wordt veroorzaakt door oppervlaktespanning. De capillaire druk is omgekeerd evenredig met de straal van de poriënkeel. Dit betekent dat er een hogere druk nodig is om een niet-bevochtigende vloeistof (bv. olie in een water-nat systeem) in een kleinere porie te persen. Dit fenomeen verklaart het bestaan van een transitiezone in een reservoir, waar de waterverzadiging geleidelijk afneemt met de hoogte boven het vrije waterniveau, omdat de koolwaterstoffen eerst de grootste poriën binnendringen.

De Noodzaak van Meerdere Technieken

Het karakteriseren van een reservoir dat zich kilometers onder de grond bevindt, is een complexe puzzel. Er is geen enkele techniek die ons alle antwoorden kan geven. Een robuust begrip van het reservoir wordt verkregen door het integreren van data van verschillende bronnen, die elk op een andere schaal en met een andere meetmethode informatie verschaffen. De belangrijkste technieken zijn directe metingen aan boorkernen (core analysis), indirecte metingen via boorgatmetingen (well logging) en dynamische metingen door middel van puttesten (well testing).

1. Directe Metingen: Boorkernanalyse (Core Analysis)

Boorkernanalyse wordt beschouwd als de 'ground truth' of de waarheid op de grond. Het is de enige methode waarbij we daadwerkelijk een stuk van het reservoirgesteente naar het laboratorium kunnen brengen voor gedetailleerde metingen. Er wordt onderscheid gemaakt tussen Routine Core Analysis (RCAL) en Special Core Analysis (SCAL).

• Routine Core Analysis (RCAL): Dit omvat de standaardmetingen die op regelmatige intervallen (bv. elke 30 cm) op kleine 'pluggen' (cilindrische monsters) uit de kern worden uitgevoerd.
  • Porositeit: Gemeten met technieken zoals de Boyle's Law methode met helium.
  • Permeabiliteit: Gemeten met een gaspermeameter, waarbij stikstofgas onder een bekend drukverschil door de plug stroomt. Vaak wordt zowel horizontale als verticale permeabiliteit gemeten om anisotropie vast te stellen.
  • Vloeistofverzadiging: De hoeveelheid olie en water in de kern wordt bepaald met methoden zoals de Dean-Stark extractie.
  RCAL levert de basisdata voor de kalibratie van andere meetmethoden.
• Special Core Analysis (SCAL): Dit zijn meer complexe, tijdrovende en duurdere metingen die op een select aantal representatieve pluggen worden uitgevoerd. SCAL meet eigenschappen die de interactie tussen gesteente en vloeistoffen beschrijven, wat essentieel is voor dynamische modellering.
  • Elektrische Eigenschappen: Meting van de cementatie-exponent (m), saturatie-exponent (n) en de tortuositeitsfactor (a). Deze parameters zijn cruciaal voor de interpretatie van resistiviteitslogs (Archie's Wet).
  • Capillaire Druk: Gemeten met methoden als kwikinjectie, poreuze plaat of centrifuge. De resulterende curve relateert de waterverzadiging aan de hoogte boven het vrije waterniveau.
  • Relatieve Permeabiliteit: Gemeten door het gelijktijdig injecteren van twee vloeistoffen (bv. olie en water) door een kernmonster en het meten van drukvallen en verzadigingen.
  • Wettability: Bepaald met methoden zoals de Amott-test of de USBM-index.

Voordelen: Directe, nauwkeurige metingen. Nadelen: Data is alleen beschikbaar op de locaties waar kernen zijn genomen (zeer beperkte ruimtelijke dekking), het proces is duur en de kern kan tijdens het boren beschadigd raken, wat de metingen kan beïnvloeden.

2. Indirecte Metingen: Boorgatmetingen (Well Logging)

Na het boren van een put wordt een set meetinstrumenten ('tools') in het boorgat neergelaten om continu de fysische eigenschappen van de gesteentelagen te meten. Dit levert een 'log' op, een grafiek van een eigenschap versus diepte.

• Gamma Ray (GR) Log: Meet de natuurlijke radioactiviteit van het gesteente. Kleien (shales) bevatten radioactieve elementen (Kalium, Thorium, Uranium) en geven een hoge GR-waarde. Schone zandstenen en carbonaten hebben een lage GR-waarde. De GR-log is de primaire indicator voor lithologie en wordt gebruikt om het volume klei (Vshale) te berekenen.
• Porosity Logs: Er zijn drie hoofdtypen die samen worden gebruikt om porositeit te bepalen.
  • Density Log: Zendt gammastralen uit en meet de terugverstrooiing, die gerelateerd is aan de elektronendichtheid van de formatie. Hieruit kan de bulkdichtheid (ρ_bulk) worden berekend. Met kennis van de korreldichtheid (ρ_matrix) en de vloeistofdichtheid (ρ_fluid) kan de porositeit worden afgeleid.
  • Neutron Log: Zendt snelle neutronen uit die vertragen door botsingen, voornamelijk met waterstofatomen. De tool meet de populatie van vertraagde neutronen, wat een maat is voor de waterstofindex, en dus porositeit (aangenomen dat de poriën gevuld zijn met water of olie).
  • Sonic Log: Meet de reistijd (Δt) van een geluidsgolf door de formatie. In een compact gesteente reist geluid sneller (lage Δt) dan in een poreus gesteente. De Wyllie time-average-vergelijking wordt gebruikt om porositeit te berekenen.
• Resistivity Log: Meet de elektrische weerstand van de formatie. Koolwaterstoffen zijn zeer resistief, terwijl zout formatiewater zeer geleidend is (lage weerstand). Door de gemeten weerstand (Rt) te vergelijken met de verwachte weerstand als het gesteente 100% met water gevuld zou zijn, kunnen we de waterverzadiging (S_w) berekenen met behulp van Archie's Wet:

S_w = [(a * R_w) / (Φ^m * R_t)]^(1/n)

Deze vergelijking is de hoeksteen van kwantitatieve log-interpretatie. De parameters a, m, en n moeten worden gekalibreerd met SCAL-data.

Voordelen: Continue data over de gehele lengte van het boorgat, relatief goedkoop vergeleken met kernen. Nadelen: Indirecte metingen met een beperkte resolutie en onderzoeksdiepte (enkele tientallen centimeters rondom het boorgat), en de interpretatie is afhankelijk van vele aannames.

3. Dynamische Metingen: Puttesten (Well Testing)

Een puttest meet de drukrespons van het reservoir als reactie op een verandering in de productiesnelheid. Bij een 'drawdown'-test wordt de put met een constante snelheid geproduceerd en de bodemdruk gemeten. Bij een 'build-up'-test wordt de put gesloten en wordt het drukherstel gevolgd. De analyse van deze drukdata kan informatie opleveren over een veel groter deel van het reservoir dan kernen of logs.

• Gemiddelde Permeabiliteit: De analyse levert de effectieve permeabiliteit (k) op over het drainagegebied van de put, wat een schaal is van tientallen tot honderden meters. Dit is een veel representatievere waarde voor de stromingscapaciteit van het reservoir dan de metingen op kernpluggen.
• Huid-effect (Skin): Een maat voor de schade of verbetering van de permeabiliteit direct rondom het boorgat. Een positieve skin duidt op schade (bv. door boorspoeling), een negatieve skin op stimulatie (bv. door 'fracturing').
• Reservoirgrenzen: De vorm van de drukrespons kan de aanwezigheid van grenzen zoals breuken of een olie-water contact onthullen.

Voordelen: Meet dynamische eigenschappen op een representatieve, grote schaal. Nadelen: Duur om uit te voeren (productieverlies), levert een gemiddelde waarde op en verliest details over kleinschalige heterogeniteit.

Integratie: De Kracht van de Combinatie

Geen enkele techniek is superieur; ze zijn complementair. De sleutel tot een succesvolle reservoirkarakterisatie ligt in de integratie. Kern-data (RCAL en SCAL) worden gebruikt om de log-interpretatie te kalibreren (bv. de porositeit-permeabiliteit relatie en Archie's parameters). De gekalibreerde logs bieden vervolgens een continue, hoog-resolutie beschrijving van de eigenschappen langs de put. Deze put-data worden vervolgens gebruikt als 'harde data' in het geostatistische reservoirmodel. Ten slotte worden de dynamische data uit puttesten gebruikt om het upscaled reservoirmodel te valideren en te kalibreren (een proces dat 'history matching' wordt genoemd). Alleen door deze hiërarchische aanpak, waarbij de verschillende schalen van meting worden gecombineerd, kan een betrouwbaar en voorspellend reservoirmodel worden gebouwd.