CN111238551B - 分布式相位敏感光时域反射仪传感系统及相位提取方法 - Google Patents

分布式相位敏感光时域反射仪传感系统及相位提取方法 Download PDF

Info

Publication number
CN111238551B
CN111238551B CN202010071038.3A CN202010071038A CN111238551B CN 111238551 B CN111238551 B CN 111238551B CN 202010071038 A CN202010071038 A CN 202010071038A CN 111238551 B CN111238551 B CN 111238551B
Authority
CN
China
Prior art keywords
signal
light
phase
disturbance
optical fiber
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.)
Active
Application number
CN202010071038.3A
Other languages
English (en)
Other versions
CN111238551A (zh
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.)
University of Science and Technology Beijing USTB
Original Assignee
University of Science and Technology Beijing USTB
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 University of Science and Technology Beijing USTB filed Critical University of Science and Technology Beijing USTB
Priority to CN202010071038.3A priority Critical patent/CN111238551B/zh
Publication of CN111238551A publication Critical patent/CN111238551A/zh
Application granted granted Critical
Publication of CN111238551B publication Critical patent/CN111238551B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01D—MEASURING 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/00—Mechanical 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/26—Mechanical 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/32—Mechanical 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/34—Mechanical 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/353—Mechanical 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/35338—Mechanical 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 other arrangements than interferometer arrangements
    • G01D5/35354—Sensor working in reflection
    • G01D5/35358—Sensor working in reflection using backscattering to detect the measured quantity
    • G01D5/35361—Sensor working in reflection using backscattering to detect the measured quantity using elastic backscattering to detect the measured quantity, e.g. using Rayleigh backscattering
    • 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

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Optical Transform (AREA)

Abstract

本发明提供一种分布式相位敏感光时域反射仪传感系统及相位提取方法，能够提高扰动位置定位和扰动相位提取的准确性。所述方法包括：在发送端，通过声光调制器将信号光调制成脉冲光，脉冲光经过环形器注入传感光纤，从传感光纤反射回的瑞利散射光通过环形器输出，其中，脉冲光在声光调制器频移的作用下成为单边带信号；在接收端，将本振光和环形器输出的瑞利散射光在第二耦合器中进行混频，其中,混频之后的光信号为单边带信号,混频之后的光信号通过单个光电探测器进行光电转换，对光电转换后的电信号进行带通滤波并取滤波后电信号中的一个边带，对提取的边带进行数字信号处理，得到光纤扰动信号的位置及相位。本发明涉及分布式光纤传感技术领域。

Description

