CN108074282B - 基于三维地质建模软件的物性参数的厚度校正方法 - Google Patents

Abstract

本发明提供一种基于三维地质建模软件的物性参数的厚度校正方法，包括：步骤1，将储层厚度数据统一格式，并进行整数编码；步骤2，将厚度属性编码正确加载到三维地质建模软件中；步骤3，利用厚度属性编码计算得到准确的储层厚度属性曲线数据；步骤4，利用建立的厚度校正的解释模型，对孔隙度、渗透率等物性参数进行厚度校正计算。该基于三维地质建模软件的物性参数的厚度校正方法，通过计算和加载，准确地构建了厚度属性曲线，解决了三维地质建模中缺少储层厚度属性曲线的问题，在三维地质建模软件中实现了储层厚度的校正，计算得到了更为准确的孔隙度、渗透率等储层属性数据，并为利用厚度属性曲线进行的其他计算提供了基础数据。

Description

基于三维地质建模软件的物性参数的厚度校正方法

技术领域

本发明属于石油勘探开发的三维地质建模领域，特别是涉及到一种基于三维地质建模软件的物性参数的厚度校正方法。

背景技术

三维地质建模是对地质体的数字化表述过程，主要应用于开发阶段的油藏静态描述，借助于地质模型可以更好地进行油田储层评价、开发方案编制与调整及开发效果预测等工作。属性参数建模是三维地质建模的重要组成部分，主要包括反映地层岩性、油气资源分布等特性的多种参数,如孔隙度、渗透率、含油饱和度等。这些属性参数均由测井曲线计算得到，而测井曲线受纵向分辨率的影响，在薄层的测井响应中包含了较多围岩的信号。因此，薄层的属性参数结果应该进行储层厚度的影响校正，储层的厚度属性是进行该项校正的基础。

Petrel软件由斯伦贝谢公司开发，是目前比较成熟的三维地质建模软件之一，利用Petrel软件可以建立起定量的三维可视化油藏模型。但在Petrel软件中并不能直接获得储层厚度属性，给薄层的属性参数校正带来不便。为此，我们发明了一种基于三维地质建模软件的物性参数的厚度校正方法，解决了以上技术问题。

发明内容

本发明的目的是提供一种能够方便快捷的物性参数的厚度校正方法。

本发明的目的可通过如下技术措施来实现：基于三维地质建模软件的物性参数的厚度校正方法，该基于三维地质建模软件的物性参数的厚度校正方法包括：步骤1，将储层厚度数据统一格式，并进行整数编码；步骤2，将厚度属性编码正确加载到三维地质建模软件中；步骤3，利用厚度属性编码计算得到准确的储层厚度属性曲线数据。步骤4，利用建立的厚度校正的解释模型，对孔隙度、渗透率等物性参数进行厚度校正计算。

本发明的目的还可通过如下技术措施来实现：

步骤1包括：

1)将储层厚度数据保留小数点后1位数字；

2)构造出1套离散厚度属性数据，其数值等于储层厚度值×10；

3)对离散数据的整数编码特征需要进行初始设置，根据储层的最大厚度确定储层编码的最大值Numer_Max：

Numer_Max＝max(T_i)

其中，Numer_Max表示所有小层厚度属性的最大值；T_i表示第i个小层的厚度值，i＝1,2,3,…,n；n为总小层数。

在步骤2中，在地质建模软件加载属性数据时，需要满足的格式为：

起始深度 属性代码

终止深度 非代码字符

准备的小层厚度属性表格式为：

井号 起始深度 终止深度 厚度属性代码

将小层数据表快速导出为建模软件可用的数据格式，具体实现过程为：

1)读取井号列数据，并作为新建文本的文件名；

2)读取同一井号的所有行，并将井号、起始深度、终止深度、厚度属性代码分别存储；

3)将起始深度和厚度属性代码写入第一行，终止深度和非代码字符写入第一行，并用空格隔开；

4)重复以上3个步骤，将所有井的厚度属性数据加载到三维地质建模软件中。

