CN116009110A - 密度曲线校正方法、装置、存储介质及电子设备 - Google Patents

Abstract

本申请涉及非常规油气勘探开发技术领域，具体涉及一种密度曲线校正方法、装置及电子设备，该方法首先根据待校正储层的测井曲线信息划分出密度曲线失真段，然后在密度曲线未失真段，分析密度测井曲线与总有机碳数据之间的经验关系，同时获取基于自然伽马能谱的分析化验总有机碳计算模型。然后在整个页岩层段，根据经验关系式以及分析化验总有机碳计算模型进行计算，便可得到反映真实页岩地层的密度曲线，即密度校正曲线。上述方式以总有机碳数据为桥梁，利用自然伽马能谱测井曲线、密度曲线以及总有机碳数据，实现对待校正储层的密度曲线进行准确校正的目的。

Description

密度曲线校正方法、装置、存储介质及电子设备

技术领域

本申请涉及非常规油气勘探开发技术领域，特别地涉及一种密度曲线校正方法、装置、存储介质及电子设备。

背景技术

在页岩气储层钻井过程中，泥页岩易吸水膨胀，井壁强度降低，井眼易扩径坍塌。同时密度测井仪器贴井壁测量，在井眼扩径垮塌处，密度曲线探测深度较浅，极易发生失真，不能反映地层的真实密度特征。而密度曲线又是进行孔隙度测井解释以及合成地震记录不可缺少的曲线。

目前主要采用以下几种方法对密度曲线进行校正。第一种，采用多元回归法。该方法在标准层内利用曲线之间的内在关系，用多元回归分析的方法把畸变失真的曲线校正到真实的测量范围。但该方法存在一定的局限性，标准层内的参考曲线与目标曲线相关性不一定满足精度要求，并且也缺乏物理意义。第二种，岩石物理模型法。该方法根据测井评价得到的各矿物含量、孔隙度和饱和度等参数，然后建立密度响应的岩石物理体积模型，该方法应用资料过多，容易在测井参数的评价过程中引入新的误差，存在不准确的问题。第三章，Gardner经验公式法。该方法利用测井数据进行统计，得到密度测井和纵波曲线之间的对应关系。然而该方法不具有普适性，很难进行推广应用。

发明内容

针对上述问题，本申请提供一种密度曲线校正方法、装置、存储介质及电子设备。

第一方面，本申请提供了一种密度曲线校正方法，所述方法包括：

获取待校正储层的测井曲线信息，并根据所述测井曲线信息确定出所述测井曲线信息中的密度曲线失真段及密度曲线未失真段；

基于所述密度曲线未失真段的密度测井曲线以及总有机碳数据建立经验关系式；

获取所述待校正储层基于自然伽马能谱的分析化验总有机碳计算模型；

基于所述经验关系式以及所述分析化验总有机碳计算模型确定与所述密度曲线失真段对应的密度校正曲线。

上述实施方式中，首先根据待校正储层的测井曲线信息划分出密度曲线失真段，然后在密度曲线未失真段，分析密度测井曲线与总有机碳数据之间的经验关系，同时获取基于自然伽马能谱的分析化验总有机碳计算模型。然后在整个页岩层段，根据经验关系式以及分析化验总有机碳计算模型进行计算，便可得到反映真实页岩地层的密度曲线，即密度校正曲线。上述方式以总有机碳数据为桥梁，利用自然伽马能谱测井曲线、密度曲线以及总有机碳数据，实现对待校正储层的密度曲线进行准确校正的目的。

根据本申请的实施例，可选的，上述密度曲线校正方法中，所述根据所述测井曲线信息确定出所述测井曲线信息中的密度曲线失真段及密度曲线未失真段的步骤，包括：

确定所述测井曲线信息中的密度曲线、去铀伽马曲线以及井径曲线；

将所述密度曲线与所述去铀伽马曲线进行比对，以获取比对结果；

根据所述比对结果以及所述井径曲线的扩径程度确定出所述测井曲线信息的密度曲线失真段及密度曲线未失真段。

上述实施方式中，测井曲线信息中的井径曲线能直接反映待校正储层中井眼扩径的程度，其中的去铀伽马曲线则反映的是地层的岩性特征。地层的密度曲线与该地层的去铀伽马曲线存在着统计性的关系，同时，去铀伽马曲线的径向探测距离相对较远，受井眼扩径的影响较小，不易失真。因此可以对比密度曲线与去铀伽马曲线之间的对应关系，然后根据比对结果再结合井径曲线的扩径程度，确定出密度曲线的失真段。

根据本申请的实施例，可选的，上述密度曲线校正方法中，所述基于所述密度曲线未失真段的密度测井曲线以及总有机碳数据建立经验关系式的步骤，包括：

基于所述密度曲线未失真段的密度测井曲线获取反映所述待校正储层真实密度信息的密度测井曲线值；

获取总有机碳数据中与所述密度曲线未失真段对应的分析化验总有机碳；

根据所述密度测井曲线值与所述分析化验总有机碳确定经验关系式。

上述实施方式中，对于页岩地层，其总有机碳的含量越高，则地层密度越小，总有机碳与地层密度两者之间具有较强的负相关性，因此可以基于所述密度曲线未失真段的密度测井曲线以及总有机碳数据建立经验关系式，经验关系可以准确的表示出待校正储层的密度测井曲线以及总有机碳数据之间的关系。

根据本申请的实施例，可选的，上述密度曲线校正方法中，所述根据所述密度测井曲线值与所述分析化验总有机碳确定经验关系式的步骤，包括：

根据所述密度测井曲线值、所述分析化验总有机碳以及第一计算式确定所述经验关系式；

所述经验关系式为：

DEN＝a*TOC+b

其中，DEN表示密度测井曲线值，TOC表示分析化验总有机碳，a表示第一经验参数，b表示第二经验参数。

根据本申请的实施例，可选的，上述密度曲线校正方法中，所述获取所述待校正储层基于自然伽马能谱的分析化验总有机碳计算模型的步骤之前，所述方法还包括：

获取所述待校正储层的自然伽马能谱测井曲线数据和分析化验总有机碳；

根据所述自然伽马能谱测井曲线数据和所述分析化验总有机碳建立分析化验总有机碳计算模型。

上述实施方式中，由于自然伽马能谱测井曲线径向探测距离相对较远，受井眼扩径影响较小，不易失真。因此，可以在待校正储层的整个页岩层段，根据自然伽马能谱测井曲线和分析化验总有机碳，建立基于自然伽马能谱的分析化验总有机碳计算模型。

根据本申请的实施例，可选的，上述密度曲线校正方法中，所述根据所述自然伽马能谱测井曲线数据和所述分析化验总有机碳建立分析化验总有机碳计算模型的步骤，包括：

确定所述自然伽马能谱测井曲线数据中的自然伽马能谱测井总伽马曲线、去铀伽马曲线、钍曲线和铀曲线；

根据所述自然伽马能谱测井总伽马曲线、所述去铀伽马曲线、钍曲线、铀曲线、所述分析化验总有机碳以及第二计算式获取所述分析化验总有机碳计算模型；

所述第二计算式为：

其中，TOC表示分析化验总有机碳，GR表示自然伽马能谱测井总伽马曲线，GR_min表示自然伽马能谱测井总伽马曲线以及去铀伽马曲线重叠时自然伽马能谱测井总伽马曲线的最小刻度，GR_max表示自然伽马能谱测井总伽马曲线以及去铀伽马曲线重叠时自然伽马能谱测井总伽马曲线的最大刻度，KTH表示去铀伽马曲线，KTH_min表示自然伽马能谱测井总伽马曲线以及去铀伽马曲线重叠时去铀伽马曲线的最小刻度，KTH_max表示自然伽马能谱测井总伽马曲线以及去铀伽马曲线重叠时去铀伽马曲线的最大刻度，A表示第三经验参数，B表示第四经验参数，TH表示钍曲线，U表示铀曲线。

根据本申请的实施例，可选的，上述密度曲线校正方法中，所述基于所述经验关系式以及所述分析化验总有机碳计算模型确定与所述密度曲线失真段对应的密度校正曲线的步骤，包括：

根据所述分析化验总有机碳计算模型计算所述密度曲线失真段对应的总有机碳含量曲线；

根据所述总有机碳含量曲线以及所述经验关系式确定与所述密度曲线失真段对应的密度校正曲线。

上述实施方式中，在密度曲线失真的井段，根据分析化验总有机碳计算模型进行计算可以得到能反映原状地层总有机碳含量的总有机碳含量曲线，然后将该总有机碳含量曲线代入基于所述经验关系式，即可以准确的得到与密度曲线失真段对应的密度校正曲线。

第二方面，本申请还提供了一种密度曲线校正装置，所述装置包括：失真段确定模块，用于获取待校正储层的测井曲线信息，并根据所述测井曲线信息确定出所述测井曲线信息中的密度曲线失真段及密度曲线未失真段；

经验关系建立模块，用于基于所述密度曲线未失真段的密度测井曲线以及总有机碳数据建立经验关系式；

分析化验总有机碳计算模型获取模块，用于获取所述待校正储层基于自然伽马能谱的分析化验总有机碳计算模型；

密度校正曲线确定模块，用于基于所述经验关系式以及所述分析化验总有机碳计算模型确定与所述密度曲线失真段对应的密度校正曲线。

根据本申请的实施例，可选的，上述密度曲线校正装置中，所述失真段确定模块包括：

曲线确定单元，用于确定所述测井曲线信息中的密度曲线、去铀伽马曲线以及井径曲线；

比对单元，用于将所述密度曲线与所述去铀伽马曲线进行比对，以获取比对结果；

失真段确定单元，用于根据所述比对结果以及所述井径曲线的扩径程度确定出所述测井曲线信息的密度曲线失真段及密度曲线未失真段。

根据本申请的实施例，可选的，上述密度曲线校正装置中，所述经验关系建立模块包括：

密度测井曲线值获取单元，用于基于所述密度曲线未失真段的密度测井曲线获取反映所述待校正储层真实密度信息的密度测井曲线值；

分析化验总有机碳获取单元，用于获取总有机碳数据中与所述密度曲线未失真段对应的分析化验总有机碳；

经验关系式确定单元，用于根据所述密度测井曲线值与所述分析化验总有机碳确定经验关系式。

根据本申请的实施例，可选的，上述密度曲线校正装置中，所述经验关系式确定单元包括：

经验关系式确定子单元，用于根据所述密度测井曲线值、所述分析化验总有机碳以及第一计算式确定所述经验关系式；

所述经验关系式为：

DEN＝a*TOC+b

其中，DEN表示密度测井曲线值，TOC表示分析化验总有机碳，a表示第一经验参数，b表示第二经验参数。

根据本申请的实施例，可选的，上述密度曲线校正装置中，所述装置还包括：

数据获取模块，用于获取所述待校正储层的自然伽马能谱测井曲线数据和分析化验总有机碳；

分析化验总有机碳计算模型建立模块，用于根据所述自然伽马能谱测井曲线数据和所述分析化验总有机碳建立分析化验总有机碳计算模型。

根据本申请的实施例，可选的，上述密度曲线校正装置中，所述分析化验总有机碳计算模型建立模块包括：

自然伽马能谱测井曲线确定单元，用于确定所述自然伽马能谱测井曲线数据中的自然伽马能谱测井总伽马曲线、去铀伽马曲线、钍曲线和铀曲线；

分析化验总有机碳计算模型获取单元，用于根据所述自然伽马能谱测井总伽马曲线、所述去铀伽马曲线、钍曲线、铀曲线、所述分析化验总有机碳以及第二计算式获取所述分析化验总有机碳计算模型；

所述第二计算式为：

其中，TOC表示分析化验总有机碳，GR表示自然伽马能谱测井总伽马曲线，GR_min表示自然伽马能谱测井总伽马曲线以及去铀伽马曲线重叠时自然伽马能谱测井总伽马曲线的最小刻度，GR_max表示自然伽马能谱测井总伽马曲线以及去铀伽马曲线重叠时自然伽马能谱测井总伽马曲线的最大刻度，KTH表示去铀伽马曲线，KTH_min表示自然伽马能谱测井总伽马曲线以及去铀伽马曲线重叠时去铀伽马曲线的最小刻度，KTH_max表示自然伽马能谱测井总伽马曲线以及去铀伽马曲线重叠时去铀伽马曲线的最大刻度，A表示第三经验参数，B表示第四经验参数，TH表示钍曲线，U表示铀曲线。

根据本申请的实施例，可选的，上述密度曲线校正装置中，所述密度校正曲线确定模块包括：

总有机碳含量曲线计算单元，用于根据所述分析化验总有机碳计算模型计算所述密度曲线失真段对应的总有机碳含量曲线；

密度校正曲线确定单元，用于根据所述总有机碳含量曲线以及所述经验关系式确定与所述密度曲线失真段对应的密度校正曲线。

第三方面，本申请提供了一种存储介质，该存储介质存储的计算机程序，可被一个或多个处理器执行，可用来实现如上述的密度曲线校正方法。

第四方面，本申请提供了一种电子设备，包括存储器和处理器，所述存储器上存储有计算机程序，该计算机程序被所述处理器执行时，执行上述的密度曲线校正方法。

与现有技术相比，上述方案中的一个或多个实施例可以具有如下优点或有益效果：

本申请提供的一种密度曲线校正方法、装置、存储介质及电子设备，该方法包括：获取待校正储层的测井曲线信息，并根据所述测井曲线信息确定出所述测井曲线信息中的密度曲线失真段及密度曲线未失真段；基于所述密度曲线未失真段的密度测井曲线以及总有机碳数据建立经验关系式；获取所述待校正储层基于自然伽马能谱的分析化验总有机碳计算模型；基于所述经验关系式以及所述分析化验总有机碳计算模型确定与所述密度曲线失真段对应的密度校正曲线。首先根据待校正储层的测井曲线信息划分出密度曲线失真段，然后在密度曲线未失真段，分析密度测井曲线与总有机碳数据之间的经验关系，同时获取基于自然伽马能谱的分析化验总有机碳计算模型。然后在整个页岩层段，根据经验关系式以及分析化验总有机碳计算模型进行计算，便可得到反映真实页岩地层的密度曲线，即密度校正曲线。上述方式以总有机碳数据为桥梁，利用自然伽马能谱测井曲线、密度曲线以及总有机碳数据，实现对待校正储层的密度曲线进行准确校正的目的。

附图说明

在下文中将基于实施例并参考附图来对本申请进行更详细的描述。

图1为本申请实施例一提供的一种密度曲线校正方法的流程示意图。

图2为本申请实施例三提供的一种去铀伽马曲线和密度曲线的对比分析图。

图3为本申请实施例三提供的一种密度曲线失真段示意图。

图4为本申请实施例三提供的一种分析密度测井曲线与总有机碳数据分析图。

图5为本申请实施例三提供的一种密度校正曲线图。

图6为本申请实施例四提供的一种密度曲线校正装置的连接框图。

图7为本申请实施例六提供的一种电子设备的连接框图。

在附图中，相同的部件使用相同的附图标记，附图并未按照实际的比例绘制。

具体实施方式

以下将结合附图及实施例来详细说明本申请的实施方式，借此对本申请如何应用技术手段来解决技术问题，并达到相应技术效果的实现过程能充分理解并据以实施。本申请实施例以及实施例中的各个特征，在不相冲突前提下可以相互结合，所形成的技术方案均在本申请的保护范围之内。

实施例一

本发明提供一种密度曲线校正方法，请参阅图1，该方法包括如下步骤：

步骤S110：获取待校正储层的测井曲线信息，并根据所述测井曲线信息确定出所述测井曲线信息中的密度曲线失真段及密度曲线未失真段。

地层密度的大小与地层岩性、孔隙和流体等因素息息相关，其中受地层岩性的影响最大。测井曲线信息中的数据可以反应待校正储层的井眼扩径的程度、地层的岩性特征等特征，因此，可以根据测井曲线信息确定出测井曲线信息中的密度曲线失真段。

根据本申请的实施例，可选的，上述密度曲线校正方法中，所述根据所述测井曲线信息确定出所述测井曲线信息中的密度曲线失真段及密度曲线未失真段的步骤S110，包括如下步骤：

步骤S111：确定所述测井曲线信息中的密度曲线、去铀伽马曲线以及井径曲线；

步骤S112：将所述密度曲线与所述去铀伽马曲线进行比对，以获取比对结果；

步骤S113：根据所述比对结果以及所述井径曲线的扩径程度确定出所述测井曲线信息的密度曲线失真段及密度曲线未失真段。

测井曲线信息中的井径曲线能直接反映待校正储层中井眼扩径的程度，其中的去铀伽马曲线则反映的是地层的岩性特征。地层的密度曲线与该地层的去铀伽马曲线存在着统计性的关系，同时，去铀伽马曲线的径向探测距离相对较远，受井眼扩径的影响较小，不易失真。因此可以对比密度曲线与去铀伽马曲线之间的对应关系，然后根据比对结果再结合井径曲线的扩径程度，即能够准确的确定出密度曲线的失真段。在确定出密度曲线的失真段之后，即可得到密度曲线的未失真段。

步骤S120：基于所述密度曲线未失真段的密度测井曲线以及总有机碳数据建立经验关系式。

对于页岩地层，其总有机碳的含量越高，则地层密度越小，总有机碳与地层密度两者之间具有较强的负相关性，因此可以基于所述密度曲线未失真段的密度测井曲线以及总有机碳数据建立经验关系式，经验关系可以准确的表示出待校正储层的密度测井曲线以及总有机碳数据之间的关系。

步骤S130：获取所述待校正储层基于自然伽马能谱的分析化验总有机碳计算模型。

由于自然伽马能谱测井曲线径向探测距离相对较远，受井眼扩径影响较小，不易失真。因此，可以在待校正储层的整个页岩层段，根据自然伽马能谱测井曲线和分析化验总有机碳，建立基于自然伽马能谱的分析化验总有机碳计算模型。

步骤S140：基于所述经验关系式以及所述分析化验总有机碳计算模型确定与所述密度曲线失真段对应的密度校正曲线。

可选的，上述密度曲线校正方法中，所述基于所述经验关系式以及所述分析化验总有机碳计算模型确定与所述密度曲线失真段对应的密度校正曲线的步骤S140，包括以下步骤：

步骤S141：根据所述分析化验总有机碳计算模型计算所述密度曲线失真段对应的总有机碳含量曲线。

步骤S142：根据所述总有机碳含量曲线以及所述经验关系式确定与所述密度曲线失真段对应的密度校正曲线。

上述实施方式中，在密度曲线失真的井段，根据分析化验总有机碳计算模型进行计算可以得到能反映原状地层总有机碳含量的总有机碳含量曲线，然后将该总有机碳含量曲线代入基于所述经验关系式，即可以准确的得到与密度曲线失真段对应的密度校正曲线。

综上所述，本申请提供一种密度曲线校正方法，包括：获取待校正储层的测井曲线信息，并根据所述测井曲线信息确定出所述测井曲线信息中的密度曲线失真段及密度曲线未失真段；基于所述密度曲线未失真段的密度测井曲线以及总有机碳数据建立经验关系式；获取所述待校正储层基于自然伽马能谱的分析化验总有机碳计算模型；基于所述经验关系式以及所述分析化验总有机碳计算模型确定与所述密度曲线失真段对应的密度校正曲线。上述实施方式中，首先根据待校正储层的测井曲线信息划分出密度曲线失真段，然后在密度曲线未失真段，分析密度测井曲线与总有机碳数据之间的经验关系，同时获取基于自然伽马能谱的分析化验总有机碳计算模型。然后在整个页岩层段，根据经验关系式以及分析化验总有机碳计算模型进行计算，便可得到反映真实页岩地层的密度曲线，即密度校正曲线。上述方式以总有机碳数据为桥梁，利用自然伽马能谱测井曲线、密度曲线以及总有机碳数据，实现对待校正储层的密度曲线进行准确校正的目的。

实施例二

在实施例一的基础上，本实施例通过具体实施案例对实施例一中的方法进行说明。

根据本申请的实施例，可选的，上述密度曲线校正方法中，所述基于所述密度曲线未失真段的密度测井曲线以及总有机碳数据建立经验关系式的步骤S120，包括如下步骤：

步骤S121：基于所述密度曲线未失真段的密度测井曲线获取反映所述待校正储层真实密度信息的密度测井曲线值。

步骤S122：获取总有机碳数据中与所述密度曲线未失真段对应的分析化验总有机碳。

步骤S123：根据所述密度测井曲线值与所述分析化验总有机碳确定经验关系式。

作为一种实施方式，可选的，上述密度曲线校正方法中，所述根据所述密度测井曲线值与所述分析化验总有机碳确定经验关系式的步骤，包括：

根据所述密度测井曲线值、所述分析化验总有机碳以及第一计算式确定所述经验关系式；

所述经验关系式为：

DEN＝a*TOC+b

其中，DEN表示密度测井曲线值，TOC表示分析化验总有机碳，a表示第一经验参数，b表示第二经验参数。

根据本申请的实施例，可选的，上述密度曲线校正方法中，所述获取所述待校正储层基于自然伽马能谱的分析化验总有机碳计算模型的步骤之前，所述方法还包括以下过程。首先获取所述待校正储层的自然伽马能谱测井曲线数据和分析化验总有机碳，然后根据所述自然伽马能谱测井曲线数据和所述分析化验总有机碳建立分析化验总有机碳计算模型。

作为一种实施方式，可选的，上述密度曲线校正方法中，所述根据所述自然伽马能谱测井曲线数据和所述分析化验总有机碳建立分析化验总有机碳计算模型的步骤，包括以下过程。首先确定所述自然伽马能谱测井曲线数据中的自然伽马能谱测井总伽马曲线、去铀伽马曲线、钍曲线和铀曲线。然后根据所述自然伽马能谱测井总伽马曲线、所述去铀伽马曲线、钍曲线、铀曲线、所述分析化验总有机碳以及第二计算式获取所述分析化验总有机碳计算模型。其中所述第二计算式为：

在第二计算式中，TOC表示分析化验总有机碳，GR表示自然伽马能谱测井总伽马曲线，GR_min表示自然伽马能谱测井总伽马曲线以及去铀伽马曲线重叠时自然伽马能谱测井总伽马曲线的最小刻度，GR_max表示自然伽马能谱测井总伽马曲线以及去铀伽马曲线重叠时自然伽马能谱测井总伽马曲线的最大刻度，KTH表示去铀伽马曲线，KTH_min表示自然伽马能谱测井总伽马曲线以及去铀伽马曲线重叠时去铀伽马曲线的最小刻度，KTH_max表示自然伽马能谱测井总伽马曲线以及去铀伽马曲线重叠时去铀伽马曲线的最大刻度，A表示第三经验参数，B表示第四经验参数，TH表示钍曲线，U表示铀曲线。

实施例三

在实施例一的基础上，本实施例通过具体实施案例对实施例一中的方法进行说明。

以A1井为例，A1井位四川盆地武陵褶皱带桑柘坪向斜，目的层段为上奥陶统五峰组-下志留统龙马溪组页岩，从下往上，依次可划分为9个小层。首先进行去铀伽马曲线和密度曲线的对比分析，如图2所示，图2中横坐标为地层的去铀伽马值，纵坐标为地层的密度值，从图2中可以看出，绝大多数样本点对应的密度值和去铀伽马值具有正相关性，但还有部分样本点具有较低的密度值，如图中黑色线框内的样本点所示。这些具有异常密度值的样本点便可被认为是密度曲线失真的井段。因此，确定所述测井曲线信息中的密度曲线、去铀伽马曲线以及井径曲线之后，将所述密度曲线与所述去铀伽马曲线进行比对，以获取比对结果。根据所述比对结果以及所述井径曲线的扩径程度确定出所述测井曲线信息的密度曲线失真段及密度曲线未失真段。例如，综合利用去铀伽马KTH-密度DEN交会图以及井径曲线可确定出A1井的密度曲线失真段，如图3所示。

图3中第一道为深度道(MD)；第二道为测井解释小层道(小层)；第三道中实线为总伽马曲线(GR)，点虚线为去铀伽马曲线(KTH)；第四道中点虚线为井径曲线(CAL)，垂直实线代表钻头直径，为8.5in，灰色填充面积反映井眼扩径的程度；第五道中实线为密度曲线(DEN)；第六道中杆状线为分析化验总有机碳TOC含量(TOC)；第七道为上述步骤确定出的密度曲线失真段。

在密度未失真的井段，分析密度测井曲线与总有机碳数据之间的关系(如图4所示)，从图4中可以看出，密度测井曲线与总有机碳数据两者之间具有较高的相关性，说明利用总有机碳数据能够较准确的得到原状地层的密度信息，相应的计算公式为DEN＝-0.0385*TOC+2.6653，即基于所述密度曲线未失真段的密度测井曲线以及总有机碳数据建立经验关系式DEN＝a*TOC+b中的a＝-0.0385，b＝2.6653。

再在整个页岩层段，依据分析化验总有机碳计算模型利用自然伽马能谱测井曲线进行总有机碳含量的计算，分析化验总有机碳计算模型中相应的经验参数A＝0.7483，B＝0.1124。在密度曲线失真段，将依据分析化验总有机碳计算模型计算得到的总有机碳含量代入经验关系中，即可得到失真段的密度曲线校正值。相应的密度校正曲线如图5所示。图中第五道实线为利用自然伽马能谱计算的总有机碳含量(TOCM)，与实验分析得到的总有机碳(杆状线)具有较好的吻合性，说明利用自然伽马能谱能准确的计算页岩储层的TOC含量，图中第六道为校正后的密度曲线，其中密度曲线失真段采用的是密度曲线校正值，密度曲线未失真段采用的是密度曲线测井值。

实施例四

请参看图6，本申请提供了一种密度曲线校正装置600，该装置包括：

失真段确定模块610，用于获取待校正储层的测井曲线信息，并根据所述测井曲线信息确定出所述测井曲线信息中的密度曲线失真段及密度曲线未失真段；

经验关系建立模块620，用于基于所述密度曲线未失真段的密度测井曲线以及总有机碳数据建立经验关系式；

分析化验总有机碳计算模型获取模块630，用于获取所述待校正储层基于自然伽马能谱的分析化验总有机碳计算模型；

密度校正曲线确定模块640，用于基于所述经验关系式以及所述分析化验总有机碳计算模型确定与所述密度曲线失真段对应的密度校正曲线。

根据本申请的实施例，可选的，上述密度曲线校正装置中，所述失真段确定模块包括：

曲线确定单元，用于确定所述测井曲线信息中的密度曲线、去铀伽马曲线以及井径曲线；

比对单元，用于将所述密度曲线与所述去铀伽马曲线进行比对，以获取比对结果；

失真段确定单元，用于根据所述比对结果以及所述井径曲线的扩径程度确定出所述测井曲线信息的密度曲线失真段及密度曲线未失真段。

根据本申请的实施例，可选的，上述密度曲线校正装置中，所述经验关系建立模块包括：

密度测井曲线值获取单元，用于基于所述密度曲线未失真段的密度测井曲线获取反映所述待校正储层真实密度信息的密度测井曲线值；

分析化验总有机碳获取单元，用于获取总有机碳数据中与所述密度曲线未失真段对应的分析化验总有机碳；

经验关系式确定单元，用于根据所述密度测井曲线值与所述分析化验总有机碳确定经验关系式。

根据本申请的实施例，可选的，上述密度曲线校正装置中，所述经验关系式确定单元包括：

经验关系式确定子单元，用于根据所述密度测井曲线值、所述分析化验总有机碳以及第一计算式确定所述经验关系式；

所述经验关系式为：

DEN＝a*TOC+b

其中，DEN表示密度测井曲线值，TOC表示分析化验总有机碳，a表示第一经验参数，b表示第二经验参数。

根据本申请的实施例，可选的，上述密度曲线校正装置中，所述装置还包括：

数据获取模块，用于获取所述待校正储层的自然伽马能谱测井曲线数据和分析化验总有机碳；

分析化验总有机碳计算模型建立模块，用于根据所述自然伽马能谱测井曲线数据和所述分析化验总有机碳建立分析化验总有机碳计算模型。

根据本申请的实施例，可选的，上述密度曲线校正装置中，所述分析化验总有机碳计算模型建立模块包括：

自然伽马能谱测井曲线确定单元，用于确定所述自然伽马能谱测井曲线数据中的自然伽马能谱测井总伽马曲线、去铀伽马曲线、钍曲线和铀曲线；

分析化验总有机碳计算模型获取单元，用于根据所述自然伽马能谱测井总伽马曲线、所述去铀伽马曲线、钍曲线、铀曲线、所述分析化验总有机碳以及第二计算式获取所述分析化验总有机碳计算模型；

所述第二计算式为：

其中，TOC表示分析化验总有机碳，GR表示自然伽马能谱测井总伽马曲线，GR_min表示自然伽马能谱测井总伽马曲线以及去铀伽马曲线重叠时自然伽马能谱测井总伽马曲线的最小刻度，GR_max表示自然伽马能谱测井总伽马曲线以及去铀伽马曲线重叠时自然伽马能谱测井总伽马曲线的最大刻度，KTH表示去铀伽马曲线，KTH_min表示自然伽马能谱测井总伽马曲线以及去铀伽马曲线重叠时去铀伽马曲线的最小刻度，KTH_max表示自然伽马能谱测井总伽马曲线以及去铀伽马曲线重叠时去铀伽马曲线的最大刻度，A表示第三经验参数，B表示第四经验参数，TH表示钍曲线，U表示铀曲线。

