Immiscible Displacement in the Interacting Capillary Bundle Model Part II. Applications of Model and Comparison of Interacting and Non-Interacting Capillary Bundle Models

In the first part of this work (Dong et al., Transport Porous Media, 59, 1–18, 2005), an interacting capillary bundle model was developed for analysing immiscible displacement processes in porous media. In this paper, the second part of the work, the model is applied to analyse the fluid dynamics of immiscible displacements. The analysis includes: (1) free spontaneous imbibition, (2) the effects of injection rate and oil–water viscosity ratio on the displacement interface profile, and (3) the effect of oil–water viscosity ratio on the relative permeability curves. Analysis of a non-interacting tube bundle model is also presented for comparison. Because pressure equilibration between the capillaries is stipulated in the interacting capillary model, it is able to reproduce the behaviour of immiscible displacement observed in porous media which cannot be modelled by using non-interacting tube bundle models.     

Influence of Relative Permeability on the Stability Characteristics of Immiscible Flow in Porous Media

Relative permeability functions for immiscible displacements in porous media show a wide range of profiles. Although, this behavior is well known, its impact on the stability of the displacement process is unexplored. Our analysis clearly demonstrates for the first time that the viscous instability characteristics of two-phase flows are governed not only by their end point values, but are strongly dependent on the actual profile of relative permeability functions. Linear stability analysis predicts the capacity of the flow to develop large scale fingers which can result in substantial bypassing of the resident fluid. It is observed that relative permeability functions attributed to drainage processes yield a more unstable displacement as compared to functions related to imbibition processes. Moreover, instability is observed to increase for those relative permeability functions which result from increased wettability of the wetting fluid. High accuracy numerical simulations show agreement with these predictions and demonstrate how large amplitude viscous fingers result in significant bypassing for certain relative permeability functions. In the nonlinear regime, the finger amplitude grows at a rate ∝ t^1/2 initially, drops to t^1/4 at a later time and finally grows ∝ t. The basic mechanisms of finger interaction, however, are not substantially influenced by relative permeability functions.     

Nonequilibrium Effects During Spontaneous Imbibition

Accurate models of multiphase flow in porous media and predictions of oil recovery require a thorough understanding of the physics of fluid flow. Current simulators assume, generally, that local capillary equilibrium is reached instantaneously during any flow mode. Consequently, capillary pressure and relative permeability curves are functions solely of water saturation. In the case of imbibition, the assumption of instantaneous local capillary equilibrium allows the balance equations to be cast in the form of a self-similar, diffusion-like problem. Li et al. [J. Petrol. Sci. Eng. 39(3) (2003), 309–326] analyzed oil production data from spontaneous countercurrent imbibition experiments and inferred that they observed the self-similar behavior expected from the mathematical equations. Others (Barenblatt et al. [Soc. Petrol. Eng. J. 8(4) (2002), 409–416]; Silin and Patzek [Transport in Porous Media 54 (2004), 297–322]) assert that local equilibirum is not reached in porous media during spontaneous imbibition and nonequilibirium effects should be taken into account. Simulations and definitive experiments are conducted at core scale in this work to reveal unequivocally nonequilbirium effects. Experimental in-situ saturation data obtained with a computerized tomography scanner illustrate significant deviation from the numerical local-equilibrium based results. The data indicates: (i) capillary imbibition is an inherently nonequilibrium process and (ii) the traditional, multi-phase, reservoir simulation equations may not well represent the true physics of the process.     

Percolation as a Basic Concept for Macroscopic Capillarity

The concepts of relative permeability and capillary pressure are crucial for the accepted traditional theory of two phase flow in porous media. Recently, a theoretical approach was introduced that does not require these concepts as input (Hilfer, Physica A, 359:119–128, 2006a; Phys. Rev. E, 73:016307, 2006b). Instead it was based on the concept of hydraulic percolation of fluid phases. This paper presents the first numerical solutions of the coupled nonlinear partial differential equations introduced in Hilfer (Phys. Rev. E, 73:016307, 2006b). Approximate numerical results for saturation profiles in one spatial dimension have been calculated. Long time limits of dynamic time-dependent profiles are compared to stationary solutions of the traditional theory. The initial and boundary conditions are chosen to model the displacement processes that occur when a closed porous column containing two immiscible fluids of different density is raised from a horizontal to a vertical position in a gravitational field. The nature of the displacement process may change locally in space and time between drainage and imbibition. The theory gives local saturations for nonpercolating trapped fluids near the endpoint of the displacement.     

A Dynamic Network Model for Two-Phase Flow in Porous Media

We present a dynamic model of immiscible two-phase flow in a network representation of a porous medium. The model is based on the governing equations describing two-phase flow in porous media, and can handle both drainage, imbibition, and steady-state displacement. Dynamic wetting layers in corners of the pore space are incorporated, with focus on modeling resistivity measurements on saturated rocks at different capillary numbers. The flow simulations are performed on a realistic network of a sandpack which is perfectly water-wet. Our numerical results show saturation profiles for imbibition in agreement with experiments. For free spontaneous imbibition we find that the imbibition rate follows the Washburn relation, i.e., the water saturation increases proportionally to the square root of time. We also reproduce rate effects in the resistivity index for drainage and imbibition.     

