[go: up one dir, main page]

CN111693513A - 一种基于激光诱导击穿光谱的在线检测方法 - Google Patents

一种基于激光诱导击穿光谱的在线检测方法 Download PDF

Info

Publication number
CN111693513A
CN111693513A CN202010731010.8A CN202010731010A CN111693513A CN 111693513 A CN111693513 A CN 111693513A CN 202010731010 A CN202010731010 A CN 202010731010A CN 111693513 A CN111693513 A CN 111693513A
Authority
CN
China
Prior art keywords
tested
sample
spectral
point standard
prediction model
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
Application number
CN202010731010.8A
Other languages
English (en)
Other versions
CN111693513B (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.)
Aerospace Information Research Institute of CAS
Original Assignee
Aerospace Information Research Institute of CAS
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 Aerospace Information Research Institute of CAS filed Critical Aerospace Information Research Institute of CAS
Priority to CN202010731010.8A priority Critical patent/CN111693513B/zh
Publication of CN111693513A publication Critical patent/CN111693513A/zh
Application granted granted Critical
Publication of CN111693513B publication Critical patent/CN111693513B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/71Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light thermally excited
    • G01N21/718Laser microanalysis, i.e. with formation of sample plasma
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Algebra (AREA)
  • Plasma & Fusion (AREA)
  • Operations Research (AREA)
  • Evolutionary Biology (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Optics & Photonics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)

Abstract

本发明公开了一种基于激光诱导击穿光谱的在线检测方法,首先预先选择或制造一块单点标准样品;在相同检测条件下测定所述单点标准样品和待测样品的光谱信息;获得待测样品中待测元素的对应光谱强度与所述单点标准样品中对应元素的光谱强度比值R;针对需要检测的前五个待测样品,根据所述光谱强度比值R和预测模型获得待测样品中待测元素含量的预测值;从需要检测的第六个待测样品起,更新所述预测模型中的矩阵效应修正系数β和其他因素影响补偿Off,并利用更新后的预测模型对新的待测样品进行检测。该方法以尽可能低的模型建立成本,随着预测次数的增多不断提升预测水平,大大提升了LIBS的应用空间。

Description

一种基于激光诱导击穿光谱的在线检测方法
技术领域
本发明涉及激光诱导击穿光谱技术领域,尤其涉及一种基于激光诱导击穿光谱的在线检测方法。
背景技术
激光诱导击穿光谱(Laser-Induced Breakdown Spectroscopy,LIBS)是基于激光和材料相互作用产生的发射光谱的一种物质成分检测技术,与传统检测技术相比,LIBS优势之一在于可以通过光学系统实现复杂、恶劣环境下的隔离式(Stand-Off)在线物质成分分析。熔融金属在线检测是LIBS该优势的重要应用。在冶金生产过程中,通过LIBS进行在线检测能够缩短冶炼时间、避免废料回炉、大大提升多个批次精密合金之间的一致性,具有重大的经济价值。
在冶金生产中,熔融合金出炉冷却后,往往会经过实验室二次检测以确保合金成分符合生产标准。二次检测往往会采用精度优于LIBS在线检测的火花光谱法(SparkOptical Emission Spectrometry)或者电感耦合等离子体光谱法(Inductively CoupledPlasma Optical Emission Spectrometry)等手段。这些检测数据一方面可以对原LIBS定量模型进行一些补充,另一方面也一定程度上反应了检测环境带来的影响。将这些检测数据合理的利用起来,能够对LIBS在线检测的定量模型进行不断的自我校正,提升其下次检测的能力。
但现有技术往往不加甄选的将收到的二次检测数据添入原模型,重新回归建模,这种方法产生的模型往往会因为少数“较差”的二次检测数据,大大降低模型的预测能力;另外,LIBS冶金在线检测应用中还需考虑建立定量模型成本的问题,常用于在线检测的定量模型往往采用数块至数十块标准样品建立,这些标准样品在模型建立过程中,需要熔融后采集其光谱信息,标准样品价格昂贵,因此在保证模型性能的前提下,如何降低模型建立成本也是急需解决的问题。
发明内容
本发明的目的是提供一种基于激光诱导击穿光谱的在线检测方法,该方法以尽可能低的模型建立成本,随着预测次数的增多不断提升预测水平,大大提升了LIBS的应用空间。
本发明的目的是通过以下技术方案实现的:
一种基于激光诱导击穿光谱的在线检测方法,所述方法包括:
步骤1、预先选择或制造一块单点标准样品;
步骤2、在相同检测条件下测定所述单点标准样品和待测样品的光谱信息;
步骤3、获得所述待测样品中待测元素的对应光谱强度与所述单点标准样品中对应元素的光谱强度比值R;
步骤4、针对需要检测的前五个待测样品,根据所述光谱强度比值R和预测模型获得待测样品中待测元素含量的预测值;其中,所述预测模型表示为:
Csp=RCstd+β(R-1)+Off
其中,Csp为待测样品中待测元素含量的预测值;Cstd为所述单点标准样品中对应元素的含量;R为光谱强度比值;β为矩阵效应修正系数;Off为其他因素影响补偿;对于前五个待测样品,β和Off为0;
步骤5、从需要检测的第六个待测样品起,更新所述预测模型中的矩阵效应修正系数β和其他因素影响补偿Off,并利用更新后的预测模型对新的待测样品进行检测。
由上述本发明提供的技术方案可以看出,上述方法以尽可能低的模型建立成本,随着预测次数的增多不断提升预测水平,大大提升了LIBS的应用空间,为真正的工业级LIBS提供了参考。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域的普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他附图。
图1为本发明实施例提供的基于激光诱导击穿光谱的在线检测方法流程示意图。
具体实施方式
下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明的保护范围。
下面将结合附图对本发明实施例作进一步地详细描述,如图1所示为本发明实施例提供的基于激光诱导击穿光谱的在线检测方法流程示意图,所述方法包括:
步骤1、预先选择或制造一块单点标准样品;
在该步骤中,对于冶金行业而言,通过生产厂商提供的待检测产品信息,获得当前生产水平下该检测产品不同批次中各个元素含量的波动范围表;
以该波动范围表中各个元素对应的含量均值或生产厂商的目标元素配比值,选择或制造一块单点标准样品。
步骤2、在相同检测条件下测定所述单点标准样品和待测样品的光谱信息;
在该步骤中,标准样品的检测环境应当与待检测样品的检测环境尽可能保持一致,包括生产环境、仪器参数、检测状态;
所述单点标准样品检测得到的光谱信息为Sstd;待检测样品测得的光谱信息为Ssp[k],其中,k表示第k个待检测样品。
步骤3、获得所述待测样品中待测元素的对应光谱强度与所述单点标准样品中对应元素的光谱强度比值R;
在该步骤中,所述光谱强度比值R可以有多种表述形式,包括待测样品中待测元素与单点标准样品中对应元素的某条谱线强度或某几条谱线强度之和的比值;
或,待测样品中待测元素若干条谱线强度与单点标准样品中对应元素若干条谱线的强度拟合斜率比,具体公式为:
Figure BDA0002603344710000031
其中,
Figure BDA0002603344710000032
代表待测样品的待测元素的第i条谱线;
Figure BDA0002603344710000033
代表单点标准样品的对应元素的第i条谱线;linefitting代表基于最小二乘法的线性拟合;R代表拟合后得到的斜率,也即光谱强度比值。
具体实现中,谱线的强度可以采用峰值强度、谱线积分面积、谱线洛伦兹拟合后的积分面积等表述方式。
步骤4、针对需要检测的前五个待测样品,根据所述光谱强度比值R和预测模型获得待测样品中待测元素含量的预测值;
其中,所述预测模型表示为:
Csp=RCstd+β(R-1)+Off
其中,Csp为待测样品中待测元素含量的预测值;Cstd为所述单点标准样品中对应元素的含量;R为光谱强度比值;β为矩阵效应修正系数;Off为其他因素影响补偿;对于前五个待测样品,β和Off为0。
步骤5、从需要检测的第六个待测样品起,更新所述预测模型中的矩阵效应修正系数β和其他因素影响补偿Off,并利用更新后的预测模型对新的待测样品进行检测。
在该步骤中,预测模型的更新过程具体为:
首先以前五个待测样品的二次检测反馈值与所述预测模型预测值之间的差为因变量y,并以前五个待测样品对应的R-1的值为自变量x,假设y=βx+Off
采用最小二乘拟合得到β和Off的估计值
Figure BDA0002603344710000041
Figure BDA0002603344710000042
构建t统计量对x的系数进行检验,具体为:
Figure BDA0002603344710000043
其中,
Figure BDA0002603344710000044
Figure BDA0002603344710000045
Figure BDA0002603344710000046
xi和yi为第i个待测样品对应的x和y;
Figure BDA0002603344710000047
是i个x的均值;
Figure BDA0002603344710000048
是经过最小二乘拟合后得到的yi的估计值;
当|t0|<ta/2,n-2时,保持原模型中的β不变,更新Off,具体为:
Figure BDA0002603344710000049
其中,v代表权重,一般可选1/2;ta/2,n-2通过查询t统计量表得到,一般的a=0.05;
当|t0|>ta/2,n-2时,同时更新原模型中的β和Off,具体为:
Figure BDA0002603344710000051
Figure BDA0002603344710000052
其中,β0
Figure BDA0002603344710000053
表示原预测模型中的值;β和Off表示更新后的值。
然后再利用上述更新后的预测模型对新的待测样品进行检测,这样随着后续预测次数的增多,就可以不断提升预测水平,大大提升了LIBS的应用空间,为真正的工业级LIBS提供了参考。
值得注意的是,本发明实施例中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
需要说明的是,除了上述场景外,所有对LIBS在线检测数据进行了离线二次验证的应用,均可应用本文所提出的在线监测方法。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明披露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求书的保护范围为准。

Claims (5)

1.一种基于激光诱导击穿光谱的在线检测方法,其特征在于,所述方法包括:
步骤1、预先选择或制造一块单点标准样品;
步骤2、在相同检测条件下测定所述单点标准样品和待测样品的光谱信息;
步骤3、获得所述待测样品中待测元素的对应光谱强度与所述单点标准样品中对应元素的光谱强度比值R;
步骤4、针对需要检测的前五个待测样品,根据所述光谱强度比值R和预测模型获得待测样品中待测元素含量的预测值;其中,所述预测模型表示为:
Csp=RCstd+β(R-1)+Off
其中,Csp为待测样品中待测元素含量的预测值;Cstd为所述单点标准样品中对应元素的含量;R为光谱强度比值;β为矩阵效应修正系数;Off为其他因素影响补偿;对于前五个待测样品,β和Off为0;
步骤5、从需要检测的第六个待测样品起,更新所述预测模型中的矩阵效应修正系数β和其他因素影响补偿Off,并利用更新后的预测模型对新的待测样品进行检测。
2.根据权利要求1所述基于激光诱导击穿光谱的在线检测方法,其特征在于,在步骤1中,具体通过生产厂商提供的待检测产品信息,获得当前生产水平下该检测产品不同批次中各个元素含量的波动范围表;
以该波动范围表中各个元素对应的含量均值或生产厂商的目标元素配比值,选择或制造一块单点标准样品。
3.根据权利要求1所述基于激光诱导击穿光谱的在线检测方法,其特征在于,在步骤2中,所述相同检测条件包括生产环境、仪器参数、检测状态;
所述单点标准样品检测得到的光谱信息为Sstd;待检测样品测得的光谱信息为Ssp[k],其中,k表示第k个待检测样品。
4.根据权利要求1所述基于激光诱导击穿光谱的在线检测方法,其特征在于,在步骤1中,在步骤3中,所述光谱强度比值R的形式包括待测样品中待测元素与单点标准样品中对应元素的某条谱线强度或某几条谱线强度之和的比值;
或,待测样品中待测元素若干条谱线强度与单点标准样品中对应元素若干条谱线的强度拟合斜率比,具体公式为:
Figure FDA0002603344700000021
其中,
Figure FDA0002603344700000022
代表待测样品的待测元素的第i条谱线;
Figure FDA0002603344700000023
代表单点标准样品的对应元素的第i条谱线;linefitting代表基于最小二乘法的线性拟合;R代表拟合后得到的斜率,也即光谱强度比值。
5.根据权利要求1所述基于激光诱导击穿光谱的在线检测方法,其特征在于,在步骤5中,预测模型的更新过程具体为:
首先以前五个待测样品的二次检测反馈值与所述预测模型预测值之间的差为因变量y,并以前五个待测样品对应的R-1的值为自变量x,假设y=βx+Off
采用最小二乘拟合得到β和Off的估计值
Figure FDA0002603344700000024
Figure FDA0002603344700000025
构建t统计量对x的系数进行0检验,具体为:
Figure FDA0002603344700000026
其中,
Figure FDA0002603344700000027
Figure FDA0002603344700000028
Figure FDA00026033447000000212
xi和yi为第i个待测样品对应的x和y;
Figure FDA0002603344700000029
是i个x的均值;
Figure FDA00026033447000000210
是经过最小二乘拟合后得到的yi的估计值;
当|t0|<ta/2,n-2时,保持原模型中的β不变,更新Off,具体为:
β=0,
Figure FDA00026033447000000211
其中,v代表权重;ta/2,n-2通过查询t统计量表得到;
当|t0|>ta/2,n-2时,同时更新原模型中的β和Off,具体为:
Figure FDA0002603344700000031
Figure FDA0002603344700000032
其中,β0
Figure FDA0002603344700000033
表示原预测模型中的值;β和Off表示更新后的值。
CN202010731010.8A 2020-07-27 2020-07-27 一种基于激光诱导击穿光谱的在线检测方法 Expired - Fee Related CN111693513B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010731010.8A CN111693513B (zh) 2020-07-27 2020-07-27 一种基于激光诱导击穿光谱的在线检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010731010.8A CN111693513B (zh) 2020-07-27 2020-07-27 一种基于激光诱导击穿光谱的在线检测方法

