CN106442472A - 一种基于等离子体温度精确测量的自由定标方法 - Google Patents

一种基于等离子体温度精确测量的自由定标方法 Download PDF

Info

Publication number
CN106442472A
CN106442472A CN201611132076.5A CN201611132076A CN106442472A CN 106442472 A CN106442472 A CN 106442472A CN 201611132076 A CN201611132076 A CN 201611132076A CN 106442472 A CN106442472 A CN 106442472A
Authority
CN
China
Prior art keywords
elements
plasma
temperature
content
plasma temperature
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
CN201611132076.5A
Other languages
English (en)
Other versions
CN106442472B (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.)
Shanxi University
Original Assignee
Shanxi University
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 Shanxi University filed Critical Shanxi University
Priority to CN201611132076.5A priority Critical patent/CN106442472B/zh
Publication of CN106442472A publication Critical patent/CN106442472A/zh
Application granted granted Critical
Publication of CN106442472B publication Critical patent/CN106442472B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/71—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light thermally excited
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K15/00—Testing or calibrating of thermometers
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/01—Arrangements or apparatus for facilitating the optical investigation

Landscapes

  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)

Abstract

本发明属于激光光谱分析与检测方法技术领域，具体涉及一种基于等离子体温度精确测量的自由定标方法。本方法主要是通过对比不同等离子体温度时已知样品中所有元素含量的自由定标预测值与实际值间的平均相对误差，选择最小平均相对误差所对应的温度值作为最佳等离子体温度。利用该最佳等离子体温度并结合自由定标法，可以实现同种激光诱导击穿光谱检测条件下相同基体未知样品中所有元素含量的精确测量，解决了前自由定标激光诱导击穿光谱定量分析过程中等离子体温度测量精度不够高而导致元素含量测量误差大的技术问题。本方法较于传统方法的优势在于提高了离子体温度的测量精度，从而实现了元素含量精确测量的目的。

Description

一种基于等离子体温度精确测量的自由定标方法
技术领域
本发明属于激光光谱分析与检测方法技术领域，具体涉及一种基于等离子体温度精确测量的自由定标方法。
背景技术
自由定标激光诱导击穿光谱属于一种激光烧蚀光谱分析技术，具有无需标准样品定标、可多元素同时分析、分析速度快、安全无辐射等优势，可广泛应用于勘矿、冶金、制药、化工、考古、珠宝鉴定、星际探测等领域。该技术的定量分析精度主要取决于对激光诱导等离子体温度的测定精度。传统的等离子体温度测定方法有Boltzmann平面法和Saha双线法，前者主要利用同一元素多条原子谱线或多条离子谱线来计算温度，后者主要利用同一元素的单条原子谱线和单条离子谱线来计算温度。利用以上方法，通过若干元素谱线计算出若干温度后将其均值作为等离子体温度。然而，该温度均值与等离子体实际温度仍有较大差异，并不能全面反映等离子体温度情况，若将其用于自由定标激光诱导击穿光谱的定量分析计算中，必然会导致各元素含量测量结果的误差较大。例如，意大利科学院等离子体研究所利用自由定标激光诱导击穿光谱技术对陨石进行了分析，其中采用基于25条Fe原子谱线的Boltzmann平面对等离子体温度进行了测量，结合该温度值对陨石中Fe、Ni、Co三种元素含量的相对测量误差分别为±14％、±11％、±20％；白俄罗斯明斯克原子分子物理研究所利用自由定标激光诱导击穿光谱技术对铜合金进行了分析，其中采用基于8条Cu原子谱线和4条Pb原子谱线的Boltzmann平面分别对等离子体温度进行了测量并求均值，结合该温度均值对铜合金中Zn、Pb、Sn三种元素含量的相对测量误差分别为±15％、±21％、±11％；奥地利约翰开普勒大学利用自由定标激光诱导击穿光谱技术对钢渣进行了分析，其中采用基于13条Ca原子谱线、16条Mn原子谱线、7条Mn离子谱线、6条Ti离子谱线的Boltzmann平面分别对等离子体温度进行了测量并求均值，结合该温度均值对Ca、Al、Mg、Si、Fe、Mn、Ti七种元素含量的相对测量误差分别为±6％、±11％、±10％、±21％、±32％、±18％、±23％。可见，基于传统等离子体温度测定方法的自由定标激光诱导击穿光谱技术对元素含量测量的相对误差为10～20％。
发明内容
本发明的目的是解决目前自由定标激光诱导击穿光谱定量分析过程中等离子体温度测量精度不够高而导致元素含量测量误差大的技术问题，提供一种基于等离子体温度精确测量的自由定标方法。
为解决上述技术问题，本发明采用的技术方案为：
一种基于等离子体温度精确测量的自由定标方法，包括以下步骤：
(1)选出已知标准样品中所有元素的原子谱线，利用本领域技术人员所公知的Boltzmann平面法分别计算出等离子体温度Ti，并给出Ti的变化范围，其中i代表元素种类；
(2)由一元三阶方程拟合求出所有元素的原子和离子配分函数，其中原子配分函数记为离子配分函数记为
(3)以ΔT＝100K为步长求出不同温度T时所有元素的相对原子数密度
其中，Min(Ti)≤T≤Max(Ti)，g为简并度，A为自发辐射系数，h为普朗克常量，v为发射光子频率，Em为上能级能量，k为Boltzmann常数；
(4)利用Stark展宽求出等离子体电子密度Ne：
其中，Δλ为谱线线宽，ω为谱线的电子碰撞参数；
(5)求出所有元素的相对离子数密度
其中，me为电子质量，En为下能级能量；
(6)求出不同温度T时所有元素含量的自由定标预测值Ci(T)：
其中，M为元素相对原子质量；
(7)求出Ci(T)的相对误差ΔCi(T)：
其中，C*为标准样品中元素真实含量；
(8)求出ΔCi(T)的归一化值
(9)求出的均值
(10)将最小时所对应的温度Topt作为最佳等离子体温度；
(11)将Topt作为同种激光诱导击穿光谱检测条件下相同基体待测样品的等离子体温度，利用本领域技术人员所公知的自由定标方法，即可实现对该待测样品中各元素含量的精确测量。
本发明采用以上技术方案，通过对比不同等离子体温度时已知样品中所有元素含量的自由定标预测值与实际值间的平均相对误差，选择最小平均相对误差所对应的温度值作为最佳等离子体温度。利用该最佳等离子体温度并结合自由定标法，可以实现同种激光诱导击穿光谱检测条件下相同基体未知样品中所有元素含量的精确测量，解决了前自由定标激光诱导击穿光谱定量分析过程中等离子体温度测量精度不够高而导致元素含量测量误差大的技术问题。
附图说明
图1是本发明所述自由定标激光诱导击穿光谱定量分析装置结构示意图；
图中，1-Nd:YAG脉冲激光器，2-镀铝反射镜，3-平凸透镜，4-样品，5-电控旋转台，6-耦合透镜，7-光纤，8-光栅光谱仪，9-ICCD，10-计算机；
图2是水泥样品等离子体光谱，其中标注了Si、Al、Ca、Fe、Mg谱线；
图3是Al元素谱线计算等离子体温度的Boltzmann平面；
图4是Ca元素相对原子数密度的数值模拟；
图5是Ca元素相对离子数密度的数值模拟；
图6是Ca元素含量的自由定标预测值CCa(T)；
图7是CCa(T)的相对误差ΔCCa(T)；
图8是ΔCCa(T)的归一化值
图9是水泥样品中Si、Al、Ca、Fe、Mg所对应的均值
具体实施方式
本实施例以水泥样品为例，结合附图对本发明所述的一种基于等离子体温度精确测量的自由定标方法进一步描述，包括以下步骤：
(1)选出标准样品中所有元素的原子谱线，利用Boltzmann平面法分别计算出等离子体温度Ti，并给出Ti的变化范围，其中i代表元素种类；选出已知标准水泥样品等离子体光谱中元素Si、Al、Ca、Fe、Mg的原子谱线(表1)，利用本领域技术人员所公知的Boltzmann平面法分别计算出等离子体温度Ti分别为：TSi＝5560K、TAL＝8360K、TCa＝12530K、TFe＝9540K、TMg＝1360K，则Ti的范围为5000K≤Ti≤15000K；图3为利用Al元素谱线计算等离子体温度的Boltzmann平面；
表1水泥样品中元素Si、Al、Ca、Fe、Mg的原子谱线及参数
(2)由一元三阶方程拟合求出所有元素的原子和离子配分函数，其中原子配分函数记为离子配分函数记为对水泥中Si、Al、Ca、Fe、Mg的原子和离子配分函数进行一元三阶阶方程拟合，其中原子配分函数分别为：
离子配分函数分别为：
(3)以ΔT＝100K为步长求出不同温度T时所有元素的相对原子数密度
其中，Min(Ti)≤T≤Max(Ti)，g为简并度，A为自发辐射系数，h为普朗克常量，v为发射光子频率，Em为上能级能量，k为Boltzmann常数；
以ΔT＝100K为步长求出5000-15000K范围内Si、Al、Ca、Fe、Mg的相对原子数密度例如，Ca元素的相对原子数密度的计算结果如图4所示；
(4)利用Stark展宽求出等离子体电子密度Ne：
其中，Δλ为谱线线宽，ω为谱线的电子碰撞参数；
利用Mg 277.827nm谱线的Stark展宽求出等离子体电子密度Ne＝1.4558×1018cm-3；
(5)求出所有元素的相对离子数密度
其中，me为电子质量，En为下能级能量；
求出Si、Al、Ca、Fe、Mg相对原子数密度所对应的相对离子数密度例如，Ca元素的相对离子数密度的数值模拟如图5所示；
(6)求出不同温度T时所有元素含量的自由定标预测值Ci(T)：
其中，M为元素相对原子质量；
求出不同温度T时Si、Al、Ca、Fe、Mg元素含量的自由定标预测值Ci(T)。例如，Ca元素含量的自由定标预测值CCa(T)如图6所示；
(7)求出Ci(T)的相对误差ΔCi(T)：
其中，C*为标准样品中元素真实含量；
求出Ci(T)的相对误差ΔCi(T)；例如，Ca元素含量预测的相对误差ΔCCa(T)如图7所示；
(8)求出ΔCi(T)的归一化值
例如，Ca元素含量预测的相对误差归一化值如图8所示；
(9)求出的均值
求出的均值结果如图9所示；
(10)将最小时所对应的温度Topt＝12700K作为最佳等离子体温度；
(11)将Topt＝12700K作为同种激光诱导击穿光谱检测条件下未知水泥样品的等离子体温度，利用自由定标方法，即可实现对该未知样品中Si、Al、Ca、Fe、Mg元素含量的相对测量误差分别为±5.8％、±12％、±0.07％、±9.4％、±12.6％。
如图1所示，本发明基于等离子体温度精确测量的自由定标方法所采用的装置Nd:YAG脉冲激光器1出射的高能窄脉冲激光经镀铝反射镜2后传播方向由水平变为竖直，再经过平凸透镜3后聚焦在样品4表面，诱导产生等离子体，等离子体荧光经过耦合透镜6汇聚至光纤7中传输至光栅光谱仪8，光栅光谱仪8对等离子体荧光按波长进行分光后投射到ICCD9上进行光电转换和模数转换，转换后的相应光谱数据传输至计算机10中进行处理。其中，电控旋转台5用于连续旋转样品4，以确保每个激光脉冲都有新的作用点。

Claims (1)

1.一种基于等离子体温度精确测量的自由定标方法，其特征在于：包括以下步骤：
(1)选出标准样品中所有元素的原子谱线，利用Boltzmann平面法分别计算出等离子体温度Ti，并给出Ti的变化范围，其中i代表元素种类；
(2)由一元三阶方程拟合求出所有元素的原子和离子配分函数，其中原子配分函数记为离子配分函数记为
(3)以ΔT＝100K为步长求出不同温度T时所有元素的相对原子数密度
N i 0 ( T ) = I i j Z i 0 ( T ) g i j A i j hv i j exp ( E m i j k T ) ,
其中，Min(Ti)≤T≤Max(Ti)，g为简并度，A为自发辐射系数，h为普朗克常量，v为发射光子频率，Em为上能级能量，k为Boltzmann常数；
(4)利用Stark展宽求出等离子体电子密度Ne：
N e = 5 × 10 15 Δλ i j ω i j ,
其中，Δλ为谱线线宽，ω为谱线的电子碰撞参数；
(5)求出所有元素的相对离子数密度
N i 1 ( T ) = 2 N i 0 ( T ) Z i 1 ( T ) N e Z i 0 ( T ) ( 2 πm e k T h 2 ) 3 / 2 exp ( - E n i j K T ) ,
其中，me为电子质量，En为下能级能量；
(6)求出不同温度T时所有元素含量的自由定标预测值Ci(T)：
C i ( T ) = ( N i 0 ( T ) + N i 1 ( T ) ) M i Σ i ( N i 0 ( T ) + N i 1 ( T ) ) M i ,
其中，M为元素相对原子质量；
(7)求出Ci(T)的相对误差ΔCi(T)：
ΔC i ( T ) = | C i * - C i ( T ) | C i * ,
其中，C*为标准样品中元素真实含量；
(8)求出ΔCi(T)的归一化值
ΔC i ( T ) ‾ = ΔC i ( T ) M a x ( ΔC i ( T ) ) ;
(9)求出的均值
Δ C ( T ) ‾ ‾ = Σ i ΔC i ( T ) ‾ i ;
(10)将最小时所对应的温度Topt作为最佳等离子体温度；
(11)将Topt作为同种激光诱导击穿光谱检测条件下相同基体待测样品的等离子体温度，利用自由定标方法，即可实现对该待测样品中各元素含量的精确测量。
CN201611132076.5A 2016-12-09 2016-12-09 一种基于等离子体温度精确测量的自由定标方法 Expired - Fee Related CN106442472B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611132076.5A CN106442472B (zh) 2016-12-09 2016-12-09 一种基于等离子体温度精确测量的自由定标方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611132076.5A CN106442472B (zh) 2016-12-09 2016-12-09 一种基于等离子体温度精确测量的自由定标方法

