基于不同含水率煤的瓦斯解吸特性推算煤层瓦斯压力方法研究.pdf

8移动扫码阅读coalwithdifferentmoisturecontentsLJSafetinCoalMines.2023,54(9):8-14.FENG Yufeng,DONG Huzi.Estimation of coal seam gas pressure based on gas desorption characteristics of安全4(9):8-14.冯玉凤，董虎子.基于不同含水率煤的瓦斯解吸特性推算煤层瓦斯压力方法研究 .煤矿Safetyin Coal Mines2023年9 月Sep.2023No.9煤砺安全Vo1.54第9 期第54 卷DOI:10.13347/ki.mkaq.2023.09.002基于不同含水率煤的瓦斯解吸特性推算煤层瓦斯压力方法研究冯玉凤，董虎子（中国矿业大学（北京）应急管理与安全工程学院，北京10 0 0 8 3）摘要：针对直接法测定瓦斯压力工艺复杂、成功率低，而间接法测定适用条件较为苛刻的问题，以河南鹤壁三矿为例，基于不同含水率煤的瓦斯解吸实验规律，引入“孙重旭”式系数,提出了能够推算不同含水率煤层瓦斯压力的新方法。结果表明：“孙重旭式”拟合不同含水率煤样的解吸曲线效果最好，确定孙重旭公式为最佳拟合公式；系数在不同的拟合时间下，可保持稳定的值，且与瓦斯压力呈现指数关系，可利用系数与瓦斯压力的关系对不同含水率煤层瓦斯压力进行推算。现场验证发现实测值与理论计算值较为接近，该解吸法推算煤层瓦斯压力是可行的。关键词：含水率；瓦斯压力；解吸规律；孙重旭式系数；指数关系中图分类号：TD712文献标志码：A文章编号：10 0 3-4 9 6 X（2 0 2 3）0 9-0 0 0 8-0 7Estimation of coal seam gas pressure based on gas desorption characteristics of coalwith different moisture contentsFENGYufeng,DONGHuzi(School of Emergency Management and Safety Engineering,China University of Mining&Technology(Beijing),Beijing 100083,China)Abstract:The direct method to measure gas pressure process is complex,the success rate is low,the applicable conditions of the in-direct method are relatively harsh.Taking Hebi No.3 Mine as an example,based on the experimental law of gas desorption of coalwith different moisture contents,the coefficient a of the“Sun Zhongxu formula is introduced,and a new method that can estimatethe gas pressure of coal seams with different moisture contents is proposed.The results show that the“Sun Zhongxu equation is thebest way to fit the desorption curve of coal samples with different moisture contents,the“Sun Zhongxu equation is determined tobe the best fitting formula.The coefficient a can maintain a stable value under different fitting time and has an exponential relation-ship with gas pressure.The relationship between the coefficient a and gas pressure can be used to estimate the gas pressure of coalseams with different moisture contents.Field verification found that the measured value is close to the theoretical calculation value,the desorption method is feasible to calculate the coal seam gas pressure.Key words:moisture content;gas pressure;desorption law;coefficient a of Sun Zhongxu formula;exponential relationship煤层瓦斯压力是瓦斯在煤层中所呈现的压力，它与煤层瓦斯含量的测定及瓦斯抽放、瓦斯涌出与突出危险性预测等瓦斯防治及治理工作息息相关。煤层瓦斯压力的快速准确测定对于制定安全可靠的瓦斯治理方案和保障矿井的安全高效开采具有重要的意义 1-2 收稿日期：2 0 2 2-0 4-0 1责任编辑：兰莹基金项目：国家自然科学基金资助项目（52 130 4 0 9）作者简介：冯玉凤（19 8 8），女，内蒙古赤峰人，工程师，硕士，从事安全工程研究和实验室安全管理工作。E-mail：f e n g y u f e n g 8 8 16 3.c o m9SafetyinCoal MinesSep.20232023年9 月No.9煤砺发全Vol.54第9 期第54 卷目前常用的瓦斯压力测定方法分为直接法和间接法2 类 3。直接法是通过向煤层内部打钻孔，再利用封孔材料和测压仪来测定煤层瓦斯压力，封孔的严密性和封孔质量的好坏直接影响数据的准确定 4-6 。间接法指先测定煤层瓦斯含量，采用朗格缪尔吸附公式计算得出。由于公式中涉及煤层瓦斯含量、吸附常数、煤的水分和灰分、煤的真假密度等多个参数，影响这些参数测定准确性的因素较多，因此间接法测定煤层瓦斯压力也存在较多不可控因素 7 。赵发军等 8 考虑了煤层瓦斯的分析组分，认为煤对不同气体的吸附能力相差较大，如果分析组分中CH4浓度偏低却得不到修正，其测定结果可能有较大偏差；此外，在利用瓦斯压力解吸式计算煤层瓦斯压力方面贺建民等 9 从热力学和地应力梯度变化的角度探讨了煤层中瓦斯压力的理论计算方法；王宏图等 10 考虑了地应力和地温的影响，提出了煤层中瓦斯压力的理论计算式，但是式中所需参数繁多，计算较为复杂。煤炭开采中，通常利用直接法从岩巷向煤层打钻测定煤层瓦斯压力是可行的，但是直接法测定煤层瓦斯压力工艺复杂，受测压封孔技术和封孔质量影响较大 1，对于不具备穿岩测压的条件的煤层就要采用间接法测定煤层瓦斯压力。间接测定煤层瓦斯压力的方法具有用时短、操作方便等优点，为不具备穿岩测压条件或含水率较大不易采取直接打钻测压的煤层提供了测定瓦斯压力的重要手段，同时可对直接法测定瓦斯压力的测压数据的可靠性进行校验。而一般的间接测压方法通常受煤层特点和间接测压方法自身特性的影响，均有一定的适用条件，因此，为了能够较准确推算含水率较大煤层的瓦斯压力，进一步探索煤层瓦斯压力的间接测定方法是十分必要的。基于此，利用孙重旭解吸函数式，推导出一种瓦斯压力指数计算式，能够准确有效的推算煤层瓦斯压力。1煤层瓦斯压力与瓦斯解吸量的关系煤层瓦斯运移是个复杂的动力学过程，经历脱附/吸附、扩散、渗流3个过程 12 。当煤层中的瓦斯向外涌出时，瓦斯先通过脱吸附作用从煤分子表面运动到煤的孔裂隙空间中，且瓦斯的脱附是瞬间完成的；孔裂隙空间中的瓦斯分子再通过裂隙通道渗流到煤体外面，这个过程需要一定的时间。从不同的研究角度上看，瓦斯运移遵从不同的规律，但大多数研究者都赞同煤粒的瓦斯扩散运动符合菲克扩散定律，杨其銮等 13 正是基于菲克扩散定律，通过数学推导得到了煤屑瓦斯扩散规律的一般表达式，并通过阳泉一矿煤样瓦斯扩散率的测定验证了该表达式的正确性：681enBt(1)Q元2nn=1式中：Q,为t时刻累计的解吸气体量，cm/g;Q。为极限瓦斯解吸量，cm/g；B=元D/d；D 为扩散系数，cm/s；d 为煤粒的半径，cm；n 为常数，n=1，2，3，;t 为解吸时间。煤的吸附曲线呈抛物线状，可用对数公式近似表示为：X=Clnp+D(2)式中：X为煤层瓦斯含量，cm/g；p 为瓦斯压力，MPa；C、D 为拟合常数。1-Mad-AdQ=(Clnp+D-Xa)(3)1+0.31Mad式中：Xa为常压解吸量，cm/g；Ma 为煤中灰分，%；A为煤中水分，%。由式(1)式(3)可得：Q,=Eg(t)(Clnp+D)6人1g(t)=1enBt元22n=1D=D-Xa1-Mad-AdE=(4)1+0.31MadQ,可看作时间的解吸函数f(t)，代入式(4)可得：P=ecEg()C(5)由式(5)可以看出，煤层瓦斯压力与煤样的解吸规律、瓦斯解吸量间有着紧密的关系。如果确定了解吸函数f(t)，就可以推算出煤层的瓦斯压力。为了进一步研究不同含水率煤层的瓦斯压力与解吸函数间的关系，从现场采取煤样通过实验室制备了不同含水率的煤样进行了解吸规律的研究。2不同含水率煤样瓦斯解吸规律实验研究2.1实验方法和过程从河南鹤壁三矿采集煤样，在实验室破碎后筛分出粒径为0.18 0.2 5mm的煤样，采用干燥箱干燥、浸泡抽真空、恒温恒湿箱培养3种方式制10SafetyinCoal MinesSep.20232023年9 月No.9Vol.54煤砺岁全第9 期第54 卷备了干燥煤样、平衡水煤样（5.0 9%）、湿煤样（3.7 8%）和原始水分煤样（1.2 2%）4 种不同含水率的煤样进行解吸实验研究。瓦斯解吸实验装置示意图如图1，装置主要由真空脱气单元、充气单元、恒温控制单元和瓦斯吸附、解吸单元4 个部分组成。系统安装了2 套解吸仪，在2 个独立回路内可以同时进行2 个煤样的瓦斯吸附-解吸平行实验，提高了实验效率，同时最大限度地确保了实验条件（温度、吸附平衡压力等)的一致性。压力表压力表减压阀充气阀-2 0 平衡阀-2 0 解吸阀-2 变卡2一样品管阀-30 4D样品管阀-4压力表样品样品仪表阀压力表。压力表。罐罐控制阀变卡1解吸仪-2D平衡阀-1D4解吸阀-1充气阀-1样品D4样品管阀-2充气管阀管阀-1气瓶真空阀样品样品充气罐罐罐2000ml百解吸仪-1负压表真空泵图1瓦斯解吸实验装置示意图Fig.1Gas desorption experiment device在装有精密压力表的耐高压煤样罐中装人一定质量实验所需的煤样，在温控单元的恒温水域设定实验所需的固定温度；检查系统仪器的气密性，在保证气密性良好的条件下对实验系统的充气罐、煤样罐及所属的管路进行抽真空，连续抽真空12 h以上；向煤样罐中充入一定压力值的甲烷气体，煤样吸附4 8 h以上，充分达到吸附平衡后，瞬间释放煤样罐压力，开始模拟煤体暴露在煤巷中的瓦斯解吸规律的实验，力求最大限度地保障模拟条件与煤矿井下实际情况的相似性，从而测定煤样的常压解吸量。2.2实验结果鹤壁矿不同含水率煤样在不同吸附平衡压力下的瓦斯解吸量变化曲线如图2。由图2 可以看出，在相同的粒度（0.18 0.2 5mm）、恒定的实验系统温度30（精度1)条件下，吸附平衡压力对不同含水率煤的瓦斯解吸规律有如下影响：对于同一种煤，不同含水率的煤样与干燥煤具有相同的瓦斯解吸规律；对于同一含水率的煤样，吸附平衡压力越大，其初始阶段的瓦斯解吸速度越大，在相同时间段内的瓦斯解吸量越大，且随着吸附平衡压力的增加，解吸曲线上移的幅度变大。目前为止，国内外众多学者通过对煤瓦斯解吸规律的大量研究，提出了许多经验或半经验公式。但是由于研究对象和角度不同，提出的经验公式既有合理的部分也有不同的局限性存在，并不适合所有煤屑的解吸规律 4 。为寻求描述适合不同含水率煤样的瓦斯解吸规律公式，通过巴雷尔式、乌斯基诺夫式等8 类具有代表性的经验公式，分别对鹤壁矿不同含水率煤样的解吸曲线前3、10、2 0、6 0、12 0 m i n 的瓦斯累计量进行拟合分析，其中吸附平衡压力为1.5MPa、含水率1.2 2%的鹤壁矿煤样不同时间段瓦斯解吸量的回归分析见表1。从表1中的拟合相关系数R来看，上述几种描述煤瓦斯解吸过程的理论和经验公式中，从整体效果来看以孙重旭式最为理想（文特式和孙重旭式表达式子不同，但在本质上是同一类公式），拟合精度最高且较稳定，随拟合时间段的延长拟合精度仍然较高，适用于水分变化较大煤层的瓦斯解吸规律的研究。对其他鹤壁矿不同含水率的煤样在不同吸附平衡压力下的实测瓦斯解吸数据的拟合效果也支持上述结论，受限于篇幅，未附其11SafetyinCoal Mines2023年9 月Sep.2023煤砺发全No.9Vol.54第9 期第54 卷0.8MPa0.8MPa201.5MPa161.5 MPa602.5 MPa602.5 MPa3.5MPa14）3.5 MPa4.5MPa124.5MPa15101086542020406080100120020406080100120t/mint/min（a）鹤壁干燥煤样（b）鹤壁含水率1.2 2%煤样0.8MPa0.8MPa141.5 MPa81.5 MPa(8.2.5 MPa2.5MPa3.5MPa712）3.5MPa4.5MPa4.5 MPa61054863422020406080100120020406080100120t/mint/min（c）鹤壁含水率3.7 8%煤样（d）鹤壁含水率5.0 9%煤样图2鹤壁矿不同含水率煤样在不同吸附平衡压力下的瓦斯解吸量变化曲线Fig.2Gas desorption volume curves at different adsorption balance pressures indifferent moisture conditions of coal in Hebi Coal Mine表1鹤壁矿煤样不同时间段瓦斯解吸量的回归分析(1.5MPa、含水率1.2 2%)Table1 Regression analysis results of the desorptionamount of coal samples collected from Hebi Coal Minein different time(1.5 MPa,moisture content 1.22%)R2公式名称回归公式形式3min 10 min 20min 60 min 120min巴雷尔式Qi=Ki Vi0.9910.9250.8720.7310.519V1文特式Qt1-K,t1-K,0.9990.9990.9990.9970.991乌斯基(1+t)I-nQ=Vo0.9320.8820.8840.897 0.889诺夫式1-nQ博特式1-Ae-At0.8710.8440.8810.8960.881QABt王佑安式Q0.9980.9940.9910.9910.9901+Bt孙重旭式Qi=ati0.9990.9990.999 0.9970.991Vo指数式Qe-bt)0.8160.5400.4310.2290.061b对数式Q,=Alnt+B0.9990.9920.9870.9840.988拟合效果的图表。同时以相同方法分别对其他变质程度的义马千秋矿（长焰煤）、山西马兰矿（1/3焦煤）、焦作古汉山矿（无烟煤)的不同含水率的煤样在不同吸附平衡压力下解吸数据进行拟合，得出各组数据拟合的效果基本相似，也符合上述结论，且拟合效果与煤样的含水率大小并无直接关系，得出的结论是适用于不同变质程度、不同含水率的煤样，说明孙重旭式作为煤样的解吸函数具有一定的普遍性，也证实了瓦斯解吸模拟实验数据的可靠性。利用孙重旭式拟合不同吸附平衡压力下鹤壁矿含水率1.2 2%煤样解吸曲线，鹤壁矿煤样孙重旭式拟合系数见表2，因其他含水率煤样解吸量拟合回归系数、i 值规律相同，不再详列。从表2 可以看出：同一含水率的煤样回归系数i值随拟合时间不同、吸附平衡压力不同其拟合值是稳定的，可以看作与解吸时间、吸附平衡压12Safetyin Coal MinesSep.20232023年9 月No.9Vol.54煤防岁全第9 期第54 卷表2鹤壁矿煤样孙重旭式拟合系数（含水率1.2 2%）Table 2 Fitting coefficients of formula created by SunZhongxu of Hebi Coal Mine(Moisture content 1.22%)时间0.8MPa1.5MPa2.5MPa3.5MPa4.5MPaminaiaiaiaai31.7340.2912.7010.2913.7810.2944.4420.2915.1490.291101.7340.2952.716.0.2903.7900.2894.4620.2925.13990.292201.741 0.2872.710 0.2923.792 0.2924.4490.2845.1420.294601.7320.2902.7080.28993.7950.2854.5320.2915.141(0.2861201.7300.2882.713 0.2903.791 0.2944.676 0.2885.1430.290力基本无关；回归系数与解吸时间无关，而与吸附平衡压力密切相关，且能体现同一含水率煤层不同瓦斯压力下煤样解吸特性的差异。因此，在井下同一种含水率煤层不同取样损失时间下都可以对现场解吸曲线进行孙重旭式拟合，得到回归系数来推算取样地点煤层瓦斯压力。不同含水率的煤样瓦斯解吸曲线符合孙重旭式，故解吸函数可写作f(t)=at，代人式(5)可得：P三eCEg()(6)由于系数与瓦斯压力密切相关，而与解吸时间无关，系数i与解吸时间、吸附平衡压力均无关，因此，在解吸时间段确定时，瓦斯压力和系数的关系可表示为：p=aela(7)式中：=e-%；BCEgo、为解吸时间t的函数，解吸时间t确定时、为常数。鹤壁矿含水率1.2 2%煤样拟合时间段为前3min、前2 0 min、6 0 m i n 及前12 0 min时，系数a与吸附平衡压力的回归分析如图3。4.54.53.53.5p=0.356 2e0.505 3ap=0.3544e0.507.0aR2=0.981 1R2=0.981 72.52.51.51.50.860.81.732.73.74.45.11.732.73.74.45.1系数系数(a)前 3 min(b)前 2 0 min4.543.53.5p=0.3550e0.5069ap=0.366 0e0.493 8a2.5R2=0.98152.5R2=0.982 51.51.50.80.81.732.73.74.45.11.732.73.74.45.1系数a系数a(c)前 6 0 min(d)前 12 0 min图3系数a与吸附平衡压力的回归分析Fig.3Regression analysis of fitting coefficient a with adsorption equilibrium pressure13Safety in Coal MinesSep.20232023年9 月煤砺发全No.9Vol.54第9 期第54 卷图3中指数关系拟合的相关系数R均较高，拟合效果较好。因此，根据煤样在不同吸附平衡压力下的解吸实验规律，可以建立系数与吸附平衡压力的指数函数公式(式(7)）。对任意取样损失时间下的解吸曲线进行孙重旭式拟合得到回归系数，即可推算出所测地点的煤层瓦斯压力。3应用验证鹤壁三矿二，煤层为石炭二叠系，煤层结构简单，无夹石层，煤质为贫瘦煤。实验室研究瓦斯解吸规律实验所用的鹤壁矿新鲜原煤煤样取自该煤层的四水平南翼辅助回风巷4 10 1工作面的2#穿层钻孔。在现场对4 10 1工作面瓦斯压力进行测定时，采用直接测压法中的被动测压法，在辅助回风巷顶板中部向上依次打3个孔径为9 4 mm的穿层钻孔，钻孔倾角4 0，穿层钻孔依次编号命名为1#3#穿层钻孔。经分析选用了2#穿层钻孔进行取样，将取到的煤屑迅速筛选出粒径为0.18 0.2 5mm的煤样，装人解吸罐中进行现场瓦斯解吸测定，现场实测瓦斯解吸量随时间变化的曲线如图4。765432020406080100120时间/min图4 王现场实测瓦斯解吸量随时间变化的曲线Fig.4Curve of the measured gas desorption amountwith time利用孙重旭式对现场解吸曲线进行拟合，可以得到拟合函数Q=1.76r0.277，R为0.9 9 1，拟合系数a=1.76，i=0.2 7 7。将=1.76代人到图3中拟合的指数式，分别得到拟合时间段为前3min、前2 0 min、前6 0 min及前12 0 min时，推算的煤层瓦斯压力依次为0.8 6 2、0.8 6 4、0.8 6 6、0.8 7 5MPa。此煤层相同埋藏深度辅助回风巷的工作面1#3#穿层钻孔所测得的瓦斯压力结果分别为0.8 8、0.92、0.8 9 MPa，现场实测值与推算值比较接近，结果表明，利用解吸法推算煤层瓦斯压力是准确可靠的。4结语1)通过对鹤壁三矿煤样解吸数据分析发现孙重旭式是描述不同含水率煤的瓦斯解吸特性的最佳经验公式，可作为煤样的解吸函数。2)对于同一含水率的煤样，孙重旭式中回归系数i值与拟合时间段和吸附平衡压力无关，始终保持稳定；回归系数与解吸时间无关，而与吸附平衡压力密切相关，且能体现同一含水率煤层不同瓦斯压力下煤样解吸特性的差异。3)基于瓦斯压力与瓦斯解吸量之间的数学关系建立孙重旭式系数与瓦斯压力之间的指数函数关系式p=ea。对于不同含水率的煤层，通过现场煤样的瓦斯解吸规律，得到孙重旭式拟合参数，即可推算所测地点的煤层瓦斯压力。4)在鹤壁三矿二，煤层含水率为1.2 2%的原煤进行了现场效果检验，验证结果表明利用该方法推算煤层瓦斯压力是可行的。参考文献(References):1王亮，孙毅民，褚鹏，等.基于时空分布特征的煤层瓦斯测压准确性研究 J.中国安全科学学报，2 0 2 1,31(2):40-47.WANG Liang,SUN Yimin,CHU Peng,et al.Study onaccuracy of coal seam gas pressure measurement basedon its spatial and temporal distribution characteristicsJ.China Safety Science Journal,2021,31(2):40-47.2邹永洺.煤层瓦斯压力快速测定技术应用研究 J.煤矿安全,2 0 19,50(4)：138-14 1.ZOU Yongming.Application research on rapid determ-ination technology of coal seam gas pressureJ.Safetyin Coal Mines,2019,50(4):138-141.3崔永杰，王兆丰，范道鹏.免封孔煤层瓦斯压力测定技术可行性研究 J.煤矿安全,2 0 2 1,52(11)：2 2-2 7.CUI Yongjie,WANG Zhaofeng,FAN Daopeng.Rapidmeasuring method and device of working face residualgas parametersJ.Safety in Coal Mines,2021,52(11):22-27.4李皓伟，刘海，王兆丰，等.近水平穿层钻孔两级变径瓦斯压力测定封孔技术 J.煤炭科学技术，2 0 17,45(10):101-104.LI Haowei,LIU Hai,WANG Zhaofeng,et al.Boreholesealing technology of gas pressure measurement in twostage variable diameter borehole with subhorizontal14SafetyinCoal MinesSep.20232023年9 月No.9煤砺发全Vol.54第9 期第54 卷passed through strata drillingJ.Coal Science andTechnology,2017,45(10):101-104.5安丰华，程远平，吴冬梅，等.基于瓦斯解吸特性推算煤层瓦斯压力的方法 J.采矿与安全工程学报，2 0 11,28(1):81-85.AN Fenghua,CHENG Yuanping,WU Dongmei,et al.Determination of coal gas pressure based on character-istics of gas desorption J.Journal of Mining&SafetyEngineering,2011,28(1):81-85.6高明忠，陈领，凡东，等.深部煤矿原位保压保瓦斯取芯原理与技术探索 J.煤炭学报，2 0 2 1，4 6(3)：885-897.GAO Mingzhong,CHEN Ling,FAN Dong,et al.Prin-ciple and technology of coring with in-situ pressure andgas maintaining in deep coal mine J.Journal of ChinaCoal Society,2021,46(3):885-897.7张东旭.外加水分对煤层瓦斯压力测定的影响分析J.矿业安全与环保,2 0 2 1,4 8（4)：118-12 2.ZHANG Dongxu.Analysis of the influence of addition-al water on the measurement of coal seam gaspressureJ.Mining Safety&Environmental Protec-tion,2021,48(4):118-122.8赵发军，温志辉，刘明举，等.分析组分对间接法测定煤层瓦斯压力的影响 J.煤炭科学技术，2 0 10,38（4)：52-54.ZHAO Fajun,WEN Zhihui,LIU Mingju,et al.Analys-is content affected to seam gas pressure measure withindirect methodJ.Coal Science and Technology,2010,38(4):52-54.9贺建民，王宏图，鲜学福，等.煤层温度和应力梯度变化对煤层瓦斯压力计算的影响 J，重庆大学学报（自然科学版），19 9 9,2 2(5)：9 5-9 8.HE Jianmin,WANG Hongtu,XIAN Xuefu,et al.Theinfluence of in-situ stress gradient and geothermal tem-perature on the calculation of coalbed gas pressureJ.Journal of Chongqing University(Natural Science Edi-tion),1999,22(5):95-98.10 王宏图，鲜学福，尹光志，等.煤矿深部开采瓦斯压力计算的解析算法 J.煤炭学报，19 9 9,2 4(3)：2 7 9-2 8 3.WANG Hongtu,XIAN Xuefu,YIN Guangzhi,et al.In-vestigation of analytic method for coalbed gas pressurein high temperature and deeply mining mines J.Journ-al of China Coal Society,1999,24(3):279-283.11赵嵘，齐黎明.井下直接法测定煤层瓦斯压力研究现状及分析 J.煤炭技术，2 0 17,36(3)：2 0 2-2 0 3.ZHAO Rong,QI Liming.Research status of directmethod of gas pressure measurement in coal seam un-der mine and its analysisJ.Coal Technology,2017,36(3):202-203.【12 贾宏福，安丰华，彭信山.巴雷尔法计算煤瓦斯扩散系数的局限性分析 J.中国安全科学学报，2 0 2 0,30(2)：73-78.JIA Hongfu,AN Fenghua,PENG Xinshan.Limitationof Barrell formula in calculating gas diffusion coeffi-cient of coalJ.China Safety Science Journal,2020,30(2):73-78.13杨其銮，王佑安.煤屑瓦斯扩散理论及其应用 J.煤炭学报,19 8 6,11(3):8 7-9 4.YANG Qiluan,WANG Youan.Theory of methane dif-fusion from coal cuttings and its applicationJ.Journalof China Coal Society,1986,11(3):87-94.14】秦跃平，史浩洋，乔廷，等.关于煤粒瓦斯解吸经验公式的探讨 J.矿业安全与环保，2 0 15,4 2(1)：10 9-111.QIN Yueping,SHI Haoyang,QIAO Ting,et al.Discus-sion on empirical formula for gas desorption of coalparticlesJ.Mining Safety&Environmental Protec-tion,2015,42(1):109-111.