分布式相位敏感光时域反射仪传感系统及相位提取方法
技术领域
本发明涉及分布式光纤传感技术领域，特别是指一种分布式相位敏感光时域反射仪传感系统及相位提取方法。
背景技术
随着科学技术及工程建设的蓬勃发展，电网中的电缆与通信网中的光缆遍布在各个角落，例如，跨海大桥。这些工程往往横跨几十甚至上百公里，单单依靠人力对这些大型工程的安防及扰动监测是远远不够的，相位敏感光时域反射仪(Φ-OTDR)技术可以代替人力实现长距离的安防及扰动监测。通过分析Φ-OTDR传感系统接收的瑞利散射光，可以实现探测距离内的扰动信号定位以及其相位提取。目前，Φ-OTDR传感系统中使用的调制器多种多样，如光纤声光调制器(AOM)、电光调制器(Electro-Optic Modulator，EOM)等将连续的激光调制成脉冲光。EOM采用同相/正交(IQ)调制，调制方式复杂，并且EOM的消光比较小，很容易造成Φ-OTDR传感系统扰动的误判，影响扰动信号的定位结果，更进一步影响扰动信号的相位提取及分析；因此，现有技术中，常用AOM与平衡光电探测器(BPD)结合，利用外差相干检测来实现扰动信号相位的提取，但是BPD在使信号功率增大一倍的同时，噪声功率同样会增加一倍，导致扰动位置定位和扰动相位提取的准确性下降。
发明内容
本发明要解决的技术问题是提供一种分布式相位敏感光时域反射仪传感系统及相位提取方法，以解决现有技术所存在的由于噪声功率增加，导致扰动位置定位和扰动相位提取的准确性下降的问题。
为解决上述技术问题，本发明实施例提供一种分布式相位敏感光时域反射仪传感系统，包括：窄线宽激光器、第一耦合器、声光调制器、环形器、第二耦合器、光电探测器和数字信号处理器；其中，
窄线宽激光器，用于发射激光并将发射的激光通过第一耦合器分成两路，其中，一路作为信号光被声光调制器调制成脉冲光并在声光调制器频移的作用下成为单边带信号，作为单边带信号的脉冲光经过环形器注入传感光纤，从传感光纤反射回的瑞利散射光通过环形器输出；另一路作为本振光发送至第二耦合器；
第二耦合器，用于对本振光和环形器输出的瑞利散射光进行混频，混频之后的光信号通过单个光电探测器进行光电转换；其中，混频之后的光信号为单边带信号；
数字信号处理器，用于对光电转换后的电信号进行带通滤波并取滤波后电信号中的一个边带，对提取的边带进行数字信号处理，得到光纤扰动信号的位置及相位。
进一步的，所述系统还包括：第一掺铒光纤放大器和第一带通滤波器；
所述第一掺铒光纤放大器，用于对作为单边带信号的脉冲光进行功率放大；
所述第一带通滤波器，用于对放大后的脉冲光进行滤波处理，滤波后的脉冲光经过环形器注入传感光纤。
进一步的，所述系统还包括：第二掺铒光纤放大器和第二带通滤波器；
所述第二掺铒光纤放大器，用于对环形器输出的瑞利散射光进行功率放大；
所述第二带通滤波器，用于对放大后的瑞利散射光进行滤波处理，滤波后的瑞利散射光输入至第二耦合器。
本发明实施例还提供一种应用于权利要求1所述传感系统的相位提取方法，包括：
在发送端，通过声光调制器将信号光调制成脉冲光，脉冲光经过环形器注入传感光纤，从传感光纤反射回的瑞利散射光通过环形器输出；其中，脉冲光在声光调制器频移的作用下成为单边带信号；
在接收端，将本振光和环形器输出的瑞利散射光在第二耦合器中进行混频，其中，混频之后的光信号为单边带信号，混频之后的光信号通过单个光电探测器进行光电转换，对光电转换后的电信号进行带通滤波并取滤波后电信号中的一个边带，对提取的边带进行数字信号处理，得到光纤扰动信号的位置及相位。
进一步的，所述脉冲光经过环形器注入传感光纤，从传感光纤反射回的瑞利散射光通过环形器输出包括：
将调制后的脉冲光通过放大、带通滤波之后经过环形器注入传感光纤，从传感光纤反射回的瑞利散射光sRBS(t)的时域表达式为：
Figure GDA0002761014220000031
其中，t为时间，As(t)为瑞利散射光的幅值，f0为窄线宽激光器的频率，fb为声光调制器的频移量，
Figure GDA0002761014220000032
为瑞利散射光的瞬时相位偏移。
进一步的，本振光为载波，本振光sLO(t)的时域表达式为：
Figure GDA0002761014220000033
其中，AL为本振光的幅值，
Figure GDA0002761014220000034
为本振光的瞬时相位偏移。
进一步的，本振光和环形器输出的经过放大、带通滤波之后的瑞利散射光在第二耦合器中混频后，输出的光信号smix(t)的时域表达式为：
Figure GDA0002761014220000035
进一步的，光电转换后得到的电信号I(t)的时域表达式为：
Figure GDA0002761014220000036
其中，AS 2(t)为信号自拍频干扰项。
进一步的，所述对光电转换后的电信号进行带通滤波并取滤波后电信号中的一个边带，对提取的边带进行数字信号处理，得到光纤扰动信号的位置及相位包括：
对光电转换后得到的电信号进行带通滤波并取滤波后电信号中的一个边带，得到瑞利散射光功率IPF(t)表达式：
Figure GDA0002761014220000037
其中，时间t＝t0+NT，N为提取的信号IPF(t)中有完整瑞利散射光周期的个数，T为瑞利散射光的周期，t0为提取的信号IPF(t)中第一个不完整周期瑞利散射光持续时间；
去除t0时间内的信号，得到t’＝NT时间内的信号，将t’时间内的信号转换为N×T的数据矩阵，其中，每行数据与传感光纤中随距离变化的瑞利散射光相对应；
以数据矩阵的行为单位，确定数据矩阵中行之间数据的功率差，通过功率差实现扰动位置定位；其中，无扰动的位置功率差为0，有扰动的位置功率差不为0，以此来实现扰动信号的定位；
以数据矩阵的行为单位，确定数据矩阵中行之间数据的相位差，将扰动位置处瑞利散射光相位差与扰动之前无扰动位置瑞利散射光相位差做差，得到光纤扰动信号相位。
进一步的，所述将扰动位置处瑞利散射光相位差与扰动之前无扰动位置瑞利散射光相位差做差，得到光纤扰动信号相位包括：
设无扰动位置对应的时间为ta+nT，扰动位置对应的时间为tb+nT，根据相位差公式：
Figure GDA0002761014220000041
分别确定扰动位置和无扰动位置的瑞利散射光相位差，并将扰动位置和无扰动位置的瑞利散射光相位差做差，得到光纤扰动信号相位
Figure GDA0002761014220000042
Figure GDA0002761014220000043
其中，t1为第一个瑞利散射光周期中扰动位置对应的时间，IPF *(t1)为瑞利散射光功率IPF(t)表达式的共轭，ta为第一个瑞利散射光周期内无扰动位置(参考位置)对应的时间，tb为第一个瑞利散射光周期内扰动位置对应的时间，n＝2…N,0<t1≤T，
Figure GDA0002761014220000044
为扰动信号引起的瑞利散射光的相位变化，
Figure GDA0002761014220000045
为本振光的瞬时相位变化。
本发明的上述技术方案的有益效果如下：
上述方案中，利用声光调制器的频移功能，使得经过声光调制器频移的单边带信号在光电探测器接收之后，与光电转换过程中的信号自拍频干扰(SSBI)不发生重叠，具有隔离于SSBI的天然保护带，以便数字信号处理器通过带通滤波滤出不含SSBI的接收信号，从而排除SSBI对接收信号的干扰；采用单光电探测器进行光电转换，能够降低光电转换过程中引入的噪声，在长探测距离、高扰动频率的场景中发挥出更加明显的优势；在数字信号处理部分只取电信号的一个边带即可实现光纤扰动信号的定位及相位的检测，这样，在提高扰动位置定位和扰动相位提取的准确性的同时，还能实现光纤扰动信号的定位及相位的快速检测。
附图说明
图1为本发明实施例提供的分布式相位敏感光时域反射仪传感系统的结构示意图；
图2为本发明实施例提供的应用于分布式相位敏感光时域反射仪传感系统的相位提取方法的流程示意图；
图3为本发明实施例提供的信号频谱示意图，其中，图3中的(a)、(b)、(c)、(d)分别对应图1中的①②③④；
图4为本发明实施例提供的数字信号处理流程示意图。
具体实施方式
为使本发明要解决的技术问题、技术方案和优点更加清楚，下面将结合附图及具体实施例进行详细描述。
本发明针对现有的由于噪声功率增加，导致扰动位置定位和扰动相位提取的准确性下降的问题，提供一种分布式相位敏感光时域反射仪传感系统及相位提取方法。
实施例一
如图1所示，本发明实施例还提供一种分布式相位敏感光时域反射仪(Φ-OTDR)传感系统，包括：窄线宽激光器、第一耦合器、声光调制器、环形器、第二耦合器、光电探测器和数字信号处理器；其中，
窄线宽激光器，用于发射激光并将发射的激光通过第一耦合器分成两路，其中，一路作为信号光被声光调制器调制成脉冲光并在声光调制器频移的作用下成为单边带信号，作为单边带信号的脉冲光经过环形器注入传感光纤，从传感光纤反射回的瑞利散射光通过环形器输出；另一路作为本振光发送至第二耦合器；
第二耦合器，用于对本振光和环形器输出的瑞利散射光进行混频，混频之后的光信号通过单个光电探测器进行光电转换；其中，混频之后的光信号为单边带信号；
数字信号处理器，用于对光电转换后的电信号进行带通滤波并取滤波后电信号中的一个边带，对提取的边带进行数字信号处理，得到光纤扰动信号的位置及相位。
本发明实施例所述的分布式相位敏感光时域反射仪传感系统，利用声光调制器的频移功能，使得经过声光调制器频移的单边带信号在光电探测器接收之后，与光电转换过程中的信号自拍频干扰(SSBI)不发生重叠，具有隔离于SSBI的天然保护带，以便数字信号处理器通过带通滤波滤出不含SSBI的接收信号，从而排除SSBI对接收信号的干扰；采用单光电探测器进行光电转换，能够降低光电转换过程中引入的噪声，在长探测距离、高扰动频率的场景中发挥出更加明显的优势；在数字信号处理部分只取电信号的一个边带即可实现光纤扰动信号的定位及相位的检测，这样，在提高扰动位置定位和扰动相位提取的准确性的同时，还能实现光纤扰动信号的定位及相位的快速检测。
在前述分布式相位敏感光时域反射仪传感系统的具体实施方式中，进一步地，所述系统还包括：第一掺铒光纤放大器和第一带通滤波器；
所述第一掺铒光纤放大器，用于对作为单边带信号的脉冲光进行功率放大；
所述第一带通滤波器，用于对放大后的脉冲光进行滤波处理，滤波后的脉冲光经过环形器注入传感光纤。
本实施例中，如图1所示，在发送端，将窄线宽激光器(laser)(1550.12nm，线宽<100Hz)发出的光通过99：1的第一耦合器(Coupler)分成两路，99％作为信号光被声光调制器(AOM)调制成脉冲光，1％作为本振光发送至第二耦合器；在AOM频移(例如，上频移)的作用下脉冲光成为单边带(SSB)信号，作为单边带信号的脉冲光经过第一掺铒光纤放大器(EDFA)进行功率放大，但是脉冲光功率不宜放大太多，以避免发生受激布里渊和受激拉曼散射；此时，EDFA由于放大引入的自发辐射噪声(ASE)不可忽视，需要1550.12nm波段的第一带通滤波器(BPF)进行滤波处理滤除带外噪声；之后，脉冲光从环形器1端口进入，从环形器2端口输出到传感光纤中，从传感光纤反射回的瑞利散射光通过环形器3端口输出；其中，由于是在实验，还需要在传感光纤某一位置放置压电陶瓷扰动器(PZT)来模拟扰动。
本实施例中，例如，可以选择单模裸光纤作为传感光纤。
在前述分布式相位敏感光时域反射仪传感系统的具体实施方式中，进一步地，所述系统还包括：第二掺铒光纤放大器和第二带通滤波器；
所述第二掺铒光纤放大器，用于对环形器输出的瑞利散射光进行功率放大；
所述第二带通滤波器，用于对放大后的瑞利散射光进行滤波处理，滤波后的瑞利散射光输入至第二耦合器。
本实施例中，如图1所示，在接收端，环形器输出的瑞利散射光通过第二掺铒光纤放大器(EDFA)进行功率放大和第二带通滤波器(BPF)进行滤波处理，滤波处理后的瑞利散射光与1％的本振光在一个50：50的2×1的第二耦合器中进行混频，混频之后的光信号经过单个光电探测器(PD)进行光电转换，然后，通过数字信号处理器对光电转换后的电信号进行带通滤波并取滤波后电信号中的一个边带，对提取的边带进行数字信号处理，得到光纤扰动信号的位置及相位差。
实施例二
如图2所示，本发明实施例提供的应用于分布式相位敏感光时域反射仪传感系统的相位提取方法，包括：
S101，在发送端，通过声光调制器将信号光调制成脉冲光，脉冲光经过环形器注入传感光纤，从传感光纤反射回的瑞利散射光通过环形器输出；其中，脉冲光在声光调制器频移的作用下成为单边带信号；
S102，在接收端，将本振光和环形器输出的瑞利散射光在第二耦合器中进行混频，其中，混频之后的光信号为单边带信号，混频之后的光信号通过单个光电探测器进行光电转换，对光电转换后的电信号进行带通滤波并取滤波后电信号中的一个边带，对提取的边带进行数字信号处理，得到光纤扰动信号的位置及相位。
本发明实施例所述的应用于分布式相位敏感光时域反射仪传感系统的相位提取方法，利用声光调制器的频移功能，使得经过声光调制器频移的单边带信号在光电探测器接收之后，与光电转换过程中的信号自拍频干扰(SSBI)不发生重叠，具有隔离于SSBI的天然保护带，以便数字信号处理器通过带通滤波滤出不含SSBI的接收信号，从而排除SSBI对接收信号的干扰；采用单光电探测器进行光电转换，能够降低光电转换过程中引入的噪声，在长探测距离、高扰动频率的场景中发挥出更加明显的优势；在数字信号处理部分只取电信号的一个边带即可实现光纤扰动信号的定位及相位的检测，这样，在提高扰动位置定位和扰动相位提取的准确性的同时，还能实现光纤扰动信号的定位及相位的快速检测。
本发明实施例所述的相位提取方法，有望成为未来分布式相位敏感光纤传感系统信噪比提升、简洁的算法处理流程且低成本的实施方案。
本实施例中，在分布式相位敏感光时域反射仪(Φ-OTDR)传感系统中声光调制器(AOM)具有两个作用：
1)将窄线宽激光器发射的连续的高相干光调制成脉冲光，其中，脉冲光的重复周期根据传感光纤的长度来设定，重复周期应大于脉冲光在传感光纤中传输时长的2倍，脉冲光的脉宽决定了Φ-OTDR系统的空间分辨率，设脉冲光的脉宽为τ，则系统空间分辨率可表示为
Figure GDA0002761014220000081
其中，c为光在真空中传输的速度，n为光在光纤中传输的折射率；
2)将调制后的脉冲光上频移或者下频移，例如，假设AOM上频移fb，则调制后的脉冲光信号与原来的脉冲光信号有fb的频差。
本实施例中，经过AOM调制以后的脉冲光信号用s0(t)来表示，其光谱如图3(a)所示，其中，B为调制信号的带宽，脉冲光信号s0(t)的光谱S0(j2πf)表达式为：
Figure GDA0002761014220000082
式(1)中，f为自变量，f表示频率，τ为脉冲光的脉宽，T为瑞利散射光的周期，fb为AOM的频移量，k为负无穷到正无穷范围内的整数,随着k取值的不同，S0(j2πf)也发生相应变化，从式(1)可以看出经AOM调制以后的信号中心频率位于fb处，信号的包络为Sa函数，内部由频率间隔为
Figure GDA0002761014220000083
的冲激函数δ(·)构成，如图3(a)所示。
本实施例中，将调制后的脉冲光通过放大、带通滤波之后经过环形器注入传感光纤，从传感光纤反射回来的瑞利散射光sRBS(t)的时域表达式为：
Figure GDA0002761014220000084
式(2)中，t为时间，f0为窄线宽激光器的频率，fb为声光调制器的频移量，As(t)为瑞利散射光的幅值，
Figure GDA0002761014220000091
为瑞利散射光的瞬时相位偏移。瑞利散射光的光谱如图3(b)所示，信号位于频率为fb位置处，调制信号的带宽为B。瑞利散射光的相位对扰动极其敏感，当传感光纤的某一处发生扰动时，该处的折射率发生变化，从而引起瑞利散射光相位的变化，通过检测瑞利散射光的相位变化来实现扰动信号的位置及功率分析。
本实施例中，本振光为载波，本振光sLO(t)的时域表达式为：
Figure GDA0002761014220000092
式(3)中，AL为本振光的幅值，
Figure GDA0002761014220000093
为本振光的瞬时相位偏移。
本实施例中，通过放大、带通滤波之后的瑞利散射光与窄线宽激光器发出的1％的本振光在50：50的第二耦合器中混频后，输出的光信号smix(t)的时域表达式为：
Figure GDA0002761014220000094
光信号smix(t)的光谱如图3(c)所示，通过观察此时光信号smix(t)的光谱，可知在光电转换之前的信号为SSB信号。在单PD对混频后的光信号smix(t)进行光电转换之后，得到的电信号I(t)的时域表达式为：
Figure GDA0002761014220000095
其中，AS 2(t)为信号自拍频干扰(SSBI)项，用图3(d)中位于基带的三角形表示其频谱，SSBI的带宽为B/2，B为调制信号的带宽。
本实施例中，使用的AOM频移一般为80MHz或200MHz，AOM频移之后的信号带宽一般不超过40MHz。由于AOM的频移作用，信号位于频率为fb的位置，fb通常为80MHz或200MHz，信号的带宽一般不超过40MHz，SSBI的带宽一般不超过20MHz，由于AOM的频移特性，经过AOM频移的SSB信号在PD接收之后，与光电转换过程中信号自拍频干扰(SSBI)不发生重叠，对信号形成了隔离于SSBI的天然保护带，这种保护带的存在使直接带通滤波就可滤出不含SSBI的接收信号，从而排除SSBI对接收信号的干扰。
因此，如图4所示，本实施例中的数字信号处理器首先需要对光电转换后得到的电信号进行带通滤波，由于声光调制器的移频作用,信号自拍频干扰项与信号项没有在同一带宽范围内,此处的滤波操作能过滤除信号自拍频干扰项,并取滤波后电信号的一个边带，得到瑞利散射光功率IPF(t)表达式：
Figure GDA0002761014220000101
其中，时间t＝t0+NT，N为提取的信号IPF(t)中有完整瑞利散射光周期的个数，T为瑞利散射光的周期，t0为提取的信号IPF(t)中第一个不完整周期瑞利散射光持续时间；
去除t0时间内的信号，得到t’＝NT时间内的信号，将t’时间内的信号转换为N×T的数据矩阵，其中，每行数据与传感光纤中随时间变化的瑞利散射光相对应；
以数据矩阵的行为单位，确定数据矩阵中行之间数据的功率差，通过功率差实现扰动位置定位；这是根据瑞利散射光的幅度和相位只与位置和光纤的折射率有关的特性，在无扰动的位置功率差为0，在有扰动的位置光纤的折射率发生变化，功率差不再为0，以此来实现扰动位置定位；
以数据矩阵的行为单位，确定数据矩阵中行之间数据(即：瑞利散射光)的相位差，瑞利散射光的相位差可用式(7)来表示：
Figure GDA0002761014220000102
其中，t1为第一个瑞利散射光周期中扰动位置对应的时间，IPF *(t1)为瑞利散射光功率IPF(t)表达式的共轭，根据式(7)取出的瑞利散射光相位差除包含瞬时相位差之外，还有残余的2πfbnT项，由光纤扰动信号引起的相位是在之前无扰动相位的基础之上变化的，所以在取扰动信号相位时需要扰动位置处瑞利散射光相位差与扰动之前无扰动位置瑞利散射光相位差做差，因此，设无扰动位置(参考位置)对应的时间为ta+nT，扰动位置对应的时间为tb+nT，其中，ta为第一个瑞利散射光周期内无扰动位置(参考位置)对应的时间，tb为第一个瑞利散射光周期内扰动位置对应的时间。根据式(7)，分别确定扰动位置和无扰动位置的瑞利散射光相位差，然后，将扰动位置和无扰动位置的瑞利散射光相位差做差之后，得到光纤扰动信号相位
Figure GDA0002761014220000103
Figure GDA0002761014220000104
式(8)中，
Figure GDA0002761014220000105
为本振光的瞬时相位变化，
Figure GDA0002761014220000106
为扰动信号引起的瑞利散射光的相位变化，在数字信号处理后可发现
Figure GDA0002761014220000107
与信号发生器施加给压电陶瓷振荡器(PZT)的信号一致，据此可通过快速傅里叶变换(FFT)，实现扰动信号的频率分析。
将扰动位置处瑞利散射光相位与扰动之前无扰动位置瑞利散射光相位做差，得到光纤扰动信号相位差。
需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。
以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明所述原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (2)

1.一种应用于分布式相位敏感光时域反射仪传感系统的相位提取方法，其特征在于，所述系统包括：窄线宽激光器、第一耦合器、声光调制器、环形器、第二耦合器、光电探测器和数字信号处理器；其中，
窄线宽激光器，用于发射激光并将发射的激光通过第一耦合器分成两路，其中，一路作为信号光被声光调制器调制成脉冲光并在声光调制器频移的作用下成为单边带信号，作为单边带信号的脉冲光经过环形器注入传感光纤，从传感光纤反射回的瑞利散射光通过环形器输出；另一路作为本振光发送至第二耦合器；
第二耦合器，用于对本振光和环形器输出的瑞利散射光进行混频，混频之后的光信号通过单个光电探测器进行光电转换；其中，混频之后的光信号为单边带信号；
数字信号处理器，用于对光电转换后的电信号进行带通滤波并取滤波后电信号中的一个边带，对提取的边带进行数字信号处理，得到光纤扰动信号的位置及相位；
所述方法包括：
在发送端，通过声光调制器将信号光调制成脉冲光，脉冲光经过环形器注入传感光纤，从传感光纤反射回的瑞利散射光通过环形器输出；其中，脉冲光在声光调制器频移的作用下成为单边带信号；
在接收端，将本振光和环形器输出的瑞利散射光在第二耦合器中进行混频，其中，混频之后的光信号为单边带信号，混频之后的光信号通过单个光电探测器进行光电转换，对光电转换后的电信号进行带通滤波并取滤波后电信号中的一个边带，对提取的边带进行数字信号处理，得到光纤扰动信号的位置及相位；
其中，所述脉冲光经过环形器注入传感光纤，从传感光纤反射回的瑞利散射光通过环形器输出包括：
将调制后的脉冲光通过放大、带通滤波之后经过环形器注入传感光纤，从传感光纤反射回的瑞利散射光sRBS(t)的时域表达式为：
Figure FDA0002761014210000011
其中，t为时间，As(t)为瑞利散射光的幅值，f0为窄线宽激光器的频率，fb为声光调制器的频移量，
Figure FDA0002761014210000021
为瑞利散射光的瞬时相位偏移；
其中，本振光为载波，本振光sLO(t)的时域表达式为：
Figure FDA0002761014210000022
其中，AL为本振光的幅值，
Figure FDA0002761014210000023
为本振光的瞬时相位偏移；
其中，本振光和环形器输出的经过放大、带通滤波之后的瑞利散射光在第二耦合器中混频后，输出的光信号smix(t)的时域表达式为：
Figure FDA0002761014210000024
其中，光电转换后得到的电信号I(t)的时域表达式为：
Figure FDA0002761014210000025
其中，AS 2(t)为信号自拍频干扰项；
其中，所述对光电转换后的电信号进行带通滤波并取滤波后电信号中的一个边带，对提取的边带进行数字信号处理，得到光纤扰动信号的位置及相位包括：
对光电转换后得到的电信号进行带通滤波并取滤波后电信号中的一个边带，得到瑞利散射光功率IPF(t)表达式：
Figure FDA0002761014210000026
其中，时间t＝t0+NT，N为提取的信号IPF(t)中有完整瑞利散射光周期的个数，T为瑞利散射光的周期，t0为提取的信号IPF(t)中第一个不完整周期瑞利散射光持续时间；
去除t0时间内的信号，得到t’＝NT时间内的信号，将t’时间内的信号转换为N×T的数据矩阵，其中，每行数据与传感光纤中随距离变化的瑞利散射光相对应；
以数据矩阵的行为单位，确定数据矩阵中行之间数据的功率差，通过功率差实现扰动位置定位；其中，无扰动的位置功率差为0，有扰动的位置功率差不为0，以此来实现扰动信号的定位；
以数据矩阵的行为单位，确定数据矩阵中行之间数据的相位差，将扰动位置处瑞利散射光相位差与扰动之前无扰动位置瑞利散射光相位差做差，得到光纤扰动信号相位。
2.根据权利要求1所述的相位提取方法，其特征在于，所述将扰动位置处瑞利散射光相位差与扰动之前无扰动位置瑞利散射光相位差做差，得到光纤扰动信号相位包括：
设无扰动位置对应的时间为ta+nT，扰动位置对应的时间为tb+nT，根据相位差公式：
Figure FDA0002761014210000031
分别确定扰动位置和无扰动位置的瑞利散射光相位差，并将扰动位置和无扰动位置的瑞利散射光相位差做差，得到光纤扰动信号相位
Figure FDA0002761014210000032
Figure FDA0002761014210000033
其中，t1为第一个瑞利散射光周期中扰动位置对应的时间，IPF *(t1)为瑞利散射光功率IPF(t)表达式的共轭，ta为第一个瑞利散射光周期内无扰动位置对应的时间，tb为第一个瑞利散射光周期内扰动位置对应的时间，n＝2…N,0<t1≤T，
Figure FDA0002761014210000034
为扰动信号引起的瑞利散射光的相位变化，
Figure FDA0002761014210000035
为本振光的瞬时相位变化。
CN202010071038.3A 2020-01-21 2020-01-21 分布式相位敏感光时域反射仪传感系统及相位提取方法 Active CN111238551B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010071038.3A CN111238551B (zh) 2020-01-21 2020-01-21 分布式相位敏感光时域反射仪传感系统及相位提取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010071038.3A CN111238551B (zh) 2020-01-21 2020-01-21 分布式相位敏感光时域反射仪传感系统及相位提取方法

Publications (2)

Publication Number Publication Date
CN111238551A CN111238551A (zh) 2020-06-05
CN111238551B true CN111238551B (zh) 2021-02-26

Family

ID=70871484

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010071038.3A Active CN111238551B (zh) 2020-01-21 2020-01-21 分布式相位敏感光时域反射仪传感系统及相位提取方法

Country Status (1)

Country Link
CN (1) CN111238551B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111609918A (zh) * 2020-06-09 2020-09-01 重庆大学 基于包络检波电路的光纤分布式振动传感系统
CN112611444B (zh) * 2020-12-30 2023-04-28 西安和其光电科技股份有限公司 一种可精确定位的分布式光纤振动监测系统及方法
CN113447772A (zh) * 2021-06-18 2021-09-28 安徽龙联智能光电有限公司 一种高压电缆局部放电在线监测系统及方法
CN113701793B (zh) * 2021-08-09 2022-05-06 南京大学 相位敏感型光时域反射系统中非硬件实现相位锁定的方法
CN113776567B (zh) * 2021-08-27 2024-06-14 武汉普赛斯电子技术有限公司 一种基于光频域反射计的偏振衰落消除装置及方法
CN114111860B (zh) * 2021-12-03 2022-08-30 北京科技大学 基于多频脉冲编码的分布式φ-otdr传感方法及系统
CN114964329B (zh) * 2022-05-22 2023-03-14 上海交通大学 双边带光频域反射仪
CN116295782B (zh) * 2023-03-08 2023-10-03 浙江信测通信股份有限公司 一种基于φ-otdr的分布式光纤振动传感系统及相位解调方法
CN116698096A (zh) * 2023-06-26 2023-09-05 国网江苏省电力有限公司南京供电分公司 一种基于双波长光源的光纤传感方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6856400B1 (en) * 2000-12-14 2005-02-15 Luna Technologies Apparatus and method for the complete characterization of optical devices including loss, birefringence and dispersion effects
CN102628698B (zh) * 2012-04-06 2015-02-18 中国科学院上海光学精密机械研究所 分布式光纤传感器及信息解调方法
CN204718622U (zh) * 2015-06-10 2015-10-21 贵州电网公司信息通信分公司 一种相干otdr装置
CN105509868B (zh) * 2015-12-16 2019-06-04 上海大学 相位敏感光时域反射光纤分布式传感系统相位计算方法
CN108683453B (zh) * 2018-04-18 2020-08-21 华中科技大学 一种基于合成单边带信号的线性直接探测方法与系统

