Effective Flux Boundary Conditions for Upscaling Porous Media Equations

We introduce a new algorithm for setting pressure boundary conditions in subgrid simulations of porous media flow. The algorithm approximates the flux in the boundary cell as the flux through a homogeneous inclusion in a homogeneous background, where the permeability of the inclusion is given by the cell permeability and the permeability of the background is given by the ambient effective permeability. With this approximation, the flux in the boundary cell scales with the cell permeability when that permeability is small, and saturates at a constant value when that permeability is large. The flux conditions provide Neumann boundary conditions for the subgrid pressure. We call these boundary conditions effective flux boundary conditions (EFBCs). We give solutions for the flux through ellipsoidal inclusions in two and three dimensions, assuming symmetric tensor permeabilities whose principal axes align with the axes of the ellipse. We then discuss the considerations involved in applying these equations to scale up problems in geological porous media. The key complications are heterogeneity, fluctuations at all length scales, and boundary conditions at finite scales.     

多孔介质渗透率的一种新分形模型

多孔介质的渗流特性是油气藏工程、地下水资源利用、高放废物深地质处置等实际工程领域的热门研究问题．基于分形理论及多孔介质由一束面积大小不等的椭圆形毛细管组成的假设,本文建立了流体在分形多孔介质中渗流时的绝对渗透率及相对渗透率的分形渗透率模型．结果表明,绝对渗透率是最大和最小孔隙面积、分形维数、形状因子ε的函数,且当ε =1时,本文模型可以简化成Yu与Cheng模型;而非饱和多孔介质的相对渗透率与饱和度和多孔介质微结构参数有关．将本文提出的渗透率分形模型预测与实验测量数据及其他模型结果进行对比,显示它们整体吻合很好．     

A theory of macrodispersion for the scale-up problem

Dispersion is the result, observable on large length scales, of events which are random on small length scales. When the length scale on which the randomness operates is not small, relative to the observations, then classical dispersion theory fails. The scale up problem refers to situations in which randomness occurs on all length scales, and for which classical dispersion theory necessarily fails. The purpose of this article is to present non-Fickian, theories of dispersion, which do not assume a scale separation between the randomness and the observed consequences, and which do not assume a single length scale.Porous media flow properties are heterogeneous on all length scales. The geological variation on length scales below the observational length scale can be regarded as unknown and unknowable, and thus as a random variable.We develop a systematic theory relating scaling behavior of the geological heterogeneity to the scaling behavior of the fluid dispersivity. Three qualitatively distinct regimes (Fickian, non-Fickian and nonrenormalizable) are found. The theory gives consistent answers within several distinct analytic approximations, and with numerical simulation of the equations of porous media flow.Comparison to field data is made. The use of Kriging to generate constrained ensembles for conditional simulation is discussed.     

Finite-Difference Approximation for Fluid-Flow Simulation and Calculation of Permeability in Porous Media

We introduce a finite-difference method to simulate pore scale steady-state creeping fluid flow in porous media. First, a geometrical approximation is invoked to describe the interstitial space of grid-based images of porous media. Subsequently, a generalized Laplace equation is derived and solved to calculate fluid pressure and velocity distributions in the interstitial space domain. We use a previously validated lattice-Boltzmann method (LBM) as ground truth for modeling comparison purposes. Our method requires on average 17 % of the CPU time used by LBM to calculate permeability in the same pore-scale distributions. After grid refinement, calculations of permeability performed from velocity distributions converge with both methods, and our modeling results differ within 6 % from those yielded by LBM. However, without grid refinement, permeability calculations differ within 20 % from those yielded by LBM for the case of high-porosity rocks and by as much as 100 % in low-porosity and highly tortuous porous media. We confirm that grid refinement is essential to secure reliable results when modeling fluid flow in porous media. Without grid refinement, permeability results obtained with our modeling method are closer to converged results than those yielded by LBM in low-porosity and highly tortuous media. However, the accuracy of the presented model decreases in pores with elongated cross sections.     

Multiscale, Multiphysics Network Modeling of Shale Matrix Gas Flows

