气相色谱法测定大白菜中百菌清农药残留的不确定度评定

摘 要：目的：评定气相色谱法测定大白菜中百菌清农药残留的不确定度。方法：参照《蔬菜和水果中有机磷、有机氯、拟除虫菊酯和氨基甲酸酯类农药多残留的测定》（NY/T 761—2008）标准测定大白菜中百菌清农药残留，并建立数学模型，分析不确定度的来源。结果：6次平行添加测定大白菜中百菌清平均含量为0.187 4 mg·kg-1，置信概率为95%时，测定的扩展不确定度为0.018 4 mg·kg-1（k=2）。结论：大白菜中百菌清农药残留测量结果引入的不确定度主要有标准曲线拟合、试剂移取、样品重复测定、标准物质、定容和称量等。其中，标准曲线拟合引入的不确定度最大，样品测量引入不确定度最小。

关键词：农药残留;不确定度;风险评估;农产品安全

Evaluation of Uncertainty in the Determination of Chlorothalonil Pesticide Residues in Chinese-Cabbage by Gas Chromatography

WANG Lianguang， WU Xiangzhu*， ZHOU Shuya， JI Qingmei， MA Peng， SHI Anxian

（Zhaotong Green Food Development Center， Zhaotong 657000， China）

Abstract： Objective： To evaluate the uncertainty evaluation of chlorothalonil pesticide residues in Chinese cabbage by gas chromatography. Method： According to the NY/T 761—2008 standard， the pesticide residues of chlorothalonil in Chinese cabbage were determined， and a mathematical model was established to analyze the source of uncertainty. Result： The average content of chlorothalonil in Chinese cabbage was 0.187 4 mg·kg-1 by 6 parallel additions. When the confidence probability was 95%， the expanded uncertainty of the determination was 0.018 4 mg·kg-1（k=2）. Conclusion： The uncertainty introduced in the measurement results of chlorothalonil pesticide residues in Chinese cabbage mainly includes standard curve fitting， reagent pipetting， repeated sample determination， standard material， constant volume and weighing. Among them， the uncertainty introduced by the standard curve fitting is the largest， and the uncertainty introduced by the sample measurement is the smallest.

Keywords： pesticide residues; uncertainty; risk assessment; agricultural product safety

不确定度为测量不确定度的简称。1927年，德国物理学家维纳海森堡提出了测量的不确定性概念，表示测量结果的分散程度，亦是分析结果准确性的可疑程度，习惯俗称测量的“不可靠程度”，用于表达分析测量结果优劣的一种方式，通常在测量结果完整表达中包括测量不确定度，检验检测机构通用要求规定检验检测出现检测临界值时也要给出不确定度[1]。不确定度分为标准不确定度和扩展不确定度两大类。标准不确定度分为A类、B类标准不确定度。A类是用统计分析方法计算，B类是以估计的标准偏差表示，合成标准不确定度即为各类A、B类标准不确定度的合成。扩展不确定度为总不确定度，是合成标准不确定度与置信因子k的倍数，置信因子k与概率发布关系为正态分布，k=2～3。其中，置信概率P=95.45%时，k=2;P=99.73%时，k=3;三角分布，k=;平均分布，k=等。

百菌清作为一种广谱、保护性杀菌剂，被广泛应用于蔬菜、水果真菌侵染的防治。昭通大白菜是当地的经济特色产业，但其霜霉病的发生严重影响其产量和质量。百菌清对大白菜霜霉病的防治有重要的作用，不仅可促进农民增收，还可为昭通产业发展提供指导。但百菌清农药的不合理使用或过量使用，使得大白菜中残留增大。为测量其在大白菜中的含量的测量不确定度，对有效监管大白菜中百菌清残留提供更加可靠的结果。本实验参考NY/T 761—2008[2]标准方法，在大白菜样品中加入一定量的百菌清标准液，通过提取、净化和使用气相色谱测定以及建立不确定度测量模型，计算其不确定度值。

1 材料与方法

1.1 材料与试剂

大白菜，购买于当地;百菌清（1 000 mg·L-1），农业农村部环境保护科研监测所（天津）;正己烷、二氯甲烷（色谱纯），赛默飞世尔科技（中国）科技有限公司;弗罗里矽柱（1 g/6 mL），天津博纳艾杰尔科技有限公司;乙腈（分析纯），国药集团化学试剂有限公司;其他试剂均为分析纯。

1.2 仪器与设备

安捷伦7890B气相色谱仪配uECD检测器，美国安捷伦公司;旋转蒸发RE-2000B，上海亚荣生化仪器厂;电子天平DT-500B，常熟市金羊砝码仪器有限公司;匀浆机IKA T18，艾卡有限公司;固相萃取装QSE-240011，得泰仪器科技有限公司;MFV-36智能氮吹仪，得泰仪器科技有限公司。

1.3 实验方法

1.3.1 样品前处理

准确称取15.0 g大白菜均质样品，添加0.2 mg·kg-1百菌清，摇匀静止30 min，加入30.0 mL乙腈，均质2 min，过滤到4～5 g NaCl的100 mL具塞量筒中，振荡摇匀，静置30 min，量取10.00 mL上清液于150 mL旋转蒸发小瓶中，在40 ℃旋转蒸发至进干后，加入2 mL正己烷，用弗罗里矽柱净化待测。

1.3.2 色谱条件

色谱柱：DB-5石英毛细管柱（30 m×0.32 mm×0.25 µm）;载气：氮气（99.999%），恒流;流量：3.0 mL·min-1;进样口温度：220 ℃;检测器温度：300 ℃;柱温：初温150 ℃保持2 min，以8 ℃·min-1升至260 ℃，保持7.25 min;进样方式：不分流进样。

1.3.3 校准标准曲线溶液配制

用正己烷稀释1 000 mg·L-1的百菌清配制成100 mg·L-1的标准储备液，继续稀释成10 mg·L的标准工作液，再用正己烷溶液配制成0.01 mg·L-1、0.05 mg·L-1、0.10 mg·L-1、0.50 mg·L-1和1.00 mg·L-1的系列百菌清校准标准溶液。

2 结果与分析

2.1 测量不确定度主要来源

测定大白菜中百菌清农药残留不确定度主要来源于样品的称量、试剂的体积受环境温度变化、量筒和移液管自带的误差、标准物质不确定度和仪器重复测定的随机误差等。

2.2 建立测量不确定度数学模型

样品中被测农药残留量以质量分数ω计，单位以mg·kg-1表示，公式为

ω=V1×A×V3×ρ/V2×AS×m（1）

式中：ρ为标准溶液中农药的质量浓度，mg·L-1;A为样品溶液中被测农药的峰面积;AS为农药标准溶液中被测农药的峰面积;V1为提取溶剂总体积，mL;V2为吸取出用于检测的提取溶液的体积，mL;V3为样品定容体积，mL;m为样品的质量，g。

2.3 测量不确定度评定

2.3.1 天平称量引入的不确定度

该天平在0～50 g称量范围的最大允许误差MPE=±0.05 g[3]，其允许误差的半峰宽为0.05，按均匀分布计算，k=。称量15.0 g样品引入的B类相对标准不确定度urel（a）=MPE/（k×15.0）=0.192 5%。

2.3.2 标准物质引入的相对标准不确定度

标准物质引入的不确定度urel（b）为标准物质证书给定的扩展不确定度以及样品中加入标准物质引入的不确定度。

（1）标准物质纯度引入的相对不确定度urel（b1）。标准物质证书给出的百菌清浓度为1 000 mg·L-1，则标准物质百菌清的扩展不确定度U（b）=1.20 mg·L-1，k=2，则B类相对标准不确定度urel（b1）=U（b）/（k×1 000）=0.060 0%。

（2）标准物质加入样品中引入的不确定度urel（b2）。加入标准物质引入的B类相对不确定度urel（b2）为稀释1 000 mg·L-1百菌清标准溶液配制成100 mg·L-1标准储备液，继续稀释成10 mg·L-1标准工作液，再移取0.3 mL工作液到样品过程中引入的不确定度。

实验室温差变化为（20±5） ℃，正己烷的膨胀系数为1.26×10 -3 mL·℃-1，按均匀分布计算，k=，3次使用1 mL移液管移取1 mL百菌清标准溶液，2次用正己烷定容于10 mL容量瓶，按均匀分布计算，k=，引入的B类相对标准不确定度urel（b2-1）为使用1 mL A级移液管移取3次标准溶液，根据《常用玻璃量具检定规程》（JJG 196—2006）规定，1 mL A级移液管的最大允许误差MPE=±0.007 mL，按三角分布计算[4]，k=，其引入的B类相对标准不确定度urel（b2-2）=MPE×3/（k×1）=0.857 4%。

使用10 mL A 级容量瓶定容标准溶液2次，根据《常用玻璃量具检定规程》（JJG 196—2006）规定，10 mL A级容量瓶的最大允许误差MPE=±0.02 mL，按三角分布计算，k=，其允许误差引入的B类相对标准不确定度urel（b2-3）=MPE×2/（k×10）=1.633 0%。

加入标准物质引入的B类相对合成不确定度urel（b2）为

（3）标准物质引入的不确定度urel（b）。标准物质引入的B类相对标准不确定度为

2.3.3 量取乙腈引入的不确定度

使用100 mL量筒取30.0 mL乙腈引入的B类相对不确定度urel（c）为量筒的允许误差、乙腈体积受实验室温度变化和重复移取乙腈引入。

100 mL量筒量取30 mL乙腈，按《常用玻璃量器检定规程》（JJG 196—2006）规定，100 mL的最大允许误差MPE=±1 mL，按均匀分布计算，k=，则允许误差引入的B类相对标准不确定度urel（c1）=MPE/（k×30.0）=1.924 5%。实验室温差变化为（20±5） ℃，乙腈的体积膨胀系数为1.37×10-3 mL·℃-1，按均匀分布计算，k=，则由乙腈受温度变化引入的B类相对标准不确定度urel（c2）为urel（c2）=1.37×10-3×5/（×30.0）=0.013 2%重复6次移取乙腈的体积分别为38.081 3 mL、30.052 8 mL、30.142 0 mL、30.032 3 mL、30.075 7 mL和30.055 6 mL，标准偏差为0.019 0 mL，则重复移取乙腈引入的A类相对标准不确定度urel（c3）为量取乙腈引入的B类相对合成不确定度urel（c）=。