Near-surface seismic traveltime tomography using a direct-push source and surface-planted geophones

Geophysics ◽  
2009 ◽  
Vol 74 (4) ◽  
pp. G17-G25 ◽  
Author(s):  
Hendrik Paasche ◽  
Ulrike Werban ◽  
Peter Dietrich
Keyword(s):  
Velocity Distribution ◽  
Seismic Velocity ◽  
Velocity Model ◽  
Unconsolidated Sediments ◽  
Alluvial Deposits ◽  
Near Surface ◽  
Traveltime Tomography ◽  
Direct Push ◽  
Refraction Seismic ◽  
Push Technology

Information about seismic velocity distribution in heterogeneous near-surface sedimentary deposits is essential for a variety of environmental and engineering geophysical applications. We have evaluated the suitability of the minimally invasive direct-push technology for near-surface seismic traveltime tomography. Geophones placed at the surface and a seismic source installed temporarily in the subsurface by direct-push technology quickly acquire reversed multioffset vertical seismic profiles (VSPs). The first-arrival traveltimes of these data were used to reconstruct the 2D seismic velocity distribution tomographically. After testing this approach on synthetic data, we applied it to field data collected over alluvial deposits in a former river floodplain. The resulting velocity model contains information about high- and low-velocity anomalies and offers a significantly deeper penetration depth than conventional refraction tomography using surface-planted sources and receivers at the investigated site. A combination of refraction seismic and direct-push data increases resolution capabilities in the unsaturated zone and enables reliable reconstruction of velocity variations in near-surface unconsolidated sediments. The final velocity model structurally matches the results of cone-penetration tests and natural gamma-radiation data acquired along the profile. The suitability of multiple rapidly acquired reverse VSP surveys for 2D tomographic velocity imaging of near-surface unconsolidated sediments was explored.

Two robust imaging methodologies for challenging environments: Wave-equation dispersion inversion of surface waves and guided waves and supervirtual interferometry + tomography for far-offset refractions

2018 ◽  
Vol 6 (4) ◽  
pp. SM27-SM37 ◽  
Author(s):  
Jing Li ◽  
Kai Lu ◽  
Sherif Hanafy ◽  
Gerard Schuster
Keyword(s):  
Wave Equation ◽  
Velocity Distribution ◽  
Surface Waves ◽  
Guided Waves ◽  
Velocity Model ◽  
Dispersion Curves ◽  
P Wave ◽  
Head Wave ◽  
Near Surface ◽  
Velocity Models

Two robust imaging technologies are reviewed that provide subsurface geologic information in challenging environments. The first one is wave-equation dispersion (WD) inversion of surface waves and guided waves (GW) for the shear-velocity (S-wave) and compressional-velocity (P-wave) models, respectively. The other method is traveltime inversion for the velocity model, in which supervirtual refraction interferometry (SVI) is used to enhance the signal-to-noise ratio of far-offset refractions. We have determined the benefits and liabilities of both methods with synthetic seismograms and field data. The benefits of WD are that (1) there is no layered-medium assumption, as there is in conventional inversion of dispersion curves. This means that 2D or 3D velocity models can be accurately estimated from data recorded by seismic surveys over rugged topography, and (2) WD mostly avoids getting stuck in local minima. The liability is that WD for surface waves is almost as expensive as full-waveform inversion (FWI) and, for Rayleigh waves, only recovers the S-velocity distribution to a depth no deeper than approximately 1/2 to 1/3 wavelength of the lowest-frequency surface wave. The limitation for GW is that, for now, it can estimate the P-velocity model by inverting the dispersion curves from GW propagating in near-surface low-velocity zones. Also, WD often requires user intervention to pick reliable dispersion curves. For SVI, the offset of usable refractions can be more than doubled, so that traveltime tomography can be used to estimate a much deeper model of the P-velocity distribution. This can provide a more effective starting velocity model for FWI. The liability is that SVI assumes head-wave first arrivals, not those from strong diving waves.

Crosscoherence-based interferometry for the retrieval of first arrivals and subsequent tomographic imaging of differential weathering

Geophysics ◽  
2019 ◽  
Vol 84 (4) ◽  
pp. Q37-Q48 ◽  
Author(s):  
Joachim Place ◽  
Deyan Draganov ◽  
Alireza Malehmir ◽  
Christopher Juhlin ◽  
Chris Wijns
Keyword(s):  
Seismic Waves ◽  
Seismic Velocity ◽  
Signal To Noise Ratio ◽  
Seismic Profile ◽  
Seismic Interferometry ◽  
Near Surface ◽  
Tomographic Images ◽  
Traveltime Tomography ◽  
Seismic Methods ◽  
Source Signal

