Depth migration anisotropy analysis in the time domain

In areas of complex geology such as the Canadian Foothills, the effects of anisotropy are apparent in seismic data and estimation of anisotropic parameters for use in seismic imaging is not a trivial task. Here we explore the applicability of common‐focus point (CFP)‐based velocity analysis to estimate anisotropic parameters for the variably tilted shale thrust sheet in the Canadian Foothills model. To avoid the inherent velocity‐depth ambiguity, we assume that the elastic properties of thrust‐sheet with respect to transverse isotropy symmetry axis are homogeneous, the reflector below the thrust‐sheet is flat, and that the anisotropy is weak. In our CFP approach to velocity analysis, for a poorly imaged reflection point, a traveltime residual is obtained as the time difference between the focusing operator for an assumed subsurface velocity model and the corresponding CFP response obtained from the reflection data. We assume that this residual is due to unknown values for anisotropy, and we perform an iterative linear inversion to obtain new model parameters that minimize the residuals. Migration of the data using parameters obtained from our inversion results in a correctly positioned and better focused reflector below the thrust sheet. For traveltime computation we use a brute force mapping scheme that takes into account weakly tilted transverse isotropy media. For inversion, the problem is set up as a generalized Newton's equation where traveltime error (differential time shift) is linearly dependent on the parameter updates. The iterative updates of parameters are obtained by a least‐squares solution of Newton's equations. The significance of this work lies in its applicability to areas where transverse isotropy layers are heterogeneous laterally, and where transverse isotropy layers are overlain by complex structures that preclude a moveout curve fitting.     

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.     

Traveltime inversion and error analysis for layered anisotropy

While tilted transverse isotropy (TTI) is a good approximation of the velocity structure for many dipping and fractured strata, it is still challenging to estimate anisotropic depth models even when the tilted angle is known. With the assumption of weak anisotropy, we present a TTI traveltime inversion approach for models consisting of several thickness-varying layers where the anisotropic parameters are constant for each layer. For each model layer the inversion variables consist of the anisotropic parameters ε and δ, the tilted angle φ of its symmetry axis, layer velocity along the symmetry axis, and thickness variation of the layer. Using this method and synthetic data, we evaluate the effects of errors in some of the model parameters on the inverted values of the other parameters in crosswell and Vertical Seismic Profile (VSP) acquisition geometry. The analyses show that the errors in the layer symmetry axes sensitively affect the inverted values of other parameters, especially δ. However, the impact of errors in δ on the inversion of other parameters is much less than the impact on δ from the errors in other parameters. Hence, a practical strategy is first to invert for the most error-tolerant parameter layer velocity, then progressively invert for ε in crosswell geometry or δ in VSP geometry.     

Migration velocity analysis for tilted transversely isotropic media

Tilted transversely isotropic formations cause serious imaging distortions in active tectonic areas (e.g., fold‐and‐thrust belts) and in subsalt exploration. Here, we introduce a methodology for P‐wave prestack depth imaging in tilted transversely isotropic media that properly accounts for the tilt of the symmetry axis as well as for spatial velocity variations. For purposes of migration velocity analysis, the model is divided into blocks with constant values of the anisotropy parameters ε and δ and linearly varying symmetry‐direction velocity V_P0 controlled by the vertical (k_z) and lateral (k_x) gradients. Since determination of tilt from P‐wave data is generally unstable, the symmetry axis is kept orthogonal to the reflectors in all trial velocity models. It is also assumed that the velocity V_P0 is either known at the top of each block or remains continuous in the vertical direction. The velocity analysis algorithm estimates the velocity gradients k_z and k_x and the anisotropy parameters ε and δ in the layer‐stripping mode using a generalized version of the method introduced by Sarkar and Tsvankin for factorized transverse isotropy with a vertical symmetry axis. Synthetic tests for several models typical in exploration (a syncline, uptilted shale layers near a salt dome and a bending shale layer) confirm that if the symmetry‐axis direction is fixed and V_P0 is known, the parameters k_z, k_x, ε and δ can be resolved from reflection data. It should be emphasized that estimation of ε in tilted transversely isotropic media requires using nonhyperbolic moveout for long offsets reaching at least twice the reflector depth. We also demonstrate that application of processing algorithms designed for a vertical symmetry axis to data from tilted transversely isotropic media may lead to significant misfocusing of reflectors and errors in parameter estimation, even when the tilt is moderate (30°). The ability of our velocity analysis algorithm to separate the anisotropy parameters from the velocity gradients can be also used in lithology discrimination and geologic interpretation of seismic data in complex areas.     

