中国激光, 2020, 47 (9): 0911001, 网络出版: 2020-09-16   

基于响应面法的激光诱导击穿光谱实验装置参数优化研究 下载: 756次

Parameter Optimization of Laser-Induced Breakdown Spectroscopy Experimental Device Based on Response Surface Methodology
作者单位
南京航空航天大学机电学院, 江苏 南京 210001
摘要
激光诱导击穿光谱技术(LIBS)具有无需制备样品、多元素同时测量和快速检测等优点,在物质检测领域得到了广泛应用,选择合适的实验装置参数是其取得良好检测效果的重要基础。在单因素实验的基础上,以激光能量、延迟时间和焦深为影响因素,以光谱信背比作为响应指标,建立了一种优化激光诱导击穿光谱实验装置参数的多因素响应面模型。研究了不同影响因素及它们之间的耦合作用对光谱质量的影响,得到的最优实验参数为激光能量114 mJ,延迟时间1.86 μs,焦深1.75 mm。最后,对响应面法所得的最优实验参数进行实验验证,得出谱线信背比均值为7.45,较单因素实验法提高了1.92%,相对标准偏差为3.16%,较单因素实验法下降了0.94%。结果表明,响应面法比单因素实验法更加有效、可靠。
Abstract
Laser-induced breakdown spectroscopy (LIBS) has been widely used in the field of material detection due to its advantages such as no sample preparation, multi-element simultaneous analysis and rapid detection. Selecting appropriate parameters of the experimental device is an important prerequisite for achieving a good detection effect. Based on the single factor experiment, a multi-factor response surface model for optimizing the parameters of the laser-induced breakdown spectroscopy experimental device is established with laser energy, delay time and depth of focus as the influencing factors and with spectral signal-to-background ratio as the response factor. The influence of different influencing factors and their coupling effects on the spectral quality are studied and the obtained optimal experimental parameters are laser energy of 114 mJ, delay time of 1.86 μs, and depth of focus of 1.75 mm. Finally, the optimal experimental parameters obtained by the response surface methodology are verified by experiments. The mean value of the signal-to-background ratio of the spectral line is 7.45, which is 1.92% higher than that by the single factor method, and the relative standard deviation is 3.16%, which is 0.94% lower than that by the single factor method. The results show that the response surface methodology is more effective and reliable than the single factor method.

1 引言

激光诱导击穿光谱技术(Laser-Induced Breakdown Spectroscopy,LIBS)是一种原子发射光谱技术[1],它利用高能脉冲激光轰击样品表面产生等离子体,然后通过对光谱仪采集的等离子体光谱信号进行光谱分析,从而达到确定物质成分的目的。LIBS具有近似无损分析、无需制备样品、多元素同时测量和快速检测等优点[2-4],在艺术品鉴定[5]、冶金[6]和环境工程[7]等领域得到了广泛应用。

激光诱导等离子体的产生和演化是个非常复杂的过程,实验参数如激光能量、延迟时间和焦深等因素会直接影响分析检测的精度与灵敏度。为了改善此问题,许多学者对实验参数优化进行了相关研究。Shaikh等[8]研究了不同激光波长对铝等离子体的激发温度与电子数密度的影响,发现长波长下的等离子体的激发温度较高,而短波长下的等离子体的谱线强度更强,电子密度更大。张雷等[9]以信背比作为参数优化的依据,发现经过优化的LIBS实验装置对煤粉中C元素测量的标准偏差由优化前的6.7%下降至1.6%,测量精度得到了较大提高。李乘等[10]以混合水溶液中选定的分析谱线作为研究对象,以LIBS信号信背比的最大化为优化标准,优化了光谱探测系统的关键实验参数。

传统的参数优化过程都是研究单因素的影响规律,即通过改变一种因素,并保持其他因素不变,研究该因素对检测结果的影响。对于不同的样品与元素,其最佳的实验参数是不同的,若通过单因素实验的方法来寻找每个实验参数的最优值,需要大量的实验测试,步骤十分繁琐,且不同实验参数之间可能存在耦合关系,每个单因素最优值的组合可能并非是当前实验条件下的最佳实验参数组合。考虑多种影响因素及各因素间的交互作用的统计学方法也许是更加有效的办法。孙兰香等[11]利用二次回归正交设计(QROD)法构建了LIBS多因素参数优化模型,结果表明,该模型可以获得比单因素优化方法相对更优的检测结果,且大大减少了实验次数,提高了检测效率。鄢然等[12]通过响应面法建立了尺寸精度与工艺参数之间的数学模型,研究了不同实验参数及其交互作用对选择性激光烧结制件成型精度的影响,并得出了最优的工艺参数。

响应面法相比于正交实验法,可以连续地对实验中的各个因素进行分析,进而找到整个区域中的最佳实验参数组合与响应值的最优值,所得结果更加可靠、合理。因此,基于上述研究,本文以分析谱线的信背比作为响应,建立了多因素响应面优化模型,研究了激光能量、延迟时间、焦深及它们之间的耦合作用对光谱质量的影响,从而得出了最优的实验参数组合,并将结果与单因素最优组合进行对比,验证了基于响应面法所建立的优化模型的有效性。

2 基本原理

响应面法[13-15](Response Surface Methodology,RSM)通过对指定设计空间内的样本点集合进行有限次数的实验设计,回归拟合出响应目标与设计变量之间的函数关系,从而获得设计变量的最优组合,使目标函数达到最优。

当求解问题涉及多目标优化时,可以采用多项式响应面模型来拟合因素与响应之间的函数关系,当实验区域远离响应面的最优区域时,宜采用响应面的一阶模型进行逼近;当实验区域接近或位于最优区域时,可用二阶多项式模型来逼近响应面。其中,一阶与二阶多项式响应面模型[16]可表示为

Y=β0+i=1nβixi,(1)

Y=β0+i=1nβixi+i=1nβiixi+i=1ni<jnβijxixj,(2)

式中:x为影响因素变量;ij为影响因素的编号;n为影响因素的数量;Y为预测响应值;β0βiβiiβij分别表示偏移项、线性偏移、二阶偏移系数和交互作用系数。本文中影响因素为激光能量、延迟时间和焦深,其中焦深表示聚焦点到样品表面的距离。

在实际光谱分析中,光谱强度增大意味着背景噪声也会增大,当背景噪声的增长速度高于谱线强度的增长速度时,分析谱线的信背比就会下降,所以本文选取信背比RSBR作为响应面模型的响应。响应指标RSBR越大,表明该组实验参数的性能越好。

3 实验部分

3.1 实验装置

实验系统原理及组成如图1所示。实验采用调Q开关的Nd∶YAG脉冲激光器(Dawa-200,镭宝光电技术有限公司,北京),波长为1064 nm,最大重复频率为20 Hz,脉冲宽度为8 ns,最大输出能量为200 mJ;光纤光谱仪(AvaSpec-ULS2048CL-4-EVO,爱万提斯科技有限公司,荷兰)的平均光谱分辨率为0.08~0.2 nm,波长探测范围为200~1080 nm。激光经过45°反射镜和平凸透镜L1(焦距为100 mm)聚焦于三维移动平台上的铝合金样品上,激发出等离子体,随后等离子体发射光经平凸透镜L2(焦距为100 mm)与L3(焦距为100 mm)被光纤耦合器采集传输至光纤光谱仪,最后利用计算机对光谱仪得出的光谱数据进行分析处理。其中,脉冲激光器与光纤光谱仪之间的时序关系由数字脉冲发生器(DG535,斯坦福公司,美国)控制。

图 1. LIBS实验系统原理图

Fig. 1. Schematic of LIBS experimental system

下载图片 查看所有图片

3.2 单因素实验设计

实验样品为一块铝合金标样,生产厂家为济南众标。基础工作参数为:脉冲激光能量100 mJ,脉冲频率5 Hz,光谱仪相对激光器Q开关延迟3 μs,积分时间30 μs,焦深1 mm。在上述基础工作参数下,通过单独调整激光能量、延迟时间和焦深,分别研究各因素对光谱分析线信背比的影响,其中单因素及其水平信息如表1所示。每组实验均在标样表面选取20个测量点,每个测量点用激光激发30次,为降低样品表面污渍的影响,取后20次测量数据的平均值作为该点的测量结果,最终在每组实验参数下获得20组实验数据。图2为不同激光能量下铝合金标样在200~500 nm波长范围内的光谱图。

表 1. 单因素实验的因素及水平

Table 1. Factors and levels for single factor experiment

FactorValue
Laser energy /mJ80,90,100,110,120
Delay time /μs1,2,3,4,5
Depth of focus /mm-1,0,1,2,3

查看所有表

图 2. 铝合金光谱图

Fig. 2. LIBS spectra of aluminum alloy

下载图片 查看所有图片

3.3 响应面实验设计

在单因素实验的基础上,利用Design-Export软件中的Box-Behnken实验设计方法进行实验方案的设计,以激光能量、延迟时间及焦深作为实验因素,以信背比RSBR作为响应,每个因素取三个水平,以(0,±1)进行编码,其中“0”代表中心点,“1”代表高水平,“-1”代表低水平。实验因素及水平如表2所示。以Al 265.25作为分析谱线,将各因素及其水平值输入软件,最终确定17组实验参数,详细的实验设计参数及结果如表3所示。

