Application of ground penetrating radar to detect tunnel lining defects based on improved full waveform inversion and reverse time migration

To obtain a physical understanding of gradient-based descent methods in full-waveform inversion (FWI), we find a connection between the FWI gradient and the image provided by reverse time migration (RTM). The gradient uses the residual data as a virtual source, and RTM uses the observed data directly as the boundary condition, so the FWI gradient is similar to a time integration of the RTM image using the residual data, which physically converts the phase of the reflectivity to the phase of the velocity. Therefore, gradient-based FWI can be connected to the classical reflectivity-to-velocity/impedance inversion (RVI). We have developed a new FWI scheme that provides a self-contained and physically intuitive derivation, which naturally establishes a connection among the amplitude-preserved RTM, the Zoeppritz equations (amplitude variation with angle inversion), and RVI, and combines them into a single framework to produce a preconditioned inversion formula. In this scheme, the relative velocity update is a phase-modified and deconvolved RTM image obtained with the residual data. Consistent with the deconvolution, the multiscale approach applies a gradually widening low-pass frequency filter to the deconvolved wavelet at early iterations, and then it uses the unfiltered deconvolved wavelet for the final iterations. Our numerical testing determined that the new method makes a significant improvement to the quality of the inversion result.

Download Full-text

Full-waveform inversion and least-squares reverse-time migration imaging of collimated ultrasonic-beam data for high-resolution wellbore integrity monitoring

Applied Physics Letters ◽

10.1063/1.5034122 ◽

2018 ◽

Vol 113 (7) ◽

pp. 071903 ◽

Cited By ~ 1

Author(s):

Yu Chen ◽

Kai Gao ◽

Eric S. Davis ◽

Dipen N. Sinha ◽

Cristian Pantea ◽

...

Keyword(s):

High Resolution ◽

Waveform Inversion ◽

Full Waveform Inversion ◽

Reverse Time ◽

Reverse Time Migration ◽

Ultrasonic Beam ◽

Wellbore Integrity ◽

Full Waveform ◽

Time Migration ◽

Migration Imaging

Download Full-text

Multisource Reverse-time Migration and Full-waveform Inversion on a GPGPU

72nd EAGE Conference and Exhibition incorporating SPE EUROPEC 2010 ◽

10.3997/2214-4609.201401241 ◽

2010 ◽

Cited By ~ 2

Author(s):

C. Boonyasiriwat ◽

G. Zhan ◽

M. Hadwiger ◽

M. Srinivasan ◽

G. T. Schuster

Keyword(s):

Waveform Inversion ◽

Full Waveform Inversion ◽

Reverse Time ◽

Reverse Time Migration ◽

Full Waveform ◽

Time Migration

Download Full-text

Efficient snapshot-free reverse time migration and computation of multiparameter gradients in full waveform inversion

Geophysics ◽

10.1190/geo2020-0606.1 ◽

2021 ◽

pp. 1-79

Author(s):

Johan O. A. Robertsson ◽

Fredrik Andersson ◽

René-Édouard Plessix

Keyword(s):

Waveform Inversion ◽

Full Waveform Inversion ◽

Reverse Time ◽

Reverse Time Migration ◽

Principle Of Superposition ◽

Model Parameter ◽

Inversion Model ◽

Full Waveform ◽

Time Migration ◽

Anisotropic Elastic

Computing images in reverse time migration and model parameter gradients from adjoint wavefields in full waveform inversion requires the correlation of a forward propagated wavefield with another reverse propagated wavefield. Although in theory only two wavefield propagations are required, one forward propagation and one reverse propagation, it requires storing the forward propagated wavefield as a function of time to carry out the correlations which is associated with significant I/O cost. Alternatively, three wavefield propagations can be carried out to reverse propagate the forward propagated wavefield in tandem with the reverse propagated wavefield. We show how highly accurate reverse time migrated images and full waveform inversion model parameter gradients for anisotropic elastic full waveform inversion can be efficiently computed without significant disk I/O using two wavefield propagations by means of the principle of superposition.

Download Full-text

Breakthrough acquisition and technologies for subsalt imaging

Geophysics ◽

10.1190/geo2010-0399.1 ◽

2011 ◽

Vol 76 (5) ◽

