三天前！国人团队又用COMSOL跑模拟了

来源公众号：COMSOL仿真交流

最新文献

# 派言派语 写在前面 #

在双碳能源转型过程中，煤炭仍然发挥着"压舱石"的作用。随着开采深度的不断增加，煤层赋存条件变得更加复杂，地应力增大、低透气性、高瓦斯压力等，煤矿瓦斯动力灾害形式更加严峻。

然而，对于煤层瓦斯赋存不均匀条件下，考虑钻孔周围塑性破坏的钻孔有效抽采半径和钻孔间距的研究较少。

鉴于此，2024年3月26日，河南理工大学安全科学与工程学院研究人员在ACS Omega上发表论文成果。

DOI: 10.1021/acsomega.4c00494

为了解决瓦斯赋存不均匀时钻孔布置的盲目性以及忽略钻孔周围的塑性破坏会低估瓦斯产量的问题，研究者建立了考虑塑性破坏的气固耦合模型。

他们利用有限元软件COMSOL Multiphysics探究了瓦斯非均匀赋存过程中有效半径的演化规律，确定了有效半径与钻孔间距的关系。

结果表明：考虑塑性破坏时，钻孔周围渗透率显著增大，且上升区域呈蝴蝶状，最大增幅达94倍。不仅瓦斯抽采流量峰值增大，而且其出现的时间也有所延迟。

最大瓦斯抽采流量增加40.7%，累计抽采量增加37%。得到了非均匀瓦斯赋存条件下有效半径与合理钻孔间距的关系。该研究可为煤层瓦斯非均匀分布条件下钻孔间距的选取提供参考。

# 内容梳理 结果讨论 #

一、模型概述

1、工作面瓦斯赋存情况概述

以某煤矿某工作面为例，该工作面煤层原始瓦斯压力在1.0~1.6MPa之间，瓦斯含量在10~20m3/t之间。工作面埋深530~660m，煤层平均厚度5.4m。

图1：某矿及某工作面瓦斯赋存状态

气体压力呈现不均匀分布，如图1所示。瓦斯分布主要分为4种状态，原始瓦斯压力分别为1.0，1.2，1.4，1.6MPa。

2、几何模型及定解条件

为研究不同赋存条件下的有效抽采半径，建立如图2所示的几何模型，长30m，宽25.4m，煤体厚度为5.4m。模型受上覆岩层重力影响，计算过程中施加16MPa的垂直主应力和10MPa的水平主应力。

图2：瓦斯单孔抽采数值模拟的几何模型

左侧设置为滚轮边界，下边界为固定支撑。模型外侧设置为不透水边界，表明在这些边界处没有气体流动。

二、结果讨论

1、考虑钻孔周围损坏对排水的影响

以煤层初始瓦斯压力为1MPa为例，模拟含瓦斯煤在塑性破坏后主要参数的变化。提取AC线上的渗透率。

考虑塑性破坏后，钻孔周围的渗透率显著增加，隆起区域呈蝶形，如图3所示。钻孔周围区域分为破坏区、塑性软化区和弹性区。

图3：抽采钻孔施工后钻孔周围的渗透率比值

2、非均匀瓦斯赋存有效半径的演化规律

在不同初始瓦斯压力下，随着抽采时间的增加，残余瓦斯压力不断降低，如图4所示。由于受煤层赋存等条件的影响，残余瓦斯压力降低的区域呈现椭圆形压降区。

在瓦斯压力的演化过程中，可以发现有效影响区的范围扩展比有效抽采区的范围扩展更快。

图4：不同抽气时间压力分布下瓦斯赋存不均匀

3、非均匀瓦斯赋存钻孔间距的确定

考虑钻孔与渗流速度之间的叠加效应，将排水钻孔周围区域依次划分为开放渗流区、过渡区和未影响区。渗流开放区对应瓦斯有效抽采区，渗流过渡区对应有效影响区，如图5所示。

图5：钻孔周围瓦斯渗流状态分布

根据防止张开渗流区域之间交叉重复影响和考虑过渡渗流区域之间相互叠加的原则，孔间距S应满足：2倍有效排水半径<s<< span="">有效影响半径。因此，不同初始瓦斯压力下钻孔的合理间距如图6所示，即图中的充填区域。

图6：瓦斯赋存不均匀情况下钻孔间距的布置范围

总之，研究者提出一种考虑塑性失效的气固耦合模型。研究了非均匀瓦斯赋存条件下有效半径的演化规律，得到了有效半径与钻孔间距的关系，可为非均匀瓦斯赋存条件下钻孔间距的选择提供参考。