Imaging Seismic Reflections

T.J.P.M. op 't Root, Timotheus Johannes Petrus Maria Op 't Root

    Research output: ThesisPhD Thesis - Research UT, graduation UT

    126 Downloads (Pure)


    The goal of reflection seismic imaging is making images of the Earth subsurface using surface measurements of reflected seismic waves. Besides the position and orientation of subsurface reflecting interfaces it is a challenge to recover the size or amplitude of the discontinuities. We investigate two methods or techniques of reflection seismic imaging in order to improve the amplitude. First we study the \textsl{one-way wave equation} with attention to the amplitude of the waves. In our study of \textsl{reverse-time migration} the amplitude refers to the amplitude of the image. Though the approach in the thesis is formal, our results have practical implications for seismic imaging algorithms, which improve or correct the amplitudes. The one-way wave equation is a $1^{\rm st}$-order equation that describes wave propagation in a predetermined direction. We derive the equation and identify the steps that determine the wave amplitude. We introduce a \textsl{symmetric square root operator} and a wave field normalization operator and show that they provide the correct amplitude. The idea to use a symmetric square root is generally applicable. Our amplitude claims are numerically verified. The one-way wave equation is an application of \textsl{pseudo-differential operators}. \textsl{Reverse-time migration} (RTM) is an imaging method that uses simulations of the source and receiver wave fields through a slowly varying estimate of the subsurface medium. The receiver wave is an in reverse time continued field that matches the measurements. An \textsl{imaging condition} transforms the fields into an image of the small scale medium contrast. We investigate the linearized inverse problem and use the RTM procedure to reconstruct the perturbation of the medium. We model the scattering event by the scattering operator, which maps the medium perturbation to the scattered wave. We propose an approximate inverse of the scattering operator, derive a novel imaging condition and show that it yields a reconstruction of the perturbation with correct amplitude. The study extensively uses the theory of \textsl{Fourier integral operators} (FIO). Besides that we globally characterize the scattering operator as a FIO, we also obtain a local expression to calculate the amplitude explicitly. The claims are confirmed and illustrated by numerical simulations.
    Original languageUndefined
    Awarding Institution
    • University of Twente
    • van Groesen, Embrecht W.C., Supervisor
    • Stolk, C.C., Advisor
    Thesis sponsors
    Award date6 Apr 2011
    Place of PublicationZutphen
    Print ISBNs978-90-365-3150-4
    Publication statusPublished - 6 Apr 2011


    • MSC-00A69
    • Reverse-time migration
    • One-way wave modeling
    • EWI-19953
    • Seismic inversion
    • IR-76619
    • METIS-277600
    • NWO-VIDI 639.032.509

    Cite this