首页 / 专利库 / 光学 / 多普勒光谱 / 一种基于激光吸收光谱层析成像技术的气流三维速度分布测量方法

一种基于激光吸收光谱层析成像技术的气流三维速度分布测量方法

阅读:265发布:2020-05-21

专利汇可以提供一种基于激光吸收光谱层析成像技术的气流三维速度分布测量方法专利检索,专利查询,专利分析的服务。并且本 发明 提供一种基于激光吸收 光谱 层析成像 技术的气流三维速度分布测量方法,分为以下步骤:1.获取多个不同 角 度的激光吸收光谱投影数据;2.将所述被测区域离散化并计算层析成像的灵敏度矩阵;3.根据激光吸收光谱多普勒频移原理及所述不同角度的激光吸收光谱投影数据建立方程组;4.求解所述非线性方程组获得三维速度分布。本发明提出的方法可以通过非 接触 的方式测量和重建气流的三维速度分布,适用性强、可靠性高,有非常广泛的应用前景。,下面是一种基于激光吸收光谱层析成像技术的气流三维速度分布测量方法专利的具体信息内容。

1.一种基于激光吸收光谱层析成像技术的气流三维速度分布测量方法:
步骤1,获取多个不同度的激光吸收光谱投影数据:选择气流当中一种已知气体分子的一条吸收谱线,利用与所述气体分子吸收谱线相对应且频率为ν的激光束,从不同的角度分别穿过所述气流流过的三维被测区域,所述气体分子与所述激光束相互作用产生了所述气体分子吸收谱线对应的激光吸收光谱投影数据,假设所述不同角度的激光束共有M条,并以标号m区别所述不同角度的激光束,将所述第m条激光束产生的投影数据记为τm(ν);
步骤2,将所述被测区域离散化并计算层析成像的灵敏度矩阵:根据层析成像的分辨率要求,将所述三维被测区域离散地划分为N个不同网格,并以标号n区分所述不同网格,根据所述不同角度激光束与所述网格的几何关系,可以计算所述第m条激光束穿过所述第n个网格的光程lmn,根据所述网格中所述气流的温度Tn、压Pn、组分浓度Xn等状态参数并结合所述光程,可计算所述第m条激光束穿过所述第n个网格对应的灵敏度矩阵分量Amn,计算公式为:
Amn=lmn×Pn×Xn×S(Tn),  (1)
其中,S(T)是所述气体分子吸收谱线对应的线强度函数,将所述灵敏度矩阵分量按照标号m、n依次排列,可获得层析成像的灵敏度矩阵:
A=[Amn]M×N;  (2)
步骤3,根据激光吸收光谱多普勒频移原理及所述不同角度的激光吸收光谱投影数据建立方程组:在所述三维被测区域所在的平面上建立空间直角坐标系O-xyz,所述第m条激光束的方向向量可以表示为:
假设所述气流在所述第n个网格中对应于所述平面直角坐标系的x、y、z方向速度分量分别是 和 则所述气流在所述第n个网格中使所述第m条激光束产生的激光吸收光谱多普勒频移为:
C
其中,ν是所述气体分子吸收谱线的中心频率,c是光速;假设所述第n个网格中由其他因素引起的频移为Δνn,所述吸收谱线因所述第m条激光束穿过所述第n个网格而产生的线型函数为:
其中,σn、γn是所述线型函数的两个未知参数,所述线型函数 的数学形式为:
假设在实际测量中有J个不同频率的测量点,并以标号j区别所述不同频率的测量点,那么所述吸收谱线因所述第m条激光束穿过所述第n个网格而产生的线型函数在第j个频率测量点νj上的线型函数值为:
根据激光吸收光谱技术的基本原理,所述吸收谱线因所述第m条激光束在所述第j个频率测量点νj上的投影值为:
以 Δνn、σn、γn为未知数,利用(8)式所述的定量关系,可建立如下形式的非线性方程组:
步骤4,求解所述非线性方程组获得三维速度分布:所述非线性方程组(9)有J×M个非线性方程,可以利用数值方法求解,将所述非线性方程组(9)化为如下等价的非线性最优化问题求解:
其中, 是在求解所述非线性最优化问题(10)迭代过程中得到的所述第m条激光束穿过所述第n个网格而产生的线型函数在第j个频率测量点νj上线型函数值的计算值,如果采用全局最优化算法,用于计算 的未知变量将逐步收敛到所述非线性方程组(9)的解,其中包括所述速度分布的三个分量 将所述速度分量合成为矢量,即可得到所述气流在所述三维速度分布。

