CN113176075A - 一种损耗与振动同步监测的相干及偏振衰落抑制方法 - Google Patents

一种损耗与振动同步监测的相干及偏振衰落抑制方法 Download PDF

Info

Publication number
CN113176075A
CN113176075A CN202110453267.6A CN202110453267A CN113176075A CN 113176075 A CN113176075 A CN 113176075A CN 202110453267 A CN202110453267 A CN 202110453267A CN 113176075 A CN113176075 A CN 113176075A
Authority
CN
China
Prior art keywords
polarization state
amplitude
distribution matrix
phase
disturbance
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
CN202110453267.6A
Other languages
English (en)
Other versions
CN113176075B (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.)
Nanjing University
Original Assignee
Nanjing 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 Nanjing University filed Critical Nanjing University
Priority to CN202110453267.6A priority Critical patent/CN113176075B/zh
Publication of CN113176075A publication Critical patent/CN113176075A/zh
Priority to US17/908,248 priority patent/US12117361B2/en
Priority to PCT/CN2022/082419 priority patent/WO2022227948A1/zh
Application granted granted Critical
Publication of CN113176075B publication Critical patent/CN113176075B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
    • G01M11/30—Testing of optical devices, constituted by fibre optics or optical waveguides
    • G01M11/31—Testing of optical devices, constituted by fibre optics or optical waveguides with a light emitter and a light receiver being disposed at the same side of a fibre or waveguide end-face, e.g. reflectometers
    • G01M11/3109—Reflectometers detecting the back-scattered light in the time-domain, e.g. OTDR
    • G01M11/3145—Details of the optoelectronics or data analysis
    • 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/02—Testing optical properties
    • 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
    • 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
    • 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
    • G01H9/006—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means using fibre optic sensors the vibrations causing a variation in the relative position of the end of a fibre and another element
    • 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/3181—Reflectometers dealing with polarisation
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Optics & Photonics (AREA)
  • Engineering & Computer Science (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Optical Communication System (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明涉及一种损耗与振动同步监测的相干及偏振衰落抑制方法，在仅使用普通单模传感光纤的前提下，基于偏振分集的测量方式，解决了Φ‑OTDR和COTDR在偏振态问题上的冲突，并且对于损耗参量的监测，不仅能够获得更好的相干衰落噪声抑制效果，还能够增强微小损耗事件的识别能力，同时对于扰动参量的监测，能够抑制相干衰落对Φ‑OTDR相位解调的影响，高保真地重构振动信号，从而有效地降低外部扰动预警的误报率。

Description

一种损耗与振动同步监测的相干及偏振衰落抑制方法
技术领域
本发明涉及一种损耗与振动同步监测的相干及偏振衰落抑制方法，属于光纤传感技术领域。
背景技术
相干光时域反射仪(Coherent Optical Time Domain Reflectometry，COTDR)采用相干探测将瑞利背向散射信号(Rayleigh Backscattering Signal，RBS)光学信号转换为某一特定中频的电学信号。通过对该中频信号进行电学的窄带滤波，抑制绝大部分的ASE噪声，从而使其能够监测多中继超长距离光纤通信线路。近年来，COTDR技术仍有持续改进，主要依托频分复用技术提高监测通道的带宽利用效率，缩短测量耗时，但就其目标参量而言仍局限在对光纤损耗参量的监测。这就造成基于COTDR的光纤链路监测系统仅能够进行事后的故障诊断与定位。事实上，在出现真正的不可逆转的损伤之前，虽然损耗值尚未出现明显的变化，但线路周边环境中的潜在威胁常通过直接接触或间接的振动信号传递，对光缆形成扰动，从而为提前预警提供依据。在扰动监测方面，相位敏感型光时域反射仪(Phase-sensitive Optical Time-Domain Reflectometry，Φ-OTDR)凭借灵敏度高和响应速度快的优点成为现阶段海缆监测技术的研究热点。
考虑到Φ-OTDR和COTDR都使用线宽窄、频率漂移小的激光器作为光源，并且在结构上外差探测Φ-OTDR与COTDR系统也具有高度的相似性。为了满足光纤监测系统对振动和损耗同步监测的需求，可以考虑将Φ--OTDR与COTDR技术进行结合。但是两种测量原理的融合，存在诸多矛盾和问题，例如：(1)Φ-OTDR系统为了能够获得较高的监测精度，相邻RBS曲线之间必须具备较高的相似性，因此需要一个频率稳定的光源，而COTDR则需要进行主动扫频和多点数字平均技术来抑制相干衰落噪声，产生对光源工作模式的需求冲突；(2)对于Φ-OTDR而言，偏振态的改变会在检测信号中引入随时间变化的偏振分量，影响振动测量精度。而COTDR中常采用扰偏器对光本振光(Optical Local Oscillator，OLO)光的偏振态进行快速扰动，以抑制由于RBS与OLO偏振态不匹配导致的偏振衰落噪声，这与Φ-OTDR希望偏振态高度稳定的要求相矛盾；(3)目前Φ-OTDR解调方法有两种：鉴幅型和鉴相型。鉴幅型解调方法根据RBS的幅度信息快速定位振动事件。由于RBS的幅度变化与振动事件引起的光纤动态应变没有明确的线性关系，因此鉴幅型解调方法通常仅能够做定性测量，无法高保真地重构扰动事件的时频特征。鉴相型解调方法利用RBS的相位变化与光纤局部区域的动态应变之间的线性关系，可以精准还原扰动事件，提升后续模式识别的准确性。但是由于相干衰落，光纤某些区域强度趋近于零，在这些区域相位解调结果很容易出现信噪比的急剧恶化导致频繁的误报警。
发明内容
本发明所要解决的技术问题是提供一种损耗与振动同步监测的相干及偏振衰落抑制方法，采用全新设计方法，在实现光纤链路损耗监测的同时，能够实现扰动事件的预警。
本发明为了解决上述技术问题采用以下技术方案：本发明设计了一种损耗与振动同步监测的相干及偏振衰落抑制方法，基于针对待测光纤采用偏振分集探测方式获得的两路正交偏振态中频信号，执行如下步骤A至步骤F：
步骤A.分别针对两路正交偏振态中频信号中的P偏振态中频信号与S偏振态中频信号，执行IQ解调，获得P偏振态中频信号所对应的P偏振态幅度分布矩阵、P偏振态相位分布矩阵，以及S偏振态中频信号所对应的S偏振态幅度分布矩阵、S偏振态相位分布矩阵，然后进入步骤B；其中，幅度分布矩阵表示信号幅度随时间、距离的分布，相位分布矩阵表示信号相位随时间、距离的分布；
步骤B.根据P偏振态幅度分布矩阵和S偏振态幅度分布矩阵，分别获得P偏振态扰动位置和S偏振态扰动位置，然后进入步骤C；
步骤C.分别针对P偏振态幅度分布矩阵和S偏振态幅度分布矩阵，按预设大小时域窗口TWsize，沿时间轴进行划分，获得P偏振态幅度分布矩阵所对应的各个TWsize时长的信号区间，以及S偏振态幅度分布矩阵所对应的各个TWsize时长的信号区间，然后进入步骤D；
步骤D.分别针对P偏振态幅度分布矩阵和S偏振态幅度分布矩阵划分后的各信号区间，计算P偏振态和S偏振态分别对应当前信号区间扰动位置前空间分辨率宽度内的第一幅度特征值、以及所在位置；同时，计算P偏振态和S偏振态分别对应当前信号区间扰动位置后空间分辨率宽度内的第二幅度特征值、以及所在位置；进而分别获得P偏振态与S偏振态各个信号区间扰动位置前和扰动位置后空间分辨率宽度内的第一幅度特征值和第二幅度特征值，以及该第一幅度特征值和第二幅度特征值所对应的位置，然后进入步骤E；
步骤E.基于P偏振态相位分布矩阵与S偏振态相位分布矩阵，分别针对P偏振态幅度分布矩阵和S偏振态幅度分布矩阵划分后的各信号区间，首先选择P偏振态和S偏振态对应当前信号区间扰动位置前空间分辨率宽度内较大的第一幅度特征值所对应的偏振态相位、以及所在位置，作为当前信号区间扰动位置前的偏振态相位；同时，选择P偏振态和S偏振态对应当前信号区间扰动位置后空间分辨率宽度内较大的第二幅度特征值所对应的偏振态相位、以及所在位置，作为当前信号区间扰动位置后的偏振态相位；然后对当前信号区间扰动位置前后所选的偏振态相位进行相位差运算，实现相位提取，获得当前信号区间扰动位置的相位信号；进而获得各个信号区间扰动位置的相位信号，实现扰动位置相位信号重构，然后进入步骤F；
步骤F.根据各信号区间扰动位置的重构相位信号，执行相位解缠绕，还原由外部扰动信号所引起的扰动位置处的光纤伸缩情况，即获得待测光纤的扰动监测结果。
作为本发明的一种优选技术方案，基于步骤A所获P偏振态幅度分布矩阵与S偏振态幅度分布矩阵，还包括执行如下步骤I至步骤II，实现对待测光纤损耗监测结果的获得：
步骤I.沿时间轴顺序，分别对P偏振态幅度分布矩阵和S偏振态幅度分布矩阵进行平方和处理，获得P偏振态幅度分布矩阵平方和结果，以及S偏振态幅度分布矩阵平方和结果，然后进入步骤II；
步骤II.针对P偏振态幅度分布矩阵平方和结果与S偏振态幅度分布矩阵平方和结果，进行累加平均处理，获得待测光纤的损耗监测结果。
作为本发明的一种优选技术方案，所述步骤D包括如下：
分别针对P偏振态幅度分布矩阵和S偏振态幅度分布矩阵划分后的各信号区间，执行步骤D1至步骤D2，获得P偏振态与S偏振态各信号区间扰动位置前和扰动位置后空间分辨率宽度内的第一幅度特征值和第二幅度特征值，以及该第一幅度特征值和第二幅度特征值所对应的位置，然后进入步骤E；
步骤D1.根据P偏振态幅度分布矩阵和S偏振态幅度分布矩阵，分别获得当前信号区间P偏振态和S偏振态扰动位置前、空间分辨率宽度范围内各空间位置的最小幅度值，并选取P偏振态和S偏振态分别所对应该各最小幅度值中的最大幅度值，作为P偏振态和S偏振态分别对应当前信号区间各自扰动位置前所对应的第一幅度特征值，分别记为Ap ra和As ra,其中下标ra表示P偏振态、S偏振态分别对应信号区间各自扰动位置前第一幅度特征值所在位置；
步骤D2.根据P偏振态幅度分布矩阵和S偏振态幅度分布矩阵，分别获得当前信号区间P偏振态和S偏振态扰动位置后、空间分辨率宽度范围内各空间位置的最小幅度值，并选取P偏振态和S偏振态分别所对应该各最小幅度值的最大幅度值，作为P偏振态和S偏振态分别对应当前信号区间各自扰动位置后所对应的第二幅度特征值，分别记为Ap rb和As rb,其中下标rb表示P偏振态、S偏振态分别对应信号区间各自扰动位置后第二幅度特征值所在位置。
作为本发明的一种优选技术方案，所述步骤E中各信号区间扰动位置前、后的偏振态相位选择包括：
(1)信号区间扰动位置前选择了P偏振态相位，扰动位置后选择了P偏振态相位；
(2)信号区间扰动位置前选择了P偏振态相位，扰动位置后选择了S偏振态相位；
(3)信号区间扰动位置前选择了S偏振态相位，扰动位置后选择了P偏振态相位；
(4)信号区间扰动位置前选择了S偏振态相位，扰动位置后选择了S偏振态相位。
本发明所述一种损耗与振动同步监测的相干及偏振衰落抑制方法，采用以上技术方案与现有技术相比，具有以下技术效果：
本发明所设计一种损耗与振动同步监测的相干及偏振衰落抑制方法，在仅使用普通单模传感光纤的前提下，基于偏振分集的测量方式，解决了Φ-OTDR和COTDR在偏振态问题上的冲突，并且对于损耗参量的监测，不仅能够获得更好的相干衰落噪声抑制效果，还能够增强微小损耗事件的识别能力，同时对于扰动参量的监测，能够抑制相干衰落对Φ-OTDR相位解调的影响，高保真地重构振动信号，从而有效地降低外部扰动预警的误报率。
附图说明
图1是本发明所设计损耗与振动同步监测的相干及偏振衰落抑制方法的系统架构示意图；
图2为本发明所设计中待测光纤连接示意图；
图3是本发明所设计损耗与振动同步监测的相干及偏振衰落抑制方法流程图；
图4是实施例中P态PZT附近2s间的幅度瀑布图；
图5是实施例中216条曲线损耗测试结果；
图6是实施例中损耗测试结果局部放大图；
图7是实施例中两态提取的相位信息及经过优选后的重建结果。
其中，1.待测光纤1，2.待测光纤2，3.实验用圆筒形柱状压电陶瓷。
具体实施方式
下面结合说明书附图对本发明的具体实施方式作进一步详细的说明。
本发明设计了一种损耗与振动同步监测的相干及偏振衰落抑制方法，该方法是基于具体所设计的系统架构去实现执行的，实际应用设计中，如图1所示，该系统具体包括激光器、光耦合器、声光调制器、光放大器、环形器、预放型偏振分集接收器(Pre-PDR)、偏振控制器、控制及处理模块、数据采集卡(DAQ)、脉冲信号发生器、以及P偏振态信号滤波装置和S偏振态滤波装置。
其中，激光器与控制及处理模块相连，激光器基于来自控制及处理模块的预设扫频速度进行扫频，产生连续光信号；激光器的光输出端对接光耦合器的输入端，由光耦合器针对所接收的连续光信号进行分束，分别获得探测光信号与本振光信号；光耦合器的本振光信号输出端对接偏振控制器的输入端，由偏振控制器针对所接收的本振光信号的偏振态进行调整；光耦合器的探测光信号输出端对接声光调制器的输入端，由声光调制器在脉冲信号发生器所输出脉冲调制信号的控制下、针对所接收到的探测光信号进行调制，获得探测光脉冲，并进行移频。
声光调制器的输出端对接光放大器的输入端，由声光调制器将探测光脉冲输送至光放大器进行放大处理；光放大器的输出端对接环形器，由光放大器将探测光脉冲输送至环形器，由环形器将探测光脉冲传输至待测光纤，并接收待测光纤返回的RBS信号；预放型偏振分集接收器(Pre-PDR)分别与环形器、偏振控制器相连接，由预放型偏振分集接收器(Pre-PDR)分别接收来自环形器的RBS信号、以及来自偏振控制器的本振光信号，并实现RBS信号中各偏振态与本振光信号的相干处理，分别获得P偏振态中频信号和S偏振态中频信号进行输出。
P偏振态信号滤波装置和S偏振态信号滤波装置结构相同，分别针对来自预放型偏振分集接收器(Pre-PDR)输出端的P偏振态中频信号和S偏振态中频信号进行滤波-放大-再滤波处理，实现偏振态中频信号的更新、并输出。
P偏振态信号滤波装置和S偏振态信号滤波装置的输出端对接数据采集卡(DAQ)的两路输入端，由数据采集卡(DAQ)对来自两路偏振态信号滤波装置输出更新后的偏振态中频信号进行采集，并传输至控制及处理模块进行数据处理，执行所述一种损耗与振动同步监测的相干及偏振衰落抑制方法；脉冲信号发生器基于与数据采集卡(DAQ)所输出各偏振态中频信号相同步的时钟信号，生成用于驱动声光调制器产生探测脉冲光的脉冲调制信号，并输送至声光调制器。
实际设计应用中，该系统中的两路偏振态信号滤波装置分别均包括第一带通滤波器、低噪声电学放大器、第二带通滤波器，两路偏振态信号滤波装置结构中，第一带通滤波器的输入端构成偏振态信号滤波装置的输入端，第一带通滤波器的输出端对接对应低噪声电学放大器的输入端，低噪声电学放大器的输出端对接对应第二带通滤波器的输入端，第二带通滤波器的输出端构成对应偏振态信号滤波装置的输出端；第一带通滤波器用于滤除所接收偏振态中频信号中的无用信号、以及宽带噪声，执行第一级滤波处理，实现对偏振态中频信号的更新；低噪声电学放大器用于针对第一带通滤波器输出的偏振态中频信号，实现放大更新；第二带通滤波器用于针对对应低噪声电学放大器输出的偏振态中频信号，滤除其中的无用信号、以及宽带噪声，执行第二级滤波处理，实现对偏振态中频信号的更新输出。
并且在实际的系统搭建中，针对光放大器，具体选用为掺铒光纤放大器；针对激光器，具体选用窄线宽可调激光器。
基于上述实际应用所搭建的系统，具体执行本发明所设计损耗与振动同步监测的相干及偏振衰落抑制方法，基于针对待测光纤采用偏振分集探测方式获得的两路正交偏振态中频信号，实际应用中，如图3所示，具体执行如下步骤A至步骤F。
步骤A.分别针对两路正交偏振态中频信号中的P偏振态中频信号与S偏振态中频信号，执行IQ解调，获得P偏振态中频信号所对应的P偏振态幅度分布矩阵、P偏振态相位分布矩阵，以及S偏振态中频信号所对应的S偏振态幅度分布矩阵、S偏振态相位分布矩阵，然后进入步骤B；其中，幅度分布矩阵表示信号幅度随时间、距离的分布，相位分布矩阵表示信号相位随时间、距离的分布。
步骤B.根据P偏振态幅度分布矩阵和S偏振态幅度分布矩阵，分别获得P偏振态扰动位置和S偏振态扰动位置，然后进入步骤C。
步骤C.分别针对P偏振态幅度分布矩阵和S偏振态幅度分布矩阵，按预设大小时域窗口TWsize，沿时间轴进行划分，获得P偏振态幅度分布矩阵所对应的各个TWsize时长的信号区间，以及S偏振态幅度分布矩阵所对应的各个TWsize时长的信号区间，然后进入步骤D。
步骤D.分别针对P偏振态幅度分布矩阵和S偏振态幅度分布矩阵划分后的各信号区间，计算P偏振态和S偏振态分别对应当前信号区间扰动位置前空间分辨率宽度内的第一幅度特征值、以及所在位置；同时，计算P偏振态和S偏振态分别对应当前信号区间扰动位置后空间分辨率宽度内的第二幅度特征值、以及所在位置；进而分别获得P偏振态与S偏振态各个信号区间扰动位置前和扰动位置后空间分辨率宽度内的第一幅度特征值和第二幅度特征值，以及该第一幅度特征值和第二幅度特征值所对应的位置，然后进入步骤E。
实际应用当中，上述步骤D具体设计，分别针对P偏振态幅度分布矩阵和S偏振态幅度分布矩阵划分后的各信号区间，执行步骤D1至步骤D2，获得P偏振态与S偏振态各信号区间扰动位置前和扰动位置后空间分辨率宽度内的第一幅度特征值和第二幅度特征值，以及该第一幅度特征值和第二幅度特征值所对应的位置，然后进入步骤E。
步骤D1.根据P偏振态幅度分布矩阵和S偏振态幅度分布矩阵，分别获得当前信号区间P偏振态和S偏振态扰动位置前、空间分辨率宽度范围内各空间位置的最小幅度值，并选取P偏振态和S偏振态分别所对应该各最小幅度值中的最大幅度值，作为P偏振态和S偏振态分别对应当前信号区间各自扰动位置前所对应的第一幅度特征值，分别记为Ap ra和As ra,其中下标ra表示P偏振态、S偏振态分别对应信号区间各自扰动位置前第一幅度特征值所在位置。
步骤D2.根据P偏振态幅度分布矩阵和S偏振态幅度分布矩阵，分别获得当前信号区间P偏振态和S偏振态扰动位置后、空间分辨率宽度范围内各空间位置的最小幅度值，并选取P偏振态和S偏振态分别所对应该各最小幅度值的最大幅度值，作为P偏振态和S偏振态分别对应当前信号区间各自扰动位置后所对应的第二幅度特征值，分别记为Ap rb和As rb,其中下标rb表示P偏振态、S偏振态分别对应信号区间各自扰动位置后第二幅度特征值所在位置。
步骤E.基于P偏振态相位分布矩阵与S偏振态相位分布矩阵，分别针对P偏振态幅度分布矩阵和S偏振态幅度分布矩阵划分后的各信号区间，首先选择P偏振态和S偏振态对应当前信号区间扰动位置前空间分辨率宽度内较大的第一幅度特征值(max{Ap ra,As ra})所对应的偏振态相位、以及所在位置，作为当前信号区间扰动位置前的偏振态相位；同时，选择P偏振态和S偏振态对应当前信号区间扰动位置后空间分辨率宽度内较大的第二幅度特征值(max{Ap rb,As rb})所对应的偏振态相位、以及所在位置，作为当前信号区间扰动位置后的偏振态相位；然后对当前信号区间扰动位置前后所选的偏振态相位进行相位差运算，实现相位提取，获得当前信号区间扰动位置的相位信号；进而获得各个信号区间扰动位置的相位信号，实现扰动位置相位信号重构，然后进入步骤F。
实际应用中，步骤E中各信号区间扰动位置前、后的偏振态相位选择包括：
(1)信号区间扰动位置前选择了P偏振态相位，扰动位置后选择了P偏振态相位；
(2)信号区间扰动位置前选择了P偏振态相位，扰动位置后选择了S偏振态相位；
(3)信号区间扰动位置前选择了S偏振态相位，扰动位置后选择了P偏振态相位；
(4)信号区间扰动位置前选择了S偏振态相位，扰动位置后选择了S偏振态相位。
为抑制相干及偏振衰落提供了多种选择，从而实现更好的相干及偏振衰落抑制效果。
步骤F.根据各信号区间扰动位置的重构相位信号，执行相位解缠绕，还原由外部扰动信号所引起的扰动位置处的光纤伸缩情况，即获得待测光纤的扰动监测结果。
实际应用中，基于步骤A所获P偏振态幅度分布矩阵与S偏振态幅度分布矩阵，还包括执行如下步骤I至步骤II，实现对待测光纤损耗监测结果的获得：
步骤I.沿时间轴顺序，分别对P偏振态幅度分布矩阵和S偏振态幅度分布矩阵进行平方和处理，获得P偏振态幅度分布矩阵平方和结果，以及S偏振态幅度分布矩阵平方和结果，然后进入步骤II；
步骤II.针对P偏振态幅度分布矩阵平方和结果与S偏振态幅度分布矩阵平方和结果，进行累加平均处理，获得待测光纤的损耗监测结果。
将本发明所设计损耗与振动同步监测的相干及偏振衰落抑制方法，应用于实际当中，以长距离海缆监测为例，结合附图对本发明的技术方案做进一步的详细说明：
如图1所示，系统采用波长为1550.12nm、线宽为3.7kHz的窄线宽可调谐NKT激光器作为光源，通过控制及处理模块控制的可编程驱动器对其光频率进行调谐。NKT激光器输出的连续光被90:10的光耦合器分成两部分，其中90％的一路作为探测光，10％的一路作为本振光。探测光经声光调制器被调制为200MHz频移的光脉冲，再利用光放大器进行功率放大。针对长距离海缆监测海上定位的需求，100m以内的空间分辨率足以满足要求。因此实验中声光调制器的调制脉冲采用1us脉宽，周期为1ms。为了确保Φ-OTDR对外部扰动的测量精度，同时抑制RBS曲线的幅度波动，根据Φ-OTDR对相邻RBS曲线相关度的要求，以及扫频速度与调制信号脉宽、周期之间的关系，实验中选取了5MHz/s的扫频速度，这足以有效地抑制相干衰落对损耗测量精度的影响，而不牺牲Φ-OTDR的扰动事件识别能力。经光放大器放大后的探测脉冲光通过环形器的2端口注入待测光纤。
基于图2所示，待测光纤1和待测光纤2分别为长度约71km和1km的单模光纤。在实验用圆筒形柱状压电陶瓷(PZT)表面缠绕了约30m光纤，用于模拟产生外界声场的振动。实验中，对PZT施加了30Hz的正弦波模拟传感光纤受到外部声场扰动。
光纤中产生的RBS光经环形器的2端口返回，并由3端口输出后与本振光被预放型偏振分集相干接收器(Pre-PDR)接收，其中本振光一路通过偏振控制器实现对其偏振态的调整，以均衡Pre-PDR输出的P态和S态的中频信号强度，使其平均强度尽量一致。相干拍频产生P态和S态两路中频信号，再分别经中心频率为200MHz的带通滤波器和低噪声放大器进一步滤波和放大，最终由数据采集卡(DAQ)进行采集，传输至控制及处理模块进行数据处理。系统采样率为1.25GSa/s。其中，第一带通滤波器和第二带通滤波器的带宽通带范围为195MHz～205MHz，低噪声放大器的增益为27dB，该偏振态信号滤波装置的作用主要是为了滤除无用的信号，同时放大光电转换产生的微弱电信号。
采集后的数据处理如图3所示，具体过程如下：
(1)分别对采集到的两路中频信号(P态和S态)分别进行IQ解调，获得两态中频信号的幅度分布矩阵Ap m×n和As m×n，以及相位分布矩阵Φp m×n和Φs m×n，其中，A和Φ均为由时间行和距离列组成的m×n列的矩阵，二者中每一行分别代表着幅度和相位随距离的分布，每一列分别代表着幅度和相位随时间的变化。
(2)根据P态幅度分布矩阵和S态幅度分布矩阵，分别获得P偏振态扰动位置和S偏振态扰动位置，如图4所示，为P态PZT附近2s间的幅度瀑布图，图中虚线区域可以清楚地看到约71.9km处检测到了一个周期性的信号。
(3)分别对P态幅度分布矩阵和S态幅度分布矩阵，按预设大小时域窗口TWsize，沿时间轴进行时域划分。考虑到海底扰动事件的低频特性，这里设置时域窗口TWsize选择为20个点，即0.02s。将两态幅度分布矩阵分别划分为首尾相接且等宽度的428个信号区间。
(4)沿被划分的428个信号区间的时间轴，比较两态幅度分布矩阵在当前时域宽度TWsize下扰动位置前的第一幅度特征值和扰动位置后的第二幅度特征值，选择两态当前信号区间扰动位置前空间分辨率宽度内较大的第一幅度特征值(max{Ap ra,As ra})所对应的偏振态相位及所在位置ra，作为扰动位置前的偏振态相位；同时，选择两态当前信号区间扰动位置后空间分辨率宽度内较大的第二幅度特征值(max{Ap rb,As rb})所对应的偏振态相位及所在位置rb，作为扰动位置后的偏振态相位。对当前信号区间扰动位置前后所选的偏振态相位进行相位差运算，实现相位提取，获得当前信号区间扰动位置的相位信号。进而获得428个信号区间扰动位置的相位信号，实现扰动位置相位信号重构。
(5)对重构的相位信号进行解缠绕，还原由外部扰动信号引起的扰动位置处的光纤伸缩情况，获得扰动监测结果，从而实现对外部扰动信号的感测。
(6)分别对两态幅度矩阵Ap m×n和As m×n沿时间轴进行平方，得到两态幅度矩阵的平方和AP和AS。
(7)对两态平方和矩阵AP和AS再进行累加平均，得到AT曲线，从而获得损耗监测结果，如图5所示。其中，图5中的(a)和(b)分别为P偏振态和S偏振态各自获得的216条RBS曲线进行平均的结果，图5中的(c)为两态RBS曲线总条数平均的改进结果，可以看到将两态的RBS曲线进一步平均后，得到的COTDR曲线更加平滑，获得了59.2dB的总动态范围，对应的单程动态范围(Single-Way Dynamic Range，SWDR)为29.6dB。图6为图5中30km～40km的局部放大图，经计算P态和S态获得的COTDR曲线的波动分别为0.24dB和0.35dB，将两态累加平均获得COTDR曲线的波动为0.09B。可见，利用本发明提供的一种损耗与振动同步监测的相干及偏振衰落抑制方法，不仅能够获得更好的相干衰落噪声抑制效果，还能够增强微小损耗事件的识别能力。
如图7所示，展示了经过本发明提供的方法优选前后的相位提取结果。为了能够凸显偏振分集幅度优选的效果，在已知施加信号频率的基础上，我们对各自鉴相结果进行带通滤波。由上至下，图7中(a)和(b)分别为P态和S态各自RBS的重构振动信号波形图，可以看到未进行本发明的优选算法之前，两态RBS提取的相位信号均存在多处畸变；图7中(c)展示了上文P态和S态RBS信号解调得到的幅度在TWsize时域内ra和rb位置的幅度总和。与相应的信号强度相比，图7中(a)和(b)中的失真位置均为幅值较弱的位置。图7中(d)对应了图7中(c)的幅值选择策略，在当前时域窗口内选择最大总幅度特征值的一方进行相位提取。其中，蓝色域表示选择了P态，红色域表示选择了S态，从图中可以看到，随着时间的推移，其相位提取一直在两个偏振态中跳转，从而获得了如图7中(e)所示的重建相位信号波形图。
可见，利用本发明提供的一种损耗与振动同步监测的相干及偏振衰落抑制方法，能够高保真地重构振动信号，并且同步实现了单端传感距离至少71.9km的损耗与振动监测。
上面结合附图对本发明的实施方式作了详细说明，但是本发明并不限于上述实施方式，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下做出各种变化。

Claims (4)

1.一种损耗与振动同步监测的相干及偏振衰落抑制方法，其特征在于：基于针对待测光纤采用偏振分集探测方式获得的两路正交偏振态中频信号，执行如下步骤A至步骤F：
步骤A.分别针对两路正交偏振态中频信号中的P偏振态中频信号与S偏振态中频信号，执行IQ解调，获得P偏振态中频信号所对应的P偏振态幅度分布矩阵、P偏振态相位分布矩阵，以及S偏振态中频信号所对应的S偏振态幅度分布矩阵、S偏振态相位分布矩阵，然后进入步骤B；其中，幅度分布矩阵表示信号幅度随时间、距离的分布，相位分布矩阵表示信号相位随时间、距离的分布；
步骤B.根据P偏振态幅度分布矩阵和S偏振态幅度分布矩阵，分别获得P偏振态扰动位置和S偏振态扰动位置，然后进入步骤C；
步骤C.分别针对P偏振态幅度分布矩阵和S偏振态幅度分布矩阵，按预设大小时域窗口TWsize，沿时间轴进行划分，获得P偏振态幅度分布矩阵所对应的各个TWsize时长的信号区间，以及S偏振态幅度分布矩阵所对应的各个TWsize时长的信号区间，然后进入步骤D；
步骤D.分别针对P偏振态幅度分布矩阵和S偏振态幅度分布矩阵划分后的各信号区间，计算P偏振态和S偏振态分别对应当前信号区间扰动位置前空间分辨率宽度内的第一幅度特征值、以及所在位置；同时，计算P偏振态和S偏振态分别对应当前信号区间扰动位置后空间分辨率宽度内的第二幅度特征值、以及所在位置；进而分别获得P偏振态与S偏振态各个信号区间扰动位置前和扰动位置后空间分辨率宽度内的第一幅度特征值和第二幅度特征值，以及该第一幅度特征值和第二幅度特征值所对应的位置，然后进入步骤E；
步骤E.基于P偏振态相位分布矩阵与S偏振态相位分布矩阵，分别针对P偏振态幅度分布矩阵和S偏振态幅度分布矩阵划分后的各信号区间，首先选择P偏振态和S偏振态对应当前信号区间扰动位置前空间分辨率宽度内较大的第一幅度特征值所对应的偏振态相位、以及所在位置，作为当前信号区间扰动位置前的偏振态相位；同时，选择P偏振态和S偏振态对应当前信号区间扰动位置后空间分辨率宽度内较大的第二幅度特征值所对应的偏振态相位、以及所在位置，作为当前信号区间扰动位置后的偏振态相位；然后对当前信号区间扰动位置前后所选的偏振态相位进行相位差运算，实现相位提取，获得当前信号区间扰动位置的相位信号；进而获得各个信号区间扰动位置的相位信号，实现扰动位置相位信号重构，然后进入步骤F；
步骤F.根据各信号区间扰动位置的重构相位信号，执行相位解缠绕，还原由外部扰动信号所引起的扰动位置处的光纤伸缩情况，即获得待测光纤的扰动监测结果。
2.根据权利要求1所述一种损耗与振动同步监测的相干及偏振衰落抑制方法，其特征在于：基于步骤A所获P偏振态幅度分布矩阵与S偏振态幅度分布矩阵，还包括执行如下步骤I至步骤II，实现对待测光纤损耗监测结果的获得：
步骤I.沿时间轴顺序，分别对P偏振态幅度分布矩阵和S偏振态幅度分布矩阵进行平方和处理，获得P偏振态幅度分布矩阵平方和结果，以及S偏振态幅度分布矩阵平方和结果，然后进入步骤II；
步骤II.针对P偏振态幅度分布矩阵平方和结果与S偏振态幅度分布矩阵平方和结果，进行累加平均处理，获得待测光纤的损耗监测结果。
3.根据权利要求1或2所述一种损耗与振动同步监测的相干及偏振衰落抑制方法，其特征在于：所述步骤D包括如下：
分别针对P偏振态幅度分布矩阵和S偏振态幅度分布矩阵划分后的各信号区间，执行步骤D1至步骤D2，获得P偏振态与S偏振态各信号区间扰动位置前和扰动位置后空间分辨率宽度内的第一幅度特征值和第二幅度特征值，以及该第一幅度特征值和第二幅度特征值所对应的位置，然后进入步骤E；
步骤D1.根据P偏振态幅度分布矩阵和S偏振态幅度分布矩阵，分别获得当前信号区间P偏振态和S偏振态扰动位置前、空间分辨率宽度范围内各空间位置的最小幅度值，并选取P偏振态和S偏振态分别所对应该各最小幅度值中的最大幅度值，作为P偏振态和S偏振态分别对应当前信号区间各自扰动位置前所对应的第一幅度特征值，分别记为Ap ra和As ra,其中下标ra表示P偏振态、S偏振态分别对应信号区间各自扰动位置前第一幅度特征值所在位置；
步骤D2.根据P偏振态幅度分布矩阵和S偏振态幅度分布矩阵，分别获得当前信号区间P偏振态和S偏振态扰动位置后、空间分辨率宽度范围内各空间位置的最小幅度值，并选取P偏振态和S偏振态分别所对应该各最小幅度值的最大幅度值，作为P偏振态和S偏振态分别对应当前信号区间各自扰动位置后所对应的第二幅度特征值，分别记为Ap rb和As rb,其中下标rb表示P偏振态、S偏振态分别对应信号区间各自扰动位置后第二幅度特征值所在位置。
4.根据权利要求1所述一种损耗与振动同步监测的相干及偏振衰落抑制方法，其特征在于，所述步骤E中各信号区间扰动位置前、后的偏振态相位选择包括：
(1)信号区间扰动位置前选择了P偏振态相位，扰动位置后选择了P偏振态相位；
(2)信号区间扰动位置前选择了P偏振态相位，扰动位置后选择了S偏振态相位；
(3)信号区间扰动位置前选择了S偏振态相位，扰动位置后选择了P偏振态相位；
(4)信号区间扰动位置前选择了S偏振态相位，扰动位置后选择了S偏振态相位。
CN202110453267.6A 2021-04-26 2021-04-26 一种损耗与振动同步监测的相干及偏振衰落抑制方法 Active CN113176075B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN202110453267.6A CN113176075B (zh) 2021-04-26 2021-04-26 一种损耗与振动同步监测的相干及偏振衰落抑制方法
US17/908,248 US12117361B2 (en) 2021-04-26 2022-03-23 Method for suppressing coherent and polarization-induced fading by simultaneous monitoring of loss and vibration
PCT/CN2022/082419 WO2022227948A1 (zh) 2021-04-26 2022-03-23 一种损耗与振动同步监测的相干及偏振衰落抑制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110453267.6A CN113176075B (zh) 2021-04-26 2021-04-26 一种损耗与振动同步监测的相干及偏振衰落抑制方法

Publications (2)

Publication Number Publication Date
CN113176075A true CN113176075A (zh) 2021-07-27
CN113176075B CN113176075B (zh) 2022-03-25

Family

ID=76926187

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110453267.6A Active CN113176075B (zh) 2021-04-26 2021-04-26 一种损耗与振动同步监测的相干及偏振衰落抑制方法

Country Status (3)

Country Link
US (1) US12117361B2 (zh)
CN (1) CN113176075B (zh)
WO (1) WO2022227948A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113916351A (zh) * 2021-10-28 2022-01-11 苏州光格科技股份有限公司 光纤振动监测系统
WO2022227948A1 (zh) * 2021-04-26 2022-11-03 南京大学 一种损耗与振动同步监测的相干及偏振衰落抑制方法

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116089782B (zh) * 2022-11-23 2023-11-21 国网甘肃省电力公司临夏供电公司 一种多导向匹配的配电网精细化线损数据系统
CN116907627B (zh) * 2023-09-13 2023-12-19 之江实验室 基于光程差辅助的大动态范围分布式相位传感方法和装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111609919A (zh) * 2020-06-09 2020-09-01 重庆大学 光纤分布式振动和损耗同时检测系统
CN112595492A (zh) * 2020-11-03 2021-04-02 烽火通信科技股份有限公司 用于损耗与扰动事件同步监测的cotdr方法及系统

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103023570B (zh) * 2011-09-26 2016-03-30 富士通株式会社 非线性补偿装置、方法和发射机
CN103259597B (zh) * 2012-02-20 2016-12-07 富士通株式会社 非线性补偿装置、方法和发射机
US9435712B2 (en) * 2014-12-29 2016-09-06 Verizon Patent And Licensing Inc. OTDR signal injection for multi-band optical power balancing
CN108023643B (zh) * 2016-10-31 2020-08-11 富士通株式会社 偏振相关损耗的估计装置、方法以及接收机
CN113176075B (zh) * 2021-04-26 2022-03-25 南京大学 一种损耗与振动同步监测的相干及偏振衰落抑制方法
CN117129023A (zh) * 2023-08-29 2023-11-28 电子科技大学 一种相敏光时域反射仪的相干及偏振衰落抑制方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111609919A (zh) * 2020-06-09 2020-09-01 重庆大学 光纤分布式振动和损耗同时检测系统
CN112595492A (zh) * 2020-11-03 2021-04-02 烽火通信科技股份有限公司 用于损耗与扰动事件同步监测的cotdr方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
何成波等: "基于Φ-OTDR的海光缆扰动监测技术原理及研究进展", 《通信技术》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022227948A1 (zh) * 2021-04-26 2022-11-03 南京大学 一种损耗与振动同步监测的相干及偏振衰落抑制方法
US20230273092A1 (en) * 2021-04-26 2023-08-31 Nanjing University Method for suppressing coherent and polarization-induced fading by simultaneous monitoring of loss and vibration
US12117361B2 (en) * 2021-04-26 2024-10-15 Nanjing University Method for suppressing coherent and polarization-induced fading by simultaneous monitoring of loss and vibration
CN113916351A (zh) * 2021-10-28 2022-01-11 苏州光格科技股份有限公司 光纤振动监测系统
CN113916351B (zh) * 2021-10-28 2024-03-12 苏州光格科技股份有限公司 光纤振动监测系统

Also Published As

Publication number Publication date
US20230273092A1 (en) 2023-08-31
WO2022227948A1 (zh) 2022-11-03
CN113176075B (zh) 2022-03-25
US12117361B2 (en) 2024-10-15

Similar Documents

Publication Publication Date Title
CN113176075B (zh) 一种损耗与振动同步监测的相干及偏振衰落抑制方法
JP7220288B2 (ja) 光ファイバセンシングシステム、方法、構造及び用途
US9909903B2 (en) Apparatus for deriving measurement signals from fibre optic distributed sensors
JP6698164B2 (ja) 周波数合成に基づいた光周波数領域反射方法及びシステム
CN108507662B (zh) 基于多波长双光脉冲的光纤分布式传感方法与装置
CN110470327A (zh) 一种光时域分析仪和分析方法
CN108981890B (zh) 基于多载波非线性调频的分布式振动传感器
US11159200B2 (en) High speed frequency hopping DAS interrogation using AOM-gated re-circulating loop and frequency-shifted receiver LO
CN113418539B (zh) 具有自由多空间分辨率的自外差式的φ-OTDR系统
CN112595492B (zh) 用于损耗与扰动事件同步监测的cotdr方法及系统
CN116592986A (zh) 一种动态应变范围可调的光纤分布式声波传感装置
JP6969506B2 (ja) 光周波数多重型コヒーレントｏｔｄｒ、試験方法、信号処理装置、及びプログラム
CN112697257A (zh) 无衰落多波长分布式声波传感系统和差分旋转矢量叠加法
CN109724529B (zh) 基于多斜坡辅助的大动态范围布里渊快速测量系统
CN111307188B (zh) 一种基于噪声调制的免扫频botda装置
CN113810098A (zh) 一种基于双边带啁啾脉冲调制的光时域反射计
US20230152151A1 (en) Colorless distributed fiber optic sensing / distributed vibration sensing
CN111609919B (zh) 光纤分布式振动和损耗同时检测系统
JP5941877B2 (ja) 光パルス試験装置及び光パルス試験方法
Hu et al. A Frequency-Scanned Φ-OTDR Implemented by a Frequency Shifting Loop Used for Distributed Strain Sensing
Xu et al. Interference Fading Suppression Solution for Phase-sensitive OTDR with Sub Signal Phase Reconstruction
Lei et al. A high fidelity Ф-OTDR system based on compact and flexible multi-frequency probe pulse modulation
Waagaard et al. Real-time phase-recording DAS in 171 km low-loss fiber
Zhang et al. Optical fiber distributed acoustic sensing based on NUFDM-NLFM
EP4446711A1 (en) Method and apparatus of distributed acoustic sensing

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