CN115639170A - 一种基于波分复用的抗光学偏折吸收谱提取方法与系统 - Google Patents
一种基于波分复用的抗光学偏折吸收谱提取方法与系统 Download PDFInfo
- Publication number
- CN115639170A CN115639170A CN202211356722.1A CN202211356722A CN115639170A CN 115639170 A CN115639170 A CN 115639170A CN 202211356722 A CN202211356722 A CN 202211356722A CN 115639170 A CN115639170 A CN 115639170A
- Authority
- CN
- China
- Prior art keywords
- laser
- absorption
- light
- optical fiber
- generation module
- 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.)
- Granted
Links
- 238000000862 absorption spectrum Methods 0.000 title claims abstract description 47
- 238000000605 extraction Methods 0.000 title claims abstract description 14
- 238000010521 absorption reaction Methods 0.000 claims abstract description 75
- 230000000694 effects Effects 0.000 claims abstract description 41
- 239000013307 optical fiber Substances 0.000 claims abstract description 40
- 230000008859 change Effects 0.000 claims abstract description 19
- 238000000034 method Methods 0.000 claims description 44
- 230000003287 optical effect Effects 0.000 claims description 17
- 238000001228 spectrum Methods 0.000 claims description 7
- 230000008878 coupling Effects 0.000 claims description 6
- 238000010168 coupling process Methods 0.000 claims description 6
- 238000005859 coupling reaction Methods 0.000 claims description 6
- 238000002835 absorbance Methods 0.000 claims description 3
- 238000013480 data collection Methods 0.000 claims description 2
- 238000000926 separation method Methods 0.000 claims 1
- 238000012360 testing method Methods 0.000 abstract description 27
- 238000002485 combustion reaction Methods 0.000 abstract description 19
- 238000003745 diagnosis Methods 0.000 abstract 1
- 239000007789 gas Substances 0.000 description 33
- 238000005259 measurement Methods 0.000 description 21
- 238000000041 tunable diode laser absorption spectroscopy Methods 0.000 description 19
- 238000001514 detection method Methods 0.000 description 12
- 238000005516 engineering process Methods 0.000 description 10
- 239000000835 fiber Substances 0.000 description 9
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 6
- 238000003384 imaging method Methods 0.000 description 6
- 238000001285 laser absorption spectroscopy Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 238000004847 absorption spectroscopy Methods 0.000 description 5
- HSFWRNGVRCDJHI-UHFFFAOYSA-N alpha-acetylene Natural products C#C HSFWRNGVRCDJHI-UHFFFAOYSA-N 0.000 description 5
- 125000002534 ethynyl group Chemical group [H]C#C* 0.000 description 5
- 230000035939 shock Effects 0.000 description 5
- 230000003595 spectral effect Effects 0.000 description 5
- 238000004611 spectroscopical analysis Methods 0.000 description 5
- 238000009529 body temperature measurement Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000011065 in-situ storage Methods 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 230000003313 weakening effect Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 229910002092 carbon dioxide Inorganic materials 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 239000002245 particle Substances 0.000 description 3
- 239000012071 phase Substances 0.000 description 3
- 238000011084 recovery Methods 0.000 description 3
- 230000002269 spontaneous effect Effects 0.000 description 3
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 description 2
- MYMOFIZGZYHOMD-UHFFFAOYSA-N Dioxygen Chemical compound O=O MYMOFIZGZYHOMD-UHFFFAOYSA-N 0.000 description 2
- 239000000654 additive Substances 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 2
- 229910002091 carbon monoxide Inorganic materials 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 229910001882 dioxygen Inorganic materials 0.000 description 2
- 238000012625 in-situ measurement Methods 0.000 description 2
- 230000005764 inhibitory process Effects 0.000 description 2
- 238000012886 linear function Methods 0.000 description 2
- 238000004476 mid-IR spectroscopy Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000002829 reductive effect Effects 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 244000145845 chattering Species 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 239000007791 liquid phase Substances 0.000 description 1
- 238000004020 luminiscence type Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 230000005180 public health Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
本发明提供了一种基于波分复用的抗光学偏折吸收谱提取方法与系统。系统包括吸收激光发生模块、参考激光发生模块、光纤耦合器、分束器、干涉仪、光纤准直镜对、波分复用器、光电探测器、数据采集模块与计算机。吸收激光发生模块发出的激光分成两束；一束由干涉仪记录波长，另一束与参考激光发生模块发出的无吸收谱段激光合束后通过光纤准直镜发出；合束激光穿过待测气体后由光纤准直镜收集到波分复用器中进行分束，再分别由两个光电探测器将光强信号转换成电压信号；利用参考激光的强度变化补偿光线偏折效应导致的吸收激光的强度变化，进而实现吸收光谱的准确提取。该发明能够克服光线偏折效应对激光吸收光谱的干扰，在燃烧诊断等动态吸收光谱测试领域有广阔应用前景。
Description
技术领域
本发明涉及激光气体吸收光谱测量系统与方法，尤其涉及一种基于波分复用的抗光学偏折吸收光谱提取方法与系统，属于激光光谱及气体参数测量技术领域。
背景技术
可调谐二极管激光吸收光谱技术(Tunable diode laser absorptionspectroscopy，TDLAS) 利用窄带激光扫描气体的吸收谱，通过获取气体分子对激光的吸收率，进而获取待测气体的浓度、温度等信息。从上世纪七十年代起，该技术被应用于燃烧场的参数测量，具有非侵入性，测量速度快，灵敏度高的特点。
目前，TDLAS的主要方法有直接吸收法(Direct absorption Spectroscopy，DAS)与波长调制法(Wavelength Modulation Spectroscopy，WMS)，两种方法的目的在于准确获取气体吸收线型函数。DAS是将透射光强和入射光强的比值直接拟合得到气体的吸收率线型函数，拟合得到的吸收率线型函数包括了待测气体的温度、浓度和压力等信息。但直接吸收法易受颗粒物浓度、激光强度波动和高压下谱线重叠等因素的影响，导致测量误差的增大，直接吸收法更适合于强吸收条件下的测量。Adam J Reid等在1981年发表在《应用物理学报》(Journal of Applied Physics)第26卷第203页的论文《可调谐二极管激光器的二次谐波检测——实验与理论的比较》(Second-harmonic detection with tunable diodelasers—comparison of experiment and theory)提出了波长调制方法，针对窄带激光扫描过程进行了高频调制，利用高频信号的高次谐波进行气体吸收率的计算与气体参数的提取，该方法有效地降低了测量系统中的背景干扰，提高了TDLAS的测量精度。但是单纯利用高次谐波信号进行气体参数的计算需要对已知气体进行标定，增加了系统的复杂度。WeiWei等在2019年发表在《应用物理B》(Applied Physics B)第125卷的第9号论文《在高压冲击管中使用波长调制光谱方法抑制非吸收干扰》(Demonstration of non-absorbinginterference rejection using wavelength modulation spectroscopy in high-pressure shock tubes)直接比较了波长调制法(Wavelength Modulation Spectroscopy，WMS)的测量值和直接吸收法的测量值，对WMS检测系统在3.5atm、8.5atm和18atm 的反射冲击压力下进行了评估，结果表明，与并行DAS测量相比，WMS检测系统的信噪比有所提高。WMS方法一般是通过在锯齿波上叠加了正弦信号的电流控制激光器，并对采集到的信号进行解调处理，获取激光吸收线型的部分参数。由于进行了调制解调的过程，WMS测得的信号的信噪比要比DAS方法高，但得到的谐波分量无法直接用于恢复吸收谱。
尽管可调谐二极管激光吸收光谱技术非常适合燃烧监测，但该方法作为一种利用视线场激光的吸收效应进行的测量技术，当燃烧场中的燃烧产物浓度过大，或者是燃烧器内的压力环境发生改变之后，光线的穿过的路线上的折射率发生变化后，激光光线在不同折射率的环境上不再是沿直线传播。由于燃烧场的高速动态环境，这种由于浓度和压力的改变引起的折射率场的改变，引起了光线传播路径的变化，使得激光照射到光电探测器上的光照强度发生了变化，这会引起非吸收性质的光强变弱。基于激光吸收光谱的测量技术的本身也是根据气体的对激光的吸收度判断气体的浓度、温度与压力等状态，当这种非吸收效应引起的光强变化耦合到了气体吸收效应引起的光强变化之中，会导致无法对气体的吸收谱进行准确的提取。H.Eric Schlosser等人在2002年发表在《国际燃烧会议论文集》(Proceedings of the Combustion Institute)上第29卷第353-360页的论文《用可调谐二极管激光吸收光谱法在全尺寸灭火试验中就地测定分子氧浓度》(In situ determinationof molecular oxygen concentrations in full-scale fire-suppression tests usingtunable diode laser absorption spectroscopy)就提出了基于激光吸收光谱在燃烧应用中的模型，该模型考虑了在火焰监测的过程中需要考虑的火焰自发辐射，以及火焰燃烧产物压力与浓度的变化给光线带来的偏折，在探测器上导致光强变化的影响。该模型中，火焰的辐射噪声被视为加法噪声，光纤偏折与颗粒散射等因素导致的噪声被视为乘法噪声。Lin Ma等在2007年发表在《应用光学》(Applied Optics)第21期第47卷第3751页的论文《温度和浓度同步成像的高光谱层析成像数值研究》(Numerical investigation ofhyperspectral tomography for simultaneous temperature and concentrationimaging)在对基于 TDLAS技术的燃烧场测量的时候添加了0.01％的噪声，用于表示基础加性随机噪声，用 5.12％噪声等级表示了光线偏折等效应带来的乘法噪声，说明了在TDLAS的测试应用中，来源于光线偏折效应的噪声是主要的噪声来源。
针对光线偏折效应的抑制与降低，大量研究人员进行了丰富的研究。直接把光线带来的噪声项作为一个待求解的参数，代入激光吸收光谱模型，从而利用采集到的数据对包含偏折效应的吸收谱模型进行统一拟合。Steven Wagner等在2009年发表在《Proceedings of the Combustion Institute》第32卷第839-846页上的论文《基于TDLAS的二维层流扩散火焰中乙炔绝对浓度的原位测量》(TDLAS-based in situ measurementof absolute acetylene concentrations in laminar 2D diffusion flames)在利用直接吸收法对乙炔浓度测量的过程中，将光线偏折与颗粒散射造成的光线衰减定义成了吸收率上的乘法噪声，并将其当作随机噪声代入了吸收谱线型的拟合中，拟合结果误差大。S.Wagner.等人在2012年发表在《应用物理B》(Applic Physic B)上第107期第585-589页的论文《非预混层流逆流火焰中乙炔绝对浓度剖面的原位TDLAS测量》(In situ TDLASmeasurement of absolute acetylene concentration profiles in a non-premixedlaminar counter-flow flame)讨论了基于高斯光束的激光在发生偏折时候的误差，通过增加光束直径，降低了光线偏折带来的影响，并通过对偏差进行建模，利用非线性拟合消除影响。T.Kamimoto等在2015 年发表在《流量测量和仪器仪表》(Flow Measurement andInstrumentation)第46卷第51- 57页的论文《温度和浓度同步成像的高光谱层析成像数值研究》(High-temperature field application of two-dimensional temperaturemeasurement technology using CT tunable diode laser absorption spectroscopy)中，在利用TDLAS技术对发动机的点火过程的水分温度与浓度变化的二维成像监测中，使用了多项式拟合技术对光线偏折效应进行了补偿。这种方式能够在一定程度上缓解光线偏折效应对吸收谱求解的影响，但是无法解决吸收效应与偏折效应同时造成光线变弱的吸收谱恢复，并且这种方法依赖多次扫描的数据量，使得测试帧率与速度下降。
从光学角度出发，采用聚光镜的方案降低光线偏折效应的方案也被广泛研究。采用聚光镜将不同方向和不同位置的激光收束到一个稳定的光斑，从而避免光线偏折效应带来的非吸收误差被广泛研究与使用。P.weibring等人在2006年发表在《应用物理B》(Applic Physic B)上第85期第207-218页的论文《超高精密中红外光谱仪II：系统描述和光谱性能》(Ultra-high-precision mid-IR spectrometer II:system description andspectroscopic performance)在使用超快的中红外的激光器进行测试的时候，由于光线偏折的影响，导致了信号的失准，采用了一个类似照相机使用的消色差透镜实现了光学汇聚，提高了光学系统对光线抖动带来的影响的抵抗能力。H.Yang.等人在2011年发表在《应用物理B》(Applic Physic B)上第104期第21-27页的论文《可调谐二极管激光吸收传感器同步水膜厚度、液相和气相的测量温度》(Tunable diode laser absorption sensor for thesimultaneous measurement of water film thickness,liquid-and vapor-phasetemperature)在进行水膜厚度的测试过程中，使用准直镜进行光线对准接收时，由于薄膜厚度发生变化导致了光线的偏转，其采用了在准直镜前加上一个小型的收光透镜，从而降低光纤偏折抖动的对光强收集的影响。X.Lin等人在2012年发表在《应用物理B》(ApplicPhysic B) 上第110期第401-409页的论文《通过耦合OES和TDLAS用于激波管波中Martian混合物燃烧后的CO浓度和温度的测量》(CO concentration and temperaturemeasurements in a shock tube for Martian mixtures by coupling OES and TDLAS)中将火焰自发辐射的光谱与结合可调谐激光吸收光谱的测试手段，对激波管中的高温高压的一氧化碳的温度与浓度进行测量，该系统中涉及到两部分的光学对准，一个是将激光器发出的用于激光吸收光谱的激光进行收集测试，另一个是将火焰的自发光进行收集，两个方式都使用了透镜的方案进行收集，其中对火焰自发光的收集中，能够很好的抑制住在复杂环境中的噪声抖动，而在进行激光吸收光谱的测试中所需的定向激光中，透镜收集到的具有方向性的激光信号受到的干扰严重，通过一系列的降噪拟合算法，使得吸收谱的恢复误差在1％左右。这种使用透镜的方案，虽然能够在一定程度上降低沿着直线传播的激光器的方向偏转，但是由于测试环境的复杂性，复杂的光学透镜组合难以有效的布置，并且进一步提升精度的空间有限。
除了增加透镜，扩大测试激光的光斑直径，也是研究人员在克服光线偏折效应是一个常用技术。Chang Liu等人在2015年发表在《光学快讯》(Optic Express)上的论文《用于快速成像温度和气体浓度的扇形波束TDLAS层析传感器的研制》(Development of afan-beam TDLAS-based tomographic sensor for rapid imaging of temperature andgas concentration)，介绍了一种利用扇形束的激光解决横向平面内的光线偏折效应的方法，但是在纵向上依然没有得到有效解决，其在浓度与温度的结果中，使用了100组的数据进行平均，用以克服光线偏折效应带来的偏差影响。Peter Fendt等人在2021年发表在《光学快讯》(Optic Express)上第25期第29卷的论文《超连续激光在快速压缩膨胀机中的Herriott气体池的增强SMF耦合多标量燃烧诊断》(Herriott cell enhanced SMF-coupledmulti-scalar combustion diagnostics in a rapid compression expansion machineby supercontinuum laser)分析了光线偏折效应主要是由燃烧场内不均匀的密度分布与折射率分布造成的，并且提到了利用单模光纤来耦合空间光，消除光线抖动效应，但是这种方式会导致耦合到光纤中的信号偏小，引起了信号的噪声上升，测试的不确定度变大。这种使用大光束，或者使用光纤耦合方式的技术都需要牺牲激光器的能量集中性，使得探测器能够响应到的光强信号下降，信号中随机噪声上升，在实际的操作中，为了降低随机噪声，需要采取更多的测试周期以提高精度。
近些年，结合波长调制技术，也发展了一些通过调制激光信号的抗偏折方法。Sebastian Burkle等人在2018年发表在《湍流燃烧》(Flow Turbulence Combust)上第101卷第139-159页的论文《使用TDLAS测量发动机缸内温度测量的IC引擎》(In-CylinderTemperature Measurements in a Motored IC Engine using TDLAS)在利用激光吸收光谱技术对发动机缸内的燃烧过程的监测中发现在燃烧室出口处的湍流化严重，使用激光光线对吸收进行测量的时候，受到的严重的气流干扰导致光线严重偏折影响了吸收谱的恢复，其在研究中发现这种抖动现象与发动机的转速有关，因此在其展望中提到了可以通过提高扫描速度来避免光线偏折效应对吸收谱的测量产生影响。Guangzhen Gao等人在2013年发表在《光学与光谱学》(Optics and Spectroscopy)上第114卷第340-346页的论文《利用1570nm附近的可调谐二极管激光吸收光谱技术在高温下同时检测CO和CO2》(Simultaneous Detection of CO and CO2 at Elevated Temperatures Using TunableDiode Laser Absorption Spectroscopy Near 1570nm)在利用激光吸收光谱技术对高温环境下的一氧化碳与二氧化碳的浓度进行测试的时候，采用了波长调制吸收谱的测试技术，并且利用了从波长调制吸收光谱的吸收信号中提取的一次谐波信号来归一化提取的二次谐波信号，从而实现了对偏折效应的免疫。但是在其实验的过程中，是通过一个静态的温度气体池实现的温度测量，在出口处也使用了一个光学透镜进行恢复。采用波长调制方法虽然能够通过调制解调的计算，将对吸收谱的求解带宽限制在窄带范围内，但是当这种光线偏折与抖动带来的乘性噪声所在频率范围与吸收谱的频率范围相近的时候，这种窄带滤波的方法无法将噪声从吸收谱的提取过程中分离。
采用一束无吸收的参考光，与进行激光吸收光谱测试的激光共光路，进而通过分光技术将两束光的抖动分开，从而补偿光纤偏折带来的光强变化，在近年来的文献中也有体现。中国专利CN 107505063 A《一种基于高频正弦校准光的激光光线偏折校正装置及方法》(专利号：201710568850.5)介绍了一种基于高频正弦校准光的激光光线偏折校正装置及方法，包括可见光激光光源、红外可调谐激光光源、光纤耦合器、准直器、光电探测器等，控制可见光激光光源与红外可调谐激光光源，使其分别以高频正弦波形式和扫描锯齿波形式输出可见光信号和红外激光信号，其中，以可见光激光器发出的可见光波段激光作为校准光，而红外可调谐激光器发出的红外波段激光作为探测光；校准光与探测光通过光纤耦合器耦合并经过准直器后入射到被测温度区域，接收端使用一个宽波段的光电探测器同时检测透射校准光和探测光信号；基于上述过程，利用校准光与探测光信号的光线偏折及激光吸收的特征来消除光线偏折的影响。这种方案需要高频正弦波的频率至少为锯齿波扫描频率的100倍以上，以保证对光线偏折效应的细节刻画，并且其在解调的过程中，需要采用一段时间内的点数进行测量，此外，该方法使用可见光作为参考光，需要同时能够探测到可见光与红外光的探测器，大带宽的探测器在使用中会受到环境中的可见光影响，给测试带了很大的噪声。中国专利CN 107560738 A《一种基于二向色分光的激光光线偏折校正装置及方法》(专利号：201710567248.X)介绍了一种利用二向色镜对可见光光源与红外激光光源的分光方案，红外激光经过光纤分束器分成两束，一束作为参考，另一束与近红外激光通过光纤耦合器后经准直器出射，出射光穿过被测火焰后照射到二向色分光片，红外激光探测器放置于反射光线路径上，近红外激光探测器放置于透射光线路径上，利用三个激光信号之间的吸收与偏折关系消除激光偏折的影响。该方法虽然保证了在到二向分光镜之前的激光的传播路径一致，但是二向色镜之后的光线路径无法控制。中国专利CN 107560754 A《一种基于调制锯齿波信号的激光光线偏折校正装置及方法》(专利号：201710564892.1)介绍了一种利用二向色镜对可见光光源与红外激光光源进行分光的方案，控制可见光激光光源和红外可调谐激光光源以相同幅值及相位的、经过高频正弦波调制的锯齿波扫描信号形式，发出可见光激光和红外波段激光；可见光激光作为校准光，红外波段激光作为探测光；校准光与探测光通过光纤耦合器耦合、准直器准直后入射到被测温度区域，接收端利用二向色分光片将校准光与探测光分离并分别入射到可见光和近红外光电探测器上进行检测。该技术方案够在一定程度上解决光线偏折的问题，但仍需要利用一段时间内的调制信号对其进行求解。
针对以上激光吸收光谱测量技术存在的问题，本发明提出了一种基于波分复用的抗光学偏折吸收光谱提取方法与系统。通过吸收激光发生模块与参考激光发生模块两个模块分别发出一束针对待测气体有吸收的窄带扫频激光与光强恒定的无吸收的参考激光信号，两束激光耦合到一根光纤中，经过光纤准直镜发出，在测试过程中完全共光路，再利用光纤准直镜将两束激光耦合到波分复用器中，波分复用器将两束激光分开后直接利用光纤法兰紧密耦合到光电探测器上；本发明充分利用光线在光纤中传播的稳定性，避免光线偏折效应引起的非吸收效应引起的光强变化对吸收谱求解的影响，结构简单，成本低，为激光吸收光谱技术在燃烧场测试中提供了抗偏折测试的新手段。
发明内容
(一)要解决的技术问题
本发明的目的在于提供一种基于波分复用的抗光学偏折的吸收光谱提取方法与系统，要解决可调谐激光吸收光谱测试技术在实际的测试环境中受光线偏折效应从而无法正确提取气体吸收光谱的问题。
(二)技术方案
本发明，一种基于波分复用的抗光学偏折吸收谱提取方法与系统的主要实现过程包括以下步骤：
步骤一：由102吸收激光发生模块产生锯齿波调制电流驱动该模块中的可调谐激光器发出在待测气体吸收中心λ1附近的谱段扫描的激光，由101参考激光发生模块产生直流电流驱动该模块中的可调谐激光器发出波长为λ2的恒定强度的激光；其中102吸收激光发生模块发出的测试激光器经过103光纤分束器被分为两束，一束经过104干涉仪形成干涉信号后根据干涉信号记录波长，另一束和101参考激光发生模块发出的直流恒定光强的激光同时输入105光纤耦合器，形成混合的激光信号，可以表示为：
I0(t)＝IS(t)+IR(t) (1)
其中，IS为吸收激光发生模块发出的激光强度，IR为参考激光发生模块发出的激光强度，驱动吸收激光发生模块发光的电流与时间t的关系可以表示如下：
iS(t)＝a0+a1t (2)
其中，a0、a1分别表示调制电流的偏置项与一阶变化项；以大小为iR的电流驱动参考激光发生模块发光。
步骤二：105光纤耦合器将吸收激光与参考激光的混合激光传入106光纤准直镜，激光从106光纤准直镜发出穿过待测气体，由107光纤准直镜接收，将激光耦合到108波分复用器中，其中，经过光学吸收后与光线偏折后以及准直镜耦合到光纤中产生损耗后的光强信号It(t)可以表示为：
It(t)＝(IS(t)·τ(t)+IR(t))·P(t) (3)
其中，τ(t)表示气体吸收率，P(t)表示由光线偏折与耦合损耗等因素导致的光强变化。
步骤三：108波分复用器将λ1附近的扫谱吸收激光与λ2的参考激光两束混合的激光分开后，分别通过光纤法兰紧密耦合到109光电探测器与110光电探测器上，两个光电探测器将光强信号转换成电压信号，由111数据采集模块采集，可分别表示为A1和A2；111 数据模块将两个光电探测器探测到的电压信号采集后上传到112计算机。
步骤四：112计算机根据A1和A2恢复的抗光线偏折效应的光强信号为：
A(t)＝A1(t)/A2(t) (4)
进一步针对A拟合基线Aref，计算吸收度τ为：
将τ与从104干涉仪中记录并由111数据采集模解算出来的波长变化信号v按照时间对应，可恢复待测气体的吸收谱α。
(三)有益效果
本发明的有益效果在于提出一种基于波分复用器的抗光学偏折的吸收光谱提取方法与系统。主要同时发出针对待测气体有吸收波段的扫谱吸收激光与针对待测气体无吸收的直流激光；利用光纤准直镜将穿过待测气体的空间测试激光耦合到光纤波分复用器中，在空间光的部分保证了两束激光完全同光路，进而使得参考激光能够有效地刻画光线偏折效应带的光强变弱，以用于对吸收激光的非吸收效应的光强变弱进行补偿。该发明结构简单，方案切实可行，方法准确可靠，为针对燃烧场的抗光线偏折的吸收谱获取提供可靠的新方案
附图说明
附图1：一种基于波分复用的抗光学偏折吸收光谱提取系统图
附图2：吸收激光与参考激光的混合信号图
附图3：耦合光线偏折效应与吸收效应的混合激光的光强变化
附图4：波分复用器分光后光电探测器探测到的吸收激光的光强信号
附图5：波分复用器分光后光电探测器探测到的参考激光的光强信号
附图6：补偿后的吸收谱
具体实施方式
参考附图1为一种基于波分复用的抗光学偏折吸收光谱提取系统图，参考附图2为吸收激光与参考激光的混合信号图。结合一个实例的操作步骤如下：
步骤一：由102吸收激光发生模块产生锯齿波调制电流驱动该模块中的可调谐激光器发出水分子在1342nm处吸收谱的扫谱激光，由101参考激光发生模块产生直流电流驱动该模块中的可调谐激光器发出1400nm的恒定强度的激光；其中102吸收激光发生模块发出的测试激光器经过103光纤分束器被分为两束，一束经过104干涉仪形成干涉信号后根据干涉信号记录波长，另一束和101参考激光发生模块发出的直流恒定光强的激光同时输入105光纤耦合器，形成混合的激光信号，驱动吸收激光发生模块发光的电流一个扫描周期T＝1ms内与时间t的关系可以表示如下：
iS(t)＝a0+a1t (1)
其中，a0＝5mA、a1＝50mA/ms分别表示为调制电流的偏置项与一阶变化项；以大小为iR＝5mA的电流驱动参考激光发生模块发光；形成混合的激光信号如附图2所示。
步骤二：105光纤耦合器将吸收激光与参考激光的混合激光传入106光纤准直镜，激光从106光纤准直镜发出穿过待测气体，由107光纤准直镜接收，将激光耦合到108波分复用器中，其中，经过光学吸收后与光线偏折后以及准直镜耦合到光纤中产生损耗后的光强信号It(t)可以表示为：
It(t)＝(IS(t)τ(t)+IR(t))·P(t) (2)
其中，τ(t)表示气体吸收率，P(t)表示光线偏折与耦合损耗等因素导致的光强变化，如附图3所示。
步骤三：108波分复用器将λ1附近的扫谱吸收激光与λ2的参考激光两束混合的激光分开后，分别通过光纤法兰紧密耦合到109光电探测器与110光电探测器上，两个光电探测器将光强信号转换成电压信号，由111数据采集模块采集，可分别表示为A1，如附图4所示，以及A2，如附图5所示；111数据模块将从两个光电探测器探测到的电压信号采集后上传到112计算机。
步骤四：112计算机选取无吸收位置数据点，拟合基线Aref，计算吸收度τ为：
将τ与从干涉仪中记录并解算出来的波长变化信号ν按照时间对应，可恢复待测气体的吸收谱α为，如附图6所示。
以上对本发明及其实施方式的描述，并不局限于此，附图中所示仅是本发明的实施方式之一。在不脱离本发明创造宗旨的情况下，不经创造地设计出与该技术方案类似的结构或实施例，均属本发明保护范围。
Claims (3)
1.一种基于波分复用的抗光学偏折吸收谱提取方法与系统，其系统组成包括系统包括吸收激光发生模块、参考激光发生模块、光纤耦合器、分束器、干涉仪、光纤准直镜对、波分复用器、光电探测器、数据采集模块与计算机；吸收激光发生模块按照锯齿波调制的模式发出包括待测气体吸收谱段的激光信号，由吸收激光发生模块产生的锯齿波电流信号驱动激光器发出的吸收激光经过分束器分光后，一束吸收激光信号由干涉仪记录波长，另一束吸收激光信号与由参考激光发生模块产生直流信号驱动的激光器发出的激光经过光纤耦合器合束之后经准直镜发出，穿过待测气体后由另一个光纤准直镜收集，由波分复用器将两束激光分开后分别由两个光电探测器同时探测，再由数据采集模块采集后，上传到计算机，利用无吸收的参考激光的强度变化补偿光线偏折效应带给吸收激光的强度变化，实现准确的吸收光谱提取。
2.按照权利要求1所述的一种基于波分复用的抗光学偏折吸收谱提取方法与系统，其特征在于，吸收激光发生模块采用可调谐近红外激光器，使用锯齿波的电流调制可调谐激光器发出在待测气体吸收谱附近的扫描激光信号；使用的参考激光发生模块采用与吸收激光发生模块所选用的近红外激光器中心波段距离50nm以上的激光信号，利用直流电流驱动可调谐激光器发出波长和强度稳定的红外激光；两束激光经过可变增益的光纤耦合器合束后，由准直镜发出；其中，吸收激光发生模块采用的锯齿波电流在一个扫描周期T内随时间t的变化可以表示为：
iS(t)＝a0+a1t (1)
其中，a0、a1分别表示为调制电流的偏置项与一阶变化项；在此调制模式下，吸收激光发生模块控制的可调谐激光器发出中心波长为待测气体的吸收峰中心，记为λ1；参考激光发生模块控制的可调谐激光器的中心发光波长控制在λ2，保持λ1与λ2之间相差50nm以上；其中，参考激光器发生模块采用的直流电流为iR；吸收激光发生模块产生的吸收激光经过分束器分成两束光；一束激光信号由干涉仪记录波长，另一束激光信号与由参考激光发生模块产生的参考激光经过可变增益的光纤耦合器合束之后经准直镜发出，由准直镜发出的激光强度可以表示为：
I0(t)＝IS(t)+IR(t) (2)
其中，IS为吸收激光发生模块发出的激光强度，IR为参考激光发生模块发出的激光强度。
3.按照权利要求1所述的一种基于波分复用的抗光学偏折吸收谱提取方法与系统，其特征在于，吸收激光信号与参考激光发生模块产生直流信号调制的激光经过光纤耦合器合束之后由准直镜发出，穿过待测气体后由另一个光纤准直镜收集，由波分复用器将两束激光分开后分别由两个光电探测器同时探测，再由数据采集模块采集后，上传到计算机，利用参考激光的强度变化补偿光线偏折效应带给吸收激光的强度变化，实现准确的吸收光谱提取；其中，吸收激光信号与参考激光发生模块产生直流信号调制的参考激光经过光纤耦合器合束之后由准直镜发出，穿过待测气体后由另一个光纤准直镜收集的光强信号可以表示为：
It(t)＝(IS(t)·τ(t)+IR(t))·P(t) (3)
其中，τ(t)表示气体吸收率，P(t)表示由光线抖动以及准直镜到波分复用器的耦合损耗导致的光强变化；由波分复用器将在λ1附近的扫谱吸收激光与波长为λ2的参考激光的合束激光分开后分别由两个光电探测器同时探测到的信号可分别表示为A1和A2；数据采集模块将光电探测器探测到的电压信号采集后上传到计算机，计算机根据A1和A2计算抗光线偏折效应的光强信号为：
A(t)＝A1(t)/A2(t) (4)
进一步针对A拟合基线Aref，计算吸收度τ为：
将τ与从干涉仪中记录并解算出来的波长变化信号ν按照时间对应，可恢复待测气体的吸收谱α。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211356722.1A CN115639170B (zh) | 2022-11-01 | 2022-11-01 | 一种基于波分复用的抗光学偏折吸收谱提取方法与系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211356722.1A CN115639170B (zh) | 2022-11-01 | 2022-11-01 | 一种基于波分复用的抗光学偏折吸收谱提取方法与系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115639170A true CN115639170A (zh) | 2023-01-24 |
| CN115639170B CN115639170B (zh) | 2024-04-12 |
Family
ID=84946759
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211356722.1A Active CN115639170B (zh) | 2022-11-01 | 2022-11-01 | 一种基于波分复用的抗光学偏折吸收谱提取方法与系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115639170B (zh) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2003025552A1 (de) * | 2001-09-18 | 2003-03-27 | Austria Wirtschaftsservice Gesellschaft M.B.H. | Verfahren zur laser-absorptionsspektroskopie und einrichtungen zur durchführung dieses verfahrens |
| WO2013023426A1 (zh) * | 2011-08-18 | 2013-02-21 | 天津大学 | 多波段混叠式内腔气体传感系统及传感方法 |
| DE102014226487A1 (de) * | 2014-12-18 | 2016-06-23 | Bruker Optik Gmbh | FTIR-Spektrometer mit Stabilisierung des Referenzlasers über eine natürliche Absorptionslinie |
| CN107560754A (zh) * | 2017-07-12 | 2018-01-09 | 北京航空航天大学 | 一种基于调制锯齿波信号的激光光线偏折校正装置及方法 |
| CN113252204A (zh) * | 2021-04-01 | 2021-08-13 | 北京航空航天大学 | 一种基于幅度调制的激光吸收光谱温度测量方法 |
| CN114544547A (zh) * | 2022-02-23 | 2022-05-27 | 北京航空航天大学 | 一种用于波长调制光谱信号在线递推解调的方法与系统 |
| CN114548172A (zh) * | 2022-02-23 | 2022-05-27 | 北京航空航天大学 | 一种波长调制激光气体吸收光谱的参数提取方法与系统 |
-
2022
- 2022-11-01 CN CN202211356722.1A patent/CN115639170B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2003025552A1 (de) * | 2001-09-18 | 2003-03-27 | Austria Wirtschaftsservice Gesellschaft M.B.H. | Verfahren zur laser-absorptionsspektroskopie und einrichtungen zur durchführung dieses verfahrens |
| WO2013023426A1 (zh) * | 2011-08-18 | 2013-02-21 | 天津大学 | 多波段混叠式内腔气体传感系统及传感方法 |
| DE102014226487A1 (de) * | 2014-12-18 | 2016-06-23 | Bruker Optik Gmbh | FTIR-Spektrometer mit Stabilisierung des Referenzlasers über eine natürliche Absorptionslinie |
| CN107560754A (zh) * | 2017-07-12 | 2018-01-09 | 北京航空航天大学 | 一种基于调制锯齿波信号的激光光线偏折校正装置及方法 |
| CN113252204A (zh) * | 2021-04-01 | 2021-08-13 | 北京航空航天大学 | 一种基于幅度调制的激光吸收光谱温度测量方法 |
| CN114544547A (zh) * | 2022-02-23 | 2022-05-27 | 北京航空航天大学 | 一种用于波长调制光谱信号在线递推解调的方法与系统 |
| CN114548172A (zh) * | 2022-02-23 | 2022-05-27 | 北京航空航天大学 | 一种波长调制激光气体吸收光谱的参数提取方法与系统 |
Non-Patent Citations (2)
| Title |
|---|
| 李宁;吕晓静;翁春生;: "基于光强与吸收率非线性同步拟合的吸收光谱测量方法", 物理学报, no. 05, 8 February 2018 (2018-02-08), pages 243 - 250 * |
| 阚瑞峰;夏晖晖;许振宇;姚路;阮俊;范雪丽;: "激光吸收光谱流场诊断技术应用研究与进展", 中国激光, no. 09, 27 June 2018 (2018-06-27), pages 67 - 82 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115639170B (zh) | 2024-04-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105277503B (zh) | 基于两种量子级联激光光谱的多组分气体同时检测装置及方法 | |
| CN104280362B (zh) | 一种高温水汽激光光谱在线检测系统 | |
| CA3025935C (en) | Photothermal interferometry apparatus and method | |
| CN104596996B (zh) | 基于空芯光纤光热效应的气体检测方法和系统 | |
| CN101308090B (zh) | 一种火场多参量激光波长调制光谱检测方法和装置 | |
| US10302563B2 (en) | Apparatus and method of gas analysis using laser light | |
| CN204924934U (zh) | 基于两种量子级联激光光谱的多组分气体同时检测装置 | |
| CN211235536U (zh) | 一种燃烧场温度与气体组分浓度测试仪 | |
| CN108225543B (zh) | 一种自混合干涉的正交检波振动测量装置及方法 | |
| KR101642473B1 (ko) | 원격검출용 주파수 및 강도 변조 레이저 흡수 분광장치 및 방법 | |
| US20060262316A1 (en) | System and method for interferometric laser photoacoustic spectroscopy | |
| CN113252204B (zh) | 一种基于幅度调制的激光吸收光谱温度测量方法 | |
| CN113447458B (zh) | 一种基于激光吸收阻抗谱的气体温度与浓度参数测量方法 | |
| CN105683727B (zh) | 测量气体的气体组分的浓度的测量方法和气体分析仪 | |
| KR20060124111A (ko) | 파장 변조 방법을 이용한 유해 가스 측정 장치 | |
| CN111579532A (zh) | 一种基于频分复用和主峰扫描的激光吸收光谱层析成像方法 | |
| CN108061722A (zh) | 一种一氧化碳浓度的检测装置及检测方法 | |
| Shen et al. | Methane near-infrared laser remote detection under non-cooperative target condition based on harmonic waveform recognition | |
| CN113916802A (zh) | 一种自动校准开路式激光气体探测装置及实现方法 | |
| Kasyutich et al. | An off-axis cavity-enhanced absorption spectrometer at 1605 nm for the 12 CO 2/13 CO 2 measurement | |
| CN114397271A (zh) | 一种光谱分析温室气体的检测装置及方法 | |
| CN112729544B (zh) | 一种扫描吸收光谱的激光波长锁定系统与方法 | |
| CN115639170B (zh) | 一种基于波分复用的抗光学偏折吸收谱提取方法与系统 | |
| CN115639171A (zh) | 一种抑制乘性噪声的激光吸收光谱提取方法与系统 | |
| CN114295583A (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 |