scholarly journals A direct inverse method for subsurface properties: The conceptual and practical benefit and added value in comparison with all current indirect methods, for example, amplitude-variation-with-offset and full-waveform inversion

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.

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.


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