CN101509872A - 一种基于回归分析的煤质在线检测分析方法 - Google Patents

Abstract

一种基于回归分析的煤质在线检测分析方法，利用一组定标煤炭样品的谱线强度，通过最小二乘法对所有元素建立两组定标曲线C⁰ _k＝a_kI⁰ _k+e_k和C⁰ _k＝A_kI⁰ _k+B_jC⁰ _k+D_kg(C⁰ _k)+E_k，对于未知煤流中待测样品进行检测时，通过各元素的激光诱导等离子特征光谱谱线强度计算各元素的一次计算质量浓度，然后根据定标曲线求取待测元素的质量浓度，或者通过对碳元素进行两次修正，求取各元素的二次计算质量浓度后再根据定标曲线求取待测元素的质量浓度。该方法利用了LIBS光谱的全部信息，避免了复杂定标方法导致的误差，减小了母体效应的影响，定标速度快，精度高，能对煤炭进行在线的全元素分析，可以为用户提供实时的煤质元素数据。

Description

一种基于回归分析的煤质在线检测分析方法

技术领域

本发明涉及一种煤炭在线检测分析方法。具体来说，方法的基本原理是激光诱导等离子光谱技术(LIBS)，使用的是基于回归分析中最小二乘法对煤炭元素进行定量分析。

背景技术

在煤矿、煤厂和电厂等用煤单位，根据需要须及时掌控皮带输送机上的煤炭中各种成分，以便适时进行调控，这就要求能实时在线地对输送带上的煤炭进行成分分析。目前普遍采用的方法大多是离线分析。这种方法采样和制样代表性差、分析速度慢、工序繁琐等缺点，不能及时反馈煤炭的信息，不能为操作人员提供实时在线的参考数据，难以适应工业生产的需求。

目前煤质在线检测中使用的技术为X射线荧光技术，中子感生瞬发γ射线分析技术和双能γ射线透射技术。但X射线荧光技术只适合于测量原子序数大于11的元素，测量精度和灵敏度不高。中子感生瞬发γ射线分析技术存在投资大、辐射危害和放射源半衰期短的缺点。而双能γ射线投射技术最大的缺点是无法全元素分析、成本较高和有安全隐患。由于这些技术本身的缺点，所以并没有得到更广泛的应用。各用煤单位急需一种精度较高，并能实现全元素分析的煤质在线快速检测方法。

近年来，激光诱导等离子光谱技术(简称LIBS)由于具有高灵敏度、无需样品预处理和实现多元素测量等优点，成为一种新的激光分析技术，在煤质在线检测上有很大的应用潜力。可是由于该技术母体效应比较明显，直接测量物质成分时精度不高，也限制了该技术在煤质在线检测中的应用。准确的定量化测量是LIBS系统在煤质在线检测中发挥作用的前提和基础。

发明内容

本发明的目的是针对目前的煤质在线分析技术存在投资大、精度不够高或不能实现全元素分析的缺点，设计了一种基于回归分析的煤质检测方法，可在使用激光感生击穿光谱系统上运用，解决了煤质在线快速分析的问题。本发明通过运用最小二乘法拟合其他元素的影响系数，同时运用自吸收函数来考虑自吸收效应的影响，既可以实现煤质的全元素分析，又可以提高测量精度。

本发明的技术方案是：

一种基于回归分析的煤质在线检测分析方法，其特征是该方法包括了如下步骤：

1)首先使用各元素质量浓度已知的一组煤炭样品进行定标，利用安装在输煤皮带上的激光诱导等离子光谱系统对煤炭样品进行在线检测：以脉冲激光器1为激发光源，从激光器出射的激光经过聚焦透镜2聚焦后作用于煤炭样品3表面，在聚焦点产生等离子体，等离子体在保护气体的氛围中进行冷却，产生的辐射光信号通过采焦透镜4被实时收集，通过光纤5并经过光谱仪6处理后转化成电信号而被计算机7采集，得到各元素质量浓度已知的一组煤炭样品的光谱谱线，进一步得到定标煤炭样品中各种元素的激光诱导等离子特征光谱谱线强度；

2)利用各元素质量浓度已知的一组定标煤炭样品的等离子特征光谱谱线强度，建立两组定标曲线；定标曲线是通过在计算机中利用回归分析中的最小二乘法拟合方法求取的，第一组定标曲线是建立元素k质量浓度

和元素k特征谱线强度

的关系式，即

$C_k^0=a_k{I}_k^0+e_k---\left(1\right)$

上式中，

为定标样品中已知的元素k的质量浓度，

为元素k的特征谱线强度，a_k，e_k为通过最小二乘法拟合确定的常数；

第二组定标曲线是建立元素k质量浓度

元素k特征谱线强度

与其它元素质量浓度

之间的关系式，即

$C_k^0=A_k{I}_k^0+\stackrel{n}{\mathrm{\Σ}}{B}_j{C}_j^0+D_{k}g\left(C_k^0\right)C_k^0+E_k---\left(2\right)$

上式中，

为定标样品中已知的元素k的质量浓度，

为元素k的特征谱线强度，

为定标样品中已知的除元素k以外的所有其它元素的质量浓度，A_k，B_j，D_k，E_k为通过最小二乘法拟合确定的常数，

为描述自吸收效应的函数， $g\left(C_k^0\right)=a^{C_k^0}-1,$ a为大于等于1且小于10的实数，a的具体值根据实验数据选用，选用原则是使元素k的第二组定标曲线公式(2)的线性拟合优度在0.9以上，对于定标煤炭样品中所有元素，都要根据公式(1)和(2)建立两组定标曲线；

3)对于各元素质量浓度未知的煤流中的待测样品进行检测，通过激光诱导等离子光谱系统得到煤流的特征光谱，进一步得到煤流中待测样品各种元素的激光诱导等离子特征光谱谱线强度；

4)根据待测样品中各种元素的激光诱导等离子特征光谱谱线强度和定标曲线求取待测元素i的质量浓度：对于待测样品中所有元素，根据第一组定标曲线公式(1)，得到

C_k＝a_kI_k+e_k            (3)

上式中，C_k为元素k质量浓度，I_k为元素k的特征谱线强度，a_k，e_k为定标过程中已确定的常数；根据测得的待测样品各元素的特征谱线强度I_k和公式(3)，可求得各元素的一次计算质量浓度值

当中包括待测元素i一次计算质量浓度值；待测元素i的二次计算质量浓度值根据第二组定标曲线公式(2)：

$C_i^2=A_i{I}_i+\stackrel{n}{\mathrm{\Σ}}{B}_j{C}_i^1+D_{i}g\left(C_i^1\right)C_i^1+E_i$ 求得，上式中I_i为待测元素i的特征谱线强度，

为样品中除待测元素i外所有其它元素的一次计算质量浓度值，A_i，B_j，D_i，E_i为在定标过程中已确定的常数，

的具体形式也在定标过程中已确定，待测元素i的二次计算质量浓度值作为待测元素i质量浓度的最终结果。

一种基于回归分析的煤质在线检测方法，其特征是：

1)首先使用各元素质量浓度已知的一组煤炭样品进行定标，利用安装在输煤皮带上的激光诱导等离子光谱系统对煤炭样品进行在线检测：以脉冲激光器1为激发光源，从激光器出射的激光经过聚焦透镜2聚焦后作用于煤炭样品3表面，在聚焦点产生等离子体，等离子体在保护气体的氛围中进行冷却，产生的辐射光信号通过采焦透镜4被实时收集，通过光纤5并经过光谱仪6处理后转化成电信号而被计算机7采集，得到各元素质量浓度已知的一组煤炭样品的光谱谱线，进一步得到定标煤炭样品中各种元素的激光诱导等离子特征光谱谱线强度；

2)利用各元素质量浓度已知的一组定标煤炭样品的等离子特征光谱谱线强度，建立两组定标曲线；定标曲线是通过在计算机中利用回归分析中的最小二乘法拟合方法求取的，第一组定标曲线是建立元素k质量浓度

和元素k特征谱线强度

的关系式，即

$C_k^0=a_k{I}_k^0+e_k---\left(4\right)$

上式中，

为定标样品中已知的元素k的质量浓度，

为元素k的特征谱线强度，a_k，e_k为通过最小二乘法拟合确定的常数；

第二组定标曲线是建立元素k质量浓度元素k特征谱线强度

与其它元素质量浓度

之间的关系式，即

$C_k^0=A_k{I}_k^0+\stackrel{n}{\mathrm{\Σ}}{B}_j{C}_j^0+D_{k}g\left(C_k^0\right)C_k^0+E_k---\left(5\right)$

上式中，

为定标样品中已知的元素k的质量浓度，

为元素k的特征谱线强度，

为定标样品中已知的除元素k以外的所有其它元素的质量浓度，A_k，B_j，D_k，E_k为通过最小二乘法拟合确定的常数，

为描述自吸收效应的函数， $g\left(C_k^0\right)=a^{C_k^0}-1,$ a为大于等于1且小于10的实数，a的具体值根据实验数据选用，选用原则是使元素k的第二组定标曲线公式(5)的线性拟合优度在0.9以上，对于定标煤炭样品中所有元素，都要根据公式(4)和(5)建立两组定标曲线；

3)对于各元素质量浓度未知的煤流中的待测样品进行检测，通过激光诱导等离子光谱系统得到煤流的特征光谱，进一步得到煤流中待测样品各种元素的激光诱导等离子特征光谱谱线强度；

4)根据待测样品中各种元素的激光诱导等离子特征光谱谱线强度和定标曲线求取待测元素i的质量浓度：考虑到煤质中主要元素为碳，故先对碳成分的测量结果做两次修正，对于待测样品中所有元素k，根据第一组定标曲线公式(4)，得到

C_k＝a_kI_k+e_k              (6)

上式中，C_j为元素k质量浓度，I_k为元素k的特征谱线强度，a_k，e_k为定标过程中已确定的常数；根据测得的待测样品各元素的特征谱线强度I_k和公式(6)，可求得各元素的一次计算质量浓度值

然后根据第二组定标曲线公式(5)：

求得碳元素的二次计算质量浓度；接着把

和

作为已知量再次代入第二组定标曲线公式(5)：

求得各元素二次计算质量浓度值

最后根据

和

以及第二组定标曲线公式(5)：

再次对C_碳进行计算，所得结果

作为碳元素质量浓度的最终结果；

对于待测元素i，基于

和

以及第二组定标曲线公式(5)：

计算所得质量浓度值为待测元素i质量浓度的最终结果。

所述的基于回归分析的煤质在线检测分析方法的步骤1)所述的保护气体包括空气，N₂，CO₂和惰性气体。

本发明具有以下优点：

本发明综合利用了LIBS光谱的全部有用信息，考虑了其它元素和自吸收效应对待测元素谱线强度的影响，避免了复杂定标方法导致的误差，定标速度快，精度高，能对皮带上煤炭进行实时在线的全元素分析；通过本发明可以提供实时准确的煤质元素数据，指导煤的燃烧或其它利用过程，提高生产效率。

附图说明

图1是本发明回归分析方法流程示意图，

图2是本发明的激光诱导等离子光谱系统结构框图。

具体实施方式

下面结合附图对本发明作进一步的说明。

如图1、2所示。

本发明提供的一种基于回归分析的煤质在线检测分析方法，其特征是该方法包括了如下步骤：

1)首先使用各元素质量浓度已知的一组煤炭样品进行定标，利用安装在输煤皮带上的激光诱导等离子光谱系统对煤炭样品进行在线检测：以脉冲激光器1为激发光源，从激光器出射的激光经过聚焦透镜2聚焦后作用于煤炭样品3表面，在聚焦点产生等离子体，等离子体在保护气体的氛围中进行冷却，产生的辐射光信号通过采焦透镜4被实时收集，通过光纤5并经过光谱仪6处理后转化成电信号而被计算机7采集，得到各元素质量浓度已知的一组煤炭样品的光谱谱线，进一步得到定标煤炭样品中各种元素的激光诱导等离子特征光谱谱线强度，所述的保护气体包括空气、N₂、CO₂或惰性气体。

2)利用各元素质量浓度已知的一组定标煤炭样品的等离子特征光谱谱线强度，建立两组定标曲线；定标曲线是通过在计算机中利用回归分析中的最小二乘法拟合方法求取的，第一组定标曲线是建立元素k质量浓度

和元素k特征谱线强度

的关系式，即

$C_k^0=a_k{I}_k^0+e_k---\left(7\right)$

上式中，

为定标煤样品中已知的元素k的质量浓度，

为元素k的特征谱线强度，a_k，e_k为通过最小二乘法拟合确定的常数；

第二组定标曲线是建立元素k质量浓度

元素k特征谱线强度

与其它元素质量浓度

之间的关系式，即

$C_k^0=A_k{I}_k^0+\stackrel{n}{\mathrm{\Σ}}{B}_j{C}_j^0+D_{k}g\left(C_k^0\right)C_k^0+E_k---\left(8\right)$

上式中，

为定标煤炭样品中已知的元素k的质量浓度，

为元素k的特征谱线强度，

为定标煤炭样品中已知的除元素k以外的所有其它元素的质量浓度，A_k，B_j，D_k，E_k为通过最小二乘法拟合确定的常数，为描述自吸收效应的函数， $g\left(C_k^0\right)=a^{C_k^0}-1,$ a为大于等于1且小于10的实数，a的具体值根据实验数据选用，选用原则是使元素k的第二组定标曲线公式(8)的线性拟合优度在0.9以上，对于定标煤炭样品中所有元素，都要根据公式(7)和(8)建立两组定标曲线；

3)对于各元素质量浓度未知的煤流中的待测样品进行检测，通过激光诱导等离子光谱系统得到煤流的特征光谱，进一步得到煤流中待测样品各种元素的激光诱导等离子特征光谱谱线强度；

4)根据待测样品中各种元素的激光诱导等离子特征光谱谱线强度和定标曲线求取待测元素i的质量浓度：对于待测样品中所有元素，根据第一组定标曲线公式(7)，得到

C_k＝a_kI_k+e_k                 (9)

上式中，C_k为元素k质量浓度，I_k为元素k的特征谱线强度，a_k，e_k为定标过程中已确定的常数；根据测得的待测样品各元素的特征谱线强度I_k和公式(9)，可求得各元素的一次计算质量浓度值

当中包括待测元素i一次计算质量浓度值；待测元素i的二次计算质量浓度值根据第二组定标曲线公式(8)：

$C_i^2=A_i{I}_i+\stackrel{n}{\mathrm{\Σ}}{B}_j{C}_j^1+D_{i}g\left(C_i^1\right)C_i^1+E_i$ 求得，上式中I_i为待测元素i的特征谱线强度，为待测样品中除待测元素i外所有其它元素的一次计算质量浓度值，A_i，B_j，D_i，E_i为在定标过程中已确定的常数，

的具体形式也在定标过程中已确定，待测元素i的二次计算质量浓度值作为待测元素i质量浓度的最终结果。

本发明提供的一种基于回归分析的煤质在线检测方法，其特征是该方法包括了如下步骤：

1)首先使用各元素质量浓度已知的一组煤炭样品进行定标，利用安装在输煤皮带上的激光诱导等离子光谱系统对煤炭样品进行在线检测：以脉冲激光器1为激发光源，从激光器出射的激光经过聚焦透镜2聚焦后作用于煤炭样品3表面，在聚焦点产生等离子体，等离子体在保护气体的氛围中进行冷却，产生的辐射光信号通过采焦透镜4被实时收集，通过光纤5并经过光谱仪6处理后转化成电信号而被计算机7采集，得到各元素质量浓度已知的一组煤炭样品的光谱谱线，进一步得到定标煤炭样品中各种元素的激光诱导等离子特征光谱谱线强度，所述的保护气体包括空气、N₂、CO₂或惰性气体。

2)利用各元素质量浓度已知的一组定标煤炭样品的等离子特征光谱谱线强度，建立两组定标曲线；定标曲线是通过在计算机中利用回归分析中的最小二乘法拟合方法求取的，第一组定标曲线是建立元素k质量浓度