Converted-wave imaging in anisotropic media: theory and case studies

Common‐conversion‐point binning associated with converted‐wave (C‐wave) processing complicates the task of parameter estimation, especially in anisotropic media. To overcome this problem, we derive new expressions for converted‐wave prestack time migration (PSTM) in anisotropic media and illustrate their applications using both 2D and 3D data examples. The converted‐wave kinematic response in inhomogeneous media with vertical transverse isotropy is separated into two parts: the response in horizontally layered vertical transverse isotrophy media and the response from a point‐scatterer. The former controls the stacking process and the latter controls the process of PSTM. The C‐wave traveltime in horizontally layered vertical transverse isotrophy media is determined by four parameters: the C‐wave stacking velocity V_C2, the vertical and effective velocity ratios γ₀ and γ_eff, and the C‐wave anisotropic parameter χ_eff. These four parameters are referred to as the C‐wave stacking velocity model. In contrast, the C‐wave diffraction time from a point‐scatterer is determined by five parameters: γ₀, V_P2, V_S2, η_eff and ζ_eff, where η_eff and ζ_eff are, respectively, the P‐ and S‐wave anisotropic parameters, and V_P2 and V_S2 are the corresponding stacking velocities. V_P2, V_S2, η_eff and ζ_eff are referred to as the C‐wave PSTM velocity model. There is a one‐to‐one analytical link between the stacking velocity model and the PSTM velocity model. There is also a simple analytical link between the C‐wave stacking velocities V_C2 and the migration velocity V_Cmig, which is in turn linked to V_P2 and V_S2. Based on the above, we have developed an interactive processing scheme to build the stacking and PSTM velocity models and to perform 2D and 3D C‐wave anisotropic PSTM. Real data applications show that the PSTM scheme substantially improves the quality of C‐wave imaging compared with the dip‐moveout scheme, and these improvements have been confirmed by drilling.     

3D P‐wave traveltime computation in transversely isotropic media using layer‐by‐layer wavefront marching

Subsurface rocks (e.g. shale) may induce seismic anisotropy, such as transverse isotropy. Traveltime computation is an essential component of depth imaging and tomography in transversely isotropic media. It is natural to compute the traveltime using the wavefront marching method. However, tracking the 3D wavefront is expensive, especially in anisotropic media. Besides, the wavefront marching method usually computes the traveltime using the eikonal equation. However, the anisotropic eikonal equation is highly non‐linear and it is challenging to solve. To address these issues, we present a layer‐by‐layer wavefront marching method to compute the P‐wave traveltime in 3D transversely isotropic media. To simplify the wavefront tracking, it uses the traveltime of the previous depth as the boundary condition to compute that of the next depth based on the wavefront marching. A strategy of traveltime computation is designed to guarantee the causality of wave propagation. To avoid solving the non‐linear eikonal equation, it updates traveltime along the expanding wavefront by Fermat's principle. To compute the traveltime using Fermat's principle, an approximate group velocity with high accuracy in transversely isotropic media is adopted to describe the ray propagation. Numerical examples on 3D vertical transverse isotropy and tilted transverse isotropy models show that the proposed method computes the traveltime with high accuracy. It can find applications in modelling and depth migration.     

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.     

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.     

Local Determination of Weak Anisotropy Parameters from Walkaway VSP qP-Wave Data in the Java Sea Region

We apply the inversion scheme of Zheng and Peník (2002) to the walkaway VSP data of Horne and Leaney (2000) collected in the Java Sea region. The goal is a local determination of parameters of the medium surrounding the borehole receiver array. The inversion scheme is based on linearized equations expressing qP-wave slowness and polarization vectors in terms of weak anisotropy (WA) parameters. It thus represents an alternative approach to Horne and Leaney (2000), who based their procedure on inversion of the Christoffel equation using a global optimization method. The presented inversion scheme is independent of structural complexities in the overburden and of the orientation of the borehole. The inverson formula is local, and has therefore potential to separate effects of anisotropy from effects of inhomogeneity. The data used are components of the slowness vector along the receiver array and polarization vectors. The inversion is performed without any assumptions concerning the remaining components of the slowness vector. The inversion is made (a) assuming arbitrary anisotropy, i.e., without any assumptions about symmetry of the medium, (b) assuming transverse isotropy with a vertical axis of symmetry and (c) assuming isotropy of the medium. Inverted are the raw data as well as data, in which weighting is used to reduce the effect of outliers. It is found that the WA parameters _z, ₁₅ and ₃₅ are considerably more stable than the parameters _x and _x. The latter two parameters are also found to be strongly correlated. Weaker correlation is also found between the mentioned two parameters and _z. The results of inversion show clearly that the studied medium is not isotropic. They also seem to indicate that the studied medium does not possess the VTI symmetry.     

P‐wave velocity distribution in basalt flows of the Enni Formation in the Faroe Islands from refraction seismic analysis

The main objective of this work is to establish the applicability of shallow surface‐seismic traveltime tomography in basalt‐covered areas. A densely sampled ～1300‐m long surface seismic profile, acquired as part of the SeiFaBa project in 2003 ( Japsen et al. 2006 ) at Glyvursnes in the Faroe Islands, served as the basis to evaluate the performance of the tomographic method in basalt‐covered areas. The profile is centred at a ～700‐m deep well. V_P, V_S and density logs, a zero‐offset VSP, downhole‐geophone recordings and geological mapping in the area provided good means of control. The inversion was performed with facilities of the Wide Angle Reflection/Refraction Profiling program package ( Ditmar et al. 1999 ). We tested many inversion sequences while varying the inversion parameters. Modelled traveltimes were verified by full‐waveform modelling. Typically an inversion sequence consists in several iterations that proceed until a satisfactory solution is reached. However, in the present case with high velocity contrasts in the subsurface we obtained the best result with two iterations: first obtaining a smooth starting model with small traveltime residuals by inverting with a high smoothing constraint and then inverting with the lowest possible smoothing constraint to allow the inversion to have the full benefit of the traveltime residuals. The tomogram gives usable velocity information for the near‐surface geology in the area but fails to reproduce the expected velocity distribution of the layered basalt flows. Based on the analysis of the tomogram and geological mapping in the area, a model was defined that correctly models first arrivals from both surface seismic data and downhole‐geophone data.     

THREE-TERM TAYLOR SERIES FOR t2 - x2-CURVES OF P- AND S-WAVES OVER LAYERED TRANSVERSELY ISOTROPIC GROUND*

The arrival-time curve of a reflection from a horizontal interface, beneath a homogeneous isotropic layer, is a hyperbola in the x - t-domain. If the subsurface is one-dimensionally inhomogeneous (horizontally layered), or if some or all of the layers are transversely isotropic with vertical axis of symmetry, the statement is no longer strictly true, though the arrival-time curves are still hyperbola-like. In the case of transverse isotropy, however, classical interpretation of these curves fails. Interval velocities calculated from t² - x²-curves do not always approximate vertical velocities and therefore cannot be used to calculate depths of reflectors. To study the relationship between velocities calculated from t² - x²-curves and the true velocities of a transversely isotropic layer, we approximate t² - x²-curves over a vertically inhomogeneous transversely isotropic medium by a three-term Taylor series and calculate expressions for these terms as a function of the elastic parameters. It is shown that both inhomogeneity and transverse isotropy affect slope and curvature of t² - x²-curves. For P-waves the effect of transverse isotropy is that the t² - x²-curves are convex upwards; for SV-waves the curves are convex downwards. For SH-waves transverse isotropy has no effect on curvature.     

Accounting for velocity anisotropy in seismic traveltime tomography: a case study from the investigation of the foundations of a Byzantine monumental building

