A numerical tour of wave propagation |

It is possible for us to split the wave field into two components: x-component and z-component (Carcione et al., 2002). Then the acoustic wave equation becomes

To absorb the boundary reflection, by adding the decaying coefficients the SPML governing equation can be specified as (Collino and Tsogka, 2001)

where and are the ABC coefficients designed to attenuate the reflection in the boundary zone, see Figure 8. There exists many forms of ABC coefficients function. In the absorbing layers, we use the following model for the damping parameter (Collino and Tsogka, 2001):

(24) |

where indicates the PML thinkness; represents the distance between current position (in PML) and PML inner boundary. is always chosen as . It is important to note that the same idea can be applied to elastic wave equation (Collino and Tsogka, 2001). The split version of wave equation is very suitable for the construction of seismic Poynting vector. A straightforward application is the angle gather extration using Poynting vector, see Section .

A numerical example of SPML using 8th order staggered finite difference scheme is given in Figure 9.

snapspml
Wavefield snap of SPML with 8th order finite difference
Figure 9. |
---|

A numerical tour of wave propagation |

2021-08-31