Frequency domain least-squares reverse time migration with a modified scattering-integral approach

Author(s):  
Jizhong Yang* ◽  
Yuzhu Liu ◽  
Liangguo Dong
Geophysics ◽  
2016 ◽  
Vol 81 (6) ◽  
pp. S497-S509 ◽  
Author(s):  
Jizhong Yang ◽  
Yuzhu Liu ◽  
Liangguo Dong

Least-squares migration (LSM) is commonly regarded as an amplitude-preserving or true amplitude migration algorithm that, compared with conventional migration, can provide migrated images with reduced migration artifacts, balanced amplitudes, and enhanced spatial resolution. Most applications of LSM are based on the constant-density assumption, which is not the case in the real earth. Consequently, the amplitude performance of LSM is not appropriate. To partially remedy this problem, we have developed a least-squares reverse time migration (LSRTM) scheme suitable for density variations in the acoustic approximation. An improved scattering-integral approach is adopted for implementation of LSRTM in the frequency domain. LSRTM images associated with velocity and density perturbations are simultaneously used to generate the simulated data, which better matches the recorded data in amplitudes. Summation of these two images provides a reflectivity model related to impedance perturbation that is in better accordance with the true one, than are the velocity and density images separately. Numerical examples based on a two-layer model and a small part of the Sigsbee2A model verify the effectiveness of our method.


2021 ◽  
Vol 1719 (1) ◽  
pp. 012030
Author(s):  
Phudit Sombutsirinun ◽  
Chaiwoot Boonyasiriwat

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.


Sign in / Sign up

Export Citation Format

Share Document