Effect of network topology on two-phase imbibition relative permeability

In a previous study Arns et al. (2004, Transport Porous Media 55, 21–46) we considered the role of topology on drainage relative permeability curves computed using network models derived from a suite of tomographic images of Fontainebleau sandstone. The present study extends the analysis to more complex imbibition displacements where the non-wetting fluid can be disconnected by snap-off as a result of swelling of wetting films in the corners of pores and throats. In contrast to the findings for drainage displacements which showed that relative permeabilities are significantly affected by network topology, the present study shows that the effect of topology on imbibition relative permeabilities depends on the level of snap-off. For strongly wetting conditions where snap-off dominates the displacement the effect of network topology is significantly smaller than for weakly wet conditions where snap-off is suppressed. For contact angles sufficiently large to completely suppress snap-off, the effect of topology on imbibition relative permeabilities is similar to that for drainage displacements. The findings are valid for random networks and for networks displaying short-range pore–throat and longer range spatial correlations.     

Oil Displacement for One-Dimensional Three-Phase Flow with Phase Change in Fractured Media

In this article, the numerical simulations for one-dimensional three-phase flows in fractured porous media are implemented. The simulation results show that oil displacement in matrix is dominated by oil–water capillary pressure only under certain conditions. When conditions are changed to decrease the amount of water entering into the fractured media from the boundary of the flow field, water in fracture may be vaporized to superheated steam. In these cases, the appearance of superheated steam in fracture rather than in matrix will decrease the fracture pressure and generate the pressure difference between matrix and fracture, which results in oil flowing from matrix to fracture. Assuming that oil is wetting to steam, the matrix steam–oil capillary pressure will decrease the matrix oil-phase pressure as the matrix steam saturation increases. After the steam–oil capillary pressure finally exceeds the pressure difference due to the appearance of superheated steam in fracture, the oil displacement in matrix will stop. It is also shown that variations of the water relative permeability curve in matrix do not result in different mechanisms for oil displacement in matrix. The simulation results suggest that the amount of liquid water supply from the boundary of flow field fundamentally influence the mechanisms for oil displacement in matrix.     

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.     

Pore-throat size correlation from capillary pressure curves

Void spaces in porous media can be considered as three-dimensional networks consisting of bulges (pores) connected by constrictions (throats). Computer simulations of drainage-imbibition processes show that the critical end points of wetting-phase and nonwetting-phase saturation, in drainage and imbibition respectively, and the form of simulated relative permeability curves all were significantly different for uncorrelated and correlated pore-throat models. Since these models were identical except for the arrangement of throats in relation to pores, the degree of pore-throat size correlation appears to be an important property influencing flow and fluid displacement. Examples of uncorrelated and correlated pore-throat structures in rocks are presented and it is shown that this property, although difficult to quantify by direct observation, can be evaluated from capillary pressure curves.     

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.     

A Discussion of the Effect of Tortuosity on the Capillary Imbibition in Porous Media

