煤岩体瓦斯水渗流耦合过程数值模型及其在矿山工程中的应用ppt课件.pptVIP

模型建立 初始条件： 域内具有一个大气压，瓦斯初始浓度3mol/m3；推进新的工作面区域时，旧区域瓦斯浓度模拟结果C0再附加一初始补充浓度Cbc （取3mol/m3）作为旧推进工作面的初始浓度 。 模型建立 时间步长： 按照非均等的积数步长增大，初始值7s,终止指7e5s，设定100个中间时间值。 计算参数： 动粘系数=1.8e-5pa·s，流体密度=1200kg/m3，扩散系数D=2e-5m2 ，瞬态时间比例系数=0.55 模型建立 1 2 3 4 5 1 3.24e-8 2 2.88e-8 3.24e-8 3 2.52e-8 2.88e-8 3.24e-8 4 2.16e-8 2.52e-8 2.88e-8 3.24e-8 5 1.8e-8 2.16e-8 2.52e-8 2.88e-8 3.24e-8 已推区域 透气率 （m2） 新推区域 随推进进行透气率变化表 模拟结果分析 推进区域1时的计算结果 下图分别是当时间为：7s，7e3s，7e4s，1.1e5s，2.5e5s，3.5e5s，7e5s瞬态时间区域1采空区瓦斯浓度分布图 模拟结果分析 Time=7s时瓦斯浓度 模拟结果分析 Time=7e3s时瓦斯浓度 模拟结果分析 Time=7e4s时瓦斯浓度 模拟结果分析 Time=1.1e5s时瓦斯浓度 模拟结果分析 Time=2.5e5s时瓦斯浓度 模拟结果分析 Time=3.5e5s时瓦斯浓度 模拟结果分析 Time=7e5s时瓦斯浓度和流线 A1 A1’ 模拟结果分析 第一步推进时A1-A1’切面浓度变化曲线 模拟结果分析 推进区域2时的计算结果 下图分别是当时间为：7s，7e3s，1.1e5s，7e5s瞬态时间区域2采空区瓦斯浓度分布图 模拟结果分析 Time=7s时瓦斯浓度 模拟结果分析 Time=7e3s时瓦斯浓度 模拟结果分析 Time=1.1e5s时瓦斯浓度 模拟结果分析 第二步推进时Time=7e5s时瓦斯浓度和流线 A2 A2’ 模拟结果分析 第二步推进时A2-A2’切面浓度变化曲线 模拟结果分析 推进区域3,4时的计算结果 模拟结果分析 第三步推进时Time=7e5s时瓦斯浓度和流线 A3 A3’ 模拟结果分析 第三步推进时A3-A3’切面浓度变化曲线 模拟结果分析 第四步推进时Time=7e5s时瓦斯浓度和流线 A4 A4’ 模拟结果分析 第四步推进时A4-A4’切面浓度变化曲线 讨论 （1）在区域1右侧不透气边界附近，随着工作面向左推进，风流作用逐渐减弱，该处瓦斯浓度不断集聚增大，每个推进步瓦斯浓度分别为7mol，12mol，17 mol，24 mol，表明该处距离通风口越远，瓦斯集聚程度越大。 讨论 （2）对于每个推进步，当时间达到7e5s（8d左右）时，风流基本上可以把本区域瓦斯浓度降低到安全范围内（接近0mol），表明风流在本区域的流速明显大于前一个区域。而前一个区域的下左半部分（占整个区域的1/4）的瓦斯浓度也能得到有效降低，而对其他区域瓦斯浓度驱散作用影响有限。 讨论 （3）图9可见，瓦斯浓度沿工作面推进方向从近到远瓦斯浓度成台阶状渐次增大，这表明通风量只能在一定范围内降低瓦斯浓度，同时也看出距离通风口越远，瓦斯补给的时间和补给量越大。 讨论 （4）虽然模型上下边界都是瓦斯补给边界，但由于模型左下边界是进风口，瓦斯随风流对流作用得到有效降低；而模型左上边界是出风口，瓦斯容易在回风隅角大量聚集，所以该处瓦斯浓度较高。 三、三维渗流耦合模型及瓦斯抽放 根据实际的三维煤层瓦斯抽放过程，建立如图1所示的长宽高为10m×10m×10m理想化的三维计算模型，模型左下部边界假设为巷道，三个瓦斯抽放孔——K1、K 2、K 3之间的距离1m，与水平面呈45°倾角展布，K 1、K 3与K 2呈16°夹角分布在孔2两侧，三个瓦斯抽放孔长度7m，瓦斯抽放孔按照“以缝代孔”[16]原则简化为定压力边界（25kPa）。 巷道所在位置 45o 16o L1 L2 K1 K3 K2 外部载荷 (1) 边界条件：四周和上下面均为零通量不透气边界；四周和下部都约束法线方向的位移，上部自由，作用有P分别为10MPa、20MPa、2MPa、0.1MPa的外部载荷，同时模型具有自重载荷。 (2) 初始条件：内部有1MPa的初始瓦斯压力，三个抽放孔的压力为0.25e5Pa。 (3) 时间步长：按照非均等的积数步长增大，初始值1s，终止指1e7s（10d左右），设定100个中间时间步长。 (4) 计算方案：模拟不同外部载荷条件下（2MPa、10MPa、20MPa、0.1MPa），瓦斯抽放效果和渗透性变化规律。 (5) 计算参数：相关参数列表1所列。 生活家饮