[发明专利]一种基于测井多参数识别低渗透储层成岩相的方法

摘要

本发明公开了一种基于测井多参数识别低渗透储层成岩相的方法，对研究区目的层取芯井段岩心进行观察分析，对铸体薄片、普通薄片、扫描电镜等资料进行研究，根据研究区目的层段内控制储层物性的胶结物类型、产状和成岩作用类型划分成岩相类型；选取标准成岩相样本，初步选择成岩相划分的测井参数，利用极值化方法，对成岩相划分的测井参数进行处理，并利用Wilks统计量方法，对能够表征成岩相类型的测井参数进行筛选；利用多元统计理论，建立成岩相的测井分类识别体系，实现低渗透储层的成岩相的测井识别。本发明的有益效果是能够利用测井资料对鄂尔多斯盆地西南缘演武地区低渗透储层的成岩相进行准确识别。

主权项

1.一种基于测井多参数识别低渗透储层成岩相的方法，其特征在于按照以下步骤进行：步骤1：对研究区目的层段取芯井段岩心进行观察分析，对铸体薄片、普通薄片和扫描电镜资料进行研究，根据研究区目的层段内控制储层物性的胶结物类型、产状和成岩作用类型划分成岩相类型；步骤2：选取标准成岩相样本，并根据研究区目的层段内测井成岩相的响应特征，初步选择成岩相划分的测井参数，利用极值化方法，对成岩相划分的测井参数进行处理，并利用Wilks统计量方法，对能够表征成岩相类型的测井参数进行筛选；步骤3：利用多元统计理论，建立成岩相的测井分类识别体系，实现低渗透储层的成岩相的测井识别；具体步骤如下：(1)根据岩心分析、物性测试、微观试验资料，将研究区目的层段内的成岩相根据控制储层物性的胶结物类型、产状和成岩作用类型分为n类；(2)根据成岩相的测井响应特征，初步选取对成岩相比较敏感的g个测井参数，参数包括自然伽马GR、岩性密度DEN、电阻率RT、补偿中子CNL、自然电位SP、补偿声波AC，中子‑密度视石灰岩孔隙度差Ф_ND参数；(3)对成岩相测井参数极值化，并求取能够反映成岩相层段内测井形态的参数S，即：v_j为v_i测井参数第j个采样点的极差转化值，N为选定沉积微相段内的采样点数，GR极差转化均值v₁，CNL极差转化均值v₂，DEN极差转化均值v₃，RT极差转化均值v₄，Ф_ND极差转化均值v₅，SP极值转化均值v₆，AC极值转化均值v₇，S为GR极差转化离差平方和；(4)利用Wilks准则，初选对成岩相比较敏感的g个测井参数、引入参数中子‑密度视石灰岩孔隙度差Ф_ND及能够反映成岩相层段内测井形态的参数S做出筛选，对判别函数无关或影响较小的变量进行剔除，构建合理的判别方程，其步骤为：先将样品从P个总体b₁、b₂、…、b_p中分别取出n₁、n₂、…、n_p样品，每个样品有g个测井参数，样品构成的观测样本组成g维矩阵，则其成分变量为：各变量的累计概率符合正态分布，则建立类内离差矩阵为：建立总离差矩阵为：统计量U为：U＝|W|/|T|   (6)将已带入的h个变量中，检验变量x^(r)是否应该被剔除，则统计量U_r此时为其判别方程为：若F值小于检验水平的临界值F′_a，就将x^(r)剔除，否则保留，最终提取m个参数作为成岩相的测井判别参数；(5)判别函数的建立，将所得p个样品，每个样品提取m个参数极值化后，总体样本的平均值为：式中为样品的变量平均值；p＝1,2,…,P；P为样品个数；i＝1,2，…,m；利用公式求取样本的协方差矩阵，公式为：i，j＝1，2，…，m根据多元统计理论，采用判别函数：式中：所得判别数值模型为：其中，F₁、F₂…F_n为权重函数；x₁，x₂…x_m为选取的m个测井参数变量；a₀₀，a₀₁…a_0m、a₁₀，a₁₁…a_1m、a_(n‑1)0，a_(n‑1)1…a_(n‑1)m分别为权重函数的判别系数；(6)同类型成岩相具有不同的判别函数值F_n(x)n＝1,2,…,p，将新的未知样本参数带入判别函数，达到最大值即为第n^*类成岩相，即：样本属于n^*类成岩相的条件概率为：式中：F_n(x)为待判定成岩相样本的参数带入各个成岩相判别公式后的函数值；为成岩相判别函数中的最大值；n为成岩相的分类数；n^*为函数值大的测井相的类型；Q_n为待判定样本属于n^*类成岩相的条件概率。