Abstract:ABSTRACTWe present a Gaussian packet migration method based on Gabor frame decomposition and asymptotic propagation of Gaussian packets. A Gaussian packet has both Gaussian‐shaped time–frequency localization and space–direction localization. Its evolution can be obtained by ray tracing and dynamic ray tracing. In this paper, we first briefly review the concept of Gaussian packets. After discussing how initial parameters affect the shape of a Gaussian packet, we then propose two Gabor‐frame‐based Gaussian packet decomposition methods that can sparsely and accurately represent seismic data. One method is the dreamlet–Gaussian packet method. Dreamlets are physical wavelets defined on an observation plane and can represent seismic data efficiently in the local time–frequency space–wavenumber domain. After decomposition, dreamlet coefficients can be easily converted to the corresponding Gaussian packet coefficients. The other method is the Gabor‐frame Gaussian beam method. In this method, a local slant stack, which is widely used in Gaussian beam migration, is combined with the Gabor frame decomposition to obtain uniform sampled horizontal slowness for each local frequency. Based on these decomposition methods, we derive a poststack depth migration method through the summation of the backpropagated Gaussian packets and the application of the imaging condition. To demonstrate the Gaussian packet evolution and migration/imaging in complex models, we show several numerical examples. We first use the evolution of a single Gaussian packet in media with different complexities to show the accuracy of Gaussian packet propagation. Then we test the point source responses in smoothed varying velocity models to show the accuracy of Gaussian packet summation. Finally, using poststack synthetic data sets of a four‐layer model and the two‐dimensional SEG/EAGE model, we demonstrate the validity and accuracy of the migration method. Compared with the more accurate but more time‐consuming one‐way wave‐equation‐based migration, such as beamlet migration, the Gaussian packet method proposed in this paper can correctly image the major structures of the complex model, especially in subsalt areas, with much higher efficiency. This shows the application potential of Gaussian packet migration in complicated areas.

Seismic modeling using the frozen Gaussian approximation

Frozen Gaussian approximation for the fractional Schrödinger equation

Frozen Gaussian Sampling for Scalar Wave Equations

Gaussian beam reconstruction of seismic data

Convergence of frozen Gaussian approximation for high frequency wave propagation

A new wave field modeling method by using Gaussian Packets

Frozen Gaussian approximation for high frequency wave propagation in periodic media

Borehole and crosshole seismic wave field forward modeling

Modeling seismic wave propagation within complex structures

A fast space-time-domain Gaussian beam migration approach using the dominant frequency approximation

Gabor‐frame‐based Gaussian Packet Migration

Seismic Characterization of Fractured Reservoirs by Focusing Gaussian Beams

Modeling the Propagation of Diffusive-Viscous Waves Using Flux Corrected Transport-Finite Difference Method

Efficient Gaussian packets representation and seismic imaging

Frozen Gaussian Series Representation of the Imaginary Time Propagator Theory and Numerical Tests

Complex topography and media implementation using hp-adaptive DG-FEM for seismic wave modeling

Seismic Data Interpolation Using Gaussian Process Regression

Double-beam Stacking to Infer Seismic Properties of Fractured Reservoirs

Modeling obliquely incident seismic wave propagation in rocks via the FDM-DEM coupling scheme: Algorithms and applications

An Efficient Immersed Free Surface Boundary Method for 3-D Scalar Seismic Waves Finite-Difference Modeling in Presence of Topography

AN EFFICIENT METHOD TO MODEL SEISMIC PROPAGATION IN DIFFUSIVE-VISCOUS MEDIA WITH DIPPING INTERFACES