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.     

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.     

Determination of permeability in fibrous porous media using the lattice Boltzmann method with application to PEM fuel cells

The lattice Boltzmann method (LBM) is used to simulate the flow through an idealized proton exchange membrane fuel cell (PEMFC) porous transport layer (PTL) geometry generated using a Monte Carlo method. Using the calculated flow field, Darcy's law is applied and the permeability is calculated. This process is applied in both through‐ and in‐plane directions of the paper as both of these permeability values are important in computational fluid dynamics models of PEMFCs. It is shown that the LBM can be used to determine permeability in a random porous media by solving the flow in the microstructure of the material. The permeability in the through‐ and in‐plane directions is shown to be different and the anisotropic nature of the geometry creates anisotropic permeability. It is also found that fiber arrangement plays a large role in the permeability of the PTL. New correlations are presented for in‐ and though‐plane permeabilities of fibrous porous media with (0.6<ε<0.8). Copyright © 2008 John Wiley & Sons, Ltd.     

A LGA model for dispersion in heterogeneous porous media

The lattice gas automaton (LGA) model proposed in the previous paper is applied to the problem of simulating dispersion and mixing in heterogeneous porous media. We demonstrate here that tracer breakthrough profiles and longitudinal dispersion coefficients can be computed for heterogeneous porous media.     

Knudsen’s Permeability Correction for Tight Porous Media

Various flow regimes including Knudsen, transition, slip and viscous flows (Darcy’s law), as applied to flow of natural gas through porous conventional rocks, tight formations and shale systems, are investigated. Data from the Mesaverde formation in the United States are used to demonstrate that the permeability correction factors range generally between 1 and 10. However, there are instances where the corrections can be between 10 and 100 for gas flow with high Knudsen number in the transition flow regime, and especially in the Knudsen’s flow regime. The results are of practical interest as gas permeability in porous media can be more complex than that of liquid. The gas permeability is influenced by slippage of gas, which is a pressure-dependent parameter, commonly referred to as Klinkenberg’s effect. This phenomenon plays a substantial role in gas flow through porous media, especially in unconventional reservoirs with low permeability, such as tight sands, coal seams, and shale formations. A higher-order permeability correlation for gas flow called Knudsen’s permeability is studied. As opposed to Klinkenberg’s correlation, which is a first-order equation, Knudsen’s correlation is a second-order approximation. Even higher-order equations can be derived based on the concept used in developing this model. A plot of permeability correction factor versus Knudsen number gives a typecurve. This typecurve can be used to generalize the permeability correction in tight porous media. We conclude that Knudsen’s permeability correlation is more accurate than Klinkenberg’s model especially for extremely tight porous media with transition and free molecular flow regimes. The results from this study indicate that Klinkenberg’s model and various extensions developed throughout the past years underestimate the permeability correction especially for the case of fluid flow with the high Knudsen number.     

Multiresolution Wavelet Scale Up of Unstable Miscible Displacements in Flow Through Heterogeneous Porous Media

We describe scale up of geological models of field-scale porous media using a new method based on the wavelet transformations. The porous media of interest contain broadly-distributed and correlated permeabilities. Wavelet transformation of the permeability field of such porous media coarsens the geological model from smallest to largest length scales, drastically reduces the number of equations to be solved, preserves the important information on the permeability field at all the relevant length scales, and yields numerical results for any fluid flow property that are as accurate as those that are obtained with the highly detailed geological model of the same porous media. To test this method, we carry out extensive computer simulations of unstable miscible displacement processes and the associated viscous fingering phenomenon in highly heterogeneous porous media, both with the fine-scale geological model and the coarsened model. Excellent agreement is found between the results of the two sets of simulations.     

裂缝性低渗透油藏流-固耦合理论与数值模拟

根据裂缝性低渗油藏的储层特征,建立适合裂缝性砂岩油藏渗流的等效连续介质模型。将渗流力学与弹塑性力学相结合,建立裂缝性低渗透油藏的流-固耦合渗流数学模型,并给出其数值解.通过数值模拟对一实际井网开发过程中孔隙度、渗透率的变化以及开发指标进行计算,并和刚性模型以及双重介质模型的计算结果进行了分析比较.     

考虑气体压缩性的多孔材料渗透率和惯性系数的测定

