《一种冲击地压的微震多参量预警方法.pdf》由会员分享，可在线阅读，更多相关《一种冲击地压的微震多参量预警方法.pdf（8页珍藏版）》请在专利查询网上搜索。
1、10申请公布号CN104153814A43申请公布日20141119CN104153814A21申请号201410284705022申请日20140623E21F17/1820060171申请人中国矿业大学地址221000江苏省徐州市泉山区大学路1号72发明人陆菜平刘广建张农刘洋左源耀74专利代理机构徐州市淮海专利事务所32205代理人华德明54发明名称一种冲击地压的微震多参量预警方法57摘要本发明公开了一种冲击地压的微震多参量预警方法，包括在矿井中安装微震监测系统；利用矿井微震监测系统，测量微震能量和频次，计算断层总面积AT、缺震B值、ZMAP值并记录该四项微震信号；将所得数据进行分析，绘制。
2、出能量、频次时间序列图；断层总面积演变图；缺震B值演变图以及ZMAP值演变图四项演变图；根据所绘制的演变图，分析是否在同一时间段内出现矿震频次降低、能量增加；断层总面积值出现高值异常；缺震B值低值异常；|Z|2四项指标，若同一时间段内四项指标中出现24项，则判断会出现冲击地压的情况。有益效果是提高预警冲击地压的准确性。51INTCL权利要求书1页说明书4页附图2页19中华人民共和国国家知识产权局12发明专利申请权利要求书1页说明书4页附图2页10申请公布号CN104153814ACN104153814A1/1页21一种冲击地压的微震多参量预警方法，其特征在于，具体步骤是A在矿井中安装微震监测系。
3、统；B利用矿井微震监测系统，测量微震能量和频次，计算断层总面积AT、缺震B值、ZMAP值并记录该四项微震信号；其中，断层总面积AT介于频度和能量之间，根据求得AT值；其中，NK为时间T至TT间隔内能级为K的微震数目，即微震弹性波能量在10K10K1之间的微震数目,K0为所统计微震的下限，K为每个微震的能级；ZMAP值通过分析平均震级样本的变化即可了解研究区域不同时段平均的变化情况，根据求得Z值；其中，为整个时间区间上对所有平均震级样本的算术平均值，为要考察的时间区段内样本平均震级样本的算术平均值；M和M分别是两样本的标准差；缺震B值根据宇津德治公式计算古登堡公式求得B值，B值的变化依赖于平均能。
4、级其中为统计时段内微震的平均能级，M0为起始能级，N为微震总数，为短期平均能级；C将所得数据进行分析，绘制出能量、频次时间序列图；断层总面积演变图；缺震B值演变图以及ZMAP值演变图四项演变图；D根据所绘制的演变图，分析是否在同一时间段内出现矿震频次降低、能量增加；断层总面积值出现高值异常；缺震B值低值异常；|Z|2四项指标，若同一时间段内四项指标中出现24项，则判断会出现冲击地压的情况。2根据权利要求1所述的一种冲击地压的微震多参量预警方法，其特征在于，D步骤中，同时出现出现矿震频次降低、能量增加；断层总面积值出现高值异常；缺震B值低值异常；|Z|2四种情况，则判断会出现冲击地压的情况。权利。
5、要求书CN104153814A1/4页3一种冲击地压的微震多参量预警方法技术领域0001本发明涉及一种冲击地压的微震多参量预警方法，属于微震监测预警技术领域。背景技术0002近年来，随着我国煤矿开采深度的逐年增加，很多矿区冲击地压问题越来越严重。煤矿冲击地压煤岩体中的压力超过强度极限，聚积在巷道或采场周围煤岩体中的能量突然释放，在井巷发生爆炸性事故，动力将煤岩抛向巷道，同时发出剧烈声响，是一种开采诱发的矿山地震，更是一种灾害，不仅造成井巷破坏、人员伤亡、地面建筑物破坏，还会引发瓦斯、煤尘爆炸以及水灾，干扰通风系统，严重威胁着煤矿的生产安全。由于这种灾害发生的时间与地域方面具有复杂性和突发性，因。
6、此，对其进行预测预警也成为世界性难题。传统的利用矿压预测冲击地压的方法局限性较大，利用微震监测的方法大多指标单一，预警准确率低；费时费力，人员需要现场记录，安全性低，经济上劣势突出。发明内容0003针对上述现有技术存在的问题，本发明提供一种冲击地压的微震多参量预警方法，用能量与频次、断层总面积、缺震、ZMAP值四个指标来作为预测冲击地压的前兆性指标。0004为了实现上述目的，本发明采用的技术方案是一种冲击地压的微震多参量预警方法，具体步骤是0005A在矿井中安装微震监测系统；0006B利用矿井微震监测系统，测量微震能量和频次，计算断层总面积AT、缺震B值、ZMAP值并记录该四项微震信号；其中，。
7、0007断层总面积AT介于频度和能量之间，根据求得AT值；其中，NK为时间T至TT间隔内能级为K的微震数目，即微震弹性波能量在10K10K1之间的微震数目,K0为所统计微震的下限，K为每个微震的能级；0008ZMAP值,简写为Z值通过分析平均震级样本的变化即可了解研究区域不同时段平均的变化情况，根据求得Z值；其中，为整个时间区间上对所有平均震级样本的算术平均值，为要考察的时间区段内样本平均震级样本的算术平均值；M和M分别是两样本的标准差；0009缺震B值根据宇津德治公式计算古登堡公式求得B说明书CN104153814A2/4页4值，B值的变化依赖于平均能级其中为统计时段内微震的平均能级，M0为。
8、起始能级，N为微震总数，为短期平均能级；0010C将所得数据进行分析，绘制出能量、频次时间序列图；断层总面积演变图；缺震B值演变图以及ZMAP值演变图四项演变图；0011D根据所绘制的演变图，分析是否在同一时间段内出现矿震频次降低、能量增加；断层总面积值出现高值异常；缺震B值低值异常；|Z|2四项指标，若同一时间段内四项指标中出现24项，则判断会出现冲击地压的情况。0012优选地，D步骤中，同时出现出现矿震频次降低、能量增加；断层总面积值出现高值异常；缺震B值低值异常；|Z|2四种情况，则判断会出现冲击地压的情况。0013本发明中，所述四项指标的基本原理分别是0014一、能量、频次0015微震。
9、监测技术当中，使用比较多的两个参数是能量和频次。微震监测系统的主要功能是对全矿范围内的矿震进行监测，自动记录矿震活动，通过分析和计算矿震能量和频次的数值。研究结果表明，频次主要反映煤岩体变形及微破裂的频次，能量主要反映煤岩体的受载程度及变形破裂强度。能量与频次的结合分析能较好的反映出煤岩体的破裂变化。0016二、断层总面积AT0017用频度和能量这两个参量评价微震活动性的高低，由于低强度事件的数目通常远多于高强度事件，所以频度实际上是由低强度事件所决定的，但总能量却是由为数不多的几个高强度事件所决定的，而低强度事件被忽略了。为了克服这一矛盾，定义一个介于频度和能量之间的AT，它既考虑到微震频度。
10、，同时也包含了微震能量。00180019式中NK为时间T至TT间隔内能级为K的微震数目，即微震弹性波能量在10K10K1之间的微震数目。显然，当L1时，AT仅表示微震数目；L10时，AT表示微震总能量。若取L45，则LGL2/3，此时上式写为00200021设微震弹性波能量E正比于微震能量E0，即EE0，因能量比例于震源体积即震源体大小的三次方，即正比于震源断层面积S的3/2次方，E0S3/2，因此，SE2/3102K/3；即L45时，AT相当于断层总面积。这就是AT的物理意义。为方便起见，将上式写为00220023式中K0为所统计微震的下限，K为每个微震的能级。0024较强微震往往发生在活动。
11、断裂带上，微震在孕育过程中积聚大量的能量，这些能量要通过裂隙发育的地方来释放，因此断层裂隙多的地方发生强矿震的可能性较大，理论上在强能量释放前，AT出现高值异常，表明在强能量释放前，微震活动性增强。0025三、ZMAP值，简写为Z值说明书CN104153814A3/4页50026通过分析平均震级样本的变化即可了解研究区域不同时段平均的变化情况002700280029式中，为整个时间区间上对所有平均震级样本的算术平均值，是一个较为稳定的量，表征研究区域的背景特征；为要考察的时间区段内样本平均震级样本的算术平均值；M和M分别是两样本的标准差。0030因为Z是从相同母体中抽取出的平均值之差，近似服从。
12、于标准正态分布，因此它具有正态分布的显著特征，即Z164，具有90的显著水平；Z196，具有95的显著水平；当Z257时，具有99的显著水平。据此，可分析平均震级变化所反映的平静特征显著性。0031当Z0时，微震发生率与背景微震发生率相同；当Z0时，该时段微震发生率下降，以此来判断强能量释放前的微震活动平静。0032无论是当Z2，还是Z2为异常临界值。0033四、缺震B值0034根据宇津德治公式，古登堡公式中的B值可为003500360037这里为统计时段内微震的平均能级，M0为起始能级，N为微震总数。由宇津公式可知，B值的变化依赖于平均能级在正常微震活动情况下，一个区域的B值基本可看作稳定的。
13、，则平均能级也应该是比较稳定的,即平均能级代表一个区域的平均矿震水平。如果一个区域在一段时间内的短期平均能级比长期平均能级偏小的话，那么这个区域就可能发生一些较大微震来补足这个长期平均能级的缺额。这就是缺震的基本含义，“缺震”意味着将要发生缺失能级的趋势。如果一个区域在一段时间内的短期平均能级比长期平均能级偏大的话，那么这个区域将不会发生一些较大微震，这就是过震的基本含义。0038本发明的有益效果是通过对能量与频次、断层总面积AT、缺震B值、ZMAP值四个指标来作为预测冲击地压的前兆性指标，从而提高预警冲击地压的准确性。附图说明0039图1为实施例中能量、频次时间序列图；说明书CN104153。
14、814A4/4页60040图2为实施例中断层总面积演变图；0041图3为实施例中ZMAP值演变图；0042图4为实施例中缺震B值演变图。具体实施方式0043下面结合实施例对本发明作进一步说明。0044本发明的具体步骤是0045A在矿井中安装微震监测系统；0046B利用矿井微震监测系统，测量微震能量和频次，计算断层总面积AT、缺震B值、ZMAP值并记录该四项微震信号；其中，0047断层总面积AT介于频度和能量之间，根据求得AT值；其中，NK为时间T至TT间隔内能级为K的微震数目，即微震弹性波能量在10K10K1之间的微震数目,K0为所统计微震的下限，K为每个微震的能级；0048ZMAP值通过分析。
15、平均震级样本的变化即可了解研究区域不同时段平均的变化情况，根据求得Z值；其中，为整个时间区间上对所有平均震级样本的算术平均值，为要考察的时间区段内样本平均震级样本的算术平均值；M和M分别是两样本的标准差，；0049缺震B值根据宇津德治公式计算古登堡公式求得B值，B值的变化依赖于平均能级其中为统计时段内微震的平均能级，M0为起始能级，N为微震总数，为短期平均能级；0050C将所得数据进行分析，绘制出能量、频次时间序列图；断层总面积演变图；缺震B值演变图以及ZMAP值演变图四项演变图；0051D根据所绘制的演变图，分析是否在同一时间段内出现矿震频次降低、能量增加；断层总面积值出现高值异常；缺震B值低值异常；|Z|2四项指标，若同一时间段内四项指标中出现24项，则判断会出现冲击地压的情况。0052D步骤中，若同时出现矿震频次降低、能量增加；断层总面积值出现高值异常；缺震B值低值异常；|Z|2四种情况，则出现冲击地压的情况会大大增加，通过这种方法提高预警冲击地压的准确性。0053实施例为某矿一起冲击地压事故，根据本发明的预警方法，从图中可以分析出，该起事故于2013年3月15日发生。具体四项分析的时间演变图见附图。说明书CN104153814A1/2页7图1图2说明书附图CN104153814A2/2页8图3图4说明书附图CN104153814A。