CN114111553B - 一种快速获取新增耕地重构土体填充层厚度的方法 - Google Patents
一种快速获取新增耕地重构土体填充层厚度的方法 Download PDFInfo
- Publication number
- CN114111553B CN114111553B CN202111413763.5A CN202111413763A CN114111553B CN 114111553 B CN114111553 B CN 114111553B CN 202111413763 A CN202111413763 A CN 202111413763A CN 114111553 B CN114111553 B CN 114111553B
- Authority
- CN
- China
- Prior art keywords
- filling layer
- thickness
- dielectric constant
- time
- radar
- 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.)
- Active
Links
- 239000002689 soil Substances 0.000 title claims abstract description 60
- 238000000034 method Methods 0.000 title claims abstract description 47
- 230000000149 penetrating effect Effects 0.000 claims abstract description 36
- 238000005259 measurement Methods 0.000 claims abstract description 32
- 238000004364 calculation method Methods 0.000 claims abstract description 18
- 238000002474 experimental method Methods 0.000 claims abstract description 8
- 230000000694 effects Effects 0.000 claims abstract description 5
- 239000000945 filler Substances 0.000 claims abstract description 5
- 238000013401 experimental design Methods 0.000 claims abstract description 4
- 238000007781 pre-processing Methods 0.000 claims abstract description 4
- 238000001514 detection method Methods 0.000 claims description 14
- 238000001914 filtration Methods 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 10
- 238000004458 analytical method Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 4
- 230000006698 induction Effects 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 2
- 238000013461 design Methods 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims description 2
- 239000010410 layer Substances 0.000 description 52
- 230000008569 process Effects 0.000 description 6
- 239000002344 surface layer Substances 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000005056 compaction Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000004062 sedimentation Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000009412 basement excavation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B7/00—Measuring arrangements characterised by the use of electric or magnetic techniques
- G01B7/02—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring length, width or thickness
- G01B7/06—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring length, width or thickness for measuring thickness
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/411—Identification of targets based on measurements of radar reflectivity
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Geophysics And Detection Of Objects (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及土体填充层厚度的测量领域，具体的是一种快速获取新增耕地重构土体填充层厚度的方法，包括以下步骤：布置实验，根据不同填充位置进行实验设计，并对填充层上下分界面位置进行标定、获取每个填充位置下地探地雷达实时的剖面数据、探测图像、对所获取的雷达数据进行预处理，获取效果明确的雷达剖面数据、根据雷达图像获取填充物所在的上、下分界面信号，根据分界面时间计算探地雷达在填充层的双程走时，根据标定深度计算填充层介电常数、结合点测数据对标定的填充层位置的介电常数进行表示，获取精确的介电常数。本发明可以直接对地下填充层厚度进行计算，不再需要从表土层开始进而类推的方式，使得深层土层厚度计算更加简便。
Description
技术领域
本发明涉及土体填充层厚度的测量领域，具体的是一种快速获取新增耕地重构土体填充层厚度的方法。
背景技术
土壤重构是土地复垦中的重要环节，在土壤复垦过程中，填充物的选择、填充厚度、覆土厚度的确定均会对土壤质量造成影响，而在复垦过程中，填充层的厚度或大或小，导致覆土层的厚度差异较大，从而造成不同地方土壤结构不同，而对于作物生长过程中，土壤结构不同带来的问题便是作物长势差异较大。而在复垦过程中，同一地区一般采用同一种填充物进行土地复垦，因此主要影响作物生长的因素便是覆土厚度及填充厚度的差异，确定好覆土厚度和填充厚度的比例，不仅可以减小土壤结构对作物生长的影响，而且可以保证经济效益的稳定。
探地雷达(GPR)作为一种新型的浅层地球物理探测技术，已广泛应用于地理探测、工程检测等领域。探地雷达具有快读、便捷、无损、连续、可重复探测、可大面积测量等特征。在土层厚度及分层探测方面，探地雷达表现出了一定的优势，可以应用在土地复垦工程中对土层厚度的验收。
在传统的探地雷达分层识别中，主要是根据人为经验判断和波形简化方式来确定分层位置，但这样做存在较大误差，对于表层土壤的探测可能影响不大，但对之后层次的土壤厚度计算精度将会减小，所以，忽略表层直接对表层之下的土层进行计算，在一定程度上可以提高深层土壤的探测精度，同时，对于只需要表层之下土层厚度的项目，可以在一定程度上减少工作量。
发明内容
为解决上述背景技术中提到的不足，本发明的目的在于提供一种快速获取新增耕地重构土体填充层厚度的方法。
本发明的目的可以通过以下技术方案实现：
一种快速获取新增耕地重构土体填充层厚度的方法，所述方法包括以下步骤：
一、布置实验，根据不同填充位置进行实验设计，并对填充层上下分界面位置进行标定；
二、获取每个填充位置下地探地雷达实时的剖面数据、探测图像；
三、对所获取的雷达数据进行预处理，获取效果明确的雷达剖面数据；
四、根据雷达图像获取填充物所在的上、下分界面信号，结合时频分析方法(短时傅里叶变换)，根据分界面时间计算探地雷达在填充层的双程走时，根据标定深度计算填充层介电常数；
五、结合点测数据对标定的填充层位置的介电常数进行表示，获取精确的介电常数；
六、基于获取的介电常数来计算填充层厚度的实时变化情况。
进一步地，所述步骤一在实验时分别设计填充层厚度相同、覆土厚度不同的实验区域，并且每个区域均标定填充层上下分界面的位置及其深度。
进一步地，所述步骤二的具体步骤如下：
1)在每个覆土区域进行位置标定，确定标定的位置和上下标定位置之间的厚度；
2)设置探地雷达参数，在不同覆土区域进行探地雷达线测数据采集，获取实时线测剖面图；
3)保持探地雷达采集参数不变，在标定位置进行点测数据采集。
进一步地，所述步骤三的要操作流程为：去直流漂移、AGG增益、背景去除、巴斯沃斯滤波，具体步骤如下：
1)对采集的雷达数据进行去除直流漂移处理，其目的是使有效信号不受漂移现象的影响；
2)去直达波处理是去除雷达数据中由于收发天线距离较近以及地面波与空气波产生的低频信号，从而消除感应现象的失真，可以在保留高频信号的同时移除不需要的低频信号；
3)背景去除处理，由于带通滤波整体上改善了信号的信噪比，但在某些特定情况下，雷达回波剖面信号中包含呈水平线的“背景”干扰，其掩埋了有效反射信号，对信号进行背景去除处理；
4)滤波处理采用方法为巴斯沃斯滤波，获取分层更加明显的探地雷达图像。
进一步地，所述步骤四的具体步骤如下：
1)根据标定的深度，结合短时傅里叶变换，获取上下分界面的时间信息，计算填充层的双程走时，具体使用以下公式计算：
t＝t1-t2，t表示填充层的双程走时，t1表示电磁波到达上分界面的时间，t2代表电磁波到达下分界面的时间；
2)根据标定深度，结合计算的时间，获取填充层的介电常数，以此来表征其余区域的填充层介电常数，介电常数获取方法使用以下公式计算：
ε＝(ct/h)2；其中，ε代表介电常数，c代表光速，一般为3.0×108m/s，h代表标定深度。
进一步地，所述步骤五中获取的介电常数以点测数据获取的介电常数为准。
进一步地，所述步骤六中进行多次点测数据采集，采用空间插值获取测线上个点的介电常数，最后进行填充层厚度的反算。
本发明的有益效果：
本发明可以直接对地下填充层厚度进行计算，越过了表层土壤，不再需要从表土层开始进而类推的方式，使得深层土层厚度计算更加简便，并且该方法适用于重构土壤及分层较为简单的土体结构，在一定程度上避免了误差的传递，使得填充层厚度的探测精度更高。
本发明未考虑表层土的影响，直接从分层位置进行计算，由于主要是依据分层位置进行时间计算，避免了由表层土壤计算的分层不清的问题。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图；
图1是本发明基于探地雷达获取重构土壤填充层厚度的示意图；
图2是本发明探地雷达图像处理流程示意图；
图3是本发明探地雷达采集到的原始图像示意图；
图4是本发明预处理之后的探地雷达图像示意图；
图5是本发明短时傅里叶方法处理后的探地雷达图像示意图；
图6是本发明填充层厚度变化曲线示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
一种快速获取新增耕地重构土体填充层厚度的方法，如图1-6所示，包括以下步骤：
一、布置实验，根据不同填充位置进行实验设计，并对填充层上下分界面位置进行标定；
其中，分别设计填充层厚度相同、覆土厚度不同的实验区域，并且每个区域均标定填充层上下分界面的位置及其深度；
由于探地雷达主要应用于表层土壤结构探测，本发明主要研究表层土壤结构进行探测，研究深度均在1m以内，即复垦土壤的开垦范围。
具体地，本发明中实验区域共有三个，通过开挖三个50×50cm的正方形坑，每个坑的开挖深度不同，分别为60cm、80cm、100cm，其中，填充厚度均为20cm，覆土厚度为40cm、60cm、80cm，填充物均为建筑垃圾，实验中使用60cm的地方作为具体分析；
然后利用探地雷达在表层进行探测，采样点设计为1026，采样间隔为0.02m。由于使用机械压实会影响土层的分层变化，所以主要通过土壤自然沉降，对土壤沉降稳定后的土层进行探测，从而减小人为压实产生的影响。
二、获取每个填充位置下探地雷达实时的剖面数据、探测图像；
其中，在图像获取过程中，需要点测数据与线测数据同时采集，在标定位置不仅获取具体线测数据，同时采集点测数据；
具体步骤如下：
1)在每个覆土区域进行位置标定，确定标定的位置和上下标定位置之间的厚度；
2)设置探地雷达参数，在不同覆土区域进行探地雷达线测数据采集，获取实时线测剖面图；
3)保持探地雷达采集参数不变，在标定位置进行点测数据采集。
具体的，分别对三个不同覆土区域进行探地雷达探测，获取每个区域的雷达剖面图。测线长为0.5m，共获取三条测线的剖面图，在标定位置使用探地雷达进行点测，具体操作方式为：保持探地雷达不动，连续向同一位置发射20条电磁波，以20条电磁波的叠加作为测点的图像，每个覆土区域均获取一个标定位置的点测数据和线测数据。
三、对所获取的雷达数据进行预处理，获取效果明确的雷达剖面数据；
对探地雷达获取的初始图像进行预处理，使分层效果更加明显，主要操作流程为：去直流漂移、AGG增益、背景去除、巴斯沃斯滤波，具体如下：
1)首先对采集的雷达数据进行去除直流漂移处理，其目的是使有效信号不受漂移现象的影响；
2)去直达波处理，主要是去除雷达数据中由于收发天线距离较近以及地面波与空气波产生的低频信号，从而消除感应现象的失真，它可以在保留高频信号的同时，移除不需要的低频信号；
3)之后进行背景去除处理，由于带通滤波整体上改善了信号的信噪比，但在某些特定情况下，雷达回波剖面信号中包含呈水平线的“背景”干扰，这些掩埋了有效反射信号，因此对信号进行背景去除处理。
4)最后进行滤波处理，主要采用方法为巴斯沃斯滤波，获取分层更加明显的探地雷达图像。
四、根据雷达图像获取填充物所在的上、下分界面信号，结合时频分析方法(短时傅里叶变换)，根据分界面时间计算探地雷达在填充层的双程走时，根据标定深度计算填充层介电常数，具体地步骤如下：
1)根据标定的深度，通过短时傅里叶方法对波形进行变换，获取上下分界面的时间信息，计算填充层的双程走时，具体使用以下公式计算：
t＝t1-t2，t表示填充层的双程走时，t1表示电磁波到达上分界面的时间，t2代表电磁波到达下分界面的时间；
2)根据标定深度，结合计算的时间，获取填充层的介电常数，以此来表征其余区域的填充层介电常数，介电常数获取方法使用以下公式计算：
ε＝(ct/h)2；其中，ε代表介电常数，c代表光速，一般为3.0×108m/s，h代表标定深度。
五、结合点测数据对标定的填充层位置的介电常数进行表示，获取较为精确的介电常数；
1)根据所得的点测数据下的介电常数作为最后的介电常数；
2)根据校正参数，对未标定位置分层位置进行双程走时计算，根据通过插值方法获取的预测的介电常数作为填充层的介电常数。
六、基于获取的介电常数来计算填充层厚度的实时变化情况。
进行多次点测数据采集，采用空间插值获取测线上个点的介电常数，最后进行填充层厚度的反算，如图6所示。
由于覆土厚度不同，在不同区域预测的介电常数有所差别，主要是由于地下结构不同、含水率的影响不同，本发明进行介电常数预测，不涉及含水率的影响。
由于本发明未设计实验，填充层厚度较为均一，同时位置大概均处于同一位置，在探地雷达图像中易于分辨，但在实际中，由于土壤结构的变化情况，可能无法精确判断分层位置，需利用单道波形进行分析，同时可以结合具体的时频分析方法进行时间判定。由于在深度计算中无法进行逐点计算，所以需要借助插值手段进行图层变化预测。
综上，本发明可以直接对地下填充层厚度进行计算，越过了表层土壤，不再需要从表土层开始进而类推的方式，使得深层土层厚度计算更加简便，并且该方法适用于重构土壤及分层较为简单的土体结构，在一定程度上避免了误差的传递，使得填充层厚度的探测精度更高。
本发明未考虑表层土的影响，直接从分层位置进行计算，由于主要是依据分层位置进行时间计算，避免了由表层土壤计算的分层不清的问题。
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。
Claims (6)
1.一种快速获取新增耕地重构土体填充层厚度的方法，其特征在于，所述方法包括以下步骤：
一、布置实验，根据不同填充位置进行实验设计，并对填充层上下分界面位置进行标定；
二、获取每个填充位置下探地雷达实时的剖面数据、探测图像；
三、对所获取的雷达数据进行预处理，获取效果明确的雷达剖面数据；
四、根据雷达图像获取填充物所在的上、下分界面信号，结合时频分析方法，根据分界面时间计算探地雷达在填充层的双程走时，根据标定深度计算填充层介电常数；
五、结合点测数据对标定的填充层位置的介电常数进行表示，获取精确的介电常数；
六、基于获取的介电常数来计算填充层厚度的实时变化情况；
所述步骤四的具体步骤如下：
1)根据标定的深度，结合短时傅里叶变换，获取上下分界面的时间信息，计算填充层的双程走时，具体使用以下公式计算：
t＝t1-t2，t表示填充层的双程走时，t1表示电磁波到达上分界面的时间，t2代表电磁波到达下分界面的时间；
2)根据标定深度，结合计算的时间，获取填充层的介电常数，以此来表征其余区域的填充层介电常数，介电常数获取方法使用以下公式计算：
ε＝(ct/h)2；其中，ε代表介电常数，c代表光速，为3.0×108m/s，h代表标定深度；
所述步骤五的具体步骤如下：
1)根据所得的点测数据下的介电常数作为最后的介电常数；
2)根据校正参数，对未标定位置分层位置进行双程走时计算，根据通过插值方法获取的预测的介电常数作为填充层的介电常数。
2.根据权利要求1所述的一种快速获取新增耕地重构土体填充层厚度的方法，其特征在于，所述步骤一在实验时分别设计填充层厚度相同、覆土厚度不同的实验区域，并且每个区域均标定填充层上下分界面的位置及其深度。
3.根据权利要求1所述的一种快速获取新增耕地重构土体填充层厚度的方法，其特征在于，所述步骤二的具体步骤如下：
1)在每个覆土区域进行位置标定，确定标定的位置和上下标定位置之间的厚度；
2)设置探地雷达参数，在不同覆土区域进行探地雷达线测数据采集，获取实时线测剖面图；
3)保持探地雷达采集参数不变，在标定位置进行点测数据采集。
4.根据权利要求1所述的一种快速获取新增耕地重构土体填充层厚度的方法，其特征在于，所述步骤三的要操作流程为：去直流漂移、AGC增益、背景去除、巴斯沃斯滤波，具体步骤如下：
1)对采集的雷达数据进行去除直流漂移处理，其目的是使有效信号不受漂移现象的影响；
2)去直达波处理是去除雷达数据中由于收发天线距离较近以及地面波与空气波产生的低频信号，从而消除感应现象的失真，可以在保留高频信号的同时移除不需要的低频信号；
3)背景去除处理，由于带通滤波整体上改善了信号的信噪比，但在某些特定情况下，雷达回波剖面信号中包含呈水平线的背景干扰，其掩埋了有效反射信号，对信号进行背景去除处理；
4)滤波处理采用方法为巴斯沃斯滤波，获取分层更加明显的探地雷达图像。
5.根据权利要求1所述的一种快速获取新增耕地重构土体填充层厚度的方法，其特征在于，所述步骤五中获取的介电常数以点测数据获取的介电常数为准。
6.根据权利要求1所述的一种快速获取新增耕地重构土体填充层厚度的方法，其特征在于，所述步骤六中进行多次点测数据采集，采用空间插值获取预测的介电常数，最后进行填充层厚度的反算。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111413763.5A CN114111553B (zh) | 2021-11-25 | 2021-11-25 | 一种快速获取新增耕地重构土体填充层厚度的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111413763.5A CN114111553B (zh) | 2021-11-25 | 2021-11-25 | 一种快速获取新增耕地重构土体填充层厚度的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114111553A CN114111553A (zh) | 2022-03-01 |
| CN114111553B true CN114111553B (zh) | 2023-12-08 |
Family
ID=80373162
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111413763.5A Active CN114111553B (zh) | 2021-11-25 | 2021-11-25 | 一种快速获取新增耕地重构土体填充层厚度的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114111553B (zh) |
Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6380747B1 (en) * | 1998-05-12 | 2002-04-30 | Jentek Sensors, Inc. | Methods for processing, optimization, calibration and display of measured dielectrometry signals using property estimation grids |
| AU2002318970B2 (en) * | 2001-07-06 | 2007-03-29 | Gecoz Pty Ltd | Method of determining salinity of an area of soil |
| US8947102B1 (en) * | 2011-08-05 | 2015-02-03 | The United States Of America As Represented By The Secretary Of Agriculture | Soil water and conductivity sensing system |
| CN106767371A (zh) * | 2016-12-27 | 2017-05-31 | 山东农业大学 | 一种探地雷达快速探测复垦土壤覆盖表土厚度的方法 |
| CN107305226A (zh) * | 2017-04-28 | 2017-10-31 | 厦门大学 | 一种层状介质介电常数和厚度同时反演算法 |
| CN108982548A (zh) * | 2018-07-20 | 2018-12-11 | 浙江大学 | 一种基于被动微波遥感数据的地表土壤水分反演方法 |
| CN109031290A (zh) * | 2018-07-27 | 2018-12-18 | 山东农业大学 | 一种土地工程的质量检测方法 |
| CN109115104A (zh) * | 2018-10-23 | 2019-01-01 | 西北农林科技大学 | 一种基于磁化率测量土壤层厚度的测定方法 |
| US10234552B1 (en) * | 2018-06-27 | 2019-03-19 | University Of South Florida | Precise infrastructure mapping using full-waveform inversion of ground penetrating radar signals |
| CN110161495A (zh) * | 2019-05-23 | 2019-08-23 | 中国地质大学（北京） | 一种土壤有效土层厚度确定方法及系统 |
| CN110579760A (zh) * | 2019-09-23 | 2019-12-17 | 中国地质大学（北京） | 一种矿区土壤物理性质的三维呈现方法及系统 |
| CN111880176A (zh) * | 2020-07-24 | 2020-11-03 | 安徽理工大学 | 一种煤矸石充填复垦地土壤含水量的测定方法 |
| CN111929678A (zh) * | 2020-07-24 | 2020-11-13 | 安徽理工大学 | 一种浅地表土壤含水量的测定方法 |
| CN113156427A (zh) * | 2021-05-08 | 2021-07-23 | 吉林大学 | 探地雷达数据的反演方法和装置 |
| CN113156423A (zh) * | 2021-04-08 | 2021-07-23 | 淮北矿业股份有限公司 | 一种基于探测雷达的复垦土壤耕层含水率预测方法 |
| CN113391283A (zh) * | 2021-05-24 | 2021-09-14 | 中国农业大学 | 基于探地雷达的土壤分层信息识别方法和装置 |
-
2021
- 2021-11-25 CN CN202111413763.5A patent/CN114111553B/zh active Active
Patent Citations (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6380747B1 (en) * | 1998-05-12 | 2002-04-30 | Jentek Sensors, Inc. | Methods for processing, optimization, calibration and display of measured dielectrometry signals using property estimation grids |
| AU2002318970B2 (en) * | 2001-07-06 | 2007-03-29 | Gecoz Pty Ltd | Method of determining salinity of an area of soil |
| US8947102B1 (en) * | 2011-08-05 | 2015-02-03 | The United States Of America As Represented By The Secretary Of Agriculture | Soil water and conductivity sensing system |
| CN106767371A (zh) * | 2016-12-27 | 2017-05-31 | 山东农业大学 | 一种探地雷达快速探测复垦土壤覆盖表土厚度的方法 |
| CN107305226A (zh) * | 2017-04-28 | 2017-10-31 | 厦门大学 | 一种层状介质介电常数和厚度同时反演算法 |
| US10234552B1 (en) * | 2018-06-27 | 2019-03-19 | University Of South Florida | Precise infrastructure mapping using full-waveform inversion of ground penetrating radar signals |
| CN108982548A (zh) * | 2018-07-20 | 2018-12-11 | 浙江大学 | 一种基于被动微波遥感数据的地表土壤水分反演方法 |
| CN109031290A (zh) * | 2018-07-27 | 2018-12-18 | 山东农业大学 | 一种土地工程的质量检测方法 |
| CN109115104A (zh) * | 2018-10-23 | 2019-01-01 | 西北农林科技大学 | 一种基于磁化率测量土壤层厚度的测定方法 |
| WO2020082475A1 (zh) * | 2018-10-23 | 2020-04-30 | 西北农林科技大学 | 一种基于磁化率测量土壤层厚度的测定方法 |
| CN110161495A (zh) * | 2019-05-23 | 2019-08-23 | 中国地质大学（北京） | 一种土壤有效土层厚度确定方法及系统 |
| CN110579760A (zh) * | 2019-09-23 | 2019-12-17 | 中国地质大学（北京） | 一种矿区土壤物理性质的三维呈现方法及系统 |
| CN111880176A (zh) * | 2020-07-24 | 2020-11-03 | 安徽理工大学 | 一种煤矸石充填复垦地土壤含水量的测定方法 |
| CN111929678A (zh) * | 2020-07-24 | 2020-11-13 | 安徽理工大学 | 一种浅地表土壤含水量的测定方法 |
| CN113156423A (zh) * | 2021-04-08 | 2021-07-23 | 淮北矿业股份有限公司 | 一种基于探测雷达的复垦土壤耕层含水率预测方法 |
| CN113156427A (zh) * | 2021-05-08 | 2021-07-23 | 吉林大学 | 探地雷达数据的反演方法和装置 |
| CN113391283A (zh) * | 2021-05-24 | 2021-09-14 | 中国农业大学 | 基于探地雷达的土壤分层信息识别方法和装置 |
Non-Patent Citations (3)
| Title |
|---|
| A novel multi-frequency inversion algorithm for the retrieval of the subsurface properties of layered soil media;Chih-hao Kuo;《2007 IEEE Antennas and Propagation Society International Symposium》;第1785-1788页 * |
| 基于探地雷达的特厚煤层厚度动态探测技术;刘万里;《煤炭学报》;第46卷(第8期);第2706-2714页 * |
| 煤矿采空区复垦土地探地雷达评价技术研究与应用;安志元;《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》(第02(2018)期);第A011-922页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114111553A (zh) | 2022-03-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Cui et al. | Modeling tree root diameter and biomass by ground-penetrating radar | |
| CN101915771B (zh) | 基于低频探地雷达地波法的测定砂质土壤含水量的方法 | |
| CN107576674A (zh) | 一种基于探地雷达测量土壤压实程度的方法 | |
| AU2011100086A4 (en) | A method for measuring biomass of tree-roots by ground-penetrating radar | |
| CN107479042B (zh) | 一种表层岩溶带空间蓄水能力的估算方法 | |
| CN113325419B (zh) | 一种路面坍塌检测方法 | |
| Rosa et al. | Determining the number of manual measurements required to improve peat thickness estimations by ground penetrating radar | |
| CN112666554A (zh) | 一种沥青路面雷达振幅特征裂缝宽度识别方法 | |
| CN111880176A (zh) | 一种煤矸石充填复垦地土壤含水量的测定方法 | |
| Pezdir et al. | Investigating peatland stratigraphy and development of the Šijec bog (Slovenia) using near-surface geophysical methods | |
| Samet et al. | Using interpolation techniques to determine the optimal profile interval in ground-penetrating radar applications | |
| CN108897050B (zh) | 一种可控源音频大地电磁测深法最佳收发距确定方法及系统 | |
| CN109541695B (zh) | 人工场源频率域电场梯度远区视电阻率快速成像方法 | |
| CN114111553B (zh) | 一种快速获取新增耕地重构土体填充层厚度的方法 | |
| CN111142165A (zh) | 一种利用探地雷达获取含水层的水位信息的方法 | |
| Ruan et al. | Ground Penetrating Radar (GPR) Identification Method for Agricultural Soil Stratification in a Typical Mollisols Area of Northeast China | |
| CN113688692A (zh) | 基于探地雷达时频尺度特征的供水管道漏损检测方法 | |
| CN111929678A (zh) | 一种浅地表土壤含水量的测定方法 | |
| Cousin et al. | From spatial-continuous electrical resistivity measurements to the soil hydraulic functioning at the field scale | |
| CN113391283A (zh) | 基于探地雷达的土壤分层信息识别方法和装置 | |
| CN113156423A (zh) | 一种基于探测雷达的复垦土壤耕层含水率预测方法 | |
| CN115638719A (zh) | 一种快速获取新增耕地重构土地覆土层厚度的方法 | |
| CN109031266A (zh) | 一种土壤湿润体体积的估算方法 | |
| CN112666552B (zh) | 一种探地雷达数据背景杂波自适应去除方法 | |
| CN113030953A (zh) | 基于探地雷达与波前干扰偏移的树干内部缺陷成像方法 |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |