Effect of fluids and frequencies on Poisson’s ratio of sandstone samples

Geophysics ◽  
2016 ◽  
Vol 81 (2) ◽  
pp. D183-D195 ◽  
Author(s):  
Lucas Pimienta ◽  
Jérôme Fortin ◽  
Yves Guéguen

Poisson’s ratio [Formula: see text] is an important parameter when interpreting measured geophysical and seismic data. For an isotropic medium, it directly relates to the ratio of P- and S-wave velocities. We have measured [Formula: see text] as a function of pressure and frequency in fluid-saturated sandstones. The method of measuring [Formula: see text] was first tested as a function of pressure and frequency using standard samples. The phase shift [Formula: see text] between radial and axial strains was also measured. For all standard samples, such as the linear viscoelastic Plexiglas, the data indicated that [Formula: see text] correlated with [Formula: see text] and related to a dissipation on [Formula: see text]. Then, [Formula: see text] and [Formula: see text] were measured as a function of pressure and frequency for two dry and fluid-saturated Fontainebleau sandstone samples. Under dry conditions, no frequency dependence and very small pressure dependence were observed. Unusual behaviors were observed under fluid-saturated conditions. In particular, [Formula: see text] of one sample indicated a frequency-dependent bell-shaped dispersion under water and glycerin saturation that correlated with peaks in [Formula: see text]. Plotting the measurements as a function of apparent frequency (i.e., normalizing by the fluid viscosity) indicated a good fit between the water- and glycerin-saturated measurements. The bell-shaped dispersion in [Formula: see text] that was observed for one particular sandstone held for all effective pressures. These variations fully correlated with the peaks of [Formula: see text] observed. Our results can be interpreted using fluid flow and effective medium theories in the case of a porous microcracked rock. Drained/undrained and relaxed/unrelaxed transitions have frequency and magnitude of variations that are consistent with the measurements. The rock sample microcrack density strongly affects this frequency dependence. The inferred [Formula: see text] ratio at low effective pressures also indicates a large frequency-dependent bell-shaped dispersion. The parameter [Formula: see text] is a clear indicator of the frequency-dependent dissipation of [Formula: see text] and relates to the attenuation of P- and S-waves.

1982 ◽  
Vol 72 (1) ◽  
pp. 129-149
Author(s):  
S. W. Roecker ◽  
B. Tucker ◽  
J. King ◽  
D. Hatzfeld

abstract Digital recordings of microearthquake codas from shallow and intermediate depth earthquakes in the Hindu Kush region of Afghanistan were used to determine the attenuation factors of the S-wave coda (Qc) and primary S waves (Qβ). An anomalously rapid decay of the coda shortly after the S-wave arrival, observed also in a study of coda in central Asia by Rautian and Khalturin (1978), seems to be due primarily to depth-dependent variations in Qc. In particular, we deduce the average Qc in the crust and uppermost mantle (<100-km depth) is approximately four times lower than the deeper mantle (<400-km depth) over a wide frequency range (0.4 to 24 Hz). Further, while Qc generally increases with frequency at any depth, the degree of frequency dependence of Qc depends on depth. Except at the highest frequency studied here (∼48 Hz), the magnitude of Qc at a particular frequency increases with depth while its frequency dependence decreases. For similar depths, determinations of Qβ and Qc agree, suggesting a common wave composition and attenuation mechanism for S waves and codas. Comparison of these determinations of Qc in Afghanistan with those in other parts of the world shows that the degree of frequency dependence of Qc correlates with the expected regional heterogeneity. Such a correlation supports the prejudice that Qc is primarily influenced by scattering and suggests that tectonic processes such as folding and faulting are instrumental in creating scattering environments.


Geophysics ◽  
2000 ◽  
Vol 65 (2) ◽  
pp. 559-564 ◽  
Author(s):  
Ran Bachrach ◽  
Jack Dvorkin ◽  
Amos M. Nur

We determined P- and S-wave velocity depth profiles in shallow, unconsolidated beach sand by analyzing three‐component surface seismic data. P- and S-wave velocity profiles were calculated from traveltime measurements of vertical and tangential component seismograms, respectively. The results reveal two discrepancies between theory and data. Whereas both velocities were found to be proportional to the pressure raised to the power of 1/6, as predicted by the Hertz‐Mindlin contact theory, the actual values of the velocities are less than half of those calculated from this theory. We attribute this discrepancy to the angularity of the sand grains. Assuming that the average radii of curvature at the grain contacts are smaller than the average radii of the grains, we modify the Hertz‐Mindlin theory accordingly. We found that the ratio of the contact radius to the grain radius is about 0.086. The second disparity is between the observed Poisson’s ratio of 0.15 and the theoretical value (0.008 for random pack of quartz spheres). This discrepancy can be reconciled by assuming slip at the grain contacts. Because slip decreases the shearing between grains, Poisson’s ratio increases.


Geophysics ◽  
2018 ◽  
Vol 83 (4) ◽  
pp. T209-T234 ◽  
Author(s):  
Jing-Bo Chen ◽  
Jian Cao

Because of its high computational cost, we needed to develop an efficient numerical scheme for the frequency-domain 3D elastic wave equation. In addition, the numerical scheme should be applicable to media with a liquid-solid interface. To address these two issues, we have developed a new average-derivative optimal 27-point scheme with arbitrary directional grid intervals and a corresponding numerical dispersion analysis for the frequency-domain 3D elastic wave equation. The novelty of this scheme is that its optimal coefficients depend on the ratio of the directional grid intervals and Poisson’s ratio. In this way, this scheme can be applied to media with a liquid-solid interface and a computational grid with arbitrary directional grid intervals. For media with a variable Poisson’s ratio, we have developed an effective and stable interpolation method for optimization coefficients. Compared with the classic 19-point scheme, this new scheme reduces the required number of grid points per wavelength for equal and unequal directional grid intervals. The reduction of the number of grid points increases as the Poisson’s ratio becomes larger. In particular, the numerical S-wave phase velocity of this new scheme becomes zero, whereas the classic 19-point scheme produces a spurious numerical S-wave phase velocity, as Poisson’s ratio reaches 0.5. We have performed numerical examples to develop the theoretical analysis.


Geophysics ◽  
1982 ◽  
Vol 47 (5) ◽  
pp. 819-824 ◽  
Author(s):  
Harsh K. Gupta ◽  
Ronald W. Ward ◽  
Tzeu‐Lie Lin

Analysis of P‐ and S‐waves from shallow microearthquakes in the vicinity of The Geysers geothermal area, California, recorded by a dense, telemetered seismic array operated by the U.S. Geological Survey (USGS) shows that these phases are easily recognized and traced on record sections to distances of 80 km. Regional average velocities for the upper crust are estimated to be [Formula: see text] and [Formula: see text] for P‐ and S‐waves, respectively. Poisson’s ratio is estimated at 23 locations using Wadati diagrams and is found to vary from 0.13 to 0.32. In general, the Poisson’s ratio is found to be lower at the locations close to the steam production zones at The Geysers and Clear Lake volcanic field to the northeast. The low Poisson ratio corresponds to a decrease in P‐wave velocity in areas of high heat flow. The decrease may be caused by fracturing of the rock and saturation with gas or steam.


1969 ◽  
Vol 59 (3) ◽  
pp. 1189-1200
Author(s):  
John R. McGinley ◽  
Don L. Anderson

abstract The unified magnitude, the ratio of the amiplitudes of S to P waves, and travel-time residuals were compiled from published data for the five Seismological Observatories, TFO, UBO, BMO, WMO and CBO. Using one of the stations as a reference, a relative measure of the above quantities was calculated for each of the other stations for each of a number of earthquakes. The stations in the Basin and Range Province are consistent with a markedly higher attentuation of P waves and a high attenuation of S relative to P when compared to the other stations. This latter observation indicates a high Poisson's ratio in the mantle under the Basin and Range. The delay times to these stations are also consistent with the high Poisson's ratio and with a low-velocity upper mantle. The ratio of the amplitudes of long-period S waves to short-period P waves varies by a factor of 4 among these stations. BMO, in eastern Oregon, has a high S/P amplitude ratio compared to other stations and a travel-time residual that is comparable to the observatories in the mid-continent. This may be another example of a seismic “window” into the upper mantle that is generated by underthrusting of the oceanic lithosphere.


2014 ◽  
Vol 6 (2) ◽  
pp. 2567-2613 ◽  
Author(s):  
P. Tong ◽  
D. Zhao ◽  
D. Yang ◽  
X. Yang ◽  
J. Chen ◽  
...  

Abstract. High-resolution 3-D P and S wave crustal velocity and Poisson's ratio models of the 1992 Landers earthquake (Mw 7.3) area are determined iteratively by a wave-equation based traveltime seismic tomography (WETST) technique as developed in the first paper. The details of data selection, synthetic arrival-time determination, and trade-off analysis of damping and smoothing parameters are presented to show the performance of this new tomographic inversion method. A total of 78 523 P wave and 46 999 S wave high-quality arrival-time data from 2041 local earthquakes recorded by 275 stations during the period of 1992–2013 is used to obtain the final tomographic models which costs around 10 000 CPU h. Checkerboard resolution tests are conducted to verify the reliability of inversion results for the chosen seismic data and the wave-equation based traveltime seismic tomography method. Significant structural heterogeneities are revealed in the crust of the 1992 Lander earthquake area which may be closely related to the local seismic activities. Strong variations of velocity and Poisson's ratio exist in the source regions of the Landers and three other strong earthquakes in this area. Most seismicity occurs in areas with high-velocity and low Poisson's ratio, which may be associated with the seismogenic layer. Pronounced low-velocity anomalies revealed in the lower crust along the Elsinore, the San Jacinto and the San Andreas faults may reflect the existence of fluids in the lower crust. The recovery of these strong heterogeneous structures are facilitated by the use of full wave equation solvers and WETST and verifies their ability in generating high-resolution tomographic models.


