使用迭代方法对采空区流态区域进行划分.pdfVIP

【科技创新论坛 】 疆 VALLE-_L 使用迭代方法对采空区流态 区域进行划分 周 林锦 (衡水学院 数学与计算机科学系 河北 衡水 050030) 摘 要 ： 矿井瓦斯是煤矿生产过程 中主要的不安全 索，对矿井瓦斯运移规律预测是瓦斯防治不 可缺少 的重要技术环节。采空区流态 区域划分是采空区瓦斯 运移规律探索最基础的 步 为此将对多孔介质非线性渗流运动方程采用改进 的迭代方法对采空区流态 区域进行划分。 关键词 ： 瓦斯运移 ；流态区域 ；渗流方程；迭代方法 中图分类号：V2113 文献标识码 ：A 文章编号：1671 7597(2011)0320195～01 1采空区内的三种流态 3，2算法步骤 步骤1 对所给的方程 ，(x)=0恰当的选择迭代函数 ()，将方程变 形为 = (x)： 步骤2 提供初值 ，迭代求 出 = Xo)： 图1 气体在采空区的流动示意图 和遗煤 的压实程度 的差_骨很大，与进 步骤3 对预先给定的精度 ￡，若 l一Xol≥￡，则以 ·替换 ，从步 风巷 、回风巷 的距离相差也很大 ，因此气体在采空区各点的渗流流速也相 骤2继续计算；若 —Xol￡，则取X．为近似根，计算结束 差较大，采空区内各处渗流雷诺数的变化也较大。国内外大量的研 究表 4改进的迭代形式 明，采空区内一股 同时存在层流流区、紊流流区和过渡流流区，但是人们 将式5改为迭代形式： ”=一 +bllVII)。。Vp (6) 通过实验又发现，在靠近开采工作面 个很小的范尉内，由于漏风风速较 其中k为迭代步数。 大，其流场是紊流 、过渡流，而在采空区的其他区域则是类似于小雷诺数 记 ：～(口+b ) 贝： Vk+= vp (7) 的层流 [1]。 其中式 (7)运用变系数达西渗透迭代 求解 非线性渗 流方程的迭 代 2多孔介质非线性渗流运动方程 式。当k=l时即为变系数达西渗流。所 以可以说非线性渗流是从达两渗流 非均质冒落采空区是 个‘紊流、层流、过渡流三者同时存在的非线性 中不断迭代收敛得到的。 渗流场 ，其大空隙介质流体运动 的雷诺数Re很大 ，符合Bachmat提 山的多 为了加快迭代的收敛速度 ，引入收敛因予以，@ ，将式 (6)改为： 孔介质非线性渗流运动方程。 1z =一q(+6 ff ) 其 中 为初始收敛因子 ，叫∈[0，l】；鸥为振荡阻尼收敛 凼子。 P．_，： (1+ )q (1) g V 引入 的依据是初始迭代之前的速度为零，第一次迭代结果远远高 式中： 为压力梯度；e为渗透率 ；g为重力加速度 ；n为多孔介质的 于最终的收敛结果，此 时 应取小值。在q 引入后 ，迭代的结果仍然是 空隙度； 为介质颗粒的形状系数；d为多孔介质平均粒径 ；g 为比流量的 波动振荡的，为降低波动振荡的幅度 ，有必要 引入 ，其作用是削弱迭 模；q为比流量；v为运动粘性系数，对 于采