Effective Diffusivities of Point-Like Molecules in Isotropic Porous Media by Monte Carlo Simulation

Monte Carlo simulations of point-like molecules in random and structured media are used to determine and characterize the effective diffusion coefficients of the molecules in the media. Simulations were carried out in 2D and 3D media. Monte Carlo simulation results in 2D and 3D media are compared with those obtained by analytical techniques. Simulation results indicate that for the structured, isotropic media the effective diffusivities can be characterized according to percolation thresholds in addition to porosity. The effective diffusivities in two isotropic media with the same porosity but different percolation thresholds can differ significantly. The effects of dimensionality on the effective diffusivities can also be significant. It is shown that in general the effective diffusion coefficients obtained from 2D simulation are not a good approximation to those of 3D, especially when the percolation thresholds of the 2D media and the 3D media are very different.     

General Anisotropic Effective Medium Theory for the Effective Permeability of Heterogeneous Reservoirs

One of the techniques to calculate the effective property of a heterogeneous medium is the effective medium theory. The present paper presents a general mathematical formulation for the effective medium approximation using a self-consistent choice of the effective permeability, to apply it to the case of a general anisotropic 2D medium and to the case of a 3D isotropic medium with randomly oriented ellipsoidal inclusions. The 2D results are compared with analytical results and with a homogenization technique with good result. The 3D correlations are used to derive percolation thresholds in two-phase systems with a large permeability contrast, which are compared to numerical results from the literature, also with good results.     

Percolation Thresholds in 2-Dimensional Prefractal Models of Porous Media*

Considerable effort has been directed towards the application of percolation theory and fractal modeling to porous media. We combine these areas of research to investigate percolation in prefractal porous media. We estimated percolation thresholds in the pore space of homogeneous random 2-dimensional prefractals as a function of the fractal scale invariance ratio b and iteration level i. The percolation thresholds for these simulations were found to increase beyond the 0.5927l... porosity expected in Bernoulli (uncorrelated) percolation networks. Percolation in prefractals occurs through large pores connected by small pores. The thresholds increase with both b (a finite size effect) and i. The results allow the prediction of the onset of percolation in models of prefractal porous media and can be used to bound modeling efforts. More fundamental applications are also possible. Only a limited range of parameters has been explored empirically but extrapolations allow the critical fractal dimension to be estimated for a large combination of b and i values. Extrapolation to infinite iterations suggests that there may be a critical fractal dimension of the solid at which the pore space percolates. The extrapolated value is close to 1.89 – the well-known fractal dimension of percolation clusters in 2-dimensional Bernoulli networks.     

Archie’s Law in Microsystems

Since 1942 Archie??s law is used every day to estimate, from electrical measurements, the quantity of oil present in oil fields. In this article, we perform the first experimental analysis of electric conductivity in well controlled models of porous media. We used microfluidic networks (called micromodels in the oil industry jargon), incorporating thousands of pores, with controlled wettability. Different electrode and pore geometries are considered. In all cases the evolution of the conductivity with the conductive fluid fraction (??saturation??) clearly reveals the presence of percolation thresholds, signaling that as the fraction of the conductive fluid decreases below some critical value, there are no more pathways involving only channels entirely filled with the conductive fluid that connect the electrodes. This behavior is observed in all cases, for all the network/electrode geometries and wetting properties we investigated, and is consequently likely to reflect a genuine behavior for microfluidic ??2D?? networks. The existing models??based on percolation theory or on mean field approach??reproduce correctly the structure of this behavior, but generally at a semi-quantitative level. The most successful case is obtained with the effective medium theory (EMT) model, with drainage and perpendicular electrodes. This outcome suggests that, despite the complexity of these systems, very simple models can describe correctly the physics of the system. Nonetheless, more precise modeling requires case-by-case studies. Our results are consistent with the current body of knowledge accumulated for decades on three-dimensional samples. The key point is that in 3D systems, owing to topological reasons, the threshold is extremely low in terms of water saturations. Archie??s law completely neglects the threshold effect. Nonetheless the percolation threshold should not be overlooked, and modeling should take this aspect systematically into account, as it is already done by several investigators.     

Effective Medium Analysis of Random Lattices

Calculations of effective conductivities of generalized random bond lattices representing porous media are compared to approximations using effective medium theory (EMT). We use numerical simulations of flow through 2D and 3D random lattice models, which allow for variable lattice densities and a lognormal distribution of local conductivities, to compare effective conductivities to effective medium approximations. We find that the analytical expressions provide good agreement to the simulations in 2D systems, but are in significant error in 3D systems when the standard deviation of the local conductivities is large.     

A 3D Computational Study of Effective Medium Methods Applied to Fractured Media

This work evaluates and improves upon existing effective medium methods for permeability upscaling in fractured media. Specifically, we are concerned with the asymmetric self-consistent, symmetric self-consistent, and differential methods. In effective medium theory, inhomogeneity is modeled as ellipsoidal inclusions embedded in the rock matrix. Fractured media correspond to the limiting case of flat ellipsoids, for which we derive a novel set of simplified formulas. The new formulas have improved numerical stability properties, and require a smaller number of input parameters. To assess their accuracy, we compare the analytical permeability predictions with three-dimensional finite-element simulations. We also compare the results with a semi-analytical method based on percolation theory and curve-fitting, which represents an alternative upscaling approach. A large number of cases is considered, with varying fracture aperture, density, matrix/fracture permeability contrast, orientation, shape, and number of fracture sets. The differential method is seen to be the best choice for sealed fractures and thin open fractures. For high-permeable, connected fractures, the semi-analytical method provides the best fit to the numerical data, whereas the differential method breaks down. The two self-consistent methods can be used for both unconnected and connected fractures, although the asymmetric method is somewhat unreliable for sealed fractures. For open fractures, the symmetric method is generally the more accurate for moderate fracture densities, but only the asymmetric method is seen to have correct asymptotic behavior. The asymmetric method is also surprisingly accurate at predicting percolation thresholds.     

Geometric and Hydrodynamic Characteristics of Three-dimensional Saturated Prefractal Porous Media Determined with Lattice Boltzmann Modeling

Fractal and prefractal geometric models have substantial potential for contributing to the analysis of flow and transport in porous media such as soils and reservoir rocks. In this study, geometric and hydrodynamic parameters of saturated 3D mass and pore–solid prefractal porous media were characterized using the lattice Boltzmann model (LBM). The percolation thresholds of the 3D prefractal porous media were inversely correlated with the fraction of micro-pore clusters and estimated as 0.36 and 0.30 for mass and pore–solid prefractal porous media, respectively. The intrinsic permeability and the dispersivity of the 3D pore–solid prefractals were larger than those of the 3D mass prefractals, presumably because of the occurrence of larger solid and pore cluster sizes in the former. The intrinsic permeability and dispersivity of both types of structure increased with increasing porosity, indicating a positive relationship between permeability and dispersivity, which is at odds with laboratory data and current theory. This discrepancy may be related to limitations of the convection dispersion equation at the relatively high porosity values employed in the present study.     

论岩体的渗透特性

从岩体结构的控渗作用出发, 可将岩体的渗流介质划分为多孔介质、准多孔介质、面状流不连续介质及脉状流不连续介质四类；论述了岩体复杂的渗透特性, 包括渗流的不均匀性、各向异性、非饱和性及渗流与变形的耦合等问题。文中还对岩体渗流分析的技术思路进行了讨论。     

An Algorithm for 3D Pore Space Reconstruction from a 2D Image Using Sequential Simulation and Gradual Deformation with the Probability Perturbation Sampler

The construction of a faithful 3D pore space model of a porous medium that could reproduce the macroscopic behavior of that medium is of great interest in various fields including medicine, material science, hydrology and petroleum engineering. A computationally efficient algorithm is developed that uses the probability perturbation method and sequential multiple-point statistics simulations to generate 3D stochastic and equiprobable representations of random porous media when only a 2D thin section image is available. By employing the probability perturbation method as a gradual deformation technique, the pore patterns of a single 2D image are deformed to generate a series of 2D stochastically simulated images. The 3D pore structure is then generated by simply stacking the 2D-simulated images. The quality of the 3D reconstruction is critically dependent on the rate of deformation and a simple general procedure for choosing this parameter is presented. Various criteria such as porosity, two-point auto-correlation function, multiple-point connectivity function, local percolation probability, absolute permeability obtained by lattice-Boltzmann method (LBM), formation factor and two-phase relative permeability calculations are used to validate the results. The method is tested on two random porous solids; Berea Sandstone and synthetic Silica, for which directly measured 3D micro-CT images are available. The stochastically reconstructed 3D pore space preserves the low- and high-order spatial statistics, the macroscopic flow properties and the microstructure of the 3D micro-CT images.     

Effects of the Correlation Length on the Hydraulic Parameters of a Fracture Network

We consider the influences of correlation length and aperture variability on the REV, the equivalent permeability of a fracture network, and the uncertainty in the equivalent permeability using a two-dimensional orthogonal bond percolation model. The percolation threshold, correlation length, effective conductivity, and coefficient of variation of the effective conductivity are investigated over statistically representative multiple realizations with Monte Carlo simulations in 2D fracture networks that have log-normally distributed individual fracture permeabilities. We show that although the aperture variability is large, the REV and the correlation length are similar near the percolation threshold. In contrast, when the fracture density is much larger than the percolation threshold they diverge as the aperture variability increases. We characterize the effects of correlation length and aperture variability on effective conductivity with a simple function. From the coefficient of variation analysis, the correlation length can be a criterion for evaluating which conceptual model is appropriate for describing the flow system for a given fracture network when aperture variability is sufficiently small. However, discrete fracture network models are recommended for flow simulation models because of the difficulty of REV estimation and the uncertainty in equivalent hydraulic parameters when aperture variability is large.     

Technical Note: Percolation Thresholds and Conductivities of a Uniaxial Anisotropic Simple-Cubic Lattice

We have investigated the percolation and transport behavior of uniaxial anisotropic networks, in which the bond occupation probability in one, perpendicular, direction is different from that in the other two, parallel, directions. The percolation threshold in the perpendicular direction depends on the bond occupation probability in the parallel directions, and vice versa. We report simulation results for these thresholds, and for the conductivity of finite-sized lattices, and some extrapolated estimates of percolation thresholds of infinite anisotropic lattices.     

On Determining the Percolation Reynolds Number and the Characteristic Linear Dimension for Ideal and Fictitious Porous Media

Various versions of representations of the percolation Reynolds number for porous media with isotropic and anisotropic flow properties are considered. The formulas are derived and the variants are analyzed with reference to model porous media with a periodic microstructure formed by systems of capillaries and packings consisting of spheres of constant diameter (ideal and fictitious porous media, respectively). A generalization of the Kozeny formula is given for determining the capillary diameter in an ideal porous medium equivalent to a fictitious medium with respect to permeability and porosity and it is shown that the capillary diameter is nonuniquely determined. Relations for recalculating values of the Reynolds number determined by means of formulas proposed earlier are given and it is shown that taking the microstructure of porous media into account, as proposed in [1, 2], makes it possible to explain the large scatter of the numerical values of the Reynolds number in processing the experimental data.     

Stiffness threshold of randomly distributed carbon nanotube networks

For carbon nanotube (CNT) networks, with increasing network density, there may be sudden changes in the properties, such as the sudden change in electrical conductivity at the electrical percolation threshold. In this paper, the change in stiffness of the CNT networks is studied and especially the existence of stiffness threshold is revealed. Two critical network densities are found to divide the stiffness behavior into three stages: zero stiffness, bending dominated and stretching dominated stages. The first critical network density is a criterion to judge whether or not the network is capable of carrying load, defined as the stiffness threshold. The second critical network density is a criterion to measure whether or not most of the CNTs in network are utilized effectively to carry load, defined as bending–stretching transitional threshold. Based on the geometric probability analysis, a theoretical methodology is set up to predict the two thresholds and explain their underlying mechanisms. The stiffness threshold is revealed to be determined by the statical determinacy of CNTs in the network, and can be estimated quantitatively by the stabilization fraction of network, a newly proposed parameter in this paper. The other threshold, bending–stretching transitional threshold, which signs the conversion of dominant deformation mode, is verified to be well evaluated by the proposed defect fraction of network. According to the theoretical analysis as well as the numerical simulation, the average intersection number on each CNT is revealed as the only dominant factor for the electrical percolation and the stiffness thresholds, it is approximately 3.7 for electrical percolation threshold, and 5.2 for the stiffness threshold of 2D networks. For 3D networks, they are 1.4 and 4.4. And it also affects the bending–stretching transitional threshold, together with the CNT aspect ratio. The average intersection number divided by the fourth root of CNT aspect ratio is found to be an invariant at the bending–stretching transitional threshold, which is 6.7 and 6.3 for 2D and 3D networks, respectively. Based on this study, a simple piecewise expression is summarized to describe the relative stiffness of CNT networks, in which the relative stiffness of networks depends on the relative network density as well as the CNT aspect ratio. This formula provides a solid theoretical foundation for the design optimization and property prediction of CNT networks.     

Advanced Study About the Permeability for Micro-Porous Structures Using the Lattice Boltzmann Method

The micro flows through two-dimensional (2D) and three-dimensional (3D) granular porous media at various Knudsen numbers are studied by using the lattice Boltzmann method. For 2D cases, the correlation between the permeability, the porosity, and the Knudsen number is derived. For 2D cases, the correlation can estimate the permeability well except for the staggered square cylinder. The permeability of the porous media, which have the inclusions of different sizes, is calculated. For 3D cases, simulations for the uniform overlapping and non-uniform non-overlapping granular media are carried out. The results are compared with the correlation of previous study. The effect of rarefaction on the permeability is also discussed.     

On the Convection in a Porous Medium with Inclined Temperature Gradient and Vertical Throughflow. Part II. Absolute and Convective Instabilities,and Spatially Amplifying Waves

This article presents a new methodology to estimate the effective permeability of random fractured media of any anisotropy containing both microfractures and a large number of long fractures crosscutting the representative volume element. The fractures are replaced by fictitious permeable materials for which the tangential permeability is deduced from a Poiseuille flow. A self-consistent scheme is proposed to derive the macroscopic permeability. On the one hand, the contribution of long fractures to the effective permeability writes by simple superposition of the fracture tangential permeabilities. On the other hand, the contribution of microfractures needs to resort to auxiliary problems requiring the computation of second-order Hill (or Eshelby) tensors related to ellipsoids embedded in an anisotropic matrix, for which a complete procedure is detailed. The effect of the microfracture normal permeability is put in evidence in the upscaling scheme and analyzed. In particular, it is shown that it must be chosen large enough to allow the connections between families. Examples are finally developed and compared to numerical simulations in the 2D case.     

A constitutive approach to 3-d nonlinear fluid flow through finite deformable porous continua

The present paper proposes a thermodynamically consistent Forchheimer-type filter law for application in macroscopic porous media theories. The constitutive flow equation is thereby capable of describing the essential nonlinearities during 3-d fluid percolation through deformable porous solids. In particular, tortuosity effects, anisotropic properties, and the indispensable influence of finite distortions of the interconnected pore space are accounted for. However, the common shape of a Darcy-type relation is retained by assigning all nonlinearities to a general permeability tensor. Finally, to show the validity and applicability of the proposed formulation, the filter law is correlated with the data of permeability experiments on a high-porosity polyurethane foam and is used in a 3-d finite element analysis to simulate the pneumatic damping properties of the material.     

A lattice model of foam flow in porous media: A percolation approach

Because fluid flow in porous media is opaque to most observational techniques simulations of the processes occurring in porous media have become important. Typical reservoir simulations treat the flow as taking place in some averaged (Darcy-scale) medium but simulations can also be carried out at the level of the network of pores and throats of the porous medium. We report the results of a pore-scale investigation of mechanisms for the alteration of mobility by foam lamella blockage in a network of these spaces and channels of porous media. Saturation and relative permeability curves are obtained using well-known power-law expressions of percolation theory and a rescaling of the percolation parameter readily permits a number of lamella-blocking mechanisms to be treated. An explanation of the shift in breakthrough gas saturation and the deformation of the shape of permeabilityvs saturation curves upon introduction of foam is provided for a variety of blocking mechanisms. The qualitatively different features seen in experimental studies of modification of gas mobility by foam can be rationalized using only two parameters which characterize the throat-size at which blockage commences and the degree of blockage.     

Gas cluster growth by solute diffusion in porous media. Experiments and automaton simulation on pore network

The study presented in this paper deals with the liquid–gas phase change by pressure decline of supersaturated CO₂ solutions in 2D porous media. The growth of the gas phase is studied experimentally and numerically as a function of supersaturation, wettability and gravity. Experiments are performed on a transparent etched network (micromodel) and simulations with a specific numerical automaton.In the experiments, the nucleation process, i.e. the occurrence of the gas bubbles, as well as the growth of these bubbles are visualised and analysed by means of a micro video camera and an image processing apparatus. The observations confirm the heterogeneous nature of nucleation and the disordered growth pattern of the gas phase. The analysis of the growth rate of a single gas cluster shows that this phenomenon is different from the compact growth of an isolated single bubble in the bulk. As previously predicted, the bubble growth by mass transfer and volume expansion in porous media is characterised by a pattern of the invasion percolation type under normal laboratory conditions.Numerical simulations of the growth pattern and the growth rate of a single gas cluster are performed with a numerical automaton. Based on a pore network modelling technique and a set of hypotheses derived from the observations, this automaton is first validated by comparing the numerical results with the experiments. Then, the automaton is used to conduct a sensitivity study. In particular, the influences of the Jakob number, pressure decline rate, Bond number, wettability and characteristics of the microstructure are investigated.     

An efficient method for discretizing 3D fractured media for subsurface flow and transport simulations

We introduce a new method to discretize inclined non‐planar two‐dimensional (2D) fractures in three‐dimensional (3D) fractured media for subsurface flow and transport simulations. The 2D fractures are represented by ellipsoids. We first discretize the fractures and generate a 2D finite element mesh for each fracture. Then, the mesh of fractures is analyzed by searching and treating critical geometric configurations. Based on that search, the method generates a quality mesh and allows for including finer grids. A solute transport problem in fractured porous media is solved to test the method. The results show that the method (i) adequately represents the fractured domain by maintaining the geometric integrity of input surfaces and geologic data, (ii) provides accurate results for both simple and complex fractured domains, (iii) is insensitive to spatial discretization, and (iv) is computationally very efficient. For inclined and vertical fractures, analytical and numerical solutions are shown to be in good agreement. The method is therefore suitable to discretize fracture networks for flow and transport simulations in fractured porous media. Copyright © 2010 John Wiley & Sons, Ltd.     

Calculating the Anisotropic Permeability of Porous Media Using the Lattice Boltzmann Method and X-ray Computed Tomography

A lattice Boltzmann (LB) method is developed in this article in a combination with X-ray computed tomography to simulate fluid flow at pore scale in order to calculate the anisotropic permeability of porous media. The binary 3D structures of porous materials were acquired by X-ray computed tomography at a resolution of a few microns, and the reconstructed 3D porous structures were then combined with the LB model to calculate their permeability tensor based on the simulated velocity field at pore scale. The flow is driven by pressure gradients imposed in different directions. Two porous media, one gas diffusion porous layer used in fuel cells industry and glass beads, were simulated. For both media, we investigated the relationship between their anisotropic permeability and porosity. The results indicate that the LB model is efficient to simulate pore-scale flow in porous media, and capable of giving a good estimate of the anisotropic permeability for both media. The calculated permeability is in good agreement with the measured date; the relationship between the permeability and porosity for the two media is well described by the Kozeny–Carman equation. For the gas diffusion layer, the simulated results showed that its permeability in one direction could be one order of magnitude higher than those in other two directions. The simulation was based on the single-relaxation time LB model, and we showed that by properly choosing the relaxation time, it could give similar results to those obtained using the multiple-relaxation time (MRT) LB method, but with only one third of the computational costs of MRTLB model.