本发明涉及煤矿开采，特别涉及一种平切煤体破裂范围确定方法。
背景技术：
1、随着我国煤矿开采逐渐向深部延伸，地质条件趋于复杂。深部煤层高地应力、高地温和高瓦斯压力的特性导致其普遍具有微孔隙、低渗透和强吸附的特征，瓦斯抽采极为困难，严重制约了矿井灾害防治和瓦斯资源利用。通过外界扰动作用，人为诱导煤体裂隙萌生、扩展及贯通，创造瓦斯流动通道是增加煤层透气性、提高瓦斯抽采效果的主要思路。水力割缝能够使煤层内应力场重新分布，煤层裂隙的数量、长度、张开度及连通面积得到增加，透气性得到明显改善。高压射流切缝造成的扰动与煤体损伤破裂密切相关。然而，射流扰动和缝槽流变影响因素众多，其中主控因素难以确定，导致不同影响因素作用下平切煤体破裂范围尚不清晰。
2、因此，亟需开发一种平切煤体破裂范围确定方法。
技术实现思路
1、本发明的目的是提供一种平切煤体破裂范围确定方法，以解决现有技术中存在的问题。
2、为实现本发明目的而采用的技术方案是这样的，一种平切煤体破裂范围确定方法，包括以下步骤，其中，步骤1)和步骤2)同时进行或不分先后顺序：
3、1)构建平切煤体破坏影响因素集合r。
4、2)提取平切扰动煤体破坏表征参数。
5、3)使用基于聚类的特征选择方法确定平切扰动下缝槽破坏主控因素。
6、4)获取现场矿井煤层基本物理力学参数，构建相似模拟物理模型。
7、5)开展高压射流平切煤体相似实验，在相似模型中布置应变检测仪，实时得到煤体应力。根据切缝过程中射流压力的变化情况，构建射流扰动阶段煤体应力分布本构模型。
8、6)根据切缝后缝槽间距和缝槽宽度的变化情况，确定含缝槽煤体变形长期演化规律，构建缝槽流变阶段煤体应力分布模型。
9、7)根据射流扰动阶段中煤体参数，包括弹性模量、粘聚力、内摩擦角，构建煤体损伤本构模型，得到切缝过程中的煤体破坏临界应力σ1，描述煤体在水射流冲击下的破坏行为。
10、8)根据缝槽流变阶段中煤体参数，包括累积应变及流变速率，构建损伤煤体流变模型，得到切缝后的煤体破坏临界应力σ2，描述煤体在流变效应下的破坏行为。
11、9)将步骤7)和步骤8)中得到的煤体应力，根据缝槽破坏临界条件，与煤体破坏临界应力进行对比，判断煤体是否发生破坏。其中，若射流扰动阶段煤体应力大于σ1，则煤体在切缝过程中发生了破坏。若缝槽流变阶段煤体应力大于σ2，则煤体在切缝后发生了破坏。
12、进一步，步骤1)中，通过对高压射流切割煤体的文献调研和现场资料分析，构建平切煤体破坏影响因素集合。所述影响因素集合包括射流扰动影响因素和缝槽流变影响因素。所述射流扰动阶段的影响因素包括冲击压力、冲击时间和冲击流量。所述缝槽流变阶段的影响因素包括缝槽宽度、缝槽间距、煤体强度和煤体含水率。
13、进一步，使用动力学软件对高压射流切割煤体进行数值模拟，改变射流扰动影响因素和缝槽流变影响因素，基于数值模拟结果变化提取平切扰动煤体破坏表征参数。所述表征参数包括射流扰动阶段破坏表征参数和缝槽流变阶段破坏表征参数。所述射流扰动阶段破坏表征参数包括峰值应力和塑性能占比。所述缝槽流变阶段破坏表征参数包括累计应变和流变速率。
14、进一步，步骤2)之后，还具有开展高压射流切割煤样力学实验，随机抽取并改变部分影响因素，根据特征参数变化验证数值模拟结果的步骤。
15、进一步，步骤3)具体包括以下子步骤：
16、3.1)使用聚类算法将集合r中的因素进行聚类，获得聚类结果ri＝{r1，r2，…，rk}。其中，聚类条件为每个集合ri的改变均能够使煤体破坏表征参数发生对应的改变。
17、3.2)计算每个集合ri到其他任意一个集合ri′的区分度ri。
18、计算公式如式(1)所示。
19、
20、式中，|ri|为集合ri的大小，为集合r1和r2在特征参数上的差异。
21、3.3)根据所有区分度计算每个集合ri的总区分度平均值meani，及针对每个表征参数的平均区分度最大值maxi和平均区分度最小值mini。
22、3.4)计算区分度参数fi。计算公式如式(2)所示。
23、
24、3.5)对集合ri按照区分度参数fi降序排列，得到在的折线图中找到急剧变化的点或拐点i0，即为从影响因素集合r中选择的破坏主控因素。
25、进一步，步骤5)中，在引入应变率和损伤因子基础上，将材料的等效应力通过式(3)进行量化表征：
26、
27、当煤体的损伤因子d＝1时煤体完全失效，等效应力如式(4)所示：
28、
29、式中，σi为煤体的实际破裂强度，hel为hugoniot弹性极限。p*为规范化静水压力，p*＝pa/phel，pa为实际压力。phel为煤体处于hugoniot弹性极限时的压力。为煤体实际等效应变率，为10s-1参考应变率。kn为第n个主控因素影响下的损伤变量，εn为第n个主控因素影响下的应变，ε0为无损状态下的应变。c为煤体常数。b和m为煤体残余强度常数。
30、进一步，步骤6)中，基于改进后的西原模型构建缝槽流变阶段煤体应力分布模型。缝槽流变阶段煤体应力分布模型由一个胡克体、一个开尔文体和一个理想粘塑性体串联而成。考虑缝槽周围煤体损伤改进的西原模型的一维流变方程表达式如式(5)所示：
31、当σ<σs时，流变方程表达为：
32、
33、当σf＞σ≥σs时，流变方程表达如式(6)所示：
34、
35、当σ≥σf时，流变方程表达如式(7)所示：
36、
37、其中，h为考虑主控因素影响下的损伤变量如式(8)所示：
38、
39、式中，hn为考虑第n个主控因素影响下的损伤变量，en为考虑第n个主控因素影响下的弹性模量。e0为无损状态下的弹性模量。
40、本发明的技术效果是毋庸置疑的：
41、a.采用调研分析、数值模拟、物理实验相结合的方法，从射流扰动阶段及缝槽流变阶段分别开展研究，并引入基于聚类的特征选择方法，从而精准确定影响煤体破裂的主控因素；
42、b.通过理论研究，对煤体应力分布本构模型进行优化，实现了平切煤体破裂条件的准确判断；
43、c.基于主控因素与煤体应力的函数关系，通过数值模拟实现了不同主控因素下平切煤体在空间、时间双尺度下的破裂范围确定；
44、d.发明对维持缝槽稳定、提升卸压效果有重要意义，避免了瓦斯运移通道堵塞，从而提高了瓦斯抽采效果。
1.一种平切煤体破裂范围确定方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的一种平切煤体破裂范围确定方法，其特征在于：步骤1)中，通过对高压射流切割煤体的文献调研和现场资料分析，构建平切煤体破裂影响因素集合；所述影响因素集合包括射流扰动阶段影响因素和缝槽流变阶段影响因素；所述射流扰动阶段的影响因素包括冲击压力、冲击时间和冲击角度；所述缝槽流变阶段的影响因素包括缝槽宽度、缝槽间距、煤体强度和煤体含水率。
3.根据权利要求1所述的一种平切煤体破裂范围确定方法，其特征在于：步骤4)中，选定一个水平的平切煤体破裂影响因素，通过监测相似模拟实验中的煤体应力，得到此水平下煤体破裂表征参数，同时为数值模拟的结果提供对比，以验证数值模拟的可靠性。所述表征参数包括射流扰动阶段破裂表征参数和缝槽流变阶段破裂表征参数；所述射流扰动阶段破裂表征参数包括峰值应力和塑性能占比；所述缝槽流变阶段破裂表征参数包括累计应变和流变速率。
4.根据权利要求1所述的一种平切煤体破裂范围确定方法，其特征在于：步骤5)中，使用仿真软件对高压射流切割煤体和割缝煤体长期流变过程分别进行数值模拟，分别改变射流扰动阶段影响因素和缝槽流变阶段影响因素的不同水平，获得各影响因素在不同水平下的破裂表征参数。
5.根据权利要求1所述的一种平切煤体破裂范围确定方法，其特征在于，步骤6)具体包括以下子步骤：
6.根据权利要求1所述的一种平切煤体破裂范围确定方法，其特征在于，步骤7)具体包括以下子步骤：
7.根据权利要求1所述的一种平切煤体破裂范围确定方法，其特征在于：步骤8)中，在引入应变率和损伤因子基础上，将材料的等效应力通过式(6)进行量化表征：
8.根据权利要求1所述的一种平切煤体破裂范围确定方法，其特征在于：步骤9)中，基于改进后的西原模型构建缝槽流变阶段煤体应力分布模型；缝槽流变阶段煤体应力分布模型由一个胡克体、一个开尔文体和一个理想粘塑性体及一个可以反映加速蠕变阶段的黏塑性模型串联而成；考虑缝槽周围煤体损伤的改进西原模型的流变方程表达式如式(8)所示：
9.根据权利要求1所述的一种平切煤体破裂范围确定方法，其特征在于：步骤12)中，通过下式判断煤体是否发生破裂：
