Mercury Displacement Method

The mercury displacement method (also called mercury injection capillary pressure, MICP, or the Purcell method) is a laboratory technique for characterizing the pore structure of reservoir rock by injecting mercury under progressively increasing pressure into a cleaned and dried core plug, measuring the volume of mercury that enters the rock at each pressure step, and using the capillary pressure equation to convert the pressure-volume data into a pore throat size distribution that describes the rock's internal network of pore connections — the fundamental controls on permeability, residual saturations, and hydrocarbon column heights in the subsurface; mercury is used because it is a non-wetting fluid on virtually all geological surfaces at room temperature, meaning it can only enter pore spaces by capillary pressure being overcome rather than by wetting-driven spontaneous imbibition, and its high interfacial tension (approximately 480 mN/m) with air creates measurable pressures across pore throats in the micron and nanometer size range; the pressure required to force mercury into a pore throat of radius r is given by the Washburn equation (P = 2σ cos θ / r, where σ is the mercury-air interfacial tension and θ is the mercury contact angle with the rock), so measuring the injection pressure directly translates to the throat radius being invaded at that pressure, generating a complete pore throat size distribution from the injection curve; MICP provides information on total porosity (from the total mercury volume injected at maximum pressure), pore throat size distribution (the shape of the capillary pressure curve), pore throat sorting (steepness of the capillary pressure curve), entry pressure (the minimum pressure for mercury to enter the largest throats), and irreducible non-wetting phase saturation (the fraction of pore space mercury cannot access even at maximum pressure, corresponding to isolated pores or very small micro-porosity).

Key Takeaways

  • Pore throat size distribution from MICP is the definitive characterization of a rock's flow capacity — permeability is controlled not by the total pore volume (porosity) but by the size and connectivity of the pore throats that connect those pores; a rock with 20% porosity but very fine pore throats (tight shale or cemented sandstone) may have nanodarcy permeability, while a rock with 15% porosity but large, well-connected throats (clean aeolian sandstone) may have darcy-level permeability; MICP resolves this distinction by revealing the pore throat size distribution directly, allowing permeability prediction from the Winland R35 method (which correlates permeability to the pore throat radius at 35% mercury saturation) or more sophisticated petrophysical transforms; the pore throat distribution also classifies reservoir rocks into rock quality types that guide hydraulic unit definition, permeability-porosity transform calibration, and reservoir simulation cell property assignment.
  • MICP capillary pressure curves convert to reservoir fluid conditions for hydrocarbon column height prediction — mercury injection pressures measured in the laboratory must be converted to reservoir capillary pressures between actual reservoir fluids (oil-water or gas-water) using the ratio of interfacial tensions and contact angles; the converted capillary pressure curve relates water saturation to the height above the free water level in the reservoir, because capillary pressure at any height equals the pressure difference between the hydrocarbon and water phases (which increases with height above the free water contact proportional to density difference); the maximum hydrocarbon column that a seal formation can trap without leaking is the entry pressure of the seal converted to reservoir conditions — if the seal's entry pressure is exceeded by the buoyancy pressure of the hydrocarbon column, hydrocarbons will migrate through the seal and the trap will not retain a commercial column; MICP seal analysis is therefore a fundamental tool in trap integrity evaluation during prospect risking for exploration well decisions.
  • The drainage-imbibition asymmetry in MICP reveals irreversible pore-scale trapping — during mercury injection (drainage), mercury invades progressively smaller pore throats as pressure increases; during mercury withdrawal (imbibition), mercury does not simply retract along the same path; instead, mercury is trapped in large pore bodies connected to the outside only through small pore throats (snap-off mechanism), creating an irreversible residual mercury saturation analogous to residual oil saturation after waterflooding; the ratio of trapped mercury volume to total injected mercury volume characterizes the rock's propensity to trap non-wetting phase, which in reservoir terms corresponds to the proportion of the hydrocarbon pore volume that will be trapped as residual oil after waterflood displacement; rocks with large pore bodies connected by narrow throats (high aspect ratio) trap more mercury (and by analogy, more residual oil) than rocks where pore body and throat sizes are similar in scale.
  • MICP is the calibration dataset for NMR pore size distributions in log interpretation — nuclear magnetic resonance (NMR) logs measure the relaxation time distribution of formation fluids, which is sensitive to pore size (fluids in small pores relax faster than fluids in large pores); to convert the NMR T2 distribution (which reflects pore body size) to a permeability estimate, a surface relaxivity parameter must be calibrated by comparing the NMR T2 distribution to the pore throat size distribution from MICP on the same core samples; this NMR-MICP calibration, established for representative samples from each major rock type in the formation, allows NMR logs to provide continuous permeability profiles throughout the well rather than only at the discrete depths where core samples were analyzed; the quality of the NMR-derived permeability profile depends entirely on the quality of the NMR-MICP calibration, making MICP data one of the most important lab datasets for wells with NMR logging programs.
  • Sample preparation profoundly affects MICP results and must be standardized for inter-sample comparison — before mercury injection, core plugs must be cleaned to remove hydrocarbons and brines from the pore system (otherwise mercury cannot access the pore throats) and dried to remove all water (water in the pore space blocks mercury access and gives incorrect pore volume measurements); cleaning methods (Dean-Stark solvent extraction, flow-through cleaning, or high-temperature retort) vary in their aggressiveness and can alter clay minerals, micro-fractures, or diagenetic cements that affect pore throat size measurements; the drying temperature affects whether clay-bound water is removed along with free water (high-temperature drying above 105°C removes clay-bound water and may cause irreversible clay shrinkage that enlarges apparent pore throats); standardization of sample preparation procedures within a laboratory program is essential for ensuring that MICP results reflect actual formation rock properties rather than artifacts of the preparation process.

