鄒德智，李 丹，張 青，于坤弘，于殿宇,*，王立琦*

（1.東北農(nóng)業(yè)大學(xué)食品學(xué)院，黑龍江 哈爾濱 150030；2.黑龍江北大荒豐威食品有限公司，黑龍江 哈爾濱 150060；3.哈爾濱商業(yè)大學(xué)計(jì)算機(jī)與信息工程學(xué)院，黑龍江 哈爾濱 150028）

太赫茲光譜法檢測(cè)面粉中過(guò)氧化苯甲酰含量

鄒德智1，李 丹1，張 青1，于坤弘2，于殿宇1,*，王立琦3,*

（1.東北農(nóng)業(yè)大學(xué)食品學(xué)院，黑龍江 哈爾濱 150030；2.黑龍江北大荒豐威食品有限公司，黑龍江 哈爾濱 150060；3.哈爾濱商業(yè)大學(xué)計(jì)算機(jī)與信息工程學(xué)院，黑龍江 哈爾濱 150028）

檢測(cè)40 個(gè)不同過(guò)氧化苯甲酰含量的面粉樣品的太赫茲時(shí)域譜，獲得其在0.2～2.0 THz頻段內(nèi)的頻域譜、折射譜及吸收譜。時(shí)域譜表明，與參考信號(hào)相比，添加過(guò)氧化苯甲酰的面粉樣品，其太赫茲波幅值出現(xiàn)了不同程度的降低，面粉樣品的頻域信號(hào)出現(xiàn)了不同程度的延遲現(xiàn)象。頻域譜表明面粉樣品中過(guò)氧化苯甲酰的含量越高，其頻域譜的幅值越大，對(duì)太赫茲波的吸收強(qiáng)度越小，過(guò)氧化苯甲酰在一定程度上降低了面粉樣品對(duì)太赫茲波的吸收。折射譜表明隨著頻率的增加，折射率逐漸降低，隨著面粉中過(guò)氧化苯甲酰的含量增加，面粉樣品的折射率隨之增大。吸收譜表明面粉樣品中的過(guò)氧化苯甲酰在1.6 THz處存在一個(gè)較為明顯的特征吸收峰。因此，太赫茲光譜可用于面粉中過(guò)氧化苯甲酰的特征識(shí)別。利用40 個(gè)面粉樣品的吸收譜數(shù)據(jù)作為輸入變量，分別采用偏最小二乘（partial least squares，PLS）、支持向量回歸（support vector regression，SVR）及最小二乘支持向量回歸（least squares support vector regression，LS-SVR）建立校正模型，其中PLS和SVR模型得到的決定系數(shù)（R2）都低于0.9，相對(duì)標(biāo)準(zhǔn)偏差（relative standard deviation，RSD）都大于10%，而LS-SVR模型的R2為0.986，RSD為4.30%，結(jié)果表明太赫茲光譜分析技術(shù)可以用于檢測(cè)面粉中過(guò)氧化苯甲酰含量。

太赫茲波譜；過(guò)氧化苯甲酰；定性；最小二乘支持向量回歸；定量；快速無(wú)損檢測(cè)

太赫茲（terahertz，THz）是指頻率在0.1～10 THz之間的電磁波，THz波譜中包含了物質(zhì)中豐富的物理和化學(xué)信息，研究樣品的THz光譜可以得到它的復(fù)介電常數(shù)、吸收系數(shù)和折射率等參數(shù)[1-5]。由于THz脈沖的波長(zhǎng)相對(duì)較長(zhǎng)，其光子能量很小，測(cè)量時(shí)不會(huì)引起有害的電離反應(yīng)，可以對(duì)樣品進(jìn)行快速無(wú)損的檢測(cè)[6-12]。近20年間，太赫茲時(shí)域光譜（terahertz tim e dom ain spectroscopy，THz-TDS）技術(shù)在食品[13-14]、中藥[15-16]、生物科學(xué)[17]等方面都有很廣泛的應(yīng)用，在今后的研究中會(huì)與X射線、近紅外等檢測(cè)分析技術(shù)互補(bǔ)結(jié)合[18-19]。

