CN104155689A - 碳酸盐岩缝洞雕刻方法及装置 - Google Patents

Abstract

本发明实施例提供了一种碳酸盐岩缝洞雕刻方法及装置，其中，该方法包括：将碳酸盐岩的工区平面图划分为多个缝洞带，并将每个缝洞带划分为多个缝洞单元，所述缝洞带为不规则的闭合区域，所述缝洞单元是不规则的闭合区域；根据预设缝洞振幅值范围，确定每个缝洞带中每个缝洞单元内的缝洞体；计算每个缝洞带中每个缝洞单元内每个缝洞体的体积和储量；对于每个缝洞单元，将该缝洞单元内缝洞体的体积累加，得到该缝洞单元的体积，将该缝洞单元内缝洞体的储量累加，得到该缝洞单元的储量。该方案可以一次性将缝洞带内的缝洞体全部雕刻出来，并量化缝洞体的体积和储量。

Description

碳酸盐岩缝洞雕刻方法及装置

技术领域

本发明涉及地球物理勘探地震资料计算机识别技术领域，特别涉及一种碳酸盐岩缝洞雕刻方法及装置。

背景技术

在地球物理勘探处理和解释过程中，对碳酸盐岩缝洞进行雕刻、量化统计并形成报表以满足勘探生产开发的需要，是碳酸盐岩区域必须解决的问题。

只有将碳酸盐岩缝洞雕刻出来，统计分布范围并量化其体积、储量才能进行下一步的勘探生产和开发。但是，目前针对碳酸盐岩这样的特殊地质目标，常规解释软件无法满足量化缝洞的体积、储量等信息需要，无法实现工业化生产。

发明内容

本发明实施例提供了一种碳酸盐岩缝洞雕刻方法，以满足量化缝洞的体积、储量的需要，该方法包括：将碳酸盐岩的工区平面图划分为多个缝洞带，并将每个缝洞带划分为多个缝洞单元，所述缝洞带为不规则的闭合区域，所述缝洞单元是不规则的闭合区域；根据预设缝洞振幅值范围，确定每个缝洞带中每个缝洞单元内的缝洞体；计算每个缝洞带中每个缝洞单元内每个缝洞体的体积和储量；对于每个缝洞单元，将该缝洞单元内缝洞体的体积累加，得到该缝洞单元的体积，将该缝洞单元内缝洞体的储量累加，得到该缝洞单元的储量。

在一个实施例中，根据预设缝洞振幅值范围，确定每个缝洞带中每个缝洞单元内的缝洞体，包括：对于每个缝洞带，获取该缝洞带的线号个数、道号个数和样点个数；将该缝洞带内每个种子点的瞬时振幅值存储在三维数组中，并为瞬时振幅值在所述预设缝洞振幅值范围内的种子点设置第一特征值，为瞬时振幅值在所述预设缝洞振幅值范围外的种子点设置第二特征值；以该缝洞带内第一个具有第一特征值的种子点为中心定义一个缝洞体；在该种子点的上方、下方、左方、右方、前方和后方与该种子点联通的种子点中，将具有第一特征值的种子点确定为属于所述一个缝洞体的种子点，并以该确定出的种子点为中心，从与该确定出的种子点上方、下方、左方、右方、前方和后方联通的种子点中，将具有第一特征值的种子点确定为属于所述一个缝洞体的种子点，直至与该确定出的种子点联通的种子点中，不包括具有第一特征值的种子点。

在一个实施例中，计算每个缝洞带中每个缝洞单元内缝洞体的体积，包括：根据工区的线线间距、道号间距和采样间隔，计算缝洞体内每个种子点的体积；累加缝洞体内每个种子点的体积，得到缝洞体的体积。

在一个实施例中，计算每个缝洞带中每个缝洞单元内缝洞体的储量，包括：根据缝洞体内每个种子点的体积和缝洞体内每个种子点的孔隙度数据值，计算缝洞体的有效体积；根据含油饱和度、密度和缝洞体的有效体积，计算缝洞体的储量。

在一个实施例中，根据预设缝洞振幅值范围，确定每个缝洞带中每个缝洞单元内的缝洞体，包括：根据用户输入的预设缝洞振幅值范围，确定每个缝洞带中每个缝洞单元内的缝洞体；所述方法还包括：确定每个缝洞带中每个缝洞单元内的缝洞体之前，对所述工区的地震资料数据进行瞬时振幅变换，得到瞬时振幅数据，所述瞬时振幅数据为用户确定所述预设缝洞振幅值范围提供依据。

在一个实施例中，还包括：输出外部文件，所述外部文件包括每个缝洞体的体积、储量、位置、面积和厚度。

