scholarly journals Near surface reconstruction using the full waveform inversion method: a resolution analysis

Author(s):  
К.Г. Гадыльшин ◽  
Д.А. Неклюдов ◽  
И.Ю. Сильвестров

Изучается влияние волны Рэлея на результат обращения полного волнового поля. Показано, что разрешающая способность борновского оператора выше в случае отсутствия рэлеевских волн в наблюдаемых данных. Предложен численный алгоритм, позволяющий анализировать влияние различных параметров обратной динамической задачи сейсмики на результат восстановления верхней части геологического разреза: геометрии системы наблюдений, частотного диапазона, параметризации упругой модели и др. В качестве демонстрационного примера показано, что для поверхностной системы наблюдений в задаче реконструкции макроскоростного строения среды наиболее приемлема параметризация упругой модели через скорости продольных и поперечных волн и плотность. The effect of the Rayleigh wave on full waveform inversion (FWI) is studied. It is shown that the Born operator resolution is better when the Rayleigh waves are not taken into account. A numerical algorithm to analyze the influence of different parameters of the inverse dynamical seismic problem (such as source/receiver geometry, frequency range, elastic model parametrization, etc.) on the near-surface geological section reconstruction is proposed. As a demonstrative example, we show that the best elastic parameters for macrovelocity reconstruction via FWI for a subsurface observation system are P- and S-velocities and density.

2017 ◽  
Vol 5 (3) ◽  
pp. SL89-SL107 ◽  
Author(s):  
Arthur B. Weglein

Direct inverse methods solve the problem of interest; in addition, they communicate whether the problem of interest is the problem that we (the seismic industry) need to be interested in. When a direct solution does not result in an improved drill success rate, we know that the problem we have chosen to solve is not the right problem — because the solution is direct and cannot be the issue. On the other hand, with an indirect method, if the result is not an improved drill success rate, then the issue can be either the chosen problem, or the particular choice within the plethora of indirect solution methods, or both. The inverse scattering series (ISS) is the only direct inversion method for a multidimensional subsurface. Solving a forward problem in an inverse sense is not equivalent to a direct inverse solution. All current methods for parameter estimation, e.g., amplitude-variation-with-offset and full-waveform inversion, are solving a forward problem in an inverse sense and are indirect inversion methods. The direct ISS method for determining earth material properties defines the precise data required and the algorithms that directly output earth mechanical properties. For an elastic model of the subsurface, the required data are a matrix of multicomponent data, and a complete set of shot records, with only primaries. With indirect methods, any data can be matched: one trace, one or several shot records, one component, multicomponent, with primaries only or primaries and multiples. Added to that are the innumerable choices of cost functions, generalized inverses, and local and global search engines. Direct and indirect parameter inversion are compared. The direct ISS method has more rapid convergence and a broader region of convergence. The difference in effectiveness increases as subsurface circumstances become more realistic and complex, in particular with band-limited noisy data.


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.


Geophysics ◽  
2021 ◽  
pp. 1-37
Author(s):  
Xinhai Hu ◽  
Wei Guoqi ◽  
Jianyong Song ◽  
Zhifang Yang ◽  
Minghui Lu ◽  
...  

Coupling factors of sources and receivers vary dramatically due to the strong heterogeneity of near surface, which are as important as the model parameters for the inversion success. We propose a full waveform inversion (FWI) scheme that corrects for variable coupling factors while updating the model parameter. A linear inversion is embedded into the scheme to estimate the source and receiver factors and compute the amplitude weights according to the acquisition geometry. After the weights are introduced in the objective function, the inversion falls into the category of separable nonlinear least-squares problems. Hence, we could use the variable projection technique widely used in source estimation problem to invert the model parameter without the knowledge of source and receiver factors. The efficacy of the inversion scheme is demonstrated with two synthetic examples and one real data test.


Geophysics ◽  
2019 ◽  
Vol 84 (1) ◽  
pp. B95-B105 ◽  
Author(s):  
Yao Wang ◽  
Richard D. Miller ◽  
Shelby L. Peterie ◽  
Steven D. Sloan ◽  
Mark L. Moran ◽  
...  

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.


2017 ◽  
Author(s):  
Yao Wang ◽  
Richard Miller ◽  
Shelby Peterie ◽  
Steven Sloan ◽  
Mark Moran ◽  
...  

Geophysics ◽  
2018 ◽  
Vol 83 (1) ◽  
pp. R1-R11 ◽  
Author(s):  
Dmitry Borisov ◽  
Ryan Modrak ◽  
Fuchun Gao ◽  
Jeroen Tromp

Full-waveform inversion (FWI) is a powerful method for estimating the earth’s material properties. We demonstrate that surface-wave-driven FWI is well-suited to recovering near-surface structures and effective at providing S-wave speed starting models for use in conventional body-wave FWI. Using a synthetic example based on the SEG Advanced Modeling phase II foothills model, we started with an envelope-based objective function to invert for shallow large-scale heterogeneities. Then we used a waveform-difference objective function to obtain a higher-resolution model. To accurately model surface waves in the presence of complex tomography, we used a spectral-element wave-propagation solver. Envelope misfit functions are found to be effective at minimizing cycle-skipping issues in surface-wave inversions, and surface waves themselves are found to be useful for constraining complex near-surface features.


Sign in / Sign up

Export Citation Format

Share Document