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.     

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.     

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.     

Simultaneous Estimation of Relative Permeability and Capillary Pressure Using Ensemble-Based History Matching Techniques

An ensemble-based technique has been developed and successfully applied to simultaneously estimate the relative permeability and capillary pressure by history matching the observed production profile. Relative permeability and capillary pressure curves are represented by using a power-law model. Then, forward simulation is performed with the initial coefficients of the power-law model, all of which are to be tuned automatically and finally determined once the observed data is assimilated completely and history matched. The newly developed technique has been validated by a synthetic coreflooding experiment with two scenarios. The endpoints are fixed for the first scenario, whereas they are completely free in the second scenario. Simultaneous estimation of relative permeability and capillary pressure has been found to improve gradually as more observation data is assimilated. There exists an excellent agreement between both the updated relative permeability and capillary pressure and their corresponding reference values, once the discrepancy between the simulated and observed production history has been minimized. Compared with coefficients of capillary pressure curve, coefficients of relative permeability curves, irreducible water saturation and residual oil saturation are found to be more sensitive to the observed data. In addition, water relative permeability is more sensitive to the observation data than either oil relative permeability or capillary pressure. It is shown from its application to a laboratory coreflooding experiment that relative permeability and capillary pressure curves can be simultaneously evaluated once all of the experimental measurements are assimilated and history matched.     

Insights into the Relationships Among Capillary Pressure,Saturation, Interfacial Area and Relative Permeability Using Pore-Network Modeling

To gain insight in relationships among capillary pressure, interfacial area, saturation, and relative permeability in two-phase flow in porous media, we have developed two types of pore-network models. The first one, called tube model, has only one element type, namely pore throats. The second one is a sphere-and-tube model with both pore bodies and pore throats. We have shown that the two models produce distinctly different curves for capillary pressure and relative permeability. In particular, we find that the tube model cannot reproduce hysteresis. We have investigated some basic issues such as effect of network size, network dimension, and different trapping assumptions in the two networks. We have also obtained curves of fluid–fluid interfacial area versus saturation. We show that the trend of relationship between interfacial area and saturation is largely influenced by trapping assumptions. Through simulating primary and scanning drainage and imbibition cycles, we have generated two surfaces fitted to capillary pressure, saturation, and interfacial area (P _c –S _w –a _nw ) points as well as to relative permeability, saturation, and interfacial area (k _r –S _w –a _nw ) points. The two fitted three-dimensional surfaces show very good correlation with the data points. We have fitted two different surfaces to P _c –S _w –a _nw points for drainage and imbibition separately. The two surfaces do not completely coincide. But, their mean absolute difference decreases with increasing overlap in the statistical distributions of pore bodies and pore throats. We have shown that interfacial area can be considered as an essential variable for diminishing or eliminating the hysteresis observed in capillary pressure–saturation (P _c –S _w ) and the relative permeability–saturation (k _r –S _w ) curves.     

Determination of Capillary Pressure Function from Resistivity Data

A model has been derived theoretically to correlate capillary pressure and resistivity index based on the fractal scaling theory. The model is simple and predicts a power law relationship between capillary pressure and resistivity index (P _c = p _e · I ^β) in a specific range of low water saturation. To verify the model, gas-water capillary pressure and resistivity were measured simultaneously at a room temperature in 14 core samples from two formations in an oil reservoir. The permeability of the core samples ranged from 0.028 to over 3000 md. The porosity ranged from less than 8 to over 30. Capillary pressure curves were measured using a semi-permeable porous-plate technique. The model was tested against the experimental data obtained in this study. The results demonstrated that the model could match the experimental data in a specific range of low water saturation. The experimental results also support the fractal scaling theory in a low water saturation range. The new model developed in this study may be deployed to determine capillary pressure from resistivity data both in laboratories and reservoirs, especially in the case in which permeability is low or it is difficult to measure capillary pressure.     

Pore Network Modelling of Electrical Resistivity and Capillary Pressure Characteristics

In order to model petrophysical properties of hydrocarbon reservoir rocks, the underlying physics occurring in realistic rock pore structures must be captured. Experimental evidence showing variations of wetting occurring within a pore, and existence of the so-called 'non-Archie' behaviour, has led to numerical models using pore shapes with crevices (for example, square, elliptic, star-like shapes, etc.). This paper presents theoretical derivations and simulation results of a new pore space network model for the prediction of petrophysical properties of reservoir rocks. The effects of key pore geometrical factors such as pore shape, pore size distribution and pore co-ordination number (pore connectivity) have been incorporated into the theoretical model. In particular, the model is used to investigate the effects of wettability and saturation history on electrical resistivity and capillary pressure characteristics. The petrophysical characteristics were simulated for reservoir rock samples. The use of the more realistic grain boundary pore (GBP) shape allows simulation of the generic behaviour of sandstone rocks, with various wetting scenarios. The predictions are in close agreement with electrical resistivity and capillary pressure characteristics observed in experiments.     

Grain Shape Effects on Permeability,Formation Factor,and Capillary Pressure from Pore-Scale Modeling

We invoke pore-scale models to evaluate grain shape effects on petrophysical properties of three-dimensional (3D) images from micro-CT scans and consolidated grain packs. Four sets of grain-packs are constructed on the basis of a new sedimentary algorithm with the following shapes: exact angular grain shapes identified from micro-CT scans, ellipsoids fitted to angular grains, and spheres with volume and surface-to-volume ratio equal to original angular grains on a grain-by-grain basis. Subsequently, a geometry-based cementation algorithm implements pore space alteration due to diagenesis. Eight micro-CT scans and 144 grain-pack images with $500 \times 500 \times 500$ voxels (the resolution units of 3D images) are analyzed in this study. Absolute permeability, formation factor, and capillary pressure are calculated for each 3D image using numerical methods and compared to available core measurements. Angular grain packs give rise to the best agreement with experimental measurements. Cement volume and its spatial distribution in the pore space significantly affect all calculated petrophysical properties. Available empirical permeability correlations for non-spherical grains underestimate permeability between 30 and 70 % for the analyzed samples. Kozeny–Carman’s predictions agree with modeled permeability for spherical grain packs but overestimate permeability for micro-CT images and non-spherical grain packs when volume-based radii are used to calculate the average grain size in a pack. We identify surface-to-volume ratio and grain shape as fundamental physical parameters that control fluid distribution and flow in porous media for equivalent porosity samples.     

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.     

Steady-State Relative Permeability Measurements of Tight and Shale Rocks Considering Capillary End Effect

Transport in Porous Media - Relative permeability (kr) data are the key factors for describing the behaviour of the multi-phase flow in porous media. During the kr measurements of low-permeability...     

A Markov Chain Model for Effective Relative Permeabilities and Capillary Pressure

We present a new method for calculating the effective two-phase parameters of one-dimensional randomly heterogeneous porous media, which avoids the timeconsuming use of simulations on explicit realizations. The procedure is based on the steady state saturation distribution. The idea is to model the local variation of saturation and saturation dependent parameters as Markov chains, in such a way that the effective parameters are given by the asymptotic expectations of the chains. We derive the exact asymptotic moment equations and solve them numerically, based on their second order approximation. The method determines the effective parameters to a high degree of accuracy, even with large variations in rock properties. In particular, the capillary limit and viscous limit effective parameters are recovered exactly. The applicability of the effective parameters in the unsteady state case is studied by comparing the displacement production profiles in heterogeneous media and their homogenized counterpart.     

A New Method for Calculating Two-Phase Relative Permeability from Resistivity Data in Porous Media

Many resistivity data from laboratory measurements and well logging are available. Papers on the relationship between resistivity and relative permeability have been few. To this end, a new method was developed to infer two-phase relative permeability from the resistivity data in a consolidated porous medium. It was found that the wetting phase relative permeability is inversely proportional to the resistivity index of a porous medium. The proposed model was verified using the experimental data in different rocks (Berea, Boise sandstone, and limestone) at different temperatures up to 300°F. The results demonstrated that the oil and water relative permeabilities calculated from the experimental resistivity data by using the model proposed in this article were close to those calculated from the capillary pressure data in the rock samples with different porosities and permeabilities. The results demonstrated that the proposed approach to calculating two-phase relative permeability from resistivity data works satisfactorily in the cases studied.     

Predicting Resistivity and Permeability of Porous Media Using Minkowski Functionals

Transport in Porous Media - Permeability and formation factor are important properties of a porous medium that only depend on pore space geometry, and it has been proposed that these transport...     

Collagen Fiber Network Infiltration: Permeability and Capillary Infiltration

In dentin restoration, collagen fiber network infiltration is an issue. Using data from the literature, we have constructed a relevant numerical geometrical model of the network. The specificity of our model is that the fibers are taken into account implicitly using a regularized Heaviside function. This function is either used to set the viscosity or to localize the contact line where capillary forces are applied. A level set technique with respect to fluid infiltration front tracking in five fiber networks using the level set method and Navier–Stokes equations with capillary terms is used to point out efficient critical infiltration parameters. A variational formulation which can be implemented in a finite elements model is proposed both for the infiltration front and the contact line. Because of lack of knowledge on fiber orientation, different configurations were tested through permeability assessment of the whole network. Fiber orientation, interfibrillar space and contact angle influence were investigated.     

Dynamic Relative Permeability and Simulation of WAG Injection Processes

