Efficient least-squares reverse time migration with the approximate asymptotic inverse operator in place of the adjoint operator

Kai Yang,Jianfeng Zhang

DOI: https://doi.org/10.1190/geo2024-0002.1

IF: 3.264

2024-10-17

Geophysics

Abstract:GEOPHYSICS, Volume 89, Issue 6, Page S391-S404, November-December 2024.

geochemistry & geophysics

What problem does this paper attempt to address?

Least-squares Reverse Time Migration in Frequency Domain Using the Adjoint-State Method

Haoran Ren,Huazhong Wang,Shengchang Chen

DOI: https://doi.org/10.1088/1742-2132/10/3/035002

2013-01-01

Journal of Geophysics and Engineering

Abstract:A new scheme is presented to implement a least-squares frequency domain reverse time migration (LS-FRTM). This scheme expresses the gradient of the misfit function with respect to the model as the product of the conjugated Green's functions and the data residuals in the frequency domain based on the adjoint state method. In the 2D case, for each frequency all the Green's functions from the shots to the reflectors and from the reflectors to the receivers which depend on the background velocity can be calculated once using the lower/upper decomposition. The pseudo-Hessian matrix which is also expressed as a function of the Green's function is used as a substitute for the approximate Hessian to amplitude compensation for the gradient. Since the linearized inversion does not update the background velocity, the Green's function needs to be calculated only once. An iteration based LS-FRTM can be implemented with high efficiency. As examples supporting our assertion, we present the results obtained by applying our method to the 2D Marmousi model.
An efficient least-squares reverse time migration in image domain

Chen ShengChang,Li DaiGuang,Jin ChengMei

DOI: https://doi.org/10.6038/cjg2022P0532

2022-01-01

Chinese Journal of Geophysics

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.
Non-Linear Least-Squares Reverse-Time Migration

G. Yao,H. Jakubowicz

DOI: https://doi.org/10.3997/2214-4609.20130474

2012-01-01

Proceedings

Abstract:Migration attempts to produce an image of the subsurface by reversing the propagation effects in seismic data. Although in principle this requires the inverse of a modelling operator, in practice the adjoint of the modelling operator is used instead. In cases where the data are subject to significant aliasing, truncation, noise, or are incomplete, the adjoint of a modelling operator is not a good approximation to the inverse, and this degrades the resolution of the final migrated image. An improved approximation to the inverse operator can be obtained using a least-squares approach. In this work we present a least-squares formulation of reverse-time migration scheme that is based on an explicit matrix representation of generalised diffraction-stack migration. Our implementation uses a modified source wavelet to perform the forward and inverse steps at each iteration, and is regularised using a roughness penalty constraint. The results show increased resolution compared to conventional reverse-time migration.
Elastic Reverse-Time Migration Without Wavefield Decomposition

Daiguang Li,Sheng‐Chang Chen

DOI: https://doi.org/10.1190/image2022-3738445.1

2022-01-01

Abstract:PreviousNext No AccessSecond International Meeting for Applied Geoscience & EnergyElastic reverse-time migration without wavefield decompositionAuthors: Daiguang LiShengchang ChenDaiguang LiZheJiang UniversitySearch for more papers by this author and Shengchang ChenZheJiang UniversitySearch for more papers by this authorhttps://doi.org/10.1190/image2022-3738445.1 SectionsAboutPDF/ePub ToolsAdd to favoritesDownload CitationsTrack CitationsPermissions ShareFacebookTwitterLinked InRedditEmail AbstractElastic reverse time migration is a linear inversion of the reflectivity of the subsurface elastic interface using primarily reflected waves and primarily converted waves given a smooth migration model with accurate elastic wave kinematics. Combining this with the noncurl characteristic of P-waves and nondivergence characteristics of S-waves, we propose an innovative elastic reverse-time migration approach without wavefield decomposition applied to heterogeneous isotropic media. In this paper, pure propagators of P- and S-wave independent of each other are derived. We obtain ideal imaging results when we apply the proposed method to elastic synthetic seismic data generated by pure P- or S-wave source, demonstrating the validity of our method. Because the innovative migration method does not include wavefield decomposition which requires a large number of complex calculations, it has a great significant computational cost advantage over the migration method based on vector decomposition. Meanwhile, when compared to the Helmholtz separation migration method, it can avoid the loss of important geological information due to phase shift and amplitude distortion.Note: This paper was accepted into the Technical Program but was not presented at IMAGE 2022 in Houston, Texas.Keywords: elastic, reverse-time migration, linear inversion, pure propagatorsPermalink: https://doi.org/10.1190/image2022-3738445.1FiguresReferencesRelatedDetails Second International Meeting for Applied Geoscience & EnergyISSN (print):1052-3812 ISSN (online):1949-4645Copyright: 2022 Pages: 3694 publication data© 2022 Published in electronic format with permission by the Society of Exploration Geophysicists and the American Association of Petroleum GeologistsPublisher:Society of Exploration Geophysicists HistoryPublished Online: 15 Aug 2022 CITATION INFORMATION Daiguang Li and Shengchang Chen, (2022), "Elastic reverse-time migration without wavefield decomposition," SEG Technical Program Expanded Abstracts : 2778-2781. https://doi.org/10.1190/image2022-3738445.1 Plain-Language Summary Keywordselastic reverse-time migrationlinear inversionpure propagatorsPDF DownloadLoading ...
Isotropic Elastic Reverse Time Migration

Chen ShengChang,Li DaiGuang,Gong XinYue

DOI: https://doi.org/10.6038/cjg2022q0415

2024-01-01

Abstract:In order to obtain the migration imaging results with clear physical meaning and geological interpretation, the current elastic reverse time migration methods require P/S waves decomposition to the extrapolated elastic wavefields in the migration process to eliminate the cross-talk artifacts generated by physically mismatched wave modes in the migration imaging results. Given a smooth model of subsurface media with accurate kinematic characteristics of seismic waves, this paper first proposes a general forward representation for migration as a linear inverse problem of seismic data. Then, by using the isotropic elastic wave equation, the P/S wave decoupled elastic wave equations suitable for the smooth model of elastic media are derived, namely a P wave elastic wave equation and a S wave elastic wave equation. Finally, a new method and new realization for the isotropic elastic reverse time migration are proposed. For the reverse time migration of isotropic elastic wave seismic data excited by the pure P (or S) wave source, a new method is proposed, which does not need the P/S waves decomposition to the extrapolated elastic wavefields and does not appear cross-talk artifacts in the migration imaging results. For the reverse time migration of isotropic elastic wave seismic data excited by P/S waves composite source, a new realization is put forward that even if the P/S waves decomposition is used to the extrapolated elastic wavefields, the cross-talk artifacts still appear in the migration imaging results.
Elastic Reverse Time Migration Without Wavefield Decomposition

Daiguang Li,Shengchang Chen

DOI: https://doi.org/10.1016/j.jappgeo.2022.104818

IF: 1.845

2022-01-01

Journal of Applied Geophysics

Abstract:In this paper, we provide a clear definition represented by a concise expression for seismic data migration that is a linear inversion of subsurface exploration targets using primary waves given a smooth migration model with accurate wave kinematics. Combining the definition with noncurl characteristic of P-wave and nondivergence characteristic of S-wave, we propose an elastic reverse-time migration approach without wavefield decomposition applied to heterogeneous isotropic media, in which decoupled P- and S-wave propagators independent of each other are derived. We obtain ideal imaging results when we apply the proposed migration approach to elastic synthetic seismic data of a simple four-layered model and the complex Marmousi2 model with pure P- or S-wave source exciting, demonstrating the validity of our method. Because the proposed method does not include wavefield decomposition in the backward extrapolation of recorded wavefield and demand less spatial derivative numbers in the forward extrapolation of source wavefield, it is capable of achieving high-efficiency computing. In the case of P/S composite source exciting, we suggest that although the recorded wavefield can be divided into P- and S-wave component, the P-wave component consists of primarily reflected waves PP plus primarily converted waves SP and the S-wave component consists of primarily reflected waves SS plus primarily converted waves PS. Because PP and SP or SS and PS cannot be further separated from each other by current decoupling methods, the crosstalk artifacts are unavoidable, which is verified by the migration imaging results of a two-layered model.
Reflected wave least squares reverse time migration with angle illumination compensation

Huamin Zhou,Shengchang Chen,Liming Zhou,Haoran Ren,Rushan Wu,Guoqiang Xiao

DOI: https://doi.org/10.1007/s11600-019-00368-1

IF: 2.293

2019-01-01

Acta Geophysica

Abstract:High-quality seismic data imaging plays an important role in the lithological interpretation of subsurface structures. However, high-quality imaging remains a challenging task. Based on the linear inversion theory of reflected wave equations, this paper proposes reflected wave least squares reverse time migration with angle illumination compensation to better balance the amplitude of seismic imaging. We use the reflected wave migration equation to unify forward and backward propagation, which helps to obtain an image with correct phase and symmetric waveform. Under the assumption that the spectrum of seismic wavefield remains unchanged, the Poynting vector method is used to efficiently calculate the propagation direction of seismic waveform and seismic illumination in the angle domain. During iteration, angle-domain illumination is used as a preconditioner to compensate for the amplitude of the iterated gradient terms based on the angle value. In this manner, we can enhance the imaging energy of steeply inclined structures. To improve the stability of linear inversion, the spatial derivative of the image is used as a regularized constraint term. Numerical tests show that the proposed method can suppress imaging noise as well as improve resolution and amplitude fidelity of the images. Furthermore, the inversed result can be used to estimate underground reflectivity, which is important for the further development of seismic inversion technology.
An efficient reverse time migration method using local nearly analytic discrete operator

JingShuang Li,Dinghui Yang,Faqi Liu

DOI: https://doi.org/10.1190/geo2012-0247.1

IF: 3.264

2013-01-01

Geophysics

Abstract:In recent years, reverse time migration (RTM), the most powerful depth imaging method, has become the preferred imaging tool in many geologic settings because of its ability to handle complex velocity models including steeply dipping interfaces. Finite difference is one of the most popular numerical methods applied in RTM in the industry. However, it often encounters a serious issue of numerical dispersion, which is typically suppressed by reducing the propagation grid sizes, resulting in large computation and memory increment. Recently, a nearly analytic discrete operator has been developed to approximate the partial differential operators, from which many antidispersion schemes have been proposed, and are confirmed to be superior to conventional algorithms in suppressing numerical dispersion. We apply an optimal nearly analytic discrete (ONAD) method to RTM to improve its accuracy and performance. Numerical results show that ONAD can be used effectively in seismic modeling and migration based on the full wave equations. This method produces little numerical dispersion and requires much less computation and memory compared to the traditional finite-difference methods such as Lax-Wendroff correction method. The reverse time migration results of the 2D Marmousi model and the Sigsbee2B data set show that ONAD can improve the computational efficiency and maintain image quality by using large extrapolation grids.
A new correlation based least-squares reverse time migration

Feipeng Li,Jinghuai Gao,Xiudi Jiang,Wenbo Sun

DOI: https://doi.org/10.1190/segam2018-2995223.1

2018-01-01

Abstract:PreviousNext No AccessSEG Technical Program Expanded Abstracts 2018A new correlation based least-squares reverse time migrationAuthors: Feipeng LiJinghuai GaoXiudi JiangWenbo SunFeipeng LiSchool of Electronic and Information Engineering, Xi’an Jiaotong University, ChinaSearch for more papers by this author, Jinghuai GaoSchool of Electronic and Information Engineering, Xi’an Jiaotong University, CHINANational Engineering Laboratory for Offshore Oil Exploration, Xi’an Jiaotong University, ChinaSearch for more papers by this author, Xiudi JiangCNOOC Research Institue, ChinaSearch for more papers by this author, and Wenbo SunCNOOC Research Institue, ChinaSearch for more papers by this authorhttps://doi.org/10.1190/segam2018-2995223.1 SectionsSupplemental MaterialAboutPDF/ePub ToolsAdd to favoritesDownload CitationsTrack CitationsPermissions ShareFacebookTwitterLinked InRedditEmail Abstract Least-squares reverse time migration (LSRTM) applies reverse time migration and demigration as migration and modeling operator to minimize the difference between the simulated and observed seismic data. We propose a new scheme for LSRTM and derive its corresponding conjugate gradient method to seek for the optimal image. By applying a correlation based objective function to LSRTM, we are able to emphasize not only on amplitude matching but also on phase information. The objective function incorporates the auto-correlation besides the cross-correlation. And it is the least-square of the two correlations. Numerical tests on synthetic models demonstrate that our new LSRTM scheme improves true amplitude imaging in the case of severely undersampled data and simultaneous source data. Presentation Date: Monday, October 15, 2018 Start Time: 1:50:00 PM Location: 207A (Anaheim Convention Center) Presentation Type: Oral Keywords: least squares, reverse time migration, inversionPermalink: https://doi.org/10.1190/segam2018-2995223.1FiguresReferencesRelatedDetails SEG Technical Program Expanded Abstracts 2018ISSN (print):1052-3812 ISSN (online):1949-4645Copyright: 2018 Pages: 5520 publication data© 2018 Published in electronic format with permission by the Society of Exploration GeophysicistsPublisher:Society of Exploration Geophysicists HistoryPublished Online: 27 Aug 2018 CITATION INFORMATION Feipeng Li, Jinghuai Gao, Xiudi Jiang, and Wenbo Sun, (2018), "A new correlation based least-squares reverse time migration," SEG Technical Program Expanded Abstracts : 4236-4240. https://doi.org/10.1190/segam2018-2995223.1 Plain-Language Summary Keywordsleast squaresreverse time migrationinversionPDF DownloadLoading ...
Frequency Domain Adjoint-state Least-squares Migration

h z wang,s c chen,h r ren

DOI: https://doi.org/10.3997/2214-4609.20130241

2013-01-01

Abstract:Least-Squares Migration based on the inversion formulation will provide the amplitude-preserved imaging if inversed-Hessian kernel are applied. The inversed-Hessian kernel reflects the forward and backward propagation of the operator. Under the assuming of infinite receiver coverage and smooth velocity model, LSM using one-way operator can provide believable imaging. However, in other cases, for the areas the operator can not covered, LSM can not compensate. To avoid the zero-illuminated cases, we apply the full-wave operator to propagate the wave in frequency domain. Frequency-domain Reverse Time Migration use efficient full wave extrapolator in the frequency domain. We combine the LSM frame and the FRTM method and implement a Least-Squares Frequency domain Reverse Time Migration. To enhance the computational efficiency, we introduce the Adjoint State method based on the local Born approximation. The algorithms are tested using the BP 2004 Benchmark model and results show improved image quality and amplitudes.
Noniterative Least-Squares Reverse Time Migration Based on an Efficient Asymptotic Hessian/point-spread Function Estimation

Peng Xu,Bo Feng,Huazhong Wang,Shaoyong Liu

DOI: https://doi.org/10.1190/geo2021-0401.1

IF: 3.264

2022-01-01

Geophysics

Abstract:Although reverse time migration (RTM) is a powerful imaging method, it uses only the adjoint of the wave-propagation operator and produces only a blurred image. Through linearized inversion, least-squares RTM (LSRTM) reduces migration artifacts, balances the amplitudes, and improves the resolution of the imaging results. Compared to data-domain LSRTM, image-domain LSRTM (ID-LSRTM) is a more practical choice because it typically requires one round of RTM and subsequently solves the deblurring problem of the conventional RTM image involving the Hessian matrix. The Hessian matrix can be effectively replaced by a set of point-spread functions (PSFs). PSFs are traditionally obtained through one round of demigration-migration; hence, the traditional PSF estimation has a computational cost that is twice that of the RTM. We have deduced an analytic expression of the Hessian/PSF by using an asymptotic Green’s function and derived an efficient estimation strategy for the Hessian/PSF. The proposed method significantly reduces the computational cost and allows more flexible spatial sampling of PSFs without interference. Based on the explicitly calculated PSFs, ID-LSRTM can be achieved by solving an optimization problem with combined total-variation and curvelet-domain [Formula: see text] regularizations. The total computational cost of the proposed LSRTM is approximately the same as that of the conventional RTM. Numerical examples of synthetic and field data find that the proposed LSRTM significantly improves image quality while being highly efficient.
Correlative Least-Squares Reverse Time Migration in Viscoelastic Media

