CN111121836A - 基于改进二次多项式拟合的布里渊频移快速准确提取方法 - Google Patents

Abstract

一种基于改进二次多项式拟合的布里渊频移快速准确提取方法，所述方法是在分布式光纤传感器所输出的测量信号中，以布里渊频移为中心截取1个布里渊增益谱的线宽Δv_B的谱信号，对该信号进行二次多项式处理方法的拟合，得到布里渊谱：g_B(v)＝av²+bv+c，则布里渊频移为

Description

基于改进二次多项式拟合的布里渊频移快速准确提取方法

技术领域

本发明涉及一种改进的布里渊频移快速提取方法，属于测量技术领域。

背景技术

分布式光纤传感器除具有一般光纤传感器抗电磁干扰、耐腐蚀以及电绝缘性好等优点外，还具有只需一次测量即可获取沿整个光纤被测场分布信息等独特优点，另外其测量精度高、定位准确、传感距离可达上百公里。因此，应用领域非常广泛。其中，基于光纤布里渊散射的分布式传感技术在温度和应变测量上所达到的测量精度、测量范围以及空间分辨率高于其它类型的分布式光纤传感技术，因此引起了国内外的广泛关注和研究。利用该技术对油气管道、大型水利水电工程结构和电力线缆等的温度和应变进行在线监测，可实现故障隐患及故障点的快速、准确定位。

在与温度和应变有关的谱特征量中，布里渊频移与温度和应变呈线性关系且最为稳定，目前基于布里渊散射的光纤分布式传感测量时绝大多数需要提取布里渊频移信息。布里渊频移的准确、快速测量对基于布里渊散射的光纤分布式传感***非常关键。整体上看布里渊频移的提取方法以拟合基于的方法为主，可以分为模型基于的处理方法和非模型基于的处理方法。由于能充分利用布里渊谱的波形数据来提取布里渊频移，整体上看模型基于的处理方法具有更高的准确性。当入射矩形脉冲宽度明显大于10ns(大于50ns)时光纤中实测布里渊谱近似满足洛伦兹函数，目前洛伦兹模型基于的拟合处理方法应用最为广泛。当然，随着脉冲宽度的增加空间分辨率会随之减少，为了提高空间分辨率常使用较小的脉冲宽度，当其值小于10ns时布里渊谱近似满足高斯分布，此时更应该采用高斯模型基于的处理方法。显然，这样根据脉冲宽度来选择处理方法的方式给布里渊频移的提取增加了一个干扰因素，迫切需要一个通用模型基于的处理方法。可以认为在矩形脉冲的整个脉宽范围内布里渊谱更满足Voigt函数，因此Voigt模型基于的处理方法具有更高的准确性。但Voigt函数并非代数表达式，只有数值解，虽然基于它的处理方法准确性最高但计算速度过慢，不是特别适合于长传感范围、高空间分辨率的场合。为了提高计算速度，采用洛伦兹和高斯模型线性组合的伪Voigt模型来近似布里渊谱，具有很好的表征效果。为了进一步提高计算速度和收敛能力，人们开展了大量工作。为了提高布里渊频移提取的准确性，有的学者采用有限元分析来改进Newton处理方法和Levenberg-Marquardt处理方法，然后采用不同信噪比和线宽的布里渊谱的计算结果验证了改进处理方法在精度提高方面的有效性。针对变量初值选择比较困难可能导致处理方法发散的问题，有的学者分别采用了布谷鸟牛顿搜索和粒子群优化来获得更加靠近最优解的初值，然后采用Levenberg-Marquardt处理方法来优化变量，能确保收敛。虽然以上处理方法提高了准确性、降低处理方法发散的概率，但是该模型所对应的是非线性目标函数，需要采用迭代优化处理方法才能得到最优解，通常情况下计算速度仍然偏慢。这在传感距离较长、空间分辨率较高时比较突出。该问题仍有待于进一步研究。

发明内容

本发明的目的在于针对现有技术之弊端，提供一种基于改进二次多项式拟合的布里渊频移快速准确提取方法，在保证布里渊频移提取精度的同时，提高处理速度。

本发明所述问题是以下述技术方案解决的：

一种基于改进二次多项式拟合的布里渊频移快速准确提取方法，所述方法针对分布式光纤传感器所输出的测量信号，以布里渊频移为中心截取1个布里渊增益谱的线宽Δv_B的谱信号，对该信号进行二次多项式处理方法的拟合，得到布里渊谱：

g_B(v)＝av²+bv+c

