Multiscale reverse-time-migration-type imaging using the dyadic parabolic decomposition of phase space

TL;DR

This paper introduces a novel multiscale reverse-time migration method using dyadic parabolic decomposition of phase space, enabling wave packet-based wave equation solutions for improved seismic imaging.

Contribution

It presents a new Fourier integral operator framework for reverse-time migration with a dyadic parabolic decomposition, enhancing numerical wavefield continuation and imaging.

Findings

Developed a wave packet-based solution for wave equations

Created algorithms for reverse time continuation from scattering data

Demonstrated improved RTM migration techniques

Abstract

We develop a representation of reverse-time migration in terms of Fourier integral operators the canonical relations of which are graphs. Through the dyadic parabolic decomposition of phase space, we obtain the solution of the wave equation with a boundary source and homogeneous initial conditions using wave packets. On this basis, we develop a numerical procedure for the reverse time continuation from the boundary of scattering data and for RTM migration. The algorithms are derived from those we recently developed for the discrete approximate evaluation of the action of Fourier integral operators and inherit from their conceptual and numerical properties.