To support the development of desalination processes for high salinity produced waters, we developed a comprehensive thermodynamic model that reliably predicts the thermodynamic properties and phase behavior of produced water. This study extends a previously developed thermodynamic model for the aqueous hexary oceanic salt system by incorporating the Sr2+ and Ba2+ ions. Based on the electrolyte nonrandom two-liquid theory, the model requires two adjustable binary interaction parameters for each water-electrolyte and electrolyte-electrolyte pair that has a common ion. The binary interaction parameters for the electrolyte-electrolyte pairs involving Sr2+ and Ba2+ ions were identified using thermodynamic and salt solubility data. The model was validated for temperatures from 273.15 K to 473.15 K and electrolyte concentrations from infinite dilution to salt saturation. Extension of the model to include HCO3−, CO32−, and CO2 is currently in progress.