Algorithms for staggered-grid computations for poroelastic, elastic, acoustic, and scalar wave equations

Heterogeneous wave equations are more complicated numerically than homogeneous wave equations, but are necessary for physical validity. A wide variety of numerical solutions of seismic wave equations is available, but most produce strong numerical artefacts and local instabilities where model parameters change rapidly. Accuracy and stability of heterogeneous equations is achieved through staggered-grid formulations. A new pseudospectral staggered-grid algorithm is developed for the poroelastic (Blot) equations. The algorithm may be reduced to handle the elastic and acoustic limits of the Blot equations. Comparisons of results from poroelastic, elastic, acoustic and scalar computations for a 2D model show that porous medium parameters may affect amplitudes significantly. The use of homogeneous wave equations for modelling of a heterogeneous medium, or of a centred rather than a staggered grid, or of simplified (e.g. acoustic) wave equations when elastic or poroelastic media are synthesized, may produce erroneous or ambiguous interpretations.