Seismic critical-angle reflectometry: A method to characterize azimuthal anisotropy?

Geophysics ◽  
2007 ◽  
Vol 72 (3) ◽  
pp. D41-D50 ◽  
Author(s):  
Martin Landrø ◽  
Ilya Tsvankin
Keyword(s):  
Critical Angle ◽  
Transverse Isotropy ◽  
Symmetry Plane ◽  
Azimuthal Anisotropy ◽  
P Wave ◽  
Head Wave ◽  
Orthorhombic Symmetry ◽  
Amplitude Curve ◽  
Diagnostic Features ◽  
Azimuthal Variation

Existing anisotropic parameter-estimation algorithms that operate with long-offset data are based on the inversion of either nonhyperbolic moveout or wide-angle amplitude-variation-with-offset (AVO) response. We show that valuable information about anisotropic reservoirs can also be obtained from the critical angle of reflected waves. To explain the behavior of the critical angle, we develop weak-anisotropy approximations for vertical transverse isotropy and then use Tsvankin’s notation to extend them to azimuthally anisotropic models of orthorhombic symmetry. The P-wave critical-angle reflection in orthorhombic media is particularly sensitive to the parameters [Formula: see text] and [Formula: see text] responsible for the symmetry-plane horizontal velocity in the high-velocity layer. The azimuthal variation of the critical angle for typical orthorhombic models can reach [Formula: see text], which translates into substantial changes in the critical offset of the reflected P-wave. The main diagnostic features of the critical-angle reflection employed in our method include the rapid amplitude increase at the critical angle and the subsequent separation of the head wave. Analysis of exact synthetic seismograms, generated with the reflectivity method, confirms that the azimuthal variation of the critical offset is detectable on wide-azimuth, long-spread data and can be qualitatively described by our linearized equations. Estimation of the critical offset from the amplitude curve of the reflected wave, however, is not straightforward. Additional complications may be caused by the overburden noise train and by the influence of errors in the overburden velocity model on the computation of the critical angle. Still, critical-angle reflectometry should help to constrain the dominant fracture directions and can be combined with other methods to reduce the uncertainty in the estimated anisotropy parameters.

3-D description of normal moveout in anisotropic inhomogeneous media

Geophysics ◽  
1998 ◽  
Vol 63 (3) ◽  
pp. 1079-1092 ◽  
Author(s):  
Vladimir Grechka ◽  
Ilya Tsvankin
Keyword(s):  
Fractured Reservoirs ◽  
Symmetry Plane ◽  
Transversely Isotropic ◽  
Inhomogeneous Media ◽  
P Wave ◽  
Orthorhombic Symmetry ◽  
Azimuthal Variation ◽  
Anisotropic Models ◽  
Transversely Isotropic Media ◽  
Isotropic Media

We present a new equation for normal‐moveout (NMO) velocity that describes azimuthally dependent reflection traveltimes of pure modes from both horizontal and dipping reflectors in arbitrary anisotropic inhomogeneous media. With the exception of anomalous areas such as those where common‐midpoint (CMP) reflection time decreases with offset, the azimuthal variation of NMO velocity represents an ellipse in the horizontal plane, with the orientation of the axes determined by the properties of the medium and the direction of the reflector normal. In general, a minimum of three azimuthal measurements is necessary to reconstruct the best‐fit ellipse and obtain NMO velocity in all azimuthal directions. This result provides a simple way to correct for the azimuthal variation in stacking velocity often observed in 3-D surveys. Even more importantly, analytic expressions for the parameters of the NMO ellipse can be used in the inversion of moveout data for the anisotropic coefficients of the medium. For homogeneous transversely isotropic media with a vertical axis of symmetry (VTI media), our equation for azimuthally dependent NMO velocity from dipping reflectors becomes a relatively simple function of phase velocity and its derivatives. We show that the zero‐dip NMO velocity Vnmo(0) and the anisotropic coefficient η are sufficient to describe the P-wave NMO velocity for any orientation of the CMP line with respect to the dip plane of the reflector. Using our formalism, Vnmo(0) and η (the only parameters needed for time processing) can be found from the dip‐dependent NMO velocity at any azimuth or, alternatively, from the azimuthally dependent NMO for a single dipping reflector. We also apply this theory to more complicated azimuthally anisotropic models with the orthorhombic symmetry used to describe fractured reservoirs. For reflections from horizontal interfaces in orthorhombic media, the axes of the normal moveout ellipse are aligned with the vertical symmetry planes. Therefore, azimuthal P-wave moveout measurements can be inverted for the orientation of the symmetry planes (typically determined by the fracture direction) and the NMO velocities within them. If the vertical velocity is known, symmetry‐plane NMO velocities make it possible to estimate two anisotropic parameters equivalent to Thomsen’s coefficient δ for transversely isotropic media.

Quartic moveout coefficient for a dipping azimuthally anisotropic layer

Geophysics ◽  
2004 ◽  
Vol 69 (3) ◽  
pp. 699-707 ◽  
Author(s):  
Andrés Pech ◽  
Ilya Tsvankin
Keyword(s):  
Reservoir Characterization ◽  
Fractured Reservoirs ◽  
Symmetry Plane ◽  
Exact Expression ◽  
Naturally Fractured Reservoirs ◽  
P Wave ◽  
Orthorhombic Symmetry ◽  
Fracture Density ◽  
Weak Anisotropy ◽  
Velocity Anisotropy

Interpretation and inversion of azimuthally varying nonhyperbolic reflection moveout requires accounting for both velocity anisotropy and subsurface structure. Here, our previously derived exact expression for the quartic moveout coefficient A4 is applied to P‐wave reflections from a dipping interface overlaid by a medium of orthorhombic symmetry. The weak‐anisotropy approximaton for the coefficient A4 in a homogeneous orthorhombic layer is controlled by the anellipticity parameters η(1), η(2), and η(3), which are responsible for time processing of P‐wave data. If the dip plane of the reflector coincides with the vertical symmetry plane [x1, x3], A4 on the dip line is proportional to the in‐plane anellipticity parameter η(2) and always changes sign for a dip of 30○. The quartic coefficient on the strike line is a function of all three η–parameters, but for mild dips it is mostly governed by η(1)—the parameter defined in the incidence plane [x2, x3]. Whereas the magnitude of the dip line A4 typically becomes small for dips exceeding 45○, the nonhyperbolic moveout on the strike line may remain significant even for subvertical reflectors. The character of the azimuthal variation of A4 depends on reflector dip and is quite sensitive to the signs and relative magnitudes of η(1), η(2), and η(3). The analytic results and numerical modeling show that the azimuthal pattern of the quartic coefficient can contain multiple lobes, with one or two azimuths of vanishing A4 between the dip and strike directions. The strong influence of the anellipticity parameters on the azimuthally varying coefficient A4 suggests that nonhyperbolic moveout recorded in wide‐azimuth surveys can help to constrain the anisotropic velocity field. Since for typical orthorhombic models that describe naturally fractured reservoirs the parameters η(1,2,3) are closely related to the fracture density and infill, the results of azimuthal nonhyperbolic moveout analysis can also be used in reservoir characterization.

scholarly journals Physical and numerical modeling of tuning and diffraction in azimuthally anisotropic media

Geophysics ◽  
2000 ◽  
Vol 65 (5) ◽  
pp. 1613-1621 ◽  
Author(s):  
Richard L. Gibson ◽  
Stephen Theophanis ◽  
M. Nafi Toksöz
Keyword(s):  
Numerical Modeling ◽  
Fractured Reservoirs ◽  
Homogeneous Model ◽  
Azimuthal Anisotropy ◽  
P Wave ◽  
Orthorhombic Symmetry ◽  
Fractured Reservoir ◽  
Sh Wave ◽  
Wave Data ◽  
Ray Methods

Fractured reservoirs are an important target for exploration and production geophysics, and the azimuthal anisotropy often associated with these reservoirs can strongly influence seismic wave propagation. We created a physical model of a fractured reservoir to simulate some of these propagation effects. The reservoir is represented by a phenolite disk that is thin with respect to the elastic wavelengths in the experiment, creating model dimensions that are representative of realistic reservoirs. Phenolite is strongly anisotropic with orthorhombic symmetry, which suggests that azimuthal amplitude versus offset (AVO) effects should be obvious in data. We acquired both SH- and P-wave data in common‐offset gathers with a near offset and a far offset and found that although the SH-wave data show clear azimuthal variations in AVO, the P-wave signals show no apparent changes with azimuth. We then applied numerical modeling to analyze the data. Because ray methods cannot model diffractions from the disk edge, we first used a ray‐Born technique to simulate variations in waveforms associated with such scattering. The synthetic seismograms reproduced variations in the SH-wave waveforms accurately, though the amplitude contrast between acquisition azimuths was overestimated. Assuming a laterally homogeneous model, we then applied ray methods to simulate tuning effects in SH- and P-wave data and confirmed that in spite of the large contrasts in elastic properties, the tuning of the P-wave reflections from the thin disk changed so there was negligible contrast in AVO with azimuth. Models of field scale reservoirs showed that the same effects could be expected for field applications.

Seismic traveltime analysis for azimuthally anisotropic media: Theory and experiment

Geophysics ◽  
1997 ◽  
Vol 62 (5) ◽  
pp. 1570-1582 ◽  
Author(s):  
Colin M. Sayers ◽  
Daniel A. Ebrom
Keyword(s):  
Model Simulation ◽  
Symmetry Plane ◽  
Transversely Isotropic ◽  
Bedding Plane ◽  
Seismic Profile ◽  
Azimuthal Anisotropy ◽  
P Wave ◽  
Principal Axes ◽  
Physical Model Simulation ◽  
Single Set

Natural fractures in reservoirs, and in the caprock overlying the reservoir, play an important role in determining fluid flow during production. The density and orientation of sets of fractures is therefore of great interest. Rocks possessing an anisotropic fabric and a preferred orientation of fractures display both polar and azimuthal anisotropy. Sedimentary rocks containing several sets of vertical fractures may be approximated as having monoclinic symmetry with symmetry plane parallel to the layers if, in the absence of fractures, the rock is transversely isotropic with symmetry axis perpendicular to the bedding plane. A nonhyperbolic traveltime equation, which can be used in the presence of azimuthally anisotropic layered media, can be obtained from an expansion of the inverse‐squared ray velocity in spherical harmonics. For a single set of aligned fractures, application of this equation to traveltime data acquired at a sufficient number of azimuths allows the strike of the fractures to be estimated. Analysis of the traveltimes measured in a physical model simulation of a reverse vertical seismic profile in an azimuthally anisotropic medium shows the medium to be orthorhombic with principal axes in agreement with those given by an independent shear‐wave experiment. In contrast to previous work, no knowledge of the orientation of the symmetry planes is required. The method is therefore applicable to P‐wave data collected at multiple azimuths using multiple offset vertical seismic profiling (VSP) techniques.

Estimation of fracture parameters from reflection seismic data—Part II: Fractured models with orthorhombic symmetry

Geophysics ◽  
2000 ◽  
Vol 65 (6) ◽  
pp. 1803-1817 ◽  
Author(s):  
Andrey Bakulin ◽  
Vladimir Grechka ◽  
Ilya Tsvankin
Keyword(s):  
Symmetry Axis ◽  
Fractured Reservoirs ◽  
Symmetry Plane ◽  
P Wave ◽  
Single Fracture ◽  
Orthorhombic Symmetry ◽  
Weak Anisotropy ◽  
Horizontal Symmetry ◽  
Vertical Symmetry Plane ◽  
Vertical Fractures

Existing geophysical and geological data indicate that orthorhombic media with a horizontal symmetry plane should be rather common for naturally fractured reservoirs. Here, we consider two orthorhombic models: one that contains parallel vertical fractures embedded in a transversely isotropic background with a vertical symmetry axis (VTI medium) and the other formed by two orthogonal sets of rotationally invariant vertical fractures in a purely isotropic host rock. Using the linear‐slip theory, we obtain simple analytic expressions for the anisotropic coefficients of effective orthorhombic media. Under the assumptions of weak anisotropy of the background medium (for the first model) and small compliances of the fractures, all effective anisotropic parameters reduce to the sum of the background values and the parameters associated with each fracture set. For the model with a single fracture system, this result allows us to eliminate the influence of the VTI background by evaluating the differences between the anisotropic parameters defined in the vertical symmetry planes. Subsequently, the fracture weaknesses, which carry information about the density and content of the fracture network, can be estimated in the same way as for fracture‐induced transverse isotropy with a horizontal symmetry axis (HTI media) examined in our previous paper (part I). The parameter estimation procedure can be based on the azimuthally dependent reflection traveltimes and prestack amplitudes of P-waves alone if an estimate of the ratio of the P- and S-wave vertical velocities is available. It is beneficial, however, to combine P-wave data with the vertical traveltimes, NMO velocities, or AVO gradients of mode‐converted (PS) waves. In each vertical symmetry plane of the model with two orthogonal fracture sets, the anisotropic parameters are largely governed by the weaknesses of the fractures orthogonal to this plane. For weak anisotropy, the fracture sets are essentially decoupled, and their parameters can be estimated by means of two independently performed HTI inversions. The input data for this model must include the vertical velocities (or reflector depth) to resolve the anisotropic coefficients in each vertical symmetry plane rather than just their differences. We also discuss several criteria that can be used to distinguish between the models with one and two fracture sets. For example, the semimajor axis of the P-wave NMO ellipse and the polarization direction of the vertically traveling fast shear wave are always parallel to each other for a single system of fractures, but they may become orthogonal in the medium with two fracture sets.

Nonhyperbolic reflection moveout for horizontal transverse isotropy

Geophysics ◽  
1998 ◽  
Vol 63 (5) ◽  
pp. 1738-1753 ◽  
Author(s):  
AbdulFattah Al‐Dajani ◽  
Ilya Tsvankin
Keyword(s):  
Transverse Isotropy ◽  
Single Layer ◽  
Azimuthal Velocity ◽  
Azimuthal Anisotropy ◽  
P Wave ◽  
Analytic Description ◽  
Pure Mode ◽  
Quartic Term ◽  
Hti Media ◽  
Vertical Transverse Isotropy

The transversely isotropic model with a horizontal axis of symmetry (HTI) has been used extensively in studies of shear‐wave splitting to describe fractured formations with a single system of parallel vertical penny‐shaped cracks. Here, we present an analytic description of longspread reflection moveout in horizontally layered HTI media with arbitrary strength of anisotropy. The hyperbolic moveout equation parameterized by the exact normal‐moveout (NMO) velocity is sufficiently accurate for P-waves on conventional‐length spreads (close to the reflector depth), although the NMO velocity is not, in general, usable for converting time to depth. However, the influence of anisotropy leads to the deviation of the moveout curve from a hyperbola with increasing spread length, even in a single‐layer model. To account for nonhyperbolic moveout, we have derived an exact expression for the azimuthally dependent quartic term of the Taylor series traveltime expansion [t2(x2)] valid for any pure mode in an HTI layer. The quartic moveout coefficient and the NMO velocity are then substituted into the nonhyperbolic moveout equation of Tsvankin and Thomsen, originally designed for vertical transverse isotropy (VTI). Numerical examples for media with both moderate and uncommonly strong nonhyperbolic moveout show that this equation accurately describes azimuthally dependent P-wave reflection traveltimes in an HTI layer, even for spread lengths twice as large as the reflector depth. In multilayered HTI media, the NMO velocity and the quartic moveout coefficient reflect the influence of layering as well as azimuthal anisotropy. We show that the conventional Dix equation for NMO velocity remains entirely valid for any azimuth in HTI media if the group‐velocity vectors (rays) for data in a common‐midpoint (CMP) gather do not deviate from the vertical incidence plane. Although this condition is not exactly satisfied in the presence of azimuthal velocity variations, rms averaging of the interval NMO velocities represents a good approximation for models with moderate azimuthal anisotropy. Furthermore, the quartic moveout coefficient for multilayered HTI media can also be calculated with acceptable accuracy using the known averaging equations for vertical transverse isotropy. This allows us to extend the nonhyperbolic moveout equation to horizontally stratified media composed of any combination of isotropic, VTI, and HTI layers. In addition to providing analytic insight into the behavior of nonhyperbolic moveout, these results can be used in modeling and inversion of reflection traveltimes in azimuthally anisotropic media.

