基于热-流-固耦合的松软煤层瓦斯抽采数值模拟研究

基于热一流-固耦合的松软煤层瓦斯抽采数值模拟研究摘要：为了探索热-流-固耦合作用下煤储层渗透率的变化规律，本研究应用ONBOLMltiphysics探讨抽采时间为30d60d90d120d时山西潞安五阳煤矿8003工作面煤层瓦斯抽采效果研究结果表明：在不同抽采时间下渗透率变化较大；随着抽采时间增加，渗透率随之变小；渗透率较大的区域瓦斯渗流速度快，瓦斯抽采效率高；随着抽采时间增加，煤体孔隙率逐步减少；距离钻孔越远，煤体孔隙率越增加用煤层气是在煤化活动过程中产生的以甲烷和一氧界气候变暖的条件下，已变成一种清洁优质资源，有效五”计划推行，煤层气开发和利用被视为占重要地位的新能源产业，得到政府大力支持。同时，煤层气开发和利用在我国能源政策中占据较大比重，在减少煤层气排放时，对环保提供较大帮助。我国煤层气可采储量较国煤层气开发的主要因素之一11。煤层气开发的主要动和温度升高等因素影响，并对瓦斯抽采效率有较大影效应的综合作用过程。随着矿井进入深部开采，原岩应煤体透气性降低导致瓦斯抽采难度越来越大12]。其中，煤的渗透率不仅与本身孔隙、裂隙和煤层尺度相关，还由于煤层的瓦斯压力、气体含量、温度和地应力等多种因素而变化，结合这些因素可以得出煤层变化特性。研究瓦斯抽采过程中煤层渗透率变化规律及对抽采的影近年来，众多学者对瓦斯抽采的热-流-固耦合条件下煤层渗透率的变化规律进行探索，张民波等³1考虑抽采时煤层温度的变化，基于煤体渗透率的各向异性建立热-流-固三场耦合模型，研究温度和瓦斯压力对有效抽采区域的影响，发现有效抽采区域面积随着温度和瓦斯压力增加而减少，随着抽采时间增加而减少；李明等[4采用耐腐蚀气液两相渗透仪，进行不同埋深不同温度无烟煤渗流试验，研究有效应力、孔隙压力和温度对渗透性的影响规律，分析热刺激下滑脱效应的变化规律；范超军等⁵建立深部煤层瓦斯抽采流-固-热耦合模型，进行深部煤层瓦斯抽采数值模拟实验，研究了渗透率、瓦斯压力和温度等深部煤层参数对瓦斯抽采的影响；郝建峰等⁶J基于不同压力梯度下温度变化规律，并进行煤体吸附/解吸瓦斯热效应实验，建立含瓦斯煤层热-流-固耦合模型，分析吸附/解吸热效应在煤与瓦斯相互作用中是否起到较大作用；舒才等考虑解吸热能项优化了温度场控制方程，基于煤层瓦斯两能态吸附热理论构建煤层瓦斯流动的热-流-固多场耦合模型，分析煤层瓦斯流动规律，并通过对比实验结果验证模型准确性和合理性。为研究煤层瓦斯压力、渗透率和孔隙率对抽采效果的影响，以山西潞安五阳煤矿8003工作面煤样为对象，采用数值试验研究抽采时渗透率和孔隙率的变化对抽1煤岩与瓦斯热-流-固耦合模型①假定随煤层温度提升，处于吸附态的瓦斯被瞬间解吸，含瓦斯煤岩体被视为各向同性的线弹性均与多孔基体；②煤岩体在因地应力产生形变过程中符合广义胡克定律；③瓦斯在煤体中的渗流遵循Darcy准则；煤体中处于吸附态的瓦斯遵循Langmuir方程；④煤体为线弹性体且变形微小可恢复，遵守广义胡克定律。71.2建立物理模型基于太沙基有效应力原理，并考虑煤岩吸附瓦斯产生的吸附膨胀应力和温度产生的变形，结合应力方程、几何方程及本构方程得到应力场方程：表示变量μ的i方向分量，第二个下标j表示对μ求j方向偏导数，第三个下标j表示对μ;求j方向偏导数；为煤体骨架热膨胀系数，1/K;ε,为煤体骨架吸附瓦斯程为：式中q为瓦斯在煤体中的渗流速度，m³/s;k为煤根据质量守恒方程，在煤岩体中的连续性瓦斯流动方程为：煤中的瓦斯含量M可表示为：朗缪尔气体常数。联立得裂隙内瓦斯压力随时间变化的表达式：吸附量，m³/kg;b为气体的吸附常数，1/MP煤层中，煤体骨架、瓦斯、水共同存在于一个单元体积系统，温度变化造成煤体中内能变化，煤体因地应力造成的体积应变发生变形做功，在瓦斯吸附解吸过程中伴随热量释放和吸收。同时，煤体与外界通过热对流和传递产生热能交换。因煤体遵守能量守恒定随温度升高，煤的瓦斯吸附量逐渐减少。根据吸附试验得到不同温度时的吸附常数，采用二次函数和线性方程拟合得到吸附常数a、b与温度T的关系：式(1)、(5)、(6)、(7)、(8)构成含瓦斯煤层热-流-固耦合模型，应力场方程中的压力、温度项改变会引起煤岩体变形；流动方程中体积应变、瓦斯压力和温度项改变会使孔隙率和渗透率的方程结果发生变化，说明瓦斯抽采效率受到煤体变形影响。1.3耦合项煤层是孔隙和裂隙结合的双重结构，瓦斯主要储存在孔隙中，大多在裂隙通道中进行瓦斯渗流。根据有效应力、温度变化和热膨胀系数等参数，得到煤体的孔隙度方程为：单位体积煤体瓦斯解吸引起的膨胀应变可表示如下：因温度变化造成煤基质内发生热膨胀效应，所以，温度变化后得到含瓦斯煤在线弹性阶段时的动态渗透率方程可表示为：1.4工程背景试验背景为山西潞安五阳煤矿8003工作面，主采3#煤层，8003工作面走向长1595m,倾斜长315.8m,平均厚度5.37m,煤层最大原始瓦斯含量18.6m³/t,平均10.76m³/t,煤层初始瓦斯压力1.23MPa。1.5几何模型整个抽采过程是在三维空间进行，设计二维几何模型进行数值模拟将问题简化，便于计算。建立长在模型内设置如图1所示4个监测点在煤层钻孔平(10#0)、(20#0)、(30#0)、(40#0)通过监测每个点不同时刻的瓦斯压力、渗透率、孔隙率和温度得出各因煤层上边界受地应力为12.1MPaf左右两侧设置辊支撑开下边界为固定约束四边界为零通量边界模型左侧边界和右侧边界为煤层壁面上下边界为壁面用煤层初始瓦斯压力1.23MPa4分别研究抽采时间为30d、60d90d和120d的率、孔隙率和煤层瓦斯压力分析对抽采效果的影响实验煤样采用五阳煤矿8003工作面模拟所需的部分参数通过实验获得少数需要复杂实验获得的参数有效热对流系数η泊松比2采用OMEOLMltiphysics软件进行数值实验用抽采时间为250d选取30~120d煤层内部的瓦斯压力分布图如图2所示可以看出随着时间推移煤层钻孔周围的瓦斯压力向外逐渐减少#抽采30d抽采90d图2不同抽采时刻下煤层瓦斯压力分布图时间/d图3不同抽采时间下瓦斯压力变化图由图3可看出随着与钻孔中心距离增加瓦斯压力减少幅度变缓P监测点从左到右随抽采时间增加瓦斯压力慢慢变小用抽采120d后开各个监测点的瓦斯压力从左到右分别为0.55MPa0.MPa0.82MPa0.91MPa#距钻孔中心位置越远开瓦斯压力越小#因可以清晰看出随着抽采时间增加瓦斯压力在煤层分布由外到内逐渐减少并持续向外扩散用抽采后监测点A处的瓦斯压力降到0.75MPa瓦斯压力减少37%抽采90d后监测点A处瓦斯压力降低到0.57MPa瓦斯压力相较第30天时减少53%抽采天时瓦斯压力减少62%通过对比同一点处的瓦斯压力变化开可以发现煤体瓦斯压力随抽采时间增加而逐渐降低对比云图可以发现开随着距钻孔距离变远开煤层瓦斯压力有显著增加用走向长度/m图4不同抽采时间下渗透率变化图在打孔过程中影响煤层的初始应力平衡扩大了裂隙和孔隙范围钻孔周围的瓦斯压力变小开煤层瓦斯压力逐渐变小开孔隙率和渗透率随之增大开瓦斯抽采量增加并通过对比监测点发现瓦斯抽采对距离钻孔较远的瓦斯压力影响较小用9天时渗透率减少33%在抽采90d后，监测点A处渗透率为0.67×10¹6m²,相比抽采30d时渗透率减少56%;透率变化可以得出，煤体渗透率随着抽采时间增加而孔钻进过程中对应力平衡的破坏造成收缩进而导致煤低；受到的影响随着与钻孔距离增加而逐渐减少，瓦斯五阳煤矿为背景建立瓦斯抽采模型在一定边界条件下，得出各监测点瓦斯压力分布图以及不同抽采时间下孔隙率和渗透率随抽采时间变化的关系图用扩大裂隙和孔隙范围，钻孔周围的瓦斯压力变小，煤层瓦斯压力逐渐变小，孔隙率和渗透率随之增大，瓦斯抽采量增加，并通过对比监测点发现瓦斯抽采对距离钻孔较远处的瓦斯压力影响较小用透率随着抽采时间增加而减少，煤层瓦斯压力降低，煤基质受到地应力作用及钻孔钻进过程中对应力平衡的破坏造成收缩，进而导致煤体孔隙率变小，渗透率减少用通过对比孔隙率，可以发现随着抽采孔形成，在孔周围形成卸压区，这些区域的瓦斯压力降低使煤体产生压缩变形导致基质收缩，造成孔隙率降低；受到的影响随着与钻孔距离增加而逐渐减少，瓦斯压力变化幅度平缓，孔隙率变化较小用通过观察图5可知，在同一时刻监测点从左到右，孔隙率逐渐增加，距离抽采孔越远孔隙率越大抽采30d时，监测点A处孔隙率为0.043抽采60d后，监测点A处孔隙率减少到0.029,孔隙率相比抽采30d时下降34%抽采90d时，监测点A处孔隙率为0.023,孔隙率相比30d时减少50%当抽采120d时，A点处孔隙率为0.021#通过对比同一点处孔隙率变化可以发现，煤体孔隙率随着抽采时间增加逐渐降低；同一时刻随着监测点从左到