Effect of fluid viscosity on elastic wave attenuation in porous rocks

Attenuation and dispersion of elastic waves in fluid-saturated rocks due to pore fluid viscosity is investigated using an idealized exactly solvable example of a system of alternating solid and viscous fluid layers. Waves in periodic layered systems at low frequencies can be studied using an asymptotic analysis of Rytov's exact dispersion equations. Since the wavelength of the shear wave in the fluid (viscous skin depth) is much smaller than the wavelength of the shear or compressional waves in the solid, the presence of viscous fluid layers requires a consideration of higher-order terms in the low-frequency asymptotic expansions. This expansion leads to asymptotic low-frequency dispersion equations. For a shear wave with the directions of propagation and of particle motion in the bedding plane, the dispersion equation yields the low-frequency attenuation (inverse quality factor) as a sum of two terms which are both proportional to frequency omega but have different dependencies on viscosity eta: one term is proportional to omega/eta the other to omegaeta. The low-frequency dispersion equation for compressional waves allows for the propagation of two waves corresponding to Biot's fast and slow waves. Attenuation of the fast wave has the same two-term structure as that of the shear wave. The slow wave is a rapidly attenuating diffusion-type wave, whose squared complex velocity again consists of two terms which scale with iomega/eta and iomegaeta. For all three waves, the terms proportional to eta are responsible for the viscoelastic phenomena (viscous shear relaxation), whereas the terms proportional to eta(-1) account for the visco-inertial (poroelastic) mechanism of Biot's type. Furthermore, the characteristic frequencies of visco-elastic omega(V) and poroelastic omega(B) attenuation mechanisms obey the relation omega(V)omega(B) = Aomega(R)(2), where omega(R) is the R resonant frequency of the layered system, and A is a dimensionless constant of order 1. This result explains why the visco-elastic and poroelastic mechanisms are usually treated separately in the context of macroscopic theories that imply omega much less than omega(R). The poroelastic mechanism dominates over the visco-elastic one when the frequency-indeperient parameter B = omega(B)/omega(V) = 12eta2/mu(s)rho(f)h(f)(2) much less than 1, and vice versa, where h(f) is the fluid layer thickness, rho(f) the fluid density, and mu(s) represents the shear modulus of the solid.