📝 摘要
超大型碳酸盐岩油田（如沙特阿拉伯的Ghawar油田和巴西盐下的Lula油田）表现出高度非均质性，这与薄层中的高渗透率区间密切相关。本文提出将局部网格加密与粗化流程相结合，通过在流动模拟中保留地质模型的静态属性和动态趋势，改善强层状油藏的模拟表征。研究分为五个主要步骤：(1) 定义常规粗网格；(2) 根据超K层和井位在常规粗网格中定义LGR；(3) 对所有方案应用粗化流程；(4) 直接在模拟模型中定义LGR，但不整合地质趋势；(5) 比较所有方案的动态响应。研究使用了基于巴西盐下和Ghawar油田信息构建的UNISIM-II-R基准模型。结果表明：将地质模型粗化为包含超K层LGR的粗网格，相比常规粗网格和仅在井周加密的LGR，能更准确地再现地质特征。对于常规粗网格或双重介质模拟中的LGR方案，应考虑拟相对渗透率曲线。直接在模拟模型中采用LGR的方法效果较差，因为将动态属性外推到粗块会偏离地质模型的动态行为。
🔍 关键词
局部网格加密；层状油藏；粗化；流动模拟；超K层；双重介质
🖥️ 模拟软件应用情况总结
| 项目 | 内容 |
|---|---|
| 所用软件 | 商业黑油模拟器（如ECLIPSE、CMG、PUMA FLOW） |
| 网格类型 | 结构化网格 |
| LGR类型 | 垂向加密（每个父网格分为3层） |
| 模拟对象 | 碳酸盐岩层状油藏（包含超K薄层） |
| 粗化方法 | 基于流动的有限差分法 |
| 比较对象 | 常规粗网格、LGR（井周/超K层/两者）、LGR SIM（常规模拟中直接加密） |
| 评价指标 | 累计产油量、累计产水量、NQDS（归一化平方距离带符号）、含水饱和度分布 |
文中明确指出：使用ECLIPSE、CMG或PUMA FLOW等商业黑油模拟器进行结构化网格的流动模拟，并采用LGR技术。
📊 研究方法与流程
| 步骤 | 内容 |
|---|---|
| 1 | 定义常规粗网格（网格尺寸100×100×8 m） |
| 2 | 根据超K层（渗透率>800 mD）和井位选择父网格 |
| 3 | 对所有方案应用粗化流程（流动法） |
| 4 | 直接在模拟模型中定义LGR（不整合地质趋势） |
| 5 | 比较所有方案与参考模型的动态响应 |
参考模型：UNISIM-II-R，网格尺寸50×50×2 m，模拟区域2000×2000×150 m。
井网：反五点井网（1口注水井，4口生产井）。
LGR方案：
-
仅井周加密
-
仅超K层加密
-
超K层+井周加密
🧪 主要结果与结论
| 结论 | 说明 |
|---|---|
| 网格角点对齐至关重要 | 若粗网格与细网格角点不匹配，LGR会失真，超K层渗透率被低估 |
| 超K层LGR效果最佳 | 井间超K层的LGR比井周LGR更显著改善动态匹配 |
| 相对渗透率外推问题 | 若将超K层的相对渗透率赋值给整个粗块，LGR效果接近于常规粗网格 |
| 拟相对渗透率必要性 | 对于常规粗网格或双重介质LGR，必须使用拟相对渗透率曲线 |
| 双重介质LGR局限性 | 商业软件无法对LGR块中的离散裂缝网络进行粗化；Oda方法仅适用于高连通裂缝网络 |
| LGR SIM方法效果差 | 常规方法（直接在模拟器中加密）因未整合地质趋势，动态匹配差 |
| 时间消耗 | LGR仅加密超K层时间适中（23秒单重介质，693秒双重介质），效果与超K+井周方案相同，效率更高 |
📈 定量评价（NQDS）
-
采用10%容差
-
Case A（网格角点不匹配）：LGR效果差，超出容差范围
-
Case B（网格角点匹配）：超K层LGR方案落在容差范围内
-
双重介质：所有模型更分散，LGR优势被掩盖
-
LGR SIM：全部超出容差范围，效果差于本文提出方法
✅ 主要结论
-
LGR必须与粗化流程整合，单纯在模拟器中加密无法反映薄层垂向非均质性。
-
超K层LGR是模拟层状油藏的最优方案，不建议仅加密井周。
-
必须为LGR中的加密块单独赋岩石类型和相对渗透率，否则效果等同于粗网格。
-
双重介质模型中，受限于裂缝粗化方法，LGR效果有限。
-
拟相对渗透率曲线对常规粗网格和双重介质LGR至关重要。
🏛️ 作者及单位信息
| 作者 | 单位 |
|---|---|
| Manuel Gomes Correia | 巴西坎皮纳斯大学（UNICAMP），石油工程系 |
| Célio Maschio | 巴西坎皮纳斯大学（UNICAMP），石油工程系 |
| Denis José Schiozer | 巴西坎皮纳斯大学（UNICAMP），石油工程系 |