Moveout analysis for tranversely isotropic media with a tilted symmetry axis

The transversely isotropic (TI) model with a tilted axis of symmetry may be typical, for instance, for sediments near the flanks of salt domes. This work is devoted to an analysis of reflection moveout from horizontal and dipping reflectors in the symmetry plane of TI media that contains the symmetry axis. While for vertical and horizontal transverse isotropy zero-offset reflections exist for the full range of dips up to 90°, this is no longer the case for intermediate axis orientations. For typical homogeneous models with a symmetry axis tilted towards the reflector, wavefront distortions make it impossible to generate specular zero-offset reflected rays from steep interfaces. The ‘missing’ dipping planes can be imaged only in vertically inhomogeneous media by using turning waves. These unusual phenomena may have serious implications in salt imaging. In non-elliptical TI media, the tilt of the symmetry axis may have a drastic influence on normal-moveout (NMO) velocity from horizontal reflectors, as well as on the dependence of NMO velocity on the ray parameter p (the ‘dip-moveout (DMO) signature’). The DMO signature retains the same character as for vertical transverse isotropy only for near-vertical and near-horizontal orientation of the symmetry axis. The behaviour of NMO velocity rapidly changes if the symmetry axis is tilted away from the vertical, with a tilt of ±20° being almost sufficient to eliminate the influence of the anisotropy on the DMO signature. For larger tilt angles and typical positive values of the difference between the anisotropic parameters ε and δ, the NMO velocity increases with p more slowly than in homogeneous isotropic media; a dependence usually caused by a vertical velocity gradient. Dip-moveout processing for a wide range of tilt angles requires application of anisotropic DMO algorithms. The strong influence of the tilt angle on P-wave moveout can be used to constrain the tilt using P-wave NMO velocity in the plane that includes the symmetry axis. However, if the azimuth of the axis is unknown, the inversion for the axis orientation cannot be performed without a 3D analysis of reflection traveltimes on lines with different azimuthal directions.     

3D description and inversion of reflection moveout of PS‐waves in anisotropic media

Common‐midpoint moveout of converted waves is generally asymmetric with respect to zero offset and cannot be described by the traveltime series t²(x²) conventionally used for pure modes. Here, we present concise parametric expressions for both common‐midpoint (CMP) and common‐conversion‐point (CCP) gathers of PS‐waves for arbitrary anisotropic, horizontally layered media above a plane dipping reflector. This analytic representation can be used to model 3D (multi‐azimuth) CMP gathers without time‐consuming two‐point ray tracing and to compute attributes of PS moveout such as the slope of the traveltime surface at zero offset and the coordinates of the moveout minimum. In addition to providing an efficient tool for forward modelling, our formalism helps to carry out joint inversion of P and PS data for transverse isotropy with a vertical symmetry axis (VTI media). If the medium above the reflector is laterally homogeneous, P‐wave reflection moveout cannot constrain the depth scale of the model needed for depth migration. Extending our previous results for a single VTI layer, we show that the interval vertical velocities of the P‐ and S‐waves (V_P0 and V_S0) and the Thomsen parameters ε and δ can be found from surface data alone by combining P‐wave moveout with the traveltimes of the converted PS(PSV)‐wave. If the data are acquired only on the dip line (i.e. in 2D), stable parameter estimation requires including the moveout of P‐ and PS‐waves from both a horizontal and a dipping interface. At the first stage of the velocity‐analysis procedure, we build an initial anisotropic model by applying a layer‐stripping algorithm to CMP moveout of P‐ and PS‐waves. To overcome the distorting influence of conversion‐point dispersal on CMP gathers, the interval VTI parameters are refined by collecting the PS data into CCP gathers and repeating the inversion. For 3D surveys with a sufficiently wide range of source–receiver azimuths, it is possible to estimate all four relevant parameters (V_P0, V_S0, ε and δ) using reflections from a single mildly dipping interface. In this case, the P‐wave NMO ellipse determined by 3D (azimuthal) velocity analysis is combined with azimuthally dependent traveltimes of the PS‐wave. On the whole, the joint inversion of P and PS data yields a VTI model suitable for depth migration of P‐waves, as well as processing (e.g. transformation to zero offset) of converted waves.     

PRACTICAL ASPECTS OF THE DETERMINATION OF 3-D STACKING VELOCITIES*

Proper stacking of three-dimensional seismic CDP-data generally requires the knowledge of normal moveout velocities in all source-receiver directions contributing to a CDP-gather. The azimuthal variation of the stacking velocities mainly depends on the dip of the seismic interfaces. For a single dipping plane a simple relation exists between the dip and the azimuthal variation of NMO-velocity. Varying strike and dip of subsequent reflectors, however, result in a complex dependency of the seismic parameters. Reliable information on the spatial distribution of the normal moveout (NMO)-velocity can be derived from a wavefront curvature estimation using a 3-D ray-tracing technique. These procedures require additional information, e.g. reflection time gradients or depth maps to show interval velocities between leading interfaces. Moreover, their application to an extended 3-D data volume is restricted by high costs. The need for a routine 3-D procedure resulted in a special data selection to create pseudo 2-D profiles and to apply existing velocity estimation routines to these profiles. At least three estimates in different directions are necessary to derive the full azimuthal velocity variation, characterized by the large and the small main axis and the orientation of the velocity ellipse. Errors are estimated by means of computer models. Stacking velocities obtained by mathematical routines (least-squares fit) and by seismic standard routines (NMO-correction and correlation) are compared. Finally, a general 3-D velocity procedure using cross-correlation of preliminarily NMO-corrected traces is proposed.     

What is DMO coverage?

'Coverage' or 'fold' is defined as the multiplicity of common-midpoint (CMP) data. For CMP stacking the coverage is consistent with the number of traces sharing a common reflection point on flat subsurface reflectors. This relationship is not true for dipping reflectors. The deficiencies of CMP stacking with respect to imaging dipping events have long been overcome by the introduction of the dip-moveout (DMO) correction. However, the concept of coverage has not yet satisfactorily been updated to a 'DMO coverage' consistent with DMO stacking. A definition of constant-velocity DMO coverage will be proposed here. A subsurface reflector will be illuminated from a given source and receiver location if the time difference between the reflector zero-offset traveltime and the NMO- and DMO-corrected traveltime of the reflection event is less than half a dominant wavelength. Due to the fact that a subsurface reflector location is determined by its zero-offset traveltime, its strike and its dip, the DMO coverage also depends on these three parameters. For every surface location, the proposed DMO coverage consists of a 3D fold distribution over reflector strike, dip and zero-offset traveltime.     

Seismic characterization of vertical fractures described as general linear-slip interfaces

Fluid flow in many hydrocarbon reservoirs is controlled by aligned fractures which make the medium anisotropic on the scale of seismic wavelength. Applying the linear‐slip theory, we investigate seismic signatures of the effective medium produced by a single set of ‘general’ vertical fractures embedded in a purely isotropic host rock. The generality of our fracture model means the allowance for coupling between the normal (to the fracture plane) stress and the tangential jump in displacement (and vice versa). Despite its low (triclinic) symmetry, the medium is described by just nine independent effective parameters and possesses several distinct features which help to identify the physical model and estimate the fracture compliances and background velocities. For example, the polarization vector of the vertically propagating fast shear wave S₁ and the semi‐major axis of the S₁‐wave normal‐moveout (NMO) ellipse from a horizontal reflector always point in the direction of the fracture strike. Moreover, for the S₁‐wave both the vertical velocity and the NMO velocity along the fractures are equal to the shear‐wave velocity in the host rock. Analysis of seismic signatures in the limit of small fracture weaknesses allows us to select the input data needed for unambiguous fracture characterization. The fracture and background parameters can be estimated using the NMO ellipses from horizontal reflectors and vertical velocities of P‐waves and two split S‐waves, combined with a portion of the P‐wave slowness surface reconstructed from multi‐azimuth walkaway vertical seismic profiling (VSP) data. The stability of the parameter‐estimation procedure is verified by performing non‐linear inversion based on the exact equations.     

Moveout approximation for horizontal transversely isotropic and vertical transversely isotropic layered medium. Part I: 1D ray propagation‡

