一种基于VG模型预测非饱和土相对渗透系数的方法技术

本发明专利技术公开了一种本发明专利技术基于VG模型预测非饱和土相对渗透系数的方法是以土样的土水特征曲线为参考状态，结合土水特征曲线分形模型、TK相对渗透系数模型与VG模型，利用分形理论，推导出一种非饱和相对渗透系数模型，然后用VG模型拟合土水特征曲线，求得进气值，进而预测土样的非饱和相对渗透系数，本预测方法克服了对渗透系数的预测以经验法为主、缺少理论支撑、模型较复杂、计算繁琐的缺点，对于非饱和土的渗流理论、流固耦合研究及工程应用具有极其重要的价值。

【技术实现步骤摘要】
一种基于VG模型预测非饱和土相对渗透系数的方法
本专利技术属于土木、岩土工程中非饱和土相对渗透系数研究
，具体涉及一种基于VG模型预测非饱和土相对渗透系数的方法。
技术介绍
非饱和土的复杂性和多变性决定了其渗透特性明显不同于饱和土，无法根据土壤的基本性质从理论上分析得出，试验难度也较大。对于土体的非饱和相对渗透系数是研究饱和/非饱和土渗流的关键参数，与饱和土相比，对非饱和土渗透系数的实验测定要困难得多，尤其是在低饱和度情况下，土中的水极难排出，因此通过实验室试验直接测试非饱和土渗透系数将耗费大量时间。而确定非饱和土的渗透系数很重要，因为非饱和土的渗透系数不仅在分析降雨对土坡稳定性影响时需要，而且在其他重要工程中，如地下水对垃圾填埋场和核废料填埋场的影响、地下污水的迁移、防洪堤渗漏等工程中都需要准确确定非饱和土的渗透系数。在实际工程(包括以上提到的工程)中所遇到的土体，不仅是处于非饱和或接近非饱和状态，而且还处于不同的应力状态下，因此，确定非饱和渗透系数是工程实践的需要，具有广阔的应用前景。
技术实现思路
本专利技术的目的就是针对上述技术的不足，提供一种既准确又简便的基于VG模型预测非饱和土相对渗透系数的方法。为实现上述目的，本专利技术所设计的基于VG模型预测非饱和土相对渗透系数的方法，包括如下步骤：1)结合TK相对渗透系数模型、土水特征曲线分形模型及VG模型提出一种预测非饱和土相对渗透系数模型，预测非饱和土相对渗透系数模型为：kr(ψ)为预测非饱和土相对渗透系数，a、m、n均为VG模型拟合参数，ψa为进气值，ψ为基质吸力；其中，VG模型拟合参数a的单位为kPa-1、其值为进气值ψa的倒数，VG模型拟合参数n与土的孔隙分布有关，VG模型拟合参数m与土体特征曲线的整体对称性有关；2)基于压力板试验测得土水特征曲线实测数据，得到测量基质吸力ψ1和测量饱和度Sr1；3)根据步骤2)中测量基质吸力ψ1和测量饱和度Sr1代入VG模型中拟合出土水特征曲线，得出测量拟合参数a0、m0、n0，根据拟合参数a与进气值ψa的关系，计算测量进气值ψa0；4)将设定的基质吸力ψ0、计算出测量进气值Ψa0和步骤3)中测量拟合参数a0,m0,n0代入公式Ⅳ中计算出非饱和土相对渗透系数，即为非饱和土相对渗透系数的预测值，其中ψ0>0。进一步地，所述步骤1)中，结合TK相对渗透系数模型、土水特征曲线分形模型及VG模型提出的预测非饱和土相对渗透系数模型具体过程如下：1a)TK相对渗透系数模型为其中kr(θ)为相对渗透系数、θ为体积含水率、ψ为基质吸力、θr为残余体积含水率、θs为饱和体积含水率；土水特征曲线分形模型为其中，Sr为饱和度、ψa为进气值、ψ为基质吸力、D为分维数；VG模型为其中，Se为有效饱和度，ψ为基质吸力，a、m、n均为VG模型拟合参数，其中VG模型拟合参数a的单位为kPa-1、其值为进气值ψa的倒数，VG模型拟合参数n与土的孔隙分布有关，VG模型拟合参数m与土体特征曲线的整体对称性有关；1b)将土水特征曲线分形模型公式Ⅱ两边同时求导，得由于体积含水率θ对应饱和度Sr或ψ，θr为残余体积含水率对应饱和度Srmin或ψd，饱和体积含水率θs对应饱和度为1或ψa，且ψd为最大基质吸力，并将公式Ⅴ代入公式Ⅰ中，得由于ψd>>ψa，所以忽略(ψa/ψd)5-D，并结合公式Ⅱ简化得到公式Ⅵ将Ⅲ式代入Ⅵ式，且有效饱和度Se与饱和度近Sr近似相等得到预测非饱和土相对渗透系数模型进一步地，所述步骤2)中，压力板试验所用仪器为体积压力板仪，气压力分别为5kPa、10kPa、25kPa、50kPa、100kPa、200kPa，绘出全部干燥曲线，完成干燥过程后试验继续沿浸湿过程进行，减小基质吸力，递减段气压分别为100kPa、50kPa、25kPa；试验结束后卸除气压，称量湿土样和烘干后土样重，根据量管起始和结束读数，计算出最后一组试样的含水量，然后反算相应于其他吸力值的质量含水量，然后推出体积含水量，最后绘制基质吸力与含水量关系曲线，即土水特征曲线。进一步地，所述步骤3)中，具体拟合过程如下：根据步骤2)中测量基质吸力ψ1和测量饱和度Sr1代入VG模型中，Sr1为VG模型的Sr、ψ1为VG模型的ψ，舍去含水量未改变的低基质吸力段数据，进而拟合出土水特征曲线，得出测量拟合参数a0、m0、n0，根据VG模型拟合参数a的值为进气值ψa的倒数，计算测量进气值ψa0，ψa0为测量拟合参数a0的倒数。与现有技术相比，本专利技术具有以下优点：本专利技术基于VG模型预测非饱和土相对渗透系数的方法是以土样的土水特征曲线(SWCC)为参考状态，结合土水特征曲线分形模型、TK相对渗透系数模型与VG模型，利用分形理论，推导出一种非饱和相对渗透系数模型，然后用VG模型拟合土水特征曲线，求得进气值，进而预测土样的非饱和相对渗透系数，本预测方法克服了对渗透系数的预测以经验法为主、缺少理论支撑、模型较复杂、计算繁琐的缺点，对于非饱和土的渗流理论、流固耦合研究及工程应用具有极其重要的价值。附图说明图1为实施例Yolo轻黏土土水特征曲线拟合图；图2为实施Yolo轻黏土相对渗透系数预测值与实测值对比图。具体实施方式下面结合附图和具体实施例对本专利技术作进一步的详细说明。一种基于VG模型预测非饱和土相对渗透系数的方法，具体方法如下：1)结合TK相对渗透系数模型(陶孔相对渗透系数模型)、土水特征曲线分形模型及VG(vangenuchten)模型提出一种预测非饱和土相对渗透系数模型，预测非饱和土相对渗透系数模型为：kr(ψ)为预测非饱和土相对渗透系数，a、m、n均为VG模型拟合参数(其中，VG模型拟合参数a的单位为kPa-1、其值为进气值ψa的倒数，VG模型拟合参数n与土的孔隙分布有关，VG模型拟合参数m与土体特征曲线的整体对称性有关)，ψa为进气值，ψ为基质吸力；结合TK相对渗透系数模型、土水特征曲线分形模型及VG模型提出的预测非饱和土相对渗透系数模型具体过程如下：1a)TK相对渗透系数模型为其中kr(θ)为相对渗透系数、θ为体积含水率、ψ为基质吸力、θr为残余体积含水率、θs为饱和体积含水率；土水特征曲线分形模型为其中，Sr为饱和度、ψa为进气值、ψ为基质吸力、D为分维数；VG模型为其中，Se为有效饱和度，ψ为基质吸力，a、m、n均为VG模型拟合参数，其中VG模型拟合参数a的单位为kPa-1、其值为进气值ψa的倒数，VG模型拟合参数n与土的孔隙分布有关，VG模型拟合参数m与土体特征曲线的整体对称性有关；1b)将土水特征曲线分形模型公式Ⅱ两边同时求导，得由于体积含水率θ对应饱和度Sr或ψ，θr为残余体积含水率对应饱和度Srmin或ψd，饱和体积含水率θs对应饱和度为1或ψa，且ψd为最大基质吸力，并将公式Ⅴ代入公式Ⅰ中，得由于ψd>>ψa，所以忽略(ψa/ψd)5-D，并结合公式Ⅱ简化得到公式Ⅵ将Ⅲ式代入Ⅵ式，且有效饱和度Se与饱和度近Sr近似相等得到预测非饱和土相对渗透系数模型2)基于压力板试验测得土水特征曲线实测数据，得到测量基质吸力ψ1和测量饱和度Sr1；压力板试验所用仪器为体积压力板仪，气压力分别为5kPa、10本文档来自技高网...

【技术保护点】
1.一种基于VG模型预测非饱和土相对渗透系数的方法，其特征在于：所述预测方法包括如下步骤：1)结合TK相对渗透系数模型、土水特征曲线分形模型及VG模型提出一种预测非饱和土相对渗透系数模型，预测非饱和土相对渗透系数模型为：

【技术特征摘要】
1.一种基于VG模型预测非饱和土相对渗透系数的方法，其特征在于：所述预测方法包括如下步骤：1)结合TK相对渗透系数模型、土水特征曲线分形模型及VG模型提出一种预测非饱和土相对渗透系数模型，预测非饱和土相对渗透系数模型为：kr(ψ)为预测非饱和土相对渗透系数，a、m、n均为VG模型拟合参数，ψa为进气值，ψ为基质吸力；其中，VG模型拟合参数a的单位为kPa-1、其值为进气值ψa的倒数，VG模型拟合参数n与土的孔隙分布有关，VG模型拟合参数m与土体特征曲线的整体对称性有关；2)基于压力板试验测得土水特征曲线实测数据，得到测量基质吸力ψ1和测量饱和度Sr1；3)根据步骤2)中测量基质吸力ψ1和测量饱和度Sr1代入VG模型中拟合出土水特征曲线，得出测量拟合参数a0、m0、n0，根据拟合参数a与进气值ψa的关系，计算测量进气值ψa0；4)将设定的基质吸力ψ0、计算出测量进气值Ψa0和步骤3)中测量拟合参数a0,m0,n0代入公式Ⅳ中计算出非饱和土相对渗透系数，即为非饱和土相对渗透系数的预测值，其中ψ0>0。2.根据权利要求1所述基于VG模型预测非饱和土相对渗透系数的方法，其特征在于：所述步骤1)中，结合TK相对渗透系数模型、土水特征曲线分形模型及VG模型提出的预测非饱和土相对渗透系数模型具体过程如下：1a)TK相对渗透系数模型为其中kr(θ)为相对渗透系数、θ为体积含水率、ψ为基质吸力、θr为残余体积含水率、θs为饱和体积含水率；土水特征曲线分形模型为其中，Sr为饱和度、ψa为进气值、ψ为基质吸力、D为分维数；VG模型为其中，Se为有效饱和度，ψ为基质吸力，a、m、n均为VG模型拟合参数，其中VG模型拟合参数a的单位为kPa-1、其值为进气...

【专利技术属性】
技术研发人员：陶高梁，罗晨晨，肖衡林，马强，万娟，朱志政，付佩，陈阳阳，彭弯，张琎炜，
申请(专利权)人：湖北工业大学，
类型：发明
国别省市：湖北,42