CN102538845A - 多点扰动定位方法、光纤分布式扰动传感器 - Google Patents

Abstract

本发明公开了一种多点扰动定位方法，包括以下步骤：S1：分别对接收到的两干涉信号进行隔直操作以去掉两干涉信号中的直流项和低频干扰项，两干涉信号中均包括扰动信号；S2：对两干涉信号分别进行移相90°的操作，并分别利用移相90°的信号除以原干涉信号，得到相除后的两干涉信号；S3：从步骤S2后的两干涉信号中分别提取三角函数内的相位信息；S4：对提取后的两干涉信号进行高通滤波操作以滤除相位缓变的干扰信号；S5：对高通滤波后的两干涉信号分别进行三角变换、频谱分析并分析得到扰动位置。还公开了一种光纤分布式扰动传感器。本发明有效地消除了偏振衰落的影响，实现了多点扰动定位功能，且传感器的光路结构简单、稳定性好、硬件成本低。

Description

多点扰动定位方法、光纤分布式扰动传感器

技术领域

本发明涉及光纤传感器技术领域，特别涉及一种多点扰动定位方法、光纤分布式扰动传感器。

背景技术

光纤传感器由于其高灵敏度，体积小，重量轻，本质安全，电绝缘性，抗电磁干扰，相对成本低，多功能性，可靠性高，硬件匹配光纤通信接口，易于组网，特别是可以实现分布式测量等优良特性，在工业、民用和军事领域具有广泛的应用。其中，光纤分布式扰动传感器在周界安防，油气管线监测，大型结构监测和通信线路监测等领域具有重要意义。

光纤分布式扰动传感器可以对传感光纤上任意一点处的扰动(时变信号)进行监测，得到扰动信号的时域波形，根据扰动事件性质进行判断，给出报警信息；同时给出扰动事件发生的空间位置信息。

目前，根据不同的工作原理，光纤分布式传感器可以分为干涉仪型、光纤光栅型、光时域反射计型，光频域反射计型以及强度调制型等传感技术。

光纤光栅型分布式传感器采用光纤光栅作为敏感元件，在一定长度的间隔之间铺设光纤光栅，通过复用技术实现准分布式传感，因此，光纤光栅型分布式扰动传感器的空间分辨率具有不连续性，且受到光纤光栅空间分布间隔的限制。同时，光纤光栅的集成基于波长复用，在一根光纤上可以复用的光纤光栅数量受到波长区间的限制，其测量长度的增加需要以增大光纤光栅间隔即降低空间分辨率为代价。除了空间分辨率和测量长度间的矛盾外，光纤光栅型分布式传感器的成本也限制了其作为分布式扰动传感器在大范围环境中的应用。

光时域反射计型分布式传感器可以用来检测外界环境中温度或压力的变化，但其响应时间较长，对于外界扰动的实时定位比较困难，不能应用于对时变信号的分布式传感，因此限制了其作为分布式扰动传感器的应用。

光频域反射计型分布式传感器，基于非线性光学效应，布里渊或拉曼散射，可以对外界的温度和压力进行传感，但其传感信号相对微弱，使得信号的检测和解调相对困难，同时其器件成本也相对较高，限制了其在长距离扰动传感中的应用。

强度调制型传感器基于单模光纤和多模光纤中的模式耦合机理，可以实现对扰动的分布式传感，但其灵敏度和精度较低，为了能够在实际中应用还需进一步解决增敏和提高精度的问题。

综上，在光纤分布式传感器中，干涉仪型分布式传感器具有实现原理简单，灵敏度高，响应速度快，硬件成本低，适于长距离传感等优良特性，已经成为光纤分布式扰动传感器的主要技术方案。

目前，干涉仪型分布式光纤扰动传感器的理论方案主要包括单萨格奈克型、双马赫-泽德型、双萨格奈克型、萨格奈克+迈克尔逊型和萨格奈克+马赫-泽德型，双波长萨格奈克型，双调制频率萨格奈克型等。

基于双马赫-泽德干涉仪的光纤分布式扰动传感器的光路原理图如图1所示。激光器发出的光经过耦合器C₁分成两束光，分别沿顺时针和逆时针方向经过由耦合器C₂、C₃和它们之间的两根敏感光纤构成的马赫-泽德干涉仪，在耦合器C₃和C₂处发生干涉并通过探测器PD₂和探测器PD₁接收干涉信号。以上光路中的光纤均为单模光纤。当扰动作用于传感臂上时，应力会引起光纤长度和传播常数的变化，从而在信号臂和参考臂上产生一个相位差的变化。当光纤长度变化ΔL，传播常数变化Δβ时，相位差

可以表示为：

逆时针和顺时针方向传播的干涉光通过PD₁和PD₂接收到的信号可以分别表示为：

其中，变量t表示时间，t₁，t₂和t₃分别是光沿着光纤L1，L2和L3的传播时间，这里忽略了光源到耦合器C₁和耦合器C₂到PD₁的距离。忽略传感光纤和传导光纤的长度差，光纤L3的长度可以近似为L1和L2的长度之和。I₁和I₂由光源输出的光强决定，K₁和K₂是干涉仪的可见度，

