请选择 进入手机版 | 继续访问电脑版
游客您好
第三方账号登陆
  • 点击联系客服

    在线时间:8:00-16:00

    客服电话

    17712677918

    电子邮件

    857020071@qq.com
  • 扫描二维码

    关注微信公众号

推荐阅读
vguangxian 版主
未知星球 | 未知职业
  • 关注0
  • 粉丝0
  • 帖子361
热议话题
Node.js 简介
2021-08-28 0
Node.js 是一个开源与跨平台的 JavaScript 运行时环境。 它是一个可用于几乎任何项目
精选帖子

基于分布式光纤传感的输电线路温度和应变快速测量方法

[复制链接]
vguangxian 发表于 2022-12-10 10:28:55 | 显示全部楼层 |阅读模式 打印 上一主题 下一主题
摘要

为了提高基于布里渊散射的分布式光纤传感的输电线路温度和应变测量的实时性,为基于分布式光纤传感的电气设备振动测量奠定基础,将斜坡法引入输电线路温度和应变测量中来。介绍了布里渊谱模型,基于谱拟合、单斜坡和双斜坡的布里渊频移计算方法,基于实测裸光纤、光纤复合架空地线和光电复合海底电缆中复合光纤的布里渊谱,系统比较了谱拟合法、单斜坡法和双斜坡法的性能。基于数值产生以及实测的光纤沿线布里渊谱研究了工作点以及信噪比对双斜坡法算得布里渊频移误差的影响规律。结果表明:双斜坡法左右2个最佳工作点为光纤沿线布里渊频移的均值加减半个线宽;随信噪比增加布里渊频移误差成指数规律下降;在高信噪比(≥30 dB)和光纤沿线布里渊频移波动程度不大(明显小于1个线宽)时斜坡法可以有接近谱拟合法的准确性;双斜坡法的准确性高于单斜坡法。较之谱拟合法,单斜坡法和双斜坡法可以显著降低谱测量时间,传感距离较短时有望实现对振动的测量且没有相位敏感光时域反射技术(phase-sensitive optical time domain reflectometer, φ-OTDR)测量结果存在不确定性的问题。


关键词 : 布里渊散射; 分布式光纤传感; 输电线路; 温度; 应变; 快速测量; 双斜坡法;

DOI:10.13336/j.1003-6520.hve.20200792


ABSTRACT

To improve real-time performance of temperature and strain measurement for transmission lines by use of Brillouin-scattering-based distributed optical fiber sensor, and to lay the foundation for vibration measurement based on distributed optical fiber sensing, the slope-assisted techniques are used to measure the temperature and strain of transmission lines. The Brillouin spectrum model is introduced. Based on the measured Brillouin spectra along bare optical fibers, optical fiber composite overhead ground wires (OPGW), and optical fiber composite submarine power cables, the performances of the spectrum fitting method, slope-assisted technique, and double slope-assisted technique are systematically compared. Based on the numerically generated and measured Brillouin spectra along optical fibers, the influences of working point and signal-to-noise ratio (SNR) on the BFS error estimated by the double slope-assisted technique is investigated. The results reveal that the optimal left and right working points of the double slope-assisted technique are vB0±ΔvB/2, where vB0 and ΔvB respectively are the mean BFS along the optical fiber and the linewidth. BFS error exponentially decreases with SNR. If SNR is high (≥30 dB) and BFS variation along optical fiber is low (significantly less than ΔvB), the slope-assisted techniques have similar accuracy with the spectrum fitting method. The measurement time of the slope-assisted techniques is much shorter than that of the spectrum fitting method. The accuracy of the double slope-assisted technique is higher than that of the slope-assisted technique. The vibration measurement for a short fiber under test can be conducted by the slope-assisted techniques and it can avoid the problem of measurement uncertainty of the φ-OTDR (phase-sensitive optical time domain reflectometer).


KEY WORDS : Brillouin scattering; distributed optical fiber sensing; transmission line; temperature; strain; fast measurement; double slope-assisted technique;


0 引言

光纤传感(optical fiber sensing, OFS)是自20世纪70年代发展起来的传感技术。与传统传感技术相比,光纤传感具有质量轻、体积小、抗电磁干扰能力强、安全性高、耐高温、便于融合进物联网等优点[1-2],在特殊环境下可完成传统传感器无法实现的任务,近些年在电气设备状态监测领域逐渐得到应用[3-7]。

光纤传感按原理基本可分为单点光纤传感[8]、准分布光纤传感[9]和全分布光纤传感[10-11]3类。基于布拉格光栅的光纤光栅传感和干涉型光纤传感器属于单点光纤传感技术,其仅可在极小范围内监测选定点附近的参量变化;在单点传感基础上,Ren L将光纤光栅传感技术的传感点连接成1条离散的传感路径,实现了准分布式光纤传感[12]。但该技术只能监测预先设置传感单元位置的状态,并非真正的全分布式传感,且随着传感单元设置数量的增加其成本、施工难度等会显著上升;随着技术的发展,全分布式光纤传感技术逐渐成熟,该类技术除具有传统光纤传感的优点外,还可实现长距离、大范围的传感[13],1根光纤即可代替成千上万个点式传感器,具有很好的应用前景。

