欢迎您访问《智慧农业(中英文)》官方网站! English

Smart Agriculture ›› 2022, Vol. 4 ›› Issue (1): 84-96.doi: 10.12133/j.smartag.SA202202010

• 专题--作物生长及其环境监测 • 上一篇    下一篇

基于探地雷达的田块尺度下不同深度土壤含水量监测

张文瀚1(), 杜克明2(), 孙彦坤1(), 刘布春2, 孙忠富2, 马浚诚2, 郑飞翔2   

  1. 1.东北农业大学 资源与环境学院,黑龙江 哈尔滨 150030
    2.中国农业科学院农业环境与可持续发展研究所,北京 100081
  • 收稿日期:2021-12-25 出版日期:2022-03-30
  • 基金资助:
    国家重点研发计划课题(2016YFD0300606);国家自然科学基金项目(31628015)
  • 作者简介:张文瀚(1997-),男,硕士研究生,研究方向为农业资源与环境监测。E-mail:2542312177@qq.com
  • 通信作者:

Monitoring Specified Depth Soil Moisture in Field Scale with Ground Penetrating Radar

ZHANG Wenhan1(), DU Keming2(), SUN Yankun1(), LIU Buchun2, SUN Zhongfu2, MA Juncheng2, ZHENG Feixiang2   

  1. 1.College of Resources and Environment, Northeast Agricultural University, Harbin 150030
    2.Institute of Environment and Sustainable Development in Agriculture, Chinese Academy of Agricultural Sciences, Beijing 100081
  • Received:2021-12-25 Online:2022-03-30

摘要:

为确定田块尺度下探地雷达对不同深度及相邻反射层间土壤含水量的反演精度、有效反演深度、最佳反演深度及最优反演模型,本研究采用1000 MHz中心频率探地雷达设备,分别在无降雨偏干旱土壤和降雨后湿润土壤两种条件下,在选定农田区域基于共中心点法采集雷达波数据,提取有效地表波与反射波数据,通过双曲线拟合法分别获取不同深度反射层雷达波的传播速度,计算得出土壤的相对介电常数,最后根据土壤体积含水量和介电常数之间的经验模型计算获得不同深度的土壤体积含水量。通过Topp、Roth、Herkelrath和Ferre四种经验模型分别进行土壤体积含水量反演测定,同时以利用烘干法获取的代表性测定土壤含水量实测值为指标进行精度验证。田间试验结果表明,1000 MHz探地雷达的有效反演深度范围为0~50 cm;土壤偏干旱和偏湿润条件的最佳反演深度分别为50 cm和40 cm;Roth模型相关性最好,决定系数R2最高为0.750,且Roth模型反演土壤含水量值最稳定,在土壤偏干旱和偏湿润条件下均方根误差(Root Mean Square Error,RMSE)平均值分别为0.0401 m3/m3和0.0335 m3/m3;相对误差(Relative Error,RE)最低为3.0%。探地雷达具备对定量深度田间土壤含水量快速、精准监测的能力,但其反演模型需根据不同土壤类型等条件进行相应参数校正。

关键词: 探地雷达, 土壤含水量监测, 共中心点法, 最佳反演深度, 最优反演模型, Roth模型

Abstract:

Ground-penetrating radar (GPR) is one of the emerging technologies for soil moisture measurement. However, the measurement accuracy is difficult to determine due to some influence factors including radar wave frequency, soil texture type, etc. The GPR equipment with 1000 MHz center frequency and the measurement method of common midpoint (CMP) were adopted in the research to collect radar wave raw data in the selected field area under arid soil and moist soil conditions. The transmitter and receiver antennas of the GPR equipment were moved 0.01 m respectively in opposite directions on each radar wave raw data collection. Therefore, a CMP radar image consisted of 100 pieces of radar wave raw data by increasing the antenna distance from 0 m to 2 m. Each radar wave raw data indicated that the radar waves were reflected in the reflective layer with different dielectric constant under the same antenna distance. And the reflected and refracted radar waves were acquired by the receiving antenna at different two-way travel time respectively, and recorded in the computer. The collection of CMP soundings aimed to determine the inversion accuracy, optimum inversion depth, effective inversion depth and optimal inversion model of soil moisture content at different depth ranges and adjacent reflective layers by GPR at field scale. The reflected and refracted radar wave data were extracted from the raw data. The velocities of the surface waves and reflected waves were obtained respectively from the line slope of the surface wave data and the hyperbolic curves fitting of the reflected wave data. In addition, the relative dielectric constant of the soil at specified depth were deduced according to the soil dielectric constant and its reflected wave velocity. Moreover, 4 different models including Topp, Roth, Herkelrath and Ferre were used to figure out the soil volumetric water content inversion. Meanwhile, the measured data of soil volumetric moisture content obtained by oven drying method were used to verify the accuracy of the inversion results. The results showed that the effective inversion depth of 1000 MHz GPR ranged from 0 to 50 cm. The best inversion depth was 50 cm in arid soil and 40 cm in moist soil. The Roth model had the best correlation and stability with the highest R2 was 0.750, the Root Mean Square Error (RMSE) was 0.0114 m3/m3 and the lowest Relative Error (RE) was 3.0%. The GPR could possess the capacity of quick, precise and non-destructive measurement of specified depth soil moisture in field scale. The inversion model of soil moisture content needs to be calibrated according to different soil conditions.

Key words: ground penetrating radar, soil water content monitoring, common midpoint, optimum inversion depth, optimum inversion model, Roth model

中图分类号: