一种基于空间矢量的砂体储层构型建模方法及系统的制作方法
【专利摘要】本发明提供一种基于空间矢量的砂体储层构型建模的方法及系统，所述方法包括：获取当前砂体的测井数据；对测井数据进行数字处理，得到测积层数据、废弃河道数据以及点坝砂体数据；对废弃河道数据进行拟合，得到平面分布的深泓线；根据深泓线以及点坝砂体数据建立废弃河道三维几何体；根据反曲曲线公式、侧积层数据确定侧积层的空间分布体；根据废弃河道三维几何体、侧积层的空间分布体以及点坝砂体数据确定矢量化的砂体储层构型模型；对矢量化的砂体储层构型模型进行网格化，得到网格化的砂体储层构型模型。本发明的模拟结果可以刻画不同规模层次的非均质体，可以快速的对模型进行调节。
【专利说明】一种基于空间矢量的砂体储层构型建模方法及系统
【技术领域】
[0001]本发明关于油田开采【技术领域】，特别是关于储层的空间建模技术，具体的讲是一种基于空间矢量的砂体储层构型建模方法及系统。
【背景技术】
[0002]三维储层建模，即建立储层特征三维分布的数字化模型，其本质是基于三维网格表征储层特征的分布，其成果是三维数据体。这一技术是上世纪80年代随着计算机技术的发展而发展起来的。基于计算机存储和显示技术，将储层三维网块化后，对各个网格(grid)赋以各自的储层参数值，并按三维空间分布位置存入计算机内，形成了三维数据体，这样就可以进行储层的三维显示，可以任意切片和切剖面(不同层位、不同方向剖面)，以及进行各种运算和分析。对于一个给定的储层，按X，Y，Z三个方向划分网格，将储层在三维空间划分成众多网格，即为三维网格化(如图13所示)。
[0003]在储层地质学的研究中，对于储层在三维空间的表征主要是应用地质统计学的理论和方法，综合多学科的信息对储层非均质的空间分布进行确定性和随机性的赋值。确定性的途径主要应用各种克里金方法，随机性的途径则主要是应用各种高斯模拟(如截断高斯模拟、序贯指示模拟)方法。按照随机模拟基本单元不同，随机建模方法可以分为基于目标的方法和基于象元的方法两类。
[0004]基于象元的方法的模拟基本单元为象元(cell)，即网格化储层模型中单个网格(grid)，借助地质统计学中的变差函数，以不同变差函数表征变量的空间变异特征，逐网格分配储层属性值。基于象元的模拟方法中，对于河流相应用较多的方法主要为序贯指示模拟。尽管序贯指示模拟应用相当广泛，然而由于它是以变差函数为依据来反映储层形态及空间分布特征，而变差函数只考虑了两点相关性。地质现象复杂性决定了对其描述应该是基于更多点信息的，而变差函数很难真实反映储层地下结构及形态，难于再现复杂地质的界面信息。
[0005]基于目标的方法的模拟基本单元为目标物体(即具有离散性质的地质特征，如沉积相、流动单元等)，通过对目标几何形态(如长、宽、厚及其之间定量关系)的定量研究，在建模过程中通过对多个网格同时赋予生目标物体而建立储层的三维模型。但是，依沉积学原理，沉积砂体的横向宽度、厚度受其沉积时的水动力条件制约，两者之间存在一定的定量统计关系，基于目标的算法对沉积因素考虑不够。
[0006]另外的学者则对基于变差函数的两点统计学储层建模的改进，提出了多点地质统计学。应用“数据事件”代替变差函数表达地质变量的空间结构性，克服传统地质统计学不能再现目标几何形态的不足。但是多点地质统计学作为一种新的建模思路和方法，在再现河道的连续性方面仍需改进，并且对各种非平稳现象(具有趋势分布的)的模拟方面仍然存在困难，尤其在对地质体的三维形态表达和模拟方面仍需改进。
[0007]储层构型建模则多采用序贯指示模拟与人机结合再处理的手段，特点是工作量大，且无法表征模型中不确定的因素。也有学者提出了二次网格加密方法、基于沉积过程的模、层次约束下的随机抽样等随机建模方法。
[0008]但各种方法仍然沿用了传统的基于象元的建模方法。即首先定义网格，之后将构型要素作为一种离散属性定义在网格之上。因此在网格化研究区时必须考虑到所定义的网格能予以刻画最关键的非均质体。对曲流河储集层构型模型而言，虽然点坝内部的侧积层厚度较小(0.1m?1.0m级)，但却具有一定规模的横向延伸范围，是造成点坝砂体剩余油富集的主要因素。在精细储层构型建模中，为了刻画此种规模的非均质体，选择之一便是细化网格，这样势必会造成网格节点个数的膨胀。虽然计算机有了很大的发展，对于一个纵向、横向上划分的网格精度还是有一定的限制，并不能采取任意细化的网格。例如，对于一个IkmX IkmX 50m的模型，如果使用5mX 5mX Im的建模网格，网格总数将达到200 X IO4个。即使这样的网格依然无法刻画垂向上小于0.5m的非均值体。因此，在建模阶段采用基于网格的方法导致网格数目急剧增加，增加了计算时间，也限制了随机模拟实现的个数。虽然已有随机建模算法对高级别构型要素，如单一河道砂体的三维建模有较好的效果，但难于满足砂体内部构型建模解剖与建模的需要。尤其对于微小级次的构型单元，人工绘制进行三维建模的工作量巨大，效率低下。
【发明内容】
[0009]为了克服现有技术存在的上述问题，本发明提供了一种基于空间矢量的砂体储层构型建模的方法及系统，对模拟地质体的定义采用了基于实数集的点、线、面所限定的体，因而模型在模拟过程中不受网格尺寸的限制，使得模拟结果可以刻画不同规模层次的非均质体，可以快速的对模型进行调节，在模拟过程中只是对地质体的空间位置和形态的定义参数进行调节。
[0010]本发明的目的之一是，提供一种基于空间矢量的砂体储层构型建模的方法，包括:获取当前砂体的测井数据；对所述的测井数据进行数字处理，得到侧积层数据、废弃河道数据以及点坝砂体数据；对所述的废弃河道数据进行拟合，得到平面分布的深泓线；根据所述的深泓线以及点坝砂体数据建立废弃河道三维几何体；根据反曲曲线公式、侧积层数据确定侧积层的空间分布体；根据所述废弃河道三维几何体、侧积层的空间分布体以及点坝砂体数据确定矢量化的砂体储层构型模型；对所述的矢量化的砂体储层构型模型进行网格化，得到网格化的砂体储层构型模型。
[0011]本发明的目的之一是，提供了一种基于空间矢量的砂体储层构型建模的系统，包括:测井数据获取装置，用于获取当前砂体的测井数据；测井数据处理装置，用于对所述的测井数据进行数字处理，得到侧积层数据、废弃河道数据以及点坝砂体数据；拟合装置，用于对所述的废弃河道数据进行拟合，得到平面分布的深泓线；三维几何体建立装置，用于根据所述的深泓线以及点坝砂体数据建立废弃河道三维几何体；空间分布体建立装置，用于根据反曲曲线公式、侧积层数据确定侧积层的空间分布体；矢量化模型确定装置，用于根据所述废弃河道三维几何体、侧积层的空间分布体以及点坝砂体数据确定矢量化的砂体储层构型模型；网格化装置，用于对所述的矢量化的砂体储层构型模型进行网格化，得到网格化的砂体储层构型模型。
[0012]本发明的有益效果在于，提供了一种基于空间矢量的砂体储层构型建模的方法及系统，与现有技术相比，该方法提出的主要目的是在储层构型建模中描述不同级别的非均质体以及准确刻画各级地质体的复杂空间形态。在模拟阶段的地质体的存储借鉴了矢量图形存储的思路，即对模拟地质体的定义采用了基于实数集的点、线、面所限定的体，因而模型在模拟过程中不受网格尺寸的限制，使得模拟结果可以刻画不同规模层次的非均质体。应用基于矢量体存储方法的另外一个优点是可以快速的对模型进行调节，因为在模拟过程中只是对地质体的空间位置和形态的定义参数进行调节，因而这样的调节方式较基于网格的调节更便捷，从而使结果既忠实条件数据点又易再现模拟目标地质体的几何形态。
[0013]为让本发明的上述和其他目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附图式，作详细说明如下。
【专利附图】
【附图说明】
[0014]为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0015]图1为本发明实施例提供的一种基于空间矢量的砂体储层构型建模的方法的流程图；
[0016]图2为图1中的步骤S104的具体流程图；
[0017]图3为图1中的步骤S105的具体流程图；
[0018]图4为图1中的步骤S106的具体流程图；
[0019]图5为图1中的步骤S107的具体流程图；
[0020]图6为图5中的步骤S503的具体流程图；
[0021]图7为本发明实施例提供的一种基于空间矢量的砂体储层构型建模的系统的示意图；
[0022]图8为本发明提供的一种基于空间矢量的砂体储层构型建模的系统中三维几何体建立装置400的结构框图；
[0023]图9为本发明提供的一种基于空间矢量的砂体储层构型建模的系统中空间分布体建立装置500的结构框图；
[0024]图10为本发明提供的一种基于空间矢量的砂体储层构型建模的系统中矢量化模型确定装置600的结构框图；
[0025]图11为本发明提供的一种基于空间矢量的砂体储层构型建模的系统中网格化装置700的结构框图；
[0026]图12为本发明提供的一种基于空间矢量的砂体储层构型建模的系统中数值填充模块703的结构框图；
[0027]图13为现有技术中的三维网格化示意图；
[0028]图14为曲流河的河曲变化形态实施方式一的示意图；
[0029]图15为曲流河的河曲变化形态实施方式二的示意图；
[0030]图16为曲流河的河曲变化形态实施方式三的示意图；
[0031]图17为废弃河道的空间矢量表征方法得到的废弃河道剖面形态示意图；
[0032]图18为废弃河道的空间矢量表征方法得到的废弃河道剖面沿着深泓线的移动的三维几何体示意图；
[0033]图19为侧积层剖面的空间矢量表达方法对应的理论反曲曲线形态示意图；
[0034]图20为侧积层剖面的空间矢量表达方法对应的平移反曲曲线所圈定的侧积层剖面示意图；
[0035]图21为将侧积层剖面对应的两个反曲曲线沿两个圆相交的圆弧线进行平行移动侧积层的空间分布体示 意图；
[0036]图22为将点坝砂体数据填充到图18的废弃河道三维几何体的区域中得到的第一填充三维几何体的示意图；
[0037]图23为矢量化的砂体储层构型模型的剖面图；
[0038]图24为依据水平密度填充后得到的矢量化的砂体储层构型模型的剖面图；
[0039]图25为网格化后的单个点坝砂体构型模拟示意图。
【具体实施方式】
[0040]下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0041]本发明提出一种基于空间矢量的砂体储层构型建模的方法，图1为该方法的具体流程图，由图1可知，所述的方法包括:
[0042]SlOl:获取当前砂体的测井数据，在具体的实施方式中，即获取每口井在当前建模层位的测井数据。
[0043]S102:对所述的测井数据进行数字处理，得到侧积层数据、废弃河道数据以及点坝砂体数据。在具体的实施方式中，对所述的测井数据进行数字处理，依据其形态确定砂体的构型类型，得到侧积层数据、废弃河道数据以及点坝砂体数据。
[0044]S103:对所述的废弃河道数据进行拟合，得到平面分布的深泓线。
[0045]虽然自然界的废弃河道形态呈现出各种形态，但是用数学方法去刻画它则主要考虑到其共性，通过设置不同的参数则可以描述不同的河曲形态。本发明中对于废弃河道的矢量表达是采用了 2部分参数。第I部分是河曲控制点，共为7个，可用p0、pl、……p6来表示。第2部分是河曲形态参数，可用α来表示。
[0046]对于深泓线的空间形态表征是采用了关键点控制下的三次样条内插法来进行表达。对于一个河曲采用了7个关键控制点(？0，？1，...，？6)和1个形态参数角α (如图14所示)。不同的河曲形态是通过设定形态控制关键点和形态参数角进行差值来确定其空间各个点的位置，图14和图15中是使用不同参数得到的两个不同形态的河曲形态例子。图14、图15、图16中的曲线(即p4、p3、p2、p1、ρ6组成的曲线)即为深泓线。
[0047]S104:根据所述的深泓线以及点坝砂体数据建立废弃河道三维几何体。废弃河道的三维几何体是通过沿着深泓线平移废弃河道的剖面来表征的。废弃河道的剖面形态中，最主要一点就是设定深泓线与剖面的交点，即河道剖面最大厚度点(图17中的a3)的相对位置。a3即为深泓线上的点。图2为步骤S104的具体流程图，由图2可知，该步骤具体包括:[0048]S201:从所述的点坝砂体数据中获取砂体厚度；
[0049]S202:将所述深泓线上的任一点作为数据点。如图17所示，选取的数据点为a3，砂体厚度即为图17中的H，即|a2a3|。
[0050]S203:根据所述的数据点以及曲河流点坝砂体经验公式确定河道宽度。当其宽深比确定后，侧积层的倾角就已经给定，而废弃河道深度H应该满足宽深比的要求。此处提及的曲河流点坝砂体经验公式包括:
[0051]①由点坝砂体厚度(hs)经压实校正求得原始砂体厚度，这一厚度也近似等于河道满岸深度(hr)。 [0052]②依据满岸深度由经验公式lgWfl.541g hr+0.83，计算得到满岸宽度(Wr)。
[0053]③利用公式Wm=7.44ffrL01计算出曲流带宽度(Wm)，并依公式Wd=6.5894ffr+0.0677计算点坝跨度(Wd)；
[0054]④根据经验规模Wl=2/3Wr计算出侧积层水平宽度(W1)，并使用公式β =33exp (-0.1ffZ/HZ)计算侧积层倾角(β )。
[0055]S204:根据所述的数据点、侧积层水平宽度、河道宽度以及砂体厚度确定废弃河道的剖面。如图17所示，数据点为a3，根据侧积层水平宽度、河道宽度可以确定出a1、a4，根据砂体厚度可以确定出a2，弧线aia2则可以使用以点&3为中心点的1/4椭圆的方法来刻画其形态，其中长轴和短轴分别为Iap3I和|a2a3|。同理对于弧线a2a4也是使用以点a3为中心点的1/4椭圆的方法来刻画其形态，其中长轴和短轴分别为|a3a4|和|a2a3|，如此即可确定出废弃河道的剖面(如图17所示)。
[0056]S205:遍历所述深泓线上的点，得到多个废弃河道的剖面。深泓线上的所有点皆可作为数据点构建剖面，如此得到多个剖面。
[0057]S206:多个废弃河道的剖面之间圈定的空间体积即为废弃河道的三维几何体。沿所模拟的河道深泓线在不同的位置计算河道空间的剖面形态后，各个面之间所圈定的空间体积即为废弃河道的三维几何体，如图18所示。在实际模拟中，河道深泓线的参数以及剖面形态等参数都可以依据地质家对于河曲形态的理解和研究区的具体问题进行设置。
[0058]由图1可知，该方法还包括:
[0059]S105:根据反曲曲线公式、侧积层数据确定侧积层的空间分布体。图3为步骤S105的具体流程图，由图3可知，该步骤具体包括:
[0060]S301:根据反曲曲线公式确定在实数范围内的反曲曲线形态。
[0061]曲流河构型模式中另一个重要的构型要素则是侧积层。在本发明中对于侧积层剖面形态描述采用了反曲曲线。反曲曲线的数学表达为
I
[0062]J ~—
1 + e
[0063]其在实数范围内的形状为图19所示。
[0064]S302:对所述的反曲曲线形态进行坐标平移或坐标旋转，得到侧积层剖面。在实际模拟中，可采用坐标平移或者坐标旋转的方法生成形态各异的侧积层剖面形态，如图20所示，是对一个反曲曲线进行平移所圈定的剖面，这个剖面构成了侧积层剖面。侧积层剖面包括着两个反曲曲线。
[0065]S303:从所述的侧积层数据中获取侧积层中心厚度、侧积层水平跨度。[0066]S304:设置两个圆心距为侧积层中心厚度、直径为侧积层水平跨度的圆，如图21中所示的两个圆，两个圆的圆心的距离Ic1C2I等于侧积层中心的厚度，两个圆的直径为侧积层水平跨度，而延展的距离则由图21中的Iap2I和|a2a3|来控制，则两个圆可以形成如图21所示的两个相交的圆弧线。
[0067]S305:将所述侧积层剖面对应的两个反曲曲线沿两个圆相交的圆弧线进行平行移动；
[0068]S306:平行移动所圈定的空间即为侧积层的空间分布体。如图21所示。模拟过程中是通过给定图21中的|aia2|和|a2a3|的分布函数从而控制侧积层侧向延伸的距离和延伸的方式。将两个反曲曲线沿着所定义好的圆弧线进行平行移动所圈定的空间体积即是侧积层的空间分布体积。
[0069]如上所述，本发明虽然模拟仍然采用了传统的基于目标体的模拟，但是并不设定模型的网格尺寸，而是直接在模型内用面和体的方法生成模拟地质体。模型只存储地质体的空间分布形态参数而不存储每个网格上的目标体的属性。故而必须用数学公式或形态参数刻画模型中的基本地质体。虽然地下的地质体形态是极其复杂的，但是这些复杂的地质体都可以用一些简单形状比如带状、球状、椭球状，楔形等经过一定的修正和调整来刻画。其中地质学家的储层构型概念地质模式是定义模拟基本单元的基础。故而本发明以现在的构型地质研究作为基础，以构型建模算法的开发为目标，对各类储层构型要素进行三维形态刻画和表征。
[0070]由图1可知，该方法还包括:
[0071]S106:根根据所述废弃河道三维几何体、侧积层的空间分布体以及点坝砂体数据确定矢量化的砂体储层构型模型。图4为步骤S106的具体流程图，由图4可知，该步骤具体包括:
[0072]S401:将所述的点坝砂体数据填充到所述废弃河道三维几何体所圈定的区域中，得到第一填充三维几何体，如图22所示，将点坝砂体数据填充到图18的废弃河道三维几何体的区域中得到第一填充三维几何体。
[0073]S402:从所述的点坝砂体数据中获取侧积层在井身上的位置、侧积层的发育角度以及侧积层的水平密度；
[0074]S403:根据所述侧积层在井身上的位置、侧积层的发育角度以及侧积层的水平密度，将所述侧积层的空间分布体填充到第一填充三维几何体，得到第二填充几何体。侧积层在测井上的位置用于表征此口井钻遇侧积层的位置(如图23示例中的LI和L2)，即必须于井点的构型要素解释类型相吻合，角度用于表征侧积层发育的角度(见图23中示例)。如图23所示，将侧积层的空间分布体分别填充到井A、井B，与点坝砂体表面相交的点为a、C。
[0075]S404:从所述的点坝砂体数据中获取侧积层的水平密度，水平密度用于表征侧积层发育的多少。
[0076]S405:根据所述的水平密度将将所述侧积层的空间分布体填充到所述第二填充几何体，得到矢量化的砂体储层构型模型。如图24所示，将侧积层的空间分布体分别填充到点坝砂体中，与点坝砂体表面的交点为b、d。图25即为网格化后的单个点坝砂体构型模拟示意图。
[0077]上述步骤利用所定义的基本构型要素地质体进行随机模拟。地下储层构型表征的核心是井间预测，而预测的关键是预知对象的分布规律或模式。虽然模型的调节只是对各个地质体的表征参数进行调节，但是约束依然是各级储层构型定量地质模式。如各个构型要素的展布规模和相互之间的切割和叠置关系。国内外进行过不少的储层地质知识库的研究，初步建立了垂向信息与侧向规模的关系，涉及了河流相、三角洲相、扇三角洲相、深水浊积相等，取得了丰富的研究成果，这些研究都是进行随机模拟的基础。故而本发明将现在的构型地质定量模式研究作为基础，以构型随机模拟算法的开发为目标，对各类储层构型要素的三维空间随机模拟方法进行研究，最终得到基于矢量表达的储层构型模型。
[0078]在定义好各个构型要素之后则是进行基于目标的随机模拟。与传统的基于目标体的模拟方法一样，认为目标中心点位置符合齐次泊松点过程，而在无井区，模拟采用非齐次的泊松点过程，从而满足了井间与井点分布具有差异的要求。模拟仍然采用迭代的算法，并且在模拟过程的必须受到各个目标体的叠置切割模式和井数据的约束。
[0079]地质构型分析所得到的单井构型解释是作为模拟阶段的硬数据，最终的构型模型必须满足井数据的约束。本发明在模拟程序中采用了优先拟合井点数据的做法。首先通过调整各个地质体的形态参数使其满足井点约束然。之后在选定的地质模式约束下进行空间其它构型要素的投放以满足目标函数。在点坝砂体内部模拟侧积层为例(如图24)，依据井数据模拟出了侧积层a和C。而依据侧积层间距等定量构型模式信息可以模拟出侧积层b和d。
[0080]由图1可知，该方法还包括:
[0081]S107:对所述的矢量化的砂体储层构型模型进行网格化，得到网格化的砂体储层构型模型。图5为步骤S107的具体流程图，由图5可知，该步骤包括:
[0082]S501:根据所述的侧积层数据确定网格尺寸，即网格必须满足刻画侧积层的最小厚度。在具体的实施方式中，如侧积层数据为0.8米厚，则确定的网格尺寸可以是0.5米。
[0083]S502:根据所述的网格尺寸对所述矢量化的砂体储层构型模型建立模拟区的网格；
[0084]S503:对所述模拟区的网格根据构型要素类型进行数值填充，得到网格化的砂体储层构型模型。图6为步骤S503的具体流程图，由图6可知，该步骤具体包括:
[0085]S601:判断所述模拟区的网格对应的构型要素类型是否仅为废弃河道；
[0086]S602:当判断为是时，将所述的网格填充为废弃河道数据；
[0087]S603:否则，判断所述模拟区的网格对应的构型要素类型是否仅为侧积层；
[0088]S604:当判断为是时，将所述的网格填充为侧积层数据；
[0089]S605:否则，判断所述模拟区的网格对应的构型要素类型是否仅为点坝砂体；
[0090]S606:当判断为是时，将所述的网格填充为点坝砂体数据；
[0091]S607:否则，确定所述网格对应的构型要素类型的比例。当网格中包括多种构型要素类型时，分别确定各个要素所占的比例，如:网格中既有点坝砂体，又有侧积层，则确定二者的比例，为2:5。
[0092]S608:获取所述比例中的最大值。
[0093]S609:确定所述最大值对应的构型要素类型；
[0094]S610:将所述的网格填充为最大值对应的构型要素类型对应的数据。
[0095]如上的实施例中，网格中既有点坝砂体，又有侧积层，二者的比例为2:5，则比例最大之为5，对应的构型要素类型为侧积层，因此该网格填充的数据为侧积层数据。如图25为单个点坝砂体构型模型网格化结果。
[0096]由于后期物性参数的建模方法仍然是基于网格的，而且后期的数值模拟也是基于离散网格数据体，故而模拟的最终结果应该是一个基于网格的模型。所以必须将矢量存储的空间地质体进行离散化，将各个储层构型属性赋于每一个网格。由于基于矢量存储的地质模型不受空间网格大小的限制，并且不再考虑模拟速度的要求，此阶段的网格定义可以采取任意细化以满足刻画不同规模的非均质地质体。故而如何对基于矢量表达的模型进行网格化处理也是本发明的内容之一。由于此阶段不早模拟，没有迭代等计算，所以可以选择能够刻画最小规模的构型要素的模型网格进行离散化，如图25为单个点坝砂体构型模型网格化结果。
[0097]如上即为本发明提供的方法，本发明还提出一种基于空间矢量的砂体储层构型建模的系统，图7为该系统的结构示意图，由图7可知，所述的系统包括:
[0098]测井数据获取装置100，用于获取当前砂体的测井数据，在具体的实施方式中，SP获取每口井在当前建模层位的测井数据。
[0099]测井数据处理装置200，用于对所述的测井数据进行数字处理，得到侧积层数据、废弃河道数据以及点坝砂体数据。在具体的实施方式中，对所述的测井数据进行数字处理，依据其形态确定砂体的构型类型，得到侧积层数据、废弃河道数据以及点坝砂体数据。此处的数字处理具体包括:依据测井相进行。
[0100]拟合装置300，用于对所述的废弃河道数据进行拟合，得到平面分布的深泓线。
[0101]虽然自然界的废弃河道形态呈现出各种形态，但是用数学方法去刻画它则主要考虑到其共性，通过设置不同的参数则可以描述不同的河曲形态。本发明中对于废弃河道的矢量表达是采用了 2部分参数。第I部分是河`曲控制点，共为7个，可用p0、pl、……p6来表示。第2部分是河曲形态参数，可用α来表示。
[0102]对于深泓线的空间形态表征是采用了关键点控制下的三次样条内插法来进行表达。对于一个河曲采用了7个关键控制点(？0，？1，...，？6)和1个形态参数角α (如图14所示)。不同的河曲形态是通过设定形态控制关键点和形态参数角进行差值来确定其空间各个点的位置，图14和图15中是使用不同参数得到的两个不同形态的河曲形态例子。图14、图15、图16中的曲线(即p4、p3、p2、p1、ρ6组成的曲线)即为深泓线。
[0103]三维几何体建立装置400，用于根据所述的深泓线以及点坝砂体数据建立废弃河道三维几何体。废弃河道的三维几何体是通过沿着深泓线平移废弃河道的剖面来表征的。废弃河道的剖面形态中，最主要一点就是设定深泓线与剖面的交点，即河道剖面最大厚度点(图17中的a3)的相对位置。a3即为深泓线上的点。图8为三维几何体建立装置400的结构框图，由图8可知，三维几何体建立装置具体包括:
[0104]砂体厚度获取模块401，用于从所述的点坝砂体数据中获取砂体厚度；
[0105]数据点选取模块402，用于将所述深泓线上的任一点作为数据点。如图17所示，选取的数据点为a3，砂体厚度即为图17中的H，即|a2a3|。
[0106]河道宽度确定模块403，用于根据所述的数据点以及曲河流点坝砂体经验公式确定凹岸水平宽度以及河道宽度。
[0107]当其宽深比确定后，侧积层的倾角就已经给定，而废弃河道深度H应该满足宽深比的要求。此处提及的曲河流点坝砂体经验公式包括:
[0108]①由点坝砂体厚度(hs)经压实校正求得原始砂体厚度，这一厚度也近似等于河道满岸深度(hr)。
[0109]②依据满岸深度由经验公式lgWfl.541g hr+0.83，计算得到满岸宽度(Wr)。
[0110]③利用公式Wm=7.44ffrL01计算出曲流带宽度(Wm)，并依公式Wd=6.5894ffr+0.0677计算点坝跨度(Wd)；
[0111]④根据经验规模Wl=2/3Wr计算出侧积层水平宽度(W1)，并使用公式β =33exp (-0.1ffZ/HZ)计算侧积层倾角(β )。
[0112]废弃河道剖面确定模块404，用于根据所述的数据点、侧积层水平宽度、河道宽度以及砂体厚度确定废弃河道的剖面。如图17所示，数据点为a3，根据侧积层水平宽度、河道宽度可以确定出apa4,根据砂体厚度可以确定出a2,弧线则可以使用以点a3为中心点的1/4椭圆的方法来刻画其形态，其中长轴和短轴分别为Iap3I和|a2a3|。同理对于弧线a2a4也是使用以点a3为中心点的1/4椭圆的方法来刻画其形态，其中长轴和短轴分别为
a3a4|和|a2a3|， 如此即可确定出废弃河道的剖面(如图17所示)。
[0113]数据点遍历模块405，用于遍历所述深泓线上的点，得到多个废弃河道的剖面。深泓线上的所有点皆可作为数据点构建剖面，如此得到多个剖面。
[0114]三维几何体确定模块406，多个废弃河道的剖面之间圈定的空间体积即为废弃河道的三维几何体。沿所模拟的河道深泓线在不同的位置计算河道空间的剖面形态后，各个面之间所圈定的空间体积即为废弃河道的三维几何体，如图18所示。在实际模拟中，河道深泓线的参数以及剖面形态等参数都可以依据地质家对于河曲形态的理解和研究区的具体问题进行设置。
[0115]由图7可知，该系统还包括:
[0116]空间分布体建立装置500，用于根据反曲曲线公式、侧积层数据确定侧积层的空间分布体。图9为空间分布体建立装置500的结构框图，由图9可知，空间分布体建立装置具体包括:
[0117]反曲曲线形态确定模块501，用于根据反曲曲线公式确定在实数范围内的反曲曲线形态。
[0118]曲流河构型模式中另一个重要的构型要素则是侧积层。在本发明中对于侧积层剖面形态描述采用了反曲曲线。反曲曲线的数学表达为
[0119]y=^~7
1 + e
[0120]其在实数范围内的形状为图19所示。
[0121]侧积层剖面确定模块502，用于对所述的反曲曲线形态进行坐标平移或坐标旋转，得到侧积层剖面。在实际模拟中，可采用坐标平移或者坐标旋转的方法生成形态各异的侧积层剖面形态，如图20所示，是对一个反曲曲线进行平移所圈定的剖面，这个剖面构成了侧积层剖面。侧积层剖面包括着两个反曲曲线。
[0122]水平跨度获取模块503，用于从所述的侧积层数据中获取侧积层中心厚度、侧积层水平跨度。
[0123]设置模块504，用于设置两个圆心距为侧积层中心厚度、直径为侧积层水平跨度的圆，
[0124]如图21中所示的两个圆，两个圆的圆心的距离Ic1C2I等于侧积层中心的厚度，两个圆的直径为侧积层水平跨度，而延展的距离则由图21中的|aia2|和|a2a3|来控制，则两个圆可以形成如图21所示的两个相交的圆弧线。
[0125]平行移动模块505，用于将所述侧积层剖面对应的两个反曲曲线沿两个圆相交的圆弧线进行平行移动；
[0126]空间分布体确定模块506，平行移动所圈定的空间即为侧积层的空间分布体。如图21所示。模拟过程中是通过给定图21中的|aia2|和|a2a3|的分布函数从而控制侧积层侧向延伸的距离和延伸的方式。将两个反曲曲线沿着所定义好的圆弧线进行平行移动所圈定的空间体积即是侧积层的空间分布体积。
[0127]如上所述，本发明虽然模拟仍然采用了传统的基于目标体的模拟，但是并不设定模型的网格尺寸，而是直接在模型内用面和体的方法生成模拟地质体。模型只存储地质体的空间分布形态参数而不存储每个网格上的目标体的属性。故而必须用数学公式或形态参数刻画模型中的基本地质体。虽然地下的地质体形态是极其复杂的，但是这些复杂的地质体都可以用一些简单形状比如带状、球状、椭球状，楔形等经过一定的修正和调整来刻画。其中地质学家的储层构型概念地质模式是定义模拟基本单元的基础。故而本发明以现在的构型地质研究作为基础，以构型建模算法的开发为目标，对各类储层构型要素进行三维形态刻画和表征。由图7可知，该系统还包括:
[0128]矢量化模型确定装置600，用于根根据所述废弃河道三维几何体、侧积层的空间分布体以及点坝砂体数据确定矢量化的砂体储层构型模型。图10为本发明提供的一种基于空间矢量的砂体储层构型建模的系统中矢量化模型确定装置600的结构框图，由图10可知，矢量化模型确定装置具体包括:
[0129]第一填充模块601，用于将所述的点坝砂体数据填充到所述废弃河道三维几何体的所圈定的区域中，得到第一填充三维几何体，如图22所示，将点坝砂体数据填充到图18的废弃河道三维几何体的区域中得到第一填充三维几何体。
[0130]角度获取模块602，用于从所述的点坝砂体数据中获取侧积层在井身上的位置、侧积层的发育角度以及侧积层的水平密度；
[0131]第二填充模块603，用于根据所述侧积层在井身上的位置、侧积层的发育角度以及侧积层的水平密度，将所述侧积层的空间分布体填充到第一填充三维几何体，得到第二填充几何体。侧积层在测井上的位置用于表征此口井钻遇侧积层的位置(如图23示例中的LI和L2)，即必须于井点的构型要素解释类型相吻合，角度用于表征侧积层发育的角度(见图23中示例)。如图23所示，将侧积层的空间分布体分别填充到井A、井B，与点坝砂体表面相交的点为a、C。
[0132]水平密度获取模块604，用于从所述的点坝砂体数据中获取侧积层的水平密度，水平密度用于表征在点坝砂体内部非井点控制处侧积层发育的可能性。
[0133]矢量化模型确定模块605，用于根据所述的水平密度将将所述侧积层的空间分布体填充到所述第二填充几何体，得到矢量化的砂体储层构型模型。如图24所示，将侧积层的空间分布体分别填充到点坝砂体中，与点坝砂体表面的交点为b、d。图25即为网格化后的单个点坝砂体构型模拟示意图。[0134]上述步骤利用所定义的基本构型要素地质体进行随机模拟。地下储层构型表征的核心是井间预测，而预测的关键是预知对象的分布规律或模式。虽然模型的调节只是对各个地质体的表征参数进行调节，但是约束依然是各级储层构型定量地质模式。如各个构型要素的展布规模和相互之间的切割和叠置关系。国内外进行过不少的储层地质知识库的研究，初步建立了垂向信息与侧向规模的关系，涉及了河流相、三角洲相、扇三角洲相、深水浊积相等，取得了丰富的研究成果，这些研究都是进行随机模拟的基础。故而本发明将现在的构型地质定量模式研究作为基础，以构型随机模拟算法的开发为目标，对各类储层构型要素的三维空间随机模拟方法进行研究，最终得到基于矢量表达的储层构型模型。
[0135]在定义好各个构型要素之后则是进行基于目标的随机模拟。与传统的基于目标体的模拟方法一样，认为目标中心点位置符合齐次泊松点过程，而在无井区，模拟采用非齐次的泊松点过程，从而满足了井间与井点分布具有差异的要求。模拟仍然采用迭代的算法，并且在模拟过程的必须受到各个目标体的叠置切割模式和井数据的约束。
[0136]地质构型分析所得到的单井构型解释是作为模拟阶段的硬数据，最终的构型模型必须满足井数据的约束。本发明在模拟程序中采用了优先拟合井点数据的做法。首先通过调整各个地质体的形态参数使其满足井点约束然。之后在选定的地质模式约束下进行空间其它构型要素的投放以满足目标函数。在点坝砂体内部模拟侧积层为例(如图24)，依据井数据模拟出了侧积层a和C。而依据侧积层间距等定量构型模式信息可以模拟出侧积层b和d。
[0137]由图7可知，该系统还包括:
[0138]网格化装置700，用于对所述的矢量化的砂体储层构型模型进行网格化，得到网格化的砂体储层构型模型。图11为网格化装置700的结构框图，由图11可知，网格化装置包括:
[0139]网格尺寸确定模块701，用于根据所述的侧积层数据确定网格尺寸。在具体的实施方式中，如侧积层数据为垂向上最厚处为0.8米，则确定的网格尺寸为0.5米。
[0140]网格建立模块702，用于根据所述的网格尺寸对所述矢量化的砂体储层构型模型建立模拟区的网格；
[0141]数值填充模块703，用于对所述模拟区的网格根据构型要素类型进行数值填充，得到网格化的砂体储层构型模型。图12为数值填充模块703的结构框图，由图12可知，数值填充模块具体包括:
[0142]第一判断单元7031，用于判断所述模拟区的网格对应的构型要素类型是否仅为废弃河道；
[0143]第一填充单元7032，用于当所述的第一判断单元判断为是时，将所述的网格填充为废弃河道数据；
[0144]第二判断单元7033，用于当所述的第一判断单元判断为否时，判断所述模拟区的网格对应的构型要素类型是否仅为侧积层；
[0145]第二填充单元7034，用于当所述的第二判断单元判断为是时，将所述的网格填充为侧积层数据；
[0146]第三判断单元7035，用于当所述的第二判断单元判断为否时，判断所述模拟区的网格对应的构型要素类型是否仅为点坝砂体；[0147]第三填充单元7036，用于当所述的第三判断单元判断为是时，将所述的网格填充为点坝砂体数据；
[0148]比例确定单元7037，用于当所述第三判断单元判断为否时，确定所述网格对应的构型要素类型的比例。当网格中包括多种构型要素类型时，分别确定各个要素所占的比例，如:网格中既有点坝砂体，又有侧积层，则确定二者的比例，为2:5。
[0149]最大值获取单元7038，用于获取所述比例中的最大值。
[0150]要素类型确定单元7039，用于确定所述最大值对应的构型要素类型；
[0151]第四填充单元7040，用于将所述的网格填充为最大值对应的构型要素类型对应的数据。
[0152]如上的实施例中，网格中既有点坝砂体，又有侧积层，二者的比例为2:5，则比例最大之为5，对应的构型要素类型为侧积层，因此该网格填充的数据为侧积层数据。如图25为单个点坝砂体构型模型网格化结果。
[0153]由于后期物性参数的建模方法仍然是基于网格的，而且后期的数值模拟也是基于离散网格数据体，故而模拟的最终结果应该是一个基于网格的模型。所以必须将矢量存储的空间地质体进行离散化，将各个储层构型属性赋于每一个网格。由于基于矢量存储的地质模型不受空间网格大小的限制，并且不再考虑模拟速度的要求，此阶段的网格定义可以采取任意细化以满足刻画不同规模的非均质地质体。故而如何对基于矢量表达的模型进行网格化处理也是本发明的内容之一。由于此阶段不早模拟，没有迭代等计算，所以可以选择能够刻画最小规模的构型要素的模型网格进行离散化，如图25为单个点坝砂体构型模型网格化结果。
[0154]下面结合具体的实施例，详细介绍本发明的技术方案。研究区井距最大约380m，最小约25m，平均井距IOOm左右。研究区目的层位为馆陶组上段，其细分为5个小层，其中本文选取S3小层为例。S3小层发育点坝、天然堤、决口扇等微相。岩性自下而上为中、细粒度砂岩渐变为细、粉砂岩至泥岩具典型正韵律特征，底部有冲刷面，河流二元结构特征明显。
[0155](I)读取工区85 口井的自然电位、自然伽马、微电极曲线等测井曲线数据、分层数据。
[0156](2)依据测井曲线特征，确定S3小层各井的储层构型类型。如曲流河河道砂体中的泥质侧积层的岩性主要包括泥岩、页岩、粉砂质泥岩及部分泥质粉砂岩。岩电标定结果显示研究区内泥质侧积层较薄，一般为0.2?0.Sm，测井曲线表现出微电极曲线明显回返，幅度差减小，自然伽马曲线见回返，自然电位曲线轻微回返或者不明显的回返。过工区的一条沉积相剖面连井对比分析。通过以上地质分析所解释的单井各级构型要素将作为构型模拟的条件数据。
[0157](3)对于本工区，根据经验公式和井数据，推断单河道河流满岸深度为7.2m，平均河流满岸宽度162m，单一侧积体水平宽度约为106m左右，侧积层的倾角约为5?10°。
[0158](4)得到本工区中使用本发明得到的单个点坝砂体构型模拟示意图，模型较好的再现了点坝内部及点坝内部侧积层的空间分布特征。
[0159]综上所述，本发明提供了一种基于空间矢量的砂体储层构型建模的方法及系统，采用的基于空间矢量的储层构型建模方法可以有效解决现有技术存在的问题。该方案提出的主要目的是在储层构型建模中描述不同级别的非均质体以及准确刻画各级地质体的复杂空间形态。在模拟阶段的地质体的存储借鉴了矢量图形存储的思路，即对模拟地质体的定义采用了基于实数集的点、线、面所限定的体。因而模型在模拟过程中不受网格尺寸的限制，使得模拟结果可以刻画不同规模层次的非均质体。应用基于矢量体存储方法的另外一个优点是可以快速的对模型进行调节，因为在模拟过程中只是对地质体的空间位置和形态的定义参数进行调节，因而这样的调节方式较基于网格的调节更便捷，从而使结果既忠实条件数据点又易再现模拟目标地质体的几何形态。
[0160]该方案依据地质家对研究工区的地质模式定义模型中的构型要素。用构型要素形态参数和数学关系来描述其空间形态及彼此之间的相对位置关系。
[0161]该方案利用所定义的构型要素进行随机模拟。虽然仍然采用了传统的基于目标体的模拟，但是并不设定模型的网格尺寸，而是直接在模型内用面和体的方法生成构型模型。模型只存储构型要素空间分布形态参数。并且模型的条件化过程仅对构型要素的形态表征参数和位置关系参数进行调节，这样避免了基于网格的方法中网格个数增加带来的计算时间挑战。
[0162]该方案对模拟的结果进行网格化。由于不再考虑模拟速度的要求，此阶段可以采取任意细化的网格以满足刻画所关注的关键构型要素，如点坝砂体中的侧积层。
[0163]本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一般计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory, ROM)或随机存储记忆体(Random AccessMemory, RAM)等。
[0164]本领域技术人员还可以了解到本发明实施例列出的各种功能是通过硬件还是软件来实现取决于特定的应用和整个系统的设计要求。本领域技术人员可以对于每种特定的应用，可以使用各种方法实现所述的功能，但这种实现不应被理解为超出本发明实施例保护的范围。
[0165]本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在【具体实施方式】及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。
【权利要求】
1.一种基于空间矢量的砂体储层构型建模的方法，其特征是，所述的方法具体包括: 获取当前砂体的测井数据； 对所述的测井数据进行数字处理，得到侧积层数据、废弃河道数据以及点坝砂体数据； 对所述的废弃河道数据进行拟合，得到平面分布的深泓线； 根据所述的深泓线以及点坝砂体数据建立废弃河道三维几何体； 根据反曲曲线公式、侧积层数据确定侧积层的空间分布体； 根据所述废弃河道三维几何体、侧积层的空间分布体以及点坝砂体数据确定矢量化的砂体储层构型模型； 对所述的矢量化的砂体储层构型模型进行网格化，得到网格化的砂体储层构型模型。
2.根据权利要求1所述的方法，其特征是，根据所述的深泓线以及点坝砂体数据建立废弃河道三维几何体具体包括: 从所述的点坝砂体数据中获取砂体厚度； 将所述深泓线上的任一点作为数据点； 根据所述的数据点以及曲河流点坝砂体经验公式确定侧积层水平宽度以及河道宽 度； 根据所述的数据点、侧积层水平宽度、河道宽度以及砂体厚度确定废弃河道的剖面； 遍历所述深泓线上的点，得到多个废弃河道的剖面； 多个废弃河道的剖面之间圈定的空间体积即为废弃河道的三维几何体。
3.根据权利要求2所述的方法，其特征是，根据反曲曲线公式、侧积层数据确定侧积层的空间分布体具体包括: 根据反曲曲线公式确定在实数范围内的反曲曲线形态； 对所述的反曲曲线形态进行坐标平移或坐标旋转，得到侧积层剖面； 从所述的侧积层数据中获取侧积层中心厚度、侧积层水平跨度； 设置两个圆心距为侧积层中心厚度、直径为侧积层水平跨度的圆； 将所述侧积层剖面对应的两个反曲曲线沿两个圆相交的圆弧线进行平行移动； 平行移动侧积层剖面所圈定的空间即为侧积层的空间分布体。
4.根据权利要求3所述的方法，其特征是，根据所述废弃河道三维几何体、侧积层的空间分布体以及点坝砂体数据确定矢量化的砂体储层构型模型网格具体包括: 将所述的点坝砂体数据填充到所述废弃河道三维几何体所圈定的区域中，得到第一填充三维几何体； 从所述的点坝砂体数据中获取侧积层在井身上的位置、侧积层的发育角度以及侧积层的水平密度； 根据所述侧积层在井身上的位置、侧积层的发育角度以及侧积层的水平密度，将所述侧积层的空间分布体填充到第一填充三维几何体，得到第二填充几何体； 从所述的点坝砂体数据中获取侧积层的水平密度； 根据所述的水平密度将将所述侧积层的空间分布体填充到所述第二填充几何体，得到矢量化的砂体储层构型模型。
5.根据权利要求4所述的方法，其特征是，对所述的矢量化的砂体储层构型模型进行网格化，得到网格化的砂体储层构型模型具体包括: 根据所述的侧积层数据确定网格尺寸； 根据所述的网格尺寸对所述矢量化的砂体储层构型模型建立模拟区的网格； 对所述模拟区的网格根据构型要素类型进行数值填充，得到网格化的砂体储层构型模型。
6.根据权利要求5所述的方法，其特征是，对所述模拟区的网格根据构型要素类型进行数值填充具体包括: 判断所述模拟区的网格对应的构型要素类型是否仅为废弃河道； 当判断为是时，将所述的网格填充为废弃河道数据； 否则，判断所述模拟区的网格对应的构型要素类型是否仅为侧积层； 当判断为是时，将所述的网格填充为侧积层数据； 否则，判断所述模拟区的网格对应的构型要素类型是否仅为点坝砂体； 当判断为是时 ，将所述的网格填充为点坝砂体数据； 否则，确定所述网格对应的构型要素类型的比例； 获取所述比例中的最大值； 确定所述最大值对应的构型要素类型； 将所述的网格填充为最大值对应的构型要素类型对应的数据。
7.一种基于空间矢量的砂体储层构型建模的系统，其特征是，所述的系统具体包括: 测井数据获取装置，用于获取当前砂体的测井数据； 测井数据处理装置，用于对所述的测井数据进行数字处理，得到侧积层数据、废弃河道数据以及点坝砂体数据； 拟合装置，用于对所述的废弃河道数据进行拟合，得到平面分布的深泓线； 三维几何体建立装置，用于根据所述的深泓线以及点坝砂体数据建立废弃河道三维几何体； 空间分布体建立装置，用于根据反曲曲线公式、侧积层数据确定侧积层的空间分布体； 矢量化模型确定装置，用于根据所述废弃河道三维几何体、侧积层的空间分布体以及点坝砂体数据确定矢量化的砂体储层构型模型； 网格化装置，用于对所述的矢量化的砂体储层构型模型进行网格化，得到网格化的砂体储层构型模型。
8.根据权利要求7所述的系统，其特征是，所述的三维几何体建立装置具体包括: 砂体厚度获取模块，用于从所述的点坝砂体数据中获取砂体厚度； 数据点选取模块，用于将所述深泓线上的任一点作为数据点； 河道宽度确定模块，用于根据所述的数据点以及曲河流点坝砂体经验公式确定侧积层水平宽度、河道宽度； 废弃河道剖面确定模块，用于根据所述的数据点、侧积层水平宽度、河道宽度以及砂体厚度确定废弃河道的剖面； 数据点遍历模块，用于遍历所述深泓线上的点，得到多个废弃河道的剖面； 三维几何体确定模块，多个废弃河道的剖面之间圈定的空间体积即为废弃河道的三维几何体。
9.根据权利要求8所述的系统，其特征是，所述的空间分布体建立装置具体包括: 反曲曲线形态确定模块，用于根据反曲曲线公式确定在实数范围内的反曲曲线形态； 侧积层剖面确定模块，用于对所述的反曲曲线形态进行坐标平移或坐标旋转，得到侧积层剂面； 水平跨度获取模块，用于从所述的侧积层数据中获取侧积层中心厚度、侧积层水平跨度； 设置模块，用于设置两个圆心距为侧积层中心厚度、直径为侧积层水平跨度的圆；平行移动模块，用于将所述侧积层剖面对应的两个反曲曲线沿两个圆相交的圆弧线进行平行移动； 空间分布体确定模块，平行移动侧积层剖面所圈定的空间即为侧积层的空间分布体。
10.根据权利要求9所述的系统，其特征是，所述的矢量化模型确定装置具体包括: 第一填充模块，用于将所述的点坝砂体数据填充到所述废弃河道三维几何体所圈定的区域中，得到第一填充三维几何体； 角度获取模块，用于从所述的点坝砂体数据中获取侧积层在井身上的位置、侧积层的发育角度以及侧积层的水平密度； 第二填充模块，用于根据所述侧积层在井`身上的位置、侧积层的发育角度以及侧积层的水平密度，将所述侧积层的空间分布体填充到第一填充三维几何体，得到第二填充几何体； 水平密度获取模块，用于从所述的点坝砂体数据中获取侧积层的水平密度； 矢量化模型确定模块，用于根据所述的水平密度将将所述侧积层的空间分布体填充到所述第二填充几何体，得到矢量化的砂体储层构型模型。
11.根据权利要求10所述的系统，其特征是，所述的网格化装置具体包括: 网格尺寸确定模块，用于根据所述的侧积层数据确定网格尺寸； 网格建立模块，用于根据所述的网格尺寸对所述矢量化的砂体储层构型模型建立模拟区的网格； 数值填充模块，用于对所述模拟区的网格根据构型要素类型进行数值填充，得到网格化的砂体储层构型模型。
12.根据权利要求11所述的系统，其特征是，所述的数值填充模块具体包括: 第一判断单元，用于判断所述模拟区的网格对应的构型要素类型是否仅为废弃河道； 第一填充单元，用于当所述的第一判断单元判断为是时，将所述的网格填充为废弃河道数据； 第二判断单元，用于当所述的第一判断单元判断为否时，判断所述模拟区的网格对应的构型要素类型是否仅为侧积层； 第二填充单元，用于当所述的第二判断单元判断为是时，将所述的网格填充为侧积层数据； 第三判断单元，用于当所述的第二判断单元判断为否时，判断所述模拟区的网格对应的构型要素类型是否仅为点坝砂体； 第三填充单元，用于当所述的第三判断单元判断为是时，将所述的网格填充为点坝砂体数据； 比例确定单元，用于当所述第三判断单元判断为否时，确定所述网格对应的构型要素类型的比例； 最大值获取单元，用于获取所述比例中的最大值； 要素类型确定单元，用于确定所述最大值对应的构型要素类型； 第四填充单元，用于将 所述的网格填充为最大值对应的构型要素类型对应的数据。
【文档编号】G06F17/50GK103699751SQ201310744975
【公开日】2014年4月2日 申请日期:2013年12月30日 优先权日:2013年12月30日
【发明者】李宇鹏, 吴胜和 申请人:中国石油大学(北京)