Estimation of interval anisotropic attenuation from reflection data

Geophysics ◽  
2009 ◽  
Vol 74 (6) ◽  
pp. A69-A74 ◽  
Author(s):  
Jyoti Behura ◽  
Ilya Tsvankin
Keyword(s):  
Attenuation Coefficient ◽  
Reservoir Characterization ◽  
Symmetry Plane ◽  
Transversely Isotropic ◽  
Spectral Ratio ◽  
P Wave ◽  
Ratio Method ◽  
Reflection Seismic ◽  
Reflection Data ◽  
Target Layer

Knowledge of interval attenuation can be highly beneficial in reservoir characterization and lithology discrimination. We combine the spectral-ratio method with velocity-independent layer stripping to develop a technique for the estimation of the interval attenuation coefficient from reflection seismic data. The layer-stripping procedure is based on identifying the reflections from the top and bottom of the target layer that share the same ray segments in the overburden. The algorithm is designed for heterogeneous, arbitrarily anisotropic target layers, but the overburden is assumed to be laterally homogeneous with a horizontal symmetry plane. Although no velocity information about the overburden is needed, interpretation of the computed anisotropic attenuation coefficient involves the phase angle in the target layer. Tests on synthetic P-wave data from layered transversely isotropic and orthorhombic media confirm the high accuracy of 2D and 3D versions of the algorithm. We also demonstrate that the interval attenuation estimates are independent of the inhomogeneity angle of the incident and reflected waves.

Estimation of shear-wave interval attenuation from mode-converted data

Geophysics ◽  
2011 ◽  
Vol 76 (6) ◽  
pp. D11-D19 ◽  
Author(s):  
Bharath Shekar ◽  
Ilya Tsvankin
Keyword(s):  
Shear Wave ◽  
Attenuation Coefficient ◽  
Reservoir Characterization ◽  
Wave Attenuation ◽  
Pure Shear ◽  
Synthetic Data ◽  
Spectral Ratio ◽  
Ratio Method ◽  
S Wave ◽  
Target Layer

Interval attenuation measurements provide valuable information for reservoir characterization and lithology discrimination. We extend the attenuation layer-stripping method of Behura and Tsvankin to mode-converted (PS) waves with the goal of estimating the S-wave interval attenuation coefficient. By identifying PP and PS events with shared ray segments and applying the [Formula: see text] method, we first perform kinematic construction of pure shear (SS) events in the target layer and overburden. Then, the modified spectral-ratio method is used to compute the effective shear-wave attenuation coefficient for the target reflection. Finally, application of the dynamic version of velocity-independent layer stripping to the constructed SS reflections yields the interval S-wave attenuation coefficient in the target layer. The attenuation coefficient estimated for a range of source-receiver offsets can be inverted for the interval attenuation parameters. The method is tested on multicomponent synthetic data generated with the anisotropic reflectivity method for layered VTI (transversely isotropic with a vertical symmetry axis) and orthorhombic media.

Physical modeling and analysis of P-wave attenuation anisotropy in transversely isotropic media

Geophysics ◽  
2007 ◽  
Vol 72 (1) ◽  
pp. D1-D7 ◽  
Author(s):  
Yaping Zhu ◽  
Ilya Tsvankin ◽  
Pawan Dewangan ◽  
Kasper van Wijk
Keyword(s):  
Attenuation Coefficient ◽  
Symmetry Axis ◽  
Wave Attenuation ◽  
Transversely Isotropic ◽  
P Wave ◽  
Velocity Variation ◽  
Ratio Method ◽  
Fracture Detection ◽  
Wide Range ◽  
Attenuation Anisotropy

Anisotropic attenuation can provide sensitive attributes for fracture detection and lithology discrimination. This paper analyzes measurements of the P-wave attenuation coefficient in a transversely isotropic sample made of phenolic material. Using the spectral-ratio method, we estimate the group (effective) attenuation coefficient of P-waves transmitted through the sample for a wide range of propagation angles (from [Formula: see text] to [Formula: see text]) with the symmetry axis. Correction for the difference between the group and phase angles and for the angular velocity variation help us to obtain the normalized phase attenuation coefficient [Formula: see text] governed by the Thomsen-style attenuation-anisotropy parameters [Formula: see text] and [Formula: see text]. Whereas the symmetry axis of the angle-dependent coefficient [Formula: see text] practically coincides with that of the velocity function, the magnitude of the attenuation anisotropy far exceeds that of the velocity anisotropy. The quality factor [Formula: see text] increases more than tenfold from the symmetry axis (slow direction) to the isotropy plane (fast direction). Inversion of the coefficient [Formula: see text] using the Christoffel equation yields large negative values of the parameters [Formula: see text] and [Formula: see text]. The robustness of our results critically depends on several factors, such as the availability of an accurate anisotropic velocity model and adequacy of the homogeneous concept of wave propagation, as well as the choice of the frequency band. The methodology discussed here can be extended to field measurements of anisotropic attenuation needed for AVO (amplitude-variation-with-offset) analysis, amplitude-preserving migration, and seismic fracture detection.

Estimation of interval anisotropy parameters using velocity-independent layer stripping

Geophysics ◽  
2009 ◽  
Vol 74 (5) ◽  
pp. WB117-WB127 ◽  
Author(s):  
Xiaoxiang Wang ◽  
Ilya Tsvankin
Keyword(s):  
Symmetry Plane ◽  
Synthetic Data ◽  
Single Layer ◽  
Transversely Isotropic ◽  
P Wave ◽  
Correlated Noise ◽  
Fractured Reservoir ◽  
Processing Parameter ◽  
Reflection Data ◽  
Error Amplification

Moveout analysis of long-spread P-wave data is widely used to estimate the key time-processing parameter [Formula: see text] in layered transversely isotropic media with a vertical symmetry axis (VTI). Inversion for interval [Formula: see text] values, however, suffers from instability caused by the trade-off between the effective moveout parameters and by subsequent error amplification during Dix-type layer stripping. We propose an alternative approach to nonhyperbolic moveout inversion based on the velocity-independent layer-stripping (VILS) method of Dewangan and Tsvankin. Also, we develop the 3D version of VILS and apply it to interval parameter estimation in orthorhombic media using wide-azimuth, long-spread data. If the overburden is laterally homogeneous and has a horizontal symmetry plane, VILS produces the exact interval traveltime-offset function in the target layer without knowledgeof the velocity field. Hence, Dix-type differentiation of moveout parameters used in existing techniques is replaced by the much more stable layer stripping of reflection traveltimes. The interval traveltimes are then inverted for the moveout parameters using the single-layer nonhyperbolic moveout equation. The superior accuracy and stability of the algorithm are illustrated on ray-traced synthetic data for typical VTI and orthorhombic models. Even small correlated noise in reflection traveltimes causes substantial distortions in the interval [Formula: see text] values computed by conventional Dix-type differentiation. In contrast, the output of VILS is insensitive to mild correlated traveltime errors. The algorithm is also tested on wide-azimuth P-wave reflection data recorded above a fractured reservoir at Rulison field in Colorado. The interval moveout parameters estimated by VILS in the shale layer above the reservoir are more plausible and less influenced by noise than those obtained by the Dix-type method.

Anisotropic geometrical-spreading correction for wide-azimuth P-wave reflections

Geophysics ◽  
2006 ◽  
Vol 71 (5) ◽  
pp. D161-D170 ◽  
Author(s):  
Xiaoxia Xu ◽  
Ilya Tsvankin
Keyword(s):  
Transversely Isotropic ◽  
Azimuthal Anisotropy ◽  
P Wave ◽  
Fractured Reservoir ◽  
Geometrical Spreading ◽  
Spreading Factor ◽  
Velocity Anisotropy ◽  
Reflection Data ◽  
Fitting Procedure ◽  
Avo Attributes

