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.     

湿润性对孔隙介质两相渗流驱替效率的影响

孔隙介质中多相渗流的驱替效率对二氧化碳封存效率和石油采收率具有决定性影响, 是实际工程调控中的一个关键指标. 湿润性是影响多相渗流驱替模式及其效率的一个重要因素. 本文通过微流体模型-显微镜-高速相机可视化实验平台, 对基于真实砂岩孔隙结构的微流体模型进行湿润性修饰, 开展了5种流量和2种湿润性的两相驱替可视化实验, 研究了湿润性对砂岩孔隙结构中两相渗流驱替模式及其效率的重要影响. 实验结果表明: 随着流速的增大, 两相渗流驱替模式由毛细指流向稳定流发生转变; 在低流速条件下, 由于毛细力的主导效应, 亲水性介质中指进的宽度和被驱替流体团簇的数目均小于疏水性介质, 而被驱替流体团簇的最大半径、平均半径和方差均大于疏水性介质. 实验结果还证实了亲水性介质中由于单支优势通道和"绕流"现象的发生, 驱替效率显著小于疏水性介质. 最后, 通过考虑接触角效应对毛细管数进行修正, 建立了考虑湿润性影响的驱替效率和毛细管数之间的统一关系式, 为不同湿润性条件下驱替效率的预测提供了一种潜在方法.      

Pore-Scale Network Modeling of Three-Phase Flow Based on Thermodynamically Consistent Threshold Capillary Pressures. I. Cusp Formation and Collapse

We present a pore-scale network model of two- and three-phase flow in disordered porous media. The model reads three-dimensional pore networks representing the pore space in different porous materials. It simulates wide range of two- and three-phase pore-scale displacements in porous media with mixed-wet wettability. The networks are composed of pores and throats with circular and angular cross sections. The model allows the presence of multiple phases in each angular pore. It uses Helmholtz free energy balance and Mayer–Stowe–Princen (MSP) method to compute threshold capillary pressures for two- and three-phase displacements (fluid configuration changes) based on pore wettability, pore geometry, interfacial tension, and initial pore fluid occupancy. In particular, it generates thermodynamically consistent threshold capillary pressures for wetting and spreading fluid layers resulting from different displacement events. Threshold capillary pressure equations are presented for various possible fluid configuration changes. By solving the equations for the most favorable displacements, we show how threshold capillary pressures and final fluid configurations may vary with wettability, shape factor, and the maximum capillary pressure reached during preceding displacement processes. A new cusp pore fluid configuration is introduced to handle the connectivity of the intermediate wetting phase at low saturations and to improve model’s predictive capabilities. Based on energy balance and geometric equations, we show that, for instance, a gas-to-oil piston-like displacement in an angular pore can result in a pore fluid configuration with no oil, with oil layers, or with oil cusps. Oil layers can then collapse to form cusps. Cusps can shrink and disappear leaving no oil behind. Different displacement mechanisms for layer and cusp formation and collapse based on the MSP analysis are implemented in the model. We introduce four different layer collapse rules. A selected collapse rule may generate different corner configuration depending on fluid occupancies of the neighboring elements and capillary pressures. A new methodology based on the MSP method is introduced to handle newly created gas/water interfaces that eliminates inconsistencies in relation between capillary pressures and pore fluid occupancies. Minimization of Helmholtz free energy for each relevant displacement enables the model to accurately determine the most favorable displacement, and hence, improve its predictive capabilities for relative permeabilities, capillary pressures, and residual saturations. The results indicate that absence of oil cusps and the previously used geometric criterion for the collapse of oil layers could yield lower residual oil saturations than the experimentally measured values in two- and three-phase systems.     

海洋地震工程流固耦合问题统一计算框架

海底地震动的模拟以及海洋工程结构的地震反应分析中,涉及到海水、饱和海床、弹性基岩、结构之间的相互耦合.传统的方法分别采用声波方程描述理想流体、Biot方程描述饱和海床、弹性波方程描述基岩和结构,分别考虑相互之间的耦合,十分不便.本文基于理想流体、固体分别为饱和多孔介质的特殊情形(孔隙率分别为1和0),由饱和多孔介质的Biot方程可退化得到理想流体的声波方程和固体的弹性波方程.然后,以饱和多孔介质方程为基础,经集中质量有限元离散,考虑不同孔隙率的饱和多孔介质之间耦合的一般情形,建立了该耦合情形的求解方法.进一步论证了该一般情形的耦合计算方法可分别退化到流体与固体、流体与饱和多孔介质、固体与饱和多孔介质之间的耦合计算,从而将流体、固体、饱和多孔介质间的耦合问题纳入到统一计算框架,并编制了相应的三维并行分析程序.以P-SV波垂直入射时,半无限层状海水-饱和海床、海水-弹性基岩、海水-饱和海床-弹性基岩三种情形的动力分析为例,采用统一计算框架结合透射边界条件进行求解,并与传递矩阵方法得到的解进行对比,验证了该统一计算框架的有效性以及并行计算的可行性.      

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.     

Flow behaviour in the presence of wettability heterogeneities

The physical processes occurring during fluid flow and displacement within porous media having wettability heterogeneities have been investigated in specially designed heterogeneous visual models. The models were packed with glass beads, areas of which were treated with a water repellent to create wettability variations. Immiscible displacement experiments show visually the effect of wettability heterogeneities on the formation of residual oil and recovery due to capillary trapping. This work demonstrates by experiment the importance of incorporating reservoir heterogeneity into pore displacement analysis, essential for the correct interpretation of core data and for directing the route for scale-up of the processes to reservoir scale.     

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     

Pressure-driven and electroosmotic non-Newtonian flows through microporous media via lattice Boltzmann method

The lattice Boltzmann method is developed to simulate the pressure-driven flow and electroosmotic flow of non-Newtonian fluids in porous media based on the representative elementary volume scale. The flow through porous media was simulated by including the porosity into the equilibrium distribution function and adding a non-Newtonian force term to the evolution equation. The non-Newtonian behavior is considered based on the Herschel–Bulkley model. The velocity results for pressure-driven non-Newtonian flow agree well with the analytical solutions. For the electroosmotic flow, the influences of porosity, solid particle diameter, power law exponent, yield stress and electric parameters are investigated. The results demonstrate that the present lattice Boltzmann model is capable of modeling non-Newtonian flow through porous media.     

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

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

A computational framework for fluid–porous structure interaction with large structural deformation

We study the effect of poroelasticity on fluid–structure interaction. More precisely, we analyze the role of fluid flow through a deformable porous matrix in the energy dissipation behavior of a poroelastic structure. For this purpose, we develop and use a nonlinear poroelastic computational model and apply it to the fluid–structure interaction simulations. We discretize the problem by means of the finite element method for the spatial approximation and using finite differences in time. The numerical discretization leads to a system of non-linear equations that are solved by Newton’s method. We adopt a moving mesh algorithm, based on the Arbitrary Lagrangian–Eulerian method to handle large deformations of the structure. To reduce the computational cost, the coupled problem of free fluid, porous media flow and solid mechanics is split among its components and solved using a partitioned approach. Numerical results show that the flow through the porous matrix is responsible for generating a hysteresis loop in the stress versus displacement diagrams of the poroelastic structure. The sensitivity of this effect with respect to the parameters of the problem is also analyzed.

     

Nonlinear fluid flow laws for orthotropic porous media

Nonlinear fluid flow laws for orthotropic porous media are written in invariant tensor form. As usual in the theory of fluid flow through porous media [1, 2], the equations contain the flow velocity up to the second power. Expressions that determine the nonlinear resistances to fluid flow are presented and it is shown that, on going over from linear to nonlinear flow laws, the asymmetry effect may manifest itself, that is, the fluid flow characteristics may differ along the same straight line in the positive and negative directions. It is shown that, as compared with the linear fluid flow law for orthotropic media when for three symmetry groups a single flow law is sufficient, in nonlinear laws the anisotropy manifestations are much more variable and each symmetry group must be described by specific equations. A system of laboratory measurements for finding the nonlinear flow characteristics for orthotropic porous media is considered.     

海洋地震工程流固耦合问题统一计算框架——————不规则界面情形

海底地震动场及海洋声场的模拟中,需要考虑复杂海床介质及海底地形的影响,涉及到海水、饱和海床、弹性基岩之间的相互耦合.传统的方法分别采用声波方程描述理想流体、Biot方程描述饱和海床、弹性波方程描述基岩,分别进行空间离散和界面耦合, 十分不便.本文基于理想流体、固体分别为饱和多孔介质的特殊情形(孔隙率分别为1和0),由饱和多孔介质的Biot方程可退化得到理想流体的声波方程和固体的弹性波方程.然后, 以饱和多孔介质方程为基础, 经集中质量有限元离散,严格考虑不同孔隙率的饱和多孔介质在不规则界面的耦合条件,通过求解法向和切向界面力的途径,建立了不同孔隙率的饱和多孔介质耦合情形的求解方法,将流体、固体、饱和多孔介质间的耦合问题纳入到统一计算框架,并编制了相应的三维并行分析程序.考虑海水--弹性基岩、海水--饱和海床--弹性基岩体系中凹陷地形情形,采用本文提出的统一计算框架, 结合透射边界条件,分析了P波入射时的动力反应, 并通过结果是否满足界面条件,验证了该统一计算框架的有效性以及并行计算的可行性.      

The Physical Origin of Interfacial Coupling in Two-Phase Flow through Porous Media

Recently developed transport equations for two-phase flow through porous media usually have a second term that has been included to account properly for interfacial coupling between the two flowing phases. The source and magnitude of such coupling is not well understood. In this study, a partition concept has been introduced into Kalaydjian's transport equations to construct modified transport equations that enable a better understanding of the role of interfacial coupling in two-phase flow through natural porous media. Using these equations, it is demonstrated that, in natural porous media, the physical origin of interfacial coupling is the capillarity of the porous medium, and not interfacial momentum transfer, as is usually assumed. The new equations are also used to show that, under conditions of steady-state flow, the magnitude of mobilities measured in a countercurrent flow experiment is the same as that measured in a cocurrent flow experiment, contrary to what has been reported previously. Moreover, the new equations are used to explicate the mechanism by which a saturation front steepens in an unstabilized displacement, and to show that the rate at which a wetting fluid is imbibed into a porous medium is controlled by the capillary coupling parameter, . Finally, it is argued that the capillary coupling parameter, , is dependent, at least in part, on porosity. Because a clear understanding of the role played by interfacial coupling is important to an improved understanding of two-phase flow through porous media, the new transport equations should prove to be effective tools for the study of such flow.     

Contribution to the Modelling of Coupled Heat and Mass Transfers on 3D Real Structure of Heterogeneous Building Materials: Application to Hemp Concrete

Hydraulic flow, electrical flow and the passage of elastic waves through porous media are all linked by electrokinetic processes. In its simplest form, the passage of elastic waves through the porous medium causes fluid to flow through that medium and that flow gives rise to an electrical streaming potential and electrical counter-current. These processes are frequency-dependent and governed by coupling coefficients which are themselves frequency-dependent. The link between fluid pressure and fluid flow is described by dynamic permeability, which is characterised by the hydraulic coupling coefficient (C_hp). The link between fluid pressure and electrical streaming potential is characterised by the streaming potential coefficient (C_sp). While the steady-state values of such coefficients are well studied and understood, their frequency dependence is not. Previous work has been confined to unconsolidated and disaggregated materials such as sands, gravels and soils. In this work, we present an apparatus for measuring the hydraulic and streaming potential coefficients of high porosity, high permeability consolidated porous media as a function of frequency. The apparatus operates in the range 1 Hz to 2 kHz with a sample of 10 mm diameter and 5–30 mm in length. The full design and validation of the apparatus are described together with the experimental protocol it uses. Initial data are presented for three samples of Boise sandstone, which present as dispersive media with the critical transition frequency of 918.3?±?99.4 Hz. The in-phase and in-quadrature components of the measured hydraulic and streaming potential coefficients have been compared to the Debye-type dispersion model as well as theoretical models based on bundles of capillary tubes and porous media. Initial results indicate that the dynamic permeability data present an extremely good fit to the capillary bundle and Debye-type dispersion models, while the streaming potential coefficient presents an extremely good fit to all of the models up to the critical transition frequency, but diverges at higher frequencies. The streaming potential coefficient data are best fitted by the Pride model and its Walker and Glover simplification. Characteristic pore size values calculated from the measured critical transition frequency fell within 1.73% of independent measures of this parameter, while the values calculated directly from the Packard model showed an underestimation by about 12%.

     

Dependency of Tortuosity and Permeability of Porous Media on Directional Distribution of Pore Voids

Single-phase fluid flow in porous media is usually direction dependent owing to the tortuosity associated with the internal structures of materials that exhibit inherent anisotropy. This article presents an approach to determine the tortuosity and permeability of porous materials using a structural measure quantifying the anisotropic distribution of pore voids. The approach uses a volume averaging method through which the macroscopic tortuosity tensor is related to both the average porosity and the directional distribution of pore spaces. The permeability tensor is derived from the macroscopic momentum balance equation of fluid in a porous medium and expressed as a function of the tortuosity tensor and the internal structure of the material. The analytical results generally agree with experimental data in the literature.     

Modeling Geometric State for Fluids in Porous Media: Evolution of the Euler Characteristic

Multiphase flow in porous media is strongly influenced by the pore-scale arrangement of fluids. Reservoir-scale constitutive relationships capture these effects in a phenomenological way, relying only on fluid saturation to characterize the macroscopic behavior. Working toward a more rigorous framework, we make use of the fact that the momentary state of such a system is uniquely characterized by the geometry of the pore-scale fluid distribution. We consider how fluids evolve as they undergo topological changes induced by pore-scale displacement events. Changes to the topology of an object are fundamentally discrete events. We describe how discontinuities arise, characterize the possible topological transformations and analyze the associated source terms based on geometric evolution equations. Geometric evolution is shown to be hierarchical in nature, with a topological source term that constrains how a structure can evolve with time. The challenge associated with predicting topological changes is addressed by constructing a universal geometric state function that predicts the possible states based on a non-dimensional relationship with two degrees of freedom. The approach is validated using fluid configurations from both capillary and viscous regimes in ten different porous media with porosity between 0.10 and 0.38. We show that the non-dimensional relationship is independent of both the material type and flow regime. We demonstrate that the state function can be used to predict history-dependent behavior associated with the evolution of the Euler characteristic during two-fluid flow.

     

Smoothed Particle Hydrodynamics Modeling of Transverse Flow in Randomly Aligned Fibrous Porous Media

The Lagrangian smoothed particle hydrodynamics (SPH) method is employed to obtain a meso-/micro-scopic pore-scale insight into the transverse flow across the randomly aligned fibrous porous media in a 2D domain. Fluid is driven by an external body force, and a square domain with periodic boundary conditions imposed at both the streamwise and transverse flow direction is assumed. The porous matrix is established by randomly embedding a certain number of fibers in the square domain. Fibers are represented by position-fixed SPH particles, which exert viscous forces upon, and contribute to the density variations of, the nearby fluid particles. An additional repulsive force, similar in form to the 12-6 Lennard-Jones potential between atoms, is introduced to consider the no-penetrating restraint prescribed by the solid pore structure. This force is initiated from the fixed solid material particle and may act on its neighboring moving fluid particles. Fluid flow is visualized by plotting the local velocity vector field; the meandering fluid flow around the porous microstructures always follow the paths of least resistance. The simulated steady-state flow field is further used to calculate the macroscopic permeability. The dimensionless permeability (normalized by the squared characteristic dimension of the fiber cross section) exhibits an exponential dependence on the porosity within the intermediate porosity range, and the derived dimensionless permeability—porosity relation is found to have only minor dependence on either the relative arrangement condition among fibers or the fiber cross section (shape or area).     

Thermodynamically Consistent Limiting Forced Convection Heat Transfer in a Asymmetrically Heated Porous Channel: An Analytical Study

Fluid transport and the associated heat transfer through porous media is of immense importance because of its numerous practical applications. In view of the widespread applications of porous media flow, the present study attempts to investigate the forced convective heat transfer in the limiting condition for the flow through porous channel. There could be many areas, where heat transfer through porous channel attain some limiting conditions, thus, the analysis of limiting convective heat transfer is far reaching. The primary aim of the present study is focused on the limiting forced convection analysis considering the flow of Newtonian fluid between two asymmetrically heated parallel plates filled with saturated porous media. Utilizing a few assumptions, which are usually employed in the literature, an analytical methodology is executed to obtain the closed-form expression of the temperature profile, and in the following the expression of the limiting Nusselt numbers. The parametric variations of the temperature profile and the Nusselt numbers in different cases have been shown highlighting the influential role of different performance indexing parameters, like Darcy number, porosity of the media, and Brinkman number of forced convective heat transfer in porous channel. In doing so, the underlying physics of the transport characteristics of heat has been delineated in a comprehensive way. Moreover, a discussion has been made regarding an important feature like the onset of point of singularity as appeared on the variation of the Nusselt number from the consideration of energy balance in the flow field, and in view of second law of thermodynamics.     

Effect of Organosilicon-Acrylic Emulsion Treatment on Wettability of Porous Media

To investigate the influence of the organosilicon-acrylic on wetting properties of porous media, contact angle tests were performed on two different sandstones. In addition, the effectiveness of the emulsion on wettability alteration of porous media was validated by capillary rise and spontaneous imbibition tests. The results of wettability tests showed that the wettability of two sandstones was altered from water-wet to gas-wet after treatment with the emulsion. The principle that the critical radius of pore throats and wettability of porous media affect liquids flow was derived analytically and verified experimentally. Coreflood results demonstrated that the latex resulted in increasing the water permeability through altering the rock wettability to gas-wetting, then decreasing the friction drag between liquids and rocks surface. Thereby, the emulsion treatment could increase the flowback rate of trapped liquids. Experimental results were in good agreement with the theoretical analysis. In conclusion, all results indicated that the emulsion could alter the wettability from water-wet to intermediate gas-wet and enhance water permeability in porous media. It was extrapolated that the emulsion had the tremendous potential to be applied in field conditions, enhancing gas productivity through the cleanup of trapped water in the vicinity of the wellbore.