Measurements of seismic attenuation and transient fluid pressure in partially saturated Berea sandstone: evidence of fluid flow on the mesoscopic scale

A novel laboratory technique is proposed to investigate wave-induced fluid flow on the mesoscopic scale as a mechanism for seismic attenuation in partially saturated rocks. This technique combines measurements of seismic attenuation in the frequency range from 1 to 100 Hz with measurements of transient fluid pressure as a response of a step stress applied on top of the sample. We used a Berea sandstone sample partially saturated with water. The laboratory results suggest that wave-induced fluid flow on the mesoscopic scale is dominant in partially saturated samples. A 3-D numerical model representing the sample was used to verify the experimental results. Biot's equations of consolidation were solved with the finite-element method. Wave-induced fluid flow on the mesoscopic scale was the only attenuation mechanism accounted for in the numerical solution. The numerically calculated transient fluid pressure reproduced the laboratory data. Moreover, the numerically calculated attenuation, superposed to the frequency-independent matrix anelasticity, reproduced the attenuation measured in the laboratory in the partially saturated sample. This experimental—numerical fit demonstrates that wave-induced fluid flow on the mesoscopic scale and matrix anelasticity are the dominant mechanisms for seismic attenuation in partially saturated Berea sandston

Digital rock physics applied to squirt flow

Meshed small and big models with triangular elements on their surfaces (.stl files) and tetrahedral elements (.mphbin)

Resolving Wave Propagation in Anisotropic Poroelastic Media Using Graphical Processing Units (GPUs)

Biot's equations describe the physics of hydromechanically coupled systems establishing the widely recognized theory of poroelasticity. This theory has a broad range of applications in Earth and biological sciences as well as in engineering. The numerical solution of Biot's equations is challenging because wave propagation and fluid pressure diffusion processes occur simultaneously but feature very different characteristic time scales. Analogous to geophysical data acquisition, high resolution and three dimensional numerical experiments lately redefined state of the art. Tackling high spatial and temporal resolution requires a high-performance computing approach. We developed a multi- graphical processing units (GPU) numerical application to resolve the anisotropic elastodynamic Biot's equations that relies on a conservative numerical scheme to simulate, in a few seconds, wave fields for spatial domains involving more than 1.5 billion grid cells. We present a comprehensive dimensional analysis reducing the number of material parameters needed for the numerical experiments from ten to four. Furthermore, the dimensional analysis emphasizes the key material parameters governing the physics of wave propagation in poroelastic media. We perform a dispersion analysis as function of dimensionless parameters leading to simple and transparent dispersion relations. We then benchmark our numerical solution against an analytical plane wave solution. Finally, we present several numerical modeling experiments, including a three-dimensional simulation of fluid injection into a poroelastic medium. We provide the Matlab, symbolic Maple, and GPU CUDA C routines to reproduce the main presented results. The high efficiency of our numerical implementation makes it readily usable to investigate three-dimensional and high-resolution scenarios of practical applications.ISSN:2169-9313ISSN:0148-0227ISSN:2169-935

Frequency-dependent AVO inversion applied to physically based models for seismic attenuation

Seismic inversion of amplitude versus offset (AVO) data in viscoelastic media can potentially provide high-resolution subsurface models of seismic velocities and attenuation from offset/angle seismic gathers. P- and S-wave quality factors (Q), whose inverse represent a measure of attenuation, depend on reservoir rock and pore fluid properties, in particular, saturation, permeability, porosity, fluid viscosity and lithology; however, these quality factors are rarely taken into account in seismic AVO inversion. For this reason, in this work, we aim to integrate quality factors derived from physically based models in AVO inversion by proposing a gradient descent optimization-based inversion technique to predict the unknown model properties (P- and S-wave velocities, the related quality factors and density). The proposed inversion minimizes the non-linear least-squares misfit with the observed data. The optimal solution is iteratively obtained by optimizing the data misfit using a second-order limited-memory quasi-Newton technique. The forward model is performed in the frequency–frequency-angle domain based on a convolution of broad-band signals and a linearized viscoelastic frequency-dependent AVO (FAVO) equation. The optimization includes the adjoint-state-based gradients with the Lagrangian formulation to improve the efficiency of the non-linear seismic FAVO inversion process. The inversion is tested on synthetic seismic data, in 1-D and 2-D, with and without noise. The sensitivity for seismic quality factors is evaluated using various rock physics models for seismic attenuation and quality factors. The results demonstrate that the proposed inversion method reliably retrieves the unknown elastic and an-elastic properties with good convergence and accuracy. The stability of the inverse solution especially seismic quality factors estimation relies on the noise level of the seismic data. We further investigate the uncertainty of the solution as a function of the variability of the initial models.Frequency-dependent AVO inversion applied to physically based models for seismic attenuationpublishedVersio