CN101762445A - 一种基于入渗时间特征参数的土壤饱和导水率测算方法 - Google Patents
一种基于入渗时间特征参数的土壤饱和导水率测算方法 Download PDFInfo
- Publication number
- CN101762445A CN101762445A CN201010011465A CN201010011465A CN101762445A CN 101762445 A CN101762445 A CN 101762445A CN 201010011465 A CN201010011465 A CN 201010011465A CN 201010011465 A CN201010011465 A CN 201010011465A CN 101762445 A CN101762445 A CN 101762445A
- Authority
- CN
- China
- Prior art keywords
- soil
- hydraulic conductivity
- infiltrating
- infiltrate
- time data
- 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.)
- Granted
Links
Images
Abstract
本发明提供一种基于入渗时间特征参数的土壤饱和导水率测算方法,该方法是在对对粘土、壤土和砂土等三种土壤质地条件下Philip-Dunne入渗数据分析的基础上,建立了土壤饱和导水率与入渗时间数据的特征参数(偏度、峰度和变异系数)的定量关系,利用该方法测算土壤饱和导水率时不需要考虑入渗前后土壤的含水率变动,只需观测Philip-Dunne入渗仪内不同水位对应的累积入渗时间,因而简化了饱和导水率的计算过程并提高了测量效率。
Description
技术领域:
本发明涉及土壤物理学领域,特别涉及一种基于入渗时间特征参数的土壤饱和导水率测算方法。
背景技术:
土壤的饱和导水率ks是土壤被水饱和时,单位水势梯度下,单位时间内通过单位面积的水量。土壤饱和导水率是土壤重要的物理性质之一。它是估计土壤非饱和导水率,计算土壤剖面中水的通量和设计灌溉、排水系统工程的一个重要土壤参数。同时,有关土壤水分传导参数的获取和掌握对于了解和治理土壤水和地表水污染也是非常重要的。因此,迅速和准确地测定饱和导水率具有十分重要的意义。土壤饱和导水率的大小,决定于土壤的性质和水的黏度,所以土壤的质地、结构、孔隙状况、盐分含量、温度、有机质含量以及土地利用方式等因素均可对其产生影响。土壤饱和导水率可在田间测定,也可采取土样在室内测定。由于获取土样过程对土样不可避免地施加一定程度的扰动,从这些扰动的土样得到的水分传导参数并不能准确地代表田间的真实情况,因而在田间进行现场入渗试验的结果更加可靠。
目前,国内外学者提出了多种测量土壤饱和导水率的方法及其计算模型,诸如单环入渗仪、双环入渗仪、Guelph入渗仪、张力盘式入渗仪、Philip-Dunne入渗仪、压力入渗仪等。上述测量方式中,以Philip-Dunne入渗仪构造最为简单,只需要一根一定直径(一般5cm左右)和一定长度(一般30cm左右)透明的塑料直管。与其它测量方法相比较,Philip-Dunne入渗方法测量饱和导水率时用水量最少且测量时间较短,再加上该设备成本低廉而非常适宜于野外大面积的饱和导水率测量,特别适合于运输不便的偏远地区土壤参数测量工作。1993年John Philip基于Green-Ampt模型基本假设提出了Philip-Dunne入渗方法,在降水头多维入渗过程中利用少量水便可快速测定饱和导水率ks和Green-Ampt模型中的湿润锋吸力Ψ。为推广Philip-Dunne入渗方法的运用,De Haro做了饱和导水率ks的灵敏度分析,Muuoz-Carpena对比了Philip-Dunne入渗和常水头井渗透,并对计算方法做了改进,得到Philip-Dunne入渗的参数计算公式:
公式①对土壤饱和导水率而言是隐式方程,求解过程十分复杂,并且需要测量入渗前后土壤体积含水率的变动。不考虑入渗前后土壤含水率的变动,基于数据统计分析理论,分析Philip-Dunne入渗过程的入渗时间数据的特征参数值与土壤饱和导水率的关系,利用入渗时间的偏度、峰度和变异系数来推求土壤饱和导水率的方法至今未见报道。
发明内容:
本发明的目的是克服上述已有技术的不足,而提供一种基于入渗时间特征参数的土壤饱和导水率测算方法,主要解决了现有方法需要测量入渗前后土壤体积含水率的变动以及计算过程冗长、复杂的问题。
为了达到上述目的,本发明是这样实现的:一种基于入渗时间特征参数的土壤饱和导水率测算方法,其特殊之处在于包括以下步骤:
(1)在待测区域选定的测试点进行Philip-Dunne入渗试验,从入渗开始记录入渗过程中Philip-Dunne入渗仪内特定水位对应的累积入渗时间数据,并测量每个测试点入渗前后的土壤含水率;
(2)利用Muuoz-Carpena公式结合第(1)步所测的每个测试点入渗前后的土壤含水率及特定水位对应的累积入渗时间数据,计算出每个测试点的土壤饱和导水率;
(3)利用第(1)步所测特定水位对应的累积入渗时间数据,计算出每个测试点入渗时间数据的特征参数,结合第(2)步计算出的每个测试点对应的土壤饱和导水率,建立土壤饱和导水率和入渗时间数据的特征参数之间的经验关系;
(4)在待测区未测位置处进行Philip-Dunne入渗试验,只需观测入渗过程中Philip-Dunne入渗仪内特定水位对应的累积入渗时间数据,并求出每个测试点的入渗时间数据的特征参数;
(5)把第(4)步求得的入渗时间数据的特征参数代入第(3)步获得的土壤饱和导水率和入渗时间数据的特征参数之间的经验关系,估算出未测位置的土壤饱和导水率。
本发明的一种基于入渗时间特征参数的土壤饱和导水率测算方法,其所述的入渗时间数据的特征参数指入渗时间数据的偏度、峰度和变异系数。
所述的入渗时间是指从入渗开始计时Philip-Dunne入渗仪内特定水位对应的入渗累积时间。
所述的特定水位必须包括中间水位及最低水位,并且所有的特定水位是一致的。
本发明的一种基于入渗时间特征参数的土壤饱和导水率测算方法,可用于野外利用Philip-Dunne入渗仪测定土壤饱和导水率参数的场合。对于Philip-Dunne入渗,建立了一个土壤饱和导水率ks与入渗时间数据的特征参数值的定量关系,二者之间存在指数函数关系,该关系的一般形式为:
ks=a·eb·s; ②
式中ks(cm/min)为饱和导水率,e≈2.71828为自然对数的底数,s为Philip-Dunne入渗过程中入渗时间的特征参数值,a、b为拟合的系数。根据具体情况也可采用其它函数关系。
本发明的一种基于入渗时间特征参数的土壤饱和导水率测算方法与已有技术相比具有突出的实质性特点和显著进步:1、利用Philip-Dunne入渗方法测量土壤饱和导水率时不需要测量入渗前后的土壤含水率数据,减少了测量环节,节约了测量时间和相应设备投入,极大提高了测量效率;2、所提出的估算方法形式简单计算方便易于应用,有效减少了计算时长,极大提高了计算效率。
附图说明:
图1是本发明的计算流程图。
具体实施方式:
为了更好的理解与实施,下面结合附图给出具体实施例详细说明本发明一种基于入渗时间特征参数的土壤饱和导水率测算方法。
实施例1:第一步,首先测得土壤饱和导水率参数ks和对应的入渗时间数据的特征参数值的定量关系。在待测区域的若干测试点利用Philip-Dunne入渗方法,采用Philip-Dunne入渗仪,即一根直径为5cm左右和长度35cm左右(但要长于30cm)的透明塑料直管,将Philip-Dunne入渗仪插入土壤地表以下10cm左右的深度处,并将管内的土壤取出,在室外对需要测量土壤饱和导水率的田块进行多个试验点的降水头多维入渗实验。Philip-Dunne入渗仪的初始水头保持为30cm,实验前后用TDR土壤水分速测系统或其他方法,测定实验开始前的土壤体积含水率的初始值θ0和水分全部入渗后对应的土壤含水率值θs,实验过程中用秒表或其他计时工具记录Philip-Dunne入渗仪中的水位在不同位置对应的累积入渗时间,从水位开始下降时开始计时,水位每变化5cm观测一次时间,直到水全部入渗停止计时。根据测定的土壤含水率及不同水位对应的累积入渗时间数据,利用公式①计算各试验点的土壤饱和导水率参数ks。对不同水位对应的累积入渗时间数据序列,利用EXCELL等常用数据分析软件可十分方便、快捷地计算出每个测试点累积入渗时间序列数据的偏度、峰度和变异系数等特征参数值。对所有测试点上的土壤饱和导水率参数ks和对应的入渗时间数据的特征参数值进行分析,通过回归分析拟合出入渗参数ks与入渗时间数据的特征参数值的关系,并做显著性检验,要求二者显著相关。对粘土、壤土和砂土等多种土壤质地条件下的研究结果表明(见表1,2,3),土壤饱和导水率参数ks和对应的入渗时间数据的特征参数值均存在指数函数关系,即公式②,二者在0.001置信水平上显著相关,表明本方法能够用于土壤饱和导水率参数ks的测算。
表1饱和导水率和入渗时间偏度之间指数函数拟合结果
表2饱和导水率和入渗时间峰度之间指数函数拟合结果
表3饱和导水率和入渗时间变异系数之间指数函数拟合结果
第二步,利用土壤饱和导水率参数ks和对应的入渗时间数据的特征参数值的定量关系方程,利用Philip-Dunne入渗仪测量数据测算饱和导水率。对于尚未进行测量的其他位置的土壤进行Philip-Dunne入渗试验,将Philip-Dunne入渗仪插入土壤地表以下10cm左右的深度处,并将管内的土壤取出,Philip-Dunne入渗仪的初始水头保持为30cm,实验过程中用秒表或其他计时工具记录Philip-Dunne入渗仪中的水位在不同位置对应的累积入渗时间,从水位开始下降时开始计时,水位每变化5cm观测一次时间。利用计算器或者EXCELL等常用软件可十分方便地计算出累计入渗时间序列数据的特征参数值,把得到的每个测点的入渗时间数据的特征参数值分别带入第一步的得到的土壤饱和导水率参数ks和对应的入渗时间数据的特征参数值的定量关系方程②,即可估算出每个测点对应的未测位置的土壤饱和导水率值。
Claims (4)
1.一种基于入渗时间特征参数的土壤饱和导水率测算方法,其特征在于包括以下步骤:
(1)在待测区域选定的测试点进行Philip-Dunne入渗试验,从入渗开始记录入渗过程中Philip-Dunne入渗仪内特定水位对应的累积入渗时间数据,并测量每个测试点入渗前后的土壤含水率;
(2)利用Muuoz-Carpena公式结合第(1)步所测的每个测试点入渗前后的土壤含水率及特定水位对应的累积入渗时间数据,计算出每个测试点的土壤饱和导水率;
(3)利用第(1)步所测特定水位对应的累积入渗时间数据,计算出每个测试点入渗时间数据的特征参数,结合第(2)步计算出的每个测试点的土壤饱和导水率,建立土壤饱和导水率和入渗时间数据的特征参数之间的经验关系;
(4)在待测区未测位置处进行Philip-Dunne入渗试验,只需观测入渗过程中Philip-Dunne入渗仪内特定水位对应的累积入渗时间数据,并求出每个测试点的入渗时间数据的特征参数;
(5)把第(4)步求得的入渗时间数据的特征参数代入第(3)步获得的土壤饱和导水率和入渗时间数据的特征参数之间的经验关系,估算出未测位置的土壤饱和导水率。
2.根据权利要求1所述的一种基于入渗时间特征参数的土壤饱和导水率测算方法,其特征在于所述的入渗时间数据的特征参数指入渗时间数据的偏度、峰度和变异系数。
3.根据权利要求1所述的一种基于入渗时间特征参数的土壤饱和导水率测算方法,其特征在于所述的入渗时间是指从入渗开始计时Philip-Dunne入渗仪内特定水位对应的入渗累积时间。
4.根据权利要求1或3所述的一种基于入渗时间特征参数的土壤饱和导水率测算方法,其特征在于所述的特定水位包括中间水位及最低水位,并且对不同测点所观测记录的特定水位是完全一致的。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010100114659A CN101762445B (zh) | 2010-01-15 | 2010-01-15 | 一种基于入渗时间特征参数的土壤饱和导水率测算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010100114659A CN101762445B (zh) | 2010-01-15 | 2010-01-15 | 一种基于入渗时间特征参数的土壤饱和导水率测算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101762445A true CN101762445A (zh) | 2010-06-30 |
CN101762445B CN101762445B (zh) | 2011-11-30 |
Family
ID=42493735
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010100114659A Expired - Fee Related CN101762445B (zh) | 2010-01-15 | 2010-01-15 | 一种基于入渗时间特征参数的土壤饱和导水率测算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101762445B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102692359A (zh) * | 2012-05-10 | 2012-09-26 | 张振华 | 一种土壤水气传导率原位测量系统与测量方法 |
CN103344747A (zh) * | 2013-07-03 | 2013-10-09 | 兰州大学 | 一种缩短粉质粘土饱和时间的技术方法 |
CN103822863A (zh) * | 2014-02-17 | 2014-05-28 | 同济大学 | 应用于现场钻孔注水测定土体渗透系数的仪器 |
CN104764934A (zh) * | 2015-04-27 | 2015-07-08 | 塔里木大学 | 土壤饱和浸提液电导率的测量方法 |
CN105547957A (zh) * | 2015-12-10 | 2016-05-04 | 中国科学院地理科学与资源研究所 | 土壤下渗湿润峰吸力和降雨初损量计算方法 |
CN106018196A (zh) * | 2016-05-30 | 2016-10-12 | 武汉大学 | 土颗粒团聚与离散程度的微观视频图像监测与评估方法 |
CN107703045A (zh) * | 2017-11-09 | 2018-02-16 | 石家庄铁道大学 | 海绵城市绿地雨水收集能力分析系统及分析方法 |
CN110595984A (zh) * | 2019-10-29 | 2019-12-20 | 兰州理工大学 | 一种测量原状土饱和导水率的圆筒入渗仪及其测量方法 |
CN112881167A (zh) * | 2020-09-10 | 2021-06-01 | 北京理工大学 | 一种微动疲劳裂纹萌生位置判定方法 |
CN116148445A (zh) * | 2023-04-19 | 2023-05-23 | 中国环境科学研究院 | 定量降水溅蚀对稻田生态系统土壤结构侵蚀影响的方法 |
CN107703045B (zh) * | 2017-11-09 | 2024-04-26 | 石家庄铁道大学 | 海绵城市绿地雨水收集能力分析系统及分析方法 |
-
2010
- 2010-01-15 CN CN2010100114659A patent/CN101762445B/zh not_active Expired - Fee Related
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102692359A (zh) * | 2012-05-10 | 2012-09-26 | 张振华 | 一种土壤水气传导率原位测量系统与测量方法 |
CN103344747A (zh) * | 2013-07-03 | 2013-10-09 | 兰州大学 | 一种缩短粉质粘土饱和时间的技术方法 |
CN103344747B (zh) * | 2013-07-03 | 2015-07-15 | 兰州大学 | 一种缩短粉质粘土饱和时间的技术方法 |
CN103822863A (zh) * | 2014-02-17 | 2014-05-28 | 同济大学 | 应用于现场钻孔注水测定土体渗透系数的仪器 |
CN104764934A (zh) * | 2015-04-27 | 2015-07-08 | 塔里木大学 | 土壤饱和浸提液电导率的测量方法 |
CN104764934B (zh) * | 2015-04-27 | 2019-04-26 | 塔里木大学 | 土壤饱和浸提液电导率的测量方法 |
CN105547957B (zh) * | 2015-12-10 | 2018-05-04 | 中国科学院地理科学与资源研究所 | 土壤下渗湿润峰吸力和降雨初损量计算方法 |
CN105547957A (zh) * | 2015-12-10 | 2016-05-04 | 中国科学院地理科学与资源研究所 | 土壤下渗湿润峰吸力和降雨初损量计算方法 |
CN106018196A (zh) * | 2016-05-30 | 2016-10-12 | 武汉大学 | 土颗粒团聚与离散程度的微观视频图像监测与评估方法 |
CN106018196B (zh) * | 2016-05-30 | 2019-04-09 | 武汉大学 | 土颗粒团聚与离散程度的微观视频图像监测与评估方法 |
CN107703045A (zh) * | 2017-11-09 | 2018-02-16 | 石家庄铁道大学 | 海绵城市绿地雨水收集能力分析系统及分析方法 |
CN107703045B (zh) * | 2017-11-09 | 2024-04-26 | 石家庄铁道大学 | 海绵城市绿地雨水收集能力分析系统及分析方法 |
CN110595984A (zh) * | 2019-10-29 | 2019-12-20 | 兰州理工大学 | 一种测量原状土饱和导水率的圆筒入渗仪及其测量方法 |
CN112881167A (zh) * | 2020-09-10 | 2021-06-01 | 北京理工大学 | 一种微动疲劳裂纹萌生位置判定方法 |
CN112881167B (zh) * | 2020-09-10 | 2022-01-04 | 北京理工大学 | 一种微动疲劳裂纹萌生位置判定方法 |
CN116148445A (zh) * | 2023-04-19 | 2023-05-23 | 中国环境科学研究院 | 定量降水溅蚀对稻田生态系统土壤结构侵蚀影响的方法 |
CN116148445B (zh) * | 2023-04-19 | 2023-07-28 | 中国环境科学研究院 | 定量降水溅蚀对稻田生态系统土壤结构侵蚀影响的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101762445B (zh) | 2011-11-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101762445B (zh) | 一种基于入渗时间特征参数的土壤饱和导水率测算方法 | |
CN101738356B (zh) | 一种测算土壤湿润锋吸力的方法 | |
Briggs et al. | Practical limitations on the use of diurnal temperature signals to quantify groundwater upwelling | |
CN103592690B (zh) | 基于电成像测井孔隙度谱信息自动识别储层裂缝的方法 | |
Tromp-van Meerveld et al. | Assessment of multi-frequency electromagnetic induction for determining soil moisture patterns at the hillslope scale | |
CN104374827B (zh) | 横观各向同性岩体原位动弹性模量的各向异性系数的测量方法 | |
CN104793249B (zh) | 一种系统检测海水入侵的方法 | |
CN104331639A (zh) | 一种土壤含水率间接获取及快速评价方法 | |
Águila et al. | Parametric and numerical analysis of the estimation of groundwater recharge from water-table fluctuations in heterogeneous unconfined aquifers | |
Siyuan et al. | Estimating impacts of water-table depth on groundwater evaporation and recharge using lysimeter measurement data and bromide tracer | |
Shi et al. | Measuring shallow water flow velocity with virtual boundary condition signal in the electrolyte tracer method | |
Langhans et al. | Development and parameterization of an infiltration model accounting for water depth and rainfall intensity | |
CN105138761B (zh) | 降雨条件下坡面地表糙率及土壤吸湿率的估算方法 | |
CN104881564B (zh) | 结构面粗糙度系数尺寸效应概率密度函数模型的构建方法 | |
CN205719870U (zh) | 原位土壤降雨入渗测量装置 | |
Liu et al. | A note on the time compression approximation | |
Hajhouji et al. | Dynamics of groundwater recharge near a semi-arid Mediterranean intermittent stream under wet and normal climate conditions | |
Zhu et al. | Measurement and simulation of nitrogen leaching loss in hillslope cropland of purple soil | |
Altafi Dadgar et al. | Transient potential groundwater recharge under surface irrigation in semiarid environment: An experimental and numerical study | |
Zeng et al. | Evaluation of hydraulic parameters obtained by different measurement methods for heterogeneous gravel soil. | |
Rucker | Inverse upscaling of hydraulic parameters during constant flux infiltration using borehole radar | |
CN109271675A (zh) | 单指数模型矿区生态临界地下水位动态预测方法及装置 | |
CN104179164A (zh) | 用标贯击数N探测长江河口、砂嘴、砂岛地层砂性土相对密度Dr的方法 | |
Zhang et al. | Surface-positioned double-ring to improve traditional infiltrometer for measuring soil infiltration | |
CN204085926U (zh) | 一种可控采样量的表层土壤采样器 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20111130 Termination date: 20120115 |