武夷茶中苯醚甲環(huán)唑的殘留檢測(cè)

李凌

福建林業(yè)職業(yè)技術(shù)學(xué)院,福建 南平 353000

中國(guó)是最早種植茶葉的國(guó)家[1]，茶葉不僅富含礦物質(zhì)元素和蛋白質(zhì)，而且還含有茶多糖以及茶多酚等功能性成分[2]。茶葉在生長(zhǎng)過程中極易發(fā)生不同種類的病蟲害，由這些病蟲害引起的一系列問題會(huì)使得茶葉的品質(zhì)下降，最終造成茶種植業(yè)嚴(yán)重虧損[3]?；瘜W(xué)方法是目前我國(guó)對(duì)茶葉病蟲害治理的主要手段，不僅對(duì)茶葉的品質(zhì)造成了影響，還會(huì)對(duì)生態(tài)環(huán)境造成不可修復(fù)的損傷[4]。

我國(guó)檢測(cè)茶葉農(nóng)藥殘留的方法主要有氣相色譜法（GC）、氣象色譜質(zhì)譜聯(lián)用法（GC-MS）以及液相色譜質(zhì)譜聯(lián)用法（LC-MS）等[5]，這些化學(xué)分析方法在重復(fù)性以及靈敏度方面具有較大的優(yōu)勢(shì)，但設(shè)備造價(jià)昂貴，樣品前處理要求較高等都是不可忽視的問題，這些缺點(diǎn)使得茶葉大規(guī)模檢測(cè)變得較為困難。研究設(shè)計(jì)一種精密度好、準(zhǔn)確度高的檢測(cè)方法，對(duì)于茶葉的大規(guī)模檢測(cè)有著深遠(yuǎn)的意義[6]。本文采用SERS技術(shù)快速分析茶葉中苯醚甲環(huán)唑農(nóng)藥殘留并結(jié)合快速溶劑前處理的方法[7]。通過乙腈提取茶葉，利用變壓吸附（PSA）和納米竹炭（NBC）去除茶葉中的有機(jī)酸和生物堿等熒光物質(zhì)。樣本的基底采用OTR202金納米顆粒作為增強(qiáng)物質(zhì)，采集苯醚甲環(huán)唑以及乙腈溶液的表面增強(qiáng)拉曼光譜，并用SERS檢測(cè)含苯醚甲環(huán)唑標(biāo)準(zhǔn)溶液的茶葉樣本，檢出農(nóng)藥殘留樣本的最低濃度，建立農(nóng)藥殘留的預(yù)測(cè)模型并用偏最小二乘法進(jìn)行數(shù)值計(jì)算，對(duì)測(cè)得的原始光譜進(jìn)行預(yù)處理，并對(duì)理論值進(jìn)行分析比對(duì)，最終得到了適合茶葉大規(guī)模檢測(cè)的方法。

1 預(yù)測(cè)模型的建立

1.1 偏最小二乘法的數(shù)學(xué)模型

偏最小二乘法即PLS，是一種常用的多因子分析方法[8]。在利用偏最小二乘法進(jìn)行數(shù)值計(jì)算時(shí)，通過交換因變量Y與自變量X的得分，使得被分析成分含量與主成分直接發(fā)生關(guān)聯(lián)。

交互驗(yàn)證均方根誤差即RMSECV，作為偏最小二乘法的評(píng)價(jià)指標(biāo)，常被用來衡量方法的準(zhǔn)確性。

RMSECV的計(jì)算公式如上所示，其中，n為樣本總數(shù)，yi為第i個(gè)樣本的實(shí)測(cè)值，為第i個(gè)樣本的預(yù)測(cè)值。由公式可知，實(shí)測(cè)值與預(yù)測(cè)值的差值越小即得到的RMSECV計(jì)算值越小，最后的效果越好。

