A Dynamic Network Model for Two-Phase Flow in Porous Media

We present a dynamic model of immiscible two-phase flow in a network representation of a porous medium. The model is based on the governing equations describing two-phase flow in porous media, and can handle both drainage, imbibition, and steady-state displacement. Dynamic wetting layers in corners of the pore space are incorporated, with focus on modeling resistivity measurements on saturated rocks at different capillary numbers. The flow simulations are performed on a realistic network of a sandpack which is perfectly water-wet. Our numerical results show saturation profiles for imbibition in agreement with experiments. For free spontaneous imbibition we find that the imbibition rate follows the Washburn relation, i.e., the water saturation increases proportionally to the square root of time. We also reproduce rate effects in the resistivity index for drainage and imbibition.     

A Network Model for Diffusion in Media with Partially Resolvable Pore Space Characteristics

A microstructure-informed meso-scale model for diffusion of foreign species in porous media is proposed. The model is intended for media where the pore geometry data acquired experimentally represent a fraction of total porosity. A cellular complex, with a cell representing the average pore neighbourhood, is used to generate 3D graphs of sites at cell centres and bonds between neighbouring cells. The novel interpretation of pore systems as graphs allows for clear separation between topology (here connectedness) and physics (here diffusion) in the mathematical formulation of transport. Further, it allows for easy introduction of dynamics into the system, i.e., local changes in topology due to other physical mechanisms, such as micro-cracking or blockage of pores. A mapping between microstructure features and graph elements is used for model construction. The mapping is based on data for clays, where the experimentally resolved pore system comprises isolated elongated pores of preferred orientation with a large volume fraction of unresolved pores. Both the resolved and the “hidden” systems are accounted for. The graph geometry is described by a principal length, the cell size in the preferred orientation, and a secondary length, the cell size out of preferred orientation. This is considered as a representation of mineralogical heterogeneity of clays. Analysis on graphs, a specialisation of the discrete exterior calculus, is used to obtain connectivity and diffusivity properties of formed networks. Since the experimental data are not sufficient to determine the principal length, upper and lower limits are determined from the limited information. Effects of the principal cell size between limits and of the secondary cell size are studied. The results are within the range of experimentally measured macroscopic (bulk) diffusivity for the material studied, including anisotropic diffusion coefficients. The variation of calculated diffusivity coefficients with principal and secondary lengths provides an explanation for the variability in experimentally measured coefficients across different clays.     

Nonlinear Coupled Mathematical Model for Solid Deformation and Gas Seepage in Fractured Media

Based on detailed investigation into the interactional physical mechanism of solid deformations and gas seepage in rock matrix and fracture, a nonlinear coupled mathematical model of solid deformation and gas seepage is put forward and the FEM model is built up to carry out numerical analysis. The coupled interaction laws between solid deformations and gas seepage in rock matrix and fractures has been emphasized in the model, which is a vital progress for coupled mathematical model of solid deformation and gas seepage of rock mass media. As an example, the methane extraction in fractured coal seam has been numerically simulated. By analyzing the simulation results, the law of methane migration and exchange in rock matrix and fractures is interpreted.     

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.     

Finite-Volume Discretisations for Flow in Fractured Porous Media

Over the last decades, finite-volume discretisations for flow in porous media have been extended to handle situations where fractures dominate the flow. Successful discretisations have been based on the discrete fracture-matrix models to yield mass conservative methods capable of explicitly incorporating the impact of fractures and their geometry. When combined with a hybrid-dimensional formulation, two central concerns are the restrictions arising from small cell sizes at fracture intersections and the coupling between fractures and matrix. Focusing on these aspects, we demonstrate how finite-volume methods can be efficiently extended to handle fractures, providing generalisations of previous work. We address the finite-volume methods applying a general hierarchical formulation, facilitating implementation with extensive code reuse and providing a natural framework for coupling of different subdomains. Furthermore, we demonstrate how a Schur complement technique may be used to obtain a robust and versatile method for fracture intersection cell elimination. We investigate the accuracy of the proposed elimination method through a series of numerical simulations in 3D and 2D. The simulations, performed on fractured domains containing permeability heterogeneity and anisotropy, also demonstrate the flexibility of the hierarchical framework.     

A (Dual) Network Model for Heat Transfer in Porous Media

Transport in Porous Media - We present a dual network model to simulate coupled single-phase flow and energy transport in porous media including conditions under which local thermal equilibrium...     

A Lattice Gas Automaton Simulation of the Nonlinear Diffusion Equation: A Model for Moisture Flow in Unsaturated Porous Media

We investigate a two-dimensional lattice gas automaton (LGA) for simulating the nonlinear diffusion equation in a random heterogeneous structure. The utilility of the LGA for computation of nonlinear diffusion arises from the fact that, the diffusion coefficient in the LGA depends on the local density of fluid particles which statistically determines the collision rate and thus, the mean free path of the particles at the microscopic scale. The LGA may therefore be used as a physical analogue to simulate moisture flow in unsaturated porous media. The capability of the LGA to account for unsaturated flow is tested through a set of numerical experiments simulating one-dimensional infiltration in a simplified semi-infinite homogenous isotropic porous material. Different mechanisms of interactions are used between the fluid and the solid phase to simulate various fluid–solid interfaces. The heterogeneous medium, initially at low density is submitted to a steep density gradient by continuously injecting fluid particles at high concentration and zero velocity along one face of the model. The propagation of the infiltration front is visualized at different time steps through concentration profiles parallel to the applied concentration gradient and the infiltration rate is measured continuously until steady-state flow is reached. The numerical results show close agreement with the classical theory of flow in unsaturated porous media. The cumulative absorption exhibits the expected t ^1/2 dependence. The evolution of the effective diffusion coefficient with the particle concentration is estimated from the measured density profiles for the various porous materials. Depending on the applied fluid–solid interactions, the macroscopic effective diffusivity may vary by more than two orders of magnitude with density.     

Network Model of Flow,Transport and Biofilm Effects in Porous Media

In this paper, we develop a network model to determine porosity and permeability changes in a porous medium as a result of changes in the amount of biomass. The biomass is in the form of biofilms. Biofilms form when certain types of bacteria reproduce, bond to surfaces, and produce extracellular polymer (EPS) filaments that link together the bacteria. The pore spaces are modeled as a system of interconnected pipes in two and three dimensions. The radii of the pipes are given by a lognormal probability distribution. Volumetric flow rates through each of the pipes, and through the medium, are determined by solving a linear system of equations, with a symmetric and positive definite matrix. Transport through the medium is modeled by upwind, explicit finite difference approximations in the individual pipes. Methods for handling the boundary conditions between pipes and for visualizing the results of numerical simulations are developed. Increases in biomass, as a result of transport and reaction, decrease the pipe radii, which decreases the permeability of the medium. Relationships between biomass accumulation and permeability and porosity reduction are presented.     

Comparison of Three FE-FV Numerical Schemes for Single- and Two-Phase Flow Simulation of Fractured Porous Media

We benchmark a family of hybrid finite element–node-centered finite volume discretization methods (FEFV) for single- and two-phase flow/transport through porous media with discrete fracture representations. Special emphasis is placed on a new method we call DFEFVM in which the mesh is split along fracture–matrix interfaces so that discontinuities in concentration or saturation can evolve rather than being suppressed by nodal averaging of these variables. The main objective is to illustrate differences among three discretization schemes suitable for discrete fracture modeling: (a) FEFVM with volumetric finite elements for both fractures and porous rock matrix, (b) FEFVM with lower dimensional finite elements for fractures and volumetric finite elements for the matrix, and (c) DFEFVM with a mesh that is split along material discontinuities. Fracture discontinuities strongly influence single- and multi-phase fluid flow. Continuum methods, when used to model transport across such interfaces, smear out concentration/saturation. We show that the new DFEFVM addresses this problem producing significantly more accurate results. Sealed and open single fractures as well as a realistic fracture geometry are used to conduct tracer and water-flooding numerical experiments. The benchmarking results also reveal the limitations/mesh refinement requirements of FE node-centered FV hybrid methods. We show that the DFEFVM method produces more accurate results even for much coarser meshes.     