在步骤3中，从厚度属性编码计算得到准确的井点储层厚度属性曲线数据的计算公式为：

T_i＝T_i/10

其中，T_i表示第i个小层的厚度值，i＝1,2,3,…,n。

在步骤4中，厚度校正的解释模型的计算公式为：

PorC＝K×Por

其中，Por为厚度校正前的孔隙度值，PorC为厚度校正的孔隙度值，K为校正系数，含义是岩心声波孔隙度与测井解释声波孔隙度比值，其解释模型为：

K＝a*exp(-T/b)+c，

a，b，c为模型系数，由地区资料回归得到。

本发明由于采取以上技术方案，其具有以下优点：1、该方法以地质研究小层数据表成果为基础，储层厚度参数来源准确可靠。2、在具体实现时，利用Office软件的VBA工具，可以快速地将小层数据表转换为建模软件可用的数据格式，并正确加载到地质建模软件中。3、加载后，利用软件自带的计算器进行简单计算就可以得到准确的储层厚度参数，为后续的研究提供重要的储层参数。该方法方便快捷，技巧性强。利用本发明，准确地构建了厚度属性曲线，从而解决了三维地质建模中缺少储层厚度属性曲线的问题，并为利用厚度属性曲线进行的计算提供了数据。弥补了Petrel地质建模软件无法直接提供储层厚度属性的不足，为储层厚度的影响校正，进而精确求取孔隙度等储层属性提供了数据。

附图说明

图1为本发明的基于三维地质建模软件的物性参数的厚度校正方法的一具体实施例的流程图；

图2为本发明的一具体实施例中厚度数据成果图。

具体实施方式

为使本发明的上述和其他目的、特征和优点能更明显易懂，下文特举出较佳实施例，并配合附图所示，作详细说明如下。

如图1所示，图1为本发明的一种基于三维地质建模软件的物性参数的厚度校正方法的流程图。

在步骤101，具体实现过程为：

1)将储层厚度数据保留小数点后1位数字；

2)构造出1套离散厚度属性数据，其数值等于储层厚度值×10；

3)对离散数据的整数编码特征需要进行初始设置，根据储层的最大厚度确定储层编码的最大值Numer_Max：

Numer_Max＝max(T_i)

其中，Numer_Max表示所有小层厚度属性的最大值；T_i表示第i个小层的厚度值，i＝1,2,3,…,n；n为总小层数。

比如储层最大厚度为2.8米，则离散数据的属性编码增加至28个。流程进入到步骤103。

在步骤102，在地质建模软件加载属性数据时，一般需要满足的格式为：

起始深度 属性代码

终止深度 -9999(或其他非代码字符)

准备的小层厚度属性表格式为：

井号 起始深度 终止深度 厚度属性代码

因此，利用Office软件的VBA功能将小层数据表快速导出为建模软件可用的数据格式。其具体实现过程为：

1)读取井号列数据，并作为新建文本的文件名；

2)读取同一井号的所有行，并将井号、起始深度、终止深度、厚度属性代码分别存储；

3)将起始深度和厚度属性代码写入第一行，终止深度和-9999写入第一行，并用空格隔开；

4)重复以上3个步骤，将所有井的厚度属性数据加载到三维地质建模软件中。

实例井X1的小层数据表如表1。

表1实例井X1的小层数据表

转换后生成名为X1.txt的文档，储层厚度参数属性，其格式为：

上例中，实例X1井2914.3-2917.1米油层处加载的厚度属性为28(代表厚度2.8米×10)，将X1.txt加载到Petrel模型中。流程进入到步骤103。

在步骤103，从厚度属性编码计算得到准确的井点储层厚度属性曲线数据的计算公式为：

T_i＝T_i/10

其中，T_i表示第i个小层的厚度值，i＝1,2,3,…,n。

利用建模软件自带的计算器功能，将厚度属性离散数据/10，得到准确的储层厚度属性曲线。

实例井X1某干层起止深度为2907.3-2908.1米，厚度为0.8米，米油层处加载的厚度属性计算方法如下：8/10，计算得到结果为0.8，代表该储层厚度属性为7.3米。

