A Thermodynamic Model for the Solubility Prediction of Barite, Calcite, Gypsum, and Anhydrite, and the Association Constant Estimation of CaSO4(0) Ion Pair up to 250 °C and 22000 psi