COMSOL多物理场耦合在干热岩地热开发中的应用

1. 干热岩开采的工程挑战与技术背景

干热岩（Hot Dry Rock）作为地热资源开发的前沿领域，其储层温度通常在150-650℃之间，埋深在3-10公里范围内。与传统地热系统不同，干热岩缺乏天然流体作为传热介质，需要通过人工压裂形成裂隙网络，并注入工作流体进行热交换。这一过程涉及三个核心物理场的强耦合作用：

• 热力学（Thermal）：岩石与流体的热量传递
• 水力学（Hydraulic）：裂隙网络中的流体流动
• 力学（Mechanical）：岩体应力变形与裂隙扩展

在COMSOL中实现这类THM（Thermo-Hydro-Mechanical）耦合问题，需要特别关注几个关键参数：

python复制# 典型干热岩参数范围
rock_density = 2600-2800  # kg/m³
thermal_conductivity = 2.0-3.5  # W/(m·K)
specific_heat = 800-1200  # J/(kg·K)
permeability = 1e-18-1e-15  # m² (压裂后)
injection_rate = 10-50  # kg/s

关键提示：实际建模时需注意温度对岩石力学参数的非线性影响，建议采用温度相关的本构模型而非固定参数。

2. COMSOL多物理场建模框架搭建

2.1 物理场接口选择策略

在COMSOL 6.1版本中，推荐使用以下模块组合：

1. 固体力学（Solid Mechanics）：处理岩石变形
2. 达西定律（Darcy's Law）：模拟裂隙流体流动
3. 非等温流动（Nonisothermal Flow）：耦合热对流/传导
4. 断裂力学（Fracture Mechanics）：可选模块用于显式模拟裂隙扩展

对于初学者，建议从简化模型入手：

• 先建立纯热流耦合（TH）模型验证热交换效率
• 再加入力学场进行完整THM耦合
• 最后引入化学场（THMC）考虑矿物溶解沉积效应

2.2 几何建模的特殊处理

干热岩储层的几何建模需要特别注意：

matlab复制% 典型几何层次结构
reservoir = createGeometry('cylinder', [0,0,0], 1000, 300); % 储层主体
fractures = generateStochasticFractures(50, 10-100); % 随机裂隙网络
injection_well = createPipe([0,0,-3000], [0,0,0], 0.15); % 注入井
production_well = createPipe([500,0,-3000], [500,0,0], 0.15); % 生产井

实际操作中建议：

1. 使用CAD导入功能处理复杂地质结构
2. 对裂隙网络采用离散裂缝网络（DFN）建模
3. 通过"形成装配体"处理井筒-储层接触关系

3. 材料属性与边界条件设置

3.1 岩石本构模型选择

花岗岩类干热岩推荐使用：

• 热弹性模型（线性阶段）
• Drucker-Prager塑性准则（非线性阶段）
• 温度相关的杨氏模量：

  code复制E(T) = E0*(1 - 0.5*(T-20)/800)  # T in ℃
  

流体属性设置要点：

• 使用"非牛顿流体"接口处理高温高压下水的性质变化
• 考虑超临界状态（>374℃, >22.1MPa）的特殊物性
• 注入流体建议采用TABULATED格式输入实验数据

3.2 关键边界条件配置

典型边界条件设置：

特别注意：生产井边界建议使用"压力约束+回流防止"组合条件，避免数值振荡。

4. 求解器配置与计算优化

4.1 多物理场耦合策略

推荐采用分步耦合方案：

1. 稳态初始应力场（仅力学）
2. 瞬态TH耦合（固定裂隙几何）
3. 全耦合THM分析（启用几何非线性）

求解器配置技巧：

java复制// 典型求解器序列
solverSequence = {
  "Stationary",  // 初始应力
  "TimeDependent1",  // TH耦合
  "TimeDependent2"   // THM全耦合
};
solverConfig = {
  "relativeTolerance": 0.01,
  "maxIterations": 50,
  "dampingFactor": 0.8  // 非线性问题必备
};

4.2 计算资源管理

针对大规模模型建议：

• 使用"集群扫描"功能进行参数化研究
• 启用"几何多重网格"加速收敛
• 临时文件存储设置：
  • Windows: C:\Users\[user]\AppData\Local\Temp\COMSOL
  • Linux: /tmp/comsol_[pid]

内存优化技巧：

1. 对对称模型使用"对称简化"
2. 裂隙网络采用"降维建模"（2D裂缝）
3. 输出时启用"时间步选择"减少存储量

5. 后处理与结果验证

5.1 关键结果可视化

必看结果项：

• 温度场时空演化（动画）
• 裂隙孔径分布云图
• 热提取率随时间变化曲线
• 储层阻抗因子（阻抗=ΔP/流量）

自定义派生值示例：

code复制热回收效率 = (T_inject - T_produce)/(T_inject - T_initial)
应力阴影区 = where(σ3 > 0.9*σ1)  // 压应力集中区

5.2 模型验证方法

建议验证步骤：

1. 与解析解对比（如Theis方程）
2. 网格独立性检验（3种网格尺寸）
3. 参数敏感性分析（Morris筛选法）
4. 现场数据对比（如有）

常见问题排查：

• 温度场异常：检查热容单位是否统一
• 不收敛问题：调整阻尼因子或时间步长
• 质量不守恒：检查裂隙连通性设置

6. 工程案例：增强型地热系统仿真

某EGS项目参数：

• 储层体积：2.5 km³
• 初始温度：235℃
• 注入流量：30 kg/s
• 运行周期：30年

模拟结果要点：

• 热突破时间：第8年
• 平均热提取率：5.2 MW
• 裂隙网络扩展：主要沿NE70°方向

优化建议：

1. 采用周期性注采策略延缓热突破
2. 在储层边缘布置补充注入井
3. 控制注入压力在临界值以下防诱发地震

在长期运行中观察到有趣的"热虹吸"现象：高温流体自动向生产井聚集，这与传统的达西流动预测存在偏差。通过引入Forchheimer方程修正，模型预测精度提高了18%。这个案例说明，干热岩模拟需要持续跟踪现场数据并迭代模型。