Hussein Muhammed , AbdelHafiz Gadelmula , Zhenchun Li
{"title":"On the steadfastness of the least-squares reverse-time migration wavefield extrapolation via 1st-order Riemannian axis finite-difference solver","authors":"Hussein Muhammed , AbdelHafiz Gadelmula , Zhenchun Li","doi":"10.1016/j.rines.2025.100121","DOIUrl":null,"url":null,"abstract":"<div><div>Exploring Earth's deep regions via Least-Squares Reverse-Time Migration (LSRTM) methods is of significant interest due to its exceptional structural-amplitude clarity. This cutting-edge seismic imaging technique is time-consuming and memory-intensive, thus wavefield extrapolation was proposed to be in the Pseudodepth domain (1<sup>st</sup>-order Riemannian coordinate system’s axis) to address these issues and to prevent oversampling/aliasing when modeling deeper subsurface zones. Stabilizing the generated Riemannian wavefield involves implementing an appropriate mapping velocity and obtaining the vertical axis operator which partially converts the finite difference solver from time to frequency domains. Each Cartesian point <span><math><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow></math></span> has a corresponding vertical-time point <span><math><mrow><mo>(</mo><msub><mrow><mi>ξ</mi></mrow><mrow><mn>1</mn></mrow></msub><mo>,</mo><msub><mrow><mi>ξ</mi></mrow><mrow><mn>2</mn></mrow></msub><mo>,</mo><msub><mrow><mi>ξ</mi></mrow><mrow><mn>3</mn></mrow></msub><mo>)</mo></mrow></math></span>, allowing interpolation of the reconstructed source wavefield through a Cartesian-to-Riemannian mapping function. Our stability and convergence analysis indicates that the spatial derivatives of the 1<sup>st</sup>-order Riemannian axis can be approximated by Fourier pseudo-spectral methods and fast-Fourier transforms using a special Gaussian-like impulse function. This function generates the source term vector-matrix within the finite-difference operator. The mapping velocity, derived as a differential form of the initial input velocity model, controls the CFL conditions of the associated Riemannian-finite difference operator. Numerical, synthetic, and seismic field data examples show that this approach is more stable and efficient in extrapolating a smooth 1<sup>st</sup>-order Riemannian axis-based finite-difference wavefield while adhering to Claerbout’s principle for locating subsurface reflectors. Additionally, choosing the appropriate sampling rate for the new vertical axis is inversely related to the maximum frequency of the impulse wavelet and directly related to the minimum velocity value in the model.</div></div>","PeriodicalId":101084,"journal":{"name":"Results in Earth Sciences","volume":"3 ","pages":"Article 100121"},"PeriodicalIF":0.0000,"publicationDate":"2025-08-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Results in Earth Sciences","FirstCategoryId":"1085","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S2211714825000639","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0
Abstract
Exploring Earth's deep regions via Least-Squares Reverse-Time Migration (LSRTM) methods is of significant interest due to its exceptional structural-amplitude clarity. This cutting-edge seismic imaging technique is time-consuming and memory-intensive, thus wavefield extrapolation was proposed to be in the Pseudodepth domain (1st-order Riemannian coordinate system’s axis) to address these issues and to prevent oversampling/aliasing when modeling deeper subsurface zones. Stabilizing the generated Riemannian wavefield involves implementing an appropriate mapping velocity and obtaining the vertical axis operator which partially converts the finite difference solver from time to frequency domains. Each Cartesian point has a corresponding vertical-time point , allowing interpolation of the reconstructed source wavefield through a Cartesian-to-Riemannian mapping function. Our stability and convergence analysis indicates that the spatial derivatives of the 1st-order Riemannian axis can be approximated by Fourier pseudo-spectral methods and fast-Fourier transforms using a special Gaussian-like impulse function. This function generates the source term vector-matrix within the finite-difference operator. The mapping velocity, derived as a differential form of the initial input velocity model, controls the CFL conditions of the associated Riemannian-finite difference operator. Numerical, synthetic, and seismic field data examples show that this approach is more stable and efficient in extrapolating a smooth 1st-order Riemannian axis-based finite-difference wavefield while adhering to Claerbout’s principle for locating subsurface reflectors. Additionally, choosing the appropriate sampling rate for the new vertical axis is inversely related to the maximum frequency of the impulse wavelet and directly related to the minimum velocity value in the model.