Bp sesmic interpretation

Seismic data interpretation
by
Dr. HARILAL, DGM (GP)
GEOPIC, ONGC, Dehradun
27th March 2012 Rajahmundry

Seismic data interpretation (2 )
What is seismic data ?

Outline
• Introduction
– Objectives
• Basic concepts : Geophysical and Geologic
• Reflection, Velocity, P- and S-waves
• Polarity, Phase, Resolution, Detectability
• Occurrence of petroleum in subsurface
• depositional environments, rock types, faults and folds
• Project and data management and overview of seismic data loading
• Structural Interpretation
– Calibration, horizon and fault correlation and map generation
• Stratigraphic interpretation
– Seismic Attributes: Physical and structural
– Multi-attribute seismic guided log property mapping
• Direct Hydrocarbon Indictors (DHI) and their validation
• Rock Physics, AVO and Inversion (Impedance and rock properties)
• Time Lapse Reservoir Monitoring (4D)
• Unconventional Reservoirs (Shale –Gas)
• Pitfalls

Interlink in Seismic API in Exploration
Data Processing
and Imaging
Data Visualization, Interpretation and
modeling, Prospects and Locations
Data Acquisition
Drilling
(Hydrocarbon and Data)
(log, VSP, Model, test data)

What is seismic interpretation?
• Interpretation is telling the geologic story contained
in seismic data.
• It is correlating the features we see in seismic data
with elements of geology as we know them.
• Seismic interpretation is a type of inverse problem:
• Given a set of data collected by experimentation and
the physical laws controlling the experiment
and then
 Develop a model of the physical system (subsurface)
being studied

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
– Structural, stratigraphic, combination
• 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, saturation, etc.).

Basic concepts in seismic
interpretation

Generation of seismic signal
Travel time (t) : infer and delineate structure, paleo structure
Amplitude and reflection character: lithology and fluid saturation
Signal + noise
source
Receiver
h
x
(V2t2/4h2 ) – (x2/4h2 ) ) = 1
Traveltime curve is hyperbola
Travel
time

Reflection of seismic waves (normal incidence)
• When a seismic wave
encounters a
boundary between
two different
materials with
different impedances,
some of the energy of
the wave will be
reflected off the
boundary, while
some of it will be
transmitted through
the boundary.
Snell’s Law
Sin θ1/V1 = Sinθ2 /V2= p
θ1
θ2
Angle of refraction
Angle of
incidence
Normal
incidence
V1 ρ1
V2
Boundary
ρ2
1122
1122
0
VV
VV
R






Reflected
P-wave = RP
Reflected
SV-wave
Transmitted
P-wave
Incident
P-wave
Transmitted
SV-wave
VP1 , VS1 , r1
VP2 , VS2 , r2
i
r
r
t
t
If q > 0°, an incident P-wave will produce both P and SV
reflected and transmitted waves. This is called mode conversion.
Reflection of seismic waves (non-normal incidence)

At non-normal incidence, the reflection coefficient depends on other
parameters also, and is estimated by the Zoeppritz equations:
1 1 2 2
1 1 2 2
2
2 2 1 2 2 11
21 1 22 2
1 1 1 2 1 1
1 2 22 2
1 1 2 2
1 1 1 1 1
sin cos sin cos
cos sin cos sin
sin2 cos2 sin2 cos2
cos2 sin2 cos2 sin2
P
S S P S PP
P S S P S
S S SP
P P P
R
R V V V VV
T V V V V
T V VV
V V V
   
   
 
   
 

   
 
  
     
   
    
   
        
 
1
1
1
1
1
sin
cos
sin2
cos2





 
 
 
 
 
  

Zoeppritz derived the amplitudes of the reflected and transmitted
waves using the conservation of stress and displacement across the
layer boundary, which gives four equations with four unknowns.
Inverting the matrix form of the Zoeppritz equations gives us the exact
amplitudes as a function of angle:
Reflection of seismic waves (non-normal incidence)

Direction of energy propagation and particle motion for
compressional (P) waves, horizontal shear (SH) waves, and
vertical shear (SV) waves.
Multi component (P-, Sv-, Sh-wave)

Data phase and polarity critically determine seismic character.
Character is more important than amplitude in directly identifying
hydrocarbons with seismic data.
Once data phase and polarity are determined, hydrocarbon character
can be predicted, and this is of major importance in analysing
prospectivity in younger sediments.
Character is also key in making an effective well tie and thus
correctly identifying seismic horizons.
Polarity and Phase

Polarity Convention
The SEG standard for causal seismic data specifies that the
onset of a compression from an explosive source is represented
by a negative number that is, by a downward deflection when
displayed graphically
Minimum phase
Zero phase
Low Impedance
High Impedance
High Impedance

Top and bottom of a gas reservoir
(low impedance zone) in (a)
American polarity and (b) European
polarity
Polarity Convention
American polarity is
described as: An increase in
impedance yields positive
amplitude normally displayed in
blue. A decrease in impedance
yields negative amplitude
normally displayed in red.
European (or Australian)
polarity is described as the
reverse, namely: An increase in
impedance yields negative
amplitude normally displayed in
red. A decrease in impedance
yields positive amplitude
normally displayed in blue.
American
European
Impedance

In an extremely simplified way seismic pulses displayed on seismic
sections can be grouped into two main types, minimum phase and
zero phase.
Phase

Phase
Zero-phase pulses consist of a central peak and two side lobes of
opposite sign and lesser amplitude. Here the boundary is
located at the central peak and not at the wavelet onset as is the
case for minimum-phase pulses.
Given the same amplitude spectrum, a zero-phase signal is
always shorter and always has greater amplitude than the
equivalent minimum-phase signal; it therefore has a greater
signal/noise ratio.
A minimum-phase pulse has its energy concentrated at its front,
and is thought to be typical of many seismic signals. The pulse is
said to be "front loaded," with its onset at the acoustic-
impedance boundary.

Subsurface features which can generate sufficiently
high reflections to be useful for interpretive assessment
of phase and polarity
Assessment of phase and polarity

Resolution
The ability to separate two features that are close
together
The resolving power of seismic data is always measured
in terms of seismic wavelength (λ=V/F)
Limit of (resolution) separability = λ/4
The predominant frequency decreases with depth because the
higher frequencies in the seismic signal are more quickly
attenuated. Wavelength increases with depth.
Resolution decreases with depth
For thinner intervals amplitude is progressively attenuated until
Limit of visibility=λ/25 is reached when reflection signal
becomes obscured by the background noise

Resolution: has both horizontal and vertical aspects
Resolution
Attenuation of noise
in data processing
Broad bandwidth by
maximum effort in
data acquisition
Horizontal
minimum size
Fresnel zone
sampling
Seismic migration
Vertical minimum
Thickness
Deconvolution
Wavelet
Can be improved through both acquisition and processing
Migration improves
horizontal resolution
Deconvolution
enhances vertical
resolution

Limit of Separability
Age of rocks Very
young
young medium old Very old
Depth Very
shallow
shallow medium deep Very deep
Velocity
(m/s)
1600 2000 3500 5000 6000
Predominant
frequency
70 50 35 25 20
Wavelength 23 40 100 200 300
Separability 6 10 25 50 75

Limit of visibility
• Factors affecting the
visibility
– Impedance contrast of
the geologic layer of
interest relative to the
embedding material
– Random and
systematic noise in the
data
– Phase of the data or
shape of seismic
wavelet
– It may be less than 1
m to more than 40 m
Limit of visibility
S/N Example Limit
Poor Water sand poor
data
λ/8
Moderate Water or oil sand
fairly good data
λ/12
High Gas sand good
data
λ/20
Outstanding Gas sand
excellent data
λ/30

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
f
tv
rf
2

where
rf = radius of the Fresnel zone.
v = average velocity.
t = two-way time in seconds.
f = dominant frequency in hertz.
Post – migration Fresnel zone size is given by /4 = V/4Fλ
On migrated section vertical and horizontal resolutions are of
the same order

Horizontal resolution
Pre migration Fresnel zone= 1130 m
Post migration Fresnel zone= 40 m
T = 2 sec
F= 25 Hz
V= 4000 m/s
Effects on Fresnel
zone size and shape
on 2-D (ellipse) and
3-D migration
(circle)

