滩羊肉冷藏期间TBARS含量的高光谱检测

2022-03-11樊奈昀刘贵珊张晶晶袁瑞瑞孙有瑞

光谱学与光谱分析 2022年3期

樊奈昀，刘贵珊，张晶晶，袁瑞瑞，孙有瑞，李月

宁夏大学食品与葡萄酒学院，宁夏银川 750021

引言

脂肪氧化是羊肉品质劣化和新鲜度下降的重要影响机制，仅次于微生物引起的变化[1]。硫代巴比妥酸反应物(thiobarbituric acid reactive substances，TBARS)被广泛用于评估肉品中脂肪氧化程度，是反映脂肪氧化程度的重要氧化参数[2]。因此，实时检测滩羊肉中TBARS含量对保障肉品安全可食具有重要意义。分光光度法是测定TBARS含量的常规方法，具有破坏性、操作繁琐、一定的危害性和污染性等缺点。

高光谱成像是一种先进的光学检测技术，能同时捕获待测样本内部成分的光谱数据和外部属性的图像信息，通过分析频谱波形和强度信息可实现肉品品质的定性或定量检测[3-4]。然而，肉的近红外光谱通常复杂且重叠，不仅包含有用的光学、物理、化学信息，还包含被测样品内部因素产生的其他光谱响应[5]。因此，更好地了解光谱信息和更准确的波段分配有助于开发稳健的校正模型[6]。二维相关光谱技术(two-dimensional correlation spectra，2DCOS)将光谱信号扩展到第二维上，在研究体系中添加一定的外扰诱发光谱信号发生动态变化，通过解析二维相关光谱谱图特征寻找与微扰相关的敏感变量[7]。目前，高光谱技术结合2DCOS分析无损检测羊肉中TBARS含量的研究鲜有报道。

以滩羊肉中TBARS含量为外界扰动，分析随外部扰动变化的光谱特征信息，确定与TBARS含量相关的特征变量，探究高光谱技术结合2DCOS无损检测滩羊肉中TBARS含量的可行性并建立分析模型，为实时快速检测TBARS含量提供参考。

1 实验部分

1.1 材料与试剂

新鲜羊肉背最长肌购买于宁夏盐池县鑫海食品有限公司，将样本存储于0～4 ℃的便携式冷藏培养箱中，2 h内运输到实验室。使用无菌刀具除去肉样表面的脂肪和结缔组织，整形切块，尺寸约为40 mm×40 mm×10 mm(长×宽×厚)，共获得羊肉样本180个。所有样本均单独真空包装，随机分为10组(n=18)，4 ℃贮藏，分别在贮藏期为1，3，5，7，9，11，14，16，18和20 d时进行测试。

1.2 仪器

可见/近红外高光谱成像系统(美国Headwall Photonics公司)和UV1902PC型紫外-可见分光光度计(上海棱光技术有限公司)。推扫式高光谱成像系统主要由6部分组成：光谱范围为400～1 000 nm的线扫描成像光谱仪，卤素灯光源，G4-232增强型EMCCD相机，VT-80精密电控位移平台，计算机和数据采集软件。光谱成像仪共有125个波段，光谱分辨率为2.8 nm，狭缝宽度为25 μm，光源为卤素灯双线光源。

1.3 TBARS含量测定[2]

准确称取5 g切碎肉样，将25 mL 20%的三氯乙酸溶液和20 mL预冷去离子水添加到肉样中均匀混合，使用均质器均质，室温下静置20 min。提取液在8 000 g条件下离心10 min后过滤，然后使用去离子水将上清液稀释至50 mL。最后，取5 mL稀释后的上清液与0.02 mol·L-1的2-硫代巴比妥酸水溶液(5 mL)混合，沸水浴反应30 min。将反应溶液通过流动水冷却至室温后，以去离子水为空白，在532 nm下读取反应液的吸光度值。按下式计算得到的TBARS值：

TBARS=A532×7.8

其中，A532是反应液在532 nm处的吸光度值。

2 结果与讨论

2.1 样本集划分

常规化学分析方法测定的TBARS含量统计结果如表1所示，TBARS值具有较大的变化范围(0.16～1.47 mg·kg-1)，有利于开发相对稳健的校正模型。采用Kennard-Stone(KS)算法按3∶1的比例选取135个样本作为校正集建立模型，其余45个样本为预测集用于校正模型性能的验证。KS算法划分样本集时首先选择具有最长欧氏距离的两个矢量进入训练库，在接下来的迭代过程中，选择拥有最大最小距离的待选样本进入训练库，以此类推得到所需的样本数量[8]。

表1 滩羊肉TBARS含量统计Table 1 Statistics of measured TBARS contents in Tan mutton

2.2 预处理方法优选

表2 不同预处理方法下TBARS含量的PLSR模型Table 2 PLSR modeling results of TBARS content with different pretreatment methods

图1 TBARS含量实测值与预测值的线性拟合图(a):Raw-PLSR;(b):SG-PLSR;(c):De-trending-PLSR;(d):SG+De-trending-PLSRFig.1 Linear fitting diagrams of measured and predicted values of TBARS contents(a):Raw-PLSR;(b):SG-PLSR;(c):De-trending-PLSR;(d):SG+De-trending-PLSR

2.3 2DCOS分析

化合物的含氢基团分子键(O—H，C—H和N—H)在受到光照射时，会造成分子键振动能的改变。通过解析波形和频谱强度信息可以实现滩羊肉中TBARS含量的定量检测。图2为以TBARS含量为微扰量，由滩羊肉的可见近红外光谱进行2DCOS分析产生的二维相关光谱及其切割光谱。如图2(b)所示；在同步二维相关光谱的对角线位置出现了4个主要的自相关峰，分别为579，699，756和867 nm。这些自相关峰的出现表明随着TBARS含量的增大，自相关峰所在波段的光谱吸收强度很容易发生变化，该处的光谱信号对外扰较敏感，是与TBARS值相关的敏感变量。与自相关峰对应的主要交叉峰位于同步二维相关光谱的非对角线位置[图2(a)]，表2总结了各交叉峰的正负性。正峰表明v1和v2处的信号强度沿相同方向变化，即一起增加或减少，负峰表明v1和v2处的信号强度沿相反方向变化。结合2DCOS分析，选择579～867 nm光谱范围内的特征变量用作为滩羊肉中TBARS含量的检测。

图2 二维相关光谱及其切割光谱(a)：同步谱图；(b)：切割光谱Fig.2 Two-dimensional correlation spectra and its slice spectra(a):Synchronous contour map;(b):Slice spectra

表2 二维相关分析产生的交叉峰符号Table 2 Signs of cross-peaks generated from two-dimensional correlation analysis

2.4 特征波长优选

采用变量组合集群分析(variable combination population analysis，VCPA)[11]、竞争性自适应加权算法(competitive adaptive reweighted sampling，CARS)[12]和连续投影算法(successive projections algorithm，SPA)[13]从全光谱数据(FS)和2DCOS分析中提取具有最多相关信息的特征波长，建立简化模型并进行比较。由图3(a)可知，当采样次数为36时，RMSECV下降到0.136 5，VCPA从FS中选择了7个特征波长，占到总波长的5.6%。CARS算法根据达尔文进化论中“适者生存”的原则淘汰变量，在重复循环采样后提取代表性波长。使用CARS算法从FS中选择了16个特征变量，占到总波长的12.8%。SPA通过RMSECV值评估波长的候选子集，然后使用消除程序除去不相关的波长。当RMSECV最小值为0.141时，SPA从FS中选择了20个特征变量，占总波长数的16%。采用VCPA，CARS 和SPA算法对2DCOS分析所确定的研究区域进行二次筛选，分别优选出了8，24和14个特征波长，具体特征波长提取结果见表3。

图3 变量选择算法提取特征波长(a):VCPA算法;(b):SPA算法Fig .3 Characteristic wavelength extracted by variable selection algorithm(a):VCPA algorithm;(b):SPA algorithm

表3 特征波长提取结果Table 3 The results of extracting characteristic wavelengths

2.5 特征波长预测性能分析

表4 不同波长选取方法所建立的PLSR模型效果对比Table 4 Comparison of PLSR models based on different wavelength extraction methods

2.6 定量分析模型的开发

基于2DCOS+CARS方法选取的关键性波长建立了光谱数据与滩羊肉中TBARS参考值之间的定量分析模型，模型函数为：Y(TBARS)=-0.15+2.99λ588-7.01λ593+7.45λ598-6.14λ603+7.06λ612-8.25λ622+2.64λ631-4.18λ636+13.91λ646-11.3λ655+12.64λ675-8.51λ684-7.81λ689+1.08λ703-2.54λ713+5.47λ727+6.62λ742+5.69λ751+2.48λ775-1.93λ780-6.95λ790+7.09λ799-3.56λ809+1.82λ819。λ表示特征波长所对应的光谱反射信息，本模型可为TBARS含量实时在线检测研究提供理论参考。