We estimate velocity anisotropy factors from seismic traveltime tomographic data and apply a correction for anisotropy in the inversion procedure to test possible improvements on the traveltime fit and the quality of the resulting tomographic images. We applied the anisotropy correction on a traveltime data set obtained from the investigation of the foundation structure of a monumental building: a Byzantine church from the 11th century AD, in Athens, Greece. Vertical transverse isotropy is represented by one axis of symmetry and one anisotropy magnitude for the entire tomographic inversion grid. We choose the vertical direction for the symmetry axis by analysing the available data set and taking into account information on the character of the foundations of the church from the literature and past excavations. The anisotropy magnitude is determined by testing a series of values of anisotropy and examining their effect on the tomographic inversion results. The best traveltime fit and image quality are obtained with an anisotropy value (V_max/V_min) of 1.6, restricted to the high velocity structures in the subsurface. We believe that this anisotropy value, which is significantly higher than the usual values reported for near‐surface geological material, is related to the fabric of the church foundations, due to the shape of the individual stone blocks and the layout of the stonework. Inversion results obtained with the correction for anisotropy indicate that both the traveltime fit and the image quality are improved, providing an enhanced reconstruction of the velocity field, especially for the high‐velocity features. Based on this enhanced and more reliable reconstruction of velocity distribution, an improved image of the subsurface material character was made possible. In particular, the pattern and state of the church foundations and possible weak ground material areas were revealed more clearly. This improved subsurface knowledge may assist in a better design of restoration measures for monumental buildings such as Byzantine churches.     

Kinematical characteristics of the factorized velocity model

The factorized velocity model that incorporates both vertical heterogeneity and constant anisotropy is one of the complicated analytical models used in seismic data processing and interpretation. In this paper, I derive the analytic equations for offset, traveltime and relative geometrical spreading for the quasi‐compressional (qP‐) waves that can be used for modelling and inversion of the traveltime parameters. I show that the presence of anelliptic anisotropy usually dominates over the vertical heterogeneity with respect to the non‐hyperbolicity of the factorized velocity model.     

VSP traveltime inversion for anisotropy in a buried layer

We present a method for calculating the anisotropy parameter of a buried layer by inverting the total traveltimes of direct arrivals travelling from a surface source to a well‐bore receiver in a vertical seismic profiling (VSP) geometry. The method assumes two‐dimensional media. The medium above the layer of interest (and separated from it by a horizontal interface) can exhibit both anisotropy and inhomogeneity. Both the depth of the interface as well as the velocity field of the overburden are assumed to be known. We assume the layer of interest to be homogeneous and elliptically anisotropic, with the anisotropy described by a single parameter χ. We solve the function describing the traveltime between source and receiver explicitly for χ. The solution is expressed in terms of known quantities, such as the source and receiver locations, and in terms of quantities expressed as functions of the single argument x_r, which is the horizontal coordinate of the refraction point on the interface. In view of Fermat's principle, the measured traveltime T possesses a stationary value or, considering direct arrivals, a minimum value, . This gives rise to a key result ‐‐ the condition that the actual anisotropy parameter . Owing to the explicit expression , this result allows a direct calculation of in the layer of interest. We perform an error analysis and show this inverse method to be stable. In particular, for horizontally layered media, a traveltime error of one millisecond results in a typical error of about 20% in the anisotropy parameter. This is almost one order of magnitude less than the error inherent in the slowness method, which uses a similar set of experimental data. We conclude by detailing possible extensions to non‐elliptical anisotropy and a non‐planar interface.     

TTI介质的相速度研究

推导了二维TTI介质的相速度表达式,并且依据推导出来的相速度表达式,模拟并分析了二维TTI介质相速度的传播快照以及TI介质相速度的传播快照;对比并分析了TTI介质和TI介质模型的相速度理论计算值的X分量特征的差异。TTI介质的相速度研究具有较高的理论研究价值和实际应用价值．     

Seismic anisotropy analysis at the Low-Noise Underground Laboratory (LSBB) of Rustrel (France)

Seismic anisotropy of a fractured karstic limestone massif in sub-parallel underground galleries is studied. As the fractures are mostly vertically oriented, the seismic properties of the massif are approximated by horizontal transverse isotropy (HTI). Several data inversion methods were applied to a seismic dataset of arrival-times of P and S-waves.The applied methods include: isotropic tomography, simple cosine function fit, homogeneous Monte-Carlo anisotropic inversion for the parameters of horizontal transverse isotropy and anisotropic tomography for tilted transversely isotropic bodies. All methods lead to the conclusion that there is indeed an anisotropy present in the rock massif and confirm the direction of maximum velocity parallel to the direction of fracturing. Strong anisotropy of about 15% is found in the studied area. Repeated measurements show variations of the P-wave parameters, but not of the S-wave parameters, which is reflecting a change in water saturation.     

