微震监测冲击地压预警指标优化方法.pdf

《微震监测冲击地压预警指标优化方法.pdf》由会员分享，可在线阅读，更多相关《微震监测冲击地压预警指标优化方法.pdf（19页完成版）》请在专利查询网上搜索。

1、(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 202010269897.3 (22)申请日 2020.04.08 (71)申请人 中煤能源研究院有限责任公司 地址 710054 陕西省西安市碑林区雁塔路 北段66号中煤能源大厦二层 申请人 乌审旗蒙大矿业有限责任公司 (72)发明人 韩刚郑忠友赵忠证张寅 王朝引赵雷王黔张振配 吕玉磊陈国华白俊杰曹京龙 解嘉豪靳俊晓刘洋张翔 宁石杰石义恒杨振周林 刘成勇孙俊彦李娟 (74)专利代理机构 西安通大专利代理有限责任 公司 61200 代理人 闵岳峰 (51)Int.Cl. G06Q。

2、 10/04(2012.01) G06Q 10/06(2012.01) (54)发明名称 一种微震监测冲击地压预警指标优化方法 (57)摘要 本发明公开了一种微震监测冲击地压预警 指标优化方法， 该方法以冲击地压矿井实际动力 显现、 强矿震事件的长期观测数据为基础， 首先 建立单一预警指标确定准则， 对单一预警指标进 行预警能效评定， 明确其各自预警效能， 再综合 各项预警指标的预警效能， 形成多参数联合预警 评价方法， 实现微震监测的多参数联合预警， 提 高微震系统对冲击地压预警的准确性。 本发明主 要适用于冲击地压矿井微震监测系统的监测预 警， 特别是微震监测预警指标优化方面。 权利要求书。

3、4页 说明书12页 附图2页 CN 111476420 A 2020.07.31 CN 111476420 A 1.一种微震监测冲击地压预警指标优化方法， 其特征在于， 包括以下步骤： (1)收集数据 收集冲击地压矿井工作面采掘期间微震监测数据、 所有动力显现以及矿震事件相关微 震监测数据； (2)选择、 拓展单一预警指标种类 初期借鉴本矿井其他采掘工作面、 邻近矿井、 行业常规单一预警指标种类， 在日常监测 预警工作中不断发掘新的预警指标； 常规单一预警指标种类选择： 单一大能量事件、 104J以 上事件频次、 每日微震事件总能量、 每日微震事件总频次以及每刀能量； (3)对每项单一预警指标。

4、在工作面采掘期间微震监测数据进行倒序排列 对目标工作面采掘期间监测到的微震事件按每项单一预警指标进行统计计算， 并将计 算结果进行倒序排列， 按横坐标为时间轴绘制倒序排列图表； (4)在倒序排列图中标识历次动力显现、 矿震事件 将收集到的本工作面历次动力显现、 矿震事件按实际发生时间表示在倒序排列图表 中； (5)按以下原则确定每项单一预警指标具体预警值 设定强冲击危险预警指标能够预警80冲击事件， 在倒序排列图表中在选定单一预 警指标具体预警值|kqi|； 设定中等冲击危险预警指标能够预警100冲击事件， 在倒序排列图表中在选定单 一预警指标具体预警值|kzi|； 设定弱等冲击危险预警指标能。

5、够预警100强震动事件， 在倒序排列图表中在选定 单一预警指标具体预警值|kri|； 设定无冲击危险预警指标为|kri|； 其中， i为第i个单一预警指标种类， q、 z、 r、 w分别代表强、 中等、 弱、 无冲击危险等级； (6)计算每项单一预警指标预警效能评价指标 预警指标预警效能评价指标包括准确率、 漏报率和误报率； (7)有效预警指标筛选 依据步骤(6)中计算结果， 排查各项单一预警指标预警效能指标是否同时具备准确率 最小、 漏报率最大、 误报率最大， 如满足以上条件， 说明该项单一预警指标预警效能差， 不适 合作为一项预警指标， 应删去， 否则保留； (8)计算每项单一预警指标预警。

6、权重 计算每项单一预警指标实际预警率 其中： 权利要求书 1/4 页 2 CN 111476420 A 2 i0,1,2,n； Eiq为第i项单一预警指标种类强冲击危险判定指标的实际预警率； Diq为第i项单一预警指标种类采用强冲击危险判定指标时在工作面采掘期间总的预警 天数； Eiz为第i项单一预警指标种类中等冲击危险判定指标的实际预警率； Diz为第i项单一预警指标种类采用中等冲击危险判定指标时在工作面采掘期间总的预 警天数； Eir为第i项单一预警指标种类弱冲击危险判定指标的实际预警率； Dir为第i项单一预警指标种类采用弱冲击危险判定指标时在工作面采掘期间总的预警 天数； S为采掘工作。

7、面全部生产天数； 选择每项单一预警指标实际预警率最小值 EminqMIN(Eiq) EminzMIN(Eiz) EminrMIN(Eir) 计算每项单一预警指标预警权重 其中： iq为第i项单一预警指标种类采用强冲击危险判定指标时预警权重； iz为第i项单一预警指标种类采用中等冲击危险判定指标时预警权重； ir为第i项单一预警指标种类采用弱冲击危险判定指标时预警权重； (9)计算当日综合动载预警值 计算当日单一预警指标的预警权重 若ki|kqi|， 表示当日第i项单一预警指标实际监测值达到强冲击危险等级， 则取 i iq； 若|kqi|ki|kzi|， 表示当日第i项单一预警指标实际监测值达到。

8、中等冲击危险等级， 则取 i iz； 若|kzi|ki|kri|， 表示当日第i项单一预警指标实际监测值达到弱冲击危险等级， 则 取 i ir； 若ki|kri|， 表示当日第i项单一预警指标实际监测值为无冲击危险等级， 则取 i0； 计算当日动载综合预警值 权利要求书 2/4 页 3 CN 111476420 A 3 (10)每日进行综合预警 按以下原则判定当日动载扰动条件下冲击危险程度： 若则判定为当日强冲击危险； 若则判定为当日中等冲击危险； 若则判定为当日弱冲击危险； 若| |0， 则判定为当日无冲击危险； (11)当日是否出现新的动力显现、 矿震事件 如当日出现新的动力显现、 矿震事。

9、件， 则将当日微震监测数据统计后返回步骤(3)， 重 新优化预警值， 直至工作面安全回采结束。 2.根据权利要求1所述的一种微震监测冲击地压预警指标优化方法， 其特征在于， 步骤 (6)中， 每项单一预警指标具体预警值的效能指标计算方法如下： 准确率： 其中： Yiq为第i项单一预警指标种类强冲击危险判定指标的预警准确率； Aiq为采用kqi作为第i项单一预警指标种类的强冲击危险判定指标时， 能够预测到的历 史冲击显现事件次数； Yiz为第i项单一预警指标种类中等冲击危险判定指标的预警准确率； Aiz为采用kqi作为第i项单一预警指标种类的强冲击危险判定指标时， 能够预测到的历 史冲击显现事件。

10、次数； Yir为第i项单一预警指标种类弱冲击危险判定指标的预警准确率； Air为采用kqi作为第i项单一预警指标种类的弱冲击危险判定指标时， 能够预测到的历 史冲击显现事件次数； B为步骤(1)中采掘工作面收集到的动力显现事件数量； 漏报率： 漏报率Liq1-Yiq 漏报率Liz1-Yiz 权利要求书 3/4 页 4 CN 111476420 A 4 漏报率Lir1-Yir Liq为第i项单一预警指标种类强冲击危险判定指标的预警漏报率； Liz为第i项单一预警指标种类中等冲击危险判定指标的预警漏报率； Lir为第i项单一预警指标种类弱冲击危险判定指标的预警漏报率； 误报率： 其中： Xiq为第。

11、i项单一预警指标种类强冲击危险判定指标的预警误报率； Ciq为采用kqi作为第i项单一预警指标种类的强冲击危险判定指标时， 工作面采掘期间 具备强冲击危险的总次数； Xiz为第i项单一预警指标种类中等冲击危险判定指标的预警误报率； Ciz为采用kqi作为第i项单一预警指标种类的中等冲击危险判定指标时， 工作面采掘期 间具备中等冲击危险的总次数； Xir为第i项单一预警指标种类弱冲击危险判定指标的预警误报率； Cir为采用kqi作为第i项单一预警指标种类的弱冲击危险判定指标时， 工作面采掘期间 具备弱冲击危险的总次数。 权利要求书 4/4 页 5 CN 111476420 A 5 一种微震监测冲。

12、击地压预警指标优化方法 技术领域 0001 本发明涉及一种微震监测冲击地压预警指标优化方法， 主要用于有冲击地压灾害 矿井实际生产过程中对冲击地压危险程度的预测预报， 用于指导煤矿的防冲工作。 背景技术 0002 冲击地压是煤矿深部开采面临的主要灾害之一， 对煤矿安全高效生产及煤矿工人 人身安全构成巨大威胁。 目前， 国内关于冲击地压发生机理、 监测预警、 防治措施等方面相 关研究不断深入， 但在实际生产过程中， 现场冲击地压发生可能性判别与预警、 预警指标优 化等方面尚需进一步研究。 0003 目前， 国内冲击地压矿井常用的冲击危险性预警方法主要有钻屑法、 煤岩体变形 测量法、 煤体应力测量。

13、法、 微震监测法等。 其中， 微震监测作为区域监测手段在冲击地压矿 井中应用广泛。 其主要用于监测覆岩裂隙发育程度， 微震事件是覆岩发生裂隙、 破断过程 中， 所集聚的弹性能瞬间释放过程中的一种物理效应。 利用微震监测系统可以较为准确的 监测和计算包括微震事件发生时间、 位置和强度的多种微震活动参数， 为实现矿山微震活 动时空变化规律进行分析提供了可能性， 采用微震活动参数来评价与预测冲击地压， 可以 是冲击地压预测方法与指标建立在一定的理论基础上， 可以提高冲击地压预测的准确性， 目前生产过程中基于微震监测系统的冲击地压预警技术为微震能量分析预测法， 即依据实 际监测到的当日微震能量最大值、。

14、 一定推进距离内释放的微震能量综合作为冲击危险性预 警指标， 当其实际监测值大于阈值时， 判定具备冲击地压发生的可能性， 根据危险程度实施 相应的解危措施。 0004 目前国内冲击地压矿井， 微震能量预警指标通常包括104J以上事件频次、 每日微 震事件总能量、 每刀能量、 每日事件总频次、 单位推进度总能量、 单位推进度总频次等， 实际 使用过程中每日将实测微震数据汇总统计， 形成监测日报， 用于冲击危险性判别。 0005 但在实际监测工作中发现， 目前的微震监测任何单一预警指标均无法实现准确预 警， 需要将各项指标综合起来实现多参数综合预警； 其次， 目前对每项预警指标实际预警效 能评价不。

15、足， 以104J以上事件频次为例， 设定每日104J以上事件为10个时判定为强冲击危 险， 那如果当日出现10个104J以上事件， 其预报准确性如何？ 有无漏报可能？ 有无误报的可 能？ 能否准确反映当日冲击危险的真实水平？ 预警指标如何进一步优化。 此类疑问为目前现 场微震监测预警工作中继续攻克难题， 也是目前冲击地压监测预警行业面临的共性难题。 0006 此外， 目前国内外对于冲击地压发生机理认识尚处于研究过程中， 对于预警前兆 信息识别及预警方法的研究尚处于摸索阶段， 不同的冲击地压矿井实际条件千差万别， 监 测预警指标需要结合各矿实际条件采用 “一矿一策” 甚至是 “一面一策” ， 如。

16、何结合各矿实际 情况形成相对普遍使用的预警指标优化方法， 是目前煤矿冲击地压防治重要研究课题之 一。 说明书 1/12 页 6 CN 111476420 A 6 发明内容 0007 本发明针对以上问题， 提供了一种微震监测冲击地压预警指标优化方法， 该方法 以冲击地压矿井实际动力显现、 强矿震事件的长期观测数据为基础， 首先建立单一预警指 标确定准则， 对单一预警指标进行预警能效评定， 明确其各自预警效能， 再综合各项预警指 标的预警效能， 形成多参数联合预警评价方法， 实现微震监测的多参数联合预警， 提高微震 系统对冲击地压预警的准确性。 本发明主要适用于冲击地压矿井微震监测系统的监测预 警。

17、， 特别是微震监测预警指标优化方面。 0008 本发明采用如下技术方案来实现的： 0009 一种微震监测冲击地压预警指标优化方法， 包括以下步骤： 0010 (1)收集数据 0011 收集冲击地压矿井工作面采掘期间微震监测数据、 所有动力显现以及矿震事件相 关微震监测数据； 0012 (2)选择、 拓展单一预警指标种类 0013 初期借鉴本矿井其他采掘工作面、 邻近矿井、 行业常规单一预警指标种类， 在日常 监测预警工作中不断发掘新的预警指标； 常规单一预警指标种类选择： 单一大能量事件、 104J以上事件频次、 每日微震事件总能量、 每日微震事件总频次以及每刀能量； 0014 (3)对每项单。

18、一预警指标在工作面采掘期间微震监测数据进行倒序排列 0015 对目标工作面采掘期间监测到的微震事件按每项单一预警指标进行统计计算， 并 将计算结果进行倒序排列， 按横坐标为时间轴绘制倒序排列图表； 0016 (4)在倒序排列图中标识历次动力显现、 矿震事件 0017 将收集到的本工作面历次动力显现、 矿震事件按实际发生时间表示在倒序排列图 表中； 0018 (5)按以下原则确定每项单一预警指标具体预警值 0019 设定强冲击危险预警指标能够预警80冲击事件， 在倒序排列图表中在选定单 一预警指标具体预警值|kqi|； 0020 设定中等冲击危险预警指标能够预警100冲击事件， 在倒序排列图表中。

19、在选 定单一预警指标具体预警值|kzi|； 0021 设定弱等冲击危险预警指标能够预警100强震动事件， 在倒序排列图表中在 选定单一预警指标具体预警值|kri|； 0022 设定无冲击危险预警指标为|kri|； 0023 其中， i为第i个单一预警指标种类， q、 z、 r、 w分别代表强、 中等、 弱、 无冲击危险等 级； 0024 (6)计算每项单一预警指标预警效能评价指标 0025 预警指标预警效能评价指标包括准确率、 漏报率和误报率； 0026 (7)有效预警指标筛选 0027 依据步骤(6)中计算结果， 排查各项单一预警指标预警效能指标是否同时具备准 确率最小、 漏报率最大、 误报。

20、率最大， 如满足以上条件， 说明该项单一预警指标预警效能差， 不适合作为一项预警指标， 应删去， 否则保留； 0028 (8)计算每项单一预警指标预警权重 说明书 2/12 页 7 CN 111476420 A 7 0029 计算每项单一预警指标实际预警率 0030 0031 0032 0033 其中： 0034 i0,1,2,n； 0035 Eiq为第i项单一预警指标种类强冲击危险判定指标的实际预警率； 0036 Diq为第i项单一预警指标种类采用强冲击危险判定指标时在工作面采掘期间总的 预警天数； 0037 Eiz为第i项单一预警指标种类中等冲击危险判定指标的实际预警率； 0038 Diz。

21、为第i项单一预警指标种类采用中等冲击危险判定指标时在工作面采掘期间总 的预警天数； 0039 Eir为第i项单一预警指标种类弱冲击危险判定指标的实际预警率； 0040 Dir为第i项单一预警指标种类采用弱冲击危险判定指标时在工作面采掘期间总的 预警天数； 0041 S为采掘工作面全部生产天数； 0042 选择每项单一预警指标实际预警率最小值 0043 EminqMIN(Eiq) 0044 EminzMIN(Eiz) 0045 EminrMIN(Eir) 0046 计算每项单一预警指标预警权重 0047 0048 0049 0050 其中： 0051 iq为第i项单一预警指标种类采用强冲击危险判。

22、定指标时预警权重； 0052 iz为第i项单一预警指标种类采用中等冲击危险判定指标时预警权重； 0053 ir为第i项单一预警指标种类采用弱冲击危险判定指标时预警权重； 0054 (9)计算当日综合动载预警值 说明书 3/12 页 8 CN 111476420 A 8 0055 计算当日单一预警指标的预警权重 0056 若ki|kqi|， 表示当日第i项单一预警指标实际监测值达到强冲击危险等级， 则取 i iq； 0057 若|kqi|ki|kzi|， 表示当日第i项单一预警指标实际监测值达到中等冲击危险 等级， 则取 i iz； 0058 若|kzi|ki|kri|， 表示当日第i项单一预警。

23、指标实际监测值达到弱冲击危险等 级， 则取 i ir； 0059 若ki|kri|， 表示当日第i项单一预警指标实际监测值为无冲击危险等级， 则取 i 0； 0060 计算当日动载综合预警值 0061 0062 (10)每日进行综合预警 0063 按以下原则判定当日动载扰动条件下冲击危险程度： 0064若则判定为当日强冲击危险； 0065若则判定为当日中等冲击危险； 0066若则判定为当日弱冲击危险； 0067 若| |0， 则判定为当日无冲击危险； 0068 (11)当日是否出现新的动力显现、 矿震事件 0069 如当日出现新的动力显现、 矿震事件， 则将当日微震监测数据统计后返回步骤 (3。

24、)， 重新优化预警值， 直至工作面安全回采结束。 0070 本发明进一步的改进在于， 步骤(6)中， 每项单一预警指标具体预警值的效能指标 计算方法如下： 0071 准确率： 0072 0073 0074 0075 其中： 0076 Yiq为第i项单一预警指标种类强冲击危险判定指标的预警准确率； 0077 Aiq为采用kqi作为第i项单一预警指标种类的强冲击危险判定指标时， 能够预测到 的历史冲击显现事件次数； 0078 Yiz为第i项单一预警指标种类中等冲击危险判定指标的预警准确率； 说明书 4/12 页 9 CN 111476420 A 9 0079 Aiz为采用kqi作为第i项单一预警指。

25、标种类的强冲击危险判定指标时， 能够预测到 的历史冲击显现事件次数； 0080 Yir为第i项单一预警指标种类弱冲击危险判定指标的预警准确率； 0081 Air为采用kqi作为第i项单一预警指标种类的弱冲击危险判定指标时， 能够预测到 的历史冲击显现事件次数； 0082 B为步骤(1)中采掘工作面收集到的动力显现事件数量； 0083 漏报率： 0084 漏报率Liq1-Yiq 0085 漏报率Liz1-Yiz 0086 漏报率Lir1-Yir 0087 Liq为第i项单一预警指标种类强冲击危险判定指标的预警漏报率； 0088 Liz为第i项单一预警指标种类中等冲击危险判定指标的预警漏报率； 0。

26、089 Lir为第i项单一预警指标种类弱冲击危险判定指标的预警漏报率； 0090 误报率： 0091 0092 0093 0094 其中： 0095 Xiq为第i项单一预警指标种类强冲击危险判定指标的预警误报率； 0096 Ciq为采用kqi作为第i项单一预警指标种类的强冲击危险判定指标时， 工作面采掘 期间具备强冲击危险的总次数； 0097 Xiz为第i项单一预警指标种类中等冲击危险判定指标的预警误报率； 0098 Ciz为采用kqi作为第i项单一预警指标种类的中等冲击危险判定指标时， 工作面采 掘期间具备中等冲击危险的总次数； 0099 Xir为第i项单一预警指标种类弱冲击危险判定指标的预。

27、警误报率； 0100 Cir为采用kqi作为第i项单一预警指标种类的弱冲击危险判定指标时， 工作面采掘 期间具备弱冲击危险的总次数。 0101 本发明与现有技术相比， 具备以下几个优点： 0102 1、 通过引入预警指标预警效能评价指标， 包括准确率、 漏报率、 误报率， 可量化评 价单一预警指标每日预报准确性， 有无漏报可能， 有无误报的可能， 使得现场监测预警效果 更加准确。 此外， 通过选择合理的准确率、 漏报率、 误报率， 控制过大的误报率会导致工作面 频繁预警， 可以避免对工作面采掘作业的过度影响。 0103 2、 依据步骤(7)的有效预警指标筛选， 可以筛选合理的预警指标类型， 删。

28、除无效预 警指标类型， 使得现场监测预警工作更加准确。 0104 3、 通过当日动载综合预警值， 可以有效反应出采掘工作面当日动载扰动条件下发 说明书 5/12 页 10 CN 111476420 A 10 生冲击危险的可能性， 对于现场安全管理指导意义明显。 本发明的对动静载叠加型冲击地 压现场监测与防治指导意义较大。 附图说明 0105 图1为本发明方法的流程原理示意图。 0106 图2为每日104J以上事件频次倒序排列及动力显现事件分布图。 0107 图3为每日微震事件总能量倒序排列及动力显现事件分布图。 0108 图4为每刀能量降序倒序排列及动力显现事件分布图。 0109 图5为每日事。

29、件总频次降序倒序排列及动力显现事件分布图。 具体实施方式 0110 下面通过具体实施例并结合附图对本发明作进一步详细说明。 0111 如图1所示， 本发明提供的一种微震监测冲击地压预警指标优化方法， 其具体步骤 如下： 0112 步骤1： 本实施例涉及的实验回采工作面采用KJ551微震监测系统， 于2019年5月投 入使用， 截止2019年12月22日， 系统在实验工作面区域共监测到微震事件23845个。 期间记 录动力显现8次， 明显震感强矿震事件1次。 0113 步骤2： 工作面回采期间借鉴本矿井邻近回采工作面微震监测预警指标， 选择单一 大能量事件、 104J以上事件频次、 每日微震事件。

30、总能量、 每日微震事件总频次、 每刀能量， 此 外实际预警工作中发现能频差异系数预警指标实际预警效果较好， 亦作为一项单一预警指 标类型。 0114 步骤3： 对每项单一预警指标在工作面采掘期间微震监测数据进行倒序排列 0115 对实验工作面回采期间监测到的微震事件按照104J以上事件频次、 每日微震事件 总能量、 每日微震事件总频次、 每刀能量等预警指标类型进行统计计算， 并将计算结果进行 倒序排列， 按横坐标为时间轴绘制倒序排列图表。 0116 步骤4： 在倒序排列图中标识历次动力显现、 矿震事件 0117 将收集到的本工作面历次动力显现、 矿震事件按实际发生时间表示在倒序排列图 表中， 。

31、结果如图2图5所示。 0118 步骤5： 按以下原则确定每项单一预警指标具体预警值 0119 设定强冲击危险预警指标能够预警80冲击事件， 在倒序排列图表中在选定单 一预警指标具体预警值|kqi|； 0120 选定8月18日对应的事件频次作为强冲击危险预警指标， 104J以上事件每日频次 16次， 取|kq1|16。 0121 选定8月01日对应的当日微震事件总能量作为强冲击危险预警指标， 即867955J， 取|kq2|8.67105J。 0122 选定9月03日对应的当日微震事件总能量作为强冲击危险预警指标， 即140836J， 取|kq3|1.40105J。 0123 选定7月29日对应。

32、的每日事件总频次作为强冲击危险预警指标， 即85个， 取|kq4| 85。 说明书 6/12 页 11 CN 111476420 A 11 0124 设定中等冲击危险预警指标能够预警100冲击事件， 在倒序排列图表中在选 定单一预警指标具体预警值|kzi|； 0125 选定9月21日对应的事件频次作为中等冲击危险预警指标， 即104J以上事件每日 频次12次， 取|kz1|12。 0126 选定9月03日对应的事件频次作为中等冲击危险预警指标， 即563345J， 取|kz2| 5.63105J。 0127 选定8月01日对应的事件频次作为中等冲击危险预警指标， 即108494J， 取|kz3。

33、| 1.08105J。 0128 选定8月01日对应的事件频次作为中等冲击危险预警指标， 取|kz4|68。 0129 设定弱等冲击危险预警指标能够预警100强震动事件， 在倒序排列图表中在 选定单一预警指标具体预警值|kri|； 0130 选定弱冲击危险预警指标能够预警100以上强震动事件， 即104J以上事件每日 频次1次， 累计强震动事件56次， 占总数的100， 取|kr1|1 0131 选定弱冲击危险预警指标能够预警100以上强震动事件， 即91341J， 取|kr2| 9.13104J。 0132 选定弱冲击危险预警指标能够预警100以上强震动事件， 即22835.25J， 取|k。

34、r3| 2.28104J。 0133 设定弱冲击危险预警指标能够预警100强震动事件， 取|kr4|29。 0134 设定无冲击危险预警指标为|kwi|kri|； 0135 取|kw1|0； 取|kw2|9.13104； 取|kw3|2.28104； 取|kw4|29； 0136 步骤6： 计算每项单一预警指标预警效能指标 0137 依据步骤6计算结果如下： 0138 表1 104J以上事件频次预警指标预警效能指标 0139 0140 表2每日微震事件总能量预警指标预警效能指标 说明书 7/12 页 12 CN 111476420 A 12 0141 0142 表3每刀能量预警指标预警效能指标。

35、 0143 0144 表4每日事件总频次指标预警效能指标 0145 0146 表5能-频差异系数预警效能指标 说明书 8/12 页 13 CN 111476420 A 13 0147 0148 步骤7： 有效预警指标筛选 0149 经排查步骤6中各项单一预警指标预警效能指标不存在同时具备准确率最小、 漏 报率最大、 误报率最大的条件， 均予以保留。 0150 汇总以上结果， 得到实验工作面微震监测系统预警指标优化结果， 如表6所示。 0151 表6实验工作面微震系统预警指标优化结果 0152 说明书 9/12 页 14 CN 111476420 A 14 0153 0154 步骤8： 计算每项。

36、单一预警指标预警权重 0155 计算每项单一预警指标实际预警率 0156 计算结果如表7所示。 0157 表7每项单一预警指标实际预警率计算结果 0158 0159 选择每项单一预警指标实际预警率最小值 说明书 10/12 页 15 CN 111476420 A 15 0160 EminqMIN(Eiq)(i0,1,2.n)0.197183 0161 EminzMIN(Eiz)(i0,1,2.n)0.211268 0162 EminrMIN(Eir)(i0,1,2.n)0.230047 0163 计算每项单一预警指标预警权重 0164 计算结果如表8所示。 0165 表8每项单一预警指标预警权。

37、重计算结果 0166 0167 步骤9： 计算当日综合动载预警值 0168 计算当日单一预警指标的预警权重 0169 该方法自2019年12月25日试用， 以2019年12月25日实验工作面微震监测数据为例 说明当日综合动载预警值计算过程。 0170 2019年12月25日， 实验工作面各单一预警指标实际监测值如下： 0171 说明书 11/12 页 16 CN 111476420 A 16 0172 0173 计算当日动载综合预警值 0174 0175 步骤10： 每日进行综合预警 0176 故判定当日动载扰动条件下冲击危险程度为中等。 说明书 12/12 页 17 CN 111476420 A 17 图1 图2 说明书附图 1/2 页 18 CN 111476420 A 18 图3 图4 图5 说明书附图 2/2 页 19 CN 111476420 A 19 。