是不是COMSOL的文献才这么爆炸？！

来源公众号：COMSOL仿真交流

瓦斯抽采

煤层静态爆破预抽瓦斯

煤层静态爆破预抽瓦斯是高强度开采矿井保证工作面安全生产的重要保证，低浓度瓦斯的王家岭矿在强力开采状态下瓦斯会大量涌出，为提高矿煤层预抽的针对性和有效性，利用COMSOL软件模拟矿井12316工作面瓦斯抽采有效半径随抽采时间的演化规律。

 - 01 - 

静态爆破膨胀致裂力学分析

将致裂剂在钻孔壁膨胀致裂简化为壁厚圆简体模型，取钻孔边缘一单元体abcd进行研究，钻孔形变后单元体ad边位移至a'd',故得出环向和径向应变为：

εφ，表示环向应变；ερ 表示径向应变；ρ表示单元体上任一点距原点的距离；u表示钻孔边缘各点沿半径方向的位移；φ表示角度。 

单元体在各应力作用下的静力平衡方程为：

致裂孔内膨胀压仅施加于孔壁内部，致裂剂每一单元点合力为 0。 

圆筒体应力表达式为：

σρ表示压应力; σφ表示拉应力；ρ1表示膨胀压。

膨胀压变化到定值，内部结构开始破坏，围岩出现损伤。致裂孔内部，在钻孔r1半径处逐渐传递到r0处时，膨胀压由ρ0变化为ρ1，ρ1可表示为：

经过致裂压应力和拉应力作用，通过迭代计算得到当静态致裂产生的膨胀压到达钻孔内壁上的膨胀压近似值约为 40 MPa，传递后的膨胀压值根据实测数据作为现场的考察依据。煤层气被认为是理想气体，符合达西定律。气体泄漏的解吸是在等温线上处理的，并遵循理想气体方程。

 - 02 - 

建立模型

数值模型

1）多物理场模型假设。

① 煤层内部应力应变服从摩尔库仑准则；②吸附瓦斯遵循朗格缪尔吸附等温线方程；③瓦斯渗流运动符合连续性方程和达西定律。 

2）工作流程。

最初确定模型的建立，在抽采瓦斯的过程中，把钻孔周围的瓦斯流动视作二维径向流动，由此在建立数值模型的过程中选取 2D状态空间模型。

其次是物理场接口的选取。模拟场景为低渗透低瓦斯厚煤层在静态爆破下的瓦斯有效抽采半径，使用一般形式偏微分方程接口和固体力学接口。 

最后，添加研究状态为稳态。在研究过程中测定数据需要在瓦斯抽采效果稳定后确定，研究稳态的情况。 

3）连续性方程。

连续性方程是瓦斯流动的基本方程，质量守恒原理的流体力学表达式：

ρg表示瓦斯密度；Vg表示瓦斯的渗流速度；Q表示单位体积煤中的瓦斯含量；t 表示时间变量。 

4）瓦斯渗流场方程。

如果认为煤层中的瓦斯流动是线性的、稳定的水平流动，则符合达西定律：

ηg表示瓦斯动力粘度系数；k表示煤体渗透性；p 表示煤层瓦斯压力。 

5）煤岩体的变形方程。

假设煤层为弹性体，则煤体变形本构方程可表示为：

G 表示拉梅常数；e 表示体积变形；εij表示体积应变；δij 表示 Kronecker 符号。

模型参数

将模型四周假设为零流量边界，煤层的左右两侧设立棍支撑，底部设置固定约束，模型侧面设置水平约束，模型顶部和钻孔边界为自由边界，在模型顶部的自由边界上设置一个 10 MPa 的竖直向下的应力表达上覆岩层载荷，并在模型的钻孔壁面施加 40 MPa的载荷表达静态爆破所产生的膨胀力。模型的初始速度场和位移场为 0。

几何模型

在数值模拟中，将钻孔附近的瓦斯流动看成径向流动，并选择垂直钻孔方向的二维断面来研究在抽采过程中瓦斯压力和浓度的变化。模型如图 2 所示，尺寸为：长度 15 m、高度 8 m，其中煤层厚度为 6m，顶板厚度为 2 m，底板厚度为 1.3 m。在煤层中部布置了一个爆破孔和两个抽采孔，爆破孔和抽采孔之间的水平距离均为 1.5 m，爆破孔和考察孔的孔径均为 0.94 mm。根据上述参数进行网格划分，如图 3所示。

 - 03 - 

模拟结果分析

煤层塑性区变形

对于煤体这种弹塑性材料屈服就是破坏，即煤体发生塑性应变的区域就是煤体受到破坏的区域。由图 4 可知，爆破孔周围塑性区分布明显大于抽采孔的塑性区，爆破孔的塑性区呈现纺锤体状分布，塑性区向左右两侧横向扩展约 0.45 m，纵向向上延伸0.729 m，向下扩展 0.740 m。抽采孔的塑形区呈蝶翼状分布，塑性区半径约 0.12 m。这表明在爆破孔壁施加膨胀力后，使爆破孔周围煤体发生了显著的塑性变形，也大大促进了煤层裂隙发育，导致该区域范围内煤体的透气性大幅提升。

瓦斯有效抽采区域分布

随抽采时间增加，煤层整体瓦斯压力均呈下降趋势。但值得注意的是，随着抽采时间的增长，钻孔周围瓦斯压力要比其他位置瓦斯压力的下降速度更快，即任意一点的瓦斯压力的下降速度与该点距钻孔中心的直线距离成正比。

PS：后台私信“瓦斯”获取完整文献