Bayesian full waveform inversion in anisotropic elastic media using the iterated extended Kalman filter

Geophysics ◽  
2020 ◽  
pp. 1-69
Author(s):  
Xingguo Huang ◽  
Kjersti Solberg Eikrem ◽  
Morten Jakobsen ◽  
Geir Naevdal

Uncertainty quantification in the context of seismic imaging is important for interpretinginverted subsurface models and updating reservoir models. The limited illumination, noisydata and poor initial model in the seismic full waveform inversion (FWI) lead to inversionuncertainties. This is particularly true for anisotropic elastic FWI, which suffers from extra parameter trade-off problems. In this work, we address the uncertainty quantificationof anisotropic elastic FWI problem in the framework of Bayesian inference. Specially, weestimate the uncertainties of the subsurface elastic parameters in the Bayesian anisotropicelastic FWI by combining the iterated extended Kalman filter with an explicit representation of the sensitivity matrix with Green’s functions. The sensitivity matrix is based onthe integral equation approach, which is also within the context of nonlinear inverse scattering theory. We give the results of numerical tests with examples for anisotropic elasticmedia. They show that the proposed Bayesian inversion method can provide reasonablereconstructed results for the elastic coefficients of the stiffness tensor and the framework issuitable for accessing the uncertainties.

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.


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.


2018 ◽  
Vol 8 (2) ◽  
Author(s):  
Katherine Flórez ◽  
Sergio Alberto Abreo Carrillo ◽  
Ana Beatriz Ramírez Silva

Full Waveform Inversion (FWI) schemes are gradually becoming more common in the oil and gas industry, as a new tool for studying complex geological zones, based on their reliability for estimating velocity models. FWI is a non-linear inversion method that iteratively estimates subsurface characteristics such as seismic velocity, starting from an initial velocity model and the preconditioned data acquired. Blended sources have been used in marine seismic acquisitions to reduce acquisition costs, reducing the number of times that the vessel needs to cross the exploration delineation trajectory. When blended or simultaneous without previous de-blending or separation, stage data are used in the reconstruction of the velocity model with the FWI method, and the computational time is reduced. However, blended data implies overlapping single shot-gathers, producing interference that affects the result of seismic approaches, such as FWI or seismic image migration. In this document, an encoding strategy is developed, which reduces the overlap areas within the blended data to improve the final velocity model with the FWI method.


Sign in / Sign up

Export Citation Format

Share Document