A Fractional Hybrid Staggered-Grid Grünwald–Letnikov Method for Numerical Simulation of Viscoelastic Seismic Wave Propagation

The accurate and efficient simulation of seismic wave energy dissipation and phase dispersion during propagation in subsurface media due to inelastic attenuation is critical for the hydrocarbon-bearing distinction and improving the quality of seismic imaging in strongly attenuating geological media. The fractional viscoelastic equation, which quantifies frequency-independent anelastic effects, has recently become a focal point in seismic exploration. We have developed a novel hybrid staggered-grid Grünwald–Letnikov (HSGGL) finite difference method for solving the fractional viscoelastic equation in the time domain. The proposed method achieves accurate and computationally efficient solutions by using a staggered grid to discretize the first-order partial derivatives of the velocity–stress equations, combined with Grünwald–Letnikov finite difference discretization for the fractional-order terms. To improve the computational efficiency, we employ a preset accuracy to truncate the difference stencil, resulting in a compact fractional-order difference scheme. A stability analysis using the eigenvalue method reveals that the proposed method confers a relaxed stability condition, providing greater flexibility in the selection of sampling intervals. The numerical experiments indicate that the HSGGL method achieves a maximum relative error of no more than 0.17% compared to the reference solution (on a finely meshed domain) while being significantly faster than the conventional global FD method (GFD). In a 500 × 500 computational domain, the computation times for the proposed methods, which meet the specified accuracy levels used, are only approximately 4.67%, 4.47%, 4.44%, and 4.42% of that of the GFD method. This indicates that the novel HSGGL method has the potential as an effective forward modeling tool for understanding complex subsurface structures by employing a fractional viscoelastic equation.