(2) |

When ignoring perturbations in the source function or in density, the adjoint (transpose) can be written in terms of material parameters:

(3) |

The above equation is intuitively more understandable when written in terms of wavefields:

(4) |

Changes in the material parameters can be estimated by crosscorrelating the backward-propagated modeled strain field and the backward-propagated data strain field. Both fields are tensor fields and what are crosscorrelated are actually the strain tensor components. In terms of prestack migration this crosscorrelation is the imaging condition. In the following I use this imaging condition to find perturbations of elastic stiffness parameters for a given starting model. One might ask: why image stiffnesses and not velocities? For isotropic media the imaging of velocities or wave potential works well, because the wave types are separable into a P wave and S wave potential and those potentials are decoupled. In the case of anisotropic media such a decoupling is only possible in a few rare cases of special symmetry. Imaging stiffness naturally evolves out of the anisotropic wave equation and wave types don't have to be decoupled in order to perform imaging. Some stiffness components, however, may not be be very well resolved, depending strongly on the acquisition geometry and source radiation pattern illuminating the subsurface. The above stated imaging condition reduces to the well-known U/D principle (Claerbout, 1986) in the case of an isotropic scalar wavefield.

12/18/1997