Least squares migration compensates for the effects of missing data, noise, and illumination by imposing various constraints on an image while ensuring the model fits the observed data. Multicomponent seismic data are well suited for least squares migration as they generally suffer from many of the same complications as single component data. This article extends least squares wave equation migration to two component elastic data in isotropic media. Forward and adjoint operators are written using Helmholtz recomposition/decomposition operators implemented in the Fourier domain, while extrapolation is carried out using a split-step operator. Poynting vectors calculated using source and receiver side P-wave potentials are used to calculate angle gathers. We regularize the inversion by dip filtering in the angle domain to reduce the effect of source/receiver sampling, noise, and PP/PS crosstalk artifacts.


Article metrics loading...

Loading full text...

Full text loading...


  1. Etgen, J.T.
    [1988] Elastic prestack migration of two component data. Stanford Exploration Project, 57, 289–306.
    [Google Scholar]
  2. Higginbotham, J. and Brown, M.
    [2009] Wave equation simulation and angle-domain imaging for acquisition testing, velocity analysis, and more. 71st EAGE Conference and Exhibition.
    [Google Scholar]
  3. Kaplan, S., Routh, P. and Sacchi, M.
    [2010] Derivation of forward and adjoint operators for least-squares shot-profile split-step migration. Geophysics, 75(6), S225–S235, doi:10.1190/1.3506146.
    https://doi.org/10.1190/1.3506146 [Google Scholar]
  4. Martin, G.S., Wiley, R. and Marfurt, K.J.
    [2006] Marmousi2: An elastic upgrade for marmousi. The Leading Edge, 25(2), 156–166, doi:10.1190/1.2172306.
    https://doi.org/10.1190/1.2172306 [Google Scholar]
  5. Nemeth, T., Wu, C. and Schuster, G.
    [1999] LeastâĂŘsquares migration of incomplete reflection data. Geophysics, 64(1), 208–221, doi:10.1190/1.1444517.
    https://doi.org/10.1190/1.1444517 [Google Scholar]
  6. Scales, J.
    [1987] Tomographic inversion via the conjugate gradient method. Geophysics, 52(2), 179–185, doi: 10.1190/1.1442293.
    https://doi.org/10.1190/1.1442293 [Google Scholar]

Data & Media loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error