表 2. 响应面实验的因素及水平

Table 2. Factors and levels for response surface methodology

FactorVariablesymbolLevel
-101
Laser energy /mJX1100110120
Delay time /μsX2123
Depth of focus /mmX3123

查看所有表

表 3. 实验设计参数与结果

Table 3. Design parameters and results of experiments

Experiment No.X1 / mJX2 /μsX3 / mmRSBR
1120236.42
2110227.12
3110227.31
4120127.02
5110316.61
6110227.42
7100216.48
8110116.85
9100236.29
10110336.45
11110227.39
12110227.27
13100326.54
14120326.87
15120217.04
16100126.35
17110136.56

查看所有表

4 实验结果与分析

4.1 单因素实验结果分析

根据3.2中的单因素实验数据,以Al 265.25作为分析谱线,激光能量、延迟时间和焦深各单因素对谱线强度及信背比的影响如图3所示。由图3可知,光谱强度随着激光能量的增大而增大,信背比先增大后减小,光谱强度峰值出现在110 mJ处,这主要是因为随着激光能量的增加,样品表面的聚焦区

域更容易获得高能量,产生的高能等离子体的温度和电子数密度也随之增加,但由于自吸收效应的增强,谱线强度在110 mJ之后的上升速度变缓,出现饱和现象;分析谱线的强度随着延迟时间的增大而不断减小,当延迟时间为2 μs左右时,信背比取得最大值,这是因为在等离子体产生初期,背景辐射的衰减速度相对于光谱信号更快,经过一定延迟时间后能够得到光谱质量更高的信号。当焦深从-1 mm到3 mm变化时,光谱强度先增大后减小,信背比在2 mm处达到最大值,这是由于透镜聚焦点位于样品表面上方时,空气在聚焦点处被击穿,激光到达样品表面时已损耗一部分能量,得到的谱线强度比较低;当聚焦点位于样品表面下方不远处时,激光在光束传播方向上的有效作用距离增大,激光接触的物质增多,谱线强度得到增强;而当聚焦点继续向下移动时,光斑变大,作用在样品表面的激光能量密度不断降低,这使得被激发物质的光谱强度降低。因此,在当前实验条件下,激光能量110 mJ、延迟时间2 μs和焦深2 mm为各单因素的最优实验参数。

图 3. 单因素实验结果。(a)激光能量;(b)延迟时间;(c)焦深

Fig. 3. Results of single factor experiment. (a) Laser energy; (b) delay time; (c) depth of focus

下载图片 查看所有图片

4.2 响应面法实验结果分析

利用Design-Expert 软件对表3中的实验结果进行多元回归方程拟合,得到的拟合回归方程为

Y=-41.83+0.79×X1+1.93×X2+2.60×X3-0.0085×X1X2-0.011×X1X3+0.033×X2X3-0.0033×X12-0.27×X22-0.41×X32(3)

进一步对回归方程进行方差分析,结果如表4所示。其中,Y代表显著,N代表不显著;X1X2表示激光能量和延迟时间的交互作用,X1X3X2X3的含义同理;PF代表相关系数的显著性,通常当P小于0. 05时,可认为该因素对目标函数的影响显著,P值越小,所对应的因素的影响越显著;失拟项表示回归方程未能拟合的部分或未考虑到的因素,失拟项的P值在不超过1的情况下,P值越大说明方程的拟合效果越好。Design-Expert软件还给出了残差的正态概率分布图、预测值与实际值分布图,分别如图4、5所示。

表 4. 方差分析

Table 4. Analysis of variance

SourceSum of squareDegree of freedomMean square errorFPSignificant
Model2.3190.2617.430.0005Y
X10.3610.3624.240.0017Y
X20.01210.0120.820.3965
X30.2010.2013.470.0080Y
X1X20.02910.0291.960.2040
X1X30.04610.0463.140.1198
X2X30.004210.00420.290.6088
X120.4710.4731.790.0008Y
X220.3110.3121.380.0024Y
X320.7110.7148.290.0002Y
Residual0.1070.015
Lack of fit0.04730.0161.130.4380N
Pure error0.05640.014

查看所有表

图 4. 残差的正态概率分布图

Fig. 4. Normal probability distribution of residuals

下载图片 查看所有图片

图 5. 预测值与实际值分布图

Fig. 5. Distribution of predicted and actual values

下载图片 查看所有图片

