基于页岩油赋存状态表征的页岩游离油饱和度预测方法与流程

1.本技术涉及石油勘探技术，尤其涉及一种基于页岩油赋存状态表征的页岩游离油饱和度预测方法。


背景技术：

2.在开采页岩中的石油之前，判断待开采地区的页岩是否满足开采条件是至关重要的。其中，可动流体饱和度是判断页岩是否值得开采的重要依据。上述可动流体饱和度主要是指页岩的游离油饱和度。页岩的游离油饱和度可以通过该页岩的岩石样品(简称岩心)的游离油饱和度类比得到。岩心的游离油饱和度越高，对应的页岩的游离油饱和度越高，说明该页岩中可开采的石油越多，满足开采条件的可能性越大。岩心的游离油饱和度越低，对应的页岩的游离油饱和度越低，说明该页岩中可开采的石油越少，满足开采条件的可能性越小。因此，获取准确的页岩岩心的游离油饱和度是至关重要的。
3.目前，现有的对页岩岩心的游离油饱和度进行预测的方法主要是基于阿尔奇公式的预测方法。具体的，在使用基于阿尔奇公式对页岩岩心的游离油饱和度进行预测时，需要将页岩岩心的电阻率作为阿尔奇公式的其中一个输入项，以获取该岩心的游离油饱和度预测结果。
4.然而，页岩的粘土矿物含量较多，导致页岩的电阻率低于实际电阻率，进而可能导致获取的页岩岩心的游离油饱和度的准确性较差。


技术实现要素：

5.本技术提供一种基于页岩油赋存状态表征的页岩游离油饱和度预测方法，以解决预测页岩岩心的游离油饱和度的准确性较差的问题。
6.第一方面，本技术提供一种基于页岩油赋存状态表征的页岩游离油饱和度预测模型的构建方法，所述页岩油赋存状态表征指的是页岩在各横向弛豫时间区间对应的游离油饱和度，所述方法包括：
7.获取n个来自所述页岩的样本岩心在m个横向弛豫时间区间对应的孔隙度；所述样本岩心在m个横向弛豫时间区间对应的孔隙度的和为所述样本岩心的总孔隙度，所述n和所述m均为大于或等于2的整数、且所述n大于或等于所述m；
8.获取n个所述样本岩心的页岩油赋存状态表征；
9.根据所述n个所述样本岩心的页岩油赋存状态表征，获取所述n个所述样本岩心的游离油饱和度；
10.根据n个所述样本岩心在m个横向弛豫时间区间对应的孔隙度，以及，n个所述样本岩心的游离油饱和度，确定每个所述横向弛豫时间区间对应的权重；
11.根据每个所述横向弛豫时间区间对应的权重，构建所述页岩游离油饱和度预测模型；所述页岩游离油饱和度预测模型用于根据待测岩心在m个横向弛豫时间区间对应的孔隙度，获取所述待测岩心的游离油饱和度；所述待测岩心所在的页岩的游离油饱和度等于
所述待测岩心的游离油饱和度。
12.可选的，所述根据n个所述样本岩心在m个横向弛豫时间区间对应的孔隙度，以及，n个所述样本岩心的游离油饱和度，确定每个所述横向弛豫时间区间对应的权重，包括：
13.构建每个所述样本岩心对应的求解权重方程；所述样本岩心对应的求解权重方程为：
[0014][0015]
其中，所述ω
i
表示第i个横向弛豫时间区间对应的权重，所述bin
i
表示所述样本岩心在所述第i个横向弛豫时间区间对应的孔隙度，所述φ表示所述样本岩心的总孔隙度，所述p表示所述样本岩心的游离油饱和度；i为大于或等于1，且小于或等于m的整数；
[0016]
根据n个所述样本岩心对应的求解权重方程，求解各所述横向弛豫时间区间对应的权重的取值。
[0017]
可选的，所述根据n个所述样本岩心对应的求解权重方程，求解各所述横向弛豫时间区间对应的权重的取值，包括：
[0018]
将所述n个样本岩心划分为至少一个样本岩心组，每个样本岩心组包括m个样本岩心，任意两个样本岩心组至少有一个样本岩心不同；
[0019]
根据每组样本岩心对应的求解权重方程，获取每组样本岩心在各横向弛豫时间区间对应的初始权重；
[0020]
将所有组样本岩心在同一横向弛豫时间区间对应的初始权重的平均值，作为该横向弛豫时间区间对应的权重。
[0021]
可选的，获取n个所述样本岩心的页岩油赋存状态表征，包括：
[0022]
针对n个所述样本岩心中的任一样本岩心，获取所述样本岩心的横向弛豫时间谱，以及，所述样本岩心对应的饱和油样本岩心的横向弛豫时间谱；所述饱和油样本岩心为对所述样本岩心进行抽真空，以及，饱和页岩原油得到的；所述样本岩心的横向弛豫时间谱用于表征所述样本岩心在各横向弛豫时间区间对应的束缚油饱和度，所述样本岩心对应的饱和油样本岩心的横向弛豫时间谱用于表征所述样本岩心在各横向弛豫时间区间对应的束缚油饱和度与游离油饱和度的和；
[0023]
根据所述样本岩心的横向弛豫时间谱，以及，所述样本岩心对应的饱和油样本岩心的横向弛豫时间谱，获取所述样本岩心的页岩油赋存状态表征。
[0024]
第二方面，本技术提供一种页岩游离油饱和度预测方法，所述方法包括：
[0025]
获取待测岩心的各横向弛豫时间区间对应的孔隙度，所述待测岩心为从所述页岩开采得到的；
[0026]
根据所述待测岩心的各横向弛豫时间区间对应的孔隙度，利用如第一方面任一项所述的页岩游离油饱和度预测模型，获取所述待测岩心的游离油饱和度；所述待测岩心所在的页岩的游离油饱和度等于所述待测岩心的游离油饱和度；
[0027]
输出所述页岩的游离油饱和度。
[0028]
可选的，在所述获取所述待测岩心的各横向弛豫时间区间对应的孔隙度之前，还包括：
[0029]
确定所述待测岩心中存在游离油。
[0030]
可选的，在所述获取所述待测岩心的游离油饱和度之后，还包括：
[0031]
根据所述待测岩心的游离油饱和度，确定所述页岩是否满足开采条件；
[0032]
输出提示信息，所述提示信息包括所述确定所述页岩是否满足开采条件的结果。
[0033]
第三方面，本技术提供一种基于页岩油赋存状态表征的页岩游离油饱和度预测模型的构建装置，所述页岩油赋存状态表征指的是页岩在各横向弛豫时间区间对应的游离油饱和度，所述装置包括：
[0034]
第一获取模块，用于获取n个来自所述页岩的样本岩心在m个横向弛豫时间区间对应的孔隙度；所述样本岩心在m个横向弛豫时间区间对应的孔隙度的和为所述样本岩心的总孔隙度，所述n和所述m均为大于或等于2的整数、且所述n大于或等于所述m；
[0035]
第二获取模块，用于获取n个所述样本岩心的页岩油赋存状态表征；
[0036]
第三获取模块，用于根据所述n个所述样本岩心的页岩油赋存状态表征，获取所述n个所述样本岩心的游离油饱和度；
[0037]
处理模块，用于根据n个所述样本岩心在m个横向弛豫时间区间对应的孔隙度，以及，n个所述样本岩心的游离油饱和度，确定每个所述横向弛豫时间区间对应的权重；
[0038]
构建模块，用于根据每个所述横向弛豫时间区间对应的权重，构建所述页岩游离油饱和度预测模型；所述页岩游离油饱和度预测模型用于根据待测岩心在m个横向弛豫时间区间对应的孔隙度，获取所述待测岩心的游离油饱和度；所述待测岩心所在的页岩的游离油饱和度等于所述待测岩心的游离油饱和度。
[0039]
第四方面，本技术提供一种页岩游离油饱和度预测装置，所述装置包括：
[0040]
第一获取模块，用于获取待测岩心的各横向弛豫时间区间对应的孔隙度，所述待测岩心为从所述页岩开采得到的；
[0041]
第二获取模块，用于根据所述待测岩心的各横向弛豫时间区间对应的孔隙度，利用如第一方面任一项所述的页岩游离油饱和度预测模型，获取所述待测岩心的游离油饱和度；所述待测岩心所在的页岩的游离油饱和度等于所述待测岩心的游离油饱和度；
[0042]
输出模块，用于输出所述页岩的游离油饱和度。
[0043]
第五方面，本技术提供一种电子设备，电子设备应用于井下页岩核磁共振测井数据的处理与解释，包括：至少一个处理器、存储器；
[0044]
所述存储器存储计算机执行指令；
[0045]
所述至少一个处理器执行所述存储器存储的计算机执行指令，使得所述电子设备执行第一方面或第二方面任一项所述的方法。
[0046]
本技术提供的基于页岩油赋存状态表征的页岩游离油饱和度预测方法，通过n个样本岩心在m个横向弛豫时间区间对应的孔隙度，以及，n个样本岩心的页岩油赋存状态表征，确定每个横向弛豫时间区间对应的权重。然后基于上述每个横向弛豫时间区间对应的权重，构建页岩游离油饱和度预测模型。通过上述方法得到的页岩游离油饱和度预测模型，不需要基于岩心的电阻率来构建页岩游离油饱和度预测模型，提高了构建页岩游离油饱和度预测模型的准确性，进而提高了使用该模型获取待测岩心的游离油饱和度的准确性，因此提高了对页岩的游离油饱和度进行表征的准确性。
附图说明
[0047]
为了更清楚地说明本技术或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图做一简单地介绍，显而易见地，下面描述中的附图是本技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
[0048]
图1为岩心的横向弛豫时间谱的示意图；
[0049]
图2为本技术提供的一种基于页岩油赋存状态表征的页岩游离油饱和度预测模型的构建方法的流程示意图；
[0050]
图3为本技术提供的一种页岩游离油饱和度预测方法的流程示意图；
[0051]
图4为本技术提供的一种基于页岩油赋存状态表征的页岩游离油饱和度预测模型的构建装置的结构示意图；
[0052]
图5为本技术提供的一种页岩游离油饱和度预测装置的结构示意图；
[0053]
图6为本技术提供的一种电子设备结构示意图。
[0054]
通过上述附图，已示出本技术明确的实施例，后文中将有更详细的描述。这些附图和文字描述并不是为了通过任何方式限制本技术构思的范围，而是通过参考特定实施例为本领域技术人员说明本技术的概念。
具体实施方式
[0055]
为使本技术的目的、技术方案和优点更加清楚，下面将结合本技术中的附图，对本技术中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本技术一部分实施例，而不是全部的实施例。基于本技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本技术保护的范围。
[0056]
下面对本技术涉及到的名词概念进行解释：
[0057]
页岩油：指的是以页岩为主的页岩层系中所含的石油。
[0058]
页岩：指的是具有连通孔隙，且允许石油、天然气等在其中储存和渗滤的页岩岩层。
[0059]
束缚油：通过石油开采，从页岩中无法开采出的石油。
[0060]
游离油：通过石油开采，从页岩中可以开采出的石油，也可以称为可流动油(或者可动油)。
[0061]
总孔隙度：指的是岩心中所有的孔隙空间体积之和与岩心总体积的比值。
[0062]
横向弛豫时间：通过核磁共振检测设备，对岩心进行核磁共振检测实验，可以得到该岩心的横向弛豫时间。岩心的横向弛豫时间可以反映岩心中存在的物质种类(束缚油或者游离油)。示例性的，以3个岩心为例，图1为岩心的横向弛豫时间谱的示意图。如图1所示，以岩心1、岩心2，或者，岩心3中的任意一个岩心对应的孔隙度分量曲线为例，在横向弛豫时间0毫秒(ms)到35ms之间的曲线与横轴(时间轴)围成的面积，表示该岩心中的束缚油所占该岩心的孔隙分量。在横向弛豫时间35毫秒之后的曲线与横轴(时间轴)围成的面积，表示该岩心中的游离油所占该岩心的孔隙分量。
[0063]
对于待开采地区的页岩，页岩的游离油饱和度越高，说明该页岩中可开采的石油越多，该待开采地区的页岩满足开采条件的可能性越大。页岩的游离油饱和度越低，说明该
页岩中可开采的石油越少，该待开采地区的页岩满足开采条件的可能性越小。因此，获取准确的页岩的游离油饱和度是至关重要的。
[0064]
通常，页岩的游离油饱和度可以通过该页岩的岩石样品(简称岩心)的游离油饱和度类比得到。也就是说，可以通过研究页岩岩心的游离油饱和度，确定该待开采地区的页岩是否满足开采条件。
[0065]
目前，现有的对页岩岩心的游离油饱和度进行预测的方法主要是基于阿尔奇公式的预测方法。具体的，在使用基于阿尔奇公式对页岩岩心的游离油饱和度进行预测时，需要将页岩岩心的电阻率作为阿尔奇公式的其中一个输入项，以获取该岩心的游离油饱和度的预测结果。
[0066]
然而，页岩的粘土矿物含量较多，导致页岩的电阻率低于实际电阻率，进而可能导致获取的页岩岩心的游离油饱和度的准确性较差。
[0067]
考虑到现有的页岩游离油饱和度预测方法存在准确性较低的问题是页岩的电阻率低于实际电阻率导致的，本技术提出了一种基于页岩游离油饱和度预测模型，以及，岩心在多个横向弛豫时间区间对应的孔隙度，获取页岩游离油饱和度的方法。其中，页岩游离油饱和度预测模型是基于页岩油赋存状态表征和多个横向弛豫时间区间对应的权重确定的。页岩油赋存状态表征指的是页岩在各横向弛豫时间区间对应的游离油饱和度。岩心中的游离油只存储在岩心的孔隙中，与岩心中除了孔隙之外的结构无关，且现有技术已经能够获取准确性较高的岩心的孔隙度。因此，根据岩心的孔隙度确定岩心的游离油饱和度，保证了确定岩心所属页岩的游离油饱和度的准确性。
[0068]
此外，通过上述方法，不需要使用岩心的电阻率来获取岩心的游离油饱和度，避免了页岩岩心的电阻率不准确带来的误差，提高了对页岩的游离油饱和度进行预测的准确性。
[0069]
下面结合具体地实施例对本技术的提供的页岩游离油饱和度预测模型的构建方法进行详细说明。具体实现时，该方法的执行主体例如可以为终端，或者，服务器等具有处理能力的电子设备。下面这几个具体的实施例可以相互结合，对于相同或相似的概念或过程可能在某些实施例不再赘述。
[0070]
图2为本技术提供的一种基于页岩油赋存状态表征的页岩游离油饱和度预测模型的构建方法的流程示意图。如图2所示，该方法包括以下步骤：
[0071]
s101、获取n个来自页岩的样本岩心在m个横向弛豫时间区间对应的孔隙度。n和m均为大于或等于2的整数，且n大于或等于m。
[0072]
上述n个样本岩心可以是从同一个测井开采得到的岩心，也可以是从不同测井开采得到的岩心。此外，上述n个样本岩心所在页岩的深度(或者称为井段)可以相同，也可以不同。
[0073]
上述m个横向弛豫时间区间指的是将一段横向弛豫时间的总时长划分为m段，得到的m个区间。示例性的，以上述一段横向弛豫时间总时长为32毫秒(ms)，m等于8为例，该8个横向弛豫时间区间对应的区间范围例如可以如下表1所示：
[0074]
表1
[0075]
横向弛豫时间区间区间范围区间1(0ms，4ms)
区间2(4ms，16ms)区间3(16ms，32ms)区间4(32ms，64ms)区间5(64ms，128ms)区间6(128ms，256ms)区间7(256ms，512ms)区间8(512ms，1024ms)
[0076]
应理解，表1仅是对各个横向弛豫时间区间对应的区间范围的示例。具体实现时，不同横向弛豫时间区间对应的时长可以不同，也可以相同。
[0077]
样本岩心在上述m个横向弛豫时间区间对应的孔隙度的和为该样本岩心的总孔隙度。示例性的，仍然以m等于8为例，样本岩心在8个横向弛豫时间区间对应的孔隙度例如可以如下表2所示：
[0078]
表2
[0079]
横向弛豫时间区间孔隙度区间1孔隙度1区间2孔隙度2区间3孔隙度3区间4孔隙度4区间5孔隙度5区间6孔隙度6区间7孔隙度7区间8孔隙度8
[0080]
上述孔隙度1到孔隙度8的和等于该样本岩心的总孔隙度。
[0081]
作为一种可能的实现方式，电子设备例如可以根据该样本岩心所在井段的测井数据，获取该样本岩心在m个横向弛豫时间区间对应的孔隙度。其中，上述测井数据例如可以包括通过核磁共振检测设备对该样本岩心进行核磁共振检测得到的测井数据。示例性的，电子设备例如可以通过应用程序接口(application programming interface，api)，或者，图形用户界面(graphical user interface，gui)获取样本岩心所在井段的测井数据。然后，电子设备例如可以通过techlog(一种测井资料处理软件的名称)软件对该测井数据进行处理，以获取该样本岩心在m个横向弛豫时间区间对应的孔隙度。
[0082]
s102、获取n个样本岩心的页岩油赋存状态表征。
[0083]
作为一种可能的实现方式，针对n个样本岩心中的任一样本岩心，电子设备例如可以根据该样本岩心和饱和油样本岩心的横向弛豫时间谱，确定该样本岩心的页岩油赋存状态表征。其中，上述饱和油样本岩心为对该样本岩心进行抽真空，以及，饱和页岩原油得到的。
[0084]
具体的，电子设备首先可以获取样本岩心的横向弛豫时间谱，以及，样本岩心对应的饱和油样本岩心的横向弛豫时间谱。其中，样本岩心的横向弛豫时间谱用于表征样本岩心在各横向弛豫时间区间对应的束缚油饱和度。样本岩心对应的饱和油样本岩心的横向弛豫时间谱用于表征该样本岩心在各横向弛豫时间区间对应的束缚油饱和度与游离油饱和
度的和。
[0085]
可选的，核磁共振检测设备可以对上述样本岩心和饱和油样本岩心进行核磁共振检测实验，得到样本岩心的横向弛豫时间谱，以及，样本岩心对应的饱和油样本岩心的横向弛豫时间谱。然后，电子设备可以接收该核磁共振检测设备发送的样本岩心的横向弛豫时间谱，以及，样本岩心对应的饱和油样本岩心的横向弛豫时间谱。或者，电子设备例如还可以通过api或者gui等获取用户输入的样本岩心的横向弛豫时间谱，以及，样本岩心对应的饱和油样本岩心的横向弛豫时间谱。
[0086]
然后，电子设备可以根据上述样本岩心的横向弛豫时间谱，以及，样本岩心对应的饱和油样本岩心的横向弛豫时间谱，获取样本岩心的页岩油赋存状态表征。
[0087]
具体的，针对任意一个横向弛豫时间区间，电子设备可以将该横向弛豫时间区间对应的束缚油饱和度与游离油饱和度的和，减去，该横向弛豫时间区间对应的束缚油饱和度，得到该横向弛豫时间区间对应的游离油饱和度。通过上述方法，电子设备可以获取该样本岩心在各横向弛豫时间区间对应的游离油饱和度，进而得到该样本岩心的页岩油赋存状态表征。
[0088]
作为另一种可能的实现方式，在获取上述n个样本岩心的游离油饱和度之前，还可以对上述n个样本岩心进行切割、制样、烘干等常规预处理。也就是说，上述n个样本岩心还可以是经过常规预处理之后的样本岩心。
[0089]
s103、获取n个样本岩心的游离油饱和度。
[0090]
样本岩心的游离油饱和度等于样本岩心在各横向弛豫时间区间对应的游离油饱和度的和。可选的，针对任意一个样本岩心，在获取该样本岩心的页岩油赋存状态表征之后，可以将该样本岩心在各横向弛豫时间区间对应的游离油饱和度，作为该样本岩心的游离油饱和度。
[0091]
s104、根据n个样本岩心在m个横向弛豫时间区间对应的孔隙度，以及，n个样本岩心的游离油饱和度，确定每个横向弛豫时间区间对应的权重。
[0092]
作为第一种可能的实现方式，当上述n等于m时，电子设备可以直接根据n个样本岩心的m个孔隙度，以及，该n个样本岩心的游离油饱和度，确定每个横向弛豫时间区间对应的权重。
[0093]
作为第二种可能的实现方式，当上述n大于m时，电子设备可以基于该n个样本岩心中的任意m个样本岩心的m个孔隙度，以及，该m个样本岩心的游离油饱和度，确定每个横向弛豫时间区间对应的权重。
[0094]
作为第三种可能的实现方式，当上述n大于m时，可以将该n个样本岩心分为多组样本岩心，其中每组样本岩心的数量为m。然后，电子设备可以基于上述多组样本岩心中各样本岩心的m个孔隙度，以及，该多组样本岩心中各样本岩心的游离油饱和度，确定每个横向弛豫时间区间对应的权重。
[0095]
s105、根据每个横向弛豫时间区间对应的权重，构建页岩游离油饱和度预测模型。
[0096]
其中，上述页岩游离油饱和度预测模型用于根据待测岩心在m个横向弛豫时间区间对应的孔隙度，获取待测岩心的游离油饱和度。该待测岩心所在的页岩的游离油饱和度等于该待测岩心的游离油饱和度。
[0097]
示例性的，上述页岩游离油饱和度预测模型例如可以通过下述公式(1)表示：
[0098][0099]
其中，ω
i
上述m个横向弛豫时间区间中第i个横向弛豫时间区间对应的权重。bin
i
表示样本岩心在该第i个横向弛豫时间区间对应的孔隙度。φ表示该样本岩心的总孔隙度。p表示该样本岩心的游离油饱和度。i为大于或等于1，且小于或等于m的整数。
[0100]
在本实施例中，通过n个样本岩心在m个横向弛豫时间区间对应的孔隙度，以及，n个样本岩心的页岩油赋存状态表征，确定每个横向弛豫时间区间对应的权重。然后基于上述每个横向弛豫时间区间对应的权重，构建页岩游离油饱和度预测模型。通过上述方法得到的页岩游离油饱和度预测模型，不需要基于岩心的电阻率来构建页岩游离油饱和度预测模型，提高了构建页岩游离油饱和度预测模型的准确性，进而提高了使用该模型获取待测岩心的游离油饱和度的准确性，因此提高了对页岩的游离油饱和度进行表征的准确性。
[0101]
下面对电子设备如何根据n个样本岩心在m个横向弛豫时间区间对应的孔隙度，以及，n个样本岩心的游离油饱和度，确定每个横向弛豫时间区间对应的权重进行详细说明。作为一种可能的实现方式，上述步骤s104可以包括以下步骤：
[0102]
步骤1、构建每个样本岩心对应的求解权重方程。
[0103]
针对任意一个样本岩心，上述样本岩心对应的求解权重方程为：
[0104][0105]
其中，该求解权重方程中各字母表示的意义与上述公式(1)相同。在该求解权重方程中，bin
i
、φ，以及，p均为已知数，ω
i
为未知数。以上述m等于8为例，此时，该求解权重方程中包括8个未知数。
[0106]
针对上述n个样本岩心中的任意一个岩心，均构建上述求解权重方程。示例性的，以n等于8为例，电子设备可以构建8个上述求解权重方程。
[0107]
步骤2、根据n个样本岩心对应的求解权重方程，求解各横向弛豫时间区间对应的权重的取值。
[0108]
作为一种可能的实现方式，当n等于m时，电子设备可以构建m个样本岩心对应的求解权重方程。然后，电子设备可以将根据该n个样本岩心对应的求解权重方程，求解得到的各横向弛豫时间区间对应的权重的取值，作为每个横向弛豫时间区间对应的权重。
[0109]
作为另一种可能的实现方式，当n大于m时，电子设备可以基于多组样本岩心组，获取各个横向弛豫时间区间对应的权重，以提高确定横向弛豫时间区间对应的权重的准确性，进而提高上述页岩游离油饱和度预测模型的准确性。
[0110]
具体的，电子设备首先可以将n个样本岩心划分为至少一个样本岩心组。其中，每个样本岩心组包括m个样本岩心，且任意两个样本岩心组至少有一个样本岩心不同。然后，电子设备可以根据每组样本岩心对应的求解权重方程，获取每组样本岩心在各横向弛豫时间区间对应的初始权重。
[0111]
在获取各组样本岩心在各横向弛豫时间区间对应的初始权重之后，电子设备可以将所有组样本岩心在同一横向弛豫时间区间对应的初始权重的平均值，作为该横向弛豫时间区间对应的权重。
[0112]
示例性的，以上述n等于10，m等于8为例，假设电子设备可以将该10个样本岩心划
分为3组样本岩心组。其中每个样本岩心组包括8个样本岩心，且任意两个样本岩心组至少有一个样本岩心不同。示例性的，上述3组样本岩心组在各横向弛豫时间区间对应的初始权重例如可以如下表3所示：
[0113]
表3
[0114][0115]
以表3所示内容为例，示例性的，电子设备可以将ω
11
、ω
21
和ω
31
三者的平均值，作为该横向弛豫时间区间1对应的权重，将ω
12
、ω
22
和ω
33
三者的平均值，作为该横向弛豫时间区间2对应的权重。参照上述方式，电子设备可以获取每个横向弛豫时间区间对应的权重。
[0116]
在确定每个横向弛豫时间区间对应的权重之后，电子设备可以根据每个横向弛豫时间区间对应的权重，构建页岩游离油饱和度预测模型。
[0117]
在得到上述页岩游离油饱和度预测模型之后，可以使用该页岩游离油饱和度预测模型，对页岩的游离油饱和度进行预测。下面结合具体地实施例对使用该页岩游离油饱和度预测模型，对页岩的游离油饱和度进行预测的方法进行详细说明。下面这几个具体的实施例可以相互结合，对于相同或相似的概念或过程可能在某些实施例不再赘述。
[0118]
具体实现时，该页岩游离油饱和度预测方法的执行主体也可以为终端或者服务器等具有处理能力的电子设备。应理解，该执行页岩游离油饱和度预测方法的电子设备，与，执行上述页岩游离油饱和度预测模型构建方法的电子设备，可以是同一设备，也可以是不同设备。
[0119]
图3为本技术提供的一种页岩游离油饱和度预测方法的流程示意图。如图3所示，该方法包括以下步骤：
[0120]
s201、获取待测岩心的各横向弛豫时间区间对应的孔隙度，该待测岩心为从页岩开采得到的。
[0121]
示例性的，电子设备可以获取用户输入的待测岩心的各横向弛豫时间区间对应的孔隙度。或者，电子设备还可以获取核磁共振检测设备发送的待测岩心的各横向弛豫时间区间对应的孔隙度。
[0122]
s202、根据待测岩心的各横向弛豫时间区间对应的孔隙度，利用如前述任一实施例提供的页岩游离油饱和度预测模型，获取待测岩心的游离油饱和度。
[0123]
以该游离油饱和度预测模型可以通过上述公式(1)表示为例，电子设备首先可以获取该待测岩心的总孔隙度。然后将该待测岩心的各横向弛豫时间区间对应的孔隙度、每个横向弛豫时间区间对应的权重，以及，该待测岩心的总孔隙度带入上述公式(1)，以获取该待测岩心的游离油饱和度。
[0124]
其中，该待测岩心所在的页岩的游离油饱和度等于该待测岩心的游离油饱和度。
[0125]
s203、输出页岩的游离油饱和度。
[0126]
示例性的，以该电子设备包括显示装置为例，电子设备在获取页岩的游离油饱和度之后，可以将该页岩的游离油饱和度输出至电子设备的显示装置。或者，电子设备还可以将该页岩的游离油饱和度输出至与该电子设备连接的显示装置。或者，电子设备还可以将该页岩的游离油饱和度发送给用户使用的终端设备。再或者，电子设备例如还可以将页岩的游离油饱和度输出至语音播报装置。
[0127]
在本实施例中，根据待测岩心的各横向弛豫时间区间对应的孔隙度，使用前述实施例构建的页岩游离油饱和度预测模型，即可获取待测岩心的游离油饱和度。该待测岩心所在的页岩的游离油饱和度等于该待测岩心的游离油饱和度。通过该方法，不需要使用岩心的电阻率来获取岩心的游离油饱和度，避免了页岩岩心的电阻率不准确带来的误差，提高了对岩心的游离油饱和度进行预测的准确性。也就是说，通过上述方法提高了确定待测岩心所在的页岩的游离油饱和度的准确性，进而可以提高确定该岩心所属页岩的游离油饱和度的准确性，因此提高了对该页岩是否满足开采条件进行判断的准确性。
[0128]
作为一种可能的实现方式，在获取待测岩心的各横向弛豫时间区间对应的孔隙度之前，电子设备还可以确定待测岩心中是否存在游离油，并在确定待测岩心中存在游离油时，执行前述步骤s201。若电子设备确定待测岩心中不存在游离油，可选的，电子设备可以不执行上述步骤s201，以节约计算资源。
[0129]
可选的，电子设备例如可以根据用户输入的指示信息，确定待测岩心中是否存在游离油。示例性的，若用户输入的指示信息包括用于指示该待测岩心中存在游离油的标识，则电子设备可以确定待测岩心中存在游离油。若用户输入的指示信息不包括用于指示该待测岩心中存在游离油的标识，则电子设备可以确定待测岩心中不存在游离油。
[0130]
具体实现时，用户例如可以将待测岩心的二维纵向弛豫
‑
横向弛豫时间谱、待测岩心对应的饱和油待测岩心的二维纵向弛豫
‑
横向弛豫时间谱，以及，标准二维核磁共振图版进行比对，以判断该待测岩心中可能存在的流体的类别。在确定该待测岩心中可能包括石油时，用户可以进一步对待测岩心的二维纵向弛豫
‑
横向弛豫时间谱、待测岩心对应的饱和油待测岩心的二维纵向弛豫
‑
横向弛豫时间谱，以及，标准二维核磁共振图版进行比对，以判断该待测岩心中是否存在游离油。若用户确定该待测岩心中存在游离油，则可以向电子设备输入用于表示该待测岩心中存在游离油的指示信息。然后，根据该指示信息，电子设备可以确定待测岩心中存在游离油。
[0131]
作为一种可能的实现方式，在获取待测岩心的游离油饱和度之后，电子设备还可以根据该待测岩心的游离油饱和度，确定页岩是否满足开采条件。然后，电子设备可以输出包括确定页岩是否满足开采条件的结果的指示信息，以使用户可以根据该指示信息及时知晓该页岩是否满足开采条件。
[0132]
可选的，在获取待测岩心的游离油饱和度之后，电子设备例如可以判断该游离油饱和度是否大于或等于预设饱和度阈值。若该游离油饱和度大于或等于预设饱和度阈值，说明该待测岩心所在的页岩的游离油饱和度大于或等于预设饱和度阈值，也就是说该页岩中存在的游离油的量较多，因此，电子设备可以确定该页岩满足开采条件。若该游离油饱和度小于预设饱和度阈值，说明该待测岩心所在的页岩的游离油饱和度小于预设饱和度阈值，也就是说该页岩中存在的游离油的量较少，因此，电子设备可以确定该页岩不满足开采
条件。
[0133]
若电子设备确定该页岩满足开采条件，示例性的，电子设备例如可以输出“页岩满足开采条件”的提示信息。若电子设备确定该页岩不满足开采条件，示例性的，电子设备例如可以输出“页岩不满足开采条件”的提示信息。
[0134]
图4为本技术提供的一种基于页岩油赋存状态表征的页岩游离油饱和度预测模型的构建装置的结构示意图。如图4所示，该装置300包括：第一获取模块301、第二获取模块302、第三获取模块303、处理模块304，以及，构建模块305。其中，
[0135]
第一获取模块301，用于获取n个来自所述页岩的样本岩心在m个横向弛豫时间区间对应的孔隙度。其中，所述样本岩心在m个横向弛豫时间区间对应的孔隙度的和为所述样本岩心的总孔隙度，所述n和所述m均为大于或等于2的整数、且所述n大于或等于所述m。
[0136]
第二获取模块302，用于获取n个所述样本岩心的页岩油赋存状态表征。
[0137]
第三获取模块303，用于根据所述n个所述样本岩心的页岩油赋存状态表征，获取所述n个所述样本岩心的游离油饱和度。
[0138]
处理模块304，用于根据n个所述样本岩心在m个横向弛豫时间区间对应的孔隙度，以及，n个所述样本岩心的游离油饱和度，确定每个所述横向弛豫时间区间对应的权重。
[0139]
构建模块305，用于根据每个所述横向弛豫时间区间对应的权重，构建所述页岩游离油饱和度预测模型。其中，所述页岩游离油饱和度预测模型用于根据待测岩心在m个横向弛豫时间区间对应的孔隙度，获取所述待测岩心的游离油饱和度；所述待测岩心所在的页岩的游离油饱和度等于所述待测岩心的游离油饱和度。
[0140]
可选的，处理模块304，具体用于构建每个所述样本岩心对应的求解权重方程；根据n个所述样本岩心对应的求解权重方程，求解各所述横向弛豫时间区间对应的权重的取值。其中，所述样本岩心对应的求解权重方程为：其中，所述ω
i
表示第i个横向弛豫时间区间对应的权重，所述bin
i
表示所述样本岩心在所述第i个横向弛豫时间区间对应的孔隙度，所述φ表示所述样本岩心的总孔隙度，所述p表示所述样本岩心的游离油饱和度。i为大于或等于1，且小于或等于m的整数。
[0141]
可选的，处理模块304，具体用于将所述n个样本岩心划分为至少一个样本岩心组；根据每组样本岩心对应的求解权重方程，获取每组样本岩心在各横向弛豫时间区间对应的初始权重；将所有组样本岩心在同一横向弛豫时间区间对应的初始权重的平均值，作为该横向弛豫时间区间对应的权重。其中，每个样本岩心组包括m个样本岩心，任意两个样本岩心组至少有一个样本岩心不同。
[0142]
可选的，第二获取模块302，具体用于针对n个所述样本岩心中的任一样本岩心，获取所述样本岩心的横向弛豫时间谱，以及，所述样本岩心对应的饱和油样本岩心的横向弛豫时间谱；根据所述样本岩心的横向弛豫时间谱，以及，所述样本岩心对应的饱和油样本岩心的横向弛豫时间谱，获取所述样本岩心的页岩油赋存状态表征。其中，所述饱和油样本岩心为对所述样本岩心进行抽真空，以及，饱和页岩原油得到的；所述样本岩心的横向弛豫时间谱用于表征所述样本岩心在各横向弛豫时间区间对应的束缚油饱和度，所述样本岩心对应的饱和油样本岩心的横向弛豫时间谱用于表征所述样本岩心在各横向弛豫时间区间对应的束缚油饱和度与游离油饱和度的和。
[0143]
本实施例提供的页岩游离油饱和度预测模型的构建装置300，用于执行前述页岩游离油饱和度预测模型的构建方法实施例，其实现原理与技术效果类似，对此不再赘述。
[0144]
图5为本技术提供的一种页岩游离油饱和度预测装置的结构示意图。如图5所示，该装置400包括：第一获取模块401、第二获取模块402，以及，输出模块403。其中，
[0145]
第一获取模块401，用于获取所述待测岩心的各横向弛豫时间区间对应的孔隙度。其中，所述待测岩心为从所述页岩开采得到的。
[0146]
第二获取模块402，用于根据所述待测岩心的各横向弛豫时间区间对应的孔隙度，利用如前述任一项实施例所述的游离油饱和度预测模型，获取所述待测岩心的游离油饱和度。其中，所述待测岩心所在的页岩的游离油饱和度等于所述待测岩心的游离油饱和度。
[0147]
输出模块403，用于输出所述页岩的游离油饱和度。
[0148]
可选的，该装置400还可以包括确定模块404，用于在所述获取所述待测岩心的各横向弛豫时间区间对应的孔隙度之前，确定所述待测岩心中存在游离油。
[0149]
可选的，输出模块403还用于在所述获取所述待测岩心的游离油饱和度之后，根据所述待测岩心的游离油饱和度，确定所述页岩是否满足开采条件；输出提示信息。其中，所述提示信息包括所述确定所述页岩是否满足开采条件的结果。
[0150]
本实施例提供的岩心游离油饱和度预测装置400，用于执行前述岩心游离油饱和度预测方法实施例，其实现原理与技术效果类似，对此不再赘述。
[0151]
图6为本技术提供的一种电子设备结构示意图。该电子设备应用于井下页岩核磁共振测井数据的处理与解释。如图6所示，该电子设备500可以包括：至少一个处理器501和存储器502。
[0152]
存储器502，用于存放程序。具体地，程序可以包括程序代码，程序代码包括计算机操作指令。
[0153]
存储器502可能包含高速ram存储器，也可能还包括非易失性存储器(non
‑
volatile memory)，例如至少一个磁盘存储器。
[0154]
处理器501用于执行存储器502存储的计算机执行指令，以实现前述方法实施例所描述的页岩游离油饱和度预测模型的构建或者页岩游离油饱和度预测方法。其中，处理器501可能是一个中央处理器(central processing unit，简称为cpu)，或者是特定集成电路(application specific integrated circuit，简称为asic)，或者是被配置成实施本技术实施例的一个或多个集成电路。
[0155]
可选的，该电子设备500还可以包括通信接口503。在具体实现上，如果通信接口503、存储器502和处理器501独立实现，则通信接口503、存储器502和处理器501可以通过总线相互连接并完成相互间的通信。总线可以是工业标准体系结构(industry standard architecture，简称为isa)总线、外部设备互连(peripheral component，简称为pci)总线或扩展工业标准体系结构(extended industry standard architecture，简称为eisa)总线等。总线可以分为地址总线、数据总线、控制总线等，但并不表示仅有一根总线或一种类型的总线。
[0156]
可选的，在具体实现上，如果通信接口503、存储器502和处理器501集成在一块芯片上实现，则通信接口503、存储器502和处理器501可以通过内部接口完成通信。
[0157]
本技术还提供了一种计算机可读存储介质，该计算机可读存储介质可以包括：u
盘、移动硬盘、只读存储器(rom，read
‑
only memory)、随机存取存储器(ram，random access memory)、磁盘或者光盘等各种可以存储程序代码的介质，具体的，该计算机可读存储介质中存储有程序指令，程序指令用于上述实施例中的方法。
[0158]
本技术还提供一种程序产品，该程序产品包括执行指令，该执行指令存储在可读存储介质中。电子设备的至少一个处理器可以从可读存储介质读取该执行指令，至少一个处理器执行该执行指令使得电子设备实施上述的各种实施方式提供的页岩游离油饱和度预测模型的构建或者页岩游离油饱和度预测方法。
[0159]
最后应说明的是：以上各实施例仅用以说明本技术的技术方案，而非对其限制；尽管参照前述各实施例对本技术进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本技术各实施例技术方案的范围。