CN114166798B - 一种基于曲率的tdlas露点测量信号增强方法 - Google Patents
一种基于曲率的tdlas露点测量信号增强方法 Download PDFInfo
- Publication number
- CN114166798B CN114166798B CN202111523155.XA CN202111523155A CN114166798B CN 114166798 B CN114166798 B CN 114166798B CN 202111523155 A CN202111523155 A CN 202111523155A CN 114166798 B CN114166798 B CN 114166798B
- Authority
- CN
- China
- Prior art keywords
- curvature
- dew point
- measurement signal
- tdlas
- point measurement
- 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
- 238000005259 measurement Methods 0.000 title claims abstract description 29
- 238000000041 tunable diode laser absorption spectroscopy Methods 0.000 title claims abstract description 19
- 238000000034 method Methods 0.000 title claims abstract description 17
- 230000002708 enhancing effect Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000000354 decomposition reaction Methods 0.000 description 8
- 238000010521 absorption reaction Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000000862 absorption spectrum Methods 0.000 description 3
- 238000009434 installation Methods 0.000 description 3
- 239000004065 semiconductor Substances 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 108010076504 Protein Sorting Signals Proteins 0.000 description 1
- 238000002835 absorbance Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000001285 laser absorption spectroscopy Methods 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
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 Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明公开了一种基于曲率的TDLAS露点测量信号增强方,首先在露点测量信号中,根据长度为N的移动窗口在第n个时刻的点的左右各取个点,得到数据序列,计算N个数值的平均值,并作为第n个时刻的数值;若数据序列是非线性分布时,首先通过公式计算第N个数据的曲率,再根据曲率通过公式对平均值进行修正;本发明针对TDLAS露点测量信号增强的问题,并考虑到数据序列的非线性特征导致滑动平均值的结果出现较大误差的问题,对数据序列的非线性特征导致滑动平均值的结果出现的较大误差进行修正,以降低噪声的影响。
Description
技术领域
本发明涉及一种增强信号的方法,具体讲是一种基于曲率的TDLAS露点测量信号增强方法。
背景技术
激光吸收光谱技术(TDLAS)具有高灵敏、连续在线、快速响应、环境适应性强等特点,是当前国际上公认的最有希望能够满足超低温风洞超低露点检测的技术手段。TDLAS技术以波长可连续调谐的窄线宽半导体激光器做光源,获得目标气体高分辨率吸收光谱来反演气体的状态参数。激光器波长可通过其有源区工作温度和驱动电流进行调谐,温度调谐响应速率相对较慢(Hz至亚Hz量级),通常采用电流调谐(最大响应速率在百kHz~MHz以上),或者采用固定激光器温度,电流调谐方式工作,波长调谐波形采用锯齿波。典型测量光路如图1所示,半导体激光器出射激光准直之后,穿过待测量气团后最终被光电探测器(PD)接收,目标气体的吸收造成激光光强衰减,衰减量遵守Lambert-Beer定律,
其中S(T)为目标气体吸收线强,对选定的吸收线仅为温度的函数;为面积归一化线型函数;P、x分别为环境静压和水汽的体积比浓度,/>即绝对分压;/>称为吸光度,描述路径上目标气体的整体吸收大小,采用锯齿波调谐激光器电流实现波长调谐的典型吸收信号及如图2所示。图2中,实线为接收的光谱强度/>,虚线为发射的光谱强度。
为了解决“倾斜”安装中存在的一些问题,采用了垂直安装方式,在垂直安装方式下,存在激光器穿过垂直于光路的玻璃窗口后的干涉问题,使测量的TDLAS信号中存在混频噪声影响测量精度。另外,随着测试对象浓度接近极限时,图2的观察信号中存在大量噪声。如何对这些噪声进行抑制,是提升测量精度的有效有段。
针对TDLAS测量信号增强问题,学者们提出了一系列方法:
[1] 奇异值分解用于可调谐二极管激光吸收光谱技术去除系统噪声;
[2] 基于经验模态分解的可调谐半导体激光吸收光谱中干涉条纹的抑制_郭心骞;
[3]基于Gabor变换的TDLAS检测信号的降噪研究_崔海滨;
在上述方法中,Gabor变换是一种人工设定标识基方法,经验模态分解、奇异值分解是一种数据驱动方法,其中,经验模态分解通过将当前信号分解为几个信号的融合,再通过线性组合重建进行去噪,分解的层次、数量等不可控;SVD分解是一种正交基方法,重建质量较低。
发明内容
针对上述问题,本发明的技术方案是:针对TDLAS露点测量信号增强的问题,并考虑到数据序列的非线性特征导致滑动平均值的结果出现较大误差的问题,提出了一种基于曲率的TDLAS露点测量信号增强方法,以降低噪声的影响。
本发明提供的一种基于曲率的TDLAS露点测量信号增强方,包括如下步骤:
首先在露点测量信号中,通过长度为N的移动窗口在第n个时刻的点的左右各取/>个点,得到数据序列/>,计算N个数值的平均值/>,并作为第n个时刻的数值;
若数据序列是非线性分布时,
首先通过公式计算第N个数据的曲率,
再根据曲率通过公式对平均值进行修正。
本发明具有如下有益效果:
本发明针对TDLAS露点测量信号增强的问题,并考虑到数据序列的非线性特征导致滑动平均值的结果出现较大误差的问题,对数据序列的非线性特征导致滑动平均值的结果出现的较大误差进行修正,以降低噪声的影响,从而增强测量信号。
附图说明
图1是TDLAS对射式测量原理示意图;
图2是典型吸收信号曲线;
图3是数据序列接近线性分布时示意图;
图4是数据序列为非线性分布时示意图;
图5是数据序列为非线性分布时表示弧长的示意图;
图6是本发明的方法流程示意图。
具体实施方式
下面将结合附图1-图6对本发明进行详细说明,对本发明实施例中的技术方案进行清楚、完整地描述。
实施例一:
一种基于曲率的TDLAS露点测量信号增强方法,以降低噪声的影响,其具体方法如下:
首先,在露点测量信号中,通过长度为N的移动窗口在第n个时刻的点的左右各取/>个点,得到数据序列/>,其中/>为中间数据;
计算N个数值的平均值,并作为第n个时刻的数值;
其中,n为整数;
其中,/>,x i 、y i 分别为通过窗口获取的N个数据中第i个数据的横坐标和纵坐标。
如图3所示,当数据序列接近线性分布时,实心圆n为测量值,通过公式和/>得到平均值为/>,对应的时刻为图中的空心圆n’,这与真实的时刻为空心正方形n相近。
实施例二:
在实施例一的基础上,当数据序列为非线性分布时,如图4所示,实心圆n为测量值,通过公式和/>得到平均值为/>,对应的时刻为图中的空心圆/>,与真实的时刻为空心正方形n相差较远,这是由于曲线弯曲引起的。
为了解决这个问题,本专利将考虑曲率的影响,首先计算N个数据的曲率,
其中,代表弧长,/>代表的就是切线方向角的变化,/>、/>和/>分别为第1个时刻、第n个时刻和第N个时刻的切线斜率。
由于数据密度较大,可以用两点之间的距离的公式来求一个近似值,如图3所示,各段线段的和为弧长,即:
其中,表示第i个时刻与第i-1个时刻的距离,/>和/>分别为第i个时刻的横坐标和纵坐标,/>和/>分别为第i-1个时刻的横坐标和纵坐标。
最后,对平均值进行修正。
其中,,/>和/>分别为常数。
本发明针对TDLAS露点测量信号增强的问题,通过计算TDLAS露点测量信号数列的平均值可以去噪,实现测量信号的增强;但是当数据序列是非线性时,即存在一定曲率时,直接通过求取平均值的话会存在一定的误差,如图4和图5,所以通过计算曲率和修正平均值,实现去噪,以增强测量信号。
Claims (1)
1.一种基于曲率的TDLAS露点测量信号增强方法,其特征在于,包括如下步骤:
在露点测量信号中,通过长度为N的移动窗口在第n个时刻的点的左右各取/>个点,得到数据序列/>,其中/>为中间数据;
计算N个数值的平均值,并作为第n个时刻的数值;
其中,n为整数;
其中,/>;
x i、y i分别为通过窗口获取的N个数据中第i个数据的横坐标和纵坐标;
考虑数据序列是非线性分布时,进行如下步骤:
S1、通过如下公式计算第N个数据的曲率,
;
其中,并代表的就是切线方向角的变化, />、/>和分别为第1个时刻、第n个时刻和第N个时刻的切线斜率,/>代表数据序列的弧长;
S2、通过如下公式对平均值进行修正,
;
其中,x n 和y n 分别表示修正后的第n个时刻的横坐标和纵坐标, , />和分别为常数;
通过如下公式计算弧长,
;
其中,表示第i个时刻与第i-1个时刻的距离,/>和/>分别为第i个时刻的横坐标和纵坐标, />和/>分别为第i-1个时刻的横坐标和纵坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111523155.XA CN114166798B (zh) | 2021-12-14 | 2021-12-14 | 一种基于曲率的tdlas露点测量信号增强方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111523155.XA CN114166798B (zh) | 2021-12-14 | 2021-12-14 | 一种基于曲率的tdlas露点测量信号增强方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114166798A CN114166798A (zh) | 2022-03-11 |
CN114166798B true CN114166798B (zh) | 2023-09-01 |
Family
ID=80486075
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111523155.XA Active CN114166798B (zh) | 2021-12-14 | 2021-12-14 | 一种基于曲率的tdlas露点测量信号增强方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114166798B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115326349B (zh) * | 2022-10-17 | 2023-02-03 | 中国科学院合肥物质科学研究院 | 一种低温风洞试验段露点测量测试的光路系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2016070686A (ja) * | 2014-09-26 | 2016-05-09 | 神栄テクノロジー株式会社 | Tdlas法による濃度測定装置。 |
CN108709871A (zh) * | 2018-08-06 | 2018-10-26 | 电子科技大学 | 一种基于tdlas的气体浓度检测方法 |
CN111829980A (zh) * | 2020-07-23 | 2020-10-27 | 安徽农业大学 | 一种基于谐波技术的线型非线性修正的检测系统及方法 |
-
2021
- 2021-12-14 CN CN202111523155.XA patent/CN114166798B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2016070686A (ja) * | 2014-09-26 | 2016-05-09 | 神栄テクノロジー株式会社 | Tdlas法による濃度測定装置。 |
CN108709871A (zh) * | 2018-08-06 | 2018-10-26 | 电子科技大学 | 一种基于tdlas的气体浓度检测方法 |
CN111829980A (zh) * | 2020-07-23 | 2020-10-27 | 安徽农业大学 | 一种基于谐波技术的线型非线性修正的检测系统及方法 |
Non-Patent Citations (1)
Title |
---|
可调谐二极管激光吸收光谱技术测量低温流场水汽露点温度;聂伟 等;光学精密工程;第26卷(第08期);1862-1869 * |
Also Published As
Publication number | Publication date |
---|---|
CN114166798A (zh) | 2022-03-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110441777A (zh) | 一种基于激光雷达的气溶胶垂直廓线的反演方法 | |
CN107389606B (zh) | 一种基于可调谐半导体激光吸收谱的光程长度分析方法 | |
CN105628649B (zh) | 气体现场监测特征吸收峰提取方法 | |
US20220026577A1 (en) | Dispersion gating-based atmospheric composition measurement laser radar | |
CN109991189B (zh) | 一种基于波数漂移修正的固定点波长调制气体浓度测量装置及其测量方法 | |
CN114166798B (zh) | 一种基于曲率的tdlas露点测量信号增强方法 | |
CN102692273B (zh) | 一种干涉型高光谱成像仪的mtf在轨检测方法 | |
CN104076345B (zh) | 一种针对纯转动拉曼激光雷达测量温度的饱和修正方法 | |
CN108956491A (zh) | 一种大气hcho遥感监测方法 | |
CN114384045A (zh) | 一种痕量气体浓度和路径长度实时检测系统与方法 | |
CN110471046B (zh) | 一种差分吸收激光雷达系统常数定标方法 | |
Joly et al. | Development of a compact CO 2 sensor open to the atmosphere and based on near-infrared laser technology at 2.68 μm | |
KR101968087B1 (ko) | 초분광 이미징 센서 관측 자료를 이용한 대기중 오존 및 이산화질소의 기주 농도 산출 장치 | |
CN108333143B (zh) | 一种基于可调谐激光吸收光谱的水汽浓度测量修正方法 | |
Pedireddi et al. | Characterization of atmospheric turbulence effects and their mitigation using wavelet-based signal processing | |
Gong et al. | Multi-points scanning method for wavelength locking in CO2 differential absorption lidar | |
Li et al. | Retrieval of the aerosol extinction coefficient of 1064nm based on high-spectral-resolution lidar | |
Zhu et al. | Improved Savitzky-Golay filtering algorithm for measuring a pharmaceutical vial’s oxygen content based on wavelength modulation spectroscopy | |
CN113390826B (zh) | 一种用于视线非均匀吸收光谱测量的积分吸光度测量误差估计方法 | |
Bilbe et al. | An improved Raman lidar system for the remote measurement of natural gas releases into the atmosphere | |
Notholt et al. | Test of in situ measurements of atmospheric aerosols and trace gases by long path transmission spectroscopy | |
Ma et al. | Noise reduction for ground-based atmospheric detection lidar: A universal method based on signal segmentation and reconstruction | |
CN108663328B (zh) | 适用于光学水质测量系统的光噪音滤除装置 | |
Yang et al. | Performance analysis of moving average filtering used in TDLAS carbon dioxide sensing system | |
Song et al. | Performance enhancement of methane detection using a novel self-adaptive mid-infrared absorption spectroscopy technique |
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 |