本发明实施例提供了一种碳酸盐岩缝洞雕刻装置，以满足量化缝洞的体积、储量的需要。该装置包括：划分模块，用于将碳酸盐岩的工区平面图划分为多个缝洞带，并将每个缝洞带划分为多个缝洞单元，所述缝洞带为不规则的闭合区域，所述缝洞单元是不规则的闭合区域；确定模块，用于根据预设缝洞振幅值范围，确定每个缝洞带中每个缝洞单元内的缝洞体；计算模块，用于计算每个缝洞带中每个缝洞单元内每个缝洞体的体积和储量；累加模块，用于对于每个缝洞单元，将该缝洞单元内缝洞体的体积累加，得到该缝洞单元的体积，将该缝洞单元内缝洞体的储量累加，得到该缝洞单元的储量。

在一个实施例中，所述确定模块，包括：获取单元，用于对于每个缝洞带，获取该缝洞带的线号个数、道号个数和样点个数；设置单元，用于将该缝洞带内每个种子点的瞬时振幅值存储在三维数组中，并为瞬时振幅值在所述预设缝洞振幅值范围内的种子点设置第一特征值，为瞬时振幅值在所述预设缝洞振幅值范围外的种子点设置第二特征值；定义单元，用于以该缝洞带内第一个具有第一特征值的种子点为中心定义一个缝洞体；雕刻单元，用于在该种子点的上方、下方、左方、右方、前方和后方与该种子点联通的种子点中，将具有第一特征值的种子点确定为属于所述一个缝洞体的种子点，并以该确定出的种子点为中心，从与该确定出的种子点上方、下方、左方、右方、前方和后方联通的种子点中，将具有第一特征值的种子点确定为属于所述一个缝洞体的种子点，直至与该确定出的种子点联通的种子点中，不包括具有第一特征值的种子点。

在一个实施例中，所述计算模块，包括：第一计算单元，用于根据工区的线线间距、道号间距和采样间隔，计算缝洞体内每个种子点的体积；第二计算单元，用于累加缝洞体内每个种子点的体积，得到缝洞体的体积。

在一个实施例中，所述计算模块，还包括：第三计算单元，用于根据缝洞体内每个种子点的体积和缝洞体内每个种子点的孔隙度数据值，计算缝洞体的有效体积；第四计算单元，用于根据含油饱和度、密度和缝洞体的有效体积，计算缝洞体的储量。

在一个实施例中，所述确定模块，具体用于根据用户输入的预设缝洞振幅值范围，确定每个缝洞带中每个缝洞单元内的缝洞体；所述装置还包括：变换模块，用于根据用户输入的预设缝洞振幅值范围，确定每个缝洞带中每个缝洞单元内的缝洞体之前，对所述工区的地震资料数据进行瞬时振幅变换，得到瞬时振幅数据，所述瞬时振幅数据为用户确定所述预设缝洞振幅值范围提供依据。

在一个实施例中，还包括：输出模块，用于输出外部文件，所述外部文件的内容包括每个缝洞体的体积、储量、位置、面积和厚度。

在本发明实施例中，通过将碳酸盐岩的工区平面图划分为多个缝洞带，并将每个缝洞带划分为多个缝洞单元，然后，根据预设缝洞振幅值范围，确定每个缝洞带中每个缝洞单元内的缝洞体，并计算每个缝洞单元内缝洞体的体积和储量，最后得到每个缝洞单元的体积和储量，实现了一次性将缝洞带内的缝洞体全部雕刻出来，并量化缝洞体的体积和储量，使得可以对缝洞体进行快速、便捷地雕刻，并使解释人员对缝洞体能够有一个更加直观的认识，以满足下一步勘探生产的需求。

附图说明

此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中：

图1是本发明实施例提供的一种碳酸盐岩缝洞雕刻方法的流程图；

图2是本发明实施例提供的一种缝洞带划分示意图；

图3是本发明实施例提供的一种缝洞单元划分示意图；

图4是本发明实施例提供的一种缝洞的联通示意图；

图5是本发明实施例提供的一种缝洞的三维显示效果图；

图6是本发明实施例提供的一种碳酸盐岩缝洞雕刻装置的结构框图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施方式和附图，对本发明做进一步详细说明。在此，本发明的示意性实施方式及其说明用于解释本发明，但并不作为对本发明的限定。

在本发明实施例中，提供了一种碳酸盐岩缝洞雕刻方法，如图1所示，该方法包括：

步骤101：将碳酸盐岩的工区平面图划分为多个缝洞带，并将每个缝洞带划分为多个缝洞单元，所述缝洞带为不规则的闭合区域，所述缝洞单元是不规则的闭合区域；

步骤102：根据预设缝洞振幅值范围，确定每个缝洞带中每个缝洞单元内的缝洞体；

步骤103：计算每个缝洞带中每个缝洞单元内每个缝洞体的体积和储量；

步骤104：对于每个缝洞单元，将该缝洞单元内缝洞体的体积累加，得到该缝洞单元的体积，将该缝洞单元内缝洞体的储量累加，得到该缝洞单元的储量。

由图1所示的流程可知，在本发明实施例中，通过将碳酸盐岩的工区平面图划分为多个缝洞带，并将每个缝洞带划分为多个缝洞单元，然后，根据预设缝洞振幅值范围，确定每个缝洞带中每个缝洞单元内的缝洞体，并计算每个缝洞单元内缝洞体的体积和储量，最后得到每个缝洞单元的体积和储量，实现了一次性将缝洞带内的缝洞体全部雕刻出来，并量化缝洞体的体积和储量，使得可以对缝洞体进行快速、便捷地雕刻，并使解释人员对缝洞体能够有一个更加直观的认识，以满足下一步勘探生产的需求。

具体实施时，可以通过以下方式将工区平面图划分为多个缝洞带，并将每个缝洞带划分为多个缝洞单元：例如，解释用户在工区平面图上用鼠标圈出闭合缝洞带，如图2所示，划分为缝洞带1和缝洞带2，这些缝洞带用来圈定需要雕刻缝洞的工区边界，表现为多个闭合的多边形，每个缝洞带的名字唯一；使用同样方法定义缝洞单元，用鼠标在缝洞带内圈出缝洞单元，如图3所示，缝洞单元指缝洞带内需要进一步统计缝洞信息的边界范围。

具体实施时，可以通过以下步骤来实现根据预设缝洞振幅值范围，确定每个缝洞带中每个缝洞单元内的缝洞体：

首先，对于每个缝洞带，获取该缝洞带的线号个数、道号个数和样点个数；例如，搜索某个指定的缝洞带所在的工区线道号最小范围(line_min，line_max，cmp_min，cmp_max)，搜索目的层所在的采样索引范围(sample_min，sample_max)，则该缝洞带的线号个数、道号个数、样点个数分别为：

line_count＝line_max–line_min+1；

cmp_count＝cmp_max–cmp_min+1；

sample_count＝sample_max–sample_min+1。

其次，将该缝洞带内每个种子点的瞬时振幅值存储在三维数组中，并为瞬时振幅值在所述预设缝洞振幅值范围内的种子点设置第一特征值，为瞬时振幅值在所述预设缝洞振幅值范围外的种子点设置第二特征值；例如，定义bool类型三维数组，以存放该缝洞带内每个种子点的瞬时振幅特征值，并将特征数组中每个特征点初始化为false：

bool flag[line_count，cmp_count，sample_count]

对每一个flag元素赋值为false；

依次读取该缝洞带范围内每个种子点的瞬时振幅数据，根据瞬时振幅值是否位于用户定义的缝洞振幅范围(即上述预设缝洞振幅值范围)内，对振幅特征数组进行赋值，假设某个样点线号为line_no，道号为cmp_no，样点号为sample_no，那么该点在所加载数据范围内的线索引、道索引、样点索引分别为：

line_idx＝line_no–line_min；

cmp_idx＝cmp_no–cmp_min；

sample_idx＝sample_no–sample_min；

该位置的特征值根据用户定义的缝洞振幅范围确定，如果该样点的瞬时振幅值在所述预设缝洞振幅值范围内，则为该样点设置的特征值为true(即第一特征值)，否则，为该样点设置的特征值为false(第二特征值)。

再次，以该缝洞带内第一个具有第一特征值的种子点为中心定义一个缝洞体；例如，分别对该缝洞带内每一个flag特征点进行判断，如果该样点的特征值为false，则忽略该特征点，否则定义一个缝洞体，并对位于(line_i，cmp_i，sample_i)位置的种子特征点执行雕刻操作。

最后，对具有第一特征值的种子点执行以下步骤进行缝洞的雕刻：在该种子点的上方、下方、左方、右方、前方和后方与该种子点联通的种子点中，将具有第一特征值的种子点确定为属于所述一个缝洞体的种子点，并以该确定出的种子点为中心，从与该确定出的种子点上方、下方、左方、右方、前方和后方联通的种子点中，将具有第一特征值的种子点确定为属于所述一个缝洞体的种子点，直至与该确定出的种子点联通的种子点中，不包括具有第一特征值的种子点。例如，首先将flag(line_i，cmp_i，sample_i)的特征值赋值为false，防止该点被再次分析，并将种子特征点(line_i，cmp_i，sample_i)的索引保存到缝洞体中，取出位于下列位置种子点的特征值：

(line_i,cmp_i,sample_i-1)

(line_i,cmp_i,sample_i+1)

(line_i,cmp_i-1,sample_i)

(line_i,cmp_i+1,sample_i)

(line_i-1,cmp_i,sample_i)

