Computation of additional components of the first-order ray approximation in isotropic media

This study shows that the use of the first-order additional components of the ray method in the seismic wave field modeling is easy and that it can bring a substantial improvement of the standard ray results obtained with the zero-order ray approximation only. For the calculation of a first-order additional component, spatial derivatives of the parameters of the medium and spatial derivatives of the zero-order ray amplitude term are necessary. The evaluation of the former derivatives is straightforward; the latter derivatives can be calculated approximately from neighboring rays by substituting the derivatives by finite differences. This allows an effective calculation of the first-order additional terms in arbitrary laterally varying layered media.The importance of the first-order additional terms is demonstrated by the study of individual higher-order terms of the ray series representing elementaryP andS elastodynamic Green functions for a homogeneous isotropic medium. The study shows clearly that the consideration of the first-order additional terms leads to a more substantial decrease of the difference between approximate and exact elementary Green functions than any other higher-order term. With this in mind, effects of the first-order additional terms on the ray synthetic seismograms for aVSP configuration are studied. It is shown that the use of the additional terms leads to such phenomena, unknown in the zero-order approximation of the ray method, like quasi-elliptical and transverse polarization of a singleP wave or longitudinal polarization of a singleS wave.