馬劭瑾 李永玉 彭彥昆 閆 帥 劉亞超

(中國農(nóng)業(yè)大學(xué)工學(xué)院， 北京 100083)

0 引言

葵花籽油含豐富的亞油酸、油酸、不飽和脂肪酸等營養(yǎng)物質(zhì)，在世界食用油中產(chǎn)量位居第二[1]。近些年來，我國葵花籽油的消費(fèi)量也不斷增加。但是葵花籽油中的高不飽和脂肪酸含量也導(dǎo)致其更易被氧化，在售賣或家庭使用過程中都很容易發(fā)生氧化變質(zhì)[2]。變質(zhì)的葵花籽油會產(chǎn)生多種自由基、氫過氧化物和有毒聚合物，容易引發(fā)多種疾病。在植物油食品安全國家標(biāo)準(zhǔn)中明確規(guī)定，食用植物油過氧化值不能超過0.25 g/(100 g)，在各地食品安全抽檢信息中植物油過氧化值超標(biāo)等頻見報道。

目前食用油過氧化值的傳統(tǒng)檢測方法主要為化學(xué)滴定法，具有操作復(fù)雜、耗時較長、需采用有毒化學(xué)試劑等缺點(diǎn)，而且不能滿足快速實(shí)時的需求。為克服傳統(tǒng)化學(xué)手段在食用油產(chǎn)銷鏈中品質(zhì)監(jiān)測的局限性，諸多學(xué)者進(jìn)行了食用油氧化指標(biāo)快速檢測方法的研究[3-8]。其中，熒光光譜以其靈敏度高、能直接反映食用油中與氧化機(jī)制直接相關(guān)的如酚類物質(zhì)和維生素E等熒光物質(zhì)的特性，在食用油氧化程度評估方面得到迅速發(fā)展。諸多研究者均基于三維熒光光譜[9-11]，或者從三維熒光光譜中提取合適的波長間隔，探究同步熒光光譜與氧化食用油的相關(guān)關(guān)系[12-20]。三維熒光光譜能夠同時反映熒光強(qiáng)度隨激發(fā)、發(fā)射波長的變化，但也因此具有較慢的掃描速度，且獲取三維熒光光譜的熒光分光光度計體積較大、價格昂貴，并不具備食用油品質(zhì)安全監(jiān)測所需的高效性與便捷性。激光誘導(dǎo)熒光(Laser induced fluorescence, LIF)光譜則可以彌補(bǔ)上述熒光檢測方法的不足，可以實(shí)現(xiàn)食用油過氧化值的現(xiàn)場實(shí)時監(jiān)測，然而基于激光誘導(dǎo)熒光對食用油過氧化值進(jìn)行定量檢測的相關(guān)研究尚未見報道。

本文針對食用油易被氧化產(chǎn)生多種自由基、氫過氧化物和有毒聚合物等特性，根據(jù)食用油在產(chǎn)銷鏈中品質(zhì)監(jiān)測的實(shí)用需求，基于激光誘導(dǎo)熒光光譜搭建食用油快速檢測系統(tǒng)，編寫相關(guān)采集控制軟件，并以葵花籽油為研究對象，確定熒光響應(yīng)時間等光譜采集參數(shù)，基于自行搭建系統(tǒng)采集不同過氧化值葵花籽油的熒光光譜，探討不同預(yù)處理對熒光光譜的校正效果，建立葵花籽油過氧化值實(shí)時定量預(yù)測模型，最終實(shí)現(xiàn)對葵花籽油過氧化值的無損、快速檢測。

1 實(shí)驗材料與方法

1.1 實(shí)驗材料

葵花籽油樣品購于物美超市，分為生產(chǎn)日期2019年11月4日和2020年6月10日兩批次。為制備不同氧化程度葵花籽油樣品，將購買的葵花籽油分別放置于60℃恒溫箱和25℃的室溫環(huán)境中保存。60℃恒溫箱中，加熱氧化樣品每隔6 h取樣一次，室溫氧化處理樣品每隔5 d取樣一次，共制備90個不同氧化程度葵花籽油樣品。

1.2 過氧化值(PV)的標(biāo)準(zhǔn)理化值測定方法

參照GB 5009.227—2016中的滴定法，每個葵花籽油樣品平行測定3次取平均值作為過氧化值的標(biāo)準(zhǔn)理化值。

2 葵花籽油熒光光譜檢測系統(tǒng)搭建

