本发明属于化学检测
技术领域：
，特别是涉及到一种单点内标法定量水产品中多氯联苯含量检测结果不确定度的评估方法。
背景技术：
：多氯联苯(polychlorinatedbiphenyls，简称pcbs)是最具有代表性的一类持久性有机污染物(pops)，具有突出的致癌、致畸、致突变的生物毒性、生物蓄积性和环境持久性，对人类健康危害极大。全球环境监测系统/食品规划部分(gems/food)中规定7种“指示性pcb”单体(pcb28、pcb52、pcb101、pcb118、pcb153、pcb138、pcb180)作为常用的监测pcbs污染情况指标，通常以水产类食品的pcbs总量分析为主要监测指标，美国食品药品管理局(fda)就建议水产品中pcbs最大允许浓度为2ug/g(湿质量)。我国最近发布的《食品安全国家标准食品中污染物限量》(gb2762-2017)中规定，水产动物及其制品中多氯联苯的限量指标为0.5mg/kg(以pcb28、pcb52、pcb101、pcb118、pcb138、pcb153和pcb180总和计)。目前水产品中pcbs含量测定数据通常由气相色谱法和气相色谱—质谱联用法获得，由于水产品基质成分复杂，前处理过程难度大，多组分仪器分析干扰多，这些因素都会对检测结果的质量产生影响。为此，有必要对检测结果质量进行定量评估，以确保监测数据的准确性与可比性，以便及时发现、准确预警、正确分析、科学解决水产品质量安全问题。我国在《检测和校准实验室能力的通用要求》就规定，“当检测方法要求；测量不确定度与检测结果的有效性或应用领域有关；客户提出要求；当测试结果处于规定指标临界值附近时，测量不确定度对判断结果符合性会产生影响时，需要对测量不确定度进行评定，并在检测报告中给出不确定度值”。测量不确定度的大小很大程度上决定了测量结果的可用性。因此，水产品pcbs检测作为一项与环境、食品安全密切相关的重要指标，其不确定度的正确评估、准确计算和科学表达对实现测量数据科学检测、准确度量、合理判定、互认共享以及量值溯源、检测过程关键因素控制、提高检测水平具有重要意义。因此现有技术当中亟需要一种新型的技术方案来解决这一问题。技术实现要素：本发明所要解决的技术问题是：提供单点定量水产品中多氯联苯含量检测结果不确定度的评定方法，通过建立不确定度计算数学模型、分析相关不确定度来源、对各分量合理的量化计算，实现样品处理过程、仪器测量过程以及定量分析过程不确定度的准确分析与合理表达。单点定量水产品中多氯联苯含量检测结果不确定度的评定方法，其特征是：包括以下步骤，且按以下步骤顺次进行，步骤一、采用气相色谱仪测定、单点内标计算水产品中多氯联苯含量；步骤二、建立水产品中多氯联苯含量检测结果不确定度的定量数学模型：式中：x—样品中pcbs含量，μg/kg；cs—标准溶液中目标物质量浓度，ng/ml；cn3—样品溶液中内标物质量浓度，ng/ml；cn2—上机标准溶液中内标物质量浓度，ng/ml；a—样品溶液中目标峰响应面积；an0—上机标准溶液中内标峰响应面积；as—上机标准溶液中目标峰响应面积；an1—样品溶液中内标峰响应面积；v1—样品最终定容体积，ml；m—待测样品称量质量，g；rec—回收率％；frep—重复性系数。步骤三、分析不确定度来源，包括a、测量重复性引入的不确定度u(frep)；b、测量样本回收率引入的不确定度c、测量质量引入的不确定度u(m)；d、测量样本最终定容体积引入的不确定度u(v1)；e、测量色谱峰面积引入的不确定度u(a)；f、测量上机标准溶液中内标物浓度引入的不确定度u(cn2)；g、测量样品溶液中内标物浓度引入的不确定度u(cn3)；h、上机标准溶液中目标物浓度cs55引入的不确定度u(cs55)；步骤四、测量不确定度各分量的评定a、测量重复性引入的不确定度u(frep)所述步骤一进行n次独立重复测量，结果为n次测定的平均值，计算公式为相对标准不确定度其中，单次测量标准偏差按下式计算：式中：xi为每一次独立测量的样品中pcbs含量；为n次独立重复测量样品中pcbs含量的均值；b、测量样本平均回收率引入的不确定度对样品回收率值rmax和最小值rmin相对于最佳回收率100％并不对称并假定矩形分布引入不确定度分量：c、测量质量引入的不确定度u(m)采用ⅰ级电子天平，实际分度值d＝0.1mg，检定分度值e＝10d＝1mg，称量重复偏差和最大误差均取最大值1.0mg，通过矩形分布，换算成标准不确定度为：测定的鱼肉样本称量质量引入的不确定度：urel(m)＝u(m)/m；d、测量样本最终定容体积引入的不确定度u(v1)综合定容所用移液器校准和温度对体积的影响并假定矩形分布得到：u(v1)＝0.0058v1ml,urel(v1)＝0.0058；e、测量色谱峰面积引入的不确定度u(a)综合气相色谱定量重复性并假定矩形分布计算可得：u(a)＝0.0173a；样品溶液中内标物峰面积an1引入的不确定度u(an1)、上机标准溶液中内标物峰面积an0引入的不确定度u(an0)、样品溶液中目标物峰面积a引入的不确定度u(a)、标准溶液目标物峰面积as引入的不确定度u(as)之间的关系为：u(a)＝u(as)＝u(an0)＝u(an1)＝0.0173a；f、测量上机标准溶液中内标物浓度引入的不确定度u(cn2)采用定值浓度标准不确定度u(cn0)的内标标准溶液配制上机内标标准溶液，得到上机标准溶液中内标物浓度引入的不确定度其中g、测量样品溶液中内标物浓度引入的不确定度u(cn3)采用定值浓度标准不确定度u(cn0)的内标标准溶液配制内标物浓度为cn3的测量样品溶液，样品溶液中内标物浓度不确定度h、上机标准溶液中目标物浓度引入的不确定度u(cs55)采用定值浓度cs0、标准不确定度u(cs0)的指示性多氯联苯混合标准溶液配制上机标准溶液，则上机标准溶液中目标物浓度不确定度其中urel(cs0)＝u(cs0)/cs0；步骤五、合成标准不确定度的计算综合所述步骤四中获得的各分量不确定度评估值，合成标准不确定度为步骤六、扩展不确定度的计算水产品样本中多氯联苯含量的测定结果为x，测量结果标准不确定度u(x)＝urel(x)×x，扩展不确定度为：u＝k×u(x)k为包含因子，测量结果在95％包含概率下，k＝2则测量结果表示为水产品样本中多氯联苯含量：x±u，k＝2。所述步骤四b、测量样本平均回收率引入的不确定度中平均回收率与100％存在显著性差异，包含在结果计算中，回收率为不确定度来源应包含在计算结果中。通过上述设计方案，本发明可以带来如下有益效果：单点定量水产品中多氯联苯含量检测结果不确定度的评定方法，通过建立不确定度计算数学模型、分析相关不确定度来源、对各分量合理的量化计算，实现样品处理过程、仪器测量过程以及定量分析过程不确定度的准确分析与合理表达；本发明可以带来的进一步的有益效果在于：1、通过建立不确定度数学评估模型，给出不确定度来源，并按类别采用a类或b类评估方法，实现了水产品多残留检测仪器分析的不确定度量化计算，填补了水产品中药物残留检测量值溯源相关研究的空白，为提供准确可靠的水产品质量安全检测数据，以便于实现更广范围的互认共享奠定基础；2、通过对水产品中pcbs残留测定的不确定度分析，将单点内标定量模式通过数学模型转化为各相关不确定度来源的相对不确定度之和，实现了单点内标定量与单点外标定量不确定度评估模式的统一，提出了单点定量不确定度评估的模块化模式；该模式既统一又简化了相关仪器分析不确定度的量化评估过程；3、校准过程是水产品质量安全仪器分析必须的定量过程，校准溶液的配制是仪器能否准确定量的关键步骤，一般都需经过标准储备液、标准工作液、标准使用液，其不确定度影响因素具有共同特征，通过建立配制过程不确定度公式，使得该步骤的不确定度评估简单化、标准化、系统化，而且给出了该步骤不确定度来源及其灵敏度，更加有利于标准溶液配制过程的质量控制。4、利用b类评估模式定量评估回收率这一不确定度来源，并将t检验作为判定方法，确定水产品药物残留检测回收率是否作为不确定度来源。该方法与人们常用的a类评估方法不同，是建立在以往方法确认或长期积累数据基础上，其适用范围更宽，评估方式更为合理，使得水产品药残检测不确定度评估更加科学、准确。具体实施方式单点定量水产品中多氯联苯含量检测结果不确定度的评定方法，测定原理为：用气相色谱-ecd法测定水产品中的多氯联苯含量，是以pcb198为定量内标，将试样水浴加热震荡提取后，经硫酸净化，采用气相色谱-电子捕获检测器测定，保留时间定性，内标法定量。单点定量水产品中多氯联苯含量检测结果不确定度的评定方法，包括以下步骤：(1)称量：称取粉碎均匀的鱼肉样品，给出称量质量m为5.0g，精确到0.1g。(2)40ng/ml(cn1)内标标准工作液的配制：将购买的浓度为2.0μg/ml(cn0)国家标准物质定级证书中制定编号的pcb198标准溶液1ml(v1.0)全部转移并定容至50ml(v50.0)容量瓶中。(3)样品中内标物的加入：在称量好的鱼肉样品中，用量程为100μl～1000μl的移液器，精密移取内标物标准溶液(cn1)0.1ml(v0.1)，混合均匀后待提取，使得样品最终定容后内标物浓度为0.8ng/ml(cn3)。(4)提取：在上述样品中等比例加入正己烷和二氯甲烷适量，于40℃水浴震荡提取2h，保持震荡速度200r/min，离心后倾出全部上清液，逐渐减压蒸馏至-0.1mpa，然后移取5ml正己烷溶解并转移至刻度管中，获得定容体积为5ml(v5.0)的提取液。(5)净化：在提取液中加入浓硫酸除去杂质，直至提取液清澈透明后，上机测定。(6)测定：注射1μl净化彻底的样品溶液至气相色谱仪gc中，得到目标物峰面积a和内标物峰面积an1。(7)80ng/ml(cs1)混合标准工作液的配制：将购买的定值浓度为2.0ug/ml(cs0)pcb混合标准溶液1ml(v1.0)全部转移并定容于25ml(v25.0)容量瓶；(8)56ng/ml(cs2)混合标准使用液的配制：用移液器(量程1000ul～5000ul)分两次移取浓度cs1标准溶液各3.5(v3.5)ml，并定容至10ml(v10.0)容量瓶中。(9)40ng/ml(cs3)、24ng/ml(cs4)、8ng/ml(cs5)混合标准使用液的配制：用移液器(量程1000ul～5000ul)移取浓度为cs1的标准溶液5ml(v5.0)、3ml(v3.0)、1ml(v1.0)并分别定容至10ml(v10.0)容量瓶中。(10)加入内标的混合标准溶液的配制：用移液器(量程100ul～1000ul)移取步骤(7)、(8)、(9)配制的不同浓度的混合标准溶液各0.8ml(v0.8)，分别与0.2ml(v0.2)内标标准溶液(cn1)混合均匀，得到最终pcb混合标准溶液浓度(cs55)分别为5.6ng/ml并且内标物pcb198最终浓度(cn2)为8ng/ml。(11)用制备的标准溶液作gc校准：注射1μl分别含有8ng/ml(cn2)pcb198和5.6ng/ml(cs55)pcb28,pcb52，pcb101,pcb118，pcb138，pcb153和pcb180混合标准溶液，得到参考标准峰响应值as和内标峰响应值an0。将上述测量结果进行计算公式和不确定度评估数学模型的建立，不确定度主要来源分析，不确定度分量的量化，合成标准不确定度的计算，扩展不确定度的计算以及测量结果的表示。一、测量结果数学模型的建立上述单点内标法定量计算单点定量水产品中多氯联苯含量结果不确定度的评定方法，对于气相色谱法测定鱼类样本中多氯联苯含量，其定量数学模型为：式中：x—样品中pcbs含量，μg/kg；cs—标准溶液中目标物质量浓度，ng/ml；cn3—样品溶液中内标物质量浓度，ng/ml；cn2—上机标准溶液中内标物质量浓度，ng/ml；a—样品溶液中目标峰响应面积；an0—上机标准溶液中内标峰响应面积；as—上机标准溶液中目标峰响应面积；an1—样品溶液中内标峰响应面积；v1—样品最终定容体积，ml；m—待测样品称量质量，g；rec—回收率％；frep—重复性系数。二、不确定度来源分析根据数学模型和检测过程分析不确定度来源包括：测量重复性引入的不确定度u(frep)；样本平均回收率引入的不确定度样本称量质量引入的不确定度um；样本最终定容体积引入的不确定度u(v1)；样品溶液中内标物峰面积引入的不确定度u(an1)；上机标准溶液中内标物峰面积引入的不确定度u(an0)；样品溶液中目标物峰面积引入的不确定度u(a)；标准溶液目标物峰面积引入的不确定度u(as)；上机标准溶液中内标物浓度引入的不确定度u(cn2)；样品溶液中内标物浓度引入的不确定度u(cn3)；上机标准溶液中目标物浓度引入的不确定度u(cs55)。三、不确定度分量的量化1、测量重复性引入的不确定度u(frep)上述单点内标法计算水产品中多氯联苯含量检测结果不确定度的评定方法，在重复性条件下，对同一样本从取样开始重复测量6次，结果是6次测定的平均。按计算测定结果见表1。表1pcbs重复性测定结果2、样本平均回收率引入的不确定度上述单点内标法计算水产品中多氯联苯含量检测结果不确定度的评定方法，样本平均回收率引入的不确定度考虑样品回收率最大值和最小值相对于最佳回收率100％并不对称并假定矩形分布，则测量样本回收率引入的不确定度分量其中为样本测定回收率的均值。回收率是否作为不确定度来源应通过显著性检测来确定，即以平均回收率是否与100％有显著性差异为判定标准。当平均回收率与100％有显著性差异时，则应被明确地包含在结果计算中，回收率应作为明确的不确定度来源。平均回收率的统计检测数据t按下式计算：计算所得t值与95％置信度下，n-1自由度的双边临界值tcrit相比较(其中n是用来评估的测试结果的数目)，假如t大于或等于tcrit值，则与1有显著差异，则应被明确地包含在结果计算中。根据相同操作条件下获得的鱼肉样品三水平回收率试验数据，计算平均回收率、平均回收率不确定度见表2。表2回收率及其不确定度结果首先，对表2每一组分的平均回收率数据进行显著性检测，检测结果见表3。表3显著性检测结果由表3显著性比较结果可知，样品平均回收率与100％均存在显著性差异，需要将平均回收率作为不确定度来源之一进行评定。每一组分平均回收率不确定度见表2。3、称量质量引入的不确定度u(m)对于精度为十万分之一且检定合格的ⅰ级电子天平，当实际分度值d＝0.1mg，检定分度值e＝10d＝1mg，在5×104≤m≤2×105范围内，最大允许称量误差为±1.0e＝±1.0×1mg＝±1mg。称量重复偏差同样满足不超过相应载荷最大允许误差绝对值的要求，取最大值1.0mg。假定为矩形分布，换算成标准不确定度为：则样本称量质量引入的不确定度urel(m)＝um/m＝1.154/5000＝0.0002308。4、样本最终定容体积引入的不确定度u(v1)综合考虑定容所移液器校准和温度对体积的影响并假定矩形分布，试样定容体积引入的不确定度u(v1)＝0.0058v1。5、峰面积引入的不确定度根据气相色谱定量重复性并假定矩形分布，色谱峰面积引入的不确定度u(a)＝u(as)＝u(an0)＝u(an1)＝0.0173a。其中，u(an1)为样品溶液中内标物峰面积不确定度、u(an0)为上机标准溶液中内标物峰面积不确定度、u(a)样品溶液中目标物峰面积不确定度、u(as)标准溶液目标物峰面积不确定度。6、上机标准溶液中内标物浓度引入的不确定度u(cn2)按含有内标物的上机标准溶液配制过程得到数学模型：则内标标准上机液引入的不确定度(此过程v1.0为购买储备液全部转移体积，忽略其不确定度)。其中，与体积相关量值的不确定度一般与定量所用器具校准和溶剂随温度变化引入的不确定度有关，综合考虑这两个因素，其中，单标线容量瓶容量允差以及移液器的允许误差，分别依据jjg646-2006移液器检定规程和jjg196-2006常用玻璃量器检定规程给定数值，假定矩形分布计算校准产生的不确定度。单标线容量瓶和移液器校准温度均为20℃，而实验室温度在20±8℃范围内变化，有机溶剂膨胀系数为1.0×10-3，由此产生的体积变化为：±(v×8×1.0×10-3)＝±0.008vml，按矩形分布有：带入后整理得到：依据国家标准物质定级证书，使用浓度为2.0ug/ml，内标标准储备液不确定度u(cn0)为0.018ug/ml，urel(cn0)＝0.018/2＝0.009。综合上述，7、样品溶液中内标物浓度引入的不确定度u(cn3)根据含有内标的样品处理过程得到样品溶液内标物浓度cn3的数学模型为则样品溶液中内标物引入不确定度为(忽略购买储备液全部转移体积v1.0不确定度)。其中，与体积相关量值的不确定度一般与定量所用器具校准和溶剂随温度变化引入的不确定度有关，综合考虑这两个因素，有同样，根据单标线容量瓶容量允差以及移液器允许误差，以及实验室温度在20±8℃范围内变化引起体积变化±(v×8×1.0×10-3)＝±0.008vml，按矩形分布有：带入整理并计算得到：8、上机标准溶液中目标物浓度引入的不确定度u(cs55)按照上机标准溶液的配制过程得到数学模型为：则目标物上机液引入的不确定度urel(cs55)：(此过程v1.0为购买储备液全部转移体积，忽略其不确定度)。综合考虑与体积相关量值一般与定量所用器具校准和溶剂随温度变化引入的不确定度有关，则并且带入整理并计算得到：其中标准储备液引入的不确定度u(cs0)根据国家标准物质定级证书给定标准值及扩展不确定度，标准不确定度为扩展不确定度除以扩展因子2，相对标准不确定度urel(cs0)＝u(cs0)/cs0，以及标准上机测定液相对标准不确定度见表4。表4标准储备液不确定度u(cs0)四、合成标准不确定度的计算表8水产品中多氯联苯含量相对测量不确定度分量表来源符号pcb28pcb52pcb101pcb118pcb138pcb153pcb180cs55urel(cs55)0.01870.01860.01880.01860.01890.01860.0187cn3urel(cn3)0.01700.01700.01700.01700.01700.01700.0170cn2urel(cn2)0.0190.0190.0190.0190.0190.0190.019aurel(a)0.01730.01730.01730.01730.01730.01730.0173an0urel(an0)0.01730.01730.01730.01730.01730.01730.0173asurel(as)0.01730.01730.01730.01730.01730.01730.0173an1urel(an1)0.01730.01730.01730.01730.01730.01730.0173v1urel(v1)0.0580.0580.0580.0580.0580.0580.058murel(m)0.000230.000230.000230.000230.000230.000230.00023recurel(rec)0.04510.04680.05340.05170.05590.06180.0917frepurel(frep)0.02690.02060.01690.01310.01670.02620.0219将表8中各不确定度分量按下式合成，气相色谱法测定、单点内标法定量的水产品中多氯联苯含量不确定度评估结果见表9。若水产品样本中多氯联苯含量的测定结果为x，则u(x)＝urel(x)×x扩展不确定度的计算：u＝k×u(x)(95％包含概率下，取包含因子k＝2)测量结果表示：水产品样本中多氯联苯含量：x±u，k＝2。表9不确定度评估结果urel(x)x/ng/gu(x)/ng/g扩展不确定度upcb280.09120.3720.03390.0678pcb520.09040.3420.03090.0618pcb1010.09330.3390.03160.0632pcb1180.09170.4170.03820.0764pcb1380.09470.3690.03490.0698pcb1530.10030.3430.03440.0688pcb1800.12020.3560.04280.0856当前第1页12