CN104112285A - 面向油田勘探与开发的沉积相图的智能绘制方法与系统 - Google Patents

Abstract

本发明实施例公开了一种面向油田勘探与开发的沉积相图的智能绘制方法与系统。所述方法包括：获取井点坐标及其对应相值，并根据河道流向和连通程度，对所述井点坐标进行旋转压缩，生成单井相平面分布图；在所述单井相平面分布图中，对所述井点进行带有约束的Delauney三角剖分，生成主连通线和辅连通线；沿着每一相的辅连通线追踪边界，生成每一相对应的边界多边形；对所述每一相对应的边界多边形进行曲线化，生成对应的相带边界；将带有所述相带边界的单井相平面分布图的坐标变换为初始坐标，生成对应的沉积相图。

Description

面向油田勘探与开发的沉积相图的智能绘制方法与系统

技术领域

本发明涉及石油勘探领域，尤其涉及一种面向油田勘探与开发的沉积相图的智能绘制方法及系统。

背景技术

在油田地质建模过程中，对沉积相进行建模是一个很重要的步骤。在实践中，对沉积相的建模是按照如下步骤进行的：

1、对测井曲线进行分析，确定整个区块的古沉积环境，并划分地层；

2、对分层后的测井曲线进行沉积相的判定；

3、综合所有测井的分析数据，对于每个地层，根据该地层所有测井位置的沉积相绘制出沉积相图。

在油田现场工作中，现有的绘制沉积相图方法往往通过人工手绘的方式进行的，因此效率较低，编辑较难，且不易进行数字化存档。而现有的由程序自动生成沉积相图的方法是按照如下步骤进行的：

1、选定要作图的地层，输入该层井点坐标数据以及对应的沉积相；

2、加入约束，自动绘图。

为了让自动成图的效果接近人工绘制的效果，需要在成图前加入一些约束。原因是，由测井曲线所分析得到的沉积相分为河道相，席状砂相，表外相等等，对应的是古地理环境中河流的不同位置产生的沉积物。因此在绘制沉积相平面图的时候需要对地理环境进行多方位的考虑。具体的约束条件有：(1)古河流流向约束，这一条件一般是由地质专家给出。(2)河道宽度约束，也是由地质专家给出。(3)相序约束，这一约束是为了保证在不同相之间交接地带的分界线走势符合地质认识。(4)河道连通程度约束，也就是距离多远的河道井点被视为同一条河道。

3、用户对自动成图结果进行编辑，查缺补漏，以达到满意的效果。

而现有的自动生成沉积相图的方法有以下两种：

（一）利用Delaunay三角网剖分绘制沉积相图，其按照以下步骤进行：

1、对地层中所有井点形成的点集进行Delaunay剖分；

2、对这个Delaunay剖分的每一条边，分析它的两个端点是否是同一个沉积相。如果不是，就将这条边的中点加入到沉积相边界点的集合中，称为边界控制点；

3、根据Delaunay三角网的三角形邻接关系，对每个相区域进行追踪。每一条边界线都是由一系列边界控制点组成，之后通过在相邻边界控制点之间加入贝塞尔曲线得到光滑的相边界线；

4、用户可以通过添加虚拟井点和修改边界线的每段贝塞尔曲线的控制点的方法进行编辑。

但是，上述的这种利用Delaunay三角网剖分绘制相图的方法，无法保证河道的流向和连通性；且取连线中点作为边界线的控制点，无法控制河宽；直接用贝塞尔曲线连接相邻的控制点，在曲度过大的地方会导致边界线相交；添加虚拟井+编辑贝塞尔曲线控制点的编辑方法效率较低。

（二）利用指示克里金插值生成沉积相图，其按照以下步骤进行：

1、在整个区域建立平面规则网格系统；

2、在规则网格系统之上进行指示克里金插值。井点处的相值作为克里金插值的控制点，河流的流向作为插值的变差函数模型的各向异性指数；

3、在取得插值结果之后，根据每个网格点的相值来追踪相边界，并进行平滑处理。

上述的这种绘制沉积相图的方法，由于使用了克里金插值的方法，所以能在一定程度上满足河道流向的约束要求。但是还有以下的一些缺点：使用指示克里金插值会造成相带区域形态不好，这种方法不能保证河道宽度，并且在没有井点控制的区域会出现孤立的相区域；且无法满足相交界地区的相序约束；这种方法由于是一次成图，基本上很难进行编辑操作。

发明内容

本发明的目的是为了克服现有技术中存在的不足，提供一种面向油田勘探与开发的沉积相图的智能绘制方法及系统，利用已知井点处的相值来生成平面沉积相图，使之满足各项地质约束，并能够让用户通过此方便的进行沉积相图的编辑。

为了达到上述目的，本发明实施例公开了一种面向油田勘探与开发的沉积相图的智能绘制方法，包括：获取井点坐标及其对应相值，并根据河道流向和连通程度，对所述井点坐标进行旋转压缩，生成单井相平面分布图；在所述单井相平面分布图中，对所述井点进行带有约束的Delauney三角剖分，生成主连通线和辅连通线；沿着每一相的辅连通线追踪边界，生成每一相对应的边界多边形；对所述每一相对应的边界多边形进行曲线化，生成对应的相带边界；将带有所述相带边界的单井相平面分布图的坐标变换为初始坐标，生成对应的沉积相图。

为了达到上述目的，本发明实施例还公开了一种面向油田勘探与开发的沉积相图的智能绘制系统，包括：单井相平面分布图生成单元，用于获取井点坐标及其对应相值，并根据河道流向和连通程度，对所述井点坐标进行旋转压缩，生成单井相平面分布图；三角剖分单元，用于在所述单井相平面分布图中，对所述井点进行带有约束的Delauney三角剖分，生成主连通线和辅连通线；边界多边形生成单元，用于沿着每一相的辅连通线追踪边界，生成每一相对应的边界多边形；相带边界生成单元，用于对所述每一相对应的边界多边形进行曲线化，生成对应的相带边界；坐标变换单元，用于将带有所述相带边界的单井相平面分布图的坐标变换为初始坐标，生成对应的沉积相图。

本发明实施例的面向油田勘探与开发的沉积相图的智能绘制方法与系统，生成的沉积相图满足流向、河宽以及连通程度等的约束，生成的沉积相图比较美观，符合现场工作要求，并且用该方法生成的沉积相图的编辑比较方便，用户通过简单的操作就可以很好地优化成图效果。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例的一种面向油田勘探与开发的沉积相图的智能绘制方法的方法流程图；

图2为本发明实施例的一种面向油田勘探与开发的沉积相图的智能绘制系统的结构示意图；

图3为图2所示实施例中的单井相平面分布图生成单元101的结构示意图；

图4为图2所示实施例中的三角剖分单元102的结构示意图；

图5为图2所示实施例中的边界多边形生成单元103的结构示意图；

图6为本发明具体实施例中的生成的单井相平面分布图；

图7为本发明具体实施例中的标记主连通线和辅连通线的示意图；

图8为本发明具体实施例中生成的边界多边形；

图9为本发明具体实施例中生成的相带边界；

图10为对本发明生成的沉积相图进行添加井点操作时的效果示意图；

图11为对本发明生成的沉积相图进行删除井点操作时的效果示意图；

图12为对本发明生成的沉积相图进行添加主连通线操作时的效果示意图；

图13为对本发明生成的沉积相图进行删除主连通线操作时的效果示意图；

图14为对本发明生成的沉积相图进行边界线微调操作时的效果示意图；

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

图1为本发明实施例的一种面向油田勘探与开发的沉积相图的智能绘制方法的方法流程图。在本实施例中，该方法包括：步骤S101，获取井点坐标及其对应相值，并根据河道流向和连通程度，对所述井点坐标进行旋转压缩，生成单井相平面分布图；步骤S102，在所述单井相平面分布图中，对所述井点进行带有约束的Delauney三角剖分，生成主连通线和辅连通线；步骤S103，沿着每一相的辅连通线追踪边界，生成每一相对应的边界多边形；步骤S104，对所述每一相对应的边界多边形进行曲线化，生成对应的相带边界；步骤S105，将带有所述相带边界的单井相平面分布图的坐标变换为初始坐标，生成对应的沉积相图。

在本实施例的步骤S101中，所述获取井点坐标及其对应相值，并根据河道流向和连通程度，对所述井点坐标进行旋转压缩，生成单井相平面分布图，包括：根据所述河道流向，将所述井点坐标进行旋转使得所述河道流向的方向为将所述井点坐标进行旋转使得所述河道流向的方向为Y轴正方向；以及沿着所述河道流向对所述单井相平面分布图进行压缩，压缩比例由所述连通程度决定。而连通程度、河道流向这些约束参数均是由地质专家设定，可以直接获取到的约束。在本实施例的步骤S101中，生成单井相平面分布图之前，如果单井相平面分布图没有外边界，则自动创建所述单井相平面分布图的边界。

在本实施例的步骤S102中，所述在所述单井相平面分布图中，对所述井点进行带有约束的Delauney三角剖分，生成主连通线和辅连通线，包括：在距离小于设定值的河道相井点间加入主连通线约束，使得所述距离小于设定值的河道相井点间有边连接，生成主连通线；根据每个相的相序标号大小，对每个相标记辅连通线。相序标号是这么设定的：河道相序标号最大，表外最小，标记为0，其余在中间按顺序排列。在创建某个相N(>=1)的边界时，将相序标号大于等于N的都视为相N，其余相视为相0。此时，三角剖分中得到的边之中，两端点均为相N的称为主连通线。两端点一个是相N一个是相0的称为辅连接线。

在本实施例的步骤S103中，述沿着每一相的辅连通线追踪边界，生成每一相对应的边界多边形，包括：在所述每一相的每一条辅连通线上根据河宽约束和边界光滑性选取控制点；将所述控制点顺序连接起来，生成边界多边形。

在本实施例的步骤S104中，对所述每一相对应的边界多边形进行曲线化，生成对应的相带边界，包括：沿着同一个旋转方向将所述边界多边形的边标记为E(1)，E(2)，…，E(n)，之后将E(k)的中点和E(k+1)的中点作为贝塞尔曲线的节点，E(k)的3/4分点和E(k+1)的1/4分点作为控制点，生成贝塞尔曲线，所述贝塞尔曲线为所述相带边界。这样做可以使得生成的边界曲线平滑，且因为是边界多边形的内切曲线，所以不会与其他曲线相交。

在本实施例中，在生成所述沉积相图后，可对所述沉积相图进行编辑，所述编辑包括：添加井点、删除井点，添加主连通线，删除主连通线以及边界线的微调等。

图2为本发明实施例的一种面向油田勘探与开发的沉积相图的智能绘制系统的结构示意图。如图所示，本实施例的智能绘制系统包括：

单井相平面分布图生成单元101，用于获取井点坐标及其对应相值，并根据河道流向和连通程度，对所述井点坐标进行旋转压缩，生成单井相平面分布图；三角剖分单元102，用于在所述单井相平面分布图中，对所述井点进行带有约束的Delauney三角剖分，生成主连通线和辅连通线；边界多边形生成单元103，用于沿着每一相的辅连通线追踪边界，生成每一相对应的边界多边形；相带边界生成单元104，用于对所述每一相对应的边界多边形进行曲线化，生成对应的相带边界；坐标变换单元105，用于将带有所述相带边界的单井相平面分布图的坐标变换为初始坐标，生成对应的沉积相图。

在本实施例中，如图3所示，所述单井相平面分布图生成单元101包括：

旋转单元1011，用于根据所述河道流向，将所述井点坐标进行旋转使得所述河道流向的方向为将所述井点坐标进行旋转使得所述河道流向的方向为Y轴正方向；以及压缩单元1012，用于沿着所述河道流向对所述单井相平面分布图进行压缩，压缩比例由所述连通程度决定。

所述单井相平面分布图生成单元101还包括：边界创建单元1013，用于在生成单井相平面分布图之前，自动创建所述单井相平面分布图的边界。

在本实施例中，如图4所示，所述三角剖分单元102包括：主连通线生成单元1021，用于在距离小于设定值的河道相井点间加入主连通线约束，使得所述距离小于设定值的河道相井点间有边连接，生成主连通线；辅连通线生成单元1022，用于根据每个相的相序标号大小，对每个相标记辅连通线。

在本实施例中，如图5所示，所述边界多边形生成单元103包括：控制点选取单元1031，用于在所述每一相的每一条辅连通线上根据河宽约束和边界光滑性选取控制点；连接单元1032，用于将所述控制点顺序连接起来，生成边界多边形。

在本实施例中，所述相带边界生成单元104对所述每一相对应的边界多边形进行曲线化，生成对应的相带边界，包括：沿着同一个旋转方向将所述边界多边形的边标记为E(1)，E(2)，…，E(n)，之后将E(k)的中点和E(k+1)的中点作为贝塞尔曲线的节点，E(k)的3/4分点和E(k+1)的1/4分点作为控制点，生成贝塞尔曲线，所述贝塞尔曲线为所述相带边界。

利用本发明实施例的面向油田勘探与开发的沉积相图的智能绘制系统生成的沉积相图，满足流向、河宽以及连通程度等的约束，且生成的沉积相图比较美观，符合现场工作要求，编辑比较方便，用户通过简单的操作就可以很好地优化成图效果。

具体实施例：

步骤1、预处理，生成单井相平面分布图

(1)、获取井点位置坐标以及对应相值，输入相序标号信息；其中，相序标号信息是如此设定的：河道相序标号最大，表外最小，标记为相0，其余在中间按顺序排列。

(2)、根据河道流向和连通程度的要求，对井点进行坐标变换。方法是先将井点坐标旋转，使得河道流向为y轴正方向。随后沿着河道流向进行压缩，压缩比例由连通程度决定。

(3)、如果没有外边界，则自动创建边界。

该步骤1生成的单井相平面分布图如图6所示。其中，实心方框A表示河道相，三角形B1、圆圈B2、空心方框C依次表示距离中心越来越远，纱厚越来越薄，直至到尖灭。

步骤2、根据约束生成三角剖分

(1)、在距离小于某值的河道相井点间加入主连通线约束，在此约束的基础上对压缩后的井点进行Delauney三角剖分。也就是说在这个三角剖分中，距离小于某值的河道相井点之间是有边连接的；

(2)、根据每个相的相序关系，从大到小对每个相标记辅连通线。

具体方法是：在创建某个相N(>=1)的边界时，将相序标号大于等于N的都视为相N，其余相视为相0。此时，三角剖分中得到的边之中，两端点均为相N的称为主连通线。两端点一个是相N一个是相0的称为辅连接线，如图7所示。

步骤3、根据相序从大到小，对每个相生成相带边界

(1)、沿着当前追踪相的辅连接线来追踪边界。每条边界线的追踪结果是一系列辅连接线的标号集；

(2)、在每一条辅连接线上面根据河宽约束和边界光滑性考虑来选取控制点；

(3)、对于每条边界，由辅连接线集合得到对应的控制点集合。将控制点顺序连接起来组成边界多边形，如图8所示；

(4)、对边界多边形进行曲线化。方法是：沿着同一个旋转方向将边界多边形的边标记为E(1),E(2)…,E(n).之后将E(k)的中点和E(k+1)的中点作为贝塞尔曲线的节点，E(k)的3/4分点和E(k+1)的1/4分点作为控制点，生成贝塞尔曲线。这样做的好处是生成的边界曲线平滑，且因为是边界多边形的内切曲线，所以不会与其他曲线相交。

该步骤3后生成的相带边界如图9所示。

步骤4、在所有相带边界生成之后，将整张图的坐标化为初始坐标。并显示在屏幕上，供用户观察和编辑。该步骤即是步骤1中的旋转压缩的反变换过程，目的是把经过旋转压缩后的坐标转换为初始坐标。

