一种含气饱和度确定方法和装置

摘要

本公开提供了一种含气饱和度确定方法和装置，包括：对分角度叠加的数据体进行射线弹性阻抗反演得到多个角度的射线弹性阻抗反演体；根据岩性和孔隙度值进行岩相划分；根据提取的井旁道数据及测井岩相进行计算得到AVO阻抗体；对AVO阻抗体的阻抗值域的划分，计算形成地震岩相体；再次进行射线弹性阻抗反演得到储层弹性参数体；在地震岩相体的控制下，对储层弹性参数体进行散点交汇分析，获得对含气饱和度敏感的弹性参数体；提取井旁道数据，进行压实分析形成干岩石骨架体；根据利用干岩石骨架体计算得到的流体体积模块和流体密度，定量计算含气饱和度。本公开实现了排除现有叠前反演含气性预测方法的多解性，准确的预测有效储层中的含气饱和度效果。

权利要求

1.一种含气饱和度确定方法，其特征在于，包括：对分角度叠加的数据体进行射线弹性阻抗反演得到多个角度的射线弹性阻抗反演体；

根据岩性和孔隙度值对多个角度的射线弹性阻抗反演体进行岩相划分，得到测井岩相；

根据提取的第一井旁道数据及测井岩相进行计算得到AVO阻抗体；

对AVO阻抗体的阻抗值域的划分，计算形成地震岩相体；

将多个射线弹性阻抗反演体进行射线弹性阻抗反演得到储层弹性参数体；

在地震岩相体的控制下，对储层弹性参数体进行散点交汇分析，获得对含气饱和度敏感的弹性参数体；

对敏感弹性参数体提取第二井旁道数据，进行压实分析形成干岩石骨架体；

根据利用干岩石骨架体计算得到的流体体积模量和流体密度，定量计算含气饱和度。

2.根据权利要求1所述的含气饱和度确定方法，其特征在于，所述根据岩性和孔隙度值对多个角度的射线弹性阻抗反演体进行岩相划分，得到测井岩相包括：从多个角度的射线弹性阻抗反演体中提取小角度射线弹性阻抗和大角度射线弹性阻抗；

提取小角度射线弹性阻抗中的孔隙度信息及大角度射线弹性阻抗中的岩性信息；

根据岩性和孔隙度值、小角度射线弹性阻抗中的孔隙度信息及大角度射线弹性阻抗中的岩性信息进行岩相划分，得到测井岩相。

3.根据权利要求2所述的含气饱和度确定方法，其特征在于，所述根据提取的第一井旁道数据及划分出的测井岩相进行计算得到AVO阻抗体包括：根据抽取的第一井旁道数据及划分出的测井岩相，进行散点交汇分析；

利用属性投影的方法对小角度射线弹性阻抗体和大角度射线弹性阻抗体进行坐标旋转，得到坐标转化公式；

根据散点交汇分析的结果及坐标转化公式对小角度射线弹性阻抗和大角度射线弹性阻抗进行计算得到AVO阻抗体。

4.根据权利要求1所述的含气饱和度确定方法，其特征在于，所述对敏感弹性参数体提取第二井旁道数据，进行压实分析形成干岩石骨架体包括：对敏感弹性参数体提取第二井旁道数据；

对第二井旁道数据进行压实分析；

根据压实分析的结果获取干岩石骨架的属性变化趋势；

根据干岩石骨架的属性变化趋势进行多属性井插值形成干岩石骨架体。

5.根据权利要求1所述的含气饱和度确定方法，其特征在于，所述对含气饱和度敏感的弹性参数体包括密度和体积模量。

6.一种含气饱和度确定装置，其特征在于，包括：第一反演模块，用于对分角度叠加的数据体进行射线弹性阻抗反演，得到多个角度的射线弹性阻抗反演体；

岩相划分模块，用于根据岩性和孔隙度值对多个角度的射线弹性阻抗反演体进行岩相划分，得到测井岩相；

阻抗体计算模块，用于根据提取的第一井旁道数据及划分出的测井岩相进行计算得到AVO阻抗体；

岩相体计算模块，用于对AVO阻抗体的阻抗值域的划分，计算形成地震岩相体；

第二反演模块，用于将多个射线弹性阻抗反演体进行射线弹性阻抗反演，得到储层弹性参数体；

交汇分析模块，用于在地震岩相体的控制下，对储层弹性参数体进行散点交汇分析，获得对含气饱和度敏感的弹性参数体；

压实分析模块，用于对敏感弹性参数体提取第二井旁道数据，进行压实分析形成干岩石骨架体；

饱和度计算模块，用于根据利用干岩石骨架体计算得到的流体体积模量和流体密度，定量计算含气饱和度。

7.根据权利要求6所述的含气饱和度确定装置，其特征在于，所述岩相划分模块包括：第一提取子模块，用于从多个角度的射线弹性阻抗反演体中提取小角度射线弹性阻抗和大角度射线弹性阻抗；

第二提取子模块，用于提取小角度射线弹性阻抗中的孔隙度信息及大角度射线弹性阻抗中的岩性信息；

划分子模块，用于根据岩性和孔隙度值、小角度射线弹性阻抗中的孔隙度信息及大角度射线弹性阻抗中的岩性信息进行岩相划分，得到测井岩相。

8.根据权利要求7所述的含气饱和度确定装置，其特征在于，所述阻抗体计算模块包括：第一分析子模块，用于根据抽取的第一井旁道数据及划分出的测井岩相，进行散点交汇分析；

坐标旋转子模块，用于利用属性投影的方法对小角度射线弹性阻抗体和大角度射线弹性阻抗体进行坐标旋转，得到坐标转化公式；

计算子模块，用于根据散点交汇分析的结果及坐标转化公式对小角度射线弹性阻抗和大角度射线弹性阻抗进行计算得到AVO阻抗体。

9.根据权利要求6所述的含气饱和度确定装置，其特征在于，所述压实分析模块包括：数据提取子模块，用于对敏感弹性参数体提取第二井旁道数据；

压实分析子模块，用于对第二井旁道数据进行压实分析；

获取子模块，用于根据压实分析的结果获取干岩石骨架的属性变化趋势；

插值子模块，用于根据干岩石骨架的属性变化趋势进行多属性井插值形成干岩石骨架体。

10.根据权利要求6所述的含气饱和度确定装置，其特征在于，所述对含气饱和度敏感的弹性参数体包括密度和体积模量。

说明书

一种含气饱和度确定方法和装置

技术领域

[0001] 本发明涉及地球物理勘探地震解释领域领域，更具体地，涉及一种基于岩相体控制的含气饱和度确定方法及系统。

背景技术

[0002] 二十世纪90年代以来，我国加到了天然气勘探力度，探明储量每年快速稳步增长，但是其探明程度依然很低，仅为6.7％。随着勘探和开发程度的增加，储层评价和含气性预测成为工作中的核心部分。如何从地震数据中获得饱和度、孔隙度等储层物性参数是进行储层孔隙流体预测的关键。目前,利用地震数据进行油气预测基本还是针对储层所含流体的定性预测。

[0003] 然而地下储层岩性的复杂性决定了不能用一种属性或是一个岩石物理方程来描述地下储层分布，通过相关预测方法得到的流体因子或是反演得到的储层物性参数都会受到非有效储层的干扰，预测精度低。

发明内容

[0004] 本公开提出了一种排除非有效储层的干扰、预测精度高的含气饱和度确定方法及系统。

[0005] 根据本公开的一方面，提出了一种含气饱和度确定方法，包括：

[0006] 对分角度叠加的数据体进行射线弹性阻抗反演得到多个角度的射线弹性阻抗反演体；

[0007] 根据岩性和孔隙度值对多个角度的射线弹性阻抗反演体进行岩相划分，得到测井岩相；

[0008] 根据提取的井旁道数据及测井岩相进行计算得到AVO阻抗体；

[0009] 对AVO阻抗体的阻抗值域的划分，计算形成地震岩相体；

[0010] 将多个射线弹性阻抗反演体进行射线弹性阻抗反演得到储层弹性参数体；

[0011] 在地震岩相体的控制下，对储层弹性参数体进行散点交汇分析，获得对含气饱和度敏感的弹性参数体；

[0012] 对敏感弹性参数体提取井旁道数据，进行压实分析形成干岩石骨架体；

[0013] 根据利用干岩石骨架体计算得到的流体体积模量和流体密度，定量计算含气饱和度。

[0014] 根据本公开的第一方面，在本公开的第一种可能的实现方式中，所述根据岩性和孔隙度值对多个角度的射线弹性阻抗反演体进行岩相划分，得到测井岩相包括：

[0015] 从多个角度的射线弹性阻抗反演体中提取小角度射线弹性阻抗和大角度射线弹性阻抗；

[0016] 提取小角度射线弹性阻抗中的孔隙度信息及大角度射线弹性阻抗中的岩性信息；

[0017] 根据岩性和孔隙度值、小角度射线弹性阻抗中的孔隙度信息及大角度射线弹性阻抗中的岩性信息进行岩相划分，得到测井岩相。

[0018] 根据本公开第一方面第一种可能的实现方式，在本公开的第二种可能的实现方式中，所述根据提取的井旁道数据及划分出的测井岩相进行计算得到AVO阻抗体包括：

[0019] 根据抽取的井旁道数据及划分出的测井岩相，进行散点交汇分析；

[0020] 利用属性投影的方法对小角度射线弹性阻抗体和大角度射线弹性阻抗体进行坐标旋转，得到坐标转化公式；

[0021] 根据散点交汇分析的结果及坐标转化公式对小角度射线弹性阻抗和大角度射线弹性阻抗进行计算得到AVO阻抗体。

[0022] 根据本公开的第一方面，在本公开的第三种可能的实现方式中，所述对敏感弹性参数体提取井旁道数据，进行压实分析形成干岩石骨架体包括：

[0023] 对敏感弹性参数体提取井旁道数据；

[0024] 对井旁道数据进行压实分析；

[0025] 根据压实分析的结果获取干岩石骨架的属性变化趋势；

[0026] 根据干岩石骨架的属性变化趋势进行多属性井插值形成干岩石骨架体。

[0027] 根据本公开的第一方面，所述对含气饱和度敏感的弹性参数体包括密度和体积模量。

[0028] 根据本公开的第二方面，一种含气饱和度确定装置，包括：

[0029] 第一反演模块，用于对分角度叠加的数据体进行射线弹性阻抗反演，得到多个角度的射线弹性阻抗反演体；

[0030] 岩相划分模块，用于根据岩性和孔隙度值对多个角度的射线弹性阻抗反演体进行岩相划分，得到测井岩相；

[0031] 阻抗体计算模块，用于根据提取的井旁道数据及划分出的测井岩相进行计算得到AVO阻抗体；

[0032] 岩相体计算模块，用于对AVO阻抗体的阻抗值域的划分，计算形成地震岩相体；

[0033] 第二反演模块，用于将多个射线弹性阻抗反演体进行射线弹性阻抗反演，得到储层弹性参数体；

[0034] 交汇分析模块，用于在地震岩相体的控制下，对储层弹性参数体进行散点交汇分析，获得对含气饱和度敏感的弹性参数体；

[0035] 压实分析模块，用于对敏感弹性参数体提取井旁道数据，进行压实分析形成干岩石骨架体；

[0036] 饱和度计算模块，用于根据利用干岩石骨架体计算得到的流体体积模量和流体密度，定量计算含气饱和度。

[0037] 根据本公开的第二方面，在第二方面第一种可能的实现方式中，所述岩相划分模块包括：

[0038] 第一提取子模块，用于从多个角度的射线弹性阻抗反演体中提取小角度射线弹性阻抗和大角度射线弹性阻抗；

[0039] 第二提取子模块，用于提取小角度射线弹性阻抗中的孔隙度信息及大角度射线弹性阻抗中的岩性信息；

[0040] 划分子模块，用于根据岩性和孔隙度值、小角度射线弹性阻抗中的孔隙度信息及大角度射线弹性阻抗中的岩性信息进行岩相划分，得到测井岩相。

[0041] 根据本公开的第二方面第一种可能的实现方式，在第二方面第二种可能的实现方式中，所述阻抗体计算模块包括：

[0042] 第一分析子模块，用于根据抽取的井旁道数据及划分出的测井岩相，进行散点交汇分析；

[0043] 坐标旋转子模块，用于利用属性投影的方法对小角度射线弹性阻抗体和大角度射线弹性阻抗体进行坐标旋转，得到坐标转化公式；

[0044] 计算子模块，用于根据散点交汇分析的结果及坐标转化公式对小角度射线弹性阻抗和大角度射线弹性阻抗进行计算得到AVO阻抗体。

[0045] 根据本公开的第二方面，在第二方面第三种可能的实现方式中，所述压实分析模块包括：

[0046] 数据提取子模块，用于对敏感弹性参数体提取井旁道数据；

[0047] 压实分析子模块，用于对井旁道数据进行压实分析；

[0048] 获取子模块，用于根据压实分析的结果获取干岩石骨架的属性变化趋势；

[0049] 插值子模块，用于根据干岩石骨架的属性变化趋势进行多属性井插值形成干岩石骨架体。

[0050] 根据本公开的第二方面，在第二方面第四种可能的实现方式中，所述对含气饱和度敏感的弹性参数体包括密度和体积模量。

[0051] 本公开的各方面通过利用射线弹性阻抗交汇形成新的AVO阻抗反演地震岩相体，预测储层有利分布，在岩相体控制下排除非储层影响，进行定量含气饱和度转化，转化过程中，利用Gassmann方程的简化式，将孔隙岩石理解为干燥骨架与孔隙流体两部分，在去除干燥部分，得到的流体密度与含气饱和度之间存在线性关系，实现了排除现有叠前反演含气性预测方法的多解性，准确的预测有效储层中的含气饱和度效果。

附图说明

[0052] 通过结合附图对本公开示例性实施方式进行更详细的描述，本公开的上述以及其它目的、特征和优势将变得更加明显，其中，在本公开示例性实施方式中，相同的参考标号通常代表相同部件。

[0053] 图1示出了本公开一示例性实施例提供的一种含气饱和度确定方法的流程示意图。

[0054] 图2示出了本公开一示例性实施例示出的步骤102的一种流程示意图。

[0055] 图3示出了本公开一示例性实施例示出的步骤103的一种流程示意图。

[0056] 图4示出了本公开一示例性实施例示出的步骤107的一种流程示意图。

[0057] 图5示出了本公开一示例性实施例提供的一种含气饱和度确定装置的结构示意图。

[0058] 图6示出了本公开一示例性实施例示出的岩相划分模块202的一种结构示意图。

[0059] 图7示出了本公开一示例性实施例示出的阻抗体计算模块203的一种结构示意图。

[0060] 图8示出了本公开一示例性实施例示出的压实分析模块207的一种结构示意图。

[0061] 图9示出了本公开涉及的岩石物理模型示意图。

具体实施方式

[0062] 下面将参照附图更详细地描述本公开的优选实施方式。虽然附图中显示了本公开的优选实施方式，然而应该理解，可以以各种形式实现本公开而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了使本公开更加透彻和完整，并且能够将本公开的范围完整地传达给本领域的技术人员。

[0063] 图1示出了根据本发明的一个实施例的含气饱和度确定方法的流程图，如图1所示，该方法包括以下步骤。

[0064] 步骤101，对分角度叠加的数据体进行射线弹性阻抗反演得到多个角度的射线弹性阻抗反演体。

[0065] 所述步骤101在本公开实施例中可具体实现为，由于多个角度的射线弹性阻抗反演体中包括小角度射线弹性阻抗和大角度射线弹性阻抗，而小角度射线弹性阻抗含有孔隙度信息，大角度射线弹性阻抗含有岩性信息，所以可以通过因此可以通过大角度射线弹性阻抗(REI25)和小角度射线弹性阻抗(REI5)区分岩相。角度的数值范围一般看实际资料的角度范围而定，例如一般叠前道集入射角范围为0～30，平均分为三份，小角度为5，中角度为15，大角度为25。

[0066] 步骤102，根据岩性和孔隙度值对多个角度的射线弹性阻抗反演体进行岩相划分，得到测井岩相。

[0067] 步骤103，根据提取的井旁道数据及测井岩相进行计算得到AVO阻抗体。

[0068] 所述步骤103在本公开实施例中可具体实现为，交汇图中散点形状，可以通过测井岩相解释获得，如可令三角形表示泥岩，正方形为干砂岩储层，圆形为砂岩储层等。

[0069] 岩相体的数值范围应该是岩相分类值域(1,2,3…n)。在本实施例中，两个角度的射线弹性阻抗通过坐标转换，形成新的阻抗体(AVOImp)，通过AVOImp与密度交汇图显示可以看出，AVO阻抗在值域上区分储层、干层及泥岩类型，通过阈值划分的方式，可以将AVO阻抗转化成岩性体(定义1—泥岩，2—干层，3—储层)。

[0070] 步骤104，对AVO阻抗体的阻抗值域的划分，计算形成地震岩相体。

[0071] 所述步骤104在本公开实施例中可具体实现为，岩相体LITH与AVO阻抗体应满足如下公式：

[0072]

[0073] value1和value2为区分泥岩的AVOImp的最小值、最大值；value1和value2为区分干层的AVOImp的最小值、最大值；value1和value2为区分储层的AVOImp的最小值、最大值。

[0074] 其中，AVO阻抗与输入的大小角度射线弹性阻抗体满足如下公式：

[0075] AVOImp＝A*X+B*Y+C

[0076] 式中，X是小角度射线弹性阻抗体，Y是大角度射线弹性阻抗体，A、B和C为拟合参数。

[0077] 步骤105，将多个射线弹性阻抗反演体进行射线弹性阻抗反演得到储层弹性参数体。

[0078] 所述步骤105在本公开实施例中可具体实现为，在岩相体的控制下，储层弹性参数可以用Gassmann方程进行描述，如图9所示，当储层为饱和岩石时：

[0079]

[0080] μs＝μd

[0081]

[0082]

[0083] ρs＝(1-φ)ρm+φρf

[0084] 式中：Ks是饱和岩石的体积模量；Kd是干燥岩石骨架的体积模量；Km是岩石固体基质的体积模量；Kf是孔隙流体的体积模量；φ是孔隙度；μs是饱和岩石的剪切模量；μd是干燥岩石骨架的剪切模量；ρs是饱和岩石的密度。

[0085] 这里Gassmann方程孔隙流体只涉及到单相流体。

[0086] 步骤106，在地震岩相体的控制下，对储层弹性参数体进行散点交汇分析，获得对含气饱和度敏感的弹性参数体。

[0087] 所述步骤106在本公开实施例中可具体实现为，所述对含气饱和度敏感的弹性参数体包括密度和体积模量，如果假设流体包含气和水两相，流体密度可写成两相流体的加权和，流体的体积模量可用Wood方程来表示：

[0088]

[0089] 或采用Brie经验流体混合方程，e为经验值，一般取2.5，值越小逼近Wood方程，越大逼近斑块饱和关系流体混合方程：

[0090] Kf＝(Kw-Kg)Sge+Kg

[0091] 流体密度ρf可写为：

[0092] ρf＝(1-Sg)ρw+Sgρg

[0093] 通过以上三个公式可知，当孔隙中含两相流体时，其中一相流体的饱和度可由流体体积模量用非线性公式来表达，流体密度用线性公式来表达。

[0094] 步骤107，对敏感弹性参数体提取井旁道数据，进行压实分析形成干岩石骨架体。

[0095] 所述步骤107在本公开实施例中可具体实现为，对敏感弹性参数体提取井旁道数据，由于Kd、ρd与孔隙度及成岩作用有关，随深度、压力趋势变化，可以通过测井曲线压实分析，形成干岩石骨架体。

[0096] 步骤108，根据利用干岩石骨架体计算得到的流体体积模量和流体密度，定量计算含气饱和度。

[0097] 通过对Gassmann方程重新整理，可以将流体密度和流体体积模量写成与饱和岩石、骨架、基质和孔隙度有关的量：

[0098]

[0099]

[0100] 从上式可以看出，要求得流体密度和流体体积模量，需要知道的物理量有：干岩石的体积模量及密度Kd、ρd，固体基质的体积模量及密度Km、ρm，以及孔隙度φ；由于孔隙度φ随深度变化，在特定物相分类中，可定义为常数；Km、ρm排除与泥质含量的关系，与孔隙度无关，不随深度变化，岩样数据测试及查表获得。

[0101] 本实施例通过利用射线弹性阻抗交汇形成新的AVO阻抗反演地震岩相体，预测储层有利分布，在岩相体控制下排除非储层影响，进行定量含气饱和度转化，转化过程中，利用Gassmann方程的简化式，将孔隙岩石理解为干燥骨架与孔隙流体两部分，在去除干燥部分，得到的流体密度与含气饱和度之间存在线性关系，实现了排除现有叠前反演含气性预测方法的多解性，准确的预测有效储层中的含气饱和度效果。

[0102] 由于小角度射线弹性阻抗含有孔隙度信息，大角度射线弹性阻抗含有岩性信息，因此可以通过大、小角度射线弹性阻抗区分岩相，如图2所示，所述步骤102包括：

[0103] 步骤1021，从多个角度的射线弹性阻抗反演体中提取小角度射线弹性阻抗和大角度射线弹性阻抗；

[0104] 步骤1022，提取小角度射线弹性阻抗中的孔隙度信息及大角度射线弹性阻抗中的岩性信息；

[0105] 步骤1023，根据岩性和孔隙度值、小角度射线弹性阻抗中的孔隙度信息及大角度射线弹性阻抗中的岩性信息进行岩相划分，得到测井岩相。

[0106] 本实施例通过小角度射线弹性阻抗中的孔隙度信息及大角度射线弹性阻抗中的岩性信息进行岩相划分，可以准确的得到测井岩相。

[0107] 为了能够在AVO阻抗上可以排除泥岩和干层的干扰，直接找到孔隙度与AVO阻抗的关系，如图3所示，所述步骤103包括：

[0108] 步骤1031，根据抽取的井旁道数据及划分出的测井岩相，进行散点交汇分析；

[0109] 步骤1032，利用属性投影的方法对小角度射线弹性阻抗体和大角度射线弹性阻抗体进行坐标旋转，得到坐标转化公式；

[0110] 步骤1033，根据散点交汇分析的结果及坐标转化公式对小角度射线弹性阻抗和大角度射线弹性阻抗进行计算得到AVO阻抗体。

[0111] 本实施例通过根据抽取的井旁道数据及划分出的测井岩相，进行散点交汇分析，利用属性投影的方法对小角度射线弹性阻抗体和大角度射线弹性阻抗体进行坐标旋转，得到坐标转化公式，根据散点交汇分析的结果及坐标转化公式对小角度射线弹性阻抗和大角度射线弹性阻抗进行计算得到AVO阻抗体，可以实现利用计算得到的AVO阻抗排除泥岩和干层的干扰，直接找到孔隙度与AVO阻抗的关系。

[0112] 为了能够去除压实作用的岩石骨架，获得反映流体的含气饱和度体，达到定量预测储层含气性目的，如图4所示，所述步骤107包括：

[0113] 步骤1071，对敏感弹性参数体提取井旁道数据；

[0114] 步骤1072，对井旁道数据进行压实分析；

[0115] 步骤1073，根据压实分析的结果获取干岩石骨架的属性变化趋势；

[0116] 步骤1074，根据干岩石骨架的属性变化趋势进行多属性井插值形成干岩石骨架体。

[0117] 本实施例通过对敏感弹性参数体提取井旁道数据，对井旁道数据进行压实分析，根据压实分析的结果获取干岩石骨架的属性变化趋势，根据干岩石骨架的属性变化趋势进行多属性井插值形成干岩石骨架体，实现了去除压实作用的岩石骨架，获得反映流体的含气饱和度体。

[0118] 另外，作为对上述各实施例的实现，本公开实施例还提供了一种含气饱和度确定装置，如图5所示，包括：第一反演模块201、岩相划分模块202、阻抗体计算模块203、岩相体计算模块204、第二反演模块205、交汇分析模块206、压实分析模块207和饱和度计算模块208。

[0119] 第一反演模块201，用于对分角度叠加的数据体进行射线弹性阻抗反演，得到多个角度的射线弹性阻抗反演体。

[0120] 岩相划分模块202，用于根据岩性和孔隙度值对多个角度的射线弹性阻抗反演体进行岩相划分，得到测井岩相。

[0121] 阻抗体计算模块203，用于根据提取的井旁道数据及划分出的测井岩相进行计算得到AVO阻抗体。

[0122] 岩相体计算模块204，用于对AVO阻抗体的阻抗值域的划分，计算形成地震岩相体。

[0123] 第二反演模块205，用于将多个射线弹性阻抗反演体进行射线弹性阻抗反演，得到储层弹性参数体。

[0124] 交汇分析模块206，用于在地震岩相体的控制下，对储层弹性参数体进行散点交汇分析，获得对含气饱和度敏感的弹性参数体，所述对含气饱和度敏感的弹性参数体包括密度和体积模量。

[0125] 压实分析模块207，用于对敏感弹性参数体提取井旁道数据，进行压实分析形成干岩石骨架体；

[0126] 饱和度计算模块208，用于根据利用干岩石骨架体计算得到的流体体积模量和流体密度，定量计算含气饱和度。

[0127] 本实施例通过利用射线弹性阻抗交汇形成新的AVO阻抗反演地震岩相体，预测储层有利分布，在岩相体控制下排除非储层影响，进行定量含气饱和度转化，转化过程中，利用Gassmann方程的简化式，将孔隙岩石理解为干燥骨架与孔隙流体两部分，在去除干燥部分，得到的流体密度与含气饱和度之间存在线性关系，实现了排除现有叠前反演含气性预测方法的多解性，准确的预测有效储层中的含气饱和度效果。

[0128] 由于小角度射线弹性阻抗含有孔隙度信息，大角度射线弹性阻抗含有岩性信息，因此可以通过大、小角度射线弹性阻抗区分岩相，如图6所示，所述岩相划分模块202包括：第一提取子模块2021、第二提取子模块2022和划分子模块2023。

[0129] 第一提取子模块2021，用于从多个角度的射线弹性阻抗反演体中提取小角度射线弹性阻抗和大角度射线弹性阻抗。

[0130] 第二提取子模块2022，用于提取小角度射线弹性阻抗中的孔隙度信息及大角度射线弹性阻抗中的岩性信息。

[0131] 划分子模块2023，用于根据岩性和孔隙度值、小角度射线弹性阻抗中的孔隙度信息及大角度射线弹性阻抗中的岩性信息进行岩相划分，得到测井岩相。

[0132] 本实施例通过小角度射线弹性阻抗中的孔隙度信息及大角度射线弹性阻抗中的岩性信息进行岩相划分，可以准确的得到测井岩相。

[0133] 为了能够在AVO阻抗上可以排除泥岩和干层的干扰，直接找到孔隙度与AVO阻抗的关系，如图7所示，所述阻抗体计算模块203包括：第一分析子模块2031、坐标旋转子模块2032和计算子模块2033。

[0134] 第一分析子模块2031，用于根据抽取的井旁道数据及划分出的测井岩相，进行散点交汇分析。

[0135] 坐标旋转子模块2032，用于利用属性投影的方法对小角度射线弹性阻抗体和大角度射线弹性阻抗体进行坐标旋转，得到坐标转化公式。

[0136] 计算子模块2033，用于根据散点交汇分析的结果及坐标转化公式对小角度射线弹性阻抗和大角度射线弹性阻抗进行计算得到AVO阻抗体。

[0137] 本实施例通过根据抽取的井旁道数据及划分出的测井岩相，进行散点交汇分析，利用属性投影的方法对小角度射线弹性阻抗体和大角度射线弹性阻抗体进行坐标旋转，得到坐标转化公式，根据散点交汇分析的结果及坐标转化公式对小角度射线弹性阻抗和大角度射线弹性阻抗进行计算得到AVO阻抗体，可以实现利用计算得到的AVO阻抗排除泥岩和干层的干扰，直接找到孔隙度与AVO阻抗的关系。

[0138] 为了能够去除压实作用的岩石骨架，获得反映流体的含气饱和度体，达到定量预测储层含气性目的，如图8所示，所述压实分析模块207包括：数据提取子模块2071、压实分析子模块2072、获取子模块2073和插值子模块2074。

[0139] 数据提取子模块2071，用于对敏感弹性参数体提取井旁道数据。

[0140] 压实分析子模块2072，用于对井旁道数据进行压实分析。

[0141] 获取子模块2073，用于根据压实分析的结果获取干岩石骨架的属性变化趋势。

[0142] 插值子模块2074，用于根据干岩石骨架的属性变化趋势进行多属性井插值形成干岩石骨架体。

[0143] 本实施例通过对敏感弹性参数体提取井旁道数据，对井旁道数据进行压实分析，根据压实分析的结果获取干岩石骨架的属性变化趋势，根据干岩石骨架的属性变化趋势进行多属性井插值形成干岩石骨架体，实现了去除压实作用的岩石骨架，获得反映流体的含气饱和度体。

[0144] 应用示例

[0145] 为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。

[0146] 川西某地区蓬莱镇组发育多套气层，展示了很好的勘探潜力。但在开发过程中，发现这种低孔低渗储层气藏气水关系复杂，不同含气丰度气层的地震响应特征差异不明显，导致常规油气检测方法难以解决本区气水预测问题，针对本区地质特点及地震资料现状，利用岩相控制下含气饱和度预测解决了当前的开发难点。

[0147] 研究区定义地震岩相综合考虑了沉积作用和成岩作用、孔隙度，并结合了测井储层评价的标准，对研究区岩相划分出三类岩相(表1)，通过反演获得地震岩相体。通过岩相体划分，将I、II类物相孔隙度定义为

[0148] 表1‐研究区岩相划分标准

[0149]岩相 储层级别 划分标准
I类 1类储层、2类储层 7II类 干层 Por<7
III类 泥岩 Por<7