Wei Zhang,Jinghuai Gao,Feipeng Li,Ying Shi,Xuan Ke

DOI: https://doi.org/10.1016/j.jappgeo.2021.104256

IF: 1.845

2021-01-01

Journal of Applied Geophysics

Abstract:Standard elastic least-squares reverse time migration (LSRTM) generally ignores the subsurface attenuation effects which can lead to blurred and unfocused migration images. To correct these attenuation effects, we develop a correlative Q-LSRTM approach in viscoelastic media. We use the viscoelastic wave equation based on the standard linear solid model to simulate the subsurface attenuation. The Born modeling and adjoint state equations in viscoelastic media can be derived by the Born approximation and adjoint state method, respectively. Instead of the conventional amplitude-matching misfit function, the proposed Q-LSRTM approach utilizes the normalized zero-lag crosscorrelation misfit function to emphasize the phase measure between the predicted and observed data. The gradients of the crosscorrelation misfit function with regard to reflectivity images of P- and S-wave impedance are estimated by a zero-lag crosscorrelation between the forward background wavefields and the reversal adjoint wavefields. In view of computational efficiency, we also derive an analytical step-length formula for the correlative Q-LSRTM approach, which only requires to solve the Born modeling operator once. Numerical examples using several 2D models have determined that the proposed Q-LSRTM can improve better image quality, compared with the standard elastic LSRTM approach. Furthermore, the crosscorrelation misfit function is less sensitive to amplitude errors.
Elastic reverse time migration using an efficient and accurate adaptive variable grid discretization method

Xiujin Yang,Jianping Huang,Yi Shen,Liang Chen

DOI: https://doi.org/10.1093/jge/gxae085

2024-08-26

Journal of Geophysics and Engineering

Abstract:The conventional reverse time migration utilizes regularly sampled computational grids to simulate wave propagation. Selecting the appropriate grid sampling is important for computational accuracy and efficiency. In general, the uniform-size grid cannot well represent the complexity of the geology. The grid may appear sparse in the low-velocity zone, especially in shallow depths where dispersion may occur. Conversely, it may appear excessively dense in the high-velocity zone, such as at greater depths or within a salt body, which results in higher computational memory and time consumption. To overcome these issues, we developed an efficient and accurate adaptive variable grid discretization method that automatically selects the vertical grid size based on the velocity, depth, and dominant frequency of the wavelet in elastic medium. Then we reformulated the elastic equations based on the adaptive variable grid by introducing a mapping relationship. To test the effectiveness, accuracy, and efficiency of the equation, we implemented it to both the forward propagation and migration of elastic wavefield. Synthetic numerical examples demonstrate that our proposed method can achieve elastic wavefield separation and no significant dispersion phenomenon. The multi-component imaging accuracy of reverse time migration is nearly equivalent to the traditional method, while significantly improving computational efficiency and saving storage space.

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

Zilong Ye,Jianping Huang,Yi Shen,Xinru Mu,Qiwei Zhan

DOI: https://doi.org/10.1093/jge/gxad054

2023-01-01

Journal of Geophysics and Engineering

