用于确定页岩气储层总孔隙度的方法、装置、介质及设备与流程

1.本发明属于非常规油气勘探开发领域，具体涉及一种用于确定页岩气储层总孔隙度的方法、装置、介质及设备。


背景技术：

2.页岩储层中孔隙度的确定是天然气资源量和经济评价的关键参数。页岩的孔隙可简单的分为无机孔隙和有机孔隙，无机孔隙包括微裂缝、黏土孔隙以及碎屑孔隙，有机孔隙主要为纳米级别的孔隙，形成于有机质生烃热演化阶段。
3.目前确定孔隙度的测井方法主要包括核磁测井法、岩心刻度法、基于密度的岩石物理体积模型法和中子密度交会法。其中，丁娱娇等(2014)指出核磁共振测井在页岩气孔隙度评价中具有不可比拟的优势，但是该方法成本较高，并不适合大范围开展测量。岩心刻度法建立在大量分析化验的基础上，由于岩石组分的复杂性，建立的模型往往过于复杂，并且缺乏相应的物理意义；同时由于不同的实验标准对页岩孔隙度测量结果影响较大，差异可达1～3倍(李军等，2017)，使得岩心刻度进行孔隙度的评价由于缺乏统一的实验标准而难以开展。ultey(2005)、mavor(2009)和sondergeld(2010)将页岩简化为无机骨架、有机质和孔隙三部分，根据密度测井的岩石物理体积响应方程进行总孔隙度的计算；但是该方法受井眼质量的影响比较大，同时无机骨架的混合密度随岩相的变化也在发生变化。考虑到单条密度测井曲线的不足，李军等(2017)利用密度测井、中子测井和声波测井，将页岩等效为干黏土、非黏土颗粒骨架、有机质以及孔隙，进而确定总孔隙度，该方法可以较为准确的获取页岩的总孔隙度，但是在页岩水平井测井过程中，由于安全性及经费的考虑，一般没有密度和中子测井曲线，只有伽马测井和声波测井。
4.在页岩测井过程中，伽马测井响应值主要受泥质含量背景值和有机质含量的影响，当有机孔隙的比例发育一致时，扣除泥质含量背景值的影响，可以近似认为伽马测井值与有机孔隙度大小是线性相关的，同时声波孔隙度主要反映的是原生孔隙度。
5.因此，可提出一种利用伽马和声波曲线，适用于页岩气储层总孔隙度的确定方法。


技术实现要素：

