Fowler DMO and time migration in transversely isotropic media

Author(s):  
John E. Anderson ◽  
Tariq A. Alkhalifah ◽  
Ilya Tsvankin
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.


Geophysics ◽  
2006 ◽  
Vol 71 (1) ◽  
pp. S29-S35 ◽  
Author(s):  
Tariq Alkhalifah

Using a newly developed nonhyperbolic offset-mid-point traveltime equation for prestack Kirchhoff time migration, instead of the conventional double-square-root (DSR) equation, results in overall better images from anisotropic data. Specifically, prestack Kirchhoff time migration for transversely isotropic media with a vertical symmetry axis (VTI media) is implemented using an analytical offset-midpoint traveltime equation that represents the equivalent of Cheop's pyramid for VTI media. It includes higher-order terms necessary to better handle anisotropy as well as vertical inhomogeneity. Application of this enhanced Kirchhoff time-migration method to the anisotropic Marmousi data set demonstrates the effectiveness of the approach. Further application of the method to field data from Trinidad results in sharper reflectivity images of the subsurface, with the faults better focused and positioned than with images obtained using isotropic methods. The superiority of the anisotropic time migration is evident in the flatness of the image gathers.


Geophysics ◽  
2000 ◽  
Vol 65 (4) ◽  
pp. 1316-1325 ◽  
Author(s):  
Tariq Alkhalifah

Prestack Kirchhoff time migration for transversely isotropic media with a vertical symmetry axis (VTI media) is implemented using an offset‐midpoint traveltime equation, Cheop’s pyramid equivalent equation for VTI media. The derivation of such an equation for VTI media requires approximations that pertain to high frequency and weak anisotropy. Yet the resultant offset‐midpoint traveltime equation for VTI media is highly accurate for even strong anisotropy. It is also strictly dependent on two parameters: NMO velocity and the anisotropy parameter, η. It reduces to the exact offset‐midpoint traveltime equation for isotropic media when η = 0. In vertically inhomogeneous media, the NMO velocity and η parameters in the offset‐midpoint traveltime equation are replaced by their effective values: the velocity is replaced by the rms velocity and η is given by a more complicated equation that includes summation of the fourth power of velocity.


Geophysics ◽  
2020 ◽  
pp. 1-79
Author(s):  
Ali Fathalian ◽  
Daniel O. Trad ◽  
Kristopher A. Innanen

Anisotropy and absorption are critical to the modeling and analysis of seismic amplitude,phase, and traveltime data. To neglect any of these phenomena, which are often bothoperating simultaneously, degrades the resolution and interpretability of migrated images.However, a full accounting of anisotropy and anelasticity is computationally complex andexpensive. One strategy for accommodating these aspects of wave propagation, while keepingcost and complexity under control, is to do so within an acoustic approximation. Weset up a procedure for solving the time-domain viscoacoustic wave equation for tilted transverselyisotropic (TTI) media, based on a standard linear solid model and, from this, developa viscoacoustic reverse time migration (Q-RTM) algorithm. In this approach, amplitudecompensation occurs within the migration process through a manipulation of attenuationand phase dispersion terms in the time domain differential equations. Specifically, theback-propagation operator is constructed by reversing the sign only of the amplitude lossoperators, but not the dispersion-related operators, a step made possible by reformulatingthe absorptive TTI equations such that the loss and dispersion operators appear separately.The scheme is tested on synthetic examples to examine the capacity of viscoacoustic RTM to correct for attenuation, and the overall stability of the procedure.


Geophysics ◽  
1993 ◽  
Vol 58 (10) ◽  
pp. 1454-1467 ◽  
Author(s):  
Ken L. Larner ◽  
Jack K. Cohen

