CN105277272A - 分布式光纤振动传感多点扰动定位算法 - Google Patents
分布式光纤振动传感多点扰动定位算法 Download PDFInfo
- Publication number
- CN105277272A CN105277272A CN201510698781.0A CN201510698781A CN105277272A CN 105277272 A CN105277272 A CN 105277272A CN 201510698781 A CN201510698781 A CN 201510698781A CN 105277272 A CN105277272 A CN 105277272A
- Authority
- CN
- China
- Prior art keywords
- optical fiber
- disturbance
- point
- cepstrum
- distributed optical
- 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
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Optical Transform (AREA)
Abstract
本发明属于分布式测量技术领域,具体为一种分布式光纤传感多点扰动定位算法。本发明方法运用倒谱分析的方法,首先对分布式光纤传感主机采集到的两路干涉信号进行相位还原解调得到两路干涉信号的相位差信号,对相位差信号进行傅里叶变换后取自然对数,最后进行傅里叶反变换得到倒谱,根据倒谱获取扰动信号的时延信息计算出多点扰动的位置信息,实现了分布式光纤传感定位系统中多点扰动同时出现的多点位置确定问题。
Description
技术领域
本发明属于分布式测量技术领域,具体涉及一种分布式光纤振动传感系统中的多点扰动定位算法。
背景技术
分布式光纤振动传感器具有获得被测量信息在空间和时间上连续分布的特点,灵敏度高,便于构成连续分布的传感网络,适合于长距离光信号的传输,在一些基础设施如通信电缆、输油管道,高压电网等以及一些具有特殊要求的安防系统中具有广阔的应用范围。
传统的分布式光纤振动传感扰动源定位方法主要有陷波点法和自相关时延法。陷波点法是通过扰动源引起的光相位差的频谱图中的频率缺失点来获得振动点的位置信息。但是该方法中的频率缺失点位置易受到其他因素的干扰使得定位精度下降,更无法处理多点同时扰动的情况。自相关时延法是对解调后的相位信号做自相关,利用不同路径光信号之间受到扰动的时延差对应于自相关函数副峰与主峰之间的位置差来实现扰动定位。但当多点扰动同时作用时,自相关函数的副峰难以分辨以致无法定位。
发明内容
本发明的目的是提出一种针对分布式光纤传感系统的多点扰动定位算法,以解决传统扰动源定位方法无法实现多点定位的问题。
本发明提出的分布式光纤传感多点扰动定位算法,具体步骤如下:
(1)首先,对分布式光纤传感主机采集到的两路干涉信号dataX和dataY进行相位还原解调得到两路干涉信号的相位差信号
(2)然后,对相位差信号进行傅里叶变换得到频谱P(ω);对频谱P(ω)取自然对数ln[P(ω)];
(3)对频谱的自然对数ln[P(ω)]进行傅里叶反变换得到倒谱C(t)=F-1{ln[P(ω)]};
(4)根据倒谱C(t),取其四个最大峰值对应的时间值从小到大分别为τ1,τ2,τ3,τ4,计算τ1+2τ2和τ2+2τ1,与τ4进行比较:
若τ1+2τ2与τ4值大小较近(一般地,两者差值小于10微秒),则将τ1和τ4根据求得两处扰动位置到线路末端的距离LA和LB;
若τ2+2τ1与τ4值大小较近(一般地,两者差值小于10微秒),则将τ2和τ4根据求得两处扰动位置到线路末端的距离LA和LB;
其中C是真空中光速,neff是光纤的有效折射率,τ是不同峰值对应的时间值,即为τ1,τ2,τ3,τ4。
与现有技术相比较,本发明的有益效果是:本算法解决了分布式光纤传感定位系统中多点扰动同时出现的多点位置确定问题,且定位算法简便快捷,定位精度高。
附图说明
图1示意了整个发明的算法流程图。
图2示意了实施例所采用的光路结构示意图。
图3示意了本发明处理得到的相位差信号示意图。
图4示意了本发明处理得到的倒谱信号示意图。
图中标号:1为光源;2为3*3耦合器;3为延迟线圈;4为2*2耦合器;5为扰动点A;6为扰动点B;7为PGC模块;8为FRM。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
实施例中分布式光纤传感器采用基于Mach-Zehnder和Sagnac混合干涉仪的结构,其结构框图如图2所示。根据光路结构,三路等能量光束中的两路光束经过光路后形成4路不同光束,分别为:
光束A:光源1→3*3耦合器2→延迟线圈3→2*2耦合器4→扰动点A5→扰动点B6→PGC模块7→FRM8→PGC模块7→扰动点B6→扰动点A5→2*2耦合器4→延迟线圈3→3*3耦合器2→光源1。
光束B:光源1→3*3耦合器2→2*2耦合器4→扰动点A5→扰动点B6→PGC模块7→FRM8→PGC模块7→扰动点B6→扰动点A5→2*2耦合器4→3*3耦合器2→光源1。
光束C:光源1→3*3耦合器2→延迟线圈3→2*2耦合器4→扰动点A5→扰动点B6→PGC模块7→FRM8→PGC模块7→扰动点B6→扰动点A5→2*2耦合器4→3*3耦合器2→光源1。
光束D:光源1→3*3耦合器2→2*2耦合器4→扰动点A5→扰动点B6→PGC模块7→FRM8→PGC模块7→扰动点B6→扰动点A5→2*2耦合器4→延迟线圈3→3*3耦合器2→光源1。
由于光纤延迟线长度远远大于激光器的相干长度,所以仅仅有光束C和D能够产生干涉,其余光束形成直流信号。这两路光束具有相同光程,在3*3耦合器2处发生干涉,分布式光纤传感主机采集到的两路干涉信号dataX和dataY。
图1示意了整个发明的算法流程图。首先对分布式光纤传感主机采集到的两路干涉信号dataX和dataY进行相位还原解调得到两路干涉信号的相位差信号当A,B两点同时施加扰动信号时,去程光路和回程光路均在A,B两点各受到扰动,一共四次扰动,四次扰动两两之间均存在一定的时间延迟,分别为τ56,2τ68,τ56+2τ68,2τ56+2τ68,τ56,τ56+2τ68,其中τ56,τ56+2τ68均重复一次,最后解调出的信号为 如图3中所示,其中s(t)为扰动信号,α是衰减因子,在计算中假设α均为1,τ为2个子信号之间的时延,n(t)为噪声。
对相位差信号进行傅里叶变换得到频谱P(ω),对频谱P(ω)取自然对数ln[P(ω)]。对频谱的自然对数ln[P(ω)]进行傅里叶反变换得到倒谱C(t)=F-1{ln[P(ω)]}。如图4中所示,根据倒谱C(t),其四个最大峰值对应的时间值从小到大分别为τ1,τ2,τ3,τ4,分别对应于时间延迟τ56,2τ68,τ56+2τ68,2τ56+2τ68,计算τ1+2τ2和τ2+2τ1,与τ4进行比较,若τ1+2τ2与τ4值大小较为相近,则将τ1和τ4根据求得两处扰动位置到线路末端的距离LA和LB,若τ2+2τ1与τ4值大小较为相近,则将τ2和τ4根据求得两处扰动位置到线路末端的距离LA和LB。
Claims (1)
1.一种分布式光纤振动传感多点扰动定位算法,其特征在于,具体步骤如下:
(1)首先,对分布式光纤传感主机采集到的两路干涉信号dataX和dataY进行相位还原解调得到两路干涉信号的相位差信号;
(2)然后,对相位差信号进行傅里叶变换得到频谱;对频谱取自然对数;
(3)对频谱的自然对数进行傅里叶反变换得到倒谱;
(4)根据倒谱,取其四个最大峰值对应的时间值从小到大分别为τ1,τ2,τ3,τ4,计算τ1+2τ2和τ2+2τ1,与τ4进行比较:
若τ1+2τ2与τ4值大小差值小于10微秒,则将τ1和τ4根据求得两处扰动位置到线路末端的距离LA和LB;
若τ2+2τ1与τ4值大小差值小于10微秒,则将τ2和τ4根据求得两处扰动位置到线路末端的距离LA和LB;
其中c是真空中光速,neff是光纤的有效折射率,τ是不同峰值对应的时间值,即为τ1,τ2,τ3,τ4。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510698781.0A CN105277272A (zh) | 2015-10-25 | 2015-10-25 | 分布式光纤振动传感多点扰动定位算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510698781.0A CN105277272A (zh) | 2015-10-25 | 2015-10-25 | 分布式光纤振动传感多点扰动定位算法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105277272A true CN105277272A (zh) | 2016-01-27 |
Family
ID=55146684
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510698781.0A Pending CN105277272A (zh) | 2015-10-25 | 2015-10-25 | 分布式光纤振动传感多点扰动定位算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105277272A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107014409A (zh) * | 2017-03-26 | 2017-08-04 | 天津大学 | 一种长距离光频域反射光纤分布式多点扰动传感方法 |
CN108303626A (zh) * | 2018-03-09 | 2018-07-20 | 广东电网有限责任公司珠海供电局 | 一种基于分布式光纤传感阵列的局部放电超声测量系统与方法 |
WO2019019735A1 (zh) * | 2017-07-27 | 2019-01-31 | 天津求实飞博科技有限公司 | 基于光频域反射长距离光纤分布式扰动传感信号处理方法 |
CN109668620A (zh) * | 2019-02-27 | 2019-04-23 | 深圳大学 | 一种基于干涉型分布式光纤振动传感器的定位方法 |
CN110057387A (zh) * | 2019-05-10 | 2019-07-26 | 南昌航空大学 | 直线型Sagnac分布式光纤传感系统的定位方法 |
CN111157102A (zh) * | 2020-01-02 | 2020-05-15 | 河海大学常州校区 | 一种分布式光纤传感系统中消除频率干扰的定位方法 |
CN111238552A (zh) * | 2020-02-27 | 2020-06-05 | 南昌航空大学 | 基于深度学习的分布式光纤传感系统扰动定位方法 |
CN111780857A (zh) * | 2020-06-05 | 2020-10-16 | 南京曦光信息科技有限公司 | 一种基于谐波累加的p-otdr系统的多点扰动定位检测方法 |
CN111780856A (zh) * | 2020-06-01 | 2020-10-16 | 哈尔滨工业大学 | 基于瑞利散射谱的相位谱分析的光纤分布式振动测量方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6834146B2 (en) * | 2002-01-03 | 2004-12-21 | The Boeing Company | Differential phase shift keyed demodulator system |
CN102538845A (zh) * | 2011-12-08 | 2012-07-04 | 北京航空航天大学 | 多点扰动定位方法、光纤分布式扰动传感器 |
CN102564476A (zh) * | 2011-10-28 | 2012-07-11 | 北京航空航天大学 | 多点扰动定位方法 |
CN102997945A (zh) * | 2011-09-16 | 2013-03-27 | 北京航空航天大学 | 光纤分布式扰动传感器的多点扰动定位方法 |
CN102997946A (zh) * | 2011-09-16 | 2013-03-27 | 北京航空航天大学 | 光纤分布式扰动传感器及其扰动定位方法 |
CN103090961A (zh) * | 2013-02-01 | 2013-05-08 | 华中科技大学 | 一种分布式光纤传感系统的扰动源定位方法 |
CN103743461A (zh) * | 2013-12-24 | 2014-04-23 | 南京邮电大学 | 分布式光纤振动传感器振源位置定位方法 |
CN103983342A (zh) * | 2014-05-26 | 2014-08-13 | 北京航天易联科技发展有限公司 | 一种确定传感光纤的震动位置的方法及光纤传感系统 |
CN104077267A (zh) * | 2014-06-27 | 2014-10-01 | 复旦大学 | 一种双重傅立叶变换扰动定位方法 |
CN104215271A (zh) * | 2014-07-30 | 2014-12-17 | 复旦大学 | 一种分布式光纤扰动监测系统中扰动位置的定位方法 |
-
2015
- 2015-10-25 CN CN201510698781.0A patent/CN105277272A/zh active Pending
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6834146B2 (en) * | 2002-01-03 | 2004-12-21 | The Boeing Company | Differential phase shift keyed demodulator system |
CN102997945A (zh) * | 2011-09-16 | 2013-03-27 | 北京航空航天大学 | 光纤分布式扰动传感器的多点扰动定位方法 |
CN102997946A (zh) * | 2011-09-16 | 2013-03-27 | 北京航空航天大学 | 光纤分布式扰动传感器及其扰动定位方法 |
CN102564476A (zh) * | 2011-10-28 | 2012-07-11 | 北京航空航天大学 | 多点扰动定位方法 |
CN102538845A (zh) * | 2011-12-08 | 2012-07-04 | 北京航空航天大学 | 多点扰动定位方法、光纤分布式扰动传感器 |
CN103090961A (zh) * | 2013-02-01 | 2013-05-08 | 华中科技大学 | 一种分布式光纤传感系统的扰动源定位方法 |
CN103743461A (zh) * | 2013-12-24 | 2014-04-23 | 南京邮电大学 | 分布式光纤振动传感器振源位置定位方法 |
CN103983342A (zh) * | 2014-05-26 | 2014-08-13 | 北京航天易联科技发展有限公司 | 一种确定传感光纤的震动位置的方法及光纤传感系统 |
CN104077267A (zh) * | 2014-06-27 | 2014-10-01 | 复旦大学 | 一种双重傅立叶变换扰动定位方法 |
CN104215271A (zh) * | 2014-07-30 | 2014-12-17 | 复旦大学 | 一种分布式光纤扰动监测系统中扰动位置的定位方法 |
Non-Patent Citations (3)
Title |
---|
杨德森等: "基于倒谱分析的水声信号被动定位时延估计算法研究", 《系统仿真学报》 * |
王志强等: "倒谱和复倒谱在被动声纳目标识别中的应用", 《声学技术》 * |
陈京惠: "基于Mach-Zehnder干涉仪分布式扰动传感器的定位方法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107014409B (zh) * | 2017-03-26 | 2019-05-28 | 天津大学 | 一种长距离光频域反射光纤分布式多点扰动传感方法 |
CN107014409A (zh) * | 2017-03-26 | 2017-08-04 | 天津大学 | 一种长距离光频域反射光纤分布式多点扰动传感方法 |
WO2019019735A1 (zh) * | 2017-07-27 | 2019-01-31 | 天津求实飞博科技有限公司 | 基于光频域反射长距离光纤分布式扰动传感信号处理方法 |
CN108303626A (zh) * | 2018-03-09 | 2018-07-20 | 广东电网有限责任公司珠海供电局 | 一种基于分布式光纤传感阵列的局部放电超声测量系统与方法 |
CN108303626B (zh) * | 2018-03-09 | 2024-04-12 | 广东电网有限责任公司珠海供电局 | 基于分布式光纤传感阵列的局部放电超声测量系统与方法 |
CN109668620A (zh) * | 2019-02-27 | 2019-04-23 | 深圳大学 | 一种基于干涉型分布式光纤振动传感器的定位方法 |
CN110057387B (zh) * | 2019-05-10 | 2021-01-19 | 南昌航空大学 | 直线型Sagnac分布式光纤传感系统的定位方法 |
CN110057387A (zh) * | 2019-05-10 | 2019-07-26 | 南昌航空大学 | 直线型Sagnac分布式光纤传感系统的定位方法 |
CN111157102A (zh) * | 2020-01-02 | 2020-05-15 | 河海大学常州校区 | 一种分布式光纤传感系统中消除频率干扰的定位方法 |
CN111238552B (zh) * | 2020-02-27 | 2021-06-22 | 南昌航空大学 | 基于深度学习的分布式光纤传感系统扰动定位方法 |
CN111238552A (zh) * | 2020-02-27 | 2020-06-05 | 南昌航空大学 | 基于深度学习的分布式光纤传感系统扰动定位方法 |
CN111780856A (zh) * | 2020-06-01 | 2020-10-16 | 哈尔滨工业大学 | 基于瑞利散射谱的相位谱分析的光纤分布式振动测量方法 |
CN111780856B (zh) * | 2020-06-01 | 2022-03-29 | 哈尔滨工业大学 | 基于瑞利散射谱的相位谱分析的光纤分布式振动测量方法 |
CN111780857A (zh) * | 2020-06-05 | 2020-10-16 | 南京曦光信息科技有限公司 | 一种基于谐波累加的p-otdr系统的多点扰动定位检测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105277272A (zh) | 分布式光纤振动传感多点扰动定位算法 | |
Zhu et al. | Characterizing thunder‐induced ground motions using fiber‐optic distributed acoustic sensing array | |
CN102879081B (zh) | 一种分布式光纤振动系统中的数据处理方法 | |
CN102636196A (zh) | 一种基于瑞利散射光谱相关系数的分布式扰动传感装置和解调方法 | |
CN104237871B (zh) | 一种基于相位补偿的时延差估计方法 | |
CN103047540B (zh) | 基于光纤传感的天然气管道泄漏监测光路系统 | |
CN103499383B (zh) | 一种提高光纤振动传感器定位精度的自相关定位方法 | |
CN107727227B (zh) | 基于φ-otdr的高压输电线路覆冰舞动监测方法 | |
CN111398901B (zh) | 一种扰动源多维空间定位系统和方法 | |
CN103090961A (zh) | 一种分布式光纤传感系统的扰动源定位方法 | |
CN105187121A (zh) | 一种通信光缆故障点地表位置定位方法及系统 | |
CN101813238B (zh) | Sagnac/Mach-Zehnder干涉仪分布光纤传感系统及其时域定位方法 | |
Pang et al. | Comparison between time shifting deviation and cross-correlation methods | |
CN101969344B (zh) | 基于光纤光弹效应的大区域声音监听系统 | |
CN103115633A (zh) | 利用相位生成载波降低干涉路径散(反)射光干扰的方法 | |
CN102419247B (zh) | 高精度反射式光纤湍流检测装置及方法 | |
CN110161458A (zh) | 一种基于光纤迈克尔逊干涉仪的声源定位系统 | |
CN110057387B (zh) | 直线型Sagnac分布式光纤传感系统的定位方法 | |
CN103148872B (zh) | 一种光纤分布式扰动传感器定位方法 | |
Gonzalez-Herraez et al. | Distributed acoustic sensing for seismic monitoring | |
CN111157102B (zh) | 一种分布式光纤传感系统中消除频率干扰的定位方法 | |
CN102589551A (zh) | 一种基于小波变换的船用光纤陀螺信号实时滤波方法 | |
Peng et al. | Method of precise sensing channel positioning forϕ-OTDR-based DAS | |
CN114136430B (zh) | 一种精准寻缆系统及方法 | |
Wu et al. | 2-Dimensional Timestamp-Based Sliding Window Cross-Correlation Localization Algorithm for Distributed Optical Fiber Vibration Sensing Interferometer |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into 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: 20160127 |