A weak-anisotropy approximation to multicomponent induction responses in cross-bedded formations

Geophysics ◽  
2006 ◽  
Vol 71 (4) ◽  
pp. F61-F66 ◽  
Author(s):  
Tsili Wang
Keyword(s):  
Apparent Resistivity ◽  
Transversely Isotropic ◽  
Anisotropy Ratio ◽  
Weak Anisotropy ◽  
Bedding Planes ◽  
Geometric Average ◽  
Induction Logging ◽  
Anisotropic Formation ◽  
The Cross

The multicomponent induction logging response to a cross-bedded formation has been modeled under a weak-anisotropy approximation. With the approximation, a cross-bedded formation can be modeled as a transversely isotropic (TI) medium. The validity of the approximation has been tested for the main (coplanar and coaxial) components of the induction response. The conditions for the weak-anisotropy approximation to be valid depend on the component of the response. For the coplanar components, the approximation is valid for an anisotropy ratio up to 2 if the relative dipping angle between the cross-bedded formation and the borehole axis is below [Formula: see text]. For the coaxial component, the approximation reduces to a previously established result that the apparent resistivity for such a component is the geometric average of the resistivities, parallel and perpendicular to the bedding planes of an anisotropic formation, respectively, if the borehole is ignored. Hence, the approximation holds for the coaxial component regardless of the anisotropy ratio.

Modeling of triaxial induction logging responses in multilayered anisotropic formations

Geophysics ◽  
2021 ◽  
pp. 1-46
Author(s):  
Yunyun Hu ◽  
Qingtao Sun
Keyword(s):  
Fourier Transform ◽  
Analytical Method ◽  
Inverse Fourier Transform ◽  
Singular Term ◽  
Transversely Isotropic ◽  
Spectral Domain ◽  
Induction Logging ◽  
Anisotropic Formation ◽  
Electric Anisotropy ◽  
Induction Logs

Triaxial induction logging tools have been widely applied to formation characterization due to its sensitivity to electric anisotropy. To model triaxial induction logs in multilayered general anisotropic formations, where the anisotropy can be arbitrary, an analytical method is applied to compute the tool responses. For the analytical method, Maxwell's equations in the spectral domain are written into a compact first-order differential equation. The equation is solved to obtain the spectral-domain fields, which are transformed to the spatial domain through the inverse Fourier transform. The singular issue for the tool located in highly deviated wells, is handled by subtracting the singular term in the spectral domain. The singularity treatment makes the integrands in the inverse Fourier transform decay faster, thus making the infinite integration computation faster. Formations with isotropic, transversely isotropic, biaxially anisotropic and general anisotropic conductivity are modeled and compared to investigate the effects of anisotropy on the tool responses. For a tool in a general anisotropic formation, all the H components are nonzero. For a tool in a vertical well in transversely isotropic and biaxially anisotropic formations, only the diagonal components of H are nonzero. For a tool located in a deviated well, the effects of tool deviation and electric anisotropy are coupled. The diagonal components are more sensitive to the electric anisotropy than the off-diagonal components, and the off-diagonal ones can clearly indicate bed boundaries.

Geometrical spreading in a layered transversely isotropic medium with vertical symmetry axis

Geophysics ◽  
2003 ◽  
Vol 68 (6) ◽  
pp. 2082-2091 ◽  
Author(s):  
Bjørn Ursin ◽  
Ketil Hokstad
Keyword(s):  
Ray Tracing ◽  
Seismic Data ◽  
Symmetry Axis ◽  
Transversely Isotropic ◽  
Analytical Approximation ◽  
Geometrical Spreading ◽  
S Wave ◽  
Weak Anisotropy ◽  
P Waves ◽  
Amplitude Versus

Compensation for geometrical spreading is important in prestack Kirchhoff migration and in amplitude versus offset/amplitude versus angle (AVO/AVA) analysis of seismic data. We present equations for the relative geometrical spreading of reflected and transmitted P‐ and S‐wave in horizontally layered transversely isotropic media with vertical symmetry axis (VTI). We show that relatively simple expressions are obtained when the geometrical spreading is expressed in terms of group velocities. In weakly anisotropic media, we obtain simple expressions also in terms of phase velocities. Also, we derive analytical equations for geometrical spreading based on the nonhyperbolic traveltime formula of Tsvankin and Thomsen, such that the geometrical spreading can be expressed in terms of the parameters used in time processing of seismic data. Comparison with numerical ray tracing demonstrates that the weak anisotropy approximation to geometrical spreading is accurate for P‐waves. It is less accurate for SV‐waves, but has qualitatively the correct form. For P waves, the nonhyperbolic equation for geometrical spreading compares favorably with ray‐tracing results for offset‐depth ratios less than five. For SV‐waves, the analytical approximation is accurate only at small offsets, and breaks down at offset‐depth ratios less than unity. The numerical results are in agreement with the range of validity for the nonhyperbolic traveltime equations.

Reflection moveout approximations for P-waves in a moderately anisotropic homogeneous tilted transverse isotropy layer

Geophysics ◽  
2017 ◽  
Vol 82 (5) ◽  
pp. C175-C185 ◽  
Author(s):  
Ivan Pšenčík ◽  
Véronique Farra
Keyword(s):  
Transverse Isotropy ◽  
Transversely Isotropic ◽  
P Wave ◽  
Weak Anisotropy ◽  
P Waves ◽  
3D Space ◽  
Axis Of Symmetry ◽  
Relative Errors ◽  
Symmetry Tests ◽  
Ti Media

We have developed approximate nonhyperbolic P-wave moveout formulas applicable to weakly or moderately anisotropic media of arbitrary anisotropy symmetry and orientation. Instead of the commonly used Taylor expansion of the square of the reflection traveltime in terms of the square of the offset, we expand the square of the reflection traveltime in terms of weak-anisotropy (WA) parameters. No acoustic approximation is used. We specify the formulas designed for anisotropy of arbitrary symmetry for the transversely isotropic (TI) media with the axis of symmetry oriented arbitrarily in the 3D space. Resulting formulas depend on three P-wave WA parameters specifying the TI symmetry and two angles specifying the orientation of the axis of symmetry. Tests of the accuracy of the more accurate of the approximate formulas indicate that maximum relative errors do not exceed 0.3% or 2.5% for weak or moderate P-wave anisotropy, respectively.

Normal moveout from dipping reflectors in anisotropic media

Geophysics ◽  
1995 ◽  
Vol 60 (1) ◽  
pp. 268-284 ◽  
Author(s):  
Ilya Tsvankin
Keyword(s):  
Symmetry Axis ◽  
Transverse Isotropy ◽  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Weak Anisotropy ◽  
Anisotropic Models ◽  
Transversely Isotropic Media ◽  
Wide Range ◽  
Isotropic Media ◽  
The Difference

Description of reflection moveout from dipping interfaces is important in developing seismic processing methods for anisotropic media, as well as in the inversion of reflection data. Here, I present a concise analytic expression for normal‐moveout (NMO) velocities valid for a wide range of homogeneous anisotropic models including transverse isotropy with a tilted in‐plane symmetry axis and symmetry planes in orthorhombic media. In transversely isotropic media, NMO velocity for quasi‐P‐waves may deviate substantially from the isotropic cosine‐of‐dip dependence used in conventional constant‐velocity dip‐moveout (DMO) algorithms. However, numerical studies of NMO velocities have revealed no apparent correlation between the conventional measures of anisotropy and errors in the cosine‐of‐dip DMO correction (“DMO errors”). The analytic treatment developed here shows that for transverse isotropy with a vertical symmetry axis, the magnitude of DMO errors is dependent primarily on the difference between Thomsen parameters ε and δ. For the most common case, ε − δ > 0, the cosine‐of‐dip–corrected moveout velocity remains significantly larger than the moveout velocity for a horizontal reflector. DMO errors at a dip of 45 degrees may exceed 20–25 percent, even for weak anisotropy. By comparing analytically derived NMO velocities with moveout velocities calculated on finite spreads, I analyze anisotropy‐induced deviations from hyperbolic moveout for dipping reflectors. For transversely isotropic media with a vertical velocity gradient and typical (positive) values of the difference ε − δ, inhomogeneity tends to reduce (sometimes significantly) the influence of anisotropy on the dip dependence of moveout velocity.

scholarly journals Seismic AVOA Inversion for Weak Anisotropy Parameters and Fracture Density in a Monoclinic Medium

2020 ◽  
Vol 10 (15) ◽  
pp. 5136 ◽  
Author(s):  
Zijian Ge ◽  
Shulin Pan ◽  
Jingye Li
Keyword(s):  
Fractured Rock ◽  
Synthetic Data ◽  
Transversely Isotropic ◽  
P Wave ◽  
Inversion Method ◽  
Fracture Density ◽  
Weak Anisotropy ◽  
Amplitude Versus Offset ◽  
Porosity And Permeability ◽  
Transverse Isotropic

In shale gas development, fracture density is an important lithologic parameter to properly characterize reservoir reconstruction, establish a fracturing scheme, and calculate porosity and permeability. The traditional methods usually assume that the fracture reservoir is one set of aligned vertical fractures, embedded in an isotropic background, and estimate some alternative parameters associated with fracture density. Thus, the low accuracy caused by this simplified model, and the intrinsic errors caused by the indirect substitution, affect the estimation of fracture density. In this paper, the fractured rock of monoclinic symmetry assumes two non-orthogonal vertical fracture sets, embedded in a transversely isotropic background. Firstly, assuming that the fracture radius, width, and orientation are known, a new form of P-wave reflection coefficient, in terms of weak anisotropy (WA) parameters and fracture density, was obtained by substituting the stiffness coefficients of vertical transverse isotropic (VTI) background, normal, and tangential fracture compliances. Then, a linear amplitude versus offset and azimuth (AVOA) inversion method, of WA parameters and fracture density, was constructed by using Bayesian theory. Tests on synthetic data showed that WA parameters, and fracture density, are stably estimated in the case of seismic data containing a moderate noise, which can provide a reliable tool in fracture prediction.

scholarly journals Assessment of 2-D resistivity structures using 1-D inversion

Geophysics ◽  
1991 ◽  
Vol 56 (6) ◽  
pp. 874-883 ◽  
Author(s):  
L. P. Beard ◽  
F. D. Morgan
Keyword(s):  
Cross Section ◽  
Cross Sections ◽  
Apparent Resistivity ◽  
Inversion Process ◽  
Square Prism ◽  
Contour Lines ◽  
The Cross ◽  
Schlumberger Array ◽  
Wenner Array

Schlumberger and Wenner array resistivity soundings over 2-D resistivity structures are interpreted using apparent resistivity pseudosections and cross‐sections constructed from 1-D inversions in order to determine the effectiveness of 1-D interpretations over such structures. Cross‐sections contoured from resistivities of inverted “layers” show distinct differences from the apparent resistivity pseudosections and may be used as interpretational aids. Contour lines in the cross‐sections locate the horizontal interfaces of the 2-D structures quite well. The vertically oriented segments of the cross‐section contours are relatively undistorted in the inversion process and are similar to the vertically oriented portions of contours in the apparent resistivity pseudosection. A simple, empirically determined formula is used to separate the sections into resistive and conductive zones and helps to define the geometry of the anomaly. In order to apply the formula, it is necessary to know whether the target is a relative conductor or a relative resistor. Except for the case of a square prism, the Schlumberger array appears to hold advantages over the Wenner in qualitatively assessing an anomaly. The primary drawback of the Wenner array is that its expanding potential electrodes create false anomalous zones and complicate interpretation. As might be expected, structures with long horizontal interfaces, i.e. those more nearly 1-D, yield the most accurate interpretations.

Quartic moveout coefficient: 3D description and application to tilted TI media

Geophysics ◽  
2003 ◽  
Vol 68 (5) ◽  
pp. 1600-1610 ◽  
Author(s):  
Andres Pech ◽  
Ilya Tsvankin ◽  
Vladimir Grechka
Keyword(s):  
Heterogeneous Media ◽  
Symmetry Axis ◽  
High Sensitivity ◽  
Transversely Isotropic ◽  
Seismic Inversion ◽  
P Wave ◽  
Individual Values ◽  
Weak Anisotropy ◽  
Essential Information ◽  
Ti Media

Nonhyperbolic (long‐spread) moveout provides essential information for a number of seismic inversion/processing applications, particularly for parameter estimation in anisotropic media. Here, we present an analytic expression for the quartic moveout coefficient A4 that controls the magnitude of nonhyperbolic moveout of pure (nonconverted) modes. Our result takes into account reflection‐point dispersal on irregular interfaces and is valid for arbitrarily anisotropic, heterogeneous media. All quantities needed to compute A4 can be evaluated during the tracing of the zero‐offset ray, so long‐spread moveout can be modeled without time‐consuming multioffset, multiazimuth ray tracing. The general equation for the quartic coefficient is then used to study azimuthally varying nonhyperbolic moveout of P‐waves in a dipping transversely isotropic (TI) layer with an arbitrary tilt ν of the symmetry axis. Assuming that the symmetry axis is confined to the dip plane, we employed the weak‐anisotropy approximation to analyze the dependence of A4 on the anisotropic parameters. The linearized expression for A4 is proportional to the anellipticity coefficient η ≈ ε − δ and does not depend on the individual values of the Thomsen parameters. Typically, the magnitude of nonhyperbolic moveout in tilted TI media above a dipping reflector is highest near the reflector strike, whereas deviations from hyperbolic moveout on the dip line are substantial only for mild dips. The azimuthal variation of the quartic coefficient is governed by the tilt ν and reflector dip φ and has a much more complicated character than the NMO–velocity ellipse. For example, if the symmetry axis is vertical (VTI media, ν = 0) and the dip φ < 30°, A4 goes to zero on two lines with different azimuths where it changes sign. If the symmetry axis is orthogonal to the reflector (this model is typical for thrust‐and‐fold belts), the strike‐line quartic coefficient is defined by the well‐known expression for a horizontal VTI layer (i.e., it is independent of dip), while the dip‐line A4 is proportional to cos4 φ and rapidly decreases with dip. The high sensitivity of the quartic moveout coefficient to the parameter η and the tilt of the symmetry axis can be exploited in the inversion of wide‐azimuth, long‐spread P‐wave data for the parameters of TI media.

Reflection and transmission responses in a layered transversely isotropic medium with horizontal symmetry axis

Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. C143-C157 ◽  
Author(s):  
Song Jin ◽  
Alexey Stovas
Keyword(s):  
Good Accuracy ◽  
Symmetry Axis ◽  
Local Property ◽  
Transversely Isotropic ◽  
Weak Anisotropy ◽  
Reflection And Transmission ◽  
Numerical Tests ◽  
Transversely Isotropic Media ◽  
Horizontal Symmetry ◽  
Isotropic Media

Seismic wave reflection and transmission (R/T) responses characterize the subsurface local property, and the widely spread anisotropy has considerable influences even at small incident angles. We have considered layered transversely isotropic media with horizontal symmetry axes (HTI), and the symmetry axes were not restricted to be aligned. With the assumption of weak contrast across the interface, linear approximations for R/T coefficients normalized by vertical energy flux are derived based on a simple layered HTI model. We also obtain the approximation with the isotropic background medium under an additional weak anisotropy assumption. Numerical tests illustrate the good accuracy of the approximations compared with the exact results.

Sign in / Sign up

Export Citation Format

Share Document