CN113536218A - 一种计算岩石孔隙体积压缩系数的方法及装置 - Google Patents

Abstract

本申请提供一种计算岩石孔隙体积压缩系数的方法和装置，通过统计资源储层的基本参数，使用该基本参数计算岩石孔隙体积压缩系数，资源储层的基本参数易于获取，降低了岩石孔隙体积压缩系数计算的难度，使该方法能够广泛应用。另外，在计算岩石孔隙体积压缩系数时，考虑了有效应力的变化，通过计算岩石在零有效应力下孔隙体积压缩系数，进而根据不同的有效应力和该岩石在零有效应力下孔隙体积压缩系数，得到不同有效应力下孔隙体积压缩系数，该方法适用于任意的储层，应用更加广泛。通过统一的方法计算不同储层的岩石孔隙体积压缩系数，使得储层弹性能量评价、油气藏储量计算和压力不稳定分析等资源勘探更加准确且易于实现，降低了勘探开发成本。

Description

一种计算岩石孔隙体积压缩系数的方法及装置

技术领域

本发明涉及石油勘探开发技术，尤其涉及一种计算岩石孔隙体积压缩系数的方法及装置。

背景技术

石油勘探是指为了寻找和查明油气资源，而利用各种勘探手段了解地下的地质状况，搞清油气层情况和产出能力的过程。钻探法是油气勘探中必须采用的重要手段，该手段用于调查、发现油气藏一直到油气藏开采的过程。储层岩石孔隙压缩系数作为油气藏工程中的一个重要参数，经常被用于储层弹性能量评价、油气藏储量计算和压力不稳定分析。

目前，确定岩石孔隙体积压缩系数的测量方法主要有实验测量和公式计算两种方法。实验测量方法需要测量不同方向应力变化下的岩石视体积和岩石骨架体积变化量，根据该变化量进一步计算岩石孔隙体积压缩系数；公式计算方法大都为根据实验室测定数据，采用统计学方法建立的经验相关式。并且目前在油气藏工程计算中，往往将岩石孔隙体积压缩系数视为一个常数。

实验测量方法耗时长，成本高，应用受限；而经验相关式测量方法由于影响岩石压缩系数的因素较多，不具有普遍的适用性，例如适用于相对高孔隙度储层的方法当用在低孔隙度储层会出现较大误差。所以，已有的岩石孔隙体积压缩系数的测量方法都具有一定的局限性，无法得到广泛的应用。

发明内容

本申请提供一种计算岩石孔隙体积压缩系数的方法和装置，解决了可在在不同有效应力变化的情况下计算资源储层的岩石孔隙体积压缩系数问题。

一方面，本申请提供一种计算岩石孔隙体积压缩系数的方法，包括：

统计资源储层的基本参数，所述基本参数包括：岩石孔隙度、岩石的泊松比、岩石视体积的杨氏模量、孔隙压力随有效应力变化的衰减系数和有效应力；

根据所述岩石孔隙度、所述岩石的泊松比和所述岩石视体积的杨氏模量，计算得到岩石骨架体积压缩系数，所述岩石骨架体积压缩系数用于表示单位压力变化时引起的岩石骨架体积的变化率；

根据所述岩石的泊松比和所述岩石视体积的杨氏模量，计算得到岩石视体积压缩系数，所述岩石视体积压缩系数用于表示单位压力变化时引起的岩石视体积的变化率；

根据所述岩石孔隙度、所述岩石骨架体积压缩系数和所述岩石视体积压缩系数，计算得到在零有效应力下岩石孔隙体积压缩系数，所述岩石孔隙体积压缩系数用于表示单位压力变化时引起的岩石孔隙体积的变化率；

根据所述孔隙压力随有效应力变化的衰减系数、所述有效应力和所述在零有效应力下岩石孔隙体积压缩系数，计算得到在所述有效应力下的岩石孔隙体积压缩系数。

可选的，所述根据所述岩石孔隙度、所述岩石的泊松比和岩石视体积的杨氏模量，计算得到岩石骨架体积压缩系数，包括：

根据如下公式计算所述岩石骨架体积压缩系数：

其中，C_s为所述岩石骨架体积压缩系数，ν是所述岩石的泊松比，E为所述岩石视体积的杨氏模量，φ为所述岩石孔隙度。

可选的，所述根据所述岩石的泊松比和所述岩石视体积的杨氏模量，计算得到岩石视体积压缩系数，包括：

根据如下公式计算所述岩石视体积压缩系数：

其中，C_b为所述岩石视体积压缩系数，v是所述岩石的泊松比，E为所述岩石视体积杨氏模量。

可选的，所述根据所述岩石孔隙度、所述岩石骨架体积压缩系数和所述岩石视体积压缩系数，计算得到在零有效应力下岩石孔隙体积压缩系数，包括：

根据如下公式计算所述在零有效应力下岩石孔隙体积压缩系数：

其中，C_p0为所述零有效应力下岩石孔隙体积压缩系数，C_b为所述岩石视体积压缩系数，C_s为所述岩石骨架体积压缩系数，φ为所述岩石孔隙度。

可选的，所述根据所述孔隙压力随有效应力变化的衰减系数、所述有效应力和所述在零有效应力下岩石孔隙体积压缩系数，计算得到在所述有效应力下的岩石孔隙体积压缩系数，包括：

根据如下公式计算在所述有效应力下的岩石孔隙体积压缩系数：

其中，C_p为在所述有效应力下岩石孔隙体积压缩系数，C_p0为所述零有效应力下岩石孔隙体积压缩系数，γ为所述孔隙压力随有效应力变化的衰减系数，Δσ_eff为所述有效应力。

本申请第二方面提供一种计算岩石孔隙体积压缩系数的装置，包括：

统计模块，用于统计资源储层的基本参数，所述基本参数包括：岩石孔隙度、岩石的泊松比、岩石视体积的杨氏模量、孔隙压力随有效应力变化的衰减系数和有效应力；

第一计算模块，用于根据所述岩石孔隙度、所述岩石的泊松比和岩石视体积的杨氏模量，计算得到岩石骨架体积压缩系数，所述岩石骨架体积压缩系数用于表示单位压力变化时引起的岩石骨架体积的变化率；

第二计算模块，用于根据所述岩石的泊松比和所述岩石视体积的杨氏模量，计算得到岩石视体积压缩系数，所述岩石视体积压缩系数用于表示单位压力变化时引起的岩石视体积的变化率；

第三计算模块，用于根据所述岩石孔隙度、所述岩石骨架体积压缩系数和所述岩石视体积压缩系数，计算得到在零有效应力下岩石孔隙体积压缩系数，所述岩石孔隙体积压缩系数用于表示单位压力变化时引起的岩石孔隙体积的变化率；

第四计算模块，用于根据所述孔隙压力随有效应力变化的衰减系数、所述有效应力和所述在零有效应力下岩石孔隙体积压缩系数，计算得到在所述有效应力下的岩石孔隙体积压缩系数。

可选的，所述第三计算模块具体用于：

根据如下公式计算所述在零有效应力下岩石孔隙体积压缩系数：

其中，C_p0为所述零有效应力下岩石孔隙体积压缩系数，C_b为所述岩石视体积压缩系数，C_s为所述岩石骨架体积压缩系数，φ为所述岩石孔隙度。

可选的，所述第四计算模块具体用于：

根据如下公式计算在所述有效应力下的岩石孔隙体积压缩系数：

其中，C_p为在所述有效应力下岩石孔隙体积压缩系数，C_p0为所述零有效应力下岩石孔隙体积压缩系数，γ为所述孔隙压力随有效应力变化的衰减系数，Δσ_eff为所述有效应力。

第三方面，本发明提供一种电子设备，包括：至少一个处理器和存储器；

所述存储器存储计算机执行指令；

所述至少一个处理器执行所述存储器存储的计算机执行指令，所述至少一个处理器执行所述存储器存储的计算机执行指令，使得所述至少一个处理器执行如本发明第一方面所述的方法。

第四方面，本发明提供一种计算机可读存储介质，所述计算机可读存储介质中存储有计算机执行指令，所述计算机执行指令被处理器执行时用于实现如本发明第一方面所述的方法。

第五方面，本发明提供一种计算机程序产品，包括计算机程序，所述计算机程序被处理器执行时，实现如本发明第一方面所述的方法。

本申请提供的计算岩石孔隙体积压缩系数的方法和装置，通过统计资源储层的基本参数，使用资源储层的基本参数计算岩石孔隙体积压缩系数，资源储层的基本参数易于获取，从而降低了岩石孔隙体积压缩系数计算的难度和成本等，使得该计算方法能够广泛应用。另外，在计算岩石孔隙体积压缩系数时，考虑了有效应力的变化，通过计算在零有效应力下岩石孔隙体积压缩系数，进而根据不同的有效应力和该在零有效应力下岩石孔隙体积压缩系数，可以得到不同有效应力下岩石孔隙体积压缩系数，该计算方式可以适用于任意的储层，使得应用更加广泛。通过统一的方法计算不同储层的岩石孔隙体积压缩系数的方法，使得储层弹性能量评价、油气藏储量计算和压力不稳定分析等资源勘探开发方法更加准确且易于实现，降低了勘探开发成本。

附图说明

此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本申请的实施例，并与说明书一起用于解释本申请的原理。

图1为本发明实施例一提供的计算岩石孔隙体积压缩系数的方法的流程示意图；

图2为本发明实施例二提供的计算岩石孔隙体积压缩系数的装置的结构示意图；

图3为本发明实施例三提供的一种电子设备的一种结构示意图。

通过上述附图，已示出本申请明确的实施例，后文中将有更详细的描述。这些附图和文字描述并不是为了通过任何方式限制本申请构思的范围，而是通过参考特定实施例为本领域技术人员说明本申请的概念。

具体实施方式

这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本申请相一致的所有实施方式。相反，它们仅是与如所附权利要求书中所详述的、本申请的一些方面相一致的装置和方法的例子。

首先对本申请所涉及的名词进行解释：

岩石孔隙度：是指岩石的孔隙体积与岩石的表面体积的比值，岩石的孔隙是指岩石中存在的未被固态物质填满的空洞，孔隙度是计算储量和评价油藏特性的一个重要指标，岩石孔隙度是一个分数，通常用符号φ表示。

岩石的泊松比：是指岩石在纵向荷载作用下产生的横向应变与纵向应变的比值，反映材料横向变形的弹性常数，岩石的泊松比是评价岩石变形特征的重要参数之一，岩石的泊松比是一个无量纲量，通常用符号ν表示。

岩石视体积的杨氏模量：是指描述固体材料抵抗形变能力的物理量，定义为理想材料有小形变时应力与相应的应变之比，应力是指单位面积上所受到的力，应变是指在外力作用下的相对形变，杨氏模量的大小标志了材料的刚性，杨氏模量越大，越不容易发生形变，岩石视体积的杨氏模量的单位是MPa，通常用符号E表示。

有效应力：是指物体由于外因(受力、湿度、温度场变化等)而变形时，在物体内各部分之间产生相互作用的内力，以抵抗这种外因的作用，并试图使物体从变形后的位置恢复到变形前的位置，这种内力称为有效应力，有效应力的单位是MPa，通常用符号Δσ_eff表示。

岩石骨架体积压缩系数：是指单位压力变化时引起的岩石骨架体积的变化率，岩石骨架体积压缩系数的单位是MPa^-1，通常用符号C_s表示。

岩石视体积压缩系数：是指单位压力变化时引起的岩石视体积的变化率，岩石视体积压缩系数的单位是MPa^-1，通常用符号C_b表示。

岩石孔隙体积压缩系数：是指单位压力变化时引起的岩石孔隙体积的变化率，岩石孔隙体积压缩系数的单位是MPa^-1，通常用符号C_p表示。

已有的通过实验方法测量岩石孔隙体积压缩系数，需要测量不同方向应力变化下的岩石视体积和岩石骨架体积变化量，耗时长，成本高，应用受限。而通过公式计算岩石孔隙体积压缩系数方法大都为经验相关式，在某些特定情况下可能不适用，如Hall提出的经验相关式适用于相对高孔隙度储层，对于相对低孔隙度储层会出现较大误差，Newman给出了一种计算低孔隙度胶结砂岩的经验相关式，对于相对高孔隙度储层会出现较大的误差。因此，有必要利用岩石力学参数，建立计算储层岩石孔隙压缩系数的广义理论公式。

目前在油气藏工程计算中，往往将岩石孔隙体积压缩系数视为一个常数。但在某些特定情况下，将岩石孔隙体积压缩系数作为常数进行计算可能会产生较大的误差，而实验方法和经验方法计算不同有效应力下的岩石孔隙压缩系数都有一定的局限性，因此有必要建立考虑应力变化的岩石孔隙压缩系数计算方法。

所以，本申请提供一种岩石孔隙体积压缩系数的方法，该方法能够基于岩石力学参数，计算岩石孔隙体积压缩系数，并且在计算时考虑有效应力的变化情况，计算所需参数少且参数易获取，并且适用于所有特性的储层，应用更加广泛。

下面以具体地实施例对本申请的技术方案以及本申请的技术方案如何解决上述技术问题进行详细说明。下面这几个具体的实施例可以相互结合，对于相同或相似的概念或过程可能在某些实施例中不再赘述。下面将结合附图，对本申请的实施例进行描述。

参考图1，图1为本发明实施例一提供的计算岩石孔隙体积压缩系数的方法的流程示意图，该方法可以由计算岩石孔隙体积压缩系数的装置执行，该装置可以是电脑、平板电脑和手机等，该方法可以包括如下步骤。

步骤101：统计资源储层的基本参数，该基本参数包括：岩石孔隙度、岩石的泊松比、岩石视体积的杨氏模量、孔隙压力随有效应力变化的衰减系数和有效应力。

该资源储层可以为油气储层，还可以是煤储层等，统计资源储层的基本参数并存储资源储层的基本参数，用于后续岩石孔隙体积压缩系数的计算，该基本参数可以以表一形式存储。

表一

可以理解，表一的各基本参数的取值只是举例说明，并不构成限定。

步骤102：根据岩石孔隙度、岩石的泊松比和岩石视体积的杨氏模量，计算得到岩石骨架体积压缩系数。

该岩石骨架体积压缩系数用于表示单位压力变化时引起的岩石骨架体积的变化率，示例性的可以通过如下公式(1)计算该岩石骨架体积压缩系数。

式中，C_s为岩石骨架体积压缩系数，ν是岩石的泊松比，E为岩石视体积的杨氏模量，φ为岩石孔隙度。

以表一所示参数的取值为例，岩石孔隙度φ值为0.02，岩石的泊松比ν的值为0.3，岩石视体积的杨氏模量E的值为552MPa，将各参数的取值带入公式(1)得到岩石骨架体积压缩系数C_s＝4.348×10^-5MPa^-1。

步骤103：根据岩石的泊松比和岩石视体积的杨氏模量，计算得到岩石视体积压缩系数。

该岩石视体积压缩系数用于表示单位压力变化时引起的岩石视体积的变化率，示例性的可以通过如下公式(2)计算该岩石视体积压缩系数。

式中，C_b为岩石视体积压缩系数，v是岩石的泊松比，E为岩石视体积杨氏模量。

以表一所示参数的取值为例，岩石孔隙度φ值为0.02，岩石视体积的杨氏模量E的值为552MPa，将各参数的取值带入公式(2)得到岩石视体积压缩系数C_b＝8.282×10^-4MPa^-1。

步骤104：根据岩石孔隙度、岩石骨架体积压缩系数和岩石视体积压缩系数，计算得到在零有效应力下岩石孔隙体积压缩系数。

该岩石孔隙体积压缩系数用于表示单位压力变化时引起的岩石孔隙体积的变化率；示例性的可以通过如下公式(3)计算岩石在零有效应力下视体积压缩系数。

式中，C_p0为零有效应力下岩石孔隙体积压缩系数，C_b为岩石视体积压缩系数，C_s为岩石骨架体积压缩系数，φ为岩石孔隙度。

以表一所示参数的取值为例，岩石孔隙度φ值为0.02，岩石骨架体积压缩系数C_s的值为4.348×10^-5MPa^-1，岩石视体积压缩系数C_b的值为8.282×10^-4MPa^-1，将各参数的取值带入公式(3)得到在零有效应力下岩石孔隙体积压缩系数C_p0＝4.354×10^-2MPa^-1。

步骤105：根据孔隙压力随有效应力变化的衰减系数、有效应力和在零有效应力下岩石孔隙体积压缩系数，计算得到在有效应力下的岩石孔隙体积压缩系数。

示例性的可以通过如下公式(4)计算在有效应力下岩石视体积压缩系数。

式中，C_p为有效应力下岩石孔隙体积压缩系数，C_p0为在零有效应力下岩石孔隙体积压缩系数，γ为孔隙压力随有效应力变化的衰减系数，Δσ_eff为有效应力。

表一所示参数的取值为例，效应力Δσ_eff的值为3.5MPa；压力随有效应力变化的衰减系数孔隙压缩系数γ的值为0.15；在零有效应力下初始孔隙压缩系数C_p0的值为4.354×10^-2MPa^-1。将各参数的取值带入公式(4)得到在有效应力为3.5MPa下的岩石孔隙体积压缩系数C_p＝3.387×10^-2MPa^-1。

可以理解，本实施例中涉及到的公式(1)-(4)只是一种实现方式，在计算岩石骨架体积压缩系数、岩石视体积压缩系数、在零有效应力下岩石孔隙体积压缩系数和在有效应力下的岩石孔隙体积压缩系数还可以采用其他公式，或者公式(1)-(4)的简单变形，本实施例不对此进行限制。

本实施例中，通过统计资源储层的基本参数，使用资源储层的基本参数计算岩石孔隙体积压缩系数，资源储层的基本参数易于获取，从而降低了岩石孔隙体积压缩系数计算的难度和成本等，使得该计算方法能够广泛应用。另外，在计算岩石孔隙体积压缩系数时，考虑了有效应力的变化，通过计算在零有效应力下岩石孔隙体积压缩系数，进而根据不同的有效应力和该在零有效应力下岩石孔隙体积压缩系数，可以得到不同有效应力下岩石孔隙体积压缩系数，该计算方式可以适用于任意的储层，使得应用更加广泛。通过统一的方法计算不同储层的岩石孔隙体积压缩系数的方法，使得储层弹性能量评价、油气藏储量计算和压力不稳定分析等资源勘探开发方法更加准确且易于实现，降低了勘探开发成本。

参考图2，图2为本发明实施例二提供的计算岩石孔隙体积压缩系数的装置的结构示意图，该装置100包括统计模块11、第一计算模块12、第二计算模块13、第三计算模块14和第四计算模块15。

统计模块11，用于统计资源储层的基本参数，基本参数包括：岩石孔隙度、岩石的泊松比、岩石视体积的杨氏模量、孔隙压力随有效应力变化的衰减系数孔隙压缩系数和有效应力；

第一计算模12，用于根据岩石孔隙度、岩石的泊松比、岩石视体积的杨氏模量，计算得到岩石骨架体积压缩系数，岩石骨架体积压缩系数用于表示单位压力变化时引起的岩石骨架体积的变化率；

第二计算模块13，用于根据岩石的泊松比和岩石视体积的杨氏模量，计算得到岩石视体积压缩系数，岩石视体积压缩系数用于表示单位压力变化时引起的岩石视体积的变化率；

