Research on gel plugging and fire prevention technology for roadway excavation in residual coal re-mining area
-
摘要:
采空区中遗留煤炭自燃现象直接影响着煤矿的安全高效开采，而采用堵漏材料封闭漏风通道是解决采空区煤自燃问题的有效方法。以运河煤矿F1301复采工作面为工程背景，针对在采空区中掘进巷道时遗留松散煤体的自燃隐患问题，制定了松散煤体凝胶堵漏防灭火技术实施方案，采用COMSOL数值模拟方法研究了不同注胶时间与压力下的凝胶扩散规律，并进行了工程应用。结果表明：注胶时间和压力与凝胶扩散范围成正比，综合考虑注胶效果与成本，选取10 MPa、140~160 min为最佳的注胶压力与时间，注胶钻孔间距设置为5 m；实施凝胶堵漏防灭火方案后，巷道两侧观测钻孔监测数据显示，标志气体CO体积分数稳定于24×10−6以下，凝胶堵漏防灭火效果良好，有效消除了采空区遗留煤炭自然发火隐患。
Abstract:The phenomenon of spontaneous combustion of coal left in the goaf directly affects the safe and efficient mining of coal mines, and the use of plugging materials to close the air leakage channel is an effective method to solve the problem of spontaneous combustion of coal in the goaf. Taking the F1301 re-mining face of Yunhe Coal Mine as the engineering background, aiming at the hidden danger of spontaneous combustion of loose coal left in the tunneling roadway in the goaf, the implementation scheme of gel plugging and fire extinguishing technology for loose coal was formulated. The COMSOL numerical simulation method was used to study the gel diffusion law under different injection time and pressures, and the engineering application was carried out. The results show that the injection time and pressure are proportional to the gel diffusion range. Considering the injection effect and cost, 10 MPa and 140-160 min are selected as the best injection pressure and time, and the distance between the injection holes is set to 5 m. After the implementation of the gel plugging fire prevention and extinguishing scheme, the monitoring data of the observation boreholes on both sides of the roadway show that the volume fraction of the iconic gas CO is stable below 24×10−6, and the gel plugging fire prevention and extinguishing effect is good, which effectively eliminates the hidden danger of spontaneous combustion of coal left in the goaf.
-
煤炭开采机械化程度不断提高，粉尘量也明显增多，在煤炭掘进和开采过程中都会产生大量的煤尘、矽尘等混合性粉尘，会给作业人员的生命健康带来巨大的危害。因此对存在粉尘的长期作业场所，粉尘治理成了越来越重视的问题。
随着分子模拟技术的发展，经常使用分子动力学模拟研究分子之间的微观相互作用[1-2]。高正阳等[3]研究了4种不同煤阶煤分子表面水分子吸附的微观机理；XU等[4]和HUANG等[5]分别研制了针对建筑工地粉尘和岩石粉尘的新型环保抑尘剂，并通过试验得出新型抑尘剂具有突出的抑尘效果；CRAWFORD等[6]研究了表面活性剂对不同煤表面特性的影响，发现添加表面活性剂可以明显降低煤的接触角，提高煤的润湿效果；ZHANG等[7]通过分子动力学模拟研究了表面活性剂对褐煤吸附行为和润湿性的影响，发现表面活性剂可以改变水与褐煤之间的相互作用力；陈璐瑶[8]研究了表面活性剂结构对SiO2粉尘润湿性，通过分子动力学模拟，得出SiO2表面基团性质是影响SiO2润湿性的根本原因；孟筠青等[9]探究不同种类表面活性剂及其复配对SiO2润湿性的影响，同时结合沉降试验，找出SiO2粉尘润湿效果最好的单体表面活性剂及复配溶液；邱玄等[10]通过分子动力学模拟了 Pluronic F-68 在萤石和石英表面的吸附行为，结果显示 Pluronic F-68在石英上的吸附能力更强。
上述研究表明,分子模拟已成为研究润湿和吸附性能的重要方法[11-12]。但研究都是针对单一性粉尘研究润湿性。为此，借助于Materials studio分子模拟软件构建模型，采用分子动力学模拟和量子化学的方法，以表面活性剂CTAB为例，研究褐煤煤尘、矽尘的润湿机理的差异性。
1. 模拟方法
1.1 模拟体系构建
采用Materials Studio 分子模拟软件，构建了水−表面活性剂−褐煤煤尘/矽尘模型。初始模型如图1所示。
选用Wender经典褐煤模型，其煤分子结构如图1（a）所示；矽尘模型为SiO2超晶胞（图1（b）），添加1.5 nm的真空层以消除周期性的影响；表面活性剂选用十六烷基三甲基溴化铵（CTAB）（图1（c））。利用MS中Amorphous Cell 模块构建各个盒子，盒子详细参数见表1。
表 1 盒子或超晶胞参数Table 1. Box or supercell parameters名称 长/nm 宽/nm 高/nm 分子数量 α、β、γ/（°） 密度/ （g·cm−3） 水盒子 2.946 0 3.241 2 1.566 4 500 90 1.000 CTAB盒子 2.946 0 3.241 2 0.239 7 5 90 1.322 煤分子盒子 2.946 0 3.241 2 1.342 6 8 90 1.350 SiO2超晶胞 2.946 0 3.241 2 2.288 3 — 90 2.200 注：α、β、γ为三维空间中与x、y、z轴的夹角。 通过build layer构成水−表面活性剂−褐煤/矽尘（SiO2）2个模型，并在体系上方都添加了0.5 nm真空层。模拟体系的构建如图2所示，模拟平衡后水的密度场分布图如图3所示。
![]()
图 3 模拟平衡后水的密度场分布图Figure 3. Density field distribution of water after simulated equilibrium1.2 计算模拟方法
1）分子动力学模拟。运用Forcite模块对褐煤分子、水分子、CTAB分子和SiO2超晶胞进行结构优化。首先采用Geometry Optimization任务中通过Smart算法进行局部结构优化。然后，利用Amorphous Cell 模块构建500个水分子盒子、5个表面活性剂盒子、8个褐煤分子盒子，通过build layer构成模型，在模型上方添加0.5 nm真空层。然后，为了得到模型全局最小和稳定的能量，在Anneal任务中采用NVT系综（温度为298 K）进行退火动力学模拟。最后，在Dynamics任务中采用NVT系综进行1 ns的分子动力学模拟，其中力场采用COMPASSⅡ，电荷采用Forcefield assigned。
2）量子化学模拟。量子力学计算在DMol3模块中完成[13]，在广义梯度近似(GGA)[14]下，利用BLYP泛函模拟电子交换和相关相互作用[15]，模拟精度为Fine，其中SCF收敛精度为0.000 001。选择所有由D极化函数（DNP）增强的电子双数值原子轨道作为基础集，并将轨道截止质量定为良好。smearing的值选择为0.005，并且在计算过程中不考虑基集叠加误差（BSSE）[16]，因为在DMol3中实现的数值基集可最小化甚至消除基集叠加误差。
2. 模拟结果
2.1 相互作用能
在模拟体系中，水分子的形态和流动的差异性主要是受表面活性剂水溶液和褐煤分子或者SiO2相互作用的结果。相互作用能可以数字化体现分子之间的相互作用的强度，如果相互作用能为负值，则表明吸附过程是自发进行的；若相互作用能为正值，则反之。相互作用能的绝对值越大时，分子之间的相互作用也越强。通过分子动力学模拟，计算体系相互作用能公式如下：
$$ {E_{{{\mathrm{inter}}} }} = {E_{{\mathrm{total}}}} - \left( {{E_{\mathrm{i}}} + {E_{{\mathrm{w + s}}}}} \right) $$ (1) 式中：$ {E_{{{\mathrm{inter}}} }} $为活性剂水溶液与煤或SiO2的相互作用能；$ {E_{{\mathrm{total}}}} $为体系内的总能量；$ {E_{\mathrm{i}}} $为体系中褐煤分子或者SiO2的能量；$ {E_{{\mathrm{w + s}}}} $为水和表面活性剂的能量。
根据式（1）计算出平衡后体系的相互作用能，结果见表2。由表2可以看出，褐煤煤尘体系（简称WD-DS）的相互作用能的绝对值比矽尘体系（简称SiO2-DS）的大，说明相同条件下，褐煤分子和活性剂水溶液的相互作用强度更强，褐煤煤尘表面吸附水的能力更强。
表 2 水−CTAB−煤尘/矽尘的相互作用能Table 2. Interaction energies of water-CTAB-coal dust/silicon dust体系 能量/（kcal·mol−1） 总能量 水和CTAB 粉尘 相互作用能 水−CTAB−煤尘 −2 739.640 2 −4 053.856 8 2 107.936 5 −793.720 4 水−CTAB−矽尘 −24 706.914 0 −4 065.075 7 −20 451.980 0 −135.598 0 注：1 kcal=4.184 kJ。 2.2 径向分布函数（RDF）
径向分布函数（RDF）被定义为特定原子在r距离内的密度与容重的比率。换句话说，特定原子的密度随与参考分子距离的变化随体积密度的变化表示RDF。所以，它可以用来演示给定分子周围的密度分布，数学表达式如下：
$$ g\left( r \right) = {{{\mathrm{d}}N} / {( {4\rho {\pi ^2}} )}} $$ (2) 式中：$ {\mathrm{d}}N $为与中心的距离为$ r $到$ r + {\mathrm{d}}r $之间的原子数目，个；$ \pi $为圆周率，取3.141 592 6；$ \rho $为该原子的平均数目密度。
根据分子动力学模拟的结果，分析得到的CTAB头基和H2O的径向分布曲线如图4所示，CTAB尾基和H2O的径向分布曲线如图5所示。
![]()
图 4 CTAB头基和H2O的径向分布曲线Figure 4. Radial distribution curves of CTAB head groups and H2O![]()
图 5 CTAB尾基和H2O的径向分布曲线Figure 5. Radial distribution curves of CTAB tail groups and H2O由图4可以看出，模拟平衡后，CTAB头基和水分子在WD-DS和SiO2-DS中的0.619 nm和0.597 nm处出现峰值，且SiO2-DS的峰值要比WD-DS的高，0.7 nm两曲线出现交点；WD-DS在0.885 nm处又出现第二峰，交点之后WD-DS的径向分布曲线一直在SiO2-DS的上方，两者缓慢下降。根据峰的位置判断3处峰是因为CTAB头基和水之间的静电力作用而形成的，虽然SiO2-DS出现的峰值要高，但WD-DS总体有2个峰值，并且交点以后WD-DS的径向分布曲线高于SiO2-DS的，说明WD-DS中头基和水的相互作用更强。
由图5可以看出，0.7 nm之前SiO2-DS的径向分布曲线略微高于WD-DS的曲线，0.7 nm之后相反。两曲线在1.4、1.6 nm处出现强峰，显然WD-DS的CTAB尾基和水分子的相互作用强于SiO2-DS。
2.3 均方位移和自扩散系数
均方位移MSD用来描述时间变化下目标粒子空间位置与初始位置的偏差程度。通过计算水分子的MSD和扩散系数D，可以进一步阐明表面活性剂加入对水分子运动的影响[17]。可以将MSD对时间的曲线拟合成直线$ y = ax + b $，然后利用以下公式计算扩散系数。
$$ {\mathrm{MSD}} = \frac{1}{N}\sum\limits_{{{i}} = 1}^N {{{{{[}}{{{r}}_i}\left( t \right){{ - }}{{{r}}_i}\left( 0 \right){{]}}}^{{2}}}} $$ (3) $$ D = \frac{a }{6} $$ (4) 式中：$ N $为粒子数；$ {r_i}\left( t \right) $为$ i $粒子在$ t $时刻的位置；$ {r_i}\left( 0 \right) $为$ i $粒子的初始位置；$ a $为拟合之后直线的斜率。
水分子的MSD拟合直线如图6所示。根据图6均方位移MSD拟合直线的斜率计算可得自扩散系数DWD=0.004 5 nm2/ps和$D_{{\mathrm{SiO}}_2 }$=0.004 2 nm2/ps，从D的取值来看，CTAB增强了褐煤煤尘和矽尘的润湿性，使得表面和水的相互作用增强，加速了水分子的运动，进而提高了表面的润湿性，但对褐煤的效果更好。
2.4 数密度
界面区域是相变现象发生的主要场所，是其本质就是从一种相态密度快速变化至另一种相态密度的区域。根据10−90准则[18-19]，界面厚度被定义为密度为其体积相密度的10%~90%的2个位置之间的距离。仅对z轴方向水分子的相对数密度进行分析，WD-DS中的水分子相对数密度如图7所示，SiO2-DS中的水分子相对数密度如图8所示。
由图7和图8可知：WD-DS中水的最大相对数密度为6.49，并依据10−90准则得出最大相对数密度的10%和90%位点之间的距离约为1.71、3.30 nm，由此可以算出活性剂水溶液和煤的界面厚度约为1.59 nm；同理，算得活性剂水溶液和SiO2的界面厚度约为0.971 nm。由此可得WD-DS的固液界面层厚度大于SiO2-DS。
3. 静电势
静电势ESP是指从无穷远处移动单位正电荷至该点时所需做的功。分析分子表面静电势可预测其容易发生吸附的位点。各个分子或晶胞静电势图如图9所示。图中：红色区域表示正电荷，该区域的电子收缩表明它很难失去电子，且可以相对容易地接受电子形成氢键；蓝色区域表示负电荷，该区域较容易提供电子形成氢键；白色区域表示电位接近于0的区域，具有相对稳定结构[20]。各个分子或晶胞的静电势对比表见表3。
![]()
图 9 各个分子或晶胞静电势图Figure 9. Diagrams of the electrostatic potential of individual molecules or crystals表 3 各个分子或晶胞的静电势对比表Table 3. Comparison of the electrostatic potential of individual molecules or crystals分子或晶胞种类 静电势/（Ha·e−1） 正 负 尾链 H2O 0.091 56 −0.065 44 — CTAB 0.203 40 — 0.021 89 WD 0.124 20 −0.070 96 — SiO2 0.568 80 −0.020 57 — CTAB头基的最大正静电势0.203 40 Ha/e，大于水分子的最大正静电势0.091 56 Ha/e，可以得出 CTAB的头基是亲水的。尾链的正静电势为0.021 89 Ha/e，小于水分子的最大负静电势0.065 44 Ha/e，因此，CTAB尾链是疏水的。
平衡后，CTAB头基一部分吸附于WD表面的负静电势位置，另一部分与水进行吸附。当头基和WD表面负静电势吸附时，头基会占据原先水分子吸附的界面位置，尾基存在疏水性，朝向水中。吸附前该位置与水分子的电势差为0.162 52 Ha/e，当头基吸附后，该位置与水的电势差变为0.092 85 Ha/e，WD的亲水性减弱。CTAB尾基一部分吸附于WD表面的负静电势位置时，吸附后电势差变为0.274 36 Ha/e, WD的亲水性增强。CTAB头基尾基吸附电势差变化是增大的，对WD分子亲水性是增强的，润湿性有所提升。同理，CTAB对SiO2的润湿作用有所提升。相比之下SiO2加入活性剂之后的电势差增幅更大，润湿效果应该更好。
通过对分子间相互作用能对比可知，WD与水分子之间的作用能更大，润湿性更好，可能主要是因为WD分子中存在大量官能团，比如醚、醇、羧酸、酚、酮等，官能团的存在一定程度上会增加氢键的形成，进而增强分子间的相互作用力，这可能是SiO2和活性剂水溶液的相互作用较弱的原因，导致了SiO2润湿性差。
4. 结 语
1）基于分子动力学模拟的相互作用能和均方位移定量分析表明：在同等质量的WD和SiO2、相同表面活性剂CTAB水溶液条件下，WD和活性剂水溶液的相互作用能为−793.720 4 kcal/mol，SiO2的为−135.598 0 kcal/mol，即CTAB对WD的润湿性较好；结合水分子自扩散系数D验证了CTAB对WD的润湿性较好，其中DWD=0.004 5 nm2/ps 和$D_{{\mathrm{SiO}}_2 } $=0.004 2 nm2/ps，D值越大，水的扩散效果越好。
2）基于分子动力学模拟的径向分布函数和相对数密度的定性分析表明：相较于尾基，CTAB头基和水峰值高、作用强，矽尘和水作用相比于WD弱；由相对数密度所得的固液界面厚度可知，WD和水的界面厚度1.59 nm大于SiO2和水的0.971 nm，所以WD中沿z轴方向润湿性效果好。
3）静电势分析结果表明：表面活性剂CTAB头基是亲水头基，尾基是疏水的。活性剂吸附使得表面发生电势差的变化是影响表面润湿性的因素，CTAB对SiO2和WD 都有润湿效果，单单考虑电势差的变化程度，SiO2的电势差变化程度要大，但是考虑到WD官能团的作用，WD和CTAB水溶液的相互作用能大，CTAB对WD润湿性更好。
-
![]()
图 5 不同注胶压力、时间下凝胶扩散形态变化云图
Figure 5. The cloud diagrams of gel diffusion morphology changes under different injection pressure and time conditions
![]()
图 7 凝胶扩散面积占比随时间变化曲线
Figure 7. The curves of the proportion of gel diffusion area with time
表 1 数值模拟方案
Table 1 Scheme of numerical simulation
序号 注胶压力/MPa 时间/min 渗透率/m2 孔隙率 1 2 180 5×10−15 0.4 2 4 180 5×10−15 0.4 3 6 180 5×10−15 0.4 4 8 180 5×10−15 0.4 5 10 180 5×10−15 0.4 6 12 180 5×10−15 0.4 
下载: 导出CSV
-
[1] 秦波涛,仲晓星,王德明,等. 煤自燃过程特性及防治技术研究进展[J]. 煤炭科学技术,2021,49(1):66−99. QIN Botao, ZHONG Xiaoxing, WANG Deming, et al. Research progress of coal spontaneous combustion process characteristics and prevention technology[J]. Coal Science and Technology, 2021, 49(1): 66−99.
[2] 张创业,王晓川,吕有厂,等. 下保护层开采扰动卸压效应及瓦斯与煤自燃协同治理技术[J/OL]. 矿业安全与环保,2024: 1-8. [2024-04-26]. http://kns.cnki.net/kcms/detail/50.1062.td.20230829.0847.001.html. ZHANG Chuangye, WANG Xiaochuan, LYU Youchang, et al. Collaborative technology of mining disturbance pressure relief effect and gas-spontaneous combustion collaborative control technology[J/OL]. Mining Safety & Environmental Protection, 2024: 1-8. [2024-04-26]. http://kns.cnki.net/kcms/detail/50.1062.td.20230829.0847.001.html.
[3] 赵亚明,杨伟,陈洋. 冲击地压矿井易自燃厚煤层回撤工作面自燃防治技术[J]. 煤矿安全,2022,53(4):81−86. ZHAO Yaming, YANG Wei, CHEN Yang. Spontaneous combustion prevention technology of withdrawal working face of spontaneous combustion thick coal seam in rock burst mine[J]. Safety in Coal Mines, 2022, 53(4): 81−86.
[4] 仲晓星,杨正杰,窦国兰,等. 地下煤火热交换系统相似准则研究[J]. 中国矿业大学学报,2015,44(6):965−970. ZHONG Xiaoxing, YANG Zhengjie, DOU Guolan, et al. Study on similarity criterion for heat exchange system of underground coal fires[J]. Journal of China University of Mining & Technology, 2015, 44(6): 965−970.
[5] 付伟,宋海洲,陈军,等. 新型矿用防灭火凝胶泡沫制备及阻化性能研究[J]. 煤炭工程,2023,55(9):127−133. FU Wei, SONG Haizhou, CHEN Jun, et al. Study on preparation and chemical inhibition properties of a new gel foam for mine fire extinguishing[J]. Coal Engineering, 2023, 55(9): 127−133.
[6] 陈鹏燕,周春山,程熙宇. CMC/AlCit/GDL防灭火凝胶的流变性能及本构方程研究[J]. 应用化工,2023,52(4):1010−1014. CHEN Pengyan, ZHOU Chunshan, CHENG Xiyu. Study on rheological ability and constitutive equation of mineral CMC/AlCit/GDL gel[J]. Application Chemical Indystry, 2023, 52(4): 1010−1014.
[7] 李洋,戴广龙,杨苗苗. 高保水矿用凝胶阻化材料的制备和性能研究[J]. 煤炭技术,2023,42(3):173−177. LI Yang, DAI Guanglong, YANG Miaomiao. Preparation and performance study of high water retention mineral hydrogel inhibitors[J]. Coal Technology, 2023, 42(3): 173−177.
[8] 段西凯,周睿,何勇军. 矿用超吸水防灭火凝胶制备及性能研究[J]. 煤矿安全,2023,54(2):71−76. DUAN Xikai, ZHOU Rui, HE Yongjun. Study on preparation and performance of mine superabsorbent fire prevention and extinguishing gel[J]. Safety in Coal Mines, 2023, 54(2): 71−76.
[9] 石政锋. 大采高煤层综采工作面采空区自燃“三带”立体分布规律研究[J]. 矿业安全与环保,2019,46(2):47−50. SHI Zhengfeng. Study on stereoscopic distribution law of spontaneous combustion“three-zone”in goaf of fully mechanized working face with large mining height[J]. Mining Safety & Environmental Protection, 2019, 46(2): 47−50.
[10] 王正帅. 复杂通风综采面采空区“三带”分布及防灭火参数研究[J]. 矿业研究与开发,2023,43(8):73−77. WANG Zhengshuai. Study on the distribution of three zones and fire prevention and extinguishing parameters in goaf of complex ventilation fully mechanized mining face[J]. Mining Research and Development, 2023, 43(8): 73−77.
[11] 张嘉麟,董子文. 基于COMSOL动态模型的采空区连续注氮模拟研究[J]. 矿业安全与环保,2019,46(2):15−19. ZHANG Jialin, DONG Ziwen. Simulation study on continuous nitrogen injection in goaf based on COMSOL dynamic model[J]. Mining Safety & Environmental Protection, 2019, 46(2): 15−19.
[12] 李尚国. 采空区氧化中心数值模拟[J]. 煤炭技术,2019,38(4):109−112. LI Shangguo. Numerical simulation of oxidation center in goaf[J]. Coal Technology, 2019, 38(4): 109−112.
[13] 齐庆杰,祁云,周新华. 供风量对采空区自燃动态影响及防灭火技术[J]. 中国安全科学学报,2019,29(4):120−126. QI Qingjie, QI Yun, ZHOU Xinhua. Dynamic effects of air supply volume on spontaneous combustion in goafs and fire prevention and extinguishing technology[J]. China Safety Science Journal, 2019, 29(4): 120−126.
[14] 李宗翔,刘宇,王政,等. 九道岭矿采空区注CO2防灭火技术数值模拟研究[J]. 煤炭科学技术,2018,46(9):153−157. LI Zongxiang, LIU Yu, WANG Zheng, et al. Study on numerical simulation of CO2 injection fire prevention and extinguishing technique in goaf of Jiudaoling mine[J]. Coal Science and Technology, 2018, 46(9): 153−157.
[15] 刘玉良,彭宝山,蒙文军,等. 浅埋深煤层采空区注氮防灭火技术[J]. 煤矿安全,2017,48(5):73−76. LIU Yuliang, PENG Baoshan, MENG Wenjun, et al. Goaf nitrogen injection fire-fighting technology of shallow buried depth coal seam[J]. Safety in Coal Mines, 2017, 48(5): 73−76.
[16] 邓军,李贝,王凯,等. 我国煤火灾害防治技术研究现状及展望[J]. 煤炭科学技术,2016,44(10):1−7. DENG Jun, LI Bei, WANG Kai, et al. Research status and outlook on prevention and control technology of coal fire disaster in China[J]. Coal Science and Technology, 2016, 44(10): 1−7.
[17] 余明高,王亮,李海涛,等. 我国煤矿防灭火材料的研究现状及发展趋势[J]. 矿业安全与环保,2022,49(4):22−36. YU Minggao, WANG Liang, LI Haitao et al. Research status and development trend of fire-extinguishing materials in Chinese coal mines[J]. Mining Safety & Environmental Protection, 2022, 49(4): 22−36.
[18] 谢军,张心静,徐书荣,等. 泡沫防灭火材料研究现状及发展趋势[J]. 煤炭技术,2023,42(1):129−132. XIE Jun, ZHANG Xinjing, XU Shurong, et al. Research status and development trend of foam fire prevention materials[J]. Coal Technology, 2023, 42(1): 129−132.
[19] 韩福志,邵和,司俊鸿,等. 矿用堵漏无机触变防灭火凝胶的流变特性研究[J]. 矿业安全与环保,2022,49(1):1−7. HAN Fuzhi, SHAO He, SI Junhong, et al. Study on rheological characteristics of inorganic thixotropic fire extinguishing gel for blocking air leakage of coal mine[J]. Mining Safety & Environmental Protection, 2022, 49(1): 1−7.
-
期刊类型引用(1)
1. 刘啸. 吸附条件下含裂隙煤体渗透率演化规律数值模拟分析. 煤矿安全. 2025(03): 12-20 . 
本站查看
其他类型引用(0)
计量
- 文章访问数: 31
- HTML全文浏览量: 3
- PDF下载量: 7
- 被引次数: 1
