[go: up one dir, main page]

CN116592922A - Branch distributed positioning system based on single-core feedback interferometer - Google Patents

Branch distributed positioning system based on single-core feedback interferometer Download PDF

Info

Publication number
CN116592922A
CN116592922A CN202310435948.9A CN202310435948A CN116592922A CN 116592922 A CN116592922 A CN 116592922A CN 202310435948 A CN202310435948 A CN 202310435948A CN 116592922 A CN116592922 A CN 116592922A
Authority
CN
China
Prior art keywords
branch
disturbance
optical cable
path
main
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.)
Granted
Application number
CN202310435948.9A
Other languages
Chinese (zh)
Other versions
CN116592922B (en
Inventor
贾波
肖倩
李楚瑞
宋昱辰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fudan University
Original Assignee
Fudan University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Fudan University filed Critical Fudan University
Priority to CN202310435948.9A priority Critical patent/CN116592922B/en
Publication of CN116592922A publication Critical patent/CN116592922A/en
Application granted granted Critical
Publication of CN116592922B publication Critical patent/CN116592922B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01DMEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D5/00Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable
    • G01D5/26Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light
    • G01D5/32Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light
    • G01D5/34Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells
    • G01D5/353Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre
    • G01D5/35306Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre using an interferometer arrangement
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01DMEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D5/00Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable
    • G01D5/26Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light
    • G01D5/32Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light
    • G01D5/34Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells
    • G01D5/353Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre
    • G01D5/35306Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre using an interferometer arrangement
    • G01D5/35322Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre using an interferometer arrangement using interferometer with one loop with several directions of circulation of the light, e.g. Sagnac interferometer
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H9/00Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
    • G01H9/004Measuring 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)
  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Instruments For Measurement Of Length By Optical Means (AREA)

Abstract

The invention belongs to the technical field of optical fiber sensing, and particularly relates to a branch distributed positioning system based on a single-core feedback interferometer. The system combines two single-core feedback interferometers with different wavelengths through wavelength division multiplexing; the induction optical cable is divided into a main path and a branch path; multiple branches can extend out of the sensing optical cable; sensing that the optical cable is disturbed by the outside to cause the phase change of light in the optical fiber; demodulating interference signals with two different wavelengths generated by the same vibration to obtain two paths of phase difference signals; a group of signals with fixed time delay can be obtained through signal processing; determining whether the disturbance occurs on the main path or the branch path by using a positioning algorithm of the time delay estimation: the invention can not only distinguish whether the disturbance is applied to the main road or a certain branch road, but also accurately position the disturbance, and can reduce the complexity and the equipment cost of the sensing system.

Description

基于单芯反馈式干涉仪的支路分布式定位系统Branch Distributed Positioning System Based on Single-Core Feedback Interferometer

技术领域technical field

本发明属于光纤传感技术领域,具体涉及支路上的分布式光纤定位系统。The invention belongs to the technical field of optical fiber sensing, in particular to a distributed optical fiber positioning system on branch roads.

背景技术Background technique

外界的扰动会引起光纤内光学性质的变化,光纤传感技术正是基于此特性。如今光纤传感技术迅猛发展,其中分布式光纤传感技术由于其灵敏度高,动态范围大,监测距离远等优点,广泛应用在油气管道,边境安全,桥梁监测等领域。External disturbances will cause changes in the optical properties of the optical fiber, and optical fiber sensing technology is based on this characteristic. Nowadays, optical fiber sensing technology is developing rapidly. Among them, distributed optical fiber sensing technology is widely used in oil and gas pipelines, border security, bridge monitoring and other fields due to its high sensitivity, large dynamic range, and long monitoring distance.

目前,分布式光纤传感系统主要基于光时域反射(OTDR)和相干干涉法,其中干涉型光纤传感器获通过检测相干光的相位变化来探测外界信号,有着更高的实时性和灵敏性。在过去的几十年里,基于光纤干涉仪的各种结构已被提出,以实现分布式传感系统,如迈克尔逊干涉仪(MI),马赫-曾德尔干涉仪(MZI)和Sagnac干涉仪(SI)、单芯反馈式干涉仪(线性Sagnac干涉仪)。为适应各种复杂的实际工程环境并提升传感器性能,一些新型的干涉仪结构和定位算法被提出。At present, the distributed optical fiber sensing system is mainly based on optical time domain reflectometry (OTDR) and coherent interferometry. Among them, the interferometric optical fiber sensor can detect external signals by detecting the phase change of coherent light, which has higher real-time and sensitivity. In the past decades, various structures based on fiber optic interferometers have been proposed to realize distributed sensing systems, such as Michelson interferometer (MI), Mach-Zehnder interferometer (MZI) and Sagnac interferometer (SI), single-core feedback interferometer (linear Sagnac interferometer). In order to adapt to various complex actual engineering environments and improve sensor performance, some new interferometer structures and positioning algorithms have been proposed.

但在许多实际的应用场景中,如对重要的光通信网络或油气管道进行监测时,通信线路和管道线路的走向可能会存在星型、树型等多支路、多分岔的情况。在多支路的场景中,如果要对支路上的扰动进行探测和定位,使用现有的技术方案,就需要搭建并同时运行多个传感系统,会大大提高系统的复杂度与稳定性,并产生更高的成本。However, in many practical application scenarios, such as monitoring important optical communication networks or oil and gas pipelines, the directions of communication lines and pipelines may have multi-branches and multi-branches such as star and tree. In a multi-branch scenario, if you want to detect and locate the disturbance on the branch, using the existing technical solutions, you need to build and run multiple sensing systems at the same time, which will greatly increase the complexity and stability of the system. and incur higher costs.

本发明提出了一种基于单芯反馈式干涉仪的支路定位系统,实现只用一套传感系统对有支路和分叉的线路进行监控,能够探测主路与支路上的扰动并定位,能够分辨出扰动作用在主路还是支路上,并能够计算出扰动作用的具体位置。The present invention proposes a branch positioning system based on a single-core feedback interferometer, which realizes the monitoring of branch and bifurcated lines with only one set of sensing systems, and can detect and locate disturbances on the main road and the branch , it can be distinguished whether the disturbance is on the main road or the branch road, and the specific location of the disturbance can be calculated.

发明内容Contents of the invention

本发明针对有多条支路的监控线路,提出一种能延伸出多条支路的传感光纤线路结构,即基于单芯反馈式干涉仪的支路分布式定位系统。Aiming at a monitoring line with multiple branches, the present invention proposes a sensing optical fiber line structure capable of extending multiple branches, that is, a branch distributed positioning system based on a single-core feedback interferometer.

本发明提供的基于单芯反馈式干涉仪的支路分布式定位系统,是将两个不同波长的单芯反馈式干涉仪通过波分复用结合在一起,系统结构如图1所示。其中,感应光缆分为主路和支路;主路利用光缆中的一根芯,支路上利用同一根光缆中的两根芯;在靠近主路的一端,支路的两根芯分别和主路的光纤连接,在远离主路的一端,支路的两根光纤尾尾相互连接。通过这种连接方式,使感应光缆上可以延伸出多条支路,一般设置为3-6个支路。主路与支路连接方式如图2所示。The branch distributed positioning system based on the single-core feedback interferometer provided by the present invention combines two single-core feedback interferometers with different wavelengths through wavelength division multiplexing. The system structure is shown in Figure 1 . Among them, the induction optical cable is divided into main road and branch road; the main road uses one core in the optical cable, and the branch road uses two cores in the same optical cable; At the end far away from the main road, the two optical fiber tails of the branch road are connected to each other. Through this connection method, multiple branches can be extended from the sensing optical cable, generally 3-6 branches. The connection mode of the main road and the branch road is shown in Figure 2.

基于上结构,感应光缆受到外界的扰动,引起光纤内光的相位变化。将由相同振动产生的具有两个不同波长的干涉信号解调出来,得到两路相位差信号。通过进一步的信号处理,能得到一组具有固定时延的信号。利用时延估计的定位算法,确定扰动是发生在主路上还是支路上:Based on the above structure, the sensing optical cable is disturbed by the outside, causing the phase change of the light in the optical fiber. The interference signals with two different wavelengths generated by the same vibration are demodulated to obtain two phase difference signals. Through further signal processing, a set of signals with fixed time delay can be obtained. Use the positioning algorithm of time delay estimation to determine whether the disturbance occurs on the main road or the branch road:

(1)当扰动作用在主路上的,利用时延估计的定位算法直接解算出扰动发生的位置;(1) When the disturbance acts on the main road, use the positioning algorithm of time delay estimation to directly solve the location where the disturbance occurs;

(2)当扰动作用在支路上,先利用时延估计的定位算法确定扰动作用的支路,再使用频谱的陷波点定位法进行定位。在支路采用陷波点进行定位的处理过程中,由于直接利用某一路相位差信号陷波点进行定位,会存在多个干扰的陷波点,影响位置的确定,因此,本发明先将所得到的两路相位差信号相减,再对相减所得信号做频谱分析,通过陷波点确定支路上的扰动位置。(2) When the disturbance acts on the branch, first use the positioning algorithm of time delay estimation to determine the branch of the disturbance, and then use the notch point positioning method of the frequency spectrum to locate. In the process of using the notch point for positioning in the branch, since the notch point of a certain phase difference signal is directly used for positioning, there will be a plurality of disturbing notch points, which will affect the determination of the position. Therefore, the present invention first The obtained two phase difference signals are subtracted, and then the spectrum analysis is performed on the subtracted signal, and the disturbance position on the branch is determined through the notch point.

具体地,本发明提供的基于单芯反馈式干涉仪的支路分布式定位系统,参见图1,包括前端模块、感应光缆和尾端模块三大部分。前端模块由一个或若干光纤器件构成的分束合束单元1、分路器2和第一延迟光纤3组成;分束合束单元1包含一个3×3耦合器,其具有两个端口1b1、1b2,分路器2具有两个输入端口2b1、2b2,一个输出端口2a;从分束合束单元1输入的光被分成两束,分别经端口1b1、1b2输出;从光纤干涉路径返回的光分别进入端口1b1、1b2,在分束合束单元1内会合,发生干涉。端口1b1经第一延迟光纤3与分路器2的输入端口2b1相连,端口1b2直接与分路器2的输入端口2b2相连,分路器2的输出端口2a与感应光缆的一端相连。感应光缆分为主路光缆9和多条支路光缆10,其中主路9利用光缆中的一根芯,支路光缆10上利用同一根光缆中的两根芯。在靠近主路9的一端,支路10的两根芯分别和主路的光纤连接,在远离主路光缆9的一端,支路光缆10的两根光纤尾尾相互连接。主路上扰动点的位置11,支路上扰动点的位置12。尾端模块由两个波分复用器4和5,第二延时光纤6,相位调制装置7和反馈装置8组成;其中,波分复用器4具有一个输入端口4a、两个输出端口4b1、4b2,波分复用器5具有两个输入端口5b1、5b2和一个输出端口5a;波分复用器4的输入端口4a与主路光缆9的输出端连接,输出端口4b1经第二延时光纤6与波分复用器5的输入端口4b1连接,波分复用器4的输出端口4b2与波分复用器5的输入端口4b2连接;波分复用器5的输出端口5a依次与相位调制装置7和反馈装置8。Specifically, the branch distributed positioning system based on the single-core feedback interferometer provided by the present invention, as shown in FIG. 1 , includes three parts: a front-end module, a sensing optical cable and a tail-end module. The front-end module is composed of a beam splitting and combining unit 1, a splitter 2, and a first delay fiber 3 composed of one or several optical fiber components; the beam splitting and combining unit 1 includes a 3×3 coupler, which has two ports 1b1, 1b2, the splitter 2 has two input ports 2b1, 2b2, and one output port 2a; the light input from the beam splitting unit 1 is divided into two beams, which are respectively output through ports 1b1, 1b2; the light returned from the optical fiber interference path Enter ports 1b1 and 1b2 respectively, meet in the beam splitting and combining unit 1, and interfere. The port 1b1 is connected to the input port 2b1 of the splitter 2 via the first delay fiber 3, the port 1b2 is directly connected to the input port 2b2 of the splitter 2, and the output port 2a of the splitter 2 is connected to one end of the induction optical cable. The induction optical cable is divided into a main optical cable 9 and a plurality of branch optical cables 10, wherein the main optical cable 9 uses one core in the optical cable, and the branch optical cables 10 use two cores in the same optical cable. At one end close to the main path 9, the two cores of the branch path 10 are respectively connected to the optical fibers of the main path, and at the end far away from the main path optical cable 9, the two optical fiber tails of the branch path optical cable 10 are connected to each other. The location of the disturbance point on the main road is 11, and the location of the disturbance point on the branch is 12. Tail-end module is made up of two wavelength division multiplexers 4 and 5, the second delay fiber 6, phase modulation device 7 and feedback device 8; Wherein, wavelength division multiplexer 4 has an input port 4a, two output ports 4b1, 4b2, the wavelength division multiplexer 5 has two input ports 5b1, 5b2 and an output port 5a; the input port 4a of the wavelength division multiplexer 4 is connected to the output end of the main road optical cable 9, and the output port 4b1 passes through the second Delay fiber 6 is connected with the input port 4b1 of the wavelength division multiplexer 5, and the output port 4b2 of the wavelength division multiplexer 4 is connected with the input port 4b2 of the wavelength division multiplexer 5; the output port 5a of the wavelength division multiplexer 5 In turn with the phase modulation device 7 and the feedback device 8 .

设本发明系统中,波长分别为λ1、λ2的光,产生两对相干光束;于是:Assume that in the system of the present invention, the light with the wavelengths λ 1 and λ 2 generates two pairs of coherent light beams; then:

波长λ1的光的干涉路径为:The interference path of light with wavelength λ1 is:

(a)1b1→3→2b1→2a→9→10→4a→4b1→6→5b1→5a→7→8→7→5a→5b1→6→4b1→4a→10→9→2a→2b2→1b2;(a) 1b1→3→2b1→2a→9→10→4a→4b1→6→5b1→5a→7→8→7→5a→5b1→6→4b1→4a→10→9→2a→2b2→1b2 ;

(b)1b2→2b2→2a→9→10→4a→4b1→6→5b1→5a→7→8→7→5a→5b1→6→4b1→4a→10→9→2a→2b1→3→1b1;(b) 1b2→2b2→2a→9→10→4a→4b1→6→5b1→5a→7→8→7→5a→5b1→6→4b1→4a→10→9→2a→2b1→3→1b1 ;

同样的,波长为λ2的干涉路径为:Similarly, the interference path with wavelength λ2 is:

(c)1b1→3→2b1→2a→9→10→4a→4b2→5b2→5a→7→8→7→5a→5b2→4b2→4a→10→9→(c) 1b1→3→2b1→2a→9→10→4a→4b2→5b2→5a→7→8→7→5a→5b2→4b2→4a→10→9→

2a→2b2→1b2;2a→2b2→1b2;

(d)1b2→2b2→2a→9→10→4a→4b2→5b2→5a→7→8→7→5a→5b2→4b2→4a→10→9→2a→2b1→3→1b1;(d) 1b2→2b2→2a→9→10→4a→4b2→5b2→5a→7→8→7→5a→5b2→4b2→4a→10→9→2a→2b1→3→1b1;

其中,“→”表示路径,标号代表光所经过的器件以及该器件的端口。Wherein, "→" represents the path, and the label represents the device through which the light passes and the port of the device.

波长为λ1的光经路径a和路径b后,在分束合束单元11中形成稳定干涉,最终被光电探测器探测到。After the light with a wavelength of λ1 passes through the path a and the path b, it forms stable interference in the beam splitting and combining unit 11, and is finally detected by the photodetector.

波长为λ2的光经路径c和路径d后,在分束合束单元11中形成稳定干涉,最终被光电探测器探测到。After the light with a wavelength of λ2 passes through the paths c and d, it forms stable interference in the beam splitting and combining unit 11, and is finally detected by the photodetector.

由于这两个干涉光路中的光波长不同,故两者之间不会发生干涉。Since the light wavelengths in the two interfering optical paths are different, no interference will occur between the two.

对于干涉型传感系统,外界扰动会引起光的相位发生变化。而干涉仪的输出光强是相位变化的非线性函数,为了获得相位信号,需要对探测到的强度信号进行解调,从而得到导致相位变化的外界扰动信号。For interferometric sensing systems, external disturbances can cause changes in the phase of light. The output light intensity of the interferometer is a nonlinear function of the phase change. In order to obtain the phase signal, it is necessary to demodulate the detected intensity signal to obtain the external disturbance signal that causes the phase change.

在本系统中,假设在支路上某处施加一个扰动信号,波长为λ1的干涉信号为:In this system, assuming that a disturbance signal is applied somewhere on the branch, the interference signal with a wavelength of λ1 is:

波长为λ2的干涉信号为:The interference signal whose wavelength is λ2 is:

其中,A1,A2,B1,B2是与输入光功率大小相关的参数,和/>是扰动引起的相位差,φcos(ωt)为相位调制器加载的相位信号,/>和/>是由分束合束单元1中3×3耦合器引起的初始相位差;通过调节相位生成载波的幅度,获得IPD的最大值和最小值,为进一步的相位解调奠定基础。另一方面,相位载波信号可以将干涉信号调制到高频,实现干涉信号和低频噪声信号的分离,为消除低频信号提供可能。Among them, A 1 , A 2 , B 1 , B 2 are parameters related to the input optical power, and /> is the phase difference caused by the disturbance, φcos(ωt) is the phase signal loaded by the phase modulator, /> and /> It is the initial phase difference caused by the 3×3 coupler in the beam splitting and combining unit 1; by adjusting the amplitude of the phase-generated carrier, the maximum and minimum values of IPD are obtained, which lays the foundation for further phase demodulation. On the other hand, the phase carrier signal can modulate the interference signal to a high frequency, realize the separation of the interference signal and the low-frequency noise signal, and provide the possibility to eliminate the low-frequency signal.

进一步通过相位还原算法,如,反三角法的解调算法,能够分别解调出信号:Further through the phase restoration algorithm, such as the demodulation algorithm of the inverse triangulation method, the signals can be demodulated separately:

之后通过低通滤波器将相位载波信号φcos(ωt)滤除,可以得到外界扰动引起的相位差Afterwards, the phase carrier signal φcos(ωt) is filtered out by a low-pass filter, and the phase difference caused by external disturbance can be obtained

和/> and />

本发明中系统使用波分复用技术,两个不同波长的两个子干涉系统能监测到同一扰动经过不同路径传播的信号。当扰动信号施加在主路光缆上时,传感器监测获得的同一外界扰动信号引起的相位差和/>即为:In the present invention, the system uses wavelength division multiplexing technology, and two sub-interference systems with two different wavelengths can monitor the signals propagated by the same disturbance through different paths. When the disturbance signal is applied to the main optical cable, the sensor monitors the phase difference caused by the same external disturbance signal and /> That is:

上式中,c是真空光速,n为纤芯折射率,L1是第一延迟线圈3的长度,L2是第二延迟线圈6的长度,L是扰动点11到尾端模块8的光纤长度(不包括第二延迟线圈6的长度)。由光弹效应造成的系数α,使得 In the above formula, c is the speed of light in vacuum, n is the core refractive index, L1 is the length of the first delay coil 3, L2 is the length of the second delay coil 6, and L is the length of the optical fiber from the disturbance point 11 to the end module 8 ( The length of the second delay coil 6 is not included). The coefficient α caused by the photoelastic effect, such that

对同一扰动信号产生的相位差和/>进行变换,令:The phase difference generated for the same disturbance signal and /> To transform, let:

得到具有固定时延差的信号Δφ1m和Δφ2mSignals Δφ 1m and Δφ 2m with fixed delay difference are obtained:

固定的时间延迟Δτ=2nL/c。对(7)和(8)进行互相关运算:Fixed time delay Δτ=2nL/c. Perform cross-correlation operations on (7) and (8):

时间延迟可以通过R(t)的峰值位置来确定,传感距离L则可由下式得出:The time delay can be determined by the peak position of R(t), and the sensing distance L can be obtained by the following formula:

当扰动信号施加在支路光缆上时,传感器监测获得的同一外界扰动信号引起的相位差和/>可以写为:When the disturbance signal is applied to the branch optical cable, the sensor monitors the phase difference caused by the same external disturbance signal and /> can be written as:

Lx是扰动点12到支路远端距离,Lz是支路的长度,Lb是支路近端到尾端模块8的距离。L x is the distance from the disturbance point 12 to the far end of the branch, L z is the length of the branch, and L b is the distance from the near end of the branch to the end module 8 .

如果继续使用在主路上的定位算法,对同一扰动信号产生的相位差进行变换,令:If you continue to use the positioning algorithm on the main road, the phase difference for the same disturbance signal and To transform, let:

得到具有固定时延差的信号Δφ1b和Δφ2bSignals Δφ 1b and Δφ 2b with fixed delay difference are obtained:

固定的时间延迟在这个式子中并不包含Lx这一项。所以当在支路上的不同位置施加扰动时,时间延迟的方法都只能计算出Lb+Lz,即支路远端到尾端模块的距离。我们无法通过该方法来计算扰动点的具体位置。但在多支路的情况下,因为每一条支路的长度和在主路上的位置不同,在不同的支路上施加扰动时,通过此方法能仅能计算出受扰动支路远端到尾端模块的距离,从而分辨出是哪一条支路收到了扰动。fixed time delay The term L x is not included in this formula. Therefore, when disturbances are applied at different positions on the branch, the time delay method can only calculate L b +L z , that is, the distance from the far end of the branch to the tail end module. We cannot use this method to calculate the specific location of the disturbance point. However, in the case of multiple branches, because the length of each branch and the position on the main road are different, when disturbances are applied to different branches, this method can only calculate the distance from the far end to the tail end of the disturbed branch. Module distance, so as to distinguish which branch has received the disturbance.

在确定了受扰动的支路后,为了更进一步确定在支路上的具体扰动位置,采用频谱陷波点法进行定位。频谱陷波点法常见于Sagnac干涉仪、单芯反馈式干涉仪的扰动定位中。根据振动频谱分析原理,任何一个振动都可以分解为一系列单一振动频率的叠加。那么当频率为ωi的扰动施加在支路光缆上时引起的相位差和/>就可以写作:After the disturbed branch is determined, in order to further determine the specific disturbing position on the branch, the notch point method of frequency spectrum is used for positioning. The spectrum notch method is commonly used in the disturbance location of Sagnac interferometer and single-core feedback interferometer. According to the principle of vibration spectrum analysis, any vibration can be decomposed into a series of superposition of single vibration frequency. Then when the disturbance with frequency ω i is applied to the branch optical cable, the phase difference caused by and /> You can write:

其中,对式(16)进一步化简:in, To further simplify formula (16):

(17)(17)

如果直接对一个波长的相位差信号进行频谱分析,根据式(17),当:If the spectrum analysis is performed directly on the phase difference signal of one wavelength, according to formula (17), when:

时,在/>的频谱上于频率ωi对应的交流量为0,频率谱上会出现陷波点。此时,有三种可能的情况: at /> The AC quantity corresponding to frequency ω i on the frequency spectrum is 0, and a notch point will appear on the frequency spectrum. At this point, there are three possible situations:

cosωi(-tx)=0,cosωi(-tb-tz)=0或 cosω i (-t x )=0, cosω i (-t b -t z )=0 or

因为不同支路对应的Lb+Lz不同,且这两个长度由应用本系统的实际光路所决定,所以当cosωi(-tb-tz)=0时引起的频谱上的陷波点会对我们造成干扰,以至于难以找到由cosωi(-tx)=0引起的陷波点,阻碍支路上扰动点的定位。Because L b + L z corresponding to different branches is different, and these two lengths are determined by the actual optical path of the system, so when cosω i (-t b -t z ) = 0, the notch on the frequency spectrum caused The points will disturb us, so that it is difficult to find the notch point caused by cosω i (-t x )=0, which hinders the location of the disturbance point on the branch.

利用本发明提出的新型的双波长的陷波点定位法,先将支路光缆上扰动引起的两波长的相位差相减,并做进一步的化简,可以得到:Utilize the novel dual-wavelength notch point positioning method proposed by the present invention, first subtract the phase difference of the two wavelengths caused by the disturbance on the branch optical cable, and do further simplification to obtain:

由式(18)可知,当时,在/>的频谱上于频率ωi对应的交流量为0,在频谱上存在明显的陷波点。这种情况可分为以下两种可能:From formula (18), it can be seen that when at /> The frequency ω i corresponds to 0 AC on the frequency spectrum, and there are obvious notch points in the frequency spectrum. This situation can be divided into the following two possibilities:

1)当cosωi(-tx)=0时,其中k为自然数;此时陷波点频率记为fn(k),则:1) When cosω i (-t x )=0, Where k is a natural number; at this time, the frequency of the notch point is recorded as f n (k), then:

与陷波点频率所对应的扰动点的位置为:The location of the disturbance point corresponding to the frequency of the notch point is:

2)当或sinωi(t2)=0时,/>或/>其中k为自然数。由于L1和L2可以取很小,那么与其对应的第一个陷波点频率fn(1)就会很大。只要选择适当的第一延迟线圈和第二延迟线圈长度,就可以避免fn (k)对fn(k)的干扰。2) when Or when sinω i (t 2 )=0, /> or /> where k is a natural number. Since L 1 and L 2 can be taken very small, the corresponding first notch frequency f n (1) will be very large. As long as the lengths of the first delay coil and the second delay coil are selected appropriately, the interference of f n (k) on f n (k) can be avoided.

因此,对光电探测器得到的两个波长的四个干涉信号做相位解调,计算出两个波长干涉仪的相位差和/>再将两路相位差相减,之后对其做傅里叶变换得到频域谱,从频域谱中找到陷波点,即可根据式(20)求得Lx,从而能确定作用在支路上的扰动点的具体位置。Therefore, the phase demodulation is performed on the four interference signals of two wavelengths obtained by the photodetector, and the phase difference of the two wavelength interferometers is calculated. and /> Then subtract the two-way phase difference, and then perform Fourier transform on it to obtain the frequency domain spectrum, find the notch point from the frequency domain spectrum, and then obtain Lx according to formula (20), so that it can be determined to act on the branch road The specific location of the disturbance point.

本发明提出了一种可以延伸出多条支路的传感光缆布局方式。传感光缆分为主路和支路,主路上利用光缆中的一根芯,支路利用光缆中的两根芯。支路与主路连接处,支路中两芯分别与主路中一芯前后相连;支路尽头,两芯相互连接。通过此连接方式,光纤中传输的光束先从主路中进入支路,在支路光缆中走过一个来回,再次从支路中回到主路。The invention proposes a sensor optical cable layout mode that can extend multiple branches. The sensing optical cable is divided into a main road and a branch road. The main road uses one core in the optical cable, and the branch road uses two cores in the optical cable. At the junction of the branch road and the main road, the two cores in the branch road are respectively connected to the front and back of one core in the main road; at the end of the branch road, the two cores are connected to each other. Through this connection method, the light beam transmitted in the optical fiber first enters the branch circuit from the main circuit, goes back and forth in the branch cable, and returns to the main circuit from the branch circuit again.

本发明利用时间延迟估计的定位算法对传感光缆主路上的扰动进行准确的传感定位,对传感光缆支路上的扰动进行辨别,能够分辨具体是哪一条支路受到了扰动。The invention utilizes a positioning algorithm of time delay estimation to accurately sense and locate the disturbance on the main path of the sensing optical cable, and distinguish the disturbance on the branch paths of the sensing optical cable, and can distinguish which branch path is disturbed.

本发明通过将光缆上同一扰动产生两波长相位差信号相减,再对这个差值信号做频谱分析,找出频谱上的陷波点来确定支路上扰动的具体位置。The invention subtracts two wavelength phase difference signals generated by the same disturbance on the optical cable, and then performs spectrum analysis on the difference signal to find notch points on the spectrum to determine the specific position of the disturbance on the branch.

本发明为具有多条支路或岔路的监测线路提供了一套新的多支路光纤传感系统,使用一套系统便可对主路和支路上的扰动进行传感定位,降低了传感系统的复杂度和设备成本。The present invention provides a new set of multi-branch optical fiber sensing system for monitoring lines with multiple branches or branch roads. Using a set of systems can sense and locate disturbances on the main road and branch roads, reducing the sensor System complexity and equipment cost.

附图说明Description of drawings

图1是本发明的一种支路定位的系统结构。Fig. 1 is a branch location system structure of the present invention.

图2是感应光缆上主路与支路连接方式示意图。Fig. 2 is a schematic diagram of the connection mode of the main circuit and the branch circuit on the induction optical cable.

图3是本发明实施例的支路定位系统的具体结构。Fig. 3 is the specific structure of the branch positioning system of the embodiment of the present invention.

图4是在实施例中,在三条支路上的10km处施加人工扰动,通过本发明中时间延迟算法测算出的定位结果。Fig. 4 is the positioning result measured and calculated by the time delay algorithm in the present invention by applying artificial disturbance at 10 km on the three branch roads in the embodiment.

图5是在实施例中,在三条支路上的10km处施加人工扰动,通过本发明中陷波点算法测算出的定位结果。Fig. 5 is the positioning result calculated by the notch point algorithm in the present invention by applying artificial disturbance at 10 km on the three branch roads in the embodiment.

图中标号:1为分束合束单元,2为光纤分路器,3为延迟光纤;1b1、1b2为分束合束单元1的输出端口,2b1、2b2为分路器2的两个输入端口,2a为分路器2的输出端口2a;4和5为两个波分复用器,4a为波分复用器4的输入端口,4b1、4b2为波分复用器4的两个输出端口,5b1、5b2为波分复用器5的两个输入端口,波分复用器5的输出端口;6为延时光纤,7为相位调制装置,8为反馈装置,如带尾纤的反射镜;9为主路光缆,10为支路光缆。11为主路上扰动点,12为支路上扰动点。13为主路光缆,14为支路光缆。15为超辐射发光二极管(SLD),16为掺铒光纤放大器,17、18、20、21分别是四个光电探测器,19和22为两个波分复用器,23为对称3X3光纤耦合器,24为对称2X2光纤耦合器,25是第一延时光纤。26为主路光缆,27是主路上延伸出的三条支路光缆,28是主路上扰动施加的位置,29是支路上扰动施加的位置。30和31是两个波分复用器,32为第二延时光纤,33为相位调制器,34为法拉第旋转镜。Numbers in the figure: 1 is the beam splitting and combining unit, 2 is the fiber splitter, 3 is the delay fiber; 1b1 and 1b2 are the output ports of the beam splitting and combining unit 1, and 2b1 and 2b2 are the two inputs of the splitter 2 port, 2a is the output port 2a of splitter 2; 4 and 5 are two wavelength division multiplexers, 4a is the input port of wavelength division multiplexer 4, 4b1, 4b2 are two wavelength division multiplexers 4 Output ports, 5b1 and 5b2 are two input ports of the wavelength division multiplexer 5, and the output port of the wavelength division multiplexer 5; 6 is a delay fiber, 7 is a phase modulation device, and 8 is a feedback device, such as a pigtail reflector; 9 is the main road optical cable, and 10 is the branch optical cable. 11 is the disturbance point on the main road, and 12 is the disturbance point on the branch road. 13 is the main road optical cable, and 14 is the branch optical cable. 15 is a superluminescent light-emitting diode (SLD), 16 is an erbium-doped fiber amplifier, 17, 18, 20, and 21 are four photodetectors respectively, 19 and 22 are two wavelength division multiplexers, and 23 is a symmetrical 3X3 fiber coupling device, 24 is a symmetrical 2X2 fiber coupler, and 25 is a first delay fiber. 26 is the optical cable for the main road, 27 is the three branch optical cables extending from the main road, 28 is the position where the disturbance is applied on the main road, and 29 is the position where the disturbance is applied on the branch road. 30 and 31 are two wavelength division multiplexers, 32 is a second delay fiber, 33 is a phase modulator, and 34 is a Faraday rotating mirror.

具体实施方式Detailed ways

下面通过具体实施例结合附图进一步介绍本发明。The present invention will be further described below through specific embodiments in conjunction with the accompanying drawings.

图1为本发明支路定位系统结构。Fig. 1 is the structure of the branch positioning system of the present invention.

图2是感应光缆上主路与支路连接方式示意图。13为主路光缆,14为主路光缆上延伸出的一条支路光缆。其中主路利用光缆中的一芯,支路上利用同一根光缆中的两根芯。在靠近主路的一端,支路的两根芯分别和主路的光纤连接,在远离主路9的一端,支路10的两根光纤尾尾相互连接。粗箭头指示了感应光缆中光束的传播路径,光束先从主路中进入支路,在支路光缆中走过一个来回,再次从支路中回到主路。Fig. 2 is a schematic diagram of the connection mode of the main circuit and the branch circuit on the induction optical cable. 13 is the main optical cable, and 14 is a branch optical cable extending from the main optical cable. One core of the optical cable is used for the main path, and two cores of the same optical cable are used for the branch path. At one end close to the main road, the two cores of the branch road are respectively connected to the optical fibers of the main road, and at the end far away from the main road 9, the two optical fiber tails of the branch road 10 are connected to each other. The thick arrow indicates the propagation path of the light beam in the induction optical cable. The light beam first enters the branch from the main road, travels back and forth in the branch optical cable, and returns to the main road from the branch again.

图3所示。15为超辐射发光二极管(SLD),16为掺铒光纤放大器,17、18、20、21分别是四个光电探测器,19和22为两个波分复用器,23为对称3X3光纤耦合器,24为对称2X2光纤耦合器,25是第一延时光纤。感应光缆可以分为主路和支路,26为主路,27是主路上延伸出的三条支路,28是主路上扰动施加的位置,29是支路上扰动施加的位置。30和31是两个波分复用器,32为第二延时光纤,33为相位调制器,34为法拉第旋转镜。Figure 3 shows. 15 is a superluminescent light-emitting diode (SLD), 16 is an erbium-doped fiber amplifier, 17, 18, 20, and 21 are four photodetectors respectively, 19 and 22 are two wavelength division multiplexers, and 23 is a symmetrical 3X3 fiber coupling device, 24 is a symmetrical 2X2 fiber coupler, and 25 is a first delay fiber. The induction optical cable can be divided into a main road and a branch road, 26 is the main road, 27 is three branch roads extending from the main road, 28 is the position of the disturbance on the main road, and 29 is the position of the disturbance on the branch road. 30 and 31 are two wavelength division multiplexers, 32 is a second delay fiber, 33 is a phase modulator, and 34 is a Faraday rotating mirror.

超辐射发光二极管15的中心波长为1550nm,带宽40nm,其相干长度在几十微米左右,作为本系统的光源,系统中的四个波分复用器19、22、30、31能将超辐射发光二极管15发出的光分为两个波长λ1=1553nm,λ2=1548nm。输出光的光功率通过16调整后通过波分复用器19被分为两束,一束经过延时光纤25和2X2光纤耦合器24,另一束直接经过2X2光纤耦合器24。两束光继续通过传感光纤后进入波分复用器30,波分复用器30将光束分为两个不同波长,相比波长为λ2的光束,λ1会多经过一段延时光纤32,之后波分复用器31将两束光合并。光束被相位调制器33调制后通过法拉第旋转镜34反射沿原路径返回,最终回到前端中的3X3光纤耦合器23中形成稳定干涉,两个波长的光再被波分复用器19和波分复用器22分开,四个光电探测器17、18、20、21采集干涉信号。使用高速数据采集卡实现A/D转换,采样率为5Mbps。使用基于Labview的处理程序,实现信号的处理算法。本实施例中,选取第一延迟线圈25为200m,第二延迟线圈32的长度为4km。主路光缆总长4km,靠近尾端的位置延伸出三条15km长的支路。The central wavelength of the superradiant light-emitting diode 15 is 1550nm, the bandwidth is 40nm, and its coherent length is about tens of microns. The light emitted by the light emitting diode 15 is divided into two wavelengths λ 1 =1553nm and λ 2 =1548nm. The optical power of the output light is adjusted by 16 and divided into two beams by the wavelength division multiplexer 19 , one beam passes through the delay fiber 25 and the 2X2 fiber coupler 24 , and the other beam directly passes through the 2X2 fiber coupler 24 . The two beams of light continue to enter the wavelength division multiplexer 30 after passing through the sensing fiber, and the wavelength division multiplexer 30 divides the light beam into two different wavelengths. Compared with the light beam with a wavelength of λ2 , λ1 will pass through a longer delay fiber 32, and then the wavelength division multiplexer 31 combines the two beams of light. After the light beam is modulated by the phase modulator 33, it is reflected by the Faraday rotating mirror 34 and returns along the original path, and finally returns to the 3X3 fiber coupler 23 in the front end to form stable interference. The multiplexer 22 separates the four photodetectors 17, 18, 20, 21 to collect interference signals. A/D conversion is realized by using a high-speed data acquisition card with a sampling rate of 5Mbps. Use the processing program based on Labview to realize the signal processing algorithm. In this embodiment, the length of the first delay coil 25 is selected to be 200 m, and the length of the second delay coil 32 is 4 km. The total length of the main optical cable is 4km, and three 15km-long branch roads are extended near the end.

用小铁锤分别在三条支路上的约10km处(由OTDR测得距离为9.98km)施加一个扰动,连续敲击20次,使用时间延迟法能够计算出扰动所在支路对应的支路远离主路一端,距离系统远端的光纤长度。敲击多次,计算出的结果如图4所示。敲击支路27a的计算结果为75000米左右,支路27b的计算结果为45000米左右,支路27c的结果为15000米左右。说明,主路上同一位置延伸出的支路,能通过时延法分辨出受扰动的支路。通过时间延迟法判断出受扰动的支路后,进一步对两路相位还原信号相减后进行频谱分析,从一阶陷波点的位置能计算出扰动作用在支路上具体的位置。多次敲击三条支路上10km处扰动的定位结果如图5所示。大部分的定位结果都为9.958km,能够确定出扰动点在支路上的具体位置,与实际位置误差为0.022km。Use a small hammer to apply a disturbance at about 10km on the three branch roads (the distance measured by OTDR is 9.98km), and knock 20 times in a row. Using the time delay method, it can be calculated that the branch where the disturbance is located is far away from the main road. At one end of the road, the length of the optical fiber from the far end of the system. Tap multiple times, and the calculated results are shown in Figure 4. The calculation result of the knocking branch 27a is about 75000 meters, the calculation result of the branch 27b is about 45000 meters, and the result of the branch 27c is about 15000 meters. It shows that the branches extending from the same position on the main road can distinguish the disturbed branch through the time delay method. After the disturbed branch is judged by the time delay method, the spectrum analysis is carried out after subtracting the two phase recovery signals, and the specific position of the disturbance on the branch can be calculated from the position of the first-order notch point. The positioning results of the disturbance at 10 km on the three branch roads after multiple knocks are shown in Figure 5. Most of the positioning results are 9.958km, which can determine the specific position of the disturbance point on the branch road, and the error with the actual position is 0.022km.

Claims (4)

1. A branch distributed positioning system based on a single-core feedback interferometer is characterized in that two single-core feedback interferometers with different wavelengths are combined together through wavelength division multiplexing; the induction optical cable is divided into a main path and a branch path; the main path uses one core in the optical cable, and the branch path uses two cores in the same optical cable; two cores of the branch are respectively connected with the optical fibers of the main road at one end close to the main road; at one end far away from the main road, two optical fiber tails of the branch road are connected with each other; by the connection mode, a plurality of branches extend out of the sensing optical cable;
based on the sensing optical cable structure, the sensing optical cable is disturbed by the outside to cause the phase change of light in the optical fiber; demodulating interference signals with two different wavelengths generated by the same vibration to obtain two paths of phase difference signals; a group of signals with fixed time delay can be obtained through signal processing; determining whether the disturbance occurs on the main path or the branch path by using a positioning algorithm of the time delay estimation:
(1) When the disturbance acts on the main road, the position where the disturbance occurs is directly solved by using a positioning algorithm of time delay estimation;
(2) When the disturbance acts on the branch, determining the branch of the disturbance by using a positioning algorithm of time delay estimation, and then positioning by using a notch point positioning method of a frequency spectrum; in the processing process of positioning the branch by adopting the notch point, the two obtained phase difference signals are subtracted, then spectrum analysis is carried out on the subtracted signals, and the disturbance position on the branch is determined through the notch point.
2. The single-core feedback interferometer-based branch distributed positioning system of claim 1, comprising three major parts, a front end module, an inductive fiber optic cable, and a tail end module; wherein:
the front end module consists of a beam splitting and combining unit (1) formed by one or a plurality of optical fiber devices, a splitter (2) and a first delay optical fiber (3); the beam splitting and combining unit (1) comprises a 3x3 coupler with two ports (1 b1, 1b 2), and the splitter (2) is provided with two input ports (2 b1, 2b 2) and an output port (2 a); light input from the beam splitting and combining unit (1) is split into two beams, and the two beams are output through two ports (1 b1, 1b 2) respectively; light returning from the optical fiber interference path enters ports (1 b1, 1b 2) respectively, and is converged in the beam splitting and combining unit (1) to generate interference; the first port (1 b 1) of the beam splitting and combining unit (1) is connected with the first input port (2 b 1) of the splitter (2) through the first delay optical fiber (3), the second port (1 b 2) of the beam splitting and combining unit (1) is directly connected with the second input port (2 b 2) of the splitter (2), and the output port (2 a) of the splitter (2) is connected with one end of the sensing optical cable;
the induction optical cable is divided into a main optical cable (9) and a plurality of branch optical cables (10), wherein one core is utilized in the main optical cable (9), and two cores in the same optical cable are utilized in the branch optical cable (10); two cores of the branch optical cable (10) are respectively connected with optical fibers of the main optical cable (9) at one end close to the main optical cable (9), and two optical fiber tails of the branch optical cable (10) are connected with each other at one end far away from the main optical cable (9); a position (11) of a disturbance point on the main optical cable, and a position (12) of a disturbance point on the branch optical cable;
the tail end module consists of two wavelength division multiplexers (4 and 5), a second delay optical fiber (6), a phase modulation device (7) and a feedback device (8); wherein the first wavelength division multiplexer (4) has one input port (4 a), two output ports (4 b1, 4b 2), and the second wavelength division multiplexer (5) has two input ports (5 b1, 5b 2) and one output port (5 a); an input port (4 a) of the first wavelength division multiplexer (4) is connected with an output end of the main optical cable (9), a first output port (4 b 1) of the first wavelength division multiplexer (4) is connected with a first input port (4 b 1) of the second wavelength division multiplexer (5) through a second delay optical fiber (6), a second output port (4 b 2) of the first wavelength division multiplexer (4) is connected with a second input port (4 b 2) of the second wavelength division multiplexer (5), and an output port (5 a) of the second wavelength division multiplexer (5) is sequentially connected with the phase modulation device (7) and the feedback device (8).
3. The single-core feedback interferometer-based branch distributed positioning system of claim 2, wherein the wavelengths are λ, respectively 1 、λ 2 Generating two pairs of coherent light beams; then:
wavelength lambda 1 The interference path of light is:
(a)1b1→3→2b1→2a→9→10→4a→4b1→6→5b1→5a→7→8→7→5a→5b1→6→4b1→4a→10→9→2a→2b2→1b2;
(b)1b2→2b2→2a→9→10→4a→4b1→6→5b1→5a→7→8→7→5a→5b1→6→4b1→4a→10→9→2a→2b1→3→1b1;
similarly, the wavelength is lambda 2 The interference path of (a) is:
(c)1b1→3→2b1→2a→9→10→4a→4b2→5b2→5a→7→8→7→5a→5b2→4b2→4a→10→9→2a→2b2→1b2;
(d)1b2→2b2→2a→9→10→4a→4b2→5b2→5a→7→8→7→5a→5b2→4b2→4a→10→9→2a→2b1→3→1b1;
wherein "→" denotes a path, and a reference numeral denotes a device through which light passes and a port of the device;
wavelength lambda 1 After passing through the path a and the path b, stable interference is formed in the beam splitting and combining unit 1, and finally the light is detected by the photoelectric detector;
wavelength lambda 2 After passing through the path c and the path d, stable interference is formed in the beam splitting and combining unit 1, and finally, the light is detected by the photodetector.
4. The single-core feedback interferometer-based branch distributed positioning system of claim 3, wherein the specific process for identifying and determining the external disturbance signal is as follows:
assuming that a disturbing signal is applied somewhere on the branch, the wavelength is lambda 1 The interference signals of (a) are:
wavelength lambda 2 The interference signals of (a) are:
wherein ,A1 ,A 2 ,B 1 ,B 2 Is a parameter related to the magnitude of the input optical power, and />Is the phase difference caused by disturbance, phi cos (ωt) is the phase signal loaded by the phase modulator,/-> and />Is an initial phase difference caused by a 3x3 coupler in the beam splitting and combining unit (1); by adjusting the amplitude of the phase-generated carrier, I is obtained PD Maximum and minimum of (2); on the other hand, phase carrier signalsModulating the interference signal to high frequency to realize the separation of the interference signal and the low-frequency noise signal;
the signals are demodulated by a phase recovery algorithm, namely a demodulation algorithm of an inverse trigonometry, respectively:
then filtering the phase carrier signal phi cos (omega t) by a low-pass filter to obtain the phase difference caused by external disturbance and />
Two sub-interference systems with different wavelengths in the system monitor signals of the same disturbance transmitted through different paths;
when disturbance signals are applied to the main optical cable, the sensor monitors the obtained phase difference caused by the same external disturbance signals and />The method comprises the following steps:
in the above formula, c is the vacuum light velocity, n is the refractive index of the fiber core, L1 is the length of the first delay coil (3), L2 is the length of the second delay coil (6), and L is the disturbance point (11) on the main optical cable toThe length of the optical fiber of the tail end module (8); coefficient alpha caused by photoelastic effect such that
Phase difference generated for the same disturbance signal and />And (3) performing transformation, and enabling:
obtaining a signal delta phi with a fixed time delay difference 1m and Δφ2m
A fixed time delay Δτ=2nl/c; performing cross-correlation operation on the steps (7) and (8):
the time delay is determined by the peak position of R (t), and the sensing distance L is given by:
when disturbance signal is applied to the branchWhen the optical cable is on, the sensor monitors the phase difference caused by the same external disturbance signal obtained and />The writing is as follows:
L x is the distance from the disturbance point (12) on the branch optical cable to the far end of the branch, L z Is the length of the branch, L b Is the distance from the proximal end of the branch to the tail end module (8);
if the positioning algorithm on the main road is continuously used, the phase difference generated for the same disturbance signal and />And (3) performing transformation, and enabling:
obtaining a signal delta phi with a fixed time delay difference 1b and Δφ2b
Fixed time delay
Under the condition of multiple branches, as the length of each branch and the position of each branch on a main road are different, when disturbance is applied to different branches, the distance from the far end of the disturbed branch to the tail end module can be calculated, so that the branch which is disturbed can be distinguished;
after determining a disturbed branch, determining a specific disturbance position on the branch, and positioning by adopting a frequency spectrum notch point method; according to the vibration spectrum analysis principle, any vibration can be decomposed into a series of superposition of single vibration frequencies; then when the frequency is omega i Phase difference caused when disturbance of (b) is applied to branch optical cable and />Writing:
wherein ,further simplifying formula (16):
the method comprises the steps of subtracting phase differences of two wavelengths caused by disturbance on a branch optical cable by adopting a double-wavelength notch point positioning method, and further simplifying to obtain:
as can be seen from formula (18), whenAt->Is spectrally at frequency omega i The corresponding alternating current quantity is 0, and obvious notch points exist on the frequency spectrum; there are two possibilities for this:
1) When cos omega i (-t x ) When the value of the sum is =0,wherein k is a natural number; the notch frequency is denoted as f n (k) Then:
the positions of the disturbance points corresponding to the notch point frequencies are as follows:
2) When (when)Or sin omega i (t 2 ) When=0,>or->Wherein k is a natural number; due to L 1 and L2 Can be made very small, then the first notch point frequency f corresponding thereto n (1) Will be very large; by selecting the appropriate lengths of the first delay coil and the second delay coil, f can be avoided n (k) For f n (k) Is a part of the interference of the (c).
CN202310435948.9A 2023-04-23 2023-04-23 Branch distributed positioning system based on single-core feedback interferometer Active CN116592922B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310435948.9A CN116592922B (en) 2023-04-23 2023-04-23 Branch distributed positioning system based on single-core feedback interferometer

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310435948.9A CN116592922B (en) 2023-04-23 2023-04-23 Branch distributed positioning system based on single-core feedback interferometer

