scholarly journals Q-least-squares reverse time migration with viscoacoustic deblurring filters

Geophysics ◽  
2017 ◽  
Vol 82 (6) ◽  
pp. S425-S438 ◽  
Author(s):  
Yuqing Chen ◽  
Gaurav Dutta ◽  
Wei Dai ◽  
Gerard T. Schuster
Keyword(s):  
Convergence Rate ◽  
Least Squares ◽  
Field Data ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Acoustic Data ◽  
Numerical Tests ◽  
Time Migration ◽  
Preconditioning Method ◽  
Number Of Iterations

Viscoacoustic least-squares reverse time migration, also denoted as Q-LSRTM, linearly inverts for the subsurface reflectivity model from lossy data. Compared with conventional migration methods, it can compensate for the amplitude loss in the migrated images due to strong subsurface attenuation and can produce reflectors that are accurately positioned in depth. However, the adjoint [Formula: see text] propagators used for backward propagating the residual data are also attenuative. Thus, the inverted images from [Formula: see text]-LSRTM with a small number of iterations are often observed to have lower resolution when compared with the benchmark acoustic LSRTM images from acoustic data. To increase the resolution and accelerate the convergence of [Formula: see text]-LSRTM, we used viscoacoustic deblurring filters as a preconditioner for [Formula: see text]-LSRTM. These filters can be estimated by matching a simulated migration image to its reference reflectivity model. Numerical tests on synthetic and field data demonstrate that [Formula: see text]-LSRTM combined with viscoacoustic deblurring filters can produce images with higher resolution and more balanced amplitudes than images from acoustic RTM, acoustic LSRTM, and [Formula: see text]-LSRTM when there is strong attenuation in the background medium. Our preconditioning method is also shown to improve the convergence rate of [Formula: see text]-LSRTM by more than 30% in some cases and significantly compensate for the lossy artifacts in RTM images.

Plane-wave least-squares reverse time migration with a preconditioned stochastic conjugate gradient method

Geophysics ◽  
2018 ◽  
Vol 83 (1) ◽  
pp. S33-S46 ◽  
Author(s):  
Chuang Li ◽  
Jianping Huang ◽  
Zhenchun Li ◽  
Rongrong Wang
Keyword(s):  
Plane Wave ◽  
Least Squares ◽  
Conjugate Gradient ◽  
Sampling Method ◽  
Singular Spectrum Analysis ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Numerical Tests ◽  
Time Migration

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.

scholarly journals Q-compensated least-squares reverse time migration using low-rank one-step wave extrapolation

Geophysics ◽  
2016 ◽  
Vol 81 (4) ◽  
pp. S271-S279 ◽  
Author(s):  
Junzhe Sun ◽  
Sergey Fomel ◽  
Tieyuan Zhu ◽  
Jingwei Hu
Keyword(s):  
Convergence Rate ◽  
Least Squares ◽  
Seismic Attenuation ◽  
Low Rank ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Attenuation Compensation ◽  
Time Migration ◽  
Rank One ◽  
One Step

Attenuation of seismic waves needs to be taken into account to improve the accuracy of seismic imaging. In viscoacoustic media, reverse time migration (RTM) can be performed with [Formula: see text]-compensation, which is also known as [Formula: see text]-RTM. Least-squares RTM (LSRTM) has also been shown to be able to compensate for attenuation through linearized inversion. However, seismic attenuation may significantly slow down the convergence rate of the least-squares iterative inversion process without proper preconditioning. We have found that incorporating attenuation compensation into LSRTM can improve the speed of convergence in attenuating media, obtaining high-quality images within the first few iterations. Based on the low-rank one-step seismic modeling operator in viscoacoustic media, we have derived its adjoint operator using nonstationary filtering theory. The proposed forward and adjoint operators can be efficiently applied to propagate viscoacoustic waves and to implement attenuation compensation. Recognizing that, in viscoacoustic media, the wave-equation Hessian may become ill-conditioned, we propose to precondition LSRTM with [Formula: see text]-compensated RTM. Numerical examples showed that the preconditioned [Formula: see text]-LSRTM method has a significantly faster convergence rate than LSRTM and thus is preferable for practical applications.

