Theory and modeling of constant-Q P- and S-waves using fractional time derivatives

Abstract

I have developed and solved the constant-[Formula: see text] model for the attenuation of P- and S-waves in the time domain using a new modeling algorithm based on fractional derivatives. The model requires time derivatives of order [Formula: see text] applied to the strain components, where [Formula: see text] and [Formula: see text], with [Formula: see text] the P-wave or S-wave quality factor. The derivatives are computed with the Grünwald-Letnikov and central-difference fractional approximations, which are extensions of the standard finite-difference operators for derivatives of integer order. The modeling uses the Fourier method to compute the spatial derivatives, and therefore can handle complex geometries and general material-property variability. I verified the results by comparison with the 2D analytical solution obtained for wave propagation in homogeneous Pierre Shale. Moreover, the modeling algorithm was used to compute synthetic seismograms in heterogeneous media corresponding to a crosswell seismic experiment.