Saturation index (SI) is an important parameter to describe the water-rock interaction state, while complex has a significant influence on the state. It is possible to improve the accuracy of SI model by correctly evaluating the content of complex. According to the basic theory of chemical thermodynamics, a FORTRAN code was designed to calculate SI value for modifying solution method of saturation index model. It is proved that this model has better stabilization and precision.