目前主要有基于瑞利散射、拉曼散射和布里渊散射的3类分布式光纤传感技术[14]。基于瑞利散射的光时域反射技术(optical time domain reflectometer, OTDR)常用于光纤断点和损耗监测,其对温度和应变不敏感,测量精度较低,且较难实现长距离的监测。基于拉曼散射的光纤传感技术可实现沿线温度的监测,但其无法监测应力且测量范围通常小于30 km,难以实现更长距离监测[15]。基于布里渊散射的分布式光纤传感技术可实现温度和应变的监测,且具有长监测距离、高空间分辨率和高测量精度的特点,其性能明显优于前2种技术[16],在电气设备的监测中具有良好的应用前景。应用布里渊光时域反射(Brillouin optical time domain reflectometer, BOTDR)技术有实现光纤复合海底电缆的应变和温度监测的报道[17-18];基于布里渊光时域反射技术和光纤光栅(fiber Bragg grating, FBG)技术相结合可实现分布式在线温度和应变光纤传感系统[19],基于该系统可实现100 km以内架空输电线路温度的在线监测。

通常来说,基于布里渊散射的传统分布式传感技术采用谱拟合法获得布里渊频移(Brillouin frequency shift, BFS)。谱拟合法为保证测量的精度,扫频点数不会选取的过少,导致测量时间为分钟级[20],甚至更长。虽然通过采用快速布里渊频移提取方法[21-22]可有效减少测量的总耗时,但该方法的实时性仍有待于进一步提高。为减少测量时间,R. Bernini和A. Minardo等在2009年首次提出了基于布里渊散射的单斜坡法(slope-assisted BOTDA, SA-BOTDA)[23],该技术仅需选取1个探测光频率即可解调出布里渊频移,极大缩短了测量时间。2014年,A. Motil提出了基于布里渊分布式光纤传感的双斜坡法(double slope-assisted BOTDA, DSA-BOTDA)[24],该方法通过测量2处探测光频率来解调布里渊频移,一定程度上牺牲了测量速度,但免除了功率波动造成的测量误差。单斜坡法和双斜坡法实现了基于布里渊散射的光纤分布式传感的快速测量,Aldo Minardo在2013年采用单斜坡法实现了铁路的在线动态监测[25];2017年,Damien Maraval将单斜坡法应用于管道振动监测[26]。以上研究显著推动了基于光纤分布式传感的温度和应变快速测量的进步。但该方面尚存在以下问题:1)2种斜坡法以及谱拟合法针对实测布里渊谱的性能比较未见文献公开报道;2)双斜坡法的左右2个工作点以及布里渊谱信噪比对布里渊频移测量准确性的影响未见报道。

本文将谱拟合法、单斜坡法和双斜坡法应用于裸光纤、光纤复合架空地线(optical fiber composite overhead ground wire, OPGW)和光电复合海缆监测,并对效果进行了比较,指出了斜坡法尤其是双斜坡法的特点和优势。在此基础上基于数值产生和实测光纤沿线的布里渊谱研究了左右两侧工作点和布里渊谱信噪比对双斜坡法提取布里渊频移准确性的影响,并确定了最佳工作点。本文工作对基于布里渊散射的分布式光纤温度和应变快速测量技术在输电线路乃至电气设备监测领域的应用具有较好的参考价值。

1 原理介绍1.1 基于布里渊散射的温度和应变测量

通常来说光纤分布式传感中入射光为矩形光脉冲,布里渊谱近似满足伪Voigt模型,表示如下

img_1.png

式中:g01和g02分别为洛伦兹和高斯模型布里渊增益峰值,g0=g01+g02为峰值增益;∆vB和vB分别为.线宽和布里渊频移;gB为增益;v为频率。

伪Voigt模型不仅能较好逼近Voigt模型,而且避免了卷积的计算,计算量不是很大,编程实现也相对容易。因此,也有大量布里渊频移提取方法是基于以上的伪Voigt模型。

谱峰对应频率称为布里渊频移,光纤温度及应变与布里渊频移呈线性关系,如式(2)所示。

img_2.png

式中:CvT和Cvε分别为温度和应变敏感系数;∆T和∆ε分别为温度和应变的变化量;vBO为基准温度和应变下的布里渊频移。

因此,基于温度和应变敏感系数、基准条件下的布里渊频移以及实测的布里渊频移即可解析光纤中的温度或应变。图1展示了布里渊谱及参数。

1.2 基于谱拟合的频移提取方法

采用非线性最小二乘拟合方法,建立目标函数如下

img_3.png

式中:E为所有频率点增益误差的平方和;vi和gBi分别为第i个频率值及对应的布里.渊增益;ei为第i个频率点增益的误差;N为扫频点数。式(3)中的gB可用式(1)计算,对应的方法是基于伪Voigt模型的谱拟合法[27]。

以上的目标函数采用Levenberg-Marquardt算法优化,对应待优化变量的修正计算式为

img_4.png

式中W=[w1, w2, …, wM]T为待优化变量构成的矩阵,即W=[g01, vB, ∆vB, g02];e=[e0, e1, …, eN-1]T为误差矢量矩阵;J为雅克比矩阵;I为M维单位矩阵;l为当前迭代次数。如果迭代后误差E增加,则λ=10λ,同时W的原始值保留,即W(l+1)=W(l)。如果迭代后误差E减少,则λ=λ/10,按照式(4)修正W。λ的初值设置为1。

wi(l)为第l次迭代后得到的第i个参数。如果式(5)满足则l+K为最终的迭代次数。

. img_5.png

K选择为6时式(5)具有不错的效果。

当目标函数趋于最小化后,得到最优化的vB即为求解得到的布里渊频移。

1.3 基于单斜坡的频移提取方法

不同于谱拟合法的完整扫频,单斜坡法只需在布里渊谱上升或下降斜坡选取1个探测光频率(称为工作点[23],其值为v0),即可解调出布里渊频移的信息,极大缩短了谱测量时间。由图2可知,在增益的单调变化区间范围内,增益与布里渊频移存在一一对应的关系,由选定探测光频率的增益值即可算得布里渊频移。

设δvB为布里渊频移的改变量,定义如下

<span class="MathJax" id="MathJax-Element-1-Frame" tabindex="0" data-mathml="δvB=vB1−vB0" role="presentation" style="box-sizing: border-box; display: inline; line-height: normal; text-indent: 0px; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: 0px; min-height: 0px; border: 0px; position: relative;">δvB=vB1−vB0δvB=vB1−vB0 (6)

式中:vB0和vB1分别为基准谱与实测谱的布里渊频移。基准谱为用于获得工作点增益与布里渊频移关系且信噪比足够高的布里渊谱。

以洛伦兹线型为例,工作点选择为布里渊频移的左侧半峰(v0=vB0-ΔvB/2)时,gB与δvB的关系及动态范围(测量允许的范围)如图2所示。

img_6.png 图1 布里渊谱及参数Fig.1 Brillouin spectrum and its parameters

img_7.png 图2 单斜坡法中gB与δvB关系Fig.2 Relationship between gB and δvB in slope-assisted technique

虽然布里渊频移处于理论动态范围内时可以采用单斜坡法计算布里渊频移,但非线性区间内灵敏度偏低或由于引入模型的原因会导致更大的误差。因此,其动态范围也比较有限。

γ为选取探测光频率的参数,探测光频率v0满足

img_8.png

根据基准谱可以确定工作点处增益与布里渊频移的关系,再结合实测条件下工作点处的增益即可获得对应的布里渊频移。研究表明,最佳的γ值约为1[28]。

单斜坡法扫频个数为1,极大地缩短了测量时间,优势十分明显。但受制于布里渊谱的线性区域,其动态范围(可应用范围)较小,超出该范围时噪声的影响变大;同时功率波动的影响给探测光频率处测得的增益值带来了一定误差,对其测量准确性产生了一定的影响。

1.4 基于双斜坡的频移提取方法

与单斜坡法不同,双斜坡法通过在基准谱上升与下降斜坡处各选择1个频率点,根据其对应增益比值的对数来解调布里渊频移信息[24],表示为

. img_9.png

式中:z为光纤位置;k为比例系数(z)为泵浦光功率;v+和v-是工作点的频率(频移),表示为:

img_10.png

img_11.png

显然γ+和γ-取值对准确性会有影响,但其最佳值未见文献研究并明确给出,需要进一步研究。式(8)可简化为

img_12.png

由式(11)可知,双斜坡法通过对光纤2个频移对应的布里渊增益进行比值取对数得到RB,消除掉泵浦光功率因子P(z),避免了功率波动带来的影响。根据式(1)和式(8)—式(11)可知,RB与ΔvB、vB、v+和v-有关。实际应用时ΔvB、v+和v-均为确定值或已知值,因此,根据RB即可计算vB。为了获得ΔvB、v+和v-,该方案也需首先进行完整的扫频以获得基准谱,且扫频范围应尽可能大,频率间隔尽可能小,叠加平均次数一般取到允许情况下的最大值。由基准谱计算得到布里渊频移vB0和线宽ΔvB,同时根据vB0、ΔvB、γ+和γ-确定探测光频率v+和v-。利用基准谱获得RB与δvB的关系,再结合实测谱的RB值即可获得δvB。

.选取初始布里渊频移vB0=10.85 GHz,g01和g02分别为0.5,ΔvB=40 MHz,γ+=γ-=1,由式(11)得到RB与δvB的关系如图3所示。

RB与δvB关系曲线中的单调区间为双斜坡法理论上的动态范围,由图3可知,为1.5ΔvB左右。如果布里渊频移的变化>1.5ΔvB(图中虚线),则有2个δvB解可以满足要求,无法确定正确的布里渊频移。但在区间的边缘,灵敏度较低,噪声的影响导致误差增大,实际上不应该使用,这会使动态范围有所减少。

虽然双斜坡法利用了2个频率点增益的比较抑制了功率波动造成的误差,但比值方式可能会进一步放大增益误差,其与单斜坡相比准确性如何需要系统验证。

布里渊谱测量时间tT如式(12)[29]所示

<span class="MathJax" id="MathJax-Element-2-Frame" tabindex="0" data-mathml="tT=N(tsf+NatT)" role="presentation" style="box-sizing: border-box; display: inline; line-height: normal; text-indent: 0px; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: 0px; min-height: 0px; border: 0px; position: relative;">tT=N(tsf+NatT)tT=N(tsf+NatT) (12)

式中:T为泵浦脉冲光的周期;Nat为信号叠加平均次数;tsf为频率扫描的切换和稳定时间。

由式(12)可知,谱测量时间与扫频点数成正比,从测量耗时的角度来看,由于谱拟合法需要测量几十个频率点,而单斜坡法和双斜坡法仅分别需要测量1和2个频率点。因此,斜坡法的测量时间可能仅为谱拟合法的几十分之一,这当然与扫频点数有关。由于无需涉及谱拟合,它们的另外一个优势是计算速度很快,一旦谱测量完成,很快就能获得光纤沿线的布里渊频移。

img_13.png 图3 双斜坡法中RB与δvB关系Fig.3 Relationship between RB and δvB in double lope-assisted method

2 输电线路监测中的应用

本章采用实测谱比较不同方法的效果。不失一般性,γ+和γ-均取值为1。

2.1 裸光纤

由于输电线路中敷设的光纤多为裸光纤,或含有护套的裸光纤。缠绕于光纤轴上的裸光纤也会承受应变,使得布里渊频移随光纤沿线逐渐变化,能较好模拟输电线路中敷设光纤的状况,而且裸光纤长度可以更加灵活选择,可以采集更多的布里渊谱去充分验证结果的可靠性。因此,本节选择其进行分析。待分析的裸光纤为1盘长度为1 km左右的G657A2型单模光纤,采用BOTDR在室温下测量布里渊谱。扫频范围为10.52~10.92 GHz,扫频间隔为1 MHz,量程为2 km,入射脉宽为200 ns,波长为1 550 nm。如果信噪比太低则双斜坡法误差太大,可能不能满足实际要求,如果信噪比很高则要求叠加平均次数很高,测量的实时性受到限制,综合考虑叠加平均次数为212。采用基于伪Voigt模型的谱拟合法得到峰值处的信噪比为32.53 dB。为了获得信噪比较高情况下的光纤沿线布里渊谱作为基准谱,叠加平均次数选择足够高的218,采用基于伪Voigt模型的谱拟合法获得峰值处的信噪比为45.12 dB,线宽和布里渊频移的均值分别为41.08 MHz和10 724.36 MHz。根据式(7)可得单斜坡法工作点选择为10 704 MHz,根据式(9)—式(10)可得双斜坡法工作点分别选择为10 704和10 745 MHz。

选择1个典型基准谱及对应基于伪Voigt模型的谱拟合法结果如图4所示,可知实测谱与拟合谱符合度较高。这是由于该信号的信噪比足够高且扫频点数也够多,此时采用基于伪Voigt模型的谱拟合法有望得到比较准确的布里渊频移,该值可以作为准确值去评估低信噪比情况下不同方法的误差。

针对叠加平均次数为212的信号,3种方法的布里渊频移计算结果与高信噪比情况(叠加平均次数为218)的比较如图5所示。采用谱拟合法算得的叠加平均次数=218情况下的布里渊频移比叠加平均次数=212情况下的整体偏大,平均大0.28 MHz,分析后认为这一波动是由于温度变化导致,故需修正,图5中的情况是修正后的结果。以谱拟合法针对叠加平均次数=218情况的计算结果为准确值(信噪比高、点数多、准确性高),3种方法的布里渊频移误差分别为7.48×10-2、0.58和0.28 MHz。此时双斜

img_14.png 图4 裸光纤高信噪比的布里渊谱及拟合结果(叠加平均次数为218,入射光脉宽为200 ns)Fig.4 Measured spectrum with high SNR for bare fiber and corresponding fitted curve(average time is 218, incident light pulse width is 200 ns)

坡法的误差为单斜坡法的48.28%。由式(12)可知,测量时间与扫频点数成正比。扫频401个点时光纤沿线20 000个谱的测量时间为541 s,故结合式(12)可知,单斜坡和双斜坡法测量时间仅分别为1.35和2.70 s。斜坡法可有效减少测量时间。

双斜坡法准确性比单斜坡法更高的原因除了其利用2个频率点的信息外,还通过比值消除了不同次测量之间的布里渊增益的波动,而后者对单斜坡法来说几乎无法直接避免。

2.2 光纤复合架空地线

待分析的光纤复合于220 kV的OPGW中,长度约为14 km。实测布里渊谱的扫频范围为10.4~11.2 GHz,扫频间隔为5 MHz。量程为20 km。入射脉宽为20 ns,波长为1 550 nm,叠加平均次数为213。采用基于伪Voigt模型的谱拟合法得到光纤沿线线宽和布里渊频移的均值分别为81.69 MHz和10 830.58 MHz,峰值处布里渊谱的信噪比为30.98 dB。因此,根据式(7)得到单斜坡法工作点选择为10 790 MHz,根据式(9)—式(10)可得双斜坡法工作点分别选择为10 790和10 870 MHz。针对典型的实测谱,采用基于伪Voigt模型谱拟合法的拟合结果如图6所示。3种方法算得的光纤沿线布里渊频移如图7所示。

由图6可知,实测谱虽然含有一定的噪声,但噪声含量并不是很高,且扫频点数众多,测量得到的布里渊谱能充分表征谱特征参数,谱拟合得到的布里渊谱与实测谱能较相符,得到的布里渊频移准确性有保障。

img_15.png 图5 不同方法算得光纤沿线的布里渊频移(裸光纤)Fig.5 Brillouin frequency shift along optical fiber estimated by different methods(bare fiber)

由图7可知,整体上看斜坡法测量得到的光纤沿线的布里渊频移与谱拟合法得到的结果基本相符,尤其是双斜坡法准确性较高。由图7中数据分析可得,单斜坡法和双斜坡法算得布里渊频移与谱

img_16.png 图6 OPGW复合光纤实测谱及对应的拟合结果Fig.6 Measured spectrum of OPGW composited fiber and corresponding fitted curve

img_17.png 图7 不同方法算得OPGW复合光纤沿线的布里渊频移Fig.7 Brillouin frequency shift along optical fiber composited in OPGW estimated by different methods

拟合法的差距分别为2.87和1.02 MHz。考虑到此时谱拟合准确性较高,将谱拟合法得到的布里渊频移作为准确值,则双斜坡法的误差为单斜坡法的35.54%。具体而言单斜坡法在0~3 km左右得到的布里渊频移与谱拟合法结果差距偏大,这是因为单斜坡法线性范围较小,进入了非线性区间,误差增大。

普通单模光纤的温度和应变的敏感系数非常稳定,分别取值为1.06 MHz/℃和0.047 7 MHz/10-6。但针对不同类型和批次的光纤基准温度下的布里渊频移会有一定的波动,由于该光纤已经复合于OPGW中无法取出后测量,故查阅文献后取基准温度下布里渊频移为10.8 GHz(即实际工程应用时在厂家复合光纤之前可实测基准温度下的布里渊频移)。考虑到OPGW沿线的温度近似为同一值,其光纤沿线布里渊频移的变化主要由应变导致,结合式(2)可得3种方法算得的光纤应变如图8所示。

本案例如果使用谱拟合法,则测量OPGW复合光纤上所有布里渊谱需要的时间为595.75 s。扫频点数为161,故结合式(12)可知,如果采用单斜坡和双斜坡法则对应的测量时间仅分别为3.61和7.22 s。

2.3 光电复合海缆

待分析的为110 kV单相高压光电复合海底电缆,实测布里渊谱的扫频范围为10.71~11.03 GHz,扫频间隔为10 MHz。入射脉宽为200 ns,波长为1 550 nm,叠加平均次数为213。采用基于伪Voigt模型的谱拟合法计算得到光纤沿线线宽和布里渊频移的均值分别为62.36 MHz和10 873.51 MHz,峰值处布里渊谱的信噪比为37.41 dB。因此,根据式(7)可得单斜坡法工作点选择为10 840 MHz,根据式(9)—式(10)可得双斜坡法工作点分别选择为10 840和10 900 MHz。3种方法算得的光纤沿线布里渊频移如图9(a)所示。对海缆中复合光纤进行了标定,结果表明在基准温度下的布里渊频移νBO=10.856 GHz,应变敏感系数为0.047 7 MHz/10-6。考虑到实际情况光纤沿线的布里渊频移主要是由应变导致,如果认为光纤沿线的温度为基准值,结合式(2)得到以上情况下不同方法求得光纤沿线的应变见图9(b)。

单斜坡法和双斜坡法算得布里渊频移与谱拟合法的差距分别为0.57和0.47 MHz。考虑到此时谱拟合准确性较高,近似忽略其误差,故双斜坡法的误差为单斜坡法的82.46%。结合图9可知,与2.1和2.2节类似,双斜坡法的布里渊频移、应变的测量准确性要(略)高于单斜坡法。

本案例如果使用谱拟合法,则测量海缆复合光纤上所有布里渊谱需要的时间为24.63 s。扫频点数为33,故结合式(12)可知,如果采用单斜坡和双斜坡法则对应的测量时间仅分别为0.75 s和1.49 s。由于温度和应变能有效反映输电线路状态,而根据本文研究结果基于光纤分布式传感能快速测量输电线路中复合光纤的温度和应变,因此能为其状态的快速评估奠定基础。当然,如果希望温度和应变同时测量,则需要再引入一个变量,比如光强,这时候图8和9(b)不再成立,但是这并不影响本文的布里渊频移快速测量方法的有效性,也可以为温度和应变的同时、快速测量提供参考。

由图5可知,对于1 km光纤,叠加平均次数Nat=212时则双斜坡法的准确性有一定的保证,这样

img_18.png 图8 不同方法算得OPGW复合光纤沿线的应变Fig.8 Strain along optical fiber composited in OPGW estimated by different methods

img_19.png 图9 不同方法算得的海缆复合光纤布里渊频移和应变Fig.9 Brillouin frequency shift and strain along optical fiber composited in submarine cable estimated by different methods

如果光纤长度仅为100 m同时保持叠加平均次数不变则双斜坡法的准确性可以进一步提高。考虑到频率扫描的切换和稳定时间tsf的典型值为300 μs[30];泵浦脉冲光的周期T=0.97 μs(光在100 m光纤上走1个来回);而双斜坡法N=2;故根据式(12)可得双斜坡法对应的谱测量时间为8.57 ms。双斜坡法计算布里渊频移时间<0.1 ms,也就是说此时双斜坡法有望能测量100 Hz级别的振动。虽然基于相位敏感光时域反射技术(phase-sensitive optical time domain reflectometer,φ-OTDR)也可以测量振动甚至高频振动,但它存在的关键问题是光纤应变与瑞利背向散射光相位关系存在模糊性[31],导致其无法准确测量振动幅值,基于布里渊散射的分布式光纤传感器的优点是可以准确测量振动幅值。

在实际应用中,如长距离OPGW等监测对象中,整根光纤往往由多段具有不同布里渊频移的短光纤组成,不同光纤段之间的布里渊频移之差可能较大。如果差别特别大则这种情况不是特别适合于斜坡法,其实也给谱拟合法有效提取布里渊频移带来了一些难度。此时可以通过适当减少入射光脉宽、增加布里渊谱的线宽[32]来提高斜坡法的可用性。