可以看出,X1X3X12X22X32为显著因素,失拟项不显著,残差的正态概率基本分布在一条直线上,实际值与预测值分布图呈现明显的线性关系,表明多元回归模型的拟合精度好,适应性强。在本实验所取的各因素水平范围内,各因素对谱线信背比的影响顺序为: X32, X12,X1, X22,X3,X1X3,X1X2,X2,X2X3

图6~8为各因素的交互作用对谱线信背比影响的三维曲面图和等高线图,其中等高线图的形状越接近椭圆形则表示因素之间的交互作用对响应值的影响越强,越接近圆形则表示影响越弱,等高线上的数字表示RSBR。可以看出,激光能量与焦深的交互作用对响应值的影响最明显,而延迟时间与焦深的交互作用对响应值的影响最弱;不同因素之间的交互作用的三维曲面均存在峰值点,这意味着响应值在当前各因素水平范围内存在最大值。经分析可知,当响应值最大时,激光能量为114 mJ,延迟时间为1.86 μs,焦深为1.75 mm,多元回归模型预测的谱线信背比为7.36。

图 6. 激光能量与延迟时间的交互作用对谱线信背比的影响。(a)响应面图;(b)等高线图

Fig. 6. Effect of interaction between laser energy and delay time on SBR of spectral line. (a) Response surface diagram; (b) contour map

下载图片 查看所有图片

图 7. 激光能量与焦深的交互作用对谱线信背比的影响。(a)响应面图;(b)等高线图

Fig. 7. Effect of interaction between laser energy and depth of focus on SBR of spectral line. (a) Response surface diagram; (b) contour map

下载图片 查看所有图片

图 8. 延迟时间与焦深的交互作用对谱线信背比的影响。(a)响应面图;(b)等高线图

Fig. 8. Effect of interaction between delay time and depth of focus on SBR of spectral line. (a) Response surface diagram; (b) contour map

下载图片 查看所有图片

为了验证响应面法用于LIBS实验装置参数优化的有效性,分别采用响应面法与单因素实验法所得的最佳实验参数进行了5次实验验证,单因素实验法所得的谱线信背比均值为7.31,相对标准偏差为4.10%;响应面法所得的谱线信背比均值为7.45,较单因素实验法提高了1.92%,相对标准偏差为3.16%,较单因素实验法下降了0.94%,与模型预测值的相对误差为1.22%。由此可见,利用响应面法优化激光诱导击穿光谱实验装置参数是可行的,且相对于传统单因素实验方法更加有效、可靠。

5 结论

基于响应面法提出了一种优化激光诱导击穿光谱实验装置参数的方法,以谱线信背比作为响应指标,以激光能量、延迟时间和焦深为影响因素,研究了不同影响因素及它们之间的耦合作用对光谱质量的影响,最终得到的最优实验参数为:激光能量114 mJ,延迟时间1.86 μs,焦深1.75 mm。对该结果进行实验验证,得出谱线信背比均值为7.45,较单因素实验法提高了1.92%,相对标准偏差为3.16%,较单因素实验法下降了0.94%。结果表明,在当前实验条件下,利用响应面法能够有效优化激光诱导击穿光谱实验装置参数,且优化效果优于传统单因素实验法。

参考文献

[1] 王金梅, 郑慧娟, 郑培超, 等. 正交再加热双脉冲激光诱导黄连等离子体的光谱特性[J]. 中国激光, 2018, 45(7): 0702006.

    Wang J M, Zheng H J, Zheng P C, et al. Spectral characteristics of coptis chinensis plasma induced by orthogonal re-heating double-pulse laser[J]. Chinese Journal of Lasers, 2018, 45(7): 0702006.

[2] Zhu Z H, Li J M, Guo Y M, et al. Accuracy improvement of boron by molecular emission with a genetic algorithm and partial least squares regression model in laser-induced breakdown spectroscopy[J]. Journal of Analytical Atomic Spectrometry, 2018, 33(2): 205-209.

[3] Rohwetter P, Yu J, Mejean G, et al. Remote LIBS with ultrashort pulses: characteristics in picosecond and femtosecond regimes[J]. Journal of Analytical Atomic Spectrometry, 2004, 19(4): 437-444.

[4] 应璐娜, 周卫东. 对比分析多种化学计量学方法在激光诱导击穿光谱土壤元素定量分析中的应用[J]. 光学学报, 2018, 38(12): 1214002.

    Ying L N, Zhou W D. Comparative analysis of multiple chemometrics methods in application of laser-induced breakdown spectroscopy for quantitative analysis of soil elements[J]. Acta Optica Sinica, 2018, 38(12): 1214002.

[5] De Giacomo A. Dell'Aglio M, De Pascale O, et al. Laser induced breakdown spectroscopy methodology for the analysis of copper-based-alloys used in ancient artworks[J]. Spectrochimica Acta Part B: Atomic Spectroscopy, 2008, 63(5): 585-590.

[6] 王琦, 陈兴龙, 王静鸽, 等. 影响激光诱导等离子体稳定性的因素研究[J]. 光学学报, 2014, 34(6): 0630002.

    Wang Q, Chen X L, Wang J G, et al. Research on factors affecting the stability of laser-induced plasmas[J]. Acta Optica Sinica, 2014, 34(6): 0630002.

[7] Kumar R, Rai A K, Alamelu D, et al. Monitoring of toxic elements present in sludge of industrial waste using CF-LIBS[J]. Environmental Monitoring and Assessment, 2013, 185(1): 171-180.

[8] Shaikh N M, Hafeez S, Rashid B, et al. Spectroscopic studies of laser induced aluminum plasma using fundamental, second and third harmonics of a Nd∶YAG laser[J]. European Physical Journal D, 2007, 44(2): 371-379.

[9] 张雷, 马维光, 闫晓娟, 等. 激光诱导击穿光谱实验装置的参数优化研究[J]. 光谱学与光谱分析, 2011, 31(9): 2355-2360.

    Zhang L, Ma W G, Yan X J, et al. Research on parameters optimization of laser-induced breakdown spectroscopy based experimental device[J]. Spectroscopy and Spectral Analysis, 2011, 31(9): 2355-2360.

[10] 李乘, 姚关心, 杨新艳, 等. 液相基质激光诱导击穿光谱的关键实验参数优化[J]. 激光与光电子学进展, 2019, 56(7): 073002.

    Li C, Yao G X, Yang X Y, et al. Key experimental parameter optimization for laser induced breakdown spectroscopy of liquid matrix[J]. Laser & Optoelectronics Progress, 2019, 56(7): 073002.

[11] 孙兰香, 辛勇, 丛智博, 等. 通过二次回归正交设计对激光诱导击穿光谱实验参数优化建模[J]. 光学学报, 2014, 34(5): 0530003.

    Sun L X, Xin Y, Cong Z B, et al. Optimization modeling for experimental parameters of laser-induced breakdown spectroscopy by quadratic regression orthogonal design[J]. Acta Optica Sinica, 2014, 34(5): 0530003.

[12] 鄢然, 李浩, 李军超, 等. 基于响应面法的聚苯乙烯粉末选择性激光烧结成型工艺参数优化[J]. 中国激光, 2019, 46(3): 0302015.

    Yan R, Li H, Li J C, et al. Process parameters optimization of polystyrene powder selective laser sintering based on response surface methodology[J]. Chinese Journal of Lasers, 2019, 46(3): 0302015.

[13] 李莉, 张赛, 何强, 等. 响应面法在试验设计与优化中的应用[J]. 实验室研究与探索, 2015, 34(8): 41-45.

    Li L, Zhang S, He Q, et al. Application of response surface methodology in experiment design and optimization[J]. Laboratory Research and Exploration, 2015, 34(8): 41-45.

[14] 伍彦伟, 刘会霞, 李品, 等. 激光透射焊接聚碳酸酯和聚苯醚实验研究[J]. 中国激光, 2015, 42(5): 0503010.

    Wu Y W, Liu H X, Li P, et al. Experimental study of laser transmission welding between polycarbonate and polyphenylene oxide[J]. Chinese Journal of Lasers, 2015, 42(5): 0503010.

[15] Algin H M. Optimised design of jet-grouted raft using response surface method[J]. Computers and Geotechnics, 2016, 74: 56-73.

[16] 刘其蒙, 黄俊, 王克鸿, 等. 基于响应面法的离体皮肤组织激光融合切口性能多元非线性回归模型[J]. 中国激光, 2018, 45(8): 0807002.

    Liu Q M, Huang J, Wang K H, et al. Multivariate nonlinear regression model of laser fusion in vitro skin tissue incision performance based on response surface methodology[J]. Chinese Journal of Lasers, 2018, 45(8): 0807002.

潘立剑, 陈蔚芳, 周晏锋, 崔榕芳. 基于响应面法的激光诱导击穿光谱实验装置参数优化研究[J]. 中国激光, 2020, 47(9): 0911001. Pan Lijian, Chen Weifang, Zhou Yanfeng, Cui Rongfang. Parameter Optimization of Laser-Induced Breakdown Spectroscopy Experimental Device Based on Response Surface Methodology[J]. Chinese Journal of Lasers, 2020, 47(9): 0911001.

本文已被 3 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!