Weak‐anisotropy approximation for P‐wave reflection coefficient at the boundary between two tilted transversely isotropic media

Existing and commonly used in industry nowadays, closed‐form approximations for a P‐wave reflection coefficient in transversely isotropic media are restricted to cases of a vertical and a horizontal transverse isotropy. However, field observations confirm the widespread presence of rock beds and fracture sets tilted with respect to a reflection boundary. These situations can be described by means of the transverse isotropy with an arbitrary orientation of the symmetry axis, known as tilted transversely isotropic media. In order to study the influence of the anisotropy parameters and the orientation of the symmetry axis on P‐wave reflection amplitudes, a linearised 3D P‐wave reflection coefficient at a planar weak‐contrast interface separating two weakly anisotropic tilted tranversely isotropic half‐spaces is derived. The approximation is a function of the incidence phase angle, the anisotropy parameters, and symmetry axes tilt and azimuth angles in both media above and below the interface. The expression takes the form of the well‐known amplitude‐versus‐offset “Shuey‐type” equation and confirms that the influence of the tilt and the azimuth of the symmetry axis on the P‐wave reflection coefficient even for a weakly anisotropic medium is strong and cannot be neglected. There are no assumptions made on the symmetry‐axis orientation angles in both half‐spaces above and below the interface. The proposed approximation can be used for inversion for the model parameters, including the orientation of the symmetry axes. Obtained amplitude‐versus‐offset attributes converge to well‐known approximations for vertical and horizontal transverse isotropic media derived by Rüger in corresponding limits. Comparison with numerical solution demonstrates good accuracy.     

A multi‐azimuth seismic refraction study for a horizontal transverse isotropic medium: physical modelling results

To investigate the characteristics of the anisotropic stratum, a multi‐azimuth seismic refraction technique is proposed in this study since the travel time anomaly of the refraction wave induced by this anisotropic stratum will be large for a far offset receiver. To simplify the problem, a two‐layer (isotropy–horizontal transverse isotropy) model is considered. A new travel time equation of the refracted P‐wave propagation in this two‐layer model is derived, which is the function of the phase and group velocities of the horizontal transverse isotropic stratum. In addition, the measured refraction wave velocity in the physical model experiment is the group velocity. The isotropic intercept time equation of a refraction wave can be directly used to estimate the thickness of the top (isotropic) layer of the two‐layer model because the contrast between the phase and group velocities of the horizontal transverse isotropic medium is seldom greater than 10% in the Earth. If the contrast between the phase and group velocities of an anisotropic medium is small, the approximated travel time equation of a refraction wave is obtained. This equation is only dependent on the group velocity of the horizontal transverse isotropic stratum. The elastic constants A₁₁, A₁₃, and A₃₃ and the Thomsen anisotropic parameter ε of the horizontal transverse isotropic stratum can be estimated using this multi‐azimuth seismic refraction technique. Furthermore, under a condition of weak anisotropy, the Thomsen anisotropic parameter δ of the horizontal transverse isotropic stratum can be estimated by this technique as well.     

Post‐stack diffraction imaging in vertical transverse isotropy media using non‐hyperbolic moveout approximations

Diffractions carry valuable information about local discontinuities and small‐scale objects in the subsurface. They are still not commonly used in the process of geological interpretation. Many diffraction imaging techniques have been developed and applied for isotropic media, whereas relatively few techniques have been developed for anisotropic media. Ignoring anisotropy can result in low‐resolution images with wrongly positioned or spurious diffractors. In this article, we suggest taking anisotropy into account in two‐dimensional post‐stack domain by considering P‐wave non‐hyperbolic diffraction traveltime approximations for vertical transverse isotropy media, previously developed for reflection seismology. The accuracy of the final images is directly connected to the accuracy of the diffraction traveltime approximations. We quantified the accuracy of six different approximations, including hyperbolic moveout approximation, by the application of a post‐stack diffraction imaging technique on two‐dimensional synthetic data examples.