Depth‐domain velocity analysis in VTI media using surface P-wave data: Is it feasible?

Geophysics ◽  
2001 ◽  
Vol 66 (3) ◽  
pp. 897-903 ◽  
Author(s):  
Yves Le Stunff ◽  
Vladimir Grechka ◽  
Ilya Tsvankin
Keyword(s):  
Vertical Velocity ◽  
Heterogeneous Media ◽  
Wave Reflection ◽  
P Wave ◽  
Model Parameters ◽  
Velocity Analysis ◽  
Depth Migration ◽  
Velocity Models ◽  
Trade Offs ◽  
Vti Media

The main difficulties in anisotropic velocity analysis and inversion using surface seismic data are associated with the multiparameter nature of the problem and inherent trade‐offs between the model parameters. For the most common anisotropic model, transverse isotropy with a vertical symmetry axis (VTI media), P-wave kinematic signatures are controlled by the vertical velocity V0 and the anisotropic parameters ε and δ. However, only two combinations of these parameters—NMO velocity from a horizontal reflector Vnmo(0) and the anellipticity coefficient η—can be determined from P-wave reflection traveltimes if the medium above the reflector is laterally homogeneous. While Vnmo(0) and η are sufficient for time‐domain imaging in VTI media, they cannot be used to resolve the vertical velocity and build velocity models needed for depth migration. Here, we demonstrate that P-wave reflection data can be inverted for all three relevant VTI parameters (V0, ε and δ) if the model contains nonhorizontal intermediate interfaces. Using anisotropic reflection tomography, we carry out parameter estimation for a two‐layer medium with a curved intermediate interface and reconstruct the correct anisotropic depth model. To explain the success of this inversion procedure, we present an analytic study of reflection traveltimes for this model and show that the information about the vertical velocity and reflector depth was contained in the reflected rays which crossed the dipping intermediate interface. The results of this work are especially encouraging because the need for depth imaging (such as prestack depth migration) arises mostly in laterally heterogeneous media. Still, we restricted this study to a relatively simple model and constrained the inversion by assuming that one of the layers is isotropic. In general, although lateral heterogeneity does create a dependence of P-wave reflection traveltimes on the vertical velocity, there is no guarantee that for more complicated models all anisotropic parameters can be resolved in a unique fashion.

Image gathers of SV-waves in homogeneous and factorized VTI media

Geophysics ◽  
2005 ◽  
Vol 70 (5) ◽  
pp. D55-D64 ◽  
Author(s):  
Ramzy M. Al-Zayer ◽  
Ilya Tsvankin
Keyword(s):  
Shear Wave ◽  
Vertical Velocity ◽  
Migration Velocity ◽  
P Wave ◽  
Model Parameters ◽  
Velocity Analysis ◽  
Wave Data ◽  
Migration Velocity Analysis ◽  
Wave Migration ◽  
Vti Media

Reflection moveout of SV-waves in transversely isotropic media with a vertical symmetry axis (VTI media) can provide valuable information about the model parameters and help to overcome the ambiguities in the inversion of P-wave data. Here, to develop a foundation for shear-wave migration velocity analysis, we study SV-wave image gathers obtained after prestack depth migration. The key issue, addressed using both approximate analytic results and Kirchhoff migration of synthetic data, is whether long-spread SV data can constrain the shear-wave vertical velocity [Formula: see text] and the depth scale of VTI models. For homogeneous media, the residual moveout of horizontal SV events on image gathers is close to hyperbolic and depends just on the NMO velocity [Formula: see text] out to offset-to-depth ratios of about 1.7. Because [Formula: see text] differs from [Formula: see text], flattening moderate-spread gathers of SV-waves does not ensure the correct depth of the migrated events. The residual moveout rapidly becomes nonhyperbolic as the offset-to-depth ratio approaches two, with the migrated depths at long offsets strongly influenced by the SV-wave anisotropy parameter σ. Although the combination of [Formula: see text] and σ is sufficient to constrain the vertical velocity [Formula: see text] and reflector depth, the tradeoff between σ and the Thomsen parameter ε on long-spread gathers causes errors in time-to-depth conversion. The residual moveout of dipping SV events is also controlled by the parameters [Formula: see text], σ, and ε, but in the presence of dip, the contributions of both σ and ε are significant even at small offsets. For factorized v(z) VTI media with a constant SV-wave vertical-velocity gradient [Formula: see text], flattening of horizontal events for a range of depths requires the correct NMO velocity at the surface, the gradient [Formula: see text], and, for long offsets, the parameters σ and ε. On the whole, the nonnegligible uncertainty in the estimation of reflector depth from SV-wave moveout highlights the need to combine P- and SV-wave data in migration velocity analysis for VTI media.

