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.33.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%