SOME EFFECTS OF FORMATION ANISOTROPY ON RESISTIVITY MEASUREMENTS IN BOREHOLES

Geophysics ◽  
1958 ◽  
Vol 23 (4) ◽  
pp. 770-794 ◽  
Author(s):  
K. S. Kunz ◽  
J. H. Moran

It is shown that a wide class of potential problems involving anisotropic media can be transformed into equivalent problems involving only isotropic media. By means of such transformations it is possible, in a large number of cases, to determine the apparent resistivities which would be observed in anisotropic formations, using electrode‐type resistivity logging devices. Discussion is given of an infinite, anisotropic medium with and without borehole, of two semi‐infinite anisotropic beds (without borehole), and of a thin isotropic bed bounded by anisotropic adjacent formations (without borehole). An interpretation chart for the normal device is presented for thick, non‐invaded, anisotropic beds penetrated by a borehole.

Geophysics ◽  
2001 ◽  
Vol 66 (1) ◽  
pp. 308-319 ◽  
Author(s):  
Ivan Pšenčík ◽  
Joe A. Dellinger

In inhomogeneous isotropic regions, S-waves can be modeled using the ray method for isotropic media. In inhomogeneous strongly anisotropic regions, the independently propagating qS1- and qS2-waves can similarly be modeled using the ray method for anisotropic media. The latter method does not work properly in inhomogenous weakly anisotropic regions, however, where the split qS-waves couple. The zeroth‐order approximation of the quasi‐isotropic (QI) approach was designed for just such inhomogeneous weakly anisotropic media, for which neither the ray method for isotropic nor anisotropic media applies. We test the ranges of validity of these three methods using two simple synthetic models. Our results show that the QI approach more than spans the gap between the ray methods: it can be used in isotropic regions (where it reduces to the ray method for isotropic media), in regions of weak anisotropy (where the ray method for anisotropic media does not work properly), and even in regions of moderately strong anisotropy (in which the qS-waves decouple and thus could be modeled using the ray method for anisotropic media). A modeling program that switches between these three methods as necessary should be valid for arbitrary‐strength anisotropy.


Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. C159-C170 ◽  
Author(s):  
Yuriy Ivanov ◽  
Alexey Stovas

Based on the rotation of a slowness surface in anisotropic media, we have derived a set of mapping operators that establishes a point-to-point correspondence for the traveltime and relative-geometric-spreading surfaces between these calculated in nonrotated and rotated media. The mapping approach allows one to efficiently obtain the aforementioned surfaces in a rotated anisotropic medium from precomputed surfaces in the nonrotated medium. The process consists of two steps: calculation of a necessary kinematic attribute in a nonrotated, e.g., orthorhombic (ORT), medium, and subsequent mapping of the obtained values to a transformed, e.g., rotated ORT, medium. The operators we obtained are applicable to anisotropic media of any type; they are 3D and are expressed through a general form of the transformation matrix. The mapping equations can be used to develop moveout and relative-geometric-spreading approximations in rotated anisotropic media from existing approximations in nonrotated media. Although our operators are derived in case of a homogeneous medium and for a one-way propagation only, we discuss their extension to vertically heterogeneous media and to reflected (and converted) waves.


Geophysics ◽  
1995 ◽  
Vol 60 (1) ◽  
pp. 268-284 ◽  
Author(s):  
Ilya Tsvankin

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.


Geophysics ◽  
2021 ◽  
pp. 1-74
Author(s):  
Iury Araújo ◽  
Murillo Nascimento ◽  
Jesse Costa ◽  
Alan Souza ◽  
Jörg Schleicher

We present a procedure to derive low-rank evolution operators in the mixed space-wavenumber domain for modeling the qP Born-scattered wavefield at perturbations of an anisotropic medium under the pseudo-acoustic approximation. To approximate the full wavefield, this scattered field is then added to the reference wavefield obtained with the corresponding low-rank evolution operator in the background medium. Being built upon a Hamiltonian formulation using the dispersion relation for qP waves, this procedure avoids pseudo-S-wave artifacts and provides a unified approach for linearizing anisotropic pseudo-acoustic evolution operators. Therefore it is immediately applicable to any arbitrary class of anisotropy. As an additional asset, the scattering operators explicitly contain the sensitivity kernels of the Born-scattered wavefield with respect to the anisotropic medium parameters. This enables direct access to important information like its offset dependence or directional characteristics as a function of the individual parameter perturbations. For our numerical tests, we specify the operators for a mildly anisotropic tilted transversly isotropic (TTI) medium. We validate our implementation in a simple model with weak contrasts and simulate reflection data in the BP TTI model to show that the procedure works in a more realistic scenario. The Born-scattering results indicate that our procedure is applicable to strongly heterogeneous anisotropic media. Moreover, we use the analytical capabilities of the kernels by means of sensitivity tests to demonstrate that using two different medium parameterizations leads to different results. The mathematical formulation of the method is such that it allows for an immediate application to least-squares migration in pseudo-acoustic anisotropic media.


Geophysics ◽  
1996 ◽  
Vol 61 (4) ◽  
pp. 947-963 ◽  
Author(s):  
Tariq Alkhalifah

Nearly all dip‐moveout correction (DMO) implementations to date assume isotropic homogeneous media. Usually, this has been acceptable considering the tremendous cost savings of homogeneous isotropic DMO and considering the difficulty of obtaining the anisotropy parameters required for effective implementation. In the presence of typical anisotropy, however, ignoring the anisotropy can yield inadequate results. Since anisotropy may introduce large deviations from hyperbolic moveout, accurate transformation to zero‐offset in anisotropic media should address such nonhyperbolic moveout behavior of reflections. Artley and Hale’s v(z) ray‐tracing‐based DMO, developed for isotropic media, provides an attractive approach to treating such problems. By using a ray‐tracing procedure crafted for anisotropic media, I modify some aspects of their DMO so that it can work for v(z) anisotropic media. DMO impulse responses in typical transversely isotropic (TI) models (such as those associated with shales) deviate substantially from the familiar elliptical shape associated with responses in homogeneous isotropic media (to the extent that triplications arise even where the medium is homogeneous). Such deviations can exceed those caused by vertical inhomogeneity, thus emphasizing the importance of taking anisotropy into account in DMO processing. For isotropic or elliptically anisotropic media, the impulse response is an ellipse; but as the key anisotropy parameter η varies, the shape of the response differs substantially from elliptical. For typical η > 0, the impulse response in TI media tends to broaden compared to the response in an isotropic homogeneous medium, a behavior opposite to that encountered in typical v(z) isotropic media, where the response tends to be squeezed. Furthermore, the amplitude distribution along the DMO operator differs significantly from that for isotropic media. Application of this anisotropic DMO to data from offshore Africa resulted in a considerably better alignment of reflections from horizontal and dipping reflectors in common‐midpoint gather than that obtained using an isotropic DMO. Even the presence of vertical inhomogeneity in this medium could not eliminate the importance of considering the shale‐induced anisotropy.


Geophysics ◽  
2013 ◽  
Vol 78 (4) ◽  
pp. T99-T112 ◽  
Author(s):  
Siyi Shao ◽  
Laura J. Pyrak-Nolte

The detection of fractures in an anisotropic medium is complicated by discrete modes that are guided or confined by fractures such as fracture interface waves. Fracture interface waves are generalized coupled Rayleigh waves whose existence and velocity in isotropic media depend on the stiffness of the fracture, frequency of the source, and shear-wave polarization. We derived the analytic solution for fracture interface waves in an orthorhombic medium and found that the existence and velocity of interface waves in anisotropic media are also affected by the orientation of a fracture relative to the layering. Laboratory measurements of fracture interface waves using ultrasonic transducers (central frequency [Formula: see text] MHz) on garolite specimens confirmed that the presence of fracture interface waves can mask the textural shear-wave anisotropy of waves propagating parallel to the layering. At low stresses, a layered medium appears almost isotropic when a fracture is oriented perpendicular to the layering, and conversely, a layered medium exhibits stronger anisotropy than the matrix for a fracture oriented parallel to the layering. The matrix shear-wave anisotropy is recovered when sufficient stress is applied to close a fracture. The theory and experimental results demonstrated that the interpretation of the presence of fractures in anisotropic material can be unambiguously interpreted if measurements are made as a function of stress, which eliminates many fractured-generated discrete modes such as fracture interface waves.


Geophysics ◽  
1994 ◽  
Vol 59 (4) ◽  
pp. 591-596 ◽  
Author(s):  
Suhas Phadke ◽  
S. Kapotas ◽  
N. Dai ◽  
Ernest R. Kanasewich

Wave propagation in transversely isotropic media is governed by the horizontal and vertical wave velocities. The quasi‐P(qP) wavefront is not an ellipse; therefore, the propagation cannot be described by the wave equation appropriate for elliptically anisotropic media. However, for a limited range of angles from the vertical, the dispersion relation for qP‐waves can be approximated by an ellipse. The horizontal velocity necessary for this approximation is different from the true horizontal velocity and depends upon the physical properties of the media. In the method described here, seismic data is migrated using a 45-degree wave equation for elliptically anisotropic media with the horizontal velocity determined by comparing the 45-degree elliptical dispersion relation and the quasi‐P‐dispersion relation. The method is demonstrated for some synthetic data sets.


Sign in / Sign up

Export Citation Format

Share Document