Capillary Pressure and Wettability Behavior of CO2 Sequestration in Coal at Elevated Pressures

SPE Journal ◽  
2008 ◽  
Vol 13 (04) ◽  
pp. 455-464 ◽  
Author(s):  
Willem-Jan Plug ◽  
Saikat Mazumder ◽  
Johannes Bruining

Summary Enhanced coalbed-methane (ECBM) recovery combines recovery of methane (CH4) from coal seams with storage of carbon dioxide (CO2). The efficiency of ECBM recovery depends on the CO2 transfer rate between the macrocleats, via the microcleats to the coal matrix. Diffusive transport of CO2 in the small cleats is enhanced when the coal is CO2-wet. Indeed, for water-wet conditions, the small fracture system is filled with water and the rate of CO2 sorption and CH4 desorption is affected by slow diffusion of CO2. This work investigates the wetting behavior of coal using capillary pressures between CO2 and water, measured continuously as a function of water saturation at in-situ conditions. To facilitate the interpretation of the coal measurements, we also obtain capillary pressure curves for unconsolidated-sand samples. For medium- and high-rank coal, the primary drainage capillary pressure curves show a water-wet behavior. Secondary forced-imbibition experiments show that the medium-rank coal becomes CO2-wet as the CO2 pressure increases. High-rank coal is CO2-wet during primary imbibition. The imbibition behavior is in agreement with contact-angle measurements. Hence, we conclude that imbibition tests provide the practically relevant data to evaluate the wetting properties of coal. Introduction Geological sequestration (Orr 2004) of CO2 is one of the viable methods to stabilize the concentration of greenhouse gases in the atmosphere and to satisfy the Kyoto protocol. The main storage options are depleted oil and gas reservoirs (Shtepani 2006; Pawar et al. 2004), deep (saline) aquifers (Kumar et al. 2005; Pruess et al. 2003; Pruess 2004), and unmineable coalbeds (Reeves 2001). Laboratory studies and recent pilot field tests (Mavor et al. 2004; Pagnier et al. 2005) demonstrate that CO2 injection has the potential to enhance CH4 production from coal seams. This technology can be used to sequester large volumes of CO2, thereby reducing emissions of industrial CO2 as a greenhouse gas (Plug 2007). The efficiency of CO2 sequestration in coal seams strongly depends on the coal type, the pressure and temperature conditions of the reservoir (Siemons et al. 2006a, 2006b), and the interfacial interactions of the coal/gas/water system (Gutierrez-Rodriguez et al. 1984; Gutierrez-Rodriguez and Aplan 1984; Orumwense 2001; Keller 1987). It can be expected that in highly fractured coal systems the wetting behavior positively influences the efficiency of ECBM recovery. It is generally accepted that the coal structure consists of the macrocleat and fracture system (>50 nm) and the coal matrix (<50 nm). The macrofracture system is initially filled with water and provides the conduits where the mass flow is dominated by Darcy flow. The coal matrix can be subdivided in mesocleats (from 2 to 50 nm), microcleats (from 0.8 to 2 nm), and the micropores (<0.8 nm). The matrix system is relatively impermeable, and the mass transfer is dominated by diffusion. After a dewatering stage, CO2 is injected and flows through the larger cleats of the coal. Subsequently, CO2 is transported through the smaller cleats and is sorbed in the matrix blocks (Siemons et al. 2006a). Depending on the wettability of coal, we can distinguish the following gas exchange mechanisms:The coal is water-wet, and CO2 and CH4 diffuse in the water-filled cleats.The coal is CO2-wet or gas-wet, and countercurrent capillary diffusion can take place.The coal is gas-wet, and binary diffusion of CO2 and CH4 occurs. Capillary diffusion finds its origin in capillary pressure (Pc) effects, where Pc is defined as the pressure difference between the nonaqueous and aqueous phase. The storage rate for CO2 is much smaller if the microcleat system is water-wet. This is because of the small CO2 molecular-diffusion coefficient (DCO2 ˜ 2 x 10-9 m2/s). For CO2-wet conditions, a faster and more efficient sorption rate is expected and the molecular diffusion is much larger (i.e., DCO2 ˜ 1.7 x 10-7 m2/s at 100 bar) (Bird et al. 1960). Therefore, we assert that the wettability of coal is important for ECBM recovery applications. For this reason, we have undertaken an experimental study to investigate the wetting properties of two different coal types under reservoir conditions, measuring the capillary pressure between CO2 and water. The dissolution properties of CO2 in water (Wiebe and Gaddy 1940), the interfacial tension between water and CO2 (Chun and Wilkinson 1995), and the CO2 sorption (Siemons et al. 2003) play important roles in the interpretation of capillary pressure experiments. The CO2, will sorb on the coal and will cause a swelling-induced permeability decrease (Mazumder et al. 2006). The higher the pressure, the more CO2 can be sorbed and the more the coal swells (Reucroft and Sethuraman 1987). The largest amount of sorption-induced swelling in intact coal is approximately 4%. It is found that the swelling for ground coal is much higher than intact coal and has been reported to be in the order of 15-20%. The swelling causes a porosity reduction, thus the water saturation decreases. In the Background section, relevant literature about the wettability of coal and the capillary pressure is summarized. The Experimental Design section describes the experimental setup we have developed to measure the capillary pressure as a function of the CO2 pressure. Furthermore, we describe the sample preparation and experimental procedure. In the Results and Discussion section, the experimental results are presented and discussed. We end with Conclusions.

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.


2007 ◽  
Vol 10 (03) ◽  
pp. 260-269 ◽  
Author(s):  
Eric P. Robertson ◽  
Richard L. Christiansen

Summary Sorption-induced strain and permeability were measured as a function of pore pressure using subbituminous coal from the Powder River basin of Wyoming, USA, and high-volatile bituminous coal from the Uinta-Piceance basin of Utah, USA. We found that for these coal samples, cleat compressibility was not constant, but variable. Calculated variable cleat-compressibility constants were found to correlate well with previously published data for other coals. Sorption-induced matrix strain (shrinkage/swelling) was measured on unconstrained samples for different gases: carbon dioxide (CO2), methane (CH4), and nitrogen (N2). During permeability tests, sorption-induced matrix shrinkage was demonstrated clearly by higher-permeability values at lower pore pressures while holding overburden pressure constant; this effect was more pronounced when gases with higher adsorption isotherms such as CO2 were used. Measured permeability data were modeled using three different permeability models that take into account sorption-induced matrix strain. We found that when the measured strain data were applied, all three models matched the measured permeability results poorly. However, by applying an experimentally derived expression to the strain data that accounts for the constraining stress of overburden pressure, pore pressure, coal type, and gas type, two of the models were greatly improved. Introduction Coal seams have the capacity to adsorb large amounts of gases because of their typically large internal surface area (30 to 300 m2/g) (Berkowitz 1985). Some gases, such as CO2, have a higher affinity for the coal surfaces than others, such as N2. Knowledge of how the adsorption or desorption of gases affects coal permeability is important not only to operations involving the production of natural gas from coalbeds but also to the design and operation of projects to sequester greenhouse gases in coalbeds (RECOPOL Workshop 2005). As reservoir pressure is lowered, gas molecules are desorbed from the matrix and travel to the cleat (natural-fracture) system, where they are conveyed to producing wells. Fluid movement in coal is controlled by diffusion in the coal matrix and described by Darcy flow in the fracture (cleat) system. Because diffusion of gases through the matrix is a much slower process than Darcy flow through the fracture (cleat) system, coal seams are treated as fractured reservoirs with respect to fluid flow. However, coalbeds are more complex than other fractured reservoirs because of their ability to adsorb (or desorb) large quantities of gas. Adsorption of gases by the internal surfaces of coal causes the coal matrix to swell, and desorption of gases causes the coal matrix to shrink. The swelling or shrinkage of coal as gas is adsorbed or desorbed is referred to as sorption-induced strain. Sorption-induced strain of the coal matrix causes a change in the width of the cleats or fractures that must be accounted for when modeling permeability changes in the system. A number of permeability-change models (Gray 1987; Sawyer et al. 1990; Seidle and Huitt 1995; Palmer and Mansoori 1998; Pekot and Reeves 2003; Shi and Durucan 2003) for coal have been proposed that attempt to account for the effect of sorption-induced strain. Accurate measurement of sorption-induced strain becomes important when modeling the effect of gas sorption on coal permeability. For this work, laboratory measurements of sorption-induced strain were made for two different coals and three gases. Permeability measurements also were made using the same coals and gases under different pressure and stress regimes. The objective of this current work is to present these data and to model the laboratory-generated permeability data using a number of permeability-change models that have been described by other researchers. This work should be of value to those who model coalbed-methane fields with reservoir simulators because these results could be incorporated into those reservoir models to improve their accuracy.


