一种基于冲击风险熵的冲击地压预警方法技术

本发明专利技术提出了一种基于冲击风险熵的冲击地压预警方法，通过微震监测手段获取微震频度和微震总能量，获取通过地音监测手段获取的班频次变化率和班能量变化率，将所述微震频度、微震总能量、班频次变化率和班能量变化率作为冲击风险熵指标；以专家打分表为基础，利用层次分析法确定各冲击风险熵指标的权重值；根据各冲击风险熵指标构建熵增熵减估算模型，再结合确定的冲击风险熵指标的权重值计算出熵增熵减估算模型的熵增值和熵减值；将计算出的熵增值和熵减值带入建立的突变模型中，判断冲击地压的风险，进而预警。本发明专利技术可结合突变理论模型对冲击地压中的冲击风险熵指标进行分析，以达到预警冲击地压危险的效果，从而预防此类事故的发生。事故的发生。事故的发生。

【技术实现步骤摘要】
一种基于冲击风险熵的冲击地压预警方法


[0001]本专利技术属于煤矿工程事故应急预警
，具体涉及一种基于冲击风险熵的冲击地压预警方法。

技术介绍

[0002]随着我国煤矿开采强度和深度的不断增加，矿山动力灾害问题越来越严重，尤其是其中的冲击地压问题更是严重的影响了安全和经济问题。近年来，因为冲击地压遇难的人数逐渐剧增，因此国家越来越关注冲击地压的预警问题。很明显，冲击地压危险已成为我国威胁煤矿安全生产的重大灾害之一。因此，准确监测矿山并预警冲击地压危险是必要且必须重视的问题。
[0003]微震与地音监测技术是利用岩体受力变形和破坏后本身发射出的微震和声波来进行监测工程岩体稳定性的技术方法。因此，微震法结合地音法可以使监测预警更准确。另外，冲击地压又具有震动性、瞬时性、复杂性以及突发性，但目前针对解决冲击地压突发性的有效方法比较有限，所以通过一种冲击风险熵对冲击地压突发性进行预警的问题亟待解决。
[0004]因此，现在针对煤矿矿山的突发冲击地压问题亟需一种有效精准的预警方法来进行监测预警，从而采用应急措施及卸压防治方法进行应对处理。

技术实现思路

