Heterogeneities of flow in stochastically generated porous media

Effect of the Heterogeneity on Sorptivity in Sandstones with High and Low Permeability in Water Imbibition Process

Capillary imbibition in unsaturated rocks is important for the exploitation of tight reservoirs, such as oil and gas reservoirs. However, the physical properties of natural rocks tend to be relatively uneven, mainly in the heterogeneity of material composition and pore space. Reservoir heterogeneity is an important factor affecting the exploitation of oil fields and other reservoirs, which can be evaluated by the pore structure tortuosity fractal dimension DT of rock. The greater the value of DT, the stronger the heterogeneity of sandstone. Two types of sandstone with high and low permeability were selected to study the effect of heterogeneity on the imbibition behavior by using high-resolution X-ray imaging and neutron radiography. Quantitative results of the wetting front position for each specimen were extracted from the neutron images. The wetting front advanced linearly with the power index of time t1/(2DT). Different values of DT were selected to estimate and discuss the effect of the heterogeneity on sorptivity. A modified L-W equation was employed to predict the sorptivity. Comparing with the experimental results, the heterogeneity plays a significant role in determining the sorptivity. The modified model provides a reference for the prediction of the sorptivity of the same types of sandstones studied in this paper.

Pore-scale statistics of flow and transport through porous media

Flow in porous media is known to be largely affected by pore morphology. In this work, we investigate the effects of pore geometry on the transport and spatial correlations of flow through porous media in two distinct pore structures arising from three-dimensional assemblies of overlapping and nonoverlapping spheres. Using high-resolution direct numerical simulations (DNS), we perform Eulerian and Lagrangian analysis of the flow and transport characteristics in porous media. We show that the Eulerian velocity distributions change from nearly exponential to Gaussian distributions as porosity increases. A stretched exponential distribution can be used to represent this behavior for a wide range of porosities. Evolution of Lagrangian velocities is studied for the uniform injection rule. Evaluation of tortuosity and trajectory length distributions of each porous medium shows that the model of overlapping spheres results in higher tortuosity and more skewed trajectory length distributions compared to the model of nonoverlapping spheres. Wider velocity distribution and higher tortuosity for overlapping spheres model give rise to non-Fickian transport while transport in nonoverlapping spheres model is found to be Fickian. Particularly, for overlapping spheres model our analysis of first-passage time distribution shows that the transport is very similar to those observed for sandstone. Finally, using three-dimensional (3D) velocity field obtained by DNS at the pore-scale, we quantitatively show that despite the randomness of pore-space, the spatially fluctuating velocity field and the 3D pore-space distribution are strongly correlated for a range of porous media from relatively homogeneous monodisperse sphere packs to Castlegate sandstone.

Pore-scale simulation of fluid flow and solute dispersion in three-dimensional porous media

In the present work fluid flow and solute transport through porous media are described by solving the governing equations at the pore scale with finite-volume discretization. Instead of solving the simplified Stokes equation (very often employed in this context) the full Navier-Stokes equation is used here. The realistic three-dimensional porous medium is created in this work by packing together, with standard ballistic physics, irregular and polydisperse objects. Emphasis is placed on numerical issues related to mesh generation and spatial discretization, which play an important role in determining the final accuracy of the finite-volume scheme and are often overlooked. The simulations performed are then analyzed in terms of velocity distributions and dispersion rates in a wider range of operating conditions, when compared with other works carried out by solving the Stokes equation. Results show that dispersion within the analyzed porous medium is adequately described by classical power laws obtained by analytic homogenization. Eventually the validity of Fickian diffusion to treat dispersion in porous media is also assessed.

Geometrical interpretation of long-time tails of first-passage time distributions in porous media with stagnant parts

Using a combined experimental-numerical approach, we study the first-passage time distributions (FPTD) of small particles in two-dimensional porous materials. The distributions in low-porosity structures show persistent long-time tails, which are independent of the Péclet number and therefore cannot be explained by the advection-diffusion equation. Instead, our results suggest that these tails are caused by stagnant, i.e., quiescent areas where particles are trapped for some time. Comparison of measured FPTD with an analytical expression for the residence time of particles, which diffuse in confined regions and are able to escape through a small pore, yields good agreement with our data.

Relationship between pore size and velocity probability distributions in stochastically generated porous media

We perform a set of detailed numerical simulations of single-phase, fully saturated flow in stochastically generated, three-dimensional pore structures with diverse porosities (ϕ) and degrees of connectivity, and analyze the probability density functions (PDFs) of the pore sizes, S, and vertical velocity components, w, which are aligned with the mean flow direction. Both of the PDFs are markedly skewed with pronounced positive tails. This feature of the velocity PDF is dictated by the pore structure and determines the shortest travel times, one of the key transport attributes that underpins the success or the failure of environmental remediation techniques. Using a maximum likelihood approach, we determine that the PDFs of S and w decay according to an exponential and a stretched exponential model, respectively. A strong correlation between the key parameters governing the decay of the upper tails of the two PDFs is found, which provides a quantitative result for this analogy that so far has been stated only qualitatively. The parameter governing the concavity of the tail of the velocity PDF varies linearly with porosity over the entire range of tested values (0.2≤ϕ≤0.6). The parameters controlling the spread of the upper tails of the PDFs of S and w appear to be linked by a power-law relationship.

Hyperbolic regions in flows through three-dimensional pore structures

Finite time Lyapunov exponents are used to determine expanding, contracting, and hyperbolic regions in computational simulations of laminar steady-state fluid flows within realistic three dimensional pore structures embedded within an impermeable matrix. These regions correspond approximately to pores where flow converges (contraction) or diverges (expansion), and to throats between pores where the flow mixes (hyperbolic). The regions are sparse and disjoint from one another, occupying only a small percentage of the pore space. Nonetheless, nearly every percolating fluid particle trajectory passes through several hyperbolic regions indicating that the effects of in-pore mixing are distributed throughout an entire pore structure. Furthermore, the observed range of fluid dynamics evidences two scales of heterogeneity within each of these flow fields. There is a larger scale that affects dispersion of fluid particle trajectories across the connected network of pores and a relatively small scale of nonuniform distributions of velocities within an individual pore.

Anisotropy of flow in stochastically generated porous media

Models of porous media are often applied to relatively small systems, which leads not only to system-size-dependent results, but also to phenomena that would be absent in larger systems. Here we investigate one such finite-size effect: anisotropy of the permeability tensor. We show that a nonzero angle between the external body force and macroscopic flux vector exists in three-dimensional periodic models of sizes commonly used in computer simulations and propose a criterion, based on the ratio of the system size to the grain size, for this phenomenon to be relevant or negligible. The finite-size anisotropy of the porous matrix induces a pressure gradient perpendicular to the axis of a porous duct and we analyze how this effect scales with the system and grain sizes. We also analyze how the size of the representative elementary volume (REV) for anisotropy compares with the REV for permeability.