Geofluids ◽  
2021 ◽  
Vol 2021 ◽  
pp. 1-15
Author(s):  
Yingfang Zhou ◽  
Dimitrios Georgios Hatzignatiou ◽  
Johan Olav Helland ◽  
Yulong Zhao ◽  
Jianchao Cai

In this work, we developed a semianalytical model to compute three-phase capillary pressure curves and associated fluid configurations for gas invasion in uniformly wet rock images. The fluid configurations and favorable capillary entry pressures are determined based on free energy minimization by combining all physically allowed three-phase arc menisci. The model was first validated against analytical solutions developed in a star-shaped pore space and subsequently employed on an SEM image of Bentheim sandstone. The simulated fluid configurations show similar oil-layer behavior as previously imaged three-phase fluid configurations. The simulated saturation path indicates that the oil-water capillary pressure can be described as a function of the water saturation only. The gas-oil capillary pressure can be represented as a function of gas saturation in the majority part of the three-phase region, while the three-phase displacements slightly reduce the accuracy of such representation. At small oil saturations, the gas-oil capillary pressure depends strongly on two-phase saturations.


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.


SPE Journal ◽  
2013 ◽  
Vol 18 (04) ◽  
pp. 634-645 ◽  
Author(s):  
Y.. Zhou ◽  
J.O.. O. Helland ◽  
E.. Jettestuen

Summary In reservoir multiphase-flow processes with high flow rates, both viscous and capillary forces determine the pore-scale fluid configurations, and significant dynamic effects could appear in the capillary pressure/saturation relation. We simulate dynamic and quasistatic capillary pressure curves for drainage and imbibition directly in scanning-electron-microscope (SEM) images of Bentheim sandstone at mixed-wet conditions by treating the identified pore spaces as tube cross sections. The phase pressures vary with length positions along the tube length but remain unique in each cross section, which leads to a nonlinear system of equations that are solved for interface positions as a function of time. The cross-sectional fluid configurations are computed accurately at any capillary pressure and wetting condition by combining free-energy minimization with a menisci-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. Dynamic capillary pressure is calculated with volume-averaged phase pressures, and dynamic capillary coefficients are obtained by computing the time derivative of saturation and the difference between the dynamic and static capillary pressure. Consistent with previously reported measurements, our results demonstrate that, for a given water saturation, simulated dynamic capillary pressure curves are at a higher capillary level than the static capillary pressure during drainage, but at a lower level during imbibition, regardless of the wetting state of the porous medium. The difference between dynamic and static capillary pressure becomes larger as the pressure step applied in the simulations is increased. The model predicts that the dynamic capillary coefficient is a function of saturation and is independent of the incremental pressure step, which is consistent with results reported in previous studies. The dynamic capillary coefficient increases with decreasing water saturation at water-wet conditions, whereas, for mixed- to oil-wet conditions, it increases with increasing water saturation. The imbibition simulations performed at mixed- to oil-wet conditions also indicate that the dynamic capillary coefficient increases with decreasing initial water saturation. The proposed modeling procedure provides insights into the extent of dynamic effects in capillary pressure curves for realistic mixed-wet pore spaces, which could contribute to the improved interpretation of core-scale experiments. The simulated capillary pressure curves obtained with the pore-scale model could also be applied in reservoir-simulation models to assess dynamic pore-scale effects on the Darcy scale.


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.


1971 ◽  
Vol 11 (01) ◽  
pp. 13-22 ◽  
Author(s):  
Ali A. Sinnokrot ◽  
H.J. Ramey ◽  
S.S. Marsden

Abstract A number of recent studies of drainage relative permeability ratio by dynamic displacement have permeability ratio by dynamic displacement have indicated temperature sensitivity. Poston et al. found that the irreducible water saturation appeared to increase significantly with temperature-level increase and speculated that capillary pressure saturation data would also change to show this effect. Although there have been capillary pressure-saturation studies which show important pressure-saturation studies which show important differences between laboratory and reservoir conditions (presumably higher temperatures), the effects have usually been attributed to adsorption and desorption of polar components from the liquid phases. There appears to be no systematic studies phases. There appears to be no systematic studies of the effect of temperature level upon capillary pressure. pressure. Equipment was constructed to permit measuring capillary pressures for simple systems at temperatures ranging from room temperature to about 350 deg. F. Drainage and imbibition capillary pressure curves were measured for three consolidated pressure curves were measured for three consolidated sandstones and one limestone sample, at either three or four temperature levels form 70 deg to 325 deg F. Fluid used were a filtered white oil and distilled water. Results for the sandstone samples were similar. The practical irreducible water saturation increased significantly as temperature was raised from 70 deg F to the maximum temperature studies - about 325 deg F. Surprisingly, the hysteresis between drainage and imbibition cycles decreased as temperature increased and was nearly absent at 300 deg F. Results indicated that the sandstone samples became markedly more water-wet as temperature level increased. Results for the limestone sample were quite different. All capillary pressure-saturation curves for the various isotherms were found to lie within the envelope of the room-temperature drainage and imbibition curves. The main objective of this study was to determine whether the supposition of Poston et al. was correct. Results are in agreement with the previous dynamic displacement work. Introduction In 1967, Poston et al reported displacement experiments on unconsolidated sands at elevated temperatures and found that the irreducible water saturation increased with temperature increase. The oil viscosity appeared to have had no real effect on their results. Although less conclusive, practical residual oil saturations (to a producing practical residual oil saturations (to a producing water-oil ratio of 100) appeared to decrease with temperature increase. The results also indicated important increases in both oil and water relative permeabilities as temperature increased. This led permeabilities as temperature increased. This led Poston et al. to suggest that temperature affects. Poston et al. to suggest that temperature affects. the sand wettability. Sessile-drop contact angle measurements indicated that the water-oil-glass contact angle decreased with temperature increase. The results of Poston et al. regarding an increase in irreducible water saturation with temperature increase and the nondependence of this finding on the viscosity ratio deserve more attention. it is a well established concept in the literature that increasing water wetness of sands is reflected in an increase in the irreducible water saturation and an increase in oil recovery efficiency. The effect on a capillary pressure-saturation curve would be to cause a shift toward increasing irreducible wetting-phase saturation. if this is the case, then the studies of McNiel and Moss and Willman et al. should give partial credit for the added oil recovery efficiency involved in hot water flooding to the effect of temperature level upon wettability. In view of the potential importance of hot fluid injection for improving oil recovery and the lack of an adequate description of the flow process and thermodynamics involved, it was decided to study the speculation of Poston et al. that capillary pressure-saturation curves should be temperature pressure-saturation curves should be temperature dependent. This study concerns the effect of temperature level upon capillary pressure-saturation relationships for consolidated porous media. SPEJ p. 13


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.


Sign in / Sign up

Export Citation Format

Share Document