Migration velocity analysis in factorized VTI media

Geophysics ◽  
2004 ◽  
Vol 69 (3) ◽  
pp. 708-718 ◽  
Author(s):  
Debashish Sarkar ◽  
Ilya Tsvankin
Keyword(s):  
Vertical Velocity ◽  
Velocity Gradient ◽  
Migration Velocity ◽  
P Wave ◽  
Velocity Analysis ◽  
Prestack Depth Migration ◽  
Depth Migration ◽  
Migration Velocity Analysis ◽  
Wave Migration ◽  
Vti Media

One of the main challenges in anisotropic velocity analysis and imaging is simultaneous estimation of velocity gradients and anisotropic parameters from reflection data. Approximating the subsurface by a factorized VTI (transversely isotropic with a vertical symmetry axis) medium provides a convenient way of building vertically and laterally heterogeneous anisotropic models for prestack depthmigration. The algorithm for P‐wave migration velocity analysis (MVA) introduced here is designed for models composed of factorized VTI layers or blocks with constant vertical and lateral gradients in the vertical velocity VP0. The anisotropic MVA method is implemented as an iterative two‐step procedure that includes prestack depth migration (imaging step) followed by an update of the medium parameters (velocity‐analysis step). The residual moveout of the migrated events, which is minimized during the parameter updates, is described by a nonhyperbolic equation whose coefficients are determined by 2D semblance scanning. For piecewise‐factorized VTI media without significant dips in the overburden, the residual moveout of P‐wave events in image gathers is governed by four effective quantities in each block: (1) the normal‐moveout velocity Vnmo at a certain point within the block, (2) the vertical velocity gradient kz, (3) the combination kx[Formula: see text] of the lateral velocity gradient kx and the anisotropic parameter δ, and (4) the anellipticity parameter η. We show that all four parameters can be estimated from the residual moveout for at least two reflectors within a block sufficiently separated in depth. Inversion for the parameter η also requires using either long‐spread data (with the maximum offset‐to‐depth ratio no less than two) from horizontal interfaces or reflections from dipping interfaces. To find the depth scale of the section and build a model for prestack depth migration using the MVA results, the vertical velocity VP0 needs to be specified for at least a single point in each block. When no borehole information about VP0 is available, a well‐focused image can often be obtained by assuming that the vertical‐velocity field is continuous across layer boundaries. A synthetic test for a three‐layer model with a syncline structure confirms the accuracy of our MVA algorithm in estimating the interval parameters Vnmo, kz, kx, and η and illustrates the influence of errors in the vertical velocity on the image quality.

Analysis of image gathers in factorized VTI media

Geophysics ◽  
2003 ◽  
Vol 68 (6) ◽  
pp. 2016-2025 ◽  
Author(s):  
Debashish Sarkar ◽  
Ilya Tsvankin
Keyword(s):  
Vertical Velocity ◽  
Velocity Gradient ◽  
Migration Velocity ◽  
P Wave ◽  
Velocity Variation ◽  
Velocity Analysis ◽  
Lateral Heterogeneity ◽  
Velocity Models ◽  
Isotropic Media ◽  
Vti Media

