Transcript Summary of SEP research

Overview of SEP research
Paul Sava
The problem
Modeling
operator
Lm  d
Seismic
image
Seismic
data
Migration
Migration
operator
mLd
*
Seismic
image
Seismic
data
Migration operator: wavefield propagation
mLd
*
• Downward continuation
– Common-azimuth migration
– Narrow-azimuth migration
• Reverse time migration
•
•
•
•
Propagation
Imaging
Amplitudes
Velocity
Downward continuation
Biondi, 2003
Commonazimuth
Narrowazimuth
Reverse-time migration
Biondi & Shan, 2002
Reverse-time
migration
Downwardcontinuation
Migration operator: angle-gathers
mLd
*
•
•
•
•
Propagation
Imaging
Amplitudes
Velocity
• Data-space
• Prucha et al
• Image-space
•
•
•
•
Sava&Fomel, Tisserant&Biondi (S-G)
Rickett&Sava (shot profile)
Rosales&Rickett (converted waves)
Biondi&Shan (reverse-time)
Angle-domain common image gathers
Prucha et. Al (1999), Sava&Fomel (2002)
Image space
ADCIG
Data-space
ADCIG
3-D angle gathers
Tisserant & Biondi (2003)
Migration operator: amplitudes
mLd
*
•
•
•
•
Propagation
Imaging
Amplitudes
Velocity
• Amplitude preserving wavefield extrapolation
• Sava & Biondi
• Amplitude corrections of extrapolation operators
• Vlad et. al
• Hi resolution imaging condition
• Valenciano & Biondi
Amplitude-preserving migration
Sava & Biondi (2002)
Kinematic
migration
Amplitude
preserving
Migration operator: velocity
mLd
*
• Traveltime-based
• Clapp
• Wavefield-based
• Sava & Biondi
•
•
•
•
Propagation
Imaging
Amplitudes
Velocity
Angle-domain traveltime tomography
Clapp (2001)
Tomography: geological constraints
Clapp (2001)
1
-.4
-.6
Wave-equation MVA
Sava & Biondi (2003)
g
x
z
g
x
z
Slowness
perturbation
Image
perturbation
Wave-equation MVA
Sava & Biondi (2003)
Multiple attenuation
mLd
*
• Multiple attenuation
– Data space
– Image space
• Multiple imaging
• Joint imaging
Multiple attenuation: data space
Guitton (2003)
Input data
Time (s)
2.6
5.
Hyperbolic
Radon
Adaptive
filtering
Pattern
recognition
Multiple attenuation: image space
Sava & Guitton (2003)
PRT
image
space
HRT
data
space
Multiple imaging
mLd
*
• Multiple attenuation
– Data space
– Image space
• Multiple imaging
• Joint imaging
Multiple imaging: shot-profile migration
Guitton (2002)
Down-going
wavefield
Impulse
Up-going
wavefield
Down-going
wavefield
Up-going
wavefield
Primaries
Primaries
Multiples
Primaries imaging
Multiples imaging
Multiple imaging: S-G migration
Shan (2003)
xD
xs
xD
xs
xU
h
x
h
t1
R1
R1
R2
t2
t2
R2
xU
Joint imaging
mLd
*
• Multiple attenuation
– Data space
– Image space
• Multiple imaging
• Joint imaging
Joint imaging
Brown (2003)
S
G
Migration
Migration
operator
mLd
*
Seismic
image
Seismic
data
Least-squares imaging
 
Inversion
operator
1
m LL Ld
*
Seismic
image
*
Seismic
data
Least-squares imaging
m  WL d
*
• Illumination compensation
• Rickett (2001)
• Prucha (2003)
•
•
•
•
Normalized migration
Matching filters
Least-squares inverse
Multiple realizations
Normalized migration
Rickett (2003)
Least-squares imaging
Guitton (2003)
•
•
•
•
m  BL d
*
L*
data
Normalized migration
Matching filters
Least-squares inverse
Multiple realizations
L*L
m1
m2
(L*L)-1
B
Find B such that
Bm2  m1
" mtrue"
Imaging with mathching filters
Guitton (2003)
Least-squares inverse
 
1
•
•
•
•
m LL Ld
*

*
Normalized migration
Matching filters
Least-squares inverse
Multiple realizations

1
m  LLε A A Ld
*
2
*
*
Regularization
Prucha (2003)
ph
z
z-ph
preconditioning
x
z
z-x
preconditioning
Regularized inversion
Prucha (2003)
Least-squares
inverse
Migration
Multiple realizations
 
1
•
•
•
•
m LL Ld
*

*
Normalized migration
Matching filters
Least-squares inverse
Multiple realizations

1
m  LLε A A Ld
*
2
*
Lm  d  0
Am  0
*
Multiple realizations: interpolation
Clapp (2002)
Lm  d  0
Am  0
Lm  d  0
Am  n
Multiple realizations: velocity
Clapp (2002)
http://sepwww.stanford.edu
• Reports (all online)
• Seplib
• Computers: 4 clusters
• 3D real data