CN103971002A - 一种非饱和土相对渗透系数的计算方法 - Google Patents
一种非饱和土相对渗透系数的计算方法 Download PDFInfo
- Publication number
- CN103971002A CN103971002A CN201410199064.9A CN201410199064A CN103971002A CN 103971002 A CN103971002 A CN 103971002A CN 201410199064 A CN201410199064 A CN 201410199064A CN 103971002 A CN103971002 A CN 103971002A
- Authority
- CN
- China
- Prior art keywords
- ref
- value
- under
- parameter
- swcc
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明公开一种非饱和土相对渗透系数的计算方法,该方法包括针对某一类型的土样,选定某初始孔隙比为参考状态,利用该参考状态下的土水特征曲线(SWCC)标定出van Genuchten SWCC方程参数a1,b1;再利用该类型土样在初始孔隙比为时的土水特征曲线,标定出参数ζ;基于标定出的参数a1,b1,ζ,可预测出该类型土样在任意初始孔隙比e0条件下的土水特征曲线;基于e0条件下的土水特征曲线,再次拟合出其相应的van Genuchten SWCC方程参数a2,b2,并代入本发明关系式进行计算,即可得出相应的非饱和土相对渗透系数值kr。采用本发明所述技术方案可有效克服在实验直接测定渗透系数上的难度,从而降低人力、物力及时间上的消耗,同时大幅度减小实际工程中的工作量,有效提高工作进度。
Description
技术领域
本发明涉及土木、水利及环境工程中非饱和土相对渗透系数的计算,其特别是一种可对不同初始孔隙比条件下的非饱和土相对渗透系数进行有效预测的计算方法。
背景技术
渗流问题是土力学的三大基本问题之一。随着理论研究的深入和实际工程的发展,非饱和土渗流问题在岩土工程及环境岩土工程中受到越来越多的重视,例如降雨引起的土质边坡饱和度及渗透性变化,并最终导致的滑坡灾害;地基基础或路基中压实土的固结及膨胀土的隆起;土质堤坝的防渗;垃圾填埋场、核废料处置库周边地下水中污染物的迁移等。
对于可压缩的饱和土,其渗透系数主要受到孔隙比的影响;对于不可压缩的非饱和土,其渗透系数是饱和度或吸力的函数。而实际工程中遇到的土体,往往既可压缩又处于非饱和状态,因此其渗透系数将同时受饱和度和孔隙比的影响,是它们共同的函数。岩土工程界较早地认识到了可压缩饱和土的渗透系数随孔隙比而变化的规律;不考虑变形的影响,研究非饱和土渗透系数随饱和度/吸力的变化规律,则起源于土壤学。而对于不同孔隙比条件下非饱和土渗透系数演化规律的研究,则相对较少,相关研究成果也很难直接应用于工程实践。
渗透系数是研究饱和/非饱和土渗流的关键参数,与饱和土相比,对非饱和土渗透系数的试验测定要困难得多,尤其是在低饱和度情况下,土中的水极难排出,因此通过实验室试验直接测试非饱和土渗透系数将耗费大量时间及人力物力。通过试验直接测定非饱和土相对渗透系数所面临的困难,严重阻碍了该问题的解决。
因此,需要提供一种非饱和土相对渗透系数的计算方法,以克服常规实验直接测定渗透系数的难度。
发明内容
本发明要解决的技术问题是提供一种预测不同初始孔隙比条件下非饱和土相对渗透系数的计算方法,从而克服在试验直接测定上所面临的困难,为获得实际工程中非饱和土渗流特性提供有益参考,进而可以大幅减小实际工程中的试验工作量,有效提高工程进度。
为解决上述技术问题,本发明采用下述技术方案
一种非饱和土相对渗透系数的计算方法,该方法包括
S1、确定待测土样类型,选定初始空隙比作为参考状态;
S2、利用方程求得参考状态下的每个饱和度对应的吸力值
S3、根据数值分析法,利用公式 其中,f(x)=[x(1-x)ζ]-1,0≤ζ≤1,和参考状态下的饱和度求得当前状态的饱和度Se值,并获得(sref,Se)曲线图;
S4、根据当前状态下的试验点和(sref,Se)曲线图比较,对参数ζ进行优化,获得最优参数ζ;
S5、保持参考状态不变,将最优参数ζ代入步骤S3中的公式,求得待测土样在任意初始孔隙比e0条件下所对应的土水特征曲线,利用该曲线标定出该条件下中的参数a和b的值;
S6、将步骤S5中标定出的参数a2和b2的值代入公式 中,计算得出任意吸力值s时,该土样的相对渗透系数值。
优选的,所述步骤S2中求得参考状态下的每个饱和度对应的吸力值的步骤包括
S21、利用参考状态的SWCC标定出方程中参数a1和b1;
S22、将参考状态下的SWCC沿有效饱和度轴进行m等分,并根据方程计算出每一个饱和度对应的吸力值
优选的,所述步骤S4中优化参数ζ的步骤包括
S41、在0至1范围内对参数ζ的值进行多次迭代,并计算出对应的当前状态下的饱和度Se值及(sref,Se)曲线图;
S42、直到(sref,Se)曲线图与初始状态下的试验点的拟合曲线相吻合,此时结束迭代,获得最优参数ζ。
本发明的有益效果如下:
采用本发明所述技术方案可有效克服在实验直接测定渗透系数上的难度,从而降低人力、物力及时间上的消耗,同时大幅度减小实际工程中的工作量,有效提高工作进度。
附图说明
下面结合附图对本发明的具体实施方式作进一步详细的说明;
图1示出一种非饱和土相对渗透系数计算方法的示意图;
图2示出Touchet粉砂壤土van Genuchten模型参数a,b标定的示意图;
图3示出Touchet粉砂壤土计算优化参数ζ选定的示意图;
图4示出e0=0.653时Touchet粉砂壤土模型预测结果与试验结果对比;
图5示出e0=0.733时Touchet粉砂壤土模型预测结果与试验结果对比;
图6示出e0=0.815时Touchet粉砂壤土模型预测结果与试验结果对比。
具体实施方式
本发明公开了一种预测不同初始孔隙比条件下非饱和土相对渗透系数的计算方法,该方法包括
确定待测土样类型,选定某初始孔隙比作为参考状态,利用该参考状态下的土水特征曲线(SWCC)标定出van Genuchten SWCC公式1参数a1,b1;
在本发明中,所提出的不同初始孔隙比下SWCC曲线的预测方程,给出了在同一吸力下不同初始孔隙比和有效饱和度的一一对应关系。所以,将利用公式1中得到的初始孔隙比土样的土水特征曲线沿有效饱和度轴等分为“m”等分,m≥100时即可满足精度要求,并利用式1计算出每一个饱和度对应的吸力值
选定初始孔隙比为当前状态,并利用该参考状态时的SWCC试验数据标定参数ζ。预先设定一个ζ的值(0≤ζ≤1),基于步骤(2)中被“m”等分出的值,利用式2计算出初始孔隙比的土样所对应的Se值,并绘制(sref,Se)曲线图,重新设定新的ζ值,并重复上述步骤,直到预测出的(sref,Se)曲线图能与初始孔隙比状态下的试验点较好地吻合,选定此时的ζ值为目标值;
其中,f(x)=[x(1-x)ζ]-1。
将标定出的最优参数ζ值代入公式2中,并仍以初始孔隙比为参考状态,在已知前提下,可计算出上述被等分的每一个值所对应的值;再利用m个点(sref,根据公式1标定出初始孔隙比条件下的vanGenuchten SWCC方程参数a2,b2;
将步骤(5)中得到的参数a2,b2代入公式3中:
计算得出任意吸力值时该土样的相对渗透系数值。
下面根据实具体施例对本发明作进一步描述:
本实施例所用土样为Touchet粉砂壤土的重塑土样,其作为一种粗粒粉土,含砂32%、粉土53%、粘土15%,土粒密度为2.599 g/cm3。对土样风干、过筛,装入制样筒进行振动,以获得不同初始孔隙比的试样。分别测试不同初始孔隙比条件下Touchet粉砂壤土的土水特征曲线干燥段、非饱和相对渗透系数。
初始孔隙比时Touchet粉砂壤土土水特征曲线干燥段试验结果如表1所示。
表1初始孔隙比时Touchet粉砂壤土土水特征曲线干燥段试验结果
首先选定初始孔隙比为参考状态,利用表1所示的该状态下的土水特征曲线(SWCC)试验点,对van Genuchten SWCC公式1中参数a1、b1进行标定,
本实施例中,参数a1和b1标定如图1所示,标定结果为a1=7.537、b1=3.153。
将利用公式1得到的参考状态时的土水特征曲线,沿有效饱和度轴等分为“m”等分,并利用式1计算出每一个饱和度对应的吸力值sref,本实施例中,取m=(0.999-0.001)/0.001;
选定初始孔隙比为参考状态,初始孔隙比为当前状态,并预先设定一个ζ的值(0≤ζ≤1)。由于本发明所提出的不同初始孔隙比下SWCC曲线的预测方程,是在同一吸力下不同初始孔隙比和有效饱和度的一一对应关系,因此,基于步骤(2)中被“m”等分出的值,利用式2计算出当前状态的Se值,并绘制(sref,Se)曲线图,重新选取ζ值后,计算出相应的饱和度值,并绘制(sref,Se)曲线图,直到绘制出的(sref,Se)曲线与初始孔隙比状态下的试验点能够较好地吻合,选定此时的ζ值为目标值,本实施例中,最优参数ζ=0.001,参数ζ的标定优化过程如图2所示,
其中,f(x)=[x(1-x)ζ]-1。
将标定出的最优参数ζ值代入公式2中,并保持初始孔隙比为参考状态不变,可计算出该土样在任意初始孔隙比e0条件下所对应的土水特征曲线,本实施例中,分别计算e0=0.653,0.733,0.815时与sref相对应的Se,得到各不同初始孔隙比条件下的(sref,Se)数据。
基于上述步骤预测出的初始孔隙比e0=0.653,e0=0.733,e0=0.815条件下相应的土水特征曲线数据(sref,Se),利用公式1标定出上述初始孔隙比条件下的van Genuchten SWCC方程参数a2,b2,标定结果如表2所示。
表2 初始孔隙比为0.653,0.733和0.815时van Genuchten SWCC方程参数a,b拟合结果
将标定的参数a2和b2代入公式3中:
计算得出该土样在初始孔隙比e0=0.653,0.733,0.815条件下的相对渗透系数值随吸力变化关系,结果如图3-5所示。通过与试验结果的比较可以看出,基于本发明方法的预测结果与试验结果吻合较好。
综上所述,通过本发明所述方案能够克服在试验直接测定上人力、物力以及时间耗费的困难,为获得实际工程中非饱和土渗流特性提供有有效可靠的方法,进而可以大幅减小实际工程中的试验工作量,有效提高工程进度。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定,对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动,这里无法对所有的实施方式予以穷举,凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。
Claims (3)
1.一种非饱和土相对渗透系数的计算方法,其特征在于:该方法包括
S1、确定待测土样类型,选定初始空隙比作为参考状态;
S2、利用方程求得参考状态下的每个饱和度对应的吸力值
S3、根据数值分析法,利用公式 其中,f(x)=[x(1-x)ζ]-1,0≤ζ≤1,和参考状态下的饱和度求得当前状态的饱和度Se值,并获得(sref,Se)曲线图;
S4、根据当前状态下的试验点和(sref,Se)曲线图比较,对参数ζ进行优化,获得最优参数ζ;
S5、保持参考状态不变,将最优参数ζ代入步骤S3中的公式,求得待测土样在任意初始孔隙比e0条件下所对应的土水特征曲线,利用该曲线标定出该条件下 中的参数a和b的值;
S6、将步骤S5中标定出的参数a2和b2的值代入公式 中,计算得出任意吸力值s时,该土样的相对渗透系数值。
2.根据权利要求1所述的一种非饱和土相对渗透系数的计算方法,其特征在于:所述步骤S2中求得参考状态下的每个饱和度对应的吸力值的步骤包括
S21、利用参考状态的SWCC标定出方程中参数a1和b1;
S22、将参考状态下的SWCC沿有效饱和度轴进行m等分,并根据方程计算出每一个饱和度对应的吸力值
3.根据权利要求1所述的一种非饱和土相对渗透系数的计算方法,其特征在于:所述步骤S4中优化参数ζ的步骤包括
S41、在0至1范围内对参数ζ的值进行多次迭代,并计算出对应的当前状态下的饱和度Se值及(sref,Se)曲线图;
S42、直到(sref,Se)曲线图与初始状态下的试验点的拟合曲线相吻合,此时结束迭代,获得最优参数ζ。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410199064.9A CN103971002A (zh) | 2014-05-12 | 2014-05-12 | 一种非饱和土相对渗透系数的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410199064.9A CN103971002A (zh) | 2014-05-12 | 2014-05-12 | 一种非饱和土相对渗透系数的计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103971002A true CN103971002A (zh) | 2014-08-06 |
Family
ID=51240490
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410199064.9A Pending CN103971002A (zh) | 2014-05-12 | 2014-05-12 | 一种非饱和土相对渗透系数的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103971002A (zh) |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104182647A (zh) * | 2014-09-03 | 2014-12-03 | 塔里木大学 | 获得土壤水分特征曲线Van Genchten模型参数的方法 |
CN104198317A (zh) * | 2014-09-03 | 2014-12-10 | 塔里木大学 | 测定土壤物理参数s的方法 |
CN105606793A (zh) * | 2016-03-07 | 2016-05-25 | 长沙理工大学 | 一种基于亚甲蓝值的细集料土水特征曲线预估方法 |
CN107063967A (zh) * | 2017-04-21 | 2017-08-18 | 湖北工业大学 | 一种基于微观孔隙通道的饱和∕非饱和渗透系数预测方法 |
CN107449706A (zh) * | 2017-06-06 | 2017-12-08 | 湖北工业大学 | 基于分形理论的变形土体饱和、非饱和渗透系数预测方法 |
CN108052783A (zh) * | 2018-01-29 | 2018-05-18 | 济南大学 | 一种基于自适应步长的非饱和土动力数值计算方法 |
CN108844881A (zh) * | 2018-08-06 | 2018-11-20 | 湖北工业大学 | 一种基于vg模型预测非饱和土相对渗透系数的方法 |
CN109238911A (zh) * | 2018-09-10 | 2019-01-18 | 湖北工业大学 | 一种基于土中水分蒸发曲线预测土-水特征曲线vg模型的方法 |
CN109325206A (zh) * | 2018-09-10 | 2019-02-12 | 柳创新 | 一种降雨径流模型参数优化方法 |
CN109342295A (zh) * | 2018-11-28 | 2019-02-15 | 东南大学 | 一种基于流固耦合分析的透水沥青路面渗透性能评价方法 |
CN110023769A (zh) * | 2016-11-08 | 2019-07-16 | 香港科技大学 | 测量非饱和土壤电阻率各向异性的电阻率测量单元 |
CN110688613A (zh) * | 2019-09-11 | 2020-01-14 | 四川大学 | 基于数值结果挖掘的确定松散颗粒材料渗透性的方法 |
CN110717693A (zh) * | 2019-10-24 | 2020-01-21 | 东南大学 | 一种污染物作用下土质防渗工程屏障渗透系数的确定方法 |
CN111220792A (zh) * | 2020-01-22 | 2020-06-02 | 长安大学 | 一种非饱和黄土水分入渗深度计算方法 |
CN111651906A (zh) * | 2020-07-10 | 2020-09-11 | 四川大学 | 一种高效的大变形尾矿固结渗透性能求解新方法 |
CN113362417A (zh) * | 2021-04-20 | 2021-09-07 | 河海大学 | 一种用拟合曲线描述裂隙率发展规律的方法 |
CN114117766A (zh) * | 2021-11-19 | 2022-03-01 | 江苏科技大学 | 一种针对vg模型的参数拟合方法 |
CN117030565A (zh) * | 2023-08-08 | 2023-11-10 | 甘肃电通电力工程设计咨询有限公司 | 一种基于土水特征曲线的非饱和土渗透系数计算方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101487787A (zh) * | 2009-02-12 | 2009-07-22 | 中国科学院武汉岩土力学研究所 | 快速确定非饱和土水力参数的方法 |
CN101813596A (zh) * | 2010-05-11 | 2010-08-25 | 中国科学院地质与地球物理研究所 | 一种测试土体非饱和气相渗透系数的方法 |
-
2014
- 2014-05-12 CN CN201410199064.9A patent/CN103971002A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101487787A (zh) * | 2009-02-12 | 2009-07-22 | 中国科学院武汉岩土力学研究所 | 快速确定非饱和土水力参数的方法 |
CN101813596A (zh) * | 2010-05-11 | 2010-08-25 | 中国科学院地质与地球物理研究所 | 一种测试土体非饱和气相渗透系数的方法 |
Non-Patent Citations (5)
Title |
---|
孙大松等: "非饱和土的渗透系数", 《水利学报》 * |
张红芬等: "非饱和土渗透系数直接试验法和间接计算法", 《人民黄河》 * |
胡冉等: "考虑变形效应的非饱和土相对渗透系数模型", 《岩石力学与工程学报》 * |
蔡国庆等: "考虑初始孔隙比影响的非饱和土相对渗透系数方程", 《HTTP://WWW.CNKI.NET/KCMS/DETAIL/32.1124.TU.20140303.1154.004.HTM1》 * |
蔡田庆等: "一种预测不同温度下非饱和土相对渗透系数的间接方法", 《岩土力学》 * |
Cited By (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104198317A (zh) * | 2014-09-03 | 2014-12-10 | 塔里木大学 | 测定土壤物理参数s的方法 |
CN104198317B (zh) * | 2014-09-03 | 2016-06-01 | 塔里木大学 | 测定土壤物理参数s的方法 |
CN104182647B (zh) * | 2014-09-03 | 2017-07-21 | 塔里木大学 | 获得土壤水分特征曲线Van Genuchten模型参数的方法 |
CN104182647A (zh) * | 2014-09-03 | 2014-12-03 | 塔里木大学 | 获得土壤水分特征曲线Van Genchten模型参数的方法 |
CN105606793A (zh) * | 2016-03-07 | 2016-05-25 | 长沙理工大学 | 一种基于亚甲蓝值的细集料土水特征曲线预估方法 |
CN105606793B (zh) * | 2016-03-07 | 2017-11-07 | 长沙理工大学 | 一种基于亚甲蓝值的细集料土水特征曲线预估方法 |
CN110023769A (zh) * | 2016-11-08 | 2019-07-16 | 香港科技大学 | 测量非饱和土壤电阻率各向异性的电阻率测量单元 |
CN107063967B (zh) * | 2017-04-21 | 2019-05-03 | 湖北工业大学 | 一种基于微观孔隙通道的饱和/非饱和渗透系数预测方法 |
CN107063967A (zh) * | 2017-04-21 | 2017-08-18 | 湖北工业大学 | 一种基于微观孔隙通道的饱和∕非饱和渗透系数预测方法 |
CN107449706A (zh) * | 2017-06-06 | 2017-12-08 | 湖北工业大学 | 基于分形理论的变形土体饱和、非饱和渗透系数预测方法 |
CN107449706B (zh) * | 2017-06-06 | 2019-11-08 | 湖北工业大学 | 基于分形理论的变形土体饱和、非饱和渗透系数预测方法 |
CN108052783A (zh) * | 2018-01-29 | 2018-05-18 | 济南大学 | 一种基于自适应步长的非饱和土动力数值计算方法 |
CN108844881B (zh) * | 2018-08-06 | 2020-08-07 | 湖北工业大学 | 一种基于vg模型预测非饱和土相对渗透系数的方法 |
CN108844881A (zh) * | 2018-08-06 | 2018-11-20 | 湖北工业大学 | 一种基于vg模型预测非饱和土相对渗透系数的方法 |
CN109325206A (zh) * | 2018-09-10 | 2019-02-12 | 柳创新 | 一种降雨径流模型参数优化方法 |
CN109238911A (zh) * | 2018-09-10 | 2019-01-18 | 湖北工业大学 | 一种基于土中水分蒸发曲线预测土-水特征曲线vg模型的方法 |
CN109325206B (zh) * | 2018-09-10 | 2023-03-24 | 柳创新 | 一种降雨径流模型参数优化方法 |
CN109342295A (zh) * | 2018-11-28 | 2019-02-15 | 东南大学 | 一种基于流固耦合分析的透水沥青路面渗透性能评价方法 |
CN109342295B (zh) * | 2018-11-28 | 2021-05-11 | 东南大学 | 一种基于流固耦合分析的透水沥青路面渗透性能评价方法 |
CN110688613A (zh) * | 2019-09-11 | 2020-01-14 | 四川大学 | 基于数值结果挖掘的确定松散颗粒材料渗透性的方法 |
CN110688613B (zh) * | 2019-09-11 | 2021-10-01 | 四川大学 | 基于数值结果挖掘的确定松散颗粒材料渗透性的方法 |
CN110717693A (zh) * | 2019-10-24 | 2020-01-21 | 东南大学 | 一种污染物作用下土质防渗工程屏障渗透系数的确定方法 |
CN110717693B (zh) * | 2019-10-24 | 2023-04-07 | 东南大学 | 一种污染物作用下土质防渗工程屏障渗透系数的确定方法 |
CN111220792A (zh) * | 2020-01-22 | 2020-06-02 | 长安大学 | 一种非饱和黄土水分入渗深度计算方法 |
CN111651906A (zh) * | 2020-07-10 | 2020-09-11 | 四川大学 | 一种高效的大变形尾矿固结渗透性能求解新方法 |
CN113362417A (zh) * | 2021-04-20 | 2021-09-07 | 河海大学 | 一种用拟合曲线描述裂隙率发展规律的方法 |
CN114117766A (zh) * | 2021-11-19 | 2022-03-01 | 江苏科技大学 | 一种针对vg模型的参数拟合方法 |
CN114117766B (zh) * | 2021-11-19 | 2024-04-26 | 江苏科技大学 | 一种针对vg模型的参数拟合方法 |
CN117030565A (zh) * | 2023-08-08 | 2023-11-10 | 甘肃电通电力工程设计咨询有限公司 | 一种基于土水特征曲线的非饱和土渗透系数计算方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103971002A (zh) | 一种非饱和土相对渗透系数的计算方法 | |
Luo et al. | Stability analysis of super-large-section tunnel in loess ground considering water infiltration caused by irrigation | |
CN107543775B (zh) | 基于分形理论确定堆石料填筑指标及现场填筑质量检测的方法 | |
CN103953023B (zh) | 一种真空联合堆载预压法沉降量的测量方法 | |
CN104794298A (zh) | 一种铀尾矿库核素对库区地下水污染风险评价的方法 | |
CN107807520A (zh) | 一种确定粗粒料填筑指标与现场质量控制的方法 | |
Fei et al. | Runoff processes and lateral transport of soil total carbon induced by water erosion in the hilly region of southern China under rainstorm conditions | |
Xin et al. | An analytical solution for predicting the transient seepage from a subsurface drainage system | |
CN104345131A (zh) | 一种稻田氮磷径流流失负荷的田间尺度估算方法 | |
Li et al. | Response of the groundwater environment to rapid urbanization in Hohhot, the provincial capital of western China | |
CN106596358A (zh) | 基于图像处理的透水混凝土集料粒径测量方法 | |
Chen et al. | Quantifying natural recharge characteristics of shallow aquifers in groundwater overexploitation zone of North China | |
Li et al. | Effect of rock exposure on runoff and sediment on karst slopes under erosive rainfall conditions | |
Wu et al. | Factors affecting the stability of soil aggregates of plinthosols in the middle reaches of the Yangtze River | |
Ma et al. | Research on the improvement of rainfall infiltration behavior of expansive soil slope by the protection of polymer waterproof coating | |
Zukowski et al. | Evaluation of field hydraulic conductivity data: Comparing spot infiltrometer test data to continuous recession data | |
Mbayaki et al. | Soil hydraulic properties of a chromic luvisol in Katumani, Kenya | |
CN104805875B (zh) | 滨海水泥土桩长期综合刚度预测方法 | |
CN101950504B (zh) | 地质力学模型分层压实风干制作与切槽埋设测试仪器方法 | |
Zhang et al. | Water infiltration and soil-water characteristics of compacted loess under applied vertical stress | |
Jabro et al. | Evaluation and comparison of five simulation models for estimating water drainage fluxes under corn | |
Rahardjo et al. | Numerical analyses for assessment of geobarrier system performance | |
Ali | Quantifying natural groundwater recharge using tracer and other techniques | |
Cen et al. | Nutrient and Stoichiometric Characteristics of Underground Fissures Soil in Karst Rocky Desertification Regions | |
Liang et al. | Rain erosion characteristics of multi-stage loess high slope of highway 1 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20140806 |