发明名称 提高差分吸收光谱在线监测灵敏度的方法
摘要 本发明涉及提高差分吸收光谱在线监测灵敏度的方法。对气体的标准吸收截面和测量光谱数据做同样特性的特征变化提取,基于特征变化的数据进行气体成分和浓度的计算。气体吸收光谱的特征变化是其标准吸收截面频域图中能谱集中的一点或一个范围,该特征变化包括一个慢变化点υL和一个快变化点νH。通过对气体标准吸收截面做通带为(υL,νH)的带通处理而获得的,通过对实测光谱数据做快变化点(νH)的平滑(低通)处理,到等效出射光强I′(λ)。本发明可以有效去除各种噪声和干扰对在线测量的影响,保留信号中对信噪比和检测灵敏度贡献最大的部分,无需试凑法即可找到信号处理的最佳分界点,最终提高差分吸收光谱方法的在线集成的精度和灵敏度。
申请公布号 CN101726464B 申请公布日期 2011.07.27
申请号 CN200910229102.X 申请日期 2009.12.10
申请人 天津大学;天津市先石光学技术有限公司 发明人 杜振辉;马艺闻;陈文亮;徐可欣
分类号 G01N21/31(2006.01)I;G01J3/42(2006.01)I 主分类号 G01N21/31(2006.01)I
代理机构 天津市北洋有限责任专利代理事务所 12201 代理人 王丽
主权项 1.一种提高差分吸收光谱在线监测灵敏度的方法,其特征是步骤如下:1)对所有待测目标气体的标准吸收光谱做频域变换,获得各个气体吸收光谱随波长变化的特性;2)确定每一种待测气体吸收光谱的频域图中能谱集中的一点或一个范围;3)选取包括所有待测气体吸收光谱的共同特征变化作为总的特征变化,所述总特征变化是包括一个慢变化点和一个快变化点的一个区间;总特征变化也可以由几个区间构成;4)对待测气体的标准吸收截面做特征提取,基于总特征变化的快变化点ν<sub>H</sub>和慢变化点ν<sub>L</sub>分别作为快慢分界点;即将气体吸收截面分解为三个部分:慢变化、特征变化和更快变化:        σ<sub>i</sub>(λ)=σ<sub>i0</sub>(λ)+σ<sub>ic</sub>(λ)+σ<sub>in</sub>(λ)        (4)其中σ<sub>ic</sub>(λ)为气体吸收截面中随波长的特征变化部分,σ<sub>in</sub>(λ)为气体吸收截面中随波长的更快变化部分,σ<sub>i0</sub>(λ)是气体吸收截面中随波长的慢变化部分;气体吸收截面的特征变化σ<sub>ic</sub>(λ)由气体吸收截面σ<sub>i</sub>(λ)做通带为(υ<sub>L</sub>,ν<sub>H</sub>)的带通滤波得到:σ<sub>ic</sub>(λ)=σ<sub>i</sub>(λ)·F(υ<sub>L</sub>,ν<sub>H</sub>),其中,F(υ<sub>L</sub>,ν<sub>H</sub>)采用公知的带通滤波算法;5)对实测光谱数据做特征提取,基于总特征变化的快变化点ν<sub>H</sub>和慢变化点ν<sub>L</sub>分别作为快慢分界点;对实测信号I(λ)做变化特征为ν<sub>H</sub>的平滑低通处理,得到等效出射光强I′(λ):           I′(λ)=I(λ)F(ν<sub>H</sub>)                (5)对实测信号I(λ)做变化特征为ν<sub>L</sub>的平滑低通处理,得到等效入射光强<img file="FSB00000474446200011.GIF" wi="132" he="58" /><maths num="0001"><![CDATA[<math><mrow><msubsup><mi>I</mi><mn>0</mn><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mi>&lambda;</mi><mo>)</mo></mrow><mo>=</mo><mi>I</mi><mrow><mo>(</mo><mi>&lambda;</mi><mo>)</mo></mrow><mi>F</mi><mrow><mo>(</mo><msub><mi>&nu;</mi><mi>L</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,F(ν)采用公知的平滑低通滤波算法;6)基于气体吸收截面的特征变化σ<sub>ic</sub>(λ)、等效入射光强<img file="FSB00000474446200013.GIF" wi="156" he="88" />和等效出射光强I′(λ),计算待测气体的浓度:<maths num="0002"><![CDATA[<math><mrow><mi>L&Sigma;</mi><mrow><mo>(</mo><msub><mi>&sigma;</mi><mi>ic</mi></msub><mrow><mo>(</mo><mi>&lambda;</mi><mo>)</mo></mrow><msub><mi>c</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>ln</mi><mrow><mo>(</mo><msubsup><mi>I</mi><mn>0</mn><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mi>&lambda;</mi><mo>)</mo></mrow><mo>/</mo><msup><mi>I</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>&lambda;</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>通过对式(7)选择多个不同波长可以构成线性方程组(8),使用最小二乘法对线性方程组(8)进行求解,可以同时得到待测气体中的多种成分的浓度<maths num="0003"><![CDATA[<math><mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>D</mi><mi>c</mi></msub><mrow><mo>(</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>D</mi><mi>c</mi></msub><mrow><mo>(</mo><msub><mi>&lambda;</mi><mn>2</mn></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>D</mi><mi>c</mi></msub><mrow><mo>(</mo><msub><mi>&lambda;</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>L&Sigma;</mi><mrow><mo>(</mo><msub><mi>&sigma;</mi><mi>ic</mi></msub><mrow><mo>(</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>)</mo></mrow></mrow><msub><mi>c</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><mi>L&Sigma;</mi><mrow><mo>(</mo><msub><mi>&sigma;</mi><mi>ic</mi></msub><mrow><mo>(</mo><msub><mi>&lambda;</mi><mn>2</mn></msub><mo>)</mo></mrow></mrow><msub><mi>c</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mi>L&Sigma;</mi><mrow><mo>(</mo><msub><mi>&sigma;</mi><mi>ic</mi></msub><mrow><mo>(</mo><msub><mi>&lambda;</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow></mrow><msub><mi>c</mi><mi>i</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow><mo>.</mo></mrow></math>]]></maths>
地址 300072 天津市南开区卫津路92号天津大学