[发明专利]一种基于震源参数的微震与爆破事件识别方法

说明书

技术领域

本发明涉及一种基于震源参数的微震与爆破事件识别方法，结果明确，适用性强，准确率高。

背景技术

微震监测作为一种有效的地压监测手段在国内外得到了广泛应用，而微震与爆破事件的区分识别是数据处理过程中的关键问题，因此对微震与爆破事件识别方法的研究具有重要意义。但微震与爆破事件波形均较为复杂，爆破事件往往会造成重复波形，其中被分离的单个波形可能会人为的被错误处理成微震事件，造成严重后果。目前主要根据一天内的爆破时间和波形重复来手动区分微震与爆破事件，但易受操作人员专业知识和经验的影响，且识别数量有限，降低了微震监测实时分析的效果。

目前针对微震与爆破事件识别的方法可分为：波形频谱分析法、多参数统计法和软计算法，这些方法通常包括特征提取和特征识别两个过程。常用的特征提取包括震级、能量、视应力、视体积、静态应力降、动态应力降和时域波形特征参数(振幅、频率)等，特征识别方法包括fisher分类法、Logistic回归法、随机森林法、神经网络法、支持向量机法和贝叶斯分类法等。常用的特征参数大多难以自动得到，波形时域特征参数是在单一尺度上分析得到的，其信息量较少，限制了微震与爆破事件的自动识别。

对于分析微震与爆破事件之间波形和波谱的不同特征，波形频谱分析法不失为一种好方法，其主要借助小波分析、小波包分析和频率切片小波分析。唐守锋等(2011)根据采样定理和Mallat算法确定了小波分解最大尺度，并提出采用小波特征能谱系数作为煤岩破裂微震信号识别的定量表征；朱权洁等(2012a)运用小波包分析对微震信号进行了5层多尺度分解，并与爆破震动信号的能量分布进行了对比分析；朱权洁等(2012b)结合小波分析与分形理论对爆破振动、岩石破裂及电磁干扰3类信号进行了5层小波包分解，并以筛选后的23个小波分形盒维数作为支持向量机识别的特征向量；赵国彦等(2015)采用频率切片小波变换对矿山岩体破裂信号和爆破振动信号特定频带能量比和相关系数进行了研究。小波分析、小波包分析和频率切片小波分析具有较好的自适应性，但很容易受到信号中相邻谐波成分的交叠影响，致使不同频带信号存在混叠。

可见现有的微震事件与爆破事件识别方法存在较大的局限，需要研究一种适用性强、准确率高的自动识别方法。

发明内容

本发明所要解决的技术问题是，针对现有技术的不足，提供一种基于震源参数的微震与爆破事件识别方法，该微震与爆破事件识别方法适用性强、准确率高。

本发明的技术解决方案如下：

一种基于震源参数的微震与爆破事件识别方法，包括以下步骤：

步骤1：获得微震事件与爆破事件的判别式：

以已确定的N组微震事件与L组爆破事件为样本，获得微震与爆破事件的判别式：

其中Y_E为微震事件判别式的因变量，Y_B为爆破事件判别式的因变量，式中，X_i，i＝1,2,…,6为6个与典型震源参数相关的值，C₀为微震事件判别式中的常数项，C_i为在微震事件判别式中对应于X_i的最优系数，B₀爆破事件判别式中的常数项，B_i为在爆破事件判别式中对应于X_i的最优系数；C₀、C_i、B₀和B_i基于样本辨识得到；N和L为整数，且N，L≥10；

步骤2：对待识别事件进行识别：

计算待识别事件对应的X_i取值，代入所述的微震事件与爆破事件的判别式，计算得到Y_E和Y_B；若Y_E大于或等于Y_B，则说明该待识别事件为微震事件，否则若Y_E小于Y_B，则说明该待识别事件为爆破事件。

所述步骤1中，针对每一组微震事件或爆破事件，进行以下操作：

(1)通过微震监测系统获取震源参数；

包括：S波与P波能量之比值E_S/E_P，地震力矩M，静应力降SSD，动应力降DSD，震源发震时间t以及力矩震级M_m；