Abstract:Seismic anisotropy exists in various types of strata and should be considered in seismic imaging schemes. Seismic imaging algorithms based on isotropic assumption neglect the effects of anisotropy on seismic data, which cause migration artifacts and waveform distortion. To correct the effects of anisotropy on seismic wave propagation, we propose an imaging algorithm that performs least-squares reverse time migration in vertical transversely isotropic acoustic media. We derive the following operators to implement this algorithm, the de-migration operator, its adjoint migration operator, and the corresponding gradient. However, an inaccurate estimated source wavelet will introduce error in the seismic simulation, and thus increase the mismatch between observed and synthetic data for least-squares reverse time migration. In addition, the noises, especially the noises with abnormal amplitudes in the seismic data, damage the inversion convergence and reduce the imaging resolution. To improve the image quality, we propose to use convolved wavefields between observed and synthetic data so that such a mismatch can be independent of the source wavelets. Also, we use the Student's t-distribution instead of L2 norm in our inversion scheme to better handle the seismic noise. Its implementation only modifies the gradient of the conventional least-square reverse time migration scheme. Our numerical tests show a clear improvement using our proposed imaging algorithm when compared with the conventional isotropic migration scheme for the anisotropic data. Also, the synthetic examples demonstrate the feasibility and effectiveness of our proposed source-independent algorithm using the Student's t-distribution.
Decoupled elastic least-squares reverse time migration and its application in tunnel geologic forward prospecting

Bin Liu,Jiansen Wang,Yuxiao Ren,Xu Guo,Lei Chen,Lanbo Liu

DOI: https://doi.org/10.1190/geo2020-0875.1

IF: 3.264

2021-12-21

Geophysics

Abstract:Accurate seismic imaging can ensure safe and efficient tunnel construction under complex geologic conditions. As a high-precision migration method, reverse time migration (RTM) has been introduced into tunnel seismic forward prospecting. However, the resolution of traditional RTM imaging results may not meet the requirements in a complex tunnel environment, which affects the interpretation of tunnel seismic forward-prospecting results. We have developed a least-squares RTM method based on the decoupled elastic wave equation in tunnels. The Born forward-modeling operator and its exact adjoint migration imaging operator are derived to ensure a stable convergence of the conjugate-gradient method. Moreover, a pseudo-Hessian-based preconditioning operator is adopted to accelerate the convergence. Numerical examples are provided to verify the efficiency of our scheme. A field test in a traffic tunnel construction site is performed to show the good application effect of the decoupled elastic least-squares RTM in practical situations.

geochemistry & geophysics
Sparse least-squares reverse time migration using 2-D undecimated wavelet transform

Feipeng Li,Jinghuai Gao

DOI: https://doi.org/10.1190/igc2018-151

2018-01-01

Abstract:PreviousNext No AccessInternational Geophysical Conference, Beijing, China, 24-27 April 2018Sparse least-squares reverse time migration using 2-D undecimated wavelet transformAuthors: Feipeng LiJinghuai GaoFeipeng LiSearch for more papers by this author and Jinghuai GaoXi'an Jiaotong UniversitySearch for more papers by this authorhttps://doi.org/10.1190/IGC2018-151 SectionsAboutPDF/ePub ToolsAdd to favoritesDownload CitationsTrack CitationsPermissions ShareFacebookTwitterLinked InRedditEmail Abstract We propose a new scheme for least-squares reverse time migration (LSRTM) under the constraint of sparseness in 2-D undecimated wavelet domain. Taking 2-D undecimated wavelet as an effective representation for the reflectivity model, we incorporate wavelet domain sparseness into LSRTM as a regularization term. To achieve sparseness, we use an iteratively reweighted least-squares algorithm (IRLS) whose weighting matrix is constructed via solution itself. Then we derive a wavelet domain conjugate gradient method to seek for the optimal image. Numerical tests show that this new scheme is capable of suppressing migration artifacts. Keywords: least-squares, reverse time migration, algorithmPermalink: https://doi.org/10.1190/IGC2018-151FiguresReferencesRelatedDetailsCited byDeep convolutional neural network and sparse least-squares migrationZhaolun Liu, Yuqing Chen, and Gerard Schuster13 June 2020 | GEOPHYSICS, Vol. 85, No. 4 International Geophysical Conference, Beijing, China, 24-27 April 2018ISSN (online):2159-6832Copyright: 2018 Pages: 1821 publication data© 2018 Published in electronic format with permission by the Society of Exploration Geophysicists and Chinese Geophysical SocietyPublisher:Society of Exploration Geophysicists HistoryPublished Online: 11 Dec 2018 CITATION INFORMATION Feipeng Li and Jinghuai Gao, (2018), "Sparse least-squares reverse time migration using 2-D undecimated wavelet transform," SEG Global Meeting Abstracts : 616-619. https://doi.org/10.1190/IGC2018-151 Plain-Language Summary Keywordsleast-squaresreverse time migrationalgorithmPDF DownloadLoading ...
Elastic least‐squares reverse time migration from topography through anisotropic tensorial elastodynamics

Tugrul Konuk,Jeffrey Shragge

DOI: https://doi.org/10.1111/1365-2478.13619

2024-10-10

Geophysical Prospecting

Abstract:Least‐squares reverse time migration is an increasingly popular technique for subsurface imaging, especially in the presence of complex geological structures. However, elastic least‐squares reverse time migration algorithms face significant practical and numerical challenges when migrating multi‐component seismic data acquired from irregular topography. Many associated issues can be avoided by abandoning the Cartesian coordinate system and migrating the data to a generalized topographic coordinate system conformal to surface topology. We introduce a generalized anisotropic elastic least‐squares reverse time migration methodology that uses the numerical solutions of tensorial elastodynamics for propagating wavefields in computational domains influenced by free‐surface topography. We define a coordinate mapping assuming unstretched vertically translated meshes that transform an irregular physical domain to a regular computational domain on which calculating numerical elastodynamics solutions is straightforward. This allows us to obtain numerical solutions of forward and adjoint elastodynamics and generate subsurface images directly in topographic coordinates using a tensorial energy‐norm imaging condition. Numerical examples demonstrate that the proposed generalized elastic least‐squares reverse time migration algorithm is suitable for generating high‐quality images with reduced artefacts and better balanced reflectivity that can accurately explain observed data acquired from topography in a medium characterized by arbitrary heterogeneity and anisotropy. Finally, the computational cost of our method is comparable to that of an equivalent Cartesian elastic least‐squares reverse time migration numerical implementation.

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

Yangkang Chen,Min Bai,Yatong Zhou,Qingchen Zhang,Yufeng Wang,Hanming Chen

DOI: https://doi.org/10.1190/geo2017-0298.1

IF: 3.264

2019-01-01

Geophysics

Abstract:Seismic migration can be formulated as an inverse problem, the model of which can be iteratively inverted via the least-squares migration framework instead of approximated by applying the adjoint operator to the observed data. Least-squares reverse time migration (LSRTM) has attracted more and more attention in modern seismic imaging workflows because of its exceptional performance in obtaining high-resolution true-amplitude seismic images and the fast development of the computational capability of modern computing architecture. However, due to a variety of reasons, e.g., insufficient shot coverage and data sampling, the image from least-squares inversion still contains a large amount of artifacts. This phenomenon results from the ill-posed nature of the inverse problem. In traditional LSRTM, the minimum least-squares energy of the model is used as a constraint to regularize the inverse problem. Considering the residual noise caused by the smoothing operator in traditional LSRTM, we regularize the model using a powerful low-rank decomposition operator, which can better suppress the migration artifacts in the image during iterative inversion. We evaluate in detail the low-rank decomposition operator and the way to apply it along the geologic structure of seismic reflectors. We comprehensively analyze the performance of our algorithm in attenuating crosstalk noise caused by simultaneous source acquisition and migration artifacts caused by insufficient space sampling via two synthetic examples and one field data example. Our results indicate that compared to the conventional smoothing operator, our low-rank decomposition operator can help obtain a cleaner LSRTM image and obtain a slightly better edge-preserving performance.
Least-Squares Reverse Time Migration for Reflection-Angle-Dependent Reflectivity.

Chuang Li,Yikang Zhang,Zhaoqi Gao,Feipeng Li,Zhen Li,Jinghuai Gao

