Capillary pressure

Development Geology Reference Manual
Series	Methods in Exploration
Part	Laboratory methods
Chapter	Capillary pressure
Author	Charles L. Vavra, John G. Kaldi, Robert M. Sneider
Link	Web page
PDF	PDF file (requires access)

Capillary pressure concepts can be used by geologists, petrophysicists, and petroleum engineers to evaluate the following:

• Reservoir rock quality
• Pay versus nonpay
• Expected fluid saturations
• Seal capacity (thickness of hydrocarbon column a seal can hold before it leaks)
• Depth of the reservoir fluid contacts
• Thickness of the transition zone
• An approximation of the recovery efficiency during primary or secondary recovery.

Evaluating capillary pressure of potential reservoir and seal rocks is important because capillarity controls the static distribution of fluids in the reservoir prior to production and remaining hydrocarbons after primary production.

Capillary pressure concepts

Capillary pressure results from interactions of forces acting within and between fluids and their bounding solids. These include both cohesive forces (surface and interfacial tension) and adhesive (liquid-solid) forces. When adhesive forces are greater than cohesive forces, the liquid is said to be wetting (Figure 1a). When cohesive forces exceed adhesive forces, the liquid is nonwetting (Figure 1b). The relative wettability of the fluids is described by the contact angle (θ), which is the angle between the solid and the fluid-fluid interface as measured through the denser fluid (Figure 1). (For information on the measurement of wettability, see the chapter on Wettabilityin Part 5.)

File:Charles-l-vavra-john-g-kaldi-robert-m-sneider capillary-pressure 1.png

Figure 1 Effects of interaction of adhesive and cohesive forces on wettability. (a) If adhesive forces are greater than the cohesive forces, the fluid spreads out on the surface and is termed wetting. (b) If cohesive forces exceed adhesive forces, the liquid beads up and is termed nonwetting. The measure of relative wettability is the contact angle (θ), which is measured through the denser phase.

File:Charles-l-vavra-john-g-kaldi-robert-m-sneider capillary-pressure 2.png

Figure 2 The wetting phase rises above the original or free surface in the capillary tube experiment until adhesive and gravitational forces balance. Capillary pressure (P_c) is the difference in pressure measured across the interface in the capillary (P_c = P_nw - P_w). This pressure results from the contrast in pressure gradients caused by the different densities of the nonwetting (ρ_nw) and wetting (ρ_w) phases (right).

If the end of a narrow capillary tube is placed in a wetting fluid, net adhesive forces draw the fluid into the tube (Figure 2). The wetting phase rises in the capillary above the original interface or free surface until adhesive and gravitational forces are balanced. Because the wetting and nonwetting fluids have different densities, they also have different pressure gradients (Figure 2). Capillary pressure (P_c) is defined as the difference in pressure across the meniscus in the capillary tube. Put another way, capillary pressure is the amount of extra pressure required to force the nonwetting phase to displace the wetting phase in the capillary. Capillary pressure can be calculated as follows:

$P_{\mathrm {c} }=(\rho _{\mathrm {w} }-\rho _{\mathrm {nw} })gh$, or

$P_{\mathrm {c} }={\frac {2\sigma \cos \theta }{r_{\mathrm {c} }}}$

where

• ρ_w = density of the wetting nonwetting fluid
• ρ_nw = density of the nonwetting fluid
• g = gravitational constant
• h = height above the free surface
• σ = interfacial tension
• θ = contact angle between the fluids and the capillary tube
• r_c = radius of the capillary

These equations show that capillary pressure increases with greater height above the free surface and with smaller capillary size.

The importance of capillary pressure in reservoir studies is that many reservoir rocks can be approximated by a bundle of capillaries, with formation water being the wetting phase and hydrocarbons the nonwetting phase. As hydrocarbons begin to migrate into a rock, displacing the pore water, the hydrocarbons first enter the pores with the largest pore throats (capillaries), leaving the wetting phase (water) in the pores with smaller throats or in small nooks and crannies (surface roughness). As the hydrocarbon column increases, the height above the surface where P_c = 0, called the free surface or free water level (FWL), becomes greater and the capillary pressure increases, allowing hydrocarbons to enter pores with smaller and smaller throats. This process continues ntil one of several things occurs: (1) generation or migration ends, (2) the trap reaches its spill point, or (3) capillary pressure is sufficient to force hydrocarbons into the seal (displacement pressure is exceeded), allowing the seal to leak.

Reservoir capillary pressure relationships can be evaluated by using either of the following:

1. Porous plate or centrifuge method, which uses the actual or simulated hydrocarbon-brine system of the reservoir to approximate the wetting properties
2. Mercury injection, which simulates the premigration wetting characteristics of the reservoir

The mercury injection method is recommended to evaluate the initial static distributions of reservoir fluids prior to production because these distributions are typically controlled by the premigration wetting characteristics. Mercury injection is also favored because it is simpler, cheaper, and less time consuming than porous plate or centrifuge methods. In addition, mercury injection can be conducted on cuttings or sidewall samples. Mercury injection capillary pressure data can also be converted to the reservoir fluid system. Although this involves a number of assumptions on wettability and the effects of cleaning and extracting the core are largely ignored, the cost and time benefits typically favor mercury injection.

Terminology and mercury injection systematics

Mercury injection capillary pressure data are acquired by injecting mercury into an evacuated, cleaned, and extracted core plug. Mercury injection pressure is increased in a stepwise manner, and the percentage of rock pore volume saturated by mercury at each step is recorded after allowing sufficient time for equilibrium to be reached. The pressure is then plotted against the mercury saturation (Figure 3), resulting in the injection curve (which is also called the drainage curve because the wetting phase is being drained from the sample).

The pressure at which mercury first enters the sample (after the mercury has filled any surface irregularities on the sample) is termed the displacement pressure (P_d). The percentage of pore volume saturated by mercury at the maximum pressure is the maximum saturation (S_max). The unsaturated pore volume at that pressure is the minimum unsaturated pore volume (u_min) (Figure 3). This is sometimes incorrectly referred to as irreducible saturation. This term is inappropriate for the air-mercury system because saturation depends on applied pressure and on the duration of the experiment (Wardlaw and Taylor, 1976).

After the maximum pressure is reached, the pressure is reduced in steps and air (the wetting phase) is allowed to imbibe into the sample. The amount of mercury expelled from the sample at each pressure is expressed as a percentage of total pore volume or bulk volume. Again, pressure is plotted against mercury saturation in the withdrawal curve (Figure 3). The volume of pore space still saturated with mercury after pressure is reduced to the minimum is called the residual mercury saturation (S_R).

File:Charles-l-vavra-john-g-kaldi-robert-m-sneider capillary-pressure 3.png

Figure 3 Mercury injection capillary pressure curve terminology.

Data from the mercury injection curve can be used to approximate the distribution of pore volume accessible by throats of given effective radii (in cm) using the following equation:

$r\mathrm {c} ={\frac {2\sigma \cos \theta }{P_{\mathrm {c} }}}$

where

• σ = interfacial tension of the air-mercury system (480 dynes/cm)
• θ = air/mercury/solid contact angle (140°)
• P_c = capillary pressure in dynes/cm² (1 psi = 69035 dynes/cm²)

Note that this equation calculates the radius of a cylindrical capillary tube; however, real pore throats have very complex geometries. Therefore, the calculated values represent the effective radii of the pore throats, which may not equal their actual dimensions. Samples dominated by throats of similar size (well sorted) have broad, flat plateaus (Figure 4). As sorting becomes poorer, the plateau steepens then disappears and the slope of the curve approaches 45° (unsorted).

Data from the mercury withdrawal (air imbibition) curve can provide information regarding the efficiency with which the nonwetting phase can be withdrawn from the pore system. Withdrawal efficiency (W_E) is defined as the ratio of the mercury saturation in the sample at minimum pressure after pressure is reduced (S_R) to the saturation at maximum pressure (S_max) (Figure 3) (Wardlaw and Taylor, 1976). Because few samples reach 100% mercury saturation at routinely available injection pressures, data are normalized by the following equation:

$W_{\mathbf {E} }={\frac {S_{\mathrm {max} }-S\mathrm {R} }{S_{\mathrm {max} }}}\times 100\%$

Withdrawal efficiency in the air-mercury system is controlled exclusively by pore geometry, which can be interpreted from the capillary pressure curves and by direct observation of the rock or an epoxy resin pore cast. Critical pore geometry factors affecting withdrawal efficiency include effective pore throat heterogeneity; the ratio of pore body to pore throat size (r_b/r_t); and the number of throats connecting each pore (coordination number). Table 1 shows the effect of these parameters on withdrawal efficiency.

Relating withdrawal efficiency in the air-mercury system to recovery efficiency in the hydrocarbon-water system is dependent on properties of the fluids as well as properties of the pore system. Fluid properties that affect recovery include viscosity, density, interfacial tension, wettability, contact angle hysteresis, and rate of displacement (Wardlaw and Cassan, 1978). Nevertheless, for a given range of fluid properties in a water wet reservoir, the same pore geometry factors that contribute to increased mercury withdrawal efficiency also increase hydrocarbon recovery efficiency.

Reservoir applications

Conversion from air-mercury to brine-hydrocarbon system

Before mercury injection data can be applied to reservoirs, the data must be converted from the air-mercury system to the brine-hydrocarbon system of the reservoir using the following relationship:

$P_{\mathbf {C_{b/hc}} }=P_{\mathbf {C_{a/m}} }\times {\frac {\sigma _{\mathrm {b/hc} }\cos \theta _{\mathrm {b/hc} }}{\sigma _{\mathrm {a/m} }\cos \theta _{\mathrm {a/m} }}}$

where

• P_Cb/hc = capillary pressure in the brine-hydrocarbon system of the reservoir
• P_Ca/m = capillary pressure in the air-mercury system
• σ_b/hc = interfacial tension of the brine-hydrocarbon system
• σ_a/m = interfacial tension of air-mercury system
• θ_b/hc = contact angle of the brine-hydrocarbon-solid system
• θ_a/m = contact angle of the air-mercury-solid system

Ideally, values for brine-hydrocarbon contact angle, surface tension, and fluid densities at reservoir temperature and pressure should be determined in the laboratory using actual reservoir fluids. However, these measurements are difficult and expensive, so approximations such as those given in Tables 2 and 3 are commonly used.

File:Charles-l-vavra-john-g-kaldi-robert-m-sneider capillary-pressure 4.png

Figure 4 Idealized mercury injection capillary pressure curve shapes. Note that all of the curves have identical displacement pressures and minimum unsaturated pore volumes, but that the saturation profiles would differ dramatically due to differences in pore throat size distributions.

File:Table charles-l-vavra-john-g-kaldi-robert-m-sneider capillary-pressure 1.png

Table 1 Pore geometry factors affecting recovery efficiency (RE)

Height above the free water level

Having converted the capillary pressure data to the hydrocarbon-brine system of the reservoir and knowing the hydrocarbon and brine densities from laboratory analysis, we can now calculate the amount of hydrocarbon column required (height above the free water level or level where P_c = 0) to attain a pressure of interest:

