基于lattice Boltzmann方法对裂隙煤体中瓦斯运移规律的模拟研究

基于lattice Boltzmann方法建立一个新的模拟裂隙煤体内瓦斯渗流的动力学模型,并利用该模型模拟二维裂隙煤体内瓦斯流动。模拟结果表明,采动压力差对瓦斯流动速度、孔隙瓦斯压力及瓦斯压力梯度都有很大的影响。在采动压力条件下,瓦斯在裂隙煤体中的流动压力具有波动性,波动幅度的大小与流场两端的压力差有关。在瓦斯压力达到峰值前的瞬时,孔隙瓦斯压力梯度很大。随着采动压力差的增大,瓦斯流动由层流过渡到紊流,同时孔隙压力出现发散,压力变化出现明显的非线性特征。基于lattice Boltzmann方法的模拟结果与用其他方法得到的瓦斯渗流规律比较吻合,表明lattice Boltzmann方法可有效模拟瓦斯在裂隙煤体中的运移规律,这为进一步探讨煤与瓦斯两相耦合机制、煤与瓦斯突出机制及瓦斯抽放方案的设计提供新的思路。     

基于lattice Boltzmann方法对裂隙煤体中瓦斯运移规律的模拟研究

基于lattice Boltzmann方法建立一个新的模拟裂隙煤体内瓦斯渗流的动力学模型,并利用该模型模拟二维裂隙煤体内瓦斯流动.模拟结果表明,采动压力差对瓦斯流动速度、孔隙瓦斯压力及瓦斯压力梯度都有很大的影响.在采动压力条件下,瓦斯在裂隙煤体中的流动压力具有波动性,波动幅度的大小与流场两端的压力差有关.在瓦斯压力达到峰值前的瞬时,孔隙瓦斯压力梯度很大.随着采动压力差的增大,瓦斯流动由层流过渡到紊流,同时孔隙压力出现发散,压力变化出现明显的非线性特征.基于lattice Boltzmann方法的模拟结果与用其他方法得到的瓦斯渗流规律比较吻合,表明lattice Boltzmann方法可有效模拟瓦斯在裂隙煤体中的运移规律,这为进一步探讨煤与瓦斯两相耦合机制、煤与瓦斯突出机制及瓦斯抽放方案的设计提供新的思路.     

瓦斯抽采过程中的煤层透气性动态演化规律与数值模拟

 为了研究煤层瓦斯抽采过程中的煤体渗透性变化规律,基于Kozeny-Carman方程,利用表面物理化学与含瓦斯煤的有效应力理论,建立考虑有效应力变化、瓦斯解吸和煤基质收缩效应的煤层渗透率动态变化模型,并结合数值模拟分析煤层瓦斯抽采过程中煤体透气性动态演化规律。研究结果表明：(1) 所建立的煤层渗透率动态演化模型能较好地描述煤层瓦斯抽采过程中的煤体透气性动态演化规律。(2) 煤体渗透率与煤体孔隙压力之间呈现出“V”字型变化趋势,低瓦斯压力阶段煤基质收缩效应占主导地位,煤层渗透率随瓦斯压力降低而增大;高瓦斯压力阶段有效应力作用占主导地位,煤层渗透率随瓦斯压力降低而减小。(3) 从煤层内部逐渐接近抽采钻孔过程中,煤层瓦斯压力较高时,煤体渗透率先减小后增加;煤层瓦斯压力较低时,煤体渗透率不断增大。研究结果可以为我国煤矿瓦斯治理和煤层瓦斯抽采提供理论支撑,具有指导性意义。     

基于扩散-渗流机理瓦斯抽采三维模拟研究

根据煤体孔裂隙多重介质特性,分析煤体内瓦斯在基质孔隙和裂隙中不同的运移机理,建立瓦斯流动方程;考虑孔隙压力、吸附膨胀效应对渗透率的综合影响,建立渗透率动态方程;以弹性力学为基础,建立煤体变形控制方程;共同构建能够描述瓦斯扩散-渗流运移机理的耦合模型。以真实煤层赋存参数为例,根据理论模型对COMSOL Multiphysics进行二次开发,建立三维瓦斯抽采模型,模拟研究不同抽采时间下的瓦斯压力、抽采有效半径变化规律;利用压降法进行现场实测考察,现场实测结果与模拟结果基本吻合,验证了扩散-渗流耦合模型及三维数值模拟的可靠性,从而为该矿实际瓦斯抽采工作提供安全可靠的理论技术指导。     

