3‐D moveout inversion in azimuthally anisotropic media with lateral velocity variation: Theory and a case study

1998 ◽  
Author(s):  
Vladimir Grechka ◽  
Ilya Tsvankin
Keyword(s):  
Anisotropic Media ◽  
Velocity Variation ◽  
Variation Theory ◽  
Lateral Velocity

3-D moveout inversion in azimuthally anisotropic media with lateral velocity variation: Theory and a case study

Geophysics ◽  
1999 ◽  
Vol 64 (4) ◽  
pp. 1202-1218 ◽  
Author(s):  
Vladimir Grechka ◽  
Ilya Tsvankin
Keyword(s):  
Correction Term ◽  
Anisotropic Media ◽  
Azimuthal Anisotropy ◽  
P Wave ◽  
Velocity Variation ◽  
Variation Theory ◽  
Fractured Reservoir ◽  
Lateral Velocity ◽  
Data Set ◽  
Zero Offset

Reflection moveout recorded over an azimuthally anisotropic medium (e.g., caused by vertical or dipping fractures) varies with the azimuth of the source‐receiver line. Normal‐moveout (NMO) velocity, responsible for the reflection traveltimes on conventional‐length spreads, forms an elliptical curve in the horizontal plane. While this result remains valid in the presence of arbitrary anisotropy and heterogeneity, the inversion of the NMO ellipse for the medium parameters has been discussed so far only for horizontally homogeneous models above a horizontal or dipping reflector. Here, we develop an analytic moveout correction for weak lateral velocity variation in horizontally layered azimuthally anisotropic media. The correction term is proportional to the curvature of the zero‐offset traveltime surface at the common midpoint and, therefore, can be estimated from surface seismic data. After the influence of lateral velocity variation on the effective NMO ellipses has been stripped, the generalized Dix equation can be used to compute the interval ellipses and evaluate the magnitude of azimuthal anisotropy (measured by P-wave NMO velocity) within the layer of interest. This methodology was applied to a 3-D “wide‐azimuth” data set acquired over a fractured reservoir in the Powder River Basin, Wyoming. The processing sequence included 3-D semblance analysis (based on the elliptical NMO equation) for a grid of common‐midpoint “supergathers,” spatial smoothing of the effective NMO ellipses and zero‐offset traveltimes, correction for lateral velocity variation, and generalized Dix differentiation. Our estimates of depth‐varying fracture trends in the survey area, based on the interval P-wave NMO ellipses, are in good agreement with the results of outcrop and borehole measurements and the rotational analysis of four‐ component S-wave data.

One-way true-amplitude migration using matrix decomposition

Geophysics ◽  
2018 ◽  
Vol 83 (5) ◽  
pp. S387-S398 ◽  
Author(s):  
Jiachun You ◽  
Ru-Shan Wu ◽  
Xuewei Liu
Keyword(s):  
Wave Equations ◽  
Matrix Decomposition ◽  
Dynamic Imaging ◽  
Velocity Variation ◽  
Lateral Velocity ◽  
Amplitude Wave ◽  
Contrast Model ◽  
Migration Method ◽  
The One ◽  
Imaging Performance

To meet the requirement of true-amplitude migration and address the shortcomings of the classic one-way wave equations on the dynamic imaging, one-way true-amplitude wave equations were developed. Migration methods, based on the Taylor or other series approximation theory, are introduced to solve the one-way true-amplitude wave equations. This leads to the main weakness of one-way true-amplitude migration for imaging the complex or strong velocity — contrast media — the limited imaging angles. To deal with this issue, we apply a matrix decomposition method to accurately calculate the square-root operator and impose the boundary conditions of the one-way true-amplitude wave equations. Our migration method and the conventional one-way true-amplitude Fourier finite-difference (FFD) migration method are used by us to test and compare the imaging performance. The impulse responses in a strong velocity-contrast model prove that our migration method works for larger imaging angles than the one-way true-amplitude FFD method. The amplitude calculations in a strong-lateral velocity variation media with one reflector and in the Marmousi model demonstrate that our migration method provides better amplitude-preserving performance and offers higher structural imaging quality than the one-way true-amplitude FFD method. We also use field data to indicate the imaging enhancement and the feasibility of our method compared with the one-way true-amplitude FFD method. Our one-way true-amplitude migration method using matrix decomposition fully exploits the features of one-way true-amplitude wave equations with less approximation, and it is capable of producing more accurate amplitude estimations and potentially wider imaging angles.

Velocity estimation in laterally varying media

Geophysics ◽  
1982 ◽  
Vol 47 (6) ◽  
pp. 884-897 ◽  
Author(s):  
Walter S. Lynn ◽  
Jon F. Claerbout
Keyword(s):  
Taylor Series Expansion ◽  
Velocity Estimation ◽  
Lateral Resolution ◽  
Velocity Variation ◽  
Lateral Velocity ◽  
Depth Migration ◽  
Derivative Term ◽  
Fourth Derivative ◽  
Zero Offset ◽  
Lateral Variations

In areas of large lateral variations in velocity, stacking velocities computed on the basis of hyperbolic moveout can differ substantially from the actual root mean square (rms) velocities. This paper addresses the problem of obtaining rms or migration velocities from stacking velocities in such areas. The first‐order difference between the stacking and the vertical rms velocities due to lateral variations in velocity are shown to be related to the second lateral derivative of the rms slowness [Formula: see text]. Approximations leading to this relation are straight raypaths and that the vertical rms slowness to a given interface can be expressed as a second‐order Taylor series expansion in the midpoint direction. Under these approximations, the effect of the first lateral derivative of the slowness on the traveltime is negligible. The linearization of the equation relating the stacking and true velocities results in a set of equations whose inversion is unstable. Stability is achieved, however, by adding a nonphysical fourth derivative term which affects only the higher spatial wavenumbers, those beyond the lateral resolution of the lateral derivative method (LDM). Thus, given the stacking velocities and the zero‐offset traveltime to a given event as a function of midpoint, the LDM provides an estimate of the true vertical rms velocity to that event with a lateral resolution of about two mute zones or cable lengths. The LDM is applicable when lateral variations of velocity greater than 2 percent occur over the mute zone. At variations of 30 percent or greater, the internal assumptions of the LDM begin to break down. Synthetic models designed to test the LDM when the different assumptions are violated show that, in all cases, the results are not seriously affected. A test of the LDM on field data having a lateral velocity variation caused by sea floor topography gives a result which is supported by depth migration.

Migration velocity analysis for TI media in the presence of quadratic lateral velocity variation

Geophysics ◽  
2012 ◽  
Vol 77 (6) ◽  
pp. U87-U96 ◽  
Author(s):  
Mamoru Takanashi ◽  
Ilya Tsvankin
Keyword(s):  
Migration Velocity ◽  
P Wave ◽  
Velocity Variation ◽  
Small Scale ◽  
Velocity Analysis ◽  
Lateral Heterogeneity ◽  
Lateral Velocity ◽  
Symmetry Direction ◽  
Migration Velocity Analysis ◽  
Ti Media

One of the most serious problems in anisotropic velocity analysis is the trade-off between anisotropy and lateral heterogeneity, especially if velocity varies on a scale smaller than the maximum offset. We have developed a P-wave MVA (migration velocity analysis) algorithm for transversely isotropic (TI) models that include layers with small-scale lateral heterogeneity. Each layer is described by constant Thomsen parameters [Formula: see text] and [Formula: see text] and the symmetry-direction velocity [Formula: see text] that varies as a quadratic function of the distance along the layer boundaries. For tilted TI media (TTI), the symmetry axis is taken orthogonal to the reflectors. We analyzed the influence of lateral heterogeneity on image gathers obtained after prestack depth migration and found that quadratic lateral velocity variation in the overburden can significantly distort the moveout of the target reflection. Consequently, medium parameters beneath the heterogeneous layer(s) are estimated with substantial error, even when borehole information (e.g., check shots or sonic logs) is available. Because residual moveout in the image gathers is highly sensitive to lateral heterogeneity in the overburden, our algorithm simultaneously inverts for the interval parameters of all layers. Synthetic tests for models with a gently dipping overburden demonstrate that if the vertical profile of the symmetry-direction velocity [Formula: see text] is known at one location, the algorithm can reconstruct the other relevant parameters of TI models. The proposed approach helps increase the robustness of anisotropic velocity model-building and enhance image quality in the presence of small-scale lateral heterogeneity in the overburden.

Plumes: Response of time migration to lateral velocity variation

Geophysics ◽  
1995 ◽  
Vol 60 (4) ◽  
pp. 1118-1127 ◽  
Author(s):  
Dimitri Bevc ◽  
James L. Black ◽  
Gopal Palacharla
Keyword(s):  
Impulse Response ◽  
Synthetic Data ◽  
Migration Velocity ◽  
Velocity Variation ◽  
Velocity Dependence ◽  
Geometrical Construction ◽  
Lateral Velocity ◽  
Time Migration ◽  
Offset Time ◽  
Zero Offset

We analyze how time migration mispositions events in the presence of lateral velocity variation by examining the impulse response of depth modeling followed by time migration. By examining this impulse response, we lay the groundwork for the development of a remedial migration operator that links time and depth migration. A simple theory by Black and Brzostowski predicted that the response of zero‐offset time migration to a point diffractor in a v(x, z) medium would be a distinctive, cusp‐shaped curve called a plume. We have constructed these plumes by migrating synthetic data using several time‐migration methods. We have also computed the shape of the plumes by two geometrical construction methods. These two geometrical methods compare well and explain the observed migration results. The plume response is strongly influenced by migration velocity. We have studied this dependency by migrating synthetic data with different velocities. The observed velocity dependence is confirmed by geometrical construction. A simple first‐order theory qualitatively explains the behavior of zero‐offset time migration, but a more complete understanding of migration velocity dependence in a v(x, z) medium requires a higher order finite‐offset theory.

Enhancements to prestack frequency‐wavenumber (f-k) migration

Geophysics ◽  
1991 ◽  
Vol 56 (1) ◽  
pp. 27-40 ◽  
Author(s):  
Z. Li ◽  
W. Lynn ◽  
R. Chambers ◽  
Ken Larner ◽  
Ray Abma
Keyword(s):  
Computational Effort ◽  
Velocity Variation ◽  
Lateral Velocity ◽  
Time Migration ◽  
Migration Velocity Analysis ◽  
Prestack Time Migration ◽  
And Migration ◽  
Definition Of ◽  
Ray Bending ◽  
Excellent Source

Prestack frequency‐wavenumber (f-k) migration is a particularly efficient method of doing both full prestack time migration and migration velocity analysis. Conventional implementations of the method, however, can encounter several drawbacks: (1) poor resolution and spatial aliasing noise caused by insufficient sampling in the offset dimension, (2) poor definition of steep events caused by insufficient sampling in the velocity dimension, and (3) inadequate handling of ray bending for steep events. All three of these problems can be mitigated with modifications to the prestack f-k algorithm. The application of linear moveout (LMO) in the offset dimension prior to migration reduces event moveout and hence increases the bandwidth of non‐spatially aliased signals. To reduce problems of interpolation for steep events, the number of constant‐velocity migrations can be economically increased by performing residual poststack migrations. Finally, migration with a dip‐dependent imaging velocity addresses the issue of ray bending and thereby improves the positioning of steep events. None of these enhancements substantially increases the computational effort of f-k migration. Prestack f-k migration possesses a limitation for which no solution is readily available. Where lateral velocity variation is modest, steep events (such as fault‐plane reflections in sediments) may not be imaged as well as by other migration approaches. This shortcoming results from the restriction that, in the prestack f-k approach, a single velocity field must serve to perform two different functions: imaging and stacking. Nevertheless, in areas of strong velocity variation and gentle to moderate dip, the detailed velocity control afforded by the prestack f-k method is an excellent source of geologic information.

Sign in / Sign up

Export Citation Format

Share Document