CONVERSION POINTS AND TRAVELTIMES OF CONVERTED WAVES IN PARALLEL DIPPING LAYERS1

For converted waves, stacking as well as AVO analysis requires a true common reflection point gather which, in this case, is also a common conversion point (CCP) gather. The coordinates of the conversion points for PS or SP waves, in a single homogeneous layer can be calculated exactly as a function of the offset, the reflector depth and the ratio v_p/v_s. An approximation of the conversion point on a dipping interface as well as for a stack of parallel dipping layers is given. Numerical tests show that the approximation can be used for offsets smaller than the depth of the reflector under consideration. The traveltime of converted waves in horizontal layers can be expanded into a power series. For small offsets a two-term truncation of the series yields a good approximation. This approximation can also be used in the case of dipping reflectors if a correction is applied to the traveltimes. This correction can be calculated from the approximated conversion point coordinates.     

NUMERICAL SOLUTION OF THE ACOUSTIC AND ELASTIC WAVE EQUATIONS BY A NEW RAPID EXPANSION METHOD1

We present a new rapid expansion method (REM) for the time integration of the acoustic wave equation and the equations of dynamic elasticity in two spatial dimensions. The method is applicable to spatial grid methods such as finite differences, finite elements or the Fourier method. It is based on a Chebyshev expansion of the formal solution to the appropriate wave equation written in operator form. The method yields machine accuracy yet it is faster than methods based on temporal differencing. Its disadvantages are that it does not apply to all types of material rheology, and it can also require much storage when many snapshots and time sections are desired. Comparisons between numerical and analytical solutions for simple acoustic and elastic problems demonstrate the high accuracy of the REM.     

PROCESSING OF PS-REFLECTION DATA APPLYING A COMMON CONVERSION-POINT STACKING TECHNIQUE1

Converted waves require special data processing as the wave paths are asymmetrical. The CMP concept is not applicable for converted PS waves, instead a sorting algorithm for a common conversion point (CCP) has to be applied. The coordinates of the conversion points in a single homogeneous layer can be calculated as a function of the offset, the reflector depth and the velocity ratio v_P/ v_s. For multilayered media, an approximation for the coordinates of the conversion points can be made. Numerical tests show that the traveltime of PS reflections can be approximated with sufficient accuracy for a certain offset range by a two-term series truncation. Therefore NMO corrections can be calculated by standard routines which use the hyperbolic approximation of the reflection traveltime curves. The CCP-stacking technique has been applied to field data which have been generated by three vertical vibrators. The in-line horizontal components have been recorded. The static corrections have been estimated from additional P- and SH-wave measurements for the source and geophone locations, respectively. The data quality has been improved by processes such as spectral balancing. A comparison with the stacked results of the corresponding P- and SH-wavefield surveys shows a good coherency of structural features in P-, SH- and PS-time sections.     

COMMON REFLECTION POINT DATA-STACKING TECHNIQUE FOR CONVERTED WAVES1

For converted waves stacking requires a true common reflection point gather which, in this case, is also a common conversion point (CCP) gather. We consider converted waves of the PS- and SP-type in a stack of horizontal layers. The coordinates of the conversion points for waves of PS- or SP-type, respectively, in a single homogeneous layer are calculated as a function of the offset, the reflector depth and the velocity ratio v_p/v_s. Knowledge of the conversion points enables us to gather the seismic traces in a common conversion point (CCP) record. Numerical tests show that the CCP coordinates in a multilayered medium can be approximated by the equations given for a single layer. In practical applications, an a priori estimate of v_p/v_s is required to obtain the CCP for a given reflector depth. A series expansion for the traveltime of converted waves as a function of the offset is presented. Numerical examples have been calculated for several truncations. For small offsets, a hyperbolic approximation can be used. For this, the rms velocity of converted waves is defined. A Dix-type formula, relating the product of the interval velocities of compressional and shear waves to the rms velocity of the converted waves, is presented.