虽然以上案例均针对输电线路,但为了监测电气设备运行状态也可以在其中敷设光纤,本文方法同样也可以监测电气设备光纤沿线的温度和应变,从而为其状态评估奠定基础。

3 影响因素研究

工作点的位置和布里渊谱信噪比会影响双斜坡法的准确性,本章对这2个关键因素的影响进行研究以期能为实际应用时提供参考。

3.1 工作点

本节分别用数值产生布里渊谱和实测布里渊谱分析γ+和γ-对双斜坡法准确性的影响。

采用式(1)的伪Voigt模型数值产生布里渊谱,线宽为40 MHz,g01和g02分别设置为0.6和0.2,扫频间隔为1 MHz。按照实测布里渊谱噪声特性对数值产生的谱添加高斯白噪声,为了与实测谱的信噪比一致,峰值处信噪比设置为32.53 dB。设线路长度为2 km,按取样分辨率为2 m产生布里渊谱。第1个测点的布里渊频移为10.712 GHz,沿着光纤到1 km时线性增加到10.72 GHz,然后又沿着光纤到2 km时布里渊频移线性减少到10.712 GHz,布里渊谱扫频范围为10.64~10.8 GHz。γ+和γ-均在0.5~2范围内变化,得到光纤沿线布里渊频移误差均值与γ+和γ-的关系如图10所示。

由图10可知,随着γ+、γ-从较小值开始增加,布里渊频移误差先减少后增加。最小误差为0.30 MHz,γ+=γ-=1时误差为0.31 MHz,此时误差基本达到或接近最小值。

实测谱如2.1节所示,γ+和γ-在0.5~2范围内变化时采用双斜坡法计算光纤沿线的布里渊频移,光纤沿线布里渊频移平均误差与γ+和γ-的关系如图11所示。最小误差为0.28 MHz,γ+=γ-=1时误差为0.29 MHz,此时误差基本达到或接近最小值。显然,仿真结果与实验结果相符,即γ+=γ-=1时准确性最高。注意,图11中双斜坡法误差计算是将叠加平均次数选择足够高的218时宽扫频范围时谱拟合法的结果作为准确值而算得。

img_20.png 图10 双斜坡法布里渊频移误差与γ+和γ-的关系(仿真谱)Fig.10 Change of error in Brillouin frequency shift estimated by double slope-assisted technique with γ+ and γ-(numerically generated spectra)

img_21.png 图11 双斜坡法布里渊频移误差与γ+和γ-的关系(实测谱)Fig.11 Change of error in Brillouin frequency shift estimated by double slope-assisted technique with γ+ and γ-(measured spectra)

3.2 信噪比

γ+和γ-均设置为1,信噪比在20~40 dB范围内变化,不考虑噪声情况下的光纤沿线布里渊频移的变化情况如图12所示,其他参数与3.1节一致。图中波动小、波动中、波动大分别对应光纤沿线最大布里渊频移与最小布里渊频移的差距为8、24、40 MHz。采用式(1)的伪Voigt模型数值产生布里渊谱,采用双斜坡法的布里渊频移计算误差如图13所示,同时图中给出了布里渊频移误差采用式(13)的拟合结果。

<span class="MathJax" id="MathJax-Element-3-Frame" tabindex="0" data-mathml="EBFS=ae−bRSN" role="presentation" style="box-sizing: border-box; display: inline; line-height: normal; text-indent: 0px; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: 0px; min-height: 0px; border: 0px; position: relative;">EBFS=ae−bRSNEBFS=ae−bRSN (13)

式中:EBFS为布里渊频移误差幅值的均值;RSN为峰值处的信噪比;a和b均为待定系数。

由图13可知,针对仿真布里渊谱,随着信噪比的增加双斜坡法算得的布里渊频移误差呈指数规律下降。另外,当信噪比>30 dB且光纤沿线布里渊频移波动不大(明显小于1个线宽)时斜坡法误差已然不大(与第2章结果相符),尤其是双斜坡法的布里渊频移误差接近于0.5 MHz,已经较小,也与谱拟合法的结果较为接近了。进一步地本文也采用实测布里渊谱进行了验证,入射脉宽设置为10、20和200 ns,叠加平均次数分别为210、212、214、216、218。采用双斜坡法的布里渊频移误差与信噪比的关系如图14所示,同时图中也给出了布里渊频移误差采用式(13)的拟合结果。

由图14可知,与仿真布里渊谱类似,针对实测布里渊谱,随着信噪比的增加双斜坡法算得的布里渊频移误差呈指数规律下降。

4 结论

1)双斜坡法左右2个最佳工作点为光纤沿线布里渊频移的均值加减半个线宽;随着信噪比的增加双斜坡法误差呈指数规律下降。

2)谱拟合法具有最高的准确性,在信噪比较高(不建议小于30 dB)且光纤沿线布里渊频移波动较小(明显小于1个线宽)时,斜坡法和双斜坡法也能达到与谱拟合法相似的准确性,尤其是双斜坡法准确性要高于单斜坡法(前者误差为后者的35.54%~82.46%);谱测量时间可以减少为典型设置下谱拟合法的几十分之一,甚至更小。

基于本文工作有望实现基于布里渊散射的光纤分布式传感的输电线路温度和应变快速测量,甚至传感距离较短时电气设备振动的测量,可为电气设备状态快速评估奠定基础。

img_22.png 图12 光纤沿线的布里渊频移(仿真)Fig.12 Brillouin frequency shift along the optical fiber(simulation)

img_23.png 图13 双斜坡法误差与信噪比的关系(仿真谱)Fig.13 Change of error in Brillouin frequency shift estimated by double slope-assisted technique with SNR(numerically generated spectra)

img_24.png 图14 双斜坡法误差与信噪比的关系(实测谱)Fig.14 Change of error in Brillouin frequency shift estimated by double slope-assisted technique with SNR(measured spectra)

参考文献

[1] 廖延彪,苑立波,田芊. 中国光纤传感40年[J]. 光学学报,2018,38(3):1-21. LIAO Yanbiao, YUAN Libo, TIAN Qian.China fiber optic sensing for 40 years[J]. Acta Optica Sinica, 2018, 38(3): 1-21.

[2] 张智娟,陈飞飞,徐志钮,等. 分布式光纤传感中布里渊动态光栅研究现状与发展趋势[J]. 半导体光电,2018,39(4):455-461. ZHANG Zhijuan, CHEN Feifei, XU Zhiniu, et al.Research status and development trends of Brillouin dynamic grating in distributed optical fiber sensing[J]. Semiconductor Optoelectronics, 2018, 39(4): 455-461.

[3] XU Z N, HU Z W, ZHAO L J, et al.Application of temperature field modeling in monitoring of optic-electric composite submarine cable with insulation degradation[J]. Measurement, 2019, 133: 479-494.

[4] 张施令,姚强,李新田. 基于法兰式光纤传感技术的高压组合电器HF气体在线检测[J]. 高电压技术,2019,45(2):402-409. ZHANG Shiling, YAO Qiang, LI Xintian.On-line monitoring of HF gas in high voltage GIS based on flange optical fiber sensing technology[J]. High Voltage Engineering, 2019, 45(2): 402-409.

[5] 阳林,史尊伟,黄吉超,等. 基于BOTDR分布式光纤传感技术的架空线路温度检测现场试验研究[J]. 高电压技术,2015,41(3):925-930. YANG Lin, SHI Zunwei, HUANG Jichao, et al.Field test research on temperature detection of overhead transmission line based on BOTDR distributed fiber sensing technology[J]. High Voltage Engineering, 2015, 41(3): 925-930.

[6] 徐志钮,胡志伟,赵丽娟,等. 光电复合海缆温度场建模分析及在在线监测中的应用[J]. 红外与激光工程,2018,47(7):0717003. XU Zhiniu, HU Zhiwei, ZHAO Lijuan, et al.Analysis of temperature field model in optic-electric composite submarine cable and its applications in online monitoring[J]. Infrared and Laser Engineering, 2018, 47(7): 0717003.

[7] 吴雪瑞,江军,汪卓玮,等. 基于微纳光纤倏逝场传感的变压器油中微水含量检测[J]. 高电压技术,2020,46(6):1955-1961. WU Xuerui, JIANG Jun, WANG Zhuowei, et al.Detection of moisture in transformer oil based on micro-nano fiber evanescent field sensing[J]. High Voltage Engineering, 2020, 46(6): 1955-1961.

[8] 姚维强,司文荣,吕佳明,等. EFPI光纤超声传感器及其潜在局放检测应用综述[J]. 高电压技术,2020,46(6):1855-1866. YAO Weiqiang, SI Wenrong, LÜ Jiaming, et al.Review on EFPI fiber-based ultrasonic sensors and its potentially partial discharge detection application[J]. High Voltage Engineering, 2020, 46(6): 1855-1866.

[9] 王永千,侯继勇,孟凡勇,等. 基于光纤光栅传感的低功耗小型化电网智能感知系统[J]. 电力信息与通信技术,2020,18(4):51-57. WANG Yongqian, HOU Jiyong, MENG Fanyong, et al.Low power consumption miniaturized power grid intelligent sensing system based on fiber Bragg grating sensing[J]. Electric Power Information and Communication Technology, 2020, 18(4): 51-57.

[10] 郭经红,梁云,陈川,等. 电力智能传感技术挑战及应用展望[J]. 电力信息与通信技术,2020,18(4):15-24. GUO Jinghong, LIANG Yun, CHEN Chuan, et al.Challenge and application prospect of power intelligent sensor technology[J]. Electric Power Information and Communication Technology, 2020, 18(4): 15-24.

[11] ZHAO L J, XU Z N, LI T.Birefringence-induced frequency shift error estimation for Brillouin dynamic grating reflected spectra[J]. Optik, 2020, 220: 165164.

[12] REN L, SONG G, CONDITT M, et al.Fiber Bragg grating displacement sensor for movement measurement of tendons and ligaments[J]. Applied Optics, 2007, 46(28): 6867-6871.

[13] 冯欣,曾周末,封皓,等. 基于小波信息熵的分布式振动传感系统的扰动评价方法[J]. 光学学报,2013,33(11):47-53. FENG Xin, ZENG Zhoumo, FENG Hao, et al.A disturbance evaluation method for distributed vibration sensing system based on wavelet information entropy[J]. Acta Optica Sinica, 2013, 33(11): 47-53.

