[发明专利]一种全井段岩性识别方法

说明书

技术领域

本发明涉及一种岩性识别方法，特别是涉及一种采用交互式多维直方图对岩性进行全井段岩性类型定量识别的方法。

背景技术

目前利用测井信息识别岩性的方法主要有常规测井资料解释法、判别分析法、神经网络算法等。常规测井资料解释法目的主要是确定地层泥质含量，解释结果仅为砂岩和泥岩两种岩性,解释结果太粗,无法反映地质非均质性情况；利用判别分析法识别地层岩性较为普遍，但对原始样本数据的代表性和影响岩性类别的指标的全面性要求较高，针对性较差，由于受到流体性质等因素的影响，可靠性也较差，无法满足储层物性的精细评价和沉积学研究的需要；神经网络算法，该方法泛化能力不强，识别结果准确性较差，且需要人工先验知识，在岩心资料较少的情况下，很难用常规方法进行岩性识别。上述三种方法岩性识别精度不够，均无法满足测井精细评价和地质研究的需要。

发明内容

本发明提供一种全井段岩性识别方法，其目的是提供一种利用岩心资料和测井资料通过数理统计学建立多维直方图模型进行岩性的分类和定量识别，最终建立精细连续的岩性剖面。

为了实现上述发明目的，本发明采取的技术方案如下：

一种全井段岩性识别方法，包括：

从多个深度位置取出多个岩心样本，确定每一岩心样本的粒度数据及其岩性类型；

确定所述多个深度位置处的多种测井参数的值，选择随着岩性样本粒度增大而增大或减小的N种测井参数作为特征参数，N≥2；

根据所述多个深度位置处岩心样本点的岩性类型及所述N种特征参数的值，确定以所述特征参数坐标的N维空间中，各N维单元格与岩性类型之间的对应关系；

根据全井段各深度位置处所述特征参数的值及所述对应关系，确定各深度位置处的代表岩性类型。

进一步地，在所述N维空间中确定所述特征参数与岩性类型之间的对应关系的方法为：

对所述多个深度位置中的每一深度位置，将该深度位置处的岩性类型确定为该深度位置处N种特征参数值所指向的N维单元格对应的一份岩性类型，由此得到所述N维空间中，一个或多个N维单元格对应的一种或多种岩性类型及每一种岩性类型的份数；

对所述N维空间中没有对应的岩性类型的N维单元格，通过就近插值法插值生成对应的一种或多种岩性类型及其份数。

进一步地，各深度位置处的代表岩性类型为此深度位置所有岩性类型中份数最多的岩性类型。

进一步地，所述岩心样本点粒度数据容量不少于五个。

进一步地，所述测井参数为自然伽马参数、密度参数、声波参数和中子参数。

另一种可替代的技术方案为：

一种交互式多维直方图全井段岩性识别方法，包括：

步骤一，从多个深度位置取出多个岩心样本，确定每一岩心样本的粒度数据，并根据岩心样本粒度数据确定岩性类型；

步骤二，用同一深度位置处的岩性类型刻度多种测井参数的值，并建立测井参数与岩性类型之间的交会图；

步骤三，选择交会图中，随着岩性样本粒度增大而增大或减小的测井参数作为特征参数，分别统计每一个特征参数中各类岩性类型出现的份数，并建立多个特征参数-岩性类型-份数之间的直方图模型；

步骤四，将步骤三中得到所有直方图模型统一表示在一个N维坐标空间中，N≥2；

步骤五，对于N维坐标空间的单元格中没有岩性类型的，查找距拟生成单个格空间距离最近的样本点单元格，采用就近插值法插值生成无样本点单元格中的各类岩性类型，并更新N维坐标空间，最终在N维坐标空间形成包括全井段岩性类型分布图；

步骤六，分别统计N维坐标空间每个单元格中各类岩性类型的数目，并计算各类岩性类型在每个单元格内出现的份数；

步骤七，将全井段的每一个深度区段所对应的特征参数值逐个在N维坐标空间中单元格寻找对应位置，并以此单元格中出现份数最多的岩性作为对应深度区段的岩性类型代表，然后对连续且岩性类型相同的特征参数值合并，形成全井段的连续岩性剖面。

优选地，所述步骤一中，岩心样本点粒度数据容量不少于五个。

优选地，所述步骤二中，测井参数为自然伽马参数、密度参数、声波参数或中子参数。

优选地，所述步骤三，直方图中每个特征参数轴根据其所代表的特征参数值范围等间距划分，形成空间单元格。

优选地，所述步骤四，N维坐标空间中每个维度间距和数值范围的确定同其所代表特征参数在直方图中范围和间距。

本发明和现有技术相比，具有如下有益效果：