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.

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.

Multiwave inversion: A key step for depth model building — Examples from the Sultanate of Oman

2021 ◽  
Vol 40 (8) ◽  
pp. 610-618
Author(s):  
Daniela Donno ◽  
Mohammad Sheryar Farooqui ◽  
Mostafa Khalil ◽  
David McCarthy ◽  
Didier Solyga ◽  
...  
Keyword(s):  
Surface Waves ◽  
Model Building ◽  
Sand Dunes ◽  
P Wave ◽  
Surface Model ◽  
Sultanate Of Oman ◽  
Surface Wave Dispersion ◽  
Low Frequencies ◽  
Near Surface ◽  
Velocity Models

The near surface in the Middle East, particularly in the Sultanate of Oman, is characterized by very shallow high-velocity carbonates and anhydrites interleaved by slow-velocity clastic layers, resulting in sharp velocity inversions in the first few hundred meters below the surface. In addition, the surface is characterized by features such as unconsolidated materials within dry riverbeds (known as “wadis”), small jebels, and sand dunes, which cause distortions in the underlying shallow and deeper seismic images. This work presents the building of a near-surface model by using multiwave inversion that jointly inverts information from P-wave first breaks and surface-wave dispersion curves. The use of surface waves in combination with first breaks captures the lateral and vertical velocity variations, especially in the shallowest parts of the near surface. This paper focuses on the analysis of two drawbacks of this technology: the picking of the input data information, which can be cumbersome and time consuming, and the limited penetration depth of surface waves at the typical frequencies of active data. To overcome these issues, an innovative workflow is proposed that combines the use of an unsupervised machine learning technique to guide the pick extraction phase and the reconstruction of ultra-low-frequency surface waves (0.5 to 1.5 Hz) through an interferometry process using information from natural and ambient noise. Deeper near-surface P- and S-wave velocity models can be obtained with multiwave inversion using these ultra-low frequencies. The integration of a near-surface model into the velocity model building workflow brings a major improvement in depth imaging from shallow to deep structures, as demonstrated on two data sets from the Sultanate of Oman.

Tomographic modeling of a cross‐borehole data set

Geophysics ◽  
1989 ◽  
Vol 54 (10) ◽  
pp. 1249-1257 ◽  
Author(s):  
Larry R. Lines ◽  
Edward D. LaFehr
Keyword(s):  
Wave Equation ◽  
Ray Tracing ◽  
Velocity Model ◽  
Velocity Estimation ◽  
P Wave ◽  
Seismic Survey ◽  
Equation Modeling ◽  
Data Set ◽  
True Solution ◽  
Velocity Models

In this paper we describe a methodology for estimating P‐wave velocities from a cross‐borehole seismic survey that uses straight‐ray tomography, ray tracing, and finite‐difference wave‐equation modeling to produce velocity models that fit the first‐break traveltimes. After a starting model is established by straight‐ray tomography, the velocity model is checked by ray tracing and wave‐equation modeling. Since the models for each procedure show consistent results and the modeled traveltimes closely match those traveltimes from the actual data, we felt our interpretation was confirmed. However, the fitting of cross‐well first break traveltimes is only a necessary validity check and is not sufficient to guarantee that the true solution has been found. Two wells were drilled through the areas that were anomalous on the derived tomogram and check‐shot velocity surveys were run. Due primarily to a lateral ambiguity in velocity estimation caused by too few near‐vertical raypaths, the check‐shot surveys did not agree with the tomogram velocities. However, subsequently the check‐shot traveltimes were used to place bounds on velocity in a constrained least‐squares procedure; the combined modeling of uphole and cross‐well rays produced an optimum velocity model which satisfies all available data.

scholarly journals Mitigating elastic effects in marine 3-D full-waveform inversion

2019 ◽  
Vol 220 (3) ◽  
pp. 2089-2104
Author(s):  
Òscar Calderón Agudo ◽  
Nuno Vieira da Silva ◽  
George Stronge ◽  
Michael Warner
Keyword(s):  
Wave Equation ◽  
Field Data ◽  
Waveform Inversion ◽  
Velocity Model ◽  
P Wave ◽  
Data Sets ◽  
Acoustic Wave Equation ◽  
Data Set ◽  
Velocity Models ◽  
Full Waveform

SUMMARY The potential of full-waveform inversion (FWI) to recover high-resolution velocity models of the subsurface has been demonstrated in the last decades with its application to field data. But in certain geological scenarios, conventional FWI using the acoustic wave equation fails in recovering accurate models due to the presence of strong elastic effects, as the acoustic wave equation only accounts for compressional waves. This becomes more critical when dealing with land data sets, in which elastic effects are generated at the source and recorded directly by the receivers. In marine settings, in which sources and receivers are typically within the water layer, elastic effects are weaker but can be observed most easily as double mode conversions and through their effect on P-wave amplitudes. Ignoring these elastic effects can have a detrimental impact on the accuracy of the recovered velocity models, even in marine data sets. Ideally, the elastic wave equation should be used to model wave propagation, and FWI should aim to recover anisotropic models of velocity for P waves (vp) and S waves (vs). However, routine three-dimensional elastic FWI is still commercially impractical due to the elevated computational cost of modelling elastic wave propagation in regions with low S-wave velocity near the seabed. Moreover, elastic FWI using local optimization methods suffers from cross-talk between different inverted parameters. This generally leads to incorrect estimation of subsurface models, requiring an estimate of vp/vs that is rarely known beforehand. Here we illustrate how neglecting elasticity during FWI for a marine field data set that contains especially strong elastic heterogeneities can lead to an incorrect estimation of the P-wave velocity model. We then demonstrate a practical approach to mitigate elastic effects in 3-D yielding improved estimates, consisting of using a global inversion algorithm to estimate a model of vp/vs, employing matching filters to remove elastic effects from the field data, and performing acoustic FWI of the resulting data set. The quality of the recovered models is assessed by exploring the continuity of the events in the migrated sections and the fit of the latter with the recovered velocity model.

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.

A method for correcting acoustic finite-difference amplitudes for elastic effects

Geophysics ◽  
2014 ◽  
Vol 79 (4) ◽  
pp. T243-T255 ◽  
Author(s):  
James W. D. Hobro ◽  
Chris H. Chapman ◽  
Johan O. A. Robertsson
Keyword(s):  
Wave Equation ◽  
Finite Difference ◽  
Effective Means ◽  
Cost Effective ◽  
P Wave ◽  
S Waves ◽  
Velocity Models ◽  
Acoustic Simulation ◽  
Wide Range ◽  
Elastic Simulation

We present a new method for correcting the amplitudes of arrivals in an acoustic finite-difference simulation for elastic effects. In this method, we selectively compute an estimate of the error incurred when the acoustic wave equation is used to approximate the behavior of the elastic wave equation. This error estimate is used to generate an effective source field in a second acoustic simulation. The result of this second simulation is then applied as a correction to the original acoustic simulation. The overall cost is approximately twice that of an acoustic simulation but substantially less than the cost of an elastic simulation. Because both simulations are acoustic, no S-waves are generated, so dispersed converted waves are avoided. We tested the characteristics of the method on a simple synthetic model designed to simulate propagation through a strong acoustic impedance contrast representative of sedimentary geology. It corrected amplitudes to high accuracy for reflected arrivals over a wide range of incidence angles. We also evaluated results from simulations on more complex models that demonstrated that the method was applicable in realistic sedimentary models containing a wide range of seismic contrasts. However, its accuracy was reduced for wide-angle reflections from very high impedance contrasts such as a shallow top-salt interface. We examined the influence of modeling at coarse grid resolutions, in which converted S-waves in the equivalent elastic simulation are dispersed. These results provide some validation for the accuracy of the method when applied using finite-difference grids designed for acoustic modeling. The method appears to offer a cost-effective means of modeling elastic amplitudes for P-wave arrivals in a useful range of velocity models. It has several potential applications in imaging and inversion.

Shallow velocity structure and poisson's ratio at the Tarzana, California, strong-motion accelerometer site

1996 ◽  
Vol 86 (6) ◽  
pp. 1704-1713 ◽  
Author(s):  
R. D. Catchings ◽  
W. H. K. Lee
Keyword(s):  
Urban Areas ◽  
Velocity Structure ◽  
Strong Motion ◽  
P Wave ◽  
S Wave ◽  
Near Surface ◽  
Velocity Models ◽  
Wave Velocities ◽  
Poisson's Ratios ◽  
Poisson’S Ratios

