Flow behaviour in the presence of wettability heterogeneities

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

Experimental displacement patterns in a 2 × 2 quadrant block with permeability and wettability heterogeneities—problems for numerical modelling

The effect of heterogeneities on miscible and immiscible flood displacements in 2D bead packs in quadrant form, 2 × 2 block heterogeneity, with either a permeability or a wettability contrast is the subject of this paper. The physical processes occurring during miscible and immiscible flow and displacement within permeability and wettability quadrant bead pack models have been studied experimentally. This geometry occurs in a number of situations relevant to hydrocarbon production: particularly faults where adjacent rocks have large permeability contrasts with rapid changes, in the laboratory with core butting, in reservoir simulation where grid blocks have different permeability and in reservoirs having near-wellbore damage problems. The model quadrants 1–4, had 1 and 4 and 2 and 3 with identical properties, either in permeability or wettability. Reported are complete unit mobility miscible displacements, then the effects of viscosity differences (mobility modifiers) and finally immiscible displacements on displacement patterns for initial linear injection. The experiments demonstrate that nodal flow occurs for both miscible and immiscible flow, but for immiscible flow there are boundary effects due to capillary pressure differences created by water saturation changes or wettability contrasts which can leave patches of isolated fluid within a quadrant. The displacement patterns for the different models and fluids change significantly with the viscosity and wettability changes, particularly for the immiscible displacements. This is due to the changing capillary pressure between the quadrant blocks as the water saturation change. These are difficult to address in numerical modelling but should be accounted for. Other effects include coupling of all physical processes governing the flow through the node and creations of microzones of trapped residual oil. Our displacement patterns can therefore be a valuable verification benchmark tool for numerical modelling and a calibration data source for those wishing to simulate the effects of capillary pressure under differing wettability conditions and for those investigating upscaling modelling procedures. However, the possible loss of physical reality when averaging must always be considered.     

Simulation of three-phase displacement experiments

Three-phase displacement experiments for a water-benzyl alcohol-decane system are simulated. Literature experimental three-phase relative permeabilities for the system are used to describe the relative permeabilities in the three-phase region for different three-phase relative permeability models. Saturation trajectories and elliptical regions are mapped in the three-phase region. Simulations are performed to model displacement experiments including breakthrough and the formation of multiple shocks. The model can be used to predict the results for other displacements. In an experiment where significant gravity segregation is present, the displacement is more accurately modeled by assuming a uniform initial condition than by using the actual vertical saturation and assuming no cross flow. It is shown how different residual saturation values can be measured in the laboratory depending on the initial saturation conditions in the core. The experimental residual saturations can be significantly different than the ‘theoretical’ or model values.     

An Experimental Method for One Dimensional Permeability Characterization of Heterogeneous Porous Media at the Core Scale

Many reservoir simulator inputs are derived from laboratory experiments. Special core analysis techniques generally assume that core samples are homogeneous. This assumption does not hold for porous media with significant heterogeneities. This paper presents a new method to characterize core scale permeability heterogeneity. The method is validated by both numerical and experimental results. The leading idea consists in injecting a high viscosity miscible fluid into a core sample saturated with a low viscosity fluid. In such conditions, the fluid displacement is expected to be piston-like. We investigate the evolution of the pressure drop as a function of time. A continuous permeability profile is estimated along flow direction from the pressure drop assuming that the core sample is a stack of infinitely thin cross sections perpendicular to flow direction. Thus, we determine a permeability value for each cross section. Numerical and laboratory experiments are carried out to validate the method. Flow simulations are performed for numerical models representing core samples to estimate the pressure drop. The selected models are sequences of plugs with constant permeabilities. In addition, laboratory displacements are conducted for both low permeability and high permeability core samples. To investigate whether there is dispersion inside the porous medium, CT scan measurements are performed during fluid displacement: the location of the front is observed at successive time intervals. The results validate the methodology developed in this paper as long as heterogeneity is one dimensional.     

Numerical Analysis of Imbibition Front Evolution in Fractured Sandstone under Capillary-Dominated Conditions

Fractures serve as primary conduits having a great impact on the migration of injected fluid into fractured permeable media. Appropriate transport properties such as relative permeability and capillary pressure are essential for successful simulation and prediction of multi-phase flow in such systems. However, the lack of a thorough understanding of the dynamics governing immiscible displacement in fractured media, limits our ability to properly represent their macroscopic transport properties. Previous experimental observations of imbibition front evolution in fractured rocks are examined in the present study using an automated history-matching approach to obtain representative relative permeability and capillary pressure curves. Predicted imbibition front evolution under different flow conditions resulted in an excellent agreement with experimental observations. Sensitivity analyses, in combination with direct experimental observation, allowed exploring the competing effects of relative permeability and capillary pressure on the development of saturation distribution and imbibing front evolution in fractured porous media. Results show that residual saturations are most sensitive to matrix relative permeability to oil, while the ratio of oil and water relative permeability, rock heterogeneity, boundary condition, and matrix–fracture capillary pressure contrast, affect displacement shape, speed, and geometry of the imbibing front.     

Modelling fluid flow in cross-bedded sections

Beadpack experiments and numerical simulations have been carried out to study flow displacements, effluent profiles and streamline patterns for layered systems with flow not parallel to the layers. The effects of layer thickness, permeability contrast, angle of layer to flow direction, mobility ratio and flood rate have been examined. Each of these parameters influence the displacement profiles, and disperse the flood front. Such real effects must be considered when subsuming reservoir heterogeneities in average reservoir parameters in simulation studies, or interpreting core tests.     

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.     

An Experimental Study of Three-Phase Trapping in Sand Packs

The trapped saturations of oil and gas are measured as functions of initial oil and gas saturation in water-wet sand packs. Analogue fluids—water, octane and air—are used at ambient conditions. Starting with a sand-pack column which has been saturated with brine, oil (octane) is injected with the column horizontal until irreducible water saturation is reached. The column is then positioned vertically and air is allowed to enter from the top of the column, while oil is allowed to drain under gravity for varying lengths of time. At this point, the column may be sliced and the fluids analyzed by gas chromatography to obtain the initial saturations. Alternatively, brine is injected through the bottom of the vertical column to trap oil and gas, before slicing the columns and measuring the trapped or residual saturations by gas chromatography and mass balance. The experiments show that in three-phase flow, the total trapped saturations of oil and gas are considerably higher than the trapped saturations reported in the literature for two-phase systems. It is found that the residual saturation of oil and gas combined could be as high as 23 %, as opposed to a maximum two-phase residual of only 14 %. For very high initial gas saturations, the residual gas saturation, up to 17 %, was also higher than for two-phase displacement. These observations are explained in terms of the competition between piston-like displacement and snap-off. It is also observed that less oil is always trapped in three-phase flow than in two-phase displacement, and the difference depends on the amount of gas present. For low and intermediate initial gas saturations, the trapped gas saturation rises linearly with initial saturation, followed by a constant residual, as seen in two-phase displacements. However, at very high initial gas saturations, the residual saturation rises again.     

Invasion Percolation on Correlated and Elongated Lattices: Implications for the Interpretation of Residual Saturations in Rock Cores

The invasion percolation model is used to investigate the effect of correlated heterogeneity on capillary dominated displacements in porous media. The breakthrough and residual saturations are shown to be strongly influenced by the correlations. Correlated heterogeneity leads to lower residual saturations than those observed in random systems and the scatter commonly observed in laboratory core measurements of the residual saturations can be attributed to the presence of such heterogeneity at the pore scale. Invasion percolation computations on elongated lattices, those with a geometry of L ^d–1 × nL where n denotes the aspect ratio, show that residual saturations for systems with correlated heterogeneity exhibit a strong dependence on aspect ratio. This effect is not considered by experimentalists who advocate the use of long (high aspect ratio) cores in order to overcome end-effects in experiments on shorter cores. A new scaling law is proposed for the residual saturations in elongated systems with correlated heterogeneity, and is confirmed by numerical simulations.     

Fluid Flow Across Capillary Heterogeneities,Experiments and Simulations

A comparative study of numerical modelling and laboratory experiments of two-phase immiscible displacements in a 33 cm × 10 × 3 cm thick cross-bedded reservoir model is reported. Dynamic two-dimensional fluid saturation development was obtained from experiments by use of a nuclear tracer imaging technique and compared to numerical predictions using a full-field black oil simulator.The laboratory cross-bedded reservoir model was a sandpack consisting of two strongly waterwet sands of different grain sizes, packed in sequential layers. The inlet and outlet sand consisted of low permeable, high capillary, sand while the central crosslayer with a dip angle of 30° was a high permeable, low capillary, sand. Results on moderate contrasts in permeability and capillary heterogeneities in the cross-bedded reservoir model at different mobility ratios and capillary number floods temporarily showed a bypass of oil, resulting in a prolonged two-phase production. The final remaining oil saturations, however, were as for isolated samples. Hence, permanently trapped oil was not observed.Simulations of waterfloods, using a commercial software package, displayed correct water breakthrough at low flow rate and unity viscosity ratio, but failed in predicting local saturation development in detail, probably due to numerical diffusion.The simulator was used to test several cases of heterogeneity contrasts, and influence from different relative permeability curves. Further, by altering the capillary pressure at the outlet, the end effects were proven important.     

Experimental investigation of capillary pressure effects on immiscible displacement in lensed and layered porous media

This paper reports the results of extensive experimental studies of the effects of well-defined heterogeneous porous media on immiscible flooding. The heterogeneities were layers and lenses, with some of the lenses being a wettability contrast. Drainage and imbibition displacements, with and without an initial residual fluid saturation, were carried out at a variety of flow rates on layered and lensed two-dimensional glass beads models of the size of a typical large core test (58×10×0.6 cm). These displacements were followed photographically and the effluent saturation profiles recorded. In most of the experiments the glass beads were water-wet, but in some the lens beads were coated with a water repellent chemical. In all experiments, the displacement fronts became highly irregular due to the different capillary pressures acting in the different areas of the models. In this paper, these displacements are fully reported and their implications for reservoir simulation and for interpretation of laboratory core tests, where the inner heterogeneities are not known, are discussed.     

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.     

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.     

Stochastic Simulations of NAPL Mass Transport in Variably Saturated Heterogeneous Porous Media

A multiphase flow and transport numerical model is developed to study the effects of porous media heterogeneities on residual NAPL mass partitioning and transport of dissolved and/or volatilized NAPL mass in variably saturated media. The results indicate the significance of porous media heterogeneity in influencing the mass transfer processes and NAPL transport in the subsurface. Among the parameters investigated in this study, the heterogeneity of the permeability field has the most significant influence on the NAPL mass partitioning and transport. In general, the heterogeneity of the porous media properties enhances the NAPL mass plume spreading in both the water phase and the gas phase while the influence on the water phase is much more significant. Overall, the porous media property heterogeneities tend to increase the accumulation of NAPL mass in the water phase. The nonequilibrium mass transfer processes result in the expected trend of decreasing the NAPL mass dissipation rate and causing long-term groundwater contamination.     

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.     

An Experimental Study on the Influence of Sub-Core Scale Heterogeneities on CO<Subscript>2</Subscript> Distribution in Reservoir Rocks