In this work, we investigate the impact of mobility changes due to flow reversals from co-current to counter-current flow on the displacement performance of water alternating gas (WAG) injection processes. In WAG processes, the injected gas will migrate toward the top of the formation while the injected water will migrate toward the bottom of the formation. The segregation of gas, oil and water phases will result in counter-current flow occurring in the vertical direction in some portions of the reservoir during the displacement process. Previous experimental and theoretical studies of counter-current flow have shown that the relative mobility of each of the phases in a porous medium is considerably less when counter-current flow prevails as compared to co-current flow settings. A reduction of the relative permeability in the vertical direction results in a dynamic anisotropy in phase mobilities. This effect has, to the best of our knowledge, not previously been considered in the modeling and simulation of WAG processes. A new flow model that accounts for flow reversals in the vertical direction has been implemented and tested in a three-phase compositional reservoir simulator. In order to investigate the impact of flow reversals, results from the new flow model are compared to cases where counter-current flow effects on the phase mobilities are ignored. A range of displacement settings, covering relevant slug sizes, have been investigated to gauge the impact of mobility reductions due to flow reversals. Significant differences, in terms of saturation distribution, producing GOR and oil recovery, are observed between the conventional flow model (ignoring mobility reductions due to counter-current flow) and the proposed new model that accounts for reductions in phase mobility during counter-current flow. Accordingly, we recommend that an explicit representation of flow transitions between co-current and counter-current flow (and the related impact on phase mobilities) should be considered to ensure accurate and optimal design of WAG injection processes.     

Effects of Organic Base Chemistry on Interfacial Tension,Wettability, and Capillary Pressure in Multiphase Subsurface Waste Systems

The presence of surfactants may have profound effects on the transport of organic contaminants in multiphase systems. It is a common practice, however, to model the subsurface migration of liquids independently of the aqueous phase composition. As such, transport in these systems may not be adequately characterized. This study investigates the impact of pH on interfacial tension, wettability, and the drainage capillary pressure–saturation relationship in air–water–quartz and oxylene–water–quartz systems containing dodecylamine, an organic base. In these systems, three mechanisms, speciation, partitioning, and sorption, are important in determining the interfacial tension and contact angle, and consequently, important in determining the capillary pressure. By adjusting the pH above and below, the pK_a of the base, the relative importance of these mechanisms was altered. Below dodecylamine's pK_a of 10.6, the base was primarily in a cationic form resulting in minimal partitioning into the nonaqueous liquid and greater sorption at the quartz surface. Above the pK_a, the base was primarily in a neutral form which did not sorb to the quartz, and, furthermore, partitioned into the organic liquid phase where its surface activity was minimized. The combination of these processes caused the capillary pressure to change in a manner consistent with porescale theory of capillarity. The utility in this approach lies in the possibility of predicting transport properties in multiphase systems while incorporating the direct effects of solution chemistry.     

Impact of Capillary Pressure, Salinity and In situ Conditions on CO2 Injection into Saline Aquifers

This article presents a numerical investigation of the combined effects of capillary pressure, salinity and in situ thermodynamic conditions on CO₂-brine-rock interactions in a saline aquifer. We demonstrate that the interrelations between capillary pressure, salinity, dissolution and drying-out affect CO₂ injectivity and storage capacity of a saline aquifer. High capillary forces require a high injection pressure for a given injection rate. Depending on salinity, the increase in injection pressure due to capillary forces can be offset by the dissolution of CO₂ in formation water and its compressibility. Higher capillary forces also reduce gravity segregation, and this gives a more homogeneous CO₂ plume which improves the dissolution of CO₂. The solubility of CO₂ in formation water decreases with increasing salinity which requires an increased injection pressure. Higher salinity and capillary pressure can even block the pores, causing an increased salt precipitation. Simulations with various pressure-temperature conditions and modified salinity and capillary pressure curves demonstrate that, with the injection pressures similar for both cold and warm basins at a given injection rate, CO₂ dissolves about 10% more in the warm basin water than in the cold basin. The increase in dissolution lowers the injection pressure compensating the disadvantage of low CO₂ density and compressibility for storage in warm basins.     

Pore Shape Effects on the Relative Permeability of Gas and Condensate in Three-Dimensional Pore Networks

Transport in Porous Media -     

A Numerical Study of Two-Phase Flow Models with Dynamic Capillary Pressure and Hysteresis

Saturation overshoot and pressure overshoot are studied by incorporating dynamic capillary pressure, capillary pressure hysteresis and hysteretic dynamic coefficient with a traditional fractional flow equation in one-dimensional space. Using the method of lines, the discretizations are constructed by applying the Castillo–Grone’s mimetic operators in the space direction and a semi-implicit integrator in the time direction. Convergence tests and conservation properties of the schemes are presented. Computed profiles capture both the saturation overshoot and pressure overshoot phenomena. Comparisons between numerical results and experiments illustrate the effectiveness and different features of the models.     

微梁与基底间的毛细粘附作用

在表面微型机械结构的制造过程中,强的毛细相互作用常常使得组成这些结构的微桥、微梁与基底粘附而导致失效。而在微尺度实验中,微桥与微梁又是微尺度材料常数和性能检测的常用的试件样式,如果实验中加载端与被检测的微尺度试件发生毛细粘附,将直接影响检测数据的准确性。本文应用微悬臂梁试件,讨论微梁与基底间的毛细粘附作用,并通过能量原理计算其粘附力的大小和试件几何尺寸、粘附面距离、粘附液体特性之间的关系。最后应用微散斑干涉,检测粘附平衡态时微桥和微梁的粘附力以及由毛细粘附所导致的弯曲变形,并与理论计算结果进行比较。