Velocity tends to increase with
depth due to compaction and
diagenetic effects.
Frequency decreases due to
attenuation of the seismic wave.
Higher frequencies are attenuated
more than lower frequencies for a
particular path length.
With increasing depth vertical
and lateral resolution decreases,
and interference effects become
more pronounced as the pulse
length increases ( due to lower
frequency).
Relationship among wavelength, frequency and velocity

Occurrence of Petroleum in Subsurface
How to locate
subsurface trap?
How to locate reservoir
rock?
How to find H/C in the
trap?
How to monitor
changes due to
production?
Gas-Oil
Gas-
Water
Oil-
Water
Original
fluid
contact

How to map structural closure, fault closure, pinchout and salt-related structures?

How to identify different types of traps (Structural and/or stratigraphic )

How to develop unconventional shale-gas reservoirs which are fine-grained
organic-rich shales of tight-porosity and ultralow-permeability and acting as
the source, seal, and the reservoir rock
Unconventional
shale-gas in the
syncline

Sedimentary rock types
• Sandstone
– Low porosity sands: very high value of velocity and density
– High porosity sand: low values of density
– High porosity or unconsolidated: low values of density
– Interfaces between low-porosity sand and shales : good reflector.
– Interfaces between high-porosity sand and shales: poor reflector
• Shale
– Important characteristic is readiness to compact after deposition
– Draping of shale sediments over sandstone and carbonate bodies
– Depth is most important factor affecting the seismic properties
• Carbonates
– At a given porosity carbonate rocks have higher density than sands and shales
– At zero porosity, carbonates have higher velocity than sands and shale
– Good reflection from sand and tight carbonates and shales and carbonates
– Secondary porosity observed in the carbonates develop complications
• Fractures lower velocity without developing significant porosity or change in
density
• Vugs can produce depression in velocity in the bulk material, but it may not appear
on sonic logs

Sedimentary Depositional Environments
Alluvial
Beach
Deltaic
Reef
Deep water
marine
Glacial
Lacustrine
Aeolian

Classification of depositional environments, based on the relative
contributions of nonmarine and marine processes.
(After : O. Catuneanu)

Continental
Glacial: not good reservoir rocks(fine material silt and
clays)
 Alluvial: generally not a reservoir rock
Aeolian: Complex heterogeneous reservoir
Fluvial:
Braded system: Deposited from a river system,
Potentially good reservoir rock. Stratigraphic traps are
not common
Meandering System: Potentially good reservoir
rock. Stratigraphic traps are common
Lacustrine:
Good source rocks:

Transitional
Deltaic : Good reservoir characteristics. Tectonic and
Stratigraphic traps are abundant.
 Delta Plain: sub aerial lowland part
Upper delta Plain: above the tidal or marine influence, migratory
distributary channel, inter-distributary
Lower delta plain: Low tide mark to limit of tidal influence, Bay
fill and abandoned distributary-fill
Delta front: Subaqueous, High energy sub environment
Sheet sands
Distributary mouth bar
longshore
Prodelta: below the effective depth of wave erosion
River dominated delta
Wave dominated delta
Tide dominated delta whare tidal range is high reversing
flow in distributaries

Marine
Shallow water marine
Tidal sand ridge: elongated sand bodies formed by tidal
currents
Linear submarine sand bars
Carbonate environment
Good reservoir characteristics with source and seal
rock potential
Deep water marine
Deposited in large body of water below the actions of
waves resulting from sediment gravity flow mechanism.
Turbidites are sediments which are transported and
deposited by density flow in deep ocean.

Alluvial fans

Braided channel
Meander channel
A braided river consists of a network of
small channels separated by small and
often temporary islands called braid bars.
Braided streams occur in rivers with high
slope and/or large sediment load. Braided
channels are also typical of environments
that dramatically decrease channel depth,
and consequently channel velocity, such as
river deltas, alluvial fans and peneplains.

River dominated delta
Mississippi river delta
Nile River delta, as seen from
Earth orbit. The Nile is an
example of a wave-dominated
delta that has the classic Greek
delta (Δ) shape after which River
deltas were named
If waves, tides and longitudinal
currents are week, and if the
volume of sediments carried by
river is high, rapid seaward
progradation takes place, and a
variety of river dominated
depositional environments develop.
The main sediment is deposited in
distributary mouth bar.
Nile river delta

Tide dominated Ganges delta
Tide-dominated deltas: where
river mouths hit the sea in areas
affected by large tidal ranges, the
delta shape can be extensively
reshaped by the twice a day flood
and ebb tidal currents moving in
and out of the river mouth. This
usually happens in bays and
estuaries where the river mouth is
protected from much wave
activity. The relentless in and out
currents of tides can sculpt the
sediment into elongate tidal bars.

Normal Faults caused by
extension or tension
Reverse Faults caused by
compression
Thrust Faults caused by
compression
Strike Slip Faults caused by
lateral shear, lateral shift
Tectonic elements (faults)

Listric (growth) faults, are generally parallel to ancient shoreline and concave seaward
in plan. Significant source of sediment is landward.
• Thickening of correlatable intervals on the downthrown side
•Often, rollover of downthrown reflections into the fault, and sometimes a rising of the
upthrown reflections towards fault
•Concavity towards the basin in the section, with angle approaching 60 degrees at the
top and attenuating to bedding-plane dip at the bottom
•Antithetic (or adjustment) faults , generally develop on the crestal part of rollover

• Introduction
– Objectives
• Basic concepts
• Project and data management and overview of
seismic data loading
– Calibration, horizon and fault correlation and map generation
• Pitfalls

Interpretation functionalities
(through workstation)
• Project and data management
– (Project and survey)
• Data conditioning (interpretive processing)
– Scaling, filtering, wavelet processing etc.
• Integration of various types of data (seismic, well etc.)
• Calibration (well to seismic tie)
– Synthetic seismogram generation and correlation with seismic
• Visualization
• Horizon and fault Correlation
• Generation of time maps
• Depth conversion and generation of depth maps
• Special studies
– Seismic Attributes
– Direct hydrocarbon indicators (DHI) and AVO
– Seismic Inversion
– Time lapse reservoir monitoring
• Integration of structure, attributes, impedance, geologic model etc.
• Prospect generation and location identification
• Reports/proposals

IEEE floating point : 1.18E-38 to 3.37E+38, and zero
IBM floating point: 5.4E-79 to 7.2E+75, and zero
8 bit Integers from -128 to +127
16 bit Integers from -32768 to +32767
32 bit Integers from -2147483648 to +2147483647, or
in scientific notation, approximately +2.147E+10
Overview of data loading
Storage concepts
Storage in lower bit (8 or 16) during interpretation to maximize
response time and minimize the storage space
Why scaling is required?

Overview of data loading
Storage in lower bit (8 or 16) during interpretation to maximize
response time and minimize the storage space
Statistical distribution of amplitude in a data volume
-128 1280
Numbers of samples in data volume
Amplitude
Most interpretive
time
May be bright spot
(DHI)
May be bright spot
(DHI)

• When data are stored in lower bits (8 bit or 16 bit) clipping should be avoided.
• Clipped data are unfit for any attribute or stratigraphic analysis
Overview of data loading (Scaling and clipping)
Clipped Normal (properly scaled)

• Introduction
– Objectives
• Basic concepts
– Calibration, horizon and fault correlation and map
generation
• Pitfalls

Data volume
Vertical section Horizontal section Surface tracking
Arbitrary
line line
Crossline Time slice Depth slice Horizon
slice
Horizon
attribute
Fault slice
Composite Display integrating 2 or more
Elements of display
Display products from a 3-D seismic data

A seismic section showing colour convention and other display elements
Elements of display
Positive amplitude Blue
Negative amplitude Red
Doublegradationalscheme

3-D Cube of
seismic data
Elements of display
Time slice
In-linecross-line

Successive time
slices depict the
anticlinal closures
Time 2430
Time 2390
Elements of display
Vertical section
and time slices

Composite Display

WATER
BOTTOM Canyon cut
3800
3300
2600
2100
1800
1600 m/s
1500
Velocity data Display

Horizon Correlation
– Identification of log markers on seismic section
– Identification of sequence boundaries
• Reflection configuration and termination pattern, Onlap,
Top lap and Down Lap
– Tracking
• Auto and manual mode
– Auto dip and correlation
• 3D voxel tracking
• Voxel detection and horizon converter
Fault correlation
– Picking on vertical and horizontal sections
– Visualizing DAT, coherency, Event similarity and
instantaneous attribute volumes
– Analyzing Dip Azimuth and other attributes
extracted along horizon
– Fault plane correlation
– Computing Heave and throw of faults

Synthetic to seismic matching before starting the interpretation. Understand
geologic elements (model) and their geophysical responses (Trace).
DT RHOB IMP GR LLD RC
Black
traces
Seismic
Red
traces
synthetic
Well to seismic tie (Synthetic Seismograms)

Overlay of synthetic and logs on
seismic sections. Here synthetic is
shown by yellow trace
Identification of log markers on seismic section through synthetic
Well to seismic tie

min
max
+/- zero crossing
-/+ zero crossing
Minimum. Autopicking along the minimum amplitude values of the event.
Zero. Autopicking along the zero crossing of the event.
Maximum. Autopicking along the maximum amplitude values of the event.
Zero (+/-). Autopicking zero amplitudes where the signal crosses from positive
to negative.
Zero (-/+). Autopicking zero amplitudes where the signal crosses from
negative to positive.
Illustration of Points Specified by Onset Type
Auto Picking

For stratigraphic and reservoir
interpretation, the correct point to
pick on the waveform is the crest.
Picking on the waveform
flank may be good for
structure map but amplitude
extracted from that point is
inadequate and may be
meaningless for stratigraphic
interpretation.
Principles of seismic
resolution tell us that all the
reflecting interfaces within an
interval equal to one-quarter
wavelength (λ/4) will
contribute to each reflection.
Thus we should identify the
relevant reflection and pick on
the crest of it. The crest of the
waveform is the correct point
to pick for zero-phase data.

Seismic reflection patterns
Seismic reflections follow impedance contrasts that coincide with
stratal surfaces, which are surfaces where depositional processes occur
at a fixed moment in geologic time (Chronostratigraphic surfaces)
The study of rock relationships within a chronostratigraphic
framework of genetically related strata bounded by surfaces of
erosion or nondeposition or their correlative conformities.
Sequence stratigraphy

Sequence boundaries may be identified from bounding
configuration and can be taken for horizon correlation
At Base: Base onlap, downlap or parallelism
At Top : Toplap, erosional truncation or parallelism to boundary
Configuration at sequence boundaries

Horizon correlation
Mounded feature
onlapdownlap

Coherency time slice at
1600 ms showing strike
slip and normal faults.
Signatures on vertical
sections are shown on
sections AA’ and BB’
Strike slip fault
(dextral)
Normal fault
B
B’
A A’
Fault correlation

Seismic section showing Normal faults
Normal faults
A A’
Fault correlation

Flexure
Time-slice at 1300 ms from dip-azimuth
volume showing faults and flexures
Fault correlation

Curvature attribute showing faults and flexure
Maximum NegativeMaximum Positive
Fault correlation

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

Map at correlated horizons
Map generation

• Introduction
– Objectives
• Basic concepts
– Calibration-synthetic seismogram, horizon and fault correlation and map generation
• Pitfalls

 Stratigraphic interpretation is based
on the identification and mapping of
changes in reflection character and
correlating it with analogues.
 A shape or pattern which is unrelated
to structure may be due to depositional,
erosional, lithologic or other features of
interest
 Horizontal sections and horizon slices
can provide a bird’s-eye-view of
ancient stratigraphy
 Seismic attributes improve the
mapping of stratigraphic features
which may be “hidden” on normal
sections
(In this lecture Stratigraphic interpretation is
covered in attribute section)
Stratigraphic interpretation
Horizon slice showing a
stratigraphic feature (Channel). Let
us see how this type of
stratigraphic features are mapped.

Links between Seismic attributes and
rock/reservoir properties
• Seismic responses are based on
Impedance contrast
• Depends upon rock properties
– Lithology, porosity, fluid saturation
– Many others
GR IMP

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
•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.

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)

Seismic attributes
Attribute Information
Time-derived Structural information
Amplitude-derived Stratigraphic and
reservoir
Frequency-derived Stratigraphic and
reservoir
Attenuation Permeability (future)

(After A. R. Brown)

