scholarly journals An Experimental Study of Velocity-Saturation Relationships in Volcanic Rocks

2015 ◽  
Vol 8 (1) ◽  
pp. 142-152 ◽  
Author(s):  
Zhidi Liu ◽  
Jingzhou Zhao

In this paper, experiments are carried out under different pressures and water saturations using core samples of volcanic rocks from the Junggar Basin in China to understand how water saturation affects P- and S-wave velocities. The results show that water saturated rocks exhibit significantly higher P- and S-wave velocities than gas saturated rocks. In addition, the P- and S-wave velocity ratio declines with increasing water saturation. Furthermore, a P- and S-wave velocity ratio vs. resistivity cross plot is created to identify gas reservoirs in the volcanic rocks in the Junggar Basin.

Geophysics ◽  
2019 ◽  
Vol 84 (6) ◽  
pp. MR205-MR222 ◽  
Author(s):  
Sheyore John Omovie ◽  
John P. Castagna

In situ P- and S-wave velocity measurements in a variety of organic-rich shales exhibit P-to-S-wave velocity ratios that are significantly lower than lithologically similar fully brine-saturated shales having low organic content. It has been hypothesized that this drop could be explained by the direct influence of kerogen on the rock frame and/or by the presence of free hydrocarbons in the pore space. The correlation of hydrocarbon saturation with total organic content in situ makes it difficult to separate these possible mechanisms using log data alone. Theoretical bounding equations, using pure kerogen as an end-member component without associated gas, indicate that kerogen reduces the P- and S-wave velocities but does not in general reduce their ratio enough to explain the observed low velocity ratio. The theoretical modeling is consistent with ultrasonic measurements on organic shale core samples that indicate no dependence of velocity ratios on the kerogen volume alone. Sonic log measurements of P- and S-wave velocities in seven organic-rich shale formations deviate significantly (typically more than 5%) from the Greenberg-Castagna empirical brine-saturated shale trend toward lower velocity ratios. In these formations, and on core measurements, Gassmann fluid substitution to 100% brine saturation yields velocity ratios consistent with the Greenberg-Castagna velocity trend for fully brine-saturated shales, despite the high organic content. These sonic and ultrasonic measurements, as well as theoretical modeling, suggest that the velocity ratio reduction in organic shales is best explained by the presence of free hydrocarbons.


Author(s):  
Lei Wang ◽  
Erik Rybacki ◽  
Audrey Bonnelye ◽  
Marco Bohnhoff ◽  
Georg Dresen

Abstract Knowledge of pressure-dependent static and dynamic moduli of porous reservoir rocks is of key importance for evaluating geological setting of a reservoir in geo-energy applications. We examined experimentally the evolution of static and dynamic bulk moduli for porous Bentheim sandstone with increasing confining pressure up to about 190 MPa under dry and water-saturated conditions. The static bulk moduli (Ks) were estimated from stress–volumetric strain curves while dynamic bulk moduli (Kd) were derived from the changes in ultrasonic P- and S- wave velocities (~ 1 MHz) along different traces, which were monitored simultaneously during the entire deformation. In conjunction with published data of other porous sandstones (Berea, Navajo and Weber sandstones), our results reveal that the ratio between dynamic and static bulk moduli (Kd/Ks) reduces rapidly from about 1.5 − 2.0 at ambient pressure to about 1.1 at high pressure under dry conditions and from about 2.0 − 4.0 to about 1.5 under water-saturated conditions, respectively. We interpret such a pressure-dependent reduction by closure of narrow (compliant) cracks, highlighting that Kd/Ks is positively correlated with the amount of narrow cracks. Above the crack closure pressure, where equant (stiff) pores dominate the void space, Kd/Ks is almost constant. The enhanced difference between dynamic and static bulk moduli under water saturation compared to dry conditions is possibly caused by high pore pressure that is locally maintained if measured using high-frequency ultrasonic wave velocities. In our experiments, the pressure dependence of dynamic bulk modulus of water-saturated Bentheim sandstone at effective pressures above 5 MPa can be roughly predicted by both the effective medium theory (Mori–Tanaka scheme) and the squirt-flow model. Static bulk moduli are found to be more sensitive to narrow cracks than dynamic bulk moduli for porous sandstones under dry and water-saturated conditions.


2021 ◽  
Vol 9 ◽  
Author(s):  
Rupeng Ma ◽  
Jing Ba ◽  
José Carcione ◽  
Maxim Lebedev ◽  
Changsheng Wang

The petrophysical properties can be proper indicators to identify oil and gas reservoirs, since the pore fluids have significant effects on the wave response. We have performed ultrasonic measurements on two sets of tight siltstones and dolomites at partial saturation. P- and S-wave velocities are obtained by the pulse transmission technique, while attenuation is calculated using the centroid-frequency shift and spectral-ratio methods. The fluid sensitivities of different properties (i.e., P- and S-wave velocities, impedances and attenuation, Poisson's ratio, density, and their combinations) are quantitatively analyzed by considering the data distribution, based on the crossplot technique. The result shows that the properties (P- to S-wave velocity and attenuation ratios, Poisson's ratio, and first to second Lamé constant ratio) with high fluid-sensitivity indicators successfully distinguish gas from oil and water, unlike oil from water. Moreover, siltstones and dolomites can be identified on the basis of data distribution areas. Ultrasonic rock-physics templates of the P- to S-wave velocity ratio vs. the product of first Lamé constant with density obtained with a poroelastic model, considering the structural heterogeneity and patchy saturation, are used to predict the saturation and porosity, which are in good agreement with the experimental data at different porosity ranges.


1968 ◽  
Vol 8 (04) ◽  
pp. 389-404 ◽  
Author(s):  
A.L. Podio ◽  
A.R. Gregory ◽  
K.E. Gray

Abstract Dynamic elastic properties of dry and water-saturated Green River shale samples were computed from compressional and shear-wave velocity measurements. P- and S-wave velocity measurements were made in three mutually perpendicular directions with respect to the bedding planes. Measurements were also made in several different directions by varying the angle between the bedding planes and the direction of propagation of the wave for angles of 0, 30, 45, 60 and 90 degrees. The oriented samples were subjected to both confining pressure and axial loads, in excess of the confining stress, in the direction of Propagation. In general, P- and S-wave velocities increased with increasing stress levels, with a corresponding increase in Young's modulus. Water saturation caused the P-wave velocity to increase and the S-wave velocity to decrease. Elastic moduli decreased upon saturation, except for Poisson's ratio, which increased, indicating some degree of weakening of the material. The samples showed a moderate degree of anisotropy; this was to be expected from the laminated nature and shallow occurence of Green River shale. Introduction This paper presents some results of an experimental determination of the elastic coefficients of anisotropic materials (in particular, finely layered rocks and minerals such as Green River shale) from measurements of dilatational and shear-ultrasonic-wave velocities. Ultrasonic techniques have been used extensively in nondestructive testing. Several methods have been proposed by McSkimmin, and some of these have been used successfully to measure ultrasonic velocities in rocks. Hughes and Cross, Wyllie et al., and Birch, developed pulse first-arrival techniques for the measurement of dilatational and shear velocities. Williams and Lamb used the method of cancellation of a traveling wave, which was later modified by Myers et al. and perfected by McSkimmin. Although this method is highly accurate, it has not been used as widely as the pulse-transmission methods recently reported by Jamieson and Hoskins, King, and Mattaboni and Schreiber. It has been common practice to use some form of crystal transducers, either quartz or ceramic, that has been cut or polarized in different directions in order to generate either compressional or shear waves. However, accurate determination of shear wave velocities has been difficult due to problems that arise in obtaining a pure shear wave from cross-polarized crystals, which usually also generate a small amount of compressional energy. As reported by Gregory, this energy can be seen as a long precursor preceding the sharp break of the shear first arrival. The need for generating pure shear waves led to interest in mode-conversion techniques, which are based upon conversion of the mode of vibration through wave reflection or refraction at a discontinuity. Arenberg showed that for certain materials and for certain angles of incidence it is possible to generate pure shear modes by reflection at a boundary. Jamieson and Hoskins used a pyrex glass-air interface for generating pure shear waves, and King used this method successfully for measuring shear-wave velocities in rocks. Gregory arrived at a similar result by refraction of a wave at an aluminum-oil interface. A plane compressional wave, traveling in the oil phase, is incident on the aluminum at an angle larger than the critical angle for compressional waves, and thereby generates a purely transverse, plane-polarized wave in the aluminum. During the last few years methods have been developed that allow the simultaneous determination of shear and compressional velocities in solids. SPEJ P. 389ˆ


Geophysics ◽  
1985 ◽  
Vol 50 (4) ◽  
pp. 571-581 ◽  
Author(s):  
J. P. Castagna ◽  
M. L. Batzle ◽  
R. L. Eastwood

New velocity data in addition to literature data derived from sonic log, seismic, and laboratory measurements are analyzed for clastic silicate rocks. These data demonstrate simple systematic relationships between compressional and shear wave velocities. For water‐saturated clastic silicate rocks, shear wave velocity is approximately linearly related to compressional wave velocity and the compressional‐to‐shear velocity ratio decreases with increasing compressional velocity. Laboratory data for dry sandstones indicate a nearly constant compressional‐to‐shear velocity ratio with rigidity approximately equal to bulk modulus. Ideal models for regular packings of spheres and cracked solids exhibit behavior similar to the observed water‐saturated and dry trends. For dry rigidity equal to dry bulk modulus, Gassmann’s equations predict velocities in close agreement with data from the water‐saturated rock.


Geophysics ◽  
1994 ◽  
Vol 59 (1) ◽  
pp. 87-92 ◽  
Author(s):  
Gary Mavko ◽  
Diane Jizba

Seismic velocity dispersionin fluid-saturated rocks appears to be dominated by tow mecahnisms: the large scale mechanism modeled by Biot, and the local flow or squirt mecahnism. The tow mechanisms can be distuinguished by the ratio of P-to S-wave dispersions, or more conbeniently, by the ratio of dynamic bulk to shear compliance dispersions derived from the wave velocities. Our formulation suggests that when local flow denominates, the dispersion of the shear compliance will be approximately 4/15 the dispersion of the compressibility. When the Biot mechanism dominates, the constant of proportionality is much smaller. Our examination of ultrasonic velocities from 40 sandstones and granites shows that most, but not all, of the samples were dominated by local flow dispersion, particularly at effective pressures below 40 MPa.


2021 ◽  
Author(s):  
Yair Gordin ◽  
Thomas Bradley ◽  
Yoav O. Rosenberg ◽  
Anat Canning ◽  
Yossef H. Hatzor ◽  
...  

Abstract The mechanical and petrophysical behavior of organic-rich carbonates (ORC) is affected significantly by burial diagenesis and the thermal maturation of their organic matter. Therefore, establishing Rock Physics (RP) relations and appropriate models can be valuable in delineating the spatial distribution of key rock properties such as the total organic carbon (TOC), porosity, water saturation, and thermal maturity in the petroleum system. These key rock properties are of most importance to evaluate during hydrocarbon exploration and production operations when establishing a detailed subsurface model is critical. High-resolution reservoir models are typically based on the inversion of seismic data to calculate the seismic layer properties such as P- and S-wave impedances (or velocities), density, Poisson's ratio, Vp/Vs ratio, etc. If velocity anisotropy data are also available, then another layer of data can be used as input for the subsurface model leading to a better understanding of the geological section. The challenge is to establish reliable geostatistical relations between these seismic layer measurements and petrophysical/geomechanical properties using well logs and laboratory measurements. In this study, we developed RP models to predict the organic richness (TOC of 1-15 wt%), porosity (7-35 %), water saturation, and thermal maturity (Tmax of 420-435⁰C) of the organic-rich carbonate sections using well logs and laboratory core measurements derived from the Ness 5 well drilled in the Golan Basin (950-1350 m). The RP models are based primarily on the modified lower Hashin-Shtrikman bounds (MLHS) and Gassmann's fluid substitution equations. These organic-rich carbonate sections are unique in their relatively low burial diagenetic stage characterized by a wide range of porosity which decreases with depth, and thermal maturation which increases with depth (from immature up to the oil window). As confirmation of the method, the levels of organic content and maturity were confirmed using Rock-Eval pyrolysis data. Following the RP analysis, horizontal (HTI) and vertical (VTI) S-wave velocity anisotropy were analyzed using cross-dipole shear well logs (based on Stoneley waves response). It was found that anisotropy, in addition to the RP analysis, can assist in delineating the organic-rich sections, microfractures, and changes in gas saturation due to thermal maturation. Specifically, increasing thermal maturation enhances VTI and azimuthal HTI S-wave velocity anisotropies, in the ductile and brittle sections, respectively. The observed relationships are quite robust based on the high-quality laboratory and log data. However, our conclusions may be limited to the early stages of maturation and burial diagenesis, as at higher maturation and diagenesis the changes in physical properties can vary significantly.


Geosciences ◽  
2020 ◽  
Vol 10 (9) ◽  
pp. 337
Author(s):  
Hanane Sghiouri El Idrissi ◽  
Abderrahim Samaouali ◽  
Younes El Rhaffari ◽  
Salah El Alami ◽  
Yves Geraud

In this work, we study the variability of the lithological composition and organic matter content of samples were taken from the different layers M, X and Y of the Timahdit oil shale in Morocco, in order to experimentally analyze the impact of this variability on petrophysical measurements. The objective of this study is to predict the properties of the layers, including their thermal conductivity, thermal diffusivity, porosity and P and S wave velocities. The results of the study of the impact of the organic matter content of the samples on the petrophysical measurements show that, regardless of the organic matter content, thermal conductivity and diffusivity remain insensitive, while P and S wave velocities decrease linearly and porosity increases with increasing organic matter content. On the other hand, the study of the organic matter variability content is consistent with the velocity ratio, so can be used as an organic matter indicator of the layers. Conductivity and thermal diffusivity are almost invariant to the variability of the organic matter.


Geophysics ◽  
2018 ◽  
Vol 83 (3) ◽  
pp. KS11-KS22 ◽  
Author(s):  
Nicola Piana Agostinetti ◽  
Alberto Malinverno

We use teleseismic P-to-S converted waves from a permanent station to estimate the uncertainties in a 1D elastic model of the shallow crust (0–7 km depth) obtained from the inversion of receiver function (RF) data. Our earth model consists of layers with a constant S-wave velocity [Formula: see text] and P- to S-wave velocity ratio ([Formula: see text]). We apply a Bayesian formulation and transdimensional Monte Carlo sampling to compute the posterior uncertainties of the earth model. The model uncertainties rely on a realistic representation of the data uncertainties, and we estimate directly from the stacking of the teleseismic data, a full-error covariance matrix. To explore the effect of the number of teleseismic events and the RF frequency content, we compare the results of inverting a single RF computed for a cut-off filter frequency of 4 Hz with the joint inversion of four RFs computed from independent ensembles in a larger pool of events for cut-off frequencies of 0.5, 1, 2, and 4 Hz. The inversion results are compared with the lithostratigraphy and sonic-log measurements from a 7 km deep borehole drilled near the seismic station. The inversion of a single RF results in larger uncertainties in the recovered [Formula: see text] profile and in the depth to seismic discontinuities compared with the multifrequency inversion. Moreover, the multifrequency inversion predicts more accurately the depth to a velocity inversion at approximately 6 km below the surface and matches more closely the borehole sonic-log data. Our results indicate that RF data can be used to map shallow (3–5 km depth) crustal interfaces with uncertainties in the order of 300–500 m, whereas uncertainties are consistently smaller (<300 m) for interfaces in the top kilometer.


Sign in / Sign up

Export Citation Format

Share Document