一种基于角点网格的高精度油气运移模拟方法与流程

本发明涉及盆地模拟和油气成藏过程模拟技术领域，特别是勘探程度较高的盆地级、凹陷级不同层次油气系统的一种基于角点网格的高精度油气运移模拟方法。

背景技术：

随着油气藏勘探开发的不断深入，对基于高精度地质模型下的油气成藏模拟提出了更高的要求。传统的直角正交网格模型在油气数值模拟过程中应用较为广泛，其特点是计算速度快，其所有网格的长宽相等，垂向连接网格面垂直。但是在有断层情况下，这种方式就不能建立精细的地质模型。角点网格模型克服了正交网格模型的缺点，可以根据控制区域的长度和宽度灵活调整各油层的网格步长，从而使断层附近的网格边线与断层的走向重合。该技术能够更精准的完成油气运移模拟，为油藏开发提供更好的分析工具。更精准的角点网格模型带来了更精确的模拟刻画，同时也模拟算法实现带来了更大的难度。目前的各类模拟软件，仅能够实现传统正交网格模型的模拟方法，如何在更精准的角点网格模型下，处理好构造作用、沉积作用、地热作用、生烃作用、排烃作用及运烃作用、聚烃作用之间的关系，将该技术运用到油藏开发中，成为了目前亟待解决的问题。

技术实现要素：

本发明的目的是为了解决以上技术问题，提供一种基于角点网格的高精度油气运移模拟方法。

本发明的目的可通过如下技术措施来实现：

一种基于角点网格的高精度油气运移模拟方法，包括：

步骤1，建立基于角点网格的地质模型和属性模型，进行前期的热史、生排烃模拟，形成模拟成果数据；

步骤2，首先计算每个单元体的中心坐标，加密体处理；其次是进行运移量分配评估；第三是进行油气运移计算；第四是进行初次运移、查找通道，将当前通道的油气运移达到体平衡；

其中：在步骤1中，对于地质模型和属性模型，对于研究工区的每套对应地层、断层都要建立，网格数统一；断层要根据实际情况进行封堵性设置，属性模型信息至少应包括岩性、孔隙度、渗透率；在步骤2中，对于局部加密的每个加密体，每个小网格要进行渗透率、孔隙度处理，然后统计到大网格三个方向的渗透率和综合孔隙度。

上述方案进一步包括：

进行前期的热史、生排烃模拟，要形成标准格式的温度文件、超压文件、生排烃文件、油气流体势文件。

在步骤2中，通过达西定律进行运移量分配评估，计算每个单元格的油粘度、油密度、液体深度压力、XYZ三个方向的流体速度。

在步骤2中，进行油气运移计算，根据地质体属性计算每个网格能容纳的烃量，计算每个网格向相邻网格的运移情况。

在步骤2中，单元格间运移，要对于每个单元格的各个方向，根据岩性计算网格间接触面积、排替压力，进行断层封堵性判断。

本发明中的基于角点网格的高精度油气运移模拟方法，采用了能够较为准确描述裂缝空间展布和特征的角点网格模型，同时运用了动态的体平衡技术，从而提高了运移模拟精度，为勘探开发提供有效依据。

附图说明

图1为本发明的一种具体实施例的流程简图；

图2是本发明的角点网格模型的应用效果图；

图3是本发明的油气运移模拟应用效果图。

具体实施方式

为使本发明的上述和其他目的、特征和优点能更明显易懂，通过典型实施例，作详细说明如下。

基于角点网格的高精度油气运移模拟方法包括：

步骤1，建立基于角点网格的地质模型和属性模型，进行运移前的热史、生烃、排烃等模拟，逐期次获取每个期次的各类成果数据（包括体数据、断层信息、生排烃信息、温度信息、超压信息、油气流体势等）和属性数据（包括岩性、孔隙度、渗透率等）；步骤2，对每个期次进行油运聚模拟、加密体处理、运移量分配、单元格间运移等运算处理，最后输出油气运移强度、运移方向和运移量。

在步骤1中，要首先建立基于角点网格的地质模型和属性模型，对于地质模型和属性模型，对于研究工区的每套对应地层、断层都要建立，网格数统一，局部可进行加密处理。断层要根据实际情况进行封堵性设置，属性信息至少应包括岩性、孔隙度、渗透率等。对建好的地质模型和属性模型，进行前期的热史、生排烃模拟，形成标准格式的温度文件、超压文件、生排烃文件、油气流体势文件，以便于在油气运移模拟的属性提取。

在步骤2中，首先计算每个单元体的中心坐标，对每个加密体的小网格处理，然后统计到大网格三个方向的渗透率和综合孔隙度；其次是通过达西定律进行运移量分配评估，计算每个单元格的油粘度、油密度、液体深度压力、XYZ三个方向的流体速度；第三是进行油气运移，根据地质体属性计算每个网格能容纳的烃量，计算每个网格可以向相邻网格的运移情况；第四是进行初次运移、查找通道，将当前通道的油气运移达到体平衡。

下面结合附图1对做进一步说明。

在步骤101，通过收集建模所需各项基础数据，应用模型建立工具，建立基于角点网格的地质模型和属性模型。流程进入到步骤102。

在步骤102，通过热史、生排烃方法等，进行运移模拟前的模拟处理，形成各期次的成果文件。流程进入到步骤103。

在步骤103，根据期次循环处理以下各步骤。获取每个期次的各类成果数据和属性数据。成果数据包括体数据、断层信息、生排烃信息、温度信息、超压信息、油气流体势等；属性数据包括岩性、孔隙度、渗透率等。流程进入到步骤104。

在步骤104，对于有局部加密的数据体，要进行加密体处理。对每个加密体，获取加密网格逻辑单元起止序号，计算XYZ方向被加密网格数，提取加密体属性（岩性、孔隙度、渗透率等），对每个加密网格进行处理，统计到大网格三个方向的渗透率及综合孔隙度。流程进入到步骤105。

在步骤105，进行运移量分配估算，对每个单元格分别计算油粘度、油密度、流体深度压力、XYZ方向的流体速度。流程进入到步骤106。

在步骤106，计算油气运移。根据地质体属性计算网格能容纳的烃量，计算每个网格向相邻单元格的运移情况，然后进行初次运移、查找通道，最终将当前通道中的油气运移到体平衡，输出该期次的油气运移强度、运移量、运移方向文件等。

其中，在单元格间运移时，首先获取角点网格和断层文件，其次依次计算向上和运移的其他方向，根据岩性，计算网格间接触面积、排替压力，判断断层封堵性，得到单元格间的运移关系。

图2是车西地区基于角点网格的现阶段三维地质模型效果图；

图3是车西地区应用该运移模拟方法形成的Ng沉积末时期运移后各地层叠加的油气分布强度效果图。

流程结束。