Research on the potential of gas resource development in abandoned mines using matter element extension theory
-
摘要:
“双碳”目标下,为了精准研判废弃矿井采动稳定区瓦斯资源开发潜力和找出其影响开发潜力的关键因素,构建了“精选−内联−互补−优化”四位一体为理念的废弃矿井瓦斯资源开发潜力评价方法,最后与工程试验对比,验证评估模型的可行性。“精选”指:首先初步选取大量的评价指标,然后通过德尔菲法和敏感性分析筛选初始指标,去除不符合要求的指标因素;“内联”指:对筛选出的指标,采用DEMATEL-ISM进行内部关联性分析,判定关键指标;“互补”指:首先采用FAHP和熵值法分别对指标进行定性和定量的赋权,然后采用博弈论对指标进行综合权重确定;“优化”指:结合综合权重,采用物元可拓理论计算其关联度,判定废弃矿井瓦斯资源开发潜力等级。结果表明:煤层埋深、煤矿面积是判定废弃矿井瓦斯资源开发潜力的关键指标因素与DEMATEL-ISM模型评价一致;根据物元可拓模型计算的关联度为K3=
0.0575 ,可以判定其开发潜力为“良好”,根据工程试验统计评估可知,13个可采煤层的瓦斯资源量为61.81×108 m3,潘一矿瓦斯资源规模较大,属于中型煤层气田,其开发潜力较大与物元可拓模型评价结论一致,符合实际情况,该模型为研究废弃矿井瓦斯资源开发潜力等级提供方向。-
关键词:
- 废弃矿井 /
- DEMATEL-ISM /
- 模糊层次分析法(FAHP) /
- 熵值法 /
- 博弈论 /
- 物元可拓模型
Abstract:Based on the background of national carbon peak and carbon neutrality, in order to accurately judge the development potential of gas resources in the mining stability area of abandoned mines and find out the key factors affecting the development potential, the author constructs the “selection-inline-complementary-optimization”. “Selecting” refers to: firstly, a large number of evaluation indexes are preliminarily selected, and then the initial indexes are screened by Delphi method and sensitivity analysis to remove the index factors that do not meet the requirements. “Inline” refers to: for the selected indicators, DEMATEL-ISM is used for internal correlation analysis to determine the key indicators. “Complementarity” refers to: firstly, FAHP and entropy method are used to determine the qualitative and quantitative weights of the indicators, and then the game theory is used to determine the comprehensive weights of the indicators. “Optimization” refers to a four-in-one evaluation method for the development potential of abandoned mine gas resources based on the concept of comprehensive weight, using matter-element extension theory to calculate its correlation degree, and determining the development potential level of abandoned mine gas resources. Finally, the feasibility of the evaluation model is verified by comparison with engineering tests. The results show that the buried depth of coal seam and the area of coal mine are the key index factors to determine the development potential of gas resources in abandoned mines, which are consistent with the evaluation of DEMATEL-ISM model. According to the correlation degree calculated by the matter-element extension model, K3 =
0.0575 , it can be judged that its development potential is “good”. According to the statistical evaluation of engineering tests, the gas resources of 13 coal seams are 61.81 × 108 m3. The gas resources of Panyi Mine are large and belong to medium-sized coalbed methane fields. Its development potential is consistent with the evaluation conclusion of the matter-element extension model, which is in line with the actual situation. The model provides a direction for studying the development potential level of gas resources in abandoned mines. -
0. 引 言
煤炭是我国的基础能源,在较长时期内仍将是我国的主体能源并发挥着国家能源安全稳定供应压舱石的作用[1]。我国浅部煤炭资源逐渐枯竭,深部资源开采成为必然趋势。随着开采深度和开采强度的增大,地应力、瓦斯压力和地温不断增加,采动影响也越来越大,深部煤炭采掘过程将伴随着更加复杂的煤与瓦斯突出灾害,严重威胁煤矿安全生产[2-5]。
在煤层群条件下,应优先采用保护层开采措施,该项技术在多煤层有效增透卸压、实现大面积消突上具有重要作用。开采保护层是被公认为最有效、最经济的区域防突措施,《防治煤与瓦斯突出细则》规定[6]:“具备开采保护层条件的突出危险区,必须开采保护层”。对于不具备保护层开采条件的单一突出煤层目前主要采用预抽煤层瓦斯区域防突措施,《防治煤与瓦斯突出细则》针对区域防突措施提出了多种区域瓦斯抽采消突方式。其中,最为常用的是穿层钻孔和顺层钻孔瓦斯抽采技术。但对强突出单一煤层主要适合采用穿层钻孔预抽,并采取人为增透措施强化瓦斯抽采消突。
煤层卸压与煤体裂缝改造是低渗煤层增透两个主要技术途径,目前单一低渗煤层增透技术主要可分为水力化技术、爆破类技术、钻扩类技术以及其他新技术。其中,水力化技术主要包括:水力压裂、水射流(磨料水射流)割缝、水力冲孔等;爆破类技术主要包括:深孔爆破、CO2相变致裂、高压空气爆破致裂、可控冲击波等;钻扩类技术主要包括:大直径钻孔、超短半径侧钻等;其他新技术包括注液氮冷加载、微波辐射、超声波、变频气动致裂、化学增透等。由此可见,对强突出煤层,卸压增透抽采消突技术更多层面考虑低渗煤层的增透抽采消突效应,对卸压消突因素考虑不足,缺少对相关增透措施的大区域连片均匀卸压消突理论及构建方式等的研究,导致穿层钻孔卸压增透预抽煤层瓦斯技术,在局部煤岩层卸压、增透强化抽采等方面虽然取得较好的效果,但仍存在消突范围小、时间长、消突区域不均衡等关键技术问题。在无保护层开采条件下,通过对煤层强制卸压,致使煤层大范围均匀卸压,形成类似于在煤层内部构建保护层的卸压开采效果,实现高效消突等,需要开展深入的理论研究、工程实践及配套技术装备研发。
笔者以焦作矿区为研究对象,以地质勘探为基础,利用底抽巷精准施工均匀、高密度穿层钻孔,通过穿层钻孔均匀合理掏煤,在单一煤层内形成“内保护层”,均匀卸压、増透构建内保护层,使其达到保护层开采的卸压消突效果。即以煤层卸压消突为目标,以内保护层消突技术参数和高效构建内保护层技术装备为依托,以系统保障体系和消突效果考察为保障,实现构建的内保护层全域穿层钻孔“布孔均匀密集、卸煤科学高效、邻孔串气相透、卸压均匀充分”,通过均匀卸压与瓦斯抽采实现区域高效消突。
1. 内保护层技术科学内涵
保护层开采是目前最为有效的区域防突手段,因此《煤矿安全规程》和《防治煤与瓦斯突出细则》中明确规定,具备开采保护层条件的突出危险区必须开采保护层作为区域防突措施。但对于不具备保护层开采条件的单一突出厚煤层,《细则》也明确指出:“利用上分层或者上区段开采后形成的卸压作用保护下分层或者下区段···”[6]。这种卸压保护作用与保护层开采的消突原理殊途同归。
虽然在分层开采过程中,顶分层开采后对底分层起到较好的卸压保护作用,但突出煤层在顶分层开采过程中仍然存在突出危险性。因此笔者借鉴单一突出煤层分层开采的卸压消突作用,提出了内保护层技术。该技术是以卸除突出发生的应力条件与增透强化抽采相结合的高效消突措施,通过合理设计穿层钻孔间距、优化不同煤层的卸压措施和合理的掏出煤量,强化卸压抽采及消突效果考察,实现大区域连片均匀卸压增透,可形成类似于在煤层内构建类似保护层的消突途径,实现全区域充分卸压、快速消突。
焦作矿区主采煤层为山西组二1,在其上部或下部赋存着一层厚度为0.2~1.2 m的构造软分层,这一构造软分层突出危险性极大。为此,结合焦作矿区煤层赋存实际情况,并在单一煤层分层开采和保护层开采的理论基础上,提出了针对单一突出煤层构建内保护层的技术思想,根据软分层厚度或硬煤层赋存条件,构建内保护层消突技术的原理分别如图1、图2所示。
第Ⅰ阶段,在煤层底板岩巷中,利用水力冲孔、气动柔性刀具、超高压水力割缝、磨料水射流等高效破煤卸压技术手段,提前将煤层顶部或底部的构造软煤掏出,也可在硬煤层中通过水力割缝等技术掏煤,减少强突出煤体的总量,以密集钻孔、强效卸煤的方式在单一煤层内部主动构建出保护层。
第Ⅱ阶段,在构建内保护层的卸压作用下,软煤分层因应变软化而产生塑性流变,硬煤因应力降低而发生膨胀变形,煤层整体孔裂隙得到发育,煤层瓦斯压力梯度降低,煤层弹性潜能和瓦斯膨胀能开始释放。
第Ⅲ阶段,煤层实现区域连片均匀卸压,煤层孔裂隙充分发育透气性增加,形成“邻孔串气”效应,利用顶板“源汇抽采”技术,煤层瓦斯压力及含量得到快速下降,煤层弹性潜能及瓦斯膨胀能降至临界突出能量以下,煤层整体实现消突。
构建内保护层技术以地质勘探为基础,以煤层卸压消突为目标,以内保护层消突技术参数和高效构建技术装备为依托,以消突效果考察和系统保障体系为保障,实现穿层钻孔的“布孔均匀密集、卸煤科学高效、邻孔串气相透、卸压均匀充分”,统筹了煤层区域高效均匀卸压和瓦斯抽采,适用于单一突出煤层的主动构建保护层的消突理念。
2. 构建内保护层消突力学判据及评判指标
2.1 煤与瓦斯突出的能量条件
煤与瓦斯突出是一种复杂的动力现象,涉及岩石力学、渗流力学、地球物理学等多领域的复杂问题,突出机理至今仍是国际性、世纪性难题。但究其根源,突出过程实际上是能量的转移、聚积和释放过程,这个过程必须满足热力学第一定律。因此从能量守恒的角度,可总体上描述突出发生的条件。对于某一范围内的煤岩体,煤与瓦斯突出的能量条件可以表示如下[7]:
$$ {{W}}_{1}+{{W}}_{2}+{{W}}_{3}+{{W}}_{4}={{A}}_{1}+{{A}}_{2}+{{A}}_{3} $$ (1) 式中:W1为突出范围内煤岩体的弹性势能;W2为突出范围内煤岩体的瓦斯内能;W3为突出范围内煤岩体的重力势能;W4为周围煤岩体对突出范围内煤岩所做的功;A1为煤岩体的破碎功;A2为煤岩体的抛出功;A3为其他能量耗散,如煤体发热、声发射等。
式(1)所表达的是突出发生必然遵守的能量条件。在突出发动中,由于突出煤体的弹性模量远小于煤层顶、底板岩石的弹性模量,导致煤层顶板、底板弹性潜能对突出做功远小于煤体弹性潜能,同时煤壁失稳后顶底板可能来不及快速下沉对煤体做功,煤体就已经被抛出了。其次,突出范围内煤岩体的重力势能 W3只在一些特殊的开采环境中,才会对突出的发生起到一定作用,因此为了便于计算和分析,可忽略顶板、底板对煤层做功及重力势能。同时也不考虑其他的如声发射、热量放散等的能量损耗[8]。近似的认为瓦斯内能和煤体弹性潜能是突出的主要能量来源。这样,发生煤与瓦斯突出的能量关系可简化为式(2):
$$ {{W}}_{1}+{{W}}_{2}={{A}}_{1}+{{A}}_{2}$$ (2) 2.2 突出过程中煤层释放的弹性势能
为了便于分析,以煤巷掘进为研究对象,将掘进工作面看作半无限体内的球形孔洞,同时兼顾地应力分布特征,建立煤与瓦斯突出应力模型,如图3所示。其中R0为原生破裂带半径;Rp为塑性区半径;L为掘进进尺。
模型将掘进工作面前方区域分为2部分,分别为靠近掘进头的塑性区和塑性区之外的弹性区。由图3可知,掘进工作面向前推进L后,部分弹性区转化成为塑性区。根据对煤样破坏过程能量演化规律可知,煤岩破坏后煤体弹性能释放且瓦斯膨胀能突变。因此,掘进工作面前方煤体由弹性区转换为塑性区后将引起能量变化,该部分即为能量释放区,如图中Vs[9-11]。
通过岩石力学的计算,在工作面前方塑性区内,由平衡微分方程和摩尔−库仑强度准则组成的方程组,以及边界条件:当r=R0时,σrp=0,可计算得出塑性区范围Rp为:
$$ {R}_{{\mathrm{p}}}={R}_{0}{\left[\frac{3({\sigma }_{0}+c\cot \;\varphi )(1-\sin \;\varphi )}{(3+\sin \;\varphi )c\cot \;\varphi }\right]}^{\tfrac{1-\sin \;\varphi }{4\sin \;\varphi }} $$ (3) 式中:σ0为地应力,MPa;φ为内摩擦角,(°);c为黏聚力,MPa。
在弹性区内,由平衡微分方程、几何方程和本构方程,以及边界条件:r =+∞时,σre=σ0,当r =Rp 时,σre=σrp。能量释放区内煤体释放的能量可表示为
$$ {W}_{1}={\int }_{{V}_{{\mathrm{s}}}}^{}{U}^{e}{\rm{d}}V={\iiint }_{{V}_{{\mathrm{s}}}}^{}{U}^{{\mathrm{e}}}\sin \;\psi {r}^{2}{\rm{d}}\psi {\rm{d}}\theta {\rm{d}}r $$ (4) 式中:Ue为单位体积煤体可释放的弹性应变能;μ为泊松比;Vs为煤体弹性释放区体积。
式(4)中各积分变量的取值范围,可通过图4中的各参数的关系求得,设掘进前,破裂带球心为O,进尺L后,新的破裂带球心为O1,OO1=L,假设原破裂带、新破裂带半径均为RP。
过O画一条直线,与x轴夹角为θA,与原塑性区、新塑性区分别相交于A'(RP,ψA),A(rA,ψA),能量释放区内煤体释放的能量可表示为
$$ \begin{array}{c} {W_1} = \dfrac{\pi }{E}\displaystyle \int _0^{{\mathrm{arccos}} \tfrac{L}{{2{R_{\mathrm{p}}}}}}{\mathrm{d}}\psi \displaystyle \int _{{R_{\mathrm{p}}}}^{L\cos \;\psi + \sqrt {R_{\mathrm{p}}^2 - {L^2} - \sin{^2}\psi } }\left[ \left( {3 - 6\mu } \right){\sigma _0} + \right.\\ \left.\dfrac{{3R_{\mathrm{p}}^2{k^2}(1 + \mu )}}{{2{r^6}}} \right]{r^2}{\rm{d}}r \end{array}$$ (5) 式中:$k={-\sigma }_{0}+c\mathrm{c}\mathrm{o}\mathrm{t}\;\varphi \left[{\left(\dfrac{{R}_{{\mathrm{p}}}}{{R}_{0}}\right)}^{\tfrac{4\sin\; \varphi }{1-\sin \;\varphi }}-1\right]$。
对式(5)进行数值积分,即可求得煤层弹性释放能。
2.3 突出过程中煤层瓦斯膨胀内能
1)突出过程中煤层瓦斯膨胀内能的计算方法。
瓦斯膨胀能是突出过程的重要动力来源[12]。目前瓦斯膨胀能确定主要基于能量计算公式进行。但由于突出实际工况的复杂性,能量公式没有统一的理论假设,导致现阶段膨胀能计算公式繁多,且进行膨胀能计算所必需的参与突出过程的瓦斯量确定也没有统一方法[13]。为此,以目前主流的膨胀能计算公式进行分析和计算。根据热力学定理及突出的实际工况,将突出过程看作为等温多变过程,膨胀内能可按下式进行计算:
$$ {W}_{2}=V{P}_{0}\mathrm{l}\mathrm{n}\left(\frac{{P}_{0}}{{P}_{{\mathrm{a}}}}\right) $$ (6) 式中:Pa为大气压,MPa;P0为突出煤层瓦斯压力,MPa;V为参与突出的瓦斯体积,m3。
2)参与突出过程中瓦斯量的计算。
突出发生时,能够释放出来参与做功的瓦斯内能包括两部分,一部分是游离态的瓦斯体积;另一部分是吸附态的瓦斯。
游离瓦斯在能量释放区范围内参与突出的游离瓦斯含量可用下式计算。
$$ V_{{\mathrm{f}}}=nV_{{\mathrm{s}}}$$ (7) 式中:Vf为参与突出过程中游离瓦斯体积,m3;n为煤体孔隙率;Vs为煤体能量释放区体积。
刘明举等[14]认为,参与突出的瓦斯量大于游离瓦斯含量,小于瓦斯总含量。李成武等[15]认为,参与突出的瓦斯量与突出煤体内游离瓦斯含量接近。胡千庭等[7]采用一定比例的游离瓦斯膨胀能来表征吸附瓦斯的能量大小。为了准确性和可行性,借鉴胡千庭的思路,将参与突出吸附瓦斯含量采用一定比例的游离瓦斯含量来表征。能量释放区内参与突出的瓦斯膨胀内能则可通过下式表示,同时通过试验测定,η=1.5~2。
$$ W_2=\eta n V_{\mathrm{s}} P_0 \ln \left(\frac{P_0}{P_{\mathrm{a}}}\right) $$ (8) 3)突出过程中煤体破碎功。
基于破碎功的基本计算原理,根据前人对全国各主要突出矿井21个煤样的冲击破碎试验得到煤体破碎功的表达式[16]为
$$ {A}_{1}=46.914{f}^{1.437}{Y}_{{\mathrm{P1}}}^{1.679}{V}_{{\mathrm{S}}}\rho $$ (9) 式中:YP1为破碎成0.2 mm以下粒度煤样质量占总煤样质量的百分比,取20;ρ为煤体密度。
4)突出过程中煤体抛出功。
煤体的抛出功就是碎煤抛出时的动能[17]:
$$ {A}_{2}=\frac{{V}_{{\mathrm{S}}}\rho {v}^{2}}{2} $$ (10) 式中:A2为煤的抛出功;v为碎煤抛出的速度,m/s,根据前苏联“红色国际工会”矿杰列卓夫卡煤层的突出实测资料,抛出碎煤的速度在17.6~55.5 m/s。
2.4 内保护层消突评判分析
在煤与瓦斯突出的历程中,地应力、瓦斯压力、煤的物理力学性质,三者共同决定了煤与瓦斯突出的危险程度。煤与瓦斯突出与三者之间的关系可以用式(11)表达[18]。
$$ {{F}}={{f}}\left(\sigma 、{{P}}、{{G}}\right) $$ (11) 式中:σ 为煤层的应力状态;P为煤层瓦斯情况;G为煤层的力学性质。
将式(5)、(8)、(9)和(10)代入式(11),得到煤与瓦斯突出的能量条件为
$$ \begin{array}{c} \dfrac{\pi}{E} \displaystyle\int_0^{\arccos \tfrac{L}{2 R_{\mathrm{p}}}} {\mathrm{d}} \psi \displaystyle\int_{R_{\mathrm{p}}}^{L \cos \psi+\sqrt{R_{\mathrm{p}}^2-L^2-\sin ^2 \psi}}\left[(3-6 \mu) \sigma_0+\right. \\ \left.\dfrac{3 R_{\mathrm{p}}^2 k^2(1+\mu)}{2 r^6}\right] r^2 {\mathrm{d}} r+\eta n V _{\mathrm{s}} P_0 \ln \left(\dfrac{P_0}{P_{\mathrm{a}}}\right)= \\ 46.914 f^{1.437} Y_{{\mathrm{P 1}}}^{1.679} V_{\mathrm{S}} \rho+\dfrac{V_{\mathrm{S}} \rho v^2}{2} \end{array} $$ (12) 根据式(12)编制计算程序,可计算出不同地应力条件下,发生煤与瓦斯突出的临界压力。
结合焦作矿区煤层基本情况:主采煤层山西组二1煤,煤层中存在一层软分层,厚度为0.2~1.2 m,无临近开采煤层,煤层厚度6~7 m。经现场实测,软分层煤体f值为0.2~0.3,而硬煤f值为0.5~2.0。根据式(12),可分别计算出软煤和硬煤在不同地应力条件下,发生煤与瓦斯突出的临界瓦斯压力。在此,选取软煤f值为0.2,硬煤f值为0.5的极限情况。软煤和硬煤相关参数取值及计算结果见表1和表2。
表 1 软、硬煤突出临界瓦斯压力计算参数值Table 1. Thecritical gas pressure values for outburst soft and hard coal类型 弹性模量/MPa 黏聚力
/MPa泊松比 内摩擦角
/(°)孔隙率
/%巷道半径
/m掘进进尺
/mη 软煤 800 0.4 0.28 30 6 3 1 2 硬煤 1200 3 0.3 35 4 3 1 2 表 2 软、硬煤不同应力条件下的突出临界瓦斯压力值Table 2. Outburst critical gas pressure values for soft and hard coals under different stress conditionsMPa f σ0 0 2 4 6 8 10 12 14 16 18 20 22 0.2 1.48 1.26 1.10 0.95 0.84 0.73 0.66 0.6 0.53 0.47 0.42 0.39 0.5 2.81 2.25 1.87 1.59 1.40 1.25 1.13 1.03 0.94 0.87 0.81 0.77 目前,焦作矿区下属各矿井开采深度已达到500 m,实测地应力为14 MPa。根据表1,软煤与硬煤在14 MPa地应力环境下,理论突出临界瓦斯压力分别为0.6 MPa和1.03 MPa。经现场实测瓦斯压力为1.62 MPa,大于理论突出临界瓦斯压力,因此煤层具有突出危险性。
图5为焦作矿区软、硬煤在不同地应力与理论突出临界瓦斯压力关系的图。由图5可知,无论是构造软煤还是硬煤,突出临界瓦斯压力均随地应力的降低而升高,降低煤层地应力可有效起到防治煤与瓦斯突出的作用。但从曲线的陡峭程度来看,煤体f值越大,地应力对临界突出瓦斯压力越敏感,即小幅度地降低煤层地应力,就可大幅提高突出临界煤层瓦斯压力,是由于硬煤弹性模量较大,煤体释放的弹性能是突出发动过程中的主要能量来源,地应力的减少可有效降低煤层自身突出总能量储备。而煤体f值越小,地应力对突出临界瓦斯压力的敏感度越低,即构造软煤降低应力的消突效果无法像硬煤那样显著,因为构造软煤的强度低,煤层瓦斯的膨胀能是突出发动过程中的主要能量来源。综上所述,硬煤的卸压消突作用要优于煤层瓦斯抽采,而软煤瓦斯抽采的消突作用要优于煤层卸压。
开采保护层是目前较为有效的区域防突技术措施,而该技术的本质是使被保护层得到充分卸压,自身的弹性储能得到有效释放,从而引起煤层发生膨胀变形。同时,《细则》中规定:“若经实际考察被保护层的最大膨胀变形量大于3‰,则检验和考察结果可适用于具有同一保护层和被保护层关系的其他区域。”因此,煤层充分卸压是消除突出的首要任务。
构建“内保护层”技术是针对单一煤层特殊的保护层开采消突技术,煤层的充分卸压同样是该技术所要达成的首要目标。同时,考虑到瓦斯膨胀内能对煤与瓦斯突出过程中的作用,尤其是软煤。结合目前焦作矿区地应力14 MPa、瓦斯压力1.62 MPa软、硬复合煤层的开采环境分析可知:
1)对于硬煤而言,在目前的开采煤层底板标高500 m,实测地应力14 MPa,可计算得出当前开采深度的理论突出临界瓦斯压力为 1.03 MPa,实测煤层瓦斯压力为1.62 MPa,大于理论突出临界瓦斯压力,因此,存在突出危险性。为实现硬煤层的有效消突,在不进行煤层瓦斯预抽的情况,煤层应力需从14 MPa降至5.9 MPa 后,卸压程度达到58%以上,理论上可以实现消除煤层突出危险性。
2)对于软煤而言,即使卸压程度达到了 100%,软煤所处应力为0,其理论突出临界瓦斯压力1.48 MPa,实际煤层瓦斯压力仍高于理论临界瓦斯压力,因此软煤无法像硬煤那样,仅靠煤层卸压就可到达消突的目的,须在对软煤进行卸压的同时,进行瓦斯抽采才能达到消突的目的。
3)针对焦作矿区特殊的复合煤层,以煤层总体预抽率为30%为前提条件,对于软分层其卸压程度至少应达到40%以上,而硬煤理论上即使不卸压也可消除突出危险,但考虑到硬煤卸压对煤层弹性潜能的释放作用,以及对煤层增透强化抽采的效果,硬煤卸压程度也应达到10%以上。
3. 内保护层构建技术参数研究
3.1 内保护层构建正交试验方案
为了指导不同地质条件下的内保护层参数构建,设计了正交试验,研究不同地质条件及不同施工参数下的内保护层构建消突效果,通过分析不同条件下卸压控制范围内的最薄弱点应力、应变变化规律,结合前文突出判定力学依据,分析各施工参数下是否消突,选择合适的目标函数,构建内包层施工效果预测模型,并结合预测模型分析得出各影响因素的权重,最终为后续现场施工提供指导。
方案中,构建了4种不同地质条件下的正交试验分析表,分别是构建上部内保护层正交模拟方案,构建下部内保护层正交模拟方案,全软煤时正交模拟方案以及全硬煤时正交模拟方案。详细参数见表3—表6。
表 3 构建上部内保护层正交试验分析(工况1)Table 3. Orthogonal test scheme for construction of upper inner protective layer (working condition 1)编号 钻孔倾角/
(°)出煤量/
(t·m−1)孔间距/
m原始应力/
MPa软煤厚度/
m1 45 0.6 5 12 0.3 2 45 1 6 10 0.9 3 30 0.6 6 16 0.6 4 60 0.8 6 14 0.3 5 30 0.4 4 10 0.3 6 30 1.0 5 14 1.2 7 60 0.4 5 16 0.9 8 30 0.8 7 12 0.9 9 75 0.8 5 10 0.6 10 45 0.8 4 16 1.2 11 60 1.0 4 12 0.6 12 45 0.4 7 14 0.6 13 75 0.4 6 12 1.2 14 60 0.6 7 10 1.2 15 75 1.0 7 16 0.3 16 75 0.6 4 14 0.9 表 4 构建下部内保护层正交试验分析(工况2)Table 4. Orthogonal test scheme for construction of lower inner protective layer (working condition 2)编号 钻孔倾角/
(°)出煤量/
(t·m−1)孔间距/
m原始应力/
MPa软煤厚度/
m1 45 0.6 5 12 0.3 2 45 1.0 6 10 0.9 3 30 0.6 6 16 0.6 4 60 0.8 6 14 0.3 5 30 0.4 4 10 0.3 6 30 1.0 5 14 1.2 7 60 0.4 5 16 0.9 8 30 0.8 7 12 0.9 9 75 0.8 5 10 0.6 10 45 0.8 4 16 1.2 11 60 1.0 4 12 0.6 12 45 0.4 7 14 0.6 13 75 0.4 6 12 1.2 14 60 0.6 7 10 1.2 15 75 1.0 7 16 0.3 16 75 0.6 4 14 0.9 表 5 全软煤正交试验正交试验方案(工况3)Table 5. Orthogonal test scheme for all-soft coal (working condition 3)编号 钻孔倾角
/(°)出煤量
/(t·m−1)孔间距
/m原始应力
/MPa1 60 0.6 7 14 2 60 1.0 5 10 3 75 0.6 6 10 4 45 0.4 5 14 5 30 0.6 5 12 6 75 1.0 4 14 7 30 1.0 7 16 8 75 0.8 6 16 9 30 0.8 6 14 10 60 0.8 4 12 11 45 0.8 7 10 12 45 1.0 6 12 13 60 0.4 6 16 14 75 0.4 7 12 15 30 0.4 4 10 16 45 0.6 4 16 表 6 全硬煤正交试验分析(工况4)Table 6. Orthogonal test scheme for all-hard coal (working condition 4)编号 钻孔倾角/(°) 原始应力/MPa 孔间距/m 割缝长度/MPa 割缝数量/个 割缝高度/cm 1 45 3.5 10 0.3 2 2 2 60 3.5 12 0.9 4 2 3 30 2.5 10 0.9 8 4 4 30 2.0 16 0.9 2 2 5 30 2.5 14 0.3 6 2 6 30 3.5 14 0.6 2 8 7 60 2.5 16 1.2 2 8 8 75 2.5 10 0.6 4 2 9 45 3.0 10 0.9 2 6 10 75 3.0 10 1.2 2 2 11 75 2.0 14 0.9 2 6 12 60 3.0 14 0.3 2 6 13 45 2.5 12 0.3 2 6 14 60 2.0 10 0.6 6 6 15 30 3.5 10 1.2 8 6 16 30 2.0 10 0.3 2 2 17 60 2.0 10 0.3 2 4 18 30 2.0 12 1.2 6 2 19 30 3.0 16 0.3 4 6 20 30 3.0 12 0.6 2 4 21 30 2.0 10 0.3 4 8 22 75 2.0 12 0.3 8 8 23 45 2.0 16 0.6 8 2 24 75 3.5 16 0.3 6 4 25 45 2.0 14 1.2 4 4 3.2 多元线性回归分析
运用SPSS软件进行多元线性回归分析,由于各个自变量的单位不一样,自变量前系数的大小并不能说明该因素的重要程度,因此需要将所有变量包括因变量都先转化为标准分,再进行线性回归,此时得到的回归系数就能反映对应自变量的重要程度。详细计算结果如下:
工况1:在构建上部内保护层时,通过改变钻孔倾角、出煤量、孔间距、原始应力及软分层厚度来进行正交分析,以卸压后软、硬分层的卸压程度、变形量为目标值。结合焦作矿区的煤层瓦斯压力及构建内保护层的消突力学判据,以煤层瓦斯抽采率 30%为前提条件,获得不同应力条件下软、硬煤的消突卸压指标进行判定。结合理论判据可知,突出危险性的消除主要通过应力卸压程度进行判定,因此根据正交试验结果,以最薄弱点的卸压程度为目标函数,以钻孔倾角、出煤量、孔间距、原始应力及软煤厚度作为因变量,同时采用多元线性回归的方法对上述数据进行拟合可得到如下拟合方程:
$$ \begin{array}{c} F_{\text {s }}= 88.35-0.106\,92 A+30.739\,5 B-9.440\,55 C- \\ 0.356\,4 D+4.306\,5 E \end{array} $$ (13) $$ \begin{array}{c} F_{\text {h}}= 22.360\,38-0.027\,06 A+7.779\,75 B-2.389\,28 C -\\ 0.090\,2 D+1.089\,917 E \end{array} $$ (14) 式中:A为钻孔倾角,(°);B为出煤量,t/m;C为钻孔间距,m;D为原始应力,MPa;E为软分层厚度,m;F为卸压程度,%。
公式(13)和(14)即为工况1条件下软、硬煤最薄弱点卸压程度的预测模型,将钻孔倾角、出煤量、孔间距、原始应力及软煤厚度各参数进行归一化处理,然后结合预测模型的斜率可计算得到各因素的影响程度,如图6所示。可以看出,各因素对软、硬煤卸压效果的影响程度排列次序是一致的,各影响因素对卸压程度影响大小依次为钻孔间距>出煤量>钻孔倾角>原始应力>软煤厚度,其中以钻孔间距和出煤量影响最为重要。
工况2:当内保护层在下部时,通过改变钻孔倾角、出煤量、孔间距、原始应力及软煤厚度来进行正交分析,以卸压后软、硬分层的卸压程度、变形量为目标值。得到如下拟合方程:
$$ \begin{array}{c} F_{\text {s }}= 89.370\,25-0.097\,39 A+32.453 B-9.412\,2 C- \\ 0.502\,15 D+4.412\,833 E \end{array} $$ (15) $$ \begin{array}{c} F_{\text {h }}= 23.117\,18-0.023\,81 A+8.300\,563 B-2.412\,24 C -\\ 0.143\,01 D+1.260\,708 E \end{array} $$ (16) 公式(15)和(16)即为工况2条件下软、硬煤最薄弱点卸压程度的预测模型,同时将钻孔倾角、出煤量、孔间距、原始应力及软煤厚度各参数进行归一化处理,然后结合预测模型的斜率可计算得到各因素的影响程度,如图7所示。可以看出,各因素对软、硬煤卸压效果的影响程度排列次序是一致的,各影响因素对卸压程度影响大小依次为钻孔间距>出煤量>原始应力>钻孔倾角>软煤厚度,其中以钻孔间距和出煤量影响最为重要。
工况3:当煤层为全软煤时,通过改变钻孔倾角、出煤量、孔间距、原始应力及软煤厚度来进行正交分析,以卸压后软、硬分层的卸压程度、变形量为目标值。得到如下拟合方程:
$$ \begin{array}{c} F_{\text {cs }}= 95.032\,89-0.115\,47 A+36.513\,63 B -\\ 10.859\,9 C-0.370\,59 D \end{array} $$ (17) 公式(17)即为工况3条件下最薄弱点卸压程度的预测模型,同时将钻孔倾角、出煤量、孔间距及原始应力等参数进行归一化处理,然后结合预测模型的斜率可计算得到各因素的影响程度,如图8所示。可以看出,各影响因素对卸压程度影响大小依次为钻孔间距>出煤量>钻孔倾角>原始应力,其中以钻孔间距和出煤量影响最为重要。
工况4:当煤层为全硬煤时,通过改变钻孔倾角、孔间距、原始应力、割缝长度、割缝数量、割缝高度等来进行正交分析,以卸压后煤层的卸压程度、变形量为目标值。得到如下拟合方程:
$$ \begin{split} F_{\text {ch }} =& 12.193\,88-0.007\,97 A - 2.081\,22 C -0.039\,63 D +\\&\qquad 2.929\,02 G+0.263\,529 H+0.168 I \end{split} $$ (18) 式中:G为割缝长度,m;H为割缝数量,个;I为割缝高度,cm;F为最薄弱点卸压程度,%。
式(18)即为工况4条件下最薄弱点卸压程度的预测模型,同时将钻孔倾角、孔间距、原始应力、割缝长度、割缝数量、割缝高度等参数进行归一化处理,然后结合预测模型的斜率可计算得到各因素的影响程度,如图9所示。可以看出,各因素对卸压程度影响大小依次为钻孔间距>割缝长度>割缝数量>割缝高度>原始应力>钻孔倾角,其中以钻孔间距和割缝长度影响最为重要。
4. 内保护层构建技术研究
在煤层构建内保护层能够对煤层进行充分卸压并提高煤层透气性,达到消突目的,通过正交试验给出了不同地质条件煤层构建内保护层的技术参数。根据煤层硬度及软、硬分层分布特征,划分了内保护层构建方法体系。在此基础上,根据内保护层构建方法体系,针对性地研制或完善了内保护层构建技术。不同地质条件煤层内保护层构建技术研发思路如图10所示。
根据煤层硬度及软、硬分层分布特征,划分了内保护层构建方法,在软分层采用水力冲孔或柔性刀具扩孔构建内保护层,在中硬分层采用超高压水射流割缝构建内保护层,在硬及中硬分层采用脉冲磨料水射流割缝构建内保护层。构建技术归纳总结如下:
1)软分层水力冲孔构建内保护层技术。水力冲孔构建内保护层适用于构造煤发育、打钻时喷孔、夹钻、瓦斯压力与瓦斯含量高的低透气性难以抽采的突出煤层和地区。根据煤层地质条件确定卸压空间形态和体积后,结合煤体力学强度,确定冲孔压力、冲孔时间及转速。依据确定的冲孔压力及钻速确定高压柱塞泵及液压马达选型。
2)软分层气动柔性刀具扩孔构建内保护层技术。气动柔性刀具破煤技术应用在全软煤层增透过程中,以高压气体替代水动力带动柔性刀具破煤,避免水力化措施在松软煤层应用时的塌孔、水锁等负面效果[19]。气动柔性刀具扩孔是一种无水化卸压增透技术,从根源上解决了塌孔、堵孔等问题,实现了松软煤层精准卸压。与水力冲孔构建内保护层方法相同,确定气体压力、流量及扩孔时间。依据确定的气体压力及流量确定高压气泵及气动螺杆泵选型。
3)中硬分层超高压水射流割缝构建内保护层技术。超高压水力割缝技术适用于构造煤不发育、坚固性系数大、难以冲孔卸压的突出煤层和地区。通过在煤层中形成连续均匀卸压缝网构建内保护层,调控封孔宽度、深度和分布控制卸压程度和范围。结合超高压水射流破煤理论,确定水射流割缝压力、流量,明确喷嘴选型、布置方式和高压泵组选型。
4)预混合脉冲磨料水射流割缝构建内保护层适用范围。在坚固性系数大于1的坚硬煤层中,超高压水射流无法形成有效缝隙,不能达到煤层均匀卸压的目的。采用脉冲磨料水射流割缝能够在坚硬煤层中形成均匀缝隙,达到构建内保护层的目的。该技术使用的设备包括预混搅拌罐、混砂电机、注浆泵、压力表、高压旋转水尾、钻杆、一体化钻头。一体化钻头包括射流喷嘴和转换阀芯,通过转换阀芯调节射流喷嘴的开关[20]。根据煤层地质条件,结合内保护层构建消突原理,确定割缝空间形态和体积后,依据煤体力学强度,确定射流压力、流量和磨料质量流量,继而确定设备选型。
5. 内保护层构建消突效果评价技术体系
对于焦作矿区单一煤层采用水力冲孔、水力割缝等卸压抽采消突技术措施,可以有效地将大量的煤体冲蚀掉,引起煤层破坏卸压,形成大面积的连通孔洞,释放地应力的同时增加瓦斯解吸有效面积和通道,起到构建煤层内保护层的消突效果。为对单一煤层构建内保护层消突技术的效果进行评价,必须建立行之有效的内保护层构建消突效果评价技术体系。
在焦煤公司下属的九里山矿、中马村矿和古汉山矿进行了构建内保护层消突技术的工程试验,根据前述理论研究及现场考察的相关数据结果,按照构建钻孔的卸压有效影响半径确定内保护层构建钻孔的间距,并严格按照构建钻孔的设计参数施工。
5.1 构建内保护层煤层变形卸压效果评价
传统水力冲孔等水力化措施通常采用单孔卸煤量来考察卸煤效果,但单孔卸煤量是一个粗略的统计参数,且容易导致钻孔之间卸煤不均匀;实际上钻孔卸煤量的多少主要取决于钻孔内煤层段的长度,钻孔施工的角度、位置、长度等参数不同,都会导致钻孔内的煤层段长度出现差异[21]。
在单一煤层内保护层构建消突过程中,采用“密集钻孔、高效卸煤、邻孔串气、均匀卸压”的技术思想构建内保护层达到煤层大面积卸压增透的效果。通过对九里山矿、中马村矿和古汉山矿二1煤层内保护层构建过程中的单孔、组孔、区段的卸煤量进行统计分析,可以看出更为科学的卸煤量统计参数应当是“卸煤率”,即“区段内平均每米煤孔段钻孔对应的卸煤量”,在孔内卸煤过程中应尽量减少各钻孔的卸煤量相对差异性,并在钻孔控制区域内达到一定的卸煤总量。
对九里山矿、中马村矿和古汉山矿二1煤层内保护层构建过程中的单孔、组孔、区段的卸煤量进行统计分析,计算了内保护层构建过程中单孔、组孔、区段控制区域的卸煤率。
九里山矿1613中间底抽巷第一区段的单孔卸煤率为0.56~0.96 t/m,组孔的卸煤率为0.60~0.83 t/m,区段的平均卸煤率为0.67 t/m,该区段内构建钻孔的卸煤总量为控制区域内煤炭储量的33‰。
中马村矿3901运输底抽巷的单孔卸煤率为0.62~1.30 t/m,组孔的卸煤率为0.70~0.80 t/m,区段的平均卸煤率为0.76 t/m,该区段内构建钻孔的卸煤总量为控制区域内煤炭储量的38‰;中马村矿3901 上中间底抽巷的单孔卸煤率为0.53~1.52 t/m,组孔的卸煤率为0.56~0.95 t/m,区段的平均卸煤率为0.75 t/m,该区段内构建钻孔的卸煤总量为控制区域内煤炭储量的37‰。
古汉山矿1606工作面底抽巷第6区段的单孔卸煤率为0.50~1.14 t/m,组孔的卸煤率为0.50~0.82 t/m,区段的平均卸煤率为0.59 t/m,该区段内构建钻孔的卸煤总量为控制区域内煤炭储量的26‰。
根据统计分析1000多个钻孔数据以及矿井在实际生产中统计的相关数据,在目前的钻孔布置方式下,确定达到煤层内保护层构建消突效果的单孔卸煤率为不小于0.60 t/m,并针对具有软分层赋存的层位应加强卸煤措施,加大卸煤量,区段内钻孔卸煤总量不得低于控制区域内煤炭储量的20‰,保证内保护层构建的效果。
5.2 构建内保护层瓦斯抽采效果评价
1)单孔瓦斯抽采情况。基于对内保护层构建后煤层瓦斯抽采钻孔抽采参数的准确监测,对比分析了普通钻孔和内保护层构建消突技术钻孔在瓦斯抽采期间单个钻孔、单组钻孔及区段钻孔的瓦斯抽采情况。
在九里山矿1613中间底抽巷对每个钻孔都进行了整个抽采生命周期的钻孔流量及瓦斯浓度的测定,选择3组具有代表性的典型钻孔进行统计分析,对比钻孔具有相同的钻孔施工参数见表7,各对比钻孔在抽采期间的瓦斯抽采浓度及瓦斯纯量测定结果如图11所示。
表 7 1613中间底抽巷单孔瓦斯抽采对比钻孔参数Table 7. Comparison parameters of single-hole gas extraction drilling for 1613 intermediate bottom extraction roadway对比组 钻孔倾角 煤段长度/m 钻孔总长度/m I 41°20′ 13 35 41°20′ 13 34 II 80°20′ 7 18 80° 7 16.5 III 6° 18 48 6° 18 49 从九里山矿1613中间底抽巷单孔瓦斯抽采情况对比可以看出,内保护层构建消突技术实施后,由于在煤层顶部形成较大的孔洞和蠕变卸压区,导致煤层被充分损伤卸压,下部原生结构煤体膨胀卸压变形,煤层中瓦斯流动的通道显著增加,考察钻孔的瓦斯抽采浓度和纯量显著提高。
与普通钻孔相比较,构建内保护层后穿层抽采钻孔的瓦斯抽采浓度在30%以上的高浓度瓦斯抽采期平均为263 d,而普通钻孔的高浓度瓦斯抽采期平均仅为33 d;构建内保护层后穿层抽采钻孔单孔瓦斯抽采纯量最高为0.38 m3/min,而普通钻孔的单孔瓦斯抽采纯量最高为0.06 m3/min;与普通钻孔相比,构建内保护层穿层抽采钻孔的高浓度瓦斯抽采期提高了8倍,单孔最高瓦斯抽采纯量提高了6.3倍。
2)组孔瓦斯抽采情况。与单孔抽采浓度对比类似,分别在九里山矿1613中间底抽巷和1605底抽巷选取了具有代表性的2组参数类似的瓦斯抽采钻孔,每组各21个钻孔,煤段总长度均为239.0 m,将抽采过程中对比组钻孔的瓦斯浓度及百米纯量的变化进行统计,结果如图12所示。
从图12中组孔的瓦斯抽采情况对比曲线变化情况可以看出,与普通钻孔相比,构建内保护层穿层钻孔能够显著提高煤体中瓦斯运移能力,使瓦斯抽采效率显著提高。在统计组孔瓦斯抽采数据中,构建内保护层穿层钻孔的组孔高浓度瓦斯抽采期达到220 d,最大百米钻孔瓦斯抽采纯量为2.99 m3/(min·hm),且有近180 d百米钻孔抽采纯量高于0.2 m3/(min·hm);而普通钻孔的组孔高浓度瓦斯抽采期仅为40 d,最大百米钻孔瓦斯抽采纯量为0.57 m3/(min·hm),仅有20 d百米钻孔抽采纯量高于0.2 m3/(min·hm)。
瓦斯抽采的450 d内,构建内保护层穿层钻孔的组孔平均百米瓦斯抽采纯量为0.285 m3/(min·hm),普通钻孔平均百米瓦斯抽采纯量为0.038 m3/(min·hm)。以上组孔的瓦斯抽采数据表明,与普通钻孔相比,构建内保护层穿层钻孔的组孔高浓度瓦斯抽采期提高了4.5倍,最大百米瓦斯抽采纯量提高了5.2倍,450 d内平均百米钻孔瓦斯抽采纯量提高了7.5倍。
3)区段钻孔瓦斯抽采情况。对九里山矿1613中间底抽巷第一区段的瓦斯抽采情况进行了统计分析,该区段在进行内保护层构建消突措施实施过程中共施工32组674个钻孔,其中煤孔总长度7918.8 m,在煤层内保护层构建过程中所有钻孔的总出煤量为5346.56 t,该区段控制区域的煤孔段平均卸煤率0.67 t/m,该区段内构建钻孔的卸煤总量为控制区域内煤炭储量的33‰,并选取了以往条件类似、煤孔总长度基本相同的普通钻孔治理区段进行了对比分析,其瓦斯抽采浓度及百米钻孔抽采纯量对比曲线分别如图13和图14所示。
从内保护层构建消突钻孔与普通钻孔区段,在抽采期间内的瓦斯纯量及百米纯量中可以看出,内保护层构建消突措施钻孔比普通钻孔能够显著提高钻孔瓦斯抽采浓度,瓦斯浓度在30%以上的高浓度瓦斯抽采期有210 d,百米瓦斯抽采纯量在0.2 m3/(min·hm)以上的时间有196 d;而普通钻孔在抽采期间,高浓度瓦斯抽采期只有32 d,百米瓦斯抽采纯量在0.2 m3/(min·hm)以上的时间只有38 d;与普通钻孔相比,构建内保护层穿层钻孔的区段钻孔高浓度瓦斯抽采期提高了6.6倍, 百米瓦斯抽采纯量在0.2 m3/(min·hm)以上的时间提高了5.2倍。
综合以上,对内保护层构建消突钻孔和普通钻孔在瓦斯抽采期间单个钻孔、单组钻孔及区段钻孔的瓦斯抽采情况进行对比分析。从所测得的钻孔瓦斯抽采负压和浓度、百米钻孔瓦斯抽采量等参数,可以得出内保护层构建消突技术实施后,由于将大量的煤体冲蚀掉,在煤层顶部形成较大的连通孔洞和卸压区,导致煤层充分卸压,下部原生结构煤体膨胀变形,煤层中瓦斯流动的通道显著增加,钻孔的瓦斯抽采浓度和纯量显著提高,显著提高了煤体中瓦斯运移能力,瓦斯抽采效率显著提升,该消突技术实施后取得了较好的瓦斯抽采效果,减少了突出煤层消突的时间。
5.3 内保护层构建技术消突效果评价体系
1)内保护层构建消突效果的日常评价指标。内保护层构建消突技术主要通过水力冲孔、柔性刀具破煤扩孔、超高压水力割缝、脉冲磨料水射流割缝等技术,达到破硬煤、割缝槽、冲孔洞等强力卸煤的目的,导致煤层破坏卸压,形成大面积的连通孔洞,增加瓦斯解吸有效面积和通道的同时释放地应力,起到构建煤层内保护层的消突效果,因此,在内保护层构建过程中,应按照理论研究及现场考察的卸压有效影响半径,确定内保护层构建钻孔的间距,并严格按照构建钻孔的设计参数施工。
结合对焦作矿区3对矿井典型工作面煤层变形卸压效果的统计分析,提出了更为科学的卸煤量统计参数——“卸煤率”,即“区段内平均每米煤孔段钻孔对应的卸煤量”。经统计九里山矿1613中间底抽巷第一区段的单孔卸煤率为0.56~0.96 t/m,组孔的卸煤率为0.60~0.83 t/m,区段的平均卸煤率为0.67 t/m,该区段内构建钻孔的卸煤总量为控制区域内煤炭储量的33‰;中马村矿3901上运输底抽巷的单孔卸煤率为0.62~1.30 t/m,组孔的卸煤率为0.70~0.80 t/m,区段的平均卸煤率为0.76 t/m,该区段内构建钻孔的卸煤总量为控制区域内煤炭储量的38‰;中马村矿3901上中间底抽巷的单孔卸煤率为0.53~1.52 t/m, 组孔的卸煤率为0.56~0.95 t/m,区段的平均卸煤率为0.75 t/m,该区段内构建钻孔的卸煤总量为控制区域内煤炭储量的37‰;古汉山矿1606工作面底抽巷第六区段的单孔卸煤率为0.50~1.14 t/m,组孔的卸煤率为0.50~0.82 t/m,区段的平均卸煤率为0.59 t/m,该区段内构建钻孔的卸煤总量为控制区域内煤炭储量的26‰。根据研究分析数据以及各矿井在实际生产中统计数据,将内保护层构建钻孔的间距应按考察研究的钻孔卸压有效影响半径确定、达到煤层内保护层构建消突效果的卸煤率不小于每米煤孔0.60 t,在实施过程中应尽量减少各钻孔的卸煤量相对差异性,且钻孔卸煤总量不得低于控制区域内煤炭储量的20‰,作为日常评价指标,并应对软分层赋存的层位加强卸煤措施,加大卸煤量,增大卸煤率,确保内保护层构建后各钻孔的连通性和煤层的全域卸压,保证内保护层构建的效果。
2)内保护层构建煤层卸压均匀性及多物理场叠加区评价指标。在单一煤层内保护层构建消突过程中,采用“密集钻孔、孔内卸煤、孔孔贯通、均匀卸压”的技术思想构建内保护层达到煤层大面积卸压增透的效果。内保护层构建是钻孔连片卸压形成的,容易在各钻孔影响范围的叠加区形成消突效果的薄弱区,根据传统钻孔经验易产生应力集中、卸压程度低、煤层渗透性差等现象。
为保证内保护层的构建效果,通过对内保护层构建消突技术实施过程中多物理场叠加区的应力、变形、透气性等参数的考察和分析,提出内保护层构建消突效果评价的辅助指标I:煤层垂向膨胀变形量不低于3‰,用来表征煤层中内保护层的形成和煤体的膨胀变形程度。
同时,通过对内保护层构建后煤层瓦斯抽采情况的进行实时跟踪观测,提出内保护层构建消突效果评价的辅助指标Ⅱ:单孔初始抽采浓度不低于30%,单孔、组孔、区段钻孔的瓦斯抽采浓度在30%以上的高浓度抽采期不低于90 d,瓦斯抽采率差异性系数低于10%,瓦斯抽采累计时间不少于设计预抽期,最终以残存瓦斯压力和瓦斯含量均小于“双六”指标进行评判。
通过对内保护层构建消突技术煤层变形卸压、瓦斯抽采效果的评价以及对构建内保护层后多物理场叠加区卸压效果的考察验证,在充分结合构建内保护层消突技术特点及现场工程实际的基础上,提出了单一煤层构建内保护层消突技术的效果评价技术体系和指标,如图15所示。
6. 结 论
1)提出了强突出单一煤层构建内保护层消突技术思想,阐释了构建内保护层消突的科学内涵。构建内保护层技术以煤层卸压消突为目标,以地质勘探为基础,以内保护层消突技术参数和高效构建技术装备为依托,以消突效果考察和系统保障体系为保障,实现穿层钻孔的“布孔均匀密集、卸煤科学高效、邻孔串气相透、卸压均匀充分”,统筹了煤层区域高效均匀卸压与瓦斯抽采,提出适用于强突出单一煤层的主动构建内保护层的消突理念。通过多个采掘工作面构建内保护层的工程实践,说明此方法是科学、安全和高效的。
2)揭示了构建内保护层消突机理及消突力学评判依据。从地应力和瓦斯压力的角度,基于发生煤与瓦斯突出的能量方程,确定了煤层开采时发生煤与瓦斯突出能量条件。明确了不同硬度煤体突出主控因素和卸压消突力学判据:硬煤卸压消突效果优于软煤,硬分层仅靠掏煤卸压到一定程度即可达到消突效果,软分层卸压配合瓦斯抽采也可达到消突;试验研究条件下,f为0.5时,当卸压程度达到58%以上,可实现煤层消突,在瓦斯预抽率15%前提下,卸压程度达到10%即可实现消突;f为0.2时,在瓦斯预抽率30%前提下,卸压程度达到40% 即可实现消突效果。
3)形成了内保护层构建方法体系。根据煤体软硬分层结构划分4种工况,研究了不同煤层赋存、钻孔间距、出煤量、软分层厚度、钻孔角度、原始应力、水力割缝参数等条件下,构建内保护层时的卸压程度及煤层变形量的主控因素及变化规律。通过正交实验模拟分析了4种工况不同内保护层构建条件下,钻孔施工参数的卸压效果,结合多元线性回归的方法,分析了不同影响因素对消突效果的影响程度,构建了钻孔控制范围内最薄弱作用点的卸压程度预测模型,形成了内保护层构建方法体系,为构建内保护层消突技术提供了施工参数确定依据。
4)建立了内保护层消突效果评价体系。针对不同瓦斯地质条件,在确保合理钻孔布置前提下,提出了采用“卸煤率”作为内保护层构建效果的日常考察指标,煤层垂向膨胀变形量、瓦斯抽采浓度和抽采差异性等为辅助评价指标;每米煤孔的卸煤率为不小于0.60 t/m,钻孔卸煤总量不得低于控制区域内煤炭储量的20‰。
-
表 1 废弃矿井采动稳定区瓦斯资源开发潜力评价初始指标体系及选取依据
Table 1 Initial index system and selection basis of gas resource development potential evaluation in abandoned mine mining stability area
目标层 准则层 指标层 选取依据 参数阈值(范围) 废弃矿井采动稳定区瓦斯资源开发潜力评价指标体系 瓦斯资源量因素 采空区遗煤瓦斯资源量 遗煤越多,瓦斯量越大 潘一矿可采煤层的瓦斯总量为61.81×
108 m3,属于中型煤层气田采空区未动用煤瓦斯资源量 未动用煤越多,瓦斯量越大 邻近层煤瓦斯资源量 邻近层煤越多,瓦斯量越大 煤矿面积 越大,遗煤越多,瓦斯含量越大 54.6717 km2煤炭采出率 越多,遗煤越少 — 煤厚 越厚,瓦斯含量越高 矿井含可采煤层13层,可采煤层平均总厚 30.38 m 资源丰度 越大,价值越大 — 稳定性因素 采空区封闭条件 越好,瓦斯越难逸散 良好 巷道支护完好程度 越好,瓦斯储存空间越大 良好 煤层埋深 越深,地应力越大,瓦斯保存条件越好 −500~−800 围岩性质 越硬,岩层透气性越差,瓦斯保存条件越好 砾岩、细砂岩、粗砂岩、粉砂岩 矿井涌水量 越大,抽采成本越多 144.1~177.7 m3/h 瓦斯抽采条件因素 废弃时间 越长,解吸瓦斯越多 6 a 煤层孔裂隙发育程度 越好,孔隙率越高,越有利于抽采 各可采煤层孔隙率为 0.7%~8.3% 煤体破坏程度 越高,渗透率越差,抽采时易堵孔 破坏程度较高 地面钻孔甲烷浓度 越大,瓦斯可利用率越高 — 环境风险 井下积水、遗煤自燃,瓦斯逸散严重,不利于开发利用 有自燃倾向性、自燃等级 Ⅱ 级 社会影响 瓦斯抽采,有利于减少瓦斯在井下聚集的浓度,减少瓦斯灾害 — 可持续性 瓦斯资源量越多、稳定性越好等,瓦斯抽采的可持续性越长 — 经济条件因素 矿井现有设备条件 越好,投入成本越低 良好 矿井巷道开拓状态 越好,抽采提升速度越快 — 瓦斯价格 越高,成本回升越快 2.5元/m3(受市场影响较大) 抽采工艺 越好,抽采效率越好 “L”型钻井抽采 管网阻力 越大,抽采率越低 — 表 2 DEMATEL-ISM模型结果
Table 2 DEMATEL-ISM Model Results
指标 影响度 被影响度 中心度 原因度 指标 影响度 被影响度 中心度 原因度 C1 0.3675 1.2837 1.6512 − 0.9162 C10 0.9944 0.4926 1.4870 0.5018 C2 0.3675 1.0597 1.4272 − 0.6922 C11 0.3852 0.2318 0.6170 0.1534 C3 0.3675 1.0597 1.4272 − 0.6922 C12 0 2.1601 2.1601 − 2.1601 C4 0.7888 0 0.7888 0.7888 C13 0.8835 0.5572 1.4407 0.3263 C5 1.3839 0.3672 1.7511 1.0167 C14 0.0483 0.5820 0.6303 − 0.5337 C6 0.4641 0.5251 0.9892 −0.0610 C15 0.3015 0.1420 0.4435 0.1595 C7 1.4192 0 1.4192 1.4192 C16 0.0639 0.3590 0.4229 − 0.2951 C8 0.4070 0.4926 0.8996 − 0.0856 C17 0.3423 0.3998 0.7421 − 0.0575 C9 1.0085 0.0622 1.0707 0.9463 C18 0.4645 0.2829 0.7474 0.1816 表 3 判断矩阵标度含义
Table 3 3Meaning of Judgment Matrix Scale
标度 含义 0.9 ri元素极端重要于rj元素 0.8 ri元素强烈重要于rj元素 0.7 ri元素明显重要于rj元素 0.6 ri元素稍微重要于rj元素 0.5 ri元素与rj元素同样重要 互补数 若ri比rj相比较判断为rij,那么rj比ai相比较判断为1−rij 表 4 判断矩阵B-G及特征向量
Table 4 Judgment matrix B-G and eigenvectors
B-G B1 B2 B3 B4 $ {G_i} $ B1 0.5000 0.5304 0.5469 0.5834 0.2917 B2 0.4696 0.5000 0.5166 0.5536 0.2583 B3 0.4531 0.4834 0.5000 0.5371 0.2417 B4 0.4166 0.4464 0.4629 0.5000 0.2083 表 5 基于FAHP的指标层对目标层的各指标权重值
Table 5 Weight values of each indicator in target layer based on FAHP indicator layer
C-G B1( 0.2917 )B2( 0.2583 )B3( 0.2417 )B4( 0.2083 )$ {G_i} $ C1 0.2500 — — — 0.0729 C2 0.2417 — — — 0.0705 C3 0.2250 — — — 0.0656 C4 0.2833 — — — 0.0826 C5 — 0.2833 — — 0.0732 C6 — 0.2500 — — 0.0646 C7 — 0.2417 — — 0.0624 C8 — 0.2250 — — 0.0581 C9 — — 0.2000 — 0.0483 C10 — — 0.1533 — 0.0371 C11 — — 0.1500 — 0.0363 C12 — — 0.1600 — 0.0387 C13 — — 0.1600 — 0.0387 C14 — — 0.1767 — 0.0427 C15 — — — 0.2833 0.0590 C16 — — — 0.2333 0.0486 C17 — — — 0.2583 0.0538 C18 — — — 0.2250 0.0469 表 6 基于Entropy的指标层对目标层的各指标权重值
Table 6 Weight values of each indicator in target layer based on Entropy’s indicator layer
C-G B1( 0.2917 )B2( 0.2583 )B3( 0.2417 )B4( 0.2083 )C1 0.8138 — — — C2 0.7918 — — — C3 0.7934 — — — C4 0.8222 — — — C5 — 0.8280 — — C6 — 0.0025 — — C7 — 0.6070 — — C8 — 0.8263 — — C9 — — 0.8615 — C10 — — 0.8210 — C11 — — 0.8147 — C12 — — 0.8201 — C13 — — 0.0025 — C14 — — 0.0025 — C15 — — — 0.3967 C16 — — — 0.8401 C17 — — — 0.6830 C18 — — — 0.8264 表 7 基于博弈论的指标层对目标层的各指标权重值
Table 7 Weight values of each indicator combination in target layer based on game theory for indicator layer
指标 FAHP权重 熵值法权重 组合权重 C1 0.0729 0.8138 0.3414 C2 0.0705 0.7918 0.3319 C3 0.0656 0.7934 0.3293 C4 0.0826 0.8222 0.3506 C5 0.0732 0.8280 0.3467 C6 0.0646 0.0025 0.0421 C7 0.0624 0.6070 0.2598 C8 0.0581 0.8263 0.3365 C9 0.0483 0.8615 0.3430 C10 0.0371 0.8210 0.3212 C11 0.0363 0.8147 0.3184 C12 0.0387 0.8201 0.3219 C13 0.0387 0.0025 0.0256 C14 0.0427 0.0025 0.0281 C15 0.0590 0.3967 0.1814 C16 0.0486 0.8401 0.3354 C17 0.0538 0.6830 0.2818 C18 0.0469 0.8264 0.3294 表 8 指标评估等级划分
Table 8 Classification of Index Evaluation Levels
评估指标 差Ⅰ 一般Ⅱ 良好Ⅲ 优秀Ⅳ 采空区遗煤瓦斯资源量C1 较差[0,60) 一般[60,70) 较高[70,80) 高[80,100] 采空区未动用煤瓦斯资源量C2 较差[0,60) 一般[[60,70) 较高[70,80) 高[80,100] 邻近层煤瓦斯资源量C3 较差[0,60) 一般[[60,70) 较高[70,80) 高[80,100] 煤矿面积C4 较差[0,60) 一般[[60,70) 较高[70,80) 高[80,100] 采空区封闭条件C5 较差[0,60) 一般[[60,70) 较高[70,80) 高[80,100] 巷道支护完好程度C6 较差[0,60) 一般[[60,70) 较高[70,80) 高[80,100] 煤层埋深C7 较差[0,60) 一般[[60,70) 较高[70,80) 高[80,100] 围岩性质C8 较差[0,60) 一般[[60,70) 较高[70,80) 高[80,100] 废弃时间C9 较差[0,60) 一般[[60,70) 较高[70,80) 高[80,100] 煤层孔裂隙发育程度C10 较差[0,60) 一般[[60,70) 较高[70,80) 高[80,100] 煤体破坏程度C11 高[0,60) 较高[60,70) 一般[70,80) 低[80,100] 地面钻孔甲烷浓度C12 少[0,60) 一般[60,70) 较多[70,80) 多[80,100] 环境风险C13 高[0,60) 较高[60,70) 一般[70,80) 低[80,100] 社会影响C14 少[0,60) 一般[60,70) 较多[70,80) 多[80,100] 矿井现有设备条件C15 少[0,60) 一般[60,70) 较多[70,80) 多[80,100] 矿井巷道开拓状态C16 较差[0,60) 一般[60,70) 较好[70,80) 好[80,100] 瓦斯价格C17 低[0,60) 一般[60,70) 较高[70,80) 高[80,100] 抽采工艺C18 较差[0,60) 一般[60,70) 较好[70,80) 好[80,100] 表 9 二级评估等级关联度
Table 9 Level 2 evaluation level correlation
因素层B1−4 指标层 权重 评估指标关联度Kj 所属等级 j=1 j=2 j=3 j=4 B1 C1 0.3414 −0.425 −0.233 0.3 −0.115 Ⅲ C2 0.3319 −0.4 −0.2 0.4 −0.143 Ⅲ C3 0.3293 −0.079 0.3 −0.167 −0.327 Ⅱ C4 0.3506 −0.525 −0.367 −0.05 0.05 Ⅳ B2 C5 0.3467 −0.350 −0.133 0.400 −0.188 Ⅲ C6 0.0421 −0.375 −0.167 0.500 −0.167 Ⅲ C7 0.2598 −0.500 −0.667 0 0 Ⅲ C8 0.3365 0.083 −0.100 −0.250 −0.357 Ⅰ B3 C9 0.3430 −0.300 −0.067 0.200 −0.222 Ⅲ C10 0.3212 −0.500 −0.667 0 0 Ⅲ C11 0.3184 −0.150 0.400 −0.105 −0.292 Ⅱ C12 0.3219 −0.125 0.500 −0.125 −0.300 Ⅱ C13 0.0256 −0.500 −0.333 0 0 Ⅲ C14 0.0281 −0.05 0.2 −0.174 −0.391 Ⅱ B4 C15 0.1814 −0.300 −0.067 0.200 −0.222 Ⅲ C16 0.3354 −0.250 0 0 −0.250 Ⅲ C17 0.2818 0.067 −0.241 −0.833 −0.353 Ⅰ C18 0.3294 −0.375 −0.167 0.500 −0.167 Ⅲ 表 10 一级评估等级关联度
Table 10 Level 1 evaluation level correlation
目标层 指标层 权重 评估指标关联度Kj 所属等级 j=1 j=2 j=3 j=4 G B1 0.2917 − 0.4879 − 0.1758 0.1627 − 0.1769 Ⅲ B2 0.2583 − 0.2391 − 0.2601 0.0756 − 0.1923 Ⅲ B3 0.2417 − 0.3657 0.0482 − 0.0100 − 0.2767 Ⅱ B4 0.2083 − 0.2429 − 0.1351 − 0.0338 − 0.2786 Ⅲ — — − 0.3431 − 0.1350 0.0575 − 0.2262 Ⅲ -
[1] 袁亮,赵毅鑫. “我国关闭/废弃矿井资源精准开发利用研究” 专题特邀主编致读者[J]. 煤炭学报,2022,47(6):2125−2126. YUAN Liang,ZHAO Yixin. Invited editor-in-chief of “research on precise development and utilization of closed/abandoned mine resources in China” to readers[J]. Journal of China Coal Society,2022,47(6):2125−2126.
[2] 袁亮. 关闭煤矿资源开发利用迫在眉睫[J]. 新能源经贸观察,2018(7):50. YUAN Liang. It is urgent to close coal mine resources and develop and utilize them[J]. Energy Outlook,2018(7):50.
[3] 袁亮. 我国煤炭资源高效回收及节能战略研究[M]. 北京:科学出版社,2017. [4] KARACAN C Ö,WARWICK P D. Assessment of coal mine methane (CMM) and abandoned mine methane (AMM) resource potential of longwall mine panels:Example from Northern Appalachian Basin,USA[J]. International Journal of Coal Geology,2019,208:37−53. doi: 10.1016/j.coal.2019.04.005
[5] YOUNGER P L. Abandoned coal mines:From environmental liabilities to low-carbon energy assets[J]. International Journal of Coal Geology,2016,164:1−2. doi: 10.1016/j.coal.2016.08.006
[6] KUBIT O E,PLUHAR C J,DE GRAFF J V. A model for prioritizing sites and reclamation methods at abandoned mines[J]. Environmental Earth Sciences,2015,73(12):7915−7931. doi: 10.1007/s12665-014-3949-3
[7] 孟召平,师修昌,刘珊珊,等. 废弃煤矿采空区煤层气资源评价模型及应用[J]. 煤炭学报,2016,41(3):537−544. MENG Zhaoping,SHI Xiuchang,LIU Shanshan,et al. Evaluation model of CBM resources in abandoned coal mine and its application[J]. Journal of China Coal Society,2016,41(3):537−544.
[8] 刘小磊,闫江伟,刘操,等. 废弃煤矿瓦斯资源估算与评价方法构建及应用[J]. 煤田地质与勘探,2022,50(4):45−51. LIU Xiaolei,YAN Jiangwei,LIU Cao,et al. Construction and application of gas resource estimation and evaluation method in abandoned coal mines[J]. Coal Geology & Exploration,2022,50(4):45−51.
[9] 胡长勤,刘坤鹏. 废弃煤矿瓦斯资源量预测模型的构建与应用[J]. 能源与环保,2020,42(9):1−6,12. HU Changqin,LIU Kunpeng. Establishment and application of prediction model for gas resources in abandoned coal mines[J]. China Energy and Environmental Protection,2020,42(9):1−6,12.
[10] 刘文革,张康顺,韩甲业,等. 废弃煤矿瓦斯开发利用技术与前景分析[J]. 中国煤层气,2016,13(6):3−6. LIU Wenge,ZHANG Kangshun,HAN Jiaye,et al. Technology and prospect analysis of AMM development and utilization[J]. China Coalbed Methane,2016,13(6):3−6.
[11] 文光才,孙海涛,李日富,等. 煤矿采动稳定区煤层气资源评估方法及其应用[J]. 煤炭学报,2018,43(1):160−167. WEN Guangcai,SUN Haitao,LI Rifu,et al. Assessment method and application of coalbed methane resources in coal mining stability area[J]. Journal of China Coal Society,2018,43(1):160−167.
[12] 彭金刚. 鸡西、鹤岗矿区封闭采空区煤层气资源评价[D]. 徐州:中国矿业大学,2017. PENG Jingang. Evaluation of coalbed methane resources in closed goaf of Jixi and Hegang mining areas[D]. Xuzhou:China University of Mining and Technology,2017.
[13] 韩保山,李健武,董敏涛. 用下降曲线估算废弃矿井煤层气资源量[J]. 中国煤田地质,2005,17(5):37−39,46. HAN Baoshan,LI Jianwu,DONG Mintao. Estimation of abandoned mine CBM resources through drop-down curves[J]. Coal Geology of China,2005,17(5):37−39,46.
[14] CHEN S,SUN L,ZHANG Y H. Research on the distributing law of the gas in the gob area based on flow-tube model[J]. Procedia Engineering,2011,26:1043−1050. doi: 10.1016/j.proeng.2011.11.2272
[15] 李树刚,程皓,潘红宇,等. 崔家沟煤矿采空区瓦斯抽采效果评价模型[J]. 西安科技大学学报,2020,40(1):11−17,87. LI Shugang,CHENG Hao,PAN Hongyu,et al. Evaluation model of gas drainage effect in goaf of Cuijiagou coal mine[J]. Journal of Xi’an University of Science and Technology,2020,40(1):11−17,87.
[16] 贾志超,黄华州,黄绍博,等. 基于AHP-熵权法的采动区煤层气开发潜力评价[J]. 煤田地质与勘探,2021,49(2):117−124. JIA Zhichao,HUANG Huazhou,HUANG Shaobo,et al. Evaluation of the development potential of the coalbed methane resources in mining area based on AHP-entropy method[J]. Coal Geology & Exploration,2021,49(2):117−124.
[17] 冯玉龙,周林元,王乾,等. 基于多层次模糊数学法的煤层气井产能综合评价模型[J]. 煤田地质与勘探,2021,49(2):125−132. FENG Yulong,ZHOU Linyuan,WANG Qian,et al. Multi-level fuzzy mathematics-based comprehensive evaluation model of CBM well productivity[J]. Coal Geology & Exploration,2021,49(2):125−132.
[18] 李丹,苏现波. 煤与煤层气资源开发全过程阶段划分及其开发效果评价[J]. 煤炭科学技术,2023,51(3):137−147. LI Dan,SU Xianbo. Stage division and development effect evaluation of whole process of coal and coalbed methane resources development[J]. Coal Science and Technology,2023,51(3):137−147.
[19] 王红杰,卞孝东,邓晓伟,等. 基于模糊数学理论的煤矿地下空间开发利用适宜性评价:以白源煤矿为例[J]. 西北地质,2021,54(4):156−170. WANG Hongjie,BIAN Xiaodong,DENG Xiaowei,et al. Suitability evaluation of underground space development and utilization of coal mine based on fuzzy mathematics theory:Taking baiyuan coal mine as an example[J]. Northwestern Geology,2021,54(4):156−170.
[20] 杨明,朱林,全星苑,等. 废弃煤矿采空区瓦斯抽采可行性评价模型构建及应用[J]. 煤田地质与勘探,2022,50(7):107−117. doi: 10.12363/issn.1001-1986.21.11.0686 YANG Ming,ZHU Lin,QUAN Xingyuan,et al. Construction and application of feasibility evaluation model for gas extraction in goafs of abandoned coal mines[J]. Coal Geology & Exploration,2022,50(7):107−117. doi: 10.12363/issn.1001-1986.21.11.0686
[21] 李开心. 废弃矿井瓦斯资源开发利用潜力分析及风险评价[D]. 焦作:河南理工大学,2022. LI Kaixin. Potential analysis and risk assessment of gas resource development and utilization in abandoned mines [D]. Jiaozuo:Henan University of Technology,2022.
[22] 张凯鑫. 阳泉矿区关闭煤矿剩余资源特征及其利用[D]. 徐州:中国矿业大学,2023. ZHANG Kaixin. Characteristics and utilization of remaining resources from closed coal mines in yangquan mining area [D]. Xuzhou:China University of Mining and Technology,2023.
[23] 周一沫. 废弃煤矿地下空间开发利用评价研究[D]. 重庆:重庆大学,2020. ZHOU Yimo. Study on evaluation of underground space development and utilization in abandoned coal mines[D]. Chongqing:Chongqing University,2020.
[24] 樊燕燕,王瑞. 物元可拓法在地铁盾构施工安全风险评估中的应用[J]. 安全与环境学报,2023,23(6):1779−1790. FAN Yanyan,WANG Rui. Application of matter-element extension method in safety risk assessment of subway shield construction[J]. Journal of Safety and Environment,2023,23(6):1779−1790.
[25] 王福生,高湛翔,董宪伟,等. 基于物元可拓理论的化工企业风险评价[J]. 安全与环境学报,2021,21(6):2401−2406. WANG Fusheng,GAO Zhanxiang,DONG Xianwei,et al. Risk assessment of the chemical enterprises based on the matter-element extension theory[J]. Journal of Safety and Environment,2021,21(6):2401−2406.
[26] 吴丹红,张美霞,张汉斌,等. 基于可拓学的地铁车站深基坑施工安全评价[J]. 安全与环境学报,2019,19(3):761−766. WU Danhong,ZHANG Meixia,ZHANG Hanbin,et al. On the safety evaluation for the deep foundation pit of the subway stations based on the theory of extenics[J]. Journal of Safety and Environment,2019,19(3):761−766.
[27] 傅宁,邵月龄. 基于物元可拓模型的机场飞行区安全风险评价及影响程度分析[J]. 安全与环境学报,2023,23(12):4247−4254. FU Ning,SHAO Yueling. Safety risk assessment and influence degree analysis of the flight area based on the matter-element extension theory[J]. Journal of Safety and Environment,2023,23(12):4247−4254.
[28] 孙中博,赵毅鑫,任建东. 废弃矿井地下空间再利用根源因素识别及评价模型构建[J]. 中国矿业大学学报,2024,53(3):585−599. SUN Zhongbo,ZHAO Yixin,REN Jiandong. Identification and evaluation model construction of root factors for reuse of underground space in abandoned mines[J]. Journal of China University of Mining and Technology,2024,53(3):585−599.
[29] 陈红,周继彪,王建军,等. 公路隧道运行环境安全评价指标与方法[J]. 长安大学学报(自然科学版),2013,33(4):54−61. CHEN Hong,ZHOU Jibiao,WANG Jianjun,et al. Safety evaluation indexes and method for traffic environment of highway tunnels[J]. Journal of Chang'an University ( Natural Science Edition),2013,33(4):54−61.
[30] 陈赟,刘湘慧,陈玉斌. 基于决策实验室分析法及解释结构模型的交通隧道施工事故致因链构建[J]. 科学技术与工程,2021,21(35):15222−15229. CHEN Yun,LIU Xianghui,CHEN Yubin. Construction of the cause chain of traffic tunnel construction accidents based on DEMATEL-ISM[J]. Science Technology and Engineering,2021,21(35):15222−15229.
-
期刊类型引用(9)
1. 穆朝元,武腾飞,张朝文,刘勇,熊俊杰. 全封闭负压连抽高效防喷装置研究及应用. 煤炭工程. 2025(01): 218-224 . 百度学术
2. 樊世星,宋金锁,程小蛟,米万升,李日军,时玉坤,乔季洁. 中低压力N_2/CO_2驱替煤层CH_4非线性渗流特征研究. 煤矿安全. 2025(04): 9-18 . 百度学术
3. 马腾,王峰,赵明广,陈森军,郭志远. 保护层下巷道掘进围岩应力演化规律研究. 煤炭技术. 2025(05): 39-44 . 百度学术
4. 刘浩东. 高瓦斯矿井综采工作面瓦斯治理技术研究. 西部探矿工程. 2024(03): 53-55 . 百度学术
5. 张晓铭. 基于定向水力掏煤的层内“保护层”强化抽采技术. 陕西煤炭. 2024(06): 133-137 . 百度学术
6. 卢庆红,汪强,高少鹏,申晟,葛家君. 月亮田煤矿选煤厂煤泥浓缩技术改造及应用. 煤炭加工与综合利用. 2024(05): 29-31+36 . 百度学术
7. 周爱星,张颖异,马钱钱. 煤矿瓦斯抽采系统新型自动排水连孔工艺研究. 煤炭与化工. 2024(07): 111-114 . 百度学术
8. 李江涛,王飞. 神经网络结合机器学习的煤与瓦斯突出量和危险性等级预测组合模型. 中国矿业. 2024(S2): 176-184 . 百度学术
9. 王宏伟,段晨烨. 王庄煤矿动压巷道顺层钻孔二次注浆提浓技术研究. 煤炭科技. 2024(05): 158-162 . 百度学术
其他类型引用(6)