scholarly journals Time-average velocity estimation through surface-wave analysis: Part 1 — S-wave velocity

Geophysics ◽  
2017 ◽  
Vol 82 (3) ◽  
pp. U49-U59 ◽  
Author(s):  
Laura Valentina Socco ◽  
Cesare Comina ◽  
Farbod Khosro Anjom
Keyword(s):  
Wave Velocity ◽  
Average Velocity ◽  
Field Data ◽  
Velocity Model ◽  
Velocity Estimation ◽  
S Wave ◽  
Simple Method ◽  
Near Surface ◽  
Velocity Models ◽  
Time Average

In some areas, the estimation of static corrections for land seismic data is a critical step of the processing workflow. It often requires the execution of additional surveys and data analyses. Surface waves (SWs) in seismic records can be processed to extract local dispersion curves (DCs) that can be used to estimate near-surface S-wave velocity models. Here we focus on the direct estimation of time-average S-wave velocity models from SW DCs without the need to invert the data. Time-average velocity directly provides the value of one-way time, given a datum plan depth. The method requires the knowledge of one 1D S-wave velocity model along the seismic line, together with the relevant DC, to estimate a relationship between SW wavelength and investigation depth on the time-average velocity model. This wavelength/depth relationship is then used to estimate all the other time-average S-wave velocity models along the line directly from the DCs by means of a data transformation. This approach removes the need for extensive data inversion and provides a simple method suitable for industrial workflows. We tested the method on synthetic and field data and found that it is possible to retrieve the time-average velocity models with uncertainties less than 10% in sites with laterally varying velocities. The error on one-way times at various depths of the datum plan retrieved by the time-average velocity models is mostly less than 5 ms for synthetic and field data.

Computing near-surface velocity models for S-wave static corrections using raypath interferometry

Geophysics ◽  
2018 ◽  
Vol 83 (3) ◽  
pp. U23-U34
Author(s):  
Raul Cova ◽  
David Henley ◽  
Kristopher A. Innanen
Keyword(s):  
Wave Velocity ◽  
Velocity Model ◽  
P Wave ◽  
Surface Velocity ◽  
S Wave ◽  
Converted Wave ◽  
Near Surface ◽  
Velocity Models ◽  
Consistent Solution ◽  
Static Corrections

A near-surface velocity model is one of the typical products generated when computing static corrections, particularly in the processing of PP data. Critically refracted waves are the input usually needed for this process. In addition, for the converted PS mode, S-wave near-surface corrections must be applied at the receiver locations. In this case, however, critically refracted S-waves are difficult to identify when using P-wave energy sources. We use the [Formula: see text]-[Formula: see text] representation of the converted-wave data to capture the intercept-time differences between receiver locations. These [Formula: see text]-differences are then used in the inversion of a near-surface S-wave velocity model. Our processing workflow provides not only a set of raypath-dependent S-wave static corrections, but also a velocity model that is based on those corrections. Our computed near-surface S-wave velocity model can be used for building migration velocity models or to initialize elastic full-waveform inversions. Our tests on synthetic and field data provided superior results to those obtained by using a surface-consistent solution.

scholarly journals Time-average velocity estimation through surface-wave analysis: Part 2 — P-wave velocity

Geophysics ◽  
2017 ◽  
Vol 82 (3) ◽  
pp. U61-U73 ◽  
Author(s):  
Laura Valentina Socco ◽  
Cesare Comina
Keyword(s):  
Wave Velocity ◽  
Poisson’S Ratio ◽  
Average Velocity ◽  
P Wave ◽  
Poisson's Ratio ◽  
S Wave ◽  
Seismic Line ◽  
Velocity Models ◽  
P Wave Velocity ◽  
Time Average

Surface waves (SWs) in seismic records can be used to extract local dispersion curves (DCs) along a seismic line. These curves can be used to estimate near-surface S-wave velocity models. If the velocity models are used to compute S-wave static corrections, the required information consists of S-wave time-average velocities that define the one-way time for a given datum plan depth. However, given the wider use of P-wave reflection seismic with respect to S-wave surveys, the estimate of P-wave time-average velocity would be more useful. We therefore focus on the possibility of also extracting time-average P-wave velocity models from SW dispersion data. We start from a known 1D S-wave velocity model along the line, with its relevant DC, and we estimate a wavelength/depth relationship for SWs. We found that this relationship is sensitive to Poisson’s ratio, and we develop a simple method for estimating an “apparent” Poisson’s ratio profile, defined as the Poisson’s ratio value that relates the time-average S-wave velocity to the time-average P-wave velocity. Hence, we transform the time-average S-wave velocity models estimated from the DCs into the time-average P-wave velocity models along the seismic line. We tested the method on synthetic and field data and found that it is possible to retrieve time-average P-wave velocity models with uncertainties mostly less than 10% in laterally varying sites and one-way traveltime for P-waves with less than 5 ms uncertainty with respect to P-wave tomography data. To our knowledge, this is the first method for reliable estimation of P-wave velocity from SW data without any a priori information or additional data.

scholarly journals The potential of seismic cross-hole tomography for geotechnical site investigation