其中a，b，c为二次多项式系数，其数值随布里渊谱的变化而不同，则布里渊频移为

上述基于改进二次多项式拟合的布里渊频移快速准确提取方法，在以布里渊频移为中心截取1个布里渊增益谱的线宽Δv_B的谱信号时，所述布里渊频移由布里渊谱峰值对应频率代替。

本发明选取峰值附近对称的一个线宽的布里渊谱、然后采用二次多项式来逼近布里渊谱进而提取布里渊频移，该方法不仅能大幅提高计算速度,计算时间仅为现有典型方法(洛伦兹、高斯、伪Voigt和Voigt模型基于的拟合方法)的1.15％、1.80％、1.51％和0.51％。而且计算准确性与典型的基于洛伦兹、高斯、伪Voigt和Voigt模型基于的拟合方法的处理方法相当，折算为温度的差距为温度误差仅分别为0.35℃、0.21℃、0.21℃和0.23℃。采用数值仿真及实测布里渊谱验证了提出处理方法的有效性。本发明为实现光纤分布式传感的快速测量提供了很好的支持。

附图说明

下面结合附图对本发明作进一步详述。

图1是基于伪Voigt和二次多项式处理方法的拟合结果；

图2是不同处理方法布里渊频移的提取结果和计算时间，实测谱，其中图2(a)为布里渊频移，图2(b)为计算时间；

图3是不同处理方法布里渊谱的拟合结果，实测谱；

图4是不同处理方法布里渊频移的提取结果，仿真谱；

图5是不同处理方法布里渊谱的拟合结果，仿真谱；

图6是布里渊频移误差与扫频范围的关系；

图7是布里渊频移误差与扫频点数的关系；

图8是布里渊频移误差与信噪比的关系；

图9是布里渊频移误差与线宽的关系；

图10是布里渊频移误差与扫频范围偏差的关系；

图11是本发明与经典处理方法的布里渊频移提取结果和计算时间，实测谱，其中图11(a)为布里渊频移，图11(b)为计算时间；

图12是本发明与经典处理方法的布里渊频移提取结果，仿真谱。

文中所用符号为：v为频率值；v_B是布里渊频移；Δv_BL和Δv_BG分别为Lorentzian型和Gaussian型布里渊谱的线宽；Δv_B为设布里渊增益谱的线宽；g₀₁和g₀₂分别为布里渊散射谱的Lorentzian和Gaussian峰值增益；E_vB为布里渊频移误差幅值的均值；g_B(v)为布里渊散射谱的增益；v为频率；a，b，c为二次多项式系数。

具体实施方式

针对实测和仿真布里渊谱的计算发现，二次多项式拟合处理方法的计算速度较现有经典模型基于处理方法有显著提高，但误差偏大。此后***研究了扫频范围、扫频点数、信噪比、线宽和扫频范围偏差对基于二次多项式的布里渊频移提取准确性的影响。根据研究结果，本发明提出了改进的二次多项式拟合处理方法，改进的处理方法不仅能大幅提高计算速度且计算准确性与经典处理方法相当。采用数值产生及实测布里渊谱验证了提出处理方法的有效性。

1处理方法原理

入射脉冲光为矩形波时会发生布里渊散射，布里渊谱近似满足如下的Voigt模型，它是洛伦兹函数与高斯函数的卷积形式

式中，v为频率值，单位为GHz；v_B是布里渊频移，单位为GHz；Δv_BL和Δv_BG分别为Lorentzian型和Gaussian型布里渊谱的线宽，单位为GHz。可以采用计算更为快速的伪Voigt模型来表征布里渊谱，设布里渊增益谱的线宽为Δv_B；g₀₁和g₀₂分别为布里渊散射谱的Lorentzian和Gaussian峰值增益，布里渊谱表示如下

设

v＝x+v_B (3)

将式(3)代入式(2)同时根据泰勒级数展开则有

式中O(x)为x的高阶无穷小。

如果x较小，如x≤Δv_B/2，忽略高阶无穷小后根据式(4)则有

根据式(3)和(5)则有

因此，在一个线宽范围内布里渊增益谱近似满足二次多项式关系。可以采用二次多项式来逼近布里渊谱，即

g_B(v)＝av²+bv+c (7)

采用二次多项式拟合后布里渊频移计算公式如下

注意，如果求得布里渊频移超出了谱的扫频范围，此时布里渊频移误差可能会超乎正常的大。则处理方法实现时选择扫频范围内满足均匀分布的随机值作为算得的布里渊频移，这样可以减少测量误差且使测量结果显得更正常。

由图1可知，伪Voigt函数能很好逼近Voigt函数，而二次多项式虽然与真实的Voigt函数从表达式上差别不小，但也有希望逼近对应的谱线。

2处理方法比较

2.1实测谱

本发明基于中电科仪器仪表有限公司生产的AV6419型光时域反射计(BOTDR,Brillouin Optical Time Domain Reflectometry)搭建了光纤布里渊谱测量***，选择了1km左右的SM 9/125μm光纤。扫频范围为10.52～10.92GHz，扫频间隔为1MHz，入射脉冲光波长为1550nm，脉冲宽度为10ns，采样分辨率为10m，叠加平均次数为2¹⁸。实验在室温下进行，但由于绕制光纤承受了应变，因此沿线布里渊频移并非恒定值。采用基于洛伦兹、高斯、伪Voigt和Voigt模型处理方法以及二次多项式拟合处理方法计算，布里渊频移的提取结果和计算时间如图2所示。选择一个典型谱5种处理方法的拟合结果如图3所示。

由图2(a)可知，基于洛伦兹、高斯、伪Voigt和Voigt模型处理方法的布里渊频移提取结果几乎一致，这与图3中基于高斯、伪Voigt和Voigt模型处理方法的拟合结果能较好逼近实测谱这一结果吻合，这些结果也基本验证了基于谱模型的布里渊频移提取的可靠性。然而，在大部分区域基于二次多项式处理方法的布里渊提取结果要小于其他处理方法的计算结果，它与以上4种处理方法差距的均值分别为0.70、0.74、0.73和0.73MHz，差距的最大值分别为2.47、2.60、2.57和2.60MHz。由于以上分析认为基于谱模型处理方法的计算结果比较可靠，故二次多项式拟合处理方法存在明显误差，这与二次多项式拟合结果与实测谱差距较大比较吻合(图3)。当然这一结论会在后续分析中得到进一步验证。由图2(b)可知，基于Voigt模型处理方法计算速度最慢，然后依次是基于洛伦兹、伪Voigt、高斯模型的处理方法，最快的是二次多项式拟合处理方法，图2(b)中5种处理方法对应的平均计算时间分别为58.97ms、37.69ms、44.72ms、131.48ms和6.76ms，即二次多项式拟合处理方法的计算时间仅为前4种处理方法的1.15％、1.80％、1.51％和0.51％。也就是说二次多项式拟合处理方法具有非常快的计算速度，但计算误差可能偏大，因此需要改正该处理方法，这也就是本发明技术方案的核心。

注意，本发明不是根据一定温度和应变下光纤布里渊谱数据来分析处理方法准确性。这是因为这种方式通常默认处理方法计算结果为准确值，而对于目前来说，二次多项式拟合处理方法的准确性本身就需要核对，故不宜采用该方式。本发明选择信噪比较高信号的原因是过大的噪声可能会掩盖布里渊频移随光纤位置的变化规律。

2.2数值产生谱

2.1节已有实测谱的计算，本节采用数值产生谱进一步计算的原因是：1)实测谱虽然可靠但布里渊频移不能足够准确获得；2)后续涉及大量数值产生谱的分析，需要验证数值产生谱分析结果的可靠性。布里渊谱根据式(1)产生，其中的A、v_B、Δv_BL和Δv_BG由Voigt模型处理方法针对2.1节实测谱计算获得。信噪比与实测谱一致，约为33dB。布里渊频移的提取结果如图4所示。5种处理方法针对一个与图3实测谱对应典型谱的拟合结果如图5所示。

由图4可知，5种处理方法算得布里渊频移及随光纤位置的变化规律与实测谱非常接近。仍然是基于洛伦兹、高斯、伪Voigt和Voigt模型处理方法的布里渊频移提取结果几乎一致，二次多项式拟合处理方法计算结果明显小于前者，这与图2非常吻合。比较图3和5可知，实测谱和仿真谱的拟合结果非常吻合。这验证了数值产生谱信号能较好模拟实测谱。由于数值产生信号的布里渊频移已知，据此可知5种处理方法对应误差幅值的均值分别为0.06、0.05、0.05、0.05和1.35MHz。这与基于实测谱分析结果基本吻合，即基于洛伦兹、高斯、伪Voigt和Voigt模型处理方法的准确性足够高，但二次多项式拟合处理方法存在显著误差。也基本上验证了数值产生谱的可靠性。

3.多项式拟合方法的影响因素

要想减少二次多项式拟合处理方法的误差，则需要研究各种因素对该处理方法准确性的影响规律。根据我们***研究发现数值产生信号时谱模型选择对布里渊频移提取结果影响不大，同时为了适当加快计算速度，因此后续分析采用洛伦兹模型产生谱。考虑到实际情况多变，分析时较之实际情况适当扩展了参数的取值范围，这不会减少本发明的可信度。

3.1扫频范围

g₀、v_B和Δv_B分别取0.9、10.7GHz、0.03GHz；信噪比设置为10、20和30dB；扫频点数为61；扫频范围在0.2Δv_B到10Δv_B范围内变化。针对每种参数组合，产生10000个谱信号，计算后得布里渊频移误差幅值的均值。二次多项式拟合处理方法算得的布里渊频移误差如图6所示。

由图6可知，扫频点数固定时当扫频范围较小时布里渊频移的计算误差较大。随着扫频范围的增加误差逐渐减少，当扫频范围为一个线宽时误差达到最小值。然后，随着扫频范围的增加误差又逐渐增大。这是因为扫频范围过小时对应的信号中没有包含足够多的谱特征，当扫频范围过大时在谱特征提取有效范围内的点数太少。因此两种情况下的误差均较大。

3.2扫频点数

扫频范围为Δv_B；扫频点数N在3到501范围内变化。其他参数与3.1节一致，二次多项式拟合处理方法算得的布里渊频移误差如图7所示。

由图7可知，随着扫频点数的增加布里渊频移误差逐渐减少。拟合发现布里渊频移误差幅值的均值E_vB满足：E_vB＝aN^b，以上3种情况下a的值分别为3.2812、0.8580和0.2677，对应b的值分别为-0.5095、-0.4722和-0.4710，拟合相对误差分别为3.11％、4.13％和4.12％。

3.3信噪比

扫频范围为Δv_B；信噪比在0到40dB范围内变化。其他参数与3.1节一致，二次多项式拟合处理方法算得的布里渊频移误差如图8所示。

由图8可知，随着信噪比(SNR,signal to noise ratio)的增加布里渊频移误差逐渐减少。当信噪比大于10dB时布里渊频移误差与信噪比近似满足指数规律变化。当信噪比小于10dB时误差的变化规律与大于10dB时不同，这是因为信噪比较低时直接采用二次多项式拟合处理方法得到的布里渊频移超出了扫频范围，处理方法实现时采用扫频范围内的随机值来作为算得的布里渊频移，这样可以减少测量误差且使测量结果显得更正常。

3.4线宽

扫频范围为Δv_B；Δv_B在0.03到1GHz范围内变化。其他参数与3.1节一致，二次多项式拟合处理方法算得的布里渊频移误差如图9所示。

由图9可知，随着线宽的增加布里渊频移误差线性增大。因此，不能通过选择过窄的入射脉冲光来提高空间分辨率。为了提高布里渊频移的准确性，应该选择适当配置使布里渊谱线宽减少。

3.5扫频范围偏差

由于待测光纤布里渊频移未能完全准确获得，扫频范围的中点与布里渊频移未必重合，二者的差距在本发明方法中称为扫频范围偏差。扫频范围偏差在0到0.3Δv_B范围内变化；扫频范围为Δv_B。其他参数与3.1节一致，二次多项式拟合处理方法算得的布里渊频移误差如图10所示。

由图10可知，随着扫频范围中点逐渐偏离布里渊频移布里渊频移误差有增加的趋势。因此，应该围绕布里渊频移为中心选择谱信号用于布里渊频移提取。

4改进的处理方法及验证

4.1改进处理方法的提出

由上述分析可知，二次多项式拟合处理方法的准确性影响较大的因素分别是扫频范围、扫频点数、信噪比、线宽和扫频范围偏差。虽然扫频点数、信噪比和线宽对二次多项式拟合处理方法的准确性影响较大，但这3个因素主要由实际测量状况决定。而用于拟合谱的扫频范围和扫频范围偏差可以调整接近最优值以减少处理方法的误差。具体来说就是不再是所有测量信号而是以布里渊频移为中心截取1个Δv_B的谱信号用于二次多项式处理方法的拟合。这样改进的二次多项式拟合处理方法的准确性能提高，而且由于待拟合点数的减少二次多项式拟合处理方法的计算速度会进一步加快。

由于布里渊频移为待测量，在处理方法执行前未知，故实际处理方法执行时以布里渊谱峰值对应频率为中心截取1个Δv_B的谱信号用于二次多项式处理方法的拟合。线宽虽然与多个因素，包括光纤密度、折射率、石英光纤材料的粘滞系数、激光器输出中心波长和入射脉冲光宽度等有关，通常实验布置时除了入射脉冲光宽度其他参数基本固定。因此，线宽的估算比较容易。

4.2实测谱

采用2.1节的实测谱进行分析，区别在于本节仅仅以增益峰值为中心截取了1个Δv_B的谱信号。Δv_B取实测结果的近似值，即0.1GHz。5种处理方法的计算结果和计算时间如图11所示。

由图11(a)可知，基于洛伦兹、高斯、伪Voigt和Voigt模型处理方法的布里渊频移提取结果几乎一致，这与2(a)吻合。比较图2(a)与11(a)可知，本发明提出的改进二次多项式拟合处理方法的计算结果与经典谱模型处理方法的计算结果非常接近，它与以上4种处理方法差距的均值分别仅为0.39、0.24、0.23和0.26MHz。考虑到温度系数为1.12MHz/℃，以上布里渊频移误差对应的温度误差仅分别为0.35℃、0.21℃、0.21℃和0.23℃。另外，4种经典处理方法的布里渊频移提取结果在待拟合谱范围调整前后差别很小，也进一步验证了以上4种处理方法的准确性。图11(b)与图2(b)非常相似，也是二次多项式拟合处理方法的计算量远小于其他处理方法。以上结果说明了改进二次多项式拟合处理方法计算量远小于其他经典处理方法，但准确性与其他处理方法相当。

4.3数值产生谱

谱参数及谱产生方法与2.2节一致，但本节以增益峰值为中心截取了0.1GHz的谱信号用于布里渊频移提取。5种处理方法的计算结果如图12所示。

比较图11(a)与12可知，针对5种处理方法实测谱和仿真谱的计算结果非常相似。针对数值产生谱，它与以上4种处理方法差距的均值分别仅为0.34、0.21、0.23和0.24MHz。以上4种经典处理方法及本发明提出的改进处理方法误差幅值的均值分别仅为0.13、0.06、0.06、0.06和0.24MHz，按照典型的温度敏感系数1.2MHz/℃来估算，温度单一因素测量时以上5种处理方法的误差分别为0.11、0.05、0.05、0.05和0.20℃。以上结果进一步验证了本发明提出的改进二次多项式拟合处理方法的有效性。

5结论

基于实测和数值产生的布里渊谱，本发明对基于二次多项式拟合处理方法的布里渊频移提取问题进行了***处理，在保证实时性的基础上显著提高了计算准确性，具体来说结论如下：

1)、当扫频点数固定时随扫频范围增加布里渊频移误差先减少到最小值后逐渐增加，最佳扫频范围为1个线宽；随扫频点数和信噪比增加布里渊频移误差分别成幂和指数规律减少；随线宽增加误差线性增大；随扫频范围偏差增加误差逐渐增大，用于特征提取的谱信号尽量选择围绕布里渊频移左右对称。

2)、二次多项式拟合处理方法的计算速度远快于基于洛伦兹、高斯、伪Voigt和Voigt模型的处理方法，但原始二次多项式拟合处理方法可能存在显著误差。通过选择合适扫频范围的改进二次多项式拟合处理方法(以布里渊谱峰值对应频率为中心截取1个Δv_B的谱信号用于二次多项式处理方法的拟合)可以在保持快速计算的基础上使其准确性与基于洛伦兹、高斯、伪Voigt和Voigt模型的处理方法相似。

3)、本发明实现了光纤分布式传感的快速测量。

Claims (2)

1.一种基于改进二次多项式拟合的布里渊频移快速准确提取方法，其特征是，所述方法是在分布式光纤传感器所输出的测量信号中，以布里渊频移为中心截取1个布里渊增益谱的线宽Δv_B的谱信号，对该信号进行二次多项式处理方法的拟合，得到布里渊谱：

g_B(v)＝av²+bv+c

其中，g_B(v)为布里渊散射谱的增益；v为频率；a，b，c为二次多项式系数，其数值随布里渊谱的变化而不同，则布里渊频移为

2.根据权利要求1所述的基于改进二次多项式拟合的布里渊频移快速准确提取方法，其特征是，在以布里渊频移为中心截取1个布里渊增益谱的线宽Δv_B的谱信号时，所述布里渊频移由布里渊谱峰值对应频率代替。

Images