Seismic attribute (extraction)
• 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 (strat
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

Seismic attributes (instantaneous)
Attribute name Mathematical definition
Reflection strength A(t) =√ { f2 (t) + h2 (t) }
Instantaneous phase θ (t) = tan-1 { h(t) / f(t) }
Instantaneous frequency w (t) = d θ / dt
Quadrature Trace h(t) is Hilbert transform of f(t),
a 90 degree phase shift of f(t)
The basis for all the complex trace attributes is the idea of a complex trace.
Complex-trace analysis separates amplitude and phase information.
Taner et al. (1979) show that a seismic trace f(t) can be considered as the real
part of a complex trace, F (t) = f ( t) + ih(t). The h(t), quadrature or
imaginary component is determined from f(t) by Hilbert Transform

Seismic attributes (instantaneous)
Name Description Applicability
Instantaneous
real amplitude
Default expression of
seismic trace data
Structural and stratigraphic, Isolate high
and low amplitude areas such as bright
and dim spots
Instantaneous
amplitude
Amplitude
envelope., Reflection
strength
Identifying bright / dim / flat spot s.
Often used to determine lateral
fluid ,lithologic, and stratigraphic
variations in reservoir
Instantaneous
quadrature
amplitude
Time domain vibration
amplitude with phase
delay of 90 degree from
instantaneous real
amplitude
The phase-delay feature is useful in
quality control of instantaneous phase,
identify some AVO anomalies from thin-
layers since they may only be observable
at specific phase
Instantaneous
phase
Angle whose tangent is
h(t) / f(t)
Enhances weak intra-reservoir events.
One of the hydrocarbon indictors
Instantaneous
frequency
Derivative of
Instantaneous phase
To estimate seismic attenuation. Drop of
high frequencies due to oil and gas
reservoir. Measure cyclicity of geologic
intervals (After Chen et al., TLE May-1997)

Seismic attributes (Apparent Polarity)
recorded
seismic data
Reflection
strength
Apparent
polarity
Apparent polarity
displays have a
blocky
appearance.
In the recorded seismic
data, the gas-oil contact
(flat horizon)
corresponds to the center
of a zone with large
amplitudes and positive
polarity. In the apparent
polarity display, this
interface is even more
prominent since the effect
of apparent polarity is to
highlight both the
dominant amplitude and
the polarity of the event.

Sweetness Sr(t) is defined as response amplitude ar(t) divided by
the square root of response frequency fr(t),
Seismic attributes (Sweetness)
Sweetness is an empirical seismic attribute designed to identify
“sweet spots,” places that are oil and gas prone. In young clastic
sediments, sweet spots are often characterized seismically by high
amplitudes and low frequencies. Sometimes reflection strength and
instantaneous frequency are used instead of response amplitude and
response frequency.
Sweetness closely resembles reflection strength. Sweetness
anomalies of most interest are those that are relatively stronger than
their corresponding reflection strength anomalies.

Thin Bed Indicator is a hybrid complex trace attribute invented by
T. Taner. It is defined as the absolute value of the instantaneous
frequency minus weighted average instantaneous frequency.
The weighted average instantaneous frequency represents the
background trend in the instantaneous frequency. Its resolution is
set by a window length. Use shorter window lengths to subtract
more detail from the instantaneous frequency, and use longer
window lengths to subtract less detail.
The spikes in instantaneous frequency are likely caused by
interference in reflected wavelets.
These can be used to indicate thin beds, pinchouts, and bedding
patterns in general.
Removing the background trend in frequency makes these spikes
stand out better.
Seismic attributes (Thin bed indicator)

Real amplitude, instantaneous amplitude (reflection strength),
phase and frequency sections
Instantaneous amplitude
Instantaneous Phase Instantaneous Frequency

Example showing 3-D visualization of derived frequency volume
depicting sand dispersal trends
Sand
dispersal
trends

Time derived Horizon attributes
• Residual
– Arithmetic difference between high precision automatically-
tracked time map and its spatially-smoothed equivalent
• Subtle faults and data collection irregularities
• Dip (dip magnitude)
– Ms/trace
– Provides structural details (faults)
– Quality control the performance of horizon tracking process
• Azimuth (direction of dip)
– Subtle faults
• Dip-Azimuth
– Combines dip and azimuth attributes onto the one display
• Subtle faults
• Curvature
– Derivative of dip and azimuth
• Subtle faults

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.
Time derived Horizon attributes (Dip)
A dip map shows the magnitude of the time gradient. It is constructed
by comparing each sample of the horizon with two adjacent samples in
orthogonal directions. A plane is fit through the three points. The plane
has a magnitude of dip in milliseconds per unit distance x 1000.
Individual dip values are not particularly meaningful. but the
relative difference between various dip values is significant which
can be effectively displayed by a simple dichromatic ramp color
map (black and white)

Where :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.
An azimuth value of 0o is
aligned with true north as
shown in the figure.
The 0o axis always coincides with North,
90o coincides with East, 180o with South,
and 270o with West.
Time derived Horizon attributes (Azimuth)
The azimuth map is closely associated with the dip map. It
shows the direction of maximum dip. Like the dip map, it is
constructed by comparing each sample of the horizon with two
adjacent samples in orthogonal directions. A plane is fit through
the three points. The plane has a direction of dip in degrees.

It provides a technique for maximizing the effectiveness of a
single structural attribute map. Its usefulness is based on the fact
that, when mapped on standard azimuth maps and standard dip
maps, faults may or may not appear, depending on the relationship
of their dip and azimuth to the dip and azimuth of the horizons.
Azimuth maps exhibit the following tendencies:
•Faults show up best where dip direction of the fault is
opposite the dip direction of the beds.
•Faults show up poorly where dip direction of fault is similar
to dip direction of the horizons.
Dip maps exhibit the following tendencies:
•Faults show up best where dip angle of fault plane is opposite
direction of beds.
•Faults show up poorly where dip angle is close to that of
horizon.
Time derived Horizon attributes (Dip Azimuth)

Time derived Horizon attributes (Curvature)
Curvature is defined as the reciprocal of the radius of a circle that is
tangent to the given curve at a point.
Thus curvature will be large for a curve that is “bent more” and will
be zero for a straight line.
Mathematically, curvature may simply be defined as a second-order
derivative of the curve.
If the radii of the circles at the point of contact on the curve are
replaced by normal vectors, it is possible to assign a sign to
curvature for different shapes
Diverging vectors on the curve are associated with anticlines,
converging vectors with synclines, parallel vectors with planar
surfaces (which have zero curvature).

Time derived Horizon attributes (Curvature attributes)
Type Definition Interpretation
Mean Avg of min and max High and low (throw of fault)
Gaussian Product of Min. max Distortion of surface
Dip Along the dip direction Rate of change of dip in maximum
dip direction
Strike Along the strike direction patterns connecting high with lows
on both sides of the main fault
Most-
Negative
greatest negative value Synclinal and bowl features.
Positive values of the most
negative indicate a dome feature
Most-
Positive
greatest positive value anticlinal and domal features.
Negative values of the most
positive indicate a bowl feature

3D display showing topography and fault systems along Panna-4 reflector
Faults
Median high
Highs within
Central Low
High

Dip –Azimuth horizon slice at UBT showing variation in magnitude and azimuth
of dip of reflector. Faults, flexures, channels and other discontinuities are
highlighted
Red ➛ yellow indicates dip to the northeast.
Red ➛ blue indicates dip to the northwest.
Green ➛ yellow indicates dip to the southeast.
Green ➛ blue indicates dip to the southwest.
structure (faults and flexures) and
stratigraphy (channels and fans)
and to the recognition of
geophysical artifacts (statics,
acquisition effects, and noise).

Time derived Window attributes
• Coherence
– Volume of discontinuity (faults and other
boundaries) free from interpretive bias
• Subtle faults
• Stratigraphic features (channel edges)
• Lithologic boundaries
– Coherence horizon slice extracted along
smoothed picked horizon
• Enhanced visibility for Channel and faults

Horizon and windowed amplitude
Horizon Amplitude
• Reservoir study
• No Contamination
• Highly dependent on
tracking accuracy
• Horizon identification
Critical
• Horizon name only
required
Windowed Amplitude
• Reconnaissance
• Subject to contamination by
other amplitudes
• More immune to tracking
errors
• Horizon interpretation less
critical
• Record required of window
definition and type of
amplitude extracted

Coherency
slice
at 1500 ms

• Discontinuity volumes
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.
Carbonate buildup
coherency
Carbonate build up
boundary
Top of Carbonate
build up

Spectral Decomposition
Spectral decomposition decomposes the seismic data with normal
frequency bandwidth into a set of sections having discrete or very
narrow bandwidth. It is based on the concept that reflections from
a thin bed have a characteristic expression in a particular
frequency domain that is indicative of temporal bed thickness.
Very effective in mapping stratigraphic features (Channels)

Spectral Decomposition at 28 Hz frequencyCoherency horizon slice

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

Seismic facies
map showing
differences in
lithofacies.
The red curve shows
cumulative
difference between
adjacent synthetic
classes and green
curve shows
correlation between
selected trace and
synthetic traces.

Seismic facies
map with
Neural
Network
technique of
a channel
showing
facies within
and on the
banks of
channel

Horizon slice
Flattened section
• Channel
mapping
through
horizon slice

Amplitude horizon slice of Pay-2 reflector extracted from PSTM reflection strength
volume, with overlay of time structure. Anomaly is grouped into sectors I and II with
gap (green polygon) at well-A.

Opacity cube showing massive buildup and slope facies (off-buildups). Here low amplitude
has been made opaque and surfaces are displayed in their color. (b) Coherency horizon slice
of buildup top showing buildup boundary and other discontinuities. (c) 3D map
simultaneously showing structural configuration and amplitude variation (in color) of
buildup top reflector. Illumination by directional light has highlighted the lowamplitude
(yellow) massive buildup in the center, and highamplitude (red) surrounding talus (off-
buildup) features.
Horizon slice

Geobody (Channel-fills) mapped through voxel
based 3-D visualization and detection method.
Panna-2
Panna Formation

Multi-attribute Seismic-guided log-property mapping
• Seismic attributes can be used for mapping of log properties
by finding the interrelationships between logs and seismic
data at tie points
– Multi-attribute linear regression
– Artificial neural network techniques
• Training the seismic to predict the reservoir parameter of
interest at the tie (well) locations
• Applying results of training to the seismic volume
• Any log property volume , original: Velocity, Density etc., or
derived: porosity, Saturation, etc. can be estimated

Original log Modelled log Analysis window
Correlation Coefficient: 94% , Average Error: 24.16
Time(ms)
Prediction of DT log through Probabilistic Neural Network (PNN)
DT (µs/m) DT (µs/m) DT (µs/m)
DT (µs/m)
1800
1850
1900
1950
1950
Multi-attribute seismic guided mapping

DT (µs/m)
DT(µs/m)
NPHI (fraction)
NPHI(fraction)
CR: 94% CR: 86%
Crossplot between actual and
predicted Sonic (DT ) logs.
Correlation Coefficient: 94% ,
Average Error: 24.16
Crossplot between actual and
predicted Neutron-Porosity
(NPHI) logs.
Correlation Coefficient: 88%,
Average Error: 0.036

CDT (µs/m)
Prediction of DT (Sonic ) log
Multi-attribute Seismic-guided log-property mapping
Low DT (high
velocity) areas
indicate good
sand zones

• Introduction
– Objectives
• Basic concepts
– Calibration-synthetic seismogram, horizon and fault correlation and map
generation
• Pitfalls

Direct hydrocarbon indicator (DHI)
A measurement that suggest the presence of
a hydrocarbon accumulation

Hydrocarbon indicators. All indicators can have causes other than
hydrocarbons; a case for hydrocarbon accumulation is stronger where
several indicators agree
DHI

DHI• Bright spot
– Water sand has lower
impedance than
embedding medium
and impedance of gas
sand is further
reduced.
– Top and base
reflections show
natural pairing
– If sand is thick
enough, flat spot or
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)
Peak on
synthetic
seismogram

