CN107340284A - 一种元素定量分析方法及装置 - Google Patents

一种元素定量分析方法及装置 Download PDF

Info

Publication number
CN107340284A
CN107340284A CN201710550399.4A CN201710550399A CN107340284A CN 107340284 A CN107340284 A CN 107340284A CN 201710550399 A CN201710550399 A CN 201710550399A CN 107340284 A CN107340284 A CN 107340284A
Authority
CN
China
Prior art keywords
spectral line
measured
intensity
sample
target
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
CN201710550399.4A
Other languages
English (en)
Other versions
CN107340284B (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.)
Wuhan Hongshi Technology Co ltd
Original Assignee
Hubei Engineering 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 Hubei Engineering University filed Critical Hubei Engineering University
Priority to CN201710550399.4A priority Critical patent/CN107340284B/zh
Publication of CN107340284A publication Critical patent/CN107340284A/zh
Application granted granted Critical
Publication of CN107340284B publication Critical patent/CN107340284B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

Landscapes

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

Abstract

本发明实施例提供一种元素定量分析方法及装置,属于光谱分析技术领域。该元素定量分析方法首先获取定标样品集中的每个样品的待测元素的浓度和待测元素的至少两条目标谱线的谱线强度,然后根据预设的每条目标谱线的权重系数和每个样品中每条目标谱线的谱线强度,计算每个样品中待测元素的加权谱线强度,最后再根据每个样品中待测元素的浓度和计算得到的每个样品中待测元素的加权谱线强度进行线性拟合,确定用于对该待测元素定量分析的定标曲线。这种元素定量分析方法采用多条谱线进行加权处理,避免了单条谱线在定量分析中,因等离子体的波动或外界干扰引起的谱线波动,而造成的分析不准确,有效的提高了系统的分析准确度和稳定性。

Description

一种元素定量分析方法及装置
技术领域
本发明涉及光谱分析技术领域,具体而言,涉及一种元素定量分析方法及装置。
背景技术
随着经济与科技的飞速发展,在许多领域诸如采矿、冶金、环保、化工、能源、文物珠宝、食品安全、生物制药等中,需要对物质成分进行定性或定量分析。激光诱导击穿光谱(Laser-induced Breakdown Spectroscopy,LIBS),是一种基于原子发射光谱的物质成分分析技术,其原理是将激光脉冲经过透镜聚焦之后照射到被测样品靶材的表面,激发产生等离子体,同时对等离子体演化过程中辐射产生的元素特征光谱进行采集,并分析出该样品中的元素种类及其含量。但是在现有LIBS中,常采用的定量分析方法容易因谱线波动或受其它谱线干扰造成分析结果不准确。
发明内容
有鉴于此,本发明的目的在于提供一种元素定量分析方法及装置,以改善上述问题。
本发明较佳实施例提供一种元素定量分析方法,所述方法包括:获取定标样品集中的每个样品的待测元素的浓度和所述待测元素的至少两条目标谱线的谱线强度;根据每条目标谱线预设的初始权重系数和每个样品中每条目标谱线的谱线强度,计算每个样品中所述待测元素的加权谱线强度;根据获取到的每个样品中所述待测元素的浓度和计算得到的每个样品中所述待测元素的加权谱线强度进行线性拟合,确定用于对该待测元素进行定量分析的第一定标曲线。
本发明另一较佳实施例提供一种元素定量分析装置,包括:谱线强度获取模块,用于获取每个样品的待测元素的浓度和所述待测元素的至少两条目标谱线的谱线强度;加权谱线强度计算模块,用于根据每条目标谱线预设的初始权重系数和每个样品中每条目标谱线的谱线强度,计算每个样品中所述待测元素的加权谱线强度;第一定标曲线确定模块,用于根据获取到的每个样品中所述待测元素的浓度和计算得到的每个样品中所述待测元素的加权谱线强度进行线性拟合,确定用于对该待测元素进行定量分析的第一定标曲线。
本发明实施例提供的元素定量分析方法及装置,首先获取每个样品中待测元素的浓度和待测元素的至少两条目标谱线的谱线强度,然后根据预设的每条目标谱线的初始权重系数和每个样品中每条目标谱线的谱线强度,计算每个样品中待测元素的加权谱线强度,最后再根据每个样品中待测元素的浓度和计算得到的每个样品中待测元素的加权谱线强度进行线性拟合,确定用于对该待测元素定量分析的定标曲线。这种元素定量分析方法采用多条谱线进行加权处理,避免了单条谱线在定量分析中,因等离子体的波动或外界干扰引起的谱线波动,而造成的分析不准确,有效的提高了系统的分析准确度和稳定性。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
图1为本发明实施例提供的一种元素定量分析方法的应用场景示意图;
图2为本发明实施例提供的一种执行元素定量分析方法的计算设备的方框示意图;
图3为本发明实施例提供的一种元素定量分析方法的流程图;
图4为本发明实施例提供的一示例中生铁样品的光谱图;
图5为本发明实施例提供的另一种元素定量分析方法的流程图;
图6为本发明实施例提供的一种元素定量分析装置的功能模块框图。
图标:100-计算设备;110-元素定量分析装置;120-存储器;130-处理器;1102-谱线强度获取模块;1104-加权谱线强度计算模块;1106-第一定标曲线确定模块;1108-参数遍历模块;1110-第二定标曲线确定模块;1112-决定系数计算模块;1114-目标定标曲线确定模块;1116-平均谱线强度计算模块;1118-浓度求取模块。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。
因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请参阅图1,是本发明实施例提供的一种元素定量分析方法的应用场景示意图。该场景中采用便携式LIBS系统获取样品的光谱,并将光谱输出至计算机进行后续光谱分析,例如本公开涉及的待测元素的定量分析。
请参阅图2,是本发明实施例提供的一种执行元素定量分析方法的计算设备100的方框示意图。计算设备可以是图1中所示的计算机,但不限于此。该计算设备包括元素定量分析装置110、存储器120和处理器130。
存储器120和处理器130之间直接或间接的电性连接,以实现数据的传输或交互。元素定量分析装置110包括至少一个可以软件或固件的形式存储于存储器120中或固化在计算设备100的操作系统中的软件功能模块。处理器130用于执行存储器120中存储的可执行模块,例如元素定量分析装置110所包括的软件功能模块及计算机程序,以根据从LIBS系统获取的光谱进行待测元素定量分析。
请参阅图3,是本发明实施例提供的一种元素定量分析方法的流程图。本实施例中的方法适用于计算设备100,以对样品中的待测元素进行定量分析。所应说明的是,本方法不以图3及下述的具体顺序为限制。下面将对图3所示的各步骤进行详细阐述。
步骤S101,获取定标样品集中的每个样品的待测元素的浓度和该待测元素的至少两条目标谱线的谱线强度。
本实施例中,定标样品集中的样品的数量为多个,各样品中待测元素的浓度可以获知,而且浓度呈一定的梯度变化。
本实施例中,作为一种实施方式,可以采用图1中所示的LIBS系统获取各样品的光谱,并针对获取到的光谱选择待测元素的两条或两条以上的目标谱线。换言之,不同样品中所选择的待测元素的目标谱线是相同的。
作为一个示例,如图4所示,是本示例中某一生铁样品的光谱图。假设待测元素为微量元素Mn(锰)。参照图4,可以选择待测元素Mn的波长分别为403.08nm、403.31nm和403.45nm的三条Mn元素谱线作为目标谱线。这样,在其他生铁样品中,也同样必须选择波长分别为403.08nm、403.31nm和403.45nm的三条Mn元素谱线作为其目标谱线。
步骤S103,根据每条目标谱线的预设的初始权重系数和每个样品中每条目标谱线的谱线强度,计算每个样品中待测元素的加权谱线强度。
目标谱线的谱线强度可以从光谱中获知,例如在上述示例的图4中,可以获知每条Mn元素的目标谱线的强度。
本实施例中,每条目标谱线对应有预设的初始权重系数,并且所有的初始权重系数之和为1。将每条目标谱线的初始权重系数和其谱线强度相乘之后,再求和,可以计算出每个样品中待测元素的加权谱线强度。
作为一个示例,分别对m个同类基体的标准样品(比如一套微合金钢标准样品,样品1#,样品2#,…,样品m#),选取其中待测元素(例如Mn元素)的n条目标谱线。以样品1#为例,其n条目标谱线的强度为[I1,1,I1,2,I1,3,...,I1,n]。为每条目标谱线预设初始的权重系数,构成的初始权重系数组为[g1,g2,g3,...,gn]。将样品1#的n条目标谱线的强度[I1,1,I1,2,I1,3,...,I1,n]与初始权重系数组为[g1,g2,g3,...,gn]对应相乘(点乘),得到该样品1#中待测元素的加权谱线强度y1,即:
[I1,1,I1,2,I1,3,...,I1,n]·[g1,g2,g3,...,gn]T=y1
那么,对于m个样品,可以通过下式中的矩阵形式表达加权谱线强度的计算:
其中,yi,i=1,2,...,m表示样品i#中待测元素的加权谱线强度。
需要说明的是,待测元素的多条目标谱线的初始权重系数可以相同,也可以不同,但总和须为1。例如,n条目标谱线的初始权重系数可以均为1/n。
步骤S105,根据获取到的每个样品中待测元素的浓度和计算得到的每个样品中待测元素的加权谱线强度进行线性拟合,确定用于对该待测元素定量分析的第一定标曲线。
仍以步骤S103中的示例为例,令m个样品的待测元素的浓度记作C=[c1,c2,c3..,cm]T。其中,ci表示样品i#中待测元素的浓度,i=1,2,3,…,m。
根据获取到的定标样品集中每个样品中待测元素的浓度和计算得到的每个样品中待测元素的加权谱线强度进行线性拟合,即,根据下面示出的各个拟合点(c1,y1),(c2,y2),(c3,y3),...,(cm,ym)进行线性拟合,得到对应该待测元素的第一定标曲线,也可以称之为初始定标曲线:
其中,a和b为拟合参数。
本实施例中提供的初始定标曲线虽然可能非最优权重系数组对应的定标曲线,但其仍然可以一定程度的避免因等离子体的波动或外界干扰引起的谱线波动所造成的分析结果失准,一定程度的提高系统的分析准确度和稳定性。
请参阅图5,是本发明实施例提供的另一种元素定量分析方法的流程图。如图5所示,本实施例中元素定量分析方法包括:
步骤S201,获取定标样品集中的每个样品的待测元素的浓度和该待测元素的至少两条目标谱线的谱线强度。
步骤S203,根据每条目标谱线的预设的初始权重系数和每个样品中每条目标谱线的谱线强度,计算每个样品中待测元素的加权谱线强度。
步骤S205,根据获取到的每个样品中待测元素的浓度和计算得到的每个样品中待测元素的加权谱线强度进行线性拟合,确定用于对该待测元素定量分析的第一定标曲线。
本实施例中步骤S201~步骤S205的具体实现方法可参照上述图3所示步骤S101~步骤S105的详细描述,在此不再赘述。
步骤S207,采用预设算法以预设值为步进单位基于所述初始权重系数组进行参数遍历。
步骤S209,根据每次参数遍历得到的权重系数组,重新计算每个样品中所述待测元素的加权谱线强度。
本实施例中,所述预设算法例如可以是,穷举法或蚂蚁算法,但不限于此。
本实施例中,所述预设值可以是0.01,但不限于此。
不失一般性地,在上述的m个样品的示例中,以穷举法为例,示例性地阐述以0.01为步进单位基于初始权重系数组[g1,g2,g3,...,gn]进行参数遍历的实现方式。
详细地,作为一种实施方式,可以将初始权重系数组分成两组,分别为变动组A和不动组B。例如,将初始权重系数组中的任意2个权重系数抽取出来构成变动组,其余的n-2个权重系数构成不动组,则共有Cn 2种分组方式。例如,变动组A可以是[g1,g2],相应地,不动组B可以是[g3,g4,...,gn]。
进行穷举遍历时,外层大循环:令A组、B组一个每次增大0.01,另一个每次减小0.01,直至减小的那一组减至0为止。
第一内层循环:先保持B组内的各权重系数不变,使A组内的某一部分权重系数每次增大0.01,相同个数的另一部分权重系数同时每次减小0.01,直至减小的那一部分减至0为止。例如,A为[g1,g2]时,可以令g1,g2中的一个(例如g1)每次增大0.01,另一个(例如g2)每次减小0.01,直至减小的那一个为0为止。
第二内层循环:将B组再按照上述分组方式进一步分成两组,其中一组为变动组A1,另一组为不动组B1。例如,B组为[g3,g4,...,gn]时,选取其中的任意2个权重系数构成变动组A1,其余的权重系数构成不动组B1,共有种组合方式。例如,n=5,则[g3,g4]可以为A1组,[g5]可以为B1组,令[g3,g4]的权重系数之和每次增大0.01,[g5]每次减小0.01,直至g5为0为止。
第三层内循环:令[g3,g4]中的一个权重系数(例如g3)每次增大0.01,另一个(例如g4)每次减小0.01,直至减小至0为止。
采用穷举法时,每循环一次(参数遍历一次)得到一组权重系数,根据每次得到的权重系数组重新计算一次每个样品中待测元素的加权谱线强度,计算方法参照上述步骤S203的描述。
步骤S211,根据每个样品中所述待测元素的浓度和上述重新计算得到的每个样品中待测元素的加权谱线强度进行线性拟合,确定用于对该待测元素定量分析的第二定标曲线。
不难理解,第二定标曲线有多个。
该步骤S211的具体实现过程可以参照上述步骤S205和步骤S105的阐述。
步骤S213,计算所述第一定标曲线和每个第二定标曲线的决定系数。
所述决定系数也称作拟合优度,表征回归直线对观测值的拟合程度。决定系数R2的值越大表示回归直线对观测值的拟合程度越好,R2的最大值为1。
分别计算出第一定标曲线和多个第二定标曲线的决定系数,并进行比较,以便选择最优的定标曲线。
当然,可以理解的是,具体实现时可以预设一变量Max,然后计算出第一定标曲线的决定系数R0 2,并将R0 2赋值给Max。在后续参数遍历时,每循环一次则计算出当前第二定标曲线的决定系数R1 2,并将R1 2与当前的Max值比较,若R1 2大于Max,则将R1 2赋值给Max并继续下一次循环,若R1 2小于或等于Max,则直接继续下一次循环,如此直至循环结束。结束后,Max值即为最大的决定系数值。
步骤S215,选取决定系数最大的定标曲线作为待测元素的目标定标曲线,并记录下此时待测元素多条目标谱线的权重系数。
本实施例中提供的目标定标曲线是最优权重系数组对应的定标曲线,其可以有效的避免因等离子体的波动或外界干扰引起的谱线波动所造成的分析结果失准,明显提高系统的分析准确度和稳定性。
进一步地,本实施例中,在获取目标定标曲线后,可以根据目标定标曲线以及对应的权重系数组计算出未知浓度的待测样品的浓度值。
详细地,计算未知浓度的待测样品浓度值的实现方式可以是:首先,获取待测样品中待测元素的多条目标谱线的谱线强度,所述待测样品中待测元素的浓度未知;然后,根据所述目标定标曲线对应的最优权重系数组和上述获取到的所述待测样品中待测元素的每条所述目标谱线的谱线强度,计算该待测样品中待测元素的加权谱线强度;最后,将计算得到的待测样品中待测元素的加权谱线强度代入目标定标曲线,求出待测样品中待测元素的浓度。
以一Mn元素浓度未知的生铁样品为示例,假设求取生铁样品中Mn元素的目标定标曲线时,选择的是波长分别为403.08nm、403.31nm和403.45nm的三条Mn元素谱线作为目标谱线。那么,在求该Mn元素浓度未知的生铁样品中Mn元素浓度时,可以通过LIBS系统获得该样品的光谱,然后从光谱中获知403.08nm、403.31nm和403.45nm的三条Mn元素目标谱线的强度。而后,根据目标定标曲线对应的最优权重系数组和目标谱线的谱线强度计算出该样品中Mn元素的加权谱线强度。最后,将加权谱线强度代入目标定标曲线的方程表达式中,求出Mn元素的浓度值。
较佳地,可以依照上述方法多次采集待测样品的光谱,然后获得多个浓度值,然后再求多个浓度值的平均值作为最终结果。
当然,作为另一种实施方式,还可以多次采集待测样品中待测元素的光谱,然后根据多次获取结果,计算每条目标谱线的平均谱线强度,再根据目标定标曲线对应的最优权重系数组和上述获取到的每条目标谱线的平均谱线强度,计算该待测样品中待测元素的加权谱线强度,最后将计算得到的待测样品中待测元素的加权谱线强度代入目标定标曲线,求出待测样品中待测元素的浓度。
本实施例提供的元素定量分析方法不仅适用于便携式物质成分分析仪,而且适用于其他激光诱导击穿光谱的元素成分定量分析仪器。
请参阅图6,是本发明实施例提供的一种元素定量分析装置110的功能模块框图。该元素定量分析装置110包括谱线强度获取模块1102、加权谱线强度计算模块1104、第一定标曲线确定模块1106、参数遍历模块1108、第二定标曲线确定模块1110、决定系数计算模块1112、目标定标曲线确定模块1114、平均谱线强度计算模块1116、浓度求取模块1118。
所述谱线强度获取模块1102,用于获取每个样品中待测元素的浓度和所述待测元素的至少两条目标谱线的谱线强度。该谱线强度获取模块1102可用于执行图3中所示的步骤S101和图5中所示的步骤S201,具体的操作方法可参照上述步骤S101和步骤S201的描述。
所述加权谱线强度计算模块1104,用于根据每条目标谱线的预设的初始权重系数和每个样品中每条目标谱线的谱线强度,计算每个样品中所述待测元素的加权谱线强度,还用于根据每次参数遍历得到的权重系数组,重新计算每个样品中所述待测元素的加权谱线强度。该加权谱线强度计算模块1104可用于执行图3中所示的步骤S103和图5中所示的步骤S203及S209,具体的操作方法可参照上述步骤S101、步骤S201及步骤S209的描述。
所述第一定标曲线确定模块1106,用于根据获取到的每个样品中所述待测元素的浓度和计算得到的每个样品中所述待测元素的加权谱线强度进行线性拟合,确定用于对该待测元素定量分析的第一定标曲线。该第一定标曲线确定模块1106可用于执行上述步骤S105以及步骤S205,具体的操作方法可参照上述步骤S105以及步骤S205的描述。
所述参数遍历模块1108,用于采用预设算法以预设值为步进单位基于所述初始权重系数组进行参数遍历。该参数遍历模块1108可用于执行图5中所示的步骤S207,具体的操作方法可参照上述步骤S207的描述。
所述第二定标曲线确定模块1110,用于根据每个样品中所述待测元素的浓度和上述重新计算得到的每个样品中所述待测元素的加权谱线强度进行线性拟合,确定用于对该待测元素定量分析的第二定标曲线。该第二定标曲线确定模块1110可用于执行图5中所示的步骤S211,具体的操作方法可参照上述步骤S211的描述。
所述决定系数计算模块1112,用于计算所述第一定标曲线和每个第二定标曲线的决定系数。该决定系数计算模块1112可用于执行图5中所示的步骤S213,具体的操作方法可参照上述步骤S213的描述。
所述目标定标曲线确定模块1114,选取决定系数最大的定标曲线作为所述待测元素的目标定标曲线。该目标定标曲线确定模块1114可用于执行图5中所示的步骤S215,具体的操作方法可参照上述步骤S215的描述。
所述平均谱线强度计算模块1116,用于根据上述多次获取结果,计算每条目标谱线的平均谱线强度。
所述浓度求取模块1118,用于将计算得到的所述待测样品中待测元素的加权谱线强度代入所述目标定标曲线,求出所述待测样品中待测元素的浓度。
以上各模块可以是由软件代码实现,此时,上述的各模块可存储于计算设备100的存储器内。以上各模块同样可以由硬件例如集成电路芯片实现。
需要说明的是,本说明书中的各个实施例均采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似的部分互相参见即可。
本发明实施例所提供的元素定量分析装置,其实现原理及产生的技术效果和前述方法实施例相同,为简要描述,装置实施例部分未提及之处,可参考前述方法实施例中相应内容。
在本申请所提供的几个实施例中,应该理解到,所揭露的装置和方法,也可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的,例如,附图中的流程图和框图显示了根据本发明的多个实施例的装置、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分,所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意,在有些作为替换的实现方式中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个连续的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合,可以用执行规定的功能或动作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的组合来实现。
所述功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器、随机存取存储器、磁碟或者光盘等各种可以存储程序代码的介质。需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。

Claims (10)

1.一种元素定量分析方法,其特征在于,所述方法包括:
获取定标样品集中的每个样品的待测元素的浓度和所述待测元素的至少两条目标谱线的谱线强度;
根据每条目标谱线预设的初始权重系数和每个样品中每条目标谱线的谱线强度,计算每个样品中所述待测元素的加权谱线强度;
根据获取到的每个样品中所述待测元素的浓度和计算得到的每个样品中所述待测元素的加权谱线强度进行线性拟合,确定用于对该待测元素进行定量分析的第一定标曲线。
2.根据权利要求1所述的元素定量分析方法,其特征在于,所述根据每条目标谱线预设的初始权重系数构成初始权重系数组,所述方法还包括:
采用预设算法以预设值为步进单位基于所述初始权重系数组进行参数遍历;
根据每次参数遍历得到的权重系数组,重新计算每个样品中所述待测元素的加权谱线强度;
根据每个样品中所述待测元素的浓度和上述重新计算得到的每个样品中所述待测元素的加权谱线强度进行线性拟合,确定用于对该待测元素定量分析的第二定标曲线;
计算所述第一定标曲线和每个第二定标曲线的决定系数;
选取决定系数最大的定标曲线作为所述待测元素的目标定标曲线。
3.根据权利要求2所述的元素定量分析方法,其特征在于,所述方法还包括:
获取待测样品中所述待测元素的每条所述目标谱线的谱线强度;
根据所述目标定标曲线对应的权重系数组和上述获取到的所述待测样品中待测元素的每条所述目标谱线的谱线强度,计算该待测样品中待测元素的加权谱线强度;
根据计算得到的所述待测样品中待测元素的加权谱线强度以及所述目标定标曲线,计算得到所述待测样品中待测元素的浓度并输出。
4.根据权利要求2所述的元素定量分析方法,其特征在于,所述方法还包括:
多次获取待测样品中所述待测元素的每条所述目标谱线的谱线强度;
根据上述多次获取结果,计算每条目标谱线的平均谱线强度;
根据所述目标定标曲线对应的权重系数组和上述获取到的每条目标谱线的平均谱线强度,计算该待测样品中待测元素的加权谱线强度;
根据计算得到的所述待测样品中待测元素的加权谱线强度以及所述目标定标曲线,计算得到所述待测样品中待测元素的浓度并输出。
5.根据权利要求2所述的元素定量分析方法,其特征在于,所述预设算法包括穷举法或蚂蚁算法,所述预设值包括0.01。
6.一种元素定量分析装置,其特征在于,包括:
谱线强度获取模块,用于获取每个样品中待测元素的浓度和所述待测元素的至少两条目标谱线的谱线强度;
加权谱线强度计算模块,用于根据每条目标谱线预设的初始权重系数和每个样品中每条目标谱线的谱线强度,计算每个样品中所述待测元素的加权谱线强度;
第一定标曲线确定模块,用于根据定标样品集中获取到的每个样品中所述待测元素的浓度和计算得到的每个样品中所述待测元素的加权谱线强度进行线性拟合,确定用于对该待测元素进行定量分析的第一定标曲线。
7.根据权利要求6所述的元素定量分析装置,其特征在于,所述根据每条目标谱线预设的初始权重系数构成初始权重系数组,所述装置还包括:
参数遍历模块,用于采用预设算法以预设值为步进单位基于所述初始权重系数组进行参数遍历;
所述加权谱线强度计算模块,还用于根据每次参数遍历得到的权重系数组,重新计算每个样品中所述待测元素的加权谱线强度;
第二定标曲线确定模块,用于根据定标样品集中每个样品中所述待测元素的浓度和上述重新计算得到的每个样品中所述待测元素的加权谱线强度进行线性拟合,确定用于对该待测元素定量分析的第二定标曲线;
决定系数计算模块,用于计算所述第一定标曲线和每个第二定标曲线的决定系数;
目标定标曲线确定模块,选取决定系数最大的定标曲线作为所述待测元素的目标定标曲线。
8.根据权利要求7所述的元素定量分析装置,其特征在于,所述谱线强度获取模块,还用于获取待测样品中所述待测元素的每条所述目标谱线的谱线强度;
所述加权谱线强度计算模块,还用于根据所述目标定标曲线对应的权重系数组和上述获取到的所述待测样品中待测元素的每条所述目标谱线的谱线强度,计算该待测样品中待测元素的加权谱线强度;
所述装置还包括:
浓度求取模块,用于根据计算得到的所述待测样品中待测元素的加权谱线强度以及所述目标定标曲线,计算得到所述待测样品中待测元素的浓度并输出。
9.根据权利要求7所述的元素定量分析装置,其特征在于,所述谱线强度获取模块,还用于多次获取待测样品中所述待测元素的每条所述目标谱线的谱线强度;
所述装置还包括:
平均谱线强度计算模块,用于根据上述多次获取结果,计算每条目标谱线的平均谱线强度;
所述加权谱线强度计算模块,还用于根据所述目标定标曲线对应的权重系数组和上述获取到的每条目标谱线的平均谱线强度,计算该待测样品中待测元素的加权谱线强度;
浓度求取模块,用于根据计算得到的所述待测样品中待测元素的加权谱线强度以及所述目标定标曲线,计算得到所述待测样品中待测元素的浓度并输出。
10.根据权利要求7所述的元素定量分析装置,其特征在于,所述预设算法包括穷举法或蚂蚁算法,所述预设值包括0.01。
CN201710550399.4A 2017-07-06 2017-07-06 一种元素定量分析方法及装置 Active CN107340284B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710550399.4A CN107340284B (zh) 2017-07-06 2017-07-06 一种元素定量分析方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710550399.4A CN107340284B (zh) 2017-07-06 2017-07-06 一种元素定量分析方法及装置

Publications (2)

Publication Number Publication Date
CN107340284A true CN107340284A (zh) 2017-11-10
CN107340284B CN107340284B (zh) 2018-04-13

Family

ID=60218512

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710550399.4A Active CN107340284B (zh) 2017-07-06 2017-07-06 一种元素定量分析方法及装置

Country Status (1)

Country Link
CN (1) CN107340284B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108572177A (zh) * 2018-05-21 2018-09-25 中北大学 基于发射光谱的焊接电弧污染物在线定温检测方法
CN114324301A (zh) * 2021-12-21 2022-04-12 杭州谱育科技发展有限公司 提高便携式libs设备检测精度的方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102262076A (zh) * 2011-07-26 2011-11-30 清华大学 基于谱线组合的激光诱导击穿光谱元素浓度测量方法
CN102313731A (zh) * 2010-07-09 2012-01-11 中国科学院沈阳自动化研究所 一种未知物组成元素含量在线检测方法
CN103411931A (zh) * 2013-07-30 2013-11-27 中国科学院上海技术物理研究所 基于加权多谱线标定的远程libs元素定量分析方法
US20140168645A1 (en) * 2012-12-13 2014-06-19 Gwangju Institute Of Science And Technology Quantitative analysis method for measuring target element in specimen using laser-induced plasma spectrum
US20140336971A1 (en) * 2013-05-08 2014-11-13 Gwangju Institute Of Science And Technology Quantitative analyzing method of cigs film using a laser induced breakdown spectroscopy

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102313731A (zh) * 2010-07-09 2012-01-11 中国科学院沈阳自动化研究所 一种未知物组成元素含量在线检测方法
CN102262076A (zh) * 2011-07-26 2011-11-30 清华大学 基于谱线组合的激光诱导击穿光谱元素浓度测量方法
US20140168645A1 (en) * 2012-12-13 2014-06-19 Gwangju Institute Of Science And Technology Quantitative analysis method for measuring target element in specimen using laser-induced plasma spectrum
US20140336971A1 (en) * 2013-05-08 2014-11-13 Gwangju Institute Of Science And Technology Quantitative analyzing method of cigs film using a laser induced breakdown spectroscopy
CN103411931A (zh) * 2013-07-30 2013-11-27 中国科学院上海技术物理研究所 基于加权多谱线标定的远程libs元素定量分析方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
KUOHU LI.ET AL: "Analytical-performance improvement of laserinduced breakdown spectroscopy for steel using multi-spectral-line calibration with an artificial neural network", 《ROYAL SOCIETY OF CHEMISTRY》 *
QINGDONG ZENG.ETAL: "Laser-induced breakdown spectroscopy using laser pulses delivered by optical fibers for analyzing Mn and Ti elements in pig iron", 《ROYAL SOCIETY OF CHEMISTRY》 *
谢承利: "激光诱导击穿光谱数据处理方法及在煤分析中的应用研究", 《中国博士学位论文全文数据库工程科技Ⅰ辑》 *
陈兴龙等: "一种多元非线性模型在激光诱导击穿光谱定量分析中的应用", 《光谱学与光谱分析》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108572177A (zh) * 2018-05-21 2018-09-25 中北大学 基于发射光谱的焊接电弧污染物在线定温检测方法
CN114324301A (zh) * 2021-12-21 2022-04-12 杭州谱育科技发展有限公司 提高便携式libs设备检测精度的方法
CN114324301B (zh) * 2021-12-21 2023-10-13 杭州谱育科技发展有限公司 提高便携式libs设备检测精度的方法

Also Published As

Publication number Publication date
CN107340284B (zh) 2018-04-13

Similar Documents

Publication Publication Date Title
JP5870851B2 (ja) 情報処理装置、情報処理方法、及びプログラム
JP6091493B2 (ja) 試料に存在する成分を決定するための分光装置と分光法
CN107340284B (zh) 一种元素定量分析方法及装置
CN106507681A (zh) 光谱定量方法、光谱定量装置及程序
Dörfer et al. Quantitative mineral analysis using Raman spectroscopy and chemometric techniques
JP6355137B2 (ja) 信号分析装置、信号分析方法及びコンピュータプログラム
Beattie Optimising reproducibility in low quality signals without smoothing; an alternative paradigm for signal processing
EP0428593A1 (en) DEVICE AND METHOD FOR MEASURING THE ACTIVITY OF RADIOACTIVE PATTERNS CONTAINING SEVERAL RADIOACTIVE ISOTOPES.
Dlask et al. Short-time fractal analysis of biological autoluminescence
Jarvis et al. Multiobjective evolutionary optimisation for surface-enhanced Raman scattering
JP2005526246A (ja) 蛍光光線のスペクトル微分画像を測定する方法と装置
JP4819680B2 (ja) ソースモデリングにおけるマルチチャンネル測定信号を用いた方法およびデバイス
US9773090B2 (en) Method computer program and system to analyze mass spectra
CN107132214B (zh) 一种基于多谱线加权的元素测量方法及装置
Rowlands et al. Denoising of spectra with no user input: a spline‐smoothing algorithm
CN113435115B (zh) 一种荧光光谱特征波长筛选方法、装置、计算机设备及可读储存介质
JP6280910B2 (ja) 分光システムの性能を測定するための方法
US20220252516A1 (en) Spectroscopic apparatus and methods for determining components present in a sample
US20230133180A1 (en) Modeling method and apparatus for model of tracing the origin of durians, and method for tracing the origin of durians
EP1105814A1 (en) Method of determining an intrinsic spectrum from a measured spectrum
JP6949034B2 (ja) 信号分析装置、信号分析方法、コンピュータプログラム、測定装置及び測定方法
US6813569B2 (en) Assembly and method for a correlator structure
CN108152262A (zh) 一种毛细管电泳核酸分析方法及系统
JP2023074746A (ja) データ生成方法、学習モデル生成方法、コンピュータプログラム、情報処理装置及び分析装置
CN118347987A (zh) 一种基于sers基底电磁增强模拟结果分析的表面增强拉曼光谱检测方法

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20230425

Address after: 430073 Room 2041, 1st Floor, Unit 1, Building 10, Guanggu Chuangye Street, Dongxin Road, Wuhan Donghu New Technology Development Zone, Wuhan, Hubei Province (formerly 3) (Wuhan Area, Free Trade Zone)

Patentee after: Wuhan Hongshi Technology Co.,Ltd.

Address before: No. 272, Jiaotong Avenue, Xiaonan District, Xiaogan City, Hubei Province, 432000

Patentee before: HUBEI ENGINEERING University