Transcript Document
Prestack depth migration in angle-domain using beamlet decomposition: Local image matrix and local AVA
Ru-Shan Wu and Ling Chen
Modeling and Imaging Laboratory/IGPP University of California, Santa Cruz ------------------------------------------------------ † Presently at Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing, China
Outline
Beamlet decomposition: Wave field in angle-domain Local image matrix and local scattering matrix Effect of acquisition aperture Local AVA: Preliminary tests Conclusion
True-reflection imaging in angle-domain
• • • Preserving the relative amplitudes of scattered waves w.r.t. incident waves.
Benefits:
Improve image (
total strength image)
quality, especially for steep reflectors.
Reduce artifacts
(angle-domain filtering)
.
Provide basis for local AVA and local inversion
True-amplitude imaging in angle-domain
• • • • • Amplitude corrections ( for ray theory see Hubral et al., Bleistein et al., Xu et al., Audebert et al., ……):
Transmission loss
( boundary reflection and scattering)
Geometric spreading
( for ray method )
Nonuniform information distribution
: Jacobian ( Beylkin determinant )
Acquisition aperture effects
( hit-count for ray method)
Intrinsic attenuation
( Anelasticity )
True-reflection imaging in angle-domain for wave-equation based methods
• • • Preserving the relative amplitudes of scattered waves w.r.t. incident waves:
Nonuniform information distribution
: Jacobian
Acquisition aperture effects
(in angle-domain) (including the geometric spreading and hit-count for ray method)
Transmission and anelastic losses
are less important, especially for small-angle reflections
Beamlet decomposition of the wavefield:
u
u z
n m
n m u z u
,
x n
~
b mn
,
m b mn
b mn
where ~
b mn
(
x
,
z
) ----- decomposition vectors (atoms),
u b mn
z x
(
n x
, ,
z m
) ----- reconstruction vectors (atoms), ----- coefficients of the decomposition atoms,
x
n m
n
m
x
----- window location, ----- local wavenumber .
Windowed plane waves
G-D frame atoms
g mn
e i
m x g
x
x n
is a Windowed Plane wave
(each beamlet is a windowed plane wave)
Local plane waves
Local plane wave
: a
superposition of windowed plane waves
of the same local wavenumber from all neighboring windows : Partial reconstruction of wavefield (mixed domain wavefield: local phase–space):
u
x
,
z
,
j
,
e i
j x
l g
x
The corresponding propagating angle:
x l
u z
x l
,
j
,
j
cos 1
j
v
Source
**
Receiver High-velocity body Local Scattering Matrix
S
(
in
,
sc
)
sc
Target area
in
Local Image Matrix (includes aperture and propagation effects)
L
(
in
,
sc
)
shallow Point scatterer Planar reflector deep
Local image matrix in homogeneous medium (total 201 shots with 176 left-hand receivers )
Local image matrix:
image condition in
beamlet domain
and
mixed domain
Forward-propagated source field:
u S
x
,
z
,
j l u z S
x l
,
j
,
g jl
Backward-propagated scattered field:
u R S
x
,
z
,
p q u z R S
x q
,
p
,
g pq
Local image matrix
:
L a
in , sc ,
x
,
z
, 0 2 cos in cos sc
S
W a S
x
, in ,
z
,
R S W a R S
x
, sc ,
z
, Where
k
0 2 cos in cos sc Serves as the Jacobian W s and W Rs are the wave fields in angle-domain by beamlet decomposition
Stacking over frequency to get the final image
I
In the local angle-domain: in , sc , Re
d
e i v
sin in sin sc
x L a
in , sc ,
x
,
z
, The final image in space domain:
I
in sc
I
in , sc ,
x
,
z
Local Reflection Analysis (AVA):
For planar reflectors:
the local image matrix can be represented as:
I
n , r ,
x
,
z
with
n
r
in in
sc sc
2 2
CDAI (common dip-angle image) gathers Sum up all reflections for a common dip angle:
I
CDAI
n
,
x
,
gather
z
r I
n
,
r
,
x
,
z
Obtain the dip-angle of the local reflectors from CDAI.
CRAI (common reflection-angle image) gathers Sum up reflected energy for a common reflection-angle for all possible dip-angles: CRAI.
Performing local AVA from CRAI gathers.
The calculation of local reflection coefficients:
R
r
,
x
,
z
1
N
n I
n
,
r
,
x
,
z
Local AVA for an oblique interface in homogeneous background
Local image matrices
at a point on the middle of dipping interface 14 ° obtained from 80 shots with a two-side receiver array (513 receivers). The dotted line corresponds to the theoretical prediction without aperture effect.
CDAI gathers for a local reflector at its central location Obtain the dip-angle of local reflectors from the CDAI gathers
Calculated reflection coefficients from CRAI gathers .
Angle-dependent reflection coefficients at the interface using 256 shots with 513 two-side detectors for the horizontal layered model with different velocity contrasts: (a) 10%; (b) 25%; (c) 50%; (d)150%
Dotted:
synthetic;
Red:
513 points two-sides
Blue:
257 points one-side;
Green:
129 points two-sides
Angle-dependent reflection coefficients at the interface obtained from LIM in case of 10% velocity contrast for the horizontal layered model
Local image matrix and the local scattering matrix
The
local image matrix
has the acquisition-aperture and propagation effects included. The purpose of the imaging/inversion is to recover the
real local scattering matrix coefficients
.
and obtain the local reflection
To achieve the true-reflection imaging, we need to
estimate the acquisition-aperture effect and apply the correction
.
Acquisition-Aperture Efficacy
(Effect of the source-receiver configuration) • • • Acquisition-Aperture Efficacy (AAE) Matrix Acquisition-Aperture Dip-response function Aperture corrections
Source
**
Receiver Overburden structures Assume scattering Coefficients as 1
sc
Target area
in
Acquisition-Aperture Efficacy: (includes propagation effects)
E
(
in
,
sc
)
Acquisition-aperture efficacy matrix
With unit impulses at both source and receivers, the local
acquisition-aperture efficacy matrix
is obtained as
:
E z
j
,
p
,
x
,
S l z
x l
,
j
,
x S
,
g
x
x l
2
q R S z
*
x q
,
p
,
x R S
,
g
x
x q
2 1 2 Where G’s are the Green’s functions in beamlet domain
Acquisition-aperture dip-response function
Acquisition-aperture dip-response
as a function of dip-angle of local interface (reflector), which reduce the AAE matrix into a vector:
E
z
n
,
x
,
r E z
n
,
r
,
x
, with
n
r
in in
sc sc
2 2
Acquisition-Aperture Dip-Response
(Acquisition Configuration Response)
*
S3
*
S2
*
S1
image by common-shot prestack G-D migration Acquisition-Dip-Response (horizontal reflector) from all the 325 shots Acquisition-Dip-Response (45
down from horizontal) from all the 325 shots
G-D beamlet prestack migration image Acquisition-Dip Response for 45 o reflectors Improved image after Directional illumination correction
Conclusion
• • Local image matrix can be obtained from the local incident and scattered plane waves based on beamlet decomposition The goal of true-reflection imaging in angle-domain is to remove the acquisition aperture effect and propagation effect through directional illumination analysis and the corresponding corrections
Conclusion
(
continued
) • • • CDAI and CRAI gathers can be deduced from local image matrices (after corrections) CDAI gathers can be used to determine the dip-angles of local reflectors CRAI gathers can be used for local AVA analysis (and further for local inversion)
Acknowlegement
We thank the support from WTOPI Research Consortium at UCSC We thank the support from DOE Project at UCSC ___________________________________________ Welcome to visit our Consortium booth #2745