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...     

An Analysis of Three-Phase Pore Occupancies and Relative Permeabilities in Porous Media with Variable Wettability

In three-phase flow, the macroscopic constitutive relations of capillary pressure and relative permeability as functions of saturation depend in a complex manner on the underlying pore occupancies. These three-phase pore occupancies depend in turn on the interfacial tensions, the pore sizes and the degree of wettability of the pores, as characterised by the cosines of the oil–water contact angles. In this work, a quasi-probabilistic approach is developed to determine three-phase pore occupancies in media where the degree of wettability varies from pore to pore. Given a set of fluid and rock properties, a simple but novel graphical representation is given of the sizes and oil–water contact angles underlying three-phase occupancies for every allowed combination of capillary pressures. The actual phase occupancies are then computed using the contact angle probability density function. Since a completely accessible porous medium is studied, saturations, capillary pressures, and relative permeabilities are uniquely related to the pore occupancies. In empirical models of three-phase relative permeability it is of central importance whether a phase relative permeability depends only on its own saturation and how this relates to the corresponding two-phase relative permeability (if at all). The new graphical representation of pore sizes and wettabilities clearly distinguishes all three-phase pore occupancies with respect to these saturation-dependencies. Different types of saturation-dependencies may occur, which are shown to appear in ternary saturation diagrams of iso-relative permeability curves as well, thus guiding empirical approaches. However, for many saturation combinations three-phase and two-phase relative permeabilities can not be linked. In view of the latter, the present model has been used to demonstrate an approach for three-phase flow modelling on the basis of the underlying pore-scale processes, in which three-phase relative permeabilities are computed only along the actual flow paths. This process-based approach is used to predict an efficient strategy for oil recovery by simultaneous water-alternating-gas (SWAG) injection.     

Analysis of Microscopic Displacement Mechanisms of a MIOR Process in Porous Media with Different Wettability

Transport in Porous Media - The wettability of the reservoir rock has an important effect on the displacement of fluids on a microscopic scale in all EOR processes, especially in the microbial...     

Two-Phase Flow in Heterogeneous Porous Media with Non-Wetting Phase Trapping

This article presents a mathematical model describing flow of two fluid phases in a heterogeneous porous medium. The medium contains disconnected inclusions embedded in the background material. The background material is characterized by higher value of the non-wetting-phase entry pressure than the inclusions, which causes non-standard behavior of the medium at the macroscopic scale. During the displacement of the non-wetting fluid by the wetting one, some portions of the non-wetting fluid become trapped in the inclusions. On the other hand, if the medium is initially saturated with the wetting phase, it starts to drain only after the capillary pressure exceeds the entry pressure of the background material. These effects cannot be represented by standard upscaling approaches based on the assumption of local equilibrium of the capillary pressure. We propose a relevant modification of the upscaled model obtained by asymptotic homogenization. The modification concerns the form of flow equations and the calculation of the effective hydraulic functions. This approach is illustrated with two numerical examples concerning oil–water and CO₂–brine flow, respectively.     

Accurate Modelling of Pore-Scale Films and Layers for Three-Phase Flow Processes in Clastic and Carbonate Rocks with Arbitrary Wettability

Three-phase flow is a key process occurring in subsurface reservoirs, for example, during $\text{ CO }_2$ sequestration and enhanced oil recovery techniques such as water alternating gas (WAG) injection. Predicting three-phase flow processes, for example, the increase in oil recovery during WAG, requires a sound understanding of the fundamental flow physics in water- to oil-wet rocks to derive physically robust flow functions, i.e. relative permeability and capillary pressure. In this study, we use pore-network modelling, a reliable and physically based simulation tool, to predict the flow functions. We have developed a new pore-scale network model for rocks with variable wettability, from water- to oil-wet. It comprises a constrained set of parameters that mimic the wetting state of a reservoir. Unlike other models, it combines three main features: (1) A novel thermodynamic criterion for formation and collapse of oil layers. The new model hence captures wetting film and layer flow of oil adequately, which affects the oil relative permeability at low oil saturation and leads to accurate prediction of residual oil. (2) Multiple displacement chains, where injection of one phase at the inlet triggers a chain of interface displacements throughout the network. This allows for the accurate modelling of the mobilisation of many disconnected phase clusters that arise, in particular, during higher order WAG floods. (3) The model takes realistic 3D pore-networks extracted from pore-space reconstruction methods and CT images as input, preserving both topology and pore shape of the sample. For water-wet systems, we have validated our model with available experimental data from core floods. For oil-wet systems, we validated our network model by comparing 2D network simulations with published data from WAG floods in oil-wet micromodels. This demonstrates the importance of film and layer flow for the continuity of the various phases during subsequent WAG cycles and for the residual oil saturations. A sensitivity analysis has been carried out with the full 3D model to predict three-phase relative permeabilities and residual oil saturations for WAG cycles under various wetting conditions with different flood end-points.     

