CN102262076B - 基于谱线组合的激光诱导击穿光谱元素浓度测量方法 - Google Patents

基于谱线组合的激光诱导击穿光谱元素浓度测量方法 Download PDF

Info

Publication number
CN102262076B
CN102262076B CN 201110210501 CN201110210501A CN102262076B CN 102262076 B CN102262076 B CN 102262076B CN 201110210501 CN201110210501 CN 201110210501 CN 201110210501 A CN201110210501 A CN 201110210501A CN 102262076 B CN102262076 B CN 102262076B
Authority
CN
China
Prior art keywords
line
spectral line
calibration
concentration
sample
Prior art date
Application number
CN 201110210501
Other languages
English (en)
Other versions
CN102262076A (zh
Inventor
王哲
李政
侯宗余
李立志
袁廷璧
Original Assignee
清华大学
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 清华大学 filed Critical 清华大学
Priority to CN 201110210501 priority Critical patent/CN102262076B/zh
Publication of CN102262076A publication Critical patent/CN102262076A/zh
Application granted granted Critical
Publication of CN102262076B publication Critical patent/CN102262076B/zh

Links

Abstract

基于谱线组合的激光诱导击穿光谱(LIBS)元素浓度测量方法,用在LIBS测量系统上。先用LIBS系统得到定标样品的光谱,选择目标元素的一条原子谱线和一条离子谱线,对这两条谱线的强度进行归一化然后加权得到组合强度,加权系数应使得定标样品的组合强度的波动最小。随后用单个或多个组合强度建立定标模型,即拟合出元素浓度和组合强度之间的函数关系。对于待测样品,先用LIBS系统得到其光谱,使用已求得的加权系数对相应的原子谱线和离子谱线的归一化后的强度进行加权得到待测样品的组合强度,然后代入定标模型中即可得到目标元素的浓度。该方法可降低等离子体参数的波动对LIBS元素浓度测量的影响并提高测量的精准度。

Description

基于谱线组合的激光诱导击穿光谱元素浓度测量方法
技术领域
[0001] 本发明涉及一种基于谱线组合的激光诱导击穿光谱元素浓度测量方法。方法的基本原理是在激光诱导等离子光谱(Laser Induced Breakdown Spectroscopy, LIBS)技术中,选择特定的一条原子谱线和一条离子谱线,对其谱线强度进行组合,得到组合强度,并利用组合强度进行定标和预测,从而降低LIBS元素浓度测量的不确定度,并提高定标优度和预测精度。
背景技术
[0002] LIBS技术的基本原理是将激光聚焦并击打在样品表面以形成高温、高电子密度的等离子体,然后用光谱仪记录等离子体发射的光谱信息用于分析样品的元素成分和浓度信息。该技术的主要优点有:1)几乎适用于各种样品(固、液、气);2)响应快,可用于实时测量;3)很少或者不用制备样品;4)对样品基本无破坏;5)能够实现全元素分析。这些优点使得激光诱导击穿技术广泛适用于很多领域。
[0003] 在LIBS元素浓度测量中,由于激光能量的波动以及激光-样品相互作用的波动,所产生的等离子体特性(形状、温度、电子密度等)也会产生波动,会导致测量的不确定度,这阻碍了 LIBS元素浓度测量精度的提高和商业化发展。因而降低信号的不确定度一直是LIBS研究的一个重要方向。
[0004]目前,传统的单变量模型都选择一条原子线的强度进行定标和预测,多变量模型,比如PLS,用到了多条谱线,但大部分仍是原子线或者未加选择地采用所有波长处的强度。大部分研究都选择原子线的强度进行定标和预测,而没有用到离子线所含的信息,主要原因是相对离子线来说,原子线强度波动较小,可得到较好的定标优度和预测精度。
[0005] 实际上,原子线和离子线强度的波动很大一部分是由等离子特性的波动造成的,而原子线和离子线强度随等离子特征参数的变化是明确的,而且原子线和离子线强度随等离子体温度和电子密度的变化有不同的响应,本文在分段归一化的基础上,根据原子线和离子线强度随温度和电子密度的不同响应,有针对性的对一对原子线和离子线进行组合,得到组合强度用于定标和预测,使得原子线和离子线强度随等离子温度和电子密度的不同响应在组合强度中相互抵消,大大降低了 LIBS元素浓度测量的不确定度,并且提高了定标优度和预测精度。
发明内容
[0006] 本发明的目的是针对目前的LIBS元素浓度测量中几乎只用到原子谱线强度,而对离子谱线却没有针对性的研究的缺陷,找到了一种利用原子谱线和离子谱线的强度进行组合的方法,从而通过组合强度来降低LIBS元素浓度测量的不确定度并提高测量的精准度。
[0007] 本发明的技术方案是:
[0008] 基于谱线组合的激光诱导击穿光谱元素浓度测量方法,其特征在于该方法包括如下步骤:
[0009] I)对于各元素浓度已知的一个定标样品,在其表面选定n个不同的位置,利用激光诱导击穿光谱测量系统在这n个不同位置进行检测,每个位置得到一幅包含各个元素原子和离子特征谱线的光谱;
[0010] 2)对于选定的不同的m个定标样品,重复步骤I),即得到每个定标样品在n个不同位置的包含各个元素原子和离子特征谱线的光谱,共mXn幅光谱;
[0011] 3)将待测量的元素定为目标元素,从定标样品的光谱中选择目标元素的一条原子特征谱线和一条离子特征谱线,求出各个定标样品各个位置的相应的原子特征谱线和离子特征谱线的强度,并对各个谱线强度进行归一化,归一化后的原子和离子谱线强度分别记为d,其中±标I代表原子、II代表离子,下标i,j代表第i个定标样品的第j个位置,其中i和j为正整数,i = 1,2,...,m,j = l,2,…,n;
[0012] 4)对/&和Zf7.进行加权得到组合强度Ii, j,/;J = a• I1uj + (1-a)!1^,a为加权系数,
确定加权系数a的方法为:在0到I之间变化a的值,直到组合强度Im随不同位置的波动最小;
[0013] 5)对于步骤3)中已选定的目标元素,重复步骤3)和4),即得到目标元素的不同的原子特征谱线和不同的离子特征谱线所组成的多种组合强度;
[0014] 6)用单个或多个组合强度建立定标模型,即用拟合的方法建立已选定的目标元素的浓度和组合强度之间的函数关系;
[0015] 7)对于目标元素浓度未知的待测样品,同样用激光诱导击穿光谱测量系统得到其单个测量位置的光谱或多个测量位置的平均光谱,求出相应的原子特征谱线强度和离子特征谱线强度,然后使用已求得的加权系数a对相应的归一化后的原子特征谱线强度/&和离子特征谱线强度进行加权,得到待测样品的组合强度Iq,把待测样品的组合强度
Ii, j代入定标模型中即得到待测样品中目标元素的浓度。
[0016] 本技术方案步骤3)中所述的谱线强度归一化方法采用全谱归一化或者分段归一化,即谱线强度除以全谱面积或特定波长范围内的面积。
[0017] 本技术方案步骤4)中所述的Ii, j的波动最小是指Ii, j的平均相对标准差达到最
小,即:
Figure CN102262076BD00041
达到最小,其中
Figure CN102262076BD00042
或者是指的最大相对标
准差达到最小,即
Figure CN102262076BD00043
达到最小。
[0018] 本技术方案步骤5)中所述的定标模型是单变量模型或多变量模型。
[0019] 本发明具有以下突出性特点:
[0020] 通过利用原子线强度和离子线强度随等离子体参数波动的不同响应,使用组合强度定标,在组合强度中部分抵消了等离子温度和电子密度对LIBS元素浓度测量的影响,从而大大降低了 LIBS元素浓度测量的不确定度,并提高了定标优度和测量精度。附图说明
[0021] 图1是本发明中LIBS测量系统的具体示意图。
[0022] 图2是本发明的技术方案示意图。
[0023] 图3是本发明的RSD降低效果图。
[0024] 图4是本发明的定标和预测的改善效果图。
具体实施方式
[0025] 下面结合附图对本发明作进一步的说明,但本发明并不止包含这一种实施方式。
[0026] 如图1、2所示。
[0027] 本发明提供的基于谱线组合的激光诱导击穿光谱元素浓度测量方法,该方法包括如下步骤:
[0028] I)对于各元素浓度已知的一个定标样品,在其表面选定n个不同的位置,利用激光诱导击穿光谱测量系统在这n个不同位置进行检测,每个位置得到一幅包含各个元素原子和离子特征谱线的光谱;
[0029] 2)对于选定的不同的m个定标样品,重复步骤I),即得到每个定标样品在n个不同位置的包含各个元素原子和离子特征谱线的光谱,共mXn幅光谱;
[0030] 3)将待测量的元素定为目标元素,从定标样品的光谱中选择目标元素的一条原子特征谱线和一条离子特征谱线,求出各个定标样品各个位置的相应的原子特征谱线和离子特征谱线的强度,并对各个谱线强度进行归一化,即谱线强度除以全谱面积或特定波长范围内的面积,归一化后的原子和离子谱线强度分别记为和,其中上标I代表原子、II代
表离子,下标i,j代表第i个定标样品的第j个位置,其中i和j为正整数,i = 1,2,…,m, j = 1,2, •••, n ;
[0031 ] 4)对I1u和g进行加权得到组合强度 a为加
权系数,确定加权系数a的方法为:在0到I之间变化a的值,直到组合强度
Figure CN102262076BD00051
不同位置的波动最小,Iq的波动最小是指Iq的平均相对标准差达到最小,即:
Figure CN102262076BD00052
达到最小,其中7 = iyj..或者是指Ii,j的最大相对标准差达到
'
最小,即
Figure CN102262076BD00053
达到最小;
[0032] 5)对于步骤3)中已选定的目标元素,重复步骤3)和4),即得到目标元素的不同的原子特征谱线和不同的离子特征谱线的多种组合强度;
[0033] 6)用单个或多个组合强度建立定标模型,即用拟合的方法建立已选定的目标元素的浓度和组合强度之间的函数关系,定标模型可以采用单变量定标模型或者多变量定标模型;
[0034] 7)对于目标元素浓度未知的待测样品,同样用激光诱导击穿光谱测量系统得到其单个测量位置的光谱或多个测量位置的平均光谱,求出相应的原子特征谱线强度和离子特征谱线强度,并使用已求得的加权系数a对相应的归一化后的原子谱线强度/&和离子谱线强度/&进行加权,得到待测样品的组合强度IM,把待测样品的组合强度Im代入定标
模型中即得到待测样品中目标元素的浓度。
[0035] 实施例:
[0036] I)使用各元素质量浓度已知的29种黄铜合金样品进行分析,其中20种作为定标样品,即m为20品,9种作为预测样品,各个样品的主要元素的质量浓度和如表I所示,此例中把铜作为目标元素。使用激光诱导击穿光谱测量系统对29种黄铜合金样品进行检测:如图1所示,以脉冲激光器I为激发光源,从激光器出射的激光经过聚焦透镜2聚焦后作用于黄铜样品3表面,在聚焦点产生等离子体,等离子体产生的辐射光信号通过采焦透镜4被实时收集,通过光纤5并经过光谱仪6处理后转化成电信号而被计算机7采集,得到各样品的光谱谱线,进一步得到各样品中铜元素的特征光谱谱线强度。其中每个样品在不同位置击打35个点,即n为35,每个点收集一条光谱,每条光谱都可求出特征光谱谱线强度,据此可求出每个样品35个点的特征光谱谱线强度均值和相对标准偏差(RSD)。
[0037] 表I黄铜合金样品的成分
[0038]
Figure CN102262076BD00061
[0039]
Figure CN102262076BD00071
[0040] *预测样品,其余为定标样品。
[0041] 2)选择铜原子261.837nm的特征谱线和铜离子201.69nm的特征谱线,选择这两条谱线的原因是其形状较规则,收到自吸收、互干扰的影响较小。对于20种定标样品,求出各个定标样品的两条特征谱线的强度并进行分段归一化,即各个强度除以该谱线所在的光谱仪波长范围内的面积,归一化后的原子线强度和离子线强度分别记为和^;其中上标I代表铜原子261.837nm.1I代表铜离子201.69mn,下标i,j代表第i个定标样品的第j个位置,其中i和j为正整数,i = 1,2,…,20,j = 1,2,…,35;
[0042] 3)对&和进行加权得到组合强度Ii, p/y =〃•/& ,a为加权系数,确定加权系数a的方法为:在0到I之间变化a的值,直到20个定标样品的组合强度的平均相对标准差达到最小,即
Figure CN102262076BD00081
达到最小,按照此方法,此例中求出a =
0.4250716 ;
[0043] 4)此例中的定标模型采用单变量线性定标模型,即通过线性拟合建立铜元素浓度和组合强度的直线关系,详见图4 ;
[0044] 5)对于9种预测样品中的每个样品,使用已求得的加权系数a = 0.4250716对相应的归一化后的原子谱线和离子谱线强度进行加权,得到待测样品的组合强度,归一化方法同样采用分段归一化法。把待测样品的组合强度代入4)中的线性定标模型中即可得到待测元素的浓度。
[0045] 图3对比了铜原子261.837nm的强度、铜离子201.69nm的强度以及两者组合强度的RSD,发现应用本发明提出的谱线组合方法可大大降低RSD。图4对比了 261.837nm的强度、铜离子201.69nm的强度以及两者组合强度的定标和预测效果,发现应用本发明的方法,也可提高定标优度和预测精度。表2总结了原子线强度、离子线强度和组合强度的定标和预测效果。
[0046] 表2本发明的改善效果
Figure CN102262076BD00082
[0048] 与上述方法相对应的在线检测设备(如图1所示),包括脉冲激光器(型号可以为Nd: YAG)、聚焦透镜2、光纤探头4、光纤5、光谱仪6、计算机7,脉冲激光器I发出的激光通过聚焦透镜2后,在样品3的表面聚焦,使得样品3被激光激发为等离子体,光纤探头4收集等离子体的发射光并通过光纤5输入光谱仪6,光谱仪6与计算机7相连,计算机7在相应的软件控制下可进行计算并实现显示或打印计算结果,同时提供相应的数据接口。

