Wave-equation reflection tomography: annihilators and sensitivity kernels

In seismic tomography, the finite frequency content of broad-band data leads to interference effects in the process of medium reconstruction, which are ignored in traditional ray theoretical implementations. Various ways of looking at these effects in the framework of transmission tomography can be found in the literature. Here,we consider inverse scattering of bodywaves to develop a method of wave-equation reflection tomography with broad-band waveform data— which in exploration seismics is identified as a method of wave-equation migration velocity analysis. In the transition from transmission to reflection tomography the usual cross correlation between modelled and observedwaveforms of a particular phase arrival is replaced by the action of operators (annihilators) to the observed broad-bandwavefields. Using the generalized screen expansion for one-way wave propagation, we develop the Fr´echet (or sensitivity) kernel, and show how it can be evaluated with an adjoint state method. We cast the reflection tomography into an optimization procedure; the kernel appears in the gradient of this procedure.We include a numerical example of evaluating the kernel in a modified Marmousi model, which illustrates the complex dependency of the kernel on frequency band and, hence, scale. In heterogeneous media the kernels reflect proper wave dynamics and do not reveal a self-similar dependence on frequency: low-frequency wave components sample preferentially the smoother parts of the model, whereas the high-frequency data are—as expected—more sensitive to the stronger heterogeneity.We develop the concept for acoustic waves but there are no inherent limitations for the extension to the fully elastic case.

[1]  A one-way wave equation for modelling seismic waveform variations due to elastic anisotropy , 2004 .

[2]  W. A. Mulder,et al.  Automatic velocity analysis by differential semblance optimization , 2002 .

[3]  M. V. de Hoop,et al.  Generalized‐Screen Approximation and Algorithm for the Scattering of Elastic Waves , 2003 .

[4]  Ru-Shan Wu,et al.  Generalization of the phase-screen approximation for the scattering of acoustic waves , 2000 .

[5]  Christiaan C. Stolk,et al.  Kinematics of shot-geophone migration , 2009 .

[6]  Maarten V. de Hoop,et al.  Generalization of the Bremmer coupling series , 1996 .

[7]  Adrianus T. de Hoop,et al.  Wave–field reciprocity and optimization in remote sensing , 2000 .

[8]  William W. Symes,et al.  Velocity inversion by differential semblance optimization , 1991 .

[9]  Guust Nolet,et al.  Fréchet kernels for finite-frequency traveltimes—I. Theory , 2000 .

[10]  M. V. Hoop,et al.  Approximate dispersion relations for qP-qSV-waves in transversely isotropic media , 2000 .

[11]  H. Chauris,et al.  Two-dimensional velocity macro model estimation from seismic reflection data by local differential semblance optimization: applications to synthetic and real data sets , 2001 .

[12]  Louis Fishman,et al.  Exact constructions of square-root Helmholtz operator symbols: The focusing quadratic profile , 2000 .

[13]  Maarten V. de Hoop,et al.  Uniform Asymptotic Expansion of the Square-Root Helmholtz Operator and the One-Way Wave Propagator , 2003, SIAM J. Appl. Math..

[14]  Lateral heterogeneity of the upper mantle determined from the travel times of multiple ScS , 1975 .

[15]  John E. Peterson,et al.  Beyond ray tomography: Wavepaths and Fresnel volumes , 1995 .

[16]  B. Ursin,et al.  Seismic velocity analysis in the scattering‐angle/azimuth domain , 2003 .

[17]  Paul Sava,et al.  Amplitude-preserved common image gathers by wave-equation migration , 2001 .

[18]  Van,et al.  Generalized Bremmer series with rational approximation for the scattering of waves in inhomogeneous media , 1996 .

[19]  L. Fishman.  Exact and operator rational approximate solutions of the Helmholtz, Weyl composition equation in underwater acoustics−The quadratic profile , 1992 .

[20]  Jon F. Claerbout,et al.  Imaging the Earth's Interior , 1985 .

[21]  Louis Fishman,et al.  Uniform high-frequency approximations of the square root Helmholtz operator symbol , 1997 .

[22]  R. Gerhard Pratt,et al.  Efficient waveform inversion and imaging: A strategy for selecting temporal frequencies , 2004 .

[23]  Christiaan C. Stolk,et al.  Smooth objective functionals for seismic velocity inversion , 2003 .

[25]  Kim B. Olsen,et al.  Frechet Kernels for Imaging Regional Earth Structure Based on Three-Dimensional Reference Models , 2005 .

[26]  Seismic inverse scattering in the `wave-equation' approach , 2001, math/0112172.

[27]  Jérôme Le Rousseau,et al.  Modeling and imaging with the scalar generalized‐screen algorithms in isotropic media , 2001 .

[28]  Louis Fishman,et al.  Derivation and application of extended parabolic wave theories. II. Path integral representations , 1984 .

[29]  Gerard T. Schuster,et al.  Wave-equation traveltime inversion , 1991 .

[30]  Sergey Fomel,et al.  Theory of differential offset continuation , 2003 .

[31]  Louis Fishman,et al.  Derivation and application of extended parabolic wave theories. I. The factorized Helmholtz equation , 1984 .

[32]  A. J. Berkhout,et al.  Angle-dependent reflectivity by means of prestack migration , 1990 .

[33]  C. Wunsch.  The Ocean Circulation Inverse Problem , 1996 .

[34]  Paul Sava,et al.  Wave‐equation migration velocity analysis , 2004 .

[35]  Jan Thorbecke,et al.  Common focus point technology , 1997 .

[36]  Gunther Uhlmann,et al.  Oscillatory integrals with singular symbols , 1981 .

[37]  T. Jordan,et al.  Mantle layering from ScS reverberations: 1. Waveform inversion of zeroth-order reverberations , 1991 .

[38]  T. Jordan,et al.  Lateral heterogeneity of the upper mantle determined from the travel times of multiple ScS , 1975 .

[39]  Maarten V. de Hoop,et al.  SEISMIC INVERSE SCATTERING IN THE DOWNWARD CONTINUATION APPROACH , 2006 .

[40]  Maarten V. de Hoop,et al.  Modeling of Seismic Data in the Downward Continuation Approach , 2005, SIAM J. Appl. Math..

[41]  Christiaan C. Stolk.  A pseudodifferential equation with damping for one-way wave propagation in inhomogeneous acoustic media , 2003 .

[42]  Qinya Liu,et al.  Tomography, Adjoint Methods, Time-Reversal, and Banana-Doughnut Kernels , 2004 .

[43]  Assessment of crustal velocity models using seismic refraction and reflection tomography , 2003 .

[44]  R. Pratt.  Seismic waveform inversion in the frequency domain; Part 1, Theory and verification in a physical scale model , 1999 .

[45]  Tariq Alkhalifah,et al.  Acoustic approximations for processing in transversely isotropic media , 1998 .

[46]  B. Biondi,et al.  Symplectic structure of wave-equation imaging: a path-integral approach based on the double-square r , 2003 .

[47]  Maarten V. de Hoop,et al.  On sensitivity kernels for ‘wave-equation’ transmission tomography , 2005 .

[48]  M. D. Hoop,et al.  Characterization and `Source-Receiver' Continuation of Seismic Reflection Data , 2006 .

[49]  Paul Sava,et al.  Angle-domain common-image gathers by wavefield continuation methods , 2003 .

[50]  A. J. Berkhout.  Pushing the limits of seismic imaging, Part I: Prestack migration in terms of double dynamic focusing , 1997 .

[51]  S. Wales,et al.  Factorization and path integration of the Helmholtz equation: Numerical algorithms , 1987 .

[52]  F. Scherbaum,et al.  Double beam imaging: Mapping lower mantle heterogeneities using combinations of source and receiver arrays , 1997 .

[53]  Thomas H. Jordan,et al.  Sensitivity of frequency-dependent traveltimes to laterally heterogeneous, anisotropic Earth structure , 1998 .

[54]  Maarten V. de Hoop,et al.  Uniform Asymptotic Expansion of the Generalized Bremmer Series , 2000, SIAM J. Appl. Math..

[55]  A one‐way wave equation for modelling variations in seismic waveforms due to elastic anisotropy , 2004 .

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

[57]  R. Hilst,et al.  Global wave-equation reflection tomography , 2003 .

[58]  Maarten V. de Hoop,et al.  Microlocal analysis of seismic inverse scattering in anisotropic elastic media , 2002 .

[59]  Louis Fishman,et al.  Derivation And Application Of Extended Parabolic Wave Theories , 1982, Optics & Photonics.