Modeling Non-Darcian Single- and Two-Phase Flow in Transparent Replicas of Rough-Walled Rock Fractures

While it is generally assumed that in the viscous flow regime, the two-phase flow relative permeabilities in fractured and porous media depend uniquely on the phase saturations, several studies have shown that for non-Darcian flows (i.e., where the inertial forces are not negligible compared with the viscous forces), the relative permeabilities not only depend on phase saturations but also on the flow regime. Experimental results on inertial single- and two-phase flows in two transparent replicas of real rough fractures are presented and modeled combining a generalization of the single-phase flow Darcy’s law with the apparent permeability concept. The experimental setup was designed to measure injected fluid flow rates, pressure drop within the fracture, and fluid saturation by image processing. For both fractures, single-phase flow experiments were modeled by means of the full cubic inertial law which allowed the determination of the intrinsic hydrodynamic parameters. Using these parameters, the apparent permeability of each fracture was calculated as a function of the Reynolds number, leading to an elegant means to compare the two fractures in terms of hydraulic behavior versus flow regime. Also, a method for determining the experimental transition flow rate between the weak inertia and the strong inertia flow regimes is proposed. Two-phase flow experiments consisted in measuring the pressure drop and the fluid saturation within the fractures, for various constant values of the liquid flow rate and for increasing values of the gas flow rate. Regardless of the explored flow regime, two-phase flow relative permeabilities were calculated as the ratio of the single phase flow pressure drop per unit length divided by the two-phase flow pressure drop per unit length, and were plotted versus the measured fluid saturation. Results confirm the dependence of the relative permeabilities on the flow regime. Also the proposed generalization of Darcy’s law shows that the relative permeabilities versus fluid saturation follow physical meaningful trends for different liquid and gas flow rates. The presented model fits correctly the liquid and gas experimental relative permeabilities as well as the fluid saturation.     

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.     

Interrelationship between Resistivity Index,Capillary Pressure and Relative Permeability

It is known that the three important parameters, resistivity, capillary pressure, and relative permeability, are all a function of fluid saturation in a porous medium. This implies that there may be a correlation among the three parameters. There have been many papers on the approach to inferring relative permeability from capillary pressure data. However, the literature on the interrelationship between resistivity index, capillary pressure, and relative permeability has been few. The models representing such relationships have been proposed in this study, including a new model correlating relative permeability and capillary pressure. Some of the models were verified using experimental data for the first time. It has been shown that the other two parameters could be determined using these models if one of the three parameters (capillary pressure, relative permeability, and resistivity) is known. Using this approach, it would be possible to quickly obtain a distribution of capillary pressure and relative permeability characteristics as a function of depth and location across an entire reservoir.     

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.     

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.     

A Multiscale Network Model for Simulating Moisture Transfer Properties of Porous Media

A multiscale network model is presented to model unsaturated moisture transfer in hygroscopic capillary-porous materials showing a broad pore-size distribution. Both capillary effects and water sorption phenomena, water vapour and liquid water transfer are considered. The multiscale approach is based on the concept of examining the porous space at different levels of magnification. The conservation of the water vapour permeability of dry material is used as scaling criterion to link the different pore scales. A macroscopic permeability is deduced from the permeabilities calculated at the different levels of magnification. Each level of magnification is modelled using an isotropic nonplanar 2D cross-squared network. The multiscale network simulates the enhancement of water vapour permeability due to capillary condensation, the hysteresis phenomenon between wetting and drying, and the steep increase of moisture permeability at the critical moisture saturation level. The calculated network permeabilities are compared with experimental data for calcium silicate and ceramic brick and a good agreement is observed.     

A Simple Model for the Variation of Permeability due to Partial Saturation in Dual Scale Porous Media

The main focus of this work is to model macroscopically the effects of partial saturation upon the permeability of dual scale fibrous media made of fiber bundles when a Newtonian viscous fluid impregnates it. A new phenomenological model is proposed to explain the discrepancies between experimental pressure results and analytical predictions based on Darcy's law. This model incorporates the essential features of relative permeability but without the necessity of measuring saturation of the liquid for its prediction. The model is very relevant for the small scale industrial systems where a liquid is forced to flow through a fibrous porous medium. It requires four parameters. Two of them are the two permeability values based on the two length scales. One length scale is of the order of magnitude of the individual fiber radius and corresponds to the permeability of the completely staurated medium, the other is of the order of magnitude of the distance between the fiber bundles and corresponds to the permeability of the partially saturated medium. The other two parameters are the lengths of the two partially saturated regions of the flow domain. The two lengths of the partially saturated region and the permeability of the fully saturated flow domain can be directly measured from the experiments. The excellent agreement between the model and the experimental results of inlet pressure profile with respect to time suggests that this model may be used to describe the variation of the permeability behind a moving front in such porous media for correct pressure prediction. It may also be used to characterize the fibrous medium by determining the two different permeabilities and the relative importance of the unsaturated portion of the flow domain for a given architecture.     

Two-Phase Filtration in a Transversally Isotropic Porous Medium: Experiment and Theory

The results of an experimental determination of the relative phase permeabilities during flow of two immiscible fluids in stratified sandstone with transversally isotropic characteristics are presented. The measurements were performed on samples oriented in three directions: along, perpendicular to and at an angle of 45° to the stratification plane. An approximate solution of the problem of steady two-phase flow toward a finite gallery in an anisotropic porous medium for arbitrary relative orientation of the gallery and the principal axes of the absolute permeability tensor is given. This solution was tested against the experimental results. The good agreement between the theoretical and experimental results makes it possible to recommend for engineering calculations both the relations between the absolute and phase permeabilities for transversally isotropic and orthotropic porous media and the approximate solution obtained.     

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.     

An Air Permeability Study of Anisotropic Glass Wool Fibrous Products

This article reports an experimental study of the permeability of fibrous porous media based on air flow. Two glass wools have been analyzed in three different flow directions to evaluate their anisotropic structure. The usual flow theories such as Darcy’s law, valid for granular materials, have been shown to apply also to these fibrous media which is a more complex structure (fiber arrangement, high porous media,…). The intrinsic permeabilities and the anisotropic factors have been determined from an air permeability measurement based on a standard test for industrial fabrics. The experimental results have been compared to different models, showing the limits of some simplistic models when applied to our highly porous structures. One model, Kyan’s model based on the flow around submerged objects, was, however, found to agree with our experimental data. This study also demonstrates the interesting use of an experimental technique, which is air permeability measurement (water permeability is not adapted), in the investigation of flow phenomena in a fibrous medium.     

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.     

Origin and quantification of coupling between relative permeabilities for two-phase flows in porous media

An extended formulation of Darcy's two-phase law is developed on the basis of Stokes' equations. It leads, through results borrowed from the thermodynamics of irreversible processes, to a matrix of relative permeabilities. Nondiagonal coefficients of this matrix are due to the viscous coupling exerted between fluid phases, while diagonal coefficients represent the contribution of both fluid phases to the total flow, as if they were alone. The coefficients of this matrix, contrary to standard relative permeabilities, do not depend on the boundary conditions imposed on two-phase flow in porous media, such as the flow rate. This formalism is validated by comparison with experimental results from tests of two-phase flow in a square cross-section capillary tube and in porous media. Coupling terms of the matrix are found to be nonnegligible compared to diagonal terms. Relationships between standard relative permeabilities and matrix coefficients are studied and lead to an experimental way to determine the new terms for two-phase flow in porous media.     

Relative Permeability Analysis of Tube Bundle Models

The analytical solution for calculating two-phase immiscible flow through a bundle of parallel capillary tubes of uniform diametral probability distribution is developed and employed to calculate the relative permeabilities of both phases. Also, expressions for calculating two-phase flow through bundles of serial tubes (tubes in which the diameter varies along the direction of flow) are obtained and utilized to study relative permeability characteristics using a lognormal tube diameter distribution. The effect of viscosity ratio on conventional relative permeability was investigated and it was found to have a significant effect for both the parallel and serial tube models. General agreement was observed between trends of relative permeability ratios found in this work and those from experimental results of Singhal et al. (1976) using porous media consisting of mixtures of Teflon powder and glass beads. It was concluded that neglecting the difference between the average pressure of the non-wetting phase and the average pressure of the wetting phase (the macro-scale capillary pressure) – a necessary assumption underlying the popular analysis methods of Johnson et al. (1959) and Jones and Roszelle (1978) – was responsible for the disparity in the relative permeability curves for various viscosity ratios. The methods therefore do not account for non-local viscous effects when applied to tube bundle models. It was contended that average pressure differences between two immiscible phases can arise from either capillary interfaces (micro-scale capillary pressures) or due to disparate pressure gradients that are maintained for a flow of two fluids of viscosity ratio that is different from unity.     

An Experimental Study of the Influence of Interfacial Tension on Water–Oil Two-Phase Relative Permeability

Adding surfactant into the displacing aqueous phase during surfactant-enhanced aquifer remediation of NAPL contamination and in chemical flooding oil recovery significantly changes interfacial tension (IFT) (σ) on water–oil interfaces within porous media. The change in IFT may have a large impact on relative permeability for the two-phase flow system. In most subsurface flow investigations, however, the influence of IFT on relative permeability has been ignored. In this article, we present an experimental study of two-phase relative- permeability behavior in the low and more realistic ranges of IFT for water–oil systems. The experimental work overcomes the limitations of the existing laboratory measurements of relative permeability (which are applicable only for high ranges of IFT (e.g., σ > 10⁻² mN/m). In particular, we have (1) developed an improved steady-state method of measuring complete water–oil relative permeability curves; (2) proven that a certain critical range of IFT exists such that IFT has little impact on relative permeability for σ greater than this range, while within the range, relative permeabilities to both water and oil phases will increase with decreasing IFT; and (3) shown that a functional correlation exists between water–oil two-phase relative permeability and IFT. In addition, this work presents such correlation formula between water–oil two-phase relative permeability and IFT. The experimental results and proposed conceptual models will be useful for quantitative studies of surfactant-enhanced aquifer remediation and chemical flooding operations in reservoirs.     

Network model evaluation of permeability and spatial correlation in a real random sphere packing

In principle, network models can replicate exactly the microstructure of porous media. In practice, however, network models have been constructed using various assumptions concerning pore structure. This paper presents a network model of a real, disordered porous medium that invokes no assumptions regarding pore structure. The calculated permeability of the model agrees well with measured permeabilities, providing a new and more rigorous confirmation of the validity of the network approach. Several assumptions commonly used in constructing network models are found to be invalid for a random packing of equal spheres. In addition, the model permits quantification of the effect of pore-scale correlation (departure from randomness) upon permeability. The effect is comparable to reported discrepancies between measured permeabilities and predictions of other network models. The implications of this finding are twofold. First, a key assumption of several theories of transport in porous media, namely that pore dimensions are randomly distributed upon a network, may be invalid for real porous systems. Second, efforts both to model and to measure pore-scale correlations could yield more accurate predictions of permeability.     

Determining equations of two-phase flows through anisotropic porous media

A new interpretation of the concept of relative phase permeability is given. Relative phase permeabilities are represented in the form of fourth-rank tensors. It is shown that in the case of anisotropic porous media functions depending not only on the saturation but also on the anisotropy parameters represented in the form of ratios of the principal values of the absolute permeability coefficient tensor correspond to the classical representation of the relative phase permeabilities. For a two-phase flow in anisotropic porous media with orthotropic and transversely-isotropic symmetry a generalized two-term Darcy’s law is analyzed. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 87–94, March–April, 1998. The work was carried out with support from the Russian Foundation for Fundamental Research (project No. 96-01-00623).     

On the use of conventional cocurrent and countercurrent effective permeabilities to estimate the four generalized permeability coefficients which arise in coupled,two-phase flow

In the case of coupled, two-phase flow of fluids in porous media, the governing equations show that there are four independent generalized permeability coefficients which have to be measured separately. In order to specify these four coefficients at a specific saturation, it is necessary to conduct two types of flow experiments. The two types of flow experiments used in this study are cocurrent and countercurrent, steady-state permeability experiments. It is shown that, by taking this approach, it is possible to define the four generalized permeability coefficients in terms of the conventional cocurrent and countercurrent effective permeabilities for each phase. It is demonstrated that a given generalized phase permeability falls about midway between the conventional, cocurrent effective permeability for that phase, and that for the countercurrent flow of the same phase. Moreover, it is suggested that the conventional effective permeability for a given phase can be interpreted as arising out of the effects of two types of viscous drag: that due to the flow of a given phase over the solid surfaces in the porous medium and that due to momentum transfer across the phase 1-phase 2 interfaces in the porous medium. The magnitude of the viscous coupling is significant, contributing at least 15% to the total conventional cocurrent effective permeability for both phases. Finally, it is shown that the nontraditional generalized permeabilities which arise out of viscous coupling effects cannot equal one another, even when the viscosity ratio is unity and the surface tension is zero.