過(guò)氧化苯甲酰是一種強(qiáng)氧化劑，它能氧化類胡蘿卜素、葉黃素等色素而提高面粉白度，促進(jìn)面粉熟化，其還原產(chǎn)物苯甲酸還可以抑制蛋白酶的活性和微生物的生長(zhǎng)，避免蛋白質(zhì)的分解，提高面食制品的貯藏性能。但是面粉中加入過(guò)氧化苯甲酰后會(huì)破壞面粉中的各種營(yíng)養(yǎng)成分，使面筋質(zhì)都受到不同程度的破壞，并且使脂肪酸氧化酸敗，過(guò)氧化苯甲酰的還原產(chǎn)物苯甲酸的攝入量過(guò)多也不利于人體健康，是限制使用的食品添加劑[20]。動(dòng)物試驗(yàn)表明，長(zhǎng)期攝入過(guò)氧化苯甲酰對(duì)試驗(yàn)鼠的肝臟、脾臟、腎臟均有較大損傷[21]。所以，快速準(zhǔn)確檢測(cè)面粉中過(guò)氧化苯甲酰的含量具有重要意義。

杜娟等[22]研究在pH 6.84的Na2HPO4-KH2PO4緩沖溶液中，以羅丹明B為光譜探針，用分光光度法測(cè)定面粉中過(guò)氧化苯甲酰，在室溫條件下，羅丹明B和過(guò)氧化苯甲酰組成的體系在550 nm波長(zhǎng)處有最大吸收，加標(biāo)回收率為97.5%～101.0%，但是小麥粉中微量的過(guò)氧化苯甲酰，在放置過(guò)程中受熱、光照或輻射等會(huì)使部分過(guò)氧化苯甲酰分子中的過(guò)氧鏈斷裂，分解生成氧氣和苯甲酸，使測(cè)定結(jié)果偏低，所以本法適用于新出廠的面粉測(cè)定。Abe-Onishi等[23]利用高效液相色譜-質(zhì)譜法分析面粉及其制品中的苯甲酸和過(guò)氧化苯甲酸的殘留量，但是該方法比較耗時(shí)且操作繁瑣。夏燚等[24]利用THz光譜的指紋識(shí)別特性檢測(cè)增白劑中吊白塊含量，獲得了在0.2～2.0 THz波段的吸收譜和折射率譜，采用偏最小二乘（partial least squares，PLS）方法對(duì)增白劑中摻雜吊白塊的含量進(jìn)行了定量分析，結(jié)果表明增白劑中的吊白塊含量與THz吸收系數(shù)具有較高的相關(guān)性。而本實(shí)驗(yàn)將深入研究利用THz光譜技術(shù)快速定性和定量檢測(cè)面粉中的過(guò)氧化苯甲酰。

本實(shí)驗(yàn)是對(duì)企業(yè)經(jīng)常在面粉中添加的增白劑（過(guò)氧化苯甲酰）的含量進(jìn)行檢測(cè)，利用THz-TDS技術(shù)，通過(guò)向面粉中添加不同量的過(guò)氧化苯甲酰，經(jīng)充分混合后壓制成40 個(gè)含有不同過(guò)氧化苯甲酰含量的面粉樣品，采用柱后衍生-高效液相色譜-熒光檢測(cè)法檢測(cè)的面粉樣品實(shí)際含量值，在室溫狀態(tài)下對(duì)40 個(gè)面粉樣品進(jìn)行光譜掃描，得到了面粉樣品的時(shí)域譜，研究通過(guò)傅里葉變換和計(jì)算得到0.2～2.0 THz頻段內(nèi)的頻域譜、折射譜和吸收譜。根據(jù)面粉樣品在THz頻段內(nèi)的吸收峰特征，對(duì)面粉樣品中的過(guò)氧化苯甲酰進(jìn)行定性分析，將40 個(gè)面粉樣品的光譜數(shù)據(jù)作為輸入變量建立PLS、支持向量回歸（support vector regression，SVR）及最小二乘支持向量回歸（least squares support vector regression，LS-SVR）模型[25-27]，選擇相對(duì)標(biāo)準(zhǔn)偏差（relative standard deviation，RSD）較小，且R2較大的模型預(yù)測(cè)面粉樣品中過(guò)氧化苯甲酰的含量，研究其與柱后衍生-高效液相色譜-熒光檢測(cè)法檢測(cè)值之間的偏差，以期為控制面粉中違規(guī)添加劑的質(zhì)量安全控制提供了一種快速無(wú)損的檢測(cè)方法。

