CN105866003A - 一种页岩气储层总孔隙度获取方法及装置 - Google Patents

Abstract

本发明属于地质勘探技术领域，公开了一种页岩气储层总孔隙度获取方法，包括：获取页岩气储层的岩石密度ρ_matrix；获取页岩气储层的干岩石体积模量K_matrix；获取页岩气储层的干岩石剪切模量μ；依据公式求得声波在页岩气储层颗粒中的传播速度Vp_matrix；依据声波在页岩气储层颗粒中的传播速度Vp_matrix，采用声波计算孔隙度方法，求得页岩气储层总孔隙度φ。本发明提供的方法能够实现高精度高效率的孔隙度测量。

Description

一种页岩气储层总孔隙度获取方法及装置

技术领域

本发明涉及技术领域，特别涉及一种页岩气储层总孔隙度获取方法及装置。

背景技术

页岩气储层总孔隙度是指岩石样品中所有孔隙体积与该岩石样品总体积的比值，或称绝对孔隙度，可用百分比表示。页岩气储层中孔隙多为纳米级孔隙，孔隙类型多样，主要分为有机孔隙、无机孔隙和微裂缝，是页岩气储存的主要空间。在页岩气的勘探开发过程中，页岩气储层总孔隙度是计算页岩气储量和制定开发方案的关键参数。页岩气储层中总孔隙度计算是页岩气储层评价的重要内容，也是优选页岩气有利钻探目标的依据。

目前对页岩气储层总孔隙度的测定，通常采用复杂工艺流程和庞大的设备群组实现。所需花费较大，时间周期较长，稳定性较差；现今还无法通过常规测井资料，低成本快速准确的获取页岩气储层总孔隙度。

发明内容

本发明提供一种页岩气储层总孔隙度获取方法及装置，解决现有技术中页岩气储层总孔隙度测定工艺复杂，周期长，稳定性差，成本高的技术问题。

为解决上述技术问题，本发明提供了一种页岩气储层总孔隙度获取方法，包括：

获取页岩气储层的岩石密度ρ_matrix；

获取页岩气储层的干岩石体积模量K_matrix；

获取页岩气储层的干岩石剪切模量μ；

依据公式求得声波在页岩气储层颗粒中的传播速度Vp_matrix；

依据声波在页岩气储层颗粒中的传播速度Vp_matrix，采用声波计算孔隙度方法，求得页岩气储层总孔隙度

进一步地，所述获取页岩气储层的岩石密度ρ_matrix具体为：

获取所述页岩气储层各矿物的密度ρ_i及各矿物含量f_i；

依据公式求得页岩气储层岩石密度ρ_matrix。

进一步地，所述获取页岩气储层的干岩石体积模量K_matrix具体为：

获取所述页岩气储层各矿物的体积模量K_i及各矿物含量f_i；

依据公式求得干岩石体积模量K_matrix。

进一步地，所述获取页岩气储层的获取页岩气储层的干岩石剪切模量μ具体为：

获取声波在饱水状态下的横波速度Vs_saturated及所述页岩气储层饱水状态下密度ρ_saturated；

依据公式求得干岩石的剪切模量μ。

进一步地，所述声波计算孔隙度方法具体包括：

获取声波在所述页岩气储层的流体中传播的转换系数C；

获取声波在饱水状态下的纵波速度V_p、声波在天然气下的传播速度V_gas；

依据公式求得所述的总孔隙度

一种页岩气储层总孔隙度获取装置，包括：

第一获取单元，获取页岩气储层的岩石密度ρ_matrix、页岩气储层的干岩石体积模量K_matrix以及页岩气储层的干岩石剪切模量μ；

第一计算单元，与所述第一获取单元相连，依据所述页岩气储层的岩石密度ρ_matrix、所述页岩气储层的干岩石体积模量K_matrix以及所述页岩气储层的干岩石剪切模量μ计算声波在页岩气储层颗粒中的传播速Vp_matrix；

声波计算孔隙度单元，与所述第一计算单元相连，获取所述声波在页岩气储层颗粒中的传播速Vp_matrix并采用声波计算孔隙度方法，计算页岩气储层总孔隙度

其中，所述第一计算单元依据公式计算声波在页岩气储层颗粒中的传播速度Vp_matrix。

进一步地，所述第一获取单元包括：