以食用油品質(zhì)安全現(xiàn)場實(shí)時監(jiān)測技術(shù)裝置研發(fā)為目的，自行搭建了基于激光誘導(dǎo)熒光光譜的葵花籽油檢測系統(tǒng)，主要包括光源部分、樣品室、光譜采集部分、供電部分以及控制部分，如圖1所示。光源部分包括紫外LED光源、濾光片與準(zhǔn)直透鏡。LED光源功率為3 W，中心波長為(365±10)nm，LED激發(fā)光源發(fā)出的紫外光經(jīng)過窄帶濾光片和準(zhǔn)直透鏡平行照射到樣品室內(nèi)的樣品。樣品室在相鄰兩側(cè)進(jìn)行開孔設(shè)計，用于進(jìn)行直角熒光檢測，并且可以保證在放入比色皿后，激發(fā)光入射和發(fā)射光出射的位置均位于比色皿中部。樣品室配有遮光蓋，可以確保在比色皿放入后進(jìn)行光譜采集時不受外界雜散光的影響。光源系統(tǒng)與樣品室采用一體化設(shè)計，從而避免比色皿放入或取出時檢測系統(tǒng)發(fā)生晃動。光譜采集部分主要包括光纖與波長范圍為200～1 100 nm的光譜儀(USB2000+型)，光譜分辨率為1 nm。激發(fā)樣品產(chǎn)生的熒光信號，經(jīng)垂直于入射光方向安裝的光纖探頭傳入光譜儀中，再由光譜儀通過USB接口將數(shù)據(jù)傳輸?shù)接嬎銠C(jī)。供電部分主要是由小型恒流LED電源驅(qū)動器與電源組成?？刂撇糠种糜跇悠肥蚁路?，分別設(shè)計光源和光譜儀基座以確保光源部分與光譜采集部分的高度以及樣品室入射光與出射光高度的一致性。為后期整機(jī)設(shè)計的內(nèi)部空間利用率考慮，光源基座的上層布置激光誘導(dǎo)熒光檢測系統(tǒng)的光源部分，下層采用中空設(shè)計，用于放置開發(fā)板?？ㄗ延瓦^氧化值檢測系統(tǒng)實(shí)物如圖2所示。

物質(zhì)的熒光強(qiáng)度會隨著照射時間的增加而下降，即存在熒光褪色效應(yīng)[21-22]?？ㄗ延蜔晒夤庾V采集時間越長，其強(qiáng)度可能越低。為減小熒光褪色效應(yīng)對實(shí)驗結(jié)果造成的不利影響，基于LabVIEW編寫了檢測系統(tǒng)熒光采集控制軟件，實(shí)現(xiàn)對光譜采集時間的控制，軟件工作流程如圖3所示。比色皿放入樣品室后，在外力作用下觸發(fā)控制部分，控制部分會發(fā)送信號給計算機(jī)開始計時并根據(jù)預(yù)先設(shè)定的時間間隔和總時長自動采集和保存光譜數(shù)據(jù)，達(dá)到設(shè)定總時長光譜采集自動停止并以.txt格式保存。通過Matlab軟件對光譜數(shù)據(jù)進(jìn)行預(yù)處理后建立模型，并將最優(yōu)過氧化值模型以.m文件格式導(dǎo)出。當(dāng)對未知葵花籽油樣品進(jìn)行預(yù)測時，利用LabVIEW調(diào)用最優(yōu)模型文件進(jìn)行實(shí)時過氧化值預(yù)測。

為確保系統(tǒng)可靠性，樣品室中放入空白比色皿采集了光譜響應(yīng)，如圖4a所示。僅在光源365 nm處存在峰值，說明背景光譜無其他信號干擾?；谧孕写罱ǖ募す庹T導(dǎo)熒光檢測系統(tǒng)采集不同氧化程度的葵花籽油熒光光譜如圖4b所示?？梢钥吹?，當(dāng)葵花籽油受到365 nm紫外光照射時，在440、495、520、680 nm附近產(chǎn)生熒光特征峰，且隨著過氧化值增加，其熒光強(qiáng)度呈現(xiàn)遞減趨勢。在400～600 nm范圍內(nèi)熒光強(qiáng)度的衰減最顯著，其中520 nm附近的熒光特征峰應(yīng)該是由葵花籽油中的維生素E產(chǎn)生的[23]。680 nm附近葉綠素的熒光光譜強(qiáng)度隨氧化程度的增加而降低，即葵花籽油中的葉綠素也會隨著氧化的進(jìn)行而被破壞。以上說明，基于自行搭建的激光誘導(dǎo)熒光檢測系統(tǒng)能夠較好地反映熒光光譜強(qiáng)度隨葵花籽油過氧化值的變化，而且完全可以檢測到國家標(biāo)準(zhǔn)限定值0.25 g/(100 g)以下的葵花籽油過氧化值。

為探究熒光褪色效應(yīng)對葵花籽油熒光強(qiáng)度的影響，基于自行搭建的系統(tǒng)每隔10 s自動采集過氧化值為0.12 g/(100 g)的葵花籽油熒光光譜，共連續(xù)采集熒光光譜7條。葵花籽油激光誘導(dǎo)熒光光譜隨時間的熒光響應(yīng)如圖5所示。1 s時在430～468 nm和480～520 nm范圍內(nèi)的光譜強(qiáng)度明顯高于10 s至60 s時的光譜強(qiáng)度。因此，在后續(xù)的研究中均將熒光光譜采集時間設(shè)定為1 s。

