CN105717006A - 一种激光粒度分析仪测量结果不确定度评定方法 - Google Patents
一种激光粒度分析仪测量结果不确定度评定方法 Download PDFInfo
- Publication number
- CN105717006A CN105717006A CN201610178470.6A CN201610178470A CN105717006A CN 105717006 A CN105717006 A CN 105717006A CN 201610178470 A CN201610178470 A CN 201610178470A CN 105717006 A CN105717006 A CN 105717006A
- Authority
- CN
- China
- Prior art keywords
- uncertainty
- particle size
- laser particle
- standard substance
- size analyzer
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000002245 particle Substances 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims abstract description 25
- 239000000126 substance Substances 0.000 claims abstract description 46
- 238000005259 measurement Methods 0.000 claims abstract description 33
- 238000004458 analytical method Methods 0.000 claims abstract description 19
- 239000004531 microgranule Substances 0.000 claims description 8
- 239000002131 composite material Substances 0.000 claims description 7
- 230000003252 repetitive effect Effects 0.000 claims description 5
- 238000000926 separation method Methods 0.000 claims description 5
- 230000009897 systematic effect Effects 0.000 claims description 4
- 239000000463 material Substances 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 abstract description 19
- 238000004364 calculation method Methods 0.000 abstract description 10
- 238000003921 particle size analysis Methods 0.000 abstract 1
- 238000010561 standard procedure Methods 0.000 abstract 1
- 238000013076 uncertainty analysis Methods 0.000 abstract 1
- 238000012360 testing method Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 239000010421 standard material Substances 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/02—Investigating particle size or size distribution
- G01N15/0205—Investigating particle size or size distribution by optical means
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明属于计算测量与统计领域,具体涉及一种激光粒度分析仪测量结果不确定度评定方法。包括以下步骤:(1)选取三个粒度不同的国家微粒标准物质,用激光粒度分析仪进行粒度分析;(2)统计影响测量结果的各项不确定度分量,由偶然误差引起的不确定度<i>U1</i>;由仪器系统误差引起的不确定度<i>U2</i>;由标准物质引入的不确定度<i>U3</i>;由激光粒度分析仪中的传感器布局引起的不确定度<i>U4</i>;(3)进行不确定度分析。本发明为评价激光粒度仪测量结果的可靠性提供了一种科学规范的方法;不同用户使用的激光粒度仪可以采用本发明提供的方法进行评估,以评价分析结果的可靠性;同时,使不同用户的分析结果具有了可比性。
Description
技术领域
本发明属于计算测量与统计领域,具体涉及一种激光粒度分析仪测量结果不确定度评定方法。
背景技术
激光粒度分析仪由于其使用方便,分析速度快,广泛用于粉体材料的粒度分析。由于激光粒度仪的应用厂家众多,各应用厂家使用的激光粒度仪生产厂家、规格不同、使用条件不同、操作人员的技术也不同,这些差异都会导致分析结果的差异,到底哪一家的仪器分析结果可靠,必须通过不确定度评估,根据不确定度评估结果判断准确与否。同时GB/T27025-2008/ISO/IEC17025:2005《检测和校准实验室能力的通用要求》也要求实验室对检测仪器设备,都应进行不确定度的评估。
由于激光粒度仪的不确定度的评估比较复杂,目前没有统一的不确定度评估方法。通常的不确定评估都是对一个点的测量不确定度进行评估,由于激光粒度仪进行粒度分析时,分析的是一个粒度范围而不是一个点,激光粒度仪分析误差分布也不是线性的,因此只对一个点进行评估是不够的;为解决这一问题,本专利提出了激光粒度仪不确定度的三点评估方法。
发明内容
本发明所要解决的技术问题是提供一种科学、完整、准确的激光粒度分析仪测量结果不确定度评定方法。
本发明所述的一种激光粒度分析仪测量结果不确定度评定方法,包括以下步骤:
1)选取三个粒度不同的国家微粒标准物质,用激光粒度分析仪进行粒度分析;
2)统计影响测量结果的各项不确定度分量,具体为:
2.1)由偶然误差引起的不确定度U1;即,由多次重复测量引起的不确定度,以均方差作为重复测量引起的不确定度,其计算公式为:
其中:Di-----标准物质单次测量值;
-----标准物质多次测量平均值;
n-----测量次数;
2.2)由仪器系统误差引起的不确定度U2;即,由仪器对标准物质多次测量结果的平均值与标准物质额定值之间的差别引起的不确定度,其计算公式为:
其中:-----标准物质多次测量平均值;
di------标准物质额定值;
2.3)由标准物质引入的不确定度U3;采用国家标准物质给定的不确定度计算方法,其计算公式为:
其中:Ub-----标准物质的质量中值不确定度;
k-----包含因子,取k=2;
2.4)由激光粒度分析仪中的传感器布局引起的不确定度U4;激光粒度分析仪中有两个传感器计数通道,以标准物质给定值所处的计数通道前后分界点所表征的粒径值之差的二分之一作为U4的计算标准,其计算公式为:
其中:d1,d2----标准物质给定值所处的计数通道前后分界点所表征的粒径值;
3)不确定度分析:
3.1)单一测量点的合成不确定度Ud50;计算公式为:
注:本文中的“合成不确定度”的含义是“当测量结果是由若干个其他量的值求得时,按其他各量的方差算得的不确定度”,与“标准合成不确定度”的含义是一样的,但目前还没有激光粒度仪不确定度的“标准合成不确定度”,所以本文采用“合成不确定度”来表示。
3.2)单一测量点的相对合成不确定度Ud;计算公式为:
其中:Ud-----单一测量点的相对合成不确定度
3.3)单一测量点的相对扩展不确定度Uk为:计算公式为:
Uk=K×Ud(7)
其中,K-----包含因子,本方法评估的不确定度以包含概率95%为目标,K取2。
3.4)仪器综合不确定度U;计算公式为:
U=(Uk1+UK2+UK3)/3(8)
其中:U-----仪器综合不确定度。
进一步,用质量中值d50作为步骤2.2)中标准物质额定值di。
与现有技术相比,本发明具有以下积极有益效果:
(1)本发明对激光粒度仪测量结果提出了不确定度的评定方法,为激光粒度仪测量结果不确定度评定建立了函数关系。
(2)本发明为评价激光粒度仪测量结果的可靠性提供了一种科学规范的方法;不同用户使用的激光粒度仪可以采用本专利提供的方法进行评估,以评价分析结果的可靠性;同时,使不同用户的分析结果具有了可比性。
具体实施方式
为了使本发明的技术方案及优点更加清楚明了,结合一下实施例,对本发明进一步详细描述,但本发明的保护范围不限于下述实施例。
实施例1
评估仪器:1号仪器
仪器名称:激光粒度分析仪
规格型号:Mastersizer2000
测试范围:1~30μm
标准物质:三个国家微粒标准物质,编号分别是:
GB(W)120002a,质量中值d50=3.6μm
GB(W)120004a,质量中值d50=9.8μm
GB(W)120006a,质量中值d50=22.4μm
粒度分析:按照仪器说明书的要求,对三个标准物质进行粒度分析,分别取每个标准物质的质量中值d50作为标准物质额定值di进行评估。
以国家微粒标准物质GB(W)120002a的评定结果为例,对测量不确定度各分量的评定如下:
(1)由偶然误差引起的不确定度U1的评定
采用激光粒度分析仪测定国家微粒标准物质的粒度组成,对三个微粒标准物质,进行10次独立的重复测量,以均方差作为重复测量引起的不确定度,利用公式(1)进行计算,其计算结果为:
(2)由仪器系统误差引起的不确定度U2的评定
系统误差可通过计算10次测量结果平均值与额定值d50之差得出,利用公式(2)进行计算,其计算结果为。
(3)由标准物质引入的不确定度U3的评定
国家微粒标准物质GB(W)120002a证书给出的质量中值不确定度为0.2,不确定因子k=2,代入公式(3)进行计算,其计算结果为:
(4)由激光粒度分析仪中的传感器布局引起的不确定度U4的评定
由激光粒度分析仪中的传感器布局引起的不确定度U4,激光粒度分析仪中有两个传感器计数通道,以标准物质给定值所处的计数通道前后分界点所表征的粒径值之差的二分之一作为计数通道不确定度,检测报告表明,对应于d50=3.6μm传感器通道布局对应d2=3.802,d1=3.311。将该数据代入公式(4),其计算结果为:
不确定度分析:
(1)合成不确定度Ud50的评定
将上述各不确定度分量代入公式(5),其计算结果为:
(2)相对合成不确定度Ud的评定
以标准物质额定值为准,评定出相对不确定度,利用公式(6)进行计算,其计算结果为:
(3)相对扩展不确定度Uk的评定,利用公式(7)进行计算,其计算结果为:
取包含因子K=2.则
Uk=9.50%×2=19.01%
其余各测量点的数据,依据上述方法进行评定,评定结果见表1。
(4)仪器综合不确定度
以三个标准物质测量点的相对扩展不确定度的平均值作为仪器的综合不确定度。不确定度的覆盖范围,在最小测量点的测量值与最大测量点的测量值之间。将三个标准物质测量点的相对扩展不确定度,代入将公式(8),其计算结果为:
U=(Uk1+UK2+UK3)/3=(19.01%+15.84%+13.93%)/3=16.26%
实施例1中三个标准物质的测量和分析结果如表1所示:
表1实施例1中三个标准物质的测量和分析结果
实施例2
评估仪器:2号仪器
仪器名称:激光粒度分析仪
规格型号:SCF-106
测试范围:1~30μm
标准物质:三个国家微粒标准物质,编号分别是:
GB(W)120002a,质量中值d50=3.6μm
GB(W)120004a,质量中值d50=9.8μm
GB(W)120006a,质量中值d50=22.4μm
粒度分析:按照激光粒度仪说明书的操作要求,对三个标准物质进行粒度分析,分别取每个标准物质的质量中值d50作为标准物质额定值di进行评估。分析结果如表2所示。具体计算方法参见实施例1。
表2实施例2中三个标准物质的测量和分析结果
验证实施
验证仪器:1号仪器、2号仪器
标准物质:GB(W)120005,质量中值d50=16.3μm
验证方法:按照仪器说明书的要求,对标准物质GB(W)120005进行粒度分析,取质量d50值作为验证数据,测试结果如表3所示。
表3统计与评定结果
仪器序号 | 1号仪器 | 2号仪器 |
额定值 | 16.3 | 16.3 |
测量数据 | 15.72 | 15.30 |
偏差 | 0.58 | 1.00 |
实例1和实例2对两台仪器的不确定度评估结果表明,1号仪器的不确定度小于2号仪器的不确定度。表3分析验证实验1号仪器的测量偏差小于2号仪器的测量偏差。两台仪器测量结果与本专利的不确定度评估结果是一致的,验证结果表明本专利提出的不确定度评估方法是合理可行的。
Claims (2)
1.一种激光粒度分析仪测量结果不确定度评定方法,其特征在于,包括以下步骤:
1)选取三个粒度不同的国家微粒标准物质,用激光粒度分析仪进行粒度分析;
2)统计影响测量结果的各项不确定度分量,具体为:
2.1)由偶然误差引起的不确定度U1;其计算公式为:
其中:Di-----标准物质单次测量值;
-----标准物质多次测量平均值;
n-----测量次数;
2.2)由仪器系统误差引起的不确定度U2;其计算公式为:
其中:-----标准物质多次测量平均值;
di-----标准物质额定值;
2.3)由标准物质引入的不确定度U3;其计算公式为:
其中:Ub-----标准物质的质量中值不确定度;
k-----包含因子,取k=2;
2.4)由激光粒度分析仪中的传感器布局引起的不确定度U4;激光粒度分析仪中有两个传感器计数通道,以标准物质给定值所处的计数通道前后分界点所表征的粒径值之差的二分之一作为U4的计算标准,计算公式为:
其中:d1,d2----标准物质给定值的计数通道分界点的粒径值;
3)不确定度分析:
3.1)单一测量点的合成不确定度Ud50;计算公式为:
3.2)单一测量点的相对合成不确定度Ud;计算公式为:
其中:Ud-----单一测量点的相对合成不确定度
3.3)单一测量点的相对扩展不确定度Uk为:计算公式为:
Uk=K×Ud(7)
其中,K-----包含因子,本方法评估的不确定度以包含概率95%为目标,K取2。
3.4)仪器综合不确定度U;计算公式为:
U=(Uk1+UK2+UK3)/3(8)
其中:U-----仪器综合不确定度。
2.根据权利要求1所述的一种激光粒度分析仪测量结果不确定度评定方法,其特征在于,用质量中值d50作为步骤2.2)中标准物质额定值di。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610178470.6A CN105717006B (zh) | 2016-03-28 | 2016-03-28 | 一种激光粒度分析仪测量结果不确定度评定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610178470.6A CN105717006B (zh) | 2016-03-28 | 2016-03-28 | 一种激光粒度分析仪测量结果不确定度评定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105717006A true CN105717006A (zh) | 2016-06-29 |
CN105717006B CN105717006B (zh) | 2018-11-09 |
Family
ID=56158302
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610178470.6A Active CN105717006B (zh) | 2016-03-28 | 2016-03-28 | 一种激光粒度分析仪测量结果不确定度评定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105717006B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106644864A (zh) * | 2016-12-05 | 2017-05-10 | 中国计量科学研究院 | Pm2.5切割器切割特性检测用标准物质的定值方法及测量系统 |
CN107271345A (zh) * | 2017-06-15 | 2017-10-20 | 中国食品药品检定研究院 | 一种两法联合标定氧气透过量标准膜及其制备方法 |
CN109696586A (zh) * | 2017-10-23 | 2019-04-30 | 中国电力科学研究院有限公司 | 一种用于测试配电自动化系统不确定度的评定方法及系统 |
CN112923967A (zh) * | 2021-05-10 | 2021-06-08 | 北京康斯特仪表科技股份有限公司 | 支持对仪表校验不确定度进行快速评定的仪表校验系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101603958A (zh) * | 2009-07-16 | 2009-12-16 | 朱家平 | 分析不确定度的方法及装置 |
CN102135412A (zh) * | 2011-01-05 | 2011-07-27 | 西安交通大学 | 一种基于不确定度分布规律的工件组合测量方法 |
CN102565857A (zh) * | 2011-12-16 | 2012-07-11 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 自动剩余动校正方法 |
-
2016
- 2016-03-28 CN CN201610178470.6A patent/CN105717006B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101603958A (zh) * | 2009-07-16 | 2009-12-16 | 朱家平 | 分析不确定度的方法及装置 |
CN102135412A (zh) * | 2011-01-05 | 2011-07-27 | 西安交通大学 | 一种基于不确定度分布规律的工件组合测量方法 |
CN102565857A (zh) * | 2011-12-16 | 2012-07-11 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 自动剩余动校正方法 |
Non-Patent Citations (2)
Title |
---|
张文阁: "激光粒度分析仪校准方法及校准结果不确定度分析", 《中国粉体技术》 * |
李士林: "激光粒度分析仪粒径示值误差的测量结果不确定度评定", 《计量与测试技术》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106644864A (zh) * | 2016-12-05 | 2017-05-10 | 中国计量科学研究院 | Pm2.5切割器切割特性检测用标准物质的定值方法及测量系统 |
CN107271345A (zh) * | 2017-06-15 | 2017-10-20 | 中国食品药品检定研究院 | 一种两法联合标定氧气透过量标准膜及其制备方法 |
CN109696586A (zh) * | 2017-10-23 | 2019-04-30 | 中国电力科学研究院有限公司 | 一种用于测试配电自动化系统不确定度的评定方法及系统 |
CN112923967A (zh) * | 2021-05-10 | 2021-06-08 | 北京康斯特仪表科技股份有限公司 | 支持对仪表校验不确定度进行快速评定的仪表校验系统 |
Also Published As
Publication number | Publication date |
---|---|
CN105717006B (zh) | 2018-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Thomas et al. | Excess clustering on large scales in the MegaZ DR7 photometric redshift survey | |
CN105528288B (zh) | 一种软件测试方法以及装置 | |
CN105717006A (zh) | 一种激光粒度分析仪测量结果不确定度评定方法 | |
CN104713845B (zh) | 基于太赫兹吸收光谱数据处理的混合物成分识别方法 | |
CN102749378B (zh) | 模拟尾气排放的汽油车简易瞬态工况排放检测标准装置 | |
CN103940767B (zh) | 基于多流形学习的气体浓度反演方法 | |
CN104569314B (zh) | 一种商品化定量检测试剂盒评价方法 | |
CN110441241A (zh) | 一种光声光谱多组分气体分析仪器的性能评价装置及方法 | |
CN103902451A (zh) | 一种智能电能表软件品质评价方法 | |
EP3410078A3 (en) | Apparatus and methods for calibrating and/or validating pulmonary function test equipment | |
CN109358095B (zh) | 一种混合气体各组分气体浓度的定量测量方法和系统 | |
CN104198416A (zh) | 一种光谱仪波长漂移引起的测量误差实时补偿方法 | |
CN104569886B (zh) | 基于时频参数标准仪器方式的信号检测设备校准方法 | |
CN202661459U (zh) | 一种模拟尾气排放的汽油车瞬态工况排放检测标准装置 | |
CN104833652A (zh) | 一种混合材料中纤维素硝酸酯含量的快速分析方法 | |
CN103995147B (zh) | 一种适用于声学多普勒流速仪的数据后处理系统与应用 | |
CN113252281B (zh) | 一种结冰云雾液滴尺寸分布的重构方法 | |
CN104316511B (zh) | 火花源原子发射光谱分析中谱线干扰的校正方法 | |
CN104502875B (zh) | 基于时频参数标准信号源方式的信号检测设备校准方法 | |
CN109222220A (zh) | 用于评测烟支卷制质量指数的方法和系统 | |
CN102023135A (zh) | 偏峰定位光谱分析方法 | |
CN103616345B (zh) | 红外碳硫分析仪的含碳量测定方法 | |
CN107132302A (zh) | 一种外标法计算毒品组分含量检验结果不确定度的评定方法 | |
CN102213674A (zh) | 近红外光谱法检测面粉中石灰类物质的方法 | |
CN103983535A (zh) | 非标准大气压下卷烟焦油量和一氧化碳量测定的修正方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |