scholarly journals Migration Using a Transversely Isotropic Medium with Symmetry Normal to the Reflector Dip

2011 ◽  
Vol 2011 ◽  
pp. 1-5 ◽  
Author(s):  
Tariq Alkhalifah ◽  
Paul Sava

A transversely isotropic (TI) model in which the tilt is constrained to be normal to the dip (DTI model) allows for simplifications in the imaging and velocity model building efforts as compared to a general TI (TTI) model. Although this model cannot be represented physically in all situations, for example, in the case of conflicting dips, it handles arbitrary reflector orientations under the assumption of symmetry axis normal to the dip. Using this assumption, we obtain efficient downward continuation algorithms compared to the general TTI ones, by utilizing the reflection features of such a model. Phase-shift migration can be easily extended to approximately handle lateral inhomogeneity using, for example, the split-step approach. This is possible because, unlike the general TTI case, the DTI model reduces to VTI for zero dip. These features enable a process in which we can extract velocity information by including tools that expose inaccuracies in the velocity model in the downward continuation process. We test this model on synthetic data corresponding to a general TTI medium and show its resilience.

Geophysics ◽  
2011 ◽  
Vol 76 (3) ◽  
pp. WA13-WA21 ◽  
Author(s):  
Mamoru Takanashi ◽  
Ilya Tsvankin

Nonhyperbolic moveout analysis plays an increasingly important role in velocity model building because it provides valuable information for anisotropic parameter estimation. However, lateral heterogeneity associated with stratigraphic lenses such as channels and reefs can significantly distort the moveout parameters, even when the structure is relatively simple. We analyze the influence of a low-velocity isotropic lens on nonhyperbolic moveout inversion for horizontally layered VTI (transversely isotropic with a vertical symmetry axis) models. Synthetic tests demonstrate that a lens can cause substantial, laterally varying errors in the normal-moveout velocity [Formula: see text] and the anellipticity parameter [Formula: see text]. The area influenced by the lens can be identified using the residual moveout after the nonhyperbolic moveout correction as well as the dependence of errors in [Formula: see text] and [Formula: see text] on spreadlength. To remove such errors in [Formula: see text] and [Formula: see text], we propose a correction algorithm designed for a lens embedded in a horizontally layered overburden. This algorithm involves estimation of the incidence angle of the ray passing through the lens for each recorded trace. With the assumption that lens-related perturbation of the raypath is negligible, the lens-induced traveltime shifts are computed from the corresponding zero-offset time distortion (i.e., from “pull-up” or “push-down” anomalies). Synthetic tests demonstrate that this algorithm substantially reduces the errors in the effective and interval parameters [Formula: see text] and [Formula: see text]. The corrected traces and reconstructed “background” values of [Formula: see text] and [Formula: see text] are suitable for anisotropic time imaging and producing a high-quality stack.


Geophysics ◽  
2010 ◽  
Vol 75 (4) ◽  
pp. D27-D36 ◽  
Author(s):  
Andrey Bakulin ◽  
Marta Woodward ◽  
Dave Nichols ◽  
Konstantin Osypov ◽  
Olga Zdraveva

Tilted transverse isotropy (TTI) is increasingly recognized as a more geologically plausible description of anisotropy in sedimentary formations than vertical transverse isotropy (VTI). Although model-building approaches for VTI media are well understood, similar approaches for TTI media are in their infancy, even when the symmetry-axis direction is assumed known. We describe a tomographic approach that builds localized anisotropic models by jointly inverting surface-seismic and well data. We present a synthetic data example of anisotropic tomography applied to a layered TTI model with a symmetry-axis tilt of 45 degrees. We demonstrate three scenarios for constraining the solution. In the first scenario, velocity along the symmetry axis is known and tomography inverts for Thomsen’s [Formula: see text] and [Formula: see text] parame-ters. In the second scenario, tomography inverts for [Formula: see text], [Formula: see text], and velocity, using surface-seismic data and vertical check-shot traveltimes. In contrast to the VTI case, both these inversions are nonunique. To combat nonuniqueness, in the third scenario, we supplement check-shot and seismic data with the [Formula: see text] profile from an offset well. This allows recovery of the correct profiles for velocity along the symmetry axis and [Formula: see text]. We conclude that TTI is more ambiguous than VTI for model building. Additional well data or rock-physics assumptions may be required to constrain the tomography and arrive at geologically plausible TTI models. Furthermore, we demonstrate that VTI models with atypical Thomsen parameters can also fit the same joint seismic and check-shot data set. In this case, although imaging with VTI models can focus the TTI data and match vertical event depths, it leads to substantial lateral mispositioning of the reflections.


Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. WB191-WB207 ◽  
Author(s):  
Yaxun Tang ◽  
Biondo Biondi

We present a new strategy for efficient wave-equation migration-velocity analysis in complex geological settings. The proposed strategy has two main steps: simulating a new data set using an initial unfocused image and performing wavefield-based tomography using this data set. We demonstrated that the new data set can be synthesized by using generalized Born wavefield modeling for a specific target region where velocities are inaccurate. We also showed that the new data set can be much smaller than the original one because of the target-oriented modeling strategy, but it contains necessary velocity information for successful velocity analysis. These interesting features make this new data set suitable for target-oriented, fast and interactive velocity model-building. We demonstrate the performance of our method on both a synthetic data set and a field data set acquired from the Gulf of Mexico, where we update the subsalt velocity in a target-oriented fashion and obtain a subsalt image with improved continuities, signal-to-noise ratio and flattened angle-domain common-image gathers.


Geophysics ◽  
2013 ◽  
Vol 78 (1) ◽  
pp. S25-S36 ◽  
Author(s):  
Ernesto V. Oropeza ◽  
George A. McMechan

An efficient Kirchhoff-style prestack depth migration, called “parsimonious” migration, was developed a decade ago for isotropic 2D and 3D media by using measured slownesses to reduce the amount of ray tracing by orders of magnitude. It is conceptually similar to “map” migration, but its implementation has some differences. We have extended this approach to 2D tilted transversely isotropic (TTI) media and illustrated it with synthetic P-wave data. Although the framework of isotropic parsimonious may be retained, the extension to TTI media requires redevelopment of each of the numerical components, calculation of the phase and group velocity for TTI media, development of a new two-point anisotropic ray tracer, and substitution of an initial-angle isotropic shooting ray-trace algorithm for an anisotropic one. The model parameterization consists of Thomsen’s parameters ([Formula: see text], [Formula: see text], [Formula: see text]) and the tilt angle of the symmetry axis of the TI medium. The parsimonious anisotropic migration algorithm is successfully applied to synthetic data from a TTI version of the Marmousi2 model. The quality of the image improves by weighting the impulse response by the calculation of the anisotropic Fresnel radius. The accuracy and speed of this migration makes it useful for anisotropic velocity model building. The elapsed computing time for 101 shots for the Marmousi2 TTI model is 35 s per shot (each with 501 traces) in 32 Opteron cores.


Geophysics ◽  
2010 ◽  
Vol 75 (3) ◽  
pp. A19-A24 ◽  
Author(s):  
Tariq Alkhalifah ◽  
Paul Sava

The computational tools for imaging in transversely isotropic media with tilted axes of symmetry (TTI) are complex and in most cases do not have an explicit closed-form representation. Developing such tools for a TTI medium with tilt constrained to be normal to the reflector dip (DTI) reduces their complexity and allows for closed-form representations. The homogeneous-case zero-offset migration in such a medium can be performed using an isotropic operator scaled by the velocity of the medium in the tilt direction. For the nonzero-offset case, the reflection angle is always equal to the incidence angle, and thus, the velocities for the source and receiver waves at the reflection point are equal and explicitly dependent on the reflection angle. This fact allows for the development of explicit representations for angle decomposition as well as moveout formulas for analysis of extended images obtained by wave-equation migration. Although setting the tilt normal to the reflector dip may not be valid everywhere (i.e., on salt flanks), it can be used in the process of velocity model building, in which such constrains are useful and typically are used.


2018 ◽  
Vol 6 (4) ◽  
pp. SM63-SM70 ◽  
Author(s):  
Tian Jun ◽  
Peng Gengxin ◽  
Junru Jiao ◽  
Grace (Yan) Yan ◽  
Xianhuai Zhu

A special challenge for land seismic exploration is estimating velocities, in part due to complex near-surface structures, and in some instances because of rugose topography over foothills. We have developed an integrated turning-ray and reflection-tomographic method to face this challenge. First, turning-ray tomography is performed to derive a near-surface velocity-depth model. Then, we combine the near-surface model with the initial-subsurface model. Taking the combined model as starting model, we go through a reflection tomographic process to build the model for imaging. During reflection tomography, the near-surface model and subsurface models are jointly updated. Our method has been successfully applied to a 2D complex synthetic data example and a 3D field data example. The results demonstrate that our method derives a very decent model even when there is no reflection information available in a few hundred meters underneath the surface. Joint tomography can lead to geologic plausible models and produce subsurface images with high fidelity.


Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE39-VE47 ◽  
Author(s):  
Jessé C. Costa ◽  
Francisco J. da Silva ◽  
Ellen N. Gomes ◽  
Jörg Schleicher ◽  
Luiz André Melo ◽  
...  

Seismic imaging in depth is limited by the accuracy of velocity model estimation. Slope tomography uses the slowness components and traveltimes of picked reflection or diffraction events for velocity model building. The unavoidable data incompleteness requires additional information to assure stability to inversion. One natural constraint for ray-based tomography is a smooth velocity model. We propose a new, reflection-angle-based kind of smoothness constraint as regularization in slope tomography and have compared its effects to three other, more conventional constraints. The effects of these constraints were evaluated through angle-domain common-image gathers, computed with wave-equation migration using the estimated velocity model. We found that the smoothness constraints have a distinct effect on the velocity model but a weaker effect on the migrated data. In numerical tests on synthetic data, the new constraint leads to geologically more consistent models.


Geophysics ◽  
2018 ◽  
Vol 83 (1) ◽  
pp. S47-S55 ◽  
Author(s):  
Parsa Bakhtiari Rad ◽  
Benjamin Schwarz ◽  
Dirk Gajewski ◽  
Claudia Vanelle

Diffraction imaging can lead to high-resolution characterization of small-scale subsurface structures. A key step of diffraction imaging and tomography is diffraction separation and enhancement, especially in the full prestack data volume. We have considered point diffractors and developed a robust and fully data-driven workflow for prestack diffraction separation based on wavefront attributes, which are determined using the common-reflection-surface (CRS) method. In the first of two steps, we apply a zero-offset-based extrapolation operator for prestack diffraction separation, which combines the robustness and stability of the zero-offset CRS processing with enhanced resolution and improved illumination of the finite-offset CRS processing. In the second step, when the finite-offset diffracted events are separated, we apply a diffraction-based time migration velocity model building that provides high-quality diffraction velocity spectra. Applications of the new workflow to 2D/3D complex synthetic data confirm the superiority of prestack diffraction separation over the poststack method as well as the high potential of diffractions for improved time imaging.


Sign in / Sign up

Export Citation Format

Share Document