CN113091947B - Ase时域相关压缩的分布式光纤拉曼温度传感装置和方法 - Google Patents

Ase时域相关压缩的分布式光纤拉曼温度传感装置和方法 Download PDF

Info

Publication number
CN113091947B
CN113091947B CN202110420282.0A CN202110420282A CN113091947B CN 113091947 B CN113091947 B CN 113091947B CN 202110420282 A CN202110420282 A CN 202110420282A CN 113091947 B CN113091947 B CN 113091947B
Authority
CN
China
Prior art keywords
ase
optical fiber
raman
time domain
temperature
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
CN202110420282.0A
Other languages
English (en)
Other versions
CN113091947A (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.)
Taiyuan University of Technology
Original Assignee
Taiyuan University of Technology
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 Taiyuan University of Technology filed Critical Taiyuan University of Technology
Priority to CN202110420282.0A priority Critical patent/CN113091947B/zh
Publication of CN113091947A publication Critical patent/CN113091947A/zh
Application granted granted Critical
Publication of CN113091947B publication Critical patent/CN113091947B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K11/00—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00
    • G01K11/32—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using changes in transmittance, scattering or luminescence in optical fibres
    • 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
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring Temperature Or Quantity Of Heat (AREA)
  • Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)

Abstract

本发明涉及分布式光纤传感领域，公开了一种ASE时域相关压缩的分布式光纤拉曼温度传感装置和方法，装置包括ASE信号发生器，其发出ASE激光脉冲信号分为两束，其中一束作为参考光束被第二探测器探测；另一束作为探测光束经波分复用器后入射到传感光纤，产生的后向拉曼散射光经波分复用器、掺铒光纤放大器后，被第一探测器探测；数据采集卡采集两个光电探测器的探测信号后发送至计算机；计算机用于对ASE拉曼后向散射信号进行时域差分信号重构，进行时域相关压缩运算，进而定位和解调光纤沿线的温度突变点信息。本发明可以达到毫米量级的空间分辨率，提高了系统的信噪比，温度精度且无需定标过程。

Description

ASE时域相关压缩的分布式光纤拉曼温度传感装置和方法
技术领域
本发明涉及分布式光纤传感领域，具体是一种基于ASE时域相关压缩的分布式光纤拉曼传感装置和方法，其空间分辨率可达毫米量级。
背景技术
分布式光纤拉曼传感系统可以连续测量传感光纤沿线的分布式温度特征信息。在分布式光纤拉曼传感系统中，传感光纤沿线的环境温度可以对光纤中的拉曼散射光进行强度调制，系统通过采集这些携带着温度信息的拉曼散射光，经过解调就可以得到传感光纤沿线的温度变化情况。分布式光纤拉曼传感系统具有环境适应能力强、抗电磁干扰、检测范围大、温度精度高等优点，因此被广泛应用于煤矿、油气管道、桥梁、建筑等温度安全监测领域。
在分布式光纤拉曼传感系统中，空间分辨率是一项主要的技术指标，可以反映测温系统中能够分辨光纤温度变化的最小长度。该系统所用探测信号为脉冲信号，定位原理为脉冲时间飞行法，这种方法采集得到的拉曼后向散射信号并非传感光纤L位置点的光强信息，而是大小为半个脉冲宽度的一段传感光纤上的光强信息叠加，这使得传统分布式拉曼传感系统空间分辨率受限于脉冲宽度。
基于此，有必要发明一种全新的温度解调方法，解决现有分布式光纤拉曼传感系统空间分辨率受限于OTDR原理，导致其难以突破1m的技术难题，以及系统解调温度前需要进行定标处理问题。
发明内容
为了解决现有分布式光纤拉曼传感系统空间分辨率受限于OTDR原理，导致其难以突破1m的技术难题，以及系统解调温度前需要进行定标处理的问题，本发明提出了一种基于ASE时域相关压缩的分布式光纤拉曼传感装置和方法，可以准确识别并解调出光纤沿线温度突变点信息且不需要定标过程，最终实现光纤沿线的毫米量级空间分辨率测量。
为了解决上述技术问题，本发明采用的技术方案为：一种基于ASE时域相关压缩的分布式光纤拉曼温度传感装置，包括ASE脉冲信号源，光纤耦合器，波分复用器，传感光纤，掺铒光纤放大器，第一光电探测器，第二光电探测器，数据采集卡和计算机；
所述ASE脉冲信号源发出的ASE脉冲信号经光纤耦合器分为两束，其中一束作为参考光束被第二探测器探测；另一束作为探测光束经波分复用器后入射到传感光纤，在传感光纤中产生的后向拉曼散射光经波分复用器、掺铒光纤放大器后，被第一探测器探测；数据采集卡采集第一光电探测器，第二光电探测器的探测信号后发送至所述计算机；
所述计算机用于对ASE拉曼后向散射信号进行时域差分信号重构，得到短尺度下的ASE拉曼后向散射信号，并通过将ASE脉冲参考信号和重构后的ASE拉曼后向散射信号进行时域相关压缩运算，进而定位和解调光纤沿线的温度突变点信息。
ASE脉冲信号源包括ASE信号发生器，隔离器，声光调制器，脉冲信号发生器，所述ASE信号发生器发出的ASE激光经隔离器后入射到声光调制器，经脉冲信号发生器驱动的声光调制器调制后输出ASE脉冲信号；
ASE信号发生器的工作波长为1550nm。
光纤耦合器用于将ASE脉冲信号分比为功率比为1∶99的参考光和探测光。
本发明还提供了一种基于ASE时域相关压缩的分布式光纤拉曼温度传感方法，采用所述的一种基于ASE时域相关压缩的分布式光纤拉曼温度传感装置实现，其特征在于，包括以下步骤：
S1、搭建传感装置，并进行ASE脉冲参考信号Iref和ASE拉曼后向散射信号Ias(L)的数据采集；Ias(L)表示传感光纤L位置处的ASE拉曼后向散射信号强度；
S2、对ASE拉曼后向散射信号进行时域差分信号重构；
S3、根据ASE脉冲参考信号和时域差分信号重构后的ASE拉曼后向散射信号，进行短尺度时域相关压缩系数的运算；
S4、寻找短尺度时域相关压缩系数的正峰值和负峰值，分别对应的延时采样点个数a0和a1，并计算温度突变位置，温度突变区域长度和温度突变区域温度信息。
所述步骤S2中，对ASE拉曼后向散射信号进行时域差分信号重构的表达式为：
Figure BDA0003027579130000021
其中，Fas(L，i)表示传感光纤L位置处的ASE拉曼后向散射信号进行时域差分信号重构后第i个数据点的时域差分信号，Ias(L)和
Figure BDA0003027579130000022
分别表示传感光纤L位置处和/>
Figure BDA0003027579130000023
位置处的后向拉曼反斯托克斯光的强度，其中，1≤i≤W·fs。
所述步骤S3中，进行短尺度时域相关压缩系数的运算公式如下：
Figure BDA0003027579130000031
其中，C(a)表示延时a个采样点时的短尺度时域相关压缩系数；N表示重构后ASE拉曼后向散射信号总采样点个数，Fas(n)表示进行时域差分信号重构后第n个采样点的时域差分信号；Iref(n+a)表示延时a个采样点的ASE脉冲参考信号；N表示重构后ASE拉曼后向散射信号总采样点个数。
温度突变位置和温度突变区域长度的计算公式分别为：
L1＝a0·c/(2n0·fs)；
ΔL＝(a1-a0)·c/(2n0·fs)；
其中，a0和a1分别表示短尺度时域相关压缩系数的正峰值和负峰值分别对应的延时采样点个数，c表示光速，n0表示光纤中的折射率，fs表示采样率，L1和ΔL分别表示温度突变位置和温变区域长度。
温度突变区域温度信息的计算公式为：
Figure BDA0003027579130000032
其中，Cpeak表示短尺度时域相关压缩系数的正峰值，W表示脉冲宽度，fs表示采样率，Pi为ASE脉冲激光第i个数据点的功率，Ka表示与拉曼反斯托克斯背向散射截面有关的系数，λa是拉曼反斯托克斯散射信号的波长，Δv为拉曼频移，h为普朗克常数，k为波尔兹曼常数，T0表示非温变区传感光纤温度，
Figure BDA0003027579130000033
表示温变位置L1处的附加损耗信息，/>
Figure BDA0003027579130000034
表示温变位置L1处的衰减信息，/>
Figure BDA0003027579130000035
α0、αas分别是入射光和反斯托克斯光在传感光纤中单位长度上的损耗系数，ΔAp(m)表示温度突变区域长度对正相关峰值的影响函数，其计算公式为：
Figure BDA0003027579130000041
其中，m表示温度突变区域的长度ΔL在时间尺度内的表示，m＝ΔL·2n0/c，ΔL表示温变区域长度，c表示光速，n0表示光纤中的折射率。
本发明与现有技术相比具有以下有益效果：本发明基于自发辐射噪声(ASE)激光信号产生拉曼后向散射信号，将其进行时域差分重构，在光纤沿线温度突变区域得到与ASE激光信号时序特征一致的拉曼后向散射信号，随后将重构的拉曼后向散射信号与ASE激光信号进行短尺度时域相关压缩处理，根据相关峰值定位和解调光纤沿线的温度突变点信息，时域差分重构方法消除了脉冲宽度对空间分辨率的限制，使得系统可以达到毫米量级的空间分辨率；短尺度时域相关压缩方法提高了系统的信噪比，提升了温度精度且无需定标过程。
附图说明
图1为本发明实施例提供的一种基于ASE时域相关压缩的分布式光纤拉曼传感装置的结构示意图；
图中：1-ASE信号发生器，2-隔离器，3-声光调制器，4-脉冲信号发生器，5-光纤耦合器，6-波分复用器，7-传感光纤，8-掺铒光纤放大器，9-第一光电探测器，10-第二光电探测器，11-高速数据采集卡，12-计算机。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例；基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例一
如图1所示，本发明实施例提供了一种基于ASE时域相关压缩的分布式光纤拉曼传感装置，其通过将ASE脉冲信号源输出的ASE脉冲激光分为两束，一束作为探测光束入射至传感光纤，在传感光纤中产生的背向ASE拉曼散射光中的反斯托克斯光被第一探测器探测，另一束作为参考光束被第二探测器探测；通过数据采集卡采集两个探测器的探测信号，然后其基于时域差分信号重构方法对后向拉曼散射信号进行重构，短尺度时域相关压缩方法定位和解调光纤沿线温度突变点温度信息。
具体地，如图1所示，本实施例的一种基于ASE时域相关压缩的分布式光纤拉曼温度传感装置包括ASE激光发生器1，隔离器2，声光调制器3，脉冲信号发生器4，光纤耦合器5，波分复用器6，传感光纤7，掺铒光纤放大器8，第一光电探测器9，第二光电探测器10，高速数据采集卡11和计算机12。
ASE激光器1的输出端与隔离器2输入端连接；隔离器2输出端与声光调制器3输入端连接，脉冲信号发生器4与声光调制器3连接，声光调制器3输出端与光纤耦合器5的a端口连接；波分复用器6、第二光电探测器10分别与光纤耦合器5的b端口、c端口连接；波分复用器6的b端口、c端口分别与传感光纤7、掺铒光纤放大器8连接；掺铒光纤放大器8的输出端与第一光电探测器9连接；第一光电探测器9、第二光电探测器10的输出端与高速数据采集卡11连接；数据采集卡11与计算机12连接。
所述ASE激光发生器1发出ASE激光信号经隔离器2后入射到声光调制器3、脉冲信号发生器4驱动的声光调制器3将ASE激光信息调制成变成脉冲信号，然后经光纤耦合器5分为两束，其中一束作为参考光束被第二探测器10探测；另一束作为探测光束经波分复用器6后入射到传感光纤7，在传感光纤7中产生的背向拉曼散射光经波分复用器6、掺铒光纤放大器8，被第一探测器9探测；通过数据采集卡11采集两个探测器的探测信号并发送至所述计算机12。
进一步地，ASE激光器1的工作波长为1550nm，脉冲激光经过分支比为1：99的光纤耦合器3后，分为功率比为1:99的参考路和探测路；探测路的ASE脉冲激光经波分复用器后进入传感光纤7，在传感光纤7中产生后向拉曼散射光从波分复用器6的1450nm(c)端口输出。
所述计算机12用于对ASE拉曼后向散射信号进行时域差分信号重构，得到短尺度下的ASE拉曼后向散射信号，并通过将ASE脉冲参考信号和重构ASE拉曼后向散射信号进行时域相关压缩运算，进而定位和解调光纤沿线的温度突变点信息。
实施例二
本实施例提供了一种基于ASE时域相关压缩的分布式光纤拉曼温度传感方法，其基于实施例的一种装置实现，通过时域差分信号重构方法对后向拉曼散射信号进行重构，短尺度时域相关压缩方法定位和解调光纤沿线温度突变点温度信息，其具体包括以下步骤：
S1、搭建传感装置，并进行ASE脉冲参考信号Iref和ASE拉曼后向散射信号Ias(L)的数据采集；Ias(L)表示传感光纤L位置处的ASE拉曼后向散射信号强度；
S2、对ASE拉曼后向散射信号进行时域差分信号重构；
S3、根据ASE脉冲参考信号和时域差分信号重构后的ASE拉曼后向散射信号，进行短尺度时域相关压缩系数的运算；
S4、寻找短尺度时域相关压缩系数的正峰值和负峰值分别对应的延时采样点个数a0和a1，并计算温度突变位置，温度突变区域长度和温度突变区域温度信息。
下面介绍本实施例的传感原理。
一、ASE脉冲参考信号、拉曼后向散射信号光强采集处理
(1)ASE拉曼后向散射信号的采集与处理。
在温度解调中，设激光脉冲宽度为W，传感光纤L位置处的后向拉曼散射信号(anti-Stokes)强度为：
Ias(L)＝Ka·λa -4·P·Ras(T)·exp[-(α0+αas)·L] (1)
式中P为脉冲激光的入射功率，Ka表示与拉曼反斯托克斯背向散射截面有关的系数，λa是拉曼反斯托克斯散射信号的波长，α0、αas分别是入射光和反斯托克斯光在传感光纤中单位长度上的损耗系数，L为传感光纤的位置，Ras(T)是反斯托克斯散射光的温度调制函数：
Figure BDA0003027579130000061
其中，Δv为拉曼频移，h为普朗克常数，k为波尔兹曼常数，T为传感光纤温度。
实际上，在分布式光纤拉曼传感系统中，由于光时域反射原理和时间飞行法定位原理，采集系统采集到传感光纤某一位置的拉曼反斯托克斯信号并不是传感光纤该点受温度调制的光强信息，而是脉宽尺度内传感光纤沿线所激发的拉曼散射光强信息累加。基于ASE信号的分布式光纤拉曼传感系统，ASE信号具有随机起伏特性。因此，将脉宽为W的ASE脉冲信号注入传感光纤，高速数据采集卡采集到传感光纤L位置处的拉曼后向散射信号强度(反斯托克斯光)可以表示为：
Figure BDA0003027579130000062
其中，Pi为ASE脉冲激光第i个数据点的功率，fs为采样率，W·fs为ASE脉冲信号在fs采样率下的单位脉冲数据点个数，Wi为i个数据点对应的脉冲宽度，Wi＝i/fs，
Figure BDA0003027579130000063
为位置L-Wi·c/2n0处的反斯托克斯散射光的温度调制函数值，表征了单位长度受温度调制的反斯托克斯光子后向散射因子，c为光在传感光纤中的传播速度，n0为传感光纤的折射率。
(2)ASE脉冲参考信号的采集与处理。
ASE信号发生器1发出的ASE脉冲激光经过光纤耦合器5的c端口，由采集卡11收集得到参考信号为:
Iref(i)＝Pi,1≤i≤W·fs； (4)
(3)基于时域差分信号重构方法重新构建拉曼后向散射信号。
采集得到的拉曼后向散射信号，其每个采样点均为整个ASE脉冲序列的光强信息叠加且受传感光纤温度突变区域调制的ASE信号不断发生变化。基于此，将后向散射信号进行重构，时域差分重构方法表达式为：
Figure BDA0003027579130000071
其中，Fas(L，i)为重构ASE拉曼后向散射信号。
设温变区温度为T1，非温变区的温度为T0，位置为L1，ΔL为温度突变区域的长度(空间尺度)，其长度大小在时间尺度内表示为m，其中m＝ΔL·2n0/c，经时域差分重构处理后，得到非温度突变区域重构ASE拉曼后向散射信号为0，温度突变点位置的重构拉曼后向散射信号Fas(L1,i)表示为：
Figure BDA0003027579130000072
其中，
Figure BDA0003027579130000073
Figure BDA0003027579130000074
是附加损耗信息，可以由公式(7)表示:
Figure BDA0003027579130000075
其中，
Figure BDA0003027579130000076
通过公式(3)可以看出，采集到的每个信号点实际为脉冲尺度长度的光纤叠加，本发明实施例中，通过公式(5)时域的差分操作，将公式(3)所示的反斯托克斯信号进行差分处理，时域差分之后的反斯托克斯信号如公式(6)所示，对比公式(3)和(6)可看出，经过时域差分处理后的反斯托克斯信号不再是信号累加，这样就消除了脉冲宽度限制的空间分辨率问题，提升了空间分辨率，然后再做自相关根据相关峰确定出温变区的位置和温度，因此，本发明的空间分辨率取决于相关峰的半高全宽。
二、利用短尺度时域相关压缩方法进行定位。
将ASE脉冲参考信号和重构后的ASE拉曼后向散射信号进行短尺度时域相关压缩运算，运算公式如下：
Figure BDA0003027579130000081
其中C(a)为短尺度时域相关压缩系数公式，N表示重构ASE拉曼后向散射信号总采样点个数，Fas(n)表示重构后的ASE拉曼后向散射信号的第n个采样点，Iref(n+a)表示延时a个采样点的ASE脉冲参考信号。
当参考信号延时a0个采样点到温度突变区域起始端点位置时相关系数呈现正峰值，因此可以确定温度突变区域位置为：
L1＝a0·c/(2n0·fs)。 (9)
当参考信号延时a1个采样点到温度突变区域另一端点位置时相关系数呈现负峰值，根据两峰值延时差值即可确定温度突变区域长度，即温度突变区域长度的计算公式为：
ΔL＝(a1-a0)·c/(2n0·fs)。 (10)
三、利用短尺度时域相关压缩正峰值解调温度。
将参考信号与温度突变位置的重构后的ASE拉曼后向散射信号进行相关短尺度时域压缩，相关正峰值表达式为：
Figure BDA0003027579130000082
为了分析温度突变区域长度对相关峰值的影响，引入参考信号自相关函数A(m)，m≤W，m为温度突变区域的空间尺度长度大小在时间尺度内表示，其表达式为：
Figure BDA0003027579130000083
定义温度突变区域长度对正相关峰值的影响函数ΔAp(m)，表示为：
Figure BDA0003027579130000084
结合公式(4)(6)，对公式(11)进行展开，并结合式(13)，可以得到：
Figure BDA0003027579130000085
结合Ras(T)＝[exp(hΔν/kT)-1]-1提取温度突变区域温度信息，其解调方程如公式(15)所示。
Figure BDA0003027579130000091
综上所述，本发明将ASE拉曼后向散射信号进行时域差分重构，在光纤沿线温度突变区域得到与ASE探测信号时序特征一致的ASE拉曼反斯托克斯信号，随后将重构的ASE拉曼后向散射信号与ASE参考信号进行短尺度时域相关压缩处理，根据相关峰值定位和解调光纤沿线的温度突变点信息。时域差分重构方法消除了脉冲宽度对空间分辨率的限制，使得系统可以达到毫米量级的空间分辨率。本发明通过短尺度时域相关压缩方法提高了系统的信噪比，提升了温度精度且无需定标过程。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。

