A Smoothing Scheme for Seismic Wave Propagation Simulation with a Mixed FEM of Diffusionized Wave Equation

In this paper, we propose a smoothing scheme for seismic wave propagation simulation. The proposed scheme is based on a diffusionized wave equation with the fourth-order spatial derivative term. So, the solution requires higher regularity in the usual weak formulation. Reducing the diffusionized wave equation to a coupled system of diffusion equations yields a mixed FEM to ease the regularity. We mathematically explain how our scheme works for smoothing. We construct a semi-implicit time integration scheme and apply it to the wave equation. This numerical experiment reveals that our scheme is effective for filtering short wavelength components in seismic wave propagation simulation.