DOI：10.3969/j.issn.10001565.2024.02.005

摘" 要：為建立菊花中29種農(nóng)藥殘留的氣相色譜-串聯(lián)質(zhì)譜測(cè)定方法及不確定度評(píng)價(jià)，采用氣相色譜-串聯(lián)質(zhì)譜內(nèi)標(biāo)法測(cè)定菊花中29種農(nóng)藥殘留的含量，并依據(jù)測(cè)量不確定度的相關(guān)評(píng)定標(biāo)準(zhǔn)與表示，通過(guò)建立數(shù)學(xué)模型，對(duì)檢測(cè)過(guò)程中產(chǎn)生的不確定度進(jìn)行分析與評(píng)價(jià)，包括標(biāo)準(zhǔn)溶液配制、標(biāo)準(zhǔn)曲線擬合、供試品測(cè)定3大方面.結(jié)果表明：菊花中29 種農(nóng)藥質(zhì)量濃度在0.002～0.1 mg/L內(nèi)線性關(guān)系良好，回歸方程決定系數(shù)R2為0.997 36～0.999 82，定量限0.002 03～1.18 963 mg/L，檢出限0.000 64～0.356 89 mg/L，添加回收實(shí)驗(yàn)中測(cè)定平均回收率為69.85%～93.26%，29種農(nóng)藥測(cè)量數(shù)值為0.04、0.06、0.1 mg/kg時(shí)，得到擴(kuò)展不確定度為0.004 787～0.058 935 mg/kg.測(cè)量結(jié)果不確定度顯示標(biāo)準(zhǔn)工作液的配制和標(biāo)準(zhǔn)曲線擬合是本實(shí)驗(yàn)測(cè)定方法不確定度的主要影響因素，該實(shí)驗(yàn)結(jié)果可以為測(cè)定菊花中農(nóng)藥殘留量的方法及不確定度評(píng)價(jià)提供參考.

關(guān)鍵詞：氣相色譜-串聯(lián)質(zhì)譜法；菊花；農(nóng)藥殘留；不確定度

中圖分類號(hào)：TS207.3""" 文獻(xiàn)標(biāo)志碼：A""" 文章編號(hào)：10001565（2024）02014610

Uncertainty evaluation of 29 pesticide residues in chrysanthemum by gas chromatography-tandem mass spectrometry

ZHAO Zhilei1，2，3， SUN Ming1，2， LIU Dongdong1，2， LI Xun1，2， MA Lihui1，2

（1.Hebei Key Laboratory of Energy Metering and Safety Testing Technology， College of Quality and Technical Supervision， Hebei University， Baoding" 071002， China； 2. National amp; Local Joint Engineering Research Center of Metrology Instrument and System， Hebei University， Baoding" 071002， China； 3.Institute of Geographical Indications，Hebei University， Baoding" 071002， China ）

Abstract： In order to establish the determination method of 29 pesticide residues in chrysanthemum by gas chromatography-tandem mass spectrometry and its uncertainty evaluation， the content of 29 pesticide residues in chrysanthemum was determined by gas chromatography-tandem mass spectrometry with internal standard method， and the uncertainty generated in the detection process was analyzed and evaluated by establishing a mathematical model according to the evaluation criteria and expression of the measurement uncertainty. It includes three aspects： preparation of standard solution， fitting of standard curve and determination of test products. The results showed that the linearity of the 29 pesticides in chrysanthemum

收稿日期：20221225;修回日期：20230606

基金項(xiàng)目：河北省重點(diǎn)研發(fā)計(jì)劃項(xiàng)目（21344801D）

第一作者：趙志磊（1977—），男，河北大學(xué)教授，博士，主要從事食品藥品質(zhì)量檢測(cè)技術(shù)研究，E-mail：15930721279@163.com

通信作者：馬力輝（1967—），男，河北大學(xué)教授，博士，主要從事檢測(cè)技術(shù)及產(chǎn)品創(chuàng)新研究.E-mail：malihui@aliyun.com

were good in the range of 0.002～0.1 mg/L， and the coefficient of determination R2 was between 0.997 36～0.999 82. The limits of quantitation were 0.002 03～1.189 63 mg/L， and the limits" of detection were 0.000 64～0.356 89 mg/L. In addition， the average recoveries were between 69.85% and 93.26%， and the extended uncertainties were between 0.004 787～0.058 935 mg/kg for the measured values of 0.04， 0.06 and 0.1 mg/kg for the 29 pesticides. The uncertainty of measurement results showed that the preparation of standard working solutions and the fitting of standard curve were the main factors affecting the uncertainty of the measurement method in this experiment. The results could provide reference for determination of pesticide residues in chrysanthemum and uncertainty evaluation.

Key words： gas chromatography-tandem mass spectrometry; the chrysanthemum; pesticide residues; the uncertainty

菊花是菊科、菊屬的多年生宿根草本植物，具有散風(fēng)清熱、消炎抗菌、保肝降壓、抗氧護(hù)眼等功效，屬藥食同源植物［1-2］.市場(chǎng)上的菊花多為人工栽培，菊花在栽培過(guò)程中病蟲(chóng)害較多，其防治主要依賴使用農(nóng)藥，不規(guī)范的農(nóng)藥使用導(dǎo)致菊花中農(nóng)藥殘留超標(biāo)現(xiàn)象時(shí)有發(fā)生［3-4］，農(nóng)殘超標(biāo)嚴(yán)重影響菊花的質(zhì)量安全.因此，監(jiān)測(cè)菊花中的農(nóng)藥殘留勢(shì)在必行，其檢測(cè)的準(zhǔn)確性對(duì)于菊花質(zhì)量監(jiān)管及保護(hù)消費(fèi)者的健康具有重要意義.不確定度的評(píng)定可以保證檢測(cè)結(jié)果更準(zhǔn)確、可靠，為檢測(cè)結(jié)果的質(zhì)量控制提供依據(jù)［5］.不同的農(nóng)藥殘留測(cè)定方法引入的不確定度也不盡相同，為了確保檢驗(yàn)結(jié)果的可信度，必須對(duì)所用測(cè)定方法及實(shí)驗(yàn)過(guò)程進(jìn)行測(cè)量不確定度的分析與評(píng)定.

氣相色譜-串聯(lián)質(zhì)譜法（GC-MS/MS）是目前較為常用的農(nóng)藥殘留測(cè)定方法之一，其能夠在存在復(fù)雜基質(zhì)的情況下對(duì)化合物進(jìn)行準(zhǔn)確的識(shí)別，對(duì)于測(cè)定食品及藥材中農(nóng)藥殘留及相應(yīng)的代謝產(chǎn)物十分有效［6-7］.韓梅等［8］ 采用GC-MS/MS法建立了檢測(cè)菊花中多種農(nóng)藥殘留的方法，為菊花中多農(nóng)藥殘留的篩査與定量分析提供參考.此外，李俊等［9］研究采用GC-MS對(duì)茶葉中噠螨靈等9種農(nóng)藥含量進(jìn)行測(cè)定，并對(duì)實(shí)驗(yàn)過(guò)程中產(chǎn)生的不確定度作出分析，得出了影響不確定度的主要因素.殷超等［10］ 對(duì)GC-MS法測(cè)定葡萄樣品中苯醚甲環(huán)唑含量進(jìn)行了不確定度評(píng)估，對(duì)客觀評(píng)價(jià)實(shí)驗(yàn)檢測(cè)結(jié)果的準(zhǔn)確性提供了參考.陽(yáng)曦等［11］ 詳細(xì)評(píng)定了QuEChERS-GC-MS法測(cè)定結(jié)球甘藍(lán)中甲基毒死蜱農(nóng)藥殘留的不確定度，為改進(jìn)和提高檢測(cè)質(zhì)量提供了科學(xué)可靠的依據(jù).目前，對(duì)于食品藥品中農(nóng)藥檢測(cè)不確定度的研究大多集中在1種或幾種農(nóng)藥的不確定度分析，對(duì)菊花中多種農(nóng)藥殘留檢測(cè)的不確定度分析鮮有報(bào)道.

本研究從2020年版《中國(guó)藥典》藥材及飲片（植物類）規(guī)定的禁用農(nóng)藥中選擇了29種農(nóng)藥作為研究對(duì)象，依據(jù)相關(guān)檢驗(yàn)規(guī)定［12，13］，建立氣相色譜-串聯(lián)質(zhì)譜法測(cè)定菊花中29 種農(nóng)藥殘留的方法，并對(duì)測(cè)定結(jié)果進(jìn)行不確定度評(píng)定及分析，為客觀評(píng)價(jià)和正確運(yùn)用菊花農(nóng)藥殘留檢測(cè)方法及結(jié)果提供參考.

1" 材料與方法

1.1" 材料與試劑

菊花購(gòu)于河北安國(guó)藥材市場(chǎng)；農(nóng)藥混合標(biāo)準(zhǔn)品（40～100 mg/L，上海安譜實(shí)驗(yàn)科技股份有限公司）；HLB固相萃取柱（涿州凱斯科生物技術(shù)有限公司）；磷酸三苯酯標(biāo)準(zhǔn)品（1 000 mg/L，上海安譜實(shí)驗(yàn)科技股份有限公司）；氯化鈉（分析純，天津市匯杭化工科技有限公司）；乙腈（色譜純，天津市科密歐化學(xué)試劑有限公司）.

1.2" 實(shí)驗(yàn)儀器

GCMS-TQ8040NX氣相色譜-質(zhì)譜聯(lián)用儀（日本島津公司）；PX225DZH電子天平（美國(guó)OHAUS公司）；FC5718R冷凍離心機(jī)（美國(guó)OHAUS公司）；濾膜（0.22 μm，津隆公司）；旋轉(zhuǎn)蒸發(fā)儀（上海葉招科技有限公司）；E0AA-HM-01多管漩渦混合儀（上海安譜實(shí)驗(yàn)科技股份有限公司）；WH-861旋渦混合器（太倉(cāng)市華利達(dá)實(shí)驗(yàn)設(shè)備有限公司）；WJX-350高速多功能粉碎機(jī)（上海緣沃工業(yè)有限公司）.

1.3" 檢測(cè)條件

1.3.1" 氣相色譜條件

氣相色譜升溫程序：初始溫度60 ℃，保持1 min.30 ℃/min升溫至120 ℃；再10 ℃/min升溫至160 ℃；以2 ℃/min升溫至230 ℃；再15 ℃/min升溫至300 ℃并保持6 min.載氣為氦氣，流速1.2 mL/min；進(jìn)樣量1.0 μL；不分流進(jìn)樣.

1.3.2" 質(zhì)譜條件

電子轟擊 （EI） 離子源；電子能量70 eV；離子源溫度250 ℃；接口溫度250 ℃；MRM采集模式；溶劑延遲時(shí)間3 min；SHIMADZU GCMS工作軟件.

1.4" 農(nóng)藥標(biāo)準(zhǔn)品溶液配制

混合農(nóng)藥標(biāo)準(zhǔn)品儲(chǔ)備液配制：吸取0.5 mL農(nóng)藥混合標(biāo)準(zhǔn)品40～100 mg/L于100 mL容量瓶中，用乙腈稀釋定容，配成質(zhì)量濃度為0.2～0.5 mg/L 的混合農(nóng)藥標(biāo)準(zhǔn)品儲(chǔ)備液，-18 ℃避光保存.

內(nèi)標(biāo)物溶液配制：準(zhǔn)確吸取25 μL磷酸三苯酯標(biāo)品（1 000 mg/L）于25 mL容量瓶中，用乙腈稀釋配制成1 mg/L內(nèi)標(biāo)溶液，-18 ℃避光保存.

農(nóng)藥標(biāo)準(zhǔn)品溶液配制：精密吸取混合農(nóng)藥標(biāo)準(zhǔn)品儲(chǔ)備液10、20、50、100、150、200 μL分別置于棕色小瓶中，用前處理方法配制的空白菊花基質(zhì)溶液逐級(jí)稀釋至1 mL，再分別加入內(nèi)標(biāo)物溶液10 μL.渦旋混勻，配成系列質(zhì)量濃度的基質(zhì)匹配農(nóng)藥標(biāo)準(zhǔn)品溶液，現(xiàn)配現(xiàn)用.

加標(biāo)農(nóng)藥溶液配制：將0.5 mL的農(nóng)藥混合標(biāo)準(zhǔn)品（40～100 mg/L）置于10 mL容量瓶中，加乙腈稀釋配制成質(zhì)量濃度為2～5 mg/L 的溶液，搖勻，-18 ℃避光保存.

1.5" 菊花供試品前處理

提?。壕芊Q取菊花供試品粉末（過(guò)三號(hào)篩）5.0 g，加氯化鈉1.0 g，搖散，加入50 mL乙腈2 500 r/min勻漿5 min，4 000 r/min離心5 min，分取上清液，沉淀再加乙腈50 mL，勻漿3 min，離心，合并2次提取的上清液，減壓濃縮至約4 mL，放冷，用乙腈稀釋至10 mL，搖勻備用.

凈化：由于菊花供試品中有機(jī)酸類、黃酮類、揮發(fā)油等小分子成分含量較高［14］，易干擾測(cè)定，且對(duì)儀器污染較大，因此在進(jìn)樣前，吸取5 mL菊花乙腈提取液過(guò)HLB固相萃取小柱，流速每秒1～2滴，收集凈化流出液，過(guò) 0.22 μm 有機(jī)濾膜后上機(jī)測(cè)定.

1.6" 數(shù)學(xué)模型的建立

菊花中29 種農(nóng)藥殘留量的結(jié)果計(jì)算

W=c×vm×f×1 0001 000，

其中：W表示菊花中目標(biāo)物農(nóng)藥的殘留量，mg/kg；c表示供試品溶液中目標(biāo)物農(nóng)藥的質(zhì)量濃度，mg/L；v表示供試品溶液的定容體積，mL；m表示菊花供試品的質(zhì)量，g；f：稀釋因子；1 000 為單位轉(zhuǎn)換系數(shù).

根據(jù)上述數(shù)學(xué)模型及測(cè)定方法，菊花中29 種農(nóng)藥殘留測(cè)定方法的不確定度主要來(lái)源于農(nóng)藥標(biāo)準(zhǔn)品溶液配制（混合農(nóng)藥標(biāo)準(zhǔn)品儲(chǔ)備液、農(nóng)藥標(biāo)準(zhǔn)品溶液）；標(biāo)準(zhǔn)曲線擬合；供試品測(cè)定（包括前處理過(guò)程中的稱量、提取、進(jìn)樣、加標(biāo)回收率及測(cè)定重復(fù)性）.

2" 結(jié)果與分析

2.1" 農(nóng)藥標(biāo)準(zhǔn)品溶液配制引入的相對(duì)不確定度

2.1.1" 混合農(nóng)藥標(biāo)準(zhǔn)品儲(chǔ)備液配制引入的相對(duì)標(biāo)準(zhǔn)不確定度urel（sr）

混合農(nóng)藥標(biāo)準(zhǔn)品儲(chǔ)備液配制的操作步驟即稀釋定容農(nóng)藥混合標(biāo)準(zhǔn)品.因此，在混合農(nóng)藥標(biāo)準(zhǔn)品儲(chǔ)備液配制的過(guò)程中，其不確定度主要受玻璃量器選用和溫度的變化因素影響.

主要選用的玻璃量器有100～1 000 μL可調(diào)移液器，用于量取500 μL農(nóng)藥混合對(duì)照品儲(chǔ)備液1次及10 mL容量瓶（A級(jí)）定容1次.參考JJG 196—2006《常用玻璃量器》［15］及文獻(xiàn)［16］，當(dāng)移液器的標(biāo)稱容量在容量誤差表中給定的2個(gè)容量點(diǎn)區(qū)間內(nèi)時(shí)，以2個(gè)容量點(diǎn)之間較大容量為標(biāo)準(zhǔn).因此，可以得到使用100～1 000 μL可調(diào)移液器的容量允差ΔV為5 μL，使用10 mL容量瓶的容量允差ΔV為0.02 mL.采用B類評(píng)定，按矩形分布計(jì)算，取k=3，量取V（mL）溶液，次數(shù)n=1，公式如下：

urel（Vsr）=n×ΔVk×V.（1）

實(shí)驗(yàn)室環(huán)境溫度通常保持（20±5） ℃（Δt=5 ℃），溫度變化會(huì)造成試劑膨脹，引入不確定度.乙腈的體積膨脹系數(shù)α=1.37×10-3/℃，按矩形分布計(jì)算，取k=3，量取V（mL）溶液，次數(shù)n = 1，公式如下：

urel（tsr）=n×α×Δt×Vk×V.（2）

合成urel（Vsr）和urel（tsr），則由混合農(nóng)藥標(biāo)準(zhǔn)品儲(chǔ)備液配制引入的相對(duì)標(biāo)準(zhǔn)不確定度計(jì)算如下：

urel（sr）=u2rel（Vsr）+u2rel（tsr）=0.007 093.

2.1.2" 農(nóng)藥標(biāo)準(zhǔn)品溶液配制引入的相對(duì)標(biāo)準(zhǔn)不確定度urel（sw）

根據(jù)1.4節(jié)的配制方法，農(nóng)藥標(biāo)準(zhǔn)品溶液主要由菊花基質(zhì)溶液、混合標(biāo)準(zhǔn)儲(chǔ)備液及內(nèi)標(biāo)溶液制取而成.加標(biāo)用農(nóng)藥溶液同混合農(nóng)藥標(biāo)準(zhǔn)品儲(chǔ)備液制備過(guò)程相同.因此，urel（sw） 來(lái)自玻璃量器引入的相對(duì)標(biāo)準(zhǔn)不確定度urel（Vsw）和溫度變化引入的相對(duì)標(biāo)準(zhǔn)不確定度urel（tsw）2項(xiàng)合成.

實(shí)驗(yàn)室環(huán)境溫度誤差Δt = ±5 ℃，乙腈的體積膨脹系數(shù)α= 1.37×10-3/℃，按矩形分布（B類評(píng)定），取k=3，則農(nóng)藥標(biāo)準(zhǔn)品溶液配制所用玻璃量器引入的相對(duì)標(biāo)準(zhǔn)不確定度urel（Vsw）和溫度變化引入的相對(duì)標(biāo)準(zhǔn)不確定度urel（tsw），依據(jù)式（1）、（2）計(jì)算，結(jié)果見(jiàn)表1.由農(nóng)藥標(biāo)準(zhǔn)品溶液引入的相對(duì)標(biāo)準(zhǔn)不確定度為

urel（sw）=u2rel（Vsw）+u2rel（tsw）=0.058 916.

2.2" 標(biāo)準(zhǔn)曲線擬合引入的相對(duì)標(biāo)準(zhǔn)不確定度urel（sc）

取1.4節(jié)配制的6種質(zhì)量濃度的農(nóng)藥標(biāo)準(zhǔn)品溶液檢測(cè)，內(nèi)標(biāo)法定量，以各種農(nóng)藥與內(nèi)標(biāo)的質(zhì)量濃度比值為橫坐標(biāo)，與內(nèi)標(biāo)物的峰面積比值為縱坐標(biāo)，繪制標(biāo)準(zhǔn)曲線.每個(gè)質(zhì)量濃度點(diǎn)重復(fù)測(cè)定3次，使用最小二乘法進(jìn)行線性擬合，得到29 種農(nóng)藥殘留的線性方程Y=aX+b.繪制標(biāo)準(zhǔn)曲線，決定系數(shù)R2為0.997 36～0.999 82，29 種農(nóng)藥在0.002～0.1 mg/L內(nèi)線性關(guān)系良好.在添加水平0.001 mg/L，3倍的信噪比為檢出限，10倍的信噪比為定量限［17］.根據(jù)測(cè)得的數(shù)據(jù)，按公式（3）、（4）計(jì)算由標(biāo)準(zhǔn)曲線擬合引入的相對(duì)標(biāo)準(zhǔn)不確定度.結(jié)果見(jiàn)表2.

usc=sa×1p+1n（X0-X）2∑ni=1（Xi-X）2，（3）

urel（sc）=uscX.（4）

式（3）中，usc為標(biāo)準(zhǔn)曲線擬合引入的標(biāo)準(zhǔn)不確定度；s 為線性回歸方程剩余標(biāo)準(zhǔn)差；a為擬合曲線的斜率；p為被測(cè)供試品的平行測(cè)量次數(shù)；n為標(biāo)準(zhǔn)曲線擬合的數(shù)據(jù)總數(shù)（標(biāo)準(zhǔn)溶液重復(fù)測(cè)定次數(shù)×標(biāo)準(zhǔn)溶液濃度點(diǎn)個(gè)數(shù)）；X0為加標(biāo)供試品通過(guò)線性回歸方程校正所得相對(duì)濃度比平均值；X為標(biāo)準(zhǔn)溶液各相對(duì)質(zhì)量濃度點(diǎn)平均值（農(nóng)藥標(biāo)準(zhǔn)溶液的平均相對(duì)質(zhì)量濃度）；Xi 為線性曲線回歸方程校正所得各點(diǎn)相對(duì)質(zhì)量濃度比（目標(biāo)物與內(nèi)標(biāo)的質(zhì)量濃度比）.

s=∑ni=1［Yi-（aXi+b）］2n-2，（5）

式（5）中：b為擬合曲線的截距；a為擬合曲線的斜率；Yi 為線性回歸方程各質(zhì)量濃度比的實(shí)際相對(duì)響應(yīng)的峰面積比；aXi+b為依據(jù)線性回歸方程計(jì)算得到的理論峰面積比.

2.3" 供試品測(cè)定引入的相對(duì)標(biāo)準(zhǔn)不確定度

2.3.1" 供試品稱量引入的相對(duì)不準(zhǔn)確度urel（sd）

稱取5.0 g菊花樣品使用的電子天平最大秤量值為220.0 g，最小秤量值為0.001 g，實(shí)際分度值0.000 01/0.000 1 g，檢定分度值0.001 g.依據(jù)JJG 1036-2008《電子天平》［18］規(guī)定該天平最大允許誤差為d=±0.05 g，則按三角分布（B類評(píng)定），取k=6，公式如下：

urel（sd）=dk×m.（6）

由供試品稱量引入的相對(duì)標(biāo)準(zhǔn)不確定度urel（sd）為0.004 082.

2.3.2" 供試品提取引入的相對(duì)標(biāo)準(zhǔn)不確定度urel（se）

在1.5節(jié)供試品提取過(guò)程中，不確定度主要玻璃量器量取、定容和環(huán)境溫度變化引入.其中，使用50 mL量筒量取50 mL乙腈2次，容量誤差為±1 mL；10 mL容量瓶（A級(jí)）定容1次，容量誤差為±0.02 mL.參照2.1.1的計(jì)算方法，由式（1）計(jì)算urel（Vse1） = 0.023 094，urel（Vse2） = 0.001 155；由式（2）計(jì)算使用量筒量取及容量瓶定容乙腈溶劑時(shí)由溫度變化引入的相對(duì)標(biāo)準(zhǔn)不確定度分別為urel（tse1） = 0.007 910，urel（tse2） = 0.003 955.

合成urel（Vse）和urel（tse），則由供試品提取引入的相對(duì)標(biāo)準(zhǔn)不確定度為

urel（se）=u2rel（Vse1）+u2rel（Vse2）+u2rel（tsel）+u2rel（tse2）=0.024 756.

2.3.3" 供試品進(jìn)樣體積引入的相對(duì)標(biāo)準(zhǔn)不確定度urel（siv）

實(shí)驗(yàn)使用氣相色譜-質(zhì)譜聯(lián)用儀進(jìn)行測(cè)定，由儀器引入的不確定度來(lái)源主要為配備的微量注射器，其容積為10 μL，相對(duì)標(biāo)準(zhǔn)偏差（RSD）為±1%，供試品進(jìn)樣量為1 μL，采用B類評(píng)定，按矩形分布（k=3），則供試品進(jìn)樣體積引入的相對(duì)不確定度為

urel（siv）=RSDk=0.005 774.（7）

2.3.4" 回收率引入的相對(duì)標(biāo)準(zhǔn)不確定度urel（r）

回收率作為測(cè)量結(jié)果不確定度的重要分量之一，能夠反應(yīng)儀器的準(zhǔn)確度和精密度，以及分析方法是否適合被測(cè)供試品，確保分析數(shù)據(jù)準(zhǔn)確、可靠.在菊花供試品中添加1.4配制的加標(biāo)用農(nóng)藥溶液進(jìn)行加樣回收實(shí)驗(yàn)，添加水平為低（0.01、0.015、0.025 mg/kg）、中（0.02、0.03、0.05 mg/kg）和高（0.04、0.06、0.1 mg/kg）3個(gè)水平，每個(gè)水平重復(fù)測(cè)定3次，計(jì)算回收率.采用A類不確定度評(píng)定，由式（8）計(jì)算回收率r的標(biāo)準(zhǔn)偏差s（r），式（9）計(jì)算回收率r引入的標(biāo)準(zhǔn)不確定度，式（10）計(jì)算回收率r引入的相對(duì)標(biāo)準(zhǔn)不確定度.如表3所示.從表3中可知：本文所用方法對(duì)菊花中29 種農(nóng)藥殘留測(cè)定的平均回收率為69.85%～93.26%，相對(duì)標(biāo)準(zhǔn)偏差為4.15%～10.58%，說(shuō)明該方法準(zhǔn)確可靠.

s（r）=∑ni=1（ri-r）2n-1，（8）

u（r）=s（r）n，（9）

urel（r）=u（r）r.（10）

2.3.5" 測(cè)定重復(fù)性引入的相對(duì)不確定度urel（rd）

在相同條件下，平行被測(cè)供試品在重復(fù)測(cè)量中測(cè)量值會(huì)存在變化.實(shí)驗(yàn)對(duì)3份同一菊花空白試樣中添加高水平的29 種混合標(biāo)準(zhǔn)儲(chǔ)備液農(nóng)藥的殘留量進(jìn)行測(cè)定.采用A類不確定度評(píng)定，按式（11）計(jì)算測(cè)量值的標(biāo)準(zhǔn)偏差s（rd），按式（12）計(jì)算平均值的標(biāo)準(zhǔn)偏差s（rd），按式（13）計(jì)算測(cè)量方法重復(fù)性引入的相對(duì)標(biāo)準(zhǔn)不確定度urel（rd），結(jié)果為0.010 651～0.072 800.

s（rd）=∑ni=1（rdi-rd）2n-1，（11）

s（rd）=s（rd）n，（12）

urel（rd）=s（rd）rd.（13）

2.4" 合成相對(duì)標(biāo)準(zhǔn)不確定度

由上述不確定度來(lái)源得到的各分量相對(duì)標(biāo)準(zhǔn)不確定度按式（14）計(jì)算，得到氣相色譜-串聯(lián)質(zhì)譜方法測(cè)定菊花中29 種農(nóng)藥殘留測(cè)定的合成相對(duì)標(biāo)準(zhǔn)不確定度urel（s），結(jié)果見(jiàn)表4.同時(shí)可以推算出菊花中29 種農(nóng)藥殘留測(cè)定各不確定度分量在合成標(biāo)準(zhǔn)不確定度中的平均占比情況，見(jiàn)圖1.

urel（s）=u2rel（sr）+u2rel（sw）+u2rel（sc）+u2rel（sd）+u2rel（se）+u2rel（siv）+u2rel（r）+u2rel（rd）.（14）

2.5" 擴(kuò)展不確定度和測(cè)量結(jié)果表示

依據(jù)JJF 1135—2005《化學(xué)分析測(cè)量不確定度評(píng)定》，采用置信概率p=95%，包含因子k=2，按式（15）計(jì)算擴(kuò)展不確定度U.

U=k×x×urel（s）.（15）

因此，采用氣相色譜-串聯(lián)質(zhì)譜法對(duì)菊花中29 種農(nóng)藥殘留量進(jìn)行測(cè)定，當(dāng)測(cè)量數(shù)值為0.04、0.06、0.1 mg/kg時(shí)，得到擴(kuò)展不確定度為0.004 787～0.058 935 mg/kg，具體結(jié)果見(jiàn)表4.

3" 結(jié)論

采用乙腈提取菊花中29 種農(nóng)藥殘留，經(jīng)HLB固相萃取小柱凈化，氣相色譜-串聯(lián)質(zhì)譜內(nèi)標(biāo)法測(cè)定農(nóng)藥的含量，29 種農(nóng)藥在0.002～0.1 mg/L內(nèi)線性關(guān)系良好，回歸方程決定系數(shù)R2為0.997 36.加標(biāo)回收實(shí)驗(yàn)的平均回收率為69.85%～93.26%，說(shuō)明該方法準(zhǔn)確度、精密度良好，滿足結(jié)果需要.本文分析了實(shí)驗(yàn)過(guò)程中影響該方法不確定度的各個(gè)分量.通過(guò)對(duì)各分量不確定度占比的統(tǒng)計(jì)可知，該方法不確定度主要來(lái)自農(nóng)藥標(biāo)準(zhǔn)品溶液的配制和標(biāo)準(zhǔn)曲線擬合，平均占比分別為66.69%和45.67%；其次為測(cè)量重復(fù)性和供試品提取分別為39.64%及28.02%，回收率次之，為10.58%，而供試品稱量及供試品進(jìn)樣對(duì)不確定度的影響較小，只有4.62%和6.54%.其中，混合農(nóng)藥標(biāo)準(zhǔn)品儲(chǔ)備液配制不確定度貢獻(xiàn)率8.03%.因此，在實(shí)際實(shí)驗(yàn)分析中應(yīng)盡量簡(jiǎn)化配制步驟，減少玻璃量器及可調(diào)移液器的使用次數(shù)或者更換精度更高的量具，提高實(shí)驗(yàn)人員的操作規(guī)范性.對(duì)于標(biāo)準(zhǔn)曲線擬合，可以適當(dāng)增加參與標(biāo)準(zhǔn)曲線繪制的系列濃度點(diǎn)的個(gè)數(shù)及標(biāo)準(zhǔn)溶液測(cè)定的次數(shù)，同時(shí)確保供試品測(cè)定的結(jié)果落在標(biāo)準(zhǔn)曲線的線性范圍內(nèi)，以減小不確定度來(lái)保證檢測(cè)結(jié)果數(shù)據(jù)的科學(xué)性、嚴(yán)謹(jǐn)性.

參" 考" 文" 獻(xiàn)：

［1］" 王德勝，黃艷梅，石巖，等.菊花化學(xué)成分及藥理作用研究進(jìn)展［J］.安徽農(nóng)業(yè)科學(xué)， 2018， 46（23）： 9-11.DOI： 10.13989/j.cnki.0517-6611.2018.23.003.

［2］" GU J， SCOTTI F， REICH E， et al. Chrysanthemum species used as food and medicine： Understanding quality differences on the global market［J］. South African Journal of Botany， 2022， 148（8）： 123-134.DOI： 10.1016/J.SAJB.2022.04.009.

［3］" 楊丹，王珊，胡云飛，等.菊花農(nóng)藥殘留研究進(jìn)展［J］.中國(guó)中藥雜志， 2021， 46（6）： 1339-1344.DOI： 10.19540/j.cnki.cjcmm.2020.1123.602.

［4］" JIANG M， ZHANG W， ZHANG T， et al. Assessing transfer of pesticide residues from chrysanthemum flowers into tea solution and associated health risks［J］. Ecotoxicology and Environmental Safety， 2020， 187（1）： 109859.DOI： 10.1016/j.ecoenv.2019.109859.

