2012年第3期 文章编号:1005-3387(2012)03-0054-58 CMC评估和表示的方法初探 胡 畅 邵鹏飞 刘国义 张(1.西安计量技术研究院,西安摘娟 张倩 100011) 710068;2.北京日月驰管理咨询有限公司,北京要:本文论述了利用Excel电子表格对校准和测量能力(CMC)评估及表示的方法和过程。借助Excel软件,对不同被 校准值的不确定度分量进行分析,并利用Excel绘制被测量与其对应测量不确定度拟合曲线。合理选择CMC评估中测量范 围分段和CMC的合理表述。 关键词:CMC;EXcel;线性拟合;非线性拟合 O 引言 随着ILAC-P14:2010(校准领域测量不确定度 的》的实施,目前CNAS认可的部分校准实验室 “校准和测量能力”(以下简称CMC)的表示方式不 插值算法以给出区间内的值的测量不确定度; 3)CMC用被测量值或参数的函数表示; 4)CMC用矩阵表示。此时,不确定度的值取决 于被测量的值以及与其相关的其他参数; 5)CMC用图形表示。此时,每个数轴应有足够 的分辨率,使得到的CMC至少有2位有效数字; 6)CMC不允许用开区问表示(例如“U<X”)。 一能满足上述文件的要求。为此CNAS下发了1 18号 文件,要求2012年5月1日前所有认可的校准实验 室完成CMC核查,未完成或提供的CMC仍不满足 要求的已获认可的校准实验室,CNAS将暂停或撤 销其校准能力的认可。科学、严谨、合理的选择 CMC的表示方式并进行评估是以对不同专业、不同 参数的测量不确定度来源、测量方法的科学分析的 结果。由于测量不确定度是针对被测量而言的,即 针对测量点,而非测量范围。这就意味着对于整个 测量范围或不同量程内每个被测点都应进行测量不 确定度评估。如缺少合理、方便、快捷的计算工具, 将极大地增加实验室的工作量。甚至会影响到 CMC的评估质量。 般情况下,CMC应该用包含概率约为95%的 扩展不确定度表示。CMC的单位应当始终与被测 量一致,或者使用与被测量的单位相关的其他单位 表示,例如用百分比表示。当CMC的单位与被测量 不一致时,应给出必要的说明。 2 CMC的概念中的几个关键点 2.1对“常规条件”的理解 此处的“常规条件”是强调实验室必须根据自 己的实际情况,而不必刻意要求最好的人员、设备以 及环境。应尽量与日常校准保持一致。这样提供客 1 CMC的概念 根据CNAS-CL07:201 1《测量不确定度的要求》 (第二次修订)中规定: 校准和测量能力(CMC)是校准实验室在常规 条件下能够提供给客户的校准和测量的能力。其应 是在常规条件下的校准中可获得的最小的测量不确 定度。应特别注意当被测量的值是一个范围时, CMC通常可以用下列一种或多种方式表示:… 户的CMC才更加真实可靠。 2.2对“最小的测量不确定度”和“可校准的最佳 性能的器具”的理解 “最小的测量不确定度”可以理解为,实验室在 对被校对象的全部准确度等级(以下简称等级)、全 部校准项目、全部校准参量和校准结果的测量不确 定度完成预评估(以下简称MU)后。针对同一测量 点,不同类型、型号的被校对象、使用不同测量方法 和(或)测量设备所有预评估结果的最小值。 若实验室尚未完成完整的预评估,可按下述步 骤完成CMC评估: 1)对所有被校对象按其准确度等级由高到低 1)CMC用整个测量范围内都适用的单一值 表示; 2)CMC用范围表示。此时,实验室应有适当的 ~54一 进行排序,并确定其测量范围。如,三等量块测量范 围为(0.5~100)mm。四等、五等量块测量范围均 为(0.5—500)mm; 2)按准确度等级高低,依次用该等级被校对象 的测量范围减去所有高于本等级被校对象已覆盖的 测量范围所得到的区间作为其CMC评估时的测量 范围。当不同等级被校对象的测量范围相同时,应 选择其中等级最高的被校对象作为CMC的评估对 象。如本文2.2.1例中,先确定三等量块为CMC评 估对象,其测量范围为(0.5—100)mm。再确定四 等量块为CMC评估对象,减去三等量块已经覆盖的 测量范围,其评估时的测量范围为(125~500)mm。 最后选择五等量块作为评估对象,因其测量范围与 四等相同。因此,只须进行预评估,而不在CMC中 体现; 3)在初步确定各等级被校对象CMC评估时的 测量范围后,还应进一步判断在其测量范围内,测量 方法或(和)测量设备是否不同。若不同,则应进一 步对其进行拆分。如在本文2.2.1例中,三等量块 在(0.5—100)mm内,方法和设备均相同,无须拆 分。但如果其测量范围改变为(0.5~500)mm,则 在(0.5~100)mm内,使用干涉仪进行校准,在(125 ~500)mm内。使用测长机进行校准;为了便于对 以后的预评估结果线性拟合,最好分为两段; 4)确定完各等级被校对象及其CMC评估时的 测量范围后,可以利用本文第4部分“CMC的评估 过程”中介绍方法进行MU和CMC评估。 上述方法评估出的CMC,虽然不是来自预评估 的最小值。但由于所选评估对象基本上为“该参量 可校准的最好性能的器具”。从某种程度上也确保 了其为该测量点的“最小的测量不确定度”。 2.3对“CMC”与“MU”的理解 根据CMC的定义,CMC的本质是一个最小的 测量不确定度。而这个最小的测量不确定度是测量 范围内各测量点实验室能够开展校准的所有被校对 象在该点预评估的最小值。因此,CMC源于科学、 完整的MU。是实验室提供给客户的最好测量能 力。这里应注意:①CMC并不是所有准确度等级、 所有型号被校对象的测量不确定度;②校准证书中 各测量点对应的测量不确定度是MU值。而非 CMC值;其一定是大于或等于该测量点对应的CMC 值;③对同一被校仪器同一参数同一被测点,MU可 以有多个。其随着准确度等级、测量原理、测量方 法、测量设备等不同而不同。但CMC却只有一个 值。即所有MU中最小的那个值。 2.4“CMC”与能力验证的关系 CNAS--GL02((能力验证结果的统计处理和能 力评价指南》附录B.4中规定“在En值中使用了由 实验室报告的测量不确定度。如果实验室报告的不 确定度大于其被认可的不确定度,它应该作出解 释” 。目前,多数能验证提供者和测量审核机构 都要求机构提供的测量不确定度不大于被认可的不 确定度(即CMC)。这样就会产生两个问题①、如能 力验证的样品的准确度等级低于CMC评估时所选 择的被校对象的准确度等级。则实验室报告的不确 定度有可能大于其被认可的CMC。因此,应尽量保 持一致。如CMC评估时在某测量范围内使用的三 等量块。则参加能力验证,最好参加三等量块的能 力验证,而不能参加五等量块的能力验证;②、在进 行CMC评估时,应充分考虑到以后能力验证的问 题,可人为合理地放大CMC,以降低风险。 2.5 CMC与被校准对象最大允许误差的关系 CMC的评估结果应对照计量标准的技术指标、 被校对象的技术指标,参考国家计量检定系统 表。一般情况下,CMC不应大于现有最佳仪器(即 可校准的最好性能的器具)在该测量点规程\规范 规定的不确定度的最大值或最大允许误差(以下简 称最大允许误差)的三分之一。 3 CMC评估中应注意的问题 1)CMC评估时选择的校准对象应是对该参量 可校准的最佳性能的器具。 2)CMC的评估是在常规条件下进行的。 3)CMC的评估结果应对照计量标准的技术指 标、被校对象的技术指标,参考国家计量检定系 统表。 ・ 4)CMC的评估应是对全部校准项目、校准参量 和校准结果的测量不确定度的评估。 5)对照国家计量检定系统表给出的同等级计 量标准设备的不确定度要求,可以检查CMC评估值 是否合理。 6)慎用CMC用范围表示的方法,一方面,因为 实验室应有适当的插值算法以给出区间内的值的测 量不确定度。而事实上,插值算法通常难以给出。 另一方面,因为容易造成范围过宽,而造成测量范围 分段不合理,用户不易使用,如表1。 一55— 表1 实验室推荐认可的实验室校准能力范围表 校准 规范代号(含 扩展不确定度 名称 参量 年号)名称 测量范围 和测量能力 (k=2) 数字 质量 数字指示秤 JJG539-1997 lOOmg一 U=0.016mg一 指示秤 检定规程 34700 kg 0.56 kg 4 CMC的评估过程 4.1确定测量范围内的评估点 尽管CMC源于对整个测量范围、全部校准项 目、校准参量和校准结果的测量不确定度的MU。 但实际上并不需要对每个等级、每种型号、每个测量 点逐一完整的预评估。可参考以下步骤确定待评估 的测量点。 1)参考CNAS公布的《校准实验室CMC表示方 式参考范例》(以下简称范例)。选取与之相一致的 测量范围分段原则、测量点选择及CMC表示方式。 2)若“范例”采用单一值表示,预评估中,可针 对不同等级的被校对象。选取其测量范围内的最 大、最小和中间值进行预评估。并对其是否为单一 值进行验证。 3)若“范例”采用函数表示,预评估中,依次判 断各等级被校对象在整个测量范围内(不用减去高 于该等级所有仪器已经覆盖的测量范围)。是否存 在测量方法或(和)测量设备不同,若存在,需进一 步拆分。若相同,则在规程\规范“计量性能要求” 中寻找该等级被校对象在该测量范围下,最大允许 误差对应的测量范围的所有分段。并对这些分段的 最大、最小和中问值进行预评估。最后利用3.3介 绍的方法得到散点图、拟合曲线和拟合方程。如本 文2.2.1例中,四等量块测量范围为0.5~500mm。 由于在此区间内存在选用测量方法和(或)测量设 备不同,须进一步分为0.5~lOOmm和125~500mm 两段。查量块检定规程,在125~500mm内,最大允 许误差对应测量范围分段为100~150mm,150~ 200mm,200~250mm,250~300mm,300~400mm, 400~500mm,取其每段的最大、最小和中间值。得 到125ram,150mm,175ram,200mm,225mm,250mm, 275mm,300mm,350ram,400mm,450mm,500mm。其 中225mm,350mm两个点在日常校准中不会出现, 可以舍弃。另外,由于“范例”中采用函数表示。因 此,可以在上述12个值中对称选择5~7个点进行 预评估和线性拟合(见表2),以减少预评估工作量。 一56— 4)若“范例”中采用范围表示或“范例”中未包 含待校准对象。可在最大允许误差测量范围各分段 尽可多的选择测量点进行预评估。并观察、分析能 否进行线性拟合、曲线拟合或多段线性拟合。若无 法拟合,则需对所有测量点进行预评估。 4.2测量不确定度分量汇总表分析 测量范围内的评估点确定后,由于同一测量范 围区间内,所使用的测量方法和测量设备基本相同。 因此,这些测量点的测量不确定度来源及分量也同 样大体相同,根据《校准和测量能力(CMC)的表示 方式应用指南》(以下简称CMC表示指南)第三部 分“对CMC评估和应用的建议”的要求,可先完成 一个典型点的测量不确定度评估,并对其测量不确 定度分量汇总表进行如下分析、归纳其关联,以及在 不确定度分量数值上的关系。已达到简化处理的效 果。当这些关系确定后,测量范围内其他各点和其 他等级被校量值的测量不确定度也就可以容易的计 算出来。 1)对所有测量不确定度来源分量按从大到小 排序。并按大小顺序排序 ≥tt:≥ ,≥…≥u 。 1 2)判断第一大分量 。是否满足÷“ ≥ 厂i——一 ^/∑(y),若满足,则“。为占绝对优势的不确定度 分量。此时,CMC的表示方式宜与“。的表达方式 一致。特别是对于计量标准最大允许误差为最大分 量时,CMC的表示方式宜与该最大允许误差的表达 方式一致,可以是单一值也可以是函数。 3)若不满足,则进一步寻找是否存在对校准结 果的不确定度贡献可忽略不计的分量,即满足 ≤ Ui1 。若存在, 以后的分量均可忽略不计。特别 上U 是在预评估中,当计量标准器具的准确度等级(或 其示值的测量不确定度)高于被校计量器具最大允 许误差的一个数量级时,计量标准器具的最大允许 误差(或其示值的测量不确定度)对校准结果的不 确定度贡献可忽略不计。 4)依据“CMC表示指南”的要求,测量重复性 和被校对象读数误差对测量不确定度的贡献存在重 复,因此,这两个分量在计算合成不确定度时,只取 其中最大的。… 5)若不存在可忽略分量,应进一步分析各分量 中是否存在不随被校准对象、测量点和测量设备不 同而变化的分量。若存在则以不同颜色予以标注。 如下分析处理: 在分析不变分量时,应特别注意区分计量标准的重 复性和被测对象的重复性。计量标准的重复性是设 备的固有属性,可以认为是一个固定值,但被测对象 1)判断所有评估点的预评估结果的绝对或相 对值是否相等或近似(即 ttmin ≤ ),若相等,进 lU 一的重复性受其准确度等级和测量点的影响,很可能 发生变化。因此,不能简单认为重复性是一成不变 的。除非重复性分量为可忽略,或在测量范围内所 有评估点计算得到的重复性均相等。 4.3预评估结果分析处理 步判断规程\规范“计量性能要求”中最大允许误 差是否存在驻点。若存在,则须增加驻点对应的测 量点进行预评估。若不存在,则用单一值表示该测 量范围内的CMC。并停止该等级仪器的预评估。 2)若不相等或近似,可进行如下分析: ①利用EXCEL散点图制作预评估结果随被测 量变化的曲线Y 。 评估结果判断 对依据本文3.1所确定被校对象测量范围内的 所有评估点进行预评估,将其结果列入表2,并进行 表2 预评估结果分析表 预评估得到的 评估点 测量不 最小二乘法线性拟合后 (mm) 确定度 的测量不确定度 ( m) (O.381xm+O.8x10 ) 125 20o 300 400 5o0 0.40 0.59 0.68 O.71 0.73 0.48 0.54 0.62 0.70 0.78 规程\规范中规定的 不确定度允许值 MPE MPE 是否满足 是否满足 是否满足 是否满足 J (O.201 ̄m+2x10 £) ( m) (≤ Ul≤ ≤峨 ≤ wm) 0.45 O.60 0.80 1.00 1.20 0.53 0.67 0.93 1.20 1.6 2.O 2.8 3.6 、/ 、/ 、/ 、/ 、/ 、/ 、/ 、/ 、/ 、/ × 、/ 、/ 、/ 、/ 、/ 、/ 、/ 、/ 、/ 1.47 4.4 ②从EXCEL图表栏中选择添加趋势线或在绘 图区中鼠标指向某一数据点来选择添加Y 的趋势线 。③依据规程\规范“计量性能要求”,计算各评 估点对应的三分之一最大允许误差眠和最大允许 误差 。 在趋势线格式“类型”项上的趋势预测/回归分析 类型(线性、对数、多项式、乘幂、指数、移动平均)时, 一④若规程\规范中给出了各评估点不确定度的 允许值。如在JJG146-2003.附录A.2中就规定了 方面要观测Y 的图形特点,如是否呈线性或近似 线性或在某些段内呈线性等;另一方面应依据行业惯 例,并尽量与规程\规范中最大允许误差或被测量不 确定度的允许值的数学模型(线性或多项式)保持一 致。选择合适的趋势预测/回归分析类型后,在趋势 线格式“选项”项上勾选“显示公式”即可获得Y,的拟 各等级量块长度测量的不确定度允许值的计算公 式,其中4等量块长度测量的不确定度允许值的计 算公式为(0.20 ̄m+2×10 L),利用该计算公式在 表2中计算出各评估点不确定度允许值 。 ⑤分别对U 、 与 、 间的大小关系进行 判断,以合理选择CMC的表达方式。(见表2) 合公式。最后,利用该拟合公式在表2中计算出各测 量点拟合后计算出的测量不确定度 。 CMC汇总图 ⑥在本文4.3.2.①已形成的散点图的基础上 增加不确定度允许值 和三分之一最大允许误差 v4 百 e 的折线图Y 和,,4。 4l}l{ I嘣 I Y3=0.002 ̄+0. 一/ 3)综合考虑a、预评估结果随评估点变化的曲 线Y。的变化趋势Ib、U。及£,2对 、 的关系;c、 行业惯例;d、CNAS公布的《校准实验室CMC表示 0.37 厂 /讧O.fl ̄x+一 1 .—・— y1 方式参考范例》及其他权威资料,选择合理的CMC 表述及测量范围分段方式。如图1中,Y.基本符合 线性,并且规程\规范中给出的不确定度允许值也是 一髂 0 50 lO0 150 200 250 300 350 4O0 450 500 550 被测量(mm) 图1测量不确定度随被测量变化趋势图 个线性模型。因此,宜选择线性模型进行表述。 另外,观察到Y 的部分数值大于Y,,即通过拟合方 一57— 程计算出的不确定度大于规程\规范中规定不确定 度允许值。但实际评定出的测量不确定度却小于规 定不确定度允许值,因此,不易采用Y 作为表述方 8 6 4 1 l 1 1 2 100mg~34700 kg的分段,显然不太合适,可以按 mg、g、kg进行分段。②同一个范围内CMC最大值 不宜超过最小值的5~10倍。③要确保CMC的范 8 6 4 2 ( 0 0 O O 式。综合上述考虑,可选择Y (即U=0.20txm+2× 10 L)作为CMC的表述方式。 4)若预评估结果随评估点变化呈规律曲线变 化,如图2。在EXCEL选择趋势预N/回归分析类 型上应选择“多项式”,并选择适当的“阶数”,以便 得最接近的拟合曲线。如图2中,拟合后得到Y= 0.0002x 一0.0079x+0.2895,即 =2x10一 L 一7.9  ̄10一 L+0.29IJ,m。 cMc汇总图 :0 0002x2-0 0079x}0 2895 言 i ~ 、 …“ 露 .9 您 O O O 0 0 O O 0 0 雕 9 8 7 j6 5 .5 4 3 2 1 O 。一 一0_22’0 0 10 20 30 40 50 60 7O 8O 90 100 被测量(mm) 图2预评估结果随评估点变化趋势图 5)通常情况下,CMC用二次方程表示较少。选 用前,一方面应考虑规程\规范中最大允许误差或不 确定度允许值是否也采用二次方程,另一方面要符 合行业惯例。如规程\规范或行业惯例中最大允许 误差或不确定度允许值均采用线性模型。当这两情 况发生时,可将二次曲线按其变化趋势,转化为多段 线性拟合(见图3)。如图3中,在测量范围l0~ 20mm内,采用单一值,得到U=0.22txm;在20~ 50mm内,采用线性拟合,得到U=0.026txm+9.4× l0 ;在50—70mm内,采用线性拟合,得到U=2× l0~L一0.501xm;在70~90mm内,采用线性拟合, 得到U=3.5×10~L一1.561.Lm。 6)若预评估结果随评估点变化呈不规律曲线 变化,如图3,且无法采用分段线性拟合。这时可考 虑采用第二种表述方式,即用范围表述。此时应对 各等级被校对象、所有测量点均开展预评估,以形成 完整的预评估报告。并从中挑选出最小值作为该测 量点的CMC。用范围表示时,应注意①测量范围分 段应尽量确保在同一计量单位。如表1中测量范围 一58— 围应与测量范围前后对应。 CMC汇总图 ” … 一 冒 厂、0. ” 一i 一九 / 器 厂 匠删 u / f 、 … 嚣 。 \./… 0 10 20 30 40 50 60 70 80 9O l0O 被测量(mm1 图3预评估结果随评估点变化趋势图 l 1 O CMC汇总罔 O 0 O 8 6 4 2 l 8 6 4 2 0 y=0.035x 5667厂… 善 i l 、 、一 …“ 越 器 y=0.02x一0. .9 咖 . 燕 y=0.0094x+0 0三 .5 A . 0.22--022 0 l0 20 3O 40 50 60 70 80 90 1OO 被测量(mm1 图4预评估结果随评估点变化趋势图 5结论和展望 CMC的评估和表示方式通常根据各专业实际 而定,并无一成不变的方法。本文的所提供的评定 方法其目的在于,通过对CMC评定流程的探索,减 少评估点,简化不确定度流程、降低运算量。以减少 实验室工作量,加快整改进度。 参考文献 [1]校准和测量能力(CMC)的表示方式应用指南[S].北 京:中国合格评定国家认可委员会,2011 [2] CNAS--GL02(能力验证结果的统计处理和能力评价 指南》[s].北京:中国合格评定国家认可委员 会,2006