Empirical Measures of Wettability in Porous Media and the Relationship between Them Derived From Pore-Scale Modelling

The wettability of a crude oil/brine/rock system is of central importance in determining the oil recovery efficiency of water displacement processes in oil reservoirs. Wettability of a rock sample has traditionally been measured using one of two experimental techniques, viz. the United States Bureau of Mines and Amott tests. The former gives the USBM index, I _USBM, and the latter yields the Amott–Harvey index, I _AH. As there is no well-established theoretical basis for either test, any relationship between the two indices remains unclear.Analytical relationships between I _AH and I _USBM for mixed-wet and fractionally-wet media have been based on a number of simplifying assumptions relating to the underlying pore-scale displacement mechanisms. This simple approach provides some guidelines regarding the influence of the distribution of oil-wet surfaces within the porous medium on I _AH and I _USBM. More detailed insight into the relationship between I _AH and I _USBM is provided by modelling the pore-scale displacement processes in a network of interconnected pores. The effects of pore size distribution, interconnectivity, displacement mechanisms, distribution of volume and of oil-wet pores within the pore space have all been investigated by means of the network model.The results of these analytical calculations and network simulations show that I _AH and I _USBM need not be identical. Moreover, the calculated indices and the relationship between them suggest explanations for some of the trends that appear in experimental data when both I _USBM and I _AH have been reported in the literature for tests with comparable fluids and solids. Such calculations should help with the design of more informative wettability tests in the future.     

Pore-Scale Simulations of Single- and Two-Phase Flow in Porous Media: Approaches and Applications

Transport in Porous Media - We present a review of pore-scale simulations of immiscible fluid transport with focus on two of the most popular approaches: lattice Boltzmann modeling for direct...     

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.     

Pore-Scale Analysis of NAPL Blob Dissolution and Mobilization in Porous Media

A pore-scale analysis of nonaqueous phase liquid (NAPL) blob dissolution and mobilization in porous media was presented. Dissolution kinetics of residual NAPLs in an otherwise water-saturated porous medium was investigated by conducting micromodel experiments. Changes in residual NAPL volume were measured from recorded video images to calculate the mass transfer coefficient, K and the lumped mass transfer rate coefficient, k. The morphological characteristics of the blobs such as specific and intrinsic area were found to be independent of water flow rate except at NAPL saturations below 2%. Dissolution process was also investigated by separating the mass transfer into zones of mobile and immobile water. The fractions of total residual NAPL perimeters in contact with mobile water and immobile water were measured and their relationship to the mass transfer coefficient was discussed. In general, residual NAPLs are removed by dissolution and mobilization. Although these two mechanisms were studied individually by others, their simultaneous occurrence was not considered. Therefore, in this study, mobilization of dissolving NAPL blobs was investigated by an analysis of the forces acting on a trapped NAPL blob. A dimensional analysis was performed to quantify the residual blob mobilization in terms of dimensionless Capillary number (Ca _I). If Ca _I is equal to or greater than the trapping number defined as , then blob mobilization is expected.     

High-Resolution Temporo-Ensemble PIV to Resolve Pore-Scale Flow in 3D-Printed Fractured Porous Media

Transport in Porous Media - Fractures are conduits that can enable fast advective transfer of (fluid, solute, reactant, particle, etc.) mass and energy. Such fast transfer can significantly affect...     

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.     

Analytical Solution to the Riemann Problem of Three-Phase Flow in Porous Media

In this paper we study one-dimensional three-phase flow through porous media of immiscible, incompressible fluids. The model uses the common multiphase flow extension of Darcys equation, and does not include gravity and capillarity effects. Under these conditions, the mathematical problem reduces to a 2 × 2 system of conservation laws whose essential features are: (1) the system is strictly hyperbolic; (2) both characteristic fields are nongenuinely nonlinear, with single, connected inflection loci. These properties, which are natural extensions of the two-phase flow model, ensure that the solution is physically sensible. We present the complete analytical solution to the Riemann problem (constant initial and injected states) in detail, and describe the characteristic waves that may arise, concluding that only nine combinations of rarefactions, shocks and rarefaction-shocks are possible. We demonstrate that assuming the saturation paths of the solution are straight lines may result in inaccurate predictions for some realistic systems. Efficient algorithms for computing the exact solution are also given, making the analytical developments presented here readily applicable to interpretation of lab displacement experiments, and implementation of streamline simulators.     

Three-Phase Flow Modelling Using Pore-Scale Capillary Pressures and Relative Permeabilities for Mixed-Wet Media at the Continuum-Scale

