摘要
为快速准确获取土壤含水率信息,便于农业精准灌溉,引入支持向量机算法(SVM)对4种不同干湿交替处理下超声波速度与土壤含水率进行拟合分析和回归训练优化,构建基于超声波速度的土壤含水率预测模型。结果显示,与传统的烘干法相比较,利用该模型在田间验证土壤含水率,平均相对误差为1.5%左右。研究结果表明,基于SVM模型构建的超声波速度-土壤含水率预测模型能够较好地描述被研究区域内土壤含水率,可为利用超声波特性实现对农田土壤水分的持续监测提供参考。
土壤水分的有效性对植被的光合及蒸腾作用起着关键作用,土壤水分不仅是作物赖以生长的重要因素之
土壤的基本物理参数如颗粒尺寸、孔隙度、含水量等,在降雨入渗、土体水分蒸散等气象条件或灌溉排水等水文循环过程中都会发生相应变化,直接导致土壤微结构的改
本研究考虑水分干湿交替对于土壤结构的影响,对4种不同程度干湿交替处理下的土壤含水率与超声波速度进行相关性试验,并采用支持向量机(SVM)的机器学习算法构建模型并选取优化参数值,对比预测集的结果精确度,旨在引入机器学习算法实现声学监测农田土壤含水率信息。
供试土样采自华南农业大学科研试验基地(23°09´48"N,113°21´44"E)常年耕作土壤耕层0~15 cm,采用直径为5 cm的取土圆柱环刀对被研究土壤区域进行取样,测量土样样本原始条件包括土壤水分以及组成成分等。取样过程中需要对表面的石块和落叶杂物等进行清理,取样之后为了保证土壤的结构不被破坏,不对土壤进行过筛研磨,以确保土壤可以更加接近自然状态下的条件,按照林业行业标准LY/T 1225—1999(土壤颗粒组成)测量被研究土壤的颗粒含量,分别为:砂粒(颗粒直径为0.05~2.00 mm)占48.7%、粉粒(颗粒直径为0.002~0.050 mm)占39.1%、黏粒(颗粒直径<0.002 mm)占12.2%。
对于一般田间耕作层土壤,土壤最大含水率为30%、土壤平均含水率为15%、土壤凋萎含水率为5
如

图1 土壤样本实物图
Fig.1 Physical picture of soil sample
(1) |
干湿交替设计:土壤不断经历水分变化,土粒细化团聚体重组,将对声学检测土壤含水率造成影响,因此需要对土壤干湿交替的变化进行考虑。分别将土壤样本自然风干,期间多次取样利用烘干法确定土壤含水率,直到水分达到终止含水率时,判断为完成1次干湿交替。接着重复上述操作,直到共完成4次干湿交替,视为结束。
在土壤由湿到干的变化过程中,将利用如

