Archie Equation

Key Takeaways

• The formation factor F = a / phi^m is the fundamental rock property that links the rock's electrical tortuosity to its pore architecture: The formation factor is the ratio of the resistivity of a rock fully saturated with brine (Ro) to the resistivity of the brine itself (Rw): F = Ro / Rw. It depends only on the pore geometry (through the cementation exponent m) and the tortuosity of the pore paths (through the constant a). In simple intergranular sandstones, m typically ranges from 1.8 to 2.2: m is lower for loosely cemented or poorly sorted sands where the pore paths are more direct, and higher for tightly cemented or vugular carbonates where the pore paths are more tortuous. The Humble formula, developed from measurements on Gulf Coast sandstones, uses a = 0.62 and m = 2.15; the Archie cementation factor commonly used for carbonates is a = 1.0 and m = 2.0. Core-derived formation factors are measured by saturating clean core plugs with synthetic brine of known Rw and measuring the core resistivity to compute F = Ro / Rw at each sample; plotting F versus phi on a log-log scale gives a straight line whose slope is m and whose intercept at phi = 1 gives a. The accuracy of the Archie water saturation calculation is controlled primarily by the accuracy of m: a 10 percent error in m propagates to a 15 to 25 percent error in computed Sw depending on the porosity level, enough to misclassify a well as a producer or a duster in thin or marginal reservoirs.
• Formation water resistivity (Rw) must be known accurately for the Archie Equation to produce reliable water saturation estimates: Rw is measured from produced water samples, pressure bomb samples of downhole formation water, or estimated from the spontaneous potential (SP) log using the electrochemical potential theory that relates the SP deflection between a sand and shale to the ratio of Rw to Rmf (the filtrate resistivity). The SP-derived Rw is less accurate than the direct water sample measurement but is available from the standard open-hole log suite in any well. Rw varies spatially across a field due to differences in formation water salinity: in the Cardium Formation of central Alberta, Rw ranges from 0.020 to 0.055 ohm-metres depending on proximity to the paleo-recharge zone of the formation, with higher Rw (lower salinity) in the northwest updip areas and lower Rw (higher salinity) in the southeast down-dip areas. Using a single representative Rw for all wells in a field introduces a systematic spatial bias in the computed Sw: wells in higher-Rw areas will appear wetter than they are if the lower regional Rw is applied, potentially leading to unjustified poor-well assessments. Rw maps calibrated to produced water chemistry are therefore essential inputs to any multi-well formation evaluation study intended to guide development drilling decisions in the WCSB conventional play areas.
• The saturation exponent n is the least well-constrained Archie parameter and has the largest sensitivity to wettability and oil saturation level: The saturation exponent n relates the resistivity index (Ir = Rt/Ro) to water saturation: Ir = Sw^(-n), or equivalently Sw = Ir^(-1/n). For a water-wet rock with simple intergranular pores fully occupied by water film, n is typically 1.8 to 2.0. For an oil-wet or mixed-wet rock, where the hydrocarbon phase bridges across pore throats and creates longer electrical pathways through the water phase, n can range from 2.5 to 8.0 or higher. A rock with n = 4.0 (moderately oil-wet) at Rt/Ro = 10 would give Sw = 10^(-1/4) = 0.562, or 56 percent water saturation; using the water-wet assumption of n = 2.0 would give Sw = 10^(-0.5) = 0.316, or 32 percent water saturation: a massive 24 percentage point error in computed Sw that would cause the interpreter to dramatically underestimate the water saturation and overestimate recoverable oil. Reservoir wettability in the WCSB Cardium varies from strongly water-wet in clean quartzose zones to mixed-wet or oil-wet in organic-rich, bitumen-bearing facies, and the appropriate n value must be measured on representative core samples using RCAL (routine core analysis) and SCAL (special core analysis) resistivity measurements at reservoir conditions rather than defaulted to 2.0 for all facies without differentiation.
• The Archie Equation fails in shaly sands because clay minerals provide an additional electrical conduction pathway independent of Rw and Sw: Clay minerals (smectite, illite, kaolinite) have excess cations on their surface (characterised by the cation exchange capacity, CEC) that remain mobile even when the pore water is at low salinity, providing a surface conductance path that reduces the apparent formation resistivity below the value predicted by Archie. In a shaly sand with Vsh greater than 10 to 15 percent, the Archie Equation applied without shale correction will compute Sw higher than the true value (because it attributes the excess conductance from clay to additional pore water rather than clay surface conductance), leading to a pessimistic formation evaluation that classifies productive zones as wet. The Waxman-Smits model (1968) and the dual-water model (Clavier et al., 1977) provide physically motivated corrections that separate the clay conductance (expressed as QV, the CEC per unit pore volume in milliequivalents per mL) from the pore-brine conductance, recovering the correct Sw in the shaly intervals. In WCSB Viking shaly sands with clay volumes of 15 to 30 percent (common in the upper Viking), applying the simple Archie equation without shale correction systematically underestimates recoverable oil by 15 to 25 percent across the productive zone, a quantitatively significant error that has historically led to decisions to bypass productive intervals or to mis-rank development well priorities.
• Archie Equation application in carbonates requires special attention to the pore type because vugular and fracture porosity have different cementation exponents from intergranular porosity: In carbonate reservoirs, the porosity system commonly includes a mixture of intercrystalline (or intergranular) matrix porosity, moldic or vugular porosity, and fracture porosity. Each pore type has a different electrical tortuosity (and therefore a different effective m value). Fracture porosity, which is essentially a flat crack of very low tortuosity, has m approximately 1.0 (nearly linear resistivity-porosity relationship); vugular porosity, in which the vugs are poorly connected isolated cavities, has m greater than 2.5 to 3.0 (the isolated pores contribute to total porosity without contributing proportionally to connectivity). The dual-porosity or triple-porosity Archie models for carbonates use a volume-weighted m that depends on the fractions of each pore type at each depth. In Leduc reef carbonates at Redwater, Alberta, where vugular porosity constitutes 25 to 40 percent of total porosity in some dolomite facies, the effective m from core formation factor measurements is 2.6 to 3.2 rather than the standard 2.0 assumed in simple Archie, and using m = 2.0 would overestimate the formation factor and underestimate the water saturation by 10 to 20 percentage points, causing an overly optimistic evaluation of water saturation that could lead to premature well abandonment when high water cuts are encountered in production.