A correlation-based misfit criterion for wave-equation traveltime tomography

Wave-equation traveltime tomography tries to obtain a subsurface velocity model from seismic data, either passive or active, that explains their traveltimes. A key step is the extraction of traveltime differences, or relative phase shifts, between observed and modelled finite-frequency waveforms. A standard approach involves a correlation of the observed and measured waveforms. When the amplitude spectra of the waveforms are identical, the maximum of the correlation is indicative of the relative phase shift. When the amplitude spectra are not identical, however, this argument is no longer valid. We propose an alternative criterion to measure the relative phase shift. This misfit criterion is a weighted norm of the correlation and is less sensitive to differences in the amplitude spectra. For practical application it is important to use a sensitivity kernel that is consistent with the way the misfit is measured. We derive this sensitivity kernel and show how it differs from the standard banana–doughnut sensitivity kernel. We illustrate the approach on a cross-well data set.

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

[2]  R. Plessix A review of the adjoint-state method for computing the gradient of a functional with geophysical applications , 2006 .

[3]  A. Tarantola,et al.  Generalized Nonlinear Inverse Problems Solved Using the Least Squares Criterion (Paper 1R1855) , 1982 .

[4]  Wim A. Mulder,et al.  FAST TRACK PAPER: The exact solution of the time-harmonic wave equation for a linear velocity profile , 2006 .

[5]  Guust Nolet,et al.  Three-dimensional waveform sensitivity kernels , 1998 .

[6]  G. Folland Harmonic Analysis in Phase Space. (AM-122), Volume 122 , 1989 .

[7]  C. Vogel Computational Methods for Inverse Problems , 1987 .

[8]  Wim A. Mulder,et al.  A comparison of seismic velocity inversion methods for layered acoustics , 2009 .

[9]  C. Pekeris,et al.  Theory of Propagation of Sound in a Half‐Space of Variable Sound Velocity under Conditions of Formation of a Shadow Zone , 1946 .

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

[11]  Guust Nolet,et al.  Three-dimensional sensitivity kernels for finite-frequency traveltimes: the banana–doughnut paradox , 1999 .

[12]  M. Cara,et al.  Waveform inversion using secondary observables , 1987 .

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

[14]  W. A. Mulder,et al.  E-15 TIME- VERSUS FREQUENCY-DOMAIN MODELLING OF SEISMIC WAVE PROPAGATION , 2002 .

[15]  Comparing Virtual Versus Real Crosswell Surveys , 2008 .

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

[17]  Detection of Wave Front Set Perturbations via Correlation: Foundation for Wave-equation Tomography , 2001, math-ph/0104003.

[18]  K. Hannabuss,et al.  HARMONIC ANALYSIS IN PHASE SPACE: (Annals of Mathematics Studies 122) , 1990 .

[19]  Thomas H. Jordan,et al.  Generalized seismological data functionals , 1992 .

[20]  Wim A. Mulder,et al.  Velocity analysis based on data correlation , 2008 .

[21]  R. Strichartz A Guide to Distribution Theory and Fourier Transforms , 1994 .

[22]  William W. Symes,et al.  Migration velocity analysis and waveform inversion , 2008 .

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

[24]  G. Folland Harmonic analysis in phase space , 1989 .