Produced Water Chemistry

Oilfield Scale Prediction — Engineering Fundamentals

Stiff-Davis for calcite, plus Templeton, Marshall-Slusher and Jacques-Bourland Ksp correlations for the three common sulfate scales.

Carbonate index

Stiff-Davis SSI

High-TDS calcite scaling index (1952).

Sulfate scales

BaSO₄ / CaSO₄ / SrSO₄

Ksp(T, I) saturation-ratio correlations.

Scaling threshold

SR > 1 / SSI > +0.5

Precipitation favored — inhibitor required.

Use this guide when you need to:

  • Predict calcite scaling with the Stiff-Davis index at high TDS.
  • Screen barite, gypsum and celestite from Ksp saturation ratios.
  • Test water-mixing incompatibility and pick an inhibitor dose.

1. Why scale forms

Scale precipitates when an ion product (e.g. [Ca²⁺][CO₃²⁻]) exceeds the temperature-and-pressure-specific solubility product (Ksp). Four mechanisms push field waters above Ksp:

  1. Temperature change. Calcite and barite have retrograde solubility — Ksp falls as T rises, so heating a saturated water causes precipitation.
  2. Pressure change. CO₂ partial pressure drops as fluid moves up a wellbore or across a choke; carbonate equilibrium shifts and bicarbonate flashes to CO₂ + carbonate.
  3. Water mixing. Formation water (often Ba-rich, SO₄-poor) mixes with seawater or injection water (SO₄-rich); the product [Ba][SO₄] spikes far above Ksp.
  4. Evaporation. Tank vents or heater treaters concentrate dissolved solids; ions move past Ksp by simple concentration.

2. Ionic strength

All scale calculations need ionic strength because activity coefficients depend on it:

I = ½ · Σ ( Ci · zi2 ), Ci in mol/L

Convert each ion's mg/L analysis to mol/L: Ci = mg/L ÷ MWi ÷ 1000. A typical Bakken brine at 80,000 mg/L Cl⁻, 32,000 Na⁺, 5,000 Ca²⁺, etc. has I ≈ 1.5–2.5 mol/L. Seawater is ~0.7 mol/L.

For oilfield brines (I > 0.1), the Langelier Saturation Index (LSI) is no longer valid because it assumes activity = concentration. Stiff-Davis was developed specifically for high-TDS water by re-fitting the K constant to ionic strength.

3. Stiff-Davis Stability Index (CaCO₃)

SSI = pH − pCa − pAlk − K(T, I)

Where pCa = −log10[Ca²⁺] mol/L, pAlk = −log10 of total alkalinity (eq/L) ≈ −log10([HCO₃⁻] + 2·[CO₃²⁻]).

The K constant captures the temperature- and ionic-strength-dependent equilibrium constants for the carbonate system at high TDS. K is tabulated in the original 1952 Stiff & Davis paper as a chart (Fig. 1). K rises steeply with ionic strength to a broad peak near I ≈ 1.8–2.0, then declines gently at very high TDS, and it falls as temperature rises. The ionic-strength dependence is strong and non-monotonic, so a simple linear-in-log(I) regression is inadequate — it flattens and inverts the I-term and over-predicts SSI by ~1 unit at high TDS. This calculator interpolates the digitized Fig. 1 chart directly:

K — T ↓ / I →00.51.02.03.0
0 °C (32 °F)2.303.053.503.883.68
50 °C (122 °F)1.882.552.883.052.85
90 °C (194 °F)1.201.651.902.021.85

For final design refer to the full Stiff-Davis chart or the Oddo-Tomson modified correlations that include HP/HT and CO₂ partial-pressure terms.

SSIInterpretation
> +1.0Severe scaling — inhibitor required
+0.5 to +1.0Moderate scaling — inhibitor recommended
0 to +0.5Mild scaling — monitor
−0.5 to 0Stable
< −0.5Corrosive (CO₂-aggressive) — corrosion inhibitor likely needed

4. Sulfate Ksp correlations

For each sulfate, the saturation ratio:

SR = [M2+] · [SO42−] / Ksp(T, I)

SR > 1 means precipitation is thermodynamically favored. The Ksp correlations capture both intrinsic solubility (T) and the activity-coefficient corrections (I) without needing separate Debye-Hückel calculations:

BaSO₄ — Templeton 1960

log Ksp = −10.03 + 0.0146·T°C − 1.95·√I + 0.59·I

Industry-standard correlation for barite. At 25°C, I = 0 → Ksp = 10⁻¹⁰·⁰³ = 9.3 × 10⁻¹¹ mol²/L² — barite is extremely insoluble. The square-root term dominates: at I = 1, Ksp rises by 10^(−1.95 + 0.59) ≈ 23×, but a typical mix-water [Ba][SO₄] product is still well above this.

CaSO₄ (gypsum) — Marshall-Slusher 1966 (simplified)

log Ksp = −4.61 + 0.005·T°C − 0.7·√I + 0.2·I

Calcium sulfate is much more soluble than barite (Ksp ≈ 2.4 × 10⁻⁵ at 25°C, I = 0). Gypsum forms below ~40°C; anhydrite (CaSO₄ anhydrous) forms above ~98°C with similar Ksp.

SrSO₄ (celestite) — Jacques-Bourland 1983

log Ksp = −6.50 + 0.005·T°C − 1.0·√I + 0.35·I

Often a "hidden" scale — Sr²⁺ is present in many formation waters at 100–500 mg/L and can precipitate alongside barite when SO₄ is introduced.

5. Water incompatibility

The classic mixing problem: a formation water rich in Ba²⁺ but poor in SO₄²⁻ meets a seawater-injection stream rich in SO₄²⁻ but poor in Ba²⁺. Each water alone is stable; the mix is dramatically supersaturated:

  • Formation: Ba²⁺ = 100 mg/L, SO₄²⁻ = 10 mg/L → [Ba][SO₄] = 7.6 × 10⁻⁸ > Ksp at I = 0 but stable at high I.
  • Seawater: Ba²⁺ = 0.02 mg/L, SO₄²⁻ = 2,700 mg/L → [Ba][SO₄] tiny, stable.
  • 50/50 mix: Ba²⁺ = 50 mg/L, SO₄²⁻ = 1,355 mg/L → [Ba][SO₄] = 5.1 × 10⁻⁶ → SR ≈ 10,000× → severe barite scaling.

Field practice: compatibility-test all mixing ratios from 10/90 to 90/10 before commingling any two waters. Re-test annually as formation water chemistry shifts with depletion.

6. Inhibitor selection

SeverityChemistryDoseMethod
Mild (SR 1–2)HEDP phosphonate5–10 ppmContinuous
Moderate (SR 2–10)DTPMP or polyacrylate10–25 ppmContinuous wellhead or downhole
Severe (SR > 10)DTPMP + polyacrylate blend25–100 ppmContinuous + periodic squeeze
Barite SR > 10Aminomethylene-phosphonate50–100 ppm + squeezeSqueeze every 6–12 months

Phosphonates work by site-blocking nucleation; polyacrylates work by adsorbing on growing crystals and distorting them. Both are dosed in ppm because they are catalytic, not stoichiometric.

7. Worked example — Permian formation water

From the calculator sample: Ca²⁺ 5,000, Na⁺ 32,000, Mg²⁺ 500, Ba²⁺ 50, Sr²⁺ 200, SO₄²⁻ 200, HCO₃⁻ 500, Cl⁻ 80,000 mg/L, pH 6.5, T = 150°F.

  • Convert to mol/L: [Ca] = 0.125, [Na] = 1.39, [Ba] = 3.64e-4, [SO₄] = 2.08e-3, [HCO₃] = 8.20e-3, [Cl] = 2.26.
  • Ionic strength (all ions, incl. Na/Mg/Sr): I ≈ 2.13 mol/L (high-TDS oilfield brine).
  • Stiff-Davis K(150°F = 65.6°C, I = 2.13) ≈ 2.64, read from the Fig. 1 chart.
  • pCa = 0.90, pAlk = 2.09, SSI = 6.5 − 0.90 − 2.09 − 2.64 = +0.87 → moderate calcite scaling.
  • BaSO₄ Ksp at 65.6°C, I = 2.13: log Ksp = −10.03 + 0.96 − 2.85 + 1.26 = −10.66 → Ksp = 2.18e-11.
  • BaSO₄ SR = (3.64e-4)(2.08e-3)/(2.18e-11) = ≈34,800 → catastrophic barite (≈30 lb BaSO₄ per 1,000 bbl).
  • Inhibitor: 25–100 ppm DTPMP phosphonate + squeeze treatment.

(The saturation ratio is enormous because barite is extraordinarily insoluble: even though the ionic-strength term raises Ksp roughly 20-fold at I ≈ 2, the ion product still exceeds Ksp by ~35,000×. The practical conclusion — barite will precipitate and an inhibitor squeeze is mandatory — is unambiguous.)

8. References

  • Stiff, H. A. & Davis, L. E. (1952). "A Method for Predicting the Tendency of Oil Field Waters to Deposit Calcium Carbonate." JPT, Sept 1952.
  • Templeton, C. C. (1960). "Solubility of barium sulfate in sodium chloride solutions from 25° to 95° C." J. Chem. Eng. Data 5(4), 514–516.
  • Marshall, W. L. & Slusher, R. (1966). "Thermodynamics of calcium sulfate dihydrate in aqueous sodium chloride solutions, 0–110°C." J. Phys. Chem. 70(12), 4015–4027.
  • Jacques, D. F. & Bourland, B. I. (1983). "A study of solubility of strontium sulfate." SPE Journal 23(2).
  • Oddo, J. E. & Tomson, M. B. (1994). "Why scale forms in the oil field and methods to predict it." SPE Production & Facilities 9(1), 47–54.
  • NACE TM0374 — Laboratory Screening Tests to Determine the Ability of Scale Inhibitors to Prevent Scale Deposition.
  • NACE SP0775 — Preparation, Installation, Analysis, and Interpretation of Corrosion Coupons in Oilfield Operations.
  • Kemmer, F. N. (1988). The NALCO Water Handbook, 2nd ed.

Frequently Asked Questions

Why use Stiff-Davis instead of the Langelier index?

The Langelier Saturation Index assumes activity equals concentration, which fails for oilfield brines where ionic strength exceeds 0.1 mol/L. Stiff-Davis was developed for high-TDS water by re-fitting the K constant to ionic strength, so it stays valid for the 1.5–2.5 mol/L brines common in the field.

What causes water-mixing scale incompatibility?

A formation water rich in Ba²⁺ but poor in SO₄²⁻ can be stable on its own, as can a SO₄-rich seawater injection stream. Mixed, the [Ba][SO₄] product spikes far above Ksp — a 50/50 mix can reach a saturation ratio near 10,000× and cause severe barite scaling. Compatibility-test all mixing ratios before commingling.

How is a scale inhibitor selected and dosed?

Dose scales with severity: ~5–10 ppm HEDP phosphonate for mild (SR 1–2), 10–25 ppm DTPMP or polyacrylate for moderate, and 25–100 ppm blends plus periodic squeeze for severe or high-SR barite. Phosphonates block nucleation and polyacrylates distort growing crystals — both are catalytic, so they are dosed in ppm rather than stoichiometrically.