CN115879289A - 一种基于水分下渗过程原位快速获取土壤水力学参数的方法 - Google Patents
一种基于水分下渗过程原位快速获取土壤水力学参数的方法 Download PDFInfo
- Publication number
- CN115879289A CN115879289A CN202211501243.4A CN202211501243A CN115879289A CN 115879289 A CN115879289 A CN 115879289A CN 202211501243 A CN202211501243 A CN 202211501243A CN 115879289 A CN115879289 A CN 115879289A
- Authority
- CN
- China
- Prior art keywords
- infiltration
- soil
- probe
- tdr
- length
- 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
- 239000002689 soil Substances 0.000 title claims abstract description 104
- 230000008595 infiltration Effects 0.000 title claims abstract description 90
- 238000001764 infiltration Methods 0.000 title claims abstract description 90
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 88
- 238000000034 method Methods 0.000 title claims abstract description 83
- 238000011065 in-situ storage Methods 0.000 title claims abstract description 18
- 239000000523 sample Substances 0.000 claims description 41
- 238000009736 wetting Methods 0.000 claims description 19
- 229920006395 saturated elastomer Polymers 0.000 claims description 14
- 238000012544 monitoring process Methods 0.000 claims description 10
- 238000012360 testing method Methods 0.000 claims description 8
- 238000005259 measurement Methods 0.000 claims description 7
- 238000009825 accumulation Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 5
- 206010016807 Fluid retention Diseases 0.000 claims description 4
- 238000012897 Levenberg–Marquardt algorithm Methods 0.000 claims description 4
- 230000001186 cumulative effect Effects 0.000 claims description 4
- 238000005520 cutting process Methods 0.000 claims description 4
- 230000000717 retained effect Effects 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 6
- 238000005457 optimization Methods 0.000 abstract description 4
- 238000009795 derivation Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 abstract description 2
- 239000012528 membrane Substances 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 244000035744 Hura crepitans Species 0.000 description 1
- 229920005372 Plexiglas® Polymers 0.000 description 1
- 238000005273 aeration Methods 0.000 description 1
- 238000007605 air drying Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- HPNSNYBUADCFDR-UHFFFAOYSA-N chromafenozide Chemical compound CC1=CC(C)=CC(C(=O)N(NC(=O)C=2C(=C3CCCOC3=CC=2)C)C(C)(C)C)=C1 HPNSNYBUADCFDR-UHFFFAOYSA-N 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 238000000227 grinding Methods 0.000 description 1
- 230000002262 irrigation Effects 0.000 description 1
- 238000003973 irrigation Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 239000004926 polymethyl methacrylate Substances 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000007873 sieving Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
- G01N33/246—Earth materials for water content
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
- G01N33/245—Earth materials for agricultural purposes
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Food Science & Technology (AREA)
- Analytical Chemistry (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Medicinal Chemistry (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明公开了一种基于水分下渗过程原位快速获取土壤水力学参数的方法,采用全新推导方法,基于Richards水流控制方程与Brooks‑Corey模型推导得出可以准确描述恒定水头边界条件下均质土壤一维水分下渗的解析解,首次刻画了饱和区随时间增长这一重要入渗特性,可以更准确地模拟土壤水分下渗过程。基于解析解的逆过程,一种参数优化估计方法被提出。该方法仅需通过时域反射仪(TDR)测量积水下渗过程中累积入渗量和湿润区长度随时间的变化数据,便可在田间原位快速获取土壤水力特性。该计算方法成本低,反演过程稳定,避免了参数多解、不收敛等问题,为大面积流域水文模拟所需的参数提供了强有力的获取手段。
Description
技术领域
本发明涉及到一种基于水分下渗过程原位快速获取土壤水力学参数的方法,属于农田水利和土壤水文测量技术领域。
背景技术
土壤水力学特性参数(土壤水分特征曲线、饱和/非饱和导水率)是模拟包气带土壤水分运动的基本输入参数。由于土壤水力学性质具有相对较大的时空变异特性,在一个地区或流域内确定有代表性的、准确的土壤水力学参数数值是一项巨大而耗时的工作。因此,开发简单、低成本、省时的方法对研究一个地区或流域内的土壤水力特性至关重要。
传统的土壤水力学参数的获取方法通常依赖于环刀采样,进行室内分析,该方法需要使用昂贵的设备(如:沙箱、吸力平板仪、压力膜仪、离心机等),通过缓慢的排水过程来获取参数,并且饱和导水率需要通过额外的定水头或变水头法单独测量,耗时费力,成本高昂。垂直入渗或水平吸渗试验是较为快速的土壤水分运动过程,常被用于间接快速获取土壤水力学特性参数。此类方法多采用数值反演途径来获取参数,存在数值解不收敛和多解问题。基于近似解的反演是当前最有潜力的发展方向,目前关于下渗问题的近似解很多,但其精度尚不足以支持土壤水力学特性参数反演预测的需求,其中一个重要的原因在于目前的入渗公式中没有考虑土壤水分剖面在下渗过程中饱和区随时间增长这一重要特性。此外,土壤结构对于水力特性特别重要,但用环刀采样原状土壤样品带回室内测量在采样和运输过程中难免破坏原状土样的结构,造成测量结果失真。因此,亟需发展一种可在田间原位直接应用的快速获取土壤水力学特性参数的方法。
发明内容
本发明的目的是克服现有技术的不足,提供一种基于水分下渗过程原位快速获取土壤水力学参数的方法。
本发明采用的技术方案为:
一种基于水分下渗过程原位快速获取土壤水力学参数的方法,其步骤包括:
(4)测量获得土壤的饱和含水量θs、初始含水量θi、滞留含水量θr;
(5)针对积水条件下田间原位一维水分下渗过程,测量获得累积入渗量I和湿润区长度zf与入渗时间t的变化关系;
(6)根据目标函数(9)优化计算待求参数n、hd和Ks
优选的,所述步骤(3)中,采用Levenberg-Marquardt算法优化目标函数(9),其算法思想在于通过估计正确的待求参数n、hd和Ks,使得目标函数(9)达到最小值,式中的估计值和/>通过以下方法获得:
上式中未知参数的定义为:
优选的,步骤(1)和(2)的参数获取方法具体为:
入渗前将TDR探针插入待测土壤,进行水头高度为hp的积水入渗试验,从入渗前即开始收集TDR探针的时域反射信号,计算出土壤初始含水量θi和饱和含水量θs,对土壤进行环刀取样,测量其滞留含水量θr,因为滞留含水量θr仅由土壤质地所决定,因此可在入渗结束后进行原位破坏性采样,室内采用压力膜仪进行15bar下的排水试验,通过称重法测定获得θr;
在积水入渗的过程中通过解译时域TDR的波形图实时获取湿润锋的时间序列数据zf~t,以及累积入渗量的时间序列数据I~t。
优选的,步骤(1)和(2)的参数获取方法更具体为:
A、在水入渗前,将长度为L的TDR探针垂直的完全插入待测土壤,直到探针首部与土壤齐平,水入渗开始前即用时域反射仪以一定频率测量收集TDR探针反射信号,持续收集到湿润锋超过探针长度,探针监测土层达到饱和状态;
B、入渗前和入渗结束时TDR所测区间没有湿润锋,从TDR反射信号波形中提取入渗前探针首部和尾部的表观位置,进而获取初始TDR探针监测土层的表观长度Lad,以及入渗结束时探针首部和尾部的表观位置,进而获取饱和时TDR探针监测土层的表观长度Las,利用电磁波传输理论计算监测土层初始平均介电常数Kad=(Lad/L)2和饱和时的平均介电常数Kas=(Las/L)2,然后利用土壤含水量与介电常数之间的关系将Kad和Kas转换为土壤初始含水量θi和饱和含水量θs;
C、在入渗开始后,湿润锋在TDR所测区间,每个时刻从湿润锋到达探针尾部前的TDR反射信号波形TDR探针监测土层的表观长度La和湿润区的表观长度Law。根据La计算得到监测土层的平均介电常数Ka=(La/L)2,利用土壤含水量与介电常数之间的关系将TDR探针监测土层的平均介电常数Ka转换为实时平均含水量θ,进而计算出土壤累积入渗量I=L*(θ-θi),得到累积入渗量的时间序列数据I~t,以及湿润区长度zf=L*(Lad-La+Law)/Lad,得到湿润区长度的时间序列数据zf~t。
优选的,步骤(C)中湿润锋和探针尾部在TDR波形中的表观位置采用二阶导数法提取,具体如下:对波形数据分别求一阶和二阶导数,在一阶导数最大值出现位置之后,确定波形数据一阶导数的第一个零点即为TDR探针首部对应的表观位置Xh;二阶导数在Xh之后出现的第二个和第三个局部最大值出现的位置分别是湿润锋对应的表观位置Xf和TDR探针尾部对应的表观位置Xe,监测土层的表观长度La为探针首部和尾部表观位置之间的距离Xe-Xh;湿润区的表观长度Law为探针首部和尾部表观位置之间的距离Xf-Xh。
本发明的有益效果如下:
与已有技术相比,本发明采用全新推导方法,利用Richards水流控制方程与Brooks-Corey模型推导出恒定水头边界条件下均质土壤一维水分下渗问题的解析解,首次详细刻画了饱和区随时间增长这一重要入渗特性,可以更准确地描述土壤水分下渗过程。本发明通过积水条件下均质土一维水分下渗解析解的逆过程,发展了一种土壤水力学参数获取方法,结合TDR技术,仅需在田间原位进行的积水入渗试验中利用时域反射仪(TDR)实时记录反射波形信息,便可解译得到累积入渗量和湿润区长度随时间的变化数据,进而获得完整的土壤水力学参数。该方法进行的参数反演过程稳定,计算成本低,避免了参数多解、不收敛等问题,为大面积高效获取流域或区域水文地球化学循环过程模型模拟所需的参数提供了强有力的手段。
附图说明
图1是本发明的计算流程图;
图2是TDR法和常规方法观测的湿润区长度变化对比图;
图3是解析解模拟和试验观测的累积入渗量变化对比图;
图4是本发明方法预测的土壤水分特征曲线与压力膜法实测结果的比较。
具体实施方式
以下结合附图,进一步说明本发明的内容,但不应理解为对本发明的限制。在不背离本发明精神和实质的情况下,对本发明方法、步骤或条件所作的修改和替换,均属于本发明的范围。若未特别指明,以下实施案例中所用的技术手段为本领域技术人员所熟知的常规手段。图1说明了本发明的计算过程。
实施例1
技术原理:基于理查德水流控制方程与Brooks-Corey模型,通量假定以及积分中值定理,获得了描述土壤水分下渗的剖面方程,再基于质量守恒原理与剖面方程推导得出完整的描述恒定水头边界条件下均质土一维水分下渗的解析模型。其中:
理查德水流控制方程及其初始与边界条件表达式如下:
式中:θ和θi分别为土壤含水量和初始土壤含水量(cm3cm-3),t为入渗时间,z为垂向坐标,K为土壤导水率(cm min-1),h为土壤基质吸力(cm),hp是土壤表层的积水深度(cm)。
所述采用的Brooks-Corey模型描述土壤非饱和导水率、土壤基质势和土壤含水量之间关系的表达式分别如下:
式中:θs和θr分别为土壤饱和含水量和滞留含水量(cm3cm-3),Ks是土壤饱和导水率(cm min-1),n是土壤孔隙尺寸分布指数,hd是土壤进气吸力值(cm),m=3n+2。
公式(4)描述了一维水分下渗的水分剖面中有效土壤水饱和度S随土壤深度z(cm)的变化关系,式(4)中zf为湿润区长度(cm),zs为饱和区长度(cm),
公式(5)描述了土壤水分剖面在下渗过程中饱和区长度zs随入渗时间t(min)的定量关系,
公式(6)描述了土壤水分剖面在下渗过程中饱和区长度zs(cm)与湿润区长度zf(cm)的定量关系,
公式(7)描述了土壤水分剖面在下渗过程中饱和区长度zs(cm)与累积入渗量I(cm)的定量关系,
式中:θi是初始含水量,Ki是土壤初始含水量下的导水率(cm min-1),其余未知参数定义为,
通过上述解析模型的逆过程,获得一种对土壤水力学参数(n,hd和Ks)的优化计算方法,所述方法步骤如下:
步骤1,针对恒定水头条件下(hp)田间原位一维水分下渗过程,测量获得累积入渗量I和湿润区长度zf随入渗时间t的变化数据;
步骤2,测量获得土壤的饱和含水量θs、初始含水量θi、滞留含水量θr,根据上述参数以及测定获得的数据集I~t和zf~t,通过目标函数(9)优化计算待求参数n,hd和Ks:
目标函数(9)采用Levenberg-Marquardt算法优化得到,其算法思想在于通过估计正确的待求参数(n,hd和Ks),使得目标函数(9)达到最小值,式中的估计值和通过以下方法获得:对于给定的一组估计参数(ne,hd,e和Ks,e),对任一入渗时间t,可通过公式(5),采用牛顿-辛普森方法计算获得饱和区长度zs的估计值/>
式中:参数定义同上。
实施例2
步骤1:采集江西鹰潭红壤(北纬28.202942°,东经116.948483°,砂粒20.9%,粉粒34.9%,粘粒44.2%),风干碾碎过2mm筛,按1.29g cm-3的容重装土到内径19cm,长30cm的有机玻璃柱中。入渗试验前,首先在土壤表面插入30cm长的时域TDR探针,通过时域TDR测定得到待测土壤的初始含水量θi=0.058cm cm-3,因初始含水量比较低,滞留含水量暂以初始含水量近似,即θr=0.058cm cm-3。实际应用时,若需更准确的滞留含水量值,可在研究区0.5m以内通过环刀采样,带回实验室内利用压力膜仪法测量获得滞留含水量。
步骤2:随后采用行业通用的马氏瓶提供恒压供水,进行水头高度为hp=1cm的积水入渗试验,并配备水位自动监测和采集系统,在积水入渗过程中通过时域TDR获取所测土壤的实时反射波形信号,通过波形分析解译获取累积入渗量和湿润区长度随时间变化的数据序列I~t和zf~t(图2和3中的观测值)。入渗试验结束后继续供水直至土柱饱和,然后用TDR测量土壤饱和含水量θs=0.513cm cm-3。
步骤3,根据已知参数θi、θs、θr和hp,以及测量得到的累积入渗量的时间序列数据I~t和湿润锋的时间序列数据zf~t,利用Levenberg-Marquardt算法优化公式(9),进而同时估计得到参数n、hd和Ks。其估计结果与标准压力膜仪法和定水头法的实测值如表1和图4所示,从表中可以看出采用本发明估计土壤水力学参数的方法是一种准确的方法。
表1基于本发明估计得到的土壤水力学参数与土壤实测水力学参数之间的比较
n | h<sub>d</sub>(cm) | K<sub>s</sub>(cm min<sup>-1</sup>) | |
实测值 | 0.17 | 1.2 | 0.079 |
估计值 | 0.1965 | 2.588 | 0.081 |
。
最后所应说明的是,以上具体实施方式仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对发明进行了详细说明,本领域的普通技术人员应当理解,可以对发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (5)
3.根据权利要求1所述的基于水分下渗过程原位快速获取土壤水力学参数的方法,其特征在于步骤(1)和(2)的参数获取方法具体为:
入渗前将TDR探针插入待测土壤,进行水头高度为hp的积水入渗试验,从入渗前即开始收集TDR探针的时域反射信号,计算出土壤初始含水量θi和饱和含水量θs,对土壤进行环刀取样,测量其滞留含水量θr;
通过解译积水入渗过程中TDR的反射波形实时获取湿润区长度的时间序列数据zf~t,以及累积入渗量的时间序列数据I~t。
4.根据权利要求3所述的基于水分下渗过程原位快速获取土壤水力学参数的方法,其特征在于步骤(1)和(2)的参数获取方法更具体为:
A、在水入渗前,将长度为L的TDR探针垂直的完全插入待测土壤,直到探针首部与土壤齐平,水入渗开始前即用时域反射仪以一定频率测量收集TDR探针反射信号,持续收集到湿润锋超过探针长度,探针监测土层达到饱和状态;
B、入渗前和入渗结束时TDR所测区间没有湿润锋,从TDR反射信号波形中提取入渗前探针首部和尾部的表观位置,进而获取初始TDR探针监测土层的表观长度Lad,以及入渗结束时探针首部和尾部的表观位置,进而获取饱和时TDR探针监测土层的表观长度Las,利用电磁波传输理论计算监测土层初始平均介电常数Kad=(Lad/L)2和饱和时的平均介电常数Kas=(Las/L)2,然后利用土壤含水量与介电常数之间的关系将Kad和Kas转换为土壤初始含水量θi和饱和含水量θs;
C、在入渗开始后,湿润锋在TDR所测区间,每个时刻从湿润锋到达探针尾部前的TDR反射信号波形中提取TDR探针监测土层的表观长度La,和湿润区的表观长度Law,根据La计算得到监测土层的平均介电常数Ka=(La/L)2,利用土壤含水量与介电常数之间的关系将TDR探针监测土层的平均介电常数Ka转换为实时平均含水量θ,进而计算出土壤累积入渗量I=L*(θ-θi),得到累积入渗量的时间序列数据I~t,以及湿润区长度zf=L*(Lad-La+Law)/Lad,得到湿润区长度的时间序列数据zf~t。
5.根据权利要求4所述的基于水分下渗过程原位快速获取土壤水力学参数的方法,其特征在于:步骤(C)中湿润锋和探针尾部在TDR波形中的表观位置采用二阶导数法提取,具体如下:对波形数据分别求一阶和二阶导数,在一阶导数最大值出现位置之后,确定波形数据一阶导数的第一个零点即为TDR探针首部对应的表观位置Xh;二阶导数在Xh之后出现的第二个和第三个局部最大值出现的位置分别是湿润锋对应的表观位置Xf和TDR探针尾部对应的表观位置Xe;监测土层的表观长度La为探针首部和尾部表观位置之间的距离Xe-Xh;湿润区的表观长度Law为探针首部和尾部表观位置之间的距离Xf-Xh。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211501243.4A CN115879289A (zh) | 2022-11-28 | 2022-11-28 | 一种基于水分下渗过程原位快速获取土壤水力学参数的方法 |
US18/242,011 US20240183838A1 (en) | 2022-11-28 | 2023-09-05 | Method for quickly acquiring soil hydraulic properties in situ based on ponded infiltration experiment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211501243.4A CN115879289A (zh) | 2022-11-28 | 2022-11-28 | 一种基于水分下渗过程原位快速获取土壤水力学参数的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115879289A true CN115879289A (zh) | 2023-03-31 |
Family
ID=85764297
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211501243.4A Pending CN115879289A (zh) | 2022-11-28 | 2022-11-28 | 一种基于水分下渗过程原位快速获取土壤水力学参数的方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US20240183838A1 (zh) |
CN (1) | CN115879289A (zh) |
-
2022
- 2022-11-28 CN CN202211501243.4A patent/CN115879289A/zh active Pending
-
2023
- 2023-09-05 US US18/242,011 patent/US20240183838A1/en active Pending
Also Published As
Publication number | Publication date |
---|---|
US20240183838A1 (en) | 2024-06-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ward et al. | Laboratory measurements of solute transport using time domain reflectometry | |
Noborio | Measurement of soil water content and electrical conductivity by time domain reflectometry: a review | |
Greco | Soil water content inverse profiling from single TDR waveforms | |
Rothe et al. | Changes in soil structure caused by the installation of time domain reflectometry probes and their influence on the measurement of soil moisture | |
CN109164033B (zh) | 一种基于电阻率法的获取非饱和土壤水力传导力的方法 | |
CN104062654B (zh) | 一种基于超宽带雷达的土壤含水量测量方法 | |
CN107621438B (zh) | 一种坡面地形演变与水蚀过程耦合的动态监测方法 | |
CN106970036B (zh) | 一种绝缘子表面污秽度的判定方法 | |
CN114397427B (zh) | 一种基于环刀土样上渗过程的土壤水力学特性预测方法 | |
CN115754221A (zh) | 一种步进式时域反射土壤有效态养分原位速测系统及方法 | |
CN112666554A (zh) | 一种沥青路面雷达振幅特征裂缝宽度识别方法 | |
CN210243553U (zh) | 一种环刀称重法土壤水分传感器标定装置 | |
CN105547957A (zh) | 土壤下渗湿润峰吸力和降雨初损量计算方法 | |
CN114755269A (zh) | 基于无损时域反射技术的黄土湿陷性原位评价方法及系统 | |
CN113468804B (zh) | 一种基于矩阵束和深度神经网络的地下管道识别方法 | |
Irmak et al. | Performance of frequency-domain reflectometer, capacitance, and psuedo-transit time-based soil water content probes in four coarse-textured soils | |
CN115879289A (zh) | 一种基于水分下渗过程原位快速获取土壤水力学参数的方法 | |
Comegna et al. | Nonreactive solute transport in variously structured soil materials as determined by laboratory-based time domain reflectometry (TDR) | |
Baker et al. | Measuring water exchange between soil and atmosphere with TDR-microlysimetry | |
CN102692359A (zh) | 一种土壤水气传导率原位测量系统与测量方法 | |
CN110836799A (zh) | 一种吸力控制型快速非饱和土试样制样控制系统及方法 | |
Shi et al. | Design and performance analysis of soil temperature and humidity sensor | |
Galagedara et al. | Measuring and modeling of direct ground wave depth penetration under transient soil moisture conditions | |
Kanemitsu et al. | Quantitative evaluation of the relative permittivity of artificial soil with altered soil types and water content | |
Tanriverdi | Using TDR in the agricultural water management |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |