蒙特卡洛进行车内空气中甲苯浓度测定不确定度评定

车内空气中甲苯浓度测定不确定度评定
―――自适应蒙特卡洛法
一、检测方法
车内空气中的甲苯通过样品采集吸附在Tenax 管上，然后样品进行热解吸，采用气相色谱质谱联用法进行定性定量分析。

二、测量模型
00000
0.37()(273.15)
273.15101.325(273.15)
m a A b a A b t c PV V PfT
t +++=
==
+
式中：
c 车内空气中甲苯浓度，3/mg m 0V 标准采样体积，L m 样品质量， g µ 0a 标准曲线斜率 0b 标准曲线截距 A 样品峰面积 V 采样体积 L t 采样温度 ℃
P 采样现场大气压 kPa f 采样流量 /L min T 采样时间 min 三、各分量PDF 及其数字特征 数字特征
备注 分量 分布 均值
方差 a 2(,)N µσ
b
2(,)N µσ
A
样品不能重复进样，不引入的不确定 t
(,)U a b
P (,)U a b f 2(,)t νµσ
T
通过检定和长期观察，引入的不确定度可以忽略
四
、各分量MCM 模拟过程取值讨论
符号约定：meanNum ，NormalNum ，tdisNum 分别代表均匀（0,1）分布，标准正态分布以及t 分布数。

20(,)a N a σ ，由于在实际标准曲线制作过程中，2σ不能准确得到，
但是统计量，符合(2)t n −分布，因此，在MCM 过程：0ˆa a =+
符号符合通用现行回归符号 同理，在MCM 过程中，0b b =+ t 符合(,)a b 均匀分布，因此，()t a b a meanNum =+−× 同理()t a b a meanNum =+−×
f 分布信息来自于检定证书，其不确定以GUM 的格式给出，依据GUM
评定规范f 符合2(,)t νµσ，结合GUM 评定经验，自由度ν较大，依据数理统计原理，自由度较大的t 分布可当正态分布处理，因此，f normalNum σµ=+ 五、本评定例中输入量取值
标准曲线数据
甲苯含量g µ 0.050.1
0.250.51 信号强度 248750
5294221408911
2788054
5424160
7327
03
0 1.8410 1.73100.9991.8410
1.7310y x R a b −−−−=×−×==×=×
通过计算：
97
2.1610 1.8410
a a tdisNum−
=+=×+×
－
33 0
6.0810 1.7310
b b tdisNum−
=+×+×
－
采样温度t＝24.8℃，依据采样器校准证书，温度计的最大误差是4.0℃，按均匀分布处理，t＝24.6+0.4meanNum℃
采样压力P＝100.3kPa，气压计的最大误差是kPa
5.0，按均匀分布处理，
P＝100.05+0.5meanNum kPa
5.0
采样流量f＝0.106/L min，依据检定证书，4
7.5100.106
f normalNum
−
=×+峰面积A＝3858910，采样时间T＝30min
六、不确定度评定
1、取2
dig
n=
2、取95%
p=，
4
4
max(100/(1),10)
10
M p
M
=−
=
3、设h=1，经第一次MCM模拟，结果为：
(1)3
(1)33
3
0.24/
() 6.310/
(,)(0.231,0.255)/
low high
c mg m
u c mg m
c c mg m
−
=
=×
=
取最短包含区间
4、经计算()
u y的容差限是43
510/
mg m
−
×
5、依据MCM评定方法要求，当h=5时，
()
2,2,2,2
c u c clow chigh
s s s s皆小于容差限。

6、经过对h M
×次计算结果的不确定度评定，
3
33
3
0.24/
() 6.410/
(,)(0.231,0.255)/(P95%,
low high
c mg m
u c mg m
c c mg m
−
=
=×
==最短包含区间）。