本研究采用分层积分方法对非均匀介质中的三维溶质运移问题进行了数值求解。忽略每层土壤和水力性质的垂直空间变化,三维溶质运移可以简化为准三维溶质迁移,它耦合了水平二维模拟和垂直一维计算。采用有限分析数值方法对导出的二维控制方程进行离散。为了保证相邻层间界面浓度和通量的连续性,采用二次函数近似计算了垂直一维浓度分布。通过对每一层的积分,可以为一列垂直单元生成一组方程组,并通过数值求解得到垂直溶质浓度分布。然后,通过求解垂直单元的所有列来获得溶质浓度场,从而通过迭代求解过程实现收敛。该模型通过已发表文献中的实例进行了验证,包括四个分析和实验案例的验证。仿真结果对比表明,该模型满足实验获得的溶质浓度时空分布。