Mode Conversion Technique Employed in Shear Wave Velocity Studies of Rock Samples Under Axial and Uniform Compression

1967 ◽  
Vol 7 (02) ◽  
pp. 136-148 ◽  
Author(s):  
A.R. Gregory

Abstract A shear wave velocity laboratory apparatus and techniques for testing rock samples under simulated subsurface conditions have been developed. In the apparatus, two electromechanical transducers operating in the frequency range 0.5 to 5.0 megahertz (MHz: megacycles per second) are mounted in contact with each end of the sample. Liquid-solid interfaces of Drakeol-aluminum are used as mode converters. In the generator transducer, there is total mode conversion from P-wave energy to plain S-wave energy, S-wave energy is converted back to P-wave energy in the motor transducer. Similar transducers without mode converters are used to measure P-wave velocities. The apparatus is designed for testing rock samples under axial or uniform loading in the pressure range 0 to 12,000 psi. The transducers have certain advantages over those used by King,1 and the measurement techniques are influenced less by subjective elements than other methods previously reported. An electronic counter-timer having a resolution of 10 nanoseconds measures the transit time of ultrasonic pulses through the sample; elastic wave velocities of most homogeneous materials can be measured with errors of less than 1 percent. S- and P-wave velocity measurements on Bandera sandstone and Solenhofen limestone are reported for the axial pressure range 0 to 6,000 psi and for the uniform pressure range 0 to 10,000 psi. The influence of liquid pore saturants on P- and S-wave velocity is investigated and found to be in broad agreement with Biot's theory. In specific areas, the measurements do not conform to theory. Velocities of samples measured under axial and uniform loading are compared and, in general, velocities measured under uniform stress are higher than those measured under axial stress. Liquid pore fluids cause increases in Poisson's ratio and the bulk modulus but reduce the rigidity modulus, Young's modulus and the bulk compressibility. INTRODUCTION Ultrasonic pulse methods for measuring the shear wave velocity of rock samples in the laboratory have been gradually improved during the last few years. Early experimental pulse techniques reported by Hughes et al.2, and by Gregory3 were beset by uncertainties in determining the first arrival of the shear wave (S-wave) energy. Much of this ambiguity was caused by the multiple modes propagated by piezoelectric crystals and by boundary conversions in the rock specimens. Shear wave velocity data obtained from the critical angle method, described by Schneider and Burton4 and used later by King and Fatt5 and by Gregory,3,6 are of limited accuracy, and interpreting results is too complicated for routine laboratory work. The mode conversion method described by Jamieson and Hoskins7 was recently used by King1 for measuring the S-wave velocities of dry and liquid-saturated rock samples. Glass-air interfaces acted as mode converters in the apparatus, and much of the compressional (P-wave) energy apparently was eliminated from the desired pure shear mode. A more detailed discussion of the current status of laboratory pulse methods applied to geological specimens is given in a review by Simmons.8

2019 ◽  
Author(s):  
Michael Behm ◽  
Feng Cheng ◽  
Anna Patterson ◽  
Gerilyn Soreghan

Abstract. The advent of cable-free nodal arrays for conventional seismic reflection and refraction experiments is changing the acquisition style for active source surveys. Instead of triggering short recording windows for each shot, the nodes are continuously recording over the entire acquisition period from the first to the last shot. The main benefit is a significant increase in geometrical and logistical flexibility. As a by-product, a significant amount of continuous data might also be collected. These data can be analysed with passive seismic methods and therefore offer the possibility to complement subsurface characterization at marginal additional cost. We present data and results from a 2.4 km long active source profile which has been recently acquired in Western Colorado (US) to characterize the structure and sedimentary infill of an over-deepened alpine valley. We show how the leftover passive data from the active source acquisition can be processed towards a shear wave velocity model with seismic interferometry. The shear wave velocity model supports the structural interpretation of the active P-wave data, and the P-to-S-wave velocity ratio provides new insights into the nature and hydrological properties of the sedimentary infill. We discuss the benefits and limitations of our workflow and conclude with recommendations for acquisition and processing of similar data sets.


Solid Earth ◽  
2019 ◽  
Vol 10 (4) ◽  
pp. 1337-1354 ◽  
Author(s):  
Michael Behm ◽  
Feng Cheng ◽  
Anna Patterson ◽  
Gerilyn S. Soreghan

Abstract. The advent of cable-free nodal arrays for conventional seismic reflection and refraction experiments is changing the acquisition style for active-source surveys. Instead of triggering short recording windows for each shot, the nodes are continuously recording over the entire acquisition period from the first to the last shot. The main benefit is a significant increase in geometrical and logistical flexibility. As a by-product, a significant amount of continuous data might also be collected. These data can be analyzed with passive seismic methods and therefore offer the possibility to complement subsurface characterization at marginal additional cost. We present data and results from a 2.4 km long active-source profile, which have recently been acquired in western Colorado (US) to characterize the structure and sedimentary infill of an over-deepened alpine valley. We show how the “leftover” passive data from the active-source acquisition can be processed towards a shear wave velocity model with seismic interferometry. The shear wave velocity model supports the structural interpretation of the active P-wave data, and the P-to-S-wave velocity ratio provides new insights into the nature and hydrological properties of the sedimentary infill. We discuss the benefits and limitations of our workflow and conclude with recommendations for the acquisition and processing of similar datasets.


2016 ◽  
Vol 59 (5) ◽  
Author(s):  
Luca Minarelli ◽  
Sara Amoroso ◽  
Gabriele Tarabusi ◽  
Marco Stefani ◽  
Gabriele Pulelli

The 2012 earthquakes sequence stroke a wide area of the alluvial plain in the Emilia-Romagna Region and triggered a new research interest on the role of the subsurface stratigraphic architecture and petrophysical property distribution in the modulation of the local seismic effects. Few direct shear wave velocity V_S data were however available below the depth of 50m. The only available V_S measurements were obtained from an anticline area, characterized by a reduced stratigraphic thickness and peculiar sedimentary facies, hardly representative of the majority of the alluvial plain subsurface. The study provides the first V_S profile available from middle-upper Quaternary successions deposited into a fast subsiding syncline area of the Apennine Foredeep Basin. The P-wave velocity V_P and the S-wave velocity V_S logs fill in the previous data gap on the geophysical parameters needed for the estimation of the local seismic response. Both V_P and V_S logs were continuously acquired to the depth of 265 m. The log records a velocity increase with depth, punctuated by sharp increases at some stratigraphic discordance surfaces. The value of 800 m/s that characterizes the “seismic bedrock”, as defined by the Italian building code [NTC 2008] was never reached at any depth. The investigated succession records a depositional evolution from deltaic-marine to alluvial plain conditions, punctuated by six glacio-eustatic depositional cycles, developed in Middle-Upper Quaternary times. The stratigraphic units described in the syncline log were correlated at a regional scale, with the thinner anticline succession of Mirandola. Correlatable units deposited into syncline and anticline areas reveal similar shear wave velocity values, supporting the regional extrapolation of the measured values.<br /><br /><br /><em></em>


1998 ◽  
Vol 41 (1) ◽  
Author(s):  
G. A. Tselentis ◽  
G. Delis

The importance of detailed knowledge of the shear-wave velocity structure of the upper geological layers was recently stressed in strong motion studies. In this work we describe an algorithm which we have developed to infer the 1D shear wave velocity structure from the inversion of multichannel surface wave dispersion data (ground-roll). Phase velocities are derived from wavenumber-frequency stacks while the inversion process is speeded up by the use of Householder transformations. Using synthetic and experimental data, we examined the applicability of the technique in deducing S-wave profiles. The comparison of the obtained results with those derived from cross-hole measurements and synthesized wave fields proved the reliability of the technique for the rapid assessment of shear wave profiles during microzonation investigations.


Geophysics ◽  
1996 ◽  
Vol 61 (2) ◽  
pp. 437-443 ◽  
Author(s):  
Ningya Cheng ◽  
Chuen Hon Cheng

Field data sets collected by an array monopole acoustic logging tool and a shear wave logging tool are processed and interpreted. The P‐ and S‐wave velocities of the formation are determined by threshold detection with cross‐correlation correction from the full waveform and the shear‐wave log, respectively. The array monopole acoustic logging data are also processed using the extended Prony’s method to estimate the borehole Stoneley wave phase velocity and attenuation as a function of frequency. The well formation between depths of 2950 and 3150 ft (899 and 960 m) can be described as an isotropic elastic medium. The inverted [Formula: see text] from the Stoneley wave phase velocity is in excellent agreement with the shear‐wave log results in this section. The well formation between the depths of 3715 and 3780 ft (1132 and 1152 m) can be described as a porous medium with shear‐wave velocity anisotropy about 10% to 20% and with the symmetry axis perpendicular to the borehole axis. The disagreement between the shear‐wave velocity from the Stoneley wave inversion and the direct shear‐wave log velocity in this section is beyond the errors in the measurements. Estimated permeabilities from low‐frequency Stoneley wave velocity and attenuation data are in good agreement with the core measurements. Also it is proven that the formation permeability is not the cause of the discrepancy. From the estimated “shear/pseudo‐Rayleigh” phase velocities in the array monopole log and the 3-D finite‐difference synthetics in the anisotropic formation, the discrepancy can best be explained as shear‐wave anisotropy.


