CN110658179B - 基于多信号叠加和伪逆法的激光拉曼气体浓度检测方法 - Google Patents

基于多信号叠加和伪逆法的激光拉曼气体浓度检测方法 Download PDF

Info

Publication number
CN110658179B
CN110658179B CN201910952887.7A CN201910952887A CN110658179B CN 110658179 B CN110658179 B CN 110658179B CN 201910952887 A CN201910952887 A CN 201910952887A CN 110658179 B CN110658179 B CN 110658179B
Authority
CN
China
Prior art keywords
gas
raman
laser
signal
pseudo
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910952887.7A
Other languages
English (en)
Other versions
CN110658179A (zh
Inventor
张雷洪
徐润初
叶华龙
张大伟
王凯民
康祎
张志晟
王朝锐
袁晓
熊锐
王阳
卞志祥
韩晓宣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Shanghai for Science and Technology
Original Assignee
University of Shanghai for Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by University of Shanghai for Science and Technology filed Critical University of Shanghai for Science and Technology
Priority to CN201910952887.7A priority Critical patent/CN110658179B/zh
Publication of CN110658179A publication Critical patent/CN110658179A/zh
Application granted granted Critical
Publication of CN110658179B publication Critical patent/CN110658179B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/65Raman scattering
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/28Investigating the spectrum
    • G01J3/44Raman spectrometry; Scattering spectrometry ; Fluorescence spectrometry

Abstract

本发明公开了一种基于多信号叠加和伪逆法的激光拉曼气体浓度检测方法,其特征在于,采用包括气体模块、激光模块、滤光模块和信号检测与处理模块的激光拉曼气体检测系统,通过透镜组和不同滤光片,让m种气体的斯托克斯散射光通过滤光模块,信号检测与处理模块检测叠加的拉曼信号并使用伪逆法重构单个气体的拉曼信号,进而计算气体浓度。本发明不改变现有的激光拉曼气体检测系统,对于低浓度气体,通过拉曼信号的叠加可以有效地增强抵抗噪声攻击的能力,从而提高对低浓度气体的浓度检测精度。

Description

基于多信号叠加和伪逆法的激光拉曼气体浓度检测方法
技术领域
本发明涉及一种气体检测方法,具体涉及一种基于多信号叠加和伪逆法的激光拉曼气体浓度检测方法,属于光学检测领域。
背景技术
随着我国经济的高速发展,石油和天然气需求日益增多,当前我国探明的石油储量约占世界总量的11%,天然气仅占0.4%,每年需进口大量石油和天然气。因此,深层油气的探勘和开采有着重要意义。在深层油气的勘探中,录井技术常被用来采集和分析油气信息。用于分析气体成分和浓度的气体检测技术是录井技术的核心之一。
目前常见的气体检测方法有红外光谱法、气象色谱法、质谱与色谱结合法和激光拉曼气体检测法。红外光谱法分析速度快,但检测精度差且无法检测对称结构非极性双原子分子。气象色谱法能同时检测多种气体且精度较高,但实时性差,检测周期长。质谱与色谱结合法分析速度快,精度高且实时性好,但设备昂贵且操作复杂。不同的气体在激光的照射下会产生特定的拉曼频移和斯托克斯散射光,因此可以根据拉曼频移和斯托克斯散射光的强度实现对混合气体组分的测定。激光拉曼气体检测方法无需接触气体,可以在各种恶劣环境下进行气体检测。但气体的分子密度远低于固体和液体,且气体的散射截面小,因此斯托克斯散射光的强度较弱,对于低浓度气体的检测误差较大。
综上所述,现有的气体检测方法,都有各自的缺点,要么检测范围小,要么设备复杂,要么精度较差。
发明内容
本发明所要解决的技术问题是:现有气体检测方法检测范围小,设备复杂,精度较差等问题。
为了解决上述问题,本发明提供了一种基于多信号叠加和伪逆法的激光拉曼气体浓度检测方法,其特征在于,采用包括气体模块、激光模块、滤光模块和信号检测与处理模块的激光拉曼气体检测系统,通过透镜组和不同滤光片,让m种气体的斯托克斯散射光通过滤光模块,信号检测与处理模块检测叠加的拉曼信号并使用伪逆法重构单个气体的拉曼信号,进而计算气体浓度。
优选地,所述激光拉曼气体检测系统中的气体模块包括通过气压阀与混合气体样品连通的过滤器;激光模块包括一对反射镜,两块反射镜之间依次设有激光光源、气体反应腔/激光谐振腔、棱镜;滤光模块包括用于对气体反应腔/激光谐振腔发出的激光进行处理的滤光片;信号检测与处理模块包括用于检测滤光片处理后的光的信号检测装置、用于将信号检测装置检测的信号进行转换的A/D转换器及用于对A/D转换器转换的信号进行数据处理的计算机信号处理系统。
混合气体样品需要经过干燥和过滤后才可导入气体反应腔。所述激光模块中激光在激光谐振腔中通过棱镜去除杂散光,在两面反射镜之间来回反射将激光的能量增强百倍。所述滤光模块位于信号检测处理模块前,由透镜组和窄带滤光片组成。透镜组对斯托克斯散射光起到准直聚焦的作用,不同的窄带滤光片起到了过滤杂散光和通过特定波长的散射光的作用。所述信号检测处理模块由光信号探测器、电信号放大电路(A/D转换器)和计算机组成,该模块的主要功能是将检测到的微弱光信号转换为电信号并放大后传输至计算机进行后期处理。
更优选地,包括以下步骤:
步骤1:使用Gaussian 09W软件计算气体的拉曼频移,进而计算理论斯托克斯散射光波长,为选择合适的滤光片提供参考;
步骤2:使用激光照射激光谐振腔,激光在激光谐振腔中来回反射增强激光的能量;
步骤3:干燥和过滤后的混合气体样品通入气体反应腔;
步骤4:根据气体特征光波长和所需叠加的拉曼信号数目,选择合适的滤光片,使信号检测与处理模块获取叠加的气体拉曼信号;
步骤5:利用伪逆法重构不同气体的拉曼信号;
步骤6:基于重构的拉曼信号,根据式1
Figure BDA0002226330890000021
获取不同气体的浓度C;式1中,xcs为单组份气体的重构拉曼信号值,E100为单组分气体体积浓度为100%时的拉曼信号值。
进一步地,所述步骤5中伪逆法使用的公式为:
x=φ+y=φT(φφT)-1y 式2;
式2中,y表示叠加的气体拉曼信号值,φ表示转换矩阵,x表示待重构的单气体拉曼信号。φ+表示矩阵的伪逆,T表示矩阵的转置。
更进一步地,所述伪逆法计算过程中,进行如下改进:
将信号的叠加数从2扩充至m,构造一个多维转换矩阵,所述转换矩阵的维数大小为m×m(m-1)/2。
对于浓度过低的气体检测,双路信号叠加的拉曼信号仍然微弱,增加叠加信号数目虽然可以提高散射光强度,增加抵抗系统噪声的能力,但是对于伪逆法而言,多个气体拉曼信号相互干扰,单气体的拉曼信号数值特征不明显,从而导致伪逆法的重构精度很差,故要进行上述改进。
进一步地,将单组份气体的体积浓度为100%时的拉曼信号值通过所述激光拉曼气体检测系统测量,测量方法包括如下步骤:
使用氩气驱除气体反应腔内的余气后,向气体反应腔内通入体积浓度99.9%的待测气体和体积浓度0.1%的氩气,使气体反应腔内压强为一个标准大气压;开启激光光源,激光光源发出的光依次通过气体反应腔、滤光模块,被信号检测与处理模块接收,信号检测与处理模块计算得到该单组份气体含量为体积浓度100%时的拉曼能量值。
更进一步地,所述待测气体与氩气通入气体反应腔内的流量为200-300mL/min。
根据本发明所提供的一种基于多信号叠加和伪逆法的激光拉曼气体检测方法,无需改变现有的激光拉曼气体检测系统,通过透镜组和不同滤光片的选择,信号检测装置可以检测叠加的气体拉曼信号,使用伪逆法重构即可获得单个气体的拉曼信号值,进而获取气体浓度。因此,本发明操作快速便捷,对低浓度气体有着较好的检测精度。
附图说明
图1为本发明提供的基于多信号叠加和伪逆法的激光拉曼气体浓度检测方法的流程图;
图2为实施例中使用的激光拉曼气体检测系统的示意图。
具体实施方式
为使本发明更明显易懂,兹以优选实施例,并配合附图作详细说明如下。
实施例
图2为激光拉曼气体浓度检测系统的示意图,其由气体模块、激光模块、滤光模块和信号检测与处理模块组成。所述气体模块包括通过气压阀与混合气体样品连通的过滤器;激光模块包括一对反射镜,两块反射镜之间依次设有激光光源、气体反应腔/激光谐振腔、棱镜;滤光模块包括用于对气体反应腔/激光谐振腔发出的激光进行处理的滤光片;信号检测与处理模块包括用于检测滤光片处理后的光的信号检测装置、用于将信号检测装置检测的信号进行转换的A/D转换器及用于对A/D转换器转换的信号进行数据处理的计算机信号处理系统。所述激光光源为波长632.8nm的氦氖激光。
一种基于多信号叠加和伪逆法的激光拉曼气体浓度检测方法,其采用上述激光拉曼气体检测系统分别检测由m种气体组成的混合气体中各气体的浓度。其流程如图1所示,包括如下步骤(本实施例中对6种气体组成的混合气体进行检测):
步骤1:通过Gaussian 09W软件计算气体的拉曼频移,进而计算气体的斯托克斯散射光波长,6种气体H2、C2H6、CH4、N2、O2、CO2的拉曼频移和斯托克斯散射光波长如表1所示。
表1不同气体的拉曼频移和特征光波长
Figure BDA0002226330890000041
步骤2:打开激光光源,将632.8nm的氦氖激光通入激光谐振腔,激光在激光谐振腔中来回反射增强激光的能量;
步骤3:将经过过滤和干燥的气体体积浓度为10%的H2、10%的C2H6、1%的CH4、1%的N2、0.5%的02和0.1%的CO2以及作为背景气的77.4%的氩气配制成待测混合气体通入气体腔,进入步骤3,由于氩气不会产生拉曼信号,所以在本实施例中作为背景气,不会影响实验结果;
步骤4:根据步骤1所获的6种气体的特征光波长,选择合适的滤光片,气体的拉曼信号叠加数目由2扩充至6,构造一个多维转换矩阵,公式如式3所示:
Figure BDA0002226330890000051
式3中,y表示叠加的气体拉曼信号的测量值,φ表示信号叠加的测量矩阵,x表示待重构的单个气体的拉曼信号;
步骤5:根据上述伪逆法的公式,重构单个气体的拉曼信号;
步骤6:根据上述气体浓度的计算公式,计算气体浓度。
本实施例的各数据如表2、3所示。
表2样品气体对应的斯托克斯散射光光子数
H<sub>2</sub> C<sub>2</sub>H<sub>6</sub> CH<sub>4</sub> N<sub>2</sub> O<sub>2</sub> CO<sub>2</sub>
光子数 4432.3 1217.4 5408.9 4220.6 4668.4 7901.3
表3不同方法检测误差对比
H<sub>2</sub> C<sub>2</sub>H<sub>6</sub> CH<sub>4</sub> N<sub>2</sub> O<sub>2</sub> CO<sub>2</sub>
实施例 0.22% 0.35% 1.93% 1.07% 9.36% 5.70%
直接检测 0.62% 1.19% 7.86% 5.82% 13.10% 27.64%
由表3可知,对于体积浓度为10%的H2和C2H6,体积浓度为1%的CH4和N2,体积浓度为0.5%的O2和体积浓度为0.1%的CO2,本发明提供的一种基于多信号叠加和伪逆法的激光拉曼气体浓度检测方法均不同程度的提高了检测精度,特别是对于低浓度的气体,该方法的检测精度有很大提高。
本发明无需改变现有的激光拉曼气体检测系统,通过透镜组和不同滤光片的选择,信号检测装置可以检测叠加的气体拉曼信号,使用伪逆法重构即可获得单个气体的拉曼信号值,进而获取气体浓度。因此,本发明操作快速便捷,对低浓度气体有着较好的检测精度。

Claims (7)

1.一种基于多信号叠加和伪逆法的激光拉曼气体浓度检测方法,其特征在于,采用包括气体模块、激光模块、滤光模块和信号检测与处理模块的激光拉曼气体检测系统,通过透镜组和不同滤光片,让m种气体的斯托克斯散射光通过滤光模块,信号检测与处理模块检测叠加的拉曼信号并使用伪逆法重构单个气体的拉曼信号,进而计算气体浓度;其中,m大于2。
2.如权利要求1所述的基于多信号叠加和伪逆法的激光拉曼气体浓度检测方法,其特征在于,所述激光拉曼气体检测系统中的气体模块包括通过气压阀与混合气体样品连通的过滤器;激光模块包括一对反射镜,两块反射镜之间依次设有激光光源、气体反应腔/激光谐振腔、棱镜;滤光模块包括用于对气体反应腔/激光谐振腔发出的激光进行处理的滤光片;信号检测与处理模块包括用于检测滤光片处理后的光的信号检测装置、用于将信号检测装置检测的信号进行转换的A/D转换器及用于对A/D转换器转换的信号进行数据处理的计算机信号处理系统。
3.如权利要求2所述的基于多信号叠加和伪逆法的激光拉曼气体浓度检测方法,其特征在于,包括以下步骤:
步骤1:使用Gaussian 09W软件计算气体的拉曼频移,进而计算理论斯托克斯散射光波长,为选择合适的滤光片提供参考;
步骤2:使用激光照射激光谐振腔,激光在激光谐振腔中来回反射增强激光的能量;
步骤3:干燥和过滤后的混合气体样品通入气体反应腔;
步骤4:根据气体特征光波长和所需叠加的拉曼信号数目,选择合适的滤光片,使信号检测与处理模块获取叠加的气体拉曼信号;
步骤5:利用伪逆法重构不同气体的拉曼信号;
步骤6:基于重构的拉曼信号,根据式1:
Figure FDA0003241472840000011
获取不同气体的浓度C;式1中,xcs为单组份气体的重构拉曼信号值,E100为单组分气体体积浓度为100%时的拉曼信号值。
4.如权利要求3所述的基于多信号叠加和伪逆法的激光拉曼气体浓度检测方法,其特征在于,所述步骤5中伪逆法使用的公式为:
x=φ+y=φT(φφT)-1y 式2;
式2中,y表示叠加的气体拉曼信号值,φ表示转换矩阵,x表示待重构的单气体拉曼信号;φ+表示矩阵的伪逆,T表示矩阵的转置。
5.如权利要求4所述的基于多信号叠加和伪逆法的激光拉曼气体浓度检测方法,其特征在于,所述伪逆法计算过程中,进行如下改进:
将信号的叠加数从2扩充至m,构造一个多维转换矩阵,所述转换矩阵的维数大小为m×m(m-1)/2。
6.如权利要求3所述的基于多信号叠加和伪逆法的激光拉曼气体浓度检测方法,其特征在于,将单组份气体的体积浓度为100%时的拉曼信号值通过所述激光拉曼气体检测系统测量,测量方法包括如下步骤:
使用氩气驱除气体反应腔内的余气后,向气体反应腔内通入体积浓度99.9%的待测气体和体积浓度0.1%的氩气,使气体反应腔内压强为一个标准大气压;开启激光光源,激光光源发出的光依次通过气体反应腔、滤光模块,被信号检测与处理模块接收,信号检测与处理模块计算得到该单组份气体含量为体积浓度100%时的拉曼能量值。
7.如权利要求6所述的基于多信号叠加和伪逆法的激光拉曼气体浓度检测方法,其特征在于,所述待测气体与氩气通入气体反应腔内的流量为200-300mL/min。
CN201910952887.7A 2019-10-09 2019-10-09 基于多信号叠加和伪逆法的激光拉曼气体浓度检测方法 Active CN110658179B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910952887.7A CN110658179B (zh) 2019-10-09 2019-10-09 基于多信号叠加和伪逆法的激光拉曼气体浓度检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910952887.7A CN110658179B (zh) 2019-10-09 2019-10-09 基于多信号叠加和伪逆法的激光拉曼气体浓度检测方法

Publications (2)

Publication Number Publication Date
CN110658179A CN110658179A (zh) 2020-01-07
CN110658179B true CN110658179B (zh) 2021-11-19

Family

ID=69040242

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910952887.7A Active CN110658179B (zh) 2019-10-09 2019-10-09 基于多信号叠加和伪逆法的激光拉曼气体浓度检测方法

Country Status (1)

Country Link
CN (1) CN110658179B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
USRE34153E (en) * 1985-09-11 1992-12-29 University Of Utah Molecular gas analysis by Raman scattering in intracavity laser configuration
CN102466520A (zh) * 2010-11-11 2012-05-23 香港纺织及成衣研发中心 多光谱成像颜色测量系统及其成像信号处理方法
CN103207015A (zh) * 2013-04-16 2013-07-17 华东师范大学 一种光谱重构方法及其光谱仪装置
CN106644071A (zh) * 2016-09-30 2017-05-10 上海理工大学 基于加权测量矩阵的光谱重构方法
CN107084789A (zh) * 2017-04-01 2017-08-22 上海理工大学 基于稀疏先验的单像素探测器光谱反射率重构方法
CN107196661A (zh) * 2017-05-16 2017-09-22 中国人民解放军国防科学技术大学 基于稀疏理论的谐波组合信号非均匀欠采样盲重构方法
CN108474755A (zh) * 2015-11-20 2018-08-31 集成动态电子解决方案公司 时间压缩感测系统

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10114052C1 (de) * 2001-03-15 2002-07-25 Hertz Inst Heinrich Funkübertragungsverfahren im Innenraumbereich zur parallelen Funkübertragung von digitalen Datenteilströmen und mobiles Funkübertragungssystem
CA2445680A1 (en) * 2001-05-02 2002-11-07 Csir Spectrometry using broadand filters with overlapping spectral ranges
US9395244B2 (en) * 2010-09-23 2016-07-19 Nanolambda, Inc. Spectrum reconstruction method for miniature spectrometers

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
USRE34153E (en) * 1985-09-11 1992-12-29 University Of Utah Molecular gas analysis by Raman scattering in intracavity laser configuration
CN102466520A (zh) * 2010-11-11 2012-05-23 香港纺织及成衣研发中心 多光谱成像颜色测量系统及其成像信号处理方法
CN103207015A (zh) * 2013-04-16 2013-07-17 华东师范大学 一种光谱重构方法及其光谱仪装置
CN108474755A (zh) * 2015-11-20 2018-08-31 集成动态电子解决方案公司 时间压缩感测系统
CN106644071A (zh) * 2016-09-30 2017-05-10 上海理工大学 基于加权测量矩阵的光谱重构方法
CN107084789A (zh) * 2017-04-01 2017-08-22 上海理工大学 基于稀疏先验的单像素探测器光谱反射率重构方法
CN107196661A (zh) * 2017-05-16 2017-09-22 中国人民解放军国防科学技术大学 基于稀疏理论的谐波组合信号非均匀欠采样盲重构方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于多信号叠加和伪逆法的激光拉曼气体检测方法;张雷洪等;《应用激光》;20191015;第39卷(第5期);第866-871页 *

Also Published As

Publication number Publication date
CN110658179A (zh) 2020-01-07

Similar Documents

Publication Publication Date Title
CN109100325B (zh) 一种基于光谱吸收率二次谐波特征提取的气体浓度测量方法
CA2599535C (en) Method to reduce background noise in a spectrum
CN107462535B (zh) 一种基于高斯多峰拟合的光谱解析方法
CN108918446B (zh) 一种超低浓度二氧化硫紫外差分特征提取算法
CN105044026A (zh) 基于双光谱吸收线和波形匹配的激光甲烷浓度测量方法
CN103543132B (zh) 一种基于小波变换的煤质特性测量方法
CN101819140A (zh) 气态单质汞浓度的连续监测装置和方法
CN110414169B (zh) 一种傅立叶红外气测录井方法及其装置
CN112504983A (zh) 一种基于浊度色度补偿的硝酸盐浓度预测方法
CN101408501A (zh) 采用近红外光谱-偏最小二乘法定量检测dna碱基的方法
CN110658179B (zh) 基于多信号叠加和伪逆法的激光拉曼气体浓度检测方法
CN101943658A (zh) 主元分析理论和光谱分析技术结合对水质连续监测的方法
Zhao et al. Relationships Between DOC and CDOM Based on the Total Carbon‐Specific Fluorescence Intensities for River Waters Across China
CN102944546A (zh) 一种钻井液分析装置及其分析方法
CN112255192B (zh) 一种基于光谱反演的多组分痕量呼吸气体协同测量方法
CN109632680A (zh) 一种基于排列熵的水体中磷的检测方法
CN203894167U (zh) 一种拉曼光谱气体检测系统
CN103163100B (zh) 一种气态元素汞检测方法
CN110346346B (zh) 一种基于压缩感知关联算法的拉曼气体检测方法
CN115753661A (zh) 一种基于紫外光谱的喷气燃料水分离指数测量处理方法及系统
CN110553988B (zh) 一种基于叠加态思想的no低浓度信号提取方法
CN103398966A (zh) 光谱仪在有机溶液中测定tmc浓度的方法
Cheng et al. Identification and determination of oil pollutants based on 3-D fluorescence spectrum combined with self-weighted alternating trilinear decomposition algorithm
CN112240870B (zh) 原油含水率的测量系统及方法
CN212364040U (zh) 一种用于高精度测量气体浓度的滤光片

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant