Theoretical and quantitative evaluation of hybrid PML-ABCs for seismic wave simulation

Abstract

A good artificial boundary treatment in a seismic wave grid-based numerical simulation can reduce the size of the computational region and increase the computational efficiency, which is becoming increasingly important for seismic migration and waveform inversion tasks requiring hundreds or thousands of simulations. Two artificial boundary techniques are commonly used: perfectly matched layers (PMLs), which exhibit the excellent absorption performance but impose a greater computational burden by using finite layers to gradually reduce wave amplitudes; and absorbing boundary conditions (ABCs), which have the high computational efficiency but are less effective in absorption because they employ the one-way wave equation at the exterior boundary. Naturally, PMLs have been combined with ABCs to reduce the number of PMLs, thus improving the computational efficiency; many studies have proposed such hybrid PMLs. Depending on the equations from which the ABCs are derived, there are two hybrid PML variants: the PML+unstretched ABC (UABC), in which the ABC is derived from a physical equation; or the PML+stretched ABC (SABC), in which the ABC is derived from the PML equation. Even though all the previous studies concluded that hybrid PMLs can improve the absorption performance, none of them quantified how many PMLs can be removed by combining the PML with the ABC compared with the pure PML. In this paper, we systematically study the absorption performance of the two hybrid PML variants. We develop a method to distinguish the artificial reflections from the PML-interior interface and those caused by the PML exterior boundary to accurately approximate the additional absorption achieved by using the UABC and the SABC. The reflection coefficients based on a theoretical derivation and numerical tests both show that the UABC amplifies most reflections and is not recommended in any situation; conversely, the SABC can always diminish reflections, but the additional absorption achieved by the SABC is relatively poor and cannot effectively reduce the number of PMLs. In contrast, we find that simply increasing the damping parameter improves absorption better than the PML+SABC. Our results show that the improvement in absorption achieved by combining the PML with either the SABC or the UABC is not better than that obtained by simply adjusting the damping profile of the PML; thus, combining the PML with the ABC is not recommended in practice.