Publications (2)

Publication Number Publication Date
CN106442472A true CN106442472A (zh) 2017-02-22
CN106442472B CN106442472B (zh) 2018-12-07

Family

ID=58216331

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611132076.5A Expired - Fee Related CN106442472B (zh) 2016-12-09 2016-12-09 一种基于等离子体温度精确测量的自由定标方法

Country Status (1)

Country Link
CN (1) CN106442472B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108572168A (zh) * 2018-04-24 2018-09-25 山西大学 基于谱线自吸收量化的激光诱导等离子体参数表征方法
CN108572177A (zh) * 2018-05-21 2018-09-25 中北大学 基于发射光谱的焊接电弧污染物在线定温检测方法
CN109932355A (zh) * 2017-12-15 2019-06-25 中石化石油工程技术服务有限公司 钻井液地层岩性-流体激光信息的分离方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100309604A1 (en) * 2005-08-05 2010-12-09 Lam Research Corporation Method and Apparatus for Chuck Thermal Calibration
CN101949852A (zh) * 2010-07-30 2011-01-19 清华大学 一种基于光谱标准化的煤质在线检测方法
CN102262076A (zh) * 2011-07-26 2011-11-30 清华大学 基于谱线组合的激光诱导击穿光谱元素浓度测量方法
CN102313731A (zh) * 2010-07-09 2012-01-11 中国科学院沈阳自动化研究所 一种未知物组成元素含量在线检测方法
CN102410993A (zh) * 2011-08-01 2012-04-11 清华大学 基于激光诱导等离子体发射光谱标准化的元素测量方法
CN102410992A (zh) * 2011-08-01 2012-04-11 清华大学 简化的激光诱导等离子体光谱标准化的元素测量方法
CN103983353A (zh) * 2014-03-22 2014-08-13 中国科学技术大学 一种利用等离子体发射光谱实现光学系统传输效率的标定方法
CN105203526A (zh) * 2015-09-22 2015-12-30 中国科学院上海技术物理研究所 免定标的远程定量激光诱导击穿光谱分析方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100309604A1 (en) * 2005-08-05 2010-12-09 Lam Research Corporation Method and Apparatus for Chuck Thermal Calibration
CN102313731A (zh) * 2010-07-09 2012-01-11 中国科学院沈阳自动化研究所 一种未知物组成元素含量在线检测方法
CN101949852A (zh) * 2010-07-30 2011-01-19 清华大学 一种基于光谱标准化的煤质在线检测方法
CN102262076A (zh) * 2011-07-26 2011-11-30 清华大学 基于谱线组合的激光诱导击穿光谱元素浓度测量方法
CN102410993A (zh) * 2011-08-01 2012-04-11 清华大学 基于激光诱导等离子体发射光谱标准化的元素测量方法
CN102410992A (zh) * 2011-08-01 2012-04-11 清华大学 简化的激光诱导等离子体光谱标准化的元素测量方法
CN103983353A (zh) * 2014-03-22 2014-08-13 中国科学技术大学 一种利用等离子体发射光谱实现光学系统传输效率的标定方法
CN105203526A (zh) * 2015-09-22 2015-12-30 中国科学院上海技术物理研究所 免定标的远程定量激光诱导击穿光谱分析方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
ILARIABORGIA等: "Self-calibrated quantitative elemental analysis by laser-induced plasma spectroscopy: application to pigment analysis", 《JOURNAL OF CULTURAL HERITAGE》 *
张雷等: "激光诱导等离子体光谱用于煤中氧的定量分析", 《量子光学学报》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109932355A (zh) * 2017-12-15 2019-06-25 中石化石油工程技术服务有限公司 钻井液地层岩性-流体激光信息的分离方法
CN109932355B (zh) * 2017-12-15 2024-01-26 中国石油化工集团有限公司 钻井液地层岩性-流体激光信息的分离方法
CN108572168A (zh) * 2018-04-24 2018-09-25 山西大学 基于谱线自吸收量化的激光诱导等离子体参数表征方法
CN108572177A (zh) * 2018-05-21 2018-09-25 中北大学 基于发射光谱的焊接电弧污染物在线定温检测方法