多孔介质材料的渗透率和惯性系数是决定多孔介质中流体流动特性的重要参数，一般需要通过实验进行测定，在测定渗透律和惯性系数量，选用气体作为工作介质可以为实验带来极大的方便，然而通常的实验都将气体看作不可压缩流体，直接根据Darcy-Forchheimer定律得到这两个参数，这种近似对实验条件如样品厚度、工作压力等提出了很多要求，本文提出了在考虑气体压缩性的情况下测定渗透率和惯性系数的方法，该方法可以大大降低实验时对样品厚度、工作压力等条件的要求。本文还根据该方法对多孔材料PVF进行了渗透率和惯性系数的测定，并对测量结果进行了验证。     

Upscaling LBM-TPM simulation approach of Darcy and non-Darcy fluid flow in deformable,heterogeneous porous media

This paper presents a numerical approach for the simulation of fluid flow through porous media by proposing a theoretical and numerical meso-to-macro multiscale framework, which combines the advantages of the lattice Boltzmann method (LBM) with the continuum Theory of Porous Media (TPM) to efficiently and accurately model fluid transport in heterogeneous porous media. In particular, LBM presents an alternative to experiments by studying the flow from a mesoscopic perspective, which in turn, allows the derivation of the material parameters needed for simulating the flow in the macroscopic TPM model. In this work, a meso-macro hierarchic upscaling scheme is applied to investigate the deformation-dependent intrinsic permeability properties and the Darcy/non-Darcy fluid flow regime. Concerning the mesoscale, the intrinsic permeability of the porous domain is computed by means of the LBM model at the first stage. Subsequently, deformation of the medium takes place in furtherance of determining the relation of the aforementioned deformation dependency. Thereupon, these findings are input into the TPM model in order to compute the primary unknown variables, where special focus is laid on the stability challenges in the compaction and near compaction states. With respect to the criteria of non-Darcy fluid flow, the conditions of its onset, i.e. the induced pressure gradient and mean fluid flow velocity, are computed as well using the LBM solver and conveyed afterwards to the macroscopic TPM model. Herein, the non-Darcy intrinsic permeability has been investigated in the TPM approach based on the Forchheimer equation. Simulations done on a synthetic porous micro-structure show that the combined framework proved to stand well between the two approaches.     

Analytical Solutions to Gas Flow Problems in Low Permeability Porous Media

In most of conventional porous media the flow of gas is basically controlled by the permeability and the contribution of gas flow due to gas diffusion is ignored. The diffusion effect may have significant impact on gas flow behavior, especially in low permeability porous media. In this study, a dual mechanism based on Darcy flow as well as diffusion is presented for the gas flow in homogeneous porous media. Then, a novel form of pseudo pressure function was defined. This study presents a set of novel analytical solutions developed for analyzing steady-state and transient gas flow through porous media including effective diffusion. The analytical solutions are obtained using the real gas pseudo pressure function that incorporates the effective diffusion. Furthermore, the conventional assumption was used for linearizing the gas flow equation. As application examples, the new analytical solutions have been used to design new laboratory and field testing method to determine the porous media parameters. The proposed laboratory analysis method is also used to analyze data from steady-state flow tests of three core plugs. Then, permeability (k) and effective diffusion coefficient (D _e) was determined; however, the new method allows one to analyze data from both transient and steady-state tests in various flow geometries.     

Pore-Scale Simulation of Shear Thinning Fluid Flow Using Lattice Boltzmann Method

The present work attempts to identify the roles of flow and geometric variables on the scaling factor which is a necessary parameter for modeling the apparent viscosity of non-Newtonian fluid in porous media. While idealizing the porous media microstructure as arrays of circular and square cylinders, the present study uses multi-relaxation time lattice Boltzmann method to conduct pore-scale simulation of shear thinning non-Newtonian fluid flow. Variation in the size and inclusion ratio of the solid cylinders generates wide range of porous media with varying porosity and permeability. The present study also used stochastic reconstruction technique to generate realistic, random porous microstructures. For each case, pore-scale fluid flow simulation enables the calculation of equivalent viscosity based on the computed shear rate within the pores. It is observed that the scaling factor has strong dependence on porosity, permeability, tortuosity and the percolation threshold, while approaching the maximum value at the percolation threshold porosity. The present investigation quantifies and proposes meaningful correlations between the scaling factor and the macroscopic properties of the porous media.     

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.     