深部低透气性首采层煤与瓦斯气固耦合模型

立足于消除煤层渗透及扩散特性对于煤与瓦斯气固耦合模型的干扰,在分析首采煤层所处应力状态特点的基础上,建立更符合煤体的孔隙裂隙二重介质特性的修正的P-M渗透率模型,提出考虑解吸–扩散效应及Klinkenberg效应的煤与瓦斯气固耦合模型,详细阐述多物理场之间的耦合作用关系。应用该模型模拟分析深部首采层顺层钻孔预抽消突过程中煤层瓦斯压力及渗透率的演化规律。模拟结果表明,Klinkenberg效应对低渗透煤层瓦斯运移的促进作用显著,并随着瓦斯压力减小促进效果增大;煤体绝对渗透率的动态变化是骨架压缩效应及基质收缩效应的竞争结果,瓦斯压力开始减小时,骨架压缩效应首先起主导作用,渗透率减小,瓦斯压力持续降低时,基质收缩效应逐渐取代其成为主导作用,渗透率增大。     

综放开采围岩活动影响下瓦斯运移规律及其控制

 博士学位论文摘要　将综放开采矿山压力新理论、覆岩结构及其活动特征与瓦斯运移聚集规律有机结合, 通过现场矿压观测、适时瓦斯监测、相似材料模拟、电液伺服试验、理论分析及有限元解算等方法, 对含瓦斯厚煤层实施综放工艺时顶煤顶板活动规律及其内瓦斯运移聚集形态进行了系统深入的定性定量分析。从分析含瓦斯煤岩体微观结构特征入手, 探讨了矿压作用下煤岩与瓦斯的宏观运动形态, 从而得出综放开采富含瓦斯煤层的安全控制理论与实践技术。论文视含瓦斯煤岩体为不连续的多相介质材料, 并深入到材料的微观机制抽象出其微结构模型, 分析了其结构形态在煤岩体变形和瓦斯运移过程中的变化形式; 定量分析了煤岩体内孔隙裂隙分布及破碎后的块度变化的分形维数值, 认为采动后支承压力等因素使煤岩体破坏, 从而影响其分维值, 煤岩的极限粒度越小,瓦斯放散解析速度越快, 瓦斯运移潜能越大, 表征其吸附特性的V an derW aals 力下降越多; 综放开采使煤岩力学介质变化显著, 支承压力作用下煤体经历初期压缩、强烈压缩、滑移与离散、流动放出4 个过程, 顶煤介质状态可由宏观连续介质转化为碎裂介质和块裂介质。基于魏家地矿和阳泉五矿的综放面及巷道矿压观测、瓦斯监测和采空区束管监测, 认为本煤层开采时支承压力使煤体孔隙率、渗透系数和瓦斯压力发生很大变化, 且支承压力变化状态与瓦斯运移变化在时间和空间上有一致性; 首次将综放面前方煤体划分为无压瓦斯自由放散区、卸压瓦斯涌出活跃区、降压瓦斯涌出变化区、升压瓦斯涌出变化区及稳压瓦斯正常涌出区; 周期来压时, 瓦斯绝对涌出量和相对涌出量分别比来压前增加了44% 和54%; 统计国内11 个矿井生产面, 覆岩关键层初次来压时, 采场瓦斯绝对涌出量是来压前的1. 21～ 3. 44 倍, 相对涌出量是初垮前的1. 21～ 6. 75 倍; 关键层失稳经历不同, 使综放采空区后方垮矸碎胀系数具有明显的分区性, 据计算碎胀系数和瓦斯浓度、氧浓度变化可将其划分为自然堆积区、载荷影响区和压实稳定区, 自然堆积区和载荷影响区是分析瓦斯流态的关键区, 其内瓦斯的紊流和过渡流态可用非线性渗流方程表示。关键层初次破断失稳的动态冲击力, 一方面克服了风流压力和漏风阻力将富集于支架上方断裂煤壁及架后采空区瓦斯压向采场, 另一方面在碰撞作用点附近使采空区垮矸空隙率降低形成模拟墙, 从而使墙体前方瓦斯迅速挤向采场; 顶板来压前工作面支承压力达到最大, 使直接顶、顶煤体剪切破坏, 煤体屈服, 卸压范围扩大, 支承压力与瓦斯压力梯度联合作用使煤体瓦斯解析并向采场涌出; 关键层作为板的破断呈“O 2X”型特征, 在其交点处首先形成导气通道, 顶板来压期间关键层及其覆岩因变形特性不同而不协调垮落, 便将离层裂隙聚集的瓦斯通过“O 2X”破断裂隙挤入采场, 因此得出结论: 综放面瓦斯大量快速涌出是矿山压力的一种显现。综放开采使垮落带和规则移动带高度增加, 为瓦斯运移聚集提供了较大活动范围。关键层与其上覆或下伏岩层间不协调变形将形成覆岩离层裂隙和破断裂隙; 在煤层采厚2. 6～ 3. 4 倍高度以下破断裂隙较发育, 其上以离层裂隙为主, 随综放开采两类裂隙的时空发展有明显的三阶段特征, 即切眼和回采面附近覆岩采动裂隙发育, 采空区中部裂隙则被重新压实; 覆岩垮落及离层高度受关键层及其结构的影响而呈动态变化, 随顶板初次来压和周期来压, 离层高度呈跳跃性变化, 从而为裂隙带内瓦斯运移聚集的剧变提供了时空条件。综放开采前期是瓦斯运移及控制的关键时期, 覆岩采动裂隙带是经破断与离层裂隙贯通后在空间形成关键层下似椭圆抛物面内外边界所包围的椭抛带(EPZ) 分布, 椭抛带层面的切割为椭圆形裂隙发育区; 关键层破断后, 裂隙带宽在初采边界处相当于初次来压步距, 在综放面上方则变化在1～ 2 倍周期来压步距之间; 采动裂隙带的发生、发展基本受制于覆岩关键层层位及其所形成砌体梁结构的变形、破断和失稳形态; 当主关键层切割椭抛带时, 采动裂隙呈椭球台状, 层面展布的椭圆形裂隙区仍将存在。表征覆岩离层特征的位移曲线主要取决于关键层断裂块度的大小, 由此可计算出关键层初次失稳前后离层裂隙的当量面积及其空隙率和渗透系数。综放开采时瓦斯涌出特点决定了其在覆岩采动裂隙带内具有升浮和扩散两种运移方式。不均衡性瓦斯涌出带, 与周围环境气体存在密度差而升浮, 在浮力作用下沿破断裂隙上升过程中不断渗入周围气体, 使涌出源瓦斯与环境气体的密度差渐减至零, 瓦斯则会漂浮在离层裂隙发育区, 瓦斯升浮高度与本煤层及邻近层瓦斯含量及涌出强度成正比; 混入矿井空气中的瓦斯(CH4) 在其浓度梯度作用下会引起气体分子的普通扩散和压强扩散, 瓦斯扩散流方向与重力压强梯度反向, 即瓦斯具有向上扩散的趋势。从理论上解释了裂隙带是瓦斯运移及聚集带, 为覆岩裂隙带内钻孔抽放、巷道排放瓦斯技术提供了科学依据。首次提出煤样全应力应变过程中渗透系数是体积应变的双值函数, 体积缩小时为2 次多项式, 体积膨胀时为5 次多项式; 渗透系数在弹塑性段急增, 峰值后仍增大, 但梯度渐缓, 最大值发生在软化段或塑性流动段, 且与最小值相差上百倍; 主应力差增大时, 渗透系数变化范围增大, 反之则小。渗透系数是影响煤层瓦斯运移的最重要指标, 而支承压力则是影响渗透系数的主导因素; 支承压力作用下综放面前方不同部位煤体渗透系数变化范围相当大, 支架上方顶煤煤体及煤壁前方5m 内渗透系数最高, 塑性变形区内, 煤层渗透系数急剧降低, 到弹性变形区则接近原始值, 两极值相差可达数10 倍甚至数百倍。认为不论原始渗透系数怎样低的煤层, 采动影响下煤层卸压后, 其内瓦斯渗流速度大增, 瓦斯涌出量也随之剧增, 为瓦斯抽排提供便利条件, 由此提出“煤层与瓦斯共采”的新概念。有限元计算表明: 均衡推进的综放面, 采用短距离循环推进则可降低煤体中因渗流场结构变化而引起的瓦斯压力较大的波动, 一定程度上可减弱综放面前方煤体中瓦斯挤压和抛出煤体(动力异常) 的危险程度。提出以采场矿压监测为主的连续危险源非接触式法预测采场瓦斯大量涌出或涌出异常, 并成功预测了打通一矿工作面突出; 魏家地矿的应用实践表明, 三巷型布置较适宜富含瓦斯倾斜厚煤层的综放开采, 高抽巷应开掘在采动裂隙带内; 预采顶分层或开采解放层即可预释放大量瓦斯又可减缓综放开采的矿压显现程度;代替采空区井的采动区井(孔) 底处于裂隙带内能够充分抽排瓦斯, 淮北局的应用充分证明了该论点; 充分监测综放面顶板来压, 可有效防止综放采空区瓦斯爆炸; 充分利用矿压显现特点且有利于综放开采防治瓦斯、煤层自燃和煤尘等的煤体注水技术在魏家地矿取得了成功。     

