CN110031082A - 分布式光纤振动监测系统事件定位方法、装置及应用 - Google Patents
分布式光纤振动监测系统事件定位方法、装置及应用 Download PDFInfo
- Publication number
- CN110031082A CN110031082A CN201910389122.7A CN201910389122A CN110031082A CN 110031082 A CN110031082 A CN 110031082A CN 201910389122 A CN201910389122 A CN 201910389122A CN 110031082 A CN110031082 A CN 110031082A
- Authority
- CN
- China
- Prior art keywords
- event
- sensor fibre
- spatial resolution
- max
- field signal
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H9/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
- G01H9/004—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means using fibre optic sensors
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本申请提供了一种分布式光纤振动监测系统事件定位方法、装置及应用,在对传感光纤中发生事件的位置进行定位时,通过对脉冲光在所述传感光纤中传输时,产生的后向瑞利散射光所对应的各干涉场信号的光功率进行比较,得到光功率最大的干涉场信号;然后,计算该光功率最大的干涉场信号的传感位置;最后,根据该传感位置与事件的位置分布规律,确定出传感光纤中发生事件的位置。由于本申请提供的定位方式通过对各干涉场信号进行筛选,再进行事件位置计算,使得其空间分辨率不再单纯的受光信号脉宽决定的,进而可以对事件的位置定位可以更为准确,提高空间分辨精度。
Description
技术领域
本申请涉及光纤传感技术,尤其设计一种分布式光纤振动监测系统事件定位方法、装置及应用。
背景技术
分布式光纤传感技术是应用光纤纵向特性进行测量的技术,因其具有全尺度连续性、网络智能化、长距离等特点,受到广泛应用与关注。
分布式光纤传感技术按照光纤内部信号性质不同,可以分为基于瑞利散射的分布式光纤传感、基于拉曼散射的分布式光纤传感和基于布里渊散射的分布式光纤传感。其中,瑞利散射属于弹性散射,相对于拉曼和布里渊散射有着更高的能量,更容易被监测到,因此,目前多利用基于瑞利散射的分布式光纤振动监测系统,对在整个光纤长度上沿光纤几何路径分布的外部物理参量进行连续的测量。
在分布式光纤振动监测系统中的一个重要指标是空间分辨率,其中,空间分辨率是指当系统同时出现两个事件,能够区分辨别出两个事件之间的最近距离,决定着系统事件漏报率的高低。光纤分布式振动监测系统的空间分辨率主要系统的光脉冲调制器的驱动信号脉宽决定,光脉冲调制器由声光调制器或者电光调制两类。其中,电光调制器(Electrooptic Modulator,EOM)的物理基础是电光效应,但是其消光比较低;声光调制器(Acousto-Optic Modulator,AOM)的物理效应为晶体的声光效应,其主要优点是消光比较大。在实际应用过程中,由于需要监测的距离少则几公里,多则上百公里,为了实现长距离高灵敏度的监测,需要光脉冲调制器具有高消光比,因此,在光纤分布式振动监测系统多数采用声光调制器。
但是,声光调制器的驱动信号的上升沿和下降沿时间较长,所以脉冲宽度一般较大,数值一般在几百纳秒,进而导致光纤分布式振动监测系统的空间分辨率交底,通常在几十米范围。因此,亟待提供一种高空间分辨率的光纤分布式振动监测系统,以满足需求在米级别量级的空间分辨率的工程需要。
发明内容
本申请提供了一种分布式光纤振动监测系统事件定位方法、装置及应用,以解决现有的光纤分布式振动监测系统空间分辨率低的问题。
根据本申请实施例的第一方面,提供了一种分布式光纤振动事件定位方法,应用于基于瑞利散射的分布式光纤振动监测系统,所述方法包括:
将脉冲光注入到所述分布式光纤振动监测系统中的传感光纤中;
获取所述脉冲光在所述传感光纤中传输时,产生的后向瑞利散射光所对应的各干涉场信号;
对各所述干涉场信号的光功率进行比较,得到光功率最大的干涉场信号;
根据接收到所述光功率最大的干涉场信号的时间点、以及所述脉冲光注入传感光纤中的时间点,计算出产生所述光功率最大的干涉场信号的传感位置;
根据所述传感位置,确定所述传感光纤中发生事件的位置。
可选地,根据所述传感位置,确定所述传感光纤中发生事件的位置,包括:
根据所述脉冲光的脉冲宽度,计算出所述分布式光纤振动监测系统的空间分辨率;
获取所述空间分辨率的奇偶性;
当所述空间分辨率为奇数时,则所述传感光纤中发生事件的位置为Lth=Lmax-(S-1)/2;
当所述空间分辨率为偶数时,则所述传感光纤中发生事件的位置为Lth=Lmax-S/2+1;
其中,Lmax为所述传感位置,S为空间分辨率。
可选地,根据所述传感位置,确定所述传感光纤中发生事件的位置,包括:
根据所述脉冲光的脉冲宽度,计算出所述分布式光纤振动监测系统的空间分辨率;
当所述空间分辨率为奇数时,则所述传感光纤中发生事件的位置分别为Lth1=Lmax-(S+1)/2、Lth2=Lmax-(S+1)/2+M-S+1;
当所述空间分辨率为偶数时,则所述传感光纤中发生事件的位置分别为Lth1=Lmax-S/2、Lth2=Lmax-S/2+M-S+1;
其中,Lmax为所述有效传感位置,Lth1为发生事件第一个位置点,Lth2为发生事件第二个位置点,S为空间分辨率M为具有扰动信息的干涉场信号的个数。
可选地,根据所述传感位置,确定所述传感光纤中发生事件的位置,包括:
根据所述脉冲光的脉冲宽度,计算出所述分布式光纤振动监测系统的空间分辨率;
获取所述空间分辨率的奇偶性;
当所述空间分辨率为奇数时,则所述传感光纤中发生事件的位置为Lth=Pmax-(S-1)/2+N;
当所述空间分辨率为偶数时,则所述传感光纤中发生事件的位置为Lth=Pmax-S/2+1+N;
其中,Lmax为所述传感位置,S为空间分辨率,N为所述脉冲光的一个脉冲宽度对应的传播距离。
可选地,根据所述传感位置,确定所述传感光纤中发生事件的位置,包括:
根据所述脉冲光的脉冲宽度,计算出所述分布式光纤振动监测系统的空间分辨率;
当所述空间分辨率为奇数时,则所述传感光纤中发生事件的位置分别为Lth1=Lmax-(S+1)/2+N、Lth2=Lmax-(S+1)/2+M-S+1+N;
当所述空间分辨率为偶数时,则所述传感光纤中发生事件的位置分别为Lth1=Lmax-S/2+N、Lth2=Lmax-S/2+M-S+1+N;
其中,Lmax为所述有效传感位置,Lth1为发生事件第一个位置点,Lth2为发生事件第二个位置点,S为空间分辨率M为具有扰动信息的干涉场信号的个数,N为所述脉冲光的一个脉冲宽度对应的传播距离。
根据本申请实施例的第二方面,提供了一种分布式光纤事件定位装置,应用于基于瑞利散射的分布式光纤振动监测系统,所述装置包括:
干涉场信号获取单元:用于获取脉冲光在所述分布式光纤振动监测系统中的传感光纤中传输时,产生的后向瑞利散射光所对应的各干涉场信号;
光功率比较单元:用于对各所述干涉场信号的光功率进行比较,得到光功率最大的干涉场信号;
传感位置计算单元:用于根据接收到所述光功率最大的干涉场信号的时间点、以及所述脉冲光注入传感光纤中的时间点,计算出产生所述光功率最大的干涉场信号的传感位置;
事件位置确定单元:用于根据所述传感位置,确定所述传感光纤中发生事件的位置。
根据本申请实施例的第三方面,提供了一种分布式光纤振动监测系统,所述系统包括本申请实施例第二方面提供的分布式光纤振动监测系统事件定位装置,还包括激光器、与所述激光器连接的声光调制器、与所述声光调制器连接的环形器、与所述环形器连接的传感光纤、以及与所述分布式光纤振动监测系统事件定位装置和环形器连接的光电探测器。
由上述实施例可见,本申请实施例提供的分布式光纤振动监测系统事件定位方法、装置及应用,在对传感光纤中发生事件的位置进行定位时,通过对脉冲光在所述传感光纤中传输时,产生的后向瑞利散射光所对应的各干涉场信号的光功率进行比较,得到光功率最大的干涉场信号;然后,计算该光功率最大的干涉场信号的传感位置;最后,根据该传感位置与事件的位置分布规律,确定出传感光纤中发生事件的位置。由于本申请提供的定位方式通过对各干涉场信号进行筛选,再进行事件位置计算,使得其空间分辨率不再单纯的受光信号脉宽决定的,进而可以对事件的位置定位可以更为准确,提高空间分辨精度。
应当理解的是,以上的一般描述和后文的细节描述仅是示例性和解释性的,并不能限制本发明。
附图说明
为了更清楚地说明本发明的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,对于本领域普通技术人员而言,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本申请实例提供的一种分布式光纤振动监测系统的基本结构示意图;
图2为本申请实例提供的一种分布式光纤振动监测系统事件定位方法的基本流程示意图;
图3为本申请实例提供的后向瑞利散射离散模型图示意图;
图4为本申请实例提供的后向瑞利散射振动监测光路示意图;
图5为本申请实例提供的在传感光纤Li处发生一个事件时,对应的后向瑞利散射离散模型图示意图;
图6为本申请实例提供的当在空间分辨率范围内发生一个事件时,不同时刻相位解调结果图;
图7为本申请实例提供的当在空间分辨率范围内发生两个事件且空间分辨率为偶数时,不同时刻相位解调结果图;
图8为本申请实例提供的当在空间分辨率范围内发生两个事件且空间分辨率为奇数时,不同时刻相位解调结果图;
图9为本申请实例提供的一种分布式光纤振动监测系统事件定位装置的基本结构示意图。
具体实施方式
这里将详细地对示例性实施例进行说明,其示例表示在附图中。下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本发明相一致的所有实施方式。相反,它们仅是与如所附权利要求书中所详述的、本发明的一些方面相一致的装置和方法的例子。
图1为本实例提供的一种分布式光纤振动监测系统的基本结构示意图。如图1所示,该系统主要包括激光器10、与激光器10连接的声光调制器(AOM,Acousto-opticalModulators)20、与声光调制器20连接的环形器40、为实现脉冲光的功率放大在声光调制器20与环形器40之间本实施例还设置有第一掺铒光纤放大器(EDFA,Erbium Doped FiberAmplifier)30,还包括与环形器40的C2端连接的传感光纤50、与环形器40的C3端连接的光电探测器(PD,Photodetector)80,为实现环形器输出信号的处理,在环形器40的C3端与光电探测器80之间还设置有第二掺铒光纤放大器60和滤波器70,最后,还包括本实施例提供的事件定位装置90。需要说明的是,在实际使用还可以对上述振动监测系统中的部分元件替换、增设一些元件或减少一部分元件。
利用上述振动监测系统,激光器10发出窄线宽连续激光,经过声光调制器20的斩波为光脉冲序列,脉冲光经过第一掺铒光纤放大器30的功率放大,进入环形器40的C1端,然后经过环形器40的C2端注入到传感光纤50,传感光纤50中的后向瑞利散射信号经过环形器40的C3端进入到第二掺铒光纤放大器60,放大后的瑞利散射信号经过滤波器70后,进入到光电探测器80,光电探测器80进行光电转化后的电信号进入事件定位装置90,事件定位装置90通过对接收的电信号进行分析,确定传感光纤50中发生事件的位置,具体的,下面将对事件定位装置90的数据处理过程进行详细介绍。
图2为本申请实例提供的一种分布式光纤振动监测系统事件定位方法的基本流程示意图。如图2所示,该方法具体包括如下步骤:
S110:获取脉冲光在所述分布式光纤振动监测系统中的传感光纤中传输时,产生的后向瑞利散射光所对应的各干涉场信号。
图3为本申请实例提供的后向瑞利散射离散模型图示意图。如图3所示,当在传感光纤中注入相干长度长的窄线宽激光时,本实施例将传感光纤中的散射点看作是一系列离散的反射镜,某个反射镜反射信号可认为是在单位散射长度ΔL范围内,随机分布的散射点的后向散射光的矢量和。
其中,单元散射长度ΔL定义为:
在公式(1)中,Sa为振动监测系统的采样率,C为真空中光速,nf为传感光纤的折射率。
如图3所示,设ΔL内有M个随机分布的瑞利散射点,偏振态相同,第p个反射镜处的光场是M个散射点场矢量和,可以表达为:
在公式(2)中,rp为第p段光纤的M个后向散射点的光场振幅矢量和,定义为第p个反射镜的反射率;φp为M个后向散射点的相位矢量和,定义为第p个反射镜的相位;am为ΔL光纤长度内第m个后向散射点的光场振幅值,Ωm为第m个后向散射点的光场相位值。
由于脉冲宽度为W的激光脉冲注入传感光纤,在不同时刻只有一段光纤中有光,即只有这段光纤被“点亮”,对应的,探测器观察到的光纤点亮长度为(q-1)ΔL,其中:
在公式(3)中,q为点亮光纤内的等效反射镜个数。
基于上述设定,距离初始端iΔL处的后向瑞利散射光的干涉场是i-q+1到i个等效反射镜的场矢量和,Li处的光强表达式如下所示:
在公式(4)中,Pk为第k个等效反射镜的光场偏振态,rk为第k个等效反射镜的光场振幅矢量,为第k个等效反射镜的相位矢量,α为光纤衰减系数,Li是第i个单位散射长度的位置,即:Li=iΔL。
由图3及公式(4)所知,Li处的后向瑞利散射光的干涉场是i-q+1到i个等效反射镜的场矢量和,即脉宽内q个等效反射镜的场矢量和。
图4为本申请实例提供的后向瑞利散射振动监测光路示意图。如图4所示,将一束光频率为f、脉冲宽度为W的高相干脉冲光在t=0时刻从环形器处注入光纤,光电探测器在t时刻得到的光场表达式如下:
在公式(5)中,ak为光场振幅,当时,矩形函数其他情况rect[(t-τk)/W]=0。τk为光纤任意第k个等效反射镜的时间延迟,τk与Lk的关系为τk=2nfLk/C=2nfkΔL/C,N为等效反射镜总个数。
图5为本申请实例提供的在传感光纤Li处发生一个事件时,对应的后向瑞利散射离散模型图示意图。如图5所示,假设在第i个等效反射镜处(即传感光纤Li处)发生一个事件、即施加一个扰动信息ΔΦ,施加在某一个等效反射镜上的扰动信号称为点扰动信号,光电探测器观察到的光脉冲前沿到达传感光纤的第j个等效反射镜的观测时间tj为:
将公式(6)和公式(3)代入公式(5)中,便可以得出,光电探测器在tj时刻得到的光场表达式如下:
在公式(7)中,ak是在tj时刻,光电探测器观测到的脉冲范围内的第k个等效反射镜的光场强度。
当光脉冲前沿到达传感光纤的第i个等效反射镜时,即j=i时,由于在传感光纤Li处发生一个事件,根据公式(7),可得光电探测器在tj时刻得到的光场表达式如下:
根据公式(8),探测到的光功率值Ib(ti)为:
根据公式(9),可得在ti时刻观测到第i个等效反射镜处有扰动信息ΔΦ。
当光脉冲前沿到达传感光纤的第i+1个等效反射镜时,即j=i+1,可得光电探测器在tj+1时刻得到的光场表达式如下:
对应的,探测到的光功率信号Ib(ti+1)为:
根据公式(11),可得在tj+1时,观测到第i+1个等效反射镜处有扰动信息ΔΦ。
以此类推,当光脉冲前沿到达传感光纤的第i+q-2个等效反射镜时,即j=i+q-2,可得光电探测器在tj+q-2时刻得到的光场表达式如下:
对应的,探测的光功率信号Ib(ti+q-2)为:
根据公式(13),可以得出在ti+q-2时刻观测到第i+q-2个等效反射镜处有扰动信息ΔΦ。
当光脉冲前沿继续前进,到达传感光纤的第i+q-1个等效反射镜时,即j=i+q-1,可得光电探测器在tj+q-1时刻得到的光场表达式如下:
对应的,探测的光功率信号Ib(ti+q-1)为:
根据公式(15),可得探测器在ti+q-1时刻观测到第i+q-1个等效反射镜处没有扰动信息ΔΦ。
综上所述,{Ib(ti),Ib(ti+1),……,Ib(ti+q-2)}中均含有扰动信息ΔΦ,{Ib(ti+q-1)……}中不含有有扰动信息ΔΦ,即第i个等效反射镜处的点扰动信息ΔΦ展宽到i个等效反射镜之后的q-2个等效反射镜区域内。
S120:对各所述干涉场信号的光功率进行比较,得到光功率最大的干涉场信号。
基于步骤S110中,探测器在不同时刻观测到各等效反射镜处扰动信息的分布规律,可以将具有ΔΦ的记为有效的干涉场信号,同时,基于光功率信号中含有ΔΦ的余弦项的个数不同,则对应的光功率值也不同,其中,含有ΔΦ的余弦项的个数越多,光功率值越大。
进而,通过对接收到的干涉场信号的光功率进行比较,便可以光功率最大的干涉场信号。
S130:根据接收到所述光功率最大的干涉场信号的时间点、以及所述脉冲光注入传感光纤中的时间点,计算出产生所述光功率最大的干涉场信号的传感位置。
设接收到光功率最大的干涉场信号的时间点为tmax,脉冲光注入传感光纤中的时间点为to,对应的,产生所述光功率最大的干涉场信号的传感位置Lmax则为:
S140:根据所述传感位置,确定所述传感光纤中发生事件的位置
根据步骤S110中,在传感光纤Li处发生事件时,第i个等效反射镜处的点扰动信息ΔΦ展宽到i个等效反射镜之后的q-2个等效反射镜区域内的分布规律,便可以根据传感位置,所述传感光纤中发生事件的位置,其中:
1)当空间分辨率S内发生一个事件时,对应的,传感光纤中发生事件的位置与光功率最大的干涉场信号的传感位置之间的对应关系为:
a、当所述空间分辨率S为奇数时,则所述传感光纤中发生事件的位置为:
Lth=Lmax-(S-1)/2 公式(17)
其中,空间分辨率P为脉冲光信号的脉宽。
b、当所述空间分辨率S为偶数时,则所述传感光纤中发生事件的位置为:
Lth=Lmax-S/2 公式(18)
2)当空间分辨率S内发生两个事件时,对应的,传感光纤中发生事件的位置与光功率最大的干涉场信号的传感位置之间的对应关系为:
a、当所述空间分辨率为奇数时,则所述传感光纤中发生事件的位置分别为
Lth1=Lmax-(S+1)/2、Lth2=Lmax-(S+1)/2+M-S+1 公式(19)
b、当所述空间分辨率为偶数时,则所述传感光纤中发生事件的位置分别为
Lth1=Lmax-S/2、Lth2=Lmax-S/2+M-S+1 公式(20)
在公式(20)中,Lth1为发生事件第一个位置点,Lth2为发生事件第二个位置点,S为空间分辨率,M为具有扰动信息的干涉场信号的个数。
需要说明的是,上述公式(17)至(20)的推导,是以当得到两个或两个光功率最大的干涉场信号时,将得到的第一个光功率最大的干涉场信号的传感位置作为有效传感位置,推导出的,若是以其他功率最大的干涉场信号的传感位置作为有效传感位置,则在上述公式的基础上进行相应的变形即可。
下面以Sa的数值为100Mbps,C=3×108m/s、nf≈1.5,ΔL=1m为例,对上述理论推导进行验证。
图6为本申请实例提供的当在空间分辨率范围内发生一个事件时,不同时刻相位解调结果图。如图6所示,空间分辨率S=5(在实际情况中S的大小为几十米级别),假设事件ΔΦ发生在在传感光纤位置i+5处,则第i+5个瑞利散射点在t=ti+5时刻,探测的光功率信号Ii+5为:
Ii+5=[Ei(ti+5)+Ei+1(ti+5)+Ei+2(ti+5)+Ei+3(ti+5)+Ei+4(ti+5)+Ei+5(ti+5)]*
×[Ei(ti+5)+Ei+1(ti+5)+Ei+2(ti+5)+Ei+3(ti+5)+Ei+4(ti+5)+Ei+5(ti+5)]
对应的,Ii+5中含有10项ΔΦ的余弦项。
同理求得ti+6时刻的光场表达式Ii+6为:
Ii+6中含有16项ΔΦ的余弦项。
依次递推,可以得出如图6所示的关系
ti+7时刻Ii+7中含有18项ΔΦ的余弦项。
ti +8时刻Ii+8中含有16项ΔΦ的余弦项。
ti+9时刻Ii+9中含有10项ΔΦ的余弦项。
ti+10时刻Ii+10中没有含有ΔΦ的余弦项。
进而,Lmax为i+7处,根据公式(17),可计算出发生事件的位置为i+5。
图7为本申请实例提供的当在空间分辨率范围内发生两个事件且空间分辨率为偶数时,不同时刻相位解调结果图。如图7所示,空间分辨率S=4,假设事件ΔΦ发生在在传感光纤位置i+4和i+6处。
对应的,ti+4时刻Ii+4中含有8项ΔΦ的余弦项。
ti+5时刻Ii+5中含有12项ΔΦ的余弦项。
ti+6时刻Ii+6中含有12项ΔΦ和4项2ΔΦ的余弦项。
ti+7时刻Ii+7中含有12项ΔΦ和4项2ΔΦ的余弦项。
ti+8时刻Ii+8中含有12项ΔΦ的余弦项。
ti+9时刻Ii+9中含有8项ΔΦ的余弦项。
进而,Lmax为i+6处,根据公式(20),可计算出发生事件的位置为i+4和i+6处。
图8为本申请实例提供的当在空间分辨率范围内发生两个事件且空间分辨率为奇数时,不同时刻相位解调结果图。如图8所示,空间分辨率S=5,假设事件ΔΦ发生在在传感光纤位置i+4和i+7处。
对应的,ti+4时刻Ii+4中含有10项ΔΦ的余弦项。
ti+5时刻Ii+5中含有16项ΔΦ的余弦项。
ti+6时刻Ii+6中含有18项ΔΦ的余弦项。
ti+7时刻Ii+7中含有18项ΔΦ和4项2ΔΦ的余弦项。
ti+8时刻Ii+8中含有18项ΔΦ和4项2ΔΦ的的余弦项。
ti +9时刻Ii+9中含有18项ΔΦ的余弦项。
ti+10时刻Ii+10中含有16项ΔΦ的余弦项。
ti+11时刻Ii+11中含有10项ΔΦ的余弦项。
进而,Lmax为i+7处,根据公式(19),可计算出发生事件的位置为i+4和i+7处。
需要说明的是,上述推导均是以光脉冲的前沿到达的时刻进行计算的,若是取光脉冲的后沿,则在上述公式的基础上加上脉冲光的一个脉冲宽度对应的传播距离即可,对应如下:
3)当空间分辨率S内发生一个事件时,对应的,传感光纤中发生事件的位置与光功率最大的干涉场信号的传感位置之间的对应关系为:
a、当所述空间分辨率S为奇数时,则所述传感光纤中发生事件的位置为:
Lth=Pmax-(S-1)/2+N 公式(21)
其中,空间分辨率P为脉冲光信号的脉宽。
b、当所述空间分辨率S为偶数时,则所述传感光纤中发生事件的位置为:
Lth=Pmax-S/2+1+N 公式(22)
4)当空间分辨率S内发生两个事件时,对应的,传感光纤中发生事件的位置与光功率最大的干涉场信号的传感位置之间的对应关系为:
a、当所述空间分辨率为奇数时,则所述传感光纤中发生事件的位置分别为
Lth1=Lmax-(S+1)/2+N、Lth2=Lmax-(S+1)/2+M-S+1+N 公式(23)
b、当所述空间分辨率为偶数时,则所述传感光纤中发生事件的位置分别为
Lth1=Lmax-S/2+N、Lth2=Lmax-S/2+M-S+1+N 公式(24)
基于与上述方法相同的发明构思,本申请还提供了一种分布式光纤事件定位装置,应用于基于瑞利散射的分布式光纤振动监测系统。图9为本申请实例提供的一种分布式光纤振动监测系统事件定位装置的基本结构示意图。如图9所示,该装置包括:
干涉场信号获取单元910:用于获取脉冲光在所述分布式光纤振动监测系统中的传感光纤中传输时,产生的后向瑞利散射光所对应的各干涉场信号;
光功率比较单元920:用于对各所述干涉场信号的光功率进行比较,得到光功率最大的干涉场信号;
传感位置计算单元930:用于根据接收到所述光功率最大的干涉场信号的时间点、以及所述脉冲光注入传感光纤中的时间点,计算出产生所述光功率最大的干涉场信号的传感位置;
事件位置确定单元940:用于根据所述传感位置,确定所述传感光纤中发生事件的位置。
本实施例提供的分布式光纤事件定位方法,在对传感光纤中发生事件的位置进行定位时,通过对脉冲光在所述传感光纤中传输时,产生的后向瑞利散射光所对应的各干涉场信号的光功率进行比较,得到光功率最大的干涉场信号;然后,计算该光功率最大的干涉场信号的传感位置;最后,根据该传感位置与事件的位置分布规律,确定出传感光纤中发生事件的位置。由于本申请提供的定位方式通过对各干涉场信号进行筛选,再进行事件位置计算,使得其空间分辨率不再单纯的受光信号脉宽决定的,进而可以对事件的位置定位可以更为准确,提高空间分辨精度,其中,可以将空间分辨率由几十米范围降低到几米的精度。
基于上述方法和装置,本实施例还提供了一种分布式光纤振动监测系统,该系统包括上述提供的所述的分布式光纤振动监测系统事件定位装置,另外还包括激光器、声光调制器、环形器、传感光纤、以及光电探测器等元件。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。
本领域技术人员在考虑说明书及实践这里发明的发明后,将容易想到本发明的其它实施方案。本申请旨在涵盖本发明的任何变型、用途或者适应性变化,这些变型、用途或者适应性变化遵循本发明的一般性原理并包括本发明未发明的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的,本发明的真正范围和精神由下面的权利要求指出。
应当理解的是,本发明并不局限于上面已经描述并在附图中示出的精确结构,并且可以在不脱离其范围进行各种修改和改变。本发明的范围仅由所附的权利要求来限制。
Claims (7)
1.一种分布式光纤振动监测系统事件定位方法,应用于基于瑞利散射的分布式光纤振动监测系统,其特征在于,所述方法包括:
获取脉冲光在所述分布式光纤振动监测系统中的传感光纤中传输时,产生的后向瑞利散射光所对应的各干涉场信号;
对各所述干涉场信号的光功率进行比较,得到光功率最大的干涉场信号;
根据接收到所述光功率最大的干涉场信号的时间点、以及所述脉冲光注入传感光纤中的时间点,计算出产生所述光功率最大的干涉场信号的传感位置;
根据所述传感位置,确定所述传感光纤中发生事件的位置。
2.根据权利要求1所述的方法,其特征在于,根据所述传感位置,确定所述传感光纤中发生事件的位置,包括:
根据所述脉冲光的脉冲宽度,计算出所述分布式光纤振动监测系统的空间分辨率;
获取所述空间分辨率的奇偶性;
当所述空间分辨率为奇数时,则所述传感光纤中发生事件的位置为Lth=Lmax-(S-1)/2;
当所述空间分辨率为偶数时,则所述传感光纤中发生事件的位置为Lth=Lmax-S/2+1;
其中,Lmax为所述传感位置,S为空间分辨率。
3.根据权利要求1所述的方法,其特征在于,根据所述传感位置,确定所述传感光纤中发生事件的位置,包括:
根据所述脉冲光的脉冲宽度,计算出所述分布式光纤振动监测系统的空间分辨率;
当所述空间分辨率为奇数时,则所述传感光纤中发生事件的位置分别为Lth1=Lmax-(S+1)/2、Lth2=Lmax-(S+1)/2+M-S+1;
当所述空间分辨率为偶数时,则所述传感光纤中发生事件的位置分别为Lth1=Lmax-S/2、Lth2=Lmax-S/2+M-S+1;
其中,Lmax为所述有效传感位置,Lth1为发生事件第一个位置点,Lth2为发生事件第二个位置点,S为空间分辨率M为具有扰动信息的干涉场信号的个数。
4.根据权利要求1或2所述的方法,其特征在于,根据所述传感位置,确定所述传感光纤中发生事件的位置,包括:
根据所述脉冲光的脉冲宽度,计算出所述分布式光纤振动监测系统的空间分辨率;
获取所述空间分辨率的奇偶性;
当所述空间分辨率为奇数时,则所述传感光纤中发生事件的位置为Lth=Pmax-(S-1)/2+N;
当所述空间分辨率为偶数时,则所述传感光纤中发生事件的位置为Lth=Pmax-S/2+1+N;
其中,Lmax为所述传感位置,S为空间分辨率,N为所述脉冲光的一个脉冲宽度对应的传播距离。
5.根据权利要求1或3所述的方法,其特征在于,根据所述传感位置,确定所述传感光纤中发生事件的位置,包括:
根据所述脉冲光的脉冲宽度,计算出所述分布式光纤振动监测系统的空间分辨率;
当所述空间分辨率为奇数时,则所述传感光纤中发生事件的位置分别为Lth1=Lmax-(S+1)/2+N、Lth2=Lmax-(S+1)/2+M-S+1+N;
当所述空间分辨率为偶数时,则所述传感光纤中发生事件的位置分别为Lth1=Lmax-S/2+N、Lth2=Lmax-S/2+M-S+1+N;
其中,Lmax为所述有效传感位置,Lth1为发生事件第一个位置点,Lth2为发生事件第二个位置点,S为空间分辨率,M为具有扰动信息的干涉场信号的个数,N为所述脉冲光的一个脉冲宽度对应的传播距离。
6.一种分布式光纤振动监测系统事件定位装置,应用于基于瑞利散射的分布式光纤振动监测系统,其特征在于,所述装置包括:
干涉场信号获取单元:用于获取脉冲光在所述分布式光纤振动监测系统中的传感光纤中传输时,产生的后向瑞利散射光所对应的各干涉场信号;
光功率比较单元:用于对各所述干涉场信号的光功率进行比较,得到光功率最大的干涉场信号;
传感位置计算单元:用于根据接收到所述光功率最大的干涉场信号的时间点、以及所述脉冲光注入传感光纤中的时间点,计算出产生所述光功率最大的干涉场信号的传感位置;
事件位置确定单元:用于根据所述传感位置,确定所述传感光纤中发生事件的位置。
7.一种分布式光纤振动监测系统,其特征在于,所述系统包括权利要求6所述的分布式光纤振动监测系统事件定位装置,还包括激光器、与所述激光器连接的声光调制器、与所述声光调制器连接的环形器、与所述环形器连接的传感光纤、以及与所述分布式光纤振动监测系统事件定位装置和环形器连接的光电探测器。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910389122.7A CN110031082B (zh) | 2019-05-10 | 2019-05-10 | 分布式光纤振动监测系统事件定位方法、装置及应用 |
AU2019204200A AU2019204200B1 (en) | 2019-05-10 | 2019-06-14 | Event positioning method, device and application in distributed fiber vibration monitoring system |
US16/443,397 US11054302B2 (en) | 2019-05-10 | 2019-06-17 | Event positioning method, device and application in distributed fiber vibration monitoring system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910389122.7A CN110031082B (zh) | 2019-05-10 | 2019-05-10 | 分布式光纤振动监测系统事件定位方法、装置及应用 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110031082A true CN110031082A (zh) | 2019-07-19 |
CN110031082B CN110031082B (zh) | 2021-06-18 |
Family
ID=67241780
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910389122.7A Active CN110031082B (zh) | 2019-05-10 | 2019-05-10 | 分布式光纤振动监测系统事件定位方法、装置及应用 |
Country Status (3)
Country | Link |
---|---|
US (1) | US11054302B2 (zh) |
CN (1) | CN110031082B (zh) |
AU (1) | AU2019204200B1 (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110470376A (zh) * | 2019-08-16 | 2019-11-19 | 天津大学 | 一种干涉分布式光纤声传感装置及其传感方法 |
CN110686166A (zh) * | 2019-10-21 | 2020-01-14 | 上海大学 | Sagnac分布光纤传感系统的离散化定位方法 |
CN111147133A (zh) * | 2019-12-24 | 2020-05-12 | 武汉理工光科股份有限公司 | 一种基于φ-OTDR的车流量实时监测系统及方法 |
CN112484837A (zh) * | 2020-11-24 | 2021-03-12 | 电子科技大学 | 光纤空间定位系统及其实现方法 |
CN114172567A (zh) * | 2020-09-10 | 2022-03-11 | 中国石油天然气股份有限公司 | 光纤中故障点的定位方法及装置 |
CN115371794A (zh) * | 2022-10-25 | 2022-11-22 | 杭州水务数智科技股份有限公司 | 基于分布式光纤传感器的地下管廊地表施工事件定位方法 |
CN117109721A (zh) * | 2023-10-25 | 2023-11-24 | 之江实验室 | Das系统中散射点位置查找装置、查找方法及解调方法 |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113310663A (zh) * | 2021-06-08 | 2021-08-27 | 国网山东省电力公司济南供电公司 | 一种分布式光纤振动监测系统 |
CN113588259B (zh) * | 2021-08-03 | 2024-05-31 | 山东中科普锐检测技术有限公司 | 一种设备振动信号标度曲线转折点检测方法及工况监测装置 |
CN114459593B (zh) * | 2022-01-25 | 2024-01-30 | 北京信维科技股份有限公司 | 一种提高光纤振动系统探测距离的方法 |
CN114964459B (zh) * | 2022-04-14 | 2023-08-04 | 同济大学 | 一种载具轮迹实时定位测速系统及方法 |
CN117029999B (zh) * | 2023-10-09 | 2024-01-30 | 山东省科学院激光研究所 | 基于脉冲调制技术的分布式声波传感系统及测量方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4743113A (en) * | 1985-08-29 | 1988-05-10 | Western Atlas International, Inc. | Optical fiber interferometer network |
CN1414283A (zh) * | 2002-12-10 | 2003-04-30 | 西安交通大学 | 基于分布式光纤传感器的油气管线泄漏智能在线监测方法 |
CN201083536Y (zh) * | 2007-04-11 | 2008-07-09 | 北京京安能科技有限公司 | 可定位的光缆振动探测装置 |
CN102628698A (zh) * | 2012-04-06 | 2012-08-08 | 中国科学院上海光学精密机械研究所 | 分布式光纤传感器及信息解调方法 |
CN105157874A (zh) * | 2015-05-20 | 2015-12-16 | 东华大学 | 一种分布式边界安防监测系统及方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6545760B1 (en) * | 1999-03-25 | 2003-04-08 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Apparatus and method for measuring strain in optical fibers using rayleigh scatter |
GB0820658D0 (en) * | 2008-11-12 | 2008-12-17 | Rogers Alan J | Directionality for distributed event location (del) |
CN101825498B (zh) * | 2010-04-13 | 2012-03-21 | 中国计量学院 | 色散与损耗光谱自校正分布式光纤拉曼温度传感器 |
WO2014012411A1 (zh) * | 2012-07-19 | 2014-01-23 | 南京大学 | 基于脉冲编码和相干探测的botda系统 |
US10638514B2 (en) * | 2016-02-04 | 2020-04-28 | Nec Corporation | Method for operating a plurality of wireless networks |
CN106404154B (zh) * | 2016-11-23 | 2022-09-20 | 山东省科学院激光研究所 | 光纤声波探测系统 |
CN108489598A (zh) * | 2018-05-21 | 2018-09-04 | 浙江杰昆科技有限公司 | 相位敏感光时域反射计相位模拟相干解调系统和方法 |
US11391622B2 (en) * | 2019-02-06 | 2022-07-19 | Nec Corporation | Optical fiber sensing system with reduced spatial resolution and noise for multiple application uses |
-
2019
- 2019-05-10 CN CN201910389122.7A patent/CN110031082B/zh active Active
- 2019-06-14 AU AU2019204200A patent/AU2019204200B1/en active Active
- 2019-06-17 US US16/443,397 patent/US11054302B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4743113A (en) * | 1985-08-29 | 1988-05-10 | Western Atlas International, Inc. | Optical fiber interferometer network |
CN1414283A (zh) * | 2002-12-10 | 2003-04-30 | 西安交通大学 | 基于分布式光纤传感器的油气管线泄漏智能在线监测方法 |
CN201083536Y (zh) * | 2007-04-11 | 2008-07-09 | 北京京安能科技有限公司 | 可定位的光缆振动探测装置 |
CN102628698A (zh) * | 2012-04-06 | 2012-08-08 | 中国科学院上海光学精密机械研究所 | 分布式光纤传感器及信息解调方法 |
CN105157874A (zh) * | 2015-05-20 | 2015-12-16 | 东华大学 | 一种分布式边界安防监测系统及方法 |
Non-Patent Citations (1)
Title |
---|
尚盈 等: "采用后向瑞利散射空间差分的周界安防分布式振动监测", 《红外与激光工程》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110470376A (zh) * | 2019-08-16 | 2019-11-19 | 天津大学 | 一种干涉分布式光纤声传感装置及其传感方法 |
CN110686166A (zh) * | 2019-10-21 | 2020-01-14 | 上海大学 | Sagnac分布光纤传感系统的离散化定位方法 |
CN110686166B (zh) * | 2019-10-21 | 2021-11-05 | 上海大学 | Sagnac分布光纤传感系统的离散化定位方法 |
CN111147133A (zh) * | 2019-12-24 | 2020-05-12 | 武汉理工光科股份有限公司 | 一种基于φ-OTDR的车流量实时监测系统及方法 |
CN111147133B (zh) * | 2019-12-24 | 2021-09-14 | 武汉理工光科股份有限公司 | 一种基于φ-OTDR的车流量实时监测系统及方法 |
CN114172567A (zh) * | 2020-09-10 | 2022-03-11 | 中国石油天然气股份有限公司 | 光纤中故障点的定位方法及装置 |
CN114172567B (zh) * | 2020-09-10 | 2023-04-25 | 中国石油天然气股份有限公司 | 光纤中故障点的定位方法及装置 |
CN112484837A (zh) * | 2020-11-24 | 2021-03-12 | 电子科技大学 | 光纤空间定位系统及其实现方法 |
CN112484837B (zh) * | 2020-11-24 | 2021-12-28 | 电子科技大学 | 光纤空间定位系统及其实现方法 |
CN115371794A (zh) * | 2022-10-25 | 2022-11-22 | 杭州水务数智科技股份有限公司 | 基于分布式光纤传感器的地下管廊地表施工事件定位方法 |
CN117109721A (zh) * | 2023-10-25 | 2023-11-24 | 之江实验室 | Das系统中散射点位置查找装置、查找方法及解调方法 |
CN117109721B (zh) * | 2023-10-25 | 2024-01-30 | 之江实验室 | Das系统中散射点位置查找装置、查找方法及解调方法 |
Also Published As
Publication number | Publication date |
---|---|
US20200355547A1 (en) | 2020-11-12 |
US11054302B2 (en) | 2021-07-06 |
AU2019204200B1 (en) | 2020-10-15 |
CN110031082B (zh) | 2021-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110031082A (zh) | 分布式光纤振动监测系统事件定位方法、装置及应用 | |
CN105547460B (zh) | 融合弱反射光栅的双脉冲相位敏感光时域反射计及其方法 | |
CN104034409B (zh) | 基于脉冲编码外调制的分布式光纤振动传感方法及其系统 | |
CN102879081B (zh) | 一种分布式光纤振动系统中的数据处理方法 | |
CN107505041A (zh) | 一种基于相位敏感光时域反射计的相位解调装置和方法 | |
CN108507662B (zh) | 基于多波长双光脉冲的光纤分布式传感方法与装置 | |
CN105157874B (zh) | 一种分布式边界安防监测系统及方法 | |
CN104568119A (zh) | 一种单光源脉冲编码的光纤振动传感系统及其传感方法 | |
CN101634571B (zh) | 光纤脉栅分布传感装置 | |
CN105181111A (zh) | 一种超弱光纤光栅阵列同φ-otdr相结合的光纤振动传感系统 | |
CN106019313A (zh) | 基于偏振双边缘的单像素探测测风激光雷达 | |
US10190940B2 (en) | Spectrally efficient optical frequency-domain reflectometry using I/Q detection | |
CN105973277A (zh) | 基于单光子探测的分布式光纤传感系统的实现装置及方法 | |
Faustov et al. | Comparison of simulated and experimental results for distributed radiation-induced absorption measurement using OFDR reflectometry | |
CN110426067A (zh) | 一种基于多芯光纤的布里渊时域分析系统和方法 | |
Liu et al. | Phase sensitive distributed vibration sensing based on ultraweak fiber Bragg grating array using double-pulse | |
JP6730295B2 (ja) | 共有リスクリンクグループを検出する方法及び装置 | |
CN212320747U (zh) | 一种线性调频的分布式光纤传感装置 | |
CN110082000A (zh) | 多参量分布式智能光纤传感系统 | |
CN107727122B (zh) | 双端探测的联合拉曼和布里渊散射的分布式光纤传感装置 | |
Costa et al. | Reaching pε/√ Hz sensitivity in a distributed optical fiber strain sensor | |
CN107167225B (zh) | 一种分布式光纤应力及振动的传感系统及其传感方法 | |
CN211452794U (zh) | 定位光缆故障点的装置 | |
CN110455495A (zh) | 一种光纤激光器模式稳定性检测装置及方法 | |
CN110375960A (zh) | 一种基于超连续谱光源otdr的装置及方法 |
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 |