第三计算模块14，用于根据岩石孔隙度、岩石骨架体积压缩系数和岩石视体积压缩系数，计算得到在零有效应力下岩石孔隙体积压缩系数，岩石孔隙体积压缩系数用于表示单位压力变化时引起的岩石孔隙体积的变化率；

第四计算模15，用于根据孔隙压力随有效应力变化的衰减系数、有效应力、和在零有效应力下初始岩石孔隙压缩系数，计算得到在有效应力下的岩石孔隙体积压缩系数。

本实施例的装置，可用于执行实施例一中的计算岩石孔隙体积压缩系数的方法，具体实现方式和技术效果类似，这里不再赘述

参考图3，图3为本发明实施例三提供的一种电子设备的一种结构示意图，如图3所示，该电子设备200包括：处理器21、存储器22，收发器23，存储器22用于存储指令，收发器23用于和其他设备通信，处理器21用于执行存储器中存储的指令，以使电子设备200执行如上述实施例一中的计算岩石孔隙体积压缩系数的方法步骤，具体实现方式和技术效果类似，这里不再赘述。

本发明实施例四提供一种计算机可读存储介质，计算机可读存储介质中存储有计算机执行指令，计算机执行指令被处理器执行时用于实现如上述实施例一中的计算岩石孔隙体积压缩系数的方法步骤，具体实现方式和技术效果类似，这里不再赘述。

本发明实施例五提供一种计算机程序产品，包括计算机程序，所述计算机程序被处理器执行时，实现如上述实施例一中所述的计算岩石孔隙体积压缩系数的方法步骤，具体实现方式和技术效果类似，这里不再赘述。

本领域技术人员在考虑说明书及实践这里公开的发明后，将容易想到本申请的其它实施方案。本申请旨在涵盖本申请的任何变型、用途或者适应性变化，这些变型、用途或者适应性变化遵循本申请的一般性原理并包括本申请未公开的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的，本申请的真正范围和精神由下面的权利要求书指出。

应当理解的是，本申请并不局限于上面已经描述并在附图中示出的精确结构，并且可以在不脱离其范围进行各种修改和改变。本申请的范围仅由所附的权利要求书来限制。

Claims (10)

1.一种计算岩石孔隙体积压缩系数的方法，其特征在于，包括：

统计资源储层的基本参数，所述基本参数包括：岩石孔隙度、岩石的泊松比、岩石视体积的杨氏模量、孔隙压力随有效应力变化的衰减系数和有效应力；

根据所述岩石孔隙度、所述岩石的泊松比和所述岩石视体积的杨氏模量，计算得到岩石骨架体积压缩系数，所述岩石骨架体积压缩系数用于表示单位压力变化时引起的岩石骨架体积的变化率；

根据所述岩石的泊松比和所述岩石视体积的杨氏模量，计算得到岩石视体积压缩系数，所述岩石视体积压缩系数用于表示单位压力变化时引起的岩石视体积的变化率；

根据所述岩石孔隙度、所述岩石骨架体积压缩系数和所述岩石视体积压缩系数，计算得到在零有效应力下岩石孔隙体积压缩系数，所述岩石孔隙体积压缩系数用于表示单位压力变化时引起的岩石孔隙体积的变化率；

根据所述孔隙压力随有效应力变化的衰减系数、所述有效应力和所述在零有效应力下岩石孔隙体积压缩系数，计算得到在所述有效应力下的岩石孔隙体积压缩系数。

2.根据权利要求1所述的方法，其特征在于，所述根据所述岩石孔隙度、所述岩石的泊松比和岩石视体积的杨氏模量，计算得到岩石骨架体积压缩系数，包括：

根据如下公式计算所述岩石骨架体积压缩系数：

其中，C_s为所述岩石骨架体积压缩系数，ν是所述岩石的泊松比，E为所述岩石视体积的杨氏模量，φ为所述岩石孔隙度。

3.根据权利要求1所述的方法，其特征在于，所述根据所述岩石的泊松比和所述岩石视体积的杨氏模量，计算得到岩石视体积压缩系数，包括：

根据如下公式计算所述岩石视体积压缩系数：

其中，C_b为所述岩石视体积压缩系数，v是所述岩石的泊松比，E为所述岩石视体积杨氏模量。

4.根据权利要求1所述的方法，其特征在于，所述根据所述岩石孔隙度、所述岩石骨架体积压缩系数和所述岩石视体积压缩系数，计算得到在零有效应力下岩石孔隙体积压缩系数，包括：

根据如下公式计算所述在零有效应力下岩石孔隙体积压缩系数：

其中，C_p0为所述零有效应力下岩石孔隙体积压缩系数，C_b为所述岩石视体积压缩系数，C_s为所述岩石骨架体积压缩系数，φ为所述岩石孔隙度。

5.根据权利要求1所述的方法，其特征在于，所述根据所述孔隙压力随有效应力变化的衰减系数、所述有效应力和所述在零有效应力下岩石孔隙体积压缩系数，计算得到在所述有效应力下岩石孔隙体积压缩系数，包括：

根据如下公式计算在所述有效应力下岩石孔隙体积压缩系数：

其中，C_p为在所述有效应力下岩石孔隙体积压缩系数，C_p0为所述零有效应力下岩石孔隙体积压缩系数，γ为所述孔隙压力随有效应力变化的衰减系数，Δσ_eff为所述有效应力。

6.一种计算岩石孔隙体积压缩系数的装置，其特征在于，包括：

统计模块，用于统计资源储层的基本参数，所述基本参数包括：岩石孔隙度、岩石的泊松比、岩石视体积的杨氏模量、孔隙压力随有效应力变化的衰减系数和有效应力；

第一计算模块，用于根据所述岩石孔隙度、所述岩石的泊松比和岩石视体积的杨氏模量，计算得到岩石骨架体积压缩系数，所述岩石骨架体积压缩系数用于表示单位压力变化时引起的岩石骨架体积的变化率；

第二计算模块，用于根据所述岩石的泊松比和所述岩石视体积的杨氏模量，计算得到岩石视体积压缩系数，所述岩石视体积压缩系数用于表示单位压力变化时引起的岩石视体积的变化率；

第三计算模块，用于根据所述岩石孔隙度、所述岩石骨架体积压缩系数和所述岩石视体积压缩系数，计算得到在零有效应力下岩石孔隙体积压缩系数，所述岩石孔隙体积压缩系数用于表示单位压力变化时引起的岩石孔隙体积的变化率；

第四计算模块，用于根据所述孔隙压力随有效应力变化的衰减系数、所述有效应力和所述在零有效应力下岩石孔隙体积压缩系数，计算得到在所述有效应力下岩石孔隙体积压缩系数。

7.根据权利要求6所述的装置，其特征在于，所述第三计算模块具体用于：

根据如下公式计算所述在零有效应力下岩石孔隙体积压缩系数：

其中，C_p0为所述零有效应力下岩石孔隙体积压缩系数，C_b为所述岩石视体积压缩系数，C_s为所述岩石骨架体积压缩系数，φ为所述岩石孔隙度。

8.根据权利要求6所述的装置，其特征在于，所述第四计算模块具体用于：

根据如下公式计算在所述有效应力下的岩石孔隙体积压缩系数：

其中，C_p为在所述有效应力下岩石孔隙体积压缩系数，C_p0为所述零有效应力下岩石孔隙体积压缩系数，γ为所述孔隙压力随有效应力变化的衰减系数，Δσ_eff为所述有效应力。

9.一种电子设备，其特征在于，包括：至少一个处理器和存储器；

所述存储器存储计算机执行指令；

所述至少一个处理器执行所述存储器存储的计算机执行指令，使得所述至少一个处理器执行如权利要求1至5任一项所述的方法。

10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质中存储有计算机执行指令，所述计算机执行指令被处理器执行时用于实现如权利要求1至5任一项所述的方法。

Images