页岩气储层的岩石密度ρ_matrix获取模块，获取所述页岩气储层中各矿物的密度ρ_i及所述页岩气储层中各矿物的含量f_i，并依据公式计算页岩气储层岩石密度ρ_matrix；

页岩气储层的干岩石体积模量K_matrix获取模块，获取所述页岩气储层各矿物的体积模量K_i及各矿物含量f_i，依据公式计算所述干岩石体积模量K_matrix；

页岩气储层的干岩石剪切模量μ获取模块，获取声波在饱水状态下的横波速度Vs_saturated及所述页岩气储层饱水状态下密度ρ_saturated，依据公式计算干岩石的剪切模量μ。

进一步地，所述第一计算单元包括：

声波在页岩气储层颗粒中的传播速Vp_matrix计算模块，分别与所述页岩气储层的岩石密度ρ_matrix获取模块、页岩气储层的干岩石体积模量K_matrix获取模块以及页岩气储层的干岩石剪切模量μ获取模块相连，获取所述页岩气储层的岩石密度ρ_matrix、所述页岩气储层的干岩石体积模量K_matrix以及所述页岩气储层的干岩石剪切模量μ计算声波在页岩气储层颗粒中的传播速Vp_matrix；

其中，所述第声波在页岩气储层颗粒中的传播速Vp_matrix计算模块依据公式计算声波在页岩气储层颗粒中的传播速度Vp_matrix。

进一步地，所述声波计算孔隙度单元包括：

第二获取单元，获取声波在所述页岩气储层的流体中传播的转换系数C、声波在饱水状态下的纵波速度V_p以及声波在天然气下的传播速度V_gas；

总孔隙度计算模块，与所述声波在页岩气储层颗粒中的传播速Vp_matrix计算模块以及所述第二获取单元相连，依据声波在页岩气储层颗粒中的传播速度Vp_matrix、声波在所述页岩气储层的流体中传播的转换系数C、声波在饱水状态下的纵波速度V_p以及声波在天然气下的传播速度V_gas，通过公式计算所述的总孔隙度

进一步地，所述第二获取单元包括：

声波在所述页岩气储层的流体中传播的转换系数C获取模块，获取声波在所述页岩气储层的流体中传播的转换系数C获取模块；

声波在饱水状态下的纵波速度V_p获取模块，获取声波在饱水状态下的纵波速度V_p；

声波在天然气下的传播速度V_gas获取模块，获取声波在天然气下的传播速度V_gas；

其中，所述声波在所述页岩气储层的流体中传播的转换系数C获取模块、所述声波在饱水状态下的纵波速度V_p获取模块以及所述声波在天然气下的传播速度V_gas获取模块分别与所述总孔隙度计算模块相连。

本申请实施例中提供的一个或多个技术方案，至少具有如下技术效果或优点：

本申请实施例中提供的页岩气储层总孔隙度获取方法及装置，通过获取声波在饱水状态下的横波速度Vs_saturated、所述页岩气储层饱水状态下密度ρ_saturated、页岩气储层各矿物的体积模量K_i、各矿物密度ρ_i及各矿物含量f_i等参数并计算出声波在页岩气储层颗粒中的传播速度Vp_matrix；从而实现声波在页岩气储层传播速度Vp_matrix的精确参数测定，并进一步结合声波在所述页岩气储层的流体中传播的转换系数C、纵波速度V_p、声波在天然气下的传播速度V_gas、计算出页岩气储层总孔隙度，大大提升了可靠性和精度。解决了现有技术中，孔隙度测量复杂低效的技术问题；改变了无法通过常规测井资料，低成本快速准确获取页岩气储层总孔隙度的现状。

附图说明

图1为本发明实施例提供的页岩气储层总孔隙度获取方法流程图。

具体实施方式

本申请实施例通过提供一种页岩气储层总孔隙度获取方法及装置，解决现有技术中页岩气储层总孔隙度测定工艺复杂，周期长，稳定性差，成本高的技术问题；达到了大大提升了孔隙度测定的精度，效率和可靠性的技术效果。

为解决上述技术问题，本申请实施例提供技术方案的总体思路如下：

一种页岩气储层总孔隙度获取方法，包括：

获取页岩气储层的岩石密度ρ_matrix；

获取页岩气储层的干岩石体积模量K_matrix；

获取页岩气储层的干岩石剪切模量μ；

