页岩储层的孔隙度确定方法、装置及计算机存储介质与流程

1.本技术实施例涉及页岩气勘探开发技术领域，特别涉及一种页岩储层的孔隙度确定方法、装置及计算机存储介质。

背景技术：

2.随着页岩气在全球能源领域的迅速发展，对页岩气的勘探开发也来越来越受到重视。由于不同地区、不同井之间页岩储层存在差异，有高产井、低产井、无产井等，而页岩储层的产能与页岩储层的孔隙度相关，因此，需要确定页岩孔隙度。
3.但是，由于有的页岩储层电阻率高，有的电阻率低，即电阻率小于一定阈值时，该页岩储层为低阻页岩储层。目前，都是通过对正常电阻率页岩储层的孔隙度确定方式来确定低阻页岩储层的孔隙度，这导致确定的低阻页岩储层的孔隙度不准确，进而导致对低阻页岩储层的产能估计不准确。

技术实现要素：

4.本技术实施例提供了一种页岩储层的孔隙度确定方法、装置及计算机存储介质，可以用于解决相关技术中对低阻页岩储层的孔隙度确定不准确的问题。
5.所述技术方案如下：
6.一方面，提供了一种页岩储层的孔隙度确定方法，所述方法包括：
7.获取针对目标页岩储层的目标单井的测井数据；
8.根据所述测井数据，依次确定所述目标页岩储层内有机碳含量与储层密度之间的第一变化关系，以及所述目标页岩储层内有机碳含量与储层体积之间的第二变化关系；
9.根据所述第一变化关系、所述第二变化关系和所述测井数据，确定消除所述目标页岩储层内有机碳含量影响后的所述目标页岩储层的孔隙度。
10.在一些实施例中，所述测井数据包括自然伽马测井曲线、电阻率曲线、声波时差曲线和密度曲线；
11.所述根据所述测井数据，依次确定所述目标页岩储层内有机碳含量与储层密度之间的第一变化关系，以及所述目标页岩储层内有机碳含量与储层体积之间的第二变化关系，包括：
12.根据所述自然伽马测井曲线、所述电阻率曲线和所述声波时差曲线，确定所述目标页岩储层中的富有机质层段对应的测井曲线段和贫有机质层段对应的测井曲线段；
13.根据所述富有机质层段对应的测井曲线段和所述贫有机质层段对应的测井曲线段，确定所述目标页岩储层的有机碳含量确定模型；
14.将所述密度曲线的倒数与所述有机碳含量模型进行线性拟合，得到所述第一变化关系；
15.根据所述有机碳含量确定模型和所述第一变化关系，确定所述第二变化关系。
16.在一些实施例中，所述根据所述自然伽马测井曲线、所述电阻率曲线和所述声波
时差曲线，确定所述目标页岩储层中的富有机质层段对应的测井曲线段和贫有机质层段对应的测井曲线段，包括：
17.将所述电阻率曲线和所述声波时差曲线进行重叠，得到重叠曲线；
18.根据所述自然伽马曲线中的非有机质层段对应的曲线所在位置，将所述重叠曲线中非有机质层段对应的曲线剔除，得到有机质重叠曲线；
19.从所述有机质重叠曲线中确定富有机质层段对应的测井曲线段和贫有机质层段对应的测井曲线段，所述富有机质层段对应的测井曲线段为所述有机质重叠曲线中所述电阻率曲线和所述声波时差曲线走势不相同且不发生重合的测井曲线段，所述贫有机质层段对应的测井曲线段为所述有机质重叠曲线中所述电阻率曲线和所述声波时差曲线走势相同且发生重合的测井曲线段。
20.在一些实施例中，所述根据所述富有机质层段对应的测井曲线段和所述贫有机质层段对应的测井曲线段，确定所述目标页岩储层的有机碳含量确定模型，包括：
21.获取所述目标页岩储层的地层电阻率和不同地层对应的有机碳含量；
22.在所述贫有机质层段对应的测井曲线段中确定声波时差基线和电阻率基线，所述声波时差基线为所述贫有机质层段对应的测井曲线段重叠部分中声波对应的基线，所述电阻率基线为所述贫有机质层段对应的测井曲线段重叠部分中电阻率对应的基线；
23.分别根据所述贫有机质层段对应的声波时差基线和电阻率基线，在所述富有机质层段对应的测井曲线段中确定参考声波时差基线和参考电阻率基线，并确定所述富有机质层段对应的测井曲线段中的参考声波时差曲线；
24.根据所述参考声波时差基线、所述参考电阻率基线、所述参考声波时差曲线和所述地层电阻率确定不同地层对应的有机质参考丰度；
25.将所述不同地层对应的有机质参考丰度和所述不同地层对应的有机碳含量进行线性拟合，得到所述有机碳含量确定模型。
26.在一些实施例中，所述根据所述有机碳含量确定模型和所述第一变化关系，确定所述第二变化关系，包括：
27.从所述第一变化关系中确定所述有机碳含量为0和100％时分别对应的密度；
28.将所述有机碳含量为0时对应的密度确定为所述目标页岩储层的岩石骨架密度，并将所述有机碳含量为100％时对应的密度确定为所述有机碳含量的密度；
29.根据所述有机碳含量密度、所述有机碳含量确定模型和不同地层对应的参考密度，确定所述第二变化关系，所述不同地层对应的参考密度通过所述密度曲线获取得到。
30.在一些实施例中，所述根据所述第一变化关系和所述第二变化关系，确定消除所述目标页岩储层内有机碳含量影响后的所述目标页岩储层的孔隙度，包括：
31.根据所述第二变化关系、所述第一变化关系中描述的有机碳含量的密度和不同地层对应的参考密度，确定消除所述有机碳含量的影响后所述目标页岩储层的不同地层的地层密度；
32.根据所述不同地层的地层密度、所述岩石骨架密度和所述目标页岩储层内流体的流体密度，确定所述不同地层的孔隙度。
33.另一方面，提供了一种页岩储层的孔隙度确定装置，所述装置包括：
34.获取模块，用于获取针对目标页岩储层的目标单井的测井数据；
35.第一确定模块，用于根据所述测井数据，依次确定所述目标页岩储层内有机碳含量与储层密度之间的第一变化关系，以及所述目标页岩储层内有机碳含量与储层体积之间的第二变化关系；
36.第二确定模块，用于根据所述第一变化关系、所述第二变化关系和所述测井数据，确定消除所述目标页岩储层内有机碳含量影响后的所述目标页岩储层的孔隙度。
37.在一些实施例中，所述测井数据包括自然伽马测井曲线、电阻率曲线、声波时差曲线和密度曲线；
38.所述第一确定模块包括：
39.第一确定子模块，用于根据所述自然伽马测井曲线、所述电阻率曲线和所述声波时差曲线，确定所述目标页岩储层中的富有机质层段对应的测井曲线段和贫有机质层段对应的测井曲线段；
40.第二确定子模块，用于根据所述富有机质层段对应的测井曲线段和所述贫有机质层段对应的测井曲线段，确定所述目标页岩储层的有机碳含量确定模型；
41.拟合子模块，用于将所述密度曲线的倒数与所述有机碳含量模型进行线性拟合，得到所述第一变化关系；
42.第三确定子模块，用于根据所述有机碳含量确定模型和所述第一变化关系，确定所述第二变化关系。
43.在一些实施例中，所述第一确定子模块用于：
44.将所述电阻率曲线和所述声波时差曲线进行重叠，得到重叠曲线；
45.根据所述自然伽马曲线中的非有机质层段对应的曲线所在位置，将所述重叠曲线中非有机质层段对应的曲线剔除，得到有机质重叠曲线；
46.从所述有机质重叠曲线中确定富有机质层段对应的测井曲线段和贫有机质层段对应的测井曲线段，所述富有机质层段对应的测井曲线段为所述有机质重叠曲线中所述电阻率曲线和所述声波时差曲线走势不相同且不发生重合的测井曲线段，所述贫有机质层段对应的测井曲线段为所述有机质重叠曲线中所述电阻率曲线和所述声波时差曲线走势相同且发生重合的测井曲线段。
47.在一些实施例中，所述第二确定子模块用于：
48.获取所述目标页岩储层的地层电阻率和不同地层对应的有机碳含量；
49.在所述贫有机质层段对应的测井曲线段中确定声波时差基线和电阻率基线，所述声波时差基线为所述贫有机质层段对应的测井曲线段重叠部分中声波对应的基线，所述电阻率基线为所述贫有机质层段对应的测井曲线段重叠部分中电阻率对应的基线；
50.分别根据所述贫有机质层段对应的声波时差基线和电阻率基线，在所述富有机质层段对应的测井曲线段中确定参考声波时差基线和参考电阻率基线，并确定所述富有机质层段对应的测井曲线段中的参考声波时差曲线；
51.根据所述参考声波时差基线、所述参考电阻率基线、所述参考声波时差曲线和所述地层电阻率确定不同地层对应的有机质参考丰度；
52.将所述不同地层对应的有机质参考丰度和所述不同地层对应的有机碳含量进行线性拟合，得到所述有机碳含量确定模型。
53.在一些实施例中，所述第三确定子模块用于：
54.从所述第一变化关系中确定所述有机碳含量为0和100％时分别对应的密度；
55.将所述有机碳含量为0时对应的密度确定为所述目标页岩储层的岩石骨架密度，并将所述有机碳含量为100％时对应的密度确定为所述有机碳含量的密度；
56.根据所述有机碳含量密度、所述有机碳含量确定模型和不同地层对应的参考密度，确定所述第二变化关系，所述不同地层对应的参考密度通过所述密度曲线获取得到。
57.在一些实施例中，所述第二确定模块包括：
58.第四确定子模块，用于根据所述第二变化关系、所述第一变化关系中描述的有机碳含量的密度和不同地层对应的参考密度，确定消除所述有机碳含量的影响后所述目标页岩储层的不同地层的地层密度；
59.第五确定子模块，用于根据所述不同地层的地层密度、所述岩石骨架密度和所述目标页岩储层内流体的流体密度，确定所述不同地层的孔隙度。
60.另一方面，提供了一种计算机可读存储介质，所述计算机可读存储介质上存储有指令，所述指令被处理器执行时实现上述页岩储层的孔隙度确定方法中的任一步骤。
61.本技术实施例提供的技术方案带来的有益效果至少包括：
62.在本技术实施例中，能够通过目标页岩储层内有机碳含量与储层密度之间的关系，以及有机碳含量与储层体积之间的关系，确定消除目标页岩储层内有机碳含量影响后的目标页岩储层的孔隙度，从而提高了确定页岩储层的孔隙度的准确性，进而提高了对低阻页岩储层的产能估计的准确性。
附图说明
63.为了更清楚地说明本技术实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
64.图1是本技术实施例提供的一种页岩储层的孔隙度确定方法流程图；
65.图2是本技术实施例提供的一种页岩储层的孔隙度确定方法流程图；
66.图3是本技术实施例提供的一种目标页岩储层内有机碳含量与储层密度之间的关系示意图；
67.图4是本技术实施例提供的一种目标页岩储层内有机碳含量与储层体积之间的关系示意图；
68.图5是本技术实施例提供的一种页岩储层的孔隙度确定装置的结构示意图；
69.图6是本技术实施例提供的一种第一确定模块的结构示意图；
70.图7是本技术实施例提供的一种第二确定模块的结构示意图；
71.图8是本技术实施例提供的一种终端的结构示意图。
具体实施方式
72.为使本技术实施例的目的、技术方案和优点更加清楚，下面将结合附图对本技术实施方式作进一步地详细描述。
73.在对本技术实施例提供的一种页岩储层的孔隙度确定方法进行详细的解释说明
之前，先对本技术实施例提供的一种应用场景进行解释说明。
74.由于页岩的孔隙度是评价页岩储层的重要物性特征，在页岩油、页岩气的勘探和开发中，弄清页岩储层中的孔隙度的大小，可以为评价页岩油和页岩气的赋存方式、流动状态和可采性等提供依据和基础。但是，由于目前都是通过对正常电阻率页岩储层的孔隙度确定方式来确定低阻页岩储层的孔隙度，这导致确定的低阻页岩储层的孔隙度不准确，进而导致对低阻页岩储层的产能估计不准确。
75.基于这样的应用场景，本技术实施例提供了一种提高确定低阻页岩储层的孔隙度准确性的页岩储层的孔隙度确定方法。
76.图1是本技术实施例提供的一种页岩储层的孔隙度确定方法流程图，该页岩储层的孔隙度确定方法可以包括如下几个步骤：
77.步骤101：获取针对目标页岩储层的目标单井的测井数据。
78.步骤102：根据该测井数据，依次确定该目标页岩储层内有机碳含量与储层密度之间的第一变化关系，以及该目标页岩储层内有机碳含量与储层体积之间的第二变化关系。
79.步骤103：根据该第一变化关系、该第二变化关系和该测井数据，确定消除该目标页岩储层内有机碳含量影响后的该目标页岩储层的孔隙度。
80.在本技术实施例中，能够通过目标页岩储层内有机碳含量与储层密度之间的关系，以及有机碳含量与储层体积之间的关系，确定消除目标页岩储层内有机碳含量影响后的目标页岩储层的孔隙度，从而提高了确定页岩储层的孔隙度的准确性，进而提高了对低阻页岩储层的产能估计的准确性。
81.在一些实施例中，该测井数据包括自然伽马测井曲线、电阻率曲线、声波时差曲线和密度曲线；
82.根据该测井数据，依次确定该目标页岩储层内有机碳含量与储层密度之间的第一变化关系，以及该目标页岩储层内有机碳含量与储层体积之间的第二变化关系，包括：
83.根据该自然伽马测井曲线、该电阻率曲线和该声波时差曲线，确定该目标页岩储层中的富有机质层段对应的测井曲线段和贫有机质层段对应的测井曲线段；
84.根据该富有机质层段对应的测井曲线段和贫有机质层段对应的测井曲线段，确定该目标页岩储层的有机碳含量确定模型；
85.将该密度曲线的倒数与该有机碳含量模型进行线性拟合，得到该第一变化关系；
86.根据该有机碳含量确定模型和该第一变化关系，确定该第二变化关系。
87.在一些实施例中，根据该自然伽马测井曲线、该电阻率曲线和该声波时差曲线，确定该目标页岩储层中的富有机质层段对应的测井曲线段和贫有机质层段对应的测井曲线段，包括：
88.将该电阻率曲线和该声波时差曲线进行重叠，得到重叠曲线；
89.根据该自然伽马曲线中的非有机质层段对应的曲线所在位置，将该重叠曲线中非有机质层段对应的曲线剔除，得到有机质重叠曲线；
90.从该有机质重叠曲线中确定富有机质层段对应的测井曲线段和贫有机质层段对应的测井曲线段，该富有机质层段对应的测井曲线段为该有机质重叠曲线中该电阻率曲线和该声波时差曲线走势不相同且不发生重合的测井曲线段，该贫有机质层段对应的测井曲线段为该有机质重叠曲线中该电阻率曲线和该声波时差曲线走势相同且发生重合的测井
曲线段。
91.在一些实施例中，根据该富有机质层段对应的测井曲线段确定该目标页岩储层的有机碳含量确定模型和该贫有机质层段对应的测井曲线段，包括：
92.获取该目标页岩储层的地层电阻率和不同地层对应的有机碳含量；
93.在该贫有机质层段对应的测井曲线段中确定声波时差基线和电阻率基线，该声波时差基线为该贫有机质层段对应的测井曲线段重叠部分中声波对应的基线，该电阻率基线为该贫有机质层段对应的测井曲线段重叠部分中电阻率对应的基线；
94.分别根据该贫有机质层段对应的声波时差基线和电阻率基线，在该富有机质层段对应的测井曲线段中确定参考声波时差基线和参考电阻率基线，并确定该富有机质层段对应的测井曲线段中的参考声波时差曲线；
95.根据该参考声波时差基线、该参考电阻率基线、该参考声波时差曲线和该地层电阻率确定不同地层对应的有机质参考丰度；
96.将该不同地层对应的有机质参考丰度和该不同地层对应的有机碳含量进行线性拟合，得到该有机碳含量确定模型。
97.在一些实施例中，根据该有机碳含量确定模型和该第一变化关系，确定该第二变化关系，包括：
98.从该第一变化关系中确定该有机碳含量为0和100％时分别对应的密度；
99.将该有机碳含量为0时对应的密度确定为该目标页岩储层的岩石骨架密度，并将该有机碳含量为100％时对应的密度确定为该有机碳含量的密度；
100.根据该有机碳含量密度、该有机碳含量确定模型和不同地层对应的参考密度，确定该第二变化关系，该不同地层对应的参考密度通过该密度曲线获取得到。
101.在一些实施例中，根据该第一变化关系和该第二变化关系，确定消除该目标页岩储层内有机碳含量影响后的该目标页岩储层的孔隙度，包括：
102.根据该第二变化关系、该第一变化关系中描述的有机碳含量的密度和不同地层对应的参考密度，确定消除该有机碳含量的影响后该目标页岩储层的不同地层的地层密度；
103.根据该不同地层的地层密度、该岩石骨架密度和该目标页岩储层内流体的流体密度，确定该不同地层的孔隙度。
104.上述所有可选技术方案，均可按照任意结合形成本技术的可选实施例，本技术实施例对此不再一一赘述。
105.图2是本技术实施例提供的一种页岩储层的孔隙度确定方法流程图，本实施例以该页岩储层的孔隙度确定方法应用于终端中进行举例说明，该页岩储层的孔隙度确定方法可以包括如下几个步骤：
106.步骤201：终端获取针对目标页岩储层的目标单井的测井数据。
107.由于单井的测井数据能与页岩储层的孔隙度存在关系，因此，为了确定目标页岩储层的孔隙度，终端能够获取针对目标页岩储层的目标单井的测井数据。该目标页岩储层能够为任一低电阻率页岩储层，该目标单井为目标页岩储层范围内的任一单井。
108.作为一种示例，终端能够在接收到获取指令时，从本地存储的文件中获取目标单井的测井数据。也即是，终端能够事先获取并存储目标单井的测井数据，并在接收到获取指令时，从本地存储的文件中获取目标单井的测井数据。或者，终端还能够在接收到获取指令
时，向其他设备发送获取请求，该获取请求中携带目标单井的单井标识；其他设备在接收到终端发生的获取请求时，根据目标单井的单井标识向终端返回目标单井的测井数据，从而使终端获取到目标单井的测井数据。
109.需要说明的是，获取指令为用户通过指定操作作用在终端上时触发，该指定操作包括语音操作、点击操作、滑动操作等等。
110.在一些实施例中，测井数据至少包括自然伽马测井曲线(gr)、电阻率曲线(rt)、声波时差曲线和密度曲线。
111.步骤202：终端根据测井数据，依次确定目标页岩储层内有机碳含量与储层密度之间的第一变化关系，以及该目标页岩储层内有机碳含量与储层体积之间的第二变化关系。
112.由于目标页岩储层中有机碳含量的多少影响到储层密度和体积，进而影响到目标页岩储层的孔隙度，因此，终端需要根据测井数据，依次确定目标页岩储层内有机碳含量与储层密度之间的第一变化关系，以及该目标页岩储层内有机碳含量与储层体积之间的第二变化关系。
113.作为一种示例，终端根据测井数据，依次确定目标页岩储层内有机碳含量与储层密度之间的第一变化关系，以及目标页岩储层内有机碳含量与储层体积之间的第二变化关系的操作至少包括：根据自然伽马测井曲线、电阻率曲线和声波时差曲线，确定目标页岩储层中的富有机质层段对应的测井曲线段和贫有机质层段对应的测井曲线段；根据富有机质层段对应的测井曲线段和贫有机质层段对应的测井曲线段，确定目标页岩储层的有机碳含量确定模型；将密度曲线的倒数与有机碳含量模型进行线性拟合，得到第一变化关系；根据有机碳含量确定模型和第一变化关系，确定第二变化关系。
114.由于影响目标页岩储层的孔隙度的层段主要为有机质层段，因此，终端需要根据自然伽马测井曲线、电阻率曲线和声波时差曲线，确定目标页岩储层中的富有机质层段对应的测井曲线段和贫有机质层段对应的测井曲线段。
115.作为一种示例，终端根据自然伽马测井曲线、电阻率曲线和声波时差曲线，确定目标页岩储层中的富有机质层段对应的测井曲线段和贫有机质层段对应的测井曲线段的操作包括：将电阻率曲线和声波时差曲线进行重叠，得到重叠曲线；根据自然伽马曲线中的非有机质层段对应的曲线所在位置，将重叠曲线中非有机质层段对应的曲线剔除，得到有机质重叠曲线；从有机质重叠曲线中确定富有机质层段对应的测井曲线段和贫有机质层段对应的测井曲线段，该富有机质层段对应的测井曲线段为有机质重叠曲线中电阻率曲线和声波时差曲线走势不相同且不发生重合的测井曲线段，贫有机质层段对应的测井曲线段为有机质重叠曲线中电阻率曲线和声波时差曲线走势相同且发生重合的测井曲线段。
116.在一些实施例中，终端将电阻率曲线和声波时差曲线进行重叠是指将带有刻度的声波时差曲线在深测向上叠加在电阻率曲线上，从而得到重叠曲线。
117.需要说明的是，终端将电阻率曲线和声波时差曲线进行重叠的操作能够参考相关技术，本技术实施例对此不再进行一一赘述。
118.由于非机质层段对目标页岩储层的孔隙度没有影响或者影响几乎可忽略，且由于自然伽马曲线能够反应不同的储层结构，因此，终端能够通过根据自然伽马曲线中的非有机质层段所在位置，将重叠曲线中非有机质层段对应的曲线剔除，得到有机质重叠曲线。
119.在一些实施例中，终端能够从自然伽马曲线中识别油层、蒸发岩等非有机质层段
对应的曲线位置，根据自然伽马曲线中非有机质层段对应的曲线位置，将重叠曲线中对应的非有机质层段对应的曲线进行剔除，得到有机质重叠曲线。
120.由于有机质重叠曲线中可能包括贫有机质层段对应的测井曲线段和富有机质层段对应的测井曲线段，因此，终端还能够从有机质重叠曲线中确定富有机质层段对应的测井曲线段和贫有机质层段对应的测井曲线段。
121.在一些实施例中，终端能够将有机质重叠曲线中电阻率曲线和声波时差曲线走势相同且重合在的测井曲线段确定为贫有机质层段(通常为非烃源岩)对应的测井曲线段。将有机质重叠曲线中电阻率曲线和声波时差曲线走势不相同且不发生重合的测井曲线段确定为富有机质层段对应的测井曲线段。
122.在一些实施例中，终端据富有机质层段对应的测井曲线段和贫有机质层段对应的测井曲线段，确定目标页岩储层的有机碳含量确定模型的操作包括：获取目标页岩储层的地层电阻率和不同地层对应的有机碳含量；在贫有机质层段对应的测井曲线段中确定声波时差基线和电阻率基线，声波时差基线为贫有机质层段对应的测井曲线段重叠部分中声波对应的基线，电阻率基线为贫有机质层段对应的测井曲线段重叠部分中电阻率对应的基线；分别根据贫有机质层段对应的声波时差基线和电阻率基线，在富有机质层段对应的测井曲线中确定参考声波时差基线和参考电阻率基线，并确定富有机质层段对应的测井曲线段中的参考声波时差曲线；根据参考声波时差基线、参考电阻率基线、参考声波时差曲线和地层电阻率确定不同地层对应的有机质参考丰度；将不同地层对应的有机质参考丰度和不同地层对应的有机碳含量进行线性拟合，得到有机碳含量确定模型。
123.由于在目标页岩储层中，不同深度地层对应的有机碳含量不同，因此，终端需要获取目标页岩储层的地层电阻率和不同地层对应的有机碳含量。
124.又由于贫有机质层段对孔隙度影响几乎可忽略不计，且贫有机质层段中的电阻率基线和声波时差基线与富有机质层段中的电阻率基线和声波时差基线相关，因此，终端能够先确定贫有机质层段中的电阻率基线和声波时差基线，然后分别根据贫有机质层段中的电阻率基线和声波时差基线，在富有机质层段对应的测井曲线段中确定对应的参考声波时差基线和参考电阻率基线。
125.在一些实施例中，终端能够与将贫有机质层段中的电阻率基线和声波时差基线分别确定为富有机质层段对应参考电阻率基线和参考声波时差基线。
126.需要说明的是，参考声波时差曲线为富有机质层段对应的测井曲线段中任一段声波时差曲线。
127.在一些实施例中，终端根据参考声波时差基线、参考电阻率基线、参考声波时差曲线和地层电阻率，能够通过下述第一公式确定不同地层对应的有机质参考丰度。
[0128][0129]
需要说明的是，在上述第一公式(1)中，index为不同地层对应的有机质参考丰度，r
t
为参考电阻率基线，δt为参考声波时差曲线，r
bl
为地层电阻率，δt
bl
为参考声波时差基线，d为常数。
[0130]
在一些实施例中，终端将不同地层对应的有机质参考丰度和不同地层对应的有机碳含量进行线性拟合，得到的有机碳含量确定模型能够通过下述图3所示的图像进行表示，
也能够通过下述第二公式所示的公式进行表示。
[0131]
toc＝a+b*index
ꢀꢀꢀ
(2)
[0132]
需要说明的是，在上述第二公式(2)中，a和b均为常数，toc为有机碳含量，index为不同地层对应的有机质参考丰度。
[0133]
为了确定目标页岩储层内有机碳含量与储层密度之间的第一变化关系，终端能够将目标页岩储层的密度曲线与有机碳含量进行拟合。即终端将密度曲线的倒数与有机碳含量模型进行线性拟合，得到第一变化关系。
[0134]
需要说明的是，有机碳含量模型与密度曲线的倒数拟合后，能够得到如图4所示的第一变化关系。
[0135]
由于随着有机碳含量和密度发生变化，储层体积也将会发生变化，因此，终端能够据有机碳含量确定模型和第一变化关系，确定第二变化关系。
[0136]
在一些实施例中，终端根据有机碳含量确定模型和第一变化关系，确定第二变化关系的操作包括：从第一变化关系中确定有机碳含量为0和100％时分别对应的密度；将有机碳含量为0时对应的密度确定为目标页岩储层的岩石骨架密度，并将有机碳含量为100％时对应的密度确定为有机碳含量的密度；根据有机碳含量密度、有机碳含量确定模型和不同地层对应的参考密度，确定第二变化关系，不同地层对应的参考密度通过密度曲线获取得到。
[0137]
由于第一变化关系为目标页岩储层内有机碳含量与储层密度之间的关系，当有机碳含量为0时，说明已不存在有机碳含量，那么目标页岩储层此时的储层密度即为目标页岩储层的岩石骨架密度，而当有机碳含量为100％时，说明不存在其他物质，因此，目标页岩储层此时的储层密度即为有机碳含量的密度。
[0138]
作为一种示例，终端根据有机碳含量密度、有机碳含量确定模型和不同地层对应的参考密度，通过下述第三公式确定第二变化关系。
[0139]vk
＝(toc*ρb*υ)/ρkꢀꢀꢀ
(3)
[0140]
需要说明的是，在上述第三公式(3)中，vk为不同地层(或称为不同深度地层)对应的体积，toc为通过上述第二公式(2)确定的有机碳含量，ρb为不同地层对应的参考密度，ρk为有机碳含量的密度(即当有机碳含量为100％时对应的密度)，υ为常数，通常取值为1.19。
[0141]
步骤203：终端根据第一变化关系和第二变化关系，确定消除目标页岩储层内有机碳含量影响后的目标页岩储层的孔隙度。
[0142]
作为一种示例，终端根据第一变化关系和第二变化关系，确定消除目标页岩储层内有机碳含量影响后的目标页岩储层的孔隙度的操作包括：根据第二变化关系、第一变化关系中描述的有机碳含量的密度和不同地层对应的参考密度，确定消除有机碳含量的影响后目标页岩储层的不同地层的地层密度；根据不同地层的地层密度、岩石骨架密度和目标页岩储层内流体的流体密度，确定不同地层的孔隙度。
[0143]
在一些实施例中，终端根据第二变化关系、第一变化关系中描述的有机碳含量的密度和不同地层对应的参考密度，通过下述第四公式确定消除有机碳含量的影响后目标页岩储层的不同地层的地层密度。
[0144]
ρ
bc
＝(ρ
b-vk*ρk)/(1-vk)
ꢀꢀꢀ
(4)
[0145]
需要说明的是，在上述第四公式(4)中，ρ
bc
为消除有机碳含量的影响后目标页岩储
层的不同地层的地层密度,vk为通过上述第三公式(第二关系)确定的不同地层对应的体积，ρb为不同地层对应的参考密度，ρk为有机碳含量的密度。
[0146]
在一些实施例中，终端根据不同地层的地层密度、岩石骨架密度和目标页岩储层内流体的流体密度，通过下述第五公式确定不同地层的孔隙度。
[0147]
φ＝(ρ
ma-ρ
bc
)/(ρ
ma-ρf)
ꢀꢀꢀ
(5)
[0148]
需要说明的是，在上述第五公式(5)中，φ为目标岩储层的不同地层的孔隙度，ρ
bc
为消除有机碳含量的影响后目标页岩储层的不同地层的地层密度，ρ
ma
为目标页岩储层的岩石骨架密度，ρf为流体密度。
[0149]
步骤204：终端通过第一提示信息提示目标页岩储层的孔隙度。
[0150]
为了使工作人员清楚的了解到目标页岩储层的孔隙度，终端能够在确定目标页岩储层的孔隙度后，通过第一提示信息提示目标页岩储层的孔隙度。
[0151]
在一些实施例中，为了使工作人员了解到测井得到的孔隙度与实际得到的孔隙度之间的差异，终端还能够将测井得到的孔隙度与通过上述方式得到的孔隙度进行对比，并将对比结果通过第二提示信息进行提示。
[0152]
需要说明的是，第一提示信息和第二提示信息均能够为文字、图像、视频等至少一种形式的提示信息。
[0153]
在本技术实施例中，终端能够通过目标页岩储层内有机碳含量与储层密度之间的关系，以及有机碳含量与储层体积之间的关系，确定消除目标页岩储层内有机碳含量影响后的目标页岩储层的孔隙度，从而提高了确定页岩储层的孔隙度的准确性，进而提高了对低阻页岩储层的产能估计的准确性。
[0154]
图5是本技术实施例提供的一种页岩储层的孔隙度确定装置的结构示意图，该页岩储层的孔隙度确定装置可以由软件、硬件或者两者的结合实现。该页岩储层的孔隙度确定装置可以包括：获取模块501、第一确定模块502和第二确定模块503。
[0155]
获取模块501，用于获取针对目标页岩储层的目标单井的测井数据；
[0156]
第一确定模块502，用于根据所述测井数据，依次确定所述目标页岩储层内有机碳含量与储层密度之间的第一变化关系，以及所述目标页岩储层内有机碳含量与储层体积之间的第二变化关系；
[0157]
第二确定模块503，用于根据所述第一变化关系、所述第二变化关系和所述测井数据，确定消除所述目标页岩储层内有机碳含量影响后的所述目标页岩储层的孔隙度。
[0158]
在一些实施例中，所述测井数据包括自然伽马测井曲线、电阻率曲线、声波时差曲线和密度曲线；
[0159]
参见图6，所述第一确定模块502包括：
[0160]
第一确定子模块5021，用于根据所述自然伽马测井曲线、所述电阻率曲线和所述声波时差曲线，确定所述目标页岩储层中的富有机质层段对应的测井曲线段和贫有机质层段对应的测井曲线段；
[0161]
第二确定子模块5022，用于根据所述富有机质层段对应的测井曲线段和所述贫有机质层段对应的测井曲线段，确定所述目标页岩储层的有机碳含量确定模型；
[0162]
拟合子模块5023，用于将所述密度曲线的倒数与所述有机碳含量模型进行线性拟合，得到所述第一变化关系；
[0163]
第三确定子模块5024，用于根据所述有机碳含量确定模型和所述第一变化关系，确定所述第二变化关系。
[0164]
在一些实施例中，所述第一确定子模块5021用于：
[0165]
将所述电阻率曲线和所述声波时差曲线进行重叠，得到重叠曲线；
[0166]
根据所述自然伽马曲线中的非有机质层段对应的曲线所在位置，将所述重叠曲线中非有机质层段对应的曲线剔除，得到有机质重叠曲线；
[0167]
从所述有机质重叠曲线中确定富有机质层段对应的测井曲线段和贫有机质层段对应的测井曲线段，所述富有机质层段对应的测井曲线段为所述有机质重叠曲线中所述电阻率曲线和所述声波时差曲线走势不相同且不发生重合的测井曲线段，所述贫有机质层段对应的测井曲线段为所述有机质重叠曲线中所述电阻率曲线和所述声波时差曲线走势相同且发生重合的测井曲线段。
[0168]
在一些实施例中，所述第二确定子模块5022用于：
[0169]
获取所述目标页岩储层的地层电阻率和不同地层对应的有机碳含量；
[0170]
在所述贫有机质层段对应的测井曲线段中确定声波时差基线和电阻率基线，所述声波时差基线为所述贫有机质层段对应的测井曲线段重叠部分中声波对应的基线，所述电阻率基线为所述贫有机质层段对应的测井曲线段重叠部分中电阻率对应的基线；
[0171]
分别根据所述贫有机质层段对应的声波时差基线和电阻率基线，在所述富有机质层段对应的测井曲线段中确定参考声波时差基线和参考电阻率基线，并确定所述富有机质层段对应的测井曲线段中的参考声波时差曲线；
[0172]
根据所述参考声波时差基线、所述参考电阻率基线、所述参考声波时差曲线和所述地层电阻率确定不同地层对应的有机质参考丰度；
[0173]
将所述不同地层对应的有机质参考丰度和所述不同地层对应的有机碳含量进行线性拟合，得到所述有机碳含量确定模型。
[0174]
在一些实施例中，所述第三确定子模块5024用于：
[0175]
从所述第一变化关系中确定所述有机碳含量为0和100％时分别对应的密度；
[0176]
将所述有机碳含量为0时对应的密度确定为所述目标页岩储层的岩石骨架密度，并将所述有机碳含量为100％时对应的密度确定为所述有机碳含量的密度；
[0177]
根据所述有机碳含量密度、所述有机碳含量确定模型和不同地层对应的参考密度，确定所述第二变化关系，所述不同地层对应的参考密度通过所述密度曲线获取得到。
[0178]
在一些实施例中，参见图7，所述第二确定模块503包括：
[0179]
第四确定子模块5031，用于根据所述第二变化关系、所述第一变化关系中描述的有机碳含量的密度和不同地层对应的参考密度，确定消除所述有机碳含量的影响后所述目标页岩储层的不同地层的地层密度；
[0180]
第五确定子模块5032，用于根据所述不同地层的地层密度、所述岩石骨架密度和所述目标页岩储层内流体的流体密度，确定所述不同地层的孔隙度。
[0181]
在本技术实施例中，终端能够通过目标页岩储层内有机碳含量与储层密度之间的关系，以及有机碳含量与储层体积之间的关系，确定消除目标页岩储层内有机碳含量影响后的目标页岩储层的孔隙度，从而提高了确定页岩储层的孔隙度的准确性，进而提高了对低阻页岩储层的产能估计的准确性。
[0182]
需要说明的是：上述实施例提供的页岩储层的孔隙度确定装置在确定页岩储层的孔隙度时，仅以上述各功能模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将装置的内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。另外，上述实施例提供的页岩储层的孔隙度确定装置与页岩储层的孔隙度确定方法实施例属于同一构思，其具体实现过程详见方法实施例，这里不再赘述。
[0183]
图8示出了本技术一个示例性实施例提供的终端800的结构框图。该终端800可以是：智能手机、平板电脑、mp3播放器(moving picture experts group audio layer iii，动态影像专家压缩标准音频层面3)、mp4(moving picture experts group audio layer iv，动态影像专家压缩标准音频层面4)播放器、笔记本电脑或台式电脑。终端800还可能被称为用户设备、便携式终端、膝上型终端、台式终端等其他名称。
[0184]
通常，终端800包括有：处理器801和存储器802。
[0185]
处理器801可以包括一个或多个处理核心，比如4核心处理器、8核心处理器等。处理器801可以采用dsp(digital signal processing，数字信号处理)、fpga(field－programmable gate array，现场可编程门阵列)、pla(programmable logic array，可编程逻辑阵列)中的至少一种硬件形式来实现。处理器801也可以包括主处理器和协处理器，主处理器是用于对在唤醒状态下的数据进行处理的处理器，也称cpu(central processing unit，中央处理器)；协处理器是用于对在待机状态下的数据进行处理的低功耗处理器。在一些实施例中，处理器801可以集成有gpu(graphics processing unit，图像处理器)，gpu用于负责显示屏所需要显示的内容的渲染和绘制。一些实施例中，处理器801还可以包括ai(artificial intelligence，人工智能)处理器，该ai处理器用于处理有关机器学习的计算操作。
[0186]
存储器802可以包括一个或多个计算机可读存储介质，该计算机可读存储介质可以是非暂态的。存储器802还可包括高速随机存取存储器，以及非易失性存储器，比如一个或多个磁盘存储设备、闪存存储设备。在一些实施例中，存储器802中的非暂态的计算机可读存储介质用于存储至少一个指令，该至少一个指令用于被处理器801所执行以实现本技术中方法实施例提供的页岩储层的孔隙度确定方法。
[0187]
在一些实施例中，终端800还可选包括有：外围设备接口803和至少一个外围设备。处理器801、存储器802和外围设备接口803之间可以通过总线或信号线相连。各个外围设备可以通过总线、信号线或电路板与外围设备接口803相连。具体地，外围设备包括：射频电路804、显示屏805、摄像头组件806、音频电路807、定位组件808和电源809中的至少一种。
[0188]
外围设备接口803可被用于将i/o(input/output，输入/输出)相关的至少一个外围设备连接到处理器801和存储器802。在一些实施例中，处理器801、存储器802和外围设备接口803被集成在同一芯片或电路板上；在一些其他实施例中，处理器801、存储器802和外围设备接口803中的任意一个或两个可以在单独的芯片或电路板上实现，本实施例对此不加以限定。
[0189]
射频电路804用于接收和发射rf(radio frequency，射频)信号，也称电磁信号。射频电路804通过电磁信号与通信网络以及其他通信设备进行通信。射频电路804将电信号转换为电磁信号进行发送，或者，将接收到的电磁信号转换为电信号。可选地，射频电路804包
括：天线系统、rf收发器、一个或多个放大器、调谐器、振荡器、数字信号处理器、编解码芯片组、用户身份模块卡等等。射频电路804可以通过至少一种无线通信协议来与其它终端进行通信。该无线通信协议包括但不限于：城域网、各代移动通信网络(2g、3g、4g及5g)、无线局域网和/或wifi(wireless fidelity，无线保真)网络。在一些实施例中，射频电路804还可以包括nfc(near field communication，近距离无线通信)有关的电路，本技术对此不加以限定。
[0190]
显示屏805用于显示ui(user interface，用户界面)。该ui可以包括图形、文本、图标、视频及其它们的任意组合。当显示屏805是触摸显示屏时，显示屏805还具有采集在显示屏805的表面或表面上方的触摸信号的能力。该触摸信号可以作为控制信号输入至处理器801进行处理。此时，显示屏805还可以用于提供虚拟按钮和/或虚拟键盘，也称软按钮和/或软键盘。在一些实施例中，显示屏805可以为一个，设置终端800的前面板；在另一些实施例中，显示屏805可以为至少两个，分别设置在终端800的不同表面或呈折叠设计；在另一些实施例中，显示屏805可以是柔性显示屏，设置在终端800的弯曲表面上或折叠面上。甚至，显示屏805还可以设置成非矩形的不规则图形，也即异形屏。显示屏805可以采用lcd(liquid crystal display，液晶显示屏)、oled(organic light-emitting diode,有机发光二极管)等材质制备。
[0191]
摄像头组件806用于采集图像或视频。可选地，摄像头组件806包括前置摄像头和后置摄像头。通常，前置摄像头设置在终端的前面板，后置摄像头设置在终端的背面。在一些实施例中，后置摄像头为至少两个，分别为主摄像头、景深摄像头、广角摄像头、长焦摄像头中的任意一种，以实现主摄像头和景深摄像头融合实现背景虚化功能、主摄像头和广角摄像头融合实现全景拍摄以及vr(virtual reality，虚拟现实)拍摄功能或者其它融合拍摄功能。在一些实施例中，摄像头组件806还可以包括闪光灯。闪光灯可以是单温闪光灯，也可以是双温闪光灯。双温闪光灯是指暖光闪光灯和冷光闪光灯的组合，可以用于不同温下的光线补偿。
[0192]
音频电路807可以包括麦克风和扬声器。麦克风用于采集用户及环境的声波，并将声波转换为电信号输入至处理器801进行处理，或者输入至射频电路804以实现语音通信。出于立体声采集或降噪的目的，麦克风可以为多个，分别设置在终端800的不同部位。麦克风还可以是阵列麦克风或全向采集型麦克风。扬声器则用于将来自处理器801或射频电路804的电信号转换为声波。扬声器可以是传统的薄膜扬声器，也可以是压电陶瓷扬声器。当扬声器是压电陶瓷扬声器时，不仅可以将电信号转换为人类可听见的声波，也可以将电信号转换为人类听不见的声波以进行测距等用途。在一些实施例中，音频电路807还可以包括耳机插孔。
[0193]
定位组件808用于定位终端800的当前地理位置，以实现导航或lbs(location based service，基于位置的服务)。定位组件808可以是基于美国的gps(global positioning system，全球定位系统)、中国的北斗系统、俄罗斯的格雷纳斯系统或欧盟的伽利略系统的定位组件。
[0194]
电源809用于为终端800中的各个组件进行供电。电源809可以是交流电、直流电、一次性电池或可充电电池。当电源809包括可充电电池时，该可充电电池可以支持有线充电或无线充电。该可充电电池还可以用于支持快充技术。
[0195]
在一些实施例中，终端800还包括有一个或多个传感器810。该一个或多个传感器810包括但不限于：加速度传感器811、陀螺仪传感器812、压力传感器813、指纹传感器814、光学传感器815以及接近传感器816。
[0196]
加速度传感器811可以检测以终端800建立的坐标系的三个坐标轴上的加速度大小。比如，加速度传感器811可以用于检测重力加速度在三个坐标轴上的分量。处理器801可以根据加速度传感器811采集的重力加速度信号，控制显示屏805以横向视图或纵向视图进行用户界面的显示。加速度传感器811还可以用于游戏或者用户的运动数据的采集。
[0197]
陀螺仪传感器812可以检测终端800的机体方向及转动角度，陀螺仪传感器812可以与加速度传感器811协同采集用户对终端800的3d动作。处理器801根据陀螺仪传感器812采集的数据，可以实现如下功能：动作感应(比如根据用户的倾斜操作来改变ui)、拍摄时的图像稳定、游戏控制以及惯性导航。
[0198]
压力传感器813可以设置在终端800的侧边框和/或显示屏805的下层。当压力传感器813设置在终端800的侧边框时，可以检测用户对终端800的握持信号，由处理器801根据压力传感器813采集的握持信号进行左右手识别或快捷操作。当压力传感器813设置在显示屏805的下层时，由处理器801根据用户对显示屏805的压力操作，实现对ui界面上的可操作性控件进行控制。可操作性控件包括按钮控件、滚动条控件、图标控件、菜单控件中的至少一种。
[0199]
指纹传感器814用于采集用户的指纹，由处理器801根据指纹传感器814采集到的指纹识别用户的身份，或者，由指纹传感器814根据采集到的指纹识别用户的身份。在识别出用户的身份为可信身份时，由处理器801授权该用户执行相关的敏感操作，该敏感操作包括解锁屏幕、查看加密信息、下载软件、支付及更改设置等。指纹传感器814可以被设置终端800的正面、背面或侧面。当终端800上设置有物理按键或厂商logo时，指纹传感器814可以与物理按键或厂商logo集成在一起。
[0200]
光学传感器815用于采集环境光强度。在一个实施例中，处理器801可以根据光学传感器815采集的环境光强度，控制显示屏805的显示亮度。具体地，当环境光强度较高时，调高显示屏805的显示亮度；当环境光强度较低时，调低显示屏805的显示亮度。在另一个实施例中，处理器801还可以根据光学传感器815采集的环境光强度，动态调整摄像头组件806的拍摄参数。
[0201]
接近传感器816，也称距离传感器，通常设置在终端800的前面板。接近传感器816用于采集用户与终端800的正面之间的距离。在一个实施例中，当接近传感器816检测到用户与终端800的正面之间的距离逐渐变小时，由处理器801控制显示屏805从亮屏状态切换为息屏状态；当接近传感器816检测到用户与终端800的正面之间的距离逐渐变大时，由处理器801控制显示屏805从息屏状态切换为亮屏状态。
[0202]
本领域技术人员可以理解，图8中示出的结构并不构成对终端800的限定，可以包括比图示更多或更少的组件，或者组合某些组件，或者采用不同的组件布置。
[0203]
本技术实施例还提供了一种非临时性计算机可读存储介质，当所述存储介质中的指令由终端的处理器执行时，使得终端能够执行上实施例提供的页岩储层的孔隙度确定方法。
[0204]
本技术实施例还提供了一种包含指令的计算机程序产品，当其在终端上运行时，
使得终端执行上述实施例提供的页岩储层的孔隙度确定方法。
[0205]
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成，也可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，上述提到的存储介质可以是只读存储器，磁盘或光盘等。
[0206]
以上所述仅为本技术实施例的较佳实施例，并不用以限制本技术实施例，凡在本技术实施例的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本技术的保护范围之内。

技术特征：

1.一种页岩储层的孔隙度确定方法，其特征在于，所述方法包括：获取针对目标页岩储层的目标单井的测井数据；根据所述测井数据，依次确定所述目标页岩储层内有机碳含量与储层密度之间的第一变化关系，以及所述目标页岩储层内有机碳含量与储层体积之间的第二变化关系；根据所述第一变化关系、所述第二变化关系和所述测井数据，确定消除所述目标页岩储层内有机碳含量影响后的所述目标页岩储层的孔隙度。2.如权利要求1所述的方法，其特征在于，所述测井数据包括自然伽马测井曲线、电阻率曲线、声波时差曲线和密度曲线；所述根据所述测井数据，依次确定所述目标页岩储层内有机碳含量与储层密度之间的第一变化关系，以及所述目标页岩储层内有机碳含量与储层体积之间的第二变化关系，包括：根据所述自然伽马测井曲线、所述电阻率曲线和所述声波时差曲线，确定所述目标页岩储层中的富有机质层段对应的测井曲线段和贫有机质层段对应的测井曲线段；根据所述富有机质层段对应的测井曲线段和所述贫有机质层段对应的测井曲线段，确定所述目标页岩储层的有机碳含量确定模型；将所述密度曲线的倒数与所述有机碳含量模型进行线性拟合，得到所述第一变化关系；根据所述有机碳含量确定模型和所述第一变化关系，确定所述第二变化关系。3.如权利要求2所述的方法，其特征在于，所述根据所述自然伽马测井曲线、所述电阻率曲线和所述声波时差曲线，确定所述目标页岩储层中的富有机质层段对应的测井曲线段和贫有机质层段对应的测井曲线段，包括：将所述电阻率曲线和所述声波时差曲线进行重叠，得到重叠曲线；根据所述自然伽马曲线中的非有机质层段对应的曲线所在位置，将所述重叠曲线中非有机质层段对应的曲线剔除，得到有机质重叠曲线；从所述有机质重叠曲线中确定富有机质层段对应的测井曲线段和贫有机质层段对应的测井曲线段，所述富有机质层段对应的测井曲线段为所述有机质重叠曲线中所述电阻率曲线和所述声波时差曲线走势不相同且不发生重合的测井曲线段，所述贫有机质层段对应的测井曲线段为所述有机质重叠曲线中所述电阻率曲线和所述声波时差曲线走势相同且发生重合的测井曲线段。4.如权利要求2所述的方法，其特征在于，所述根据所述富有机质层段对应的测井曲线段和所述贫有机质层段对应的测井曲线段，确定所述目标页岩储层的有机碳含量确定模型，包括：获取所述目标页岩储层的地层电阻率和不同地层对应的有机碳含量；在所述贫有机质层段对应的测井曲线段中确定声波时差基线和电阻率基线，所述声波时差基线为所述贫有机质层段对应的测井曲线段重叠部分中声波对应的基线，所述电阻率基线为所述贫有机质层段对应的测井曲线段重叠部分中电阻率对应的基线；分别根据所述贫有机质层段对应的声波时差基线和电阻率基线，在所述富有机质层段对应的测井曲线段中确定对应的参考声波时差基线和参考电阻率基线，并确定所述富有机质层段对应的测井曲线段中的参考声波时差曲线；
根据所述参考声波时差基线、所述参考电阻率基线、所述参考声波时差曲线和所述地层电阻率确定不同地层对应的有机质参考丰度；将所述不同地层对应的有机质参考丰度和所述不同地层对应的有机碳含量进行线性拟合，得到所述有机碳含量确定模型。5.如权利要求2所述的方法，其特征在于，所述根据所述有机碳含量确定模型和所述第一变化关系，确定所述第二变化关系，包括：从所述第一变化关系中确定所述有机碳含量为0和100％时分别对应的密度；将所述有机碳含量为0时对应的密度确定为所述目标页岩储层的岩石骨架密度，并将所述有机碳含量为100％时对应的密度确定为所述有机碳含量的密度；根据所述有机碳含量密度、所述有机碳含量确定模型和不同地层对应的参考密度，确定所述第二变化关系，所述不同地层对应的参考密度通过所述密度曲线获取得到。6.如权利要求1所述的方法，其特征在于，所述根据所述第一变化关系和所述第二变化关系，确定消除所述目标页岩储层内有机碳含量影响后的所述目标页岩储层的孔隙度，包括：根据所述第二变化关系、所述第一变化关系中描述的有机碳含量的密度和不同地层对应的参考密度，确定消除所述有机碳含量的影响后所述目标页岩储层的不同地层的地层密度；根据所述不同地层的地层密度、所述岩石骨架密度和所述目标页岩储层内流体的流体密度，确定所述不同地层的孔隙度。7.一种页岩储层的孔隙度确定装置，其特征在于，所述装置包括：获取模块，用于获取针对目标页岩储层的目标单井的测井数据；第一确定模块，用于根据所述测井数据，依次确定所述目标页岩储层内有机碳含量与储层密度之间的第一变化关系，以及所述目标页岩储层内有机碳含量与储层体积之间的第二变化关系；第二确定模块，用于根据所述第一变化关系、所述第二变化关系和所述测井数据，确定消除所述目标页岩储层内有机碳含量影响后的所述目标页岩储层的孔隙度。8.如权利要求7所述的装置，其特征在于，所述测井数据包括自然伽马测井曲线、电阻率曲线、声波时差曲线和密度曲线；所述第一确定模块包括：第一确定子模块，用于根据所述自然伽马测井曲线、所述电阻率曲线和所述声波时差曲线，确定所述目标页岩储层中的富有机质层段对应的测井曲线段和贫有机质层段对应的测井曲线段；第二确定子模块，用于根据所述富有机质层段对应的测井曲线段和所述贫有机质层段对应的测井曲线段，确定所述目标页岩储层的有机碳含量确定模型；拟合子模块，用于将所述密度曲线的倒数与所述有机碳含量模型进行线性拟合，得到所述第一变化关系；第三确定子模块，用于根据所述有机碳含量确定模型和所述第一变化关系，确定所述第二变化关系。9.如权利要求8所述的装置，其特征在于，所述第一确定子模块用于：
将所述电阻率曲线和所述声波时差曲线进行重叠，得到重叠曲线；根据所述自然伽马曲线中的非有机质层段对应的曲线所在位置，将所述重叠曲线中非有机质层段对应的曲线剔除，得到有机质重叠曲线；从所述有机质重叠曲线中确定富有机质层段对应的测井曲线段和贫有机质层段对应的测井曲线段，所述富有机质层段对应的测井曲线段为所述有机质重叠曲线中所述电阻率曲线和所述声波时差曲线走势不相同且不发生重合的测井曲线段，所述贫有机质层段对应的测井曲线段为所述有机质重叠曲线中所述电阻率曲线和所述声波时差曲线走势相同且发生重合的测井曲线段。10.如权利要求8所述的装置，其特征在于，所述第二确定子模块用于：获取所述目标页岩储层的地层电阻率和不同地层对应的有机碳含量；在所述贫有机质层段对应的测井曲线段中确定声波时差基线和电阻率基线，所述声波时差基线为所述贫有机质层段对应的测井曲线段重叠部分中声波对应的基线，所述电阻率基线为所述贫有机质层段对应的测井曲线段重叠部分中电阻率对应的基线；分别根据所述贫有机质层段对应的声波时差基线和电阻率基线，在所述富有机质层段对应的测井曲线段中确定参考声波时差基线和参考电阻率基线，并确定所述富有机质层段对应的测井曲线段中的参考声波时差曲线；根据所述参考声波时差基线、所述参考电阻率曲线、所述参考声波时差曲线和所述地层电阻率确定不同地层对应的有机质参考丰度；将所述不同地层对应的有机质参考丰度和所述不同地层对应的有机碳含量进行线性拟合，得到所述有机碳含量确定模型。11.如权利要求8所述的装置，其特征在于，所述第三确定子模块用于：从所述第一变化关系中确定所述有机碳含量为0和100％时分别对应的密度；将所述有机碳含量为0时对应的密度确定为所述目标页岩储层的岩石骨架密度，并将所述有机碳含量为100％时对应的密度确定为所述有机碳含量的密度；根据所述有机碳含量密度、所述有机碳含量确定模型和不同地层对应的参考密度，确定所述第二变化关系，所述不同地层对应的参考密度通过所述密度曲线获取得到。12.如权利要求7所述的装置，其特征在于，所述第二确定模块包括：第四确定子模块，用于根据所述第二变化关系、所述第一变化关系中描述的有机碳含量的密度和不同地层对应的参考密度，确定消除所述有机碳含量的影响后所述目标页岩储层的不同地层的地层密度；第五确定子模块，用于根据所述不同地层的地层密度、所述岩石骨架密度和所述目标页岩储层内流体的流体密度，确定所述不同地层的孔隙度。13.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有指令，所述指令被处理器执行时实现上述权利要求1至权利要求6中的任一项权利要求所述的方法的步骤。

技术总结

本申请实施例公开了一种页岩储层的孔隙度确定方法、装置及计算机存储介质，属于页岩气勘探开发技术领域。该方法包括：获取针对目标页岩储层的目标单井的测井数据；根据测井数据，依次确定目标页岩储层内有机碳含量与储层密度之间的第一变化关系，以及目标页岩储层内有机碳含量与储层体积之间的第二变化关系；根据第一变化关系、第二变化关系和测井数据，确定消除目标页岩储层内有机碳含量影响后的目标页岩储层的孔隙度。本申请实施例通过目标页岩储层内有机碳含量与储层密度的关系及有机碳含量与储层体积的关系，确定消除目标页岩储层内有机碳含量影响后的目标页岩储层的孔隙度，从而提高了确定页岩储层的孔隙度的准确性。性。性。