New Integrated Capillary Pressure and Initial Water Saturation Modeling

2020 ◽  
Author(s):  
Yizhu Liao
2005 ◽  
Vol 127 (3) ◽  
pp. 240-247 ◽  
Author(s):  
D. Brant Bennion ◽  
F. Brent Thomas

Very low in situ permeability gas reservoirs (Kgas<0.1mD) are very common and represent a major portion of the current exploitation market for unconventional gas production. Many of these reservoirs exist regionally in Canada and the United States and also on a worldwide basis. A considerable fraction of these formations appear to exist in a state of noncapillary equilibrium (abnormally low initial water saturation given the pore geometry and capillary pressure characteristics of the rock). These reservoirs have many unique challenges associated with the drilling and completion practices required in order to obtain economic production rates. Formation damage mechanisms affecting these very low permeability gas reservoirs, with a particular emphasis on relative permeability and capillary pressure effects (phase trapping) will be discussed in this article. Examples of reservoirs prone to these types of problems will be reviewed, and techniques which can be used to minimize the impact of formation damage on the productivity of tight gas reservoirs of this type will be presented.


SPE Journal ◽  
2012 ◽  
Vol 18 (02) ◽  
pp. 296-308 ◽  
Author(s):  
Y.. Zhou ◽  
J.O.. O. Helland ◽  
D.G.. G. Hatzignatiou

Summary It has been demonstrated experimentally that Leverett's J-function yields almost unique dimensionless drainage capillary pressure curves in relatively homogeneous rocks at strongly water-wet conditions, whereas for imbibition at mixed-wet conditions, it does not work satisfactorily because the permeability dependency on capillary pressure has been reported to be weak. The purpose of this study is to formulate a new dimensionless capillary pressure function for mixed-wet conditions on the basis of pore-scale modeling, which could overcome these restrictions. We simulate drainage, wettability alteration, and imbibition in 2D rock images by use of a semianalytical pore-scale model that represents the identified pore spaces as cross sections of straight capillary tubes. The fluid configurations occurring during drainage and imbibition in the highly irregular pore spaces are modeled at any capillary pressure and wetting condition by combining the free-energy minimization with an arc meniscus (AM)-determining procedure that identifies the intersections of two circles moving in opposite directions along the pore boundary. Circle rotation at pinned contact lines accounts for mixed-wet conditions. Capillary pressure curves for imbibition are simulated for different mixed-wet conditions in Bentheim sandstone samples, and the results are scaled by a newly proposed improved J-function that accounts for differences in formation wettability induced by different initial water saturations after primary drainage. At the end of primary drainage, oil-wet-pore wall segments are connected by many water-wet corners and constrictions that remain occupied by water. The novel dimensionless capillary pressure expression accounts for these conditions by introducing an effective contact angle that depends on the initial water saturation and is related to the wetting property measured at the core scale by means of a wettability index. The accuracy of the proposed J-function is tested on 36 imbibition capillary pressure curves for different mixed-wet conditions that are simulated with the semianalytical model in scanning-electron-microscope (SEM) images of Bentheim sandstone. The simulated imbibition capillary pressure curves and the reproduced curves, based on the proposed J-function, are in good agreement for the mixed-wet conditions considered in this study. The detailed behavior is explained by analyzing the fluid displacements occurring in the pore spaces. It is demonstrated that the proposed J-function could be applied to mixed-wet conditions to generate a family of curves describing different wetting states induced by assigning different wetting properties on the solid surfaces or by varying the initial water saturation after primary drainage. The variability of formation wettability and permeability could be described more accurately in reservoir-simulation models by means of the proposed J-function, and hence the opportunity arises for improved evaluation of core-sample laboratory experiments and reservoir performance.


2021 ◽  
Vol 11 (4) ◽  
pp. 1577-1595
Author(s):  
Rasoul Ranjbar-Karami ◽  
Parisa Tavoosi Iraj ◽  
Hamzeh Mehrabi

AbstractKnowledge of initial fluids saturation has great importance in hydrocarbon reservoir analysis and modelling. Distribution of initial water saturation (Swi) in 3D models dictates the original oil in place (STOIIP), which consequently influences reserve estimation and dynamic modelling. Calculation of initial water saturation in heterogeneous carbonate reservoirs always is a challenging task, because these reservoirs have complex depositional and diagenetic history with a complex pore network. This paper aims to model the initial water saturation in a pore facies framework, in a heterogeneous carbonate reservoir. Petrographic studies were accomplished to define depositional facies, diagenetic features and pore types. Accordingly, isolated pores are dominant in the upper parts, while the lower intervals contain more interconnected interparticle pore types. Generally, in the upper and middle parts of the reservoir, diagenetic alterations such as cementation and compaction decreased the primary reservoir potential. However, in the lower interval, which mainly includes high-energy shoal facies, high reservoir quality was formed by primary interparticle pores and secondary dissolution moulds and vugs. Using huge number of primary drainage mercury injection capillary pressure tests, we evaluate the ability of FZI, r35Winland, r35Pittman, FZI* and Lucia’s petrophysical classes in definition of rock types. Results show that recently introduced rock typing method is an efficient way to classify samples into petrophysical rock types with same pore characteristics. Moreover, as in this study MICP data were available from every one meter of reservoir interval, results show that using FZI* method much more representative sample can be selected for SCAL laboratory tests, in case of limitation in number of SCAL tests samples. Integration of petrographic analyses with routine (RCAL) and special (SCAL) core data resulted in recognition of four pore facies in the studied reservoir. Finally, in order to model initial water saturation, capillary pressure data were averaged in each pore facies which was defined by FZI* method and using a nonlinear curve fitting approach, fitting parameters (M and C) were extracted. Finally, relationship between fitting parameters and porosity in core samples was used to model initial water saturation in wells and between wells. As permeability prediction and reservoir rock typing are challenging tasks, findings of this study help to model initial water saturation using log-derived porosity.


1999 ◽  
Vol 2 (01) ◽  
pp. 25-36 ◽  
Author(s):  
A.B. Dixit ◽  
S.R. McDougall ◽  
K.S. Sorbie ◽  
J.S. Buckley

Summary The wettability of a crude oil/brine/rock system influences both the form of petrophysical parameters (e.g., Pc and krw/kro) and the structure and distribution of remaining oil after secondary recovery. This latter issue is of central importance for improved oil recovery since it represents the "target" oil for any IOR process. In the present study, we have developed a three-dimensional network model to derive capillary pressure curves from nonuniformly wetted (mixed and fractionally wet) systems. The model initially considers primary drainage and the aging process leading to wettability alterations. This is then followed by simulations of spontaneous water imbibition, forced water drive, spontaneous oil imbibition and forced oil drive—i.e., we consider a complete flooding sequence characteristic of wettability experiments. The model takes into account many pore level flow phenomena such as film flow along wetting phase clusters, trapping of wetting and nonwetting phases by snapoff and bypassing. We also consider realistic variations in advancing and receding contact angles. There is a discussion of the effects of additional parameters such as the fraction of oil-wet pores, mean coordination number and pore size distribution upon fractionally and mixed wet capillary pressure curves. Moreover, we calculate Amott oil and water indices using the simulated curves. Results indicate that oil recovery via water imbibition in weakly water-wet cores can often exceed that obtained from strongly water-wet samples. Such an effect has been observed experimentally in the past. The basic physics governing this enhancement in spontaneous water imbibition can be explained using the concept of a capillarity surface. Based on these theoretical calculations, we propose a general "regime based" theory of wettability classification and analysis. We classify a range of experimentally observed and apparently inconsistent waterflood recovery trends into various regimes, depending upon the structure of the underlying oil- and water-wet pore clusters and the distribution of contact angles. Using this approach, numerous published experimental Amott indices and waterflood data from a variety of core/crude oil/brine systems are analyzed. Introduction In crude oil/brine/rock (COBR) systems, pore level displacements of oil and brine and hence the corresponding petrophysical flow parameters (e.g., Pc and krw/kro) describing these displacements are governed by the pore geometry, topology and wettability of the system. A number of excellent review papers are available that describe experimental investigations of the effect of wettability on capillary pressure and oil-water relative permeability curves.1–5 In COBR systems, wettability alterations depend upon the mineralogical composition of the rock, pH and/or composition of the brine, crude oil composition, initial water saturation, reservoir temperature, etc.6–12 Therefore, in recent years, interest in restoring the wettability of reservoir core using crude oil and formation brine has greatly increased.3,4,13,14 In this approach, cleaned reservoir core is first saturated with brine and then oil flooded to initial water saturation using crude oil. The core containing crude oil and brine is then aged to alter its wettability state. Wettability measurements, such as Amott and USBM tests, and waterflood experiments are then typically conducted on the aged core. This entire process broadly mimics the actual flow sequences in the reservoir; consequently, the wettability alterations are more realistic than those achieved using chemical treatment methods. During the aging process, wettability may be altered to vastly different degrees depending upon many factors, including those mentioned above. In addition, aging time, thickness of existing water films and wetting film disjoining pressure isotherms also play important roles. Hence, the final wettability state of a re-conditioned core will generally be case specific.


2014 ◽  
Vol 2014 ◽  
pp. 1-12 ◽  
Author(s):  
Olugbenga Falode ◽  
Edo Manuel

An understanding of the mechanisms by which oil is displaced from porous media requires the knowledge of the role of wettability and capillary forces in the displacement process. The determination of representative capillary pressure (Pc) data and wettability index of a reservoir rock is needed for the prediction of the fluids distribution in the reservoir: the initial water saturation and the volume of reserves. This study shows how wettability alteration of an initially water-wet reservoir rock to oil-wet affects the properties that govern multiphase flow in porous media, that is, capillary pressure, relative permeability, and irreducible saturation. Initial water-wet reservoir core samples with porosities ranging from 23 to 33%, absolute air permeability of 50 to 233 md, and initial brine saturation of 63 to 87% were first tested as water-wet samples under air-brine system. This yielded irreducible wetting phase saturation of 19 to 21%. The samples were later tested after modifying their wettability to oil-wet using a surfactant obtained from glycerophtalic paint; and the results yielded irreducible wetting phase saturation of 25 to 34%. From the results of these experiments, changing the wettability of the samples to oil-wet improved the recovery of the wetting phase.


1964 ◽  
Vol 4 (03) ◽  
pp. 195-202 ◽  
Author(s):  
P.M. Blair

Abstract This paper presents numerical solutions of the equations describing the imbibition of water and the countercurrent flow of oil in porous rocks. The imbibition process is of practical importance in recovering oil from heterogeneous formations and has been studied principally by experimental means. Calculations were made for imbibition of water into both linear and radial systems. Imbibition in the linear systems was allowed to take place through one open, or permeable, face of the porous medium studied. In the radial system, water was imbibed inward from the outer radius. The effects on rate of imbibition of varying the capillary pressure and relative permeability curves, oil viscosity and the initial water saturation were computed. For each case studied, the rate of water imbibition and the saturation and pressure profiles were calculated as functions of time. The results of these calculations indicate that, for the porous medium studied, the time required to imbibe a fixed volume of water of a certain viscosity is approximately proportional to the square root of the viscosity of the reservoir oil whenever the oil viscosity is greater than the water viscosity. Results are also presented illustrating the effects on rate of imbibition of the other variables studied. Introduction The process of imbibition, or spontaneous flow of fluids in porous media under the influence of capillary pressure gradient s, occurs wherever there exist in permeable rock capillary pressure gradients which are not exactly balanced by opposing pressure gradients (such as those resulting from the influence of gravity). The importance of such capillary movement in the displacement of oil by water or gas was recognized in early investigations and described by Leverett, Lewis and True in 1942. Methods advanced by these authors for studying the process using dynamically scaled models were rendered more general and flexible by the research of later workers. The influence of capillary forces in laboratory water floods has also been discussed by several authors. While imbibition plays a very important role in the recovery of oil from normal reservoirs, Brownscombe and Dyes pointed out that imbibition might be the dominant displacement process in water flooding reservoirs characterized by drastic variations in permeability, such as in fractured- matrix reservoirs. In water-wet, fractured-matrix reservoirs, water will be imbibed from fractures into the matrix with a countercurrent expulsion of oil into the fractures. If the imbibition occurs at a sufficiently rapid rate, a very successful water flood can result; if the imbibition proceeds slowly the project might not be economically attractive. Scaled-model studies have demonstrated the vital importance of imbibition in secondary recovery in fractured reservoirs. It is therefore important in the evaluation of waterflooding prospects to develop a thorough understanding of the quantitative relationships of the factors which control the rapidity of capillary imbibition. The imbibition process serves reservoir engineers in still another important way by providing a technique for studying the wettability of reservoir core samples. Such experiments are usually conducted by observing the rate of expulsion of oil or water from core samples submerged in the appropriate fluid. Several papers have been published on the experimental techniques involved. Although Handy has recently published a method for calculating capillary pressures from experiments with gas-saturated cores, it has not yet been possible to deduce quantitative information regarding water-oil relative permeability and capillary pressure characteristics of the rock from the experimental results. Thus a technique is needed for studying the quantitative dependence of imbibition rate on oil and water viscosity, initial water saturation, relative permeability-saturation, and capillary pressure-saturation relations. The development of such information, including saturation and pressure profiles by laboratory experiments, would be very difficult. SPEJ P. 195ˆ


SPE Journal ◽  
2016 ◽  
Vol 21 (01) ◽  
pp. 152-169 ◽  
Author(s):  
Y.. Zhou ◽  
J. O. Helland ◽  
D. G. Hatzignatiou

Summary In this study, we present a three-phase, mixed-wet capillary bundle model with cross sections obtained from a segmented 2D rock image, and apply it to simulate gas-invasion processes directly on images of Bentheim sandstone after two-phase saturation histories consisting of primary drainage, wettability alteration, and imbibition. We calculate three-phase capillary pressure curves, corresponding fluid configurations, and saturation paths for the gas-invasion processes and study the effects of mixed wettability and saturation history by varying the initial water saturation after primary drainage and simulating gas invasion from different water saturations after imbibition. In this model, geometrically allowed gas/oil, oil/water, and gas/water interfaces are determined in the pore cross sections by moving two circles in opposite directions along the pore/solid boundary for each of the three fluid pairs separately. These circles form the contact angle with the pore walls at their front arcs. For each fluid pair, circle intersections determine the geometrically allowed interfaces. The physically valid three-phase fluid configurations are determined by combining these interfaces systematically in all permissible ways, and then the three-phase capillary entry pressures for each valid interface combination are calculated consistently on the basis of free-energy minimization. The valid configuration change is given by the displacement with the most favorable (the smallest) gas/oil capillary entry pressure. The simulation results show that three-phase oil/water and gas/oil capillary pressure curves are functions of two saturations at mixed wettability conditions. We also find that oil layers exist in a larger gas/oil capillary pressure range for mixed-wet conditions than for water-wet conditions, even though a nonspreading oil is considered. Simulation results obtained in sandstone rock sample images show that gas-invasion paths may cross each other at mixed-wet conditions. This is possible because the pores have different and highly complex, irregular shapes, in which simultaneous bulk-gas and oil-layer invasion into water-filled pores occur frequently. The initial water saturation at the end of primary drainage has a significant effect on the gas-invasion processes after imbibition. Small initial water saturations yield more-oil-wet behavior, whereas large initial water saturations show more-water-wet behavior. However, in both cases, the three-phase capillary pressure curves must be described by a function of two saturations. For mixed-wet conditions, in which some pores are water-wet and other pores are oil-wet, the gas/oil capillary pressure curves can be grouped into two curve bundles that represent the two wetting states. Finally, the results obtained in this work demonstrate that it is important to describe the pore geometry accurately when computing the three-phase capillary pressure and related saturation paths in mixed-wet rock.


2019 ◽  
Vol 89 ◽  
pp. 02002
Author(s):  
Stefano Pruno ◽  
Hans Erik Rodvelt ◽  
Ola Skjæveland

