1. 项目背景与核心目标
这个模型聚焦于非常规天然气开采中的关键科学问题——二氧化碳驱替煤层气过程中的多物理场耦合现象。在深层煤层气开发中,传统水力压裂技术面临水资源消耗大、环境污染等挑战,而CO2-ECBM(二氧化碳增强煤层气采收)技术既能提高甲烷产量,又能实现碳封存,成为近年来的研究热点。
我们团队通过COMSOL Multiphysics构建的热-流-固全耦合模型,能够精确模拟:
- 注入超临界CO2与吸附态CH4的竞争吸附过程
- 煤体基质收缩/膨胀引发的渗透率动态变化
- 温度场变化对吸附/解吸平衡的影响
- 多场耦合作用下的气体运移规律
2. 模型构建的关键技术路线
2.1 多物理场耦合机制解析
模型包含三个核心耦合模块:
- 热力学模块:采用非等温达西定律描述气体流动,考虑焦耳-汤姆逊效应导致的温度变化
- 化学吸附模块:基于扩展的Langmuir方程,建立CO2/CH4二元竞争吸附模型
- 力学模块:引入Palmer-Mansoori方程表征有效应力与基质应变的关系
关键参数:CO2吸附热(~32 kJ/mol)显著高于CH4(~16 kJ/mol),这是驱替效率差异的主因
2.2 几何建模与网格划分技巧
采用二维轴对称模型简化计算:
- 煤层厚度设定为5m(典型中阶煤参数)
- 注气井与生产井间距100m
- 使用极细化边界层网格(最小单元尺寸0.01m)捕捉近井地带的高梯度变化
- 网格独立性验证显示当单元数超过15万时,结果收敛误差<2%
3. 材料属性与边界条件设置
3.1 煤岩参数实验测定
通过实验室测试获取关键参数:
| 参数 |
数值 |
测试方法 |
| 初始孔隙度 |
0.05-0.12 |
氦气孔隙度仪 |
| 初始渗透率 |
0.1-1mD |
脉冲衰减法 |
| 弹性模量 |
2-5GPa |
三轴压缩试验 |
| 吸附常数b(CO2) |
0.85MPa⁻¹ |
高压吸附仪 |
3.2 边界条件实现方案
-
注入端:设置混合气体质量流量边界
- CO2纯度>99%
- 注入压力梯度8-12MPa/km
- 温度保持318K(典型储层条件)
-
生产端:采用压力边界条件
4. 求解器配置与计算优化
4.1 多物理场耦合求解策略
采用分离式求解器分步迭代:
- 先求解流体流动与传热(瞬态分析)
- 将孔隙压力场映射到固体力学接口
- 更新渗透率张量后返回第一步迭代
经验提示:设置阻尼因子0.7可有效改善强非线性问题的收敛性
4.2 高性能计算配置
- 使用MUMPS直接求解器处理刚度矩阵
- 在64核服务器上并行计算,典型案例耗时约6-8小时
- 存储方案:每10个时间步保存一次完整解,关键变量(如渗透率)每步保存
5. 典型模拟结果与分析
5.1 驱替前沿演化特征
通过浓度场云图可见:
- CO2突破时间约45天(注入速率5t/d)
- 驱替效率呈现明显各向异性(水平/垂向波及系数比≈3:1)
- 甲烷产量曲线出现双峰特征(对应基质孔隙与裂隙系统产出)
5.2 渗透率动态变化规律
渗透率演化呈现三阶段特征:
- 初始上升期(0-30天):CH4解吸导致基质收缩
- 快速下降期(30-90天):CO2吸附引发膨胀
- 稳定期(>90天):裂隙系统重新压实
6. 现场应用验证与参数反演
6.1 山西沁水盆地案例对比
将模拟结果与现场监测数据对比:
| 指标 |
模拟值 |
实测值 |
误差 |
| 累计产气量(1年) |
4.2×10⁶m³ |
3.8×10⁶m³ |
+10.5% |
| 井口CO2突破时间 |
68天 |
73天 |
-6.8% |
| 储层压力降幅 |
34% |
31% |
+9.7% |
6.2 敏感性分析关键发现
通过Morris筛选法识别出主导因素:
- 初始裂隙渗透率(灵敏度指数0.62)
- CO2吸附常数(0.51)
- 弹性模量(0.39)
7. 工程优化建议与扩展应用
基于模拟结果提出改进方案:
- 采用脉冲注气模式(注3天停2天)可提高波及效率18%
- 注入前酸处理能改善近井地带渗透性
- 结合微震监测优化注采井网布置
该模型框架经适当修改后,还可应用于:
- 页岩气藏CO2驱替开发
- 天然气水合物开采模拟
- 地热储层改造评估