1. 项目背景与工程意义
在非常规天然气开采领域，注气驱替技术已成为提高甲烷采收率的重要手段。这项技术通过向储层注入二氧化碳或氮气等气体，利用压力驱动和组分置换效应将吸附在煤岩或页岩中的甲烷置换出来。然而传统数值模拟往往只考虑流体流动（Flow）与组分传输（Mass Transfer），忽略了热力学（Thermal）与固体力学（Solid Mechanics）的耦合效应，导致预测结果与实际工况存在显著偏差。
THM（Thermo-Hydro-Mechanical）耦合分析的核心价值在于：
- 温度变化影响气体吸附/解吸平衡（Langmuir方程中的参数对温度敏感）
- 注气过程引发的热应力会改变储层孔隙结构和渗透率
- 流体压力与固体变形存在双向耦合关系（有效应力原理）
- 多场耦合效应共同决定了甲烷解吸-扩散-渗流的全过程
以某煤层气田的实际数据为例，仅考虑流固耦合（HM）时预测的甲烷采收率比THM全耦合模型高估约18%，这种偏差在长期注气项目中会累积成巨大的经济差异。Comsol Multiphysics凭借其真正的多物理场耦合求解能力，成为解决此类复杂问题的理想工具。
2. 模型构建关键步骤
2.1 几何建模与材料定义
采用二维轴对称模型简化实际井筒-储层系统：
matlab复制% COMSOL几何脚本示例
geom = ModelUtil.create('Model');
geom.feature.create('cyl', 'Cylinder');
geom.feature('cyl').set('r', '5[m]'); % 储层半径
geom.feature('cyl').set('h', '50[m]'); % 储层高度
geom.feature('cyl').set('pos', [0,0]);
geom.run;
材料属性定义需特别注意：
- 煤岩：采用Drucker-Prager塑性模型，渗透率设为应变张量的指数函数
- 甲烷：使用Peng-Robinson状态方程，考虑临界温度(Tc=190.56K)和压力(pc=4.59MPa)
- 注入气体(CO₂)：设置与甲烷不同的吸附特性（吸附热差约10kJ/mol）
2.2 物理场耦合设置
核心耦合关系通过PDE接口实现：
-
热场耦合：
- 能量方程包含对流项(ρCp*u·∇T)和热源项（吸附热Q_ads）
- 热膨胀系数设为各向异性（垂向α_v=1.2e-5/K，水平α_h=8e-6/K）
-
流场耦合：
math复制\frac{\partial(φρ)}{\partial t} + ∇·(ρu) = Q_m其中孔隙率φ与体积应变ε_v的关系：
math复制φ = φ_0 + (1-φ_0)(ε_v - ε_{v0}) -
固体力学耦合：
采用增量形式的本构关系：math复制Δσ = C:(Δε - αΔTI - Δε_p)其中塑性应变增量Δε_p通过DP准则计算
2.3 边界条件设置
典型边界条件配置：
| 边界类型 | 热场 | 流场 | 固体力学 |
|---|---|---|---|
| 注入井壁 | 恒温T_inj | 流量边界q_inj | 固定位移 |
| 生产井壁 | 热对流系数h | 压力边界p_prod | 自由位移 |
| 外边界 | 地温梯度 | 无流动 | 弹性支撑 |
特别注意：移动网格（ALE）接口需启用几何非线性选项，并设置适当的网格平滑参数（Laplacian系数建议0.3-0.7）
3. 求解器配置技巧
3.1 多物理场求解策略
推荐采用分离式求解器（Segregated Solver）分步迭代：
- 先求解热-流耦合（瞬态分析，时间步长自适应）
- 将热流结果作为载荷传递到固体力学（准静态分析）
- 更新几何和材料参数后返回第1步
关键参数设置：
matlab复制solver = ModelUtil.create('SolverConfig');
solver.set('tlist', 'range(0,10,3600)'); % 1小时步长
solver.set('rtol', 1e-4);
solver.set('atol', 1e-6);
solver.set('maxiter', 50);
3.2 非线性收敛控制
常见问题及解决方案：
- 振荡发散：启用阻尼因子（建议初始值0.7）
- 收敛慢：采用增量加载（分阶段施加注气压力）
- 网格畸变：启用自动重划分（remeshing）阈值设为0.85
实测表明，当注入压力超过8MPa时，需要将固体力学模块的几何非线性选项从"小应变"切换为"有限应变"。
4. 后处理与结果分析
4.1 关键指标提取
创建以下派生变量监控过程：
- 甲烷采收率：
math复制η = \frac{\int Γ_{prod} c_{CH4} u·n dΓ}{\int Ω_0 ρ_{CH4} dΩ} - 储层平均渗透率变化：
math复制k_{avg} = \frac{1}{V}\int k(ε_v) dV - 热影响半径：
math复制R_T = max(||x||), where |T(x)-T_0|>5K
4.2 典型结果可视化
- 流线图叠加温度云图（显示热锋面推进）
- 应力张量的第一不变量分布（识别潜在裂缝区域）
- 甲烷摩尔分数动画（观察置换前沿形态）
某案例的模拟结果显示：
- 注气30天后热影响半径达12.7m
- 近井地带渗透率提升约40%（由于剪切膨胀效应）
- 甲烷解吸速率在50℃时比30℃高2.3倍
5. 模型验证与实验对比
采用某煤层气井的现场数据进行验证：
| 参数 | 实测值 | 模拟值 | 相对误差 |
|---|---|---|---|
| 日产气量(m³/d) | 2350 | 2210 | 6.0% |
| 井底温度(℃) | 48.2 | 46.7 | 3.1% |
| 套压(MPa) | 3.2 | 3.4 | 6.3% |
误差主要来源于：
- 实际储层的非均质性未在模型中完全体现
- 忽略了地层水的相变影响
- 边界条件简化带来的系统误差
6. 工程应用建议
基于数百次模拟的经验总结：
-
注气参数优化：
- 最佳注入温度区间：50-70℃（过高会导致煤基质收缩过度）
- 压力提升速率控制：<0.5MPa/h（避免诱发剪切破坏）
-
监测方案设计：
- 在径向距离5m、10m处布置温度-压力传感器
- 定期进行声波测井监测体积应变
-
经济性评估：
- 当采收率提升<8%时，THM耦合分析可能不具经济性
- 对于埋深>800m的储层，热应力贡献可达总应变的30%
实际工程中曾遇到一个典型案例：某项目初期仅进行HM分析，导致注气井间距设计过大。后期THM耦合分析发现热传导效应使有效驱替半径缩小15%，及时调整井距后单井产量提高22%。