一种基于近红外光谱分析技术的植物蛋白饮料中主要成分含量快速检测的方法与流程

本发明属于食品质量与安全快速检测技术领域，具体涉及一种利用近红外光谱分析技术对植物蛋白饮料中蛋白质、脂肪及可溶性固形物含量的快速检测方法。

背景技术：

植物蛋白饮料，在当前市场中，由于其营养、保健、天然、口感优良等优点而备受消费者喜爱。蛋白质、脂肪及可溶性固形物作为植物蛋白饮料中的三大主要成分，直接影响着饮料的风味、营养及稳定性。因此，对植物蛋白饮料中的蛋白质、脂肪及可溶性固形物含量进行及时控制十分必要，对保证饮料品质起到重要的作用。

目前植物蛋白饮料中蛋白质、脂肪及可溶性固形物的测定分别为凯氏定氮法、酸水解法及折光计法，但是，这些方法普遍存在检测过程耗时、操作较为复杂等缺陷，且对检测人员的技术水平要求较高，无法实现植物蛋白饮料中主要成分含量的大批量快速定量检测，无法及时对产品质量及稳定性进行控制，不适合企业用于样品的快速检测。

近红外光谱（nearinfraredspectroscopy，nir）技术是近年来迅速发展起来的快速无损检测技术，几乎可用于所有与含氢基团有关的样品物化性质分析，具有操作简单、样品无需前处理、易于实现快速分析的优势，已经在食品、药品、石油化工、农牧等领域广泛应用，市场前景好。但在国内至今为止尚未见有应用近红外光谱分析技术对植物蛋白饮料中成分进行快速检测的相关报道。

技术实现要素：

本发明目的在于提供一种基于近红外光谱分析技术的植物蛋白饮料中蛋白质、脂肪及可溶性固形物含量的快速检测方法，该方法相比现有检测方法所存在的操作复杂、费力费时等缺陷，具有快速、方便等优点。

本发明所提供的植物蛋白饮料中主要成分含量快速检测方法，包括如下步骤：

1）收集不同种类、不同品牌的植物蛋白饮料样品，随机选取若干作为样品集，利用近红外光谱仪采集样品集的近红外光谱图；

2）检测得到样品集植物蛋白饮料样品中蛋白质、脂肪及可溶性固形物含量的实测值；

3）选取建模样品，按照2:1的比例，采取kennard-stone（k-s）法进行样本集划分，分为校正集和验证集，并将步骤1）中得到的校正集中植物蛋白饮料样品的近红外光谱数据分别与步骤2）中得到的蛋白质、脂肪及可溶性固形物含量的实测值相对应，通过化学计量学软件建立酿酒葡萄样品中酒石酸及苹果酸含量的定量校正模型，并用验证集样品对所建定量校正模型进行外部验证，最终得到植物蛋白饮料中主要成分含量近红外校正预测模型；

4）取待测样品，按步骤1）中的光谱测量条件采集待测样品的近红外光谱数据，导入校正预测模型中，经模型运算，即可得到未知植物蛋白饮料样品中的蛋白质、脂肪及可溶性固形物含量。

上述检测方法中，步骤1）中，所述收集的植物蛋白饮料样品数至少为800，优选为800~1200。所述采集样品集中植物蛋白饮料样品的近红外光谱图的方法如下：测量前将样品放置于室温下统一静置30min，采用透反射方式进行样品的光谱采集，每一份样品体积为25ml，样品测量温度为40摄氏度，样品倒入样品杯，放入透反射盖，并用透反射盖赶走气泡，近红外光透过样品后在透反射盖上发生漫反射，漫反射光又透过样品进入检测器中，扫描次数为32次，波长范围为4000-10000cm^-1，透反射盖下样品厚度固定为0.3mm，每个样品重复装样采集2次光谱，光谱以吸光度log(1/r)形式存储。

所述植物蛋白饮料样品涵盖不同产地、不同生产日期、不同型号规格的样品。

所述近红外光谱仪为nirmasterm54p傅里叶变换近红外光谱仪（瑞士步琪有限公司）。

所述化学计量学软件是由瑞士步琪有限公司提供。

上述检测方法中，步骤2）中，所述检测的方法分别为gb/t5009.5-2010《食品中蛋白质的测定》的凯氏定氮法、gb/t5009.6-2003《食品中脂肪的测定》的酸水解法以及gb/t12143-2008《饮料通用分析方法》的折光计法。

上述检测方法中，步骤3）中，所述对所建立的近红外定量校正模型进行外部验证的步骤：将步骤1）得到的验证集植物蛋白饮料样品的近红外光谱导入所建立的定量校正模型中，得到验证集中样品蛋白质、脂肪及可溶性固形物含量的预测值，将其与步骤2）中所述验证集中植物蛋白饮料样品的蛋白质、脂肪及可溶性固形物含量的实测值相比，进行模型准确度的检验，若预测值（%）与实测值（%）差值的绝对值与实测值（%）之比在设定的范围内，则所述定量校正模型可用；反之，则需重复步骤3），优化建模条件直至所述定量校正模型可用。

所述定量校正模型依次采用如下步骤而建立得到：

a、光谱预处理方法：所述光谱预处理方法选自如下至少一种：多元散射校正（msc）、savitzky-golay卷积平滑、savitzky-golay一阶导数、小波变换（wt）和ncl；

b、样品集划分方法：所述样品集划分方法选自如下任意一种：随机法、kennard-stone（k-s）法、spxy法；

c、变量压缩方法：所述变量压缩方法选自如下任意一种：cars法、无信息变量消除法、组合区间偏最小二乘法、遗传算法；

d、化学计量学建模方法：所述化学计量学建模方法选自如下任意一种：偏最小二乘法（pls）、主成份回归（pcr）、最小二乘支持向量机（ls-svm）或人工神经网络法（ann）。

所述定量校正模型具体依次采用如下步骤而建立得到：

a、光谱预处理方法：ncl；

b、样品集划分方法：kennard-stone（k-s）法；

c、变量压缩方法：组合区间偏最小二乘法、遗传算法；

d、化学计量学建模方法：偏最小二乘法（pls）。

上述检测方法中，步骤3）中,所述变量压缩方法的实现步骤：将植物蛋白饮料4000-10000cm-1全光谱分为k个子区间（k=10~40，间隔5），在不同子区间数下，分别就不同组合数（1~4）进行组合区间偏最小二乘（sipls）计算，再以交互验证标准差rmsecv最小为标准将sipls筛选出的建模波段使用遗传算法进行计算，分别筛选出植物蛋白饮料的蛋白质、脂肪及可溶性固形物的最佳建模变量。

上述检测方法中，步骤3）中,所述变量压缩方法中遗传算法的运行参数设置为：初始种群大小80，变异概率pm=0.01，交叉概率pc=0.5，最大因子数10，遗传迭代次数120次，以rmsecv值确定出最佳的建模变量。

上述检测方法中，步骤3）中,所述最佳建模变量为：

1）经sipls筛选出的植物蛋白饮料各指标建模波段分别为：蛋白质（4404~4800、5604~6000、6804~7200、9204~9600cm^-1）；脂肪（4244~4480、4724~4960、5684~5920、8804~9040cm^-1）；可溶性固形物（4000~4480、4964~5200、7604~7840cm^-1）；

2）经遗传算法筛选出的植物蛋白饮料各指标最佳建模变量分别为：蛋白质（4404、4460、4500、4528、4552、4580、4596、4612、4624、4628、4672、4680、4700、4716、5608、5616、5620、5624、5640、5644、5676、5680、5688、5700、5708、5720、5732、5752、5768、5772、5780、5788、5792、5796、5800、5804、5808、5816、5820、5824、5828、5840、5844、5852、5864、5936、5964、6816、6964、7112、7148cm^-1）；脂肪（4272、4292、4312、4316、4320、4324、4328、4332、4344、4348、4352、4356、4380、4384、4404、4408、4432、4436、4464、4800、4832、4928、4936、5692、5716、5724、5728、5732、5736、5740、5744、5748、5756、5788、5800、5804、5808、5812、5816、5828、5832、5836、5844、5852、5856、5864、5868、5872、5876、5892、5920、8868、8936cm^-1）；可溶性固形物（4096、4228、4232、4244、4248、4288、4292、4300、4340、4368、4372、4376、4380、4384、4388、4392、4396、4400、4404、4408、4412、4416、4420、4428、4432、4436、4440、4444、4452、4456、4460、4464、4468、4472、4476、4480、4980、4988、5036、5080、7620、7672cm^-1）。

本发明为植物蛋白饮料蛋白质、脂肪及可溶性固形物含量的测定提供一种全新的快速、准确、简单的检测方法，对植物蛋白饮料的品质及稳定性控制以及整个行业意义重大。

与现有技术相比，本发明的技术优势体现在如下方面：（1）本发明操作步骤简单，操作过程简洁，为企业节省了传统检测所需的人力物力成本；（2）适用于植物蛋白饮料主要成分含量的快速高效分析，具有检测速度快、分析效率高、稳定性及重现性好的优势，采用本发明可于40秒内完成植物蛋白饮料中蛋白质、脂肪及可溶性固形物含量的测定，降低了检测成本同时大大提高检测效率；（3）在检测过程中不使用任何化学试剂，也不对环境产生任何污染，提供了一种绿色可靠的植物蛋白饮料中蛋白质、脂肪及可溶性固形物含量的分析方法。

附图说明

图1为建模用的植物蛋白饮料蛋白质、脂肪及可溶性固形物含量的光谱图。

图2、图3、图4为植物蛋白饮料蛋白质、脂肪及可溶性固形物含量近红外预测模型效果图。

具体实施方式

下面通过具体实施例对本发明的方法进行说明，但本发明并不局限于此，凡在本发明的精神和原则之内所做的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

下述实施例中所述实验方法，如无特殊说明，均为常规方法；所述试剂和材料，如无特殊说明，均可从商业途径获得。

实施例1、基于近红外光谱分析技术的植物蛋白饮料主要成分含量快速检测方法

具体操作步骤如下：

（1）收集样品和采集样品近红外光谱：收集不同种类、不同品牌的植物蛋白饮料样品，按照2:1的比例，采取kennard-stone（k-s）法进行样本集划分，分为校正样品集和验证样品集，最终选择786份样品组成校正样品集，388份样品为验证样品集。样品集涵盖不同产地、不同生产日期、不同型号规格的样品；

采用透反射方式进行样品的光谱采集，样品倒入样品杯，放入透反射盖，并用透反射盖赶走气泡，近红外光透过样品后在透反射盖上发生漫反射，漫反射光又透过样品进入检测器中，扫描次数为32次，波长范围为4000-10000cm^-1。透反射盖下样品厚度固定为0.3mm，每个样品重复装样采集2次光谱。收集所得2条光谱以吸光度log(1/r)形式均加入植物蛋白饮料样品的近红外光谱图库，最终所获校正样品集的光谱如图1所示。

（2）化学法测定植物蛋白饮料样品中的蛋白质、脂肪及可溶性固形物含量：分别依据gb/t5009.5-2010《食品中蛋白质的测定》的凯氏定氮法、gb/t5009.6-2003《食品中脂肪的测定》的酸水解法以及gb/t12143-2008《饮料通用分析方法》的折光计法测得样品中蛋白质、脂肪及可溶性固形物含量的实测值，构建校正样品集与验证样品集的近红外模型的参比值数据集。

（3）近红外定量校正模型的建立和验证样品集对近红外定量校正模型进行检验评价：将步骤（1）所得校正样品集的近红外光谱图分别与步骤（2）中所得的校正样品集中样品的蛋白质、脂肪及可溶性固形物含量的实测值一一对应，通过化学计量学软件分别建立植物蛋白饮料中蛋白质、脂肪及可溶性固形物含量的定量校正模型，所述定量校正模型依次采用如下步骤而建立得到：

a、光谱预处理方法：ncl；

b、样品集划分方法：kennard-stone（k-s）法；

c、变量压缩方法：组合区间偏最小二乘法、遗传算法；

d、化学计量学建模方法：偏最小二乘法。

其中变量压缩方法通过如下步骤实现：将植物蛋白饮料4000-10000cm-1全光谱分为k个子区间（k=10~40，间隔5），在不同子区间数下，分别就不同组合数（1~4）进行组合区间偏最小二乘（sipls）计算，再以交互验证标准差rmsecv最小为标准将sipls筛选出的建模波段使用遗传算法进行计算，分别筛选出植物蛋白饮料的蛋白质、脂肪及可溶性固形物的最佳建模变量。

遗传算法运行参数设置为：初始种群大小50，变异概率pm=0.01，交叉概率pc=0.5，最大因子数10，遗传迭代次数120次，以rmsecv值确定出最佳的建模变量。

经sipls初步筛选出的植物蛋白饮料各指标建模波段分别为：蛋白质（4404~4800、5604~6000、6804~7200、9204~9600cm^-1）；脂肪（4244~4480、4724~4960、5684~5920、8804~9040cm^-1）；可溶性固形物（4000~4480、4964~5200、7604~7840cm^-1）。

经遗传算法筛选出的植物蛋白饮料各指标最佳建模变量分别为：蛋白质（4404、4460、4500、4528、4552、4580、4596、4612、4624、4628、4672、4680、4700、4716、5608、5616、5620、5624、5640、5644、5676、5680、5688、5700、5708、5720、5732、5752、5768、5772、5780、5788、5792、5796、5800、5804、5808、5816、5820、5824、5828、5840、5844、5852、5864、5936、5964、6816、6964、7112、7148cm^-1）；脂肪（4272、4292、4312、4316、4320、4324、4328、4332、4344、4348、4352、4356、4380、4384、4404、4408、4432、4436、4464、4800、4832、4928、4936、5692、5716、5724、5728、5732、5736、5740、5744、5748、5756、5788、5800、5804、5808、5812、5816、5828、5832、5836、5844、5852、5856、5864、5868、5872、5876、5892、5920、8868、8936cm^-1）；可溶性固形物（4096、4228、4232、4244、4248、4288、4292、4300、4340、4368、4372、4376、4380、4384、4388、4392、4396、4400、4404、4408、4412、4416、4420、4428、4432、4436、4440、4444、4452、4456、4460、4464、4468、4472、4476、4480、4980、4988、5036、5080、7620、7672cm^-1）。

相应的植物蛋白饮料蛋白质、脂肪及可溶性固形物含量近红外预测模型效果图分别如图2、图3及图4所示，从图可知植物蛋白饮料蛋白质、脂肪及可溶性固形物的近红外预测值与实测值均匀分布在对角线两侧，说明模型能够对样品蛋白质、脂肪及可溶性固形物含量进行预测。将验证样品集中样品按照与步骤（1）中建模样品集中样品相同的处理方式进行处理，将步骤（1）得到的验证集样品的近红外光谱导入所建立的定量校正模型中，得到验证样品集中蛋白质、脂肪及可溶性固形物含量的预测值，将其与步骤（2）中所述验证集植物蛋白饮料中蛋白质、脂肪及可溶性固形物含量的实测值相比，进行模型准确度的检验，若预测值（%）与实测值（%）差值的绝对值与实测值（%）之比在国标要求范围之内，则所述定量校正模型可用；反之，则需重复步骤（3），优化建模条件直至所述定量校正模型可用。定量校正模型验证结果如表1所示。可以看到，模型验证结果良好，准确度满足需求。

表1验证集近红外蛋白质、脂肪及可溶性固形物结果分析

（4）待测样品测试：将50个未参与建模的植物蛋白饮料样品组成预测样品集，分别依据gb/t5009.5-2010《食品中蛋白质的测定》的凯氏定氮法、gb/t5009.6-2003《食品中脂肪的测定》的酸水解法以及gb/t12143-2008《饮料通用分析方法》的折光计法测得样品中蛋白质、脂肪及可溶性固形物含量，作为实测值。同步骤（1）相同的近红外光谱扫描条件下进行近红外光谱扫描，将所得光谱数据输入步骤（3）所得植物蛋白饮料中蛋白质、脂肪及可溶性固形物含量的定量校正模型中，得到待测植物蛋白饮料中蛋白质、脂肪及可溶性固形物含量的预测值，根据国家标准，植物蛋白饮料中蛋白质、脂肪及可溶性固形物模型检测结果如下表2所示，从表2可知：蛋白质、脂肪及可溶性固形物指标分别有98%、100%、100%的测试结果满足国家标准，说明本发明在植物蛋白饮料中蛋白质、脂肪及可溶性固形物的含量测定中具有良好的应用效果。

表2植物蛋白饮料蛋白质、脂肪及可溶性固形物模型检测结果