Geophysics ◽  
2000 ◽  
Vol 65 (1) ◽  
pp. 35-45
Author(s):  
Jarrod C. Dunne ◽  
Greg Beresford ◽  
Brian L. N Kennett

We developed guidelines for building a detailed elastic depth model by using an elastic synthetic seismogram that matched both prestack and stacked marine seismic data from the Gippsland Basin (Australia). Recomputing this synthetic for systematic variations upon the depth model provided insight into how each part of the model affected the synthetic. This led to the identification of parameters in the depth model that have only a minor influence upon the synthetic and suggested methods for estimating the parameters that are important. The depth coverage of the logging run is of prime importance because highly reflective layering in the overburden can generate noise events that interfere with deeper events. A depth sampling interval of 1 m for the P-wave velocity model is a useful lower limit for modeling the transmission response and thus maintaining accuracy in the tie over a large time interval. The sea‐floor model has a strong influence on mode conversion and surface multiples and can be built using a checkshot survey or by testing different trend curves. When an S-wave velocity log is unavailable, it can be replaced using the P-wave velocity model and estimates of the Poisson ratio for each significant geological formation. Missing densities can be replaced using Gardner’s equation, although separate substitutions are required for layers known to have exceptionally high or low densities. Linear events in the elastic synthetic are sensitive to the choice of inelastic attenuation values in the water layer and sea‐floor sediments, while a simple inelastic attenuation model for the consolidated sediments is often adequate. The usefulness of a 1-D depth model is limited by misties resulting from complex 3-D structures and the validity of the measurements obtained in the logging run. The importance of such mis‐ties can be judged, and allowed for in an interpretation, by recomputing the elastic synthetic after perturbing the depth model to simulate the key uncertainties. Taking the next step beyond using simplistic modeling techniques requires extra effort to achieve a satisfactory tie to each part of a prestack seismic record. This is rewarded by the greater confidence that can then be held in the stacked synthetic tie and applications such as noise identification, data processing benchmarking, AVO analysis, and inversion.


Geophysics ◽  
2009 ◽  
Vol 74 (3) ◽  
pp. E135-E147 ◽  
Author(s):  
Gregor T. Baechle ◽  
Gregor P. Eberli ◽  
Ralf J. Weger ◽  
Jose Luis Massaferro

To assess saturation effects on acoustic properties in carbonates, we measure ultrasonic velocity on 38 limestone samples whose porosity ranges from 5% to 30% under dry and water-saturated conditions. Complete saturation of the pore space with water causes an increase and decrease in compressional- and shear-wave velocity as well as significant changes in the shear moduli. Compressional velocities of most water-saturated samples are up to [Formula: see text] higher than the velocities of the dry samples. Some show no change, and a few even show a decrease in velocity. Shear-wave velocity [Formula: see text] generally decreases, but nine samples show an increase of up to [Formula: see text]. Water saturation decreases the shear modulus by up to [Formula: see text] in some samples and increases it by up to [Formula: see text] in others. The average increase in the shear modulus with water saturation is [Formula: see text]; the average decrease is [Formula: see text]. The [Formula: see text] ratio shows an overall increase with water saturation. In particular, rocks displaying shear weakening have distinctly higher [Formula: see text] ratios. Grainstone samples with high amounts of microporosity and interparticle macro-pores preferentially show shear weakening, whereas recrystallized limestones are prone to increase shear strengths with water saturation. The observed shear weakening indicates that a rock-fluid interaction occurs with water saturation, which violates one of the assumptions in Gassmann’s theory. We find a positive correlation between changes in shear modulus and the inability of Gassmann’s theory to predict velocities of water-saturated samples at high frequencies. Velocities of water-saturated samples predicted by Gassmann’s equation often exceed measured values by as much as [Formula: see text] for samples exhibiting shear weakening. In samples showing shear strengthening, Gassmann-predicted velocity values are as much as [Formula: see text] lower than measured values. In 66% of samples, Gassmann-predicted velocities show a misfit to measured water-saturated P-wave velocities. This discrepancy between measured and Gassmann-predicted velocity is not caused solely by velocity dispersion but also by rock-fluid interaction related to the pore structure of carbonates. Thus, a pore analysis should be conducted to assess shear-moduli changes and the resultant uncertainty for amplitude variation with offset analyses and velocity prediction using Gassmann’s theory.


