[发明专利]一种土水特征曲线获取方法、系统、设备及存储介质

说明书

本发明公开了一种土水特征曲线获取方法、系统、设备及存储介质，建立累积百分数与颗粒粒径R的函数关系F(R)；确定任意颗粒粒径Ri所属的组分[Ri‑△R，Ri+△R]；根据S1和S2的结果，计算每个组分对应的土粒数量nsubgt;i/subgt;根据土粒数量nsubgt;i/subgt;，计算任意颗粒粒径Ri对应的任意孔隙半径ri；计算任意孔隙半径ri小于或等于某一孔隙半径r的累积百分数，根据累计百分数计算任意孔隙半径ri对应的土体的体积含水率θsubgt;vi/subgt;；根据任意孔隙半径ri计算土体的基质吸力；通过任意孔隙半径ri建立基质吸力与饱和度之间的关系，得到土水特征曲线。消除了随机粒组的划分进而引起的计算结果的差异，可以预测不同土体、不同颗粒级配、及不同压实度的连续性的土水特征曲线。

技术领域

本发明属于土水特征研究领域，涉及一种土水特征曲线获取方法、系统、设备及存储介质。

背景技术

在水文学中，土水特征曲线(SWCC)是定义非饱和土中土壤吸力和含水率之间关系最基本的概念之一，也是描述非饱和土水、温度和溶质运移建模的基本工具。通过室内试验获得完整的SWCC通常是非常耗时耗力的，因此从基本的非饱和土性质中建立SWCC数学模型是目前最可行的方案。现有的数学模型一般通过两种方式得到：经验模型(Van Genuchten和FredlundXing模型)和理论模型(Chen、MohammadiMeskini-Vishkaee模型)。经验模型往往是通过拟合许多离散的试验数据，形成一个连续的函数得到。这类模型通过许多拟合参数，往往可以很好的拟合试验数据。然而这些经验模型依赖于模型校准和试验数据的准确性，而且诸多拟合参数在物理上是没有意义的。因此许多学者提出了基于土壤物质属性，如矿物含量、类型，土壤的密度和孔隙形状大小的预测土水特征曲线。与经验模型相比，该类模型能普遍反映土壤性质对SWCC的影响，且不依赖大量的试验时间和数据。其中，一种考虑颗粒级配(GSD)对土水特征曲线预测的模型Arya and Paris模型，以下简称AP模型，已被广泛运用。AP模型是基于颗粒级配和孔隙条件对不同土样的土水特征曲线进行预测的概念和数学模型。虽然AP模型已经被许多学者所采用，但是在使用的过程中发现了许多问题，比如利用AP模型预测土体的基质吸力的大小与土样的质量密切相关，土样的质量越大，预测值则越大；AP模型的结果与颗粒级配粒组的划分密切相关，目前没有确定一个合理且统一的划分方式；AP模型没有实现土水特征曲线的连续化处理；从而会对得到土水特征曲线产生结果偏差，并且土水特征曲线精度较差。

发明内容

本发明的目的在于克服上述现有技术的缺点，提供一种土水特征曲线获取方法、系统、设备及存储介质，消除了随机粒组的划分进而引起的计算结果的差异，可以预测不同土体、不同颗粒级配、及不同压实度的连续性的土水特征曲线。

为达到上述目的，本发明采用以下技术方案予以实现：

一种土水特征曲线获取方法，包括以下过程：

S1、建立累积百分数与颗粒粒径R的函数关系F(R)；

S2、确定任意颗粒粒径Ri所属的组分[Ri-△R，Ri+△R]；

S3、根据S1和S2的结果，计算每个组分对应的土粒数量n_i；

S4、根据土粒数量n_i，计算任意颗粒粒径Ri对应的任意孔隙半径ri；

S5、计算任意孔隙半径ri小于或等于某一孔隙半径r的累积百分数；

S6、根据任意孔隙半径ri计算土体的基质吸力；

S7、根据累积百分数计算任意孔隙半径ri对应的土体的体积含水率θ_vi；

S8、根据基质吸力和体积含水率θ_vi，建立基质吸力与饱和度之间的关系，得到土水特征曲线。