在经过以上4步以后，沉积相图即自动生成。生成以后，用户可以根据成图效果对沉积相图进行编辑。编辑操作包括两类。

(1)、添加/删除井点，添加/删除主连通线。编辑结果会通过对Delaunay剖分进行局部更新的方式体现出来。编辑之后的结果返回步骤2，然后从步骤2开始重新运行，局部更新三角剖分，并重新追踪相带边界曲线。

a)添加井点：找到新加顶点所在三角形，用新加入的点剖分此三角形，其效果如图10所示。

b)删除井点：找到所有包含被删除点的三角形，这些三角形组成一个多边形。在清除这个多边形里面所有三角形，然后对这个多边形重新进行Delaunay剖分，其效果如图11所示。

c)添加主连通线：找到主连通线通过的三角形。这些三角形组成的多边形被添加进来的主连通线分为两个部分。对这两个部分进行Delaunay剖分，其效果如图12所示。

d)删除主连通线：分为两种情况。如果被删除这条主连通线是在编辑过程中后添加的，那么只需删除这条线，并将添加它时所影响到的三角形恢复即可。如果不是编辑过程中添加的，那么可以将主连通线标记为无效的主连通线，并在之后的处理中将其视为辅连接线，其效果如图13所示。

(2)、边界线的微调。用户可以通过调整辅连接线上贝塞尔曲线节点的位置来对边界曲线进行微调。微调结果直接反映在成果图中，其效果如图14所示。注意如果先进行微调后进行局部更新，被更新的区域中之前的微调结果会消失。

本发明实施例的面向油田勘探与开发的沉积相图的智能绘制方法与系统，生成的沉积相图满足流向、河宽以及连通程度等的约束，生成的沉积相图比较美观，符合现场工作要求，并且用该方法生成的沉积相图的编辑比较方便，用户通过简单的操作就可以很好地优化成图效果。

本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质（包括但不限于磁盘存储器、CD-ROM、光学存储器等）上实施的计算机程序产品的形式。

本发明是参照根据本发明实施例的方法、设备（系统）、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (13)

1.一种面向油田勘探与开发的沉积相图的智能绘制方法，其特征在于，所述方法包括：

获取井点坐标及其对应相值，并根据河道流向和连通程度，对所述井点坐标进行旋转压缩，生成单井相平面分布图；

在所述单井相平面分布图中，对所述井点进行带有约束的Delauney三角剖分，生成主连通线和辅连通线；

沿着每一相的辅连通线追踪边界，生成每一相对应的边界多边形；

对所述每一相对应的边界多边形进行曲线化，生成对应的相带边界；

将带有所述相带边界的单井相平面分布图的坐标变换为初始坐标，生成对应的沉积相图。

2.根据权利要求1所述的面向油田勘探与开发的沉积相图的智能绘制方法，其特征在于，所述获取井点坐标及其对应相值，并根据河道流向和连通程度，对所述井点坐标进行旋转压缩，生成单井相平面分布图，包括：

根据所述河道流向，将所述井点坐标进行旋转使得所述河道流向的方向为将所述井点坐标进行旋转使得所述河道流向的方向为Y轴正方向；以及

沿着所述河道流向对所述单井相平面分布图进行压缩，压缩比例由所述连通程度决定。

3.根据权利要求1所述的面向油田勘探与开发的沉积相图的智能绘制方法，其特征在于，生成单井相平面分布图之前，还包括自动创建所述单井相平面分布图的边界。

4.根据权利要求1所述的面向油田勘探与开发的沉积相图的智能绘制方法，其特征在于，所述在所述单井相平面分布图中，对所述井点进行带有约束的Delauney三角剖分，生成主连通线和辅连通线，包括：

在距离小于设定值的河道相井点间加入主连通线约束，使得所述距离小于设定值的河道相井点间有边连接，生成主连通线；

根据每个相的相序标号大小，对每个相标记辅连通线。

5.根据权利要求1所述的面向油田勘探与开发的沉积相图的智能绘制方法，其特征在于，所述沿着每一相的辅连通线追踪边界，生成每一相对应的边界多边形，包括：

在所述每一相的每一条辅连通线上根据河宽约束和边界光滑性选取控制点；

将所述控制点顺序连接起来，生成边界多边形。

6.根据权利要求1所述的面向油田勘探与开发的沉积相图的智能绘制方法，其特征在于，对所述每一相对应的边界多边形进行曲线化，生成对应的相带边界，包括：沿着同一个旋转方向将所述边界多边形的边标记为E(1)，E(2)，…，E(n)，之后将E(k)的中点和E(k+1)的中点作为贝塞尔曲线的节点，E(k)的3/4分点和E(k+1)的1/4分点作为控制点，生成贝塞尔曲线，所述贝塞尔曲线为所述相带边界。

7.根据权利要求1所述的面向油田勘探与开发的沉积相图的智能绘制方法，其特征在于，在生成所述沉积相图后，可对所述沉积相图进行编辑，所述编辑包括：添加井点、删除井点，添加主连通线，删除主连通线以及边界线的微调。

8.一种面向油田勘探与开发的沉积相图的智能绘制系统，其特征在于，所述智能绘制系统包括：

单井相平面分布图生成单元，用于获取井点坐标及其对应相值，并根据河道流向和连通程度，对所述井点坐标进行旋转压缩，生成单井相平面分布图；

三角剖分单元，用于在所述单井相平面分布图中，对所述井点进行带有约束的Delauney三角剖分，生成主连通线和辅连通线；

边界多边形生成单元，用于沿着每一相的辅连通线追踪边界，生成每一相对应的边界多边形；

相带边界生成单元，用于对所述每一相对应的边界多边形进行曲线化，生成对应的相带边界；

坐标变换单元，用于将带有所述相带边界的单井相平面分布图的坐标变换为初始坐标，生成对应的沉积相图。

9.根据权利要求8所述的面向油田勘探与开发的沉积相图的智能绘制系统，其特征在于，所述单井相平面分布图生成单元包括：

旋转单元，用于根据所述河道流向，将所述井点坐标进行旋转使得所述河道流向的方向为将所述井点坐标进行旋转使得所述河道流向的方向为Y轴正方向；以及

压缩单元，用于沿着所述河道流向对所述单井相平面分布图进行压缩，压缩比例由所述连通程度决定。

10.根据权利要求8所述的面向油田勘探与开发的沉积相图的智能绘制系统，其特征在于，所述单井相平面分布图生成单元还包括：

边界创建单元，用于在生成单井相平面分布图之前，自动创建所述单井相平面分布图的边界。

11.根据权利要求8所述的面向油田勘探与开发的沉积相图的智能绘制系统，其特征在于，所述三角剖分单元包括：

主连通线生成单元，用于在距离小于设定值的河道相井点间加入主连通线约束，使得所述距离小于设定值的河道相井点间有边连接，生成主连通线；

辅连通线生成单元，用于根据每个相的相序标号大小，对每个相标记辅连通线。

12.根据权利要求8所述的面向油田勘探与开发的沉积相图的智能绘制系统，其特征在于，所述边界多边形生成单元包括：

控制点选取单元，用于在所述每一相的每一条辅连通线上根据河宽约束和边界光滑性选取控制点；

连接单元，用于将所述控制点顺序连接起来，生成边界多边形。

13.根据权利要求8所述的面向油田勘探与开发的沉积相图的智能绘制系统，其特征在于，所述相带边界生成单元对所述每一相对应的边界多边形进行曲线化，生成对应的相带边界，包括：沿着同一个旋转方向将所述边界多边形的边标记为E(1)，E(2)，…，E(n)，之后将E(k)的中点和E(k+1)的中点作为贝塞尔曲线的节点，E(k)的3/4分点和E(k+1)的1/4分点作为控制点，生成贝塞尔曲线，所述贝塞尔曲线为所述相带边界。