Surface and near-surface effects on seismic waves—theory and borehole seismometer results

1987 ◽  
Vol 77 (4) ◽  
pp. 1168-1196
Author(s):  
P. M. Shearer ◽  
J. A. Orcutt

Abstract A simple plane wave model is adequate to explain many surface versus borehole seismometer data sets. Using such a model, we present a series of examples which demonstrate the effects of the free-surface, near-surface velocity gradients, and low impedance surface layers on the amplitudes of upcoming body waves. In some cases, these amplitudes are predictable from simple free-surface and impedance contrast expressions. However, in other cases these expressions are an unreliable guide to the complete response, and the full plane wave calculation must be performed. Large surface amplifications are possible, even without focusing due to lateral heterogeneities or nonlinear effects. Both surface and borehole seismometer site responses are almost always frequency-dependent. Ocean bottom versus borehole seismic data from the 1983 Ngendei Seismic Experiment in the southwest Pacific are consistent with both a simple plane wave model and a more complete synthetic seismogram calculation. The borehole seismic response to upcoming P waves is reduced at high frequencies because of interference between the upgoing P wave and downgoing P and SV waves reflected from the sediment-basement interface. However, because of much lower borehole noise levels, the borehole seismometer enjoys a P-wave signal-to-noise advantage of 3 to 7 dB over nearby ocean bottom instruments.

1990 ◽  
Vol 80 (6B) ◽  
pp. 2032-2052 ◽  
Author(s):  
D. C. Jepsen ◽  
B. L. N. Kennett

Abstract Both phased array techniques for single-component sensors and vectorial analysis of three-component recordings can provide estimates of the azimuth and slowness of seismic phases. However, a combination of these approaches provides a more powerful tool to estimate the propagation characteristics of different seismic phases at regional distances. Conventional approaches to the analysis of three-component seismic records endeavor to exploit the apparent angles of propagation in horizontal and vertical planes as well as the polarization of the waves. The basic assumption is that for a given time window there is a dominant wavetype (e.g., a P wave) traveling in a particular direction arriving at the seismic station. By testing a range of characteristics of the three-component records, a set of rules can be established for classifying much of the seismogram in terms of wavetype and direction. It is, however, difficult to recognize SH waves in the presence of other wavetypes. Problems also arise when more than one signal (in either wavetype or direction) arrive in the same window. The stability and robustness of the classification scheme is much improved when records from an array of three-component sensors are combined. For a set of three-component instruments forming part of a larger array, it is possible to estimate the slowness and azimuth of arrivals from the main array and then extract the relative proportions of the current P-, SV-, and SH-wave contributions to the seismogram. This form of wavetype decomposition depends on a model of near-surface propagation. A convenient choice for hard-rock sites is to include just the effect of the free surface, which generates a frequency-independent operation on the three-component seismograms and which is not very sensitive to surface velocities. This approach generates good estimates of the character of the S wavefield, because the phase distortion of SV induced by the free surface can be removed. The method has been successfully applied to regional seismograms recorded at the medium aperture Warramunga array in northern Australia, and the two small arrays NORESS and ARCESS in Norway, which were designed for studies of regional phases. The new wavefield decomposition scheme provides results in which the relative proportions of P, SV, and SH waves as a function of time can be compared without the distortion imposed by free surface amplification. Such information can provide a useful adjunct to existing measures of signal character used in source discrimination.


2021 ◽  
Author(s):  
Nicola Piana Agostinetti ◽  
Alberto Villa ◽  
Gilberto Saccorotti

