[发明专利]基于椭球体双重孔隙模型的岩石纵波速度预测方法

权利要求书

1.一种基于椭球体双重孔隙模型的岩石纵波速度预测方法，包括以下步骤：步骤S10，通过获取测井数据和实验观测数据，得到岩石物理参数，并生成岩石的干骨架模型和孔隙流体模型；

步骤S20，以岩石的干骨架模型和孔隙流体模型为基础，建立双重孔隙介质椭球斑块饱和模型，所述双重孔隙介质椭球斑块饱和模型包括三层具有不同孔隙度和不同流体饱和的椭球壳区域模型，每层球壳代表具有不同孔隙特征和不同流体特征的区域，并计算孔隙介质的势能/动能，以及内嵌入体中流体的动能和耗散方程，推导出拉格朗日方程组，并求取纵横波速度；

步骤S30，根据平面波分析方法，获得波动方程的频散关系，并得到纵波速度频散和衰减计算公式。

2.如权利要求1所述的基于椭球体双重孔隙模型的岩石纵波速度预测方法，其特征在于，所述岩石物理参数包括岩石的矿物成分、矿物体积比率、渗透率、孔隙率、泥质含量，以及孔隙中的流体数据包括流体的密度、粘性、弹性模量。

3.如权利要求2所述的基于椭球体双重孔隙模型的岩石纵波速度预测方法，其特征在于，岩石物理参数由目标区域的地质报告、测井资料、岩芯切片、地层的温度、压力、矿化度信息获得。

4.如权利要求3所述的基于椭球体双重孔隙模型的岩石纵波速度预测方法，其特征在于，所述岩石本身的等效弹性模量采用Voigt-Reuss-Hill平均模型计算：

MVRH=12(Σi=1nviMi+1/Σi=1nviMi);]]>

其中：

M_VRH：岩石本身的等效弹性模量；

v_i：第i种矿物的体积率；

M_i：第i种矿物的弹性模量；

n：岩石中矿物的总种类数；

采用Pride半经验性的公式来计算干骨架的等效体积模量和剪切模量:

Kb=1-φ1+cφKs,μb=1-φ1+c'φμs;]]>

K_b,μ_b:岩石干骨架的体积模量和剪切模量；

K_s,μ_s:岩石基质的体积模量和剪切模量；

φ:岩石的孔隙度；

c，c′:为经验性参数，与岩石的固结程度有关。