A flat spot associated with a bright spot offshore Indonesia. The
flat spot is fairly flat, but actually dipping slightly to the left - that
is, against structure.
DHI

DHI
• Dim spot
– Water sand has
higher impedance
than embedding
medium and water is
replaced by gas
impedance is reduced
– Contrast is reduced at
upper and lower
boundaries and
reservoir is seen as
dim spot
– Flat spot can be
expected at the point
where the dimming
occurs.
– More common in
deeper sandstone
reservoirs where shale
impedance is lower
than sandstone
impedance

DHI
• Polarity reversal
– Water sand is of
higher impedance than
enclosing medium and
gas sand impedance is
lower than enclosing
medium
– The polarity of
reflections from water-
sand-shale interface
and gas-sand-shale
interface have opposite
sign and thus polarity
reversal
– Flat spot from GWC
may show bright
amplitude .
– More common in
medium depth range
sandstone reservoirs

DHI
Low impedance
reservoir embedded
in high impedance
High impedance
reservoir embedded
in medium
impedance
High impedance
reservoir
embedded in Low
impedance
Hydrocarbon fluid
effects in seismic data

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

RMS-Amplitude anomaly maps showing envisaged pay sands
Amplitude slice at 24 ms w.r.t.
Daman4 reflector
Pay-1
Pay-3
Bright Spot

Flat Spot

Occurrence of DHI
Depthofmaximumburial
• Nature of DHI (bright, Dim
or Phase change) depends on
relative impedance of
hydrocarbon sand, water
sand and shales
• Impedance increases with
depth and age of rock
• Visibility of indicators
decreases with depth
• There is a cut-off below
which no hydrocarbon
observations will be
possible
• DHI may occur at great
depth for younger rock
• Older rocks may show
indicators at shallower
depth.

Gas sand, water sand and
shale acoustic impedances
all increase with depth and
age at different rates. The
crossover points define
whether the hydrocarbon
indication is bright spot,
phase change or dim spot.
Overpressure moves the
shale line to lower
impedance, thus increasing
the effect of dim spot and
decreasing the effect of
bright spot
Occurrence of DHI

Validation of DHI
• AVO
– In many practical cases gas sands show an increase
of amplitude with offset
• AVO attribute Cross-plot
• LMR inversion
• Elastic inversion
– Many difficulties of theoretical and practical nature
• Data is pre-stack hence lower S/N
• Multidimensional, many mode of display
Amplitude anomalies observed in a seismic event can be due to
lithology, over pressures and the presence of hydrocarbons

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

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

• Introduction
– Objectives
• Basic concepts
• Rock Physics, AVO and Inversion (Impedance and
rock properties)
• Pitfalls

Rock physics
•Bridges seismic data and reservoir properties and parameters
•Seismic data are commonly analysed for
 lithology
 porosity
 pore fluids
Saturation
• Technologies based on rock physics are:
4-D seismic reservoir monitoring
Seismic lithology discrimination
Direct hydrocarbon detection with “bright-spot”
Angle dependent reflectivity (AVO) analyses

Rock physics
Factors influencing seismic properties of sedimentary rocks
(with increasing importance from top to bottom)
Rock properties Fluid properties Environment
Compaction Viscosity Frequency
Consolidation history Density Stress history
Age Wettability Depositional environment
Cementation Fluid composition Temperature
Texture Phase Reservoir process
Bulk density Fluid type Production history
Clay content Gas-oil,, gas-water ratio Layer geometry
Anisotropy Saturation Net reservoir pressure
Fractures ,Porosity,
Lithology Pore shape

VP VS
ρ
SW = 50%
VP VS
ρ
SW = 100%
Input Logs Output Logs
(1) Fluid Substitution
The basic use of the Biot-Gassmann equations is to “substitute” or replace
the fluids in a set of target layers with another set of fluids.
In this case, VP, VS, and ρ must all be known for the input logs, along with
the fluid content (SW). Generally all three logs are changed within the
target zone.
Using the Biot-Gassmann Equations

VP ρ
SW = 50%
VP VS
ρ
Input Logs Output Logs
(2) Calculating Vs
The second use of the Biot-Gassmann equations is to calculate a VS
curve, which has not been measured in the well. Either KDRY is assumed
known or the mudrock equation is assumed to hold for wet sands.
In this case, VP and ρ must both be known, along with the fluid content
(SW). The VP and ρ logs are unchanged, and a new VS log is created.
SW = 50%
Using the Biot-Gassmann Equations

AVO
Amplitude variation with angle/offset (AVA/AVO)
The variation in amplitude of a seismic reflection with
angle of incidence or source-geophone distance. It depends
on changes in velocity, density, and Poisson’s
ratio. It is often used as a hydrocarbon gas indicator
because gas generally decreases Poisson’s ratio and
often increases amplitude with incident angle/offset.
The input data for this process consists of the following
elements:
• One or more well logs.
• A pre-stack seismic volume. This volume has usually
been processed to the final CDP gather stage.

In seismic gather reflection coefficient at an incidence angle θ
is given by: R(θ ) = R0 + BSin2θ
Where: Ro: RC at zero offset and R (θ) : RC at angle θ
B is a gradient term which produces the AVO effect. It is
dependent on changes in density, ρ, P-wave velocity, Vp, and S-
wave velocity, Vs.
Principle of AVO Analysis
Reflector
S3 S2 S1 S0R0 R1 R2 R3

 222
sintanCsinBA)(R 
where:
p
P
2
P
S
S
S
2
P
S
p
P
p
P
V
V
2
1
C
V
V
2
V
V
V
V
4
V
V
2
1
B
V
V
2
1
A




























A more intuitive, but totally equivalent, form was derived by Wiggins.
He separated the equation into three reflection terms, each weaker
than the previous term:
Wiggins’ Version of the Aki-Richards Equation

Principle of AVO Analysis
AVO response depends on properties of P-wave velocity (VP), S-wave
velocity (VS), and density (ρ) in a porous reservoir rock. This involves
the matrix material, the porosity, and the fluids filling the pores.
R(θ ) = A + BSin2θ +Ctan2θ *Sin2θ
Poisson’s Ratio
K = the bulk modulus, or the reciprocal
of compressibility, µ = The shear
modulus, modulus of rigidity

-0.2
-0.1
0
0.1
0.2
0.3
0.4
0.5
0 1 2 3 4 5 6 7 8 9 10
Poisson'sRatio
Vp/Vs
Vp/Vs vs Poisson's Ratio
Gas Case Wet Case
If VP/VS = 2, then  = 0
If VP/VS = 1.5, then  = 0.1 (Gas Case)
If VP/VS = 2, then  = 1/3 (Wet Case)
If VP/VS = , then  = 0.5 (VS = 0)

Offset
+A
-A
+B
- B
sin2q
Time
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.
Two basic attribute volumes are
produced, Intercept (A) and gradient (B)
The pick amplitudes are extracted at
all times, two of which are shown.
A(θ) = Ao + BSin2θ
Estimating intercept (A) and gradient (B)

AVO Workshop 135
Intercept: A
Gradient: B
AVO Attributes are
used to analyze large
volumes of seismic
data, looking for
hydrocarbon
anomalies.

•Rutherford and Williams (1989) derived the following
classification scheme for AVO anomalies, with further
modifications by Ross and Kinman (1995) and Castagna (1997):
• Class 1: High impedance gas sand with decreasing AVO
• Class 2: Near-zero impedance contrast
• Class 2p: Same as 2, with polarity change
• Class 3: Low impedance gas sand with increasing AVO
• Class 4: Low impedance sand with decreasing AVO
Rutherford/Williams Classification

The generic AVO curves at the top of the gas sand
Class-II
Class-IV
Class-IV gas sand
amplitude decreases with
increasing offset
θ

Anomaly Class 1: High Acoustic Impedance Sandstone
Sandstone Class 1 has relatively high impedance than its seal layer,
which usually is shale. Interface between shale and this sandstone
will result relative high positive coefficient reflection (R0).
Usually this type sandstone is found in offshore area. This
sandstone is a mature sandstone which have been moderately to
highly compacted.
Stack section (Dimpout on productive zone) CDP Gather

Anomaly Class 2: Near-Zero Acoustic Impedance
Sandstone
Class 2 Sandstone has almost the same Acoustic Impedance as
the sealing rock.
This sandstone is a compacted and moderately consolidated
sandstone.
Gradient of sandstone class 2 usually has big magnitude, but
generally it’s smaller than the magnitude of sandstone class 1
Reflectivity of sandstone class 2 on small offset is zero.
This is often blurred due to the presence of noise on our seismic data
The reflectivity suddenly emerge on bigger offset, that is when the
reflection amplitude is located on a higher level than the noise.

Anomaly Class 2: Near-Zero Acoustic Impedance
Sandstone
Polarity change happen if the RC0 is positive
but usually it’s undetected, because it happen on the near offset
where the signal level is under the noise
Sandstone class 2 might and might not be related to amplitude
anomaly on stack data.
If the angular range is available, so the amplitude will rise as the
offset increasing, the amplitude anomaly on stacked data is apparent.

Anomaly Class 3: Low Acoustic Impedance Sandstone
Sandstone class 3 has lower acoustic impedance than the seal
rock.
Usually this sandstone is less compacted and unconsolidated.
On seismic stack data, sandstone class 3 has big amplitude anomaly
and reflectivity in the whole offset.
Usually, the gradient is significant enough but it has lower
magnitude than the sandstone class 1 and 2.
The RC at normal incidence angle is always negative.
In some conditions, relatively small change of amplitude to offset
can cause detection difficulties because the presence of tuning
thickness, attenuation, recording array, and decreasing of signal-to-
noise to the offset ratio.

Anomaly Class 4:
The fourth class of gas-sandstone anomaly associate with the
reflection coefficient which becoming more positive as the offset
increases, but the magnitude decreased as the offset increases.
Sandstone class 4 often emerged when the porous sandstone,
bounded by lithology which has high seismic wave velocity, such
as hard shale (e.g : siliceous or calcareous), siltstone, tightly
cemented sand, or carbonate.
Sandstone class 4, fall in quadrant II, and has big negative Intercept
(A) and positive Gradient (B).
Sandstone class 4 is the bright spot, but the reflection magnitude
decreases as the offset increases.

Class Relative
Impedance
Quadrant A B Remark
I Higher than the
sealing layer
IV +  Reflection Coefficient (and
magnitude) decreases as the
offset increases
II Almost the
same with the
sealing layer
III or IV   Reflection Magnitude decreases
or increases to offset and the
polarity reversal can happen
III Lower than the
sealing layer
III   Reflection Magnitude increases
to offset
IV Lower than the
sealing layer
II  + Reflection Magnitude decreases
to offset
Reflection coefficient of gas-sandstone top against the
offset for 4 classes of AVO anomaly

AVO classification and exploration risk
AVO
Class
Comparison
to full offset
Risk Impact Remark
I Not much High High Can be generated at major
lithologic boundaries
II enhanced medium highest
III Not much Low Low-to-
moderate
Can be used when
conventional amplitude
analysis is ambiguous
IV nil high low Difficult to produce in a
uniform lithologic section

Gradient (B)
Intercept (A)
“Wet” Trend






 2
V
V
s
p
Base III
Base IIBase II P
Top IV
Top III
Top II Top II P
Base IV
Top I
Base I
Crossplot
Showing
Anomalies
Intercept-gradient cross-plot
Vp /Vs =2

AVO attributes
(1) AVO Product : A*B
(2) Scaled Poisson’s Ratio
Change : A+B
(3) Shear Reflectivity : A-B
(4) Fluid factor
+A
-A
+B
- B
sin2q
The raw A and B attribute volumes are rarely used in that form.
Instead, other AVO attributes are usually calculated from them.

Near Stack Far Stack Fluid factor
AVO attribute Fluid Factor
High fluid factor vale for a envisaged Class-III Gas Sand

AVO attributes λρ, µρ
Modulus Mu (µ) is the capability of rocks to resist
shear strain due to the shear stress
Shear stress changes the shape but retain the volume
µ scaled by Rho (ρ) can indicate the rock type:
 Shear-able rocks such as shale or un-compacted
sands have lower µρ
 whereas rigid rocks such as compacted sands or
tight carbonate associate with bigger µρ

Incompressibility modulus Lambda (λ) scaled with density ρ can indicate the
presence of fluid as it’s value reflects the fluid resistance to the change of
volume due to the compressional stress.
Rocks and fluids can not be compressed easily means that their Lambda is
relatively high.
On the other hand, gas is compressible; therefore the presence of
gas in rocks will decrease its Lambda.
The small λρ indicate a porous rock.
For the same porosity, rocks filled by gas have smaller λρ than
rocks filled by water or oil.
Conversely, tight rocks have bigger λρ.

Vs= √ (µ ∕ρ )
P-Impedance: (ZP )2 = (ρVp)2 = (λ +2µ)ρ
S-Impedance: (Zs )2 = (ρVs)2 =µρ
λρ= (ZP )2 -2 (Zs )2

LMR Example

VP/VS (VP/VS)2 s l+2m m l l/m
Shale 2.25 5.1 0.38 20.37 4.035 12.3 3.1
Gas Sand 1.71 2.9 0.24 18.53 6.314 5.9 0.9
Change /
Average (%)
-27 -55 -45 -9.5 44 -70 -110
Goodway et al. (1997)
Petrophysical analysis indicating that l/mis the most sensitive to
variations in rock properties going from shale into gas sand.
The Power of the LMR Method

Simultaneous Inversion
Simultaneous Inversion inverts for ZP, ZS, and possibly Density using pre-stack
angle gathers as input. The benefit of this procedure is that it allows
constraints to be imposed between these variables. This can stabilize the
results and reduce the non-uniqueness problem.
1 2 3( )PP P S DR c R c R c R   
2
1
2 2
2
2 2 2
3
1 tan
8 sin
1
tan 2 sin
2
S
P
c
c
c
V
V

 
  

 
 
  

1
2
1
2
.
P
P
P
S
S
S
D
V
R
V
V
R
V
R






  
  
 
  
  
 


where:

Simultaneous Inversion
The interpretation of Simultaneous Inversion volumes is similar to
other AVO Inversion results:

Inversion

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.

Acoustic impedance is used to produce more accurate and
detailed structural and stratigraphic interpretations than can be
obtained from seismic (or seismic attribute) interpretation.
In many geological environments acoustic impedance has a
strong relationship to petrophysical properties such as porosity,
lithology, and fluid saturation.
Moreover, the acoustic impedance models are more readily
understood (versus seismic attributes) because it is layer
property unlike the seismic amplitude which is interface
property.
Acoustic Impedance

Cross plot between GR and Impedance in the depth range 3120-
3150 m. Colour bar is coded with respect to lithology (Blue: shales,
red and pink: Sands.
3150 m
3140 m
3130 m
3120 m
Acoustic Impedance-Lithology

Acoustic Impedance-Saturation
Water-sands
Gas-sands
Sg

Crossplots of P-wave velocity and density, color-coded by (a) gamma-ray
measurements
GR

Acoustic Impedance
The definition of the zero-offset reflection coefficient, shown in the
figure above. R0 , the reflection coefficient, is the amplitude of the
seismic peak shown and represents relative impedance contrast.
1122
1122
0
VV
VV
R





Seismic
raypath
Interface at
depth = d
r1 V1
r2 V2
t
Reflection at time
t = 2d/V1
GeologyGeology SeismicSeismic
Surface
Seismic
Wavelet
Shale
Gas Sand

Impedance
Acoustic
Shear Elastic
Reflectivity
Wavelet
Seismic trace
The common forward model for all inversions
General Forward Model for Inversion

ImpedanceReflectivity
Inverse
Wavelet
Seismic
Inversion tries to reverse the forward model
Inverse Model
Acoustic
Shear Elastic

Elastic Impedance
(Patrick Connolly, Bp, TLE-1999)
The elastic impedance (EI) is a generalization of acoustic impedance
for variable incidence angle.
EI provides a consistent and absolute framework to calibrate and
invert nonzero-offset seismic data just as AI does for zero-offset data.
EI, an approximation derived from a linearization of the Zoeppritz
equations.
As might be expected, EI is a function of P-wave velocity, S-wave
velocity, density, and incidence angle.
Elastic impedance is the latest improvement in a continuing process
of integrating seismic and well data for reservoir characterization.

Elastic Impedance
(Patrick Connolly, BP, TLE-1999)
To relate EI to seismic, the stacked data must be some form of
angle stack rather than a constant range of offsets.
An EI log provides an absolute frame of reference and so can
also calibrate the inverted data to any desired rock property
with which it correlates.
EI was initially developed by BP in the early 1990s to help
exploration and development in the Atlantic Margins province,
west of the Shetlands, where Tertiary reservoirs are typified by
class II and class III AVO responses.

Comparison of an AI curve with a 30° EI curve
The 30° elastic-impedance log, EI(30), is broadly similar in appearance
to the acoustic-impedance log
Elastic Impedance

Elastic Impedance
Elastic ImpedanceAcoustic Impedance
Ref. SEG 2000 expanded abstract, Milos Savic, ARCO British Ltd., Bruce VerWest, ARCO British Ltd., Ron
Masters, ARCO Exploration, Arcangelo Sena, ARCO Exploration, and Dean Gingrich, ARCO Alaska Inc
Comparison between acoustic and elastic impedance inversion.
Lithology (GR Log) is better correlated with elastic impedance
GRAI GREI

Elastic Impedance
EI (elastic impedance) is analogous to acoustic impedance, such
that reflectivity can be derived for any incidence angle
R(θ) = (EI2 – EI1)/ (EI2 + EI1)
R(0) = (AI2 – AI1)/ (AI2 + AI1)
K=
At θ = 0 EI = AI = VPρ
Elastic impedance = acoustic impedance

Extended Elastic Impedance (EEI)
It is an extension of the Elastic impedance and
can provide a good approximation to:
 α/β ratio (P-wave velocity / S-wave velocity)
acoustic impedance: αρ
shear impedance: βρ
bulk modulus: K
Lame’s parameter: µ, λ
shear modulus (modulus of rigidity): µ
It can be optimised as a fluid or lithology
discriminator.
Having established a desired value of χ = sin 2
θ an equivalent seismic section can be
constructed using conventional impedance
inversion and AVO processing techniques.
Fluid and lithology impedance can be estimated
K=
E=
σ =

Extended Elastic Impedance (EEI)
Comparison of the optimised EEI curves with the target
Sw and gamma-ray logs

Workflow (Poststack Inversion)
Seismic
Data
Constraints
Geological
Model
Combine and
Invert
Optimum
Section
Final
Inversion

Types of inversion
Type Significance Remark
Recursive Traditional Band Limited Post stack
Sparse Spike Constrained to produce few events Post stack
Coloured Modern derivative of Recursive
Inversion
Post stack
Model Based Iteratively updates a layered
initial model
Post stack
Elastic Enhancement for AVO data Pre-stack
LMR Enhancement for AVO data Pre-stack
Simultaneous Enhancement for AVO data Pre-stack

Band limited recursive inversion
It is simplest and earliest form of Inversion
Starting at the first layer , the impedance of each successive layer
is determined by recursively applying the formula
Ii+1- Ii
Ii+1+ Ii
Ri = 1+ Ri
1- Ri
Ii+1 = Ii
1+ Ri
1- Ri
In = I1 ∏
Recursive inversion produces a result which is band limited to
the same frequency range of the input seismic trace

175
Input Seismic
Recursive Inversion
Recursive Inversion
produces a result
which is
bandlimited to the
same frequency
range as the input
seismic data.
Note the loss of
high frequency
detail, as compared
with the well logs.
Bandlimited Inversion

Issues in Recursive Inversion:
1. The wavelet is ignored. This means that the input seismic
data must be zero phase.
2. Even if the seismic is zero-phase, side-lobes from the
actual wavelet will be interpreted by the algorithm as
lithologic variations.
3. The inversion result is bandlimited to the frequency range
of the seismic data.
4. The scaling of the seismic trace to reflectivity is critical to
get the proper range of impedance changes.
Band limited recursive inversion

Colored inversion
Colored Inversion is modification of Recursive Inversion
In this process, There is a single operator, O, which is applied to
the seismic trace to transform it directly into inversion result:
I= O*S
The operator, O, is defined in the frequency domain
Operator phase is -90 ◦
Amplitude spectra of acoustic impedance of wells are plotted on
log-log scale. The fitted straight line represents “desired” output
impedance. The average seismic spectrum is calculated from
seismic traces around the wells
From two above spectra, the operator spectra is calculate. This has
the effect of shaping the impedance spectrum within seismic band

The amplitude spectrum of the operator is derived this way:
Using a set of wells
from the area, the
amplitude spectra of
the acoustic
impedance for all the
wells are plotted on a
log-log scale.
As predicted by
theory, we can fit a
straight line which
represents the
“desired” output
impedance spectrum.
Log(Frequency)
Log(Impedance)
Amplitude Spectrum of Acoustic Impedance
Colored Inversion

Then, using a set of seismic traces
from around the wells, the average
seismic spectrum is calculated.
Spectrum of Seismic Data
Frequency (Hz)
Operator Spectrum
From the two preceding spectra,
the operator spectrum is
calculated. This has the effect of
shaping the seismic spectrum to
the impedance spectrum within
the seismic band.
Colored Inversion

Input seismic
Colored Inversion
Colored Inversion
produces a result
very similar to
Recursive Inversion.
One difference is
that, in the original
implementation, the
scale is relative
Acoustic Impedance,
with positive and
negative values.
-3000
0
+3000
Colored Inversion

Issues in Colored Inversion:
(1) Very little dependence on the initial model, except to
determine the general impedance trend.
(2) Very fast to apply.
(3) Very simple with few user parameters.
(4) Assumes the data is zero-phase.
(5) In the initial implementation, the method produced a relative
impedance result, although we now have an option to add
back the low frequency trend.
Colored Inversion

Sparse Spike Inversion
It assumes that the actual reflectivity can be thought of a series of
large spikes embedded in the background of small spikes. Only
large spikes are considered meaningful.
Sparse Spike Inversion builds reflectivity sequence one spike at a
time. Spikes are added until the trace is modeled accurately enough.
It produces a broad-band high frequency results.
•Sparse spike Inversion puts events only where the seismic demands
•It attempts to produce simplest possible model consists with the
seismic data
•It often produces fever events than are known to be geologically true
(very thin beds may be missing)
•It may be less dependent on the initial guess model than Model
Based Inversion

Input Seismic
Sparse Spike
Inversion
produces a
broad-band,
high frequency
result.

Constrained Sparse Spike Inversion (CSSI)
The Constrained Sparse Spike Inversion (CSSI) algorithm
produces high quality acoustic impedance volumes from full or
post-stack seismic data. The four output volumes are:
Full Bandwidth Impedance
Band-Limited Impedance
Reflectivity model
Low frequency component
Fugro Jason

Model Based Inversion
It assume that seismic trace, S, and the wavelet, W, are known.
Noise is random and uncorrelated with the signal. It solve for the
reflectivity, R, which satisfies this equation.
It produces a broad-band high frequency results.
Model based Inversion starts with the equation for convolution model
S = W*R + Noise Seismic = wavelet convolved
with reflectivity plus noise
•High frequency details may be coming from initial guess model,
not from the seismic model data
•Errors in the estimated wavelet will cause errors in the inversion
•Seismic does not have to be zero-phase, as long as the wavelet has
the same phase as the seismic
•There is a non-uniqueness problem, as with all inversion.

Step 1:
The initial background model for Model Based Inversion is formed by
blocking an impedance log from a well:
The user specifies the layer size
in milliseconds.
All the layers are originally set
to the same size (in time).

Step 2:
Using the blocked model, and the known wavelet, a synthetic trace is
calculated.
This is compared with the
actual seismic trace.
By analyzing the errors or
“misfit” between synthetic and
real trace, each of the layers is
modified in thickness and
amplitude to reduce the error.
This is repeated through a
series of iterations.
Synthetic Seismic

Input Seismic
Model Based
Inversion
produces a broad-
band, high
frequency result.
The problem is
that the high
frequency detail
may be coming
from the initial
guess model, and
not from the
seismic data.

Differentiating different lithology types through acoustic
impedance inversion

Acoustic impedance
within -20 to 0 ms
window w.r.t.
correlated reflector
showing pay
mapping (gas
saturation) through
inversion. The
channel is filled with
high impedance
shales and not
interesting for
hydrocarbon
accumulation
Channel
Low
impedance
Pay anomaly

Low impedance layer enclosed in high impedance material
Impedanceamplitude

Outline
• Introduction
– Objectives
• Basic concepts
• 4-D Monitoring
• Pitfalls

Time Lapse Reservoir Monitoring (4D)

• Study of two or more 3-D surveys over the same reservoir
for observing changes with time caused by petroleum
production
• Results are often displayed as difference sections or maps
• 4-D surveys are mostly studied for secondary recovery
• Based on acoustic impedance changes due to common
production process
• Easier to monitor
– Steam injection for heavy oil and tar
– Water flood for light oil
– Gas injection
– Oil production involving gas cap expansion

• Critical parameters
– Frame elastic properties of reservoir rocks
– Contrasts in pore fluid compressibility
– Nature of recovery process
– Reservoir parameters
• Depth, pressure, Temperature etc.
Cases of possible high contrasts in pore fluid compressibility
Reservoir fluids Changes from/to Reservoir fluids
Liquid (water, oil) Gas (hydrocarbon gas, steam)
Oil/water Co2 (liquid or gas)
Live oil Water/brine
Oil (live or dead) High salinity brine
Live oil Dead Oil
Low temperature oil High temperature oil

• Good candidate reservoirs for monitoring
– Reservoir with weak rocks
• Unconsolidated or poorly consolidated sands
– Reservoirs undergoing large pore fluid
compressibility changes
• Thermal recovery
– Reservoirs undergoing large rock compressibility
changes
• Thermal recovery, fracturing
– Reservoirs undergoing large temperature changes
• Thermal recovery, water injection

• Time-lapse seismic images can identify bypassed oil to be
targeted for infill drilling, and add major reserves to
production to extend a field’s economic life.
• 4-D seismic can monitor the progress of costly injected fluid
fronts (water, gas, steam,CO2, etc.) and optimizing injection
programs.
• 4-D seismic can map reservoir compartmentalization and the
fluid-flow properties of faults (sealing versus leaking), which
can be extremely useful for optimal design of production
facilities and well paths in complex reservoir flow systems.

• Introduction
– Objectives
• Basic concepts
• 4-D Monitoring
• Pitfalls

What is unconventional gas-shale-systems
• Organic-rich shale formations
– source rock as well as the reservoir
– Very tight porosity and ultralow permeability
• 10 t0 100 nano darcies
• The gas is stored in the form of
– Free gas in limited pore spaces
– Adsorbed gas on the surface of organic
material
• Profitable production requires
– horizontal well drilling
– multistage hydraulic fracturing and stimulation
• Typical shale-gas reservoirs may exhibit
– a net thickness of 15 to 200 m
– porosity of 2–8%
– total organic carbon (TOC) of 1–14%
– depths ranging from 300 m to 4000 m
Core photo of gas shale -
about 50% clay, 50%
quartz plus calcite, 10 -
15% total porosity, 3 - 6%
effective porosity, < 0.001

Inverting Seismic and log data to get rock properties
Basic concepts of Characterization
Lithology inversion
Geology
Elastic parameters
Vp
Vs
Density
Seismic
modeling
Seismic data
Acoustic/elastic inversion
Lithology
Mineralogy
TOC
Porosity
saturation
fracture
Pore structure
permeability
Rock
physics
modeling
The Rock
Travel time
Amplitude
Frequency
Phase
Other pre and
post stack
attributes

Challenges and method in shale-gas exploration
• Exploration stage: finding depth, thickness and
areal extent of promising shales
– Structural and stratigraphic interpretation of seismic
data
• Assessment of source potential:
– Geochemical, sedimentological, log evaluation and
basin modelling methods
– The factors which make good source rocks (TOC,
maturity etc.) also influence the seismic response and
hence can be inferred from seismic attributes
– If shallower sequences are producing oil and gas,
underlying shales (source rock) may have in-situ
hydrocarbon saturation

Challenges and methods in shale gas exploration
• Locating prospective intervals (sweet spots):
distinguishing between reservoir and non-reservoir shales
– Log evaluation methods
– Seismic attributes, Impedance, inversion for log property
volumes, AVO etc.
• Identifying target zones favourable for drilling and
stimulation: mapping natural fractures, stress orientation,
geomechanical properties
– Seismic attributes
– Wide azimuth and multi-component seismic (3D/3C)
• Optimisation of production: activation of fracture
network
– Microseismic imaging
– 4-D seismic surveys

1. Organic-rich sediments have a higher resistivity than organic lean
sediments
2. Organic-rich rocks - increase in sonic transit time and an increase
in resistivity
3. Organic-rich rocks can be relatively highly radioactive
(higher gamma-ray reading than ordinary shales and limestones)
4. The ΔLogR cross plot is a quick and easy determination of
resource
Log response of Shale-gas reservoir
Barnett Thin Section and Log Rock Types
1.Black shales 1.Clay-rich; separation of neutron-density logs
2.Calcareous black shale 2.High Gamma Ray –130-140 API
3.Phosphaticblack shales 3.Highest Gamma Ray –150+ API
4.Limey grainstones 4.Lower GR 100-120
5.Dolomitic black shales 5. Lower GR 100-120; neutron and density stack
or cross over

DT/LLD GR RHOB/NPHI
130
615
DT (µs/m)
1.0 100.0
LLD (Ώm)
source
interval
Dominantly sands
Dominantlyshales
Source shale
Depth (m)
Possible source interval identified from gamma-ray (GR), density (RHOB), Neutron
porosity( NPHI), transit time (DT), and relatively high resistivity (LLD)
Shale
Sand
Coal
source
DT-Resistivity overlay method of source rock identification

Blue: Low DT (high velocity/ low resistivity, representing lean shales
Green-yellow-red: Medium DT (medium velocity)/ medium resistivity shale-
reservoirs
Pink: high DT (low velocity)/high resistivity coals.
DT (Sonic) Map Resistivity Map
Sweet Spot Thickness about 12 m

• Introduction
– Objectives
• Basic concepts
• Pitfalls

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

• 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
Pitfalls in Seismic Interpretation

DHI
Seismic section with
dry hole and good well
penetrate high
amplitudes, but only
the latter
with red-over blue
(trough-over-peak) is
representative of low-
impedance prospective
sand. The
blue-over-red
amplitude represents a
poor prospect.

DHI-fizz water
Typical effect of gas saturation on P-velocity of
rocks under shallow conditions.

The difference in amplitude behavior between the near- and far-trace
stacks in (a) implies the changes in elastic properties shown in (b).
These properties are not known with certainty, but even if they were,
there might be multiple geologic models (c).
Hypothetical example of Pitfalls in DHI

The difference in
amplitude behavior
between the near- and
far-trace stacks in (a)
implies the changes in
elastic properties shown
in (b). These properties
are not known with
certainty, but even if they
were, there might be
multiple geologic models
(c) that have those
properties. All models
explain the observation,
but only one contains
hydrocarbons.
C. Possible
geologic models.
All models explain
the observation, but
only one contains the
hydrocarbons.
Hypothetical example of Pitfalls in DHI
DHI: Low impedance sands
Pitfall: Low gas saturation
Pitfall: anomalous lithologic variation

Bp sesmic interpretation

More Related Content

What's hot (20)

Similar to Bp sesmic interpretation (20)

Recently uploaded (20)

Bp sesmic interpretation