1 材料與方法

1.1 材料與試劑

面粉：從黑龍江北大荒豐威食品有限公司扦取沒(méi)有添加任何添加劑的特一粉。

過(guò)氧化苯甲酰（化學(xué)純）、H3PO4（純度85%）、氨水（純度25%）、NH4C l、冰乙酸（均為分析純）國(guó)藥集團(tuán)化學(xué)試劑有限公司；甲醇、乙醇（均為色譜純）迪馬科技公司；氯化血紅素（分析純） 美國(guó)Sigm a公司；對(duì)羥基苯乙酸（純度98%） 美國(guó)Acros公司。

1.2 儀器與設(shè)備

FW-4型壓片機(jī) 天津天光光學(xué)儀器有限公司；QL-901旋渦混合器 海門市麒麟醫(yī)用儀器廠；R968電子數(shù)顯外徑千分尺 上海量具刃具廠；Z-2型THz-TDS系統(tǒng) 美國(guó)Zomega公司；A lliance e2695型高效液相色譜儀（配有2475型熒光檢測(cè)器） 美國(guó)Waters公司；PINNACLE PCX型柱后衍生器 美國(guó)Pickering公司；S40K型pH計(jì) 瑞士Mettler公司。

1.3 方法

1.3.1 面粉樣品制備

稱取10 kg面粉，添加一定量的過(guò)氧化苯甲酰，充分?jǐn)嚢枋蛊浠旌暇鶆颍苽涑鲆环莶缓羞^(guò)氧化苯甲酰的空白樣品和添加不同量的過(guò)氧化苯甲酰面粉大樣品共40 份，并編排成1～40號(hào)面粉大樣品。

分別從40 個(gè)面粉大樣品中稱取一定量的樣品，利用柱后衍生-高效液相色譜-熒光檢測(cè)法測(cè)定面粉樣品中過(guò)氧化苯甲酰的含量[28]，分別從40 個(gè)面粉大樣品中扦取0.15 g的樣品與對(duì)THz波無(wú)吸收作用的聚乙烯粉0.05 g混合均勻，用5 t的重力壓成直徑為13 mm，厚度為1 mm的質(zhì)地均勻、兩端面平行的面粉薄片，編成1～40號(hào)面粉樣品。

1.3.2 面粉樣品的THz光譜采集及光學(xué)參數(shù)提取

在首都師范大學(xué)“太赫茲波譜與成像重點(diǎn)實(shí)驗(yàn)室”和中央民族大學(xué)理學(xué)院共同組建的“太赫茲重點(diǎn)實(shí)驗(yàn)室”中對(duì)40 個(gè)面粉樣品進(jìn)行THz光譜采集，THz-TDS系統(tǒng)其重要組成部件是美國(guó)Coherent公司的Vitesse-800-5型鈦藍(lán)寶石飛秒激光器，以及美國(guó)Zomega公司的Z-2型太赫茲系統(tǒng)；以ZnTe為探測(cè)晶體，以摻鈦藍(lán)寶石激光器作為激發(fā)光源，重復(fù)頻率為82 MHz，脈沖寬度為100 fs，中心波長(zhǎng)為800 nm，它的平均功率在960 mW。在探測(cè)過(guò)程中，系統(tǒng)信噪比為1 000，譜分辨率好于40 GHz。在采集過(guò)程中，將產(chǎn)生THz信號(hào)的砷化鎵（GaAs）、待測(cè)面粉樣品和探測(cè)晶體碲化鋅（ZnTe）組成的光路密封在充有氮?dú)獾南潴w內(nèi)，為了防止空氣中水蒸氣對(duì)THz信號(hào)的影響，需要將水分對(duì)實(shí)驗(yàn)的影響降到最低，保持箱內(nèi)的相對(duì)濕度小于2%，溫度為21 ℃，產(chǎn)生的THz脈沖經(jīng)過(guò)2 個(gè)離軸拋物面鏡聚焦后準(zhǔn)直入射到壓成的面粉樣品上，透過(guò)面粉樣品并與探測(cè)光共線后再次聚焦入射到探測(cè)元件ZnTe晶體上，產(chǎn)生線性電光效應(yīng)，并改變探測(cè)脈沖的偏振態(tài)，即可得到載有面粉樣品信息的THz脈沖電場(chǎng)的大小及變化信號(hào)，將信號(hào)送入鎖相放大器進(jìn)行放大，分別通過(guò)測(cè)量氮?dú)夂兔娣蹣悠返玫絽⒖夹盘?hào)和面粉樣品信號(hào)，取3 次測(cè)量平均值。

1.4 數(shù)據(jù)處理

將時(shí)域信號(hào)進(jìn)行快速傅里葉變換得到參考和面粉樣品的頻域信號(hào)，采用Dorney[29]和Dubillaret[30]等提出的基于菲涅爾公式的數(shù)據(jù)處理模型處理數(shù)據(jù)，得到樣品的折射率與吸收系數(shù)。

2 結(jié)果與分析

2.1 面粉樣品的制備

利用柱后衍生-高效液相色譜-熒光檢測(cè)法測(cè)得面粉樣品中過(guò)氧化苯甲酰的含量，編成1～40號(hào)面粉樣品如表1所示。

表1 40 個(gè)面粉樣品中過(guò)氧化苯甲酰的含量Tab le 1 Content of benzoyl peroxide in 40 flour sam p les

2.2 參考信號(hào)與樣品信號(hào)的時(shí)域譜

利用THz光譜系統(tǒng)掃描40 個(gè)面粉樣品，得到參考信號(hào)與樣品信號(hào)的時(shí)域譜，如圖1所示。

圖1 參考信號(hào)與樣品信號(hào)的時(shí)域譜Fig. 1 Time domain spectra of reference signal and samp le signal

從圖1可以看出，由于不同過(guò)氧化苯甲酰含量的面粉樣品對(duì)THz波的吸收程度不同，衰減的程度與延遲的時(shí)間長(zhǎng)短也有所不同，與參考信號(hào)相比，含有不同含量的過(guò)氧化苯甲酰的面粉樣品的幅值出現(xiàn)了不同程度的降低，面粉樣品信號(hào)產(chǎn)生的時(shí)間出現(xiàn)了不同程度的延遲現(xiàn)象。

2.3 參考信號(hào)與樣品信號(hào)的頻域譜

對(duì)40 個(gè)面粉樣品的時(shí)域信號(hào)進(jìn)行傅里葉變換，得到參考信號(hào)與樣品信號(hào)的頻域譜，如圖2所示。在0.2～2.0 THz波段內(nèi)，參考樣品和面粉樣本的頻域譜的波形基本相同，差別不大，參考信號(hào)的振幅最高，面粉樣品在該THz頻段的幅值均有不同程度的降低，過(guò)氧化苯甲酰的含量越高，其頻域譜的幅值越大，對(duì)THz波的吸收強(qiáng)度越小，可以看出，面粉中的過(guò)氧化苯甲酰在一定程度上降低了面粉樣品對(duì)THz波的吸收。

圖2 參考信號(hào)與樣品信號(hào)的頻域譜Fig. 2 Frequency domain spectra of reference signal and sample signal

2.4 面粉樣品的THz折射譜

將經(jīng)過(guò)傅里葉變換得到的相位差代入到菲涅爾公式的數(shù)據(jù)處理模型中，經(jīng)計(jì)算得到面粉樣品的THz折射率，40 個(gè)面粉樣品的THz折射譜如圖3所示。

圖3 面粉樣品的THz折射譜Fig. 3 Terahertz refraction spectra of fl our sam p les

由圖3可知，面粉中不同含量的過(guò)氧化苯甲酰面粉樣品的折射率譜中，在0.5～2.0 THz波段內(nèi)，隨著頻率的增加，面粉樣品的折射率呈現(xiàn)先減小再平緩增大，最后再減小的趨勢(shì)。過(guò)氧化苯甲酰面粉樣品的折射率范圍為1.40～1.75，平均折射率為1.58。從整體上看，隨著面粉中過(guò)氧化苯甲酰的含量增加，面粉樣品的折射率隨之增大，面粉樣品的折射譜的變化趨勢(shì)基本相同，均隨頻率的增大而減小。因?yàn)槊娣鄣臉悠沸盘?hào)都在它的平均折射率附近變化，沒(méi)有出現(xiàn)較為明顯的吸收峰與反常色散現(xiàn)象，這可能是因?yàn)槊娣蹣悠贩肿娱g的相互作用比較弱。

