RTM using effective boundary saving: A staggered grid GPU implementation |

To combine the absorbing effects into the acoustic equation, CPML boundary condition is such a nice way that we merely need to combine two convolutional terms into the above equations:

where , are the convolutional terms of and ; , are the convolutional terms of and . These convolutional terms can be computed via the following relation:

where and . In the absorbing layers, the damping parameter we used is (Collino and Tsogka, 2001):

(15) |

where indicates the PML thickness; represent the distance between current position (in PML) and PML inner boundary. is always chosen as . For more details about the derivation of CPML, the interested readers are referred to Collino and Tsogka (2001) and Komatitsch and Martin (2007). The implementation of CPML boundary condition is easy to carry out: in each iteration the wavefield extrapolation is performed according to the first equation in (13); it follows by adding the convolutional terms in terms of (14).

RTM using effective boundary saving: A staggered grid GPU implementation |

2021-08-31