Role of the inhomogeneity angle in anisotropic attenuation analysis

Geophysics ◽  
2009 ◽  
Vol 74 (5) ◽  
pp. WB177-WB191 ◽  
Author(s):  
Jyoti Behura ◽  
Ilya Tsvankin

The inhomogeneity angle (the angle between the real and imaginary parts of the wave vector) is seldom taken into account in estimating attenuation coefficients from seismic data. Wave propagation through the subsurface, however, can result in relatively large inhomogeneity angles [Formula: see text], especially for models with significant attenuation contrasts across layer boundaries. Here we study the influence of the angle [Formula: see text] on phase and group attenuation in arbitrarily anisotropic media using the first-order perturbation theory verified by exact numerical modeling. Application of the spectral-ratio method to transmitted or reflected waves yields the normalized group attenuation coefficient [Formula: see text], which is responsible for amplitude decay along seismic rays. Our analytic solutions show that for a wide range of inhomogeneity angles, the coefficient [Formula: see text] is close to the normalized phase attenuation coefficient [Formula: see text] computed for [Formula: see text] [Formula: see text]. The coefficient[Formula: see text] can be inverted directly for the attenuation-anisotropy parameters, so no knowledge of the inhomogeneity angle is required for attenuation analysis of seismic data. This conclusion remains valid even for uncommonly high attenuation with the quality factor [Formula: see text] less than 10 and strong velocity and attenuation anisotropy. However, the relationship between group and phase attenuation coefficients becomes more complicated for relatively large inhomogeneity angles approaching so-called ‘‘forbidden directions.’’ We also demonstrate that the velocity function remains practically independent of attenuation for a wide range of small and moderate angles [Formula: see text]. In principle, estimation of the attenuation-anisotropy parameters from the coefficient [Formula: see text] requires computation of the phase angle, which depends on the anisotropic velocity field. For moderately anisotropic models, however, the difference between the phase and group directions should not significantly distort the results of attenuation analysis.

Geophysics ◽  
2007 ◽  
Vol 72 (1) ◽  
pp. D9-D19 ◽  
Author(s):  
Yaping Zhu ◽  
Ilya Tsvankin

Orthorhombic models are often used in the interpretation of azimuthally varying seismic signatures recorded over fractured reservoirs. Here, we develop an analytic framework for describing the attenuation coefficients in orthorhombic media with orthorhombic attenuation (i.e., the symmetry of both the real and imaginary parts of the stiffness tensor is identical) under the assumption of homogeneous wave propagation. The analogous form of the Christoffel equation in the symmetry planes of orthorhombic and VTI (transversely isotropic with a vertical symmetry axis) media helps to obtain the symmetry-plane attenuation coefficients by adapting the existing VTI equations. To take full advantage of this equivalence with transverse isotropy, we introduce a parameter set similar to the VTI attenuation-anisotropy parameters [Formula: see text], [Formula: see text], and [Formula: see text]. This notation, based on the same principle as Tsvankin’s velocity-anisotropy parameters for orthorhombic media, leads to concise linearized equations for thesymmetry-plane attenuation coefficients of all three modes (P, [Formula: see text], and [Formula: see text]).The attenuation-anisotropy parameters also allow us to simplify the P-wave attenuation coefficient [Formula: see text] outside the symmetry planes under the assumptions of small attenuation and weak velocity and attenuation anisotropy. The approximate coefficient [Formula: see text] has the same form as the linearized P-wave phase-velocity function, with the velocity parameters [Formula: see text] and [Formula: see text] replaced by the attenuation parameters [Formula: see text] and [Formula: see text]. The exact attenuation coefficient, however, also depends on the velocity-anisotropy parameters, while the body-wave velocities are almost unperturbed by the presence of attenuation. The reduction in the number of parameters responsible for the P-wave attenuation and the simple approximation for the coefficient [Formula: see text] provide a basis for inverting P-wave attenuation measurements from orthorhombic media. The attenuation processing must be preceded by anisotropic velocity analysis that can be performed (in the absence of pronounced velocity dispersion) using existing algorithms for nonattenuative media.


Geophysics ◽  
2007 ◽  
Vol 72 (1) ◽  
pp. D1-D7 ◽  
Author(s):  
Yaping Zhu ◽  
Ilya Tsvankin ◽  
Pawan Dewangan ◽  
Kasper van Wijk

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.


Geophysics ◽  
1990 ◽  
Vol 55 (8) ◽  
pp. 998-1014 ◽  
Author(s):  
T. Klimentos ◽  
C. McCann

Anelastic attenuation is the process by which rocks convert compressional waves into heat and thereby modify the amplitude and phase of the waves. Understanding the causes of compressional wave attenuation is important in the acquisition, processing, and interpretation of high‐resolution seismic data, and in deducing the physical properties of rocks from seismic data. We have measured the attenuation coefficients of compressional waves in 42 sandstones at a confining pressure of 40 MPa (equivalent to a depth of burial of about 1.5 km) in a frequency range from 0.5 to 1.5 MHz. The compressional wave measurements were made using a pulse‐echo method in which the sample (5 cm diameter, 1.8 cm to 3.5 cm long) was sandwiched between perspex (lucite) buffer rods inside the high‐pressure rig. The attenuation of the sample was estimated from the logarithmic spectral ratio of the signals (corrected for beam spreading) reflected from the top and base of the sample. The results show that for these samples, compressional wave attenuation (α, dB/cm) at 1 MHz and 40 MPa is related to clay content (C, percent) and porosity (ϕ, percent) by α=0.0315ϕ+0.241C−0.132 with a correlation coefficient of 0.88. The relationship between attenuation and permeability is less well defined: Those samples with permeabilities less than 50 md have high attenuation coefficients (generally greater than 1 dB/cm) while those with permeabilities greater than 50 md have low attenuation coefficients (generally less than 1 dB/cm) at 1 MHz at 40 MPa. These experimental data can be accounted for by modifications of the Biot theory and by consideration of the Sewell/Urick theory of compressional wave attenuation in porous, fluid‐saturated media.


2018 ◽  
Author(s):  
Marius Kriegerowski ◽  
Simone Cesca ◽  
Matthias Ohrnberger ◽  
Torsten Dahm ◽  
Frank Krüger

Abstract. We develop an amplitude spectral ratio method for event couples from clustered earthquakes to estimate seismic wave attenuation (Q−1) in the source volume. The method allows to study attenuation within the source region of earthquake swarms or aftershocks at depth, independent of wave path and attenuation between source region and surface station. We exploit the high frequency slope of phase spectra using multitaper spectral estimates. The method is tested using simulated full wavefield seismograms affected by recorded noise and finite source rupture. The synthetic tests verify the approach and show that solutions are independent of focal mechanisms, but also show that seismic noise may broaden the scatter of results. We apply the event couple spectral ratio method to North-West Bohemia, Czech Republic, a region characterized by the persistent occurrence of earthquake swarms in a confined source region at mid-crustal depth. Our method indicates a strong anomaly of high attenuation in the source region of the swarm with an averaged attenuation factor of Qp 


2016 ◽  
Vol 1 (1) ◽  
Author(s):  
Muhammad Farid

<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>


2018 ◽  
Vol 106 (4) ◽  
pp. 337-344 ◽  
Author(s):  
Moustafa A. Hilal ◽  
Mohamed F. Attallah

AbstractIn the present work, materials, namely mud, sand mud, ferruginous sandstone and sandstone with different densities are used to evaluate the effect of chemical composition and moisture content on the self-attenuation coefficient factor at γ-energy range from 59.5 to 1332.5 keV. The results revealed that the attenuation coefficient increases with increasing the moisture content until the material saturate with moisture. The average value of increasing linear attenuation coefficients based on increasing moisture content are 14.3%, 16.0%, 18.2%, 28.1% and 24.8% at γ-energies 59.5, 356.0, 661.7, 1173.4 and 1332.5 keV, respectively. Chemical composition of material affected on the values of attenuation, i.e. the elements with high density have high attenuation coefficient. Significant effect of self-attenuation correction factor was observed at low γ-energies up to 500 KeV. Application of the calculated correction of environmental sample with low radioactivity content has been carried out.


Solid Earth ◽  
2019 ◽  
Vol 10 (1) ◽  
pp. 317-328
Author(s):  
Marius Kriegerowski ◽  
Simone Cesca ◽  
Matthias Ohrnberger ◽  
Torsten Dahm ◽  
Frank Krüger

Abstract. We develop an amplitude spectral ratio method for event couples from clustered earthquakes to estimate seismic wave attenuation (Q−1) in the source volume. The method allows to study attenuation within the source region of earthquake swarms or aftershocks at depth, independent of wave path and attenuation between source region and surface station. We exploit the high-frequency slope of phase spectra using multitaper spectral estimates. The method is tested using simulated full wave-field seismograms affected by recorded noise and finite source rupture. The synthetic tests verify the approach and show that solutions are independent of focal mechanisms but also show that seismic noise may broaden the scatter of results. We apply the event couple spectral ratio method to northwest Bohemia, Czech Republic, a region characterized by the persistent occurrence of earthquake swarms in a confined source region at mid-crustal depth. Our method indicates a strong anomaly of high attenuation in the source region of the swarm with an averaged attenuation factor of Qp<100. The application to S phases fails due to scattered P-phase energy interfering with S phases. The Qp anomaly supports the common hypothesis of highly fractured and fluid saturated rocks in the source region of the swarms in northwest Bohemia. However, high temperatures in a small volume around the swarms cannot be excluded to explain our observations.


Geophysics ◽  
2011 ◽  
Vol 76 (6) ◽  
pp. D11-D19 ◽  
Author(s):  
Bharath Shekar ◽  
Ilya Tsvankin

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.


2016 ◽  
Vol 1 ◽  
Author(s):  
Muhammad Farid

<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>


Geophysics ◽  
2019 ◽  
Vol 84 (6) ◽  
pp. M25-M36 ◽  
Author(s):  
Mengqiang Pang ◽  
Jing Ba ◽  
José M. Carcione ◽  
Stefano Picotti ◽  
Jian Zhou ◽  
...  

Rock-physics templates establish a link between seismic properties (e.g., velocity, density, impedance, and attenuation) and reservoir properties such as porosity, fluid saturation, permeability, and clay content. We focus on templates based on attenuation (seismic [Formula: see text] or quality factor), which are highly affected by those properties, and we consider carbonate reservoirs that constitute 60% of the world oil reserves and a potential for additional gas reserves. The seismic properties are described with mesoscopic-loss models, such as the White model of patchy saturation and the double double-porosity model, which include frame and fluid heterogeneities. We have performed ultrasonic experiments, and we estimate the attenuation of the samples and the reservoir by using the spectral ratio method and the improved frequency-shift method. Then, multiscale calibrations of the templates are performed by using laboratory, well log, and seismic data. On this basis, reservoir porosity and fluid saturation are quantitatively evaluated. We first apply the templates to ultrasonic data of limestone using the White model. Then, we consider seismic data of a carbonate gas reservoir of MX work area in the Sichuan Basin, southwest China. A survey line in the area is selected to detect the reservoir by using the templates. The results indicate that the estimated porosity and saturation are consistent with well-log data and actual gas production results. The methodology indicates that the microstructural characteristics of a high-quality reservoir can effectively be predicted using seismic [Formula: see text].


Sign in / Sign up

Export Citation Format

Share Document