用于动态应变测量的快速分布式布里渊光纤传感 下载: 1330次特邀综述
1 引言
自20世纪60年代以来,伴随着激光器和光纤的发明及应用,光纤传感技术也快速发展起来。目前,光纤传感技术已广泛应用于土木结构健康监测、现代工业控制及**安全等诸多领域。与传统电阻式、电容式等电子传感技术相比,光纤传感技术具有很明显的优势:1)光纤既是损耗很小的传光介质,又是灵敏的传感器件,可以认为是由成千上万个传感单元串联而成,因此能够进行长距离和分布式的传感测量;2)光纤易于接入全光网络,而且质量轻、体积小、便于布设;3)光纤具有抗电磁干扰、耐腐蚀的特性,能够在恶劣的环境下测量;4)光纤传感器具有分辨率高、灵敏度高、测量范围宽等优点[1]。由于光纤传感器具有以上优点,因此在现代化的工程项目中,都离不开各式各样的光纤传感器。值得注意的是,在现代工业及基础设施的结构健康监测中,急需一种可以实时监测应变或温度信息的分布式超快测量技术,因而光纤传感器可以应用于此。
自20世纪90年代以来,基于布里渊散射的分布式光纤传感器被广泛报道[2-5]。通常地,两束光波(抽运光和探测光)相向注入待测光纤(FUT),当两束光波的频率失谐量接近FUT的布里渊频移(BFS)时,由于受激布里渊散射(SBS)效应[6],光能量从高频光波转移到低频光波。布里渊增益谱(BGS)可以通过扫频获得,再通过洛伦兹曲线或高斯曲线拟合获得光纤的BFS。值得注意的是,由于BFS和光纤的应变及温度呈线性关系[2, 4,7-8],因此可以通过解调BFS来计算出对应的应变或温度。
基于布里渊散射的光纤传感方案主要包括:布里渊光时域分析(BOTDA)技术[9-14]、布里渊光时域反射(BOTDR)技术[15-17]、布里渊光相关域分析(BOCDA)技术[5, 7, 18-22]、布里渊光相关域反射(BOCDR)技术[5, 23-25]、布里渊光频域分析(BOFDA)技术[26-27]和布里渊光频域反射(BOFDR)技术[28]。采用反射技术的光纤传感方案,其布里渊信号较弱,信噪比(SNR)差;采用光相关域技术的光纤传感方案,通常采用点式传感,需要通过扫描相关峰(传感点)的位置[5, 23, 25]才能获得分布式的应变或温度测量,这样会极大地限制动态应变的采样率;采用光频域技术的光纤传感方案,光波调制及解调算法较复杂。
本文主要针对传统的BOTDA系统测量时间较长,仅适用于静态或缓慢变化应变测量的局限,从理论上分析BOTDA系统,并综述近些年BOTDA系统在快速分布式应变测量方面的进展。
2 基本原理
2.1 受激布里渊散射
光束通过某一介质,当传输介质光学特性不均匀时,就会产生光散射现象[6]。当入射光功率超过受激布里渊散射阈值时,入射光能量就会大量转移给后向传播的散射光(斯托克斯光),由于存在多普勒效应,斯托克斯光有一个BFS。对于普通的单模光纤,BFS[10, 29-30]的公式为
式中
式中
2.2 BOTDA系统的物理模型
图 1. 传统的BOTDA系统原理示意图。 (a)时域关系; (b)频域关系
Fig. 1. Schematic of classical BOTDA system. (a) Time-domain relationship; (b) frequency-domain relationship
传统的BOTDA系统(
麦克斯韦波动方程可以描述抽运光和探测光的光场,Navier-Stokes方程给出光纤介质声波场的变化。在忽略横向场变化并采用慢变包络近似的情况下,BOTDA系统可由以下三波耦合方程组描述为
式中
式中
2.3 动态分布式应变测量的限制因素
由于传统的BOTDA系统的分布式应变或温度的测量时间通常需要数秒到数分钟,因此该系统适合于静态或缓慢的应变测量,应变的最大振动采样率为
由(6)式可知,影响测量时间的限制因素[35]主要有以下4个方面。
1)FUT的长度
2)平均次数
3)有限的光波频率切换时间
4)扫频个数
限制因素1是固有限制,决定了BOTDA系统的动态测量极限值,目前还没有突破的方案。限制因素2是普遍的限制因素,也可以通过增加抽运光功率和探测光功率以增加SNR,降低对平均次数的要求;另外,对于单模光纤,可以通过两束正交偏振的探测光合成的布里渊信号来消除偏振噪声,从而降低平均次数。限制因素3中的光波频率切换时间主要受限于电学微波源或者任意波形发生器(AWG)的性能,最近有报道指出基于高性能AWG的捷变频技术(FAT)可以有效地压缩光波频率切换时间。对于限制因素4,较多的扫频个数会增加采集的BGS的精度或增大应变变化范围,但是会限制采集时间,最近有报道指出利用斜坡法可以增大频率间隔,有效地减少扫频个数。
3 用于动态应变测量的快速BOTDA的分类
到目前为止,为了提升BOTDA系统的动态分布式传感性能,研究人员提出很多改进方案:偏振补偿技术、光学捷变频技术、斜坡法、光学啁啾链技术、光学频率梳技术等,相对应的动态传感系统分别介绍如下。
3.1 基于偏振补偿技术的快速BOTDA系统
由于SBS效应对抽运光和探测光的偏振态十分敏感,因此使得BOTDA系统中会产生偏振噪声。对于普通的单模光纤,通常需要使用扰偏器和数千次平均来减小系统的偏振噪声,这极大地增加了测量时间,限制了应变振动的采样率。针对限制因素2,2014年,López-Gil等[36-37]提出一种利用平衡探测技术[38]和正交的双边带(DSB)探测光(斯托克斯-反斯托克斯光)来消除偏振噪声的技术方案。基于平衡探测和偏振消除的BOTDA系统的实验装置如
图 2. 基于平衡探测和偏振消除的BOTDA系统的实验装置
Fig. 2. Experimental setup of BOTDA system with balanced detection and polarization noise elimination
图 3. 基于平衡探测和偏振消除的BOTDA系统的实验结果。 (a)单探测器获得的布里渊增益信号(蓝色)、单探测器获得的布里渊衰减信号(红色)、平衡探测器获得的合成布里渊信号(绿色); (b) 25 km处三种信号的局部放大图
Fig. 3. Experimental results of BOTDA system with balanced detection and polarization noise elimination. (a) Brillouin gain trace obtained using single-detector (blue), Brillouin attenuation trace obtained using single-detector (red), and synthetic Brillouin trace obtained using balanced-detector (green); (b) local large diagram of the three traces around 25 km
随后,如
图 4. 基于偏振补偿的应变快速采集测量实验结果
Fig. 4. Experimental results of fast-acquisition measurement of the strain based on polarization compensation
3.2 基于光学捷变频技术的快速BOTDA系统
针对限制因素3,2012年Peled等[40]在传统BOTDA系统的基础上,提出了基于光学FAT的快速BOTDA系统,采用高性能的电学数字式AWG,实现了光波频率的快速切换。首先,将100个扫频光波段所对应的电信号波形以首尾串联的形式提前写入AWG的内存中,通常电信号频率在FUTBFS(约为11 GHz)附近。然后,内存中的波形转化为电信号输出,同时驱动调制器加载到单频率光波,在最小工作点状态,形成一阶USB和一阶LSB。最后,通过滤出一阶USB或一阶LSB,即可获得100个首尾串联的扫频光波段。这样,由于AWG输出电学信号的频率转换时间小于纳秒数量级,因此对应的光波频率切换时间也小于纳秒数量级。与其他因素相比,光波频率切换时间
基于FAT的快速BOTDA系统原理示意图如
图 5. 基于FAT的快速BOTDA系统原理示意图。 (a)探测光经过FAT调制; (b)抽运脉冲光经过FAT调制
Fig. 5. Schematic of fast BOTDA system based on FAT. (a) Probe wave modulated by FAT; (b) pump pulse modulated by FAT
为了降低对AWG高带宽(>10 GHz)和高空间分辨率的实际应用需求,2014年,本课题组[12]结合差分双脉冲技术[13, 41]和二阶边带调制技术,实现了高空间分辨率的快速BOTDA系统。如
随后,二阶LSB探测光波经过FAT技术调制为扫频光波段。抽运光调制为52 ns与50 ns的脉冲对,对应200 mm的高空间分辨率。然而由于每个光波段持续时间最短为2倍2
为了进一步降低对AWG带宽的要求,本课题组又提出了基于双调制方案的捷变频技术[42]。如
图 7. 振动信号的测量。 (a)测量所得光纤振动截面BGS随时间的演进图; (b)应变的振动波形; (c)应变的振动波形功率谱
Fig. 7. Measurements of the vibration signals. (a) Evolution of the measured BGS of the vibrated section of fiber; (b) vibration waveforms of the strain; (c) power spectra of the vibration waveforms of the strain
图 8. 基于双调制方案的捷变频技术示意图
Fig. 8. Schematic of the frequency-agile technique based on dual-modulation scheme
3.3 基于斜坡法的快速BOTDA系统
针对限制因素4,2009年,Bernini等[43]利用BGS的单侧斜坡,搭建了斜坡辅助BOTDA(SA-BOTDA)系统,其原理如
图 9. SA-BOTDA系统原理示意图。 (a)斜坡法操作原理; (b)斜坡法和曲线拟合法应变振动波形测量的比较; (c)不同振动频率下谱域中的动态应变
Fig. 9. Schematic of SA-BOTDA system. (a) Operation principle of slope-assisted method; (b) comparison of strain vibration waveforms measured by slope-assisted method and curve fitting method; (c) dynamic strain in the spectral domain for different vibration frequencies
对于SA-BOTDA系统而言,其振动采样率只受限于光纤长度和平均次数,而且扫频个数
3.3.1 增加斜坡数量
2011年,Peled等[44]针对光纤的BFS分布不均匀问题,提出了一种任意BFS分布的SA-BOTDA系统。首先,用传统的BOTDA系统测量静态时分布式的BGS,通过曲线拟合获得分布式的BFS;然后,根据BFS的分布情况,通过AWG对探测光的频率随时间进行调制;最后,通过同步抽运光与探测光,使得二者频率差值始终位于BGS峰值的3 dB位置。这样光纤中每个位置都可以进行基于斜坡法的分布式的应变动态测量。实验中采集到了高达400 Hz的振动信号。随后,采用BGS双斜坡进行了BOTDA系统的动态测量。定义
图 10. 多斜坡F-BOTDA系统原理示意图。 (a)多斜坡F-BOTDA的抽运光和探测光的频域关系,以及测量范围扩展的基本思路; (b)多斜坡F-BOTDA的抽运光和探测光的时域关系
Fig. 10. Schematic of multi-slope-assisted F-BOTDA system. (a) Frequency-domain relationship of pump and probe waves of multi-slope-assisted F-BOTDA, and basic idea of the extension of measurement range; (b) time-domain relationship of pump and probe waves of multi-slope-assisted F-BOTDA
2016年,本课题组利用双斜坡解调算法和捷变频技术相结合的技术方案,实现了多斜坡F-BOTDA系统[3],其原理如
由(8)式可知,
基于多斜坡 F-BOTDA所得BFS应变振动波形的实验测量结果如
图 11. 多斜坡法(蓝色实线)和传统曲线拟合法(红色虚线)所得BFS应变振动波形的实验测量结果,分别对应不同的扫频间隔。 (a) 20 MHz; (b) 40 MHz; (c) 60 MHz; (d) 80 MHz
Fig. 11. Experimental measurement results of BFS strain vibration waveforms obtained via multi-slope-assisted method (blue solid curve) and traditional curve-fitting method (red dotted curve) in varios frequency-intervals. (a) 20, (b) 40, (c) 60, and (d) 80 MHz
3.3.2 扩大斜坡的单调范围
2012年,Urricelqui等[45]首次提出利用BPSS和射频(RF)解调算法实现动态布里渊光纤传感的技术方案。2014年,Tu等[46]通过外差探测方法和同相正交(IQ)解调算法获得了分布式的BGS和BPSS,搭建了矢量BOTDA系统,并对比了二者在单斜坡BOTDA系统中的动态应变传感。2016年,本课题组根据矢量SBS,提出一种具有大频率范围的单调斜坡的
当抽运光为矩形脉冲光时,存在关系式为
式中
由(4)式、(5)式、(9)式和(10)式可知,
当抽运光为连续光时,(9)式可以化简为
此时,KS是一条直线,斜率为2/Δ
图 12. ΔνB=30 MHz时的归一化谱仿真图,抽运脉冲光脉宽从1 ns到连续光。 (a)归一化的BGS; (b)归一化的BPSS; (c)归一化的KS
Fig. 12. Simulation diagram of normalized spectra with the pump pulse width from 1 ns to continuous wave (ΔνB=30 MHz). (a) normalized BGS; (b) normalized BPSS; (c) normalized KS
基于KS的多斜坡BOTDA系统示意图如
图 13. 多频率方案抽运光和探测光的时域关系
Fig. 13. Time-domain relationship of pump wave and probe wave for multi-slope-assisted method
BFS应变振动波形的实验测量结果如
图 14. BFS应变振动波形的实验测量结果,其中两种参考波形分别为BGS曲线拟合所得参考波形(黑色散点)和KS曲线拟合所得参考波形(蓝色散点),多斜坡F-BOTDA系统解调所得应变振动波形(红色实线)分别对应不同的扫频间隔。 (a) 80 MHz; (b) 120 MHz; (c) 160 MHz; (d) 180 MHz
Fig. 14. Experimental results of BFS strain vibration waveforms, in which two kinds of reference strain vibration waveforms are curve-fitted based on the BGS (black scatter) and KS (blue scatter) respectively, meanwhile the strain vibration waveforms(red solid curve) are demodulated via F-BOTDA system in various frequency-intervals. (a) 80, (b) 120, (c) 160, and (d) 180 MHz
3.4 基于光学啁啾链技术的快速BOTDA系统
在传统的FAT方案中,只有扫频光波段和对应的抽运脉冲光相互作用后才能获得分布式的BGS,无法实现由单发抽运脉冲光获得分布式的BGS。最近,本课题组提出了一种基于光学啁啾链(OCC)的单发测量的动态布里渊光纤传感技术方案[47]。OCC-BOTDA系统原理示意图如
由(11)式可知,该技术方案的应变振动采样率只受限于光纤长度和平均次数。
图 15. OCC-BOTDA系统原理示意图。 (a)时域关系; (b)频域关系
Fig. 15. Schematic of OCC-BOTDA system. (a) Time-domain relationship; (b) frequency-domain relationship
图 16. OCC-BOTDA系统对BGS的测量。 (a)单发抽运脉冲光即可测量分布式的BGS; (b)不同应变的改变量时的BGS
Fig. 16. Measurement of BGS via OCC-BOTDA system. (a) Distributed BGS measured by single shot of pump pulse wave; (b) various BGS for different strain changes
OCC-BOTDA系统对BGS的测量如
OCC-BOTDA系统对机械撞击的实验测量结果如
图 17. OCC-BOTDA系统对机械撞击的实验测量结果。(a)经过二次互相关算法处理后的BGS演进图; (b)解调出的应变振动波形
Fig. 17. Experimental results of mechanical shocks via OCC-BOTDA system. (a) Evolution of BGS processed by twice correlation algorithm; (b) demodulated strain vibration waveforms
为了显示OCC-BOTDA系统的快速采集特性,本课题组通过将相邻OCC的扫频频率差设置为20 MHz来模拟一个快速的应变过程。光纤长度为10 m,应变振动采样率为6.25 MHz。OCC-BOTDA系统对快速开关事件的实验测量结果如
图 18. OCC-BOTDA系统对快速开关事件的实验测量结果:算法处理后的BGS演进图。 (a)平均200次; (b)未加平均; (c)解调出的BFS的改变量与时间的关系
Fig. 18. Experimental results of the switch event via OCC-BOTDA system: evolution of BGS processed by twice correlation algorithm for various conditions. (a) 200 times averaging; (b) no averaging; (c) demodulated BFS changes over time
3.5 基于光学频率梳的动态布里渊光纤传感系统
针对限制因素4,也可以在频域直接产生多条光谱线来实现单发测量,该系统称为自扫频BOTDA系统(SF-BOTDA)。SF-BOTDA系统的演变与原理示意图如
图 19. SF-BOTDA系统的演变与原理示意图。 (a)基于光学频率梳的BOTDA系统;(b)SF-BOTDA系统
Fig. 19. Schematic of SF-BOTDA system. (a) BOTDA system based on optical frequency comb; (b) SF-BOTDA system
最近,Fang等[51]提出了一种基于新型的单发BOTDA快速测量系统技术方案。经过正交频分复用技术(OFDM)调制产生的探测光,用来采集分布式的BGS,同时采用相干探测技术增强SNR。该技术方案既不需要扫频,又不需要扰偏,更不需要平均。单发BOTDA系统原理示意图如
图 20. 单发BOTDA系统原理示意图。 (a)频域中的OFDM信号的子载波; (b)双偏振双边带OFDM探测光;(c)抽运光和探测光的相互作用及布里渊谱型解调过程;(d)SBS过程中,BGS和BLS的示意图
Fig. 20. Schematic of single-shot BOTDA system. (a) Subcarriers of OFDM signal in frequency domain; (b) dual-polarized double-sideband OFDM probe; (c) pump and probe interaction and Brillouin spectrum extraction; (d) BGS and BLS after SBS process
图 21. 根据静态温度的测量实验结果重构的布里渊谱。 (a)x偏振态; (b)y偏振态; (c)双偏振相结合的情况; (d)在(c)中A、B、C三点的测量数据及其洛伦兹拟合曲线
Fig. 21. Reconstructed Brillouin spectrograms based on the measurement experimental results of static temperature. (a) x-polarization; (b) y-polarization; (c) combined dual-polarization; (d) measured data and Lorentzian fitting curves of marker A, B, C in (c)
4 结论
与传统的电学传感器相比,光纤传感具有很突出的优势,尤其是BOTDA系统。在现代工业中,对动态分布式应变测量的需求在不断增强,然而传统的BOTDA系统受到光纤长度、平均次数、光波频率切换时间和扫频个数等因素的限制,测量时间需要数秒到数分钟,因此被认为只适合于静态或者缓慢变化应变的测量。
从BOTDA的原理出发,总结了近期利用BOTDA系统在快速分布式应变测量方面的研究进展。光纤既是传感器又是传光介质,光纤长度对应变振动采样率的限制是固有的。对于传感光纤,使用保偏光纤时,SNR较高,平均次数少,但价格昂贵,不适宜长距离传感;使用普通单模光纤时,偏振补偿技术可以取代扰偏器,有效地消除偏振噪声,并减少平均次数,但该方法需要分别测量两个偏振方向的布里渊信号,因而偏振态切换时间、频率切换时间和扫频个数同样限制应变振动采样率。对于光学FAT,利用高性能的AWG产生捷变频信号,使得光波频率切换时间可以忽略不计,然而扫频个数仍然限制着应变振动采样率。对于单斜坡法,应变振动采样率只受限于光纤长度,然而应变动态范围受限于BGS的斜坡范围;对于多斜坡法,虽然增加了动态范围,但是应变振动采样率和斜坡个数成反比。值得注意的是,与BGS和BPSS两种传感曲线相比,KS传感曲线呈现出更加优异的性能:更大的动态范围和免疫于抽运脉冲光的功率波动。对于OCC技术,应变振动采样率只受限于光纤长度,而且系统的空间分辨率和动态范围可以通过啁啾信号预编写来实现智能化(或个性化)配置。对于光学频率梳技术,虽然应变振动采样率也只受限于光纤长度,但是其空间分辨率受限于频率梳间隔和FFT的时间窗,这极大地限制了该方法的实际应用价值。随着限制因素的逐个突破,结合多个新技术搭建而成的快速BOTDA系统会具有更好的性能。例如,结合偏振补偿技术、KS传感曲线和OCC在单模光纤的长距离快速应变测量中具有潜在的优势。另外,在保持快速应变测量的同时,如何有效地提高系统的空间分辨率也将会成为重要的研究方向。
综上所述,快速BOTDA系统受到国内外的广泛研究,多种改进方案的实用化在**边境安全、现代工业控制、土木机构监测、灾害监测、航空航天器安全监测和海洋探测等领域具有广阔的应用前景。
[5] HotateK. Recent achievements in BOCDA/BOCDR[C]//Proceedings of 2014 IEEE Sensors, 2014: 142- 145.
HotateK. Recent achievements in BOCDA/BOCDR[C]//Proceedings of 2014 IEEE Sensors, 2014: 142- 145.
[6] Boyd RW. Nonlinear optics[M]. 3rd ed. Pittsburgh: Academic Press, 2008: 429- 471.
Boyd RW. Nonlinear optics[M]. 3rd ed. Pittsburgh: Academic Press, 2008: 429- 471.
[18] YamauchiT, HotateK. Distributed and dynamic strain measurement by BOCDA with time-division pump-probe generation scheme[C]. Conference on Lasers and Electro-Optics, 2004: CWA57.
YamauchiT, HotateK. Distributed and dynamic strain measurement by BOCDA with time-division pump-probe generation scheme[C]. Conference on Lasers and Electro-Optics, 2004: CWA57.
[21] Ong S SL, HotateK. Dynamic strain measurement at 50 Hz using a Brillouin optical correlation domain analysis based on fiber optic sensor[C]. 5 th Pacific Rim Conference on Lasers and Electro-Optics , 2003: 7993564.
Ong S SL, HotateK. Dynamic strain measurement at 50 Hz using a Brillouin optical correlation domain analysis based on fiber optic sensor[C]. 5 th Pacific Rim Conference on Lasers and Electro-Optics , 2003: 7993564.
[30] Agrawal GP. Nonlinear fiber optics[M]. 4th ed. Amsterdam: Elsevier Academic Press, 2007.
Agrawal GP. Nonlinear fiber optics[M]. 4th ed. Amsterdam: Elsevier Academic Press, 2007.
[31] Dong YK, Zhang HY, Zhou DP, et al. Characterization of Brillouin gratings in optical fibers and their applications[M]. London: Intech Publisher, 2012: 115- 136.
Dong YK, Zhang HY, Zhou DP, et al. Characterization of Brillouin gratings in optical fibers and their applications[M]. London: Intech Publisher, 2012: 115- 136.
[35] 王本章. 基于光学捷变频的动态分布式布里渊光纤传感技术研究[D]. 哈尔滨: 哈尔滨工业大学, 2016.
王本章. 基于光学捷变频的动态分布式布里渊光纤传感技术研究[D]. 哈尔滨: 哈尔滨工业大学, 2016.
Wang BZ. Dynamic distributed Brillouin optical fiber sensing based on optical frequency-agile technology[D]. Harbin: Harbin Institute of Technology, 2016.
Wang BZ. Dynamic distributed Brillouin optical fiber sensing based on optical frequency-agile technology[D]. Harbin: Harbin Institute of Technology, 2016.
[36] López-GilA, Domínguez-LópezA, Martín-LópezS, et al. Simple method for the elimination of polarization noise in BOTDA using balanced detection of orthogonally polarized Stokes and anti-Stokes probe sidebands[C]. 23 rd International Conference on Optical Fiber Sensors , 2014, 9157: 91573U.
López-GilA, Domínguez-LópezA, Martín-LópezS, et al. Simple method for the elimination of polarization noise in BOTDA using balanced detection of orthogonally polarized Stokes and anti-Stokes probe sidebands[C]. 23 rd International Conference on Optical Fiber Sensors , 2014, 9157: 91573U.
[39] UrricelquiJ, López-FernandinoF, SaguesM, et al. Polarization diversity for Brillouin distributed fiber sensors based on a double orthogonal pump[C]. 23 rd International Conference on Optical Fiber Sensors , 2014, 9157: 91576A.
UrricelquiJ, López-FernandinoF, SaguesM, et al. Polarization diversity for Brillouin distributed fiber sensors based on a double orthogonal pump[C]. 23 rd International Conference on Optical Fiber Sensors , 2014, 9157: 91576A.
[47] Zhou D W, Dong Y K, Wang B Z, et al. Single-shot BOTDA based on optical chirp chain probe wave for distributed ultra-fast measurement[J]. Light: Science & Applications, 2017.
Zhou D W, Dong Y K, Wang B Z, et al. Single-shot BOTDA based on optical chirp chain probe wave for distributed ultra-fast measurement[J]. Light: Science & Applications, 2017.
Article Outline
周登望, 王本章, 巴德欣, 徐金龙, 徐鹏柏, 姜桃飞, 张东昱, 李惠, 董永康. 用于动态应变测量的快速分布式布里渊光纤传感[J]. 光学学报, 2018, 38(3): 0328005. Zhou Dengwang, Wang Benzhang, Ba Dexin, Xu Jinlong, Xu Pengbai, Jiang Taofei, Zhang Dongyu, Li Hui, Dong Yongkang. Fast Distributed Brillouin Optical Fiber Sensing for Dynamic Strain Measurement[J]. Acta Optica Sinica, 2018, 38(3): 0328005.