A New Mathematical Model for Force Gravity Drainage in Fractured Porous Media

In force gas/oil gravity drainage process in fractured porous media, gas is flowing in both matrix and fractures leading to produce a finite gas pressure gradient. Consequently, viscous force plays an important role for displacing matrix oil toward fractures in addition to gravity force that is required to be modeled appropriately. A new analytical model for estimation of steady state oil saturation distribution with assumption of fixed gas pressure gradient throughout the matrix is presented. Moreover, based on some results of this analytical model a different numerical formulation is developed to predict the performance of oil production process. Comparison of the results obtained from this numerical model with the results of a conventional simulator demonstrates that the newly developed model can be applied with satisfactory accuracy. Numerical simulations show that the viscous displacement in fractured porous media can reduce the capillary threshold height, and thus it suggests the force gravity drainage as a favorable production mechanism when the matrix length is close to the threshold height.     

Macroscopic Two-phase Flow in Porous Media Assuming the Diffuse-interface Model at Pore Level

The paper presents a model for two-phase flow, where liquid and gas are treated as one fluid with variable density. A one-component fluid and the diffuse-interface model for two-phase flow are assumed at pore level. The wetting properties of the fluid are described by the Cahn theory. Macroscopic equations are deduced in the framework of the Marle formalism. It is shown that two-phase flow in porous media can be described by the Cahn–Hilliard equation for the mass density. The concept of relative permeability is not needed. For non-neutral wetting, it is shown that a capillary pressure exists but that it is not a function of state. Two numerical illustrations are presented, one of them showing that the model is, at least in a simple steady-state situation, compatible with the generalized two-continuum model.     

A Spatially Non-Local Model for Flow in Porous Media

A general mathematical model of steady-state transport driven by spatially non-local driving potential differences is proposed. The porous medium is considered to be a network of short-, medium-, and long-range interstitial channels with impermeable walls and at a continuum of length scales, and the flow rate in each channel is assumed to be linear with respect to the pressure difference between its ends. The flow rate in the model is thus a functional of the non-local driving pressure differences. As special cases, the model reduces to familiar forms of transport equations that are commonly used. An important situation arises when the phenomenon is almost, but not quite, locally dependent. The one-dimensional form of the model discussed here can be extended to multiple dimensions, temporal non-locality, and to heat, mass, and momentum transfer.     

A Subgrid-Scale Model for Turbulent Flow in Porous Media

Transport in Porous Media - Given the analogy between the filtered equations of large eddy simulation and volume-averaged Navier–Stokes equations in porous media, a subgrid-scale model is...     

Oil Displacement for One-Dimensional Three-Phase Flow with Phase Change in Fractured Media

In this article, the numerical simulations for one-dimensional three-phase flows in fractured porous media are implemented. The simulation results show that oil displacement in matrix is dominated by oil–water capillary pressure only under certain conditions. When conditions are changed to decrease the amount of water entering into the fractured media from the boundary of the flow field, water in fracture may be vaporized to superheated steam. In these cases, the appearance of superheated steam in fracture rather than in matrix will decrease the fracture pressure and generate the pressure difference between matrix and fracture, which results in oil flowing from matrix to fracture. Assuming that oil is wetting to steam, the matrix steam–oil capillary pressure will decrease the matrix oil-phase pressure as the matrix steam saturation increases. After the steam–oil capillary pressure finally exceeds the pressure difference due to the appearance of superheated steam in fracture, the oil displacement in matrix will stop. It is also shown that variations of the water relative permeability curve in matrix do not result in different mechanisms for oil displacement in matrix. The simulation results suggest that the amount of liquid water supply from the boundary of flow field fundamentally influence the mechanisms for oil displacement in matrix.     

Flow in Fractured Porous Media: A Review of Conceptual Models and Discretization Approaches

Transport in Porous Media - The last decade has seen a strong increase of research into flows in fractured porous media, mainly related to subsurface processes but also in materials science and...     

A Two-Dimensional Network Simulator for Two-Phase Flow in Porous Media

