Guided wave tomography based on least-squares reverse-time migration

2019 ◽  
Vol 19 (4) ◽  
pp. 1237-1249 ◽  
Author(s):  
Jiaze He ◽  
Daniel C Rocha ◽  
Paul Sava

A key to successful damage diagnostics and quantification is damage imaging through ultrasonic guided wave tomography. We propose the implementation of least-squares reverse-time migration in a circular array for damage imaging in an aluminum plate. The theory of least-squares reverse-time migration is formulated for guided wave applications along with the summary of an efficient optimization algorithm: the conjugate gradient method. Numerical simulation and laboratory experiments are used to evaluate its performance with a circular array setup. In order to improve the data processing efficiency, the concept of using a limited number of actuators but a relatively large number of sensors is tested. Studies are conducted on three numerical cases, including a rectangular-shaped damage site, a complex-shaped damage site, and six other damage sites varying in size. As an inversion-based method, least-squares reverse-time migration shows significantly improved shape reconstruction with the amplitude quantification capability, compared to conventional reverse-time migration. Our experimental data are generated by piezoelectric wafers as actuators, measured by a scanning laser Doppler vibrometer to form a circular array on an aluminum plate, with a rectangular notch located in the inner region of the array. The damage images using experimental data show consistency in both the simulations using Born scattering and in altered material properties in the damaged region. According to the comparison, least-squares reverse-time migration for guided wave tomography is a promising technology to provide high-resolution large area damage imaging for plate-like structures.

2019 ◽  
Vol 28 (6) ◽  
pp. 065010 ◽  
Author(s):  
Jiaze He ◽  
Daniel C Rocha ◽  
Patrick E Leser ◽  
Paul Sava ◽  
William P Leser

Geophysics ◽  
2021 ◽  
pp. 1-73
Author(s):  
Milad Farshad ◽  
Hervé Chauris

Elastic least-squares reverse time migration is the state-of-the-art linear imaging technique to retrieve high-resolution quantitative subsurface images. A successful application requires many migration/modeling cycles. To accelerate the convergence rate, various pseudoinverse Born operators have been proposed, providing quantitative results within a single iteration, while having roughly the same computational cost as reverse time migration. However, these are based on the acoustic approximation, leading to possible inaccurate amplitude predictions as well as the ignorance of S-wave effects. To solve this problem, we extend the pseudoinverse Born operator from acoustic to elastic media to account for the elastic amplitudes of PP reflections and provide an estimate of physical density, P- and S-wave impedance models. We restrict the extension to marine environment, with the recording of pressure waves at the receiver positions. Firstly, we replace the acoustic Green's functions by their elastic version, without modifying the structure of the original pseudoinverse Born operator. We then apply a Radon transform to the results of the first step to calculate the angle-dependent response. Finally, we simultaneously invert for the physical parameters using a weighted least-squares method. Through numerical experiments, we first illustrate the consequences of acoustic approximation on elastic data, leading to inaccurate parameter inversion as well as to artificial reflector inclusion. Then we demonstrate that our method can simultaneously invert for elastic parameters in the presence of complex uncorrelated structures, inaccurate background models, and Gaussian noisy data.


Geophysics ◽  
2021 ◽  
pp. 1-65
Author(s):  
Yingming Qu ◽  
Yixin Wang ◽  
Zhenchun Li ◽  
Chang Liu

Seismic wave attenuation caused by subsurface viscoelasticity reduces the quality of migration and the reliability of interpretation. A variety of Q-compensated migration methods have been developed based on the second-order viscoacoustic quasidifferential equations. However, these second-order wave-equation-based methods are difficult to handle with density perturbation and surface topography. In addition, the staggered grid scheme, which has an advantage over the collocated grid scheme because of its reduced numerical dispersion and enhanced stability, works in first-order wave-equation-based methods. We have developed a Q least-squares reverse time migration method based on the first-order viscoacoustic quasidifferential equations by deriving Q-compensated forward-propagated operators, Q-compensated adjoint operators, and Q-attenuated Born modeling operators. Besides, our method using curvilinear grids is available even when the attenuating medium has surface topography and can conduct Q-compensated migration with density perturbation. The results of numerical tests on two synthetic and a field data sets indicate that our method improves the imaging quality with iterations and produces better imaging results with clearer structures, higher signal-to-noise ratio, higher resolution, and more balanced amplitude by correcting the energy loss and phase distortion caused by Q attenuation. It also suppresses the scattering and diffracted noise caused by the surface topography.


Geophysics ◽  
2018 ◽  
Vol 83 (1) ◽  
pp. S33-S46 ◽  
Author(s):  
Chuang Li ◽  
Jianping Huang ◽  
Zhenchun Li ◽  
Rongrong Wang

This study derives a preconditioned stochastic conjugate gradient (CG) method that combines stochastic optimization with singular spectrum analysis (SSA) denoising to improve the efficiency and image quality of plane-wave least-squares reverse time migration (PLSRTM). This method reduces the computational costs of PLSRTM by applying a controlled group-sampling method to a sufficiently large number of plane-wave sections and accelerates the convergence using a hybrid of stochastic descent (SD) iteration and CG iteration. However, the group sampling also produces aliasing artifacts in the migration results. We use SSA denoising as a preconditioner to remove the artifacts. Moreover, we implement the preconditioning on the take-off angle-domain common-image gathers (CIGs) for better results. We conduct numerical tests using the Marmousi model and Sigsbee2A salt model and compare the results of this method with those of the SD method and the CG method. The results demonstrate that our method efficiently eliminates the artifacts and produces high-quality images and CIGs.


Sign in / Sign up

Export Citation Format

Share Document