[发明专利]一种冲击矿压的微震多维信息综合时序预警方法

权利要求书

1.一种冲击矿压的微震多维信息综合时序预警方法，其特征在于：包括如下步骤：

（1）对矿震数据进行分区分级筛选；

（2）建立微震多维信息指标体系，包括频次时序指标W₁、空间时序指标W₂和强度时序指标W₃，所述频次时序指标W₁包括频次指标W₁₁，所述空间时序指标W₂包括震源集中程度指标W₂₁，所述强度时序指标W₃包括最大应力当量指标W₃₁和平均应力当量指标W₃₂；

（3）根据下式计算微震多维信息综合时序预警指标值W：

W=W₁+W₂+W₃，W₁=ω₁₁·W₁₁，

W₂=ω₂₁·W₂₁，W₃=ω₃₁·W₃₁+ω₃₂·W₃₂

其中，ω₁₁,ω₂₂,ω₃₁,ω₃₂为相应指标的预测信度，且满足ω₁₁+ω₂₂+ω₃₁+ω₃₂=1，λ_ij(t)为相应指标在监测时段t内的异常隶属度，取值范围为0-1；

（4）根据微震多维信息综合时序预警指标值W确定危险等级和对应危险等级的防治对策；微震多维信息综合时序预警指标值W越大危险等级越高。

2.根据权利要求1所述的冲击矿压的微震多维信息综合时序预警方法，其特征在于：所述步骤（1）中，采用任意多边形分区方法对矿震数据进行分区筛选。

3.根据权利要求1所述的冲击矿压的微震多维信息综合时序预警方法，其特征在于：所述步骤（1）中，根据矿震能量对矿震数据进行分级筛选，即矿震能量上、下限分别以Gurtenberg-Richter幂率曲线高能量端和低能量端偏离幂率的拐点作为能量上、下限的分界线，所述拐点根据下式计算识别：

能量上限分界线根据下式计算识别：

Ri-High=Σj=i-1n(xj-x‾)(yj-y‾)Σj=i-1n(xj-x‾)2Σj=i-1n(yj-y‾)2]]>

其中，i=2,3,…,n，R_i-High为Gurtenberg-Richter幂率曲线中横坐标lgE_i对应的相关系数值，参与计算的x序列为lgE_i-1,lgE_i,lgE_i+1,…,lgE_n，y序列为lgN_i-1,lgN_i,lgN_i+1,…,lgN_n，为x序列的样本均值，为y序列样本均值；

当Rk-High=min{Rn-12,Rn+12,···,Rk-High,···,Rn}]]>时，能量上限的分界线为Gurtenberg-Richter幂率曲线中纵坐标R_k-High对应横坐标值lgE_k中的能量E_k；

能量下限分界线根据下式计算识别：

Ri-Low=Σj=1i(xj-x‾)(yj-y‾)Σj=1i(xj-x‾)2Σj=1i(yj-y‾)2]]>

其中，i=2,3,…,n，R_i-Low为Gurtenberg-Richter幂率曲线中横坐标lgE_i对应的相关系数值，参与计算的x序列为lgE₁,…,lgE_i-1,lgE_i，y序列为lgN₁,…,lgN_i-1,lgN_i，为x序列的样本均值，为y序列样本均值；

当时，能量下限的分界线为Gurtenberg-Richter幂率曲线中纵坐标R_k-Low对应横坐标值lgE_k中的能量E_k。