Claims (3)

1.基于谱线组合的激光诱导击穿光谱元素浓度测量方法,其特征在于该方法包括如下步骤: 1)对于各元素浓度已知的一个定标样品,在其表面选定n个不同的位置,利用激光诱导击穿光谱测量系统在这n个不同位置进行检测,每个位置得到一幅包含各个元素原子和离子特征谱线的光谱; 2)对于选定的不同的m个定标样品,重复步骤1),即得到每个定标样品在n个不同位置的包含各个元素原子和离子特征谱线的光谱,共mXn幅光谱; 3)将待测量的元素定为目标元素,从定标样品的光谱中选择目标元素的一条原子特征谱线和一条离子特征谱线,求出各个定标样品各个位置的相应的原子特征谱线和离子特征谱线的强度,并对各个谱线强度进行归一化,归一化后的原子和离子谱线强度分别记为和/&,其中上标I代表原子、II代表离子,下标i,j代表第i个定标样品的第j个位置,其中i和j为正整数,i = I, 2,…,m, j=l, 2,…,n ; 4)对进行加权得到组合强度Ii, j,
Figure CN102262076BC00021
a为加权系数,确定加权系数a的方法为:在O到I之间变化a的值,直到组合强度I。随不同位置的波动最小;Ii j的波动最小是指Ii ^的平均相对标准差达到最小,即
Figure CN102262076BC00022
达到最小, 其中
Figure CN102262076BC00023
或者是指的最大相对标准差达到最小,即
Figure CN102262076BC00024
达 到最小; 5)对于步骤3)中已选定的目标元素,重复步骤3)和4),即得到目标元素的不同的原子特征谱线和不同的离子特征谱线所组成的多种组合强度; 6)用单个或多个组合强度建立定标模型,即用拟合的方法建立已选定的目标元素的浓度和组合强度之间的函数关系; 7)对于目标元素浓度未知的待测样品,使用已求得的加权系数a对相应的归一化后的原子特征谱线强度和离子特征谱线强度进行加权,得到待测样品的组合强度Iu,把待测样品的组合强度Iu代入定标模型中即得到待测样品中目标元素的浓度。
2.按照权利要求1所述的基于谱线组合的激光诱导击穿光谱元素浓度测量方法,其特征在于:步骤3)中所述的谱线强度归一化方法采用全谱归一化或者分段归一化,即谱线强度除以全谱面积或特定波长范围内的面积。
3.按照权利要求1所述的基于谱线组合的激光诱导击穿光谱元素浓度测量方法,其特征还在于:步骤6)中所述的定标模型是单变量模型或多变量模型。
CN 201110210501 2011-07-26 2011-07-26 基于谱线组合的激光诱导击穿光谱元素浓度测量方法 CN102262076B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110210501 CN102262076B (zh) 2011-07-26 2011-07-26 基于谱线组合的激光诱导击穿光谱元素浓度测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110210501 CN102262076B (zh) 2011-07-26 2011-07-26 基于谱线组合的激光诱导击穿光谱元素浓度测量方法

