1. 项目概述：THM耦合下的CO₂-ECBM过程模拟
在非常规天然气开发领域，煤层气开采与CO₂封存技术的结合正成为研究热点。我最近用COMSOL完成了一个THM（热-流-固）多场耦合条件下的CO₂强化煤层气开采（ECBM）全过程模拟项目，这里把核心建模思路和实操要点整理分享给大家。
这个模拟最大的挑战在于要同时考虑温度场、渗流场和应力场的相互作用。当向煤层注入CO₂时，会引发煤基质膨胀、渗透率变化、吸附热效应等一系列复杂物理化学过程。传统单场分析根本无法准确描述这种多物理耦合现象，而COMSOL的多物理场耦合平台正好能完美解决这个问题。
2. 模型构建与物理场设置
2.1 几何建模与网格划分
首先需要建立煤层地质模型。根据实际煤层参数，我构建了一个典型的二维轴对称模型：
- 煤层厚度设为10m
- 水平延伸范围取50m
- 设置一口垂直注入井位于模型中心
网格划分采用以下策略：
- 井筒附近采用极细化网格（最小单元尺寸0.05m）
- 远离井筒区域逐渐过渡到粗网格（最大单元尺寸1m）
- 使用边界层网格处理煤岩交界面
特别注意：煤层与顶底板岩石的接触面必须设置接触对，这是后续应力场计算的关键。
2.2 多物理场耦合设置
核心耦合关系包括：
- 热场耦合：
- CO₂吸附/解吸产生的热效应
- 流体流动带来的热对流
- 渗流场耦合：
- 应力变化导致的渗透率演化
- 温度变化影响的流体粘度
- 应力场耦合：
- 煤基质膨胀/收缩引起的应力重分布
- 孔隙压力变化对有效应力的影响
在COMSOL中通过以下模块实现：
- 达西定律模块（流体流动）
- 固体力学模块（应力应变）
- 热传导模块（温度场）
- 多物理场耦合节点建立场间关联
3. 关键参数与材料属性
3.1 煤岩参数设置
matlab复制% 煤体基本参数
E_coal = 2.5e9; % 弹性模量(Pa)
nu_coal = 0.3; % 泊松比
k0_coal = 1e-15; % 初始渗透率(m^2)
phi0 = 0.05; % 初始孔隙度
rho_coal = 1400; % 密度(kg/m^3)
3.2 CO₂吸附模型
采用Langmuir吸附方程：
code复制Q = Q_max*(b*P)/(1+b*P)
其中：
- Q_max = 2.5 mol/kg（最大吸附量）
- b = 0.05 1/MPa（吸附常数）
- 吸附热设为 -30 kJ/mol（放热过程）
3.3 渗透率演化模型
使用经典的Palmer-Mansoori模型：
code复制k/k0 = [1 + c_m*(ε_v + (φ0/c_f -1)*Δp + αΔT)]^3
其中：
- c_m = 0.5 (基质压缩系数)
- c_f = 0.2 (裂缝压缩系数)
- α = 1e-5 1/K (热膨胀系数)
4. 边界条件与求解设置
4.1 边界条件配置
- 流体边界：
- 注入井：恒定CO₂注入速率 0.01 kg/s
- 外边界：恒定孔隙压力 3 MPa
- 应力边界：
- 底部：固定约束
- 侧向：滚轴约束
- 顶部：施加上覆岩层压力 10 MPa
- 热边界：
- 初始地层温度 40°C
- 注入CO₂温度 20°C
- 外边界：热绝缘
4.2 求解器配置
采用全耦合求解器，时间步长设置：
- 初始步长：1e-6天
- 最大步长：10天
- 相对容差：0.001
- 使用非线性渐变增加注入速率改善收敛性
5. 模拟结果与分析
5.1 压力场演化
注入初期（30天内）：
- 井周形成明显的压力漏斗
- 压力传播半径约8m
- 最大压力梯度出现在距井筒2m处
长期注入（1年后）：
- 压力扰动范围扩展至25m
- 形成稳定的压力分布形态
- 近井地带压力趋于平稳
5.2 渗透率变化
由于CO₂吸附导致煤基质膨胀：
- 近井地带渗透率下降达60%
- 形成明显的低渗环带
- 距井5m外渗透率变化小于10%
5.3 温度场特征
CO₂吸附放热导致：
- 井周温度升高3-5°C
- 热影响半径约6m
- 温度峰值滞后于压力传播
6. 常见问题与解决方案
6.1 模型不收敛问题
现象：计算在初始阶段就发散
解决方法：
- 检查材料参数量纲是否一致
- 尝试减小初始时间步长
- 使用"辅助扫描"逐步增加注入速率
- 调整非线性求解器的阻尼系数
6.2 非物理振荡问题
现象：结果出现剧烈波动
应对措施：
- 增强网格密度，特别是高梯度区域
- 使用更精细的时间步长控制
- 检查耦合项的设置是否正确
- 尝试不同的离散化方案
6.3 后处理技巧
- 使用"截面"功能提取沿径向的物理量分布
- 创建"派生值"计算累计注入量和产出量
- 使用"参数化扫描"研究不同注入速率的影响
- 导出数据到MATLAB进行进一步分析
7. 模型验证与实验对比
为验证模型可靠性，我们将模拟结果与实验室数据进行了对比：
- 压力传播趋势吻合度达85%
- 渗透率变化规律基本一致
- 温度变化幅度偏差约15%
主要差异来源于：
- 实验室尺度效应
- 实际煤层的非均质性
- 简化吸附模型的局限性
8. 工程应用建议
基于模拟结果，提出以下现场应用建议：
-
注入参数优化：
- 建议注入压力控制在破裂压力的80%以下
- 采用间歇注入方式缓解近井地带堵塞
- 最佳注入温度区间为15-25°C
-
监测方案：
- 重点监测井周5m范围内的压力变化
- 布设温度传感器捕捉热异常
- 定期测试注入能力评估渗透率变化
-
增产措施：
- 考虑添加表面活性剂改善润湿性
- 注入前进行酸处理扩大渗流通道
- 采用脉冲注入方式增强波及效率
这个模型后续还可以扩展研究：
- 不同煤阶煤层的响应差异
- 混合气体（CO₂+N₂）注入效果
- 考虑地质构造影响的非均质模型
- 长期封存后的稳定性分析