Adaptive-coefficient finite difference frequency domain method for time fractional diffusive-viscous wave equation arising in geophysics

The diffusive-viscous wave (DVW) equation is a widely used model to describe frequency dependent attenuation of seismic wave in fluid-saturated porous medium. In this paper taking power law frequency dependent attenuation into account, we first introduce a modified DVW equation (time fractional DVW equation) in which the first order temporal derivative of viscous term is replaced with a fractional order temporal derivative. In consideration of that most of the existing numerical simulations for seismic wave equations are based on time domain methods and truncation with some specific boundary conditions, we incorporate the absorbing boundary condition as complex-frequency-shifted (CFS) perfectly matched layer (PML) into the time fractional DVW equation, and then develop an adaptive-coefficient (AC) finite difference frequency domain (FDFD) method for numerical simulation. The corresponding analytical solution for homogeneous time fractional DVW equation is provided for model validation, and the effectiveness of the developed AC FDFD method is verified by some numerical examples including the homogeneous model and the layered model. Numerical results show that AC FDFD method is more accurate than the traditional 2nd-order FDFD method for numerical modeling of time fractional DVW equation with CFS PML absorbing boundary condition, while requiring similar computational costs.