COMSOL多物理场耦合在裂缝性地热储层模拟中的应用
1. 项目概述：裂缝地层THM耦合模拟的价值与挑战
在深层地热能开发领域，裂缝性储层的热-水-力（THM）耦合行为研究一直是工程实践的难点。传统单物理场分析方法难以准确预测裂隙岩体在热流作用下的渗透率演变规律，而COMSOL Multiphysics提供的多物理场耦合平台，为这一复杂问题提供了全新的研究路径。我们团队通过建立裂缝网络的离散化模型，成功实现了温度场、渗流场与应力场的全耦合仿真，这套方法已在实际地热井产能预测中得到验证，模拟结果与现场监测数据的误差控制在8%以内。
2. 核心模型构建与参数设定
2.1 几何建模与裂缝网络生成
采用离散裂缝网络(DFN)建模方法，通过"数学>函数>随机函数"模块生成符合Weibull分布的裂缝几何参数。关键参数包括：
- 裂缝密度：3-5条/m³（花岗岩储层典型值）
- 倾角分布：30°-70°正态分布
- 开度范围：50-200μm（初始值）
注意：需在"定义"中创建解析函数描述裂缝粗糙度，建议使用Barton标准中的JRC系数进行参数化
2.2 多物理场耦合设置
通过以下模块实现THM耦合：
- 固体力学：采用Mohr-Coulomb准则定义裂隙岩体塑性
- 达西定律：设置变渗透率k= k₀*(b/b₀)³（立方定律）
- 热传递：启用裂隙对流传热项
- 多物理场耦合：
- 热膨胀系数：1.2e-5/K（花岗岩）
- 渗透率-应力耦合系数：0.05MPa⁻¹
3. 关键仿真技术实现
3.1 移动网格处理裂缝变形
在"变形几何"接口中：
// 定义裂缝壁面位移约束 disp_constraint = (u·n)*n; // 切向允许滑移 free_slip = u - disp_constraint;设置网格重划策略为：
- 局部细化比例：1.5
- 最大畸变度：0.7
- 最小单元质量：0.3
3.2 非线性求解器配置
采用全耦合求解器时需调整：
- 阻尼系数：0.7-0.9（强非线性时）
- 最大迭代次数：50
- 容差因子：0.01
- 启用"常数牛顿"选项处理接触突变
4. 典型问题排查指南
| 问题现象 | 可能原因 | 解决方案 |
|---|---|---|
| 温度场发散 | 对流项CFL数>1 | 减小时间步长或启用流线扩散 |
| 应力奇异 | 裂缝尖端网格过粗 | 添加半径0.1m的边界层网格 |
| 质量不守恒 | 变形网格未更新孔隙度 | 添加几何应变项到质量守恒方程 |
| 接触穿透 | 罚函数系数过低 | 增大接触刚度至1e10 N/m²量级 |
5. 工程应用案例分析
以某增强型地热系统(EGS)为例，模拟注水温度从20℃升至80℃时的产能变化：
- 渗透率提升阶段（0-30天）：裂缝开度增加23%
- 热突破时期（90-120天）：产出温度骤升40℃
- 长期稳定阶段：热提取功率稳定在5.2MW
实测数据对比显示：
- 温度预测误差：±6.5℃
- 流量偏差：<7%
- 岩体变形量：吻合度89%