1.一种岩石地层压实系数和有效孔隙体积压缩系数计算方法，其特征在于，包括以下步骤：S1选定研究区域，确定研究对象的储层类型，并采集若干样品；

S2采用PoroPDP‑200型覆压孔隙度渗透率测量仪测试岩石地层压实系数，得到不同有效应力下一系列岩样品的岩石地层压实系数Cpc和孔隙体积离散值；

S3根据岩样品孔隙体积与有效应力的离散测试值，对各个样品孔隙体积和有效应力σm的非线性乘幂关系进行拟合；

S4将孔隙体积转换为孔隙度，并根据S3的结果对岩样孔隙度、有效应力和岩石地层压实系数进行拟合，得到公式：Cpc＝f(φ,σm)

式中，地层压实系数Cpc、σm为有效应力，φ为孔隙度；

S5根据组成岩石样品的各矿物体积百分数fi，采用Voigt‑Reuss‑Hill平均模量模型，计算待预测碳酸盐岩地层的岩石基质压缩系数Cs；

S6将S4中的公式与Biot‑Willis有效应力定理联立求解得到岩石孔隙度、有效应力、地层压实系数；

S7根据岩石地层压实系数和岩石基质压缩系数，计算得到有效孔隙体积压缩系数Cpp，其中用到的公式为：Cpp＝Cpc‑Cs。

2.根据权利要求1所述的岩石地层压实系数和有效孔隙体积压缩系数计算方法，其特征在于：所述步骤S3中，先将每个岩样有效应力测试范围进行扩展，然后通过实测和扩展数据拟合得到岩样孔隙体积、有效应力和岩石地层压实系数的三维曲面关系式，并由决定系数验证该三维曲面关系式的吻合效果。

3.根据权利要求1所述的岩石地层压实系数和有效孔隙体积压缩系数计算方法，其特征在于：所述步骤S5中Voigt‑Reuss‑Hill平均模量模型计算公式如下：

4.根据权利要求1所述的岩石地层压实系数和有效孔隙体积压缩系数计算方法，其特征在于：所述步骤S6中的方程组为：式中，σv为围压，为物理模拟实验给定先验值；p′为孔隙压力，为物理模拟实验给定先验值；α为有效应力系数；Cs为岩石基质压缩系数；

S6‑1：由该公式整理得到非线性方程组；所述未知数包括：地层压实系数、有效应力和孔隙度；

S6‑2：通过粒子群算法搜索一个接近于解的初值；

S6‑3：根据初值，采用非线性约束优化方法得到地层压实系数、有效应力和孔隙度的准确解。

5.根据权利要求1所述的岩石地层压实系数和有效孔隙体积压缩系数计算方法，其特征在于：所述步骤S4中的拟合方式采用非线性最小二乘法。