Claims (6)

1.一种基于ASE时域相关压缩的分布式光纤拉曼温度传感装置，其特征在于，包括ASE脉冲信号源，光纤耦合器（5），波分复用器（6），传感光纤（7），掺铒光纤放大器（8），第一光电探测器（9），第二光电探测器（10），数据采集卡（11）和计算机（12）；
所述ASE脉冲信号源发出的ASE脉冲信号经光纤耦合器（5）分为两束，其中一束作为参考光束被第二光电探测器（10）探测；另一束作为探测光束经波分复用器（6）后入射到传感光纤（7），在传感光纤（7）中产生的后向拉曼散射光经波分复用器（6）、掺铒光纤放大器（8）后，被第一光电探测器（9）探测；数据采集卡（11）采集第一光电探测器（9），第二光电探测器（10）的探测信号后发送至所述计算机（12）；
所述计算机（12）用于对ASE拉曼后向散射信号进行时域差分信号重构，得到短尺度下的ASE拉曼后向散射信号，并通过将ASE脉冲参考信号和重构后的ASE拉曼后向散射信号进行时域相关压缩运算，进而定位和解调光纤沿线的温度突变点信息；
对ASE拉曼后向散射信号进行时域差分信号重构的表达式为：
Figure QLYQS_1
；
其中，F as（L,i）表示传感光纤L位置处的ASE拉曼后向散射信号进行时域差分信号重构后第i个数据点的时域差分信号，
Figure QLYQS_2
和/>
Figure QLYQS_3
分别表示传感光纤L位置处和/>
Figure QLYQS_4
位置处的后向拉曼反斯托克斯光的强度，其中，/>
Figure QLYQS_5
，W表示脉冲宽度，f s表示采样率，c表示光速，n0表示光纤中的折射率；
进行短尺度时域相关压缩系数的运算公式如下：
Figure QLYQS_6
其中，C（a）表示延时a个采样点时的短尺度时域相关压缩系数；N表示重构后ASE拉曼后向散射信号总采样点个数，F as（n）表示进行时域差分信号重构后第n个采样点的时域差分信号；I ref (n+a)表示延时a个采样点的ASE脉冲参考信号；N表示重构后ASE拉曼后向散射信号总采样点个数。
2.根据权利要求1所述的基于ASE时域相关压缩的分布式光纤拉曼温度传感装置，其特征在于，ASE脉冲信号源包括ASE信号发生器（1），隔离器（2），声光调制器（3），脉冲信号发生器（4），所述ASE信号发生器（1）发出的ASE激光经隔离器（2）后入射到声光调制器（3），经脉冲信号发生器（4）驱动的声光调制器（3）调制后输出ASE脉冲信号；
ASE信号发生器（1）的工作波长为1550nm。
3.根据权利要求1所述的基于ASE时域相关压缩的分布式光纤拉曼温度传感装置，其特征在于，光纤耦合器（5）用于将ASE脉冲信号分比为功率比为1：99的参考光和探测光。
4.一种基于ASE时域相关压缩的分布式光纤拉曼温度传感方法，采用权利要求1所述的一种基于ASE时域相关压缩的分布式光纤拉曼温度传感装置实现，其特征在于，包括以下步骤：
S1、搭建传感装置，并进行ASE脉冲参考信号Iref和ASE拉曼后向散射信号Ias（L）的数据采集；Ias（L）表示传感光纤L位置处的ASE拉曼后向散射信号强度；
S2、对ASE拉曼后向散射信号进行时域差分信号重构；
S3、根据ASE脉冲参考信号和时域差分信号重构后的ASE拉曼后向散射信号，进行短尺度时域相关压缩系数的运算；
S4、寻找短尺度时域相关压缩系数的正峰值和负峰值，分别对应的延时采样点个数a0和a1，并计算温度突变位置，温度突变区域长度和温度突变区域温度信息。
5.根据权利要求4所述的一种基于ASE时域相关压缩的分布式光纤拉曼温度传感方法，其特征在于，温度突变位置和温度突变区域长度的计算公式分别为：
Figure QLYQS_7
；
Figure QLYQS_8
；
其中，a0和a1分别表示短尺度时域相关压缩系数的正峰值和负峰值分别对应的延时采样点个数，c表示光速，n0表示光纤中的折射率，f s表示采样率，L1和ΔL分别表示温度突变位置和温变区域长度。
6.根据权利要求4所述的一种基于ASE时域相关压缩的分布式光纤拉曼温度传感方法，其特征在于，温度突变区域温度信息的计算公式为：
Figure QLYQS_9
；
其中，Cpeak表示短尺度时域相关压缩系数的正峰值，W表示脉冲宽度，f s表示采样率，P i 为ASE脉冲激光第i个数据点的功率，K a 表示与拉曼反斯托克斯背向散射截面有关的系数，λa是拉曼反斯托克斯散射信号的波长，Δν为拉曼频移，h为普朗克常数，k为波尔兹曼常数，T 0 表示非温变区传感光纤温度，Δφ(L 1)表示温变位置L 1处的附加损耗信息，φ(L 1)表示温变位置L 1 处的衰减信息，
Figure QLYQS_10
，α 0 、α as 分别是入射光和反斯托克斯光在传感光纤中单位长度上的损耗系数，/>
Figure QLYQS_11
表示温度突变区域长度对正相关峰值的影响函数，其计算公式为：
Figure QLYQS_12
其中，m表示温度突变区域的长度ΔL在时间尺度内的表示，
Figure QLYQS_13
，ΔL表示温变区域长度，c表示光速，n0表示光纤中的折射率。
CN202110420282.0A 2021-04-19 2021-04-19 Ase时域相关压缩的分布式光纤拉曼温度传感装置和方法 Active CN113091947B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110420282.0A CN113091947B (zh) 2021-04-19 2021-04-19 Ase时域相关压缩的分布式光纤拉曼温度传感装置和方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110420282.0A CN113091947B (zh) 2021-04-19 2021-04-19 Ase时域相关压缩的分布式光纤拉曼温度传感装置和方法

