Data Reconstruction Inversion: An Augmented Lagrangian Based Full-Waveform Inversion in the Time Domain

Author(s):  
A. Gholami ◽  
H. Aghamiry ◽  
S. Operto





Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE135-VE144 ◽  
Author(s):  
Denes Vigh ◽  
E. William Starr

Prestack depth migration has been used for decades to derive velocity distributions in depth. Numerous tools and methodologies have been developed to reach this goal. Exploration in geologically more complex areas exceeds the abilities of existing methods. New data-acquisition and data-processing methods are required to answer these new challenges effectively. The recently introduced wide-azimuth data acquisition method offers better illumination and noise attenuation as well as an opportunity to more accurately determine velocities for imaging. One of the most advanced tools for depth imaging is full-waveform inversion. Prestack seismic full-waveform inversion is very challenging because of the nonlinearity and nonuniqueness of the solution. Combined with multiple iterations of forward modeling and residual wavefield back propagation, the method is computer intensive, especially for 3D projects. We studied a time-domain, plane-wave implementation of 3D waveform inversion. We found that plane-wave gathers are an attractive input to waveform inversion with dramatically reduced computer run times compared to traditional shot-gather approaches. The study was conducted on two synthetic data sets — Marmousi2 and SMAART Pluto 1.5 — and a field data set. The results showed that a velocity field can be reconstructed well using a multiscale time-domain implementation of waveform inversion. Although the time-domain solution does not take advantage of wavenumber redundancy, the method is feasible on current computer architectures for 3D surveys. The inverted velocity volume produces a quality image for exploration geologists by using numerous iterations of waveform inversion.





Geophysics ◽  
2021 ◽  
pp. 1-57
Author(s):  
Ali Gholami ◽  
Hossein S. Aghamiry ◽  
Stéphane Operto

The search space of Full Waveform Inversion (FWI) can be extended via a relaxation of the wave equation to increase the linear regime of the inversion. This wave equation relaxation is implemented by solving jointly (in a least-squares sense) the wave equation weighted by a penalty parameter and the observation equation such that the reconstructed wavefields closely match the data, hence preventing cycle skipping at receivers. Then, the subsurface parameters are updated by minimizing the temporal and spatial source extension generated by the wave-equation relaxation to push back the data-assimilated wavefields toward the physics.This extended formulation of FWI has been efficiently implemented in the frequency domain with the augmented Lagrangian method where the overdetermined systems of the data-assimilated wavefields can be solved separately for each frequency with linear algebra methods and the sensitivity of the optimization to the penalty parameter is mitigated through the action of the Lagrange multipliers.Applying this method in the time domain is however hampered by two main issues: the computation of data-assimilated wavefields with explicit time-stepping schemes and the storage of the Lagrange multipliers capturing the history of the source residuals in the state space.These two issues are solved by recognizing that the source residuals on the right-hand side of the extended wave equation, when formulated in a form suitable for explicit time stepping, are related to the extended data residuals through an adjoint equation.This relationship first allows us to relate the extended data residuals to the reduced data residuals through a normal equation in the data space. Once the extended data residuals have been estimated by solving (exactly or approximately) this normal equation, the data-assimilated wavefields are computed with explicit time stepping schemes by cascading an adjoint and a forward simulation.



Sensors ◽  
2018 ◽  
Vol 18 (9) ◽  
pp. 3114 ◽  
Author(s):  
Sixin Liu ◽  
Xintong Liu ◽  
Xu Meng ◽  
Lei Fu ◽  
Qi Lu ◽  
...  

Xiuyan Jade, produced in Xiuyan County, Liaoning Province, China is one of the four famous jade in China. King Jade, which is deemed the largest jade body of the world, was broken out from a hill. The local government planned to build a tourism site based on the jade culture there. The purpose of the investigation was to evaluate the stability of subsurface foundation, and the possible positions of mined-out zones to prevent the further rolling of the jade body. Cross-hole radar tomography is the key technique in the investigation. Conventional travel time and attenuation tomography based on ray tracing theory cannot provide high-resolution images because only a fraction of the measured information is used in the inversion. Full-waveform inversion (FWI) can provide high-resolution permittivity and conductivity images because it utilizes all the information provided by the radar signals. We deduce the gradient expression of the time-domain FWI with respect to the permittivity and conductivity using a method that is different from that of the previous work and realize the FWI algorithm that can simultaneously update the permittivity and conductivity by using the conjugate gradient method. Inverted results from synthetic data show that time-domain FWI can significantly improve the resolution compared with the ray-based tomogram methods. FWI can distinguish targets that are as small as one-half to one-third wavelength and the inverted physical values are closer to the real ones than those provided by the ray tracing method. We use the FWI algorithm to the field data measured at Xiuyan jade mine. Both the inverted permittivity and conductivity can comparably delineate four mined-out zones, which exhibit low-permittivity and low-conductivity characteristics. Furthermore, the locations of the interpreted mined-out zones are in good agreement with the existing mining channels recorded by geological data.



Sign in / Sign up

Export Citation Format

Share Document