A Parametric Model for Predicting Relative Permeability-Saturation-Capillary Pressure Relationships of Oil–Water Systems in Porous Media with Mixed Wettability

A parametric two-phase, oil–water relative permeability/capillary pressure model for petroleum engineering and environmental applications is developed for porous media in which the smaller pores are strongly water-wet and the larger pores tend to be intermediate- or oil-wet. A saturation index, which can vary from 0 to 1, is used to distinguish those pores that are strongly water-wet from those that have intermediate- or oil-wet characteristics. The capillary pressure submodel is capable of describing main-drainage and hysteretic saturation-path saturations for positive and negative oil–water capillary pressures. At high oil–water capillary pressures, an asymptote is approached as the water saturation approaches the residual water saturation. At low oil–water capillary pressures (i.e. negative), another asymptote is approached as the oil saturation approaches the residual oil saturation. Hysteresis in capillary pressure relations, including water entrapment, is modeled. Relative permeabilities are predicted using parameters that describe main-drainage capillary pressure relations and accounting for how water and oil are distributed throughout the pore spaces of a porous medium with mixed wettability. The capillary pressure submodel is tested against published experimental data, and an example of how to use the relative permeability/capillary pressure model for a hypothetical saturation-path scenario involving several imbibition and drainage paths is given. Features of the model are also explained. Results suggest that the proposed model is capable of predicting relative permeability/capillary pressure characteristics of porous media mixed wettability.     

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.     

Pore-by-Pore Modelling,Validation and Prediction of Waterflooding in Oil-Wet Rocks Using Dynamic Synchrotron Data

We predict waterflood displacement on a pore-by-pore basis using pore network modelling. The pore structure is captured by a high-resolution image. We then use an energy balance applied to images of the displacement to assign an average contact angle, and then modify the local pore-scale contact angles in the model about this mean to match the observed displacement sequence. Two waterflooding experiments on oil-wet rocks are analysed where the displacement sequence was imaged using time-resolved synchrotron imaging. In both cases the capillary pressure in the model matches the experimentally obtained values derived from the measured interfacial curvature. We then predict relative permeability for the full saturation range. Using the optimised contact angles distributed randomly in space has little effect on the predicted capillary pressures and relative permeabilities, indicating that spatial correlation in wettability is not significant in these oil-wet samples. The calibrated model can be used to predict properties outside the range of conditions considered in the experiment.

     

Pore-scale Network Modeling of Three-Phase Flow Based on Thermodynamically Consistent Threshold Capillary Pressures. II. Results

We use the model described in Zolfaghari and Piri (Transp Porous Media, 2016) to predict two- and three-phase relative permeabilities and residual saturations for different saturation histories. The results are rigorously validated against their experimentally measured counterparts available in the literature. We show the relevance of thermodynamically consistent threshold capillary pressures and presence of oil cusps for significantly improving the predictive capabilities of the model at low oil saturations. We study systems with wetting and spreading oil layers and cusps. Three independent experimental data sets representing different rock samples and fluid systems are investigated in this work. Different disordered networks are used to represent the pore spaces in which different sets of experiments were performed, i.e., Berea, Bentheimer, and reservoir sandstones. All three-phase equilibrium interfacial tensions used for the simulation of three-phase experimental data are measured and used in the model’s validation. We use a fixed set of parameters, i.e., the input network (to represent the pore space) and contact angles (to represent the wettability state), for all experiments belonging to a data set. Incorporation of the MSP method for capillary pressure calculations and cusp analysis significantly improves the agreement between the model’s predictions of relative permeabilities and residual oil saturations with experimental data.     

A Steady-State Upscaling Approach for Immiscible Two-Phase Flow

The paper presents a model for computing rate-dependent effective capillary pressure and relative permeabilities for two-phase flow, in 2 and 3 space-dimensions. The model is based on solving the equations for immiscible two-phase flow at steady-state, accounting for viscous and capillary forces, at a given external pressure drop. The computational performance of the steady-state model and its accuracy is evaluated through comparison with a commercial simulator ECLIPSE. The properties of the rate-dependent effective relative permeabilities are studied by way of computations using the developed steady-state model. Examples presented show the dependence of the effective relative permeabilities and capillary pressures, which incorporate the effects of fine scale wettability heterogeneity, on the external pressure drop, and thereby on the dimensionless macro-scale capillary number. The effective relative permeabilities converge towards the viscous limit functions as the capillary number tends to infinity. Special cases, when the effective relative permeabilities are rate-invariant, are also studied. The applicability of the steady-state upscaling algorithm in dynamic displacement situations is validated by comparing fine-gridded simulations in heterogeneous reservoirs against their homogenized counterparts. It is concluded that the steady-state upscaling method is able to accurately predict the dynamic behavior of a heterogeneous reservoir, including small scale heterogeneities in both the absolute permeability and the wettability.     

Three-Phase Capillary Pressure and Relative Permeability Relationships in Mixed-Wet Systems

A simple process-based model of three-phase displacement cycles for both spreading and non-spreading oils in a mixed-wet capillary bundle model is presented. All possible pore filling sequences are determined analytically and it is found that the number of pore occupancies that are permitted on physical grounds is actually quite restricted. For typical non-spreading gas/oil/water systems, only two important cases need to be considered to see all types of allowed qualitative behaviour for non-spreading oils. These two cases correspond to whether water or gas is the intermediate-wetting phase in oil-wet pores as determined by the corresponding contact angles, that is, cos ^o _gw > 0 or cos ^o _gw < 0, respectively. Analysis of the derived pore occupancies leads to the establishment of a number of relationships showing the phase dependencies of three-phase capillary pressures and relative permeabilities in mixed-wet systems. It is shown that different relationships hold in different regions of the ternary diagram and the morphology of these regions is discussed in terms of various rock/fluid properties. Up to three distinct phase-dependency regions may appear for a non-spreading oil and this reduces to two for a spreading oil. In each region, we find that only one phase may be specified as being the intermediate-wetting phase and it is only the relative permeability of this phase and the capillary pressure between the two remaining phases that depend upon more than one saturation. Given the simplicity of the model, a remarkable variety of behaviour is predicted. Moreover, the emergent three-phase saturation-dependency regions developed in this paper should prove useful in: (a) guiding improved empirical approaches of how two-phase data should be combined to obtain the corresponding three-phase capillary pressures and relative permeabilities; and (b) determining particular displacement sequences that require additional investigation using a more complete process-based 3D pore-scale network model.     

Two-Phase Pore-Network Modelling: Existence of Oil Layers During Water Invasion

Pore-network modelling is commonly used to predict capillary pressure and relative permeability functions for multi-phase flow simulations. These functions strongly depend on the presence of fluid films and layers in pore corners. Recently, van Dijke and Sorbie (J. Coll. Int. Sci. 293:455–463, 2006) obtained the new thermodynamically derived criterion for oil layers existence in the pore corners with non-uniform wettability caused by ageing. This criterion is consistent with the thermodynamically derived capillary entry pressures for other water invasion displacements and it is more restrictive than the previously used geometrical layer collapse criterion. The thermodynamic criterion has been included in a newly developed two-phase flow pore network model, as well as two versions of the geometrical criterion. The network model takes as input networks extracted from pore space reconstruction methods or CT images. Furthermore, a new n-cornered star shape characterization technique has been implemented, based on shape factor and dimensionless hydraulic radius as input parameters. For two unstructured networks, derived from a Berea sandstone sample, oil residuals have been estimated for different wettability scenarios, by varying the contact angles in oil-filled pores after ageing from weakly to strongly oil-wet. Simulation of primary drainage, ageing and water invasion show that the thermodynamical oil layer existence criterion gives more realistic oil residual saturations compared to the geometrical criteria. Additionally, a sensitivity analysis has been carried out of oil residuals with respect to end-point capillary pressures. For strongly oil-wet cases residuals increase strongly with increasing end-point capillary pressures, contrary to intermediate oil-wet cases.     

Relative Permeability Analysis of Tube Bundle Models,Including Capillary Pressure

