基于畸变点检测的光纤传感信号相位错误修正方法和装置与流程

1.本发明涉及相位错误抑制技术领域，尤其是涉及基于畸变点检测的光纤传感信号相位错误修正方法和装置。


背景技术：

2.基于光时域反射仪（otdr）的分布式声波传感（das）技术已经被广泛应用于管道健康监测、轨道交通监测、周界入侵监测和水下环境监测等各种领域。其利用激光器向光纤中注入窄线宽的探测脉冲，由于光纤中折射率呈不均匀分布而形成随机的散射点，因此光脉冲在向前传播过程中会产生背向瑞利散射信号（rbs）。外界的扰动会导致光纤的折射率和长度发生变化，从而引起 rbs 的强度和相位发生变化。其中，rbs 的相位延迟的变化量与扰动的强度成正比关系，故可实现定量监测及事件识别的相位敏感光时域反射仪（φ-otdr）被广泛应用。
3.尽管 φ-otdr 具有灵敏度高、传感范围长、不受电磁干扰和分布式传感等多种优点，但是其也面临一些挑战。其中普遍存在且影响不容忽视的难题之一为衰落问题。在φ-otdr 系统中，每一时刻返回的 rbs 其实是光纤上某一段长度为距离分辨率的光纤段上所有散射点的 rbs 的叠加，而 rbs 具有随机性，导致测量得到的 rbs 也具有随机性，即呈瑞利分布。另一方面，散射信号的偏振态也具有随机性，这两个因素共同作用，导致 rbs 的强度会随机上下摆动，因此会随机存在多个 rbs 强度值极小的位置。在这些位置处，信号功率极小，故系统噪声（主要为光电探测中的强度噪声）会严重干扰相位解调的过程，导致相位解调错误，对于常用的相位解调方法，解调结果为卷绕的相位，还需利用相位解卷绕算法来得到真实相位，而相位解调错误会导致相位解卷绕错误，产生接近于2π弧度的错误跳变，甚至造成大幅相位误差积累。相位错误会干扰后续目标定位、识别等任务的正确执行，因此如何切实有效的修正甚至消除这些相位错误是亟待解决的问题。
4.目前已有的相位噪声抑制方法主要可分为两类：硬件增强方法与信号处理方法。在硬件层面，可以通过对探测脉冲进行编码、同时利用不同频率的探测脉冲、使用散射增强光纤等方法抑制衰落噪声引起的相位错误，但这类方法会增加系统的硬件成本及复杂度。在信号处理层面，研究学者提出了一些基于机器学习及深度学习的去噪方法，但这类方法普遍需要配对的干净信号和含噪信号，而在现实环境中难以获取比较干净的传感信号，且通过仿真生成干净数据又需要对实验环境进行建模，难度较大且难以保证与现实条件相符。此外，也有研究学者提出了基于变化点检测的相位错误抑制方法，但其处理规则单一，只适用于简单信号。


技术实现要素：

5.本发明的目的就是为了克服上述现有技术存在实现难度高、或处理规则单一，只适用于简单信号的缺陷而提供一种基于畸变点检测的光纤传感信号相位错误修正方法和装置，设计更加鲁棒的相位错误修正方法，使其可通过信号处理的方式处理复杂信号中可
能的相位错误。
6.本发明的目的可以通过以下技术方案来实现：一种基于畸变点检测的光纤传感信号相位错误修正方法，包括以下步骤：获取光纤传感数据，进行解调，得到相位数据，该相位数据包括光纤上多个位置处的相位变化信号；对所述相位数据中的各个位置信号分别进行相位解卷绕，然后进行重采样，获取各个位置的相位解卷绕数据；根据所述相位解卷绕数据的信号变化量，确定畸变点；根据畸变点的分布确定畸变点的类型，从而进行对应的修正处理。
7.进一步地，所述畸变点的确定过程包括：对所述相位解卷绕数据中的各个信号点进行差分计算，得到差分信号，根据差分信号中各个数值的离散程度，确定相位解卷绕数据中的畸变点。
8.进一步地，确定的所述畸变点类型包括尖峰畸变点，该尖峰畸变点的确定过程包括：在所有确定的畸变点中搜索畸变点对，该畸变点对的搜索过程为：若存在畸变点s
unwrap
(i)以及畸变点s
unwrap
(j+1)，满足时刻与时刻的时间间隔低于预设的尖峰畸变相邻阈值，则(s
unwrap
(i)，s
unwrap
(j))构成畸变点对；若存在畸变点对(s
unwrap
(i)，s
unwrap
(j))，满足s
unwrap
(i-1)信号值与s
unwrap
(j+1)信号值的差值的绝对值小于预设的尖峰畸变阈值，则为尖峰畸变点对。
9.进一步地，若存在两组尖峰畸变点对之间的时间间隔低于预设的尖峰相邻阈值，则仅根据前一组尖峰畸变点进行修正处理。
10.进一步地，对尖峰畸变点的修正处理过程包括：对待处理的尖峰畸变点对(s
unwrap
(i)，s
unwrap
(j))，设置相位解卷绕数据中时刻和j时刻之间的信号值均为s
unwrap
(i-1)信号值与s
unwrap
(j+1)信号值的均值。
11.进一步地，对尖峰畸变点的修正处理过程中，还实时循环分析处理后的光纤传感数据，搜索尖峰畸变点，若搜索出的尖峰畸变点的数量为0，或者搜索出的尖峰畸变点的数量不再继续减少的次数大于预设的稳定阈值，则完成对尖峰畸变点的修正处理，否则继续进行对尖峰畸变点的修正处理。
12.进一步地，确定的所述畸变点类型包括连续台阶畸变点，该连续台阶畸变点的确定过程包括：在所有确定的畸变点中搜索畸变点对(s
unwrap
(i)，s
unwrap
(j))，畸变点对的搜索过程为：若两个畸变点s
unwrap
(i)和s
unwrap
(j)的时间间隔低于预设的台阶畸变相邻阈值，则构成畸变点对；若存在畸变点对(s
unwrap
(i)，s
unwrap
(j))，满足时刻处的差分信号与时刻处的差分信号的乘积大于零，且s
unwrap
(j)信号值与s
unwrap
(i)信号值的差值的绝对值大于预设的台阶畸变阈值，则为台阶畸变点对；判断确定的台阶畸变点对是否存在相邻的台阶畸变点对，若存在，则将相邻的台阶畸变点对依次合并，以合并后的台阶畸变点集合中的首末端点作为连续台阶畸变点对，若不存在，则直接作为连续台阶畸变点对。
13.进一步地，对连续台阶畸变点的修正处理过程包括：对待处理的连续台阶畸变点对(s
unwrap
(i)，s
unwrap
(j))，获取所述相位解卷绕数据在i时刻的前向预测值，以及在j-1时刻的后向预测值，设置相位解卷绕数据中i位置和j位置之间的信号均为时刻的前向预测值；若j-1时刻的后向预测值与i时刻的前向预测值的差值的绝对值不低于预设的第一整体平移调整阈值，则根据j-1时刻的后向预测值与i时刻的前向预测值的差值对所述相位解卷绕数据在j时刻以后的信号值整体进行赋值调整。
14.进一步地，对连续台阶畸变点的修正处理过程中，还实时循环分析处理后的光纤传感数据，搜索连续台阶畸变点对，若搜索出的连续台阶畸变点对的数量为0，或者搜索出的连续台阶畸变点对的数量不再继续减少的次数大于预设的稳定阈值，则完成对连续台阶畸变点对的修正处理，否则继续进行对连续台阶畸变点对的修正处理。
15.进一步地，确定的所述畸变点类型包括第一单次畸变点，该第一单次畸变点的确定过程包括：在所有确定的畸变点中搜索直接相邻的畸变点对(s
unwrap
(i)，s
unwrap
(i+1))，并且畸变点对满足j+1时刻处的差分信号与时刻处的差分信号的乘积小于零；对搜索出的畸变点对进行分组，使得每组中均包含尽可能多的两两直接相邻的畸变点对且其中的畸变点对保持原始顺序，各组之间保持互斥；从所述畸变点对数大于2的各组畸变点中剔除所有处于偶数位或者奇数位的畸变点对，剩余畸变点对中的第一个畸变点作为第一单次畸变点。
16.进一步地，对第一单次畸变点的修正处理过程包括：对每一个第一单次畸变点s
unwrap
(i)，获取相位解卷绕数据中在j-w时刻到i-1时刻之间，以及i+1时刻到j+w时刻之间的信号值，并进行求均值计算，得到信号均值，设置所述第一单次畸变点s
unwrap
(i)的信号值为所述信号均值，w为窗口宽度。
17.进一步地，对第一单次畸变点的修正处理过程中，还实时循环分析处理后的光纤传感数据，搜索第一单次畸变点，若搜索出的第一单次畸变点的数量为0，或者搜索出的第一单次畸变点的数量不再继续减少的次数大于预设的稳定阈值，则完成对第一单次畸变点的修正处理，否则继续进行对第一单次畸变点的修正处理。
18.进一步地，确定的所述畸变点类型包括第二单次畸变点，该第二单次畸变点的确定过程包括：在所有确定的畸变点中搜索畸变点，满足相位解卷绕数据中在该畸变点之后时刻的信号值保持与该畸变点的差值在预设的稳定阈值范围以内，并维持预设的第一稳定时间，将搜索出的畸变点作为第二单次畸变点。
19.进一步地，对第二单次畸变点的修正处理过程包括：对待处理的第二单次畸变点s
unwrap
(j)，获取所述相位解卷绕数据在j时刻的前向预测值，以及在j时刻的后向预测值，若j时刻的后向预测值与j时刻的前向预测值的差值的绝对值不低于预设的第二整体平移调整阈值，则根据j时刻的后向预测值与j时刻的前向预测值的差值对所述相位解卷绕数据在j时刻
以后的信号值整体进行赋值调整。
20.进一步地，对第二单次畸变点的修正处理过程中，还实时循环分析处理后的光纤传感数据，搜索第二单次畸变点，若搜索出的第二单次畸变点的数量为0，或者搜索出的第二单次畸变点的数量不再继续减少的次数大于预设的稳定阈值，则完成对第二单次畸变点的修正处理，否则继续进行对第二单次畸变点的修正处理。
21.进一步地，所述方法还包括：对修正处理后的相位解卷绕数据进行下采样，然后拟合为目标函数，通过该目标函数计算在各个时间点的值，得到相位解卷绕数据的趋势曲线，根据所述趋势曲线对修正处理后的相位解卷绕数据进行信号调整，得到最终的输出相位信号。
22.本发明还提供一种基于畸变点检测的光纤传感信号相位错误修正装置，包括：相位数据获取模块；用于获取光纤传感数据，进行解调，得到相位数据，该相位数据包括光纤上多个位置处的相位变化信号；相位解卷绕模块；用于对所述相位数据中的各个位置信号分别进行相位解卷绕，然后进行重采样，获取各个位置的相位解卷绕数据；畸变点处理模块；用于根据所述相位解卷绕数据的信号变化量，确定畸变点；根据畸变点的分布确定畸变点的类型，从而进行对应的修正处理。
23.本发明还提供一种基于畸变点检测的光纤传感信号相位错误修正装置，包括存储器和一个或多个处理器，所述存储器中存储有可执行代码，所述一个或多个处理器执行所述可执行代码时，用于实现如上所述的基于畸变点检测的光纤传感信号相位错误修正方法。
24.本发明还提供一种计算机可读存储介质，其上存储有程序，该程序被处理器执行时，实现如上所述的基于畸变点检测的光纤传感信号相位错误修正方法。
25.与现有技术相比，本发明具有以下优点：（1）本发明方案可有效消除相位敏感光时域反射仪系统中因相干衰落等原因引起的解卷绕后相位信号中的相位错误。
26.（2）本发明可有效提高相位信号的信噪比，提升信号质量。
27.（3）本发明处理速率快，可部署在系统端，实现对传感信号的实时处理。
[0028] 附图说明
[0029]
图1为本发明实施例中提供的一种基于畸变点检测的光纤传感信号相位错误修正方法的流程示意图；图2为本发明实施例中提供的另一种基于畸变点检测的光纤传感信号相位错误修正方法的流程示意图；图3为本发明实施例中提供的尖峰畸变点的示意图；图4为本发明实施例中提供的连续台阶畸变点的示意图；图5为本发明实施例中提供的单次畸变点的示意图；图6为本发明实施例中提供的原始相位信号示意图；图7为本发明实施例中提供的经解卷绕处理后的相位信号示意图；
图8为本发明实施例中提供的采用本发明方法处理后的相位信号示意图；图9为本发明实施例中提供的一种基于畸变点检测的光纤传感信号相位错误修正装置的结构框图。
实施方式
[0030]
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。
[0031]
因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0032]
应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。
[0033]
需要说明的是，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本技术的描述中，“多个”的含义是两个或两个以上，除非另有明确具体的限定。
实施例
[0034]
图1为本发明一实施例提供的一种基于畸变点检测的光纤传感信号相位错误修正方法的流程示意图。如图1所示，本发明实施例中的基于畸变点检测的光纤传感信号相位错误修正方法可包括如下步骤：s1：获取光纤传感数据，进行解调，得到相位数据，该相位数据包括光纤上多个位置处的相位延迟变化信号；s2：对相位数据中的各个位置信号分别进行相位解卷绕，然后进行重采样，获取各个位置的相位解卷绕数据；s3：根据相位解卷绕数据的信号变化量，确定畸变点；s4：根据畸变点的分布确定畸变点的类型，从而进行对应的修正处理。
[0035]
本方案设计了更加鲁棒的相位错误抑制方法，使其可通过信号处理的方式，基于畸变点检测策略定位并有效处理复杂信号中可能的相位错误，对进一步扩展φ-otdr 的适用场景并提高系统的稳定性具有重要意义。
[0036]
下面对上述各步骤进行具体说明。
[0037]
s1：获取光纤传感数据，进行解调，得到相位数据，该相位数据包括光纤上多个位置处的相位延迟变化信号；本实施例中针对φ-otdr系统的传感数据进行解调，得到相位数据s，其中包含传感光纤上相互间距为距离分辨率r的z个位置处的背向瑞利散射信号的相位延迟变化信号，当信号采样率为f，信号采样时长为t秒时，采样点数为m=f*t，s 是大小为m
×
z的矩阵。
[0038]
s2：对相位数据中的各个位置信号分别进行相位解卷绕，然后进行重采样，获取各个位置的相位解卷绕数据；本实施例中对步骤s1得到的各个位置处的相位传感数据即矩阵s中的每一列数据分别进行相位解卷绕，然后进行重采样，得到s
unwrap
，其长度为n。
[0039]
需要说明的是，上述重采样的目的是使得相位数据中各个采样点数据的时间间隔在预设的标准范围以内，该标准范围优选为0.1ms~10ms，防止采样点的时间间隔过短或过长，并方便后续进行针对性的数据处理。
[0040]
步骤s3和步骤s4分别用于确定畸变点的位置及其类型，从而进行对应的修正处理，本实施例中，提供以下几种实现方式。
[0041]
实现方式一确定的畸变点类型包括尖峰畸变点，该尖峰畸变点的确定过程包括：s301：对相位解卷绕数据中的各个信号点进行差分计算，得到差分信号，根据差分信号中各个数值的离散程度，确定相位解卷绕数据中的畸变点。
[0042]
可选的，上述差分计算可以为，将相位解卷绕数据中各个时刻的信号点依次与前一时刻的信号点相减，得到差分计算结果。
[0043]
可选的，根据差分信号中各个数值的离散程度，确定相位解卷绕数据中的畸变点，可以为将差分信号中的各个数值与相位解卷绕数据的标准差对比，从而判断离散程度，确定畸变点。
[0044]
本实施例中，对经前述步骤处理后得到的s
unwrap
进行差分运算，即对于光纤上第c个位置处的相位解卷绕信号s
unwrap
(:，c)，c=1,2,
…
,m，（为便于表示，后续内容中s
unwrap
(:，c)将简写为s
unwrap
）求差分后得到信号d，其长度为n-1，且d(i)=s
unwrap
(i+1)-s
unwrap
(i)，i=1,
…
,n-1。选取某一阈值δ，当d 中某一元素d(i)的绝对值|d(i)|大于δ时，则s
unwrap
中第i+1个相位信号发生了畸变，据此可以找出s
unwrap
中所有的畸变点。其中，经解卷绕算法处理后，，故0《δ《π，考虑到外界振动声波干扰光纤时将持续一段时间，即使该声波信号持续时间非常短，也将远大于探测信号的周期，因此正常情况下相位不会发生单点畸变情况，所以当某一个相位值的变化量大于信号的标准差时，可以判定其为畸变点，故std(s
unwrap
)《δ《π。
[0045]
s302：在所有确定的畸变点中搜索畸变点对(s
unwrap
(i)，s
unwrap
(j))，畸变点对的搜索过程为：若两个畸变点的时间间隔低于预设的尖峰畸变相邻阈值δ1，则构成畸变点对；尖峰畸变点表现为信号从正常值跳跃至一异常值后立即跳回正常值或保持几步异常值后再跳回正常值附近；因此通过获取相邻的畸变点构成的畸变点对的方式可以搜索出具有持续异常过程的畸变点，并且当j=i 时，s
unwrap
(i)和s
unwrap
(j+1)为a中直接相邻的一对畸变点，(s
unwrap
(i)，s
unwrap
(j))将表示单独的一点s
unwrap
(i)。
[0046]
尖峰畸变相邻阈值δ1的取值可以根据相位解卷绕数据的采样点时间间距进行对应调整，以获得尽可能多的近似相邻畸变点。
[0047]
本实施例中，从所有的畸变点中从前向后依次挑选出所有的相邻的畸变点对，其中当a中的两个畸变点的间距不大于
阈值δ1+1时，便可认为其相邻，当δ1=0时，即当j=i 时，s
unwrap
(i)和s
unwrap
(j+1)为a中直接相邻的一对畸变点，(s
unwrap
(i)，s
unwrap
(j))将表示单独的一点s
unwrap
(i)（对应于图3的上半部分），为找出近似相邻的畸变点对（对应于图 3的下半部分），可适当增大δ1的取值，经验证δ1取值在10左右时便可找出信号中绝大多数的近似相邻畸变点。
[0048]
s303：若存在畸变点对(s
unwrap
(i)，s
unwrap
(j))，满足s
unwrap
(i-1)信号值与s
unwrap
(j+1)信号值差值的绝对值小于预设的尖峰畸变阈值，则为尖峰畸变点对。
[0049]
可选的，若存在两组尖峰畸变点对之间的时间间隔低于预设的尖峰相邻阈值，则仅根据前一组尖峰畸变点进行修正处理，后面的尖峰畸变点会自行恢复正常。
[0050]
本实施例中，尖峰畸变点可通过对b进行筛选来得到，即，其中，阈值δ2用于限制该畸变点对前后相邻信号值的差值应在0附近，且当b 中存在前后相邻的两组尖峰畸变点对(s
unwrap
(i)，s
unwrap
(j))和(s
unwrap
(j+1)，s
unwrap
(k))时，c 中将只包含(s
unwrap
(i)，s
unwrap
(j))。至此，c 中即包含了信号中所有的尖峰畸变点。
[0051]
s304：循环分析光纤传感数据，根据上述步骤搜索尖峰畸变点，若搜索出的尖峰畸变点的数量为0，或者搜索出的尖峰畸变点的数量不再继续减少的次数大于预设的稳定阈值，则完成对尖峰畸变点的修正处理，否则继续进行对尖峰畸变点的修正处理。
[0052]
本实施例中，实时判断c 中尖峰畸变点的个数是否为 0，如果为 0，则完成对尖峰畸变点的处理，可进行其它类型畸变点的处理，或结束对畸变点的处理，如果不为 0，则判断其数量不再继续减少的次数是否大于阈值m，如果大于，则完成对尖峰畸变点的处理，可进行其它类型畸变点的处理，或结束对畸变点的处理，否则跳至步骤s305。其中阈值m用于防止处理程序陷入死循环，其取值过小容易导致对畸变点的处理不彻底，但取值过大又会减慢处理程序的运行速率，经验证10次循环后基本可以解决畸变点出现的情形，且处理速率较快，故推荐的取值为10左右。
[0053]
s305：对尖峰畸变点进行修正处理过程，具体包括：对待处理的尖峰畸变点对(s
unwrap
(i)，s
unwrap
(j))，设置相位解卷绕数据中时刻和时刻之间的信号值均为s
unwrap
(i-1)信号值与s
unwrap
(j+1)信号值的均值。
[0054]
本实施例中，对c 中的每一对尖峰畸变点(s
unwrap
(i)，s
unwrap
(j))，执行
[0055]
然后跳转至步骤s301。
[0056]
实现方式二确定的畸变点类型包括连续台阶畸变点，该连续台阶畸变点的确定过程包括：s311：对相位解卷绕数据中的各个信号点进行差分计算，得到差分信号，根据差分信号中各个数值的离散程度，确定相位解卷绕数据中的畸变点，与步骤s301大体相同。
[0057]
本实施例中，对经上述步骤处理后的相位信号s
unwrap
进行差分运算得到信号d，其长度为t*f-1且。选取某一阈值 δ（其取值与步骤s301中相同），当d 中某一元素d(i)的绝对值|d(i)|大于δ时，便认为s
unwrap
(i+1)中对应时间点处
的相位信号发生了畸变，据此可以找出s
unwrap
中所有的畸变点。
[0058]
s312：在所有确定的畸变点中搜索畸变点对(s
unwrap
(i)，s
unwrap
(j))，畸变点对的搜索过程为：若两个畸变点的时间间隔低于预设的台阶畸变相邻阈值δ3，则构成畸变点对；连续台阶畸变点表现为信号从正常值跳跃至一异常值后，立即或保持一定步数后再沿相同方向跳跃至更远的异常值，然后停止跳跃或继续此跳跃过程。
[0059]
由于连续台阶畸变点需产生连续变化，即连续分布有多个畸变点，因此上述畸变点对的搜索过程需使得畸变点对包含两个不同且相邻的畸变点。
[0060]
本实施例中，从所有的畸变点中从前向后依次挑选出所有的相邻畸变点对，其中当d中的两个畸变点的间距不大于阈值δ3时，便可认为其相邻，当δ3=1时，s
unwrap
(i)和s
unwrap
(j)为d中直接相邻的一对畸变点（对应于图4的上半部分），为找出近似相邻的畸变点对（对应于图4的下半部分），可适当增大δ3的取值，经验证δ3取值在10左右时便可以找出信号中绝大多数的近似相邻畸变点。
[0061]
s313：若存在畸变点对(s
unwrap
(i)，s
unwrap
(j))，满足i-1时刻处的差分信号与j-1时刻处的差分信号的乘积大于零，且s
unwrap
(j)信号值与s
unwrap
(i)信号值的差值的绝对值大于预设的台阶畸变阈值δ4，则为台阶畸变点对。
[0062]
其中，满足i-1时刻处的差分信号与j-1时刻处的差分信号的乘积大于零的目的是确定畸变点沿相同方向进行跳跃，台阶畸变阈值δ4的作用是确保第二个畸变点较第一个畸变点有明显变化，实现呈阶梯状变化。
[0063]
本实施例中，对连续台阶畸变点进行定位需对e 进行筛选，得到 ，其中为确保第二个畸变点较第一个畸变点呈阶梯状变化，阈值δ4取值应大于0，同时考虑到信号值在第一个畸变点后可能存在小幅的变化，故阈值δ4的取值应小于选取畸变点的阈值δ，综合考虑，其推荐取值范围可为0.5≤δ4《δ。
[0064]
s314：判断确定的台阶畸变点对中是否存在相邻的台阶畸变点对，若存在，则将相邻的台阶畸变点对依次合并，以合并后的台阶畸变点集合中的首末端点作为连续台阶畸变点对，若不存在，则直接作为连续台阶畸变点对。
[0065]
对于连续台阶畸变点，其处理过程可根据整体变化的端值进行处理，因此进行上述合并操作，减少畸变点修正次数，提高处理效率。
[0066]
本实施例中，合并f 中直接相邻的畸变点对，以合并后的台阶畸变点集合中的首末端点作为连续台阶畸变点对，即进行得到
[0067]
至此，g 中即包含了信号中所有的连续台阶畸变点。
[0068]
s315：循环分析光纤传感数据，根据上述步骤搜索连续台阶畸变点对，若搜索出的连续台阶畸变点对的数量为0，或者搜索出的连续台阶畸变点对的数量不在继续减少的次数大于预设的稳定阈值，则完成对连续台阶畸变点对的修正处理，否则继续进行对连续台
阶畸变点对的修正处理。
[0069]
本实施例中，判断g 中台阶畸变点的个数是否为 0，如果为 0，则完成对连续台阶畸变点对的处理，可进行其它类型畸变点的处理，或结束对畸变点的处理，如果不为0，则判断其数量不再继续减少的次数是否大于阈值s，如果大于，则完成对连续台阶畸变点对的处理，可进行其它类型畸变点的处理，或结束对畸变点的处理，否则跳至步骤s316。其中阈值s用于防止处理程序陷入死循环，其取值过小容易导致对畸变点的处理不彻底，但取值过大又会减慢处理程序的运行速率，经验证10次循环后基本可以解决畸变点出现的情形，且处理速率较快，故推荐的取值为10左右。
[0070]
s316：对连续台阶畸变点进行修正处理，过程具体为：对待处理的连续台阶畸变点对(s
unwrap
(i)，s
unwrap
(j))，获取所述相位解卷绕数据在i时刻的前向预测值，以及在j-1时刻的后向预测值，设置相位解卷绕数据中i位置和j位置之间的信号均为时刻的前向预测值；若j-1时刻的后向预测值与i时刻的前向预测值的差值的绝对值不低于预设的第一整体平移调整阈值，则根据j-1时刻的后向预测值与i时刻的前向预测值的差值对所述相位解卷绕数据在j时刻以后的信号值整体进行赋值调整。
[0071]
首先对连续台阶畸变点对进行前向预测以及后向预测，从而对连续台阶畸变点对对应的时间范围进行重新赋值，对于连续台阶畸变后的相位信号值均大幅偏移的情况，进行了整体平移调整。
[0072]
根据时刻的后向预测值与时刻的前向预测值的差值对相位解卷绕数据在时刻以后的信号值整体进行赋值调整，不局限于以计算的差值作为时刻以后的信号值的调整量，也可以对计算的差值进行调整，如增加对差值乘以自定义系数等。
[0073]
上述对时刻以后的信号值整体赋值，可以将时刻以后的部分或全部信号值整体调整相同的量，也可以赋予不同值，以达到连续台阶畸变优化效果。
[0074]
上述前向预测值以及后向预测值的预测方法为通过前向正常数据或后向正常数据进行数据预测，可以通过现有方法实现，不做唯一限定。
[0075]
本实施例中，对g 中的每一对连续台阶畸变点(s
unwrap
(i)，s
unwrap
(j))，进行如下处理：
[0076][0077]
其中 为信号在时间点i 的前向预测值，具体地，首先从信号点s
unwrap
(i-1)开始向前截取长度为l
1 的信号片段，然后对该信号片段进行拟合，通过优化得到目标方程后再对时间点i 的信号点的值进行预测。同理，为信号在时间点j-1 的后向预测值，具体地，首先从信号点s
unwrap
(j)开始向后截取长度为l
2 的信号片段，然后对该信号片段进行拟合，通过优化得到目标方程后再对时间点j-1 的信号点的值进行预测。阈值δ5用于判断第二个畸变点后的信号值是否需要进行整体平移调整，其取值可依据阈值δ来确定，考虑
到前向预测值可能会较s
unwrap
(i)更接近于，故阈值δ5的取值应小于选取畸变点的阈值δ，经验证δ5取值在0.5δ左右时可以对连续台阶畸变点起到较好的检测效果。
[0078]
实现方式三确定的畸变点类型包括第一单次畸变点，该第一单次畸变点的确定过程包括：s321：对相位解卷绕数据中的各个信号点进行差分计算，得到差分信号，根据差分信号中各个数值的离散程度，确定相位解卷绕数据中的畸变点，与步骤s301大体相同。
[0079]
本实施例中，对经上述步骤处理后的相位信号s
unwrap
进行差分运算得到信号d，其长度为n-1，且 。选取某一阈值 δ（其取值与步骤三中相同），当d 中某一元素d(i)的绝对值|d(i)|大于 δ 时，便认为s
unwrap
(i+1)中对应时间点处的相位信号发生了畸变，据此可以找出s
unwrap
中所有的畸变点。
[0080]
s322：在所有确定的畸变点中搜索直接相邻的畸变点对(s
unwrap
(i)，s
unwrap
(i+1))，并且畸变点对满足i-1时刻处的差分信号与时刻处的差分信号的乘积小于零；对搜索出的畸变点进行分组，使得每组中均包含尽可能多的两两直接相邻的畸变点对且其中的畸变点对保持原始顺序，各组之间保持互斥；从所述畸变点对数大于2的各组畸变点中剔除所有处于偶数（或奇数）位置的畸变点，剩余畸变点对中的第一个畸变点作为第一单次畸变点。
[0081]
第一单次畸变点表现为信号从某一正常值沿一方向跳跃至一异常值，然后立即沿反方向跳跃至另一异常值，如图5的上半部分所示。
[0082]
因此，本方法搜索直接相邻的畸变点对(s
unwrap
(i)，s
unwrap
(i+1))，并且畸变点对满足i-1时刻处的差分信号与时刻处的差分信号的乘积小于零；两畸变点之间的时间间隔为重采样的单位时间长度时，即为直接相邻的畸变点；限制i-1时刻处的差分信号与时刻处的差分信号的乘积小于零，能保证时刻的畸变方向与i+1时刻的畸变方向相反，实现沿一方向跳跃至一异常值，然后立即沿反方向跳跃至另一异常值。
[0083]
对搜索出的多个畸变点构成的畸变点集合，依次等时间间隔剔除其中的畸变点，可获得多个第一单次畸变点，并对畸变情况均能进行处理。
[0084]
本实施例中，从所有的畸变点中从前向后依次挑选出所有直接相邻的畸变点对，i 中如果存在直接相邻的畸变点，则应依次等间隔剔除其中的点。至此，i 中将包含信号中的第一单次畸变点。
[0085]
s323：循环分析光纤传感数据，根据上述步骤搜索第一单次畸变点，若搜索出的第一单次畸变点的数量为0，或者搜索出的第一单次畸变点的数量不在继续减少的次数大于预设的稳定阈值，则完成对第一单次畸变点的修正处理，否则继续进行对第一单次畸变点的修正处理。
[0086]
本实施例中，判断h 中的畸变点的个数是否为 0，如果为 0，则完成对第一单次畸变点的处理，可进行其它类型畸变点的处理，或结束对畸变点的处理，如果不为 0，则判断其数量不再继续减少的次数是否大于阈值u，如果大于，则完成对第一单次畸变点的处理，可进行其它类型畸变点的处理，或结束对畸变点的处理，否则跳至步骤s324。其中阈值u用
于防止处理程序陷入死循环，其取值过小容易导致对畸变点的处理不彻底，但取值过大又会减慢处理程序的运行速率，经验证10次循环后基本可以解决畸变点出现的情形，且处理速率较快，故推荐的取值为10左右。
[0087]
s324：对第一单次畸变点进行修正处理，过程包括：对每一个第一单次畸变点s
unwrap
(i)，获取相位解卷绕数据中在i-w时刻到i-1时刻之间，以及i+1时刻到j+w时刻之间的信号值，并进行求均值计算，得到信号均值，设置第一单次畸变点s
unwrap
(i)的信号值为信号均值，w为窗口宽度。
[0088]
可通过调整窗口宽度w，实现获取最优的信号均值。
[0089]
本实施例中，对i 中的每一畸变点s
unwrap
(i)进行如下处理：
[0090]
其中参数w 用于控制窗口的宽度，mean() 表示求均值运算。
[0091]
实现方式四确定的畸变点类型包括第二单次畸变点，该第二单次畸变点的确定过程包括：s331：对相位解卷绕数据中的各个信号点进行差分计算，得到差分信号，根据差分信号中各个数值的离散程度，确定相位解卷绕数据中的畸变点，与步骤s301大体相同。
[0092]
本实施例中，对经上述步骤处理后的相位信号s
unwrap
进行差分运算得到信号d，其长度为n-1，且 。选取某一阈值 δ（其取值与步骤s301中相同），当d 中某一元素d(i)的绝对值|d(i)|大于 δ 时，便认为s
unwrap
(i+1)中对应时间点处的相位信号发生了畸变，据此可以找出 中所有的畸变点。
[0093]
s332：在所有确定的畸变点中搜索畸变点，满足相位解卷绕数据中在该畸变点之后时刻的信号值保持与该畸变点的差值在预设的稳定阈值范围以内，并维持预设的第一稳定时间，将搜索出的畸变点作为第二单次畸变点。
[0094]
信号中存在第二单次畸变点，这类畸变点表现为信号从正常值跳跃至一异常值后便保持在该值附近一段时间而不再继续跳变（附图5下半部分），而这类畸变点将包含在集合h-i 中。
[0095]
s333：循环分析光纤传感数据，根据上述步骤搜索第二单次畸变点，若搜索出的第二单次畸变点的数量为0，或者搜索出的第二单次畸变点的数量不在继续减少的次数大于预设的稳定阈值，则完成对第二单次畸变点的修正处理，否则继续进行对第二单次畸变点的修正处理。
[0096]
本实施例中，判断h 中的畸变点的个数是否为 0，如果为 0，则完成对第二单次畸变点的处理，可进行其它类型畸变点的处理，或结束对畸变点的处理，如果不为 0，则判断其数量不再继续减少的次数是否大于阈值u，如果大于，则完成对第二单次畸变点的处理，可进行其它类型畸变点的处理，或结束对畸变点的处理，否则跳至步骤s334。其中阈值u用于防止处理程序陷入死循环，其取值过小容易导致对畸变点的处理不彻底，但取值过大又会减慢处理程序的运行速率，经验证10次循环后基本可以解决畸变点出现的情形，且处理速率较快，故推荐的取值为10左右。
[0097]
s334：对第二单次畸变点进行修正处理，过程包括：对待处理的第二单次畸变点s
unwrap
(j)，获取所述相位解卷绕数据在j时刻的前向
预测值，以及在j时刻的后向预测值，若j时刻的后向预测值与j时刻的前向预测值的差值的绝对值不低于预设的第二整体平移调整阈值，则根据j时刻的后向预测值与j时刻的前向预测值的差值对所述相位解卷绕数据在j时刻以后的信号值整体进行赋值调整。
[0098]
第二单次畸变点从正常值跳跃至一异常值后便保持在该值附近一段时间而不再继续跳变，因此需对跳变后的相位信号值，进行整体平移调整。
[0099]
根据时刻的后向预测值与时刻的前向预测值的差值对相位解卷绕数据在时刻以后的信号值整体进行赋值调整，不局限于以计算的差值作为时刻以后的信号值的整体调整量，也可以对计算的差值进行调整，如增加对差值乘以自定义系数等。
[0100]
上述对时刻以后的信号值整体赋值，可以将时刻以后的部分或全部信号值调整相同量，也可以赋予不用值，以达到第二单次畸变优化效果。
[0101]
上述前向预测值以及后向预测值的预测方法为通过前向正常数据或后向正常数据进行数据预测，可以通过现有方法实现，不做唯一限定。
[0102]
本实施例中，对h-i 中的每一个畸变点s
unwrap
(j)，进行如下处理：
[0103]
其中 为信号在时间点j 的前向预测值，具体地，首先从信号点s
unwrap
(j-1)开始向前截取长度为l
3 的信号片段，然后对该信号片段进行拟合，通过优化得到目标方程后再对时间点j 的信号点的值进行预测。同理，为信号在时间点j 的后向预测值，具体地，首先从信号点s
unwrap
(j+1) 开始向后截取长度为l
4 的信号片段，然后对该信号片段进行拟合，通过优化得到目标方程后再对时间点j 的信号点的值进行预测。阈值δ6用于判断第二个畸变点后的信号值是否需要进行整体平移调整，由于只发生了一次畸变，其取值可接近于阈值δ。
[0104]
可选的，方法还包括：对修正处理后的相位解卷绕数据进行下采样，然后拟合为目标函数，通过该目标函数计算在各个时间点的值，得到相位解卷绕数据的趋势曲线，根据趋势曲线对修正处理后的相位解卷绕数据进行信号调整，得到最终的输出相位信号。
[0105]
经上述步骤处理后的数据可能存在稳定状态下的相位值不为0的情况，为使其回归0值附近，需进行去趋势操作。
[0106]
本实施例的具体过程如下：选取某一步长s 对经上述步骤处理后的s
unwrap
(i)进行下采样得到 ，然后对其进行拟合得到目标函数g
fit
后，再计算出g
fit
在各时间点1,2,3，
…
n 处的值，即可得到s
unwrap
的趋势曲线 。其中步长s 可通过分析相位信号中有用信号的频率范围来确定。最终，去除趋势后的相位信号即为最终的输出相位信号：。
[0107]
将上述可选的实施方式以及对畸变点处理的各个实现方式进行任意组合可以得
circuit hardware description language）与verilog。本领域技术人员也应该清楚，只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中，就可以很容易得到实现该逻辑方法流程的硬件电路。
[0113]
控制器可以按任何适当的方式实现，例如，控制器可以采取例如微处理器或处理器以及存储可由该（微）处理器执行的计算机可读程序代码（例如软件或固件）的计算机可读介质、逻辑门、开关、专用集成电路（application specific integrated circuit，asic）、可编程逻辑控制器和嵌入微控制器的形式，控制器的例子包括但不限于以下微控制器：arc 625d、atmel at91sam、microchip pic18f26k20 以及silicone labs c8051f320，存储器控制器还可以被实现为存储器的控制逻辑的一部分。本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
[0114]
上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的，计算机例如可以为个人计算机、膝上型计算机、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
[0115]
为了描述的方便，描述以上装置时以功能分为各种单元分别描述。当然，在实施本发明时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
[0116]
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质（包括但不限于磁盘存储器、cd-rom、光学存储器等）上实施的计算机程序产品的形式。
[0117]
本发明是参照根据本发明实施例的方法、设备（系统）、和计算机程序产品的流程图和／或方框图来描述的。应理解可由计算机程序指令实现流程图和／或方框图中的每一流程和／或方框、以及流程图和／或方框图中的流程和／或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和／或方框图一个方框或多个方框中指定的功能的装置。
[0118]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和／或方框图一个方框或多个方框中指定的功能。
[0119]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和／或方框图一
个方框或多个方框中指定的功能的步骤。
[0120]
在一个典型的配置中，计算设备包括一个或多个处理器(cpu)、输入/输出接口、网络接口和内存。
[0121]
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(ram)和/或非易失性内存等形式，如只读存储器(rom)或闪存(flash ram)。内存是计算机可读介质的示例。
[0122]
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(pram)、静态随机存取存储器(sram)、动态随机存取存储器(dram)、其他类型的随机存取存储器(ram)、只读存储器(rom)、电可擦除可编程只读存储器(eeprom)、快闪记忆体或其他内存技术、只读光盘只读存储器(cd-rom)、数字多功能光盘(dvd)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
[0123]
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个
……”
限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
[0124]
本领域技术人员应明白，本发明的实施例可提供为方法、系统或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质（包括但不限于磁盘存储器、cd-rom、光学存储器等）上实施的计算机程序产品的形式。
[0125]
本发明可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本发明，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
[0126]
本发明中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
[0127]
以上详细描述了本发明的较佳具体实施例。应当理解，本领域的普通技术人员无需创造性劳动就可以根据本发明的构思做出诸多修改和变化。因此，凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案，皆应在由权利要求书所确定的保护范围内。