光纤光栅位移传感器实时监测解调系统 下载: 1149次
1 引言
电学类传感器因有零点漂移、易被腐蚀和受电磁干扰严重等缺陷而无法较好地应用于某些恶劣环境中。光纤布拉格光栅(FBG)传感器具有无零点漂移、灵敏度高、抗环境干扰能力强和可长期监测等优势,弥补了电类传感器的不足,被广泛用于恶劣环境中温度、压力、应力、震动和位移等物理量的测量。例如光纤光栅位移传感器[1]是一种能够长期测量大型结构之间相对位移的传感器,在水坝、桥梁、工业与民用建筑等结构的裂缝或接缝开合度的测量中有较大应用前景。FBG的中心波长会随待测物理量的变化而变化,因此只要能检测到其中心波长的改变量,就能得到待测物理量的改变情况。FBG的中心波长检测技术有非平衡马赫-曾德尔(M-Z)干涉法、匹配光栅法、边缘滤波器法、可调谐法布里-珀罗(F-P)滤波器法等。
F-P滤波器以体积小、响应速度快、信噪比和分辨率高、稳定性好等特性[2],被广泛应用于FBG传感器的波长解调方案中[3],其基本原理是F-P滤波器的驱动电压和其透射峰的波长基本呈线性变化,当透射光强达到最大时,根据此时的驱动电压即可反推出入射光波长,但由于其压电陶瓷(PZT)元件具有非线性和迟滞性,并存在一定的温度漂移[4-5],因此其驱动电压和透射波长并不能用简单严格的线性关系表达,需要实时校准波长和电压的关系。常用的方法有3种:1) 引入1个FBG作为波长参考[6],其缺点是提供的校准波长参考点太少;2) 采用气体吸收谱线参考法[7],可以提供多个吸收谱线的波长参考点,但整个吸收谱线覆盖光波长范围较窄,谱线并非均匀分布,且各个峰的幅值强度差别较大,波长校正的算法复杂[7];3) 引入1个固定的参考F-P标准具[8-9],提供间隔均匀、强度相同,且覆盖光波长范围宽的多条参考波长线,但是谱线容易受到温度漂移的影响。
针对F-P滤波器压电陶瓷的驱动电压和透射波长的非严格线性关系,以及温度漂移造成的解调精度差等问题,基于可调谐F-P滤波器解调技术的原理、方法及其关键技术,同时综合考虑成本等多方面的因素,采用引入放置在恒温箱的高精度参考FBG,标定F-P标准具的某个峰的波长,根据F-P标准具的自由光谱特性可得到其余参考峰值波长,从而可进一步得到光纤传感器的波长偏移量,同时搭建了光纤光栅位移传感器实时监测解调系统平台,在数据处理中滤波算法采用改进的小波阈值去噪算法[10-11],针对原始阈值函数和阈值的不足,使用新的阈值函数和阈值,改善了原始信号,比较多种寻峰算法(如质心法、倒数法、多项式拟合法和半峰法)的优劣[12-14],最后采用指数修正高斯拟合算法[15]对提取的峰值信号进行二次过滤处理,去除假峰,修正部分峰的位置等,较好地提高了波长解调的准确度。
2 数据采集处理
2.1 光纤光栅位移传感器解调系统
根据宽带光源的光谱范围1525~1565 nm,F-P滤波器的扫描电压范围必须使其投射波长范围涵盖该宽带光源的范围。通过实验寻找到2个极值电压9.5 V(投射波长为1565 nm)和18.5 V(投射波长为1525 nm),所以锯齿波的扫描电压只要覆盖9.0~18.0 V即可。最终采用高性能的16位模数(AD)转换芯片AD6541进行模数转换,其转换建立时间为1 μs,满足系统要求。系统框架如
整个锯齿波产生过程为数字信号处理器(型号为DSPF2812)通过串行外接接口(SPI接)每隔5 μs送1个自加1的AD数值给数模转换(DAC)模块(型号为AD6541),通过数模(DA)转换产生对应的0~5 V线性电平,再通过模拟放大电路(1个比例运算放大器加1个加法器)将0~5 V电平搬迁到8~20 V,从而满足系统需求。理论上对应的F-P滤波器的扫描精度可以达到0.7 nm。
2.2 小波阈值去噪
光电探测器转换FBG反射谱信号存在大量噪声信号,而小波变换因具有低熵性、多分辨率、去相关性、选基灵活性等特性,在去噪方面得到了广泛应用,优于常见的均值滤波和中值滤波等算法。
设s(t)为采集到的原始信号,n(t)为服从
N(0,σ2)的高斯白噪声(σ为标准差),将高斯白噪声加到原始信号,得到含有噪声的信号f(t),可表示为f(t)=s(t)+n(t),对f(t)作离散小波变换
式中ωs(j,k)和ωn(j,k)分别为原始信号和噪声信号在第j层上的小波系数,J表示小波变换最大分解层,K表示信号长度。
由(1)式可知,原始信号能量大多表现在有限的小波系数上,且幅值大;噪声能量主要体现在全部小波域内,且幅值小。所以选择1个恰当的阈值,对小波系数作量化处理,再对量化后的系数进行重构,即可得到去除噪声的信号,流程如
常见的传统阈值函数包括硬阈值函数和软阈值函数,分别为
(2)式为硬阈值函数,(3)式为软阈值函数,硬阈值函数在±γ处具有不连续性,导致重构信号在不连续点附近震荡;软阈值函数虽有较好的连续性,但是待估计的小波系数与带噪的小波系数存在恒定偏差,也会给重构信号带来误差。
因此,在文献[
9-10]的研究基础上,有必要重新构造1个新的阈值函数,其特点是在阈值点处连续,且当ωj,k增大时,
由(4)式可知,当
阈值的选取常采用γ=σ
通过以上分析,使用新的阈值函数和阈值,同时采用Sym8小波基,分解层次到第5层时数据源的去噪平滑效果达到最好,如
图 4. 小波阈值去噪效果。(a)原始信号;(b)小波阈值去噪处理信号;(c)某个光栅通道某个峰值点未处理前的局部放大;(d)峰值点经过去噪后的放大
Fig. 4. Effects of wavelet threshold denoising. (a) Original signal; (b) wavelet threshold denosing signal; (c) partial magnification of some peak point in a grating channel for original signal; (d) partial magnification of the peak point for denoising signal
2.3 高斯拟合寻峰算法
对比文献[ 11-13]中质心法、倒数法、多项式拟合、半峰法和高斯拟合等多种算法的分析结果可知,高斯拟合效果最好,但需要进行大量的数据运算。结合本系统并非大规模的分布式系统,普通的计算机也满足数据的大量运算,因此采用高斯拟合作为寻峰算法的基本原理。
高斯曲线模型与传感采样得到的FBG 反射光谱信号整体相似,可表示为
式中G0为反射谱光功率的峰值,λs为反射谱光功率达到峰值时所对应的波长值,Δλ为反射谱的3 dB带宽值。将(5)式两边取对数,即可得到中心波长值点。这种方法具有较好的抗噪性能,是一种相对稳定的拟合办法。
本系统首先对源数据进行小波阈值平滑去噪,然后采用高斯拟合算法进行峰值粗提取,得到的峰值序列如
图 5. 理想峰型和双峰型。(a)理想峰型;(b)双峰型
Fig. 5. Ideal peak shape and double peak shape. (a) Ideal peak shape; (b) double peak shape
图 6. 无效峰型和左偏峰型。(a)无效峰型;(b)左偏峰型
Fig. 6. Invalid peak shape and left peakshape. (a) Invalid peak shape; (d) left peak shape
图 7. 高斯模型示意图。(a)标准高斯曲线;(b)非对称高斯曲线
Fig. 7. Schematic diagrams of Gaussian models. (a) Standard Gaussian curve; (b) asymmetric Gaussian curve
称为左偏峰值波形,其具有明显的非对称特征,与标准的高斯模型相比,峰顶位置偏向左边,如
最终通过高斯拟合寻峰算法初步得到峰值序列点,然后再对峰值序列点进行二次去除无效峰,修正部分峰值点位置,可以得到比较精确的峰值位置,效果如
经过该寻峰处理,获得峰值点的序号及位置,对
图 8. 高斯拟合寻峰算法效果。(a)整体经过高斯拟合和二次处理以后得到的峰值点位置;(b)其中某个峰的局部放大图
Fig. 8. Effects of Gaussian fitting algorithm. (a) Total peak points after Gaussian fitting and the second processing; (b) partial magnification of some peak point
于FBG传感通道的信号,其峰值点可表示为{p0,p1,…,pN},其中N表示光栅通道中光纤光栅的数目,p0表示参考光栅的位置信息,并且参考光栅的波长是已知的,记为λp0。F-P标准具通道中的峰值点信息记为{d0,d1,…,dM},其中M表示光栅通道中峰值点数量。通过对比分析光栅通道和F-P标准具通道的峰值点位置,假设参考光栅峰值点位于F-P标准具的某2个峰dk和dk+1之间,这2个峰的波长记为λdk和λdk+1。由于F-P标准具的峰值点波长是未知的,所以需要利用参考光栅对F-P标准具进行校准。在F-P标准具的光谱波形中,每个峰值点之间可近似看作波长随采样点线性变化,由此可得
F-P标准具中2个峰之间间隔1个自由光谱范围,对标准具可用频率间隔来表示自由光谱范围,即
频率f与波长λ之间满足λ=c/f,其中c为真空中的光速。联立(6)式、(7)式可得F-P标准具中序号为k的峰值点波长λdk。再根据(7)式可得F-P标准具任意1个峰值点波长:
对光栅通道中待测光栅pi,波长记为λpi,假设其位于F-P标准具中dl和dl+1、波长为λdl和λdl+1的2个峰值之间,同样根据波长近似线性变化的关系可得
从(8)式、(9)式即可得到光栅通道其他待测光栅的波长,从而根据其波长改变量和物理量灵敏度的对应关系得到其传感物理量的变化量。
3 位移传感器动态显示实验结果
实际搭建的光纤光栅位移传感器的硬件解调系统如
表 1. 测量结果
Table 1. Measurement results
|
光纤光栅位移传感器的最大量程为30 mm,选择0~25 mm量程,测量精度为0.01 mm的螺旋测微仪作为对比数据。位移传感器厂商提供的位移的具体计算公式为
式中λ1为波长较短的光栅(记为光栅1)的动态测量值,λ2为波长较长的光栅的动态测量值,λ10是光栅1的静态(0位移时)值,λ20是光栅2的静态(0位移时)值,η为传感器的灵敏度,本文为227 pm/mm,是采用螺旋测微仪测得的作为对比的数据,该传感器测量所得结果与螺旋测微仪测量结果的误差值E=X-
实验中缓慢将位移传感器从初始状态拉大到最大量程,选取其中部分数据点与螺旋测微仪所测得数据作对比得到
4 结论
本文采用F-P标准具和放置在恒温箱的高精度参考光栅校准F-P滤波器,设计搭建了测量光纤布拉格光栅位移传感器的波长解调系统,原始数据序列先通过改进的小波阈值去噪算法(使用新的阈值函数和阈值 )进行平滑滤波处理,改善原始信号,然后采用高斯拟合算法的核心思想得到粗峰值序列,再进行二次处理,去掉无效峰,修正部分峰的位置,保证提取的峰值序列可靠准确,最后对实际的光纤光栅位移传感器进行实验,使用螺旋测微仪进行对比,二者测量的绝对误差相对稳定在0.25 mm左右,在工业与民用建筑等结构上裂缝或接缝开合度的长期监测中有一定的应用前景。
[1] 田晓丹, 张会新, 刘文怡, 等. 温度补偿型光纤光栅位移传感器[J]. 电子器件, 2016, 39(3): 576-580.
[2] Ma Y, Wang C, Yang Y, et al. High resolution and wide scale fiber Bragg grating sensor[J]. Optics & Laser Technology, 2013, 50(2): 107-111.
[3] 杨刚, 许国良, 涂郭结, 等. 基于频谱分区的高精度光纤光栅波长解调系统[J]. 中国激光, 2015, 42(4): 0405001.
[5] 盛文娟, 张慧, 杨宁, 等. 基于光纤光栅传感系统的光纤法布里-珀罗滤波器热致非线性研究[J]. 激光与光电子学进展, 2017, 54(4): 042301.
[6] 王鹏, 赵虹, 刘杰, 等. 基于可调谐F-P滤波器的FBG波长解调系统的动态实时校准方法[J]. 光学学报, 2015, 35(8): 0806006.
[7] 汪金辉, 许雪梅, 丁家峰, 等. 基于法布里-珀罗标准具和多光栅校准的光纤布喇格光栅波长解调系统[J]. 光子学报, 2016, 45(6): 0606003.
[8] 梁霄, 刘铁根, 刘琨, 等. 一种可调谐光滤波器非线性实时标定方法研究[J]. 中国激光, 2010, 37(6): 1445-1449.
[9] 江俊峰, 何盼, 刘铁根, 等. 基于复合波长参考的温度稳定光纤光栅传感解调研究[J]. 光学学报, 2015, 35(10): 1006005.
[10] 江虹, 苏阳, 李进, 等. 一种小波阈值去噪在FBG解调中的改进方法[J]. 光通信研究, 2016( 2): 40- 43.
JiangH, SuY, LiJ, et al. An improved method of wavelet threshold de-noising applied to FBG demodulation[J]. Study on Optical Communications, 2016( 2): 40- 43.
[11] 曹京京, 胡辽林, 赵瑞. 一种改进小波阈值函数的光纤光栅传感信号去噪方法[J]. 传感技术学报, 2015, 28(4): 521-525.
[12] 陈志军, 白剑, 吴祖堂, 等. 光纤布喇格光栅反射谱寻峰算法优化及比较[J]. 光子学报, 2017, 44(11): 1112001.
[13] 胡正文, 庞成鑫, 程冯宇, 等. Gaussian-LM算法在光纤布拉格光栅反射光谱寻峰中的应用[J]. 激光与光电子学进展, 2017, 54(1): 013001.
[14] 胡颖, 莫文琴, 柴大虎. 基于F-P滤波器的光纤光栅解调中寻峰算法分析[J]. 光通信技术, 2014, 38(12): 17-20.
[15] 陈勇, 杨雪, 刘焕淋, 等. 指数修正高斯拟合寻峰算法处理FBG传感信号[J]. 光谱学与光谱分析, 2016, 36(5): 1526-1531.
尹建璟, 许雪梅, 丁一鹏. 光纤光栅位移传感器实时监测解调系统[J]. 激光与光电子学进展, 2018, 55(1): 010601. Yin Jianjing, Xu Xuemei, Ding Yipeng. Real-Time Monitoring and Demodulation System of Fiber Bragg Grating Displacement Sensor[J]. Laser & Optoelectronics Progress, 2018, 55(1): 010601.