(line_i+1,cmp_i,sample_i)

上面点分别表示种子特征点(line_i，cmp_i，sample_i)上、下、左、右、前、后六个点，即与种子特征点(line_i，cmp_i，sample_i)上方、下方、左方、右方、前方和后方联通的种子点(如图4所示)，将这些种子点中特征值为true的种子点，确定为同一个缝洞体，将这些种子点中特征值为false的种子点，确定为不是同一个缝洞体，并以确定出的种子点为中心，从与该确定出的种子点联通的种子点中，将具有第一特征值(即true)的种子点确定为属于一个缝洞体的种子点，直至与该确定出的种子点联通的种子点中，不包括具有第一特征值的种子点，该一个缝洞体雕刻完毕。继续判断缝洞带内其他种子点的特征值，并以具有第一特征值的种子点为中心定义另一个缝洞体，对该具有第一特征值的种子点进行缝洞雕刻的操作。

具体的，根据预设缝洞振幅值范围，确定每个缝洞带中每个缝洞单元内的缝洞体，包括：根据用户输入的预设缝洞振幅值范围，确定每个缝洞带中每个缝洞单元内的缝洞体，为了便于用户给出预设缝洞振幅值范围，在本实施例中，根据用户输入的预设缝洞振幅值范围，确定每个缝洞带中每个缝洞单元内的缝洞体之前，还包括：对所述工区的地震资料数据进行瞬时振幅变换，得到瞬时振幅数据，所述瞬时振幅数据为用户确定所述预设缝洞振幅值范围提供依据。因为缝洞在地震资料中的表现是一个个的最大或最小振幅范围，即缝洞的振幅特征为大于某个振幅值和小于某个振幅值两个区间，这样的振幅值范围很难让用户指定所雕刻的振幅值范围，所以在做缝洞雕刻之前，对地震资料做瞬时振幅变换。

首先，对地震数据进行Hilbert变换，假设地震振幅值为s(t)：

$s_r\left(t\right)=H\left\{s\right\}=h\left(t\right)\×s\left(t\right)={\∫}_{-\∞}^{\∞}s\left(\mathrm{\τ}\right)h(t-\mathrm{\τ})\mathrm{d\τ}=\frac{1}{\mathrm{\π}}{\∫}_{-\∞}^{\∞}\frac{s\left(\mathrm{\τ}\right)}{t-\mathrm{\τ}}\mathrm{d\τ}$

则瞬时振幅计算公式如下：

$\hat{A}\left(t\right)=\sqrt{s_r^2\left(t\right)+s^2\left(t\right)},$

对地震资料做了瞬时振幅变换后，其最大、最小振幅区间都变成了某个最大振幅区间，这样的振幅区间方便用户指定某个振幅范围(即预设缝洞振幅值范围)，比如大于某个振幅值的样点特征。

在雕刻出缝洞体之后，还需要量化缝洞体的体积，在本实施例中，计算每个缝洞带中每个缝洞单元内缝洞体的体积，包括：根据缝工区的线线间距、道号间距和采样间隔，计算缝洞体内每个种子点的体积；累加缝洞体内每个种子点的体积，得到缝洞体的体积。具体的，根据工区线线间距(line_distance)道号间距(cmp_distance)采样间隔(sampel_interval)，计算每个有效点(即缝洞体内每个种子点)的体积V_i，其中line_distance，cmp_distance单位为米，如果是深度域类型数据sample_interval单位为米，而时间域类型数据则需要将时间间隔变换到深度域：

其中velocity指当前目地层的速度，于某个特定目的层velocity可以认为是一个常量，该参数由用户指定。

V_i＝line_distance×cmp_distance×sample_distance

那么每个缝洞的体积则为：V＝V_i×缝洞体内种子点的个数。在雕刻出缝洞体之后，还需要量化缝洞体的体积，在本实施例中，计算每个缝洞带中每个缝洞单元内缝洞体的储量，包括：根据缝洞体内每个种子点的体积和缝洞体内每个种子点的孔隙度数据值，计算缝洞体的有效体积；根据含油饱和度、密度和缝洞体的有效体积，计算缝洞体的储量。例如，对缝洞体内的每个种子点读取相应位置的孔隙度数据值pori，该参数可以由外部作为孔隙度数据体输入，根据孔隙度数据值计算缝洞体的有效体积：V_valid＝∫V_i×por_i，

那么缝洞体的储量则为：R＝V_valid×含油饱和度×密度，其中，含油饱和度和密度可以通过测井获得，由用户输入即可。

