CN107918837A - 一种果蔬类食品安全风险预测方法 - Google Patents

Abstract

本发明涉及一种果蔬类食品安全风险预测方法，其主要技术特点是：对食品供应链的各个环节进行分析，找出各个环节的关键控制点；分析各环节关键控制点并将其作为量化指标，建立预测模型；初始化预测模型参数，根据预测模型的可见状态序列及系统的真实状态，对预测模型参数进行训练；利用预测模型对食品供应链风险等级进行风险评估和风险值计算。本发明利用预测模型对食品质量安全风险进行评价，考虑各个阶段的动态性，能够实时反应风险状态，采用定量和定性结合的评价方法，可以更精确地描述供应链风险的大小，有助于决策者及时采取应对措施，同时本发明具有很强的扩展性，可以应用到食品安全的其他领域。

Description

一种果蔬类食品安全风险预测方法

技术领域

本发明属于食品安全与计算机数据建模技术领域，尤其是一种果蔬类食品安全风险预测方法。

背景技术

频发的食品安全事件不仅严重威胁群众的身体健康，而且给社会带来了不安定因素，引发信任危机，影响经济发展。近年来食品安全问题的预防和治理受到高度关注，其中，风险评估是风险预警和治理的基础，受到格外重视。风险评估就是确认安全风险及其大小的过程，即利用适当的风险评估工具和方法，确定资产风险等级和优先控制顺序。食品安全风险评估的目的是为了全面掌握供应链风险的整体状况，从而为实施有效的风险管理措施提供决策支持。

目前常用的风险评估方法主要有贝叶斯网络法、德尔菲法、灰色理论法、SVM、人工神经网络、主成分分析法和层次分析法，下面对这些风险评估方法分别进行说明：

贝叶斯网络法是一种基于概率推理的图形化概率网络，能在有限的不确定信息条件下进行学习和推理，在智能化系统中应用广泛；但是贝叶斯网络构造繁琐，实际应用时还需反复交叉不断完善，易用性不够好。

德尔菲法是典型的定性分析方法，它能够充分发挥各位专家的不同意见，评定过程公正；但是实施过程比较复杂，花费时间较长，而且受专家主观因素影响较大，使得结果的精确度不够。

灰色系统是处于白色系统和黑箱系统之间的中间系统，是对不确定系统的研究方法；由于这种方法的指标权重和分辨系数需要人为设定，使评估结果人为干预性较大。

SVM在小样本的情况下具有较好的分类和泛化能力，但是不适合对于大样本数据的分析处理，与此同时，经典的SVM只给出了二类分类算法，解决多分类问题存在困难。

人工神经网络虽然具有高度的自适应能力，但是学习速度较慢，而且算法陷入局部极值的可能性较大。

层次分析法是一种将定量和定性相结合的多目标决策方法，是在风险分析或决策评价过程中应用最为广泛的方法；然而层次分析法存在判断矩阵一致性检验困难，缺乏科学性、未考虑评估者判断具有模糊性等缺点,同时，也不能根据系统状态的变化实时改变评价结果。

发明内容

本发明的目的在于克服现有技术的不足，提出一种果蔬类食品安全风险预测方法，解决现有风险评估方法缺乏预测性、实时性差等问题。

本发明解决其技术问题是采取以下技术方案实现的：

一种果蔬类食品安全风险预测方法，包括以下步骤：

步骤1、对食品供应链的各个环节进行分析，找出各个环节的关键控制点；

步骤2、分析各环节关键控制点并将其作为量化指标，建立预测模型；

步骤3、初始化预测模型参数，根据预测模型的可见状态序列及系统的真实状态，对预测模型参数进行训练；

步骤4、利用预测模型对食品供应链风险等级进行风险评估和风险值计算。

进一步，在步骤1中，找出各个环节的关键控制点的方法为：根据HACCP体系，并结合食品供应链的各环节，对各个环节中关键点的生物性、物理性和化学性的潜在危害进行分析，确定关键控制点。

进一步，在步骤2中，建立预测模型的方法包括以下步骤：

步骤201、确定食品供应链的各个系统可见状态；

步骤202、在食品供应链每个环节建立评估预测三元组参数模型，包括状态转移概率矩阵、观测向量的概率矩阵和最初状态分布矢量。

进一步，在步骤201中，各个系统可见状态包含系统的所有信息，并且在当前状态下的观察是独立的。

进一步，在步骤202中，每个环节的安全状态概率分布就是下一环节的初始状态概率分布。

进一步，在步骤3中，可见状态序列是指供应链各个环节检测到的数据；所述系统的真实状态为食品供应链各环节的真实风险值。

进一步，在步骤4中，风险值计算方法为：

其中，R_t表示在t时刻该环节所处的总体风险值，α_t(i)为t时刻供应链处在安全状态S_i的概率，N是安全状态的数目，c(i)是与状态S(i)关联的开销。

本发明的优点和积极效果是：

本发明利用预测模型对食品质量安全风险进行评价，考虑各个阶段的动态性，能够实时反应风险状态，采用定量和定性结合的评价方法，可以更精确地描述供应链风险的大小，有助于决策者及时采取应对措施。本发明的风险评估模型除了对果蔬类进行评估外，还具有很强的扩展性，可以应用到食品安全的其他领域。

附图说明

图1是本发明预测三大基本算法之间的关系示意图；

图2是本发明的预测模型的状态转移模型。

具体实施方式

以下结合附图对本发明实施例做进一步详述。

一种果蔬类食品安全风险预测方法，包括以下步骤：

步骤1、对食品供应链的各个环节进行分析，找出各个环节的关键控制点，其方法为：根据HACCP体系，并结合食品供应链的各环节，对各个环节中关键点的生物性、物理性和化学性的潜在危害进行分析，确定关键控制点。

步骤2、分析各环节关键控制点并将其作为量化指标，建立预测模型，其建立方法如下：

步骤201、确定食品供应链的各个系统可见状态，各个系统可见状态包含系统的所有信息，并且在当前状态下的观察是独立的。

步骤202、在食品供应链每个环节建立评估预测三元组参数模型，包括状态转移概率矩阵、观测向量的概率矩阵和最初状态分布矢量。

如图1所示，本步骤建立的预测模型为：

λ＝(N,M,π,P,Q) (1)

也可以简单表示为：

λ＝(N,M,π) (2)

其中：N为预测中状态的个数；M为预测中对应的观测值个数；π是初始状态的分布矢量，

π＝(π₁,π₂,…,π_N),

P为状态转移概率矩阵，表示如下：

Q是观察向量的概率矩阵，表示如下：

Q＝(q_jk)_N×M

1≤j≤N,1≤k≤M； (3)

Q_jk表示在状态j的情况下，观察状态k出现的概率，即：

Q_jk＝p(v_k|P_j)

1≤j≤N,1≤k≤M； (4)

在上述预测模型中，每个环节的安全状态概率分布就是下一环节的初始状态概率分布。

步骤3、初始化预测模型参数，根据预测模型的可见状态序列及系统的真实状态，对预测模型参数进行训练。其中，可见状态序列是指供应链各个环节检测到的数据，系统的真实状态为食品供应链各环节的真实风险值。

步骤4、利用预测模型对食品供应链风险等级进行风险评估和风险值计算，该风险值的计算公式如下：

式中，R_t表示在t时刻该环节所处的总体风险值，α_t(i)为t时刻供应链处在安全状态S_i的概率，N是安全状态的数目，c(i)是与状态S(i)关联的开销。

下面将预测模型应用在果蔬类领域作为本发明的实施例来进一步说明：

(1)对食品供应链的各个环节进行分析，找出各个环节的关键控制点

影响食品安全的危害因子主要分为三类：生物性因子、物理性因子和化学性危害因子。果蔬类在从农场到餐桌的过程经历了四个环节：种植、加工、储藏运输和消费。根据HACCP体系并结合果蔬类供应链的四个环节，对各个环节中的生物性、物理性和化学性的潜在危害进行分析，确定关键控制点。

由于食品一般都难以长久保持并容易受到微生物污染，食品供应链上任何环节的潜在危险因素都会影响食品质量安全。果蔬类加工环节是整个供应链最重要的部分，负责将收集到的原料加工完成，该过程对参与者的要求较高；储运环节涉及加工和物流企业，该环节受上一环节的影响巨大，同样的，消费环节的风险因素也基本来源供应链的前几个阶段。因此，生产供应链的各个环节紧密联系，要对每一个环节的检查监管都格外重视。最后一个环节的风险值即为整条供应链的风险大小。

(2)建立预测模型

食品安全风险评估会受很多因素的影响，假设影响因素为：X_N＝{X_N1,X_N2,…X_Ni}。其中X_N(1≤N≤4)为风险评估的观测值，X_Ni为第N个环节的第i个影响因素，通过一定的方法反映出风险评估值和其影响因素之间的关系，数学模型表达式为：

Y_N＝f(X_Ni), (5)

评估预测模型(预测)是序列数据处理和统计学习的一种重要概率模型，已经成功应用到风险评估领域。预测有两个随机过程：可见的状态序列(y_t；t＝1,2,…)，本文指供应链各个环节检测到的数据；系统的真实状态，本文指供应链各环节的真实风险值。供应链每个环节的评估预测三元组参数模型包括状态转移概率矩阵P，观测概率矩阵Q和最初状态分布π，字母表示为λ＝(P,Q,π)。除了第一环节——果蔬种植，每个环节的参数模型都受上一环节的影响。该环节最终的风险概率直接影响下一环节的初始状态概率分布。例如，果蔬种植环节通过预测模型计算，最终得出各个安全状态出现的概率分布为δ＝{δ₁,δ₂,…,δ_i}(i为真实状态的个数)，那么其下一环节即生产环节的初始状态分布概率π＝δ。

我们假设供应链共有五个状态，A₁-A₅，其中A₁表示正常安全状态，A₅表示重大危险状况，A₂、A₃、A₄表示危险等级逐级加深。如果用概率表示，则各个状态表示的供应链风险等级的概率分别为：

A₁：正常状态，设定供应链受威胁的概率为P(A)＝0；

A₂：低风险，设定供应链受威胁的概率为：

0＜P(A)＜0.2

A₃：中级风险，设定供应链受威胁的概率为0.2＜P(A)＜0.5；

A₄：中高级风险，设定供应链受威胁的概率为0.5＜P(A)＜0.8；

A₅：高风险，设定供应链受威胁的概率为0.8＜P(A)＜1。

这里所涉及的预测为离散预测，在建立离散的预测模型系统时，可以假定：系统状态包含系统的所有信息，并且在当前状态下的观察是独立的。则预测模型的状态转移模型可以表示为图2所示。从一个节点移动到另一个节点表示：系统在源节点显示的状态，会转变成目的节点的状态，该模型图是一个完全连接图，表明任意安全状态都有转变为其他任意安全状态的可能.

(3)实验

根据上述原理分析，为了简化实验过程，以供应链的最后一个环节为例，利用预测模型供应链进行风险评估。由于供应链各个环节紧密相连，环环相接，因此，最后一个环节的风险即可反应整条供应链的风险等级。

假设系统的可见状态有4个：V₁-V₄；隐藏状态我们假设有五个：A₁-A₅。上一环节的各个安全状态出现的概率分布为δ＝{0.56,0.12,0.19,0.05,0.08}，这个概率分布就是消费环节的初始状态概率分布。

设定状态转移矩阵P如下：

观测向量概率矩阵Q可以设定如下：