We present a pore network model to determine the permeability of shale gas matrix. Contrary to the conventional reservoirs, where permeability is only a function of topology and morphology of the pores, the permeability in shale depends on pressure as well. In addition to traditional viscous flow of Hagen–Poiseuille or Darcy type, we included slip flow and Knudsen diffusion in our network model to simulate gas flow in shale systems that contain pores on both micrometer and nanometer scales. This is the first network model in 3D that combines pores with nanometer and micrometer sizes with different flow physics mechanisms on both scales. Our results showed that estimated apparent permeability is significantly higher when the additional physical phenomena are considered, especially at lower pressures and in networks where nanopores dominate. We performed sensitivity analyses on three different network models with equal porosity; constant cross-section model (CCM), enlarged cross-section model (ECM) and shrunk length model (SLM). For the porous systems with variable pore sizes, the apparent permeability is highly dependent on the fraction of nanopores and the pores’ connectivity. The overall permeability in each model decreased as the fraction of nanopores increased.     

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.     

A Simple Model for the Variation of Permeability due to Partial Saturation in Dual Scale Porous Media

The main focus of this work is to model macroscopically the effects of partial saturation upon the permeability of dual scale fibrous media made of fiber bundles when a Newtonian viscous fluid impregnates it. A new phenomenological model is proposed to explain the discrepancies between experimental pressure results and analytical predictions based on Darcy's law. This model incorporates the essential features of relative permeability but without the necessity of measuring saturation of the liquid for its prediction. The model is very relevant for the small scale industrial systems where a liquid is forced to flow through a fibrous porous medium. It requires four parameters. Two of them are the two permeability values based on the two length scales. One length scale is of the order of magnitude of the individual fiber radius and corresponds to the permeability of the completely staurated medium, the other is of the order of magnitude of the distance between the fiber bundles and corresponds to the permeability of the partially saturated medium. The other two parameters are the lengths of the two partially saturated regions of the flow domain. The two lengths of the partially saturated region and the permeability of the fully saturated flow domain can be directly measured from the experiments. The excellent agreement between the model and the experimental results of inlet pressure profile with respect to time suggests that this model may be used to describe the variation of the permeability behind a moving front in such porous media for correct pressure prediction. It may also be used to characterize the fibrous medium by determining the two different permeabilities and the relative importance of the unsaturated portion of the flow domain for a given architecture.     

A LGA model for fluid flow in heterogeneous porous media

A lattice gas automaton (LGA) model is proposed to simulate fluid flow in heterogeneous porous media. Permeability fields are created by distributing scatterers (solids, grains) within the fluid flow field. These scatterers act as obstacles to flow. The loss in momentum of the fluid is directly related to the permeability of the lattice gas model. It is shown that by varying the probability of occurrence of solid nodes, the permeability of the porous medium can be changed over several orders of magnitude. To simulate fluid flow in heterogeneous permeability fields, isotropic, anisotropic, random, and correlated permeability fields are generated. The lattice gas model developed here is then used to obtain the effective permeability as well as the local fluid flow field. The method presented here can be used to simulate fluid flow in arbitrarily complex heterogeneous porous media.     

A Semi-Analytical Model for Two Phase Immiscible Flow in Porous Media Honouring Capillary Pressure

This article describes a semi-analytical model for two-phase immiscible flow in porous media. The model incorporates the effect of capillary pressure gradient on fluid displacement. It also includes a correction to the capillarity-free Buckley–Leverett saturation profile for the stabilized-zone around the displacement front and the end-effects near the core outlet. The model is valid for both drainage and imbibition oil–water displacements in porous media with different wettability conditions. A stepwise procedure is presented to derive relative permeabilities from coreflood displacements using the proposed semi-analytical model. The procedure can be utilized for both before and after breakthrough data and hence is capable to generate a continuous relative permeability curve unlike other analytical/semi-analytical approaches. The model predictions are compared with numerical simulations and laboratory experiments. The comparison shows that the model predictions for drainage process agree well with the numerical simulations for different capillary numbers, whereas there is mismatch between the relative permeability derived using the Johnson–Bossler–Naumann (JBN) method and the simulations. The coreflood experiments carried out on a Berea sandstone core suggest that the proposed model works better than the JBN method for a drainage process in strongly wet rocks. Both methods give similar results for imbibition processes.     

Numerical Simulations of Flows in a Heterogeneous Porous Medium

Numerical simulations are employed to investigate the fluid flow and pressure loss in a heterogeneous block within a composite porous medium. The mean permeability of the heterogeneous block is seen to affect the overall effective flux significantly. The heterogeneous parameters of the permeability field, such as the correlation length and variance, affect it quite differently. Because of the channelling effects, the effective flux depends strongly on the realization of the permeability for larger correlation length. Under a specific permeability field, higher effective flux results from smaller variances. The influences of the inertial factor are found to be insignificant within the range of practical interests.     

