Stokesian Dynamics Simulation of Suspension Flow in Porous Media

Transport in Porous Media - Suspension flow through porous medium was studied using the Stokesian dynamics simulation method. Stokesian dynamics is an efficient tool to carry out numerical...     

Numerical Simulation of Two-Phase Inertial Flow in Heterogeneous Porous Media

In this study, non-Darcy inertial two-phase incompressible and non-stationary flow in heterogeneous porous media is analyzed using numerical simulations. For the purpose, a 3D numerical tool was fully developed using a finite volume formulation, although for clarity, results are presented in 1D and 2D configurations only. Since a formalized theoretical model confirmed by experimental data is still lacking, our study is based on the widely used generalized Darcy–Forchheimer model. First, a validation is performed by comparing numerical results of the saturation front kinetics with a semi-analytical solution inspired from the Buckley–Leverett model extended to take into account inertia. Second, we highlight the importance of inertial terms on the evolution of saturation fronts as a function of a suitable Reynolds number. Saturation fields are shown to have a structure markedly different from the classical case without inertia, especially for heterogeneous media, thereby, emphasizing the necessity of a more complete model than the classical generalized Darcy’s one when inertial effects are not negligible.     

膨胀性非牛顿流体多孔介质流动的模拟分析

在表征体元尺度采用格子Boltzmann方法分析膨胀性非牛顿流体在多孔介质中的流动,基于二阶矩模型在演化方程中引入表征介质阻力的作用力项,求解描述渗流模型的广义Navier-Stokes方程．采用局部法计算形变速率张量,通过循环迭代得到非牛顿粘度和松弛时间．对多孔介质的Poiseuille流动进行分析,通过比较发现结果与孔隙尺度的解析解十分吻合,并且收敛较快,表明方法合理有效．分析了渗透率和幂律指数对速度和压力降的影响,研究结果表明,膨胀性流体的多孔介质流动不符合达西规律,压力降的增加幅度小于渗透率的减小幅度．当无量纲渗透率Da小于10-5时,流道中的速度呈现均匀分布,并且速度分布随着幂律指数的减小趋于平滑．压力降随着幂律指数的增加而增加,Da越大幂律指数对压力降的影响越明显．     

Flow in Random Porous Media

Flow in a porous medium with a random hydraulic conductivity tensor K(x) is analyzed when the mean conductivity tensor (x) is a non-constant function of position x. The results are a non-local expression for the mean flux vector (x) in terms of the gradient of the mean hydraulic head (x), an integrodifferential equation for (x), and expressions for the two point covariance functions of q(x) and (x). When K(x) is a Gaussian random function, the joint probability distribution of the functions q(x) and (x) is determined.     

Effective Rheology of Two-Phase Flow in Three-Dimensional Porous Media: Experiment and Simulation

We present an experimental and numerical study of immiscible two-phase flow of Newtonian fluids in three-dimensional (3D) porous media to find the relationship between the volumetric flow rate (Q) and the total pressure difference (\(\Delta P\)) in the steady state. We show that in the regime where capillary forces compete with the viscous forces, the distribution of capillary barriers at the interfaces effectively creates a yield threshold (\(P_t\)), making the fluids reminiscent of a Bingham viscoplastic fluid in the porous medium. In this regime, Q depends quadratically on an excess pressure drop (\(\Delta P-P_t\)). While increasing the flow rate, there is a transition, beyond which the overall flow is Newtonian and the relationship is linear. In our experiments, we build a model porous medium using a column of glass beads transporting two fluids, deionized water and air. For the numerical study, reconstructed 3D pore networks from real core samples are considered and the transport of wetting and non-wetting fluids through the network is modeled by tracking the fluid interfaces with time. We find agreement between our numerical and experimental results. Our results match with the mean-field results reported earlier.     

Upscaled Unstructured Computational Grids for Efficient Simulation of Flow in Fractured Porous Media

