CN117387673B - 一种基于参考气室的光纤光栅解调方法及设备 - Google Patents
一种基于参考气室的光纤光栅解调方法及设备 Download PDFInfo
- Publication number
- CN117387673B CN117387673B CN202311675037.XA CN202311675037A CN117387673B CN 117387673 B CN117387673 B CN 117387673B CN 202311675037 A CN202311675037 A CN 202311675037A CN 117387673 B CN117387673 B CN 117387673B
- Authority
- CN
- China
- Prior art keywords
- array
- wavelength
- spectrum
- acetylene
- grating
- 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
- 239000000835 fiber Substances 0.000 title claims abstract description 85
- 238000000034 method Methods 0.000 title claims abstract description 49
- HSFWRNGVRCDJHI-UHFFFAOYSA-N alpha-acetylene Natural products C#C HSFWRNGVRCDJHI-UHFFFAOYSA-N 0.000 claims abstract description 99
- 125000002534 ethynyl group Chemical group [H]C#C* 0.000 claims abstract description 99
- 238000010521 absorption reaction Methods 0.000 claims abstract description 87
- 238000001228 spectrum Methods 0.000 claims abstract description 86
- 238000005070 sampling Methods 0.000 claims abstract description 25
- 238000000411 transmission spectrum Methods 0.000 claims abstract description 17
- 238000012545 processing Methods 0.000 claims abstract description 11
- 238000003491 array Methods 0.000 claims description 20
- 238000012216 screening Methods 0.000 claims description 14
- 238000010606 normalization Methods 0.000 claims description 10
- 239000007789 gas Substances 0.000 description 17
- LELOWRISYMNNSU-UHFFFAOYSA-N hydrogen cyanide Chemical compound N#C LELOWRISYMNNSU-UHFFFAOYSA-N 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 239000013307 optical fiber Substances 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 238000000862 absorption spectrum Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000000985 reflectance spectrum Methods 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 231100000331 toxic Toxicity 0.000 description 1
- 230000002588 toxic effect Effects 0.000 description 1
- 239000002341 toxic gas Substances 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D5/00—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable
- G01D5/26—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light
- G01D5/32—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light
- G01D5/34—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells
- G01D5/353—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre
-
- 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
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明涉及光栅解调技术领域，具体为一种基于参考气室的光纤光栅解调方法及设备，包括采集乙炔透射光谱信号，通过拟合透射光谱基线获取乙炔差分光谱；通过采样率与扫描信号周期计算扫描信号周期的采样点个数，根据乙炔气体吸收峰波长确定扫描信号周期范围内的乙炔吸收峰个数，并建立对应乙炔吸收峰的位置索引变量数组；基于乙炔差分光谱与乙炔吸收峰个数，确定位置索引变量数组对应的绝对波长；采集光纤光栅反射谱处理得到归一化光栅光谱数组，通过光栅寻峰法得到光纤光栅中心波长对应的索引值数组与索引值；将光纤光栅中心波长的索引值代入至位置索引变量数组中，确定光纤光栅的中心波长。
Description
技术领域
本发明涉及光栅解调技术领域，具体为一种基于参考气室的光纤光栅解调方法及设备。
背景技术
现代社会工业生产和基础建筑设施中，温度、应力是两个相当重要的物理量，生产设备的工作温度变化及关键运行部位的应力变化能够反映其运行状态是否正常；基础建筑结构的环境温度的变化或关键支撑结构应力的变化则影响其结构损伤和使用寿命。随着国民经济的迅速发展，许多现代工业如电力系统、石油化工等的生产规模不断扩大而政府对基础建筑设施如桥梁、隧道、大坝的投入建设也显著增加，这些生产与工程领域的建设在给社会和国民带来的经济效益、生活便利的同时，也存在人为失误破坏、环境变迁影响导致生产设备损坏和建筑安全失效的风险。
光纤光栅传感技术作为光纤传感技术的一种，凭借其自身独特的优点，在温度和应力检测领域得到了迅速发展。光纤光栅波长解调方法的研究一直是光纤光栅传感领域研究的热点，太原理工大学硕士论文《基于扫描激光器的光纤光栅解调仪研究》中采用40nm扫描激光器模组作为光源，通过控制激光器脉冲调节波长和捕获上升索引实现波长的解调，该方案采用完全成熟的激光器模组，硬件成本高，并且绝对波长完全依靠激光器性能本身，当激光器波长发生偏移时，由于没有绝对波长参考，便容易造成波长解调不准确。中国专利文献CN108106645A提供了《基于氰化氢吸收波长参考的光纤光栅传感解调装置与方法》中提到采用氰化氢气体作为参考气体为光纤光栅解调提供绝对波长参考，但是由于解调系统采用ASE宽带光源、PF腔调节以及剧毒气体氰化氢做波长参考，使得该系统成本较大，批量化、产业化困难。
发明内容
为解决上述的问题，本发明一方面提供了一种基于参考气室的光纤光栅解调方法，包括以下步骤：
S1、采集乙炔透射光谱信号，通过拟合透射光谱基线获取乙炔差分光谱；
S2、通过采样率与扫描信号周期计算扫描信号周期的采样点个数，根据乙炔气体吸收峰波长确定扫描波长范围内包含的乙炔吸收峰个数，并建立对应乙炔吸收峰的位置索引变量数组；
S3、基于乙炔差分光谱与乙炔吸收峰个数，确定位置索引变量数组对应的绝对波长；
S4、采集光纤光栅反射谱处理得到归一化光栅光谱数组，通过光栅寻峰法得到光纤光栅中心波长对应的索引值数组与索引值；
S5、将光纤光栅中心波长的索引值代入至位置索引变量数组中，确定光纤光栅的中心波长。
作为具体的实施方式，所述S3中基于乙炔差分光谱与乙炔吸收峰个数，确定位置索引变量数组对应的绝对波长，具体方法为：
S3.1、将获取的乙炔差分光谱进行最大最小值归一化得到差分归一化光谱数组；
S3.2、将差分归一化光谱数组与预设阈值进行比较，确定相邻差分数组；将所述相邻差分数组中各元素与预设阈值进行比较处理，得到潜在吸收峰两翼的位置数组；
S3.3、基于乙炔吸收峰波长和扫描的最小步长，对位置数组进行筛选处理确定位置索引变量数组的位置索引值；
S3.4、将位置索引值与乙炔吸收峰波长数组进行线性拟合，计算索引值数组各位置索引值对应的波长，并将各段拟合得到的波长拼接得到位置索引变量数组各元素点对应的绝对波长。
所述S3.2中通过将差分归一化光谱数组与预设阈值进行比较，确定相邻差分数组，具体为：
将差分归一化光谱数组各元素与预设阈值进行比较，若大于预设阈值，差分归一化光谱数组的数值保留；若小于等于预设阈值，则置零；
将保留的差分归一化光谱数组做相邻元素差分处理，得到相邻差分数组，并对所述相邻差分数组各元素数值取绝对值。
所述S3.3中基于乙炔吸收峰个数和扫描的最小步长，对位置数组进行筛选处理确定位置索引变量数组的位置索引值，具体包括：
S3.3.1、根据相邻乙炔吸收峰的波长间隔，结合扫描波长范围与采样点个数，预估步长分辨率对位置数组进行筛选处理，确定完整吸收峰的两翼；
S3.3.2、根据确定的完整吸收峰的两翼，基于步长，计算两翼的平均值作为当前吸收峰的位置索引值，得到索引值数组。
作为具体的实施方式，所述S3.3之后还包括判断索引值数组中位置索引值的个数与乙炔吸收峰个数是否相等，
若相等，则表示确定得到完整的乙炔差分光谱对应的位置索引值；
若不相等，则改变预设阈值并重复步骤S3.2-步骤S3.3。
作为具体的实施方式，所述扫描波长范围为1531.19-1535.87nm。
作为具体的实施方式，所述S4中采集光纤光栅光谱反射谱处理得到归一化光栅光谱数组，具体包括：通过采集光纤光栅反射谱，经高斯卷积后得到光栅光谱数组；对所述光栅光谱数组进行最大最小值归一化得到归一化光栅光谱数组。
所述S4中通过光栅寻峰法得到光纤光栅中心波长对应的索引数组与索引值，具体方法为：
S4.1、通过将归一化光栅光谱数组与固定阈值进行比较，确定相邻差分数组；遍历所述相邻差分数组中的元素与固定阈值进行比较处理，得到潜在吸收峰两翼的位置数组；
S4.2、基于扫描的最小步长，对位置数组进行筛选处理确定光纤光栅中心波长对应的索引数组与索引值。
所述光纤光栅解调方法还包括所述S5在确定光纤光栅的中心波长后，将光纤光栅的中心波长与初始波长比对确定波长的变化量，根据确定的波长的变化量确定光纤光栅阵列变化的物理量。
本发明另一方面提供了一种基于参考气室的光纤光栅解调设备，包括处理器和存储器，其中，所述处理器执行所述存储器中保存的程序数据时实现如上所述的基于参考气室的光纤光栅解调方法。
有益效果在于：
（1）本申请通过寻气体吸收峰算法结合乙炔气体在1532nm附近的吸收特征，自动改变判定阈值根据吸收峰个数分析出指定吸收峰波长对应的采集信号的位置索引，为波长-索引分段拟合提供依据；
（3）本申请通过分段线性拟合将相邻乙炔吸收峰的波长与位置索引作为线性拟合的两点，通过遍历波长与位置索引值，进而得到了整个采样点数与波长的一一对应关系，进一步的将光纤光栅的解调波长与采样精度强相关，并且乙炔气体吸收峰波长不受温度、压强等环境因素影响，实现以位置索引求解的光栅波长稳定性好、准确度高；
（3）本申请采用最大最小值归一化，消除设备、通道、接头等因素对乙炔吸收谱和光纤光栅反射谱的强度影响；
（4）采用乙炔气体参考气室得到了绝对波长，其波长的获取不受环境温度、工作时长、环境湿度、压强等外界因素影响降低了光纤光栅解调仪获取标准波长的难度，提高了解调准确度，此外，相比于以含有剧毒且易燃易爆的氰化氢参考气室，乙炔参考气室不仅成本便宜，还大大提高了设备安全性。
附图说明
图1为基于参考气室的光纤光栅解调方法的流程图；
图2为光纤光栅解调方法的步骤S3.1-S3.2的流程图；
图3为光纤光栅解调方法步骤S3.3的流程图；
图4为乙炔气体参考气室的位置索引波长图。
具体实施方式
下面将结合附图更详细地描述本公开的示例性实施方式。
实施例1，本实施例提供了一种基于参考气室的光纤光栅解调方法，应用于光纤光栅解调系统，所述光纤光栅解调系统包括VCSEL激光器、激光器驱动模块、一分八耦合器、一分二耦合器、光纤光栅阵列、乙炔参考气室、光电探测器阵列和数据采集处理器。其中，该光纤光栅解调系统在进行工作时，由STM32F407为核心处理器的数据采集处理器产生频率为10Hz，幅值为1V的正锯齿波传输到激光器驱动模块，所述激光器驱动模块将该信号源进行放大，并加载到VCSEL激光器的波长调制引脚产生周期为10Hz的扫描波长，所述扫描波长的波长范围为1531.19-1535.87nm；所述VCSEL激光器出射光经光纤耦合到一分八耦合器上，所述一分八耦合器的输出通道一经过光纤耦合到密封乙炔气体的参考气室，所述乙炔参考气室的透射光信号经光纤传输至光电探测器阵列得到乙炔透射谱电压信号。所述一分八耦合器的其他通道经过光纤耦合至一分二耦合器的通道一，所述一分二耦合器的公共端连接光纤光栅阵列，出射光经布拉格反射，光纤光栅阵列的反射谱再次经过一分二耦合器的通道二被光电探测器检测，得到光栅的反射谱电压信号。所述数据采集处理器同步采集的乙炔透射光谱信号和光纤光栅反射谱信号进行光栅解调。
所述基于参考气室的光纤光栅解调方法的具体步骤如下：
S1、采集乙炔透射光谱信号，通过拟合透射光谱基线获取乙炔差分光谱。
通过乙炔透射光谱信号将未发生吸收的波形拟合得到透射光谱基线，将透射光谱基线与乙炔透射光谱信号做差，获取得到乙炔差分光谱。
S2、通过采样率与扫描信号周期计算扫描信号周期的采样点个数，根据乙炔气体吸收峰波长确定扫描波长范围内包含的乙炔吸收峰个数，并建立对应乙炔吸收峰的位置索引变量数组。
在进行光纤光栅解调处理时，载入信号的ADC采样率fs与扫描信号周期f，根据采样率与扫描信号周期计算得到一个完整周期的采样点个数为N=fs/f；同时根据乙炔气体吸收峰波长建立乙炔吸收峰波长数组。如下：
[1535.394,1534.744,1534.097,1533.467,1532.828,1532.205]；
根据乙炔吸收峰波长数组确定扫描波长范围内的乙炔吸收峰个数，并建立对应乙炔吸收峰的位置索引变量数组，所述位置索引变量数组的个数为N。
S3、基于乙炔差分光谱与乙炔吸收峰个数，确定位置索引变量数组对应的绝对波长。
其具体方法包括：
S3.1、将获取的乙炔差分光谱进行最大最小值归一化得到差分归一化光谱数组；
将乙炔差分光谱的最大吸收峰缩放到绝对数值1，基线缩放至0附近，并根据归一化公式得到差分归一化光谱数组，其具体公式如下：
；
其中，X norm表示差分归一化光谱数值，无量纲；X表示乙炔差分光谱的电压值；X min表示乙炔差分光谱电压值的最小值；X max表示乙炔差分光谱电压值的最大值；
通过对乙炔差分光谱进行最大最小值归一化，能够消除设备、通道、接头等因素对乙炔吸收透射谱的强度影响。
S3.2、将差分归一化光谱数组与预设阈值进行比较，确定相邻差分数组；将所述相邻差分数组中各元素与预设阈值进行比较处理，得到潜在吸收峰两翼的位置数组。
其中，所述通过将差分归一化光谱数组与预设阈值进行比较，确定相邻差分数组，具体为：
将差分归一化光谱数组各元素与预设阈值进行比较，若大于预设阈值，差分归一化光谱数组的数值保留；若小于等于预设阈值，则置零；
将保留的差分归一化光谱数组做相邻元素差分处理，得到相邻差分数组，并对所述相邻差分数组各元素数值取绝对值；
传统的峰值索引解调是通过监测最大值的索引确定，但是由于采样率和噪声限制直接获取最大值使得解调误差较大，因此更加常见的方法是3db带宽法也即取最大值的半高宽对应的索引再求均值；本发明中通过设定阈值比较的方式，将“半高宽”的索引值确定，以便后续确定是否为真实峰的索引。分别将所述相邻差分数组中各元素与预设阈值进行比较，将大于该预设阈值的元素保留，并作为乙炔吸收峰对应的位置索引，建立潜在吸收峰两翼的位置数组。
S3.3、基于乙炔吸收峰波长和扫描最小步长，对位置数组进行筛选处理确定位置索引变量数组的位置索引值，具体包括：
S3.3.1、根据相邻乙炔吸收峰的波长间隔，结合扫描波长范围与采样点个数，预估步长分辨率对位置数组进行筛选处理，确定完整吸收峰的两翼；
根据乙炔气体在1532附近的吸收谱线可以得出大吸收峰峰值的半高全宽约为30pm，小吸收峰峰值的半高宽约为40pm，因此确定单个吸收峰的两翼索引差值不会太大。
根据乙炔吸收峰之间的0.6nm波长间隔，只要不把前一个吸收峰的左侧索引与后一个吸收峰的右侧索引判定是一个吸收峰的两翼均能满足要求，因此结合激光器扫描的总波长范围与采样点个数共同预估采样步长分辨率Δλ，确定筛选。所述筛选的位置阈值为10*6*Δλ。
对位置数组做相邻元素差分处理，若小于位置阈值，则确定所述位置数组的相邻两个元素为一个完整吸收峰的两翼。
S3.3.2、根据确定的完整吸收峰的两翼，基于步长，计算两翼的平均值作为当前吸收峰的位置索引值，得到索引值数组。
同时，步骤S3.3还包括判断索引值数组中位置索引值的个数与乙炔吸收峰个数是否相等，
若相等，则表示确定得到完整的乙炔差分光谱对应的位置索引值；
若不相等，则改变预设阈值并重复步骤S3.2-步骤S3.3。
S3.4、将位置索引值与乙炔吸收峰波长数组进行线性拟合，计算索引值数组各位置索引值对应的波长，并将各段拟合得到的波长拼接得到位置索引变量数组各元素点对应的绝对波长。
将得到的吸收峰峰值的位置索引值与乙炔吸收峰波长形成数值对，每相邻两组数值对做线性拟合，将各段拟合得到的波长拼接得到位置索引变量数组各元素点对应的绝对波长。
S4、采集光纤光栅反射谱处理得到归一化光栅光谱数组，通过光栅寻峰法得到光纤光栅中心波长对应的索引值数组与索引值。
通过采集光纤光栅反射谱，经高斯卷积后得到光栅光谱数组，所述高斯函数系数根据气体吸收峰的个数、采样率和扫描周期等参数进行确定；并对所述光栅光谱数组进行最大最小值归一化得到归一化光栅光谱数组。
在得到归一化光栅光谱数组之后，通过光栅寻峰法得到光纤光栅中心波长对应的索引值数组与索引值，所述光栅寻峰法与步骤S3乙炔寻峰方法类似，具体方法为：
S4.1、通过将归一化光栅光谱数组与固定阈值（因光纤光栅反射谱的信噪比较好，即阈值可以采用固定阈值）进行比较，所述固定阈值可以手动设置更改，确定相邻差分数组；遍历所述相邻差分数组中的元素与固定阈值进行比较处理，得到潜在吸收峰两翼的位置数组；
S4.2、基于扫描的最小步长，对位置数组进行筛选处理确定光纤光栅中心波长对应的索引数组与索引值。
S5、将光纤光栅中心波长的索引值代入至位置索引变量数组中，确定光纤光栅的中心波长。
在确定光纤光栅的中心波长后，将光纤光栅的中心波长与初始波长比对确定波长的变化量，根据确定的波长的变化量确定光纤光栅阵列变化的物理量，实现光纤光栅波长信号的解调。
其次，本申请还提供了一种基于参考气室的光纤光栅解调设备，包括处理器和存储器，其中，所述处理器执行所述存储器中保存的程序数据时实现如上所述的基于参考气室的光纤光栅解调方法。
实施例2，参见图1-图3，本实施例提供了一种基于参考气室的光纤光栅解调方法，包括：
S1、采集乙炔透射光谱信号，并拟合透射光谱基线获取乙炔差分光谱one_period_diff。
S2、载入相关参数，包括光纤光栅信号的ADC采样率fs=12500，扫描信号周期f_scan=10hz,乙炔差分光谱one_period_diff，光纤光栅反射谱one_grating_reflection，以及乙炔吸收峰波长数组real_wavelength；
real_wavelength=[1535.394,1534.744,1534.097,1533.467,1532.828,1532.205]；
通过当前的ADC采样率fs与扫描信号周期f_scan计算扫描信号周期的采样点个数为N=fs/f_scan，根据乙炔气体吸收峰波长确定扫描波长范围内包含的乙炔吸收峰个数real_waveleng_num = len(real_wavelength)，并建立对应乙炔吸收峰的位置索引变量数组one_x=[0:N]，数组长度从0到N-1，共N个数。
S3、基于乙炔差分光谱one_period_diff与乙炔吸收峰个数real_waveleng_num，确定位置索引变量数组one_x对应的绝对波长。
S3.1、对乙炔差分光谱one_period_diff进行最大最小值归一化得到差分归一化光谱数组one_diff_norm；
分别得到one_period_diff的最小值min和最大值max，
one_diff_norm=（one_period_diff-min）/(max-min)。
S3.2、通过将差分归一化光谱数组与预设阈值进行比较，确定相邻差分数组；将所述相邻差分数组中各元素与预设阈值进行比较处理，得到潜在吸收峰两翼的位置数组；
给定预设阈值threshold_value=0.5-0.05*j，for j=0；n＜10；j++；
将one_diff_norm复制给high_threshold，遍历数组high_threshold，对数组中各元素与threshold_value分别进行比较，若大于预设阈值，差分归一化光谱数组的数值保留；若小于等于预设阈值，则置零；
将保留的high_threshold做相邻元素差分处理，得到相邻差分数组high_threshold_slope，并对所述相邻差分数组各元素数值取绝对值；
high_threshold_slope[n]=high_threshold[n+1]-high_threshold[n],
遍历所述相邻差分数组中各元素与预设阈值进行比较，将大于该预设阈值的元素进行保留，并作为乙炔吸收峰对应的位置索引，建立潜在吸收峰两翼的位置数组slope_max_point。
S3.3、基于乙炔吸收峰个数和扫描的最小步长，对位置数组进行筛选处理确定位置索引变量数组的位置索引值。
S3.3.1、根据乙炔吸收峰之间的0.6nm波长间隔，结合激光器扫描的总波长范围与采样点个数共同预估采样步长分辨率Δλ，对位置数组slope_max_point进行筛选，只有满足slope_max_point[k+1]-slope_max_point[k]＜int(10*6*Δλ),for k=0;k＜<len(slope_max_point)；k++的元素为一个完整的吸收峰的两翼，并添加到数组main_peak_index_wings中。
S3.3.2、将数组main_peak_index_wings的元素，以步长为2，计算两翼的平均值并取值作为当前吸收峰的位置索引值data_main_peak_index；
data_main_peak_index[n]=int(0.5*(main_peak_index_wings[2*n]+main_peak_index_wings[2*n+1])),得到索引值数组main_peak_index；
并判断索引值数组main_peak_index的位置索引值data_main_peak_index的个数与乙炔吸收峰个数real_waveleng_num是否相等，
若相等，则表示确定得到完整的乙炔差分光谱对应的位置索引值；
若不相等，则改变预设阈值threshold_value并重复步骤S3.2-步骤S3.3。
S3.4、将位置索引值与乙炔吸收峰进行线性拟合，得到位置索引值对应的绝对波长；
创建波长数组all_area_cal_wavelength用来存放拟合峰值，对索引值数组进行取点main_peak_index[n],main_peak_index[n+1]构成拟合的X，对乙炔吸收峰波长数组取点real_wavelength[n],real_wavelength[n+1]构成拟合的Y，进行线性拟合，得到常数项a0，一次项a1；通过从位置索引变量数组one_x中截取从main_peak_index[n]到main_peak_index[n+1]的序列作为a1*x+a0的x，代入并得到当前序列对应的计算波长are_cal_wavelength；
将每一段拟合得到的计算波长拼接到波长数组all_area_cal_wavelength中，便得到了位置索引变量数组one_x每个元素点对应的绝对波长，如图4所示。
S4、通过对光纤光栅反射谱one_grating_reflection进行高斯卷积，得到光栅光谱数组one_grating_reflection_gauss；并对one_grating_reflection_gauss做最大最小值归一化，得到归一化光栅光谱数组one_grating_reflection_gauss_renorm，并且光纤光栅反射谱的信噪比较好，即光栅寻峰法的预设阈值采用固定阈值，光栅透射两翼的固定阈值根据同一通道光栅个数进行设置；
通过光栅寻峰法得到光纤光栅中心波长对应的索引值数组grating_peak_index，所述光栅寻峰法与步骤S3乙炔寻峰算法类似，在此不再赘述。
S5、将得到的光纤光栅中心波长的索引值数组grating_peak_index带入至波长数组all_area_cal_wavelength中，确定光纤光栅的中心波长。
Claims (4)
1.一种基于参考气室的光纤光栅解调方法，其特征在于，包括以下步骤：
S1、采集乙炔透射光谱信号，通过拟合透射光谱基线获取乙炔差分光谱；
S2、通过采样率与扫描信号周期计算扫描信号周期的采样点个数，根据乙炔气体吸收峰波长确定扫描波长范围内包含的乙炔吸收峰个数，并建立对应乙炔吸收峰的位置索引变量数组；
S3、基于乙炔差分光谱与乙炔吸收峰个数，确定位置索引变量数组对应的绝对波长，具体方法为：
S3.1、将获取的乙炔差分光谱进行最大最小值归一化得到差分归一化光谱数组；
S3.2、将差分归一化光谱数组与预设阈值进行比较，确定相邻差分数组；将所述相邻差分数组中各元素与预设阈值进行比较处理，得到潜在吸收峰两翼的位置数组；
通过将差分归一化光谱数组与预设阈值进行比较，确定相邻差分数组，具体为：将差分归一化光谱数组各元素与预设阈值进行比较，若大于预设阈值，差分归一化光谱数组的数值保留；若小于等于预设阈值，则置零；将保留的差分归一化光谱数组做相邻元素差分处理，得到相邻差分数组，并对所述相邻差分数组各元素数值取绝对值；
S3.3、基于乙炔吸收峰波长和扫描的最小步长，对位置数组进行筛选处理确定位置索引变量数组的位置索引值，具体包括：
S3.3.1、根据相邻乙炔吸收峰的波长间隔，结合扫描波长范围与采样点个数，预估步长分辨率对位置数组进行筛选处理，确定完整吸收峰的两翼；
S3.3.2、根据确定的完整吸收峰的两翼，基于步长，计算两翼的平均值作为当前吸收峰的位置索引值，得到索引值数组；
所述S3.3之后还包括判断索引值数组中位置索引值的个数与乙炔吸收峰个数是否相等，若相等，则表示确定得到完整的乙炔差分光谱对应的位置索引值；
若不相等，则改变预设阈值并重复步骤S3.2-步骤S3.3；
S3.4、将位置索引值与乙炔吸收峰波长数组进行线性拟合，计算索引值数组各位置索引值对应的波长，并将各段拟合得到的波长拼接得到位置索引变量数组各元素点对应的绝对波长；
S4、采集光纤光栅反射谱处理得到归一化光栅光谱数组，通过光栅寻峰法得到光纤光栅中心波长对应的索引值数组与索引值；
所述采集光纤光栅反射谱处理得到归一化光栅光谱数组，具体包括：通过采集光纤光栅反射谱，经高斯卷积后得到光栅光谱数组；对所述光栅光谱数组进行最大最小值归一化得到归一化光栅光谱数组；
所述通过光栅寻峰法得到光纤光栅中心波长对应的索引数组与索引值，具体方法为：
S4.1、通过将归一化光栅光谱数组与固定阈值进行比较，确定相邻差分数组；遍历所述相邻差分数组中的元素与固定阈值进行比较处理，得到潜在吸收峰两翼的位置数组；
S4.2、基于扫描的最小步长，对位置数组进行筛选处理确定光纤光栅中心波长对应的索引数组与索引值；
S5、将光纤光栅中心波长的索引值代入至位置索引变量数组中，确定光纤光栅的中心波长。
2.根据权利要求1所述的基于参考气室的光纤光栅解调方法，其特征在于，所述扫描波长范围为1531.19-1535.87nm。
3.根据权利要求1所述的基于参考气室的光纤光栅解调方法，其特征在于，还包括所述S5在确定光纤光栅的中心波长后，将光纤光栅的中心波长与初始波长比对确定波长的变化量，根据确定的波长的变化量确定光纤光栅阵列变化的物理量。
4.一种基于参考气室的光纤光栅解调设备，其特征在于，包括处理器和存储器，其中，所述处理器执行所述存储器中保存的程序数据时实现如所述权利要求1-3中任一项所述的基于参考气室的光纤光栅解调方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311675037.XA CN117387673B (zh) | 2023-12-08 | 2023-12-08 | 一种基于参考气室的光纤光栅解调方法及设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311675037.XA CN117387673B (zh) | 2023-12-08 | 2023-12-08 | 一种基于参考气室的光纤光栅解调方法及设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117387673A CN117387673A (zh) | 2024-01-12 |
| CN117387673B true CN117387673B (zh) | 2024-02-23 |
Family
ID=89437672
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311675037.XA Active CN117387673B (zh) | 2023-12-08 | 2023-12-08 | 一种基于参考气室的光纤光栅解调方法及设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117387673B (zh) |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2377497C1 (ru) * | 2008-07-15 | 2009-12-27 | Федеральное государственное унитарное предприятие "ВСЕРОССИЙСКИЙ НАУЧНО-ИССЛЕДОВАТЕЛЬСКИЙ ИНСТИТУТ ОПТИКО-ФИЗИЧЕСКИХ ИЗМЕРЕНИЙ" (ФГУП "ВНИИОФИ") | Устройство для измерения деформаций на основе квазираспределенных волоконно-оптических датчиков на брэгговских решетках |
| CN101799331A (zh) * | 2010-03-16 | 2010-08-11 | 山东大学 | 基于乙炔气体吸收峰的dfb激光器波长解调方法和系统 |
| CN104215605A (zh) * | 2014-02-12 | 2014-12-17 | 中科融通物联科技无锡有限公司 | 一种快速的光纤光栅解调系统的动态定标方法 |
| CN104931081A (zh) * | 2015-06-10 | 2015-09-23 | 天津大学 | 基于复合波长参考的光纤光栅传感解调装置及方法 |
| CN106248121A (zh) * | 2016-08-11 | 2016-12-21 | 天津大学 | 环境变温下波动抑制的光纤光栅传感解调装置与解调方法 |
| CN108106645A (zh) * | 2017-12-19 | 2018-06-01 | 天津大学 | 基于氰化氢吸收波长参考的光纤光栅传感解调装置与方法 |
| CN110736708A (zh) * | 2019-09-23 | 2020-01-31 | 天津大学 | 一种变温环境下基于恢复的fbg高精度解调装置及解调方法 |
| CN114235018A (zh) * | 2021-12-09 | 2022-03-25 | 山东微感光电子有限公司 | 一种温度自适应的fbg解调方法及系统 |
| CN114527488A (zh) * | 2022-04-25 | 2022-05-24 | 中国科学技术大学 | 1548nm多普勒测风激光雷达发射激光波长绝对锁定装置及方法 |
| CN115839916A (zh) * | 2022-12-20 | 2023-03-24 | 山东星冉信息科技有限公司 | 一种免疫激光器温度变化的气体检测报警系统及方法 |
| CN116930121A (zh) * | 2023-08-22 | 2023-10-24 | 山东微感光电子有限公司 | 基于vcsel和乙炔吸收峰的光纤光栅解调方法与装置 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8358417B2 (en) * | 2010-10-21 | 2013-01-22 | Spectrasensors, Inc. | Spectrometer with validation cell |
| CN113686368A (zh) * | 2021-08-24 | 2021-11-23 | 中北大学 | 用于获得相位变化参数的光纤传感器解调系统和解调方法 |
-
2023
- 2023-12-08 CN CN202311675037.XA patent/CN117387673B/zh active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2377497C1 (ru) * | 2008-07-15 | 2009-12-27 | Федеральное государственное унитарное предприятие "ВСЕРОССИЙСКИЙ НАУЧНО-ИССЛЕДОВАТЕЛЬСКИЙ ИНСТИТУТ ОПТИКО-ФИЗИЧЕСКИХ ИЗМЕРЕНИЙ" (ФГУП "ВНИИОФИ") | Устройство для измерения деформаций на основе квазираспределенных волоконно-оптических датчиков на брэгговских решетках |
| CN101799331A (zh) * | 2010-03-16 | 2010-08-11 | 山东大学 | 基于乙炔气体吸收峰的dfb激光器波长解调方法和系统 |
| CN104215605A (zh) * | 2014-02-12 | 2014-12-17 | 中科融通物联科技无锡有限公司 | 一种快速的光纤光栅解调系统的动态定标方法 |
| CN104931081A (zh) * | 2015-06-10 | 2015-09-23 | 天津大学 | 基于复合波长参考的光纤光栅传感解调装置及方法 |
| CN106248121A (zh) * | 2016-08-11 | 2016-12-21 | 天津大学 | 环境变温下波动抑制的光纤光栅传感解调装置与解调方法 |
| CN108106645A (zh) * | 2017-12-19 | 2018-06-01 | 天津大学 | 基于氰化氢吸收波长参考的光纤光栅传感解调装置与方法 |
| CN110736708A (zh) * | 2019-09-23 | 2020-01-31 | 天津大学 | 一种变温环境下基于恢复的fbg高精度解调装置及解调方法 |
| CN114235018A (zh) * | 2021-12-09 | 2022-03-25 | 山东微感光电子有限公司 | 一种温度自适应的fbg解调方法及系统 |
| CN114527488A (zh) * | 2022-04-25 | 2022-05-24 | 中国科学技术大学 | 1548nm多普勒测风激光雷达发射激光波长绝对锁定装置及方法 |
| CN115839916A (zh) * | 2022-12-20 | 2023-03-24 | 山东星冉信息科技有限公司 | 一种免疫激光器温度变化的气体检测报警系统及方法 |
| CN116930121A (zh) * | 2023-08-22 | 2023-10-24 | 山东微感光电子有限公司 | 基于vcsel和乙炔吸收峰的光纤光栅解调方法与装置 |
Non-Patent Citations (6)
| Title |
|---|
| High resolution phase demodulation at 1530nm using low-pressure acetylene absorption lines;Diaz, J.;Proceedings of SPIE;20210101;第11700卷;全文 * |
| 基于光纤传感的航天器结构应变参数精确测量;杨宁;申景诗;张建德;徐安宏;;航天器环境工程;20201025(05);全文 * |
| 基于气室的光纤光栅反射谱高精度解调方法研究;鞠婉秋;中国优秀硕士学位论文全文数据库信息科技辑;20220315(第03（2022）期);全文 * |
| 基于直接吸收光谱的高浓度气体检测研究及甲烷气体检测样机开发;张志峰;中国优秀硕士学位论文全文数据库工程科技Ⅰ辑;20230215(第02(2023)期);全文 * |
| 波长扫描及多次平均法提高气体传感信号信噪比的理论与实验研究;齐洁;董小鹏;郑俊达;林淦斌;;江苏现代计量;20110228(01);全文 * |
| 高速高分辨率光纤布拉格光栅传感系统的解调技术;王拥军;刘永超;张靖涛;杨昭怡;王智;;中国激光;20130210(02);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117387673A (zh) | 2024-01-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20210025945A1 (en) | Fiber optic sensing apparatus, system, and method for state of charge measurement in energy storage devices | |
| CN106352905A (zh) | 一种光纤光栅解调系统及解调仪 | |
| CN105758434A (zh) | 基于线阵InGaAs扫描FBG反射谱的传感解调方法 | |
| CN117387673B (zh) | 一种基于参考气室的光纤光栅解调方法及设备 | |
| CN106017533A (zh) | 一种快速调谐实时校准光纤光栅解调装置及工作方法 | |
| GB2163251A (en) | Infrared gas detector | |
| US5949533A (en) | Optical fiber distortion measurement device and method thereof | |
| CN105823497B (zh) | 一种基于信号自相关匹配的光纤光栅反射谱解调算法 | |
| CN116105776A (zh) | 光纤法珀传感器腔长高精度原位标定与高速实时解调的方法与装置 | |
| JP4607923B2 (ja) | ファイバブラッググレーティング素子反射光波長計測処理装置及び処理方法 | |
| Beam et al. | An electroactive fiber optic chip for spectroelectrochemical characterization of ultra-thin redox-active films | |
| CN111521206B (zh) | 光纤光栅边缘叠加滤波的波长解调方法 | |
| CN111537129B (zh) | 一种高速大动态范围光纤光栅解调系统及解调方法 | |
| CN113406017A (zh) | 高集成度表面等离子体共振传感器系统 | |
| CN108303196B (zh) | 全同光纤光栅数字化解调及区域温度监测方法及其系统 | |
| CN113916828B (zh) | 一种基于tdlas技术的湿度测量方法及系统 | |
| JP4080763B2 (ja) | ファイバブラッググレーティング素子反射光波長計測処理装置 | |
| CN113624267B (zh) | 一种基于边缘滤波的光纤光栅中心波长解调系统及解调仪 | |
| CN217358427U (zh) | 一种基于布拉格光纤光栅传感器的结构健康检测装置 | |
| CN113790743B (zh) | 用于光纤光栅传感系统的波长解调算法 | |
| CN111678594B (zh) | 一种激光功率测试仪响应线性的对数校准方法 | |
| CN211505199U (zh) | 紫外doas差分吸收光谱技术专用光纤光谱仪 | |
| CN117686008B (zh) | 一种基于图像处理的光纤布拉格光栅信号解调系统与方法 | |
| CN113310564B (zh) | 一种油井套管的振动参量和温度参量的测量系统及方法 | |
| CN117740728A (zh) | 基于l-m非线性拟合的tdlas气体测量方法及系统 |
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 |