时间分数阶的非饱和渗流数值分析及其应用

Numerical Analysis of Time Fractional-Order Unsaturated Flow and Its Application

非饱和渗流过程的数值模拟对土质边坡稳定性分析、地下污染物迁移模拟等众多领域有着重要的意义.Richards方程由于其普遍适用性被广泛地应用,然而Richards方程所描述的渗流过程并未考虑在自然环境和实验中存在的反常扩散现象.针对这一问题,该文结合Caputo导数得到了具有更广泛渗流意义的时间分数阶Richards方程,采用有限差分法得到其离散格式并采用Picard法迭代求解,以及对分数阶参数和土水特征曲线进行了敏感性分析.最后,结合土柱入渗实验数据,比较了不同土水特征曲线下时间分数阶Richards方程得到的数值解.结果表明,VGM模型的时间分数阶Richards方程与实测数据具有更好的拟合效果,能够更好地描述地下水在非饱和土中的渗流过程.

Numerical simulation of the unsaturated flow process is of great significance to many fields such as soil slope stability analysis and migration simulation of underground pollutants.Generally,it is widely used due to the universal applicability of the Richards equation,but the seepage process described by the Richards equation does not involve the anomalous diffusion phenomenon in natural environment and experiments.To address this problem,the Caputo derivative was applied to obtain the time fractional-order Richards equation with broader seepage significance.Then the finite difference method was used to get the discretization scheme and the Picard method was chosen to solve it iteratively,and the sensitivity analysis of the fractional parameters and soil-water characteristic curves was carried out.Finally,combined with the experimental data of soil column infiltration,the numerical solutions obtained from the time fractional-order Richards equation under different soil-water characteristic curves were compared.The results show that,the time fractional-order Richards equation of the VGM model has better fitting effects for the measured data and can better describe the seepage process of groundwater in unsaturated soil.