Issues With Gas And Water Relative Permeability In Low-Permeability Sandstones
Alan P. Byrnes
Kansas Geological Survey
Gas and water relative permeability can be effectively modeled in many porous media using the modified Corey (1954) equations:
krg = (1 – (Sw-Swc,g)/(1-Sgc-Swc,g))p (1-((Sw-Swc,g)/(1-Swc,g))q) (1)
krw = ((Sw-Swc)/(1-Swc))r (2)
where Sw= water saturation, Sgc = critical gas saturation (expressed as fraction gas saturation), Swc,g = critical water saturation for gas equation (expressed as fraction water saturation), Swc = critical water saturation, and p, q, and r are exponents reflecting pore size distribution and architecture. These variables change with such lithologic variables as grain size, sorting, and volume, type, and distribution, of diagenetic clay and secondary porosity. Beyond influences of lithologic variables within a homogeneous rock, lithologic anisotropy such as bedding planes parallel or perpendicular to flow will affect relative permeability. Corey and Rathjens (1956) showed that bedding planes perpendicular to flow tend to increase Sgc and that parallel bedding tends to decrease Sgc and cause inflections in both the krg and krw curves. In higher permeability rocks (i.e., k > 10 md) the measurement and operational definitions of Sgc and Swc can be clearly defined within a few saturation percent. However, as permeability decreases and rocks move more into the transition zone interval of the capillary pressure curve, the nature of both becomes more complex with implications for modeling using relative permeability.
Based on comparison of wireline log-measured water saturations and the lack of significant water production for low-permeability reservoirs, it is generally thought that many low-permeability sandstones are at or near “irreducible” water saturation (Swi) or critical water saturation (Swc). “Irreducible” saturation can be operationally defined as the saturation at which further increase in capillary pressure does not result in a “significant” decrease in water saturation. Critical water saturation can be operationally defined as the saturation at which water is immobile or water flow is negligible on the time scale of importance for the evaluation of flow properties. Critical water saturation is also often experimentally defined as the saturation at which the ratio of the nonwetting phase flow to water (wetting phase) flow is greater than 1000 (i.e. water flow represents less than 0.001 of total flow). Many western low-permeability reservoirs produce gas water-free or at very low water production rates indicating that water saturations are at or near Swc or that the water relative permeability (krw) is low. Critical gas saturation (Sgc) represents the saturation below which the gas phase is discontinuous and therefore does not flow. Experimentally this is defined as the saturation at which a threshold pressure achieves first detectable gas flow.
For higher permeability rocks Swi and Swc can be similar and for clean sandstones are often between 10-20%. Because of this low water saturation the effective gas permeability at or near “irreducible” water saturation (keg,Swi) is within 75-100% of the absolute single-phase permeability (ki). However, in low-permeability rocks, keg,Swi can be significantly less than ki because water occupies critical pore space even at or near Swc. Byrnes (2003) compiled low-permeability sandstone relative permeability data and showed that with relative permeability referenced to the in situ Klinkenberg gas permeability (kik) the gas relative permeability curve could be modeled using the following parameters:
Swc,g = 0.16 + 0.053*logkik (where if kik < 0.001 md then Swc,g = 0) (3)
Sgc = 0.15 - 0.05*logkik (4)
p = 1.7 (5)
q = 2 (6)
where Swc,g and Sgc are expressed in fractions and kik is expressed in md. In Figure 1 the bounding Corey-equation curves represent curves with Swc,g and Sgc values for rocks with permeability of 0.001 md and 1 md, which approximate the range of permeabilities of the samples for which relative permeability curves are shown. Equation 3 models Swc,g = 0.16 for kik = 1 md and approaches Swc,g = 0 at kik = 0.001 md. Conversely, Sgc = 0.15 for kik = 1 md and approaches Sgc = 0.30 at kik = 0.001 md. It is important to note that Swc,g calculated using Equation 3 is relevant only to the influence of water saturation on krg and does not correspond to critical water saturation relative to water flow (i.e., Swc). Swc,g defines the water saturation below which krg remains ~ 1. In contrast, Swc defines the water saturation below which krw is approximately zero. The Sgc and Swc terms are defined separately for each phase to allow calculation of relative permeability curves at saturations below those for which flow is unmeasurable or negligible.
Assuming the Corey equations presented here are approximately correct, Equation 4 can be interpreted to indicate that Sgc increases with decreasing absolute permeability. Conventionally, for higher permeability rocks, Sgc ranges from ~5-15% and is not highly dependent on permeability. All the data presented here can be interpreted to indicate that Sgc is greater than this in low-permeability sandstones, unless the samples have fractures or high-permeability laminae, which would provide flow paths with negligible saturation change. It is possible to postulate two aspects of pore architecture that explain the relationship of increasing Sgc with decreasing ki. In higher permeability rocks Sgc might be lower because the fraction of pores needed to be occupied by gas to achieve a percolation threshold (i.e. a single connective path through the pore network) is low since there are many pores with a significant fraction of the total pore volume that do not need to be gas saturated to achieve the percolation threshold. Conversely, in low-permeability rocks with thin sheet-like pores interconnecting a limited number of large pore bodies, that represent the majority of the pore space, if a connective path encountered just a few pore bodies, the filling of these would represent a greater fraction of the pore volume and consequently a higher Sgc. In addition, lower permeability rocks tend to have greater bedding complexity. As Corey and Rathjens (1956) reported, this increased heterogeneity could also be responsible for increasing Sgc with decreasing kik. Alternately, it is possible to model the observed krg curves with a similar constant low value of Sgc for all kik but with the exponent, p, changing with kik. This would imply that pore architecture and size distribution change with decreasing permeability.
To test the above models, threshold mercury injection capillary pressure measurements were performed coupled with electrical resistivity measurements to determine Sgc. This experimental methodology allowed the observation of the threshold capillary pressure and saturation at which a connective path was established as observed by an abrupt increase (several orders of magnitude) in electrical conductivity across the sample as mercury reached the end of the core providing a conductive path from one end of the core to the other. This electrical conductivity technique is highly sensitive to connectivity. These measurements confirm that Snwc (Sgc) increases with decreasing kik but Snwc also varied as a function of rock lithology, which was not handled in equation 4. This implies that both Sgc and p change with decreasing permeability and in response to lithology.
Equation 3 implies that krg at any given Sw increases with increasing absolute permeability, that is, the krg curves shift up to higher values of krg as kik increases or, alternately, krg curves shift to higher values of Sw as kik increases. The first explanation can be considered to represent the idea that in higher permeability rock at any given saturation gas occupies larger pores in a higher permeability sample compared to a lower-permeability sample and therefore the gas relative permeability is greater. The second explanation can be considered to represent the condition that in a high permeability sample higher water saturation is possible for the same krg because the water occupies smaller pores that are less influential or are inconsequential to flow. Both mechanisms are valid. Higher permeability samples can exhibit krg values at low-moderate Sw (~0-0.2) that are 75-100% of the dry gas permeability kik at Sw=0. This can result from water occupying pore space at Sw < Swc that is inconsequential to flow and therefore when gas occupies these same pores its relative additional contribution to gas flow is also negligible. Equation 3 predicts Swc>0.20 for higher permeability rocks (kik>10 md) with Swc approaching zero with decreasing kik down to 0.001 md and below. Equation 4 predicts Sgc<0.1 for higher permeability rocks and approaches Sgc = 0.30 for rocks of kik=0.001 md. Taken together these trends imply that as absolute permeability decreases, and pores become progressively more sheet-like, any water in the pore space interferes with gas flow and that greater gas saturation is needed to establish a connective path. These two conditions are consistent. Perhaps most important is that values for Sw at Sgc in lower-permeability rocks begin to approach saturations present in the reservoir. This would indicate that very low-permeability rocks have gas saturation that it is nearly immobile on the time frame relevant to commercial production. Commercially viable gas production from intervals with these properties requires the presence of a high permeability “carrier” channel (ie., fracture or high permeability bed) within close proximity to decrease flow distance and therefore time.
Jones and Owens (1981) first reported that water permeability is progressively less than Klinkenberg gas permeability with decreasing permeability for kik < 1 md. Ward and Morrow (1987) proposed the modification of the Corey water relative permeability equation 2, shown above, to calculate water relative permeability with values presented relative to kik by using the ratio of the kw/kik:
krw = ((Sw-Swc)/(1-Swc))4 (kw/kik) (7)
Because of the experimental difficulty of measuring krw in very low permeability rocks little work is published. In addition, because of the significant difference in the mobility (k/m; permeability/viscosity) of gas and water, gas/water flow ratios of 1000/1 are achieved at high water saturations (i.e., high Swc). At a limit Swc~ Swi, but definition of Swi is problematic in low-permeability rocks because rocks with kik< 0.1 md are still in the transition zone for gas column heights of hundreds of feet. Curves generated using higher, flow-measured, Swc values are useful for short-term flow prediction but can deviate from observed flow over longer time periods. As with krg, krw curves can be modeled either with a constant exponent, r, and variable Swc that increases with decreases kik, or with a constant Swc and changing value of r with decreasing kik. High pressure laboratory gas-water drainage experiments indicate that water is mobile to low saturations at high displacement pressures. This indicates that although increasing Swc with decreasing kik effectively models short-term flow, models involving r as a function of kik may more accurately represent water relative permeability behavior.
Byrnes, A.P., 2003, “Aspects of permeability, capillary pressure, and relative permeability properties and distribution in low-permeability rocks important to evaluation, damage, and stimulation”, Proceedings Rocky Mountain Association of Geologists – Petroleum Systems and Reservoirs of Southwest Wyoming Symposium, Denver, Colorado, September 19, 2003, 12 p.
Corey, A.T., 1954, “The interrelation between gas and oil relative permeabilities”, Producers Monthly, vol. 19, no. 1, November, p. 38-41.
Corey, A.T., and Rathjens, C.H., 1956, “effect of stratification on relative permeability”, J. of Petroleum Technology, Transactions, Am. Inst. Of Mining Engrs., Technical Note 393, December, p. 69-71.
Jones, F.O., and Owens, W. W., 1980, “A laboratory study of low-permeability gas sands”, Journal of Petroleum Technology, vol. 32, no. 9, p. 1631-1644.
Ward, J.S., and Morrow, N.R., 1987, “Capillary pressure and gas relative permeabilities of low permeability sandstone”, Soc. of Petroleum Engineers Formation Evaluation, Sept., p. 345-356.
Figure 1. Gas relative permeability curves for various western sandstones from Byrnes (2003). The bounding black curves represent curves constructed using equations 3 through 6 for rocks of 0.001 md and 1 md, representing the approximate upper and low limits of the samples for which curves are shown. Sgc = 0.3 and 0.15 for kik = 0.001 md and 1 md, respectively. Bounding red curves represent Eq. 3 with Sgc = 0.1 and p = 2.8 for kik=0.001md and p=2 for kik =1 md, respectively.