6.本发明的目的是提出一种利用伽马和声波曲线，适用于页岩气储层总孔隙度的确定方法。
7.本发明提供一种用于确定页岩气储层总孔隙度的方法，包括：分别获得页岩样品和目标页岩的等效总孔隙度；获得页岩样品的实验总孔隙度；将所述页岩样品的等效总孔隙度和实验总孔隙度做拟合，获得等效总孔隙度和实验总孔隙度的拟合关系；基于所述等效总孔隙度和实验总孔隙度的拟合关系和目标页岩的等效总孔隙度，获得所述目标页岩的实验总孔隙度。
8.可选的，通过以下步骤获得页岩的等效总孔隙度：对页岩的声波时差测井曲线进行有机质含量校正，获得等效声波时差；将所述等效声波时差进行归一化处理，获得等效声
波时差的归一化处理后的值；将页岩的自然伽马曲线进行归一化处理，获得自然伽马曲线归一化处理后的值；基于所述等效声波时差的归一化处理后的值和自然伽马曲线归一化处理后的值，获得等效总孔隙度。
9.可选的，所述等效声波时差为：
10.δt'＝δt-vk(δt
k-δt
ma
)
11.其中，δt'为等效声波时差，δt为页岩的测井声波时差值，vk为有机骨架体积百分数，δtk和δt
ma
分别为有机质骨架和无机骨架的声波时差值。
12.可选的，所述等效声波时差的归一化处理后的值为：
[0013][0014]
其中：φ
ac
为等效声波时差的归一化后的值，反映的是无机孔隙度的相对大小；δt'为等效声波时差；δt'
max
和δt'
min
分别为等效声波时差的上下限值。
[0015]
可选的，所述自然伽马曲线归一化处理后的值为：
[0016][0017]
其中：φ
gr
为自然伽马曲线归一化处理后的值，反映的是有机孔隙度的相对大小；gr为页岩的自然伽马测井值，gr
max
和gr
min
分别为自然伽马测井值的上、下限值。
[0018]
可选的，采用下述公式计算所述等效总孔隙度：
[0019]
φ
ag
＝φ
ac
+φ
gr
[0020]
其中，φ
ag
为等效总孔隙度。
[0021]
可选的，所述等效总孔隙度和实验总孔隙度的拟合关系为：
[0022]
φ
t
＝aφ
ag
+b
[0023]
其中，φ
t
为实验总孔隙度，φ
ag
为等效总孔隙度，a和b为常数。
[0024]
本发明还提供一种电子设备，所述电子设备包括：存储器，存储有可执行指令；处理器，所述处理器运行所述存储器中的所述可执行指令，以实现上述用于确定页岩气储层总孔隙度的方法。
[0025]
本发明还提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现上述用于确定页岩气储层总孔隙度的方法。
[0026]
本发明还提供一种用于确定页岩气储层总孔隙度的装置，包括：等效总孔隙度获取模块，分别获得页岩样品和目标页岩的等效总孔隙度；页岩样品的实验总孔隙度获取模块，获得页岩样品的实验总孔隙度；拟合关系确定模块，将所述页岩样品的等效总孔隙度和实验总孔隙度做拟合，获得等效总孔隙度和实验总孔隙度的拟合关系；目标页岩的实验总孔隙度获取模块，基于所述等效总孔隙度和实验总孔隙度的拟合关系和目标页岩的等效总孔隙度，获得所述目标页岩的实验总孔隙度。
[0027]
可选的，所述等效声波时差为：
[0028]
δt'＝δt-vk(δt
k-δt
ma
)
[0029]
其中，δt'为等效声波时差，δt为页岩的测井声波时差值，vk为有机骨架体积百分数，δtk和δt
ma
分别为有机质骨架和无机骨架的声波时差值。
[0030]
可选的，所述等效声波时差的归一化处理后的值为：
[0031][0032]
其中：φ
ac
为等效声波时差的归一化后的值，反映的是无机孔隙度的相对大小；δt'为等效声波时差；δt'
max
和δt'
min
分别为等效声波时差的上下限值。
[0033]
可选的，所述自然伽马曲线归一化处理后的值为：
[0034][0035]
其中：φ
gr
为自然伽马曲线归一化处理后的值，反映的是有机孔隙度的相对大小；gr为页岩的自然伽马测井值，gr
max
和gr
min
分别为自然伽马测井值的上、下限值。
[0036]
可选的，采用下述公式计算所述等效总孔隙度：
[0037]
φ
ag
＝φ
ac
+φ
gr
[0038]
其中，φ
ag
为等效总孔隙度。
[0039]
可选的，所述等效总孔隙度和实验总孔隙度的拟合关系为：
[0040]
φ
t
＝aφ
ag
+b
[0041]
其中，φ
t
为实验总孔隙度，φ
ag
为等效总孔隙度，a和b为常数。
[0042]
本发明的有益效果在于：本发明的用于确定页岩气储层总孔隙度的方法利用自然伽马曲线和声波曲线这两条较为常见的曲线，简单快捷准确地计算实验总孔隙度，尤其适用于测井难度比较大、测井资料比较少(无中子测井和密度测井)的页岩气水平井中。
[0043]
本发明具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施例中将是显而易见的，或者将在并入本文中的附图和随后的具体实施例中进行详细陈述，这些附图和具体实施例共同用于解释本发明的特定原理。
附图说明
[0044]
通过结合附图对本发明示例性实施方式进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施方式中，相同的参考标号通常代表相同部件。
[0045]
图1示出了根据本发明的一个实施例的一种用于确定页岩气储层总孔隙度的方法的流程图。
[0046]
图2示出了四川盆地页岩气储层实验总孔隙度和等效总孔隙度关系图。
[0047]
图3示出了四川盆地直井a1计算的实验总孔隙度与岩心分析的总孔隙度对比图。
[0048]
图4示出了四川盆地水平井b1计算的实验总孔隙度与利用中子密度交会计算的总孔隙度对比图。
[0049]
图5示出了根据本发明的一个实施例的一种用于确定页岩气储层总孔隙度的装置的框图。
[0050]
附图标记说明
[0051]
102、等效总孔隙度获取模块；104、页岩样品的实验总孔隙度获取模块；106、拟合关系确定模块；108、目标页岩的实验总孔隙度获取模块。
具体实施方式
[0052]
下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明的优选实施方式，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。
[0053]
本发明提供一种用于确定页岩气储层总孔隙度的方法，包括：分别获得页岩样品和目标页岩的等效总孔隙度；获得页岩样品的实验总孔隙度；将页岩样品的等效总孔隙度和实验总孔隙度做拟合，获得等效总孔隙度和实验总孔隙度的拟合关系；基于等效总孔隙度和实验总孔隙度的拟合关系和目标页岩的等效总孔隙度，获得目标页岩的实验总孔隙度。
[0054]
具体的，首先将声波测井进行有机质校正并进行归一化，以消除有机质和无机骨架的共同影响，得到反映无机孔隙相对大小的量；其次将伽马测井进行归一化，以消除无机骨架的影响，得到反映有机孔隙相对大小的量；然后将归一化后两个孔隙度相加，获得等效总孔隙度，并建立等效总孔隙度与岩心分析的实验总孔隙度的定量关系，最终基于该定量关系和目标页岩的等效总孔隙度，实现目标页岩的实验总孔隙度的定量计算。
[0055]
根据示例性的实施方式，用于确定页岩气储层总孔隙度的方法利用自然伽马曲线和声波曲线这两条较为常见的曲线，简单快捷准确地计算实验总孔隙度，尤其适用于测井难度比较大、测井资料比较少(无中子测井和密度测井)的页岩气水平井中。
[0056]
作为可选方案，通过以下步骤获得页岩的等效总孔隙度：对页岩的声波时差测井曲线进行有机质含量校正，获得等效声波时差；将等效声波时差进行归一化处理，获得等效声波时差的归一化处理后的值；将页岩的自然伽马曲线进行归一化处理，获得自然伽马曲线归一化处理后的值；基于等效声波时差的归一化处理后的值和自然伽马曲线归一化处理后的值，获得等效总孔隙度。
[0057]
作为可选方案，等效声波时差为：
[0058]
δt'＝δt-vk(δt
k-δt
ma
)
[0059]
其中，δt'为等效声波时差，δt为页岩的测井声波时差值，vk为有机骨架体积百分数，δtk和δt
ma
分别为有机质骨架和无机骨架的声波时差值。
[0060]
作为可选方案，等效声波时差的归一化处理后的值为：
[0061][0062]
其中：φ
ac
为等效声波时差的归一化后的值，反映的是无机孔隙度的相对大小；δt'为等效声波时差；δt'
max
和δt'
min
分别为等效声波时差的上下限值。
[0063]
具体的，将声波曲线进行有机质校正，得到与该页岩具有同等孔隙度、相同流体和无机骨架的等效声波时差δt'，然后将等效声波时差δt'进行归一化，以消除岩性的影响：
[0064]
δt'＝δt-vk(δt
k-δt
ma
)
ꢀꢀꢀꢀꢀꢀ
(1)
[0065][0066]
式中：φ
ac
为等效声波时差归一化后的值，无量纲，反映的是无机孔隙度的相对大小；vk有机骨架体积百分数，v/v；δtk和δt
ma
分别为有机质骨架和无机骨架的声波时差值，us/ft；δt为页岩的测井声波时差值，us/ft；δt'
max
和δt'
min
分别为等效声波时差的上下
限值，us/ft。
[0067]
作为可选方案，自然伽马曲线归一化处理后的值为：
[0068][0069]
其中：φ
gr
为自然伽马曲线归一化处理后的值，反映的是有机孔隙度的相对大小；gr为页岩的自然伽马测井值，gr
max
和gr
min
分别为自然伽马测井值的上、下限值。
[0070]
具体的，将自然伽马曲线进行归一化，以消除岩性对自然伽马的影响，得到反映有机孔隙度相对大小的值：
[0071][0072]
式中：φ
gr
为自然伽马曲线归一化后的值，无量纲，反映的是有机孔隙度的相对大小。
[0073]
作为可选方案，采用下述公式计算等效总孔隙度：
[0074]
φ
ag
＝φ
ac
+φ
gr
[0075]
其中，φ
ag
为等效总孔隙度。
[0076]
具体的，将等效声波曲线和伽马曲线归一化后的结果进行相加，得到反映等效总孔隙度相对大小的曲线φ
ag
：
[0077]
φ
ag
＝φ
ac
+φ
gr
ꢀꢀꢀꢀꢀꢀꢀꢀ
(4)。
[0078]
作为可选方案，等效总孔隙度和实验总孔隙度的拟合关系为：
[0079]
φ
t
＝aφ
ag
+b
[0080]
其中，φ
t
为实验总孔隙度，φ
ag
为等效总孔隙度，a和b为常数。
[0081]
具体的，对比φ
ag
曲线与岩心分析实验总孔隙度，φ
ag
与岩心分析实验孔隙度存在线性正相关关系。即φ
ag
越大，相应的页岩实验总孔隙度也就越大。
[0082]
φ
t
＝aφ
ag
+b
ꢀꢀꢀꢀꢀꢀꢀ
(5)
[0083]
经实验标定后可得：a＝7.08，b＝4.38。
[0084]
本发明还提供一种电子设备，电子设备包括：存储器，存储有可执行指令；处理器，处理器运行存储器中的可执行指令，以实现上述用于确定页岩气储层总孔隙度的方法。
[0085]
本发明还提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现上述用于确定页岩气储层总孔隙度的方法。
[0086]
本发明还提供一种用于确定页岩气储层总孔隙度的装置，包括：等效总孔隙度获取模块，分别获得页岩样品和目标页岩的等效总孔隙度；页岩样品的实验总孔隙度获取模块，获得页岩样品的实验总孔隙度；拟合关系确定模块，将页岩样品的等效总孔隙度和实验总孔隙度做拟合，获得等效总孔隙度和实验总孔隙度的拟合关系；目标页岩的实验总孔隙度获取模块，基于等效总孔隙度和实验总孔隙度的拟合关系和目标页岩的等效总孔隙度，获得目标页岩的实验总孔隙度。
[0087]
具体的，首先将声波测井进行有机质校正并进行归一化，以消除有机质和无机骨架的共同影响，得到反映无机孔隙相对大小的量；其次将伽马测井进行归一化，以消除无机骨架的影响，得到反映有机孔隙相对大小的量；然后将归一化后两个孔隙度相加，获得等效总孔隙度，并建立等效总孔隙度与岩心分析的实验总孔隙度的定量关系，最终基于该定量
关系和目标页岩的等效总孔隙度，实现目标页岩的实验总孔隙度的定量计算。
[0088]
根据示例性的实施方式，用于确定页岩气储层总孔隙度的装置利用自然伽马曲线和声波曲线这两种较为常见的两条曲线，简单快捷准确地计算实验总孔隙度，尤其适用于测井难度比较大、测井资料比较少(无中子测井和密度测井)的页岩气水平井中。
[0089]
作为可选方案，通过以下步骤获得页岩的等效总孔隙度：对页岩的声波时差测井曲线进行有机质含量校正，获得等效声波时差；将等效声波时差进行归一化处理，获得等效声波时差的归一化处理后的值；将页岩的自然伽马曲线进行归一化处理，获得自然伽马曲线归一化处理后的值；基于等效声波时差的归一化处理后的值和自然伽马曲线归一化处理后的值，获得等效总孔隙度。
[0090]
作为可选方案，等效声波时差为：
[0091]
δt'＝δt-vk(δt
k-δt
ma
)
[0092]
其中，δt'为等效声波时差，δt为页岩的测井声波时差值，vk为有机骨架体积百分数，δtk和δt
ma
分别为有机质骨架和无机骨架的声波时差值。
[0093]
作为可选方案，等效声波时差的归一化处理后的值为：
[0094][0095]
其中：φ
ac
为等效声波时差的归一化后的值，反映的是无机孔隙度的相对大小；δt'为等效声波时差；δt'
max
和δt'
min
分别为等效声波时差的上下限值。
[0096]
具体的，将声波曲线进行有机质校正，得到与该页岩具有同等孔隙度、相同流体和无机骨架的等效声波时差δt'，然后将等效声波时差δt'进行归一化，以消除岩性的影响：
[0097]
δt'＝δt-vk(δt
k-δt
ma
)
ꢀꢀꢀꢀꢀꢀꢀ
(1)
[0098][0099]
式中：φ
ac
为等效声波时差归一化后的值，无量纲，反映的是无机孔隙度的相对大小；vk有机骨架体积百分数，v/v；δtk和δt
ma
分别为有机质骨架和无机骨架的声波时差值，us/ft；δt为页岩的测井声波时差值，us/ft；δt'
max
和δt'
min
分别为等效声波时差的上下限值，us/ft。
[0100]
作为可选方案，自然伽马曲线归一化处理后的值为：
[0101][0102]
其中：φ
gr
为自然伽马曲线归一化处理后的值，反映的是有机孔隙度的相对大小；gr为页岩的自然伽马测井值，gr
max
和gr
min
分别为自然伽马测井值的上、下限值。
[0103]
具体的，将自然伽马曲线进行归一化，以消除岩性对自然伽马的影响，得到反映有机孔隙度相对大小的值：
[0104][0105]
式中：φ
gr
为自然伽马曲线归一化后的值，无量纲，反映的是有机孔隙度的相对大小。
[0106]
作为可选方案，采用下述公式计算等效总孔隙度：
[0107]
φ
ag
＝φ
ac
+φ
gr
[0108]
其中，φ
ag
为等效总孔隙度。
[0109]
具体的，将等效声波曲线和伽马曲线归一化后的结果进行相加，得到反映等效总孔隙度相对大小的曲线φ
ag
：
[0110]
φ
ag
＝φ
ac
+φ
gr
ꢀꢀꢀꢀꢀꢀꢀꢀ
(4)。
[0111]
作为可选方案，等效总孔隙度和实验总孔隙度的拟合关系为：
[0112]
φ
t
＝aφ
ag
+b
[0113]
其中，φ
t
为实验总孔隙度，φ
ag
为等效总孔隙度，a和b为常数。
[0114]
具体的，对比φ
ag
曲线与岩心分析实验总孔隙度，φ
ag
与岩心分析实验孔隙度存在线性正相关关系。即φ
ag
越大，相应的页岩实验总孔隙度也就越大。
[0115]
φ
t
＝aφ
ag
+b
ꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
[0116]
经实验标定后可得：a＝7.08，b＝4.38。
[0117]
实施例一
[0118]
图1示出了根据本发明的一个实施例的一种用于确定页岩气储层总孔隙度的方法的流程图。图2示出了四川盆地页岩气储层实验总孔隙度和等效总孔隙度关系图。图3示出了四川盆地直井a1计算的实验总孔隙度与岩心分析的总孔隙度对比图。图4示出了四川盆地水平井b1计算的实验总孔隙度与利用中子密度交会计算的总孔隙度对比图。
[0119]
结合图1、图2、图3和图4所示，该用于确定页岩气储层总孔隙度的方法，包括：
[0120]
步骤1：分别获得页岩样品和目标页岩的等效总孔隙度；
[0121]
其中，通过以下步骤获得页岩的等效总孔隙度：对页岩的声波时差测井曲线进行有机质含量校正，获得等效声波时差；将等效声波时差进行归一化处理，获得等效声波时差的归一化处理后的值；将页岩的自然伽马曲线进行归一化处理，获得自然伽马曲线归一化处理后的值；基于等效声波时差的归一化处理后的值和自然伽马曲线归一化处理后的值，获得等效总孔隙度。
[0122]
其中，等效声波时差为：
[0123]
δt'＝δt-vk(δt
k-δt
ma
)
[0124]
其中，δt'为等效声波时差，δt为页岩的测井声波时差值，vk为有机骨架体积百分数，δtk和δt
ma
分别为有机质骨架和无机骨架的声波时差值。
[0125]
其中，等效声波时差的归一化处理后的值为：
[0126][0127]
其中：φ
ac
为等效声波时差的归一化后的值，反映的是无机孔隙度的相对大小；δt'为等效声波时差；δt'
max
和δt'
min
分别为等效声波时差的上下限值。
[0128]
其中，自然伽马曲线归一化处理后的值为：
[0129][0130]
其中：φ
gr
为自然伽马曲线归一化处理后的值，反映的是有机孔隙度的相对大小；gr为页岩的自然伽马测井值，gr
max
和gr
min
分别为自然伽马测井值的上、下限值。
[0131]
其中，采用下述公式计算等效总孔隙度：
[0132]
φ
ag
＝φ
ac
+φ
gr
[0133]
其中，φ
ag
为等效总孔隙度。
[0134]
步骤2：获得页岩样品的实验总孔隙度；
[0135]
步骤3：将页岩样品的等效总孔隙度和实验总孔隙度做拟合，获得等效总孔隙度和实验总孔隙度的拟合关系；
[0136]
其中，等效总孔隙度和实验总孔隙度的拟合关系为：
[0137]
φ
t
＝aφ
ag
+b
[0138]
其中，φ
t
为实验总孔隙度，φ
ag
为等效总孔隙度，a和b为常数。
[0139]
步骤4：基于等效总孔隙度和实验总孔隙度的拟合关系和目标页岩的等效总孔隙度，获得目标页岩的实验总孔隙度。
[0140]
以四川盆地的直井a1井为例，对直井a1的总孔隙度进行了解释。如图3所示，图中第一道实线为测井声波时差曲线，虚线为自然伽马曲线；图中第二道虚线为自然电位曲线，实线为井径曲线；第三道为toc质量百分含量曲线；第四道为深度道。
[0141]
对声波时差测井曲线ac进行有机质含量校正，得到不受有机质含量变化影响的等效声波时差曲线，如图中第五道中实线所示；对等效声波时差曲线进行归一化，得到反映无机孔隙度相对大小的曲线，如图中第六道中的实线所示；对自然伽马曲线进行归一化，得到反映有机孔隙度相对大小的曲线，如图中第六道中的虚线所示。然后将归一化后的两条曲线进行相加，得到反映总孔隙度相对大小的曲线，获得等效总孔隙度和实验总孔隙度的拟合关系，最后得到目标页岩的实验总孔隙度曲线，如图中第七道中的实线所示，且第七道中以杆状线表示气测法测量得到的页岩总孔隙度。与岩心分析对比表明，利用本发明提出的页岩总孔隙度的计算方法具有较好的评价精度。
[0142]
以四川盆地的水平井b1井为例，进行了实验总孔隙度的解释。如图4所示，图中第一道虚线为声波时差曲线，实线为自然伽马曲线；第二道虚线为井径曲线；第三道虚线为中子曲线，实线为密度曲线；第四道为toc含量曲线；第五道为深度道。对该井孔隙度的解释可参考直井a1井的解释过程，图中第八道实线为利用新方法计算的页岩实验总孔隙度，虚线为利用中子密度交会的方法计算的页岩总孔隙度，两者具有较好的一致性。
[0143]
从导眼井和水平井的测井评价结果对比中可知，利用本发明提出的总孔隙度的评价方法，不仅具有较高的精度，还具有较强的适用性。
[0144]
实施例二
[0145]
图5示出了根据本发明的一个实施例的一种用于确定页岩气储层总孔隙度的装置的框图。
[0146]
如图5所示，该用于确定页岩气储层总孔隙度的装置，包括：
[0147]
等效总孔隙度获取模块102，分别获得页岩样品和目标页岩的等效总孔隙度；
[0148]
页岩样品的实验总孔隙度获取模块104，获得页岩样品的实验总孔隙度；
[0149]
拟合关系确定模块106，将页岩样品的等效总孔隙度和实验总孔隙度做拟合，获得等效总孔隙度和实验总孔隙度的拟合关系；
[0150]
目标页岩的实验总孔隙度获取模块108，基于等效总孔隙度和实验总孔隙度的拟合关系和目标页岩的等效总孔隙度，获得目标页岩的实验总孔隙度。
[0151]
其中，通过以下步骤获得页岩的等效总孔隙度：对页岩的声波时差测井曲线进行有机质含量校正，获得等效声波时差；将等效声波时差进行归一化处理，获得等效声波时差的归一化处理后的值；将页岩的自然伽马曲线进行归一化处理，获得自然伽马曲线归一化处理后的值；基于等效声波时差的归一化处理后的值和自然伽马曲线归一化处理后的值，获得等效总孔隙度。
[0152]
其中，等效声波时差为：
[0153]
δt'＝δt-vk(δt
k-δt
ma
)
[0154]
其中，δt'为等效声波时差，δ
t
为页岩的测井声波时差值，vk为有机骨架体积百分数，δtk和δt
ma
分别为有机质骨架和无机骨架的声波时差值。
[0155]
其中，等效声波时差的归一化处理后的值为：
[0156][0157]
其中：φ
ac
为等效声波时差的归一化后的值，反映的是无机孔隙度的相对大小；δt'为等效声波时差；δt'
max
和δt'
min
分别为等效声波时差的上下限值。
[0158]
其中，自然伽马曲线归一化处理后的值为：
[0159][0160]
其中：φ
gr
为自然伽马曲线归一化处理后的值，反映的是有机孔隙度的相对大小；gr为页岩的自然伽马测井值，gr
max
和gr
min
分别为自然伽马测井值的上、下限值。
[0161]
其中，采用下述公式计算等效总孔隙度：
[0162]
φ
ag
＝φ
ac
+φ
gr
[0163]
其中，φ
ag
为等效总孔隙度。
[0164]
其中，等效总孔隙度和实验总孔隙度的拟合关系为：
[0165]
φ
t
＝aφ
ag
+b
[0166]
其中，φ
t
为实验总孔隙度，φ
ag
为等效总孔隙度，a和b为常数。
[0167]
实施例三
[0168]
本公开提供一种电子设备包括，该电子设备包括：存储器，存储有可执行指令；处理器，处理器运行存储器中的可执行指令，以实现上述用于确定页岩气储层总孔隙度的方法。
[0169]
根据本公开实施例的电子设备包括存储器和处理器。
[0170]
该存储器用于存储非暂时性计算机可读指令。具体地，存储器可以包括一个或多个计算机程序产品，该计算机程序产品可以包括各种形式的计算机可读存储介质，例如易失性存储器和/或非易失性存储器。该易失性存储器例如可以包括随机存取存储器(ram)和/或高速缓冲存储器(cache)等。该非易失性存储器例如可以包括只读存储器(rom)、硬盘、闪存等。
[0171]
该处理器可以是中央处理单元(cpu)或者具有数据处理能力和/或指令执行能力的其它形式的处理单元，并且可以控制电子设备中的其它组件以执行期望的功能。在本公开的一个实施例中，该处理器用于运行该存储器中存储的该计算机可读指令。
[0172]
本领域技术人员应能理解，为了解决如何获得良好用户体验效果的技术问题，本实施例中也可以包括诸如通信总线、接口等公知的结构，这些公知的结构也应包含在本公开的保护范围之内。
[0173]
有关本实施例的详细说明可以参考前述各实施例中的相应说明，在此不再赘述。
[0174]
实施例四
[0175]
本公开提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现上述用于确定页岩气储层总孔隙度的方法。
[0176]
根据本公开实施例的计算机可读存储介质，其上存储有非暂时性计算机可读指令。当该非暂时性计算机可读指令由处理器运行时，执行前述的本公开各实施例方法的全部或部分步骤。
[0177]
上述计算机可读存储介质包括但不限于：光存储介质(例如：cd－rom和dvd)、磁光存储介质(例如：mo)、磁存储介质(例如：磁带或移动硬盘)、具有内置的可重写非易失性存储器的媒体(例如：存储卡)和具有内置rom的媒体(例如：rom盒)。
[0178]
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。