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

权利要求书

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^*类成岩相的条件概率。