Publications (2)

Publication Number Publication Date
CN111693513A true CN111693513A (zh) 2020-09-22
CN111693513B CN111693513B (zh) 2023-03-31

Family

ID=72486741

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010731010.8A Expired - Fee Related CN111693513B (zh) 2020-07-27 2020-07-27 一种基于激光诱导击穿光谱的在线检测方法

Country Status (1)

Country Link
CN (1) CN111693513B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114739980A (zh) * 2022-06-10 2022-07-12 合肥金星智控科技股份有限公司 元素信息的预测方法、装置、设备及介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103234944A (zh) * 2013-04-17 2013-08-07 清华大学 一种基于主导因素结合偏最小二乘法的煤质特性分析方法
CN103792214A (zh) * 2014-01-24 2014-05-14 国电科学技术研究院 一种提高钢铁中碳元素含量测量精度的方法
US20190079019A1 (en) * 2017-08-15 2019-03-14 Pablo Sobron Non-linear methods for quantitative elemental analysis and mineral classification using laser-induced breakdown spectroscopy (LIBS)
CN109781711A (zh) * 2019-02-21 2019-05-21 华中科技大学 一种基于单标样定标的激光诱导击穿光谱定量分析方法
CN109799195A (zh) * 2019-01-22 2019-05-24 上海交通大学 一种激光诱导击穿光谱的高精度定量分析方法
CN110018294A (zh) * 2019-05-14 2019-07-16 北京农业信息技术研究中心 土壤重金属检测值校正方法、装置和计算机存储介质

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103234944A (zh) * 2013-04-17 2013-08-07 清华大学 一种基于主导因素结合偏最小二乘法的煤质特性分析方法
CN103792214A (zh) * 2014-01-24 2014-05-14 国电科学技术研究院 一种提高钢铁中碳元素含量测量精度的方法
US20190079019A1 (en) * 2017-08-15 2019-03-14 Pablo Sobron Non-linear methods for quantitative elemental analysis and mineral classification using laser-induced breakdown spectroscopy (LIBS)
CN109799195A (zh) * 2019-01-22 2019-05-24 上海交通大学 一种激光诱导击穿光谱的高精度定量分析方法
CN109781711A (zh) * 2019-02-21 2019-05-21 华中科技大学 一种基于单标样定标的激光诱导击穿光谱定量分析方法
CN110018294A (zh) * 2019-05-14 2019-07-16 北京农业信息技术研究中心 土壤重金属检测值校正方法、装置和计算机存储介质

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114739980A (zh) * 2022-06-10 2022-07-12 合肥金星智控科技股份有限公司 元素信息的预测方法、装置、设备及介质
WO2023236387A1 (zh) * 2022-06-10 2023-12-14 合肥金星智控科技股份有限公司 元素信息的预测方法、装置、设备及介质

Also Published As

Publication number Publication date
CN111693513B (zh) 2023-03-31

Similar Documents

Publication Publication Date Title
CN107703097B (zh) 利用近红外光谱仪构建快速预测原油性质的模型的方法
CN104697966A (zh) 一种基于最小二乘支持向量机算法结合激光诱导击穿光谱的钢铁中铬和锰定量分析方法
WO2018184262A1 (zh) 用于激光诱导击穿光谱采集的阶梯光谱仪动态校正方法
CN103808695B (zh) 一种基于激光诱导击穿光谱技术检测铁矿石全铁的方法
CN105241866B (zh) 利用荧光褪色效应消除拉曼光谱中基线干扰的方法
CN105717094B (zh) 一种基于大数据库辨识的金属元素含量分析方法
CN113588597A (zh) 一种提升炉渣分析精度的方法
CN104569019A (zh) 铜阳极泥中金属元素含量的测定方法
CN111693513A (zh) 一种基于激光诱导击穿光谱的在线检测方法
CN104730043A (zh) 一种基于偏最小二乘法的墨水中重金属测定方法
CN107941758B (zh) 一种基于光谱预处理的煤质热值测量方法
CN119104542B (zh) 一种基于标靶校准的煤质检测方法、装置、设备及介质
CN118225752B (zh) 一种煤质检测方法、装置、电子设备和存储介质
CN101105456A (zh) 一种合金钢样品夹杂含量快速测定的分析方法
CN114636687A (zh) 基于深度迁移学习的小样本煤质特性分析系统及方法
CN117990615A (zh) 一种实时灵敏直测染料中间体浓度的色度学方法
CN105717093B (zh) 一种基于大数据库辨识的水泥特性分析方法
CN113155879A (zh) 一种测定萤石中二氧化硅和氟化钙含量的方法
CN106841175A (zh) 一种炼钢精炼炉的钢水成分分析系统
CN108267422B (zh) 基于近红外光谱分析的异常样本剔除法
CN116482079A (zh) 基于激光诱导击穿光谱的检测方法及系统
CN116858870A (zh) 一种铝锌合金中成分含量的检测方法
CN107328760B (zh) 一种针对合金样品表面蒸汽对lips测量改进的方法
CN103472007A (zh) 一种合金中元素的检测方法
CN111189866A (zh) 一种x荧光光谱仪监控直读光谱仪的检测方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20230331