[发明专利]一种获取测井密度曲线的方法

说明书

技术领域

本发明涉及测井密度曲线的处理技术，具体涉及一种获取测井密度曲线的方法。

背景技术

密度测井就其方法的物理基础而言，是属于放射性测井，在测量时需要使用到伽马射线源，一般是铯¹³⁷(₅₅Cs¹³⁷)或钴⁶⁰(₂₇Co⁶⁰)，较高的测井作业要求导致在20世纪90年代以前，我国油田开发的密度测井施工作业不多，因此，也导致大部分油井都缺少密度资料，而补充这部分油井的密度资料对石油勘探科技研发是很有意义的工作。如果能够仅仅通过计算获得接近于实际的测井密度曲线，将大大减少科技研发的人力物力，并提高石油勘探的效率。

目前，在进行测井解释和地震工作中，可以通过地层速度曲线(声波时差曲线)来求取密度，Gardner曾根据大量实际资料得到了一个由速度推算密度的经验表达式：

ρ＝0.31V^0.25

V为地层速度。在Gardner表达式中，密度与速度表现为指数关系的正相关，当有测井速度的情况下，可以利用Gardner表达式来求取密度。这是地震解释和反演工作中，为制作合成地震记录经常采用的方法。但由于该表达式是根据不同岩性速度和密度的关系，统计得到的一般规律性关系，因此，其计算误差往往较大，在某些地区甚至完全不符合适实际测量结果。比如我国东部高孔高渗的疏松砂泥岩储层，受差异压实作用的影响，往往导致砂泥岩速度重叠，速度与岩性相关程度低，图1为大庆杏树岗油田X5-33-721井砂岩段与泥岩段声波时差的分布范围，可以看出，两种岩性的声波时差分布在(320～380)μs/m之间重合，反映出速度与岩性相关度较差。而实际测量的密度却与岩性具有很高的相关程度，能够较好地区分砂泥岩性，如图2为X5-33-721井实测曲线所示，其中，左数第一道曲线是自然伽马曲线，第二道是自然电位曲线，根据这两条曲线可以解释得到该井的岩性分布，即左数第四道的岩性解释结论。左数第五道和第六道曲线分别为密度曲线和声波时差曲线，从曲线形态可以看出，密度曲线与岩性具有很高的相关程度，砂岩密度较低，泥岩密度高。通过对比发现，表示地层速度的声波时差曲线与岩性无明显相关关系，这种条件下密度与速度之间并不符合指数正相关关系，因此，依据Gardner表达式由速度计算的密度也是不符合真实的地质认识和规律。

发明内容

针对现有技术使用单一速度曲线来计算密度的不足，本发明的目的在于提供一种获取测井密度曲线方法，通过有效提取多种测井曲线内所包含的岩性及物性特征，采用小波重构的信息融合方式，得到更准确更接近实际的密度曲线。

本发明实施例提供了一种获取测井密度曲线的方法，所述方法包括：

获取目标钻井的自然电位曲线和声波时差曲线；

基于小波变换Mallat算法将所述自然电位曲线和所述声波时差曲线进行曲线融合计算，得到所述目标钻井的密度系数曲线；

根据所述目标钻井所在区域的密度数据和所述密度系数曲线计算得到所述目标钻井的密度曲线。

可选的，在本发明一实施例中，所述曲线融合计算，包括：

将所述自然电位曲线和声波速度曲线分别当作小波变换Mallat算法中的近似部分和细节部分，采用Morlet小波系数进行重构计算，所述声波速度曲线为所述声波时差曲线的倒数。

可选的，在本发明一实施例中，所述基于小波变换Mallat算法将所述自然电位曲线和所述声波时差曲线进行曲线融合计算还包括：剔除所述自然电位曲线在测量及记录过程中产生的异常值和无效值。

可选的，在本发明一实施例中，在所述剔除所述自然电位曲线在测量及记录过程中产生的异常值和无效值之后，还包括：

对所述自然电位曲线进行泥岩基线漂移，表达式为：

SSP＝SP-SP_泥

其中，SSP为计算得到的静自然电位曲线，SP为原始自然电位曲线，SP_泥为泥岩基线。

可选的，在本发明一实施例中，在所述对所述自然电位曲线进行泥岩基线漂移之后，还包括：

对所述自然电位曲线进行归一化计算，表达式为：