说明书全文

一种基于激光吸收光谱层析成像技术的气流三维速度分布测

量方法

技术领域

[0001] 本发明涉及一种基于激光吸收光谱层析成像技术的气流三维速度分布测量方法,主要内容为利用激光吸收光谱层析成像技术对气流速度在三维空间的分布进行测量和重建的方法,涉及激光吸收光谱测速技术、层析成像技术等领域。

背景技术

[0002] 激光吸收光谱测量技术是一种用于测量气体环境状态参数的重要技术,广泛应用于工业过程检测、燃烧检测、大气环境监测等各个领域。研究人员利用激光吸收光谱技术开展了单一路径上气体温度、组分浓度、气流速度等方面的研究。结合层析成像技术(亦称断层扫描成像技术),研究人员还实现了对气体温度、组分浓度等分布参数的测量和重建,然而一直未有利用激光吸收光谱对气流速度分布进行测量和重建方面的报道,本发明则为利用激光吸收光谱对气流三维速度分布进行重建提供了一种切实可行的方法。
[0003] 激光吸收光谱技术所依赖的基本物理原理是Beer-Lambert定律,该定律指出,频率为v入射光强为Iin(v)的激光经过一段长度为L的气体环境时,特定的气体分子会对光强产生吸收作用,使得出射光强发生衰减,其吸收率τ(v)为:
[0004]
[0005] 其中l是激光与物质相互作用的位置,T是温度,S(T)是线强度函数,P是压,X是产生吸收作用的气体分子对应的组分浓度, 是吸收光谱线型函数。(1)式在形式上表达了气体状态参数对路径的积分关系,吸收率τ(v)在数学上被认为是路径L上的投影。这样的关系与层析成像技术的数学形式相同,因此可以将激光吸收光谱技术与层析成像技术结合,实现对三维空间内气体状态参数分布的测量和重建。
[0006] 根据单一视激光吸收光谱速度测量原理,当激光穿过被测气流时,气流速度会使得吸收光谱线型函数 的中心波长产生多普勒频移作用,导致线型函数发生变化,线型函数的变化则进一步由(1)式的关系传递到吸收率上,引起吸收率发生变化。而(1)式同时也表明,由气流速度引起的吸收率变化也是与路径积分有关的,因此也可以结合层析成像技术对气流速度的三维分布进行测量和重建。
[0007] 基于以上背景,本发明提供了一种基于激光吸收光谱层析成像技术的气流三维速度分布测量方法,通过获得多个角度的激光吸收光谱吸收率测量值,结合层析成像技术,实现对气流三维流速分布的测量和重建。

发明内容

