CN109918698B - 确定冲击地压临界能量的方法、装置、设备和存储介质 - Google Patents

Abstract

本申请公开了一种确定冲击地压临界能量的方法、装置、设备和存储介质。方法包括以下步骤：基于煤体损伤本构模型得到煤体破坏损失能量，基于所述煤体破坏损失能量和破碎煤体临界抛出能量，确定发生冲击地压的临界能量。

Description

确定冲击地压临界能量的方法、装置、设备和存储介质

技术领域

本公开一般涉及地质灾害防治技术领域，尤其涉及一种确定冲击地压临界能量的方法、装置、设备和存储介质。

背景技术

随着煤矿开采深度和开采强度的越来越大，矿井动力灾害发生的频次也随之增多，对于冲击地压矿井，冲击地压动力灾害是煤矿安全生产的重中之重。冲击地压是在煤岩体能量达到或超过冲击地压临界能量时，煤岩体聚积的弹性能突然释放，将煤岩体抛出，伴随着强烈声响，对井下设备、井下工程空间造成破坏，甚至对井下工作人员造成伤亡。此外，冲击地压发生还会对其他矿井灾害产生影响。由于冲击地压发生机理复杂，又受地质条件和开采条件影响较大，因而对冲击地压研究是保证煤矿安全生产的重要课题之一。

目前国内外学术界关于冲击地压发生机理研究主要有刚度理论、强度理论、冲击倾向理论，随着对冲击地压发生机理研究更加深入，又进一步提出了三准则理论、剪切滑移理论、″三因素″理论和变形失稳理论等，对煤矿冲击地压防治和保障安全生产起到了举足轻重的作用。(1)刚度理论：压力试验机出现以后，Petukhov和Cook及Digest最先发现冲击地压可以用岩石试件在刚度较小的柔性试验机上的动态破坏来描述，进而提出了冲击地压刚度理论。Cook在提出能量理论之后，又认识到冲击地压是由于结构失稳造成的。又利用大理岩进一步进行实验。Salaman、Brady和Petukhov等也对多个矿柱的冲击地压问题进行大量的研究。(2)强度理论：强度理论观点为当煤岩体承受的载荷到达其强度极限时，煤岩体就会开始破坏。夹持煤体经典理论是由布霍依诺提出的。若煤体-围岩交界面处以及煤体本身达到极限平衡条件，就达到了冲击地压发生的强度条件。但冲击地压还和煤岩体突然破坏有关，有时煤岩破坏并不意味着冲击地压一定发生。(3)冲击倾向理论：冲击倾向理论指出：若煤岩体的冲击倾向度K_E≥K_EC时，则煤岩体就可能发生冲击破坏。(4)三准则理论：李玉生指出强度准则是煤体破坏准则，而能量准则和冲击倾向准则是突然破坏准则，提出了冲击地压发生机理可由强度准则、能量准则及冲击倾向性准则来描述，当三个准则全部满足时，将发生冲击地压。(5)“三因素”理论：齐庆新等认为冲击地压发生受三个因素影响，煤岩体内在因素、开采工程活动力源因素和结构因素。据此提出了冲击地压发生″三因素″理论。(6)变形失稳理论：章梦涛认为冲击地压是煤岩体在进入了峰值强度后出现的应***化现象，基于有限元法，建立了数学模型，对冲击地压进行定量研究。

综上所述，对于冲击地压机理研究国内外学者都从不同的角度进行分析论证，针对不同的地质动力条件也给出了***的描述，取得了很多宝贵的成果与经验。

目前，关于煤岩体冲击破坏的能量理论分析研究较多，但这些能量理论大多数都是从同一出发点进行研究的，即：煤岩体冲击破坏过程时所释放的能量与其破坏所消耗的各种能量之间的大小关系，煤岩体破坏后剩余的能量为冲击地压发生提供了能量。到目前为止现有技术还不能全面揭示冲击地压发生规律和机理，在矿井施工过程中，即使对地应力进行监测，由于无法给出发生冲击地压的相关依据，因此无法实现对是否会发生冲击地压的相关预测。

发明内容

鉴于现有技术中的上述缺陷或不足，期望提供能预测是否会发生冲击地压的一种确定冲击地压临界能量的方案。

第一方面，本申请实施例提供了一种确定冲击地压临界能量的方法，包括以下步骤：

基于煤体损伤本构模型得到煤体破坏损失能量，基于所述煤体破坏损失能量和破碎煤体临界抛出能量，确定发生冲击地压的临界能量。

所述发生冲击地压的临界能量等于所述煤体破坏损失能量和所述破碎煤体临界抛出能量之和。

通过相应的能量密度来表达相应的所述煤体破坏损失能量和所述破碎煤体临界抛出能量，可以通过相应的理论公式计算能量密度，进而反映相应的煤体破坏损失能量和破碎煤体抛出能量。

所述煤体损伤本构模型的公式为：

Y为损伤耗能率；

σ为应力；

E为弹性模量；

D为损伤变量，表示单位体积微裂纹数；

D_E为在位置E完全断裂破坏时的损伤值；

D_D为在位置D完全断裂破坏时的损伤值；

U_D为煤岩体损伤消耗能量。

所述煤体损伤本构模型中的损伤变量通过所述矿井中煤体的应力-应变试验关系确定。其中，建立煤岩体损伤本构模型拟合曲线、损伤变量曲线关系、使用Matlab软件计算煤矿煤岩体破坏损失的能量为本领域的现有技术。根据煤岩体微裂纹的扩展、内部损伤的发展演化及分布密切相关性，建立煤岩体损伤模型，建立模型时可以认为：

(1)煤岩体损伤是在应力主轴的方向上，弹性与损伤的耦合下进行的，视为损伤张量主对角矩阵；

(2)煤岩体损伤演化应该与应力或应变存在幂函数关系，且损伤变量可以通过内拉应变宏观物理量来反映。

所述破碎煤体临界抛出能量采用以下公式：

U_F：破碎煤体临界抛出能量；

ρ：破碎煤体平均密度；

v：破碎煤体向自由空间抛出的平均初速度。

所述破碎煤体临界抛出能量所对应的破碎煤体的抛出速度为8-12m/s。当抛出速度v≤1m/s时不发生冲击地压，而当v≥10m/s时，则冲击地压发生具有较高的可能性。一般情况下，可以将抛出速度设定为8-12m/s，来计算临界能量公式。

第二方面，本申请实施例提供了一种确定冲击地压临界能量的装置，所述装置包括：

获取模块，用于获取煤体损伤本构模型和破碎煤体临界抛出能量的相关参数；

计算模块，用于基于所述煤体损伤本构模型得到煤体破坏损失能量，基于所述煤体破坏损失能量和所述破碎煤体临界抛出能量，确定发生冲击地压的临界能量。

第三方面，本申请实施例提供了一种确定冲击地压临界能量的设备，包括：

一个或多个处理器；

存储器，用于存储一个或多个程序；

当所述一个或多个程序被所述一个或多个处理器执行时，使得所述一个或多个处理器执行任一项所述的确定冲击地压临界能量的方法。

第四方面，本申请实施例提供了一种存储有计算机程序的计算机可读存储介质，该程序被处理器执行时，实现任一项所述的确定冲击地压临界能量的方法。

本申请实施例提供的确定冲击地压临界能量的方案，基于冲击地压能量理论，从冲击地压发生的根源入手，研究出了冲击地压与煤岩体能量的关系，提出原岩能量为冲击地压发生的主要能量、决定了发生冲击地压的临界能量条件。损伤力学能够根据煤岩体的岩性，得到媒体破坏时所损失的能量；而破碎媒体抛出能量则根据媒体被抛出时速度的不同，对应不同的抛出能量及其冲击大小，进而可以预测不同冲击地压的强度。当抛出速度较小时，不会发生冲击地压；当抛出速度达到临界抛出速度时，则会发生冲击地压。该方法能全面揭示冲击地压发生规律和机理，可以定量的给出发生冲击地压的临界能量，在矿井施工过程中，依据该冲击地压临界能量，可以实现对是否会发生冲击地压的相关预测，在矿井施工期做好冲击地压防治的前瞻性工作，从根源上防治冲击地压。

附图说明

通过阅读参照以下附图所作的对非限制性实施例所作的详细描述，本申请的其它特征、目的和优点将会变得更明显：

图1示出了本申请实施例中工作面前方弹性潜能变化示意图；

图2示出了煤岩体变形破坏应力-应变曲线图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整的描述。显然，所描述的实施例是本发明一部分而不是全部的实施例。

需要说明的是，在不冲突的情况下，通常在此附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。可以理解的是，此处所描述的具体实施例仅仅用于解释相关发明，而非对该发明的限定。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

如背景技术中所说，目前，关于煤岩体冲击破坏的能量理论分析研究较多，但这些能量理论大多数都是从同一出发点进行研究的，即：煤岩体冲击破坏过程时所释放的能量与其破坏所消耗的各种能量之间的大小关系，煤岩体破坏后剩余的能量为冲击地压发生提供了能量。然而，现有技术中，并没有对冲击地压孕育、发展、发生过程中能量之间关系的相关研究，尤其并未考虑到原岩应力对冲击地压的影响。

本申请的理论研究成果如下：

矿井不具备发生矿井冲击地压的区域地质动力条件，与开采活动引起的采动应力相互耦合，从而导致应力的增加和能量的积聚，达到发生矿井冲击地压的条件，从而诱发矿井冲击地压。在这种动态平衡过程中，能量首先集中在煤壁附近，当能量值达到煤体的软化强度后，煤体发生软化变形，使集中能量向煤体深部转移，达到新的能量平衡后。当集中能量值超过煤体的损伤耗散能时，煤体发生破坏，可能发生冲击地压。

工作面回采过程中，能量聚集区逐步推进，导致煤岩体破碎，能量释放。定义煤岩体在初始地应力场作用下积聚的能量为原岩能量，用U_B表示；煤岩体损伤消耗的能量为U_D；冲击地压发生的临界能量为U_L；开采条件下，随工作面推进煤岩体总能量为U_T。

在工作面回采过程中，煤岩体内总能量U_T是一个不断变化的变量。煤岩体总能量U_T与煤岩体损伤消耗的能量为U_D之差为U_T-D，它的变化可能会出现如下的三种可能性：

①当U_T-D＞0，说明煤岩体内积聚的能量多于煤岩体损伤消耗的能量，U_T-D随工作面推进不断增加；

②当U_T-D＝0，说明煤岩体内积聚的能量与煤岩体损伤消耗的能量相等，能量的积聚与耗散处于平衡状态；

③当U_T-D＜0，说明煤岩体损伤消耗的能量多于岩体内积聚的能量，U_T-D随工作面推进不断降低。

因此，工作面前方煤岩体内能量的变化可由附图1所示模型来表示，其中U_T-D为U_T(煤岩体总能量)与U_D(煤岩体损伤消耗的能量)之差，U_L为冲击地压的临界能量。

当U_T-D＞0时，说明煤岩体内所积聚的能量大于耗散的能量，能量处于不断增加的状态。但只要U_T-D＜U_L，都不会出现冲击地压。U_T-D-U_L之差越大，危险性越大。当U_T-D＝U_L时，如再有能量增加，就会发生冲击地压。如果U_T-D＞U_L，而且不及时采取解危措施来释放能量，则在工作面回采过程中，其能量进一步增加，就可能发生冲击地压。

当U_T-D＜0时，说明煤岩体内所积聚的能量小于煤岩体内耗散的能量，煤岩体能量逐渐释放，能量不断减少，U_T-D-U_L之差越来越大，危险性也就越来越小，这种情况下，不采用任何解危措施，也不会出现冲击地压危险。

Grady和Kipp(1980)提出了岩体内部单位体积微裂纹数N与应变ε在单轴拉伸应力状态下，满足双参数(k和m)Weibull分布规律，如公式(1)。

根据煤岩体微裂纹的扩展、内部损伤的发展演化及分布密切相关性，建立煤岩体损伤模型，建立模型时可以认为：

(1)煤岩体损伤是在应力主轴的方向上，弹性与损伤的耦合下进行的，视为损伤张量主对角矩阵；

(2)煤岩体损伤演化应该与应力或应变存在幂函数关系，且损伤变量可以通过内拉应变宏观物理量来反映；

一维煤岩体损伤模型建立时假定：岩石在单轴压力状态下，张应变引起内部微裂纹扩展，因此建立微裂纹数与张应变之间的关系，用公式(2)表示。

N＝∝mεⁿ (2)

式中：N-单位体积的微裂纹数；

ε-轴向压应变，在单轴作用下，张应变与压应变数值相等；

m，n-煤岩体参数。

用损伤变量D表示单位体积微裂纹数，用公式(3)表示。

式中：

-有效应力；σ-应力。

设：

则：

损伤煤岩体本构关系公式(7)。

E为弹性模量。

故煤岩体损伤变量的演化方程为公式(8)。

受损煤岩体弹性应变能释放率dΦ_e和损伤耗能率Y可以表示为公式(9)、公式(10)。

根据有效应力概念和应变等价原理，可得受损煤岩体损伤释放率公式(11)。

对于各向同性损伤，损伤变量为标量。

则：

设：

则偏应力张量为：S_ij＝σ-σ_HI (14)

Misses等价应力为公式(15)。

由公式(12)得到公式(16)。

式中：S_t-三轴应力因子，反映了三轴应力比的影响。

其中

在单轴应力下，

定义损伤等价应力为：

则：

煤岩体随着外载的增加，损伤将不断发展，直到岩石完全断裂，从理论上讲当D＝1时，认为煤岩体完全断裂破坏，但很多试验表明，当D＜1时材料已完全破坏，如图2位置E，设在位置E完全断裂破坏时的损伤值为D_E，则对应的损伤能量释放率为Y_E，位置D煤岩体内部损伤迅速增长，导致应***化加剧，煤岩体内部开始破坏，则煤岩体损伤消耗的能量为公式(21)。

下述实施例中，木城涧煤矿煤体损伤消耗的能量通过推导公式计算，根据木城涧煤矿煤体应力-应变试验曲线采用Origin9.0软件分段拟合了峰值前后的应力-应变关系式。

通过应力-应变曲线建立与D的关系如下式所示。

根据应力-应变曲线，当应变为0.0117时为煤体破坏的起点，煤体内部损伤迅速增长，煤体内部开始破坏，当应变为0.014时为煤体已完全破坏，故按公式(23)，使用Matlab软件计算木城涧煤矿煤体破坏损失能量为1.45×10⁵J/m³。

冲击地压发生时，一方面煤体损伤消耗的能量U_D，另一方面破碎煤体抛出能量U_F，煤体损伤消耗的能量与破碎煤体抛出能量之和定义为冲击地压临界能量，用U_L表示，计算公式(25)。通常物体运动以质点的速度作为衡量标准，破碎煤体向自由空间抛出的的平均初速度(v)决定了煤体产生的动能。当v≤1m/s时不发生冲击地压，而当v≥10m/s时，则冲击地压发生具有较高的可能性。通过临界能量计算公式(25)，可得能量密度为2.41×10⁵J/m³。

U_L＝U_D+U_F (25)

本申请中，确定冲击地压临界能量的目的，即为将U_T-D与冲击地压临界能量值进行对比，当U_T-D＜冲击地压临界能量时，可以确定为无冲击地压的危险；当U_T-D＞冲击地压临界能量时，可以确定存在冲击地压的危险。可以理解的是，由于冲击地压临界能量的确定与破碎煤体临界抛出能量有关，在某些要求非常严格的情况下，不允许存在有破碎煤体抛出的情况，此时可以将破碎煤体临界抛出能量设为零。而在本申请的下述实施例中，将临界抛出能量设为10m/s，下例所说的无冲击地压和弱冲击地压，均理解为不存在冲击地压的风险；而中等冲击地压和强冲击地压，理解为属于存在冲击地压风险的情况。

本申请实施例提供的一种确定冲击地压临界能量的方法，包括以下步骤：

基于煤体损伤本构模型得到煤体破坏损失能量，基于所述煤体破坏损失能量和破碎煤体临界抛出能量，确定发生冲击地压的临界能量。

所述发生冲击地压的临界能量等于所述煤体破坏损失能量和所述破碎煤体临界抛出能量之和。

通过相应的能量密度来表达相应的所述煤体破坏损失能量和所述破碎煤体临界抛出能量。

所述煤体损伤本构模型的公式为：

Y为损伤耗能率；

σ为应力；

E为弹性模量；

D为损伤变量，表示单位体积微裂纹数；

D_E为在位置E完全断裂破坏时的损伤值；

D_D为在位置D完全断裂破坏时的损伤值；

U_D为煤岩体损伤消耗能量。

所述煤体损伤本构模型中的损伤变量通过所述矿井中煤体的应力-应变试验关系确定。

所述破碎煤体临界抛出能量采用以下公式：

U_F：破碎煤体临界抛出能量；

ρ：破碎煤体平均密度；

v：破碎煤体向自由空间抛出的平均初速度。

所述破碎煤体临界抛出能量所对应的破碎煤体的抛出速度为8-12m/s。

第二方面，本申请实施例提供了一种确定冲击地压临界能量的装置，所述装置包括：

获取模块，用于获取煤体损伤本构模型和破碎煤体临界抛出能量的相关参数；

计算模块，用于基于所述煤体损伤本构模型得到煤体破坏损失能量，基于所述煤体破坏损失能量和所述破碎煤体临界抛出能量，确定发生冲击地压的临界能量。

第三方面，本申请实施例提供了一种确定冲击地压临界能量的设备，包括；

一个或多个处理器；

存储器，用于存储一个或多个程序；

当所述一个或多个程序被所述一个或多个处理器执行时，使得所述一个或多个处理器执行任一项所述的确定冲击地压临界能量的方法。

第四方面，本申请实施例提供了一种存储有计算机程序的计算机可读存储介质，该程序被处理器执行时，实现任一项所述的确定冲击地压临界能量的方法。

实施例1

1区域构造特征及能量场特征

1.1煤田区域构造及其活动性特征

(1)冲击地压能量来源及煤岩体破坏能量研究

宏观分析了地质构造、新构造运动和现代构造应力对冲击地压能量孕育作用，利用地质动力区划方法对冲击地压危险性进行评价；微观分析了煤岩体变形破坏时能量转移与耗散情况，基于冲击地压现象形成的特点，确定冲击地压发生临界能量密度。

(2)井田构造运动与原岩能量分析

基于京西煤田构造应力场、能量场及地壳应变能特征，分析了木城涧井田内构造形态及运动方式对冲击地压的控制作用，利用地质区划方法划分活动断裂，分析其对冲击地压控制作用，计算出自然地质条件下的原岩能量，基于冲击地压临界能量密度对能量区划分，划分三槽、二槽煤层的高能量区和释放区，分析高能量区对冲击地压的控制作用，进而建立原岩能量与冲击地压强度关系。

(3)开采条件下采动能量分析

基于自然条件下的原岩能量分别特征，应用FLAC3D数值模拟软件，对三槽、二槽煤层原岩能量积聚区工作面进行模拟开采，分析开采条件下采动活动对工作面前方、煤柱区能量分布特征的影响，对冲击地压强度进行分类。基于工作面前方、煤柱区能量分布特征和冲击地压强度分类，进而建立在开采条件下煤岩体能量与冲击地压强度的关系，针对不同冲击地压强度区域采取分等级的防治技术措施。

2构造应力场与能量场特征

课题组近年来通过对京西矿区木城涧煤矿、大安山煤矿、大台井、长沟峪煤矿等课题的研究，地应力测量及能量数值结果表明(表3.1)，京西地区最大主应力方位为北东-东西，能量密度最大值为3.25×10⁵J/m³，能量密度最小值为0.44×10⁵J/m3，能量密度平均值为1.51×10⁵J/m³。研究表明构造运动对京西煤田能量分布影响较大，能量大小不同，由于构造运动导致有些区域能量积聚，而有些区域能量释放。通过赵德安、景锋等统计分析了我国大量地应力测量数据，根据地应力算出能量密度，得到根据式计算得到的京西煤田能量密度与全国能量密度的对比图，表明京西煤田能量密度高于全国平均水平，特别是在埋深700m以上，京西煤田的能量密度更加显著，京西煤田能量密度分布不均匀，表明京西煤田煤岩体承受更为强烈的非均匀能量环境，据此将对木城涧井田区域进行能量区划分。

3井田位置与构造概况

3.1井田位置及交通

木城涧井田位于北京西部山区，隶属于京西煤田，矿区地跨房山、门头沟两个区。

3.2煤系地层及顶底板岩性

(1)煤系地层

侏罗纪含煤地层为窑坡组，最大厚度可超过720m。沉积物岩性、岩相变化大，旋回复杂，由砂岩、粉砂岩、泥质岩、煤层、夹砾岩、凝灰质砂岩等组成。

(2)煤层顶底板

井田内煤层顶板裂隙及小断层发育，顶板较破碎，顶底板起伏不平，对支护要求较高。

3.3井田地质构造特征

木城涧侏罗纪井田地处京西庙安岭-髫髻山向斜南西段南东翼。井田内发育一系列次级褶皱构造，主要轴迹方位呈NE56°～70°、NE40°～45°、NE5°～15°、NE15°～30°、NW330°～350°的褶皱。井田内褶曲构造对冲击地压能量的积聚起到重要控制作用。

木城涧煤矿井田内原岩能量受地质构造的控制，不同类型、不同时期所形成的地质构造对原岩能量的影响不同，同一构造的不同位置煤岩体能量分布情况也有很大差异，分析木城涧煤矿井田构造对原岩能量分布影响程度，为下面计算井田能量区能量分布至关重要。

在构造运动的作用下，断层的端部积聚大量的弹性能，为冲击地压的发生提供动力条件。当该区域受到采掘工作的扰动时，煤岩体内的应力和积聚的弹性能急剧释放，破坏煤体，在能量膨胀的作用下，将煤体抛向采掘空间，发生冲击地压。对木城涧煤矿发生的冲击地压事故统计分析表明，在木城涧煤矿近五年发生的20次冲击地压事故中，几乎全部位于断层附近，或受到井田内断层的影响。

通过以上分析，木城涧煤矿煤体受褶曲构造、断裂构造、构造复合、联合、组合等地质构造的控制，矿井冲击地压多发生在地质构造带。3次冲击地压在牛道沟向斜附近、17次冲击地压在北港沟逆断层和盲沟同沉积正断层附近。

3.4京西地区断块划分

地质动力区划工作应该遵循从总体到局部的原则，通过逐步缩小划分范围，分析井田尺度或者更小的范围的断块构造，建立板块构造与原岩能量之间的联系，岩体能量在井田划分断块图为背景下进行计算，因此，分析井田活动断裂是计算原岩能量的基础。

①I级断块构造划分

在1∶250万比例的地形图上查明I级断块构造，I级断裂共14条。

②II级断块构造划分

在1∶100万比例的地形图上查明的II级断块构造，在研究区域内共划出活动断裂40条。

③III级断块构造划分

在1∶20万的地形图上，以木城涧煤矿为中心，查明III级断块构造，划分出的III级断裂共有12条。

④IV级断块构造划分

在1∶5万比例的地形图上查明的IV级断块构造，在研究区域内共划出活动断裂29条。

⑤V级断块构造划分

在1∶1万比例的地形图上查明的V级断块构造，在研究区域内共划出活动断裂21条。

4自然条件下原岩能量区划分

4.1井田地应力分布规律

应力场的研究对了解构造活动过程有明显的重要性。根据已知断裂来估算原岩体能量分布。基于板块学说的地质动力区划方法在分析原岩能量方面具有显著代表性。区域构造和岩体应力状态间的内在关系是通过理论或数值分析、反演、回算和模拟来揭示的。

结合木城涧煤矿应力异常发生的实际区域，以及井下实际的生产技术条件。最终确定地应力测试区域定在+450m水平2-4石门的3槽围岩区域，+450m水平一共布置了3个现场测试点。

针对木城涧煤矿进行地应力测试，得到3个地应力测孔计算的最大主应力值、中间主应力值、最小主应力值以及方位角与倾角等结参数。

4.2自然条件下原岩应力计算及划分

(1)顶板岩性分析及划分

通过对木城涧煤矿钻孔资料进行分析处理，最终得到木城涧煤矿三槽和二槽煤层顶板的岩性分布图。

(2)建立计算模型、网格划分

通过I～V级断块图构建木城涧煤矿区域现代构造运动的格架，选择了V级断块图来形成模型。

(3)参数定义

需要定义的参数主要包括以下几种：地应力、岩体和断裂力学参数、岩体力学参数、断裂几何参数等。地应力测试结果获得上述参数。岩性分布依据井田钻孔数据确定。

4.3能量密度计算及能量区图形输出

根据计算所得应力数据，反算能量密度值，用等值线图方式将能量密度显示出来。

4.4原岩能量区划分及其对冲击地压影响分析

自然地质条件下能量区划分原则依据冲击地压发生临界能量密度条件与原岩能量密度比较划分高能量区、低能量区，超过冲击地压临界能量密度条件的原岩能量密度的区域为高能量区，低能量区是按无冲击地压能量密度条件比原岩能量密度小的区域称为低能量区。

(1)煤层能量区划分

在岩体能量密度计算的基础上，将木城涧煤矿三槽煤层顶板能量密度划分为高能量区、低能量区。20次冲击地压12次在高能量区，水平最大主应力值在27.5MPa～31.5MPa，能量密度值在1.67×10⁵J/m³～2.77×10⁵J/m³，1次在水平最大主应力值为26MPa、能量密度值为2.37×10⁵J/m³，1次在水平最大主应力值为25MPa、能量密度值为2.39×10⁵J/m³，1次在水平最大主应力值为24MPa、能量密度值为2.43×10⁵J/m³，3次在水平最大主应力值为23MPa、能量密度值为2.53×10⁵J/m³。

5开采条件下采动能量分析

5.1原岩能量积聚区采动能量密度分析

5.1.1三槽西五壁模型建立

为了全面***地反映+250m水平西一石门三槽西五壁回采工作面在回采过程中工作面前方和煤柱内的能量密度分布规律，以地质条件和开采技术条件为背景，利用FLAC3D大数值模拟软件建立模型，模型长、宽、高分别为375m、300m、203m。

计算模型边界条件为X轴、Y轴两端边界分别施加沿X轴、Y轴的约束，X方向、Y轴位移都为零，模型顶部为自由边界，底部边界固定。计算模型边界载荷条件为X轴、Y轴、Z轴方向分别施加19.9MPa、45.0MPa、22.6MPa的梯度应力。

5.1.2三槽西五壁采动能量密度分析

计算得到三槽西五壁工作面回采30m时能量密度分布，具体得到最大能量密度，最小能量密度。沿工作面走向中部位置能量密度分布曲线等。工作面中部能量密度峰值为5.48×10⁵J/m³，属于中等冲击地压强度，能量密度峰值距工作面煤壁距离为5.3m。

5.2煤岩体能量特征与冲击地压强度关系

5.2.1冲击地压强度分类

(1)无冲击地压危险性

在工作面开采过程中，煤体的能量小于煤体损伤消耗的能量时，不会发生冲击地压，根据煤体损伤本构模型计算木城涧煤矿煤体破坏损失能量为1.45×10⁵J/m³，因此，当煤体在开采条件下，能量密度小于1.45×10⁵J/m³，确定为无冲击地压危险。

(2)弱冲击地压危险性界定值确定

依据弱冲击地压的现象描述，煤或岩石向已采空间抛出，但破坏性不是很大，对支架、机器和设备基本上没有破坏；围岩产生震动，伴有很大声响；产生煤尘，在瓦斯煤层中可能有大量瓦斯涌出，破碎煤体向自由空间抛出的平均初速度v＜10m/s。通过计算公式

可得当抛出的的平均初速度v＝1Om/s时，能量密度为2.41×10⁵J/m³。因此，当煤体在开采条件下，能量密度在1.45×10⁵J/m³～2.41×10⁵J/m³之间时，确定为弱冲击地压危险。

(3)中等冲击地压强度界定值确定

依据中等冲击地压的现象描述，部分煤或岩石急剧破碎，大量向已采空间抛出，破碎煤体向自由空间抛出的平均初速度v≥10m/s。

由冲击地压***的能量与震级的关系式，可以确定面波震级M_L＝2时，冲击地压***的尺度半径R计算公式：

得出冲击地压***尺度半径为1.83m，通过微震能量在10⁸J时，通过公式5.3计算得出煤体能量密度。下面针对木城涧煤矿2008年5月三槽东一壁已发生的冲击地压实际情况，依据冲击地压***尺度半径的计算公式(5.2)，其中，σ₁＝28.7MPa，σ₂＝23.1MPa，σ₃＝15.9MPa，E＝1828MPa，μ＝0.25，γ＝40000KN/m3，H＝500，得出冲击地压***尺度半径为1.83m，通过微震能量在10⁸J时，通过公式5.3计算得出煤体能量密度为1.3×10⁶J/m³。

因此，当煤体在开采条件下，能量密度在2.41×10⁵J/m³～1.3×10⁶J/m³之间时，确定为中等冲击地压危险。

R＝1.83m

(4)强冲击地压强度界定值确定

依据强冲击地压的现象描述，大部分煤或岩石急剧破碎，大量向已采空间抛出，出现支架折损、设备移动和围岩震动，面波震级M_L在2级以上，伴有巨大声响，形成大量煤尘和产生冲击波，微震能量在10⁸J以上。因此，当煤体在开采条件下，能量密度大于1.3×10⁶J/m3之间时，确定为强冲击地压危险。

采动条件下煤岩体能量特征与冲击地压强度关系

工作面开采工程活动，引起围岩、上覆岩层力学性质及赋存状态变化，破坏了自然状态下的应力平衡，导致应力状态重新分布，达到新的平衡状态。自然地质动力条件下原岩能量与开采活动引起的采动能量相互耦合作用，导致能量的积聚，达到发生矿井冲击地压的能量条件，从而诱发矿井冲击地压。

木城涧井田目前工程活动已进入了冲击地压危险强度区域。在地质及工程允许条件下，选择区域性、局部解危措施来降低危险程度，目的是使高能量区煤岩体能量得以释放，达到解除或降低危险程度。在工作面回采的同时，采取相应的局部检测措施，一方面是为了降低开采活动区域的危险程度；另一方面为检验解危措施的有效性。根据相应的局部检测结果，决定是否继续进行开采工程活动或是否继续采取冲击地压解危措施。目前木城涧煤矿采取的解危措施有钻孔卸压、煤层注水、卸压***等措施。

以上描述仅为本申请的较佳实施例以及对所运用技术原理的说明。本领域技术人员应当理解，本申请中所涉及的发明范围，并不限于上述技术特征的特定组合而成的技术方案，同时也应涵盖在不脱离所述发明构思的情况下，由上述技术特征或其等同特征进行任意组合而形成的其它技术方案。例如上述特征与本申请中公开的(但不限于)具有类似功能的技术特征进行互相替换而形成的技术方案。

Claims (5)

1.一种确定冲击地压临界能量的方法，其特征在于，包括以下步骤：

基于煤体损伤本构模型得到煤体破坏损失能量，基于所述煤体破坏损失能量和破碎煤体临界抛出能量，确定发生冲击地压的临界能量；

所述发生冲击地压的临界能量等于所述煤体破坏损失能量和所述破碎煤体临界抛出能量之和；

通过相应的能量密度来表达相应的所述煤体破坏损失能量和所述破碎煤体临界抛出能量；

所述煤体损伤本构模型的公式为：

Y为损伤耗能率；

σ为应力；

E为弹性模量；

D为损伤变量，表示单位体积微裂纹数；

DE为在位置E完全断裂破坏时的损伤值；

DD为在位置D完全断裂破坏时的损伤值；

UD为煤岩体损伤消耗能量；

所述煤体损伤本构模型中的损伤变量通过矿井中煤体的应力-应变试验关系确定；

所述破碎煤体临界抛出能量采用以下公式：

UF：破碎煤体临界抛出能量；

ρ：破碎煤体平均密度；

ν：破碎煤体向自由空间抛出的平均初速度。

2.根据权利要求1所述的确定冲击地压临界能量的方法，其特征在于，所述破碎煤体临界抛出能量所对应的破碎煤体的抛出速度为8-12m/s。

3.一种确定冲击地压临界能量的装置，其特征在于，所述装置包括：

获取模块，用于获取煤体损伤本构模型和破碎煤体临界抛出能量的相关参数；

计算模块，用于基于所述煤体损伤本构模型得到煤体破坏损失能量，基于所述煤体破坏损失能量和所述破碎煤体临界抛出能量，确定发生冲击地压的临界能量；

所述发生冲击地压的临界能量等于所述煤体破坏损失能量和所述破碎煤体临界抛出能量之和；

通过相应的能量密度来表达相应的所述煤体破坏损失能量和所述破碎煤体临界抛出能量；

所述煤体损伤本构模型的公式为：

Y为损伤耗能率；

σ为应力；

E为弹性模量；

D为损伤变量，表示单位体积微裂纹数；

DE为在位置E完全断裂破坏时的损伤值；

DD为在位置D完全断裂破坏时的损伤值；

UD为煤岩体损伤消耗能量；

所述煤体损伤本构模型中的损伤变量通过矿井中煤体的应力-应变试验关系确定；

所述破碎煤体临界抛出能量采用以下公式：

UF：破碎煤体临界抛出能量；

ρ：破碎煤体平均密度；

ν：破碎煤体向自由空间抛出的平均初速度。

4.一种确定冲击地压临界能量的设备，其特征在于，包括：

一个或多个处理器；

存储器，用于存储一个或多个程序；

当所述一个或多个程序被所述一个或多个处理器执行时，使得所述一个或多个处理器执行如权利要求1-2任一所述的确定冲击地压临界能量的方法。

5.一种存储有计算机程序的计算机可读存储介质，其特征在于，该程序被处理器执行时，实现如权利要求1-2任一所述的确定冲击地压临界能量的方法。

Images