Session 2 Rock Physics Analysis - Pre-Upscaling For Seismic
Session 2 Rock Physics Analysis - Pre-Upscaling For Seismic
Session 2 Rock Physics Analysis - Pre-Upscaling For Seismic
ROCK PHYSICS
Study of the physical properties and their relationships of sedimentary rocks.
Specifically , in terms of how seismic properties of velocity, impedance and
attenuation behave and relate to other petrophysical properties (elastic
behaviour) and reservoir parameters.
Advent of the LSS and ARRAY SONIC has enabled the FULL
WAVEFORM to be recorded.
• ARRAY SONIC designed to analyse complete composite
waveform by combining an array of receivers and dual
transmitters/receivers. The tool combines the real-time
measurement of the Compressional, Stonely and Mud Wave
Velocities from which the Vs is derived.
ROCK PHYSICS THEORY
In the beginning ……
which
which ultimately
P transmitted
P refracted
P reflected S transmitted
S reflected
Incident Energy
Layer 1 Layer 2
SHEAR WAVES Wavefront that is generated as a result of the primary
wave front in a direction that is perpendicular to the
Note the beautiful shear
direction of particle displacemnt, and generally has
spray pattern a velocity 2/3 of Vp for a given medium.
Vp/Vs= 1.6 for clean porous media.
P transmitted
P refracted
S transmitted
P reflected
S reflected
Incident Energy
Layer 1 Layer 2
Seismic Reflection schematic
Incident energy S reflected
1 0.1 P reflected
Layer 1
0.5
Vp, Vs, s2, r2 0.2
P transmitted
S transmitted
Seismic Reflection schematic
Incident energy S reflected
1 0.05 P reflected
Layer 1
0.4
Vp, Vs, s2, r2 ... 0.25
P transmitted
S transmitted
ELASTIC ROCK PROPERTIES (Mechanical & Seismic)
Comp Stress Comp Strain
Index of Confinement
Bulk Modulus
Comp Stress
K = C Stress
C Strain
K = RHOB.(Vp2-4/3Vs2)
Relates change in Volume with Stress-Strain Ratio under simple
hydrostatic pressure
Shear Stress Shear -Lateral
Strain
Index of Distortion
Shear Modulus
G = S Stress
S Strain
G = RHOB.Vs2
Normal Strain
The amount by which a material
increases along one axis when
Radial Strain
compressed along the other.
SUMMARY of ROCK PHYSICS KEY ELEMENTS
Help bridging the gap between Rock Physics and Seismic is the fact that Elastic
Moduli do not have to be in terms of Mechanical Properties (of stress and strain)
but
measured in
Travel time frames are governed by the Mechanical Properties of the Rock
where
Porosity vs P-Impedance
12
P - IMPEDANCE (km/s.g/cc)
10
Legend
AWate r
8
AClst/Sh
AA/A Clst/s
AA W ater
6
AA G as
AA Clst/Sh
4 A Sand
0 5 10 15 20 25 30 35
Porosity - (%)
P-Impedance vs VpVs
21
Legend
20
DingoClst
19 Dingo/Brig?
x10 -1
Barrow G p
18
VP_VS
AWate r
AClst/Sh
17
AA/A Clst/s
16 AA W ater
AA G as
15
AA Clst/Sh
14 A Sand
4 5 6 7 8 9 10 11
.
P-IMPEDANCE
EFFECTS of LITHOLOGY and TEXTURE
13
12
11 B Gp Slt/C ls t+
D in go C ls t
10
9
Impedance
8
D u puy /D ingo
7
B Gp Sa nds Ma rdie
6 Mu dero ng
5
W inda lia
4 Legend
3 W etSn d
2 c ls t/s ls t
0.00 .05 .10 .15 .20 .25 .30 .35 .40 .45 .50
Poisson's Ratio
F igu re 15
EFFECTS of LITHOLOGY and PORE FLUIDS
DTVPVS vs DTCO
vs VpVs Ratio
2.30
2.20
2.10 O 'ba nk
2.00
C r av Sply
1.90
VPVS
1.80
W et Ss t
1.70
Carbonates
1.60
S lty Sn dGa s
Legend
1.50
C ln Ga s W etSn d
1.40
Ga s Sn d
1.30 c ls t/s ls t
40 50 60 70 80 90 10 0 11 0 12 0 13 0 14 0
DT CO (us/ft)
EFFECTS of BULK DENSITY
RHOB vs Vp
60 00
55 00 SAND
DingoClst
50 00
Dingo/Brig?
45 00 Barrow G p
VP m/s
AWate r
40 00
AClst/Sh
35 00 AA/A Clst/s
AA W ater
30 00
AA G as
25 00
AA Clst/Sh
20 00 A Sand
1.8 0 1.9 0 2.0 0 2.1 0 2.2 0 2.3 0 2.4 0 2.5 0 2.6 0 2.7 0 2.8 0 2.9 0 3.0 0
RHOB - g/cc
In general, seismic properties increase with bulk density, although wide scatter
can occur.
EFFECT of POROSITY
Velocity Vp vs Vs
4000
3500
3000
2500 MungarooW et
Vs - m/s
2000 MungarooG as
MungarooO 'B nk
1500 Dingo
1000
Legend
W etSnd
500
GasSnd
0 cls t/slst
2000 2200 2400 2600 2800 3000 3200 3400 3600 3800 4000 4200 4400 4600 4800 5000
Vp - m/s
SUMMARY of ROCK PHYSICS PROPERTIES
hence
Unconsolidated high porosity rocks are less rigid
& more compressible
than
Consolidated low porosity rocks
2 2
Sand RHOB Vp(dsi) Vs(dsi) Vs Vp Vp/Vs Impedance PR ShearMod BulkMod Young'sMod
Unit g/c
3
m/s ft/s m/s ft/s km/s km/s m/s
3
g/c /m/s v/v m- K- E-
3 3 3
m.g/s.cm m.g/s.cm m.g/s.cm
Soft Clay 2.35 1100 3608.9 550 1804.5 0.30 1.21 2.00 2.59 0.33 0.71 1.90 1.90
Clay 2.5 2100 6889.8 900 2952.8 0.81 4.41 2.33 5.25 0.39 2.03 8.33 5.62
Loose Sand 2.65 2743 8999.3 1714 5623.4 2.94 7.52 1.60 7.27 0.18 7.79 9.56 18.37
Sand 2.66 3505 11499.3 2100 6889.8 4.41 12.29 1.67 9.32 0.22 11.73 17.04 28.62
Limestone 2.71 3962 12998.7 2133 6998.0 4.55 15.70 1.86 10.74 0.30 12.33 26.10 31.96
Dolomites 2.87 4572 15000.0 2286 7500.0 5.23 20.90 2.00 13.12 0.33 15.00 40.00 39.99
MarineShale 3 5180 17000.0 2590 8497.4 6.71 26.83 2.00 15.54 0.33 20.12 53.67 53.67
Granite 2.67 5640 18503.9 2870 9416.0 8.24 31.81 1.97 15.06 0.33 21.99 55.62 58.29
FLUID SUBSTITUTION PLOT
Input Parameters
TDDepth Mineral Volume Fractions
Res Depth Qtz Fluid Properties
FmTemp Calcite Salinity Previous Fluid Output
GG Dolomite Oil Gravity Brine Vp, R, K Previous FluidMix
Fm Press Clay GOR Oil Vp, R, K & SolidOutput
Hmin Gas Grav Gas Vp, R, K & Fluid Vp, R, K
Gas Index in Brine Log Sw Solid R, K, G &
Gas Index in Oil Log So Log Porosity
Pore Press & Temp Log Sg
Output
Solid R Brine Vp, R, K Fluid Vp Rock Vp, R, K, G
Solid Vp Oil Vp, R, K Fluid R RockVs
Solid Vs Gas Vp, R, K Fluid K RockPR
Solid K
Solid G
Solid PR
EFFECTS of VALID SHEAR
EFFECT of ERRONEOUS SHEAR
POISSON’S RATIO vs Vs SENSITIVITIES
230
210
190 Dtp130
Unconsol
Sand
D
170
Delta T Shear - m/s
D
D
150
Sand
D
D
130
Silt D
D
110 Shale
D
90 Dtp60
Gas Sand Water Sand Silt/Shale Sand Siltstone/Shale
70
50
0.0000 0.0500 0.1000 0.1500 0.2000 0.2500 0.3000 0.3500 0.4000 0.4500 0.5000
Poissons Ratio
VELOCITY RATIO vs Vs SENSITIVITIES
Idealised Vp/Vs vs Dt Shear for range of Dtp
260
240 Unconsol Sand Silt
Shale
Sand
220
200
Dtp130
180
DtShear - m/s
160
140
120
100
80 Dtp60
Clean Coarse to Shaley Siltstone
60 Fine Sand Sand /Shale
40
20
0
1.4 1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 2.4
Vp / Vs Ratio
ROCK PHYSICS PLOT
ROCK PHYSICS PLOT
ROCK PHYSICS PLOT
Mungaroo AA Sand
Marginal MarineSequence
Fluvial Sequence
GWC @ 3280.0mRT
ROCK PHYSICS PLOT
Reference Aquifer
Sand Anomaly
Gas Anomaly
Stack
Gradient
AI-S PLOT
Reference Aquifer
(ref values on the stack)
Reservoir
Anomaly
Fluid and porosity sweep on A1-S plot
Rx PROPERTY UPSCALING WORKFLOW REVIEW
Other effects not covered by this discussion
ANISOTROPY
were isotropic rocks have 2 independent Elastic Constants
INVERSION
inversion for rock and fluid properties via multi-attibute
analysis / acoustic impedance inversion
CONCLUSIONS
thereby
Raw Curves
Gamma Ray
Deep, Medium, Shallow Induction Resistivity
Density, Neutron & PEF Curves
Compressional Sonic
OBDT- Oil Based Dipmeter
VSP
Processed Curves
STC Processing DSI - Dipole Shear, PR & VS/VP
Startigraphic / Structural Dips
VSP Processing
PROCESSED LOG INPUT required for ROCK PHYSICS
MAIN OUTPUTS
. Total Porosity
. Effective Porosity
. Water Saturation
. MobileSaturation
. Immobile Saturation
. Invasion Profile
. Sandstone Volume
. ShaleVolume
. Limestone Volume
. Dolomite Volume
. Heavy Mineral Volume
. Net Sand Flag
. Net Pay Flag
MAIN LOG OUTPUT for ROCK PHYSICS
MAIN OUTPUTS
.Velocity/Sonic
-Compressional
-Shear
. Elasstic Moduli
-Bulk Modulus
-Shear Modulua
-Poisson’s Ratio
-AcousticImpedance
EFFECTS of LITHOLOGY and PORE FLUIDS
Porosity vs P-Impedance
12
P - IMPEDANCE (km/s.g/cc)
10
Legend
AWate r
8
AClst/Sh
AA/A Clst/s
AA W ater
6
AA G as
AA Clst/Sh
4 A Sand
0 5 10 15 20 25 30 35
Porosity - (%)
P-Impedance vs VpVs
21
Legend
20
DingoClst
19 Dingo/Brig?
x10 -1
Barrow G p
18
VP_VS
AWate r
AClst/Sh
17
AA/A Clst/s
16 AA W ater
AA G as
15
AA Clst/Sh
14 A Sand
4 5 6 7 8 9 10 11
.
P-IMPEDANCE
EFFECTS of LITHOLOGY and PORE FLUIDS
VPVS vs DTCO
2.30
2.20
2.10 O 'ba nk
2.00
C r av Sply
1.90
VPVS
1.80
W et Ss t
1.70
1.60
S lty Sn dGa s
Legend
1.50
C ln Ga s W etSn d
1.40
Ga s Sn d
1.30 c ls t/s ls t
40 50 60 70 80 90 10 0 11 0 12 0 13 0 14 0
DT CO (us/ft)
EFFECTS of LITHOLOGY & PORE FLUIDS
14
12
Acoustic Impedance
10
D in goC ls t
6
Legend
Mu ngar ooG as
W etSn d
4
Ga s Sn d
2 c ls t/s ls t
0.00 .05 .10 .15 .20 .25 .30 .35 .40 .45 .50
PRmonopole
KB:27.4m WD 1100m f :\xl\f lis...JLRoche
MODELLING SHEAR WAVES
Rock consists of
– Mineral Components
– Porosity
– Fluid
• Each constituent contributes to the Mechanical behaviour of the
rock as a whole
ROCK consists of Solid Rock Constituents
QTZ
CLAY
LST Vp, Vs, R, K, G, PR
DOL
HMIN
POROSITY & SAT Effected by Por & Sw
FLUID BRINE Effected by Salinity,
OIL PP, Tmp, GG, GOR,
GAS Oil Grav, Gas Grav, Sw
THE UPSCALE PROCESS -FLUID PROPERTIES SETUP
Answer
2 main factors are:
Rock Properties (Lithology)
Presence of hydrocarbons (Fluid)
Amplitude Analysis WorkFlow
Modelling Seismic
Arrive at interpretation
Data Summary for UpScale Modelling
• Compressional velocity (Vp)*
• Shear velocity (Vs)
• Density (r) *
• Poisson’s ratio (s)
• Saturation (%)
• Clay content (%)
• Thickness (m, ft)
• Depth (m, ft)
• Fluid properties (salinity, GOR, API…)
* Note: if some info is missing, there is a large rock properties database from which
probable values can derived
The Gassmann Equation
(1 - K d / K m ) 2
K* = Kd +
1- Kd
+ - 2
Kf Km Km
G* = Gd r * = r d + r f
K=bulk modulus : Km- solid rx; Kd-dry frame rx; Kf-fluid; K*-saturated rx
G=shear modulus : Gm- solid rx; Gd-dry frame rx; Gf-fluid; G*-saturated rx