An improved model for the calculation of CO2 solubility in aqueous solutions containing Na+, K+, Ca2+, Mg2+, Cl-, and SO42-

An improved model is presented for the calculation of the solubility of carbon dioxide in aqueous solutions containing Na+, K+, Cal(+), Mg2+, Cl-, and SO42- in a wide temperature-pressure-ionic strength range (from 273 to 533 K, from 0 to 2000 bar, and from 0 to 4.5 molality of salts) with experimental accuracy. The improvements over the previous model [Duan, Z. and Sun, R., 2003. An improved model calculating CO2 solubility in pure water and aqueous NaCl solutions from 273 to 533K and from 0 to 2000 bar. Chemical Geology, 193: 257-271] include: (1) By developing a non-iterative equation to replace the original equation of state in the calculation of CO2 fugacity coefficients, the new model is at least twenty times computationally faster and can be easily adapted to numerical reaction-flow simulator for such applications as CO, sequestration and (2) By fitting to the new solubility data, the new model improved the accuracy below 288 K from 6% to about 3% of uncertainty but still retains the high accuracy of the original model above 288 K. We comprehensively evaluate all experimental CO, solubility data. Compared with these data, this model not only reproduces all the reliable data used for the parameterization but also predicts the data that were not used in the parameterization. In order to facilitate the application to CO2 sequestration, we also predicted CO2 solubility in seawater at two-phase coexistence (vapor-liquid or liquid-liquid) and at three-phase coexistence (CO2 hydrate-liquid water-vapor CO2 [or liquid CO2]). The improved model is programmed and can be downloaded from the website http://www.geochem-model.org/programs.htm. (c) 2005 Elsevier B.V All rights reserved.