Calculation of slowness vectors from ray directions for qP-, qSV-, and qSH-waves in tilted transversely isotropic media

Geophysics ◽  
2018 ◽  
Vol 83 (4) ◽  
pp. C153-C160 ◽  
Author(s):  
Lijiao Zhang ◽  
Bing Zhou
Keyword(s):  
Ray Tracing ◽  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Seismic Wave Propagation ◽  
Previous Method ◽  
New Methods ◽  
Slowness Vector ◽  
Isotropic Media ◽  
Tti Media ◽  
Wave Modes

Kinematic ray tracing is an effective way to simulate the seismic wave propagation in isotropic and anisotropic media. It is essential to know the ray velocity when tracing seismic rays. But in anisotropic media, the ray velocity is a function of the direction of the slowness vector instead of the ray direction and it often deviates from the phase velocity. In this case, it causes a critical problem for ray tracing, which is how to calculate the ray velocity from a known ray direction. If we could calculate the phase slowness vector from ray directions, the ray velocity could be computed. We have evaluated a previous method in the first place. Then, we developed two new methods to solve two existing problems of the previous method: (1) It leads to complex and multiple solutions of the slowness vector and (2) it mixes up the qP- and qSV-wave modes. Our first method solves the two problems by applying eigenvalues to separate the wave modes and decrease the two unknowns ([Formula: see text] and [Formula: see text]) to only one unknown in two equations. Our second method is based on the general relationship between the slowness and ray-velocity vectors and shows that only one unknown is involved in one equation for tilted transversely isotropic (TTI) media. After obtaining the slowness vector, the ray velocity can be computed easily. A 2D model is designed to test the feasibility of the new methods. Using the results for the model, we found that the presented approaches were applicable for ray tracing in TTI media.

scholarly journals Transformation to zero offset in transversely isotropic media

Geophysics ◽  
1996 ◽  
Vol 61 (4) ◽  
pp. 947-963 ◽  
Author(s):  
Tariq Alkhalifah
Keyword(s):  
Ray Tracing ◽  
Impulse Response ◽  
Homogeneous Medium ◽  
Anisotropy Parameter ◽  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Amplitude Distribution ◽  
Isotropic Media ◽  
Zero Offset ◽  
Vertical Inhomogeneity

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.

Effective and efficient approaches for calculating seismic ray velocity and attenuation in viscoelastic anisotropic media

Geophysics ◽  
2020 ◽  
Vol 86 (1) ◽  
pp. C19-C35
Author(s):  
Jianlu Wu ◽  
Bing Zhou ◽  
Xingwang Li ◽  
Youcef Bouzidi
Keyword(s):  
Velocity Vector ◽  
Hamiltonian Function ◽  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Quality Factors ◽  
Approximate Methods ◽  
Slowness Vector ◽  
Homogeneous Complex ◽  
Isotropic Media ◽  
Complex Ray

In viscoelastic anisotropic media, the elastic moduli, slowness vector, phase, and ray velocity are all complex-valued quantities in the frequency domain. Solving the complex eikonal equation becomes computationally complex and time-consuming. We have developed two approximate methods to effectively calculate the ray velocity vector, attenuation, and quality factor in viscoelastic transversely isotropic media with a vertical symmetry axis (VTI) and in orthorhombic (ORT) anisotropy. The first method is based on the perturbation theory (PER) under the assumption of a homogeneous complex ray vector, which is obtained by applying the elastic background and viscoelastic perturbations to the real and imaginary components of the modulus tensor, respectively. The perturbations of the slowness vectors of the three wave modes (qP, qSV, and qSH) are determined through the vanishing Hamiltonian function. The second method is derived by applying a real slowness direction (RSD) to the inhomogeneous complex slowness vector and then approximately calculating the complex ray velocity vector with the condition of the homogeneous complex vector. The numerical results verify that the two approaches can produce accurate ray velocity vector, attenuation, and quality factors of the qP-wave in viscoelastic VTI and ORT media. The RSD method can yield high accuracies of ray velocity for the qSV- and qSH-wave in viscoelastic VTI models even at triplication of the qSV wavefronts, as well as qS1 and qS2 in a weak ORT medium ([Formula: see text] > 20), except for near the cusp of the qS1 wavefronts (errors approximately 6%) where the PER has more than 10% error.

Gaussian beam depth migration for anisotropic media

Geophysics ◽  
1995 ◽  
Vol 60 (5) ◽  
pp. 1474-1484 ◽  
Author(s):  
Tariq Alkhalifah
Keyword(s):  
Gaussian Beam ◽  
Ray Tracing ◽  
Anisotropy Parameter ◽  
Synthetic Data ◽  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Inhomogeneous Media ◽  
Computational Effort ◽  
Isotropic Media ◽  
Gaussian Beam Migration

Gaussian beam migration (GBM), as it is implemented today, efficiently handles isotropic inhomogeneous media. The approach is based on the solution of the wave equation in ray‐centered coordinates. Here, I extend the method to work for 2-D migration in generally anisotropic inhomogeneous media. Extension of the Gaussian‐beam method from isotropic to anisotropic media involves modification of the kinematics and dynamics in the required ray tracing. While the accuracy of the paraxial expansion for anisotropic media is comparable to that for isotropic media, ray tracing in anisotropic media is much slower than that in isotropic media. However, because ray tracing is just a small portion of the computation in GBM, the increased computational effort in general anisotropic GBM is typically only about 40%. Application of this method to synthetic examples shows successful migration in inhomogeneous, transversely isotropic media for reflector dips up to and beyond 90°. Further applications to synthetic data of layered anisotropic media show the importance of applying the proper smoothing to the velocity field used in the migration. Also, tests with synthetic data show that the quality of anisotropic migration of steep events in a medium with velocity increasing with depth is much more sensitive to the Thomsen anisotropy parameter ε than to the parameter δ. Thus, a good estimate of ε is needed to apply anisotropic migration with confidence.

Reverse time migration in tilted transversely isotropic media

2020 ◽  
Vol 38 (2) ◽  
Author(s):  
Razec Cezar Sampaio Pinto da Silva Torres ◽  
Leandro Di Bartolo
Keyword(s):  
Seismic Anisotropy ◽  
Synthetic Data ◽  
Transversely Isotropic ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Computer Clusters ◽  
Time Migration ◽  
Transversely Isotropic Media ◽  
Isotropic Media ◽  
Tti Media

ABSTRACT. Reverse time migration (RTM) is one of the most powerful methods used to generate images of the subsurface. The RTM was proposed in the early 1980s, but only recently it has been routinely used in exploratory projects involving complex geology – Brazilian pre-salt, for example. Because the method uses the two-way wave equation, RTM is able to correctly image any kind of geological environment (simple or complex), including those with anisotropy. On the other hand, RTM is computationally expensive and requires the use of computer clusters. This paper proposes to investigate the influence of anisotropy on seismic imaging through the application of RTM for tilted transversely isotropic (TTI) media in pre-stack synthetic data. This work presents in detail how to implement RTM for TTI media, addressing the main issues and specific details, e.g., the computational resources required. A couple of simple models results are presented, including the application to a BP TTI 2007 benchmark model.Keywords: finite differences, wave numerical modeling, seismic anisotropy. Migração reversa no tempo em meios transversalmente isotrópicos inclinadosRESUMO. A migração reversa no tempo (RTM) é um dos mais poderosos métodos utilizados para gerar imagens da subsuperfície. A RTM foi proposta no início da década de 80, mas apenas recentemente tem sido rotineiramente utilizada em projetos exploratórios envolvendo geologia complexa, em especial no pré-sal brasileiro. Por ser um método que utiliza a equação completa da onda, qualquer configuração do meio geológico pode ser corretamente tratada, em especial na presença de anisotropia. Por outro lado, a RTM é dispendiosa computacionalmente e requer o uso de clusters de computadores por parte da indústria. Este artigo apresenta em detalhes uma implementação da RTM para meios transversalmente isotrópicos inclinados (TTI), abordando as principais dificuldades na sua implementação, além dos recursos computacionais exigidos. O algoritmo desenvolvido é aplicado a casos simples e a um benchmark padrão, conhecido como BP TTI 2007.Palavras-chave: diferenças finitas, modelagem numérica de ondas, anisotropia sísmica.

Transverse isotropy of the upper mantle in the vicinity of Pacific fracture zones

1969 ◽  
Vol 59 (1) ◽  
pp. 59-72
Author(s):  
Robert S. Crosson ◽  
Nikolas I. Christensen
Keyword(s):  
Upper Mantle ◽  
Symmetry Axis ◽  
Transverse Isotropy ◽  
Transversely Isotropic ◽  
Seismic Wave Propagation ◽  
Fracture Zones ◽  
Mantle Material ◽  
Transversely Isotropic Media ◽  
Horizontal Symmetry ◽  
Isotropic Media

Abstract Several recent investigations suggest that portions of the Earth's upper mantle behave anisotropically to seismic wave propagation. Since several types of anisotropy can produce azimuthal variations in Pn velocities, it is of particular geophysical interest to provide a framework for the recognition of the form or forms of anisotropy most likely to be manifest in the upper mantle. In this paper upper mantle material is assumed to possess the elastic properties of transversely isotropic media. Equations are presented which relate azimuthal variations in Pn velocities to the direction and angle of tilt of the symmetry axis of a transversely isotropic upper mantle. It is shown that the velocity data of Raitt and Shor taken near the Mendocino and Molokai fracture zones can be adequately explained by the assumption of transverse isotropy with a nearly horizontal symmetry axis.

Normal moveout from dipping reflectors in anisotropic media

Geophysics ◽  
1995 ◽  
Vol 60 (1) ◽  
pp. 268-284 ◽  
Author(s):  
Ilya Tsvankin
Keyword(s):  
Symmetry Axis ◽  
Transverse Isotropy ◽  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Weak Anisotropy ◽  
Anisotropic Models ◽  
Transversely Isotropic Media ◽  
Wide Range ◽  
Isotropic Media ◽  
The Difference

Description of reflection moveout from dipping interfaces is important in developing seismic processing methods for anisotropic media, as well as in the inversion of reflection data. Here, I present a concise analytic expression for normal‐moveout (NMO) velocities valid for a wide range of homogeneous anisotropic models including transverse isotropy with a tilted in‐plane symmetry axis and symmetry planes in orthorhombic media. In transversely isotropic media, NMO velocity for quasi‐P‐waves may deviate substantially from the isotropic cosine‐of‐dip dependence used in conventional constant‐velocity dip‐moveout (DMO) algorithms. However, numerical studies of NMO velocities have revealed no apparent correlation between the conventional measures of anisotropy and errors in the cosine‐of‐dip DMO correction (“DMO errors”). The analytic treatment developed here shows that for transverse isotropy with a vertical symmetry axis, the magnitude of DMO errors is dependent primarily on the difference between Thomsen parameters ε and δ. For the most common case, ε − δ > 0, the cosine‐of‐dip–corrected moveout velocity remains significantly larger than the moveout velocity for a horizontal reflector. DMO errors at a dip of 45 degrees may exceed 20–25 percent, even for weak anisotropy. By comparing analytically derived NMO velocities with moveout velocities calculated on finite spreads, I analyze anisotropy‐induced deviations from hyperbolic moveout for dipping reflectors. For transversely isotropic media with a vertical velocity gradient and typical (positive) values of the difference ε − δ, inhomogeneity tends to reduce (sometimes significantly) the influence of anisotropy on the dip dependence of moveout velocity.

Migration of P‐wave reflection data in transversely isotropic media

Geophysics ◽  
1994 ◽  
Vol 59 (4) ◽  
pp. 591-596 ◽  
Author(s):  
Suhas Phadke ◽  
S. Kapotas ◽  
N. Dai ◽  
Ernest R. Kanasewich
Keyword(s):  
Wave Equation ◽  
Dispersion Relation ◽  
Synthetic Data ◽  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Limited Range ◽  
Horizontal Velocity ◽  
P Wave ◽  
Transversely Isotropic Media ◽  
Isotropic Media

Wave propagation in transversely isotropic media is governed by the horizontal and vertical wave velocities. The quasi‐P(qP) wavefront is not an ellipse; therefore, the propagation cannot be described by the wave equation appropriate for elliptically anisotropic media. However, for a limited range of angles from the vertical, the dispersion relation for qP‐waves can be approximated by an ellipse. The horizontal velocity necessary for this approximation is different from the true horizontal velocity and depends upon the physical properties of the media. In the method described here, seismic data is migrated using a 45-degree wave equation for elliptically anisotropic media with the horizontal velocity determined by comparing the 45-degree elliptical dispersion relation and the quasi‐P‐dispersion relation. The method is demonstrated for some synthetic data sets.

Prestack Gaussian-beam depth migration in anisotropic media

Geophysics ◽  
2007 ◽  
Vol 72 (3) ◽  
pp. S133-S138 ◽  
Author(s):  
Tianfei Zhu ◽  
Samuel H. Gray ◽  
Daoliu Wang
Keyword(s):  
Gaussian Beam ◽  
Ray Tracing ◽  
Synthetic Data ◽  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Elastic Parameters ◽  
Depth Migration ◽  
Kirchhoff Migration ◽  
Beam Depth ◽  
Dynamic Ray Tracing

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.

scholarly journals 3-D two‐point ray tracing for heterogeneous, weakly transversely isotropic media

Geophysics ◽  
1996 ◽  
Vol 61 (6) ◽  
pp. 1883-1894 ◽  
Author(s):  
Vladimir Y. Grechka ◽  
George A. McMechan
Keyword(s):  
Ray Tracing ◽  
Conjugate Gradient Method ◽  
Conjugate Gradient ◽  
Gradient Method ◽  
Heterogeneous Media ◽  
Chebyshev Approximation ◽  
Transversely Isotropic ◽  
Transversely Isotropic Media ◽  
Isotropic Media ◽  
Derivatives Of

A two‐point ray‐tracing technique for 3-D smoothly heterogeneous, weakly transversely isotropic media is based on Fermat’s principle and takes advantage of global Chebyshev approximation of both the model and curved rays. This approximation gives explicit relations for derivatives of traveltime with respect to ray parameters and allows use of the rapidly converging conjugate gradient method to compute traveltimes. The method is fast because, for most smoothly heterogeneous media, approximation of rays by only a few polynomials and a few conjugate gradient iterations provide excellent precision in traveltime calculation.

Sign in / Sign up

Export Citation Format

Share Document