CN113723207B - 一种基于直方图距离的声发射信号突变检测方法 - Google Patents

Abstract

本发明提供一种基于直方图距离的声发射信号突变检测方法，所述包括以下步骤：(1)采集原始声发射信号；(2)分割声发射信号，将原始声发射信号以k点为分界点分为两个窗口段；(3)计算两个窗口段的相对频率直方图；(4)获取两个窗口段的相对频率直方图距离；(5)确定声发射信号突变点，移动两相邻窗口段，搜索直方图距离最大的点，将该点的对应时间视为信号突变时间。本发明能准确、快速识别声发射信号发生突变的时间，解决了现有声发射信号突变检测方法检测精度低、稳定性差问题，适用于故障诊断、地质学等多种领域。

Description

一种基于直方图距离的声发射信号突变检测方法

技术领域

本发明涉及信号检测技术领域，尤其涉及一种基于直方图距离的声发射信号突变检测方法。

背景技术

突变信号中的突变点经常携带重要的信息，是信号重要的特征之一。在机械故障诊断领域，信号突变点往往表示故障的产生。在地质学领域，地震信号在P波、S波到达时会产生突变。然而，在实际环境中，由于背景噪声过大，或者信号本身比较微弱，很难对信号突变时间做出准确的判断。这对信号的处理和后续工作带来了影响。因此，准确识别信号突变时间具有重要意义。

现有的突变信号检测的方法包括基于赤池信息准则(AIC)方法，时窗能量法等。然而，赤池信息方法对低信噪比信号处理效果不理想，并且突变点难以确定，容易产生误判；时窗能量法窗口选择对信号处理影响过大，同样也不适用于低信噪比情况。

发明内容

为了解决不同信噪比环境下突变信号检测精度低、稳定性差等问题，本发明提出一种基于直方图距离的声发射信号突变检测方法。本方法具有效率高、判断准确的特点，克服了窗口选择带来的问题，可应用于多种领域和不同信噪比的信号，例如机械加工监测和地震信号拾取。

本发明提供一种基于直方图距离的声发射信号突变检测方法，所述包括以下步骤：

(1)采集信号：以符合香农采样定理的采样频率来采集时域声发射信号；

(2)分割信号：对含有n个元素的原始离散信号进行分割，将原始信号分割成任意两相邻窗口段，记为interval1和interval2：interval1包含从第1个到第k个元素，interval2包含从第(k+1)个到第n个元素；

(3)计算窗口的相对频率直方图：定义信号时间序列直方图并计算相邻两段信号直方图；

(4)获取两窗口段相对频率直方图距离：使用巴氏系数定义直方图距离，使其满足度量公理，计算相邻两窗口段的直方图距离；

(5)突变点确定：根据d_B(h^A,h^B)的定义，移动两个相邻窗口段interval1和interval2，搜索直方图距离最大的点，并将该点的对应时间视为声发射信号的突变时间点。

步骤(3)中计算窗口的相对频率直方图包括以下步骤：

(3.1)定义直方图：定义为窗口直方图，其中/>

(3.2)计算直方图：根据公式，计算生成两个窗口段的相对频率直方图：其中/>为第i个直方中样本的概率，/>为属于[x_i,x_i+1]区间的元素数量，n是离散信号中的元素数量。

步骤(4)中获取两窗口段相对频率直方图距离包括以下步骤：

(4.1)定义直方图距离：其中B(h^A,h^B)是直方图h^A与直方图h^B之间巴氏系数，记为：/>d_B(h^A,h^B)为直方图h^A与直方图h^B的直方图距离，其满足度量公理，即满足正定性、对称性和三角不等式；

(4.2)直方图距离判定突变点：d_B介于0到1之间，当两个直方图趋近相同时，d_B为趋近0；信号突变点前后两直方图距离最大。

本发明的有益成果为：

(1)本发明改进了算法流程，通过移动窗口段来代替传统的窗口分段方式，有效解决了传统算法由于窗口选择带来的问题，方法更具有实用性。

(2)本发明利用断铅实验和多组地震波的P相识别评估本方法的效果，并与现有技术相比较。测试结果表明，本发明相对于现有技术识别精度更高、稳定性更好、识别速度更快。本发明可以运用于机械加工、地质学等多种领域和不同信噪比下声发射信号突变检测。

附图说明

图1为本发明基于直方图距离的声发射信号突变检测方法的流程示意图。

图2为本发明基于直方图距离的声发射信号突变检测方法的原始信号分割示意图。

图3为本发明基于直方图距离的声发射信号突变检测方法检测误差及分布曲线。

图4为本发明基于直方图距离的声发射信号突变检测方法断铅实验的检测结果。

图5为本发明基于直方图距离的声发射信号突变检测方法的检测结果与人工标记的P波到达时间的偏移量。

图6为本发明基于直方图距离的突变信号检测方法地震波的P相识别实验的检测结果。

图7为现有技术基于AIC的突变信号检测方法误差。

图8为现有技术人工标记断铅实验的检测结果。

图9为现有技术基于AIC的突变信号检测方法断铅实验的检测结果。

图10为现有技术基于AIC的突变信号检测方法的检测结果与人工标记的P波到达时间的偏移量。

图11为现有技术人工标记地震波的P相识别实验的检测结果。

具体实施方式

本发明提供一种基于直方图距离的声发射信号突变检测方法，图1为为本发明基于直方图距离的声发射信号突变检测方法的流程示意图。下面结合图1，描述本发明方法的具体实施方式。

识别信号发生突变时间的难度为：低信噪比信号处理效果不理想；突变点难以确定，容易产生误判；窗口选择对信号处理影响过大。为解决以上问题，本发明提出一种基于直方图距离的声发射信号突变检测方法，用于判断信号发生突变的准确时间，可以识别振动等故障信号、地震信号等到达时间。尤其在低信噪比情况下，本发明依然有着较高的精确度和稳定性，并且可以避免窗口选择带来的问题。

本发明提供一种基于直方图距离的声发射信号突变检测方法，所述包括如下步骤：

(1)采集信号：以符合香农采样定理的采样频率采集的时域信号为对象；

(2)分割信号：如图2所示对原始离散信号进分割，将原始信号分割成任意两相邻窗口段，记为interval1和interval2。interval1包含从第1个到第k个元素，interval2包含从第(k+1)个到第n个元素。通过此分段方式，避免了窗口长度等参数的影响；

(3)计算窗口的相对频率直方图：定义信号时间序列直方图并计算相邻两段信号直方图；

(4)获取两窗口段相对频率直方图距离：使用巴氏系数定义直方图距离，使其满足度量公理，计算相邻两窗口段的直方图距离；

(6)突变点确定：根据d_B(h^A,h^B)的定义，移动两个相邻窗口段interval1和interval2，搜索直方图距离最大的点，并将该点的对应时间视为信号突变时间。

步骤(3)中计算窗口的相对频率直方图包括以下步骤：

(3.1)定义直方图：定义为窗口直方图，其中/>

(3.2)计算直方图：根据公式：计算生成两个窗口段的相对频率直方图：其中/>为第i个直方中样本的概率，/>为属于[x_i,x_i+1]区间的元素数量，n是离散时间序列中的元素数量。

步骤(4)中获取两窗口段相对频率直方图距离包括以下步骤：

(4.1)定义直方图距离：其中B(h^A,h^B)是直方图h^A与直方图h^B之间巴氏系数，记为：/>d_B(h^A,h^B)为直方图h^A与直方图h^B的直方图距离，其满足度量公理，即满足正定性、对称性和三角不等式；

(4.2)直方图距离判定突变点：d_B介于0到1之间，当两个直方图趋近相同时，d_B为趋近0；信号突变点前后两直方图距离最大。

本发明提出的一种基于直方图距离的声发射信号突变检测方法可以准确的识别突变点，不容易误判。实验验证如下：

(1)断铅实验验证：铅芯断裂时会产生突变的声发射信号，通过与人工检测结果比较(人工检测结果，定义为分析人员在其常规分析中确定的结果，经经验认为是正确的)，从而验证本发明所述的方法识别不同信噪比信号突变时间的精度。将直径0.3mm，长约30mm的铅笔芯固定在混凝土梁的表面。使铅笔芯断裂在三个不同点。铅芯断裂释放的能量由10个工作频率范围为200-850kHz的压电式声发射传感器收集。原始信号使用26dB增益放大。采样频率为1MHz，阈值为31dB。铅芯断裂测试生成了100组声发射信号。经过上述5个步骤，得到100组声信号的突变时间点。将其与人工检测的突变时间点比较。

自动检测结果与手动检测结果的误差如图3所示。信噪比是声发射信号与噪声信号的比值，信噪比越小，噪声相对越大。密度曲线为相应的经验分布，提供了对均值的估计。偏移量为自动检测结果的误差，均值为偏移量数据集的平均大小，标准差用来衡量偏移量数据点的离散程度，均值和标准差越小，误差越小，检测结果越稳定，越精确。当信噪比大于10dB时，误差值大多集中在均值附近。当信噪比低于10dB时，自动检测结果和手动检测结果之间的分散程度高。高信噪比组(信噪比大于10dB)中误差小于20μs的百分比为94％。低信噪比组(信噪比小于10dB)中小于20μs的误差百分比为76.5％。

基于赤池信息准则(AIC)的算法是目前常用的突变信号自动识别算法。将本发明所述的检测方法与AIC方法作比较。如图7所示，AIC方法在高信噪比(信噪比大于10dB)时结果与本发明所述的方法相似，但在低信噪比(信噪比小于10dB)时结果更为分散，表明在低信噪比情况下赤池信息方法不容易准确识别低信噪比信号。使用本发明的检测方法和AIC方法的检测结果如图4、图8和图9所示：图8为人工标记断铅实验的检测结果。信号在时间为0.1920秒时发生突变，图4为本发明所示方法对断铅实验的检测结果，信号在时间为0.1920秒时发生突变，与人工标记的结果一样。图9为AIC方法对对断铅实验的检测结果，AIC值最小点的对应时间视为AIC方法检测到的信号突变时间，可以看出在手动选取点处(0.1920秒)的AIC值发生明显下降，但此后AIC值开始波动，最小值点出现在手动选取点的后面，信号在0.2580秒处发生突变。AIC方法生成的突变时间产生了延迟。相比之下，本发明所述的方法准确地反映信号突变，突变点为直方图距离的最大值。

(2)地震波的P相识别验证：地震信号在P波到达时会发生突变，通过与地震台站测得的P波到达时间比较即可验证本发明所述方法识别信号突变时间的效果。由于地震信号垂直分量上P相最强，这里使用地震信号的垂直分量数据。地震信号从IRIS(IncorporatedResearch Institutions for Seismology，美国地震学研究联合会)地震台站获取，采样频率为100Hz。选择300组不同波形特性的信号，其涵盖了广泛的信噪比。通过本发明所述的方法步骤计算得到信号突变时间。

图5显示了本发明所示方法的检测结果与人工标记的P波到达时间的时间差(偏移量)。图10显示了基于AIC方法的检测结果与人工标记的P波到达时间的时间差。可以看出本发明所述方法的结果偏移分布比AIC方法的要窄得多，即本发明所述的方法误差更小。对于本发明所述的方法，±2秒内的偏移量数占数据集的45％。对于AIC方法，落在同一范围内偏移数的百分比为39.7％。相比之下，本发明所述的方法有更高的精度。图11为人工标记人地震波的P相识别实验的检测结果，在时间为21:39:26时P波到达，信号发生突变。图6为本发明所示方法P相识别实验的检测结果，在时间为21:39:26时P波到达，信号发生突变，这与人工标记结果一样。Z相为垂直分量上的地震信号。对于不同的波形特性，本发明所述的方法可以成功地检测到P的到达，第一个峰值即本发明所述方法认为的突变时间，这与人工标记的P波到达时间精确一致。

上述的两个验证性实验中，断铅实验是最常用的声发射相关检测和验证实验，而声发射技术广泛运用于无损检测领域。地震波P相识别是地震信号后处理，震源定位的基础，在地质学领域有着重要意义。本发明所述的方法历经多次实验，可信度高；本发明所述的方法在与现有常用的AIC方法比较中，展现了精度高的特点。本发明所述的方法历经广泛分布的信噪比信号测试，尤其对于低信噪比声发射信号突变的识别，本发明所述的方法具有准确率高的特点，适用范围广。综上，本发明所述的一种基于直方图距离的声发射信号突变检测方法具有很高的实用性，可广泛地运用于机械加工、地质学等领域。

本说明书实施例所述的内容仅仅是对发明构思的实现形式的列举，本发明的保护范围不应当被视为仅限于实施例所陈述的具体形式，本发明的保护范围也及于本领域技术人员根据本发明构思所能够想到的等同技术手段。

Claims (1)

1.一种基于直方图距离的声发射信号突变检测方法，其特征在于，包括以下步骤：

(1)采集信号：以符合香农采样定理的采样频率来采集时域声发射信号；

(2)分割信号：对含有n个元素的原始离散信号进行分割，将原始信号分割成任意两相邻窗口段，记为interval1和interval2：interval1包含从第1个到第k个元素，interval2包含从第(k+1)个到第n个元素；

(3)计算窗口的相对频率直方图：定义信号时间序列直方图并计算相邻两段信号直方图。包括以下步骤：

(3.1)定义直方图：定义为窗口直方图，其中/>

(3.2)计算直方图：根据公式，计算生成两个窗口段的相对频率直方图：其中为第i个直方中样本的概率，/>为属于[x_i,x_i+1]区间的元素数量，n是离散信号中的元素数量；

(4)获取两窗口段相对频率直方图距离：使用巴氏系数定义直方图距离，使其满足度量公理，计算相邻两窗口段的直方图距离。包括以下步骤：

(4.1)定义直方图距离：其中B(h^A,h^B)是直方图h^A与直方图h^B之间巴氏系数，记为：/>d_B(h^A,h^B)为直方图h^A与直方图h^B的直方图距离，其满足度量公理，即满足正定性、对称性和三角不等式；

(4.2)通过直方图距离判定突变点：d_B介于0到1之间，当两个直方图趋近相同时，d_B为趋近0；信号突变点前后两直方图距离最大；

(5)突变点确定：根据d_B(h^A,h^B)的定义，移动两个相邻窗口段interval1和interval2，搜索直方图距离最大的点，并将该点的对应时间视为声发射信号的突变时间点。