线性回归标准曲线法不确定度检验检疫.docx
《线性回归标准曲线法不确定度检验检疫.docx》由会员分享，可在线阅读，更多相关《线性回归标准曲线法不确定度检验检疫.docx（10页珍藏版）》请在冰豆网上搜索。
线性回归标准曲线法不确定度检验检疫
仪器分析中线性回归标准曲线法分析结果不确定度评估
一、前言
对测试方法制定不确定度评估程序是ISO/IEC17025对实验室的要求[1]，也是检验工作的需要。
由ISO等7个国际组织联合发布的《测量不确定度表达指南》[2]采用当前国际通行的观点和方法，使涉及测量的技术领域和部门可以用统一的准则对测量结果及其质量进行评定、表示和比较，满足了不同学科之间交往的需要[3]。
采用《测量不确定度表达指南》对测试结果不确定度进行评估，也是检验工作同国际标准接轨的需要。
线性回归标准曲线法是仪器分析中最常用的方法，这类仪器包括原子吸收分光光度计、发射光谱仪、分光光度计、气相（液相）色谱仪等。
这类分析测定结果的不确定度都有相似的来源，可概括为仪器精密度、标准物质不确定度及溶液制备过程中带来的不确定度等。
因此，可用相似的方法对它们进行评估。
本文以ICP-AES法测定钢铁中磷为例，推导了仪器分析中线性回归标准曲线法测定不确定度的计算方法，并提供了计算过程所需的各参数的采集和计算方法，评估了标准不确定度、自由度和扩展不确定度的数值。
二、测定过程和数学模型
仪器分析中线性回归标准曲线测定方法，利用被测物质相应的信号强度与其浓度成正比关系，通过测定已知浓度的溶液（即标准溶液）的信号强度，回归出浓度-信号强度标准曲线，从标准曲线上得到被测定溶液信号强度相应的浓度。
计算过程的数学模型如下：
用yi和yt分别表示标准溶液和被测溶液的信号线强度，以xi和xt分别表示第i个标准溶液和被测样品溶液的浓度，i=1~n，n表示标准溶液个数，则：
（1）
其中，
（2）
（3）
（1）式也可表示成：
（4）
把式
（2）、（3）代入式（4）得：
（5）
式（5）表明了被测量xt与输入量x1,x2...xn和y1,y2...yn、yt的函数关系，可简写成：
由上式可知，样品溶液浓度测定结果不确定度可分成标准溶液浓度不确定度分量及其信号强度不确定度分量和被测定溶液信号强度不确定度分量，其中标准溶液浓度不确定度分量可由标准样品标称含量不确定度和配制过程引入的不确定度合成得到，而信号强度不确定度分量是由仪器测量的误差引起的，可从仪器的精密度数据得到。
如果标准溶液是从各自不同的母标准物质配制，并且在配制过程使不同的器皿，则各标准溶液之间浓度估计值之间无关，因此xt各输入量估计值之间无关，由文献[2]5.1款式10，得样品溶液浓度的合成标准不确定度如下：
（6）
由式（5）得传播系数的表达式如下：
（7）
（8）
（9）
我们以测定钢铁中磷的浓度为例，选择5个标准钢铁样品，按SN/T0750-1999《进出口碳钢、低合金中铝、砷、铬、钴、铜、磷、钼、镍、硅、锡、钛、钒含量的测定-电感耦合等离子体原子发射光谱（ICP-AES）法》的方法制备标准溶液和被测样品溶液。
在ICP-AES上选用P213.6nm做分析谱线，分别测定标准溶液和被测溶液的分析谱线强度，用最小二乘线性回归法做出浓度-分析谱线强度标准曲线，从标准曲线上得到与被溶液分析谱线强度相应的浓度，再把被测溶液浓度转换成样品的百分比含量。
测试过程的各种量如表1所示。
表1钢铁中磷的测定有关量值
样品编号
标称含量（%）
称样量（g）
溶液浓度（ug/ml）
谱线强度（cps）
BH0507-9A
0.011
0.4755
0.523
57.6
30#第39号
0.022
0.4864
1.07
139.6
GB01202
0.029
0.5167
1.39
190
GB01203
0.038
0.5053
1.92
266.7
BH0506-3A
0.063
0.4524
2.58
351.9
被测样品
0.4805
2.30
316.7
三、各量值不确定度的计算
如前所述，仪器分析中标准曲线法测定结果不确定度来源包括仪器精密度、标准物质不确定度及制备各种溶液过程中引入的不确定度，这些不确定度来源可从仪器说明书及其它有关文件中得到，通过这些数据及合适的合成方法，可得到测定结果的不确定度。
我们以上面测定钢铁中磷的例子，各种量值不确定度计算如下：
1.标准溶液浓度不确定度u（xi）
标准溶液溶度由下式得出：
（10）
其中，xi为标准溶液浓度（单位为ug/ml），ci为标准样品质量比浓度（单位为1），mi为标样质量（单位为ug），vi为溶液体积（单位为ml）。
由此得出标准溶液浓度不确定度来自三个分量：
第一，钢铁标样标值的不确定度uc（xi）；第二，称量过程引入的不确定度um（xi）；第三，定容过程引入的不确定度uv（xi）。
由文献[3]公式10得：
（11）
式中传播系数：
（12）
（13）
（14）
（1）标值的标准不确定度uc（xi）
标准钢铁采用国家标准，其相对扩展不确定度为1%，按正态分布95%置信概率计算，标准不确定度：
uc（xi）=0.01/1.96ci
=0.0051ci
国标钢样由国家标准研究所提供，其标不准确定度可信度很高，自由度：
vc=∞
（2）称量过程引入的不确定度um（xi）
称量不确定度来自两个方面，第一，天平校正产生的不确定度，按检定证书上给出的在95%置信概率时为±0.1mg，换算成标准偏差为u1=0.1/1.96=0.052mg自由度为∞；第二称量变动性，根据历史记录，在50g以内，变动性标准偏差为u2=0.07mg，自由度为∞。
样品称量过程中，样品质量m=mG-mT得出，其中，mG是样品及容器总质量，mT是容器质量，由此得出样品质量不确定度的四个分量：
第一，系统误差导致的毛重不确定度uG1=0.052mg；第二，系统误差导致的皮重不确定度ut1=0.052mg；第三，由随机误差导致的毛重不确定度uG2=0.007mg；第四，由随机误差导致的皮重不确定度uT2=0.007mg。
由于毛重和皮重由同一天平称出，因此，系统误差导致的毛重和皮重值不确定度是强相关的，即相关系数为1，而随机误差导致的不确定度是无关的，由此得：
由文献[2]式（G2），自由度为：
（3）定容过程引入的不确定度uv（xi）[4]
使用100ml容量瓶配制标准和样品溶液，其不确定度包括三部分：
第一，容量瓶体积的不确定度，按制造商给定为±0.10ml，按均匀分布换算成标准偏差为：
，自由度为：
∞；第二，充满液体至容量瓶刻度的变动性，可通过重复称量进行统计，例如重复10次统计出标准偏差为0.012ml，自由度：
10-1=9；第三，容量瓶和溶液温度与校正时温度不同引起的体积不确定度，假设差3︒C，水体积膨胀系数为2.1X10-4/︒C，则95%置信概率时体积变化的区间为±100X3X2.1X10-4=0.063ml，转换成标准偏差：
0.063/1.96=0.032ml，自由为：
∞。
以上三项合成得出：
自由度为：
把5点标准溶液体积、称样量和标称含量代入式（12）、（13）、（14），得各量值的不确定度传播系数。
对标准溶液浓度不确定度有贡献的量、其不确定度及传播系数如表2所示。
把表2中有关量值代入式（11），得第i个标准溶液不确定度u（xi），其自由度v（xi）=∞，列于表3中。
表2标准溶液浓度不确定度u（xi）计算的有关量值
量
量值
标准不确定度
自由度
传播系数
不确定度分量（ug/ml）
c1
0.00011
0.00000056
∞
4755ug/ml
0.00269
c2
0.00022
0.00000112
∞
4864ug/ml
0.005466
c3
0.000269
0.00000137
∞
5167ug/ml
0.00709
c4
0.00038
0.00000194
∞
5053ug/ml
0.00979
c5
0.00063
0.00000321
∞
4524ug/ml
0.0145
m1
475500ug
10ug
∞
0.0000011ml-1
0.000011
m2
486400ug
10ug
∞
0.0000022ml-1
0.000022
m3
516700ug
10ug
∞
0.0000029ml-1
0.000030
m4
505300ug
10ug
∞
0.0000038ml-1
0.000038
m5
452400ug
10ug
∞
0.0000063ml-1
0.000063
v1
100ml
0.067ml
∞
-0.00523ug/ml
-0.00035
v2
100ml
0.067ml
∞
-0.0107ug/ml
-0.00072
v3
100ml
0.067ml
∞
-0.0150ug/ml
-0.00093
v4
100ml
0.067ml
∞
-0.0192ug/ml
-0.00129
v5
100ml
0.067ml
∞
-0.0285ug/ml
-0.00191
表3标准溶液浓度不确定度u（xi）及自由度v（xi）
量xi
量值
标准不确定度u（xi）
自由度v（xi）
x1
0.523ug/ml
0.00269ug/ml
∞
x2
1.07ug/ml
0.00551ug/ml
∞
x3
1.39ug/ml
0.00716ug/ml
∞
x4
1.92ug/ml
0.00989ug/ml
∞
x5
2.58ug/ml
0.0147ug/ml
∞
2.分析谱线强度的标准不确定度
分析谱线强度标准不确定度来源于ICP-AES仪器本身的长时间精密度。
以我们采用的美国PE公司OPTIMA3300RL为例，其仪器说明书上列明的短时间精密度为10次连续重复测量结果的相对标准偏差不超过1%，长时间精密度是每10分钟重复测量一次相同溶液，4小时内得到数据的相对极差小于10%，即相对标准偏差小于：
我们认为取长时间稳定性做仪器不确定度计算依据较为可信，其标准不确定度：
u（yi）=0.0202yi
把5点标准溶液的分析谱线强度代入上式，得其标准不确定度，如表4所示。
自由度：
v（yi）=4X60/10-1=23。
相似地，可得样品溶液分析谱线的强度不确定度：
u（yt）=0.0202yt,
自由度：
v（yt）=23
如表4所示。
表4溶液分析谱线强度不确定度u（yi）和自由度v（yi）
量
量值
标准不确定度u（yi）
自由度v（yi）
y1
57.6cps
1.16cps
23
y2
139.6cps
2.82cps
23
y3
190.0cps
3.84cps
23
y4
266.7cps
5.39cps
23
y5
351.9cps
7.10cps
23
yt
316.7cps
6.40cps
23
四、样品溶液浓度合成标准不确定度的计算
计算完上述有关量值，把它们分别代入式（7）、（8）、（9），得出各不确定度分量的传播系数。
以上述例，结果如表5所示。
表5被测样品溶液浓度不确定度u（xt）计算的有关量值
量
量值
标准
不确定度
自由度
传播系数
不确定度分量（ug/ml）
x1
0.523ug/ml
0.00269ug/ml
∞
-0.0445
-0.0012
x2
1.07ug/ml
0.00551ug/ml
∞
0.0904
0.000498
x3
1.39ug/ml
0.00716ug/ml
∞
0.165
0.00118
x4
1.92ug/ml
0.00989ug/ml
∞
0.301
0.00297
x5
2.85ug/ml
0.0147ug/ml
∞
0.488
0.00717
y1
57.6cps
1.16cps
23
0.000356ug/ml/cps
0.00414
y2
139.6cps
2.82cps
23
-0.000630ug/ml/cps
-0.00178
y3
190.0cps
3.84cps
23
-0.00120ug/ml/cps
-0.00461
y4
266.7cps
5.39cps
23
-0.00215ug/ml/cps
-0.0116
y5
351.9cps
7.11cps
23
-0.00334ug/ml/cps
-0.0237
yt
316.7cps
6.33cps
23
0.00696ug/ml/cps
0.0445
把表5有关量值代入式（6），得样品溶液浓度的合成标准不确定度：
u（xt）=0.0525ug/ml
自由度：
五、被测样品含量不确定度的计算
上面得出样品溶液浓度的不确定度，据此可计算出样品百分含量的不确定度。
钢铁样品的磷含量与样品溶液的关系式：
其中，ct为样品的质量分数含量（单位为1），xt为样品溶液浓度（单位为ug/ml），mt为样品质量（单位为ug，标准不确定度与标样相同，为10ug）。
代入表1的有关数据，可得：
ct=0.000479
由文献[2]公式10得：
（14）
把表1、2、3中有关量值代入上式，得：
自由度：
取置信概率p=95，自由度v=41，查t分布表得：
k=tp（v）=t0.95（41）=2.02
扩展不确定度：
U=tp（v）uc=2.02X0.000011=0.0022%
所以，例中被测样品中磷含量可表示成
ct=0.0479±0.0022（%）（其中扩展不确定包含系数k=tp（v）=t0.95（41）=2.02）。
参考文献
[1]ISO/IEC17025，校准和测试实验室能力的通用要求，北京：
中国国家进出口商品检验实验室认可委员会秘书处，1999年，1-81
[2]肖明耀等译，测量不确定度表达指南，北京：
中国计量出版社出版，1997年，1-164
[3]JJF1059-1999，测量不确定评定与表示，北京：
中国计量出版社出版，1999年，1-11
[4]韩永志，某些分析测定不确定度的计算示例，理化检验-化学分册，2000年，第36卷第9期，424-425
升级会员 