LC-MS/MS同時測定茭白中阿維菌素和甲氨基阿維菌素苯甲酸鹽殘留量的不確定度評估

楊德毅，吾建祥*，劉莉，刁銀軍，朱麗燕，虞冰，馬婧妤

(1.金華市農產品質量綜合監督檢測中心，浙江 金華 321007； 2.金華職業技術學院，浙江 金華 321017；3.金華市農業科學研究院，浙江 金華 321007)

阿維菌素(avermectin，ABA)是1976年由日本北里大學學者Takahito Omur等和美國Merck公司首先開發研制的殺蟲劑，是由灰色鏈霉菌發酵產生的一類十六元大環內酯化合物；甲氨基阿維菌素苯甲酸鹽(emamectin benzoate，EMA)是美國Merck公司1984年在阿維菌素B1基礎上合成的一種高效抗生素殺蟲劑。其殺蟲活性是阿維菌素的103倍，殺蟲譜更廣[1]。阿維菌素和甲氨基阿維菌素苯甲酸鹽在茭白生產中常用于防治綠飛虱、螟蟲、蚜蟲等蟲害[2-4]。GB 2763—2016中茭白中的阿維菌素和甲氨基阿維菌素苯甲酸鹽最大殘留限量未作規定，但近年來研究表明，該2種農藥若不合理使用，對生態環境和人類健康存在較大風險[5-6]。本研究參照NY/T 761—2008和GB/T 20769—2008前處理方法和檢測方法，建立了液相色譜-串聯質譜(LC-MS/MS)同時測定茭白中阿維菌素和甲氨基阿維菌素苯甲酸鹽殘留量的分析方法；根據JJF 1059.1—2012《測量不確定度評定與表示》[7]、CNAS—GL06《化學分析中不確定度的評估指南》[8]中規定的程序與方法和相關文獻[9-13]，對該方法進行不確定度評估，分析不確定度的主要來源，尋求影響不確定度的最大因素，并通過相應的措施減少不確定度對結果帶來的影響，為檢測機構質量控制提供科學、準確的依據，為用LC-MS/MS測定其他農藥殘留的不確定度評估提供參考。

1 材料與方法

1.1 材料

測試材料為茭白，系金華市農業科學研究院基地種植的本地品種茭白。

試劑有乙腈(色譜純)、甲醇(色譜純)、乙酸銨(優級純)；甲醇中阿維菌素標準溶液(100 mg·L-1±0.6%)、甲醇中甲氨基阿維菌素甲酸鹽標準溶液(100 mg·L-1±0.6%)，購于農業部環境保護科研監測所；NH2氨基萃取柱，天津博納艾杰爾科技有限公司生產。

所用儀器與設備有液相色譜-串聯四級桿質譜儀(美國Waters公司)，T25均質器(德國IKA公司)，S-EVAP氮氣吹干儀(美國Organomation公司)，固相萃取裝置(美國Supelco公司)，電子天平(瑞士Mettler Toledo公司)，超純水儀(美國Milli-pore公司)，渦旋振蕩器(美國Thermo公司)，移液器(德國艾本德公司)及實驗室常規儀器。

1.2 方法

1.2.1 樣品前處理

將茭白樣品切碎，充分混勻，采用四分法縮分后，經食品加工器搗碎，裝入封口容器中，貼上標簽，-20 ℃保存，待測。

準確稱取25.0 g樣品(精確至0.01 g)于250 mL燒杯中，加入50 mL乙腈和50 mL超純水，勻質器(轉速1 300 r·min-1)充分勻漿2 min后經濾紙過濾至裝有10.0 g氯化鈉的具塞量筒中，濾液劇烈振蕩1 min后，靜置20 min，用移液器轉移10 mL上清液至小燒杯中，氮吹至近干，待凈化[14]。

