Stochastic theory of dynamic permeability in poroelastic media

Dynamic Relative Permeabilities for Partially Saturated Porous Media Accounting for Viscous Coupling Effects: An Analytical Solution

We present an analytical model to compute frequency-dependent relative permeability functions for partially saturated porous media accounting for viscous coupling effects. For this, we consider the oscillatory motion of two immiscible fluid phases and solve the Navier–Stokes equations at the pore scale using suitable interface conditions between fluids. These calculations are combined with the generalized two-phase flow Darcy equations to obtain the corresponding upscaled macroscopic fluxes. By means of an analog pore model consisting of a bundle of cylindrical capillaries in which pore fluids are distributed in a concentric manner, we find closed analytical expressions for the complex-valued and frequency- and saturation-dependent relative permeability functions. These expressions allow for a direct assessment of viscous coupling effects on oscillatory flow for all frequencies and saturations. Our results show that viscous coupling effects significantly affect flow characteristics in the viscous and inertial regimes. Dynamic relative permeabilities are affected by the pore fluid densities and viscosities. Moreover, viscous coupling effects may induce two critical frequencies in the dynamic relative permeability curves, a characteristic that cannot be addressed by extending the classic dynamic permeability definition to partially saturated scenarios using effective fluids. The theoretical derivations and results presented in this work have implications for the estimation and interpretation of seismic and seismoelectric responses of partially saturated porous media.

Wettability effect on wave propagation in saturated porous medium

Micro-fluid mechanics studies have revealed that fluid slip on the boundary of a flow channel is a quite common phenomenon. In the context of a fluid-saturated porous medium, this implies that the fluid slippage increases with the increase of the hydrophobicity, which is the non-wetting degree. Previous studies find that wettability of the pore surface is strongly related to the slippage, which is characterized by slip length. To accurately predict acoustical properties of a fluid-saturated porous medium for different wettability conditions, the slippage of the wave-induced flow has to be taken into account. This paper introduces the slip length as a proxy for wettability into the calculation of the viscous correction factor, dynamic permeability, and dynamic tortuosity of the Biot theory for elastic waves in a porous medium. It demonstrates that, under different wettability conditions, elastic waves in a saturated porous medium have different phase velocity and attenuation. Specifically, it finds that increasing hydrophobicity yields a higher phase velocity and attenuation peak in a high-frequency range.

Seismic wave attenuation and dispersion due to wave-induced fluid flow in rocks with strong permeability fluctuations

Oscillatory fluid movements in heterogeneous porous rocks induced by seismic waves cause dissipation of wave field energy. The resulting seismic signature depends not only on the rock compressibility distribution, but also on a statistically averaged permeability. This so-called equivalent seismic permeability does not, however, coincide with the respective equivalent flow permeability. While this issue has been analyzed for one-dimensional (1D) media, the corresponding two-dimensional (2D) and three-dimensional (3D) cases remain unexplored. In this work, this topic is analyzed for 2D random medium realizations having strong permeability fluctuations. With this objective, oscillatory compressibility simulations based on the quasi-static poroelasticity equations are performed. Numerical analysis shows that strong permeability fluctuations diminish the magnitude of attenuation and velocity dispersion due to fluid flow, while the frequency range where these effects are significant gets broader. By comparing the acoustic responses obtained using different permeability averages, it is also shown that at very low frequencies the equivalent seismic permeability is similar to the equivalent flow permeability, while for very high frequencies this parameter approaches the arithmetic average of the permeability field. These seemingly generic findings have potentially important implications with regard to the estimation of equivalent flow permeability from seismic data.

Frequency-dependent effective hydraulic conductivity of strongly heterogeneous media

The determination of the transport properties of heterogeneous porous rocks, such as an effective hydraulic conductivity, arises in a range of geoscience problems, from groundwater flow analysis to hydrocarbon reservoir modeling. In the presence of formation-scale heterogeneities, nonstationary flows, induced by pumping tests or propagating elastic waves, entail localized pressure diffusion processes with a characteristic frequency depending on the pressure diffusivity and size of the heterogeneity. Then, on a macroscale, a homogeneous equivalent medium exists, which has a frequency-dependent effective conductivity. The frequency dependence of the conductivity can be analyzed with Biot's equations of poroelasticity. In the quasistatic frequency regime of this framework, the slow compressional wave is a proxy for pressure diffusion processes. This slow compressional wave is associated with the out-of-phase motion of the fluid and solid phase, thereby creating a relative fluid-solid displacement vector field. Decoupling of the poroelasticity equations gives a diffusion equation for the fluid-solid displacement field valid in a poroelastic medium with spatial fluctuations in hydraulic conductivity. Then, an effective conductivity is found by a Green's function approach followed by a strong-contrast perturbation theory suggested earlier in the context of random dielectrics. This theory leads to closed-form expressions for the frequency-dependent effective conductivity as a function of the one- and two-point probability functions of the conductivity fluctuations. In one dimension, these expressions are consistent with exact solutions in both low- and high-frequency limits for arbitrary conductivity contrast. In 3D, the low-frequency limit depends on the details of the microstructure. However, the derived approximation for the effective conductivity is consistent with the Hashin-Shtrikman bounds.