2019 ◽  
Vol 92 ◽  
pp. 18006
Author(s):  
Yannick Choy Hing Ng ◽  
William Danovan ◽  
Taeseo Ku
Keyword(s):  
Wave Velocity ◽  
Site Investigation ◽  
Velocity Model ◽  
P Wave ◽  
Reconstruction Technique ◽  
Sufficient Information ◽  
Eastern Region ◽  
S Wave ◽  
Near Surface ◽  
Soil Layers

Seismic cross-hole tomography has been commonly used in oil and gas exploration and the mining industry for the detection of precious resources. For near-surface geotechnical site investigation, this geophysical method is relatively new and can be used to supplement traditional methods such as the standard penetration test, coring and sampling, thus improving the effectiveness of site characterization. This paper presents a case study which was carried out on a reclaimed land in the Eastern region of Singapore. A seismic cross-hole test was performed by generating both compressional waves and shear waves into the ground. The signals were interpreted by using first-arrival travel time wave tomography and the arrival times were subsequently inverted using Simultaneous Iterative Reconstruction Technique (SIRT). A comparison with the borehole logging data indicated that P-wave velocity model cannot provide sufficient information about the soil layers, especially when the ground water table is near the surface. The S-wave velocity model seemed to agree quite well with the variation in the SPT-N value and could identify to a certain extent the interface between the different soil layers. Finally, P-wave and S-wave velocities are used to compute the Poisson's ratio distribution which gave a good indication of the degree of saturation of the soil.

Full-waveform matching of VP and VS models from surface waves

2019 ◽  
Vol 218 (3) ◽  
pp. 1873-1891 ◽  
Author(s):  
Farbod Khosro Anjom ◽  
Daniela Teodor ◽  
Cesare Comina ◽  
Romain Brossier ◽  
Jean Virieux ◽  
...  
Keyword(s):  
Surface Waves ◽  
Wave Velocity ◽  
Real Data ◽  
Test Site ◽  
P Wave ◽  
S Wave ◽  
Surface Wave Dispersion ◽  
Data Set ◽  
Near Surface ◽  
Velocity Models

SUMMARY The analysis of surface wave dispersion curves (DCs) is widely used for near-surface S-wave velocity (VS) reconstruction. However, a comprehensive characterization of the near-surface requires also the estimation of P-wave velocity (VP). We focus on the estimation of both VS and VP models from surface waves using a direct data transform approach. We estimate a relationship between the wavelength of the fundamental mode of surface waves and the investigation depth and we use it to directly transform the DCs into VS and VP models in laterally varying sites. We apply the workflow to a real data set acquired on a known test site. The accuracy of such reconstruction is validated by a waveform comparison between field data and synthetic data obtained by performing elastic numerical simulations on the estimated VP and VS models. The uncertainties on the estimated velocity models are also computed.

PS Energy Imaging Condition for Microseismic Data - Part 1: Theory and Applications in 3D Isotropic Media

Geophysics ◽  
2020 ◽  
pp. 1-79
Author(s):  
Can Oren ◽  
Jeffrey Shragge
Keyword(s):  
Wave Velocity ◽  
Model Building ◽  
Velocity Model ◽  
Microseismic Monitoring ◽  
Elastic Model ◽  
S Wave ◽  
Imaging Condition ◽  
Velocity Models ◽  
Monitoring Programs ◽  
Imaging Conditions

Accurately estimating event locations is of significant importance in microseismic investigations because this information greatly contributes to the overall success of hydraulic fracturing monitoring programs. Full-wavefield time-reverse imaging (TRI) using one or more wave-equation imaging conditions offers an effective methodology for locating surface-recorded microseismic events. To be most beneficial in microseismic monitoring programs, though, the TRI procedure requires using accurate subsurface models that account for elastic media effects. We develop a novel microseismic (extended) PS energy imaging condition that explicitly incorporates the stiffness tensor and exhibits heightened sensitivity to isotropic elastic model perturbations compared to existing imaging conditions. Numerical experiments demonstrate the sensitivity of microseismic TRI results to perturbations in P- and S-wave velocity models. Zero-lag and extended microseismic source images computed at selected subsurface locations yields useful information about 3D P- and S-wave velocity model accuracy. Thus, we assert that these image volumes potentially can serve as the input into microseismic elastic velocity model building algorithms.

Crosshole reflection imaging with ground-penetrating radar data: Applications in near-surface sedimentary settings

Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. H61-H69
Author(s):  
Niklas Allroggen ◽  
Stéphane Garambois ◽  
Guy Sénéchal ◽  
Dominique Rousset ◽  
Jens Tronicke
Keyword(s):  
Ground Penetrating Radar ◽  
Field Data ◽  
Electromagnetic Property ◽  
Velocity Model ◽  
Radar Data ◽  
Near Surface ◽  
Velocity Models ◽  
Reflection Image ◽  
Detailed Interpretation ◽  
Ground Penetrating

Crosshole ground-penetrating radar (GPR) is applied in areas that require a very detailed subsurface characterization. Analysis of such data typically relies on tomographic inversion approaches providing an image of subsurface parameters. We have developed an approach for processing the reflected energy in crosshole GPR data and applied it on GPR data acquired in different sedimentary settings. Our approach includes muting of the first arrivals, separating the up- and the downgoing wavefield components, and backpropagating the reflected energy by a generalized Kirchhoff migration scheme. We obtain a reflection image that contains information on the location of electromagnetic property contrasts, thus outlining subsurface architecture in the interborehole plane. In combination with velocity models derived from different tomographic approaches, these images allow for a more detailed interpretation of subsurface structures without the need to acquire additional field data. In particular, a combined interpretation of the reflection image and the tomographic velocity model improves the ability to locate layer boundaries and to distinguish different subsurface units. To support our interpretations of our field data examples, we compare our crosshole reflection results with independent information, including borehole logs and surface GPR data.

Joint wave-equation traveltime inversion of diving/direct and reflected waves for the P- and S-wave velocity macromodel building

Geophysics ◽  
2021 ◽  
pp. 1-145
Author(s):  
Zhiming Ren ◽  
Qianzong Bao ◽  
Bingluo Gu
Keyword(s):  
Wave Velocity ◽  
Velocity Model ◽  
S Wave ◽  
Direct Wave ◽  
Reflected Waves ◽  
Traveltime Inversion ◽  
Long Wavelength ◽  
Velocity Models ◽  
Reflected Wave ◽  
Background Models

Full waveform inversion (FWI) suffers from the local minima problem and requires a sufficiently accurate starting model to converge to the correct solution. Wave-equation traveltime inversion (WETI) is an effective tool to retrieve the long-wavelength components of the velocity model. We develop a joint diving/direct and reflected wave WETI (JDRWETI) method to build the P- and S-wave velocity macromodels. We estimate the traveltime shifts of seismic events (diving/direct waves, PP and PS reflections) through the dynamic warping scheme and construct a misfit function using both the time shifts of diving/direct and reflected waves. We derive the adjoint wave equations and the gradients with respect to the background models based on the joint misfit function. We apply the kernel decomposition scheme to extract the kernel of the diving/direct wave and the tomography kernels of PP and PS reflections. For an explosive source, the kernels of diving/direct wave and PP reflections and the kernel of PS reflections are used to compute the P- and S-wave gradients of the background models, respectively. We implement JDRWETI by a two-stage inversion workflow: first invert the P- and S-wave velocity models using the P-wave gradients and then improve the S-wave velocity model using the S-wave gradients. Numerical tests on synthetic and field datasets reveal that the JDRWETI method successfully recovers the long-wavelength components of P- and S-wave velocity models, which can be used for an initial model for the subsequent elastic FWI. Moreover, the proposed JDRWETI method prevails over the existing reflection WETI method and the cascaded diving/direct and reflected wave WETI method, especially when large velocity errors are present in the shallow part of the starting models. The JDRWETI method with the two-stage inversion workflow can give rise to reasonable inversion results even for the model with different P- and S-wave velocity structures.

Elastic reflection waveform inversion based on the decomposition of sensitivity kernels

Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. R235-R250 ◽  
Author(s):  
Zhiming Ren ◽  
Zhenchun Li ◽  
Bingluo Gu
Keyword(s):  
Wave Velocity ◽  
Waveform Inversion ◽  
Velocity Model ◽  
P Wave ◽  
Background Model ◽  
Reverse Time ◽  
S Wave ◽  
Velocity Models ◽  
Time Migration ◽  
Starting Model

Full-waveform inversion (FWI) has the potential to obtain an accurate velocity model. Nevertheless, it depends strongly on the low-frequency data and the initial model. When the starting model is far from the real model, FWI tends to converge to a local minimum. Based on a scale separation of the model (into the background model and reflectivity model), reflection waveform inversion (RWI) can separate out the tomography term in the conventional FWI kernel and invert for the long-wavelength components of the velocity model by smearing the reflected wave residuals along the transmission (or “rabbit-ear”) paths. We have developed a new elastic RWI method to build the P- and S-wave velocity macromodels. Our method exploits a traveltime-based misfit function to highlight the contribution of tomography terms in the sensitivity kernels and a sensitivity kernel decomposition scheme based on the P- and S-wave separation to suppress the high-wavenumber artifacts caused by the crosstalk of different wave modes. Numerical examples reveal that the gradients of the background models become sufficiently smooth owing to the decomposition of sensitivity kernels and the traveltime-based misfit function. We implement our elastic RWI in an alternating way. At each loop, the reflectivity model is generated by elastic least-squares reverse time migration, and then the background model is updated using the separated traveltime kernels. Our RWI method has been successfully applied in synthetic and real reflection seismic data. Inversion results demonstrate that the proposed method can retrieve preferable low-wavenumber components of the P- and S-wave velocity models, which are reliable to serve as a starting model for conventional elastic FWI. Also, our method with a two-stage inversion workflow, first updating the P-wave velocity using the PP kernels and then updating the S-wave velocity using the PS kernels, is feasible and robust even when P- and S-wave velocities have different structures.

Sign in / Sign up

Export Citation Format

Share Document