[14] 于婷婷. 分布式光纤传感系统测试平台的设计[D]. 北京:北京邮电大学,2019. YU Tingting.Design of test platform for distributed optical fiber sensing system[D]. Beijing, China: Beijing University of Posts and Telecommunications, 2019.

[15] 陶一鸣. 基于布里渊效应的分布式光纤传感器技术研究[D]. 北京:北京邮电大学,2017. TAO Yiming.Research on distributed fiber optic sensor based on Brillouin effect[D]. Beijing, China: Beijing University of Posts and Telecommunications, 2017.

[16] 刘文哲. 基于布里渊效应的分布式传感信号处理关键算法的研究[D]. 秦皇岛:燕山大学,2015. LIU Wenzhe.Research on the key algorithm of distributed sensing signal processing based on Brillouin effect[D]. Qinhuangdao, China: Yanshan University, 2015.

[17] 李永倩,赵丽娟,杨志,等. 基于BOTDR的海缆3D立体监测系统设计与实现[J]. 仪器仪表学报,2014,35(5):1029-1036. LI Yongqian, ZHAO Lijuan, YANG Zhi, et al.Design and realization of the submarine cable three-dimensional monitoring system based on BOTDR[J]. Chinese Journal of Scientific Instrument, 2014, 35(5): 1029-1036.

[18] ZHAO L J, LI Y Q, XU Z N, et al.On-line monitoring system of 110 kV submarine cable based on BOTDR[J]. Sensors and Actuators A: Physical, 2014, 216: 28-35.

[19] LUO J B, HAO Y P, YE Q, et al.Development of optical fiber sensors based on Brillouin scattering and FBG for on-line monitoring in overhead transmission lines[J]. Journal of Lightwave Technology, 2013, 31(10): 1559-1565.

[20] MOTIL A, BERGMAN A, TUR M.State of the art of Brillouin fiber-optic distributed sensing[J]. Optics Laser Technology, 2016, 78: 81-103.

[21] 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.

[22] 徐志钮,胡宇航,赵丽娟,等. 基于改进二次多项式拟合的布里渊频移快速高精度提取算法[J]. 光谱学与光谱分析,2020,40(3):842-848. XU Zhiniu, HU Yuhang, ZHAO Lijuan, et al.Fast and highly accurate Brillouin frequency shift extracted algorithm based on modified quadratic polynomial fit[J]. Spectroscopy and Spectral Analysis, 2020, 40(3): 842-848.

[23] BERNINI R, MINARDO A, ZENI L.Dynamic strain measurement in optical fibers by stimulated Brillouin scattering[J]. Optical Letters, 2009, 34(17): 2613-2615.

[24] MOTIL A, DANON O, PELED Y, et al.Pump-power-independent double slope-assisted distributed and fast Brillouin fiber-optic sensor[J]. IEEE Photonics Technology Letters, 2014, 26(8): 797-800.

[25] MINARDO A, PORCARO G, GIANNETTA D, et al.Real-time monitoring of railway traffic using slope-assisted Brillouin distributed sensors[J]. Applied Optics, 2013, 52(16): 3770-3776.

[26] MARAVAL D, GABET R, JAOUEN Y, et al.Dynamic optical fiber sensing with Brillouin optical time domain reflectometry: application to pipeline vibration monitoring[J]. Journal of Lightwave Technology, 2017, 35(16): 3296-3302.

[27] 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.

[28] 徐志钮,胡宇航,赵丽娟,等. 基于单斜坡法的光电复合海缆温度、应变快速测量方法[J]. 电力自动化设备,2020,40(5):202-208. XU Zhiniu, HU Yuhang, ZHAO Lijuan, et al.Rapid temperature and strain measurement method for optic-electric composite submarine cable based on slope-assisted method[J]. Electric Power Automation Equipment, 2020, 40(5): 202-208.

[29] 涂晓波. 光纤布里渊矢量分布式传感技术研究[D]. 长沙:国防科学技术大学,2015. TU Xiaobo.Distributed fiber sensing technique employing vector Brillouin optical time domain analyzer[D]. Changsha, China: National University of Defense Technology, 2015.

[30] FARAHANI A M, CASTILLO-GUERRA E.Acceleration of measurements in BOTDA sensors using adaptive linear prediction[J]. IEEE Sensors Journal, 2013, 13(1): 263-272.

[31] 李永倩,马立. φ-OTDR振动传感系统中的相位模糊问题[J]. 红外与激光工程,2015,44(9):2800-2804. LI Yongqian, MA Li.Phase ambiguity problem in φ-OTDR vibration sensing system[J]. Infrared and Laser Engineering, 2015, 44(9): 2800-2804.

[32] CUI Q, PAMUKCU S, XIAO W, et al.Truly distributed fiber vibration sensor using pulse base BOTDA with wide dynamic range[J]. IEEE Photonics Technology Letters, 2011, 23(24): 1887-1889.







上一篇:光纤传感器在电力系统的应用
下一篇:分布式光纤传感技术应用电力光缆网监测
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

免责声明
icitu.com 信息谷-icit软硬件产品营销与行业智慧化解决方案实施供需对接空间,信息人家园-信息人之间互助互惠互利的资源整合空间,icit人团结起来帮助人们沟通与建设,icit unite to help people communicating&constructing。
0512-68181818
关注我们
  • 访问移动手机版
  • 官方微信公众号

Powered by Discuz! X3.4© 2001-2013 Comsenz Inc. 优康设计 ( 苏ICP备10066271号-2 )|网站地图