Abstract. We use PoroTOMO experimental data to compare the performance of Distributed Acoustic Sensing (DAS) and geophone data in executing standard exploration and monitoring activities. The PoroTOMO experiment consists of two "seismic systems": (a) a 8.6 km long optical fibre cable deployed across the Brady geothermal field and covering an area of 1.5 x 0.5 km with 100 m long segments, and (b) an array of 238 co-located geophones with an average spacing of 60 m. The PoroTOMO experiment recorded continuous seismic data between March 10th and March 25th 2016. During such period, a ML 4.3 regional event occurred in the southwest, about 150 km away from the geothermal field, together with several microseismic local events related to the geothermal activity. The seismic waves generated from such seismic events have been used as input data in this study. For the exploration tasks, we compare the propagation of the ML 4.3 event across the geothermal field in both seismic systems in term of relative time-delay, for a number of configurations and segments. Defined the propagation, we analyse and compare the amplitude and the signal-to-noise ratio (SNR) of the P-wave in the two systems at high resolution. For testing the potential in monitoring local seismicity, we first perform an analysis of the geophone data for locating a microseismic event, based on expert opinion. Then, we a adopt different workflow for the automatic location of the same microseismic event using DAS data. For testing the potential in monitoring distant event, data from the regional earthquake are used for retrieving both the propagation direction and apparent velocity of the wavefield, using a standard plane-wave-fitting approach. Our results indicate that: (1) at a local scale, the seismic P-waves propagation and their characteristics (i.e. SNR and amplitude) along a single cable segment are robustly consistent with recordings from co-located geophones (delay-times δt ∼ 0.3 over 400 m for both seismic systems) ; (2) the interpretation of seismic wave propagation across multiple separated segments is less clear, due to the heavy contamination of scattering sources and local velocity heterogeneities; nonetheless, results from the plane-wave fitting still indicate the possibility for a consistent detection and location of the event; (3) at high-resolution (10 m), large amplitude variations along the fibre cable seem to robustly correlate with near surface geology; (4) automatic monitoring of microseismicity can be performed with DAS recordings with results comparable to manual analysis of geophone recordings (i.e. maximum horizontal error on event location around 70 m for both geophones and DAS data) ; and (5) DAS data pre-conditioning (e.g., temporal sub-sampling and channel-stacking) and dedicated processing techniques are strictly necessary for making any real-time monitoring procedure feasible and trustable.


Geophysics ◽  
2013 ◽  
Vol 78 (1) ◽  
pp. U1-U8 ◽  
Author(s):  
Benoit de Cacqueray ◽  
Philippe Roux ◽  
Michel Campillo ◽  
Stefan Catheline

We tested a small-scale experiment that is dedicated to the study of the wave separation algorithm and to the velocity variations monitoring problem itself. It handles the case in which velocity variations at depth are hidden by near-surface velocity fluctuations. Using an acquisition system that combines an array of sources and an array of receivers, coupled with controlled velocity variations, we tested the ability of beam-forming techniques to track velocity variations separately for body waves and surface waves. After wave separation through double beam forming, the arrival time variations of the different waves were measured through the phase difference between the extracted wavelets. Finally, a method was tested to estimate near-surface velocity variations using surface waves or shallow reflection and compute a correction to isolate target velocity variations at depth.


2020 ◽  
Vol 221 (1) ◽  
pp. 683-691 ◽  
Author(s):  
F Brenguier ◽  
R Courbis ◽  
A Mordret ◽  
X Campman ◽  
P Boué ◽  
...  

SUMMARY Unveiling the mechanisms of earthquake and volcanic eruption preparation requires improving our ability to monitor the rock mass response to transient stress perturbations at depth. The standard passive monitoring seismic interferometry technique based on coda waves is robust but recovering accurate and properly localized P- and S-wave velocity temporal anomalies at depth is intrinsically limited by the complexity of scattered, diffracted waves. In order to mitigate this limitation, we propose a complementary, novel, passive seismic monitoring approach based on detecting weak temporal changes of velocities of ballistic waves recovered from seismic noise correlations. This new technique requires dense arrays of seismic sensors in order to circumvent the bias linked to the intrinsic high sensitivity of ballistic waves recovered from noise correlations to changes in the noise source properties. In this work we use a dense network of 417 seismometers in the Groningen area of the Netherlands, one of Europe's largest gas fields. Over the course of 1 month our results show a 1.5 per cent apparent velocity increase of the P wave refracted at the basement of the 700-m-thick sedimentary cover. We interpret this unexpected high value of velocity increase for the refracted wave as being induced by a loading effect associated with rainfall activity and possibly canal drainage at surface. We also observe a 0.25 per cent velocity decrease for the direct P-wave travelling in the near-surface sediments and conclude that it might be partially biased by changes in time in the noise source properties even though it appears to be consistent with complementary results based on ballistic surface waves presented in a companion paper and interpreted as a pore pressure diffusion effect following a strong rainfall episode. The perspective of applying this new technique to detect continuous localized variations of seismic velocity perturbations at a few kilometres depth paves the way for improved in situ earthquake, volcano and producing reservoir monitoring.


Geophysics ◽  
2021 ◽  
pp. 1-69
Author(s):  
Daniel Wehner ◽  
Filipe Borges ◽  
Martin Landrø

Monitoring the shear modulus of formations around boreholes is of interest for various applications, ranging from near-surface investigation to reservoir monitoring. Downhole logging tools and borehole seismic are common techniques applied to measure and characterize formation properties. These methods rely on transmitted and reflected waves to retrieve the rock properties. Wave modes travelling along the interface between the well and the formation, such as tube waves, are often considered as noise. However, tube waves are less attenuated than body waves, and contain information about the shear modulus of the formation surrounding the well. Hence, a potential use of this interface wave is of interest. As tube-wave properties depend on several parameters, e.g. well geometry, highly accurate measurements should be performed for use in inferring rock properties. We study the feasibility of tube-wave measurements as a monitoring method. Different experiments are conducted using a hydrophone array in two boreholes, with depths of 30 m and 95 m. The experiments are used to investigate how accurate the tube-wave velocity can be measured, and which parameters have most impact on the measurements. Our results suggest that it is hard to estimate the absolute shear modulus of the geological formation using tube-wave velocities only. However, it seems feasible to use them to monitor changes of the shear modulus, depending on the borehole set up and geological formation. The tube-wave monitoring can be used as a first step method to determine the depth along the well where changes occur before more accurate measurements are performed in a second step.


1990 ◽  
Vol 80 (5) ◽  
pp. 1245-1271 ◽  
Author(s):  
Y.-G. Li ◽  
P. C. Leary

Abstract Two instances of fault zone trapped seismic waves have been observed. At an active normal fault in crystalline rock near Oroville in northern California, trapped waves were excited with a surface source and recorded at five near-fault borehole depths with an oriented three-component borehole seismic sonde. At Parkfield, California, a borehole seismometer at Middle Mountain recorded at least two instances of the fundamental and first higher mode seismic waves of the San Andreas fault zone. At Oroville recorded particle motions indicate the presence of both Love and Rayleigh normal modes. The Love-wave dispersion relation derived for an idealized wave guide with velocity structure determined by body-wave travel-time inversion yields seismograms of the fundamental mode that are consistent with the observed Love-wave amplitude and frequency. Applying a similar velocity model to the Parkfield observations, we obtain a good fit to the trapped wavefield amplitude, frequency, dispersion, and mode time separation for an asymmetric San Andreas fault zone structure—a high-velocity half-space to the southwest, a low-velocity fault zone, a transition zone containing the borehole seismometer, and an intermediate velocity half-space to the northeast. In the Parkfield borehole seismic data set, the locations (depth and offset normal to fault plane) of natural seismic events which do or do not excite trapped waves are roughly consistent with our model of the low velocity zone. We conclude that it is feasible to obtain near-surface borehole records of fault zone trapped waves. Because trapped modes are excited only by events close to the fault zone proper—thereby fixing these events in space relative to the fault—a wider investigation of possible trapped mode waveforms recorded by a borehole seismic network could lead to a much refined body-wave/tomographic velocity model of the fault and to a weighting of events as a function of offset from the fault plane. It is an open question whether near-surface sensors exist in a stable enough seismic environment to use trapped modes as an earth monitoring device.


1990 ◽  
Vol 80 (4) ◽  
pp. 807-831 ◽  
Author(s):  
Jon B. Fletcher ◽  
Tom Fumal ◽  
Hsi-Ping Liu ◽  
Linda C. Carroll