Because events in image gathers generated after prestack depth migration are sensitive to the velocity field, they are often used in migration velocity analysis for isotropic media. Here, we present an analytic and numerical study of P‐wave image gathers in transversely isotropic media with a vertical symmetry axis (VTI) and establish the conditions for flattening such events and positioning them at the true reflector depth. Application of the weak‐anisotropy approximation leads to concise expressions for reflections in image gathers from homogeneous and factorized v(z) media in terms of the VTI parameters and the vertical velocity gradient kz. Flattening events in image gathers for any reflector dip requires accurate values of the zero‐dip NMO velocity at the surface [Vnmo (z = 0)], the gradient kz, and the anellipticity coefficient η. For a fixed error in Vnmo and kz, the magnitude of residual moveout of events in image gathers decreases with dip, while the moveout caused by an error in η initially increases for moderate dips but then decreases as dips approach 90°. Flat events in image gathers in VTI media, however, do not guarantee the correct depth scale of the model because reflector depth depends on the vertical migration velocity. For factorized v(x, z) media with a linear velocity variation in both the x‐ and z‐directions, the moveout on image gathers is controlled by Vnmo (x = z = 0), kz, η, and a combination of the horizontal velocity gradient kx and the Thomsen parameter δ (specifically, kx[Formula: see text]). If too large a value of any of these four quantities is used in migration, reflections in the image gathers curve downward (i.e., they are undercorrected; the inferred depth increases with offset), while a negative error results in overcorrection. Lateral heterogeneity tends to increase the sensitivity of moveout of events in image gathers to the parameter η, and errors in η may lead to measurable residual moveout of horizontal events in v(x, z) media even for offset‐to‐depth ratios close to unity. These results provide a basis for extending to VTI media conventional velocity analysis methods operating with image gathers. Although P‐wave traveltimes alone cannot be used to separate anisotropy from lateral heterogeneity (i.e., kx is coupled to δ), moveout of events in image gathers does constrain the vertical gradient kz. Hence, it may be possible to build VTI velocity models in depth by supplementing reflection data with minimal a priori information, such as the vertical velocity at the top of the factorized VTI layer.

Migration moveout analysis and depth focusing

Geophysics ◽  
1993 ◽  
Vol 58 (1) ◽  
pp. 91-100 ◽  
Author(s):  
Claude F. Lafond ◽  
Alan R. Levander
Keyword(s):  
Heterogeneous Media ◽  
A Priori ◽  
Complex Structure ◽  
Tomographic Reconstruction ◽  
Synthetic Data ◽  
Real Data ◽  
Velocity Model ◽  
Velocity Analysis ◽  
Data Set ◽  
Velocity Models

Prestack depth migration still suffers from the problems associated with building appropriate velocity models. The two main after‐migration, before‐stack velocity analysis techniques currently used, depth focusing and residual moveout correction, have found good use in many applications but have also shown their limitations in the case of very complex structures. To address this issue, we have extended the residual moveout analysis technique to the general case of heterogeneous velocity fields and steep dips, while keeping the algorithm robust enough to be of practical use on real data. Our method is not based on analytic expressions for the moveouts and requires no a priori knowledge of the model, but instead uses geometrical ray tracing in heterogeneous media, layer‐stripping migration, and local wavefront analysis to compute residual velocity corrections. These corrections are back projected into the velocity model along raypaths in a way that is similar to tomographic reconstruction. While this approach is more general than existing migration velocity analysis implementations, it is also much more computer intensive and is best used locally around a particularly complex structure. We demonstrate the technique using synthetic data from a model with strong velocity gradients and then apply it to a marine data set to improve the positioning of a major fault.

Prestack depth migration of an Alberta Foothills data set—The Husky experience

Geophysics ◽  
1998 ◽  
Vol 63 (2) ◽  
pp. 392-398 ◽  
Author(s):  
W.-J. Wu ◽  
L. Lines ◽  
A. Burton ◽  
H.-X. Lu ◽  
J. Zhu ◽  
...  
Keyword(s):  
Velocity Analysis ◽  
Reverse Time ◽  
Prestack Depth Migration ◽  
Data Set ◽  
Depth Migration ◽  
Geological Interpretation ◽  
Prestack Migration ◽  
Velocity Models ◽  
Depth Images ◽  
Migration Velocity Analysis

We produce depth images for an Alberta Foothills line by iteratively using a number of migration and velocity analysis techniques. In imaging steeply dipping layers of a foothills data set, it is apparent that thrust belt geology can violate the conventional assumptions of elevation datum corrections and common midpoint (CMP) stacking. To circumvent these problems, we use migration from topography in which we perform prestack depth migration on the data using correct source and receiver elevations. Migration from topography produces enhanced images of steep shallow reflectors when compared to conventional processing. In addition to migration from topography, we couple prestack depth migration with the continuous adjustment of velocity depth models. A number of criteria are used in doing this. These criteria require that our velocity estimates produce a focused image and that migrated depths in common image gathers be independent of source‐receiver offset. Velocity models are estimated by a series of iterative and interpretive steps involving prestack migration velocity analysis and structural interpretation. Overlays of velocity models on depth migrations should generally show consistency between velocity boundaries and reflection depths. Our preferred seismic depth section has been produced by using prestack reverse‐time depth migration coupled with careful geological interpretation.

Waveform inversion for microseismic velocity analysis and event location in VTI media