最后，为了将雕刻结果呈现给用户，并便于供其他软件打开，在本实施例中，上述碳酸盐岩缝洞雕刻方法，还包括：输出外部文件，所述外部文件包括每个缝洞体的体积、储量、位置、面积和厚度，该外部文件可以是报表的形式，如下表1所示。其中，缝洞体的位置是缝洞体中心的种子点所在的位置，缝洞体的面积是缝洞体从上向下投影的面积。此外，还可以使用三维可视化显示雕刻结果(如图5所示)，以更加直观的方式观察缝洞的位置、体积等信息。

表1

以下以具体示例详细描述上述碳酸盐岩缝洞雕刻方法，该方法包括以下步骤：

步骤1：对所述工区的地震资料数据进行瞬时振幅变换，得到瞬时振幅数据；因为缝洞在地震资料中的表现是一个个的最大或最小振幅范围，即缝洞的振幅特征为大于某个振幅值和小于某个振幅值两个区间，这样的振幅值范围很难让用户指定所雕刻的振幅值范围，所以在做缝洞雕刻之前，对地震资料做瞬时振幅变换。

首先，对地震数据进行Hilbert变换，假设地震振幅值为s(t)：

$s_r\left(t\right)=H\left\{s\right\}=h\left(t\right)\×s\left(t\right)={\∫}_{-\∞}^{\∞}s\left(\mathrm{\τ}\right)h(t-\mathrm{\τ})\mathrm{d\τ}=\frac{1}{\mathrm{\π}}{\∫}_{-\∞}^{\∞}\frac{s\left(\mathrm{\τ}\right)}{t-\mathrm{\τ}}\mathrm{d\τ}$

则瞬时振幅计算公式如下：

$\hat{A}\left(t\right)=\sqrt{s_r^2\left(t\right)+s^2\left(t\right)}$

对地震资料做了瞬时变换后，其最大、最小振幅区间都变成了某个最大振幅区间，这样的振幅区间方便用户指定某个振幅范围(即预设缝洞振幅值范围)，比如大于某个振幅值的样点特征。

步骤2：将碳酸盐岩的工区平面图划分为多个缝洞带，并将每个缝洞带划分为多个缝洞单元；解释用户在平面图用鼠标圈出闭合缝洞带，如图2所示，划分为缝洞带1和缝洞带2，这些缝洞带用来圈定需要雕刻缝洞的工区边界，表现为多个闭合的多边形，每个缝洞带的名字唯一，使用同样方法定义缝洞单元，用鼠标在缝洞带内圈出缝洞单元，如图3所示，缝洞单元指缝洞带内需要进一步统计缝洞信息的边界范围。

步骤3：对于每个缝洞带，获取该缝洞带的线号个数、道号个数和样点个数；搜索某个指定的缝洞带所在的工区线道号最小范围(line_min，line_max，cmp_min，cmp_max)，搜索目的层所在的采样索引范围(sample_min，sample_max)，则该缝洞带的线号个数、道号个数、样点个数分别为：

line_count＝line_max–line_min+1；

cmp_count＝cmp_max–cmp_min+1；

sample_count＝sample_max–sample_min+1；

定义bool类型三维数组，以将缝洞带内每个种子点的瞬时振幅值存储在三维数组中，并将特征数组每个特征点初始化为false：

bool flag[line_count,cmp_count,sample_count]

对每一个flag元素赋值为false。

步骤4：为瞬时振幅值在所述预设缝洞振幅值范围内的种子点设置第一特征值，为瞬时振幅值在所述预设缝洞振幅值范围外的种子点设置第二特征值；依次读取该缝洞带范围内每个种子点的瞬时振幅数据，根据瞬时振幅值是否位于用户定义的缝洞振幅范围(即上述预设缝洞振幅值范围)内，对振幅特征数组进行赋值，假设某个样点线号为line_no，道号为cmp_no，样点号为sample_no，那么该点在所加载数据范围内的线索引、道索引、样点索引分别为：

line_idx＝line_no–line_min；

cmp_idx＝cmp_no–cmp_min；

sample_idx＝sample_no–sample_min；

该位置的特征值根据用户定义的缝洞振幅范围确定，如果该样点的瞬时振幅值在所述预设缝洞振幅值范围内，则为该样点设置的特征值为true(即第一特征值)，否则，为该样点设置的特征值为false(第二特征值)；

步骤5：以该缝洞带内第一个具有第一特征值的种子点为中心定义一个缝洞体；分别对每一个flag特征点进行判断，如果该样点为false，则忽略该特征点，否则定义一个缝洞体，并对位于(line_i，cmp_i，sample_i)位置的种子特征点执行步骤6的雕刻操作。

步骤6：首先将flag(line_i，cmp_i，sample_i))的特征值赋值为false，防止该点被再次分析，并将种子特征点(line_i，cmp_i，sample_i)的索引保存到缝洞体中，取出位于下列位置种子点的特征值：

(line_i,cmp_i,sample_i-1)

(line_i,cmp_i,sample_i+1)

(line_i,cmp_i-1,sample_i)

(line_i,cmp_i+1,sample_i)

(line_i-1,cmp_i,sample_i)

(line_i+1,cmp_i,sample_i)

上面点分别表示种子特征点(line_i，cmp_i，sample_i)上、下、左、右、前、后六个点，即与种子特征点(line_i，cmp_i，sample_i)联通的种子点(如图4所示)，将这些种子点中特征值为true的种子点，确定为同一个缝洞体，将这些种子点中特征值为false的种子点，确定为不是同一个缝洞体，并以确定出的种子点为中心执行步骤6，从与该确定出的种子点联通的种子点中，将具有第一特征值(即true)的种子点确定为属于一个缝洞体的种子点，直至与该确定出的种子点联通的种子点中，不包括具有第一特征值的种子点，该一个缝洞体雕刻完毕。继续判断缝洞带内其他种子点的特征值，并以具有第一特征值的种子点为中心定义另一个缝洞体，对该具有第一特征值的种子点执行步骤6的缝洞雕刻操作。

步骤7：计算缝洞体的体积；根据工区线线间距(line_distance)、道号间距(cmp_distance)和采样间隔(sampel_interval)，计算每个有效点(即缝洞体内每个种子点)的体积Vi，其中line_distance，cmp_distance单位为米，如果是深度域类型数据sample_interval单位为米，而时间域类型数据则需要将时间间隔变换到深度域：

其中velocity指当前目的层的速度，于某个特定目的层velocity可以认为是一个常量，该参数由用户指定。

V_i＝line_distance×cmp_distance×sample_distance，

那么每个缝洞的体积则为：V＝V_i×缝洞体内种子点的个数。

步骤8：计算缝洞体的储量；对缝洞体内的每个点读取相应位置的孔隙度数据值pori，该参数可以由外部作为孔隙度数据体输入，根据孔隙度数据值计算缝洞体的有效体积：

V_valid＝∫V_i×pori，

那么缝洞体的储量则为：

R＝V_valid×含油饱和度×密度，其中，含油饱和度和密度可以通过测井获得，由用户输入即可。

步骤9：将该缝洞单元内缝洞体的体积累加，得到该缝洞单元的体积，将该缝洞单元内缝洞体的储量累加，得到该缝洞单元的储量；对缝洞带内的缝洞体依次进行判断，如果落在某个缝洞单元(闭合多边形区域)，则将体积累加到缝洞单元的体积中。

V_unit＝∫V_i

同样，储量也统计到缝洞单元储量中,

V_unit＝∫R_i

为了将雕刻结果呈现给用户，并便于供其他软件打开，上述碳酸盐岩缝洞雕刻方法，可以输出外部文件，改外部文件的内容包括每个缝洞体的体积、储量、位置、面积和厚度，该外部文件可以是报表的形式，如上表1所示。其中，缝洞体的位置是缝洞体中心的种子点所在的位置，缝洞体的面积是缝洞体从上向下投影的面积。此外，还可以使用三维可视化显示雕刻结果(如图5所示)，以更加直观的方式观察缝洞的位置、体积等信息。

基于同一发明构思，本发明实施例中还提供了一种碳酸盐岩缝洞雕刻装置，如下面的实施例所述。由于碳酸盐岩缝洞雕刻装置解决问题的原理与碳酸盐岩缝洞雕刻方法相似，因此碳酸盐岩缝洞雕刻装置的实施可以参见碳酸盐岩缝洞雕刻方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。

图6是本发明实施例的碳酸盐岩缝洞雕刻装置的一种结构框图，如图6所示，包括：划分模块601、确定模块602、计算模块603和累加模块604，下面对该结构进行说明。

划分模块601，用于将碳酸盐岩的工区平面图划分为多个缝洞带，并将每个缝洞带划分为多个缝洞单元，所述缝洞带为不规则的闭合区域，所述缝洞单元是不规则的闭合区域；

确定模块602，与划分模块601连接，用于根据预设缝洞振幅值范围，确定每个缝洞带中每个缝洞单元内的缝洞体；

计算模块603，与确定模块602连接，用于计算每个缝洞带中每个缝洞单元内每个缝洞体的体积和储量；

累加模块604，与计算模块603连接，用于对于每个缝洞单元，将该缝洞单元内缝洞体的体积累加，得到该缝洞单元的体积，将该缝洞单元内缝洞体的储量累加，得到该缝洞单元的储量。

