Common-Reflection-Surface Stack for Converted Waves

The finite-offset (FO) common-reflection-surface (CRS) stack has been shown to be able to handle not only P-P or S-S but also arbitrarily converted reflections. It can provide different stack sections such as common-offset (CO), common-midpoint (CMP) and common-shot (CS) sections with significantly increased signal-to-noise ratio from the multi-coverage pre-stack seismic data in a data-driven way. It is our purpose in this paper to demonstrate the performance of the FO CRS stack on data involving converted waves in inhomogeneous layered media. In order to do this we apply the FO CRS stack for common-offset to a synthetic seismic data set involving P-P as well as P-S converted primary reflections. We show that the FO CRS stack yields convincing improvement of the image quality in the presence of noisy data and successfully extracts kinematic wavefield attributes useful for further analyses. The extracted emergence angle information is used to achieve a complete separation of the wavefield into its P-P and P-S wave components, given the FO CRS stacked horizontal and vertical component sections.     

Model-independent Travel-time Attributes for 2-D, Finite-offset Multicoverage Reflections

— In this paper, we provide a 5-parameter stacking formula to transform 2-D prestack data into a particular common-offset section. This requires the knowledge of the near-surface velocity only and it is expected that ray theory holds to describe primary reflections. The earth model can be arbitrarily inhomogeneous. The new stacking approach can be viewed as a generalization of the 3-parameter common-reflection-surface (CRS) stack, by which 2-D multicoverage data are stacked into a simulated zero-offset section. The new 5-parameter formula can handle P-P, P-S and S-S reflections.     

Constant velocity migration in the various guises of plane-wave theory

The various analytic schemes for performing a wavefield extrapolation or seismic migration from measurements upon a planar surface within a constant velocity medium are inherently related to each other. All schemes can be derived from a simple plane-wave representation of the recorded wavefield. One scheme that is very easy to conceive is based on the Radon transform. It enables one to perform a wavefield extrapolation or seismic migration by a filtered projection and a back projection of the recorded wavefield. This reveals that the theory of seismic migration as well as the theory of seismic tomography are very closely related to each other.     

A Wave-Front Curvature Approach to Computing Ray Amplitudes in Inhomogeneous Media with Curved Interfaces

A system of three ordinary non-linear first order differential equations is proposed for the computation of the geometrical spreading of the wave front of a seismic body wave in a three-dimensional medium. The variables of the system are the parameters which provide a second order approximation of the wave front.     

True-amplitude migration of 2D synthetic data1

True-amplitude (TA) migration, which is a Kirchhoff-type modified weighted diffraction stack, recovers (possibly) complex angle-dependent reflection coefficients which are important for amplitude-versus-offset (AVO) inversion. The method can be implemented using existing prestack or post-stack Kirchhoff migration and fast Green's function computation programs. Here, it is applied to synthetic single-shot and constant-offset seismic data that include post-critical reflections (complex reflection coefficients) and caustics. Comparisons of the amplitudes of the TA migration image with theoretical reflection coefficients show that the (possibly complex) angle-dependent reflection coefficients are correctly estimated.     

On the computation of shear waves by the ray method

Summary. The ray series solution of the elastodynamic equation of motion for shear waves propagating through a laterally inhomogeneous three-dimensional medium can be simplified by the use of a particular coordinate system that accompanies the wave front along the ray of investigation. The system is entirely determined by parameters that are obtainable from the ray. The transport equations for the principal shear wave components are then no longer coupled, but reduce to the same type of equation which determines the principal compressional wave component.