預(yù)測(cè)均方根誤差即RMSEP，作為樣本預(yù)測(cè)能力的評(píng)價(jià)標(biāo)準(zhǔn)，其數(shù)值大小與預(yù)測(cè)能力緊密相關(guān)，數(shù)值越大，預(yù)測(cè)能力越低。

上式為預(yù)測(cè)均方根誤差的計(jì)算公式，其中，yi和^分別為集合中第i個(gè)樣本的實(shí)測(cè)值與預(yù)測(cè)值。

相關(guān)系數(shù)r用于考察樣本實(shí)測(cè)值與預(yù)測(cè)值之間的相關(guān)性，當(dāng)預(yù)測(cè)值與實(shí)測(cè)值之間的相關(guān)性越好時(shí)，r值越接近1。

2 試驗(yàn)方法

2.1 實(shí)驗(yàn)材料與儀器

苯醚甲環(huán)唑標(biāo)準(zhǔn)品（99%，Stanford Chemical），乙酸乙酯（色譜純，天津福晨化工公司），乙腈（色譜純，百靈威科技有限公司），PSA（北京迪馬科技有有限公司），納米竹炭（100 nm，安吉滄海精細(xì)化工有限公司），OTR202和OTR103（Opto Trace Technologies），鮮茶葉和干茶葉均取自福建省武夷山市的武夷茶。

K-Sens-532高靈敏度激光拉曼光譜儀（上海復(fù)享光學(xué)有限公司），Agilent 7890A氣相色譜質(zhì)譜聯(lián)用儀（美國(guó)安捷倫科技有限公司），JOYN-DCY-24S型水浴吹氮儀（上海喬躍電子有限公司），STLH-32低速離心機(jī)（天津美特斯試驗(yàn)機(jī)廠）。

2.2 干茶葉樣本制備

（1）將沒有農(nóng)藥殘留（陰性）的武夷茶干茶葉平均分成20份，每份重量保持在15 g左右。

（2）用色譜純的乙腈將苯醚甲環(huán)唑的標(biāo)準(zhǔn)溶液稀釋成20份不同濃度的溶液，每份溶液的體積為50 mL。濃度為0 mg/mL，編號(hào)為1；濃度為5 mg/mL，編號(hào)為2；濃度為10 mg/mL，編號(hào)為3；濃度范圍15～27 mg/mL內(nèi)間隔為3 mg/mL，編號(hào)為4～8；濃度范圍28～31 mg/mL內(nèi)間隔為1 mg/mL，編號(hào)為9～12；濃度范圍32～42 mg/mL內(nèi)間隔為2 mg/mL，編號(hào)為13～18；濃度45 mg/mL，編號(hào)為19，濃度55 mg/mL，編號(hào)為20。

（3）將配置好的溶液均勻噴灑到干茶葉上，得到20個(gè)樣本，干燥后粉碎備用。

2.3 干茶葉樣品前處理的優(yōu)化

PSA用量?jī)?yōu)化：將含有苯醚甲環(huán)唑農(nóng)殘的茶葉提取液均分為3份，分別裝到25 mg、50 mg、75 mg PSA的離心管中，劇烈振蕩5 min后，以4500 r/min的轉(zhuǎn)速離心20 min，將上清液過濾備用。

NBC用量?jī)?yōu)化：經(jīng)苯醚甲環(huán)唑浸泡過的茶葉粉碎后取上層液，分別裝到25 mg、30 mg、35 mg的離心管中，劇烈振蕩5 min，以4500 r/min的轉(zhuǎn)速離心20 min，取上層液，過濾待測(cè)。

2.4 光譜數(shù)據(jù)采集

拉曼光譜：功率250 mW，波長(zhǎng)設(shè)置為780 nm，積分時(shí)間30 s，分辨率為5 cm-1，積分5次求平均值。

液樣：樣品池中分別加入0.5 mL的OTR202、0.05 mL待測(cè)樣和0.5 mLOTR103并混合均勻，靜置15 min后采集拉曼信號(hào)。

農(nóng)藥殘留氣質(zhì)聯(lián)用：OV-101毛細(xì)管柱，程序升溫，柱子初溫80℃，保持5 min，升溫速率10℃/min，終溫250℃，保持15 min；氣化室溫度260℃；載氣為氦氣，流速1.5 mL/min，進(jìn)樣量1 μL。

3 武夷茶中殘留的苯醚甲環(huán)唑SERS檢測(cè)

3.1 苯醚甲環(huán)唑的普通拉曼光譜和表面增強(qiáng)拉曼光譜

由圖1苯醚甲環(huán)唑標(biāo)準(zhǔn)溶液與乙腈溶液的SERS及普通拉曼光譜圖可知，苯醚甲環(huán)唑與乙腈溶液的普通拉曼光譜圖即b和d，在曲線走勢(shì)及出峰位置上大致相同，說明信號(hào)峰并非農(nóng)藥殘留產(chǎn)生。苯醚甲環(huán)唑標(biāo)準(zhǔn)溶液的SERS（a）在加入表面增強(qiáng)劑OTR202以及OTR103之后，信號(hào)得到增強(qiáng)，出現(xiàn)特征峰。而乙腈溶液的SERS(c)并沒有出現(xiàn)特征峰。

圖1 拉曼光譜圖Fig.1 Raman spectrum

圖2 不同濃度苯醚甲環(huán)唑溶液的SERS圖Fig.2 SERS diagram of different difenoconazole concentrations

圖2為經(jīng)過預(yù)處理之后的不同濃度苯醚甲環(huán)唑溶液的SERS圖，由上至下濃度分別為55 mg/mL，45 mg/mL，32 mg/mL，28 mg/mL，17 mg/mL，10 mg/mL，5 mg/mL，0 mg/mL。觀察SERS信號(hào)可知，SERS信號(hào)強(qiáng)度的大小受溶液濃度大小的影響很明顯，當(dāng)溶液濃度變小時(shí)，SERS信號(hào)強(qiáng)度也變小，當(dāng)溶液濃度為55 mg/mL，45 mg/mL，32 mg/mL時(shí)，曲線都出現(xiàn)幾個(gè)明顯的特征峰，當(dāng)濃度降低時(shí)，特征峰已不再明顯，當(dāng)濃度大小為5 mg/mL時(shí)，特征峰已不再出現(xiàn)，由此可知，SERS能夠檢出苯醚甲環(huán)唑特征溶液的極限濃度為10 mg/mL。

3.2 試驗(yàn)條件的優(yōu)化分析

不同的提取劑效果不同，在檢測(cè)武夷茶中苯醚甲環(huán)唑含量時(shí)通常用到的提取劑有乙醚、乙腈以及二氯甲烷。乙醚和二氯甲烷在提取農(nóng)殘時(shí)會(huì)產(chǎn)生大量的共萃物，效果較差，而在使用乙腈作為提取劑時(shí)則不會(huì)出現(xiàn)這些情況。在使用SERS方法檢測(cè)茶葉中苯醚甲環(huán)唑的含量時(shí)，由于苯醚甲環(huán)唑是弱極性物質(zhì)，因此在提取時(shí)不需加水，防止茶葉中的水溶性物質(zhì)析出。

圖3 苯醚甲環(huán)唑提取液與標(biāo)準(zhǔn)溶液的SERS圖Fig.3 SERS diagrams of difenoconazole extraction solution and standard solution