Compensation for geometrical spreading along a raypath is one of the key steps in AVO (amplitude-variation-with-offset) analysis, in particular, for wide-azimuth surveys. Here, we propose an efficient methodology to correct long-spread, wide-azimuth reflection data for geometrical spreading in stratified azimuthally anisotropic media. The P-wave geometrical-spreading factor is expressed through the reflection traveltime described by a nonhyperbolic moveout equation that has the same form as in VTI (transversely isotropic with a vertical symmetry axis) media. The adapted VTI equation is parameterized by the normal-moveout (NMO) ellipse and the azimuthally varying anellipticity parameter [Formula: see text]. To estimate the moveout parameters, we apply a 3D nonhyperbolic semblance algorithm of Vasconcelos and Tsvankin that operates simultaneously with traces at all offsets andazimuths. The estimated moveout parameters are used as the input in our geometrical-spreading computation. Numerical tests for models composed of orthorhombic layers with strong, depth-varying velocity anisotropy confirm the high accuracy of our travetime-fitting procedure and, therefore, of the geometrical-spreading correction. Because our algorithm is based entirely on the kinematics of reflection arrivals, it can be incorporated readily into the processing flow of azimuthal AVO analysis. In combination with the nonhyperbolic moveout inversion, we apply our method to wide-azimuth P-wave data collected at the Weyburn field in Canada. The geometrical-spreading factor for the reflection from the top of the fractured reservoir is clearly influenced by azimuthal anisotropy in the overburden, which should cause distortions in the azimuthal AVO attributes. This case study confirms that the azimuthal variation of the geometrical-spreading factor often is comparable to or exceeds that of the reflection coefficient.

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.

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.

scholarly journals Microseismic Wave Measurements to Detect Landslides in Bengkulu Shore with Attenuation Coefficient and Shear Strain Indicator

2016 ◽  
Vol 1 (1) ◽  
Author(s):  
Muhammad Farid
Keyword(s):  
Shear Strain ◽  
Attenuation Coefficient ◽  
Peak Ground Acceleration ◽  
Amplification Factor ◽  
Spectral Ratio ◽  
Coastal Areas ◽  
Ratio Method ◽  
Ground Acceleration ◽  
Wave Measurements ◽  
Microseismic Data

<p>It has been detected that the condition of landslides that occurred in Bengkulu Shore can change the position of the shoreline. This research aimed to: (1) calculate of shear strain (γ) and attenuation coefficient (ά) value  based on microseismic data in coastal areas that experienced landslides; (2) determine the correlation between levels of landslides with  shear strain  and attenuation coefficient value (3) determine the correlation between the shear strain and attenuation coefficient value. Microseismic data were processed and analyzed quantitatively using the Horizontal to Vertical Spectral Ratio method (HVSR) to obtain the ground vibrations resonance frequency (<em>f<sub>o</sub></em>) and amplification factor (<em>A</em>). Shear strain value was calculated from the of <em>f<sub>o</sub></em>, <em>A</em> and Peak Ground Acceleration (<em>α<sub>max</sub></em>) value. Peak Ground Acceleration value was calculated based on 100-year period of recorded earthquake data.  Attenuation coefficient was calculated based on the equation (2). The results of study showed that the value of shear strain in the coastal areas varied from 1.0 × 10<sup>-4</sup> to 3.6 × 10<sup>-3</sup>,  in accordance with the conditions of landslides. The attenuation coefficient value varied from 0.005 to 0.020.  Level of landslides that occurred varied from moderate, to very severe. There was a tendency that the more severe the landslide level,  the greater the shear strain and attenuation coefficient value were.</p>

Estimation of Q from surface seismic reflection data

Geophysics ◽  
1998 ◽  
Vol 63 (6) ◽  
pp. 2120-2128 ◽  
Author(s):  
Rahul Dasgupta ◽  
Roger A. Clark
Keyword(s):  
Attenuation Factor ◽  
Spectral Ratio ◽  
Ratio Method ◽  
Amplitude Analysis ◽  
Gas Reservoir ◽  
Seismic Reflection Data ◽  
Reflection Data ◽  
Anelastic Attenuation ◽  
Zero Offset ◽  
Nmo Stretch

Reliable estimates of the anelastic attenuation factor, Q, are desirable for improved resolution through inverse Q deconvolution and to facilitate amplitude analysis. Q is a useful petrophysical parameter itself, yet Q is rarely measured. Estimates must currently be made from borehole seismology. This paper presents a simple technique for determining Q from conventional surface seismic common midpoint (CMP) gathers. It is essentially the classic spectral ratio method applied on a trace‐by‐trace basis to a designatured and NMO stretch‐corrected CMP gather. The variation of apparent Q versus offset (QVO) is extrapolated to give a zero‐offset Q estimate. Studies on synthetics suggest that, for reasonable data quality (S/N ratios better than 3:1, shallow (<5°) dips, and stacking velocity accuracy <5%), source‐to‐reflector average Q is recoverable to within some 3% and Q for a specific interval (depending on its two‐way time thickness and depth) is recoverable to 15–20%. Three case studies are reported. First, Q versus offset and vertical seismic profiling (VSP) Q estimates for a southern North Sea line were in close agreement, validating the method. For Chalk, Mushelkalk‐Keuper, and Bunter‐Zechstein, Q was estimated as 130 ± 15, 47 ± 8, and 156 ± 18, respectively. Next, two alternative lithological interpretations of a structure seen in a frontier area were discriminated between when Q estimates of 680 to 820 were obtained (compared to some 130–170 in the overlying units), favoring a metamorphic/crystalline lithology rather than (prospective) sediments. This was later confirmed by drilling. Third, a profile of Q estimates along a 200-ms-thick interval, known to include a gas reservoir, showed a clear and systematic reduction in Q to a low of 50 ± 11, coincident with the maximum reservoir thickness, from some 90–105 outside the reservoir. Q for the reservoir interval itself was estimated at 17 ± 7.

An effective anisotropy parameter in transversely isotropic media

Geophysics ◽  
1987 ◽  
Vol 52 (12) ◽  
pp. 1654-1664 ◽  
Author(s):  
N. C. Banik
Keyword(s):  
Transverse Isotropy ◽  
Anisotropy Parameter ◽  
Transversely Isotropic ◽  
P Wave ◽  
Weak Anisotropy ◽  
Reflection Seismic ◽  
The North ◽  
Reflection Amplitude ◽  
Transversely Isotropic Media ◽  
Isotropic Media

An interesting physical meaning is presented for the anisotropy parameter δ, previously introduced by Thomsen to describe weak anisotropy in transversely isotropic media. Roughly, δ is the difference between the P-wave and SV-wave anisotropies of the medium. The observed systematic depth errors in the North Sea are reexamined in view of the new interpretation of the moveout velocity through δ. The changes in δ at an interface adequately describe the effects of transverse isotropy on the P-wave reflection amplitude, The reflection coefficient expression is linearized in terms of changes in elastic parameters. The linearized expression clearly shows that it is the variation of δ at the interface that gives the anisotropic effects at small incidence angles. Thus, δ effectively describes both the moveout velocity and the reflection amplitude variation, two very important pieces of information in reflection seismic prospecting, in the presence of transverse isotropy.

Accurate interval Q-factor estimation from VSP data

Geophysics ◽  
2012 ◽  
Vol 77 (3) ◽  
pp. WA149-WA156 ◽  
Author(s):  
E. Blias
Keyword(s):  
Real Data ◽  
Spectral Ratio ◽  
Ratio Method ◽  
Seismic Reflection Data ◽  
Reflection Data ◽  
Seismic Quality Factor ◽  
Vsp Data ◽  
Considerable Impact ◽  
The Difference ◽  
Factor Estimation

Inelastic attenuation, quantified by [Formula: see text], the seismic quality factor, has considerable impact on surface seismic reflection data. A new method for interval [Formula: see text]-factor estimation using near-offset VSP data was based on an objective function minimization measuring the difference between cumulative [Formula: see text] estimates and those calculated through interval [Formula: see text]. To calculate interval [Formula: see text], we used all receiver pairs that provided reasonable [Formula: see text] values. To estimate [Formula: see text] between two receiver levels, we used the equation that links amplitudes at different levels and could provide more accurate [Formula: see text] values than the spectral-ratio method. To improve interval [Formula: see text] estimates, which rely on traveltimes, we used a high-accuracy approach in the frequency domain to determine time shifts. Application of this method to real data demonstrated reasonable correspondence between [Formula: see text] estimates and log data.

Sign in / Sign up

Export Citation Format

Share Document