Approximations for Far Ray-Tracing through Layered Anisotropic Media with Velocity Gradients

Author(s):  
P. R. Williamson ◽  
B. Raynaud
1991 ◽  
Author(s):  
J. Costa ◽  
M. Schoenberg ◽  
D. Miller

Geophysics ◽  
2018 ◽  
Vol 83 (5) ◽  
pp. T281-T289 ◽  
Author(s):  
Qianru Xu ◽  
Weijian Mao

We have developed a fast ray-tracing method for multiple layered inhomogeneous anisotropic media, based on the generalized Snell’s law. Realistic geologic structures continuously varying with embedded discontinuities are parameterized by adopting cubic B-splines with nonuniformly spaced nodes. Because the anisotropic characteristic is often closely related to the interface configuration, this model parameterization scheme containing the natural inclination of the corresponding layer is particularly suitable for tilted transverse isotropic models whose symmetry axis is generally perpendicular to the direction of the layers. With this model parameterization, the first- and second-order spatial derivatives of the velocity within the interfaces can be effectively obtained, which facilitates the amplitude computation in dynamic ray tracing. By using complex initial conditions for the dynamic ray system and taking the multipath effect into consideration, our method is applicable to Gaussian beam migration. Numerical experiments of our method have been used to verify its effectiveness, practicability, and efficiency in memory storage and computation.


Geophysics ◽  
1996 ◽  
Vol 61 (4) ◽  
pp. 947-963 ◽  
Author(s):  
Tariq Alkhalifah

Nearly all dip‐moveout correction (DMO) implementations to date assume isotropic homogeneous media. Usually, this has been acceptable considering the tremendous cost savings of homogeneous isotropic DMO and considering the difficulty of obtaining the anisotropy parameters required for effective implementation. In the presence of typical anisotropy, however, ignoring the anisotropy can yield inadequate results. Since anisotropy may introduce large deviations from hyperbolic moveout, accurate transformation to zero‐offset in anisotropic media should address such nonhyperbolic moveout behavior of reflections. Artley and Hale’s v(z) ray‐tracing‐based DMO, developed for isotropic media, provides an attractive approach to treating such problems. By using a ray‐tracing procedure crafted for anisotropic media, I modify some aspects of their DMO so that it can work for v(z) anisotropic media. DMO impulse responses in typical transversely isotropic (TI) models (such as those associated with shales) deviate substantially from the familiar elliptical shape associated with responses in homogeneous isotropic media (to the extent that triplications arise even where the medium is homogeneous). Such deviations can exceed those caused by vertical inhomogeneity, thus emphasizing the importance of taking anisotropy into account in DMO processing. For isotropic or elliptically anisotropic media, the impulse response is an ellipse; but as the key anisotropy parameter η varies, the shape of the response differs substantially from elliptical. For typical η > 0, the impulse response in TI media tends to broaden compared to the response in an isotropic homogeneous medium, a behavior opposite to that encountered in typical v(z) isotropic media, where the response tends to be squeezed. Furthermore, the amplitude distribution along the DMO operator differs significantly from that for isotropic media. Application of this anisotropic DMO to data from offshore Africa resulted in a considerably better alignment of reflections from horizontal and dipping reflectors in common‐midpoint gather than that obtained using an isotropic DMO. Even the presence of vertical inhomogeneity in this medium could not eliminate the importance of considering the shale‐induced anisotropy.


Geophysics ◽  
2020 ◽  
Vol 85 (6) ◽  
pp. T331-T342
Author(s):  
Xing-Wang Li ◽  
Bing Zhou ◽  
Chao-Ying Bai ◽  
Jian-Lu Wu

In a viscoelastic anisotropic medium, velocity anisotropy and wave energy attenuation occur and are often observed in seismic data applications. Numerical investigation of seismic wave propagation in complex viscoelastic anisotropic media is very helpful in understanding seismic data and reconstructing subsurface structures. Seismic ray tracing is an effective means to study the propagation characteristics of high-frequency seismic waves. Unfortunately, most seismic ray-tracing methods and traveltime tomographic inversion algorithms only deal with elastic media and ignore the effect of viscoelasticity on the seismic raypath. We have developed a method to find the complex ray velocity that gives the seismic ray speed and attenuation in an arbitrary viscoelastic anisotropic medium, and we incorporate them with the modified shortest-path method to determine the raypath and calculate the real and imaginary traveltime (wave energy attenuation) simultaneously. We determine that the complex ray-tracing method is applicable to arbitrary 2D/3D viscoelastic anisotropic media in a complex geologic model and the computational errors of the real and imaginary traveltime are less than 0.36% and 0.59%, respectively. The numerical examples verify that the new method is an effective and powerful tool for accomplishing seismic complex ray tracing in heterogeneous viscoelastic anisotropic media.


Geophysics ◽  
2007 ◽  
Vol 72 (3) ◽  
pp. S133-S138 ◽  
Author(s):  
Tianfei Zhu ◽  
Samuel H. Gray ◽  
Daoliu Wang

Gaussian-beam depth migration is a useful alternative to Kirchhoff and wave-equation migrations. It overcomes the limitations of Kirchhoff migration in imaging multipathing arrivals, while retaining its efficiency and its capability of imaging steep dips with turning waves. Extension of this migration method to anisotropic media has, however, been hampered by the difficulties in traditional kinematic and dynamic ray-tracing systems in inhomogeneous, anisotropic media. Formulated in terms of elastic parameters, the traditional anisotropic ray-tracing systems aredifficult to implement and inefficient for computation, especially for the dynamic ray-tracing system. They may also result inambiguity in specifying elastic parameters for a given medium.To overcome these difficulties, we have reformulated the ray-tracing systems in terms of phase velocity.These reformulated systems are simple and especially useful for general transversely isotropic and weak orthorhombic media, because the phase velocities for these two types of media can be computed with simple analytic expressions. These two types of media also represent the majority of anisotropy observed in sedimentary rocks. Based on these newly developed ray-tracing systems, we have extended prestack Gaussian-beam depth migration to general transversely isotropic media. Test results with synthetic data show that our anisotropic, prestack Gaussian-beam migration is accurate and efficient. It produces images superior to those generated by anisotropic, prestack Kirchhoff migration.


Geophysics ◽  
2005 ◽  
Vol 70 (6) ◽  
pp. D65-D75 ◽  
Author(s):  
Ivan Pšenčík ◽  
Véronique Farra

We propose approximate ray-tracing equations for qP-waves propagating in smooth, inhomogeneous, weakly anisotropic media. For their derivation, we use perturbation theory, in which deviations of anisotropy from isotropy are considered to be the first-order quantities. The proposed ray-tracing equations and corresponding traveltimes are of the first order. Accuracy of the traveltimes can be increased by calculating a secondorder correction along first-order rays. The first-order ray-tracing equations for qP-waves propagating in a general weakly anisotropic medium depend on only 15 weak-anisotropy parameters (generalization of Thomsen’s parameters). The equations are thus considerably simpler than the exact ray-tracing equations. For higher-symmetry anisotropic media the equations differ only slightly from equations for isotropic media. They can thus substitute for the traditional isotropic ray tracers used in seismic processing. For vanishing anisotropy, the first-order ray-tracing equations reduce to standard, exact ray-tracing equations for isotropic media. Numerical tests for configuration and models used in seismic prospecting indicate negligible dependence of accuracy of calculated traveltimes on inhomogeneity of the medium. For anisotropy of about 8%, considered in the examples presented, the relative errors of the traveltimes, including the second-order correction, are well under 0.05%; for anisotropy of about 20%, they do not exceed 0.3%.


Sign in / Sign up

Export Citation Format

Share Document