测量不确定度评定举例
- 格式:docx
- 大小:146.43 KB
- 文档页数:21
相对频率偏差的测量不确定度评定1. 测量方法相对频率偏差:参考频标:铯原子频率标准5071A 被检频标:铷原子频率标准 频标比对器:PO7D 2. 测量结果测量10次,数据如下:oox f f f y -=)(τ3. 测量不确定度来源(1)铯原子频标不准引入的不确定度1u铯原子频标检定证书给出其频率准确度为5×10-13, 按B 类方法进行不确定度评定。
视其为均匀分布,包含因子3=k ,则有:13131109.23/105--⨯=⨯=u(2)铯原子频标不稳引入的不确定度2u测量相对频率偏差的取样时间为100s 。
铯原子频标检定证书给出其100s 频率稳定度为4.9×10-13,按A 类方法进行评定,k=1,则有:132109.4-⨯=u(3)频标比对器引入的不确定度3u频标比对器检定证书100s 比对不确定度为1.2×10-13,按A 类方法进行不确定度评定,k=1,则有:133102.1-⨯=u(4)测量重复性引入的不确定度4u实验标准偏差)(x s n1212109.11)()(-=⨯=--=∑n y yx s ni i in对于平均值,重复性测量引入的不确定度为:13124100.610/109.1--⨯=⨯=u3. 合成标准不确定度c u相对频率偏差测量结果的不确定度分量如下表:以上各不确定度分量互相独立各不相关,可得合成标准不确定度c u :21321321321324232221)100.6()102.1()109.4()109.2(----⨯+⨯+⨯+⨯=++++=u u u u u c 13104.8-⨯= 4. 扩展不确定度取k=2, 则扩展不确定度: 12102-⨯=U 5. 结论相对频率偏差:11100.7-⨯ 不确定度: 12102-⨯ (k=2)频率稳定度的测量不确定度评定1. 测量方法参考频标:高稳晶振8607 被检频标:铷原子频率标准 频标比对器:5120A 2. 测量结果3.不确定度来源(1) 参考频标引入的不确定度测量频率稳定度时使用的参考源为高稳晶振8607,根据其检定证书,其1 s 频率稳定度为7.2E-14,按B 类方法进行评定,k=1,则有:141102.7-⨯=u(2) 测量装置引入的不确定度测量装置使用5120,实测1 s 比对不确定度为1.19E-13,按A 类方法进行不确定度评定,k=1,则有:1321019.1-⨯=u(3) 有限次测量引入的不确定度按A 类方法进行有限次测量不确定度的评定。
测量不确定度评定的方法以及实例1.标准不确定度方法:U =sqrt(∑(xi-x̅)^2/(n-1))其中,xi表示测量值,x̅表示测量值的平均值,n表示测量次数。
标准不确定度包含随机误差和系统误差等。
例如,对一组长度进行测量,测得的数据为10.2、10.3、10.1、10.2、10.3,计算平均值为10.22,标准差为0.069、则标准不确定度为0.069/√5≈0.031,即U=0.0312.扩展不确定度方法:扩展不确定度是在标准不确定度的基础上,考虑到误差的正态分布,对标准不确定度进行扩展得到的结果,通常以U'表示。
其计算公式如下:U'=kU其中,k表示不确定度的覆盖因子,代表了误差分布的概率密度曲线下的面积,一般取k=2例如,对上述例子中的长度进行测量,标准不确定度为0.031,取k=2,则扩展不确定度为0.031×2=0.062,即U'=0.0623.组合不确定度方法:4.直接测量法:直接测量法是通过多次测量同一物理量,统计测得值的离散程度来评估测量的不确定度。
该方法适用于一些简单的测量,如长度、质量等物理量的测量。
例如,对一些小球的直径进行测量,测得的数据为2.51 cm、2.49 cm、2.52 cm、2.50 cm,计算平均值为2.505 cm,标准差为0.013 cm。
则标准不确定度为0.013/√4≈0.007 cm,即U=0.0075.间接测量法:间接测量法是通过已知物理量之间的数学关系,求解未知物理量的方法来评估测量的不确定度。
该方法适用于一些复杂的测量,如测量速度、加速度等物理量的测量。
例如,测量物体的速度v,则有v=S/t,其中S为位移,t为时间。
若S的不确定度为U_S,t的不确定度为U_t,则根据误差传递法则,计算得到v的不确定度为U_v = sqrt(U_S^2 + (U_t * (∂v/∂t))^2 )。
总之,测量不确定度评定的方法包括标准不确定度方法、扩展不确定度方法、组合不确定度方法、直接测量法和间接测量法。
标准不确定度B类评定的举例:(例1)校准证书上给出标称值为1000g的不锈钢标准砝码质量m s的校准值为,且校准不确定度为24g(按三倍标准偏差计),求砝码的标准不确定度。
评定:a =U =24g k=3则砝码的标准不确定度为u B(m s)= 24g/3 =8g(例2)校准证书上说明标称值为10的标准电阻,在23℃时的校准值为,扩展不确定度为90,置信水平为99%,求电阻的相对标准不确定度。
评定:由校准证书的信息知道:a =U99=90,P =;p241假设为正态分布,查表得到k=;则电阻校准值的标准不确定度为:u B(R S)=90/=35相对标准不确定度为:u B(R S)/ R S=×10-6。
(例3)手册给出了纯铜在20℃时线热膨胀系数20(Cu)为×10-6℃-1,并说明此值的误差不超过×10-6℃-1,求20(Cu)的标准不确定度。
评定:根据手册,a =×10-6℃-1,依据经验假设为等概率地落在区间内,即均匀分布,查表得,铜的线热膨胀系k3数的标准不确定度为:u (20)=×10-6℃-1/ =×10-6℃-1(例4) 由数字电压表的仪器说明书得知,该电压表的最大允许误差为(14×10-6×读数+2×10-6×量程),在10 V 量程上测1 V 时,测量10次,其平均值作为测量结果, V = V ,求电压表仪器的标准不确定度。
评定:电压表最大允许误差的模为区间的半宽度: a =(14×10-6× +2×10-6×10 V )=33×10-6 V=33 V 。
设在区间内为均匀分布,查表得到 。
则:电压表仪器的标准不确定度为: u (V )= 33 V/3=19 V[案例]:某法定计量技术机构为要评定被测量Y 的测量结果y 的合成标准不确定度u c (y )时,y 的输入量中,有碳元素C 的原子量,通过资料查出C 的原子量Ar (C )为:Ar (C )=±。
第一节有关术语的定义3.量值 value of a quantity一般由一个数乘以丈量单位所表示的特定量的大小。
例: 5.34m 或 534cm, 15kg, 10s,- 40℃。
注:对于不可以由一个乘以丈量单位所表示的量,能够参照商定参照标尺,或参照丈量程序,或二者参照的方式表示。
4.〔量的〕真值 rtue value〔of a quantity〕与给定的特定量定义一致的值。
注:(1)量的真值只有经过完美的丈量才有可能获取。
(2)真值按其天性是不确立的。
(3)与给定的特定量定义一致的值不必定只有一个。
5.〔量的〕商定真值 conventional true value〔of a quantity〕对于给定目的拥有适合不确立度的、给予特定量的值,有时该值是商定采纳的。
例: a) 在给定地址,取由参照标准复现而给予该量的值人作为给定真值。
b) 常数委员会 (CODATA)1986年介绍的阿伏加得罗常数值 6.0221367 × 1023mol-1。
注:(1)商定真值有时称为指定值、最正确预计值、商定值或参照值。
(2)经常用某量的多次丈量结果来确立商定真值。
13.影响量 influence quantity不是被丈量但对丈量结果有影响的量。
例: a) 用来丈量长度的千分尺的温度;b)沟通电位差幅值丈量中的频次;c)丈量人体血液样品血红蛋浓度时的胆红素的浓度。
14.丈量结果 result of a measurement由丈量所获取的给予被丈量的值。
注:(1)在给出丈量结果时,应说明它是示值、示修正丈量结果或已修正丈量结果,还应表示它能否为几个值的均匀。
(2)在丈量结果的完好表述中应包含丈量不确立度,必需时还应说明有关影响量的取值范围。
15.〔丈量仪器的〕示值 indication〔of a measuring instrument〕丈量仪器所给出的量的值。
注:(1)由显示器读出的值可称为直接示值,将它乘以仪器常数即为示值。
标准不确定度A类评定的实例【案例】对一等活塞压力计的活塞有效面积检定中,在各种压力下,测得10次活塞有效面积与标准活塞面积之比l (由l的测量结果乘标准活塞面积就得到被检活塞的有效面积)如下:0。
250670 0。
250673 0.250670 0。
250671 0.250675 0。
250671 0。
250675 0.250670 0。
250673 0。
250670问l 的测量结果及其A 类标准不确定度。
【案例分析】由于n =10, l 的测量结果为l ,计算如下∑===n i i .l n l 125067201 由贝塞尔公式求单次测量值的实验标准差()612100521-=⨯=--=∑.n l l )l (s n i i由于测量结果以10次测量值的平均值给出,由测量重复性导致的测量结果l 的A 类标准不确定度为610630-=⨯=.)l (u n )l (s A 【案例】对某一几何量进行连续4次测量,得到测量值:0。
250mm 0.236mm 0.213mm 0。
220mm ,求单次测量值的实验标准差。
【案例分析】由于测量次数较少,用极差法求实验标准差.)()(i i x u CR x s ==式中,R-—重复测量中最大值与最小值之差;极差系数c及自由度ν可查表3-2表3-2极差系数c及自由度ν查表得c n =2.06mm ../mm )..()x (u CR )x (s i i 018006221302500=-=== 2)测量过程的A 类标准不确定度评定对一个测量过程或计量标准,如果采用核查标准进行长期核查,使测量过程处于统计控制状态,则该测量过程的实验标准偏差为合并样本标准偏差S P 。
若每次核查时测量次数n 相同,每次核查时的样本标准偏差为Si ,共核查k 次,则合并样本标准偏差S P 为k s s ki ip ∑==12此时S P 的自由度ν=(n —1)k .则在此测量过程中,测量结果的A 类标准不确定度为 n S A P u '=式中的n '为本次获得测量结果时的测量次数。
不确定度评估实例1、测量问题本次评定实验以物资(商品)检验所游标卡尺09059为测试量具,用游标卡尺测量结构长度270mm的长度ι。
已知卡尺的最大误差为1mm。
用6次测量的平均值作为测量结果。
卡尺的温度效应、弹性效应及其他不确定度来源均忽略不计。
2、数学模型卡尺上得到的读数χ即为测量结果,故得被测长度ι=χ。
但除了读数χ可能引入测量不确定度外,卡尺刻度误差对测量结果也会有影响。
由于卡尺的校准证书未给出其示值误差,因此只能根据其最大允许误差来估计它对测量结果的影响。
若卡尺刻度误差对测量结果的影响διS,则数学模型可以表示为ι=χ+διS式中διS的数学期望值为零,即Ε(διS)=0,但需考虑其不确定度,即μ(διS)≠0。
数学模型是相对的,即使对于同样的被测量,当要求的测量准确度不同时,需要考虑的测量不确定度来源也会有相应的增减,因此数学模型也会不同。
3、测量不确定度分量本测量共有两个不确定度分量,由读数的重复性引入的不确定度μ(χ)和卡尺刻度误差所引起的不确定度μ(διS)。
⑴读数χ的不确定度,μ1(ι)=μ(χ)6次测量结果分别为270、3mm270、1mm270mm271、4mm269、8mm271、2mm则6次测量结果的平均值为==270、47mm平均值的实验标准差为 s()==0、074mm故μ1(ι)=μ()=s()=0、074mm⑵卡尺误差引入的不确定度, μ2(ι)=μ(διS)由于证书未给出卡尺的示值误差,故卡尺刻度误差引入的不确定度由卡尺的最大允许误差得到。
已知卡尺的最大误差为1mm,并以矩形分布估计,于是μ2(ι)=μ(διS)==0、577mm下表给出不确定度分量汇总表符号栏中u1=s1 意为用实验标准s来表示不确定度,言外之意是该不确定度分量有A类评定得到的。
反之,对于未标u=s的不确定度分量,则表示是由B 类评定得到的。
这是经常采用的标明A类评定和B类评定不确定度分量的方法之一。
第一节有关术语的定义3.量值value of a quantity 一般由一个数乘以测量单位所表示的特定量的大小。
例:5.34m或534cn, 15kg, 10s,—40C。
注:对于不能由一个乘以测量单位所表示的量,可以参照约定参考标尺,或参照测量程序,或两者参照的方式表示。
4.〔量的〕真值rtue value 〔of a quantity 〕与给定的特定量定义一致的值。
注:(1)量的真值只有通过完善的测量才有可能获得。
(2)真值按其本性是不确定的。
(3)与给定的特定量定义一致的值不一定只有一个5.〔量的〕约定真值conventional true valueof a quantity 〕对于给定目的具有适当不确定度的、赋予特定量的值,有时该值是约定采用的。
例:a) 在给定地点,取由参考标准复现而赋予该量的值人作为给定真值b)常数委员会(CODATA)1986年推荐的xx加得xx数值 6.0221367 x 1023mol-1。
注:(1)约定真值有时称为指定值、最佳估计值、约定值或参考值。
(2)常常用某量的多次测量结果来确定约定真值。
1 3.影响量influence quantity 不是被测量但对测量结果有影响的量。
例:a)用来测量xx的千分尺的温度;b) 交流电位差幅值测量中的频率;c) 测量人体血液样品血红蛋浓度时的胆红素的浓度。
14.测量结果result of a measurement由测量所得到的赋予被测量的值。
注:(1)在给出测量结果时,应说明它是示值、示xx 测量结果或已xx 测量结果,还应表明它是否为几个值的平均(2)在测量结果的完整表述中应包括测量不确定度,必要时还应说明有关影响量的取值范围。
15.〔测量仪器的〕示值indication 〔of a measuring instrument 〕测量仪器所给出的量的值注:(1)由显示器读出的值可称为直接示值,将它乘以仪器常数即为示值。
一、力学测量应用实例用拉力试验机测量金属试件拉伸强度。
已知试件的标准直径mm d 10=,断裂时拉力为40kN 。
拉力试验机的量程为200kN ,分度值为0.5kN ,示值误差为F %1+,示值误差的不确定度为0.2%F 。
试件直径用千分尺测量,其示值误差为m μ3+。
求拉伸强度的测量不确定度。
2.1 数学模型 24d FA F R m π==m R — 拉伸强度 (Mpa )A — 试件截面积 (2mm )d — 试件直径 (mm )F — 拉力 (N )2.2 不确定度传播律)(4)()(222d u F u R u rel rel m rel c +=2.3 求相对标准不确定度分量)(d u rel2.3.1 千分尺示值误差导致的不确定度 )(1d u以均匀分布估计 m d u μ73.133)(1==2.3.2 由操作者引起的测量不确定度)(2d u经验估计,该测量误差在m μ10+范围内,以均匀分布估计, m d u μ77.5310)(2==以上二者合成 m d u μ02.677.573.1)(22=+=以上相对不确定度表示: %06.01010*02.6)(3==-d u rel2.4 求拉力F 的测量不确定度 )(F u rel2.4.1 拉力机的示值误差引入的测量不确定度)(1F u由于仪器说明书未说明置信概率,故取2=k%5.0%1)(1==k F u2.4.2 拉力机校准的不确定度)(2F u这是由上一级标准器对拉力机校准时产生的不确定度,即拉力机示值误差的不确定度,校准证书亦未给出置信概率,故取2=k%1.0%2.0)(2==k F u2.4.3 拉力机读数不准产生的不确定度)(3F u人工读数可以估计到刻度的五分之一,即0.1kN ,读数误差的不确定度可按均匀分布估计,3=k %144.03401.0)(3==F u以上三者合成 %53.0)144.0(%)1.0(%)5.0()(222=++=F u rel2.5 合成标准不确定度c u %543.0%)06.0(4%)53.0()(4)()(2222=+=+=d u F u R u rel rel m rel c 223.5094mm N d F R m ==π 28.2%543.0*3.509)(mmN R u R u m rel c m c === 2.6 扩展不确定度 U取包含因子 2=k26.58.2*2mm N ku U c ===2.7 测量结果报告 2)6.53.509(mm N R m +=……二、 电学测量应用实例用数学电压表测量电压9次,得到平均值V v 928571.0=,标准偏差V v s μ36)(=。
丈量不确立度评定实例一.体积丈量不确立度计算1.丈量方法直接丈量圆柱体的直径 D 和高度 h,由函数关系是计算出圆柱体的体积v D 2 4由分度值为 0.01mm 的测微仪重复 6 次丈量直径 D 和高度 h,测得数据见下表。
表:丈量数据i123456D i / mm10.07510.08510.09510.06510.08510.080 h i / mm10.10510.11510.11510.11010.11010.115计算: D 10.080 mm, h 10.110 mmV D2 h 806.8 mm3 42.不确立度评定剖析丈量方法可知,体积 V 的丈量不确立度影响要素主要有直径和高度的重复丈量惹起的不确立度 u1, u2和测微仪示值偏差惹起的不确立度 u3。
剖析其特色,可知不确立度 u1,u2应采纳A类评定方法,而不确立度 u3采纳B类评定方法。
①.直径 D 的重复性丈量惹起的不确立度重量直径 D 的 6 次丈量均匀值的标准差:s D0.0048 mm直径 D 偏差传达系数:V D hD2直径 D 的重复性丈量惹起的不确立度重量:u1Vs D 0.77mm3 D② .高度 h 的重复性丈量惹起的不确立度重量高度 h 的 6 次丈量均匀值的标准差:s h0.0026 mm高度 h 的偏差传达系数:V D 2h4高度 h 的重复性丈量惹起的不确立度重量:u2Vs h 0.21mm3 h③测微仪示值偏差惹起的不确立度重量由说明书获取测微仪的示值偏差范围0.005mm ,按均匀散布,示值的标准不确立度0.005u q0.00293由示值偏差惹起的直径丈量的不确立度u3D V u qD由示值偏差惹起的高度丈量的不确立度u 3hVhu q由示值偏差惹起的体积丈量的不确立度重量221.04 mm 3u 3u3 Du3h3. 合成不确立度评定u c u 12u 22u 321.3 mm 34. 扩展不确立度评定当置信因子 k 3时,体积丈量的扩展不确立度为Uku c 3 1.3 3.9 mm 35.体积丈量结果报告V V U806.8 3.9 mm 3考虑到有效数字的观点,体积丈量的结果应为V807 4 mm 3二.伏安法电阻丈量不确立度计算1.丈量方法:经过丈量电阻两头电压和所经过的电流,计算被测电阻。
欢迎阅读测量不确定度评定举例A.3.1 量块的校准通过这个例子说明如何建立数学模型及进行不确定度的评定;并通过此例说明如何将相关的输入量经过适当处理后使输入量间不相关,这样简化了合成标准不确定度的计算。
最后说明对于非线性测量式中:L—被校量块长度;L s—标准量块在20℃时的长度,由标准量块的校准证书给出;?—被校量块的热膨胀系数;?s—标准量块的热膨胀系数;?—被校量块的温度与20℃参考温度的差值;?s —标准量块的温度与20℃参考温度的差值。
在上述测量模型中,由于被校量块与标准量块处于同一温度环境中,所以?与?s 是相关的量;两个量块采用同样的材料,?与?s 也是相关的量。
为避免相关,设被校量块与标准量块的温度差为??,??= ?-?s ;他们的热膨胀系数差为??,??= ?-?s ;将?s = ?-?? 和 ?=??+?s 代入式(A.1),由此,数学模型可改写成:34s 量结果的不确定度没有影响。
合成标准不确定度公式可写成(A.5):)()()()()(22222222θαδαδθu l u l d u l u l u s s s s c +++= (A.5)4).标准不确定度分量的评定○1标准量块的校准引入的标准不确定度u (l s )标准量块的校准证书给出:校准值为l s =50.000623mm ,U = 0.075?m (k =3),有效自由度为?eff (l s )=18。
则标准量块校准引入的标准不确定度为:u (L s )=0.075/3=25nm , ?eff (L s )=18 ○2测得的长度差引入的不确定度u (d ) a. 用对两个量块的长度差进行25次独立重复观测,用贝塞尔公c. 由以上分析得到长度差引入的标准不确定度分量u (d )为: 8.97.85.4)()()(2222=+=+=d u d u d u nm 自由度?eff (d )为:○3膨胀系数差值引入的标准不确定度u (??)估计两个量块的膨胀系数之差在?1×10-6℃-1区间内,假设在区间内为均匀分布,则标准不确定度为: u (??)=1×10-6℃-1/3=0.58×10-6℃-1自由度:估计u (??)的不可靠程度⎦⎤⎢⎣⎡∆)()(ααδδu u 为10%,计算得到?(??)=50%)10(12=- u(△)= 0.5℃/2 =0.35℃ θ的标准不确定度可由下式得到:u(θ)= 41.035.02.0)()(2222=+=∆+u u θ℃ 由于c 4 = c θ=0=-=∂∂θδθs l f, 这个不确定度对l 的不确定度不引入一阶的贡献, 然而它具有二阶贡献.○6 热膨胀系数引入的标准不确定度u (αS ) 标准量块的热膨胀系数给定为αS =11.5×10-6℃-1, 具有一个矩形分布的不确定度,其界限为?2×10-6℃-1, 则标准不确定度为: u (αS )= 2×10-6℃-1/3 = 1.2×10-6℃-1 由于c 3 = c αs =0=-=∂∂θδαs Sl f, 这个不确定度对L 的不确定度不?eff (l )=3.172)6.16(50)9.2(12)8.9(18)25()32(44444=+++ 取?eff (l )=17 6)确定扩展不确定度要求包含概率P 为0.99,由?eff (l )=17,查表得:t0.99(17)=2.90,取k99= t0.99(17)=2.90,扩展不确定度U99= k99u c(l)= 2.90,×32nm=93nm。
气相色谱法测定绝缘油溶解气体含量测量不确定度的评定(供参考)一、概述1.1 目的评定绝缘油溶解气体含量测量结果的不确定度。
1.2 依据的技术标准GB/T 17623-1998《绝缘油中溶解气体组分含量的气相色谱测定法》。
1.3 使用的仪器设备(1) 气相色谱分析仪HP5890,经检定合格。
(2) 多功能全自动振荡仪ZHQ701,经检定合格,允差±1℃,分辨力0.1℃。
(3) 经检验合格注射器,在20℃时,体积100mL±0.5mL;体积5mL±0.05mL;体积1mL±0.02mL。
1.4 测量原理气相色谱分析原理是利用样品中各组分,在色谱柱中的气相和固定相之间的分配及吸附系数不同,由载气把绝缘油中溶解气体一氧化碳、二氧化碳、甲烷、乙烷、乙烯、乙炔、氢气带入色谱柱中进行分离,并经过电导和氢火焰检测器进行检测,采用外标法进行定性、定量分析。
1.5 测量程序(1) 校准。
采用国家计量部门授权单位配制的甲烷标准气体。
进样器为1mL玻璃注射器,采用外标气体的绝对校正因子定性分析。
(2) 油样处理。
用100mL玻璃注射器A,取40mL油样并用胶帽密封,并用5mL玻璃注射器向A中注入5mL氮气。
将注入氮气的注射器A放入振荡器中振荡脱气,在50℃下,连续振荡20分钟,静止10分钟。
(3) 油样测试。
然后用5mL玻璃注射器将振荡脱出的气体样品取出,在相同的色谱条件下,进样量与标准甲烷气体相同,对样品进行测定,仪器显示谱图及测量结果。
气体含量测定过程如下。
1.6 不确定度评定结果的应用符合上述条件或十分接近上述条件的同类测量结果,一般可以直接使用本不确定度评定测量结果。
二、 数学模型和不确定度传播律2.1 根据GB/T 17623-1998《绝缘油中溶解气体组分含量的气相色谱测定法》试验方法,绝缘油中溶解气体含量C 的表示式为S s=⨯hC C h μL/L (1) 式中,C ——被测绝缘油中溶解气体甲烷含量,μL/L ;C S ——标准气体中甲烷含量,μL/L ; h ——被测气体中甲烷的峰高A ; h s ——标准气体中甲烷的峰高A 。
2021 June第测量结果的不确定度评定实例分析刘海利中国石化销售股份有限公司油品技术研究所以GB/T 261—2008《闪点的测定 宾斯基-马丁闭口杯法》测量车用柴油闭口闪点为例,按照JJF 1059.1—2012《测量不确定度与表示》要求进行检测实验室测量不确定度评定,通过对实验室测量结果的不确定度评定,实现测量结果不确定度规范与正确表达,进而提升实验室测量结果质量。
作者简介:刘海利,硕士,高级工程师,现主要从事油品质量管理与应用研究工作。
E-mail:liuhaili119@163.com测量不确定度是表征检测和校准实验室测量结果的质量参数,对于一定的测量结果而言,它的不确定度值越小,其质量就越高,使用价值也越高;反之则低。
在CNAS-CL01:2018《检测和校准实验室能力认可准则》中,要求实验室应制定与检测工作相适应的测量不确定度评定程序,对每一项有数值要求的结果进行测量不确定度评定。
因此,测量不确定度评定在检测和校准实验室认可中是一项不可缺少的重要工作[1]。
JJF 1059.1—2012《测量不确定度评定与表示》是评定不确定度最常用、最基本的方法[2]。
闭口闪点是轻质油品运输、储存和使用安全的重要指标,本文以GB/T 261—2008《闪点的测定 宾斯基-马丁闭口杯法》测量车用柴油闭口闪点不确定度为例,阐述测量闭口闪点不确定度步骤,为实验室开展所有测量项目结果的不确定度评定提供参考,提高实验室检测能力。
Teat and Appraisal测试与评定8282三期83一2021 June第各不确定度分量的评定重复性测量引入的标准不确定度分量u 1(T c )车用柴油闭口闪点测量时,试样量、加热速率、搅拌速率、试验过程中温度计深入位置、温度计读数、压力表读数等随机因素带来的不确定度,一并列入重复性测量不确定度分量中进行评定。
试验用温度计修正值∆T =0.0 ℃,压力表修正值∆p =0.1 kPa,在重复性试验条件下,对同一试样独立重复测量10次,结果见表1。
测量不确定度评定举例A.3.1 量块的校准通过这个例子说明如何建立数学模型及进行不确定度的评定;并通过此例说明如何将相关的输入量经过适当处理后使输入量间不相关,这样简化了合成标准不确定度的计算。
最后说明对于非线性测量函数考虑高阶项后测量不确定度的评定结果。
1).校准方法标称值为50mm 的被校量块,通过与相同长度的标准量块比较,由比较仪上读出两个量块的长度差d ,被校量块长度的校准值L 为标准量块长度L s 与长度差d 之和。
即:L=L s +d实测时,d 取5次读数的平均值d ,d =0.000215mm ,标准量块长度L s 由校准证书给出,其校准值L s =50.000623mm 。
2)测量模型长度差d 在考虑到影响量后为:d =L (1+?? )-L s (1+?s ?s ) 所以被校量的测量模型为:此模型为非线性函数,可将此式按泰勒级数展开:L =ΛΛ+-++)(θαθαs s s s L d L忽略高次项后得到近似的线性函数式: )(θαθα-++=s s s s L d L L()式中:L —被校量块长度;L s —标准量块在20℃时的长度,由标准量块的校准证书给出; ? —被校量块的热膨胀系数; ?s —标准量块的热膨胀系数;? —被校量块的温度与20℃参考温度的差值; ?s —标准量块的温度与20℃参考温度的差值。
在上述测量模型中,由于被校量块与标准量块处于同一温度环境中,所以?与?s 是相关的量;两个量块采用同样的材料,?与?s 也是相关的量。
为避免相关,设被校量块与标准量块的温度差为??,??= ?-?s ;他们的热膨胀系数差为??,??= ?-?s ;将?s = ?-?? 和 ?=??+?s 代入式(),由此,数学模型可改写成:= ][θαδαθδs s s l d l +-+ () 测量模型中输入量??与?s 以及??与?不相关了。
特别要注意:在此式中的??和??是近似为零的,但他们的不确定度不为零,在不确定度评定中要考虑。
由于??和??是近似为零,所以被测量的估计值可以由下式得到:L =L s +d () 3).测量不确定度分析 根据测量模型,即: l = ][θαδαθδs s s l d l +-+由于各输入量间不相关,所以合成标准不确定度的计算公式为:)()()()()()()(222222222222θδαδθαδδθαθαu c u c u c u c d u c l u c l u s d s s c s +++++= () 式中灵敏系数为:1)(11=+-=∂∂==θαδαθδs ss l fc c , 由此可见,灵敏系数c 3和c 4为零,也就是说明?s 及? 的不确定度对测量结果的不确定度没有影响。
合成标准不确定度公式可写成: )()()()()(22222222θαδαδθu l u ld u l u l u s s s s c +++=()4).标准不确定度分量的评定○1标准量块的校准引入的标准不确定度u (l s ) 标准量块的校准证书给出:校准值为l s =50.000623mm ,U = 0.075?m (k =3),有效自由度为?eff (l s )=18。
则标准量块校准引入的标准不确定度为: u (L s )=3=25nm , ?eff (L s )=18 ○2测得的长度差引入的不确定度u (d ) a. 用对两个量块的长度差进行25次独立重复观测,用贝塞尔公式计算的实验标准偏差为s (d )=13nm ;本次比较时仅测5次,取5次测量的算术平均值为被校量块的长度,所以读数观测的重复性引入的标准不确定度u (d )是平均值的实验标准偏差为s (d )8.55/13/)()()(====n d s d s d u nm由于s (d )是通过25次测量得到,所以u (d )的自由度?1=25-1=24。
b. 由比较仪示值不准引起长度差测量的不确定度u B (d):由比较仪的校准证书给出最大允许误差为±0.015?m,有效期内的检定证书证明该比较仪的示值误差合格.则由比较仪示值不准引起长度差测量的标准不确定度用B 类评定,可能值区间的半宽度a 为0.015?m ,设在区间内呈均匀分布, 取包含因子k = 3。
标准不确定度u B (d)为: u (d )=0.015?m /3= 8,7 nm按下式估计其自由度: 2)()(21-⎥⎦⎤⎢⎣⎡∆≈i i i x u x u ν假设评定u B (d)的不可靠程度达25%, 计算得到8≈i ν c. 由以上分析得到长度差引入的标准不确定度分量u (d )为: 8.97.85.4)()()(2222=+=+=d u d u d u nm 自由度?eff (d )为:○3膨胀系数差值引入的标准不确定度u (??) 估计两个量块的膨胀系数之差在?1×10-6℃-1区间内,假设在区间内为均匀分布,则标准不确定度为: u (??)=1×10-6℃-1/3=×10-6℃-1自由度:估计u (??)的不可靠程度⎥⎦⎤⎢⎣⎡∆)()(ααδδu u 为10%,计算得到?(??)=50%)10(212=- ○4量块温度差引入的标准不确定度u (??) 希望被校量块与标准量块处于同一温度,但实际存在温度差异,温度差估计以等概率落在?0.05℃区间内,则标准不确定度为: u (??)=3=0.029℃估计u (??)只有50%的可靠性,计算得到自由度为:?(??)=2%)50(212=-○5量块温度偏差引入的标准不确定度u (θ) 报告给出的测试台温度为?℃, 在热作用下温度的近似周期性变化的幅度为0.5℃. 平均温度的偏差值为:1.0209.19-=-=θ(℃)由于测试台的平均温度的不确定度引起的θ 的标准不确定度为: u (θ)=0.2 ℃而温度随时间周期变化形成U 形的分布(即反正弦分布),则: u(△)= 0.5℃/2 =0.35℃θ的标准不确定度可由下式得到:u(θ)= 41.035.02.0)()(2222=+=∆+u u θ℃ 由于c 4 = c θ=0=-=∂∂θδθs l f, 这个不确定度对l 的不确定度不引入一阶的贡献, 然而它具有二阶贡献.○6 热膨胀系数引入的标准不确定度u (αS ) 标准量块的热膨胀系数给定为αS =×10-6℃-1, 具有一个矩形分布的不确定度,其界限为?2×10-6℃-1, 则标准不确定度为: u (αS )= 2×10-6℃-1/3 = ×10-6℃-1 由于c 3 = c αs =0=-=∂∂θδαs Sl f, 这个不确定度对L 的不确定度不引入一阶的贡献, 然而它具有二阶贡献. 5)计算合成标准不确定度的 ○1计算灵敏系数 由标准量块的校准证书得到L s =50.000623mm ,被校量块与参考温度20℃之差估计为-0. 1℃,标准量块的热膨胀系数?s 为×10-6℃-1,由这些信息计算得到: c 1=1,c 2=1, c 3=0, c 4=0,c 5=-l s ? = -50.000623mm×(-0.1℃)=-5.0000623mm℃, c 6=-l s ?s = -50.000623mm××10-6℃-1=×10-4mm ℃-1○2计算合成标准不确定度 = 32 nm ○3u (l )的自由度: ?eff (l )=3.172)6.16(50)9.2(12)8.9(18)25()32(44444=+++ 取?eff (l )=17 6)确定扩展不确定度要求包含概率P 为,由?eff (l )=17,查表得: (17)=,取k 99= (17)=,扩展不确定度U 99= k 99u c (l )= ,×32nm=93nm 。
7)校准结果:l =l s +d =50.000623mm+ 0.000215mm =50.000838mm U 99= 93nm (?eff =17) 或l =(?)mm其中?号后的值是扩展不确定度U 99,由u c =32nm 乘包含因子k =得到,k 是由自由度?=17,包含概率p =时查t 分布值表得到,由该扩展不确定度所包含的区间具有包含概率为。
量块校准时标准不确定度分量汇总见表表 量块校准时标准不确定度分量汇总表可见,不确定度的主要分量显然是标准量块的不确定度u (l s )= 25nm 。
注: 用蒙特卡洛法(MCM)验证, 得到:传播输出量分布的标准偏差 u(l)=36nm, 最小包含区间的半宽度U 99=94nm, 与本规范的结果基本一致.说明本规范的方法评定不确定度基本可信的. 8)考虑二阶项时不确定度的评定前面所进行的不确定度的评定是不完全的,实际上在本案例中,测量模型存在着明显的非线性,在泰勒级数展开中的高阶项不可忽略。
在合成标准不确定度评定中,有两项明显的不可忽略的二阶项对u c (l )有贡献:)()(222θδαu u l s =(0.05m )× ×10-6℃-1) ×(0.41℃)= )()(222θδαu u l s s =(0.05m )× ×10-6℃-1) ×(0.029℃)=考虑二阶项后的合成标准不确定度:u′c(l)=227.12+=34nm32+7.11扩展不确定度U99(l)= 99nm (k =,?eff =16,P =或相对扩展不确定度U99/l=×10-6A.3.2温度计的校准这个例子说明用最小二乘法获得线性校准曲线时,如何用校准曲线的截距、斜率和他们的估计方差与协方差,计算由校准曲线获得的预期修正值及其标准不确定度。
测量问题温度计是用与已知的参考温度相比较的方法校准的。
相应的已知参考温度为t R,k,其温度范围为21℃到27℃。
进行了n=11次比较,温度计的温度读数为t k,温度计读数的修正值为b k=t R,k-t k。
根据测得的修正值b k和测得的温度t k,用最小二乘法拟合成直线得到温度计修正值的线性校准曲线b(t)为:b(t)=y1+y2(t-t0)式中, y1为校准曲线的截距,y2为校准曲线的斜率,t0是所选择的参考温度;y1和y2是两个待测定的输出量。
一旦找到y1和y2以及它们的方差和协方差,式可用于预示温度计对任意一个温度值t的修正值和最小二乘法拟合引入的标准不确定度。
by10k图A3.2.1最小二乘法拟合的示意图最小二乘法拟合根据最小二乘法和的假设条件,输出量y 1和y 2及它们的估计方差和协方差是在残差平方和Q 最小时得到:Q=21201[()]nk k k s b y y t t ==---∑这就导出了y 1和y 2的以下各公式,它们的实验方差s 2(y 1)和s 2(y 2),以及它们的估计的相关系数r (y 1,y 2)=s (y 1,y 2)/ s (y 1) s (y 2),其中s (y 1,y 2)是估计的协方差()()()()21kkk kkb b y Dθθθ-=∑∑∑∑ (A.7a)()()Db b n y k k k k ∑∑∑-=θθ22221()k s s y Dθ=∑ (A.7c)222()s s y n D=r (y 1,y 2)= θ 22[()]2kk b b t sn -=-∑ (A.7f)=∑∑-=-22)()(t t n n k k θθ (A.7g) 式中,k =1,2,…,n ; ?k =t k -t 0; n k /)(∑=θθ; n t t k /)(∑=;[b k -b (t k )]是在t k 温度时测得或观测到的修正值b k 与拟合曲线b (t )=y 1+y 2(t-t 0)上在t k 时预示的修正值b (t k )之间的差值;估计方差s 2是总的拟合的不确定度的度量;其中因子n -2反映了由n 次观测确定二个参数y 1和y 2时,s 2的自由度为?=n -2。