Publications (2)

Publication Number Publication Date
CN116592922A true CN116592922A (en) 2023-08-15
CN116592922B CN116592922B (en) 2024-10-01

Family

ID=87598267

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310435948.9A Active CN116592922B (en) 2023-04-23 2023-04-23 Branch distributed positioning system based on single-core feedback interferometer

Country Status (1)

Country Link
CN (1) CN116592922B (en)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5402231A (en) * 1992-08-24 1995-03-28 Mcdonnell Douglas Corporation Distributed sagnac sensor systems
CN101578506A (en) * 2007-01-10 2009-11-11 住友电气工业株式会社 Method for measuring temperature distribution of object, and sensor unit
CN102519492A (en) * 2011-12-26 2012-06-27 复旦大学 Distributed single core feedback interference optical path structure having low background light
CN103047540A (en) * 2011-09-14 2013-04-17 中国石油天然气集团公司 Natural gas pipe leakage monitoring optical path system based on optical fiber sensing
CN108168686A (en) * 2018-03-14 2018-06-15 中国科学院半导体研究所 Dual wavelength distribution type fiber-optic sound sensor-based system
US20180292240A1 (en) * 2015-12-25 2018-10-11 Tianjin University Distributed optical fiber disturbance positioning system based on the asymmetric dual mach-zehnder interference, and positioning method thereof
CN109141490A (en) * 2018-09-20 2019-01-04 天津理工大学 A kind of fibre-optical sensing device and demodulation method of disturbance waveform and position measurement simultaneously
CN113324568A (en) * 2021-05-21 2021-08-31 复旦大学 Distributed optical fiber sensing positioning system based on asymmetric fusion interferometer
CN113324569A (en) * 2021-05-21 2021-08-31 复旦大学 Dual-wavelength unidirectional ring type distributed optical fiber sensing and positioning system

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5402231A (en) * 1992-08-24 1995-03-28 Mcdonnell Douglas Corporation Distributed sagnac sensor systems
CN101578506A (en) * 2007-01-10 2009-11-11 住友电气工业株式会社 Method for measuring temperature distribution of object, and sensor unit
CN103047540A (en) * 2011-09-14 2013-04-17 中国石油天然气集团公司 Natural gas pipe leakage monitoring optical path system based on optical fiber sensing
CN102519492A (en) * 2011-12-26 2012-06-27 复旦大学 Distributed single core feedback interference optical path structure having low background light
US20180292240A1 (en) * 2015-12-25 2018-10-11 Tianjin University Distributed optical fiber disturbance positioning system based on the asymmetric dual mach-zehnder interference, and positioning method thereof
CN108168686A (en) * 2018-03-14 2018-06-15 中国科学院半导体研究所 Dual wavelength distribution type fiber-optic sound sensor-based system
CN109141490A (en) * 2018-09-20 2019-01-04 天津理工大学 A kind of fibre-optical sensing device and demodulation method of disturbance waveform and position measurement simultaneously
CN113324568A (en) * 2021-05-21 2021-08-31 复旦大学 Distributed optical fiber sensing positioning system based on asymmetric fusion interferometer
CN113324569A (en) * 2021-05-21 2021-08-31 复旦大学 Dual-wavelength unidirectional ring type distributed optical fiber sensing and positioning system

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
GAO HAO ET.: "A long-haul and high-accuracy fiber length measurement technology based on Mach–Zehnder and Sagnac hybrid interferometer", MEASUREMENT SCIENCE AND TECHNOLOGY, 26 July 2020 (2020-07-26), pages 1 - 8 *
YUCHEN SONG ET AL.: "Fiber Optic Sensor Based on Linear Sagnac Interferometer for Branch Localization", JOURNAL OF LIGHTWAVE TECHNOLOGY, vol. 41, no. 20, 15 October 2023 (2023-10-15), pages 6621 - 6630 *
卞庞: "利用波分复用技术消除瑞利散射噪声的双波长sagnac干涉仪的研究", 中国博士学位论文全文数据库(电子期刊)工程科技Ⅱ辑, no. 01, 15 January 2016 (2016-01-15), pages 030 - 4 *

Also Published As

Publication number Publication date
CN116592922B (en) 2024-10-01

Similar Documents

Publication Publication Date Title
Sun et al. Distributed fiber-optic vibration sensor using a ring Mach-Zehnder interferometer
CN110864714B (en) Distributed sensing system based on Michelson-Sagnac fiber optic interferometer
CN103196584B (en) Measurement method for temperature and stress in fiber and Brillouin optical time domain reflectometer
CN102506913B (en) Interference type optical fiber distribution disturbance sensor and disturbance location method thereof
CN102064884B (en) Long-distance distributed optical fiber positioning interference structure based on wavelength division multiplexing (WDM)
KR101000974B1 (en) Measurement method of color dispersion characteristics of optical waveguide sample using interference fringe measurement system
CN109141490B (en) Optical fiber sensing device and demodulation method for simultaneously measuring disturbance waveform and position
CN102538846A (en) Phase-sensitive optical time domain reflectometer type optical fiber distributed disturbing sensor of double-arm pulse optical interference
CN113324568B (en) Distributed optical fiber sensor positioning system based on asymmetric fusion interferometer
CN105973450B (en) Optical fiber Fizeau interference arrays distributed vibration sensing system and method
CN110940364B (en) A distributed two-arm sensing system based on Michelson interferometer
CN102564476B (en) Multipoint disturbance positioning method
CN118010145A (en) Distributed vibration sensing system based on multi-core optical fiber
CN116499573A (en) Single-core bidirectional optical fiber distributed vibration detection system based on forward light sensing
Sun et al. Distributed disturbance sensor based on a novel Mach-Zehnder interferometer with a fiber-loop
CN117968822B (en) A multi-point vibration positioning method based on forward distributed sensing and phase demodulation
CN107314823B (en) The Method for Phase Difference Measurement and device of interferometric optical fiber sensor
CN116592922B (en) Branch distributed positioning system based on single-core feedback interferometer
CN113804412A (en) Optical fiber device micro-chromaticity dispersion measuring device based on ring light path structure
CN111157102A (en) A Positioning Method for Eliminating Frequency Interference in Distributed Optical Fiber Sensing System
CN102496231B (en) Long-distance Trunk Safety Optical Fiber Wavelength Division Multiplexing Early Warning System
CN103487067B (en) Utilize wavelength-division multiplex technique to subdue the method and system of interference path scattering, reverberation interference
Sun et al. Distributed fiber-optic sensor with a ring Mach-Zehnder interferometer
CN111157100B (en) All-fiber sensing positioning system and positioning method based on feedback
CN114111855B (en) Distributed optical fiber sensing positioning system based on two-way Michelson interferometer

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
GR01 Patent grant
GR01 Patent grant