Compressing the memory variables in constant-Q viscoelastic wave propagation via an improved sum-of-exponentials approximation

TL;DR

This paper introduces an improved sum-of-exponentials approximation method to efficiently simulate viscoelastic wave propagation with minimal memory requirements, enhancing seismic modeling accuracy.

Contribution

It presents a nearly optimal SOE approximation for the fractional derivative, reducing memory variables and establishing a mathematical link to rheological models.

Findings

Accurately captures amplitude and phase changes in wave simulations.

Reduces memory storage requirements for fractional wave equations.

Demonstrates effectiveness in 3D homogeneous and inhomogeneous media.

Abstract

Earth introduces strong attenuation and dispersion to propagating waves. The time-fractional wave equation with very small fractional exponent, based on Kjartansson's constant-Q theory, is widely recognized in the field of geophysics as a reliable model for frequency-independent Q anelastic behavior. Nonetheless, the numerical resolution of this equation poses considerable challenges due to the requirement of storing a complete time history of wavefields. To address this computational challenge, we present a novel approach: a nearly optimal sum-of-exponentials (SOE) approximation to the Caputo fractional derivative with very small fractional exponent, utilizing the machinery of generalized Gaussian quadrature. This method minimizes the number of memory variables needed to approximate the power attenuation law within a specified error tolerance. We establish a mathematical equivalence between this SOE approximation and the continuous fractional stress-strain relationship, relating it to the generalized Maxwell body model. Furthermore, we prove an improved SOE approximation error bound to thoroughly assess the ability of rheological models to replicate the power attenuation law. Numerical simulations on constant-Q viscoacoustic equation in 3D homogeneous media and variable-order P- and S- viscoelastic wave equations in 3D inhomogeneous media are performed. These simulations demonstrate that our proposed technique accurately captures changes in amplitude and phase resulting from material anelasticity. This advancement provides a significant step towards the practical usage of the time-fractional wave equation in seismic inversion.