Fractional order viscous sound wave equation numerical simulation method, system and equipment

The invention discloses a fractional order viscous sound wave equation numerical simulation method, system and device, and relates to the technical field of seismic wave propagation numerical simulation, and the method comprises the steps: substituting a wave equation plane wave solution into a fractional order Laplacian operator finite difference formula to obtain a generation function of the corresponding finite difference formula; replacing the regular grid finite difference coefficient with the optimized regular grid finite difference coefficient, and performing inverse Fourier transform on the generated function to obtain an optimized difference coefficient; and based on the optimized difference coefficient of the finite difference formula and the fractional order Laplacian operator finite difference formula, numerical discretization is performed on a preset fractional order Laplacian viscoacoustic wave equation to obtain a fractional order Laplacian viscoacoustic wave equation finite difference wave fie