Effective viscoelastic representation of gas-hydrate bearing sediments from finite-element harmonic experiments

We present a novel numerical upscaling technique for modeling the wave response of gas-hydrate bearing sediments composed of a rock frame, gas-hydrate and water, where the hydrate consists of ice-like lattice of water molecules with methane trapped inside. These sediments are highly heterogeneous at mesoscopic scales, much smaller than the wavelength but much larger than the pore size, inducing substantial seismic wave attenuation and dispersion due to mode conversions. The proposed numerical upscaling procedure simulates the wave-induced fluid-flow loss mechanism by computing an average effective viscoelastic medium having the same behavior of the original sediment. The method determines the complex stiffness coefficients associated with the viscoelastic medium by solving numerically boundary value problems formulated in the space-frequency domain, representing compressibility and shear experiments. The procedure is applied to composite media with regions of different amounts of hydrate with patchy or periodic-layer distributions, which define an anisotropic effective viscoelastic medium, respectively. The examples demonstrate that variations in hydrate content induce strong attenuation and dispersion effects on seismic waves due to the mesoscopic loss mechanism.