[0005]本专利技术所要解决的技术问题在于针对上述现有技术的不足，提供了一种基于冲击风险熵的冲击地压预警方法，该方法选择微震监测手段结合地音监测手段选定冲击风险熵指标，可结合突变理论模型对冲击地压中的冲击风险熵指标进行分析，以达到预警冲击地压危险的效果，从而预防此类事故的发生。
[0006]为解决上述技术问题，本专利技术采用的技术方案是：一种基于冲击风险熵的冲击地压预警方法，其特征在于，该方法包括：
[0007]S1、通过微震监测手段获取微震频度和微震总能量，获取通过地音监测手段获取的班频次变化率和班能量变化率，将所述微震频度、微震总能量、班频次变化率和班能量变化率作为冲击风险熵指标；
[0008]S2、以专家打分表为基础，利用层次分析法确定各冲击风险熵指标的权重值；
[0009]S3、根据各冲击风险熵指标构建熵增熵减估算模型，再结合确定的冲击风险熵指标的权重值计算出熵增熵减估算模型的熵增值和熵减值；
[0010]S4、将计算出的熵增值和熵减值带入建立的突变模型中，判断冲击地压的风险，进而预警；
[0011]所述突变模型的判别公式为：
[0012]Δ＝8u3+27v2，其中u表示熵增值，v表示熵减值。
[0013]优选地，S2中以专家打分表为基础，利用层次分析法确定各冲击风险熵指标的权
重值，具体包括：
[0014]S201、构建冲击地压预警层次分析模型，确定目标层、中间层和最底层；
[0015]S202、以层次分析模型为基础设计专家打分表，收集专家打分表，提取其相关的数据，建立微震频度、微震总能量、班频次变化率和班能量变化率的判断矩阵；
[0016]S203、求解微震频度、微震总能量、班频次变化率和班能量变化率两两判断矩阵的特征根，得出权重向量；
[0017]S204、对判断矩阵进行一致性检验，采用三个指标进行检验：一致性指标CI、一致性平均随机指标RI、相对一致性指标CR，若符合一致性检验条件，则接受该权重向量；若不符合，则调整标度值重新计算；
[0018]S205、计算微震频度、微震总能量、班频次变化率和班能量变化率的最终权重值。
[0019]优选地，S3中根据所述冲击风险熵指标构建熵增熵减估算模型，再结合确定的冲击风险熵指标的权重值计算出熵增熵减估算模型的熵增值和熵减值，具体包括：
[0020]S301、建立熵增熵减矩阵：
[0021]分别构建微震频度、微震总能量、班频次变化率和班能量变化率这四个冲击风险熵指标n个时刻的熵增矩阵B
z
和熵减矩阵B
j
；
[0022][0023]所述熵增矩阵B
z
中b
z,ki
表示第i个风险熵指标第k时刻的风险熵增值，i＝1，2，3，4；
[0024][0025]所述熵减矩阵B
j
中b
j,ki
表示第i个风险熵指标中第k时刻的风险熵减值，j＝1，2，3，4；
[0026]S302、计算熵增熵减值；
[0027]u＝G
u
×
W
T
＝E
×
B
z
×
W
T
；
[0028]v＝T
v
×
W
T
＝E
×
B
j
×
W
T
；
[0029]其中，u表示熵增值，v表示熵减值，W表示冲击风险熵指标权重矩阵，W
T
为W的转置矩阵；E表示1
×
n的矩阵，G
u
表示冲击风险熵指标的熵增值之和矩阵，T
v
表示冲击风险熵指标
的熵减值之和矩阵；
[0030]G
u
＝E
×
B
z
＝[z1,z2,z3,z4]；
[0031]T
v
＝E
×
B
j
＝[j1,j2,j3,j4]；
[0032]其中，z1～z4表示四个冲击风险熵指标所有时刻的熵增值之和，j1～j4表示四个冲击风险熵指标所有时刻的熵减值之和；
[0033]W＝[ω1,ω2,ω3,ω4][0034]其中，ω1～ω4表示微震频度、微震总能量、班频次变化率和班能量变化率四个冲击风险熵指标的依次对应的权重值。
[0035]优选地，将计算出熵增熵减估算模型的熵增值和熵减值带入建立的突变模型中，判断冲击地压的风险，进而预警，具体包括：
[0036]以判别公式Δ作为突变模型，将计算出熵增熵减估算模型的熵增值和熵减值带入判别公式Δ，若Δ>0，说明冲击地压的检测处于稳定状态，暂时不会突发危险；若Δ<0，说明受载煤岩系统风险发生突变，有严重的冲击风险会发生，应实施有效的防治方法，做好提前的应急工作；若Δ＝0，说明风险达到临界状态，随时会因为一些外在不利因素的加入导致风险发生，应保持持续监测，并适当的采用防治方法。
[0037]本专利技术与现有技术相比具有以下优点：
[0038]1、本专利技术微震法结合地音法可以使监测预警更准确，能够解决冲击地压突发性，可结合突变理论模型对冲击地压中的冲击风险熵指标进行分析，以达到预警冲击地压危险的效果，从而预防此类事故的发生
[0039]下面通过附图和实施例对本专利技术的技术方案作进一步的详细说明。
附图说明
[0040]图1为本专利技术实施例1公开的一种基于冲击风险熵的冲击地压预警方法的流程示意图。
具体实施方式
[0041]实施例1
[0042]如图1所示，本专利技术实施例的一种基于冲击风险熵的冲击本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于冲击风险熵的冲击地压预警方法，其特征在于，该方法包括：S1、通过微震监测手段获取微震频度和微震总能量，获取通过地音监测手段获取的班频次变化率和班能量变化率，将所述微震频度、微震总能量、班频次变化率和班能量变化率作为冲击风险熵指标；S2、以专家打分表为基础，利用层次分析法确定各冲击风险熵指标的权重值；S3、根据各冲击风险熵指标构建熵增熵减估算模型，再结合确定的冲击风险熵指标的权重值计算出熵增熵减估算模型的熵增值和熵减值；S4、将计算出的熵增值和熵减值带入建立的突变模型中，判断冲击地压的风险，进而预警；所述突变模型的判别公式为：Δ＝8u3+27v2，其中u表示熵增值，v表示熵减值。2.根据权利要求1所述的一种基于冲击风险熵的冲击地压预警方法，其特征在于，S2中以专家打分表为基础，利用层次分析法确定各冲击风险熵指标的权重值，具体包括：S201、构建冲击地压预警层次分析模型，确定目标层、中间层和最底层；S202、以层次分析模型为基础设计专家打分表，收集专家打分表，提取其相关的数据，建立微震频度、微震总能量、班频次变化率和班能量变化率的判断矩阵；S203、求解微震频度、微震总能量、班频次变化率和班能量变化率两两判断矩阵的特征根，得出权重向量；S204、对判断矩阵进行一致性检验，采用三个指标进行检验：一致性指标CI、一致性平均随机指标RI、相对一致性指标CR，若符合一致性检验条件，则接受该权重向量；若不符合，则调整标度值重新计算；S205、计算微震频度、微震总能量、班频次变化率和班能量变化率的最终权重值。3.根据权利要求1所述的一种基于冲击风险熵的冲击地压预警方法，其特征在于，S3中根据所述冲击风险熵指标构建熵增熵减估算模型，再结合确定的冲击风险熵指标的权重值计算出熵增熵减估算模型的熵增值和熵减值，具体包括：S301、建立熵增熵减矩阵：分别构建微震频度、微震总能量、班频次变化率和班能量变化率这四个冲击风险熵指标n个时刻的熵增矩阵B
z
和熵减矩阵B
j
；所述熵增矩阵B
z
中b
z,ki
表示第i个风险熵指标第k时刻的风险熵增值，i＝1，2，3，4；
所述熵减矩阵B
j<...

【专利技术属性】
技术研发人员：金佩剑，鲁鑫，李杭，许金融，葛素伶，娄胜阳，范鹏震，崔永祥，
申请(专利权)人：吉林建筑大学，
类型：发明
国别省市：