具体实施方式

下面结合附图来具体描述本发明的优选实施例，其中，附图构成本申请一部分，并与本发明的实施例一起用于阐释本发明的原理，并非用于限定本发明的范围。

请参照图1，图1为本发明提供的利用测井资料计算页岩油储层含油量的方法的一实施例的结构示意图，所述利用测井资料计算页岩油储层含油量的方法包括：

S100、获取页岩油储层的干酪根体积分数，本实施例中，通过测井资料得到页岩油储层的干酪根体积分数；

S200、获取页岩油储层的总有机碳(TOC)体积分数，页岩油储层的总有机碳质量分数可通过测井资料获取，此为现有技术，对此不再赘述，再将有机碳质量分数经过转换得到有机碳体积分数；

S300、根据所述页岩油储层的干酪根体积分数及总有机碳体积分数，得到所述页岩油储层的含油质量分数，由于总有机碳由干酪根及储层孔隙中的石油组成，因此，在获取总有机碳的体积分数及干酪根的体积分数后，可得到页岩油储层的含油体积分数，再将页岩油储层的含油体积分数经过转换得到页岩油储层的含油质量分数。

本发明通过页岩油储层的干酪根体积分数及总有机碳体积分数，得到页岩油储层的含油质量分数，由于干酪根体积分数及总有机碳体积分数均可通过测井资料获取，因此，通过本方法可利用测井资料计算页岩油储层含油量，同时，本方法计算过程简单且合乎物理推导论证过程，可靠性良好，有利于促进页岩油储层勘探开发。

本发明中，获取步骤S100中页岩油储层的干酪根体积分数的方法有两种，一种是通过核磁共振测井与密度测井的结合获得，另一种是通过核磁共振测井与中子测井的结合获得，两种方法原理类似，在具体应用时，既可选择二者之一计算页岩油储层的干酪根体积分数，也可以分别通过两种方法计算得到页岩油储层的干酪根体积分数，再取平均值，下面分别对这两种方法进行介绍。

(一)通过核磁共振测井与密度测井的结合获取干酪根体积分数

请参照图2，在一优选的实施例中，所述步骤S100具体包括：

S101、获取页岩储层的核磁共振测井孔隙度值；

S102、获取页岩储层的密度测井值；

S103、获取页岩储层中岩石矿物骨架密度值、流体密度值及干酪根密度值；

S104、根据所述页岩储层的核磁共振测井孔隙度值、密度测井值、岩石矿物骨架密度值、流体密度值及干酪根密度值，得到所述页岩油储层的干酪根体积分数，具体为：

其中，为所述页岩储层的干酪根体积分数，DEN为页岩储层的密度测井值，φ_NMR为核磁共振测井孔隙度值，ρ_ma为页岩储层中岩石矿物骨架密度值，ρ_f为页岩储层中流体密度值，ρ_kero为页岩储层中的干酪根密度值，k₁为第一标定系数。

公式(1)的推导过程如下：

请参照图3，图3为页岩油储层的体积模型(其中，φ_DEN为密度测井孔隙度，φ_N为中子测井孔隙度，φ_NMR为核磁共振测井孔隙度)，根据对页岩油储层地层的分析，地层体积模型的组分可以分为孔隙(Pores)、固体干酪根(Kerogen)和岩石矿物骨架(Mineral Matrix)三部分，从而可得：

V_kero+V_ma+V_pore＝1 (2)

其中V_kero为干酪根的体积分数，V_ma为岩石矿物骨架的体积分数，V_pore为孔隙度，孔隙中存在流体，因此可以被核磁共振测井方法识别，而干酪根对核磁共振测井的响应贡献较小，可以忽略。同时，因为干酪根密度及含氢指数与流体接近，其密度测井响应及中子测井响应与孔隙流体较为接近。因此，其密度测井的响应方程可以写作以下形式：

DEN＝V_kero·ρ_kero+V_ma·ρ_ma+V_pore·ρ_f (3)

其中，DEN是密度测井值，ρ_ma为页岩储层中岩石矿物骨架密度值，ρ_f为页岩储层中流体密度值，ρ_kero为页岩储层中的干酪根密度值，V_kero为干酪根体积分数，V_ma为岩石矿物骨架的体积分数，V_pore为孔隙度。

同时，核磁共振测井得到的孔隙度φ_NMR与公式(3)中的v_pore相等，即：

φ_NMR＝V_pore (4)

将公式(2)及公式(4)代入公式(3)，可得：

由于在参数设定等环节可能存在一定的误差，因此，引入第一标定系数k₁来对公式(5)进行修正，即可得到公式(1)，需要指出的是：公式(1)中的与公式(5)中的V_kero的物理含义相同。

具体标定方法是：利用公式(5)通过测井测量得到V_kero的值，再对相同的储层(同一口井，同一深度)通过岩心测量得到从而可计算得到第一标定系数k₁的值。

(二)通过核磁共振测井与中子测井的结合获取干酪根体积分数

请参照图4，在一优选的实施例中，所述步骤S100具体包括：

S111、获取页岩储层的核磁共振测井孔隙度值；

S112、获取页岩储层的中子测井孔隙度值；

S113、根据所述页岩储层的核磁共振测井孔隙度值及中子测井孔隙度值得到所述页岩油储层的干酪根体积分数，具体为：

其中，为所述页岩储层的干酪根体积分数，φ_N为密度测井孔隙度值，φ_NMR为核磁共振测井孔隙度值，k₂为第二标定系数。

公式(6)的推导过程如下：

由于孔隙中存在流体，因此可以被核磁共振测井方法识别，而干酪根对核磁共振测井的响应贡献较小，可以忽略。同时，因为干酪根密度及含氢指数与流体接近，其中子测井响应与孔隙流体较为接近，因此，通过中子测井方法获得的孔隙度实际上等于干酪根的体积分数与真实的孔隙度之和，而核磁共振测井得到的孔隙度则为真实的孔隙度，因此，核磁共振测井计算的孔隙度和密度测井计算的孔隙度的差值等于干酪根的体积分数，类似的，为减小误差，引入第二标定系数k₂来对公式进行修正，即可得到公式(6)，第二标定系数k₂的确定方法与第一标定系数k₁类似，对此不再赘述。

本发明中，请参照图5，获取步骤S200中页岩油储层的总有机碳体积分数的方法，具体包括如下步骤：

S201、获取页岩油储层的总有机碳质量分数，目前利用测井资料计算总有机碳质量分数的有效方法有很多种，例如△logR法，双泥质含量法(dual-V_sh method)等，因此在本发明中不再加以赘述，计算总有机碳质量分数的方法可根据资料情况和需求进行选择；

S202、获取页岩油储层的有机质的平均密度；

S203、根据所述页岩油储层的总有机碳质量分数及有机质的平均密度，得到页岩油储层的总有机碳体积分数，具体为：

其中，V_TOC为页岩油储层的总有机碳体积分数，TOC为页岩油储层的总有机碳质量分数，DEN为页岩油储层的密度测井值，ρ_TOC为页岩油储层的有机质的平均密度，其中，ρ_TOC可根据实验数据或者利用经验进行该值的选择。

本发明中，请参照图6，步骤S300中根据所述页岩油储层的干酪根体积分数及总有机碳体积分数，得到所述页岩油储层的含油质量分数，具体包括如下步骤：

S301、根据所述页岩油储层的干酪根体积分数及总有机碳体积分数，获取所述页岩油储层的含油体积分数，具体为：

V_o＝V_TOC-V_kero (8)

其中，V_o为页岩油储层的含油体积分数，V_TOC为页岩油储层的干酪根体积分数，V_kero为页岩油储层的干酪根体积分数。

S302、获取所述页岩油储层内石油的密度；

S303、根据所述页岩油储层的含油体积分数及所述页岩油储层内石油的密度，得到所述页岩油储层的含油质量分数，具体为：

其中C_o为页岩油储层的含油质量分数；ρ_o为油的密度,可根据实验数据或者利用经验进行该值的选择，V_o为页岩油储层的含油体积分数，ρ_o为页岩油储层内石油的密度，DEN为页岩油储层的密度测井值。

请参照图8，本发明还提供了一种电子设备，所述电子设备可以是移动终端、桌上型计算机、笔记本、掌上电脑及服务器等计算设备。该电子设备包括处理器10、存储器20及显示器30。图8仅示出了该电子设备的部分组件，但是应理解的是，并不要求实施所有示出的组件，可以替代的实施更多或者更少的组件。