［5］" 馬力輝，崔炎，李迅，等.HPLC法測(cè)定遠(yuǎn)志中黃曲霉毒素B1含量的不確定度評(píng)定 ［J］.河北大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2023， 43 （5）： 484-491.DOI： 10.3969/j.issn.1000-1565.2023.05.006.

［6］" SAITO-SHIDA S， SAITO M， NEMOTO S， et al. GC-MS/MS method for determining quizalofop ethyl， quizalofop tefuryl， and their metabolites in foods［J］. Journal of Food Composition and Analysis， 2023， 115（1）： 105011.DOI： 10.1016/J.JFCA.2022.105011.

［7］" TANKIEWICZ M， BERG A. Improvement of the QuEChERS method coupled with GC-MS/MS for the determination of pesticide residues in fresh fruit and vegetables［J］. Microchemical Journal， 2022， 181（10）： 107794.DOI： 10.1016/J.MICROC.2022.107794.

［8］" 韓梅，侯雪，邱世婷，等.加速溶劑萃取-氣相色譜-串聯(lián)質(zhì)譜法測(cè)定菊花中多種農(nóng)藥殘留［J］.食品安全質(zhì)量檢測(cè)學(xué)報(bào)， 2019， 10（10）： 3065-3077.

［9］" 李俊，蔡滔，代虹鏡，等.基于氣相色譜質(zhì)譜法測(cè)定茶葉中噠螨靈等9種農(nóng)藥殘留量的測(cè)量不確定度評(píng)定［J］.食品工業(yè)科技， 2022， 43（14）： 311-322.DOI： 10.13386/j.issn.1002-0306.2021090073.

［10］" 殷超，陳榮玲，鹿毅.氣相色譜-質(zhì)譜法測(cè)定葡萄中苯醚甲環(huán)唑的不確定度評(píng)定［J］.食品安全質(zhì)量檢測(cè)學(xué)報(bào)， 2020， 11（18）： 6412-6416.DOI： 10.19812/j.cnki.jfsq11-5956/ts.2020.18.031.

［11］" 陽(yáng)曦，陳慧斐，向世杰.QuEChERS-氣相色譜-串聯(lián)質(zhì)譜法測(cè)定結(jié)球甘藍(lán)中甲基毒死蜱農(nóng)藥殘留不確定度評(píng)估［J］.食品安全質(zhì)量檢測(cè)學(xué)報(bào)， 2021， 12（3）： 1136-1141.DOI： 10.19812/j.cnki.jfsq11-5956/ts.2021.03.056.

［12］" 全國(guó)法制計(jì)量管理計(jì)量技術(shù)委員會(huì).測(cè)量不確定度評(píng)定與表：JJF 1059.1-2012 ［S］.北京： 中國(guó)計(jì)量出版社， 2012： 9-26.

［13］" 全國(guó)物理化學(xué)計(jì)量技術(shù)委員會(huì).化學(xué)分析測(cè)量不確定度評(píng)定： JJF 1135—2005 ［S］.北京： 中國(guó)計(jì)量出版社， 2005：1-12.

［14］" CHEN L， LIU Y， HUANG X， et al. Comparison of chemical constituents and pharmacological effects of different varieties of chrysanthemum Flos in China［J］. Chemistry amp; Biodiversity， 2021， 18（8）： e2100206.DOI： 10.1002/CBDV.202100206.

［15］" 全國(guó)流量容器計(jì)量技術(shù)委員會(huì).常用玻璃量器：JJG 196—2006 ［S］.北京： 中國(guó)計(jì)量出版社， 2007： 7-8.

［16］" 姚毅珍.移液器標(biāo)稱容量及測(cè)量誤差的允許標(biāo)準(zhǔn)［J］.中國(guó)計(jì)量， 2012（7）： 82-83.DOI： 10.16569/j.cnki.cn11-3720/t.2012.07.048.

［17］" PANO-FARIAS NS， CEBALLOS-MAGAA SG， MUNIZ-VALENCIA R， et al. Validation and assessment of matrix effect and uncertainty of a gas chromatography coupled to mass spectrometry method for pesticides in papaya and avocado samples［J］. Journal of Food and Drug Analysis， 2017， 25（3）： 501-509.DOI： 10.1016/j.jfda.2016.09.005.

［18］" 全國(guó)質(zhì)量、密度計(jì)量技術(shù)委員會(huì).電子天平： JJG 1036—2008 ［S］.北京： 中國(guó)計(jì)量出版社， 2007：1-13.

（責(zé)任編輯：梁俊紅）