和元素k特征谱线强度

的关系式，即

$C_k^0=a_k{I}_k^0+e_k---\left(10\right)$

上式中，

为定标样品中已知的元素k的质量浓度，

为元素k的特征谱线强度，a_k，e_k为通过最小二乘法拟合确定的常数；

第二组定标曲线是建立元素k质量浓度

元素k特征谱线强度

与其它元素质量浓度

之间的关系式，即

$C_k^0=A_k{I}_k^0+\stackrel{n}{\mathrm{\Σ}}{B}_j{C}_j^0+D_{k}g\left(C_k^0\right)C_k^0+E_k---\left(11\right)$

上式中，

为定标样品中已知的元素k的质量浓度，

为元素k的特征谱线强度，

为定标样品中已知的除元素k以外的所有其它元素的质量浓度，A_k，B_j，D_k，E_k为通过最小二乘法拟合确定的常数，

为描述自吸收效应的函数， $g\left(C_k^0\right)=a^{C_k^0}-1,$ a为大于等于1且小于10的实数，a的具体值根据实验数据选用，选用原则是使元素k的第二组定标曲线公式(11)的线性拟合优度在0.9以上，对于定标煤炭样品中所有元素，都要根据公式(10)和(11)建立两组定标曲线；

3)对于各元素质量浓度未知的煤流中的待测样品进行检测，通过激光诱导等离子光谱系统得到煤流的特征光谱，进一步得到煤流中待测样品各种元素的激光诱导等离子特征光谱谱线强度；

4)根据待测样品中各种元素的激光诱导等离子特征光谱谱线强度和定标曲线求取待测元素i的质量浓度：考虑到煤质中主要元素为碳，碳元素对其它元素的特征谱线强度相干性较大，所以碳元素的质量浓度测量结果会对精度有较大影响，故先对碳成分的测量结果做两次修正，对于待测样品中所有元素，根据第一组定标曲线公式(10)，得到

C_k＝a_kI_k+e_k                (12)

上式中，C_k为元素k质量浓度，I_k为元素k的特征谱线强度，a_k，e_k为定标过程中已确定的常数；根据测得的待测样品各元素的特征谱线强度I_k和公式(12)，可求得各元素的一次计算质量浓度值

然后根据第二组定标曲线公式(11)：

求得碳元素的二次计算质量浓度；接着把

和

作为已知量再次代入第二组定标曲线公式(11)：

求得各元素二次计算质量浓度值最后根据

和以及第二组定标曲线公式(11)：

再次对C_碳进行计算，所得结果

作为碳元素质量浓度的最终结果；

对于待测元素i，基于和

以及第二组定标曲线公式(11)：

计算所得质量浓度值为待测元素i质量浓度的最终结果。由于对碳元素质量浓度的测量结果修正了两次，所以待测元素i质量浓度的最终结果会更为比直接计算更为准确。

实施例1：

一种基于回归分析的煤质在线检测分析方法，其特征是

1)首先使用各元素质量浓度已知的一组无烟煤煤炭样品进行定标，种煤炭样品的主要元素的质量浓度分别为样品1含C 92.27％，H 1.14％，O 4.72％，N 0.88％，S 0.54％；样品2含C 95.76％，H 1.16％，O 2.37％，N 0.88％，S 0.32％；样品3含C 94.90％，H 1.18％，O 2.55％，N0.72％，S 0.65％；样品4：C 91.59％，H 4.04％，O 2.59％，N 1.46％，S 0.32％，因为痕量元素质量浓度很小，为方便说明，此处不考虑痕量元素对待测元素信号的影响，当然也可根据实际测量需要考虑痕量元素。把五种煤炭样品依次放在输煤皮带上，利用安装在输煤皮带上的激光诱导等离子光谱系统对煤炭样品进行在线检测，如图1所示：以脉冲激光器1为激发光源，从激光器出射的激光经过聚焦透镜2聚焦后作用于煤炭样品3表面，在聚焦点产生等离子体，等离子体在空气的氛围中进行冷却，产生的辐射光信号通过采焦透镜4被实时收集，通过光纤5并经过光谱仪6处理后转化成电信号而被计算机7采集，得到各元素质量浓度已知的一组煤炭样品的特征光谱，进一步得到煤炭样品中各种元素的激光诱导等离子特征光谱谱线强度；以N为待测元素，其它元素只考虑含量相对较多的C、H、O和S。

2)利用定标煤炭样品中N元素的等离子光谱谱线强度

建立两组定标曲线；定标曲线是通过在计算机中利用回归分析中的最小二乘法拟合方法求取的，第一组定标曲线是建立元素N质量浓度

和元素N特征谱线强度

的关系式，即

$C_N^0=a_N{I}_N^0+e_N---\left(13\right)$

上式中，为定标样品中已知的元素N的质量浓度，

为元素N的特征谱线强度，根据实验数据利用最小二乘法拟合得a_N＝517.62，e_N＝-2.55；同理可利用定标煤炭样品中各元素的特征谱线强度，

$C_C^0=a_C{I}_C^0+e_C---\left(14\right)$

$C_H^0=a_H{I}_H^0+e_H---\left(15\right)$

$C_O^0=a_O{I}_O^0+e_O---\left(16\right)$

$C_S^0=a_S{I}_S^0+e_S---\left(17\right)$

公式(14)，(15)，(16)，(17)中为C元素的特征谱线强度，

为H元素的特征谱线强度，

为O元素的特征谱线强度，

为N元素的特征谱线强度，

为N元素的特征谱线强度，

为定标煤炭样品中N元素的质量浓度，

为定标煤炭样品中C元素的质量浓度，

为定标煤炭样品中H元素的质量浓度，

为定标煤炭样品中0的质量浓度，

为定标样品中N的质量浓度，以上均为已知量。公式(14)，(15)，(16)，(17)中各未知数a_C，e_C，a_H，e_H，a_O，e_O，a_S，e_S的值通过最小二乘法拟合所得。

第二组定标曲线是建立元素N质量浓度

元素N特征谱线强度

和其它元素质量浓度

之间的关系式，即

$C_N^0=A_N{I}_N^0+B_C^N{C}_C^0+B_H^N{C}_H^0+B_O^N{C}_O^0+B_S^N{C}_S^0+D_{N}g\left(C_N^0\right)C_N^0+E_N---\left(18\right)$

上式中

为定标煤炭样品中N的质量浓度，

为N的特征谱线强度，为定标煤炭样品中C的质量浓度，

为定标煤炭样品中H的质量浓度，

为定标煤炭样品中0的质量浓度，

为定标煤炭样品中S的质量浓度，根据实验数据，取 $g\left(C_N^0\right)=0$ 可使公式(18)拟合优度达到0.96，利用最小二乘法拟合得A_N＝586.48， $B_C^N=-4.14,$ $B_H^N=-0.32,$ $B_O^N=-0.65,$ $B_C^N=-0.48,$ E_N＝-0.24；同理根据定标煤炭样品中各元素的特征谱线强度，

$C_C^0=A_C{I}_C^0+B_H^C{C}_H^0+B_O^C{C}_O^0+B_N^C{C}_N^0+B_S^C{C}_S^0+D_{C}g\left(C_C^0\right)C_C^0+E_C---\left(19\right)$

$C_H^0=A_H{I}_H^0+B_C^H{C}_C^0+B_O^H{C}_O^0+B_N^H{C}_N^0+B_S^H{C}_S^0+E_H---\left(20\right)$

$C_O^0=A_O{I}_O^0+B_C^O{C}_C^0+B_H^O{C}_H^0+B_N^O{C}_N^0+B_S^O{C}_S^0+E_O---\left(21\right)$

$C_S^0=A_S{I}_S^0+B_C^S{C}_C^0+B_H^S{C}_H^0+B_O^S{C}_O^0+B_N^S{C}_N^0+E_S---\left(22\right)$

公式(19)，(20)，(21)，(22)中

为C元素的特征谱线强度，

为H元素的特征谱线强度，

为0元素的特征谱线强度，

为N元素的特征谱线强度，

为S元素的特征谱线强度，

为定标煤炭样品中S元素的质量浓度，

为定标煤炭样品中C元素的质量浓度，

为定标煤炭样品中H元素的质量浓度，

为定标煤炭样品中0的质量浓度，

为定标样品中N的质量浓度，以上均为已知量。公式(19)，(20)，(21)，(22)中各未知数A_C，

D_C，E_C，A_H，

E_H，A_O，

E_O，A_S，

E_S的值通过最小二乘法拟合所得。根据实验数据，取 $g\left(C_C^0\right)=e^{C_C^0}-1$ 可使公式(19)的拟合优度达到0.98，取 $g\left(C_H^0\right)=g\left(C_O^0\right)=g\left(C_S^0\right)=0,$ 可使公式(20)，(21)，(22)的拟合优度均达到0.97以上。

3)对于各元素质量浓度未知的煤流中待测样品进行检测，通过安装在皮带上的激光诱导等离子光谱系统得到煤流的特征光谱，求得煤流中待测样品C，H，O，N，S元素的激光诱导等离子光谱谱线强度分别为I_C，I_H，I_O，I_N，I_S；

4)根据待测样品中各种元素的激光诱导等离子特征光谱谱线强度和定标曲线求取待测元素S的质量浓度：对于待测样品中所有元素k，根据第一组定标曲线公式(14)，(15)，(16)，(17)，即根据 $C_C^1=a_C{I}_C+e_C,$ $C_H^1=a_H{I}_H+e_H,$ $C_O^1=a_O{I}_O+e_O,$ $C_S^1=a_S{I}_s+e_S$ 计算得到C，H，O，N元素的一次计算浓度值分别为

把这些值和I_N代入第二组定标曲线公式(18)，可求得待测样品中N元素的质量浓度C_N为1.07％。

5)利用计算机显示或输出结果。

为了验证方法的准确性，当使用各元素质量浓度分别为C 93.74％，H 3.37％，O 1.43％，N 1.09％，S 0.37％的无烟煤作为待测样品，经过检测后得到其N元素的质量浓度是1.09％，测量相对误差为1.1％，可见这种方法精度较高，可以符合生产需要。

实施例2：

1)首先使用各元素质量浓度已知的一组无烟煤煤炭样品进行定标，种煤炭样品的主要元素的质量浓度分别为样品1含C 92.27％，H 1.14％，O 4.72％，N 0.88％，S 0.54％；样品2含C 95.76％，H 1.16％，O 2.37％，N 0.88％，S 0.32％；样品3含C 94.90％，H 1.18％，O 2.55％，N0.72％，S 0.65％；样品4：C 91.59％，H 4.04％，O 2.59％，N 1.46％，S 0.32％，因为痕量元素质量浓度很小，为方便说明，此处不考虑痕量元素对待测元素信号的影响，当然也可根据实际测量需要考虑痕量元素。把五种煤炭样品依次放在输煤皮带上，利用安装在输煤皮带上的激光诱导等离子光谱系统对煤炭样品进行在线检测，如图1所示：以脉冲激光器1为激发光源，从激光器出射的激光经过聚焦透镜2聚焦后作用于煤炭样品3表面，在聚焦点产生等离子体，等离子体在空气的氛围中进行冷却，产生的辐射光信号通过采焦透镜4被实时收集，通过光纤5并经过光谱仪6处理后转化成电信号而被计算机7采集，得到各元素质量浓度已知的一组煤炭样品的特征光谱，进一步得到煤炭样品中各种元素的激光诱导等离子特征光谱谱线强度；以N为待测元素，其它元素只考虑含量相对较多的C、H、O和S。

2)利用定标煤炭样品中N元素的等离子光谱谱线强度

建立两组定标曲线；定标曲线是通过在计算机中利用回归分析中的最小二乘法拟合方法求取的，第一组定标曲线是建立元素N质量浓度

和元素N特征谱线强度

的关系式，即

$C_N^0=a_N{I}_N^0+e_N---\left(23\right)$

上式中，

为定标样品中已知的元素N的质量浓度，

为元素N的特征谱线强度，根据实验数据利用最小二乘法拟合得a_N＝517.62，e_N＝-2.55；同理可利用定标煤炭样品中各元素的特征谱线强度，

$C_C^0=a_C{I}_C^0+e_C---\left(24\right)$

$C_H^0=a_H{I}_H^0+e_H---\left(25\right)$

$C_O^0=a_O{I}_O^0+e_O---\left(26\right)$

$C_S^0=a_S{I}_S^0+e_S---\left(27\right)$

公式(24)，(25)，(26)，(27)中

为C元素的特征谱线强度，

为H元素的特征谱线强度，为O元素的特征谱线强度，

为N元素的特征谱线强度，

为N元素的特征谱线强度，

为定标煤炭样品中N元素的质量浓度，

为定标煤炭样品中C元素的质量浓度，

为定标煤炭样品中H元素的质量浓度，

为定标煤炭样品中O的质量浓度，

为定标样品中N的质量浓度，以上均为已知量。公式(24)，(25)，(26)，(27)中各未知数a_C，e_C，a_H，e_H，a_O，e_O，a_S，e_S的值通过最小二乘法拟合所得。

第二组定标曲线是建立元素N质量浓度

元素N特征谱线强度

和其它元素质量浓度之间的关系式，即

$C_N^0=A_N{I}_N^0+B_C^N{C}_C^0+B_H^N{C}_H^0+B_O^N{C}_O^0+B_S^N{C}_S^0+D_{N}g\left(C_N^0\right)C_N^0+E_N---\left(28\right)$

上式中为定标煤炭样品中N的质量浓度，

为N的特征谱线强度，

为定标煤炭样品中C的质量浓度，

为定标煤炭样品中H的质量浓度，

为定标煤炭样品中O的质量浓度，

为定标煤炭样品中S的质量浓度，根据实验数据，取 $g\left(C_N^0\right)=0$ 可使公式(28)拟合优度达到0.96，利用最小二乘法拟合得A_N＝586.48， $B_C^N=-4.14,$ $B_H^N=-0.32,$ $B_O^N=-0.65,$ $B_C^N=-0.48,$ E_N＝-0.24；同理根据定标煤炭样品中各元素的特征谱线强度，

$C_C^0=A_C{I}_C^0+B_H^C{C}_H^0+B_O^C{C}_O^0+B_N^C{C}_N^0+B_S^C{C}_S^0+D_{C}g\left(C_C^0\right)C_C^0+E_C---\left(29\right)$

$C_H^0=A_H{I}_H^0+B_C^H{C}_C^0+B_O^H{C}_O^0+B_N^H{C}_N^0+B_S^H{C}_S^0+E_H---\left(30\right)$

$C_O^0=A_O{I}_O^0+B_C^O{C}_C^0+B_H^O{C}_H^0+B_N^O{C}_N^0+B_S^O{C}_S^0+E_O---\left(31\right)$

$C_S^0=A_S{I}_S^0+B_C^S{C}_C^0+B_H^S{C}_H^0+B_O^S{C}_O^0+B_N^S{C}_N^0+E_S---\left(23\right)$

公式(29)，(30)，(31)，(32)中

为C元素的特征谱线强度，

为H元素的特征谱线强度，

为O元素的特征谱线强度，

为N元素的特征谱线强度，

为S元素的特征谱线强度，

为定标煤炭样品中S元素的质量浓度，

为定标煤炭样品中C元素的质量浓度，为定标煤炭样品中H元素的质量浓度，

为定标煤炭样品中0的质量浓度，为定标样品中N的质量浓度，以上均为已知量。公式(29)，(30)，(31)，(32)中各未知数A_C，

D_C，E_C，A_H，

E_H，A_O，

E_O，A_S，

E_S的值通过最小二乘法拟合所得。根据实验数据，取 $g\left(C_C^0\right)=e^{C_C^0}-1$ 可使公式(29)的拟合优度达到0.98，取 $g\left(C_H^0\right)=g\left(C_O^0\right)=g\left(C_S^0\right)=0,$ 可使公式(30)，(31)，(32)的拟合优度均达到0.97以上。

3)对于各元素质量浓度未知的煤流中待测样品进行检测，通过安装在皮带上的激光诱导等离子光谱系统得到煤流的特征光谱，求得煤流中待测样品C，H，O，N，S元素的激光诱导等离子光谱谱线强度分别为I_C，I_H，I_O，I_N，I_S；

4)根据待测样品中各种元素的激光诱导等离子特征光谱谱线强度和定标曲线求取待测元素S的质量浓度：对于待测样品中所有元素，根据第一组定标曲线公式(23)，(24)，(25)，(26)，(27)，即根据 $C_C^1=a_C{I}_C+e_C,$ $C_H^1=a_H{I}_H+e_H,$ $C_O^1=a_O{I}_O+e_O,$ $C_N^1=a_N{I}_N+e_N,$ $C_S^1=a_S{I}_s+e_S,$ 计算得到C，H，O，N，S元素的一次计算浓度值分别为

把以上值代入第二组定标曲线公式(29)，可求得待测样品中C元素的二次计算质量质量浓度

再根据

和公式(30)，(31)，(28)，(32)求取待测样品中各元素的二次计算质量浓度

接着把

代入(29)求取待测样品中C元素的最终计算质量浓度为

基于和

以及公式(28)计算所得的N元素质量浓度值为N元素质量浓度的最终结果。所得结果为N元素质量浓度为1.10％。

5)利用计算机显示或输出结果。

为了验证方法的准确性，当使用各元素质量浓度分别为C 93.74％，H 3.37％，O 1.43％，N 1.09％，S 0.37％的无烟煤作为待测样品，经过检测后得到其N元素的质量浓度是1.10％，测量相对误差为0.9％，可见这种方法精度较高，可以符合生产需要。

建立定标曲线时运用的最小二乘回归分析方法的原理为：

回归分析是研究一个变量关于另一个(些)变量的依赖关系的计算方法和理论，其目的在于通过后者的已知或设定值，去估计和预测前者的值。在建立定标曲线时，由于定标煤炭样品各元素质量浓度已知，把待测元素的浓度值作为被解释变量，把待测元素的特征谱线强度的函数、其它元素浓度的函数以及待测元素的自吸收函数作为解释变量，运用最小二乘法拟合曲线，最终求得定标曲线中各影响系数的值。而最小二乘法的原理就是使被解释变量的估计值与观测值在总体上最为接近，即估计值和观测值两者之差的平方和最小。

与上述方法相对应的在线检测设备(如图2所示)，包括脉冲激光器(型号可以为Nd:γAG)、聚焦透镜2、采集透镜4、光纤5、光谱仪6、计算机7，其特征是脉冲激光器1安装在聚焦透镜2的上部，聚焦透镜2位于样品3的上方，采集透镜4位于样品的侧面。样品从聚焦透镜2下部通过。采集透镜4通过光纤5与光谱仪6的输入段相连，光谱仪6的输出端与计算机7相连。计算机7在相应的软件控制下可进行计算并实现显示或打印计算结果，同时提供相应的数据接口。

Claims (4)

1.一种基于回归分析的煤质在线检测分析方法，其特征是该方法包括了如下步骤：

1)首先使用各元素质量浓度已知的一组煤炭样品进行定标，利用安装在输煤皮带上的激光诱导等离子光谱系统对煤炭样品进行在线检测：以脉冲激光器(1)为激发光源，从激光器出射的激光经过聚焦透镜(2)聚焦后作用于煤炭样品(3)表面，在聚焦点产生等离子体，等离子体在保护气体的氛围中进行冷却，产生的辐射光信号通过采焦透镜(4)被实时收集，通过光纤(5)并经过光谱仪(6)处理后转化成电信号而被计算机(7)采集，得到各元素质量浓度已知的一组煤炭样品的光谱谱线，进一步得到定标煤炭样品中各种元素的激光诱导等离子特征光谱谱线强度；

2)利用各元素质量浓度已知的一组定标煤炭样品的等离子特征光谱谱线强度，建立两组定标曲线；定标曲线是通过在计算机中利用回归分析中的最小二乘法拟合方法求取的，第一组定标曲线是建立元素k质量浓度

和元素k特征谱线强度

的关系式，即

$C_k^0=a_k{I}_k^0+e_k---\left(1\right)$

上式中，

为定标煤样品中已知的元素k的质量浓度，

为元素k的特征谱线强度，a_k，e_k为通过最小二乘法拟合确定的常数；

第二组定标曲线是建立元素k质量浓度

元素k特征谱线强度

与其它元素质量浓度

之间的关系式，即

$C_k^0=A_k{I}_k^0+\stackrel{n}{\mathrm{\Σ}}{B}_j{C}_j^0+D_{k}g\left(C_k^0\right)C_k^0+E_k---\left(2\right)$

上式中，

为定标煤炭样品中已知的元素k的质量浓度，

为元素k的特征谱线强度，

为定标煤炭样品中已知的除元素k以外的所有其它元素的质量浓度，A_k，B_j，D_k，E_k为通过最小二乘法拟合确定的常数，为描述自吸收效应的函数， $g\left(C_k^0\right)=a^{C_k^0}-1,$ a为大于等于1且小于10的实数，a的具体值根据实验数据选用，选用原则是使元素k的第二组定标曲线公式(2)的线性拟合优度在0.9以上，对于定标煤炭样品中所有元素，都要根据公式(1)和(2)建立两组定标曲线；

3)对于各元素质量浓度未知的煤流中的待测样品进行检测，通过激光诱导等离子光谱系统得到煤流的特征光谱，进一步得到煤流中待测样品各种元素的激光诱导等离子特征光谱谱线强度；

4)根据待测样品中各种元素的激光诱导等离子特征光谱谱线强度和定标曲线求取待测元素i的质量浓度：对于待测样品中所有元素，根据第一组定标曲线公式(1)，得到

C_k＝a_kI_k+e_k            (3)

上式中，C_k为元素k质量浓度，I_k为元素k的特征谱线强度，a_k，e_k为定标过程中已确定的常数；根据测得的待测样品各元素的特征谱线强度I_k和公式(3)，可求得各元素的一次计算质量浓度值

当中包括待测元素i一次计算质量浓度值；待测元素i的二次计算质量浓度值根据第二组定标曲线公式(2)：

$C_i^2=A_i{I}_i+\stackrel{n}{\mathrm{\Σ}}{B}_j{C}_j^1+D_{i}g\left(C_i^1\right)C_i^1+E_i$ 求得，上式中I_i为待测元素i的特征谱线强度，

为待测样品中除待测元素i外所有其它元素的一次计算质量浓度值，A_i，B_j，D_i，E_i为在定标过程中已确定的常数，

的具体形式也在定标过程中已确定，待测元素i的二次计算质量浓度值作为待测元素i质量浓度的最终结果。

2.根据权利要求1所述的基于回归分析的煤质在线检测分析方法，其特征在于：步骤1)中所述的保护气体包括空气、N₂、CO₂或惰性气体。

3.一种基于回归分析的煤质在线检测方法，其特征是该方法包括了如下步骤：

1)首先使用各元素质量浓度已知的一组煤炭样品进行定标，利用安装在输煤皮带上的激光诱导等离子光谱系统对煤炭样品进行在线检测：以脉冲激光器(1)为激发光源，从激光器出射的激光经过聚焦透镜(2)聚焦后作用于煤炭样品(3)表面，在聚焦点产生等离子体，等离子体在保护气体的氛围中进行冷却，产生的辐射光信号通过采焦透镜(4)被实时收集，通过光纤(5)并经过光谱仪(6)处理后转化成电信号而被计算机(7)采集，得到各元素质量浓度已知的一组煤炭样品的光谱谱线，进一步得到定标煤炭样品中各种元素的激光诱导等离子特征光谱谱线强度；

2)利用各元素质量浓度已知的一组定标煤炭样品的等离子特征光谱谱线强度，建立两组定标曲线；定标曲线是通过在计算机中利用回归分析中的最小二乘法拟合方法求取的，第一组定标曲线是建立元素k质量浓度

和元素k特征谱线强度

的关系式，即

$C_k^0=a_k{I}_k^0+e_k---\left(4\right)$

上式中，

为定标样品中已知的元素k的质量浓度，

为元素k的特征谱线强度，a_k，e_k为通过最小二乘法拟合确定的常数；

第二组定标曲线是建立元素k质量浓度

元素k特征谱线强度

与其它元素质量浓度

之间的关系式，即

$C_k^0=A_k{I}_k^0+\stackrel{n}{\mathrm{\Σ}}{B}_j{C}_j^0+D_{k}g\left(C_k^0\right)C_k^0+E_k---\left(5\right)$

上式中，

为定标样品中已知的元素k的质量浓度，

为元素k的特征谱线强度，

为定标样品中已知的除元素k以外的所有其它元素的质量浓度，A_k，B_j，D_k，E_k为通过最小二乘法拟合确定的常数，

为描述自吸收效应的函数， $g\left(C_k^0\right)=a^{C_k^0}-1,$ a为大于等于1且小于10的实数，a的具体值根据实验数据选用，选用原则是使元素k的第二组定标曲线公式(2)的线性拟合优度在0.9以上，对于定标煤炭样品中所有元素，都要根据公式(4)和(5)建立两组定标曲线；

3)对于各元素质量浓度未知的煤流中的待测样品进行检测，通过激光诱导等离子光谱系统得到煤流的特征光谱，进一步得到煤流中待测样品各种元素的激光诱导等离子特征光谱谱线强度；

4)根据待测样品中各种元素的激光诱导等离子特征光谱谱线强度和定标曲线求取待测元素i的质量浓度：考虑到煤质中主要元素为碳，故先对碳成分的测量结果做两次修正，对于待测样品中所有元素，根据第一组定标曲线公式(4)，得到

C_k＝a_kI_k+e_k          (6)

上式中，C_k为元素k质量浓度，I_k为元素k的特征谱线强度，a_k，e_k为定标过程中已确定的常数；根据测得的待测样品各元素的特征谱线强度I_k和公式(6)，可求得各元素的一次计算质量浓度值

然后根据第二组定标曲线公式(5)：

求得碳元素的二次计算质量浓度；接着把

作为已知量再次代入第二组定标曲线公式(5)：

求得各元素二次计算质量浓度值

最后根据

和

以及第二组定标曲线公式(5)：再次对C_碳进行计算，所得结果

作为碳元素质量浓度的最终结果；

对于待测元素i，基于

和

以及第二组定标曲线公式(5)：

计算所得质量浓度值为待测元素i质量浓度的最终结果。

4.根据权利要求3所述的基于回归分析的煤质在线检测分析方法，其特征在于：步骤1)所述的保护气体包括空气，N₂，CO₂和惰性气体。

Images