一种矿井工作面瓦斯涌出量预测的方法与流程

本发明属于煤矿矿井瓦斯治理技术领域，尤其涉及一种矿井工作面瓦斯涌出量预测的方法。

背景技术

工作面瓦斯预测是矿井瓦斯预测的基础。目前矿井瓦斯涌出量预测大多采用分源预测法和统计分析法。其中分源预测法主要针对薄及中厚煤层不分层开采，对于厚煤层开采的瓦斯涌出量预测采用瓦斯涌出系数进行计算，由于具体开采条件的差异，相关的系数往往需要进行实际测定，在没有进行测定前会造成预测数据的不精确。而统计预测法则要求被预测的回采工作面在开采方法、煤层赋存条件、瓦斯地质条件与样本工作面相同或相似。在厚煤层开采时，尤其是用放顶煤开采时，煤体放出的不均匀性、不同煤厚和顶板矿山压力对煤体破断影响的差异性、以及不同开采条件回采率不同等原因，目前的预测方法对放顶煤开采时的预测普遍存在困难。

技术实现要素：

本发明为了解决现有技术中的不足之处，提供一种可实时监控锚杆受力状况、提高安全支护监测、及时提醒工人排出安全隐患的矿井工作面瓦斯涌出量预测的方法。

为解决上述技术问题，本发明采用如下技术方案：一种矿井工作面瓦斯涌出量预测的方法，包括以下步骤，

（1）、根据矿山压力对工作面的作用效果，把工作面前方煤体由后向前依次划分为第一初始压缩及塑性变形区a、第一煤体的强烈压缩及破坏区b和第一顶板回转作用区c；

（2）、矿井下的区域运输巷平行于区域回风巷，区域运输巷位于区域回风巷的左侧，区域运输巷的右侧帮由左向右依次划分为第二初始压缩及塑性变形区a1、第二煤体的强烈压缩及破坏区b1和第二顶板回转作用区c1，区域回风巷的左侧帮由右向左依次划分为第三初始压缩及塑性变形区a2、第三煤体的强烈压缩及破坏区b2和第三顶板回转作用区c2；

（3）、对步骤（1）中工作面前方煤体瓦斯绝对涌出量qa；

（4）、对步骤（2）中区域运输巷的右侧帮和区域回风巷的左侧帮煤体的瓦斯绝对涌出量q1+2；

（5）、最后得到矿井工作面处瓦斯涌出量q=qa+q1+2。

步骤（3）中qa的确定过程为：

第一初始压缩及塑性变形区a与第一煤体的强烈压缩及破坏区b相邻，其分界线是由煤体各层位破坏支承压力点连成的，其截距及该煤壁至分界线与煤层底板界面交点的水平距离x0可以用下式表示：

其中x0——煤壁至曲线与煤层底板界面交点的水平距离，m；

m——煤层厚度，m；

λ——侧压系数，，μ为煤体的泊松比；

φ——煤体的内摩擦角，°；

k——支承压力峰值的应力集中系数；

h——煤层采深，m；

γ——上覆岩层平均容重，；

c——煤体的内聚力，mpa；

px——工作面支架对煤壁的支护反力形成的支护强度，mpa；；

ft——支架工作阻力，mn；

α——支架立柱倾角，°；

b——支架宽度，m。

第一初始压缩及塑性变形区a与第一煤体的强烈压缩及破坏区b的分界线的曲线方程式可以表示为：

x=0.0784h²+0.376h+x0

则由此可以计算破碎区宽度lp为：

lp=0.0784m²+0.376m+x0

1）第一初始压缩及塑性变形区a煤体体积的确定

由于第一煤体的强烈压缩及破坏区b为分界线下方x0至lp所围成的部分；设工作面长度为d米，煤层厚度、即采高为m米，则第一初始压缩及塑性变形区a的煤体体积va为：

2）第一初始压缩及塑性变形区a的煤体瓦斯涌出量

设工作面截深为e米，采煤机割煤速度为v米/分钟，根据现场统计，工作面割一刀煤的时间l/v一般不会小于15分钟，当15≤l/v≤30时，破碎煤体的瓦斯涌出量为总瓦斯含量的60%；当31≤l/v≤60时，破碎煤体的瓦斯涌出量为总瓦斯含量的85%；61≤l/v时，破碎煤体的瓦斯涌出量为总瓦斯含量的95%；

由此可以计算当工作面截深为e米、采煤机割煤速度为v米/分钟、原始相对瓦斯含量为q立方米/吨时，割煤后a区瓦斯绝对涌出量：

式中：

k——瓦斯释放均匀系数，根据瓦斯赋存的均匀程度可取1.0-1.5；

——瓦斯释放程度系数，当15≤l/v≤30时，取0.6；当31≤l/v≤60时，取0.85；61≤l/v时，取0.95。

步骤（4）中区域运输巷的右侧帮和区域回风巷的左侧帮煤体的瓦斯绝对涌出量q的具体确定过程为：

第二初始压缩及塑性变形区a1和第二煤体的强烈压缩及破坏区b1之间的分界线与第三初始压缩及塑性变形区a2和第三煤体的强烈压缩及破坏区b2之间的分界线左右对称布置，

其截距及该煤壁至分界线与煤层底板界面交点的水平距离x1可以用下式表示

px1——煤壁前方的超前压力

根据分界线的曲线方程式：

x=0.0784h²+0.376h+x1

则由此可以计算破碎区宽度lp1为：

lp1=0.0784m²+0.376m+x1

超前压力影响下的第二初始压缩及塑性变形区a1和第三初始压缩及塑性变形区a2的破碎情况与沿工作面方向是相似的，因为超前压力影响有一定的范围，取破碎的范围近似为20m；

则超前压力影响下的两帮的破碎区的体积

原始相对瓦斯含量为q立方米/吨时，在超前压力的影响下，两帮煤壁的瓦斯涌出量按50%来计算，则煤壁两帮的瓦斯涌出量q1+2为：。

采用上述技术方案，本发明根据矿山压力对工作面的作用效果，把工作面前方、区域运输巷的右侧帮及区域回风巷左侧帮均划分为初始压缩及塑性变形区、煤体的强烈压缩及破坏区和顶板回转作用区。而煤体的强烈压缩及破坏区和顶板回转作用区处于应力增高区，由于应力的升高致使煤体内压力升高，从而使该部分煤体中部分游离态的瓦斯会变成吸附态瓦斯，不利于瓦斯的释放。因为超前压力的影响，使初始压缩及塑性变形区与工作面前方两帮煤壁更加的破碎，从而释放出一定量的瓦斯。所以工作面瓦斯的释放可以看做是这两部分释放的，因此算出初始压缩及塑性变形区的瓦斯量与区域运输巷的右侧帮及区域回风巷左侧帮的释放量就可以近似的表达工作面的瓦斯含量。

综上所述，本发明通过将工作面前方煤层划分为不同的几个区域，通过计算破碎区域区域的瓦斯涌出含量，从而得到工作面的瓦斯涌出量。通过计算得到的数据，用来指导矿井采取更加经济合适的措施预防瓦斯，从而提高矿井安全生产水平。

附图说明

图1是工作面前方煤体沿平行于巷道长度方向垂直截面示意图；

图2是沿巷道宽度方向垂直截面示意图。

具体实施方式

本发明的一种矿井工作面瓦斯涌出量预测的方法，包括以下步骤，

（1）、如图1所示，根据矿山压力对工作面的作用效果，把工作面前方煤体由后向前依次划分为第一初始压缩及塑性变形区a、第一煤体的强烈压缩及破坏区b和第一顶板回转作用区c；

（2）、如图2所示，矿井下的区域运输巷1平行于区域回风巷2，区域运输巷1位于区域回风巷2的左侧，区域运输巷1的右侧帮由左向右依次划分为第二初始压缩及塑性变形区a1、第二煤体的强烈压缩及破坏区b1和第二顶板回转作用区c1，区域回风巷2的左侧帮由右向左依次划分为第三初始压缩及塑性变形区a2、第三煤体的强烈压缩及破坏区b2和第三顶板回转作用区c2；

（3）、对步骤（1）中工作面前方煤体瓦斯绝对涌出量qa；

（4）、对步骤（2）中区域运输巷1的右侧帮和区域回风巷2的左侧帮煤体的瓦斯绝对涌出量q1+2；

（5）、最后得到矿井工作面处瓦斯涌出量q=qa+q1+2。

步骤（3）中qa的确定过程为：

如图1所示，第一初始压缩及塑性变形区a与第一煤体的强烈压缩及破坏区b相邻，其分界线是由煤体各层位破坏支承压力点连成的，其截距及该煤壁至分界线与煤层底板界面交点的水平距离x0可以用下式表示：

其中x0——煤壁至曲线与煤层底板界面交点的水平距离，m；

m——煤层厚度，m；

λ——侧压系数，，μ为煤体的泊松比；

φ——煤体的内摩擦角，°；

k——支承压力峰值的应力集中系数；

h——煤层采深，m；

γ——上覆岩层平均容重，；

c——煤体的内聚力，mpa；

px——工作面支架对煤壁的支护反力形成的支护强度，mpa；；

ft——支架工作阻力，mn；

α——支架立柱倾角，°；

b——支架宽度，m。

第一初始压缩及塑性变形区a与第一煤体的强烈压缩及破坏区b的分界线的曲线方程式可以表示为：

x=0.0784h²+0.376h+x0

则由此可以计算破碎区宽度lp为：

lp=0.0784m²+0.376m+x0

1）第一初始压缩及塑性变形区a煤体体积的确定

由于第一煤体的强烈压缩及破坏区b为分界线下方x0至lp所围成的部分；设工作面长度为d米，煤层厚度、即采高为m米，则第一初始压缩及塑性变形区a的煤体体积va为：

2）第一初始压缩及塑性变形区a的煤体瓦斯涌出量

设工作面截深为e米，采煤机割煤速度为v米/分钟，根据现场统计，工作面割一刀煤的时间l/v一般不会小于15分钟，当15≤l/v≤30时，破碎煤体的瓦斯涌出量为总瓦斯含量的60%；当31≤l/v≤60时，破碎煤体的瓦斯涌出量为总瓦斯含量的85%；61≤l/v时，破碎煤体的瓦斯涌出量为总瓦斯含量的95%；

由此可以计算当工作面截深为e米、采煤机割煤速度为v米/分钟、原始相对瓦斯含量为q立方米/吨时，割煤后a区瓦斯绝对涌出量：

式中：

k——瓦斯释放均匀系数，根据瓦斯赋存的均匀程度可取1.0-1.5；

——瓦斯释放程度系数，当15≤l/v≤30时，取0.6；当31≤l/v≤60时，取0.85；61≤l/v时，取0.95。

如图2所示，步骤（4）中区域运输巷1的右侧帮和区域回风巷2的左侧帮煤体的瓦斯绝对涌出量q的具体确定过程为：

第二初始压缩及塑性变形区a1和第二煤体的强烈压缩及破坏区b1之间的分界线与第三初始压缩及塑性变形区a2和第三煤体的强烈压缩及破坏区b2之间的分界线左右对称布置，

其截距及该煤壁至分界线与煤层底板界面交点的水平距离x1可以用下式表示

px1——煤壁前方的超前压力

根据分界线的曲线方程式：

x=0.0784h²+0.376h+x1

则由此可以计算破碎区宽度lp1为：

lp1=0.0784m²+0.376m+x1

超前压力影响下的第二初始压缩及塑性变形区a1和第三初始压缩及塑性变形区a2的破碎情况与沿工作面方向是相似的，因为超前压力影响有一定的范围，取破碎的范围近似为20m；

则超前压力影响下的两帮的破碎区的体积

原始相对瓦斯含量为q立方米/吨时，在超前压力的影响下，两帮煤壁的瓦斯涌出量按50%来计算，则煤壁两帮的瓦斯涌出量q1+2为：。

本实施例并非对本发明的形状、材料、结构等作任何形式上的限制，凡是依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰，均属于本发明技术方案的保护范围。