CN113237431A - 一种提升ofdr***分布式空间分辨率的测量方法 - Google Patents

一种提升ofdr***分布式空间分辨率的测量方法 Download PDF

Info

Publication number
CN113237431A
CN113237431A CN202110491617.8A CN202110491617A CN113237431A CN 113237431 A CN113237431 A CN 113237431A CN 202110491617 A CN202110491617 A CN 202110491617A CN 113237431 A CN113237431 A CN 113237431A
Authority
CN
China
Prior art keywords
spatial resolution
signal
optical fiber
image
ofdr
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
CN202110491617.8A
Other languages
English (en)
Other versions
CN113237431B (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.)
Shandong University
Original Assignee
Shandong 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 Shandong University filed Critical Shandong University
Priority to CN202110491617.8A priority Critical patent/CN113237431B/zh
Publication of CN113237431A publication Critical patent/CN113237431A/zh
Application granted granted Critical
Publication of CN113237431B publication Critical patent/CN113237431B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • G01B11/18—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge using photoelastic elements

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明公开了一种提升OFDR***分布式空间分辨率的测量方法，首先将传统的一维互相关处理结果信号转化为二维图像信号，在二维图像的基础上通过高斯滤波去噪算法对图像进行去噪处理，再将处理后的图像进行下一步运算，最终得到高空间分辨率的结果。本发明所提出的高斯滤波去噪算法的OFDR传感***可以提高测量***的空间分辨率，使其在航天航空、机器设备等高精度监测领域具有更大优势和更广阔的应用。

Description

