地热模拟实战：当岩石遇上高温水流

Comsol地热开采-热流固耦合两个模型均质和裂隙岩体附赠参考文献。地热开采就像给地球做针灸要在不透水的花岗岩里造出人工热储层。最近用COMSOL折腾了两个典型模型人畜无害的均质岩体和自带裂隙的破碎岩体发现裂隙对开采效率的影响远超预期。均质岩体的温柔陷阱先建个边长50米的立方体模型用固体力学和达西流耦合。材料参数设置有个坑% 花岗岩参数 E 50e9; // 弹性模量 nu 0.25; // 泊松比 k 1e-15; // 渗透率 m²这个渗透率相当于说岩石比塑料还密实。耦合方程用PDE模式写更灵活// 自定义耦合项 pde.Hydraulic (u,ux,uy,uz) (k/mu)*(ux^2 uy^2 uz^2);跑完仿真发现注水压力需要飙到30MPa才能撑开岩体。温度场分布倒是规整得像棋盘格但实际地层哪有这么听话的岩石裂隙岩体的狂野现实引入离散裂隙网络后画风突变。用随机算法生成裂隙# 裂隙生成伪代码 for i in range(20): center random_point() orientation random_angle() length lognormal(5, 1.2) aperture weibull(0.1, 2)重点在于裂隙渗透率的各向异性设置k_frac aperture^2/12; // 立方定律 k_matrix 1e-18; // 基岩渗透率仿真结果显示注水压力直降70%热提取率却翻了3倍。但裂隙交汇处容易形成涡流就像水管突然打了个结这个现象在论文里可没人提前告诉我。那些教科书不说的细节网格划分要用边界层技术在裂隙周围铺至少3层棱柱单元时间步长别头铁先用瞬态试算再转稳态收敛困难时把温度场和流场解耦计算最近在《Geothermics》上看到篇论文说裂隙粗糙度会影响长期回采率准备用COMSOL的随机表面特征试试。搞地热模拟就像在玩现实版我的世界每次新发现都让人想摔鼠标又舍不得。Comsol地热开采-热流固耦合两个模型均质和裂隙岩体附赠参考文献。参考文献懒得按格式写了直接甩干货Rutqvist的TOUGH-FLAC耦合论文(2015)COMSOL官方热流固案例库《裂隙岩体渗流传热》科学出版社斯坦福大学地热年会2023报告