Improved subsalt images with least-squares reverse time migration

2017 ◽  
Vol 5 (3) ◽  
pp. SN25-SN32 ◽  
Author(s):  
Ping Wang ◽  
Shouting Huang ◽  
Ming Wang
Keyword(s):  
Least Squares ◽  
Field Data ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Lateral Velocity ◽  
Image Domain ◽  
Time Migration ◽  
Iterative Inversion ◽  
Velocity Variations ◽  
Subsalt Imaging

Complex overburdens often distort reservoir images in terms of structural positioning, stratigraphic resolution, and amplitude fidelity. One prime example of a complex overburden is in the deepwater Gulf of Mexico, where thick and irregular layers of remobilized (i.e., allochthonous) salt are situated above prospective reservoir intervals. The highly variant salt layers create large lateral velocity variations that distort wave propagation and the illumination of deeper reservoir targets. In subsalt imaging, tools such as reflection tomography, full-waveform inversion, and detailed salt interpretation are needed to derive a high-resolution velocity model that captures the lateral velocity variations. Once a velocity field is obtained, reverse time migration (RTM) can be applied to restore structural positioning of events below and around the salt. However, RTM by nature is unable to fully recover the reflectivity for desired amplitudes and resolution. This shortcoming is well-recognized by the imaging community, and it has propelled the emergence of least-squares RTM (LSRTM) in recent years. We have investigated how current LSRTM methods perform on subsalt images. First, we compared the formulation of data-domain versus image-domain least-squares migration, as well as methods using single-iteration approximation versus iterative inversion. Then, we examined the resulting subsalt images of several LSRTM methods applied on the synthetic and field data. Among our tests, we found that image-domain single-iteration LSRTM methods, including an extension of an approximate inverse Hessian method in the curvelet domain, not only compensated for amplitude loss due to poor illumination caused by complex salt bodies, but it also produced subsalt images with fewer migration artifacts in the field data. In contrast, an iterative inversion method showed its potential for broadening the bandwidth in the subsalt, but it was less effective in reducing migration artifacts and noise. Based on our understanding, we evaluated the current state of LSRTM for subsalt imaging.

Pure qP-wave least-squares reverse time migration in vertically transverse isotropic media and its application to field data

2020 ◽  
Vol 17 (2) ◽  
pp. 208-220
Author(s):  
Jian-Ping Huang ◽  
Xin-Ru Mu ◽  
Zhen-Chun Li ◽  
Qing-Yang Li ◽  
Shuang-Qi Yuan ◽  
...  
Keyword(s):  
Least Squares ◽  
Field Data ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Isotropic Media ◽  
Transverse Isotropic

scholarly journals Least-squares reverse time migration with local Radon-based preconditioning

Geophysics ◽  
2017 ◽  
Vol 82 (2) ◽  
pp. S75-S84 ◽  
Author(s):  
Gaurav Dutta ◽  
Matteo Giboli ◽  
Cyril Agut ◽  
Paul Williamson ◽  
Gerard T. Schuster
Keyword(s):  
Least Squares ◽  
Radon Transform ◽  
Field Data ◽  
Signal To Noise Ratio ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Inverse Mapping ◽  
Time Migration ◽  
Discrete Radon Transform ◽  
L2 Norm

Least-squares migration (LSM) can produce images with better balanced amplitudes and fewer artifacts than standard migration. The conventional objective function used for LSM minimizes the L2-norm of the data residual between the predicted and the observed data. However, for field-data applications in which the recorded data are noisy and undersampled, the conventional formulation of LSM fails to provide the desired uplift in the quality of the inverted image. We have developed a least-squares reverse time migration (LSRTM) method using local Radon-based preconditioning to overcome the low signal-to-noise ratio (S/N) problem of noisy or severely undersampled data. A high-resolution local Radon transform of the reflectivity is used, and sparseness constraints are imposed on the inverted reflectivity in the local Radon domain. The sparseness constraint is that the inverted reflectivity is sparse in the Radon domain and each location of the subsurface is represented by a limited number of geologic dips. The forward and the inverse mapping of the reflectivity to the local Radon domain and vice versa is done through 3D Fourier-based discrete Radon transform operators. The weights for the preconditioning are chosen to be varying locally based on the relative amplitudes of the local dips or assigned using quantile measures. Numerical tests on synthetic and field data validate the effectiveness of our approach in producing images with good S/N and fewer aliasing artifacts when compared with standard RTM or standard LSRTM.

Least-squares reverse time migration using controlled-order multiple reflections

Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. S347-S357 ◽  
Author(s):  
Yike Liu ◽  
Xuejian Liu ◽  
Are Osen ◽  
Yu Shao ◽  
Hao Hu ◽  
...  
Keyword(s):  
Least Squares ◽  
Seismic Data ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Multiple Reflections ◽  
Numerical Tests ◽  
Elimination Algorithm ◽  
Reflection Data ◽  
Time Migration ◽  
Least Squares Migration

Reverse time migration (RTM) using multiples generates inherent crosstalk artifacts due to the interference among multiples of different orders. We have developed a method to remove such crosstalk. This approach first separates the recorded seismic data into primary reflections and multiples using the surface-related multiples elimination algorithm and then isolates the multiples into different orders. We can take any specified, say the [Formula: see text]th, order of multiples data as the incident wave and the next higher order multiples data, ([Formula: see text])th order, as the corresponding primary reflection data for imaging. We have applied the least-squares migration scheme to these two successive orders of multiples. Our method is denoted as least-squares RTM using controlled-order multiples (LSRTM-CM). Our numerical tests demonstrated that LSRTM-CM can significantly improve imaging quality compared with straightforward seismic imaging using multiples without multiples separation.

A wavefield-separation-based elastic least-squares reverse time migration

Geophysics ◽  
2018 ◽  
Vol 83 (3) ◽  
pp. S279-S297 ◽  
Author(s):  
Bingluo Gu ◽  
Zhenchun Li ◽  
Jianguang Han
Keyword(s):  
Least Squares ◽  
Synthetic Data ◽  
P Wave ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
S Wave ◽  
S Waves ◽  
Numerical Tests ◽  
Time Migration

Elastic least-squares reverse time migration (ELSRTM) has the potential to provide improved subsurface reflectivity estimation. Compared with elastic RTM (ERTM), ELSRTM can produce images with higher spatial resolution, more balanced amplitudes, and fewer artifacts. However, the crosstalk between P- and S-waves can significantly degrade the imaging quality of ELSRTM. We have developed an ELSRTM method to suppress the crosstalk artifacts. This method includes three crucial points. The first is that the forward and backward wavefields are extrapolated based on the separated elastic velocity-stress equation of P- and S-waves. The second is that the separated vector P- and S-wave residuals are migrated to form reflectivity images of Lamé constants [Formula: see text] and [Formula: see text] independently. The third is that the reflectivity images of [Formula: see text] and [Formula: see text] are obtained by the vector P-wave wavefields achieved in the backward extrapolation of the separated vector P-wave residuals and the vector S-wave wavefields achieved in the backward extrapolation of the separated vector S-wave residuals, respectively. Numerical tests with synthetic data demonstrate that our ELSRTM method can produce images free of crosstalk artifacts. Compared with ELSRTM based on the coupled wavefields, our ELSRTM method has better convergence and higher accuracy.

From acoustic to elastic inverse extended Born modeling: a first insight in the marine environment

Geophysics ◽  
2021 ◽  
pp. 1-73
Author(s):  
Milad Farshad ◽  
Hervé Chauris
Keyword(s):  
Least Squares ◽  
Marine Environment ◽  
Computational Cost ◽  
Physical Parameters ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
S Wave ◽  
Time Migration ◽  
Physical Density ◽  
Acoustic Approximation

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.

Sign in / Sign up

Export Citation Format

Share Document