由圖3苯醚甲環(huán)唑提取液SERS圖以及苯醚甲環(huán)唑標(biāo)準(zhǔn)溶液SERS圖可知，在苯醚甲環(huán)唑的提取液（a）中并未發(fā)現(xiàn)目標(biāo)物質(zhì)的信號(hào)峰，說明檢測(cè)結(jié)果可能被茶葉中的熒光物質(zhì)干擾，使得目標(biāo)峰并未出現(xiàn)。PSA作為目前市場(chǎng)上使用較多的凈化劑，具有較好的凈化效果，能夠?qū)悠分械挠袡C(jī)酸等物質(zhì)有效的去除。在不同PSA用量下，苯醚甲環(huán)唑的拉曼特征峰并未出現(xiàn)，說明拉曼信號(hào)受到茶葉基質(zhì)的干擾較大。隨著PSA用量由25 mg一直升到75 mg時(shí)，譜圖的基線逐漸趨于平穩(wěn)，漂移現(xiàn)象減弱，當(dāng)用量為75 mg時(shí)，其譜圖與用量為50 mg時(shí)譜圖基本一致，說明PSA的最佳用量為50 mg。

相比于普通竹炭，納米竹炭具有多孔性等特性，吸附力較強(qiáng)且表面積更大。隨著NBC用量的增加，基線漂移現(xiàn)象逐漸減弱，當(dāng)用量為25 mg及30 mg時(shí)，并未出現(xiàn)苯醚甲環(huán)唑的拉曼信號(hào)峰，而用量為35 mg時(shí)出，出現(xiàn)了苯醚甲環(huán)唑的特征峰。結(jié)果表明，NBC的最佳用量為35 mg。

3.3 茶葉中苯醚甲環(huán)唑農(nóng)藥殘留的SERS檢測(cè)分析與預(yù)測(cè)

取經(jīng)過處理過的茶葉樣品凈化液，放置到氮吹儀中經(jīng)氮吹后用乙腈定容，定容完成后用氣質(zhì)聯(lián)用儀進(jìn)行真實(shí)值檢測(cè)，如表1茶葉中苯醚甲環(huán)唑含量氣質(zhì)聯(lián)用檢測(cè)結(jié)果所示。

表1 茶葉中苯醚甲環(huán)唑含量氣質(zhì)聯(lián)用檢測(cè)結(jié)果Table 1 GC-MS test results of difenoconazole in tea

在對(duì)采集到的光譜數(shù)據(jù)進(jìn)行分析時(shí)，首先采用MSC（多元散射校正）、SNV（標(biāo)準(zhǔn)正態(tài)變量變換）以及一階二階導(dǎo)數(shù)進(jìn)行預(yù)處理，在對(duì)處理之后的數(shù)據(jù)建立偏最小二乘（PLS）模型。由表2不同預(yù)處理方法建立的PLS模型結(jié)果可知，經(jīng)過MSC和SNV預(yù)處理后所建的模型優(yōu)于一階二階導(dǎo)數(shù)，而一二階導(dǎo)數(shù)處理后建立的模型又優(yōu)于原始光譜，且SNV預(yù)處理方法效果較好。

表2 不同預(yù)處理方法建立的PLS模型結(jié)果Table 2 Results of PLS models built by different pretreatment methods

在驗(yàn)證方法的準(zhǔn)確性時(shí)要對(duì)不同濃度的苯醚甲環(huán)唑茶葉樣本進(jìn)行檢測(cè)，茶葉樣本經(jīng)前處理之后，采集三次拉曼信號(hào)并求得平均值，經(jīng)SNV預(yù)處理后，采用PLS模型進(jìn)行樣本預(yù)測(cè)并與農(nóng)藥殘留氣質(zhì)結(jié)果相比較，如圖3預(yù)測(cè)值與真實(shí)值所示，發(fā)現(xiàn)茶葉樣本中農(nóng)藥殘留的預(yù)測(cè)回收率在90～105%之間，真實(shí)值與預(yù)測(cè)值誤差在-10～5%之間，誤差值在可允許的范圍之內(nèi)，與化學(xué)分析方法結(jié)果相差不大，說明茶葉中苯醚甲環(huán)唑的含量檢測(cè)也可以通過該方法來完成。

表3 預(yù)測(cè)值和真實(shí)值Table 3 Predicted and true values