We investigate a two-dimensional network simulator that model the dynamics of drainage dominated flow where film flow can be neglected. We present a new method for simulating the temporal evolution of the pressure due to capillary and viscous forces in the displacement process. To model the dynamics, we let the local capillary pressure change as if the menisci move in and out of hour-glass shaped tubes. Furthermore, a method has been developed to allow simultaneous flow of two liquids into one tube. The model is suitable to simulate different time dependencies in two-phase drainage displacements. In this paper, we simulate the temporal evolution of the fluid pressures and analyze the time dependence of the front between the two liquids. The front width was found to be consistent with a scaling relation w t h(t/t_s). The dynamical exponent, , describing the front width evolution as function of time, was estimated to = 1.0. The results are compared to experimental data of Frette and co-workers.     

Review of Steady-State Two-Phase Flow in Porous Media: Independent Variables,Universal Energy Efficiency Map,Critical Flow Conditions,Effective Characterization of Flow and Pore Network

In many applications of two-phase flow in porous media, a wetting phase is used to displace through a network of pore conduits as much as possible of a non-wetting phase, residing in situ. The energy efficiency of this physical process may be assessed by the ratio of the flow rate of the non-wetting phase over the total mechanical power externally provided and irreversibly dissipated within the process. Fractional flow analysis, extensive simulations implementing the DeProF mechanistic model, as well as a recent retrospective examination of laboratory studies have revealed universal systematic trends of the energy efficiency in terms of the actual independent variables of the process, namely the capillary number, Ca, and the flow rate ratio, r. These trends can be cast into an energy efficiency map over the (Ca, r) domain of independent variables. The map is universal for all types of non-wetting/wetting phase porous medium systems. It demarcates the efficiency of steady-state two-phase flow processes in terms of pertinent system parameters. The map can be used as a tool for designing more efficient processes, as well as for the normative characterization of two-phase flows, as to the predominance of capillary or viscous effects. This concept is based on the existence of a unique locus of critical flow conditions, for which the energy efficiency takes locally maximum values. The locus shape depends on the physicochemical characteristics of the non-wetting phase/wetting phase/porous medium system, and it shows a significant mutation as the externally imposed flow conditions change the type of flow, from capillary- to viscosity-dominated. The locus can be approached by an S-type functional form in terms of the capillary number and the system properties (viscosity ratio, wettability, pore network geometry, etc.), suggesting that formative criteria can be derived for flow characterization in any system. A new, extended definition of the capillary number is also proposed that effectively takes into account the critical properties of all the system constituents. When loci of critical flow conditions pertaining to processes with different viscosity ratio in the same pore network, are expressed in terms of this true-to-mechanism capillary number, they collapse into a unique locus. In this context, a new methodology for the effective characterization of pore networks is proposed.     

A Model for Flow and Deformation in Unsaturated Swelling Porous Media

A thermomechanical theory for multiphase transport in unsaturated swelling porous media is developed on the basis of Hybrid Mixture Theory (saturated systems can also be modeled as a special case of this general theory). The aim is to comprehensively and non-empirically describe the effect of viscoelastic deformation on fluid transport (and vice versa) for swelling porous materials. Three phases are considered in the system: the swelling solid matrix s, liquid l, and air a. The Coleman–Noll procedure is used to obtain the restrictions on the form of the constitutive equations. The form of Darcy’s law for the fluid phase, which takes into account both Fickian and non-Fickian transport, is slightly different from the forms obtained by other researchers though all the terms have been included. When the fluid phases interact with the swelling solid porous matrix, deformation occurs. Viscoelastic large deformation of the solid matrix is investigated. A simple form of differential-integral equation is obtained for the fluid transport under isothermal conditions, which can be coupled with the deformation of the solid matrix to solve for transport in an unsaturated system. The modeling theory thus developed, which involves two-way coupling of the viscoelastic solid deformation and fluid transport, can be applied to study the processing of biopolymers, for example, soaking of foodstuffs and stress-crack predictions. Moreover, extension and modification of this modeling theory can be applied to study a vast variety of problems, such as drying of gels, consolidation of clays, drug delivery, and absorption of liquids in diapers.