一种提升OFDR***分布式空间分辨率的测量方法
技术领域
本发明涉及一种提升OFDR***分布式空间分辨率的测量方法，应用于应变测量领域，提高了***的空间分辨率，属于光纤传感探测的技术领域。
背景技术
分布式光纤传感中光纤作为传感介质的同时也是测量的传输介质，利用光纤中光波的传输特性，包括拉曼散射、瑞利散射和布里渊散射，实现了沿光纤长度方向对外界环境进行实时监控。分布式光纤传感技术具有抗电磁干扰能力强，结构相对简单，空间分辨率高，传感距离长等优势。基于其上述优势，该技术正被逐渐应用在越来越多的领域，例如检测大桥安全，土木工程的检测，隧道等地下的火灾警报，勘测地质等，在社会建设中发挥着重要的影响。作为分布式光纤传感***的代表，光频域反射技术(OFDR)具有重量轻、体积小、灵敏度高、抗电磁干扰性强和空间分辨率高等优点，可以连续测量沿光纤距离上的应变、振动、温度等外界的物理量变化。近些年，随着OFDR技术的发展，也实现了形状传感和声学传感的应用。
OFDR原理为：可调谐激光光源发出的线性扫频光经过耦合器分为两束，一束进入待测光纤，待测光纤的后向瑞利散射光返回形成信号光，与另一束参考光发生拍频干涉。对拍频信号采集并进行快速傅里叶变换处理，就可以得到沿传感光纤构建的距离域信息。在OFDR***测量中，需要采集一次没有外界影响的参考信号和光纤受影响的测试信号，将测试信号与参考信号进行互相关计算得到外界信息的变化。OFDR***具有高空间分辨率的特点，其***空间分辨率可以达到毫米量级，因此其在航天航空等高精度监测领域具有非常重要的应用。然而，当测量空间分辨率提高时，参考信号与测试信号的互相关相关性会极大降低，从而导致互相关结果产生多峰和假峰，得不到正确的结果。因此，如何有效提高OFDR***的空间分辨率是一个十分重要的研究方向。
发明内容
本发明的目的是提供一种提升OFDR***分布式空间分辨率的测量方法，此方法可以有效提高***的空间分辨率。其技术方案为，
一种提升OFDR***分布式空间分辨率的测量方法，包括以下步骤，
S1.分别采集两次信号，一次为不包含应变信息的信号，为参考信号；另一次为包含应变信息的信号，为测试信号；
S2.将参考信号和测试信号通过快速傅里叶变化映射到距离域上，并取窗口大小C将信号划分为N等分，其中窗口大小C决定了***的空间分辨率；
S3.对第一份参考信号和测量信号的局部距离域信息进行快速逆傅里叶变化。
S4.对快速逆傅里叶变换后的参考信号和测量信号通过互相关计算得到一维互相关结果；
S5.重复步骤S3-S4，得到光纤每一对应位置的互相关结果，将得到的所有一维互相关结果在光纤距离上重构成二维图像信号，在二维图像的基础上通过高斯滤波去噪算法后对图像进行去噪处理；
S6.将通过高斯滤波去噪算法处理后的图像重新解构得到光纤各位置的光谱偏移量，从而能够得到高空间分辨率下的测量结果，提高测量的准确性。
优选的，利用高斯滤波对步骤S5整幅重构的二维图像进行加权平均，每一个像素点由其本身和邻域内的其他像素值经过加权平均后得到。
优选的，步骤S6中光谱偏移量获得方法为，
用一个特定的模板去扫描图像中的每一个像素，用模板确定的邻域内像素的加权平均灰度值去替代模板中心像素点的值；选择gaussian高斯低通滤波器，使用fspecial函数，选用的模板为一个M*M矩阵，标准偏差为m；将步骤S5重构的二维图像信号通过高斯滤波处理，处理后的图像重新分解到光纤的对应位置，通过寻找主峰的偏移得到对应光纤位置光谱的偏移量。
有益效果
1)本发明所提出的高斯滤波去噪算法的OFDR传感***可以提高测量***的空间分辨率，使其在航天航空、机器设备等高精度监测领域具有更大优势和更广阔的应用。
2)本发明所提出的基于高斯滤波去噪算法的OFDR传感***通过对二维图像信息的去噪处理，不仅能够提高***空间分辨率，而且能够有效去除测量结果的异常值，提高测量的准确性。
附图说明
图1基为于高斯滤波去噪算法的高空间分辨率OFDR***处理流程。
图2为OFDR***原理图。
其中1-为可调谐激光器；2-为耦合器一；3-为耦合器二；4-为环形器；5-为马赫-曾德干涉仪；6为偏振控制器一；7-为偏振控制器二；8-为耦合器三；9-为偏振分束器；10-为平衡探测器；11-为采集卡；12-传感光纤；13-菲尼尔环。
图3为传感光纤10.1-10.7m处受100με未使用此方法的结果图，空间分辨率为0.4mm。
图4为传感光纤10.1-10.7m处受100με使用此方法的结果图，空间分辨率为0.4mm。
具体实施方式
下面结合附图1-4和具体实施例对技术作进一步说明，以助于理解本发明的内容。
一种提升OFDR***分布式空间分辨率的测量方法，包括以下步骤，
S1.分别采集两次信号，一次为不包含应变信息的信号，为参考信号；另一次为包含应变信息的信号，为测试信号。
S2.将参考信号与测试信号在距离域上按一定的窗口大小C分割成N等份；
S3.将参考信号与测试信号各自的每一份距离域信息使用快速逆傅里叶变换；
S4.将逆傅里叶变换后的参考信号与测试信号进行互相关运算，可以得出一维互相关结果；
S5.重复步骤S3-S4，得到光纤每一对应位置的互相关结果，将得到的所有一维互相关结果在光纤距离上重构成二维图像信号；
S6.高斯滤波是一种线性平滑滤波，适用于消除高斯噪声。
其原理为对于步骤S5中重构成的二维图像，对整幅重构的二维图像进行加权平均，每一个像素点由其本身和邻域内的其他像素值经过加权平均后得到。用一个特定的模板去扫描图像中的每一个像素，用模板确定的邻域内像素的加权平均灰度值去替代模板中心像素点的值。在matlab软件中使用fspecial函数，选择gaussian高斯低通滤波器，选用的模板为一个100*100矩阵，标准偏差为2，将步骤S5重构的二维图像信号通过高斯滤波处理，处理后的图像重新分解到光纤的对应位置，通过寻找主峰的偏移得到对应光纤位置光谱的偏移量。
S7.在二维图像的基础上通过高斯滤波去噪算法后对图像进行去噪处理，将通过高斯滤波去噪算法处理后的图像重新解构得到光纤各位置的光谱偏移量，从而能够得到高空间分辨率下的测量结果，提高测量的准确性。
图3为传感光纤10.1-10.7m处受100με未使用此方法的结果图，空间分辨率为0.4mm，可以看到其结果有许多异常值，无法得出沿光纤长度正确的应变分布结果。
图4为传感光纤10.1-10.7m处受100με使用此方法的结果图，此方法很好的消除了异常值，提高了***分辨率，得到了0.4mm空间分辨率的正确应变分布结果。
图2是OFDR***原理图。一种基于距离域补偿的OFDR传感***，包括：可调谐激光源的连续激光输出由耦合器一2(10/90光耦合器)分成两份，10％入射到一个非平衡的马赫-曾德干涉仪5，为采集卡11提供一个触发信号，其余部分的光进入耦合器二3；之后耦合器二3(1/99光耦合器)分成两个部分，其中1％的输出通过偏振控制器一6进行调整，使“p”和“s”光分量具有相同的功率，99％通过环形器4和偏振控制器二7进入到传感光纤12探测，菲尼尔环13用来抑制光纤末端菲尼尔反射。；然后瑞利散射信号与1％激光输出从耦合器三8(50/50光耦合器)相结合得到的干涉信号，通过偏振分束器分解成“p”和“s”分量；最后“p”和“s”光由采集卡采集。
以上所述仅为本申请的优选实施例而已，并不用于限制本申请，对于本领域的技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。

Claims (3)

1.一种提升OFDR***分布式空间分辨率的测量方法，其特征在于，包括以下步骤，
S1.分别采集两次信号，一次为不包含应变信息的信号，为参考信号；另一次为包含应变信息的信号，为测试信号；
S2.将参考信号和测试信号通过快速傅里叶变化映射到距离域上，并取窗口大小C将信号划分为N等分，其中窗口大小C决定了***的空间分辨率；
S3.对第一份参考信号和测量信号的局部距离域信息进行快速逆傅里叶变化；
S4.对快速逆傅里叶变换后的参考信号和测量信号通过互相关计算得到一维互相关结果；
S5.重复步骤S3-S4，得到光纤每一对应位置的互相关结果，将得到的所有一维互相关结果在光纤距离上重构成二维图像信号，在二维图像的基础上通过高斯滤波去噪算法后对图像进行去噪处理；
S6.将通过高斯滤波去噪算法处理后的图像重新解构来得到光纤各位置的光谱偏移量，从而能够得到高空间分辨率下的测量结果，提高测量的准确性。
2.根据权利要求1所述的一种提升OFDR***分布式空间分辨率的测量方法，其特征在于，利用高斯滤波对步骤S5整幅重构的二维图像进行加权平均，每一个像素点由其本身和邻域内的其他像素值经过加权平均后得到。
3.根据权利要求1所述的一种提升OFDR***分布式空间分辨率的测量方法，其特征在于，步骤S6中光谱偏移量获得方法为，
用一个特定的模板去扫描图像中的每一个像素，用模板确定的邻域内像素的加权平均灰度值去替代模板中心像素点的值；选择gaussian高斯低通滤波器，使用fspecial函数，选用的模板为一个M*M矩阵，标准偏差为m；将步骤S5重构的二维图像信号通过高斯滤波处理，处理后的图像重新分解到光纤的对应位置，通过寻找主峰的偏移得到对应光纤位置光谱的偏移量。
CN202110491617.8A 2021-05-06 2021-05-06 一种提升ofdr***分布式空间分辨率的测量方法 Active CN113237431B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110491617.8A CN113237431B (zh) 2021-05-06 2021-05-06 一种提升ofdr***分布式空间分辨率的测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110491617.8A CN113237431B (zh) 2021-05-06 2021-05-06 一种提升ofdr***分布式空间分辨率的测量方法

Publications (2)

Publication Number Publication Date
CN113237431A true CN113237431A (zh) 2021-08-10
CN113237431B CN113237431B (zh) 2022-03-18

Family

ID=77132115

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110491617.8A Active CN113237431B (zh) 2021-05-06 2021-05-06 一种提升ofdr***分布式空间分辨率的测量方法

Country Status (1)

Country Link
CN (1) CN113237431B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114199514A (zh) * 2021-12-07 2022-03-18 天津大学 基于光频域反射分布式传感的假峰消除方法
CN114719955A (zh) * 2022-04-15 2022-07-08 北京理工大学 光频域反射仪分布式光纤测量中的相干衰落噪声抑制方法

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102322880A (zh) * 2011-08-18 2012-01-18 天津大学 偏振敏感的分布式光频域反射扰动传感装置和解调方法
WO2013010046A1 (en) * 2011-07-14 2013-01-17 Insight Photonic Solutions, Inc. System and method for improved resolution, higher scan speeds and reduced processing time in scans involving swept-wavelength interferometry
CN104155619A (zh) * 2014-08-20 2014-11-19 天津大学 基于磁致伸缩分布式光频域反射磁场传感装置和解调方法
US20140362367A1 (en) * 2013-06-10 2014-12-11 General Photonics Corporation Distributed fiber bend and stress measurement for determining optical fiber reliability by multi-wavelength optical reflectometry
CN106441386A (zh) * 2016-09-29 2017-02-22 西南交通大学 基于分布式光纤传感***的数据处理方法及装置
CN106872071A (zh) * 2016-12-30 2017-06-20 江苏骏龙光电科技股份有限公司 一种基于光频域反射技术的测温方法
CN107270952A (zh) * 2017-07-27 2017-10-20 天津求实飞博科技有限公司 基于光频域反射长距离光纤分布式扰动传感信号处理方法
CN110261799A (zh) * 2019-06-19 2019-09-20 中国科学院合肥物质科学研究院 利用分布式光纤传感技术的高温超导磁体失超检测***
CN110579177A (zh) * 2019-07-30 2019-12-17 天津大学 基于相对相位变化的光频域反射分布式传感解调方法
CN111854918A (zh) * 2020-07-01 2020-10-30 国网河南省电力公司焦作供电公司 一种提高分布式光纤检测空间分辨率的子带处理方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013010046A1 (en) * 2011-07-14 2013-01-17 Insight Photonic Solutions, Inc. System and method for improved resolution, higher scan speeds and reduced processing time in scans involving swept-wavelength interferometry
CN102322880A (zh) * 2011-08-18 2012-01-18 天津大学 偏振敏感的分布式光频域反射扰动传感装置和解调方法
US20140362367A1 (en) * 2013-06-10 2014-12-11 General Photonics Corporation Distributed fiber bend and stress measurement for determining optical fiber reliability by multi-wavelength optical reflectometry
CN104155619A (zh) * 2014-08-20 2014-11-19 天津大学 基于磁致伸缩分布式光频域反射磁场传感装置和解调方法
CN106441386A (zh) * 2016-09-29 2017-02-22 西南交通大学 基于分布式光纤传感***的数据处理方法及装置
CN106872071A (zh) * 2016-12-30 2017-06-20 江苏骏龙光电科技股份有限公司 一种基于光频域反射技术的测温方法
CN107270952A (zh) * 2017-07-27 2017-10-20 天津求实飞博科技有限公司 基于光频域反射长距离光纤分布式扰动传感信号处理方法
CN110261799A (zh) * 2019-06-19 2019-09-20 中国科学院合肥物质科学研究院 利用分布式光纤传感技术的高温超导磁体失超检测***
CN110579177A (zh) * 2019-07-30 2019-12-17 天津大学 基于相对相位变化的光频域反射分布式传感解调方法
CN111854918A (zh) * 2020-07-01 2020-10-30 国网河南省电力公司焦作供电公司 一种提高分布式光纤检测空间分辨率的子带处理方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
SHUAI QU.ET AL.: "Distributed sparse sinal sensing based on compressive sensing OFDR", 《OPTICS LETTERS》 *
刘浩等: "混凝土结构裂缝高空间分辨率应变测量技术探测模拟试验", 《苏州科技大学学报（工程技术版）》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114199514A (zh) * 2021-12-07 2022-03-18 天津大学 基于光频域反射分布式传感的假峰消除方法
CN114199514B (zh) * 2021-12-07 2023-05-30 天津大学 基于光频域反射分布式传感的假峰消除方法
CN114719955A (zh) * 2022-04-15 2022-07-08 北京理工大学 光频域反射仪分布式光纤测量中的相干衰落噪声抑制方法
CN114719955B (zh) * 2022-04-15 2023-04-07 北京理工大学 光频域反射仪分布式光纤测量中的相干衰落噪声抑制方法

