Abstract:Two-way wave equation based reverse time migration is powerful in imaging complex area. However, its result usually contains imaging noise which poses challenge to seismic interpretation. The imaging amplitude would be distorted by simple noise filter. Wavefield separation imaging condition is an effective way to suppress the imaging noise and enhance the imaging quality. Conventionally, the wavefield separation is achieved in frequency wavenumber domain which cannot be conveniently implemented in time domain reverse time migration. We propose the analytical wavefield extrapolation method which can efficiently separate wavefield into its directional components during the time domain extrapolation. We use this method in reverse time migration and separate both the source and receiver wavefield. And the imaging noise and signal can be separated after applying the imaging condition to the separated wavefield. The wavefield propagation direction is usually defined by temporal and spatial Fourier transform. In the Fourier domain, the wavefield propagation direction is defined by the sign of frequency and spatial wavenumber. In order to efficiently apply the wavefield separation imaging condition in reverse time migration, we extend the analytical time signal into the wavefield and call it analytical time wavefield. Since the analytical time wavefield contains only the positive frequency component, the wavefield propagation direction can be defined by the sign of spatial wavenumber. To avoid the I/O cost in generating the analytical wavefield at every imaging point, we propose an analytical wavefield propagation equation based on the linear relation between the source term and wavefield in wave equation. We solve the proposed equation by finite difference method. Then we separate the source and receiver wavefield into their up and down going components and apply the imaging condition to the separated wavefield. Four imaging components (i. e., imaging component from up going source and down going receiver wavefield, imaging component from down going source and down going receiver wavefield, imaging component from down going source and up going receiver wavefield and imaging component from up going source and up going receiver wavefield) are effectively separated. Then the imaging noise and signal are separated. We use the syncline model and real data to test the proposed method. Numerical example on syncline model shows that the proposed method can effectively separate imaging component from down going source and up going receiver wavefield, imaging component from up going source and down going receiver wavefield, imaging component from down going source and down going receiver wavefield and imaging component from up going source and up going receiver wavefield. When the migration velocity is not accurate, the correlation imaging condition and equivalent wavefield separation imaging condition would generate imaging artifact in central area of the syncline. And the proposed method can eliminate the artifact. Numerical example on real data shows that this method can generate imaging result free from low frequency noise. Wavefield separation imaging condition can effectively separate the imaging noise and signal in reverse time migration and enhance the imaging quality. The proposed analytical wavefield extrapolation method can Separate the source and receiver wavefield into their up and down going components. The imaging noise can be suppressed by applying the imaging condition to the separated wavefield. We want to further use this method to generate common angle imaging gather in reverse time migration and provide input data for AVA and migration velocity inversion.

Reverse Time Migration Based on Normalized Wavefield Decomposition Imaging Condition

An Effective Imaging Condition for Reverse-Time Migration Using Wavefield Decomposition

An efficient least-squares reverse time migration in image domain

Reflected wave least squares reverse time migration with angle illumination compensation

Reverse Time Migration of Crosswell GPR Data Based on Wavefield Decomposition

Angle-Weighted Reverse Time Migration with Wavefield Decomposition Based on the Optical Flow Vector

Reverse time migration using analytical time wavefield extrapolation and separation

Reverse-time migration using one-way wavefield imaging condition

Reverse Time Migration: A Prospect of Seismic Imaging Methodology

Passive multiple reverse time migration imaging based on wave decomposition and normalized imaging conditions

A causal imaging condition for reverse time migration using the Discrete Hilbert transform and its efficient implementation on GPU

Reverse Time Migration for Extended Obstacles: Electromagnetic Waves

Reverse Time Migration for Extended Obstacles: Acoustic Waves

Mesh-free Radial-Basis-function-generated Finite Differences and Their Application in Reverse Time Migration

Reverse Time Migration Based on the Pseudo-Space-Domain First-Order Velocity-Stress Acoustic Wave Equation

Wavefield decomposition in arbitrary direction and an imaging condition based on stratigraphic dip

Joint Least-Squares Reverse Time Migration of Primary and Prismatic Waves

Waveform Imaging Based on Linear Forward Representations for Scalar Wave Seismic Data

Reverse Time Migration with Source Wavefield Reconstruction Strategy

2-D and 3-D Image-Domain Least-Squares Reverse Time Migration Through Point Spread Functions and Excitation-Amplitude Imaging Condition.

GAN-enhanced directional seismic wavefield decomposition and its application in reverse-time migration