SUMMARY
Equations are derived for the point-scatterer response of the crosswell migration and traveltime tomography operators. These formulas are used to estimate the limits of spatial resolution in reflection migration images and traveltime tomograms. In particular, for a crosswell geometry with borehole length L, well offset 2x0, source wavelength δ, and a centred point scatterer I show the following.
(1) The vertical resolution Δzmig of the reflection migration image is equal to 2δx0/L under the far-field approximation. Under the Fresnel approximation, Δzmig≈2δx0/L-2δ2/L, which says that images become better resolved with an increase in aperture and a decrease in wavelength and well offset.
(2) The horizontal resolution Δzmig of the migration image is equal to 16δx02/L2. The lateral resolution of the migrated image is worse than the vertical resolution by a factor of 8x0/L (where x0/L> 1 under the far-field approximation).
(3) The vertical resolution Δztomo of the traveltime velocity tomogram is proportional to √λx0. This estimate agrees with that of a previous study. However, the tomographic image of the slowness perturbation behaves as a non-local cosine function along the depth axis, whereas the migration image behaves as a localized squared sine function in the depth coordinate. This is consistent with the empirical observation that interface boundaries are more sharply resolved by migration than by traveltime tomography.
(4) The horizontal resolution of the slowness image in a traveltime tomogram is equal to (4x0/L}√3x0δ/4, a factor more than 4x0/L worse than the vertical resolution.
(5) For Ns sources and Ng geophones, the dynamic range of the migrated image is proportional to NsNg. The dynamic range of the slowness tomogram is proportional to √NgNs.
Many of the estimates for the resolution limits have simple geometrical interpretations. For example, the minimum vertical (horizontal) resolution in a migrated reflectivity section corresponds to the minimum vertical (horizontal) stretch that a migrated wavelet undergoes in going from the time domain to the depth domain. In addition, the minimum vertical (horizontal) resolution in a traveltime tomogram corresponds to the minimum vertical (horizontal) width of the wavepath at the scatterer location
[1]
Gerard T. Schuster,et al.
Wavepath eikonal traveltime inversion: Theory
,
1993
.
[2]
Robert H. Stolt,et al.
Seismic migration : theory and practice
,
1986
.
[3]
S Hanish,et al.
Acoustic Fluctuation Workshop, Held at Naval Research Laboratory, Washington, DC on February 22-23 1978. Technical Review, Editorial Summary, Synopsis and Papers. Volume 1.
,
1978
.
[4]
M. Toksöz,et al.
Diffraction tomography and multisource holography applied to seismic imaging
,
1987
.
[5]
Bernard D. Steinberg,et al.
Microwave Imaging Techniques
,
1991
.
[6]
M. Worthington,et al.
Resolution limits in ray tomography due to wave behavior: Numerical experiments
,
1993
.
[7]
G. Schuster,et al.
Acoustic wave-equation traveltime and waveform inversion of crosshole seismic data
,
1995
.
[8]
Marta Woodward,et al.
Wave-equation tomography
,
1992
.
[9]
M. Soumekh.
Fourier Array Imaging
,
1994
.
[10]
P. Williamson,et al.
A guide to the limits of resolution imposed by scattering in ray tomography
,
1991
.
[11]
I. S. Gradshteyn,et al.
Table of Integrals, Series, and Products
,
1976
.
[12]
W. C. Elmore,et al.
Physics of Waves
,
1969
.
[13]
S. Hanish.
A treatise on acoustic radiation
,
1989
.
[14]
J. Rector,et al.
Characterization of resolution and uniqueness in crosswell direct‐arrival traveltime tomography using the Fourier projection slice theorem
,
1994
.