$h={\frac {P_{\mathbf {h} }}{0.433(\rho _{\mathrm {b} }-\rho _{\mathrm {hc} }}}$

where

• h = height above the free water level (in ft)
• P_c = hydrocarbon-brine capillary pressure (in psi)
• ρ_b = specific density of brine at ambient conditions (in gm/cm³)
• ρ_hc = specific density of hydrocarbons at ambient conditions (in gm/cm³)
• 0.433 = the gradient of pure water at ambient conditions

This information can be used to compare expected fluid saturations of different rock types at given levels in the reservoir (Figure 5). (For information on the application of capillary pressure data to evaluating fluid contacts, see the chapter on Fluid contacts in Part 6.)

Evaluating seal capacity

If the capillary pressure data from the suspected seal are available or can be estimated, the maximum hydrocarbon column each rock type could have before the seal begins leaking can be calculated by the following:

$h_{\mathrm {max} }={\frac {P_{\mathbf {dS} }-P_{\mathrm {dR} }}{0.433(\rho _{\mathrm {w} }-\rho _{\mathrm {hc} }}}$

where

• h_max = height of the hydrocarbon column (in feet)
• P_dS = brine-hydrocarbon displacement pressure of the seal (in psi)
• P_dR = brine-hydrocarbon displacement pressure of the reservoir rock (in psi)

Additional applications

Capillary pressure data can also be applied to help distinguish reservoir from nonreservoir and pay from nonpay (see Effective pay determination). Several workers have attempted to correlate capillary pressure data and brine or air permeabilities. Purcell related capillary pressures empirically to air permeability through the graphical integral of the curve of mercury saturation versus reciprocal capillary pressure squared. Swanson (1981) proposed a simple nomograph whose application improved estimation of brine permeability from capillary pressure measurements on sidewall cores and ditch cuttings.

Another type of mercury test involves injecting mercury to a saturation less than the maximum, withdrawing the mercury to some residual wetting phase saturation, and then reinjecting the mercury. This process, repeated several times to progressively higher maximum pressures, produces hysteresis loops. These loops, wherein mercury is partially withdrawn and then reinjected, can be used to investigate withdrawal efficiency at various initial saturations (Morrow, 1970; Melrose and Brandner, 1974; Wardlaw and Taylor, 1976; Wardlaw and Cassan, 1978; Wardlaw et al., 1988). Results suggest that the higher the initial saturation of the nonwetting phase, the greater the withdrawal efficiency. Porosimetry uses hysteresis loops to interpret pore body and pore throat size distributions and their patial arrangement (Dullien and Dhawan, 1974; Wardlaw et al., 1988). Pore sizes have also been evaluated with rate-controlled mercury injection (Yuan and Swanson, 1986).

It should be noted that although these more specialized procedures are quite informative, especially for approximating hydrocarbon recovery efficiencies, they are also relatively labor intensive and expensive when compared to routine mercury injection tests.

File:Table charles-l-vavra-john-g-kaldi-robert-m-sneider capillary-pressure 2.png

Table 2 Typical contact angle and surface tension values

File:Table charles-l-vavra-john-g-kaldi-robert-m-sneider capillary-pressure 3.png

Table 3 Typical fluid density ranges

File:Charles-l-vavra-john-g-kaldi-robert-m-sneider capillary-pressure 5.png

Figure 5 Effect of capillary pressure (left) on water saturation (right). At any given height above the free water level, water saturations vary widely among rock types (A-E) due to diffferences in capillarity. For example, at 50 ft above free water level, water saturations vary from 18% (rock type A) to 95% (rock type E). A well drilled into an interbedded sequence of these rock types would show multiple oil-water contacts and a highly irregular vertical saturation profile. Note also the wide transition zone in rock type B caused by poor sorting of the pore throats.

See also

• Core description
• Introduction to laboratory methods
• Relative permeability
• Paleontology
• Wettability
• Oil and condensate analysis
• Oilfield water analysis
• Permeability
• Thin section analysis
• Overview of routine core analysis
• Core-log transformations and porosity-permeability relationships
• Rock-water reaction
• Sem, xrd, cl, and xf methods

External links

find literature about Capillary pressure

• Original content in Datapages
• PDF file in Datapages
• Dictionary:Capillary pressure in Sheriff's Encyclopedic Dictionary

feet from west line