tdlas气体浓度反演的状态空间方法 state space approach for gas concentration recovery of tdlas.pdfVIP

第20卷第8期 强 激 光 与 粒 子 束 V01．20，No．8 LASER 2008年8月 HIGHPOWER ANDPARTICLEBEAMS Aug．，2008 文章编号： 1001-4322(2008)08—126l—04 TDLAS气体浓度反演的状态空间方法。 周晓巍， 刘华锋 (浙江大学现代光学仪器国家重点实验室．杭州310027) 摘要：针对基于可调谐激光二极管吸收光谱技术(TDLAS)的气体浓度测量系统，建立了气体状态空 间模型，并将卡尔曼滤波算法应用于TDLAS浓度反演中．实验表明卡尔曼滤波可以很好地消除测量噪声和 模型误差对实验结果的影响．与最小二乘拟合的方法相比．在相同信噪比下反演精度可以提高1倍以上． 关键词：可调谐激光二极管吸收光谱技术I气体浓度测量；状态空问f卡尔曼滤波I最／b--乘法 中图分类号：TN911．23ITN247 文献标志码：A 随着钢铁冶金、石油化工、发电等行业的快速发展，对过程的控制显得尤为重要。过程气体浓度的在线测 量是过程控制的一种重要手段，传统的红外光谱在线气体分析仪无法实时、连续地进行测量，无法满足过程控 描特定气体的吸收峰，通过得到的吸收光谱的特性结合最小二乘法反演出气体浓度[2-73。最小二乘法的缺点是 在反演求解过程中假设系统数据(如压强、温度等)精确已知。而实际情况是系统往往会受到各种因素的扰动， 因此造成该反演算法的鲁棒性不够[1]。此外在实际应用中，由于气体对人射光的吸收很小，光强的相对变化很 小，所以测量得到的吸收光谱信号很容易受到各种噪声和误差的干扰，在低信噪比的情况下，反演得到的结果 可靠性值得商榷。为了抑制各种噪声，TDLAS技术中发展了幅度调制、扫描积分、频率调制下的二次谐波检 测法等技术[6]，但是这些技术对硬件的要求比较高，系统比较复杂。 为了解决这些矛盾，本文把气体浓度随时间的变化视为状态方程，把TDLAS中测量所得的直接吸收光谱 建模成测量方程，建立了气体状态空间模型，并利用卡尔曼滤波阻伽进行求解。仿真实验证明，状态空间方法 应用于TDLAS浓度反演，具有鲁棒性强、精度高的特点，实现了利用简易系统的高精度测量。 l 方 法 1．1 TDLAS测量气体浓度的基本原理 TDLAS的基本原理是Beer-Lambert定律[1] J。／zo=exp[一茚(T)g(u)cL] (1) 线的线强度，它表示该谱线的吸收强度，只与温度有关，单位为cm·mol～；P为待测气体的分压，单位Pa；L 在P，S(T)，g(1，)，L都已知的情况下，测量出J．和Jo就可以反演出待测气体浓度c。 1．2 TDLAS的状态空间浓度反演理论 1．2．1气体状态方程 状态空间方程由两个方程组成，一个是表示状态变化的方程，称为状态方程，另一个是对状态变化进行间 接测量的测量方程。气体状态方程就是描述气体的浓度状态随时间变化的特征方程 —dc_(t一)=声(￡)f(￡)+r(t)硼(￡) (2) Qr 噪声的转移系数。假设气体的浓度是不随时间变化的，并且对式(2)离散化后可以得到气体状态方程 c(k)=c(k—1)+w(h一1) (3) ·收稿日期：2008-01-02f修订日期：2008-07—04 基金项目：浙江省重点科技计划资助课题(2005C219) 作者简介：周晓囊(1986一)，男，本科生，光电信息工程专业，skywingT@gmail．tom． 万方数据 强 激 光 与 粒 子 柬 第20卷 其中七时刻对应的扫描频率为vt，即c(志)=f(¨)。 1．2．2气体测量方程 L+”(v)，令H(D=pS(T)g(v)，于是得到