Interpretation
Interpretation
Interpretation
Training on
Seismic data interpretation
by
Drilling
(Hydrocarbon and Data)
(log, VSP, Model, test data)
Training on Seismic data interpretation (3)
Data interpretation
D
GEOPHYSICAL
A
T (Seismic and Non-seismic)
A
GEOLOGICAL
T
Y
P
E PETROPHYSICAL
Software
• Interpretation/Processing software
– Landmark Graphic Corporation
• Interpretation, Interpretive processing and mapping
– Geoquest, Schlumberger
• Interpretation, Interpretive processing and mapping
– Paradigm
• 3D visualization, Interpretation, Interpretive processing,
Mapping and special studies
– Hampson Russell/Jason
• Special studies, AVO and Inversion
– GX Technology corporation
• Geophysical modelling
– Miscellaneous
• Printing/plotting, ZEH
Training on Seismic data interpretation (5)
Objective
• Prospect generation and identification of suitable locations
for drilling by interpreting subsurface geo-data
– Petroleum system
• Reservoir rock
– Porous and permeable sandstone, limestone
– Any other rock type forming trap, e.g., fractured shales
• Source rock
– Shale/carbonates
• Cap rock
– Shale/carbonates
• Trap
• Migration
• Reservoir characterization
– Estimation of reservoir parameters
• Area, thickness, porosity, saturation etc.
• The primary goal of seismic interpretation is to make maps that
provide geologic information (reservoir depth structure,
thickness, porosity, etc.).
Training on Seismic data interpretation (7)
Pre-requisite of interpretation
• Knowledge of geology and geophysical processes
– Type of data and type of information which can be extracted
• Poststack / Prestack, Time (PSTM), Depth (PSDM) Domain
• Objective of interpretation and amenability of seismic data
• Elements of seismic trace data
– Modes of display
– Amplitude, time and frequency
– Role of colours and colour bar
• Contrasting (map) and gradational (sections) colour scheme
• Polarity and phase conventions
• Resolution and detectability
• Seismic to well tie
– Character based matching between synthetic and seismic
– Depth to time conversion (T-D curves)
• Check shots, VSP, Synthetics
Training on Seismic data interpretation (8)
1
2
3
4
Marine reflection seismic acquisition system
Seismic waves are reflected back from subsurface interfaces
(boundaries) having impedance contrast. Recorded and processed
seismic data reveals the subsurface geologic model (structure and
other intrinsic properties) on interpretation.
Training on Seismic data interpretation (12)
Polarity Convention
Maximum
Reflection from top of wedge
Interference
Amplitude decreases
No interference
l 100
Wavelet 40 hertz
1 wavelength
Time
Resolution of the reflections from the top and
bottom of a bed is dependent on the interaction of
closely spaced wavelets
Training on Seismic data interpretation (18)
Limit of Separability
Age of rocks Very young medium old Very old
young
Depth Very shallow medium deep Very deep
shallow
Velocity 1600 2000 3500 5000 6000
(m/s)
Predominant 70 50 35 25 20
frequency
Wavelength 23 40 100 200 300
Separability 6 10 25 50 75
Training on Seismic data interpretation (19)
Limit of visibility
• Factors affecting the
visibility Limit of visibility
– Impedance contrast of S/N Example Limit
the geologic layer of
interest relative to the
Poor Water sand poor λ/8
data
embedding material
– Random and
Moderate Wayer or oil sand λ/12
fairly good data
systematic noise in the
data
– Phase of the data or
High Gas sand good λ/20
data
shape of seismic
wavelet Outstanding Gas sand λ/30
– It may be less than 1 excellent data
m to more than 40 m
Training on Seismic data interpretation (20)
Horizontal resolution
On an un-migrated section, horizontal resolution is determined
by the size of the Fresnel zone. The magnitude of Fresnel zones
can be approximated from the relationship
where
v t
rf rf = radius of the Fresnel zone.
v = average velocity.
2 f
t = two-way time in seconds.
f = dominant frequency in hertz.
Cap rock
oil
Petroleum System
Petroleum (23)
Training on Seismic data interpretation System
DEPOSITIONAL ENVIRONMENTS
Glacial
Alluvial Lacustrine
Deltaic Deep water
marine
Aeolian
Beach
Continental: Alluvial, Aeolian, Fluvial, Lacustrine
Transitional: Deltaic, Tidal, Lagoonal, Beach
Marine: Shallow water marine, Deep water marine Reef
Training on Seismic data interpretation (28)
Braided channel
Meander channel
Training on Seismic data interpretation (29)
Map
projection,
Name measurement,
datum
information
and bounding
CRS coordinates
are required
for creation of
project data
base.
Map projection
Training on Seismic data interpretation (34)
Training on Seismic data interpretation (35)
STRUCTURAL INTERPRETATION
Elements of display
Elements of display
Time slice
cross-line In-line
3-D Cube of
seismic data
Training on Seismic data interpretation (39)
Elements of display
Time 2390
Successive time
slices depict the
anticlinal closures
Time 2430
Vertical section
and time slices
Training on Seismic data interpretation (40)
Black
traces
DT RHOB IMP GR LLD RC
Seismic
Red
traces
synthetic
Overlay of synthetic
and logs on seismic
sections. Here synthetic
is shown by yellow trace
Fault correlation
We expect to detect fault from alignment of event termination.
Alignment of termination can be seen from vertical sections,
horizontal sections and slices of attributes.
Map-based fault interpretation offers new possibilities
for establishing fault trends. A horizon is displayed in
Map View and fault segments are picked where the
horizon’s z values change abruptly. The map-based fault
segment is a continuous line marking the fault’s
intersection with that horizon. These segments will
appear as fault intersection points on the vertical
seismic sections and can guide to interpret fault cuts and
correlate segments from line to line.
Training on Seismic data interpretation (52)
Normal faults
Upthrown block
Horizon is displayed in map view. Major faults can be picked from the
abrupt change in values
Training on Seismic data interpretation (57)
Strikeslip faults
B B’
Seismic section showing strikeslip faults. Change in seismic
facies across the vertical fault is distinctly seen on the section
Training on Seismic data interpretation (58)
B’
Map generation
– Base map generation
• Depicting seismic
• Well location
• Cultural data
• Block boundaries
• Scale, coordinates and legends
– Griding
• Various method with faults and without faults
– Incorporating heaves and throw
– Contouring
– Over lay of attributes
– Map analysis
Training on Seismic data interpretation (60)
Volume LKM
Umrat 97,000
M-Tapti 12,000
Umrat
S-Tapti 24,000
M-Tapti C25 2,900
C22 34,700
C37 38,000
CA-CD 9,300
S-Tapti C23
C37 C23 85000
CA-CD
Location Map
Training on Seismic data interpretation (61)
Seismic attributes
Reflection seismic data helps to recognize and characterize
stratigraphic entities in two ways-first via seismic responses, and
second via their intrinsic seismic properties (e.g., reflectivity,
velocity, acoustic impedance). Both have limitations so the best
approach is to use both of them, if possible.
Seismic attributes are important for stratigraphic interpretation and
reservoir characterization.
Training on Seismic data interpretation (65)
Seismic attributes
• Attributes are derivatives of basic seismic
measurements/Information
– Seismic attributes extract information from seismic data that is
otherwise hidden in the data
– These information can be used for predicting, characterizing,
and monitoring hydrocarbon reservoirs
• Basic information
– Time
– Amplitude
– Frequency
– Attenuation
– Phase
• Most attributes are derived from normal stacked and
migrated data volume
• Can be derived from Pre-stack data (AVO)
Training on Seismic data interpretation (66)
Seismic attributes
Attribute Information
Time-derived Structural information
(After A. R. Brown)
Training on Seismic data interpretation (68)
Seismic attributes
• Post-stack attributes can be extracted along one horizon
or summed over window
– Window attributes provide the concept of formation attributes
– Computation window may be constant flat time interval (stat slice)
– Constant time hung from one structurally-independent horizon
depicting proper reservoir interval
– Window between two structurally-independent horizons (e.g. , top
and base of reservoir)
– Within window gross or selection attributes
– Hybrid attributes :combination of amplitude and frequency
Training on Seismic data interpretation (69)
Real amplitude
Instantaneous amplitude
Instantaneous frequency
Instantaneous phase
dt/dx is the dip in the x direction, and dt/dy is the dip in the y
direction, with x and y as real-world coordinates.
TWT map showing polygonal faults. With help of colour and light
auto correlated horizons depict all the details
Training on Seismic data interpretation (86)
residual Map
Training on Seismic data interpretation (88)
Dip Map
Training on Seismic data interpretation (89)
Dip-azimuth Map
Training on Seismic data interpretation (90)
Coherence slice
Training on Seismic data interpretation (93)
Coherence slice
Training on Seismic data interpretation (94)
Coherency
slice
at 1500 ms
Training on Seismic data interpretation (95)
Top of Carbonate
build up
• Discontinuity volumes
Carbonate buildup
coherency
Carbonate build up
boundary
Stratigrahic features, boundaries and faults are better depicted by
coherency/discontinuity volumes derived from seismic volumes. In this
example carbonate buildup boundary is depicted by horizon slice from
coherency volume.
Training on Seismic data interpretation (96)
Frequency-derived attributes
• Instantaneous frequency
– Stratigraphic and reservoir studies
• Geologic frequency is simply the number of reflecting layers
within unit thickness of rocks
• Geologic layering is expressed in seismic frequency
measurements
• More layering may be indicative of more sands
• Dominant frequencies
– First dominant – properties of wavelet
– Second dominant- geologic frequency-most useful
Training on Seismic data interpretation (100)
Hybrid attributes
• Seismic character is mixture of amplitude
and frequency
• Hybrid attributes are good descriptor of
seismic character
– Wave shape
• Seismic facies (Depositional facies-channel systems)
– Arc length
• Total length of seismic trace over a time window
• Shows depositional facies
Training on Seismic data interpretation (101)
Seismic facies
map showing
differences in
lithofacies at
well C-37-1(gas)
and C-37-2(dry).
The Analysis
interval is -24 to
+24 ms with
reference to C-
37-1 Mahuva
pay.
Training on Seismic data interpretation (102)
C-37-2
Seismic facies
map with
Neural
Network
technique of
a channel
showing
facies within
and on the
banks of
channel
Training on Seismic data interpretation (104)
Mapping of
channel through
3-D visualization
techniques.
Opacity
rendering of
channel in the
western portion
of the area.
Selected range of
voxels was made
opaque and
others were kept
transparent.
Training on Seismic data interpretation (105)
B3500 Channel
Pay anomaly
Sweetness:
Amplitude/ sqr root of
frequency
Effective in sandstone
reservoirs
Training on Seismic data interpretation (106)
Sand
dispersal
trends
• Channel
Horizon slice mapping
through
horizon slice
Flattened section
Training on Seismic data interpretation (109)
Channel
l
Channels Coherency
brought out
by horizon
and
coherency
slice
Horizon slice
Training on Seismic data interpretation (110)
Thickness (ms)
DHI or HCI
• Bright spot
– Water sand has lower
DHI
impedance than
embedding medium
and impedance of gas
sand is further
reduced.
– Top and base
reflections show
natural pairing Peak on
– If sand is thick synthetic
enough, flat spot or seismogram
fluid contact
reflection should be
visible between gas
sand and water sand
– Flat spot will have
opposite polarity than
bright spot at top
– More common in
shallower sandstone
reservoirs ( Mio-Plio)
Training on Seismic data interpretation (118)
DHI
Low impedance
reservoir embedded
in high impedance
High impedance
reservoir
embedded in Low
impedance
Hydrocarbon fluid
effects in seismic data
Training on Seismic data interpretation (122)
DHI
• Frequency
– Gas reservoir
attenuate high
frequencies more than
do rocks without gas
saturation
– Low instantaneous
frequency below a
suspected reservoir
can be a good
indicator of gas
– It is rather unreliable
indicator as frequency
lowering can often be
caused by broadening
of gas reservoir
reflections
Training on Seismic data interpretation (123)
Pay-1
Pay-3
(2537.5 m) Pay-4
Low
impedance Peak on
synthetic
seismogram
Log and seismic response of the Pay. Pay sands are characterized by low impedance
generating strong peak (blue) on normal polarity seismic section
Training on Seismic data interpretation (129)
Validation of DHI
• AVO
– In many practical cases gas sands show an
increase of amplitude with offset
– Many difficulties of theoretical and practical
nature
• Data is pre-stack hence lower S/N
• Multidimensional, many mode of display
Training on Seismic data interpretation (130)
Validation of DHI
• Shear wave
– Shear wave amplitude in conjunction with P-wave
amplitude can validate bright spot
– For “true” bright spot P-wave show high
amplitude but correlative S-wave amplitude is not
high
– For “false” bright spot P- and S-wave both have
high amplitude
– P-wave dim spot would correlate on an S-wave
higher amplitude reflection
Training on Seismic data interpretation (131)
Schematic illustration of porosity-depth trends for sands and shales. Both the sand and shale trends
can vary significantly due to composition, texture, pore fluids, temperature, and pressure gradients.
Hence, no attempt is done to assign absolute scales. However, there are a few rules of thumb: (1) The
depositional porosity of shales is normally higher than that of sands. (2) The porosity gradient with
depth is steeper for shales than for sands during mechanical compaction (i.e. at shallow depths). (3)
The porosity gradient with depth will be steeper for sands than for shales during chemical compaction
(i.e. quartz cementation of sands normally occurs at greater burial depth,beyond 2-3 km).
Training on Seismic data interpretation (132)
Occurrence of DHI
• Nature of DHI (bright, Dim
or Phase change) depends on
relative impedance of
hydrocarbon sand, water
Pitfalls in DHI
Exploration prospects based on a sound geologic model
and supported by seismic amplitude anomalies are highly
prospective and are usually assigned a high probability of
success. However, a fraction of such prospects, perhaps
10-30%, result in dry holes. Postdrill appraisal can usually
assign these results to one or more of the following
factors:
• Unusually strong lithologic variations
• Fizz water and low gas saturation
• Superposition of seismic reflections, and tuning effects
• Contamination of the seismic signal by multiples or other
undesired energy
Training on Seismic data interpretation (135)
G L I Sandstones
R L M Well H drilled here (trough (red)
D P
M1 Daman Top
M2 Daman1
M3
High amplitude Daman2
shales
Well-C (peak (blue) Well-D
High
impedance
sand
Impedance
log
In the area high amplitudes are deceptive. The drilled well confirmed thick
sandstones (16 m) but no hydrocarbon accumulation. The high amplitude
was caused by high impedance of water bearing sandstones
Training on Seismic data interpretation (137)
Validation through AVO modelling
AVO curve at top of water sand
in well H. AVO anomaly is not
seen. (Shear wave log not
recorded)
Poisson’s Ratio
0.5
0.4
0.3 If VP/VS = 2, then = 0
Poisson's Ratio
0.2
If VP/VS = 1.5, then = 0.1 (Gas Case)
0.1
0 If VP/VS = 2, then = 1/3 (Wet Case)
Stack
Gather
Increase of
amplitude with
offset
Training on Seismic data interpretation (142)
Principle of AVO Analysis
With increase
in gas
P-wave saturation, P-
wave velocity
drops
dramatically,
Velocity (m/s)
but S-wave
velocity only
S-wave increases
slightly.
Water Saturation
Training on Seismic data interpretation (143)
Estimating intercept (A) and gradient (B)
Offset
+A
+B
sin2q
Time -B
-A
The Aki-Richards equation predicts a
linear relationship between these
amplitudes and sin2θ.
Regression curves are then
calculated, to give A and B values for
each time sample. A(θ) = Ao + BSin2θ
Training on Seismic data interpretation (145)
Estimating intercept (A) and gradient (B)
Gradient (B)
Base II P Base II
Base I
Base III
Top IV
Intercept (A)
Base IV
Crossplot
Showing “Wet” Trend
Anomalies Vp
Top II Top II P 2
Vs
Training on Seismic data interpretation (147)
Rutherford/Williams Classification
•Rutherford and Williams (1989) derived the following
classification scheme for AVO anomalies, with further
modifications by Ross and Kinman (1995) and Castagna (1997):
AVO attributes
(1) AVO Product : A*B
(2) Scaled Poisson’s Ratio Change : A+B
(3) Shear Reflectivity : A-B
The AVO product shows a positive response at the top and base
of the reservoir For Class-III Anomaly
GC26-B(shift
Near
Far
Daman 3
Daman4
Increase in amplitude
Training on Seismic data interpretation (151)
Near
Far
Daman 3
Daman 4
Mahuva top
Inversion
• Inversion is the process of extracting, from the
seismic data, the underlying geology which gave rise
to that seismic.
• Traditionally, inversion has been applied to post-
stack seismic data, with the aim of extracting
acoustic impedance volumes.
• Recently, inversion has been extended to pre-stack
seismic data, with the aim of extracting both acoustic
and shear impedance volumes. This allows the
calculation of pore fluids.
• Another recent development is to use inversion
results to directly predict lithologic parameters such
as porosity and water saturation.
Training on Seismic data interpretation (154)
Inversion Non-Uniqueness
• All inversion algorithms suffer from “non-
uniqueness”.
• There is more than one possible geological model
consistent with the seismic data.
• The only way to decide\between the possibilities is to
use other information, not present in the seismic
data.
• This other information is usually provided in two
ways:
– The initial guess model
– constraints on how far the final result may deviate from the
initial guess
• The final result always depends on the “other
information” as well as the seismic data
Training on Seismic data interpretation (155)
Acoustic Impedance
Geology Seismic
Surface
Wavelet
Inverse Model
Inverse
Wavelet
Limestone
Sandstone
Basement
What is Pitfall ?
• Any concealed danger or trap for an
unsuspecting person
Webster
• A pitfall is a situation or circumstance in which
you can easily make a mistake if you proceed
with erroneous assumptions or observations or
incorrectly apply problem-solving technique
Training on Seismic data interpretation (167)
• Why Pitfall ?
– Earth crust is not a well-ordered system and
that the seismic events resulting therefrom are
even less well-ordered.
– Some pitfalls are dug by nature or some are of
our own doing
Training on Seismic data interpretation (168)
DHI
Training on Seismic data interpretation (169)
DHI-fizz water
Training on Seismic data interpretation (170)