中国地质学会岩矿测试技术专业委员会、国家地质实验测试中心主办

用质控图和稳健统计-迭代法评估环境检测实验室测量不确定度

狄一安, 孙海容, 孙培琴, 任立军, 刘岩, 周昊, 王婧瑞, 李斯明, 李玉武. 用质控图和稳健统计-迭代法评估环境检测实验室测量不确定度[J]. 岩矿测试, 2014, 33(1): 57-66.
引用本文: 狄一安, 孙海容, 孙培琴, 任立军, 刘岩, 周昊, 王婧瑞, 李斯明, 李玉武. 用质控图和稳健统计-迭代法评估环境检测实验室测量不确定度[J]. 岩矿测试, 2014, 33(1): 57-66.
Yi-an DI, Hai-rong SUN, Pei-qin SUN, Li-jun REN, Yan LIU, Hao ZHOU, Jing-rui WANG, Si-ming LI, Yu-wu LI. Evaluation of Measurement Uncertainty in an Environmental Test Laboratory by Quality Assurance, Control Charting and Robust Statistics[J]. Rock and Mineral Analysis, 2014, 33(1): 57-66.
Citation: Yi-an DI, Hai-rong SUN, Pei-qin SUN, Li-jun REN, Yan LIU, Hao ZHOU, Jing-rui WANG, Si-ming LI, Yu-wu LI. Evaluation of Measurement Uncertainty in an Environmental Test Laboratory by Quality Assurance, Control Charting and Robust Statistics[J]. Rock and Mineral Analysis, 2014, 33(1): 57-66.

用质控图和稳健统计-迭代法评估环境检测实验室测量不确定度

  • 基金项目:
    国家重大科学仪器设备开发专项(2011YQ14014708, 2011YQ17006506);中国合格评定国家认可委员会科技项目(2011CNAS11)
详细信息
    作者简介: 狄一安,高级工程师,长期从事环境监测及实验室质量管理工作。E-mail: dya_62@hotmail.com
    通讯作者: 李玉武,博士,研究员,从事大气颗粒物化学组分表征及来源解析、化学计量学研究。E-mail: liyuwu@cneac.com
  • 中图分类号: X83;O212.1

Evaluation of Measurement Uncertainty in an Environmental Test Laboratory by Quality Assurance, Control Charting and Robust Statistics

More Information
  • 基于实验室长期积累的质控数据评估测量不确定度的方法具有广泛应用前景,但常见的质控图法只能处理单一浓度,而处理多浓度水平的线性校准法建立模型时需要成套、完整的质控数据,不利于基层实验室的应用。稳健统计是指不用识别、剔除离群值,直接应用全部测量数据,将离群值对统计分析结果影响降低到最小的统计分析方法。本文尝试用回收率将不同浓度数据归一化,然后用质控图方法处理。如果存在离群数据时,可用稳健统计法计算期间精密度sR。利用本实验室积累的5套和其他实验室提供的19套环境检测领域常规项目质控数据验证了新方法的可行性。验证结果表明,对单一浓度数据,不经任何处理,稳健统计-迭代法可得到与质控图法基本相符的结果,sR′(相对值)平均偏差为0.15%。对于多浓度水平数据,经归一化后,质控图法、稳健统计-迭代法与线性校准法的结果平均偏差分别为0.43%和0.20%,质控图法与稳健统计-迭代法的结果平均偏差为0.26%,三种方法计算结果基本相符;稳健统计-迭代法更接近于线性校准法计算结果,且方法原理简单,计算步骤明显简化,适用于线性校准法比例模型数据的处理。
  • 加载中
  • 表 1  作者实验室不同检测项目质控样品标称值及测定结果

    Table 1.  Measured items of standard sample and their analysis results in author′s laboratory

    检测项目质控样品标称值不同时间测定结果
    氨氮(mg/L)2.55±0.102.57,2.54,2.57
    0.501±0.0270.524,0.507,0.508,0.497,0.504 0.491,0.504,0.517,0.490 0.513,0.493,0.503,0.523,0.491
    0.778±0.0420.756,0.776,0.764,0.763,0.772 0.766,0.752,0.774,0.757 0.798,0.780,0.764,0.806,0.775
    8.75±0.358.61
    1.22±0.061.22,1.11
    0.425±0.0250.435
    总磷(mg/L) 1.46±0.051.46,1.46,1.45,1.46,1.48 1.47,1.48,1.44,1.44,1.43
    0.539±0.0170.527,0.537,0.536,0.529 0.531,0.540,0.534,0.542
    0.356±0.0210.352,0.361
    COD(mg/L)148±7145,146,145,151,150,145,145
    112±6115,114,111,112,111
    99.9±5.0102,104,100
    61.0±4.361.4,65.0,63.7,61.8
    76.1±5.380.4,78.3,74.6,76.9
    64.3±4.466.5,64.6
    73.5±4.475.4,72.9
    土壤铅(mg/kg)23±321.4,22.4,22.3
    22.6±1.723.9,23.9,22.1,24.0
    27±225.9,26.0,27.2
    30±531.5,28.5,28.9,29.9 32.3,32.1,30.3,33.0
    58±562.6,62.3,60.0
    98±691.8
    314±13319,324,309
    552±29533
    土壤铜(mg/kg)21±220.4
    24.3±1.223.9,23.6,23.3,22.9
    26.3±1.724.3,26.6,25.8,26.1,25.0,25.8
    22.6±1.322.8
    32±231.3
    40±339.0,39.7,39.3,38.9
    144±6140,141,139,139,147
    注: 5套数据编号分别为12~16。
    下载: 导出CSV

    表 2  其他实验室提供的质控数据

    Table 2.  Quality control data from other laboratories collected by author

    编号检测项目数据个数质控样浓度单位及范围数据来源不确定度评估方法
    1COD30500 (mg/L)澳实分析检测(上海)有限公司质控图法
    2三氯乙烯25100 (mg/L)上海市环境监测中心质控图法
    3苯并[a]芘275.0 (mg/kg)澳实分析检测(上海)有限公司质控图法
    4土壤锌2568±8 (mg/kg)上海市环境监测中心质控图法
    5挥发酚450.163 (mg/L)鞍山市环境监测中心站质控图法
    6氟化物301.20 (mg/L)光大水务(济南)有限公司质控图法
    7土壤总氮500.130%±0.010%中国科学院南京地理与湖泊研究所质控图法
    8土壤总氮500.072%±0.009%中国科学院南京地理与湖泊研究所质控图法
    9亚硝酸盐氮230.0500 (mg/L)上海市供水调度监测中心水质监测站质控图法
    10BOD30186 (mg/L)哈尔滨市环境监测中心站质控图法
    11总油3092.0 (mg/L)北京生态岛科技有限责任公司实验室质控图法
    17BOD5022.8~151 (mg/L)上海市环境监测中心线性校准法
    18氯离子4540~500 (mg/L)中石化中原油田环保监测总站线性校准法
    19COD3425~500 (mg/L)澳实分析检测(上海)有限公司线性校准法
    20氨氮341~10 (mg/L)苏州吴中供水有限公司化验中心线性校准法
    21341~10 (μg/L)苏州市自来水公司水质检测中心线性校准法
    22TOC3412.5~200 (mg/L)上海海洋大学船舶压载水检测实验室线性校准法
    23苯并[a]芘2696~169 (μg/kg)上海市环境监测中心线性校准法
    24总烃3442.9~214 (mg/m3)上海市仪表电子工业环境监测站线性校准法
    下载: 导出CSV

    表 3  COD的质控图法计算示例

    Table 3.  Example for evaluation of measurement uncertainty by quality control charting for COD

    序号检测数据s式计算结果 MR式计算结果
    xi|MR|升序排列wi(s)pi1-pn+1-iAiwi(MR)pi1-pn+1-iAi
    11.0459-0.9666-1.630.05140.9858-7.223 -1.460.0720 0.9752 -6.328
    21.01200.03390.9717-1.380.08330.9820-19.502 -1.240.1077 0.9697 -17.181
    31.01400.00200.9730-1.320.09320.9567-27.570 -1.180.1183 0.9376 -24.546
    40.99200.02200.9780-1.080.14030.9356-32.942 -0.970.1670 0.9131 -29.631
    51.00600.01400.9800-0.980.16290.8885-36.073 -0.880.1894 0.8624 -32.826
    60.98000.02590.9800-0.980.16290.8716-42.538 -0.880.1894 0.8451 -38.814
    71.00600.02590.9807-0.950.17110.8672-49.203 -0.850.1975 0.8407 -44.967
    81.03190.02590.9820-0.890.18730.7431-45.510 -0.800.2132 0.7206 -42.308
    90.97800.05390.9820-0.890.18730.7111-49.586 -0.800.2132 0.6909 -46.233
    101.02400.04590.9840-0.790.21420.6472-49.066 -0.710.2391 0.6324 -46.201
    110.98400.03990.9840-0.790.21470.6396-53.741 -0.710.2396 0.6256 -50.637
    121.00400.02000.9846-0.760.22240.6396-58.046 -0.680.2469 0.6256 -54.765
    131.04390.03990.9920-0.410.34260.6057-50.044 -0.360.3583 0.5948 -48.249
    140.98000.06390.9923-0.390.34740.6057-53.670 -0.350.3627 0.5948 -51.780
    151.00780.02780.9949-0.270.39420.5682-51.351 -0.240.4050 0.5611 -50.095
    160.99610.01180.9961-0.210.41700.5411-51.265 -0.190.4255 0.5369 -50.347
    171.00780.01180.9961-0.210.41820.4918-51.104 -0.180.4267 0.4926 -50.501
    181.02350.01570.9974-0.140.44260.4426-48.985 -0.130.4485 0.4485 -48.893
    190.97170.05181.0000-0.020.49180.4182-46.301 -0.020.4926 0.4267 -46.777
    200.99740.02571.0026 0.100.54110.4170-44.991 0.090.5369 0.4255 -45.878
    210.98200.01541.0040 0.170.56820.3942-43.725 0.150.5611 0.4050 -44.977
    220.98070.00131.0060 0.270.60570.3474-39.915 0.240.5948 0.3627 -41.706
    230.99230.01161.0060 0.270.60570.3426-41.440 0.240.5948 0.3583 -43.336
    240.98460.00771.0078 0.360.63960.2224-32.828 0.320.6256 0.2469 -35.376
    250.96660.01801.0078 0.360.63960.2147-33.740 0.320.6256 0.2396 -36.405
    260.99490.02831.0083 0.380.64720.2142-34.488 0.340.6324 0.2391 -37.305
    270.97300.02191.0120 0.560.71110.1873-29.061 0.500.6909 0.2132 -32.305
    281.02570.05271.0140 0.650.74310.1873-27.739 0.580.7206 0.2132 -31.211
    291.00260.02311.0235 1.110.86720.1711-18.814 1.000.8407 0.1975 -22.431
    300.98200.02061.0240 1.130.87160.1629-18.600 1.020.8451 0.1894 -22.322
    311.03600.05401.0257 1.220.88850.1629-18.059 1.090.8624 0.1894 -21.839
    320.99610.03981.0319 1.520.93560.1403-13.721 1.360.9131 0.1670 -17.236
    330.98400.01211.0360 1.710.95670.0932-9.230 1.540.9376 0.1183 -12.369
    341.00000.01601.0439 2.100.98200.0833-7.044 1.880.9697 0.1077 -9.690
    351.00830.00831.0459 2.190.98580.0514-4.627 -1.460.0720 0.0720 -6.889
    平均值1.0000.0261-- A(i)=-1241.7-A(i)=-1236.0
    标准偏差0.02070.0232=sR′- A(s) =0.478-A(MR)=0.315
    数据量(n ) 35--- A*(s) =0.490-A*(MR) =0.322
    下载: 导出CSV

    表 4  稳健统计-迭代法(方法1)计算示例

    Table 4.  Example for robust analysis-algorithm A (method 1)

    序号第1轮第2轮第3轮第4轮第5轮 第6轮
    10.96660.96660.96660.96690.96710.9671
    20.97170.97170.97170.97170.97170.9717
    30.97300.97300.97300.97300.97300.9730
    40.97800.97800.97800.97800.97800.9780
    50.98000.98000.98000.98000.98000.9800
    60.98000.98000.98000.98000.98000.9800
    70.98070.98070.98070.98070.98070.9807
    80.98200.98200.98200.98200.98200.9820
    90.98200.98200.98200.98200.98200.9820
    100.98400.98400.98400.98400.98400.9840
    110.98400.98400.98400.98400.98400.9840
    120.98460.98460.98460.98460.98460.9846
    130.99200.99200.99200.99200.99200.9920
    140.99230.99230.99230.99230.99230.9923
    150.99490.99490.99490.99490.99490.9949
    160.99610.99610.99610.99610.99610.9961
    170.99610.99610.99610.99610.99610.9961
    180.99740.99740.99740.99740.99740.9974
    191.00001.00001.00001.00001.00001.0000
    201.00261.00261.00261.00261.00261.0026
    211.00401.00401.00401.00401.00401.0040
    221.00601.00601.00601.00601.00601.0060
    231.00601.00601.00601.00601.00601.0060
    241.00781.00781.00781.00781.00781.0078
    251.00781.00781.00781.00781.00781.0078
    261.00831.00831.00831.00831.00831.0083
    271.01201.01201.01201.01201.01201.0120
    281.01401.01401.01401.01401.01401.0140
    291.02351.02351.02351.02351.02351.0235
    301.02401.02401.02401.02401.02401.0240
    311.02571.02571.02571.02571.02571.0257
    321.03191.03191.03191.03191.03191.0319
    331.03601.03571.03331.03251.03221.0321
    341.04391.03571.03331.03251.03221.0321
    351.04591.03571.03331.03251.03221.0321
    平均值1.00040.99990.99970.99960.99960.9996
    标准偏差s0.02070.01970.01930.01910.01910.0191
    sR′=1.134×s0.02350.02230.02190.02170.02160.0216
    1.5×sR′0.03530.03340.03280.03260.03250.0324
    -1.5×sR′0.96510.96650.96690.96710.96710.9672
    +1.5×sR′1.03571.03331.03251.03221.03211.0320
    下载: 导出CSV

    表 5  稳健统计-迭代法(方法2)计算示例

    Table 5.  Example for robust analysis-algorithm A (method 2)

    编号x′i-中位值第1轮第2轮第3轮第4轮第5轮
    10.030820.96660.96760.96740.96730.9672
    20.025680.97170.97170.97170.97170.9717
    30.024390.97300.97300.97300.97300.9730
    40.019360.97800.97800.97800.97800.9780
    50.017360.98000.98000.98000.98000.9800
    60.017360.98000.98000.98000.98000.9800
    70.016680.98070.98070.98070.98070.9807
    80.015390.98200.98200.98200.98200.9820
    90.015390.98200.98200.98200.98200.9820
    100.013400.98400.98400.98400.98400.9840
    110.013370.98400.98400.98400.98400.9840
    120.012820.98460.98460.98460.98460.9846
    130.005380.99200.99200.99200.99200.9920
    140.005110.99230.99230.99230.99230.9923
    150.002540.99490.99490.99490.99490.9949
    160.001320.99610.99610.99610.99610.9961
    170.001260.99610.99610.99610.99610.9961
    180.000030.99740.99740.99740.99740.9974
    190.002601.00001.00001.00001.00001.0000
    200.005171.00261.00261.00261.00261.0026
    210.006591.00401.00401.00401.00401.0040
    220.008591.00601.00601.00601.00601.0060
    230.008591.00601.00601.00601.00601.0060
    240.010441.00781.00781.00781.00781.0078
    250.010441.00781.00781.00781.00781.0078
    260.010861.00831.00831.00831.00831.0083
    270.014581.01201.01201.01201.01201.0120
    280.016571.01401.01401.01401.01401.0140
    290.026131.02351.02351.02351.02351.0235
    300.026551.02401.02401.02401.02401.0240
    310.028311.02571.02571.02571.02571.0257
    320.034541.02981.03111.03161.03181.0319
    330.038591.02981.03111.03161.03181.0319
    340.046511.02981.03111.03161.03181.0319
    350.048511.02981.03111.03161.03181.0319
    平均值-0.99930.99950.99960.99960.9996
    标准偏差s-0.01870.01890.01900.01900.0190
    sR′=1.134×s-0.02120.02140.02150.02160.0216
    1.5×sR′0.032420.03170.03210.03230.03230.0324
    -1.5×sR′0.965000.96760.96740.96730.96720.9672
    +1.5×sR′1.029851.03111.03161.03181.03191.0320
    注: 测量数据中位值为0.9974,标准偏差估计值s0=0.01458,s0×1.483=0.0216。
    下载: 导出CSV

    表 6  单一浓度数据归一化前后质控图法及稳健统计-迭代法结果比较

    Table 6.  Comparison of results from quality control charting and robust analysis (algorithm A) based on the same concentration data before and after normalization processing

    数据编号检测项目数据个数 质控图法稳健统计-迭代法可疑数据备注
    ssR′A*(s) A*(MR)sR′
    1COD (mg/L) 30498.95.636.640.270.68498.96.390 归一化前结果
    300.99770.01130.01330.270.680.99770.01280归一化后结果
    30498.95.656.650.270.68498.96.400归一化后结果
    2三氯乙烯(mg/L) 2510412130.310.36105121归一化前结果
    251.040.120.130.310.361.050.121归一化后结果
    3苯并[a]芘(mg/L) 124.990.1410.1450.460.43--0sR′合并前结果
    154.980.1250.1490.410.40--0sR′合并前结果
    274.99-0.145----0sR′合并后结果
    270.9980.0260.0290.720.610.9980.0290归一化后结果
    274.990.1300.1440.720.614.990.1460归一化后结果
    4土壤Zn(mg/kg) 2566.94.04.00.580.5966.94.50归一化前结果
    250.9840.0590.0580.580.590.9830.0660归一化后结果
    2566.94.03.90.580.5966.84.50归一化后结果
    注: 表中测量数据取自CNAS组织的用 “top-down” 技术评估不确定度培训班教材和学员提交的报告。第3套数据前3行为数据归一化前结果。
    下载: 导出CSV

    表 7  单一浓度数据质控图法及稳健统计-迭代法结果比较

    Table 7.  Comparison of results from quality control charting and robust analysis (algorithm A) based on the same concentration data

    数据编号检测项目数据个数 质控图法稳健统计-迭代法可疑数据
    ssR′A*(s) A*(MR)sR′
    5挥发酚(mg/L)450.1630.00270.00290.640.580.1630.00311
    6氟化物(mg/L)301.200.0200.0240.440.491.200.0220
    7土壤总氮(%)500.1290.00370.00340.400.640.1290.00363
    8土壤总氮(%)500.0700.00230.00220.860.850.0700.00223
    9亚硝酸盐氮(mg/L)230.05020.00110.00110.550.580.05020.00122
    10BOD(mg/L)301885.66.60.460.481886.20
    11总油(mg/L)3092.01.21.20.480.5892.01.40
    注: 表中测量数据取自CNAS组织的用“top-down”技术评估不确定度培训班教材和学员提交的报告。
    下载: 导出CSV

    表 8  不同浓度水平质控数据归一化后质控图与稳健统计-迭代法结果比较

    Table 8.  Comparison of results from quality control charting and robust analysis (algorithm A) based on the quality control samples analysis data of different concentration after normalization processing

    数据编号检测项目数据个数 质控图法稳健统计-迭代法RMS可疑数据
    ssR′A*(s)A*(MR)sR′
    12氨氮(mg/L)351.0000.0210.0230.490.321.0000.0220.0232
    13总磷(mg/L)281.0010.0150.0100.302.321.0010.0160.0150
    14COD(mg/L)271.0110.0240.0230.370.511.0110.0260.0260
    15土壤Pb(mg/kg)261.0120.0510.0490.600.711.0110.0570.0510
    16土壤Cu(mg/kg)220.9770.0220.0220.440.430.9770.0220.0321
    下载: 导出CSV

    表 9  不同浓度水平质控数据归一化后三种计算方法结果比较

    Table 9.  Comparison of results from three calculation methods based on the analysis data of different concentration quality control samples after normalization processing

    数据编号检测项目数据个数 质控图法稳健统计-迭代法线性校准法可疑数据
    ssR′A*(s)A*(MR)sR′
    17BOD(mg/L)501.0200.0360.0380.380.421.0190.0350.0302
    18氯离子(mg/L)451.0040.0120.0100.690.811.0040.0120.0120
    19COD(mg/L)341.0090.0210.0270.320.651.0080.0220.0201
    20氨氮(mg/L)34~361.0010.0110.0101.041.311.0010.0130.0140
    21Cd(mg/L)341.0000.0130.0130.460.471.0000.0140.0130
    22TOC(mg/L)340.9960.0110.00751.601.960.9980.00800.0102
    23苯并[a]芘(mg/kg)260.8550.0390.0430.200.270.8560.0400.0371
    24总烃(mg/L)340.9970.0370.0430.410.470.9990.0400.0381
    注: 表中测量数据取自CNAS组织的不确定度培训班学员实习报告。
    下载: 导出CSV
  • [1]

    [2]

    National Association of Testing Authorities, Australia (NATA). Technical Note 33, Guidelines for Estimating and Reporting Measurement Uncertainty of Chemical Test Results [R].

    [3]

    辽宁出入境检验检疫局,中国标准化研究院.GB/Z 22553—2010,利用重复性、再现性和正确度的估计值评估测量不确定度的指南[S].北京:中国标准出版社, 2010.

    [4]

    辽宁出入境检验检疫局,中国合格评定国家认可委员会.GB/T 27407—2010,实验室质量控制-利用统计质量保证和控制图技术评价分析测量系统的性能[S].北京:中国标准出版社, 2010.

    [5]

    辽宁出入境检验检疫局,中国合格评定国家认可中心.GB/T 22554—2010,基于标准样品的线性校准[S].北京:中国标准出版社, 2010.

    [6]

    辽宁出入境检验检疫局,中国合格评定国家认可中心,中国计量科学研究院.GB/T 27411—2012,检测实验室中常用不确定度评定方法与表示[S].北京:中国标准出版社, 2013.

    [7]

    http://www.cnki.com.cn/Article/CJFDTOTAL-WJFX201201003.htm

    [8]

    http://www.cnki.com.cn/Article/CJFDTOTAL-YKCS201202032.htm

    [9]

    http://www.cnki.com.cn/Article/CJFDTOTAL-FXCQ201101015.htm

    [10]

    http://www.cnki.com.cn/Article/CJFDTOTAL-XDSJ2007S1003.htm

    [11]

    http://www.cnki.com.cn/Article/CJFDTOTAL-JILA200208039.htm

    [12]

    Eurachem Workshop.Recent Developments in Measure-ment Uncertainty [R].2011.

    [13]

    The New Draft Eurachem/CITAC Guide: A Summary of Structure and Changes. Quantifying Uncertainty in Analytical Measurement(MU指南第3版) [R].2011.

    [14]

    Van der Veen A M H.Revision of the Guide to the Expression of Uncertainty in Measurement [R].2011.

    [15]

    Magnusson B, Naykki T, Hovind H. NORDTEST Report TR537: Handbook for Calculation of Measurement Uncertainty in Environmental Laboratories [R]. 2003.

    [16]

    European Federation of National Associations of Meas-urement, Testing and Analytical Laboratories(EUROLAB).Measurement Uncertainty Revisited: Alternative Approaches to Uncertainty Evaluation[R].2007.

    [17]

    Ramsey M H, Ellison S L R. Eurachem/EUROLAB/CITAC/Nordtest/AMC Guide: Measurement Uncertainty Arising from Sampling: A Guide to Methods and Approaches[S].Eurachem ISBN 978 0 948926 26 6(http://www.eurachem.org/guides/UfS_2007.pdf).2007.

    [18]

    doi: 10.3969/j.issn.0254-5357.2012.01.029

    [19]

    ISO 13528—2005, Statistical Methods for Use in Profi-ciency Testing by Inter-laboratory Comparisons [S].

    [20]

    Analytical Methods Committee. Robust Statistics: A Method of Coping with Outliers [R].AMC Technical Brief.2001.

    [21]

    http://www.cnki.com.cn/Article/CJFDTOTAL-LHJW201002004.htm

    [22]

    http://www.cnki.com.cn/Article/CJFDTOTAL-IAOB201104003.htm

    [23]

    中国合格评定国家认可委员会.CNACNAS T0476低合金钢中化学成分分析能力验证计划结果报告(附录D) [R].2010.

    [24]

    中国实验室国家认可委员会.环境检测领域基于质控数据评定测量不确定度指南[S].2013.

  • 加载中

(9)

计量
  • 文章访问数:  1065
  • PDF下载数:  6
  • 施引文献:  0
出版历程
收稿日期:  2013-06-03
录用日期:  2013-06-21

目录