于 涛刘 晔毛肖杰

（1.煤科集团沈阳研究院有限公司，辽宁省沈阳市，113122；2.北京科技大学土木与环境工程学院，北京市海淀区，100083；3.昆明理工大学，云南省昆明市，650093）

煤田露头火区多孔介质渗流运动的数值模拟研究

于 涛1刘 晔2毛肖杰3

（1.煤科集团沈阳研究院有限公司，辽宁省沈阳市，113122；2.北京科技大学土木与环境工程学院，北京市海淀区，100083；3.昆明理工大学，云南省昆明市，650093）

煤田露头煤层自燃使得煤田露头发生塌陷，为了合理地治理煤田露头火区，研究了影响煤田露头火区的重要参数，应用Fluent数值模拟软件针对不同煤层埋藏深度、孔隙率、渗透率进行了宏观数值模拟。得到了煤田露头多孔介质速度场、温度场变化规律，并根据孔隙率与渗透率对火区温度分布的变化规律提出治理方案。

煤田露头火区 塌陷区 孔隙率 渗透率 数值模拟

煤田煤层露头在自然通风的影响下发生自燃，随着露头煤层燃烧深度的加深，煤田火区地表塌陷。煤田露头火区煤层燃烧后的塌陷区上覆岩层在垂直方向上产生冒落带及其充填区、裂隙带和弯曲下沉带，“三带”高度可按照岩层控制理论来确定。煤田露头火区多孔介质的孔隙率及渗透率影响到火区的温度场及渗流速度场分布规律，进而影响煤田露头煤层自燃的结果。因此，研究煤田露头火区多孔介质的孔隙率及渗透率的关系，有利于科学分析煤田露头火区燃烧过程和防治原理，对做好煤田火区防治工作，保护煤炭资源和生态环境具有重大的意义。

1　煤田露头火区多孔介质气体渗流形式的分析

煤田露头由露头煤层自燃产生的塌陷形成塌陷“三带”，“三带”煤岩、土层的破碎构成了多处进风口和出风口。为了研究方便，本文根据相关理论进行合理简化，将煤田露头火区简化为一主要的进风口和出风口。即新鲜自然风流由进风口通过冒落带及其充填区构成的多孔介质通道，再经过塌陷裂隙带和弯曲下沉带的裂隙及孔隙，最终由出风口排到地表，其中部分气体由冒落带及其充填区的孔隙和裂隙渗入燃烧区，在火源区域与煤层发生氧化燃烧，随后燃烧的烟气沿岩层裂隙、孔隙逸出地表，其过程如图1所示。塌陷 “三带”由于孔隙率和渗透率的差异，得到的是不连续的渗流形式，本文为研究孔隙率及渗透率的关系，将其非均匀区忽略，认为风流在塌陷 “三带”渗流运动过程中是连续的。

图1　风流在煤田露头火区渗流运动路线示意图

风流在煤田露头火区内可能存在两种流体，一种是煤自燃生成气体，另一种是外界进入的新鲜风流。且存在两种流体同时流动的情况，煤自燃生产的气体和新鲜风流一起在塌陷区 “三带”中裂隙及孔隙形成的回风通道做渗流运动，从而形成了煤田露头火区通风系统。

2　煤田露头火区多孔介质孔隙率与渗透率关系分析

煤田露头火区多孔介质的孔隙率和渗透率是物性参数中重要的参数，也是影响风流在火区内渗流速度的重要因素之一。确定孔隙率与渗透率的关系对研究煤田露头火区内渗流速度场、温度场分布规律及燃烧区火势大小变化规律提供一定的理论支持。本文主要介绍裂隙带与弯曲下沉带、冒落带及其充填区研究孔隙率与渗透率关系的结果。

（1）在煤田露头火区塌陷 “三带”中，裂隙带与弯曲下沉带主要是由松散岩土及岩块构成的多孔介质组成的，根据试验得到的孔隙率及渗透率关系曲线，孔隙率与渗透率关系函数基本符合幂函数：

式中：K——渗透率，m2；

φ——孔隙率。

（2）在煤田露头火区塌陷 “三带”中，冒落带及其充填区的自燃产物为完全燃烧的余煤及破碎岩块组成的多孔介质，新鲜空气和烟气主要以层流和分子扩散方式运动，根据相关渗透率研究成果，冒落带及其充填区孔隙率与渗透率关系可表示为：

式中：K0——初始渗透率，取1×10-9m2。

3　煤田露头火区数值模拟模型

煤田露头煤层自燃过程是自主完成的，不同孔隙率与渗透率则对火区气体渗流运动及火灾形势有不同影响，进而影响露头煤层的燃烧。根据矿区火区地质资料，进行合理简化后建立物理模型，并根据相关资料对煤田露头火区作如下假设：

（1）由于露头煤层燃烧使得地表塌陷，导致露头火区的进出口较多，本文将其合理简化为一个主要进风口和回风口；

（2）假设气体在塌陷区 “三带”多孔性介质中为稳定流动；

（3）根据相关试验得假设气体在塌陷区多孔性介质中的运动为层流且线性运动；

（4）不考虑液体影响，气体为理想气体；

（5）假设煤田露头火区塌陷区 “三带”中的裂隙带及弯曲下沉带为一整体且孔隙率和渗透率是不变的。

根据桌子山煤田骆驼山矿区煤样及相关样本测得，其煤层塌陷区的孔隙率在0.1～0.4之间。取孔隙率为0.2和0.4两种情况对煤田露头火区进行数值模拟。

煤层塌陷区的孔隙率在0.2和0.4时，由式（1）计算得裂隙带与弯曲下沉带对应的渗透率为5.104×10-11m2和4.0832×10-10m2；由式（2）计算得冒落带及其充填区对应的渗透率为5.1867 ×10-11m2和7.3767×10-10m2。

煤田露头煤层自燃时，假设煤氧化反应在区域块内进行，则煤岩露头自燃时3个区域的组分i的渗流方程如下。

露头煤层燃烧区内孔隙中渗流方程：

煤田露头塌陷区冒落带及其充填区区内孔隙中渗流方程：

煤田露头塌陷区裂隙带及弯曲下沉带区内孔隙中渗流方程：

式中：w——渗流气体的组分i的质量分数；

P——压力，Pa；

g——重力加速度，N/kg；

μ——动力粘性系数；

ρ——密度，kg/m3；

γi——区域内孔隙中组分i与裂隙中组分i的交换量；

Si——气体组分i的源项。

4　煤田露头火区数值模拟结果及分析

根据上述相关理论分别对两种孔隙率进行数值模拟，并对两种孔隙率在煤层距露头地面燃烧区深度（埋藏深度）L为60 m和70 m位置处分布进行模拟，得到对应的煤田露头火区渗流速度和温度场分布规律如图2、3所示。

图2　煤田露头火区渗流速度分布云图

通过图2火区渗流速度分布云图可以看出，煤田露头火区进、出风口的风速较高，且随着火区距煤层露头距离的增加，气体渗流速度逐渐减小；在靠近燃烧区时渗流速度几乎为0，此时气体主要以分子扩散运动方式进入火区。此外，在不同的煤层埋藏深度、不同孔隙率下，其渗流速度、分布范围不同。渗流速度随孔隙率增大而增加，渗流速度分布范围也随之增大，但随着煤层深度的加深，渗流速度会随之减小。对比图2中（a）、（b）或（c）、（d）可得在煤层相同埋藏深度时，不同孔隙率，气体在冒落带及其充填区和裂隙带、弯曲下沉带的渗流速度场分布不同。孔隙率与渗透率是影响气体在煤田露头火区渗流速度的主要因素，说明孔隙率越大在单位时间渗入燃烧区的氧气量增大。由此可以验证孔隙率与渗透率关系在不同塌陷地带的关系式也不相同，即在冒落带及其充填区孔隙率与渗透率的关系满足式（2）；而在裂隙带及弯曲下沉带孔隙率与渗透率的关系式满足式（1）。这为研究煤田露头火区煤层燃烧速度、气体渗流速度的关系及煤田露头火区火灾治理提供了一些理论支持。

由图3可以明显的看出在温度场分布范围由大到小排列顺序为图3（d）＞（c）＞（b）＞（a）。从而看出在相同火源温度下不同孔隙率、不同煤层埋藏深度，温度场的分布影响范围不同，煤层火势也不同。由图3（a）、（b）或（c）、（d）可知，孔隙率由0.2到0.4时高温范围明显变大。在入口风速不变的情况下孔隙率增大，温度场的分布范围也增大，说明了单位时间渗入燃烧区的氧气量增大，煤层的火势也增大。孔隙率增加，一方面增加了反应面积，直接加快化学反应；另一方面增强了高温气体与固体骨架间的对流换热面积，间接加快了化学反应。因此，孔隙率增加对煤自燃（燃烧）起到了促进作用。由此在煤田露头火区火灾治理时可以考虑在减小孔隙率方面做工作，如往塌陷区及滑体裂隙注浆封堵煤田露头火区原生裂隙和次生裂隙，从而减小塌陷的孔隙率和渗透率，进而减小或隔绝煤体与氧气的接触，实现灭火目的。

图3　煤田露头火区温度分布云图

5　结论

（1）得到了煤田露头火区时气体渗流运动的主要影响因素关系。分析塌陷区孔隙率与渗透率关系式，确定孔隙率与渗透率关系在不同塌陷地带的关系式也不相同，即在裂隙带及弯曲下沉带和冒落带及其充填区的孔隙率及渗透率关系分别满足式（1）和（2）。

（2）模拟结果表明，在不同的煤层埋藏深度、不同孔隙率下，其渗流速度、分布范围不同。渗流速度随孔隙率增大而增加，渗流速度分布范围也随之增大，但随着煤层埋藏深度的加深，渗流速度会随之减小。

（3）随着孔隙率的增大单位时间渗入燃烧区的氧气量增大，煤层燃烧火势也随之增大，由此在煤田露头火区火灾治理时可以考虑在减小孔隙率方面做工作，为煤田露头火区火灾治理提供技术路线。

［1］ 曾强，王德明，蔡忠勇.煤田火区裂隙场及其透气率分布特征［J］.煤炭学报，2010（10）

［2］ 孙若山，王海燕，周心权.煤田自燃表征系统的灰关联模型及应用［J］.中国煤炭，2011（4）

［3］ 梁运涛，张腾飞等.采空区孔隙率非均质模型及其流场分布模拟［J］.煤炭学报，2009（9）

［4］ 王海燕，周心权，张红军等.煤田露头自燃的渗流-热动力耦合模型及应用 ［J］.北京科技大学学报，2010（2）

［5］ 王铁利.渗透率与孔隙度、束缚水、饱和度的关系［J］.煤炭技术，2010（1）

［6］ 姚秀云，张凤莲，赵鸿儒.岩石物性综合测定——砂、泥岩孔隙度与深度及渗透率关系的定量研究［J］.石油地球物理勘探，1989（5）

［7］ 何学秋，王恩元，林海燕.孔隙气体对煤体变形及蚀损作用机理［J］.中国矿业大学学报，1996（1）

［8］ Yuan L，Smith A C.Numerical study on effects of coal properties on spontaneous heating in longwall gob areas［J］.Fuel，2008（87）

［9］ 段品佳，王芝银.煤岩孔隙率与渗透率变化规律试验研究［J］.地下空间与工程学报，2013（6）

［10］ 但汉成，李亮等.变形多孔介质流固耦合非线性渗流模型及解析［J］.长江科学院院报，2009（6）

［11］ 李唐山.煤田露头自燃气体渗流和传热传质规律研究［D］.北京：中国矿业大学（北京），2006

［12］ 陈国新.煤田露头自燃双重介质渗流规律及发展趋势预测［D］.北京：中国矿业大学（北京），2008

（责任编辑 张艳华）

Numerical modeling study on the seepage of porous media of coalfield outcrop fire zone

Yu Tao1，Liu Ye2，Mao Xiaojie3
（1.Shenyang Research Institute of China Coal Technology&Engineering Group，Shenyang，Liaoning 113122，China；2.School of Civil&Environment Engineering of University of Science&Technology Beijing，Haidian，Beijing 100083，China；3.Kunming University of Science&Technology，Kunming，Yunnan 650093，China）

Coal spontaneous combustion of coalfield outcrop made coalfield outcrop subsided. In order to reasonably manage the coalfield outcrop fire zone，the important parameters affecting coalfield outcrop fire area were studied and the macroscopic numerical modeling was made for different buried depth of coal seam，porosity and permeability by Fluent simulation software.The variation law of velocity field and temperature field of porous medium of coalfield outcrop was obtained，and the treatment scheme was put forward according to the effect of the porosity and permeability on the temperature distribution in the fire area.

coalfield outcrop fire zone，subsidence area，porosity，permeability，numerical modeling

TD752.2

A

于涛（1982-），女，汉族，辽宁营口人，工程师，研究方向：矿井通风与安全、煤与瓦斯鉴定，瓦斯涌出量预测等。