Abstract:The sparse expression of seismic data is important to signal processing, migration and parameter inversion. The conventional linear Radon transform is generalized to invert for both the local plane-wave and shot/receiver ray parameters simultaneously, which can compress the prestack seismic data in the local plane-wave domain. Consequently, it will be convenient to the subsequent ray-based imaging and inversion. Based on the combination of shot/receiver ray parameters, a new plane-wave prediction model is presented which can predict the local plane-wave within the spatial window of local shot/receiver center. The basic assumption is the local linear property of the wave-front. The seismic traces within the spatial window of local shot/receiver center are time-shifted and stacked according to the shot/receiver ray parameters and spatial distance, and the seismic phase with the given shot/receiver ray parameter can be separated naturally. Since the presented time-shift and stack procedure relies on the propagation direction of local plane-wave, we call it characteristic wave decomposition (CWD). Besides, the inversion-based CWD under the framework of compressed sensing (CS) theory can invert for both the local plane-wave and shot/receiver ray parameters simultaneously. Due to the pre-estimated ray-parameters for shot/receiver center, the loop for shot/receiver ray-parameters is avoided. Therefore the efficiency of conventional ray-based migration can be increased remarkably, and the ellipsoid migration response in constant media can be simplified to a point response. The new imaging scheme is called the characteristic wave imaging (CWI). First, the proposed CWD method is tested with the Marmousi synthetic data. Numerical examples show that the seismic phases can be separated completely. Meanwhile, even if the original seismic traces are contaminated by severe random noise, the numerical results demonstrate the robustness of CWD method. Then, the CWI method is tested with the Waxian synthetic data. Due to the high S/N ratio, the seismic data are compressed considerably. Compared with Kirchhoff PSDM method, the imaging efficiency is increased over one order of magnitude. Due to the pre-estimated ray-parameters, the angle domain CIGs can be produced directly. Finally, the CWD method is applied to the separation of first-arrival waveform, serving as the input data for wave-equation based transmission traveltime tomography. The proposed CWD method can compress the prestack seismic data in the local plane-wave domain. Under the framework of compressed sensing theory, both the local plane-wave and shot/receiver ray parameters can be inverted simultaneously. However, the local linear property of wave-front is the basic assumption of CWD. Therefore, it is not suitable for complex topography. Besides, in case of sparse shot/receiver sampling, it is difficult to invert for the shot/receiver ray-parameter simultaneously. The double beam forming procedure in CWD will be degenerated to the classical beam forming method. The presented CWI scheme can map the local plane-wave into the model space with the pre-estimated shot/receiver ray-parameters, so that the imaging efficiency is increased and the migration noise is attenuated. Meanwhile, the angle domain CIGs can be produced easily. Combined with the angle-domain tomography theory, it may be beneficial to the automation of velocity model building procedure.

Technology and application of beam ray prestack imaging method in foothill areas

A Method of Plane Wave Depth Migration for Pre-Stack Seismic Data

Prestack Depth Migration Method Based on Quasi-Linear Born Approximation

Generalized Synthesis of Seismic Data and Its Migration

Reflected wave least squares reverse time migration with angle illumination compensation

Steeply-dipping structures migration based on the wavefield vertical-and horizontal-extrapolation

High-Resolution Anisotropic Prestack Kirchhoff Dynamic Focused Beam Migration

Wave Theory Based Seismic Survey Design

Wave Equation Migration for 3-D VSP Data

Beam-ray gather stacking for attenuating noncoherent noises and PSTM from rugged topography

An efficient 3D wave-equation pre-stack time migration for high-density and wide-azimuth data

Study on the absorption tomography with pre-stack seismic reflection data based on ray theory

Characteristic wavefield decomposition, imaging and inversion with prestack seismic data

Finite-difference pre-stack time migration of 2D offset plane wave

Local Angle Domain Target Oriented Illumination Analysis and Imaging Using Beamlets

VSP COMMON SHOT PRE-STACK DEPTH MIGRATION WITH HYBRID METHOD

Efficient pre-stack depth migration method using characteristic Gaussian packets

Common offset 3D prestack seismic data anti-aliasing beamforming, migration and imaging

Target-oriented beamlet migration based on Gabor-Daubechies frame decomposition

A 3D plane-wave beamforming method for borehole azimuthal acoustic reflection imaging

A New Stacking Technique in Ambient Noise Tomography