根据本申请的实施例，可选的，上述密度曲线校正装置中，所述密度校正曲线确定模块包括：

总有机碳含量曲线计算单元，用于根据所述分析化验总有机碳计算模型计算所述密度曲线失真段对应的总有机碳含量曲线；

密度校正曲线确定单元，用于根据所述总有机碳含量曲线以及所述经验关系式确定与所述密度曲线失真段对应的密度校正曲线。

综上所述，本申请提供一种密度曲线校正装置，包括：失真段确定模块610，用于获取待校正储层的测井曲线信息，并根据所述测井曲线信息确定出所述测井曲线信息中的密度曲线失真段及密度曲线未失真段；经验关系建立模块620，用于基于所述密度曲线未失真段的密度测井曲线以及总有机碳数据建立经验关系式；分析化验总有机碳计算模型获取模块630，用于获取所述待校正储层基于自然伽马能谱的分析化验总有机碳计算模型；密度校正曲线确定模块640，用于基于所述经验关系式以及所述分析化验总有机碳计算模型确定与所述密度曲线失真段对应的密度校正曲线。首先根据待校正储层的测井曲线信息划分出密度曲线失真段，然后在密度曲线未失真段，分析密度测井曲线与总有机碳数据之间的经验关系，同时获取基于自然伽马能谱的分析化验总有机碳计算模型。然后在整个页岩层段，根据经验关系式以及分析化验总有机碳计算模型进行计算，便可得到反映真实页岩地层的密度曲线，即密度校正曲线。上述方式以总有机碳数据为桥梁，利用自然伽马能谱测井曲线、密度曲线以及总有机碳数据，实现对待校正储层的密度曲线进行准确校正的目的。

实施例五