Given the sensitivity of imaging accuracy to the velocity used in migration, migration founded (as in practice) on the erroneous assumption that a medium is isotropic can be expected to be inaccurate for steep reflectors. Here, we estimate errors in interpreted reflection time and lateral position as a function of reflector dip for transversely isotropic models in which the axis of symmetry is vertical and the medium velocity varies linearly with depth. We limit consideration to media in which ratios of the various elastic moduli are independent of depth. Tests with reflector dips up to 120 degrees on a variety of anisotropic media show errors that are tens of wavelengths for dips beyond 90 degrees when the medium (unrealistically) is homogeneous. For a given anisotropy, the errors are smaller for inhomogeneous media; the larger the velocity gradient, the smaller the errors. For gradients that are representative of the subsurface, lateral‐position errors tend to be minor for dips less than about 60 degrees, growing to two to five wavelengths as dip passes beyond 90 degrees. These errors depend on reflector depth and average velocity to the reflector only through their ratio, i.e., migrated reflection time. Migration error, which is found to be unrelated to the ratio of horizontal to vertical velocity, is such that reflections with later migrated reflection times tend to be more severely overmigrated than are those with earlier times. Over a large range of dips, migration errors that arise when anisotropy is ignored but inhomogeneity is honored tend to be considerably smaller than those encountered when inhomogeneity is ignored in migrating data from isotropic, inhomogeneous media.


Geophysics ◽  
1998 ◽  
Vol 63 (2) ◽  
pp. 623-631 ◽  
Author(s):  
Tariq Alkhalifah

When transversely isotropic (VTI) media with vertical symmetry axes are characterized using the zero‐dip normal moveout (NMO) velocity [[Formula: see text]] and the anisotropy parameter ηinstead of Thomsen’s parameters, time‐related processing [moveout correction, dip moveout (DMO), and time migration] become nearly independent of the vertical P- and S-wave velocities ([Formula: see text] and [Formula: see text], respectively). The independence on [Formula: see text] and [Formula: see text] is well within the limits of seismic accuracy, even for relatively strong anisotropy. The dependency on [Formula: see text] and [Formula: see text] reduces even further as the ratio [Formula: see text] decreases. In fact, for [Formula: see text], all time‐related processing depends exactly on only [Formula: see text] and η. This fortunate dependence on two parameters is demonstrated here through analytical derivations of time‐related processing equations in terms of [Formula: see text] and η. The time‐migration dispersion relation, the NMO velocity for dipping events, and the ray‐tracing equations extracted by setting [Formula: see text] (i.e., by considering VTI as acoustic) not only depend solely on [Formula: see text] and η but are much simpler than the counterpart expressions for elastic media. Errors attributed to this use of the acoustic assumption are small and may be neglected. Therefore, as in isotropic media, the acoustic model arising from setting [Formula: see text], although not exactly true for VTI media, can serve as a useful approximation to the elastic model for the kinematics of P-wave data. This approximation can boost the efficiency of imaging and DMO programs for VTI media as well as simplify their description.


Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. S383-S398 ◽  
Author(s):  
Chenlong Wang ◽  
Jiubing Cheng ◽  
Børge Arntsen

Recording P- and S-wave modes acquires more information related to rock properties of the earth’s interior. Elastic migration, as a part of multicomponent seismic data processing, potentially offers a great improvement over conventional acoustic migration to create a spatial image of some medium properties. In the framework of elastic reverse time migration, we have developed new scalar and vector imaging conditions assisted by efficient polarization-based mode decoupling to avoid crosstalk among the different wave modes for isotropic and transversely isotropic media. For the scalar imaging, we corrected polarity reversal of zero-lag PS images using the local angular attributes on the fly of angle-domain imaging. For the vector imaging, we naturally used the polarization information in the decoupled single-mode vector fields to automatically avoid the polarity reversal and to estimate the local angular attributes for angle-domain imaging. Examples of increasing complexity in 2D and 3D cases found that the proposed approaches can be used to obtain a physically interpretable image and angle-domain common-image gather at an acceptable computational cost. Decoupling and imaging the 3D S-waves involves some complexity, which has not been addressed in the literature. For this reason, we also attempted at illustrating the physical contents of the two separated S-wave modes and their contribution to seismic full-wave imaging.


Sign in / Sign up

Export Citation Format

Share Document