一种煤岩动力灾害危险的电磁辐射和地音监测预警方法与流程

本发明涉及一种煤岩动力灾害危险的监测预警指标和预警方法，具体是一种基于电磁辐射和地音监测技术的煤岩动力灾害危险性预警指标和预警方法，主要用于煤矿冲击地压和煤与瓦斯突出等动力灾害的实时监测和定量化预警。

背景技术：

煤岩动力灾害主要包含冲击地压、煤与瓦斯突出等，受采动作用影响常发生在高应力集中区，其破坏性大，常造成井巷破坏、人员伤亡和设备损坏。由于煤炭资源长期开采，浅部资源逐渐枯竭，我国煤矿已逐步进入深部开采，动力灾害的频次、强度及其破坏性均呈上升趋势。准确高效的监测预警是煤矿提前采取解危措施、确保安全生产的前提和保障。

目前，煤岩动力灾害的监测预警方法主要包括钻屑量监测、电磁辐射监测、地音监测、煤体应力监测、微震监测等。钻屑量监测工程量大、成本高、监测范围小、获得的信息量有限，现阶段已逐渐成为一种辅助的监测手段；电磁辐射监测可对煤岩应力集中以及微破裂、大破裂进行监测，监测预警指标主要为能量和脉冲，预警效果有限；地音监测主要用于对动力灾害孕育发展初期的高频低能微破裂进行监测，监测预警指标主要为能量和频次，预警效果有待提高；煤体应力监测属于点监测，监测范围较小，钻孔应力计与煤体钻孔耦合不佳时起不到很好的预警作用；微震监测可实现对煤岩破裂的时序监测和对破裂源的空间定位，监测范围广，但微震主要用于对低频高能的大破裂进行监测，难以捕捉到动力灾害孕育初期的微破裂事件。

电磁辐射和地音监测数据量庞大，对微破裂事件十分敏感，能反映出众多细节信息，是较为理想的动力灾害监测预警手段。然而，由于缺乏对电磁辐射和地音前兆参量的深入细致研究，仅简单使用能量、脉冲数、频次等原始统计参数进行监测预警，导致误报、漏报出现的次数较多，已影响到矿井安全生产。可见，有必要通过对电磁辐射和地音前兆参量的深入分析，提出一种新的基于电磁辐射和地音监测的煤岩动力灾害危险性预警指标和预警方法，以达到提高煤岩动力灾害监测预警可靠性的目的。

本发明针对以上问题，提供一种煤岩动力灾害危险的电磁辐射和地音监测预警指标和预警方法，以便实现对煤岩动力灾害孕育发展过程的实时监测和动力灾害的可靠预警。

技术实现要素：

为实现上述目的，本发明提供如下技术方案：一种煤岩动力灾害危险的电磁辐射和地音监测预警方法，其特征在于：该方法是在煤体和岩体上分别布置电磁辐射和地音传感器，实时采集电磁辐射的能量和脉冲数以及地音信号的能量和频次作为原始数据，其具体包括如下步骤：

步骤(1)根据电磁辐射和地音监测系统实时采集到的监测数据计算研究参量在一段时间内的加权平均值p(t)；

步骤(2)计算研究参量加权平均值p(t)的偏差值d(t)；

步骤(3)计算工作面正常开采期间研究参量偏差值绝对值|d(t)|的平均值并将作为偏差阈值；

步骤(4)计算研究参量偏差值d(t)一天内大于偏差阈值的次数ds；

步骤(5)对ds进行归一化处理得到监测预警指标ε；

步骤(6)根据监测预警指标和预警方法对工作面的动力灾害危险状态实时预警，确定动力灾害的危险等级。

进一步，作为优选，在所述步骤(1)中，所述的研究参量为电磁辐射能量、电磁辐射脉冲数、地音能量或地音频次中的一种或者几种，研究参量加权平均值p(t)为一段时间内研究参量的累积和除以这段时间的时间窗长度。

进一步，作为优选，所述一段时间取10分钟。

进一步，作为优选，在所述步骤(2)中所述的研究参量偏差值d(t)由下式计算：其中其中tt表示时间间隔，这一时间间隔与t时刻相关，是t时刻之前一段时间到t时刻之间的时间间隔；n为tt时间间隔内的研究参量加权平均值p(t)的个数。

进一步，作为优选，t时刻之前一段时间取为24小时。

进一步，作为优选，在所述步骤(3)中，所述的正常回采期间是指工作面不受断层、褶曲等地质构造影响且回采期间未出现顶板来压、动力显现等异常情况的一段时间，这段时间可以是一个月或者数月。

进一步，作为优选，在所述步骤(5)中所述的监测预警指标ε的取值由下式计算：其中ds-max为一段时间内ds的最大值，ds-min为一段时间内ds的最小值。

进一步，作为优选，该一段时间取当前时刻至前一次动力显现时刻之间的时间段。

进一步，作为优选，在所述步骤(6)中，所述的动力灾害的危险等级分为无危险、弱危险、中等危险和强危险四级，危险等级的确定方法和依据为：ε＜0.5时为无危险；0.5≤ε＜0.65时为弱危险；0.65≤ε＜0.8时为中等危险；0.8≤ε时为强危险。

