Reparameterized full-waveform inversion using deep neural networks

Geophysics ◽  
2020 ◽  
Vol 86 (1) ◽  
pp. V1-V13 ◽  
Author(s):  
Qinglong He ◽  
Yanfei Wang

Full-waveform inversion (FWI) is a powerful method for providing a high-resolution description of the subsurface. However, the misfit function of the conventional FWI method (metric [Formula: see text]-norm) is usually dominated by spurious local minima owing to its nonlinearity and ill-posedness. In addition, FWI requires intensive wavefield computation to evaluate the gradient and step length. We have considered a general inversion method using a deep neural network (DNN) for the FWI problem. This deep-learning inversion method reparameterizes physical parameters using the weights of a DNN, such that the inversion amounts to reconstructing these weights. One advantage of this deep-learning inversion method is that it can serve as an iterative regularization method, benefiting from the representation of the network. Thus, it is suitable to solve ill-posed nonlinear inverse problems. Furthermore, this method possesses good computational efficiency because it only requires first-order derivatives. In addition, it can easily be accelerated by using multiple graphics processing units and central processing units, for weight updating and forward modeling. Synthetic experiments, based on the Marmousi2, 2004 BP, and a metal ore model, are used to show the numerical performance of the deep-learning inversion method. Comprehensive comparisons with a total-variation regularized FWI are presented to show the ability of our method to recover sharp boundaries. Our numerical results indicate that this deep-learning inversion approach is effective, efficient, and can capture salient features of the model.

Author(s):  
Marcos Bernal-Romero ◽  
Ursula Iturrarán-Viveros

Summary Full-Waveform Inversion (FWI) is a procedure based on the minimization of a misfit (or cost) function applied to the difference between synthetic waveforms and real seismic traces that derives high-resolution velocity models. This is achieved through the iterative adjustment of the velocity model and/or some other physical parameters of the Earth’s subsurface, which generally implies large computational effort. In order to minimize this cost function we explore the use of Adaptive Gradient Optimization (AGO), a variant of Stochastic Gradient Descent (SGD) methods, combining them with a dynamic simultaneous sources strategy that allow us to reduce the computational cost involved in this process. AGO methods are computationally efficient, have little memory requirements and have the capability of adapting the step-length according to the optimization process’ evolution. Since a precise calibration of the step-length is needed to ensure efficiency, the AGOs are well-suited for this task because they are able to adapt the step-length according to the optimization’s development. In this work, we propose a simple non-linear relationship that allows an adjustment of the step-length with respect to the frequencies used in the multiscale FWI, avoiding the line-search strategy’s high computational burden. Additionally, the application of this new step-length rule into the AGO methods with a dynamic simultaneous sources strategy, allow us to concurrently accelerate and significantly improve the FWI’s numerical performance and results. We compare the performance and final results of seven AGO methods, using two different FWI misfit functionals (based on L1 and L2 norms) applied to estimate the final velocity models of two benchmark acoustic models: the Marmousi and the Canadian overthrust BP velocity models.


Symmetry ◽  
2021 ◽  
Vol 13 (2) ◽  
pp. 260
Author(s):  
Meng Suo ◽  
Dong Zhang ◽  
Yan Yang

Inspired by the large number of applications for symmetric nonlinear equations, an improved full waveform inversion algorithm is proposed in this paper in order to quantitatively measure the bone density and realize the early diagnosis of osteoporosis. The isotropic elastic wave equation is used to simulate ultrasonic propagation between bone and soft tissue, and the Gauss–Newton algorithm based on symmetric nonlinear equations is applied to solve the optimal solution in the inversion. In addition, the authors use several strategies including the frequency-grid multiscale method, the envelope inversion and the new joint velocity–density inversion to improve the result of conventional full-waveform inversion method. The effects of various inversion settings are also tested to find a balanced way of keeping good accuracy and high computational efficiency. Numerical inversion experiments showed that the improved full waveform inversion (FWI) method proposed in this paper shows superior inversion results as it can detect small velocity–density changes in bones, and the relative error of the numerical model is within 10%. This method can also avoid interference from small amounts of noise and satisfy the high precision requirements for quantitative ultrasound measurements of bone.


