发明名称 一种基于超光谱的高温气体二维瞬态温度场、浓度场分布同时测量的方法
摘要 本发明公开了一种基于超光谱的高温气体二维瞬态温度场、浓度场分布同时检测的方法,属于激光吸收光谱技术领域。该方法通过在待测高温区域布置激光阵列,对待测区域进行超光谱扫描,利用水蒸气特征谱线的光谱吸收来进行区域温度场和水蒸气浓度场的测量,同时利用其他气体的特征光谱吸收,来进行气体浓度场分布的测量。测量中对待测区域进行网格离散化,通过在网格各行、各列分别布置激光束,对待测气体的特征谱线进行宽光谱扫描,利用扫描得到的光谱吸收率借助智能寻优算法实现对温度场、气体组分浓度场的反演。
申请公布号 CN103557886B 申请公布日期 2016.01.13
申请号 CN201310520260.7 申请日期 2013.10.29
申请人 东南大学 发明人 周宾;李根;李剑;肖睿
分类号 G01D21/02(2006.01)I 主分类号 G01D21/02(2006.01)I
代理机构 江苏永衡昭辉律师事务所 32250 代理人 王斌
主权项 一种基于超光谱的高温气体二维瞬态温度场、浓度场分布同时检测的方法,其特征是,该方法包括以下步骤:步骤一、根据气体组分浓度和温度的空间分辨率要求,对待测二维区域进行网格离散化,划分成M×N的网格,每个网格分别对应一个温度和浓度待测值,总的待测值数目为2×M×N;步骤二、根据上述待测值的数目,从HITRAN数据库中选择I条相应气体的特征吸收谱线,其中<img file="FDA0000800526220000013.GIF" wi="507" he="80" />并提取各谱线的中心波长λ<sub>0</sub>、参考温度下线强度S(T<sub>0</sub>)、跃迁低态能级E″、分子配分函数Q(T)的系数[a,b,c,d];步骤三、由主时钟控制函数发生器产生周期性的高频驱动信号,加载到由傅里叶域锁模激光器组成的激光器系统上,以调制激光器的输出频率,激光器系统产生的激光束通过单模光纤传输,经复用器输出M+N+2个信道的激光束,其中M条激光束实现对待测区域每行进行扫描,N条激光束对待测区域每列进行扫描,经气体吸收后由光电探测器可得到相应的M个行透射光强,记为I<sub>i,r</sub>和N个列投射光强,记为I<sub>j,r</sub>,从复用器发出的另外两条激光束,其中一条激光束直接进入光电探测器用来记录参考光强I<sub>0</sub>,另一条激光束耦合到马赫‑曾德尔干涉仪来监视波长扫描;步骤四、将采集得到的行透射光强I<sub>i,r</sub>和列透射光强I<sub>j,r</sub>(r=1,2…,I)代入下式:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>r</mi></mrow></msub><mo>=</mo><mo>-</mo><mi>l</mi><mi>n</mi><mfrac><msub><mi>I</mi><mrow><mi>i</mi><mo>,</mo><mi>r</mi></mrow></msub><msub><mi>I</mi><mn>0</mn></msub></mfrac><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>M</mi><mo>;</mo><mi>r</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>I</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000800526220000011.GIF" wi="870" he="145" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>p</mi><mrow><mi>j</mi><mo>,</mo><mi>r</mi></mrow></msub><mo>=</mo><mo>-</mo><mi>l</mi><mi>n</mi><mfrac><msub><mi>I</mi><mrow><mi>j</mi><mo>,</mo><mi>r</mi></mrow></msub><msub><mi>I</mi><mn>0</mn></msub></mfrac><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>...</mo><mo>,</mo><mi>N</mi><mo>;</mo><mi>r</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>I</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000800526220000012.GIF" wi="884" he="146" /></maths>其中,I<sub>0</sub>为参考光强,p<sub>i,r</sub>为第r条谱线下,第i行的吸收率;p<sub>j,r</sub>为第r条谱线下,第j列的吸收率;步骤五、建立待测区域的温度场T、气体浓度场X的M×N×2个方程:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>A</mi><mrow><mi>i</mi><mo>,</mo><mi>r</mi></mrow></msub><mo>=</mo><mi>P</mi><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>I</mi></munderover><mi>S</mi><mrow><mo>(</mo><msub><mi>T</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>,</mo><msub><mi>&lambda;</mi><mi>r</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>X</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&CenterDot;</mo><mi>&phi;</mi><mrow><mo>(</mo><msub><mi>&lambda;</mi><mi>k</mi></msub><mo>-</mo><msub><mi>&lambda;</mi><mi>r</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>&Delta;</mi><mi>l</mi></mrow></mtd><mtd><mrow><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>...</mo><mo>,</mo><mi>M</mi><mo>;</mo><mi>r</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>...</mo><mo>,</mo><mi>I</mi></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000800526220000021.GIF" wi="1437" he="151" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>A</mi><mrow><mi>j</mi><mo>,</mo><mi>r</mi></mrow></msub><mo>=</mo><mi>P</mi><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>I</mi></munderover><mi>S</mi><mrow><mo>(</mo><msub><mi>T</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>,</mo><msub><mi>&lambda;</mi><mi>r</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>X</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&CenterDot;</mo><mi>&phi;</mi><mrow><mo>(</mo><msub><mi>&lambda;</mi><mi>k</mi></msub><mo>-</mo><msub><mi>&lambda;</mi><mi>r</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>&Delta;</mi><mi>d</mi></mrow></mtd><mtd><mrow><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>...</mo><mo>,</mo><mi>N</mi><mo>;</mo><mi>r</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>...</mo><mo>,</mo><mi>I</mi></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000800526220000022.GIF" wi="1451" he="140" /></maths>其中,A<sub>i,r</sub>为第r条谱线下第i行的吸收率计算值;A<sub>j,r</sub>为第r条谱线下第j列的吸收率计算值,P为气体压力,S(T<sub>i,j</sub>,λ<sub>r</sub>)为线强在波长为λ<sub>r</sub>、温度为T<sub>i,j</sub>时的值,X<sub>i,j</sub>为气体摩尔分数,φ(λ<sub>k</sub>‑λ<sub>r</sub>)为第r条谱线的线型函数在波长λ<sub>k</sub>处的值,Δl、Δd分别为网格的行间距和列间距;步骤六、利用智能寻优算法以下式(5)表示的相对误差函数为目标函数进行寻优计算:设定一精度值,当目标函数值小于该精度值时计算收敛,即可求得待测区域温度场T、浓度场X的分布:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>D</mi><mrow><mo>(</mo><msup><mi>T</mi><mrow><mi>r</mi><mi>e</mi><mi>c</mi></mrow></msup><mo>,</mo><msup><mi>X</mi><mrow><mi>r</mi><mi>e</mi><mi>c</mi></mrow></msup><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>I</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mfrac><msup><mrow><mo>(</mo><msub><mi>A</mi><mrow><mi>i</mi><mo>,</mo><mi>r</mi></mrow></msub><mo>-</mo><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>r</mi></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup><mrow><msup><msub><mi>A</mi><mrow><mi>i</mi><mo>,</mo><mi>r</mi></mrow></msub><mn>2</mn></msup></mrow></mfrac><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>I</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mfrac><msup><mrow><mo>(</mo><msub><mi>A</mi><mrow><mi>j</mi><mo>,</mo><mi>r</mi></mrow></msub><mo>-</mo><msub><mi>p</mi><mrow><mi>j</mi><mo>,</mo><mi>r</mi></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup><mrow><msup><msub><mi>A</mi><mrow><mi>j</mi><mo>,</mo><mi>r</mi></mrow></msub><mn>2</mn></msup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000800526220000023.GIF" wi="1219" he="181" /></maths>其中,T<sup>rec</sup>、X<sup>rec</sup>为算法重建出的温度场、浓度场,D(T<sup>rec</sup>,T<sup>rec</sup>)表示由温度场、浓度场计算值得到吸收率值和实际测量值之间的相对误差。
地址 210096 江苏省南京市四牌楼2号