本发明涉及基于定点调查和区域普查的食叶类农业虫害发生风险的评估技术，属于人工智能研究技术领域。
背景技术：
有害生物风险分析工作始于19世纪中期，最初是针对有害生物在欧洲农作物上的猖獗流行而开展起来的，当时仅对植物上可能携带的有害生物进行简单评估，仅考虑植物和有害生物的个别生物学特性和个别风险因素。自20世纪30年代开始，随着人们对有害生物风险的认识逐渐深入，与有害生物定殖可能性有关的区域适生性问题开始引起重视；20世纪70-80年代，研究的注意力转向有害生物传入和扩散可能性、潜在的经济影响及环境影响等综合因素的评估；20世纪90年代以后，有害生物风险分析逐步进入规范化阶段。为协调世界各国的有害生物风险分析，联合国粮农组织(fao)/国际植物保护公约(ippc)制定了有害生物风险分析准则以及检疫性有害生物风险分析标准。到目前为止，各国进行的有害生物风险分析工作，均是以定性分析方法为主，相关定量分析仍处于尝试阶段。
食叶害虫是以叶片为食的害虫，主要危害健康植物，以幼虫取食叶片，常咬成缺口或仅留叶脉，甚至全部吃光，是当前危害农、林业生产的主要害虫。食叶虫害发生的主要特点是虫口密度大而集中，在农业生产中，一旦发生食叶虫害往往会造成作物的大面积减产，使农户们蒙受重大经济损失。鉴于此类农业虫害的巨大危害性，有必要加强虫害风险防控方面的研究，以减低虫害发生率及危害程度。
目前，国家相关部门或省市尺度上关于此类害虫的分布及危害程度的调查周期长，时效性不够，数据集成度不高。而田间调查尽管详实，但是区域代表性不够，不易推广应用。因此，急需一种针对食叶类农业虫害的风险评估方法，初步实现预报预防的目的，以降低此类农业虫害的危害性。
技术实现要素：
本发明的目的是提供一种农业虫害的风险评估技术，为有效开展有害生物灾害的预防管控提供技术支撑。
本发明的食叶类农业虫害风险评估方法，包括如下步骤：
步骤1、建立风险评估模型公式：r＝d×e×p；
其中r：害虫灾害风险；d：致灾因子的破坏力；e：生物致灾因子中的暴露度；p：孕灾环境。
步骤2、选定待评估风险区域，计算生物致灾因子中的暴露度。
步骤2.1、分别统计待评估风险区域内的当年播种面积、近20年最大播种面积、近20年最小播种面积，并根据如下公式计算：
cr为农作物播种面积指数，将计算获得的cr值查表得到对应标准值cr标。
表1农作物播种面积指数划分标准表
步骤2.2、分别统计待评估风险区域内的当年粮食预估产量、近20年最大粮食总产量、近20年最小粮食总产量，并根据如下公式计算：
gr为粮食产量指数，将计算获得的gr值查表得到对应标准值gr标，
表2粮食总产量指数划分标准表
根据公式e＝(cr标+gr标)/2计算得到生物致灾因子中的暴露度e。
步骤3、致灾因子的破坏力计算；
步骤3.1在待评估风险区域内选定n个采样的样方，每个样方的面积均为1m2，n个采样样方的面积总和占待评估风险区域总面积的0.5％-1.0％。
步骤3.2、食叶类害虫发生指数。
对待评估风险区域内选定的所有样方中的植株计数，总数记为受调查的作物植株总数，并统计有虫害的作物植株数，根据公式计算。
oc＝(有虫害的作物植株数/受调查的作物植株总数)×100％
oc为食叶类害虫发生指数。
将计算获得的oc值查表得到发生级别oc发生，
表3害虫发生等级划分标准表
步骤3.3、食叶类害虫危害指数。
在每个样方内分别选择3～5株病株农田作物，将选出的所有病株农田作物作为病株样本，在病株样本上取等数量的健康叶片和病害叶片，采用叶面积仪扫描并求和得到健康叶片的叶片总面积以及病害叶片的叶片总面积。
根据公式计算：
da为食叶类害虫危害指数；
将计算获得的da值查表得到危害级别da危害，
表4害虫危害等级划分标准表
步骤3.4、虫害发生的单株产量损失率。
在每个样方内分别选取3～5株病株农田作物和3～5株未受害农田作物，对每个选取的作物收割并记录收割获得的作物量，分别计算病株农田作物的单株平均作物量以及未受害农田作物的单株平均作物量，根据如下公式计算：
lo为虫害发生的单株产量损失率。
将计算获得的lo值查表得到危害级别lo损失，
表5单株产量损失率等级划分标准表
步骤3.5、根据公式d＝(oc发生+da危害+lo损失)/3计算得到致灾因子的破坏力d。
步骤4、孕灾环境计算。
记待评估时间段内的预报日均气温为te，预报日降水量为pr，预报日相对湿度为hu，预报日平均风速为wi。
表6影响虫害发生的气象条件等级划分表
查表根据预报日均气温te、预报日降水量pr、预报日相对湿度hu、预报日平均风速wi所在区间，分别得到对应的等级，记为气温级别te等级，降水量级别pr等级，湿度级别hu等级，风速级别wi等级。
根据公式p＝0.3×te等级+0.3×hu等级+0.2×pr等级+0.2×wi等级计算得到孕灾气象风险值p。
步骤5、将步骤2、3、4获得的数据代入公式r＝d×e×p，计算得到害虫灾害风险r。
根据计算得到的r，判断害虫灾害风险的状态，r>80为高风险状态，r<40为低风险状态，r＝[40,80]为中度风险状态。根据风险程度即可确定应采取的风险应对措施，有效降低虫害造成的损失。
本发明的技术方案具有如下有益效果。
虫害发生风险评估被认为能够用来预测未来虫害的发生不利影响或者是评估过去某种因素导致虫害发生的可能性。本专利通过对致灾各要素的田间调查、采样、分析，可以初步确定虫害的发生区域和发生情况，再结合孕灾环境是否适合虫情发展，承载体是否支撑虫害大面积发生，预估虫害的发展动态。通过风险区判定和风险等级的计算，最终确定虫害的不同风险等级。在验证本发明方法过程中发现，采用过往数据进行验证时，无论是以乡镇一级还是以地市一级作为评价对象，本方法确定的虫害中高风险区域，其孕灾环境都适宜虫害的发生发展。进一步查阅虫害防治相关资料对比发现，该类地区无论是虫害发生等级还是受灾损失程度亦较高。表明本发明的方法可以较为准确的判定虫害的发生发展趋势与危害影响程度。
本发明的方法可以帮助相关植保部门了解和预测虫害发生的影响因素，有利于虫害防治的决策制定。在虫害可能爆发发生前，根据风险评估结果，对风险较高的区域及时采取科学有效的措施以控制虫情，将灾害损失减到最低。
附图说明
图1为本发明的食叶类农业虫害风险评估方法的模型图。
具体实施方式
本实施例的食叶类农业虫害风险评估方法，包括如下步骤：
本发明的食叶类农业虫害风险评估方法，包括如下步骤：
步骤1、建立风险评估模型公式：r＝d×e×p；
其中r：害虫灾害风险；d：致灾因子的破坏力；e：生物致灾因子中的暴露度；p：孕灾环境。
步骤2、选定待评估风险区域，计算生物致灾因子中的暴露度。
步骤2.1、分别统计待评估风险区域内的当年播种面积、近20年最大播种面积、近20年最小播种面积，并根据如下公式计算：
cr为农作物播种面积指数，将计算获得的cr值查表得到对应标准值cr标。
表1农作物播种面积指数划分标准表
步骤2.2、分别统计待评估风险区域内的当年粮食预估产量、近20年最大粮食总产量、近20年最小粮食总产量，并根据如下公式计算：
gr为粮食产量指数，将计算获得的gr值查表得到对应标准值gr标，
表2粮食总产量指数划分标准表
根据公式e＝(cr标+gr标)/2计算得到生物致灾因子中的暴露度e。
步骤3、致灾因子的破坏力计算；
步骤3.1在待评估风险区域内选定n个采样的样方，每个样方的面积均为1m2，n个采样样方的面积总和占待评估风险区域总面积的0.5％-1.0％。
步骤3.2、食叶类害虫发生指数。
对待评估风险区域内选定的所有样方中的植株计数，总数记为受调查的作物植株总数，并统计有虫害的作物植株数，根据公式计算。
oc＝(有虫害的作物植株数/受调查的作物植株总数)×100％
oc为食叶类害虫发生指数。
将计算获得的oc值查表得到发生级别oc发生，
表3害虫发生等级划分标准表
步骤3.3、食叶类害虫危害指数。
在每个样方内分别选择3～5株病株农田作物，将选出的所有病株农田作物作为病株样本，在病株样本上取等数量的健康叶片和病害叶片，采用叶面积仪扫描并求和得到健康叶片的叶片总面积以及病害叶片的叶片总面积。
根据公式计算：
da为食叶类害虫危害指数；
将计算获得的da值查表得到危害级别da危害，
表4害虫危害等级划分标准表
步骤3.4、虫害发生的单株产量损失率。
在每个样方内分别选取3～5株病株农田作物和3～5株未受害农田作物，对每个选取的作物收割并记录收割获得的作物量，分别计算病株农田作物的单株平均作物量以及未受害农田作物的单株平均作物量，根据如下公式计算：
lo为虫害发生的单株产量损失率。
将计算获得的lo值查表得到危害级别lo损失，
表5单株产量损失率等级划分标准表
步骤3.5、根据公式d＝(oc发生+da危害+lo损失)/3计算得到致灾因子的破坏力d。
步骤4、孕灾环境计算。
记待评估时间段内的预报日均气温为te，预报日降水量为pr，预报日相对湿度为hu，预报日平均风速为wi。
表6影响虫害发生的气象条件等级划分表
查表根据预报日均气温te、预报日降水量pr、预报日相对湿度hu、预报日平均风速wi所在区间，分别得到对应的等级，记为气温级别te等级，降水量级别pr等级，湿度级别hu等级，风速级别wi等级。
根据公式p＝0.3×te等级+0.3×hu等级+0.2×pr等级+0.2×wi等级计算得到孕灾气象风险值p。
步骤5、将步骤2、3、4获得的数据代入公式r＝d×e×p，计算得到害虫灾害风险r。