[发明专利]一种农产品产地安全预警方法

权利要求书

1.一种农产品安全加密调查方法，其特征在于，所述方法包括：

(1)超标概率空间分布

(1.1)建立产地环境数据库

(1.2)农产品历史超标率计算

基于产地环境数据库农产品重金属监测信息，依据《食品安全国家标准食品中污染物限量》判定农产品是否超标，并分别计算特定重金属监测指标的超标率p_i；

(1.3)土壤监测指标空间分布

结合待预警区域的区域特征进行空间插值计算，针对土壤监测指标分别进行空间插值，获取对应指标空间分布图；

(1.4)栅格数据获取

依据土壤监测指标插值结果，获取待预警区域栅格数量及各栅格内的所有土壤监测指标含量；

(1.5)超标概率区域划定

(1.5.1)基于待预警区域农产品超标率p_i与各栅格内的土壤监测指标含量x，通过QDA模型获取各栅格内的农产品超标概率q，QDA模型具体如下：

其中，μ为所有土壤监测指标均向量；∑为所有土壤监测指标协方差矩阵；当y＝0时模型结果为q₀，当y＝1时模型结果为q₁，农产品超标概率q计算公式具体如下：

(1.5.2)计算所有栅格内多种重金属农产品超标概率，选取多种重金属农产品超标概率在各栅格内的最大值，按照5-20％的概率梯度，从低到高将待预警区域划分为20-5个概率区域；

(2)加密调查

(2.1)点位布设

分别计算各个概率区域加密布点数量M_q，并采取均匀布点法、梅花形布点法、棋盘式布点法、蛇形布点法、对角线布点法、网格布点法中的一种或多种结合进行点位布设，点位数量计算公式具体如下：

M_q＝T(S/h²-t²Cv²/m²)+t²Cv²/m²

其中，M_q为q中位值概率区域内加密布设点位数量；T为区域概率系数，通过一维正态分布公式计算得到；S为区域面积，h为区域内各项重金属的最小变程，t为T检验的通过阈值，m为容忍误差，Cv为变异系数，且Cv＝σ/u，σ和u分别为概率区域内所有土壤重金属监测指标的标准差和平均值；

所述各概率区域内重金属变程通过半变异函数模型获得；

(2.2)土壤及农产品数据获取

分别获取各个概率区域的加密布点点位地理信息，进行协同样品采集及实验室检测工作，获得土壤及农产品监测指标检测结果，并上传存储至产地环境数据库。

2.如权利要求1所述的一种农产品安全加密调查方法，其特征在于，所述土壤监测指标包括：Cd、Hg、As、Pb、Cr、pH、CEC、SOM、Cu、Zn、Ni、Mg、Fe、Mn、Ca、Si、Se中的一种或多种；农产品监测指标包括：Cd、Hg、As、Pb、Cr中的一种或多种；所述加密布点存储信息包括：点位地理信息、采样时间、农产品类型、检测时间、检测结果中的一种或多种；收集待预警区域的历史土壤监测信息、背景值和农产品监测信息，建立产地环境数据库；所述土壤监测指标包括：Cd、Hg、As、Pb、Cr、pH、CEC、SOM中的一种或多种；背景值及农产品监测指标包括：Cd、Hg、As、Pb、Cr中的一种或多种。

3.如权利要求1所述的一种农产品安全加密调查方法，其特征在于，所述p_i的计算公式具体如下：

p_i＝n_i/N

其中，n_i为特定重金属监测指标i的超标点位数量，N为待预警区域内总点位数量；

所述空间插值方法包括但不限于：克里金插值、反距离权重插值、多项式插值、自然邻域法、样条函数法、趋势面法。

4.如权利要求1所述的一种农产品安全加密调查方法，其特征在于，所述概率梯度为20％时，从低到高将待预警区域划分为0-20％、20-40％、40-60％、60-80％、80-100％五个概率区域，所述五个区域分别对应名称为低概率区域、较低概率区域、中等概率区域、较高概率区域、高概率区域。