CN104458783B - 基于限量管控应用的分析仪器测量结果自动判定方法 - Google Patents

基于限量管控应用的分析仪器测量结果自动判定方法 Download PDF

Info

Publication number
CN104458783B
CN104458783B CN201410620775.9A CN201410620775A CN104458783B CN 104458783 B CN104458783 B CN 104458783B CN 201410620775 A CN201410620775 A CN 201410620775A CN 104458783 B CN104458783 B CN 104458783B
Authority
CN
China
Prior art keywords
sample
physical quantity
detection
analytical instrument
surveyed
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.)
Expired - Fee Related
Application number
CN201410620775.9A
Other languages
English (en)
Other versions
CN104458783A (zh
Inventor
殷家祥
刘小东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
SHENZHEN UNIQUE METRICAL TECHNOLOGY CO LTD
Original Assignee
SHENZHEN UNIQUE METRICAL TECHNOLOGY CO LTD
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by SHENZHEN UNIQUE METRICAL TECHNOLOGY CO LTD filed Critical SHENZHEN UNIQUE METRICAL TECHNOLOGY CO LTD
Priority to CN201410620775.9A priority Critical patent/CN104458783B/zh
Publication of CN104458783A publication Critical patent/CN104458783A/zh
Application granted granted Critical
Publication of CN104458783B publication Critical patent/CN104458783B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Automatic Analysis And Handling Materials Therefor (AREA)

Abstract

本发明公开了一种基于限量管控应用的分析仪器测量结果自动判定方法,该方法采用能够得出每次检测的标准差的算法,得到样品的检测含量的同时可以得到含量的标准差,并基于得到的标准差与置信度进行结合,通过与样品中受限成分物理量的限值进行比较判断,进行样品合格或超标的判定。与现有技术相比,本发明有效解决了分析仪器的测量测量值落在限值附近,因随机误差存在,而导致的判断是否合格存在的风险问题。

Description

基于限量管控应用的分析仪器测量结果自动判定方法
技术领域:
本发明涉及分析仪器测量结果判定方法,具体涉及的是一种基于限量管控应用的分析仪器测量结果自动判定方法。
背景技术:
众所周知,分析仪器的测量结果总是带有系统误差和随机误差,通过合理的定量算法可以尽可能减小甚至消除系统误差(后面的描述都是基于不存在系统误差的情况),但随机误差必然会存在,只会由于具体方法的不同存在随机误差大小的问题。
目前限量管控在分析仪器的测量结果判断中应用较为广泛。限量管控是指规定了检测结果某个限值,含量高于限值即为不合格,含量低于限值为合格。而分析仪器的测量结果在采用限量管控应用进行判断是否合格时,如果测量值落在限值附近,而由于随机误差的存在,可能会不同的两次测量给出合格和不合格两种结论。而且对于不同的样品情况,随机误差的幅度又可能是不确定的,那么给合格性判定就会带来很大的困扰甚至风险。因此对于限量管控应用,只是简单地根据仪器的检测结果与限量值的对比来判定合格或超标,是有风险的。
为了避免上述问题,可在限量管控应用的基础上,给出合格性判断的置信水平和受限成分的含量检测的标准差(常用测量结果的标准差来表示随机误差的幅度)来对分析仪器的检测结果进行判断。比如某种检测仪器,在设定的条件下检测样品中某种受限成分的含量,限值为C0=10ppm,而根据仪器生产商的声明,受限成分的含量检测的标准差为σ=0.5ppm。假设某次检测结果为c=9.3ppm,那么对样品做出合格判定吗?如果判定合格,那么如果再一次测量,很可能得到高于10ppm的结果,又会判定为不合格。如此判定,很显然该仪器所作的合格性判定结果会被质疑。如果给出合格性判定的置信水平,比如按95%置信水平判断或99%置信水平判断,那么对于某一次的检测结果c=9.3ppm,其代表的含义是有95%的可能真实含量为c-1.97σ到c+1.97σ之间,有99%的可能为c-2.65σ到c+2.65σ之间,将c=9.3ppm和σ=0.5ppm带入,即真实含量有95%的可能为8.3到10.3之间,有99%的可能为8.0到10.6之间。这样并不能肯定地判定样品为合格,因为其真实值可能真的超过10ppm。相应的,如果检测结果低于C0-1.97σ即低于9.0,那么在95%的置信水平上样品是合格的,而如果检测结果低于C0-2.65σ即低于8.7,那么在99%的置信水平上样品是合格的。如果检测结果高于C0+1.97σ即低于11.0,那么在95%的置信水平上样品是超标的,而如果检测结果低于C0+2.65σ即低于11.3,那么在99%的置信水平上样品是超标的。否则,不能给出是否合格或超标的确定结论。在这种情况下,要么用精密度更好的检测仪器来进一步确定是否合格,要么,通过增加检测次数,缩小随机误差。
另外通常的情况下,标准差的获得是通过对样品进行n次测量,得到样本标准差Sn-1,而当n足够大时,Sn-1近似等于σ。这也是仪器生产商通常声明的检测标准差的来源。但这是不严谨的,一般情况下标准差是会随含量不同而变化的,并且因为不可能穷举用户的应用,检测用户样品时的标准差很可能与声明的不同,所以如果简单地采用上述解决方案,仍可能会得出错误的合格性判断。另外,实际的应用中,也许由于干扰检测的因素的存在,可能使检测结果的标准差严重偏离声明值,这时再以声明值为基础去做合格性判断,会给出风险极大的判断结果。比如用EDXRF来检测塑胶材料里面的有害元素铅和汞的含量(欧盟的RoHS指令,限量1000ppm),通常厂家会给出数ppm的标准差,通常情况下不会存在风险。但如果样品中含有溴(Br)基阻燃剂,由于Br的谱线很强,使得铅和汞的检测结果的标准差急剧增大到上百甚至几百ppm,这时如果仍然依据声明的数ppm的标准差去做合格性判断,具有巨大的风险。
发明内容:
为此,本发明的目的在于提供一种基于限量管控应用的分析仪器测量结果自动判定方法,以解决分析仪器的测量测量值落在限值附近,因随机误差存在,而导致的判断是否合格存在的风险问题。
为实现上述目的,本发明主要采用如下技术方案:
一种基于限量管控应用的分析仪器测量结果自动判定方法,包括步骤:
S10、设定样品中受限成分物理量的限值为C0
S20、获取样品的可测物理量值C1和该可测物理量值的标准差σ1
S30、通过标准样品对分析仪器进行标定,获取可测物理量值C1与不可测物理量值C2之间的线性关系公式;
S40、对未知样品进行检测,采用上述线性关系公式计算未知样品中不可测物理量值C2和该不可测物理量值的标准差σ2
S50、对未知样品中不可测物理量值C2进行判断,如果C2小于C0-f*σ2,则判定该未知样品合格;如果C2大于C0+f*σ2,则判定该未知样品不合格,其中f为置信度因子,95%置信度因子f=1.97,99%置信度因子f=2.65。
进一步地,步骤S50中如果C2介于C0-f*σ2与C0+f*σ2之间,则本次检测不能给出合格或超标的明确判断,此时可以通过多次检测,用平均值和平均值的标准差来进行判断,或采用精确度更高的检测方法给出确定性判断,为了在尽可能短的时间,给出可靠地合格性判断,可采用如下步骤:
S60、将用户能够接受的检测时间T0分成n段,以1/n*T0为检测时间段依次对不可测物理量值C2进行检测判断,如果当前检测时间段可以确定判断结果合格或不合格,则结束;否则,进行下一检测时间段继续进行检测,且将完成检测后的结果与前面的检测结果进行平均,根据平均结果判断合格或不合格,并依次类推,直至获得判断结果。
进一步地,步骤S20包括:
通过标准差算法获取分析仪器测量样品结果的可测物理量C1和该可测物理量的标准差σ1
进一步地,所述标准差算法为多变量非线性最小二乘拟合算法。
进一步地,步骤S30中所述标准样品为国家标准物质是标准样品或其物理量经过更高精密度的检测方法定值。
进一步地,步骤S30中所述线性关系公式为C2=a*C1+b,其中a和b为最小二乘拟合得到的参数。
进一步地,所述不可测物理量为用户需要了解的数据,其不可通过分析仪器检测获取。
本发明基于限量管控应用的分析仪器测量结果自动判定方法,采用能够得出每次检测的标准差的算法,得到样品的检测含量的同时可以得到含量的标准差,并基于得到的标准差与置信度进行结合,通过与样品中受限成分物理量的限值进行比较判断,进行样品合格或超标的判定。与现有技术相比,本发明有效解决了分析仪器的测量测量值落在限值附近,因随机误差存在,而导致的判断是否合格存在的风险问题。
附图说明:
图1为基于限量管控应用的分析仪器测量结果自动判定方法的流程图。
具体实施方式:
为阐述本发明的思想及目的,下面将结合附图和具体实施例对本发明做进一步的说明。
本发明实施例提供的基于限量管控应用的分析仪器测量结果判定方法,适用于任何分析仪器,其核心思想在于:在现有限量管控应用的基础上,给出合格性判断的置信水平,采用能够得出每次检测的标准差的算法,得到样品的检测含量的同时得到受限成分含量的标准差,然后基于得到的标准差与置信度进行结合计算,通过与样品中受限成分物理量的限值进行比较判断,从而解决分析仪器测量测量值落在限值附近时的样品合格或超标的判定。
请参见图1所示,图1为基于限量管控应用的分析仪器测量结果自动判定方法的流程图。本发明实施例提供的是一种基于限量管控应用的分析仪器测量结果自动判定方法,其主要包括如下步骤:
S10、设定样品中受限成分物理量的限值为C0
S20、获取样品的可测物理量值C1和该可测物理量值的标准差σ1
其中通过标准差算法获取分析仪器测量样品结果的可测物理量C1和该可测物理量的标准差σ1
S30、通过标准样品对分析仪器进行标定,获取可测物理量值C1与不可测物理量值C2之间的线性关系公式;
标定是通过对标准样品的检测,得到仪器检测的物理量(荧光强度、吸光度、发光强度等)与被测量的关系。只有经过标定得到的这个关系,仪器才能对未知样品进行检测。
标准样品为国家标准物质是标准样品或其物理量经过更高精密度的检测方法定值。而对于仪器进行标定则一般通过仪器制造商使用经过准确定量的系列标准样品实现,其对仪器进行出厂标定后交给用户使用,用户也可以使用仪器自带的标定功能和系列标准样品自行标定仪器。
S40、对未知样品进行检测,采用上述线性关系公式计算未知样品中不可测物理量值C2和该不可测物理量值的标准差σ2
S50、对未知样品中不可测物理量值C2进行判断,如果C2小于C0-f*σ2,则判定该未知样品合格;如果C2大于C0+f*σ2,则判定该未知样品不合格,如果C2介于C0-f*σ2与C0+f*σ2之间,则本次检测不能给出合格或超标的明确判断。可以通过多次检测,用平均值和平均值的标准差来进行判断,或采用精确度更高的检测方法给出确定性判断。采用S60方式,可以在尽可能短的时间,给出可靠地合格性判断。
其中f为置信度因子,95%置信度因子f=1.97,99%置信度因子f=2.65。
所述不可测物理量为用户需要了解的数据,其不可通过分析仪器检测获取。
S60、将用户能够接受的检测时间T0分成n段,以1/n*T0为检测时间段依次对不可测物理量值C2进行检测判断,如果当前检测时间段可以确定判断结果合格或不合格,则结束;否则,进行下一检测时间段继续进行检测,且将完成检测后的结果与前面的检测结果进行平均,根据平均结果判断合格或不合格,并依次类推,直至获得判断结果。
下面本发明以X射线荧光能谱仪对粮食中镉(Cd)污染快速检测为例,进行详细介绍说明。
首先采用Levenberg-Marquardt非线性拟合算法(LM算法),来得到Cd:Ka峰的强度,同时得到峰强度的标准差。
其中非线性拟合算法(LM算法)也称多变量非线性最小二乘拟合算法,需指定用来拟合实验数据的函数,函数中含有多个需确定的参数,通过LM拟合算法,得到使函数与实验数据误差最小的参数。LM算法的优点是可以同时得到参数的标准差。
之后仪器经过标准样品标定,得到Cd:Ka线的强度与Cd含量的线性关系C=a*I+b(C为Cd含量,I为Cd:Ka谱线的强度,a和b为最小二乘拟合得到的参数)。
接着检测未知样品,由于采用LM拟合算法计算Cd:Ka的强度的同时得到了Cd:Ka的强度的标准差σ1。通过C=a*I+b得到Cd的含量C,通过σ=a*σ1得到含量的标准差σ。
最后判断未知样品中Cd的含量:如果含量C小于C0-f*σ,则判定样品确定合格;如果含量C大于C0+f*σ,则判定样品超标。否则样品暂时不能判定否超标(其中C0为样品中Cd含量的限值;f为置信度因子,根据系统选择的置信度取值,95%置信度f=1.97,99%置信度f=2.65)。
如果Cd的含量介于C0-f*σ和C0+f*σ之间,则暂时不能判定样品中Cd的含量是否超标。此时可根据本发明所述的方法,采用分段检测来进一步判断。比如将用户能够接受的检测时间10分钟分成5段,以每2分钟为一个检测时间段,依次对Cd的含量进行检测判断,如果在前2分钟内可以确定判断结果合格或不合格,则结束;否则进行下一时间段2分钟的检测判断,将完成检测后的结果与前面的检测结果进行平均,根据平均结果判断合格或不合格,且并依次类推,直至获得判断结果。采用这种方式的好处是,对于含量显著低于限值或高于限值的结果,可以在第一时间段检测完成时,就可以给出合格或超标的准确判定以节省检测时间;而检测结果越接近限值,才需要更多次数的检测才能给出确定的判断。
如果判断结果是暂时不能判定是否超标,则需要进行下一段检测,完成后与前面的检测结果进行平均,标准差用平均标准差除以由于次数多了,标准差减小了,上一次不能判定是否超标,这一次就可能能够确定判断了。如果样品Cd含量很低或超标较多,一般一段测量就能得到确定的判断。含量越接近限量,需要的检测段数越多。
以上是对本发明所提供的基于限量管控应用的分析仪器测量结果自动判定方法进行了详细的介绍,本文中应用了具体个例对本发明的结构原理及实施方式进行了阐述,以上实施例只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。

Claims (3)

1.一种基于限量管控应用的分析仪器测量结果自动判定方法,其特征在于,包括步骤:
S10、设定样品中受限成分物理量的限值为C0
S20、获取样品的可测物理量值C1和该可测物理量值的标准差σ1
S30、通过标准样品对分析仪器进行标定,获取可测物理量值C1与不可测物理量值C2之间的线性关系公式,其中,所述标准样品为国家标准物质或国家标准样品或其物理量经过更高精密度的检测方法定值,所述线性关系公式为C2=a*C1+b,a和b为最小二乘拟合得到的参数;
S40、对未知样品进行检测,采用上述线性关系公式计算未知样品中不可测物理量值C2和该不可测物理量值的标准差σ2,其中σ2=a*σ1
S50、对未知样品中不可测物理量值C2进行判断,如果C2小于C0-f*σ2,则判定该未知样品合格;如果C2大于C0+f*σ2,则判定该未知样品不合格;如果C2介于C0-f*σ2与C0+f*σ2之间,则本次检测不能给出合格或超标的明确判断,则进入步骤S60;其中f为置信度因子,95%置信度因子f=1.97,99%置信度因子f=2.65;
S60、将用户能够接受的检测时间T0分成n段,以1/n*T0为检测时间段依次对不可测物理量值C2进行检测判断,如果当前检测时间段可以确定判断结果合格或不合格,则结束;否则,进行下一检测时间段继续进行检测,且将完成检测后的结果与前面的检测结果进行平均,根据平均结果判断合格或不合格,并依次类推,直至获得判断结果。
2.如权利要求1所述的基于限量管控应用的分析仪器测量结果自动判定方法,其特征在于,步骤S20包括:
通过标准差算法获取分析仪器测量样品结果的可测物理量C1和该可测物理量的标准差σ1
3.如权利要求2所述的基于限量管控应用的分析仪器测量结果自动判定方法,其特征在于,
所述标准差算法为多变量非线性最小二乘拟合算法。
CN201410620775.9A 2014-11-05 2014-11-05 基于限量管控应用的分析仪器测量结果自动判定方法 Expired - Fee Related CN104458783B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410620775.9A CN104458783B (zh) 2014-11-05 2014-11-05 基于限量管控应用的分析仪器测量结果自动判定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410620775.9A CN104458783B (zh) 2014-11-05 2014-11-05 基于限量管控应用的分析仪器测量结果自动判定方法

Publications (2)

Publication Number Publication Date
CN104458783A CN104458783A (zh) 2015-03-25
CN104458783B true CN104458783B (zh) 2017-08-29

Family

ID=52905183

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410620775.9A Expired - Fee Related CN104458783B (zh) 2014-11-05 2014-11-05 基于限量管控应用的分析仪器测量结果自动判定方法

Country Status (1)

Country Link
CN (1) CN104458783B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1640460A1 (en) * 2004-09-27 2006-03-29 Hitachi High-Technologies Corporation Method and apparatus for analyzing nucleic acid amplification
CN102680494A (zh) * 2012-05-24 2012-09-19 江南大学 基于机器视觉的抛光金属弧状面瑕疵实时检测方法
CN102735642A (zh) * 2012-07-12 2012-10-17 河北大学 一种快速无损鉴别初榨橄榄油和油橄榄果渣油的方法
CN103424523A (zh) * 2013-08-01 2013-12-04 浙江工商大学 黄豆品质检测装置及检测方法
CN103837312A (zh) * 2012-11-22 2014-06-04 赛龙通信技术(深圳)有限公司 气密性自动测试系统及测试方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1640460A1 (en) * 2004-09-27 2006-03-29 Hitachi High-Technologies Corporation Method and apparatus for analyzing nucleic acid amplification
CN102680494A (zh) * 2012-05-24 2012-09-19 江南大学 基于机器视觉的抛光金属弧状面瑕疵实时检测方法
CN102735642A (zh) * 2012-07-12 2012-10-17 河北大学 一种快速无损鉴别初榨橄榄油和油橄榄果渣油的方法
CN103837312A (zh) * 2012-11-22 2014-06-04 赛龙通信技术(深圳)有限公司 气密性自动测试系统及测试方法
CN103424523A (zh) * 2013-08-01 2013-12-04 浙江工商大学 黄豆品质检测装置及检测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
建筑结构强度无损检测区间评定方法;陈海彬 等;《应用基础与工程科学学报》;20050630;第13卷(第2期);第216-218页 *
食用海产品中亚硝酸盐含量的分析;张建云 等;《南京工业大学学报(自然科学版)》;20111130;第33卷(第6期);第107-110页 *

Also Published As

Publication number Publication date
CN104458783A (zh) 2015-03-25

Similar Documents

Publication Publication Date Title
CN108806218A (zh) 一种可燃气体监测数据异常原因的判断方法及装置
CN106441632A (zh) 一种温度计的检测方法及温度计的检测装置
CN1947003B (zh) 荧光x射线分析方法及荧光x射线分析装置
CN103675022A (zh) Tds 检测方法、检测装置和检测系统
CN104458630B (zh) 一种紫外差分气体分析仪的数据处理方法及系统
CN104502306A (zh) 基于变量重要性的近红外光谱波长选择方法
CN104458783B (zh) 基于限量管控应用的分析仪器测量结果自动判定方法
CN111985163B (zh) 一种铷钟寿命预测方法及装置
CN102726833A (zh) 一种卷烟机烟支空头探测仪器可靠性判断方法
CN107451408B (zh) 一种燃煤发热量数据校核方法
CN103575757A (zh) 荧光x射线分析方法和荧光x射线分析装置
CN105427029A (zh) 缺陷检测的校验系统及校验方法
CN114459523B (zh) 一种在线质量检测仪器的标定预警方法
CN111157938B (zh) 一种评价自动化检定系统计量过程能力的方法和系统
CN105719181A (zh) 风险等级评估方法及装置
CN112011125B (zh) 一种具有增塑剂迁移值的聚丙烯膜标准样品及其制备方法
CN114137471A (zh) 基于大数据分析的单相电能表检定流水线在线期间核查方法
CN107239256A (zh) 基于综合评价的彩票行业随机序列的随机性检测方法
CN109407042B (zh) 一种智能电表的校验方法
JP5751333B2 (ja) 管理装置、管理方法、プログラムおよび記録媒体
CN103257206A (zh) 有机磷和氨基甲酸酯类农药残留快速检测仪检测结果准确性的评价方法
CN104217765A (zh) 闪存芯片操作时间的测量方法
CN113790822B (zh) 地面测量温度数据异常检测方法、装置及可读存储介质
CN109580580B (zh) 一种物质成分检测的方法、装置和电子设备
CN107273656A (zh) 一种基于分子光谱分析的定量模型评估方法和系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CP02 Change in the address of a patent holder

Address after: 518000 Guangdong province Shenzhen city Nanshan District high tech park, North Long Hill Road No. 6 Yizhongli Science Park Building 2, two floor

Patentee after: SHENZHEN UNIQUE METRICAL TECHNOLOGY Co.,Ltd.

Address before: 518000 Guangdong city of Shenzhen province Nanshan District high tech Park North Lang Mountain Road No. 6 Yizhongli Science Park No. 2 building two floor

Patentee before: SHENZHEN UNIQUE METRICAL TECHNOLOGY Co.,Ltd.

CP02 Change in the address of a patent holder
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170829

CF01 Termination of patent right due to non-payment of annual fee