Anisotropy in subsurface geological models is primarily caused by two factors: sedimentation in shale/sand layers and fractures. The sedimentation factor is mainly modelled by vertical transverse isotropy (VTI), whereas the fractures are modelled by a horizontal transversely isotropic medium (HTI). In this paper we study hyperbolic and non‐hyperbolic normal reflection moveout for a package of HTI/VTI layers, considering arbitrary azimuthal orientation of the symmetry axis at each HTI layer. We consider a local 1D medium, whose properties change vertically, with flat interfaces between the layers. In this case, the horizontal slowness is preserved; thus, the azimuth of the phase velocity is the same for all layers of the package. In general, however, the azimuth of the ray velocity differs from the azimuth of the phase velocity. The ray azimuth depends on the layer properties and may be different for each layer. In this case, the use of the Dix equation requires projection of the moveout velocity of each layer on the phase plane. We derive an accurate equation for hyperbolic and high‐order terms of the normal moveout, relating the traveltime to the surface offset, or alternatively, to the subsurface reflection angle. We relate the azimuth of the surface offset to its magnitude (or to the reflection angle), considering short and long offsets. We compare the derived approximations with analytical ray tracing.     

Interval anisotropic parameter estimation from walkaway vertical seismic profiling data

This paper presents a new explicit method for the estimation of layered vertical transverse isotropic (VTI) anisotropic parameters from walkaway VSP data. This method is based on Dix‐type normal moveout (NMO) inversion. To estimate interval anisotropic parameters above a receiver array, the method uses time arrivals of surface‐related double‐reflected downgoing waves. A three‐term NMO approximation function is used to estimate NMO velocity and a non‐hyperbolic parameter. Assuming the vertical velocity is known from zero‐offset VSP data, Dix‐type inversion is applied to estimate the layered Thomsen anisotropic parameters ?, δ above the receivers array. Model results show reasonable accuracy for estimates through Dix‐type inversion. Results also show that in many cases we can neglect the influence of the velocity gradient on anisotropy estimates. First breaks are used to estimate anisotropic parameters within the walkaway receiver interval. Analytical uncertainty analysis is performed to NMO parameter estimates. Its conclusions are confirmed by modelling.     

三维倾斜界面PS转换波CMP道集时距及参数估计

在PS转换波资料处理过程中,往往需要联合P波资料提供相应的模型.在实际应用中存在P波和PS转换波层位对比困难.本文仅利用PS转换波数据,通过三维倾斜界面PS转换波CMP道集精确时距关系推导了近似时距解析表达式;分析了PS波的精确与近似时距关系随测线方位、界面倾角与倾向的变化规律及其拟合误差;并讨论了近似时距关系的三个时距参数随方位的变化特征;理论上给出描述时距的三维倾斜界面倾角、倾向、深度、纵波速度和横波速度这5个独立参数的估计方法,并通过理论模拟数据证明了该方法的可行性.     

Non-hyperbolic moveout inversion of wide-azimuth P-wave data for orthorhombic media

The azimuthally varying non‐hyperbolic moveout of P‐waves in orthorhombic media can provide valuable information for characterization of fractured reservoirs and seismic processing. Here, we present a technique to invert long‐spread, wide‐azimuth P‐wave data for the orientation of the vertical symmetry planes and five key moveout parameters: the symmetry‐plane NMO velocities, V⁽¹⁾_nmo and V⁽²⁾_nmo , and the anellipticity parameters, η⁽¹⁾, η⁽²⁾ and η⁽³⁾ . The inversion algorithm is based on a coherence operator that computes the semblance for the full range of offsets and azimuths using a generalized version of the Alkhalifah–Tsvankin non‐hyperbolic moveout equation. The moveout equation provides a close approximation to the reflection traveltimes in layered anisotropic media with a uniform orientation of the vertical symmetry planes. Numerical tests on noise‐contaminated data for a single orthorhombic layer show that the best‐constrained parameters are the azimuth ? of one of the symmetry planes and the velocities V⁽¹⁾_nmo and V⁽²⁾_nmo , while the resolution in η⁽¹⁾ and η⁽²⁾ is somewhat compromised by the trade‐off between the quadratic and quartic moveout terms. The largest uncertainty is observed in the parameter η⁽³⁾ , which influences only long‐spread moveout in off‐symmetry directions. For stratified orthorhombic models with depth‐dependent symmetry‐plane azimuths, the moveout equation has to be modified by allowing the orientation of the effective NMO ellipse to differ from the principal azimuthal direction of the effective quartic moveout term. The algorithm was successfully tested on wide‐azimuth P‐wave reflections recorded at the Weyburn Field in Canada. Taking azimuthal anisotropy into account increased the semblance values for most long‐offset reflection events in the overburden, which indicates that fracturing is not limited to the reservoir level. The inverted symmetry‐plane directions are close to the azimuths of the off‐trend fracture sets determined from borehole data and shear‐wave splitting analysis. The effective moveout parameters estimated by our algorithm provide input for P‐wave time imaging and geometrical‐spreading correction in layered orthorhombic media.     

Migration to multiple offset and velocity analysis1

The stacking velocity best characterizes the normal moveout curves in a common-mid-point gather, while the migration velocity characterizes the diffraction curves in a zero-offset section as well as in a common-midpoint gather. For horizontally layered media, the two velocity types coincide due to the conformance of the normal and the image ray. In the case of dipping subsurface structures, stacking velocities depend on the dip of the reflector and relate to normal rays, but with a dip-dependent lateral smear of the reflection point. After dip-moveout correction, the stacking velocities are reduced while the reflection-point smear vanishes, focusing the rays on the common reflection points. For homogeneous media the dip-moveout correction is independent of the actual velocity and can be applied as a dip-moveout correction to multiple offset before velocity analysis. Migration to multiple offset is a prestack, time-migration technique, which presents data sets which mimic high-fold, bin-centre adjusted, common-midpoint gathers. This method is independent of velocity and can migrate any 2D or 3D data set with arbitrary acquisition geometry. The gathers generated can be analysed for normal-moveout velocities using traditional methods such as the interpretation of multivelocity-function stacks. These stacks, however, are equivalent to multi-velocity-function time migrations and the derived velocities are migration velocities.     

A SIMPLE EFFICIENT METHOD OF DIP-MOVEOUT CORRECTION1

Much of the success of modern seismic data processing derives from the use of the stacking process. Unfortunately, as is well known, conventional normal moveout correction (NMO) introduces mispositioning of data, and hence mis-stacking, when dip is present. Dip moveout correction (DMO) is a technique that converts non-zero-offset seismic data after NMO to true zero-offset locations and reflection times, irrespective of dip. The combination of NMO and DMO followed by post-stack time migration is equivalent to, but can be implemented much more efficiently than, full time migration before stack. In this paper we consider the frequency-wavenumber DMO algorithm developed by Hale. Our analysis centres on the result that, for a given dip, the combination of NMO at migration velocity and DMO is equivalent to NMO at the appropriate, dip-dependent, stacking velocity. This perspective on DMO leads to computationally efficient methods for applying Hale DMO and also provides interesting insights on the nature of both DMO and conventional stacking.     

Characterization of dipping fractures in a transverselyisotropic background

Although it is believed that natural fracture sets predominantly have near‐vertical orientation, oblique stresses and some other mechanisms may tilt fractures away from the vertical. Here, we examine an effective medium produced by a single system of obliquely dipping rotationally invariant fractures embedded in a transversely isotropic with a vertical symmetry axis (VTI) background rock. This model is monoclinic with a vertical symmetry plane that coincides with the dip plane of the fractures. Multicomponent seismic data acquired over such a medium possess several distinct features that make it possible to estimate the fracture orientation. For example, the vertically propagating fast shear wave (and the fast converted PS‐wave) is typically polarized in the direction of the fracture strike. The normal‐moveout (NMO) ellipses of horizontal reflection events are co‐orientated with the dip and strike directions of the fractures, which provides an independent estimate of the fracture azimuth. However, the polarization vector of the slow shear wave at vertical incidence does not lie in the horizontal plane – an unusual phenomenon that can be used to evaluate fracture dip. Also, for oblique fractures the shear‐wave splitting coefficient at vertical incidence becomes dependent on fracture infill (saturation). A complete medium‐characterization procedure includes estimating the fracture compliances and orientation (dip and azimuth), as well as the Thomsen parameters of the VTI background. We demonstrate that both the fracture and background parameters can be obtained from multicomponent wide‐azimuth data using the vertical velocities and NMO ellipses of PP‐waves and two split SS‐waves (or the traveltimes of PS‐waves) reflected from horizontal interfaces. Numerical tests corroborate the accuracy and stability of the inversion algorithm based on the exact expressions for the vertical and NMO velocities.     

P‐wave stacking‐velocity tomography for VTI media

A major complication caused by anisotropy in velocity analysis and imaging is the uncertainty in estimating the vertical velocity and depth scale of the model from surface data. For laterally homogeneous VTI (transversely isotropic with a vertical symmetry axis) media above the target reflector, P‐wave moveout has to be combined with other information (e.g. borehole data or converted waves) to build velocity models for depth imaging. The presence of lateral heterogeneity in the overburden creates the dependence of P‐wave reflection data on all three relevant parameters (the vertical velocity V_P0 and the Thomsen coefficients ε and δ) and, therefore, may help to determine the depth scale of the velocity field. Here, we propose a tomographic algorithm designed to invert NMO ellipses (obtained from azimuthally varying stacking velocities) and zero‐offset traveltimes of P‐waves for the parameters of homogeneous VTI layers separated by either plane dipping or curved interfaces. For plane non‐intersecting layer boundaries, the interval parameters cannot be recovered from P‐wave moveout in a unique way. Nonetheless, if the reflectors have sufficiently different azimuths, a priori knowledge of any single interval parameter makes it possible to reconstruct the whole model in depth. For example, the parameter estimation becomes unique if the subsurface layer is known to be isotropic. In the case of 2D inversion on the dip line of co‐orientated reflectors, it is necessary to specify one parameter (e.g. the vertical velocity) per layer. Despite the higher complexity of models with curved interfaces, the increased angle coverage of reflected rays helps to resolve the trade‐offs between the medium parameters. Singular value decomposition (SVD) shows that in the presence of sufficient interface curvature all parameters needed for anisotropic depth processing can be obtained solely from conventional‐spread P‐wave moveout. By performing tests on noise‐contaminated data we demonstrate that the tomographic inversion procedure reconstructs both the interfaces and the VTI parameters with high accuracy. Both SVD analysis and moveout inversion are implemented using an efficient modelling technique based on the theory of NMO‐velocity surfaces generalized for wave propagation through curved interfaces.     

双参数展开CRP叠加和速度分析方法研究

椭圆展开共反射点(CRP)方法可以获得比常规倾角时差校正(DMO)方法更近似的零偏移距时间剖面和相应CRP速度场.大量研究和实践证实,在非均质性较弱的地区,该方法取得的成果显著.但由于该方法没有考虑速度的横向变化和转换波等情况,当地下介质存在较强非均质性时,该方法不再准确,需要引进反映速度横向变化的双参数(上行波与下行波的平均速度和速度比)进行改进.本文详细推导了引入双参数后的叠加和速度分析算法,并通过数值模型和地震资料处理证实,修正后的算法可以更好地解决地质复杂地区速度建模和叠加成像问题.     

Anisotropic velocity analysis1

Results from walkaway VSP and shale laboratory experiments show that shale anisotropy can be significantly anelliptic. Heterogeneity and anellipticity both lead to non-hyperbolic moveout curves and the resulting ambiguity in velocity analysis is investigated for the case of a factorizable anisotropic medium with a linear dependence of velocity on depth. More information can be obtained if there are several reflectors. The method of Dellinger et al. for anisotropic velocity analysis in layered transversely isotropic media is examined and is shown to be restricted to media having relatively small anellipticity. A new scheme, based on an expansion of the inverse-squared group velocity in spherical harmonics, is presented. This scheme can be used for larger anellipticity, and is applicable for horizontal layers having monoclinic symmetry with the symmetry plane parallel to the layers. The method is applied to invert the results of anisotropic ray tracing on a model Sand/shale sequence. For transversely isotropic media with small anisotropy, the scheme reduces to the method of Byun et al. and Byun and Corrigan. The expansion in spherical harmonics allows the P-phase slowness surface of each layer to be determined in analytic form from the layer parameters obtained by inversion without the need to assume that the anisotropy is weak.     

Moveout approximation for vertical seismic profile geometry in a 2D model with anisotropic layers

I introduce a new explicit form of vertical seismic profile (VSP) traveltime approximation for a 2D model with non‐horizontal boundaries and anisotropic layers. The goal of the new approximation is to dramatically decrease the cost of time calculations by reducing the number of calculated rays in a complex multi‐layered anisotropic model for VSP walkaway data with many sources. This traveltime approximation extends the generalized moveout approximation proposed by Fomel and Stovas. The new equation is designed for borehole seismic geometry where the receivers are placed in a well while the sources are on the surface. For this, the time‐offset function is presented as a sum of odd and even functions. Coefficients in this approximation are determined by calculating the traveltime and its first‐ and second‐order derivatives at five specific rays. Once these coefficients are determined, the traveltimes at other rays are calculated by this approximation. Testing this new approximation on a 2D anisotropic model with dipping boundaries shows its very high accuracy for offsets three times the reflector depths. The new approximation can be used for 2D anisotropic models with tilted symmetry axes for practical VSP geometry calculations. The new explicit approximation eliminates the need of massive ray tracing in a complicated velocity model for multi‐source VSP surveys. This method is designed not for NMO correction but for replacing conventional ray tracing for time calculations.     

Synthetic Experiments Dealing with Anisotropic Model Building

It has been shown in the past that the interval-NMO velocity and the non-ellipticity parameter largely control the P-wave reflection time moveout of VTI media. To invert for these two parameters, one needs either reasonably large offsets, or some structure in the subsurface in combination with relatively mild lateral velocity variation.This paper deals with a simulation of an inversion approach, building on the assumption that accurately measured V _NMO, as defined by small offset asymptotics for a particular reflector, were available. Instead of such measurements we take synthetically computed data. First, an isotropic model is constructed which explains these V _NMO. Subsequently, residual moveout in common image gathers is modelled by ray tracing (replacing real data), along with its sensitivity for changes in the interval-NMO velocity and the non-ellipticity parameter under the constraint that V _NMO is preserved. This enables iterative updating of the non-ellipticity parameter and the interval-NMO velocity in a layer that can be laterally inhomogeneous.This approach is successfully applied for a mildly dipping reflector at the bottom of a layer with laterally varying medium parameters. With the exact V _NMO assumed to be given, lateral inhomogeneity and anisotropy can be distinguished for such a situation. However, for another example with a homogeneous VTI layer overlying a curved reflector with dip up to 30°, there appears to be an ambiguity which can be understood by theoretical analysis. Consistently with existing theory using the NMO-ellipse, the presented approach is successfully applied to the latter example if V _NMO in the strike direction is combined with residual moveout in dip direction.     

Reflection points and surface points1

A simple expression ties the midpoint of a surface spread to reflection points on a dipping plane. If we use two coordinate systems, an unprimed one with a z-axis perpendicular to the surface and a primed one with a z-axis perpendicular to the reflector, we have where θ is the dip angle, φ is the profile angle, X is the source-to-receiver separation, and D is the depth of the reflector. The reflection point is (x, y_p, D) and the surface midpoint is (x_c, y_c, 0). Using the expression, I show that if complete azimuthal coverage is required at a CMP position, the reflection points lie on an ellipse. Similarly, a fixed reflection point generates a circle of surface midpoints. A circle of CMP positions for fixed θ and φ becomes an ellipse of reflection points and a circle of reflection points becomes an ellipse of midpoints. A user can easily find the shape and location of the reflection area generated by a surface aperture.     

Reflection moveout inversion in azimuthally anisotropic media: accuracy, limitations and acquisition

Parameter estimation from the elliptical variations in the normal-moveout (NMO) velocity in azimuthally anisotropic media is sensitive to the angular separation between the survey lines in 2D, or equivalently, the source-to-receiver azimuth in 3D, and to the set of azimuths used in the inversion procedure. The accuracy in estimating the orientation of an NMO ellipse, in particular the parameter α, is also sensitive to the magnitude of anisotropy. On the other hand, the accuracy in estimating the semi-axes of the NMO-velocity ellipse is about the same for any magnitude of anisotropy.   To invert for the NMO ellipse parameters at least three NMO-velocity measurements along distinct azimuth directions are needed. In order to maximize the accuracy and stability in parameter estimation, it is best to have the azimuths for the three source-to-receiver directions 60° apart. Having more than three distinct source-to-receiver azimuths (e.g. full azimuthal coverage) provides a useful data redundancy that enhances the quality of the estimates.   In order to maximize quality in the inversion process, it is recommended to design the seismic data acquisition such that it contains small sectors (≤10°) with adequate fold and offset distribution.   Using three NMO-velocity measurements, 60° apart, an azimuthally anisotropic layer overlain by an azimuthally isotropic overburden (as might occur for fractured reservoirs) should have a relative thickness (in time) with respect to the total thickness at least equal to the ratio of the error in the NMO (stacking) velocity to the interval anisotropy of the fractured layer. Coverage along more than three azimuths, however, improves this limitation, which is imposed by Dix differentiation, by at most 50%, depending on the number of observations (NMO velocities) that enter the inversion procedure.