Tunnel detection at Yuma Proving Ground, Arizona, USA — Part 2: 3D full-waveform inversion experiments

Geophysics ◽  
2019 ◽  
Vol 84 (1) ◽  
pp. B107-B120 ◽  
Author(s):  
James A. Smith ◽  
Dmitry Borisov ◽  
Harley Cudney ◽  
Richard D. Miller ◽  
Ryan Modrak ◽  
...  
Keyword(s):  
Waveform Inversion ◽  
Full Waveform Inversion ◽  
3D Image ◽  
Initial Model ◽  
3D Inversion ◽  
Data Set ◽  
Tunnel Detection ◽  
Full Waveform ◽  
Wave Methods ◽  
Reduced Data

We have applied time-domain 3D elastic full-waveform inversion (FWI) to a known tunnel constructed 10 m below the surface with no distinguishing surface expressions. Multicomponent inversion experiments that use an initial model estimated from surface wave methods suggest that the vertical sources and the combination of vertical and longitudinal receivers result in the clearest image of the tunnel. We obtain an approximate 3D image of the tunnel using 24 vertical sources and 720 vertical and 720 longitudinal receivers. We find that increasing the number of vertical sources to 216 does not significantly improve the details of the tunnel. Further experimentation indicates that we can detect the tunnel using a reduced data set of 6–10 vertical sources and 216 vertical sensors. In addition, calculating the [Formula: see text] ratio helps remove artifacts from the inverted model and highlights the location of the tunnel. We compare the 3D inversion results with the 2D FWI results for the same tunnel that we previously evaluated. The variety of successful inversion experiments suggest that FWI is capable of imaging shallow tunnels in desert geologies.

Tunnel detection at Yuma Proving Ground, Arizona, USA — Part 1: 2D full-waveform inversion experiment

Geophysics ◽  
2019 ◽  
Vol 84 (1) ◽  
pp. B95-B105 ◽  
Author(s):  
Yao Wang ◽  
Richard D. Miller ◽  
Shelby L. Peterie ◽  
Steven D. Sloan ◽  
Mark L. Moran ◽  
...  
Keyword(s):  
Surface Wave ◽  
Waveform Inversion ◽  
Velocity Profiles ◽  
Infinite Length ◽  
Full Waveform Inversion ◽  
S Wave ◽  
Near Surface ◽  
Tunnel Detection ◽  
Full Waveform ◽  
Wave Methods

We have applied time domain 2D full-waveform inversion (FWI) to detect a known 10 m deep wood-framed tunnel at Yuma Proving Ground, Arizona. The acquired seismic data consist of a series of 2D survey lines that are perpendicular to the long axis of the tunnel. With the use of an initial model estimated from surface wave methods, a void-detection-oriented FWI workflow was applied. A straightforward [Formula: see text] quotient masking method was used to reduce the inversion artifacts and improve confidence in identifying anomalies that possess a high [Formula: see text] ratio. Using near-surface FWI, [Formula: see text] and [Formula: see text] velocity profiles were obtained with void anomalies that are easily interpreted. The inverted velocity profiles depict the tunnel as a low-velocity anomaly at the correct location and depth. A comparison of the observed and simulated waveforms demonstrates the reliability of inverted models. Because the known tunnel has a uniform shape and for our purposes an infinite length, we apply 1D interpolation to the inverted [Formula: see text] profiles to generate a pseudo 3D (2.5D) volume. Based on this research, we conclude the following: (1) FWI is effective in near-surface tunnel detection when high resolution is necessary. (2) Surface-wave methods can provide accurate initial S-wave velocity [Formula: see text] models for near-surface 2D FWI.

A new parameter set for anisotropic multiparameter full-waveform inversion and application to a North Sea data set

Geophysics ◽  
2016 ◽  
Vol 81 (4) ◽  
pp. U25-U38 ◽  
Author(s):  
Nuno V. da Silva ◽  
Andrew Ratcliffe ◽  
Vetle Vinje ◽  
Graham Conroy
Keyword(s):  
North Sea ◽  
Waveform Inversion ◽  
Real Data ◽  
Model Parameters ◽  
Full Waveform Inversion ◽  
Radiation Patterns ◽  
Data Set ◽  
Full Waveform ◽  
Seismic Image ◽  
Central North Sea

Parameterization lies at the center of anisotropic full-waveform inversion (FWI) with multiparameter updates. This is because FWI aims to update the long and short wavelengths of the perturbations. Thus, it is important that the parameterization accommodates this. Recently, there has been an intensive effort to determine the optimal parameterization, centering the fundamental discussion mainly on the analysis of radiation patterns for each one of these parameterizations, and aiming to determine which is best suited for multiparameter inversion. We have developed a new parameterization in the scope of FWI, based on the concept of kinematically equivalent media, as originally proposed in other areas of seismic data analysis. Our analysis is also based on radiation patterns, as well as the relation between the perturbation of this set of parameters and perturbation in traveltime. The radiation pattern reveals that this parameterization combines some of the characteristics of parameterizations with one velocity and two Thomsen’s parameters and parameterizations using two velocities and one Thomsen’s parameter. The study of perturbation of traveltime with perturbation of model parameters shows that the new parameterization is less ambiguous when relating these quantities in comparison with other more commonly used parameterizations. We have concluded that our new parameterization is well-suited for inverting diving waves, which are of paramount importance to carry out practical FWI successfully. We have demonstrated that the new parameterization produces good inversion results with synthetic and real data examples. In the latter case of the real data example from the Central North Sea, the inverted models show good agreement with the geologic structures, leading to an improvement of the seismic image and flatness of the common image gathers.

scholarly journals Normalized nonzero-lag crosscorrelation elastic full-waveform inversion

Geophysics ◽  
2019 ◽  
Vol 84 (1) ◽  
pp. R1-R10 ◽  
Author(s):  
Zhendong Zhang ◽  
Tariq Alkhalifah ◽  
Zedong Wu ◽  
Yike Liu ◽  
Bin He ◽  
...  
Keyword(s):  
Field Data ◽  
Extreme Values ◽  
Time Lag ◽  
Waveform Inversion ◽  
Optimal Time ◽  
Full Waveform Inversion ◽  
Data Set ◽  
Velocity Models ◽  
Full Waveform ◽  
The Earth

Full-waveform inversion (FWI) is an attractive technique due to its ability to build high-resolution velocity models. Conventional amplitude-matching FWI approaches remain challenging because the simplified computational physics used does not fully represent all wave phenomena in the earth. Because the earth is attenuating, a sample-by-sample fitting of the amplitude may not be feasible in practice. We have developed a normalized nonzero-lag crosscorrelataion-based elastic FWI algorithm to maximize the similarity of the calculated and observed data. We use the first-order elastic-wave equation to simulate the propagation of seismic waves in the earth. Our proposed objective function emphasizes the matching of the phases of the events in the calculated and observed data, and thus, it is more immune to inaccuracies in the initial model and the difference between the true and modeled physics. The normalization term can compensate the energy loss in the far offsets because of geometric spreading and avoid a bias in estimation toward extreme values in the observed data. We develop a polynomial-type weighting function and evaluate an approach to determine the optimal time lag. We use a synthetic elastic Marmousi model and the BigSky field data set to verify the effectiveness of the proposed method. To suppress the short-wavelength artifacts in the estimated S-wave velocity and noise in the field data, we apply a Laplacian regularization and a total variation constraint on the synthetic and field data examples, respectively.

Tackling cycle skipping in full-waveform inversion with intermediate data

Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. R411-R427 ◽  
Author(s):  
Gang Yao ◽  
Nuno V. da Silva ◽  
Michael Warner ◽  
Di Wu ◽  
Chenhao Yang
Keyword(s):  
Objective Function ◽  
Waveform Inversion ◽  
Synthetic Data ◽  
Velocity Model ◽  
Full Waveform Inversion ◽  
Data Set ◽  
Intermediate Data ◽  
Full Waveform ◽  
L2 Norm ◽  
Recorded Data

Full-waveform inversion (FWI) is a promising technique for recovering the earth models for exploration geophysics and global seismology. FWI is generally formulated as the minimization of an objective function, defined as the L2-norm of the data residuals. The nonconvex nature of this objective function is one of the main obstacles for the successful application of FWI. A key manifestation of this nonconvexity is cycle skipping, which happens if the predicted data are more than half a cycle away from the recorded data. We have developed the concept of intermediate data for tackling cycle skipping. This intermediate data set is created to sit between predicted and recorded data, and it is less than half a cycle away from the predicted data. Inverting the intermediate data rather than the cycle-skipped recorded data can then circumvent cycle skipping. We applied this concept to invert cycle-skipped first arrivals. First, we picked up the first breaks of the predicted data and the recorded data. Second, we linearly scaled down the time difference between the two first breaks of each shot into a series of time shifts, the maximum of which was less than half a cycle, for each trace in this shot. Third, we moved the predicted data with the corresponding time shifts to create the intermediate data. Finally, we inverted the intermediate data rather than the recorded data. Because the intermediate data are not cycle-skipped and contain the traveltime information of the recorded data, FWI with intermediate data updates the background velocity model in the correct direction. Thus, it produces a background velocity model accurate enough for carrying out conventional FWI to rebuild the intermediate- and short-wavelength components of the velocity model. Our numerical examples using synthetic data validate the intermediate-data concept for tackling cycle skipping and demonstrate its effectiveness for the application to first arrivals.

scholarly journals Semiglobal viscoacoustic full-waveform inversion

Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. R271-R293 ◽  
Author(s):  
Nuno V. da Silva ◽  
Gang Yao ◽  
Michael Warner
Keyword(s):  
Physical Properties ◽  
Wave Velocity ◽  
Seismic Data ◽  
Waveform Inversion ◽  
P Wave ◽  
Full Waveform Inversion ◽  
Data Set ◽  
Intrinsic Attenuation ◽  
Full Waveform ◽  
P Wave Velocity

Full-waveform inversion deals with estimating physical properties of the earth’s subsurface by matching simulated to recorded seismic data. Intrinsic attenuation in the medium leads to the dispersion of propagating waves and the absorption of energy — media with this type of rheology are not perfectly elastic. Accounting for that effect is necessary to simulate wave propagation in realistic geologic media, leading to the need to estimate intrinsic attenuation from the seismic data. That increases the complexity of the constitutive laws leading to additional issues related to the ill-posed nature of the inverse problem. In particular, the joint estimation of several physical properties increases the null space of the parameter space, leading to a larger domain of ambiguity and increasing the number of different models that can equally well explain the data. We have evaluated a method for the joint inversion of velocity and intrinsic attenuation using semiglobal inversion; this combines quantum particle-swarm optimization for the estimation of the intrinsic attenuation with nested gradient-descent iterations for the estimation of the P-wave velocity. This approach takes advantage of the fact that some physical properties, and in particular the intrinsic attenuation, can be represented using a reduced basis, substantially decreasing the dimension of the search space. We determine the feasibility of the method and its robustness to ambiguity with 2D synthetic examples. The 3D inversion of a field data set for a geologic medium with transversely isotropic anisotropy in velocity indicates the feasibility of the method for inverting large-scale real seismic data and improving the data fitting. The principal benefits of the semiglobal multiparameter inversion are the recovery of the intrinsic attenuation from the data and the recovery of the true undispersed infinite-frequency P-wave velocity, while mitigating ambiguity between the estimated parameters.

Application of multiparameter full-waveform inversion in Central Luconia Basin, Sarawak

2016 ◽  
Vol 4 (4) ◽  
pp. SU17-SU24 ◽  
Author(s):  
Vanessa Goh ◽  
Kjetil Halleland ◽  
René-Édouard Plessix ◽  
Alexandre Stopin
Keyword(s):  
Model Building ◽  
Waveform Inversion ◽  
Velocity Model ◽  
Thin Layers ◽  
Earth Model ◽  
Full Waveform Inversion ◽  
Data Set ◽  
Full Waveform ◽  
Velocity Variations ◽  
Central Luconia

Reducing velocity inaccuracy in complex settings is of paramount importance for limiting structural uncertainties, therefore helping the geologic interpretation and reservoir characterization. Shallow velocity variations due, for instance, to gas accumulations or carbonate reefs, are a common issue offshore Malaysia. These velocity variations are difficult to image through standard reflection-based velocity model building. We have applied full-waveform inversion (FWI) to better characterize the upper part of the earth model for a shallow-water field, located in the Central Luconia Basin offshore Sarawak. We have inverted a narrow-azimuth data set with a maximum inline offset of 4.4 km. Thanks to dedicated broadband preprocessing of the data set, we could enhance the signal-to-noise ratio in the 2.5–10 Hz frequency band. We then applied a multiparameter FWI to estimate the background normal moveout velocity and the [Formula: see text]-parameter. Full-waveform inversion together with broadband data processing has helped to better define the faults and resolve the thin layers in the shallow clastic section. The improvements in the velocity model brought by FWI lead to an improved image of the structural closure and flanks. Moreover, the increased velocity resolution helps in distinguishing between two different geologic interpretations.

Sign in / Sign up

Export Citation Format

Share Document