一种进口牛肉产地混合溯源方法及系统与流程

本发明涉及产地溯源领域，具体涉及一种进口牛肉产地混合溯源方法及系统。

背景技术：

随着社会和经济的发展，经济全球化，运输条件愈加便利。食品市场的全球化使消费者更加关注食品的产地，与此同时，牛肉作为高附加值的肉类，标明牛肉产地的标签标识错误、标签造假、冒用地理标志所导致的牛肉食品安全事件也常见报端，增加了公众对食品原产地认定的关注。我国《商标法》(修订)规定，原产地证明商标受法律保护，食品安全是牛肉市场规则不可或缺的一部分,每一个市场参与者都应该尊法守法，注意食品安全。同时重视牛肉的食品安全是维护牛肉市场的正常秩序和公平竞争的需要、是维护消费者的利益的要求，是实现牛肉商品价值的需要。

除此之外，疯牛病、牛口蹄疫等牛传染病的出现也会严重威胁牛肉产品的安全。当这些疫情爆发时，迅速准确地找到疫情源头并采取有效措施控制疫情的传播，有利于确保牛肉的食品安全、也有利于保障人们的身体健康。我国政府是人民意旨的执行者和人民利益的捍卫者，其宗旨是为人民服务，其原则是对人民负责，政府十分重视食品安全问题。牛肉作为市场流通性大的肉类，其食品安全问题关系到社会的各个方面，必须引起重视，防止出现食品安全事故。

目前，对牛肉的溯源多利用稳定同位素、近红外、电子舌和矿物元素等生物技术，以找出牛肉原产地，排查可能会出现病害的牛肉，预防食源性致病菌的传播。

公开号为cn109187130a的发明专利公开了一种基于矿物元素判别模型的牦牛肉溯源方法，其特征在于，包括：(1)采集牦牛肉样品，预处理，得到脱脂牦牛肉样品；(2)检测脱脂牦牛肉样品中矿物元素se、rb、ti含量；(3)将上述矿物元素含量代入判别模型，根据计算结果，判断出牦牛肉样品的产地。

上述溯源方法虽然能对牛肉进行溯源，但其主要针对国内牛肉进行溯源，如牛肉出自河南县、大通县、或碌曲县等，而不能对进口牛肉进行溯源。随着经济社会的发展，越来越多的牛肉依靠进口，且进口牛肉覆盖范围广，当出现疯牛病等时，其影响更大。且进口牛肉的特性与国内牛肉完全不同，因此现有的针对国内牛肉的溯源方法不适用于进口牛肉。此外，现有的牛肉溯源方法都是采用单一的生物技术进行溯源，溯源精度低。因此，如何为进口牛肉提供有效的溯源是本领域亟待解决的问题。

技术实现要素：

本发明的目的是针对现有技术的缺陷，提供了一种进口牛肉产地混合溯源方法及系统。生成基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型，综合四个鉴别模型的结果及排序位置判断进口牛肉的产地。本发明针对进口牛肉运输时间长、运输环境复杂等问题，提出了进口牛肉产地混合溯源的方法和系统，提高了进口牛肉的溯源精度，同时适应了现有的进口牛肉使用范围广的需求。

为了实现以上目的，本发明采用以下技术方案：

一种进口牛肉产地混合溯源方法，包括步骤：

s1、采集不同国家的牛肉样品，对牛肉样品进行粉碎、干燥、脱脂处理，脱脂完成后，对牛肉样本进行研磨，利用固定大小的筛板进行过滤，再次进行干燥，得到标准的进口牛肉样品；

s2、利用红外光谱仪采集所述标准的进口牛肉样品近红外光谱，所采集扫描的波数范围是8000-5000cm^-1，分辨率为4cm^-1，扫描的温度维持在25℃、湿度控制在保持稳定，每个样本的光谱扫描三次，对三次所采集的光谱数据求取平均值得到每个进口牛肉样本的光谱数据；根据牛肉样品近红外光谱的方差及标准差，剔除奇异点样本，得到最终的进口牛肉样品近红外光谱数据；

s3、利用同位素比率质谱仪检测所述标准的进口牛肉样品中δ¹³c、δ²h、δ¹⁵n的值，每个样本利用同位素比率质谱仪检测三次，对三次所采集的同位素质谱数据求取平均值得到每个进口牛肉样本的同位素质谱数据；根据牛肉样品同位素质谱数据的方差及标准差，剔除奇异点样本，得到最终的进口牛肉样品同位素质谱数据；

s4、利用等离子质谱仪检测所述标准的进口牛肉样品中se、rb和ti的含量，每个样本利用等离子质谱仪检测三次，对三次所采集的矿物元素数据求取平均值得到每个进口牛肉样本的矿物元素数据；根据牛肉样品同矿物元素数据的方差及标准差，剔除奇异点样本，得到最终的进口牛肉样品矿物元素数据；

s5、利用氨基酸分析仪检测所述标准的进口牛肉样品中asp、thr、ser、glu、gly、val、ile、leu、tyr、phe、lys、his12种氨基酸含量，每个样本利用氨基酸分析仪检测三次，对三次所采集的氨基酸含量求取平均值得到每个进口牛肉样本的氨基酸含量数据；根据牛肉样品同矿物元素数据的方差及标准差，剔除奇异点样本，得到最终的进口牛肉样品氨基酸数据；

s6、建立偏最小二乘法鉴别模型，对所述进口牛肉样品红外光谱数据、同位素质谱数据、矿物元素数据及氨基酸数据进行分割，选取1/n个样本数据作为测试集，其余样本数据为训练集，分别基于红外光谱数据、同位素质谱数据、矿物元素数据及氨基酸数据，不断对鉴别模型进行训练，得到基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型；

s7、交叉验证基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型的性能；

s8、对所述基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型按照性能进行排序，根据性能排名赋予鉴别模型相应的权重；

s9、分别利用基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型对进口牛肉进行检测，当判断进口牛肉不为相应国家的牛肉时，检测结果为0，否则检测结果为1；

s10、计算利用基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型得到的检测结果分别与相应的鉴别模型的加权和，将所述加权和与设定阈值进行对比，判断进口牛肉的产地。

进一步地，所述对牛肉样品进行粉碎、干燥具体为：

将进口牛肉样品切成丁，放入槽型混合机中在避光的环境中大约粉碎两个小时；牛肉粉碎后，将牛肉放入干燥室彻底干燥24小时。

进一步地，所述对牛肉样品进行粉碎、干燥具体为：所述脱脂处理具体为：采用静置法对进口牛肉样品的进行脱脂。

进一步地，所述剔除奇异点样本具体为：

计算每个产地进口牛肉样品的检测数据含量的平均值，当样本的方差和标准差超过设定阈值，将该样本从样品集中剔除；所述检测数据包括红外光谱数据、同位素质谱数据、矿物元素数据及氨基酸数据。

进一步地，所述步骤s7具体为：

将样本数据集随机划分为k个子集，一个子集作为验证集，其余的k-1组子集作为训练集；轮流将k个子集作为验证集，交叉重复k次，得到k次的结果，并用k次结果的平均值作为分类器或模型的性能指标。

进一步地，所述根据性能排名赋予鉴别模型相应的权重具体为：

假设基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型的权重分别为ω1、ω2、ω3、ω4，则：

ω1+ω2+ω3+ω4＝1

其中，基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型的性能越好，权重越大。

进一步地，所述步骤s10具体为：

通过基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型进行混合检测，假设其分别的检测结果为r1、r2、r3、r4，则鉴别模型最终的检测结果为：

r＝r1*ω1+r2*ω2+r3*ω3+r4*ω4

当r大于设定阈值，则进口牛肉为相应的国家进口，否则不为相应的国家进口。

本发明还提出一种进口牛肉产地混合溯源系统，包括：

预处理模块，用于采集不同国家的牛肉样品，对牛肉样品进行粉碎、干燥、脱脂处理，脱脂完成后，对牛肉样本进行研磨，利用固定大小的筛板进行过滤，再次进行干燥，得到标准的进口牛肉样品；

近红外光谱采集模块，用于利用红外光谱仪采集所述标准的进口牛肉样品近红外光谱，所采集扫描的波数范围是8000-5000cm^-1，分辨率为4cm^-1，扫描的温度维持在25℃、湿度控制在保持稳定，每个样本的光谱扫描三次，对三次所采集的光谱数据求取平均值得到每个进口牛肉样本的光谱数据；根据牛肉样品近红外光谱的方差及标准差，剔除奇异点样本，得到最终的进口牛肉样品近红外光谱数据；

同位素采集模块，用于利用同位素比率质谱仪检测所述标准的进口牛肉样品中δ¹³c、δ²h、δ¹⁵n的值，每个样本利用同位素比率质谱仪检测三次，对三次所采集的同位素质谱数据求取平均值得到每个进口牛肉样本的同位素质谱数据；根据牛肉样品同位素质谱数据的方差及标准差，剔除奇异点样本，得到最终的进口牛肉样品同位素质谱数据；

矿物元素采集模块，用于利用等离子质谱仪检测所述标准的进口牛肉样品中se、rb和ti的含量，每个样本利用等离子质谱仪检测三次，对三次所采集的矿物元素数据求取平均值得到每个进口牛肉样本的矿物元素数据；根据牛肉样品同矿物元素数据的方差及标准差，剔除奇异点样本，得到最终的进口牛肉样品矿物元素数据；

氨基酸采集模块，用于利用氨基酸分析仪检测所述标准的进口牛肉样品中asp、thr、ser、glu、gly、val、ile、leu、tyr、phe、lys、his12种氨基酸含量，每个样本利用氨基酸分析仪检测三次，对三次所采集的氨基酸含量求取平均值得到每个进口牛肉样本的氨基酸含量数据；根据牛肉样品同矿物元素数据的方差及标准差，剔除奇异点样本，得到最终的进口牛肉样品氨基酸数据；

模型训练模型，建立偏最小二乘法鉴别模型，对所述进口牛肉样品红外光谱数据、同位素质谱数据、矿物元素数据及氨基酸数据进行分割，选取1/n个样本数据作为测试集，其余样本数据为训练集，分别基于红外光谱数据、同位素质谱数据、矿物元素数据及氨基酸数据，不断对鉴别模型进行训练，得到基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型；

验证评价模块，用于交叉验证基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型的性能；

排序模块，用于对所述基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型按照性能进行排序，根据性能排名赋予鉴别模型相应的权重；

检测模块，用于分别利用基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型对进口牛肉进行检测，当判断进口牛肉不为相应国家的牛肉时，检测结果为0，否则检测结果为1；

综合溯源模块，用于计算利用基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型得到的检测结果分别与相应的鉴别模型的加权和，将所述加权和与设定阈值进行对比，判断进口牛肉的产地。

进一步地，所述根据性能排名赋予鉴别模型相应的权重具体为：

假设基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型的权重分别为ω1、ω2、ω3、ω4，则：

ω1+ω2+ω3+ω4＝1

其中，基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型的性能越好，权重越大。

进一步地，所述综合溯源模块包括：

通过基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型进行混合检测，假设其分别的检测结果为r1、r2、r3、r4，则鉴别模型最终的检测结果为：

r＝r1*ω1+r2*ω2+r3*ω3+r4*ω4

当r大于设定阈值，则进口牛肉为相应的国家进口，否则不为相应的国家进口。

本发明与现有技术相比，具有如下效果：

(1)本发明深入分析进口牛肉的特点，针对性地提出的进口牛肉产地混合溯源方法及系统，能够针对进口牛肉进行检测，适用范围光，具有极强的应用价值；

(2)由于进口牛肉在进口过程中，所花费的时间长，且运输的过程中环境差异大，现有的溯源技术中稳定同位素技术、矿物元素、氨基酸等存在受测试环境影响大、近红外光谱技术比较依赖数据库，需大量样品的来源和数量才能一定程度上提高溯源的精度等问题，如果直接将现有的溯源技术应用于进口牛肉，溯源的精度低。本发明结合红外光谱数据、同位素质谱数据、矿物元素数据、氨基酸混合对进口牛肉的产地进行溯源，通过四种溯源方式相结合，克服进口牛肉在运输过程中的影响，提高溯源的精度；

(3)本发明根据性能对四种溯源方式进行排序，并根据排序结果赋予不同溯源方式相应的的权重，能够根据溯源模型的性能调整其对综合溯源结果的影响，充分发挥各溯源方式的优势，而不是将各溯源方式的简单混合；

(4)本发明先对进口牛肉粉碎后再进行干燥，能够提高牛肉的干燥效果，同时采用静置法对进口牛肉样品的进行脱脂，能够快速制得大量的脱脂样品，操作简单，操作效率高；

(5)本发明通过多次采集每个进口牛肉样本的光谱数据、同位素数据、矿物元素数据、氨基酸数据，能够避免单次光谱数据采集过程中的误差；

(6)本发明对采集的每个进口牛肉样本的光谱数据、同位素数据、矿物元素数据、氨基酸数据，计算相应的方差和标准差，剔除相应的奇异点样本，避免采集环境等对数据的影响。

附图说明

图1是实施例一提供的一种进口牛肉产地混合溯源方法流程图；

图2是实施例二提供的一种进口牛肉产地混合溯源系统结构图。

具体实施方式

以下通过特定的具体实例说明本发明的实施方式，本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用，本说明书中的各项细节也可以基于不同观点与应用，在没有背离本发明的精神下进行各种修饰或改变。需说明的是，在不冲突的情况下，以下实施例及实施例中的特征可以相互组合。

需要说明的是，以下实施例中所提供的图示仅以示意方式说明本发明的基本构想，遂图式中仅显示与本发明中有关的组件而非按照实际实施时的组件数目、形状及尺寸绘制，其实际实施时各组件的型态、数量及比例可为一种随意的改变，且其组件布局型态也可能更为复杂。

下面结合附图和具体实施例对本发明作进一步说明，但不作为本发明的限定。

实施例一

如图1所示，本实施例提出了一种进口牛肉产地混合溯源方法，包括：

s1、采集不同国家的牛肉样品，对牛肉样品进行粉碎、干燥、脱脂处理，脱脂完成后，对牛肉样本进行研磨，利用固定大小的筛板进行过滤，再次进行干燥，得到标准的进口牛肉样品；

目前我国仅允许澳大利亚、新西兰、乌拉圭、阿根廷、加拿大、哥斯达黎加、智利、匈牙利等8个国家的牛肉进口。美国、巴西、日本等其他国家的牛肉被禁止进口销售。如果是名单之外的肉制品，可能经非法途径入境，或者是假冒的。因此，本发明对澳大利亚、新西兰、乌拉圭、阿根廷、加拿大、哥斯达黎加、智利、匈牙利、美国、巴西、日本这11个国家的牛肉进行采集，每个国家的牛肉采集200份，每份500g。

为了更准确地溯源进口牛肉，需要对牛肉样品进行预处理。首先将样品切成丁，放入槽型混合机中进行粉碎，在避光的环境中大约粉碎两个小时。牛肉粉碎后，进行干燥处理，可以将牛肉放入干燥室彻底干燥24小时。本发明先对进口牛肉粉碎后再进行干燥，能够提高牛肉的干燥效果。

现有技术中，进行进口牛肉脱脂主要是在样品中加入石油醚，利用索氏提取器进行提取，然后回收提取液中的石油醚,得到粗脂肪样品。但索氏提取器每次制备的脱脂样品量很少，且操作繁琐。本发明为了得到大量的脱脂样品，用于以后的检测分析，采用静置法对进口牛肉样品的进行脱脂。

s2、利用红外光谱仪采集所述标准的进口牛肉样品近红外光谱，所采集扫描的波数范围是8000-5000cm^-1，分辨率为4cm^-1，扫描的温度维持在25℃、湿度控制在保持稳定，每个样本的光谱扫描三次，对三次所采集的光谱数据求取平均值得到每个进口牛肉样本的光谱数据；根据牛肉样品近红外光谱的方差及标准差，剔除奇异点样本，得到最终的进口牛肉样品近红外光谱数据；

近红外光是波长在可见光与中红外光间的一段电磁波，其波长范围为780～2526nm。近红外光谱能反映红外射线与进口牛肉样品中物质的相互作用，在自然界中，每个分子的组成和结构是不同的，其内部不同官能团或化学键如o-h、c-h、n-h和s-h基团对红色外光的吸收作用均不同，这些氢基团再吸收近红外光中的部分能量后，会受激发从而发生跃迁，所以在近红外光谱上所显现的曲线的位置和吸收强度也不同，从而在近红外光谱仪中产生近红外光谱，并得到每个样品特征性的近红外吸收光谱。

近红外光谱反映样品中有机物的组分和含量的信息。不同产地的进口牛肉样品受品种、产地环境、加工方法、运输方式等因素影响,也受所处的生长环境、气候、土壤、水质等因素影响，导致食品的主要化学成分(如蛋白质、脂肪、水分等)的结构和含量存在一定差异，其有机物成分具有显著差异。除此之外，由于基因型、饲料类型、饲养方式、个体代谢等不同，各产地的动物机体中蛋白质、脂肪、水分等物质成分并不相同。在近红外光谱上，这些差异得到反映并由此得到每个产地的样品特征性的近红外吸收光谱。

本发明利用红外光谱仪采集所述标准的牛肉样品近红外光谱，取标准的进口牛肉样品中1～3g进行光谱采集，所采集扫描的波数范围是8000-5000cm^-1，分辨率为4cm^-1，扫描的温度维持在25℃、湿度控制在保持稳定，每个样本的光谱扫描三次，对三次所采集的光谱数据求取平均值得到每个进口牛肉样本的光谱数据。本发明通过多次采集每个进口牛肉样本的光谱数据，能够避免单次光谱数据采集过程中的误差。

此外，由于进口牛肉受环境影响较大，因此，不同的扫描环境对光谱数据的影响不同，因此，本发明在对不同的进口牛肉样本进行多次扫描时，其扫描的环境维持恒定。

在进口牛肉样本制作的过程中，可能存在被污染或者近红外光谱在采集的过程中受设备、环境等因素的影响，由此导致的进口牛肉样本数据不准确的问题。因此，本发明在利用样本数据进行进口牛肉产地溯源前，对样本过程中出现的奇异点样本进行剔除，避免奇异点样本对溯源结果的影响。

本发明采用牛肉样品近红外光谱的方差及标准差，剔除奇异点样本。同一产地的进口牛肉其近红外光谱通常表现出相似的特征，不同进口牛肉样本采集的光谱间差异较小。因此，本发明利用不同牛肉样品近红外光谱的方差及标准差进行比较，当某一近红外光谱的方差及标准差超过一定的阈值时，说明该样本偏离同一产地其它样本的幅度大，该样本很可能为奇异点样本，进行剔除。

因此，本发明计算每个产地进口牛肉200份样品的红外光谱平均值，进一步计算每个进口牛肉样品的方差及标准差。当该样品的方差及标准差超过设定阈值，将该样品从样品集中剔除。

s3、利用同位素比率质谱仪检测所述标准的进口牛肉样品中δ¹³c、δ²h、δ¹⁵n的值，每个样本利用同位素比率质谱仪检测三次，对三次所采集的同位素质谱数据求取平均值得到每个进口牛肉样本的同位素质谱数据；根据牛肉样品同位素质谱数据的方差及标准差，剔除奇异点样本，得到最终的进口牛肉样品同位素质谱数据；

稳定同位素技术的进行产地溯源的原理是利用同位素的自然分馏效应来鉴定目标样品的产地。因气温、日照、土壤、食料、空气质量等差异，致使目标样品中某元素的同位素丰度与不同产地的其他样品存在显著差异，以此可以比较准确的区分并进行动物源性农产品的产地溯源。具体来说，通过稳定同位素中¹³c和¹²c的比例可以表征饲料种类，其比值与饲料中c3、c4植物所占比例密切相关；而同位素中¹⁵n和¹⁴n的比例则受很多因素影响，主要取决于营养水平，不但与饲料种类密切相关，也可以指示土壤、气候、农业施肥中的差异，甚至与海洋和陆地植物在饲料中的占比有所联系；除此之外，同位素中¹⁸o和¹⁶o的比例和²h和¹h的比例与产地当地的气候、地形、水的蒸发、浓缩、沉降有关；同位素中³⁴s和³²s的比例与微生物作用和海洋因素有关。

由于不同地区来源的各国进口牛肉所喂养的饲料、饮用水、土壤、气候等完全不同。尤其是碳、氮、氢同位素呈极显著相关性，δ¹³c、δ²h、δ¹⁵n作为本发明区牛肉的主要指标，其中δ¹³c反映出饲料中c3和c4植物的比例，而δ²h组成与饮用水有关，δ¹⁵n反应进口牛肉所在国土壤条件。

检测所述标准的进口牛肉样品中δ¹³c的值：将每个标准的进口牛肉样品取2～4g放入锡箔杯中，通过自动进样器进入元素分析仪，将进口牛肉样品燃烧，以将其转化成纯净的co2和n2，co2经过稀释器的稀释，最终将co2送入同位素比率质谱仪进行检测。其中，燃烧炉的温度为1200℃，还原炉的温度为600℃。

检测所述标准的进口牛肉样品中δ²h的值：每个标准的进口牛肉样品取2～4g放入银杯，平衡72小时，按顺序放入自动进样器中，通过自动进样器进入元素分析仪，将进口牛肉样品裂解成co和h2，最终将h2送入同位素比率质谱仪进行检测。其中，裂解的温度为1450℃。

检测所述标准的进口牛肉样品中δ¹⁵n的值：将每个标准的进口牛肉样品取2～4g放入锡箔杯中，通过自动进样器进入元素分析仪，将进口牛肉样品燃烧，以将其转化成纯净的co2和n2，最终将n2送入同位素比率质谱仪进行检测。其中，燃烧炉的温度为1200℃，还原炉的温度为600℃。

每个样本进口牛肉样品中δ¹³c、δ²h、δ¹⁵n的值均进行三次采集，对三次所采集的数据求取平均值得到每个进口牛肉样本中δ¹³c、δ²h、δ¹⁵n的值。本发明通过多次采集每个进口牛肉样本的中δ¹³c、δ²h、δ¹⁵n的值，能够避免单次数据采集过程中的误差。同位素比率质谱仪

与近红外光谱采集一样，δ¹³c、δ²h、δ¹⁵n的值也可能受环境、采集的条件等因素影响，导致的进口牛肉样本数据不准确的问题。因此，本发明在利用样本数据进行进口牛肉产地溯源前，对样本过程中出现的奇异点样本进行剔除，避免奇异点样本对溯源结果的影响。

本发明分别根据牛肉样品δ¹³c、δ²h、δ¹⁵n的方差及标准差，剔除奇异点样本。只要δ¹³c、δ²h、δ¹⁵n任意一个值的方差及标准差超过一定的阈值时，说明该样本偏离同一产地其它样本的幅度大，该样本很可能为奇异点样本，进行剔除。

因此，本发明计算每个产地进口牛肉200份样品的δ¹³c、δ²h、δ¹⁵n平均值，进一步计算每个进口牛肉样品的方差及标准差。当该样品的方差及标准差超过设定阈值，将该样品从样品集中剔除。

s4、利用等离子质谱仪检测所述标准的进口牛肉样品中se、rb和ti的含量，每个样本利用等离子质谱仪检测三次，对三次所采集的矿物元素数据求取平均值得到每个进口牛肉样本的矿物元素数据；根据牛肉样品同矿物元素数据的方差及标准差，剔除奇异点样本，得到最终的进口牛肉样品矿物元素数据；

牛肉中的各元素含量与产品产地的水源、土壤等条件息息相关，且不同的地区各自有其特征性的元素组成，从而能够建立具有不同产地进口牛肉的矿物元素指纹图谱，实现对进口牛肉的较为准确的产地溯源。矿物元素指纹图谱技术依据进口牛肉样品中常量元素(钙、磷、镁、钾、钠、氯和硫等)和矿物元素(铁、铜、锰、锌、碘、硒和铬等)的测定数值，通过统计分析并选取有显著差异的元素建立模型，从而能够更准确的鉴别进口牛肉的产地。

不同地域土壤中痕量元素的含量与组成有其典型特，例如美国土壤含量高于欧洲。土壤、植物及进口牛肉中se含量呈极显著相关关系。从高se区和低se区来的牛肉中含量有很大差异，因此本发明对进口牛肉中的se含量进行测定，以区别不同国家的进口牛肉。北美牛肉中se含量大约是瑞士牛肉的两倍，而巴西牛肉中se的含量更高。rb含量与岩石特性有关，一般花岗岩和片麻岩地区的土壤、水、植物中含量表现最高。牛肉对rb的累积程度高，因此，本发明还对进口牛肉中的rb含量进行测定，以区别不同国家的进口牛肉。不同喂养环境的牛肉其ti含量差异较大，因此，本发明还对进口牛肉中的ti含量进行测定，以区别不同国家的进口牛肉。

利用等离子质谱仪检测所述标准的进口牛肉样品中se、rb和ti的含量具体为：将每个标准的进口牛肉样品取2～4g，放入消解管中，通过浓硝酸预消解2小时，然后利用双氧水消毒1个小时，最后放入微波消解仪进行消解。将消解好的样品利用离子质谱仪对se、rb和ti进行采集，最后使用外标法对采集的se、rb和ti进行定量，测得se、rb和ti的含量。

每个样本进口牛肉样品中se、rb和ti的含量值均进行三次采集，对三次所采集的数据求取平均值得到每个进口牛肉样本中se、rb和ti的含量的值。本发明通过多次采集每个进口牛肉样本的中se、rb和ti的含量值，能够避免单次数据采集过程中的误差。

se、rb和ti的含量值也可能受环境、采集的条件等因素影响，导致的进口牛肉样本数据不准确的问题。因此，本发明在利用样本数据进行进口牛肉产地溯源前，对样本过程中出现的奇异点样本进行剔除，避免奇异点样本对溯源结果的影响。

本发明分别根据牛肉样品se、rb和ti的含量的方差及标准差，剔除奇异点样本。只要se、rb和ti的含量任意一个值的方差及标准差超过一定的阈值时，说明该样本偏离同一产地其它样本的幅度大，该样本很可能为奇异点样本，进行剔除。

因此，本发明计算每个产地进口牛肉200份样品的se、rb和ti的含量平均值，进一步计算每个进口牛肉样品的方差及标准差。当该样品的方差及标准差超过设定阈值，将该样品从样品集中剔除。

s5、利用氨基酸分析仪检测所述标准的进口牛肉样品中asp、thr、ser、glu、gly、val、ile、leu、tyr、phe、lys、his12种氨基酸含量，每个样本利用氨基酸分析仪检测三次，对三次所采集的氨基酸含量求取平均值得到每个进口牛肉样本的氨基酸含量数据；根据牛肉样品同矿物元素数据的方差及标准差，剔除奇异点样本，得到最终的进口牛肉样品氨基酸数据；

牛肉氨基酸含量与品种、性别、年龄、肌肉部位及养殖环境等因素有关，而对相同的品种和年龄，肌肉氨基酸含量受到环境影响十分显著。产地鉴别分析主要是探寻表征不同地域来源牛肉的特异性指标。因此，本发明引入氨基酸含量的测定，以区分不同国家的进口牛肉。asp、thr、ser、glu、gly、val、ile、leu、tyr、phe、lys、his12种氨基酸是牛肉中主要的氨基酸种类，因此，为了全面分析各国进口牛肉的差异，本发明对asp、thr、ser、glu、gly、val、ile、leu、tyr、phe、lys、his12种氨基酸进行检测。

具体地，本发明采用l-8800氨基酸分析仪进行氨基酸测定。首先将所述标准的进口牛肉样本取2～4g放入水解管中，加入6mol/l的盐酸50ml，将水解管抽真空、密封，将其在110℃的恒温环境中水解24小时，并在冷却后进行过滤。调节滤液ph值至中性，定容值125ml，取滤液与0.02mol/l的盐酸按1:1进行混合，将混合液通过微孔滤膜进行过滤，然后采用l-8800氨基酸分析仪测定混合液中asp、thr、ser、glu、gly、val、ile、leu、tyr、phe、lys、his12种氨基酸的含量。

氨基酸的含量值也可能受环境、采集的条件等因素影响，导致的进口牛肉样本数据不准确的问题。因此，本发明在利用样本数据进行进口牛肉产地溯源前，对样本过程中出现的奇异点样本进行剔除，避免奇异点样本对溯源结果的影响。

本发明分别根据牛肉样品asp、thr、ser、glu、gly、val、ile、leu、tyr、phe、lys、his12种氨基酸含量的方差及标准差，剔除奇异点样本。只要asp、thr、ser、glu、gly、val、ile、leu、tyr、phe、lys、his12种氨基酸含量任意一个值的方差及标准差超过一定的阈值时，说明该样本偏离同一产地其它样本的幅度大，该样本很可能为奇异点样本，进行剔除。

因此，本发明计算每个产地进口牛肉200份样品的asp、thr、ser、glu、gly、val、ile、leu、tyr、phe、lys、his12种氨基酸含量平均值，进一步计算每个进口牛肉样品的方差及标准差。当该样品的方差及标准差超过设定阈值，将该样品从样品集中剔除。

值得注意的是，步骤s2、s3、s4、s5的先后顺序不作限定，对进口牛肉样品红外光谱数据、同位素质谱数据、矿物元素数据及氨基酸数据可以根据任意顺序进行先后采集，也可以同时进行采集。

s6、建立偏最小二乘法鉴别模型，对所述进口牛肉样品红外光谱数据、同位素质谱数据、矿物元素数据及氨基酸数据进行分割，选取1/n个样本数据作为测试集，其余样本数据为训练集，分别基于红外光谱数据、同位素质谱数据、矿物元素数据及氨基酸数据，不断对鉴别模型进行训练，得到基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型；

由于进口牛肉在进口过程中，所花费的时间长，且运输的过程中环境差异大，现有的溯源技术存在稳定同位素技术、矿物元素、氨基酸等受测试环境影响大、近红外光谱技术比较依赖数据库，需大量样品的来源和数量才能一定程度上提高溯源的精度等问题，如果直接将现有的溯源技术应用于进口牛肉，溯源的精度低。因此，本发明结合红外光谱数据、同位素质谱数据、矿物元素数据、氨基酸混合对进口牛肉的产地进行溯源，通过四种溯源方式相结合，克服进口牛肉在运输过程中的影响，提高溯源的精度。

本发明在采集到进口牛肉样品红外光谱数据、同位素质谱数据、矿物元素数据及氨基酸数据后，分别对样品数据进行处理，也就是说，本发明基于进口牛肉样品红外光谱数据、同位素质谱数据、矿物元素数据及氨基酸数据分别训练得到相应的溯源模型。

s7、交叉验证基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型的性能；

本发明将个样本数据集随机划分为k个子集(一般是均分)，将一个子集作为验证集，其余的k-1组子集作为训练集；轮流将k个子集作为验证集，交叉重复k次，得到k次的结果，并用k次结果的平均值作为分类器或模型的性能指标。

也就是说，本发明轮流对红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型进行交叉验证，获取各模型的性能。本发明通过溯源的正确率来指示模型的性能。

s8、对所述基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型按照性能进行排序，根据性能排名赋予鉴别模型相应的权重；

本发明采用正确率来指示模型的性能，因此，正确率越高，鉴别模型的性能越好，溯源精度越高。本发明通过四种鉴别模型对进口牛肉进行综合溯源，因此，为不同的鉴别模型赋予不同的权重，权重代表该鉴别模型对最终溯源结果的影响。权重越高，其对溯源结果影响越大。因此，排名越高的鉴别模型其权重越大，假设基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型的权重分别为ω1、ω2、ω3、ω4，则：

ω1+ω2+ω3+ω4＝1

s9、分别利用基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型对进口牛肉进行检测，当判断进口牛肉不为相应国家的牛肉时，检测结果为0，否则检测结果为1；

取待测未知过度的牛肉按照上述步骤s1、s2、s3、s4，测定近红外特征光谱数据、稳定同位素质谱数据、矿物元素数据和氨基酸数据，将测得数据分别代入基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型，若预测结果为0，则判断待测进口牛肉不为相应国家的进口牛肉；若预测结果为1，则判断待测样品为相应国家的进口牛肉。

s10、计算利用基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型得到的检测结果分别与相应的鉴别模型的加权和，将所述加权和与设定阈值进行对比，判断进口牛肉的产地。

本发明通过基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型进行混合检测，假设其分别的检测结果为r1、r2、r3、r4，则鉴别模型最终的检测结果为：

r＝r1*ω1+r2*ω2+r3*ω3+r4*ω4

将最终的检测结果与设定的阈值进行比较，假设设置阈值为0.75，当计算的检测结果大于设定阈值时，表面进口牛肉为相应的国家进口，否则不为相应的国家进口。

实施例二

如图2所示，本实施例提出了一种进口牛肉产地混合溯源系统，包括：

预处理模块，用于采集不同国家的牛肉样品，对牛肉样品进行粉碎、干燥、脱脂处理，脱脂完成后，对牛肉样本进行研磨，利用固定大小的筛板进行过滤，再次进行干燥，得到标准的进口牛肉样品；