This article presents the results of CO₂/brine two-phase flow experiments in rocks at reservoir conditions. X-ray CT scanning is used to determine CO₂ saturation at a fine scale with a resolution of a few pore volumes and provide 3D porosity and saturation maps that can be use to correlate CO₂ saturations and rock properties. The study highlights the strong influence of sub-core scale heterogeneities on the spatial distribution of CO₂ at steady state and provides useful relative permeability data on a sample originated from an actual storage site (CO2CRC-Otway project, Victoria, South-West Australia). Two different samples tested, although different in nature, present strong heterogeneities, but differ in the detail of the connectivity of high porosity layers. In both samples, the results of the investigations show that sub-core scale heterogeneities control the sweep efficiency and may cause channeling through the porous medium. In one of the samples, CO₂ saturation appears uncorrelated to porosity close to the outlet end of the core. This observation is understood as a result of the position and the orientation of high porosity layers with respect to the inlet face of the core. Finally, in the operating conditions of the two experiments, the saturation maps demonstrate that gravity does not play a major role since no detectable buoyancy driven flow is observed.     

The Development of Appropriate Upscaling Procedures

Permeability upscaling should be carried out with careful attention to the nature of rock heterogeneities. While there are many large-scale features which must be taken into account, there are also important heterogeneities at the small-scale. Many sedimentary structures contain laminae at the mm–cm scale, and beds at the m-scale, which give rise to strong contrasts in permeability. We use a 2D model derived from a photo-panel of an aeolian outcrop, along with permeability measurements from a North Sea oil field, to demonstrate the effects of small-scale heterogeneity. This model is similar in size to a typical cell of a reservoir geological model. We take imaginary probe and core plug measurements from the model, average them, and compare these with the effective permeability for the model computed from a finite difference flow calculation. Although this procedure is standard practice, we show that it can lead to biased estimates of the permeabilities used in flow simulation. As an alternative we suggest using models of representative beds, and performing flow simulation to calculate effective permeabilities for both single-phase and two-phase flow.     

A Fickian diffusion model for the spreading of liquid plumes infiltrating in heterogeneous media

Infiltration of water and non-aqueous phase liquids (NAPLs) in the vadose zone gives rise to complex two- and three-phase immiscible displacement processes. Physical and numerical experiments have shown that ever-present small-scale heterogeneities will cause a lateral broadening of the descending liquid plumes. This behavior of liquid plumes infiltrating in the vadose zone may be similar to the familiar transversal dispersion of solute plumes in single-phase flow. Noting this analogy we introduce a mathematical model for ‘phase dispersion’ in multiphase flow as a Fickian diffusion process. It is shown that the driving force for phase dispersion is the gradient of relative permeability, and that addition of a phase-dispersive term to the governing equations for multiphase flow is equivalent to an effective capillary pressure which is proportional to the logarithm of the relative permeability of the infiltrating liquid phase. The relationship between heterogeneity-induced phase dispersion and capillary and numerical dispersion effects is established. High-resolution numerical simulation experiments in heterogeneous media show that plume spreading tends to be diffusive, supporting the proposed convection-dispersion model. Finite difference discretization of the phase-dispersive flux is discussed, and an illustrative application to NAPL infiltration from a localized source is presented. It is found that a small amount of phase dispersion can completely alter the behavior of an infiltrating NAPL plume, and that neglect of phase-dispersive processes may lead to unrealistic predictions of NAPL behavior in the vadose zone.     

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.     

On the Validity of Darcy's Law for Stable High-Concentration Displacements in Granular Porous Media

Recently, it has been suggested that Darcy's Law might not be applicable for modelling miscible, density-dependent flow in porous media. To investigate this, three sets of careful laboratory column experiments were performed on coarse and medium sands, consisting of upward displacement of water by sodium chloride solutions with concentrations ranging from 5 to 200g/l. Data on salt concentrations and water pressures were collected in horizontal transects along the flow direction. Salt concentration data were also collected in the influent and exit lines. The experimental data were analysed using a simplified approach based on Darcy's Law alone, applied with the assumption of a sharp interface. Darcy's Law was used to estimate porous medium permeability by fitting predictions to experimental data. Consistent estimates of permeability were obtained for each set of experiments. The results indicate that Darcy's Law adequately describes high concentration displacements through saturated coarse- and medium-grained porous media.