[0008] 针对非均匀的流动气体环境,本发明提供一种基于激光吸收光谱层析成像技术的气流三维速度分布测量方法,发明所采取的方案如下:
[0009] 步骤1,获取多个不同角度的激光吸收光谱投影数据:选择气流当中一种已知气体分子的一条吸收谱线,利用与所述气体分子吸收谱线相对应且频率为v的激光束,从不同的角度分别穿过所述气流流过的三维被测区域,所述气体分子与所述激光束相互作用产生了所述气体分子吸收谱线对应的激光吸收光谱投影数据,假设所述不同角度的激光束共有M条,并以标号m区别所述不同角度的激光束,将所述第m条激光束产生的投影数据记为τm(v);
[0010] 步骤2,将所述被测区域离散化并计算层析成像的灵敏度矩阵:根据层析成像的分辨率要求,将所述三维被测区域离散地划分为N个不同网格,并以标号n区分所述不同网格,根据所述不同角度激光束与所述网格的几何关系,可以计算所述第m条激光束穿过所述第n个网格的光程lmn,根据所述网格中所述气流的温度Tn、压力Pn、组分浓度Xn等状态参数并结合所述光程,可计算所述第m条激光束穿过所述第n个网格对应的灵敏度矩阵分量Amn,计算公式为:
[0011] Amn=lmn×Pn×Xn×S(Tn),  (2)
[0012] 其中,S(T)是所述气体分子吸收谱线对应的线强度函数,将所述灵敏度矩阵分量按照标号m、n依次排列,可获得层析成像的灵敏度矩阵:
[0013] A=[Amn]M×N;  (3)
[0014] 步骤3,根据激光吸收光谱多普勒频移原理及所述不同角度的激光吸收光谱投影数据建立方程组:在所述三维被测区域所在的平面上建立空间立体直角坐标系O-xyz,所述第m条激光束的方向向量可以表示为:
[0015]
[0016] 假设所述气流在所述第n个网格中对应于所述平面直角坐标系的x、y、z方向速度分量分别是 和 则所述气流在所述第n个网格中使所述第m条激光束产生的激光吸收光谱多普勒频移为:
[0017]
[0018] 其中,vC是所述气体分子吸收谱线的中心频率,c是光速;假设所述第n个网格中由其他因素引起的频移为Δvn,所述吸收谱线因所述第m条激光束穿过所述第n个网格而产生的线型函数为:
[0019]
[0020] 其中,σn、γn是所述线型函数的两个未知参数,所述线型函数 的数学形式为:
[0021]
[0022] 假设在实际测量中有J个不同频率的测量点,并以标号j区别所述不同频率的测量点,那么所述吸收谱线因所述第m条激光束穿过所述第n个网格而产生的线型函数在第j个频率测量点vj上的线型函数值为:
[0023]
[0024] 根据激光吸收光谱技术的基本原理,所述吸收谱线因所述第m条激光束在所述第j个频率测量点vj上的投影值为:
[0025]
[0026] 以 Δvn、σn、γn为未知数,利用(9)式所述的定量关系,可建立如下形式的非线性方程组:
[0027]
[0028] 步骤4,求解所述非线性方程组获得三维速度分布:所述非线性方程组(10)有J×M个非线性方程,可以利用数值方法求解,将所述非线性方程组(10)化为如下等价的非线性最优化问题求解:
[0029]
[0030] 其中, 是在求解所述非线性最优化问题(11)迭代过程中得到的所述第m条激光束穿过所述第n个网格而产生的线型函数在第j个频率测量点vj上线型函数值的计算值,如果采用全局最优化算法,用于计算 的未知变量将逐步收敛到所述非线性方程组(10)的解,其中包括所述速度分布的三个分量 将所述速度分量合成为矢量,即可得到所述气流在所述三维速度分布。
[0031] 利用所述气流三维速度分布测量方法,重建图3中流线所示的三维气流速度分布,得到图4中流线所示的重建结果。在图3、4中,灰度的变化表示速度的大小,单位是[米/秒]。比较两幅图像可知,利用所述气流二维速度分布测量方法得到的气流三维速度分布不论大小还是方向都与给定分布一致,说明所述气流三维速度分布测量方法是有效的。
附图说明
[0032] 图1中箭头代表覆盖立方体被测区域12个激光投影角度的方向向量;
[0033] 图2中箭头代表(1,0,0)T方向上排列的平行激光束,共计10×10=100条;
[0034] 图3以流线的方式给定了立方体内三维速度分布,其灰度代表速度的大小;
[0035] 图4为利用所述气流三维速度分布测量方法重建的三维气流速度分布。

