Monitoring and experimental analysis of distributed optical fiber pipeline leakage based on wavelet packet energy spectrum
-
摘要: 基于Mach-Zehnder/Sagnac混合干涉仪原理，利用分布式光纤对城市燃气管道进行泄漏监测模拟实验研究。对测量光纤长4 km和6 km时的干涉信号利用小波包能量谱提取法和绝对距离法相结合的方法诊断管道泄漏与否，再结合频谱图进行泄漏定位，最后对系统的虚警率(false alarm rate，FAR)进行了分析。结果表明，与传统直接通过频谱图识别零点频率的方法相比，虚警率降低了8.475%，能够更准确地识别燃气管道泄漏与否，从而提高系统监测及定位的可靠性。Abstract: Based on the principle of Mach-Zehnder and Sagnac hybrid interferometer, the simulated experiment of urban gas pipeline leakage monitoring was carried out using distributed optical fiber sensor. Firstly, the method of wavelet packet energy spectrum extraction and the absolute distance method were combined to process the interference signals caused by pipeline leakage when measuring optic fiber length are 4 km and 6 km. Then the gas pipeline leakage was located according with the frequency spectrograms. Finally, the false alarm rate(FAR) of the distributed optical fiber system was analyzed.The results show that the FAR of proposed method reduces by 8.475%, compared with the traditional zero frequency spectrum localization method, it can more accurately identify the gas pipeline leak or not, thereby enhancing the reliability of system monitoring and positioning.
-
引言
随着高能高功率激光器研制及其应用的不断发展，光学元件损伤几率越来越高。激光损伤阈值是评价光学元件抗损伤能力的重要技术指标，是影响高能高功率激光输出功率和能量水平的关键因素。脉冲激光导致的损伤通常是由于缺陷、膜层应力、表面污染物或者光吸收导致的热效应引起。在一定条件下，光学元件损伤阈值测量与材料性能和激光参数，特别是与激光波长、光斑尺寸和脉冲宽度等有关，影响测量结果的因素很多[1-2]。
国际上对激光损伤阈值测量装置及其测量方法进行了大量研究[3-8]。国内相关研究单位，也参考ISO标准和国内外同行的研究成果，研制了一系列激光损伤阈值测量装置[9-14]，在激光损伤阈值的测量上进行了有益的探索。但总的来说，研制的测量装置的工程化和标准化程度不高，测量设备条件不同，测试程序不统一，人为因素、环境因素和器件因素造成测量值存在较大偏差，测量数据的复现性和可比性比较差。
为了开展光学元件损伤特性及测试方法研究，我们研制了一套1 064 nm、532 nm激光损伤测试装置，该装置由计算机集成控制，具有良好操作性和运行参数控制能力，可实现1-on-1,S-on-1, R-on-1三种测试方法下的光学元件激光损伤测试。
1 激光损伤阈值测量方法
光学元件激光损伤阈值是指在激光损伤试验中，统计分析得到的光学元件表面最大不损伤的激光能量密度或峰值功率密度值，简称为零概率损伤阈值。目前国际上对于激光损伤阈值的测量，基本遵循ISO21254标准的要求。同时为适应我国的技术条件，在ISO21254标准基础上，我国也制定了相应的GB∕T 16601-2017[15]。上述标准中提出的光学元件脉冲激光损伤阈值测量方法有3种：1-on-1、S-on-1和R-on-1。
1.1 1-on-1测量方法
1） 采用不同能量密度的脉冲激光对样品进行辐照，每个能量密度辐照一组测试点，每组测试点一般不少于10个，每个测试点只照射1次。测试点之间的距离至少是入射光斑直径的4倍，以避免测试点之间的相互重叠的干扰。激光损伤探测判别系统得到不同能量密度下测试点的损伤几率。
2） 以能量密度为横坐标，损伤几率为纵坐标，对损伤几率-能量密度进行线型拟合，拟合得到的直线与横坐标的交点处的能量密度就是零概率损伤阈值，1-on-1的测量原理如图1所示。
1.2 S-on-1测量方法
1） 采用不同能量密度的脉冲激光对样品进行辐照，每个能量密度辐照1组测试点，每组测试点一般不少于10个，每个测试点照射指定的脉冲数（比如1 000次，如果发生损伤则停止照射，移至下一个测试点）。测试点之间的距离与1-on-1的规定相同。激光损伤探测判定系统得到不同能量密度下测试点的损伤几率。
2） 和1-on-1测试方法相同，拟合得到零概率损伤阈值，S-on-1的测量原理如图2所示。
1.3 R-on-1测量方法
1） 采用不同能量密度的脉冲激光对样品进行辐照，以能量密度递增的方式辐照1个测试点，如果发现该测试点损伤， 移至下一个测试点。测试点之间的距离与1-on-1的规定相同。激光损伤探测判定系统得到不同能量密度下测试点的损伤几率。
2） 在数据处理时通常采用两种方法，一种和1-on-1测试方法相同，将数据点按照能量密度分类统计损伤概率后拟合得到零概率损伤阈值，一种是计算每个位置刚损伤时的能量密度和前一个未损伤能量密度的平均值。第二种方法可得到较多的数据点，实际工作中采用较多。该测量方法如图3所示。
2 激光损伤阈值自动测量装置
双波长光学元件激光损伤阈值自动测量装置，原理及实物如图4和图5所示。该装置主要包括：脉冲激光光源系统、激光光束诊断系统、激光损伤探测系统以及控制与处理系统组成。测试过程中，将样品测试点放置于激光光路中，计算机控制脉冲激光光源系统产生不同能量密度的脉冲激光，使用不同能量密度或功率密度的单个激光脉冲对其进行辐照。对于每个选定的脉冲能量（或功率）至少选取10个测试点进行辐照并记录。使用光束诊断系统记录下每个测试点实际使用的脉冲能量密度，同时记录下样品经激光辐照后的损伤情况（损伤或未损伤）。主控软件根据获得的辐照数据，构建一幅损伤概率随能量密度变化的关系图。损伤概率数据线性外推到零损伤概率即可得到损伤阈值。
2.1 脉冲激光光源系统
激光器采用北京镭宝公司的泵浦Nd:YAG激光器，输出基频波长为1.064 μm，通过倍频晶体，输出532 nm波长激光。输出模式为TEM00，调Q后脉宽为10 ns，最大脉冲能量为2 J（1.064 μm），工作频率1 Hz～10 Hz可调。
损伤阈值测试中的可变衰减器大多采用偏振衰减的方式，由1/2波片和薄膜偏振片（或偏振分光棱镜）组成，可变衰减器衰减范围可实现0.1%～95%。损伤阈值测试中，通常要对激光能量实现大范围衰减，衰减范围要达到103量级，即最小透过率一般要小于0.1%。因此，可变衰减器最小透过率大小对于衰减范围至关重要，主要影响因素有3个方面：1） 是入射激光偏振度；2） 是偏振器消光比；3） 是半波片位相延迟偏差。对于大能量激光器来说，存在热退偏效应，常常在激光器和可变衰减器之间增加了一个薄膜偏振片，增加激光的线偏振度。通过实验验证，本装置光源系统和可变衰减器组合后可在待测元件表面获得能量密度0.1 J/cm2～100 J/cm2的脉冲激光。
2.2 光束诊断系统
激光光束诊断系统由楔形分束镜、能量计、快速响应光电探测器、示波器、光斑空域分布测量系统等组成，以获得辐照脉冲激光的能量、光斑大小和脉宽，如图6所示。楔形分束镜1的分束比可以通过能量计进行测量。损伤阈值测量过程中，激光光束诊断系统的能量计测量得到的脉冲激光的能量值通过USB接口传输到数据采集和处理系统，计算机的主控软件将分束后的反射光的能量乘以分束比，就可以计算出辐照到待测样品测量点的脉冲激光能量。通过光路调校，光斑分析仪CCD靶面与样品的待测面对称共轭，在多路同步延时电路的控制下，光斑分析仪采集得到的脉冲激光光斑的大小就是辐照到待测样品测量点的光斑大小。时域分布测量系统由快速光电探测器和1 G带宽的示波器组成，可以满足ns级短脉冲激光脉宽精密测量的要求。
2.3 在线损伤探测系统
在线损伤探测方法主要有2种：1） 是散射光探测；2） 是在线显微成像。散射光探测可以直接探测测试激光的散射信号，也可以探测与测试激光叠加照射在测试点的另一独立激光束的散射信号。基于分离光源的散射光测量装置，需要一个指向稳定性很高且强度波动很小的激光器作为辐射光源。在线显微成像方法可以在辐照过程中直接检查样品表面状态，通过比较辐照前后样品表面图像判断损伤，具有可靠性高、直观且适合自动测量的优点，在损伤阈值自动测试领域被广泛采用。本文中装置通过照明光源与在线显微相机结合进行损伤自动判断，如图7所示。
2.4 控制与处理系统
控制与处理系统是损伤阈值自动测量装置的核心，由计算机和控制箱组成。计算机通过图8所示测试软件对激光器快门、可变衰减器衰减比例、待测元件夹持机构运动等进行控制，同时控制能量计、光斑分布测量相机、在线显微相机等测量组件的数据采集，并进行分析和计算处理。
由于短脉冲激光的脉宽只有10 ns左右，要准确捕获光斑分布和脉冲波形，需要保证激光器出光时间与触发光斑分布测量相机和示波器触发信号在时序上精确匹配。本装置中所采用激光器调Q信号与出光时间为1 μs，光斑测量相机外触发信号与实际曝光之间延迟43 μs，若采用调Q信号去触发光斑测量相机，则无法捕获到脉冲光信号，因此本装置通过自研的多路同步延迟电路同步触发激光器输出、光斑采集和脉冲波形的采集。
3 测量结果及分析
3.1 测量结果
实验样品基片为K9玻璃，单面镀制1 064 nm铝反射膜和增透膜，样品尺寸50 mm×50 mm，入射角3º，光束有效面积0.001 cm2，有效脉冲宽度9.5 ns，每个能量密度下辐照20个点，能量密度台阶数设为15。激光器能量稳定性优于1%，可变衰减器范围从0.1%～97%。图9(a)为1 064 nm增透膜的损伤概率曲线，损伤阈值为27.09 J/cm2，图9(b)为1 064 nm铝反射膜的损伤概率曲线，损伤阈值为3.21 J/cm2。图10为测试过程中在线观测到的反射膜层的损伤形貌。如图10所示，铝反射膜损伤形貌主要表现为孔洞和膜层剥离，增透膜损伤主要集中在后表面，小孔洞聚集。
![]()
图 10 1 064 nm激光辐照下铝反射膜和增透膜损伤形貌Figure 10. Damage morphology of aluminum reflection films and antireflection films under 1 064 nm laser irradiation采用1 064 nm和532 nm激光，分别获得了1-on-1、S-on-1、R-on-1方式下对增透膜和铝反射膜零概率损伤阈值的测试结果，如表1所示。
表 1 损伤阈值测量结果Table 1. Measurement results of LIDT波长/nm 样品名称 LIDT/(J/cm2) 1 on 1 100 on 1 R on 1 1 064 铝反射膜 3.21 —— 2.97 增透膜 27.09 —— 25.93 532 铝反射膜 0.56 —— 0.32 增透膜 6.42 4.55 6.38 3.2 测量结果不确定度分析
激光损伤阈值的测量结果受到能量密度测试过程的影响。同时, 用能量密度区间内有限次数测量的损伤几率代替平均能量密度处的损伤几率，来确定能量密度的几率数据点的计算以及对这些点采用线性拟合的方式，也影响到测量不确定度。根据GBT 16601.2-2017以及ISO 21254-2中测试结果评估流程，对于每一个能量密度水平Hi，观测到的损伤概率Pi，由式（1）算出：
$$ {P_i} = \frac{{n_i^{\rm{d}}}}{{n_i^{\rm{d}} + n_i^{{\rm{nd}}}}} $$ (1) 式中：nid+nind是产生nid个测试点损伤而辐照的测试点总数。
用式（2）或式（3）估算数据集中每个损伤概率Pi的不确定度σi；
$$ {\sigma _i} = \sqrt {\frac{{n_i^{{\rm{nd}}}}}{{n_i^{\rm{d}}\left( {n_i^{\rm{d}} + n_i^{{\rm{nd}}}} \right)}} + \varepsilon _f^2}, \;\;\;\;\;当{n^{\rm{d}}} > 0 $$ (2) $$ {\sigma _i} = \sqrt {\varepsilon _f^2}, \;\;\;\;当{n^{\rm{d}}} = 0 $$ (3) 式中：εf为所测能量密度的不确定度分量。能量密度由能量、分束比、光斑有效面积共同决定，其不确定度分析结果如表2所示。
表 2 能量密度不确定度分析Table 2. Uncertainty analysis of laser energy density测量不确定度 激光能量密度/% 激光能量计测量相对不确定度(u1) 1.25 分光比测量相对不确定度(u2) 0.18 光斑有效面积测量[16](u3) 1.1 能量密度计算相对合成不确定度(u4) 1.67 获得不同能量密度Hi下的损伤概率Pi数据后，对测量的损伤概率曲线进行线性拟合，可以得到的斜率m和截距b：
$$ m = \frac{1}{\varDelta }\left[ {\left( {\mathop \sum \nolimits_i \frac{1}{{\sigma _i^2}}} \right)\left( {\mathop \sum \nolimits_i \frac{{{P_i}{H_i}}}{{\sigma _i^2}}} \right) - \left( {\mathop \sum \nolimits_i \frac{{{H_i}}}{{\sigma _i^2}}} \right)\left( {\mathop \sum \nolimits_i \frac{{{P_i}}}{{\sigma _i^2}}} \right)} \right] $$ (4) $$b = \frac{1}{\varDelta }\left[ {\left( {\mathop \sum \nolimits_i \frac{{H_i^2}}{{\sigma _i^2}}} \right)\left( {\mathop \sum \nolimits_i \frac{{{P_i}}}{{\sigma _i^2}}} \right) - \left( {\mathop \sum \nolimits_i \frac{{{H_i}}}{{\sigma _i^2}}} \right)\left( {\mathop \sum \nolimits_i \frac{{{P_i}{H_i}}}{{\sigma _i^2}}} \right)} \right] $$ (5) 其中
$$\varDelta = \left( {\mathop \sum \nolimits_i \frac{1}{{\sigma _i^2}}} \right)\left( {\mathop \sum \nolimits_i \frac{{H_i^2}}{{\sigma _i^2}}} \right) - {\left( {\mathop \sum \nolimits_i \frac{{{H_i}}}{{\sigma _i^2}}} \right)^2}$$ (6) $$ {H_{{\rm{th}}}} = - \frac{b}{m} $$ (7) 同时利用最小二乘法通过式（4）和式（5）可得到拟合线性方程P=mH+b或者H=P/m−b/m，根据零概率损伤阈值定义，用式（7）确定零概率损伤阈值Hth。以上节中1 064 nm波长下增透膜和铝反射膜的1 on 1测试结果为例，计算前面图9（a）增透膜损伤概率曲线对应的拟合方程为P=0.022Hth−0.596，损伤阈值Hth为27.09 J/cm2，图9（b）铝反射膜损伤概率曲线对应的拟合方程为P=0.162Hth−0.519 7，损伤阈值Hth为3.21 J/cm2。
根据标准不确定度的A类评定方法[17]，H的标准不确定度为它的估计标准偏差，同时根据标准偏差传递规律，损伤阈值的不确定度σth为
$$ {\sigma _{{\rm{th}}}} = - \frac{1}{{{m^2}}}\sqrt {{b^2}\sigma _m^2 + {m^2}\sigma _b^2} $$ (8) 式中：
$$ {\sigma _b} = \sqrt {\frac{1}{\varDelta } \sum \frac{{H_i^2}}{{\sigma _i^2}}} $$ (9) $$ {\sigma _m} = \sqrt {\frac{1}{\varDelta } \sum \frac{1}{{\sigma _i^2}}} $$ (10) 通过式（10），上节中1 064 nm增透膜和铝反射膜损伤阈值的不确定度σth分别为1.06 J/cm2、0.18 J/cm2，相对不确定度分别为3.91%和5.61%。不确定度大小主要与数据拟合直线的斜率大小以及数据拟合引起的标准偏差差异相关。由此可见，光学元件激光损伤阈值测量过程中，能量密度带来的不确定影响较小，不确定度主要来源于损伤几率及损伤概率数据拟合带来的影响。若采用线性拟合，数据拟合不确定度与光学材料损伤几率测量数据斜率及线性程度有很大关系。
4 结论
双波长激光损伤阈值自动测量装置通过计算机集成控制，利用可变衰减器、同步控制电路和分束共轭方法对调Q脉冲激光（10 ns）不同能量密度光斑进行了准确捕捉、测量。装置测量程序自动化程度高，可实现1-on-1、r-on-1、s-on-1三种方式的激光损伤阈值自动测量，为开展光学元件的激光损伤研究提供良好的损伤实验平台。
-
![]()
图 2 分布式光纤燃气管道泄漏检测装置
Figure 2. Distributed optical fiber pipeline leak detection device
![]()
图 3 测量光纤4 km、6 km无泄漏信号频谱图
Figure 3. Signal frequency spectrograms when measuring optical fiber of 4 km and 6 km without leakage
![]()
图 4 管道不同气压和泄漏状态的小波包能量谱图。
Figure 4. Wavelet packet energy spectrograms of pipeline with different pressures and leakage states
表 1 检测信号小波包分解的频带划分
Table 1 Frequency band division of wavelet packet decomposition
序号 1 2 3 4 5 6 7 8 频谱范围/kHz 0~12.5 12.5~25 25~37.5 37.5~50 50~62.5 62.5~75 75~87.5 87.5~100 
下载: 导出CSV
-
[1] 任峰, 何仁洋, 刘长征, 等.城市燃气埋地钢质管道泄漏检测定位技术[J].管道技术与设备, 2016(2):13-16. doi: 10.3969/j.issn.1004-9614.2016.02.005 REN Feng, HE Renyang, LIU Changzheng, et al. Leakage detection and positioning technology of urban gas buried steel pipeline[J]. Pipeline Technique and Equipment, 2016 (2): 13-16. doi: 10.3969/j.issn.1004-9614.2016.02.005
[2] 张晓威, 刘锦昆, 陈同彦, 等.基于分布式光纤传感器的管道泄漏监测试验研究[J].水利与建筑工程学报, 2016, 14(03):1-6. doi: 10.3969/j.issn.1672-1144.2016.03.001 ZHANG Xiaowei, LIU Jinkun, CHEN Tongyan, et al. Experimental study on pipeline leak monitoring based on distributed optical fiber sensor[J]. Journal of Water Conservancy and Building Engineering, 2016, 14 (03): 1-6. doi: 10.3969/j.issn.1672-1144.2016.03.001
[3] SPAMMER S J, SWART P L, BOOYSEN A. Interferometric distributed optical-fiber sensor[J]. Applied Optics, 1996, 35(22):4522-4525. doi: 10.1364/AO.35.004522
[4] HUANG S C, LIN W W, TSAI M T, et al. Fiber optic in-line distributed sensor for detection and localization of the pipeline leaks[J]. Sensors and Actuators A: Physical, 2007, 135(2): 570-579. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=5b273f7554f28674e439361e847f44f8
[5] 胡正松, 杨其华, 乔波.干涉分布式光纤水下长输气管道泄漏检测系统设计[J].激光与光电子学进展, 2012, 49(7):69-73. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=jgygdzxjz201207010 HU Zhengsong, YANG Qihua, QIAO Bo. Design of interference distributed fiber-optic underwater long gas pipeline leakage detection system[J]. Laser and Optoelectronics Progress, 2012, 49(7): 69-73. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=jgygdzxjz201207010
[6] HUANG Yue, WANG Qiang, SHI Lilian, et al. Underwater gas pipeline leakage source localization by distributed fiber-optic sensing based on particle swarm optimization tuning of the support vector machine[J]. Applied Optics, 2016, 55(2):242. doi: 10.1364/AO.55.000242
[7] 崔晓燕.城市燃气管道泄漏检测技术的发展现状探究[J].中国石油石化, 2016(22):67-68. http://www.cnki.com.cn/Article/CJFDTotal-SHZG201622036.htm CUI Xiaoyan. Research on the development of leakage detection technology for urban gas pipeline[J]. China Petroleum and Petrochemical, 2016 (22):67-68. http://www.cnki.com.cn/Article/CJFDTotal-SHZG201622036.htm
[8] 黄悦, 王强, 韩玲娟, 等.基于支持向量机的分布式光纤泄漏检测定位系统及实验分析[J].应用光学, 2015, 36(3):424-429. http://d.old.wanfangdata.com.cn/Periodical/yygx201503016 HUANG Yue, WANG Qiang, HAN Lingjuan, et al. Distributed optical fiber leak detection positioning system based on support vector machines and experimental analysis[J]. Journal of Applied Optics, 2015, 36 (3): 424-429. http://d.old.wanfangdata.com.cn/Periodical/yygx201503016
[9] 章仁杰, 杨其华, 王强, 等.混合干涉型分布式光纤天然气管道泄漏检测及定位方法研究[J].中国计量学院学报, 2013, 24(02):114-118. doi: 10.3969/j.issn.1004-1540.2013.02.003 ZHANG Renjie, YANG Qihua, WANG Qiang, et al. Study on leak detection and location method of hybrid interferometric distributed optical fiber gas pipeline[J]. Journal of China Jiliang University, 2013, 24 (02): 114-118. doi: 10.3969/j.issn.1004-1540.2013.02.003
[10] 韩玲娟, 王强, 杨其华, 等.基于分布式光纤传感的水下输气管道泄漏检测与定位分析[J].传感技术学报, 2015(7):1097-1102. doi: 10.3969/j.issn.1004-1699.2015.07.026 HAN Lingjuan, WANG Qiang, YANG Qihua, et al. Leak detection and location analysis of underwater gas pipeline based on distributed optical fiber sensing[J]. Journal of Sensing Technology, 2015 (7): 1097-1102. doi: 10.3969/j.issn.1004-1699.2015.07.026
[11] 谢娜, 崔广宇.测试性虚警率指标分配技术研究[J].测控技术, 2014, 33(4):141-144. doi: 10.3969/j.issn.1000-8829.2014.04.038 XIE na, CUI Guangyu. Research on index allocation of testability false alarm rate[J]. Measurement and Control Technology, 2014, 33 (4): 141-144. doi: 10.3969/j.issn.1000-8829.2014.04.038
[12] 王超, 黄彩萍, 艾德米, 等.基于统计小波包能量分布的结构健康状况预警[J].土木工程与管理学报, 2016, 33(06):12-15. doi: 10.3969/j.issn.2095-0985.2016.06.003 WANG Chao, HUANG Caiping, AI demi, et al. Prediction of structural health status based on statistical wavelet packet energy distribution[J]. Journal of Civil Engineering and Management, 2016, 33 (06): 12-15. doi: 10.3969/j.issn.2095-0985.2016.06.003
[13] 王玮.基于聚类分析的建筑施工安全事故统计分析[J].蚌埠学院学报, 2017, 6(01):34-37. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=bbxyxb201701009 WANG Wei. Statistical analysis of construction safety accidents based on cluster analysis[J]. Journal of Bengbu College, 2017, 6 (01): 34-37. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=bbxyxb201701009
[14] 廖志芳, 李鹏, 刘克准, 等.数据聚类分析新方法研究[J].计算机工程与应用, 2009, 45(10):147-150. doi: 10.3778/j.issn.1002-8331.2009.10.044 LIAO Zhifang, LI Peng, LIU Kezhun, et al. Research on new methods of data clustering analysis[J]. Computer Engineering and Applications, 2009, 45 (10): 147-150. doi: 10.3778/j.issn.1002-8331.2009.10.044
-
期刊类型引用(1)
1. 刘宝怡，张强，蒋世磊，孙国斌，杨国锋. 基于薄膜损伤阈值测试的聚焦光学系统设计. 光学与光电技术. 2025(01): 9-14 . 
百度学术
其他类型引用(1)
计量
- 文章访问数: 762
- HTML全文浏览量: 174
- PDF下载量: 91
- 被引次数: 2
陕公网安备 61011302001501号 