DYNAMIC BEHAVIOR

The kinematics of the wavefield are not affected by the introduction of a short operator to compute the model derivatives nor by using equivalent medium theory in the evaluation of the spatial operator. I introduced these changes to better simulate the dynamics of the wave-propagation problem. A simple model with horizontal layers was used to compare the dual-operator approach with a traditional finite-difference algorithm, and with an elastic modeling algorithm based on the Haskell-Thomson propagator-matrix method (Aki and Richards, 1980), which is restricted to horizontally-layered models. To test the stability and accuracy of these three schemes in extreme conditions, a water layer was placed at the top of the model.