Stochastic Effects of 2D Random Arrays of Cylinders on Rarefied Gas Permeability Using the Lattice Boltzmann Method

A large set of 2D random arrays of circular cylinders is generated to perform a statistical study on rarefied gas flow through micro-porous media. The flow regimes in this work lie for Knudsen numbers (Kn) ranging from the continuum to the transition regimes. Arrays are built by randomly placing cylinders with constant diameter with a uniform distribution without overlapping, and are generated for three target porosities. Fluid flow is assumed to be incompressible and isothermal. A modified lattice Boltzmann model is adopted to account for discrete effects, with slip-velocity boundary conditions and a Kn-dependent multi-relaxation time collision operator. The apparent permeability is modeled with Darcy’s law with a Klinkenberg-type relationship and compared with existing correlations. Velocity fields highlight the increasing contribution of fluid flow through small pores with increasing Kn. Numerical results show that porous media randomness leads to an uncertainty on rarefied gas permeability calculation despite the same structural characteristics and may not strictly follow a specific correlation. The influence of a local collision operator based on a local Kn instead of a global one in the numerical model is also studied. Results show that the permeability in rarefied regimes undergoes significant deviation when applying the local collision operator compared to the global one. These differences could result from a more accurate capture of the pore-scale behavior with a local Kn. Thus, it emphasizes the sensitivity of the model and the apparent permeability calculation to the appropriate definition of Kn.

     

低渗透多孔介质中的非线性渗流理论简

文中论述了低渗透性多孔介质中非线性渗流理论的几个问题,阐明了渗流流体的性质,指出了多孔介质对流体通过的选择性,提出了新的非线性渗流方程,用实验资料对其进行了验证,分析了该方程演变功能,表明它可以描述各种渗流规律.该方程的各项参数都可从实验中直接得到,应用方便,并且参数的物理意义明确.     

Fluid flow through heterogeneous porous media in a rotating square channel

An analytical three-dimensional solution to the fluid flow problem through heterogeneous porous media in a rotating square channel is presented. The permeability of the fluid saturated porous domain varies in the vertical direction, thus affecting the imposed main flow in the channel. As a result of Coriolis acceleration, secondary circulation in a plane perpendicular to the main flow direction is created. A particular example of a monotonic distribution of the permeability function is analyzed leading to a single vortex secondary circulation. Nevertheless, multiple vortex secondary flow solutions are possible depending on the particular variation of the permeability in the vertical direction. No secondary motion is expected for isothermal flows in homogeneous porous media.     

LBM Simulation of Viscous Fingering Phenomenon in Immiscible Displacement of Two Fluids in Porous Media

This article presents the lattice Boltzmann simulation of viscous fingering phenomenon in immiscible displacement of two fluids in porous media. Such phenomenon generally takes place when a less viscous fluid is used to displace a more viscous fluid, and it can be found in many industrial fields. Dimensionless quantities, such as capillary number, Bond number and viscosity ratio between displaced fluid and displacing fluid are introduced to illustrate the effects of capillary force, viscous force, and gravity on the fluid behaviour. The surface wettability, which has an impact on the finger pattern, is also considered in the simulation. The numerical procedure is validated against the experiment about viscous fingering in a Hele-Shaw cell. The displacement efficiency is investigated using the parameter, areal sweep efficiency. The present simulation shows an additional evidence to demonstrate that the lattice Boltzmann method is a useful method for simulating some multiphase flow problems in porous media.     

基于物理界面的双重介质有限裂隙渗流模型

基于连续介质或者离散裂隙假设,含裂隙的多孔介质渗流问题有多种数学力学模型。受物理界面的启发,提出一种新的有限裂隙连续介质力学模型,可以为宏观裂隙-多孔介质内的流体输运问题等提供近似计算方案。该模型属于一类双重介质模型,将曲面上低维度的流场转化为三维空间的流场,并且与连续的多孔介质的流场耦合,在数学上表示为统一的输运控制方程和初始边界条件。这个近似模型为不方便实施高维度-低维度耦合求解的数值计算方法提供新的模拟思路,如光滑粒子流体动力学等无网格粒子类方法。     

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     

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

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