Building tilted transversely isotropic depth models using localized anisotropic tomography with well information

Geophysics ◽  
2010 ◽  
Vol 75 (4) ◽  
pp. D27-D36 ◽  
Author(s):  
Andrey Bakulin ◽  
Marta Woodward ◽  
Dave Nichols ◽  
Konstantin Osypov ◽  
Olga Zdraveva
Keyword(s):  
Seismic Data ◽  
Model Building ◽  
Symmetry Axis ◽  
Transverse Isotropy ◽  
Synthetic Data ◽  
Rock Physics ◽  
Transversely Isotropic ◽  
Data Set ◽  
Anisotropic Models ◽  
Well Data

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.

Localized anisotropic tomography with well information in VTI media

Geophysics ◽  
2010 ◽  
Vol 75 (5) ◽  
pp. D37-D45 ◽  
Author(s):  
Andrey Bakulin ◽  
Marta Woodward ◽  
Dave Nichols ◽  
Konstantin Osypov ◽  
Olga Zdraveva
Keyword(s):  
Vertical Velocity ◽  
Seismic Data ◽  
Model Building ◽  
Synthetic Data ◽  
Resolution Limit ◽  
Data Set ◽  
Anisotropic Models ◽  
Velocity Models ◽  
First Case ◽  
Well Data

We develop a concept of localized seismic grid tomography constrained by well information and apply it to building vertically transversely isotropic (VTI) velocity models in depth. The goal is to use a highly automated migration velocity analysis to build anisotropic models that combine optimal image focusing with accurate depth positioning in one step. We localize tomography to a limited volume around the well and jointly invert the surface seismic and well data. Well information is propagated into the local volume by using the method of preconditioning, whereby model updates are shaped to follow geologic layers with spatial smoothing constraints. We analyze our concept with a synthetic data example of anisotropic tomography applied to a 1D VTI model. We demonstrate four cases of introducing additionalinformation. In the first case, vertical velocity is assumed to be known, and the tomography inverts only for Thomsen’s [Formula: see text] and [Formula: see text] profiles using surface seismic data alone. In the second case, tomography simultaneously inverts for all three VTI parameters, including vertical velocity, using a joint data set that consists of surface seismic data and vertical check-shot traveltimes. In the third and fourth cases, sparse depth markers and walkaway vertical seismic profiling (VSP) are used, respectively, to supplement the seismic data. For all four examples, tomography reliably recovers the anisotropic velocity field up to a vertical resolution comparable to that of the well data. Even though walkaway VSP has the additional dimension of angle or offset, it offers no further increase in this resolution limit. Anisotropic tomography with well constraints has multiple advantages over other approaches and deserves a place in the portfolio of model-building tools.

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.

scholarly journals Skeletonized wave-equation inversion in vertical symmetry axis media without too much math

2017 ◽  
Vol 5 (3) ◽  
pp. SO21-SO30 ◽  
Author(s):  
Shihang Feng ◽  
Gerard T. Schuster
Keyword(s):  
Wave Equation ◽  
Symmetry Axis ◽  
Waveform Inversion ◽  
Synthetic Data ◽  
Full Waveform Inversion ◽  
Data Set ◽  
Traveltime Inversion ◽  
Anisotropic Models ◽  
Full Waveform ◽  
Starting Model

We have developed a tutorial for skeletonized inversion of pseudo-acoustic anisotropic vertical symmetry axis (VTI) data. We first invert for the anisotropic models using wave-equation traveltime inversion. Here, the skeletonized data are the traveltimes of transmitted and/or reflected arrivals that lead to simpler misfit functions and more robust convergence compared with full-waveform inversion. This provides a good starting model for waveform inversion. The effectiveness of this procedure is illustrated with synthetic data examples and a marine data set recorded in the Gulf of Mexico.

Velocity analysis for tilted transversely isotropic media: A physical modeling example

Geophysics ◽  
2001 ◽  
Vol 66 (3) ◽  
pp. 904-910 ◽  
Author(s):  
Vladimir Grechka ◽  
Andres Pech ◽  
Ilya Tsvankin ◽  
Baoniu Han
Keyword(s):  
Physical Modeling ◽  
Symmetry Axis ◽  
Transverse Isotropy ◽  
Transversely Isotropic ◽  
P Wave ◽  
Velocity Analysis ◽  
Isotropic Layer ◽  
Anisotropic Parameter ◽  
Data Set ◽  
Symmetry Direction

Transverse isotropy with a tilted symmetry axis (TTI media) has been recognized as a common feature of shale formations in overthrust areas, such as the Canadian Foothills. Since TTI layers cause serious problems in conventional imaging, it is important to be able to reconstruct the velocity model suitable for anisotropic depth migration. Here, we discuss the results of anisotropic parameter estimation on a physical‐modeling data set. The model represents a simplified version of a typical overthrust section from the Alberta Foothills, with a horizontal reflector overlaid by a bending transversely isotropic layer. Assuming that the TTI layer is homogeneous and the symmetry axis stays perpendicular to its boundaries, we invert P-wave normal‐moveout (NMO) velocities and zero‐offset traveltimes for the symmetry‐direction velocity V0 and the anisotropic parameters ε and δ. The coefficient ε is obtained using the traveltimes of a wave that crosses a dipping TTI block and reflects from the bottom of the model. The inversion for ε is based on analytic expressions for NMO velocity in media with intermediate dipping interfaces. Our estimates of both anisotropic coefficients are close to their actual values. The errors in the inversion, which are associated primarily with the uncertainties in picking the NMO velocities and traveltimes, can be reduced by a straighforward modification of the acquisition geometry. It should be emphasized that the moveout inversion also gives an accurate estimate of the thickness of the TTI layer, thus reconstructing the correct depth scale of the section. Although the physical model used here was relatively simple, our results demonstrate the principal feasibility of anisotropic velocity analysis and imaging in overthrust areas. The main problems in anisotropic processing for TTI models are likely to be caused by the lateral variation of the velocity field and overall structural complexity.

Multiparameter joint tomography for TTI model building

Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. WB183-WB190 ◽  
Author(s):  
Chaoguang Zhou ◽  
Junru Jiao ◽  
Sonny Lin ◽  
John Sherwood ◽  
Sverre Brandsberg-Dahl
Keyword(s):  
Model Building ◽  
Symmetry Axis ◽  
Anisotropy Parameter ◽  
Limited Range ◽  
Single Parameter ◽  
Symmetry Direction ◽  
Anisotropic Models ◽  
Axis Direction ◽  
Anisotropy Parameters ◽  
Well Data

Model building for tilted transversely isotropic media has commonly been performed by a single parameter tomography that updates the velocity in the symmetry direction, while the orientation of the symmetry axis and Thomsen parameters [Formula: see text] and [Formula: see text] are typically estimated from the migration stack and well data. Unfortunately, well data are often not available. In addition, when they are available, their lateral sampling is typically very sparse and their vertical sampling usually spans only a limited range of depths. In order to obtain spatially varying anisotropic models, with or without well data, we developed a multiparameter joint tomographic approach that simultaneously inverts for the velocity in the symmetry axis direction, [Formula: see text] and [Formula: see text]. We derived a set of reflection tomography equations for slowness in the symmetry axis direction and Thomsen parameters [Formula: see text] and [Formula: see text]. In order to address the nonuniqueness of the tomography, we developed a regularization strategy that uses an independent regularization operator and regularization factor for each individual anisotropy parameter. Synthetic tests found that ambiguity exists between the anisotropy parameters and that velocity has a better resolution than [Formula: see text] and [Formula: see text]. They also confirmed that joint tomography provides a better data fit than single parameter tomography. The field example was used to test a way to incorporate the sonic data in the model building process and limit the tomographic updates on certain anisotropy parameters by adjusting the regularization.

Inversion of ultrasonic data for transversely isotropic media

Geophysics ◽  
2017 ◽  
Vol 82 (1) ◽  
pp. C1-C7 ◽  
Author(s):  
Yevhen Kovalyshen ◽  
Joel Sarout ◽  
Jeremie Dautriat
Keyword(s):  
Numerical Algorithm ◽  
Seismic Data ◽  
Symmetry Axis ◽  
Rock Sample ◽  
Transverse Isotropy ◽  
Transversely Isotropic ◽  
P Wave ◽  
Wave Velocities ◽  
Transversely Isotropic Media ◽  
Isotropic Media

We have developed a new numerical algorithm for inversion of ultrasonic data in transversely isotropic media. This algorithm is able to determine from the measured P-wave velocities the orientation of the symmetry axis of a rock sample and the Thomsen’s parameters, only assuming transverse isotropy. The inversion of ultrasonic data acquired on natural and potentially heterogeneous shale samples produced reasonable results. In addition, the algorithm was successfully tested on ultrasonic data acquired on synthetic samples with predefined orientations of the symmetry axis. An additional outcome of the algorithm is a simple approximation of Thomsen’s formulation, which can be effectively used for interpretation of seismic data in transversely isotropic media.

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
Keyword(s):  
Model Building ◽  
Symmetry Axis ◽  
Synthetic Data ◽  
Transversely Isotropic ◽  
Velocity Model ◽  
Downward Continuation ◽  
Velocity Model Building ◽  
Lateral Inhomogeneity ◽  
Velocity Information ◽  
Model Phase

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.

Wave-equation migration velocity analysis for VTI models

Geophysics ◽  
2014 ◽  
Vol 79 (3) ◽  
pp. WA59-WA68 ◽  
Author(s):  
Yunyue Li ◽  
Biondo Biondi ◽  
Robert Clapp ◽  
Dave Nichols
Keyword(s):  
Wave Equation ◽  
Synthetic Data ◽  
Rock Physics ◽  
Migration Velocity ◽  
Velocity Analysis ◽  
Sources Of Information ◽  
Data Set ◽  
Anisotropic Models ◽  
Migration Velocity Analysis ◽  
Power Maximization

Anisotropic models are needed for wave simulation and inversion where a complex geologic environment exists. We extended the theory of wave equation migration velocity analysis to build vertical transverse isotropic models. Because of the ambiguity between depth and [Formula: see text] in the acoustic regime, we assumed [Formula: see text] can be accurately obtained from other sources of information, and inverted for the NMO slowness and the anellipticity parameter [Formula: see text]. We combined the differential semblance optimization objective function with the stacking power maximization to evaluate the focusing of the prestack image in the subsurface-offset domain. To regularize the multiparameter inversion, we built a framework to adapt the geologic and the rock physics information to guide the updates in NMO slowness and [Formula: see text]. This regularization step was crucial to stabilize the inversion and to produce geologically meaningful results. We tested the proposed approach on a synthetic data set and a 2D Gulf of Mexico data set starting with a fairly good initial anisotropic model. The inversion results revealed shallow anomalies collocated in NMO velocity and [Formula: see text] and improved the continuity and the resolution of the final stacked images.

Estimating frequency-dependent attenuation quality factor values from prestack surface seismic data

Geophysics ◽  
2017 ◽  
Vol 82 (1) ◽  
pp. O11-O22 ◽  
Author(s):  
James Beckwith ◽  
Roger Clark ◽  
Linda Hodgson
Keyword(s):  
Quality Factor ◽  
Seismic Data ◽  
Synthetic Data ◽  
Rock Physics ◽  
Analytic Form ◽  
Ocean Bottom ◽  
Data Set ◽  
Frequency Dependent ◽  
Seismic Quality Factor ◽  
Frequency Independent

The intrinsic seismic quality factor [Formula: see text] is known from poroelastic rock-physics theory to be frequency dependent, even within typical bandwidths of individual surface- and borehole-based surveys in which measurement methods usually deliver frequency-independent [Formula: see text]. Thus, measuring frequency-dependent [Formula: see text] instead offers better characterization of seismic properties and moreover a potential step toward estimating permeability directly from seismic data. Therefore, we have introduced a method to measure frequency-dependent [Formula: see text] from pairs of reflections in prestack [Formula: see text]-[Formula: see text] domain surface seismic data — a data type that, unlike a vertical seismic profile, offers useful areal coverage. Although, in principle, any analytic form with a manageable number of parameters could be prescribed, the frequency dependence of [Formula: see text] is modeled as a power law, [Formula: see text]. Inversion is done with a simple grid search over coefficient ([Formula: see text]) and exponent [Formula: see text], seeking a minimum [Formula: see text]-norm. We have found, using a numerical experiment and a synthetic data set, that it is robust and also accurate down to a signal-to-noise ratio of approximately 0.65. Then, [Formula: see text] is estimated for some 955 [Formula: see text] superbins of a 3D prestack ocean bottom cable data set over the Kinnoull field, central North Sea. All combinations of eight prominent reflections between Top Miocene and Base Cretaceous were treated together to give some 21,000 frequency-dependent and (for comparison) constant-[Formula: see text] results. The median coefficient ([Formula: see text]) and exponent [Formula: see text] were 0.0074 and 0.06, respectively, with sharply peaked distributions (excess kurtosis [Formula: see text]). Outlier, strongly frequency-dependent results, given by [Formula: see text], coincide with low-frequency “shadows” under amplitude anomalies, adversely affecting the spectra of reflections. The inferred frequency-dependent [Formula: see text] at 32.5 Hz, the center of the available bandwidth, is not statistically different from the frequency-independent [Formula: see text], 181 with a standard error from the distribution of one, derived from the same data. Hence for this data set, a constant-[Formula: see text] assumption would in fact be adequate. However, our method has the ability to measure stable estimates of [Formula: see text].

Regularization and datuming of seismic data by weighted, damped least squares

Geophysics ◽  
2006 ◽  
Vol 71 (5) ◽  
pp. U67-U76 ◽  
Author(s):  
Robert J. Ferguson
Keyword(s):  
Least Squares ◽  
Seismic Data ◽  
Synthetic Data ◽  
Real Data ◽  
Structural Data ◽  
Irregular Sampling ◽  
Data Set ◽  
Near Surface ◽  
Damped Least Squares ◽  
Wavefield Extrapolation

The possibility of improving regularization/datuming of seismic data is investigated by treating wavefield extrapolation as an inversion problem. Weighted, damped least squares is then used to produce the regularized/datumed wavefield. Regularization/datuming is extremely costly because of computing the Hessian, so an efficient approximation is introduced. Approximation is achieved by computing a limited number of diagonals in the operators involved. Real and synthetic data examples demonstrate the utility of this approach. For synthetic data, regularization/datuming is demonstrated for large extrapolation distances using a highly irregular recording array. Without approximation, regularization/datuming returns a regularized wavefield with reduced operator artifacts when compared to a nonregularizing method such as generalized phase shift plus interpolation (PSPI). Approximate regularization/datuming returns a regularized wavefield for approximately two orders of magnitude less in cost; but it is dip limited, though in a controllable way, compared to the full method. The Foothills structural data set, a freely available data set from the Rocky Mountains of Canada, demonstrates application to real data. The data have highly irregular sampling along the shot coordinate, and they suffer from significant near-surface effects. Approximate regularization/datuming returns common receiver data that are superior in appearance compared to conventional datuming.

Sign in / Sign up

Export Citation Format

Share Document