CN109085138B - 一种基于调频连续波的分布式气体传感装置和方法 - Google Patents
一种基于调频连续波的分布式气体传感装置和方法 Download PDFInfo
- Publication number
- CN109085138B CN109085138B CN201811101751.7A CN201811101751A CN109085138B CN 109085138 B CN109085138 B CN 109085138B CN 201811101751 A CN201811101751 A CN 201811101751A CN 109085138 B CN109085138 B CN 109085138B
- Authority
- CN
- China
- Prior art keywords
- gas
- light
- distributed
- distributed gas
- coupler
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 22
- 230000035559 beat frequency Effects 0.000 claims abstract description 44
- 238000010521 absorption reaction Methods 0.000 claims abstract description 42
- 230000003287 optical effect Effects 0.000 claims abstract description 37
- 230000001427 coherent effect Effects 0.000 claims abstract description 33
- 238000012360 testing method Methods 0.000 claims abstract description 33
- 238000004364 calculation method Methods 0.000 claims abstract description 21
- 239000013307 optical fiber Substances 0.000 claims description 49
- 238000001228 spectrum Methods 0.000 claims description 18
- 238000000862 absorption spectrum Methods 0.000 claims description 14
- 239000000835 fiber Substances 0.000 claims description 14
- 230000010287 polarization Effects 0.000 claims description 10
- 229910000831 Steel Inorganic materials 0.000 claims description 9
- 239000010959 steel Substances 0.000 claims description 9
- 239000004038 photonic crystal Substances 0.000 claims description 8
- 230000009466 transformation Effects 0.000 claims description 7
- 239000000523 sample Substances 0.000 claims description 5
- 239000013078 crystal Substances 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims description 3
- 238000002834 transmittance Methods 0.000 claims description 3
- 230000001186 cumulative effect Effects 0.000 claims description 2
- 238000001514 detection method Methods 0.000 abstract description 12
- 238000005259 measurement Methods 0.000 abstract description 6
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000005070 sampling Methods 0.000 abstract description 2
- 239000007789 gas Substances 0.000 description 162
- 238000012544 monitoring process Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 3
- 230000003993 interaction Effects 0.000 description 3
- 230000001105 regulatory effect Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 241000894007 species Species 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 241000282414 Homo sapiens Species 0.000 description 1
- CBENFWSGALASAD-UHFFFAOYSA-N Ozone Chemical compound [O-][O+]=O CBENFWSGALASAD-UHFFFAOYSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 235000013405 beer Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000840 electrochemical analysis Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000005305 interferometry Methods 0.000 description 1
- 210000001503 joint Anatomy 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004949 mass spectrometry Methods 0.000 description 1
- 230000020477 pH reduction Effects 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/39—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using tunable lasers
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Analytical Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Optics & Photonics (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
一种基于调频连续波的分布式气体传感装置和方法,涉及气体浓度分布式测量技术,为了解决现有气体传感方法仅适用于单点检测,且取样多、实时性差的问题。可调谐光源输出的激光分别入射至辅助干涉仪和主干涉仪;辅助干涉仪,用于产生时钟信号,并发送给数据采集卡的时钟端;主干涉仪的光路中包括参考光和测试光,测试光依次进入级联的多个分布式气体装置,分布式气体装置内装有待测气体,测试光经多级反射及气体的吸收后反射光与参考光发生拍频相干;数据采集卡,用于对主干涉仪产生的拍频相干信号进行采集,并发送给计算模块;计算模块,用于根据拍频相干信号计算气体浓度。本发明适用于测量气体浓度。
Description
技术领域
本发明涉及气体浓度分布式测量技术。
背景技术
大气环境的监测对于自然社会人类生存发展有着密切的联系。全球经济社会日益发展,科技迅猛进步,随之带来的环境问题也慢慢突显出来,温室效应、臭氧层破坏、酸化雨、光化学烟雾等,亟待开发和利用合理的科学手段对大气环境进行监测,防止危及社会与人民日常生活。地下综合管廊监测、隧道气体泄露监控、大气污染物排放管理、油气运输设施安全等场合都与气体监测息息相关。对于环境复杂、空气流动性差的密闭场合,不明来源的气体泄露对于维护工作人员的人身安全将构成威胁,除此以外,可燃性、易燃易爆性气体的实时监测更关乎到生命与财产的安全。
现有的气体传感方法主要以化学分析为主,包括质谱法,电化学分析方法,仅适用于单点检测,且取样多、实时性差。
发明内容
本发明的目的是为了解决现有气体传感方法仅适用于单点检测,且取样多、实时性差的问题,从而提供一种基于调频连续波的分布式气体传感装置和方法。
本发明所述的一种基于调频连续波的分布式气体传感装置,包括可调谐光源1、辅助干涉仪、主干涉仪、数据采集卡10和计算模块11;
可调谐光源1输出的激光分别入射至辅助干涉仪和主干涉仪;
辅助干涉仪,用于产生时钟信号,并发送给数据采集卡10的时钟端;
主干涉仪的光路中包括参考光和测试光,测试光依次进入级联的多个分布式气体装置8,分布式气体装置8内装有待测气体,测试光经多级反射及气体的吸收后反射光与参考光发生拍频相干;
数据采集卡10,用于对主干涉仪产生的拍频相干信号进行采集,并发送给计算模块11;
计算模块11,用于根据拍频相干信号计算气体浓度。
优选的是,可调谐光源1输出的激光经第一时延光纤3-1后由第一耦合器2-1分成两路光,一路光入射至辅助干涉仪,另一路光入射至主干涉仪;
辅助干涉仪包括第四耦合器2-4、第一法拉第旋转镜6-1、第二法拉第旋转镜6-2和第二时延光纤3-2;
第一耦合器2-1输出的一路光经第四耦合器2-4分成两路,一路光经第一法拉第旋转镜6-1反射后耦合进入第四耦合器2-4,另一路光经第二时延光纤3-2后入射至第二法拉第旋转镜6-2,再经第二法拉第旋转镜6-2反射后耦合进入第四耦合器2-4,第四耦合器2-4输出拍频相干信号,然后经TTL时钟发生器7产生方波信号,该方波信号即为时钟信号。
优选的是,主干涉仪包括第二耦合器2-2、第三耦合器2-3、偏振控制器4、光纤环形器5、多个分布式气体装置8和平衡探测器9;
第一耦合器2-1输出的另一路光经第二耦合器2-2分成两路光,一路光作为参考光经偏振控制器4后入射至第三耦合器2-3的一个输入端口,另一路光作为测试光由光纤环形器5的1端口进入,由光纤环形器5的2端口出射至级联的多个分布式气体装置8,在分布式气体装置8内经多级反射以及气体的吸收后反射光重新进入光纤环行器5的2端口并从3端口耦合进入第三耦合器2-3的另一个输入端口,参考光与测试光发生拍频相干并由第三耦合器2-3的2个输出端输出至平衡探测器9,平衡探测器9的输出端连接数据采集卡10的采集信号输入端。
优选的是,分布式气体装置8包括机械钢体气体容器8-1、第一光纤准直器8-2、第二光纤准直器8-3、进气口阀门8-4和出气口阀门8-5;
机械钢体气体容器8-1的两端分别设有第一光纤准直器8-2和第二光纤准直器8-3,靠近机械钢体气体容器8-1的两端分别设有进气口阀门8-4和出气口阀门8-5。
优选的是,分布式气体装置8包括空心光子晶体光纤和2个三通气阀;
2个三通气阀中各有1个口分别用于进气和出气;
空心光子晶体光纤的两端各通过1个三通气阀的另两个口与光路中的光纤耦合。
优选的是,计算模块11根据拍频相干信号计算气体浓度,具体为:
计算模块11对接收到的拍频相干信号进行傅里叶变换,实现由频域谱到时域谱的变换,根据反射峰间的光程差得到吸收光程;
选取时域谱中每个分布式气体装置8末端的信息数据段进行逆傅里叶变换,再对幅值数据作平方处理,得到气体的吸收光谱强度信息;
根据吸收光谱强度信息和吸收光程得到每一个分布式气体装置8内待测气体的浓度累计值;
将当前分布式气体装置8的浓度累计值减去其前一个分布式气体装置8的浓度累计值,得到当前分布式气体装置8内待测气体的浓度。
本发明所述的一种基于调频连续波的分布式气体传感方法,该方法包括:
步骤一、开启可调谐光源1线性扫频输出,可调谐光源1输出触发信号触发数据采集卡10作数据采集准备;
步骤二、辅助干涉仪产生时钟信号,传输至数据采集卡10,使数据采集卡10对主干涉仪的拍频信号进行等频率间隔采集;
主干涉仪的测试光由光纤环行器5的2端口导入级联的多个分布式气体装置8,在分布式气体装置8内经多级反射以及气体的吸收后反射光重新进入光纤环行器5;调节偏振控制器4使得从光纤环行器5的3端口导出的光与参考光发生拍频相干,然后由平衡探测器9传输至数据采集卡10;
步骤三、数据采集卡10采集拍频相干信号;
步骤四、计算模块11对拍频相干信号进行保存,并对拍频相干信号进行傅里叶变换,实现由频域谱到时域谱的变换,根据反射峰间的光程差得到吸收光程;
步骤五、选取时域谱中每个分布式气体装置8末端的信息数据段进行逆傅里叶变换,再对幅值数据作平方处理,得到气体的吸收光谱强度信息;
步骤六、根据气体的吸收光谱强度信息和吸收光程得到每一个分布式气体装置8内待测气体的浓度累计值;
步骤七、将当前分布式气体装置8的浓度累计值减去其前一个分布式气体装置8的浓度累计值,得到当前分布式气体装置8内待测气体的浓度。
优选的是,步骤四中,拍频相干信号以电流表征,电流i(ω)为:
Cm为常数、IR和IP分别为参考光和测试光的强度,ω为可调谐光源的瞬时角频率,τm为第m个反射点的测试光相对于参考光的延时时间、Tm(ω)为第m个反射点对应的分布式气体装置8内气体的光谱透射率;
对公式2进行傅里叶变换得到傅里叶变换后的电流
其中,ω0为开启可调谐光源1的中心角频率,Δω为角频率调谐范围,
Δν为频率调谐范围,τ为探测光相对于参考光的延迟时间;
第i个反射点到第j个反射点的吸收光程为:
i,j=1,2,3......M,M为反射点的总数目,n为气体的折射率。
优选的是,步骤五中,对公式3中每个分布式气体装置8末端的信息数据段进行逆傅里叶变换,得到:
再进行平方得到,
优选的是,步骤六具体为:
其中,ρ为气体的浓度累计值、M为气体的摩尔质量、l为吸收光程、S为气体的吸收线强、Φ为标准大气压下气体的吸收线型、NA为阿伏伽德罗常数,
其中,Cm为常数、IR和IP分别为参考光和测试光的强度、为气体的吸收光谱强度信息。
本发明构造出气体与光相互作用的结构,通过设置可调谐光源的调谐波段,满足并覆盖待测气体的吸收波段,通过相干检测调频连续波测试光的回波反射信号,解调出各个分布式气体装置中气体的浓度,达到气体浓度的定量检测的目的,实现分布式气体传感。
本发明的有益效果:
本发明通过改变可调谐光源的波长调谐范围,使其适应覆盖待测气体的吸收波段,即可以实现对目标气体浓度的检测,该发明对商用的半导体激光光源具有普遍适用性和兼容性;本发明可以实现长距离(km量级)上的气体浓度分布式检测,较现有气体传感装置,具有测量距离长、分布式测量、定量检测的优点。
附图说明
图1是具体实施方式一所述的一种基于调频连续波的分布式气体传感装置的结构示意图;
图2是具体实施方式二所述的一种基于调频连续波的分布式气体传感装置的结构示意图;
图3是具体实施方式二中的分布式气体装置的结构示意图;
图4是实施方式四所述的一种基于调频连续波的分布式气体传感方法的流程图。
具体实施方式
具体实施方式一:本实施方式所述的一种基于调频连续波的分布式气体传感装置,包括可调谐光源1、辅助干涉仪、主干涉仪、数据采集卡10和计算模块11;
采用的光纤为普通单模光纤,波长在1550nm附近的可调谐激光器作为光源。
可调谐光源1输出的激光分别入射至辅助干涉仪和主干涉仪;
辅助干涉仪,用于产生时钟信号,并发送给数据采集卡10的时钟端;
主干涉仪的光路中包括参考光和测试光,测试光依次进入级联的多个分布式气体装置8,分布式气体装置8内装有待测气体,测试光经多级反射及气体的吸收后反射光与参考光发生拍频相干;
数据采集卡10,用于对主干涉仪产生的拍频相干信号进行采集,并发送给计算模块11;
计算模块11,用于根据拍频相干信号计算气体浓度。
本实施方式中,可调谐光源1输出的激光经第一时延光纤3-1后由第一耦合器2-1分成两路光,一路光入射至辅助干涉仪,另一路光入射至主干涉仪;
辅助干涉仪包括第四耦合器2-4、第一法拉第旋转镜6-1、第二法拉第旋转镜6-2和第二时延光纤3-2(10m~3000m);
第一耦合器2-1输出的一路光经50:50的2×2第四耦合器2-4分成两路,一路光经第一法拉第旋转镜6-1反射后耦合进入第四耦合器2-4,另一路光经第二时延光纤3-2后入射至第二法拉第旋转镜6-2,再经第二法拉第旋转镜6-2反射后耦合进入第四耦合器2-4,第四耦合器2-4输出拍频相干信号,然后经TTL时钟发生器7产生方波信号,该方波信号即为时钟信号。
本实施方式中,主干涉仪包括第二耦合器2-2、第三耦合器2-3、偏振控制器4、光纤环形器5、多个分布式气体装置8和平衡探测器9;
第一耦合器2-1输出的另一路光经1:99的1×2的第二耦合器2-2分成两路光,1%的一路光作为参考光的偏振态通过偏振控制器4调整后入射至第三耦合器2-3的一个输入端口,99%的另一路光作为测试光由光纤环形器5的1端口进入,由光纤环形器5的2端口依次入射至级联的多个分布式气体装置8,在分布式气体装置8内经多级反射以及气体的吸收后反射光重新进入光纤环行器5的2端口并从3端口耦合进入第三耦合器2-3的另一个输入端口,参考光与测试光发生拍频相干并由第三耦合器2-3的2个输出端输出至平衡探测器9,平衡探测器9探测并转换为电压信号,数据采集卡10将模拟电压转换成数字电压值。
计算模块11根据拍频相干信号计算气体浓度,具体为:
计算模块11对接收到的拍频相干信号进行傅里叶变换,实现由频域谱到时域谱的变换,根据反射峰间的光程差得到吸收光程;
选取时域谱中每个分布式气体装置8末端的信息数据段进行逆傅里叶变换,再对幅值数据作平方处理,得到气体的吸收光谱强度信息;
根据吸收光谱强度信息和吸收光程得到每一个分布式气体装置8内待测气体的浓度累计值;
将当前分布式气体装置8的浓度累计值减去其前一个分布式气体装置8的浓度累计值,得到当前分布式气体装置8内待测气体的浓度。
相比于化学检测手段,光学方法更具有高分辨率,高灵敏度,实时性等诸多优点,利用可调谐光源1在近红外等波段上的成熟技术优势以及光学器件,将单模光纤作为气体传感的信息搭载媒介,借助光纤的波导结构特性,并构造及利用气体与光相互作用的结构(机械式级联气池或空心光子晶体光纤),并利用一定的对接方式实现级联,实现分布式长距离气体传感,同时具备实时性强、传感分辨率高等优点。
将气体充入机械式级联气池或空心光子晶体光纤,通过设置***光源的波段调谐范围使其覆盖感兴趣的样品气体吸收波段,实现光与气体的相互作用。通过解调调频连续波的回波信号,通过一定的算法即可解调出在某个气池或者某段空心光子晶体光纤中的气体浓度,利用光纤对分布在各个监测点的气池或者空心光子晶体光纤进行级联,形成一次性定量检测各个位置的气体浓度的分布式气体传感***,超越其他现有气体传感器的单点测量的局限性和能力范围。
具体实施方式二:结合图1和图2具体说明本实施方式,本实施方式是对具体实施方式一所述的一种基于调频连续波的分布式气体传感装置作进一步说明,本实施方式中,分布式气体装置8包括机械钢体气体容器8-1、第一光纤准直器8-2、第二光纤准直器8-3、进气口阀门8-4和出气口阀门8-5;
机械钢体气体容器8-1的两端分别设有第一光纤准直器8-2和第二光纤准直器8-3,靠近机械钢体气体容器8-1的两端分别设有进气口阀门8-4和出气口阀门8-5。
本实施方式中,分布式气体装置8采用机械式级联气池实现。一个准直器作为发射端,另一个准直器作为接收端,气体通过进气口阀门/出气口阀门进入/流出机械级联气池中,整个结构在常压下成密封状态。
图2中。为气体分子,为准直光束。
具体实施方式三:结合图3具体说明本实施方式,本实施方式是对具体实施方式一所述的一种基于调频连续波的分布式气体传感装置作进一步说明,本实施方式中,分布式气体装置8包括空心光子晶体光纤和2个三通气阀;
2个三通气阀中各有1个口分别用于进气和出气;
空心光子晶体光纤的两端各通过1个三通气阀的另两个口与光路中的光纤耦合。两端衔接处在常压下成密封状态。
本实施方式中,分布式气体装置8采用空心光子晶体光纤实现。
具体实施方式四:结合图4具体说明本实施方式,本实施方式所述的一种基于调频连续波的分布式气体传感方法,该方法基于具体实施方式二所述的一种基于调频连续波的分布式气体传感装置实现,该方法包括:
步骤一、开启可调谐光源1线性扫频输出,可调谐光源1输出触发信号触发数据采集卡10作数据采集准备;
步骤二、辅助干涉仪产生时钟信号,传输至数据采集卡10,使数据采集卡10对主干涉仪的拍频信号进行等频率间隔采集;
主干涉仪的测试光由光纤环行器5的2端口导入级联的多个分布式气体装置8,在分布式气体装置8内经多级反射以及气体的吸收后反射光重新进入光纤环行器5;调节偏振控制器4使得从光纤环行器5的3端口导出的光与参考光发生拍频相干,然后由平衡探测器9传输至数据采集卡10;
步骤三、数据采集卡10采集拍频相干信号;
步骤四、计算模块11对拍频相干信号进行保存,并对拍频相干信号进行快速傅里叶变换,实现由频域谱到时域谱的变换,根据反射峰间的光程差得到吸收光程;
步骤五、选取时域谱中每个分布式气体装置8末端的信息数据段进行逆傅里叶变换,再对幅值数据作平方处理,得到级联气池叠加后的气体的吸收光谱强度信息;
步骤六、根据气体的吸收光谱强度信息和吸收光程得到每一个分布式气体装置8内待测气体的浓度累计值;
步骤七、将当前分布式气体装置8的浓度累计值减去其前一个分布式气体装置8的浓度累计值,得到当前分布式气体装置8内待测气体的浓度。
通过一次性计算得到分布式气体装置中各级级联气池的叠加后的气体浓度后,为获得到当前环境下某个级联气池中目标的浓度,需要将该环境下该级联气池叠加后的浓度减去其前一个级联气池中的叠加后的浓度得到(第一个级联气池的浓度则为真实目标环境下的气体浓度)。
本实施方式基于调频连续波(FMCW)干涉测量法,交流的平衡探测器9用于探测主干涉仪上的拍频信号的干涉项,其探测器电流可以表示为:
其中ω是可调谐光源的瞬时角频率,τ是主干涉仪中探测光相对于参考光的延迟时间,C为包含探测器灵敏度的常数、特定反射点的反射率以及光学损耗系数(除气体吸收引起的光学衰减外)的常数,IR和IP分别为参考光和测试光的强度,根据BeerLambert定律T(ω,τ)=e-α(ω)l(τ)表述,T(ω,τ)为级联气池内气体的光谱透射率,其中α(ω)和l(τ)分别表示待测气体的吸收系数和吸收光程。
多个级联气池级联后形成探测臂,由于探测臂上级联气池中每个准直器的发射端和接受端上都会发生菲涅尔反射,以及连接光纤中发生的瑞利后向散射,探测器得到的探测电流可以表示成:
其中m表示沿探测臂的第m个反射点。
通过傅里叶变换获得与不同反射位置相关的T(ω):
其中角频率调谐范围Δω=2πΔν设置了积分极限。根据sinc函数,可确定在τ=τm反射位置处,/>存在最大值。同时,根据sinc函数,其时域分辨率为Δτ=1/Δν,对应的空间分辨率为Δl=c/2ngΔν,其中c为光在真空中的速度,ng为光传播的介质群折射率。由于该***的测试光从起点发射经实际反射点或反射面反射后作一次往返式传播,与参考光发生相干干涉,其相对时延所对应的光程是起点与反射点实际距离的两倍,因此吸收光程也是起点与反射点实际距离的两倍,最终的吸收光程为:
i,j=1,2,3......M,M为反射点的总数目,n为气体的折射率。
通过式4可知,基于调频连续波技术,可以实现各个反射点的位置测定,并通过两个反射点之间的距离来确定最终的吸收光程,为实现气体传感提供了有效的技术路径。
对公式3中每个分布式气体装置8末端的信息数据段进行逆傅里叶变换,得到:
则透射函数Tm(ω)可以表示为:
为推算某气体浓度ρ(g/cm3),需求得该气体的分子数密度N(molecule/cm3),已知气体的摩尔质量M(g/mol)和阿伏伽德罗常数NA(molecule/mol),气体浓度可表示为:
由HITRAN数据库获取气体的吸收线强S以及标准大气压下气体的吸收线型Φ(ν),在获知具体透射函数T(ω)以及吸收光程l(τ)的情况下,最终气体浓度简化为:
对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神或基本特征的情况下,能够以其他的具体形式实现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。
Claims (9)
1.一种基于调频连续波的分布式气体传感装置,包括可调谐光源(1)、辅助干涉仪、主干涉仪、数据采集卡(10)和计算模块(11);
可调谐光源(1)输出的激光分别入射至辅助干涉仪和主干涉仪;
辅助干涉仪,用于产生时钟信号,并发送给数据采集卡(10)的时钟端;
主干涉仪的光路中包括参考光和测试光,测试光依次进入级联的多个分布式气体装置(8),分布式气体装置(8)内装有待测气体,测试光经多级反射及气体的吸收后反射光与参考光发生拍频相干;
数据采集卡(10),用于对主干涉仪产生的拍频相干信号进行采集,并发送给计算模块(11);
计算模块(11),用于根据拍频相干信号计算气体浓度;
其特征在于,计算模块(11)根据拍频相干信号计算气体浓度,具体为:
计算模块(11)对接收到的拍频相干信号进行傅里叶变换,实现由频域谱到时域谱的变换,根据反射峰间的光程差得到吸收光程;
选取时域谱中每个分布式气体装置(8)末端的信息数据段进行逆傅里叶变换,再对幅值数据作平方处理,得到气体的吸收光谱强度信息;
根据吸收光谱强度信息和吸收光程得到每一个分布式气体装置(8)作为反射点的浓度累计值;
将当前分布式气体装置(8)作为反射点的浓度累计值减去其前一个分布式气体装置(8)作为反射点的浓度累计值,得到当前分布式气体装置(8)内待测气体的浓度。
2.根据权利要求1所述的一种基于调频连续波的分布式气体传感装置,其特征在于,可调谐光源(1)输出的激光经第一时延光纤(3-1)后由第一耦合器(2-1)分成两路光,一路光入射至辅助干涉仪,另一路光入射至主干涉仪;
辅助干涉仪包括第四耦合器(2-4)、第一法拉第旋转镜(6-1)、第二法拉第旋转镜(6-2)和第二时延光纤(3-2);
第一耦合器(2-1)输出的一路光经第四耦合器(2-4)分成两路,一路光经第一法拉第旋转镜(6-1)反射后耦合进入第四耦合器(2-4),另一路光经第二时延光纤(3-2)后入射至第二法拉第旋转镜(6-2),再经第二法拉第旋转镜(6-2)反射后耦合进入第四耦合器(2-4),第四耦合器(2-4)输出拍频相干信号,然后经TTL时钟发生器(7)产生方波信号,该方波信号即为时钟信号。
3.根据权利要求2所述的一种基于调频连续波的分布式气体传感装置,其特征在于,主干涉仪包括第二耦合器(2-2)、第三耦合器(2-3)、偏振控制器(4)、光纤环形器(5)、多个分布式气体装置(8)和平衡探测器(9);
第一耦合器(2-1)输出的另一路光经第二耦合器(2-2)分成两路光,一路光作为参考光经偏振控制器(4)后入射至第三耦合器(2-3)的一个输入端口,另一路光作为测试光由光纤环形器(5)的1端口进入,由光纤环形器(5)的2端口出射至级联的多个分布式气体装置(8),在分布式气体装置(8)内经多级反射以及气体的吸收后反射光重新进入光纤环形器(5)的2端口并从3端口耦合进入第三耦合器(2-3)的另一个输入端口,参考光与测试光发生拍频相干并由第三耦合器(2-3)的2个输出端输出至平衡探测器(9),平衡探测器(9)的输出端连接数据采集卡(10)的采集信号输入端。
4.根据权利要求1所述的一种基于调频连续波的分布式气体传感装置,其特征在于,分布式气体装置(8)包括机械钢体气体容器(8-1)、第一光纤准直器(8-2)、第二光纤准直器(8-3)、进气口阀门(8-4)和出气口阀门(8-5);
机械钢体气体容器(8-1)的两端分别设有第一光纤准直器(8-2)和第二光纤准直器(8-3),靠近机械钢体气体容器(8-1)的两端分别设有进气口阀门(8-4)和出气口阀门(8-5)。
5.根据权利要求1所述的一种基于调频连续波的分布式气体传感装置,其特征在于,分布式气体装置(8)包括空心光子晶体光纤和2个三通气阀;
2个三通气阀中各有1个口分别用于进气和出气;
空心光子晶体光纤的两端各通过1个三通气阀的另两个口与光路中的光纤耦合。
6.一种基于调频连续波的分布式气体传感方法,其特征在于,该方法包括:
步骤一、开启可调谐光源(1)线性扫频输出,可调谐光源(1)输出触发信号触发数据采集卡(10)作数据采集准备;
步骤二、辅助干涉仪产生时钟信号,传输至数据采集卡(10),使数据采集卡(10)对主干涉仪的拍频信号进行等频率间隔采集;
主干涉仪的测试光由光纤环形器(5)的2端口导入级联的多个分布式气体装置(8),在分布式气体装置(8)内经多级反射以及气体的吸收后反射光重新进入光纤环形器(5);调节偏振控制器(4)使得从光纤环形器(5)的3端口导出的光与参考光发生拍频相干,然后由平衡探测器(9)传输至数据采集卡(10);
步骤三、数据采集卡(10)采集拍频相干信号;
步骤四、计算模块(11)对拍频相干信号进行保存,并对拍频相干信号进行傅里叶变换,实现由频域谱到时域谱的变换,根据反射峰间的光程差得到吸收光程;
步骤五、选取时域谱中每个分布式气体装置(8)末端的信息数据段进行逆傅里叶变换,再对幅值数据作平方处理,得到气体的吸收光谱强度信息;
步骤六、根据气体的吸收光谱强度信息和吸收光程得到每一个分布式气体装置(8)内待测气体的浓度累计值;
步骤七、将当前分布式气体装置(8)的浓度累计值减去其前一个分布式气体装置(8)的浓度累计值,得到当前分布式气体装置(8)内待测气体的浓度;
该方法基于权利要求3所述的一种基于调频连续波的分布式气体传感装置实现。
7.根据权利要求6所述的一种基于调频连续波的分布式气体传感方法,其特征在于,步骤四中,拍频相干信号以电流表征,电流i(ω)为:
Cm为常数、IR和IP分别为参考光和测试光的强度,ω为可调谐光源的瞬时角频率,τm为第m个反射点的测试光相对于参考光的延时时间、Tm(ω)为第m个反射点对应的分布式气体装置(8)内气体的光谱透射率;
对公式(2)进行傅里叶变换得到傅里叶变换后的电流
其中,ω0为开启可调谐光源(1)的中心角频率,Δω为角频率调谐范围,Δv为频率调谐范围,τ为探测光相对于参考光的延迟时间;
第i个反射点到第j个反射点的吸收光程lij(τ)为:
i,j=1,2,3......M,M为反射点的总数目,n为气体的折射率。
8.根据权利要求7所述的一种基于调频连续波的分布式气体传感方法,其特征在于,步骤五中,对公式3中每个分布式气体装置(8)末端的信息数据段进行逆傅里叶变换,得到:
再进行平方得到,
9.根据权利要求8所述的一种基于调频连续波的分布式气体传感方法,其特征在于,步骤六具体为:
其中,ρ为气体的浓度累计值、M为气体的摩尔质量、l为吸收光程、S为气体的吸收线强、Φ为标准大气压下气体的吸收线型、NA为阿伏伽德罗常数,
其中,Cm为常数、IR和IP分别为参考光和测试光的强度、为气体的吸收光谱强度信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811101751.7A CN109085138B (zh) | 2018-09-20 | 2018-09-20 | 一种基于调频连续波的分布式气体传感装置和方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811101751.7A CN109085138B (zh) | 2018-09-20 | 2018-09-20 | 一种基于调频连续波的分布式气体传感装置和方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109085138A CN109085138A (zh) | 2018-12-25 |
CN109085138B true CN109085138B (zh) | 2024-02-13 |
Family
ID=64842048
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811101751.7A Active CN109085138B (zh) | 2018-09-20 | 2018-09-20 | 一种基于调频连续波的分布式气体传感装置和方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109085138B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109696418B (zh) * | 2019-02-27 | 2021-08-10 | 南京科远智慧科技集团股份有限公司 | 一种应用于低信噪比下的tdlas直接吸收法的外差平衡探测***及探测方法 |
CN110907384A (zh) * | 2019-12-10 | 2020-03-24 | 中国海洋大学 | 一种基于傅里叶域光学相干吸收光谱技术的气体检测***及其工作方法 |
CN113804300B (zh) * | 2021-07-20 | 2023-12-26 | 广东工业大学 | 一种基于光频域最优移频干涉的分布式偏振串音测量装置 |
CN114002185B (zh) * | 2021-11-17 | 2023-08-04 | 哈尔滨工业大学 | 基于光学调频连续波的多点色散光谱测量装置和方法 |
CN114384045B (zh) * | 2022-01-10 | 2024-01-12 | 北京航空航天大学 | 一种痕量气体浓度和路径长度实时检测***与方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104483282A (zh) * | 2014-12-31 | 2015-04-01 | 郑州光力科技股份有限公司 | 自校准气室及光纤气体传感器检测*** |
-
2018
- 2018-09-20 CN CN201811101751.7A patent/CN109085138B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104483282A (zh) * | 2014-12-31 | 2015-04-01 | 郑州光力科技股份有限公司 | 自校准气室及光纤气体传感器检测*** |
Non-Patent Citations (2)
Title |
---|
Simultaneous measurement of gas absorption spectra and optical path lengths in a multipass cell by FMCW interferometry;Xiutao Lou, et. Al.;《Optics Letters》;20180615;第43卷(第12期);第2872-2875页 * |
Xiutao Lou, et. Al..Simultaneous measurement of gas absorption spectra and optical path lengths in a multipass cell by FMCW interferometry.《Optics Letters》.2018,第43卷(第12期),第2872-2875页. * |
Also Published As
Publication number | Publication date |
---|---|
CN109085138A (zh) | 2018-12-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109085138B (zh) | 一种基于调频连续波的分布式气体传感装置和方法 | |
Dobler et al. | Atmospheric CO 2 column measurements with an airborne intensity-modulated continuous wave 1.57 μm fiber laser lidar | |
CN104280362B (zh) | 一种高温水汽激光光谱在线检测*** | |
US7884937B2 (en) | Airborne tunable mid-IR laser gas-correlation sensor | |
CN102798411B (zh) | 基于布里渊散射的分布式光纤传感测量***及测量方法 | |
CN100483107C (zh) | 分布式微结构光纤瓦斯传感***及传感方法 | |
US20210349011A1 (en) | Laser radar system apparatus for multi-wavelength measurement of atmospheric carbon dioxide concentration and vertical aerosol profile | |
CN103698298A (zh) | 基于短腔腔增强关联光谱技术测量气体浓度的装置及采用该装置测量气体浓度的方法 | |
CN105067563B (zh) | 开放空间气体平均浓度测量装置及测量方法 | |
CN103713293A (zh) | 全光纤多普勒激光雷达风场探测***及探测方法 | |
US10228490B2 (en) | Hub and spoke system for detecting and locating gas leaks | |
WO2019112459A1 (ru) | Способ дистанционного измерения концентрации газов в атмосфере | |
CN108548644A (zh) | 一种基于光纤氧气传感的联通石油储罐泄漏监测装置 | |
RU2736178C1 (ru) | Способ и устройство для автономного дистанционного определения концентрации атмосферных газовых составляющих | |
Wang et al. | Vertical wind profiling with fiber-Mach-Zehnder-interferometer-based incoherent Doppler lidar | |
Cézard et al. | Performance evaluation of a dual fringe-imaging Michelson interferometer for air parameter measurements with a 355 nm Rayleigh-Mie lidar | |
Leonard et al. | Remote sensing of ocean physical properties: a comparison of Raman and Brillouin techniques | |
CN115826088A (zh) | 中高层大气垂直风廓线的激光外差光谱测量方法及*** | |
Matson et al. | Infrared Fiber Optic Sensors for the Remote Detection of Hydrocarbons Operating in the 3.3 to 3.36 Micron Region | |
CN108593597A (zh) | 基于光纤fp腔型探头的天然气泄漏预警监控装置及方法 | |
RU2598694C2 (ru) | Устройство и способ измерения концентрации газообразных веществ | |
CN115343233B (zh) | 一种开放路径上痕量气体浓度的实时测量方法与装置 | |
CN208283270U (zh) | 基于光纤fp腔型探头的天然气泄漏预警监控装置 | |
Vasil’chenko et al. | Spectral system for measuring gaseous atmospheric components with a fiber-optic tracking system, and certain analysis results of atmospheric spectra | |
Wang et al. | Design and experimental verification of a novel Mie Doppler wind lidar based on all-fiber Mach-Zehnder frequency discriminator |
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 | ||
CB03 | Change of inventor or designer information |
Inventor after: Lou Xiutao Inventor after: Dong Yongkang Inventor after: Chen Chen Inventor before: Dong Yongkang Inventor before: Lou Xiutao Inventor before: Chen Chen |
|
CB03 | Change of inventor or designer information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |