Split PML (SPML) for acoustics

Next: Nonsplit Convolutional-PML (CPML) for Up: Perfectly Matched Layer (PML) Previous: Perfectly Matched Layer (PML)

Split PML (SPML) for acoustics

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

$\begin{equation*}\left\{ \begin{split}p=&p_x +p_z\\ \frac{\partial p_x}{\partial... ...}{\partial t}=& \frac{\partial p}{\partial z} \end{split} \right.\end{equation*}$

(22)

To absorb the boundary reflection, by adding the decaying coefficients

the SPML governing equation can be specified as (Collino and Tsogka, 2001)

$\begin{equation*}\left\{ \begin{split}p=p_x +p_z&\\ \frac{\partial p_x}{\partial... ...l t}+d(z)v_z &= \frac{\partial p}{\partial z} \end{split} \right.\end{equation*}$

(23)

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):

$\displaystyle d(u)=d_0(\frac{u}{L})^2, d_0=-\frac{3v}{2L}\ln(R), u=x,z$

(24)

where

indicates the PML thinkness;

represents the distance between current position (in PML) and PML inner boundary.

is always chosen as $10^{-3}\sim 10^{-6}$ . 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 Figure 9. Wavefield snap of SPML with 8th order finite difference

Next: Nonsplit Convolutional-PML (CPML) for Up: Perfectly Matched Layer (PML) Previous: Perfectly Matched Layer (PML)

2021-08-31