Exhumation of crust exposes rocks to weathering agents that weaken the rocks’ mechanical strength. Weakened rocks will have lower seismic velocity than intact rocks and can therefore be mapped using seismic methods. However, if the rocks are heavily weathered, they will attenuate controlled-source seismic waves to such a degree that the recorded wavefield would become dominated by ambient noise and/or surface waves. Therefore, we have examined the structure of differential weathering by first-break traveltime tomography over a seismic profile extending approximately 3.5 km and acquired at a mining site in Zambia using explosive sources and a source based on the swept-impact seismic technique (SIST). Seismic interferometry has been tested for the retrieval of supervirtual first arrivals masked by uncorrelated noise. However, use of crosscorrelation in the retrieval process makes the method vulnerable to changes in the source signal (explosives and SIST). Thus, we have developed a crosscoherence-based seismic-interferometry method to tackle this shortcoming. We investigate the method’s efficiency in retrieving first arrivals and, simultaneously, correctly handling variations in the source signal. Our results illustrate the superiority of the crosscoherence- over crosscorrelation-based method for retrieval of the first arrivals, especially in alleviating spurious ringyness and in terms of the signal-to-noise ratio. These benefits are observable in the greater penetration depth and the improved resolution of the tomography sections. The tomographic images indicate isolated bodies of higher velocities, which may be interpreted as fresh rocks embedded into a heavily weathered regolith, providing a conspicuous example of differential weathering. Our study advances the potential of seismic methods for providing better images of the near surface (the critical zone).

3D seismic traveltime tomography imaging of the shallow subsurface at the C O2 SINK project site, Ketzin, Germany

Geophysics ◽  
2009 ◽  
Vol 74 (1) ◽  
pp. G1-G15 ◽  
Author(s):  
Sawasdee Yordkayhun ◽  
Ari Tryggvason ◽  
Ben Norden ◽  
Christopher Juhlin ◽  
Björn Bergman
Keyword(s):  
Velocity Structure ◽  
Velocity Model ◽  
Geological Structure ◽  
Seismic Survey ◽  
Storage Site ◽  
Near Surface ◽  
Shallow Subsurface ◽  
Traveltime Tomography ◽  
Reflection Seismic ◽  
Reflection Data

A 3D reflection seismic survey was performed in 2005 at the Ketzin carbon dioxide [Formula: see text] pilot geological-storage site (the [Formula: see text] project) near Berlin, Germany, to image the geological structure of the site to depths of about [Formula: see text]. Because of the acquisition geometry, frequency limitations of the source, and artefacts of the data processing, detailed structures shallower than about [Formula: see text] were unclear. To obtain structural images of the shallow subsurface, we applied 3D traveltime tomography to data near the top of the Ketzin anticline, where faulting is present. Understanding the shallow subsurface structure is important for long-term monitoring aspects of the project after [Formula: see text] has been injected into a saline aquifer at about [Formula: see text] depth. We used a 3D traveltime tomography algorithm based on a combination ofsolving for 3D velocity structure and static corrections in the inversion process to account for artefacts in the velocity structure because of smearing effects from the unconsolidated cover. The resulting velocity model shows low velocities of [Formula: see text] in the uppermost shallow subsurface of the study area. The velocity reaches about [Formula: see text] at a depth of [Formula: see text]. This coincides approximately with the boundary between Quaternary units, which contain the near-surface freshwater reservoir and the Tertiary clay aquitard. Correlation of tomographic images with a similarity attribute slice at [Formula: see text] (about [Formula: see text] depth) indicates that at least one east-west striking fault zone observed in the reflection data might extend into the Tertiary unit. The more detailed images of the shallow subsurface from this study provided valuable information on this potentially risky area.

scholarly journals Ray-tracing traveltime tomography versus wave-equation traveltime inversion for near-surface seismic land data

2017 ◽  
Vol 5 (3) ◽  
pp. SO11-SO19
Author(s):  
Lei Fu ◽  
Sherif M. Hanafy
Keyword(s):  
Wave Equation ◽  
Ray Tracing ◽  
Water Table ◽  
Seismic Data ◽  
Velocity Model ◽  
Water Table Depth ◽  
Seismic Survey ◽  
Near Surface ◽  
Traveltime Tomography ◽  
Initial Starting

