特厚煤层综放工作面分层开采瓦斯渗透规律研究.docxVIP

特厚煤层综放工作面分层开采瓦斯渗透规律研究 在几年前，在煤层下方附近的煤层中发现了肥沃的砖瓦。同时，由于煤层开采的影响，这两座煤层之间的岩石存在裂缝。瓦斯流体原有的均衡状态被改变,瓦斯从下层煤层通过中间岩层贯穿裂隙向上层采空区和工作面流动,使采煤工作面瓦斯浓度增大,对采煤工作面的安全推进构成威胁。王济凯、张英华等建立了采空区渗流数学模型,揭示采空区漏风流速、压力分布规律;李晓平等通过数值分析研究了岩层水侵识别不稳定渗流理论;赵阳升、赵洪宝等通过对薄煤层开采时采空区内瓦斯运移进行模拟研究,获得采空区瓦斯分布特性规律;许满贵等利用Fluent数值模拟研究了采空区瓦斯运移规律,优化了上隅角瓦斯治理的参数;范宏伟采用Fluent软件对U型和U+L型通风方式下综放采空区瓦斯流动进行了数值模拟,研究瓦斯流动分布规律。在综放面有距离较近的邻近煤层时,其采空区瓦斯运移是否受邻近煤层瓦斯渗透的影响及影响程度,在工程实践中是一个不可回避的问题。 本研究采用数值方法,结合现场探测数据,分析中间岩层裂隙受采动影响的变化规律,着重研究邻近煤层高瓦斯渗透对上部煤层采空区瓦斯运移的影响及邻近煤层瓦斯涌出规律。总结邻近煤层瓦斯渗透对综放面瓦斯运移的影响规律,为工程实践中考虑邻近煤层瓦斯渗透因素提供了理论参考。 1 煤岩渗透率预测 在采煤工作面正常推进过程中,煤层底板岩层具有采前应力升高,采后应力降低及恢复3个阶段,并且这3个阶段随采煤工作面推进而重复出现。 煤矿井下采矿工作会使底板所受应力重新分布,使底板岩层的渗透性增高或者降低,表现为卸压区渗透性增高,应力集中带透气性降低。因而,采煤工作会使底板岩层渗透率发生变化,具体表现为下部采空区的瓦斯涌出量忽小忽大。文献通过实验室的模拟实验,得到卸载过程中渗透率κ与有效应力σ之间的关系符合幂函数方程: κ=κ1σ-C。 (1) 其中:κ—煤(岩)体渗透率,m2;κ1—煤(岩)体初始渗透率,试验测定值为0.733×10-22,m2;σ—有效应力,MPa;C—岩层孔隙压缩系数,无纲量。 本文所研究煤岩层特征与文献所采用的岩样性质基本一致,因此,可以借鉴式(1)进行岩层渗透率的预测。首先运用UDEC(universal distinct element code)软件模拟工作面开采过程中底板应力分布,然后将卸载区的应力值代入式(1)计算出卸载区底板岩层渗透率。 2 煤层应力分布规律 根据现场实例建立UDEC模型,模型高×宽=300 m×300 m,下部邻近煤层采厚5 m,开采方式为采13 m留5 m煤柱的小条带式开采;上下煤层间距为26 m;上部煤层采厚3 m,采用长壁开采。 首先,对下部煤层进行条带开采,然后,对上层煤层采煤工作面从切眼推进20 m和200 m两种情况进行模拟,得出开采后上下煤层间岩层应力分布规律(图1)。工程实践中,综放面推进20 m时,顶板未发生垮落现象,此阶段为模拟初采阶段;采煤面推进至200 m时,综放面支架后方10 m内未垮落,10 m外发生垮落现象,在模拟中垮落带以弱化材料充填,将综放面后方10 m范围以外由底板向上10 m范围划分为一个单元,并降低该单元岩性,使其岩性低于相邻区域岩层岩性。 由图1可见,上下煤层间的岩层应力原始值为4.0 MPa,下部煤层条带开采及上部煤层长壁开采后,采空区内,下部煤层顶板大约6.0 m范围应力降低为2.0 MPa,初采阶段上部煤层底板5.0 m范围应力降低为2.0 MPa,3.0 m范围应力甚至出现正值(拉应力);正常开采阶段,上部煤层底板5.0 m范围(局部地段的底板15 m范围)应力降低为2.0 MPa。上述卸载区破坏范围与现场实测破坏范围基本一致,说明上述模拟过程是正确的。 将上述应力值分别代入式(1),求得有效应力值为4.0 MPa岩层的渗透率为0.47×10-23m2,有效应力值为2.0 MPa的岩层的渗透率为1.86×10-23m2。 3 gambit建模与fluent值模拟 3.1 数值计算模型 根据相邻煤层岩层结构(图2),对煤层工作面推进距离20 m,200 m两种情况分别建立Gambit数值计算模型(图3(a),(b)),依据上下煤层间岩层渗透率变化规律,将其分为上、中、下三层。在综放面正常推进阶段,根据垮落矸石压实程度的不同,将上部煤层采空区分为“三带”:自然堆积区、载荷影响区、压实稳定区,相关参数见表1和表2。 3.2 工作面流场模型 将Gambit软件建立的数值计算模型(图3)导入Fluent模拟软件中并进行求解。模拟采用k-e湍流模型。该模型是一种可处理气体流动的湍流状态的流动模型,该湍流模型用于求解处于模型湍流核心区的流动,取C1=1.44,C2=1.92。对于模型近壁区的流动,采用标准壁面函数进行描述,对标准k-e