近红外光谱法检测炸油中基值的研究

近红外光谱法检测炸油中基值的研究

一家独特的冷冻食品生产迅速，口感独特，深受人们喜爱。然而,长时间高温油炸会引起煎炸油产生一系列的热氧化反应、聚合反应和水解反应,生成醛酮等羰基化合物。这些产物不仅影响油炸食品的色泽和风味,而且对人体健康产生危害［1］。因此,煎炸油的安全直接影响到煎炸食品的质量与安全。羰基值是煎炸油氧化的重要安全性指标之一［2－3］,传统方法一般采用2,4-二硝基苯肼法(GB/T5009.37—2003)进行羰基值测定,操作复杂、用时较长且消耗大量有机溶剂。近红外光谱(NIR)具有测定快速准确、无需样品预处理和不使用任何化学试剂等优点,已广泛应用于煎炸过程中的油脂氧化检测。国外许多学者利用近红外快速检测油脂煎炸过程中酸值、极性值和过氧化值等［4－9］,而国内研究则主要集中在油炸过程中羰基值的变化规律、影响因素以及与色泽、黏度等其他指标的相关性方面［10－15］,利用近红外光谱快速检测煎炸油羰基值以及利用模型对煎炸过程的羰基值监控均未见报道。本研究以煎炸油为对象,利用近红外光谱建立煎炸油煎炸过程中的羰基值测定模型,对模型进行验证并利用模型对煎炸油煎炸过程中的羰基值进行监控。1材料和方法1.1马铃薯、玉米油、食用调和油花旗起酥油(耐煎炸油)、大豆油(一级,金龙鱼)、菜籽油(三级,金龙鱼)、玉米油(金龙鱼)、食用调和油(金龙鱼):西安益海嘉里油脂工业有限公司;马铃薯:市购;2,4-二硝基苯肼:分析纯,上海研域生物科技有限公司。1.2仪器、分析和光度计MPA-TM傅里叶变换近红外光谱仪(BRUKEROPTICS):德国布鲁克光学仪器公司,检测器为In-GaAs,分析软件为OMNIC7.3和TQAnalyst7.2(美国尼高力公司);UV-1240紫外可见分光光度计:日本岛津公司;WY型油炸锅:煎炸最高温度300℃,广州市唯利安西厨设备制造有限公司。1.3方法1.3.1马铃薯油样品的制备从5种油中任意选取2种油,按质量比1∶1、1∶2、1∶3混合,随机获得样品24个,然后将油样置于105℃烘箱中进行加速氧化,根据氧化时间不同,取原油、氧化程度较轻和氧化程度较高3个阶段的油样;将7kg煎炸油倒入煎炸锅内,加热到油温(180±5)℃,切成片状的马铃薯片以1.5kg/h的速率进行煎炸,每30min取油样1次,每次取100mL油样。每天不间断煎炸6h,共煎炸18h,整个过程不加新油。最终共得108个样品,密封、铝箔纸包裹,置于4℃冰箱中。薯片制作工艺参见文献。1.3.2基值的测量羰基值测定:参照GB/T5009.37—2003。1.3.3温度设定和测量条件用傅里叶近红外光谱仪的液体腔(玻璃管内径为0.5mm)采集油样光谱,波数范围:12000~4000cm－1,光谱分辨率:4cm－1,扫描次数:32次(约30s),液体腔温度设定为60℃。扫描前水浴,使油样温度控制在60℃［6］。尽快完成光谱采集及羰基值测定,以减少环境因素等的影响,提高建模精度。1.3.4最优预处理方法确定光谱预处理采用一阶求导、二阶求导、Norris导数平滑、标准正交变换(SNV)及多元散射校正(MSC)。利用主成分分析确定最佳主成分数,比较不同预处理方法的建模效果确定最优预处理方法,采用PLS法建立模型并进行验证。108个样品被随机分为校正集和验证集,其中校正集样品76个,羰基值范围为0.53~44.17mmol/kg;验证集为32个,羰基值范围为1.98~43.51mmol/kg。其中模型评价指标:相关系数(R)、校正标准偏差(RMSEC)及预测标准偏差(RMSEP)。1.3.5精密度与准确度本试验采用t检验来判断近红外光谱法替代2,4-二硝基苯肼法是否可行,采用平均相对误差及相对标准偏差(RSD)分别对羰基值的验证集样品进行准确度和精密度检验。相关指标的计算公式如下:式中:n为样品测量次数;xi为样品测定值;x珋为n次测定结果的平均值。1.3.6近红外模型法使用耐煎炸油进行煎炸,煎炸过程同文献。每1h采集1个样品,取100ml油样,冷却至60℃,直接利用近红外模型进行羰基值测定,同时利用传统方法测定羰基值,对两者结果进行相关性分析。每天连续煎炸6h,当羰基值超过25mmol/kg时,停止煎炸,整个过程不加新油。2结果与分析2.1化学测定值的分析建立校正模型前,收集有代表性的样品尤为重要［17］。样品包含按不同比例混合的不同氧化程度的混合油及不同煎炸阶段的煎炸油,食用油羰基值的化学测定值分析结果见表1。由表1可以看出,校正集样品的羰基值范围为0.53~44.17mmol/kg,验证集样品羰基值的范围为1.98~43.51mmol/kg,《食用油煎炸过程中的卫生标准》(GB7102.1—2003)规定食用油羰基值≤25mmol/kg,所采集的样品覆盖了符合标准的油和不符合标准的油,样品具有较强的代表性。2.2c—光谱分析由于原始光谱数据受随机噪音、样品不均等影响,使得光谱中包含干扰信息,因此必须对光谱进行预处理来提取有效光谱信息［18］,如求导、微分、平滑、归一化等。食用油原光谱及预处理光谱见图1。由图1a可得,在8280、7074、5801~5795、5678、4663cm－1处有较强的吸收峰。羰基化合物的吸收谱带很宽,基频、一级倍频、二级倍频和三级倍频的吸收谱带分别为:1950~1550cm－1、3824~2818cm－1、5680~4043cm－1、7550~5166cm－1［19］,由此可得7074、5801~5795、5678、4663cm－1处的吸收可能与羰基的三级倍频、二级倍频有关,其中在4695cm－1附近的吸收峰与C=O伸缩振动的合频吸收有关［20］。取代基的同位素效应、质量变化、取代基团、共振等及它们之间的相互作用都会影响羰基化合物的吸收谱带位置［19］。8280cm－1处的波峰可能与C—H伸缩振动的二级倍频吸收有关［21］。图1b~图1d分别是原始光谱经过一阶求导、二阶求导及二阶求导+SF(9,6)处理后的图谱,其中SF(9,6)表示Savitzky-Golay多项式平滑处理,移动窗口宽度为9,多项式次数为6。不同的预处理方法对光谱信息的处理结果不同。光谱的低阶微分处理对噪声影响敏感性较低,微分处理(图1b和图1c)可去除部分线性或接近线性的背景、噪声光谱对目标光谱的影响［22］,强化谱带特征,提高分辨率和灵敏度;原始光谱经二阶导数和Savitzky-Golay平滑处理(图1d)后信噪比提高,可能是由于二阶求导与平滑相结合不仅消除了背景噪音的影响,而且降低了由于微分处理引起的光谱噪音增大的影响,同时,合适的平滑点数进一步提高了信噪比［18］,最终达到改善光谱性状的目的［23］。2.3求导条件下的校正集模型建立主成分分析的目的是将数据降维,以除去化学信息中相互重叠部分。通过交互检验,得出不同主成分数下的预测残差平方和(PRESS)。最小的PRESS值对应的因子数作为模型的主成分数。在9739~6274cm－1范围内,采用SNV和二阶求导条件下建立校正集模型,模型主成分分析结果见图2。由图2可知,主成分数为1时,全光谱即12000~4000cm－1贡献率为53.003%,选择范围即9739~6274cm－1的贡献率为68.958%,校正模型的RMSEC值是8.26mmol/kg;主成分数为9时,全光谱贡献率为97.104%,选择范围光谱贡献率为99.404%,校正模型的RMSEC值是4.08mmol/kg;主成分为10时,校正模型的RMSEC值最小,为3.99mmol/kg。校正模型采用10个主成分数建模能代替更多的光谱信息,校正模型中采用10个主成分。2.4阶导数和nf+二阶求导的光谱预处理效果在确定主成分的基础上,对光谱进行一阶求导、二阶求导、MSC、Norris平滑、Savitzky-Golay预处理,以消除样本不均匀、光散射、基线漂移等因素的影响。表2为食用油羰基值光谱不同预处理方法建模效果。由表2可知,在9739~6274cm－1范围内,对原始光谱分别进行单一的一阶求导或二阶求导,均可提高光谱特性,表明导数处理可以去除背景噪音,强化谱带特性。但当同时采用二阶导数与NF平滑处理,R降低,可能因为平滑过度,部分有用的信息被过滤［24］。采用NF+二阶求导处理光谱,效果较差,而在同样波段下,采用SNV(SF或MSC)+二阶求导,或者单一采用二阶求导,效果均较好,说明在此波数范围内,采用NF+二阶求导的光谱预处理效果低于其它四种处理方法。在9739~6274cm－1范围内,采用SNV和二阶求导进行光谱预处理,模型的相关系数最大(R=0.9844),模型效果最好,故将其确定为最优建模条件。2.5化学测定和近红外模型预测基值分布最优建模条件下建模,即在9739~6274cm－1波数范围内采用SNV和二阶求导,得到校正集模型并对模型进行验证。校正集和验证集的化学测定羰基值和近红外模型预测羰基值分布见图3和图4。由图3、图4可以看出,校正集中的化学测定羰基值和近红外模型预测羰基值相关系数R为0.9844,校正误差均方根RMSEC为1.47mmol/kg;验证集中的化学测定羰基值和近红外模型预测羰基值相关系数R为0.9790,验证误差均方根RMSEC为1.73mmol/kg。模型相关性高,表明羰基值模型可用。2.6近红外光谱法测定基值的准确度本试验对预测集值和实测值进行配对t检验,其t值为1.93,小于临界值t31，0．05=2.04,表明可以利用近红外光谱法快速测定羰基值。随机取5个样品测定其平均相对误差,相对误差范围为1.2%~4.3%,均小于5%,表明模型具有较高的准确度;采用模型法平行测定样品羰基值5次,其RSD为1.53%,表明模型具有较高的精密度,即稳定性较好。2.7实测值与预测值对比使用耐煎炸油进行煎炸,每小时采集1个样品,待样品冷却至60℃,直接利用近红外模型测定羰基值,对煎炸油煎炸过程进行监控。当羰基值超过国家标准(25mmol/kg)时停止煎炸。同时,对各阶段的油样采用化学法进行测定。2种方法所测羰基值相关性见图5。由图5可知,实测值与预测值决定系数R2高达0.9952,线性相关性良好,模型可用于煎炸油煎炸过程监控。在实际应用中,可与油脂企业的侧线型自动采样分析系统相结合,获得实时油样,通过近红外获得光谱,利用模型对油脂的羰基值进行实时监控［19］。同时,煎炸13h后,煎炸油的羰基值为29.53mmol/kg,明显超出国家标准(25mmol/kg),停止煎炸。2种方法测得的羰基值均随煎炸时间增长而线性增加,这与陈锋亮等［13］采用传统方法测得的煎炸过程中的羰基值变化规律相一致,进一步表明近红外光谱技术可用于