對(duì)樣本進(jìn)行t檢驗(yàn)，得到的結(jié)果如表4預(yù)測(cè)值與參考值樣本t檢驗(yàn)結(jié)果所示，t取值為-0.435，取其絕對(duì)值發(fā)現(xiàn)遠(yuǎn)遠(yuǎn)小于t0.05,3=3.2565，則表明建立的方法可行，預(yù)測(cè)值與參考值之間相差不大。

表4 預(yù)測(cè)值與參考值樣本t檢驗(yàn)結(jié)果Table 4 Results of forecast and reference values t test

4 結(jié)論

本文利用表面增強(qiáng)拉曼光譜（SERS）技術(shù)，探究了苯醚甲環(huán)唑農(nóng)藥殘留在武夷茶中含量的光譜檢測(cè)方法。將苯醚甲環(huán)唑作為目標(biāo)農(nóng)藥，將農(nóng)藥的標(biāo)準(zhǔn)溶液噴灑在茶葉樣本上，分析苯醚甲環(huán)唑標(biāo)準(zhǔn)溶液SERS的拉曼光譜峰。以乙腈為提取劑，優(yōu)化NBC和PSA的量，尋找到最佳的用量。以苯醚甲環(huán)唑作為目標(biāo)物質(zhì)采用偏最小二乘法建立武夷茶中農(nóng)藥殘留的預(yù)測(cè)模型，模型預(yù)測(cè)回收率為92%～105%，相對(duì)誤差最大值為5%，在允許的誤差范圍之內(nèi)。最終的t檢驗(yàn)結(jié)果表明利用偏最小二乘法建立的預(yù)測(cè)模型能夠很好的預(yù)測(cè)目標(biāo)農(nóng)藥的含量，與參考值之間的差值可忽略，滿足大規(guī)模苯醚甲環(huán)唑農(nóng)藥殘留的檢測(cè)要求。對(duì)苯醚甲環(huán)唑標(biāo)準(zhǔn)溶液的最低檢出濃度為10 mg/mL，農(nóng)藥殘留量滿足國(guó)家規(guī)定的殘留限量(GB 2763-2014，10 mg/mL)要求。

[1]文 洋,譚 君,譚 斌,等.茶葉中農(nóng)藥殘留檢測(cè)方法研究進(jìn)展[J].食品安全質(zhì)量檢測(cè)學(xué)報(bào),2017,8(3):826-831

[2]孫玉珍,周 群.液相色譜法對(duì)茶葉中農(nóng)藥殘留量的檢測(cè)方法研究[J].福建茶葉,2016,38(11):10-11

[3]曹培啟.茶葉中農(nóng)藥殘留量的氣相色譜-質(zhì)譜聯(lián)用檢測(cè)方法[J].食品安全導(dǎo)刊,2016(33):121

[4]王書濤,苑媛媛,王玉田,等.基于三維熒光與GA-RBF神經(jīng)網(wǎng)絡(luò)對(duì)茶葉中氯菊酯農(nóng)藥殘留的檢測(cè)[J].發(fā)光學(xué)報(bào),2016,37(10):1267-1274

[5]黃錦殷.茶葉農(nóng)藥殘留檢測(cè)的前處理新技術(shù)進(jìn)展[J].現(xiàn)代食品,2016,7(13):41-42

[6]吳民富.茶葉農(nóng)藥殘留檢測(cè)前處理新技術(shù)研究進(jìn)展[J].廣東茶業(yè),2015(3):12-14

[7]陸小磊,周衛(wèi)龍.中國(guó)茶葉農(nóng)藥殘留檢測(cè)方法標(biāo)準(zhǔn)概述[J].食品安全質(zhì)量檢測(cè)學(xué)報(bào),2015,6(5):1548-1553

[8]宋秀華,李傳榮,許景偉,等.靜態(tài)頂空與固相微萃取法測(cè)定元寶楓葉片揮發(fā)物成分比較[J].山東農(nóng)業(yè)大學(xué)學(xué)報(bào):自然科學(xué)版,2014,45(5):670-674