[发明专利]一种网格划分方法及系统

说明书

技术领域

本发明涉及数值试井分析技术领域，更具体地，特别是涉及一种网格划分方法及系统。

背景技术

油气井的井底瞬态压力数值模拟技术的应用之一是进行井底压力的试井分析，而试井分析是认识油藏、进行油藏评价和生产动态监测的重要分析手段。试井分析技术按渗流方程求解的方式，可分为解析试井分析方法与数值试井分析方法。

数值试井方法中最重要的步骤为网格划分，所划分的网络的数量与分布直接影响油藏数值模拟的计算效率、计算精度和计算时间；且网格对复杂求解区域的适应能力和逼近程度以及网格生产的自动化和智能化程度等，直接关系到数值模拟和数值结果的好坏。

现有的PEBI网格划分方法当井在边界或断层附近时，网格划分往往不成功，需要多次参数调节方能画出正确的网格，而且往往导致网格数量大，进而导致计算效率低，计算时间长。

发明内容

有鉴于此，本发明实施例提供一种网格划分方法及系统，以降低网格数量，进一步提高计算效率，缩短计算时间，降低用户负担。技术方案如下：

基于本发明的一方面，公开了一种网格划分方法，应用于试井分析领域，包括：创建油藏三维模型的水平二维面，其特征在于，还包括：

对所述水平二维面进行点模块和线模块划分，并在所述点模块和所述线模块中布点；

将所述水平二维面划分为至少一个独立区域，对所述独立区域按照径向规律布点，所述径向规律为：依据到独立区域中井点的距离值，增加网格尺度，相邻点的网格尺度增加值逐渐增大，且到独立区域中井点的距离值相同的点的网格尺度相同，所述独立区域的个数与水平二维面中井的数目相同；

去除所述独立区域、点模块和线模块中相互重合的点；

对所述独立区域中剩余的点与点模块和线模块中的点的进行Delaunay三角剖分和Voronoi网格剖分，得出网格；

删除所述网格中的无效网格，将删除无效网格后的所有网格作为二维可变垂直平分PEBI非结构网格；

对二维可变PEBI非结构网格中的每个网格进行垂向拓展，拓展后的所有网格作为三维可变PEBI非结构网格。

优选地，对所述水平二维面进行点模块和线模块划分，并在所述点模块和所述线模块中布点具体为：

将所述水平二维面的边界线的顶点，所述水平二维面的断层线的起点、转折点和终点，断层线与断层线的交点，断层线与边界线的交点作为点模块，所述点模块的大小预先设定，且所述点模块至少包括一条关联边；

将所述水平二维面中的边界线和断层线作为线模块，所述线模块的大小预先设定；

当所述点模块包括至少两条关联边，在所述关联边两侧对称位置各布置一个点，当所述点模块包括一个关联边，在所述关联边两侧对称位置各布置一个点，同时以所述点模块中心对称布点；

获取线模块各处的网格尺度，所述线模块按照轴对称布点，两对称点之间的距离等于该处的网格尺度。

优选地，定义水平二维面中点到直井的距离为点到直井的井点的距离，点到水平井的距离为点到水平井的两个端点的距离的较小值，点到裂缝井的距离为点到裂缝井的两个端点的距离的较小值；

所述将所述水平二维面划分为至少一个独立区域包括：

选取水平二维面中直井的井点，水平井和裂缝井的两个端点；

获取水平二维面中各个点分别到直井、水平井和裂缝井的距离；

选取水平二维面中的任意一个直井，将水平二维面中到选取的直井的距离小于到其他直井、水平井和裂缝井的距离的点集所覆盖区域划分为选取的直井的独立区域；

选取水平二维面中的任意一个水平井，将水平二维面中到选取的水平井的距离小于到其他水平井、直井和裂缝井的距离的点集所覆盖区域划分为选取的水平井的独立区域；

选取水平二维面中的任意一个裂缝井，将水平二维面中到选取的裂缝井的距离小于到其他裂缝井、直井和水平井的距离的点集所覆盖区域划分为选取的裂缝井的独立区域。

优选地，将所述水平二维面划分为至少一个独立区域包括：

选取水平二维面中直井的井点，水平井和裂缝井的两个端点，将井点和端点作为划分点；

对所述划分点进行三角剖分，获取三角网格；

对所述三角网格的各个边做垂直平分线，将与井点相邻的垂直平分线相交形成的区域划分为该井点的独立区域。

优选地，对所述独立区域按照径向规律布点包括：