一种基于极限逼近的负压波信号拐点定位方法及系统

1.本发明涉及管道泄漏检测技术领域，具体涉及一种基于极限逼近的负压波信号拐点定位方法及系统。


背景技术：

2.管道泄漏之后，由于泄漏点处流体损失，会产生一种负压波并以一定的速度向泄漏点的上下游传播，引起管道内部压力的瞬间降低。由于管道内部压力在正常运营状态下通常处于稳定状态，因此负压波的传播会导致管道内部压力或者管壁环向应变的时程曲线出现拐点，即负压波信号拐点。
3.通过计算不同传感器测量到拐点的时间差，就能计算管道泄漏点的位置。这种方法由于原理简单，易于操作，是目前广泛使用的一种管道泄漏定位方法。但是这种方法对拐点定位的精度有较高的要求，否则就会导致泄漏定位误差较大。而在实际工程中，由于受到外界各种因素影响，测得的数据往往呈现震荡趋势，一些常用的拐点定位方法如二次求导法就不再适用处理这种具有震荡趋势的数据，原因在于对具有震荡趋势的数据进行拐点定位时，需要进行线性拟合，操作繁琐且误差较大。
4.目前关于负压波信号的拐点定位方法还较少，主要有阈值法和奇异点定位方法。阈值法是设定一个阈值，超过阈值便视为拐点；这种方法虽然比较简单，但过于保守，拐点定位误差非常大。奇异点定位法的原理是信号的突变用小波的奇异点进行表征，确定奇异点的位置就能确定拐点的位置。然而多尺度小波变换的误差会随着尺度的增加而迅速增大，并且还具有工作量大，计算时间长，不能直接定位拐点位置等缺点。
5.因此，如何精确、迅速地定位负压波信号拐点具有重要的意义。


技术实现要素：

6.本发明目的在于提供一种基于极限逼近的负压波信号拐点定位方法及系统，以解决现有技术中所存在的一个或多个技术问题，至少提供一种有益的选择或创造条件。
7.为了实现上述目的，本发明提供以下技术方案：
8.一种基于极限逼近的负压波信号拐点定位方法，所述方法包括以下步骤：
9.接收传感器采集到的检测信号，所述检测信号包括检测信号值随时间的对应关系；
10.根据所述检测信号建立信号时程曲线；
11.若确定管道发生泄漏，则将所述信号时程曲线划分为平稳段、待测段和下降段；
12.在所述平稳段选取任一点作为起始点a，在所述下降段选取任一点作为结束点b，分别确定起始点a和结束点b所对应的坐标；
13.将a点与b点之间的所有点作为待测点，分别求解每个待测点的第一斜率值和第二斜率值；其中，所述第一斜率值为待测点与a点连线的斜率值，所述第二斜率值为待测点与b点连线的斜率值；
14.确定每个待测点的斜率差，所述斜率差根据第二斜率值减去第一斜率值得到；
15.确定斜率差最大的待测点，将斜率差最大的待测点作为负压波信号拐点。
16.进一步，所述检测信号包括管道泄漏所产生的压力信号或者管壁的环向应变信号。
17.进一步，所述若确定管道发生泄漏，则将所述信号时程曲线划分为平稳段、待测段和下降段，包括：
18.确定第一阈值，所述第一阈值为管道正常状态下多个检测信号值的平均值；
19.按时间顺序从信号时程曲线中依次选取一段连续数量的检测信号值，当确定所述连续数量的检测信号值均低于第一阈值时，则将该连续数量的检测信号值中的第一个检测信号值所对应的时间作为第一时间，将第一时间之前的信号时程曲线划分为平稳段；
20.按时间顺序继续从信号时程曲线中依次选取一段连续数量的检测信号值，当确定所述连续数量的检测信号值均低于第二阈值时，则将该连续数量的检测信号值中的最后一个检测信号值所对应的时间作为第二时间，将第二时间之后的信号时程曲线划分为下降段；其中，所述第二阈值小于第一阈值；
21.将第一时间和第二时间之间的信号时程曲线划分为待测段。
22.一种计算机可读存储介质，所述计算机可读存储介质上存储有基于极限逼近的负压波信号拐点定位程序，所述基于极限逼近的负压波信号拐点定位程序被处理器执行时实现如上述任意一项所述的基于极限逼近的负压波信号拐点定位方法的步骤。
23.一种基于极限逼近的负压波信号拐点定位系统，所述终端包括：
24.至少一个处理器；
25.至少一个存储器，用于存储至少一个程序；
26.当所述至少一个程序被所述至少一个处理器执行，使得所述至少一个处理器实现上述任一项所述的基于极限逼近的负压波信号拐点定位方法。
27.本发明的有益效果是：本发明公开一种基于极限逼近的负压波信号拐点定位方法及系统，本发明仅通过求解斜率差极值的方式就能计算负压波拐点的位置。本发明具有原理简单，计算精度以及效率高，适合实时自动化监测等优点。
附图说明
28.为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
29.图1是本发明实施例中基于极限逼近的负压波信号拐点定位方法的流程示意图；
30.图2是本发明实施例中信号时程曲线的示意图。
具体实施方式
31.以下将结合实施例和附图对本技术的构思、具体结构及产生的技术效果进行清楚、完整的描述，以充分地理解本技术的目的、方案和效果。需要说明的是，在不冲突的情况下，本技术中的实施例及实施例中的特征可以相互组合。
32.参考图1，如图1所示为本技术实施例提供的一种基于极限逼近的负压波信号拐点定位方法，所述方法包括以下步骤：
33.步骤s100、接收传感器采集到的检测信号，所述检测信号包括检测信号值随时间的对应关系；
34.在一个优选的实施例中，所述检测信号包括管道泄漏所产生的压力信号或者管壁的环向应变信号；通过在管道的不同位置设置传感器，可以采集到管道上不同位置的检测信号。
35.步骤s200、根据所述检测信号建立信号时程曲线；
36.可以理解，传感器采集到的检测信号为离散信号，为便于后续处理，将离散的检测信号通过连线，得到信号时程曲线。
37.步骤s300、若确定管道发生泄漏，则将所述信号时程曲线划分为平稳段、待测段和下降段；
38.步骤s400、在所述平稳段选取任一点作为起始点a，在所述下降段选取任一点作为结束点b，分别确定起始点a和结束点b所对应的坐标；
39.本实施例中，所述起始点a所对应的坐标包括起始点a的时间点和信号值，所述结束点b所对应的坐标包括结束点b的时间点和信号值。
40.本实施例中，选择平稳段中第2500个点作为a点，从检测信号数据中获得a点对应的时刻值t
a
与信号值y
a
。从下降段取一点b，从检测信号数据中获得b点对应的时刻点t
b
与信号值y
b
。
41.步骤s500、将a点与b点之间的所有点作为待测点，分别求解每个待测点的第一斜率值和第二斜率值；其中，所述第一斜率值为待测点与a点连线的斜率值，所述第二斜率值为待测点与b点连线的斜率值；
42.在一实施例中，设所取的待测点坐标为(t
i
,y
i
)，a点坐标为(t
a
,y
a
),b点坐标为(t
b
,y
b
)；
43.则第一斜率值通过公式计算得到，第二斜率值通过公式计算得到。
44.步骤s600、确定每个待测点的斜率差，所述斜率差根据第二斜率值减去第一斜率值得到；
45.在一实施例中，待测点的斜率差通过n
i
＝|k
bi
|
‑
|k
ai
|公式计算得到，其中，n
i
为第i个待测点的斜率差，k
ai
为第i个待测点的第一斜率值，k
bi
为第i个待测点的第二斜率值，1≤i≤n，n为a点与b点之间的总点数；
46.步骤s700、确定斜率差最大的待测点，将斜率差最大的待测点作为负压波信号拐点。
47.本实施例中，通过对检测信号进行分析，将斜率差最大的待测点作为所述信号时程曲线中的拐点，将信号时程曲线中拐点对应的时刻点作为负压波信号拐点出现的时间。
48.在一实施例中，对每个待测点的斜率差按大小进行排序，得到最大值n
i,max
对应的时刻点t
i,max
，即为拐点出现时间。
49.可以理解，本实施例中，将所述斜率差最大的待测点所在的位置作为拐点所在的
位置，由图2中的去噪信号可以看出，当取得的时刻点从a点逼近拐点时，所取的待测点与a点连线斜率接近零，而该待测点与b点连线斜率的绝对值不断增大；当所取的待测点到达拐点时，拐点与a点的连线接近于平行，即此时第一斜率值k
ai
接近0，而拐点与b点连线的第二斜率值k
bi
达到最大；当所取时刻点经过拐点后，待测点的第一斜率值的绝对值不断增大，第二斜率值的变化较小。
50.因此，利用这一规律，计算a、b两点间所有时刻点的待测点的斜率差，找到最大斜率差n
i,max
与其对应的时刻点t
i,max
，便可获得拐点出现的时间t
i,max
。
51.在一个优选的实施例中，所述步骤s300包括：
52.步骤s310、确定第一阈值，所述第一阈值为管道正常状态下多个检测信号值的平均值；
53.在管道发生泄漏前采集管道正常状态下的多个检测信号值，计算多个检测信号值的平均值，将该均值设置为第一阈值；
54.步骤s320、按时间顺序从信号时程曲线中依次选取一段连续数量的检测信号值，当确定所述连续数量的检测信号值均低于第一阈值时，则将该连续数量的检测信号值中的第一个检测信号值所对应的时间作为第一时间，将第一时间之前的信号时程曲线划分为平稳段；
55.可以理解，若没有连续的、一定数量的检测信号低于第一阈值，则认为该检测信号处于平稳段；当确定所述连续数量的检测信号值均低于第一阈值时，则确定管道发生泄漏；
56.步骤s330、按时间顺序继续从信号时程曲线中依次选取一段连续数量的检测信号值，当确定所述连续数量的检测信号值均低于第二阈值时，则将该连续数量的检测信号值中的最后一个检测信号值所对应的时间作为第二时间，将第二时间之后的信号时程曲线划分为下降段；其中，所述第二阈值小于第一阈值；
57.步骤s340、将第一时间和第二时间之间的信号时程曲线划分为待测段。
58.在本实施例中，根据管道正常运行时检测信号值的波动程度，将第二阈值设置为低于第一阈值的某一值，如果出现连续的、一定数量的信号低于第二阈值，则判定为出现负压波，低于第二阈值的信号为下降段信号；需要说明的是，第二阈值的大小根据实际情况进行设置，一般地，第二阈值相比第一阈值越小，则待测段的时间段越长，需要执行更多的计算来确定拐点时间，但是可以更准确的确定拐点出现时间；因此，第二阈值的大小需要在保证拐点出现时间准确率的前提下进行设定；在一实施例中，第一阈值的大小设定为121.217，相应的，第二阈值的大小设置为100；在平稳段中选择第2500个点作为起始点a，在下降段中选择第二阈值之后的第3000个点作为结束点b。
59.与图1的方法相对应，本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质上存储有基于极限逼近的负压波信号拐点定位程序，所述基于极限逼近的负压波信号拐点定位程序被处理器执行时实现如上述任意一实施例所述的基于极限逼近的负压波信号拐点定位方法的步骤。
60.与图1的方法相对应，本发明实施例还提供一种基于极限逼近的负压波信号拐点定位系统，所述系统包括：
61.至少一个处理器；
62.至少一个存储器，用于存储至少一个程序；
63.当所述至少一个程序被所述至少一个处理器执行，使得所述至少一个处理器实现上述任一实施例所述的基于极限逼近的负压波信号拐点定位方法。
64.上述方法实施例中的内容均适用于本系统实施例中，本系统实施例所具体实现的功能与上述方法实施例相同，并且达到的有益效果与上述方法实施例所达到的有益效果也相同。
65.所述处理器可以是中央处理单元(central
‑
processing
‑
unit，cpu)，还可以是其他通用处理器、数字信号处理器(digital
‑
signal
‑
processor，dsp)、专用集成电路(application
‑
specific
‑
integrated
‑
circuit，asic)、现场可编程门阵列(field
‑
programmable
‑
gate
‑
arr ay，fpga)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等，所述处理器是所述基于极限逼近的负压波信号拐点定位系统的控制中心，利用各种接口和线路连接整个基于极限逼近的负压波信号拐点定位系统可运行装置的各个部分。
66.所述存储器可用于存储所述计算机程序和/或模块，所述处理器通过运行或执行存储在所述存储器内的计算机程序和/或模块，以及调用存储在存储器内的数据，实现所述基于极限逼近的负压波信号拐点定位系统的各种功能。所述存储器可主要包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需的应用程序(比如声音播放功能、图像播放功能等)等；存储数据区可存储根据手机的使用所创建的数据(比如音频数据、电话本等)等。此外，存储器可以包括高速随机存取存储器，还可以包括非易失性存储器，例如硬盘、内存、插接式硬盘，智能存储卡(smart
‑
media
‑
card，smc)，安全数字(secure
‑
digital，sd)卡，闪存卡(flash
‑
card)、至少一个磁盘存储器件、闪存器件、或其他易失性固态存储器件。
67.尽管本技术的描述已经相当详尽且特别对几个所述实施例进行了描述，但其并非旨在局限于任何这些细节或实施例或任何特殊实施例，而是应当将其视作是通过参考所附权利要求，考虑到现有技术为这些权利要求提供广义的可能性解释，从而有效地涵盖本技术的预定范围。此外，上文以发明人可预见的实施例对本技术进行描述，其目的是为了提供有用的描述，而那些目前尚未预见的对本技术的非实质性改动仍可代表本技术的等效改动。