Elastic least-squares reverse time migration

Geophysics ◽  
2017 ◽  
Vol 82 (4) ◽  
pp. S315-S325 ◽  
Author(s):  
Yuting Duan ◽  
Antoine Guitton ◽  
Paul Sava
Keyword(s):  
Least Squares ◽  
Reservoir Characterization ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Imaging Condition ◽  
Time Migration ◽  
Energy Focusing ◽  
Polarity Reversals ◽  
New Perturbation ◽  
Least Squares Migration

Least-squares migration can produce images with improved resolution and reduced migration artifacts, compared with conventional imaging. We have developed a method for elastic least-squares reverse time migration (LSRTM) based on a new perturbation imaging condition that yields scalar images of squared P- and S-velocity perturbations. These perturbation images do not suffer from polarity reversals that are common for more conventional elastic imaging methods. We use 2D synthetic and field-data examples to demonstrate the proposed LSRTM algorithm using the perturbation imaging condition. Our results show that elastic LSRTM improves the energy focusing and illumination of the elastic images and it attenuates artifacts resulting, for instance, from sparseness in the wavefield sampling and crosstalk of the P- and S-modes. Compared with RTM images, the LSRTM images provide more accurate relative amplitude information that is useful for reservoir characterization.

An exact adjoint operation pair in time extrapolation and its application in least-squares reverse-time migration

Geophysics ◽  
2009 ◽  
Vol 74 (5) ◽  
pp. H27-H33 ◽  
Author(s):  
Jun Ji
Keyword(s):  
Least Squares ◽  
Adjoint Operator ◽  
Synthetic Data ◽  
Forward Modeling ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Product Test ◽  
Least Squares Migration ◽  
Time Extrapolation

To reduce the migration artifacts arising from incomplete data or inaccurate operators instead of migrating data with the adjoint of the forward-modeling operator, a least-squares migration often is considered. Least-squares migration requires a forward-modeling operator and its adjoint. In a derivation of the mathematically correct adjoint operator to a given forward-time-extrapolation modeling operator, the exact adjoint of the derived operator is obtained by formulating an explicit matrix equation for the forward operation and transposing it. The programs that implement the exact adjoint operator pair are verified by the dot-product test. The derived exact adjoint operator turns out to differ from the conventional reverse-time-migration (RTM) operator, an implementation of wavefield extrapolation backward in time. Examples with synthetic data show that migration using the exact adjoint operator gives similar results for a conventional RTM operator and that least-squares RTM is quite successful in reducing most migration artifacts. The least-squares solution using the exact adjoint pair produces a model that fits the data better than one using a conventional RTM operator pair.

A new imaging condition for least-squares reverse time migration with local similarity

2018 ◽  
Author(s):  
Young Seo Kim ◽  
Ali Almomin ◽  
Woodon Jeong ◽  
Constantine Tsingas
Keyword(s):  
Least Squares ◽  
Local Similarity ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Imaging Condition ◽  
Time Migration

scholarly journals An adaptive unstructured mesh based solution to topography least-squares reverse-time imaging

2019 ◽  
Author(s):  
Qiancheng Liu ◽  
Jianfeng Zhang
Keyword(s):  
Least Squares ◽  
Unstructured Mesh ◽  
Voronoi Tessellation ◽  
Grid Method ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Imaging Condition ◽  
Time Migration ◽  
Rugged Topography ◽  
Backward Propagation

Abstract. Least-squares reverse-time migration (LSRTM) attempts to invert for the broadband-wavenumber reflectivity image by minimizing the residual between observed and predicted seismograms via linearized inversion. However, rugged topography poses a challenge in front of LSRTM. To tackle this issue, we present an unstructured mesh-based solution to topography LSRTM. As to the forward/adjoint modeling operators in LSRTM, we take a so-called unstructured mesh-based “grid method”. Before solving the two-way wave equation with the grid method, we prepare for it a velocity-adaptive unstructured mesh using a Delaunay Triangulation plus Centroidal Voronoi Tessellation (DT-CVT) algorithm. The rugged topography acts as constraint boundaries during mesh generation. Then, by using the adjoint method, we put the observed seismograms to the receivers on the topography for backward propagation to produce the gradient through the cross-correlation imaging condition. We seek the inverted image using the conjugate gradient method during linearized inversion to linearly reduce the data misfit function. Through the 2D SEG Foothill synthetic dataset, we see that our method can handle the LSRTM from rugged topography.

Least-squares migration with primary- and multiple-guided weighting matrices

Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. S171-S185 ◽  
Author(s):  
Chuang Li ◽  
Jianping Huang ◽  
Zhenchun Li ◽  
Han Yu ◽  
Rongrong Wang
Keyword(s):  
Free Surface ◽  
Least Squares ◽  
Seismic Data ◽  
Synthetic Data ◽  
Velocity Model ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Multiple Reflections ◽  
Time Migration ◽  
Least Squares Migration

Least-squares migration (LSM) of seismic data is supposed to produce images of subsurface structures with better quality than standard migration if we have an accurate migration velocity model. However, LSM suffers from data mismatch problems and migration artifacts when noise pollutes the recorded profiles. This study has developed a reweighted least-squares reverse time migration (RWLSRTM) method to overcome the problems caused by such noise. We first verify that spiky noise and free-surface multiples lead to the mismatch problems and should be eliminated from the data residual. The primary- and multiple-guided weighting matrices are then derived for RWLSRTM to reduce the noise in the data residual. The weighting matrices impose constraints on the data residual such that spiky noise and free-surface multiple reflections are reduced whereas primary reflections are preserved. The weights for spiky noise and multiple reflections are controlled by a dynamic threshold parameter decreasing with iterations for better results. Finally, we use an iteratively reweighted least-squares algorithm to minimize the weighted data residual. We conduct numerical tests using the synthetic data and compared the results of this method with the results of standard LSRTM. The results suggest that RWLSRTM is more robust than standard LSRTM when the seismic data contain spiky noise and multiple reflections. Moreover, our method not only suppresses the migration artifacts, but it also accelerates the convergence.

Least-squares migration of multisource data with a deblurring filter

Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. R135-R146 ◽  
Author(s):  
Wei Dai ◽  
Xin Wang ◽  
Gerard T. Schuster
Keyword(s):  
Wave Equation ◽  
Least Squares ◽  
Computational Efficiency ◽  
Computational Cost ◽  
Processing Technique ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Migration Method ◽  
Least Squares Migration

Least-squares migration (LSM) has been shown to be able to produce high-quality migration images, but its computational cost is considered to be too high for practical imaging. We have developed a multisource least-squares migration algorithm (MLSM) to increase the computational efficiency by using the blended sources processing technique. To expedite convergence, a multisource deblurring filter is used as a preconditioner to reduce the data residual. This MLSM algorithm is applicable with Kirchhoff migration, wave-equation migration, or reverse time migration, and the gain in computational efficiency depends on the choice of migration method. Numerical results with Kirchhoff LSM on the 2D SEG/EAGE salt model show that an accurate image is obtained by migrating a supergather of 320 phase-encoded shots. When the encoding functions are the same for every iteration, the input/output cost of MLSM is reduced by 320 times. Empirical results show that the crosstalk noise introduced by blended sources is more effectively reduced when the encoding functions are changed at every iteration. The analysis of signal-to-noise ratio (S/N) suggests that not too many iterations are needed to enhance the S/N to an acceptable level. Therefore, when implemented with wave-equation migration or reverse time migration methods, the MLSM algorithm can be more efficient than the conventional migration method.

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.

Sign in / Sign up

Export Citation Format

Share Document