基于QuEChERS液相色譜-串聯質譜法測定茶葉中除蟲脲殘留量的不確定度評定

武 源，覃慧麗

除蟲脲又名敵滅靈，是一種苯甲酰胺類殺蟲劑，具有對人畜低毒、對環境無公害等特點，其已被廣泛用于茶園害蟲的防治[1-5]。茶葉產業關乎我國農業農村發展，茶葉質量安全影響我國茶葉對外貿易，近年來，農藥殘留已成為影響茶葉質量安全的重要因素之一，除蟲脲及其降解產物殘留問題日益突出[6-8]，嚴重影響到茶葉質量安全。我國國家標準GB 2763—2019《食品安全國家標準食品中農藥最大殘留限量》規定茶葉中除蟲脲的最大殘留限量為20 mg/kg，日本和歐盟等進口國家或地區均對茶葉制定了嚴格的農藥殘留限量[9-10]，日本規定除蟲脲在茶葉中最大殘留限量為20 mg/kg[11]。由于茶葉樣品基質復雜，含有大量的生物堿、色素、咖啡因、茶多酚等干擾物質，使得農殘檢測難度增大。為確保分析結果的準確性，對檢測結果進行不確定度評定是必要的，不確定度評定對于檢測過程的質量把控具有重要意義[12]。

目前除蟲脲在液相色譜-串聯質譜（LC-MS/MS）法的不確定度判定研究基本空缺。本研究依據文獻[3-16]，對測定茶葉中除蟲脲殘留量進行了分析與評定，明確了QuEChERS液相色譜-串聯質譜法測定茶葉中除蟲脲殘留量不確定的主要影響因素，為評價除蟲脲殘留量測定結果的準確性和可靠性提供技術依據。

1 材料與方法

1.1 材料與試劑

提取試劑：4 g MgSO4，1 g NaCl，0.5 g三水合二檸檬酸二鈉，1 g檸檬酸鈉及兩顆陶瓷均質子，美國Agilent公司；凈化試劑管：1 200 mg MgSO4，400 mg PSA，400 mg C18，200 mg GCB，美國Agilent公司。

除蟲脲對照品（純度99.65%），德國Dr. Ehrenstorfe公司；三苯基磷酸酯對照品（TPP，純度99.6%），德國Dr.Ehrenstorfe公司；甲醇、乙腈（均為色譜純），德國Merck公司；甲酸、乙酸銨（均為色譜純），德國Fluka公司。

1.2 儀器與設備

Xevo TQ-S 液 質 聯 用 儀、ACQUITY UPLCBEH C18色譜柱（2.1×50 mm，1.7 μm），美國WATERS公司；XP26微量電子天平，瑞士Mettler Toledo公司；IKA VORTEX 4 basic渦旋振蕩器，德國IKA公司；TTL-DC Ⅱ型氮吹儀，北京同泰聯科技發展有限公司；SIGMA 3-3OK離心機，德國SIGMA公司；Milli-Q超純水設備，美國Millipore公司。

1.3 方法

1.3.1 標準溶液的配制

陰性基質溶液：稱取不含目標物的茶葉按前處理步驟處理即得。

除蟲脲標準儲備液（500 μg/mL）：精確稱取5 mg除蟲脲對照品于10 mL的容量瓶中，用乙腈溶解并定容至刻線。

TPP內標儲備液（500 μg/mL）：精確稱取5 mg內標物對照品TPP于10 mL的容量瓶中，用乙腈溶解并定容至刻線。

除蟲脲標準工作液（500 ng/mL）：精確吸取0.02 mL的除蟲脲標準儲備液于20 mL容量瓶中，用乙腈定容至刻線。

TPP內標工作液（5 μg/mL）：精確吸取0.2 mL的內標儲備液于20 mL容量瓶中，用乙腈溶解定容至刻線。

基質標準工作液：用移液槍分別精確移取0.01、0.02、0.05、0.10、0.20、0.50 mL的除蟲脲標準工作液于10 mL的氮吹管中，加入0.01 mL內標工作液，氮吹近干，準確加入1.0 mL陰性基質溶液，渦旋復溶，即得基質標準工作液。以系列標準溶液中除蟲脲與內標TPP質量濃度的比值（C除蟲脲/CTPP）為橫坐標，以除蟲脲與內標TPP定量離子峰面積的比值（A除蟲脲/ATPP）為縱坐標，繪制基質標準曲線。

1.3.2 提取

茶葉粉碎后進一步充分混勻，準確稱取2 g（精確至0.01 g）茶葉置于50 mL具塞離心管中，加10 mL水渦旋混勻，靜置30 min。離心管中加入0.10 mL TPP內標工作液，10 mL乙腈，加蓋渦旋振蕩30 s后加入提取試劑，立即劇烈振蕩1 min，防止結塊。以轉速4 000 r/min離心5 min，待凈化。

1.3.3 凈化

取8 mL乙腈提取液轉移至15 mL凈化試劑管中，蓋緊離心管，充分渦旋振蕩1 min，以轉速4 000 r/min 離心5 min，凈化液過0.22 μm濾膜后，取1 mL氮吹近干，加入1 mL水-甲醇（9∶1，V/V）渦旋混勻后上機分析。

1.3.4 LC-MS/MS條件1.

3.4.1 液相條件

1.3.4.2 質譜條件

電噴霧離子源（ESI）；多反應監測（MRM）采集模式，正離子掃描。目標物和內標的保留時間、定量離子對和定性離子對見表1。

表1 化合物保留時間及離子對

1.4 數學模型的建立

茶葉中目標物殘留量X按公式（1）計算：

式中：X為茶葉中除蟲脲的殘留量，mg/kg；C0為從基質標準工作曲線得到的茶葉溶液中除蟲脲的質量濃度，ng/mL；V為茶葉經前處理后的稀釋體積，mL；m為稱取茶葉質量，g。

2 結果與分析

2.1 稱量標準物質及配制基質標準曲線時產生的不確定度μr（2.1）

2.1.1 除蟲脲純度產生的不確定度μr(P)

除蟲脲的純度為99.65%，查閱除蟲脲對照品標準證書，其不確定度U95為0.30%，在95%置信水平下，k=2，則標準不確定度μr(P)=0.15%，相對標準不確定度μr(P)=0.001 5。

2.1.2 除蟲脲稱量產生的不確定度μr(M)

精密稱量除蟲脲對照品5.757 mg，用乙腈稀釋并定容至10 mL，得到除蟲脲標準儲備液。微量電子分析天平校準不確定度為主要分量，根據分析天平檢定證書，本試驗所用天平的最大允許誤差為0.005 mg，按照均勻分布計算則除蟲脲標不確定度μr(M)=0.002 9 mg，相對標準不確定度μr(M)=0.000 5。

2.1.3 配制基質標準曲線時產生的不確定度μr(V)

根據JJG196—2006《常用玻璃量器》[17]檢定規程的A級最大允許誤差（以下玻璃量器均按A級進行評定）及廣西壯族自治區計量檢測研究院提供的實驗室移液槍的檢定證書，按照式（2）計算量器容量允差產生的相對標準不確定度μr(V校準)：

式中：a校準為玻璃器具校準引起的體積偏差；k取平均分布，即

實驗室環境溫度波動范圍為（20±5）℃，20℃乙腈的體積膨脹系數為1.37×10-3℃-1，根據式（3）計算室溫波動引起的相對標準不確定度μr(V溫度)：

式中：a溫度為實驗環境溫度波動帶來的溶劑體積偏差；k取平均分布，即

分別移取0.01、0.02、0.05、0.10、0.20、0.50 mL的除蟲脲標準工作液于10 mL的氮吹管中，加入0.01 mL內標工作液，氮吹近干，用移液槍準確加入1 mL陰性基質溶液，充分渦旋。實驗中配制基質標準曲線時移取和定容步驟，用到不同規格的移液槍和容量瓶，均按平均分布計算各分量合成配制基質標準曲線時產生的不確定度μr(V)（表2）。

基質標準曲線稀釋過程由上述19次移液，2次定容完成，配制基質標準曲線時產生的不確定度μr(V)=0.021。

表2 量器引入相對標準不確定度

由此，稱量標準物質及配制基質標準曲線時產生的不確定度μr（2.1）：

2.2 標準曲線擬合過程產生的不確定度μr（2.2）

配制除蟲脲6個質量濃度分別為5.734、11.47、28.67、57.34、114.7、286.7 ng/mL的基質標準曲線進行重復測定3次。由于實驗全程加入同一瓶內標使用液，內標濃度一致，無需考慮內標使用液濃度引入的不確定度，只需考慮內標使用體積引入的不確定度，為方便計算實驗將內標TPP質量濃度設為1 ng/mL。即除蟲脲與內標TPP質量濃度比值x分別為5.734、11.47、28.67、57.34、114.7、286.7，得到除蟲脲與內標TPP定量離子豐度比值y，以x為線性橫坐標，y為線性縱坐標，用最小二乘法進行線性擬合，得到標準線性方程y=kx＋b及相關系數R2，結果詳見表3。