2021 ◽  
Author(s):  
Yao Liu ◽  
Baodi Liu ◽  
Jianping Huang ◽  
Jun Wang ◽  
Honglong Chen ◽  
...  

Geophysics ◽  
2021 ◽  
Vol 86 (1) ◽  
pp. R31-R44
Author(s):  
Bin Liu ◽  
Senlin Yang ◽  
Yuxiao Ren ◽  
Xinji Xu ◽  
Peng Jiang ◽  
...  

Velocity model inversion is one of the most important tasks in seismic exploration. Full-waveform inversion (FWI) can obtain the highest resolution in traditional velocity inversion methods, but it heavily depends on initial models and is computationally expensive. In recent years, a large number of deep-learning (DL)-based velocity model inversion methods have been proposed. One critical component in those DL-based methods is a large training set containing different velocity models. We have developed a method to construct a realistic structural model for the DL network. Our compressional-wave velocity model building method for creating dense-layer/fault/salt body models can automatically construct a large number of models without much human effort, which is very meaningful for DL networks. Moreover, to improve the inversion result on these realistic structural models, instead of only using the common-shot gather, we also extract features from the common-receiver gather as well. Through a large number of realistic structural models, reasonable data acquisition methods, and appropriate network setups, a more generalized result can be obtained through our proposed inversion framework, which has been demonstrated to be effective on the independent testing data set. The results of dense-layer models, fault models, and salt body models that we compared and analyzed demonstrate the reliability of our method and also provide practical guidelines for choosing optimal inversion strategies in realistic situations.


2015 ◽  
Author(s):  
Haishan Li* ◽  
Wuyang Yang ◽  
Enli Wang ◽  
Xueshan Yong

2020 ◽  
Vol 23 (4) ◽  
pp. 347-358
Author(s):  
Boyoung Kim ◽  
Jun Won Kang ◽  
Yeong-Tae Choi ◽  
Seung Yup Jang

Geophysics ◽  
2021 ◽  
pp. 1-74
Author(s):  
Zhen Zhou ◽  
Anja Klotzsche ◽  
Jessica Schmäck ◽  
Harry Vereecken ◽  
Jan van der Kruk

Detailed characterization of aquifers is critical and challenging due to the existence of heterogeneous small-scale high-contrast layers. For an improved characterization of subsurface hydrological characteristics, crosshole ground penetrating radar (GPR) and Cone Penetration Test (CPT) measurements are performed. In comparison to the CPT approach that can only provide 1D high resolution data along vertical profiles, crosshole GPR enables measuring 2D cross-sections between two boreholes. Generally, a standard inversion method for GPR data is the ray-based approach that considers only a small amount of information and can therefore only provide limited resolution. In the last decade, full-waveform inversion (FWI) of crosshole GPR data in time domain has matured, and provides inversion results with higher resolution by exploiting the full recorded waveform information. However, the FWI results are limited due to complex underground structures and the non-linear nature of the method. A new approach that uses CPT data in the inversion process is applied to enhance the resolution of the final relative permittivity FWI results by updating the effective source wavelet. The updated effective source wavelet possesses a priori CPT information and a larger bandwidth. Using the same starting models, a synthetic model comparison between the conventional and updated FWI results demonstrates that the updated FWI method provides reliable and more consistent structures. To test the method, five experimental GPR cross-section results are analyzed with the standard FWI and the new proposed updated approach. Both synthetic and experimental results indicate the potential of improving the reconstruction of subsurface aquifer structures by combining conventional 2D FWI results and 1D CPT data.


Sign in / Sign up

Export Citation Format

Share Document