张莎 刘木华 陈金印 赵进辉
摘要：由于采后处理过程中脐橙保鲜剂抑霉唑易通过果皮渗进果肉中残留，不慎食用后会对人体产生危害。因此，本研究探索一种基于表面增强拉曼光谱技术（Surface-Enhanced Raman Spectroscopy ，SERS）的脐橙果皮中抑霉唑残留的快速检测方法。首先对SERS 检测条件进行优化，分别确定了最优的检测条件为反应时间2 min ，金胶加入量400µL ，NaBr作为电解质溶液且加入量为25µL 。基于以上最优检测条件，以自适应迭代惩罚最小二乘法（Adaptive Iterative Reweighted Penalized Least Squares ，air PLS）、air PLS+归一化、 air PLS+基线校正、air PLS+一阶导数、air PLS+标准正态变量（Standard Normal Distribution ，SNV）和 air PLS+多元散射校正（Multiplicative Scatter Correction ，MSC）处理后的6组光谱数据为研究对象，分别采用这 6种光谱预处理法建立支持向量回归（Support Vector Regression ，SVR）模型并对预测性能进行比较后发现， air PLS 方法所建立模型的预测集相关系数（Coefficient of the Determinant for the Prediction Set ，RP ）最大，预测集均方根误差（Root-Mean-Square Error of Prediction ，RMSEP）最小。对光谱数据进行主成分分析（Prin⁃cipal Component Analysis ，PCA）特征提取，选择前7个主成分得分作为 SVR 预测模型的输入值。采用 SVR、多元线性回归（Multiple Linear Regression ，MLR）和偏最小二乘回归（Partial Least Squares Regression， PLSR）三种建模方法分析比较其对应的预测性能，其中SVR 模型的预测集RP 可高达0.9156，预测集RMSEP 为4.8407 mg/kg ，相对标准偏差（Relative Standard Deviation ，RPD）为2.3103，表明基于 SVR 算法对脐橙表面抑霉唑残留的预测值越接近实测值，越能有效提高模型预测准确性。试验结果表明，利用SERS 结合PCA 及SVR 建模，可实现对脐橙果皮中抑霉唑残留的快速检测。
关键词：脐橙;抑霉唑;表面增强拉曼光谱;支持向量回归;多元线性回归;偏最小二乘回归
中图分类号：
O433.4文献标志码：
A文章编号：202106-SA002
引用格式：张莎， 刘木华， 陈金印， 赵进辉. 采用表面增强拉曼光谱技术快速检测脐橙果皮中抑霉唑残留[J].智慧农业（中英文）， 2021， 3（4）：42-52.
ZHANG Sha， LIU Muhua， CHEN Jinyin， ZHAO Jinhui. Rapid detection of imazalil residues in navel orange peel using surface-enhanced Raman spectroscopy[J]. Smart Agriculture， 2021， 3（4）：42-52.（in Chinese with English abstract）
1 引言
赣南脐橙作为中国赣州地理标志性产品，带动了当地经济的发展，深受当地政府重视，该产地也成为了全国最大的脐橙生产区[1]。脐橙在贮藏过程中易受病原菌侵害，导致其腐烂，造成经济损失[2]。目前，采用化学保鲜剂是柑橘类水果储存中病害防治相对成熟、效果较好的保鲜方式。抑霉唑（Imazalil，IMZ）和咪鲜胺等化学保鲜剂主要用来防治脐橙的青霉病和绿霉病[2，3]，但该类保鲜剂长期使用容易产生耐药性，进而引起保鲜效率的折损。不少商贩为了达到保鲜效果，在贮藏过程中加大保鲜剂的使用剂量，造成脐橙中抑霉唑过量残留。而且，即便是少量的抑霉唑残留也会使动物和人类的肝脏损伤，甚至可能致癌[4， 5]。为此，中国规定了脐橙中抑霉唑保鲜剂残留的最大限量为5 mg/kg[6]。为了确保脐橙的食用安全品质，检测脐橙果皮中抑霉唑的残留量具有重要的实际意义。
目前常用的抑霉唑残留检测方法主要有气相色谱法（Gas Chromatography ，GC）[7， 8]、高效液相色谱法（High Performance Liquid Chromatography ，HPLC）[9， 10]、质谱法（Mass Spectrometry， MS）[7， 9]和超高效液相色谱法（Ultra-Performance Liquid Chromatography，UPLC）[10]等。以上方法主要优点在于再现性好、精确性和敏感性高，但所需要的仪器设备较昂贵[11]，试验操作繁琐，耗时长，无法实现对抑霉唑残留量的快速检测，在实际应用中受到一定的限制。因此，有必要探索一种可以快速检测脐橙果皮中抑霉唑残留的方法。
表面增强拉曼光谱技术（Surface-Enhanced Raman Spectroscopy ，SERS）是一種当入射激光照射到金或银等粗糙金属表面时，在其表面会激发出等离子体并在粗糙金属粒子的覆盖区域内产生共振，使其表面产生比普通拉曼信号强4～11个数量级的信号强度的技术，该方法具有方便快捷、灵敏度高、稳定性好、操作简单和可适用于现场分析的优点[12， 13]。Zhao 等[11]采用 SERS 技术制备了柔性衬底，对苹果表面抑霉唑残留进行原位检测，在苹果上检测到的最低浓度可达0.073 mg/kg 。该研究表明，SERS方法具有较高的灵敏度。抑霉唑在水果采后保鲜的过程起到非常重要的作用。为避免由于脐橙果皮抑霉唑残留问题对人体造成伤害，采用 SERS技术检测脐橙果皮中抑霉唑残留具有重要意义。本研究主要以脐橙为对象，采用自适应迭代惩罚最小二乘法（Adaptive Iterative Reweighted Penalized LeastSquares，air PLS）对含不同浓度的抑霉唑的脐橙果皮的 SERS光谱进行预处理，并选择前7个主成分得分值作为支持向量回归（Support VectorRegression，SVR）模型的输入值，建立 SVR 预测模型对脐橙果皮中的抑霉唑残留量进行预测。
2 材料与方法
2.1 仪器与设备
便携式拉曼光谱仪系统主要包括QE65Pro型拉曼光谱仪（美国海洋光学公司）、785 nm 激光器、光纤和采样附件等; DXRTM 显微拉曼光谱仪（美国赛默飞世尔科技有限公司）;磁力搅拌器（LC-DMS-H ，上海力辰邦西仪器科技有限公司）;数控超声波清洗器（KQ-500DE 型，昆山市超声仪器有限公司）;电子天平（FA2004，上海舜宇恒平科学仪器有限公司）;实验室超纯水机（美国TeLedyne水质公司）。
2.2 材料与试剂
本研究所用脐橙购于江西农业大学水果店。研究所用试剂包括三水合四氯金酸（HAuCl4·3H2O≥49.0%，西格玛奥德里奇贸易有限公司），抑霉唑乳油（500 g/L ，江苏龙灯化学有限公司），抑霉唑标准品（IMZ ，纯度≥99%，阿拉丁），柠檬酸三钠（C6H5Na3O7·2H2O ，99%，西陇化工股份有限公司），溴化钙（CaBr2，≥98%，阿拉丁），溴化钠（NaBr，99%，阿拉丁），氯化钙（CaCl2，分析纯，纯度98%），氯化钠（Na‐Cl ，99.5%，西陇化工股份有限公司），以及实验室自制超纯水。
2.3 样本制备
（1）抑霉唑原液配制。取100µL 的抑霉唑乳油（500 g/L）置于100 mL 容量瓶中，加入62.5 mL 的超纯水均匀混合，得到800 mg/L的抑霉唑原液。
（2）抑霉唑工作溶液制备。取不同体积的乳油（500 g/L）用超纯水稀释成不同质量浓度的抑霉唑工作溶液（250～500 mg/L）。
（3）样本制备。将脐橙的果皮用超纯水洗净，晾干后将果皮切成若干个面积为1.5 cm×1.5 cm 的脐橙小块，对其进行称重（a ，单位：g）。
在条件优化试验中，将切好的脐橙小块浸泡于抑霉唑溶液（800 mg/L）中，90 s 后捞出晾干，再对其称重（b ，单位：g），得到含40.34 mg/kg 抑霉唑的脐橙果皮样本。
在定量试验中，将切好的65个脐橙小块浸泡于抑霉唑溶液（250～500 mg/L）中，90 s后捞出晾干，再对其称重（b ，单位：
g），得到含抑霉唑的脐橙果皮样本（1.43～40.81 mg/kg）。未浸泡在抑霉唑中的脐橙果皮样本为空白脐橙果皮样本。其中，在条件优化试验部分与定量试验部分中，a 为脐橙果皮的初始重量，b 为浸泡抑霉唑溶液且晾干后的重量; b-a可得到脐橙果皮中抑霉唑的重量，g;设 c 为抑霉唑溶液的浓度， mg/kg，则（b-a）×c为脐橙果皮抑霉唑实际吸收抑霉唑的重量，mg;（b-a）×c/b则为每块脐橙果皮上实际吸收的抑霉唑的濃度，mg/kg。
（4）金胶制备。取3 mL 浓度为1%的 HAuCL4溶液于150 mL 的烧杯中，再慢慢加入47 mL超纯水，混合均匀，将混合液放至磁力搅拌器中加热到沸腾。沸腾后立即加入2 mL 浓度为1%的柠檬酸三钠溶液，并搅拌8 min 后，冷却至室温形成金胶溶液备用[14， 15]。
（5）电解质溶液制备。取0.1 g CaBr2、0.55 g
CaCl2、5.14 g NaBr和2.9 g NaCl分别置于棕色的 50 mL容量瓶中，先加入少量水溶解，再定容至 50 mL ，可分别得到 0.01 mol/L 的 CaBr2溶液、0.1 mol/L CaCl2 溶液、 1 mol/L NaBr和1 mol/L NaCl溶液。
2.4 拉曼光谱采集
抑霉唑标准品拉曼光谱采集。取适量的抑霉唑标准品于载玻片上，用 DXRTM 显微拉曼光谱仪采集其拉曼光谱。本研究主要分析400～1800 cm-1 波段的拉曼光谱。
定性试验中脐橙果皮 SERS光谱采集。将400µL纳米金胶和25µL NaBr溶液充分混合后喷到含40.34 mg/kg 抑霉唑的脐橙小块上，待金胶与抑霉唑分子反应2 min后，采集其 SERS光谱。每个优化条件下，设置5个平行样，取其平均值为该条件下的拉曼光谱。
定量试验中脐橙果皮 SERS光谱采集。将400µL纳米金胶和25µL NaBr溶液充分混合后喷到含抑霉唑的脐橙小块（1.43～40.81 mg/kg）上，采集脐橙的 SERS光谱。每个浓度下的脐橙样本采集5条拉曼光谱作为平行样，取其平均值作为该样本的拉曼光谱。
其中，DXRTM显微拉曼光谱仪的参数设置为：激光能量5 mW，激发波长780nm，物镜的放大倍数选择10倍，采集曝光时间10 s ，预览采集时间10 s ，样品曝光10次，背景曝光16次。便携式拉曼光谱仪的参数设置为：积分时间为60 s ，平均次数为2，平滑度为1。
2.5 SERS检测试验方法
2.5.1 定性检测
通过比较空白脐橙果皮和含抑霉唑的脐橙果皮（40.34 mg/kg）的 SERS 光谱，得到 983和1047 cm-1 波段作为特征峰鉴别脐橙果皮中是否残留抑霉唑。因此，本研究通过比较983和1047 cm-1 波段处的 SERS信号强度来确定最佳的检测条件。
为研究不同反应时间对 SERS 强度的影响，
将400µL 金胶与25µL 的NaBr溶液混合均匀后喷洒到含40.34 mg/kg 抑霉唑的脐橙样本上，分别采集反应时间为 0、2、4、6和 8 min 时的SERS 光谱。通过比较983和 1047 cm-1 波段处的SERS信号强度与反应时间的关系来确定抑霉唑分子与金胶的最佳反应时间。
为研究不同金胶加入量对 SERS 强度的影响，分别将80、240、400、560和 720µL 的金胶与25µL 的NaBr溶液混合均匀后喷洒到含40.34 mg/kg抑霉唑的脐橙样本上，反应2 min后采集其 SERS光谱。通过比较983和 1047 cm-1 波段处的 SERS信号强度与金胶加入量的关系来确定最佳的金胶量。
为研究不同电解质溶液的加入量对 SERS强度的影响，将400µL 金胶分别与25µL 的 CaBr2溶液、CaCl2溶液、NaBr溶液和 NaCl 溶液混合均匀后喷洒到含40.34 mg/kg 抑霉唑的脐橙样本上，反应2 min 后采集其 SERS 光谱。通过比较983和 1047 cm-1 波段处的 SERS信号强度与电解质溶液（CaBr2、CaCl2、NaBr和 NaCl）的直方图来确定最佳的电解质溶液。
为了研究不同NaBr加入量对 SERS 强度的影响，将400µL 金胶分别与0 、10、15、20、25、30µL 的NaBr溶液混合均匀后喷洒到含 40.34 mg/kg抑霉唑的脐橙样本上，反应2 min后采集其 SERS光谱。通过比较983和 1047 cm-1 波段处的 SERS信号强度与NaBr加入量的折线图来确定最佳的电解质溶液。
2.5.2 定量检测
将400µL 金胶与 25µL 的NaBr溶液混合均匀后喷洒到含抑霉唑的脐橙样本（1.43～ 40.81 mg/kg）上，反应2 min后采集其拉曼光谱。得到65个含抑霉唑的脐橙样本，从中随机挑选35个建立预测模型，剩余的30个脐橙样本则用作预测集，用来评估该模型的预测性能。
2.6 抑霉唑理论拉曼光谱计算
为对比理论数据结果和试验结果之间的一致性，采用密度泛函理论（Density Functional Theory ，DFT）计算方法对抑霉唑分子的振动归属进行理论运算分析。
抑霉唑的3D分子结构采用Gauss View 5.0构造，Gaussian 09W软件计算理论抑霉唑分子的拉曼光谱频率，在计算时采用 DFT 方法中的B3LYP/6-31G （d ，p）基组。C 、H 、Cl 、N和 O原子也采用6-31G （d ，p）基组，在 Gaussian 09中显示计算结果。对拉曼位移用该机组相对应的校正系数（0.961）校正。其中基组为该电子体系内轨道的数学描述，“6”解释为内层轨道的函数数量;而一个基组中含3个高斯函数（Gauss‐ian Type Orbital ， GTO），每一个 Slater 函数（Slater Type Orbital ， STO）由 6个 GTO 展开，“31”解释为两组 STO 描述的价层的轨道;“d”和“p”分别是对分子中的C 、Cl 、N 、O原子和H原子做的极化函数[16]。
3 结果与分析
3.1 抑霉唑标准品的拉曼光谱分析
抑霉唑的分子式为 C14H14Cl2N2O ，由1 个烯丙氧基连接1个二氯苯基官能团组成。优化后得到的抑霉唑分子结构如图1 所示，主要包括 C、C-H 、C=C-H 、C-O 、N=C-N 、C-N 、N-C-H 和C-N-C等基团。
图2 （a）为理论计算的抑霉唑拉曼光谱，图2 （b）为实际测得的抑霉唑标准品的拉曼光谱。通过对比两者的特征峰归属可发现，1176和1378 cm-1 波段处特征峰的拉曼信号强度较弱。这可能是由于这两个位移处分别在苯环上有微弱的面外弯曲和面内弯曲振动所引起。结合表1和图2 可得，抑霉唑标准品在690 cm-1 波段处的特征峰是由于苯环弯曲振动、CH2面内摇摆、C-C 面外弯曲、二氯苯基面外弯曲共同作用产生的;983 cm-1 波段处的振动可能是由于苯环弯曲振动、二氯苯基上 C=N-C 和 C-N-C 与 C=C-N 的弯曲振动共同作用产生;1041 cm-1 处的峰是由C-H剪式振动、C=C-H 面外弯曲振动和CH2面内弯曲振动共同作用产生。图2中理论计算光谱与标准品的拉曼光谱的特征峰匹配可发现，最小偏差为0，最大偏差为10。这可能是由于实际测量拉曼光谱的标准品为固体粉末，其之间存在相互作用力，而理论计算时的形态为理想的气态单分子，忽略了分子间作用力。基于此，理论计算的抑霉唑拉曼光谱很好地印证了实测的抑霉唑拉曼光谱的正确性，可以将实测的抑霉唑标准品拉曼光谱作为后续分析脐橙表面抑霉唑残留的参考依据。
3.2 脐橙果皮中抑霉唑的SERS光谱分析
由图 3可知，含抑霉唑的脐橙样本在983与1047 cm-1 波段处有区别于空白脐橙样本的拉曼特征峰。而抑霉唑标准品拉曼光谱在983与1041 cm-1 波段处均有特征峰，并在1041 cm-1 波段处有6 cm-1 的偏差，这可能是由于脐橙本身较强的背景峰影响导致的。因此，可以将这两处特征峰作为鉴别脐橙果皮中是否含抑霉唑的依据。
3.3 样本SERS强度检测条件优化
由图4可知脐橙与金胶和NaBr溶液混合物的反应时间不同，SERS信号强度不同。随着反应时间的增长，983和 1047 cm-1 处的 SERS 强度均呈现出持续下降的趋势。产生这种现象的原因可能是刚开始吸附时，脐橙果皮的抑霉唑与金胶及NaBr溶液充分混合，使得该活性位点处的 SERS 信号得以迅速增强[17]。吸附时间过长会使金胶纳米粒子凝聚过度从而产生沉淀，降低抑霉唑与金胶结合的效率[18]。但是由于0 min 的反应时间在实际操作中比较困难，故本研究选择2min作为最优的反应时间。
SERS 主要是通过分子吸附在金和银等粗糙金属的表面，从而使局部表面产生等离激元共振效应，对 SERS 强度起到增强的效果[19]。由图5 可看出，随着金胶加入量的增大，983与 1047 cm-1 波段处的拉曼特征峰强度呈先减小后增大再次减小的趋势。出现这种现象可能的原因是当金胶的量增大到一定程度时，有更多的抑霉唑分子吸附在纳米金胶颗粒的表面，使得983 cm-1 的特征峰在560µL 时的 SERS 信号达到最强，1047 cm-1 的特征峰在400µL时的 SERS信号达到最强。而在金胶加入量为400µL时983 cm-1 波段处的 SERS強度为次强点，仅次于560µL 处的 SERS 强度。综合考虑，本研究选择最优金胶加入量为400µL。
电解质溶液对 SERS强度的增强效果取决于离子对金胶表面状态的影响及离子对金胶纳米粒子的吸附能力[20]。图6 可看出，金胶与 CaBr2的混合溶液在983 cm-1 波段处的 SERS 强度达到最强，金胶与NaBr的混合溶液在1047 cm-1 波段处的 SERS强度最强。出现此现象的原因可能是由于氯离子的盐析能力大于溴离子的盐析能力，导致溴离子与金属粒子的表面络合能力强于氯离子，结果表明含溴离子的化合物（NaBr和 Ca‐Br2）对 SERS强度的增强效果明显大于含氯离子的化合物（NaCl 和 CaCl2）对 SERS 强度的增强效果[21]。其次，NaBr在983 cm-1 特征峰处的SERS强度仅次于CaBr2。综合考虑，本研究选择NaBr作为最优的电解质溶液。
在纳米金胶溶液中加入不同量的NaBr，可能会对 SERS信号产生影响，因此有必要对电解质溶液的加入量进行优化[22]。从图7可看出，当NaBr加入量为 10µL 时， 1047 cm-1 波段处的 SERS信号强度达到局部最高。当NaBr加入量增大到25µL时，983 cm-1 波段处的 SERS强度达到最高点，而1047 cm-1 波段处虽未达到最强，却仅次于最强点。可能的原因是当卤族元素的量高到足以聚合时，使局部的电磁场增强，从而增强 SERS信号。而当加入过量NaBr后，会加剧金纳米粒子的聚集从而发生沉淀[21]，进而降低 SERS 信号强度。因此本试验确定了NaBr溶液的最佳加入量为25µL。
3.4 光谱预处理
为减少光谱分析过程中受到噪声、基线漂移等因素干扰，在建立模型之前有必要对原始光谱进行光谱预处理[23]。目前常用的光谱预处理方法主要有 air PLS 、归一化、基线校正、一阶导数、二阶导数、标准正态变量（Standard Normal Distribution ，SNV）和多元散射校正（Multiplicative Scatter Correction ，MSC）等[23，24]。其中基线校正、一阶导数和二阶导数的作用主要是扣除漂移或背景对结果的影响，在改变了原始数据相关性的变化趋势的同时，也有效提高了某些位置光谱的相关性。在光谱数据分析时，样本的误差表现在无规则分布或样本颗粒大小各异，而MSC 和 SNV可以消除该误差对光谱分析结果的影响，还能够提高原始光谱与预处理后光谱之间的相关性。归一化能够消除由于尺寸差异过大对分析结果带来的误差[23]。因此，本研究以自适应迭代惩罚最小二乘法（Adaptive Iterative Re‐weighted Penalized Least Squares ，air PLS）、airPLS+归一化、air PLS+基线校正、air PLS+一阶导数、air PLS+SNV 和 air PLS+MSC 处理后的6组光谱数据为研究对象，分别采用这6组光谱预处理法建立 SVR模型并进行预测性能的比较。
由表2 可知，air PLS 、air PLS+SNV 、air‐PLS+MSC 三组光谱预处理方法所建立模型的训练集相关系数（Coefficient of Determination forthe Training Set ，RT ）的值均达到0.96以上。且三者的训练集均方根误差（Root-Mean-SquareError of Training set ，RMSEC）的值相差不大，稳定在3.7附近。然后对该三种方法的预测集相关系数（Coefficient of the Determinant for the Pre‐diction Set ，RP ）进行比较，air PLS 方法建立模型预测的预测集RP 为最大，预测集均方根误差（Root-Mean-Square Error of Prediction ，RMSEP）为最小。因此，本研究最终选择 air PLS 作为脐橙表面抑霉唑残留的 SVR预测模型的光谱预处理方法。
3.5 PCA特征提取
建立模型之前对光谱数据进行主成分分析（Principal Component Analysis ， PCA）特征提取。一般来说，主成分分析时主成分的方差贡献率达到90%以上即可作为主要光谱信息的特征[25]。由表3可知，前7个主成分的累计方差贡献率已达91.94%。因此，选择前7个主成分的得分值作为 SVR预测模型的输入值。
3.6 预测模型的建立以及预测性能分析
为评价模型的预测精度和效果，采用了 SVR 、多元线性回归（Multiple Linear Regression ，MLR）和偏最小二乘回歸（Partial Least Squares Regression ，PLSR）三种建模方法进行分析比较。其中，SVR 预测模型的类型为epsilon SVR ，核函数为径向基核函数，epsilon 值为0.1，惩罚因子C值为10。
通常选择 RP 、 RMSEP 和相对标准偏差（Relative Standard Deviation ，RPD）的值来评价模型性能。其中，RP 和RMSEP分别用来解释真实值与预测值的相关性和精确度的能力。RP 越大，RMSEP越小则说明该模型的拟合程度越好，准确性越高[25]。RPD 值则是对模型预测能力的评估，RPD 值越大（不小于1），说明预测效果越好，其预测结果越接近实测值[25]。
由表4可得，基于 SVR建立模型预测的RP 为 SVR 、MLR 、PLSR三者中的最大，RMSEP为最小。因此，基于 SVR算法对脐橙表面抑霉唑残留的预测值越接近实测值，能有效提高模型预测准确性。其中，最低检出限为1.43 mg/kg （国家于《GB 2763-2019食品安全国家标准食品中农药最大残留限量》中规定抑霉唑在脐橙中残留最大限量为5 mg/kg）。
为进一步做预测分析，在随机挑选35个样本建立模型之后，剩余的30个样本作为预测集进行验证，即脐橙果皮抑霉唑残留浓度分别为38.93、30.45、21.53、14.63、9.24和 4.05 mg/kg的样本，实测的真实值与预测值之间的关系如图8，RP 为0.9156，RMSEP为4.8407 mg/kg。
4 結论
本研究探索了SERS技术在检测脐橙果皮中抑霉唑农药残留的可行性。首先通过对比983和1047 cm-1 两处波段的 SERS强度，进行抑霉唑残留的 SERS检测条件优化，分别确定了最优的检测条件（2 min 的反应时间，400µL 的金胶加入量，NaBr作为电解质溶液且添加量为25µL）。在最优条件下，对含抑霉唑的脐橙果皮拉曼光谱采用6组不同的光谱预处理方法，分别对这6组光谱所建立模型的预测性能进行对比，最终采用 air PLS 的光谱预处理方法。经过 PCA 的特征提取，选择前7个主成分得分作为脐橙果皮抑霉唑残留的 SVR预测模型的输入值。最后 SVR 预测模型的RP 可达0.9156，RMSEP 为4.8407 mg/kg， RPD 为2.3103。由于本研究所采用的计算实际残留浓度的方法为称重法，挥发时间不一致容易造成称量误差，进而影响最终实际残留浓度的计算。因此实际残留浓度的计算方法还有待进一步改善。研究结果表明，SERS 技术可实现脐橙果皮抑霉唑残留的快速检测。
参考文献：
[1]唐剑鸿. 新媒体视域下赣南脐橙网络营销策略优化研究[J].食品研究与开发， 2021， 42（8）：229-230. TANG J. Research on the optimization of Gannan navel orange network marketing strategy from the perspective of new media[J]. Food Research and Development， 2021， 42（8）：229-230.
[2]王文军， 曾凯芳， 刘晓佳， 等. 不同保鲜剂对柑橘果实贮藏品质的影响[J].食品与机械 ， 2017， 33（4）：110-116.
WANG W， ZENG K， LIU X， et al. Effects of different preservatives on storage quality of ; citrus fruits[J]. Food and Machinery， 2017， 33（4）：110-116.
[3]陈婷， 蔡艳， 段凯文， 等. 脐橙贮藏保鲜研究进展[J].农产品加工， 2019（20）：83-85.
CHEN T， CAI Y， DUAN K， et al. Research progress of navel orange storage and preservation[J]. Agricultural Products Processing， 2019（20）：83-85.
[4] JIN C， LUO T， FU Z， et al， Chronic exposure of miceto low doses of imazalil induces hepatotoxicity at the physiological， biochemical， and levelstranscriptomic[J]. Environmental Toxicology， 2018， 33（6）：650-658.
[5] VALCKE M， BOURGAULT M H， ROCHETTE L， etal. Human health risk assessment on the consumption of fruits and vegetables containing residual pesticides：
A cancer and non-cancer risk/benefit perspective[J]. Environ International， 2017， 108：63-74.
[6]中华人民共和国农业农村部. 食品安全国家标准食品中农药最大残留限量：
GB 2763—2019[S].北京：中国标准出版社， 2019.
Ministry of Agriculture and Rural Affairs of the Peo‐ple"s Republic of China. National food safety standard—Maximum residue limits of pesticides in food：
GB2763-2019[S]. Beijing：
Standards Press of China， 2019.
[7]谢建军， 陈捷， 李菊， 等. 改良QuEChERS法结合气相色谱串联质谱测定果蔬中20种杀菌剂[J].食品安全质量检测学报， 2013， 4（1）：82-88.
XIE J， CHEN J， LI J， et al. Determination of 20 fungi‐cides in fruits and vegetables by improved QuEChERSmethod combined with gas chromatography tandemmass spectrometry[J]. Journal of Food Safety and Qual‐ity Inspection， 2013， 4（1）：82-88.
[8]潘新明， 潘守奇， 于金玲， 等. 气相色谱法同时检测蔬菜中的7种不同性质农药残留[J].现代仪器， 2012， 18（2）：37-38.
PAN X， PAN S， YU J， et al. Simultaneous determina‐tion of seven pesticide residues in vegetables by gaschromatography[J]. Modern Instrument， 2012， 18（2）：37-38.
[9] LEHOTAY S J， DE KOK A， HIEMSTRA M， et al. Val‐idation of a fast and easy method for the determinationof residues from 229 pesticides in fruits and vegetablesusing gas and liquid chromatography and mass spectro‐metric detection[J]. Journal of AOAC International，2005， 88（2）：595-614.
[10] DANIEL D， LAGO C LDO. Determination of multi‐class pesticides residues in corn by QuEChERS andcapillary electrophoresis tandem mass spectrometry[J].Food Analytical Methods， 2019， 12（7）：1684-1692.
[11] ZHAO H， HASI W， LI N， et al. In situ analysis of pesti‐cide residues on the surface of agricultural products viasurface-enhanced Raman spectroscopy using a flexibleAu@Ag –PDMS substrate[J]. New Journal of Chemis‐try， 2019， 43（33）：13075-13082.
[12] 郭晓湲， 排尔哈提·亚生， 刘晨阳， 等. 拉曼光谱技术的发展及其在生物医学领域中的应用[J].福州大学学报（自然科学版）， 2021， 49（1）：135-142.
GUO X， PAIERHATI Y， LIU C， et al. Development ofRaman spectroscopy and its application in biomedicalfield[J]. Journal of Fuzhou University （Natural ScienceEdition）， 2021， 49（1）：135-142.
[13] 袁雯雯， 沈健. 表面增强拉曼光谱定量分析技术几种方法的探讨[J].化学世界， 2021， 62（4）：193-200.
YUAN W， SHEN J. Discussion on several methods forquantitative analysis of surface enhanced Raman spec‐troscopy[J]. Chemical World， 2021， 62（4）：193-200.
[14] DEVAL V， KUMAR A， GUPTA V， et al. Molecularstructure （monomeric and dimeric） and hydrogen bonds in 5-benzyl 2-thiohydantoin studied by FT-IR and FT-Raman spectroscopy and DFT calculations[J]. Spectrochim Acta A Mol BiomolSpectrosc， 2014， 132：15-26.
[15] 郭红青， 刘木华， 袁海超， 等. 表面增强拉曼光谱技术快速检测鸭肉中的土霉素[J].食品安全质量检测学报， 2017， 8（1）：169-176.
GUO H， LIU M， YUAN H， et al. Rapid detection of oxytetracycline in duck meat by surface enhanced Raman spectroscopy[J]. Journal of Food Safety and Quality Inspection， 2017， 8（1）：169-176.
[16] 黃双根， 王晓， 吴燕， 等. SERS 技术的小白菜中西维因农药残留检测[J].光谱学与光谱分析， 2019， 39（1）：130-136.
HUANG S， WANG X， WU Y， et al. Determination of carbaryl pesticide residues in Chinese cabbage by SERS[J]. Spectroscopy and Spectral Analysis， 2019， 39（1）：130-136.
[17] 王婷， 刘木华， 袁海超， 等. 鸡肉中丙酸睾酮残留表面增强拉曼光谱检测条件的优化[J].食品研究与开发， 2020， 41（2）：135-139， 164.
WANG T， LIU M， YUAN H， et al. Optimization of surface enhanced Raman spectroscopy for the determination of testosterone propionate residue in chicken[J]. Food Research and Development， 2020， 41（2）：
135-139， 164.
[18] ALEXANDER K D， HAMPTON M J， ZHANG S， etal. A high-throughput method for controlled hot-spot fabrication in SERS-active gold nanoparticle dimer arrays[J]. Journal of Raman Spectroscopy， 2009， 40（12）：2171-2175.
[19] 朱越洲， 张月皎， 李剑锋， 等. 表面增强拉曼光谱：应用和发展[J].应用化学， 2018， 35（9）：984-992.
ZHU Y， ZHANG Y， LI J， et al. Surface enhanced Ra‐man spectroscopy：
Application and development[J].Applied Chemistry， 2018， 35（9）：984-992.
[20] XIE L， LU J， CHEN G， et al. Key role of direct adsorp‐tion on SERS Sensitivity：
Synergistic effect among tar‐get， aggregating agent， and surface with Au or Ag col‐loid as SERS substrate[J]. The Journal of PhysicalChemistry Letters， 2020， 11（3）：1022-1029.
[21] GUO H， HUANG Q， LENG W， et al， Bromide ion-functionalized nanoprobes for sensitive and reliable pHmeasurement by surface-enhanced Raman spectrosco‐py[J]. Analyst （London）， 2019， 144（24）：7326-7335.
[22] YASEEN T， SUN D， PU H， et al. Detection of ometho‐ate residues in peach with surface-enhanced Ramanspectroscopy[J]. Food Analytical Methods， 2018， 11（9）：2518-2527.
[23] 第五鵬瑶， 卞希慧， 王姿方， 等. 光谱预处理方法选择研究[J].光谱学与光谱分析， 2019， 39（9）：2800-2806.
DIWU P， BIAN X， WANG Z， et al. Selection of spec‐tral pretreatment methods[J]. Spectroscopy and Spec‐tral Analysis， 2019， 39（9）：2800-2806.
[24] 李淑娟， 卞希慧， 李倩， 等. 基于近红外光谱的四元调和食用油定量分析[J].天津科技大学学报， 2018， 33（3）：18-24.
LI S， BIAN X， LI Q， et al. Quantitative analysis of fourelement mixed edible oil based on NIR [J]. Journal ofTianjin University of Science and Technology， 2018， 33（3）：18-24.
[25] 徐夕博 ， 吕建树 ， 吴泉源 ， 等. 基于 PCA-MLR 和PCA-BPN 的莱州湾南岸滨海平原土壤有机质高光谱预测研究[J].光谱学与光谱分析 ， 2018， 38（8）：2556-2562.
XU X， LYU J， WU Q， et al. Hyperspectral predictionof soil organic matter in coastal plain of the SouthBank of Laizhou Bay based on PCA MLR and PCABPN[J]. Spectroscopy and Spectral Analysis， 2018， 38（8）：2556-2562.
Rapid Detection of Imazalil Residues in Navel Orange Peel Using Surface-Enhanced Raman Spectroscopy
ZHANG Sha1，2， LIU Muhua1，2， CHEN Jinyin2， ZHAO Jinhui1，2*
（1. College of Engineering， Jiangxi Agricultural University/Key Laboratory of Modern Agricultural Equipment in JiangxiProvince， Nanchang 330045， China;2. Jiangxi Provincial Collaborative Innovation Center of Key Technologies and Qualityand Safety in Post-Harvest Processing of Fruits and Vegetables， Nanchang 330045， China）
Abstract：
Imazalil， a preservative for navel orange in the process of postharvest processing， is easy to seep into the flesh through the peel and produce residues in the flesh， which is vulnerable to cause endanger to human body if it was eaten accidentally. Base on this， a fast detection method of imazalil residues in navel orange peel ，namely surface-enhanced Raman spectroscopy （SERS） was proposed in this study. Firstly， the SERS detection conditions of imazalil residues in navel orange peel were optimized， and the optimal detection conditions were determined as follows：
Reaction time of 2 min， gold colloid of 400µL， NaBr as electrolyte solution， NaBr dosage of 25µL. Based on the above optimal conditions， 6 groups of spectral data processed by adaptive iterative penalized least squares （air PLS）， air PLS combination with normalization， air PLS combination with baseline correction， air PLS combination with first derivative， air PLS combination with standard normal distribution （SNV）， air PLS combination with multiplicative scatter correction （MSC） were used to establish support vector regression （SVR） models and compare the models prediction performance. And air PLS method was selected as the spectral pretreatment method， because the value of correlation coefficient computed value of prediction set （RP） is the largest， and the value of root mean square error calculated value of the prediction set （RMSEP） is the smallest. Then， principal component analysis （PCA） was used to extract the features from spectral data， and the first seven principal component scores were selected as the input values of SVR prediction model. SVR， multiple linear regression （MLR） and partial least squares regression （PLSR） were used to analyze and compare the prediction performances. The RP value of prediction set of SVR prediction model could reach 0.9156， the RMSEP value of their prediction set was 4.8407 mg/kg， and the relative standard deviation computation value （RPD） was 2.3103， which indicated that the closer the predicted value of imazalil residue on navel orange surface based on SVR algorithm was to the measured value， the more effective the prediction accuracy of the model could be. The above data indicated that the speedy detection of imazalil residues in navel orange peel could be emploied by SERS coupled with PCA and SVR modeling method.