所述存储器20在一些实施例中可以是所述电子设备的内部存储单元，例如电子设备的硬盘或内存。所述存储器20在另一些实施例中也可以是所述电子设备的外部存储设备，例如所述电子设备上配备的插接式硬盘、智能存储卡(Smart Media Card,SMC)、安全数字(Secure Digital,SD)卡、闪存卡(Flash Card)等。进一步地，所述存储器20还可以既包括电子设备的内部存储单元也包括外部存储设备。所述存储器20用于存储安装于所述电子设备的应用软件及各类数据，例如所述安装电子设备的程序代码等。所述存储器20还可以用于暂时地存储已经输出或者将要输出的数据。在一实施例中，存储器20上存储有利用测井资料计算页岩油储层含油量的方法的程序40，该利用测井资料计算页岩油储层含油量的方法的程序40可被处理器10所执行，从而实现本申请各实施例的利用测井资料计算页岩油储层含油量的方法。

所述处理器10在一些实施例中可以是一中央处理器(Central Processing Unit,CPU)、微处理器或其他数据处理芯片，用于运行所述存储器20中存储的程序代码或处理数据，例如执行所述利用测井资料计算页岩油储层含油量的方法等。

所述显示器30在一些实施例中可以是LED显示器、液晶显示器、触控式液晶显示器以及OLED(Organic Light-Emitting Diode，有机发光二极管)触摸器等。所述显示器30用于显示在所述电子设备的信息以及用于显示可视化的用户界面。所述电子设备的部件10-30通过系统总线相互通信。

在一实施例中，当处理器10执行所述存储器20中利用测井资料计算页岩油储层含油量的方法的程序40时实现如上述实施例所述的利用测井资料计算页岩油储层含油量的方法中的步骤，由于上文已对利用测井资料计算页岩油储层含油量的方法进行详细描述，在此不再赘述。

本发明还提供了一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器进行加载，以本发明提供的利用测井资料计算页岩油储层含油量的方法中的步骤。

为了验证本发明提供的技术方案的准确性及可靠性，选取一页岩油区块进行研究，如图7所示，利用本发明提出的方法，依次计算每个深度点的TOC和干酪根含量，最终可以获得连续井段的含油量值。图7中右数第2道为干酪根计算含量，最后一道为含油量计算结果与岩心结果的比较，曲线为计算值，棒状为岩心值，计算结果显示通过本发明提供技术方案计算得到的含油质量分数数据点与通过岩心资料获取的含油质量分数数据点吻合较好，证明了本发明提供的技术方案准确性及可靠性较高，能满足实际勘探开发的需求。

本发明提供的技术方案的有益效果包括：

(1)通过页岩油储层的干酪根体积分数及总有机碳体积分数，得到页岩油储层的含油质量分数，由于干酪根体积分数及总有机碳体积分数均可通过测井资料获取，因此，通过本方法可利用测井资料计算页岩油储层含油量，同时，本方法计算过程简单且合乎物理推导论证过程，可靠性良好，有利于促进页岩油储层勘探开发；

(2)通过核磁共振测井与密度测井的结合获取干酪根体积分数，由于核磁共振测井对干酪根不响应，因此，核磁共振测量得到的孔隙度为页岩油储层的真实孔隙度，同时，密度测井对干酪根响应，将核磁共振测量得到的真实孔隙度代入密度测井计算公式中，即可得到本发明中计算干酪根体积分数的第一种方法；

(3)通过核磁共振测井与中子测井的结合获取干酪根体积分数，由于核磁共振测量得到的孔隙度为页岩油储层的真实孔隙度，而中子测井方法获得的孔隙度等于干酪根的体积分数与真实的孔隙度之和，因此，核磁共振测井计算的孔隙度和密度测井计算的孔隙度的差值等于干酪根的体积分数，即可得到本发明中计算干酪根体积分数的第二种方法；

(4)本发明中，利用岩心测试获取得到干酪根体积分数对通过测井方法得到的干酪根体积分数进行标定，从而可提高干酪根体积分数的计算结果的准确性。

以上所述仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。