Synthesis of multicomponent quasi-P and converted quasi-P-S seismograms for intersecting fracture systems

Geophysics ◽  
2000 ◽  
Vol 65 (4) ◽  
pp. 1261-1271 ◽  
Author(s):  
Andrey A. Ortega ◽  
George A. McMechan
Keyword(s):  
Shear Wave ◽  
Transverse Isotropy ◽  
Normal Incidence ◽  
Shear Wave Splitting ◽  
P Wave ◽  
Orthorhombic Symmetry ◽  
S Wave ◽  
Stiffness Constant ◽  
Wave Splitting ◽  
Vertical Fractures

Dynamic ray shooting with interpolation is an economical way of computing approximate Green’s functions in 3-D heterogeneous anisotropic media. The amplitudes, traveltimes, and polarizations of the reflected rays arriving at the surface are interpolated to synthesize three‐component seismograms at the desired recording points. The algorithm is applied to investigate kinematic quasi-P-wave propagation and converted quasi-P-S-wave splitting variations produced in reflections from the bottom of a layer containing two sets of intersecting dry vertical fractures as a function of the angle between the fracture sets and of the intensity of fracturing. An analytical expression is derived for the stiffness constant C16 that extends Hudson’s second‐order scattering theory to include tetragonal-2 symmetry systems. At any offset, the amount of splitting in nonorthogonal (orthorhombic symmetry) intersecting fracture sets is larger than in orthogonal (tetragonal-1 symmetry) systems, and it increases nonlinearly as a function of the intensity of fracturing as offset increases. Such effects should be visible in field data, provided that the dominant frequency is sufficiently high and the offset is sufficiently large. The amount of shear‐wave splitting at vertical incidence increases nonlinearly as a function of the intensity of fracturing and increases nonlinearly from zero in the transition from tetragonal-1 anisotropy through orthorhombic to horizontal transverse isotropy; the latter corresponds to the two crack systems degenerating to one. The zero shear‐wave splitting corresponds to a singularity, at which the vertical velocities of the two quasi‐shear waves converge to a single value that is both predicted theoretically and illustrated numerically. For the particular case of vertical fractures, there is no P-to-S conversion of vertically propagating (zero‐offset) waves. If the fractures are not vertical, the normal incidence P-to-S reflection coefficient is not zero and thus is a potential diagnostic of fracture orientation.

Moveout inversion of P-wave data for horizontal transverse isotropy

Geophysics ◽  
1999 ◽  
Vol 64 (4) ◽  
pp. 1219-1229 ◽  
Author(s):  
Pedro Contreras ◽  
Vladimir Grechka ◽  
Ilya Tsvankin
Keyword(s):  
Vertical Velocity ◽  
Symmetry Axis ◽  
Transverse Isotropy ◽  
P Wave ◽  
Orthorhombic Symmetry ◽  
Anisotropic Parameter ◽  
Data Set ◽  
Velocity Models ◽  
Wave Data ◽  
Hti Media

The transversely isotropic model with a horizontal symmetry axis (HTI media) has been extensively used in seismological studies of fractured reservoirs. In this paper, a parameter‐estimation technique originally developed by Grechka and Tsvankin for the more general orthorhombic media is applied to horizontal transverse isotropy. Our methodology is based on the inversion of azimuthally‐dependent P-wave normal‐moveout (NMO) velocities from horizontal and dipping reflectors. If the NMO velocity of a given reflection event is plotted in each azimuthal direction, it forms an ellipse determined by three combinations of medium parameters. The NMO ellipse from a horizontal reflector in HTI media can be inverted for the azimuth β of the symmetry axis, the vertical velocity [Formula: see text], and the Thomsen‐type anisotropic parameter δ(V). We describe a technique for obtaining the remaining (for P-waves) anisotropic parameter η(V) (or ε(V)) from the NMO ellipse corresponding to a dipping reflector of arbitrary azimuth. The interval parameters of vertically inhomogeneous HTI media are recovered using the generalized Dix equation that operates with NMO ellipses for horizontal and dipping events. High accuracy of our method is confirmed by inverting a synthetic multiazimuth P-wave data set generated by ray tracing for a layered HTI medium with depth‐varying orientation of the symmetry axis. Although estimation of η(V) can be carried out by the algorithm developed for orthorhombic media, for more stable results the HTI model has to be used from the outset of the inversion procedure. It should be emphasized that P-wave conventional‐spread moveout data provide enough information to distinguish between HTI and lower‐symmetry models. We show that if the medium has the orthorhombic symmetry and is sufficiently different from HTI, the best‐fit HTI model cannot match the NMO ellipses for both a horizontal and a dipping event. The anisotropic coefficients responsible for P-wave moveout can be combined to estimate the crack density and predict whether the cracks are fluid‐filled or dry. A unique feature of the HTI model that distinguishes it from both vertical transverse isotropy and orthorhombic media is that moveout inversion provides not just zero‐dip NMO velocities and anisotropic coefficients, but also the true vertical velocity. As a result, reflection P-wave data acquired over HTI formations can be used to build velocity models in depth and perform anisotropic depth processing.

P-wave reflection-moveout approximation for horizontally layered media of arbitrary moderate anisotropy

Geophysics ◽  
2020 ◽  
Vol 85 (2) ◽  
pp. C61-C70 ◽  
Author(s):  
Véronique Farra ◽  
Ivan Pšenčík
Keyword(s):  
Coordinate System ◽  
Transverse Isotropy ◽  
Local Coordinate System ◽  
Surface Profile ◽  
Taylor Series Expansion ◽  
Layered Media ◽  
P Wave ◽  
Orthorhombic Symmetry ◽  
Weak Anisotropy ◽  
Cartesian Coordinate

We present an approximate nonhyperbolic P-wave moveout formula applicable to horizontally layered media of moderate or weak anisotropy of arbitrary symmetry and orientation. Anisotropy symmetry and its orientation may differ from layer to layer. Instead of commonly used Taylor-series expansion of the square of the reflection traveltime in terms of the square of the offset, we use the weak-anisotropy approximation, in which the square of the reflection traveltime is expanded in terms of weak-anisotropy (WA) parameters. The resulting formula is simple, and it provides a transparent relation between the traveltimes and WA parameters. Along an arbitrarily chosen single surface profile, it depends, in each layer, on the thickness of the layer, on the reference P-wave velocity used for the construction of reference rays, and on three WA parameters specified in the Cartesian coordinate system related to the profile. In each layer, these three “profile” WA parameters depend on “local” WA parameters specifying anisotropy of a given layer in a local coordinate system and on directional cosines specifying the orientation of the local coordinate system with respect to the profile one. The number of local P-wave WA parameters may vary from three for transverse isotropy or six for orthorhombic symmetry to nine for triclinic symmetry. Our tests of the accuracy indicate that the maximum relative traveltime errors do not exceed 0.5% or 2.5% for weak or moderate P-wave anisotropy, respectively.

Special section on azimuthal dependence of P-wave seismic signatures—Introduction

Geophysics ◽  
1999 ◽  
Vol 64 (4) ◽  
pp. 1139-1142 ◽  
Author(s):  
Ilya Tsvankin, ◽  
Heloise B. Lynn,
Keyword(s):  
Special Section ◽  
Azimuthal Anisotropy ◽  
P Wave ◽  
Special Issue ◽  
Azimuthal Dependence ◽  
Azimuthal Variation ◽  
Main Motivation ◽  
Development Scenarios ◽  
The Right ◽  
Lateral Heterogeneities

This special issue is based on papers presented at the post‐convention SEG workshop on azimuthal dependence of P-wave signatures held in Dallas in 1997. The main motivation for analyzing the azimuthal variation of seismic traveltimes, amplitudes, attenuation, etc. is to obtain reliable information about azimuthal anisotropy in the subsurface. Another potential application of multiazimuth techniques is in finding and mapping desirable lateral heterogeneities that could “masquerade” as azimuthal anisotropy. The last topic has not yet been fully discussed (and is not addressed in the special issue), but several exploration and development scenarios contain oriented lateral heterogeneities (sand channels, etc.) that at the right scale length could be highly visible in properly processed wide‐azimuth 3-D data.

Sign in / Sign up

Export Citation Format

Share Document