摘要: 利用煤层进行CO2地质封存成为了学术研究的热点之一。煤具有典型的双孔结构特征,在煤层中注入CO2既可有效地封存温室气体,还可驱替和提高煤层气产量,有广阔的应用前景。但是在煤层中封存CO2的过程中涉及到煤的变形、气体的吸附、渗流与扩散等一系列力学问题,极大地制约了这项技术的实施。根据多孔介质弹性力学理论,综合利用理论分析和数值模拟手段系统研究了CO2煤层封存中的煤体变形、气体吸附、气体渗流和气体扩散等多物理场耦合系统,主要得到以下结论:(1)根据裂隙和孔隙系统之间相互作用,建立了适合各自特点的渗透率变化模型。根据多孔介质弹性理论的有效应力原理建立了包含气体流动和吸附作用影响的煤的变形方程。根据气体在煤中的流动和吸附储存的特点分别建立了裂隙和孔隙中的气体流动控制方程。控制方程中包含了固体变形、流体渗流、流体吸附等多物理场耦合作用。(2)考虑到煤中残余水分对气体吸附的影响,利用Ettinger模型,建立了包含湿度效应的裂隙和孔隙系统的渗透率变化模型,并将湿度效应引入双重孔隙介质多物理场耦合模型。(3)利用COMSOL Multiphysics结合双重孔隙介质模型,对CO2煤层封存进行了模拟,并对相关影响因素进行了分析。煤层的裂隙发育程度直接决定了CO2注入的速度。煤层中裂隙越多,裂隙的初始渗透率越大,受到注气影响时变化越大。当水平方向上地应力的值不相等时,气体渗透具有方向性。气体在沿着最大应力方向上的渗透速度较快。当煤层中含有残余水分时,气体注入速度则更快,但是煤层对CO2的总封存量越小。(4)利用图像处理法和数值生成法研究了煤层中裂隙分布问题。图像处理法包括:数字图像处理技术和断层X射线扫描成像技术。前者通过处理平面图像获得初始裂隙的二维空间分布,后者通过一组CT扫描图片构造出岩样裂隙分布的三维模型。数值生成方法利用Weibull分布模型构建了岩样中裂隙的分布,讨论了Weibull均匀性指数对裂隙分布的影响。(5)将双重孔隙介质的多物理场耦合模型应用于CO2-ECBM中,建立了在二元气体竞争吸附过程中,裂隙的裂隙率和渗透率变化模型、孔隙中的孔隙率和渗透率变化模型。根据多孔介质弹性理论的有效应力原理和多元气体竞争吸附引起煤的变形关系,推导了煤层耦合变形方程。根据裂隙和孔隙中流体的质量守恒方程,动力扩散方程,Fick定律分别建立了包含二元气体竞争吸附效应的裂隙和孔隙中气体对流扩散耦合方程。(6)以沁水盆地煤层气增产试验为工程背景,将建立的CO2-ECBM模型应用于现场模拟,分别对不注气抽采煤层气、不同注气压力下抽采煤层气等情况进行了模拟研究。模拟结果表明,在煤层中注入CO2能有效地提高煤层气的采出率。