Discrete fracture modeling (DFM) is currently the most promising approach for modeling of naturally fractured reservoirs and simulation of multiphase fluid flow therein. In contrast with the classical double-porosity/double permeability models, in the DFM approach all the interactions and fluid flow in and between the fractures and within the matrix are modeled in a unified manner, using the same computational grid. There is no need for computing the shape factors, which are crucial to the accuracy of the double-porosity models. We have exploited this concept in order to develop a new method for the generation of unstructured computational grids. In the new approach the geological model (GM) of the reservoir is first generated, using square or cubic grid blocks. The GM is then upscaled using a method based on the multiresolution wavelet transformations that we recently developed. The upscaled grid contains a distribution of the square or cubic blocks of various sizes. A map of the blocks’ centers is then used with an optimized Delauney triangulation method and the advancing-front technique, in order to generate the final unstructured triangulated grid suitable for use in any general reservoir simulator with any number of fluid phases. The new method also includes an algorithm for generating fractures that, contrary to the previous methods, does not require modifying their paths due to the complexities that may arise in spatial distribution of the grid blocks. It also includes an effective partitioning of the simulation domain that results in large savings in the computation times. The speed-up in the computations with the new upscaled unstructured grid is about three orders of magnitude over that for the initial GM. Simulation of waterflooding indicates that the agreement between the results obtained with the GM and the upscaled unstructured grid is excellent. The method is equally applicable to the simulations of multiphase flow in unfractured, but highly heterogeneous, reservoirs.     

Rarefied Pure Gas Transport in Non-isothermal Porous Media: Validation and Tests of the Model

Viscous flow, effusion, and thermal transpiration are the main gas transport modalities for a rarefied gas in a macro-porous medium. They have been well quantified only in the case of simple geometries. This paper presents a numerical method based on the homogenization of kinetic equations producing effective transport properties (permeability, Knudsen diffusivity, thermal transpiration ratio) in any porous medium sample, as described, e.g. by a digitized 3D image. The homogenization procedure—neglecting the effect of gas density gradients on heat transfer through the solid—leads to closure problems in for the obtention of effective properties; they are then simplified using a Galerkin method based on a 21-element basis set. The kinetic equations are then discretized in space with a finite- volume scheme. The method is validated against experimental data in the case of a closed test tube. It shows to be coherent with past approaches of thermal transpiration. Then, it is applied to several 3D images of increasing complexity. Another validation is brought by comparison with other distinct numerical approaches for the evaluation of the Darcian permeability tensor and of the Knudsen diffusion tensor. Results show that thermal transpiration has to be described by an effective transport tensor which is distinct from the other tensors.     

High-Velocity Laminar and Turbulent Flow in Porous Media

We model high-velocity flow in porous media with the multiple scale homogenization technique and basic fluid mechanics. Momentum and mechanical energy theorems are derived. In idealized porous media inviscid irrotational flow in the pores and wall boundary layers give a pressure loss with a power of 3/2 in average velocity. This model has support from flow in simple model media. In complex media the flow separates from the solid surface. Pressure loss effects of flow separation, wall and free shear layers, pressure drag, flow tube velocity and developing flow are discussed by using phenomenological arguments. We propose that the square pressure loss in the laminar Forchheimer equation is caused by development of strong localized dissipation zones around flow separation, that is, in the viscous boundary layer in triple decks. For turbulent flow, the resulting pressure loss due to average dissipation is a power 2 term in velocity.     

Nonequilibrium Effects and Multiphase Flow in Porous Media

In this paper we develop a more general formulation for transient multiphase flow in porous media based on physics observed in core-scale and micromodel experiments. We account for non-equilibrium effects by considering redistribution time and treat saturation by evolving locally moving time-averages of the saturation. Several families of models arise from approximations to the general formulation with various degrees of accuracy. The classical Buckley-Leverett and Barenblatt expressions are special cases of these families. We explore the behaviors of a number of special cases arising from the proposed general formulation using established and novel numerical schemes that provide nonlinear physics-based preconditioning. The agreement observed between numerical and experimental results demonstrates the consistency of the proposed abstraction.     

Combined Free and Forced Convection in Vertical Channels of Porous Media

In this paper we investigate the combined free and forced convection of a fully developed Newtonian fluid within a vertical channel composed of porous media when viscous dissipation effects are taken into consideration. The flow is analysed in the region of a first critical Rayleigh number in order to interpret the multiple-valued solutions and discuss their validity. The governing fourth-order, ordinary differential equation, which contains the Darcy and the viscous dissipation terms, is solved analytically using perturbation techniques and numerically using D02HBF NAG Library. A detailed investigation of the governing O.D.E. is performed on both clear fluid and porous medium for various values of the viscous dissipation parameter, , when the wall temperature decreases linearly with height, and the pressure gradient is both above and below its hydrostatic value. Although mathematically the results in all cases show that there are two solution branches, producing four possible solutions, the study of the velocity and buoyancy profiles together with the Darcy effect indicate that only one of the two solutions at any value of the Rayleigh number appears to be physically acceptable. It is shown that the effect of the Darcy number decreases as the critical Rayleigh numbers increase.     

