CN107664571A - 一种光时域反射仪曲线数据事件检测的方法 - Google Patents
一种光时域反射仪曲线数据事件检测的方法 Download PDFInfo
- Publication number
- CN107664571A CN107664571A CN201610611889.6A CN201610611889A CN107664571A CN 107664571 A CN107664571 A CN 107664571A CN 201610611889 A CN201610611889 A CN 201610611889A CN 107664571 A CN107664571 A CN 107664571A
- Authority
- CN
- China
- Prior art keywords
- high frequency
- frequency coefficient
- wavelet
- amplitude
- wavelet decomposition
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
- G01M11/30—Testing of optical devices, constituted by fibre optics or optical waveguides
- G01M11/31—Testing of optical devices, constituted by fibre optics or optical waveguides with a light emitter and a light receiver being disposed at the same side of a fibre or waveguide end-face, e.g. reflectometers
- G01M11/3109—Reflectometers detecting the back-scattered light in the time-domain, e.g. OTDR
- G01M11/3145—Details of the optoelectronics or data analysis
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Optics & Photonics (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明揭示了一种光时域反射仪(OTDR)曲线数据事件检测的方法。小波分析对信号进行多尺度分析,可以聚焦到信号的任意细节,通过不同层级的小波分解,信号分解为不同层级的低频和高频部分,设置阈值将幅度大的高频系数保留出来,基于这些高频系数的位置,进行事件类型识别。本发明揭示的方法,能有效检测OTDR曲线数据中的事件。
Description
技术领域
本发明涉及光时域反射仪(OTDR)技术领域,有关曲线数据事件检测的方法。
背景技术
光时域反射仪(OTDR)是光在光纤传输过程中的瑞利散射和菲涅尔反射所产生的背向散射而制成的光电一体化仪表,OTDR通过向光纤中发射光脉冲,然后从OTDR接收返回的信号。当光脉冲在光纤内传输时,在光纤前端、尾端、光纤之间的连接处、光纤弯曲处、光纤的破损处产生散射、反射,其中小部分的散射、反射返回到OTDR中,通过OTDR中的信号采集器,获取到光纤内不同时间片散射、反射信息,基于散射、反射信号和返回信号之间的时间差,再确定光在光纤中的传播速度,即能求出采集到的数据对应于光纤中的距离位置。目前,OTDR作为一种非破坏性的光纤测量仪,在通信光纤的施工、维护、故障检测等方面得到广泛的应用。它可测定光纤首段、尾端、反射和非反射(衰减)等事件,也能测定光纤损耗、光纤长度等信息,是一种多功能的光纤测量仪器。
发明内容
本发明提供一种对OTDR曲线数据事件检测的方法。用时域频率域转换模块进行事件检测,转换模块中的小波分析对信号进行多尺度分析,可以聚焦到信号的任意细节,通过不同层级的小波分解,信号分解为不同层级的低频和高频系数,设置阈值将幅度大的高频系数保留出来,这些高频系数即为事件所在位置,再进行事件类型识别。本发明揭示的方法,能有效检测OTDR曲线数据中的事件。
为了实现上述目的,本发明提出如下技术方案:一种OTDR曲线数据事件检测的方法,其内容包括如下步骤:
S1,OTDR数据用小波变换进行分解,分解为不同层级的低频和高频系数;
S2,对分解所得的高频系数,设置阈值,分离出幅度大的高频系数;
S3,根据分离出的高频系数对应位置,进行事件类型识别。
所述S1中用小波变换,使用Haar小波基,
Haar尺度函数定义为: ,
Haar小波基函数定义为: 。
所述S1中用小波变换,使用Daubechies小波基,
Daubechies尺度函数定义为:,
Daubechies小波基函数定义为:,
,其中:和为滤波器系数,只有2M个非零项。
所述S1中用不同层级小波分解,小波分解层级根据一个脉宽时间内采集到的离散点数量决定,小波分解层级计算公式:,w表示脉宽,s表示采样率,L表示小波分解层级。
所述S2中分离出幅度大的高频系数,小波分解后的任一层级,存在大于该阈值的高频系数时,将这些高频系数及所在层级信息保留出来。
所述S2中分离出幅度大的高频系数,小波分解后的任一层级,不存在大于该阈值的高频系数,将所有层级的高频系数进行叠加,找出大于该阈值的叠加值。
所述S3中分离出的幅度大的高频系数对应位置,幅度大的高频系数来自小波分解的不同层级,再进行事件识别。
所述S3中分离出的幅度大的高频系数对应位置,幅度大的高频系数来自小波分解后不同层级的累加,再进行事件识别。
本发明提出了一种OTDR曲线数据事件检测方法的装置,装置中的时域频率域转换模块,OTDR数据能进行时域和频率域相互转换,将时域中的OTDR事件转换到频率域进行检测。OTDR数据通过小波分解为多层级的低频和高频系数,可以聚焦到信号的任意细节,能有效对弱信号事件进行检测,设置阈值将幅度大的高频系数提取出来,这些高频系数即为事件所在位置,再进行事件类型识别。
与现有技术相比,本发明揭示的方法,能有效检测OTDR曲线数据中的事件。
附图说明
图1是原始OTDR曲线数据。
图2是小波分解1级高频系数。
图3是小波分解2级高频系数。
图4是小波分解3级高频系数。
图5是小波分解4级高频系数。
图6是小波分解5级高频系数。
图7是小波分解6级高频系数。
图8是小波分解7级高频系数。
图9是小波分解8级高频系数。
图10是小波分解多层级高频系数累加。
图11是OTDR事件检测结果。
具体实施方式
下面将结合本发明的附图,对本发明实施的技术方案进行详细描述。
本发明提出的OTDR曲线数据事件检测的方法。用小波变换模块对信号进行处理,小波分析可对信号进行多尺度分析,可以聚焦到信号的任意细节,通过不同层级的小波分解,信号可分解为不同层级的低频和高频部分,事件位置在高频中有体现,设置阈值将幅度大的高频系数保留出来,这些高频系数即为事件所在位置,再进行事件类型识别。本发明揭示的方法,能有效检测OTDR曲线数据中的事件,其具体方法包括:
S1,OTDR数据用小波变换进行不同层级小波分解,分解为低频和高频系数;
S2,对分解所得的高频系数,设置阈值,分离出幅度大的高频系数;
S3,根据分离出的高频系数对应位置,进行事件类型识别。
S1,小波变换使用Haar小波基,
Haar尺度函数定义为: ,
Haar小波基函数定义为: 。
S1,小波变换,或使用Daubechies小波基
Daubechies尺度函数定义为:,
Daubechies小波基函数定义为:,
,其中:和为滤波器系数,只有2M个非零项。
S1,用不同层级小波分解,小波分解层级根据一个脉宽时间内采集到的离散点数量决定,小波分解层级计算公式:,表示脉宽,表示采样率,L表示小波分解层级。参数L表示对原始OTDR进行小波分解的层级数,分解后,可得共L级高频系数部分。
OTDR事件在不同层级低频和高频系数中表现为不同的离散点数量,不同层级的离散点对应原始曲线数据的步长也不同。基于上述小波分解层级计算公式:,分解所得的共L级高频系数部分,事件在D1中表现的离散点数量是D2的2倍, D1中一个离散点步长长度为D2的1/2。
S2,分解所得高频系数,设置阈值,分离出幅度大的高频系数。OTDR事件在不同层级的高频系数中表现为:事件奇异性显著,会在不同层级的高频系数中都有体现。事件奇异性较小,低层级的高频系数表现较弱,在高层级的高频系数适当加强。事件奇异性非常微弱,任一层级的高频系数都表现较弱,需将所有层级的高频系数进行累加,才会有体现。
S3,根据分离出的幅度大的高频系数对应位置,结合原始OTDR曲线数据、分解后的各层级低频数据进行事件类型识别,完成OTDR曲线数据的事件检测。
如图1所示为原始OTDR曲线数据,通过本发明方法进行事件检测的结果如图11所示,图中序号1事件为“起始”,序号5事件为“结束”,序号2、3、4分别为检测到的“衰减”、“反射”、“衰减”事件。
对原始数据进行小波分解时,事件在高频系数中会有体现,如图2-图9,共8幅图,分别表示了小波分解后1级-8级的高频系数,序号2、3事件随着小波分解级数的上升,在高频系数中表现越明显,而对于序号4事件,1级-8级的高频系数都表现较弱,将所有1级-8级的高频系数进行累加,得到图10,序号4号事件才有体现,从而实现了事件的检测。
以上所述实施例仅阐述了本发明的几种实施方式,但并不能因此而解读为对本发明专利范围的限制。应当指出的是,对于本领域的技术人员来说,在不脱离本发明构思的前提下,还可以做出若干替换和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围以所附权利要求为准。
Claims (8)
1.一种OTDR曲线数据检测的方法,其特征是包括:
S1,OTDR数据用小波变换分解为不同层级低频和高频系数;
S2,对高频系数部分,设置阈值,分离出幅度大的高频系数;
S3,根据分离出的高频系数对应位置,进行事件类型识别。
2.根据权利要求1所述的方法,其特征是,所述S1中用小波变换,采用Haar小波基,
Haar尺度函数定义为: ,
Haar小波基函数定义为: 。
3.根据权利要求1所述的方法,其特征是,所述S1中用小波变换,采用Daubechies小波基,Daubechies尺度函数定义为:,
Daubechies小波基函数定义为:,
,其中:和为滤波器系数,只有2M个非零项。
4.根据权利要求1至3任一权利要求所述方法,其特征是,所述S1中用不同层级小波分解,小波分解层级根据一个脉宽时间内采集到的离散点数量决定,小波分解层级计算公式:,w表示脉宽,s表示采样率,L表示小波分解层级。
5.根据权利要求1至4任一权利要求所述方法,其特征是,所述S2中分离出幅度大的高频系数,小波分解后的任一层级,存在大于该阈值的高频系数时,将这些高频系数及所在层级保留出来。
6.根据权利要求1至4任一权利要求所述方法,其特征是,所述S2中分离出幅度大的高频系数,小波分解后的任一层级,不存在大于该阈值的高频系数,将所有层级的高频系数进行叠加,找出大于该阈值的叠加值。
7.根据权利要求1、2、3、4、5任一权利要求所述方法,其特征是,所述S3中分离出的幅度大的高频系数对应位置,幅度大的高频系数来自小波分解的不同层级,再进行事件识别。
8.根据权利要求1、2、3、4、6,任一权利要求所述方法,其特征是,所述S3中分离出的幅度大的高频系数对应位置,幅度大的高频系数来自小波分解后不同层级的累加,再进行事件识别。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610611889.6A CN107664571A (zh) | 2016-07-31 | 2016-07-31 | 一种光时域反射仪曲线数据事件检测的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610611889.6A CN107664571A (zh) | 2016-07-31 | 2016-07-31 | 一种光时域反射仪曲线数据事件检测的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107664571A true CN107664571A (zh) | 2018-02-06 |
Family
ID=61114540
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610611889.6A Pending CN107664571A (zh) | 2016-07-31 | 2016-07-31 | 一种光时域反射仪曲线数据事件检测的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107664571A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109861746A (zh) * | 2018-12-17 | 2019-06-07 | 中博信息技术研究院有限公司 | 一种基于小波变换动态降噪的otdr曲线数据分析方法 |
CN112051030A (zh) * | 2019-06-07 | 2020-12-08 | 爱斯福公司 | 重复otdr测量检测 |
US11879802B2 (en) | 2020-10-22 | 2024-01-23 | Exfo Inc. | Testing optical fiber link continuity using OTDR backscattering patterns |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2116164A1 (en) * | 1994-02-22 | 1995-08-23 | Anton Strifler | An optical fiber testing and monitoring system with a bi-directional communication link |
CN103196465A (zh) * | 2013-04-11 | 2013-07-10 | 电子科技大学 | 一种相敏光时域反射仪传感信号噪声分离及信号提取方法 |
CN103226028A (zh) * | 2013-05-02 | 2013-07-31 | 电子科技大学 | 一种相敏光时域反射仪扰动信号检测与识别方法 |
CN104052541A (zh) * | 2014-05-06 | 2014-09-17 | 上海温光自动化技术有限公司 | 光纤检测的显示方法、系统及光时域反射仪 |
-
2016
- 2016-07-31 CN CN201610611889.6A patent/CN107664571A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2116164A1 (en) * | 1994-02-22 | 1995-08-23 | Anton Strifler | An optical fiber testing and monitoring system with a bi-directional communication link |
CN103196465A (zh) * | 2013-04-11 | 2013-07-10 | 电子科技大学 | 一种相敏光时域反射仪传感信号噪声分离及信号提取方法 |
CN103226028A (zh) * | 2013-05-02 | 2013-07-31 | 电子科技大学 | 一种相敏光时域反射仪扰动信号检测与识别方法 |
CN104052541A (zh) * | 2014-05-06 | 2014-09-17 | 上海温光自动化技术有限公司 | 光纤检测的显示方法、系统及光时域反射仪 |
Non-Patent Citations (3)
Title |
---|
张伟: "小波算法在光纤测量中的应用", 《电子科技》 * |
李建中等: "基于ф-OTDR和POTDR结合的分布式光纤微扰传感系统", 《光子学报》 * |
陈永其等: "基于小波包分析的瞬时信号侦测研究", 《通信对抗》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109861746A (zh) * | 2018-12-17 | 2019-06-07 | 中博信息技术研究院有限公司 | 一种基于小波变换动态降噪的otdr曲线数据分析方法 |
CN112051030A (zh) * | 2019-06-07 | 2020-12-08 | 爱斯福公司 | 重复otdr测量检测 |
CN112051030B (zh) * | 2019-06-07 | 2023-11-28 | 爱斯福公司 | 重复otdr测量检测 |
US11879802B2 (en) | 2020-10-22 | 2024-01-23 | Exfo Inc. | Testing optical fiber link continuity using OTDR backscattering patterns |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106225907B (zh) | 一种基于φ-otdr技术的光纤振动识别系统及方法 | |
Aktas et al. | Deep learning based multi-threat classification for phase-OTDR fiber optic distributed acoustic sensing applications | |
CA2734717C (en) | Fibre optic acoustic sensing | |
CN105067104B (zh) | 一种复合型光纤传感系统和传感方法 | |
CN103115693B (zh) | 一种分布式光纤拉曼测温系统 | |
CN104217513B (zh) | 提高相位敏感光时域反射计识别入侵事件准确率的方法 | |
CN104240455A (zh) | 一种分布式光纤管线安全预警系统中的扰动事件识别方法 | |
CN108132092A (zh) | 基于自适应均值去阈值光纤振动事件识别方法 | |
CN107664571A (zh) | 一种光时域反射仪曲线数据事件检测的方法 | |
CN104729667A (zh) | 一种分布式光纤振动传感系统中的扰动类型识别方法 | |
CN103196465B (zh) | 一种相敏光时域反射仪传感信号噪声分离及信号提取方法 | |
CN103226028B (zh) | 一种相敏光时域反射仪扰动信号检测与识别方法 | |
JP2020052030A (ja) | マンホール位置特定方法及びマンホール位置特定システム | |
CN103344989B (zh) | 可控震源地震记录中脉冲噪声干扰的分析方法 | |
CN109323751A (zh) | 一种脉冲编码的分布式光纤振动传感方法及装置 | |
CN103308943A (zh) | 一种海洋地震资料处理中层间多次波衰减的方法及装置 | |
CN110335430A (zh) | 基于深度学习的管道安全监测系统、方法和装置 | |
Sha et al. | PIG tracking utilizing fiber optic distributed vibration sensor and YOLO | |
CN104301037A (zh) | 一种实现无源光网络光纤故障自动告警的方法及系统 | |
US20230152543A1 (en) | Impulse signal detection for buried cable protection using distributed fiber optic sensing | |
Wu et al. | Performance enhancement of ROTDR using deep convolutional neural networks | |
CN110147716A (zh) | 一种基于频域与空域相结合的sar图像内波区域检测方法 | |
CN108317402B (zh) | 一种分布式光纤传感定位预警分析方法 | |
CN115265750A (zh) | 一种光纤分布式声波传感系统及方法 | |
Zhao et al. | Coupled noise reduction in distributed acoustic sensing seismic data based on convolutional neural network |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180206 |