嵌入式离散裂缝模型的多尺度模拟有限差分模拟方法.pdfVIP

嵌入式离散裂缝模型的多尺度模拟有限差分模拟方法 张庆福，黄朝琴*，姚军，王月英，李阳 （中国石油大学（华东）石油工程学院，山东青岛，266580 ） 摘要：天然裂缝性油藏和人工压裂油藏内裂缝形态多样，分布复杂，传统的离散裂缝模型采用非结 构化网格进行网格划分，导致网格划分过程复杂，计算量大。嵌入式离散裂缝模型划分网格时不需 要考虑油藏内的裂缝形态，只需对基岩系统进行简单的网格剖分，可以大大降低网格划分的复杂度， 从而能够提高计算效率。然而，在油藏级别的数值模拟和人工压裂裂缝下的产能分析中，仍然存在 计算量巨大，模拟时间过长的问题。本文提出嵌入式离散裂缝模型的多尺度数值计算格式，使用多 尺度模拟有限差分法研究嵌入式离散裂缝模型渗流问题。该方法在粗网格上求解流动方程，因而计 算效率高，并通过在粗网格上求解局部流动问题计算多尺度基函数，多尺度基函数可以捕捉裂缝与 基岩间的相互关系，反映单元内的非均质性，因此这种方法在保证计算精度的同时极大地提高了计 算效率，数值结果表明这是一种有效的裂缝性油藏数值模拟方法。 关键词：裂缝性油藏，多尺度模拟有限差分法，嵌入式离散裂缝模型，油藏数值模拟 中图分类号：TE319 文献标识码：A 第一作者：张庆福，男，1990年5月生，研究方向为油藏数值模拟和裂缝性介质渗流理论。 通讯作者：黄朝琴（1981.8- ），男，博士，研究方向复杂介质渗流理论及油藏数值模拟，Email ：emc.group.upc@ 。 基金项目：国家自然科学基金51234007,）、山东省自然科学基金（ZR2014EEQ010 ）中央高校基本科研业务费专项资 金（编号：15CX05037A）资助项目、中国石油大学创新工程资助项目 （YCXJ2016015 ）。 1 引言 多孔介质中精确有效的流动模拟对于了解地下资源的动态，实施有效管理有着重要的意义，例 如地下水，地热能源及页岩气开发项目等。这些大规模的地质构造，除了固有的非均质性，往往还 包含复杂的天然裂缝。据不完全统计，世界上裂缝性油藏约占已探明总储量的一半；在我国，已探 [1] 明的裂缝性油藏地质储量超过40亿吨，占总探明储量的28% 以上 。而致密油气和页岩气等非常规油 [2] 气资源，需压裂后方可商业开采 ，压裂后也为裂缝性油气藏。 裂缝对流动形态有着重要影响，因此裂缝的描述直接影响到最终模拟结果的正确性。目前，主 [3,4] 要采用双重介质、等效连续介质和离散裂缝模型作为裂缝性介质的流动模型。双重介质模型 是目 前广泛应用的流动模型，当裂缝性油藏内裂缝发育程度较高时，其模拟效果较好，当存在数条控制 着流体流动方向和规模的大裂缝时，其计算结果误差较大。离散裂缝模型对介质中的裂缝予以显示 表征，由于其能够精确刻画流体在裂缝中的流动，该模型越来越受到人们的重视。现有的离散裂缝 数值模型都是基于匹配型网格，即把裂缝视为基岩内边界。由于裂缝分布复杂，往往基于非结构网 格划分，其剖分过程繁琐，计算量大。尤其当裂缝相距较近时网格剖分质量较差，导致计算误差。 对此，Lee[5-8]等提出了嵌入式离散裂缝模型，该模型对基岩直接进行结构化网格划分，然后将裂缝 嵌入基岩网格系统中，并根据裂缝与基岩的相交情况形成裂缝网格，避免了上述复杂的非结构化网 格剖分过程，大幅降低了计算复杂度。 但是在规模较大的裂缝性油藏中，裂缝具有强烈多尺度性，为了得到精细的解，我们需要进行 较精细的网格划分，即便采用嵌入式离散裂缝模型，仍然可达数百万甚至数亿个网格单元。如果基 于传统的数值计算方法求解此模型，巨大的计算量将超出当今计算机的计算能力。当前尺度升级法 [9]被广泛用来降低计算量，但是此方法无法充分利用小尺度精细数据信息，模拟精度不高。为此， [10-13] 需要一种既有尺度升级法的计算效率，又有较高计算精度方法，多尺度方法 应运而生。 多尺度方法在粗网格上求解控制方程，在细网格上求解局部流动方程获得多尺度基函数，以捕 捉小尺度流动特征