Abstract:The results of least-squares migration have higher fidelity and spatial resolution than that of conventional migration. In this paper, an efficient method for least-squares reverse time migration in imaging domain (IDLSRTM) is proposed by using spatial de-convolution. Given a smooth velocity model with accurate kinematics of seismic wave used for migration (only the migration based on scalar wave equation is considered), the conventional reverse time migration (RTM) can be viewed as one kind of rough solution to the linear inversion problem of subsurface reflectivity. In image domain, the conventional reverse time migration results are a spatial-variant integral of the least-squares reverse time migration results and Hessian matrix. Different column vectors in the Hessian matrix can be regarded as the point spread functions of different imaging points in imaging domain. Because of the smoothness of migration velocity model and the multifold of seismic data acquisition system, the point spread functions of imaging points in a certain local region have well similarity. Therefore, the conventional reverse time migration results can be approximately represented as a spatial-invariant convolution between the least-square reverse time migration results and the Hessian matrix in the local region, and then a least-square reverse time migration method in image domain based on the local region spatial-invariant de-convolution is proposed. The spatial-invariant de-convolution can be efficiently and stably realized by using the pseudo-generalized inverse approach to spatial-invariant de-convolution in wavenumber domain. The size of local region is a trade-off between the efficiency and the effect of the least-square reverse time migration. The method proposed in this paper is applied to synthetic and field seismic data, ideal results are obtained.

An efficient least-squares reverse time migration in image domain

Elastic Reverse Time Migration Without Wavefield Decomposition

Reflected wave least squares reverse time migration with angle illumination compensation

Isotropic Elastic Reverse Time Migration

Reverse-time Migration Based on Generalized Synthesis of Seismic Data in F Requency-Space Domain

3-D Image-Domain Least-Squares Reverse Time Migration with L1 Norm Constraint and Total Variation Regularization

Least squares migration of seismic data with reflection wave equation

Substituting Smoothing with Low-Rank Decomposition - Applications to Least-Squares Reverse Time Migration of Simultaneous Source and Incomplete Seismic Data

2-D and 3-D Q-Compensated Image-Domain Least-Squares Reverse Time Migration Through the Hybrid Point Spread Functions and the Hybrid Deblurring Filter

Source-independent Least-Squares Reverse Time Migration in Vertical Transversely Isotropic Media Based on the Student's T-Distribution

Least-squares Reverse Time Migration in the Presence of Velocity Errors

Non-Linear Least-Squares Reverse-Time Migration

2-D and 3-D Image-Domain Least-Squares Reverse Time Migration Through Point Spread Functions and Excitation-Amplitude Imaging Condition.

Least-squares Reverse Time Migration with Sparse Regularization in the 2D Wavelet Domain

Noniterative Least-Squares Reverse Time Migration Based on an Efficient Asymptotic Hessian/point-spread Function Estimation

Least-squares Reverse Time Migration for Blended Data with a Local Low-rank Constraint

An efficient reverse time migration method using local nearly analytic discrete operator

Least-Squares Reverse Time Migration for Reflection-Angle-Dependent Reflectivity.

Least-Squares Reverse Time Migration in the Presence of Density Variations

Mesh-free Least-Squares Reverse-Time Migration

3-D Q-Compensated Image-Domain Least-Squares Reverse Time Migration Through Point Spread Functions