论文部分内容阅读
摘 要：为掌握上邻近保护层开采时Y型通风工作面采空区的瓦斯分布规律，并控制采空区的瓦斯涌出，文中首先针对平煤六矿戊8-32010工作面建立了采空区瓦斯涌出理论模型，依据现场情况设置瓦斯涌出条件并设置了机巷风量与风巷风量配比分别为1∶3，1∶2和1∶1.5时的入口风速，建立了瓦斯-空气混合气体在采空区里运移的方程组，利用FLUENT软件进行了这3种风量配比条件下的采空区瓦斯分布规律的数值试验，将此结果与现场工业性试验结果进行了对比分析，确定了适合该工作面的风量分配，最后利用数值模拟分析了留巷钻孔法的抽采效果。结果表明：越往Y型通风采空区以深，瓦斯浓度越高，且采空区尾部风巷侧体积分数高于尾部留巷侧;风巷风量越大，工作面瓦斯浓度越小，风巷風量越小，工作面和留巷瓦斯浓度越大;在留巷向采空区尾部施工抽采钻孔可以实现瓦斯的高效高浓抽采。
关键词：安全科学与工程;瓦斯分布规律;FLUENT数值模型;上邻近保护层开采;风量分配比例
中图分类号：TD 712
文献标志码：A
文章编号：1672-9315（2021）01-0015-08
DOI：10.13800/j.cnki.xakjdxxb.2021.0103
Gas distribution law and control of Y-shaped ventilated goaf
in mining upper adjacent protective layer
——A case study of Fifth8-32010 working face in Pingmei No.6 Mine
XUE Junhua1，MA Qian1，LI Yanhe2，LI Hongbiao2，YUAN Zhandong2
（1.College of Safety Science and Engineering，Xi’an University of
Science and Technology，Xi’an 710054，China;
2.No.6 Mine of Pingdingshan Tian’an Coal Mining Co.，Ltd.，Pingdingshan 467000，China）
Abstract：In order to grasp the gas distribution law in the goaf of the Y-shaped ventilation working face during the mining of the upper adjacent protective layer，and control the gas emission in the goaf，this article first established
theoretical model
of the gas emission in the goaf for the Fifth8-32010 working face of Pingdingshan Coal Mine.Gas emission conditions were set according to the site situation and the inlet wind speed was set as follows：the ratio of the air volume in the engine lane to the air volume in the wind tunnel is 1∶3，1∶2，and 1∶1.5.A set of equations of gas-air mixture transition in goaf using FLUENT software were established.The numerical test was carried out of the gas distribution law in the goaf under the conditions of the three air volume ratios for the equations moving in the goaf.The results are compared and analyzed with the ones of the on-site industrial test.The air volume distribution suitable for the working face is calculated，and finally the drainage effect of the retaining borehole method is analyzed by numerical simulation.The results show that：the deeper the Y-shaped ventilated goaf，the higher the gas concentration，and the volume fraction at the tail side of the goaf is higher than that at the tail side of the roadway;the greater the air volume in the wind road，the lower the gas concentration in the working face.The smaller the air volume，the greater the concentration of gas in the working face and in the retained lane;the construction of the drainage borehole at the end of the retained lane to the goaf can achieve efficient and high-concentration gas drainage.Key words：safety science and engineering;gas distribution law;FLUENT numerical model;mining of upper adjacent protective layer;air volume distribution ratio 0 引 言
高瓦斯煤层在井工开采过程中，工作面一直面临着采空区和上隅角的瓦斯积聚问题[1-2]，尤其是近距离上保护层开采工作面，由于开采扰动后下邻近层卸压瓦斯窜入采空区，导致采空区和上隅角瓦斯浓度居高不下，严重妨碍着矿井安全生产。Y型通风系统可以消除上隅角瓦斯积聚，解决上隅角瓦斯易于超限的问题，是一条实现安全、高效、经济一体化生产的有效途径[3]。
虽然Y型通风系统解决了上隅角瓦斯超限问题，但同时改变了采空区瓦斯浓度场[4-5]，所以在施工过程中不能套用以前U型通风的瓦斯浓度场分布规律以及瓦斯治理技术[6-7]。只有摸清楚瓦斯浓度的分布规律，找出适合该工作面的系统风量及2条进风巷的风量分配，确定出高浓度瓦斯区域，在留巷内向此区域施工抽采钻孔才是解决采空区瓦斯浓度高的科学方法[8-9]。
为此，众多学者从理论研究、数值模拟、实验室试验和现场试验等方面进行了大量研究，获得了Y型通风采空区瓦斯浓度场分布规律[10]，偏Y型通风采空区瓦斯涌出规律[11]，开发了Y型通风风控瓦斯技术、留巷钻孔法抽采瓦斯技术等[12]。但是由于理论研究将现场复杂的条件做了理想简化及假设，一般不能很好的指导现场试验[13-14]。数值模拟试验在理论研究的基础上也做了一些条件假设，会与现场的试验有些差异，但是其适用性较强[15-17]。只有结合理论研究和数值模拟试验，才可以较好的指导现场工业性实验。
为了研究平顶山矿区煤层群特征下的上邻近保护层开采Y型通风采空区瓦斯浓度场分布规律及控制技术，文中基于理论研究，建立了采空区瓦斯涌出模型，并导入到FLUENT数值模型之中，研究了上邻近保护层Y型通风采空区瓦斯分布规律，并以进风巷风量分配为自变量，研究了其对瓦斯分布规律的影响，结合现场试验，确定最合适的进风巷风量分配，并在现场进行了验证，对比分析了导致数值模擬试验和现场试验结果之间差异的原因，最后利用留巷钻孔法对抽采效果进行考察，以期能为同类型煤层条件下的采空区瓦斯治理提供一定的参考依据。
1 矿井及工作面概况
平顶山天安煤业股份有限公司六矿主采煤层从上至下分别为丁5-6煤层、戊8煤层和戊9-10煤层，丁5-6煤层和戊8煤层平均层间距71 m，戊8煤层和戊9-10煤层平均层间距5 m，经煤层突出危险性鉴定后，决定将戊8煤层作为首采层。
戊8-32010工作面位于戊8煤层，平均煤厚2.2 m，设计可采长度2 300 m，倾向长度220 m，标高-633～-768 m，埋深873～938 m，通风系统结合高位瓦斯抽采巷及机巷切顶留巷布置，采用两进一回Y型通风方式，机巷、风巷进风，机巷留巷和机巷高抽巷回风，煤层为近水平煤层，遗煤平均厚度为0.2 m，可解吸瓦斯含量为4.3 m3/t，煤壁瓦斯初始涌出速率为0.003 56 m3/（m2·min）。工作面布置示意图如图1所示。
2 采空区瓦斯浓度分布规律及抽采数值模拟
2.1 条件假设
1）采空区为冒落矸石和岩块构成的带有空隙的多孔介质，并将多孔介质视为各向同性。
2）采空区气体不可压缩，不产生化学反应，其流动近似为稳态流动、等温过程。
3）忽略瓦斯的浮升效应，瓦斯充满整个采空区。
4）采空区瓦斯涌出源来自邻近层和遗煤。
5）遗煤和邻近层产生的瓦斯都是瞬间释放的，忽略解吸时间。
2.2 模型建立及网格划分
模拟对象为戊8-32010工作面，工作面长度220 m，采空区深度300 m，据现场考察，竖向裂隙带高度稳定在57 m，所以采空区高度设置为57 m，工作面宽度10 m，工作面高度4 m，风巷、机巷长度
20 m，宽度4 m，高度3.5 m，建成的模型如图2所示。
风巷、机巷和留巷网格间距为1 m×1 m，采空区模块网格尺寸为2 m×2 m，网格划分结果如图3所示。
2.3 采空区瓦斯涌出源模型
2.3.1 采空区空隙率分布模型
采空区空隙率与渗透率分布规律符合“O”形圈覆岩垮落形态，采空区空隙率与渗透率分布的数学表达如公式（1）和（2）[18]。
式中，n（x，y，z）为空隙率;L为工作面长度，m;k（x，y，z）为渗透率，m2;Dp为采空区冒落岩块平均粒径，取250 mm[19].
2.3.2 采空区遗煤瓦斯放散模型
采空区不同位置处的瓦斯涌出速率qc（y）可用公式组（3）求得[20]。
式中 a为瓦斯涌出初始强度，m3/min;b为瓦斯涌出衰减系数，min-1;y为遗煤位置到工作面的距离，m;Vx为工作面推进速度，m/min;Qt为单位体积遗煤可解吸瓦斯放散体积，m3/m3;T为解吸时间，取288 000 min（即200 d）.
通过上述方程组结合工作面遗煤放散条件，即可求出遗煤瓦斯涌出速率。
2.3.3 下邻近层瓦斯涌出模型
戊8-32010工作面开采时，戊9-10卸压瓦斯涌入工作面采空区，下邻近层瓦斯涌出量可由公式（4）求得，其中工作面煤壁绝对瓦斯涌出量参数，可由公式（5）求得，瓦斯排放率可由公式（6）求得，最终得出下邻近层的瓦斯涌出量[11]。
（6）
式中 Qs为戊9-10煤层量瓦斯涌出量，m3/min;mi 为戊9-10煤层厚度，3.1 m;Qb 为戊8-32010工作面绝对瓦斯涌出量，m3/min;hi 为邻近层层间距，m;其余薄煤层瓦斯涌出量忽略不计;ηs为戊9-10煤层的瓦斯排放率;L为戊8-32010工作面长度，m;h为开采厚度，2.2 m;S为开采层与邻近层层间距，13 m;q0为工作面煤壁刚暴露时单位面积煤壁的瓦斯涌出强度，m3/（m2·min）;V为工作面采煤机割煤时期平均牵引速度，m/min，实测值为10.08 m/min;Tc为工作面正常生产班期间，采煤机完成一次割煤周期所用时间，min，实测值为44 min。 2.4 采空区瓦斯运移方程组选择
选择的气体组分为甲烷-空气混合气体，选择湍流方程、连续性方程和动量守恒方程作为混合气体在采空区内的流动方程。假设混气扩散符合Fick
定律，其中混合气体在采空区中的流动方程
的统一表达式如式（7）所示，Fick扩散定律的表达式如式（8）所示[21]。
式中 ρ为气体密度，kg/m3;为通用变量;t为时间，s;
Γ为广义扩散系数;S为广义源项。
式中 r为极坐标半径，m;D为扩散系数，m2/s;c为瓦斯浓度，mol/m3。
2.5 边界条件及相关参数设置
为研究系统风量为2 830 m3/min时，两进风巷进风配比对工作面及留巷的瓦斯浓度的影响，共设置了2个进风入口和一个出口。2个进风入口，分别表示戊8工作面机巷和风巷，工作面机巷风量与机巷风量之比分别为：1∶3，1∶2和1∶1.5，一个出口设置为自然出流。
2.6 抽采条件设置
在留巷向里距离工作面切顶线10 m处增加抽放汇对采空区瓦斯进行抽采，抽采管直径为500 mm，钻孔高度为巷道底部2 m处。
3 数值模拟结果及分析
3.1 瓦斯浓度分布规律实验结果及分析
在总风量2 830 m3/min的基础上，对工作面机、风巷风量分配进行调整试验，分别对机巷、风巷风量比例1∶3，1∶2和1∶1.53种状态对工作面瓦斯质量分数的影响进行研究，结果如图4所示。
由图4可得，配风比例为1∶3时，通过工作面的风量较大，上下隅角压差增大，风巷端口与留巷压差也随之增大，导致工作面漏风严重，造成采空区防灭火难度加大;配风量为1∶1.5时，采空区和留巷处瓦斯浓度区域面积明显增大，基本可以看出增大风巷风量比例，工作面漏风加剧，减少风巷风量，采空区和留巷内瓦斯浓度会增大，而当采用配风比例1∶2时，一方面可以令工作面漏风情况得以控制，另一方面也可以缓解采空区和留巷内瓦斯积聚，所以配风比例为1∶2时更适合于戊8-32010工作面。
在Y型通风条件下，瓦斯向采空区深部不断积聚，且采空区尾部风巷侧瓦斯体积分数高于尾部留巷侧，由卸压钻孔法抽采卸压瓦斯机理可知[12]，采空区和覆岩裂隙区相连通，在留巷后部向采空区布置瓦斯抽采钻孔，可以持续高浓度抽采瓦斯资源。
3.2 瓦斯抽采实验结果及分析
采用配风比例为1∶2的通风模式建立戊8-32010工作面采空区抽采模型并布置尾抽钻孔，抽采模型如图5所示，抽采结果如图6所示。
由图6可以看出，利用尾抽钻孔抽采瓦斯后，采空区瓦斯体积分数明显降低，说明在留巷深部向采空区布置钻孔可以高效抽采瓦斯。
4 讨 论
本次数值实验结果以瓦斯质量分数显示，现场工业性实验结果以瓦斯体积分数显示。假设采空区瓦斯质量分数为cg，体积浓度为c′g，则质量分数和体积浓度存在以下关系[22]
（9）
式中 Ma为空气摩尔质量，kg/mol;Mg为瓦斯摩尔质量，kg/mol。
在总风量2 830 m3/min的基础上，对工作面机、风巷风量分配进行调整试验，分别对机巷、风巷风量比例
1∶3，1∶2和1∶1.5这3种状态时，对工作面3个测点和留巷内2个测点共5个测点的瓦斯浓度进行对比分析，测点布置如图7所示，现场试验测点与对应测试位置见表1。
先将现场试验和数值模拟实验5个测点位置的瓦斯体积分数统计、并绘图，瓦斯体积浓度对比结果如图8所示。
由图8（a）可得，机、风巷风量配比1∶3时工作面瓦斯浓度在0.16%～0.3%之间，留巷段瓦斯浓度在0.32%～0.69%之间;配比为1∶1.5时工作面瓦斯浓度在0.24%～0.5%之间，留巷段瓦斯浓度在0.51%～0.64%之間。由图8（b）可得，机、风巷风量配比1∶3时工作面瓦斯浓度在0.02%～
0.05%之间，留巷段瓦斯浓度在0.05%～0.3%之间;配比为1∶1.5时工作面瓦斯浓度在0.05%～0.07%之间，留巷段瓦斯浓度在0.04%～0.28%之间。
综上可得，风巷风量比例越小，工作面瓦斯浓度越高，留巷段瓦斯浓度越高;风巷风量比例越大，工作面瓦斯浓度越低，现场试验和数值模拟试验采空区瓦斯浓度的变化趋势基本一致。但工作面处和留巷内瓦斯浓度值差异较大，造成这一差异的主要原因是数值模拟中无法控制工作面向采空区内漏风和采空区向留巷内漏风的过程，后续笔者会就这一问题做进一步研究。
最后在留巷段向里距离工作面切顶线10 m处施工了瓦斯尾抽钻孔，抽采了采空区瓦斯，在此期间，工作面和采空区瓦斯都处于安全值以内[23-24]，相比于上个相邻U型工作面推进速度，本工作面推进速度明显提高，如图9为日产量与抽采浓度之间的关系图。
从图9可以看出，利用尾抽抽采采空区尾部瓦斯，瓦斯浓度都在50%以上，说明在尾部布置瓦斯抽采效果较好。在抽采效果考察后期，日产量有短时间的相对减小，这是由于此时煤层上覆岩层破断垮落，裂隙带初步形成，将尾抽为主的抽采方式转化为高抽巷裂隙带瓦斯为主抽采方式，一部分瓦斯通过裂隙带进行抽采，暂时降低了日产量。
5 结 论
1）在采空区瓦斯涌出和下邻近层瓦斯涌出强度一定时，近距离上保护层Y型通风的瓦斯易于积聚在采空区尾部，越靠近尾部，瓦斯体积分数增加速率越大，采空区尾部留巷侧瓦斯体积分数低于尾部风巷侧，在留巷尾部一定位置处布置瓦斯抽采钻孔能高效抽采瓦斯，并且此浓度场与单煤层Y型通风瓦斯浓度场分布规律基本一致。
2）增大风巷风量比例，会导致采空区瓦斯大量涌入留巷;减小风巷风量比例，会导致工作面瓦斯增加，因此改变2条进风巷配风量，工作面和留巷的瓦斯浓度会随之变化。 参考文献（References）：
[1] 姜福兴，孔令海，刘春刚.特厚煤层综放采场瓦斯运移规律[J].煤炭学报，2011，36（3）：407-411.
JIANG Fuxing，KONG Linghai，LIU Chungang.Gas emission laws of fully mechanized sublevel caving mining in extra thick coal seam[J].Journal of China Coal Society，2011，36（3）：407-411.
[2]BALUSU R，DEGUCHI G，HOLLAND R，et al.Goaf gas flow mechanics and development of gas and Sponcom control strategies at a highly gassy mine[J].Coal and Safety，2015，20：35-45.
[3]贾廷贵，娄和壮，闻文.不同配风条件下“Y”型通风采场瓦斯涌出规律[J].辽宁工程技术大学学报（自然科学版），2019，38（1）：14-18.
JIA Tinggui，LOU Hezhuang，WEN wen.Gas migration law of Y-type ventilation stope under different air distribution conditions[J].Journal of Liaoning Technical University（Natural Science），2019，38（1）：14-18.
[4]袁亮，薛俊华.低透气性煤层群无煤柱煤与瓦斯共采关键技术[J].煤炭科学技术，2013，41（1）：5-11.
YUAN Liang，XUE Junhua.Key technology of pillarless coal and gas simultaneous mining in low permeability seam group[J].Coal Science and Technology，2013，41（1）：5-11.
[5]戚良锋.Y型通风方式下采空区瓦斯流场数值模拟研究[D].淮南：安徽理工大学，2009.
QI Liangfeng.Research on gas field in gob with Y-type ventilation system by numerical simulation[D].Huainan：
Anhui University of Science and Technology，2009.
[6]陈胜利.综放工作面瓦斯涌出特征实测及分析[J].科学技术与工程，2012，27（12）：6892-6896.
CHEN Shengli.Measurement and analysis for gas emission characters in fully mechanized caving face[J].Science Technology and Engineering，2012，27（12）：6892-6896.
[7]徐玉胜，李春元，张勇，等.大采高采场瓦斯治理模型分析及通风系统优化研究[J].煤炭科学技术，2019，47（3）：142-149.
XU Yusheng，LI Chunyuan，ZHANG Yong，et al.Study on gas control model analysis and ventilation system optimization in large mining height mines[J].Coal Science and Technology，2019，47（3）：142-149.
[8]周一力.“Y+高抽巷”工作面采空区瓦斯与氧气浓度场分布规律及其在灾害防治中的应用[D].徐州：中国矿业大学，2019.
ZHOU Yili.Distribution of gas and oxygen concentration field in the goaf of “Y+High Pumping Alley”working face and its applications in disaster prevention[D].Xuzhou：China University of Miing and Technology，2019.
[9]常海雷，刘彦青，张浪，等.“U+I”型通风形式下综放工作面采空区遗煤及支架上方煤壁瓦斯放散与运移规律研究[J].矿业安全与环保，2017，44（2）：39-44.
CHANG Hailei，LIU Yanqing，ZHANG Lang，et al.Study on gas emission and migration law of residual coal in gob and coal wall above support in fully mechanized caving face with “U+I” shape ventilation system[J].Mining Safety & Environmental Protection，2017，44（2）：39-44.
[10]張欢，赵洪宝，李晓白，等.分段留巷Y型通风采场瓦斯与流场分布三维实测与重构[J].中国安全生产科学技术，2018，14（7）：107-114.
ZHANG Huan，ZHAO Hongbao，LI Xiaobai，et al.Three-dimensional measurement and reconstruction of flow field and gas distribution for mining face and goaf of Y type ventilation[J].Journal of Safety Science and Technology，2018，14（7）：107-114. [11]赵灿，张浪，刘彦青.偏Y型通风下采空区瓦斯涌出规律及超限治理研究[J].煤炭科学技术，2019，47（4）：127-133.
ZHAO Can，ZHANG Lang，LIU Yanqing.Study on gob gas emission laws and control of gas exceeding limit under inclined type Y-ventilation[J].Coal Science and Technology，2019，47（4）：127-133.
[12]袁亮.留巷钻孔法煤与瓦斯共采技术[J].煤炭学报，2008，33（8）：898-902.
YUAN Liang.The technique of coal mining and gas extraction by roadway retaining and borehole drilling[J].Journal of China Coal Society，2008，33（8）：898-902.
[13]章梦涛，王景琰.采场空气流动状况的数学模型和数值方法[J].煤炭学报，1983，8（3）：46-54.
ZHANG Mengtao，WANG Jingyan.The mathematical model and numerical method of the air flow in a stope[J].Journal of China Coal Society，1983，8（3）：46-54.
[14]董海波，童敏明，张丽苹，等.基于实测数据的工作面瓦斯分布场重建[J].采矿与安全工程学报，2012，29（1）：144-149.
DONG Haibo，TONG Minming，ZHANG Liping，et al.Reconstruction of gas distribution field in working face based on measured data[J].Journal of Mining and Safety Engineering，2012，29（1）：144-149.
[15]蒋曙光，张人伟.综放采场流场数学模型及数值计算[J].煤炭学报，1998，23（3）：258-261.
JIANG Shuguang，ZHANG Renwei.Mathematical model and numerical calculation of the flow field in fully mechanized caving field[J].Journal of China Coal Society，1998，23（3）：258-261.
[16]胡千庭，梁运培，刘见中.采空区瓦斯流动规律的CFD模拟[J].煤炭学报，2007，32（7）：719-723.
HU Qianting，LIANG Yunpei，LIU Jianzhong.CFD simulation of goaf gas flow patterns[J].Journal of China Coal Society，2007，32（7）：719-723.
[17]WIDIATMOJO A，SASAKI K，WIDODO N P，et al.Numerical simulation to evaluate gas diffusion of turbulent flow in mine ventilation system[J].International Journal of Mining Science and Technology，2013，23（3）：349-355.
[18]钱鸣高，许家林.覆岩采动裂隙分布的“O”形圈特征研究[J].煤炭学报，1998，23（5）：3-5.
QIAN Minggao，XV Jialin.Study on the O shape circle distribution characteristics of mining induced fractures in the overlaying strata[J].Journal of China Coal Society，1998，23（5）：3-5.
[19]李生舟.采动覆岩裂隙场演化及瓦斯运移规律研究[D].重庆：重庆大学，2012.
LI Shengzhou.Study on mining-induced fractures field evolution and gas migration rule[D].Chongqing：Chongqing University，2012.
[20]刘彦青，张浪，李伟，等.移动坐标下采空区流场及遗煤瓦斯浓度场分布[J].辽宁工程技术大学学报（自然科学版），2017，36（3）：243-248.
LIU Yanqing，ZHANG Lang，LI Wei，et al.Air-flowing field and gas-concentration field of the coal left in goaf based on the moving coordinate[J].Journal of Liaoning Technical University（Natural Science），2017，36（3）：243-248.
[21]羅振敏，郝苗，苏彬，等.采空区瓦斯运移规律实验及数值模拟[J].西安科技大学学报，2020，40（1）：31-39.
LUO Zhenmin，HAO Miao，SU Bin，et al.Experiments and numerical simulation research on gas migration in goaf[J].Journal of Xi’an University of Science and Technology，2020，40（1）：31-39. [22]林海飞.综放开采覆岩裂隙演化与卸压瓦斯运移规律及工程应用[D].西安：西安科技大学，2009.
LIN Haifei.Study on the law of mining-induced fracture evolution of overlying strata and relieved methane delivery and its engineering application in fully-mechanized top coal caving[D].Xi’an：Xi’an University of Science and Technology，2009.
[23]李春元，张勇，李佳，等.采空区瓦斯宏观流动通道的高位钻孔抽采技术[J].采矿与安全工程学报，2017，34（2）：391-397.
LI Chunyuan，ZHANG Yong，LI Jia，et al.Highly-located boreholes drainage technology of gas macroscopic flow channel in goaf[J].Journal of Mining & Safety Engineering，2017，34（2）：391-397.
[24]潘丰茂，王兆丰，陈滔，等.“Y”型通风工作面邻近层卸压瓦斯高效抽采技术[J].煤矿安全，2016，47（10）：52-55.
PAN Fengmao，WANG Zhaofeng，CHEN Tao，et al.Efficient drainage technology of released gas from adjacent layer at working face of Y ventilation[J].Safety in Coal Mines，2016，47（10）：52-55.
收稿日期：2020-08-16 責任编辑：杨泉林
基金项目：
国家自然科学基金资助项目（5197042023）
第一作者：薛俊华，男，江苏泰州人，教授，博导，E-mail：xuejunhua2003@163.com
通信作者：李延河，男，河南修武人，高级工程师，E-mail：z17775207879@126.com
关键词：安全科学与工程;瓦斯分布规律;FLUENT数值模型;上邻近保护层开采;风量分配比例
中图分类号：TD 712
文献标志码：A
文章编号：1672-9315（2021）01-0015-08
DOI：10.13800/j.cnki.xakjdxxb.2021.0103
Gas distribution law and control of Y-shaped ventilated goaf
in mining upper adjacent protective layer
——A case study of Fifth8-32010 working face in Pingmei No.6 Mine
XUE Junhua1，MA Qian1，LI Yanhe2，LI Hongbiao2，YUAN Zhandong2
（1.College of Safety Science and Engineering，Xi’an University of
Science and Technology，Xi’an 710054，China;
2.No.6 Mine of Pingdingshan Tian’an Coal Mining Co.，Ltd.，Pingdingshan 467000，China）
Abstract：In order to grasp the gas distribution law in the goaf of the Y-shaped ventilation working face during the mining of the upper adjacent protective layer，and control the gas emission in the goaf，this article first established
theoretical model
of the gas emission in the goaf for the Fifth8-32010 working face of Pingdingshan Coal Mine.Gas emission conditions were set according to the site situation and the inlet wind speed was set as follows：the ratio of the air volume in the engine lane to the air volume in the wind tunnel is 1∶3，1∶2，and 1∶1.5.A set of equations of gas-air mixture transition in goaf using FLUENT software were established.The numerical test was carried out of the gas distribution law in the goaf under the conditions of the three air volume ratios for the equations moving in the goaf.The results are compared and analyzed with the ones of the on-site industrial test.The air volume distribution suitable for the working face is calculated，and finally the drainage effect of the retaining borehole method is analyzed by numerical simulation.The results show that：the deeper the Y-shaped ventilated goaf，the higher the gas concentration，and the volume fraction at the tail side of the goaf is higher than that at the tail side of the roadway;the greater the air volume in the wind road，the lower the gas concentration in the working face.The smaller the air volume，the greater the concentration of gas in the working face and in the retained lane;the construction of the drainage borehole at the end of the retained lane to the goaf can achieve efficient and high-concentration gas drainage.Key words：safety science and engineering;gas distribution law;FLUENT numerical model;mining of upper adjacent protective layer;air volume distribution ratio 0 引 言
高瓦斯煤层在井工开采过程中，工作面一直面临着采空区和上隅角的瓦斯积聚问题[1-2]，尤其是近距离上保护层开采工作面，由于开采扰动后下邻近层卸压瓦斯窜入采空区，导致采空区和上隅角瓦斯浓度居高不下，严重妨碍着矿井安全生产。Y型通风系统可以消除上隅角瓦斯积聚，解决上隅角瓦斯易于超限的问题，是一条实现安全、高效、经济一体化生产的有效途径[3]。
虽然Y型通风系统解决了上隅角瓦斯超限问题，但同时改变了采空区瓦斯浓度场[4-5]，所以在施工过程中不能套用以前U型通风的瓦斯浓度场分布规律以及瓦斯治理技术[6-7]。只有摸清楚瓦斯浓度的分布规律，找出适合该工作面的系统风量及2条进风巷的风量分配，确定出高浓度瓦斯区域，在留巷内向此区域施工抽采钻孔才是解决采空区瓦斯浓度高的科学方法[8-9]。
为此，众多学者从理论研究、数值模拟、实验室试验和现场试验等方面进行了大量研究，获得了Y型通风采空区瓦斯浓度场分布规律[10]，偏Y型通风采空区瓦斯涌出规律[11]，开发了Y型通风风控瓦斯技术、留巷钻孔法抽采瓦斯技术等[12]。但是由于理论研究将现场复杂的条件做了理想简化及假设，一般不能很好的指导现场试验[13-14]。数值模拟试验在理论研究的基础上也做了一些条件假设，会与现场的试验有些差异，但是其适用性较强[15-17]。只有结合理论研究和数值模拟试验，才可以较好的指导现场工业性实验。
为了研究平顶山矿区煤层群特征下的上邻近保护层开采Y型通风采空区瓦斯浓度场分布规律及控制技术，文中基于理论研究，建立了采空区瓦斯涌出模型，并导入到FLUENT数值模型之中，研究了上邻近保护层Y型通风采空区瓦斯分布规律，并以进风巷风量分配为自变量，研究了其对瓦斯分布规律的影响，结合现场试验，确定最合适的进风巷风量分配，并在现场进行了验证，对比分析了导致数值模擬试验和现场试验结果之间差异的原因，最后利用留巷钻孔法对抽采效果进行考察，以期能为同类型煤层条件下的采空区瓦斯治理提供一定的参考依据。
1 矿井及工作面概况
平顶山天安煤业股份有限公司六矿主采煤层从上至下分别为丁5-6煤层、戊8煤层和戊9-10煤层，丁5-6煤层和戊8煤层平均层间距71 m，戊8煤层和戊9-10煤层平均层间距5 m，经煤层突出危险性鉴定后，决定将戊8煤层作为首采层。
戊8-32010工作面位于戊8煤层，平均煤厚2.2 m，设计可采长度2 300 m，倾向长度220 m，标高-633～-768 m，埋深873～938 m，通风系统结合高位瓦斯抽采巷及机巷切顶留巷布置，采用两进一回Y型通风方式，机巷、风巷进风，机巷留巷和机巷高抽巷回风，煤层为近水平煤层，遗煤平均厚度为0.2 m，可解吸瓦斯含量为4.3 m3/t，煤壁瓦斯初始涌出速率为0.003 56 m3/（m2·min）。工作面布置示意图如图1所示。
2 采空区瓦斯浓度分布规律及抽采数值模拟
2.1 条件假设
1）采空区为冒落矸石和岩块构成的带有空隙的多孔介质，并将多孔介质视为各向同性。
2）采空区气体不可压缩，不产生化学反应，其流动近似为稳态流动、等温过程。
3）忽略瓦斯的浮升效应，瓦斯充满整个采空区。
4）采空区瓦斯涌出源来自邻近层和遗煤。
5）遗煤和邻近层产生的瓦斯都是瞬间释放的，忽略解吸时间。
2.2 模型建立及网格划分
模拟对象为戊8-32010工作面，工作面长度220 m，采空区深度300 m，据现场考察，竖向裂隙带高度稳定在57 m，所以采空区高度设置为57 m，工作面宽度10 m，工作面高度4 m，风巷、机巷长度
20 m，宽度4 m，高度3.5 m，建成的模型如图2所示。
风巷、机巷和留巷网格间距为1 m×1 m，采空区模块网格尺寸为2 m×2 m，网格划分结果如图3所示。
2.3 采空区瓦斯涌出源模型
2.3.1 采空区空隙率分布模型
采空区空隙率与渗透率分布规律符合“O”形圈覆岩垮落形态，采空区空隙率与渗透率分布的数学表达如公式（1）和（2）[18]。
式中，n（x，y，z）为空隙率;L为工作面长度，m;k（x，y，z）为渗透率，m2;Dp为采空区冒落岩块平均粒径，取250 mm[19].
2.3.2 采空区遗煤瓦斯放散模型
采空区不同位置处的瓦斯涌出速率qc（y）可用公式组（3）求得[20]。
式中 a为瓦斯涌出初始强度，m3/min;b为瓦斯涌出衰减系数，min-1;y为遗煤位置到工作面的距离，m;Vx为工作面推进速度，m/min;Qt为单位体积遗煤可解吸瓦斯放散体积，m3/m3;T为解吸时间，取288 000 min（即200 d）.
通过上述方程组结合工作面遗煤放散条件，即可求出遗煤瓦斯涌出速率。
2.3.3 下邻近层瓦斯涌出模型
戊8-32010工作面开采时，戊9-10卸压瓦斯涌入工作面采空区，下邻近层瓦斯涌出量可由公式（4）求得，其中工作面煤壁绝对瓦斯涌出量参数，可由公式（5）求得，瓦斯排放率可由公式（6）求得，最终得出下邻近层的瓦斯涌出量[11]。
（6）
式中 Qs为戊9-10煤层量瓦斯涌出量，m3/min;mi 为戊9-10煤层厚度，3.1 m;Qb 为戊8-32010工作面绝对瓦斯涌出量，m3/min;hi 为邻近层层间距，m;其余薄煤层瓦斯涌出量忽略不计;ηs为戊9-10煤层的瓦斯排放率;L为戊8-32010工作面长度，m;h为开采厚度，2.2 m;S为开采层与邻近层层间距，13 m;q0为工作面煤壁刚暴露时单位面积煤壁的瓦斯涌出强度，m3/（m2·min）;V为工作面采煤机割煤时期平均牵引速度，m/min，实测值为10.08 m/min;Tc为工作面正常生产班期间，采煤机完成一次割煤周期所用时间，min，实测值为44 min。 2.4 采空区瓦斯运移方程组选择
选择的气体组分为甲烷-空气混合气体，选择湍流方程、连续性方程和动量守恒方程作为混合气体在采空区内的流动方程。假设混气扩散符合Fick
定律，其中混合气体在采空区中的流动方程
的统一表达式如式（7）所示，Fick扩散定律的表达式如式（8）所示[21]。
式中 ρ为气体密度，kg/m3;为通用变量;t为时间，s;
Γ为广义扩散系数;S为广义源项。
式中 r为极坐标半径，m;D为扩散系数，m2/s;c为瓦斯浓度，mol/m3。
2.5 边界条件及相关参数设置
为研究系统风量为2 830 m3/min时，两进风巷进风配比对工作面及留巷的瓦斯浓度的影响，共设置了2个进风入口和一个出口。2个进风入口，分别表示戊8工作面机巷和风巷，工作面机巷风量与机巷风量之比分别为：1∶3，1∶2和1∶1.5，一个出口设置为自然出流。
2.6 抽采条件设置
在留巷向里距离工作面切顶线10 m处增加抽放汇对采空区瓦斯进行抽采，抽采管直径为500 mm，钻孔高度为巷道底部2 m处。
3 数值模拟结果及分析
3.1 瓦斯浓度分布规律实验结果及分析
在总风量2 830 m3/min的基础上，对工作面机、风巷风量分配进行调整试验，分别对机巷、风巷风量比例1∶3，1∶2和1∶1.53种状态对工作面瓦斯质量分数的影响进行研究，结果如图4所示。
由图4可得，配风比例为1∶3时，通过工作面的风量较大，上下隅角压差增大，风巷端口与留巷压差也随之增大，导致工作面漏风严重，造成采空区防灭火难度加大;配风量为1∶1.5时，采空区和留巷处瓦斯浓度区域面积明显增大，基本可以看出增大风巷风量比例，工作面漏风加剧，减少风巷风量，采空区和留巷内瓦斯浓度会增大，而当采用配风比例1∶2时，一方面可以令工作面漏风情况得以控制，另一方面也可以缓解采空区和留巷内瓦斯积聚，所以配风比例为1∶2时更适合于戊8-32010工作面。
在Y型通风条件下，瓦斯向采空区深部不断积聚，且采空区尾部风巷侧瓦斯体积分数高于尾部留巷侧，由卸压钻孔法抽采卸压瓦斯机理可知[12]，采空区和覆岩裂隙区相连通，在留巷后部向采空区布置瓦斯抽采钻孔，可以持续高浓度抽采瓦斯资源。
3.2 瓦斯抽采实验结果及分析
采用配风比例为1∶2的通风模式建立戊8-32010工作面采空区抽采模型并布置尾抽钻孔，抽采模型如图5所示，抽采结果如图6所示。
由图6可以看出，利用尾抽钻孔抽采瓦斯后，采空区瓦斯体积分数明显降低，说明在留巷深部向采空区布置钻孔可以高效抽采瓦斯。
4 讨 论
本次数值实验结果以瓦斯质量分数显示，现场工业性实验结果以瓦斯体积分数显示。假设采空区瓦斯质量分数为cg，体积浓度为c′g，则质量分数和体积浓度存在以下关系[22]
（9）
式中 Ma为空气摩尔质量，kg/mol;Mg为瓦斯摩尔质量，kg/mol。
在总风量2 830 m3/min的基础上，对工作面机、风巷风量分配进行调整试验，分别对机巷、风巷风量比例
1∶3，1∶2和1∶1.5这3种状态时，对工作面3个测点和留巷内2个测点共5个测点的瓦斯浓度进行对比分析，测点布置如图7所示，现场试验测点与对应测试位置见表1。
先将现场试验和数值模拟实验5个测点位置的瓦斯体积分数统计、并绘图，瓦斯体积浓度对比结果如图8所示。
由图8（a）可得，机、风巷风量配比1∶3时工作面瓦斯浓度在0.16%～0.3%之间，留巷段瓦斯浓度在0.32%～0.69%之间;配比为1∶1.5时工作面瓦斯浓度在0.24%～0.5%之间，留巷段瓦斯浓度在0.51%～0.64%之間。由图8（b）可得，机、风巷风量配比1∶3时工作面瓦斯浓度在0.02%～
0.05%之间，留巷段瓦斯浓度在0.05%～0.3%之间;配比为1∶1.5时工作面瓦斯浓度在0.05%～0.07%之间，留巷段瓦斯浓度在0.04%～0.28%之间。
综上可得，风巷风量比例越小，工作面瓦斯浓度越高，留巷段瓦斯浓度越高;风巷风量比例越大，工作面瓦斯浓度越低，现场试验和数值模拟试验采空区瓦斯浓度的变化趋势基本一致。但工作面处和留巷内瓦斯浓度值差异较大，造成这一差异的主要原因是数值模拟中无法控制工作面向采空区内漏风和采空区向留巷内漏风的过程，后续笔者会就这一问题做进一步研究。
最后在留巷段向里距离工作面切顶线10 m处施工了瓦斯尾抽钻孔，抽采了采空区瓦斯，在此期间，工作面和采空区瓦斯都处于安全值以内[23-24]，相比于上个相邻U型工作面推进速度，本工作面推进速度明显提高，如图9为日产量与抽采浓度之间的关系图。
从图9可以看出，利用尾抽抽采采空区尾部瓦斯，瓦斯浓度都在50%以上，说明在尾部布置瓦斯抽采效果较好。在抽采效果考察后期，日产量有短时间的相对减小，这是由于此时煤层上覆岩层破断垮落，裂隙带初步形成，将尾抽为主的抽采方式转化为高抽巷裂隙带瓦斯为主抽采方式，一部分瓦斯通过裂隙带进行抽采，暂时降低了日产量。
5 结 论
1）在采空区瓦斯涌出和下邻近层瓦斯涌出强度一定时，近距离上保护层Y型通风的瓦斯易于积聚在采空区尾部，越靠近尾部，瓦斯体积分数增加速率越大，采空区尾部留巷侧瓦斯体积分数低于尾部风巷侧，在留巷尾部一定位置处布置瓦斯抽采钻孔能高效抽采瓦斯，并且此浓度场与单煤层Y型通风瓦斯浓度场分布规律基本一致。
2）增大风巷风量比例，会导致采空区瓦斯大量涌入留巷;减小风巷风量比例，会导致工作面瓦斯增加，因此改变2条进风巷配风量，工作面和留巷的瓦斯浓度会随之变化。 参考文献（References）：
[1] 姜福兴，孔令海，刘春刚.特厚煤层综放采场瓦斯运移规律[J].煤炭学报，2011，36（3）：407-411.
JIANG Fuxing，KONG Linghai，LIU Chungang.Gas emission laws of fully mechanized sublevel caving mining in extra thick coal seam[J].Journal of China Coal Society，2011，36（3）：407-411.
[2]BALUSU R，DEGUCHI G，HOLLAND R，et al.Goaf gas flow mechanics and development of gas and Sponcom control strategies at a highly gassy mine[J].Coal and Safety，2015，20：35-45.
[3]贾廷贵，娄和壮，闻文.不同配风条件下“Y”型通风采场瓦斯涌出规律[J].辽宁工程技术大学学报（自然科学版），2019，38（1）：14-18.
JIA Tinggui，LOU Hezhuang，WEN wen.Gas migration law of Y-type ventilation stope under different air distribution conditions[J].Journal of Liaoning Technical University（Natural Science），2019，38（1）：14-18.
[4]袁亮，薛俊华.低透气性煤层群无煤柱煤与瓦斯共采关键技术[J].煤炭科学技术，2013，41（1）：5-11.
YUAN Liang，XUE Junhua.Key technology of pillarless coal and gas simultaneous mining in low permeability seam group[J].Coal Science and Technology，2013，41（1）：5-11.
[5]戚良锋.Y型通风方式下采空区瓦斯流场数值模拟研究[D].淮南：安徽理工大学，2009.
QI Liangfeng.Research on gas field in gob with Y-type ventilation system by numerical simulation[D].Huainan：
Anhui University of Science and Technology，2009.
[6]陈胜利.综放工作面瓦斯涌出特征实测及分析[J].科学技术与工程，2012，27（12）：6892-6896.
CHEN Shengli.Measurement and analysis for gas emission characters in fully mechanized caving face[J].Science Technology and Engineering，2012，27（12）：6892-6896.
[7]徐玉胜，李春元，张勇，等.大采高采场瓦斯治理模型分析及通风系统优化研究[J].煤炭科学技术，2019，47（3）：142-149.
XU Yusheng，LI Chunyuan，ZHANG Yong，et al.Study on gas control model analysis and ventilation system optimization in large mining height mines[J].Coal Science and Technology，2019，47（3）：142-149.
[8]周一力.“Y+高抽巷”工作面采空区瓦斯与氧气浓度场分布规律及其在灾害防治中的应用[D].徐州：中国矿业大学，2019.
ZHOU Yili.Distribution of gas and oxygen concentration field in the goaf of “Y+High Pumping Alley”working face and its applications in disaster prevention[D].Xuzhou：China University of Miing and Technology，2019.
[9]常海雷，刘彦青，张浪，等.“U+I”型通风形式下综放工作面采空区遗煤及支架上方煤壁瓦斯放散与运移规律研究[J].矿业安全与环保，2017，44（2）：39-44.
CHANG Hailei，LIU Yanqing，ZHANG Lang，et al.Study on gas emission and migration law of residual coal in gob and coal wall above support in fully mechanized caving face with “U+I” shape ventilation system[J].Mining Safety & Environmental Protection，2017，44（2）：39-44.
[10]張欢，赵洪宝，李晓白，等.分段留巷Y型通风采场瓦斯与流场分布三维实测与重构[J].中国安全生产科学技术，2018，14（7）：107-114.
ZHANG Huan，ZHAO Hongbao，LI Xiaobai，et al.Three-dimensional measurement and reconstruction of flow field and gas distribution for mining face and goaf of Y type ventilation[J].Journal of Safety Science and Technology，2018，14（7）：107-114. [11]赵灿，张浪，刘彦青.偏Y型通风下采空区瓦斯涌出规律及超限治理研究[J].煤炭科学技术，2019，47（4）：127-133.
ZHAO Can，ZHANG Lang，LIU Yanqing.Study on gob gas emission laws and control of gas exceeding limit under inclined type Y-ventilation[J].Coal Science and Technology，2019，47（4）：127-133.
[12]袁亮.留巷钻孔法煤与瓦斯共采技术[J].煤炭学报，2008，33（8）：898-902.
YUAN Liang.The technique of coal mining and gas extraction by roadway retaining and borehole drilling[J].Journal of China Coal Society，2008，33（8）：898-902.
[13]章梦涛，王景琰.采场空气流动状况的数学模型和数值方法[J].煤炭学报，1983，8（3）：46-54.
ZHANG Mengtao，WANG Jingyan.The mathematical model and numerical method of the air flow in a stope[J].Journal of China Coal Society，1983，8（3）：46-54.
[14]董海波，童敏明，张丽苹，等.基于实测数据的工作面瓦斯分布场重建[J].采矿与安全工程学报，2012，29（1）：144-149.
DONG Haibo，TONG Minming，ZHANG Liping，et al.Reconstruction of gas distribution field in working face based on measured data[J].Journal of Mining and Safety Engineering，2012，29（1）：144-149.
[15]蒋曙光，张人伟.综放采场流场数学模型及数值计算[J].煤炭学报，1998，23（3）：258-261.
JIANG Shuguang，ZHANG Renwei.Mathematical model and numerical calculation of the flow field in fully mechanized caving field[J].Journal of China Coal Society，1998，23（3）：258-261.
[16]胡千庭，梁运培，刘见中.采空区瓦斯流动规律的CFD模拟[J].煤炭学报，2007，32（7）：719-723.
HU Qianting，LIANG Yunpei，LIU Jianzhong.CFD simulation of goaf gas flow patterns[J].Journal of China Coal Society，2007，32（7）：719-723.
[17]WIDIATMOJO A，SASAKI K，WIDODO N P，et al.Numerical simulation to evaluate gas diffusion of turbulent flow in mine ventilation system[J].International Journal of Mining Science and Technology，2013，23（3）：349-355.
[18]钱鸣高，许家林.覆岩采动裂隙分布的“O”形圈特征研究[J].煤炭学报，1998，23（5）：3-5.
QIAN Minggao，XV Jialin.Study on the O shape circle distribution characteristics of mining induced fractures in the overlaying strata[J].Journal of China Coal Society，1998，23（5）：3-5.
[19]李生舟.采动覆岩裂隙场演化及瓦斯运移规律研究[D].重庆：重庆大学，2012.
LI Shengzhou.Study on mining-induced fractures field evolution and gas migration rule[D].Chongqing：Chongqing University，2012.
[20]刘彦青，张浪，李伟，等.移动坐标下采空区流场及遗煤瓦斯浓度场分布[J].辽宁工程技术大学学报（自然科学版），2017，36（3）：243-248.
LIU Yanqing，ZHANG Lang，LI Wei，et al.Air-flowing field and gas-concentration field of the coal left in goaf based on the moving coordinate[J].Journal of Liaoning Technical University（Natural Science），2017，36（3）：243-248.
[21]羅振敏，郝苗，苏彬，等.采空区瓦斯运移规律实验及数值模拟[J].西安科技大学学报，2020，40（1）：31-39.
LUO Zhenmin，HAO Miao，SU Bin，et al.Experiments and numerical simulation research on gas migration in goaf[J].Journal of Xi’an University of Science and Technology，2020，40（1）：31-39. [22]林海飞.综放开采覆岩裂隙演化与卸压瓦斯运移规律及工程应用[D].西安：西安科技大学，2009.
LIN Haifei.Study on the law of mining-induced fracture evolution of overlying strata and relieved methane delivery and its engineering application in fully-mechanized top coal caving[D].Xi’an：Xi’an University of Science and Technology，2009.
[23]李春元，张勇，李佳，等.采空区瓦斯宏观流动通道的高位钻孔抽采技术[J].采矿与安全工程学报，2017，34（2）：391-397.
LI Chunyuan，ZHANG Yong，LI Jia，et al.Highly-located boreholes drainage technology of gas macroscopic flow channel in goaf[J].Journal of Mining & Safety Engineering，2017，34（2）：391-397.
[24]潘丰茂，王兆丰，陈滔，等.“Y”型通风工作面邻近层卸压瓦斯高效抽采技术[J].煤矿安全，2016，47（10）：52-55.
PAN Fengmao，WANG Zhaofeng，CHEN Tao，et al.Efficient drainage technology of released gas from adjacent layer at working face of Y ventilation[J].Safety in Coal Mines，2016，47（10）：52-55.
收稿日期：2020-08-16 責任编辑：杨泉林
基金项目：
国家自然科学基金资助项目（5197042023）
第一作者：薛俊华，男，江苏泰州人，教授，博导，E-mail：xuejunhua2003@163.com
通信作者：李延河，男，河南修武人，高级工程师，E-mail：z17775207879@126.com