图2 超声波检测装置示意图
Fig.2 Schematic diagram of ultrasonic testing device
1.力值显示控制仪Force value display control instrument;2.非金属超声波检测仪Nonmetallic ultrasonic detector;3.支架Support;4.土样Soil sample;5.换能器安装盒Transducer mounting box;6.滑轨Slide rail;7.压力传感器Pressure sensor;8.超声换能器连接线Ultrasonic transducer cable;9.步进电机驱动器Stepper motor driver;10.电源变压器Power transformer;11.步进电机控制器Stepper motor controller.
将非金属超声检测仪的发射与接收换能器安装在换能器的安装盒中,调整到合适的高度正好与土样容器两侧的圆孔对正,接着涂抹超声波耦合剂在土壤表面使换能器可以紧贴土壤两侧。换能器安装盒旁边是压力传感器,首先设置压力传感器的最大限度为20 N,接着启动步进电机驱动装置的丝杆移动,带动着压力传感器推动换能器压紧土壤,直到此时的压力值已达设置值为止,这里步进电机由电源变压器驱动。记录此时土壤样本的有效长度d,接着启动非金属超声波检测仪,分别通过发射和接收换能器来完成超声波检测,重复进行10次检测,记录通过土壤样本的声速值。
试验数据采取了SVM算法模型对数据集进行训练并得到模型优化值,再利用测试集检验模型预测土壤含水率的精确度。
1)试验数据归一化处理。为避免因输入、输出数据的数量级差别较大导致训练过程不稳定以及预测误差较大等情况。首先对插值后的超声波速度与土壤含水率数据进行归一化处
(2) |
2)SVM算法样本集构造。所有归一化的数据,选取土壤样本由湿-干过程中超声波速度值作为输入数据,对应的土壤含水率作为输出数据。初始土壤含水率为15%、20%、25%和30%的土壤样本编号分别记为T1~T4,每个样本数据集总数分别有108个。按25︰2的比例划分为训练集和测试集,即分别随机选取100个超声波速度数据作为训练样本集,8个超声波速度数据作为测试样本集。
3)SVM模型参数选择及优化。在模型建立的过程中,SVM选择适用性比较强的高斯核函数(RBF)为核函数,形式如
(3) |
编号 Number | 惩罚系数C Penalty coefficient | 内核参数g Kernelparameters | 损失函数P Loss function | 核函数类型参数 Kernel function type arguments |
---|---|---|---|---|
T1 | 1.41 | 0.35 | 0.01 | 0.82 |
T2 | 45.25 | 11.31 | 0.01 | 0.54 |
T3 | 2.83 | 2.83 | 0.01 | 0.71 |
T4 | 64.00 | 0.50 | 0.01 | 0.88 |
4)数据分析方法及算法评价指标。所测声学和土壤含水率数据采用Matlab 2019软件进行整理统计。使用SPSS 24.0的皮尔逊相关系数算法计算超声波速度和土壤含水率的相关系数,并进行显著性检验。利用SVM算法得到土壤含水率的预测结果之后,选择能够反映出测量精密度的平均均方误差(MRE)以及能够衡量回归模型整体拟合度的决定系数
超声波速度对不同初始土壤含水率组的敏感程度分析结果显示,T1组的超声波速度和土壤含水率极显著负相关(P<0.001,r=-0.985),T2组的超声波速度和土壤含水率极显著负相关(P<0.001,r=-0.992),T3组的超声波速度和土壤含水率极显著负相关(P<0.001,r=-0.959),T4组的超声波速度和土壤含水率极显著负相关(P<0.001,r=-0.962)。
初始含水率分别为15%、20%、25%、30%的干湿交替处理的土壤样本采用SVM模型训练集测试的实测值与预测值结果见

