致密储层中的空间结构的相对含量确定方法、装置及介质与流程

1.本发明属于地球物理勘探技术领域，具体涉及一种致密储层中的空间结构的相对含量确定方法及系统。


背景技术：

2.在致密储层中，由于埋深或者强烈的成岩作用等因素，导致储层孔隙度低，储集空间复杂。brie(1985)等人指出，在声波和地震频段内，岩性变化带来的影响，与孔隙结构(形状和大小)带来的影响相比较小，因此在相同孔隙度下速度变化是由孔隙形状变化引起的这一假设是合理的。储集空间的孔隙结构由孔隙纵横比表示(定义为孔隙的短轴与长轴之比，值域为0-1)，孔隙结构对致密储层发育有较大的影响，例如扁长的裂缝孔隙(纵横比远小于1)比球形孔隙(纵横比接近于1)更能增加储层的渗透能力。
3.以地震资料为基础，在含裂缝致密储层中研究裂缝性相对发育情况已经成为现阶段致密储层描述的关键技术之一。目前准确表征致密储层不同孔隙类型尤其是裂缝的发育强度具有一定的难度，在缺乏成像测井资料的情况下，通过地震岩石物理分析研究获得不同孔隙类型尤其是裂缝孔隙的发育情况对储层的勘探开发具有重要意义。
4.kumar和han(2005)分析了孔隙形状对岩石弹性特征的影响，提出了不同类型孔隙的分类及其所占比例的估测方法，但该方法只能同时反演两种孔隙，而忽略第三种孔隙；xu等(2009)将原始的砂泥岩的xu-white模型推广到碳酸盐岩储层中，反演出了“类裂缝”孔隙和刚性孔隙所占的比例；sain等(2008)提出将改进的xu-white模型反演的孔隙类型进一步细分，包括微孔隙、粒间孔、粒内孔、溶洞、铸模孔、晶间孔和裂缝等，这种构想的可行性和可操作性须进一步验证；孙赞东等(2014)基于纵横波速度反演不同孔隙类型的孔隙形状和相对含量，但其孔隙纵横比的划分须在整段储层上整体考虑。
5.综上，准确表征不同孔隙类型的孔隙结构及其相对含量具有一定的难度，描述储层孔隙类型主要是靠成像测井资料和取芯资料，测试成本较高且不易获得。因此，特别需要一种成本低且能准确表征不同孔隙类型的孔隙结构及其相对含量的方法。


技术实现要素：

6.本发明的目的是提出一种致密储层中的空间结构的相对含量确定方法及系统，解决现有直井体积压裂方法没有将径向井与二氧化碳压裂相结合的问题。
7.为了实现上述目的，本发明提供一种致密储层中的空间结构的相对含量确定方法，包括：步骤1：获取致密碎屑岩储层各矿物组分含量、弹性参数和总孔隙度；步骤2：基于所述各矿物组分含量、弹性参数和总孔隙度，分别计算空间结构全部为粒间孔隙时的岩石等效纵波速度、空间结构全部为次生溶孔时的岩石等效纵波速度和空间结构全部为裂缝时的岩石等效纵波速度；步骤3：基于初始孔隙纵横比，计算给定孔隙纵横比条件下的岩石等效纵波速度；步骤4：分别将空间结构全部为粒间孔隙时的岩石等效纵波速度、空间结构全部为球形溶孔时的岩石等效纵波速度和空间结构全部为裂缝时的岩石等效纵波速度与所
述给定孔隙纵横比条件下的岩石等效纵波速度进行比较；步骤5：根据比较结果，调整所述初始孔隙纵横比，重复执行步骤3-步骤4，直到所述给定孔隙纵横比条件下的岩石等效纵波速度与实际测量的岩石等效纵波速度的误差小于第一预设阈值，将调整后的孔隙纵横比作为致密储层孔隙类型的空间结构参数；步骤6：设置粒间孔隙初始孔隙度、次生溶孔初始孔隙度和裂缝初始孔隙度；步骤7：基于所述粒间孔隙初始孔隙度、次生溶孔初始孔隙度、裂缝初始孔隙度和总孔隙度，获得该孔隙度下的岩石等效纵波速度；步骤8：将该孔隙度下的岩石等效纵波速度与实际测量的岩石等效纵波速度进行比较；步骤9：根据比较结果，调整所述粒间孔隙初始孔隙度、次生溶孔初始孔隙度和裂缝初始孔隙度，重复执行步骤7-步骤8，直到所述比较结果小于第二预设阈值，将调整后的粒间孔隙初始孔隙度、次生溶孔初始孔隙度和裂缝初始孔隙度作为孔隙类型的空间结构的相对含量。
8.可选的，利用时间平均方程计算空间结构全部为粒间孔隙时的岩石等效纵波速度；利用hashin-shtrikman上限公式计算空间结构全部为次生溶孔时的岩石等效纵波速度；利用hashin-shtrikman下限公式计算空间结构全部为裂缝时的岩石等效纵波速度。
9.可选的，利用微分等效介质模型和gassmann方程计算给定孔隙纵横比条件下的岩石等效纵波速度。
10.可选的，根据扫描电镜资料，获取孔隙类型的微观结构参数；基于所述孔隙类型的微观结构参数设定初始孔隙纵横比。
11.可选的，根据测井信息、取芯资料和地质地化信息，获取致密碎屑岩储层各矿物组分含量、弹性参数和总孔隙度。
12.第二方面，本发明还提出了一种电子设备，所述电子设备包括：存储器，存储有可执行指令；处理器，所述处理器运行所述存储器中的所述可执行指令，以实现上述致密储层中的空间结构的相对含量确定方法。
13.第二方面，本发明还提出了一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现上述致密储层中的空间结构的相对含量确定方法。
14.第三方面，本发明还提出了一种致密储层中的空间结构的相对含量确定装置，包括：第一获取模块，获取致密碎屑岩储层各矿物组分含量、弹性参数和总孔隙度；第一计算模块，基于所述各矿物组分含量、弹性参数和总孔隙度，分别计算空间结构全部为粒间孔隙时的岩石等效纵波速度、空间结构全部为次生溶孔时的岩石等效纵波速度和空间结构全部为裂缝时的岩石等效纵波速度；第二计算模块，基于初始孔隙纵横比，计算给定孔隙纵横比条件下的岩石等效纵波速度；第一比较模块，分别将空间结构全部为粒间孔隙时的岩石等效纵波速度、空间结构全部为球形溶孔时的岩石等效纵波速度和空间结构全部为裂缝时的岩石等效纵波速度与所述给定孔隙纵横比条件下的岩石等效纵波速度比较；孔隙类型的空间结构确定模块：根据比较结果，调整所述初始孔隙纵横比，直到所述给定孔隙纵横比条件下的岩石等效纵波速度与实际测量的岩石等效纵波速度的误差小于第一预设阈值，则调整后的孔隙纵横比作为致密储层孔隙类型的空间结构参数；设置模块，设置粒间孔隙初始孔隙度、次生溶孔初始孔隙度和裂缝初始孔隙度；第二获取模块，基于所述粒间孔隙初始孔隙度、次生溶孔初始孔隙度、裂缝初始孔隙度和总孔隙度，获得该孔隙度下的岩石等效纵波速度；第二比较模块，将该孔隙度下的岩石等效纵波速度与实际测量的岩石等效纵波速度进
行比较；空间结构相对含量确定模块，根据比较结果，调整所述粒间孔隙初始孔隙度、次生溶孔初始孔隙度和裂缝初始孔隙度，直到所述比较结果小于第二预设阈值，则获取调整后的粒间孔隙初始孔隙度、次生溶孔初始孔隙度和裂缝初始孔隙度作为孔隙类型的空间结构的相对含量。
15.可选的，利用时间平均方程计算空间结构全部为粒间孔隙时的岩石等效纵波速度；利用hashin-shtrikman上限公式计算空间结构全部为次生溶孔时的岩石等效纵波速度；利用hashin-shtrikman下限公式计算空间结构全部为裂缝时的岩石等效纵波速度。
16.可选的，利用微分等效介质模型和gassmann方程计算给定孔隙纵横比条件下的岩石等效纵波速度。
17.可选的，根据扫描电镜资料，获取孔隙类型的微观结构参数；基于所述孔隙类型的微观结构参数设定初始孔隙纵横比。
18.可选的，根据测井信息、取芯资料和地质地化信息，获取致密碎屑岩储层各矿物组分含量、弹性参数和总孔隙度。
19.本发明的有益效果在于：本发明的致密储层中的空间结构的相对含量确定方法通过模拟优化得到不同孔隙类型的形状，即不同孔隙类型的空间结构参数，以及通过模拟优化得到不同孔隙类型的相对含量，即孔隙度，可以较为准确地获得纵横波速度等信息，较为准确地预测了裂缝的发育位置，对致密油气的勘探开发具有重要的指导意义，且本发明测试成本低。
20.本发明具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施例中将是显而易见的，或者将在并入本文中的附图和随后的具体实施例中进行详细陈述，这些附图和具体实施例共同用于解释本发明的特定原理。
附图说明
21.通过结合附图对本发明示例性实施方式进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施方式中，相同的参考标号通常代表相同部件。
22.图1示出了根据本发明的一个实施例的一种致密储层中的空间结构的相对含量确定方法的流程图。
23.图2示出了根据本发明的一个实施例的一种致密储层中的空间结构的相对含量确定方法的well-a井的测井信息。
24.图3示出了根据本发明的一个实施例的一种致密储层中的空间结构的相对含量确定方法的各种孔隙类型的孔隙度结构参数。
25.图4示出了根据本发明的一个实施例的一种致密储层中的空间结构的相对含量确定方法的各种孔隙类型的孔隙度大小。
26.图5示出了根据本发明的一个实施例的一种致密储层中的空间结构的well-a井部分层段成像测井解释的裂缝孔隙度和岩石物理模拟的裂缝孔隙度图。
27.图6示出了根据本发明的一个实施例的一种致密储层中的空间结构的孔隙类型的空间结构参数及相对含量。
28.图7示出了根据本发明的一个实施例的一种致密储层中的空间结构的相对含量确
定装置的框图。
29.附图标记说明
30.102、第一获取模块；104、第一计算模块；106、第二计算模块；108、第一比较模块；110、孔隙类型的空间结构确定模块；112、设置模块；114、第二获取模块；116、第二比较模块；118、空间结构相对含量确定模块。
具体实施方式
31.下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明的优选实施方式，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。
32.本发明提供一种致密储层中的空间结构的相对含量确定方法，包括：步骤1：获取致密碎屑岩储层各矿物组分含量、弹性参数和总孔隙度；步骤2：基于各矿物组分含量、弹性参数和总孔隙度，分别计算空间结构全部为粒间孔隙时的岩石等效纵波速度、空间结构全部为次生溶孔时的岩石等效纵波速度和空间结构全部为裂缝时的岩石等效纵波速度；步骤3：基于初始孔隙纵横比，计算给定孔隙纵横比条件下的岩石等效纵波速度；步骤4：分别将空间结构全部为粒间孔隙时的岩石等效纵波速度、空间结构全部为球形溶孔时的岩石等效纵波速度和空间结构全部为裂缝时的岩石等效纵波速度与给定孔隙纵横比条件下的岩石等效纵波速度进行比较；步骤5：根据比较结果，调整初始孔隙纵横比，重复执行步骤3-步骤4，直到给定孔隙纵横比条件下的岩石等效纵波速度与实际测量的岩石等效纵波速度的误差小于第一预设阈值，将调整后的孔隙纵横比作为致密储层孔隙类型的空间结构参数；步骤6：设置粒间孔隙初始孔隙度、次生溶孔初始孔隙度和裂缝初始孔隙度；步骤7：基于粒间孔隙初始孔隙度、次生溶孔初始孔隙度、裂缝初始孔隙度和总孔隙度，获得该孔隙度下的岩石等效纵波速度；步骤8：将该孔隙度下的岩石等效纵波速度与实际测量的岩石等效纵波速度进行比较；步骤9：根据比较结果，调整粒间孔隙初始孔隙度、次生溶孔初始孔隙度和裂缝初始孔隙度，重复执行步骤7-步骤8，直到比较结果小于第二预设阈值，将调整后的粒间孔隙初始孔隙度、次生溶孔初始孔隙度和裂缝初始孔隙度作为孔隙类型的空间结构的相对含量。
33.具体的，获得致密碎屑岩储层各矿物组分的含量、弹性参数和总孔隙等；基于获得的各矿物组分的体积含量和弹性参数、总孔隙度等信息，分别获得储集空间全部为残余粒间孔隙时的纵波速度、全部为次生溶孔时的纵波速度、全部为裂缝时的纵波速度；根据获得的总孔隙度φ，分别给定初始孔隙纵横比(粒间孔隙为0.1-0.5，次生溶孔为0.7-1，裂缝为0.005-0.1)，分别利用dem模型和gassmann方程计算给定孔隙纵横比条件下的岩石等效纵波速度v
p,dem
；将定孔隙纵横比条件下的岩石等效纵波速度v
p,dem
分别与时间平均方程、hashin-shtrikman上下限计算的等效纵波速度结果对比，不断调整给定初始孔隙纵横比(孔隙空间结构参数)，直到使得给定孔隙纵横比条件下的岩石等效纵波速度v
p,dem
与测井中实测岩石等效纵波速度误差小于设定的第一预设阈值，则此时的给定初始孔隙纵横比为孔隙空间结构参数；
34.基于获得的各矿物组分的体积含量和弹性参数和总孔隙度等信息，结合孔隙空间
结构参数，分别给定残余粒间孔隙初始孔隙度、次生溶孔初始孔隙度和裂缝的初始孔隙度，计算得到该孔隙度下的岩石等效纵波速度；
35.将该孔隙度下的岩石等效纵波速度与测井中实测岩石纵波等效速度对比，不断调整不同孔隙类型的相对比例(φ
粒间孔隙
φ
次生溶孔
φ
裂缝
＝φ
总
，且不同类型的孔隙大于零并小于总孔隙度)，直到该孔隙度下的岩石等效纵波速度与测井实测岩石纵波等效速度小于一定的误差，则此时可以确定不同孔隙类型的孔隙度。
36.根据示例性的实施方式，致密储层中的空间结构的相对含量确定方法通过模拟优化得到不同孔隙类型的形状，即不同孔隙类型的空间结构参数，以及通过模拟优化得到不同孔隙类型的相对含量，即孔隙度，可以较为准确地获得纵横波速度等信息，较为准确地预测了裂缝的发育位置，对致密油气的勘探开发具有重要的指导意义，且本发明测试成本低。
37.作为可选方案，利用时间平均方程计算空间结构全部为粒间孔隙时的岩石等效纵波速度；利用hashin-shtrikman上限公式计算空间结构全部为次生溶孔时的岩石等效纵波速度；利用hashin-shtrikman下限公式计算空间结构全部为裂缝时的岩石等效纵波速度。
38.具体的，根据各矿物组分的体积含量和弹性参数和总孔隙度，利用时间平均方程计算储集空间全部为残余粒间孔隙，利用hashin-shtrikman上限公式，计算储集空间全部为次生溶孔纵波速度，利用hashin-shtrikman下限公式，计算储集空间全部为裂缝时的纵波速度。
39.作为可选方案，利用微分等效介质模型和gassmann方程计算给定孔隙纵横比条件下的岩石等效纵波速度。
40.具体的，将致密碎屑岩储集空间分为残余粒间孔隙、次生溶孔和裂缝。区别这些孔隙的原则是：粒间孔隙的孔隙纵横比约为0.1-0.5，次生溶孔的孔隙纵横比约为0.7-1，裂缝及裂隙的孔隙纵横比约为0.005-0.1。以dem模型和gassmann流体替换公式计算给定孔隙纵横比条件下的岩石等效纵波速度，属于现有技术，这里不再赘述。
41.(1)微分等效介质模型(dem)
42.norris(1985)等提出的微分等效介质模型通过将包含物分成体积无限小的n份，逐渐向固体基质中添加，直到其含量达到给定体积分数为止的方法对双相介质进行模拟。
[0043][0044][0045]
其中，初始条件为k
*
(0)＝k1，μ
*
(0)＝μ1；k1、k2、k
*
分别表示相1(岩石基质相)、相2(包含物相)以及模型计算得到的等效介质的体积模量；μ1、μ2、μ
*
则分别表示对应的剪切模量；y表示相2(包含物相)的体积分数；p、q是与相2相关的几何参数，用来表征包含物形状对等效模型的影响。
[0046]
(2)gassmann方程
[0047]
gassmann方程是目前最为常用的流体替换(通过干岩石速度预测饱和流体岩石速度)模型，该模型给出了低频假设下饱和流体的岩石等效弹性模量的计算公式：
[0048][0049]
μ
sat
=μ
dry
ꢀꢀ
(4)
[0050]
式中，k
sat
、k
dry
、km、kf分别为饱和岩石、干岩石骨架、混合基质及等效流体的体积模量；μ
sat
、μ
dry
分别是饱和岩石和干岩石骨架的剪切模量；φ是岩石总孔隙度。
[0051]
作为可选方案，根据扫描电镜资料，获取孔隙类型的微观结构参数；基于孔隙类型的微观结构参数设定初始孔隙纵横比。
[0052]
作为可选方案，根据测井信息、取芯资料和地质地化信息，获取致密碎屑岩储层各矿物组分含量、弹性参数和总孔隙度。
[0053]
具体的，结合测井信息、取芯资料及地质地化信息，得到致密碎屑岩储层各矿物组分的含量、弹性参数和总孔隙度。
[0054]
本发明还提供一种电子设备，电子设备包括：存储器，存储有可执行指令；处理器，处理器运行存储器中的可执行指令，以实现上述致密储层中的空间结构的相对含量确定方法。
[0055]
第二方面，本发明还还提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现上述致密储层中的空间结构的相对含量确定方法。
[0056]
本发明还提供一种致密储层中的空间结构的相对含量确定装置，包括：第一获取模块，获取致密碎屑岩储层各矿物组分含量、弹性参数和总孔隙度；第一计算模块，基于各矿物组分含量、弹性参数和总孔隙度，分别计算空间结构全部为粒间孔隙时的岩石等效纵波速度、空间结构全部为次生溶孔时的岩石等效纵波速度和空间结构全部为裂缝时的岩石等效纵波速度；第二计算模块，基于初始孔隙纵横比，计算给定孔隙纵横比条件下的岩石等效纵波速度；第一比较模块，分别将空间结构全部为粒间孔隙时的岩石等效纵波速度、空间结构全部为球形溶孔时的岩石等效纵波速度和空间结构全部为裂缝时的岩石等效纵波速度与给定孔隙纵横比条件下的岩石等效纵波速度比较；孔隙类型的空间结构确定模块：根据比较结果，调整初始孔隙纵横比，直到给定孔隙纵横比条件下的岩石等效纵波速度与实际测量的岩石等效纵波速度的误差小于第一预设阈值，则调整后的孔隙纵横比作为致密储层孔隙类型的空间结构参数；设置模块，设置粒间孔隙初始孔隙度、次生溶孔初始孔隙度和裂缝初始孔隙度；第二获取模块，基于粒间孔隙初始孔隙度、次生溶孔初始孔隙度、裂缝初始孔隙度和总孔隙度，获得该孔隙度下的岩石等效纵波速度；第二比较模块，将该孔隙度下的岩石等效纵波速度与实际测量的岩石等效纵波速度进行比较；空间结构相对含量确定模块，根据比较结果，调整粒间孔隙初始孔隙度、次生溶孔初始孔隙度和裂缝初始孔隙度，直到比较结果小于第二预设阈值，则获取调整后的粒间孔隙初始孔隙度、次生溶孔初始孔隙度和裂缝初始孔隙度作为孔隙类型的空间结构的相对含量。
[0057]
实施例一
[0058]
图1示出了根据本发明的一个实施例的一种致密储层中的空间结构的相对含量确定方法的流程图。图2示出了根据本发明的一个实施例的一种致密储层中的空间结构的相对含量确定方法的well-a井的测井信息。图3示出了根据本发明的一个实施例的一种致密
储层中的空间结构的相对含量确定方法的各种孔隙类型的孔隙度结构参数。图4示出了根据本发明的一个实施例的一种致密储层中的空间结构的相对含量确定方法的各种孔隙类型的孔隙度大小。图5示出了根据本发明的一个实施例的一种致密储层中的空间结构的well-a井部分层段成像测井解释的裂缝孔隙度和岩石物理模拟的裂缝孔隙度图。图6示出了根据本发明的一个实施例的一种致密储层中的空间结构的孔隙类型的空间结构参数及相对含量。
[0059]
结合图1、图2、图3、图4、图5和图6所示，该致密储层中的空间结构的相对含量确定方法，包括：
[0060]
步骤1：获取致密碎屑岩储层各矿物组分含量、弹性参数和总孔隙度；
[0061]
其中，根据测井信息、取芯资料和地质地化信息，获取致密碎屑岩储层各矿物组分含量、弹性参数和总孔隙度。
[0062]
步骤2：基于各矿物组分含量、弹性参数和总孔隙度，分别计算空间结构全部为粒间孔隙时的岩石等效纵波速度、空间结构全部为次生溶孔时的岩石等效纵波速度和空间结构全部为裂缝时的岩石等效纵波速度；
[0063]
其中，利用时间平均方程计算空间结构全部为粒间孔隙时的岩石等效纵波速度；利用hashin-shtrikman上限公式计算空间结构全部为次生溶孔时的岩石等效纵波速度；利用hashin-shtrikman下限公式计算空间结构全部为裂缝时的岩石等效纵波速度。
[0064]
步骤3：基于初始孔隙纵横比，计算给定孔隙纵横比条件下的岩石等效纵波速度；
[0065]
其中，利用微分等效介质模型和gassmann方程计算给定孔隙纵横比条件下的岩石等效纵波速度。
[0066]
其中，根据扫描电镜资料，获取孔隙类型的微观结构参数；基于孔隙类型的微观结构参数设定初始孔隙纵横比。
[0067]
步骤4：分别将空间结构全部为粒间孔隙时的岩石等效纵波速度、空间结构全部为球形溶孔时的岩石等效纵波速度和空间结构全部为裂缝时的岩石等效纵波速度与给定孔隙纵横比条件下的岩石等效纵波速度进行比较；
[0068]
步骤5：根据比较结果，调整初始孔隙纵横比，重复执行步骤3-步骤4，直到给定孔隙纵横比条件下的岩石等效纵波速度与实际测量的岩石等效纵波速度的误差小于第一预设阈值，将调整后的孔隙纵横比作为致密储层孔隙类型的空间结构参数；
[0069]
步骤6：设置粒间孔隙初始孔隙度、次生溶孔初始孔隙度和裂缝初始孔隙度；
[0070]
步骤7：基于粒间孔隙初始孔隙度、次生溶孔初始孔隙度、裂缝初始孔隙度和总孔隙度，获得该孔隙度下的岩石等效纵波速度；
[0071]
步骤8：将该孔隙度下的岩石等效纵波速度与实际测量的岩石等效纵波速度进行比较；
[0072]
步骤9：根据比较结果，调整粒间孔隙初始孔隙度、次生溶孔初始孔隙度和裂缝初始孔隙度，重复执行步骤7-步骤8，直到比较结果小于第二预设阈值，将调整后的粒间孔隙初始孔隙度、次生溶孔初始孔隙度和裂缝初始孔隙度作为孔隙类型的空间结构的相对含量。
[0073]
以四川盆地某致密碎屑岩储层的一口井well-a为例，图2所示为该井的测井信息，包括自然伽马曲线、含水饱和度、孔隙度、渗透率和声波时差曲线。利用已知的测井基础信
息，根据本发明提出的致密储层中不同孔隙的空间结构及其相对含量的确定方法，反演该井不同孔隙类型的孔隙结构和相对含量。
[0074]
图3中(a)为粒间孔隙的空间结构参数(孔隙纵横比)，图3中(b)为次生溶孔的空间结构参数，图3中(c)为裂缝的空间结构参数，可以看出粒间孔隙的孔隙纵横比在0.1近变化，而次生溶孔的孔隙纵横比较大，有的接近于圆孔，裂缝的孔隙纵横比较小，这与我们的地质认识是一致的。
[0075]
图4为反演得到的粒间孔隙、次生溶孔和裂缝的相对含量及总的孔隙度，可以看出储层的孔隙大部分为粒间孔隙和次生溶孔，裂缝孔隙只在局部较为发育，其他地方裂缝孔隙度很小，这与图2中的储层渗透率曲线和声波时差曲线的异常区域相符合，同时验证了本发明方法的有效性。
[0076]
图5为该井成像测井解释的裂缝强度因子和部分成像测井的资料与模拟裂缝孔隙度的对比，其解释的裂缝发育区与本发明的岩石物理模拟方法结果一致。
[0077]
图6为根据本方法确定的不同孔隙的空间结构参数及相对含量，计算得到的该致密碎屑岩裂缝性储层岩石的纵、横波速度，其中实线为估测值，虚线为测井真实值。从图中可以看出，利用反演参数计算得到的岩石纵横波速度与测井真实值之间吻合较好，误差在10％以内，满足应用要求。
[0078]
实施例二
[0079]
图7示出了根据本发明的一个实施例的一种致密储层中的空间结构的相对含量确定装置的框图。
[0080]
如图7所示，该致密储层中的空间结构的相对含量确定装置，包括：
[0081]
第一获取模块102，获取致密碎屑岩储层各矿物组分含量、弹性参数和总孔隙度；
[0082]
其中，根据测井信息、取芯资料和地质地化信息，获取致密碎屑岩储层各矿物组分含量、弹性参数和总孔隙度。
[0083]
第一计算模块104，基于各矿物组分含量、弹性参数和总孔隙度，分别计算空间结构全部为粒间孔隙时的岩石等效纵波速度、空间结构全部为次生溶孔时的岩石等效纵波速度和空间结构全部为裂缝时的岩石等效纵波速度；
[0084]
其中，利用时间平均方程计算空间结构全部为粒间孔隙时的岩石等效纵波速度；利用hashin-shtrikman上限公式计算空间结构全部为次生溶孔时的岩石等效纵波速度；利用hashin-shtrikman下限公式计算空间结构全部为裂缝时的岩石等效纵波速度。
[0085]
第二计算模块106，基于初始孔隙纵横比，计算给定孔隙纵横比条件下的岩石等效纵波速度；
[0086]
其中，利用微分等效介质模型和gassmann方程计算给定孔隙纵横比条件下的岩石等效纵波速度。
[0087]
其中，根据扫描电镜资料，获取孔隙类型的微观结构参数；基于孔隙类型的微观结构参数设定初始孔隙纵横比。
[0088]
第一比较模块108，分别将空间结构全部为粒间孔隙时的岩石等效纵波速度、空间结构全部为球形溶孔时的岩石等效纵波速度和空间结构全部为裂缝时的岩石等效纵波速度与给定孔隙纵横比条件下的岩石等效纵波速度比较；
[0089]
孔隙类型的空间结构确定模块110：根据比较结果，调整初始孔隙纵横比，直到给
定孔隙纵横比条件下的岩石等效纵波速度与实际测量的岩石等效纵波速度的误差小于第一预设阈值，则调整后的孔隙纵横比作为致密储层孔隙类型的空间结构参数；
[0090]
设置模块112，设置粒间孔隙初始孔隙度、次生溶孔初始孔隙度和裂缝初始孔隙度；
[0091]
第二获取模块114，基于粒间孔隙初始孔隙度、次生溶孔初始孔隙度、裂缝初始孔隙度和总孔隙度，获得该孔隙度下的岩石等效纵波速度；
[0092]
第二比较模块116，将该孔隙度下的岩石等效纵波速度与实际测量的岩石等效纵波速度进行比较；
[0093]
空间结构相对含量确定模块118，根据比较结果，调整粒间孔隙初始孔隙度、次生溶孔初始孔隙度和裂缝初始孔隙度，直到比较结果小于第二预设阈值，则获取调整后的粒间孔隙初始孔隙度、次生溶孔初始孔隙度和裂缝初始孔隙度作为孔隙类型的空间结构的相对含量。
[0094]
实施例三
[0095]
本公开提供一种电子设备包括，该电子设备包括：存储器，存储有可执行指令；处理器，处理器运行存储器中的可执行指令，以实现上述致密储层中的空间结构的相对含量确定方法。
[0096]
根据本公开实施例的电子设备包括存储器和处理器。
[0097]
该存储器用于存储非暂时性计算机可读指令。具体地，存储器可以包括一个或多个计算机程序产品，该计算机程序产品可以包括各种形式的计算机可读存储介质，例如易失性存储器和/或非易失性存储器。该易失性存储器例如可以包括随机存取存储器(ram)和/或高速缓冲存储器(cache)等。该非易失性存储器例如可以包括只读存储器(rom)、硬盘、闪存等。
[0098]
该处理器可以是中央处理单元(cpu)或者具有数据处理能力和/或指令执行能力的其它形式的处理单元，并且可以控制电子设备中的其它组件以执行期望的功能。在本公开的一个实施例中，该处理器用于运行该存储器中存储的该计算机可读指令。
[0099]
本领域技术人员应能理解，为了解决如何获得良好用户体验效果的技术问题，本实施例中也可以包括诸如通信总线、接口等公知的结构，这些公知的结构也应包含在本公开的保护范围之内。
[0100]
有关本实施例的详细说明可以参考前述各实施例中的相应说明，在此不再赘述。
[0101]
实施例四
[0102]
本公开提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现上述致密储层中的空间结构的相对含量确定方法。
[0103]
根据本公开实施例的计算机可读存储介质，其上存储有非暂时性计算机可读指令。当该非暂时性计算机可读指令由处理器运行时，执行前述的本公开各实施例方法的全部或部分步骤。
[0104]
上述计算机可读存储介质包括但不限于：光存储介质(例如：cd－rom和dvd)、磁光存储介质(例如：mo)、磁存储介质(例如：磁带或移动硬盘)、具有内置的可重写非易失性存储器的媒体(例如：存储卡)和具有内置rom的媒体(例如：rom盒)。
[0105]
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也
不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。