scholarly journals Subsonic near-surface P-velocity and low S-velocity observations using propagator inversion

Geophysics ◽  
2005 ◽  
Vol 70 (4) ◽  
pp. R15-R23 ◽  
Author(s):  
Robbert van Vossen ◽  
Andrew Curtis ◽  
Jeannot Trampert

Detailed knowledge of near-surface P- and S-wave velocities is important for processing and interpreting multicomponent land seismic data because (1) the entire wavefield passes through and is influenced by the near-surface soil conditions, (2) both source repeatability and receiver coupling also depend on these conditions, and (3) near-surface P- and S-wave velocities are required for wavefield decomposition and demultiple methods. However, it is often difficult to measure these velocities with conventional techniques because sensitivity to shallow-wave velocities is low and because of the presence of sharp velocity contrasts or gradients close to the earth's free surface. We demonstrate that these near-surface P- and S-wave velocities can be obtained using a propagator inversion. This approach requires data recorded by at least one multicomponent geophone at the surface and an additional multicomponent geophone at depth. The propagator between them then contains all information on the medium parameters governing wave propagation between the geophones at the surface and at depth. Hence, inverting the propagator gives local estimates for these parameters. This technique has been applied to data acquired in Zeist, the Netherlands. The near-surface sediments at this site are unconsolidated sands with a thin vegetation soil on top, and the sediments considered are located above the groundwater table. A buried geophone was positioned 1.05 m beneath receivers on the surface. Propagator inversion yielded low near-surface velocities, namely, 270 ± 15 m/s for the compressional-wave velocity, which is well below the sound velocity in air, and 150 ± 9 m/s for the shear velocity. Existing methods designed for imaging deeper structures cannot resolve these shallow material properties. Furthermore, velocities usually increase rapidly with depth close to the earth's surface because of increasing confining pressure. We suspect that for this reason, subsonic near-surface P-wave velocities are not commonly observed.

1996 ◽  
Vol 86 (6) ◽  
pp. 1704-1713 ◽  
Author(s):  
R. D. Catchings ◽  
W. H. K. Lee

Abstract The 17 January 1994, Northridge, California, earthquake produced strong ground shaking at the Cedar Hills Nursery (referred to here as the Tarzana site) within the city of Tarzana, California, approximately 6 km from the epicenter of the mainshock. Although the Tarzana site is on a hill and is a rock site, accelerations of approximately 1.78 g horizontally and 1.2 g vertically at the Tarzana site are among the highest ever instrumentally recorded for an earthquake. To investigate possible site effects at the Tarzana site, we used explosive-source seismic refraction data to determine the shallow (<70 m) P-and S-wave velocity structure. Our seismic velocity models for the Tarzana site indicate that the local velocity structure may have contributed significantly to the observed shaking. P-wave velocities range from 0.9 to 1.65 km/sec, and S-wave velocities range from 0.20 and 0.6 km/sec for the upper 70 m. We also found evidence for a local S-wave low-velocity zone (LVZ) beneath the top of the hill. The LVZ underlies a CDMG strong-motion recording site at depths between 25 and 60 m below ground surface (BGS). Our velocity model is consistent with the near-surface (<30 m) P- and S-wave velocities and Poisson's ratios measured in a nearby (<30 m) borehole. High Poisson's ratios (0.477 to 0.494) and S-wave attenuation within the LVZ suggest that the LVZ may be composed of highly saturated shales of the Modelo Formation. Because the lateral dimensions of the LVZ approximately correspond to the areas of strongest shaking, we suggest that the highly saturated zone may have contributed to localized strong shaking. Rock sites are generally considered to be ideal locations for site response in urban areas; however, localized, highly saturated rock sites may be a hazard in urban areas that requires further investigation.


Geophysics ◽  
1990 ◽  
Vol 55 (4) ◽  
pp. 470-479 ◽  
Author(s):  
D. F. Winterstein ◽  
B. N. P. Paulsson