When regions of three-phase flow arise in an oil reservoir, each of the flow parameters, i.e. capillary pressures and relative permeabilities, are generally functions of two phase saturations and depend on the wettability state. The idea of this work is to generate consistent pore-scale based three-phase capillary pressures and relative permeabilities. These are then used as input to a 1-D continuum core- or reservoir-scale simulator. The pore-scale model comprises a bundle of cylindrical capillary tubes, which has a distribution of radii and a prescribed wettability state. Contrary to a full pore-network model, the bundle model allows us to obtain the flow functions for the saturations produced at the continuum-scale iteratively. Hence, the complex dependencies of relative permeability and capillary pressure on saturation are directly taken care of. Simulations of gas injection are performed for different initial water and oil saturations, with and without capillary pressures, to demonstrate how the wettability state, incorporated in the pore-scale based flow functions, affects the continuum-scale displacement patterns and saturation profiles. In general, wettability has a major impact on the displacements, even when capillary pressure is suppressed. Moreover, displacement paths produced at the pore-scale and at the continuum-scale models are similar, but they never completely coincide.     

Relative Permeabilities for Strictly Hyperbolic Models of Three-Phase Flow in Porous Media

Traditional mathematical models of multiphase flow in porous media use a straightforward extension of Darcys equation. The key element of these models is the appropriate formulation of the relative permeability functions. It is well known that for one-dimensional flow of three immiscible incompressible fluids, when capillarity is neglected, most relative permeability models used today give rise to regions in the saturation space with elliptic behavior (the so-called elliptic regions). We believe that this behavior is not physical, but rather the result of an incomplete mathematical model. In this paper we identify necessary conditions that must be satisfied by the relative permeability functions, so that the system of equations describing three-phase flow is strictly hyperbolic everywhere in the saturation triangle. These conditions seem to be in good agreement with pore-scale physics and experimental data.     

Pore-Scale Simulation of Fluid Flow Through Deformable Porous Media Using Immersed Boundary Coupled Lattice Boltzmann Method

Transport in Porous Media - The interaction between the fluid flow and the deformable porous media is crucial in the applications of adsorption/absorption. The immersed boundary coupled lattice...     

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.     

Capillary Heterogeneity Trapping and Crossflow in Layered Porous Media

We examine the effect of capillary and viscous forces on the displacement of one fluid by a second, immiscible fluid across and along parallel layers of contrasting porosity, and relative permeability, as well as previously explored contrasts in absolute permeability and capillary pressure. We consider displacements with wetting, intermediate-wetting and non-wetting injected phases. Flow is characterized using six independent dimensionless numbers and a dimensionless storage efficiency, which is numerically equivalent to the recovery efficiency. Results are directly applicable to geologic carbon storage and hydrocarbon production. We predict how the capillary–viscous force balance influences storage efficiency as a function of a small number of key dimensionless parameters, and provide a framework to support mechanistic interpretations of complex field or experimental data, and numerical model predictions, through the use of simple dimensionless models. When flow is directed across layers, we find that capillary heterogeneity traps the non-wetting phase, regardless of whether it is the injected or displaced phase. However, minimal trapping occurs when the injected phase is intermediate-wetting or when high-permeability layers contain a smaller moveable volume of fluid than low-permeability layers. A dimensionless capillary-to-viscous number defined using the layer thickness rather than the more commonly used system length is most relevant to predict capillary heterogeneity trapping. When flow is directed along layers, we show that, regardless of wettability, increasing capillary crossflow reduces the distance between the leading edges of the injected phase in each layer and increases storage efficiency. This may be counter-intuitive when the injected phase is non-wetting. Crossflow has a significant impact on storage efficiency only when high-permeability layers contain a smaller moveable volume of fluid than low-permeability layers. In that case, capillary heterogeneity traps the wetting phase, regardless of whether it is the injected or displaced phase.     

A Pore-Scale Investigation of Low-Salinity Waterflooding in Porous Media: Uniformly Wetted Systems

The potential of low-salinity (LS) water injection as an oil recovery technique has been the source of much recent debate within the petroleum industry. Evidence from both laboratory and field-level studies has indicated significant benefits compared to conventional high-salinity (HS) waterflooding, but many conflicting results have also been reported and, to date, the underlying mechanisms remain poorly understood. In this paper, we aim to address this uncertainty by developing a novel, steady-state pore network model in which LS brine displaces oil from a HS-bearing network. The model allows systematic investigation of the crude oil/brine/rock parameter space, with the goal of identifying features that may be critical to the production of incremental oil following LS brine injection. By coupling the displacement model to a salinity-tracking tracer algorithm, and assuming that a reduction of water salinity within the pore network leads to localised wettability alteration, substantial perturbations to standard pore filling sequences are predicted. The results clearly point to two principal effects of dynamic contact angle modification at the pore scale: a “pore sequence” effect, characterised by an alteration to the distribution of displaced pore sizes, and a “sweep efficiency” effect, demonstrated by a change in the overall fraction of pores invaded. Our study indicates that any LS effect will depend on the relative (scenario-dependent) influence of each mechanism, where factors such as the initial wettability state of the system and the pore size distribution of the underlying network are found to play crucial roles. In addition, we highlight the important role played by end-point capillary pressure in determining LS efficacy.