一种岩相概率分布模型的建立方法及系统

摘要

本申请提供一种岩相概率分布模型的建立方法及系统，其中，所述方法包括：载入已建立的三维构造模型；在所述三维构造模型中导入预设岩相数据，并对所述预设岩相数据进行粗化处理，得到粗化岩相数据；将预设波阻抗数据体采样至所述三维构造模型中，得到所述三维构造模型对应的波阻抗数据；根据所述波阻抗数据以及所述粗化岩相数据，确定各个波阻抗值对应的指定岩相数据的概率值，以得到所述指定岩相数据的概率分布数据；分别载入所述粗化岩相数据中各个岩相数据的概率分布数据，以在同一个直方图窗口中展示所述粗化岩相数据的概率分布模型。本申请提供的技术方案，能够结合波阻抗和岩相数据来建立岩相概率分布模型。

权利要求

1.一种岩相概率分布模型的建立方法，其特征在于，所述方法包括：

载入已建立的三维构造模型；

在所述三维构造模型中导入预设岩相数据，并对所述预设岩相数据进行粗化处理，得到粗化岩相数据；所述粗化岩相数据中包括三种岩相数据；

将预设波阻抗数据体采样至所述三维构造模型中，得到所述三维构造模型对应的波阻抗数据；

根据所述波阻抗数据以及所述粗化岩相数据，确定各个波阻抗值对应的指定岩相数据的概率值，以得到所述指定岩相数据的概率分布数据；

分别载入所述粗化岩相数据中各个岩相数据的概率分布数据，以在同一个直方图窗口中展示所述粗化岩相数据的概率分布模型；

其中，确定各个波阻抗值对应的指定岩相数据的概率值，包括：构建所述粗化岩相数据中指定岩相数据所对应的波阻抗数据的直方图；基于所述直方图，生成所述指定岩相数据对应的波阻抗分布函数；依次从所述波阻抗分布函数中提取所述粗化岩相数据中三种岩相数据的波阻抗数据对，所述波阻抗数据对中包括波阻抗值和对应的网格数量；根据提取的所述波阻抗数据对，计算所述指定岩相数据所对应的各个波阻抗值的网格数量在指定波阻抗值总网格数量中的比值，并将所述比值作为各个波阻抗值对应的概率值，得到各种岩相数据的概率分布数据。

2.根据权利要求1所述的方法，其特征在于，在载入已建立的三维构造模型之后，所述方法还包括：将所述三维构造模型中的所有层设置为相同的网格厚度。

3.根据权利要求1所述的方法，其特征在于，在构建所述粗化岩相数据中指定岩相数据所对应的波阻抗数据的直方图时，所述方法还包括：为所述直方图设置数值滤波器，并将所述数值滤波器应用于所述指定岩相数据中。

4.根据权利要求1所述的方法，其特征在于，所述方法还包括：

当所述指定岩相数据中的波阻抗值小于或者等于第一阈值时，将所述波阻抗值对应的第一岩相的概率值设置为1，将所述波阻抗值对应的第二岩相的概率值设置为0，将所述波阻抗值对应的第三岩相的概率值设置为0；

当所述指定岩相数据中的波阻抗值大于或者等于第二阈值时，将所述波阻抗值对应的第三岩相的概率值设置为1将所述波阻抗值对应的第一岩相的概率值设置为0，将所述波阻抗值对应的第二岩相的概率值设置为0。

5.一种岩相概率分布模型的建立系统，其特征在于，所述系统包括：

模型载入单元，用于载入已建立的三维构造模型；

粗化处理单元，用于在所述三维构造模型中导入预设岩相数据，并对所述预设岩相数据进行粗化处理，得到粗化岩相数据；所述粗化岩相数据中包括三种岩相数据；

波阻抗数据采样单元，用于将预设波阻抗数据体采样至所述三维构造模型中，得到所述三维构造模型对应的波阻抗数据；

概率分布数据确定单元，用于根据所述波阻抗数据以及所述粗化岩相数据，确定各个波阻抗值对应的指定岩相数据的概率值，以得到所述指定岩相数据的概率分布数据；

概率分布模型展示单元，用于分别载入所述粗化岩相数据中各个岩相数据的概率分布数据，以在同一个直方图窗口中展示所述粗化岩相数据的概率分布模型；

其中，所述概率分布数据确定单元，包括：直方图构建模块，用于构建所述粗化岩相数据中指定岩相数据所对应的波阻抗数据的直方图；基于所述直方图，生成所述指定岩相数据对应的波阻抗分布函数；数据对提取模块，用于依次从所述波阻抗分布函数中提取所述粗化岩相数据中三种岩相数据的波阻抗数据对，所述波阻抗数据对中包括波阻抗值和对应的网格数量；概率值计算模块，用于根据提取的所述波阻抗数据对，计算所述指定岩相数据所对应的各个波阻抗值的网格数量在指定波阻抗值总网格数量中的比值，并将所述比值作为各个波阻抗值对应的概率值，得到各种岩相数据的概率分布数据。

6.根据权利要求5所述的系统，其特征在于，所述系统还包括：

第一设置模块，用于当波阻抗值小于或者等于第一阈值时，将所述波阻抗值对应的指定岩相的概率值设置为第一常数；

第二设置模块，用于当波阻抗值大于或者等于第二阈值时，将所述波阻抗值对应的指定岩相的概率值设置为第二常数。

说明书

一种岩相概率分布模型的建立方法及系统

技术领域

[0001] 本申请涉及地质建模技术领域，特别涉及一种岩相概率分布模型的建立方法及系统。

背景技术

[0002] 岩相建模对储集层研究具有重要意义，相模型的合理性直接影响储集层参数预测的精度。地震资料横向分辨率高、横向信息丰富可以降低相模型的不确定性。目前采用的岩相建模方法是利用波阻抗和单井岩相资料建立岩相的概率分布模型，再在这个相概率分布模型的约束下，结合井点资料和地质认识，建立岩相模型。然而目前还没有公开资料介绍利用波阻抗和岩相资料建立岩相概率分布模型的详细方法。

[0003] 应该注意，上面对技术背景的介绍只是为了方便对本申请的技术方案进行清楚、完整的说明，并方便本领域技术人员的理解而阐述的。不能仅仅因为这些方案在本申请的背景技术部分进行了阐述而认为上述技术方案为本领域技术人员所公知。

发明内容

[0004] 本申请实施方式的目的在于提供一种岩相概率分布模型的建立方法及系统，能够结合波阻抗和岩相数据来建立岩相概率分布模型。

[0005] 为实现上述目的，本申请一方面提供一种岩相概率分布模型的建立方法，所述方法包括：载入已建立的三维构造模型；在所述三维构造模型中导入预设岩相数据，并对所述预设岩相数据进行粗化处理，得到粗化岩相数据；所述粗化岩相数据中包括三种岩相数据；将预设波阻抗数据体采样至所述三维构造模型中，得到所述三维构造模型对应的波阻抗数据；根据所述波阻抗数据以及所述粗化岩相数据，确定各个波阻抗值对应的指定岩相数据的概率值，以得到所述指定岩相数据的概率分布数据；分别载入所述粗化岩相数据中各个岩相数据的概率分布数据，以在同一个直方图窗口中展示所述粗化岩相数据的概率分布模型。

[0006] 进一步地，在载入已建立的三维构造模型之后，所述方法还包括：将所述三维构造模型中的所有层设置为相同的网格厚度。

[0007] 进一步地，确定各个波阻抗值对应的指定岩相数据的概率值包括：构建所述粗化岩相数据中指定岩相数据所对应的波阻抗数据的直方图；提取所述直方图中的波阻抗数据对，所述波阻抗数据对中包括波阻抗值和对应的网格数量；根据提取的所述波阻抗数据对，计算各个波阻抗值对应的指定岩相数据的概率值，得到所述指定岩相数据的概率分布数据。

[0008] 进一步地，在构建所述粗化岩相数据中指定岩相数据所对应的波阻抗数据的直方图时，所述方法还包括：为所述直方图设置数值滤波器，并将所述数值滤波器应用于所述指定岩相数据中。

[0009] 进一步地，计算各个波阻抗值对应的指定岩相数据的概率值包括：计算所述指定岩相数据所对应的各个波阻抗值的网格数量在指定波阻抗值总网格数量中的比值，并将所述比值作为各个波阻抗值对应的概率值。

[0010] 进一步地，所述方法还包括：当所述指定岩相数据对应的波阻抗值小于或者等于第一阈值时，将所述波阻抗值对应的所述指定岩相数据的概率值设置为第一常数；当所述指定岩相数据对应的波阻抗值大于或者等于第二阈值时，将所述波阻抗值对应的所述指定岩相数据的概率值设置为第二常数。

[0011] 进一步地，当波阻抗值小于或者等于第一阈值时，将第一岩相的概率值设置为1，将第二岩相和第三岩相中0的概率值设置为0；当波阻抗值大于或者等于第二阈值时，将第一岩相0的概率值设置为0，将第二岩相的概率值设置为0，将第三岩相中的概率值设置为1。

[0012] 为实现上述目的，本申请另一方面还提供一种岩相概率分布模型的建立系统，所述系统包括：模型载入单元，用于载入已建立的三维构造模型；粗化处理单元，用于在所述三维构造模型中导入预设岩相数据，并对所述预设岩相数据进行粗化处理，得到粗化岩相数据；所述粗化岩相数据中包括三种岩相数据；波阻抗数据采样单元，用于将预设波阻抗数据体采样至所述三维构造模型中，得到所述三维构造模型对应的波阻抗数据；概率分布数据确定单元，用于根据所述波阻抗数据以及所述粗化岩相数据，确定各个波阻抗值对应的指定岩相数据的概率值，以得到所述指定岩相数据的概率分布数据；概率分布模型展示单元，用于分别载入所述粗化岩相数据中各个岩相数据的概率分布数据，以在同一个直方图窗口中展示所述粗化岩相数据的概率分布模型。

[0013] 进一步地，所述概率分布数据确定单元包括：直方图构建模块，用于构建所述粗化岩相数据中指定岩相数据所对应的波阻抗数据的直方图；数据对提取模块，用于提取所述直方图中的波阻抗数据对，所述波阻抗数据对中包括波阻抗值和对应的网格数量；概率值计算模块，用于根据提取的所述波阻抗数据对，计算各个波阻抗值对应的指定岩相数据的概率值，得到所述指定岩相数据的概率分布数据。

[0014] 进一步地，所述系统还包括：第一设置模块，用于当波阻抗值小于或者等于第一阈值时，将所述波阻抗值对应的指定岩相的概率值设置为第一常数；第二设置模块，用于当波阻抗值大于或者等于第二阈值时，将所述波阻抗值对应的指定岩相的概率值设置为第二常数。

[0015] 由上可见，本申请通过将岩相数据和波阻抗数据加载至预先建立的三维构造模型中，从而可以结合岩相数据和波阻抗数据来建立岩相概率分布模型。具体地，可以确定指定波阻抗值对应的各种岩相概率值，以得到所述各个波阻抗值对应的岩相的概率分布数据。然后分别载入各个岩相数据的概率分布数据。本申请提供的技术方案，能够快速准确地建立岩相概率分布模型。

[0016] 参照后文的说明和附图，详细公开了本申请的特定实施方式，指明了本申请的原理可以被采用的方式。应该理解，本申请的实施方式在范围上并不因而受到限制。在所附权利要求的精神和条款的范围内，本申请的实施方式包括许多改变、修改和等同。

[0017] 针对一种实施方式描述和/或示出的特征可以以相同或类似的方式在一个或更多个其它实施方式中使用，与其它实施方式中的特征相组合，或替代其它实施方式中的特征。

[0018] 应该强调，术语“包括/包含”在本文使用时指特征、整件、步骤或组件的存在，但并不排除一个或更多个其它特征、整件、步骤或组件的存在或附加。

附图说明

[0019] 所包括的附图用来提供对本申请实施方式的进一步的理解，其构成了说明书的一部分，用于例示本申请的实施方式，并与文字描述一起来阐释本申请的原理。显而易见地，下面描述中的附图仅仅是本申请的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。在附图中：

[0020] 图1为本申请实施方式中岩相概率分布模型的建立方法流程图；

[0021] 图2为本申请实施方式中概率分布数据的示意图；

[0022] 图3为本申请实施方式中概率分布模型的示意图；

[0023] 图4为本申请实施方式中岩相概率分布模型的建立系统结构示意图。

具体实施方式

[0024] 为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施方式中的附图，对本申请实施方式中的技术方案进行清楚、完整地描述，显然，所描述的实施方式仅仅是本申请一部分实施方式，而不是全部的实施方式。基于本申请中的实施方式，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施方式，都应当属于本申请保护的范围。

[0025] 请参阅图1，本申请提供的一种岩相概率分布模型的建立方法，包括以下步骤。

[0026] S1：载入已建立的三维构造模型。

[0027] S2：在所述三维构造模型中导入预设岩相数据，并对所述预设岩相数据进行粗化处理，得到粗化岩相数据；所述粗化岩相数据中包括三种岩相数据。

[0028] S3：将预设波阻抗数据体采样至所述三维构造模型中，得到所述三维构造模型对应的波阻抗数据。

[0029] S4：根据所述波阻抗数据以及所述粗化岩相数据，确定各个波阻抗值对应的指定岩相数据的概率值，以得到所述指定岩相数据的概率分布数据。

[0030] S5：分别载入所述粗化岩相数据及对应的各个岩相数据的概率分布数据，以在同一个直方图窗口中展示所述粗化岩相数据的概率分布模型。

[0031] 在本实施方式中，上述步骤可以在软件“Petrel”中执行。具体地，在“Petrel”中可以加载已经建立好的三维构造模型，该三维构造模型可以体现目的地层的结构。

[0032] 在本实施方式中，在加载了所述三维构造模型之后，可以将所述三维构造模型中的所有层设置为相同的网格厚度。具体地，可以通过软件中的“Layering”控件，将所有层的网格厚度设置为5米。

[0033] 在本实施方式中，在加载了三维构造模型之后，可以继续导入已经准备好的岩相数据。该岩相数据可以进行粗化处理，从而得到粗化岩相数据。此外，还可以导入已经准备好的波阻抗数据体。该波阻抗数据体可以进行采样设置。具体地，在软件中可以利用地震重采样(Seismic resampling)的方式，将该波阻抗数据体采样至所述三维构造模型中，得到所述三维构造模型对应的波阻抗数据。

[0034] 在本实施方式中，粗化岩相数据可以包括high、inter、low这三种岩相数据。在处理过程中，可以分别对这三种岩相数据进行处理。以inter岩相数据为例，在软件中可以打开直方图窗口，该直方图窗口中可以展示波阻抗数据与网格数量的对应关系。在实际应用中，可以在直方图窗口的属性栏中，为所述直方图设置数值滤波器(value filter)，并将所述数值滤波器应用于指定岩相数据(inter)中，从而可以构建所述粗化岩相数据中指定岩相数据与所述波阻抗数据之间的直方图。所述直方图中，横坐标表示波阻抗数据，纵坐标表示指定岩相数据的网格数量。

[0035] 在本实施方式中，构建了所述直方图之后，可以基于所述直方图，生成所述指定岩相数据对应的波阻抗分布函数。具体地，可以通过软件中的“Distribution function”(分布函数)功能按键，生成所述指定岩相数据对应的波阻抗分布函数。

[0036] 在本实施方式中，可以从所述波阻抗分布函数中提取出数据对，所述数据对中包括波阻抗值和对应的网格数量。具体地，可以通过软件中的“Spreadsheet”功能按键，弹出该波阻抗分布函数的数据对。然后，可以将这些数据对复制到excel表格中。其中，一列可以表示inter岩相数据对应的波阻抗值，另一列可以表示对应的网格数量。

[0037] 在本实施方式中，通过重复上述的步骤，可以依次提取出high、inter以及low这三种岩相数据的数据对，这些数据对最终可以汇总于excel表格中。根据提取的所述数据对，从而可以计算各个波阻抗值对应的各种岩相数据的概率值，从而得到各种岩相数据的概率分布数据。具体地，可以计算所述指定岩相数据所对应的各个波阻抗值的网格数量在指定波阻抗值总网格数量中的比值，并将所述比值作为各个波阻抗值对应的概率值。举个例子，波阻抗为1时，岩相high有1个，岩相inter有1个，岩相low有1个，那么各自的概率为三分之一，波阻抗为2时，岩相high有1个，岩相inter有2个，岩相low有3个，那么各自的概率为六分之一、三分之一和二分之一。在图2中共有六栏数据，每栏数据中包含了两组数据，其中，前三栏分别包含了不同岩相数据对应的波阻抗值(AI)和网格数量(个数)，而后三栏分别包含了不同岩相数据对应的波阻抗值(AI)和概率。

[0038] 在本实施方式中，针对不同大小的波阻抗值，可以设定一些人为指定的概率值。具体地，当所述波阻抗值小于或者等于第一阈值时，将所述波阻抗值对应的指定岩相数据的概率值设置为第一常数；当所述波阻抗值大于或者等于第二阈值时，将所述波阻抗值对应的指定岩相数据的概率值设置为第二常数；其中，所述第一阈值小于所述第二阈值。请参阅图2，举例来说，当波阻抗值小于或者等于12520时，将第一岩相(high)的概率值设置为1，将第二岩相(inter)和第三岩相(low)的概率值设置为0；当波阻抗值大于或者等于18040时，将第一岩相的概率值设置为0，将第二岩相的概率值设置为0，将第三岩相中的概率值设置为1。

[0039] 这样，通过上述的处理，每种岩相数据便可以得到各自的概率分布数据。在实际应用中，每种岩相数据对应的概率分布数据可以各自保存为文本文件，那么，分别载入所述粗化岩相数据中各个岩相数据的概率分布数据，从而可以在同一个直方图窗口中展示所述粗化岩相数据对应的概率分布模型。具体地，请参阅图3，三个文本文件可以被加载至软件中，然后在直方图窗口，便可以显示由三种曲线构成的概率分布模型。

[0040] 请参阅图4，本申请还提供一种岩相概率分布模型的建立系统，所述系统包括：

[0041] 模型载入单元100，用于载入已建立的三维构造模型；

[0042] 粗化处理单元200，用于在所述三维构造模型中导入预设岩相数据，并对所述预设岩相数据进行粗化处理，得到粗化岩相数据；所述粗化岩相数据中包括三种岩相数据；

[0043] 波阻抗数据采样单元300，用于将预设波阻抗数据体采样至所述三维构造模型中，得到所述三维构造模型对应的波阻抗数据；

[0044] 概率分布数据确定单元400，用于根据所述波阻抗数据以及所述粗化岩相数据，确定各个波阻抗值对应的指定岩相数据的概率值，以得到所述指定岩相数据的概率分布数据；

[0045] 概率分布模型展示单元500，用于分别载入所述粗化岩相数据的概率分布数据，以在同一个直方图窗口中展示所述粗化岩相数据的概率分布模型。

[0046] 在本实施方式中，所述概率分布数据确定单元包括：

[0047] 直方图构建模块，用于构建所述粗化岩相数据中指定岩相数据所对应的波阻抗数据的直方图；

[0048] 数据对提取模块，用于提取所述直方图中的波阻抗数据对，所述波阻抗数据对中包括波阻抗值和对应的网格数量；

[0049] 概率值计算模块，用于根据提取的所述波阻抗数据对，计算各个波阻抗值对应的指定岩相数据的概率值，得到所述指定岩相数据的概率分布数据。

[0050] 在本实施方式中，所述系统还包括：

[0051] 第一设置模块，用于当波阻抗值小于或者等于第一阈值时，将所述波阻抗值对应的指定岩相的概率值设置为第一常数。

[0052] 第二设置模块，用于当波阻抗值大于或者等于第二阈值时，将所述波阻抗值对应的指定岩相的概率值设置为第二常数。

[0053] 由上可见，本申请通过将岩相数据和波阻抗数据加载至预先建立的三维构造模型中，从而可以结合岩相数据和波阻抗数据来建立岩相概率分布模型。具体地，可以确定指定波阻抗值对应的各种岩相概率值，以得到所述各个波阻抗值对应的岩相的概率分布数据。然后分别载入各个岩相数据的概率分布数据。本申请提供的技术方案，能够快速准确地建立岩相概率分布模型。

[0054] 上面对本申请的各种实施方式的描述以描述的目的提供给本领域技术人员。其不旨在是穷举的、或者不旨在将本发明限制于单个公开的实施方式。如上所述，本申请的各种替代和变化对于上述技术所属领域技术人员而言将是显而易见的。因此，虽然已经具体讨论了一些另选的实施方式，但是其它实施方式将是显而易见的，或者本领域技术人员相对容易得出。本申请旨在包括在此已经讨论过的本发明的所有替代、修改、和变化，以及落在上述申请的精神和范围内的其它实施方式。

[0055] 本说明书中的各个实施方式均采用递进的方式描述，各个实施方式之间相同相似的部分互相参见即可，每个实施方式重点说明的都是与其他实施方式的不同之处。

[0056] 虽然通过实施方式描绘了本申请，本领域普通技术人员知道，本申请有许多变形和变化而不脱离本申请的精神，希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。