一种识别高盐富钾区注水开发油藏碎屑岩岩性的方法
【技术领域】
[0001]本发明属于储层评价技术领域，涉及一种识别高盐富钾区注水开发油藏碎肩岩岩 性的方法。
【背景技术】
[0002] 测井岩性识别是准确求取储层参数，进行储层评价工作的基础和关键。注水开发 油藏碎肩岩岩性一般通过自然电位(SP)、自然伽马(GR)、声波时差(AC)、密度(DEN)和电阻 率(RT)等测井资料进行识别。
[0003] 常规的识别方法主要包括交会图法、神经网络法和曲线重叠法等，但在高盐富钾 地区存在识别率低，不能连续识别地层岩性的缺陷。滩海地区地质条件复杂，且具有高盐富 钾的特点，盐水泥浆基液使得自然电位反向出现正异常，电阻率值比正常情况下降低从而 无法真实反映原状地层特性;地层中存在高钾矿物，又使得自然伽马曲线识别精度降低。因 此单条测井曲线不能有效反映碎肩岩岩性特征，常规方法也无法准确区分不同类岩性。
[0004] 为了更好地识别高盐富钾区注水开发油藏碎肩岩岩性并开展后续储层评价工作， 放大和提取各类测井参数的差异，综合多参数多资料进行岩性识别就变得非常有意义。
【发明内容】
[0005] 本发明的目的是提供一种识别高盐富钾区注水开发油藏碎肩岩岩性的方法，基于 数理统计原理，建立岩性识别方程，构建差值识别函数，拟合岩性识别综合曲线，达到准确 识别高盐富钾区注水开发油藏碎肩岩岩性的目的。
[0006] 本发明所采用的技术方案是，一种识别高盐富钾区注水开发油藏碎肩岩岩性的方 法，按照以下步骤进行：
[0007] 步骤1，统计已知岩性的岩心分析化验资料并将岩性分类;采用数理统计方法，剔 除非相关变量，选取对岩性反映比较敏感的η个测井参数；
[0008] 步骤2,基于数理统计原理，建立岩性识别权重函数F，计算出的权重函数根据η个 测井参数对岩性的贡献率赋予不同的函数特征值，提取与各组有最大可能多重相关的、函 数特征值最大的两个权重函数，即第一权重函数和第二权重函数：
[0009] Fi = a〇〇+a〇i · χιΗ-----i~a〇n · xn (1)
[0010] F2 = ai〇+aii · xi+."+ain · xn (2)
[0011] 其中，Fi、F2为第一和第二权重函数；χ^χπη为选取的n个测井参数变量;a00， a〇i…a〇n和ai〇，a2〇…ain分别对应两个权重函数的判别系数；
[0012] 步骤3,将给出的第一权重函数和第二权重函数用于计算各样品在降维后空间坐 标系中的坐标值，然后计算出各类别在低维空间中的中心坐标，根据各样品点距离各类别 重心的距离构造出每种岩性的识别函数，即：
[0013] fi(X) =bi〇+bii · xi+---+bin · xn (3)
[0014]其中，h为第i类岩性识别函数;b^bn···!^为第i类岩性识别函数的判别系数；
[0015] 步骤4,对于两种岩性的判别，将砂岩岩性识别函数减去泥岩岩性识别函数，得出 差值识别函数，该函数即为两种岩性的综合识别曲线，记为P(x);即：
[0016] p⑴=A f = f砂-f泥=Δ bo+Δ bi · χι+···+Δ bn · Xn (4)
[0017] 其中，P(x)为两种岩性的综合识别曲线;f砂为砂岩岩性识别函数，f泥为泥岩岩性识 另IJ函数;Χ1，χ 2···χΑ选取的η个测井参数变量；△ bo，Δ br·· Δ bn为两种岩性综合识别曲线的 判别系数；
[0018] 把待判测井数据体的观测值Χ=(χ(1)χ(2)···χ(η))代入上述差值识别函数计算出岩 性综合识别曲线值Poo，若Poo 2 〇则为砂岩岩性，相反综合曲线值Poo < 〇则为泥岩岩性； [0019]步骤5,对于Μ类岩性的判别，把样品的观测值Χ=(χ (1)χ(2)···χ(η))代入式(3)中，分 别求得fi( x)，f2(x)"_fM(x)，取最大值记为fg( x)=max{fi(x)}，1 < i SM;由判别准则知，岩性识别 函数值大为相应岩性的概率高，则认为待判测井数据体XefgW对应的岩性所属类型，同时 岩性综合识别曲线Q(x)赋值为m;对Μ类不同岩性识别结果赋值，从而拟合出Μ类岩性的综合 识别曲线Q(x);即：
[0021]其中，Quo为Μ类岩性的综合识别曲线;f1(x)，f2(x)…fM(x)为Μ类岩性识别函数;f g(x) 为代入待判数据体后Μ类岩性识别函数的最大值;?Μ,πινΠΗ为Μ类岩性综合识别曲线赋值， nu的赋值区间为-Μ < mi < Μ。
[0022]进一步的，所述步骤1中，η个测井参数包括自然电位SP、自然伽马GR、声波时差AC、 岩性密度DEN、电阻率RT。
[0023] 本发明的有益效果是克服常规岩性识别方法的不足，提供一种识别高盐富钾区注 水开发油藏碎肩岩岩性的方法，即充分利用测井、岩心分析化验等资料，以数理统计原理为 基础，建立岩性识别函数，并拟合单井岩性识别综合曲线，利用该曲线识别岩性效果较好， 自检符合率高达90%。本发明能很好地考虑多种测井参数，进而拟合出一条准确的岩性综 合识别曲线，计算结果合理且准确率高，实现对高盐富钾区注水开发油藏全井段碎肩岩岩 性的连续判别，可为类似地区的岩性识别提供借鉴。
【附图说明】
[0024] 图1是本发明实施例中研究区两种岩性的综合识别曲线。
[0025]图2是本发明实施例中研究区四种岩性的综合识别曲线。
【具体实施方式】
[0026] 为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对 本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并 不用于限定本发明。
[0027] -种识别高盐富钾区注水开发油藏碎肩岩岩性的方法，具体按照以下步骤进行： [0028]步骤1，统计已知岩性的岩心分析化验资料并将岩性分类。采用数理统计方法，剔 除非相关变量，选取对岩性反映比较敏感的η个测井参数（自然电位(SP)、自然伽马(GR)、声 波时差(AC)、岩性密度(DEN)、电阻率(RT)…）；
[0029]步骤2，基于数理统计原理，建立岩性识别权重函数F(包括Fi、F2、Fr··)，计算出的 权重函数可以根据η个测井参数对岩性的贡献率赋予不同的函数特征值，提取与各组有最 大可能多重相关的、函数特征值最大的两个权重函数，即第一权重函数和第二权重函数： [0030] Fi = aoo+aoi · xi+---+a〇n · xn (1)
[0031] F2 = aio+aii · xi+."+ain · Xn (2)
[0032] 其中，Fi、F2为第一和第二权重函数；X1，x 2…xn*选取的n个测井参数变量；aoo， aoi…a〇n和aio，a2〇…ain分别对应两个权重函数的判别系数。
[0033] 步骤3,将给出的第一权重函数和第二权重函数用于计算各样品在降维后空间坐 标系中的坐标值，然后计算出各类别在低维空间中的中心坐标，根据各样品点距离各类别 重心的距离构造出每种岩性的识别函数，即：
[0034] fi(X) =bi〇+bii · χι+...+bin · Xn (3)
[0035] 其中，fi为第i类岩性识别函数;幻，^_心为选取的11个测井参数变量;1^，1^1···!^ 为第i类岩性识别函数的判别系数。
[0036] 步骤4，对于两种岩性的判别，将砂岩岩性识别函数减去泥岩岩性识别函数，得
[0037] 出差值识别函数，该函数即为两种岩性的综合识别曲线，记为P〇〇;即：
[0038] p⑴=A f = f砂-f泥=Δ bo+Δ bi · χι+···+Δ bn · Xn (4)
[0039] 其中，P(x)为两种岩性的综合识别曲线;f砂为砂岩岩性识别函数，f泥为泥岩岩性识 另IJ函数;Χ1，χ 2···χΑ选取的η个测井参数变量；△ bo，Δ br·· Δ bn为两种岩性综合识别曲线的 判别系数。
[0040] 把待判测井数据体的观测值X=(x(1)x(2)…x(n))代入上述差值识