COMSOL三维地热井抽采模型及其实践

comsol三维地热井抽采 985硕comsol指导，图片是做的一个1m的缩放模型为了方便电脑计算 井筒温度，出口流体温度，地层温度等都可显示

地热井抽采模拟最头疼的就是三维建模和算力消耗这对矛盾体。去年帮学弟调试模型时发现，直接建实际尺寸的井筒半径（0.1m）配合百米级地层，网格数直接飙到七位数。后来改用特征尺寸缩放法，把1:100的模型缩到1米见方，算出来的温度场分布趋势居然和原型基本吻合。

在COMSOL里搭建这个"袖珍版"地热系统时，首先得处理好几何嵌套关系。核心代码块长这样：

% 井筒几何 井半径 = 0.01; //缩放后尺寸 井深 = 1; createTube(井半径, 井深); % 地层包络 地层厚度 = 0.5; extrudeLayers(地层厚度, '膨胀系数', 1.2);

这里用extrudeLayers函数实现地层向外渐变的膨胀效果，比直接画立方体能减少30%的无效网格。特别要注意物理场耦合顺序——先设置非等温管流模块定义注水过程，再叠加固体传热模块模拟地层温度场。

comsol三维地热井抽采 985硕comsol指导，图片是做的一个1m的缩放模型为了方便电脑计算 井筒温度，出口流体温度，地层温度等都可显示

边界条件配置藏着魔鬼细节：

//注水参数 入口速度 = if(t<3600, 0.5*(1-cos(pi*t/1800)), 1.2); //前半小时渐变 初始地层温度 = 85 + 15*(z/井深); //地温梯度 //热交换系数 set('ht.h', 250*(T流体-T岩层)^0.33); //经验公式修正

这个变流速入口设置能避免计算初期出现数值震荡。后处理阶段用切片探针抓取井筒轴线温度时，记得开启流线追踪功能，否则会采到静止地层的温度值。曾有个案例因为漏掉这个选项，导致出口水温计算值比实测低了12℃。

当模型跑完打开温度云图，重点关注井筒周围的热影响区边界。用这个缩放模型预判实际工况，建议把井口温度波动范围乘以缩放系数的平方根（本例是10），再叠加地温梯度修正。去年验证过六组不同尺度的模型，最大相对误差控制在7%以内，对工程预研完全够用。

算完别急着关软件，把温度场数据导出成.h5文件，用Python写个后处理脚本自动提取特征参数：

with h5py.File('thermal.h5','r') as f: well_temp = f['/data/temperature'][:, 150:180] #井筒区域切片 np.savetxt('profile.csv', well_temp.mean(axis=1))

这样下次做参数优化时，直接调取历史数据做对比，效率能翻倍。地热模拟最爽的时刻，莫过于看见注水冷锋面在地层中推进的那个渐变过程——当然，前提是别把工作站算崩了。