The analytical equations for calculating two-phase flow, including local capillary pressures, are developed for the bundle of parallel capillary tubes model. The flow equations that are derived were used to calculate dynamic immiscible displacements of oil by water under the constraint of a constant overall pressure drop across the tube bundle. Expressions for averaged fluid pressure gradients and total flow rates are developed, and relative permeabilities are calculated directly from the two-phase form of Darcy's law. The effects of pressure drop and viscosity ratio on the relative permeabilities are discussed. Capillary pressure as a function of water saturation was delineated for several cases and compared to a steady-state mercury-injection drainage type of capillary pressure profile. The bundle of serial tubes model (a model containing tubes whose diameters change randomly at periodic intervals along the direction of flow), including local Young-Laplace capillary pressures, was analyzed with respect to obtaining relative permeabilities and macroscopic capillary pressures. Relative permeabilities for the bundle of parallel tubes model were seen to be significantly affected by altering the overall pressure drop and the viscosity ratio; relative permeabilities for the bundle of serial tubes were seen to be relatively insensitive to viscosity ratio and pressure, and were consistently X-like in profile. This work also considers the standard Leverett (1941) type of capillary pressure versus saturation profile, where drainage of a wetting phase is completed in a step-wise steady fashion; it was delineated for both tube bundle models. Although the expected increase in capillary pressure at low wetting-phase saturation was produced, comparison of the primary-drainage capillary pressure curves with the pseudo-capillary pressure profiles, that are computed directly using the averaged pressures during the displacements, revealed inconsistencies between the two definitions of capillary pressure.     

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.     

A Semi-Analytical Model for Two Phase Immiscible Flow in Porous Media Honouring Capillary Pressure

This article describes a semi-analytical model for two-phase immiscible flow in porous media. The model incorporates the effect of capillary pressure gradient on fluid displacement. It also includes a correction to the capillarity-free Buckley–Leverett saturation profile for the stabilized-zone around the displacement front and the end-effects near the core outlet. The model is valid for both drainage and imbibition oil–water displacements in porous media with different wettability conditions. A stepwise procedure is presented to derive relative permeabilities from coreflood displacements using the proposed semi-analytical model. The procedure can be utilized for both before and after breakthrough data and hence is capable to generate a continuous relative permeability curve unlike other analytical/semi-analytical approaches. The model predictions are compared with numerical simulations and laboratory experiments. The comparison shows that the model predictions for drainage process agree well with the numerical simulations for different capillary numbers, whereas there is mismatch between the relative permeability derived using the Johnson–Bossler–Naumann (JBN) method and the simulations. The coreflood experiments carried out on a Berea sandstone core suggest that the proposed model works better than the JBN method for a drainage process in strongly wet rocks. Both methods give similar results for imbibition processes.     

Percolation model of two-phase equilibrium flow in a porous medium with microheterogeneous wettability

The percolation model of two-phase flow described in [1, 2] is used as a basis for examining the problem of the behavior of the characteristics of two-phase equilibrium flow in a porous medium when the capillaries have a radius distribution and differ with respect to the wettability properties of their surfaces. Analytic expressions describing the dependence of the relative phase permeability coefficients on the saturation of the medium by the displacing phase and the microinhomogeneous wettability parameters are obtained. A qualitative comparison shows the theoretical results to be consistent with the data of a direct numerical computer calculation of a grid model [3]. The effect of the microinhomogeneity parameters and the form of the capillary radius distribution function on the phase permeabilities is analyzed within the framework of the approach developed.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 86–93, September–October, 1989.     

Pore Scale Modeling of Rate Effects in Imbibition

We use pore scale network modeling to study the effects of flow rate and contact angle on imbibition relative permeabilities. The model accounts for flow in wetting layers that occupy roughness or crevices in the pore space. Viscous forces are accounted for by solving for the wetting phase pressure and assuming a fixed conductance in wetting layers. Three-dimensional simulations model granular media, whereas two-dimensional runs represent fracture flow.We identify five generic types of displacement pattern as we vary capillary number, contact angle, and initial wetting phase saturation: flat frontal advance, dendritic frontal advance, bond percolation, compact cluster growth, and ramified cluster growth. Using phase diagrams we quantify the range of physical properties under which each regime is observed. The work explains apparently inconsistent experimental measurements of relative permeability in granular media and fractures.     

A Numerical Model of Tracer Transport in a Non-isothermal Two-Phase Flow System for {\text{ CO}}_2 Geological Storage Characterization

For the purpose of characterizing geologically stored $\text{ CO}_{2}Air sparging is an in situ soil/groundwater remediation technology, which involves the injection of pressurized air through air sparging well below the zone of contamination. To investigate the rate-dependent flow properties during multistep air sparging, a rule-based dynamic two-phase flow model was developed and applied to a 3D pore network which is employed to characterize the void structure of porous media. The simulated dynamic two-phase flow at the pore scale or microscale was translated into functional relationships at the continuum-scale of capillary pressure?Csaturation (P ^c?CS) and relative permeability??saturation (K _r?CS) relationships. A significant contribution from the air injection pressure step and duration time of each air injection pressure on both of the above relationships was observed during the multistep air sparging tests. It is observed from the simulation that at a given matric potential, larger amount of water is retained during transient flow than that during steady flow. Shorter the duration of each air injection pressure step, there is higher fraction of retained water. The relative air/water permeability values are also greatly affected by the pressure step. With large air injection pressure step, the air/water relative permeability is much higher than that with a smaller air injection pressure step at the same water saturation level. However, the impact of pressure step on relative permeability is not consistent for flows with different capillary numbers (N _ca). When compared with relative air permeability, relative water permeability has a higher scatter. It was further observed that the dynamic effects on the relative permeability curve are more apparent for networks with larger pore sizes than that with smaller pore sizes. In addition, the effect of pore size on relative water permeability is higher than that on relative air permeability.     

Generalized relative permeability coefficients during steady-state two-phase flow in porous media,and correlation with the flow mechanisms

A parametric experimental investigation of the coupling effects during steady-state two-phase flow in porous media was carried out using a large model pore network of the chamber-and-throat type, etched in glass. The wetting phase saturation,S ₁, the capillary number,Ca, and the viscosity ratio,k, were changed systematically, whereas the wettability (contact angleθ _e ), the coalescence factorCo, and the geometrical and topological parameters were kept constant. The fluid flow rate and the pressure drop were measured independently for each fluid. During each experiment, the pore-scale flow mechanisms were observed and videorecorded, and the mean water saturation was determined with image analysis. Conventional relative permeability, as well as generalized relative permeability coefficients (with the viscous coupling terms taken explicitly into account) were determined with a new method that is based on a B-spline functional representation combined with standard constrained optimization techniques. A simple relationship between the conventional relative permeabilities and the generalized relative permeability coefficients is established based on several experimental sets. The viscous coupling (off-diagonal) coefficients are found to be comparable in magnitude to the direct (diagonal) coefficients over board ranges of the flow parameter values. The off-diagonal coefficients (k _rij /Μ _j ) are found to be unequal, and this is explained by the fact that, in the class of flows under consideration, microscopic reversibility does not hold and thus the Onsager-Casimir reciprocal relation does not apply. Thecoupling indices are introduced here; they are defined so that the magnitude of each coupling index is the measure of the contribution of the coupling effects to the flow rate of the corresponding fluid. A correlation of the coupling indices with the underlying flow mechanisms and the pertinent flow parameters is established.     

Fluid flow properties for different classes of intermediate wettability as studied by network modelling

Most clastic reservoirs display an intermediate type of wettability. Intermediate wettability covers several local wetting configurations like fractional wet and mixed-wet where the oil-wet sites could either be in the large or smaller pores. Clastic reservoirs show a large variation in fluid flow properties. A classical invasion–percolation network simulator is used to investigate properties of different intermediate wet situations. Variation in wetting properties like contact angles, process dependent contact angles, contact angle distribution, and fraction of oil wet sites are investigated. The fluid flow properties analysed in particular are residual oil saturation and normalized endpoint relative permeability. Results from network modelling have been compared to reservoir core analysis data. The network models applied are at the capillary limit, while the core flood results are clearly viscous influenced. Even though network modelling does not cover all the physics involved in fluid displacement processes, results show that data from simulations are sufficient to present trends in fluid flow properties which are comparable to experimental data.     

Rapoport-leas two-phase flow model for anisotropic porous media

The Rapoport-Leas mathematical model of two-phase flow is generalized to include the case of anisotropic porous media. The formula for the capillary pressure, which specifies the relationship between the phase pressures, contains a scalar function of a vector argument. In order to determine the scalar function, the capillary pressure tensor and the tensor inverse to the tensor of characteristic linear dimensions are introduced. The capillary pressure is determined by the contraction of the second-rank tensors with a unit vector collinear to the phase pressure gradients, also assumed to be collinear. It is shown that the saturation function introduced for isotropic porous media (Leverett function) can be generalized to include anisotropic media and is now determined by a fourth-rank tensor. Generalized expressions for the Leverett and relative phase permeability functions are given for orthotropic and transversely isotropic media with account for the hysteresis of the phase permeabilities and capillary pressure.     