Geophysics ◽  
2017 ◽  
Vol 82 (4) ◽  
pp. WA95-WA103 ◽  
Author(s):  
Oscar Jarillo Michel ◽  
Ilya Tsvankin
Keyword(s):  
Moment Tensor ◽  
Source Parameters ◽  
Waveform Inversion ◽  
Transversely Isotropic ◽  
Velocity Analysis ◽  
Origin Time ◽  
Velocity Models ◽  
Trade Offs ◽  
Receiver Array ◽  
Event Location

Waveform inversion (WI), which has been extensively used in reflection seismology, could provide improved velocity models and event locations for microseismic surveys. Here, we develop an elastic WI algorithm for anisotropic media designed to estimate the 2D velocity field along with the source parameters (location, origin time, and moment tensor) from microseismic data. The gradient of the objective function is obtained with the adjoint-state method, which requires just two modeling simulations at each iteration. In the current implementation the source coordinates and velocity parameters are estimated sequentially at each stage of the inversion to minimize trade-offs and improve the convergence. Synthetic examples illustrate the accuracy of the inversion for layered VTI (transversely isotropic with a vertical symmetry axis) media, as well as the sensitivity of the velocity-analysis results to noise, the length of the receiver array, errors in the initial model, and variability in the moment tensor of the recorded events.

Elastic reflection traveltime inversion with decoupled wave equation

Geophysics ◽  
2018 ◽  
Vol 83 (5) ◽  
pp. R463-R474 ◽  
Author(s):  
Guanchao Wang ◽  
Shangxu Wang ◽  
Jianyong Song ◽  
Chunhui Dong ◽  
Mingqiang Zhang
Keyword(s):  
Wave Equation ◽  
Wave Velocity ◽  
Waveform Inversion ◽  
P Wave ◽  
Model Parameters ◽  
Local Minima ◽  
S Wave ◽  
Traveltime Inversion ◽  
Velocity Models ◽  
Elastic Reflection

Elastic full-waveform inversion (FWI) updates high-resolution model parameters by minimizing the residuals of multicomponent seismic records between the field and model data. FWI suffers from the potential to converge to local minima and more serious nonlinearity than acoustic FWI mainly due to the absence of low frequencies in seismograms and the extended model domain (P- and S-velocities). Reflection waveform inversion can relax the nonlinearity by relying on the tomographic components, which can be used to update the low-wavenumber components of the model. Hence, we have developed an elastic reflection traveltime inversion (ERTI) approach to update the low-wavenumber component of the velocity models for the P- and S-waves. In our ERTI algorithm, we took the P- and S-wave impedance perturbations as elastic reflectivity to generate reflections and a weighted crosscorrelation as the misfit function. Moreover, considering the higher wavenumbers (lower velocity value) of the S-wave velocity compared with the P-wave case, optimizing the low-wavenumber components for the S-wave velocity is even more crucial in preventing the elastic FWI from converging to local minima. We have evaluated an equivalent decoupled velocity-stress wave equation to ERTI to reduce the coupling effects of different wave modes and to improve the inversion result of ERTI, especially for the S-wave velocity. The subsequent application on the Sigsbee2A model demonstrates that our ERTI method with the decoupled wave equation can efficiently update the low-wavenumber parts of the model and improve the precision of the S-wave velocity.

Processing‐induced anisotropy

Geophysics ◽  
2002 ◽  
Vol 67 (6) ◽  
pp. 1920-1928 ◽  
Author(s):  
Vladimir Grechka ◽  
Ilya Tsvankin
Keyword(s):  
Seismic Data ◽  
Symmetry Axis ◽  
Homogeneous Medium ◽  
Transversely Isotropic ◽  
P Wave ◽  
Well Logs ◽  
Velocity Analysis ◽  
Induced Anisotropy ◽  
Vertical Heterogeneity ◽  
Vti Media