有了预测的三元组模型，就可以对相关数据进行解码分析。在此给定一组观察序列{O₂,O₁,O₃,O₄,O₄,O₃,O₂,O₂,O₁,O₁,O₄,O₃,}，通过维特比算法，可以评估出系统对应的真实状态{T₂,T₁,T₃,T₄,T₄,T₃,T₂,T₂,T₁,T₁,T₄,T₃,}。

由结果可知，系统在t＝5,6,9,11时风险级别较大，在t＝1,3,7时刻，风险级别较小。假设在A₁—A₅不同状态值下供应链对应的开销为C＝{0,5,10,15,20}。根据维特比算法，系统在t时刻处在各状态的概率分别为:

δ_t＝{0.1547,0.6073,0.0645,0.1043,0.0692}。

参考风险评价模型的典型代表OCTAVE给出一个计算风险值的公式

R_t＝∑α_t(i)*c(i)。 (6)

其中R_t表示在t时刻该环节所处的总体风险值，α_t(i)为t时刻供应链处在安全状态S_i的概率，N是安全状态的数目，c(i)是与状态S(i)关联的开销。

根据式(1)，我们可以计算出该供应链的风险值为：R＝6.0495。

可以利用上述方法对不同时刻进行风险评估与开销分析，有助于生产者和决策者更全面了解供应链的风险态势，及时查出问题发生的缘由，并根据已有流程或者规定采取相应措施。若发现供应链存在的新的问题，就需要相关人员对问题所在进行深入研究，完善风险规避的方法并建立解决问题新的流程。

需要强调的是，本发明所述的实施例是说明性的，而不是限定性的，因此本发明包括并不限于具体实施方式中所述的实施例，凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式，同样属于本发明保护的范围。

Claims (7)

1.一种果蔬类食品安全风险预测方法，其特征在于包括以下步骤：

步骤1、对食品供应链的各个环节进行分析，找出各个环节的关键控制点；

步骤2、分析各环节关键控制点并将其作为量化指标，建立预测模型；

步骤3、初始化预测模型参数，根据预测模型的可见状态序列及系统的真实状态，对预测模型参数进行训练；

步骤4、利用预测模型对食品供应链风险等级进行风险评估和风险值计算。

2.根据权利要求1所述的一种果蔬类食品安全风险预测方法一，其特征在于：在步骤1中，找出各个环节的关键控制点的方法为：根据HACCP体系，并结合食品供应链的各环节，对各个环节中关键点的生物性、物理性和化学性的潜在危害进行分析，确定关键控制点。

3.根据权利要求1或2所述的一种果蔬类食品安全风险预测方法，其特征在于：在步骤2中，建立预测模型的方法包括以下步骤：

步骤201、确定食品供应链的各个系统可见状态；

步骤202、在食品供应链每个环节建立评估预测三元组参数模型，包括状态转移概率矩阵、观测向量的概率矩阵和最初状态分布矢量。

4.根据权利要求3所述的一种果蔬类食品安全风险预测方法，其特征在于：在步骤201中，各个系统可见状态包含系统的所有信息，并且在当前状态下的观察是独立的。

5.根据权利要求3所述的一种果蔬类食品安全风险预测方法，其特征在于：在步骤202中，每个环节的安全状态概率分布就是下一环节的初始状态概率分布。

6.根据权利要求1所述的一种果蔬类食品安全风险预测方法，其特征在于：在步骤3中，可见状态序列是指供应链各个环节检测到的数据；所述系统的真实状态为食品供应链各环节的真实风险值。

7.根据权利要求1所述的一种果蔬类食品安全风险预测方法，其特征在于：在步骤4中，风险值计算方法为：

<mrow> <msub> <mi>R</mi> <mi>t</mi> </msub> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mi>i</mi> <mi>N</mi> </munderover> <msub> <mi>&amp;alpha;</mi> <mi>t</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>*</mo> <mi>c</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> </mrow>

其中，R_t表示在t时刻该环节所处的总体风险值，α_t(i)为t时刻供应链处在安全状态S_i的概率，N是安全状态的数目，c(i)是与状态S(i)关联的开销。