建立了有自由面的三维渗流饱和-非饱和计算模型,该模型包含水头和饱和度两个变量,用二维、三维结果的对比证明了模型的合理性与准确性。通过分析数据执行时的相关性,算法结构的固有串行性及有限元计算的耗时分布,利用OpenMP并行编程接口,实现了渗流求解的并行化。对不同规模的问题进行了并行数值模拟并测定了加速比,通过分析可知加速比与处理器的个数呈线性分布,并且随着问题规模的增加而逐渐提高。