Publications (2)

Publication Number Publication Date
CN102262076A CN102262076A (zh) 2011-11-30
CN102262076B true CN102262076B (zh) 2013-05-22

Family

ID=45008787

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110210501 CN102262076B (zh) 2011-07-26 2011-07-26 基于谱线组合的激光诱导击穿光谱元素浓度测量方法

Country Status (1)

Country Link
CN (1) CN102262076B (zh)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102798618B (zh) * 2012-07-16 2014-08-13 国电燃料有限公司 基于偏振降噪的脉冲激光粉状物质元素含量测量方法
CN103411931B (zh) * 2013-07-30 2015-07-29 中国科学院上海技术物理研究所 基于加权多谱线标定的远程libs元素定量分析方法
WO2015017623A1 (en) * 2013-08-02 2015-02-05 Tsi, Inc. High speed spectroscopic sensor assembly and system
CN104730041B (zh) * 2013-12-20 2017-03-22 武汉新瑞达激光工程有限责任公司 一种提高激光探针塑料识别精度的方法及其装置
CN103792214B (zh) * 2014-01-24 2016-07-13 国电科学技术研究院 一种提高钢铁中碳元素含量测量精度的方法
CN103792215B (zh) * 2014-01-24 2016-08-17 国电科学技术研究院 一种快速测量钢铁中碳元素含量的方法
CN104374752B (zh) * 2014-11-17 2017-01-11 浙江大学 基于共线激光诱导击穿光谱的作物营养元素快速检测方法
CN105092540B (zh) * 2015-06-16 2017-09-19 江西农业大学 一种食用油中重金属铅含量的快速高精度检测方法
CN105717093B (zh) * 2016-01-29 2018-08-28 清华大学 一种基于大数据库辨识的水泥特性分析方法
CN105718749B (zh) * 2016-01-29 2018-08-28 清华大学 一种基于大数据库辨识的煤质特性分析方法
CN105717094B (zh) * 2016-01-29 2018-04-10 清华大学 一种基于大数据库辨识的金属元素含量分析方法
CN105973886B (zh) * 2016-05-30 2019-03-29 田发益 一种基于可见光光谱确定浓度值的方法及装置
CN106442472B (zh) * 2016-12-09 2018-12-07 山西大学 一种基于等离子体温度精确测量的自由定标方法
CN107340284B (zh) * 2017-07-06 2018-04-13 湖北工程学院 一种元素定量分析方法及装置
CN107132214B (zh) * 2017-07-06 2018-03-27 湖北工程学院 一种基于多谱线加权的元素测量方法及装置
CN109813702A (zh) * 2019-01-21 2019-05-28 柳州职业技术学院 一种基于激光诱导击穿光谱技术定量分析茶叶中Cr元素的方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101509872A (zh) * 2009-03-20 2009-08-19 清华大学 一种基于回归分析的煤质在线检测分析方法
CN101915753A (zh) * 2010-07-30 2010-12-15 浙江师范大学 基于遗传神经网络的激光诱导击穿光谱定量分析方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8582089B2 (en) * 2006-06-09 2013-11-12 Chemimage Corporation System and method for combined raman, SWIR and LIBS detection

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101509872A (zh) * 2009-03-20 2009-08-19 清华大学 一种基于回归分析的煤质在线检测分析方法
CN101915753A (zh) * 2010-07-30 2010-12-15 浙江师范大学 基于遗传神经网络的激光诱导击穿光谱定量分析方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
James D.Winefordner,et al.Comparing several atomic spectrometric methods to the super stars: special emphasis on laser induced breakdown spectrometry, LIBS, a future super star.《Journal of Analytical Atomic Spectrometry》.2004,(第9期), *
提高激光诱导等离子光谱仪测量精度;王哲,等;《清华大学学报》;20100815;第50卷(第8期);第1263-1266页 *
王哲,等.提高激光诱导等离子光谱仪测量精度.《清华大学学报》.2010,第50卷(第8期),