In the past decades, there was considerable controversy over the Lucas–Washburn (LW) equation widely applied in capillary imbibition kinetics. Many experimental results showed that the time exponent of the LW equation is less than 0.5. Based on the tortuous capillary model and fractal geometry, the effect of tortuosity on the capillary imbibition in wetting porous media is discussed in this article. The average height growth of wetting liquid in porous media driven by capillary force following the [`(L)] _s(t) ~ t^1/2D_T{\overline L _{\rm {s}}(t)\sim t^{1/{2D_{\rm {T}}}}} law is obtained (here D _T is the fractal dimension for tortuosity, which represents the heterogeneity of flow in porous media). The LW law turns out to be the special case when the straight capillary tube (D _T = 1) is assumed. The predictions by the present model for the time exponent for capillary imbibition in porous media are compared with available experimental data, and the present model can reproduce approximately the global trend of variation of the time exponent with porosity changing.     

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.     

Co-history Matching: A Way Forward for Estimating Representative Saturation Functions

Core-scale experiments and analyses would often lead to estimation of saturation functions (relative permeability and capillary pressure). However, despite previous attempts on developing analytical and numerical methods, the estimated flow functions may not be representative of coreflood experiments when it comes to predicting similar experiments due to non-uniqueness issues of inverse problems. In this work, a novel approach was developed for estimation of relative permeability and capillary pressure simultaneously using the results of “multiple” corefloods together, which is called “co-history matching.” To examine this methodology, a synthetic (numerical) model was considered using core properties obtained from pore network model. The outcome was satisfactorily similar to original saturation functions. Also, two real coreflood experiments were performed where water at high and low rates were injected under reservoir conditions (live fluid systems) using a carbonate reservoir core. The results indicated that the profiles of oil recovery and differential pressure (dP) would be significantly affected by injection rate scenarios in non-water wet systems. The outcome of co-history matching could indicate that, one set of relative permeability and capillary pressure curves can reproduce the experimental data for all corefloods.     

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.     

Fluid transfer between tubes in interacting capillary bundle models

The interacting capillary bundle model proposed by Dong et al. [Dong, M., Dullien, F.A.L., Zhou, J.: Trans. Porous Media 31, 213–237 (1998); Dong, M., Dullien, F.A.L., Dai, L., Li, D.: Trans. Porous Media 59, 1–18 (2005); Dong, M., Dullien, F.A.L., Dai, L., Li, D.: Trans. Porous Media 63, 289–304 (2006)] has simulated correctly various aspects of immiscible displacement in porous media, such as oil production histories at different viscosity ratios, the effects of water injection rate and of the oil–water viscosity ratio on the shape of the displacement front and the independence of relative permeabilities of the viscosity ratio. In the interacting capillary bundle model pressure equilibrium was assumed at any distance x measured along the bundle. Interaction between the capillaries also results in transfer of fluids across the capillaries. In the first part of this paper the process of fluid transfer between two capillaries is analysed and an algebraic expression for this flow is derived. Consistency with the assumption of pressure equilibration requires that all transfer must take place at the positions of the oil/water menisci in the tubes without any pressure drop. It is shown that fluid transfer between the tubes has no effect on the predictions obtained with the model. In the second part of the paper the interacting tube bundle model is made more realistic by assuming fluid transfer between the tubes all along the single phase flow regions across a uniform resistance, resulting in pressure differences throughout the single phase regions between the fluids present in the different tubes. The results of numerical simulations obtained with this improved interacting capillary bundle model show only small differences in the positions of the displacement front as compared with the predictions of the idealized model.     

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.     

Effect of fluids properties on non-equilibrium capillarity effects: Dynamic pore-network modeling

We have developed a Dynamic Pore-network model for Simulating Two-phase flow in porous media (DYPOSIT). The model is applicable to both drainage and imbibition processes. Employing improved numerical and geometrical features in the model facilitate a physically-based pore-scale simulator. This computational tool is employed to perform several numerical experiments (primary and main drainage, main imbibition) to investigate the current capillarity theory. Traditional two-phase flow formulations state that the pressure difference between the two phase is equal to the capillary pressure, which is assumed to be a function of saturation only. Many theoretical and experimental studies have shown that this assumption is invalid and the pressure difference between the two fluids is not only equal to the capillary pressure but is also related to the variation of saturation with time in the domain; this is referred to as the non-equilibrium capillarity effect. To date, non-equilibrium capillarity effect has been investigated mainly under drainage. In this study, we analyze the non-equilibrium capillarity theory under drainage and imbibition as a function of saturation, viscosity ratio, and effective viscosity. Other aspects of the dynamics of two-phase flow such as trapping and saturation profile are also studied.     

Evaluation of Matrix-Fracture Transfer Functions for Counter-Current Capillary Imbibition

Different functions describing matrix-fracture transfer were tested for counter-current capillary imbibition interaction. The recovery curves obtained from capillary imbibition experiments were used to fit the transfer functions. The exponential coefficients yielding the best fit to the experimental data were obtained and correlated to the effective parameters such as viscosity, IFT, matrix length and diameter, matrix permeability and porosity, and wettability using multivariable regression analysis. In order to obtain the recovery curves, experiments were conducted on Berea sandstone and Indiana limestone samples. Cylindrical samples with different shape factors were obtained by cutting the plugs 1, 2.5, and 5 cm in diameter and 2.5, 5, and 10 cm in length. All sides were coated with epoxy except one end. More than fifty static imbibition experiments were carried out on vertically and horizontally situated samples where the imbibition took place upward and lateral directions, respectively. Brine–air, brine–kerosene, brine–mineral oil, and surfactant solution–mineral oil pairs were used as fluids. For many matrix shape factors (especially longer and small diameter ones), dividing the recovery curve into three parts were needed as the early, intermediate, and late times, which are typically distinguished by the time required for the imbibition front to reach the closed boundary at the end of the core. Correlations among the exponential coefficients and rock/fluid properties were developed. It was observed that different rock/fluid properties and transfer mechanisms (capillary imbibition and gravity drainage) govern the process for each part. Hence, the analyses done in this study were useful not only for developing explicit transfer functions but also identifying the physics of the counter-current imbibition recovery.     

多孔介质自发渗吸研究进展

自发渗吸是发生在多孔介质里的一种常见自然现象, 存在于众多工程应用和自然科学领域, 多孔介质 中自发渗吸的基本静力学和动力学问题已成为当前国际研究的热点课题之一. 本文综述了传统理论研究中的 Lucas{Washburn(LW) 模型, Terzaghi 模型, Handy 模型, Mattax 和Kyte 无因次时间标度模型, Aronofsky 归 一化采收率标度模型以及近十年最新研究进展, 分析了渗吸机理判别参数研究, 简述了数值模拟研究及渗吸率 影响机理的实验研究现状, 总结了基于分形理论研究多孔介质自发渗吸的最新进展, 并展望了多孔介质以及裂 缝性双重多孔介质中牛顿流体和非牛顿流体自发渗吸研究的方向和课题.      

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.