基于相似匹配方法的光纤布里渊频移提取准确性影响因素分析 下载: 1159次
1 引言
布里渊分布式传感技术能够实现对光纤不同位置处温度和应变信息的测量,具有测量精度高、传感距离长和空间分辨率高等优势,已被广泛应用于光缆、光纤复合海底电缆和大型建筑物等健康状况的在线监测。其主要利用布里渊频移(BFS)与光纤所受的温度和应变之间的线性关系得到相应的温度和应变信息,因此温度和应变的测量精度与布里渊频移的准确性有着密切的关系。
目前,关于布里渊频移提取的主流算法有基于最小二乘拟合的算法和基于模式匹配的算法。前者主要涉及谱模型、优化算法及初值选择方面的研究[1-4]。肖尚辉等[1]研究了基于Levenberg-Marquardt(L-M)算法的布里渊散射谱最优化拟合问题,推导了按照洛伦兹型和高斯型线性加权组合的pseudo-Voigt模型;张有迪等[2]利用小波去噪结合L-M算法调节权值后向传输网络,并对布里渊散射谱进行特征提取,提高了系统在低信噪比情况下的拟合精度;郝蕴琦等[3]研究了洛伦兹曲线拟合自发布里渊散射分段傅里叶变换后每个空间分辨率范围内频域信号的准确性;李晓娟等[4]研究了多波长布里渊传感系统中的散射谱特性,利用pseudo-Voigt函数和L-M优化算法对布里渊散射叠加谱进行特征提取。该算法具有较高的准确性[5],但计算速度偏慢。可以采用优化光路设计[6-7]、斜坡法[8-9]或基于线性最小二乘的拟合算法[10-11]加快谱测量和布里渊频移提取的速度,但线性最小二乘拟合算法对扫频范围关于布里渊频移的对称性要求较高。与以上算法不同,基于模式匹配的算法采用相关性指标来确定布里渊频移,主要包括互相关算法[12-13]和相似匹配方法[14],有望抑制扫频范围偏差对布里渊频移提取准确性的影响。Farahani等[12]提出了基于互相关卷积的布里渊频移提取方法,对服从洛伦兹分布的有噪信号与理想洛伦兹信号进行卷积运算,得到在峰值附近几乎呈理想洛伦兹分布的信号,改善了信噪比,但其准确度受到扫频数据服从洛伦兹分布的程度的限制。为了进一步提高准确性,尚秋峰等[13]提出采用高阶矩质心计算来改进基于互相关卷积的布里渊频移提取方法,该算法避免了迭代过程,实时性很好。现有方法多从模型角度出发,期望通过采用更符合布里渊散射谱数据分布的模型使模拟结果更好逼近实测布里渊谱,进而能更加准确地获得布里渊频移。Wang等[14]另辟蹊径,不是从模型角度出发而是根据谱相似性来确定布里渊频移的改变量,进而提出了最大相似匹配方法(SMM),比较受到温度和应变影响前后的布里渊频谱,当相似性达到最高时提取得到对应的频移改变量,该算法的突出优点是频移提取时无需确定谱模型,在布里渊谱较为复杂的情况下也具有很好的适应性。布里渊线宽(简称线宽)、信噪比(SNR)、扫频范围、扫频间隔等参数均对该算法的准确性有影响,但目前没有给出以上参数对该算法影响的报道,这影响了该方法的推广应用。
为了解决该问题,本文在基于相似匹配方法的布里渊频移提取算法的基础上系统研究了线宽、信噪比、扫频间隔、扫频范围对布里渊频移提取准确性的影响规律,本文结果可为分布式光纤传感中布里渊频移的准确提取提供参考。
2 基于相似匹配方法的布里渊频移提取算法
2.1 布里渊谱及参数
考虑到光纤中声波的振幅按exp(t/τ)衰减,其中τ为声子寿命,当脉冲宽度大于50 ns时实测布里渊散射谱近似满足洛伦兹模型,如
式中:g0为布里渊增益峰值;vB为布里渊频移;ΔvB为线宽;v为光纤入射光频率。
布里渊频移vB可表示为
式中:n为光纤的折射率;VA为光纤中的声速;λp为泵浦光波长。若取n=1.45、VA=5945 m/s,对于普通单模光纤,当入射光波长为1550 nm时,根据(2)式计算得到的布里渊频移约为11.12 GHz。
线宽ΔvB可表示为[15]
式中:c为光纤中的光速;Γ为声波阻尼参数;vp为泵浦光频率。对于普通单模光纤,布里渊谱自然线宽约为35 MHz[16],实际测量得到的线宽还会受到入射光脉宽和泵浦功率的影响[17]。
随着注入传感光纤的调制光脉冲越来越窄,布里渊谱变得越来越宽,其形状也逐渐由洛伦兹型向高斯型过渡,这时的谱可表示为洛伦兹型谱和高斯型谱的线性组合,并近似满足pseudo-Voigt模型[18],具体可表示为
式中:g01和g02分别表示洛伦兹模型和高斯模型的增益峰值。
当光纤支持传输多个声学模式时,每个声学模式对应形成一个布里渊峰,此时的布里渊谱等于每个声学模式对应的布里渊谱的总和[19]。在一些纤芯经过特殊处理的光纤中,如光子晶体光纤和大有效面积光纤,由于是每个布里渊谱的叠加,总的布里渊谱为展宽的多峰值谱[20-21]。针对多峰值谱的计算公式较为复杂,这也给基于模型的布里渊频移提取算法增加了难度。
2.2 相似匹配方法
为减小由布里渊谱的模型选择不合适引起的误差,相似匹配方法从布里渊谱形状的相似性来考虑[22],以相关系数作为相似性的评价标准,根据温度或应变改变前后布里渊谱的相似性来获得布里渊频移的改变量。该方法的具体实现为分别获得光纤处于自然状态时的散射谱和受到温度或应
变作用后的散射谱,文中分别称为参考谱(SRS)和探测谱(SDS),将参考谱的不同部分与探测谱比较,当获得的相关系数最大时,得到的频率差即为所求的布里渊频移改变量[14],其原理如
当光纤处于自然状态时,获得的参考谱可以表示为
式中:f0为参考谱的起始频率;M为参考谱的扫频点数;F0为参考谱的扫频间隔。
当光纤受到温度或应变作用时,获得的探测谱可以表示为
式中:f为探测谱的起始频率;N为探测谱的扫频点数;F为探测谱的扫频间隔。
从频率f0开始依次截取参考谱中扫频间隔为F的N个扫频点,并计算获得参考谱与探测谱的一个相关系数。从f0开始以F0递增,这样可以获得一组相关系数,其中第j个相关系数可以表示为
式中:
设j=J时相关系数达到最大,可以得到布里渊频移的改变量Δf:
2.3 算法有效性验证
为了验证基于相似匹配方法的布里渊频移提取算法的有效性,本节分别采用基于pseudo-Voigt模型和相似匹配方法的算法对实测谱进行布里渊频移计算。pseudo-Voigt模型能较好逼近布里渊谱,无论是谱拟合结果还是布里渊频移计算结果都比较可靠,因此以pseudo-Voigt模型算法的结果作为准确值去验证相似匹配方法的准确性。通常的温升实验是将一小段光纤置于恒温箱或恒温槽中,可用于验证的谱数量较少。为了提高验证结果的可靠性,将一整盘光纤置入恒温环境中测量不同叠加平均次数下光纤沿线的布里渊谱,以两种情况下布里渊频移改变量随光纤沿线的变化作为参考对象,比较pseudo-Voigt模型算法和相似匹配方法的差异。由于频移本身基本不会影响算法的准确性,从定量误差的角度来说分析的结果与温度变化的结果一致,但谱信号数量的增加使得验证更充分。
基于布里渊光时域反射计(AV6419,中电科仪器仪表有限公司)搭建了光纤布里渊谱测量系统,所用光纤为1 km左右的普通单模光纤(G657.A2)。扫频范围为10.52~10.92 GHz,扫频间隔为1 MHz,入射脉冲光波长为1550 nm,脉冲宽度为10 ns,叠加平均次数分别为214和218,采用pseudo-Voigt模型算法得到一个线宽内的信噪比分别为27.48 dB和32.74 dB。为了测量更多的谱,将采样分辨率设置为0.1 m。实验在室温下进行,但由于绕制光纤承受了应变,布里渊频移会随光纤位置变化。参考谱和探测谱的扫频范围分别选择为线宽的4倍和2倍,后续仿真分析设置相同。两种情况下典型的布里渊谱如
图 3. 实测布里渊谱及频移计算结果。(a)计算用的布里渊谱;(b)叠加平均次数为218时模型方法算得的光纤沿线布里渊频移;(c)叠加平均次数为214时两种方法布里渊频移计算结果的差距
Fig. 3. Measured Brillouin spectrum and calculated frequency shift results. (a) Brillouin spectra for calculation; (b) BFS along the fiber by the model method when the average superposition is 218; (c) difference in calculated BFS between the two methods when the average superposition is 214
3 影响因素分析
不同扫频参数下,实际布里渊谱测量时谱信号差异很大,非常有必要研究它们的取值对算法性能的影响。不失一般性,采用(1)式所示的洛伦兹模型产生单峰布里渊谱,考虑到实际单模光纤布里渊谱测量时相关参数的设置以及相似匹配方法的特性,设置参考谱和探测谱的线宽均为35 MHz。参考谱的信噪比为30 dB,扫频间隔为1 MHz,扫频范围为4倍的线宽;探测谱的信噪比为20 dB,扫频间隔为3 MHz,扫频范围为2倍的线宽。分析某一参数的影响时,只改变该参数的值,其他参数保持不变。实际扫频过程中,扫频点未必能与实际布里渊频移重合,这会明显影响算法的准确性,但现有研究容易忽略这一点,本文为了模拟真实情况考虑了这部分误差。另外,本文设置的参考谱的信噪比一般不低于探测谱的信噪比,原因是探测谱需要与参考谱进行匹配,一次实测的参考谱可以反复使用,增加叠加平均次数可提高参考谱的信噪比,从而既可以有效提高布里渊频移提取[23]的准确性又不会对后续测量的实时性造成影响。
3.1 线宽
在实际应用中,为了提高空间分辨率需要减小泵浦脉冲的宽度,这会引起线宽的展宽。由2.1节可知,普通单模光纤的线宽约为35 MHz,本文用布里渊光时域反射计对 G657.A2型普通单模光纤进行测量,实测结果表明当脉宽在10~200 ns范围内变化时线宽在40~100 MHz范围内变化,考虑到光纤受到温度或应变影响后线宽可能存在波动,故仿真时设置参考谱和探测谱的线宽为30~100 MHz,其他参数与第3节第1段一致。不同探测谱点数下,线宽对频移提取准确性的影响结果如
图 4. 不同探测谱点数时,线宽对频移提取准确性的影响
Fig. 4. Influence of linewidth on frequency shift extractionaccuracy under different sampling point number of SDS
由
式中:ΔvB为参考谱和探测谱的线宽,单位为MHz,p为布里渊频移误差与线宽间的比例系数,参考谱和探测谱的线宽同步变化;EBFS表示布里渊频移误差,单位为MHz。
表 1. 线宽对频移提取准确性影响的拟合结果
Table 1. Fitting results of influence of linewidth on frequency shift extraction accuracy
|
由
3.2 仅探测谱线宽改变
温度和应变可能会影响探测谱的线宽,导致其与参考谱线宽未必一致,需要研究其对测量准确性的影响。文献[
24]中的研究结果表明温度在-40~80 ℃范围内变化时线宽波动程度约为±15%(以20 ℃时的线宽作为基准值),常规应变变化范围内线宽波动程度更小,文献[
16]的结果也非常相似。文献[
17]研究认为泵浦功率的起伏也会对布里渊增益谱的线宽产生影响,分析数据发现功率变化范围内线宽波动明显小于±15%。因此,后续计算中线宽的波动程度取±15%。设置参考谱的线宽为35 MHz,探测谱的线宽变化范围为30~40 MHz,线宽变化时其扫频范围始终为2倍线宽,其他参数与第3节第1段一致,探测谱线宽与参考谱线宽不一致对算法的影响如
由
图 5. 不同探测谱信噪比下,探测谱线宽对算法准确性的影响
Fig. 5. Influence of linewidth of SDS on algorithm accuracy under different SNR of SDS
3.3 参考谱的信噪比
参考谱信噪比过低则产生的误差较大。分析其信噪比对频移提取准确性的影响规律时,为使结果覆盖范围更广从而提高分析的可靠性,设置参考谱的信噪比为20 dB~40 dB,探测谱的信噪比范围为0 dB~20 dB,其他参数与第3节第1段相同。参考谱信噪比对频移提取准确性的影响结果如
图 6. 不同探测谱信噪比下,参考谱信噪比对频移提取准确性的影响
Fig. 6. Influence of SNR of SRS on frequency shift extraction accuracy under different SNR of SDS
由
式中:RSN为参考谱的信噪比,单位为dB;EBFS表示布里渊频移误差,单位为MHz。按照(11)式拟合的结果如
表 2. 参考谱信噪比对频移提取准确性影响的拟合结果
Table 2. Fitting results of influence of SNR of SRS on frequency shift extraction accuracy
|
随着参考谱信噪比的提高,误差逐渐减小,同时误差的变化幅度也逐渐减小,当参考谱的信噪比达到30 dB时,误差基本保持不变。在扫频时,通常利用时域叠加平均技术来改善信噪比,一般叠加平均次数增加到原来的4倍时信噪比提高6 dB,但增加平均次数测量时间也会随之增加。因此,参考谱的信噪比选择30 dB可保持较高的频移提取准确性。
3.4 探测谱的信噪比
考虑到3.3节中参考谱信噪比达到30 dB可得到较高的准确性,将参考谱信噪比分别设置为30 dB、35 dB和40 dB,为了系统地分析探测谱信噪比对频移提取准确性的影响,仿真时设置其范围为0 dB~30 dB,其他参数与第3节第1段一致。参考(11)式对布里渊频移误差进行拟合,
图 7. 不同参考谱信噪比下,探测谱信噪比对频移提取准确性的影响
Fig. 7. Influence of SNR of SDS on frequency shift extraction accuracy under different SNR of SRS
由
表 3. 探测谱信噪比对频移提取准确性影响的拟合结果
Table 3. Fitting results of influence of SNR of SDS on frequency shift extraction accuracy
|
由
3.5 参考谱的扫频间隔
考虑到探测谱扫频间隔需为参考谱扫频间隔的整数倍,结合谱测量的实际情况,不失一般性,设置探测谱的扫频间隔为5 MHz,参考谱的扫频间隔为探测谱的1/10~1,其他参数与第3节第1段一致。参考谱扫频间隔对布里渊频移提取准确性的影响如
由
图 8. 不同探测谱信噪比下,参考谱扫频间隔对频移提取准确性的影响
Fig. 8. Influence of sweeping step of SRS on frequency shift extraction accuracy under different SNR of SDS
3.6 探测谱的扫频间隔
设置参考谱的扫频间隔为1 MHz,探测谱的扫频间隔分别为参考谱的1~10倍,其他参数与第3节第1段一致。探测谱扫频间隔对布里渊频移提取准确性的影响如
由
图 9. 不同探测谱信噪比下,探测谱扫频间隔对频移提取准确性的影响。(a)频移提取误差;(b)归一化频移提取误差
Fig. 9. Influence of sweeping step of SDS on frequency shift extraction accuracy under different SNR of SDS.(a) Error of BFS; (b) normalized error of BFS
3.7 参考谱的扫频范围
利用相似匹配方法提取布里渊频移时,参考谱的扫频范围应包含探测谱的扫频范围,由于探测谱的扫频范围设置为2倍线宽,所以将参考谱的扫频范围设置为3~15倍的线宽,其他参数与第3节第1段一致。注意,参考谱扫频范围变化过程中其扫频间隔一直保持不变。参考谱扫频范围对频移提取准确性的影响结果如
图 10. 不同探测谱信噪比下,参考谱扫频范围对频移提取准确性的影响。(a)频移提取误差;(b)算法运行时间
Fig. 10. Influence of sweeping range of SRS on frequency shift extraction accuracy under different SNR of SDS. (a) Error of BFS; (b) algorithm runtime
由
3.8 探测谱的扫频范围
设置参考谱扫频范围为4倍的线宽,探测谱扫频范围为0.5~3倍的线宽,其中探测谱扫频范围变化时保持扫频点数不变,其他参数与第3节第1段一致,探测谱扫频范围对布里渊频移提取准确性的影响如
由
图 11. 探测谱扫频范围对频移提取准确性的影响。(a)探测谱关于频移对称;(b)探测谱不完全对称
Fig. 11. Influence of sweeping range of SDS on frequency shift extraction accuracy. (a) SDS symmetrical about vB; (b) SDS not symmetrical about vB
4 结论
在验证相似匹配方法有效性的基础上,本文研究了线宽、信噪比、扫频间隔、扫频范围对基于该方法的布里渊频移提取准确性的影响规律。在保持扫频点数不变的情况下,当参考谱和探测谱的扫频范围与线宽比值固定时,频移误差与线宽成线性关系。排除线宽对频移误差影响后仅仅由探测谱与参考谱线宽不一致导致的布里渊频移误差较小。频移误差随参考谱和探测谱信噪比的增加成指数规律减小。随着扫频间隔的增加,频移误差有增加的趋势。扫频间隔不变时,随参考谱扫频范围增加误差有增加的趋势,在包含探测谱范围的基础上,参考谱可适当选择较小的扫频范围。扫频点数不变时,探测谱扫频范围为2倍线宽时频移误差最小。参考谱和探测谱的扫频范围可分别选择为4倍和2倍线宽。以上参数可根据实际要求适当调整。研究结果可为布里渊频移的准确提取提供参考。
[1] 肖尚辉, 李立. 一种新的光纤布里渊传感散射谱拟合方法[J]. 光学技术, 2009, 35(6): 897-900, 904.
Xiao S H, Li L. New fitting method for Brillouin-based scattering spectrum of fibre-optic distributed sensing systems[J]. Optical Technique, 2009, 35(6): 897-900, 904.
[2] 张有迪, 李嘉琪, 孟钏楠, 等. 布里渊散射谱拟合的混合优化算法[J]. 强激光与粒子束, 2015, 27(9): 091013.
[3] 郝蕴琦, 叶青, 蔡海文. 基于BOTDR布里渊散射的洛伦兹拟合准确度分析[J]. 激光杂志, 2016, 37(5): 15-17.
Hao Y Q, Ye Q, Cai H W. Lorentz fitting analysis of Brillouin scattering in BOTDR sensing system[J]. Laser Journal, 2016, 37(5): 15-17.
[4] 李晓娟, 李永倩, 胡智奇, 等. 多波长BOTDR系统中布里渊散射谱的特征提取[J]. 光电子·激光, 2016, 27(12): 1286-1292.
Li X J, Li Y Q, Hu Z Q, et al. Features extraction of Brillouin scattering spectra in multi-wavelength BOTDR system[J]. Journal of Optoelectronics·Laser, 2016, 27(12): 1286-1292.
[5] Xu Z N, Zhao L J. Estimation of error in Brillouin frequency shift in distributed fiber sensor[J]. IEEE Sensors Journal, 2020, 20(4): 1829-1837.
[6] Zhou D W, Dong Y K, Wang B Z, et al. Single-shot BOTDA based on an optical chirp chain probe wave for distributed ultrafast measurement[J]. Light: Science & Applications, 2018, 7: 32.
[7] Wang B Z, Hua Z J, Pang C, et al. Fast Brillouin optical time-domain reflectometry based on the frequency-agile technique[J]. Journal of Lightwave Technology, 2020, 38(4): 946-952.
[8] Bernini R, Minardo A, Zeni L. Dynamic strain measurement in optical fibers by stimulated Brillouin scattering[J]. Optics Letters, 2009, 34(17): 2613-2615.
[9] 徐志钮, 胡宇航, 赵丽娟, 等. 基于单斜坡法的光电复合海缆温度、应变快速测量方法[J]. 电力自动化设备, 2020, 40(5): 202-208.
Xu Z N, Hu Y H, Zhao L J, et al. A rapid temperature and strain measurement method for optic-electric composite submarine cable based on slope-assisted technique[J]. Electric Power Automation Equipment, 2020, 40(5): 202-208.
[10] Xu Z N, Zhao L J. Accurate and ultra-fast estimation of Brillouin frequency shift for distributed fiber sensors[J]. Sensors and Actuators A: Physical, 2020, 303: 111822.
[11] 徐志钮, 胡宇航, 赵丽娟, 等. 基于改进二次多项式拟合的布里渊频移快速高精度提取算法[J]. 光谱学与光谱分析, 2020, 40(3): 842-848.
[12] Farahani M A, Castillo-Guerra E, Colpitts B G. Accurate estimation of Brillouin frequency shift in Brillouin optical time domain analysis sensors using cross correlation[J]. Optics Letters, 2011, 36(21): 4275-4277.
[13] 尚秋峰, 胡雨婷, 刘薇. 基于互相关卷积与高阶矩质心计算的布里渊散射谱特征提取[J]. 中国激光, 2017, 44(11): 1106011.
[14] Wang F, Zhan W W, Lu Y G, et al. Determining the change of Brillouin frequency shift by using the similarity matching method[J]. Journal of Lightwave Technology, 2015, 33(19): 4101-4108.
[15] 安琪. 矢量受激布里渊散射特性及其应用的理论与实验研究[D]. 北京: 华北电力大学, 2017: 21- 22.
AnQ. Theoretical and experimental research on the characteristics and applications of vector stimulated Brillouin scattering[D]. Beijing: North China Electric Power University, 2017: 21- 22.
[16] Bao X Y, Smith J, Brown A W. Temperature and strain measurements using the power, line-width, shape, and frequency shift of the Brillouin loss spectrum[J]. Proceedings of SPIE, 2002, 4920: 311-322.
[18] Zhang Y J, Yu C J, Fu X H, et al. An improved Newton algorithm based on finite element analysis for extracting the Brillouin scattering spectrum features[J]. Measurement, 2014, 51: 310-314.
[19] 潘宇航, 路元刚, 彭楗钦, 等. 光子晶体光纤的布里渊增益谱特性[J]. 光学学报, 2019, 39(6): 0619001.
[20] 梁浩, 张旭苹, 李新华, 等. 布里渊背向散射光谱数据拟合算法设计与实现[J]. 光子学报, 2009, 38(4): 875-879.
[22] Yoo J C, Han T H. Fast normalized cross-correlation[J]. Circuits, Systems and Signal Processing, 2009, 28(6): 819-843.
[23] 尚秋峰, 胡雨婷, 刘薇. BOTDA传感系统的布里渊频移提取方法研究进展[J]. 半导体光电, 2017, 38(5): 633-638, 669.
[24] Nikles M, Thevenaz L, Robert P A. Brillouin gain spectrum characterization in single-mode optical fibers[J]. Journal of Lightwave Technology, 1997, 15(10): 1842-1851.
Article Outline
赵丽娟, 王贺晴, 徐志钮, 车小霞. 基于相似匹配方法的光纤布里渊频移提取准确性影响因素分析[J]. 中国激光, 2020, 47(5): 0506003. Lijuan Zhao, Heqing Wang, Zhiniu Xu, Xiaoxia Che. Analysis of Factors Affecting Accuracy of Brillouin Frequency Shift Extraction Based on Similarity Matching[J]. Chinese Journal of Lasers, 2020, 47(5): 0506003.