Fast Facts

The Purcell method of mercury injection capillary pressure measurement was published by W.R. Purcell in 1949 in the Transactions of the American Institute of Mining and Metallurgical Engineers — the same journal where many of the foundational reservoir engineering papers of the 1940s-1950s appeared. Purcell recognized that the capillary pressure curve contained information not just about fluid saturations but about the pore geometry itself, and derived the relationship between the integral under the capillary pressure curve and formation permeability that later became the basis for numerous permeability prediction methods. This 1949 insight remains the foundation of pore throat characterization by mercury injection over 75 years later.

What Is the Mercury Displacement Method?

The mercury displacement method is how geologists and petrophysicists map the internal plumbing of a reservoir rock without cutting it open. By forcing mercury — a fluid that will not spontaneously wet rock surfaces and must be pushed into pores by pressure — into a core sample at progressively higher pressures, and measuring how much mercury enters at each step, laboratories generate a detailed picture of the pore throat size distribution that controls everything from permeability to hydrocarbon column heights. The pressure at which mercury enters a pore throat directly reveals how large that throat is. Do it systematically from low pressure to high, and you've characterized the entire range of pore connections in the rock.

The mercury displacement method is also called mercury injection capillary pressure (MICP), the Purcell method, or high-pressure mercury porosimetry. Related terms include capillary pressure (the physical principle measured), pore throat size (the primary output of MICP), Washburn equation (the pressure-to-pore-size conversion), Winland R35 (the permeability prediction method using MICP), residual saturation (the trapping captured by imbibition MICP), seal capacity (the maximum hydrocarbon column height predicted from MICP), NMR log (the tool calibrated using MICP data), hydraulic flow unit (the rock classification derived from MICP), and permeability (the property predicted from pore throat data).

Why MICP Is Irreplaceable in Reservoir Rock Characterization

You cannot measure a pore throat with a ruler or a wireline log. Pore throats in reservoir sandstones are typically 1-100 microns in diameter — far smaller than any mechanical probe, and far too small for direct optical observation in routine thin section work. Yet pore throat size controls permeability, determines whether a trap can hold an economic hydrocarbon column, governs how much residual oil will be left behind after a waterflood, and provides the calibration that makes NMR logs quantitative rather than qualitative. MICP is the measurement that makes all of these connections possible, translating pressure into geometry at the scale where reservoir fluid behavior is actually determined. It's a 75-year-old technique that remains, fundamentally, one of the most information-rich measurements in reservoir rock physics.