scholarly journals Rayleigh-to-shear wave conversion at the tunnel face — From 3D-FD modeling to ahead-of-drill exploration

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.

Geophysics ◽  
2002 ◽  
Vol 67 (6) ◽  
pp. 1798-1809 ◽  
Author(s):  
Lorenzo Petronio ◽  
Flavio Poletto

The tunnel boring machine (TBM) is used extensively to mechanically excavate tunnels. To optimize the mechanical drilling and work safely, an estimate of the geology to be drilled is necessary. We consider using the elastic waves produced by the TBM cutting wheel to obtain seismic‐while‐drilling (SWD) information for predicting the geology ahead of the drilling front. This method uses accelerometers mounted on the TBM together with geophones located along and outside the tunnel, similar to the technique successfully used to drill oil and geothermal wells. Study of noise and the resolution of the signal produced by the large‐diameter cutting head shows that nonstationary noise separation can be achieved by locating sensors at the front and rear ends of the tunnel. The (higher) resolution in front of the TBM is limited by pilot delays, while the (lower) lateral resolution is limited by the radial dimension of the TBM. Analysis of seismic data acquired in a field test shows that P‐ and S‐wave arrivals have a wide frequency band and high amplitude in seismic traces measured 700 m away from the drilling front. In comparison with SWD applications in wells, tunnel SWD technology has the advantage of allowing direct access to the tunnel front, which makes it easy to connect the TBM reference sensors for while‐drilling monitoring. This method can be successfully applied without interfering with drilling activity to monitor tunnel excavation continuously, reduce risks, and optimize drilling.


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.


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.


1967 ◽  
Vol 57 (2) ◽  
pp. 191-198
Author(s):  
J. Cl. De Bremaecker

abstract Realistic artificial seismograms may be computed by considering body waves as sums of normal or leaking modes of surface waves: the S wave and those arriving after S may be considered as sums of higher normal modes of Rayleigh waves (RiN) and Love waves (LiN); in this case the apparent velocity c < βn. Earlier arrivals are generally due to the first kind of leaking modes of Rayleigh waves (RiL1) for which βn < c < αn. Deep reflections in seismic prospecting are RiL2 for which c > αn. Synthetic seismograms can be computed by double Fourier transformation in those two last cases. Alternately the method of residues followed by a single Fourier (or Laplace) transformation may be used in all cases. Earth-stretching approximations should give excellent results for SH waves and may give satisfactory results for P and SV waves.


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.


Geophysics ◽  
1969 ◽  
Vol 34 (5) ◽  
pp. 696-712 ◽  
Author(s):  
Ker C. Thomson ◽  
Thomas J. Ahrens ◽  
M. Nafi Toksöz

The occasional existence of very pronounced, anomalous, horizontally polarized seismic waves from underground nuclear bomb blasts has been reported by several investigators. In order to further understanding of this phenomenon and the processes of mechanical radiation from explosions, particularly in prestressed media, a model study has been undertaken. Experimental apparatus has been developed which permits the generation and propagation of body waves from explosions in transparent plate models prestressed to various two‐dimensional stress configurations. High‐speed framing camera sequences are presented showing the explosion process and the resulting plate compressional and shear wave propagation in prestressed models. These are compared to theoretical calculations of isochromatic and [Formula: see text] isoclinic fringe patterns associated with the wave propagation in stressfree plates and plates prestressed in tension and shear. The following distinctive optical phenomena were predicted theoretically and observed in the high‐speed photoelastic patterns: a [Formula: see text] discontinuity between P and S wave isoclinics for the unstressed case; a tendency for the isoclinics to broaden and envelope the isochromatics in regions where the P and S waves are superimposed; development of serrations in the dynamic isoclinics in the presence of a prestressing field (yielding a pseudo‐isochromatic appearance to isoclinics when viewed monochromatically); and finally, a general similarity between the dynamic optical effects in media under tensile and shear prestress.


Author(s):  
Hao Wang ◽  
Ning Li ◽  
Caizhi Wang ◽  
Hongliang Wu ◽  
Peng Liu ◽  
...  

Abstract In the process of dipole-source acoustic far-detection logging, the azimuth of the fracture outside the borehole can be determined with the assumption that the SH–SH wave is stronger than the SV–SV wave. However, in slow formations, the considerable borehole modulation highly complicates the dipole-source radiation of SH and SV waves. A 3D finite-difference time-domain method is used to investigate the responses of the dipole-source reflected shear wave (S–S) in slow formations and explain the relationships between the azimuth characteristics of the S–S wave and the source–receiver offset and the dip angle of the fracture outside the borehole. Results indicate that the SH–SH and SV–SV waves cannot be effectively distinguished by amplitude at some offset ranges under low- and high-fracture dip angle conditions, and the offset ranges are related to formation properties and fracture dip angle. In these cases, the fracture azimuth determined by the amplitude of the S–S wave not only has a $180^\circ $ uncertainty but may also have a $90^\circ $ difference from the actual value. Under these situations, the P–P, S–P and S–S waves can be combined to solve the problem of the $90^\circ $ difference in the azimuth determination of fractures outside the borehole, especially for a low-dip-angle fracture.


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.


2021 ◽  
Author(s):  
Itzhak Lior ◽  
Anthony Sladen ◽  
Diego Mercerat ◽  
Jean-Paul Ampuero ◽  
Diane Rivet ◽  
...  

<p>The use of Distributed Acoustic Sensing (DAS) presents unique advantages for earthquake monitoring compared with standard seismic networks: spatially dense measurements adapted for harsh environments and designed for remote operation. However, the ability to determine earthquake source parameters using DAS is yet to be fully established. In particular, resolving the magnitude and stress drop, is a fundamental objective for seismic monitoring and earthquake early warning. To apply existing methods for source parameter estimation to DAS signals, they must first be converted from strain to ground motions. This conversion can be achieved using the waves’ apparent phase velocity, which varies for different seismic phases ranging from fast body-waves to slow surface- and scattered-waves. To facilitate this conversion and improve its reliability, an algorithm for slowness determination is presented, based on the local slant-stack transform. This approach yields a unique slowness value at each time instance of a DAS time-series. The ability to convert strain-rate signals to ground accelerations is validated using simulated data and applied to several earthquakes recorded by dark fibers of three ocean-bottom telecommunication cables in the Mediterranean Sea. The conversion emphasizes fast body-waves compared to slow scattered-waves and ambient noise, and is robust even in the presence of correlated noise and varying wave propagation directions. Good agreement is found between source parameters determined using converted DAS waveforms and on-land seismometers for both P- and S-wave records. The demonstrated ability to resolve source parameters using P-waves on horizontal ocean-bottom fibers is key for the implementation of DAS based earthquake early warning, which will significantly improve hazard mitigation capabilities for offshore and tsunami earthquakes.</p>


1967 ◽  
Vol 57 (1) ◽  
pp. 55-81
Author(s):  
E. J. Douze

abstract This report consists of a summary of the studies conducted on the subject of short-period (6.0-0.3 sec period) noise over a period of approximately three years. Information from deep-hole and surface arrays was used in an attempt to determine the types of waves of which the noise is composed. The theoretical behavior of higher-mode Rayleigh waves and of body waves as measured by surface and deep-hole arrays is described. Both surface and body waves are shown to exist in the noise. Surface waves generally predominate at the longer periods (of the period range discussed) while body waves appear at the shorter periods at quiet sites. Not all the data could be interpreted to define the wave types present.


Sign in / Sign up

Export Citation Format

Share Document