上述濃縮樣品用2.0 mL二氯甲烷+甲醇(95 mL+5 mL)溶解后過NH2萃取柱。該小柱先用4.0 mL二氯甲烷+甲醇(95 mL+5 mL)預淋條件化，當溶劑液面到達柱吸附層表面時，立即倒入樣品溶液，用15 mL尖底刻度離心管接收洗脫液，用4 mL二氯甲烷+甲醇(95 mL+5 mL)潤洗燒杯后淋洗NH2萃取柱，并重復1次。將盛有淋洗液的刻度離心管置于水浴(50 ℃)氮吹儀上氮吹蒸發至近干，用甲醇在刻度離心管內定容至5.0 mL，旋渦混勻，移入2 mL自動進樣器樣品瓶中，待測[15]。

1.2.2 混合標準液的配制

分別用1 mL的單標線吸量管各吸取1 mL阿維菌素標準溶液和甲氨基阿維菌素甲酸鹽標準溶液于10 mL容量瓶中，用甲醇定容，配制成10 mg·L-1的混合標準溶液。

用1 mL的單標線吸量管和10 mL容量瓶逐級稀釋上述混合標準溶液分別得到1 mg·L-1、100 μg·L-1、10 μg·L-1混合標準溶液；吸取5 mL 1 mg·L-1混合標準溶液定容至10 mL容量瓶中，得到500 μg·L-1混合標準溶液；吸取1 mL 500 μg·L-1混合標準溶液定容至25 mL容量瓶中，得到20 μg·L-1混合標準溶液；用1 mL的單標線吸量管和10 mL容量瓶逐級稀釋500 μg·L-1混合標準溶液分別得到50 μg·L-1、5 μg·L-1混合標準溶液。

按上述方法，得到質量濃度分別為5、10、20、50、100 μg·L-1的混合標準工作液。

1.2.3 儀器條件

色譜柱：Phenomenex Luna(150 mm×2.0 mm，5 μm)；柱溫：35 ℃；樣品室溫度：10 ℃；進樣體積：5 μL；流動相A：甲醇，流動相B：0.1%乙酸銨水溶液；流速：0.2 mL·min-1，梯度洗脫條件見表1。

表1 梯度洗脫條件

質譜掃描方式：多反應監測(MRM)；離子模式(Ionization Mode)：電噴霧正離子源ESI+；毛細管電壓(Capillary)：3.00 kV；離子源溫度(Source Temperature)：150 ℃；錐孔反吹氣流量(Cone Gas Flow)：40 L·h-1；脫溶劑氣溫度(Desolvation Temperature)：500 ℃；脫溶劑氣流量(Desolvation Gas Flow)：800 L·h-1；母離子m/z：886.2；子離子m/z：158.2(定量離子)、126.1；錐孔電壓：80、80 V；碰撞能量：50、62 eV。

1.2.4 數學模型的建立

根據測量的原理建立測量值的數學模型，進行不確定度來源的分析。采用標準曲線法定量，數學模型：

X=cV1V3÷V2÷m。

式中：X為樣品中被測組分的含量(μg·kg-1)；c為從標準工作曲線得到的試樣溶液中被測組分的濃度(μg·L-1)；V1為提取溶劑總體積(mL)；V2為吸取用于檢測的提取溶液體積(mL)；V3為試樣溶液最終的定容體積(mL)；m為試樣的質量(g)。

1.2.5 不確定度分量的主要來源

根據該方法測量過程，測定茭白樣品中阿維菌素和甲氨基阿維菌素甲酸鹽殘留量的不確定度的來源主要包括標準品的引入、標準溶液的配制、標準曲線的擬合、樣品稱量、試樣轉移和定容、測量重復性、回收率引入等。

2 結果與分析

2.1 單一因素不確定度

2.1.1 標準品引入的不確定度

根據標準物質證書，阿維菌素和甲氨基阿維菌素甲酸鹽的純度誤差均為±0.6%；按B類評定取矩形分布，由標準品引入的相對不確定度(urel(sp))：

2.1.2 標準溶液配制過程中引入的不確定度

玻璃量器引入的不確定度。根據JJG 196—2006《常用玻璃量器》[16]檢定規程規定的最大允許誤差，按以下公式計算標準溶液配制過程中玻璃量具允差引起的相對不確定度(表2)，配制過程中所用的玻璃量具均為A級，標準溶液配制過程中引入的不確定度(urel(sv))：

表2 標準溶液配制過程中玻璃量器引入的 不確定度

根據表1中的數據合成標準溶液稀釋過程中玻璃量器與移液器引入的相對不確定度為：

urel(sv)=

0.011 1。

由環境溫度變化引入的不確定度(urel(st))。標準溶液在稀釋過程中的溶劑為甲醇，實驗溫度為(20±5)℃范圍內時(溫度波動為5 ℃)，甲醇的膨脹系數為0.001 19 ℃-1，玻璃量器的膨脹系數遠小于液體的膨脹系數，忽略不計，根據以下公式計算環境溫度變化引起的相對不確定度，計算結果見表3。

根據表3中的數據合成標準溶液在配制過程中由溫度變化引入的相對不確定度：

表3 標準溶液配制過程中溫度變化引入的 不確定度

urel(st)=

0.013 7。

合成標準溶液配制過程中引入的不確定度(urel(svt))。根據標準溶液配制過程，2 種農藥的標準溶液在配制過程中引入的相對不確定度相同，合成 2 種農藥標準溶液配制過程由玻璃量具和溫度引入的相對不確定度為：

2.1.3 標準曲線擬合引入的相對不確定度

通過對質量濃度分別為5、10、20、50、100 μg·L-1的阿維菌素和甲氨基阿維菌素甲酸鹽混合標準溶液測定1次，采用最小二乘法擬合標準溶液的質量濃度與峰面積曲線，各數據結果詳見表4。

同時按照上述方法對在田間施用過2種農藥的陽性茭白樣品進行3次平行測定，測得兩種樣品溶液平均值c0(ABA)=25.15 μg·L-1、c0(EMA)=108.6 μg·L-1(根據1.2.4節公式計算陽性樣品平均殘留量X0，結果見表5)，根據以下公式對標準曲線擬合引入的標準不確定度進行計算，結果見表4。

表4 校準曲線擬合數據及不確定度

表5 陽性樣品測定結果

2.1.4 試樣質量稱量引入的不確定度

稱量時稱取25.0 g(精確到0.1 g)的試樣，根據鑒定證書天平最大允許誤差為±0.05 g，按B類評定取矩形分布，由試樣稱量引入的相對不確定度(urel(dm))：

2.1.5 試樣轉移和定容引入的不確定度

在樣品前處理過程中溫度變化范圍為(20±5)℃，用移液器轉移10 mL乙腈溶液(膨脹系數為0.001 37 ℃-1)，最終用甲醇在15 mL尖底刻度離心管定容至5 mL，根據JJG 646—2006《移液器檢定規程》[17]，10 mL移液器的最大允許誤差為±0.6%，根據JJG 10—2005《專用玻璃量器檢定規程》[18]10 mL尖底刻度離心管最大允差為±0.2 mL，尖底刻度離心管、移液器及溫度變化取矩形分布，按上述2.1.2方法分別計算移液器、玻璃量具引入的相對不確定度urel(dv)、溫度引入的相對不確定度urel(dt)及合成試樣轉移和定容引入相對不確定度(urel(dvt))，結果見表6。

由表6數據合成試樣轉移和定容過程引入的相對不確定度：

表6 前處理過程中移液器和玻璃量具的 相對不確定度

2.1.6 由測量重復性引入的不確定度

在田間施藥后取回的陽性茭白樣品經3次重復測定的結果見表5，按照農藥殘留取樣及測定的要求，對樣品充分混勻、加大樣品的制樣量及科學縮分方式保障試樣的均勻性和檢測結果的代表性。同時按A類評定對其不確定度(urel(df))進行計算，計算公式：

計算得出由測量重復性引入的相對標準不確定度為：urel(df(ABA))=0.016 97；urel(df(EMA))=0.003 446。

2.1.7 回收率引入的不確定度

在空白試樣中添加 5 組不同的混合標準溶液進行測定，按 A 類進行評定，由回收率引入的標準不確定度(urel(dr))根據如下公式計算，該方法的回收率及標準不確定度結果見表5。

按照以下計算公式對回收率進行顯著性檢驗，對照t檢驗臨界值分布表(f=n-1=4)得雙側臨界值t(0.05,4)為2.776。按照以上t值公式計算檢驗值t(表7)均大于2.776，因此與回收率100%有顯著差異，在計算結果時必須采用回收率進行修正。修正后陽性茭白樣品中阿維菌素和甲氨基阿維菌素甲酸鹽殘留量為：

表7 回收率和不確定度

因此由回收率引起的相對不確定度為：

2.1.8 進樣體積引入的不確定度

LC-MS/MS 配備的微量注射器體積為 10 μL，相對標準偏差(RSD)為±1%，樣品上機進樣量為 5 μL，按均勻分布，由體積引入的相對不確定度(urel(ev))：

2.2 復合因素不確定度

2.2.1 合成標準不確定度

LC-MS/MS 法同時測定茭白中阿維菌素和甲氨基阿維菌素苯甲酸鹽殘留的相對不確定度分量見表8，相對合成不確定度(urel(c))計算公式：

表8 相對不確定度的分量及相對合成不確定度值

注：SP，標準品純度；SVT，配制標準溶液；SW，標準曲線擬合；DM，樣品稱量； DVT，試樣轉移和定容；DF，測量重復性；DR，回收率；EV，進樣體積。

分別計算阿維菌素和甲氨基阿維菌素甲酸鹽的相對合成不確定度為：

urel(c(ABA))=0.0516 00；urel(c(EMA))=0.038 700。

根據合成標準不確定度的計算公式uc(X)=urel(X)×X0，其中X0為表5中實測陽性茭白樣品中殘留量的平均值，分別計算得到阿維菌素和甲氨基阿維菌素甲酸鹽的合成標準不確定度為：uc(X(ABA))=1.298；uc(X(EMA))=4.203。

2.2.2 擴展不確定度

依據JJF 1135—2005《化學分析測量不確定度評定》[19]，在95%的置信水平下，取包含因子k=2，擴展不確定度u=k×uc(X)，陽性樣品中阿維菌素和甲氨基阿維菌素甲酸鹽殘留量及不確定度評定結果見表9。

表9 不確定度的評定結果

2.3 影響不確定度因素

從各相對不確定度分量值(表8)看出，用LC-MS/MS法同時測定茭白中阿維菌素和甲氨基阿維菌素苯甲酸鹽殘留量產生的相對不確定度范圍分別在0.001 155～0.034 620和0.001 155～0.027 300，影響相對不確定度主要因素為配制系列標準溶液、標準曲線擬合、試樣轉移和定容等(圖1)。

圖1 測定茭白中阿維菌素和甲氨基阿維菌素 苯甲酸鹽不確定度來源

3 小結與討論

通過對茭白中阿維菌素和甲氨基阿維菌素苯甲酸鹽殘留量測定過程中引入的各相對不確定度評定分析，結果表明，測定茭白中阿維菌素殘留量的過程中對不確定度貢獻較大的因素依次為標準曲線擬合、試樣轉移和定容、標準溶液配制；測定茭白中甲氨基阿維菌素苯甲酸鹽的過程中對不確定度貢獻較大的因素依次為試樣轉移和定容、標準溶液配制、標準曲線擬合。其中在試樣轉移和定容過程中相對不確定度較大原因主要為定容過程中刻度離心管最大允許誤差較大引入。因此，在實際檢測過程中，可通過改進前處理等方法、選擇高精密度的量器、提高檢測人員業務水平和操作的規范性、選擇高純度的標準品等方法減小測量數據的不確定度，提高檢測結果的準確性。