2.5 面粉樣品的THz吸收譜

將頻域譜中的幅值和計(jì)算得到的折射率代入到菲涅爾公式的數(shù)據(jù)處理模型中，經(jīng)計(jì)算得到面粉樣品的THz吸收系數(shù)，40 個(gè)面粉樣品的THz吸收譜如圖4所示。

由圖4可知，在0.4～2.0 THz頻域波段范圍內(nèi)，出現(xiàn)了較顯著的吸收峰，吸收系數(shù)在0～1.2 cm-1范圍內(nèi)呈非線性增長(zhǎng)的趨勢(shì)，由于面粉樣品的散射，含有不同過(guò)氧化苯甲酰含量的面粉樣品的吸收系數(shù)均隨頻率的增加而增加。且過(guò)氧化苯甲酰的含量越高，衰減的程度越低，吸收系數(shù)越小，說(shuō)明過(guò)氧化苯甲酰的存在降低了面粉樣品對(duì)THz波的吸收程度。面粉樣品在1.6 THz處都出現(xiàn)了第1個(gè)較為明顯的特征吸收峰，可作為其在THz波段的指紋譜用于物質(zhì)識(shí)別，面粉樣品的吸收譜存在部分交叉的現(xiàn)象，在1.6 THz后面粉樣品的色散現(xiàn)象也相對(duì)減弱。

圖4 面粉樣品的THz吸收譜Fig. 4 Terahertz absorp tion spectra of fl our sam p les

2.6 面粉樣品的定量分析

利用THz吸收譜數(shù)據(jù)和面粉樣品中過(guò)氧化苯甲酰的含量分別建立校正模型，從40 個(gè)面粉樣本中按一定間隔抽取30 個(gè)面粉樣品作為校正集來(lái)建立模型，分別采用PLS、SVR及LS-SVR 3 種方法建模，用其余10 個(gè)面粉樣品作為預(yù)測(cè)集來(lái)驗(yàn)證模型，模型的準(zhǔn)確性由R2、預(yù)測(cè)均方根誤差（root mean square error of prediction，RMSEP）及RSD 3個(gè)參數(shù)判定，預(yù)測(cè)結(jié)果見(jiàn)表2。

表2 基于不同模型的預(yù)測(cè)結(jié)果Tab le 2 Prediction results based on different models

由表2可以看出，通過(guò)與用高效液相色譜-熒光檢測(cè)法測(cè)定的面粉樣品實(shí)際值相比，PLS和SVR的R2都低于0.9，擬合度較低，RSD都大于10%，偏差較大。而LS-SVR的R2為0.986，表明面粉樣品中過(guò)氧化苯甲酰的含量與THz吸收系數(shù)具有較高的相關(guān)性，RMSEP為5.64，RSD為4.30%，低于5%，其準(zhǔn)確性較好，明顯低于前兩種模型的偏差，說(shuō)明模型的預(yù)測(cè)效果較好。

2.7 面粉樣品的實(shí)際值與預(yù)測(cè)值

從圖5可以看出，用預(yù)測(cè)效果較好的LS-SVR方法建模，得到混合物中吊白塊的實(shí)際值與LS-SVR得到的預(yù)測(cè)值的線性擬合結(jié)果，可以看出10 個(gè)預(yù)測(cè)集的面粉樣品在實(shí)際值標(biāo)準(zhǔn)直線附近輕微波動(dòng)，預(yù)測(cè)值與實(shí)際值基本一致，表明預(yù)測(cè)值與實(shí)際值擬合度較高，定量分析產(chǎn)生的誤差主要來(lái)源于面粉樣品的濃度配制誤差、壓片的厚度均勻性不一致和THz系統(tǒng)的信噪比。

圖5 LS-SVR模型預(yù)測(cè)值與實(shí)際值Fig. 5 Linear fi tting of actual values versus predicted values from the LS-SVR model

3 結(jié) 論

