发明名称 一种基于回归分析的煤质在线检测分析方法
摘要 一种基于回归分析的煤质在线检测分析方法,利用一组定标煤炭样品的谱线强度,通过最小二乘法对所有元素建立两组定标曲线 <mrow> <msubsup> <mi>C</mi> <mi>k</mi> <mn>0</mn> </msubsup> <mo>=</mo> <msub> <mi>a</mi> <mi>k</mi> </msub> <msubsup> <mi>I</mi> <mi>k</mi> <mn>0</mn> </msubsup> <mo>+</mo> <msub> <mi>e</mi> <mi>k</mi> </msub> </mrow>和 <mrow> <msubsup> <mi>C</mi> <mi>k</mi> <mn>0</mn> </msubsup> <mo>=</mo> <msub> <mi>A</mi> <mi>k</mi> </msub> <msubsup> <mi>I</mi> <mi>k</mi> <mn>0</mn> </msubsup> <mo>+</mo> <mover> <mi>&Sigma;</mi> <mi>n</mi> </mover> <msub> <mi>B</mi> <mi>j</mi> </msub> <msubsup> <mi>C</mi> <mi>j</mi> <mn>0</mn> </msubsup> <mo>+</mo> <msub> <mi>D</mi> <mi>k</mi> </msub> <mi>g</mi> <mrow> <mo>(</mo> <msubsup> <mi>C</mi> <mi>k</mi> <mn>0</mn> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>E</mi> <mi>k</mi> </msub> </mrow>,对于未知煤流中待测样品进行检测时,通过各元素的激光诱导等离子特征光谱谱线强度计算各元素的一次计算质量浓度,然后根据定标曲线求取待测元素的质量浓度,或者通过对碳元素进行两次修正,求取各元素的二次计算质量浓度后再根据定标曲线求取待测元素的质量浓度。该方法利用了LIBS光谱的全部信息,避免了复杂定标方法导致的误差,减小了母体效应的影响,定标速度快,精度高,能对煤炭进行在线的全元素分析,可以为用户提供实时的煤质元素数据。
申请公布号 CN101509872B 申请公布日期 2010.08.18
申请号 CN200910119884.1 申请日期 2009.03.20
申请人 清华大学 发明人 李政;王哲;冯杰
分类号 G01N21/64(2006.01)I 主分类号 G01N21/64(2006.01)I
代理机构 北京鸿元知识产权代理有限公司 11327 代理人 邸更岩
主权项 1.一种基于回归分析的煤质在线检测分析方法,其特征是该方法包括了如下步骤:1)首先使用各元素质量浓度已知的一组煤炭样品进行定标,利用安装在输煤皮带上的激光诱导等离子光谱系统对煤炭样品进行在线检测:以脉冲激光器(1)为激发光源,从激光器出射的激光经过聚焦透镜(2)聚焦后作用于煤炭样品(3)表面,在聚焦点产生等离子体,等离子体在保护气体的氛围中进行冷却,产生的辐射光信号通过采集透镜(4)被实时收集,通过光纤(5)并经过光谱仪(6)处理后转化成电信号而被计算机(7)采集,得到各元素质量浓度已知的一组煤炭样品的光谱谱线,进一步得到定标煤炭样品中各种元素的激光诱导等离子特征光谱谱线强度;2)利用各元素质量浓度已知的一组定标煤炭样品的等离子特征光谱谱线强度,建立两组定标曲线;定标曲线是通过在计算机中利用回归分析中的最小二乘法拟合方法求取的,第一组定标曲线是建立元素k质量浓度C<sub>k</sub><sup>0</sup>和元素k特征谱线强度I<sub>k</sub><sup>0</sup>的关系式,即<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>C</mi><mi>k</mi><mn>0</mn></msubsup><mo>=</mo><msub><mi>a</mi><mi>k</mi></msub><msubsup><mi>I</mi><mi>k</mi><mn>0</mn></msubsup><mo>+</mo><msub><mi>e</mi><mi>k</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>上式中,C<sub>k</sub><sup>0</sup>为定标煤样品中已知的元素k的质量浓度,I<sub>k</sub><sup>0</sup>为元素k的特征谱线强度,a<sub>k</sub>,e<sub>k</sub>为通过最小二乘法拟合确定的常数;第二组定标曲线是建立元素k质量浓度C<sub>k</sub><sup>0</sup>,元素k特征谱线强度I<sub>k</sub><sup>0</sup>与其它元素质量浓度C<sub>j</sub><sup>0</sup>之间的关系式,即<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>C</mi><mi>k</mi><mn>0</mn></msubsup><mo>=</mo><msub><mi>A</mi><mi>k</mi></msub><msubsup><mi>I</mi><mi>k</mi><mn>0</mn></msubsup><mo>+</mo><mover><mi>&Sigma;</mi><mi>n</mi></mover><msub><mi>B</mi><mi>j</mi></msub><msubsup><mi>C</mi><mi>j</mi><mn>0</mn></msubsup><mo>+</mo><msub><mi>D</mi><mi>k</mi></msub><mi>g</mi><mrow><mo>(</mo><msubsup><mi>C</mi><mi>k</mi><mn>0</mn></msubsup><mo>)</mo></mrow><msubsup><mi>C</mi><mi>k</mi><mn>0</mn></msubsup><mo>+</mo><msub><mi>E</mi><mi>k</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>上式中,C<sub>k</sub><sup>0</sup>为定标煤炭样品中已知的元素k的质量浓度,I<sub>k</sub><sup>0</sup>为元素k的特征谱线强度,C<sub>j</sub><sup>0</sup>为定标煤炭样品中已知的除元素k以外的所有其它元素的质量浓度,A<sub>k</sub>,B<sub>j</sub>,D<sub>k</sub>,E<sub>k</sub>为通过最小二乘法拟合确定的常数,g(C<sub>k</sub><sup>0</sup>)为描述自吸收效应的函数,<img file="FSB00000076275600013.GIF" wi="327" he="72" />a为大于等于1且小于10的实数,a的具体值根据实验数据选用,选用原则是使元素k的第二组定标曲线公式(2)的线性拟合优度在0.9以上,对于定标煤炭样品中所有元素,都要根据公式(1)和(2)建立两组定标曲线;3)对于各元素质量浓度未知的煤流中的待测样品进行检测,通过激光诱导等离子光谱系统得到煤流的特征光谱,进一步得到煤流中待测样品各种元素的激光诱导等离子特征光谱谱线强度;4)根据待测样品中各种元素的激光诱导等离子特征光谱谱线强度和定标曲线求取待测元素i的质量浓度:对于待测样品中所有元素,根据第一组定标曲线公式(1),得到C<sub>k</sub>=a<sub>k</sub>I<sub>k</sub>+e<sub>k</sub>    (3)上式中,C<sub>k</sub>为元素k质量浓度,I<sub>k</sub>为元素k的特征谱线强度,a<sub>k</sub>,e<sub>k</sub>为定标过程中已确定的常数;根据测得的待测样品各元素的特征谱线强度I<sub>k</sub>和公式(3),可求得各元素的一次计算质量浓度值C<sub>1</sub><sup>1</sup>,C<sub>2</sub><sup>1</sup>,...,C<sub>n</sub><sup>1</sup>,当中包括待测元素i一次计算质量浓度值;待测元素i的二次计算质量浓度值根据第二组定标曲线公式(2):<maths num="0003"><![CDATA[<math><mrow><msubsup><mi>C</mi><mi>i</mi><mn>2</mn></msubsup><mo>=</mo><msub><mi>A</mi><mi>i</mi></msub><msub><mi>I</mi><mi>i</mi></msub><mo>+</mo><mover><mi>&Sigma;</mi><mi>n</mi></mover><msub><mi>B</mi><mi>j</mi></msub><msubsup><mi>C</mi><mi>j</mi><mn>1</mn></msubsup><mo>+</mo><msub><mi>D</mi><mi>i</mi></msub><mi>g</mi><mrow><mo>(</mo><msubsup><mi>C</mi><mi>i</mi><mn>1</mn></msubsup><mo>)</mo></mrow><msubsup><mi>C</mi><mi>i</mi><mn>1</mn></msubsup><mo>+</mo><msub><mi>E</mi><mi>i</mi></msub></mrow></math>]]></maths>求得,上式中I为待测元素i的特征谱线强度,C<sub>j</sub><sup>1</sup>为待测样品中除待测元素i外所有其它元素的一次计算质量浓度值,A<sub>i</sub>,B<sub>j</sub>,D<sub>i</sub>,E<sub>i</sub>为在定标过程中已确定的常数,g(C<sub>i</sub><sup>1</sup>)的具体形式也在定标过程中已确定,待测元素i的二次计算质量浓度值作为待测元素i质量浓度的最终结果。
地址 100084 北京市100084信箱82分箱清华大学专利办公室