带数值模拟报告.doc

带数值模拟报告———————————————————————————————— 作者：———————————————————————————————— 日期： 1 个人收集整理 勿做商业用途5 1121（3）采空区“三带”数值模拟采空区自然发火防治的首要技术前提是要搞清采空区自燃“三带”的分布范围，它是采空区注氮、堵漏风等技术措施的主要技术依据因此采空区“三带”的划分显得尤为重要.5 1 1121(3)工作面概况1121(3）工作面位于东一采区13—1煤层第二块段，煤层呈块状及粉末状，煤层产状变化较大，煤层倾向40°左右，倾角13～42°，煤厚30~6.0m，平均煤厚44m.工作面平均长130m，宽525～605m，高32~4.0m上风巷标高：-448.0m，下顺槽标高：—566.0m工作面顶、底板岩性如表5-1所示该工作面南侧有1111（3）工作面，与1111（3）工作面下顺槽留设10m煤柱，与1121（1）工作面相邻，高差均在80～100m左右。

其中1111(3)已回采完毕，1121（1）正在回采工作面采用区内后退走向长壁式，一次采全高综合机械化采煤法，平均推进速度为93m/月根据《精查地质报告》,本矿井恒温带深度为30m，温度168℃,地温梯度23℃/100m工作面实际温度在26~30℃工作面实际供风量为1800 m3/min该煤层具有自然发火性，自然发火期3～6个月.瓦斯相对涌出量为2.52m³/t，绝对瓦斯涌出量为884m³/min,矿压为21—24mpa表5—1 工作面顶、底板岩性岩石种类备注老顶细砂岩厚度为0～7.2/5.6m，灰白色、细粒为主直接顶泥岩伪顶炭质泥岩厚度为0～04/0.2m，灰黑色，主要为炭质泥岩直接底泥岩老底细砂岩、中细砂岩厚度为10.0~180/16.5m,浅灰～灰白色,坚硬5 2 采空区“三带”划分指标工作面正常生产时,采空区自燃“三带"处于一个动态的稳定状态采空区自燃“三带”主要指散热带、氧化带和窒息带散热带内由于冒落不充分,漏风流较大，采空区遗煤气化产生的热量不能积聚，一般不会发生自然发火氧化带内漏风风速适当，具有热量积聚的条件，O2体积分数又能满足氧化需求，因此最容易发生自然发火窒息带内由于漏风难以到达，O2体积分数往往难以满足氧化需求，一般也不会发生自然发火。

三带”是客观存在的，但如何划分具有一定的困难由于探测手段和方法的局限，想要定量地准确划分是难以做到的目前，一些研究者提出确定划分“三带”的指标有漏风风速（v)、采空区氧浓度和温升速率3种：1、采空区漏风风速(v).从理论上说，漏风风速相对较好.因为它可以体现氧浓度分布、氧化生热与散热的平衡关系.B.MMaebcka(1946）等研究表明，采空区及煤柱的漏风强度在时容易自然发火；有的研究者认为不会导致自燃的极限风速低于；封闭采空区密闭墙漏风压差在300Pa、漏风强度在时容易自燃发火.由于这些参数都是在一定条件下取得的,因此有一定局限性，但对研究自燃问题是有一定的参考价值2、采空区氧浓度(c）分布.采用氧浓度指标不能划分散热带和自燃带.因为在自燃带中氧浓度也有可能达20％以上因为氧化初期，耗氧量是很小的.对于划分自燃带和窒息带的指标有的研究者认为是5％～6％;有的认为是8%,即氧浓度c〈8％为窒息带．C〉8％为自燃带.实际上，氧化速度随氧浓度降低而减小，究竟氧化速度小到何值才算之是“窒息”是没有确切标淮；3、采空区遗煤温升速度（dt〉1OC/d为自燃带）也是一个划分“三带”的指标,但目前应用不是很广泛。

按照采空区内漏风风速指标划分可知,散热带:采空区内漏风风速大于0.24m/min的区域;氧化带：漏风风速在024～01m/min之间的区域；窒息带:漏风风速小于01m/min的区域考虑到1121(3）工作面的推进速度快，采空区浮煤还没有足够的氧化时间，氧气浓度随工作面的推进,变化并不明显因此，本章主要采用专业的流体力学模拟软件FLUENT对1121(3)工作面采空区进行流场模拟，并根据采空区内漏风风速划分指标来确定1121(3）工作面采空区“三带”的位置及范围最后改变工作面的供风量，得到工作面不同供风条件下的采空区“三带”位置及范围5.3 物理模型的建立根据潘北矿1121(3)工作面现场实际提供的相关资料,为了使模拟的结果更能够反映工作面及采空区的实际流场分布情况，把上、下顺槽，工作面和采空区三者结合起来一起作为采场整体来进行数值模拟考虑到现场实际的上、下顺槽及采空区的走向很长，而我们主要是关注工作面采空区“三带”的位置范围，因此，在对“三带”模拟结果影响不大的前提下，模型进行了简化上、下顺槽的长度都取20m,采空区走向取150m，其他模型尺寸完全和现场实际一致根据潘北矿1121（3）工作面采空区实际情况，建立复杂的采空区三维物理模型，如图5—1所示,模型尺寸如表5—2所示。

表5-2 物理模型尺寸模型名称长/m宽/m高/m备注工作面130.05.636长：指倾向长采空区150.013003.6长：指走向长上顺槽20.032下顺槽2043.0图5-1 采空区物理模型示意图上顺槽采空区150m130m20m5.6m下顺槽5.4数学模型的建立一、连续性方程的建立根据质量守衡原理直接可以得到连续性方程为： 其中：ui，uj，uk为单元体的平均流速，它和孔隙中的平均流速关系为：N——孔隙率.以下数学模型中，除特别指明外,所指速度都是指单元体的平均流速二、动量方程的建立在惯性（非加速）坐标系中i方向上的动量守恒方程为：式中:p——静压，Pa;-—应力张量,Pa;且应力张量由下式给出: ——包含了其它的模型相关源项，如多孔介质和自定义源项；采空区内是属于多孔介质流体区域，可由下式得出:其中D和C是规定的矩阵在多孔介质单元中，动量损失对于压力梯度有贡献,压降和流体速度(或速度方阵）成比例对于简单的均匀多孔介质:其中a是渗透系数，为内部阻力因子，简单的指定D和C分别为对角阵和，其它项为零通过多孔介质的层流流动中,压降和速度成比例，常数可以考虑为零5.5 FLUENT软件简介FLUENT是美国FLUENT公司推出的用于模拟具有复杂外形的流体流动以及热传导的计算机软件。

它提供了很好的网格灵活性，可以使用结构化网格和非结构化网格,例如二维三角形或四边形网格、三维四面体∕六面体∕金字塔形网格来解决具有复杂外形的流动,也可以用混合型非结构网格它还允许用户根据解的具体情况对网格进行修改(细化或粗化）FLUENT软件采用有限容积法，可选择一阶迎风、二阶迎风、中心差分、幂指数以及QUICK等格式速度和压力耦合采用SIMPLE或SIMPLEC算法，代数方程的求解可以采用多重网格及最小残差法FLUENT是用C语言写的，具有很大的灵活性与能力.除此之外,为了执行高效，易于交互控制，以及灵活的适应各种机器与操作系统,Fluent使用client∕server结构，因此它允许同时在用户桌面工作站和强有力的服务器上分离地运行程序.在FLUENT中，解的计算与显示可以通过交互界面、菜单界面来完成用户界面是用Scheme语言及LISP dialect写成的高级用户可以通过写菜单宏及菜单函数自定义及优化界面.FLUENT软件包括以下几个部分：(1）FLUENT 解法器；（2）prePDF,用于模拟PDF 燃烧过程；(3）GAMBIT，网格生成几何图形模拟以及网格生成的预处理程序;(4)TGrid，从已有边界网格中生成体网格的附加前处理程序；（5）Filters（Translators)，转换其它程序生成的网格，用于FLUENT计算.可以与ANSYS、I—DEAS、NASTRAN、PATRAN等程序进行数据交换。

在FLUENT中可以用GAMBIT产生所需的几何结构以及网格，也可以在已知边界网格（由GAMBIT或者第三方CAD/CAE软件产生的）中用TGrid产生三角网格，四面体网格或者混合网格也可以用其他软件产生FLUENT所需要的网格，比如ANSYS（Swanson Analysis Systems.Inc)、I—DEAS (SDRC)，或者MSC/ARIES,MSC/PATRAN以及MSC/NASTRAN大多数CAD/CAE软件都可以产生上述格式的网格FLUENT软件解算器对于可压缩与不可压缩流动,稳态和非稳态流动，无粘流、层流及湍流,牛顿流体及非牛顿流体，对流换热（包括自然对流和混合对流），导热与对流换热耦合，辐射换热，惯性坐标系和非惯性坐标系下的流动，多运动坐标系下的流动，化学组分混合与反应,多孔介质流动,两相流，复杂表面形状下的自由面流动等问题都可以进行较好得模拟因此，FLUENT软件在航天、机械、电子、汽车工业、材料处理、建筑设计和火灾研究以及暖通空调等领域都得到了广泛的应用5.6 采空区“三带”模拟条件模拟采用隐式求解的求解方法模拟需要的其它已知条件如下：1、工作面的推进速度取3。

1m/d；工作面供风量为2、采空区倾向上的风阻变化系数在倾斜方向上，由于采空区的上、下巷道附近顶板受实体煤的支撑，而采空区中部的上覆岩层无支撑，所以顶板岩石冒落的下沉量也不同一般是采空区中部覆岩的下沉量大于其上、下边界覆岩的下沉量因此,在倾斜方向上,一般采空区中部较上、下巷附近的孔隙率小，风阻大令k(y）表示风阻变化系数，则可得：H1——采空区中部覆岩的下沉量，m；H2——上、下边界覆岩的下沉量，m；——采空区倾向长,m;根据矿井实际条件及相关资料，H1/H2取130.3、采空区漏风风阻的计算由采空区漏风风阻的计算公式为:式中：a，b——经验系数，取决于顶板冒落岩石性质,其值如表5-3所示；c——冒落岩石的压实系数，它与工作面推进速度(m/d）有关，按式(54)计算: 式中：x—-采空区内距工作面距离，m；l——采空区滤流分支的长度，m；s——采空区滤流分支的截面积，m2；-—工作面推进速度,m/d；考虑采空区倾向的风阻变化,则可得采空区漏风风阻的计算公式为：由采空区走向长为150m，工作面长度为130m，工作面的推进速度取：31m/da，b系数，根据实际顶板岩性：a=0。

04，b=40,则可以得到： 表5—3 不同岩性冒落岩石的a、b值冒落岩石种类ab松软粘土岩，页岩101中硬粘土页岩硬粘土页岩、砂岩砂岩、石灰岩4、模拟基本假设1）、采空区内遗煤、冒落岩石与空气等混合物视为各向同性的均匀一致的多孔介质；2）、上、下顺槽在同一水平面，不考虑位能的影响；3)、粘性阻力系数和内部阻力系数在z方向不发生变化5.7 采空区“三带”模拟。