本实施例还提供一种计算机可读存储介质，如闪存、硬盘、多媒体卡、卡型存储器(例如，SD或DX存储器等)、随机访问存储器(RAM)、静态随机访问存储器(SRAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、可编程只读存储器(PROM)、磁性存储器、磁盘、光盘、服务器、App应用商城等等，其上存储有计算机程序，所述计算机程序被处理器执行时可以实现如下方法步骤：

步骤S110：获取待校正储层的测井曲线信息，并根据所述测井曲线信息确定出所述测井曲线信息中的密度曲线失真段及密度曲线未失真段；

步骤S120：基于所述密度曲线未失真段的密度测井曲线以及总有机碳数据建立经验关系式；

步骤S130：获取所述待校正储层基于自然伽马能谱的分析化验总有机碳计算模型；

步骤S140：基于所述经验关系式以及所述分析化验总有机碳计算模型确定与所述密度曲线失真段对应的密度校正曲线。

上述实施方式中，首先根据待校正储层的测井曲线信息划分出密度曲线失真段，然后在密度曲线未失真段，分析密度测井曲线与总有机碳数据之间的经验关系，同时获取基于自然伽马能谱的分析化验总有机碳计算模型。然后在整个页岩层段，根据经验关系式以及分析化验总有机碳计算模型进行计算，便可得到反映真实页岩地层的密度曲线，即密度校正曲线。上述方式以总有机碳数据为桥梁，利用自然伽马能谱测井曲线、密度曲线以及总有机碳数据，实现对待校正储层的密度曲线进行准确校正的目的。

可选的，上述密度曲线校正方法中，所述根据所述测井曲线信息确定出所述测井曲线信息中的密度曲线失真段及密度曲线未失真段的步骤，包括：

确定所述测井曲线信息中的密度曲线、去铀伽马曲线以及井径曲线；

将所述密度曲线与所述去铀伽马曲线进行比对，以获取比对结果；

根据所述比对结果以及所述井径曲线的扩径程度确定出所述测井曲线信息的密度曲线失真段及密度曲线未失真段。

上述实施方式中，测井曲线信息中的井径曲线能直接反映待校正储层中井眼扩径的程度，其中的去铀伽马曲线则反映的是地层的岩性特征。地层的密度曲线与该地层的去铀伽马曲线存在着统计性的关系，同时，去铀伽马曲线的径向探测距离相对较远，受井眼扩径的影响较小，不易失真。因此可以对比密度曲线与去铀伽马曲线之间的对应关系，然后根据比对结果再结合井径曲线的扩径程度，确定出密度曲线的失真段。

可选的，上述密度曲线校正方法中，所述基于所述密度曲线未失真段的密度测井曲线以及总有机碳数据建立经验关系式的步骤，包括：

基于所述密度曲线未失真段的密度测井曲线获取反映所述待校正储层真实密度信息的密度测井曲线值；

获取总有机碳数据中与所述密度曲线未失真段对应的分析化验总有机碳；

根据所述密度测井曲线值与所述分析化验总有机碳确定经验关系式。

上述实施方式中，对于页岩地层，其总有机碳的含量越高，则地层密度越小，总有机碳与地层密度两者之间具有较强的负相关性，因此可以基于所述密度曲线未失真段的密度测井曲线以及总有机碳数据建立经验关系式，经验关系可以准确的表示出待校正储层的密度测井曲线以及总有机碳数据之间的关系。

可选的，上述密度曲线校正方法中，所述根据所述密度测井曲线值与所述分析化验总有机碳确定经验关系式的步骤，包括：

根据所述密度测井曲线值、所述分析化验总有机碳以及第一计算式确定所述经验关系式；

所述经验关系式为：

DEN＝a*TOC+b

其中，DEN表示密度测井曲线值，TOC表示分析化验总有机碳，a表示第一经验参数，b表示第二经验参数。

可选的，上述密度曲线校正方法中，所述获取所述待校正储层基于自然伽马能谱的分析化验总有机碳计算模型的步骤之前，所述方法还包括：

获取所述待校正储层的自然伽马能谱测井曲线数据和分析化验总有机碳；

根据所述自然伽马能谱测井曲线数据和所述分析化验总有机碳建立分析化验总有机碳计算模型。

上述实施方式中，由于自然伽马能谱测井曲线径向探测距离相对较远，受井眼扩径影响较小，不易失真。因此，可以在待校正储层的整个页岩层段，根据自然伽马能谱测井曲线和分析化验总有机碳，建立基于自然伽马能谱的分析化验总有机碳计算模型。

可选的，上述密度曲线校正方法中，所述根据所述自然伽马能谱测井曲线数据和所述分析化验总有机碳建立分析化验总有机碳计算模型的步骤，包括：

确定所述自然伽马能谱测井曲线数据中的自然伽马能谱测井总伽马曲线、去铀伽马曲线、钍曲线和铀曲线；

根据所述自然伽马能谱测井总伽马曲线、所述去铀伽马曲线、钍曲线、铀曲线、所述分析化验总有机碳以及第二计算式获取所述分析化验总有机碳计算模型；

所述第二计算式为：

其中，TOC表示分析化验总有机碳，GR表示自然伽马能谱测井总伽马曲线，GR_min表示自然伽马能谱测井总伽马曲线以及去铀伽马曲线重叠时自然伽马能谱测井总伽马曲线的最小刻度，GR_max表示自然伽马能谱测井总伽马曲线以及去铀伽马曲线重叠时自然伽马能谱测井总伽马曲线的最大刻度，KTH表示去铀伽马曲线，KTH_min表示自然伽马能谱测井总伽马曲线以及去铀伽马曲线重叠时去铀伽马曲线的最小刻度，KTH_max表示自然伽马能谱测井总伽马曲线以及去铀伽马曲线重叠时去铀伽马曲线的最大刻度，A表示第三经验参数，B表示第四经验参数，TH表示钍曲线，U表示铀曲线。

可选的，上述密度曲线校正方法中，所述基于所述经验关系式以及所述分析化验总有机碳计算模型确定与所述密度曲线失真段对应的密度校正曲线的步骤，包括：

根据所述分析化验总有机碳计算模型计算所述密度曲线失真段对应的总有机碳含量曲线；

根据所述总有机碳含量曲线以及所述经验关系式确定与所述密度曲线失真段对应的密度校正曲线。

上述实施方式中，在密度曲线失真的井段，根据分析化验总有机碳计算模型进行计算可以得到能反映原状地层总有机碳含量的总有机碳含量曲线，然后将该总有机碳含量曲线代入基于所述经验关系式，即可以准确的得到与密度曲线失真段对应的密度校正曲线。

上述方法步骤的具体实施例过程可参见前述实施例，本实施例在此不再重复赘述。

实施例六

本申请实施例提供了一种电子设备，该电子设备可以是手机、电脑或平板电脑等，包括存储器和处理器，所述存储器上存储有计算器程序，该计算机程序被处理器执行时实现如实施例一中所述的密度曲线校正方法。可以理解，如图7所示，该电子设备700还可以包括：处理器701，存储器702，多媒体组件703，输入/输出(I/O)接口704，以及通信组件705。

其中，处理器701用于执行如实施例一中的密度曲线校正方法中的全部或部分步骤。存储器702用于存储各种类型的数据，这些数据例如可以包括电子设备中的任何应用程序或方法的指令，以及应用程序相关的数据。

处理器701可以是专用集成电路(Application Specific Integrated Circuit，简称ASIC)、数字信号处理器(Digital Signal Processor，简称DSP)、数字信号处理设备(Digital Signal Processing Device，简称DSPD)、可编程逻辑器件(Programmable LogicDevice，简称PLD)、现场可编程门阵列(Field Programmable Gate Array，简称FPGA)、控制器、微控制器、微处理器或其他电子元件实现，用于执行上述实施例一中的密度曲线校正方法。

存储器702可以由任何类型的易失性或非易失性存储设备或者它们的组合实现，例如静态随机存取存储器(Static Random Access Memory，简称SRAM)，电可擦除可编程只读存储器(Electrically Erasable Programmable Read-Only Memory，简称EEPROM)，可擦除可编程只读存储器(Erasable Programmable Read-Only Memory，简称EPROM)，可编程只读存储器(Programmable Read-OnlyMemory，简称PROM)，只读存储器(Read-Only Memory，简称ROM)，磁存储器，快闪存储器，磁盘或光盘。

多媒体组件703可以包括屏幕和音频组件，该屏幕可以是触摸屏，音频组件用于输出和/或输入音频信号。例如，音频组件可以包括一个麦克风，麦克风用于接收外部音频信号。所接收的音频信号可以被进一步存储在存储器或通过通信组件发送。音频组件还包括至少一个扬声器，用于输出音频信号。

I/O接口704为处理器701和其他接口模块之间提供接口，上述其他接口模块可以是键盘，鼠标，按钮等。这些按钮可以是虚拟按钮或者实体按钮。

通信组件705用于该电子设备700与其他设备之间进行有线或无线通信。无线通信，例如Wi-Fi，蓝牙，近场通信(Near Field Communication，简称NFC)，2G、3G或4G，或它们中的一种或几种的组合，因此相应的该通信组件705可以包括：Wi-Fi模块，蓝牙模块，NFC模块。

综上，本申请提供的一种密度曲线校正方法、装置、存储介质及电子设备，该方法包括：获取待校正储层的测井曲线信息，并根据所述测井曲线信息确定出所述测井曲线信息中的密度曲线失真段及密度曲线未失真段；基于所述密度曲线未失真段的密度测井曲线以及总有机碳数据建立经验关系式；获取所述待校正储层基于自然伽马能谱的分析化验总有机碳计算模型；基于所述经验关系式以及所述分析化验总有机碳计算模型确定与所述密度曲线失真段对应的密度校正曲线。首先根据待校正储层的测井曲线信息划分出密度曲线失真段，然后在密度曲线未失真段，分析密度测井曲线与总有机碳数据之间的经验关系，同时获取基于自然伽马能谱的分析化验总有机碳计算模型。然后在整个页岩层段，根据经验关系式以及分析化验总有机碳计算模型进行计算，便可得到反映真实页岩地层的密度曲线，即密度校正曲线。上述方式以总有机碳数据为桥梁，利用自然伽马能谱测井曲线、密度曲线以及总有机碳数据，实现对待校正储层的密度曲线进行准确校正的目的。

在本申请实施例所提供的几个实施例中，应该理解到，所揭露的***和方法，也可以通过其它的方式实现。以上所描述的***和方法实施例仅仅是示意性的。

需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。

虽然本申请所揭露的实施方式如上，但所述的内容只是为了便于理解本申请而采用的实施方式，并非用以限定本申请。任何本申请所属技术领域内的技术人员，在不脱离本申请所揭露的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本申请的专利保护范围，仍须以所附的权利要求书所界定的范围为准。

Claims (10)

1.一种密度曲线校正方法，其特征在于，所述方法包括：

获取待校正储层的测井曲线信息，并根据所述测井曲线信息确定出所述测井曲线信息中的密度曲线失真段及密度曲线未失真段；

基于所述密度曲线未失真段的密度测井曲线以及总有机碳数据建立经验关系式；

获取所述待校正储层基于自然伽马能谱的分析化验总有机碳计算模型；

基于所述经验关系式以及所述分析化验总有机碳计算模型确定与所述密度曲线失真段对应的密度校正曲线。

2.根据权利要求1所述的方法，其特征在于，所述根据所述测井曲线信息确定出所述测井曲线信息中的密度曲线失真段及密度曲线未失真段的步骤，包括：

确定所述测井曲线信息中的密度曲线、去铀伽马曲线以及井径曲线；

将所述密度曲线与所述去铀伽马曲线进行比对，以获取比对结果；

根据所述比对结果以及所述井径曲线的扩径程度确定出所述测井曲线信息的密度曲线失真段及密度曲线未失真段。

3.根据权利要求1所述的方法，其特征在于，所述基于所述密度曲线未失真段的密度测井曲线以及总有机碳数据建立经验关系式的步骤，包括：

基于所述密度曲线未失真段的密度测井曲线获取反映所述待校正储层真实密度信息的密度测井曲线值；

获取总有机碳数据中与所述密度曲线未失真段对应的分析化验总有机碳；

根据所述密度测井曲线值与所述分析化验总有机碳确定经验关系式。

4.根据权利要求3所述的方法，其特征在于，所述根据所述密度测井曲线值与所述分析化验总有机碳确定经验关系式的步骤，包括：

根据所述密度测井曲线值、所述分析化验总有机碳以及第一计算式确定所述经验关系式；

所述经验关系式为：

DEN＝a*TOC+b

其中，DEN表示密度测井曲线值，TOC表示分析化验总有机碳，a表示第一经验参数，b表示第二经验参数。

5.根据权利要求1所述的方法，其特征在于，所述获取所述待校正储层基于自然伽马能谱的分析化验总有机碳计算模型的步骤之前，所述方法还包括：

获取所述待校正储层的自然伽马能谱测井曲线数据和分析化验总有机碳；

根据所述自然伽马能谱测井曲线数据和所述分析化验总有机碳建立分析化验总有机碳计算模型。

6.根据权利要求5所述的方法，其特征在于，所述根据所述自然伽马能谱测井曲线数据和所述分析化验总有机碳建立分析化验总有机碳计算模型的步骤，包括：

确定所述自然伽马能谱测井曲线数据中的自然伽马能谱测井总伽马曲线、去铀伽马曲线、钍曲线和铀曲线；

根据所述自然伽马能谱测井总伽马曲线、所述去铀伽马曲线、钍曲线、铀曲线、所述分析化验总有机碳以及第二计算式获取所述分析化验总有机碳计算模型；

所述第二计算式为：

其中，TOC表示分析化验总有机碳，GR表示自然伽马能谱测井总伽马曲线，GR_min表示自然伽马能谱测井总伽马曲线以及去铀伽马曲线重叠时自然伽马能谱测井总伽马曲线的最小刻度，GR_max表示自然伽马能谱测井总伽马曲线以及去铀伽马曲线重叠时自然伽马能谱测井总伽马曲线的最大刻度，KTH表示去铀伽马曲线，KTH_min表示自然伽马能谱测井总伽马曲线以及去铀伽马曲线重叠时去铀伽马曲线的最小刻度，KTH_max表示自然伽马能谱测井总伽马曲线以及去铀伽马曲线重叠时去铀伽马曲线的最大刻度，A表示第三经验参数，B表示第四经验参数，TH表示钍曲线，U表示铀曲线。

7.根据权利要求1所述的方法，其特征在于，所述基于所述经验关系式以及所述分析化验总有机碳计算模型确定与所述密度曲线失真段对应的密度校正曲线的步骤，包括：

根据所述分析化验总有机碳计算模型计算所述密度曲线失真段对应的总有机碳含量曲线；

根据所述总有机碳含量曲线以及所述经验关系式确定与所述密度曲线失真段对应的密度校正曲线。

8.一种密度曲线校正装置，其特征在于，所述装置包括：

失真段确定模块，用于获取待校正储层的测井曲线信息，并根据所述测井曲线信息确定出所述测井曲线信息中的密度曲线失真段及密度曲线未失真段；

经验关系建立模块，用于基于所述密度曲线未失真段的密度测井曲线以及总有机碳数据建立经验关系式；

分析化验总有机碳计算模型获取模块，用于获取所述待校正储层基于自然伽马能谱的分析化验总有机碳计算模型；

密度校正曲线确定模块，用于基于所述经验关系式以及所述分析化验总有机碳计算模型确定与所述密度曲线失真段对应的密度校正曲线。

9.一种存储介质，其特征在于，该存储介质存储的计算机程序，在被一个或多个处理器执行时，用来实现如权利要求1-7中任意一项所述的方法。

10.一种电子设备，其特征在于，包括存储器和处理器，所述存储器上存储有计算机程序，该计算机程序被所述处理器执行时，执行如权利要求1-7中任意一项所述的方法。

Images