Also Published As

Publication number Publication date
CN102262076A (zh) 2011-11-30

Similar Documents

Publication Publication Date Title
Wang et al. Laser-induced breakdown spectroscopy in China
Lazic et al. Self-absorption model in quantitative laser induced breakdown spectroscopy measurements on soils and sediments
Yaroshchyk et al. Comparison of principal components regression, partial least squares regression, multi-block partial least squares regression, and serial partial least squares regression algorithms for the analysis of Fe in iron ore using LIBS
Sirven et al. Qualitative and quantitative investigation of chromium-polluted soils by laser-induced breakdown spectroscopy combined with neural networks analysis
CN102128815B (zh) 检测时间及位置可控的激光诱导击穿光谱检测装置
Nicolodelli et al. Quantification of total carbon in soil using laser-induced breakdown spectroscopy: a method to correct interference lines
CN103534578B (zh) 确定样本中存在的成分的光谱装置以及方法
US9285272B2 (en) Dual source system and method
CN102879353B (zh) 近红外检测花生中蛋白质组分含量的方法
US7663749B2 (en) Method and system to measure the concentration of constituent elements in an inhomogeneous material using LIBS
CN102735677B (zh) 一种具有普适性的表面增强拉曼光谱定量分析方法
US20120033212A1 (en) Laser induced breakdown spectroscopy instrumentation for real-time elemental analysis
Li et al. A simplified spectrum standardization method for laser-induced breakdown spectroscopy measurements
CN104251846B (zh) 一种结合判别分析的激光诱导击穿光谱定量分析方法
CN201096521Y (zh) 非接触式等离子体温度和电子密度测量装置
Wang et al. A multivariate model based on dominant factor for laser-induced breakdown spectroscopy measurements
CN103234944B (zh) 一种基于主导因素结合偏最小二乘法的煤质特性分析方法
CN101655459B (zh) 一种激光诱导击穿光谱元素谱线逆向提取方法
Zhang et al. A novel approach for the quantitative analysis of multiple elements in steel based on laser-induced breakdown spectroscopy (LIBS) and random forest regression (RFR)
CN105388123B (zh) 一种由近红外光谱预测原油性质的方法
Golik et al. Determination of iron in water solution by time-resolved femtosecond laser-induced breakdown spectroscopy
JP5983779B2 (ja) ガス吸収分光装置及びガス吸収分光方法
CN101441165A (zh) 光谱仪的波长校准方法
CN101509872B (zh) 一种基于回归分析的煤质在线检测分析方法
CN102053083B (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
C14 Grant of patent or utility model
GR01 Patent grant