CN114046900A - 一种分布式光纤测温系统自适应滤波方法及装置 - Google Patents
一种分布式光纤测温系统自适应滤波方法及装置 Download PDFInfo
- Publication number
- CN114046900A CN114046900A CN202111333819.6A CN202111333819A CN114046900A CN 114046900 A CN114046900 A CN 114046900A CN 202111333819 A CN202111333819 A CN 202111333819A CN 114046900 A CN114046900 A CN 114046900A
- Authority
- CN
- China
- Prior art keywords
- filtering
- optical fiber
- temperature
- light intensity
- counter
- 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.)
- Pending
Links
- 238000001914 filtration Methods 0.000 title claims abstract description 173
- 239000013307 optical fiber Substances 0.000 title claims abstract description 103
- 238000009529 body temperature measurement Methods 0.000 title claims abstract description 59
- 238000000034 method Methods 0.000 title claims abstract description 57
- 238000005070 sampling Methods 0.000 claims abstract description 47
- 238000012545 processing Methods 0.000 claims abstract description 41
- 230000003044 adaptive effect Effects 0.000 claims abstract description 39
- 238000004364 calculation method Methods 0.000 claims abstract description 39
- 238000009825 accumulation Methods 0.000 claims abstract description 24
- 230000008859 change Effects 0.000 claims description 76
- 238000001069 Raman spectroscopy Methods 0.000 claims description 6
- 230000001186 cumulative effect Effects 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 4
- 230000035945 sensitivity Effects 0.000 description 6
- 230000000694 effects Effects 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000000835 fiber Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000000253 optical time-domain reflectometry Methods 0.000 description 2
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000010292 electrical insulation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
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
- G01K11/324—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 using Raman scattering
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Temperature Or Quantity Of Heat (AREA)
Abstract
一种分布式光纤测温系统自适应滤波方法及装置,该方法包括：对斯托克斯光强、反斯托克斯光强原始采样值进行第一级累加平均滤波；计算测温光纤的原始温度计算值；针对测温光纤指定空间位置连续时间域上的原始温度计算值进行第二级自适应低通滤波处理，更新滤波参数；针对测温光纤不同空间位置上的原始温度计算值进行第三级自适应低通滤波处理，更新滤波参数。该装置用于执行该方法。本发明自适应滤波参数的引入，比传统滤波方法更快速地跟踪测量温度的变化趋势，同时实现良好的噪声处理，减少分布式光纤测温系统中具有零均随机噪声造成的干扰。
Description
技术领域
本发明涉及分布式系统工业控制技术领域，尤其涉及一种分布式光纤测温系统自适应滤波方法及装置。
背景技术
分布式光纤测温系统采用光时域反射(OTDR)原理和拉曼(Raman)散射效应对温度的敏感实现温度监测，具备抗电磁干扰、电绝缘性好、耐腐蚀、灵敏度高以及可以实现大范围的周界监测等优点，广泛应用于电力电缆、管廊/电力遂道、电缆沟、开关柜等场景，实现电缆及电缆接头温度、通道环境温度、开关柜重要部件温度实时监控。当测温光纤所处环境温度越高，其后向拉曼散射光的强度也越大，其中反斯托克斯(anti-Stokes)光较斯托克斯(Stokes)光对温度变化更为敏感。利用这个现象，通过计算反斯托克斯(anti-Stokes)光和斯托克斯(Stokes)光的强度比值，就可以计算出反射点的温度，这就是分布式光纤测温系统的基本原理。
拉曼散射为激光光子与光纤中分子晶格振动声子的非线性相互作用，后散射信号较弱。通过相应的增益技术能够使光信号增强，但是携温信号中夹杂的噪声信号同时也会被放大，所以必须采用一种方法对光电检测器输出信号进行噪声滤除。
现有分布式光纤测温通常采用累加平均滤波、滑动平均滤波、滑动中值滤波等噪声处理算法，上述算法处理时，需要一定时间窗口的信号统计、累加，造成信号窗口内处理数据较多，处理时间长，算法占用硬件资源偏大等问题，最终造成温度数据检测响应速度降低。
发明内容
基于现有技术的上述情况，本发明的目的在于提供一种分布式光纤测温系统自适应滤波方法及装置，能够实现滤波参数的自适应调整，改善滤波效果的同时，确保测温系统的响应灵敏性。
本发明的第一方面提供了一种分布式光纤测温系统自适应滤波方法，包括如下步骤：
对斯托克斯光强、反斯托克斯光强原始采样值进行第一级累加平均滤波；
经过第一级累加平均滤波处理后的斯托克斯光强Vs(T)、反斯托克斯光强Vas(T)的计算值经过数据对齐、去除基线后，计算测温光纤的原始温度计算值；
针对测温光纤指定空间位置连续时间域上的原始温度计算值进行第二级自适应低通滤波处理，更新滤波参数；
针对测温光纤不同空间位置上的原始温度计算值进行第三级自适应低通滤波处理，更新滤波参数。
进一步的，所述对斯托克斯光强、反斯托克斯光强原始采样值进行第一级累加平均滤波的步骤包括：
(1)根据分布式光纤测温系统原始采样率和光纤长度确定重复采样累加频率，对斯托克斯光强、反斯托克斯光强进行重复采样，得到斯托克斯光强、反斯托克斯光强采样数列，采样数列长度为n；
(2)分布式光纤测温系统将重复采样X次得到的斯托克斯光强、反斯托克斯光强采样数列进行累加求平均值；
(3)平均得到的斯托克斯光强、反斯托克斯光强采样数列即第一级累加平均滤波输出结果。
进一步的，测温光纤的原始温度计算值按照下式进行计算：
式中，T是光纤沿线某位置的待测温度，T0是标定位置处的标定温度，h是普朗克常数，Δv是拉曼频移，kB是玻尔兹曼常数，Vas(T)是光纤沿线待测温度下的反斯托克斯光强第一级累加平均值，Vs(T)是光纤沿线待测温度下的斯托克斯光强第一级累加平均值，Vas(T0)是根据标定温度T0经衰减计算的反斯托克斯光强参考值，Vs(T0)是根据标定温度T0经衰减计算的斯托克斯光强参考值。
进一步的，所述针对测温光纤指定空间位置连续时间域上的原始温度计算值进行第二级自适应低通滤波处理，更新滤波参数的步骤包括：
(1)设置温度数据变化阈值Th，滤波参数计数器i＝0；
(2)测温光纤指定空间位置连续时间域上的原始温度计算值作为输入值，减去上次滤波输出结果得到差值及变化方向；
(3)根据数值变化方向及差值大小更新滤波参数计数器；
(4)根据滤波参数计数器i更新自适应滤波参数α，其方法如下：设置数据连续变化计数低门槛M、数据连续变化计数高门槛N，M＜N；若计数器i≥0且计数器i＜低门槛M时，滤波参数增量＝0，即滤波器继续采用原始参数α0；若计数器i≥低门槛M且计数器i＜高门槛N时，则滤波器参数增量Δα与计数器同趋势变步长增大；若计数器i≥高门槛N时，则滤波器参数增量Δα达到最大值并保持不变。
进一步的，所述根据数值变化方向及差值大小更新滤波参数计数器包括：
若数值变化方向一致，则判断数值变化差值是否超过变化阈值Th，如果超过则计数器i加1；若数值变化方向与上次不一致或变化差值未超过变化阈值Th，则将计数器i归零。
进一步的，所述针对测温光纤不同空间位置上的原始温度计算值进行第三级自适应低通滤波处理，更新滤波参数的步骤包括：
(1)设置温度数据变化阈值Th，滤波参数计数器i＝0；
(2)测温光纤不同空间位置上的原始温度计算值作为输入值，减去上次滤波输出结果得到差值及变化方向；
(3)根据数值变化方向及差值大小更新滤波参数计数器；
(4)根据滤波参数计数器i更新自适应滤波参数α，其方法如下：设置数据连续变化计数低门槛M、数据连续变化计数高门槛N，M＜N；若计数器i≥0且计数器i＜低门槛M时，滤波参数增量＝0，即滤波器继续采用原始参数α0；若计数器i≥低门槛M且计数器i＜高门槛N时，则滤波器参数增量Δα与计数器同趋势变步长增大；若计数器i≥高门槛N时，则滤波器参数增量Δα达到最大值并保持不变。
进一步的，所述根据数值变化方向及差值大小更新滤波参数计数器包括：
若数值变化方向一致，则判断数值变化差值是否超过变化阈值Th，如果超过则计数器i加1；若数值变化方向与上次不一致或变化差值未超过变化阈值Th，则将计数器i归零。
本发明的第二方面提供了一种分布式光纤测温系统自适应滤波装置，包括：
第一级滤波处理模块，对斯托克斯光强、反斯托克斯光强原始采样值进行第一级累加平均滤波；
测温光纤的原始温度计算模块，经过第一级累加平均滤波处理后的斯托克斯光强Vs(T)、反斯托克斯光强Vas(T)的计算值经过数据对齐、去除基线后，计算测温光纤的原始温度计算值；
第二级滤波处理模块，针对测温光纤指定空间位置连续时间域上的原始温度计算值进行第二级自适应低通滤波处理，更新滤波参数；
第三级滤波处理模块，针对测温光纤不同空间位置上的原始温度计算值进行第三级自适应低通滤波处理，更新滤波参数；
所述滤波装置用于执行如权利要求1-7任一项所述的滤波方法。
本发明的第三方面提供了一种分布式光纤测温系统自适应滤波系统，该系统包括：
存储器以及一个或多个处理器；
其中，所述存储器与所述一个或多个处理器通信连接，所述存储器中存储有可被所述一个或多个处理器执行的指令，所述指令被所述一个或多个处理器执行，以使所述一个或多个处理器用于执行如前所述的方法。
本发明的第四方面提供了一种计算机可读存储介质，其上存储有计算机可执行指令，当所述计算机可执行指令被计算装置执行时，可操作来执行如前所述的方法。
综上所述，本发明提供一种分布式光纤测温系统自适应滤波方法及装置,该方法包括：对斯托克斯光强、反斯托克斯光强原始采样值进行第一级累加平均滤波；经过第一级累加平均滤波处理后的斯托克斯光强Vs(T)、反斯托克斯光强Vas(T)的计算值经过数据对齐、去除基线后，计算测温光纤的原始温度计算值；针对测温光纤指定空间位置连续时间域上的原始温度计算值进行第二级自适应低通滤波处理，更新滤波参数；针对测温光纤不同空间位置上的原始温度计算值进行第三级自适应低通滤波处理，更新滤波参数。该装置用于执行该方法。与现有技术相比，本发明自适应滤波参数的引入，能够比传统滤波方法更快速地跟踪测量温度的变化趋势，同时实现良好的噪声处理，减少分布式光纤测温系统中具有零均随机噪声造成的干扰。
附图说明
图1为本发明实施例中的时域滤波、空域滤波采样原理示意图；
图2是本发明实施例中的分布式光纤测温系统自适应滤波方法流程示意图；
图3是本发明实施例中的第二级自适应滤波方法流程示意图；
图4是本发明实施例与现有技术滤波后的结果对比示意图；
图5为本发明实施例中的分布式光纤测温系统自适应滤波装置结构框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明了，下面结合具体实施方式并参照附图，对本发明进一步详细说明。应该理解，这些描述只是示例性的，而并非要限制本发明的范围。此外，在以下说明中，省略了对公知结构和技术的描述，以避免不必要地混淆本发明的概念。
本发明的第一方面提供了一种分布式光纤测温系统自适应滤波方法，适用于分布式光纤测温系统空域、时域的多级自适应滤波噪声检测、处理方法，能够实现滤波参数的自适应调整，改善滤波效果的同时，确保测温系统的响应灵敏性。时域滤波、空域滤波如图1所示，时域滤波表示对整条测温光纤的同一空间位置处不同采样时刻的采样值进行滤波；空域滤波表示对整条测温光纤不同空间位置处相同采样时刻的采样值进行滤波。
本发明的分布式光纤测温系统自适应滤波方法，如图2所示，包括如下步骤：
步骤S100，对斯托克斯光强、反斯托克斯光强原始采样值进行第一级累加平均滤波。
步骤S200，经过第一级累加平均滤波处理后的斯托克斯光强Vs(T)、反斯托克斯光强Vas(T)的计算值经过数据对齐、去除基线后，计算测温光纤的原始温度计算值。
步骤S300，针对测温光纤指定空间位置连续时间域上的原始温度计算值进行第二级自适应低通滤波处理，更新滤波参数。
步骤S400，针对测温光纤不同空间位置上的原始温度计算值进行第三级自适应低通滤波处理，更新滤波参数。
具体的，在步骤S100中，对光电检测器输出的两路信号(斯托克斯、反斯托克斯)原始采样值进行第一级累加平均滤波，以系统原始采样率100M，累加频率5kHz为例，完成一次深度1000次的累加平均滤波耗时0.2秒，对应系统温度解调频率为5Hz；第一级累加平均滤波的步骤包括：
(1)根据分布式光纤测温系统原始采样率(例如100MS/s)和光纤长度(例如16km)确定重复采样累加频率(例如5kHz)，对斯托克斯光强、反斯托克斯光强进行重复采样，得到斯托克斯光强、反斯托克斯光强采样数列，单次采样时长不超过0.2ms(即
)，采样数列长度为n；
(2)以累加深度X＝1000为例，分布式光纤测温系统将重复采样X次得到的斯托克斯光强、反斯托克斯光强采样数列进行累加求平均值；
(3)平均得到的斯托克斯光强、反斯托克斯光强采样数列即第一级累加平均滤波输出结果。
具体的，在步骤S200中，第一级累加平均滤波处理后的斯托克斯光强Vs(T)、反斯托克斯光强Vas(T)计算值经过数据对齐、去除基线后，通过下面的温度解调计算公式得到整条测温光纤的原始温度计算值；
式中，T是光纤沿线某位置的待测温度，T0是标定位置处的标定温度，h是普朗克常数(6.62606896×10-34J·s)，Δv是拉曼频移(约1.32×1013Hz)，kB是玻尔兹曼常数(1.3806488×10-23J/K),Vas(T)是光纤沿线待测温度下的反斯托克斯光强第一级累加平均值，Vs(T)是光纤沿线待测温度下的斯托克斯光强第一级累加平均值，Vas(T0)是根据标定温度T0经衰减计算的反斯托克斯光强参考值，Vs(T0)是根据标定温度T0经衰减计算的斯托克斯光强参考值。
具体的，在步骤S300中，测温光纤固定空间位置上的温度测量值受环境温度和光电放大器(APD)、数据采集卡等器件引入的零均值随机噪声影响，上述噪声导致温度测量值围绕实际温度附近随机上下波动。常规低通滤波方法的公式如下：
Yn＝αXn+(1-α)Yn-1 (2)
常规低通滤波方法无法兼顾灵敏度和平稳度，为达到数据相同趋势连续变化时滤波结果能够快速响应，数据围绕某个固定值小幅度波动变化时滤波结果输出尽量平稳，针对测温光纤指定空间位置连续时间域上的原始温度计算值进行第二级的自适应低通滤波，即对不同采样循环上相同采样时刻对应的原始温度计算数据进行低通滤波处理。本发明提出的自适应低通滤波基本实现方法是，如图3所示：
(1)设置温度数据变化阈值Th，滤波参数计数器i＝0；
(2)测温光纤指定空间位置连续时间域上的原始温度计算值作为输入值，减去上次滤波输出结果得到差值及变化方向；
(3)根据数值变化方向及差值大小更新滤波参数计数器，其方法如下：若数值变化方向一致，则判断数值变化差值是否超过变化阈值Th，如果超过则计数器i加1；若数值变化方向与上次不一致或变化差值未超过变化阈值Th，则将计数器i归零(i＝0)；
(4)根据滤波参数计数器i更新自适应滤波参数α，其方法如下：设置数据连续变化计数低门槛M(例如M＝3)、数据连续变化计数高门槛N(例如N＝9)，M＜N；若计数器i≥0且计数器i＜低门槛M时，滤波参数增量＝0，即滤波器继续采用原始参数α0；若计数器i≥低门槛M且计数器i＜高门槛N时，则滤波器参数增量Δα与计数器同趋势变步长增大(例如Δα＝i×变化步长0.1)；若计数器i≥高门槛N时，则滤波器参数增量Δα达到最大值(1-α0)并保持不变。
具体的，在步骤S400中，测温光纤所处的环境温度大体一致，相邻空间位置处的温度一般不会出现剧烈变化，若被测区域存在明显温升时，其对应的光纤温度测量数据相比前后被测区域变化明显，因此噪声处理仍采用自适应低通滤波，兼顾灵敏度和平稳度。
因此解算得到测温光纤原始温度计算值后，对光纤不同位置空间域上的测量温度数值进行第三级自适应低通滤波，滤波器参数自适应更新方法与第二步基本相同，区别在于第二步的输入值为同一空间位置不同时刻的时域值，本步骤的输入值为同一采样周期输出的相邻空间位置上的空域值。具体包括：
(1)设置温度数据变化阈值Th，滤波参数计数器i＝0；
(2)测温光纤不同空间位置上的原始温度计算值作为输入值，减去上次滤波输出结果得到差值及变化方向；
(3)根据数值变化方向及差值大小更新滤波参数计数器，其方法如下：若数值变化方向一致，则判断数值变化差值是否超过变化阈值Th，如果超过则计数器i加1；若数值变化方向与上次不一致或变化差值未超过变化阈值Th，则将计数器i归零(i＝0)；
(4)根据滤波参数计数器i更新自适应滤波参数α，其方法如下：设置数据连续变化计数低门槛M(例如M＝3)、数据连续变化计数高门槛N(例如N＝9)，M＜N；若计数器i≥0且计数器i＜低门槛M时，滤波参数增量＝0，即滤波器继续采用原始参数α0；若计数器i≥低门槛M且计数器i＜高门槛N时，则滤波器参数增量Δα与计数器同趋势变步长增大(例如Δα＝i×变化步长0.1)；若计数器i≥高门槛N时，则滤波器参数增量Δα达到最大值(1-α0)并保持不变。
如图4所示，与现有技术的滤波相比，本发明科学合理的自适应滤波参数的引入，能够比传统滤波方法更快速地跟踪测量温度的变化趋势，同时实现良好的噪声处理，减少分布式光纤测温系统中具有零均随机噪声造成的干扰。
本发明的第二方面提供了一种分布式光纤测温系统自适应滤波装置，如图5所示，包括：第一级滤波处理模块，对斯托克斯光强、反斯托克斯光强原始采样值进行第一级累加平均滤波；测温光纤的原始温度计算模块，经过第一级累加平均滤波处理后的斯托克斯光强Vs(T)、反斯托克斯光强Vas(T)的计算值经过数据对齐、去除基线后，计算测温光纤的原始温度计算值；第二级滤波处理模块，针对测温光纤指定空间位置连续时间域上的原始温度计算值进行第二级自适应低通滤波处理，更新滤波参数；第三级滤波处理模块，针对测温光纤不同空间位置上的原始温度计算值进行第三级自适应低通滤波处理，更新滤波参数。该滤波装置用于执行如前所述的滤波方法，具体不再赘述。
本发明的第三方面提供了一种分布式光纤测温系统自适应滤波系统，该系统包括：存储器以及一个或多个处理器；其中，所述存储器与所述一个或多个处理器通信连接，所述存储器中存储有可被所述一个或多个处理器执行的指令，所述指令被所述一个或多个处理器执行，以使所述一个或多个处理器用于执行如前所述的方法。
本发明的第四方面提供了一种计算机可读存储介质，其上存储有计算机可执行指令，当所述计算机可执行指令被计算装置执行时，可操作来执行如前所述的方法。
综上所述，本发明提供一种分布式光纤测温系统自适应滤波方法及装置,该方法包括：对斯托克斯光强、反斯托克斯光强原始采样值进行第一级累加平均滤波；经过第一级累加平均滤波处理后的斯托克斯光强Vs(T)、反斯托克斯光强Vas(T)的计算值经过数据对齐、去除基线后，计算测温光纤的原始温度计算值；针对测温光纤指定空间位置连续时间域上的原始温度计算值进行第二级自适应低通滤波处理，更新滤波参数；针对测温光纤不同空间位置上的原始温度计算值进行第三级自适应低通滤波处理，更新滤波参数。该装置用于执行该方法。与现有技术相比，本发明自适应滤波参数的引入，能够比传统滤波方法更快速地跟踪测量温度的变化趋势，同时实现良好的噪声处理，减少分布式光纤测温系统中具有零均随机噪声造成的干扰。
应当理解的是，本发明的上述具体实施方式仅仅用于示例性说明或解释本发明的原理，而不构成对本发明的限制。因此，在不偏离本发明的精神和范围的情况下所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。此外，本发明所附权利要求旨在涵盖落入所附权利要求范围和边界、或者这种范围和边界的等同形式内的全部变化和修改例。
Claims (10)
1.一种分布式光纤测温系统自适应滤波方法，其特征在于，包括如下步骤：
对斯托克斯光强、反斯托克斯光强原始采样值进行第一级累加平均滤波；
经过第一级累加平均滤波处理后的斯托克斯光强Vs(T)、反斯托克斯光强Vas(T)的计算值经过数据对齐、去除基线后，计算测温光纤的原始温度计算值；
针对测温光纤指定空间位置连续时间域上的原始温度计算值进行第二级自适应低通滤波处理，更新滤波参数；
针对测温光纤不同空间位置上的原始温度计算值进行第三级自适应低通滤波处理，更新滤波参数。
2.根据权利要求1所述的分布式光纤测温系统自适应滤波方法，其特征在于，所述对斯托克斯光强、反斯托克斯光强原始采样值进行第一级累加平均滤波的步骤包括：
(1)根据分布式光纤测温系统原始采样率和光纤长度确定重复采样累加频率，对斯托克斯光强、反斯托克斯光强进行重复采样，得到斯托克斯光强、反斯托克斯光强采样数列，采样数列长度为n；
(2)分布式光纤测温系统将重复采样X次得到的斯托克斯光强、反斯托克斯光强采样数列进行累加求平均值；
(3)平均得到的斯托克斯光强、反斯托克斯光强采样数列即第一级累加平均滤波输出结果。
3.根据权利要求2所述的分布式光纤测温系统自适应滤波方法，其特征在于，测温光纤的原始温度计算值按照下式进行计算：
式中，T是光纤沿线某位置的待测温度，T0是标定位置处的标定温度，h是普朗克常数，Δv是拉曼频移，kB是玻尔兹曼常数，Vas(T)是光纤沿线待测温度下的反斯托克斯光强第一级累加平均值，Vs(T)是光纤沿线待测温度下的斯托克斯光强第一级累加平均值，Vas(T0)是根据标定温度T0经衰减计算的反斯托克斯光强参考值，Vs(T0)是根据标定温度T0经衰减计算的斯托克斯光强参考值。
4.根据权利要求1-3任一项所述的分布式光纤测温系统自适应滤波方法，其特征在于，所述针对测温光纤指定空间位置连续时间域上的原始温度计算值进行第二级自适应低通滤波处理，更新滤波参数的步骤包括：
(1)设置温度数据变化阈值Th，滤波参数计数器i＝0；
(2)测温光纤指定空间位置连续时间域上的原始温度计算值作为输入值，减去上次滤波输出结果得到差值及变化方向；
(3)根据数值变化方向及差值大小更新滤波参数计数器；
(4)根据滤波参数计数器i更新自适应滤波参数α，其方法如下：设置数据连续变化计数低门槛M、数据连续变化计数高门槛N，M＜N；若计数器i≥0且计数器i＜低门槛M时，滤波参数增量＝0，即滤波器继续采用原始参数α0；若计数器i≥低门槛M且计数器i＜高门槛N时，则滤波器参数增量Δα与计数器同趋势变步长增大；若计数器i≥高门槛N时，则滤波器参数增量Δα达到最大值并保持不变。
5.根据权利要求4所述的分布式光纤测温系统自适应滤波方法，其特征在于，所述根据数值变化方向及差值大小更新滤波参数计数器包括：
若数值变化方向一致，则判断数值变化差值是否超过变化阈值Th，如果超过则计数器i加1；若数值变化方向与上次不一致或变化差值未超过变化阈值Th，则将计数器i归零。
6.根据权利要求1-5任一项所述的分布式光纤测温系统自适应滤波方法，其特征在于，所述针对测温光纤不同空间位置上的原始温度计算值进行第三级自适应低通滤波处理，更新滤波参数的步骤包括：
(1)设置温度数据变化阈值Th，滤波参数计数器i＝0；
(2)测温光纤不同空间位置上的原始温度计算值作为输入值，减去上次滤波输出结果得到差值及变化方向；
(3)根据数值变化方向及差值大小更新滤波参数计数器；
(4)根据滤波参数计数器i更新自适应滤波参数α，其方法如下：设置数据连续变化计数低门槛M、数据连续变化计数高门槛N，M＜N；若计数器i≥0且计数器i＜低门槛M时，滤波参数增量＝0，即滤波器继续采用原始参数α0；若计数器i≥低门槛M且计数器i＜高门槛N时，则滤波器参数增量Δα与计数器同趋势变步长增大；若计数器i≥高门槛N时，则滤波器参数增量Δα达到最大值并保持不变。
7.根据权利要求6所述的分布式光纤测温系统自适应滤波方法，其特征在于，所述根据数值变化方向及差值大小更新滤波参数计数器包括：
若数值变化方向一致，则判断数值变化差值是否超过变化阈值Th，如果超过则计数器i加1；若数值变化方向与上次不一致或变化差值未超过变化阈值Th，则将计数器i归零。
8.一种分布式光纤测温系统自适应滤波装置，其特征在于，包括：
第一级滤波处理模块，对斯托克斯光强、反斯托克斯光强原始采样值进行第一级累加平均滤波；
测温光纤的原始温度计算模块，经过第一级累加平均滤波处理后的斯托克斯光强Vs(T)、反斯托克斯光强Vas(T)的计算值经过数据对齐、去除基线后，计算测温光纤的原始温度计算值；
第二级滤波处理模块，针对测温光纤指定空间位置连续时间域上的原始温度计算值进行第二级自适应低通滤波处理，更新滤波参数；
第三级滤波处理模块，针对测温光纤不同空间位置上的原始温度计算值进行第三级自适应低通滤波处理，更新滤波参数；
所述滤波装置用于执行如权利要求1-7任一项所述的滤波方法。
9.一种分布式光纤测温系统自适应滤波系统，其特征在于，该系统包括：
存储器以及一个或多个处理器；
其中，所述存储器与所述一个或多个处理器通信连接，所述存储器中存储有可被所述一个或多个处理器执行的指令，所述指令被所述一个或多个处理器执行，以使所述一个或多个处理器用于执行权利要求1-7任一项所述的方法。
10.一种计算机可读存储介质，其上存储有计算机可执行指令，当所述计算机可执行指令被计算装置执行时，可操作来执行权利要求1-7任一项所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111333819.6A CN114046900A (zh) | 2021-11-11 | 2021-11-11 | 一种分布式光纤测温系统自适应滤波方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111333819.6A CN114046900A (zh) | 2021-11-11 | 2021-11-11 | 一种分布式光纤测温系统自适应滤波方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114046900A true CN114046900A (zh) | 2022-02-15 |
Family
ID=80208359
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111333819.6A Pending CN114046900A (zh) | 2021-11-11 | 2021-11-11 | 一种分布式光纤测温系统自适应滤波方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114046900A (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH04148832A (ja) * | 1990-10-11 | 1992-05-21 | Asahi Glass Co Ltd | 温度測定方法及び分布型光ファイバー温度センサー |
| CN101943615A (zh) * | 2010-09-29 | 2011-01-12 | 山东大学 | 基于拉曼光反射的测温装置及方法 |
| CN104613321A (zh) * | 2015-01-30 | 2015-05-13 | 武汉工程大学 | 基于分布式光纤测温的核电厂管道泄漏探测装置及方法 |
| CN107044890A (zh) * | 2017-01-09 | 2017-08-15 | 辽宁达能电气股份有限公司 | 长距离自适应的荧光光纤测温装置及方法 |
| CN209605977U (zh) * | 2019-01-02 | 2019-11-08 | 山西天然气有限公司 | 一种含自适应滤波电路的分布式光纤测温系统 |
| CN111289142A (zh) * | 2020-03-18 | 2020-06-16 | 太原理工大学 | 一种分布式光纤拉曼测温系统信号自适应干扰抵消方法 |
-
2021
- 2021-11-11 CN CN202111333819.6A patent/CN114046900A/zh active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH04148832A (ja) * | 1990-10-11 | 1992-05-21 | Asahi Glass Co Ltd | 温度測定方法及び分布型光ファイバー温度センサー |
| CN101943615A (zh) * | 2010-09-29 | 2011-01-12 | 山东大学 | 基于拉曼光反射的测温装置及方法 |
| CN104613321A (zh) * | 2015-01-30 | 2015-05-13 | 武汉工程大学 | 基于分布式光纤测温的核电厂管道泄漏探测装置及方法 |
| CN107044890A (zh) * | 2017-01-09 | 2017-08-15 | 辽宁达能电气股份有限公司 | 长距离自适应的荧光光纤测温装置及方法 |
| CN209605977U (zh) * | 2019-01-02 | 2019-11-08 | 山西天然气有限公司 | 一种含自适应滤波电路的分布式光纤测温系统 |
| CN111289142A (zh) * | 2020-03-18 | 2020-06-16 | 太原理工大学 | 一种分布式光纤拉曼测温系统信号自适应干扰抵消方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| WO2020034558A1 (zh) | 一种面向分布式光纤拉曼传感器的高精度温度解调方法 | |
| CN105157874B (zh) | 一种分布式边界安防监测系统及方法 | |
| CN108414113B (zh) | 用多点温度离散系数预测光纤温度的火灾报警系统及方法 | |
| CN106353003B (zh) | 一种分布式光纤测温方法及系统 | |
| CN104596670A (zh) | 一种解决分布式光纤拉曼温度传感系统温度漂移的方法 | |
| CN110686765B (zh) | 一种基于φ-otdr的输电线路外破监测方法 | |
| CN109752720B (zh) | 基于分布式振动传感器振动信息定位方法及系统 | |
| CN109029770B (zh) | 基于环路解调的分布式光纤拉曼温度及应变解调方法 | |
| CN113820039A (zh) | 一种分布式光纤测温与防外破监测设备及防外破预警方法 | |
| CN114046900A (zh) | 一种分布式光纤测温系统自适应滤波方法及装置 | |
| CN110987038B (zh) | 一种分布式光纤扰动检测装置 | |
| CN113654683A (zh) | 一种分布式光纤测温系统校准方法及装置 | |
| CN111006787B (zh) | 基于差分温度补偿的分布式光纤拉曼双端温度解调方法 | |
| Saxena et al. | Empirical mode decomposition based dynamic error correction in SS covered 62.5/125 µm optical fiber based distributed temperature sensor | |
| CN113091783B (zh) | 基于二级布里渊散射的高灵敏传感装置及方法 | |
| CN112880865B (zh) | 超长距离的高空间分辨率拉曼光纤双参量传感系统和方法 | |
| RU123518U1 (ru) | Волоконно-оптическое устройство акустического мониторинга протяженных объектов | |
| CN108020345A (zh) | 一种基于单模光纤拉曼散射效应的分布式光纤测温装置 | |
| CN103033284A (zh) | 一种基于fft的分布式光纤测温系统降噪方法 | |
| CN108760078B (zh) | 分布式光纤测温方法和系统 | |
| CN112697303A (zh) | 面向智能电网的分布式光纤传感系统和检测方法 | |
| CN112268608A (zh) | 一种强干扰dvs的高精度解调算法 | |
| CN217211121U (zh) | 自适应调整的分布式光纤振动传感器系统 | |
| CN103557877B (zh) | 基于多个edfa在线放大的长距离botdr系统 | |
| CN217084005U (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 |