Modeling and mechanism of propagation of underground coal smoldering fires with thermal buoyancy driven
-
摘要:
地下煤火蔓延模型与机制研究对我国煤炭资源绿色安全开采具有重要意义。目前地下煤火阴燃正向蔓延模型主要存在问题：简化的煤燃烧单步反应忽略了水分蒸发和热解吸热反应对蔓延的重要影响作用；缺乏实验数据验证模型的有效性。为此，构建了包含水分蒸发、煤热解和碳氧化三步反应体系的地下煤火阴燃正向蔓延数学模型；采用COMSOL Multiphysics有限元软件数值计算了该理论模型；开展了不同地裂缝渗透性条件下热浮力驱动地下煤火阴燃正向蔓延实验，进行对比实验和数值计算。结果表明：构建的模型揭示了地下煤火阴燃蔓延供氧控制机制，不仅能预测地下煤火阴燃正向蔓延高温区域的温度和蔓延速率，而且还能合理地预测阴燃多步反应速率以及氧气、煤、碳和灰分质量分数的时空演化。
Abstract:Modeling and understanding of underground coal fire is of significance for green-and-safe mining of coal resources in China. Two issues remain unsolved in forward-propagation models of underground coal smoldering fires: one-step coal oxidation has been too simple so that endothermic reactions such as water evaporation and pyrolysis that have important influence on smoldering propagation of underground coal fires have been not considered; and models have not been validated via experimental data. In this work, a three-step chemical reaction scheme including water evaporation, coal pyrolysis, and char oxidation was employed and mathematic model in terms of forward propagation of underground coal smoldering fires was established. The theoretical model was numerically calculated using COMSOL Multiphysics finite element software. The theoretical model was numerically calculated using COMSOL Multiphysics finite element software, and the forward smoldering spread experiment of underground coal fire driven by thermal buoyancy under different ground crack permeability conditions was carried out. Comparison experiment and numerical calculation were carried out. Results show that the proposed model not only owns good capability to predict the peak temperature and velocity of smoldering propagation of underground coal fire, but also can inversely plot multiple step chemical reactions as well as spatiotemporal evolutions of oxygen and solid species mass fractions. In addition, the controlling mechanism of smoldering propagation of underground coal fires was further revealed.
-
煤矿粉尘是煤矿开采、运输及储存过程中伴随煤和岩石破碎而产生的混合性粉尘，主要含有煤尘、岩尘以及少量其他物质[1]。煤矿粉尘可分为呼吸性粉尘和非呼吸性粉尘。一般说来，小于5 μm的尘粒能到达并沉积于肺泡中，故称为呼吸性粉尘，是引起尘肺的主要尘粒[2-3]。在《国家职业病防治规划（2016—2020年）》《安全生产“十三五”规划》及《“健康中国2030”规划纲要》中，粉尘职业危害防治已被列入国家战略与重大需求[4]。口罩等防护用具作为煤矿环境劳动者个体防护的最后一道防线，其滤料的研发具有重要的研究意义。
防护用器具对气溶胶的过滤理论以纤维层空气过滤理论为基础，经过十几年的发展，纤维过滤理论已逐渐趋于成熟[5]。纤维过滤过程中，捕集颗粒物的过滤作用有：扩散效应、拦截效应、惯性效应、重力效应、和静电吸附[6]。与传统过滤介质相比，静电纺丝制备的纳米纤维膜具有优异的过滤性能、纤维直径小、孔隙率高、体积比显著的表面积压降[7-8]。目前，许多种高分子量聚合物，包括聚丙烯腈（PAN）、聚苯乙烯（PS）、聚乙烯醇（PVA）和聚丙烯（PP）等[9]，已通过静电纺丝加工成纳米纤维，用于空气过滤器。聚对苯二甲酸乙二酯（PET）是世界上规模最大、应用最广泛的合成纤维，是由对苯二甲酸（TA）和乙二醇（EG）脱水酯化得到的半结晶热塑性饱和聚酯[10-12]。虽然静电纺丝技术已有多年历史，但其在聚酯方面的应用还比较少。李静等[13]运用PET与含氟基团改性的SiO2共混进行静电纺丝，构建具有微纳复合结构的低表面能表面，具有优良的超疏水性；ZANDER等[14]使用塑料瓶回收并重新聚合PET，但所制备电纺过滤膜对500 nm以下的颗粒去除效果微乎其微。因此，亟须制备一种对粉尘同时具备高过滤效率和低阻力的滤料。此外，作为一种个体防护材料，呼吸过程中细菌的繁殖对人类健康有很大的影响，当前很少将抗菌性静电纺丝纳米纤维膜作为煤矿井下个体防护材料。为此，以静电纺丝技术为基础，将聚对苯二甲酸乙二醇酯和银纳米颗粒共混，以制备煤矿粉尘个体防护滤料（PET/Ag）；利用扫描电子显微镜（SEM）对纤维形貌、分布以及表面孔隙进行表征，利用X射线光电子能谱（XPS）对纤维表面的化学性质及组成进行分析；实验研究了PET/Ag纳米纤维膜作为个体防护材料的可行性，测试了其纤维膜的防水透湿性能、过滤性能以及抗菌性能。
1. 制备方法及实验测试
1）试剂。选用的试剂为：聚对苯二甲酸乙二醇酯（PET，99%）、二氯甲烷（DCM，99.5%，化学纯）、三氟乙酸（TFA，分析纯）、AgNO3溶液（0.1 mol/L）、N,N-二甲基甲酰胺（DMF，99.5%，分析纯）、聚乙烯吡咯烷酮（PVP，分析纯）。
2）纳米纤维膜的制备。将PET颗粒溶解于三氟乙酸和二氯甲烷的混合溶剂中（V（TFA）∶V（DCM）=4∶1），在室温下搅拌溶解成均一的PET纺丝溶液。通过液相还原法[15]制备银纳米粒子（AgNPs）。在pH=8的条件下，DMF作为还原剂，聚乙烯吡咯烷酮（PVP）作为保护剂，在DMF中煮沸AgNO3制备银纳米颗粒。将纳米银在去离子水中超声溶解1 h，与PET颗粒共混制备纺丝液，纳米银的含量为0.4%，然后在室温下磁力搅拌10 h，获得质量分数分别为8%、10%、12%的均一纺丝溶液。使用HZ-12静电纺丝机进行纳米纤维膜的制备。将制备的纺丝溶液置于5 mL注射器中，并使用注射泵以1 mL/h的速率推出。针头与20 kV高压电源相连。在静电纺丝过程中，针头与收集器之间的距离保持在15 cm，滚筒转速为100 r/min，温度保持在（25±2）℃，相对湿度为（35±5）%。控制纺丝时间为1、2、3 h，得到不同厚度的PET/Ag纳米纤维膜。
3）性能表征。通过扫描电子显微镜（SEM）（JSM-6360LV）观察纳米纤维膜的形态。将待测样品用导电胶带粘贴在样品台上，经喷金处理后，放入扫描电镜样品室，在10 kV加速电压下，观察纤维膜的微观形貌。在SEM图像中选择100根纳米纤维，使用图像处理软件（Image-Pro Plus）绘制纤维直径分布直方图，并计算纳米纤维平均直径；通过射线光电子能谱（XPS）（ESCALAB 250）测定电子的结合能，对纤维膜样品表面的化学性质及组成进行分析，了解其元素组成。通过水接触角测试仪（Dataphysics DCAT21）在室温和45%湿度条件下进行测试，水滴在纤维膜表面稳定后测定其接触角大小。
4）透湿性能测试。使用WPT-301水蒸气透过率测试仪测量材料的润湿性。在38 ℃、20% RH的条件下对每个样品进行3次重复测量。水蒸气透过率（WVP）由式（1）计算：
$$ {\mathrm{WVP}} = \frac{{\Delta w}}{t} \cdot \frac{x}{{A \cdot \Delta p}} $$ (1) 式中：Δw为质量变化，g；t为时间；x为膜的厚度，m；A为膜样品面积，m2；Δp 为38 ℃下通过膜的蒸汽压差，kPa。
5）过滤性能测试。使用LZC-K1型滤料综合测试仪对纳米纤维膜进行过滤性能测试，该仪器在计算纤维膜的过滤效率的同时还可通过传感器得出气流时的压阻值。过滤效率和品质因子由式（2）和式（3）计算。所得数据为3次测量的平均值和标准偏差。
$$ \eta = \frac{{N_1 - N_2}}{{N_1}} \times 100\% $$ (2) $$ {\mathrm{QF}} = \frac{{ - \ln (1 - \eta )}}{{\Delta p_1}} $$ (3) 式中：η为过滤效率；N1、N2分别为在通过纳米纤维膜之前和之后的气溶胶颗粒的数量；QF为品质因子；Δp1为相应压阻值。
6）抑菌性能测试。选择大肠杆菌（ATCC 25922）和金黄色葡萄球菌（ATCC 25923）作为实验微生物。首先，将制备的营养肉汁琼脂倒在无菌培养皿上，待其凝固，将大肠杆菌和金黄色葡萄球菌培养液（各0.5 mL）均匀涂在琼脂平板上；然后将纤维膜均匀切割成直径为1 cm的圆形，将掺杂AgNPs的样品和对照样品（PET膜）轻轻地放置在固化琼脂板上；最后将培养板在37 ℃下培养18 h，通过1个清晰的细菌抑菌圈直径来鉴定和评估抗菌活性，所得抑菌圈直径由图像处理软件进行量化。实验重复进行3次。
2. 实验结果讨论
2.1 纳米纤维膜的形态与结构
不同质量分数纺丝液所纺PET/Ag纳米纤维膜的SEM电镜图片如图1所示；通过射线光电子能谱（XPS）测定电子的结合能对PET和PET/Ag纳米纤维膜表面的元素组成进行分析得出的结果如图2所示。
![]()
图 1 不同质量分数的PET/Ag纤维膜SEM图像及纤维直径分布图Figure 1. SEM images and fiber diameter distribution of PET/Ag membranes with different mass fractions![]()
图 2 PET和PET/Ag纳米纤维膜的XPS图谱Figure 2. XPS images of PET and PET/Ag nanofibrous membranes由图1可知：通过调节纺丝液的质量分数，可以精确地控制纳米纤维膜的纤维径，而纤维直径的形态几乎不受影响，纤维间自由搭接，没有出现粘连现象，纤维间空隙较小，孔隙数目多，结构明显，单个纤维粗细均匀，没有出现串珠结构，纤维间无规则交织排列在一起，形成错综排列的网状结构[16-17]。纤维直径分布图显示了随机选取的100根纤维的直径大小分布，在1 mL/h的纺丝速率下，纳米纤维平均直径从质量分数8%时的298 nm增加到质量分数10%时的424 nm，到12%时的836 nm。纺丝液质量分数为8%的纤维膜的直径较小，主要集中在250~300 nm，纤维间空隙较小；纺丝液质量分数为10%的纤维膜的直径主要集中在400~450 nm；纺丝液质量分数为12%纤维膜的直径较大，主要集中在700~800 nm，纤维间空隙也相对较大。
由图2可以看出：PET纳米纤维膜表面的主要由C、O元素的峰非常明显，峰值分别对应284.6、531.5 eV；而PET/Ag的测试结果在284.6、367.8、531.5 eV处观察到3个峰，分别对应于C、Ag和O元素的峰， 367.8 eV处峰的出现说明Ag+还原成AgNPs并成功负载到PET膜上。相比之下，没有AgNPs负载的PET膜中没有出现Ag 3d峰。
2.2 纳米纤维膜的疏水性表征
质量分数为8%、10%、12%的PET/Ag纺丝液纺丝2 h后纤维膜表面水接触角变化情况如图3所示。
![]()
图 3 不同质量分数的PET/Ag纤维膜表面水接触角Figure 3. Water contact angle on the surface of PET/Ag fiber membranes with different mass fractions由图3可知：对不同质量分数的PET/Ag膜，水接触角可以在较长时间内保持不变，表现为较好的疏水性；质量分数为8%、10%、12%PET/Ag纤维膜表面水接触角分别为133.88°、131.51°、132.04°。水接触角测试结果表明PET/Ag 膜具有较强的疏水性。
2.3 透湿性分析
随着电纺时间的增加，PET/Ag纤维无规则交织排布，纤维的密集度不断增加，纤维间的空气流通通道愈加曲折，孔隙则不断减小[18]。不同质量分数的纺丝液在不同纺丝时间下纳米纤维膜的水蒸气透过量如图4所示。
![]()
图 4 不同质量分数的纺丝液在不同纺丝时间下纳米纤维膜的水蒸气透过量Figure 4. The water vapor permeability of nanofiber membranes with different mass fractions of spinning solution under different spinning time纺丝液质量分数影响纤维直径大小，纺丝时间影响纳米纤维膜的厚度，而纤维直径和厚度都影响了纳米纤维膜的水蒸气流通路径。随着纺丝液质量分数的增加，PET/Ag纳米纤维膜的水蒸气透过量呈增加趋势，这可能是由于纺丝液质量分数影响了纤维直径的大小，较大的纤维直径有利于形成较大的孔隙，使水蒸气流通更加通畅；随着纺丝时间的增加，PET/Ag纳米纤维膜的水蒸气透过量呈减小趋势，这可能是由于纺丝时间影响了纳米纤维膜的厚度，较大的厚度会使纤维间孔隙通路变长，水蒸气流通路径更加曲折，阻力更大[19]。从PET/Ag纳米纤维膜的透湿性能来看，纺丝液质量分数为12%、纺丝1 h的纳米纤维膜透湿性最好，高达
9727.55 g/（m2·d），纺丝液浓度为8%、纺丝3 h的纳米纤维膜透湿性最差，为2069.82 g/（m2·d）。2.4 过滤性能测试分析
煤矿中大多数煤尘颗粒的直径在1~10 μm之间，LZC-K1型滤料综合测试仪能够产生粒径范围为0.3~10 μm的NaCl气溶胶颗粒，具有更强的穿透性[20-21]。因此，以NaCl颗粒作为测试介质，在85 L/min空气流量下更严苛地反映滤料的过滤效率。实验分别测试的不同纺丝时间下的PET/Ag纳米纤维膜对空气中颗粒物的过滤效果如图5所示；不同纺丝时间下的PET/Ag纳米纤维膜在过滤过程中的压阻值如图6所示。
![]()
图 5 PET/Ag纳米纤维膜在不同纺丝时间下的过滤效率Figure 5. Filtration efficiency of PET/Ag nanofibrous membranes at different spinning time![]()
图 6 PET/Ag纳米纤维膜在不同纺丝时间下的压阻Figure 6. Pressure drop of PET/Ag nanofibrous membranes at different spinning time由图5可以看出：纺丝液质量分数越高，PET/Ag纳米纤维膜的过滤效率越好，且随着纺丝时间的延长，不同质量分数的纺丝液所纺纳米纤维膜的颗粒物过滤效率呈增加趋势；其中，纺丝液质量分数为10%和12%的PET/Ag纳米纤维膜在纺丝3 h后的过滤效率都可达99.99%，纺丝液质量分数为8%的PET/Ag纳米纤维膜过滤效率为98.75%。结果表明：在一定时间内，增加纺丝时间，纳米纤维膜不断堆积，厚度随之增加，纤维之间气流通道更加曲折，纤维膜对颗粒物的过滤效果有所提升。因此，通过控制静电纺丝时间可以得到不同厚度，对颗粒物有不同过滤效果的材料。
由图6可以看出：纺丝液质量分数越高，PET/Ag纳米纤维膜的压阻值越高，且随着纺丝时间的延长，不同纺丝液浓度所纺纳米纤维膜的压阻值呈增加趋势；当纺丝时间为1 h时，8%、10%、12% PET纤维膜的压阻分别为78.2、85.1、90.2 Pa；当纺丝时间增加到3 h时，3种PET/Ag纳米纤维膜的压阻值分别为159.3、175.6、200.1 Pa。
结合不同PET/Ag纳米纤维膜的过滤效率，通过公式计算相对应的品质因子。通过计算可知：纺丝液质量分数为12%，纺丝时间为3 h时，PET/Ag纳米纤维膜的品质因子最低，为0.013 Pa−1；纺丝液质量分数为10%，纺丝时间为2 h时，PET/Ag纳米纤维膜的品质因子最高，为0.041 Pa−1，该条件下纳米纤维膜的水蒸气透过量为同时7 101.37 g/（m2·d），透湿性较好。
纺丝液质量分数为10%，纺丝时间为2 h时PET/Ag纳米纤维膜在15~85 L/min不同空气流量下的过滤效率、压阻和品质因子见表1。
表 1 不同空气流量下PET/Ag纳米纤维膜的过滤效率、压阻和品质因子Table 1. Filtration efficiency, pressure drop and quality factor of PET/Ag nanofiber membrane under different air flow conditions空气流量/（L·min−1） 过滤效率/% 压阻/Pa 品质因数/Pa−1 15 99.99 32.5 0.084 30 99.42 49.8 0.075 45 98.19 60.6 0.069 60 96.78 89.4 0.056 85 95.14 90.7 0.041 由表1可以看出：随着空气流量增加，流速增大，纳米纤维对颗粒物的拦截、扩散作用减弱，纳米纤维膜的过滤效率逐渐降低；并且流速增强时压力差随之增大，压阻随之升高，品质因子也呈降低趋势，综合性能随流量的增加表现得越来越差，即流量越大，口罩的防护性能和舒适性均随流量增加而变差。
2.5 抑菌性能分析
利用琼脂扩散法评价所制备滤料对微生物的作用效果如图7所示。
![]()
图 7 PET和PET/Ag纳米纤维膜对大肠杆菌和金黄色葡萄球菌的抑菌区Figure 7. Bacteriostatic zones of PET and PET/Ag nanofiber membrane against E.coli and S.aureus由图7（a）可以看出，在37 ℃下培养18 h后，PET纳米纤维膜对大肠杆菌和金黄色葡萄球菌没有显示出任何抗菌活性，说明PET膜不具备抗菌性能。由图7（b）可以看出共混静电纺丝法制备的PET/Ag纳米纤维膜对2种细菌都显示出明显的抗菌活性，其中，对大肠杆菌的抑菌区直径最大，为18.12 mm，对金黄色葡萄球菌的抑菌区为13.25 mm；这可能是由于金黄色葡萄球菌的细胞壁较厚，蛋白质含量较少，因此其抑菌区直径要比大肠杆菌小。结果表明：PET/Ag纳米纤维膜对大肠杆菌和金黄色葡萄球菌细菌都具备明显的抗菌活性，作为空气过滤材料能够有效抑制细菌繁殖，在个人防护领域具有广阔的应用前景。
3. 结 语
1）SEM图显示，PET/Ag电纺纤维较细，纤维形貌良好，纤维之间无规则排列在一起，没有出现粘连现象，纤维膜孔隙结构明显，单个纤维粗细均匀，没有出现串珠结构。通过调节纺丝液的质量分数，可以精确地控制纳米纤维膜的纤维直径，而纤维直径的形态几乎不受影响，纤维间自由搭接无规则交织排列在一起，形成错综排列的网状结构。XPS表明Ag+还原成AgNPs并成功负载到PET膜上。
2）水接触角测量表明PET/Ag纳米纤维膜表现为较好的疏水性。从PET/Ag纳米纤维膜的透湿性能来看，纺丝液质量分数和纺丝时间都影响了纳米纤维膜的水蒸气流通路径。纺丝液质量分数为12%、纺丝1 h的纳米纤维膜透湿性最好，高达9 727.55 g/（m2·d）。
3）增加纺丝时间，纳米纤维膜不断堆积，厚度随之增加，纤维之间气流通道更加曲折，纤维膜对颗粒物的过滤效果有所提升。因此，通过控制静电纺丝时间可以得到不同厚度，对颗粒物有不同过滤效果的材料。
4）随空气流量从15 L/min增加到85 L/min， PET/Ag纳米纤维对颗粒物的拦截、扩散作用减弱，纳米纤维膜的过滤效率逐渐降低。并且流速增强时压力差增大，压阻随之升高，综合性能随流量的增加表现得越来越差。
5）PET/Ag纳米纤维膜对大肠杆菌和金黄色葡萄球菌细菌都具备明显的抗菌活性，作为空气过滤材料能够有效抑制细菌繁殖，在个人防护领域具有广阔的应用前景。
-
![]()
图 1 热浮力驱动地裂缝气流作用下地下煤火阴燃实验示意图
Figure 1. Schematic diagram of underground coal smoldering fire experiment
![]()
图 2 热浮力驱动作用下地下煤火阴燃正向蔓延模型验证
Figure 2. Validation of the proposed model against with experimental data
![]()
图 3 不同泡沫孔隙率条件下供氧速率
Figure 3. Oxygen supply rate under different foam porosity conditions
![]()
图 4 泡沫孔隙率为0.9522条件下t=6 h（第1列）、18 h（第2列）和30 h（第3列）的温度场及氧气体积分数、煤、碳、灰分质量分数分布
Figure 4. Cloud diagrams of temperature field, oxygen, coal, carbon and ash mass fraction distribution (from top to bottom) of t=6 h (the first column), 18 h (the second column) and 30 h (the third column) with condition of foam porosity of 0.952 2
![]()
图 5 泡沫孔隙率为0.9522条件下中心轴线3个不同位置的反应速率
Figure 5. Chemical reaction rates of three different positions at the central axis with condition of foam porosity of 0.952 2
![]()
图 6 不同泡沫孔隙率条件下中心轴线y=8 cm位置水分蒸发（rw）、煤热解（rp）和=碳氧化（ro）反应速率
Figure 6. Reaction rates of water evaporation(rw), coal pyrolysis(rp) and carbon oxidation(ro) at y=8 cm with different porosity of foam
表 1 煤阴燃反应动力学参数、化学计量系数和反应热
Table 1 Kinetic parameters, stoichiometric coefficients and reaction heat of coal smoldering reaction
反应 指前因子
（lgAj）活化能
（Ej）/(kJ·mol−1)反应级数
（ni）化学计量
系数
（vj）反应热
（ΔH）/(MJ·kg−1)水分蒸发 8.120 0 67.800 0 2.3702 0.200 2.260 00 煤热解 10.4136 194.6619 9.2711 0.731 0.335 00 碳氧化 8.0743 166.2853 2.1104 0.209 28.35132 
下载: 导出CSV
-
[1] STRACHER G B. Coal fires burning around the world: a gobal catastrophe[J]. International Journal of Coal Geology, 2004, 59(1/2): 1−6.
[2] 张建民,管海晏,曹代勇,等. 中国地下煤火研究与治理 [M]. 北京:煤炭工业出版社,2008. [3] SONG Z, KUENZER C. Coal fires in China over the last decade: A comprehensive review[J]. International Journal of Coal Geology, 2014, 133: 72−99.
[4] 邓军,李贝,王凯,等. 我国煤火灾害防治技术研究现状及展望[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 Scienceand Technology, 2016, 44(10): 1−7.
[5] 梁运涛,侯贤军,罗海珠,等. 我国煤矿火灾防治现状及发展对策[J]. 煤炭科学技术,2016,44(6):1−6. LIANG Yuntao, HOU Xianjun, LUO Haizhu, et al. Development countermeasures and current situation of coal mine fire prevention & extinguishing in China[J]. Coal Science and Technology, 2016, 44(6): 1−6.
[6] 秦波涛,仲晓星,王德明,等. 煤自燃过程特性及防治技术研究进展[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.
[7] SHI B, SU H, LI J, et al. Clean power generation from the intractable natural coalfield fires: turn harm into benefit[J]. Scientific Reports, 2017, 7(1): 5302. doi: 10.1038/s41598-017-05622-4
[8] OHLEMILLER T J, LUCCA D A. An experimental comparison of forward and reverse smolder propagation in permeable fuel beds[J]. Combust Flame, 1983, 54(1): 131−147.
[9] REIN G. Smouldering combustion phenomena in science and technology[J]. International Review of Chemical Engineering, 2009(1): 3−18.
[10] 刘晓晨,武建军,刘明,等. 气象条件对地下煤火浅层温度的影响[J]. 煤炭学报,2011,36(8):1334−1341. LIU Xiaochen, WU Jianjun, LIU Ming, et al. Effect of local meteorological conditions on shallow subsurfacetemperature above underground coal fire[J]. Journal of China Coal Society, 2011, 36(8): 1334−1341.
[11] SONG Z, ZHU H, XU J, et al. Effects of atmospheric pressure fluctuations on hill-side coal fires and surface anomalies[J]. International Journal of Mining Science and Technology, 2015, 25(6): 1037−1044. doi: 10.1016/j.ijmst.2015.09.024
[12] SONG Z, ZHU H, TAN B, et al. Numerical study on effects of air leakages from abandoned galleries on hill-side coal fires[J]. Fire Safety Journal, 2014, 69: 99−110. doi: 10.1016/j.firesaf.2014.08.011
[13] 曾强,常心坦. 新疆煤田火区火风压模式研究及其应用[J]. 煤炭学报,2007,32(9):955−958. ZENG Qiang, CHANG Xintan. Study on the model of fire-heating airflow and its application to Xinjiang coal-field fires[J]. Journal of China Coal Society, 2007, 32(9): 955−958.
[14] 邓军,文虎,张辛亥,等. 煤田火灾防治理论与技术 [M]. 徐州:中国矿业大学出版社,2014. [15] IDE S T, ORR JR F M. Comparison of methods to estimate the rate of CO2 emissions and coal consumption from a coal fire near Durango, CO[J]. International Journal of Coal Geology, 2011, 86(1): 95−107. doi: 10.1016/j.coal.2010.12.005
[16] KREVOR S C M, IDE T, BENSON S M, et al. Real-time tracking of CO2 injected into a subsurface coal fire through high-frequency measurements of the 13CO2 signature[J]. Environmental Science & Technology, 2011, 45(9): 4176−4186.
[17] SONG Z, HUANG X, KUENZER C, et al. Chimney effect induced by smoldering fire in a U-shaped porous channel: A governing mechanism of the persistent underground coal fires[J]. Process Safety and Environmental Protection, 2020, 136: 136−147. doi: 10.1016/j.psep.2020.01.029
[18] HUANG J, BRUINING J, WOLF K-H. Modeling of gas flow and temperature fields in underground coal fires[J]. Fire Safety Journal, 2001, 36(5): 477−489. doi: 10.1016/S0379-7112(01)00003-0
[19] WOLF K-H, BRUINING H. Modelling the interaction between underground coal fires and their roof rocks[J]. Fuel, 2007, 86(17/18): 2761−2777.
[20] WESSLING S, KESSELS W, SCHMIDT M, et al. Investigating dynamic underground coal fires by means of numerical simulation[J]. Geophysical Journal International, 2008, 172(1): 439−454. doi: 10.1111/j.1365-246X.2007.03568.x
[21] 宋泽阳,朱红青,徐纪元,等. 地下煤火高温阶段贫氧不完全燃烧耗氧速率的计算[J]. 煤炭学报,2014,39(12):2439−2445. SONG Zeyang, ZHU Hongqing, XU Jiyuan, et al. An approach to calculate oxygen consumption rate of underground coal fires with lean oxygen concentration and incomplete combustion at high temperature[J]. Journal of China Coal Society, 2014, 39(12): 2439−2445.
[22] WANG Y, LI X, GUO Z. Prediction of self-ignition fire propagatioin and coal loss in an inclined seam[J]. Heat Transfer Research, 2018, 49(9): 827−845. doi: 10.1615/HeatTransRes.2018020029
[23] TANG Y, ZHONG X, LI G, et al. Simulation of dynamic temperature evolution in an underground coal fire area based on an optimised thermal–hydraulic–chemical model[J]. Combustion Theory and Modelling, 2019, 23(1): 127−146. doi: 10.1080/13647830.2018.1492742
[24] BUSTAMANTE RúA M O, DAZA ARAGóN A J, BUSTAMANTE BAENA P. A study of fire propagation in coal seam with numerical simulation of heat transfer and chemical reaction rate in mining field[J]. International Journal of Mining Science and Technology, 2019, 29(6): 873−879. doi: 10.1016/j.ijmst.2019.09.003
[25] 曾强. 新疆煤田火灾烟气流动与传热特性的研究[J]. 煤炭学报,2009,34(6):792−796. ZENG Qiang. Study on the regularity of gas/smoke flow and the characteristics of heat-transfer with spontaneous combustion in Xinjiang coalfield[J]. Journal of China Coal Society, 2009, 34(6): 792−796
[26] 曹代勇,樊新杰,吴查查,等. 内蒙古乌达煤田火区相关裂隙研究[J]. 煤炭学报,2009,34(8):1009−1014. CAO Daiyong, FAN Xinjie, WU Chacha, et al. Study on the fractures related with coalfield fire area in Wuda coalfield, Inner Mongolia[J]. Journal of China Coal Society, 2009, 34(8): 1009−1014.
[27] KUENZER C, STRACHER G B. Geomorphology of coal seam fires[J]. Geomorphology, 2012, 138(1): 209−222. doi: 10.1016/j.geomorph.2011.09.004
[28] ROY N C, Nakamura Y. Investigation of unsteady behaviors of forward and opposed flow combustion of solid fuel[J]. Combustion & Flame, 2016, 163: 517-531.
[29] SONG Z, FAN H, JIANG J, et al. Insight into effects of pore diffusion on smoldering kinetics of coal using a 4-step chemical reaction model[J]. Journal of Loss Prevention in the Process Industries, 2017, 48: 312−319. doi: 10.1016/j.jlp.2017.04.020
[30] HUANG X, REIN G. Smouldering combustion of peat in wildfires: inverse modelling of the drying and the thermal and oxidative decomposition kinetics[J]. Combust Flame, 2014, 161(6): 1633−1644. doi: 10.1016/j.combustflame.2013.12.013
[31] REIN G, LAUTENBERGER C, FERNANDEZ-PELLO A C, et al. Application of genetic algorithms and thermogravimetry to determine the kinetics of polyurethane foam in smoldering combustion[J]. Combust Flame, 2006, 146(1): 95−108.
[32] SCHLöMER S, TESCHNER M, POGGENBURG J, et al. Gas and temperature monitoring of a spontaneous coal seam fire in Wuda Coal Mining Area [C]//Beijing: Tsinghua University Press and Springer Verlag, 2008: 277-305.
[33] SCHAUMANN G. Geophysical investigation of Wuda coal mining area, Inner Mongolia: electromagnetics and magnetics for coal fire detection [C]//Beijing: Tsinghua University Press and Springer Verlag, 2008: 336-349.
[34] HUANG X, REIN G. Downward spread of smouldering peat fire: the role of moisture, density and oxygen supply[J]. International Journal of Wildland Fire, 2017, 26(11): 907−918. doi: 10.1071/WF16198
-
期刊类型引用(3)
1. 郭庆彪，谢扬，汪锋，郑美楠. 废弃采空区封存CO_2地表形变特征模拟研究. 煤矿安全. 2024(05): 1-10 . 
本站查看
2. 邹永洺. 煤层深部开采瓦斯赋存规律与临界采深研究. 煤. 2024(10): 10-13+55 . 百度学术
3. 许慧敏. 温度和压力对煤样瓦斯吸附解吸的影响规律研究. 煤. 2024(12): 5-8 . 百度学术
其他类型引用(1)
计量
- 文章访问数: 57
- HTML全文浏览量: 15
- PDF下载量: 10
- 被引次数: 4