是信号臂和参考臂的臂长差引起的初相差。由式(2)和式(3)可得I₁(t)和I₂(t)的传播时间差为：

τ＝t₂+t₃-t₁        (4)

根据传播时间差τ可以得到扰动点的位置：

L₂＝cτ/2n          (5)

上式中，L₂是扰动点到耦合器C₃的距离，即光纤L2的长度，c是真空中的光速，n是光纤的折射率。计算I₁(t)和I₂(t)的互相关函数，根据相关函数的极大值可以求得τ值。根据τ值可确定两个检测信号之间的时间差并反演出扰动的位置，即实现了定位功能。

现有的双马赫泽德干涉仪型光纤分布式扰动传感器的缺点在于：由于长距离传感降低成本的需要，传感器均采用单模光纤。由于单模光纤本身固有的本征双折射和外界随机因素导致的诱导双折射，单模光纤中传输的光波的偏振态会发生随机变化，从而使得发生干涉时，传感臂和参考臂处于相同振动方向的光矢量(电场矢量)分量的幅值发生随机变化，从而使干涉仪输出信号的幅值发生变化，这将引起两路干涉信号输出波形相关性严重降低，导致基于互相关时延的定位算法计算结果错误，从而引起较大的定位误差。特别地，当两臂光波偏振态正交时，将不能发生干涉，干涉仪输出信号的幅值为0，传感器失效，该问题严重影响了传感器的可靠性。使用保偏光纤替代单模光纤作为敏感元件虽然在短距离可以有效地保持传输光的偏振态，但大大地增加了系统的硬件成本，也限制了该方案的实际应用。

发明内容

(一)要解决的技术问题

本发明要解决的技术问题是：如何有效地消除光纤分布式扰动传感器固有的偏振衰落并实现光纤分布式扰动传感器的多点扰动定位。

(二)技术方案

为解决上述技术问题，本发明提供了一种多点扰动定位方法，包括以下步骤：

S1：分别对接收到的两干涉信号进行隔直操作以去掉两干涉信号中的直流项和低频干扰项，所述两干涉信号中均包括扰动信号；

S2：对所述两干涉信号分别进行移相90°的操作，并分别利用移相90°的信号除以原干涉信号，得到相除后的两干涉信号；

S3：从步骤S2之后的两干涉信号中分别提取出三角函数内的相位信息；

S4：对提取后的两干涉信号进行高通滤波操作以滤除相位缓变的干扰信号；

S5：对高通滤波后的两干涉信号分别进行三角变换、频谱分析并分类讨论以得到扰动位置。

其中，所述步骤S1中，隔直操作后还包括步骤：分别获取两干涉信号的光强与可见度信息并消除干涉信号可见度变化，并对所述两干涉信号经过放大、滤波调理，以抑制噪声和干扰。

其中，所述步骤S2中采用希尔伯特变换将所述两干涉信号移相90°。

其中，所述步骤S3中采用三角函数相位提取算法或PGC调制算法从两干涉信号中分别提取出三角函数内的相位信息。

其中，所述步骤S5中分析得到扰动位置的过程具体如下：

通过经过频谱分析的两干涉信号相位谱的比对寻找多点扰动中不同频率成分分量的频率值，再根据这些频率值对应于经过频谱分析的两干涉信号的幅值谱的商联立方程组，从而反解出多点扰动的位置信息。

本发明还提供了一种光纤分布式扰动传感器，包括：光源、第一分光耦合设备、光接收装置、第一迈克尔逊干涉仪及第二迈克尔逊干涉仪，所述光源通过光纤连接所述第一分光耦合设备，所述第一分光耦合设备分别通过光纤分别连接所述第一迈克尔逊干涉仪和第二迈克尔逊干涉仪，所述第一迈克尔逊干涉仪和第二迈克尔逊干涉仪通过光纤分别连接光接收装置，所述光接收装置用于接收光源发出的光通过第一迈克尔逊干涉仪和第二迈克尔逊干涉仪后返回的光信号。

其中，所述第一迈克尔逊干涉仪包括：光纤延迟环、第一法拉第旋光镜、第二法拉第旋光镜及与所述第一分光耦合设备连接的第二分光耦合设备，所述第二分光耦合设备的一侧通过光纤分别连接第一法拉第旋光镜和第二法拉第旋光镜，另一侧连接所述光接收装置，所述第二分光耦合设备用于使第一法拉第旋光镜反射的光信号和第二法拉第旋光镜反射的光信号干涉，并将干涉后的信号传输给所述光接收装置，光纤延迟环位于所述第一迈克尔逊干涉仪传感臂的光纤上；

