Reverse-time migration-based reflection tomography using teleseismic free surface multiples
Name
Van der Hilst - Burdick 2014.pdf
Size
2.07 MB
Format
Adobe PDF
Checksum (MD5)
e79c59402c6669adb8e1c80aca2b0723
Author(s)
Burdick, S.
de Hoop, Maarten
Wang, S.
van der Hilst, Robert D
Date Issued
November 2013
Journal
Geophysical Journal International
Publisher
Oxford University Press on behalf of The Royal Astronomical Society
Citation
Burdick, S., M. V. de Hoop, S. Wang, and R. D. v. d. Hilst. “Reverse-Time Migration-Based Reflection Tomography Using Teleseismic Free Surface Multiples.” Geophysical Journal International 196, no. 2 (November 7, 2013): 996–1017.
Version
Author's final manuscript
Abstract
Converted and multiply reflected phases from teleseismic events are routinely used to create structural images of the crust–mantle boundary (Moho) and the elasticity contrasts within the crust and upper mantle. The accuracy of these images is to a large extent determined by the background velocity model used to propagate these phases to depth. In order to improve estimates of 3-D velocity variations and, hence, improve imaging, we develop a method of reverse-time migration-based reflection tomography for use with wavefields from teleseismic earthquakes recorded at broad-band seismograph arrays. Reflection tomography makes use of data redundancy—that is, the ability to generate numerous structural images of the subsurface with different parts of the wavefield. In exploration seismology (where it is known as migration velocity analysis) reflection tomography typically involves the generation of an extended image (e.g. offset- or angle-gathers), and the fitness of the background model is evaluated through the application of image-domain annihilators. In regional-scale passive source seismology, however, annihilation-based methods are inadequate because the sparse and irregular distribution of teleseismic sources is not likely to produce illumination over a sufficient range of angles. To overcome this problem we turn towards a source-indexed moveout scheme. Instead of extended image annihilation, we determine the success of the tomographic velocity model by cross correlating images produced with multiply scattered waves from different teleseismic sources. The optimal velocity model is the one that minimizes correlation power between windowed images away from zero depth shift. We base our inversion scheme on the seismic adjoint method and a conjugate gradient solver. For each image pair, the update direction is determined by correlations between downgoing wavefields with upgoing adjoint wavefields for both images. The sensitivity kernels used in this method is similar to those found in other forms of adjoint tomography, but their shapes are controlled by the spatial distribution of the error function. We present the method and a proof-of-concept with 2-D synthetic data.
MIT Department
Massachusetts Institute of Technology. Department of Earth, Atmospheric, and Planetary Sciences
Terms of Use
Creative Commons Attribution-Noncommercial-Share Alike
Persistent DSpace Link
DOI of Published Version
https://doi.org/10.1093/gji/ggt428