Capillary pressure and resistivity index spontaneous imbibition experiments by the porous plate method, in a core holder at elevated temperature and net reservoir stress, are both difficult and time consuming special core analysis measurements to perform. In this type of experiment, low capillary forces act against a low permeable porous plate and only one face of the cylindrical core sample is in capillary contact with the fluid saturated porous plate. In this paper, core samples having different lithology, petrophysical properties and wettability are analysed by counter current spontaneous imbibition, starting at initial water saturation (Swi), at net confining pressure and elevated temperature. Synthetic brine is used as the wetting phase and Isopar L mineral oil as the non-wetting phase. This methodology is applied to investigate and evaluate how to obtain more reliable, more efficient and faster saturation imbibition data combined with electrical measurements, during spontaneous imbibition measurements. Resistivity index (RI), saturation exponent (n) (by single saturation equilibrium point using Archie’s second law RI= Sw-n) and wettability information using representative fluids, confining pressure and temperature are also obtained by applying this specific counter current imbibition technique and improved procedural approach.


2016 ◽  
Vol 8 (1) ◽  
Author(s):  
István Nemes

AbstractThe main focus of the paper is to introduce a new approach at studying and modelling the relationship of initial water saturation profile and capillarity in water-wet hydrocarbon reservoirs, and describe the available measurement methods and possible applications. As a side track it aims to highlight a set of derivable parameters of mercury capillary curves using the Thomeer-method. Since the widely used mercury capillary pressure curves themselves can lead to over-, or underestimations regarding in-place and technical volumes and misinterpreted reservoir behaviour, the need for a proper capillary curve is reasonable. Combining the results of mercury and centrifuge capillary curves could yield a capillary curve preserving the strengths of both methods, while overcoming their weaknesses. Mercury injection capillary curves were normalized by using the irreducible water saturations derived from centrifuge capillary pressure measurements of the same core plug, and this new, combined capillary curve was applied for engineering calculations in order to make comparisons with other approaches. The most significant benefit of this approach is, that all of the measured data needed for a valid drainage capillary pressure curve represents the very same sample piece.


SPE Journal ◽  
2013 ◽  
Vol 19 (01) ◽  
pp. 88-100 ◽  
Author(s):  
Y.. Zhou ◽  
J.O.. O. Helland ◽  
D.G.. G. Hatzignatiou

Summary We simulate transient behavior of viscous- and capillary-dominated water invasion at mixed-wet conditions directly in scanning-electron-microscope (SEM) images of Bentheim sandstone by treating the pore spaces as cross sections of straight tubes. Initial conditions are established by drainage and wettability alteration. Constant rate or differential pressure is imposed along the tube bundle. The phase pressures vary with positions along the tube length but remain unique in each cross section, consistent with 1D core-scale models. This leads to a nonlinear system of equations that are solved for the interface positions as a function of time. The cross-sectional fluid configurations are computed accurately at any capillary pressure and wetting condition by a semianalytical model that is based on free-energy minimization. The fluid conductances are estimated by newly derived explicit expressions that are shown to be in agreement with numerical computations performed directly on the cross-sectional fluid configurations. An SEM image of Bentheim sandstone is taken as input to the developed model for simulating the evolution of saturation profiles during waterfloods for different flow rates and several mixed-wet conditions, which are established with various initial water saturations and contact angles. It is demonstrated that the simulated saturation profiles depend strongly on initial water saturation at mixed-wet conditions. The saturation profiles exhibit increasingly gradual behavior in time as the contact angle, defined on the oil-wet solid surfaces, increases or the initial water saturation decreases. Front menisci associated with positive capillary pressures promote oil displacement by water, whereas for large and negative capillary pressures at small flow rates, oil displaces water because the associated front menisci retract. This results in the development of pronounced gradual saturation fronts at mixed-wet conditions. The waterfloods simulated at conditions established with a large initial water saturation and small contact angle on the oil-wet solid surfaces exhibit sharp Buckley-Leverett saturation profiles for high flow rates because the capillary pressure is small and less important. The shape of the saturation profiles is interpreted on the basis of the simulated capillary pressure curves and the corresponding fluid configurations occurring in the rock image.


Sign in / Sign up

Export Citation Format

Share Document