Geomodelling For Reservoir Engineers 2020
Geomodelling For Reservoir Engineers 2020
Geomodelling For Reservoir Engineers 2020
1
11/5/2020
Online Training
❑ Certificates will be sent after the course to those who attend 75%
of the course or higher
❑ Will take a few days to process
❑ “Raw” unedited video will be sent 2 days after the course via e-mail
❑ Fully edited video will be available on the website
❑ Can take over a week to process
2
11/5/2020
Online Training
❑ Previous online training courses are available here:
https://www.cmgl.ca/training/online/previous/
Online Training
❑ Duration of the course:
❑ 9:00am to 4:00pm MDT
❑ 1 hour break ~12:00pm to 1:00pm MDT
❑ Two 15 minute breaks
3
11/5/2020
Online Training
❑ Online Agenda
❑ PowerPoint Presentation
❑ Tutorials
4
11/5/2020
Course Agenda
❑ Basic Concepts
▪ Geological modelling
▪ Geostatistics
▪ Upscaling
▪ Uncertainties in reservoir modelling
❑ Exercise
CMG Workflow
5
11/5/2020
Geological Modelling
Basic Concepts
What is Geomodelling?
❑ All mathematical methods used to produce a geological model of a subsurface
object
6
11/5/2020
What is Geomodelling?
Primary Objective:
▪ Predict the spatial variation of geological variables
Variable
▪ Properties of the geological subsurface that exhibit spatial variability and can be
measured in terms of numerical values
Spatial Variation
▪ When a quantity that is measured at different spatial locations exhibits values that differ
across the locations
7
11/5/2020
Background
❑ Modern reservoir characterization started due to:
▪ Deficiency of oil recovery techniques
❖ inadequate reservoir description
▪ Aim → predict inter-well property distributions (𝜑, K)
Static Data
Dynamic Data
8
11/5/2020
▪ In geological model construction, soft data can be used as trend or co-variance data to
decrease the geological uncertainty
9
11/5/2020
Seismic
Interpretation
Petrophysical Model
Petrophysical &
Well Correlation
Geostatistical Study
Outcrop/Analogue
❑ Outcrop is a visible exposure of bedrock on the surface of the earth
❑ Does not cover the majority of the Earth’s surface
▪ In most places the bedrock is covered by soil and vegetation
❑ Rock may be exposed due to erosion or tectonic uplift
❑ Outcrops are visually inspected and analyzed to get an idea of the subsurface
formation and bedding
10
11/5/2020
Outcrop/Analogue
❑ Analogue data is made up from nearby fields with similar rock
characteristics and distributions
Seismic Data
❑ Found by estimating the properties of the Earth's subsurface from seismic
testing
❑ Shock waves are sent into the earth and the data is generated by the sound
waves
❑ Sound waves reflect off the rock formations and are captured by geophones
on the surface
▪ A geophone is a device that converts ground movement (velocity) into voltage
▪ Deviation of the measured voltage from the base line is called the seismic
response
❑ Acoustic properties of the rock can help determine the type of rock and
geological structure
❑ Seismic Interpretation is the extraction of subsurface geologic information
from seismic data
11
11/5/2020
Seismic Data
❑ Aim of seismic interpretation is to produce structural maps that reflect
spatial variation in depth of certain geological layers
❑ Seismic data can be used to infer:
▪ Structure
❖ Faults
❖ Fractures
▪ Depositional Settings
❖ Channels
❖ Turbidite system
▪ Sand Continuity/Variation
▪ Bed Thickness
▪ Depositional Facies
▪ Reservoir Properties (Φ)
▪ Fluid Type Identification
12
11/5/2020
Structure Model
❑ Horizons picked on seismic interpretation and well data
▪ In absence of seismic data, possible to be based off well data only
❑ Only faults with sufficient throw that may impact reservoir performance continuity, fluid
distribution etc. should be incorporated
13
11/5/2020
14
11/5/2020
Stratigraphy Model
❑ Use info from well log and core correlation along with seismic depth conversion to determine
stratigraphy
❑ Subdivide reservoir in zones based on sub-surface horizons and seams
Geostatistical Study
❑ Statistical method accounting for spatial variability
❑ Goal:
▪ Utilize wide variety of data, in different scales and accuracies, to construct reservoir models
▪ Models are able to represent geological heterogeneities
▪ Quantifying uncertainties by producing many equiprobable models
15
11/5/2020
Facies Model
❑ Facies:
▪ A mass of rock that can be recognized by its composition, structures or fossil content
and mapped on the basis of those characteristics
▪ Within a given facies the reservoir properties can vary significantly
❖ Leads to a further subdivision known as Flow Units
Flow Units
▪ Regions in the sedimentary sequence that are judged to control the movement of
injected and produced fluids within the reservoir
▪ Log Facies:
▪ Separate the reservoir into Φ classes possibly linked to Vshl or other logs curves
❖ Calibrated to core properties or depositional facies
❖ This model may reflect production characteristics
▪ Depositional Facies:
▪ Rarely link to logs perfectly
▪ Dependant on core interpretations → arguable
▪ Too many facies = a complicated answer
16
11/5/2020
Petrophysical Model
❑ Petrophysical data needs to be upscaled to the geomodel
❑ Properties from log interpretation:
▪ Volume shale
▪ Effective Porosity
▪ NTG
▪ SW
❑ Log data should honor core analysis
❑ Geomodel should honor log data
❑ Fluid contacts (GOC, WOC) also need to be identified
Upscaling
❑ Geological characterizations typically contain on the order of 107 - 108 cells
▪ Impractical for flow simulation
▪ Long simulation run times
❑ Upscaling involves reducing the amount of cells by amalgamating the data
❑ Smaller sets of characteristics to represent the most significant aspect of the
reservoir
17
11/5/2020
Flow Simulation
❑ Grid block model of petroleum reservoirs where each block represents a
local part of the reservoir
❑ Within a grid block the properties are uniform
▪ Porosity, permeability, relative permeability, etc.
❑ Model incorporates:
▪ Reservoir fluids (PVT)
▪ Reservoir description (porosities, permeability, etc.)
▪ Property distribution in space
▪ Individual well pressure and flowrates with time
❑ Purpose is estimation of field performance (e.g., oil recovery) under one
or more producing schemes
18
11/5/2020
Geostatistics
Basic Concepts
Geostatistics
❑ What is Geostatistics?
▪ Data integration
19
11/5/2020
Geostatistics
❑ Why Geostatistics?
▪ Uncertainty analysis
Geostatistic Concepts
❑ Deterministic
▪ All data known beforehand
▪ No representation of uncertainty
▪ Always reproducible given same starting conditions
❑ Probabilistic (Stochastic)
▪ Element of chance involved
▪ Know likelihood something will happen, but don’t know when
▪ A random probability distribution or pattern that may be analyzed statistically, but may
not be predicted precisely
20
11/5/2020
❑ Based on sampling, we can still describe its behavior in terms of statistics and
probabilities
❑ Stochastic simulation:
▪ Generation of hypothetical data (realizations) with certain probabilities by
providing (pseudo) random input values
21
11/5/2020
22
11/5/2020
Skewed Normal
Geostatistic Terms
❑ Variance
▪ Numerical value indicating how widely individuals in a group vary
▪ The average of the squared differences from the mean
𝑁
1
σ² = (𝑥𝑖 − µ)²
𝑁
𝑖=1
N = number of samples
𝑥𝑖 = sample value
µ = mean
23
11/5/2020
Geostatistic Terms
❑ Standard Deviation
▪ Used to quantify amount of variation/dispersion for a set of data values
▪ Square root of the variance
𝑁
1
σ= σ² = (𝑥𝑖 − µ)²
𝑁
𝑖=1
Geostatistic Terms
❑ De-clustering
▪ Assignment of weights whereby redundant data are given less weight in the
construction of a representative probability distribution or the calculation of summary of
statistics
❖ Used to reduce possible bias in the data histogram
❖ Reduces the influence of over sampling in sweet spots (wells drill in good quality reservoir)
❖ Clustered data will have lesser weights than isolated data points
❖ If not used, the data are all equally weighted
24
11/5/2020
Geostatistic Terms
❑ Normal Scores Transformation
Variogram
▪ Is a function of the data variance against
the distance between the data locations
25
11/5/2020
Variogram (terms)
❑ Sill
▪ Semivariance value at which the variogram levels off
❑ Range
▪ Lag distance at which the semivariogram reaches the sill value
▪ No correlation of data beyond the range
❑ Lag
▪ Separation distance when comparing pairs of data value
❑ Nugget
▪ Semivariogram value at the origin
▪ Variability at distances smaller than typical sample spacing, including measurement error
▪ Noise level
Variogram
Sill
Semivariance
Range
Variogram Model
Nugget
Lag Distance
26
11/5/2020
Variogram – Lag
▪ Pairs of values with respect to point 55
113 99 90
103 for Lag 0:
97
112
101 92 55 62
106 91
55 53
86 Lag 2
100
104 99 74 55 78
67
96 Lag 1 111 55 98
62 Lag 0 105
109
63
107 97 82
53
55
102 ▪ Pairs for Lag 1:
109 98 94
78 88
55 99
77 55 86
95 102
68 69
87 55 67
76
78 55 74
79 84
77 55 …
Precipitation
27
11/5/2020
28
11/5/2020
Variogram Models
Isotropic Anisotropic
29
11/5/2020
❑ Simulation
▪ Use when variability observed in the data needs to be reproduced
▪ Can generate multiple results
▪ Important for uncertainty analysis
ESTIMATIONS
A. Inverse Distance Estimation
B. Ordinary Kriging Estimation
C. Ordinary Kriging with secondary variable
D. Trend Estimation
SIMULATIONS
E. Gaussian Geostatistical Simulation.
F. Gaussian Geostatistical with secondary variable
30
11/5/2020
Estimation - Kriging
❑ Weighted-average mapping technique that uses the variogram as the
weighting function
❑ Locally accurate and smooth (underestimates variability)
❑ Appropriate for visualizing trends
❑ Inappropriate for flow simulation where extreme values are important
❑ Does not assess for global uncertainty
▪ Produces one map –> Deterministic
❑ Good for Tops and Thickness
31
11/5/2020
Estimation - Kriging
❑ Basic Principles
▪ Search made around the area of the point being estimated
▪ Samples are assigned weights that reflect their spatial variability using the variogram
model
▪ A weighted average is calculated to estimate the value of the point
▪ Kriging minimizes the variance of the estimation error
▪ Provides the best average distribution
❖ High values will be underestimated
❖ Low values will be overestimated
❖ Average will be as predicted
❖ Honors data points, but does not reproduce global statistics
➢ Histogram, variance, and covariance
λ2=1/3 λ2=1/4
u u
λ3=1/3 λ3=1/2
λ1=1/3 λ1=1/4
Inverse Kriging
Distance
▪ ID gives same weight for ▪ Kriging gives shared weight to the two clustered
equal distances data points
32
11/5/2020
Geostatistical Simulation
❑ Gaussian Geostatistical Simulation
▪ Used to generate multiple versions (realizations) of grid values
▪ All realizations reproduce:
❖ Data
❖ Histogram
❖ Variogram of the data
▪ Appropriate for flow simulation
▪ Contrary to Kriging, adds an element of randomness
▪ Accounts for the uncertainty associated with some reservoir properties
▪ Small scale variability may mask large-scale trends
Geostatistical Simulation
❑ Gaussian Geostatistical Simulation
▪ Averaging multiple realizations will generate results very close to Kriging
▪ However, the average behavior of multiple realizations may be different
from the kriged model
33
11/5/2020
Estimation vs Simulation
34
11/5/2020
35
11/5/2020
Need to
Investigate
No Uncertainties? Yes
Deterministic Stochastic
Model Model
Upscaling
36
11/5/2020
Upscaling
❑ What is upscaling?
▪ Process in which a very heterogeneous region of reservoir rock
described with a large amount of fine grid cells is replaced by an
equivalent less heterogeneous region (coarse gird cells)
▪ Assign effective static properties to coarse scale cells
▪ Essentially an averaging procedure
▪ Coarse grid model should maintain the same storage and transport
properties of the reservoir rock of the fine scale model
Upscaling
❑ Why Upscale?
▪ Infeasible to run flow simulations directly on the geological model
▪ Reduce CPU time for flow simulation
▪ Make fine grid simulation practical
❖ Geological model ~ 10-100 million grid cells
▪ Reduce number of cells, but preserve accuracy in prediction
▪ Will be the starting point for prediction in the absence of historical
data
▪ Will be the starting point for history matching if historical data is
available
37
11/5/2020
Upscaling
Upscaling
❑ How QC and Verify:
▪ Visual QC
▪ Initial Volume check
❖ Pore Volume (POVO)
❖ Hydrocarbon Pore Volume (HC POVO)
▪ Spatial Check
❖ Histogram
❖ Mean
❖ Variance and Standard Deviation
▪ History Matching through flow simulation
38
11/5/2020
Upscaling
❑ Difficulties of upscaling:
▪ Bottleneck in the workflow
▪ Loss of details
❖ Dominated by loss of vertical resolution
▪ Lack of robustness
❖ Applicability to models with different global boundary conditions and well
locations
▪ How well does the upscaled model represent the fine scale
geological model?
❖ Won’t know fully until flow simulation is preformed
❖ Can be time consuming
39
11/5/2020
Uncertainty
❑ Includes:
▪ Measurement accuracy
▪ Incomplete or missing data
▪ Computation approximations
▪ Stochastic systems
Geophysical Uncertainties
❑ Acquisition
❑ Processing and interpretation
❑ Migration
❑ Time to depth conversion
❑ Horizon picking
❑ Fault position
❑ Well ties
40
11/5/2020
Geological Uncertainties
❑ Sedimentary depositional environments
❑ Rock type and associated heterogeneities
❑ Spatial distribution
❑ Sediment supply and particle size
❑ Geological markers
Petrophysical Uncertainties
❑ Net reservoir thickness
❑ Volume of shale
❑ Porosity
❑ Permeability
❑ Water saturation
❑ Fluid contact locations (GOC & WOC)
41
11/5/2020
Builder Results
Simulator
42
11/5/2020
43
11/5/2020
Maps
44
11/5/2020
Geostatistics in Builder
❑ Geostatistics on 3D Properties (Reservoir Menu)
▪ Object Modeling (Geological Bodies)
▪ Gaussian Geostatistical Simulation
45
11/5/2020
Property Statistics
Histogram based on
Color Scale intervals
Cross-Plots
❑ Scatter point plotting tool
❑ Property vs. property filtered plotting
46
11/5/2020
Cross-Plots
Ruler Tool
47
11/5/2020
Enter formula
Up-layering in Builder
❑ Options in CMG
▪ Can be done manually or automatically in Builder
❖ Note: Builder can only upscale in the vertical direction (up-layer)
▪ Manual:
❖ Builder provides mean and standard deviation of selected property
❖ Helps decide which layers to combine
❖ Guild property is typically porosity, permeability K, or combination
❖ Manually choose number of layers after upscale
❖ Manually choose which layers to combine
▪ Automatic:
❖ Builder can suggest optimal layers to combine given the desired number of layers and property
❖ Method based on “residue optimization” (SPE paper #57273)
❖ A table of residues is produced from analysis of all possible layering combinations
❖ Can determine number of simulation layers needed based on table
48
11/5/2020
Up-layering in Builder
❑ Can amalgamate layers
▪ Property variability calculations allow Builder to aid in up-layering process
▪ Builder can suggest the optimal combination of the layers
❖ Method based on “residue optimization”
❖ SPE paper #57273
49
11/5/2020
Rescue/RESQML Format
50
11/5/2020
Builder: Geostatistics
51
11/5/2020
Exercise
Geological Modelling using Builder
52