目前我国仅允许澳大利亚、新西兰、乌拉圭、阿根廷、加拿大、哥斯达黎加、智利、匈牙利等8个国家的牛肉进口。美国、巴西、日本等其他国家的牛肉被禁止进口销售。如果是名单之外的肉制品，可能经非法途径入境，或者是假冒的。因此，本发明对澳大利亚、新西兰、乌拉圭、阿根廷、加拿大、哥斯达黎加、智利、匈牙利、美国、巴西、日本这11个国家的牛肉进行采集，每个国家的牛肉采集200份，每份500g。

为了更准确地溯源进口牛肉，需要对牛肉样品进行预处理。首先将样品切成丁，放入槽型混合机中进行粉碎，在避光的环境中大约粉碎两个小时。牛肉粉碎后，进行干燥处理，可以将牛肉放入干燥室彻底干燥24小时。本发明先对进口牛肉粉碎后再进行干燥，能够提高牛肉的干燥效果。

现有技术中，进行进口牛肉脱脂主要是在样品中加入石油醚，利用索氏提取器进行提取，然后回收提取液中的石油醚,得到粗脂肪样品。但索氏提取器每次制备的脱脂样品量很少，且操作繁琐。本发明为了得到大量的脱脂样品，用于以后的检测分析，采用静置法对进口牛肉样品的进行脱脂。

近红外光谱采集模块，用于利用红外光谱仪采集所述标准的进口牛肉样品近红外光谱，所采集扫描的波数范围是8000-5000cm^-1，分辨率为4cm^-1，扫描的温度维持在25℃、湿度控制在保持稳定，每个样本的光谱扫描三次，对三次所采集的光谱数据求取平均值得到每个进口牛肉样本的光谱数据；根据牛肉样品近红外光谱的方差及标准差，剔除奇异点样本，得到最终的进口牛肉样品近红外光谱数据；

近红外光是波长在可见光与中红外光间的一段电磁波，其波长范围为780～2526nm。近红外光谱能反映红外射线与进口牛肉样品中物质的相互作用，在自然界中，每个分子的组成和结构是不同的，其内部不同官能团或化学键如o-h、c-h、n-h和s-h基团对红色外光的吸收作用均不同，这些氢基团再吸收近红外光中的部分能量后，会受激发从而发生跃迁，所以在近红外光谱上所显现的曲线的位置和吸收强度也不同，从而在近红外光谱仪中产生近红外光谱，并得到每个样品特征性的近红外吸收光谱。

近红外光谱反映样品中有机物的组分和含量的信息。不同产地的进口牛肉样品受品种、产地环境、加工方法、运输方式等因素影响,也受所处的生长环境、气候、土壤、水质等因素影响，导致食品的主要化学成分(如蛋白质、脂肪、水分等)的结构和含量存在一定差异，其有机物成分具有显著差异。除此之外，由于基因型、饲料类型、饲养方式、个体代谢等不同，各产地的动物机体中蛋白质、脂肪、水分等物质成分并不相同。在近红外光谱上，这些差异得到反映并由此得到每个产地的样品特征性的近红外吸收光谱。

本发明利用红外光谱仪采集所述标准的牛肉样品近红外光谱，取标准的进口牛肉样品中1～3g进行光谱采集，所采集扫描的波数范围是8000-5000cm^-1，分辨率为4cm^-1，扫描的温度维持在25℃、湿度控制在保持稳定，每个样本的光谱扫描三次，对三次所采集的光谱数据求取平均值得到每个进口牛肉样本的光谱数据。本发明通过多次采集每个进口牛肉样本的光谱数据，能够避免单次光谱数据采集过程中的误差。

此外，由于进口牛肉受环境影响较大，因此，不同的扫描环境对光谱数据的影响不同，因此，本发明在对不同的进口牛肉样本进行多次扫描时，其扫描的环境维持恒定。

在进口牛肉样本制作的过程中，可能存在被污染或者近红外光谱在采集的过程中受设备、环境等因素的影响，由此导致的进口牛肉样本数据不准确的问题。因此，本发明在利用样本数据进行进口牛肉产地溯源前，对样本过程中出现的奇异点样本进行剔除，避免奇异点样本对溯源结果的影响。

本发明采用牛肉样品近红外光谱的方差及标准差，剔除奇异点样本。同一产地的进口牛肉其近红外光谱通常表现出相似的特征，不同进口牛肉样本采集的光谱间差异较小。因此，本发明利用不同牛肉样品近红外光谱的方差及标准差进行比较，当某一近红外光谱的方差及标准差超过一定的阈值时，说明该样本偏离同一产地其它样本的幅度大，该样本很可能为奇异点样本，进行剔除。

因此，本发明计算每个产地进口牛肉200份样品的红外光谱平均值，进一步计算每个进口牛肉样品的方差及标准差。当该样品的方差及标准差超过设定阈值，将该样品从样品集中剔除。

同位素采集模块，用于利用同位素比率质谱仪检测所述标准的进口牛肉样品中δ¹³c、δ²h、δ¹⁵n的值，每个样本利用同位素比率质谱仪检测三次，对三次所采集的同位素质谱数据求取平均值得到每个进口牛肉样本的同位素质谱数据；根据牛肉样品同位素质谱数据的方差及标准差，剔除奇异点样本，得到最终的进口牛肉样品同位素质谱数据；

稳定同位素技术的进行产地溯源的原理是利用同位素的自然分馏效应来鉴定目标样品的产地。因气温、日照、土壤、食料、空气质量等差异，致使目标样品中某元素的同位素丰度与不同产地的其他样品存在显著差异，以此可以比较准确的区分并进行动物源性农产品的产地溯源。具体来说，通过稳定同位素中¹³c和¹²c的比例可以表征饲料种类，其比值与饲料中c3、c4植物所占比例密切相关；而同位素中¹⁵n和¹⁴n的比例则受很多因素影响，主要取决于营养水平，不但与饲料种类密切相关，也可以指示土壤、气候、农业施肥中的差异，甚至与海洋和陆地植物在饲料中的占比有所联系；除此之外，同位素中¹⁸o和¹⁶o的比例和²h和¹h的比例与产地当地的气候、地形、水的蒸发、浓缩、沉降有关；同位素中³⁴s和³²s的比例与微生物作用和海洋因素有关。

由于不同地区来源的各国进口牛肉所喂养的饲料、饮用水、土壤、气候等完全不同。尤其是碳、氮、氢同位素呈极显著相关性，δ¹³c、δ²h、δ¹⁵n作为本发明区牛肉的主要指标，其中δ¹³c反映出饲料中c3和c4植物的比例，而δ²h组成与饮用水有关，δ¹⁵n反应进口牛肉所在国土壤条件。

检测所述标准的进口牛肉样品中δ¹³c的值：将每个标准的进口牛肉样品取2～4g放入锡箔杯中，通过自动进样器进入元素分析仪，将进口牛肉样品燃烧，以将其转化成纯净的co2和n2，co2经过稀释器的稀释，最终将co2送入同位素比率质谱仪进行检测。其中，燃烧炉的温度为1200℃，还原炉的温度为600℃。

检测所述标准的进口牛肉样品中δ²h的值：每个标准的进口牛肉样品取2～4g放入银杯，平衡72小时，按顺序放入自动进样器中，通过自动进样器进入元素分析仪，将进口牛肉样品裂解成co和h2，最终将h2送入同位素比率质谱仪进行检测。其中，裂解的温度为1450℃。

检测所述标准的进口牛肉样品中δ¹⁵n的值：将每个标准的进口牛肉样品取2～4g放入锡箔杯中，通过自动进样器进入元素分析仪，将进口牛肉样品燃烧，以将其转化成纯净的co2和n2，最终将n2送入同位素比率质谱仪进行检测。其中，燃烧炉的温度为1200℃，还原炉的温度为600℃。

每个样本进口牛肉样品中δ¹³c、δ²h、δ¹⁵n的值均进行三次采集，对三次所采集的数据求取平均值得到每个进口牛肉样本中δ¹³c、δ²h、δ¹⁵n的值。本发明通过多次采集每个进口牛肉样本的中δ¹³c、δ²h、δ¹⁵n的值，能够避免单次数据采集过程中的误差。

与近红外光谱采集一样，δ¹³c、δ²h、δ¹⁵n的值也可能受环境、采集的条件等因素影响，导致的进口牛肉样本数据不准确的问题。因此，本发明在利用样本数据进行进口牛肉产地溯源前，对样本过程中出现的奇异点样本进行剔除，避免奇异点样本对溯源结果的影响。

本发明分别根据牛肉样品δ¹³c、δ²h、δ¹⁵n的方差及标准差，剔除奇异点样本。只要δ¹³c、δ²h、δ¹⁵n任意一个值的方差及标准差超过一定的阈值时，说明该样本偏离同一产地其它样本的幅度大，该样本很可能为奇异点样本，进行剔除。

因此，本发明计算每个产地进口牛肉200份样品的δ¹³c、δ²h、δ¹⁵n平均值，进一步计算每个进口牛肉样品的方差及标准差。当该样品的方差及标准差超过设定阈值，将该样品从样品集中剔除。

矿物元素采集模块，用于利用等离子质谱仪检测所述标准的进口牛肉样品中se、rb和ti的含量，每个样本利用等离子质谱仪检测三次，对三次所采集的矿物元素数据求取平均值得到每个进口牛肉样本的矿物元素数据；根据牛肉样品同矿物元素数据的方差及标准差，剔除奇异点样本，得到最终的进口牛肉样品矿物元素数据；

牛肉中的各元素含量与产品产地的水源、土壤等条件息息相关，且不同的地区各自有其特征性的元素组成，从而能够建立具有不同产地进口牛肉的矿物元素指纹图谱，实现对进口牛肉的较为准确的产地溯源。矿物元素指纹图谱技术依据进口牛肉样品中常量元素(钙、磷、镁、钾、钠、氯和硫等)和矿物元素(铁、铜、锰、锌、碘、硒和铬等)的测定数值，通过统计分析并选取有显著差异的元素建立模型，从而能够更准确的鉴别进口牛肉的产地。

不同地域土壤中痕量元素的含量与组成有其典型特，例如美国土壤含量高于欧洲。土壤、植物及进口牛肉中se含量呈极显著相关关系。从高se区和低se区来的牛肉中含量有很大差异，因此本发明对进口牛肉中的se含量进行测定，以区别不同国家的进口牛肉。北美牛肉中se含量大约是瑞士牛肉的两倍，而巴西牛肉中se的含量更高。rb含量与岩石特性有关，一般花岗岩和片麻岩地区的土壤、水、植物中含量表现最高。牛肉对rb的累积程度高，因此，本发明还对进口牛肉中的rb含量进行测定，以区别不同国家的进口牛肉。不同喂养环境的牛肉其ti含量差异较大，因此，本发明还对进口牛肉中的ti含量进行测定，以区别不同国家的进口牛肉。

利用等离子质谱仪检测所述标准的进口牛肉样品中se、rb和ti的含量具体为：将每个标准的进口牛肉样品取2～4g，放入消解管中，在恒温环境中通过浓硝酸预消解2小时，然后利用双氧水消毒1个小时，最后放入微波消解仪进行消解。将消解好的样品利用离子质谱仪对se、rb和ti进行采集，最后使用外标法对采集的se、rb和ti进行定量，测得se、rb和ti的含量。

每个样本进口牛肉样品中se、rb和ti的含量值均进行三次采集，对三次所采集的数据求取平均值得到每个进口牛肉样本中se、rb和ti的含量的值。本发明通过多次采集每个进口牛肉样本的中se、rb和ti的含量值，能够避免单次数据采集过程中的误差。

se、rb和ti的含量值也可能受环境、采集的条件等因素影响，导致的进口牛肉样本数据不准确的问题。因此，本发明在利用样本数据进行进口牛肉产地溯源前，对样本过程中出现的奇异点样本进行剔除，避免奇异点样本对溯源结果的影响。

本发明分别根据牛肉样品se、rb和ti的含量的方差及标准差，剔除奇异点样本。只要se、rb和ti的含量任意一个值的方差及标准差超过一定的阈值时，说明该样本偏离同一产地其它样本的幅度大，该样本很可能为奇异点样本，进行剔除。

因此，本发明计算每个产地进口牛肉200份样品的se、rb和ti的含量平均值，进一步计算每个进口牛肉样品的方差及标准差。当该样品的方差及标准差超过设定阈值，将该样品从样品集中剔除。

氨基酸采集模块，用于利用氨基酸分析仪检测所述标准的进口牛肉样品中asp、thr、ser、glu、gly、val、ile、leu、tyr、phe、lys、his12种氨基酸含量，每个样本利用氨基酸分析仪检测三次，对三次所采集的氨基酸含量求取平均值得到每个进口牛肉样本的氨基酸含量数据；根据牛肉样品同矿物元素数据的方差及标准差，剔除奇异点样本，得到最终的进口牛肉样品氨基酸数据；

牛肉氨基酸含量与品种、性别、年龄、肌肉部位及养殖环境等因素有关，而对相同的品种和年龄，肌肉氨基酸含量受到环境影响十分显著。产地鉴别分析主要是探寻表征不同地域来源牛肉的特异性指标。因此，本发明引入氨基酸含量的测定，以区分不同国家的进口牛肉。asp、thr、ser、glu、gly、val、ile、leu、tyr、phe、lys、his12种氨基酸是牛肉中主要的氨基酸种类，因此，为了全面分析各国进口牛肉的差异，本发明对asp、thr、ser、glu、gly、val、ile、leu、tyr、phe、lys、his12种氨基酸进行检测。

具体地，本发明采用l-8800氨基酸分析仪进行氨基酸测定。首先将所述标准的进口牛肉样本取2～4g放入水解管中，加入6mol/l的盐酸50ml，将水解管抽真空、密封，将其在110℃的恒温环境中水解24小时，并在冷却后进行过滤。调节滤液ph值至中性，定容值125ml，取滤液与0.02mol/l的盐酸按1:1进行混合，将混合液通过微孔滤膜进行过滤，然后采用l-8800氨基酸分析仪测定混合液中asp、thr、ser、glu、gly、val、ile、leu、tyr、phe、lys、his12种氨基酸的含量。

氨基酸的含量值也可能受环境、采集的条件等因素影响，导致的进口牛肉样本数据不准确的问题。因此，本发明在利用样本数据进行进口牛肉产地溯源前，对样本过程中出现的奇异点样本进行剔除，避免奇异点样本对溯源结果的影响。

本发明分别根据牛肉样品asp、thr、ser、glu、gly、val、ile、leu、tyr、phe、lys、his12种氨基酸含量的方差及标准差，剔除奇异点样本。只要asp、thr、ser、glu、gly、val、ile、leu、tyr、phe、lys、his12种氨基酸含量任意一个值的方差及标准差超过一定的阈值时，说明该样本偏离同一产地其它样本的幅度大，该样本很可能为奇异点样本，进行剔除。

因此，本发明计算每个产地进口牛肉200份样品的asp、thr、ser、glu、gly、val、ile、leu、tyr、phe、lys、his12种氨基酸含量平均值，进一步计算每个进口牛肉样品的方差及标准差。当该样品的方差及标准差超过设定阈值，将该样品从样品集中剔除。

模型训练模型，建立偏最小二乘法鉴别模型，对所述进口牛肉样品红外光谱数据、同位素质谱数据、矿物元素数据及氨基酸数据进行分割，选取1/n个样本数据作为测试集，其余样本数据为训练集，分别基于红外光谱数据、同位素质谱数据、矿物元素数据及氨基酸数据，不断对鉴别模型进行训练，得到基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型；

由于进口牛肉在进口过程中，所花费的时间长，且运输的过程中环境差异大，现有的溯源技术存在稳定同位素技术、矿物元素、氨基酸等受测试环境影响大、近红外光谱技术比较依赖数据库，需大量样品的来源和数量才能一定程度上提高溯源的精度等问题，如果直接将现有的溯源技术应用于进口牛肉，溯源的精度低。因此，本发明结合红外光谱数据、同位素质谱数据、矿物元素数据、氨基酸混合对进口牛肉的产地进行溯源，通过四种溯源方式相结合，克服进口牛肉在运输过程中的影响，提高溯源的精度。

本发明在采集到进口牛肉样品红外光谱数据、同位素质谱数据、矿物元素数据及氨基酸数据后，分别对样品数据进行处理，也就是说，本发明基于进口牛肉样品红外光谱数据、同位素质谱数据、矿物元素数据及氨基酸数据分别训练得到相应的溯源模型。

验证评价模块，用于交叉验证基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型的性能；

本发明将个样本数据集随机划分为k个子集(一般是均分)，将一个子集作为验证集，其余的k-1组子集作为训练集；轮流将k个子集作为验证集，交叉重复k次，得到k次的结果，并用k次结果的平均值作为分类器或模型的性能指标。

也就是说，本发明轮流对红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型进行交叉验证，获取各模型的性能。本发明通过溯源的正确率来指示模型的性能。

排序模块，用于对所述基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型按照性能进行排序，根据性能排名赋予鉴别模型相应的权重；

本发明采用正确率来指示模型的性能，因此，正确率越高，鉴别模型的性能越好，溯源精度越高。本发明通过四种鉴别模型对进口牛肉进行综合溯源，因此，为不同的鉴别模型赋予不同的权重，权重代表该鉴别模型对最终溯源结果的影响。权重越高，其对溯源结果影响越大。因此，排名越高的鉴别模型其权重越大，假设基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型的权重分别为ω1、ω2、ω3、ω4，则：

ω1+ω2+ω3+ω4＝1

检测模块，用于分别利用基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型对进口牛肉进行检测，当判断进口牛肉不为相应国家的牛肉时，检测结果为0，否则检测结果为1；

取待测未知过度的牛肉按照上述步骤s1、s2、s3、s4，测定近红外特征光谱数据、稳定同位素质谱数据、矿物元素数据和氨基酸数据，将测得数据分别代入基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型，若预测结果为0，则判断待测进口牛肉不为相应国家的进口牛肉；若预测结果为1，则判断待测样品为相应国家的进口牛肉。

综合溯源模块，用于计算利用基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型得到的检测结果分别与相应的鉴别模型的加权和，将所述加权和与设定阈值进行对比，判断进口牛肉的产地。

本发明通过基于红外光谱数据的鉴别模型、基于同位素质谱数据的鉴别模型、基于矿物元素数据的鉴别模型、基于氨基酸数据的鉴别模型进行混合检测，假设其分别的检测结果为r1、r2、r3、r4，则鉴别模型最终的检测结果为：

r＝r1*ω1+r2*ω2+r3*ω3+r4*ω4

将最终的检测结果与设定的阈值进行比较，假设设置阈值为0.75，当计算的检测结果大于设定阈值时，表面进口牛肉为相应的国家进口，否则不为相应的国家进口。

由此可知，本发明提出进口牛肉产地混合溯源方法及系统，能够针对进口牛肉进行检测，适用范围光，具有极强的应用价值。通过四种溯源方式相结合，克服进口牛肉在运输过程中的影响，提高溯源的精度；同时根据性能对四种溯源方式进行排序，并根据排序结果赋予不同溯源方式相应的的权重，能够根据溯源模型的性能调整其对综合溯源结果的影响，充分发挥各溯源方式的优势，而不是将各溯源方式的简单混合。此外，先对进口牛肉粉碎后再进行干燥，能够提高牛肉的干燥效果，同时采用静置法对进口牛肉样品的进行脱脂，能够快速制得大量的脱脂样品，操作简单，操作效率高。本发明通过多次采集每个进口牛肉样本的光谱数据、同位素数据、矿物元素数据、氨基酸数据，能够避免单次光谱数据采集过程中的误差。本发明对采集的每个进口牛肉样本的光谱数据、同位素数据、矿物元素数据、氨基酸数据，计算相应的方差和标准差，剔除相应的奇异点样本，避免采集环境等对数据的影响。

注意，上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解，本发明不限于这里所述的特定实施例，对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此，虽然通过以上实施例对本发明进行了较为详细的说明，但是本发明不仅仅限于以上实施例，在不脱离本发明构思的情况下，还可以包括更多其他等效实施例，而本发明的范围由所附的权利要求范围决定。