利用THz光譜快速檢測(cè)了不同過(guò)氧化苯甲酰含量的面粉，獲得面粉樣品的THz時(shí)域譜、頻域譜、折射譜及吸收譜。在0.2～2.0 THz頻段內(nèi)，在1.6 THz處有明顯的特征吸收峰，可用于面粉中過(guò)氧化苯甲酰的特征識(shí)別，以及面粉中過(guò)氧化苯甲酰的定性分析。在時(shí)域譜中，與參考信號(hào)相比，不同含量的過(guò)氧化苯甲酰面粉樣品幅值出現(xiàn)了不同程度的降低；隨著面粉樣品中過(guò)氧化苯甲酰的含量越高，頻域譜的幅值越大，折射率越大，吸收系數(shù)越小，利用PLS、SVR及LS-SVR建立過(guò)氧化苯甲酰含量與面粉樣品吸收系數(shù)之間的定量分析模型，其中PLS和SVR的R2都較小，RMSEP和RSD都較大。而LS-SVR的R2較大，RMSEP較小，RSD低于5%，能夠更加準(zhǔn)確地預(yù)測(cè)面粉樣品中過(guò)氧化苯甲酰的含量，由于面粉中的添加劑種類較多，對(duì)檢測(cè)的過(guò)程具有一定的干擾，因此首先使用THz時(shí)域光譜技術(shù)對(duì)面粉中的過(guò)氧化苯甲酰進(jìn)行定性分析及定量檢測(cè)，為快速無(wú)損檢測(cè)面粉中違規(guī)添加劑提供理論基礎(chǔ)。

[1] 劉銳, 顧春明, 賀莉蓉, 等. ZnTe晶體中光整流產(chǎn)生的THz輻射及其電光探測(cè)研究[J]. 物理學(xué)報(bào), 2004, 53(4): 1217-1222. DOI:10.3321/j.issn:1000-3290.2004.04.045.

[2] 胡穎, 張存林, 沈京玲, 等. (100)M gO和(100)LaA lO3高溫超導(dǎo)基片材料THz時(shí)域光譜研究[J]. 物理學(xué)報(bào), 2004, 53(6): 156-160.DOI:10.3321/j.issn:1000-3290.2004.06.030.

[3] 馬士華, 施宇蕾, 嚴(yán)偉, 等. 酪氨酸的太赫茲頻譜研究[J]. 光譜學(xué)與光譜分析, 2007, 27(9): 1665-1668.

[4] 涂閃, 張文濤, 熊顯名, 等. 基于太赫茲時(shí)域光譜系統(tǒng)的轉(zhuǎn)基因棉花種子主成分特性分析[J]. 光子學(xué)報(bào), 2015, 44(4): 176-181.DOI:10.3788/gzxb20154404.0430001.

[5] 曹丙花. 基于太赫茲時(shí)域光譜的檢測(cè)技術(shù)研究[D]. 杭州: 浙江大學(xué), 2009.

[6] 沈飛, 應(yīng)義斌. 太赫茲光譜和成像技術(shù)在食品安全檢測(cè)中的應(yīng)用[J].光譜學(xué)與光譜分析, 2009, 29(6): 1445-1449. DOI:10.3964/j.is sn.1000-0593(2009)06-1445-05.

[7] 王迎新, 康克軍, 陳志強(qiáng), 等. 利用太赫茲光譜定量分析固體混合物[J].清華大學(xué)學(xué)報(bào)(自然科學(xué)版), 2009, 49(2): 161-164. DOI:10.3321/j.is sn:1000-0054.2009.02.001.

[8] 郭瀾濤, 牧凱軍, 鄧朝, 等. 太赫茲波譜與成像技術(shù)[J]. 紅外與激光工程, 2013, 42(1): 51-56. DOI:10.3969/j.issn.1007-2276.2013.01.010.

[9] MARKELZ A G, KNAB J R, JING Y C, et al. Protein dynam ical transition in terahertz dielectric response[J]. Chemical Physics Letters,2007, 442(4/5/6): 413-417. DOI:10.1016/j.cplett.2007.05.080.

[10] STRACHAN C J, RADES T, NEWNHAM D A, et al. Using terahertz pulsed spectroscopy to study crystallinity of pharmaceutical m aterials[J]. Chem ical Physics Letters, 2004, 390: 20-24.DOI:10.1016/j.cplett.2004.03.117.

[11] TAKAHASHI M, ISHIKAWA Y, NISHIZAWA J I, et al. Lowfrequency vibrational modes of ribofl avin and related compounds[J].Chem ical Physics Letters, 2005, 401(4): 475-482. DOI:10.1016/j.cplett.2004.11.108.

