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

To reconstruct the modeled source wavefield in backward steps rather than read the stored history from the disk, one can reuse the same template by exchanging the role of and , that is,

We conduct the modeling (and the backward propagation in the same way due to template reuse):

where , and are and , respectively. When the modeling is finished, only the last two wave snaps ( and ) as well as the saved boundaries are required to do the backward time recursion.

As you see, RTM begs for an accurate reconstruction before applying the imaging condition using the backward propagated wavefield. The velocity model is typically extended with sponge absorbing boundary condition (ABC) (Cerjan et al., 1985) or PML and its variants (Komatitsch and Martin, 2007) to a larger size. In Figure 1, the original model size is extended to . In between is the artificial boundary ( ). Actually, the wavefield we intend to reconstruct is not the part in extended artificial boundary but the part in the original model zone . We can reduce the boundary load further (from whole to part of it ) depending on the required grids in finite difference scheme, as long as we can maintain the correctness of wavefield in . We do not care about the correctness of the wavefield neither in nor in the effective zone (i.e. the wavefield in ). Furthermore, we only need to compute the imaging condition in the zone , no concern with the part in .

fig1
Extend the model size with artificial boundary.
indicates the original model size (
).
is the extended model size
.
is the effective boundary area.
Figure 1. |
---|

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

2021-08-31