步骤104，利用建立的厚度校正的解释模型，对孔隙度、渗透率等物性参数进行厚度校正计算。厚度校正的解释模型一般计算公式为：

PorC＝K×Por

其中，Por为厚度校正前的孔隙度值，PorC为厚度校正的孔隙度值，K为校正系数，含义是岩心声波孔隙度与测井解释声波孔隙度比值，其解释模型一般为：

K＝a*exp(-T/b)+c，

a，b，c为模型系数，由地区资料回归得到。

图2为本发明的一具体实施例中一种基于三维地质建模软件的物性参数的厚度校正方法成果示意图。研究区a＝-0.10775，b＝0.48648，c＝0.99914。图中左起第1道为深度道，第2道为自然伽马(GR)、自然电位(SP)曲线道，第3道为声波时差(AC)曲线道，第4道为4米梯度电阻率(R4)曲线道，第5道为感应电导率(COND)曲线道，第6道为解释结论道，第7道为储层厚度(T)曲线道，第8、9道分别为厚度校正前(POR0)、厚度校正后(PORC)孔隙度曲线道。

第1层干层起止深度为2907.3-2908.1米，厚度为0.8米，第2层干层起止深度为2909.4-2910.4米，厚度为1.0米，第1层油层起止深度为2914.3-2917.1米，厚度为2.8米。在第7道储层厚度曲线与储层厚度相符，厚度属性构建及加载正确。通过对厚度校正前、后孔隙度曲线的对比，薄层的孔隙度得到了校正。

上述各实施例仅用于说明本发明，凡是在本发明技术方案的基础上进行的等同变换和改进，均不应排除在本发明的保护范围之外。

Claims (1)

1.基于三维地质建模软件的物性参数的厚度校正方法，其特征在于，该基于三维地质建模软件的物性参数的厚度校正方法包括：

步骤1，将储层厚度数据统一格式，并进行整数编码；

步骤2，将厚度属性编码正确加载到三维地质建模软件中；

步骤3，利用厚度属性编码计算得到准确的储层厚度属性曲线数据；

步骤4，利用建立的厚度校正的解释模型，对孔隙度、渗透率等物性参数进行厚度校正计算；

步骤1包括：

1)将储层厚度数据保留小数点后1位数字；

2)构造出1套离散厚度属性数据，其数值等于储层厚度值×10；

3)对离散数据的整数编码特征需要进行初始设置，根据储层的最大厚度确定储层编码的最大值Numer_Max：

Numer_Max＝max(T_i)

其中，Numer_Max表示所有小层厚度属性的最大值；T_i表示第i个小层的厚度值，i＝1,2,3,…,n；n为总小层数；

在步骤2中，在地质建模软件加载属性数据时，需要满足的格式为：

起始深度属性代码；

终止深度非代码字符；

准备的小层厚度属性表格式为：

井号起始深度，终止深度厚度属性代码；

将小层数据表快速导出为建模软件可用的数据格式，具体实现过程为：

1)读取井号列数据，并作为新建文本的文件名；

2)读取同一井号的所有行，并将井号、起始深度、终止深度、厚度属性代码分别存储；

3)将起始深度和厚度属性代码写入第一行，终止深度和非代码字符写入第一行，并用空格隔开；

4)重复以上3个步骤，将所有井的厚度属性数据加载到三维地质建模软件中；

在步骤3中，从厚度属性编码计算得到准确的井点储层厚度属性曲线数据的计算公式为：T_i＝T_i/10；

其中，T_i表示第i个小层的厚度值，i＝1,2,3,…,n；

在步骤4中，厚度校正的解释模型的计算公式为：

PorC＝K×Por

其中，Por为厚度校正前的孔隙度值，PorC为厚度校正的孔隙度值，K为校正系数，含义是岩心声波孔隙度与测井解释声波孔隙度比值，其解释模型为：

K＝a*exp(-T/b)+c；

a，b，c为模型系数，由地区资料回归得到。

Images