Use of a Downscaling Procedure for Macroscopic Heat Source Modeling in Porous Media

Many porous media such as rocks have mesoscale inhomogeneities. The characteristic sizes of such inhomogeneities are much larger than the pore size but much less than the characteristic scale of the problem, such as the length of the sample on which measurements are taken. In this paper, we have solved the one-particle problem for depolarization of an ellipsoidal particle located in a porous medium with electrokinetic effect. To calculate the effective physical properties of a porous medium with many ellipsoidal inclusions, we have applied the effective field method. The application of this method allows us to take into account the texture of an inhomogeneous medium. The analysis performed has shown that three effective properties of inhomogeneous media (permeability, electroosmotic coupling coefficient and electrical conductivity) are not completely independent variables. General theory is illustrated by calculations of the effective properties of media containing spherical and spheroidal inclusions.     

An Explicit Physics-Based Model for the Transverse Permeability of Multi-Material Dual Porosity Fibrous Media

A simple, explicit model for the transverse permeability of porous media comprised of alternating layers of tows of varying permeability has been proposed. The development of the model is based on results of numerical simulations as well as an earlier developed analytical model for the permeability of such systems (Markicevic and Papathanasiou, 2003). Based on the above, as well as on physical arguments regarding the flow across such dual porosity media, we formulate an explicit model in terms of dimensionless permeabilities. Extensive validation of the model is carried out through numerical simulations using the CFD package FIDAP. A good agreement between the permeabilities predicted by the simple explicit model and those calculated numerically is found. From this model, the form of a master curve for the permeability of such system has been deduced.     

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.     

考虑热流固耦合作用的多孔介质孔隙尺度两相流动模拟

为了研究深层油气资源在岩石多孔介质内的运移过程, 使用一种基于Darcy-Brinkman-Biot的流固耦合数值方法, 结合传热模型, 完成了Duhamel-Neumann热弹性应力的计算, 实现了在孔隙模拟多孔介质内的考虑热流固耦合作用的两相流动过程. 模型通过求解Navier-Stokes方程完成对孔隙空间内多相流体的计算, 通过求解Darcy方程完成流体在岩石固体颗粒内的计算, 二者通过以动能方式耦合的形式, 计算出岩石固体颗粒质点的位移, 从而实现了流固耦合计算. 在此基础上, 加入传热模型考虑温度场对两相渗流过程的影响. 温度场通过以产生热弹性应力的形式作用于岩石固体颗粒, 总体上实现热流固耦合过程. 基于数值模型, 模拟油水两相流体在二维多孔介质模型内受热流固耦合作用的流动过程. 研究结果表明: 热应力与流固耦合作用产生的应力方向相反, 使得总应力比单独考虑流固耦合作用下的应力小; 温度的增加使得模型孔隙度增加, 但当注入温差达到150 K后, 孔隙度不再有明显增加; 温度的增加使得水相的相对渗流能力增加, 等渗点左移.      

Numerical and analytical study to estimate the effect of two length scales upon the permeability of a fibrous porous medium

A fibrous porous medium with two length scales is modeled as a bed of porous cylinders aligned perpendicular to the flow of viscous fluid. The flow behavior is described using Stokes and Darcy flow equations in the regions around (higher length scale) and within the cylinders (lower length scale) respectively. The typical ratio of higher and lower length-scale regions enable us to invoke lubrication approximation and simplify the equations to develop a closed form solution for the overall permeability of this dual-scale porous medium. A parametric analysis is performed to explore the dependence of permeability on factors such as the volumetric ratio of higher and lower length-scale regions, permeability and size of inclusions in the smaller length-scale region. The analytical model is compared with the numerical results and the trend is compared with the experiments.     

Thermodynamic Automaton Simulations of Fluid Flow and Diffusion in Porous Media

A thermodynamic automaton model of fluid flow in porous media is presented. The model is a nonrelativistic version of a Lorentz invariant lattice gas model constructed by Udey et al. (1998). In the previous model it was shown that the energy momentum tensor and the relativistic Boltzman equation can be rigorously derived from the collision and propagation rules. In the present paper we demonstrate that this nonrelativistic model can be used to accurately simulate well known results involving single phase flow and diffusion in porous media. The simulation results show that (1) one-phase flow simulations in porous media are consistent with Darcy's law; (2) the apparent diffusion coefficient decreases with a decrease in permeability; (3) small scale heterogeneity does not affect diffusion significantly in the cases considered.     

On the Permeability of Fractal Tube Bundles

The permeability of a porous medium is strongly affected by its local geometry and connectivity, the size distribution of the solid inclusions, and the pores available for flow. Since direct measurements of the permeability are time consuming and require experiments that are not always possible, the reliable theoretical assessment of the permeability based on the medium structural characteristics alone is of importance. When the porosity approaches unity, the permeability?Cporosity relationships represented by the Kozeny?CCarman equations and Archie??s law predict that permeability tends to infinity and thus they yield unrealistic results if specific area of the porous media does not tend to zero. The aim of this article is the evaluation of the relationships between porosity and permeability for a set of fractal models with porosity approaching unity and a finite permeability. It is shown that the tube bundles generated by finite iterations of the corresponding geometric fractals can be used to model porous media where the permeability?Cporosity relationships are derived analytically. Several examples of the tube bundles are constructed, and the relevance of the derived permeability?Cporosity relationships is discussed in connection with the permeability measurements of highly porous metal foams reported in the literature.     

Layered Thermohaline Convection in Hypersaline Geothermal Systems

Thermohaline convection occurs in hypersaline geothermal systems due to thermal and salinity effects on liquid density. Because of its importance in oceanography, thermohaline convection in viscous liquids has received more attention than thermohaline convection in porous media. The fingered and layered convection patterns observed in viscous liquid thermohaline convection have been hypothesized to occur also in porous media. However, the extension of convective dynamics from viscous liquid systems to porous media systems is complicated by the presence of the solid matrix in porous media. The solid grains cause thermal retardation, hydrodynamic dispersion, and permeability effects. We present simulations of thermohaline convection in model systems based on the Salton Sea Geothermal System, California, that serve to point out the general dynamics of porous media thermohaline convection in the diffusive regime, and the effects of porosity and permeability, in particular. We use the TOUGH2 simulator with residual formulation and fully coupled solution technique for solving the strongly coupled equations governing thermohaline convection in porous media. We incorporate a model for brine density that takes into account the effects of NaCl and CaCl2. Simulations show that in forced convection, the increased pore velocity and thermal retardation in low-porosity regions enhances brine transport relative to heat transport. In thermohaline convection, the heat and brine transport are strongly coupled and enhanced transport of brine over heat cannot occur because buoyancy caused by heat and brine together drive the flow. Random permeability heterogeneity has a limited effect if the scale of flow is much larger than the scale of permeability heterogeneity. For the system studied here, layered thermohaline convection persists for more than one million years for a variety of initial conditions. Our simulations suggest that layered thermohaline convection is possible in hypersaline geothermal systems provided the vertical permeability is smaller than the horizontal permeability, as is likely in sedimentary basins such as the Salton Trough. Layered thermohaline convection can explain many of the observations made at the Salton Sea Geothermal System over the years.     

Evaluation of the Darcy's law performance for two-fluid flow hydrodynamics in a particle debris bed using a lattice-Boltzmann model

In the present paper, multiphase flow dynamics in a porous medium are analyzed by employing the lattice-Boltzmann modeling approach. A two-dimensional formulation of a lattice-Boltzmann model, using a D2Q9 scheme, is used. Results of the FlowLab code simulation for single phase flow in porous media and for two-phase flow in a channel are compared with analytical solutions. Excellent agreement is obtained. Additionally, fluid-fluid interaction and fluid-solid interaction (wettability) are modeled and examined. Calculations are performed to simulate two-fluid dynamics in porous media, in a wide range of physical parameters of porous media and flowing fluids. It is shown that the model is capable of determining the minimum body force needed for the nonwetting fluid to percolate through the porous medium. Dependence of the force on the pore size, and geometry, as well as on the saturation of the nonwetting fluid is predicted by the model. In these simulations, the results obtained for the relative permeability coefficients indicate the validity of the reciprocity for the two coupling terms in the modified Darcy's law equations. Implication of the simulation results on two-fluid flow hydrodynamics in a decay-heated particle debris bed is discussed. Received on 1 December 1999     

Deletion of nondominant effects in modeling transport in porous media

A methodology for eliminating nondominant effects in models that describe transport phenomena in porous media is presented. The methodology is based on the introduction of dimensionless numbers and on a proper evaluation of the order of magnitude of terms. These dimensionless numbers are redefined as characteristics of transport and transformation phenomena in porous media. It is shown that different time scales and different length scales may have to be employed for different variables. A method for evaluating the order of magnitude of the error of prediction when terms are deleted, is presented.