A later-arrival based inversion scheme to recover diffractors and reflectors

Seismic tomography techniques based on first break time inversion are a convenient way to recover velocity fields. However, these techniques only use a small fraction of the available seismogram data. Migration uses entire reflection seismograms but stops short of an actual inversion, using a summing/averaging method. The result is a smeared image of diffractor points, which makes it difficult to image weaker events. In this paper we develop a seismic inversion technique to recover diffractors and reflectors from full-waveform seismograms by direct inversion, rather than a migration-based approach. We choose a simple, idealised modelling scheme to highlight the basic advantages and limitations inherent in this class of inversion techniques, and derive several important and generic guidelines applicable to the design and development of full-waveform-based inversion algorithms. We propose a modelling scheme for diffractions, reflections and refractions in an arbitrary velocity field. By using a sourcewavelet based estimation method, we linearise the forward model and are able to obtain a diffractor/reflector field recovery scheme that is based on standard iterative matrix inversion. Our investigation reveals that the algorithm provides vastly superior resolution of the structure compared to time-of-flight tomography and requires only a surprisingly small number of sources and receivers for successful operation. However, it is highly prone to noise capture, where noise is any system unknown. This translates into inaccuracies in the macro-velocity field estimate and the seismic wave amplitudes. Our experiments show that the background velocity field must be known with an accuracy of better than 1%, suggesting that this class of linearised inversion scheme is unsuitable for field data sets. Nevertheless, the scheme highlights the pitfalls of any full-waveform inversion approach, where without careful calibration, wavefield dynamics can easily be mis-inverted to produce spurious structure.

[1]  Shunhua Cao,et al.  Finite-difference solution of the eikonal equation using an efficient, first-arrival, wavefront tracking scheme , 1994 .

[2]  Reinhold Ludwig,et al.  Finite-Element Formulation of Acoustic Scattering Phenomena with Absorbing Boundary-Condition in the Frequency-Domain , 1993 .

[3]  Leiv-J. Gelius Limited-view diffraction tomography in a nonuniform background , 1995 .

[4]  David C. Henley,et al.  Full‐waveform processing and interpretation of kilohertz cross‐well seismic data , 1993 .

[5]  G. Schuster,et al.  Acoustic wave-equation traveltime and waveform inversion of crosshole seismic data , 1995 .

[6]  T. A. Dickens Study of Diffraction Tomography For Layered Backgrounds , 1992 .

[7]  P. Williamson,et al.  A guide to the limits of resolution imposed by scattering in ray tomography , 1991 .

[8]  Gerard T. Schuster,et al.  Acoustic Wave Equation Traveltime And Waveform Inversion of Crosshole Seismic Data , 1993 .

[9]  Z. M. Song,et al.  Frequency-domain acoustic-wave modeling and inversion of crosshole data; Part II, Inversion method, synthetic experiments and real-data results , 1995 .

[10]  R. Pratt,et al.  Combining wave-equation imaging with traveltime tomography to form high-resolution images from crosshole data , 1991 .

[11]  R. G. Pratt,et al.  Two‐dimensional velocity models from wide‐angle seismic data by wavefield inversion , 1996 .

[12]  R. Pratt,et al.  The application of diffraction tomography to cross-hole seismic data , 1988 .

[13]  Mehmet Ali C.Tura Application of Diffraction Tomography to Fracture Detection , 1991 .

[14]  Jack K. Cohen,et al.  Two and one‐half dimensional Born inversion with an arbitrary reference , 1987 .

[15]  Jean Brac,et al.  Can we image complex structures with first‐arrival traveltime? , 1993 .

[16]  Steven Carroll,et al.  Techniques to improve imaging of seismic tomographic data , 1997 .

[17]  Zhou Bing,et al.  Iterative algorithm for the damped minimum norm, least-squares and constrained problem in seismic tomography , 1992 .

[18]  Robert L. Coultrip High-accuracy wavefront tracing traveltime calculation , 1993 .

[19]  Robert R. Stewart,et al.  Exploration seismic tomography : fundamentals , 1991 .

[20]  Robert H. Stolt,et al.  Seismic migration : theory and practice , 1986 .

[21]  Thomas Gruber,et al.  Precision analysis of first-break times in grid models , 1998 .

[22]  M. Toksöz,et al.  Diffraction tomography and multisource holography applied to seismic imaging , 1987 .

[23]  Reinhold Ludwig,et al.  Nonlinear Diffractive Inverse Scattering for Multiple-Scattering in Inhomogeneous Acoustic Background Media , 1995 .

[24]  W. Rodi,et al.  Nonlinear waveform tomography applied to crosshole seismic data , 1996 .

[25]  L. Johnson,et al.  Application of diffraction tomography to fracture detection , 1992 .

[26]  Shunhua Cao,et al.  Relative-error-based non-linear inversion: application to seismic traveltime tomography , 1995 .

[27]  Gerard T. Schuster,et al.  Wavepath eikonal traveltime inversion: Theory , 1993 .

[28]  David M. Pai,et al.  Crosshole seismic using vertical eigenstates , 1990 .

[29]  Shunhua Cao,et al.  Calculation of the seismic first-break time field and its ray path distribution using a minimum traveltime tree algorithm , 1993 .

[30]  Peter Mora Inversion = Migration + Tomography , 1988, Shell Conference.

[31]  Marta Woodward,et al.  Wave-equation tomography , 1992 .

[32]  William H. Press,et al.  Numerical recipes , 1990 .

[33]  M. Worthington,et al.  An introduction to geophysical tomography , 1984 .

[34]  T. J. Moser Migration using the shortest-path method , 1994 .

[35]  Dan Loewenthal,et al.  Estimation of a point source wavelet in elastic media , 1993 .