Thermodynamic modeling of calcium carbonate scale precipitation: aqueous Na + -Ca 2+ -Cl – -HCO 3 – -CO 3 2– -CO 2 system
To allow for accurate calculations of calcium carbonate scaling in highly saline produced waters, we present a comprehensive thermodynamic model based on the electrolyte nonrandom two-liquid (eNRTL) activity coefficient equation for the aqueous Na ⁺ -Ca ²⁺ -Cl – -HCO 3 – -CO 3 ²– -CO 2 system. The eNRTL binary interaction parameters for the H 2 O:(Na ⁺ -CO 3 2– ) pair, the H 2 O:(Na ⁺ -HCO 3 – ) pair, the (Na⁺-Cl–):(Na ⁺ -CO 3 2– ) pair, and the (Na ⁺ -Cl – ):(Na ⁺ -HCO 3 – ) pair are identified in this work via the regression of thermodynamic, calorimetric, and phase equilibria experimental data. The binary interaction parameters associated with the H 2 O:(Na ⁺ -Cl – ) pair, the CO 2 :(Na ⁺ -Cl – ) pair, the H 2 O:(Ca ²⁺ -Cl – ) pair, and the (Na ⁺ -Cl – ):(Ca ²⁺ -Cl – ) pair are retrieved from the literature. The remaining binary interaction parameters are retrieved from Aspen Plus or set to zero. In addition, the solubility product constants are identified for Na 2 CO 3 ·10H 2 O (s) , Na 2 CO 3 · 7H 2 O (s) , Na 2 CO 3 ·H 2 O (s) , Na 2 CO 3 ·NaHCO 3 2H 2 O (s) , Na 2 CO 3 ·3NaHCO 3 (s) , and CaCO 3(s) via regression of solubility data. Here, the model is capable of accurately calculating all phase equilibria and calorimetric properties at temperatures up to 473.15 K and salt concentrations up to saturation.