Estimation of Dynamic Relative Permeability and Capillary Pressure from Countercurrent Imbibition Experiments

Direct laboratory measurements of in situ water-phase saturation history are used to estimate relative permeability and capillary pressure functions. The magnitude of so-called nonequilibrium effects during spontaneous imbibition is quantified and, if significant, these effects are incorporated within the estimation technique. The primary constraint employed is that curves must increase or decrease monotonically; otherwise, no predetermined functionality is assumed. The technique is demonstrated using water saturation profile histories obtained for diatomite (a low-permeability and high-porosity rock). Results indicate that nonequilibrium effects detected at laboratory scale in low-permeability rocks influence the estimation of unsteady-state relative permeability and capillary pressure.     

A Mechanistic Model of Gas-Condensate Flow in Pores

Recent experimental results reported in the literature indicate that the relative permeability of gas-condensate systems increases with rate (velocity) at some conditions. To gain a better understanding of the nature of the flow and the prevailing mechanisms resulting in such behaviour flow visualisation experiments have been performed, using high pressure micromodels. The observed flow behaviour at the pore level has been employed to develop a mechanistic model describing the coupled flow of gas and condensate phases. The results of the model simulating the observed simultaneous flow of gas and condensate phases have been compared with reported core experimental results. Most features of the reported rate effect are predictable by the developed single pore model, nevertheless, its extension to include multiple pore interaction is recommended.     

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.     

近井地带凝析油聚集机理研究综述

凝析气藏衰竭开发过程中,当压力降到露点压力以下时,凝析油会在近井地带迅速聚集,严重损害凝析气井的生产能力.从实验和理论上综述凝析油在近井地带的聚集机理,通过研究凝析气藏衰竭开发过程中出现在近井地带的凝析油气流动的不同流型和相态特征,评价近井地带可动凝析油和不可动凝析油对气体流动性的影响.准确认识凝析油气的聚集特征和聚集机理对凝析气藏的高效开发具有重要的指导意义.      

Gas Condensate Relative Permeability for Well Calculations

This paper addresses several issues related to the modeling and experimental design of relative permeabilities used for simulating gas condensate well deliverability. Based on the properties of compositional flow equations, we make use of the fact that relative permeability ratio k _rg/k _ro is a purely thermodynamic variable, replacing saturation, when flow is steady-state. The key relation defining steady-state flow in gas condensate wells is relative permeability k _rg as a function of k _rg/k _ro. Consequently, determination of saturation and k _r as a function of saturation is not important for this specific calculation. Once the k _rg=f(k _rg/k _ro) relationship is experimentally established and correlated with capillary number (N _c), accurate modeling of condensate blockage is possible. A generalized model is developed for relative permeability as the function of k _rg/k _ro and capillary number. This model enables us to link the immiscible or rock curves with miscible or 'straight-line curves by a transition function dependent on the capillary number. This model is also extended to the case of high-rate, inertial gas flow within the steady-state condensate blockage regionand locally at the wellbore. We have paid particular attention to the effect of hysteresis on the relation k _rg=f(k _rg/k _ro), based on our observation that many repeated cycles of partial/complete imbibition and drainage occur in the near-well region during the life of a gas condensate well. Finally, the composite effect of condensate blockage is handled using a Muskat pseudopressure model, where relative permeabilities are corrected for the positive effect of capillary number dependence and the negative effect of inertial high velocity flow. Special steady-state experimental procedures have been developed to measure k _rg as a function of k _rg/k _ro and N _c. Saturations, though they can be measured, are not necessary. An approach for fitting steady-state gas condensate relative permeability data has been developed and used for modeling relative permeability curves.     

Effect of Network Topology on Relative Permeability

We consider the role of topology on drainage relative permeabilities derived from network models. We describe the topological properties of rock networks derived from a suite of tomographic images of Fontainbleau sandstone (Lindquist et al., 2000, J. Geophys. Res. 105B, 21508). All rock networks display a broad distribution of coordination number and the presence of long-range topological bonds. We show the importance of accurately reproducing sample topology when deriving relative permeability curves from the model networks. Comparisons between the relative permeability curves for the rock networks and those computed on a regular cubic lattice with identical geometric characteristics (pore and throat size distributions) show poor agreement. Relative permeabilities computed on regular lattices and on diluted lattices with a similar average coordination number to the rock networks also display poor agreement. We find that relative permeability curves computed on stochastic networks which honour the full coordination number distribution of the rock networks produce reasonable agreement with the rock networks. We show that random and regular lattices with the same coordination number distribution produce similar relative permeabilities and that the introduction of longer-range topological bonds has only a small effect. We show that relative permeabilities for networks exhibiting pore–throat size correlations and sizes up to the core-scale still exhibit a significant dependence on network topology. The results show the importance of incorporating realistic 3D topologies in network models for predicting multiphase flow properties.     

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.     

The Effect of Wettability on Three-Phase Relative Permeability

We study three-phase flow in water-wet, oil-wet, and fractionally-wet sandpacks. We use CT scanning to measure directly the oil and water relative permeabilites for three-phase gravity drainage. In an analogue experiment, we measure pressure gradients in the gas phase to determine the gas relative permeability. Thus we find all three relative permeabilities as a function of saturation. We find that the gas relative permeability is approximately half as much in a oil-wet medium than in an water-wet medium at the same gas saturation. The water relative permeability in the water-wet medium and the oil relative permeability in the oil-wet medium are similar. In the water-wet medium the oil relative permeability scales as k _roS _o ⁴ for S _o>S _or, where S _or is the waterflood residual oil saturation. With octane as the oil phase, k _roS _o ² for S _o<S _or, while with decane as the oil phase, k _ro falls sharply for S _o<S _or. The water relative permeability in the oil-wet medium resembles the oil relative permeability in the water-wet medium for a non-spreading oil such as decane. These observations can be explained in terms of wetting, spreading, and the pore scale configurations of fluid.     

Relative Permeability of Gas for Unconventional Reservoirs

Relative permeability of gas gains great significance in exploring unconventional gas. This paper developed a universal relative permeability model of gas, which is applicable for unconventional gas reservoirs such as coal, tight sandstone and shale. The model consists of the absolute relative permeability of gas and the gas slippage permeability. In the proposed model, the effects of water saturation and mean pore pressure on gas slippage permeability are taken into account. Subsequently, the evaluation of the model with existing model is done and then the validation of the model is made with data of tight sandstones, coals and shales from published literatures. The modeling results illustrate that a strong power-law relationship between relative permeability of gas and water saturation and the contribution of gas slippage permeability to relative permeability is determined by water saturation and mean pore pressure simultaneously. Furthermore, a sensitivity analysis of the impact of the parameters in the model is conducted and their effects are discussed.     

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.     

Relative permeability relations: A key factor for a drying model

In the modelling of heat, mass and momentum transfer phenomena which occur in a capillary porous medium during drying, the liquid and gas flows are usually described by the generalised Darcy laws. Nevertheless, the question of how to determine experimentally the relative permeability relations remains unanswered for most materials that consist of water and humid air, and as a result, arbitrary functions are used in the drying codes. In this paper, the emphasis is on deducing from both numerical and experimental studies a method for estimating pertinent relations for these key parameters. In the first part, the sensitivity of liquid velocity and, consequently, of drying kinetics in the variation of the relative permeabilities is investigated numerically by testing various forms. It is concluded that in order to predict a realistic liquid velocity behaviour, relative permeabilities can be linked to a measurable quantity: the capillary pressure. An estimation technique, based on simulations coupled with experimental measurements of capillary pressure, together with moisture content kinetics obtained for low or middle temperature convective drying, is deduced. In the second part, the proposed methodology is applied to pine wood. It is shown that the obtained relations provide closer representation of physical reality than those commonly used.     

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.     

多孔介质渗透率的一种新分形模型

多孔介质的渗流特性是油气藏工程、地下水资源利用、高放废物深地质处置等实际工程领域的热门研究问题．基于分形理论及多孔介质由一束面积大小不等的椭圆形毛细管组成的假设,本文建立了流体在分形多孔介质中渗流时的绝对渗透率及相对渗透率的分形渗透率模型．结果表明,绝对渗透率是最大和最小孔隙面积、分形维数、形状因子ε的函数,且当ε =1时,本文模型可以简化成Yu与Cheng模型;而非饱和多孔介质的相对渗透率与饱和度和多孔介质微结构参数有关．将本文提出的渗透率分形模型预测与实验测量数据及其他模型结果进行对比,显示它们整体吻合很好．     

Relative permeabilities of a near critical binary fluid

Relative permeabilities were measured at very low interfacial tensions (IFT) for two-phase mixtures of methanol and hexane flowing through Clashach sandstone. These two components pass from a two- to a single-phase system as the temperature is increased above the critical solution temperature (CST). The interfacial tension between the coexisting phases approaches zero as the solution reaches miscibility. The phase behaviour of methanol and hexane mixtures has been well characterised allowing the calculation of relative permeabilities, saturations and capillary numbers. Flow data are reported for four different temperatures in the two-phase region (i.e., four values of IFT and capillary number). The capillary desaturation curve (CDC) for the strongly wetting methanol rich phase is also presented. In addition to the novel technique presented for measurement of relative permeability, the results indicate that relative permeabilities approach straight line functions very near the critical point. Furthermore, desaturation of the wetting phase was found to be dependent on the capillary number which, in turn, depends on the location of the mixture on the fluid phase diagram and the proximity to the critical temperature.     

General Anisotropic Effective Medium Theory for the Effective Permeability of Heterogeneous Reservoirs

One of the techniques to calculate the effective property of a heterogeneous medium is the effective medium theory. The present paper presents a general mathematical formulation for the effective medium approximation using a self-consistent choice of the effective permeability, to apply it to the case of a general anisotropic 2D medium and to the case of a 3D isotropic medium with randomly oriented ellipsoidal inclusions. The 2D results are compared with analytical results and with a homogenization technique with good result. The 3D correlations are used to derive percolation thresholds in two-phase systems with a large permeability contrast, which are compared to numerical results from the literature, also with good results.     

Effects of Relative Permeability History Dependence on Two-Phase Flow in Porous Media

Hysteresis phenomena in multi-phase flow in porous media has been recognized by many researchers and widely believed to have significant effects on the flow. In an attempt to account for these effects, a theoretical model for history-dependent relative permeabilities is considered. This model is incorporated into 1-D two-phase nondiffusive flow system and the corresponding flow is predicted. Flow history is observed to have a notable impact on the saturation profile and fluids breakthrough.     

Microvisual Study of Multiphase Gas Condensate Flow in Porous Media

Gas condensate reservoirs constitute a significant portion of hydrocarbon reserves worldwide. The liquid drop-out in these reservoirs may lead to recovery problems such as near wellbore permeability impairment and uncertainty in the actual location of the target condensate. Such technical issues can be addressed through improved understanding of the formation of condensate and the multiphase flow of gas/condensate/water in the reservoir as characterized by relative permeability curves. The appropriate relative permeability curves in turn can be used in reservoir simulators to assist in optimization of field development. This paper reports results of experiments conducted in micromodels, in support of possible core flow tests, using reservoir fluids under reservoir conditions. In particular, visualizations of condensate formation with and without connate water are presented and the differences between the two cases as well as the possible implications for the relative permeability measurements are discussed. Furthermore, the flow of gas and condensate at different force ratios (capillary and Bond numbers) are presented. It is postulated that a single dimensionless number may not be sufficient to characterize the multiphase flow in gas condensate reservoirs. The physical mechanisms occurring under various field conditions are examined in the light of these observations.     

Pore-Level Modeling of Gas and Condensate Flow in Two- and Three-Dimensional Pore Networks: Pore Size Distribution Effects on the Relative Permeability of Gas and Condensate

We present a mechanistic model of retrograde condensation processes in two- and three-dimensional capillary tube networks under gravitational forces. Condensate filling-emptying cycles in pore segments and gas connection–isolation cycles are included. With the pore-level distribution of gas and condensate in hand, we determine their corresponding relative permeabilities. Details of pore space and displacement are subsumed in pore conductances. Solving for the pressure field in each phase, we find a single effective conductance for each phase as a function of condensate saturation. Along with the effective conductance for the saturated network, the relative permeability for each phase is calculated. Our model porous media are two- and three-dimensional regular networks of pore segments with distributed size and square cross-section. With a Monte Carlo sampling we find the optimum network size to avoid size effects and then we investigate the effect of network dimensionality and pore size distribution on the relative permeabilities of gas and condensate.     

Effective Relative Permeabilities and Capillary Pressure for One-Dimensional Heterogeneous Media

The paper presents an analytical construction of effective two-phase parameters for one-dimensional heterogeneous porous media, and studies their properties. We base the computation of effective parameters on analytical solutions for steady-state saturation distributions. Special care has to be taken with respect to saturation and pressure discontinuities at the interface between different rocks. The ensuing effective relative permeabilities and effective capillary pressure will be functions of rate, flow direction, fluid viscosities, and spatial scale of the heterogeneities.The applicability of the effective parameters in dynamic displacement situations is studied by comparing fine-gridded simulations in heterogeneous media with simulations in their homogeneous (effective) counterparts. Performance is quite satisfactory, even with strong fronts present. Also, we report computations studying the applicability of capillary limit parameters outside the strict limit.     

Measurement and Modelling of Gas Condensate Flow Around Rock Perforation

In the petroleum production industry wells are mostly cased and perforated in the producing formation. Perforation characteristics such as size, length, number of perforation tunnels and their arrangements as well as fluid and rock properties determine fluid flow behaviour in the wellbore region, hence, well productivity. Flow of gas and condensate around a perforation tunnel (including the damaged zone) has been studied by performing steady-state core experiments and simulating the results numerically, using a finite element modelling approach. The model allows for the changes in fluid properties and accounts for the coupling of the two phases and the inertial effect using a fractional flow based correlation. The results indicated that different sets of thickness-permeability (h−k) values obtained from matching single-phase flow performance could be assigned to the damaged zone around perforation to represent the two-phase flow performance. The status of the tip of the perforation for two extreme cases of totally closed and fully open was investigated and found to have a minimal effect on the performance of the system.