本发明专利技术公开了一种识别高盐富钾区注水开发油藏碎屑岩岩性的方法，统计已知岩性的岩心分析化验资料并将岩性分类，选取对岩性反映比较敏感的n个测井参数；基于数理统计原理，建立不同岩性识别方程，再将单井未取心段测井数据体分别代入各类岩性识别方程，根据识别准则，识别方程中函数值最大者则为相应岩性的概率最高，该岩性则确定为相应的岩性。通过拟合综合判别曲线识别高盐富钾区注水开发油藏碎屑岩岩性，实现对多类岩性的准确判别，可建立连续的地层岩性剖面，自检符合率达90％，识别效果较好，对类似地区岩性识别具有一定的借鉴意义。
【技术实现步骤摘要】
本专利技术属于储层评价
，涉及一种识别高盐富钾区注水开发油藏碎肩岩岩 性的方法。
技术介绍
测井岩性识别是准确求取储层参数，进行储层评价工作的基础和关键。注水开发 油藏碎肩岩岩性一般通过自然电位(SP)、自然伽马(GR)、声波时差(AC)、密度(DEN)和电阻 率(RT)等测井资料进行识别。 常规的识别方法主要包括交会图法、神经网络法和曲线重叠法等，但在高盐富钾 地区存在识别率低，不能连续识别地层岩性的缺陷。滩海地区地质条件复杂，且具有高盐富 钾的特点，盐水泥浆基液使得自然电位反向出现正异常，电阻率值比正常情况下降低从而 无法真实反映原状地层特性;地层中存在高钾矿物，又使得自然伽马曲线识别精度降低。因 此单条测井曲线不能有效反映碎肩岩岩性特征，常规方法也无法准确区分不同类岩性。 为了更好地识别高盐富钾区注水开发油藏碎肩岩岩性并开展后续储层评价工作， 放大和提取各类测井参数的差异，综合多参数多资料进行岩性识别就变得非常有意义。
技术实现思路
本专利技术的目的是提供一种识别高盐富钾区注水开发油藏碎肩岩岩性的方法，基于 数理统计原理，建立岩性识别方程，构建差值识别函数，拟合岩性识别综合曲线，达到准确 识别高盐富钾区注水开发油藏碎肩岩岩性的目的。 本专利技术所采用的技术方案是，一种识别高盐富钾区注水开发油藏碎肩岩岩性的方 法，按照以下步骤进行： 步骤1，统计已知岩性的岩心分析化验资料并将岩性分类;采用数理统计方法，剔 除非相关变量，选取对岩性反映比较敏感的η个测井参数； 步骤2,基于数理统计原理，建立岩性识别权重函数F，计算出的权重函数根据η个 测井参数对岩性的贡献率赋予不同的函数特征值，提取与各组有最大可能多重相关的、函 数特征值最大的两个权重函数，即第一权重函数和第二权重函数： Fi = a〇〇+a〇i · χιΗ-----i~a〇n · xn (1) F2 = ai〇+aii · xi+."+ain · xn (2) 其中，Fi、F2为第一和第二权重函数；χ^χπη为选取的n个测井参数变量;a00， a〇i…a〇n和ai〇，a2〇…ain分别对应两个权重函数的判别系数； 步骤3,将给出的第一权重函数和第二权重函数用于计算各样品在降维后空间坐 标系中的坐标值，然后计算出各类别在低维空间中的中心坐标，根据各样品点距离各类别 重心的距离构造出每种岩性的识别函数，即： fi(X) =bi〇+bii · xi+---+bin · xn (3)其中，h为第i类岩性识别函数;b^bn···!^为第i类岩性识别函数的判别系数； 步骤4,对于两种岩性的判别，将砂岩岩性识别函数减去泥岩岩性识别函数，得出 差值识别函数，该函数即为两种岩性的综合识别曲线，记为P(x);即： p⑴=A f = f砂-f泥=Δ bo+Δ bi · χι+···+Δ bn · Xn (4) 其中，P(x)为两种岩性的综合识别曲线;f砂为砂岩岩性识别函数，f泥为泥岩岩性识 另IJ函数;Χ1，χ 2···χΑ选取的η个测井参数变量；△ bo，Δ br·· Δ bn为两种岩性综合识别曲线的 判别系数； 把待判测井数据体的观测值Χ=(χ(1)χ(2)···χ(η))代入上述差值识别函数计算出岩 性综合识别曲线值Poo，若Poo 2 〇则为砂岩岩性，相反综合曲线值Poo < 〇则为泥岩岩性； 步骤5,对于Μ类岩性的判别，把样品的观测值Χ=(χ (1)χ(2)···χ(η))代入式(3)中，分 别求得fi( x)，f2(x)"_fM(x)，取最大值记为fg( x)=max{fi(x)}，1 < i SM;由判别准则知，岩性识别 函数值大为相应岩性的概率高，则认为待判测井数据体XefgW对应的岩性所属类型，同时 岩性综合识别曲线Q(x)赋值为m;对Μ类不同岩性识别结果赋值，从而拟合出Μ类岩性的综合 识别曲线Q(x);即：其中，Quo为Μ类岩性的综合识别曲线;f1(x)，f2(x)…fM(x)为Μ类岩性识别函数;f g(x) 为代入待判数据体后Μ类岩性识别函数的最大值;?Μ,πινΠΗ为Μ类岩性综合识别曲线赋值， nu的赋值区间为-Μ < mi < Μ。进一步的，所述步骤1中，η个测井参数包括自然电位SP、自然伽马GR、声波时差AC、 岩性密度DEN、电阻率RT。 本专利技术的有益效果是克服常规岩性识别方法的不足，提供一种识别高盐富钾区注 水开发油藏碎肩岩岩性的方法，即充分利用测井、岩心分析化验等资料，以数理统计原理为 基础，建立岩性识别函数，并拟合单井岩性识别综合曲线，利用该曲线识别岩性效果较好， 自检符合率高达90%。本专利技术能很好地考虑多种测井参数，进而拟合出一条准确的岩性综 合识别曲线，计算结果合理且准确率高，实现对高盐富钾区注水开发油藏全井段碎肩岩岩 性的连续判别，可为类似地区的岩性识别提供借鉴。【附图说明】 图1是本专利技术实施例中研究区两种岩性的综合识别曲线。图2是本专利技术实施例中研究区四种岩性的综合识别曲线。【具体实施方式】 为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对 本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本专利技术，并 不用于限定本专利技术。 -种识别高盐富钾区注水开发油藏碎肩岩岩性的方法，具体按照以下步骤进行： 步骤1，统计已知岩性的岩心分析化验资料并将岩性分类。采用数理统计方法，剔 除非相关变量，选取对岩性反映比较敏感的η个测井参数（自然电位(SP)、自然伽马(GR)、声 波时差(AC)、岩性密度(DEN)、电阻率(RT)…）；步骤2，基于数理统计原理，建立岩性识别权重函数F(包括Fi、F2、Fr··)，计算出的 权重函数可以根据η个测井参数对岩性的贡献率赋予不同的函数特征值，提取与各组有最 大可能多重相关的、函数特征值最大的两个权重函数，即第一权重函数和第二权重函数： Fi = aoo+aoi · xi+---+a〇n · xn (1) F2 = aio+aii · xi+."+ain · Xn (2) 其中，Fi、F2为第一和第二权重函数；X1，x 2…xn*选取的n个测井参数变量；aoo， aoi…a〇n和aio，a2〇…ain分别对应两个权重函数的判别系数。 步骤3,将给出的第一权重函数和第二权重函数用于计算各样品在降维后空间坐 标系中的坐标值，然后计算出各类别在低维空间中的中心坐标，根据各样品点距离各类别 重心的距离构造出每种岩性的识别函数，即： fi(X) =bi〇+bii · χι+...+bin · Xn (3) 其中，fi为第i类岩性识别函数;幻，^_心为选取的11个测井参数变量;1^，1^1···!^ 为第i类岩性识别函数的判别系数。 步骤4，对于两种岩性的判别，将砂岩岩性识别函数减去泥岩岩性识别函数，得 出差值识别函数，该函数即为两种岩性的综合识别曲线，记为P〇〇;即： p⑴=A f = f砂-f泥=Δ bo+Δ bi · χι+···+Δ bn · Xn (4) 其中，P(x)为两种岩性的综合识别曲线;f砂为砂岩岩性识别函数，f泥为泥岩岩性识 另IJ函数;Χ1，χ 2···χΑ选取的η个测井参数变量；△ bo，Δ br·· Δ b本文档来自技高网...
【技术保护点】
一种识别高盐富钾区注水开发油藏碎屑岩岩性的方法，其特征在于，按照以下步骤进行：步骤1，统计已知岩性的岩心分析化验资料并将岩性分类；采用数理统计方法，剔除非相关变量，选取对岩性反映比较敏感的n个测井参数；步骤2，基于数理统计原理，建立岩性识别权重函数F，计算出的权重函数根据n个测井参数对岩性的贡献率赋予不同的函数特征值，提取与各组有最大可能多重相关的、函数特征值最大的两个权重函数，即第一权重函数和第二权重函数：F1＝a00+a01·x1+…+a0n·xn (1)F2＝a10+a11·x1+…+a1n·xn (2)其中，F1、F2为第一和第二权重函数；x1，x2…xn为选取的n个测井参数变量；a00，a01…a0n和a10，a20…a1n分别对应两个权重函数的判别系数；步骤3，将给出的第一权重函数和第二权重函数用于计算各样品在降维后空间坐标系中的坐标值，然后计算出各类别在低维空间中的中心坐标，根据各样品点距离各类别重心的距离构造出每种岩性的识别函数，即：fi(X)＝bi0+bi1·x1+…+bin·xn (3)其中，fi为第i类岩性识别函数；bi0，bi1…bin为第i类岩性识别函数的判别系数；步骤4，对于两种岩性的判别，将砂岩岩性识别函数减去泥岩岩性识别函数，得出差值识别函数，该函数即为两种岩性的综合识别曲线，记为P(X)；即：P(X)＝Δf＝f砂‑f泥＝Δb0+Δb1·x1+…+Δbn·xn (4)其中，P(X)为两种岩性的综合识别曲线；f砂为砂岩岩性识别函数，f泥为泥岩岩性识别函数；x1，x2…xn为选取的n个测井参数变量；Δb0，Δb1…Δbn为两种岩性综合识别曲线的判别系数；把待判测井数据体的观测值X＝(x(1)x(2)…x(n))代入上述差值识别函数计算出岩性综合识别曲线值P(X)，若P(X)≥0则为砂岩岩性，相反综合曲线值P(X)≤0则为泥岩岩性；步骤5，对于M类岩性的判别，把样品的观测值X＝(x(1)x(2)…x(n))代入式(3)中，分别求得f1(x),f2(x)…fM(x)，取最大值记为fg(x)＝max{fi(x)}，1≤i≤M；由判别准则知，岩性识别函数值大为相应岩性的概率高，则认为待判测井数据体X∈fg(x)对应的岩性所属类型，同时岩性综合识别曲线Q(X)赋值为m；对M类不同岩性识别结果赋值，从而拟合出M类岩性的综合识别曲线Q(X)；即：Q(X)=m1····················fg(x)=f1(x)m2···················fg(x)=f2(x)m3····················fg(x)=f3(x)·....................................··.....................................··.....................................·mi·······&Cen...
【技术特征摘要】
【专利技术属性】
技术研发人员：杨少春，牛海瑞，温雅茹，
申请(专利权)人：中国石油大学华东，
类型：发明
国别省市：山东;37
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。