Abstract The 17 January 1994, Northridge, California, earthquake produced strong ground shaking at the Cedar Hills Nursery (referred to here as the Tarzana site) within the city of Tarzana, California, approximately 6 km from the epicenter of the mainshock. Although the Tarzana site is on a hill and is a rock site, accelerations of approximately 1.78 g horizontally and 1.2 g vertically at the Tarzana site are among the highest ever instrumentally recorded for an earthquake. To investigate possible site effects at the Tarzana site, we used explosive-source seismic refraction data to determine the shallow (<70 m) P-and S-wave velocity structure. Our seismic velocity models for the Tarzana site indicate that the local velocity structure may have contributed significantly to the observed shaking. P-wave velocities range from 0.9 to 1.65 km/sec, and S-wave velocities range from 0.20 and 0.6 km/sec for the upper 70 m. We also found evidence for a local S-wave low-velocity zone (LVZ) beneath the top of the hill. The LVZ underlies a CDMG strong-motion recording site at depths between 25 and 60 m below ground surface (BGS). Our velocity model is consistent with the near-surface (<30 m) P- and S-wave velocities and Poisson's ratios measured in a nearby (<30 m) borehole. High Poisson's ratios (0.477 to 0.494) and S-wave attenuation within the LVZ suggest that the LVZ may be composed of highly saturated shales of the Modelo Formation. Because the lateral dimensions of the LVZ approximately correspond to the areas of strongest shaking, we suggest that the highly saturated zone may have contributed to localized strong shaking. Rock sites are generally considered to be ideal locations for site response in urban areas; however, localized, highly saturated rock sites may be a hazard in urban areas that requires further investigation.

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.

Random objective waveform inversion of surface waves

Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. EN49-EN61
Author(s):  
Yudi Pan ◽  
Lingli Gao
Keyword(s):  
Objective Function ◽  
Surface Waves ◽  
Waveform Inversion ◽  
Velocity Model ◽  
Initial Model ◽  
S Wave ◽  
Data Set ◽  
Near Surface ◽  
Synthetic Test ◽  
Ground Penetrating

Full-waveform inversion (FWI) of surface waves is becoming increasingly popular among shallow-seismic methods. Due to a huge amount of data and the high nonlinearity of the objective function, FWI usually requires heavy computational costs and may converge toward a local minimum. To mitigate these problems, we have reformulated FWI under a multiobjective framework and adopted a random objective waveform inversion (ROWI) method for surface-wave characterization. Three different measure functions were used, whereas the combination of one measure function with one shot independently provided one of the [Formula: see text] objective functions ([Formula: see text] is the total number of shots). We have randomly chose and optimized one objective function at each iteration. We performed a synthetic test to compare the performance of the ROWI and conventional FWI approaches, which showed that the convergence of ROWI is faster and more robust compared with conventional FWI approaches. We also applied ROWI to a field data set acquired in Rheinstetten, Germany. ROWI successfully reconstructed the main geologic feature, a refilled trench, in the final result. The comparison between the ROWI result and a migrated ground-penetrating radar profile further proved the effectiveness of ROWI in reconstructing the near-surface S-wave velocity model. We also ran the same field example by using a poor initial model. In this case, conventional FWI failed whereas ROWI still reconstructed the subsurface model to a fairly good level, which highlighted the relatively low dependency of ROWI on the initial model.

Inversion of seismic refraction and reflection data for building long-wavelength velocity models

Geophysics ◽  
2015 ◽  
Vol 80 (2) ◽  
pp. R81-R93 ◽  
Author(s):  
Haiyang Wang ◽  
Satish C. Singh ◽  
Francois Audebert ◽  
Henri Calandra
Keyword(s):  
Wave Equation ◽  
Short Wavelength ◽  
Velocity Model ◽  
Density Model ◽  
Data Set ◽  
Reflected Waves ◽  
Long Wavelength ◽  
Velocity Models ◽  
Computation Efficiency ◽  
Refracted Waves

Long-wavelength velocity model building is a nonlinear process. It has traditionally been achieved without appealing to wave-equation-based approaches for combined refracted and reflected waves. We developed a cascaded wave-equation tomography method in the data domain, taking advantage of the information contained in the reflected and refracted waves. The objective function was the traveltime residual that maximized the crosscorrelation function between real and synthetic data. To alleviate the nonlinearity of the inversion problem, refracted waves were initially used to provide vertical constraints on the velocity model, and reflected waves were then included to provide lateral constraints. The use of reflected waves required scale separation. We separated the long- and short-wavelength subsurface structures into velocity and density models, respectively. The velocity model update was restricted to long wavelengths during the wave-equation tomography, whereas the density model was used to absorb all the short-wavelength impedance contrasts. To improve the computation efficiency, the density model was converted into the zero-offset traveltime domain, where it was invariant to changes of the long-wavelength velocity model. After the wave-equation tomography has derived an optimized long-wavelength velocity model, full-waveform inversion was used to invert all the data to retrieve the short-wavelength velocity structures. We developed our method in two synthetic tests and then applied it to a marine field data set. We evaluated the results of the use of refracted and reflected waves, which was critical for accurately building the long-wavelength velocity model. We showed that our wave-equation tomography strategy was robust for the real data application.

Sign in / Sign up

Export Citation Format

Share Document