一种测井数据岩性层数值归约方法技术

本发明专利技术涉及一种测井数据岩性层数值归约方法，包括以下步骤：1)采用五点二次法对原始测井数据进行平滑处理，得到滤波后的平滑测井数据；2)采用活度分层法对平滑测井数据进行岩性分层，计算平滑测井数据的活度值，并根据活度截止值得到各个地层的分界点，相邻分界点形成一个岩性地层，进而得到岩性地层分层数据；3)根据岩性地层分层数据对原始测井数据进行分段，在每一段原始测井数据上，根据测井数据形态和测井仪器的纵向分辨率，结合岩性分层数据，提取各岩性地层的特征值数据；4)对所有岩性地层的特征值进行数值归约，将特征值相同的相邻岩性层归为一类，生成测井数据岩性层数值归约数据表。本发明专利技术可以广泛应用于油气勘探、开发测井储层评价、精细油藏描述等领域。

【技术实现步骤摘要】

本专利技术涉及一种石油测井曲线数据预处理方法，特别是关于一种测井数据岩性层数值归约方法。
技术介绍
石油测井数据岩性层数值归约(也称测井数据岩性层分层取值)是石油测井储层评价和储量计算不可缺少的重要环节。复杂的数据分析和数据挖掘技术也要求对海量石油测井数据进行归约表示。经过归约处理后的石油测井数据集要比原有的石油测井数据集小很多，而且能够基本保持原有石油测井数据的完整性。这样就使得原来针对海量石油测井数据无法实现的数据处理得以进行，并且产生相同或几乎相同的分析结果。目前，测井数据岩性层数值归约方法包括两种，一种是人工分层取值，另一种是自动分层取值。其中人工分层取值通常采用的方法为首先确定取值深度，然后从深度曲线上人工读值，读值时根据人对取值深度处相邻深度曲线值大小的综合认识进行的，当人的认识不同时读取的具体数值也会不同，该方法劳动强度大，且读值依赖于人的经验，无法保证取值精度。自动分层取值分为自动分层和自动取值两步，其中，自动分层方法有微分法、活度法、句法分析法和小波分析法等；自动取值方法大多是基于曲线形态进行取值，但是其并没有考虑测井方法本身的技术参数，例如探头的源距、间距等。
技术实现思路
针对上述问题，本专利技术的目的是提供一种充分考虑测井方法本身技术参数，效率高、取值精度高的测井数据岩性层数值归约方法。为实现上述目的，本专利技术采取以下技术方案：一种测井数据岩性层数值归约方法，包括以下步骤：1)采用五点二次法对原始测井数据进行平滑处理，得到滤波后的平滑测井数据；2)采用活度分层法对平滑测井数据进行岩性分层，计算所述平滑测井数据的活度值，并根据活度截止值得到各个地层的分界点，相邻分界点形成一个岩性地层，进而得到岩性地层分层数据；3)根据岩性地层分层数据对原始测井数据进行分段，在每一段原始测井数据上，根据测井数据形态和测井仪器的纵向分辨率，结合岩性分层数据，提取各岩性地层的特征值数据；4)对所有岩性地层的特征值进行数值归约，将特征值相同的相邻岩性层归为一类，生成测井数据岩性层数值归约数据表。所述步骤2)中，所述平滑测井数据中第i个滑动平均值的活度值Ei为：Ei=Σk=i-N/2i+N/2[yk-y‾i]2,]]>其中，N为窗长L内平滑测井数据的存取样点数，L为给定活度分层窗长，且L∈N+，为滑动平均值yi前、后各N/2范围内所有滑动平均值的平均值，即为：y‾i=1NΣk=i-(N/2)i+(N/2)yk.]]>所述步骤3)中，提取各岩性地层特征值数据的规则为：①当h<＝d时，取该岩性地层测井数据的极值点作为特征值，其中h表示所述步骤2)中得到的各岩性地层的厚度，d表示测井仪器的纵向分辨率；②当h>d时，根据该岩性地层测井数据的形态采用不同的取值规则，测井数据形态分为三类：a、第一类：该岩性地层测井数据的形态特征为其只有一个极值点，此时取该极值点的值作为该岩性地层的特征值；b、第二类：该岩性地层测井数据的形态特征为其第一个极值点及最后一个极值点同为极大或极小值点，且第一个极值点与最后一个极值点差值大于该层测井数据的最大值与最小值的差值的1/3，此时取第一个极值点与最后一个极值点中较大或较小的值作为该岩性地层的特征值；c、第三类：该岩性地层测井数据的形态特征为第一个极值点与最后一个极值点同为极大或极小值，或者呈现相反的极性，且各极值点在一个值上下摆动，也就是各个极值点的两两差值都小于该段测井数据的最大值与最小值差值的1/3；则取第一个极值点与最后一个极值点之间所有极值点的平均值作为该岩性地层的特征值。本专利技术由于采取以上技术方案，其具有以下优点：1、本专利技术由于采用五点二次法平滑滤波和活度分层相结合的方法对测井曲线进行自动分层，获得了效率高、满足精度要求的分层取值效果。2、本专利技术由于针对不同的测井仪器和不同的测井曲线特征采用不同的取值原则，有效的保证了对测井曲线分层取值的精度，能够得到真实的地层测井值。本专利技术由于将将滤波、分层、取值三个技术环节配套应用，可以获得理想的岩性层测井值，因而本专利技术可以广泛应用于油气勘探、开发测井储层评价、精细油藏描述等领域。附图说明图1是本专利技术方法流程示意图，图2是本专利技术涉及到的三类测井曲线形态，其中横坐标是测井数据取值，纵坐标是测井数据的深度值，竖线为判断极值点差值的基准线。具体实施方式下面结合附图和实施例对本专利技术进行详细的描述。如图1所示，本专利技术测井数据岩性层数值归约方法包括以下步骤：1)采用五点二次法对原始测井数据进行平滑处理，得到滤波后的平滑测井数据。对原始测井数据进行采样，得到测井值xi，其中i＝1,2,...,n，n为总采样点数；则第i个采样点测井值xi的滑动平均值yi为：yi=135[-3(xi-2+xi+2)+12(xi-1+xi+1)+17xi],(i=3,4,...,n-3,n-2)---(1)]]>2)采用活度分层法对平滑测井数据进行岩性分层，根据活度截止值得到各个地层的分界点，相邻分界点形成一个岩性地层，进而得到岩性地层分层数据。首先根据平滑测井数据计算活度值，第i个滑动平均值yi的活度值Ei为：Ei=Σk=i-N/2i+N/2[yk-y‾i]2---(2)]]>其中，N为窗长L内平滑测井数据的存取样点数，L为给定活度分层窗长，且L∈N+，为滑动平均值yi前、后各N/2范围内所有滑动平均值的平均值，即为：y‾i=1NΣk=i-(N/2)i+(N/2)yk---(3)]]>根据实际应用中对分层精度的要求，分段设置活度截止值ζ(ζ>＝0)，根据活度截止值对离散测井曲线进行分层，得到各个地层之间的分界点集合O为：O＝{i|Ei>Ei-1且Ei>Ei+1且Ei>＝ζ,i＝1,.....n本文档来自技高网...

【技术保护点】
一种测井数据岩性层数值归约方法，包括以下步骤：1)采用五点二次法对原始测井数据进行平滑处理，得到滤波后的平滑测井数据；2)采用活度分层法对平滑测井数据进行岩性分层，计算所述平滑测井数据的活度值，并根据活度截止值得到各个地层的分界点，相邻分界点形成一个岩性地层，进而得到岩性地层分层数据；3)根据岩性地层分层数据对原始测井数据进行分段，在每一段原始测井数据上，根据测井数据形态和测井仪器的纵向分辨率，结合岩性分层数据，提取各岩性地层的特征值数据；4)对所有岩性地层的特征值进行数值归约，将特征值相同的相邻岩性层归为一类，生成测井数据岩性层数值归约数据表。

【技术特征摘要】
1.一种测井数据岩性层数值归约方法，包括以下步骤：
1)采用五点二次法对原始测井数据进行平滑处理，得到滤波后的平滑测井数据；
2)采用活度分层法对平滑测井数据进行岩性分层，计算所述平滑测井数据的活度
值，并根据活度截止值得到各个地层的分界点，相邻分界点形成一个岩性地层，进而
得到岩性地层分层数据；
3)根据岩性地层分层数据对原始测井数据进行分段，在每一段原始测井数据上，
根据测井数据形态和测井仪器的纵向分辨率，结合岩性分层数据，提取各岩性地层的
特征值数据；
4)对所有岩性地层的特征值进行数值归约，将特征值相同的相邻岩性层归为一类，
生成测井数据岩性层数值归约数据表。
2.如权利要求1所述的一种测井数据岩性层数值归约方法，其特征在于：所述步
骤2)中，所述平滑测井数据中第i个滑动平均值的活度值Ei为：
Ei=Σk=i-N/2i+N/2[yk-y‾i]2,]]>其中，N为窗长L内平滑测井数据的存取样点数，L为给定活度分层窗长，且
L∈N+，为滑动平均值yi前、后各N/2范围内所有滑动平均值的平均值，即为：
y‾i=1NΣk=i-(N/2)i+(N/2)yk.]]>3...

【专利技术属性】
技术研发人员：张金淼，乔悦东，李洪奇，朱丽萍，朱振宇，郝振江，王建花，糜芳，阴平，
申请(专利权)人：中国海洋石油总公司，中海油研究总院，中国石油大学北京，
类型：发明
国别省市：北京;11