Also Published As

Publication number Publication date
CN106442472B (zh) 2018-12-07

Similar Documents

Publication Publication Date Title
Nicolodelli et al. Recent advances and future trends in LIBS applications to agricultural materials and their food derivatives: An overview of developments in the last decade (2010–2019). Part I. Soils and fertilizers
Jolivet et al. Review of the recent advances and applications of LIBS-based imaging
CN102313731B (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)
Li et al. Analytical-performance improvement of laser-induced breakdown spectroscopy for steel using multi-spectral-line calibration with an artificial neural network
CN104251846B (zh) 一种结合判别分析的激光诱导击穿光谱定量分析方法
Cavalcanti et al. One-point calibration for calibration-free laser-induced breakdown spectroscopy quantitative analysis
Andrade et al. Classical univariate calibration and partial least squares for quantitative analysis of brass samples by laser-induced breakdown spectroscopy
Khater Laser-induced breakdown spectroscopy for light elements detection in steel: State of the art
Zhang et al. A plasma-image-assisted method for matrix effect correction in laser-induced breakdown spectroscopy
CN102262076B (zh) 基于谱线组合的激光诱导击穿光谱元素浓度测量方法
Sarkar et al. Evaluation of the prediction precision capability of partial least squares regression approach for analysis of high alloy steel by laser induced breakdown spectroscopy
CN102410993B (zh) 基于激光诱导等离子体发射光谱标准化的元素测量方法
US10024802B2 (en) Method for laser-induced breakdown spectroscopy and calibration
CN104697966A (zh) 一种基于最小二乘支持向量机算法结合激光诱导击穿光谱的钢铁中铬和锰定量分析方法
CN102749309B (zh) 冶金成分快速分析的激光诱导击穿光谱仪和分析方法
CN106442472A (zh) 一种基于等离子体温度精确测量的自由定标方法
CN112834485B (zh) 激光诱导击穿光谱元素定量分析的一种非定标方法
Pagnotta et al. Comparison of brass alloys composition by laser-induced breakdown spectroscopy and self-organizing maps
CN102410992B (zh) 简化的激光诱导等离子体光谱标准化的元素测量方法
CN104897626B (zh) 火星主次元素同步分析的定量激光诱导击穿光谱方法
Yin et al. Cement raw material quality analysis using laser-induced breakdown spectroscopy
Bonta et al. Metal analysis in polymers using tandem LA-ICP-MS/LIBS: eliminating matrix effects using multivariate calibration
CN106442474B (zh) 一种基于偏最小二乘法的水泥生料三率值测量方法
Yu et al. Response surface methodology for optimizing LIBS testing parameters: A case to conduct the elemental contents analysis in soil

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
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: 20181207

Termination date: 20211209