👨🎓 作者及单位
-
作者:Matteo Mario Picone(马泰奥·马里奥·皮科内)
-
学位:工程硕士(Master of Engineering)
-
专业:化学与石油工程
-
单位:加拿大卡尔加里大学(University of Calgary)
-
完成年份:2014年3月
📖 中文摘要(核心内容)
研究背景
油藏数值模拟是油藏工程中常用的预测工具,但受限于计算资源,精细地质模型(数百万网格)无法直接用于动态模拟,必须通过尺度升级(upscaling) 将粗化。在SAGD(蒸汽辅助重力泄油) 热采过程中,温度梯度大,需要精细网格才能准确追踪热前缘。然而,常规做法是垂直于井筒方向精细网格(1m),平行于井筒方向粗网格(25-100m)。但当储层存在非均质性(如页岩夹层)时,平行于井筒方向也会产生大温度梯度,导致粗网格模型与精细网格模型预测结果出现显著差异。
研究目标
本研究旨在:
-
量化网格尺寸和储层非均质性对SAGD模拟结果的影响
-
提出动态尺度升级方法,通过调整渗透率和热导率乘子,使粗网格模型匹配精细网格模型
-
建立基于SAGD产能指数(SPI) 的经验关联式,用于预测不同储层质量下的乘子值
研究方法
-
使用CMG STARS进行热采模拟
-
使用CMOST进行历史拟合,优化渗透率和热导率乘子
-
构建两类模型:
-
二元模型:砂岩(高渗)与页岩(零渗)的离散分布,页岩含量1-20%
-
岩相模型:更真实地模拟不同沉积环境(均匀、向上变粗、向上变细、河道)
-
-
网格尺寸变化:垂直井筒方向(i)2m或4m,平行井筒方向(j)25m、50m或100m
主要结论
-
均质模型中网格尺寸影响较小,差异主要体现在SAGD初期(蒸汽腔快速生长阶段)
-
非均质性对结果的影响远大于网格尺寸本身,相同页岩含量但不同分布位置,误差差异可达3-5%
-
动态尺度升级有效:通过调整渗透率和热导率乘子,粗网格模型与精细网格模型的累计产油、累计产水、累计注汽误差可降低约50-85%
-
SPI与乘子值存在良好相关性:
-
j方向(平行井筒) 乘子变化最大,与SPI相关性最强
-
i方向(垂直井筒) 乘子变化较小,可用平均值代替
-
-
计算效率显著提升:粗网格模型运行时间数分钟至数十分钟,精细网格模型数小时至数天数
🖥️ CMG软件应用情况总结
| 项目 | 说明 |
|---|---|
| 使用软件 | CMG STARS(热采模拟器)、CMG Builder(建模)、CMOST(历史拟合与优化)、Results(后处理) |
| 主要用途 | SAGD过程模拟、动态尺度升级参数优化、SPI计算 |
| 模型规模 | 粗网格:2,250网格(如4×100×1 m),精细网格:1×1×1 m |
| 模拟周期 | 10年(包含90天循环预热 + SAGD生产) |
| 优化方法 | CMG DECE(Designed Exploration Controlled Evolution)优化算法 |
| 目标函数 | 累计产油、累计产水、累计注汽的加权平均误差 |
| 优化参数 | PERM_i, PERM_j(渗透率乘子),TH_i, TH_j(热导率乘子) |
| SPI工具 | CMG内置的SAGD产能指数,用于评估地质模型的连通性 |
| 运行时间对比 | 粗网格:<1小时,精细网格:20-150小时 |
🧾 主要结论
网格尺寸与非均质性影响
| 模型类型 | 页岩含量 | 粗网格误差(4×100×1 m) | 精细网格误差(2×25×1 m) |
|---|---|---|---|
| 均质 | 0% | 4.77% | 1.10% |
| 非均质 | 10% | 8.46% | 3.41% |
-
非均质性使误差放大,且与页岩含量呈正相关
-
相同页岩含量但不同分布位置,误差差异可达3-5%
动态尺度升级效果
| 模型 | 原始误差 | 优化后误差 | 改进幅度 |
|---|---|---|---|
| 5%页岩(二元) | 7.43% | 1.82% | 75% |
| 10%页岩(二元) | 12.82% | 6.73% | 47% |
| 岩相模型(均匀) | 9.90% | 3.34% | 66% |
| 岩相模型(向上变细) | 8.11% | 2.66% | 67% |
SPI与乘子关联式
| 参数 | 关联式 | R² |
|---|---|---|
| PERM_j | -2.701×SPI + 8.724 | 较高 |
| TH_j | -6.274×SPI + 15.059 | 较高 |
| PERM_i | -0.185×SPI + 1.442 | 较低 |
| TH_i | 0.088×SPI + 0.999 | 较低 |
经济意义
-
单井对(1/3长度)误差10% → 36,750 m³沥青(900m全井)
-
50井对(5个井组)→ 183.75万立方米沥青的预测差异
🔮 未来工作建议
-
单参数敏感性分析:分别测试PERM和TH的独立影响
-
方向性参数:为每种岩相类型分配独立的渗透率/热导率乘子
-
对流换热研究:评估蒸汽腔内对流对热传递的贡献(目前认为<10%)
-
非均质位置敏感性:量化页岩夹层与井筒相对位置的影响
-
动态网格工具:对比动态尺度升级与动态网格加密的效果
📚 关键词(中文)
SAGD、尺度升级、动态尺度升级、非均质性、CMG STARS、CMOST、渗透率乘子、热导率乘子、SAGD产能指数、历史拟合、网格敏感性
