COMSOL案例：非均质储层地热能群井抽采的实践

COMSOL案例非均质储层的地热能群井抽采 适用于做地热能开采模型为非均质模型利用地质统计学模拟得到储层的非均质性加载到comsol中表征渗透率的非均质性。地热能开采最头疼的就是储层的非均质性。想象一下你打下去的井可能同时遇到透水层和隔水层就像在巧克力脆皮冰淇淋里找果仁——完全看运气。这时候地质统计学建模就派上用场了它能帮我们把随机分布的渗透率场具象化。用Python生成渗透率场的时候我习惯用高斯随机场模拟。举个栗子import numpy as np from scipy.ndimage import gaussian_filter np.random.seed(114514) # 玄学种子 raw_field np.random.randn(100,100) smooth_field gaussian_filter(raw_field, sigma5) perm_field 1e-13 * np.exp(3*smooth_field) # 渗透率对数正态分布这段代码里有个小技巧高斯滤波的sigma参数控制着渗透率场的关联长度实际操作中得结合现场地质勘探数据调整。导出的txt文件可以直接喂给COMSOL的插值函数注意坐标系对齐别翻车。COMSOL案例非均质储层的地热能群井抽采 适用于做地热能开采模型为非均质模型利用地质统计学模拟得到储层的非均质性加载到comsol中表征渗透率的非均质性。在COMSOL里设置多物理场耦合时建议先单独验证达西流模块。见过有人直接把地质模型导入就开跑结果算出来的压力场像抽象画——记得检查渗透率张量的方向性设置。群井布局建议用参数化扫描来优化特别是注采井间距这个参数对热突破时间影响能差出数量级。当看到后处理中出现这种温度云图配图示意说明热锋面推进出现了各向异性。这时候该翻回去检查地质统计模型的变差函数设置可能水平方向的渗透率相关性设得过高了。有个经验公式最优注采井距≈0.8×热影响半径×sqrt(渗透率各向异性比)不过具体项目还是得做参数敏感性分析。最后唠叨一句非均质模型特别吃网格质量建议在渗透率突变区域手动加密网格。见过有人用自动生成网格跑出负渗透率的笑话这种时候就该祭出comsol的边界层网格功能了。