Full-waveform inversion of land seismic data tends to get stuck in a local minimum associated with the waveform misfit function. This problem can be partly mitigated by using an initial velocity model that is close to the true velocity model. This initial starting model can be obtained by inverting traveltimes with ray-tracing traveltime tomography (RT) or wave-equation traveltime (WT) inversion. We have found that WT can provide a more accurate tomogram than RT by inverting the first-arrival traveltimes, and empirical tests suggest that RT is more sensitive to the additive noise in the input data than WT. We present two examples of applying WT and RT to land seismic data acquired in western Saudi Arabia. One of the seismic experiments investigated the water-table depth, and the other one attempted to detect the location of a buried fault. The seismic land data were inverted by WT and RT to generate the P-velocity tomograms, from which we can clearly identify the water table depth along the seismic survey line in the first example and the fault location in the second example.

scholarly journals Characterization of seasonal glacial seismicity from a single-station on-ice record at Holtedahlfonna, Svalbard

2019 ◽  
Vol 60 (79) ◽  
pp. 23-36 ◽  
Author(s):  
Andreas Köhler ◽  
Valerie Maupin ◽  
Christopher Nuth ◽  
Ward van Pelt
Keyword(s):  
Seismic Velocity ◽  
Temporal Distribution ◽  
Velocity Model ◽  
Ice Dynamics ◽  
Near Surface ◽  
Waveform Modeling ◽  
Single Station ◽  
Glacier Runoff ◽  
Fast Flow ◽  
Seismic Background

ABSTRACTGlacial seismicity provides important insights into glacier dynamic processes. We study the temporal distribution of cryogenic seismic signals (icequakes) at Holtedahlfonna, Svalbard, between April and August 2016 using a single three-component sensor. We investigate sources of observed icequakes using polarization analysis and waveform modeling. Processes responsible for five icequake categories are suggested, incorporating observations of previous studies into our interpretation. We infer that the most dominant icequake type is generated by surface crevasse opening through hydrofracturing. Secondly, bursts of high-frequency signals are presumably caused by repeated near-surface crevassing due to high strain rates during glacier fast-flow episodes. Furthermore, signals related to resonance in water-filled cracks, fracturing or settling events in dry firn or snow before the melt season, and processes at the glacier bed are observed. Amplitude of seismic background noise is clearly related to glacier runoff. We process ambient seismic noise to invert horizontal-to-vertical spectral ratios for a sub-surface seismic velocity model used to model icequake signals. Our study shows that a single seismic sensor provides useful information about seasonal ice dynamics in case deployment of a network is not feasible.

scholarly journals Joint interpretation of magnetotelluric, seismic, and well-log data in Hontomín (Spain)

Solid Earth ◽  
2016 ◽  
Vol 7 (3) ◽  
pp. 943-958 ◽  
Author(s):  
Xènia Ogaya ◽  
Juan Alcalde ◽  
Ignacio Marzán ◽  
Juanjo Ledo ◽  
Pilar Queralt ◽  
...  
Keyword(s):  
Seismic Data ◽  
Velocity Model ◽  
Unconsolidated Sediments ◽  
Well Log ◽  
Storage Site ◽  
Near Surface ◽  
Log Data ◽  
Velocity Models ◽  
Joint Interpretation

Abstract. Hontomín (N of Spain) hosts the first Spanish CO2 storage pilot plant. The subsurface characterization of the site included the acquisition of a 3-D seismic reflection and a circumscribed 3-D magnetotelluric (MT) survey. This paper addresses the combination of the seismic and MT results, together with the available well-log data, in order to achieve a better characterization of the Hontomín subsurface. We compare the structural model obtained from the interpretation of the seismic data with the geoelectrical model resulting from the MT data. The models correlate well in the surroundings of the CO2 injection area with the major structural differences observed related to the presence of faults. The combination of the two methods allowed a more detailed characterization of the faults, defining their geometry, and fluid flow characteristics, which are key for the risk assessment of the storage site. Moreover, we use the well-log data of the existing wells to derive resistivity–velocity relationships for the subsurface and compute a 3-D velocity model of the site using the 3-D resistivity model as a reference. The derived velocity model is compared to both the predicted and logged velocity in the injection and monitoring wells, for an overall assessment of the computed resistivity–velocity relationships. The major differences observed are explained by the different resolution of the compared geophysical methods. Finally, the derived velocity model for the near surface is compared with the velocity model used for the static corrections in the seismic data. The results allowed extracting information about the characteristics of the shallow unconsolidated sediments, suggesting possible clay and water content variations. The good correlation of the velocity models derived from the resistivity–velocity relationships and the well-log data demonstrate the potential of the combination of the two methods for characterizing the subsurface, in terms of its physical properties (velocity, resistivity) and structural/reservoir characteristics. This work explores the compatibility of the seismic and magnetotelluric methods across scales highlighting the importance of joint interpretation in near surface and reservoir characterization.