Steady-State Source Flow in Heterogeneous Porous Media

The flow of fluids in heterogeneous porous media is modelled by regarding the hydraulic conductivity as a stationary random space function. The flow variables, the pressure head and velocity field are random functions as well and we are interested primarily in calculating their mean values. The latter had been intensively studied in the past for flows uniform in the average. It has been shown that the average Darcy's law, which relates the mean pressure head gradient to the mean velocity, is given by a local linear relationship. As a result, the mean head and velocity satisfy the local flow equations in a fictitious homogeneous medium of effective conductivity. However, recent analysis has shown that for nonuniform flows the effective Darcy's law is determined by a nonlocal relationship of a convolution type. Hence, the average flow equations for the mean head are expressed as a linear integro-differential operator. Due to the linearity of the problem, it is useful to derive the mean head distribution for a flow by a source of unit discharge. This distribution represents a fundamental solution of the average flow equations and is called the mean Green function G _d (x). The mean head G _d(x) is derived here at first order in the logconductivity variance for an arbitrary correlation function (x) and for any dimensionality d of the flow. It is obtained as a product of the solution G _d ⁽⁰⁾(x) for source flow in unbounded domain of the mean conductivity K _A and the correction _d (x) which depends on the medium heterogeneous structure. The correction _d is evaluated for a few cases of interest.Simple one-quadrature expressions of _d are derived for isotropic two- and three-dimensional media. The quadratures can be calculated analytically after specifying (x) and closed form expressions are derived for exponential and Gaussian correlations. The flow toward a source in a three-dimensional heterogeneous medium of axisymmetric anisotropy is studied in detail by deriving ₃ as function of the distance from the source x and of the azimuthal angle . Its dependence on x, on the particular (x) and on the anisotropy ratio is illustrated in the plane of isotropy (=0) and along the anisotropy axis ( = /2).The head factor k ^* is defined as a ratio of the head in the homogeneous medium to the mean head, k ^*=G _d ⁽⁰⁾/G _d= _d ^–1. It is shown that for isotropic conductivity and for any dimensionality of the flow the medium behaves as a one-dimensional and as an effective one close and far from the source, respectively, that is, lim _x0 k ^*(x) = K _H/K _A and lim _x k ^*(x) = K ^efu/K _A, where K _A and K _H are the arithmetic and harmonic conductivity means and K ^efu is the effective conductivity for uniform flow. For axisymmetric heterogeneity the far-distance limit depends on the direction. Thus, in the coordinate system of (x) principal directions the limit values of k ^* are obtained as . These values differ from the corresponding components of the effective conductivities tensor for uniform flow for = 0 and /2, respectively. The results of the study are applied to solving the problem of the dipole well flow. The dependence of the mean head drop between the injection and production chambers on the anisotropy of the conductivity and the distance between the chambers is analyzed.     

Macroscopic Flow Potentials in Swelling Porous Media

In swelling porous media, the potential for flow is much more than pressure, and derivations for flow equations have yielded a variety of equations. In this article, we show that the macroscopic flow potentials are the electro-chemical potentials of the components of the fluid and that other forms of flow equations, such as those derived through mixture theory or homogenization, are a result of particular forms of the chemical potentials of the species. It is also shown that depending upon whether one is considering the pressure of a liquid in a reservoir in electro-chemical equilibrium with the swelling porous media, or the pressure of the vicinal liquid within the swelling porous media, a critical pressure gradient threshold exists or does not.     

Modelling Foamy Oil Flow in Porous Media

This paper describes a dynamic model for the simulation of foamy oil flow in porous media. The model includes expressions for the rate processes of nucleation, bubble growth and disengagement of dispersed gas bubbles from the oil. The model is used to simulate experimental results pertaining to primary depletion tests conducted in a sand pack. Using the model to interpret experimental results indicated that, although the lifetimes of supersaturation and dispersed gas bubbles may be short, supersaturated conditions are likely to exist, and dispersed gas bubbles are likely to be present during the entire production period, as long as the pressure continues to decline at a high rate. The model developed in this paper gave better agreement with experimental data than other proposed models. The effect of foamy oil flow increases as the rate of pressure decline increases.     

An Orthorhombic Lattice Boltzmann Model for Pore-Scale Simulation of Fluid Flow in Porous Media

One application of the lattice Boltzmann equation (LBE) models is in combination with tomography to simulate pore-scale flow and transport processes in porous media. Most LBE models in the literature are based on cubic lattice, and if the voxels in a tomography image are not cubic or cannot be divided into cubes due to computational limitations, these models will lose most of their advantages. How to deal with such images is, hence, an interest in use of the LBE model to simulate pore-scale processes. In this paper, we present an orthorhombic LBE model based on the single-relaxation time approach with the relaxation parameter varying with lattice directions. The equilibrium distribution functions in the standard LBE model were modified to correct the anisotropy induced by the non-cubic lattice, and the calculations of the fluid density and momentum were also redefined in order to maintain the conservation of mass and momentum during the collision. We tested the model against analytical solution for fluid flow in a tube, and against the standard cubic-based LBE model for fluid flow in a duct with an island inside. The model was then applied to simulate fluid flow in a 3D image in attempts to analyse the errors if the voxels in the image are not cubic but are assumed to be cubic.     

Direct Numerical Simulation of Pore-Scale Trapping Events During Capillary-Dominated Two-Phase Flow in Porous Media

Transport in Porous Media - This study focuses on direct numerical simulation of imbibition, displacement of the non-wetting phase by the wetting phase, through water-wet carbonate rocks. We...     

Analysis of Multiphase Non-Darcy Flow in Porous Media

Recent laboratory studies and analyses (Lai et al. Presented at the 2009 Rocky Mountain Petroleum Technology Conference, 14–16 April, Denver, CO, 2009) have shown that the Barree and Conway model is able to describe the entire range of relationships between flow rate and potential gradient from low- to high-flow rates through porous media. A Buckley and Leverett type analytical solution is derived for non-Darcy displacement of immiscible fluids in porous media, in which non-Darcy flow is described using the Barree and Conway model. The comparison between Forchheimer and Barree and Conway non-Darcy models is discussed. We also present a general mathematical and numerical model for incorporating the Barree and Conway model in a general reservoir simulator to simulate multiphase non-Darcy flow in porous media. As an application example, we use the analytical solution to verify the numerical solution for and to obtain some insight into one-dimensional non-Darcy displacement of two immiscible fluids with the Barree and Conway model. The results show how non-Darcy displacement is controlled not only by relative permeability, but also by non-Darcy coefficients, characteristic length, and injection rates. Overall, this study provides an analysis approach for modeling multiphase non-Darcy flow in reservoirs according to the Barree and Conway model.     

Investigation of Post-Darcy Flow in Thin Porous Media

Transport in Porous Media - We present numerical simulations of post-Darcy flow in thin porous medium: one consisting of staggered arrangements of circular cylinders and one random distribution of...     

A Criterion for Non-Darcy Flow in Porous Media

Non-Darcy behavior is important for describing fluid flow in porous media in situations where high velocity occurs. A criterion to identify the beginning of non-Darcy flow is needed. Two types of criteria, the Reynolds number and the Forchheimer number, have been used in the past for identifying the beginning of non-Darcy flow. Because each of these criteria has different versions of definitions, consistent results cannot be achieved. Based on a review of previous work, the Forchheimer number is revised and recommended here as a criterion for identifying non-Darcy flow in porous media. Physically, this revised Forchheimer number has the advantage of clear meaning and wide applicability. It equals the ratio of pressure drop caused by liquid–solid interactions to that by viscous resistance. It is directly related to the non-Darcy effect. Forchheimer numbers are experimentally determined for nitrogen flow in Dakota sandstone, Indiana limestone and Berea sandstone at flowrates varying four orders of magnitude. These results indicate that superficial velocity in the rocks increases non-linearly with the Forchheimer number. The critical Forchheimer number for non-Darcy flow is expressed in terms of the critical non-Darcy effect. Considering a 10% non-Darcy effect, the critical Forchheimer number would be 0.11.