一种获取测井密度曲线的方法技术

本发明专利技术实施例提供了一种获取测井密度曲线的方法，所述获取测井密度曲线方法的具体流程包括：获取目标钻井的自然电位曲线和声波时差曲线；基于小波变换Mallat算法将所述自然电位曲线和所述声波时差曲线进行曲线融合计算，得到所述目标钻井的密度系数曲线；根据所述目标钻井所在区域的密度数据和密度系数曲线计算得到所述目标钻井的密度曲线。本发明专利技术实施例提供的获取测井密度曲线的方法通过有效地提取自然电位曲线和声波时差曲线内分别所包含的岩性及物性特征，采用小波重构的信息融合方式，有效地融合两者，计算得到的测井密度曲线精度更高，更接近实际的密度曲线。

【技术实现步骤摘要】

本专利技术涉及测井密度曲线的处理技术，具体涉及一种获取测井密度曲线的方法。
技术介绍
密度测井就其方法的物理基础而言，是属于放射性测井，在测量时需要使用到伽马射线源，一般是铯137(55Cs137)或钴60(27Co60)，较高的测井作业要求导致在20世纪90年代以前，我国油田开发的密度测井施工作业不多，因此，也导致大部分油井都缺少密度资料，而补充这部分油井的密度资料对石油勘探科技研发是很有意义的工作。如果能够仅仅通过计算获得接近于实际的测井密度曲线，将大大减少科技研发的人力物力，并提高石油勘探的效率。目前，在进行测井解释和地震工作中，可以通过地层速度曲线(声波时差曲线)来求取密度，Gardner曾根据大量实际资料得到了一个由速度推算密度的经验表达式：ρ＝0.31V0.25V为地层速度。在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泥为泥岩基线。可选的，在本专利技术一实施例中，在所述对所述自然电位曲线进行泥岩基线漂移之后，还包括：对所述自然电位曲线进行归一化计算，表达式为：SSP1=SSP-SSPminSSPmax-SSPmin]]>其中，SSP1为归一化的静自然电位曲线，SSPmin为SSP的最小值，SSPmax为SSP的最大值。可选的，在本专利技术一实施例中，在所述基于小波变换Mallat算法将所述自然电位曲线和所述声波时差曲线进行曲线融合计算，还包括：剔除所述声波时差曲线在测量及记录过程中产生的异常值和无效值。可选的，在本专利技术一实施例中，在所述剔除所述声波时差曲线在测量及记录过程中产生的异常值和无效值之后，还包括：对所述声波时差曲线进行低通滤波，去除部分异常高频振荡信号。可选的，在本专利技术一实施例中，在所述对所述声波时差曲线进行低通滤波之后，还包括：对所述声波时差曲线求倒数得到所述声波速度曲线，表达式为：DT=1AC]]>其中，DT为声波速度曲线，AC为声波时差曲线；对所述声波速度曲线进行归一化计算，表达式为：DT1=DT-DTminDTmax-DTmin]]>其中，DT1为归一化的声波速度曲线，DTmin为DT的最小值，DTmax为DT的最大值。可选的，在本专利技术一实施例中，所述根据所述目标钻井所在区域的密度数据以及所述密度系数曲线计算得到所述目标钻井的密度曲线，包括：对所述密度系数曲线归一化，表达式为：DEN_WAV1=DEN_WAV_DEN_WAVminDEN_WAVmax-DEN_WAVmin]]>其中，DEN_WAV为密度系数曲线，DEN_WAV1为归一化的密度系数曲线，DEN_WAVmin为DEN_WAV的最小值，DEN_WAVmax为DEN_WAV的最大值；根据所述目标钻井所在区域的密度数据和所述密度系数曲线，计算得到所述目标钻井的密度曲线，表达式为：DEN＝Amin+(Amax-Amin)×DEN_WAV1其中，DEN为目标钻井的密度曲线，Amin和Amax为目标钻井所在区域地层密度的最小值和最大值。本专利技术实施例提供的获得测井密度曲线的方法，实际是在低频的岩性特征(自然电位)的基础上补充一个高频的孔隙流体特征(声波速度)，有效地融合两者，与传统的Gardner表达式相比，计算得到的结果更为接近实际，精度更高。该方法能够实现多口井的自动批量化处理，节省科技人员宝贵的人力。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1为大庆杏树岗油田X5-33-721井砂岩本文档来自技高网...

【技术保护点】
一种获取测井密度曲线的方法，其特征在于，所述方法包括：获取目标钻井的自然电位曲线和声波时差曲线；基于小波变换Mallat算法将所述自然电位曲线和所述声波时差曲线进行曲线融合计算，得到所述目标钻井的密度系数曲线；根据所述目标钻井所在区域的密度数据以及所述密度系数曲线计算得到所述目标钻井的密度曲线。

【技术特征摘要】
1.一种获取测井密度曲线的方法，其特征在于，所述方法包括：
获取目标钻井的自然电位曲线和声波时差曲线；
基于小波变换Mallat算法将所述自然电位曲线和所述声波时差曲线进行曲线融
合计算，得到所述目标钻井的密度系数曲线；
根据所述目标钻井所在区域的密度数据以及所述密度系数曲线计算得到所述目
标钻井的密度曲线。
2.根据权利要求1所述的获取测井密度曲线的方法，其特征在于，所述曲线融
合计算，包括：
将所述自然电位曲线和声波速度曲线分别当作小波变换Mallat算法中的近似部
分和细节部分，采用Morlet小波系数对所述自然电位曲线和所述声波速度曲线进行
重构计算，所述声波速度曲线为所述声波时差曲线的倒数。
3.根据权利要求1所述的获取测井密度曲线的方法，其特征在于，所述基于小
波变换Mallat算法将所述自然电位曲线和所述声波时差曲线进行曲线融合计算还包
括：剔除所述自然电位曲线在测量及记录过程中产生的异常值和无效值。
4.根据权利要求3所述的获取测井密度曲线的方法，其特征在于，在所述剔除
所述自然电位曲线在测量及记录过程中产生的异常值和无效值之后，还包括：
对所述自然电位曲线进行泥岩基线漂移，表达式为：
SSP＝SP-SP泥其中，SSP为计算得到的静自然电位曲线，SP为原始自然电位曲线，SP泥为泥
岩基线。
5.根据权利要求4所述的获取测井密度曲线的方法，其特征在于，在所述对所
述自然电位曲线进行泥岩基线漂移之后，还包括：
对所述静自然电位曲线进行归一化计算，表达式为：
SSP1=SSP-SSPminSSPmax-SSPmin]]>其中，SSP1为归一化的静自然电位曲线，SSPmin为SSP的最小值，SSPmax为SSP
的最大值。
6.根据权利要求1所述的获取测井密度曲线的方法，其特征在于，所述基于...

【专利技术属性】
技术研发人员：胡水清，刘文岭，宋新民，石成方，田昌炳，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：北京;11