[发明专利]一种微地震事件检测方法及系统

权利要求书

1.一种微地震事件检测方法，其特征在于，包括：

获取利用一阶线性回归函数和克里金插值方法对已标记微地震波动的微震检测序列进行训练得到的时间窗最优长度L；其中L为奇数；

获取按时间顺序采集的待检测的微地震数据序列；

建立第一时间窗和第二时间窗，所述第一时间窗和所述第二时间窗的长度均为时间窗最优长度L；所述第二时间窗紧随所述第一时间窗；

对所述微地震数据序列的左侧和右侧进行数据填充，使所述第一时间窗和所述第二时间窗从所述第一时间窗的中心对准所述微地震数据序列的第一个数据滑动至所述第一时间窗的中心对准所述微地震数据序列的最后一个数据时，所述第一时间窗和所述第二时间窗内均填充满L个数据；

从所述第一时间窗的中心对准填充前的所述微地震数据序列的第一个数据滑动所述第一时间窗和所述第二时间窗，直至所述第一时间窗的中心对准填充前的所述微地震数据序列的最后一个数据，每滑动一个数据的步长，对所述第一时间窗内的数据与所述第二时间窗内的数据进行卡方检验，从而确定所述微地震数据序列中的各个数据处是否发生微地震波动；

记录所有发生微地震波动的数据在所述微地震数据序列中的序号；

所述对所述微地震数据序列的左侧和右侧进行数据填充，具体包括：

将所述微地震数据序列的左侧补充至少个数据；左侧补充的至少个数据均与所述微地震数据序列最左侧的数据相同；

将所述微地震数据序列的右侧补充至少个数据；右侧补充的至少个数据均与所述微地震数据序列最右侧的数据相同。

2.根据权利要求1所述的一种微地震事件检测方法，其特征在于，所述每滑动一个数据的步长，对所述第一时间窗内的数据与所述第二时间窗内的数据进行卡方检验，从而确定所述微地震数据序列中的各个数据处是否发生微地震波动，具体包括：

获取利用一阶线性回归函数和克里金插值方法对已标记微地震波动的微震检测序列进行训练得到的子序列最优个数K；

求取所述第一时间窗内所有数据的最大值和最小值，得到取值区间；

将所述取值区间划分为K段，得到K段小区间；

计算所述第一时间窗内的数据介于每个所述小区间内的数量F⁽ⁿ⁾_P(k)和所述第二时间窗内的数据介于每个所述小区间内的数量F⁽ⁿ⁾_Q(k)；其中n表示位于所述第一时间窗中心的数据在所述微地震数据序列中的序号；k为小区间序号，k＝1,2,…,K；

计算数量差值比率：

根据参数K和置信度α，通过查询卡方表获取卡方值χ²(α,K-1)；

判断所述数量差值比率是否大于或等于所述卡方值，得到第一判断结果；

若所述第一判断结果表示所述数量差值比率大于或等于所述卡方值，则确定在所述微地震数据序列的第n个数据处发生微地震波动；若所述第一判断结果表示所述数量差值比率小于所述卡方值，则确定在所述微地震数据序列的第n个数据处未发生微地震波动。

3.根据权利要求2所述的一种微地震事件检测方法，其特征在于，利用一阶线性回归函数和克里金插值方法对已标记微地震波动的微震检测序列进行训练的过程为：

获取已标记微地震波动的微震检测序列的微震波动的实际位置；

以子序列个数为横坐标，时间窗长度为纵坐标建立二维坐标系，根据时间窗长度的取值范围和子序列个数的取值范围确定所述二维坐标系内的取值区域；

在所述取值区域内随机选取预设数量的点作为插值点，以每个插值点对应的纵坐标作为时间窗建立时第一时间窗和第二时间窗的长度，以每个插值点对应的横坐标作为取值区间的划分个数，计算每个插值点所对应的微震波动预测位置；

根据所述微震波动的实际位置计算每个插值点所对应的微震波动预测位置的正确率，得到每个插值点的预测正确率；

将所述取值区域进行网格划分，获取各个网格点作为被插值点；

利用每个插值点的预测正确率对每个被插值点进行线性拟合，得到拟合值；

求取所有拟合值中的最大值对应的插值点，得到最优插值点；所述最优插值点对应的横坐标即为子序列最优个数，所述最优插值点对应的纵坐标即为时间窗最优长度。