Relative permeabilities from two- and three-dimensional pore-scale network modelling

We present a computer study of two-phase flow in a porous medium. The porous medium is represented by an isotropic network of up to 80 000 randomly placed nodes connected by thin tubes. We then simulate two-fluid displacements in this network and are able to demonstrate the effects of viscous and capillary forces. We use the local average flow rates and pressures to calculate effective saturation dependent relative pemeabilities, fractional flows and capillary pressures. Using a radial Buckley-Leverett theory, the mean saturation profile can be inferred from the solution of the fractional flow equation, which is consistent with the computed saturation. We show that the relative permeability may be a function of both viscosity ratio and capillary number.     

Elliptic Regions and Stable Solutions for Three-Phase flow in Porous Media

In the limit of zero capillary pressure, solutions to the equations governing three-phase flow, obtained using common empirical relative permeability models, exhibit complex wavespeeds for certain saturation values (elliptic regions) that result in unstable and non-unique solutions. We analyze a simple but physically realizable pore-scale model: a bundle of cylindrical capillary tubes, to investigate whether the presence of these elliptic regions is an artifact of using unphysical relative permeabilities. Without gravity, the model does not yield elliptic regions unless the most non-wetting phase is the most viscous and the most wetting phase is the least viscous. With gravity, the model yields elliptic regions for any combination of viscosities, and these regions occupy a significant fraction of the saturation space. We then present converged, stable numerical solutions for one-dimensional flow, which include capillary pressure. These demonstrate that, even when capillary forces are small relative to viscous forces, they have a significant effect on solutions which cross or enter the elliptic region. We conclude that elliptic regions can occur for a physically realizable model of a porous medium, and that capillary pressure should be included explicitly in three-phase numerical simulators to obtain stable, physically meaningful solutions which reproduce the correct sequence of saturation changes.     

Experimental study of two-phase flow in three sandstones. II. Capillary pressure curve measurement and relative permeability pore space capillary models

By means of the porous plate method and mercury porosimetry intrusion tests, capillary pressure curves of three different sandstones were measured. The testing results have been exploited jointly with three relative permeability models of the pore space capillary type (Burdine’s model type), these models are widely used and in rather distinct fields. To do so, capillary pressure has been correlated to saturation degree using six of the most popular relations encountered in the literature. Model predictions were systematically compared to the experimentally measured relative permeabilities presented in the first part of this work. Comparison indicated that the studied models underestimate the water relative permeability and over-estimate that of the non-wetting phase. Moreover, this modeling proves to be unable to locate the significant points that are the limits of fields of saturation where the variation of the relative permeabilities becomes consequent. We also showed that, if pore structure is modeled as a “bundle of capillary tubes”, model predications are independent of the capillary pressure curve measuring method.     

Relative Permeability Calculations from Two-Phase Flow Simulations Directly on Digital Images of Porous Rocks

We present results from a systematic study of relative permeability functions derived from two-phase lattice Boltzmann (LB) simulations on X-ray microtomography pore space images of Bentheimer and Berea sandstone. The simulations mimic both unsteady- and steady-state experiments for measuring relative permeability. For steady-state flow, we reproduce drainage and imbibition relative permeability curves that are in good agreement with available experimental steady-state data. Relative permeabilities from unsteady-state displacements are derived by explicit calculations using the Johnson, Bossler and Naumann method with input from simulated production and pressure profiles. We find that the nonwetting phase relative permeability for drainage is over-predicted compared to the steady-state data. This is due to transient dynamic effects causing viscous instabilities. Thus, the calculated unsteady-state relative permeabilities for the drainage is fundamentally different from the steady-state situation where transient effects have vanished. These effects have a larger impact on the invading nonwetting fluid than the defending wetting fluid. Unsteady-state imbibition relative permeabilities are comparable to the steady-state ones. However, the appearance of a piston-like front disguises most of the displacement and data can only be determined for a restricted range of saturations. Relative permeabilities derived from unsteady-state displacements exhibit clear rate effects, and residual saturations depend strongly on the capillary number. We conclude that the LB method can provide a versatile tool to compute multiphase flow properties from pore space images and to explore the effects of imposed flow and fluid conditions on these properties. Also, dynamic effects are properly captured by the method, giving the opportunity to examine differences between steady and unsteady-state setups.