水-岩作用模拟是环境水文地球化学研究领域的一个重要方面。本文根据质量作用定律、质量守恒定律和化学热力学的基本原理，应用平衡常数法建立了地下水水质组份存在形式（即单一离子、络阴离子和络合物）及其浓度分布的数学模型。在此基础上，建立了旨在反映水-岩（或矿物）系列间反应状态的饱和指数SI的数学模型，采用Newton-Raphson迭代法进行求解。在计算中，分别考虑络合物的影响和不考虑络合物的影响。结果表明，基于考虑络合物的影响所求得的SI值比较符合实际。此外，还对饱和指数的灵敏度作了较深入的分析，应用所建的灵敏度函数对地下水质分析资料各变量作了误差传递的计算，得出在所有变量中pH值的误差对SI的“反应”最为灵敏。上述数学模型的具体求解是通过用Fortrain 77语言所编写的程序得以实现的。工程实例应用表明，本文所建的数学模型思路正确、计算结果也较为合理。