Sinkhole detection with 3D full seismic waveform tomography

Geophysics ◽  
2020 ◽  
Vol 85 (5) ◽  
pp. B169-B179
Author(s):  
Majid Mirzanejad ◽  
Khiem T. Tran ◽  
Michael McVay ◽  
David Horhota ◽  
Scott J. Wasman

Sinkhole collapse may result in significant property damage and even loss of life. Early detection of sinkhole attributes (buried voids, raveling zones) is critical to limit the cost of remediation. One of the most promising ways to obtain subsurface imaging is 3D seismic full-waveform inversion. For demonstration, a recently developed 3D Gauss-Newton full-waveform inversion (3D GN-FWI) method is used to detect buried voids, raveling soils, and characterize variable subsurface soil/rock layering. It is based on a finite-difference solution of 3D elastic wave equations and Gauss-Newton optimization. The method is tested first on a data set constructed from the numerical simulation of a challenging synthetic model and subsequently on field data collected from two separate test sites in Florida. For the field tests, receivers and sources are placed in uniform 2D surface grids to acquire the seismic wavefields, which then are inverted to extract the 3D subsurface velocity structures. The inverted synthetic results suggest that the approach is viable for detecting voids and characterizing layering. The field seismic results reveal that the 3D waveform analysis identified a known manmade void (plastic culvert), unknown natural voids, raveling, as well as laterally variable soil/rock layering including rock pinnacles. The results are confirmed later by standard penetration tests, including depth to bedrock, two buried voids, and a raveling soil zone. Our study provides insight into the application of the 3D seismic FWI technique as a powerful tool in detecting shallow voids and other localized subsurface features.

2021 ◽  
Vol 64 (1) ◽  
Author(s):  
Zhanyuan Liang ◽  
Guochen Wu ◽  
Xiaoyu Zhang ◽  
Qingyang Li

Reflection full-waveform inversion (RFWI) can reduce the nonlinearity of inversion providing an accurate initial velocity model for full-waveform inversion (FWI) through the tomographic components (low-wavenumber). However, elastic-wave reflection full-waveform inversion (ERFWI) is more vulnerable to the problem of local minimum due to the complicated multi-component wavefield. Our algorithm first divides kernels of ERFWI into four subkernels based on the theory of decoupled elastic-wave equations. Then we try to construct the tomographic components of ERFWI with only single-component wavefields, similarly to acoustic inversions. However, the S-wave velocity is still vulnerable to the coupling effects because of P-wave components contained in the S-wavefield in an inhomogeneous medium. Therefore we develop a method for decoupling elastic- wave equations in an inhomogeneous medium, which is applied to the decomposition of kernels in ERFWI. The new decoupled system can improve the accuracy of S-wavefield and hence further reduces the high-wavenumber crosstalk in the subkernel of S-wave velocity after kernels are decomposed. The numerical examples of Sigsbee2A model demonstrate that our ERFWI method with decoupled elastic-wave equations can efficiently recover the low-wavenumber components of the model and improve the precision of the S-wave velocity.


Geophysics ◽  
2016 ◽  
Vol 81 (4) ◽  
pp. U25-U38 ◽  
Author(s):  
Nuno V. da Silva ◽  
Andrew Ratcliffe ◽  
Vetle Vinje ◽  
Graham Conroy

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.


Geophysics ◽  
2019 ◽  
Vol 84 (1) ◽  
pp. R1-R10 ◽  
Author(s):  
Zhendong Zhang ◽  
Tariq Alkhalifah ◽  
Zedong Wu ◽  
Yike Liu ◽  
Bin He ◽  
...  

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.


2018 ◽  
Vol 37 (2) ◽  
pp. 142-145 ◽  
Author(s):  
Philipp Witte ◽  
Mathias Louboutin ◽  
Keegan Lensink ◽  
Michael Lange ◽  
Navjot Kukreja ◽  
...  

This tutorial is the third part of a full-waveform inversion (FWI) tutorial series with a step-by-step walkthrough of setting up forward and adjoint wave equations and building a basic FWI inversion framework. For discretizing and solving wave equations, we use Devito ( http://www.opesci.org/devito-public ), a Python-based domain-specific language for automated generation of finite-difference code ( Lange et al., 2016 ). The first two parts of this tutorial ( Louboutin et al., 2017 , 2018 ) demonstrated how to solve the acoustic wave equation for modeling seismic shot records and how to compute the gradient of the FWI objective function using the adjoint-state method. With these two key ingredients, we will now build an inversion framework that can be used to minimize the FWI least-squares objective function.


Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. R411-R427 ◽  
Author(s):  
Gang Yao ◽  
Nuno V. da Silva ◽  
Michael Warner ◽  
Di Wu ◽  
Chenhao Yang

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.


Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. R271-R293 ◽  
Author(s):  
Nuno V. da Silva ◽  
Gang Yao ◽  
Michael Warner

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.


2020 ◽  
Author(s):  
Anselme F.E. Borgeaud ◽  
Maria Koroni ◽  
Frédéric Deschamps

<div>This study presents a new approach for investigating the structure of the core-mantle boundary (CMB)</div><div>topography based on full-waveforms and adjoint methods. We compute intermediate period (10-20 seconds)</div><div>spectral-element seismograms using existing models of core-mantle boundary topography and we analyse the</div><div>sensitivity of relevant seismic phases. Our study adds new information about effects of CMB structure on</div><div>exact synthetics and observable traveltimes of seismic body waves by means of sensitivity kernels. It also</div><div>highlights the difficulty of imaging the boundary due to the strong trade-off between velocity and topography</div><div>variations, addressed by many previous investigators.</div><div> </div><div>Given the significance of CMB and its importance for many disciplines in geophysical research, there have</div><div>been many studies trying to understand and geographically map the variations of topography and velocity</div><div>above this boundary. However, the vast mantle area wherein seismic waves travel before and after they</div><div>interact with the CMB makes the identification of desired seismic phases somehow difficult. In addition, the</div><div>observable traveltimes can be hard to interpret as a result of the boundary’s topography only, due to the</div><div>approximate inverse methods and limited modelling methodologies. Despite considerable progress made the</div><div>past years, there is still a necessity for improving the understanding of effects of core-mantle boundary and</div><div>D″ structure on recorded waveforms.</div><div> </div><div>For our analyses, we perform comparisons between time shifts due to topography made on full-waveform</div><div>synthetics to ray theoretical predictions in order to assess methods usually deployed for imaging CMB.</div><div>Then, we calculate the corresponding sensitivity kernel for time windows around the theoretical arrival of</div><div>each phase. We focus on diffracted, core reflected and refracted <em>P</em> and <em>S</em> waves. The sensitivity kernels</div><div>depict the finite-frequency nature of these waves and possible contributions from other phases unpredictable</div><div>by ray theory. Results show that for most phases ray theory performs acceptably with some accuracy loss,</div><div>however comparisons of the effect of velocity variations to topography on traveltimes are discouraging due</div><div>to the low sensitivity to the latter.</div><div> </div><div>The conclusions drawn by our traveltime and sensitivity analyses are twofold. Firstly, using spectral-</div><div>element waveforms, the seismic phases which are frequently found in literature can be thoroughly investigated</div><div>and better understood, since their traveltime sensitivity through mantle and core is explicitly shown. The</div><div>full-waveform analysis allows us to assess the usability of phases which are informative for core-mantle</div><div>boundary structure and its topography. Secondly, we propose that using the analysed phases simultaneously</div><div>in a full-waveform inversion scheme will improve imaging of the CMB, while also allowing to jointly invert</div><div>for velocity variations along the D″ layer, which is generally poorly understood. From this study, we want</div><div>to promote advanced techniques of full-waveform inversion for improving CMB and lower mantle models.</div>


2016 ◽  
Vol 4 (4) ◽  
pp. SU17-SU24 ◽  
Author(s):  
Vanessa Goh ◽  
Kjetil Halleland ◽  
René-Édouard Plessix ◽  
Alexandre Stopin

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.


2018 ◽  
Vol 37 (1) ◽  
pp. 69-72 ◽  
Author(s):  
Mathias Louboutin ◽  
Philipp Witte ◽  
Michael Lange ◽  
Navjot Kukreja ◽  
Fabio Luporini ◽  
...  

This is the second part of a three-part tutorial series on full-waveform inversion (FWI) in which we provide a step-by-step walk through of setting up forward and adjoint wave equation solvers and an optimization framework for inversion. In Part 1 ( Louboutin et al., 2017 ), we showed how to use Devito ( http://www.opesci.org/devito-public ) to set up and solve acoustic wave equations with (impulsive) seismic sources and sample wavefields at the receiver locations to forward model shot records. Here in Part 2, we will discuss how to set up and solve adjoint wave equations with Devito and, from that, how we can calculate gradients and function values of the FWI objective function.


Sign in / Sign up

Export Citation Format

Share Document