在一个实施例中，所述确定模块602，包括：获取单元，用于对于每个缝洞带，获取该缝洞带的线号个数、道号个数和样点个数；设置单元，与获取单元连接，用于将该缝洞带内每个种子点的瞬时振幅值存储在三维数组中，并为瞬时振幅值在所述预设缝洞振幅值范围内的种子点设置第一特征值，为瞬时振幅值在所述预设缝洞振幅值范围外的种子点设置第二特征值；定义单元，与设置单元连接，用于以该缝洞带内第一个具有第一特征值的种子点为中心定义一个缝洞体；雕刻单元，与定义单元连接，用于在该种子点的上方、下方、左方、右方、前方和后方与该种子点联通的种子点中，将具有第一特征值的种子点确定为属于所述一个缝洞体的种子点，并以该确定出的种子点为中心，从与该确定出的种子点上方、下方、左方、右方、前方和后方联通的种子点中，将具有第一特征值的种子点确定为属于所述一个缝洞体的种子点，直至与该确定出的种子点联通的种子点中，不包括具有第一特征值的种子点。

在一个实施例中，所述计算模块603，包括：第一计算单元，用于根据工区的线线间距、道号间距和采样间隔，计算缝洞体内每个种子点的体积；第二计算单元，与第一计算单元连接，用于累加缝洞体内每个种子点的体积，得到缝洞体的体积。

在一个实施例中，所述计算模块603，还包括：第三计算单元，用于根据缝洞体内每个种子点的体积和缝洞体内每个种子点的孔隙度数据值，计算缝洞体的有效体积；第四计算单元，与第三计算单元连接，用于根据含油饱和度、密度和缝洞体的有效体积，计算缝洞体的储量。

在一个实施例中，所述确定模块602，具体用于根据用户输入的预设缝洞振幅值范围，确定每个缝洞带中每个缝洞单元内的缝洞体；所述装置还包括：变换模块，用于根据用户输入的预设缝洞振幅值范围，确定每个缝洞带中每个缝洞单元内的缝洞体之前，对所述工区的地震资料数据进行瞬时振幅变换，得到瞬时振幅数据，所述瞬时振幅数据为用户确定所述预设缝洞振幅值范围提供依据。

在一个实施例中，还包括：输出模块，用于输出外部文件，所述外部文件的内容包括每个缝洞体的体积、储量、位置、面积和厚度。

在本发明实施例中，通过将碳酸盐岩的工区平面图划分为多个缝洞带，并将每个缝洞带划分为多个缝洞单元，然后，根据预设缝洞振幅值范围，确定每个缝洞带中每个缝洞单元内的缝洞体，并计算每个缝洞单元内缝洞体的体积和储量，最后得到每个缝洞单元的体积和储量，实现了一次性将缝洞带内的缝洞体全部雕刻出来，并量化缝洞体的体积和储量，使得可以对缝洞体进行快速、便捷地雕刻，并使解释人员对缝洞体能够有一个更加直观的认识，以满足下一步勘探生产的需求。

显然，本领域的技术人员应该明白，上述的本发明实施例的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，并且在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明实施例不限制于任何特定的硬件和软件结合。

以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (12)

1.一种碳酸盐岩缝洞雕刻方法，其特征在于，包括：

将碳酸盐岩的工区平面图划分为多个缝洞带，并将每个缝洞带划分为多个缝洞单元，所述缝洞带为不规则的闭合区域，所述缝洞单元是不规则的闭合区域；

根据预设缝洞振幅值范围，确定每个缝洞带中每个缝洞单元内的缝洞体；

计算每个缝洞带中每个缝洞单元内每个缝洞体的体积和储量；

对于每个缝洞单元，将该缝洞单元内缝洞体的体积累加，得到该缝洞单元的体积，将该缝洞单元内缝洞体的储量累加，得到该缝洞单元的储量。

2.如权利要求1所述碳酸盐岩缝洞雕刻方法，其特征在于，根据预设缝洞振幅值范围，确定每个缝洞带中每个缝洞单元内的缝洞体，包括：

对于每个缝洞带，获取该缝洞带的线号个数、道号个数和样点个数；

将该缝洞带内每个种子点的瞬时振幅值存储在三维数组中，并为瞬时振幅值在所述预设缝洞振幅值范围内的种子点设置第一特征值，为瞬时振幅值在所述预设缝洞振幅值范围外的种子点设置第二特征值；

以该缝洞带内第一个具有第一特征值的种子点为中心定义一个缝洞体；

在该种子点的上方、下方、左方、右方、前方和后方与该种子点联通的种子点中，将具有第一特征值的种子点确定为属于所述一个缝洞体的种子点，并以该确定出的种子点为中心，从与该确定出的种子点上方、下方、左方、右方、前方和后方联通的种子点中，将具有第一特征值的种子点确定为属于所述一个缝洞体的种子点，直至与该确定出的种子点联通的种子点中，不包括具有第一特征值的种子点。

3.如权利要求2所述碳酸盐岩缝洞雕刻方法，其特征在于，计算每个缝洞带中每个缝洞单元内缝洞体的体积，包括：

根据工区的线线间距、道号间距和采样间隔，计算缝洞体内每个种子点的体积；

累加缝洞体内每个种子点的体积，得到缝洞体的体积。

4.如权利要求3所述碳酸盐岩缝洞雕刻方法，其特征在于，计算每个缝洞带中每个缝洞单元内缝洞体的储量，包括：

根据缝洞体内每个种子点的体积和缝洞体内每个种子点的孔隙度数据值，计算缝洞体的有效体积；

根据含油饱和度、密度和缝洞体的有效体积，计算缝洞体的储量。

5.如权利要求1至4中任一项所述碳酸盐岩缝洞雕刻方法，其特征在于，根据预设缝洞振幅值范围，确定每个缝洞带中每个缝洞单元内的缝洞体，包括：根据用户输入的预设缝洞振幅值范围，确定每个缝洞带中每个缝洞单元内的缝洞体；

所述方法还包括：

确定每个缝洞带中每个缝洞单元内的缝洞体之前，对所述工区的地震资料数据进行瞬时振幅变换，得到瞬时振幅数据，所述瞬时振幅数据为用户确定所述预设缝洞振幅值范围提供依据。

6.如权利要求1至4中任一项所述碳酸盐岩缝洞雕刻方法，其特征在于，还包括：

输出外部文件，所述外部文件包括每个缝洞体的体积、储量、位置、面积和厚度。

7.一种碳酸盐岩缝洞雕刻装置，其特征在于，包括：

划分模块，用于将碳酸盐岩的工区平面图划分为多个缝洞带，并将每个缝洞带划分为多个缝洞单元，所述缝洞带为不规则的闭合区域，所述缝洞单元是不规则的闭合区域；

确定模块，用于根据预设缝洞振幅值范围，确定每个缝洞带中每个缝洞单元内的缝洞体；

计算模块，用于计算每个缝洞带中每个缝洞单元内每个缝洞体的体积和储量；

累加模块，用于对于每个缝洞单元，将该缝洞单元内缝洞体的体积累加，得到该缝洞单元的体积，将该缝洞单元内缝洞体的储量累加，得到该缝洞单元的储量。

8.如权利要求7所述碳酸盐岩缝洞雕刻装置，其特征在于，所述确定模块，包括：

获取单元，用于对于每个缝洞带，获取该缝洞带的线号个数、道号个数和样点个数；

设置单元，用于将该缝洞带内每个种子点的瞬时振幅值存储在三维数组中，并为瞬时振幅值在所述预设缝洞振幅值范围内的种子点设置第一特征值，为瞬时振幅值在所述预设缝洞振幅值范围外的种子点设置第二特征值；

定义单元，用于以该缝洞带内第一个具有第一特征值的种子点为中心定义一个缝洞体；

雕刻单元，用于在该种子点的上方、下方、左方、右方、前方和后方与该种子点联通的种子点中，将具有第一特征值的种子点确定为属于所述一个缝洞体的种子点，并以该确定出的种子点为中心，从与该确定出的种子点上方、下方、左方、右方、前方和后方联通的种子点中，将具有第一特征值的种子点确定为属于所述一个缝洞体的种子点，直至与该确定出的种子点联通的种子点中，不包括具有第一特征值的种子点。

9.如权利要求8所述碳酸盐岩缝洞雕刻装置，其特征在于，所述计算模块，包括：

第一计算单元，用于根据工区的线线间距、道号间距和采样间隔，计算缝洞体内每个种子点的体积；

第二计算单元，用于累加缝洞体内每个种子点的体积，得到缝洞体的体积。

10.如权利要求9所述碳酸盐岩缝洞雕刻装置，其特征在于，所述计算模块，还包括：

第三计算单元，用于根据缝洞体内每个种子点的体积和缝洞体内每个种子点的孔隙度数据值，计算缝洞体的有效体积；

第四计算单元，用于根据含油饱和度、密度和缝洞体的有效体积，计算缝洞体的储量。

11.如权利要求7至10中任一项所述碳酸盐岩缝洞雕刻装置，其特征在于，所述确定模块，具体用于根据用户输入的预设缝洞振幅值范围，确定每个缝洞带中每个缝洞单元内的缝洞体；

所述装置还包括：

变换模块，用于根据用户输入的预设缝洞振幅值范围，确定每个缝洞带中每个缝洞单元内的缝洞体之前，对所述工区的地震资料数据进行瞬时振幅变换，得到瞬时振幅数据，所述瞬时振幅数据为用户确定所述预设缝洞振幅值范围提供依据。

12.如权利要求7至10中任一项所述碳酸盐岩缝洞雕刻装置，其特征在于，还包括：

输出模块，用于输出外部文件，所述外部文件的内容包括每个缝洞体的体积、储量、位置、面积和厚度。