Also Published As

Publication number Publication date
CN113237431B (zh) 2022-03-18

Similar Documents

Publication Publication Date Title
CN113237431B (zh) 一种提升ofdr***分布式空间分辨率的测量方法
CN113188461B (zh) 一种高空间分辨率下的ofdr大应变测量方法
CN102384799B (zh) 基于布里渊分布式光纤传感***相干检测方案的扫频及数据处理方法
Liu et al. Efficient common-mode noise suppression for fiber-optic interferometric sensor using heterodyne demodulation
CN103940363B (zh) 基于小波互相关技术的高精度光纤应变低频传感解调方法
CN116818085B (zh) 基于弱光栅的分布式光纤声波传感***及其测量方法
Chen et al. Denoising method of the Φ-OTDR system based on EMD-PCC
CN111780856B (zh) 基于瑞利散射谱的相位谱分析的光纤分布式振动测量方法
CN113218320B (zh) 一种基于距离域补偿的ofdr大应变测量方法
CN105066899B (zh) 一种互参考的光纤激光静态应变传感解调***
CN109556754A (zh) 光纤应变和温度测定装置以及光纤应变和温度测定方法
Bai et al. Noise reduction method of Φ-OTDR system based on EMD-TFPF algorithm
CN115165071A (zh) 一种基于机器学习与相位分段差分实现振动精确定位的方法和***
CN110806259A (zh) 一种用于光纤传感高频扰动定位与检测的装置
CN112880711B (zh) 一种基于双脉冲调制的分布式光纤传感方法及***
CN111579048B (zh) 一种基于滑动时间选通的ofdr***振动检测方法
CN105115528B (zh) 基于参考光纤激光器的高精度静态应变拍频解调***
CN112229501A (zh) 一种基于ofdr的振动信号自动检测装置和方法
CN115790814B (zh) 光纤振动检测***及其方法
CN110987151A (zh) 一种通讯光缆状态实时监测***
Pu et al. Acceleration of the frequency-shift demodulation in phase-sensitive OTDR
Rossi et al. Low-noise φ-OTDR employing nonlinear optical preamplification for distributed acoustic sensing
Li et al. Spatio-temporal joint oversampling-downsampling technique for ultra-high resolution fiber optic distributed acoustic sensing
Hou et al. Improved path imbalance measurement of a fiber-optic interferometer based on frequency scanning interferometry
Kiyozumi et al. Modified expression for spatial resolution in optical correlation-domain reflectometry

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