RADIATION OF BODY WAVES FROM NEAR‐SURFACE EXPLOSIVE SOURCES

Geophysics ◽  
1966 ◽  
Vol 31 (6) ◽  
pp. 1057-1065 ◽  
Author(s):  
I. N. Gupta ◽  
C. Kisslinger

Amplitude distributions obtained from field observations of the azimuthal distribution of motion from cratering shots near a vertical face in a limestone section yielded data on radiation into a half‐space. These effects have been approximately reproduced in the laboratory by means of two‐dimensional seismic models. Small chemical explosions were fired on or near the edge of a large plexiglas sheet and the radiation of both P and S waves observed. Shots on the edge of the model sheet produce P and S radiation patterns expected from a normal downward impulse on the free surface. The radiation patterns from cratering shots may be qualitatively explained by the combined action on the free surface of a normal downward stress and a pair of horizontal stresses (dipole without moment) at the source point. The observed data are not sufficient for verifying theoretical S wave distributions. Observations of SV amplitudes from nuclear explosions could yield useful information concerning the relation between the angle at which the waves leave the source and the distance at which the wave emerges.

Geology ◽  
2019 ◽  
Vol 47 (12) ◽  
pp. 1176-1180 ◽  
Author(s):  
Erika Swanson ◽  
Aviva Sussman ◽  
Jennifer Wilson

Abstract Fractures within the earth control rock strength and fluid flow, but their dynamic nature is not well understood. As part of a series of underground chemical explosions in granite in Nevada, we collected and analyzed microfracture density data sets prior to, and following, individual explosions. Our work shows an ∼4-fold increase in both open and filled microfractures following the explosions. Based on the timing of core retrieval, filling of some new fractures occurs in as little as 6 wk after fracture opening under shallow (<100 m) crustal conditions. These results suggest that near-surface fractures may fill quite rapidly, potentially changing permeability on time scales relevant to oil, gas, and geothermal energy production; carbon sequestration; seismic cycles; and radionuclide migration from nuclear waste storage and underground nuclear explosions.


Author(s):  
Maryam Safarshahi ◽  
Igor B. Morozov

ABSTRACT Empirical models of geometrical-, Q-, t-star, and kappa-type attenuation of seismic waves and ground-motion prediction equations (GMPEs) are viewed as cases of a common empirical standard model describing variation of wave amplitudes with time and frequency. Compared with existing parametric and nonparametric approaches, several new features are included in this model: (1) flexible empirical parameterization with possible nonmonotonous time or distance dependencies; (2) joint inversion for time or distance and frequency dependencies, source spectra, site responses, kappas, and Q; (3) additional constraints removing spurious correlations of model parameters and data residuals with source–receiver distances and frequencies; (4) possible kappa terms for sources as well as for receivers; (5) orientation-independent horizontal- and three-component amplitudes; and (6) adaptive filtering to reduce noise effects. The approach is applied to local and regional S-wave amplitudes in southeastern Iran. Comparisons with previous studies show that conventional attenuation models often contain method-specific biases caused by limited parameterizations of frequency-independent amplitude decays and assumptions about the models, such as smoothness of amplitude variations. Without such assumptions, the frequency-independent spreading of S waves is much faster than inferred by conventional modeling. For example, transverse-component amplitudes decrease with travel time t as about t−1.8 at distances closer than 90 km and as t−2.5 beyond 115 km. The rapid amplitude decay at larger distances could be caused by scattering within the near surface. From about 90 to 115 km distances, the amplitude increases by a factor of about 3, which could be due to reflections from the Moho and within the crust. With more accurate geometrical-spreading and kappa models, the Q factor for the study area is frequency independent and exceeds 2000. The frequency-independent and Q-type attenuation for vertical-component and multicomponent amplitudes is somewhat weaker than for the horizontal components. These observations appear to be general and likely apply to other areas.


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.


1984 ◽  
Vol 74 (2) ◽  
pp. 361-376
Author(s):  
John Boatwright ◽  
Jon B. Fletcher

Abstract Seventy-three digitally recorded body waves from nine multiply recorded small earthquakes in Monticello, South Carolina, are analyzed to estimate the energy radiated in P and S waves. Assuming Qα = Qβ = 300, the body-wave spectra are corrected for attenuation in the frequency domain, and the velocity power spectra are integrated over frequency to estimate the radiated energy flux. Focal mechanisms determined for the events by fitting the observed displacement pulse areas are used to correct for the radiation patterns. Averaging the results from the nine events gives 27.3 ± 3.3 for the ratio of the S-wave energy to the P-wave energy using 0.5 〈Fi〉 as a lower bound for the radiation pattern corrections, and 23.7 ± 3.0 using no correction for the focal mechanisms. The average shift between the P-wave corner frequency and the S-wave corner frequency, 1.24 ± 0.22, gives the ratio 13.7 ± 7.3. The substantially higher values obtained from the integral technique implies that the P waves in this data set are depleted in energy relative to the S waves. Cursory inspection of the body-wave arrivals suggests that this enervation results from an anomalous site response at two of the stations. Using the ratio of the P-wave moments to the S-wave moments to correct the two integral estimates gives 16.7 and 14.4 for the ratio of the S-wave energy to the P-wave energy.


Geophysics ◽  
2011 ◽  
Vol 76 (3) ◽  
pp. S103-S113 ◽  
Author(s):  
Robert Sun ◽  
George A. McMechan ◽  
Han-Hsiang Chuang

The reflected P- and S-waves in elastic displacement component data recorded at the earth’s surface are separated by reverse-time (downward) extrapolation of the data in an elastic computational model, followed by calculations to give divergence (dilatation) and curl (rotation) at a selected reference depth. The surface data are then reconstructed by separate forward-time (upward) scalar extrapolations, from the reference depth, of the magnitude of the divergence and curl wavefields, and extraction of the separated P- and S-waves, respectively, at the top of the models. A P-wave amplitude will change by a factor that is inversely proportional to the P-velocity when it is transformed from displacement to divergence, and an S-wave amplitude will change by a factor that is inversely proportional to the S-velocity when it is transformed from displacement to curl. Consequently, the ratio of the P- to the S-wave amplitude (the P-S amplitude ratio) in the form of divergence and curl (postseparation) is different from that in the (preseparation) displacement form. This distortion can be eliminated by multiplying the separated S-wave (curl) by a relative balancing factor (which is the S- to P-velocity ratio); thus, the postseparation P-S amplitude ratio can be returned to that in the preseparation data. The absolute P- and S-wave amplitudes are also recoverable by multiplying them by a factor that depends on frequency, on the P-velocity α, and on the unit of α and is location-dependent if the near-surface P-velocity is not constant.


Geophysics ◽  
2007 ◽  
Vol 72 (6) ◽  
pp. T67-T79 ◽  
Author(s):  
Thomas Bohlen ◽  
Ulrich Lorang ◽  
Wolfgang Rabbel ◽  
Christof Müller ◽  
Rüdiger Giese ◽  
...  

For safe tunnel excavation, it is important to predict lithologic and structural heterogeneities ahead of construction. Conventional tunnel seismic prediction systems utilize body waves (P- and S-waves) that are directly generated at the tunnel walls or near the cutter head of the tunnel boring machine (TBM). We propose a new prediction strategy that has been discovered by 3D elastic finite-difference (FD) modeling: Rayleigh waves arriving at the front face of the tunnel are converted into high-amplitude S-waves propagating further ahead. Reflected or backscattered S-waves are converted back into Rayleigh waves which can be recorded along the sidewalls. We name these waves RSSR waves. In our approach, the front face acts as an S-wave transceiver. One technical advantage is that both the sources and the receivers may be placed behind the cutter head of the TBM. The modeling reveals that the RSSR waves exhibit significantly higher amplitudes than the directly reflected body waves. The excavation damage zone causes dispersion of the RSSR wave leading to multimodal reflection response. For the detection of geologic interfaces ahead, RSSR waves recorded along the sidewalls are corrected for dispersion and stacked. From the arrival times, the distance to the S-S reflection point can be estimated. A recurrent application, while the tunnel approaches the interface, allows one to quantify the orientation of the reflecting interfaces as well. Our approach has been verified successfully in a field experiment at the Piora adit of the Gotthard base tunnel. The distance to the Piora fault zone estimated from stacked RSSR events agrees well with the information obtained by geologic surveying and exploratory drilling.


2018 ◽  
Vol 29 ◽  
pp. 00019
Author(s):  
Katarzyna Hubicka ◽  
Jakub Sokolowski

Seismic event consists of surface waves and body waves. Due to the fact that the body waves are faster (P-waves) and more energetic (S-waves) in literature the problem of their analysis is taken more often. The most universal information that is received from the recorded wave is its moment of arrival. When this information is obtained from at least four seismometers in different locations, the epicentre of the particular event can be estimated [1]. Since the recorded body waves may overlap in signal, the problem of wave onset moment is considered more often for faster P-wave than S-wave. This however does not mean that the issue of S-wave arrival time is not taken at all. As the process of manual picking is time-consuming, methods of automatic detection are recommended (these however may be less accurate). In this paper four recently developed methods estimating S-wave arrival are compared: the method operating on empirical mode decomposition and Teager-Kaiser operator [2], the modification of STA/LTA algorithm [3], the method using a nearest neighbour-based approach [4] and the algorithm operating on characteristic of signals’ second moments. The methods will be also compared to wellknown algorithm based on the autoregressive model [5]. The algorithms will be tested in terms of their S-wave arrival identification accuracy on real data originating from International Research Institutions for Seismology (IRIS) database.


Geophysics ◽  
1984 ◽  
Vol 49 (6) ◽  
pp. 822-825 ◽  
Author(s):  
Walter L. Pilant

Plane‐wave directivity patterns for both P- and S-waves approaching a free surface are well known (Knopoff et al., 1957, Figure 3–5). These have been shown to apply in a reciprocal manner to time‐harmonic S-waves emanating from vertical and horizontal sources (Miller and Pursey, 1954; Cherry, 1962) in both two‐dimensional (2-D) and three‐dimensional (3-D) cases. Knopoff and Gilbert (1959) showed that the plane‐wave directivity patterns also apply to the first motions seen in the impulsive‐source case (3-D) and Pilant (1979, sec. 9–6) showed that they held in the equivalent 2-D problem. Theoretical expressions for these patterns are given by Pilant (ibid) as [Formula: see text] and [Formula: see text] where [Formula: see text] is measured from the vertical and the positive z-axis is into the medium. The x-axis lies along the free surface and the quantity [Formula: see text]. For angles greater than critical [Formula: see text], the proper expression for the square root is given by [Formula: see text] Thus for angles of incidence (or take‐off) greater than [Formula: see text], both [Formula: see text] and [Formula: see text] become complex numbers and lead to phase‐shift induced waveform changes as the S-waves interact with the free surface. The functions [Formula: see text] and [Formula: see text] are shown in Figure 1 for the angular range 34–37 degrees which includes the angle [Formula: see text] degrees. For this example, [Formula: see text] corresponding to a Poisson’s ratio equal to one‐quarter. The null in [Formula: see text] and the maximum in [Formula: see text] are clearly seen.


Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. D101-D116
Author(s):  
Julius K. von Ketelhodt ◽  
Musa S. D. Manzi ◽  
Raymond J. Durrheim ◽  
Thomas Fechner

Joint P- and S-wave measurements for tomographic cross-borehole analysis can offer more reliable interpretational insight concerning lithologic and geotechnical parameter variations compared with P-wave measurements on their own. However, anisotropy can have a large influence on S-wave measurements, with the S-wave splitting into two modes. We have developed an inversion for parameters of transversely isotropic with a vertical symmetry axis (VTI) media. Our inversion is based on the traveltime perturbation equation, using cross-gradient constraints to ensure structural similarity for the resulting VTI parameters. We first determine the inversion on a synthetic data set consisting of P-waves and vertically and horizontally polarized S-waves. Subsequently, we evaluate inversion results for a data set comprising jointly measured P-waves and vertically and horizontally polarized S-waves that were acquired in a near-surface ([Formula: see text]) aquifer environment (the Safira research site, Germany). The inverted models indicate that the anisotropy parameters [Formula: see text] and [Formula: see text] are close to zero, with no P-wave anisotropy present. A high [Formula: see text] ratio of up to nine causes considerable SV-wave anisotropy despite the low magnitudes for [Formula: see text] and [Formula: see text]. The SH-wave anisotropy parameter [Formula: see text] is estimated to be between 0.05 and 0.15 in the clay and lignite seams. The S-wave splitting is confirmed by polarization analysis prior to the inversion. The results suggest that S-wave anisotropy may be more severe than P-wave anisotropy in near-surface environments and should be taken into account when interpreting cross-borehole S-wave data.


Sign in / Sign up

Export Citation Format

Share Document