Seismic wave propagation in porous rocks that are saturated with a liquid exhibits significant dispersion and attenuation due to fluid flow at the pore scale, so-called squirt flow. This phenomenon takes place in compliant flat pores such as microcracks and grain contacts that are connected to stiffer isometric pores. Accurate quantitative description is crucial for inverting rock and fluid properties from seismic attributes such as attenuation. Up to now, many analytical models for squirt flow were proposed based on simplified geometries of the pore space. These models were either not compared with a numerical solution or showed poor accuracy. We present a new analytical model for squirt flow which is validated against a three-dimensional numerical solution for a simple pore geometry that has been classically used to explain squirt flow; that is why we refer to it as classical geometry. The pore space is represented by a flat cylindrical (penny-shaped) pore whose curved edge is fully connected to a toroidal (stiff) pore. Compared with correct numerical solutions, our analytical model provides very accurate predictions for the attenuation and dispersion across the whole frequency range. This includes correct low-and high-frequency limits of the stiffness modulus, the characteristic frequency, and the shape of the dispersion and attenuation curves. In a companion paper (Part 2), we extend our analytical model to more complex pore geometries. We provide as supplementary material Matlab and symbolic Maple routines to reproduce our main results.

We explore the impact of roughness in crack walls on the P-wave modulus dispersion and attenuation caused by squirt flow. For that, we numerically simulate oscillatory relaxation tests on models having interconnected cracks with both simple and intricate aperture distributions. Their viscoelastic responses are compared with those of models containing planar cracks but having the same hydraulic aperture as the rough wall cracks. In the absence of contact areas between crack walls, we found that three apertures affect the P-wave modulus dispersion and attenuation: the arithmetic mean, minimum and hydraulic apertures. We show that the arithmetic mean of the crack apertures controls the effective P-wave modulus at the low- and high-frequency limits, thus representing the mechanical aperture. The minimum aperture of the cracks tends to dominate the energy dissipation process, and consequently, the characteristic frequency. An increase in the confining pressure is emulated by uniformly reducing the crack apertures, which allows for the occurrence of contact areas. The contact area density and distribution play a dominant role in the stiffness of the model and, in this scenario, the arithmetic mean is not representative of the mechanical aperture. On the other hand, for a low percentage of minimum aperture or in presence of contact areas, the hydraulic aperture tends to control de characteristic frequency. Analysing the local energy dissipation, we can more specifically visualise that a different aperture controls the energy dissipation process at each frequency, which means that a frequency-dependent hydraulic aperture might describe the squirt flow process in cracks with rough walls.

Understanding the properties of cracked rocks is of great importance in scenarios involving CO 2 geological sequestration, nuclear waste disposal, geothermal energy , and hydrocarbon exploration and production. Developing noninvasive detecting and monitoring methods for such geological formations is crucial. Many studies show that seismic waves exhibit strong dispersion and attenuation across a broad frequency range due to fluid flow at the pore scale known as squirt flow. Nevertheless, how and to what extent squirt flow affects seismic waves is still a matter of investigation. To fully understand its angle-and frequency-dependent behavior for specific geometries, appropriate numerical simulations are needed. We perform a three-dimensional numerical study of the fluid-solid deformation at the pore scale based on coupled Lamé-Navier and Navier-Stokes linear quasistatic equations. We show that seismic wave velocities exhibit strong azimuth-, angle-and frequency-dependent behavior due to squirt flow between interconnected cracks. Furthermore, the overall anisotropy of a medium mainly increases due to squirt flow, but in some specific planes the anisotropy can locally decrease. We analyze the Thomsen-type anisotropic parameters and adopt another scalar parameter which can be used to measure the anisotropy strength of a model with any elastic symmetry. This work significantly clarifies the impact of squirt flow on seismic wave anisotropy in three dimensions and can potentially be used to improve the geophysical monitoring and surveying of fluid-filled cracked porous zones in the subsurface.

Biot's equations describe the physics of hydro-mechanically 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 re-defined state of the art. Tackling high spatial and temporal resolution requires a high-performance computing approach. We developed a multi-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.