图3 SVM算法模型训练集对土壤含水率的预测结果
Fig.3 Estimation results of soil moisture content training set based on SVM algorithm model with different initial moisture content
A:初始含水率为15% Initial moisture content of 15%;B:初始含水率为20% Initial moisture content of 20%;C:初始含水率为25% Initial moisture content of 25%;D:初始含水率为30% Initial moisture content of 30%.
编号Number | 决定系数 | 平均相对误差MRE | 样本数量Sample size |
---|---|---|---|
T1 | 0.996 | 0.001 393 9 | 100 |
T2 | 0.992 | 0.002 950 4 | 100 |
T3 | 0.995 | 0.001 812 8 | 100 |
T4 | 0.991 | 0.002 902 6 | 100 |
由4个不同初始含水率条件下的超声波速度-含水率SVM算法模型训练集对土壤含水率预测结果可知,基于SVM算法的超声波声速模型对土壤含水率具有较好的识别效果(

图4 SVM算法模型测试集对土壤含水率的预测结果
Fig.4 Estimation results of soil moisture content test set based on SVM algorithm model with different initial moisture content
A:初始含水率为15% Initial moisture content of 15%;B:初始含水率为20% Initial moisture content of 20%;C:初始含水率为25% Initial moisture content of 25%;D:初始含水率为30% Initial moisture content of 30%.
笔者所在课题组提出了考虑干湿交替条件的超声波速度-土壤含水率二项式回归分析模型,具体参数如文献[

图5 2种方法预测结果平均相对误差图
Fig.5 Average relative error graph of prediction results of two methods
本研究以耕作层土壤为研究对象,考虑了干湿交替处理下土壤含水率变化对超声波速度的影响,并通过利用机器学习的SVM算法进行数据处理。引入了支持向量机的回归分析方法对4种不同干湿交替超声波速度与土壤含水率进行回归训练,构建了基于SVM优化算法的超声波速度的土壤含水率预测模型。将土壤含水率预测模型应用于不同墒情需求的土壤试验地进行含水率估测,并得到运用支持向量机模型对土壤含水率的预测结果的平均相对误差均小于1.5%。与运用二项式模型对土壤含水率的预测结果相比更加准确。本研究结果表明,土壤含水率与声学表征参数之间并非是完全的二次式关系,利用输入指标变量的二项式组合无法很好地描述输入和输出变量的依赖关系。而研究提出的SVM算法模型较为准确地反映了土壤含水率与超声波速度之间的非线性关系,表明利用支持向量机模型构建的基于超声波速度的土壤含水率预测网络模型更能够描述研究区域内不同墒情需求的土壤含水率。
该模型考虑了不同干湿交替状态对土壤结构产生的效应,在表征基于不同墒情要求的土壤含水率与超声波速度的关联性上更加精准,可为利用超声学信息对农田土壤含水率的持续监测提供参考。但同时也存在一些不足,忽略了土壤其他参数对声学测量的影响,导致准确度下降。在接下来的试验中,将考虑多个土壤参数对声参量测量的影响,以此修正土壤含水率的声学模型,实现更加精准统一的土壤水分声学检测模型。
参考文献 References
刘济明,李佳,文爱华,等.米槁幼苗光合色素与光合特征对干旱胁迫的响应[J].江苏农业科学,2019,47(9):171-174.LIU J M, LI J,WEN A H,et al.Responses of photosynthetic pigment and photosynthetic characteristics of Cinnamomum migao to drought stress[J].Jiangsu agricultural sciences,2019,47(9):171-174 (in Chinese). [百度学术]
梁飞,李云霞,关新元,等.长期大田滴灌对新疆农田水盐环境影响的研究进展[J].华中农业大学学报,2023,42(5):1-11.LIANG F,LI Y X,GUAN X Y,et al. Progress of effects of long-term field drip irrigation on water and salt balance of farmland in Xinjiang[J]. Journal of Huazhong Agricultural University,2023,42(5):1-11 (in Chinese with English abstract). [百度学术]
SANTI E,PALOSCIA S,PETTINATO S,et al.Application of artificial neural networks for the soil moisture retrieval from active and passive microwave spaceborne sensors[J].International journal of applied earth observation and geoinformation,2016,48:61-73. [百度学术]
李君,潘丽敏,潘翠敏,等.干湿交替下基于超声波速度的土壤含水率估算模型[J].农业工程学报,2020,36(8):114-119.LI J,PAN L M,PAN C M,et al.Model for estimating soil moisture content based on ultrasound velocity under dry-wet alternating conditions[J].Transactions of the CSAE,2020,36(8):114-119 (in Chinese with English abstract). [百度学术]
张国祥,申丽霞,郭云梅.微润灌溉条件下土壤质地对水分入渗的影响[J].灌溉排水学报,2016,35(7):35-39.ZHANG G X,SHEN L X,GUO Y M.Effect of soil structure on water infiltration under moistube irrigation[J].Journal of irrigation and drainage,2016,35(7):35-39 (in Chinese with English abstract). [百度学术]
XU Y,DUAN J L,JIANG R,et al.Study on the detection of soil water content based on the pulsed acoustic wave (PAW) method[J].IEEE access,2021,9:15731-15743. [百度学术]
WEI B X,HUANG Z.Assessment of change in acoustic wave velocity of compacted expansive soil through experiments[J].Scientia iranica,2017,24(1):136-142. [百度学术]
PAN L M,CHEN Y Y,XU Y,et al.A model for soil moisture content prediction based on the change in ultrasonic velocity and bulk density of tillage soil under alternating drying and wetting conditions[J/OL].Measurement,2022,189:110504[2022-12-20].https://doi.org/10.1016/j.measurement.2021.110504. [百度学术]
XU Y,LI J,DUAN J L,et al.Soil water content detection based on acoustic method and improved Brutsaert’s model[J/OL].Geoderma,2020,359:114003[2022-12-20].https://doi.org/10.1016/j.geoderma.2019.114003. [百度学术]
李君,徐岩,姜锐,等.超声波土壤含水量检测装置的模型建立与验证[J].农业工程学报,2017,33(13):127-133.LI J,XU Y,JIANG R,et al.Establishment and verification of model for ultrasonic soil water content detector[J].Transactions of the CSAE,2017,33(13):127-133 (in Chinese with English abstract). [百度学术]
李志国,白瑞涵,刘旭进,等.基于支持向量机上海地区土体物理力学指标相关性研究[J].地震科学进展,2023,53(2):66-76.LI Z G,BAI R H,LIU X J,et al.Correlations between physical and mechanical property indexes of Shanghai soil based on support vector machine[J].Progress in earthquake sciences,2023,53(2):66-76 (in Chinese with English abstract). [百度学术]
李铁牛,胡宾鑫,李化坤,等.基于改进支持向量机的微震初至波到时自动拾取方法[J].工矿自动化,2023,49(3):63-69.LI T N,HU B X,LI H K,et al.Automatic picking method of microseismic first arrival time based on improved support vector machine[J].Journal of mine automation,2023,49(3):63-69 (in Chinese with English abstract). [百度学术]
王辰璇,陈莉.小波-智能优化支持向量机模型及其对农业可持续发展的评价[J].农业工程学报,2023,39(4):208-216.WANG C X,CHEN L.Wavelet-intelligent optimization model using support vector machine and its application for evaluating sustainable agricultural development[J].Transactions of the CSAE,2023,39(4):208-216 (in Chinese with English abstract). [百度学术]
马玉莹,雷廷武,张心平,等.体积置换法直接测量土壤质量含水率及土壤容重[J].农业工程学报,2013,29(9):86-93.MA Y Y,LEI T W,ZHANG X P,et al.Volume replacement method for direct measurement of soil moisture and bulk density[J].Transactions of the CSAE,2013,29(9):86-93 (in Chinese with English abstract). [百度学术]
邹文安,姜波,顾李华.土壤水分常数的测定[J].水文,2015,35(4):62-66.ZOU W A,JIANG B,GU L H.Measurement of soil moisture constants[J].Journal of China hydrology,2015,35(4):62-66 (in Chinese with English abstract). [百度学术]
CHEN J,WANG H,YAO Y P.Determination of optimum moisture content of a clay using nonlinear ultrasonic technique[J].Japanese geotechnical society special publication,2015,1(2):1-5. [百度学术]
黄传琴,邵明安.干湿交替过程中土壤胀缩特征的实验研究[J].土壤通报,2008,39(6):1243-1247.HUANG C Q, SHAO M A.Experimental study on soil shrinking and swelling characteristics during the alternative drying and wetting processes[J].Chinese journal of soil science,2008,39(6):1243-1247 (in Chinese with English abstract). [百度学术]
WYLLIE M R J,GREGORY A R,GARDNER G H F.An experimental investigation of factors affecting elastic wave velocities in porous media[J].Geophysics,1958,23(3):459-493. [百度学术]
王浩,董振振,赵景波,等.基于改进自适应加权融合算法的土壤湿度监测研究[J].中国农机化学报,2019,40(1):152-155.WANG H,DONG Z Z,ZHAO J B,et al.Soil moisture monitoring based on improved adaptive weighted fusion algorithm[J].Journal of Chinese agricultural mechanization,2019,40(1):152-155 (in Chinese with English abstract). [百度学术]
仝玉婷.最小二乘支持向量回归机的算法研究[D].金华:浙江师范大学,2018.TONG Y T.Study of least squares support vector regression[D].Jinhua:Zhejiang Normal University,2018 (in Chinese with English abstract). [百度学术]