Energies ◽  
2021 ◽  
Vol 14 (16) ◽  
pp. 5007
Author(s):  
Stian Rørheim ◽  
Mohammad Hossain Bhuiyan ◽  
Andreas Bauer ◽  
Pierre Rolf Cerasi

Carbon capture and storage (CCS) by geological sequestration comprises a permeable formation (reservoir) for CO2 storage topped by an impermeable formation (caprock). Time-lapse (4D) seismic is used to map CO2 movement in the subsurface: CO2 migration into the caprock might change its properties and thus impact its integrity. Simultaneous forced-oscillation and pulse-transmission measurements are combined to quantify Young’s modulus and Poisson’s ratio as well as P- and S-wave velocity changes in the absence and in the presence of CO2 at constant seismic and ultrasonic frequencies. This combination is the laboratory proxy to 4D seismic because rock properties are monitored over time. It also improves the understanding of frequency-dependent (dispersive) properties needed for comparing in-situ and laboratory measurements. To verify our method, Draupne Shale is monitored during three consecutive fluid exposure phases. This shale appears to be resilient to CO2 exposure as its integrity is neither compromised by notable Young’s modulus and Poisson’s ratio nor P- and S-wave velocity changes. No significant changes in Young’s modulus and Poisson’s ratio seismic dispersion are observed. This absence of notable changes in rock properties is attributed to Draupne being a calcite-poor shale resilient to acidic CO2-bearing brine that may be a suitable candidate for CCS.


2022 ◽  
Author(s):  
◽  
Syuhada, Syuhada

<p>We investigate the mechanisms of seismic anisotropy and attenuation (1/Q) beneath the southernmost part of the Taupo Volcanic Zone (TVZ) by computing variations in S-wave attenuation factors with the direction of wave polarization. We rotate pairs of horizontal components in steps of 22.5◦ from 0◦ to 67.5◦ and into the radial and transverse directions to search for the optimal separation of the attenuation curves and thereby determine an anisotropy symmetry system. The frequency dependence of Q for the rotated S-waves is estimated by means of the non-parametric generalized inversion technique (GIT) of Castro et al. (1990) using shallow earthquakes (< 40 km depth) recorded by GeoNet within 100 km of Mt. Ruapehu. To analyze the effects on computed attenuation properties of source locations, we divide our dataset into two groups: a “TVZ” group containing earthquakes within the TVZ in a distance range of 5–55 km and a “non-TVZ” group containing earthquakes outside the TVZ in a distance range of 5–50 km. To measure Q, we compute the spectral amplitude decay with distance in terms of empirical functions at 20 separate frequencies in the frequency bands 2–10 Hz and 2– 12 Hz for the TVZ and non-TVZ datasets respectively. We construct homogeneous and two-layer Q models for the TVZ dataset based on characteristic features of the attenuation function, while for outside TVZ we only analyse a homogeneous Q model. The homogeneous Q models obtained for the two datasets indicate that S-waves are more attenuated within the TVZ than outside. The homogeneous Q model for the TVZ dataset reveals that the S-wave is anisotropic at high frequencies ( f > 6 Hz) along N–S/E– W directions with the relation QSE ( f ) = (6.15±1.22) f (1.73±0.12) and QSN ( f ) = (4.14± 1.26) f (2.06±0.14), while the non-TVZ dataset shows a weak frequency dependence of attenuation anisotropy at low frequencies in NE–SW/SE–NW directions giving the power law function QSNE ( f ) = (50.93±1.18) f (0.20±0.10) and QSSE ( f ) = (22.60±1.10) f (0.53±0.06). Here, the uncertainty estimates are 95% confidence intervals. To investigate the variation of attenuation anisotropy with depth within the TVZ, we first calculate Q along propagation paths (< 25 km, which corresponds to a maximum turning point depth of 9 km ) and then using paths of 25–55 km length. Small attenuation anisotropy with low attenuation in the N–S direction for the upper crust of TVZ may be related to heterogenous structure as reported by previous studies. Attenuation anisotropy in the northwest direction yielding lower attenuation inferred for the deeper crust suggests the presence of connected melt aligned with the extension direction of TVZ .</p>


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.


Sign in / Sign up

Export Citation Format

Share Document