CN114674352B - 基于瑞利散射光谱非相似性的分布式扰动传感和解调方法 - Google Patents
基于瑞利散射光谱非相似性的分布式扰动传感和解调方法 Download PDFInfo
- Publication number
- CN114674352B CN114674352B CN202210334748.XA CN202210334748A CN114674352B CN 114674352 B CN114674352 B CN 114674352B CN 202210334748 A CN202210334748 A CN 202210334748A CN 114674352 B CN114674352 B CN 114674352B
- Authority
- CN
- China
- Prior art keywords
- disturbance
- curve
- lnred
- domain
- rayleigh scattering
- 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
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000001228 spectrum Methods 0.000 title claims abstract description 37
- 239000013307 optical fiber Substances 0.000 claims abstract description 66
- 230000003287 optical effect Effects 0.000 claims abstract description 33
- 238000005259 measurement Methods 0.000 claims abstract description 19
- 238000001914 filtration Methods 0.000 claims abstract description 17
- 238000006243 chemical reaction Methods 0.000 claims description 12
- 230000008859 change Effects 0.000 claims description 10
- 230000009467 reduction Effects 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 7
- 238000010606 normalization Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 4
- 230000035559 beat frequency Effects 0.000 abstract description 5
- 230000010287 polarization Effects 0.000 description 24
- 239000000835 fiber Substances 0.000 description 12
- 238000005516 engineering process Methods 0.000 description 8
- 238000002168 optical frequency-domain reflectometry Methods 0.000 description 8
- 230000001427 coherent effect Effects 0.000 description 7
- 238000012360 testing method Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000001919 Rayleigh scattering spectroscopy Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D5/00—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable
- G01D5/26—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light
- G01D5/32—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light
- G01D5/34—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells
- G01D5/353—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre
- G01D5/35306—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre using an interferometer arrangement
- G01D5/35332—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre using an interferometer arrangement using other interferometers
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Light Guides In General And Applications Therefor (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明涉及一种基于瑞利散射光谱非相似性的分布式扰动传感和解调方法，包括下列步骤：在待测光纤上没有施加扰动和施加扰动时，分别利用光频域反射系统测量得到一组波长域拍频信号；利用去斜滤波进行非线性相位补偿得到参考光谱和测量光谱；利用快速傅里叶变换从波长域转换到距离域；利用滑动窗进行分段；利用快速逆傅里叶变换获得局部参考瑞利散射光谱和局部测量瑞利散射光谱；利用小波去噪降噪处理；计算非相似度；得到相对欧氏距离随滑动窗位置的归一化曲线；低通滤波；定位扰动点。
Description
技术领域
本发明属于分布式光纤传感仪器技术领域，具体涉及一种基于瑞利散射光谱非相似性的分布式扰动传 感和解调方法。
背景技术
分布式的光纤扰动传感方法在诸多领域有着广泛的应用，比如在入侵探测、油气管道监测、通信光缆 的监测等。光频域反射(Optical Frequency Domain Reflectometry,OFDR)技术，作为分布式光纤传感的一种， 具有抗电磁干扰、强适应性、定位精度高、灵敏度高等优点，其在扰动传感方面有着很大的优势。基于瑞 利散射光谱分析的OFDR技术以及被广泛用于短距离扰动的测量中，并且具有很高的测量灵敏度和空间分 辨率。但是在将其应用与长距离扰动的测量中，由于光纤中背向瑞利散射本身就很微弱，再加上环境随机 噪声的影响也逐渐增大，对瑞利散射光谱信号会产生很大影响，使得在长距离下进行分布式的扰动传感并 要保证较高的空间分辨率就变得困难。在现有基于OFDR技术实现长距离光纤扰动分布式传感方法中，刘 铁根等人[1]将不同时刻的外界扰动的瑞利散射光谱进行比较，根据互相关系数和相关图的噪声水平判断传 感光纤上是否存在扰动，该方法的测试距离可达120km-150km，但其并未提及空间分辨率的情况；杜阳等 人[2]在光频域上利用互相关分析的方法对信号扰动状态下的相似水平作评价，可以实现92.03km光纤末端 双点扰动的定位，其空间分辨率只能达到26.7m；丁振扬等人[3]通过在光频域计算瑞利散射光谱的频移定 位第一个扰动点，在距离域寻找瑞利散射信号的凹陷波谷定位第二个扰动点，该方法需要同时在波长域和 距离域作处理，定位算法复杂，且空间分辨率只能达到13.3m。
光纤中的瑞利散射是由于光纤中随机波动的折射率产生的。对于一段光纤，其瑞利散射的幅度或折射 率是随距离随机变化的，但这种变化又是相对静态的特征信息，瑞利散射可以理解为光纤的“指纹”。通过 OFDR的拍频干涉信号得到这种空间分布的“指纹”，然后通过这些指纹信息发生的变化就可以得到沿光纤 上的扰动信息。当在光纤上某点处产生扰动时，将会导致该点处光纤折射率发生变化，由于光弹效应，此 处产生的背向瑞利散射光的相位也会发生变化，从而引起背向瑞利散射光强度发生变化，在瑞利散射光谱 上会出现频移。因此，无扰动的参考瑞利散射光谱和有扰动的测量瑞利散射光谱在扰动点处的相似度会很 低，即非相似度很高。通过计算参考瑞利散射光谱和测量瑞利散射光谱的非相似度得到非相似度变化大的 点就可以定位为扰动点。当光纤上存在两个及以上扰动点时，后扰动点处的背向瑞利散射同样会受到前扰 动点处的扰动影响，因此后扰动点处的非相似度会在前扰动点处的非相似度基础上出现一个阶跃的变化， 非相似度曲线在后扰动点会出现幅度的叠加。因此在理论上，通过非相似度变化可以定位两点及以上的扰 动。
由于光纤中的瑞利散射非常微弱，以至于非常容易被一些随机噪声淹没掉，因此在长距离光纤分布式 扰动传感中，随机噪声对于扰动点定位的影响是非常大的，利用一些去噪算法对其进行处理是非常有必要 的。同时如果取光纤某一段“指纹”信息，其包含数个甚至上百个数据点，这样也可以降低随机噪声对扰动 点定位的影响。
参考文献：
[1]丁振扬，刘铁根，刘琨，江俊峰，杨迪，一种长距离光频域反射光纤分布式多点扰动传感方法，2019.5.28， 中国，ZL201710185731.1
[2]杜阳.基于OFDR的分布式光纤多参量传感方法研究[D].天津大学,2016.
[3]刘铁根，刘琨，丁振扬，江俊峰，杜阳，李定杰，一种基于瑞利散射光谱相关系数的分布式扰动传感 装置和解调方法，2014.9，中国，ZL 201210100999.8
[4]丁振扬.几种改进OFDR性能方法的提出及验证[D].天津大学,2013.
发明内容
本发明提供了一种基于瑞利散射光谱非相似性的分布式扰动传感和解调方法，本发明克服了现有基于 光频域反射的长距离光纤扰动传感方法中，对多点扰动位置无法区分、传感信号噪声较大、空间分辨率不 高以及定位算法复杂等问题，尽可能针对扰动对传感信号造成的本质影响提出测量多点扰动的新型算法， 进而实现多点扰动的区分、传感信号噪声的降低、空间分辨率的提高以及定位算法的简化，该方法的步骤 如下：
一种基于瑞利散射光谱非相似性的分布式扰动传感和解调方法，包括下列步骤：
第一步，在待测光纤上没有施加扰动时，利用光频域反射系统测量得到一组波长域拍频信号，作为参 考信号Sr；施加扰动后，利用光频域反射系统测量得到一组波长域拍频信号，作为测量信号Sm。
第二步，利用去斜滤波对Sr进行非线性相位补偿得到线性拍频信号RS，即为参考光谱；利用去斜滤 波对Sm进行非线性相位补偿得到线性拍频信号MS，即为测量光谱；
第三步，利用快速傅里叶变换将RS从波长域转换到距离域，得到距离域信号Rr；利用快速傅里叶变 换将MS从波长域转换到距离域，得到距离域信号Rm；
第四步，利用大小为N的滑动窗对Rr进行分段得到Rri；利用滑动窗对Rm进行分段得到Rmi，其中 i∈[1,M]，M为段数，N和M之间满足M＝采样点数/N；
第五步，利用快速逆傅里叶变换将Rri反变换回波长域，转换为局部参考瑞利散射光谱，记为 LRSi＝(lrs1,lrs2,…,lrsN)；利用快速逆傅里叶变换将Rmi反变换回波长域，转换为局部测量瑞利散射光 谱，记为LMSi＝(lms1,lms2,…,lmsN)；
第六步，利用小波去噪对LRSi进行降噪得到DLRSi＝(dlrs1,dlrs2,…,dlrsN)；
利用小波去噪对LMSi进行降噪得到DLMSi＝(dlms1,dlms2,…,dlmsN)；
第七步，利用相对欧氏距离衡量DLRSi和DLMSi的非相似度，相对欧氏距离公式为
相对欧氏距离值越大说明DLRSi和DLMSi之间的非相似度越大，相似度越小；
第八步，M段信号均完成第五-七步处理后，得到相对欧氏距离随滑动窗位置的变化曲线，记为曲线A， 横坐标为窗i，纵坐标为相对欧式距离(RED)；再对曲线A进行归一化，得到相对欧氏距离归一化曲线， 记为曲线B，横坐标为窗i，纵坐标为归一化相对欧氏距离(NRED)；
第九步，对曲线B进行低通滤波，滤除毛刺，得到相对欧氏距离归一化曲线的骨架，记为曲线C，横 坐标为窗i，纵坐标为低通滤波后归一化相对欧氏距离(LNRED)，C曲线上数据点为LNRED(i)；
第十步，在曲线C上定位扰动点；曲线C上判断扰动点的条件有如下三点：(1)选取LNRED(i)值不 小于阈值a1即满足LNRED(i)≥a1的所有数据点，阈值a1的选取需要避免出现误报的情况；(2)将一段连 续的满足前后点LNRED(i)差值不小于阈值a2即满足LNRED(i)-LNRED(i-1)≥a2的数据点归为一 组，分组进行处理，阈值a2的选取需要避免出现漏报和误报的情况；(3)在某一组数据点中，将该组数据第 一点横坐标记为x，该组数据最后一点横坐标记为y，若同时满足LNRED(y)-LNRED(x)≥a3， LNRED(y)-max[LNRED(i)]≥a4,i＜x两个条件，则x点即可确定为扰动点；阈值a3、a4的选取需要 避免出现漏报和误报的情况；
第十一步，将曲线C横坐标i转换为距离z即可实现扰动点的定位，距离z满足
其中，L为待测光纤总长度。
本专利提出一种基于光频域反射技术实现长距离光纤多点扰动传感方法。通过去斜滤波算法补偿光源 非线性调谐效应，通过小波去噪降低了随机噪声对长距离光纤扰动定位的影响，根据光纤中扰动点背向瑞 利散射光强度变化，在瑞利散射光谱上会出现频移，参考光谱和测量光谱在扰动点处非相似度变化大的特 征，区分多点扰动。
本发明提供的技术方案的有益效果是：
1、实现了光频域反射中光源非线性调谐效应的补偿。
2、降低了随机噪声对长距离光纤扰动定位的影响。
3、实现了长距离光纤多点扰动定位。
4、大幅提升了基于光频域反射技术定位长距离光纤中多点扰动的空间分辨率。
附图说明
图1为基于光频域反射技术实现92.03km光纤多点扰动传感方法的流程图；
图2为运用光频域反射系统测量92.03km长度光纤扰动的传感采集信息装置示意图；
图3为背向瑞利散射距离域信息(a)参考信号(b)测量信号；
图4为待测光纤末端的距离域信息差异(a)参考信号(b)测量信号；
图5为RED随滑动窗位置的变化曲线；
图6为多点扰动定位图(a)整体图(b)局部放大图。
附图中，各标号所代表的部件列表如下：
1：可调谐激光器； 2：1:99保偏光分束器；
3：第一环形器； 4：第一50:50耦合器；
5：延迟光纤； 6：第一法拉第转镜；
7：第二法拉第转镜； 8：平衡探测器；
9：50:50保偏耦合器； 10：参考臂；
11：测试臂； 12：第二环形器；
13：第一扰动源； 14：第二扰动源；
15：待测光纤； 16：偏振分集相干接收模块；
17：采集卡； 18：计算机；
19：附加干涉仪； 20：主干涉仪。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面对本发明实施方式做进一步详细描述。
本发明是通过小波去噪降低随机噪声对长距离光纤扰动定位的影响，根据光纤中发生扰动处背向瑞利 散射光强发生变化，以待测光纤同一位置在无扰动和有扰动两种情况下得到的局部瑞利散射光谱之间的非 相似度变化大为依据，在整段光纤的非相似度曲线中找寻非相似度变化大的点，从而定位多点扰动。
N为滑动窗取到的局部瑞利散射光谱的传感数据点，其值等于滑动窗的宽度，ΔZ是每个数据点的空 间分辨率，计算公式为
ΔZ＝c/2nΔF (2)
其中，c为真空中光速，n为光纤的折射率，ΔF为激光器的调谐范围，在本发明中，ΔF取为3GHz。
而空间分辨率计算公式为
ΔX＝N·ΔZ (3)
通过本发明实施例提出的基于光频域反射技术实现92.03km光纤多点扰动传感方法可以大幅减小滑动窗的 宽度，将其由800个点减小至150个点，空间分辨率由26.7米提升至4.95米。
实施例一：
本实例包括基于光频域反射系统分布式光纤测量装置
基于光频域反射分布式光纤测量系统，包括：可调谐激光器1、1:99光分束器2、采集卡17、计算机 18、附加干涉仪19和主干涉仪20。
附加干涉仪19为偏振不敏感迈克尔逊干涉结构，用于获得光源的非线性相位信息，包括：第一环形 器3、第一50:50耦合器4、延迟光纤5、第一法拉第转镜6、第二法拉第转镜7和平衡探测器8。
主干涉仪20为全保偏结构，包括：50:50保偏耦合器9、参考臂10、测试臂11、第二环形器12、第 一扰动源13、第二扰动源14、待测光纤15和偏振分集相干接收模块16。
可调谐激光器1与1:99光分束器2的a端口相连；1:99光分束器2的b端口即1％分光口与第一环形 器3的a端口相连；1:99光分束器2的c端口即99％分光口与50:50保偏耦合器9的a端口相连；第一环 形器3的b端口与第一50:50耦合器4的a端口相连；第一环形器3的c端口与平衡探测器8的输入端相 连；第一50:50耦合器4的b端口与平衡探测器8的输入端相连；第一50:50耦合器4的c端口通过延迟 光纤5与第一法拉第转镜6相连；第一50:50耦合器4的d端口与第二法拉第转镜7相连；平衡探测器8 的输出端与采集卡17的某一通道相连；50:50保偏耦合器9的c端口即50％分端口通过参考臂10与偏振 分集相干接收模块16的a端口相连；50:50保偏耦合器9的d端口即50％分端口通过测试臂11与第二环 形器12的a端口相连；第二环形器12的b端口与待测光纤15相连；第二环形器12的c端口与偏振分集 相干接收模块16的b端口相连；偏振分集相干接收模块16的c端口与采集卡17的某一通道相连；采集 卡17的输出端与计算机18的输入端相连。
装置工作时，可调谐激光器1的出射光由1:99保偏光分束器2的a端口进入，从1:99保偏光分束器2 的b端口即1％分光口进入环形器3的a端口，经环形器3的b端口，光从第一50:50耦合器4的a端口进 入，50％从c端口出射、50％从d端口出射，来自第一50:50耦合器4的c端口的出射光经延迟光纤5被第 一法拉第转镜6反射，返回第一50:50耦合器4的c端口,来自第一50:50耦合器4的d端口的出射光被第 二法拉第转镜7反射，返回第一50:50耦合器4的d端口，两束光在第一50:50耦合器4中发生干涉，从a、 b端口输出，第一50:50耦合器4的a端口的光由第一环形器3的b端口输入，c端口输出，该输出光与第 一50:50耦合器4的b端口的输出光在平衡探测器8中进行外差检测转换为模拟电信号。从1:99保偏光分 束器2的c端口即99％分光口输出光由50:50保偏耦合器9的a端口输入，c端口50％输出光进入参考臂10，d端口50％输出光进入测试臂11，测试臂11中的光由第二环形器12的a端口输入，从第二环形器12 的b端口进入待测光纤15，而待测光纤的背向瑞利散射光从第二环形器12的b端口进入，c端口输出， 该光与50:50保偏耦合器9的c端口的输出光在偏振分集相干接收模块16中发生干涉，并转换为模拟电信 号。采集卡17将采集到的模拟电信号转换为数字电信号传输至计算机18。
可调谐激光器1用于为光频域反射系统提供光源，其光频能够进行线性扫描。
第一环形器3防止附加干涉仪中50:50耦合器4的b端口反射光进入激光器。
50:50耦合器4用于分光与光干涉。
第一法拉第转镜6和第二法拉第转镜7用于为附加干涉仪提供反射，且能够消除附加干涉仪的偏振衰 落现象。
延迟光纤5用于实现非等臂的拍频干涉，能够根据拍频和延迟光纤长度得到光频。
平衡探测器8用于采集50:50耦合器4从b端口的出射光及第一环形器3从c端口的出射光，即附加 干涉仪的时域干涉拍频信号。
偏振分集相干接收模块16含偏振分集和平衡探测模块，将参考光与信号光的两个偏振态进行分别相 干，采用两路高速低噪声平衡探测器进行分别接收，消除偏振衰落问题。
采集卡17：将干涉仪产生的波长域拍频信号经模数转换后的数字信号通过USB、PCIE等总线方式传 输到计算机。
计算机18：对采集卡模数转换后的干涉信号进行数据处理，实现光纤的分布式测量。
实施例二：
本实例提供了一种基于光频域反射技术实现92.03km光纤多点扰动传感的新型算法，步骤如下：
第一步，在92.03km待测光纤上没有施加扰动时，在光频域反射仪的主干涉仪中由光纤背向瑞利散射 形成一组波长域拍频干涉信号，作为参考信号Sr；在92.03km待测光纤的90.21km和91.39km处施加两 个扰动源，在光频域反射仪的主干涉仪中由光纤背向瑞利散射形成一组波长域拍频干涉信号，作为测量信 号Sm。
第二步，利用去斜滤波算法对Sr进行非线性相位的补偿得到参考光谱RS；利用去斜滤波算法对Sm进行非线性相位的补偿得到测量光谱MS。去斜滤波算法是一种滤波算法，可以用于进行非线性相位的补 偿，本实施例选用的是文献[2]提出的优化去斜滤波算法，此算法是在文献[4]提出的去斜滤波算法的基础上 的进一步优化。
第三步，利用快速傅里叶变换将RS从波长域转换到距离域，得到距离域信号Rr，如图3(a)所示；利 用快速傅里叶变换将MS从波长域转换到距离域，得到距离域信号Rm，如图3(b)所示。Rr和Rm两组距 离域信号在光纤末端的差异如图4所示。
第四步，利用大小为N的滑动窗对Rr进行分段得到Rri；利用滑动窗对Rm进行分段得到Rmi，其中 i∈[1,M]，M为段数，N和M之间满足M＝采样点数/N。本专利中采样点数＝2761000，N＝150， 则M＝18406。N值的选取需要平衡空间分辨率和误报的情况。
第五步，利用快速逆傅里叶变换将Rri反变换回波长域，转换为局部参考瑞利散射光谱，记为 LRSi＝(lrs1,lrs2,…,lrsN)；利用快速逆傅里叶变换将Rmi反变换回波长域，转换为局部测量瑞利散射光 谱，记为LMSi＝(lms1,lms2,…,lmsN)。
第六步，利用小波去噪对LRSi进行降噪得到DLRSi＝(dlrs1,dlrs2,…,dlrsN)；利用小波去噪对 LMSi进行降噪得到DLMSi＝(dlms1,dlms2,…,dlmsN)。小波函数和分解层数的选取需要同时兼顾降噪 效果以及信号降噪后的完整性。本实施例的小波基函数选为“db2”，小波分解层数选为2，提取第一层的细 节系数来估算噪声的标准偏差确定阈值，选取硬阈值函数，每层都采用同一个阈值进行处理。
第七步，利用相对欧氏距离公式(1)衡量DLRSi和DLMSi的非相似度。
第八步，M段信号均完成第五-七步处理后，得到相对欧氏距离随滑动窗位置的变化曲线，记为曲线A， 横坐标为窗i，纵坐标为相对欧氏距离(RED)，如图5所示。再将曲线A进行归一化，得到相对欧氏距离 归一化曲线，记为曲线B，横坐标为窗i，纵坐标为归一化相对欧氏距离(NRED)。
第九步，对曲线B进行低通滤波，滤除高频干扰，得到相对欧氏距离归一化曲线的骨架，记为曲线C， 横坐标为窗i，纵坐标为低通滤波后归一化相对欧氏距离(LNRED)，C曲线上数据点为LNRED(i)。低 通滤波器的采样频率保持和采集卡的采样频率相同，低截止频率依据所需信号频带选取，要在确保将毛刺 滤除的同时保留信号变化的整体趋势。本实施例低通滤波器的采样频率设为100M，低截止频率设为0.5M。
第十步，在曲线C上定位扰动点。曲线C上判断扰动点的条件有如下三点：(1)由于在上述第八步中 做了归一化处理，放大了数据之间的差异，首先需要选取LNRED(i)值不小于阈值a1＝0.4即满足 LNRED(i)≥0.4的所有数据点。(2)将一段连续的满足前后点LNRED(i)差值不小于阈值a2＝0.00011 即满足LNRED(i)-LNRED(i-1)≥0.00011的数据点归为一组，分组进行处理。(3)在某一组数据点中， 将该组数据第一点的横坐标记为x，该组数据最后一点的横坐标记为y，如果同时满足 LNRED(y)-LNRED(x)≥a3＝0.04，LNRED(y)-max[LNRED(i)]≥a4＝0.01,i＜x两个条件，则 x点即可确定为扰动点。
上述多个阈值的选取，需要根据具体应用场景和实验进行标定。阈值a1的选取需要避免出现误报的情况。 阈值a2的选取需要避免出现漏报和误报的情况。阈值a3、a4的选取需要避免出现漏报和误报的情况。
第十一步，将曲线C横坐标i转换为距离z即可实现扰动点的定位，如图6所示。本专利中待测光纤 总长度为92.03km。
实施例三
下面结合具体的试验对实施例1-2中的光纤分布式测量系统和测量方法进行可行性验证，详见下文描 述：
本发明实施例验证实验所采用待测光纤15为50km和42.03km单模光纤经法兰盘相接，测量目的为待 测光纤末端两点扰动的定位。
在实验中，我们使用OFDR系统，可调谐激光器的线宽约为1kHz，中心波长为1550nm，光源的输出 功率约为10Mw，可调谐激光器的调谐速率为40GHz/s，调谐范围为3GHz。整个扫描时间为0.075s。采集 卡的采样率为100MS/s。
两组距离域信号Rr和Rm末端差异如图4所示，参考组中光纤末端存在一个明显的反射峰，而测量组 中光纤末端由于扰动的影响导致反射峰能量发生了扩散，本专利正是利用这一特征通过计算两组信号的非 相似度，从而实现扰动的定位。
图6为利用本专利所提出方法在光纤末端定位的两个扰动点。图中曲线首端LNRED迅速上升是由于 在此处的光纤长度还很短，随机噪声对相似度的影响还很小，非相似度保持在一个较低的水平，而在光纤 末端，非相似度保持在一个较高的水平，如图5所示，归一化处理会放大数据之间的差异，从而导致图6 中曲线首端出现骤升的现象，并非是由扰动引起，因此本专利在判断LNRED变化大的条件中加入了“选 取LNRED大于0.4的所有数据点”，目的就是为了排除曲线首端幅值骤升的影响。而图6曲线中部出现的 LNRED骤降是由于两组信号在50km处的法兰盘都存在反射峰，从而两组信号在该位置处的相似度变大。
该分布式测量系统和测量方法实现了光纤分布式测量，以4.95米的空间分辨率实现了92.03km光纤 末端多点扰动的定位。
综上所述，本专利实现了基于光频域反射技术对92.03km光纤末端多点扰动的定位，极大地提高了基 于长距离光频域反射光纤分布式扰动传感的空间分辨率。
本发明实施例对各器件的型号除做特殊说明的以外，其他器件的型号不做限制，只要能完成上述功能 的器件均可。
本领域技术人员可以理解附图只是一个优选实施例的示意图，上述本发明实施例序号仅仅为了描述， 不代表实施例的优劣。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任 何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (1)
1.一种基于瑞利散射光谱非相似性的分布式扰动传感和解调方法，包括下列步骤：
第一步，在待测光纤上没有施加扰动时，利用光频域反射系统测量得到一组波长域拍频信号，作为参考信号Sr；施加扰动后，利用光频域反射系统测量得到一组波长域拍频信号，作为测量信号Sm；
第二步，利用去斜滤波对Sr进行非线性相位补偿得到线性拍频信号RS，即为参考光谱；利用去斜滤波对Sm进行非线性相位补偿得到线性拍频信号MS，即为测量光谱；
第三步，利用快速傅里叶变换将RS从波长域转换到距离域，得到距离域信号Rr；利用快速傅里叶变换将MS从波长域转换到距离域，得到距离域信号Rm；
第四步，利用大小为N的滑动窗对Rr进行分段得到Rri；利用滑动窗对Rm进行分段得到Rmi，其中i∈[1,M]，M为段数，N和M之间满足M＝采样点数/N；
第五步，利用快速逆傅里叶变换将Rri反变换回波长域，转换为局部参考瑞利散射光谱，记为LRSi＝(lrs1,lrs2,…,lrsN)；利用快速逆傅里叶变换将Rmi反变换回波长域，转换为局部测量瑞利散射光谱，记为LMSi＝(lms1,lms2,…,lmsN)；
第六步，利用小波去噪对LRSi进行降噪得到DLRSi＝(dlrs1,dlrs2,…,dlrsN)；
利用小波去噪对LMSi进行降噪得到DLMSi＝(dlms1,dlms2,…,dlmsN)；
第七步，利用相对欧氏距离衡量DLRSi和DLMSi的非相似度，相对欧氏距离公式为
相对欧氏距离值越大说明DLRSi和DLMSi之间的非相似度越大，相似度越小；
第八步，M段信号均完成第五-七步处理后，得到相对欧氏距离随滑动窗位置的变化曲线，记为曲线A，横坐标为窗i，纵坐标为相对欧式距离(RED)；再对曲线A进行归一化，得到相对欧氏距离归一化曲线，记为曲线B，横坐标为窗i，纵坐标为归一化相对欧氏距离(NRED)；
第九步，对曲线B进行低通滤波，滤除毛刺，得到相对欧氏距离归一化曲线的骨架，记为曲线C，横坐标为窗i，纵坐标为低通滤波后归一化相对欧氏距离(LNRED)，C曲线上数据点为LNRED(i)；
第十步，在曲线C上定位扰动点；曲线C上判断扰动点的条件有如下三点：(1)选取LNRED(i)值不小于阈值a1即满足LNRED(i)≥a1的所有数据点，阈值a1的选取需要避免出现误报的情况；(2)将一段连续的满足前后点LNRED(i)差值不小于阈值a2即满足LNRED(i)-LNRED(i-1)≥a2的数据点归为一组，分组进行处理，阈值a2的选取需要避免出现漏报和误报的情况；(3)在某一组数据点中，将该组数据第一点横坐标记为x，该组数据最后一点横坐标记为y，若同时满足LNRED(y)-LNRED(x)≥a3，LNRED(y)-max[LNRED(i)]≥a4，i＜x两个条件，则x点即可确定为扰动点；阈值a3、a4的选取需要避免出现漏报和误报的情况；
第十一步，将曲线C横坐标i转换为距离z即可实现扰动点的定位，距离z满足
其中，L为待测光纤总长度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210334748.XA CN114674352B (zh) | 2022-03-31 | 2022-03-31 | 基于瑞利散射光谱非相似性的分布式扰动传感和解调方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210334748.XA CN114674352B (zh) | 2022-03-31 | 2022-03-31 | 基于瑞利散射光谱非相似性的分布式扰动传感和解调方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114674352A CN114674352A (zh) | 2022-06-28 |
| CN114674352B true CN114674352B (zh) | 2023-09-15 |
Family
ID=82076538
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210334748.XA Active CN114674352B (zh) | 2022-03-31 | 2022-03-31 | 基于瑞利散射光谱非相似性的分布式扰动传感和解调方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114674352B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115790815B (zh) * | 2023-01-17 | 2023-05-16 | 常熟理工学院 | 一种分布式光纤传感系统的扰动快速鉴别方法及系统 |
| CN116952150A (zh) * | 2023-07-04 | 2023-10-27 | 哈尔滨工业大学 | 应变值测量方法、ofdr应变测量方法 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2015165394A1 (zh) * | 2013-05-31 | 2015-11-05 | 欧普图斯(苏州)光学纳米科技有限公司 | 面向多行业检测的激光拉曼光谱智能化辨识方法及系统 |
| JP2016053525A (ja) * | 2014-09-03 | 2016-04-14 | 日本電信電話株式会社 | 光ファイバの温度・歪み分布測定方法および装置 |
| CN107657242A (zh) * | 2017-10-10 | 2018-02-02 | 湖南师范大学 | 一种大地电磁信噪辨识及分离方法 |
| CN108875684A (zh) * | 2018-06-29 | 2018-11-23 | 电子科技大学 | 基于光纤感测时频图处理的列车运行状态参数估计方法 |
| CN109100788A (zh) * | 2018-07-06 | 2018-12-28 | 东北石油大学 | 地震数据非局部均值去噪方法 |
| CN110763913A (zh) * | 2019-10-14 | 2020-02-07 | 南京信息工程大学 | 一种基于信号分段分类的导数谱平滑处理方法 |
| CN112422546A (zh) * | 2020-11-10 | 2021-02-26 | 昆明理工大学 | 一种基于变邻域算法和模糊聚类的网络异常检测方法 |
-
2022
- 2022-03-31 CN CN202210334748.XA patent/CN114674352B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2015165394A1 (zh) * | 2013-05-31 | 2015-11-05 | 欧普图斯(苏州)光学纳米科技有限公司 | 面向多行业检测的激光拉曼光谱智能化辨识方法及系统 |
| JP2016053525A (ja) * | 2014-09-03 | 2016-04-14 | 日本電信電話株式会社 | 光ファイバの温度・歪み分布測定方法および装置 |
| CN107657242A (zh) * | 2017-10-10 | 2018-02-02 | 湖南师范大学 | 一种大地电磁信噪辨识及分离方法 |
| CN108875684A (zh) * | 2018-06-29 | 2018-11-23 | 电子科技大学 | 基于光纤感测时频图处理的列车运行状态参数估计方法 |
| CN109100788A (zh) * | 2018-07-06 | 2018-12-28 | 东北石油大学 | 地震数据非局部均值去噪方法 |
| CN110763913A (zh) * | 2019-10-14 | 2020-02-07 | 南京信息工程大学 | 一种基于信号分段分类的导数谱平滑处理方法 |
| CN112422546A (zh) * | 2020-11-10 | 2021-02-26 | 昆明理工大学 | 一种基于变邻域算法和模糊聚类的网络异常检测方法 |
Non-Patent Citations (3)
| Title |
|---|
| 基于OFDR的光纤瑞利散射多特征信息分布式传感研究;杨迪;信息科技辑;全文 * |
| 基于光频域反射技术的光纤连续分布式定位应变传感;刘琨;冯博文;刘铁根;江俊峰;杜阳;;中国激光(05);全文 * |
| 基于小波系数的近红外光谱局部建模方法与应用研究;石雪;蔡文生;邵学广;分析化学(008);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114674352A (zh) | 2022-06-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114674352B (zh) | 基于瑞利散射光谱非相似性的分布式扰动传感和解调方法 | |
| CN106679790B (zh) | 一种提高分布式光纤振动传感灵敏度的互相关解调方法 | |
| CN107402082B (zh) | 一种布里渊散射信号处理方法及其分布式光纤传感系统 | |
| CN108955734B (zh) | 一种光纤f-p温度/压力复合传感器的腔长解调方法 | |
| CN102879081A (zh) | 一种分布式光纤振动系统中的数据处理方法 | |
| WO2019019735A1 (zh) | 基于光频域反射长距离光纤分布式扰动传感信号处理方法 | |
| CN113124994B (zh) | 一种基于弱信噪比信号的分布式振动传感解调方法 | |
| CN103499383A (zh) | 一种提高光纤振动传感器定位精度的自相关定位方法 | |
| CN104077267A (zh) | 一种双重傅立叶变换扰动定位方法 | |
| Ma et al. | A study on noise reduction of Φ-OTDR system based on VSS-NLMS algorithm | |
| Bai et al. | Noise reduction method of Φ-OTDR system based on EMD-TFPF algorithm | |
| Ye et al. | A positioning algorithm realized multilateration for distributed fiber‐optic sensor | |
| CN107014409B (zh) | 一种长距离光频域反射光纤分布式多点扰动传感方法 | |
| CN111157102B (zh) | 一种分布式光纤传感系统中消除频率干扰的定位方法 | |
| CN112229501A (zh) | 一种基于ofdr的振动信号自动检测装置和方法 | |
| Li et al. | SNR enhancement of far-end disturbances on distributed sensor based on phase-sensitive optical time-domain reflectometry | |
| CN112179475B (zh) | 一种多源混叠的分布式光纤振动传感信号的分离方法 | |
| CN115752696A (zh) | 基于前向传输的分布式光纤振动传感系统及信号处理方法 | |
| CN112697303B (zh) | 面向智能电网的分布式光纤传感系统和检测方法 | |
| CN114861777A (zh) | 基于光频率反射技术的局部谱匹配方法、装置和存储介质 | |
| CN112683418A (zh) | 一种光纤分布式测温的拉曼散射光双路解调方法 | |
| CN113155267A (zh) | 一种基于二次相关的ofdr系统振动检测方法、系统、存储介质和终端 | |
| Rossi et al. | Low-noise φ-OTDR employing nonlinear optical preamplification for distributed acoustic sensing | |
| CN113607075A (zh) | 一种基于自适应时延估计的光频域反射应变解调方法 | |
| CN112268608A (zh) | 一种强干扰dvs的高精度解调算法 |
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 |