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

1.本发明属于煤矿工程事故应急预警技术领域，具体涉及一种基于冲击风险熵的冲击地压预警方法。


背景技术：

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


技术实现要素：

5.本发明所要解决的技术问题在于针对上述现有技术的不足，提供了一种基于冲击风险熵的冲击地压预警方法，该方法选择微震监测手段结合地音监测手段选定冲击风险熵指标，可结合突变理论模型对冲击地压中的冲击风险熵指标进行分析，以达到预警冲击地压危险的效果，从而预防此类事故的发生。
6.为解决上述技术问题，本发明采用的技术方案是：一种基于冲击风险熵的冲击地压预警方法，其特征在于，该方法包括：
7.s1、通过微震监测手段获取微震频度和微震总能量，获取通过地音监测手段获取的班频次变化率和班能量变化率，将所述微震频度、微震总能量、班频次变化率和班能量变化率作为冲击风险熵指标；
8.s2、以专家打分表为基础，利用层次分析法确定各冲击风险熵指标的权重值；
9.s3、根据各冲击风险熵指标构建熵增熵减估算模型，再结合确定的冲击风险熵指标的权重值计算出熵增熵减估算模型的熵增值和熵减值；
10.s4、将计算出的熵增值和熵减值带入建立的突变模型中，判断冲击地压的风险，进而预警；
11.所述突变模型的判别公式为：
12.δ＝8u3+27v2，其中u表示熵增值，v表示熵减值。
13.优选地，s2中以专家打分表为基础，利用层次分析法确定各冲击风险熵指标的权
重值，具体包括：
14.s201、构建冲击地压预警层次分析模型，确定目标层、中间层和最底层；
15.s202、以层次分析模型为基础设计专家打分表，收集专家打分表，提取其相关的数据，建立微震频度、微震总能量、班频次变化率和班能量变化率的判断矩阵；
16.s203、求解微震频度、微震总能量、班频次变化率和班能量变化率两两判断矩阵的特征根，得出权重向量；
17.s204、对判断矩阵进行一致性检验，采用三个指标进行检验：一致性指标ci、一致性平均随机指标ri、相对一致性指标cr，若符合一致性检验条件，则接受该权重向量；若不符合，则调整标度值重新计算；
18.s205、计算微震频度、微震总能量、班频次变化率和班能量变化率的最终权重值。
19.优选地，s3中根据所述冲击风险熵指标构建熵增熵减估算模型，再结合确定的冲击风险熵指标的权重值计算出熵增熵减估算模型的熵增值和熵减值，具体包括：
20.s301、建立熵增熵减矩阵：
21.分别构建微震频度、微震总能量、班频次变化率和班能量变化率这四个冲击风险熵指标n个时刻的熵增矩阵bz和熵减矩阵bj；
[0022][0023]
所述熵增矩阵bz中b
z,ki
表示第i个风险熵指标第k时刻的风险熵增值，i＝1，2，3，4；
[0024][0025]
所述熵减矩阵bj中b
j,ki
表示第i个风险熵指标中第k时刻的风险熵减值，j＝1，2，3，4；
[0026]
s302、计算熵增熵减值；
[0027]
u＝gu×wt
＝e
×bz
×wt
；
[0028]
v＝tv×wt
＝e
×bj
×wt
；
[0029]
其中，u表示熵增值，v表示熵减值，w表示冲击风险熵指标权重矩阵，w
t
为w的转置矩阵；e表示1
×
n的矩阵，gu表示冲击风险熵指标的熵增值之和矩阵，tv表示冲击风险熵指标
的熵减值之和矩阵；
[0030]gu
＝e
×bz
＝[z1,z2,z3,z4]；
[0031]
tv＝e
×bj
＝[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所示，本发明实施例的一种基于冲击风险熵的冲击地压预警方法，该方法包括：
[0043]
s1、通过微震监测手段获取微震频度和微震总能量，获取通过地音监测手段获取的班频次变化率和班能量变化率，将所述微震频度、微震总能量、班频次变化率和班能量变化率作为冲击风险熵指标；
[0044]
s2、以专家打分表为基础，利用层次分析法确定各冲击风险熵指标的权重值；
[0045]
s3、根据各冲击风险熵指标构建熵增熵减估算模型，再结合确定的冲击风险熵指标的权重值计算出熵增熵减估算模型的熵增值和熵减值；
[0046]
s4、将计算出的熵增值和熵减值带入建立的突变模型中，判断冲击地压的风险，进而预警。
[0047]
本实施例中采用微震及地音(声发射)监测预警是为了能够保证监测预警的可靠性和准确性。
[0048]
本实施例中，采用微震及地音(声发射)监测预警是为了能够保证监测预警的可靠性和准确性，所述微震频度、微震总能量、班频次变化率和班能量变化率通过以下方法得到：
[0049]
①
微震频度za；
[0050]
经仪器测量得出数据。
[0051]
②
微震总能量ze。
[0052]
经仪器测量得出数据。
[0053]
s202、对于地音(声发射)监测手段，有如下冲击风险熵指标：
[0054]
①
班频次变化率aa；
[0055]
按下式计算得出：
[0056][0057][0058]
式中:
[0059]
ni———当天第i班的地音频次,“三八制”时,i≤3,“四六制”时,i≤4；
[0060]
———前10天第i班地音频次的平均值。
[0061]
②
班能量变化率ae。
[0062]
按下式计算得出：
[0063][0064][0065]
式中:
[0066]ei
———当天第i班的地音能量,单位为焦耳(j),“三八制”时,i≤3,“四六制”时,i≤4；
[0067]
———前10天第i班地音能量的平均值,单位为焦耳(j)。
[0068]
本实施例中，s2中以专家打分表为基础，利用层次分析法确定各冲击风险熵指标的权重值，具体包括：
[0069]
s201、构建冲击地压预警层次分析模型，确定目标层、中间层和最底层；
[0070]
目标层：判断是否有冲击地压风险
[0071]
中间层；所有冲击风险熵指标；所有可以造成冲击地压发生的因素，且可以作为冲击地压预警指标，即为冲击风险熵指标。
[0072]
最底层：备选的冲击风险熵指标，备选的冲击风险熵指标指从所有的冲击风险熵中选中的上述四个(微震频度、微震总能量、班频次变化率和班能量变化率)做实验要用到的冲击风险熵指标。
[0073]
s202、以层次分析模型为基础设计专家打分表，收集专家打分表，提取其相关的数据，建立微震频度、微震总能量、班频次变化率和班能量变化率的判断矩阵；
[0074]
表1专家打分表
[0075][0076]
s203、求解由微震频度、微震总能量、班频次变化率和班能量变化率构造的四阶矩阵的特征根，得出权重向量；
[0077]
s204、对判断矩阵进行一致性检验，采用三个指标进行检验：一致性指标ci、一致性平均随机指标ri、相对一致性指标cr，若符合一致性检验条件，则接受该权重向量；若不符合，则调整标度值重新计算；
[0078]
所述一致性指标ci通过下式得出，
[0079][0080]
式中:λmax为微震频度、微震总能量、班频次变化率和班能量变化率的判断矩阵最大特征值；n为判断矩阵阶数。
[0081]
表2一致性平均随机指标ri
[0082][0083]
相对一致性指标cr通过下式(7)得到，
[0084]
cr＝ci/ri
[0085]
当cr＜0.1时，认为矩阵的一致性符合要求，可以接受所计算权重向量。若cr＞0.1时，认为矩阵一致性不符合要求，需要对矩阵中的值进行分析后重新给出各因素标度值，再进行验算，直至符合一致性要求；
[0086]
s205、计算微震频度、微震总能量、班频次变化率和班能量变化率的最终权重值。
[0087]
(1)通过下式计算汇总专家意见的各指标权重初值，
[0088]
[0089]
式中:
[0090]
为第a个指标权重的计算值，a＝1，2，3，4；μ
ab
为第a个指标对应于第b个专家的权重值，b＝1，2，3
…
k，其中，k为专家数量；
[0091]
(2)通过下式归一化各指标权重得到最终的权重值，
[0092][0093]
式中：
[0094]
ωa为第a个指标最终权重值。
[0095]
本实施例中，s3中根据所述冲击风险熵指标构建熵增熵减估算模型，再结合确定的冲击风险熵指标的权重值计算出熵增熵减估算模型的熵增值和熵减值，具体包括：
[0096]
s301、建立熵增熵减矩阵：
[0097]
分别构建微震频度、微震总能量、班频次变化率和班能量变化率这四个冲击风险熵指标n个时刻的熵增矩阵bz和熵减矩阵bj；
[0098][0099]
所述熵增矩阵bz中b
z,ki
表示第i个风险熵指标第k时刻的风险熵增值，i＝1，2，3，4；
[0100][0101]
所述熵减矩阵bj中b
j,ki
表示第i个风险熵指标中第k时刻的风险熵减值，j＝1，2，3，4；
[0102]
s302、计算熵增熵减值；
[0103]
u＝gu×wt
＝e
×bz
×wt
；
[0104]
v＝tv×wt
＝e
×bj
×wt
；
[0105]
其中，u表示熵增值，v表示熵减值，w表示冲击风险熵指标权重矩阵，w
t
为w的转置矩阵；e表示1
×
n的矩阵，gu表示冲击风险熵指标的熵增值之和矩阵，tv表示冲击风险熵指标的熵减值之和矩阵；
[0106]gu
＝e
×bz
＝[z1,z2,z3,z4]；
[0107]
tv＝e
×bj
＝[j1,j2,j3,j4]；
[0108]
其中，z1～z4表示四个冲击风险熵指标所有时刻的熵增值之和，j1～j4表示四个冲击风险熵指标所有时刻的熵减值之和；
[0109]
w＝[ω1,ω2,ω3,ω4]
[0110]
其中，ω1～ω4表示微震频度、微震总能量、班频次变化率和班能量变化率四个冲击风险熵指标的依次对应的权重值。
[0111]
本实施例中，将计算出熵增熵减估算模型的熵增值和熵减值带入建立的突变模型中，判断冲击地压的风险，进而预警，具体包括：
[0112]
要建立突变模型首先要判断变量数量来选择构建哪一种突变模型，根据控制变量与状态变量的数量选择突变模型中的尖点突变模型进行建立，模型如下：
[0113]
f(x)＝x4+ux2+vx
[0114]
对其求一阶导并另f'(x)＝0，得到突变模型平衡曲面方程：
[0115]
4x3+2ux+v＝0
[0116]
再次求导并另f”(x)＝0，可得控制变量平面方程：
[0117]
12x2+2u＝0
[0118]
联立f'(x)、f”(x)可得到分叉集方程，即最终判别公式δ：
[0119]
δ＝8u3+27v2；其中u表示熵增值，v表示熵减值
[0120]
以判别公式δ作为突变模型，将计算出熵增熵减估算模型的熵增值和熵减值带入判别公式δ，若δ》0，说明冲击地压的检测处于稳定状态，暂时不会突发危险；若δ《0，说明受载煤岩系统风险发生突变，有严重的冲击风险会发生，应实施有效的防治方法，做好提前的应急工作；若δ＝0，说明风险达到临界状态，随时会因为一些外在不利因素的加入导致风险发生，应保持持续监测，并适当的采用防治方法。
[0121]
以上所述，仅是本发明的较佳实施例，并非对本发明作任何限制。凡是根据发明技术实质对以上实施例所作的任何简单修改、变更以及等效变化，均仍属于本发明技术方案的保护范围内。