Transcript Overview of Seismic Interferometry and Applications in

Seismic Interferometry: Instead of using just primary
arrivals, you also use the multiples for a wider view
Overview of Seismic Interferometry
and Applications in Exploration
Gerard Schuster
KAUST & University of Utah
Outline
•What is Seismic Interferometry?
•Applications
•VSP->SSP (surface seismic profile)
•VSP->SWP (single well profile)
•SSP->SSP
•Conclusions
SELECTIVE HISTORY SEISMIC INTERFEROMETRY
1968
!
Claerbout V(z)+passive
redatum
1970s Berryhill model-based redatum
1980s Cole+Claerbout V(x,y,z)+passive?
1990s Scherbaum earthquake V(z)+passive
1999
Rickett+Claerbout V(z) Helioseismology
Daylight Imaging, passive
Utah: Stationary Phase Theory, SSP, and VSP
Seismic Interferometric imaging, deterministic
2002-04 Wapenaar Recip. Thm. Correlation Type
Gerstoft + others Surface Wave Interferometry
Snieder Stationary Phase Redatuming
Shell Virtual Sources:Calvert+Bakulin
2001
SELECTIVE HISTORY SEISMIC INTERFEROMETRY
!
redatum
Passive Reservoir
Earthquakes
Shell, Draganov,
Wapenaar, Snieder, Polleto
Miranda, etc
Nowack, Sheng, Curtis etc
Engineering Xwell
Minato, Onishi, Matsuoka etc
Surface waves
Shapiro, Derode, Larose, Dong,
Xue, Halliday, Curtis,
Van Mannen, Robertsson,
Volcanoes+Coda
Snieder, Scales, Gret et al
Gerstoft, Sabra, Kepler, Roux,
He, Ritzwoller, Campillo etc
Model Tank
Scales, Malcolm etc
Interpolation
Sheng, Curry, Berkhout, Wang,
Dong, Hanafy, Cao, etc
VSP
Yu, Calvert, Bakulin, He,
Jiang, Hornby, Xiao, Willis, Lu,
Toksoz, Campman etc
Extrapolation
Dong, Hanafy, Cao, etc
EM
Slob, Wapenaar, Snieder
Theory: Acoustic, EM, Elastic, Potential
Fink, Wapenaar, Snieder, Papanicolaou, Blomgren,
Slob, Thorbeck, van der Neut etc
Refractions
Boise State Univ, Dong
Exploration
Curry, Guitton, Shragg, Yu,
Artman
What is Seismic Interferometry?
Answer: Redatums data by correlation of trace pairs and
stacking the result for different shot positions
iwtxB iw(t xB+ tBz+ tzB)
e
e
s
*
G(B|x) G(B|x)
VSP
Assume a
VSP experiment
z
Phase of Common
Raypath Cancels
x
= e
G(B|B)
=>
B
Point Source Response
with src at B and rec at B
SSP
F.S. multiple
direct
B
=
iw( t Bz+ tzB )
virtual
primary
virtual
source
z
A
No need to know src. location
z•
•No need to know src excitation time
•Redatum source closer to target
What is Seismic Interferometry?
Answer: Redatums data by correlation
correlation of trace pairs and
stacking the result for different shot positions
stacking
iwtxB iw(t xB+ tBz+ tzB)
x
x
x
e
e
G(B|x)* G(B|x)
~
=
~
z
z•
Phase of Common
Raypath Cancels
iw( t Bz+ tzB )
= e
G(B|B)
Point Source Response
with src at B and rec at B
z
A
No need to know src. location
•No need to know src excitation time
•Redatum source closer to target
Reciprocity Correlation Equation
2D Reflection Data
x
k =~
G(x|B)* G(x|A)
~
SSP
VSP
VSP
B
G(A|B)
B
A
A
x
x
A
Old Multiples Become
New Primaries!
•No need to know VSP rec location at x Phase of Common Raypath Cancels
•No need to know receiver statics
Reciprocity Correlation Equation
2D Reflection Data
{G(B|x)*
k =* G(A|B)
2
G(x|B)*
G(x|A)
*
n
G(A|x)
d
x
=
G(A|B)
G(B|A)
}
G(A|x)
G(B|x)
x
(Wapenaar, 2004)
S
well
Finite aperture leads to incomplete G(B|A)
B
B
A
A
x
Problems: Finite source aperture
No attenuation
1-way+ far-field approx.
x
A
Old Multiples Become
New Primaries!
Muting, Least squares or MDD
Atten. Compensation
Summary
•Seismic Interferometry:
x
G(x|B)* G(x|A)
~
~
imaginary
Im[G(A|B)]
x
x
A
G(A|x)
B
G(B|x)
k
A
B
G(A|B)
•Merits: Eliminates need for src location, excitation time,
some statics. Moves rec./srcs closer to target , no velocity
model needed (unlike Berryhill).
•Challenges: Finite aperture and noise, attenuation,
acoustic & farfield approximations
Outline
•What is Seismic Interferometry?
•Applications
•VSP->SSP (surface seismic profile)
•VSP->SWP (single well profile)
•SSP->SSP
•Conclusions
Implementation
VSP
k
x
VSP
SSP
G(A|x)* G(B|x) =
Im[G(A|B)]
1. FK Filter up and downgoing waves
2. Correlation: f(A,B,x) = G(A|x)* G(B|x)
3. Summation: k
f(A,B,x) =
Im[G(A|B)]
x
4. Migration: M(x) = Mig(G(A|B))
A B
A B
x
A B
x
Challenge: Finite Receiver Aperture = Partial Reconstruction
x
3D SEG Salt Model Test
(He, 2006)
VSP Multiples Migration
Stack of 6 receiver gathers
( Courtesy of P/GSI: ~¼ million traces, ~3 GB memory, ~4 hours on a PC )
(He, 2006)
Marine 3D VSP Field Data
Application
BP 3D VSP Survey Geometry (36 recs)
~ 11 km
1.6 km
4.0 km
3 km
(He et al., 2007)
VSP->SSP Summary
VSP
k
VSP
SSP
!
x
G(A|x)* G(B|x) =
A B
Im[G(A|B)]
A B
x
A B
x
x
Key Point #1: Every Bounce Pt on Surface Acts a New Virtual Source
Key Point #2: Kills Receiver Statics
Key Point #3: Redatuming = Huge Increase Illumination area
Key Point #4: Liabilities: Finite Aperture noise, attenuation, loss amplitudes fidelity
Outline
•What is Seismic Interferometry?
•Applications
•VSP->SSP (surface seismic profile)
•VSP->SWP (single well profile)
•SSP->SSP
•Conclusions
Motivation
Problem: Overburden+statics defocus VSP migration
Solution: VSP -> SWP Transform (Calvert, Bakulin)
VSP
VSP
SWP
Redatum sources below overburden
Local VSP migration
VSP Geometry
1500
Reflection
wavefield
Depth
(m)
0
3500
0
Offset (m)
Time (s)
1000
(He , 2006)
3
VSP Geometry
1500
Reflection
wavefield
Depth
(m)
superresolution
0
3500
0
China
Offset (m)
Time (s)
1000
(He , 2006)
3
VSP Salt Flank Imaging
(Hornby & Yu, 2006)
120 shots
Overburden
?
98 geophones
Poor image of flank by standard migration
Interferometric Migration
Result
0
2000 ft
VSP->SWP Summary
!
1. Redatum sources below overburden
2. Local VSP migration
3. Kills Source Statics and no need to know src location or excitation time
4. Super-resolution
5. Instead of redatuming receivers to surface, we
redatum sources to depth.
Outline
•What is Seismic Interferometry?
•Applications
•VSP->SSP (surface seismic profile)
•VSP->SWP (single well profile)
•SSP->SSP
•Conclusions
Surface Wave Interferometry
G(A|x)*
G(B|x)
G(B|A)
x
x
A
BB
Surface Wave Interferometry
G(A|x)* G(B|x)
= G(B|A)
x
A
B
Surface Wave Interferometry
S-velocity distribution, surface wave predic.+elimination
x
G(A|x)* G(B|x)
= G(B|A)
x
A
B
Shear velocity
Yao (2009)
3x3 Classification Matrix
out
in
SSP
VSP
SWP
SSP SSP
SSP
SSP
VSP
SSP
SWP
VSP VSP
SSP
VSP
VSP
VSP
SWP
SWP SWP
SSP
SWP
VSP
SWP
SWP
Summary
•Seismic Interferometry:
x
G(x|B)* G(x|A)
~
~
Im[G(A|B)]
x
x
A
G(A|x)
B
G(B|x)
k
A
B
G(A|B)
•Merits: Eliminates need for src location, excitation time, some statics.
Moves rec./srcs closer to target , no velocity model needed (unlike
Berryhill).
•Challenges: Finite aperture and noise, attenuation, acoustic & farfield
approximations , amplitude fidelity
•Killer Apps in Earthquake: Surface wave interferometry
•Killer Apps in Exploration: Passive reservoir monitoring? OBS? EM? VSP
Thanks
•UTAM sponsors
•Min Zhou, Chaiwoot Boonyasiriwat, Ge Zhan