[12] XIANG W F, WANG X, LIU Y, et al. Density detection of aligned nanow ire arrays using terahertz time-domain spectroscopy[J].Nanoscale Research Letters, 2016, 11(1): 337-341. DOI:10.1186/s11671-016-1551-1.

[13] YANG C N, TIAN L, ZHAO K N. Spectroscopic studies on the edible flavoring in terahertz range[J]. Acta Photonica Sinica, 2012, 41(5):627-630. DOI:10.3788/gzxb20124105.0627

[14] 朱莉. 基于太赫茲波的食品添加劑檢測(cè)技術(shù)研究[D]. 杭州: 浙江大學(xué), 2008.

[15] 周永軍, 劉勁松, 王可嘉, 等. 基于太赫茲譜分析中藥材鑒別[J].光譜學(xué)與光譜分析, 2014, 34(7): 1840-1843. DOI:10.3964/j.is sn.1000-0593(2014)07-1840-04.

[16] 張平. 中草藥的太赫茲光譜鑒別[D]. 北京: 首都師范大學(xué), 2008.

[17] 劉盛綱, 鐘任斌. 太赫茲科學(xué)技術(shù)及其應(yīng)用的新發(fā)展[J].電子科技大學(xué)學(xué)報(bào), 2009, 38(5): 481-486. DOI:10.3969/j.issn.1001-0548.2009.05.001.

[18] 李明亮, 常天英, 崔洪亮. 基于太赫茲時(shí)域光譜分析紅層泥巖風(fēng)化樣本[J]. 光譜學(xué)與光譜分析, 2016, 36(4): 929-933. DOI:10.3964/j.is sn.1000-0593(2016)04-0929-05.

[19] 趙容嬌, 何金龍, 李璟, 等. L-和DL-福多司坦的太赫茲光譜分析[J].物理化學(xué)學(xué)報(bào), 2011, 27(12): 2743-2748. DOI:10.3866/PKU.WHXB20112743.

[20] 張志勇. 基于近紅外光譜法檢測(cè)面粉中過(guò)氧化苯甲酰添加量[D]. 天津:天津大學(xué), 2012.

[21] JIA X, WU Y, LIU P. Effects of flour bleaching agent on m ice liver antioxidant status and ATPases[J]. Environmental Toxicology and Pharmaco logy, 2011, 31(3): 479-484. DOI:10.1016/j.etap.2011.03.009.

[22] 杜娟, 孫丹. 羅丹明B光譜探針測(cè)定面粉中過(guò)氧化苯甲酰[J]. 食品科學(xué), 2009, 30(10): 156-158. DOI:10.3321/j.issn:1002-6630.2009.10.031.

[23] ABE-ONISHI Y, YOMOTA C, SUGIMOTO N, et al. Determ ination o f benzoy l peroxide and benzoic acid in w heat flour by highperformance liquid chromatography and its identification by highperformance liquid chromatography-mass spectrometry[J]. Journal o f Chromatography A, 2004, 1040(2): 209-214. DOI:10.1016/j.chroma.2004.03.059.

[24] 夏燚, 杜勇, 張慧麗, 等. 增白劑中吊白塊含量的太赫茲光譜定性與定量檢測(cè)[J]. 中國(guó)糧油學(xué)報(bào), 2015, 30(2): 103-106.

[25] 于雷, 洪永勝, 耿雷, 等. 基于偏最小二乘回歸的土壤有機(jī)質(zhì)含量高光譜估算[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2015, 31(14): 103-109. DOI:10.11975/j.issn.1002-6819.2015.14.015.

[26] 任順, 于海業(yè), 周麗娜, 等. 基于支持向量機(jī)的葉綠素?zé)晒忸A(yù)測(cè)光能利用效率研究[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2015, 46(4): 273-276.DOI:10.6041/j.issn.1000-1298.2015.04.040.

[27] 畢艷蘭, 任小娜, 彭丹, 等. 粒子群最小二乘支持向量機(jī)結(jié)合偏最小二乘法用于芝麻油質(zhì)量的鑒別[J]. 分析化學(xué), 2013, 41(9): 1366-1372. DOI:10.3724/SP.J.1096.2013.21110.

[28] 馬明, 肖玲艷, 盛亞紅, 等. 柱后衍生-高效液相色譜-熒光檢測(cè)法測(cè)定面粉中的過(guò)氧化苯甲酰[J]. 分析化學(xué), 2012, 40(3): 466-469.DOI:10.3724/SP.J.1096.2012.10835.

[29] DORNEY T D, BARANIUK R G, M ITTLEMAN D M. Material parameter estimation w ith terahertz time-domain spectroscopy[J].Journal of the Optical Society of America A Optics Image Science and Vision, 2001, 18(7): 1562-1571. DOI:10.1364/JOSAA.18.001562.

[30] DUV ILLARET L, GARET F, COUTAZ J L. H igh ly p recise determ ination of optical constants and sample thickness in terahertz time-domain spectroscopy[J]. Applied Optics, 1999, 38(2): 409-415.DOI:10.1364/AO.38.000409.

Determ ination of the Content of Benzoyl Peroxide in Flour Based on Terahertz Spectroscopy

ZOU Dezhi1, LI Dan1, ZHANG Qing1, YU Kunhong2, YU Dianyu1,*, WANG Liqi3,*
(1. College of Food Science, Northeast Agricultural University, Harbin 150030, China;2. Heilongjiang Beidahuang Fengwei Food Co. Ltd., Harbin 150060, China;3. School of Computer and Information Engineering, Harbin University of Commerce, Harbin 150028, China)

Terahertz time domain spectra for 40 flour samples w ith different contents of benzoyl peroxide were measured,and frequency domain spectra, refractive spectra and absorption spectra were obtained in the range of 0.2–2.0 THz. The time domain spectra indicated that the terahertz wave amplitude value of the flour samples w ith benzoyl peroxide was reduced as compared to reference signal, and the frequency domain signal was delayed to different degrees. The frequency domain spectra indicated larger amplitude and lower intensity of terahertz wave absorption for flour samples w ith higher content of benzoyl peroxide, demonstrating that benzoyl peroxide could reduce terahertz wave absorption in flour. The refractive spectra indicated that refractive index decreased gradually w ith the increase in frequency, but it increased w ith increasing content of benzoyl peroxide in flour. The absorption spectra indicated that benzoyl peroxide in flour samples had a significant characteristic absorption peak at 1.6 THz. Therefore, terahertz spectrum can be used for the identification of benzoyl peroxide in flour. Partial least squares (PLS), support vector regression (SVR) and least squares support vector regression (LS-SVR) were used separately to establish a calibration model by using the absorption spectral data of 40 flour samples as input variables. The determ ination coefficients (R2) of both the PLS and SVR models were lower than 0.9, w ith relative standard deviations (RSD) greater than 10%, while the R2of the LS-SVR model was 0.986, w ith RSD of 4.30%,indicating that terahertz spectroscopy can be used to detect the content of benzoyl peroxide in flour.

terahertz spectroscopy; benzoyl peroxide; qualitative; least squares support vector regression; quantitative; rapid nondestructive detection

2017-05-31

“十三五”國(guó)家重點(diǎn)研發(fā)計(jì)劃重點(diǎn)專項(xiàng)（2016YD0401402）；國(guó)家自然科學(xué)基金面上項(xiàng)目（31571880）

鄒德智（1992—），男，碩士研究生，研究方向?yàn)榧Z食、油脂及植物蛋白工程。E-mail：1393600121@qq.com

*通信作者：于殿宇（1964—），男，教授，博士，研究方向?yàn)榧Z油精深加工技術(shù)。E-mail：dyyu2000@126.com

王立琦（1966—），女，教授，博士，研究方向?yàn)槭称钒踩焖贆z測(cè)。E-mail：hsdw lq@163.com

10.7506/spkx1002-6630-201722044

TS211.7

A

1002-6630（2017）22-0298-05

鄒德智, 李丹, 張青, 等. 太赫茲光譜法檢測(cè)面粉中過(guò)氧化苯甲酰含量[J]. 食品科學(xué), 2017, 38(22): 298-302.

DOI:10.7506/spkx1002-6630-201722044. http://www.spkx.net.cn

ZOU Dezhi, LI Dan, ZHANG Qing, et al. Determ ination of the content of benzoyl peroxide in flour based on terahertz spectroscopy[J]. Food Science, 2017, 38(22): 298-302. (in Chinese w ith English abstract) DOI:10.7506/spkx1002-6630-201722044. http://www.spkx.net.cn