Analytic study of the effective parameters for determination of the NMO velocity function in transversely isotropic media

Author(s):  
Jack K. Cohen
Geophysics ◽  
1997 ◽  
Vol 62 (6) ◽  
pp. 1855-1866 ◽  
Author(s):  
Jack K. Cohen

In their studies of transversely isotropic media with a vertical symmetry axis (VTI media), Alkhalifah and Tsvankin observed that, to a high numerical accuracy, the normal moveout (NMO) velocity for dipping reflectors as a function of ray parameter p depends mainly on just two parameters, each of which can be determined from surface P‐wave observations. They substantiated this result by using the weak‐anisotropy approximation and exploited it to develop a time‐domain processing sequence that takes into account vertical transverse isotropy. In this study, the two‐parameter Alkhalifah‐Tsvankin result was further examined analytically. It was found that although there is (as these authors already observed) some dependence on the remaining parameters of the problem, this dependence is weak, especially in the practically important regimes of weak to moderately strong transverse isotropy and small ray parameter. In each of these regimes, an analytic solution is derived for the anisotropy parameter η required for time‐domain P‐wave imaging in VTI media. In the case of elliptical anisotropy (η = 0), NMO velocity expressed through p is fully controlled just by the zero‐dip NMO velocity—one of the Alkhalifah‐ Tsvankin parameters. The two‐parameter representation of NMO velocity also was shown to be exact in another limit—that of the zero shear‐wave vertical velociy. The analytic results derived here are based on new representations for both the P‐wave phase velocity and normal moveout velocity in terms of the ray parameter, with explicit expressions given for the cases of vanishing onaxis shear speed, weak to moderate transverse isotropy, and small to moderate ray parameter. Using these formulas, I have rederived and, in some cases, extended in a uniform manner various results of Tsvankin, Alkhalifah, and others. Examples include second‐order expansions in the anisotropy parameters for both the P‐wave phase‐velocity function and NMO‐velocity function, as well as expansions in powers of the ray parameter for both of these functions. I have checked these expansions against the corresponding exact functions for several choices of the anisotropy parameters.


Geophysics ◽  
2014 ◽  
Vol 79 (1) ◽  
pp. C19-C26 ◽  
Author(s):  
Alexey Stovas ◽  
Tariq Alkhalifah

Moveout approximations play a very important role in seismic modeling, inversion, and scanning for parameters in complex media. We developed a scheme to map one-way moveout approximations for transversely isotropic media with a vertical axis of symmetry (VTI), which is widely available, to the tilted case (TTI) by introducing the effective tilt angle. As a result, we obtained highly accurate TTI moveout equations analogous with their VTI counterparts. Our analysis showed that the most accurate approximation is obtained from the mapping of generalized approximation. The new moveout approximations allow for, as the examples demonstrate, accurate description of moveout in the TTI case even for vertical heterogeneity. The proposed moveout approximations can be easily used for inversion in a layered TTI medium because the parameters of these approximations explicitly depend on corresponding effective parameters in a layered VTI medium.


Geophysics ◽  
1998 ◽  
Vol 63 (1) ◽  
pp. 275-278 ◽  
Author(s):  
Jack K. Cohen

Tsvankin (1995) derived an expression for the normal‐moveout (NMO) velocity from horizontal or dipping reflectors that is valid in symmetry planes of anisotropic media. This equation is written in terms of the phase velocity and its derivatives with respect to the phase angle with vertical. The phase angle, however, is not a seismic observable, while the ray parameter p is. Indeed, Alkhalifah and Tsvankin (1995) pointed out the importance of representing the NMO velocity function as a function of ray parameter and obtained a number of significant results via an indirect use of this variable. Subsequently, Cohen (1997) directly transformed the Tsvankin expression for the moveout velocity obtaining an expression involving the phase velocity expressed as a function of ray parameter. Applying this result to the normal moveout of P-waves in transversely isotropic media with a vertical symmetry axis (VTI media), he was able to get various approximations for the ray‐parameter form of the moveout function and thus give additional analytic insight into the success of the velocity‐analysis methodology Alkhalifah and Tsvankin.


Geophysics ◽  
1996 ◽  
Vol 61 (3) ◽  
pp. 835-845 ◽  
Author(s):  
John Anderson ◽  
Tariq Alkhalifah ◽  
Ilya Tsvankin

The main advantage of Fowler’s dip‐moveout (DMO) method is the ability to perform velocity analysis along with the DMO removal. This feature of Fowler DMO becomes even more attractive in anisotropic media, where imaging methods are hampered by the difficulty in reconstructing the velocity field from surface data. We have devised a Fowler‐type DMO algorithm for transversely isotropic media using the analytic expression for normal‐moveout velocity. The parameter‐estimation procedure is based on the results of Alkhalifah and Tsvankin showing that in transversely isotropic media with a vertical axis of symmetry (VTI) the P‐wave normal‐moveout (NMO) velocity as a function of ray parameter can be described fully by just two coefficients: the zero‐dip NMO velocity [Formula: see text] and the anisotropic parameter η (η reduces to the difference between Thomsen parameters ε and δ in the limit of weak anisotropy). In this extension of Fowler DMO, resampling in the frequency‐wavenumber domain makes it possible to obtain the values of [Formula: see text] and η by inspecting zero‐offset (stacked) panels for different pairs of the two parameters. Since most of the computing time is spent on generating constant‐velocity stacks, the added computational effort caused by the presence of anisotropy is relatively minor. Synthetic and field‐data examples demonstrate that the isotropic Fowler DMO technique fails to generate an accurate zero‐offset section and to obtain the zero‐dip NMO velocity for nonelliptical VTI models. In contrast, this anisotropic algorithm allows one to find the values of the parameters [Formula: see text] and η (sufficient to perform time migration as well) and to correct for the influence of transverse isotropy in the DMO processing. When combined with poststack F-K Stolt migration, this method represents a complete inversion‐processing sequence capable of recovering the effective parameters of transversely isotropic media and producing migrated images for the best‐fit homogeneous anisotropic model.


2020 ◽  
Vol 38 (2) ◽  
Author(s):  
Razec Cezar Sampaio Pinto da Silva Torres ◽  
Leandro Di Bartolo

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.


1995 ◽  
Vol 26 (2-3) ◽  
pp. 431-436 ◽  
Author(s):  
Patrick N.(Jr). Okoye ◽  
N. F. Uren ◽  
W. Waluyo

Sign in / Sign up

Export Citation Format

Share Document