发明名称 一种基于主导因素结合偏最小二乘法的煤质在线检测方法
摘要 一种基于主导因素结合偏最小二乘法的煤质在线检测方法,利用一组定标煤炭样品的谱线强度,先考虑自吸收效应,通过最小二乘法建立待测元素和其对应特征谱线的关系式,然后考虑元素互干扰效应,利用其他元素特征谱线强度修正残差;最后通过PLS方法修正前两步所得的残差,最终所得的关系式为该方法把LIBS光谱的传统定量方法和PLS结合起来,在一定程度上改进了PLS不能处理非线性效应和预测精度差的缺点。该方法可以利用主要特征谱线解释了元素浓度的定量信息,可以引入非线性形式,而且在利用谱线中的大量信息去修正定标曲线拟合的残差的同时又避免引入过多的噪音信号,减小了母体效应的影响,能达到较高的测量精度。
申请公布号 CN102004097A 申请公布日期 2011.04.06
申请号 CN201010536981.3 申请日期 2010.11.09
申请人 清华大学 发明人 李政;王哲;冯杰
分类号 G01N21/71(2006.01)I 主分类号 G01N21/71(2006.01)I
代理机构 北京鸿元知识产权代理有限公司 11327 代理人 邸更岩
主权项 1.一种基于主导因素结合偏最小二乘法的煤质在线检测方法,其特征是该方法包括了如下步骤:1)首先使用各元素质量浓度已知的一组煤炭样品作为定标样品,利用安装在输煤皮带上方的激光诱导等离子光谱系统对定标样品进行检测,得到该组定标样品的光谱谱线,即得到了定标样品中各种元素的激光诱导等离子特征光谱谱线强度;2)利用各元素质量浓度已知的一组定标样品的等离子特征光谱谱线强度,建立定标曲线;定标曲线是通过在计算机中利用回归分析中的最小二乘法拟合方法求取的,第一步是建立定标样品中某种元素x的质量浓度<img file="FDA0000031337070000011.GIF" wi="51" he="56" />和该元素的特征谱线强度<img file="FDA0000031337070000012.GIF" wi="40" he="56" />的关系式,即<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>C</mi><mi>x</mi><mn>0</mn></msubsup><mo>=</mo><mi>f</mi><mrow><mo>(</mo><msubsup><mi>I</mi><mi>x</mi><mn>0</mn></msubsup><mo>)</mo></mrow><mo>=</mo><msub><mi>C</mi><mn>0</mn></msub><mi>ln</mi><mrow><mo>(</mo><mfrac><mrow><mi>b</mi><msub><mi>C</mi><mn>0</mn></msub></mrow><mrow><mi>a</mi><mo>+</mo><mi>b</mi><msub><mi>C</mi><mn>0</mn></msub><mo>-</mo><msubsup><mi>I</mi><mi>x</mi><mn>0</mn></msubsup></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>上式中,<img file="FDA0000031337070000014.GIF" wi="51" he="56" />为定标样品中已知的某种元素x的质量浓度,<img file="FDA0000031337070000015.GIF" wi="40" he="56" />为元素x的特征谱线强度,<img file="FDA0000031337070000016.GIF" wi="115" he="56" />为描述自吸收效应的函数;C<sub>0</sub>为饱和质量浓度,根据拟合结果确定,C<sub>0</sub>取值范围为0.4-0.7,a、b为通过最小二乘法拟合确定的常数;3)以定标样品中元素x的质量浓度减去拟合公式<img file="FDA0000031337070000017.GIF" wi="115" he="56" />的计算结果,得到第一步拟合的残差<img file="FDA0000031337070000018.GIF" wi="76" he="56" />作为第二步修正的目标值;4)建立<img file="FDA0000031337070000019.GIF" wi="52" he="56" />与定标样品中其他元素特征谱线强度之间的关系式,即<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>E</mi><mn>1</mn><mn>0</mn></msubsup><mo>=</mo><msub><mi>g</mi><mi>l</mi></msub><mrow><mo>(</mo><msub><mi>I</mi><mi>l</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>上式中,I<sub>l</sub>是相关系数最大的元素l的特征谱线强度,具体做法为首先计算<img file="FDA00000313370700000111.GIF" wi="52" he="56" />与其他元素特征谱线强度的相关系数,然后根据相关系数的值,取相关系数最大的元素的特征谱线强度为I<sub>l</sub>;g<sub>l</sub>(I<sub>l</sub>)是描述元素互干扰效应的互干扰函数,选用二次、三次或四次多项式,原则是根据公式(2)的拟合优度选择;5)以<img file="FDA00000313370700000112.GIF" wi="52" he="56" />减去拟合公式g<sub>l</sub>(I<sub>l</sub>)的计算结果,得到第二步修正的残差<img file="FDA00000313370700000113.GIF" wi="76" he="56" />作为第三步修正的目标值,然后利用全谱区间的等离子体特征光谱谱线强度,使用偏最小二乘法修正残差<img file="FDA00000313370700000114.GIF" wi="75" he="56" />具体做法是建立<img file="FDA00000313370700000115.GIF" wi="52" he="56" />和全谱区间的谱线强度I<sub>w</sub>的关系式,如下所示:<maths num="0003"><![CDATA[<math><mrow><msubsup><mi>E</mi><mn>2</mn><mn>0</mn></msubsup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>w</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>d</mi><mi>w</mi></msub><msub><mi>I</mi><mi>w</mi></msub><mo>+</mo><msub><mi>b</mi><mi>x</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>上式中,I<sub>w</sub>是波长为w的光谱谱线强度,d<sub>w</sub>、b<sub>k</sub>为通过偏最小二乘法回归拟合确定的常数;根据公式(1)、(2)和(3)最终建立的定标样品中的元素x的定标曲线关系式如下所示:<maths num="0004"><![CDATA[<math><mrow><msubsup><mi>C</mi><mi>x</mi><mn>0</mn></msubsup><mo>=</mo><mi>f</mi><mrow><mo>(</mo><msubsup><mi>I</mi><mi>x</mi><mn>0</mn></msubsup><mo>)</mo></mrow><mo>+</mo><msub><mi>g</mi><mi>l</mi></msub><mrow><mo>(</mo><msub><mi>I</mi><mi>l</mi></msub><mo>)</mo></mrow><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>w</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>d</mi><mi>w</mi></msub><msub><mi>I</mi><mi>w</mi></msub><mo>+</mo><msub><mi>b</mi><mi>x</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>6)对于待测元素质量浓度未知的煤流中的待测样品进行在线检测,通过激光诱导等离子光谱系统得到煤流中的待测样品的特征光谱,即得到煤流中待测样品各种元素的激光诱导等离子特征光谱谱线强度;7)根据待测样品中的元素x的特征谱线强度<img file="FDA0000031337070000023.GIF" wi="40" he="56" />和待测样品中的元素l的特征谱线强度I<sub>l</sub>,以及待测样品特征光谱全谱区间的谱线强度I<sub>w</sub>,代入元素x对应的定标曲线关系式(4),求取元素x在待测样品中的质量浓度;如果要测定待测样品中其他元素的质量浓度,重复以上步骤1)到步骤6),即得到待测样品中其他元素的质量浓度。
地址 100084 北京市海淀区北京市100084信箱82分箱清华大学专利办公室