深松作业能够打破犁底层,耕层土壤松散结构易于蓄水、脱盐碱、疏通空气、增强热交换等[1-3], 有利于增强土壤蓄水保墒能力, 促进农作物根系下扎[4-6],有利于土壤养分转化利用,有效降低耕地土壤容重和紧实度, 提高土壤蓄水保墒能力, 增强抗旱抗涝能力,促进农作物增产增收[7-8]。
近年来,国家加快推进深松装备应用普及, 全国深松机拥有量持续增加,各地研制并推广高精度、信息化农机深松远程监测装备,通过超声波、角度传感器等设备测量作业深度来准确掌握机械深松作业质量[9-11],并利用信息化手段监测作业质量与面积, 促进了“互联网+”农机技术应用,为农业生产主体和农机化主管部门决策提供了科学依据,显著提升了农机作业管理信息化水平和效率[12-14]。通过信息技术手段对深松作业进行实时监测和远程监管取得了明显成效,为深松作业监管提供了量化依据[15-17]。
作业质量评价是农业机械化生产的重要组成部分,可为农业机械选型和配套提供科学依据,提高农业机械使用效率和经济效益。现有深松作业质量研究聚焦于作业面积、作业功耗、任务分配优化、效益评价等。刘阳春等[18]利用距离幅宽方法实现了深松作业面积计算,满足了深松作业远程监测和管理的基本需求;白晋[19]利用网格标记法建立深松作业漏耕面积计算模型,农田实验结果表明,模型最大相对误差为5.49%,能够满足深松作业系统要求;李亚硕等[20]基于半监督BP_Adaboost方法对农机效益进行了评估,并推荐最优农机机具组合,增强农机作业能力,提高农机作业效益;王培等[21]针对深松作业场景,基于线性回归分析了拖拉机功率与深松作业效率关系,建立了深松作业效率估算模型;郑侃等[22]运用离散元仿真分析了深松、旋耕作业次序对其作业质量及功耗的影响,建立了功耗回归方程,并分析了作业质量综合分析;马明龙[23]提出了基于倾角传感器与三点悬挂的几何深度算法和基于面积格点覆盖的面积测算方法,实时深度和面积测算误差均不超过5%,对研究农机监测技术提供了新思路。
现有深松作业质量评价方法无法全面覆盖作业质量、效率、成本等评价要素,也难以实现更为科学的作业评价,大部分研究仅把农机作业面积、达标率、平均深度等作为深松作业质量指标[19-20]。目前信息化监管系统仅以作业面积和深度达标率作为深松作业质量的评判依据,缺乏对整个农机作业过程的综合评价。深松作业质量是深松作业的重要考核指标,包含农机作业过程中数量、重漏等作业关键信息,本文分析深松作业面积、作业效率、作业速度、作业历史重叠、作业时间利用率等,以实现深松作业质量综合评价。
深松作业是一项关系到土壤理化性质、作物生长等多方面的农业生产技术,深松作业质量直接影响作物产量和质量;农业机械作业质量评价应遵循全面性原则,从作业精度、效率和成本等多个维度进行综合评级[24-26]。
目前利用信息化手段可以对深松作业过程中的多项指标进行精准监控,综合反映深松作业精度、效率和成本。通过计算特征表达作业精度、效率和成本。通过对深松作业中的监控数据进行分析,建立农机作业质量综合评价方法。该方法对深松作业地块、作业面积、达标面积、作业时长和总时长等进行融合分析,自动计算深松作业特征向量,具体包括数据采集、数据获取及预处理、农机作业质量特征向量计算和深松作业质量综合评价4个步骤。
数据采集使用北京市农林科学院农业智能装备技术研究中心研发的农机智能监测终端,用于农机作业全程机械化作业监控。农机智能监测终端主要包括全球导航卫星系统(Global navigation satellite system,GNSS)模块、无线传输模块、作业状态传感器模块和作业深度传感器模块,如图1所示。其中GNSS模块定位精度为2.5 m圆概率误差(Circular error probable,CEP),速度精度为0.1 m/s;利用姿态传感器,建立基于角度关系的耕深检测模型,作业深度检测均方根误差小于0.64 cm[3];采集农机作业空间运行轨迹数据和作业数据,包括UTC时间、经纬度、海拔、速度、方位角、作业深度等信息;无线传输模块将采集的空间运行轨迹信息和作业信息无线传输至服务器。
图1 农机智能监测终端
Fig.1 Intelligent monitor terminal of agricultural machinery
1.监测主机(包含无线传输模块) 2.GNSS定位模块 3.作业状态传感器
农机作业状态有工作状态和非作业状态,工作状态是农机在田块内的实际工作状态,非作业状态包括停歇和地头转弯等状态。农机作业状态和作业深度有直接关系,当作业深度大于0为工作状态,当作业深度小于等于0时为非作业状态。从数据库中获取1次在农田地块作业A中数据记录集合P,即
P={p1,p2,…,pi}
pi为第i个记录点数据,包括GNSS时间、经纬度、作业深度、作业状态等信息,并按照GNSS时间排序,根据深度计算作业状态。
采用通用横轴墨卡托投影(Universal transverse Mercator projection,UTM)将经纬度坐标转换为平面坐标,为数据平面分析计算做准备。
为了方便长度计算,需要对数据集P中的pi经纬度进行投影转换,数据统一在平面坐标系下,获得相应平面坐标为
pi={t, lon, lat, d, o, x, y}
式中 x——平面坐标横坐标
y——平面坐标纵坐标
t——GNSS时间
lon——经度
lat——纬度
d——作业深度
o——作业状态
农机作业质量特征向量是指农机在农田矢量地块B中深松作业时农机作业面积、效率、覆盖率、深度达标率、作业时间利用率等特征,用于深松作业质量综合评价。
(1)农机作业面积
按照时间顺序连接集合P内作业状态为田间作业状态的所有点,得到农机作业路径长度L,以农机深松机具幅宽W的一半为缓冲半径,对L进行缓冲操作,得到矢量面A,A大小为农机作业面积,计算式为
(1)
(2)农机作业效率
农机作业效率F是评估农机作业效果的核心指标,是指单位时间完成的作业面积,计算式为
(2)
式中 Aarea——农机作业面积,hm2
T——农机作业时长,h
(3)农机作业覆盖率
对地块B和作业矢量面A进行矢量减操作,得到作业覆盖区域面积,农机作业覆盖率S是作业覆盖区域面积占地块面积百分比,计算式为
(3)
式中 Barea——矢量地块B面积,hm2
(4)深度达标率
深度达标率H是指深松作业达标深度区域面积(超过预设达标深度区域,达标深度是区域土壤性状决定,一般不低于25 cm)占作业面积的百分比,计算式为
(4)
式中 Carea——作业深度高于预设达标深度区域面积,hm2
(5)农机作业时间利用率
农机作业时间利用率U是指农机作业时间占农机总运行时间(包括作业、转运、停歇等)百分比,计算式为
(5)
式中 Tw——农机作业时间,h
Tt——农机总运行时间,h
深松作业质量是农机作业质量特征的综合表现,通过归一化对农机作业向量特征建立特征向量权重矩阵,计算深松作业质量指数,实现作业质量综合评价。农机作业质量特征向量是指农机在农田矢量地块B中深松作业时农机作业面积、效率、覆盖率、深度达标率、作业时间利用率等特征,用于深松作业质量综合评价。
(1)农机作业质量特征归一化
归一化农机作业质量特征是将不同特征之间的数值范围统一到同一尺度上,消除不同特征间的量纲影响,使用最小-最大归一化法,统计历史农机作业向量特征,获取每个特征向量最大值和最小值,对特征向量x进行归一化处理得x′,即
(6)
式中 xmin——特征向量x最小值
xmax——特征向量x最大值
将每台农机作业面积A、农机作业效率F、农机作业覆盖率S、农机作业深度达标率H、农机作业时间利用率U等组成农机作业矩阵X,即
X=[A F S H U]
(7)
根据式(6)对农机作业矩阵进行归一化处理,得到归一化后特征向量
X′=(A′,F′,S′,H′,U′)
(8)
(2)深松作业质量指数计算
为了评价农机作业综合质量,建立深松作业质量权重,所有特征向量权重组成矩阵权重Y
Y=[wA wF wS wH wU]
(9)
其中
wA+wF+wS+wH+wU=1
式中 wA——A权重 wF——F权重
wS——S权重 wH——H权重
wU——U权重
权重依照各特征向量属性贡献程度使用德尔菲法来确定,通过匿名方式征询有关专家意见,对专家意见进行统计、处理、分析和归纳,客观地综合多数专家经验与主观判断,经过意见征询、反馈和调整获得最终权重[26]。
深松作业质量指数(AMQI)是表征深松作业质量的综合评价指数,是特征向量和权重矩阵转置矩阵乘积,计算式为
PAMQI=X′YT
(10)
黑龙江省红星农场位于小兴安岭南麓,位于黑龙江省北安市境内,隶属黑龙江农垦总局北安分局(北纬48°02′~48°17′,东经126°47′~127°15′)如图2所示,属寒温带大陆性季风气候,四季分明,全年平均降水量555.3 mm,年大于等于10℃活动积温约为2 293.1℃,农场耕地2.73×104 hm2,多为丘陵漫岗,土质肥沃,宜于耕作,并极有利于农作物生长,机械总动力2.7×104 kW,机械化率达99%以上。
图2 研究区概况
Fig.2 Overview of study area
为了验证深松作业质量指数,2022年1月1日—12月31日,在红星农场约翰迪尔、福格森、凯斯、维美德等33台拖拉机上安装车载智能终端,采集深松作业数据,累计采集457车日数据,面积总计1.47×104 hm2,每台车年度作业累计深松作业面积为118.30~2 251.54 hm2,作业时长为28.23~441.94 h,如表1所示。
表1 农机作业数据
Tab.1 Statistics of agricultural machinery
序号地块面积/hm2作业面积/ hm2达标面积/hm2作业时长/h总时长/h1234.12220.82220.3048.1560.562450.65433.14431.8291.76113.043616.62583.20581.46131.49158.074298.46271.88271.3566.1978.895171.49150.90114.3835.2141.876128.35128.43123.8131.0235.317151.32136.21135.9438.7046.558237.82271.03270.5784.1396.079152.39159.73152.9440.4950.5510146.88165.17161.5339.7149.3611279.76269.81269.2056.4367.4012477.34418.72416.6987.90105.3713455.27413.70408.5989.61106.4114432.74392.49391.9481.6597.2215198.42231.45230.2567.6885.1616172.82204.51204.1954.4269.4017424.37465.10463.17109.80134.0718112.18133.28133.1235.2740.33191530.191778.461768.06330.64385.3220126.65118.30117.7828.2333.9121419.74377.58361.16100.22111.65221276.141480.321477.36281.96288.59232028.412251.542233.75441.94482.54241098.701281.971253.78264.25301.2625132.69147.92147.7637.1042.8426312.75333.34330.2691.09109.8127166.15159.21137.6040.7747.4828320.38289.68289.0964.1174.0729355.09286.74279.3571.5488.4030339.49335.84333.2862.4682.8331472.03453.92448.04115.38128.6332142.92145.58145.3435.6040.4033242.33240.06238.9548.5356.46
通过农机作业质量特征向量计算及归一化,实现农机作业质量综合评价,根据式(1)~(5)计算农机作业面积A、农机作业效率F、农机作业覆盖率S、农机作业深度达标率H、农机作业时间利用率U等,并根据式(6)进行特征向量归一化。
使用专家打分法,确定各个特征向量权重,发送调查问卷给咨询专家,专家包括农业科研院所、政府监管部门、农业服务组织和一线农技推广人员等;根据专家打分结果,并求平均值获得各个特征向量权重,wA、wF 、wY 、wH 、wU 分别为0.2、0.1、0.3、0.3、0.1。农机作业覆盖率和作业达标率更受关注,其次为农机作业面积,最后为农机作业效率和时间利用率。
根据式(10)获得深松作业质量指数(AMQI),来表征深松作业综合质量,结果如表2所示。从表2可看出,AMQI为0.15~0.90,平均值为0.56,方差为0.14,其中序号22和序号5分别为最大值和最小值。
表2 农机特征向量归一化结果
Tab.2 Results of feature vector normalization for agricultural machinery operations
序号作业面积作业覆盖率深度达标率作业效率时间利用率AMQI10.050.360.990.630.180.5220.150.400.990.690.260.5430.220.360.990.560.350.5440.070.271.000.410.380.4750.020.190.000.490.390.1560.000.510.860.430.560.5170.010.241.000.140.350.4280.070.871.000.000.550.6390.020.630.830.340.210.50100.020.830.910.430.230.59110.070.411.000.720.370.55120.140.180.980.710.360.49130.140.270.950.650.400.50140.130.261.000.730.380.52150.050.940.980.090.180.62160.040.991.000.250.140.64170.160.760.990.470.290.63180.011.001.000.260.540.68190.780.930.981.000.470.88200.000.330.990.450.350.48210.120.240.820.250.640.43220.640.931.000.941.000.90231.000.790.970.870.730.89240.550.940.910.760.550.80250.010.811.000.350.500.63260.100.680.970.200.340.57270.020.400.440.320.470.33280.080.251.000.600.500.50290.080.000.900.360.250.35300.100.480.971.000.000.56310.160.400.950.330.640.54320.010.551.000.400.570.57330.060.480.990.800.470.58
将最大值和最小值导入雷达图中,展现数据集的多个特征向量,如图3所示。从图3可看出,AMQI最大值(序号22)作业覆盖率、深度达标率、作业效率和时间利用率等特征向量占据优势,仅作业面积处于中等;AMQI最小值(序号5)作业面积、覆盖率、达标率等特征向量都处于低位,而作业效率、时间利用率处于中等;结合专家打分表,作业覆盖率和深度达标率所占权重最大,而作业效率、时间利用率所占权重最小。
图3 AMQI极值雷达图
Fig.3 Radar map of extreme value of AMQI
当前深松作业质量以达标面积为主,从表1、2可以看出,序号23达标作业面积(2 233.75 hm2)最高,但深度达标率(99.21%)、时间利用率(91.59%)、归一化作业覆盖率(0.79)、作业效率(5.09 hm2/h)均低于序号22。
采用模型计算序号23、序号22的AMQI分别为0.89、0.90,序号22总体作业质量较高,也说明AMQI指数评价法不仅涵盖了作业量(达标面积),也考虑了作业深度(达标率)、作业覆盖率(遗漏)、作业效率、运营成本(时间利用率)等多个维度的综合分析,较传统评价方法更为科学全面。
利用车载智能终端采集深松作业数据,计算作业面积、作业效率、作业覆盖率、作业深度达标率、作业时间利用率等深松作业特征向量并进行归一化处理,融合了研究区域作业水平,将量纲表达式变为无量纲表达式,能够方便地进行比较和加权,从而构建作业质量指数评价综合作业质量,与达标面积评价法相比,评级维度更多,且评价更加全面和科学,从实验结果来看,归一化特征向量反映农机特征大小分布,AMQI也反映农机的深松作业综合质量,是农机作业面积、作业效率、作业覆盖、作业深度达标率、作业时间利用率等的综合表现,可以为农机作业补贴提供数据支撑。
[1] 郭志军,佟金,周志立,等.深松技术研究现状与展望[J].农业工程学报, 2001,17(6):169-174.
GUO Zhijun,TONG Jin,ZHOU Zhili,et al.Review of subsoiling techniques and their applications[J].Transactions of the CSAE,2001,17(6):169-174.(in Chinese)
[2] 刘学丽.机械化深松整地作业技术规范与应用[J].中国农机装备,2024(3):29-31.
[3] 尹彦鑫,王成,孟志军,等.悬挂式深松机耕整地耕深检测方法研究[J].农业机械学报,2018,49(4):68-74.
YIN Yanxin,WANG Cheng,MENG Zhijun,et al.Operation quality measurement method for tilling depth of suspended subsoiler[J]. Transactions of the Chinese Society for Agricultural Machinery,2018,49(4):68-74.(in Chinese)
[4] 李洪文,陈君达,李问盈.保护性耕作条件下深松技术研究[J].农业机械学报,2000,31(6):42-45.
LI Hongwen,CHEN Junda,LI Wenying.Study on subsoiling technique for conservation tillage field[J]. Transactions of the Chinese Society for Agricultural Machinery,2000,31(6):42-45.(in Chinese)
[5] 王慧杰,郝建平,冯瑞云,等.微孔深松耕降低土壤紧实度提高棉花产量与种籽品质[J].农业工程学报,2015,31(8):7-14.
WANG Huijie,HAO Jianping,FENG Ruiyun,et al.Microhole subsoiling decreasing soil compaction, and improving yield and seed quality of cotton[J]. Transactions of the CSAE,2015,31(8):7-14.(in Chinese)
[6] 李荣,侯贤清.深松条件下不同地表覆盖对马铃薯产量及水分利用效率的影响[J].农业工程学报,2015,31(20):115-123.
LI Rong,HOU Xianqing.Effects of different ground surface mulch under subsoiling on potato yield and water use efficiency[J].Transactions of the CSAE,2015,31(20):115-123.(in Chinese)
[7] 黄玉祥,杭程光,李伟,等.深松作业效果试验及评价方法研究[J].西北农林科技大学学报(自然科学版),2015, 43(11):228-234.
HUANG Yuxiang,HANG Chengguang,LI Wei,et al.Subsoiling test and evaluation methodology of tillage quality[J]. Journal of Northwest A&F University(Natural Science Edition),2015,43(11):228-234.(in Chinese)
[8] 张强,晋清源.模拟深松条件下深施磷肥对旱地冬小麦根系生长和水肥效应的影响[J].农业工程学报,2000,16(7): 57-60.
ZHANG Qiang,JIN Qingyuan.Effects of phosphorous fertilizer deep application on winter wheat root and soil water under soil loosening conditions[J].Transactions of the CSAE,2000,16(7):57-60.(in Chinese)
[9] 张祺睿.深松作业深度监测系统的设计与试验[D].广州:华南农业大学,2021.
ZHANG Qirui.Design and experiment of monitoring system for subsoiling depth[D].Guangzhou:South China Agricultural University,2021.(in Chinese)
[10] 李海涛.农用拖拉机深松深度作业监测系统设计——基于STM32U575[J].农机化研究,2024,46(9):128-131,137.
LI Haitao.Design of monitoring system for subsoiling depth operation of agricultural tractor based on STM32U575[J].Journal of Agricultral Mechanization Research,2024,46(9):128-131,137.(in Chinese)
[11] 赵世卿,吴永峰,李光林,等.丘陵山区农机深松作业深度和面积远程监测系统[J].西南大学学报(自然科学版),2018,40(12):22-29.
ZHAO Shiqing,WU Yongfeng,LI Guanglin,et al.A remote monitoring system for depth and area of agricultural machinery subsoiling in hilly regions[J].Journal of Southwest University (Natural Science Edition),2018,40(12):22-29.(in Chinese)
[12] WANG Z,CHEN L,LIU Y.Design and implementation of agricultural machinery monitoring and scheduling system[J]. Computer Engineering,2010,36(11):232-234,237.
[13] 冀福华.农机田间作业大数据处理关键技术研究及平台构建[D].北京:中国农业机械化科学研究院,2021.
JI Fuhua.Research on key technology and platform construction of agricultural machinery field operation big data processing[D].Beijing:Chinese Academy of Agricultural Mechanization Sciences,2021.(in Chinese)
[14] 吴才聪,陈瑛,杨卫中,等.基于北斗的农机作业大数据系统构建[J].农业工程学报,2022,38(5):1-8.
WU Caicong, CHEN Ying, YANG Weizhong,et al.Construction of big data system of agricultural machinery based on BeiDou[J].Transactions of the CSAE,2022,38(5):1-8.(in Chinese)
[15] 史国滨,王熙.基于ASP.NET的农机监控Web GIS系统性能优化[J].安徽农业科学,2011,39(5):2821-2823.
SHI Guobin,WANG Xi.Optimization of performance of agricultural machinery monitoring Web GIS system based on ASP.NET[J].Journal of Anhui Agricultural Science,2011,39(5):2821-2823.(in Chinese)
[16] 孟志军,尹彦鑫,罗长海,等.农机深松作业远程监测系统设计与实现[J].农业工程技术,2018,38(18):34-37.
[17] 翟长远,杨硕,王秀,等.农机装备智能测控技术研究现状与展望[J].农业机械学报,2022,53(4):1-20.
ZHAI Changyuan,YANG Shuo,WANG Xiu,et al.Status and prospect of intelligent measurement and control technology for agricultural equipment[J].Transactions of the Chinese Society for Agricultural Machinery,2022,53(4):1-20.(in Chinese)
[18] 刘阳春,苑严伟,张俊宁,等. 深松作业远程管理系统设计与试验[J].农业机械学报,2016,47(增刊):43-48.
LIU Yangchun,YUAN Yanwei,ZHANG Junning,et al.Design and experiment of remote management system for subsoiler[J]. Transactions of the Chinese Society for Agricultural Machinery,2016,47(Supp.):43-48.(in Chinese)
[19] 白晋.农机深松作业智能监测系统[D].太原:中北大学,2023.
BAI Jin.Intelligent monitoring system of agricultural machinery subsoiling operation[D].Taiyuan:North University of China,2023.(in Chinese)
[20] 李亚硕,赵博,徐名汉,等.基于半监督BP_Adaboost的农机作业效益评估[J].农业工程学报,2023,39(23):67-74.
LI Yashuo, ZHAO Bo, XU Minghan,et al.Evaluating operation benefit of agricultural machinery using semi-supervised BP_Adaboost[J].Transactions of the CSAE,2023,39(23):67-74.(in Chinese)
[21] 王培,孟志军,安晓飞,等.拖拉机功率与深松作业效率关系研究[J].农业机械学报,2019,50(增刊):87-90.
WANG Pei,MENG Zhijun,AN Xiaofei,et al.Relationship between agricultural machinery power and agricultural machinery subsoiling operation[J].Transactions of the Chinese Society for Agricultural Machinery,2019,50(Supp.):87-90.(in Chinese)
[22] 郑侃,何进,李洪文,等.作业次序对深松旋耕联合作业机作业质量及功耗的影响[J].农业工程学报, 2017, 33(21):52-60.
ZHENG Kan,HE Jin,LI Hongwen,et al.Influence of working order on working quality and power consumption of subsoiling and rotary tillage combined machine[J].Transactions of the CSAE,2017,33(21):52-60.(in Chinese)
[23] 马明龙.机械深松作业质量监测技术研发[D].杨凌:西北农林科技大学,2017.
MA Minglong.Research and development of mechanized subsoiling operation quality monitoring technology[D]. Yangling:Northwest A&F University,2017.(in Chinese)
[24] 金慧芳,史东梅,陈正发,等.基于聚类及PCA分析的红壤坡耕地耕层土壤质量评价指标[J].农业工程学报,2018, 34(7):155-164.
JIN Huifang,SHI Dongmei,CHEN Zhengfa,et al.Evaluation indicators of cultivated layer soil quality for red soil slope farmland based on cluster and PCA analysis[J].Transactions of the CSAE,2018,34(7):155-164.(in Chinese)
[25] 王政.基于客户满意度的玉米收割机作业调度系统开发[D].杨凌:西北农林科技大学,2023.
WANG Zheng.Development of a corn harvester operation scheduling system based on customer satisfaction[D]. Yangling:Northwest A&F University,2023.(in Chinese)
[26] 孙琳岚.农业科技发展水平测度指标与方法研究[D].合肥:安徽农业大学,2011.
SUN Linlan.Study on indicators system and method for the level of development of agricultural science and technology[D].Hefei:Anhui Agricultural University,2011. (in Chinese)