Shaly Sand Water Saturation Guide
A practical comparison of Indonesian, Archie, and Simandoux water-saturation calculations for resistivity logs, with input checks for shale volume, shale resistivity, porosity, and water resistivity.
By PetroCalcHub Editorial Team | Updated 2026-07-31
linked formula references
verified source links
Engineering context
Water saturation from resistivity is an interpretation, not a direct measurement. Archie works well as a clean-formation reference because it relates formation resistivity to water resistivity, porosity, and saturation exponents without assigning conductivity to the rock matrix. In shaly rock, conductive clay-bound water can lower measured resistivity and make an uncorrected Archie result look wetter than the pore system actually is.
Indonesian and Simandoux equations introduce shale conductivity in different empirical forms. Neither is automatically superior. A defensible workflow estimates the same depth-matched inputs, runs more than one plausible model, and treats the difference as information about model uncertainty rather than selecting whichever result produces the most hydrocarbon.
Practical workflow
- 1
Confirm whether the interval is sufficiently clean for Archie or needs an explicit shale-conductivity model.
- 2
Depth-match effective porosity, shale volume, true resistivity, water resistivity, and representative shale resistivity.
- 3
Run Indonesian and a comparison model with the same exponents and unit basis.
- 4
Reconcile model spread with core, capillary-pressure, production, and local log-interpretation evidence.
Build a clean-formation reference first
Archie saturation uses formation factor and a saturation exponent. The formation factor is commonly expressed with tortuosity factor a, porosity exponent m, and porosity. The saturation exponent n controls how resistivity changes as hydrocarbon replaces conductive water. Defaulting a, m, and n to familiar values is useful for screening but should not be confused with local calibration.
Running Archie first shows how strongly the measured Rt differs from a clean-rock interpretation. The comparison also catches unit and depth errors before shale terms are introduced. If the interval is demonstrably clean, adding a shale correction can create complexity without improving the physical model.
Understand what the shale inputs represent
Vsh is a volumetric interpretation that depends on the selected log method and end points. Gamma ray, neutron-density separation, spectral logs, and multimineral solvers can produce different shale volumes. Rsh should represent the conductive shale component relevant to the interval and may vary with clay type, salinity, compaction, and bed geometry.
Effective porosity generally excludes shale-bound water volume, while total porosity includes more of it. Mixing a total-porosity curve with a formula calibrated on effective porosity can bias saturation. The model input label should therefore be tied to the petrophysical workflow, not only to a numeric curve mnemonic.
Use model comparison as quality control
As Vsh approaches zero, a shaly-sand model should approach clean-formation behavior on a consistent parameter basis. Large disagreement at very low Vsh can signal an implementation, exponent, porosity, or resistivity-basis problem. At higher Vsh, disagreement is expected because the models mix water and shale conductivity differently.
Results below zero or above one are diagnostic flags. Clipping them silently hides a problem. Review invasion, shoulder-bed effects, thin beds, anisotropy, hydrocarbons with unusual conductivity, pyrite, conductive minerals, bad-hole response, and the temperature basis of Rw before forcing a physical range.
Calibrate to reservoir evidence
Core electrical measurements, Dean-Stark saturation, capillary pressure, formation-water samples, mineralogy, NMR, dielectric logs, pressure data, and production behavior can constrain the interpretation. Each dataset has scale and sampling limitations, so calibration should preserve uncertainty instead of reducing the interval to one unsupported saturation curve.
For volumetrics, carry low, base, and high saturation cases when model and parameter uncertainty is material. A transparent range tied to geological and measurement evidence is more useful than a highly precise number from an uncalibrated empirical equation.
Worked Indonesian-versus-Archie check
Use Rt = 20 ohm m, Vsh = 0.20, Rsh = 2.5 ohm m, effective porosity = 0.20, a = 1, m = 2, Rw = 0.05 ohm m, and n = 2.
- 1. Formation responsesqrt(1 / Rt) = sqrt(1 / 20) = 0.2236
The Indonesian expression starts from the square root of measured formation conductivity.
- 2. Shale contributionShale term = 0.20^(1 - 0.5 x 0.20) / sqrt(2.5)
This term represents the empirical shale-conductivity contribution.
- 3. Formation-water contributionWater term = sqrt(0.20^2 / (1 x 0.05))
Effective porosity and water resistivity define the clean-water contribution.
- 4. Indonesian saturationImplemented result: Sw = 0.2144, or 21.44 percent
Because n equals 2, the final exponent is one for this parameter set.
Result
The implemented Indonesian equation returns Sw = 0.2144. A clean Archie comparison returns Sw = sqrt(0.05 / (20 x 0.20^2)) = 0.25.
Interpretation
The Indonesian result is about 21.4 percent water saturation for the stated inputs. Using Archie alone on the same Rt, Rw, porosity, a, m, and n gives 25 percent, so the shale correction changes the result by about 3.6 saturation units. That spread should be investigated, not automatically accepted as improved accuracy.
Method selection guide
Match the calculation method to the physical question and the evidence available.
| Condition | Use | Why |
|---|---|---|
| Clean, clay-free water-wet formation | Archie as the baseline model | Archie attributes formation conductivity to conductive pore water and does not include a shale-conductivity term. |
| Dispersed or structural shale contribution | Indonesian and Simandoux comparison | Both introduce shale conductivity, but their empirical mixing assumptions respond differently to Vsh and Rsh. |
| Laminated, low-resistivity, or complex mineralogy | Laminated or multimineral interpretation calibrated to core | A single homogeneous shaly-sand equation may not represent parallel conductive laminations or multiple clay distributions. |
Before using the result
- Depth-match porosity, Vsh, Rt, Rw, and Rsh before calculating saturation.
- Enter porosity, shale volume, and saturation as fractions rather than percentages.
- Correct Rw and measured resistivity to compatible temperature conditions.
- Document whether porosity is total or effective and use the model's required basis.
- Select Rsh from representative wet shale rather than an isolated minimum-resistivity point.
- Compare model spread with core, capillary pressure, mineralogy, and production evidence.
Start with these calculators
These links keep the guide close to the working calculator flow.