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.

SIMULTANEOUS SEPARATION, INTERPOLATION AND TUBE WAVE SUPPRESSION OF VERTICAL SEISMIC PROFILING USING MATCHING PURSUIT BASED SPARSE BEAM FORMING

Separation Method for Multi-Source Blended Seismic Data

Data fusion for resolution improvement by combining seismic data with logging data

High Dimensional Seismic Data Interpolation with Weighted Matching Pursuit Based on Compressed Sensing

Deep Learning-Based P- and S-Wave Separation for Multicomponent Vertical Seismic Profiling.

Wave Field Separation of Vertical Seismic Profiling in Curvelet Domain

Simultaneous Interpolation and Deblending of 3-D Seismic Data by Iterative Thresholding

VSP wave field separation: An optimization method based on block relaxation and singular value thresholding

Structure-adapted Multichannel Matching Pursuit for Seismic Trace Decomposition

Beam forming by Radon spectrum constraint and its application

Three-Dimensional Curvelet-Based WPOCS Reconstruction, Doppler-Shift Correction, and Noise Attenuation of Jittered Sampled Marine Vibrator Data

Synchronous interpolation and denoising in simultaneous-source data separation

Multichannel Seismic Attenuation Compensation and Interpolation with Curvelet Sparse Constraint of Frequency-Wavenumber Spectrum

Automatic Time Picking for Weak Seismic Phase in the Strong Noise and Interference Environment: An Hybrid Method Based on Array Similarity

Simultaneous Seismic Data Interpolation and Denoising Based on Nonsubsampled Contourlet Transform Integrating With Two-Step Iterative Log Thresholding Algorithm

A novel high-resolution imaging method based on sparsity in the time domain and spectral fitting in the frequency domain

VSP Wavefield Separation Based on GCT and Clustering Algorithm Combining with Phase Space Features

Reverse time migration using analytical time wavefield extrapolation and separation

Wavefield separation using irreversible-migration filtering

Simultaneous Reconstruction and Denoising of Extremely Sparse 5-D Seismic Data by a Simple and Effective Method

Separation of magnetotelluric signals based on refined composite multiscale dispersion entropy and orthogonal matching pursuit