[发明专利]一种利用边缘探测的微震事件检测方法和系统

权利要求书

1.一种利用边缘探测的微震事件检测方法，其特征在于，包括：

步骤101获取按时间顺序采集的信号序列S；

步骤102对信号序列S进行移动平均滤波，具体为：所用滤波公式为：

y_n＝avg[s_n,M₀]

y_n为移动平均滤波后的信号序列Y的第n个元素，

avg[s_n,M₀]表示对s_n进行移动平均滤波，

s_n为所述信号序列S的第n个元素，

n＝1,2,…,N为元素下标，

N为所述信号序列S的长度，

为移动平均滤波参数，

SNR为所述信号序列S的信噪比，

对*下取整，

*表示任一自变量；

步骤103求取N个M₁点最大值和N个M₂点最大值，具体为：第k个M₁点最大值记为第k个M₂点最大值记为所用求取公式为：

其中：

为移动平均滤波之后的信号序列Y的第k-M₁个元素，

为移动平均滤波之后的信号序列Y的第k-M₁+1个元素，

为移动平均滤波之后的信号序列Y的第k+M₁-1个元素，

为移动平均滤波之后的信号序列Y的第k+M₁个元素，

为移动平均滤波之后的信号序列Y的第k-M₂个元素，

为移动平均滤波之后的信号序列Y的第k-M₂+1个元素，

为移动平均滤波之后的信号序列Y的第k+M₂-1个元素，

为移动平均滤波之后的信号序列Y的第k+M₂个元素，

为第一窗口长度，

为第二窗口长度，

k＝1,2,…,N为窗口序号；

步骤104求取N个边缘探测值，具体为：第k个边缘探测值记为H_k，所用求取公式为：

其中：

α为边缘探测因子，所用求取公式为：

步骤105求取微震事件判断阈值，具体为：微震事件判断阈值记为ε₀，所用求取公式为：

其中：

σ为信号序列S的均方差；

步骤106检测微震事件，具体为：如果第k个窗口的边缘探测值H_k满足判断条件|H_k|≥ε₀，则在所述信号序列S的第k点处，检测到微震事件；否则未检测到微震事件。

2.一种利用边缘探测的微震事件检测系统，其特征在于，包括：

模块201获取按时间顺序采集的信号序列S；

模块202对信号序列S进行移动平均滤波，具体为：所用滤波公式为：

y_n＝avg[s_n,M₀]

y_n为移动平均滤波后的信号序列Y的第n个元素，

avg[s_n,M₀]表示对s_n进行移动平均滤波，

s_n为所述信号序列S的第n个元素，

n＝1,2,…,N为元素下标，

N为所述信号序列S的长度，

为移动平均滤波参数，

SNR为所述信号序列S的信噪比，

对*下取整，

*表示任一自变量；

模块203求取N个M₁点最大值和N个M₂点最大值，具体为：第k个M₁点最大值记为第k个M₂点最大值记为所用求取公式为：

其中：

为移动平均滤波之后的信号序列Y的第k-M₁个元素，

为移动平均滤波之后的信号序列Y的第k-M₁+1个元素，

为移动平均滤波之后的信号序列Y的第k+M₁-1个元素，

为移动平均滤波之后的信号序列Y的第k+M₁个元素，

为移动平均滤波之后的信号序列Y的第k-M₂个元素，

为移动平均滤波之后的信号序列Y的第k-M₂+1个元素，

为移动平均滤波之后的信号序列Y的第k+M₂-1个元素，

为移动平均滤波之后的信号序列Y的第k+M₂个元素，

为第一窗口长度，

为第二窗口长度，

k＝1,2,…,N为窗口序号；

模块204求取N个边缘探测值，具体为：第k个边缘探测值记为H_k，所用求取公式为：

其中：

α为边缘探测因子，所用求取公式为：

模块205求取微震事件判断阈值，具体为：微震事件判断阈值记为ε₀，所用求取公式为：

其中：

σ为信号序列S的均方差；

模块206检测微震事件，具体为：如果第k个窗口的边缘探测值H_k满足判断条件|H_k|≥ε₀，则在所述信号序列S的第k点处，检测到微震事件；否则未检测到微震事件。