Least‐squares migration has been shown to improve image quality compared to the conventional migration method, but its computational cost is often too high to be practical. In this paper, we develop two numerical schemes to implement least‐squares migration with the reverse time migration method and the blended source processing technique to increase computation efficiency. By iterative migration of supergathers, which consist in a sum of many phase‐encoded shots, the image quality is enhanced and the crosstalk noise associated with the encoded shots is reduced. Numerical tests on 2D HESS VTI data show that the multisource least‐squares reverse time migration (LSRTM) algorithm suppresses migration artefacts, balances the amplitudes, improves image resolution and reduces crosstalk noise associated with the blended shot gathers. For this example, the multisource LSRTM is about three times faster than the conventional RTM method. For the 3D example of the SEG/EAGE salt model, with a comparable computational cost, multisource LSRTM produces images with more accurate amplitudes, better spatial resolution and fewer migration artefacts compared to conventional RTM. The empirical results suggest that multisource LSRTM can produce more accurate reflectivity images than conventional RTM does with a similar or less computational cost. The caveat is that the LSRTM image is sensitive to large errors in the migration velocity model.

Reverse time migration (RTM) exhibits great advantages over other imaging methods because it is based on computing numerical solutions to a two-way wave equation. It does not suffer from dip limitation like one-way downward continuation techniques do, thus enabling overturned reflections to be imaged. As well as correctly handling multipathing, RTM has the potential to image internal multiples when the boundaries responsible for generating the multiples are present in the model. In isotropic media, one can use a scalar acoustic wave equation for RTM of pressure data. In anisotropic media, P- and SV-waves are coupled together so, formally, elastic wave equations must be used for RTM. A new wave equation for P-waves is proposed in tilted transversely isotropic (TTI) media that can be solved as part of an acoustic anisotropic RTM algorithm, using standard explicit finite differencing. If the shear velocity along the axis of symmetry is set to zero, stable numerical solutions can be computed for media with a vertical axis of symmetry and [Formula: see text] not less than [Formula: see text]. In TTI media with rapid variations in the direction of the axis of symmetry, setting the shear velocity along the axis of symmetry to zero can cause numerical solutions to become unstable. A solution to this problem is proposed that involves using a small amount of nonzero shear velocity. The amount of shear velocity added is chosen to remove triplications from the SV wavefront and to minimize the anisotropic term of the SV reflection coefficient. We show modeling and high-quality RTM results in complex TTI media using this equation.

Reverse time migration (RTM) images reflectors by using time-extrapolation modeling codes to synthesize source and receiver wavefields in the subsurface. Asymptotic analysis of wave propagation in transversely isotropic (TI) media yields a dispersion relation describing coupled P- and SV-wave modes. This dispersion relation can be converted into a fourth-order scalar partial differential equation (PDE). Increased computational efficiency can be achieved using equivalent coupled second-order PDEs. Analysis of the corresponding dispersion relations as matrix eigenvalue systems allows one to characterize all possible coupled linear second-order systems equivalent to a given linear fourth-order PDE and to determine which ones yield optimally efficient finite-difference implementations. Setting the shear velocity along the axis of symmetry to zero yields a simpler approximate TI wave equation that is more efficient to implement. This simpler approximation, however, can become unstable for some plausible combinations of anisotropic parameters. The same eigensystem analysis can be applied using finite vertical shear velocity to obtain solutions that avoid these instability problems.

No abstract