Elastic Born modeling in an ocean-bottom node acquisition scenario |

The elastic isotropic wave equation in index notation reads:

where are the normal stresses, are the transverse stresses, is the source function in direction , is the spatial source location operating at time , is density and is the particle velocity in direction . The stresses are propagated using the stress-displacement relation:

where and are the Lame elastic constants.

In the staggered time grid methodology for elastic propagation (Virieux, 1986) the stresses and particle velocities are always half a time step apart. Therefore equation 2 and equation 3 are solved in alternation during the propagation.

For elastic Born modeling, these equations must be linearized. Beylkin and Burridge (1990) show a full derivation of linearized scattering for an elastic solid. Using to denote the specific volume, we have three models for the elastic isotropic case:

where , and are the smooth models, and , and are the perturbed models.

We also have the incident and scattered stress and particle velocity fields:

(5) |

The incident stress and particle velocity fields need to propagated with the smooth models so that they do not generate reflections. At each time step, the perturbed models must be multiplied by these incident fields, and then injected as an additive source function into the scattered fields. As a result of the staggered time grid, the injection must be done alternately into the stress and particle velocity fields. The scattered fields themselves must also be propagated with the smooth models, so that they do not generate any additional scattering.

Elastic Born modeling in an ocean-bottom node acquisition scenario |

2011-09-13