[0150] 在岩相计算基础上，采用四维交汇岩石物理分析，即密度(横坐标)、纵横波速度比(纵坐标)、孔隙度(颜色)以及地震岩相(符号)四种参数进行交汇分析，通过交会分析可以知道，在有利储层区，在一定孔隙度控制下密度和含气饱和度呈线性关系。

[0151] 利用密度参数作为我们本次预测含气丰度的主要参数。在此基础上将流体密度通过Gassman方程变形式转化为含气饱和度。图中红、黄色代表气层，绿色代表气水区，蓝色为水层，白色为储层不发育段，预测结果与各测试井结果基本吻合，有效验证了方法的适用性。

[0152] 本公开可以是系统、方法和/或计算机程序产品。计算机程序产品可以包括计算机可读存储介质，其上载有用于使处理器实现本公开的各个方面的计算机可读程序指令。

[0153] 计算机可读存储介质可以是可以保持和存储由指令执行设备使用的指令的有形设备。计算机可读存储介质例如可以是――但不限于――电存储设备、磁存储设备、光存储设备、电磁存储设备、半导体存储设备或者上述的任意合适的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括：便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、静态随机存取存储器(SRAM)、便携式压缩盘只读存储器(CD-ROM)、数字多功能盘(DVD)、记忆棒、软盘、机械编码设备、例如其上存储有指令的打孔卡或凹槽内凸起结构、以及上述的任意合适的组合。这里所使用的计算机可读存储介质不被解释为瞬时信号本身，诸如无线电波或者其他自由传播的电磁波、通过波导或其他传输媒介传播的电磁波(例如，通过光纤电缆的光脉冲)、或者通过电线传输的电信号。

[0154] 这里所描述的计算机可读程序指令可以从计算机可读存储介质下载到各个计算/处理设备，或者通过网络、例如因特网、局域网、广域网和/或无线网下载到外部计算机或外部存储设备。网络可以包括铜传输电缆、光纤传输、无线传输、路由器、防火墙、交换机、网关计算机和/或边缘服务器。每个计算/处理设备中的网络适配卡或者网络接口从网络接收计算机可读程序指令，并转发该计算机可读程序指令，以供存储在各个计算/处理设备中的计算机可读存储介质中。

[0155] 用于执行本公开操作的计算机程序指令可以是汇编指令、指令集架构(ISA)指令、机器指令、机器相关指令、微代码、固件指令、状态设置数据、或者以一种或多种编程语言的任意组合编写的源代码或目标代码，所述编程语言包括面向对象的编程语言—诸如Smalltalk、C++等，以及常规的过程式编程语言—诸如“C”语言或类似的编程语言。计算机可读程序指令可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络—包括局域网(LAN)或广域网(WAN)—连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。在一些实施例中，通过利用计算机可读程序指令的状态信息来个性化定制电子电路，例如可编程逻辑电路、现场可编程门阵列(FPGA)或可编程逻辑阵列(PLA)，该电子电路可以执行计算机可读程序指令，从而实现本公开的各个方面。

[0156] 这里参照根据本公开实施例的方法、装置(系统)和计算机程序产品的流程图和/或框图描述了本公开的各个方面。应当理解，流程图和/或框图的每个方框以及流程图和/或框图中各方框的组合，都可以由计算机可读程序指令实现。

[0157] 这些计算机可读程序指令可以提供给通用计算机、专用计算机或其它可编程数据处理装置的处理器，从而生产出一种机器，使得这些指令在通过计算机或其它可编程数据处理装置的处理器执行时，产生了实现流程图和/或框图中的一个或多个方框中规定的功能/动作的装置。也可以把这些计算机可读程序指令存储在计算机可读存储介质中，这些指令使得计算机、可编程数据处理装置和/或其他设备以特定方式工作，从而，存储有指令的计算机可读介质则包括一个制造品，其包括实现流程图和/或框图中的一个或多个方框中规定的功能/动作的各个方面的指令。

[0158] 也可以把计算机可读程序指令加载到计算机、其它可编程数据处理装置、或其它设备上，使得在计算机、其它可编程数据处理装置或其它设备上执行一系列操作步骤，以产生计算机实现的过程，从而使得在计算机、其它可编程数据处理装置、或其它设备上执行的指令实现流程图和/或框图中的一个或多个方框中规定的功能/动作。

[0159] 附图中的流程图和框图显示了根据本公开的多个实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或指令的一部分，所述模块、程序段或指令的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。在有些作为替换的实现中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。

[0160] 以上已经描述了本公开的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场中的技术的技术改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。