Abstract To investigate near-surface site effects in granite rock, we drilled 300-m-deep boreholes at two sites which are collocated with stations from the digital array at Anza, California. The first borehole was sited at station KNW (Keenwild fire station), which is located along a ridge line about 8.7 km east of the San Jacinto Fault zone. Station PFO (Piñon Flat Observatory), chosen for the second site, is another 6 km further to the east of station KNW and is located on a gently sloping hillside. We logged each borehole for P- and S-wave velocities, as well as for crack density and orientation. P waves were generated by striking a plate with a hammer at the surface. A tool consisting of weighted anvils driven by compressed air against end plates along a 3.5-m beam was used to generate shear waves. Signals were recorded downhole with a three-component sensor package at 2.5-m intervals from the surface to 50 m depth, and at 5-m intervals from 50 m depth to the bottom of the hole. Velocities were determined by differencing the measured arrival times of first arrivals or peaks over each interval in depth. Travel times were computed for the first breaks at shallow depths, however, below about 100 m depth, times were computed for the first peaks rather than for first breaks since the first arrival was no longer clearly distinguishable. The KNW site yielded a shear velocity of 1.9 km/sec by only 30 m in depth and reached close to 2.6 km/sec at the bottom of the hole. P-wave velocities at KNW were also high at 5.4 km/sec starting at 120 m depth. The PFO site had similar but slightly higher shear-wave velocities. The bottom-hole shear-wave velocity reached 3.0 km/sec, and the P-wave velocity was 5.4 km/sec. Shear-wave attenuation was computed using both the pulse rise time and spectral ratio methods. At station KNW, attenuation was significant only in an interval between 17.5 and approximately 40 m in depth. Over the top 50 m, attenuation corresponding to a Q of about 8 was obtained. A total T* of 0.004 sec was measured for this interval. Pulse rise times also increased rapidly in this zone. The spectral ratio data for station PFO yields two peaks in attenuation above 50 m. Similar to the attenuation found for station KNW, the peak in attenuation corresponds to a Q of about 11, averaged over the top 50 m. Spectra of the seismic pulses produced by the hammer give good signal between 20 to 80 Hz. Significant motion perpendicular to the polarizations of the first shear-wave arrival was recorded within a few meters of the surface. Apparently, the rock structure is sufficiently complicated that body waves are being converted (SH to SV at oblique incidence) very close to the surface. The presence of these elliptical particle motions within a mere few m of the pure shear-wave source suggests that the detection of polarizations perpendicular to the main shear arrival at a single location at the surface is not, by itself, a good method for detecting shearwave splitting within the upper few tens of kilometers of the earth's crust. Crack densities and orientations were determined from televiewer records. These records showed cracks with a preferred direction at station KNW and of a greater density than at station PFO. At station PFO, crack densities were smaller and more diffuse in orientation.


1971 ◽  
Vol 61 (1) ◽  
pp. 109-145 ◽  
Author(s):  
J. R. Murphy ◽  
A. H. Davis ◽  
N. L. Weaver

abstract Frequency-dependent amplification of seismic waves by near-surface low-velocity layers is a well-known phenomenon. This phenomenon was examined from both the analytic and experimental viewpoints for body waves (P, SV, SH). Groundmotion data, recorded in conjunction with the underground nuclear testing program at the Nevada Test Site, are used to provide experimental validation of the analytic models. Experimental amplification factors are derived from these data for a variety of recording-site near-surface geological configurations (alluvium, mine tailings, fill) and a wide dynamic range of ground-motion intensity (10−5 to 100 g). The variability in the mean amplification observed at a site for repeated measurements is described statistically. This analysis shows that, although the amplification at a given site varies on the average by a factor of about 1.4 across the frequency band of interest, from detonation to detonation, the frequency and magnitude of the dominant amplification are fairly consistent. The quality of the comparisons of the observed and calculated amplification indicates that the available linear analytic models are capable of describing the major features of the frequency-dependent amplification observed for this wide range of groundmotion intensity and recording-site geology.


Geophysics ◽  
2018 ◽  
Vol 83 (4) ◽  
pp. R335-R344 ◽  
Author(s):  
Lu Liu ◽  
Yan Wu ◽  
Bowen Guo ◽  
Song Han ◽  
Yi Luo

Accurate estimation of near-surface velocity is a key step for imaging deeper targets. We have developed a new workflow to invert complex early arrivals in land seismic data for near-surface velocities. This workflow is composed of two methods: source-domain full traveltime inversion (FTI) and early arrival waveform inversion (EWI). Source-domain FTI automatically generates the background velocity that kinematically matches the reconstructed plane-wave sources from early arrivals with true plane-wave sources. This method does not require picking first arrivals for inversion, which is one of the most challenging and labor-intensive steps in ray-based first-arrival traveltime tomography, especially when the subsurface medium contains low-velocity zones that cause shingled multivalue arrivals. Moreover, unlike the conventional Born-based method, source-domain FTI can determine if the initial velocity is slower or faster than the true one according to the gradient sign. In addition, the computational cost is reduced considerably by using the one-way wave equation to extrapolate the plane-wave Green’s function. The source-domain FTI tomogram is then used as the starting model for EWI to obtain the short-wavelength component associated with the velocity model. We tested the workflow on two synthetic and one onshore filed data sets. The results demonstrate that source-domain FTI generates reasonable background velocities for EWI even though the first arrivals are shingled, and that this workflow can produce a high-resolution near-surface velocity model.


Sign in / Sign up

Export Citation Format

Share Document