DOI: https://doi.org/10.1109/tgrs.2022.3218113

IF: 8.2

2022-01-01

IEEE Transactions on Geoscience and Remote Sensing

Abstract:Least-squares reverse time migration (LSRTM) can estimate high-quality reflectivity of subsurface medium from seismic data. However, the subsurface reflectivity depends on reflection angles, and its variations over reflection angles are extremely important because they can be used to estimate subsurface physical properties for seismic interpretation. We present a new formulation of the LSRTM method that can estimate reflection-angle-dependent reflectivity from seismic data. We derive a forward modeling operator which predicts the reflection data without calculating the reflection angles, and verify that it approximately equals to the reflection-angle-dependent wave-equation-based Kirchhoff modeling operator under the assumption that the velocity perturbation is small and the reflection angle is smaller than the critical angle. Based on the proposed modeling operator associated with the adjoint of the angle-dependent wave-equation-based Kirchhoff modeling operator, we reformulate LSRTM as an inverse problem to invert for reflection-angle-dependent reflectivity using a preconditioned conjugate gradient algorithm. The algorithm uses a low-rank filter as the preconditioner to attenuate migration artifacts. Imaging tests on synthetic and field seismic data are used to verify validity and superiority of the proposed method. The tests illustrate that the proposed method can produce the reflection-angle-dependent reflectivity with much higher signal-to-noise ratio (SNR), resolution and amplitude fidelity than reverse time migration. Compared with conventional LSRTM, it can produce more focused stacked image when the migration velocity contains errors. Moreover, conventional LSRTM only produces the angle-independent reflectivity, whereas the proposed method has the feasibility to produce the reflection-angle-dependent reflectivity.
Least-squares Reverse Time Migration for Blended Data with a Local Low-rank Constraint

K. Xiang,X. Chen,H. Chen,Y. Chen

DOI: https://doi.org/10.3997/2214-4609.201601200

2016-01-01

Abstract:Summary The simultaneous-source shooting technique can accelerate field acquisition and improve spatial sampling but will cause strong interferences in the recorded data. Direct imaging of blended simultaneous-source data has been demonstrated to be an promising research field since there can be no need for the separation for blended sources before the subsequent traditional processing and imaging. The key issue in direct imaging of blended data is the strong artifacts in the migrated image. Although the least-squares migration can help reduce some artifacts, there are still residual artifacts in the image. Those artifacts mainly appears on the shallow part of the image and appears as spatially incoherent noise. We propose to apply the singular spectrum analysis (SSA) operator to attenuate the such artifacts during least-squares inversion. Considering that global SSA cannot deal with over-complicated data well, we propose to use local SSA in order to better remove noise and preserve steeply dipping components. The local SSA operator corresponds to a local low-rank constraint applied in the inversion process. The migration operator used in the study is the reverse time migration (RTM) operator. We use the Marmousi model example to show the superior performance of the proposed algorithm.

Efficient least-squares reverse time migration with the approximate asymptotic inverse operator in place of the adjoint operator

Least-squares Reverse Time Migration in Frequency Domain Using the Adjoint-State Method

An efficient least-squares reverse time migration in image domain

Non-Linear Least-Squares Reverse-Time Migration

Elastic Reverse-Time Migration Without Wavefield Decomposition

Isotropic Elastic Reverse Time Migration

Elastic Reverse Time Migration Without Wavefield Decomposition

Reflected wave least squares reverse time migration with angle illumination compensation

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

A new correlation based least-squares reverse time migration

Frequency Domain Adjoint-state Least-squares Migration

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

Correlative Least-Squares Reverse Time Migration in Viscoelastic Media

Elastic reverse time migration using an efficient and accurate adaptive variable grid discretization method

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

Decoupled elastic least-squares reverse time migration and its application in tunnel geologic forward prospecting

Sparse least-squares reverse time migration using 2-D undecimated wavelet transform

Elastic least‐squares reverse time migration from topography through anisotropic tensorial elastodynamics

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

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

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