Publications (2)

Publication Number Publication Date
CN113091947A CN113091947A (zh) 2021-07-09
CN113091947B true CN113091947B (zh) 2023-06-30

Family

ID=76678601

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110420282.0A Active CN113091947B (zh) 2021-04-19 2021-04-19 Ase时域相关压缩的分布式光纤拉曼温度传感装置和方法

Country Status (1)

Country Link
CN (1) CN113091947B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113776691B (zh) * 2021-08-24 2023-05-26 太原理工大学 面向煤田采空区火源钻孔探测的分布式光纤温度传感方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB9307660D0 (en) * 1993-04-14 1993-06-02 York Ltd Distributed sensing
CN108760080A (zh) * 2018-05-18 2018-11-06 太原理工大学 一种基于ase噪声的分布式光纤拉曼测温装置及方法
CN110608760A (zh) * 2019-07-05 2019-12-24 中国地质大学（武汉） 相位敏感光时域反射传感系统扰动探测信噪比提高方法
CN112378432A (zh) * 2020-10-29 2021-02-19 太原理工大学 基于噪声匹配滤波的分布式光纤拉曼传感装置和方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101825498B (zh) * 2010-04-13 2012-03-21 中国计量学院 色散与损耗光谱自校正分布式光纤拉曼温度传感器

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB9307660D0 (en) * 1993-04-14 1993-06-02 York Ltd Distributed sensing
CN108760080A (zh) * 2018-05-18 2018-11-06 太原理工大学 一种基于ase噪声的分布式光纤拉曼测温装置及方法
CN110608760A (zh) * 2019-07-05 2019-12-24 中国地质大学（武汉） 相位敏感光时域反射传感系统扰动探测信噪比提高方法
CN112378432A (zh) * 2020-10-29 2021-02-19 太原理工大学 基于噪声匹配滤波的分布式光纤拉曼传感装置和方法

Also Published As

Publication number Publication date
CN113091947A (zh) 2021-07-09

Similar Documents

Publication Publication Date Title
CN107917738B (zh) 一种同时测量温度、应变和振动的分布式光纤传感系统
CN112378432B (zh) 基于噪声匹配滤波的分布式光纤拉曼传感装置和方法
CN108303197B (zh) 基于后向散射增强光纤分布式温度和应变双参量传感装置及其解调方法
CN110132329B (zh) 应力、温度和振动复合检测光纤传感器及信号处理方法
CN112378431B (zh) 基于宽频混沌激光的分布式光纤拉曼传感方法
JP2018146371A (ja) 温度・歪センシング装置及び温度・歪センシング方法
CN111879436B (zh) 基于双脉冲调制的分布式光纤拉曼温度解调装置及方法
CN107843357B (zh) 基于拉曼散射的分布式光纤温度及应变检测方法
US20100014071A1 (en) Frequency-scanned optical time domain reflectometry
EP0225023A2 (en) Temperature measurement
CN105067104A (zh) 一种复合型光纤传感系统和传感方法
CN111896137B (zh) 一种厘米级空间分辨率的分布式光纤拉曼传感装置及方法
CN113091947B (zh) Ase时域相关压缩的分布式光纤拉曼温度传感装置和方法
CN113091946B (zh) 混沌分布式光纤拉曼温度传感装置及方法
CN109556754A (zh) 光纤应变和温度测定装置以及光纤应变和温度测定方法
CN117232558A (zh) 一种局部受激布里渊散射光时域反射系统及数据处理方法
CN111638025B (zh) 面向隧道渗漏水监测的分布式光纤拉曼传感装置及方法
EP4206624A1 (en) Vibration distribution measuring device and method of same
CN108344432B (zh) 获取布里渊光纤分布式传感系统的传感信息的方法
CN112033572B (zh) 一种光时域反射计中补偿弯曲损耗及斯托克斯光与反斯托克斯光衰减的方法
CN110440837B (zh) 一种多参量光纤同步传感采集仪和传感采集方法
CN112697303B (zh) 面向智能电网的分布式光纤传感系统和检测方法
CN219914342U (zh) 一种长距离单多模混合光纤工程应变监测系统
CN117213659A (zh) 基于布里渊和瑞利散射的混沌botdr装置及方法
CN113566860B (zh) 基于布里渊和拉曼三阶联合放大的混沌拉曼光纤传感装置

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