  The general framework of LSRTM consists of two steps; the first one is
generating the RTM image and the second is applying the Least-Squares
Migration, however, the convergence of both operations consumes a lot of time
to extract the final Least-Squares Reverse-Time Migration image and moreover
generates oversampling when simulating the data. Applying Reverse-Time
Migration to seismic data will generate results with some migration artifacts
depending on the applied imaging conditions. To overcome this dilemma, the
Least-Squares Reverse-Time Migration is applied to the migrated section through
Born modeling and Conjugate Gradient algorithm. Vertical transverse isotropy
(VTI) media yielded as the velocity decreases with depth which distorts the
Reverse-Time Migration results significantly. This problem can be overcome by
applying the Least-Squares Reverse-Time Migration in either the Cartesian or
pseudo-depth domains by applying a proper wavefield extrapolator. Extrapolation
of Least-Squares Reverse-Time Migration reconstructed wavefield using the 2D
constant-density acoustic wave equation transformed into Riemannian domain
treats the oversampling effect of seismic signals by making even sampling and
allows more amplitude to be recovered in the final migrated image...At a
reduced cost, the Finite Difference Riemannian wavefield extrapolator acts on
the Born modelled seismic data, producing accurately similar results to the
classical LSRTM, yet some amplitude differences are appeared due to various
implementation issues and oversampling effect in the latter. The results
support that the domain transformation strategy effectively reduces the
computational time without affecting the accuracy of the conventional LSRTM
results.
