中国·永利集团(3044am-YL认证)·Holdings Limited

服务热线:
86-0755-23229824
您当前所在位置: 首页>>应用案例>>能源环保
录井专用型激光诱导击穿光谱仪测定岩屑中的8种元素

image.png

1、引言

岩屑录井工作是油气勘探和开发的基本过程之一,取不同井深的岩屑样品,分析其元素组成及含量,可以为油气开采和评估提供第一手的参考信息,对油气勘探具有重要意义激光诱导击穿光谱(LIBS)技术具有实时在线多元素同时检测样品无需复杂前处理等优点与目前岩屑录井中最常用的X射线荧光分析技术(XRF)相比,LIBS技术的检测过程更为快速便捷,没有XRF仪的辐射风险,并可分析XRF不能很好识别的轻元素作为近年来迅速发展的元素分析技术,LIBS技术的特点与油气勘探领域的岩屑录井需求十分契合,具有很好的应用前景和推广潜力但由于LIBS岩屑样品分析中基体效应谱线干扰和自吸收效应等问题,采用传统的一元定标法和内定标法对样品的元素进行定量分析通常不能满足实际应用的需要获得更为精准的元素定量分析结果一直是LIBS技术在油气勘探领域应用推广所面临的挑战之一近年来,在LIBS数据分析研究中,基于多变量统计学习的偏最小二乘回归(PLSR)遗传算法(GA)人工神经网络(ANN)支持向量机回归(SVR)得到更多运用Sarkar等将PLSR运用于钢铁样品中CoMoCrNi元素定量分析邹孝恒等将遗传算法与偏最小二乘法结合(GA-PLS)用于土壤元素定量分析,结果表明将遗传算法作为一种谱线选择预处理方法可以改善PLS模型的预测能力Haddad等通过手持式LIBS仪器采集谱图,运用ANN算法对土壤进行分类和元素含量定量分析,也取得了较好的效果张莹等在标样较少的情况下,将SVR算法运用于钢铁样品中的CrNi元素定量分析,对比了以峰强度和峰面积作为模型输入变量的模型效果,发现以峰面积作为变量获得的模型效果更好胡丽等将SVR算法结合LIBS技术运用于水中Pb元素预测,测试集的相对标准偏差和相对误差平均值分别为45%121%

目前,LIBS技术已广泛应用于铝合金钢铁土壤水体食品和煤炭等元素的定量分析,而针对钻井岩屑的元素定量研究还较少本研究选择在小样本非线性学习方面具有优势独特的SVR算法,建立优化8种元素(SiAlCaKMgMnTiFe)定量模型针对实际应用需求开发专用型仪器,并用于实践工作,是LIBS技术推广的关键

2.实验部分

2.1仪器

测量仪器采用Nd:YAG激光器,激光波长1064nm,最大发射频率1Hz,最大脉冲能量为100mJ光谱检测范围为200800nm,分辨率为01nm仪器内集成的延时器模块可以为光谱仪提供05200μs的外触发延时信号激光束通过平凸透镜聚焦于岩屑样品表面,产生的等离子体发出包含元素特性的原子发射谱线,经过合适的延时时间后通过光纤探头收集由光谱仪光电转换完成数据采集样品置于仪器三维样品台上,依据软件提供的参考图像,通过旋转三维平移台改变样品的检测位置经过优化,本研究采用的延时时间为4μs,激光能量为100mJ,激光频率为1Hz

image.png 

1仪器原理示意图

2.2样品及样品制备

18个岩屑样品来自西南某油气录井现场,基本覆盖应用现象的岩性种类,编号为YW1YW18将样品送至国家地质实验测试中心,样品玻璃熔片制样法处理后,采用XRF技术进行检测,得到各个样品中SiAlCaKMgMnTiFe元素浓度及岩性信息(1)由于PDC钻头广泛应用,钻井取出的岩屑样品十分细碎,甚至呈粉末状,处于较为一致状态为了使样品处于更为一致的状态,减少或消除环境及基体差异的影响,取15g各干燥的岩屑样品,置于玛瑙钵中研磨,利用压片机在20MPa压强下保持20s,压制成直径约30mm薄厚均匀的块状压片样品,用于谱图采集考虑到更大的样品量所获得的模型更为可靠,在YW1YW18中随机选择两种样品,按一定比例混合,经过称量搅拌研磨和压片后,另外配制了42个样品,编号为A1C42,按比例计算的成分含量见表11中样品编号标有*号的样品为测试集,其余样品作为训练集由于每个样品重复采集6次,训练集和测试集分别包含240120LIBS谱图

1岩屑样品元素质量分数(wt%)

image.pngimage.png 

image.png 

3.结果与讨论

3.1 SVR建模简介

支持向量机(SVM)最早作为一种二分类算法由贝尔实验室的Vapnik等在1995年提出,而后逐步运用于回归分析(SVR)选取合适的训练集和测试集样本数据,确立输入变量和因变量后,寻找最佳参数的SVR回归模型,可以归纳为以下步骤:

(1) 归一化对训练集和测试集数据进行归一化处理

(2) 选择核函数依据数据特点选择核函数,在本研究中选用径向基核函数(RBF),其表达式为

image.png

式中,xj为核函数中心点在高维空间的位置,xi为输入输入变量xixj绝对值为高维空间两点间的欧式距离γ为标准高斯方程中的参数

(3) 参数优化在训练集内,通过网格搜寻(Gridsearch)的方式,寻找模型最佳的关键参数惩罚因子CRBF核的关键参数γ评估指标为十折交叉验证方式得到的平方相关系数R2和交叉验证均方根误差(RMSECV)

(4) 验证将确立的最佳参数建立的模型用于测试集数据预测,并评估其预测效果评估指标为预测均方根误差(RMSEP)和平均百分比误差(MPE)考虑到算法与仪器的集成,算法采用C#语言编写,用开源的LIBSVM算法库(Library)二次开发,更为详尽算法推导过程以及参数寻优过程介绍可以参见文献

3.2光谱数据评估和预处理

由典型的岩屑样品LIBS谱图(2)可见,获得的岩屑样品光谱具有丰富的元素特征谱线,且具有较好的信噪比

image.png 

2典型岩屑样品谱图

以相对标准差(RSD)为指标,对YW3YW4YW8YW16样品分别采集12张谱图,考察对SiI28816FeII27557KI76649谱线的峰强度由表2可知,获得的数据重复性和稳定性较好此外,为建立更稳定的模型,对原始谱图各个光谱强度变量减去相应基线值进行基线校正,随后对谱图进行五点平滑滤波降低噪声,最后为使参数寻优计算更容易正确的收敛,将LIBS谱图的所有谱线强度值除以强度最大值,对所有谱图数据进行归一化处理

 

 

 

 

2部分样品部分特征谱线峰强度相对标准差

image.png 

3.3 SVR模型输入变量的选择

岩屑样品元素组成丰富,特征谱线丰富,同时也包含了交叠干扰无用冗余的信息输入变量的合理选择,是模型预测效果的关键因素分别以全部光谱数据特征谱线聚集数据段部分优选特征谱线峰面积和峰强度分别作为输入变量,建立模型,考察获得的模型效果

在获取样品LIBS谱图后,通过谱图与NIST原子光谱数据库对比,选择信号强且干扰较少的一些特征谱线,分别以峰强度与峰面积作为变量,建立模型,所选择的特征谱线见表3考虑到选择部分信号强干扰少的谱线,不可避免地忽略一些其它相关信息,选择所需分析元素特征谱线富集的8个光谱数据段(238239nm250253nm256276nm278290nm308310nm390399nm402404nm765770nm)作为输入变量,建立模型此外,将200800nm范围内所有波长的信号强度数据作为输入变量,建立了相应对比模型SiMnFe为例,最终各SVR模型输入变量与训练结果见表4,以光谱数据段及全谱所有数据作为变量,虽然包含了更多特征信息,但由于无法很好地排除冗余无用信息干扰,整体效果较差,尤其对于含量较少特征谱线强度较弱的元素效果尤为明显以部分优选特征谱线峰面积作为输入变量的效果最好,包含了峰高峰宽形状等更全面信息,优于峰强度作为输入变量的模型效果

 

 

3支持向量机回归模型特征光谱变量选择

image.png 

4 SiMnFe不同输入变量的支持向量机回归模型训练结果

image.png 

3.4定量分析结果

根据3.3节的分析,以表3中的特征光谱峰面积为输入变量依照3.1节的步骤,以RBF为核函数,交叉验证网格搜索方式得到最优的参数Cγ后,建立相应模型对测试集中8种元素的含量进行预测以横坐标为预测值,纵坐标为样品XRF检测值,各元素模型的最佳参数及预测效果见图38种元素预测的平均百分比误差(MPE)分别为Si568%Al722%Ca745%K976%Mg879%Mn119%Ti114%Fe104%,满足录井工作的要求

3.5应用测试

依据录井工作实际需求,将建立并优化好的各元素SVR定量模型集成于软件中,并对仪器硬件参数软件功能及操作界面进行修改,将仪器改造成为录井专用型号在录井现场按照岩屑录井现场工作流程,对岩屑样品研磨压片采集谱图预测定量结果,最后仪器软件以报表形式进行存档整个过程中,研磨压片过程耗时12min,可考虑多样品同时压片,提高效率,而从采集谱图到计算定量结果耗时仅为数秒测试结果表明,集成了算法的整套仪器系统可以满足岩屑录井工作需求

4结论

将自行研制的仪器结合SVR算法应用于岩屑样品中SiAlCaKMgMnTiFeLIBS定量分析,对LIBS谱图数据进行基线扣除平滑滤波及归一化处理后,考察了不同输入变量情况下SVR模型的训练效果结果表明,以峰面积为输入变量,得到的训练效果最好,各元素的平均百分比误差在5.68%11.9%之间本研究针对岩屑录井工作的实际应用需求,优化仪器硬件参数,并将研究确立的定量模型集成于仪器中,使算法结合仪器成为专用型仪器,可用于应用现场此外,本研究中岩屑标样岩性种类元素种类仍不全面,且对于含量较低的元素预测效果还有提升空间,提高LIBS信号的稳定性信噪比,建立更大岩屑样本数据库及优化建立基于此数据库的更多元素种类的定量模型将是下一步研究的重点

image.png 

3 SiAlCaKMgMnTiFe8种元素预测效果图


推荐

LSpec-LIBS400分体式LIBS激光诱导击穿光谱系统

LSpec-LIBS400分体式LIBS激光诱导击穿光谱系统主要由激光器、高分辨率光谱仪、LIBS光路收集探头、XYZ样品窗、触发延迟控制器等组成,由于其采用分体式搭建方式,可扩展非常灵活,非常适合科研研究、LIBS光学应用实验、光学应用中心等用户,可以非常方便的灵活选用配置激光器和光谱仪。

image.png 

LSpec-LIBS800 LIBS激光诱导击穿光谱系统

LSpec-LIBS800分LIBS激光诱导击穿光谱系统主要由激光器、高分辨率光谱仪、LIBS光路收集探头、XYZ样品窗、触发延迟控制器等组成,由于其采用一体化,可扩展非常灵活,非常适合科研研究、LIBS光学应用实验、光学应用中心等用户,可以非常方便的灵活选用配置激光器和光谱仪。

image.png 

手持式LIBS激光诱导击穿光谱系统

手持式LIBS激光诱导击穿光谱系统附带免费的客户端软件,它提供了强大的工具来下载测试数据、编辑牌号库、生成测试报告和编译校准文件。对于高级用户,它提供了非常有用的附加功能来满足用户的专业要求。

image.png


Copyright © 2020 All Rights Reserved 3044am永利集团(深圳) 有限公司·版权所有 备案号:粤ICP备18141551号
XML 地图