考虑煤体孔–裂隙介尺度特征的水相渗流演化模型研究

研究煤体孔裂隙结构介尺度特征下的水相渗流特性对探究煤体孔裂隙结构的跨尺度渗流具有重要意义,为研究煤体孔隙、裂隙双重结构介尺度特征下的水相渗流规律,利用表征体单元的概念界定出水在煤体孔裂隙结构中渗流的介尺度下界,基于煤体孔隙结构和裂隙结构水相渗流模型引入裂隙和孔隙结构贡献率修正系数,构建煤体孔裂隙结构跨尺度的水相渗流模型,通过数值分析和注水实验对模型分析与验证并进一步确定模型的适用范围。研究结果表明,煤体的裂隙贡献率修正系数随着孔裂隙率的增加而增大,孔隙结构贡献率修正系数随着孔裂隙率的增加而减小,且二者均呈现出幂函数的变化规律;模型渗流速度与实验渗流速度均随着孔裂隙率的增加而增大,高孔裂隙率时模型计算结果与实验结果一致性较差,低孔裂隙率时模型的计算结果与实验结果一致性较高;实验过程中应力的施加对高孔裂隙率煤体中水的渗流产生较大的影响,低孔裂隙率煤体中水的渗流受应力影响不大;孔隙率小于3%时的数据变异系数均小于15%,平均变异系数为9.70%15%,表明模型的预测值与实际的实验数据能够很好地吻合,验证了模型的正确性,所建立的预测模型适用条件为介尺度和低孔裂隙率。     

气体压力加卸载过程中无烟煤变形及渗透特性的试验研究

利用自主研发的含瓦斯煤热流固耦合三轴伺服渗流装置,以无烟煤型煤试件为研究对象,进行不同轴压、围压条件下气体压力加卸载过程中渗流试验研究,模拟不同煤层深度,以探讨煤变形及瓦斯运移演化规律.研究结果表明:(1)在加载过程中,煤应变量减小,吸附瓦斯产生较大的膨胀变形,呈现线性关系,在卸载过程中,煤应变呈增大趋势,煤逐渐被压缩.随轴压、围压增大,下降单位气体压力引起的煤应变升高量降低,应变响应程度减小.(2)在加载过程中,随气体压力升高,渗透率先减小后增大趋势,煤渗透率呈类似“V”型变化趋势,气体压力在1.2 MPa左右存在明显的拐点,体现煤孔隙扩张的程度和吸附瓦斯层增厚程度影响,依赖于吸附作用或有效应力占主导地位.在卸载过程中,随着气体压力降低,煤渗透率呈先减小后增大趋势,渗透率增大且变化速度加快,主要依赖有效应力作用或基质收缩的主导地位差异.(3)随有效应力的增大,煤渗透率呈先减小后增大的趋势.煤渗透率随有效应力增大呈对数函数或指数函数关系.(4)气体压力具有典型二阶段特征,同时渗透率与体积应变具有密切关系,体现出有效应力、吸附膨胀与煤基质收缩同时对裂隙等内部结构的影响.     

采气–采煤阶段煤岩渗透率演化机制研究

为模拟煤矿现场先采气后采煤的作业过程,利用含瓦斯煤热–流–固耦合三轴伺服渗流装置,开展孔隙压力减小的煤岩渗流试验和全应力–应变–渗流试验。通过推导温度升高时煤岩裂隙宽度变化的表达式,进一步构建考虑温度–应力的煤岩渗透率模型,探讨温度与应力作用下煤岩瓦斯渗流演化机制。新建的煤岩渗透率模型包括有效应力、吸附/解吸、热膨胀及热裂四部分影响因素,并使用损伤变量表征裂纹扩容过程中产生的基质膨胀效应(热裂),结果表明:(1)当外应力恒定时,不同温度下渗透率随孔隙压力减小先略有减小而后迅速增大;当孔隙压力恒定时,渗透率随温度增大整体呈先减小后增大的趋势。(2)在煤岩全应力-应变–渗流试验过程中,渗透率随轴向应力的增大呈先减小后增大的趋势;煤岩弹性模量及峰值强度与温度之间呈负相关关系。(3)新建渗透率模型的计算值和实测值基本一致,可较好表征渗透率随孔隙压力及有效应力的演化规律。(4)基于内膨胀应力的定义,探讨温度与应力作用过程中内膨胀变形对渗透率的贡献。在温度恒定时,渗透率随内膨胀因子的增大而减小;煤岩裂隙宽度与损伤相关,在温度突变系数增大过程中渗透率随之减小。     

幂率型裂隙分布煤层渗流场与变形应力场耦合模型及数值模拟

煤体的微观结构对其产气量有着显著的影响,然而,现阶段的煤层气开采模型都未将此因素考虑在内。为了定量分析煤层微观结构对宏观渗流过程的影响,基于幂律分形渗透模型,构建一种能够同时考虑煤层微观结构与多物理场因素相互作用的裂隙–孔隙渗流模型。模型将多物理场效应(包括气体吸附–解吸效应、气体渗流诱导储层形变效应、煤层基质形变、裂隙–基质相互作用等)定义为储层有效应力的函数,进而作用于储层孔隙率与微观结构。分析煤体的主要结构参数对宏观渗透率的影响,包括：(1)裂隙长度幂律指数α;(2)裂隙长度最值比r;(3)最大裂隙长度l,并分析煤层形变及渗透率的时空演化趋势。模拟结果表明：煤体结构参数对煤体渗透率有着显著的影响。当孔隙率等其他煤体参数保持不变时,储层宏观渗透率与裂隙长度最值比r、最大裂隙长度l成正比,与裂隙长度幂律指数α成反比;此外,与裂隙长度幂律指数、裂隙长度最值比相比,煤层最大裂隙长度对煤层渗透率有着更为显著的影响。     

Analysis of coupled gas flow and deformation process with desorption and Klinkenberg effects in coal seams

Coupled gas flow and solid deformation in porous media has received considerable attention because of its importance in pneumatic test analysis, contaminant transport, and gas outbursts during coal mining. Gas flow in porous media is quite different from liquid flow due to the large gas compressibility and pressure-dependent effective permeability. The dependence of gas pressure and gas desorption on gas permeability has a significant effect on gas flow, but has been ignored in most previous studies. Moreover, solid deformation has a direct impact on the porosity, which also leads to desorption or sorption of methane in the coal seam. In this study, a coupled mathematical model for solid deformation and gas flow is proposed and is implemented using a finite element method. The numerical code is used to solve the gas flow equation with Klinkenberg effect, and is validated by comparison with available analytical solutions. Then, it is used to simulate the coupled process during gas migration in a deformable coal seam. The numerical results indicate that the desorption and Klinkenberg effects and mechanical process effect make a significant contribution to gas flow in the coal seam. Without considering the desorption and Klinkenberg effects and the coupling action of mechanical process, the gas pressure in the coal seam would be underestimated.     

变饱和裂隙多孔介质的双渗透率数值模型与参数分析

针对裂隙多孔介质中的变饱和流动与运移,分析了双渗透率数值模型的参数构成与计算特征。双渗透率模型由分别代表裂隙等效连续介质和孔隙介质的两组并列体积单元集合构成,用平均裂隙宽度和平均裂隙间距表达裂隙等效连续介质的参数,用裂隙单元与孔隙单元之间的有效作用面积、有效作用距离和有效渗透系数描述裂隙-孔隙相对流动和运移。裂隙-孔隙相互作用参数取决于概念模型,并与流体动态相关。计算分析表明,裂隙单元与孔隙单元之间的毛细力梯度与重力的共同作用决定裂隙流动、孔隙流动和裂隙-孔隙相对流动的相对大小,流体和溶质的动态受到裂隙宽度变化的显著影响。     

功率超声影响的煤中甲烷气促解规律实验研究

 功率超声激励促解煤层甲烷气是一种不受储层地质条件和气源特性限制,具有普遍应用价值的增采技术。通过CT观测实验,对超声波作用下煤样不同尺度裂隙发展规律进行深入分析,从微观上揭示功率超声促解机制。在对比促解实验测定结果的基础上,研究功率超声作用煤样解吸量变化规律。研究结果表明：CT观测实验很好地证明了超声的机械震碎作用;在功率超声激励作用下,煤样裂隙条数显著增多,贯通裂隙增多且单条裂纹最大宽度显著增大;超声波作用后煤中甲烷气的解吸量有显著提高,建立超声声强参数影响的煤中甲烷气解吸量随解吸时间变化的修正公式。由修正公式拟合结果可知：煤中甲烷气饱和吸附后的解吸量 和 均随孔隙压力的增加而增大,相同孔隙压力下的甲烷气饱和吸附解吸量 大于无超声作用时的解吸量 ,其原因在于在功率超声的声场强度影响下,煤质点的动能和位能增加,煤表面甚至基质内部吸附的甲烷更容易脱附变成游离态,从而达到煤层甲烷气促解作用。     

煤层瓦斯运移机制的关键参数表征

 煤层的瓦斯扩散系数、浓度流动系数、初始运移强度系数和衰减系数,是煤体孔隙结构和煤质特性的力学表征。为研究煤层瓦斯运移机制,量化煤层瓦斯运移能力,依据Fick扩散理论和质量守恒方程,建立了描述煤体瓦斯浓度与扩散速率的计算模型,采用变量分离法进行数学求解,并通过数据迭代方法获取煤体内瓦斯的扩散系数和表面浓度流动系数;通过室内4种煤样瓦斯运移实验数据比较发现,煤体瓦斯质量增量与运移时间成负指数关系,且随煤阶升高而增大;瓦斯运移速度和衰减系数取决于煤的吸附能力和煤质组分,且煤种之间差异明显;渗透率与扩散系数成线性关系,与流动系数成二次函数关系,且随煤阶升高整体呈增加趋势。     

Numerical investigation of coal and gas outbursts in underground collieries

Coal and gas outbursts are a complex catastrophic unstable phenomenon that involve the ejection of large volumes of coal, and are often accompanied by gas, such as methane, carbon dioxide or a mixture of the two. Coal and gas outbursts are prevalent in deep and gassy mines where face advance rates are rapid, and where gas drainage is either poor or absent. The occurrence of progressively larger coal and gas outbursts, and the potential for the catastrophic collapse of coal pillars, is of increasing importance as mining is extended deeper in seams rich in methane and other hydrocarbons. A unique coupled gas flow and solid deformation numerical model, viz., RFPA^2D-GasFlow, has been developed and is applied to simulate the evolutionary process of such catastrophic coal failures in underground collieries. The finite element model, which incorporates the physics of gas flow in the coal seam, the physics of coal deformation and instantaneous failure, and the cross-couplings between them, is proposed. The model also incorporates small-scale variability in deformation modulus and strength of the coal and surrounding rock. The variability in modulus and strength is distributed via a fine-scale resolution model according to the Weibull distribution, where the distribution parameter determines the level of heterogeneity. This numerical model is applied to simulate the whole process of coal and gas outbursts, including stress concentration, coal fracturing, gas pressure-driven expansion, and outburst. The instantaneous outburst process and associated stress fields, gas pressure gradients and displacement vectors are presented step by step. The numerical simulations indicate that the instantaneous outburst is a complex phenomenon involving interactions between gas pressure, stress and the physico-mechanical properties of the coal, and it can occur under a variety of conditions. Successful numerical simulation of the whole coal and gas outburst process provides the basis for identifying the outburst mechanisms, parameterizing the causative processes, and to defining potential precursors of failure.     

地球物理场中煤层气渗流控制方程及其数值解

 为探索地球物理场中原地煤层气运移能力对煤层气储集和富集能力的影响,以地应力场、地温场中煤层气连续性方程、气体状态方程、吸附方程、渗流方程为基础,建立了应力、温度影响下的煤层气渗流控制方程。方程体现了地应力和地温对煤层气压力、含量、渗透率和孔隙率的影响,其中,应力和温度通过影响煤层气压力影响吸附量,通过影响煤层气压力和孔隙率影响游离量;温度还通过影响吸附常数b影响吸附量;不同的应力、温度组合条件下,渗透率的变化机制不同。通过Kaiser声发射原岩应力测试实验、不同温度下煤的甲烷等温吸附实验、不同温度及有效应力下煤体中甲烷渗流实验以及煤的孔隙率、工业分析等实验,研究应力、温度影响下的煤层气渗流特征。不同温度下煤的甲烷等温吸附实验表明,吸附常数a随温度变化不明显,b随温度升高而下降;不同温度、不同有效应力条件下煤的甲烷渗流实验表明,小有效应力条件下,煤体中甲烷渗透率随温度升高而升高;大有效应力条件下,渗透率随温度升高而下降。以实验数据和原始地质资料为基础,采用有限差分法,进行了地球物理场中原地煤层气渗流运移能力的一维、二维数值模拟。计算表明：研究区现今原地煤层气渗流运移导致的煤层气散失甚微,低渗煤层具有良好的储集和富集能力,但不利于后期开采,卸除地应力和升高温度是提高煤层气抽采率的有效途径。     

煤层气越流固-气耦合模型及可视化模拟研究

应用煤岩体变形与煤层气（瓦斯）越流相互作用的新观点，研究了双层系统煤层气越流运移的机理。在地球物理场作用下，双层系统煤层气越流实质上是可压缩性气体在各向异性且非均质的孔隙与裂隙双重可变形介质中的渗透与扩散的混合非稳定流动过程。依据这一新认识，建立了双层系统煤层气越流与煤岩弹性变形的固-气耦合数学模型，并应用该耦合模型于邻近层煤层气涌出的可视化数值模拟实例中。经实测结果验证，该固．气耦合模型是符合实际的；同时，应用三维数值场计算机模拟技术，直观地再现了煤层气越流场与煤岩体变形场耦合效应的宏观变化特征，为地下煤层气越流场中邻近层煤层气涌出的预测和控制提供了可行的新方法。     

页岩气藏流固耦合渗流模型及有限元求解

 页岩气渗流模型是页岩气藏动态分析和数值模拟的基础。将裂缝性页岩气藏视为基质孔隙–裂缝双重介质,同时考虑岩石骨架变形对气体渗流场的影响,建立页岩气藏流固耦合渗流模型。模型假设基质孔隙内作克努森流动,裂缝中作达西渗流,综合考虑页岩气壁面滑脱流动与孔内扩散作用、吸附与脱附、应力敏感性等渗流机制。采用有限元法离散控制方程及全隐式耦合求解方法,编制计算机程序。考虑真实页岩参数取值,利用该模型进行算例分析。结果表明,页岩气藏压力下降速率小于常规裂缝性气藏压力下降速率;裂缝渗透率是影响裂缝渗流压力衰减的主要因素,需考虑页岩裂缝导流能力与基质产气速率的匹配关系;原始地层压力越小,裂缝渗流压力衰减越慢。所建模型可为页岩气藏模拟器开发及动态分析提供理论基础。     

煤/瓦斯突出过程中煤介质局部化破坏的损伤机理

提出了一种煤／ 瓦斯突出过程中煤介质失效破坏局部化的损伤机理的数值分析模型。模型的理论基础是将损伤力学与煤层中瓦斯气迁移的理论通过局部化处理相结合。模型中考虑了气体流动,瓦斯吸附,煤介质的变形及煤层中损伤的发展与传播间的相互影响。该模型被用于有限元分析。通过数值模拟阐述了以局部化损伤模型为基础的煤／ 瓦斯突出机理。