Crosshole and vertical seismic profile (VST) data made possible accurate characterization of the elastic properties, including noticeable velocity anisotropy, of a near‐surface late Tertiary shale formation. Shear‐wave splitting was obvious in both crosshole and VSP data. In crosshole data, two orthologonally polarrized shear (S) waves arrived 19 ms in the uppermost 246 ft (75 m). Vertically traveling S waves of the VSP separated about 10 ms in the uppermost 300 ft (90 m) but remained at nearly constant separation below that level. A transversely isotropic model, which incorporates a rapid increase in S-wave velocities with depth but slow increase in P-wave velocities, closely fits the data over most of the measured interval. Elastic constants of the transvesely isotropic model show spherical P- and [Formula: see text]wave velocity surfaces but an ellipsoidal [Formula: see text]wave surface with a ratio of major to minor axes of 1.15. The magnitude of this S-wave anisotropy is consistent with and lends credence to S-wave anisotropy magnitudes deduced less directly from data of many sedimentary basins.


Geophysics ◽  
1995 ◽  
Vol 60 (6) ◽  
pp. 1627-1633 ◽  
Author(s):  
Bart W. Tichelaar ◽  
Klaas W. van Luik

Borehole sonic waveforms are commonly acquired to produce logs of subsurface compressional and shear wave velocities. To this purpose, modern borehole sonic tools are usually equipped with various types of acoustic sources, i.e., monopole and dipole sources. While the dipole source has been specifically developed for measuring shear wave velocities, we found that the dipole source has an advantage over the monopole source when determining compressional wave velocities in a very slow formation consisting of unconsolidated sands with a porosity of about 35% and a shear wave velocity of about 465 m/s. In this formation, the recorded compressional refracted waves suffer from interference with another wavefield component identified as a leaky P‐wave, which hampers the determination of compressional wave velocities in the sands. For the dipole source, separation of the compressional refracted wave from the recorded waveforms is accomplished through bandpass filtering since the wavefield components appear as two distinctly separate contributions to the frequency spectrum: a compressional refracted wave centered at a frequency of 6.5 kHz and a leaky P‐wave centered at 1.3 kHz. For the monopole source, the frequency spectra of the various waveform components have considerable overlap. It is therefore not obvious what passband to choose to separate the compressional refracted wave from the monopole waveforms. The compressional wave velocity obtained for the sands from the dipole compressional refracted wave is about 2150 m/s. Phase velocities obtained for the dispersive leaky P‐wave excited by the dipole source range from 1800 m/s at 1.0 kHz to 1630 m/s at 1.6 kHz. It appears that the dipole source has an advantage over the monopole source for the data recorded in this very slow formation when separating the compressional refracted wave from the recorded waveforms to determine formation compressional wave velocities.


Geophysics ◽  
2002 ◽  
Vol 67 (2) ◽  
pp. 405-412 ◽  
Author(s):  
Manika Prasad

Shallow water flows and over‐pressured zones are a major hazard in deepwater drilling projects. Their detection prior to drilling would save millions of dollars in lost drilling costs. I have investigated the sensitivity of seismic methods for this purpose. Using P‐wave information alone can be ambiguous, because a drop in P‐wave velocity (Vp) can be caused both by overpressure and by presence of gas. The ratio of P‐wave velocity to S‐wave velocity (Vp/Vs), which increases with overpressure and decreases with gas saturation, can help differentiate between the two cases. Since P‐wave velocity in a suspension is slightly below that of the suspending fluid and Vs=0, Vp/Vs and Poisson's ratio must increase exponentially as a load‐bearing sediment approaches a state of suspension. On the other hand, presence of gas will also decrease Vp but Vs will remain unaffected and Vp/Vs will decrease. Analyses of ultrasonic P‐ and S‐wave velocities in sands show that the Vp/Vs ratio, especially at low effective pressures, decreases rapidly with pressure. At very low pressures, Vp/Vs values can be as large as 100 and higher. Above pressures greater than 2 MPa, it plateaus and does not change much with pressure. There is significant change in signal amplitudes and frequency of shear waves below 1 MPa. The current ultrasonic data shows that Vp/Vs values can be invaluable indicators of low differential pressures.


Geophysics ◽  
2022 ◽  
pp. 1-49
Author(s):  
Yu-Tai Wu

Beishih Village of Hsinhua Township in southern Taiwan is a unique location for studying soil liquefaction. Soil liquefaction was observed at the same site after earthquakes in 1946, 2010, and 2016, each of which had a Richter magnitude greater than six. This recurrence provides an opportunity for analyzing soil condition variations resulting from soil liquefaction. Seismic data sets were collected in 2011, 2014, 2016, and 2017. We used seismic refraction tomography and the multichannel analysis of surface waves to estimate P- and S-wave velocities. In S-wave velocity profiles, low shear velocity zones were located beneath sand volcanoes shortly after two earthquakes and disappeared 4 years after a 2010 earthquake. However, the P-wave velocity is less sensitive to soil condition changes, possibly because groundwater obscures the effect of soil liquefaction on velocity profiles. In addition, we used seismic wave velocities to determine the importance of soil properties such as Poisson’s ratio, shear modulus, and porosity to identify the cause of the low shear velocity zone. Notably, although porosity decreased after soil grain rearrangement, sand and clay mixing increased the Poisson’s ratio, reducing the shear modulus of the soil. In addition, a soil layer between 2 and 7 m and a deeper layer below 10 m that resulted in sand volcanoes were both liquefied. We also considered how the evaluation of soil liquefaction potential could be affected by long-term variations in soil conditions and changes resulting from liquefaction. The factor of safety was used to evaluate the liquefaction potential of the site. The results revealed that the assessment conducted long after the earthquake underestimated risk because the soil developed shear strength after the earthquake.


Geophysics ◽  
2007 ◽  
Vol 72 (1) ◽  
pp. E1-E13 ◽  
Author(s):  
Michael A. Zimmer ◽  
Manika Prasad ◽  
Gary Mavko ◽  
Amos Nur

Knowledge of the pressure dependences of seismic velocities in unconsolidated sands is necessary for the remote prediction of effective pressures and for the projection of velocities to unsampled locations within shallow sand layers. We have measured the compressional- and shear-wave velocities and bulk, shear, and P-wave moduli at pressures from [Formula: see text] in a series of unconsolidated granular samples including dry and water-saturated natural sands and dry synthetic sand and glass-bead samples. The shear-wave velocities in these samples demonstrate an average pressure dependence approximately proportional to the fourth root of the effective pressure [Formula: see text], as commonly observed at lower pressures. For the compressional-wave velocities, theexponent in the pressure dependence of individual dry samples is consistently less than the exponent for the shear-wave velocity of the same sample, averaging 0.23 for the dry sands and 0.20 for the glass-bead samples. These pressure dependences are generally consistent over the entire pressure range measured. A comparison of the empirical results to theoretical predictions based on Hertz-Mindlin effective-medium models demonstrates that the theoretical models vastly overpredict the shear moduli of the dry granular frame unless the contacts are assumed to have no tangential stiffness. The models also predict a lower pressure exponent for the moduli and velocities [Formula: see text] than is generally observed in the data. We attribute this discrepancy in part to the inability of the models to account for decreases in the amount of slip or grain rotation occurring at grain-to-grain contacts with increasing pressure.


1989 ◽  
Vol 5 (3) ◽  
pp. 513-537 ◽  
Author(s):  
James F. Gibbs ◽  
Edward F. Roth

Shear (S)- and compressional (P)- wave velocities were measured to a depth of 195 m in a borehole near the San Andreas fault where a recurrence of a moderate Parkfield earthquake is predicted. S-wave velocities determined from orthogonal directions of the S-wave source show velocity differences of approximately 20 percent. An average shear-wave Q of 4 was determined in relatively unconsolidated sands and gravels of the Paso Robles Formation in the depth interval 57.5-102.5 m.


1982 ◽  
Vol 19 (8) ◽  
pp. 1535-1547 ◽  
Author(s):  
C. Wright

