发明名称 一种基于波长调制光谱技术的气体参数在线测量方法
摘要 一种基于波长调制光谱技术的气体参数在线测量方法,属于可调谐激光二极管吸收光谱技术领域。该方法基于波长调制光谱技术利用锁相放大器输出的奇数次谐波X轴和Y轴信号拟合出气体吸收率函数,并利用一次谐波背景信号对各次谐波X轴和Y轴信号进行归一化处理,以消除背景信号、激光强度、调制系数等因素的影响,进而通过气体吸收率函数直接测量气体的温度、浓度、压力和光谱常数,解决了目前波长调制光谱技术需要标定实验测量温度和浓度以及不能测量气体压力和光谱常数的问题,拓宽了波长调制光谱技术的应用范围。
申请公布号 CN102680020A 申请公布日期 2012.09.19
申请号 CN201210152917.4 申请日期 2012.05.16
申请人 清华大学 发明人 丁艳军;彭志敏;周佩丽
分类号 G01D21/02(2006.01)I;G01N21/39(2006.01)I 主分类号 G01D21/02(2006.01)I
代理机构 北京鸿元知识产权代理有限公司 11327 代理人 邸更岩
主权项 1.一种基于波长调制光谱技术的气体参数在线测量方法,其特征在于该方法包括如下步骤:1)根据待测气体种类,从美国高分辨率光谱数据库中选取相应的吸收光谱谱线,其中心频率为v<sub>0</sub>;2)以可调谐半导体激光器(5)为光源,调节激光控制器(4)的温度及电流,使可调谐半导体激光器(5)的输出频率稳定在中心频率v<sub>0</sub>处,并用波长计(6)进行标定和监测;3)将信号发生器(1)产生的低频锯齿波和锁相放大器(2)产生的高频正弦波经过加法器(3)叠加后输入到激光控制器(4),驱动激光器产生的激光在吸收光谱谱线频率处发生扫描和调制,则激光瞬时频率v和激光强度I<sub>0</sub>用公式(1)表示:<maths num="0001"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><mi>v</mi><mo>=</mo><mover><mi>v</mi><mo>&OverBar;</mo></mover><mo>+</mo><mi>a</mi><mi>cos</mi><mrow><mo>(</mo><mi>&omega;t</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>I</mi><mn>0</mn></msub><mo>=</mo><msub><mover><mi>I</mi><mo>&OverBar;</mo></mover><mn>0</mn></msub><mo>+</mo><mi>&Delta;</mi><mi>I</mi><mi>cos</mi><mrow><mo>(</mo><mi>&omega;t</mi><mo>+</mo><mi>&psi;</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中:a[cm<sup>-1</sup>]为频率调制幅度,定义调制系数m=a/γ,γ[cm<sup>-1</sup>]为谱线半高宽,ω为调制信号的角频率,<img file="FDA00001646228400012.GIF" wi="25" he="50" />为激光频率的平均值,<img file="FDA00001646228400013.GIF" wi="42" he="52" />为激光强度的平均值,ΔI为强度调制幅度,ψ为频率调制和强度调制之间的相位差;4)将激光器(5)发出的激光准直后直接由光电探测器(8)接收,然后分两路,一路输入数字示波器(9)中记录激光强度随时间的变化过程,另一路输入到锁相放大器(2)中进行一次谐波检测,锁相放大器(2)检测到的一次谐波背景信号S<sub>1-back</sub>输入到计算机数据采集与处理系统(10)中;5)将激光器(5)发出的激光准直后经过气体介质(7)由光电探测器(8)接收,然后分两路,一路输入数字示波器(9)中记录激光强度随时间的变化过程,另一路输入到锁相放大器(2)中检测奇数次谐波信号,锁相放大器(2)检测到的奇数次谐波X轴信号X<sub>2k-1</sub>和Y轴信号Y<sub>2k-1</sub>输入到计算机数据采集与处理系统(10)中;6)计算机数据采集与处理系统中采集到的奇数次谐波X轴信号X<sub>2k-1</sub>和Y轴信号Y<sub>2k-1</sub>代入下式,得到函数fun<sub>2k-1</sub>:<maths num="0002"><![CDATA[<math><mrow><msub><mi>fun</mi><mrow><mn>2</mn><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>X</mi><mrow><mn>2</mn><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>&CenterDot;</mo><mi>sin</mi><mi>&beta;</mi><mo>-</mo><msub><mi>Y</mi><mrow><mn>2</mn><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>&CenterDot;</mo><mi>cos</mi><mi>&beta;</mi></mrow><msub><mi>S</mi><mrow><mn>1</mn><mo>-</mo><mi>back</mi></mrow></msub></mfrac><mo>,</mo><mi>k</mi><mo>=</mo><mn>1,2</mn><mo>.</mo><mo>.</mo><mo>.</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中:β为锁相放大器参考信号与输入信号之间的相位差;7)将fun<sub>2k-1</sub>代入下式,得到函数Fun<sub>k</sub>:<maths num="0003"><![CDATA[<math><mrow><msub><mi>Fun</mi><mi>k</mi></msub><mo>=</mo><msub><mi>fun</mi><mn>1</mn></msub><mo>-</mo><msub><mi>fun</mi><mn>3</mn></msub><mo>+</mo><msub><mi>fun</mi><mn>5</mn></msub><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><msup><mrow><mo>(</mo><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>fun</mi><mrow><mn>2</mn><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><msup><mrow><mo>(</mo><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>fun</mi><mrow><mn>2</mn><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>8)将Fun<sub>k</sub>代入下式,即得到气体吸收率函数<img file="FDA00001646228400016.GIF" wi="124" he="79" /><maths num="0004"><![CDATA[<math><mrow><mi>&alpha;</mi><mrow><mo>(</mo><mover><mi>v</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mi>ln</mi><mrow><mo>(</mo><mfrac><msub><mi>Fun</mi><mi>k</mi></msub><mrow><mi>sin</mi><mi>&psi;</mi></mrow></mfrac><mi></mi><mo>)</mo></mrow><msub><mo>|</mo><mrow><mi>k</mi><mo>&RightArrow;</mo><mo>&infin;</mo></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>9)根据下式,按照常规的直接吸收法即得到待测气体温度、浓度、压力和光谱常数:<img file="FDA00001646228400022.GIF" wi="1352" he="80" />式中:P[atm]为气体总压,C为气体浓度,L[cm]为激光吸收光程,S(T)[cm<sup>-2</sup>atm<sup>-1</sup>]为谱线的线强度,<img file="FDA00001646228400023.GIF" wi="98" he="71" />为分子吸收线型函数,且<img file="FDA00001646228400024.GIF" wi="326" he="81" />
地址 100084 北京市海淀区100084信箱82分箱清华大学专利办公室