3 結(jié)果與討論

3.1 葵花籽油過氧化值理化分析

通過不同氧化處理制備得到80個葵花籽油樣品并測定其過氧化值標(biāo)準(zhǔn)理化值。80個樣品過氧化值范圍為0.03～1.19 g/(100 g)，平均值為0.39 g/(100 g)，中位數(shù)為0.30 g/(100 g)。

3.2 葵花籽油熒光光譜預(yù)處理及分析

采集光譜前先將自行搭建的食用油激光誘導(dǎo)熒光檢測系統(tǒng)預(yù)熱30 min，積分時間設(shè)定為100 ms，設(shè)置光譜采集時間為1 s，每個樣品采集3次光譜，取3次的平均值作為該樣品的原始光譜。80個葵花籽油樣品的原始光譜如圖6a所示。

因原始熒光光譜噪聲較大，采用SG(Savitzky-Golay)濾波對其進(jìn)行平滑?？ㄗ延偷恼扯葧S著氧化程度的加深而逐漸增加，可能導(dǎo)致各個樣品之間的光程產(chǎn)生差異。此外，不同葵花籽油樣品的光譜基線會發(fā)生不同程度的漂移。為減小上述干擾對建模結(jié)果造成的不利影響，分別對原始熒光光譜進(jìn)行多元散射校正(Multiplication scattering correction, MSC)、光路長度估計與校正[24](Optical path-length estimation and correction, OPLEC)、標(biāo)準(zhǔn)正態(tài)變量變換(Standard normal variable transform, SNV)預(yù)處理。MSC預(yù)處理以所有樣品光譜的平均光譜作為理想值，將每個光譜相對于平均光譜的變化分為偏移(斜率)變化和平移(截距)變化，根據(jù)此來校正原始光譜，校正后的光譜如圖6b所示?？梢钥闯觯?jīng)過MSC預(yù)處理后的光譜較好地消除了不同葵花籽油光譜的基線漂移。OPLEC利用正交投影的性質(zhì)，將光散射效應(yīng)從原始光譜中分離出來，僅保留與化學(xué)成分有關(guān)的光譜變化。從圖6c可以看出，熒光光譜在450～700 nm范圍內(nèi)被分為等間隔的幾段。SNV預(yù)處理是在各個波長點(diǎn)的光譜值滿足正態(tài)分布的前提下進(jìn)行的，對原始光譜數(shù)據(jù)進(jìn)行了標(biāo)準(zhǔn)化，處理后的熒光光譜變化趨勢與MSC類似，如圖6d所示。

利用CARS[25](競爭自適應(yīng)重加權(quán))算法對葵花籽油原始光譜的波長進(jìn)行篩選，提取標(biāo)準(zhǔn)化權(quán)重最大的前6個，其由大到小排列為440、580、495、680、520、458 nm。其中，440、680、520 nm處對應(yīng)的熒光特征物質(zhì)為共軛氫過氧化物、葉綠素和維生素E(或它的衍生物)?？梢钥闯鲈谘趸^程中，葵花籽油中共軛氫過氧化物產(chǎn)生的熒光信號與過氧化值最為相關(guān)，而葉綠素和維生素E的含量也間接反映過氧化值的變化。

3.3 葵花籽油過氧化值預(yù)測模型建立

利用SPXY(Sample set partitioning based on joint X-Y distance)算法，將80個葵花籽油樣品按3∶1劃分為校正集和預(yù)測集，經(jīng)不同的預(yù)處理光譜并結(jié)合CARS篩選波長后，選取合適的潛變量個數(shù)建立PLS回歸模型。原始光譜和3種預(yù)處理后的光譜建模得到的預(yù)測均方根誤差隨潛變量個數(shù)的變化如圖7所示。原始光譜、MSC和OPLEC預(yù)處理建模選取的最佳潛變量個數(shù)分別為10、7、9。經(jīng)SNV預(yù)處理后的熒光光譜共線性加劇，最終選取的潛變量個數(shù)為5。

表1 不同預(yù)處理方法的過氧化值模型結(jié)果Tab.1 Results of PV model built on different pre-processing methods

3.4 葵花籽油過氧化值快速檢測系統(tǒng)外部驗證

選取10個未參與建模的葵花籽油樣品對過氧化值定量預(yù)測模型進(jìn)行外部驗證。過氧化值預(yù)測值與理化測定值的R2為0.968 1，均方根誤差為0.041 1 g/(100 g)，如圖9a所示。過氧化值在0.30～0.80 g/(100 g)范圍內(nèi)的預(yù)測殘差較小，在低于國標(biāo)限0.25 g/(100 g)的殘差絕對值均在0.08 g/(100 g)以下，如圖9b所示。結(jié)果顯示，利用激光誘導(dǎo)熒光技術(shù)結(jié)合PLS回歸模型能夠以較高的精度和準(zhǔn)確性實(shí)現(xiàn)對葵花籽油過氧化值的快速檢測。

4 結(jié)束語