Seismological experiments have been undertaken at a test site near Chalk River, Ontario that consists of crystalline rocks covered by glacial sediments. Near-surface P and S wave velocity and amplitude variations have been measured along profiles less than 2 km in length. The P and S wave velocities were generally in the range 4.5–5.6 and 2.9–3.2 km/s, respectively. These results are consistent with propagation through fractured gneiss and monzonite, which form the bulk of the rock body. The P wave velocity falls below 5.0 km/s in a region where there is a major fault and in an area of high electrical conductivity; such velocity minima are therefore associated with fracture systems. For some paths, the P and 5 wave velocities were in the ranges 6.2–6.6 and 3.7–4.1 km/s, respectively, showing the presence of thin sheets of gabbro. Temporal changes in P travel times of up to 1.4% over a 12 h period were observed where the sediment cover was thickest. The cause may be changes in the water table. The absence of polarized SH arrivals from specially designed shear wave sources indicates the inhomogeneity of the test site. A Q value of 243 ± 53 for P waves was derived over one relatively homogeneous profile of about 600 m length. P wave velocity minima measured between depths of 25 and 250 m in a borehole correlate well with the distribution of fractures inferred from optical examination of borehole cores, laboratory measurements of seismic velocities, and tube wave studies.


Geophysics ◽  
1994 ◽  
Vol 59 (1) ◽  
pp. 131-139 ◽  
Author(s):  
M. Boulfoul ◽  
D. R. Watts

Instantaneous rotations are combined with f-k filtering to extract coherent S‐wave events from multicomponent shot records recorded by British Institutions Reflection Profiling Syndicate (BIRPS) Weardale Integrated S‐wave and P‐wave analysis (WISPA) experiment. This experiment was an attempt to measure the Poisson’s ratio of the lower crest by measuring P‐wave and S‐wave velocities. The multihole explosive source technique did generate S‐waves although not of opposite polarization. Attempts to produce stacks of the S‐wave data are unsuccessful because S‐wave splitting in the near surface produced random polarizations from receiver group to receiver group. The delay between the split wavelets varies but is commonly between 20 to 40 ms for 10 Hz wavelets. Dix hyperbola are produced on shot records after instantaneous rotations are followed by f-k filtering. To extract the instantaneous polarization, the traces are shifted back by the length of a moving window over which the calculation is performed. The instantaneous polarization direction is computed from the shifted data using the maximum eigenvector of the covariance matrix over the computation window. Split S‐waves are separated by the instantaneous rotation of the unshifted traces to the directions of the maximum eigenvectors determined for each position of the moving window. F-K filtering is required because of the presence of mode converted S‐waves and S‐waves produced by the explosive source near the time of detonation. Examples from synthetic data show that the method of instantaneous rotations will completely separate split S‐waves if the length of the moving window over which the calculation is performed is the length of the combined split wavelets. Separation may be achieved on synthetic data for wavelet delays as small as two sample intervals.


Geophysics ◽  
2018 ◽  
Vol 83 (6) ◽  
pp. B339-B353 ◽  
Author(s):  
Bojan Brodic ◽  
Alireza Malehmir ◽  
André Pugin ◽  
Georgiana Maries

We deployed a newly developed 3C microelectromechanical system-based seismic land streamer over porous glacial sediments to delineate water table and bedrock in Southwestern Finland. The seismic source used was a 500 kg vertical impact drop hammer. We analyzed the SH-wave component and interpreted it together with previously analyzed P-wave component data. In addition to this, we examined the land streamer’s potential for multichannel analysis of surface waves and delineated the site’s stratigraphy with surface-wave-derived S-wave velocities and [Formula: see text] ratios along the entire profile. These S-wave velocities and [Formula: see text] ratios complement the interpretation conducted previously on P-wave stacked section. Peculiarly, although the seismic source used is of a vertical-type nature, the data inspection indicated clear bedrock reflection on the horizontal components, particularly the transverse component. This observation led us to scrutinize the horizontal component data through side-by-side inspection of the shot records of all the three components and particle motion analysis to confirm the S-wave nature of the reflection. Using the apparent moveout velocity of the reflection, as well as the known depth to bedrock based on drilling, we used finite-difference synthetic modeling to further verify its nature. Compared with the P-wave seismic section, bedrock is relatively well delineated on the transverse component S-wave section. Some structures connected to the kettle holes and other stratigraphic units imaged on the P-wave results were also notable on the S-wave section, and particularly on the surface-wave derived S-wave velocity model and [Formula: see text] ratios. Our results indicate that P-, SV-, and SH-wave energy is generated simultaneously at the source location itself. This study demonstrates the potential of 3C seismic for characterization and delineation of the near-surface seismics.


Sign in / Sign up

Export Citation Format

Share Document