基質標準曲線擬合的標準偏差：

基質標準曲線擬合的相對標準不確定度：

式中：n為基質標準工作液測定的點數，n=6；p為樣品平行測定的個數，p=6；yi為各基質標準工作液濃度對應的響應比值（Ai目標物/Ai內標）；x0為樣品測得濃度均值，以回歸標準曲線校正所得含量均值計，為97.70 ng/mL（表4）；xˉ為標準曲線各點濃度的平均值，為84.10 ng/mL（表3）；b為直線截距，為0.007 3。

表3 基質標曲擬合數據

注：把xi代入直線方程得理論y值，把yi代入直線方程得理論x值。

因此，基質標曲擬合的相對標準不確定度μr(2.2)=0.012。

2.3 茶葉試樣前處理產生的不確定度μr（2.3）

試樣稱量、稀釋時產生的不確定度是未知試樣前處理過程常見的不確定度。實驗全程加入同一瓶內標使用液，內標濃度一致，由于考慮了內標使用體積產生的不確定度，即試樣稀釋體積不影響到定量結果，實驗無需再考慮試樣前處理稀釋體積引起的不確定度。

2.3.1 茶葉試樣稱量產生的不確定度μr(m)

檢定天平檢定證書，其最大允許誤差為0.5 mg，已知稱樣2 g，取矩形分布，則標準不確定度μ(m)=0.289 mg，相對標準不確定度為μr(m)=0.000 1。

2.3.2 茶葉試樣加入內標體積產生的不確定μr(v)

前處理過程在茶葉試樣中加入0.1 mL內標TPP使用液，根據移液槍檢定證書容量允差及溫度波動引起的體積變化（表2），相對標準不確定度μr(v)=0.009。

由此，未知試樣前處理時產生的相對標準不確定度μr(2.3)：

2.4 茶葉試樣測定重復性產生的不確定度μr（2.4）

準確稱取6份茶葉試樣，經提取、凈化步驟，進樣分析。試樣平行測量結果見表4。

表4 試樣測定重復性結果

試樣測定均值Sˉ=0.479 mg/kg（n=6），標準偏差S樣品=0.006 2 mg/kg，標準不確定度mg/kg，相對標準不確定度

2.5 茶葉試樣加標回收產生的不確定度μr（2.5）

稱取陰性基質隨行加標，考察回收率產生的不確定度。陰性基質加標回收實驗進行6次測定（n=6），除蟲脲加標量為114.7 ng，回收率計算結果見表5。

表5 茶葉中除蟲脲加標回收率結果

加收率均值Sˉ=95.2%（n=6），回收率標準偏差S回收=1.4%，標準不確定度相對標準不確定度

使用t檢驗法檢驗平均回收率與回收率期望值100%的差異：

查表得，t0.95（5）=2.57＜8.00，顯示回收率結果與100%有顯著性差異，表明樣品測定需要用回收率進行結果校正，校正后的結果為0.503 mg/kg。

2.6 不確定度的合成

除蟲脲各分量的相對標準不確定度見表6，則除蟲脲的相對合成標準不確定度：

表6 茶葉中除蟲脲各分量的相對標準不確定度統計

折算成除蟲脲合成標準不確定度為：

2.7 擴展不確定度及結果

在P=95%置信水平，測量結果擴展不確定度取擴展因子k=2，則茶葉中除蟲脲殘留量測定的擴展不確定度為：U(C)=kμ(C)=0.028 mg/kg。因此，茶葉試樣中除蟲脲的殘留量表示為（0.503±0.028）mg/kg（k=2）。

3 結論與討論

通過液相色譜-串聯質譜法測定茶葉中除蟲脲殘留量的不確定度，得出該測定方法的擴展不確定度為0.028 mg/kg。從表6可見，稱量標準物質及配制基質標準曲線時產生的不確定度μr(2.1)分量對合成不確定度的貢獻最大，標準曲線擬合過程產生的不確定度μr(2.2)次之。其中，在合成μr(2.1)眾多的小分量中，配制基質標準曲線時引入的不確定度μr(V)，影響最為顯著。因此，配制基質標準曲線時使用允差范圍小的移液槍，保持恒溫恒濕的實驗環境，注意操作人員的熟練水平，把控好內標加入體積的重復性，皆有助于減小配制基質標準曲線引入的不確定度。此外，盡可能提高曲線擬合的線性，例如標準曲線每個點測量多次，測量值取平均值進行擬合，線性范圍內多取幾個濃度點，并確保樣品測定結果在曲線范圍的中高段等，可降低測定結果的不確定度。