CN112985639B - 基于去啁啾和时域子啁啾脉冲提取的分布式光纤传感方法 - Google Patents
基于去啁啾和时域子啁啾脉冲提取的分布式光纤传感方法 Download PDFInfo
- Publication number
- CN112985639B CN112985639B CN202110166123.2A CN202110166123A CN112985639B CN 112985639 B CN112985639 B CN 112985639B CN 202110166123 A CN202110166123 A CN 202110166123A CN 112985639 B CN112985639 B CN 112985639B
- Authority
- CN
- China
- Prior art keywords
- chirp
- rayleigh scattering
- optical fiber
- signal
- pattern
- 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
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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Optical Transform (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了基于去啁啾和时域子啁啾脉冲提取的分布式光纤传感方法,涉及光纤传感测量技术领域,解决现有分布式光纤传感***受限于接收机带宽而难以实现大测量范围、高空间分辨率测量的技术问题。本发明具体为使用去啁啾技术获取大频率范围的瑞利散射信号;使用时域子啁啾脉冲提取算法将获取到的大频率范围的瑞利散射信号转换为大应变范围瑞利散射图样,光纤在未受扰动时得到瑞利散射参考图样,受扰动后得到瑞利散射测量图样;然后通过延迟估计算法,分析瑞利散射参考图样和瑞利散射测量图样在不同位置下频率轴的移动信息,定量的解调出不同位置的扰动信息,本申请适用于动态和静态传感、对空间分辨率牺牲较小、测量范围大、测量精度高。
Description
技术领域
本发明涉及光纤传感测量技术领域,具体是指基于去啁啾和时域子啁啾脉冲提取的分布式光纤传感方法。
背景技术
分布式光纤传感***目前已广泛用于地质探矿、结构健康监测、轨道交通在线监测等领域。目前基于分布式光纤传感的定量测量有几种典型的方案:1)相位解调型:通过解调光时\频域反射计瑞利散射信号的相位信息,得到沿光纤温度或应变变化的定量信息。但这种方法受到干涉衰落影响非常严重,且测量范围由于三角函数的周期性而被限制在2π之内。2)基于啁啾脉冲相位敏感型光时域反射计:通过计算扰动引起的啁啾脉冲时移来来定量计算光纤扰动信息。这种方法不受干涉衰落影响,但是其空间分辨率取决于脉冲宽度,所以和传感距离、测量灵敏度相互限制。此外,这种方法需要很高的采样率来保证量化误差。3)基于瑞利散射图样解调。这种方法可以实现很高的测量精度,没有干涉衰落,且测量范围取决于扫频范围,从而可以通过增大扫频范围来提升。但传统的基于相干光时域反射计的方法需要多次扫频来获取瑞利散射图样,使得其测量速度极为受限。基于子啁啾脉冲提取的相干光时域反射计可以避免扫频过程实现瑞利散射图样的单次获取。但由于电学设备如探测器带宽,ADC采样率,以及后端数字信号处理器的处理能力等的限制,过大的频率范围会使得接收机后端压力过大,甚至无法处理。
发明内容
基于以上问题,本发明提供了基于去啁啾和时域子啁啾脉冲提取的分布式光纤传感方法,适用于动态和静态传感、空间分辨率高、测量范围大、测量精度高。
为解决以上技术问题,本发明采用的技术方案如下:
基于去啁啾和时域子啁啾脉冲提取的分布式光纤传感方法,包括如下步骤:
步骤一、使用去啁啾技术获取大频率范围的瑞利散射信号;
步骤二、使用时域子啁啾脉冲提取算法将获取到的大频率范围的瑞利散射信号转换为大应变范围瑞利散射图样;其中在光纤未受扰动时,获取的瑞利散射图样叫瑞利散射参考图样;在光纤受到扰动时,获取的瑞利散射图样叫瑞利散射测量图样;
步骤三、通过延迟估计算法,分析瑞利散射参考图样和瑞利散射测量图样在不同位置下频率轴的移动信息,定量的解调出不同位置的扰动信息。
进一步,所述使用去啁啾技术获取大频率范围的瑞利散射信号的过程包括:设置扫频光源输出大频率范围的探测啁啾信号,并将其注入光纤,得到反向瑞利散射信号;将得到的大频率范围瑞利散射信号与另一个本地啁啾信号混频,得到去啁啾之后的信号。
进一步,所述另一个本地啁啾信号,为啁啾斜率与探测啁啾信号的啁啾斜率相同的信号,或为与啁啾斜率不同的信号。
进一步,使用时域子啁啾脉冲提取算法将获取到的大频率范围的瑞利散射信号转换为大应变范围瑞利散射图样,其具体过程为:将去啁啾之后的信号用一系列的时域窗分割,然后按照时域窗的顺序重组为矩阵;将该矩阵转换到距离域;取距离域的强度值或幅度值,即可得到大应变范围瑞利散射图样。
进一步,所述时域窗为矩形窗或汉宁窗或汉明窗或高斯窗或切比雪夫窗,所述相邻两个滑动窗之间有重叠。
进一步,所述通过延迟估计算法来获取扰动信息,其过程为:计算瑞利散射参考图样和瑞利散射测量图样在不同位置下频率轴的互相关曲线或加权互相关曲线;得到的互相关曲线的最大值对应的横坐标即为频率轴的移动信息;根据瑞利散射温度或者应变对应的频移关系,即可定量地得到沿光纤的温度或者应变信息。
进一步,所述加权互相关算法,权重可以由Roth处理器或相干平滑变换或相位变换或Eckart滤波器或HT处理器得到。
进一步,所述通过延迟估计算法来获取扰动信息,其过程为:计算瑞利散射参考图样和瑞利散射测量图样在不同位置下频率轴的最小均方误差,得到的最小均方误差曲线的最小值对应的横坐标即为频率轴的移动信息;根据瑞利散射温度或应变对应的频移关系,即可定量地得到沿光纤的温度或应变变化信息。
进一步,通过插值的方式来减小估计过程中的量化误差。
与现有技术相比,本发明的有益效果是:本发明通过去啁啾技术可以极大地降低***对接收机的需求。通过子啁啾脉冲提取技术,获取大频率范围的瑞利散射图样,从而实现动态和静态传感、空间分辨率高、测量范围大、测量精度高的分布式声波传感,且不会受到瑞利散射的相干衰落的影响。
附图说明
图1为本实施例的流程图;
图2为实施例1中去啁啾技术的示意图;
图3为实施例1中对得到的瑞利散射信号加一系列的时域窗;
图4为实施例1中得到的瑞利散射图样。
具体实施方式
下面结合附图对本发明作进一步的说明。本发明的实施方式包括但不限于下列实施例。
实施例1
如图1所示,通过去啁啾技术获取大扫频范围啁啾信号的瑞利散射信号,使用时域子啁啾脉冲提取算法获取瑞利散射图样,以相同的方法,当光纤未受扰动时,获取瑞利散射参考图样;当光纤受到扰动时,获取瑞利散射测量图样,然后利用延迟估计算法分析不同位置下频率轴的移动情况,即可获取沿光纤的扰动信息,具体方法如下:
步骤一、获取大扫频范围啁啾信号的瑞利散射信号;
首先设置扫频光源的扫频范围Bp、起始频率f0和扫频斜率k,将该啁啾信号注入到光纤,其瑞利散射信号如图2的矩形框部分所示,使用一个本地啁啾信号,如图2的线条所示与瑞利散射信号拍频,从图2中可以看出,拍频后的信号带宽为Bmea,其值远小于啁啾信号的扫频范围,从而使得接收机的压力被极大地减小。
步骤二、将去啁啾技术得到的信号通过时域子啁啾脉冲提取算法,重构成与频率、距离相关的瑞利散射图样,并获取瑞利散射参考图样和瑞利散射测量图样;
如图3所示,使用一系列的时域窗分割去啁啾得到的信号,然后按照时域窗的顺序重组为矩阵;将该矩阵转换到距离域;取距离域的强度值或幅度值,即可得到大应变范围瑞利散射图样,其中,在光纤未受到扰动时,得到的为瑞利散射参考图样在光纤受到扰动时,得到的为瑞利散射测量图样
其中,延迟估计算法为通过计算瑞利散射参考图样和瑞利散射测量图样在不同位置下频率轴的互相关曲线或加权互相关曲线;得到的互相关曲线的最大值对应的横坐标即为频率轴的移动信息;根据瑞利散射温度或者应变对应的频移关系,即可定量地得到沿光纤的温度或者应变信息,另外,在加权互相关算法中,权重可以由Roth处理器或相干平滑变换或相位变换或Eckart滤波器或HT处理器得到。另外,可通过插值的方式来减小估计过程中的量化误差。
实施例2
本实施例中采用另一种方式获取扰动信息,具体为:
计算瑞利散射参考图样和瑞利散射测量图样在不同位置下频率轴的最小均方误差,得到的最小均方误差曲线的最小值对应的横坐标即为频率轴的移动信息;根据瑞利散射温度或者应变对应的频移关系,即可定量地得到沿光纤的温度或者应变信息,另外,可通过插值的方式来减小估计过程中的量化误差。
如上即为本发明的实施例。上述实施例以及实施例中的具体参数仅是为了清楚表述发明人的发明验证过程,并非用以限制本发明的专利保护范围,本发明的专利保护范围仍然以其权利要求书为准,凡是运用本发明的说明书及附图内容所作的等同结构变化,同理均应包含在本发明的保护范围内。
Claims (8)
1.基于去啁啾和时域子啁啾脉冲提取的分布式光纤传感方法,其特征在于,包括如下步骤:
步骤一、使用去啁啾技术获取大频率范围的瑞利散射信号;
步骤二、使用时域子啁啾脉冲提取算法将获取到的大频率范围的瑞利散射信号转换为大应变范围的瑞利散射图样;其中在光纤未受扰动时,获取的瑞利散射图样叫瑞利散射参考图样;在光纤受到扰动时,获取的瑞利散射图样叫瑞利散射测量图样;
步骤三、通过延迟估计算法,分析瑞利散射参考图样和瑞利散射测量图样在不同位置下频率轴的移动信息,定量的解调出不同位置的扰动信息;
使用时域子啁啾脉冲提取算法将获取到的大频率范围的瑞利散射信号转换为大应变范围瑞利散射图样,其具体过程为:将去啁啾之后的信号用一系列的时域窗分割,然后按照时域窗的顺序重组为矩阵;将该矩阵转换到距离域;取距离域的强度值或幅度值,即可得到大应变范围瑞利散射图样。
2.根据权利要求1所述的基于去啁啾和时域子啁啾脉冲提取的分布式光纤传感方法,其特征在于:所述使用去啁啾技术获取大频率范围的瑞利散射信号的过程包括:设置扫频光源输出大频率范围的探测啁啾信号,并将其注入光纤,得到反向瑞利散射信号;将得到的大频率范围瑞利散射信号与另一个本地啁啾信号混频,得到去啁啾之后的信号。
3.根据权利要求2所述的基于去啁啾和时域子啁啾脉冲提取的分布式光纤传感方法,所述另一个本地啁啾信号,为啁啾斜率与探测啁啾信号的啁啾斜率相同的信号,或为与啁啾斜率不同的信号。
4.根据权利要求1所述的基于去啁啾和时域子啁啾脉冲提取的分布式光纤传感方法,其特征在于:所述时域窗为矩形窗或汉宁窗或汉明窗或高斯窗或切比雪夫窗。
5.根据权利要求1所述的基于去啁啾和时域子啁啾脉冲提取的分布式光纤传感方法,其特征在于:所述通过延迟估计算法来获取扰动信息,其过程为:计算瑞利散射参考图样和瑞利散射测量图样在不同位置下频率轴的互相关曲线或加权互相关曲线;得到的互相关曲线的最大值对应的横坐标即为频率轴的移动信息;根据瑞利散射温度或者应变对应的频移关系,即可定量地得到沿光纤的温度或者应变信息。
6.根据权利要求5所述的基于去啁啾和时域子啁啾脉冲提取的分布式光纤传感方法,其特征在于:加权互相关算法,权重可以由Roth处理器或相干平滑变换或相位变换或Eckart滤波器或HT处理器得到。
7.根据权利要求1所述的基于去啁啾和时域子啁啾脉冲提取的分布式光纤传感方法,其特征在于:所述通过延迟估计算法来获取扰动信息,其过程为:计算瑞利散射参考图样和瑞利散射测量图样在不同位置下频率轴的最小均方误差,得到的最小均方误差曲线的最小值对应的横坐标即为频率轴的移动信息;根据瑞利散射温度或者应变对应的频移关系,即可定量地得到沿光纤的温度或应变信息。
8.根据权利要求5或7所述的基于去啁啾和时域子啁啾脉冲提取的分布式光纤传感方法,其特征在于:通过插值的方式来减小估计过程中的量化误差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110166123.2A CN112985639B (zh) | 2021-02-06 | 2021-02-06 | 基于去啁啾和时域子啁啾脉冲提取的分布式光纤传感方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110166123.2A CN112985639B (zh) | 2021-02-06 | 2021-02-06 | 基于去啁啾和时域子啁啾脉冲提取的分布式光纤传感方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112985639A CN112985639A (zh) | 2021-06-18 |
CN112985639B true CN112985639B (zh) | 2022-02-15 |
Family
ID=76348609
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110166123.2A Active CN112985639B (zh) | 2021-02-06 | 2021-02-06 | 基于去啁啾和时域子啁啾脉冲提取的分布式光纤传感方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112985639B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113810098B (zh) * | 2021-07-20 | 2022-09-16 | 广东工业大学 | 一种基于双边带啁啾脉冲调制的光时域反射计 |
CN113607074B (zh) * | 2021-07-28 | 2023-08-01 | 天津大学前沿技术研究院有限公司 | 一种光频域反射中基于重叠窗口的应变传感解调方法 |
CN113776566A (zh) * | 2021-08-06 | 2021-12-10 | 电子科技大学 | 基于子脉冲提取算法的分布式温度应变传感方法 |
CN113810099B (zh) * | 2021-08-20 | 2022-10-14 | 广东工业大学 | 一种基于非对称双边带啁啾脉冲调制的光时域反射计 |
CN113465640B (zh) * | 2021-08-23 | 2024-05-17 | 中国科学院半导体研究所 | 基于啁啾脉冲的光纤温度压力传感*** |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107289978B (zh) * | 2017-06-09 | 2019-05-07 | 南京大学 | 一种基于potdr的测扰动的***及方法 |
CN108286992B (zh) * | 2018-01-06 | 2020-04-17 | 天津大学 | 基于数字双啁啾脉冲调制的分布式光纤声传感装置及方法 |
CN110132331B (zh) * | 2019-05-22 | 2021-05-07 | 电子科技大学 | 一种基于子啁啾脉冲提取的cotdr测量方法及其*** |
CN111829673B (zh) * | 2020-06-09 | 2021-08-03 | 中国地质大学(武汉) | 一种对称线性啁啾脉冲探测otdr的测量*** |
CN111609875B (zh) * | 2020-06-10 | 2021-12-28 | 电子科技大学 | 基于啁啾连续光的数字域可调分布式光纤传感***及方法 |
CN111854815B (zh) * | 2020-08-05 | 2022-04-01 | 国网四川省电力公司经济技术研究院 | 基于分数阶傅里叶变换的高速分布式光纤传感***及方法 |
-
2021
- 2021-02-06 CN CN202110166123.2A patent/CN112985639B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN112985639A (zh) | 2021-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112985639B (zh) | 基于去啁啾和时域子啁啾脉冲提取的分布式光纤传感方法 | |
US11193751B2 (en) | Method and apparatus for motion compensation in interferometric sensing systems | |
CN109084905B (zh) | 相位敏感型光时域反射***中降低探测死区概率的方法 | |
CN110132329B (zh) | 应力、温度和振动复合检测光纤传感器及信号处理方法 | |
US20160252651A1 (en) | Trace downsampling of distributed acoustic sensor data | |
CN108955734B (zh) | 一种光纤f-p温度/压力复合传感器的腔长解调方法 | |
Chen et al. | Denoising of distributed acoustic sensing seismic data using an integrated framework | |
CN103412137B (zh) | 旋转因子中和测速方法和装置 | |
Wang et al. | Rapid response DAS denoising method based on deep learning | |
Bai et al. | Noise reduction method of Φ-OTDR system based on EMD-TFPF algorithm | |
Zhang et al. | Anti-noise frequency estimation performance of Hanning-windowed energy centrobaric method for optical coherence velocimeter | |
CN115824378A (zh) | 高频响分布式光纤声波传感器的振动检测方法 | |
CN113514833B (zh) | 一种基于海浪图像的海面任意点波向的反演方法 | |
Zhang et al. | Enhancement for Φ-OTDR performance by using narrow linewidth light source and signal processing | |
CN114623921B (zh) | 分布式光纤振动传感的解调方法 | |
CN108490255A (zh) | 一种基于两级fft的数字宽带高精度测频方法及*** | |
CN111123286A (zh) | 一种基于自校准全光纤多普勒电缆长度测试方法及装置 | |
CN116243242A (zh) | 基于ф-otdr的自适应滤波一体化降噪高精度3d定位方法 | |
CN112179475B (zh) | 一种多源混叠的分布式光纤振动传感信号的分离方法 | |
CN107783084B (zh) | 恒虚警检测及数据处理平台 | |
CN115752696A (zh) | 基于前向传输的分布式光纤振动传感***及信号处理方法 | |
CN115345216A (zh) | 一种融合先验信息的fmcw雷达干扰消除方法 | |
CN109084816B (zh) | 一种光纤传感信号的抗偏振衰落报警算法 | |
CN207515900U (zh) | 基于脉冲编码和边沿滤波法的光纤光栅振动测量*** | |
CN117251676A (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 |