Abstract:The traditional ray method don't perform well in the caustic and shadow region because its image accuracy is low. While the 2D Gaussian beam prestack depth migration in the frequency domain is less efficient because it needs calculation of Green function of all frequencies at imaging points. To solve these problems, we propose a novel method of 2D Gaussian beam migration of common-shot records in the time domain. In Gaussian beam migration, the source and recorded wavefield are both based on the Green functions staked by Gaussian beam from the source point and the receive point. Thus we can use the source wavefield and the recorded wavefield to generate the migrated image. Amplitude and exponential terms are both complex numbers, as Green function is complex and the traveltime expressed by the exponential term is called complex traveltime. This complex traveltime at the image point is determined by the second-order Taylor expansion of the nearest point's traveltime in the center ray. With the image condition, the migration formula is a triple integral of the horizontal slowness, vertical slowness and frequency, and the integrand is also a complex number. Focused on the real part of the integrand, the triple integral is reduced to a double integral by transforming the integral variable time to frequency using the Fourier transform. In the triple integral with the real part of the integrand, we change the sequence, and transform the frequency to the time by the Fourier transform to get the result that is two-fold integral of the horizontal slowness and the vertical slowness at the image time. The image time in the formula is the summation of the complex traveltime from the source point to the image point and that from the receiver point to the image point. Because the amplitude is a complex number, so the result with the imaginary of the amplitude is needed to transform the product of the Green function and the seismic data with the Hilbert transform. Then the triple integral becomes a double integral of the horizon slowness and the vertical slowness, and the image formula becomes the summation of two double integrals. The validity of the new method is proved on a complex numerical model. Using the Marmousi data set, the imaging result and computational efficiency of Gaussian beam migration in the time domain is compared with the method in the frequency domain. We can find that the image results of these two methods are similar, while the migration algorithm in the time domain takes less computing time. Then we compare the imaging result of Gaussian beam migration with that of Kirchhoff migration on the Sigsbee model. The migration result of Gaussian beam migration is clearer, especially for the region under the high-velocity geologic body. Gaussian beam migration section can show clear interfaces and even faults under this body which appears as migration illusion in the Kirchhoff migration. The result compared with the migrated imaging result based on the ray method indicates that Gaussian beam migration in the time domain has a higher accuracy. In this article, we present a novel method of 2D common-shot Gaussian beam migration in the time domain by employing the Fourier transform and the Hilbert transform. The ability of this method to deal with complicated models has been validated using the Marmousi data set and Sigsbee data set. Furthermore, because of holding the similar integral kernels, this method can be applied to other types of Gaussian beam migration straightforwardly, such as common-offset, common-receiver, common-midpoint and even 3D conditions.

Azimuth–opening Angle Domain Imaging in 3D Gaussian Beam Depth Migration

Common-azimuth-opening angle imaging gathers from 3D Gaussian Beam migration

A computational method for wide-azimuth 3D dip-angle gathers using Gaussian beam migration

Z-99 MIGRATION AMPLITUDE CORRECTIONS IN ANGLE-DOMAIN USING BEAMLET DECOMPOSITION

Angle Domain Common Image Gather for Wave Equation

Wave Equation Migration in Angle Domain

Local angle-domain Gaussian beam prestack depth migration in a TI medium

Acquisition Aperture Correction in Angle‐domain and True‐amplitude Imaging for Wave Equation Migration

Angle-domain common-image gathers from plane-wave least-squares reverse time migration

Angle Gathers from Reverse Time Migration Using Analytic Wavefield Propagation and Decomposition in the Time Domain

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

An Effective Acoustic Impedance Imaging Based on a Broadband Gaussian Beam Migration

Angle-domain Common Imaging Gather Extraction Via Kirchhoff Prestack Depth Migration Based on a Traveltime Table in Transversely Isotropic Media

2-D Gaussian-beam migration of common-shot records in time domain

Azimuth-preserved Local Angle Domain Imaging and Its Application to Prestack Time Migration

Azimuth-angle domain Kirchhoff pre-stack time migration

3-D Wave Equation Common Azimuth Prestack Migration

The Azimuth Angle Domain Common Image Gathers Extraction by Kirchhoff Migration

Cauchy Wave Field and Its Application in Angle-Domain Imaging Condition

A True-Amplitude Imaging Method Based on Gaussian Beam Migration and Demigration

Three-dimensional Angle-Domain Double-Square-root Migration in VTI Media for the Large-Scale Wide-Azimuth Seismic Data