与现有技术相比，本发明的有益效果是：

本发明提供的一种煤岩动力灾害危险的电磁辐射和地音监测预警指标和预警方法，其通过对电磁辐射和地音监测设备监测到的原始数据进行深入分析，求出特定时间窗内研究参量加权平均值的偏差值以及偏差阈值，并得到每日研究参量偏差值出现大于偏差阈值的次数ds，并将ds作归一化处理，得到新的参量ε作为监测预警指标，通过ε的实时值来实时确定监测范围的煤岩动力灾害危险状态。该预警指标和预警方法物理意义明确、易于编程，可实现对动力灾害危险的实时和定量化预警，为矿井安全生产起到积极作用。

附图说明

图1是本发明的实施流程图；

图2是实施例的由地音能量计算的一段时间的加权平均值pde(t)；

图3是实施例的由地音能量计算的偏差值dde(t)；

图4是实施例的由地音能量计算的的偏差值在一天内大于偏差阈值的次数dsde；

图5是实施例的由地音能量计算的归一化监测预警指标ε。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

请参阅图1-5，本发明提供一种技术方案：一种煤岩动力灾害危险的电磁辐射和地音监测预警方法，在煤体和岩体上分别布置电磁辐射和地音传感器，实时采集电磁辐射的能量和脉冲次数以及地音信号的能量和频次作为原始数据。其特征在于：包括如下步骤：

步骤(1)根据电磁辐射和地音监测系统实时采集到的监测数据计算研究参量在一段时间内的加权平均值p(t)，所述的研究参量可以为电磁辐射能量、电磁辐射脉冲数、地音能量、地音频次，研究参量加权平均值p(t)为一段时间(可以取10分钟)内研究参量的累积和除以这段时间的时间窗长度。

步骤(2)计算研究参量加权平均值p(t)的偏差值d(t)，所述的偏差值d(t)由下式计算：其中其中tt表示时间间隔，这一时间间隔与t时刻相关，是t时刻之前一段时间(可以取24小时)到t时刻之间的时间间隔，n为tt时间间隔内的研究参量加权平均值p(t)的个数。

步骤(3)计算工作面正常开采期间研究参量偏差值绝对值|d(t)|的平均值并将作为偏差阈值，所述的正常回采期间是指工作面不受断层、褶曲等地质构造影响且回采期间未出现顶板来压、动力显现等异常情况的一段时间，这段时间可以是一个月或者数月。

步骤(4)计算研究参量偏差值d(t)一天内大于偏差阈值的次数ds。

步骤(5)对ds进行归一化处理得到监测预警指标ε，所述的监测预警指标ε的取值由下式计算：其中ds-max为一段时间(一般取当前时刻至前一次动力显现时刻之间的时间段)内ds的最大值，ds-min为一段时间(一般取当前时刻至前一次动力显现时刻之间的时间段)内ds的最小值。

步骤(6)根据监测预警指标和预警方法对工作面的动力灾害危险状态实时预警，确定动力灾害的危险等级，所述的动力灾害的危险等级分为无危险、弱危险、中等危险和强危险四级，危险等级的确定方法和依据为：ε＜0.5时为无危险；0.5≤ε＜0.65时为弱危险；0.65≤ε＜0.8时为中等危险；0.8≤ε时为强危险。

下面结合附图对本发明的一个实施例作进一步的描述：

本实施例以某冲击地压危险工作面的地音监测系统监测到的2016年8月8日至2017年1月25日的地音能量作为原始数据，利用本发明的方法对原始数据处理分析，得到本发明的监测预警指标，并利用预警方法判定该工作面的冲击危险状态，来说明本发明的实施。

选取地音能量作为研究参量，计算每10分钟的地音能量加权平均值pde(t)，计算结果如图2所示。

计算地音能量加权平均值pde(t)的偏差值dde(t)，dde(t)的取值由下式计算：其中tt表示时间间隔，是t时刻之前24小时到t时刻之间的时间间隔，n为tt时间间隔内的地音能量加权平均值的个数，根据图2的结果计算dde(t)，得到计算结果如图3所示。

选取工作面无异常现象的开采期间2016年8月8日至2016年9月8日期间的偏差值的绝对值|dde(t)|来计算其平均值结果为：作为偏差阈值。

计算一天内偏差值dde(t)大于偏差阈值的次数dsde，计算结果如图4所示。

选取图4中2016年11月24日(冲击显现)之后的时间窗来对dsde进行归一化处理，运用公式来计算εde并作为监测预警指标，εde的计算结果如图5所示。

根据监测预警指标εde和冲击危险判定依据来确定冲击危险等级，危险等级的确定方法和依据为：εde＜0.5时为无危险；0.5≤εde＜0.65时为弱危险；0.65≤εde＜0.8时为中等危险；εde≥0.8时为强危险，2016年11月24日至2017年1月25日的工作面冲击危险等级如图5所示。

尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由所附权利要求及其等同物限定。