煤炭试验方法不确定度评定.docx
《煤炭试验方法不确定度评定.docx》由会员分享,可在线阅读,更多相关《煤炭试验方法不确定度评定.docx(28页珍藏版)》请在冰豆网上搜索。
煤炭试验方法不确定度评定
1范围
本标准规定了煤炭检验不确定度评定的术语和定义、不确定度的主要来源、不确定度评定的基本程序、标准不确定度的评定和报告与表示。
本标准适用于煤炭检验中不确定度的评定。
2规范性引用文件
下列文件中的条款通过本标准的引用而成为本标准的条款。
凡是注日期的引用文件,其随后所有的修改单(不包括勘误的内容)或修订版均不适用于本标准,然而,鼓励根据本标准达成协议的各方研究是否可使用这些文件的最新版本。
凡是不注日期的引用文件,其最新版本适用于本标准。
GB/T483煤炭分析试验方法一般规定
JJF1059测量不确定度评定与表示
3术语和定义
GB/T483-1998和JJF1001-1998规定的术语和定义及下列术语和定义适用于本标准。
3.1可测量的量measurablequality
现象、物体或物质可定性区别和定量确定的属性。
例如:
煤炭或煤灰中各种元素的含量;煤炭各种特征属性:
焦化特性指标和气化特性指标,等。
3.2量值valueofaquantity
一般由一个数乘以测量单位所表示的特定量的大小。
例如:
-50°C,10%(m/m),5.34mm,5kg,100μg/mL。
3.3真值truevalue
与给定的特定量定义一致的值。
3.4被测量measurand
作为测量对象的特定量。
3.5实验标准(偏)差experimentalstandarddeviation
对同一被测量作n次测量,表征结果分散性的量s可按下式(贝塞尔公式)算出:
式中:
x是第i次测量结果;
是n次测量的算术平均值。
3.6不确定度uncertainty
表征合理地赋予被测量之值的分散性,与测量结果相联系的参数。
注:
此参数可以是诸如标准差或其倍数,或说明了置信水准的半宽度。
3.7标准不确定度 standarduncertainty
以标准差表示的测量不确定度。
3.8相对标准不确定度relativestandarduncertainty
标准不确定度除以被测量之值或参考值所得的商。
3.9标准不确定度的A类评定 typeAevaluationofstandarduncertainty
用对系列测量数据进行统计分析的方法来评定标准不确定度。
注:
不确定度的A类评定有时也称为A类不确定度评定。
3.10标准不确定度的B类评定 typeBevaluationofstandarduncertainty
用不同于对系列测量数据进行统计分析的方法来评定标准不确定度。
注:
不确定度的B类评定有时也称为B类不确定度评定。
3.11合成标准不确定度 combinedstandarduncertainty
当测量结果是由若干个其他量的值求得时,按其他各量的方差和协方差算得的标准不确定度。
注:
它是测量结果标准差的估计值。
3.12扩展(或展伸)不确定度 expandeduncertainty
确定测量结果区间的量,合理赋予被测量之值分布的大部分可望含于此区间。
3.13包含因子 coveragefactor
为求得扩展不确定度,对合成标准不确定度所乘之数字因子。
3.14重复性限repeatabilitylimit
在同一试验室中,由同一操作者,用同一台仪器,对同一分析试验煤样,于短期内所作的重复测定,所得结果间的差值(在95%概率下)的临界值。
3.15再现性临界差reproducibilitycriticallimit
在不同试验室中,对从煤样缩制最后阶段的同一煤样中分取出来的、具有代表性的部分所作的重复测定,所得结果的平均值间的差值(在95%概率下)的临界值。
4不确定度的主要来源
煤炭检验中可能导致不确定度的来源如下:
a)被测量定义的不完整性;
b)样品的代表性:
采样、制样、样品的污染、样品的均匀性等引入的样品代表性不充分和不一致;
c)样品称量和溶液定容时的不确定度;
d)化学试剂、药品的纯度,以及基体效应引入的误差;
e)读取数据时的人为偏移;
f)测量仪器的计量性能的局限性,如灵敏度、分辨率、稳定性、噪声、测量精密度、计量检定或校准等引入的测量结果的变动;
g)测量标准或标准物质的不确定度;
h)校准曲线拟合度的不确定度;
i)引用的数据或其他参数的不确定度;
j)测量方法和测量程序中的局限性和假设;
k)在相同条件下被测量在重复测量中的变化;
l)对于那些尚未认识到的系统效应,显然是不可能在不确定度评定中予以考虑的,但它可能导致测量的不确定度。
上述不确定度的来源可能相关,例如,第k项可能与前面各项有关。
5测量不确定度评定程序
对煤炭分析试验进行测量不确定度评定,一般需经过如下步骤:
a)被测量说明;
b)不确定度来源分析;
c)标准不确定度的A类评定;
d)标准不确定度的B类评定;
e)合成标准不确定度的评定;
f)扩展不确定度的评定;
g)不确定度的报告与表示。
5.1被测量说明
在说明被测量时一般应包括以下内容:
a)试验过程描述;
b)被测量的影响参数;
c)被测量与其影响参数(输入量)的关系,如果可能,给出被测量y与影响参数xi的关系式:
y=f(x1,x2,x3……);
d)已知的系统性影响的校正值;
e)操作说明书等文件中给出的有关信息。
5.2不确定度来源分析
根据所测定项目具体情况,识别出不确定度源。
5.3不确定度的评定
不确定度评定包括标准不确定度的A类评定、标准不确定度的B类评定、合成标准不确定度的评定、扩展不确定度的评定。
注:
在不确定评定中,需剔除测量结果中的异常值,异常值的剔除应通过对数据的适当检验进行。
不确定度的A类评定
5.3.1.1一组数据单次测量值和算术平均值的标准不确定度
在重复性或再现性条件下对被测量x进行了n次独立测量,有x1,x2,……,xn,其算术平均值
、单次测量的标准偏差s(x)和平均值的标准偏差s(
)分别按式
(1)、式
(2)、式(3)计算:
…………………………………………………
(1)
……………………………………
(2)
……………………………(3)
式中:
xi——第i次测量结果;
n——测量次数。
单次测量的标准不确定度u(x)按式(4)计算:
u(x)=s(x)…………………………………………………(4)
平均值的标准不确定度u(
)按式式(5)计算:
u(
)=
……………………………………………(5)
在煤炭检验中,除特别要求者外,每项分析项目对同一煤样进行两次测定,此时,n=2。
5.3.1.2多组数据单次测量值和总算术平均值的标准不确定度
在重复性或再现性条件下对被测量x进行了m组测量,每组测量进行n次独立测量。
将第j组测量中,被测量记为xj,即有xj1,xi2,……,xjn,其平均值为
。
对m组这样的被测量,则被测量的总平均值
、单次测量的标准偏差s(x)和总平均值的标准偏差s(
)分别按式(6)、式(7)、式(8)计算:
…………………………………………(6)
………………………(7)
……………………………(8)
单次测量的标准不确定度u(x)和总平均值的标准不确定度u(
)分别为:
u(x)=s(x)
u(
)=
5.3.1.3校准曲线拟合参数引入的标准不确定度
当被测量x的估计值是由试验数据从最小二乘法拟合的校准曲线上得到时,校准曲线拟合参数引入的标准不确定度按下述步骤进行评定。
1)建立校准曲线
当用n个标准物质进行校准曲线的拟合并对每个标准物质进行m次测定时,假定对第i个(i=1~m)校准标准物质(其标准值为xcal)每次测量的响应值为Zcal。
以xcal为横坐标,以Zcal为纵坐标(用最小二乘法)拟合校准曲线如式(9):
Z=Ax+B…………………………(9)
式中:
Z-测量的响应值;
x-与向应值Z线性相关的变量(被测物的浓度或含量等);
A—校准曲线的斜率;
B—校准曲线的截距。
2)曲线上查得值
不确定度评定
假定对样品进行k次测量,测量响应值为Z,从校准曲线上查得被测量x,其平均值为
。
按照式(10)计算校准曲线的残余标准偏差s(残余):
……………(10)
按照式(11)计算校准曲线拟合引入的
的标准不确定度
:
…………(11)
式中:
不确定度的B类评定
需要量化的不确定度分量如下:
a)对测量结果有较大影响的,未包括在A类不确定度评定中的分量;
b)通常,影响因素较少时,占不确定度最大分量
以下的因素可不予考虑。
5.3.2.1利用标准物质标定仪器引入得不确定度
假定用标准物质标定仪器,重复测定n次,按如下方法进行B类不确定度的评估:
1)计算标准值不确定度引入的分量:
当不确定度以标准差表示时,直接引用,Si=不确定度;
当不确定度以展伸不确定度表示时:
包含因子由测量数据的分布类型和置信概率决定,包含因子Kp与分布类型和置信概率P的关系见表1和表2。
煤炭检验中,标准物质测定一般为正态分布,置信概率取95%。
2)计算标准物质的重复性测量不确定度:
式中:
xi——标准物质的单次测定值;
——标准物质n次测定值的平均值;
n——标准物质的测量次数。
3)按下式计算标准物质标定仪器引入的不确定度:
5.3.2.2根据极限值和误差分布知识进行评估
式中:
——误差极限值
表1正态分布下置信概率P与包含因子kP间的关系
P(%)
50
68.27
90
95
95.45
99
99.73
kP
0.67
1
1.645
1.960
2
2.576
3
表2常用分布与kP的关系
分布类型
P(%)
kp
正态
95
1.96
三角
100
梯形
100
矩形(均匀)
100
反正弦
100
两点
100
1
*ß为梯形的上底和下底之比。
注:
在缺乏信息的情况下,一般估计为矩形分布是较合理的。
但如果xi是重复性条件下的几个测量值的算术平均值,则可估计为正态分布。
5.3.2.3按制造部门的说明书,校准证书、或其他资料给出的不确定度进行评估
1)若制造部门的说明书,校准证书、或其他资料给出的估计值xi的不确定度U(xi)是标准偏差s(xi),指明了包含因子k,则估计值的标准不确定度u(xi)=s(xi)/k。
2)若制造部门的说明书,校准证书、或其他资料给出了估计值xi在一定置信概率P的置信区间的半宽U(xi),则估计值的标准不确定度u(xi)=
。
3)当测量仪器检定证书上给出准确度等级时,可按检定规程所规定的该等级的测量不确定度或最大允许误差大小,按5.3.2.4.1或5.3.2.4.2进行评定。
4)对于数字显示式测量仪器,如其分辨力为δx,则由此带来的标准不确定度为u(x)=0.29(δx)。
5.3.2.4当某不确定度影响因素的相应试验方法国家标准或行业标准中,明确指出该因素的测量重复性限r时,则该因素的标准不确定度按照式(12)计算:
u(xi)=
………………(12)
5.3.2.5煤炭检验中常见的几种B类不确定度因素的评定
5.3.2.5.1质量标准不确定度
煤炭检验中,质量称量标准不确定度um主要考虑称量的变动性和天平校准产生的不确定度。
通常以称量误差Δ(如:
0.0002g)为误差极限值,并按下式计算um:
………………………………(13)
5.3.2.5.2体积标准不确定度
体积标准不确定度主要考虑容器体积校准、定容重复性和温度差异三个方面。
1)容器体积校准不确定度uV
2)定容重复性不确定度Sr
通过对容器进行n次定容和称量,以称量值为x,按5.3.1.1中的公式
(2)计算定容重复性标准偏差,即为定容重复性不确定度Sr。
3)温度差异的不确定度
式中:
D——体积膨胀系数,水的体积膨胀系数是2.1×10-4ºC-1;
Δt——试验温度与校准温度之差,取绝对值;
V——定容的体积。
4)体积合成标准不确定度
按照式(14)计算体积标准不确定度uL:
……………………(14)
式中:
uL——体积标准不确定度;
Sr——定容重复性标准偏差;
5.3.2.5.3试剂纯度
试剂纯度的标准不确定度uP,通常为均匀分布,按照式(15)计算:
……………………(15)
例:
铬的纯度为0.9999,其标准不确定度为:
合成标准不确定度
5.3.3.1各不确定度分量彼此独立
合成标准不确定度uc按式(16)进行评定:
………………(16)
式中:
y—测量值;
L—影响测量的独立参数的个数;
ci—各参数的灵敏度系数,其值如式(17)所示:
……………………………(17)
xi—影响测量的独立参数。
合成标准不确定度在大多数情况下可简化为以下条款的形式进行评定。
5.3.3.1.1若独立参数对测量值的影响仅涉及量的加或减,如y=(x1+x2-x3+…),可用式(18)计算合成标准不确定度:
…………(18)
5.3.3.1.2若独立参数对测量值的影响仅涉及量的乘或除,如y=x1/(x2×x3×…)或y=(x1×x2×x3×…),可用式(19)计算合成标准不确定度:
………(19)
5.3.3.1.3若独立参数对测量值的影响涉及量的加、减、乘或除,则按照计算的先后原则,将原有的数学模型分解成仅包括上述5.3.3.1.1和5.3.3.1.2中一种形式的单元,按5.3.3.1.1和5.3.3.1.2计算每个单元的不确定度,再将这些不确定度合成为合成标准不确定度。
例:
y=(x1+x2)×x3,则先按5.3.3.1.1计算p、q的合成标准不确定度为:
再按5.3.3.1.2计算y的合成标准不确定度为:
5.3.3.1.4当不确定度诸影响因素具有不同量纲时,可计算出各影响因素的相对不确定度,再按公式(19)计算合成不确定度。
注:
煤炭检验中的大多数分析项目,其各不确定度影响因素具有不同的量纲,此时按相对不确定度方式计算合成不确定度较为简单、方便。
5.3.3.2不确定度分量彼此具有相关性
5.3.3.2.1不确定度分量中有两个分量彼此具有相关时,按下式计算合成标准不确定度:
……(20)
式中:
u(x2)和u(x3)——两个彼此具有相关性的不确定度分量。
5.3.3.2.2若更多的不确定度分量彼此相关,应采用专门为相关的随机变量的数据处理设计的统计方法来分析测量值,其相关性影响的完整处理已超出本标准的范围。
扩展不确定度
5.3.4.1扩展不确定度U为合成标准不确定度uc(y)与一个包含因子k的乘积,如下式:
U=k·uc(y)
k的取值与置信概率有关。
一般在置信概率约为95%时,取k=2;置信概率约为99%时,取k=3。
5.3.4.2利用标准方法中重复性限r或再现性临界差R评估报告结果的不确定度
使用煤炭分析试验方法的国家标准或行业标准,则当报告的结果为同一化验室两个重复性测量值的平均值时,该结果的扩展不确定度可按照下式计算:
…………………………………………………………………………(21)
式中:
r——标准方法规定的重复性限。
报告结果为两个不同化验室再现性测定值的平均值时,该结果的扩展不确定度可按照下式计算:
………………………………………………………………………(22)
式中:
R——标准方法规定的再现性临界差。
6测量不确定度的报告与表示
6.1完整的不确定度报告一般应包括以下内容:
a)试验方法、试验原理和试验过程的说明;
b)不确定度源的识别;
c)标准不确定度的A类评定;
d)标准不确定度的B类评定;
e)合成标准不确定度的评定;
f)扩展不确定度的评定(若有)。
6.2除非另有要求,测量结果的不确定度应用扩展不确定度U报告,采用以下表示方式:
y±U
并说明k值的大小。
必要时,给出完整的不确定度报告。
例:
校准容器的体积为V=100.03mL,合成标准不确定度u(V)=0.03mL,则表示为:
V=(100.03±0.06)
K=2
6.3在日常的例行检验中,虽无明确的不确定度报告,但试验严格遵照标准方法规定的测量程序,仪器设备均经检定处于合格状态,则测量结果的不确定度也可表示为“以重复性限表示的测量结果的不确定度符合有关标准的要求。
”
6.4当用合成标准不确定度uc(y)报告测量结果的不确定度时应把测量结果和标准不确定度分开来表示。
例:
7.2中的例子可表示为:
容器体积:
V=100.03mL
标准不确定度:
uc(V)=0.03mL
6.5不确定度通常最多保留两位有效数字,测量结果的位数与不确定度位数平齐。
附录A
(资料性附录)
测量不确定度评定举例
A.1煤中灰分测量不确定度评定
某煤样的灰分测定数据如下:
灰皿质量,g17.265516.0562
试样质量,g1.00001.0044
加热后(灰皿+试样)质量,g17.352816.1438
灰的质量,g0.08730.0876
灰分产率,%8.738.70
平均Aad,%8.72
A1.1被测量说明
A1.1.1测量程序
称量灰皿质量称取1g左右煤样,称准至0.0002g815℃±10℃灼烧至质量恒定称量(灰皿+试样)质量计算灰分含量Aad
A1..12被测量
式中:
m1-试样灼烧后灰的质量,g;
m-称取的试样质量,g。
Aad-试样的灰分含量,%。
A1.2不确定度源分析
灰分测定值不确定度主要来自:
1)A类不确定度(测量重复性)
2)试样称量
3)质量损失称量
4)质量恒定
用函数形式表示为:
y(Aad)=f(m,m1,质量恒定,重复性)
A1.3不确定度组分的评定
A1.3.1A类不确定度分量(测量重复性)
A1.3.2B类不确定度
a)试样质量m
试样称量误差0.0002g,按均匀分布转化为标准不确定度:
u(m)=0.0002/
=0.0001154g
b)质量损失m1
灼烧后称量误差0.0002g,按均匀分布转化为标准不确定度:
u(m1)=0.0002/
=0.0001154g
c)质量恒定
本例灰分测量值8.72%,不需进行检查性灼烧,可不考虑质量恒定的不确定度。
A1.3.3合成标准不确定度
合成标准不确定度各组分:
不确定度组分量值标准不确定度相对标准不确定度
u(x)u(x)/x
试样称量1.0000g0.0001154g0.000115
灼烧后称量0.0873g0.0001154g0.000132
重复性8.72(%)0.0224(%)0.00256
标准合成不确定度:
相对不确定度:
u’(Aad)=
=0.0026
绝对不确定度:
u(c)=8.72×0.0026=0.023(%)
A1.3扩展不确定度
取扩展因子k=2,
灰分测量值Aad的扩展不确定度:
U(Vad)=0.023×2=0.046=0.05
A1.4结果表述
该煤样的灰分测量结果为:
(8.72±0.05)%
A2煤的发热量测定结果的不确定度评定
用恒温式氧弹量热量计测定某煤样的发热量含量,所得试验数据如下:
1)热容量标定
苯甲酸标准热值:
Q=26459J/g,点火热:
25J
序号苯甲酸质量(g)测量温升(K)冷却校正值(K)测得热容量E(J/K)
10.96432.42010.022310463
21.02792.57950.026010458
30.80202.01820.017210453
40.98572.47820.020910462
50.98242.46860.021910458
平均值0.95252.39290.021710459
标准差3.96
2)发热量测定
仪器热容量:
E=10459J/K,点火热:
25J
试样质量(g)测量温升(K)冷却校正值(K)测得发热量Qb,ad(J/g)
10.99841.83090.013619298
21.00631.84560.015819321
平均值1.00241.83820.014719310
全硫St,ad=3.33%
高位发热量Qgr,ad=18973J/g
A2.1被测量说明
A2.1.1测量程序
1)仪器热容量标定
称量已知热值的标准苯甲酸片→置于浸没在一定量水中的氧弹内燃烧→测量燃烧前后的水温→根据水温的升高计算仪器的热容量。
-\
2)煤的发热量的测定
称取1g左右的试验煤样→置于浸没在一定量水中的氧弹内燃烧→测量燃烧前后的水温→根据水温的升高计算试样的弹筒发热量→从弹筒发热量中扣除硫酸和硝酸校正热计算高位发热量。
A2.1.2被测量
1)热容量E
式中:
Q——苯甲酸标准热值,J/g;
mb——称取的苯甲酸片质量,g;
qi———点火热,J;
qn——硝酸形成热,Jqn=0.0015*Q*m。
2)弹筒发热量
:
式中:
E——热量计的热容量,J/K;
qi——点火热,J;
q2——添加物如包纸等产生的总热量,J;
m——试样质量,g;
to——点火时的内筒温度,K;
tn——终点时的内筒温度,K。
3)高位发热量Qgr,ad:
Qgr,ad=Qb,ad-(94.1Sb,ad+αQb,ad)
式中:
Qgr,ad——空气干燥煤样的恒容高位发热量,J/g;
Qb,ad——分析试样的弹筒发热量,J/g;
Sb,ad——由弹筒洗液测得的煤的含硫量,%;
当全硫含量St,ad低4.00%时,或发热量大于14.60MJ/kg时,用全硫St,ad代替Sb,ad;
94.1——煤样中每1.00%硫的校正值,J;
α——硝酸生成热校正系数:
当Qb≤16.70MJ/kg,α=0.0010;
当16.70MJ/kg<Qb≤25.10MJ/kg,α=0.0012;
当Qb>25.10MJ/kg,α=0.0016。
A2.2不确定度来源分析
A2.2.1热容量的影响参数有:
苯甲酸标准热值,苯甲酸片的质量,测量温升,冷却校正值(恒温式量热仪),内筒水量,点火热:
y(E)=f(Q,mb,Δt,C,qn,,qi)
A2.2.2弹筒发热量的影响参数有:
热容量,试样称量,测量温升,冷却校正值(恒温式量热仪),内筒水量,点火热:
y(Qb,ad)=f(E,Δt,C,m,qi)
热容量标定和发热量测定中采用同样的点火热,它对弹筒发热量测定结果的影响基本抵消,可不