CN114856705A - 一种基于冲击风险熵的冲击地压预警方法 - Google Patents
一种基于冲击风险熵的冲击地压预警方法 Download PDFInfo
- Publication number
- CN114856705A CN114856705A CN202210486131.XA CN202210486131A CN114856705A CN 114856705 A CN114856705 A CN 114856705A CN 202210486131 A CN202210486131 A CN 202210486131A CN 114856705 A CN114856705 A CN 114856705A
- Authority
- CN
- China
- Prior art keywords
- entropy
- microseismic
- index
- risk
- change rate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 239000011435 rock Substances 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 title claims abstract description 35
- 230000008859 change Effects 0.000 claims abstract description 56
- 238000012544 monitoring process Methods 0.000 claims abstract description 22
- 230000003247 decreasing effect Effects 0.000 claims abstract description 18
- 230000035772 mutation Effects 0.000 claims abstract description 17
- 230000008569 process Effects 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 44
- 238000004458 analytical method Methods 0.000 claims description 7
- 239000003245 coal Substances 0.000 claims description 7
- 230000009467 reduction Effects 0.000 claims description 6
- 230000002265 prevention Effects 0.000 claims description 4
- 230000002411 adverse Effects 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
- E21F17/00—Methods or devices for use in mines or tunnels, not covered elsewhere
- E21F17/18—Special adaptations of signalling or alarm devices
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Mining & Mineral Resources (AREA)
- Life Sciences & Earth Sciences (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Bioinformatics & Computational Biology (AREA)
- Operations Research (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Probability & Statistics with Applications (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geology (AREA)
- Computing Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提出了一种基于冲击风险熵的冲击地压预警方法,通过微震监测手段获取微震频度和微震总能量,获取通过地音监测手段获取的班频次变化率和班能量变化率,将所述微震频度、微震总能量、班频次变化率和班能量变化率作为冲击风险熵指标;以专家打分表为基础,利用层次分析法确定各冲击风险熵指标的权重值;根据各冲击风险熵指标构建熵增熵减估算模型,再结合确定的冲击风险熵指标的权重值计算出熵增熵减估算模型的熵增值和熵减值;将计算出的熵增值和熵减值带入建立的突变模型中,判断冲击地压的风险,进而预警。本发明可结合突变理论模型对冲击地压中的冲击风险熵指标进行分析,以达到预警冲击地压危险的效果,从而预防此类事故的发生。
Description
技术领域
本发明属于煤矿工程事故应急预警技术领域,具体涉及一种基于冲击风险熵的冲击地压预警方法。
背景技术
随着我国煤矿开采强度和深度的不断增加,矿山动力灾害问题越来越严重,尤其是其中的冲击地压问题更是严重的影响了安全和经济问题。近年来,因为冲击地压遇难的人数逐渐剧增,因此国家越来越关注冲击地压的预警问题。很明显,冲击地压危险已成为我国威胁煤矿安全生产的重大灾害之一。因此,准确监测矿山并预警冲击地压危险是必要且必须重视的问题。
微震与地音监测技术是利用岩体受力变形和破坏后本身发射出的微震和声波来进行监测工程岩体稳定性的技术方法。因此,微震法结合地音法可以使监测预警更准确。另外,冲击地压又具有震动性、瞬时性、复杂性以及突发性,但目前针对解决冲击地压突发性的有效方法比较有限,所以通过一种冲击风险熵对冲击地压突发性进行预警的问题亟待解决。
因此,现在针对煤矿矿山的突发冲击地压问题亟需一种有效精准的预警方法来进行监测预警,从而采用应急措施及卸压防治方法进行应对处理。
发明内容
本发明所要解决的技术问题在于针对上述现有技术的不足,提供了一种基于冲击风险熵的冲击地压预警方法,该方法选择微震监测手段结合地音监测手段选定冲击风险熵指标,可结合突变理论模型对冲击地压中的冲击风险熵指标进行分析,以达到预警冲击地压危险的效果,从而预防此类事故的发生。
为解决上述技术问题,本发明采用的技术方案是:一种基于冲击风险熵的冲击地压预警方法,其特征在于,该方法包括:
S1、通过微震监测手段获取微震频度和微震总能量,获取通过地音监测手段获取的班频次变化率和班能量变化率,将所述微震频度、微震总能量、班频次变化率和班能量变化率作为冲击风险熵指标;
S2、以专家打分表为基础,利用层次分析法确定各冲击风险熵指标的权重值;
S3、根据各冲击风险熵指标构建熵增熵减估算模型,再结合确定的冲击风险熵指标的权重值计算出熵增熵减估算模型的熵增值和熵减值;
S4、将计算出的熵增值和熵减值带入建立的突变模型中,判断冲击地压的风险,进而预警;
所述突变模型的判别公式为:
Δ=8u3+27v2,其中u表示熵增值,v表示熵减值。
优选地,S2中以专家打分表为基础,利用层次分析法确定各冲击风险熵指标的权重值,具体包括:
S201、构建冲击地压预警层次分析模型,确定目标层、中间层和最底层;
S202、以层次分析模型为基础设计专家打分表,收集专家打分表,提取其相关的数据,建立微震频度、微震总能量、班频次变化率和班能量变化率的判断矩阵;
S203、求解微震频度、微震总能量、班频次变化率和班能量变化率两两判断矩阵的特征根,得出权重向量;
S204、对判断矩阵进行一致性检验,采用三个指标进行检验:一致性指标CI、一致性平均随机指标RI、相对一致性指标CR,若符合一致性检验条件,则接受该权重向量;若不符合,则调整标度值重新计算;
S205、计算微震频度、微震总能量、班频次变化率和班能量变化率的最终权重值。
优选地,S3中根据所述冲击风险熵指标构建熵增熵减估算模型,再结合确定的冲击风险熵指标的权重值计算出熵增熵减估算模型的熵增值和熵减值,具体包括:
S301、建立熵增熵减矩阵:
分别构建微震频度、微震总能量、班频次变化率和班能量变化率这四个冲击风险熵指标n个时刻的熵增矩阵Bz和熵减矩阵Bj;
所述熵增矩阵Bz中bz,ki表示第i个风险熵指标第k时刻的风险熵增值,i=1,2,3,4;
所述熵减矩阵Bj中bj,ki表示第i个风险熵指标中第k时刻的风险熵减值,j=1,2,3,4;
S302、计算熵增熵减值;
u=Gu×WT=E×Bz×WT;
v=Tv×WT=E×Bj×WT;
其中,u表示熵增值,v表示熵减值,W表示冲击风险熵指标权重矩阵,WT为W的转置矩阵;E表示1×n的矩阵,Gu表示冲击风险熵指标的熵增值之和矩阵,Tv表示冲击风险熵指标的熵减值之和矩阵;
Gu=E×Bz=[z1,z2,z3,z4];
Tv=E×Bj=[j1,j2,j3,j4];
其中,z1~z4表示四个冲击风险熵指标所有时刻的熵增值之和,j1~j4表示四个冲击风险熵指标所有时刻的熵减值之和;
W=[ω1,ω2,ω3,ω4]
其中,ω1~ω4表示微震频度、微震总能量、班频次变化率和班能量变化率四个冲击风险熵指标的依次对应的权重值。
优选地,将计算出熵增熵减估算模型的熵增值和熵减值带入建立的突变模型中,判断冲击地压的风险,进而预警,具体包括:
以判别公式Δ作为突变模型,将计算出熵增熵减估算模型的熵增值和熵减值带入判别公式Δ,若Δ>0,说明冲击地压的检测处于稳定状态,暂时不会突发危险;若Δ<0,说明受载煤岩系统风险发生突变,有严重的冲击风险会发生,应实施有效的防治方法,做好提前的应急工作;若Δ=0,说明风险达到临界状态,随时会因为一些外在不利因素的加入导致风险发生,应保持持续监测,并适当的采用防治方法。
本发明与现有技术相比具有以下优点:
1、本发明微震法结合地音法可以使监测预警更准确,能够解决冲击地压突发性,可结合突变理论模型对冲击地压中的冲击风险熵指标进行分析,以达到预警冲击地压危险的效果,从而预防此类事故的发生
下面通过附图和实施例对本发明的技术方案作进一步的详细说明。
附图说明
图1为本发明实施例1公开的一种基于冲击风险熵的冲击地压预警方法的流程示意图。
具体实施方式
实施例1
如图1所示,本发明实施例的一种基于冲击风险熵的冲击地压预警方法,该方法包括:
S1、通过微震监测手段获取微震频度和微震总能量,获取通过地音监测手段获取的班频次变化率和班能量变化率,将所述微震频度、微震总能量、班频次变化率和班能量变化率作为冲击风险熵指标;
S2、以专家打分表为基础,利用层次分析法确定各冲击风险熵指标的权重值;
S3、根据各冲击风险熵指标构建熵增熵减估算模型,再结合确定的冲击风险熵指标的权重值计算出熵增熵减估算模型的熵增值和熵减值;
S4、将计算出的熵增值和熵减值带入建立的突变模型中,判断冲击地压的风险,进而预警。
本实施例中采用微震及地音(声发射)监测预警是为了能够保证监测预警的可靠性和准确性。
本实施例中,采用微震及地音(声发射)监测预警是为了能够保证监测预警的可靠性和准确性,所述微震频度、微震总能量、班频次变化率和班能量变化率通过以下方法得到:
①微震频度Za;
经仪器测量得出数据。
②微震总能量Ze。
经仪器测量得出数据。
S202、对于地音(声发射)监测手段,有如下冲击风险熵指标:
①班频次变化率Aa;
按下式计算得出:
式中:
Ni———当天第i班的地音频次,“三八制”时,i≤3,“四六制”时,i≤4;
②班能量变化率Ae。
按下式计算得出:
式中:
Ei———当天第i班的地音能量,单位为焦耳(J),“三八制”时,i≤3,“四六制”时,i≤4;
本实施例中,S2中以专家打分表为基础,利用层次分析法确定各冲击风险熵指标的权重值,具体包括:
S201、构建冲击地压预警层次分析模型,确定目标层、中间层和最底层;
目标层:判断是否有冲击地压风险
中间层;所有冲击风险熵指标;所有可以造成冲击地压发生的因素,且可以作为冲击地压预警指标,即为冲击风险熵指标。
最底层:备选的冲击风险熵指标,备选的冲击风险熵指标指从所有的冲击风险熵中选中的上述四个(微震频度、微震总能量、班频次变化率和班能量变化率)做实验要用到的冲击风险熵指标。
S202、以层次分析模型为基础设计专家打分表,收集专家打分表,提取其相关的数据,建立微震频度、微震总能量、班频次变化率和班能量变化率的判断矩阵;
表1专家打分表
S203、求解由微震频度、微震总能量、班频次变化率和班能量变化率构造的四阶矩阵的特征根,得出权重向量;
S204、对判断矩阵进行一致性检验,采用三个指标进行检验:一致性指标CI、一致性平均随机指标RI、相对一致性指标CR,若符合一致性检验条件,则接受该权重向量;若不符合,则调整标度值重新计算;
所述一致性指标CI通过下式得出,
式中:λmax为微震频度、微震总能量、班频次变化率和班能量变化率的判断矩阵最大特征值;n为判断矩阵阶数。
表2一致性平均随机指标RI
相对一致性指标CR通过下式(7)得到,
CR=CI/RI
当CR<0.1时,认为矩阵的一致性符合要求,可以接受所计算权重向量。若CR>0.1时,认为矩阵一致性不符合要求,需要对矩阵中的值进行分析后重新给出各因素标度值,再进行验算,直至符合一致性要求;
S205、计算微震频度、微震总能量、班频次变化率和班能量变化率的最终权重值。
(1)通过下式计算汇总专家意见的各指标权重初值,
式中:
(2)通过下式归一化各指标权重得到最终的权重值,
式中:
ωa为第a个指标最终权重值。
本实施例中,S3中根据所述冲击风险熵指标构建熵增熵减估算模型,再结合确定的冲击风险熵指标的权重值计算出熵增熵减估算模型的熵增值和熵减值,具体包括:
S301、建立熵增熵减矩阵:
分别构建微震频度、微震总能量、班频次变化率和班能量变化率这四个冲击风险熵指标n个时刻的熵增矩阵Bz和熵减矩阵Bj;
所述熵增矩阵Bz中bz,ki表示第i个风险熵指标第k时刻的风险熵增值,i=1,2,3,4;
所述熵减矩阵Bj中bj,ki表示第i个风险熵指标中第k时刻的风险熵减值,j=1,2,3,4;
S302、计算熵增熵减值;
u=Gu×WT=E×Bz×WT;
v=Tv×WT=E×Bj×WT;
其中,u表示熵增值,v表示熵减值,W表示冲击风险熵指标权重矩阵,WT为W的转置矩阵;E表示1×n的矩阵,Gu表示冲击风险熵指标的熵增值之和矩阵,Tv表示冲击风险熵指标的熵减值之和矩阵;
Gu=E×Bz=[z1,z2,z3,z4];
Tv=E×Bj=[j1,j2,j3,j4];
其中,z1~z4表示四个冲击风险熵指标所有时刻的熵增值之和,j1~j4表示四个冲击风险熵指标所有时刻的熵减值之和;
W=[ω1,ω2,ω3,ω4]
其中,ω1~ω4表示微震频度、微震总能量、班频次变化率和班能量变化率四个冲击风险熵指标的依次对应的权重值。
本实施例中,将计算出熵增熵减估算模型的熵增值和熵减值带入建立的突变模型中,判断冲击地压的风险,进而预警,具体包括:
要建立突变模型首先要判断变量数量来选择构建哪一种突变模型,根据控制变量与状态变量的数量选择突变模型中的尖点突变模型进行建立,模型如下:
F(x)=x4+ux2+vx
对其求一阶导并另F'(x)=0,得到突变模型平衡曲面方程:
4x3+2ux+v=0
再次求导并另F”(x)=0,可得控制变量平面方程:
12x2+2u=0
联立F'(x)、F”(x)可得到分叉集方程,即最终判别公式Δ:
Δ=8u3+27v2;其中u表示熵增值,v表示熵减值
以判别公式Δ作为突变模型,将计算出熵增熵减估算模型的熵增值和熵减值带入判别公式Δ,若Δ>0,说明冲击地压的检测处于稳定状态,暂时不会突发危险;若Δ<0,说明受载煤岩系统风险发生突变,有严重的冲击风险会发生,应实施有效的防治方法,做好提前的应急工作;若Δ=0,说明风险达到临界状态,随时会因为一些外在不利因素的加入导致风险发生,应保持持续监测,并适当的采用防治方法。
以上所述,仅是本发明的较佳实施例,并非对本发明作任何限制。凡是根据发明技术实质对以上实施例所作的任何简单修改、变更以及等效变化,均仍属于本发明技术方案的保护范围内。
Claims (4)
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个时刻的熵增矩阵Bz和熵减矩阵Bj;
所述熵增矩阵Bz中bz,ki表示第i个风险熵指标第k时刻的风险熵增值,i=1,2,3,4;
所述熵减矩阵Bj中bj,ki表示第i个风险熵指标中第k时刻的风险熵减值,j=1,2,3,4;
S302、计算熵增熵减值;
u=Gu×WT=E×Bz×WT;
v=Tv×WT=E×Bj×WT;
其中,u表示熵增值,v表示熵减值,W表示冲击风险熵指标权重矩阵,WT为W的转置矩阵;E表示1×n的矩阵,Gu表示冲击风险熵指标的熵增值之和矩阵,Tv表示冲击风险熵指标的熵减值之和矩阵;
Gu=E×Bz=[z1,z2,z3,z4];
Tv=E×Bj=[j1,j2,j3,j4];
其中,z1~z4表示四个冲击风险熵指标所有时刻的熵增值之和,j1~j4表示四个冲击风险熵指标所有时刻的熵减值之和;
W=[ω1,ω2,ω3,ω4]
其中,ω1~ω4表示微震频度、微震总能量、班频次变化率和班能量变化率四个冲击风险熵指标的依次对应的权重值。
4.根据权利要求1所述的一种基于冲击风险熵的冲击地压预警方法,其特征在于,将计算出熵增熵减估算模型的熵增值和熵减值带入建立的突变模型中,判断冲击地压的风险,进而预警,具体包括:
以判别公式Δ作为突变模型,将计算出熵增熵减估算模型的熵增值和熵减值带入判别公式Δ,若Δ>0,说明冲击地压的检测处于稳定状态,暂时不会突发危险;若Δ<0,说明受载煤岩系统风险发生突变,有严重的冲击风险会发生,应实施有效的防治方法,做好提前的应急工作;若Δ=0,说明风险达到临界状态,随时会因为一些外在不利因素的加入导致风险发生,应保持持续监测,并适当的采用防治方法。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210486131.XA CN114856705A (zh) | 2022-05-06 | 2022-05-06 | 一种基于冲击风险熵的冲击地压预警方法 |
ZA2022/11205A ZA202211205B (en) | 2022-05-06 | 2022-10-13 | Method for forewarning rock burst based on impact risk entropy |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210486131.XA CN114856705A (zh) | 2022-05-06 | 2022-05-06 | 一种基于冲击风险熵的冲击地压预警方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114856705A true CN114856705A (zh) | 2022-08-05 |
Family
ID=82635429
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210486131.XA Pending CN114856705A (zh) | 2022-05-06 | 2022-05-06 | 一种基于冲击风险熵的冲击地压预警方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN114856705A (zh) |
ZA (1) | ZA202211205B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115390130A (zh) * | 2022-08-29 | 2022-11-25 | 吉林建筑大学 | 一种煤矿开采大能量微震事件预测方法和装置 |
-
2022
- 2022-05-06 CN CN202210486131.XA patent/CN114856705A/zh active Pending
- 2022-10-13 ZA ZA2022/11205A patent/ZA202211205B/en unknown
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115390130A (zh) * | 2022-08-29 | 2022-11-25 | 吉林建筑大学 | 一种煤矿开采大能量微震事件预测方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
ZA202211205B (en) | 2022-10-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sun et al. | Correlational fractal characterisation of stress and acoustic emission during coal and rock failure under multilevel dynamic loading | |
CN103291364B (zh) | 一种冲击矿压的微震多维信息综合时序预警方法 | |
CN114810211B (zh) | 一种基于矿震群震动波能量衰减特性的冲击地压危险预测方法 | |
CN107563092A (zh) | 一种矿山动力灾害的全息预警方法 | |
CN105930624B (zh) | 一种用于挡土墙结构的稳定性预警方法 | |
CN113958366B (zh) | 基于震动-应力双场监测的冲击危险性动态量化预警方法 | |
CN111963243A (zh) | 一种基于动静组合应力分析的冲击地压危险监测预警方法 | |
CN112307601B (zh) | 一种复杂储层可压裂性评价方法 | |
CN114856705A (zh) | 一种基于冲击风险熵的冲击地压预警方法 | |
CN106202737B (zh) | 一种致密储层体积压裂适应区筛选方法及装置 | |
CN112031872A (zh) | 一种冲击地压全息预警方法及装置、存储介质、计算设备 | |
CN113821977A (zh) | 一种用于tbm隧道施工的岩爆风险评估系统及方法 | |
CN111538071B (zh) | 陡倾层状岩体洞室群高边墙位移定量预测方法 | |
CN114294062B (zh) | 一种冲击地压时空动态综合预警方法 | |
CN114662327A (zh) | 一种用于auv舱内泄露监测的辅助方法 | |
Kia et al. | Assessment the effective ground motion parameters on seismic performance of R/C buildings using artificial neural network | |
CN114705549A (zh) | 一种易泥化矿岩失稳突变的预警方法 | |
CN110008520B (zh) | 基于位移响应协方差参数和贝叶斯融合的结构损伤识别方法 | |
CN112883646B (zh) | 联合机器学习和土力学模型的建筑沉降量提取方法、系统及装置 | |
CN112329255A (zh) | 基于倾向度和未确知测度的岩爆预测方法 | |
Byrne et al. | Seismic liquefaction: centrifuge and numerical modeling | |
CN103294926B (zh) | 一种冲击矿压时序监测预警方法 | |
CN115266348B (zh) | 一种基于动静载叠加试验的岩爆倾向性等级判断方法 | |
CN114063152B (zh) | 一种基于矿震统计特征的冲击地压主控因素确定方法 | |
CN112417374B (zh) | 大采高综采工作面煤壁片帮的预测方法和装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |