P- and S-wave attenuation logs from monopole sonic data

Geophysics ◽  
2000 ◽  
Vol 65 (3) ◽  
pp. 755-765 ◽  
Author(s):  
Xinhua Sun ◽  
Xiaoming Tang ◽  
C. H. (Arthur) Cheng ◽  
L. Neil Frazer
Keyword(s):  
Wave Attenuation ◽  
Fluid Velocity ◽  
P Wave ◽  
Reservoir Rock ◽  
Rock Properties ◽  
S Wave ◽  
Intrinsic Attenuation ◽  
Modified Method ◽  
Receiver Array ◽  
Attenuation Profiles

In this paper, a modification of an existing method for estimating relative P-wave attenuation is proposed. By generating synthetic waveforms without attenuation, the variation of geometrical spreading related to changes in formation properties with depth can be accounted for. With the modified method, reliable P- and S-wave attenuation logs can be extracted from monopole array acoustic waveform log data. Synthetic tests show that the P- and S-wave attenuation values estimated from synthetic waveforms agree well with their respective model values. In‐situ P- and S-wave attenuation profiles provide valuable information about reservoir rock properties. Field data processing results show that this method gives robust estimates of intrinsic attenuation. The attenuation profiles calculated independently from each waveform of an eight‐receiver array are consistent with one another. In fast formations where S-wave velocity exceeds the borehole fluid velocity, both P-wave attenuation ([Formula: see text]) and S-wave attenuation ([Formula: see text]) profiles can be obtained. P- and S-wave attenuation profiles and their comparisons are presented for three reservoirs. Their correlations with formation lithology, permeability, and fractures are also presented.

scholarly journals Variational Bayesian inversion of seismic attributes jointly for geologic facies and petrophysical rock properties

Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. MR213-MR233 ◽  
Author(s):  
Muhammad Atif Nawaz ◽  
Andrew Curtis ◽  
Mohammad Sadegh Shahraeeni ◽  
Constantin Gerea
Keyword(s):  
Gaussian Mixture ◽  
Seismic Attributes ◽  
Joint Estimation ◽  
Training Data ◽  
P Wave ◽  
Rock Properties ◽  
Inversion Method ◽  
Computationally Efficient ◽  
S Wave ◽  
Well Data

Seismic attributes (derived quantities) such as P-wave and S-wave impedances and P-wave to S-wave velocity ratios may be used to classify subsurface volume of rock into geologic facies (distinct lithology-fluid classes) using pattern recognition methods. Seismic attributes may also be used to estimate subsurface petrophysical rock properties such as porosity, mineral composition, and pore-fluid saturations. Both of these estimation processes are conventionally carried out independent of each other and involve considerable uncertainties, which may be reduced significantly by a joint estimation process. We have developed an efficient probabilistic inversion method for joint estimation of geologic facies and petrophysical rock properties. Seismic attributes and petrophysical properties are jointly modeled using a Gaussian mixture distribution whose parameters are initialized by unsupervised learning using well-log data. Rock-physics models may be used in our method to augment the training data if the existing well data are limited; however, this is not required if sufficient well data are available. The inverse problem is solved using the Bayesian paradigm that models uncertainties in the form of probability distributions. Probabilistic inference is performed using variational optimization, which is a computationally efficient deterministic alternative to the commonly used sampling-based stochastic inference methods. With the help of a real data application from the North Sea, we find that our method is computationally efficient, honors expected spatial correlations of geologic facies, allows reliable detection of convergence, and provides full probabilistic results without stochastic sampling of the posterior distribution.

scholarly journals A method of geo-acoustic parameter inversion in shallow sea by the Bayesian theory and the acoustic pressure field

2019 ◽  
Vol 283 ◽  
pp. 06003
Author(s):  
Guangxue Zheng ◽  
Hanhao Zhu ◽  
Jun Zhu
Keyword(s):  
Sound Pressure ◽  
Pressure Field ◽  
Wave Attenuation ◽  
Acoustic Parameter ◽  
Bayesian Theory ◽  
P Wave ◽  
Shallow Sea ◽  
S Wave ◽  
Wave Velocities ◽  
Parameter Inversion

A method of geo-acoustic parameter inversion based on the Bayesian theory is proposed for the acquisition of acoustic parameters in shallow sea with the elastic seabed. Firstly, the theoretical prediction value of the sound pressure field is calculated by the fast field method (FFM). According to the Bayesian theory, we establish the misfit function between the measured sound pressure field and the theoretical pressure field. It is under the assumption of Gaussian data errors which are in line with the likelihood function. Finally, the posterior probability density (PPD) of parameters is given as the result of inversion. Our research is conducted in the light of Metropolis sample rules. Apart from numerical simulations, a scaled model experiment has been taken in the laboratory tank. The results of numerical simulations and tank experiments show that sound pressure field calculated by the result of inversion is consistent with the measured sound pressure field. Besides, s-wave velocities, p-wave velocities and seafloor density have fewer uncertainties and are more sensitive to complex sound pressure than s-wave attenuation and p-wave attenuation. The received signals calculated by inversion results are keeping with received signals in the experiment which verify the effectiveness of this method.

Rock-physics modeling of ultrasonic P- and S-wave attenuation in partially frozen brine and unconsolidated sand systems and comparison with laboratory measurements

Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. MR153-MR171 ◽  
Author(s):  
Linsen Zhan ◽  
Jun Matsushima
Keyword(s):  
Ultrasonic Wave ◽  
Wave Attenuation ◽  
Freezing Point ◽  
Rock Physics ◽  
Ultrasonic Waves ◽  
P Wave ◽  
Dominant Factor ◽  
S Wave ◽  
High Attenuation ◽  
Unconsolidated Sand

The nonintuitive observation of the simultaneous high velocity and high attenuation of ultrasonic waves near the freezing point of brine was previously measured in partially frozen systems. However, previous studies could not fully elucidate the attenuation variation of ultrasonic wave propagation in a partially frozen system. We have investigated the potential attenuation mechanisms responsible for previously obtained laboratory results by modeling ultrasonic wave transmission in two different partially frozen systems: partially frozen brine (two phases composed of ice and unfrozen brine) and unconsolidated sand (three phases composed of ice, unfrozen brine, and sand). We adopted two different rock-physics models: an effective medium model for partially frozen brine and a three-phase extension of the Biot model for partially frozen unconsolidated sand. For partially frozen brine, our rock-physics study indicated that squirt flow caused by unfrozen brine inclusions in porous ice could be responsible for high P-wave attenuation around the freezing point. Decreasing P-wave attenuation below the freezing point can be explained by the gradual decrease of squirt flow due to the gradual depletion of unfrozen brine. For partially frozen unconsolidated sand, our rock-physics study implied that squirt flow between ice grains is a dominant factor for P-wave attenuation around the freezing point. With decreasing temperature lower than the freezing point, the friction between ice and sand grains becomes more dominant for P-wave attenuation because the decreasing amount of unfrozen brine reduces squirt flow between ice grains, whereas the generation of ice increases the friction. The increasing friction between ice and sand grains caused by ice formation is possibly responsible for increasing the S-wave attenuation at decreasing temperatures. Then, further generation of ice with further cooling reduces the elastic contrast between ice and sand grains, hindering their relative motion; thus, reducing the P- and S-wave attenuation.

Useful approximations for converted‐wave AVO

Geophysics ◽  
2001 ◽  
Vol 66 (6) ◽  
pp. 1721-1734 ◽  
Author(s):  
Antonio C. B. Ramos ◽  
John P. Castagna
Keyword(s):  
Wave Velocity ◽  
Random Noise ◽  
P Wave ◽  
Rock Properties ◽  
Negative Slope ◽  
High Contrast ◽  
Fractional Change ◽  
S Wave ◽  
Converted Wave ◽  
Low Contrast

Converted‐wave amplitude versus offset (AVO) behavior may be fit with a cubic relationship between reflection coefficient and ray parameter. Attributes extracted using this form can be directly related to elastic parameters with low‐contrast or high‐contrast approximations to the Zoeppritz equations. The high‐contrast approximation has the advantage of greater accuracy; the low‐contrast approximation is analytically simpler. The two coefficients of the low‐contrast approximation are a function of the average ratio of compressional‐to‐shear‐wave velocity (α/β) and the fractional changes in S‐wave velocity and density (Δβ/β and Δρ/ρ). Because of its simplicity, the low‐contrast approximation is subject to errors, particularly for large positive contrasts in P‐wave velocity associated with negative contrasts in S‐wave velocity. However, for incidence angles up to 40° and models confined to |Δβ/β| < 0.25, the errors in both coefficients are relatively small. Converted‐wave AVO crossplotting of the coefficients of the low‐contrast approximation is a useful interpretation technique. The background trend in this case has a negative slope and an intercept proportional to the α/β ratio and the fractional change in S‐wave velocity. For constant α/β ratio, an attribute trace formed by the weighted sum of the coefficients of the low‐contrast approximation provides useful estimates of the fractional change in S‐wave velocity and density. Using synthetic examples, we investigate the sensitivity of these parameters to random noise. Integrated P‐wave and converted‐wave analysis may improve estimation of rock properties by combining extracted attributes to yield fractional contrasts in P‐wave and S‐wave velocities and density. Together, these parameters may provide improved direct hydrocarbon indication and can potentially be used to identify anomalies caused by low gas saturations.

Seismic wave attenuation and modulus dispersion in sandstones

Geophysics ◽  
2016 ◽  
Vol 81 (3) ◽  
pp. D211-D231 ◽  
Author(s):  
James W. Spencer ◽  
Jacob Shine
Keyword(s):  
Wave Attenuation ◽  
High Viscosity ◽  
P Wave ◽  
S Wave ◽  
Local Flow ◽  
Shear Moduli ◽  
Low Viscosity ◽  
Confining Pressures ◽  
Wave Induced ◽  
Dispersion And Attenuation

We have conducted laboratory experiments over the 1–200 Hz band to examine the effects of viscosity and permeability on modulus dispersion and attenuation in sandstones and also to examine the effects of partial gas or oil saturation on velocities and attenuations. Our results have indicated that bulk modulus values with low-viscosity fluids are close to the values predicted using Gassmann’s first equation, but, with increasing frequency and viscosity, the bulk and shear moduli progressively deviate from the values predicted by Gassmann’s equations. The shear moduli increase up to 1 GPa (or approximately 10%) with high-viscosity fluids. The P- and S-wave attenuations ([Formula: see text] and [Formula: see text]) and modulus dispersion with different fluids are indicative of stress relaxations that to the first order are scaling with frequency times viscosity. By fitting Cole-Cole distributions to the scaled modulus and attenuation data, we have found that there are similar P-wave, shear and bulk relaxations, and attenuation peaks in each of the five sandstones studied. The modulus defects range from 11% to 15% in Berea sandstone to 16% to 26% in the other sandstones, but these would be reduced at higher confining pressures. The relaxations shift to lower frequencies as the viscosity increased, but they do not show the dependence on permeability predicted by mesoscopic wave-induced fluid flow (WIFF) theories. Results from other experiments having patchy saturation with liquid [Formula: see text] and high-modulus fluids are consistent with mesoscopic WIFF theories. We have concluded that the modulus dispersion and attenuations ([Formula: see text] and [Formula: see text]) in saturated sandstones are caused by a pore-scale, local-flow mechanism operating near grain contacts.

The effects of pore‐fluid salinity on ultrasonic wave propagation in sandstones

Geophysics ◽  
1998 ◽  
Vol 63 (3) ◽  
pp. 928-934 ◽  
Author(s):  
Simon M. Jones ◽  
Clive McCann ◽  
Timothy R. Astin ◽  
Jeremy Sothcott
Keyword(s):  
Wave Propagation ◽  
Ultrasonic Wave ◽  
Wave Attenuation ◽  
Clay Content ◽  
Seismic Wave Propagation ◽  
Ultrasonic Pulse ◽  
Pore Fluid ◽  
P Wave ◽  
S Wave ◽  
Shear Moduli

Petrophysical interpretation of increasingly refined seismic data from subsurface formations requires a more fundamental understanding of seismic wave propagation in sedimentary rocks. We consider the variation of ultrasonic wave velocity and attenuation in sandstones with pore‐fluid salinity and show that wave propagation is modified in proportion to the clay content of the rock and the salinity of the pore fluid. Using an ultrasonic pulse reflection technique (590–890 kHz), we have measured the P-wave and S-wave velocities and attenuations of 15 saturated sandstones with variable effective pressure (5–60 MPa) and pore‐fluid salinity (0.0–3.4 M). In clean sandstones, there was close agreement between experimental and Biot model values of [Formula: see text], but they diverged progressively in rocks containing more than 5% clay. However, this effect is small: [Formula: see text] changed by only 0.6% per molar change in salinity for a rock with a clay content of 29%. The variation of [Formula: see text] with brine molarity exhibited Biot behavior in some samples but not in others; there was no obvious relationship with clay content. P-wave attenuation was independent of pore‐fluid salinity, while S-wave attenuation was weakly dependent. The velocity data suggest the frame bulk and shear moduli of sandstones are altered by changes in the pore‐fluid salinity. One possible mechanism is the formation damage caused by clay swelling and migration of fines in low‐molarity electrolytes. The absence of variation between the attenuation in water‐saturated and brine‐saturated samples indicates the attenuation mechanism is relatively unaffected by changes in the frame moduli.

Converted-wave model building and imaging based on common-focus-point methodology

Geophysics ◽  
2020 ◽  
Vol 85 (6) ◽  
pp. U139-U149
Author(s):  
Hongwei Liu ◽  
Mustafa Naser Al-Ali ◽  
Yi Luo
Keyword(s):  
Elastic Wave ◽  
Wave Velocity ◽  
Seismic Data ◽  
Model Building ◽  
P Wave ◽  
Rock Properties ◽  
S Wave ◽  
Converted Wave ◽  
Seismic Images ◽  
Focus Point

Seismic images can be viewed as photographs for underground rocks. These images can be generated from different reflections of elastic waves with different rock properties. Although the dominant seismic data processing is still based on the acoustic wave assumption, elastic wave processing and imaging have become increasingly popular in recent years. A major challenge in elastic wave processing is shear-wave (S-wave) velocity model building. For this reason, we have developed a sequence of procedures for estimating seismic S-wave velocities and the subsequent generation of seismic images using converted waves. We have two main essential new supporting techniques. The first technique is the decoupling of the S-wave information by generating common-focus-point gathers via application of the compressional-wave (P-wave) velocity on the converted seismic data. The second technique is to assume one common VP/ VS ratio to approximate two types of ratios, namely, the ratio of the average earth layer velocity and the ratio of the stacking velocity. The benefit is that we reduce two unknown ratios into one, so it can be easily scanned and picked in practice. The PS-wave images produced by this technology could be aligned with the PP-wave images such that both can be produced in the same coordinate system. The registration between the PP and PS images provides cross-validation of the migrated structures and a better estimation of underground rock and fluid properties. The S-wave velocity, computed from the picked optimal ratio, can be used not only for generating the PS-wave images, but also to ensure well registration between the converted-wave and P-wave images.

Attenuation anisotropy in the linear-slip model: Interpretation of physical modeling data

Geophysics ◽  
2009 ◽  
Vol 74 (5) ◽  
pp. WB165-WB176 ◽  
Author(s):  
Tatiana Chichinina ◽  
Irina Obolentseva ◽  
Leonid Gik ◽  
Boris Bobrov ◽  
Gerardo Ronquillo-Jarillo
Keyword(s):  
Symmetry Axis ◽  
Wave Attenuation ◽  
Joint Inversion ◽  
P Wave ◽  
Slip Model ◽  
S Wave ◽  
The Real ◽  
Complex Valued ◽  
Stack Model ◽  
Linear Slip Model

This study attempts to validate a mathematical formalism of introducing attenuation into Schoenberg’s linear slip model. This formalism is based on replacing the real-valued weaknesses by complex-valued ones. During an ultrasonic experiment, performed at a central frequency of [Formula: see text] on a plate-stack model with [Formula: see text]-thick Plexiglas™ plates, the velocity and attenuation (inverse of the quality factor [Formula: see text]) of P-, SH-, and SV-waves are measured in directions from 25° to 90° with the symmetry axis for dry and oil-saturated models and loading uniaxial pressures of 2 and [Formula: see text]. The velocity and attenuation data are fitted by the derived theoretical functions. The values of the real and imaginary parts of the complex-valued weaknesses are estimated. Thereal parts of the weaknesses, which have a clear physical meaning (they affect the weakening of the material), are three times larger for the dry model than for the oil-saturated one. The imaginary parts of the weaknesses are responsible for attenuation; their values are an order of magnitude smaller than the real parts. The derived expressions for angle-dependent velocities and attenuations can be used to distinguish between dry and oil-saturated fractures. In particular, the P-wave attenuation function in the symmetry-axis direction (normal to fracture planes) is different in dry and saturated media. The experiment shows that the plate-stack model is inhomogeneous because of the nonuniform pressure distribution, which degrades the experimental results and creates difficulties in the inversion for the complex-valued weaknesses — particularly in joint inversion of P- and S-wave data.

Sign in / Sign up

Export Citation Format

Share Document