[发明专利]确定测井曲线齿化率的方法及装置

说明书

本发明提供了一种确定测井曲线齿化率的方法及装置。所述方法包括：对砂体的测井曲线进行数据标准化，将数据标准化后的测井曲线形态校正为箱形；根据形态校正后的测井曲线，利用K均值聚类算法确定正齿、负齿及基线的重心值；根据预设的识别门槛值与所述正齿、负齿及基线的重心值，确定正齿齿数及负齿齿数；利用所述正齿齿数及所述负齿齿数，确定所述测井曲线的齿化率。本发明可以准确的确定测井曲线的齿化率，很好的解决通过测井曲线形态齿化现象来判断储层内部结构变化的问题，真实反映储层优劣变化。

技术领域

本发明涉及岩性油藏通过测井曲线识别砂体内部结构变化的技术领域，尤指一种确定测井曲线齿化率的方法及装置。

背景技术

测井曲线的齿化程度表征通常有两种方式，一是纯数学算法表征，例如方差、极差，前者原理是每个样本值与全体样本值的平均数之差的平方值的平均数，也就是齿的幅度，而极差是用来表示统计资料中的变异量数，其最大值与最小值之间的差距，即最大值减最小值后所得之数据，也就是最大齿与最小齿之间的差值，此外就是频率，单位深度范围内存在多少个齿，也是评价齿发育程度的一个重要参数。二是地质学者张广权，李浩《一种利用测井曲线齿化率刻画河道的新方法》中提出的地质上描述齿化率公式，认为齿化是因为岩性变化造成的，并也从幅度、分布密度、分布频率三方面来综合评价，结合权重分析得出齿化率综合表征公式。

但以上方法均存在具体问题，无法准确描述齿化率或没有描述所有类型的齿化，从而使其实用性大打折扣，例如，纯粹的数学方法，没有考虑到齿化的地质成因，也没有进行归类，只是从数学公式中统计，因此，其得出的平均值无法代表地质上的岩性基值，也就无法正确表征齿化幅度，更难以给出“齿”的准确门槛值。而在张广权等提出的地质上描述齿化率公式中，只是指出了一种由岩性变化导致的齿化现象，而对由层理变化，或物性变化导致的齿化现象没有考虑进去，导致其算法具有一定的局限性，而且，其对基线的选择为GR伽马曲线(一种放射性测井曲线系列，常用于区分不同岩性及泥质含量)最小值，也就是岩性最粗部分，而以这部分为基值，则齿的形成未必是差的储层，只是相对最好的部分较差，不能真实反映储层优劣变化，尚有较大的改进空间。

发明内容

为了解决上述问题，本发明实施例提供一种确定测井曲线齿化率的方法，所述方法包括：

对砂体的测井曲线进行数据标准化，将数据标准化后的测井曲线形态校正为箱形；

根据形态校正后的测井曲线，利用K均值聚类算法确定正齿、负齿及基线的重心值；

根据预设的识别门槛值与所述正齿、负齿及基线的重心值，确定正齿齿数及负齿齿数；

利用所述正齿齿数及所述负齿齿数，确定所述测井曲线的齿化率。

可选的，在本发明一实施例中，所述根据预设的识别门槛值与所述正齿、负齿及基线的重心值，确定正齿齿数及负齿齿数包括：分别求取正齿重心值及负齿重心值与基线重心值的差值，将两差值的绝对值分别与预设的识别门槛值进行比较，以确定正齿齿数及负齿齿数。

可选的，在本发明一实施例中，所述利用所述正齿齿数及所述负齿齿数，确定所述测井曲线的齿化率包括：利用所述正齿齿数及所述负齿齿数，确定正齿频率、正齿占比、负齿频率及负齿占比；根据所述正齿频率、正齿占比、负齿频率及负齿占比，确定所述测井曲线的齿化率。

可选的，在本发明一实施例中，所述根据所述正齿频率、正齿占比、负齿频率及负齿占比，确定所述测井曲线的齿化率包括：根据所述正齿频率、正齿占比、负齿频率及负齿占比，利用齿化率公式确定所述测井曲线的齿化率，所述齿化率公式为：

其中，T_正频为正齿频率，T_正比为正齿占比，T_负频为负齿频率，T_负比为负齿占比。

本发明实施例还提供一种确定测井曲线齿化率的装置，所述装置包括：