依据公式求得声波在页岩气储层颗粒中的传播速度Vp_matrix；

依据声波在页岩气储层颗粒中的传播速度Vp_matrix，采用声波计算孔隙度方法，求得页岩气储层总孔隙度

通过上述内容可以看出，页岩气储层总孔隙度获取方法通过获取：声波在饱水状态下的横波速度Vs_saturated、所述页岩气储层饱水状态下密度ρ_saturated、页岩气储层各矿物的体积模量K_i、各矿物密度ρ_i及各矿物含量f_i等参数，逐步计算出页岩气储层总孔隙度；解决了现今无法通过常规测井资料，低成本快速准确获取页岩气储层总孔隙度，也没有标准计算公式，存在的页岩气开发方案制定难度大的技术问题，有利于制定出较优的页岩气开发方案的技术效果。

为了更好的理解上述技术方案，下面将结合说明书附图以及具体的实施方式对上述技术方案进行详细说明，应当理解本发明实施例以及实施例中的具体特征是对本申请技术方案的详细的说明，而不是对本申请技术方案的限定，在不冲突的情况下，本申请实施例以及实施例中的技术特征可以相互组合。

本发明实施例提供的一种页岩气储层总孔隙度获取方法，包括：

获取页岩气储层的岩石密度ρ_matrix；

获取页岩气储层的干岩石体积模量K_matrix；

获取页岩气储层的干岩石剪切模量μ；

依据公式求得声波在页岩气储层颗粒中的传播速度Vp_matrix；

依据声波在页岩气储层颗粒中的传播速度Vp_matrix，采用声波计算孔隙度方法，求得页岩气储层总孔隙度

下面将分别介绍所述方法。

一种页岩气储层总孔隙度获取方法，包括：

获取页岩气储层的岩石密度ρ_matrix；

获取页岩气储层的干岩石体积模量K_matrix；

获取页岩气储层的干岩石剪切模量μ；

依据公式求得声波在页岩气储层颗粒中的传播速度Vp_matrix；

依据声波在页岩气储层颗粒中的传播速度Vp_matrix，采用声波计算孔隙度方法，求得页岩气储层总孔隙度

所述获取页岩气储层的岩石密度ρ_matrix具体为：

获取所述页岩气储层各矿物的密度ρ_i及各矿物含量f_i；

依据公式求得页岩气储层岩石密度ρ_matrix。

所述获取页岩气储层的干岩石体积模量K_matrix具体为：

获取所述页岩气储层各矿物的体积模量K_i及各矿物含量f_i；

依据公式求得干岩石体积模量K_matrix。

所述获取页岩气储层的获取页岩气储层的干岩石剪切模量μ具体为：

获取声波在饱水状态下的横波速度Vs_saturated及所述页岩气储层饱水状态下密度ρ_saturated；

依据公式求得干岩石的剪切模量μ。

所述声波计算孔隙度方法具体包括：

获取声波在所述页岩气储层的流体中传播的转换系数C；

获取声波在饱水状态下的纵波速度V_p、声波在天然气下的传播速度V_gas；

依据公式求得所述的总孔隙度

一种页岩气储层总孔隙度获取装置，包括：

第一获取单元，获取页岩气储层的岩石密度ρ_matrix、页岩气储层的干岩石体积模量K_matrix以及页岩气储层的干岩石剪切模量μ；

第一计算单元，与所述第一获取单元相连，依据所述页岩气储层的岩石密度ρ_matrix、所述页岩气储层的干岩石体积模量K_matrix以及所述页岩气储层的干岩石剪切模量μ计算声波在页岩气储层颗粒中的传播速Vp_matrix；

声波计算孔隙度单元，与所述第一计算单元相连，获取所述声波在页岩气储层颗粒中的传播速Vp_matrix并采用声波计算孔隙度方法，计算页岩气储层总孔隙度

其中，所述第一计算单元依据公式计算声波在页岩气储层颗粒中的传播速度Vp_matrix。

所述第一获取单元包括：

页岩气储层的岩石密度ρ_matrix获取模块，获取所述页岩气储层中各矿物的密度ρ_i及所述页岩气储层中各矿物的含量f_i，并依据公式计算页岩气储层岩石密度ρ_matrix；

