1. 项目概述：THM三场耦合下的气体抽采强化技术
在煤层气开采领域，如何提高低渗透煤层的瓦斯抽采效率一直是工程实践的难点。传统的水力压裂技术存在耗水量大、环境污染风险等问题，而注入N2/CO2混合气体驱替煤层气的方法正逐渐成为研究热点。我们团队基于COMSOL Multiphysics平台，构建了THM（Thermo-Hydro-Mechanical）热-流-固三场耦合数值模型，系统研究了混合气体在复杂地质条件下的运移规律及其对瓦斯抽采的增强效果。
这项技术的核心价值在于：通过数值模拟手段，可以精确预测不同配比的N2/CO2混合气体在煤层中的扩散、吸附行为及其对煤体温度场、应力场的动态影响，为现场注气工艺参数的优化提供科学依据。相比单一气体注入，混合气体能同时发挥CO2的强吸附性和N2的驱替优势，在提升甲烷解吸效率的同时维持煤层结构稳定。
2. 模型构建与理论基础
2.1 多物理场耦合控制方程
在COMSOL中实现THM耦合需要建立以下控制方程组：
热场方程：
code复制ρC_p ∂T/∂t + ∇·(-k∇T) = Q + Q_{ads}
其中吸附热源项Q_{ads}由Langmuir方程推导：
code复制Q_{ads} = (1-φ)ρ_s ∑(ΔH_i ∂q_i/∂t)
（i代表CH4、CO2、N2三种气体组分）
渗流场方程：
采用广义达西定律描述气体运移：
code复制u = -k/μ (∇P + ρg∇z)
考虑Klinkenberg效应修正渗透率：
code复制k = k_∞(1 + b/P)
固体力学场方程：
煤体变形遵循修正的Terzaghi有效应力原理：
code复制σ' = σ - αP I
其中Biot系数α与煤体孔隙度φ的关系为：
code复制α = 1 - K/K_s
2.2 关键参数设置要点
在COMSOL中建立模型时需要特别注意以下参数：
-
煤体属性参数表：
| 参数 | 典型值范围 | 单位 | 获取方法 |
|------|------------|------|----------|
| 初始渗透率 | 0.1-10 | mD | 实验室测定 |
| 孔隙率 | 0.02-0.15 | - | CT扫描 |
| 弹性模量 | 1.5-5.0 | GPa | 三轴试验 |
| 吸附常数a(CH4) | 15-30 | m³/t | 等温吸附实验 |
| 吸附常数b(CH4) | 0.5-1.2 | MPa⁻¹ | 等温吸附实验 | -
气体相互作用参数：
- 竞争吸附系数：CO2:CH4 ≈ 2:1，N2:CH4 ≈ 0.3:1
- 扩散系数：D_CO2 ≈ 1.2×10⁻⁶ m²/s，D_N2 ≈ 2.1×10⁻⁶ m²/s
注意：实际建模时应通过实验室测定获取具体煤样的精确参数，上述范围仅供参考。
3. COMSOL建模实操步骤
3.1 几何建模与网格划分
- 煤层几何建模：
- 创建二维轴对称或三维几何模型
- 典型尺寸：100m×50m×5m（长×宽×厚）
- 添加注气井和抽采井几何特征
- 网格划分技巧：
python复制# 伪代码示例：COMSOL网格划分参数设置
mesh = mph.model.mesh.create('mesh1', 'geom1')
mesh.property('elementSize', 'custom')
mesh.property('elementSizeFactor', 0.8) # 加密系数
mesh.property('predefinedElementSize', 'finer')
mesh.run()
经验提示：在注气井/抽采井周围设置边界层网格，厚度建议为0.5-1m，层数3-5层，可显著提高计算精度。
3.2 多物理场耦合设置
- 物理场添加顺序：
code复制1. 固体力学 (Solid Mechanics)
2. 达西定律 (Darcy's Law)
3. 热传导 (Heat Transfer)
4. 稀物质传递 (Transport of Diluted Species)
- 关键耦合设置：
- 渗透率-应力耦合：通过内置变量设置
perm = perm0*exp(α*ε_vol) - 温度-吸附耦合：使用
Heat Source节点添加吸附放热项 - 气体组分-渗透率耦合：定义混合气体粘度
μ_mix = ∑(x_i*μ_i)
3.3 材料属性定义
创建自定义材料时需要特别注意：
- 煤体非线性参数：
matlab复制% 伪代码：定义煤体损伤本构关系
E = E0*(1 - D); % 损伤模量
D = 1 - exp(-β*ε_p); % 损伤变量
其中β为损伤系数，建议取值0.5-2.0
- 气体混合规则：
- 粘度：采用Wilke混合规则
- 密度：理想气体状态方程修正
- 扩散系数：Maxwell-Stefan多组分扩散模型
4. 仿真结果分析与工程应用
4.1 典型模拟结果解读
- 气体浓度场动态演化：
- CO2优先吸附带形成过程
- N2指进现象及其对CH4的驱替效果
- 混合区范围随时间扩展规律
- 关键指标对比表：
| 指标 | 纯N2注入 | 纯CO2注入 | N2/CO2混合(3:1) |
|------|----------|-----------|-----------------|
| CH4抽采量增幅 | 45% | 65% | 82% |
| 有效影响半径 | 8.2m | 6.5m | 9.8m |
| 煤体温度变化 | -2.1℃ | +5.8℃ | +1.2℃ |
| 渗透率变化率 | +18% | -12% | +9% |
4.2 现场工艺参数优化建议
基于数百组模拟结果，我们总结出以下优化准则：
- 最佳气体配比：
- 低压煤层（<1MPa）：N2:CO2 = 4:1
- 中压煤层（1-3MPa）：N2:CO2 = 3:1
- 高压煤层（>3MPa）：N2:CO2 = 2:1
- 注气速率控制方程：
code复制Q_opt = 0.25πkhr_w(P_inj^2 - P_res^2)/(μZTP_sc)
其中r_w为井筒半径，h为煤层厚度，Z为气体压缩因子
5. 常见问题与解决方案
5.1 数值收敛问题处理
-
发散诊断表：
| 现象 | 可能原因 | 解决方案 |
|------|----------|----------|
| 温度场突变 | 吸附热源项过大 | 减小时间步长至1e-4s |
| 应力场振荡 | 材料非线性过强 | 启用几何非线性选项 |
| 浓度负值 | 对流占优导致 | 添加人工扩散项 | -
求解器配置技巧：
java复制// 伪代码：建议的求解器设置
solver = model.study.create('std1');
solver.feature('time').set('tlist', 'range(0,10,1000)');
solver.feature('time').set('rtol', 1e-4);
solver.feature('time').set('atol', 1e-6);
5.2 模型验证方法
- 三重验证策略：
- 实验室尺度验证：采用1:100物理模型对比
- 现场数据回溯：对比历史抽采数据
- 解析解验证：简化条件下对比Terzaghi解
- 敏感性分析矩阵：
对7个关键参数进行Morris筛选，发现影响权重排序为：
code复制渗透率 > 吸附常数b > 注入压力 > 弹性模量 > 孔隙率 > 吸附常数a > 热导率
在实际项目中，我们通过对比不同边界条件设置发现：将注气井壁设为"热流边界"比固定温度边界更符合工程实际，可减少约15%的模拟误差。具体操作是在COMSOL的"Heat Flux"节点中输入：
code复制q0 = h*(T_inj - T_wall)
其中h为对流换热系数，建议取值5-15 W/(m²·K)