pp. WB41-WB51 ◽

Cited By ~ 42

Author(s):

Denes Vigh ◽

Jerry Kapoor ◽

Nick Moldoveanu ◽

Hongyan Li

Keyword(s):

Velocity Field ◽

Inversion Layer ◽

Waveform Inversion ◽

Full Waveform Inversion ◽

Reverse Time ◽

Reverse Time Migration ◽

Data Set ◽

Low Frequencies ◽

Full Waveform ◽

Time Migration

The recently introduced method of wide-azimuth data acquisition offers better illumination, noise attenuation, and lower frequencies to more accurately determine a velocity field for imaging. For the field data experiment to demonstrate the technologies, we used a Gulf of Mexico (GOM) wide-azimuth data set that allows us to take advantage of possible low frequencies, relatively large crossline offsets, and increased illumination. The input data was processed through true 3D azimuthal surface-related multiple elimination (SRME) with zero-phasing and debubble. Eliminating the surface-related multiples aids the velocity determination and helps uncover the subsalt sediments at the final imaging stage. After the initial velocity derivation, which was constrained to wells and geology, full-waveform inversion (FWI) was used to further update the velocity field to achieve an enhanced image. The methodology used follows the top-down approach where suprasalt sediment model is developed followed by the top of salt, salt flanks, base of salt, and finished with a limited subsalt update. To approximate the observed data by using an acoustic inversion procedure, the propagator includes effects of attenuation, anisotropy, acquisition source, and receiver depth. The geological environment is salt related, which implies that the observed data is highly elastic, even though it is input to an acoustic full waveform inversion. To use the proper constraints for the inversion, layer-stripping method is used to develop the high-resolution velocity field. The inversion stages were carefully quality controlled through gather displays to ensure the kinematics were honored. We then demonstrated the benefit of the FWI velocity field by comparing the images derived with the traditional ray-based tomographic velocity field versus the velocity field derived by FWI performing reverse time migration to produce these images. Finally, the images were compared at key well locations to evaluate the robustness of the workflow.

Download Full-text

3D high-resolution seismic imaging of the iron-oxide deposits in Ludvika (Sweden) using full-waveform inversion and reverse-time migration

10.5194/se-2021-122 ◽

2021 ◽

Author(s):

Brij Singh ◽

Michał Malinowski ◽

Andrzej Górszczyk ◽

Alireza Malehmir ◽

Stefan Buske ◽

...

Keyword(s):

Waveform Inversion ◽

Velocity Model ◽

Mine Planning ◽

Full Waveform Inversion ◽

Reverse Time ◽

Reverse Time Migration ◽

Depth Imaging ◽

Velocity Models ◽

Full Waveform ◽

Time Migration

Abstract. A sparse 3D seismic survey was acquired over the Blötberget iron-oxide deposits of the Ludvika Mines in south-central Sweden. The main aim of the survey was to delineate the deeper extension of the mineralisation and to better understand its 3D nature and associated fault systems for mine planning purposes. To obtain a high-quality seismic image in depth, we applied time-domain 3D acoustic full-waveform inversion (FWI) to build a high-resolution P-wave velocity model. This model was subsequently used for pre-stack depth imaging with reverse time migration (RTM) to produce the complementary reflectivity section. We developed a data preprocessing workflow and inversion strategy for the successful implementation of FWI in the hardrock environment. We obtained a high-fidelity velocity model using FWI and assessed its robustness. We extensively tested and optimised the parameters associated with the RTM method for subsequent depth imaging using different velocity models: a constant velocity model, a model built using first-arrival traveltime tomography and a velocity model derived by FWI. We compare our RTM results with a priori data available in the area. We conclude that, from all tested velocity models, the FWI velocity model in combination with the subsequent RTM step, provided the most focussed image of the mineralisation and we successfully mapped its 3D geometrical nature. In particular, a major reflector interpreted as a cross-cutting fault, which is restricting the deeper extension of the mineralisation with depth, and several other fault structures which were earlier not imaged were also delineated. We believe that a thorough analysis of the depth images derived with the combined FWIRTM approach that we presented here can provide more details which will help with better estimation of areas with high mineralization, better mine planning and safety measures.

Download Full-text