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
关键提示:实际建模时需注意温度对岩石力学参数的非线性影响,建议采用温度相关的本构模型而非固定参数。
需要模型API调用? 免费领10W Token,多模型网关一键接入 Claude、DeepSeek 等主流模型。
2. COMSOL多物理场建模框架搭建
2.1 物理场接口选择策略
在COMSOL 6.1版本中,推荐使用以下模块组合:
- 固体力学(Solid Mechanics):处理岩石变形
- 达西定律(Darcy's Law):模拟裂隙流体流动
- 非等温流动(Nonisothermal Flow):耦合热对流/传导
- 断裂力学(Fracture Mechanics):可选模块用于显式模拟裂隙扩展
对于初学者,建议从简化模型入手:
- 先建立纯热流耦合(TH)模型验证热交换效率
- 再加入力学场进行完整THM耦合
- 最后引入化学场(THMC)考虑矿物溶解沉积效应
2.2 几何建模的特殊处理
干热岩储层的几何建模需要特别注意:
matlab复制% 典型几何层次结构
reservoir = createGeometry('cylinder', [0,0,0], 1000, 300)
