Anisotropic attenuation of stratified viscoelastic media

An important cause of seismic anisotropic attenuation is the interbedding of thin viscoelastic layers. However, much less attention has been devoted to layer‐induced anisotropic attenuation. Here, we derive a group of unified weighted average forms for effective attenuation from a binary isotropic, transversely isotropic‐ and orthorhombic‐layered medium in the zero‐frequency limit by using the Backus averaging/upscaling method and analyse the influence of interval parameters on effective attenuation. Besides the corresponding interval attenuation and the real part of stiffness, the contrast in the real part of the complex stiffness is also a key factor influencing effective attenuation. A simple linear approximation can be obtained to calculate effective attenuation if the contrast in the real part of stiffness is very small. In a viscoelastic medium, attenuation anisotropy and velocity anisotropy may have different orientations of symmetry planes, and the symmetry class of the former is not lower than that of the latter. We define a group of more general attenuation‐anisotropy parameters to characterize not only the anisotropic attenuation with different symmetry classes from the anisotropic velocity but also the elastic case. Numerical tests reveal the influence of interval attenuation anisotropy, interval velocity anisotropy and the contrast in the real part of stiffness on effective attenuation anisotropy. Types of effective attenuation anisotropy for interval orthorhombic attenuation and interval transversely isotropic attenuation with a vertical symmetry (vertical transversely isotropic attenuation) are controlled only by the interval attenuation anisotropy. A type of effective attenuation anisotropy for interval TI attenuation with a horizontal symmetry (horizontal transversely isotropic attenuation) is controlled by the interval attenuation anisotropy and the contrast in the real part of stiffness. The type of effective attenuation anisotropy for interval isotropic attenuation is controlled by all three factors. The magnitude of effective attenuation anisotropy is positively correlated with the contrast in the real part of the stiffness. Effective attenuation even in isotropic layers with identical isotropic attenuation is anisotropic if the contrast in the real part of stiffness is non‐zero. In addition, if the contrast in the real part of stiffness is very small, a simple linear approximation also can be performed to calculate effective attenuation‐anisotropy parameters for interval anisotropic attenuation.     

Effective ellipsoidal models for wavefield extrapolation in tilted orthorhombic media

Wavefield computations using the ellipsoidally anisotropic extrapolation operator offer significant cost reduction compared to that for the orthorhombic case, especially when the symmetry planes are tilted and/or rotated. However, ellipsoidal anisotropy does not provide accurate wavefield representation or imaging for media of orthorhombic symmetry. Therefore, we propose the use of ‘effective ellipsoidally anisotropic’ models that correctly capture the kinematic behaviour of wavefields for tilted orthorhombic (TOR) media. We compute effective velocities for the ellipsoidally anisotropic medium using kinematic high-frequency representation of the TOR wavefield, obtained by solving the TOR eikonal equation. The effective model allows us to use the cheaper ellipsoidally anisotropic wave extrapolation operators. Although the effective models are obtained by kinematic matching using high-frequency asymptotic, the resulting wavefield contains most of the critical wavefield components, including frequency dependency and caustics, if present, with reasonable accuracy. The proposed methodology offers a much better cost versus accuracy trade-off for wavefield computations in TOR media, particularly for media of low to moderate anisotropic strength. Furthermore, the computed wavefield solution is free from shear-wave artefacts as opposed to the conventional finite-difference based TOR wave extrapolation scheme. We demonstrate applicability and usefulness of our formulation through numerical tests on synthetic TOR models.     

Kinematic parameters of pure‐ and converted‐mode waves in elastic orthorhombic media

I derive the kinematic properties of single‐mode P, S1, and S2 waves as well as converted PS1, PS2, and S1S2 waves in elastic orthorhombic media including vertical velocity, two normal moveout velocities defined in vertical symmetry planes, and three anelliptic parameters (two of them are defined in vertical symmetry plane and one parameter is the cross‐term one). I show that the azimuthal dependence of normal moveout velocity and anellipticity is different in phase and group domains. The effects on‐vertical‐axis singularity and on‐vertical‐axis triplication are considered for pure‐mode S1 and S2 waves and converted‐mode S1S2 waves. The conditions and properties of on‐vertical‐axis triplication are defined in terms of kinematic parameters. The results are illustrated in four homogeneous orthorhombic models and one multilayered orthorhombic model with no variation in azimuthal orientation for all the layers.     

Low-frequency layer-induced dispersion in a weak-contrast vertically heterogeneous orthorhombic medium

In this paper, we consider wave propagation in a periodically layered medium with orthorhombic symmetry. The weak-contrast approximation is utilized to derive the low-frequency dispersion in effective properties for P, S1 and S2 waves. We show that the dispersion term for all effective properties is controlled by the second-order contrasts in elastic properties from the layers. We also compute the sensitivity matrices for second- and fourth-order coefficients from eigenvalues of frequency-dependent system matrix associated with kinematic parameters for individual wave modes.     

Validity of the long-wave approximation in periodically layered media

In seismic modelling, a stack of thin layers is often replaced by an effective equivalent anisotropic homogeneous slab. For waves with finite wavelength, this is an approximation, and the error thus introduced can be quantified by considering the relative error in the phase velocity between the layer stack and the effective medium. For periodic layering, the relative phase-velocity error can be expressed in closed form as a function of wavelength, reflection coefficients and layer thicknesses. By comparing the relative phase-velocity error with laboratory measurements and numerical simulations, we find that the difference in seismic response between a periodic layer stack and an equivalent effective medium depends not only on wavelength, but it also depends significantly on reflection coefficients and the ratio between layer thicknesses. For a 1% relative error in the phase velocity, and if all layers have the same thickness measured in vertical traveltime, we find that the wavelength must be larger than approximately three times the layer period for a reflection coefficient of 0.1, but this increases to 13 times the layer period for a reflection coefficient of 0.9, which is highly unrealistic in a geological setting.     

OnP F and other phases in the early part of seismogram

Summary It has been found that when seismic energy propagates along the surface of the homogeneous crust beside usual Rayleigh waves, it produces certain instability in layers through which it propagates. In the light of this instability, a type of motion corresponding to longitudinal wave will be prominent in horizontal component compared to the vertical component; while transverse wave will be prominent in the vertical component but weak in the horizontal component, a contradiction with the existing knowledge. This has been identified withP _F phase. On taking the medium of propagation as slightly heterogeneous which allows existence of low velocity layer, a few larger number of such instabilities have been found. Velocity equation for Rayleigh waves for such media reveals existence of different velocities corresponding to vertical and horizontal components. Table for these velocities has been furnished.     

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.     

Reflection and transmission responses of layered transversely isotropic viscoelastic media

We consider a layered heterogeneous viscoelastic transversely isotropic medium with a vertical symmetry axis (a viscoelastic TIV medium) and parameters that depend on depth only. This takes into account intrinsic attenuation, anisotropy and thin layering. The seismic wavefield is decomposed into up- and downgoing waves scaled by the vertical energy flux. This scaling gives important symmetry relationships for both reflection and transmission (R/T) responses. For a stack of homogeneous layers, the exact reflection response can be computed in a numerically stable way by a simple layer-recursive algorithm. We derive exact plane-wave R/T coefficients and several linear and quadratic approximations between two viscoelastic TIV media, as functions of the real-valued horizontal slowness. The approximations are valid for pre- and post-critical values of horizontal slowness provided that the proper complex square roots are used when computing the vertical slowness. Numerical examples demonstrate that the quadratic approximations can be used for large differences in medium parameters, while the linear approximations can be used for small differences. For weak anisotropy it is sufficient to use an isotropic background medium, while for strong anisotropy it is necessary to use a weak TIV or TIV background medium. We also extend the O'Doherty–Anstey formula to the P- and SV-wave transmission responses of a stack of viscoelastic TIV layers, taking into account intrinsic attenuation, anisotropy and thin layering.     

Preserved travel‐time smoothing in orthorhombic media

Certain degree of smoothness of velocity models is required for most ray‐based migration and tomography. Applying conventional smoothing in model parameters results in offset‐dependent travel‐time errors for reflected events, which can be large even for small contrasts in model parameters between the layers. This causes the shift in both the depth and residual moveout of the migrated images. To overcome this problem in transversely isotropic medium with a vertical symmetry axis, the preserved travel‐time smoothing method was proposed earlier. We extend this method for orthorhombic media with and without azimuthal variation between the layers. We illustrate this method for a single interface between two orthorhombic layers and show that the smoothing‐driven errors in travel time are very small for practical application.     

Out‐of‐plane geometrical spreading in anisotropic media

Two-dimensional seismic processing is successful in media with little structural and velocity variation in the direction perpendicular to the plane defined by the acquisition direction and the vertical axis. If the subsurface is anisotropic, an additional limitation is that this plane is a plane of symmetry. Kinematic ray propagation can be considered as a two-dimensional process in this type of medium. However, two-dimensional processing in a true-amplitude sense requires out-of-plane amplitude corrections in addition to compensation for in-plane amplitude variation. We provide formulae for the out-of-plane geometrical spreading for P- and S-waves in transversely isotropic and orthorhombic media. These are extensions of well-known isotropic formulae.
For isotropic and transversely isotropic media, the ray propagation is independent of the azimuthal angle. The azimuthal direction is defined with respect to a possibly tilted axis of symmetry. The out-of-plane spreading correction can then be calculated by integrating quantities which describe in-plane kinematics along in-plane rays. If, in addition, the medium varies only along the vertical direction and has a vertical axis of symmetry, no ray tracing need be carried out. All quantities affecting the out-of-plane geometrical spreading can be derived from traveltime information available at the observation surface.
Orthorhombic media possess no rotational symmetry and the out-of-plane geometrical spreading includes parameters which, even in principle, are not invertible from in-plane experiments. The exact and approximate formulae derived for P- and S-waves are nevertheless useful for modelling purposes.     

Perturbation‐based moveout approximations in anisotropic media

The moveout approximations play an important role in seismic data processing. The standard hyperbolic moveout approximation is based on an elliptical background model with two velocities: vertical and normal moveout. We propose a new set of moveout approximations based on a perturbation series in terms of anellipticity parameters using the alternative elliptical background model defined by vertical and horizontal velocities. We start with a transversely isotropic medium with a vertical symmetry axis. Then, we extend this approach to a homogeneous orthorhombic medium. To define the perturbation coefficients for a new background, we solve the eikonal equation with horizontal velocities in transversely isotropic medium with a vertical symmetry axis and orthorhombic media. To stabilise the perturbation series and improve the accuracy, the Shanks transform is applied for all the cases. We select different parameterisations for both velocities and anellipticity parameters for an orthorhombic model. From the comparison in traveltime error, the new moveout approximations result in better accuracy comparing with the standard perturbation‐based methods and other approximations.     

Azimuthal Seismic Amplitude Variation with Offset and Azimuth Inversion in Weakly Anisotropic Media with Orthorhombic Symmetry

Seismic amplitude variation with offset and azimuth (AVOaz) inversion is well known as a popular and pragmatic tool utilized to estimate fracture parameters. A single set of vertical fractures aligned along a preferred horizontal direction embedded in a horizontally layered medium can be considered as an effective long-wavelength orthorhombic medium. Estimation of Thomsen’s weak-anisotropy (WA) parameters and fracture weaknesses plays an important role in characterizing the orthorhombic anisotropy in a weakly anisotropic medium. Our goal is to demonstrate an orthorhombic anisotropic AVOaz inversion approach to describe the orthorhombic anisotropy utilizing the observable wide-azimuth seismic reflection data in a fractured reservoir with the assumption of orthorhombic symmetry. Combining Thomsen’s WA theory and linear-slip model, we first derive a perturbation in stiffness matrix of a weakly anisotropic medium with orthorhombic symmetry under the assumption of small WA parameters and fracture weaknesses. Using the perturbation matrix and scattering function, we then derive an expression for linearized PP-wave reflection coefficient in terms of P- and S-wave moduli, density, Thomsen’s WA parameters, and fracture weaknesses in such an orthorhombic medium, which avoids the complicated nonlinear relationship between the orthorhombic anisotropy and azimuthal seismic reflection data. Incorporating azimuthal seismic data and Bayesian inversion theory, the maximum a posteriori solutions of Thomsen’s WA parameters and fracture weaknesses in a weakly anisotropic medium with orthorhombic symmetry are reasonably estimated with the constraints of Cauchy a priori probability distribution and smooth initial models of model parameters to enhance the inversion resolution and the nonlinear iteratively reweighted least squares strategy. The synthetic examples containing a moderate noise demonstrate the feasibility of the derived orthorhombic anisotropic AVOaz inversion method, and the real data illustrate the inversion stabilities of orthorhombic anisotropy in a fractured reservoir.     

Slowness-domain kinematical characteristics for horizontally layered orthorhombic media. Part I: Critical slowness match

Kinematical characteristics of reflected waves in anisotropic elastic media play an important role in the seismic imaging workflow. Considering compressional and converted waves, we derive new, azimuthally dependent, slowness-domain approximations for the kinematical characteristics of reflected waves (radial and transverse offsets, intercept time and traveltime) for layered orthorhombic media with varying azimuth of the vertical symmetry planes. The proposed method can be considered an extension of the well-known ‘generalized moveout approximation’ in the slowness domain, from azimuthally isotropic to azimuthally anisotropic models. For each slowness azimuth, the approximations hold for a wide angle range, combining power series coefficients in the vicinity of both the normal-incidence ray and an additional wide-angle ray. We consider two cases for the wide-angle ray: a ‘critical slowness match’ and a ‘pre-critical slowness match’ studied in Parts I and II of this work, respectively. For the critical slowness match, the approximations are valid within the entire slowness range, up to the critical slowness. For the ‘pre-critical slowness match’, the approximations are valid only within the bounded slowness range; however, the accuracy within the defined range is higher. The critical slowness match is particularly effective when the subsurface model includes a dominant high-velocity layer where, for nearly critical slowness values, the propagation in this layer is almost horizontal. Comparing the approximated kinematical characteristics with those computed by numerical ray tracing, we demonstrate high accuracy.     

THE INSENSITIVITY OF REFLECTED SH WAVES TO ANISOTROPY IN AN UNDERLYING LAYERED MEDIUM1

Propagation in the plane of mirror symmetry of a monoclinic medium, with displacement normal to the plane, is the most general circumstance in anisotropic media for which pure shear-wave propagation can occur at all angles. Because the pure shear mode is uncoupled from the other two modes, its slowness surface in the plane is an ellipse. When the mirror symmetry plane is vertical the pure shear waves in this plane are SH waves and the elliptical SH sheet of the slowness surface is, in general, tilted with respect to the vertical axis. Consider a half-space of such a monoclinic medium, called medium M, overlain by a half-space of isotropic medium I with plane SH waves incident on medium M propagating in the vertical symmetry plane of M. Contrary to the appearance of a lack of symmetry about the vertical axis due to the tilt of the SH-wave slowness ellipse, the reflection and transmission coefficients are symmetrical functions of the angle of incidence, and further, there exists an isotropic medium E with uniquely determined density and shear speed which gives exactly the same reflection and transmission coefficients underlying medium J as does monoclinic medium M. This means that the underlying monoclinic medium M can be replaced by isotropic medium E without changing the reflection and transmission coefficients for all values of the angle of incidence. Thus no set of SH seismic experiments performed in the isotropic medium in the symmetry plane of the underlying half-space can reveal anything about the monoclinic anisotropy of that underlying half-space. Moreover, even when the underlying monoclinic half-space is stratified, there exists a stratified isotropic half-space that gives the identical reflection coefficient as the stratified monoclinic half-space for all angles of incidence and all frequencies.     

APPROXIMATIONS TO SHEAR-WAVE VELOCITY AND MOVEOUT EQUATIONS IN ANISOTROPIC MEDIA1

Backus and Crampin derived analytical equations for estimating approximate phase-velocity variations in symmetry planes in weakly anisotropic media, where the coefficients of the equations are linear combinations of the elastic constants. We examine the application of similar equations to group-velocity variations in off-symmetry planes, where the coefficients of the equations are derived numerically. We estimate the accuracy of these equations over a range of anisotropic materials with transverse isotropy with both vertical and horizontal symmetry axes, and with combinations of transverse isotropy yielding orthorhombic symmetry. These modified equations are good approximations for up to 17% shear-wave anisotropy for propagations in symmetry planes for all waves in all symmetry systems examined, but are valid only for lower shear-wave anisotropy (up to 11%) in off-symmetry planes. We also obtain analytical moveout equations for the reflection of qP-, qSH-, and qSV- waves from a single interface for off-symmetry planes in anisotropic symmetry. The moveout equation consists of two terms: a hyperbolic moveout and a residual moveout, where the residual moveout is proportional to the degree of anisotropy and the spread length of the acquisition geometry. Numerical moveout curves are computed for a range of anisotropic materials to verify the analytical moveout equations.     

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.     

The kinematics of S waves in acoustic orthorhombic media

Orthorhombic models are often used in the seismic industry nowadays to describe azimuthal and polar anisotropy and reasonably realistic in capturing the features of the earth interior. It is challenging to handle so many model parameters in the seismic data processing. In order to reduce the number of the parameters for P wave, the acoustic orthorhombic medium is proposed by setting all on-axis S wave velocities to zero. However, due to the coupled behaviour for P and S waves in the orthorhombic model, the ‘S wave artefacts’ are still remained in the acoustic orthorhombic model, which kinematics needs to be defined and analysed. In this paper, we analyse the behaviour of S wave in acoustic orthorhombic media. By analysis of the slowness surface in acoustic orthorhombic media, we define the S waves (or S wave artefacts) that are more complicated in shape comparing to the one propagating in an acoustic transversely isotropic medium with a vertical symmetry axis. The kinematic properties of these waves are defined and analysed in both phase and group domain. The caustics, amplitude and the multi-layered case for S wave in acoustic orthorhombic model are also discussed. It is shown that there are two waves propagating in this acoustic orthorhombic medium. One of these waves is similar to the one propagating in acoustic vertical symmetry axis media, whereas another one has a very complicated shape consisting of two crossing surfaces.     

Low‐frequency wave propagation in periodically layered media

In this paper, we consider wave propagation in a layered medium. Using the Baker‐Campbell‐Hausdorff series, we expand the logarithm of a propagator matrix in series of frequency. The series coefficients allow us to extend the effective Backus medium for low frequencies. The proposed technique is applied to vertical propagation in a periodically layered and binary medium as well as for a gradient medium. The velocity dispersion equations are derived for these media. We also consider the layered medium with monoclinic anisotropy. We illustrate the accuracy of the proposed method on synthetic and well‐log data.     

忽略TTI介质对称轴倾角的可行性

假设横向各向同性(TI)介质的对称轴是垂直的(VTI)或者水平的(HTI)能给实际资料处理带来便利,然而实际TI介质的对称轴往往是倾斜的(TTI),忽略对称轴倾角可能给各向异性参数提取和成像带来偏差,因此需要研究是否能、以及什么条件下能忽略TTI介质对称轴倾角.本文通过理论研究和数值分析研究了与TTI介质弹性性质最接近的VTI介质(OAVTI)的弹性常数和各向异性参数与原TTI介质的弹性常数和各向异性参数之间的联系与差别.结果表明:OAVTI介质各向异性参数与原TTI介质各向异性参数之间的差别可统一表示成F(α₀/β₀,ε,δ,γ)ξ²的形式,其中F(α₀/β₀,ε,δ,γ)是无量纲各向异性参数(ε, δ, γ)的线性函数,ξ是对称轴倾角;ξ的大小对各参数的误差起主导作用,一般不建议忽略20°~25°以上的对称轴倾角;当ξ较小时,即使是对强各向异性的TTI介质作VTI近似,引起的P波各向异性参数误差也很小,因此在纵波资料处理中忽略TTI介质对称轴倾角通常是可行的;即使在小ξ条件下,倾斜对称轴对SV波也有显著影响,因此在转换波资料处理中,不建议忽略TTI介质的对称轴倾角.本文的研究为分析忽略TTI介质对称轴倾角的可行性提供了理论依据和简便的判据.     

Moveout approximation for horizontal transversely isotropic and vertical transversely isotropic layered medium. Part II: effective model‡

We use residual moveouts measured along continuous full azimuth reflection angle gathers, in order to obtain effective horizontal transversely isotropic model parameters. The angle gathers are generated through a special angle domain imaging system, for a wide range of reflection angles and full range of phase velocity azimuths. The estimation of the effective model parameters is performed in two stages. First, the background horizontal transversely isotropic (HTI)/vertical transversely isotropic (VTI) layered model is used, along with the values of reflection angles, for converting the measured residual moveouts (or traveltime errors) into azimuthally dependent normal moveout (NMO) velocities. Then we apply a digital Fourier transform to convert the NMO velocities into azimuthal wavenumber domain, in order to obtain the effective HTI model parameters: vertical time, vertical compression velocity, Thomsen parameter delta and the azimuth of the medium axis of symmetry. The method also provides a reliability criterion of the HTI assumption. The criterion shows whether the medium possesses the HTI type of symmetry, or whether the azimuthal dependence of the residual traveltime indicates to a more complex azimuthal anisotropy. The effective model used in this approach is defined for a 1D structure with a set of HTI, VTI and isotropic layers (with at least one HTI layer). We describe and analyse the reduction of a multi‐layer structure into an equivalent effective HTI model. The equivalent model yields the same NMO velocity and the same offset azimuth on the Earth's surface as the original layered structure, for any azimuth of the phase velocity. The effective model approximates the kinematics of an HTI/VTI layered structure using only a few parameters. Under the hyperbolic approximation, the proposed effective model is exact.