Also Published As

Publication number Publication date
CN111238551A (zh) 2020-06-05

Similar Documents

Publication Publication Date Title
CN111238551B (zh) 分布式相位敏感光时域反射仪传感系统及相位提取方法
JP6698164B2 (ja) 周波数合成に基づいた光周波数領域反射方法及びシステム
US10756820B2 (en) Phase analog coherent demodulation system and method for phase sensitive optical time domain reflectometry
CN105490738B (zh) 基于频率合成的光频域反射方法及系统
CN107607135B (zh) 一种混沌布里渊光时域/相干域融合分析装置及方法
CN106768277B (zh) 一种分布式光纤振动传感装置的解调方法
WO2015010581A1 (zh) 一种基于探测频率编码的相干光时域反射仪
CN110470327A (zh) 一种光时域分析仪和分析方法
CN104792343A (zh) 一种单端结构动态测量的布里渊光纤传感系统与传感方法
CN102571200A (zh) 多频探测光相干光时域反射仪方法和装置
CN108827447B (zh) 一种异频双脉冲cotdr传感装置和方法
CN111721338B (zh) 一种泵浦光高低频交替调制的布里渊光时域分析系统
CN112697180B (zh) 一种温度和振动同时测量的融合型分布式光纤传感系统及方法
JP2019161246A (ja) デジタルコヒーレント伝送システム
CN205453695U (zh) 基于频率合成的光频域反射装置
CN112880711B (zh) 一种基于双脉冲调制的分布式光纤传感方法及系统
CN112118043B (zh) 一种基于psk脉冲编码的cotdr测量方法及其装置
WO2024098766A1 (zh) 一种同频共用光纤通信感知一体化系统
CN111637910B (zh) 时域差分高速混沌布里渊光相干域监测装置及方法
CN106656314B (zh) 一种监测光通信网络色散的方法及装置
CN109323750B (zh) 分布式光纤振动传感系统及相位解调方法
CN114623920B (zh) 一种φ-otdr型分布式光纤声波传感系统及信号解调方法
CN113390449A (zh) 一种光纤传感装置
JP2019060666A (ja) 光ファイバセンシングシステム及び光ファイバセンシング方法
JP3152314B2 (ja) 後方散乱光の測定方法およびその装置

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