免标定调制光谱的实现方法
【专利摘要】本发明基于可调谐二极管激光吸收光谱，采用适应恶劣环境的波长调制技术，提供了一种免标定调制光谱的实现方法，该方法包含数字锁相信号解调步骤、调制光谱仿真步骤和测量参数提取步骤。其具体内容是，利用数字锁相算法对光电信号解调得到其谐波信号；通过准确的理论计算得到谐波信号的仿真数据库；结合谐波信号和仿真数据库提取流场速度、温度和组分浓度。
【专利说明】免标定调制光谱的实现方法
【技术领域】
[0001]本发明属于光学测量【技术领域】，涉及信号处理与数值仿真等方法，是针对气相分子的吸收光谱测量开发的通用、可靠测量方法。
【背景技术】
[0002]光学测量技术在科研实验、环境监测和燃烧场诊断等领域获得了越来越广泛的应用。吸收光谱在理论上需要同时测量吸收前和吸收后的光强信息才能获得分子的吸收信息。在实际应用中，往往采用参考光束或函数拟合等方法确定吸收前的光强，即光谱基线。这些方法难以获得准确的光谱基线信息。此外，调制方法虽然能明显的提高测量信号的信噪比，但是调制光谱的信号强度与调制深度、谱线线型和强度调制等多方面的因素相关，因此调制光谱的定量测量也一直存在困难。
[0003]美国斯坦福大学K.H.Lyle, R.K.Hanson等人，曾利用调制光谱来增强航空发动机进气道氧气吸收光谱测量的信噪比。由于未在理论上解决调制光谱强度与线型函数间的定量对应关系，此处仅用光谱峰的中心位置来测量速度，浓度的测量则是通过直接吸收光谱实现的。此后，G.B.Rieker从理论上分析了波长调制光谱测量的原理，给出了谐波信号的表达式，着重分析了测量可能的误差来源。研究显示利用If信号归一化2f信号之后，借助准确的理论仿真可实现任意压力、浓度和温度条件下的测量。验证了水分子和二氧化碳浓度和燃烧场温度的测量结果，显示波长调制光谱较之直接吸收光谱信噪比提高约4倍。该文献并未提出详细的实现方法，仅讨论了调制光谱测量的可行性。
[0004]南京理工大学李宁等人，采用仿真的方法得到谐波信号理论值，利用二次谐波峰谷比来消除初始光强和探测器增益对信号强度的影响。因为波峰和波谷在时间上不是同时出现的，当测量过程中初始光强波动时，波峰、波谷对应的初始光强不相同，不能消除其对信号的影响。另外，该方法必需扫描得到完整的谐波信号才能提取波峰和波谷信息，对于谱线加宽严重或者出现重叠时，往往难以实现。清华大学热能工程系车璐等人，研究了利用奇数次谐波信号重建吸收谱线型。通过对吸收线型进行傅里叶级数和泰勒级数分析，得到线型函数与奇数次谐波的关系，并在实验室条件下验证了氨气的浓度测量。该方法中，谐波信号的强度仍与初始光强相关。初始光强虽然可以预先测量，但是在干扰严重的环境下，诸多非吸收因素导致的衰减致使初始光强在测量过程中出现明显的变化，此时测量结果将出现严重的偏差。因此该方法并未真正实现免标定测量，只能在实验室条件下应用。
[0005]本发明基于可调谐半导体激光吸收光谱(tunable diode laser absorptionspectroscopy, TDLAS)技术，提供一种免标定调制光谱测量的实现方法，不仅利用了调制光谱的优势，而且从理论上解决了光谱基线提取的问题，可实现真正意义上的免标定测量。
【发明内容】
[0006]本发明提供了一种免标定调制光谱的测量方法，实现基于可调谐半导体激光吸收光谱技术的真正意义上的免标定测量。包括:[0007]步骤一、数字锁相信号解调步骤；
[0008]步骤二、调制光谱仿真步骤；
[0009]步骤三、测量参数提取步骤；
[0010]其特征在于:
[0011]步骤一的数字锁相信号解调步骤包括:(1)通过数据采集设备读取探测光电信号、采样率和数据长度L ; (2)根据调制频率f、采样率和数据长度L，生成解调用的参考信号，包括一倍频正弦和余弦信号，以及二倍频正弦和余弦信号；(3)将一倍频正弦、余弦信号和二倍频正弦、余弦信号分别与载波信号相乘，得到四路混频信号；(4)根据采样频率、低通截止频率和滤波阶次，创建低通滤波算法；(5)将得到的四路混频信号分别经过数字滤波器运算，得到一次谐波信号的X分量(Xlf)和Y分量(Ylf)，以及二次谐波的X分量(X2f)和Y分量(Y2f)，并分别计算R量；
[0012]步骤二的调制光谱仿真步骤包括:(I)根据一定条件下的谱线参数进行数值积分计算，得到线型函数；(2)根据实测条件下的激光参数和计算得到的线型函数，计算谐波信号；(3)通过计算不同谱线参数对应的谐波信号，得到仿真数据库；
[0013]步骤三的测量参数提取步骤包括:(1)提取谐波信号的多普勒频移量，并计算流场速度；(2)通过提取谐波信号的最大值来获得谐波信号的高度值；(3)通过迭代算法，对比仿真谐波信号的高度值和实验谐波信号的高度值，求解流场的温度和组分浓度。
[0014]本发明利用独特的数据处理方法和准确的数值仿真，解决了调制光谱定量测量的问题。由于算法功能齐全，简化了测量硬件的需求，具有如下优点:
[0015](I)数字锁相信号解调部分不仅替代了锁相放大器，而且拓展了其功能。数字锁相算法实现了载波与参考信号混频和低通滤波等锁相放大器的全部功能，还有硬件不能实现的功能。数字锁相算法不受调制频率的限制，原则上可以实现任意频率的信号解调；算法可以无限制的扩展通道数量，解调更高阶次的谐波信号。同时利用零相位数字滤波器算法解决了低通滤波带来的相位延迟问题。
[0016](2)实现了波长调制光谱的免标定测量。调制光谱仿真部分，利用一次谐波归信号一化二次谐波信号，消除初始光强对谐波信号强度的影响，在原理上实现了免标定测量。该部分从吸收光谱的基本方程出发，通过准确的理论模型计算得到谐波信号，为实验测量提供对比值。
[0017](3)创建了一种独特的算法，解决三维空间搜索解的问题。通过迭代来求解温度、浓度与谐波信号高度比之间的对应关系，降低了仿真计算量。
【专利附图】
【附图说明】
[0018]图1为本发明的免标定调制光谱的实现方法的流程图；
[0019]图2为本发明的线型函数和谐波信号仿真示意图；
[0020]图3为本发明的仿真谐波信号及高度比与温度的关系示意图；
[0021]图4为本发明的谐波信号多普勒频移提取算法示意图；
[0022]图5为本发明的温度、浓度迭代求解算法流程图。
具体实施方案[0023]结合附图和实施案例对本发明的免标定调制光谱的实现方法作进一步详细描述。图1给出了本发明的流程图。
[0024]本发明免标定调制光谱的实现方法，实施步骤如下:
[0025]步骤一、数字锁相信号解调
[0026]免标定调制光谱的实现方法的第一步，是利用数字锁相算法处理探测信号，共分为5个任务，实施方案如下:
[0027]①、采集测量信号
[0028]测量信号采集，由数据采集设备(如数据采集卡)读取探测光电信号、采样率和数据长度L。探测光电信号是经过调制的高频载波信号，包含了被测流场的信息。
[0029]②、生成解调信号
[0030]该任务根据调制频率f、采样率和数据长度L，生成解调用的参考信号，包括一倍频正弦和余弦信号，以及二倍频正弦和余弦信号。参考信号的采样率和数据长度与探测信号相同，信号频率等于调制频率。
[0031]③、计算混频信号
[0032]数字混频计算，相当于锁相放大器中的混频过程，即将载波信号与参考信号相乘。为了求解一次谐波信号的X分量(Xlf)和Y分量(Ylf)，以及二次谐波的X分量(X2f)和Y分量(Y2f)，需要将一倍频正弦、余弦信号和二倍频正弦、余弦信号分别与载波信号相乘，得到四路混频信号。
[0033]④、生成低通滤波器
[0034]根据采用频率、低通截止频率和滤波阶次，创建低通滤波算法，即低通滤波器。本方法采用巴特沃斯(Butterworth)滤波器，因为其幅度响应曲线在通带内具有最平坦的特性，且在通带和阻带内幅度响应曲线是单调变化的。
[0035]⑤、计算谐波信号
[0036]将得到的四路混频信号分别经过数字滤波器运算，得到一次谐波信号的X分量(Xlf)和Y分量(Ylf)，以及二次谐波的X分量(X2f)和Y分量(Y2f)。并分别计算R量。引入采用零相位数字滤波算法，将消除谐波信号引滤波而产生的相位延迟，具体做法是将混频信号滤波运算、反向、再次滤波运算、再次反向。
[0037]步骤二、调制光谱仿真
[0038]①、计算线型函数
[0039]气相分子的吸收光谱线型函数通常以Voigt线型函数描述，Voigt线型函数是Gauss线型函数和Lorentz线型函数的卷积,没有解析解,只能通过数值积分计算。
[0040]②、仿真谐波信号
[0041]谐波信号仿真是数值仿真部分的核心，其目的是根据理论计算特定压力、浓度、温度条件下，某一谱线的谐波信号。谐波信号不仅与被测环境有关，还与谱线参数和激光参数有关，建立准确的理论模型，是温度与浓度测量的基础。
[0042]激光穿过距离为L的均匀介质时，光强随传播距离呈指数衰减，可由Beer-Lambert 定律表不:
[0043]In —] = -KL = -PSu(Tmv)ZL( I )[0044]其中It和Ici分别是激光入射和透射光强，kv是光谱吸收系数，P是气体总压，X是待测气体组分浓度，Φ是线型函数，V是激光频率，Svi(T)是所用谱线Vi在温度T时的谱线强度。通常通过调制DFB激光器注入电流以改变激光的频率，同时激光器发射光强也随之改变，即方程(I)中的Itl为非定值。
[0045]在扫描谱线的同时以高频信号调制激光波长，测量其谐波信号，即波长调制光谱。激光器的注入电流受到正弦调制，激光频率可以写成:
[0046]
【权利要求】
1.一种免标定调制光谱的实现方法，包括: 步骤一、数字锁相信号解调步骤； 步骤二、调制光谱仿真步骤； 步骤三、测量参数提取步骤； 其特征在于: 步骤一的数字锁相信号解调步骤包括: (1)通过数据采集设备读取探测光电信号、采样率和数据长度L； (2)根据调制频率f、采样率和数据长度L，生成解调用的参考信号，包括一倍频正弦和余弦信号，以及二倍频正弦和余弦信号； (3)将一倍频正弦、余弦信号和二倍频正弦、余弦信号分别与载波信号相乘，得到四路混频信号； (4)根据采样频率、低通截止频率和滤波阶次，创建低通滤波算法； (5)将得到的四路混频信号分别经过数字滤波器运算，得到一次谐波信号的X分量(Xlf)和Y分量(Ylf)，以及二次谐波的X分量(X2f)和Y分量(Y2f)，并分别计算R量； 步骤二的调制光谱仿真步骤包括: (1)根据一定条件下的谱线参数进行数值积分计算，得到线型函数； (2)根据实测条件下的激光参数和计算得到的线型函数，计算谐波信号； (3)通过计算不同谱线参数对应的谐波信号，得到仿真数据库； 步骤三的测量参数提取步骤包括: (1)提取谐波信号的多普勒频移量，并计算流场速度； (2)通过提取谐波信号的最大值来获得谐波信号的高度值； (3)通过迭代算法，对比仿真谐波信号的高度值和实验谐波信号的高度值，求解流场的温度和组分浓度。
【文档编号】G01N21/39GK103852443SQ201410101717
【公开日】2014年6月11日 申请日期:2014年3月19日 优先权日:2014年3月19日
【发明者】洪延姬, 王广宇, 潘虎, 赵文涛 申请人:中国人民解放军装备学院