具体实施方式

[0036] 在本实施例中,给定一个立方体区域内的状态参数分布,并将整个区域划分为5×10×10的网格,利用平行束激光吸收光谱层析成像方法进行计算,给出速度分布的重建结果。
[0037] 下面结合附图对本发明的具体实施方式作进一步说明。
[0038] 步骤1,选择中心频率在7495.5cm-1附近的H2O分子谱线进行研究,根据HITRAN数据库获得所述谱线的特征参数,计算所述谱线的线强度函数S(T)。
[0039] 步骤2,给定一个立方体区域内的温度、压力、蒸气浓度、气流速度、非多普勒频移、线型函数等,将所述立方体区域均匀地划分为5×10×10的网格,所述状态参数在每个网格中是均匀的。在所述立方体区域所在的平面上建立空间直角坐标系,并以立方体相邻三边所在的直线作为坐标系的三个坐标轴。
[0040] 以12个投影角度的平行激光束对该立方体进行探测,所述12个投影角度的方向向量如图1中所示,所述方向向量的坐标如表1所示。
[0041] 表1.12个激光投影角度的方向向量
[0042]
[0043] 每个角度有10×10=100条平行的激光束,因此整个测量过程使用的激光束共计1200条,图2中箭头展示了(1,0,0)T方向上平行激光束的排列情况。所述第m条激光束穿过所述第n个网格的光程为lmn,根据所述网格中所述气流的温度Tn、压力Pn、组分浓度Xn等状态参数并结合所述光程,可计算所述第m条激光束穿过所述第n个网格对应的灵敏度矩阵分量Amn,计算公式为:
[0044] Amn=lmn×Pn×Xn×S(Tn),  (12)
[0045] 将所述灵敏度矩阵分量按照标号m、n依次排列,可获得层析成像的灵敏度矩阵:
[0046] A=[Amn]1200×500  (13)
[0047] 步骤3,所述第m条激光束的方向向量可以表示为:
[0048]
[0049] 如图3所示,所述气流在所述第n个网格中对应于所述空间直角坐标系的x、y方向速度分量分别是 和 则所述气流在所述第n个网格中使所述第m条激光束产生的激光吸收光谱多普勒频移为:
[0050]
[0051] 其中,vC是所述气体分子吸收谱线的中心频率,c是光速;所述给定的第n个网格中非多普勒效应引起的频移为Δvn,所述吸收谱线因所述第m条激光束穿过所述第n个网格而产生的线型函数为:
[0052]
[0053] 其中,σn、γn是所述给定的线型函数的参数,所述线型函数 的数学形式为:
[0054]
[0055] 在7495.5cm-1附近的0.4cm-1范围内,选择40个不同频率的测量点,那么所述吸收谱线因所述第m条激光束穿过所述第n个网格而产生的线型函数在第j个频率测量点vj上的线型函数值为:
[0056]
[0057] 根据激光吸收光谱技术的基本原理,所述吸收谱线因所述第m条激光束在所述第j个频率测量点vj上的投影值为:
[0058]
[0059] 步骤4,在仿真计算过程中,则以 Δvn、σn、γn为未知数,根据步骤3中的推演过程,可建立如下形式的非线性方程组:
[0060]
[0061] 步骤4,所述非线性方程组(20)48 000个非线性方程,可以利用数值方法求解,多将所述非线性方程组(20)化为如下等价的非线性最优化问题求解:
[0062]
[0063] 其中, 是在求解所述非线性最优化问题(21)迭代过程中,得到的所述第m条激光束穿过所述第n个网格而产生的线型函数在第j个频率测量点vj上线型函数值的计算值,如果采用全局最优化算法,用于计算 的未知变量将逐步收敛到所述非线性方程组(20)的解,其中包括所述速度分布的三维维分量 将所述速度分布的三维分量进行矢量合成,即可得到所述气流在所述三维立方体内的速度分布,如图4所示。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