1974 ◽  
Vol 64 (2) ◽  
pp. 355-374
Author(s):  
K. L. Kaila ◽  
V. G. Krishna ◽  
Hari Narain

abstract The upper mantle shear-wave velocity structure in the Japan region has been determined from S travel times of 101 earthquakes with focal depths varying from 40 to 600 km, using a new analytical method given by Kaila (1969). In southwestern Japan, the S velocity obtained as 4.35 km/sec at a 40-km depth remains almost constant to a depth of about 170 km. The shear velocity in northeastern Japan increases linearly from 4.42 km/sec at a depth of 45 km to 4.62 km/sec at a depth of 145-km. For central Japan, the S velocity determined as 4.41 km/sec at a 40-km depth increases linearly to 4.55 km/sec at a 170-km depth, followed by a slight decrease in the velocity gradient, with velocity still increasing linearly to 4.68 km/sec at a depth of 345 km. At this transition depth, there is a first-order velocity discontinuity, the velocity increasing from 4.68 to 4.92 km/sec. Below this depth, velocity again increases linearly from 4.92 to 5.04 km/sec at a depth of 600 km. The shear velocities at depths between 440 to 640 km in Japan are found to be extremely low in comparison to those of Jeffreys (1939), Gutenberg (1959), and Arnold (1967). These low S velocities can explain satisfactorily the late S arrivals from shallow earthquakes between Δ = 20° to 30° as observed in the Japanese region. Graphs have been drawn to show the variation with depth of Δ*, the epicentral distance to the inflection point, Δ1, Δ2, (Δ2 − Δ1), ptrue =∂T/∂Δ, and αS = (T − pΔ) at the inflection point as obtained from the S-wave travel-time analysis.


Author(s):  
Jiayan Tan ◽  
Charles A. Langston ◽  
Sidao Ni

ABSTRACT Ambient noise cross-correlations, used to obtain fundamental-mode Rayleigh-wave group velocity estimates, and teleseismic P-wave receiver functions are jointly modeled to obtain a 3D shear-wave velocity model for the crust and upper mantle of Oklahoma. Broadband data from 82 stations of EarthScope Transportable Array, the U.S. National Seismic Network, and the Oklahoma Geological Survey are used. The period range for surface-wave ambient noise Green’s functions is from 4.5 to 30.5 s constraining shear-wave velocity to a depth of 50 km. We also compute high-frequency receiver functions at these stations from 214 teleseismic earthquakes to constrain individual 1D velocity models inferred from the surface-wave tomography. Receiver functions reveal Ps conversions from the Moho, intracrustal interfaces, and shallow sedimentary basins. Shallow low-velocity zones in the model correlate with the large sedimentary basins of Oklahoma. The velocity model significantly improves the agreement of synthetic and observed seismograms from the 6 November 2011 Mw 5.7 Prague, Oklahoma earthquake suggesting that it can be used to improve earthquake location and moment tensor inversion of local and regional earthquakes.


Geophysics ◽  
2014 ◽  
Vol 79 (4) ◽  
pp. T233-T241 ◽  
Author(s):  
Christopher S. Sherman ◽  
James Rector ◽  
Steven Glaser

The Born and Rytov approximation, radiative transfer theory, and other related techniques are commonly used to model features of wave propagation through heterogeneous geologic media such as scattering, attenuation, and pulse-broadening. However, due to the underlying assumptions about the scattering direction and the reference Green’s function, these methods overlook important features of the wavefield such as mode conversion and near-field term coupling. These effects are particularly important within the predicted S-wave nodes of a seismic source, so we analyzed the problem of wave propagation beneath a vertical-point force on the surface of a heterogeneous, elastic half space. To do this, we generated a suite of 3D synthetic heterogeneous geologic models using fractal statistics and simulated the wave propagation using the finite-difference method. We derived an estimate for the effective source radiation patterns, and we used these to compare the results of the models. Our numerical results showed that, due to a combination of mode conversion and near-source coupling effects, S-wave energy on the order of 10% of the P-wave energy is generated within the shear-radiation node. In some cases, this S-wave energy may occur as a coherent pulse and may be used to enhance seismic imaging.


Sign in / Sign up

Export Citation Format

Share Document