Early arrival waveform tomography on near-surface refraction data

Geophysics ◽  
2006 ◽  
Vol 71 (4) ◽  
pp. U47-U57 ◽  
Author(s):  
Jianming Sheng ◽  
Alan Leeds ◽  
Maike Buddensiek ◽  
Gerard T. Schuster
Keyword(s):  
Ground Truth ◽  
Velocity Model ◽  
Frequency Method ◽  
Near Surface ◽  
Good Convergence ◽  
Traveltime Tomography ◽  
Early Arrival ◽  
Synthetic Test ◽  
Waveform Tomography ◽  
Tomography Method

We develop a waveform-tomography method for estimating the velocity distribution that minimizes the waveform misfit between the predicted and observed early arrivals in space-time seismograms. By fitting the waveforms of early arrivals, early arrival waveform tomography (EWT) naturally takes into account more general wave-propagation effects compared to the high-frequency method of traveltime tomography, meaning that EWT can estimate a wider range of slowness wavenumbers. Another benefit of EWT is more reliable convergence compared to full-waveform tomography, because an early-arrival misfit function contains fewer local minima. Synthetic test results verify that the waveform tomogram is much more accurate than the traveltime tomogram and that this algorithm has good convergence properties. For marine data from the Gulf of Mexico, the statics problem caused by shallow, gassy muds was attacked by using EWT to obtain a more accurate velocity model. Using the waveform tomogram to correct for statics, the stacked section was significantly improved compared to using the normal move-out (NMO) velocity, and moderately improved compared to using the traveltime tomogram. Inverting high-resolution land data from Mapleton, Utah, showed an EWT velocity tomogram that was more consistent with the ground truth (trench log) than the traveltime tomogram. Our results suggest that EWT can provide supplemental, shorter-wavelength information compared to the traveltime tomogram for both shallow and moderately deep seismic data.

scholarly journals A P-wave velocity model of the upper crust of the Sannio region (Southern Apennines, Italy)

1998 ◽  
Vol 41 (4) ◽  
Author(s):  
G. Iannaccone ◽  
L. Improta ◽  
P. Capuano ◽  
A. Zollo ◽  
G. Biella ◽  
...  
Keyword(s):  
Wave Velocity ◽  
Seismic Velocity ◽  
Carbonate Platform ◽  
Velocity Model ◽  
P Wave ◽  
Seismic Velocities ◽  
Upper Crust ◽  
Near Surface ◽  
P Wave Velocity ◽  
Apulia Platform

This paper describes the results of a seismic refraction profile conducted in October 1992 in the Sannio region, Southern Italy, to obtain a detailed P-wave velocity model of the upper crust. The profile, 75 km long, extended parallel to the Apenninic chain in a region frequently damaged in historical time by strong earthquakes. Six shots were fired at five sites and recorded by a number of seismic stations ranging from 41 to 71 with a spacing of 1-2 km along the recording line. We used a two-dimensional raytracing technique to model travel times and amplitudes of first and second arrivals. The obtained P-wave velocity model has a shallow structure with strong lateral variations in the southern portion of the profile. Near surface sediments of the Tertiary age are characterized by seismic velocities in the 3.0-4.1 km/s range. In the northern part of the profile these deposits overlie a layer with a velocity of 4.8 km/s that has been interpreted as a Mesozoic sedimentary succession. A high velocity body, corresponding to the limestones of the Western Carbonate Platform with a velocity of 6 km/s, characterizes the southernmost part of the profile at shallow depths. At a depth of about 4 km the model becomes laterally homogeneous showing a continuous layer with a thickness in the 3-4 km range and a velocity of 6 km/s corresponding to the Meso-Cenozoic limestone succession of the Apulia Carbonate Platform. This platform appears to be layered, as indicated by an increase in seismic velocity from 6 to 6.7 km/s at depths in the 6-8 km range, that has been interpreted as a lithological transition from limestones to Triassic dolomites and anhydrites of the Burano formation. A lower P-wave velocity of about 5.0-5.5 km/s is hypothesized at the bottom of the Apulia Platform at depths ranging from 10 km down to 12.5 km; these low velocities could be related to Permo-Triassic siliciclastic deposits of the Verrucano sequence drilled at the bottom of the Apulia Platform in the Apulia Foreland.

Sign in / Sign up

Export Citation Format

Share Document