非饱和土渗流问题在大坝渗流、降雨入渗、垃圾填埋场内部污染物扩散等领域有着广泛的应用。对非饱和渗流 Richards 方程的数值求解，一直以来多是采用有限单元法计算基质吸力水头为自变量的方程进行的。