考虑爆破因素的煤与瓦斯流固耦合作用规律研究.docxVIP

考虑爆破因素的煤与瓦斯流固耦合作用规律研究 煤层瓦斯流动机理 在一定浓度的范围内，一定浓度的砖瓦在特定条件下燃烧和爆炸。如果浓度过高，则会导致人员窒息和死亡，或加剧煤炭和瓦片的突出事故。掘进面人员和设备较为集中,该处的瓦斯防治工作尤为重要,煤矿瓦斯的研究及预防备受各方关注。 研究煤层瓦斯流动机理是掌握瓦斯流动规律的前提,以往多是单纯地研究掘进面应力、能量或者瓦斯分布和演化特征,很少对应力与瓦斯的耦合效应进行分析。本文分别建立了机掘和炮掘2种采掘工艺下,煤岩应力场和瓦斯流动场的耦合作用模型,利用comsol模拟了掘进面煤岩应力场与瓦斯流动场的耦合作用规律,以期为掘进工作面瓦斯动力灾害防治提供理论依据,使理论分析和研究更符合工程实际。 1 煤岩介质瓦斯 为了建立本文的流固耦合数学模型,需要遵循有关物理定律,并做出如下假设:煤岩体系小变形;煤岩为多孔介质;煤岩介质为单相的瓦斯所饱和,瓦斯气体储存于煤层中分别呈游离态和吸附态,煤层瓦斯含量遵循朗格缪尔方程;游离煤层气渗流运动及煤层变形运动的惯性力、煤层气的体积力可以忽略不计;流体在煤岩孔隙中的流动服从达西定律;渗流可以按照等温过程处理;煤层瓦斯按照理想状态气体处理。 1.1 机掘工艺的应力场方程 由瓦斯及煤岩骨架所组成的多相介质受到外力作用以后,煤岩固体骨架将产生变形,瓦斯同时产生流动或原来的流动状况发生改变。煤岩固体骨架的变形是在外部载荷和瓦斯压力的共同作用下发生的。 多孔介质固体骨架的有效应力表达式: 式中p为孔隙压力,MPa; h为孔隙度; 本构方程: 式中G为剪切模量,Pa; g为泊松比; a为biot压缩系数。 将式(1)~式(3)代入式(4),得到机掘工艺的应力场方程: 动载条件下的平衡方程为在静载条件下的平衡方程右端增加一项密度与速度偏导数的乘积外载荷项,表达式为 式中vi为煤岩运动速度(m/s); d为岩石密度(kg/m3)。 1.2 气体密度和压力关系 假设渗流过程中温度保持不变,即等温情况下,渗流控制方程由质量守恒方程和达西定律组成。 煤岩多孔介质由单相瓦斯气体所饱和,瓦斯以自由态和吸附态赋存于孔隙中,当计算单位体积煤岩中瓦斯含量时,其是自由态的瓦斯含量与吸附态的瓦斯含量的和。煤中吸附瓦斯的含量可通过Langmuir方程计算: 式中m为瓦斯气体的含量(kg/m3); h为煤岩多孔介质的孔隙度; dg为瓦斯气体的密度(kg/m3); dga为一定外界压力作用下瓦斯气体的密度(kg/m3),dga=UPa,其中Pa为一个大气压; VL为Langmuir体积参数(m3/kg); 式中qg为Darcy流速度(m/s); kg为渗透系数(m2); _g为动力粘度系数(Pa·s)。 理想气体状态方程反映了气体密度和压力之间的关系: 式中dg为气体密度(kg/m3); P为孔隙压力(Pa); U为气体的压缩系数(kg·m3·Pa-1),U=Mg/RT,其中Mg为气体摩尔质量(kg/kmol),R为瓦斯气体常数(J·kmol-1·K-1),T为气体温度(K)。 在等温条件下,气体在多孔介质中的流动遵循气体质量守恒方程: 式中m为瓦斯气体含量(kg/m3); 将达西定律、瓦斯气体含量方程、理想气体状态方程代入质量守恒方程: 1.3 煤层渗流场边界 (1)流场定解条件: (1)初始条件为,其中p0为煤层中始瓦斯压力(Pa); (2)渗流场边界条件为p|边界=pi,其中pi为各边界瓦斯压力(Pa)。 (2)应力场定解条件: (1)初始条件为; 2 煤-瓦斯耦合模型 模拟建立的物理模型为三维立体模型,由于前文建立的数学模型无法利用软件中自带的简单模块进行求解,需要利用PDE,General form模型进行建模,通过调整偏微分方程的各项系数,将煤与瓦斯气固耦合关系、瓦斯气体的可压缩性、煤岩多孔介质中瓦斯气体的吸附、解吸特性等因素引入到方程的求解中。为了使所选用的PDE模块与Darcy定律模块在计算过程中相一致,在初始设置中,将Comsol Multiphysics PDE的因变量由原来表示位移因素的u修改为表示孔隙压力的p,这样在subdomain setting设置中使流固耦合控制方程和拟稳态扩散方程进行耦合,以便在求解器中进行耦合计算。 2.1 模型建立与监测 选取淮南矿业(集团)有限责任公司潘一矿11-2煤层掘进工作面进行不同掘进工艺条件下的煤岩瓦斯流固耦合效应的数值模拟研究。利用comsol软件建立的模型长×宽×高(x×y×z)=15 m×15 m×5 m,巷道长×宽×高(x×y×z)=10 m×2 m×2m,煤层初始孔隙压力1.7 MPa。在掘进面前方煤体距掘进面1 m和2 m处分别设1号与2号监测点。本文主要研究煤与瓦斯的耦合作用,而煤层顶底板为不透气层,故模拟中没有体现。几