页岩气储层的干岩石体积模量K_matrix获取模块，获取所述页岩气储层各矿物的体积模量K_i及各矿物含量f_i，依据公式计算所述干岩石体积模量K_matrix；

页岩气储层的干岩石剪切模量μ获取模块，获取声波在饱水状态下的横波速度Vs_saturated及所述页岩气储层饱水状态下密度ρ_saturated，依据公式计算干岩石的剪切模量μ。

所述第一计算单元包括：

声波在页岩气储层颗粒中的传播速Vp_matrix计算模块，分别与所述页岩气储层的岩石密度ρ_matrix获取模块、页岩气储层的干岩石体积模量K_matrix获取模块以及页岩气储层的干岩石剪切模量μ获取模块相连，获取所述页岩气储层的岩石密度ρ_matrix、所述页岩气储层的干岩石体积模量K_matrix以及所述页岩气储层的干岩石剪切模量μ计算声波在页岩气储层颗粒中的传播速Vp_matrix；

其中，所述第声波在页岩气储层颗粒中的传播速Vp_matrix计算模块依据公式计算声波在页岩气储层颗粒中的传播速度Vp_matrix。

所述声波计算孔隙度单元包括：

第二获取单元，获取声波在所述页岩气储层的流体中传播的转换系数C、声波在饱水状态下的纵波速度V_p以及声波在天然气下的传播速度V_gas；

总孔隙度计算模块，与所述声波在页岩气储层颗粒中的传播速Vp_matrix计算模块以及所述第二获取单元相连，依据声波在页岩气储层颗粒中的传播速度Vp_matrix、声波在所述页岩气储层的流体中传播的转换系数C、声波在饱水状态下的纵波速度V_p以及声波在天然气下的传播速度V_gas，通过公式计算所述的总孔隙度

所述第二获取单元包括：

声波在所述页岩气储层的流体中传播的转换系数C获取模块，获取声波在所述页岩气储层的流体中传播的转换系数C获取模块；

声波在饱水状态下的纵波速度V_p获取模块，获取声波在饱水状态下的纵波速度V_p；

声波在天然气下的传播速度V_gas获取模块，获取声波在天然气下的传播速度V_gas；

其中，所述声波在所述页岩气储层的流体中传播的转换系数C获取模块、所述声波在饱水状态下的纵波速度V_p获取模块以及所述声波在天然气下的传播速度V_gas获取模块分别与所述总孔隙度计算模块相连。

本实施例还提供两个具体的实施方案。

实施例一

参见图1，在本实施例中，提供了一种利用测井数据计算页岩气储层总孔隙度的方法，包括：

步骤S101，获取所述页岩气储层各矿物的体积模量K_i及各矿物含量f_i；

步骤S102，获取所述页岩气储层各矿物密度ρ_i及各矿物含量f_i；

步骤S103，获取声波在所述页岩气储层饱水状态下的横波速度Vs_saturated、所述页岩气储层饱水状态下密度ρ_saturated；

步骤S104，获取声波在所述页岩气储层的流体中传播的转换系数C、声波在所述页岩气储层饱水状态下的纵波速度V_p、声波在天然气下的传播速度V_gas；

步骤S105，计算出页岩气储层各矿物体积模量K_i及各矿物含量f_i乘积的和求得干岩石体积模量K_matrix；

步骤S106，计算出页岩气储层各矿物密度ρ_i与各矿物含量f_i乘积的和求得页岩气储层干岩石密度ρ_matrix；

步骤S107，计算出Vs_saturated的平方与ρ_saturated的乘积求得干岩石的剪切模量μ；

步骤S108，通过K_matrix、ρ_matrix和μ的值，运用三者的关系计算出声波在页岩气储层颗粒中的传播速度Vp_matrix；

步骤S109，通过C、V_gas、计算出的Vp_matrix和V_p的关系求得所述的总孔隙度

下面详细介绍所述页岩气储层总孔隙度的计算方法的计算原理和计算步骤。

首先，介绍所述页岩气储层总孔隙度的计算方法的计算原理。

本申请提供的计算方法的总体思路是：利用体积模量、剪切模量、密度、矿物组份的关系，推导出声波在页岩气储层干岩石颗粒中的传播速度Vp_matrix，再利用Vp_matrix与声波在所述页岩气储层的流体中传播的转换系数C、纵波速度V_p、声波在天然气中的传播速度V_gas之间的关系，计算出总孔隙度

具体来讲，对于声波在页岩气储层颗粒中的传播速度Vp_matrix的计算，根据Gassmann模型，可知声波在页岩气储层颗粒中的传播速度Vp_matrix与干岩石体积模量K_matrix、干岩石密度ρ_matrix和干岩石剪切模量有着转换关系，即

$V_{P_{matrix}}=\sqrt{\frac{K_{\mathrm{ma}trix}+\frac{4}{3}\mathrm{\μ}}{{\mathrm{\ρ}}_{matrix}}}$

而根据Voigt模型可知，K_matrix＝f₁K₁+f₂K₂+…，ρ_matrix＝f₁ρ₁+f₂ρ₂+…，其中K_i、f_i、ρ_i分别是页岩气储层各矿物体积模量、矿物组份和密度；

而根据剪切模量的定义可知，其中Vs_saturated是声波在所述页岩气储层饱水状态下的横波速度，单位是m/s、ρ_saturated是所述页岩气储层饱水状态下密度；单位是g/cm³；

计算干岩石体积模量K_matrix，我们就要获取所述页岩气储层各矿物的体积模量K_i及各矿物含量f_i；

页岩气储层各矿物的体积模量K_i是通过对所述页岩气储层进行测试，确定了所述页岩气储层矿物组成后，查表获得的，而各矿物含量f_i是通过测试手段获得的，例如页岩气储层矿物类型主要为石英、伊利石、绿泥石、伊蒙混层，通过X衍射技术，我们可以获得相应矿物的百分含量，再通过查找矿物体积模量表，找出相应矿物的体积模量值；

计算干岩石密度ρ_matrix，我们就要获取所述页岩气储层各矿物密度ρ_i及各矿物含量f_i；

页岩气储层各矿物的密度ρ_i及各矿物含量f_i是通过对所述页岩气储层进行测试，确定了所述页岩气储层矿物组成后，查表获得的，例如通过X衍射技术，我们可以获得储层的矿物种类，主要为石英、伊利石、绿泥石、伊蒙混层，再查找矿物密度表，找出相应矿物的密度值；

计算干岩石的剪切模量μ，我们就要获取声波在所述页岩气储层饱水状态下的横波速度Vs_saturated、所述页岩气储层饱水状态下密度ρ_saturated；

声波在所述页岩气储层饱水状态下的横波速度Vs_saturated及页岩气储层饱水状态下密度ρ_saturated是通过常规测井数据获得的；

在计算出干岩石体积模量K_matrix、干岩石密度ρ_matrix和干岩石的剪切模量μ后，我们利用Gassmann模型计算得到声波在页岩气储层颗粒中的传播速度Vp_matrix:

$V_{P_{\mathrm{ma}trix}}=\sqrt{\frac{K_{matrix}+\frac{4}{3}\mathrm{\μ}}{{\mathrm{\ρ}}_{matrix}}}$

将计算得到的Vp_matrix与声波在所述页岩气储层的流体中传播的转换系数C、纵波速度V_p、声波在天然气中的传播速度V_gas联立方程，计算得出所述页岩气储层总孔隙度：

而声波在所述页岩气储层的流体中传播的转换系数C是根据页岩气储层所在深度、温度和压力的不同进行调节变化的，例如储层深度在2230m-2414.5m，压力系数在1.3-1.6左右时对应的C值就在2.5-2.9之间变化，纵波速度V_p可以通过常规测井资料获得，声波在天然气中的传播速度V_gas我们认为是定值340m/s。

为了更好的使用本申请提供的页岩气储层总孔隙度的计算方法，申请人对所述计算方法进行了误差分析，本申请提供方法误差的主要来源有以下几点：

1、由于测试手段和测试精度的原因，页岩气储层矿物含量无法准确获得，计算时，矿物含量小于1％的矿物没有参与计算，但因为这部分矿物含量较少，对干岩石体积模量和密度的计算影响较小，虽会对计算结果带来一定的误差，但误差较小；

2、声波在所述页岩气储层的流体中传播的转换系数C无法实测，取值的影响因素较多，客观取值难度较大，具有不确定性，影响因素包括了储层所在深度、温度和压力等条件，会给计算结果带来一定的误差，但是根据现有数据，通过计算的总孔隙度与某些研究区已有的实测总孔隙度进行相关性比较，可以较为准确的获得C值，故该部分误差也较小。

由上可知，本申请提供的页岩气储层总孔隙度的计算方法，计算参数均能通过可以计算或检测获得，误差也在可以控制的范围内，该方法为计算页岩气储量和制定页岩气开发方案提供了关键参数，降低了页岩气开发方案的制定难度，利于制定出较优的页岩气开发方案。

基于同一发明构思，本申请还提供了实施例二。

实施例二

在本实施例中，提供了一种利用测井数据计算页岩气储层总孔隙度的装置，所述装置包括：

第一获取单元，用于获取所述页岩气储层各矿物的密度ρ_i、各矿物体积模量K_i及各矿物含量f_i、页岩气储层饱水状态下密度ρ_saturated、获取声波在页岩气储层饱水状态下的横波速度Vs_saturated；

第二获取单元，用于获取所述声波在所述页岩气储层的流体中传播的转换系数C、纵波速度V_p、声波在天然气下的传播速度V_gas；

声波在颗粒中传播速度计算单元，用于计算出K_matrix与4/3倍μ的和除以ρ_matrix后的算数平方根求得声波在页岩气储层颗粒中的传播速度Vp_matrix；

总孔隙度计算单元，用于计算出所述页岩气储层总孔隙度

在本申请实施例中，所述第一获取单元还用于，获取所述页岩气储层干岩石体积模量，用于计算出K_i与f_i乘积的和求得干岩石体积模量K_matrix的值；

在本申请实施例中，所述第一获取单元还用于，获取所述页岩气储层干岩石剪切模量，计算出Vs_saturated的平方与ρ_saturated的乘积求得干岩石剪切模量μ的值；

在本申请实施例中，所述第一获取单元还用于，获取所述页岩气储层干岩石密度，计算出ρ_i与f_i乘积的和求得页岩气储层干岩石密度ρ_matrix的值。

本申请实施例中提供的一个或多个技术方案，至少具有如下技术效果或优点：

本申请实施例中提供的页岩气储层总孔隙度获取方法及装置，通过获取声波在饱水状态下的横波速度Vs_saturated、所述页岩气储层饱水状态下密度ρ_saturated、页岩气储层各矿物的体积模量K_i、各矿物密度ρ_i及各矿物含量f_i等参数并计算出声波在页岩气储层颗粒中的传播速度Vp_matrix；从而实现声波在页岩气储层传播速度Vp_matrix的精确参数测定，并进一步结合声波在所述页岩气储层的流体中传播的转换系数C、纵波速度V_p、声波在天然气下的传播速度V_gas、计算出页岩气储层总孔隙度，大大提升了可靠性和精度。解决了现有技术中，孔隙度测量复杂低效的技术问题；改变了无法通过常规测井资料，低成本快速准确获取页岩气储层总孔隙度的现状。

最后所应说明的是，以上具体实施方式仅用以说明本发明的技术方案而非限制，尽管参照实例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。

Claims (10)

1.一种页岩气储层总孔隙度获取方法，其特征在于，包括：

获取页岩气储层的岩石密度ρ_matrix；

获取页岩气储层的干岩石体积模量K_matrix；

获取页岩气储层的干岩石剪切模量μ；

依据公式求得声波在页岩气储层颗粒中的传播速度Vp_matrix；

依据声波在页岩气储层颗粒中的传播速度Vp_matrix，采用声波计算孔隙度方法，求得页岩气储层总孔隙度

2.如权利要求1所述的页岩气储层总孔隙度获取方法，其特征在于，所述获取页岩气储层的岩石密度ρ_matrix具体为：

获取所述页岩气储层各矿物的密度ρ_i及各矿物含量f_i；

依据公式求得页岩气储层岩石密度ρ_matrix。

3.如权利要求1所述的页岩气储层总孔隙度获取方法，其特征在于，所述获取页岩气储层的干岩石体积模量K_matrix具体为：

获取所述页岩气储层各矿物的体积模量K_i及各矿物含量f_i；

依据公式求得干岩石体积模量K_matrix。

4.如权利要求1所述的页岩气储层总孔隙度获取方法，其特征在于，所述获取页岩气储层的获取页岩气储层的干岩石剪切模量μ具体为：

获取声波在饱水状态下的横波速度Vs_saturated及所述页岩气储层饱水状态下密度ρ_saturated；

依据公式求得干岩石的剪切模量μ。

5.如权利要求1～4任一项所述的页岩气储层总孔隙度获取方法，其特征在于，所述声波计算孔隙度方法具体包括：

获取声波在所述页岩气储层的流体中传播的转换系数C；

获取声波在饱水状态下的纵波速度V_p、声波在天然气下的传播速度V_gas；

依据公式求得所述的总孔隙度

6.一种页岩气储层总孔隙度获取装置，其特征在于，包括：

第一获取单元，获取页岩气储层的岩石密度ρ_matrix、页岩气储层的干岩石体积模量K_matrix以及页岩气储层的干岩石剪切模量μ；

第一计算单元，与所述第一获取单元相连，依据所述页岩气储层的岩石密度ρ_matrix、所述页岩气储层的干岩石体积模量K_matrix以及所述页岩气储层的干岩石剪切模量μ计算声波在页岩气储层颗粒中的传播速Vp_matrix；

声波计算孔隙度单元，与所述第一计算单元相连，获取所述声波在页岩气储层颗粒中的传播速Vp_matrix并采用声波计算孔隙度方法，计算页岩气储层总孔隙度

其中，所述第一计算单元依据公式计算声波在页岩气储层颗粒中的传播速度Vp_matrix。

7.如权利要求6所述的页岩气储层总孔隙度获取装置，其特征在于，所述第一获取单元包括：

页岩气储层的岩石密度ρ_matrix获取模块，获取所述页岩气储层中各矿物的密度ρ_i及所述页岩气储层中各矿物的含量f_i，并依据公式计算页岩气储层岩石密度ρ_matrix；

页岩气储层的干岩石体积模量K_matrix获取模块，获取所述页岩气储层各矿物的体积模量K_i及各矿物含量f_i，依据公式计算所述干岩石体积模量K_matrix；

页岩气储层的干岩石剪切模量μ获取模块，获取声波在饱水状态下的横波速度Vs_saturated及所述页岩气储层饱水状态下密度ρ_saturated，依据公式计算干岩石的剪切模量μ。

8.如权利要求7所述的页岩气储层总孔隙度获取装置，其特征在于，所述第一计算单元包括：

声波在页岩气储层颗粒中的传播速Vp_matrix计算模块，分别与所述页岩气储层的岩石密度ρ_matrix获取模块、页岩气储层的干岩石体积模量K_matrix获取模块以及页岩气储层的干岩石剪切模量μ获取模块相连，获取所述页岩气储层的岩石密度ρ_matrix、所述页岩气储层的干岩石体积模量K_matrix以及所述页岩气储层的干岩石剪切模量μ计算声波在页岩气储层颗粒中的传播速Vp_matrix；

其中，所述第声波在页岩气储层颗粒中的传播速Vp_matrix计算模块依据公式计算声波在页岩气储层颗粒中的传播速度Vp_matrix。

9.如权利要求8所述的页岩气储层总孔隙度获取装置，其特征在于，所述声波计算孔隙度单元包括：

第二获取单元，获取声波在所述页岩气储层的流体中传播的转换系数C、声波在饱水状态下的纵波速度V_p以及声波在天然气下的传播速度V_gas；

总孔隙度计算模块，与所述声波在页岩气储层颗粒中的传播速Vp_matrix计算模块以及所述第二获取单元相连，依据声波在页岩气储层颗粒中的传播速度Vp_matrix、声波在所述页岩气储层的流体中传播的转换系数C、声波在饱水状态下的纵波速度V_p以及声波在天然气下的传播速度V_gas，通过公式计算所述的总孔隙度

10.如权利要求9所述的页岩气储层总孔隙度获取装置，其特征在于，所述第二获取单元包括：

声波在所述页岩气储层的流体中传播的转换系数C获取模块，获取声波在所述页岩气储层的流体中传播的转换系数C获取模块；

声波在饱水状态下的纵波速度V_p获取模块，获取声波在饱水状态下的纵波速度V_p；

声波在天然气下的传播速度V_gas获取模块，获取声波在天然气下的传播速度V_gas；

其中，所述声波在所述页岩气储层的流体中传播的转换系数C获取模块、所述声波在饱水状态下的纵波速度V_p获取模块以及所述声波在天然气下的传播速度V_gas获取模块分别与所述总孔隙度计算模块相连。