Processing of seismic data is often performed under the assumption that the velocity distribution in the subsurface can be approximated by a macromodel composed of isotropic homogeneous layers or blocks. Despite being physically unrealistic, such models are believed to be sufficient for describing the kinematics of reflection arrivals. In this paper, we examine the distortions in normal‐moveout (NMO) velocities caused by the intralayer vertical heterogeneity unaccounted for in velocity analysis. To match P‐wave moveout measurements from a horizontal or a dipping reflector overlaid by a vertically heterogeneous isotropic medium, the effective homogeneous overburden has to be anisotropic. This apparent anisotropy is caused not only by velocity monotonically increasing with depth, but also by random velocity variations similar to those routinely observed in well logs. Assuming that the effective homogeneous medium is transversely isotropic with a vertical symmetry axis (VTI), we express the VTI parameters through the actual depth‐dependent isotropic velocity function. If the reflector is horizontal, combining the NMO and vertical velocities always results in nonnegative values of Thomsen's coefficient δ. For a dipping reflector, the inversion of the P‐wave NMO ellipse yields a nonnegative Alkhalifah‐Tsvankin coefficient η that increases with dip. The values of η obtained by two other methods (2‐D dip‐moveout inversion and nonhyperbolic moveout analysis) are also nonnegative but generally differ from that needed to fit the NMO ellipse. For truly anisotropic (VTI) media, the influence of vertical heterogeneity above the reflector can lead to a bias toward positive δ and η estimates in velocity analysis.

Fast and accurate dynamic raytracing in heterogeneous media

1990 ◽  
Vol 80 (5) ◽  
pp. 1284-1296
Author(s):  
Claude F. Lafond ◽  
Alan R. Levander
Keyword(s):  
Shooting Method ◽  
Heterogeneous Media ◽  
Velocity Model ◽  
Travel Times ◽  
Complex Velocity ◽  
Geometrical Spreading ◽  
Imaging Condition ◽  
Depth Migration ◽  
Velocity Models ◽  
Complex Velocity Model

Abstract We have developed a fast and accurate dynamic raytracing method for 2.5-D heterogeneous media based on the kinematic algorithm proposed by Langan et al. (1985). This algorithm divides the model into cells of constant slowness gradient, and the positions, directions, and travel times of the rays are expressed as polynomials of the travel path length, accurate to the second other in the gradient. This method is efficient because of the use of simple polynomials at each raytracing step. We derived similar polynomial expressions for the dynamic raytracing quantities by integrating the raytracing system and expanding the solutions to the second order in the gradient. This new algorithm efficiently computes the geometrical spreading, amplitude, and wavefront curvature on individual rays. The two-point raytracing problem is solved by the shooting method using the geometrical spreading. Paraxial corrections based on the wavefront curvature improve the accuracy of the travel time and amplitude at a given receiver. The computational results for two simple velocity models are compared with those obtained with the SEIS83 seismic modeling package (Cerveny and Psencik, 1984); this new method is accurate for both travel times and amplitudes while being significantly faster. We present a complex velocity model that shows that the algorithm allows for realistic models and easily computes rays in structures that pose difficulties for conventional methods. The method can be extended to raytracing in 3-D heterogeneous media and can be used as a support for a Gaussian beam algorithm. It is also suitable for computing the Green's function and imaging condition needed for prestack depth migration.

Migration velocity analysis from locally coherent events in 2‐D laterally heterogeneous media, Part II: Applications on synthetic and real data

Geophysics ◽  
2002 ◽  
Vol 67 (4) ◽  
pp. 1213-1224 ◽  
Author(s):  
Hervé Chauris ◽  
Mark S. Noble ◽  
Gilles Lambaré ◽  
Pascal Podvin
Keyword(s):  
Heterogeneous Media ◽  
Real Data ◽  
Migration Velocity ◽  
Velocity Analysis ◽  
Efficient Computation ◽  
Data Set ◽  
Seismic Reflection Data ◽  
Velocity Models ◽  
Migration Velocity Analysis ◽  
Paraxial Ray

We demonstrate a method for estimating 2‐D velocity models from synthetic and real seismic reflection data in the framework of migration velocity analysis (MVA). No assumption is required on the reflector geometry or on the unknown background velocity field, provided that the data only contain primary reflections/diffractions. In the prestack depth‐migrated volume, locations where the reflectivity exhibits local coherency are automatically picked without interpretation in two panels: common image gathers (CIGs) and common offset gathers (COGs). They are characterized by both their positions and two slopes. The velocity is estimated by minimizing all slopes picked in the CIGs. We test the applicability of the method on a real data set, showing the possibility of an efficient inversion using (1) the migration of selected CIGs and COGs, (2) automatic picking on prior uncorrelated locally coherent events, (3) efficient computation of the gradient of the cost function via paraxial ray tracing from the picked events to the surface, and (4) a gradient‐type optimization algorithm for convergence.

Sign in / Sign up

Export Citation Format

Share Document