所述第二迈克尔逊干涉仪包括：第三法拉第旋光镜、第四法拉第旋光镜及与所述第一分光耦合设备连接的第三分光耦合设备，所述第三分光耦合设备的一侧通过光纤分别连接第三法拉第旋光镜和第四法拉第旋光镜，另一侧连接所述光接收装置，所述第三分光耦合设备用于使第三法拉第旋光镜反射的光信号和第四法拉第旋光镜反射的光信号干涉，并将干涉后的信号传输给所述光接收装置。

其中，所述光接收装置包括：第一光电探测器和第二光电探测器，所述第一光电探测器通过光纤连接所述第二分光耦合设备，所述第二光电探测器通过光纤连接所述第三分光耦合设备。

其中，所述光接收装置包括：第四分光耦合设备、第五分光耦合设备、第一相位调制器、第二相位调制器、第一光电探测器、第二光电探测器、第三光电探测器和第四光电探测器，

所述第二分光耦合设备通过光纤连接所述第四分光耦合设备，所述第四分光耦合设备的一端通过光纤连接所述第一光电探测器，另一端通过光纤连接第一相位调制器，所述第一相位调制器通过光纤连接所述第三光电探测器；

所述第三分光耦合设备通过光纤连接所述第五分光耦合设备，所述第五分光耦合设备的一端通过光纤连接所述第二光电探测器，另一端通过光纤连接第二相位调制器，所述第二相位调制器通过光纤连接所述第四光电探测器。

其中，所述光接收装置包括：第四分光耦合设备、第五分光耦合设备、第一相位调制器、第二相位调制器、第一光电探测器和第二光电探测器，且所述第一光电探测器和第二光电探测器均为双通道光电探测器；

所述第二分光耦合设备通过光纤连接所述第四分光耦合设备，所述第四分光耦合设备的一端通过光纤连接所述第一光电探测器的一个通道，另一端通过光纤连接第一相位调制器，所述第一相位调制器通过光纤连接所述第一光电探测器的另一个通道；

所述第三分光耦合设备通过光纤连接所述第五分光耦合设备，所述第五分光耦合设备的一端通过光纤连接所述第二光电探测器的一个通道，另一端通过光纤连接第二相位调制器，所述第二相位调制器通过光纤连接所述第二光电探测器的另一个通道。

(三)有益效果

本发明具有如下有益效果：

1、基于迈克尔逊干涉仪结构，采用了法拉第旋光镜技术，有效地消除了偏振衰落的影响；

2、采用双迈克尔逊光路结构，具有光路结构简单，硬件成本低等优势；

3、通过对扰动信号的预处理消除了在定位计算过程中可能受到的光功率波动和信号偏振诱导衰落等因素引起的干涉信号可见度的变化，从而也间接消除了可见度变化引起的可能的传感器定位失效的问题；

4、在提取出三角函数内的信号后进行了高通滤波，滤除了相位漂移信号，从而消除了相位缓变的影响；

5、通过移相90°，三角变换，相位提取算法，频域谱分析等手段成功地提取出了扰动的位置信息，实现了双迈克尔逊型光路的多点扰动定位功能。

附图说明

图1是现有的基于双马赫-泽德干涉仪型光纤分布式扰动传感器的光路原理图；

图2是本发明实施例的一种光纤分布式扰动传感器的光路结构图；

图3是本发明实施例的另一种光纤分布式扰动传感器的光路结构图；

图4是本发明实施例的另一种光纤分布式扰动传感器的光路结构图；

图5是基于图2的传感器的多点扰动定位方法流程图。

具体实施方式

下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明，但不用来限制本发明的范围。

实施例1

如图2所示，为本实施例的光纤分布式扰动传感器的光路结构图，包括：光源、光接收装置(第一光电探测器PD₁、第二光电探测器PD₂)、耦合器C₁、第一迈克尔逊干涉仪及第二迈克尔逊干涉仪。两迈克尔逊干涉仪结构相同，第一迈克尔逊干涉仪包括：第一法拉第旋光镜FRM₁、第二法拉第旋光镜FRM₂、光纤延迟环D及耦合器C₂，第二迈克尔逊干涉仪包括：第三法拉第旋光镜FRM₃、第四法拉第旋光镜FRM₄及耦合器C₃，光源采用激光器Laser，Laser通过光纤连接耦合器C₁，耦合器C₁分别通过光纤连接第一迈克尔逊干涉仪和第二迈克尔逊干涉仪，第一迈克尔逊干涉仪和第二迈克尔逊干涉仪通过光纤分别连接第一光电探测器PD₁和第二光电探测器PD₂。具体地，耦合器C₂的一侧通过光纤分别连接第一法拉第旋光镜FRM₁和第二法拉第旋光镜FRM₂，另一侧通过光纤连接耦合器C₁，同时还连接第一光电探测器PD₁，耦合器C₂用于使第一法拉第旋光镜FRM₁反射的光信号和第二法拉第旋光镜FRM₂反射的光信号干涉，并将干涉后的信号传输给所述第一光电探测器PD₁。耦合器C₃的一侧通过光纤分别连接第三法拉第旋光镜FRM₃和第四法拉第旋光镜FRM₄，另一侧通过光纤连接耦合器C₁，同时还连接第二光电探测器PD₂，耦合器C₃用于使第三法拉第旋光镜FRM₃反射的光信号和第四法拉第旋光镜FRM₄反射的光信号干涉，并将干涉后的信号传输给第二光电探测器PD₂。

图2中，激光器Laser发出的光经过耦合器C₁分成两束光，并分别注入两个迈克尔逊干涉仪。迈克尔逊干涉仪中的光波分别沿着参考臂和传感臂传播，在光纤的末端遇到法拉第旋光镜(FRM₁、FRM₂、FRM₃、FRM₄)发生反射，并最终分别在耦合器C₂和C₃处发生干涉。在实际工程应用中，参考臂光纤L₂、L₄盘起，传感臂光纤L₁、L₃布置于同一光缆中。光纤延迟环D位于第一迈克尔逊干涉仪的传感臂上，在传感臂上的位置可以位于第一迈克尔逊干涉仪传感臂的首端(即靠近第二分光耦合设备)，也可以位于传感臂的末端(即靠近第一法拉第旋光镜)，均可以实现定位功能。

实施例2

如图3所示，本实施例中的光纤分布式扰动传感器与实施例2不同的是光接收装置还包括：耦合器C₄、耦合器C₅、第一相位调制器P₁、第二相位调制器P₂、第三光电探测器PD₃和第四光电探测器PD₄。且耦合器C₂不直接连接PD₁，耦合器C₃不直接连接PD₂。

耦合器C₂通过光纤连接耦合器C₄，耦合器C₄的一端通过光纤连接第一光电探测器PD₁，另一端通过光纤连接第一相位调制器P₁，第一相位调制器P₁通过光纤连接所述第三光电探测器PD₃；

耦合器C₃通过光纤连接耦合器C₅，耦合器C₅的一端通过光纤连接第二光电探测器PD₂，另一端通过光纤连接第二相位调制器P₂，第二相位调制器P₂通过光纤连接所述第四光电探测器PD₄。

耦合器C₄将耦合器C₂传送过来的信号分成两个信号，一个通过光纤直接送入第一光电探测器PD₁，另一个通过第一相位调制器P₁进行移相90°后再送入第三光电探测器PD₃。

实施例3

如图4所示，本实施例中的光纤分布式扰动传感器与实施例2的结构基本相同，不同的是第一光电探测器PD₁和第二光电探测器PD₂均为双通道光电探测器，因此不需要第三光电探测器PD₃和第四光电探测器PD₄。将耦合器C₂出来的光纤接入PD₁的一个通道，将P₁出来的光纤分别接入PD₁的另一个通道；将耦合器C₃出来的光纤接入PD₂的一个通道，将P₂出来的光纤分别接入PD₂的另一个通道。

上述实施例1、实施例2和实施例3中的耦合器均为2×2耦合器，实现分光和耦合作用，也可以用其他具有分光和耦合作用的器件来代替，如：环形器、分束器、半反半透膜等。

实施例4

本实施例提供了一种基于上述实施例1、实施例2或实施例3的光纤分布式扰动传感器的多点扰动定位方法。该方案多点定位的原理框图如图5所示，以两点同时扰动为例，当两点扰动f₁(t)和f₂(t)同时发生时，光纤长度和传播常数会发生变化，从而引起干涉仪中的相位变化。根据光纤传感理论，输出相位信息的变化正比于扰动信号，有：

其中，t表示时间，B是与扰动相位对应的比例因子，和

分别是两个扰动引起的相位差变化。

设激光光源发出的光为：E_in，先考虑受振后往耦合器方向传播的光。从耦合器C₂输出到PD₁的两束干涉光为：

$E_{11}=\frac{1}{4}{E}_{\mathrm{in}}\exp \left\{j\right[\frac{2\mathrm{\πn}\left(2L\right)}{\mathrm{\λ}}-\frac{\mathrm{\π}}{2}\·2\left]\right\}---\left(8\right)$

$E_{12}=\frac{1}{4}{E}_{\mathrm{in}}\exp \left\{j\right[\frac{2\mathrm{\πn}\left(2(L+\mathrm{\Δ}{L}_1+D)\right)}{\mathrm{\λ}}+{\mathrm{Bf}}_1(t-\frac{{\mathrm{nz}}_1}{c})+{\mathrm{Bf}}_2(t-\frac{{\mathrm{nz}}_2}{c})-\frac{\mathrm{\π}}{2}\·2\left]\right\}---\left(9\right)$

从耦合器C₃输出到PD₂的两束干涉光为：

$E_{21}=\frac{1}{4}{E}_{\mathrm{in}}\exp \left\{j\right[\frac{2\mathrm{\πn}\left(2L\right)}{\mathrm{\λ}}-\frac{\mathrm{\π}}{2}\·2\left]\right\}---\left(10\right)$

$E_{22}=\frac{1}{4}{E}_{\mathrm{in}}\exp \left\{j\right[\frac{2\mathrm{\πn}\left(2(L+\mathrm{\Δ}{L}_2)\right)}{\mathrm{\λ}}+{\mathrm{Bf}}_1(t-\frac{{\mathrm{nz}}_1}{c})+{\mathrm{Bf}}_2(t-\frac{{\mathrm{nz}}_2}{c})-\frac{\mathrm{\π}}{2}\·2\left]\right\}---\left(11\right)$

其中，λ为两束干涉光的波长，L为参考臂光纤L₂、L₄长度，ΔL₁、ΔL₂为传感光纤L₁、L₃与参考臂光纤L₂、L₄的长度差，D为延迟光纤长度；忽略耦合器C₂和C₃之间的距离，z₁和z₁分别为两个扰动作用点到耦合器C₂的距离，c是真空中的光速，n为单模光纤的折射率，Bf₁(t)和Bf₂(t)分别为两个扰动引入的相位调制，(π/2)·m为光波m次(式(8)～(11)中为2次)经过2×2耦合器的跨接耦合所带来的相移，j是虚数单位。

根据干涉理论，PD₁和PD₂接收到的干涉光强分别为：

$I_{11}=\frac{1}{8}{I}_0+\frac{1}{8}{I}_0{K}_1\cos \{\frac{4\mathrm{\πn\Δ}{L}_1}{\mathrm{\λ}}+\frac{4\mathrm{\πnD}}{\mathrm{\λ}}+{\mathrm{Bf}}_1(t-\frac{{\mathrm{nz}}_1}{c})+{\mathrm{Bf}}_2(t-\frac{{\mathrm{nz}}_2}{c})\}---\left(12\right)$

$I_{21}=\frac{1}{8}{I}_0+\frac{1}{8}{I}_0{K}_2\cos \{\frac{4\mathrm{\πn\Δ}{L}_2}{\mathrm{\λ}}+{\mathrm{Bf}}_1(t-\frac{{\mathrm{nz}}_1}{c})+{\mathrm{Bf}}_2(t-\frac{{\mathrm{nz}}_2}{c})\}---\left(13\right)$

其中，I₀取决于激光器的输出功率，而K₁和K₂分别表示两个迈克尔逊干涉仪的可见度。

再考虑受振时往法拉第旋光镜FRM方向传播的光，同理可得PD₁和PD₂接收到的干涉光强为：

$I_{12}=\frac{1}{8}{I}_0+\frac{1}{8}{I}_0{K}_1\cos$ (14)

$\{\frac{4\mathrm{\πn\Δ}{L}_1}{\mathrm{\λ}}+\frac{4\mathrm{\πnD}}{\mathrm{\λ}}+{\mathrm{Bf}}_1(t-n\frac{2(L+\mathrm{\Δ}{L}_1+D)-z_1}{c})+{\mathrm{Bf}}_2(t-n\frac{2(L+\mathrm{\Δ}{L}_1+D)-z_2}{c})\}$

$I_{22}=\frac{1}{8}{I}_0+\frac{1}{8}{I}_0{K}_2\cos$ (15)

$\{\frac{4\mathrm{\πn\Δ}{L}_2}{\mathrm{\λ}}+{\mathrm{Bf}}_1(t-n\frac{2(L+\mathrm{\Δ}{L}_2)-z_1}{c})+{\mathrm{Bf}}_2(t-n\frac{2(L+\mathrm{\Δ}{L}_2)-z_2}{c})\}$

干涉光强叠加可得PD₁和PD₂接收到的实际干涉光强为：

$I_1\left(t\right)=\frac{1}{4}{I}_0+\frac{1}{8}{I}_0{K}_1\left\{\cos \right[\frac{4\mathrm{\πn\Δ}{L}_1}{\mathrm{\λ}}+\frac{4\mathrm{\πnD}}{\mathrm{\λ}}+{\mathrm{Bf}}_1(t-\frac{{\mathrm{nz}}_1}{c})+{\mathrm{Bf}}_2(t-\frac{{\mathrm{nz}}_2}{c})]$ (16)

$+\cos [\frac{4\mathrm{\πn\Δ}{L}_1}{\mathrm{\λ}}++\frac{2\mathrm{\πnD}}{\mathrm{\λ}}+{\mathrm{Bf}}_1(t-\frac{2(L+\mathrm{\Δ}{L}_1+D)-z_1}{c}n)+{\mathrm{Bf}}_2(t-\frac{2(L+\mathrm{\Δ}{L}_1+D)-z_2}{c}n)]\}$

$I_2\left(t\right)=\frac{1}{4}{I}_0+\frac{1}{8}{I}_0{K}_2\left\{\cos \right[\frac{4\mathrm{\πn\Δ}{L}_2}{\mathrm{\λ}}+{\mathrm{Bf}}_1(t-\frac{{\mathrm{nz}}_1}{c})+{\mathrm{Bf}}_2(t-\frac{{\mathrm{nz}}_2}{c})]$ (17)

$+\cos [\frac{4\mathrm{\πn\Δ}{L}_2}{\mathrm{\λ}}+{\mathrm{Bf}}_1(t-\frac{n(2L+\mathrm{\Δ}{L}_2-z_1)}{c})+{\mathrm{Bf}}_2(t-\frac{n(2L+\mathrm{\Δ}{L}_2-z_2)}{c})]\}$

由于传感光纤和参考光纤的长度差相对于臂长差来说很小，据此化简式(16)和(17)，可得：

(18)

(19)

其中，和

分别表示两个迈克尔逊干涉仪臂长差引起的初相差。

下面对探测器接收到的信号进行预处理。首先滤除直流项，再通过分段求取峰峰值的计算方法求取I₀K₁和I₀K₂，来消除因光功率波动和信号偏振诱导衰落引起的干涉信号可见度的变化，可得：

(20)

(21)

本实施例中，也可以对隔直后的信号采取抗偏振衰落技术或光功率稳定控制技术来消除干涉信号可见度的变化，并通过求取峰峰值或者其它方式求取光强信息，从而得到式(20)和(21)。或者：对隔直后的信号采取抗偏振衰落技术或光功率稳定控制技术来消除干涉信号可见度的变化，并保留光强信息，则在计算式(24)和(25)的过程中可以直接约去光强信息。优选地，在隔直后对干涉信号经过放大和滤波等调理，以抑制噪声和干扰。

若采用上述实施例1的光路，则对(20)和(21)式分别做希尔伯特变换得到式(22)和(23)，即对信号I′₁(t)和I′₂(t)移相了90°。若采用上述实施例2或实施例3的光路，则实施例2中的第三光电探测器和第四光电探测器(或者实施例3中第一探测器连接P₁的通道和第二探测器连接P₂的通道)的输出分别为式(22)和(23)：

(22)

(23)

将式(22)除以(20)，式(23)除以(21)，可得：

(24)

(25)

根据万能三角公式对式(24)和(25)进行变换，可得：

(26)

(27)

对式(26)和(27)应用三角函数相位提取算法提取两干涉信号中分别提取出三角函数内的相位信息，可得：

(28)

(29)

本实施例中，也可以不经过万能三角公式对式(24)和(25)变换，直接通过三角函数相位提取算法对式(24)和(25)进行操作以得到式(28)和(29)。另外，也可以采用相位产生载波(Phase generated carrier，PGC)调制算法或其它相关算法从两干涉信号中分别提取出三角函数内的相位信息。

由于

和

是缓变信号，因此采用高通滤波以滤除相位缓变的干扰信号，经过高通滤波后，可得：

${I_3}^{\′\′}\left(t\right)={\mathrm{Bf}}_1(t-\frac{z_1}{c}n)+{\mathrm{Bf}}_2(t-\frac{z_2}{c}n)+{\mathrm{Bf}}_1(t-\frac{2(L+D)-z_1}{c}n)+{\mathrm{Bf}}_2(t-\frac{2(L+D)-z_2}{c}n)$

(30)

${I_4}^{\′\′}\left(t\right)={\mathrm{Bf}}_1(t-\frac{2L-z_1}{c}\·n)+{\mathrm{Bf}}_2(t-\frac{2L-z_2}{c}\·n)+{\mathrm{Bf}}_1(t-\frac{{\mathrm{nz}}_1}{c})+{\mathrm{Bf}}_2(t-\frac{{\mathrm{nz}}_2}{c})---\left(31\right)$

通常，假设两个扰动信号分别为f₁(t)和f₂(t)：

式中，A_i、A_j′，ω_i、ω_j′和

分别表示两个扰动信号的幅值、角频率和初相位，N和M分别表示两个扰动对应频率成分分量的数量。式(32)和(33)分别代入式(30)和(31)，可得：

(34)

(35)

根据和差积化公式，可得：

(36)

(37)

对式(36)和(37)进一步化简，提取两个扰动中相同频率成分分量，可得：

(38)

(39)

式中，A_k′″，A_k″″，ω_k″，

和分别表示两个扰动中相同频率成分分量在两个探测通道处理后的幅值(A_k′″、A_k″″)、角频率(ω_k″)和初相位

其中，P是两个扰动中相同频率成分的数量。

对式(38)和(39)进行频谱分析，经过简单的分类讨论就可以反解出两个扰动的位置信息z₁和z₂。先通过经过频谱分析的两路信号相位谱的比对寻找多点扰动中不同频率成分分量的频率值，再根据这些频率值对应于经过频谱分析的两路信号的幅值谱的商联立方程组，从而反解出多点扰动的位置信息。具体讨论过程如下：

观察式(38)和(39)，可发现：I₅(t)和I₆(t)的频谱对应的横轴频率点均相同，而频率点不外乎分为三种情况，即为ω_i、ω′_j或ω_k″。区分出ω_i、ω′_j与ω_k″的方法就是通过比较I₅(t)和I₆(t)相同频率点的相位谱。可以看出，相位谱相差为nωD/c时(其中ω代表对应的频率点)时，对应的频率点为ω_i或ω′_j。取该频率点进行运算。当ω＝ω_i或ω＝ω_j′时，可用I₅(t)和I₆(t)对应频率点的幅值谱作商，可约去B和A_i或A′_j，得到未知数为z₁和z₂的方程组，则可以解得z₁和z₂的位置信息。

将两点扰动的定位算法推广到三点和三点以上的扰动，可以实现基于双迈克尔逊干涉仪的光纤分布式扰动传感器多点同时扰动的检测和定位。

光路结构中采用的法拉第旋光镜技术可以解决偏振衰落问题，其基本原理如下：

单模光纤的双折射效应一般可用椭圆延迟器来描述，用琼斯矢量矩阵表示为：

$\stackrel{\→}{R}=\frac{\mathrm{\α}}{d}\left(\begin{array}{cc}a& {-b}^{*}\\ b& a^{*}\end{array}\right)---\left(40\right)$

其中，^*表示共轭，α为光纤损耗，d＝(a·a^*+b·b^*)，a、b与单模光纤双折射特性有关。

当光波从同一段单模光纤的另一端注入时，引入的双折射用琼斯矢量矩阵可表示为：

$\stackrel{\←}{R}=\frac{\mathrm{\α}}{d}\left(\begin{array}{cc}a& -b\\ b^{*}& a^{*}\end{array}\right)---\left(41\right)$

旋转角为θ的法拉第旋光镜的琼斯矢量矩阵表示为：

$\left(\begin{array}{cc}\cos \mathrm{\θ}& -\sin \mathrm{\θ}\\ \sin \mathrm{\θ}& \cos \mathrm{\θ}\end{array}\right)---\left(42\right)$

若光波反向经过同一法拉第旋光镜，其琼斯矢量矩阵表示为：

$\left(\begin{array}{cc}\cos \mathrm{\θ}& \sin \mathrm{\θ}\\ -\sin \mathrm{\θ}& \cos \mathrm{\θ}\end{array}\right)---\left(43\right)$

则往返经过此法拉第旋光镜的琼斯矩阵表示为：

$T=u\·\left(\begin{array}{cc}\cos \mathrm{\θ}& -\sin \mathrm{\θ}\\ \sin \mathrm{\θ}& \cos \mathrm{\θ}\end{array}\right)\left(\begin{array}{cc}-1& 0\\ 0& 1\end{array}\right)\left(\begin{array}{cc}\cos \mathrm{\θ}& \sin \mathrm{\θ}\\ -\sin \mathrm{\θ}& \cos \mathrm{\θ}\end{array}\right)---\left(44\right)$

其中，u为标量幅度项，用于体现反射镜的损耗和反射系数。化简后，可得到：

$T=u\·\left(\begin{array}{cc}-\cos 2\mathrm{\θ}& -\sin 2\mathrm{\θ}\\ -\sin 2\mathrm{\θ}& \cos 2\mathrm{\θ}\end{array}\right)---\left(45\right)$

当旋转角为45°时，有：

$T=u\·\left(\begin{array}{cc}0& -1\\ -1& 0\end{array}\right)---\left(46\right)$

当光波从单模光纤的一端注入并且经过传输和反射后返回光波注入端时，单模光纤和法拉第旋光镜的作用采用琼斯矢量矩阵表示为：

$R^{\′}=\stackrel{\←}{R}T\stackrel{\→}{R}={\mathrm{\α}}^{2}u\left[\begin{array}{cc}0& -1\\ -1& 0\end{array}\right]---\left(47\right)$

由于R′是常量矩阵，则经过法拉第旋光镜的处理，返回单模光纤注入端的光波偏振态与单模光纤的双折射效应无关，且注入端的光波偏振态与返回注入端的光波偏振态一致。

所以，将旋转角为45°的法拉第旋光镜置于信号臂和参考臂的末端，当光源发出的光经过耦合器分束后，分别进入信号臂和参考臂，在法拉第旋光镜的作用下，返回耦合器的光波偏振态与分束前的光波偏振态一致，两光束发生完全干涉，偏振衰落得以消除。

以上实施方式仅用于说明本发明，而并非对本发明的限制，有关技术领域的普通技术人员，在不脱离本发明的精神和范围的情况下，还可以做出各种变化和变型，因此所有等同的技术方案也属于本发明的范畴，本发明的专利保护范围应由权利要求限定。

Claims (10)

1.一种多点扰动定位方法，其特征在于，包括以下步骤：

S1：分别对接收到的两干涉信号进行隔直操作以去掉两干涉信号中的直流项和低频干扰项，所述两干涉信号中均包括扰动信号；

S2：对所述两干涉信号分别进行移相90°的操作，并分别利用移相90°的信号除以原干涉信号，得到相除后的两干涉信号；

S3：从步骤S2之后的两干涉信号中分别提取出三角函数内的相位信息；

S4：对提取后的两干涉信号进行高通滤波操作以滤除相位缓变的干扰信号；

S5：对高通滤波后的两干涉信号分别进行三角变换、频谱分析并分类讨论以得到扰动位置。

2.如权利要求1所述的多点扰动定位方法，其特征在于，所述步骤S1中，隔直操作后还包括步骤：分别获取两干涉信号的光强与可见度信息并消除干涉信号可见度变化，并对所述两干涉信号经过放大、滤波调理，以抑制噪声和干扰。

3.如权利要求1所述的多点扰动定位方法，其特征在于，所述步骤S2中采用希尔伯特变换将所述两干涉信号移相90°。

4.如权利要求1所述的多点扰动定位方法，其特征在于，所述步骤S3中采用三角函数相位提取算法或PGC调制算法从两干涉信号中分别提取出三角函数内的相位信息。

5.如权利要求1～4中任一项所述的多点扰动定位方法，其特征在于，所述步骤S5中分类讨论以得到扰动位置的过程具体如下：

通过经过频谱分析的两干涉信号相位谱的比对寻找多点扰动中不同频率成分分量的频率值，再根据这些频率值对应于经过频谱分析的两干涉信号的幅值谱的商联立方程组，从而反解出多点扰动的位置信息。

6.一种光纤分布式扰动传感器，其特征在于，包括：光源、第一分光耦合设备、光接收装置、第一迈克尔逊干涉仪及第二迈克尔逊干涉仪，所述光源通过光纤连接所述第一分光耦合设备，所述第一分光耦合设备分别通过光纤分别连接所述第一迈克尔逊干涉仪和第二迈克尔逊干涉仪，所述第一迈克尔逊干涉仪和第二迈克尔逊干涉仪通过光纤分别连接光接收装置，所述光接收装置用于接收光源发出的光通过第一迈克尔逊干涉仪和第二迈克尔逊干涉仪后返回的光信号。

7.如权利要求6所述的光纤分布式扰动传感器，其特征在于，所述第一迈克尔逊干涉仪包括：光纤延迟环、第一法拉第旋光镜、第二法拉第旋光镜及与所述第一分光耦合设备连接的第二分光耦合设备，所述第二分光耦合设备的一侧通过光纤分别连接第一法拉第旋光镜和第二法拉第旋光镜，另一侧连接所述光接收装置，所述第二分光耦合设备用于使第一法拉第旋光镜反射的光信号和第二法拉第旋光镜反射的光信号干涉，并将干涉后的信号传输给所述光接收装置，光纤延迟环位于所述第一迈克尔逊干涉仪传感臂的光纤上；

所述第二迈克尔逊干涉仪包括：第三法拉第旋光镜、第四法拉第旋光镜及与所述第一分光耦合设备连接的第三分光耦合设备，所述第三分光耦合设备的一侧通过光纤分别连接第三法拉第旋光镜和第四法拉第旋光镜，另一侧连接所述光接收装置，所述第三分光耦合设备用于使第三法拉第旋光镜反射的光信号和第四法拉第旋光镜反射的光信号干涉，并将干涉后的信号传输给所述光接收装置。

8.如权利要求7所述的光纤分布式扰动传感器，其特征在于，所述光接收装置包括：第一光电探测器和第二光电探测器，所述第一光电探测器通过光纤连接所述第二分光耦合设备，所述第二光电探测器通过光纤连接所述第三分光耦合设备。

9.如权利要求7所述的光纤分布式扰动传感器，其特征在于，所述光接收装置包括：第四分光耦合设备、第五分光耦合设备、第一相位调制器、第二相位调制器、第一光电探测器、第二光电探测器、第三光电探测器和第四光电探测器，

所述第二分光耦合设备通过光纤连接所述第四分光耦合设备，所述第四分光耦合设备的一端通过光纤连接所述第一光电探测器，另一端通过光纤连接第一相位调制器，所述第一相位调制器通过光纤连接所述第三光电探测器；

所述第三分光耦合设备通过光纤连接所述第五分光耦合设备，所述第五分光耦合设备的一端通过光纤连接所述第二光电探测器，另一端通过光纤连接第二相位调制器，所述第二相位调制器通过光纤连接所述第四光电探测器。

10.如权利要求7所述的光纤分布式扰动传感器，其特征在于，所述光接收装置包括：第四分光耦合设备、第五分光耦合设备、第一相位调制器、第二相位调制器、第一光电探测器和第二光电探测器，且所述第一光电探测器和第二光电探测器均为双通道光电探测器；

所述第二分光耦合设备通过光纤连接所述第四分光耦合设备，所述第四分光耦合设备的一端通过光纤连接所述第一光电探测器的一个通道，另一端通过光纤连接第一相位调制器，所述第一相位调制器通过光纤连接所述第一光电探测器的另一个通道；

所述第三分光耦合设备通过光纤连接所述第五分光耦合设备，所述第五分光耦合设备的一端通过光纤连接所述第二光电探测器的一个通道，另一端通过光纤连接第二相位调制器，所述第二相位调制器通过光纤连接所述第二光电探测器的另一个通道。

Images