Dual Reverse Time Migration: Efficient Least-squares Imaging with Multiple Scattering in Dual Space
Abstract
Least-squares reverse time migration (LSRTM) estimates a broadband reflectivity image by iteratively minimizing the data misfit. Despite its superior imaging performance over conventional reverse time migration (RTM), LSRTM is fundamentally limited by the first-order Born approximation and incurs a high computational cost due to the iterative solution of a source-dependent model-space Hessian. We propose dual RTM, a seismic imaging framework that reformulates the inverse problem in the dual space, where the optimization variable is the scattering information of all orders, rather than being limited to the first-order scattering of the Born approximation. Unlike LSRTM, which directly estimates the reflectivity model, dual RTM reconstructs the missing scattering data (dual variable) through a fixed-point iteration governed by the full wave equation, enabling the recovery of higher-order scattering and multiple reflections without repeated background velocity updates. Furthermore, dual RTM replaces the source-dependent model-space Hessian with a source-independent data-space Hessian. A hybrid time--frequency formulation enables this Hessian to be assembled and inverted once in the Fourier domain as an offline preprocessing step, reducing the online computational cost of each dual iteration to that of a single RTM application. Numerical experiments on challenging synthetic benchmarks demonstrate that dual RTM suppresses acquisition footprints, improves amplitude fidelity, and enhances deep imaging where primary reflections alone provide insufficient illumination. These results establish dual RTM as an efficient bridge between conventional RTM/LSRTM and nonlinear waveform inversion for high-resolution large-scale seismic imaging.
Related
- Hydroxy artificial water channels enable ultraselective reverse osmosis membranes for N-nitrosamine removal in water reuse
- Assessment of fate of contaminants of emerging concern in forward osmosis-reverse osmosis-anaerobic digestion hybrid: A literature-based approach
- Polyamide Reverse Osmosis Membranes Are Less Crosslinked Than Previously Thought
- Changing rounds into squares or combining stripes? In-depth analysis of Fritillaria flowers and the color-changing pattern of Sarcophaga flies sheds light on the diversity and morphogenesis of checkerboard patterns in Eukaryotes
- Elevated fluid pressure in compression facilitates flexural reverse faulting
BibTeX
@misc{sonbolestan2026dual,
title = {{Dual Reverse Time Migration: Efficient Least-squares Imaging with Multiple Scattering in Dual Space}},
author = {Mahdi Sonbolestan and Toktam Zand and Ali Gholami},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.25995},
eprint = {2610.25995},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.25995},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Dual Reverse Time Migration: Efficient Least-squares Imaging with Multiple Scattering in Dual Space
AU - Mahdi Sonbolestan
AU - Toktam Zand
AU - Ali Gholami
PY - 2026
DA - 2026/10/19
PB - arXiv
DO - 10.48550/arxiv.2610.25995
UR - https://arxiv.org/abs/2610.25995
AB - Least-squares reverse time migration (LSRTM) estimates a broadband reflectivity image by iteratively minimizing the data misfit. Despite its superior imaging performance over conventional reverse time migration (RTM), LSRTM is fundamentally limited by the first-order Born approximation and incurs a high computational cost due to the iterative solution of a source-dependent model-space Hessian. We propose dual RTM, a seismic imaging framework that reformulates the inverse problem in the dual space, where the optimization variable is the scattering information of all orders, rather than being limited to the first-order scattering of the Born approximation. Unlike LSRTM, which directly estimates the reflectivity model, dual RTM reconstructs the missing scattering data (dual variable) through a fixed-point iteration governed by the full wave equation, enabling the recovery of higher-order scattering and multiple reflections without repeated background velocity updates. Furthermore, dual RTM replaces the source-dependent model-space Hessian with a source-independent data-space Hessian. A hybrid time--frequency formulation enables this Hessian to be assembled and inverted once in the Fourier domain as an offline preprocessing step, reducing the online computational cost of each dual iteration to that of a single RTM application. Numerical experiments on challenging synthetic benchmarks demonstrate that dual RTM suppresses acquisition footprints, improves amplitude fidelity, and enhances deep imaging where primary reflections alone provide insufficient illumination. These results establish dual RTM as an efficient bridge between conventional RTM/LSRTM and nonlinear waveform inversion for high-resolution large-scale seismic imaging.
ER -