CN108088548A - 分布式光纤振动传感器高精度定位方法 - Google Patents
分布式光纤振动传感器高精度定位方法 Download PDFInfo
- Publication number
- CN108088548A CN108088548A CN201711190942.0A CN201711190942A CN108088548A CN 108088548 A CN108088548 A CN 108088548A CN 201711190942 A CN201711190942 A CN 201711190942A CN 108088548 A CN108088548 A CN 108088548A
- Authority
- CN
- China
- Prior art keywords
- curve
- optical fiber
- fiber vibration
- vibration sensor
- rayleigh scattering
- 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
- 239000013307 optical fiber Substances 0.000 title claims abstract description 26
- 238000000034 method Methods 0.000 title claims abstract description 18
- 238000000253 optical time-domain reflectometry Methods 0.000 claims abstract description 8
- 239000000835 fiber Substances 0.000 claims description 7
- 230000003321 amplification Effects 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 238000001914 filtration Methods 0.000 claims description 2
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 2
- 230000005693 optoelectronics Effects 0.000 claims 1
- 238000000926 separation method Methods 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H9/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
- G01H9/004—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means using fibre optic sensors
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
分布式光纤振动传感器高精度定位方法，属于光纤传感技术领域。本发明基于φ‑OTDR光纤振动传感原理，设计了一种基于瑞利散射曲线的时域方差的定位算法，该算法克服了分离平均法容易在分离处产生干扰尖峰，以及移动平均法运算量较大的问题。实验结果表明该算法可以提高定位精度，更容易找到事件点。
Description
技术领域
本发明属于光纤传感技术领域，尤其涉及一种分布式光纤振动传感系统中的事件点定位方法。
背景技术
φ—OTDR分布式光纤振动传感系统，采用光纤作为传感单元，终端为全光结构，通过测量发射脉冲到接收到干涉信号的时间差，即可实现对振动点的定位。适合大范围的安全和健康监测，比如大型工事、桥梁隧道、环境与公路工程、大型展览馆与会场、输油气管道等。随着我国油气管线输送项目建设的加快，国内油气管线对分布式光纤振动传感器的测量长度大幅提高，测量精度和定位准确等参数也要相应提高。
φ—OTDR分布式光纤振动传感系统是一种常见的光纤振动传感系统，向传感光纤中输入脉冲光，脉冲信号在光纤中传输会发生瑞利散射，产生瑞利后向散射信号，当外界有振动作用在光纤上某处时，由于弹光效应，该处的光纤参数会发生变化，导致瑞利散射曲线对应位置处的光强发生变化，通过信号处理，即可实现对振动事件的探测和定位，扰动位置可由下式计算：
式中：c为真空中光速△t为发出脉冲到收到瑞利曲线发生变化的时间，n为光纤纤芯折射率。
φ-OTDR型光纤振动传感系统定位时常用平均法滤除随机噪声的影响，一般采用的平均法有分离平均算法和移动平均算法。分离平均算法将连续的曲线分割开了，可能在两条分离平均法处理后的曲线间产生额外的尖峰。移动平均法虽然解决了这个问题，但是计算量比较大，难以实现较高实时性的定位功能。
发明内容
本发明所要解决的问题是：针对分离平均法和移动平均法存在的不足，设计一种新的算法，实现精确的定位。
本发明的解决方案是：设计了一种基于瑞利散射曲线时域方差的定位算法，步骤如下：
（1）向φ-OTDR系统中按一定的重复频率输入脉冲信号，脉冲信号在传感光纤中发生瑞利散射，产生后向瑞利散射信号，瑞利后向散射光具有强相干性会发生自干涉，送往光电探测器进行光电转化，经过放大、滤波后得到干涉光强I，，其中A为后向散射光波振幅，λ为光波波长，n为光纤纤芯折射率，W为脉宽；
(2) 采集多条瑞利后向散射曲线，将这些曲线按时域平均，得到一条平均曲线;
(3) 计算各条瑞利后向散射曲线与平均曲线的差值得到差值曲线;
(4) 在差值曲线中每次取若干条计算它们的时域方差，得到方差曲线。将得到的方差曲线取平均值，方差平均值曲线中的尖峰点即对应着振动所在位置。
本发明和现有技术相比的优点在于：本文提出了一种基于时域方差的定位算法，该算法克服了分离平均法容易在分离处产生干扰尖峰，以及移动平均法运算量较大的问题。
附图说明
图1为本发明的流程框图；
图2为φ-OTDR光纤振动传感系统示意图；
图3为采集到的瑞利后向散射曲线；
图4为经过差值运算后得到的差值曲线；
图5为对差值曲线计算时域方差后得到的方差曲线。
具体实施方式
如图2示，为φ-OTDR光纤振动传感系统示意图，激光器选用窄线宽激光器，输出功率为10mw,波长为1550nm的窄线宽激光、选用声光调制器为古奇（Gooch＆Housego）光纤耦合声光调制器、放大器选用上海拜安EDFA、探测器选用Thorlabs公司PDB410C光电探测器、测试用的光纤为单模光纤G652D，长度为6.5Km。采用的脉冲光的重复频率为1 kHz，脉冲宽度为300 ns，以20 MHz的采样率共采集了200条瑞利散射曲线。其中某条瑞利散射曲线如图3示，具体算法如下：
（1）采集了M（M=200）条瑞利散射曲线，将M条曲线写作t={t1, t2, t3, …, tM}。每次依次从M条曲线中取N条用于计算方差，设置一个步进因子n，表示每次取曲线的起始点相距n条。则最终可以获得K条方差曲线，其中K=[(M-N]/n]+1；
（2）首先计算M条瑞利散射曲线的时域平均值ta：
（3）然后将M条瑞利散射曲线减去平均值，得到M条差值曲线△t：
（4）设每一条△t均有Q个点，将差值曲线△t中的Q个点记作△t ={pi,1, pi,2, pi,3,…, pi,Q}，。接着，每次步进n条曲线，依次从M条曲线中取N条计算时域方差，得到K条方差曲线：
式中，，表示第i个点处的N条曲线的平均值；
（5）最终，将K条方差曲线叠加平均，得到最终的方差平均值曲线：
只要得到方差平均值曲线T中的峰值，即可得到对应的作用在光纤上的振动事件位置。
本发明的流程框图如图1所示。
本发明在4.5 km处敲击光纤模拟振动干扰，方差算法中取N = 6，n =3。得到如图4所示的差值曲线及图5所示方差曲线。由图4可以看出虽然在瑞利散射曲线图的4.5 km附近有尖峰，但是曲线的其他地方也有尖峰，容易对最终的定位结果造成干扰，由图5可以看出采用方差曲线可以对4.5公里处的振动点精确定位。
上面结合附图对本发明进行了示例性描述，显然本发明具体实现并不受上述方式的限制，只要采用了本发明的方法构思和技术方案进行的各种非实质性的改进，或未经改进将本发明的构思和技术方案直接应用于其它场合的，均在本发明的保护范围之内。
Claims (4)
1.分布式光纤振动传感器高精度定位方法，其特征在于：采用φ-OTDR光纤振动传感系统，向传感光纤中注入脉冲光，在光电探测器处就会接收到一组后向瑞利散射曲线，该散射曲线与光纤各点位置是对应的，后向瑞利散射光具有极强的相干性，因此光电探测器在环形器处探测的光实际上是脉冲宽度内后向瑞利散射光干涉的结果。
2.根据权利要求1所述的分布式光纤振动传感器高精度定位方法，其特征在于：后向瑞利散射经过光电探测器光电转换、放大及滤波后，得到输出光强信号。
3.根据权利要求1或2所述的分布式光纤振动传感器高精度定位方法，其特征在于：通过采集多条瑞利散射曲线，将这些曲线按时域平均，得到一条平均曲线，计算各条瑞利散射曲线与平均曲线的差值得到差值曲线，在差值曲线中每次取若干条计算它们的时域方差，得到方差曲线，将得到的方差曲线取平均值，方差平均值曲线中的尖峰点即对应着振动所在位置，根据尖峰推算出△t 。
4.根据权利要求1-3所述的分布式光纤振动传感器高精度定位方法，其特征在于：利用公式得到振动点位置，△t为发出脉冲到收到瑞利曲线发生变化的时间，c为真空中光速，n 为光纤纤芯折射率。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711190942.0A CN108088548A (zh) | 2017-11-24 | 2017-11-24 | 分布式光纤振动传感器高精度定位方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711190942.0A CN108088548A (zh) | 2017-11-24 | 2017-11-24 | 分布式光纤振动传感器高精度定位方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108088548A true CN108088548A (zh) | 2018-05-29 |
Family
ID=62172944
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711190942.0A Pending CN108088548A (zh) | 2017-11-24 | 2017-11-24 | 分布式光纤振动传感器高精度定位方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108088548A (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2020041972A (ja) * | 2018-09-13 | 2020-03-19 | 沖電気工業株式会社 | 振動検知光ファイバセンサ及び振動検知方法 |
| CN110954143A (zh) * | 2018-09-27 | 2020-04-03 | 冲电气工业株式会社 | 光纤传感器装置和光纤传感器系统 |
| JP2020153704A (ja) * | 2019-03-18 | 2020-09-24 | 沖電気工業株式会社 | 振動検知光ファイバセンサ及び振動検知方法 |
| CN118013401A (zh) * | 2024-04-10 | 2024-05-10 | 宁波联河光子技术有限公司 | 一种基于das的带式输送机振动误报抑制方法 |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5498867A (en) * | 1993-11-22 | 1996-03-12 | Sachio Uehara | Wavelength-division multiplex digital optical position sensor |
| EP2208016A1 (fr) * | 2007-11-09 | 2010-07-21 | Thales | Capteur optique pour mesurer la deformation au cours du temps d'une structure plane deformable |
| CN102280001A (zh) * | 2011-07-29 | 2011-12-14 | 电子科技大学 | 基于φ-otdr的分布式光纤围栏入侵检测与定位方法 |
| CN103196465A (zh) * | 2013-04-11 | 2013-07-10 | 电子科技大学 | 一种相敏光时域反射仪传感信号噪声分离及信号提取方法 |
| CN104132693A (zh) * | 2014-08-06 | 2014-11-05 | 电子科技大学 | 相位otdr系统中振动信号位置和频率的同时提取方法 |
| CN104198030A (zh) * | 2014-08-29 | 2014-12-10 | 中国石油天然气股份有限公司 | 基于相干瑞利散射的多路振动检测方法及其检测系统 |
| CN104964699A (zh) * | 2015-05-22 | 2015-10-07 | 北京交通大学 | 基于φ-OTDR光纤分布式扰动传感器的扰动判断方法和装置 |
| CN105783763A (zh) * | 2016-05-17 | 2016-07-20 | 安徽师范大学 | 一种高精度动态光纤应变传感装置及其传感方法 |
| CN106052842A (zh) * | 2016-08-05 | 2016-10-26 | 上海交通大学 | 可消衰落噪声的分布式光纤振动传感系统及其解调方法 |
-
2017
- 2017-11-24 CN CN201711190942.0A patent/CN108088548A/zh active Pending
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5498867A (en) * | 1993-11-22 | 1996-03-12 | Sachio Uehara | Wavelength-division multiplex digital optical position sensor |
| EP2208016A1 (fr) * | 2007-11-09 | 2010-07-21 | Thales | Capteur optique pour mesurer la deformation au cours du temps d'une structure plane deformable |
| CN102280001A (zh) * | 2011-07-29 | 2011-12-14 | 电子科技大学 | 基于φ-otdr的分布式光纤围栏入侵检测与定位方法 |
| CN103196465A (zh) * | 2013-04-11 | 2013-07-10 | 电子科技大学 | 一种相敏光时域反射仪传感信号噪声分离及信号提取方法 |
| CN104132693A (zh) * | 2014-08-06 | 2014-11-05 | 电子科技大学 | 相位otdr系统中振动信号位置和频率的同时提取方法 |
| CN104198030A (zh) * | 2014-08-29 | 2014-12-10 | 中国石油天然气股份有限公司 | 基于相干瑞利散射的多路振动检测方法及其检测系统 |
| CN104964699A (zh) * | 2015-05-22 | 2015-10-07 | 北京交通大学 | 基于φ-OTDR光纤分布式扰动传感器的扰动判断方法和装置 |
| CN105783763A (zh) * | 2016-05-17 | 2016-07-20 | 安徽师范大学 | 一种高精度动态光纤应变传感装置及其传感方法 |
| CN106052842A (zh) * | 2016-08-05 | 2016-10-26 | 上海交通大学 | 可消衰落噪声的分布式光纤振动传感系统及其解调方法 |
Non-Patent Citations (1)
| Title |
|---|
| 王鹏: "《Φ-OTDR的分布式光纤扰动传感系统阈值设定算法研究》", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2020041972A (ja) * | 2018-09-13 | 2020-03-19 | 沖電気工業株式会社 | 振動検知光ファイバセンサ及び振動検知方法 |
| WO2020054143A1 (ja) * | 2018-09-13 | 2020-03-19 | 沖電気工業株式会社 | 振動検知光ファイバセンサ及び振動検知方法 |
| JP7077887B2 (ja) | 2018-09-13 | 2022-05-31 | 沖電気工業株式会社 | 振動検知光ファイバセンサ及び振動検知方法 |
| CN110954143A (zh) * | 2018-09-27 | 2020-04-03 | 冲电气工业株式会社 | 光纤传感器装置和光纤传感器系统 |
| CN110954143B (zh) * | 2018-09-27 | 2022-03-08 | 冲电气工业株式会社 | 光纤传感器装置和光纤传感器系统 |
| JP2020153704A (ja) * | 2019-03-18 | 2020-09-24 | 沖電気工業株式会社 | 振動検知光ファイバセンサ及び振動検知方法 |
| JP7238507B2 (ja) | 2019-03-18 | 2023-03-14 | 沖電気工業株式会社 | 振動検知光ファイバセンサ及び振動検知方法 |
| CN118013401A (zh) * | 2024-04-10 | 2024-05-10 | 宁波联河光子技术有限公司 | 一种基于das的带式输送机振动误报抑制方法 |
| CN118013401B (zh) * | 2024-04-10 | 2024-07-26 | 宁波联河光子技术有限公司 | 一种基于das的带式输送机振动误报抑制方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| WO2018076551A1 (zh) | 能够实现全相位解调的分布式光纤传感系统及其测量方法 | |
| CN107101658B (zh) | 相位敏感光时域反射分布式光纤传感系统快速定位方法 | |
| CN101255951B (zh) | 提高油气管道泄漏和入侵检测的分布式光纤传感器性能的方法 | |
| CN101488805B (zh) | 光纤扰动探测方法及装置 | |
| CN109596205B (zh) | 一种基于延时光纤的双脉冲光纤振动传感方法 | |
| CN203561437U (zh) | 随机位置点光纤分布式声波传感装置 | |
| CN110031082B (zh) | 分布式光纤振动监测系统事件定位方法、装置及应用 | |
| CN103954349B (zh) | 一种分布式光纤振动传感系统的侧向定位方法 | |
| CN107664541A (zh) | 一种分布式光纤振动和温度融合传感系统及方法 | |
| CN107917738A (zh) | 一种同时测量温度、应变和振动的分布式光纤传感系统 | |
| CN105509868A (zh) | 相位敏感光时域反射光纤分布式传感系统相位计算方法 | |
| CN108088548A (zh) | 分布式光纤振动传感器高精度定位方法 | |
| CN106248119A (zh) | 一种分布式超高速扰动定量检测方法及装置 | |
| CN105004355B (zh) | 一种基于偏振态正交光脉冲对的φ‑otdr增强方法及装置 | |
| CN107505041A (zh) | 一种基于相位敏感光时域反射计的相位解调装置和方法 | |
| CN105157874B (zh) | 一种分布式边界安防监测系统及方法 | |
| CN201876324U (zh) | 分布式光纤拉曼温度传感器的双光源光路结构 | |
| CN104568119A (zh) | 一种单光源脉冲编码的光纤振动传感系统及其传感方法 | |
| CN101441092A (zh) | 基于相干光时域反射的周界防护传感定位系统 | |
| CN110518969B (zh) | 一种光缆振动的定位装置及方法 | |
| CN103900623A (zh) | 基于双声光调制器的光时域反射仪及其共模抑制方法 | |
| CN111811637B (zh) | 一种多径信息融合的车辆振动识别装置及判决方法 | |
| CN103017887A (zh) | 一种光纤振动传感系统及其检测方法 | |
| CN110501062A (zh) | 一种分布式光纤声音传感及定位系统 | |
| CN107664515B (zh) | 一种光纤周界传感定位系统及入侵扰动点位置的判定方法 |
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: 20180529 |