CN104504754B - 一种油气储层多点统计建模的方法及装置 - Google Patents

Abstract

本发明提供了一种多点统计建模的方法及装置，涉及油气储层地质建模技术领域，该方法包括：根据预先设置的训练图像，建立搜索树；生成多点统计模型的三维网格系统；在三维网格系统中，随机选取一个节点作为目标节点；所述目标节点和一个具有原始数据的节点或者和一个已生成模拟结果的节点相邻接；在目标节点不具有模拟结果且不具有原始数据时，以目标节点为中心，以预先设置的搜索模板尺寸为半径，在一椭圆邻域内查找具有条件数据的节点；在未查找到具有条件数据的节点时，在三维网格系统中继续随机选取一个节点作为目标节点。本发明能够解决Snesim算法在选取各个节点时，选择是没有任何规则地，当仅利用测井数据建模时，容易造成模拟结果出现不连续性的问题。

Description

一种油气储层多点统计建模的方法及装置

技术领域

本发明涉及油气储层地质建模技术领域，尤其涉及一种油气储层多点统计建模的方法及装置。

背景技术

当前，在众多储层地质建模方法中，多点统计建模(Multiple Point Simulation，简称MPS)成为储层地质建模的主流核心算法并被广泛使用。多点统计建模在理论上和应用上，极大地推动了储层地质建模的发展。多点统计建模以地质统计学为基础，利用蒙得卡洛方法来完成空间沉积微相信息的模拟。在建模时，为了模拟空间每个被模拟点的信息，需要两方面的输入信息：①被模拟点周围的各个已知点的位置和它们相应的信息；②被模拟点和周围各个已知点之间的先验的空间结构模型。多点统计建模的应用相当广泛，特别适用于河流相沉积的储层的建模。

多点统计建模方法中以Snesim(Single Normal Equation Simulation)算法的出现最早、应用最为广泛。它较好的解决了重复扫描图像的问题，从而使得多点统计建模能利用计算机迅速实施成为现实。在实际应用中，这个算法在忠实于条件数据的条件下能较好的再现储层结构特征。Snesim算法的一般原理是，选取区域中的各个点进行模拟时，在选取各个点时，选择是没有任何规则地，是随机选取的。当仅利用测井数据建模时，因此就容易造成模拟结果出现不连续性。

发明内容

本发明实施例提供一种多点统计建模的方法及装置，以解决现有技术中Snesim算法在选取各个节点时，选择是没有任何规则地，是随机进行的，当仅利用测井数据建模时，容易造成模拟结果出现不连续性的问题。

为了达到上述目的，本发明采用如下技术方案：

一种多点统计建模的方法，包括：

根据预先设置的训练图像，针对所模拟的三维储层区域，建立搜索树；

生成多点统计模型的三维网格系统；所述三维网格系统中包括各节点，其中包括具有原始数据的节点，也包括将被模拟的节点；

在三维网格系统中，随机选取一个节点作为目标节点；所述目标节点和一个具有原始数据的节点或者和一个已生成模拟结果的节点相邻接；

判断目标节点是否具有模拟结果以及是否具有原始数据；

若目标节点不具有模拟结果且不具有原始数据，以目标节点为中心，以预先设置的搜索模板尺寸为半径，在一椭圆邻域内查找具有条件数据的节点，也即其上具有原始数据的节点，或其上有模拟结果的节点；

判断在所述椭圆邻域内是否查找到具有条件数据的节点；

若未查找到具有条件数据的节点，在三维网格系统中，继续随机选取一个节点作为目标节点。

进一步的，该多点统计建模的方法，还包括：

若查找到具有条件数据的节点，在所述搜索树中统计以该目标节点为中心，以邻域存在的数据为条件的数据事件，在所述训练图像中进行扫描和逐步比对，获得相同的数据构型的重复次数；

判断所述重复次数是否大于一预先设置的最小重复次数；

若所述重复次数大于所述预先设置的最小重复次数，通过一搜索模板对训练图像进行扫描获取得到有关沉积微相的局部条件概率，从该局部条件概率的分布中抽样获得一个模拟值，将模拟值作为目标节点处的模拟结果。

进一步的，该多点统计建模的方法，还包括：

若所述重复次数小于或等于所述预先设置的最小重复次数，在已知的搜索模板中去掉离被目标节点距离最远的一个条件数据，形成更新后的搜索模板。

具体的，所述生成多点统计模型的三维网格系统，包括：

根据油气储层的几何形态建立所述油气储层相应的三维网格系统；所述三维网络系统中的各节点均具有三维坐标信息，所述三维坐标信息包括节点的X轴坐标、Y轴坐标和Z轴坐标。

具体的，所述目标节点和一个具有原始数据的节点或者和一个已生成模拟结果的节点相邻接，包括：

所述目标节点的X轴坐标、Y轴坐标和Z轴坐标分别对应于所述具有原始数据的节点的X轴坐标、Y轴坐标和Z轴坐标，且和所述具有原始数据的节点的X轴坐标、Y轴坐标和Z轴坐标之中的任意一个、两个或三个坐标的数值分别相邻；

或者，所述目标节点的X轴坐标、Y轴坐标和Z轴坐标分别对应于所述已生成模拟结果的节点的X轴坐标、Y轴坐标和Z轴坐标，且和所述已生成模拟结果的节点的X轴坐标、Y轴坐标和Z轴坐标中的任意一个、两个或三个坐标的数值分别相邻。

此外，所述原始数据为油气储层若干个节点处的测井数据和地震数据获取得到的微相数据。

一种多点统计建模的装置，包括：

搜索树建立单元，用于根据预先设置的训练图像，针对所模拟的三维储层区域，建立搜索树；

三维网格系统生成单元，用于生成多点统计模型的三维网格系统；所述三维网格系统中包括各节点；

目标节点选取单元，用于在三维网格系统中，随机选取一个节点作为目标节点；所述目标节点和一个具有原始数据的节点或者和一个已生成模拟结果的节点相邻接；

第一判断单元，用于判断目标节点是否具有模拟结果以及是否具有原始数据；

条件数据查找单元，用于在所述目标节点不具有模拟结果且不具有原始数据时，以所述目标节点为中心，以预先设置的搜索模板尺寸为半径，在一椭圆邻域内查找具有条件数据的节点；

第二判断单元，用于判断在所述椭圆邻域内是否查找到具有条件数据的节点；

所述目标节点选取单元，还用于在未查找到具有条件数据的节点时，在三维网格系统中，继续随机选取一个节点作为目标节点。

进一步的，所述多点统计建模的装置，还包括：

重复次数查找单元，用于在查找到具有条件数据的节点时，在所述搜索树中统计以该目标节点为中心，以邻域存在的数据为条件的数据事件，在所述训练图像中进行扫描和逐步比对，获得相同的数据构型的重复次数；

第三判断单元，用于判断所述重复次数是否大于一预先设置的最小重复次数；

局部条件概率获取单元，用于在所述重复次数大于所述预先设置的最小重复次数时，通过一搜索模板对训练图像进行扫描获取得到有关沉积微相的局部条件概率，从该局部条件概率的分布中抽样获得一个模拟值，将模拟值作为目标节点处的模拟结果。

进一步的，所述多点统计建模的装置，还包括：

更新单元，用于在所述重复次数小于或等于预先设置的最小重复次数时，在已知的搜索模板中去掉离被目标节点距离最远的一个条件数据，形成更新后的搜索模板。

此外，所述三维网格系统生成单元，具体用于：

根据油气储层的几何形态建立所述油气储层相应的三维网格系统；所述三维网络系统中的各节点均具有三维坐标信息，所述三维坐标信息包括节点的X轴坐标、Y轴坐标和Z轴坐标。

具体的，所述目标节点选取单元所选取的目标节点和一个具有原始数据的节点或者和一个已生成模拟结果的节点相邻接，包括：

所述目标节点的X轴坐标、Y轴坐标和Z轴坐标分别对应于所述具有原始数据的节点的X轴坐标、Y轴坐标和Z轴坐标，且和所述具有原始数据的节点的X轴坐标、Y轴坐标和Z轴坐标之中的任意一个、两个或三个坐标的数值分别相邻；

或者，所述目标节点的X轴坐标、Y轴坐标和Z轴坐标分别对应于所述已生成模拟结果的节点的X轴坐标、Y轴坐标和Z轴坐标，且和所述已生成模拟结果的节点的X轴坐标、Y轴坐标和Z轴坐标中的任意一个、两个或三个坐标的数值分别相邻。

另外，所述目标节点选取单元中的原始数据为油气储层若干个节点处的测井数据和地震数据获取得到的微相数据。

本发明实施例提供的多点统计建模的方法及装置，能够根据所述搜索树在所述三维网格系统中，随意选取一节点作为目标节点；所述目标节点和一个具有原始数据的节点或者和一个已生成模拟结果的节点相邻接，进而在所述目标节点不具有模拟结果且不具有原始数据时，以所述目标节点为中心，以预先设置的搜索模板尺寸为半径，在一椭圆邻域内查找具有条件数据的节点，并在未查找到具有条件数据的节点时，在三维网格系统中，继续随机选取一个节点作为目标节点，从而得到的模拟结果的连续性相对于现有技术较好，避免了现有技术中Snesim算法在选取各个节点时，选择是没有任何规则地，是随机进行的，当仅利用测井数据建模时，容易造成模拟结果出现不连续性的问题。

附图说明

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

图1为本发明实施例中的多点统计建模的方法的流程图一；

图2为本发明实施例中的多点统计建模的方法的流程图二；

图3为本发明实施例中的三维网络系统的示意图；

图4为本发明实施例中的训练图像和搜索模板的示意图；

图5为本发明实施例中的通过多点统计建模的方法进行模拟的示意图；

图6为本发明实施例中的通过多点统计建模的方法进行模拟和Snesim方法进行模拟的效果图一；

图7为本发明实施例中的通过多点统计建模的方法进行模拟和Snesim方法进行模拟的效果图二；

图8为本发明实施例中的通过多点统计建模的方法进行模拟和Snesim方法进行模拟的效果图三；

图9为本发明实施例中的多点统计建模的装置的结构示意图一；

图10为本发明实施例中的多点统计建模的装置的结构示意图二。

具体实施方式

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

如图1所示，本发明实施例提供一种多点统计建模的方法，包括：

步骤101、根据预先设置的训练图像，针对所模拟的三维储层区域，建立搜索树。

步骤102、生成多点统计模型的三维网格系统。

其中，该三维网格系统中包括各节点，即为各网格块的中心。

步骤103、在三维网格系统中，随机选取一个节点作为目标节点。

其中，该目标节点和一个具有原始数据的节点或者和一个已生成模拟结果的节点相邻接。其中，该目标节点为一模拟值的空间载体。

步骤104、判断目标节点是否具有模拟结果以及是否具有原始数据。

其中，该原始数据为油气储层的测井数据或地震数据获取得到的微相数据。

步骤105、若目标节点不具有模拟结果且不具备原始数据，以目标节点为中心，以预先设置的搜索模板尺寸为半径，在一椭圆邻域内查找具有条件数据的节点。

其中，该条件数据表示其位置处具有原始数据，或有模拟结果。该原始数据是进行模拟的前提条件和输入数据。

步骤106、判断在椭圆邻域内是否查找到具有条件数据的节点。

步骤107、若未查找到具有条件数据的节点，在三维网格系统中，继续随机选取一个节点作为目标节点。.

本发明实施例提供的多点统计建模的方法，能够在三维网格系统中，随机选取一节点作为目标节点；目标节点和一个具有原始数据的节点或者和一个已生成模拟结果的节点相邻接，进而在目标节点不具有模拟结果,同时有不具有原始数据时，以目标节点为中心，以预先设置的搜索模板尺寸为半径，在一椭圆邻域内查找具有条件数据的节点，并在未查找到具有条件数据的节点时能够在在三维网格系统中，继续随机选取一个节点作为目标节点，直至建模完成，从而得到的模拟结果的连续性相对于现有技术较好，避免了现有技术中Snesim算法在选取各个节点时，选择是没有任何规则地，是随机进行的，当仅利用测井数据建模时，容易造成模拟结果出现不连续性的问题。

其中，原始数据为油气储层若干个节点处的测井数据和地震数据获取得到的微相数据。

为了使本领域的技术人员更好的了解本发明，下面列举一个更为详细的实施例，如图2所示，本发明实施例提供的一种多点统计建模的方法，包括：

步骤201、根据预先设置的训练图像，针对所模拟的三维储层区域，建立搜索树。

其中，通过搜索树的数据结构可以加速模拟过程，只要对训练图像进行一次遍历搜索就可以生成该搜索树.将所有可能的概率模式存储在搜索树中，生成模拟图像时可以直接从搜索树上获得该模式的结果数据,从而加快了模拟速度。

训练图象是多点统计模拟算法中的重要的概念。训练图像(Training Image)是地质概念模型的量化表示，地质概念模型体现为在所研究的实际油藏中存在的的地质结构。训练图像可以通过地质研究人员根据他们的知识、和油藏解释，或者由露头数据等来产生。训练图像如同在传统的两点地质统计学中的变差函数一样，表现出沉积微相空间分布的多点模式。多点统计模拟算法从训练图像中捕捉微相的分布模式，并且锚定它们到不同的油藏观测数据(即原始数据)。使用训练图像能够再现复杂的或非线性的地质结构，并且使得模拟所得的几何形态与观测数据相符。

因为训练图像是纯地质的概念模型，其是从油田开发中的条件数据或者其它特殊位置油藏信息中提取出来的。一个训练图像能够反映所研究油藏中微相空间分布的形状和空间结构。

实际中应用的训练图像是定义在三维空间中的。图4给出多点模拟的一个图示，解释了怎样通过扫描训练图像，以便计算微相的局部条件分布。模拟前，需要定义对训练图像进行扫描的一个模板。例如图4中假设被模拟的网格中的位置u是正在被模拟的一个网格。在这个中心为u的模板里(椭圆区域)，有4个数据值(这4个数据值也是条件数据)。其中，两个是砂岩位置(黑色像元)和两个泥岩位置(白色像元)。所有的四个值连同它们的几何结构叫做一个数据事件。接着，这个数据事件用于对描述河道的训练图像进行扫描和逐步比对，以便推断位置u的砂岩概率。在图4的训练图像中发现有4个重复的数据事件，其中模板中心位置u是砂岩的数据事件有三个，是泥岩的数据事件有一个。因此，这个网格u中的为砂岩的概率是0.75，为泥岩的概率是0.25，并利用这个概率分布，可以给象元u一个模拟值。而网格u是砂岩或泥岩均有可能，但是其是砂岩的机会更大一些，因为砂岩比泥岩有更高的概率(0.75>0.25)。利用蒙地卡罗方法对于在网格u处，所获得的微相的概率分布，进行抽样，在这个实例中在被模拟的点处可以确定为砂岩的概率比较大。蒙地卡罗方法具体确定到底是砂岩，还是泥岩。这个模拟值可以加到条件数据组(即是通过模拟获得的微相数值，或是原始数据)，以便在以后的模拟过程中约束其他还未曾模拟的象元。接下来，这个模拟移到针对另一个网格位置。继续进行这个连续过程，直到三维网格系统中的的所有节点都被模拟到，就产生了这个区域的多点统计模拟。

步骤202、根据油气储层的几何形态建立油气储层相应的三维网格系统。

其中，三维网络系统中的各节点均具有三维坐标信息，三维坐标信息包括节点的X轴坐标、Y轴坐标和Z轴坐标。

如图3所示，对于三维网格系统，其坐标系可以为(X,Y,Z：1≤i≤M,1≤j≤N,1≤k≤10)中，对于(X,Y)平面上的每一个节点(i＝i₀,j＝j₀),都有10个节点相对应：(i＝i₀,j＝j₀,1≤k≤10)。i＝i₀,j＝j₀，k从1变化到10，分别代表着三维网格系统从顶界到底界的10片。

如图3所示，一般来说，X方向的网格数可以取为100，Y方向也可以取为100。这里，Z方向的网格数可以取为10，也可以取成20，或30。整个三维网格系统建模的目的就是要将油气储层范围内部的每个网格上模拟得出相应的值，直到将每一个网格处都经过了模拟，即被赋予一个值以后，整个建模才算完成。其中，在建模之前，有的网格块处已经被赋予了相应的值(即输入数据，或称原始数据)，例如，在井轴穿过的那些网格块处就会赋予相应的测井数据。然后，井轴以外的那些节点处则需要通过建模赋予相应的值，这些值就是由建模生成的模拟结果。

在该三维网络系统中，网格数取得越多，模拟的精度就越高，同时计算量也就越大。如果被模拟的储层范围呈现为不规则的三维几何体，那么它的边界就应该呈现为不规则的曲线(或曲面)。图3所示的储层边界比较简单，呈现为一个长方体，但实际的储层边界较为复杂，此处不在储层边界的描述问题上进行赘述。

步骤203、在三维网格系统中，随机选取一个节点作为目标节点。

该目标节点的X轴坐标、Y轴坐标和Z轴坐标分别对应于具有原始数据的节点的X轴坐标、Y轴坐标和Z轴坐标，且和具有原始数据的节点的X轴坐标、Y轴坐标和Z轴坐标之中的任意一个、两个或三个坐标的数值分别相邻。例如，具有原始数据的节点的X轴坐标、Y轴坐标和Z轴坐标为(3，4，5)，则该目标节点的X轴坐标、Y轴坐标和Z轴坐标可以为(3，4，6)，(4，4，5)，(4，5，5)等，但不仅局限于此。

该原始数据为通过油气储层的测井数据和地震数据获取得到的微相数据。

步骤204、判断目标节点是否具有模拟结果，并判断是否具有原始数据。

若目标节点具有模拟结果，执行步骤205；否则，若该目标节点不具有模拟结果，执行步骤207。若该目标节点具有原始数据，则就转到步骤203的入口，并执行步骤203。

步骤205、判断多点统计建模是否已经完成。

若多点统计建模已经完成，则执行步骤206；否则，若该多点统计建模未完成，返回执行步骤203。

步骤206、输出多点统计建模结果，确定多点统计建模完成。

步骤207、以目标节点为中心，以预先设置的搜索模板尺寸为半径，在一椭圆邻域内查找具有条件数据的节点。

其中，该搜索模板是建立在数据空间的三维网格系统上的。该数据空间就是在三维空间中的、具有明确边界的一个区域。在搜索模板内，若干个网格上具有确定的微相数据(即条件数据)。该种微相数据是由测井、地震等手段提供的原始数据，或是前若干步模拟所产生的结果。储层地质建模的目的就是为了利用这些原始的微相数据模拟出这个数据空间内部所有网格上的微相数据。在图4的左上部的图示中，椭圆区域的内部就代表了一个模板。实际上，一个模板是以三维空间中一个椭球面内部的若干个数据点组成，其椭球面的中心是需要进行模拟的点。然而，整个模板包含着两类空间网格的数据，第一类是指在相应椭球内的，原始的输入的测井提供的微相数据。第二类是包含着和前几次模拟结果所得到的，若干个网格内的微相数据。在图4中的椭圆内部有4个数据点，其中两个是代表砂岩，另外两个代表泥岩。这两类数据就是之前的数据事件，即若干个网格的位置，及其上的微相数据。在图4中，椭圆所包围的数据事件则代表着4个网格块的位置，以及网格块上的微相数据(即代表着是砂岩或者泥岩)。

需要说明的是，图4仅是示意性的简化示意图，相应的搜索模板就简化为一个椭圆。而真实的搜索模板是一个三维空间中的椭球，并且其X方向的长轴,，Y方向的短轴的跨度都有几十个网格节点之多。这种三维模板的几何尺寸可以被预先设置。

步骤208、判断在椭圆邻域内是否查找到具有条件数据的节点。

若未查找到具有条件数据的节点，返回执行步骤205；否则，若查找到具有条件数据的节点，执行步骤209。

返回执行步骤205，表示将继续重新选取一个目标节点。

步骤209、在该搜索树中统计以该目标节点为中心，以邻域存在的数据为条件的数据事件，在训练图像中进行扫描和逐步比对，获得相同的数据构型的重复次数。

步骤210、判断重复次数是否大于一预先设置的最小重复次数。

如果判断该重复次数不大于该最小重复次数，执行步骤211；如果判断到重复次数大于该最小重复次数，执行步骤212。

步骤211、在已知的搜索模板中去掉离被目标节点距离最远的一个条件数据，形成更新后的搜索模板。之后返回执行步骤209。

步骤212、通过该搜索模板对训练图像进行扫描获取得到有关沉积微相的局部条件概率，从该局部条件概率的分布中抽样获得一个模拟值，将模拟值作为目标节点处的模拟结果。之后返回执行步骤205。

例如图5所示的X:Y:Z＝35：51：2的区域内，通过上述步骤201-212的算法进行模拟的各点次序。图5左边的示意图表示纵向上的自顶面的第一片，右边的图件表示自顶面的第二片(Z方向)。灰色区域表示该网格块的模拟结果是泥岩，黑色区域表示的模拟结果是砂岩。如图所示，第13次模拟的结果紧靠着第12次模拟结果，第14次模拟的结果紧靠着第13次模拟结果，第15次模拟的结果紧靠着第14次模拟结果，每次模拟的结果均与其上一次模拟的结果相连，可见通过步骤201至步骤212的方式进行的模拟，模拟结果较为连续。在图5左面图中的左上部分，列出了第8个模拟模块的模拟结果。其横坐标为(31 51)，是和其他模拟结果都不紧靠的，所以其模拟的结果不符合步骤201-步骤212的方法。

下面列举一个本发明中多点统计建模的方法(即图6、7、8中的增长算法模拟结果对应的方法)与传统的Snesim算法的结果比较实例，在使用两种方法模拟时，X、Y、Z方向的三维网格系统的网格数分别为X、Y、Z轴方向100(网格)、100(网格)、10(网格)。训练图像X、Y、Z方向的网格数为108(网格)、160(网格)、10(网格)。硬数据为工区内的测井数据(158口井)。训练图像内共3种沉积相:泥岩、河道、心滩，它们的含量分别为：0.74、0.22、0.04。搜索模板的尺寸设置为：15(网格)*15(网格)*3(网格)。

以下的图6、图7、图8，表示了Snesim方法和增长方法的模拟结果和相应的训练图象(分别为Z方向的第2片，第9片，第10片的结果)。对比各片两种方法的模拟结果，本发明中多点统计建模的方法所得的模拟结果所呈现出的河道在整个研究区的连续性良好，而Snesim方法的模拟结果中河道的连续性会有时会出现明显的中断。这说明了增长方法能更好的展现储层形态和分布规律，更适合于储层随机建模。

本发明实施例提供的多点统计建模的方法，能够在所述三维网格系统中，随意选取一节点作为目标节点；所述目标节点和一个具有原始数据的节点或者和一个已生成模拟结果的节点相邻接，进而在所述目标节点不具有模拟结果且不具有原始数据时，以所述目标节点为中心，以预先设置的搜索模板尺寸为半径，在一椭圆邻域内查找具有条件数据的节点，并在未查找到具有条件数据的节点时，在三维网格系统中，继续随机选取一个节点作为目标节点，从而得到的模拟结果的连续性相对于现有技术较好，避免了现有技术中Snesim算法在选取各个节点时，选择是没有任何规则地，是随机进行的，当仅利用测井数据建模时，容易造成模拟结果出现不连续性的问题。

对应于上述图1和图2的方法实施例，本发明还提供一种多点统计建模的装置，如图9所示，包括：

搜索树建立单元31，可以根据预先设置的训练图像，针对所模拟的三维储层区域，建立搜索树。

三维网格系统生成单元32，可以生成多点统计模型的三维网格系统。

其中，该三维网格系统中包括各节点。

目标节点选取单元33，可以在三维网格系统中，随机选取一个节点作为目标节点。

其中，该目标节点和一个具有原始数据的节点或者和一个已生成模拟结果的节点相邻接。

第一判断单元34，可以判断目标节点是否具有模拟结果以及是否具有原始数据。

条件数据查找单元35，可以在目标节点不具有模拟结果且不具有原始数据时，以目标节点为中心，以预先设置的搜索模板尺寸为半径，在一椭圆邻域内查找具有条件数据的节点。

第二判断单元36，可以判断在椭圆邻域内是否查找到具有条件数据的节点。

所述目标节点选取单元33，还用于在未查找到具有条件数据的节点时，在三维网格系统中，继续随机选取一个节点作为目标节点。

进一步的，如图10所示，多点统计建模的装置，还包括：

重复次数查找单元38，可以在查找到具有条件数据的节点时，在搜索树中统计以该目标节点为中心，以邻域存在的数据为条件的数据事件在训练图像中进行扫描和逐步比对，获得相同的数据构型的重复次数。

第三判断单元39，可以判断重复次数是否大于一预先设置的最小重复次数。

局部条件概率获取单元40，可以在重复次数大于预先设置的最小重复次数时，通过一搜索模板对训练图像进行扫描获取得到有关沉积微相的局部条件概率，从该局部条件概率的分布中抽样获得一个模拟值，将模拟值作为目标节点处的模拟结果。

进一步的，如图10所示，多点统计建模的装置，还包括：

更新单元41，可以在数据事件的重复次数小于或等于预先设置的最小重复次数时，可以去掉离被目标节点距离最远的一个条件数据，形成更新后的搜索模板。

此外，该三维网格系统生成单元32，具体可以根据油气储层的几何形态建立油气储层相应的三维网格系统；三维网络系统中的各节点均具有三维坐标信息，三维坐标信息包括节点的X轴坐标、Y轴坐标和Z轴坐标。

具体的，该目标节点选取单元33所选取的目标节点和一个具有原始数据的节点或者和一个已生成模拟结果的节点相邻接，包括：

该目标节点的X轴坐标、Y轴坐标和Z轴坐标，与具有原始数据的节点的X轴坐标、Y轴坐标和Z轴坐标之中的任意一个、两个或三个坐标的数值分别相邻。

或者，目标节点的X轴坐标、Y轴坐标和Z轴，与已生成模拟结果的节点的X轴坐标、Y轴坐标和Z轴坐标中的任意一个、两个或三个坐标的数值分别相邻。

另外，目标节点选取单元33中的原始数据为通过油气储层的测井数据和地震数据获取得到的微相数据。

值得说明的是，本发明实施例提供的多点统计建模的装置的具体实现方式可以参见图1、图2的方法实施例，此处不再赘述。

本发明实施例提供的多点统计建模的装置，能够根据所述搜索树在所述三维网格系统中，随意选取一节点作为目标节点；所述目标节点和一个具有原始数据的节点或者和一个已生成模拟结果的节点相邻接，进而在所述目标节点不具有模拟结果且不具有原始数据时，以所述目标节点为中心，以预先设置的搜索模板尺寸为半径，在一椭圆邻域内查找具有条件数据的节点，并在未查找到具有条件数据的节点时，在三维网格系统中，继续随机选取一个节点作为目标节点，从而得到的模拟结果的连续性相对于现有技术较好，避免了现有技术中Snesim算法在选取各个节点时，选择是没有任何规则地，是随机进行的，当仅利用测井数据建模时，容易造成模拟结果出现不连续性的问题。

本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。

Claims (12)

1.一种油气储层多点统计建模的方法，其特征在于，包括：

根据预先设置的训练图像，针对所模拟的三维储层区域，建立搜索树；

生成多点统计模型的三维网格系统；所述三维网格系统中包括各节点；

在三维网格系统中，随机选取一个节点作为目标节点；所述目标节点和一个具有原始数据的节点或者和一个已生成模拟结果的节点相邻接；

判断目标节点是否具有模拟结果以及是否具有原始数据；

若目标节点不具有模拟结果且不具有原始数据，以目标节点为中心，以预先设置的搜索模板尺寸为半径，在一椭圆邻域内查找具有条件数据的节点；

判断在所述椭圆邻域内是否查找到具有条件数据的节点；

若未查找到具有条件数据的节点，在三维网格系统中，继续随机选取一个节点作为目标节点。

2.根据权利要求1所述的油气储层多点统计建模的方法，其特征在于，还包括：

若查找到具有条件数据的节点，在所述搜索树中统计以该目标节点为中心，以邻域存在的数据为条件的数据事件，在所述训练图像中进行扫描和逐步比对，获得相同的数据构型的重复次数；

判断所述重复次数是否大于一预先设置的最小重复次数；

若所述重复次数大于所述预先设置的最小重复次数，通过一搜索模板对训练图像进行扫描获取得到有关沉积微相的局部条件概率，从该局部条件概率的分布中抽样获得一个模拟值，将模拟值作为目标节点处的模拟结果。

3.根据权利要求2所述的油气储层多点统计建模的方法，其特征在于，还包括：

若所述重复次数小于或等于预先设置的最小重复次数，在已知的搜索模板中去掉离被目标节点距离最远的一个条件数据，形成更新后的搜索模板。

4.根据权利要求3所述的油气储层多点统计建模的方法，其特征在于，所述生成多点统计模型的三维网格系统，包括：

根据油气储层的几何形态建立所述油气储层相应的三维网格系统；所述三维网络系统中的各节点均具有三维坐标信息，所述三维坐标信息包括节点的X轴坐标、Y轴坐标和Z轴坐标。

5.根据权利要求4所述的油气储层多点统计建模的方法，其特征在于，所述目标节点和一个具有原始数据的节点或者和一个已生成模拟结果的节点相邻接，包括：

所述目标节点的X轴坐标、Y轴坐标和Z轴坐标分别对应于所述具有原始数据的节点的X轴坐标、Y轴坐标和Z轴坐标，且和所述具有原始数据的节点的X轴坐标、Y轴坐标和Z轴坐标之中的任意一个、两个或三个坐标的数值分别相邻；

或者，所述目标节点的X轴坐标、Y轴坐标和Z轴坐标分别对应于所述已生成模拟结果的节点的X轴坐标、Y轴坐标和Z轴坐标，且和所述已生成模拟结果的节点的X轴坐标、Y轴坐标和Z轴坐标中的任意一个、两个或三个坐标的数值分别相邻。

6.根据权利要求5所述的油气储层多点统计建模的方法，其特征在于，所述原始数据为油气储层若干个节点处的测井数据和地震数据获取得到的微相数据。

7.一种油气储层多点统计建模的装置，其特征在于，包括：

搜索树建立单元，用于根据预先设置的训练图像，针对所模拟的三维储层区域，建立搜索树；

三维网格系统生成单元，用于生成多点统计模型的三维网格系统；所述三维网格系统中包括各节点；

目标节点选取单元，用于在三维网格系统中，随机选取一个节点作为目标节点；所述目标节点和一个具有原始数据的节点或者和一个已生成模拟结果的节点相邻接；

第一判断单元，用于判断目标节点是否具有模拟结果以及是否具有原始数据；

条件数据查找单元，用于在所述目标节点不具有模拟结果且不具有原始数据时，以所述目标节点为中心，以预先设置的搜索模板尺寸为半径，在一椭圆邻域内查找具有条件数据的节点；

第二判断单元，用于判断在所述椭圆邻域内是否查找到具有条件数据的节点；

所述目标节点选取单元，还用于在未查找到具有条件数据的节点时，在三维网格系统中，继续随机选取一个节点作为目标节点。

8.根据权利要求7所述的油气储层多点统计建模的装置，其特征在于，还包括：

重复次数查找单元，用于在查找到具有条件数据的节点时，在所述搜索树中统计以该目标节点为中心，以邻域存在的数据为条件的数据事件，在所述训练图像中进行扫描和逐步比对，获得相同的数据构型的重复次数；

第三判断单元，用于判断所述重复次数是否大于一预先设置的最小重复次数；

局部条件概率获取单元，用于在所述重复次数大于所述预先设置的最小重复次数时，通过一搜索模板对训练图像进行扫描获取得到有关沉积微相的局部条件概率，从该局部条件概率的分布中抽样获得一个模拟值，将模拟值作为目标节点处的模拟结果。

9.根据权利要求8所述的油气储层多点统计建模的装置，其特征在于，还包括：

更新单元，用于在所述重复次数小于或等于预先设置的最小重复次数时，在已知的搜索模板中去掉离被目标节点距离最远的一个条件数据，形成更新后的搜索模板。

10.根据权利要求9所述的油气储层多点统计建模的装置，其特征在于，所述三维网格系统生成单元，具体用于：

根据油气储层的几何形态建立所述油气储层相应的三维网格系统；所述三维网络系统中的各节点均具有三维坐标信息，所述三维坐标信息包括节点的X轴坐标、Y轴坐标和Z轴坐标。

11.根据权利要求10所述的油气储层多点统计建模的装置，其特征在于，所述目标节点选取单元所选取的目标节点和一个具有原始数据的节点或者和一个已生成模拟结果的节点相邻接，包括：

所述目标节点的X轴坐标、Y轴坐标和Z轴坐标分别对应于所述具有原始数据的节点的X轴坐标、Y轴坐标和Z轴坐标，且和所述具有原始数据的节点的X轴坐标、Y轴坐标和Z轴坐标之中的任意一个、两个或三个坐标的数值分别相邻；

或者，所述目标节点的X轴坐标、Y轴坐标和Z轴坐标分别对应于所述已生成模拟结果的节点的X轴坐标、Y轴坐标和Z轴坐标，且和所述已生成模拟结果的节点的X轴坐标、Y轴坐标和Z轴坐标中的任意一个、两个或三个坐标的数值分别相邻。

12.根据权利要求11所述的油气储层多点统计建模的装置，其特征在于，所述目标节点选取单元中的原始数据为油气储层若干个节点处的测井数据和地震数据获取得到的微相数据。