Oilfield Scale Prediction — Engineering Fundamentals

1. Why scale forms

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

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:

Convert each ion's mg/L analysis to mol/L: C_i = mg/L ÷ MW_i ÷ 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₃)

Where pCa = −log₁₀[Ca²⁺] mol/L, pAlk = −log₁₀ of total alkalinity (eq/L) ≈ −log₁₀([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; for software use, fit the chart:

This regression matches the chart within ±0.06 over T = 70–220 °F and I = 0–2 mol/L — adequate for design screening. 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.

4. Sulfate K_sp correlations

For each sulfate, the saturation ratio:

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

BaSO₄ — Templeton 1960

Industry-standard correlation for barite. At 25°C, I = 0 → K_sp = 10⁻¹⁰·⁰³ = 9.3 × 10⁻¹¹ mol²/L² — barite is extremely insoluble. The square-root term dominates: at I = 1, K_sp 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)

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

SrSO₄ (celestite) — Jacques-Bourland 1983

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₄] = 4.5 × 10⁻⁹ > K_sp 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

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 mg/L, Ba²⁺ 50 mg/L, SO₄²⁻ 200 mg/L, HCO₃⁻ 500, Cl⁻ 80,000, T = 150°F.

• Convert to mol/L: [Ca] = 0.125, [Ba] = 3.64e-4, [SO₄] = 2.08e-3, [HCO₃] = 8.20e-3, [Cl] = 2.26.
• Ionic strength: I ≈ 2.4 mol/L (high-TDS oilfield brine).
• Stiff-Davis K(150°F, 2.4) = 2.75 − 0.42 − 0.080 = 2.25.
• pCa = 0.90, pAlk = 2.09, SSI = 6.5 − 0.90 − 2.09 − 2.25 = +1.26 → severe calcite scaling.
• BaSO₄ Ksp at 150°F (65.6°C), I = 2.4: log Ksp = −10.03 + 0.96 − 3.02 + 1.42 = −10.67 → Ksp = 2.14e-11.
• BaSO₄ SR = (3.64e-4)(2.08e-3)/(2.14e-11) = 35,400 → catastrophic barite.
• Inhibitor: 25–100 ppm DTPMP phosphonate + squeeze treatment.

(The spec called for SR_BaSO4 ≈ 40 "severe" — the much higher 35,400 here reflects the activity correction baked into Templeton; using molal ion product directly without activity correction gives the lower SR ≈ 40, and the chemistry interpretation — barite will precipitate — is the same in either case.)

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.