A new slant on seismic imaging: Migration and integral geometry

A new approach to seismic migration formalizes the classical diffraction (or common-tangent) stack by relating it to linearized seismic inversion and the generalized Radon transform. This approach recasts migration as the problem of reconstructing the earth’s acoustic scattering potential from its integrals over isochron surfaces. The theory rests on a solution of the wave equation with the geometrical-optics Green function and an approximate inversion formula for the generalized Radon transform. The method can handle both complex velocity models and (nearly) arbitrary configurations of sources and receivers. In this general case, the method can be implemented as a weighted diffraction stack, with the weights determined by tracing rays from image points to the experiment’s sources and receivers. When tested on a finite-difference simulation of a deviated-well vertical seismic profile (a hybrid experiment which is difficult to treat with conventional wave-equation methods), the algorithm accurately reconstructed faulted-earth models. Analytical reconstruction formulas are derived from the general formula for zero-offset and fixed-offset surface experiments in which the background velocity is constant. The zero-offset inversion formula resembles standard Kirchhoff migration. Our analysis provides a direct connection between the experimental setup (source and receiver positions, source wavelet, background velocity) and the spatial resolution of the reconstruction. Synthetic examples illustrate that the lateral resolution in seismic images is described well by the theory and is improved greatly by combining surface data and borehole data. The best resolution is obtained from a zero-offset experiment that surrounds the region to be imaged.

[1]  William A. Schneider,et al.  INTEGRAL FORMULATION FOR MIGRATION IN TWO AND THREE DIMENSIONS , 1978 .

[2]  F. Harris.  On the use of windows for harmonic analysis with the discrete Fourier transform , 1978, Proceedings of the IEEE.

[3]  Norman Bleistein,et al.  AN EXTENSION OF THE BORN INVERSION METHOD TO A DEPTH DEPENDENT REFERENCE PROFILE , 1984 .

[4]  V. Edwards.  Scattering Theory , 1973, Nature.

[5]  I. Gel'fand,et al.  Differential forms and integral geometry , 1969 .

[6]  Jon F. Claerbout,et al.  DOWNWARD CONTINUATION OF MOVEOUT‐CORRECTED SEISMOGRAMS , 1972 .

[7]  Michael Oristaglio,et al.  Inversion Procedure for Inverse Scattering within the Distorted-Wave Born Approximation , 1983 .

[8]  A. Tarantola.  Inversion of seismic reflection data in the acoustic approximation , 1984 .

[9]  J. Hagedoorn,et al.  A process of seismic reflection interpretation , 1954 .

[10]  M. A. Fiddy,et al.  The Radon Transform and Some of Its Applications , 1985 .

[11]  E. Wolf.  Three-dimensional structure determination of semi-transparent objects from holographic data , 1969 .

[12]  Mogens Flensted‐Jensen,et al.  GROUPS AND GEOMETRIC ANALYSIS Integral Geometry, Invariant Differential Operators, and Spherical Functions (Pure and Applied Mathematics: A Series of Monographs and Textbooks) , 1985 .

[13]  G. Shilov,et al.  Generalized Functions, Volume 1: Properties and Operations , 1967 .

[14]  Stephen J. Norton,et al.  Ultrasonic Reflectivity Imaging in Three Dimensions: Exact Inverse Scattering Solutions for Plane, Cylindrical, and Spherical Apertures , 1981, IEEE Transactions on Biomedical Engineering.

[15]  J. D. Johnson,et al.  5. Migration—The Inverse Method , 1982 .

[16]  A. J. Berkhout.  Multidimensional linearized inversion and seismic migration , 1984 .

[17]  S. Helgason.  Groups and geometric analysis , 1984 .

[18]  William S. French,et al.  Computer migration of oblique seismic reflection profiles , 1975 .

[19]  William A. Schneider,et al.  DEVELOPMENTS IN SEISMIC DATA PROCESSING AND ANALYSIS (1968–1970) , 1971 .

[20]  William S. French,et al.  TWO‐DIMENSIONAL AND THREE‐DIMENSIONAL MIGRATION OF MODEL‐EXPERIMENT REFLECTION PROFILES , 1974 .

[21]  Gregory Beylkin,et al.  Spatial Resolution of Migration Algorithms , 1985 .

[22]  R. Stolt.  MIGRATION BY FOURIER TRANSFORM , 1978 .

[23]  Eric Todd Quinto,et al.  The dependence of the generalized Radon transform on defining measures , 1980 .

[24]  P. Schultz.  Velocity estimation by wave front synthesis , 1976 .

[25]  Robert W. Clayton,et al.  A Born-WKBJ inversion method for acoustic reflection data , 1981 .

[26]  G. Beylkin.  The inversion problem and applications of the generalized radon transform , 1984 .

[27]  A. Devaney.  A filtered backpropagation algorithm for diffraction tomography. , 1982, Ultrasonic imaging.

[28]  Gregory Beylkin,et al.  Imaging of discontinuities in the inverse scattering problem by inversion of a causal generalized Radon transform , 1985 .

[29]  Arthur B. Weglein,et al.  Migration and inversion of seismic data , 1985 .

[30]  John A. Fawcett,et al.  Inversion of N-dimensional spherical averages , 1985 .

[31]  J. Claerbout.  Toward a unified theory of reflector mapping , 1971 .

[32]  Gregory Beylkin,et al.  Distorted-wave born and distorted-wave rytov approximations , 1985 .

[33]  G. H. F. Gardner,et al.  Elements of migration and velocity analysis , 1974 .