专利名称：一种网格划分方法及系统的制作方法
技术领域：
本发明涉及数值试井分析技术领域，更具体地，特别是涉及一种网格划分方法及系统。
背景技术：
油气井的井底瞬态压力数值模拟技术的应用之一是进行井底压力的试井分析，而试井分析是认识油藏、进行油藏评价和生产动态监测的重要分析手段。试井分析技术按渗流方程求解的方式，可分为解析试井分析方法与数值试井分析方法数值试井方法中最重要的步骤为网格划分，所划分的网络的数量与分布直接影响油藏数值模拟的计算效率、计算精度和计算时间；且网格对复杂求解区域的适应能力和逼近程度以及网格生产的自动化和智能化程度等，直接关系到数值模拟和数值结果的好坏。现有的PEBI网格划分方法当井在边界或断层附近时，网格划分往往不成功，需要多次参数调节方能画出正确的网格，而且往往导致网格数量大，进而导致计算效率低，计算时间长。
发明内容
有鉴于此，本发明实施例提供一种网格划分方法及系统，以降低网格数量，进一步提高计算效率，缩短计算时间，降低用户负担。技术方案如下基于本发明的一方面，公开了一种网格划分方法，应用于试井分析领域，包括创建油藏三维模型的水平二维面，其特征在于，还包括对所述水平二维面进行点模块和线模块划分，并在所述点模块和所述线模块中布占.
^ \\\ 将所述水平二维面划分为至少一个独立区域，对所述独立区域按照径向规律布点，所述径向规律为依据到独立区域中井点的距离值，增加网格尺度，相邻点的网格尺度增加值逐渐增大，且到独立区域中井点的距离值相同的点的网格尺度相同，所述独立区域的个数与水平二维面中井的数目相同；去除所述独立区域、点模块和线模块中相互重合的点；对所述独立区域中剩余的点与点模块和线模块中的点的进行Delaimay三角剖分和Voronoi网格剖分，得出网格；删除所述网格中的无效网格，将删除无效网格后的所有网格作为二维可变垂直平分PEBI非结构网格；对二维可变PEBI非结构网格中的每个网格进行垂向拓展，拓展后的所有网格作为三维可变PEBI非结构网格。优选地，对所述水平二维面进行点模块和线模块划分，并在所述点模块和所述线模块中布点具体为将所述水平二维面的边界线的顶点，所述水平二维面的断层线的起点、转折点和终点，断层线与断层线的交点，断层线与边界线的交点作为点模块，所述点模块的大小预先设定，且所述点模块至少包括一条关联边；将所述水平二维面中的边界线和断层线作为线模块，所述线模块的大小预先设定；当所述点模块包括至少两条关联边，在所述关联边两侧对称位置各布置一个点， 当所述点模块包括一个关联边，在所述关联边两侧对称位置各布置一个点，同时以所述点模块中心对称布点；获取线模块各处的网格尺度，所述线模块按照轴对称布点，两对称点之间的距离等于该处的网格尺度。优选地，定义水平二维面中点到直井的距离为点到直井的井点的距离，点到水平井的距离为点到水平井的两个端点的距离的较小值，点到裂缝井的距离为点到裂缝井的两个端点的距离的较小值；所述将所述水平二维面划分为至少一个独立区域包括选取水平二维面中直井的井点，水平井和裂缝井的两个端点；获取水平二维面中各个点分别到直井、水平井和裂缝井的距离；选取水平二维面中的任意一个直井，将水平二维面中到选取的直井的距离小于到其他直井、水平井和裂缝井的距离的点集所覆盖区域划分为选取的直井的独立区域；选取水平二维面中的任意一个水平井，将水平二维面中到选取的水平井的距离小于到其他水平井、直井和裂缝井的距离的点集所覆盖区域划分为选取的水平井的独立区域；选取水平二维面中的任意一个裂缝井，将水平二维面中到选取的裂缝井的距离小于到其他裂缝井、直井和水平井的距离的点集所覆盖区域划分为选取的裂缝井的独立区域。优选地，将所述水平二维面划分为至少一个独立区域包括选取水平二维面中直井的井点，水平井和裂缝井的两个端点，将井点和端点作为划分点；对所述划分点进行三角剖分，获取三角网格；对所述三角网格的各个边做垂直平分线，将与井点相邻的垂直平分线相交形成的区域划分为该井点的独立区域。优选地，对所述独立区域按照径向规律布点包括选取独立区域中的直井的独立区域，以直井的独立区域中直井的井点为圆心进行圆模块网格划分，对所述圆模块按照径向规律布点，所述圆模块包含在该直井的独立区域， 且所述圆模块的大小预先设定；获取独立区域中的水平井的独立区域，将水平井的独立区域中水平井的两个端点作为椭圆的两个焦点进行椭圆模块网格划分，对所述椭圆模块两端区域按径向规律布点， 椭圆模块的其他区域按规则网格布点，所述椭圆模块包含在水平井独立区域内，且所述椭圆模块的大小预先设定；选择独立区域中的裂缝井的独立区域，以裂缝井的独立区域中的两个端点作为椭圆的两个焦点进行椭圆模块网格划分，对所述椭圆模块两端区域按径向规律布点，椭圆模块的其他区域按规则网格布点，所述椭圆模块包含在裂缝井独立区域内，且所述椭圆模块的大小预先设定；独立区域中的其他区域，参考径向规律进行布点；当所述圆模块和椭圆模块相交时，缩小相交模块的大小直至模块互不相交。基于本发明的另一方面，公开了一种网格划分系统，应用于试井分析领域，包括 面创建模块，用于创建油藏三维模型的水平二维面，其特征在于，还包括划分布点模块，用于对所述水平二维面进行点模块和线模块划分，并在所述点模块和所述线模块中布点；独立区域划分布点模块，用于将所述水平二维面划分为至少一个独立区域，对所述独立区域按照径向规律布点，所述径向规律为依据到独立区域中井点的距离值，增加网格尺度，相邻点的网格尺度增加值逐渐增大，且到独立区域中井点的距离值相同的点的网格尺度相同，所述独立区域的个数与水平二维面中井的数目相同；点去除模块，用于去除所述独立区域、点模块、线模块、圆模块和椭圆模块中相互重合的点；剖分模块，用于对所述独立区域中剩余的点与点模块、线模块、圆模块和椭圆模块中的点的进行Delaunay三角剖分和Voronoi网格剖分，得出网格；二维网格获取模块，用于删除所述网格中的无效网格，将删除无效网格后的所有网格作为二维可变垂直平分PEBI非结构网格；三维网格获取模块，用于对二维可变PEBI非结构网格中格的每个网格进行垂直拓展，拓展后的所有网格作为三维可变PEBI非结构网格。优选地，所述划分布点模块包括点模块划分单元，用于将所述水平二维面的边界线的顶点，所述水平二维面的断层线的起点、转折点和终点，断层线与断层线的交点，断层线与边界线的交点作为点模块， 所述点模块的大小预先设定，且所述点模块至少包括一条关联边；线模块划分单元，用于将所述水平二维面中的边界线和断层线作为线模块，所述线模块的大小预先设定；点模块布点单元，用于当所述点模块包括至少两条关联边，在所述关联边两侧对称位置各布置一个点，当所述点模块包括一个关联边，在所述关联边两侧对称位置各布置一个点，同时以所述点模块中心对称布点；线模块布点单元，用于获取线模块各处的网格尺度，所述线模块按照轴对称布点， 两对称点之间的距离等于该处的网格尺度。优选地，定义水平二维面中点到直井的距离为点到直井的井点的距离，点到水平井的距离为点到水平井的两个端点的距离中数值小的距离，点到裂缝井的距离为点到裂缝井的两个端点的距离中数值小的距离；所述独立区域划分布点模块包括点选取单元，用于选取水平二维面中直井的井点，水平井和裂缝井的两个端点；距离获取单元，用于获取水平二维面中各个点分别到直井、水平井和裂缝井的距罔；直井独立区域划分单元，用于选取水平二维面中的任意一个直井，将水平二维面中到选取的直井的距离小于到其他直井、水平井和裂缝井的距离的点集所覆盖区域划分为选取的直井的独立区域；水平井独立区域划分单元，用于选取水平二维面中的任意一个水平井，将水平二维面中到选取的水平井的距离小于到其他水平井、直井和裂缝井的距离的点集所覆盖区域划分为选取的水平井的独立区域；裂缝井独立区域划分单元，用于选取水平二维面中的任意一个裂缝井，将水平二维面中到选取的裂缝井的距离小于到其他裂缝井、直井和水平井的距离的点集所覆盖区域划分为选取的裂缝井的独立区域。
优选地，所述独立区域划分布点模块包括划分点选取单元，用于选取水平二维面中直井的井点，水平井和裂缝井的两个端点，将井点和端点作为划分点；三角网格获取单元，用于对所述划分点进行三角剖分，获取三角网格；独立区域划分单元，用于对所述三角网格的各个边做垂直平分线，将与井点相邻的垂直平分线相交形成的区域划分为该井点的独立区域。优选地，所述独立区域划分布点模块还包括圆模块布点单元，用于选取独立区域中的直井的独立区域，以直井的独立区域中直井的井点为圆心进行圆模块划分，对所述圆模块按照径向规律布点，所述圆模块包含在直井独立区域，且所述圆模块的大小预先设定；第一椭圆模块布点单元，用于获取独立区域中的水平井的独立区域，将水平井的独立区域中水平井的两个端点作为椭圆的两个交点进行椭圆模块，对所述椭圆模块两端区域按径向规律布点，其他区域按规则网格布点，所述椭圆模块包含在水平井独立区域内，且所述椭圆模块的大小预先设定；第二椭圆模块布点单元，用于选择独立区域中的裂缝井的独立区域，以裂缝井的独立区域中裂缝井的两个端点为椭圆的两个焦点进行椭圆模块划分，对所述椭圆模块两端区域按径向规律布点，其他区域按规则网格布点，所述椭圆模块包含在裂缝井独立区域内， 且所述椭圆模块的大小预先设定；其他区域布点单元，用于在独立区域中的其他区域，参考径向规律进行布点。应用上述技术方案，对独立区域布点时，按照径向规律布点，所述径向规律为依据至IJ独立区域中井点的距离值，增加网格尺度，相邻点的网格尺度增加值逐渐增大，且到独立区域中井点的距离值相同的点的网格尺度相同。采用上述布点方式，降低了独立区域点的个数，进而在对点进行网络划分时，进一步降低网络数量，提高计算效率，缩短计算时间，降低用户负担。发明人经过多次实验发现，当油藏面积为1000m*1000m时，采用均勻PEBI划分时， 划分的网格数量为4598个，计算时间为31. 328秒，而采用本发明提供的方法，划分的网格数量为584个，计算时间为4. 531秒，相对于现有的均勻PEBI划分，计算效率明显提高。由于本发明实施例提供的网格划分方法计算时间的缩短，计算效率提高，所以采用上述网格划分方法进一步降低用户负担。
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明中记载的一些实施例，对于本领域普通技术人员来讲，还可以根据这些附图获得其他的附图。图1为本发明实施例公开的网格划分方法的一种流程图；图2为油藏三维模型的水平二维面的示意图；图3为图2所示的水平二维面经过点模块和线模块划分后的示意图；图4为点模块的一种布点示意图；图5为点模块的另一种布点示意图；图6为线模块的布点示意图；图7为本发明实施例公开的网格划分方法中独立区域划分的一种流程图；图8为图3所示的水平二维面经过独立区域划分后的示意图；图9为本发明实施例公开的网格划分方法中独立区域划分的另一种流程图；图10为图2所示的水平二维面的独立区域划分示意图；图11为本发明实施例公开的网格划分方法中独立区域布点的流程图；图12为圆模块的布点示意图；图13为椭圆模块的布点示意图；图14为独立区域的布点示意图；图15为图2所示的水平二维面的布点示意图；图16为对图11所示的点进行剖分的示意图；图17为二维可变PEBI非结构网格示意图；图18为三维可变PEBI非结构网格示意图；图19为本发明实施例公开的网格划分系统的结构示意图；图20为本发明实施例公开的网格划分系统中划分布点的结构示意图；图21为本发明实施例公开的网格划分系统中独立区域划分布点模块的一种结构示意图；图22为本发明实施例公开的网格划分系统中独立区域划分布点模块的另一种结构示意图；图23为本发明实施例公开的网格划分系统中独立区域划分布点模块的再一种结构示意图。
具体实施例方式发明人经过研究发现，目前，用于网格划分的是PEBI非结构网格和三角网格在进行网格划分时，所选用的网格尺度为均勻网格尺度，进而导致网格数量增加，进一步导致计算效率低，计算时间长。为了解决上述问题，本发明实施例公开一种网格划分方法，对独立区域布点时，按照径向规律布点。其中径向规律为依据到独立区域中井点的距离值，增加网格尺度，相邻点的网格尺度增加值逐渐增大，且到独立区域中井点的距离值相同的点的网格尺度相同。 即对独立区域布点时，采用逐渐递增的网格尺度进行布点，降低了独立区域点的个数，进一步在对点进行网络划分时，降低网络数量，提高计算效率，缩短计算时间，降低用户负担。为了使本技术领域的人员更好地理解本发明中的技术方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员所获得的所有其他实施例，都应当属于本发明保护的范围请参阅图1，图1为本发明实施例公开的网格划分方法的流程图，包括SlOl 创建油藏三维模型的水平二维面。S102 对水平二维面进行点模块和线模块划分，并在点模块和线模块中布点。其中对水平二维面进行点模块和线模块划分具体为将所述水平二维面的边界线的顶点，所述水平二维面的断层线的起点、转折点和终点，断层线与断层线的交点，断层线与边界线的交点作为点模块，所述点模块的大小预先设定，且所述点模块至少包括一条关联边；将所述水平二维面中的边界线和断层线作为线模块，所述线模块的大小预先设定。下面以一个具体的油藏三维模型的水平二维面具体说明。图2是一个油藏三维模型的水平二维面的示意图，该示意图所示的油藏具有一个水平井，三个垂直井，示意图的边界为油藏的边界，区域中的线段为断层。对图2所示的油藏三维模型的水平二维面的示意图，按照上述点模块和线模块划分方式进行划分，得出的划分示意图如图3所示，其中图中线段为线模块，线模块的个数为14。线段上的点为点模块，点模块的大小，即点模块的直径是预先由用户依据油藏三维模型设定的，各个点模块的直径可以相同，也可以不同，对此不加以阐述。图3中共有16个点模块，与每个点模块相连的线模块为该点模块的关联边。上述在点模块和线模块中布点具体为当点模块包括至少两条关联边，在关联边两侧对称位置各布置一个点，如图4所示，图4为点模块包括五个关联边，每个关联边两侧对称布点。图5为点模块包括一个关联边时的布点示意图。当点模块包括一个关联边，在关联边两侧对称位置各布置一个点，同时以点模块中心对称布点。获取线模块各处的网格尺度，线模块按照轴对称布点，两对称点之间的距离等于该处的网格尺度，如图6所示，图6为线模块布点示意图。线模块各处的网格尺度与该处到井点的距离有关。线模块上某处到井点的距离小，该处的网格尺度小；线模块上某处到井点的距离大，该处的网格尺度大。S103 将所述水平二维面划分为至少一个独立区域，对所述独立区域按照径向规律布点，所述径向规律为依据到独立区域中井点的距离值，增加网格尺度，相邻点的网格尺度增加值逐渐增大，且到独立区域中井点的距离值相同的点的网格尺度相同，所述独立区域的个数与水平二维面中井的数目相同。其中独立区域的划分有两种方式，一种是按照距离划分；一种是按照三角剖分划分。下面对两种方式一一详细说明。首先介绍按照距离划分。定义水平二维面中点到直井的距离为点到直井的井点的距离，点到水平井的距离为点到水平井的两个端点的距离中的较小值，点到裂缝井的距离为点到裂缝井的两个端点的距离中的较小值。其中水平井的两个端点为水平井的井点，裂缝井的两个端点为裂缝井的井点。步骤S103的流程图如图7所示，包括S1031 选取水平二维面中直井的井点，水平井和裂缝井的两个端点。S1032:获取水平二维面中各个点分别到直井、水平井和裂缝井的距离。S1033 选取水平二维面中的任意一个直井，将水平二维面中到选取的直井的距离小于到其他直井、水平井和裂缝井的距离的点集所覆盖区域划分为选取的直井的独立区域。如图3所示的水平二维面，选取最下面的直井，计算水平二维面中点到该直井的距离， 获取到该直井的距离小于到图3中其他两个直井和一个水平井的距离的点，将所有点构成的点集所覆盖的区域作为选取的最下面的直井的独立区域。S1034 选取水平二维面中的任意一个水平井，将水平二维面中到选取的水平井的距离小于到其他水平井、直井和裂缝井的距离的点集所覆盖区域划分为选取的水平井的独立区域。 S1035 选取水平二维面中的任意一个裂缝井，将水平二维面中到选取的裂缝井的距离小于到其他裂缝井、直井和水平井的距离的点集所覆盖区域划分为选取的裂缝井的独立区域。图3经过图7所示的独立区域划分流程图之后，得到如图8所示的水平二维面的划分后的示意图。其中虚线相交形成的区域为独立区域。图8所示的水平二维面中独立区域与线模块相交，对此不需要做处理。上述步骤S103还可以采用图9所示的独立区域划分流程图，包括S1036:选取水平二维面中直井的井点，水平井和裂缝井的两个端点，将井点和端点作为划分点；S1037 对所述划分点进行三角剖分，获取三角网格；S1038:对所述三角网格的各个边做垂直平分线，将与井点相邻的垂直平分线相交形成的区域划分为该井点的独立区域。图10为图2所示的水平二维面经过图9所示的独立区域划分流程图划分后得出的独立区域示意图。其中点为直井井点和水平井端点，实线为对点进行划分后得到的三角网格，虚线相交形成的区域为独立区域。上述步骤S103对独立区域按照径向规律布点流程图请参阅图11，包括S1039 选取独立区域中的直井的独立区域，以直井的独立区域中直井的井点为圆心进行圆模块网格划分，对圆模块按照径向规律布点，圆模块包含在直井独立区域，且圆模块的大小预先设定。对圆模块按照径向规律布点具体为以递增网格尺度作为半径将圆模块划分为多个同心圆，且将多个同心圆按照某个角度值进行均勻划分或者将多个同心圆按照多个角度值进行划分，在划分后的每个子区域内布设一个点。具体示意图请参阅图12，图12为圆模块按照径向规律布点的示意图。S1040:获取独立区域中的水平井的独立区域，将水平井的独立区域中水平井的两个端点作为椭圆的两个焦点进行椭圆模块网格划分，对椭圆模块两端区域按径向规律布点，椭圆模块的其他区域按规则网格布点。椭圆模块包含在水平井独立区域内，且椭圆模块的大小预先设定。对椭圆模块按照径向规律布点具体为以递增网格尺度作为半径将椭圆模块划分为多个同焦点椭圆，且将多个同焦点椭圆的两端区域按照某个角度值进行均勻划分或者将多个同焦点椭圆按照多个角度值进行划分，将多个同焦点椭圆的其他区域进行划分，在划分后的每个子区域内布设一个点。其中同焦点椭圆为具有相同焦点的椭圆。上述对椭圆模块布点的示意图请参阅图13，图13为椭圆模块按照径向规律布点的示意图。S1041 选择独立区域中的裂缝井的独立区域，以裂缝井的独立区域中裂缝井的两个端点为椭圆的两个焦点进行椭圆模块网格划分，对椭圆模块两端区域按径向规律布点， 椭圆模块的其他区域按规则网格布点。椭圆模块包含在裂缝井独立区域内，且椭圆模块的大小预先设定。对裂缝井独立区域布点请参阅步骤S1039和图13的具体阐述。S1042 独立区域中的其他区域，参考径向规律进行布点。其中其他区域指除圆模块和椭圆模块覆盖区域之外的其他区域上述图2所示的水平二维面的独立区域的布点示意图请参阅图14，从图14中可以看出，独立区域依据径向规律布点，能够更加直观的表明井点分部情况，能够更加直观的分辨井点区域、边界区域和断层区域。需要说明的是独立区域中的圆模块和椭圆模块，以及椭圆模块之间是不能相交的，当上述模块相交，则需要缩小相交模块至模块互不相交。S104 去除所述独立区域、点模块和线模块中相互重合的点。请参阅图15，图15为图2所示的水平二维面的布点示意图。S105 对所述独立区域中剩余的点与点模块和线模块中的点的进行Delaimay三角剖分和Voronoi网格剖分，得出网格。其中=DeIaunay三角剖分和Voronoi网格剖分使用现有的Delaimay三角剖分和Voronoi网格剖分，对此不加以阐述。步骤S105对图15所示的布点示意图中的点进行Delaimay三角剖分和Voronoi网格剖分后获取的网格示意图请参阅图16所示。S106:删除所述网格中的无效网格，将删除无效网格后的所有网格作为二维可变 PEBI非结构网格。从图16中可以看出网格中还有无效网格，如边界外的网格，因此，需要将无效网格删除，删除后的网格为二维可变PEBI非结构网格，如图17所示。S107 对二维可变PEBI非结构网格中格的每个网格进行垂向拓展，拓展后的所有网格作为三维可变PEBI非结构网格，如图18所示。应用上述技术方案，对独立区域布点时，按照径向规律布点，所述径向规律为依据到独立区域中井点的距离值，增加网格尺度，相邻点的网格尺度增加值逐渐增大，且到独立区域中井点的距离值相同的点的网格尺度相同。采用上述布点方式，降低了独立区域点的个数，进而在对点进行网络划分时，进一步降低网络数量，提高计算效率，缩短计算时间，降低用户负担。发明人经过多次实验发现，当油藏面积为1000m*1000m时，采用均勻PEBI划分时， 划分的网格数量为4598个，计算时间为31. 328秒，而采用本发明提供的方法，划分的网格数量为584个，计算时间为4. 531秒，相对于现有的均勻PEBI划分，计算效率明显提高。由于本发明实施例提供的网格划分方法计算时间的缩短，计算效率提高，所以采用上述网格划分方法进一步降低用户负担。需要说明的是当油藏面积小于300m*300m，采用均勻PEBI网格划分，而油藏面积大于300m*300m，采用本发明实施例提供的方法。并且发明人发现，油藏面积越大，计算时间相对于采用均勻PEBI网格划分花费的计算时间越少。本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分可以互相参见，每个实施例重点说明的都是与其他实施例的不同之处。通过以上的方法实施例的描述，所属领域的技术人员可以清楚地了解到本发明可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括只读存储器(ROM)、随机存取存储器(RAM)、磁碟或者光盘等各种可以存储程序代码的介质。相应于上面的方法实施例，本发明实施例还公开一种网格划分系统，应用于试井分析领域，结构示意图如图19所示，包括面创建模块11、划分布点模块12、独立区域划分布点模块13、点去除模块14、剖分模块15、二维网格获取模块16和三维网格获取模块17。 其中面创建模块11，用于创建油藏三维模型的水平二维面。划分布点模块12，用于对所述水平二维面进行点模块和线模块划分，并在所述点模块和所述线模块中布点。划分布点模块12的结构示意图请参阅图20，包括点模块划分单元121和线模块划分单元122。点模块划分单元121，用于将所述水平二维面的边界线的顶点，所述水平二维面的断层线的起点、转折点和终点，断层线与断层线的交点，断层线与边界线的交点作为点模块，所述点模块的大小预先设定，且所述点模块至少包括一条关联边。线模块划分单元122，用于将所述水平二维面中的边界线和断层线作为线模块，所述线模块的大小预先设定。上述点模块划分单元121和线模块划分单元122对水平二维面的划分请参阅图3， 图3为水平二维面经过点模块和线模块划分后的示意图。图3中共有16个点模块，与每个点模块相连的线模块为该点模块的关联边。线模块的个数为14。上述图20所示的划分布点模块12还包括点模块布点单元123，用于当所述点模块包括至少两条关联边，在所述关联边两侧对称位置各布置一个点，当所述点模块包括一个关联边，在所述关联边两侧对称位置各布置一个点，同时以所述点模块中心对称布点，如图4所示。线模块布点单元124，用于获取线模块各处的网格尺度，所述线模块按照轴对称布点，两对称点之间的距离等于该处的网格尺度，如图6所示。图6为线模块布点示意图。线模块各处的网格尺度与该处到井点的距离有关。线模块上某处到井点的距离小，该处的网格尺度小；线模块上某处到井点的距离大，该处的网格尺度大。独立区域划分布点模块13，用于将所述水平二维面划分为至少一个独立区域，对所述独立区域按照径向规律布点，所述径向规律为依据到独立区域中井点的距离值，增加网格尺度，相邻点的网格尺度增加值逐渐增大，且到独立区域中井点的距离值相同的点的网格尺度相同，所述独立区域的个数与水平二维面中井的数目相同。
独立区域划分布点模块13对独立区域的划分可以按照距离划分或者按照三角剖分划分。按照距离划分时，定义水平二维面中点到直井的距离为点到直井的井点的距离，点到水平井的距离为点到水平井的两个端点的距离的较小值，点到裂缝井的距离为点到裂缝井的两个端点的距离的较小值。其中水平井的两个端点为水平井的井点，裂缝井的两个端点为裂缝井的井点。独立区域划分布点模块13的结构示意图请参阅图21，包括点选取单元131、距离获取单元132、直井独立区域划分单元133、水平井独立区域划分单元134和裂缝井独立区域划分单元135。其中点选取单元131，用于选取水平二维面中直井的井点，水平井和裂缝井的两个端
点ο距离获取单元132，用于获取水平二维面中各个点分别到直井、水平井和裂缝井的距离。直井独立区域划分单元133，用于选取水平二维面中的任意一个直井，将水平二维面中到选取的直井的距离小于到其他直井、水平井和裂缝井的距离的点集所覆盖区域划分为选取的直井的独立区域。如图3所示的水平二维面，选取最下面的直井，计算水平二维面中点到该直井的距离，获取到该直井的距离小于到图3中其他两个直井和一个水平井的距离的点，将所有点构成的点集所覆盖的区域作为选取的最下面的直井的独立区域。水平井独立区域划分单元134，用于选取水平二维面中的任意一个水平井，将水平二维面中到选取的水平井的距离小于到其他水平井、直井和裂缝井的距离的点集所覆盖区域划分为选取的水平井的独立区域。裂缝井独立区域划分单元135，用于选取水平二维面中的任意一个裂缝井，将水平二维面中到选取的裂缝井的距离小于到其他裂缝井、直井和水平井的距离的点集所覆盖区域划分为选取的裂缝井的独立区域。上述独立区域划分布点模块13还可以按照三角剖分划分，其结构示意图如图22 所示，包括划分点选取单元136、三角网格获取单元137和独立区域划分单元138。其中划分点选取单元136，用于选取水平二维面中直井的井点，水平井和裂缝井的两个端点，将井点和端点作为划分点。三角网格获取单元137，用于对所述划分点进行三角剖分，获取三角网格。独立区域划分单元138，用于对所述三角网格的各个边做垂直平分线，将与井点相邻的垂直平分线相交形成的区域划分为该井点的独立区域。独立区域划分布点模块13对油藏三维模型的水平二维面划分后的示意图请参阅图8或图10。独立区域划分布点模块13在划分水平二维面后，对各个独立区域进行布点， 独立区域划分布点模块13在图21或22的基础上还包括圆模块布点单元139、第一椭圆模块布点单元140、第二椭圆模块布点单元141和其他区域布点单元142，如图23所示。图 23是在图21所示的独立区域划分布点模块13基础上增加其他单元。其中圆模块布点单元139，用于选取独立区域中的直井的独立区域，以直井的独立区域中直井的井点为圆心进行圆模块网格划分，对所述圆模块按照径向规律布点，所述圆模块包含在直井独立区域，且所述圆模块的大小预先设定。对圆模块按照径向规律布点具体为以递增网格尺度作为半径将圆模块划分为多个同心圆，且将多个同心圆按照某个角度值进行均勻划分或者将多个同心圆按照多个角度值进行划分，在划分后的每个子区域内布设一个点。具体示意图请参阅图12，图12为圆模块按照径向规律布点的示意图。第一椭圆模块布点单元140，用于获取独立区域中的水平井的独立区域，将水平井的独立区域中水平井的两个端点作为椭圆的两个交点进行椭圆模块网格划分，对所述椭圆模块两端区域按径向规律布点，其他区域按规则网格布点，所述椭圆模块包含在水平井独立区域内，且所述椭圆模块的大小预先设定。对椭圆模块按照径向规律布点具体为以递增网格尺度作为半径将椭圆模块划分为多个同焦点椭圆，且将多个同焦点椭圆的两端区域按照某个角度值进行均勻划分或者将多个同焦点椭圆按照多个角度值进行划分，将多个同焦点椭圆的其他区域进行划分，在划分后的每个子区域内布设一个点。其中同焦点椭圆为具有相同焦点的椭圆。上述对椭圆模块布点的示意图请参阅图13，图13为椭圆模块按照径向规律布点的示意图。第二椭圆模块布点单元141，用于选择独立区域中的裂缝井的独立区域，以裂缝井的独立区域中裂缝井的两个端点为椭圆的两个焦点进行椭圆模块网格划分，对所述椭圆模块两端区域按径向规律布点，其他区域按规则网格布点，所述椭圆模块包含在裂缝井独立区域内，且所述椭圆模块的大小预先设定。其他区域布点单元142，用于在独立区域中的其他区域，参考径向规律进行布点。 其中其他区域指除圆模块和椭圆模块覆盖区域之外的其他区域。上述图2所示的水平二维面的独立区域的布点示意图请参阅图14，从图14中可以看出，独立区域依据径向规律布点，能够更加直观的表明井点分部情况，能够更加直观的分辨井点区域、边界区域和断层区域。需要说明的是独立区域中的圆模块和椭圆模块，以及椭圆模块之间是不能相交的，当上述模块相交，则需要缩小相交模块至模块互不相交。点去除模块14，用于去除所述独立区域、点模块、线模块、圆模块和椭圆模块中相互重合的点。请参阅图15，图15为图2所示的水平二维面的布点示意图。剖分模块15，用于对所述独立区域中剩余的点与点模块、线模块、圆模块和椭圆模块中的点的进行Delaunay三角剖分和Voronoi网格剖分，得出网格。其中=DeIaunay三角剖分和Vorono i网格剖分使用现有的De 1 aunay三角剖分和Vorono i网格剖分，对此不加以阐述。网格示意图请参阅图16所示。二维网格获取模块16，用于删除所述网格中的无效网格，将删除无效网格后的所有网格作为二维可变垂直平分PEBI非结构网格。从图16中可以看出网格中还有无效网格，如边界外的网格，因此，需要将无效网格删除，删除后的网格为二维可变PEBI非结构网格，如图17所示。三维网格获取模块17，用于对二维可变PEBI非结构网格中格的每个网格进行垂直拓展，拓展后的所有网格作为三维可变PEBI非结构网格，如图18所示。本发明实施例提供的网格划分系统可以参照上述方法实施例完成对网格的划分， 因此，应用该系统，降低了独立区域点的个数，进而在对点进行网络划分时，进一步降低网络数量，提高计算效率，缩短计算时间，降低用户负担。对于设备实施例而言，由于其基本相应于方法实施例，所以相关之处参见方法实施例的部分说明即可。以上所描述的设备实施例仅仅是示意性的，其中所述作为分离部件
16说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性劳动的情况下，即可以理解并实施。在本申请所提供的几个实施例中，应该理解到，所揭露的设备和方法，在没有超过本申请的精神和范围内，可以通过其他的方式实现。当前的实施例只是一种示范性的例子， 不应该作为限制，所给出的具体内容不应该限制本申请的目的。例如，所述单元或子单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或多个子单元结合一起。另外，多个单元可以或组件可以结合或者可以集成到另一个设备，或一些特征可以忽略，或不执行。以上所述仅是本发明的具体实施方式
，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
权利要求
1.一种网格划分方法，应用于试井分析领域，包括创建油藏三维模型的水平二维面， 其特征在于，还包括对所述水平二维面进行点模块和线模块划分，并在所述点模块和所述线模块中布点； 将所述水平二维面划分为至少一个独立区域，对所述独立区域按照径向规律布点，所述径向规律为依据到独立区域中井点的距离值，增加网格尺度，相邻点的网格尺度增加值逐渐增大，且到独立区域中井点的距离值相同的点的网格尺度相同，所述独立区域的个数与水平二维面中井的数目相同；去除所述独立区域、点模块和线模块中相互重合的点；对所述独立区域中剩余的点与点模块和线模块中的点的进行Delaunay三角剖分和 Voronoi网格剖分，得出网格；删除所述网格中的无效网格，将删除无效网格后的所有网格作为二维可变垂直平分 PEBI非结构网格；对二维可变PEBI非结构网格中的每个网格进行垂向拓展，拓展后的所有网格作为三维可变PEBI非结构网格。
2.根据权利要求1所述的网格划分方法，其特征在于，对所述水平二维面进行点模块和线模块划分，并在所述点模块和所述线模块中布点具体为将所述水平二维面的边界线的顶点，所述水平二维面的断层线的起点、转折点和终点， 断层线与断层线的交点，断层线与边界线的交点作为点模块，所述点模块的大小预先设定， 且所述点模块至少包括一条关联边；将所述水平二维面中的边界线和断层线作为线模块，所述线模块的大小预先设定； 当所述点模块包括至少两条关联边，在所述关联边两侧对称位置各布置一个点，当所述点模块包括一个关联边，在所述关联边两侧对称位置各布置一个点，同时以所述点模块中心对称布点；获取线模块各处的网格尺度，所述线模块按照轴对称布点，两对称点之间的距离等于该处的网格尺度。
3.根据权利要求1所述的网格划分方法，其特征在于，定义水平二维面中点到直井的距离为点到直井的井点的距离，点到水平井的距离为点到水平井的两个端点的距离的较小值，点到裂缝井的距离为点到裂缝井的两个端点的距离的较小值；所述将所述水平二维面划分为至少一个独立区域包括 选取水平二维面中直井的井点，水平井和裂缝井的两个端点； 获取水平二维面中各个点分别到直井、水平井和裂缝井的距离； 选取水平二维面中的任意一个直井，将水平二维面中到选取的直井的距离小于到其他直井、水平井和裂缝井的距离的点集所覆盖区域划分为选取的直井的独立区域；选取水平二维面中的任意一个水平井，将水平二维面中到选取的水平井的距离小于到其他水平井、直井和裂缝井的距离的点集所覆盖区域划分为选取的水平井的独立区域；选取水平二维面中的任意一个裂缝井，将水平二维面中到选取的裂缝井的距离小于到其他裂缝井、直井和水平井的距离的点集所覆盖区域划分为选取的裂缝井的独立区域。
4.根据权利要求1所述的网格划分方法，其特征在于，将所述水平二维面划分为至少一个独立区域包括选取水平二维面中直井的井点，水平井和裂缝井的两个端点，将井点和端点作为划分点对所述划分点进行三角剖分，获取三角网格；对所述三角网格的各个边做垂直平分线，将与井点相邻的垂直平分线相交形成的区域划分为该井点的独立区域。
5.根据权利要求1至4任意一项所述的网格划分方法，其特征在于，对所述独立区域按照径向规律布点包括选取独立区域中的直井的独立区域，以直井的独立区域中直井的井点为圆心进行圆模块网格划分，对所述圆模块按照径向规律布点，所述圆模块包含在该直井的独立区域，且所述圆模块的大小预先设定；获取独立区域中的水平井的独立区域，将水平井的独立区域中水平井的两个端点作为椭圆的两个焦点进行椭圆模块网格划分，对所述椭圆模块两端区域按径向规律布点，椭圆模块的其他区域按规则网格布点，所述椭圆模块包含在水平井独立区域内，且所述椭圆模块的大小预先设定；选择独立区域中的裂缝井的独立区域，以裂缝井的独立区域中的两个端点作为椭圆的两个焦点进行椭圆模块网格划分，对所述椭圆模块两端区域按径向规律布点，椭圆模块的其他区域按规则网格布点，所述椭圆模块包含在裂缝井独立区域内，且所述椭圆模块的大小预先设定；独立区域中的其他区域，参考径向规律进行布点；当所述圆模块和椭圆模块相交时，缩小相交模块的大小直至模块互不相交。
6.一种网格划分系统，应用于试井分析领域，包括面创建模块，用于创建油藏三维模型的水平二维面，其特征在于，还包括划分布点模块，用于对所述水平二维面进行点模块和线模块划分，并在所述点模块和所述线模块中布点；独立区域划分布点模块，用于将所述水平二维面划分为至少一个独立区域，对所述独立区域按照径向规律布点，所述径向规律为依据到独立区域中井点的距离值，增加网格尺度，相邻点的网格尺度增加值逐渐增大，且到独立区域中井点的距离值相同的点的网格尺度相同，所述独立区域的个数与水平二维面中井的数目相同；点去除模块，用于去除所述独立区域、点模块、线模块、圆模块和椭圆模块中相互重合的点；剖分模块，用于对所述独立区域中剩余的点与点模块、线模块、圆模块和椭圆模块中的点的进行Delaunay三角剖分和Voronoi网格剖分，得出网格；二维网格获取模块，用于删除所述网格中的无效网格，将删除无效网格后的所有网格作为二维可变垂直平分PEBI非结构网格；三维网格获取模块，用于对二维可变PEBI非结构网格中格的每个网格进行垂直拓展， 拓展后的所有网格作为三维可变PEBI非结构网格。
7.根据权利要求6所述的网格划分系统，其特征在于，所述划分布点模块包括点模块划分单元，用于将所述水平二维面的边界线的顶点，所述水平二维面的断层线的起点、转折点和终点，断层线与断层线的交点，断层线与边界线的交点作为点模块，所述点模块的大小预先设定，且所述点模块至少包括一条关联边；线模块划分单元，用于将所述水平二维面中的边界线和断层线作为线模块，所述线模块的大小预先设定；点模块布点单元，用于当所述点模块包括至少两条关联边，在所述关联边两侧对称位置各布置一个点，当所述点模块包括一个关联边，在所述关联边两侧对称位置各布置一个点，同时以所述点模块中心对称布点；线模块布点单元，用于获取线模块各处的网格尺度，所述线模块按照轴对称布点，两对称点之间的距离等于该处的网格尺度。
8.根据权利要求6所述的网格划分系统，其特征在于，定义水平二维面中点到直井的距离为点到直井的井点的距离，点到水平井的距离为点到水平井的两个端点的距离中数值小的距离，点到裂缝井的距离为点到裂缝井的两个端点的距离中数值小的距离；所述独立区域划分布点模块包括点选取单元，用于选取水平二维面中直井的井点，水平井和裂缝井的两个端点； 距离获取单元，用于获取水平二维面中各个点分别到直井、水平井和裂缝井的距离； 直井独立区域划分单元，用于选取水平二维面中的任意一个直井，将水平二维面中到选取的直井的距离小于到其他直井、水平井和裂缝井的距离的点集所覆盖区域划分为选取的直井的独立区域；水平井独立区域划分单元，用于选取水平二维面中的任意一个水平井，将水平二维面中到选取的水平井的距离小于到其他水平井、直井和裂缝井的距离的点集所覆盖区域划分为选取的水平井的独立区域；裂缝井独立区域划分单元，用于选取水平二维面中的任意一个裂缝井，将水平二维面中到选取的裂缝井的距离小于到其他裂缝井、直井和水平井的距离的点集所覆盖区域划分为选取的裂缝井的独立区域。
9.根据权利要求6所述的网格划分系统，其特征在于，所述独立区域划分布点模块包括划分点选取单元，用于选取水平二维面中直井的井点，水平井和裂缝井的两个端点，将井点和端点作为划分点；三角网格获取单元，用于对所述划分点进行三角剖分，获取三角网格； 独立区域划分单元，用于对所述三角网格的各个边做垂直平分线，将与井点相邻的垂直平分线相交形成的区域划分为该井点的独立区域。
10.根据权利要求8或9所述的网格划分系统，其特征在于，所述独立区域划分布点模块还包括圆模块布点单元，用于选取独立区域中的直井的独立区域，以直井的独立区域中直井的井点为圆心进行圆模块划分，对所述圆模块按照径向规律布点，所述圆模块包含在直井独立区域，且所述圆模块的大小预先设定；第一椭圆模块布点单元，用于获取独立区域中的水平井的独立区域，将水平井的独立区域中水平井的两个端点作为椭圆的两个交点进行椭圆模块，对所述椭圆模块两端区域按径向规律布点，其他区域按规则网格布点，所述椭圆模块包含在水平井独立区域内，且所述椭圆模块的大小预先设定；第二椭圆模块布点单元，用于选择独立区域中的裂缝井的独立区域，以裂缝井的独立区域中裂缝井的两个端点为椭圆的两个焦点进行椭圆模块划分，对所述椭圆模块两端区域按径向规律布点，其他区域按规则网格布点，所述椭圆模块包含在裂缝井独立区域内，且所述椭圆模块的大小预先设定； 其他区域布点单元，用于在独立区域中的其他区域，参考径向规律进行布点。
全文摘要
一种网格划分方法，包括进行点模块和线模块划分，并在点模块和所述线模块中布点；将水平二维面划分为至少一个独立区域，对独立区域按照径向规律布点；去除相互重合的点；对点划分得出网格；将删除无效网格后的所有网格为二维可变PEBI非结构网格；对每个网格垂向拓展，拓展后的网格为三维可变PEBI非结构网格。应用本方案，对独立区域布点时，按照径向规律布点，径向规律为依据到独立区域中井点的距离值，增加网格尺度，相邻点的网格尺度增加值逐渐增大，且到独立区域中井点的距离值相同的点的网格尺度相同。采用上述布点方式，降低独立区域点的个数，在对点进行网络划分时，降低网络数量，提高计算效率，缩短计算时间，降低用户负担。
文档编号G06F17/50GK102243680SQ20111020479
公开日2011年11月16日 申请日期2011年7月21日 优先权日2011年7月21日
发明者卢德唐, 张龙军, 李道伦, 查文舒 申请人:中国科学技术大学