白龙山煤矿10201工作面合理挡风帘长度确定

崔传发

(郑州信息科技职业学院 建筑工程学院，河南 郑州 450000)

0 引言

采空区漏风是工作面安全回采的重要隐患之一[1]。采空区漏风不仅会导致工作面有效风量降低、上隅角瓦斯浓度增加、温度上升，同时也是引发采空区煤自燃的必要条件，严重危及井下作业人员的生命安全[2-3]。设置挡风帘具有成本低、易操作等特点，成为最常见的防治采空区漏风的方法。王小龙[4]、邢纪伟等[5]针对上隅角瓦斯超限问题，提出通过设置挡风帘来减小工作面漏风，从而降低上隅角瓦斯浓度。陈亚忠[6]通过现场实测吊挂挡风帘前后的工作面风量和瓦斯体积分数，发现挂帘后二者降幅分别为88.5%，76.2%。王云龙等[7]利用有限元法对采空区进行数值模拟，得出设置挡风帘可减小采空区自燃危险区域的宽度。高建良等[8]研究了挡风帘长度对采空区瓦斯和自燃“三带”分布的影响，结果表明，随着挡风帘长度增加，采空区瓦斯浓度增加，采空区进风侧自燃氧化带逐渐向工作面移动，采空区中部自燃氧化带的宽度明显减小。张学博等[9]研究发现挡风帘可使“U+L”型通风工作面采空区瓦斯浓度上升，采煤机附近瓦斯浓度降低。

白龙山煤矿10201工作面漏风严重且上隅角瓦斯浓度偏高，本文以该工作面为研究背景，利用Fluent软件模拟不同挡风帘长度下采空区流场，分析挡风帘长度对工作面漏风、工作面瓦斯浓度、采空区瓦斯分布及自燃氧化带范围的影响规律，从而确定工作面合理挡风帘长度，并对挂帘前后工作面风量及上隅角瓦斯和一氧化碳浓度进行现场考察。

1 采空区流场数学模型

1.1 基本假设

模型中流体视为不可压缩流体；采空区视为各向同性的多孔介质；采空区气体流动符合达西定律，扩散运动遵循菲克定律；忽视采空区温度和非定常因素对流场的影响。

1.2 采空区气体流动控制方程

(1)质量守恒方程。当不考虑流体密度的变化时，多孔介质渗流场质量守恒方程为

(1)

式中：ρ为采空区气体密度，kg/m3；vi为i(i=x，y，z)方向渗流场达西速度分量，m/s；Sm为质量源项，kg/(m3·s)。

(2)动量守恒方程。定常流动动量守恒方程在一个惯性参考系(没有加速度)内可表示为

(2)

式中：μ为流体的动力黏度，N·s/m2；p为采空区内静压力，Pa；Si为i方向动量源项，N/m3。

(3)组分输运方程。对于定常流场，采空区气体组分输运方程为

(3)

式中：ωj为组分j的质量分数；Dj为组分j的扩散系数，m2/s；Fj为组分j的生产率，1/s。

1.3 多孔介质模型

在Fluent中，使用以经验假设为主的流动阻力即动量源项来表征多孔介质属性，对于各向同性多孔介质，动量源项包括黏性损失项(式(4)中等号右边第一项)和惯性损失项(式(4)中等号右边第二项)。

(4)

(5)

(6)

式中：α为渗透性系数，m2；C为惯性阻力系数，1/m；v为渗流场达西速度，m/s；dp为介质调和平均粒径，m；n为孔隙率。

2 采空区物理模型

2.1 基本假设

在数值模拟中忽略工作面的采煤机、刮板输送机等设备，对工作面、采空区进行以下简化：忽略矿井周期来压等特殊情况，只考虑通风、高抽巷抽采及挡风帘对采空区漏风规律及自燃“三带”的影响；在不影响工程精度及基本规律的前提下，为了建模方便，将采空区简化为梯台体，将进风巷、回风巷、高抽巷空间视为长方体；将液压支架视为只有惯性损失项的多孔介质；假设已进入采空区部分的高抽巷未发生坍塌；认为设置挡风帘后，挡风帘部分不存在漏风。

2.2 物理模型

白龙山煤矿10201工作面开采的C7+8号煤层属于高瓦斯自燃煤层，工作面内煤层平均有益厚度为4.65 m，夹矸平均厚度为0.57 m，煤层厚度稳定，可采性指数为1。工作面倾向平均长219.2 m，设计风量1 200 m3/min，在回风巷上方25 m内错30 m位置施工有3 m×3 m高抽巷用于工作面卸压瓦斯治理，结合埋管抽采治理上隅角瓦斯，高抽巷抽采流量为48 m3/min，埋管抽采流量为36 m3/min。根据采煤工作面实际情况，利用SCDM软件建立回采工作面采空区物理模型，参数见表1，并利用ICEM对模型进行结构化网格划分，如图1所示。

表1 模型参数

图1 采空区物理模型及网格划分

2.3 边界条件及参数设置

(1)边界类型设置。进风巷入口边界条件为速度入口，风速设为1.33 m/s；回风巷、高抽巷和埋管抽采出口边界条件为自由出口；工作面与采空区交界面分为2个部分，设置挡风帘部分为固壁边界，未设置挡风帘部分为内部边界；其余边界设为默认固壁边界。

(2)多孔介质参数设置。根据采空区“O”形圈理论[10]，可近似拟合出采空区垮落碎胀系数[11]：

KP=KP,min+(KP,max-KP,min)×

exp(-a1l1(1-exp(-ξ1)a0l0))

(7)

式中：KP,min为压实后的垮落碎胀系数，取1.15；KP,max为初始垮落碎胀系数，取1.5；a0，a1分别为距离固壁和工作面的衰减率，分别取0.268，0.036 8；l0，l1分别为采空区任意点距固壁和工作面边界的距离，m；ξ1为控制“O”形圈模型分布形态的调整数(应小于1)，取0.233。

采空区孔隙率为

(8)

(3)源项设置。工作面瓦斯源项主要包括工作面瓦斯涌出和采空区遗煤瓦斯涌出，分别为14.7，24.9 m3/min。假设工作面和采空区瓦斯为均匀涌出，采空区3 m以上岩层无瓦斯涌出和氧气消耗，瓦斯密度取0.716 7 kg/m3，则工作面瓦斯质量源项为3.326×10-5kg/(m3·s)，采空区瓦斯质量源项为2.048×10-6kg/(m3·s)。

在不考虑温度变化时，氧气源项即遗煤耗氧速率。

ε=-ρO2λ(c-c0)

(9)

式中：ε为遗煤耗氧速率，kg/(m3·s)；ρO2为氧气密度，取1.429 kg/m3；λ为与温度和煤样有关的实验常数，取7.436×10-6s-1；c为氧气体积分数；c0为遗煤氧化所需最低氧气体积分数，取2%。

3 数值模拟结果分析

3.1 工作面风量及瓦斯浓度变化

挡风帘长度对工作面风量的影响如图2所示，其中ΔQ为漏风量，“+”代表漏入工作面，“-”代表漏出工作面。可看出在工作面距进风巷0～170 m范围内工作面风量逐渐减小，在170 m至工作面回风口范围内，采空区向工作面涌出瓦斯风流，工作面风量逐渐增大；在工作面距进风巷0～80 m范围内，工作面风量随挡风帘长度增加而逐渐增加；在工作面距进风巷80 m之后，挡风帘长度对工作面风量的影响较小；在工作面距进风巷170 m处工作面风量最小，不同挡风帘长度下最大漏风量均为125 m3/min，进回风巷风量均相差83 m3/min。

图2 不同挡风帘长度下工作面风量变化曲线

挡风帘长度对工作面瓦斯浓度的影响如图3所示。可看出工作面进风侧设置挡风帘可有效降低工作面瓦斯浓度；在工作面距进风巷0～190 m范围内，随着挡风帘长度增加，工作面瓦斯浓度逐渐下降；在工作面距进风巷190～220 m范围内，是否设置挡风帘对该区域瓦斯浓度影响较大，而挡风帘长度对瓦斯浓度影响较小；当挡风帘长度为0即不设置挡风帘时，工作面上隅角最大瓦斯体积分数为0.42%，设置5 m挡风帘后上隅角最大瓦斯体积分数降至0.37%，降幅达12%左右，当挡风帘长度增加至10，15，20 m时，挡风帘长度对上隅角最大瓦斯浓度影响不大。

图3 不同挡风帘长度下工作面瓦斯浓度变化曲线

3.2 采空区瓦斯分布

不同挡风帘长度下采空区瓦斯分布云图如图4所示。从走向上来看，采空区瓦斯浓度整体上随采空区深度增加逐渐上升，由于高抽巷和上隅角埋管的抽采作用，在采空区回风侧存在一个瓦斯聚集区，瓦斯体积分数达30%左右；从倾向上来看，采空区回风侧瓦斯浓度整体上大于采空区进风侧；挡风帘长度对采空区瓦斯分布规律的影响较小。

(a)挡风帘长度0 m

不同挡风帘长度下采空区进风侧和回风侧瓦斯浓度变化曲线如图5所示。由采空区进风侧瓦斯浓度变化曲线可看出，在采空区深度为0～140 m范围内，采空区进风侧瓦斯浓度随挡风帘长度增加逐渐升高，这是由于挡风帘可有效降低进风侧漏风量，漏风量越小对采空区瓦斯的稀释作用就越弱；在采空区深度为140～220 m范围内，采空区进风侧瓦斯浓度随挡风帘长度增加呈先升高后降低的趋势，当挡风帘长度为10 m时，该区域瓦斯体积分数最高可达23%左右。由采空区回风侧瓦斯浓度变化曲线可看出，在采空区深度为0～10 m范围内，采空区回风侧瓦斯浓度因上隅角埋管抽采作用而快速降低，受挡风帘长度影响较小；在采空区深度为10～130 m范围内，挡风帘长度为0时，采空区回风侧瓦斯浓度最高，而当挡风帘长度分别为5，10，15，20 m时，采空区回风侧瓦斯浓度随挡风帘长度增加呈先升高后降低的趋势，挡风帘长度为20 m时，采空区回风侧瓦斯浓度最小；在采空区深度为130～170 m范围内，采空区回风侧瓦斯浓度与挡风帘长度并无明显关系；在采空区深度为170～220 m范围内，采空区回风侧瓦斯浓度随挡风帘长度增加呈先升高后降低的趋势，当挡风帘长度为10 m时采空区回风侧瓦斯浓度最大，当挡风帘长度为0时采空区回风侧瓦斯浓度最小。

图5 不同挡风帘长度下采空区两侧瓦斯浓度变化曲线

3.3 采空区自燃氧化带

采空区自燃“三带”划分指标有采空区内部漏风风速、氧气体积分数和测点升温特征3种，其中氧气体积分数指标因其数据准确、可操作性强、代价低等优点而被广泛使用，因此本文选择氧气体积分数指标来划分采空区自燃“三带”，即氧气体积分数大于18%为散热带，氧气体积分数在10%～18%之间为氧化带，氧气体积分低于10%为窒息带[12-13]。

不同挡风帘长度下采空区氧气分布云图如图6所示。可看出采空区自燃氧化带即氧气体积分数为10%～18%区域呈条带状；采空区进风侧和回风侧自燃氧化带向采空区深部凸起，且采空区回风侧自燃氧化带宽度大于进风侧。

(a)挡风帘长度0 m

不同挡风帘长度下采空区进风侧和回风侧氧气浓度变化曲线如图7所示。根据氧气体积分数指标，可划分出不同挡风帘长度下采空区进风侧和回风侧自燃氧化带范围，见表2。

图7 不同挡风帘长度下采空区两侧氧气浓度变化曲线

表2 不同挡风帘长度下采空区自燃氧化带划分

从表2可看出，采空区回风侧自燃氧化带宽度整体大于采空区进风侧，这是由于高抽巷和上隅角埋管抽采的扰动效应使得采空区回风侧漏风流局部加速，自燃氧化带宽度有所增加；采空区回风侧进入自燃氧化带的时间较采空区进风侧早；在进风侧设置挡风帘使得采空区进风侧自燃氧化带宽度变大，而采空区回风侧自燃氧化带宽度减小；当挡风帘长度分别为5，10，15，20 m时，采空区进回风侧自燃氧化带宽度随挡风帘长度增加呈先增大后减小的趋势，且均在挡风帘长度为15 m时自燃氧化带宽度最大；随着挡风帘长度增加，采空区进风侧自燃氧化带逐渐向工作面靠近，即进入自燃氧化带时间变早，而挡风帘长度对采空区回风侧自燃氧化带的影响主要体现在进入自燃氧化带的位置，但并无明显规律。

4 现场应用效果

根据数值模拟结果发现，挡风帘对工作面风量和瓦斯浓度影响较大，对采空区瓦斯浓度和自燃氧化带的影响较小。考虑到挡风帘越长则实施难度越大，在保证有效降低上隅角瓦斯浓度的前提下，尽量减小挡风帘长度，确定在白龙山煤矿10201工作面进风侧布置5 m挡风帘，略大于进风巷宽度，如图8所示。

图8 白龙山煤矿10201工作面挡风帘布置

为验证白龙山煤矿10201工作面进风侧挡风帘的应用效果，在挡风帘布置前后，对工作面进行风量测定和上隅角瓦斯及一氧化碳浓度监测。

风量测定共布置10个测点，测定结果如图9所示。可看出在布置5 m挡风帘后，工作面距进风巷0～86 m范围内风量有明显增加，工作面距进风巷约5 m处风量增加60 m3/min；在工作面距进风巷86 m之后风量变化不大；在工作面距进风巷167 m处工作面风量最小。

图9 挂帘前后风量对比

从2018年10月5日至11月4日，对工作面上隅角瓦斯及一氧化碳浓度进行为期30 d的监测，监测值均取早中晚三班最大值，结果如图10所示。可看出在设置挡风帘前，瓦斯体积分数平均值为0.602%，一氧化碳体积分数平均值为2.46%；于2018年10月15日布置挡风帘后，瓦斯体积分数平均值为0.521%，降幅达13.5%，一氧化碳体积分数平均值为2.26%，降幅为8.1%。

图10 上隅角瓦斯和一氧化碳体积分数变化曲线

5 结论

(1)在工作面距进风巷0～80 m范围内，工作面风量随挡风帘长度增加而逐渐增加，而在工作面距进风巷80 m之后，挡风帘长度对工作面风量的影响较小；当不设置挡风帘时，工作面上隅角最大瓦斯体积分数为0.42%，而设置5 m挡风帘后上隅角最大瓦斯体积分数降至0.37%，降幅达12%左右，当挡风帘长度增加至10，15，20 m时，挡风帘长度对上隅角最大瓦斯浓度影响不大。

(2)挡风帘可降低采空区回风侧浅部和中部(采空区深度小于130 m)的瓦斯浓度，进而减小采空区瓦斯向工作面涌出，而对于采空区进风侧和回风侧深部(采空区深度大于170 m)区域，挡风帘会使瓦斯浓度上升。

(3)在进风侧设置挡风帘会使采空区进风侧自燃氧化带宽度变大，而采空区回风侧自燃氧化带宽度减小；随着挡风帘长度增加，采空区进风侧自燃氧化带逐渐向工作面靠近，即进入自燃氧化带时间变早，而挡风帘长度对采空区回风侧自燃氧化带的影响主要体现在进入自燃氧化带的位置，但并无明显规律。

(4)在白龙山煤矿10201工作面进风侧布置5 m挡风帘，增加了工作面有效风量，瓦斯体积分数平均值为0.521%，降幅达13.5%，一氧化碳体积分数平均值为2.26%，降幅为8.1%，降低了上隅角瓦斯超限和采空区自然发火的危险性。