PML absorbing boundary condition for structure-preserving seismic wave modeling

Abstract

Some efficient structure preserving algorithms have been proposed for high precision seismic exploration and detection in recent years. However, feasible artificial boundary conditions to truncate unbounded medium are never discussed in structure preserving seismic modeling. We pay our emphasis on application of perfectly matched layer absorbing boundary condition (PML ABC) to fit the structure preserving algorithms. In the domain of interest, a typical structure preserving algorithm is presented to discretize the conservative dynamic system, while in the PML region, traditional central difference is applied to solve the dissipative system. In the adjacent region of the interior and exterior domains, the temporary values of wave field between two time steps are approximated by interpolation using the present and the updated value of the wave fields in the PML region. A numerical experiment is adopted to demonstrate compatibility of the algorithms and no visible reflections of outgoing waves occur on the edge of interior domain.