Abstract

SUMMARY We carry out high-frequency analyses of Claerbout's double-square-root equation and its (numerical) solution procedures in heterogeneous media. We show that the double-square-root equation generates the adjoint of the single-scattering modelling operator upon substituting the leading term of the generalized Bremmer series for the background Green function. This adjoint operator yields the process of ‘wave-equation’ imaging. We finally decompose the wave-equation imaging process into common image point gathers in accordance with the characteristic strips in the wavefront set of the data.

Highlights

  • Directional wavefield decomposition is a tool for analysing and computing the propagation of waves in configurations with a certain global directionality, such as waveguiding structures

  • The method consists of three main steps: (1) decomposing the field into two constituents and propagating them upward or downward along a principal or preferred direction, (2) computing the interaction of the counterpropagating constituents and (3) recomposing the constituents into observables at the positions of interest. (The preferred direction is ‘vertical’, whereas the directions orthogonal to the preferred direction are referred to as the ‘lateral’ directions.) The method allows one to ‘trace’ wave constituents in a global coordinate system in a given medium, and distinguish constituents that have been scattered along the up/down direction a different number of times

  • Computational aspects of the Bremmer series can be found in Van Stralen et al (1998) and Le Rousseau & De Hoop (2001a)

Read more

Summary

INTRODUCTION

Directional wavefield decomposition is a tool for analysing and computing the propagation of waves in configurations with a certain global directionality, such as waveguiding structures. If the medium of the configuration were laterally homogeneous, the directional decomposition would become an algebraic operation in the lateral Fourier or wavenumber domain (see, e.g., Kennett 1985) In such a medium, the phase-shift method (Gazdag & Sguazzero 1984) is amongst the fastest (and accurate) one-way algorithms. The measurements of the sources and the receivers are propagated into the subsurface in a presumed background medium back to the time when the scattering took place. This continuation is governed by the double-square-root equation. Throughout we will make use of microlocal techniques, the fundamental concepts of which are introduced in Appendix A

ONE-WAY WAVE EQUATIONS
The reduced system of equations
The coupled system of one-way wave equations
TROTTER-PRODUCT ONE-WAY WAVE PROPAGATOR
The product integral
R 2π dω
The canonical relation
Propagation of singularities
Square-root Hamiltonian system
Amplitude of the propagator kernel
TROTTER-PRODUCT DOUBLE-SQUARE-ROOT PROPAGATOR
The imaging kernel
The isochrone and its cotangent directions
Angle-gather imaging condition
NUMERICAL EXAMPLE
DISCUSSION
The tangent and cotangent bundles
Symplectic and Lagrangian manifolds
Phase functions
Oscillatory integrals and Fourier integral operators
Full Text
Paper version not known

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.