[发明专利]一种基于记忆模拟退火的基质模量估算方法及电子设备

权利要求书

1.一种基于记忆模拟退火的基质模量估算方法，其特征在于，包括：

通过测井数据获得岩层的饱和度、密度、孔隙度、泥质含量，根据岩层属性查表获得岩层的流体组分模量；

由所述岩层的饱和度、密度、以及流体组分模量得到岩层中流体的模量和密度；由所述岩层的孔隙度和泥质含量得到岩层中干岩石体积和剪切模量；

由所述岩层中流体的模量和密度、干岩石体积的剪切模量、以及基质模量和密度初值通过Gassmann方程计算岩层的体积和剪切模量，其中所述基质模量和密度初值根据岩层属性查表获得；

由所述岩层的泥质含量、基质模量、以及密度初值得到岩石重构密度；

由所述岩层的体积和剪切模量、以及岩石重构密度计算纵横波速度；

根据所述计算纵横波速度和实测纵横波速度建立目标函数，得到基质模量最优值；

所述由所述岩层的孔隙度和泥质含量得到岩层中干岩石体积和剪切模量包括：

由Xu-White模型计算所述岩层中干岩石体积和剪切模量，包括下述公式：

其中，K'为孔隙内含物的体积模量；μ'为孔隙内含物的剪切模量；K为有效体积模量；μ为有效剪切模量；K_m和μ_m分别为骨架体积模量和剪切模量；α_l为孔隙纵横比；上述公式(1)需要φ/α1；

F₆＝1+A[1+g-R(θ+g)]+B(1-θ)(3-4R)

F₉＝A[g(R-1)-Rθ]+Bθ(3-4R)；

所述Gassmann方程为：

其中，K是所求饱和岩石体积模量；φ是孔隙度；K_b是干岩石的体积模量；K_s是岩石基质的体积模量；K_f是孔隙流体的体积模量；

所述由所述岩层的泥质含量、基质模量、以及密度初值得到岩石重构密度包括：由所述岩层的泥质含量、基质模量、以及密度初值建立VRH模型，计算得到岩石骨架模量；由所述岩石骨架模量计算得到岩石重构密度；

所述VRH模型为：

其中，M_h是Voigt上限和Reuss下限的算术平均，M_V是Voigt弹性模量、M_R是Reuss弹性模量，f_i和M_i分别是第i个组成成分的体积分量和模量；

所述由所述岩层的体积和剪切模量、以及岩石重构密度计算纵横波速度包括：由公式(8)和(9)得到所述纵横波速度：

ρ＝∑X_iρ_i (9)

ρ_i为第i种物质的密度，X_i为第i种物质在整体中的百分含量。

2.根据权利要求1所述的基于记忆模拟退火的基质模量估算方法，其特征在于，所述根据所述计算纵横波速度和实测纵横波速度建立目标函数，得到基质模量最优值包括：

根据所述计算纵横波速度和实测纵横波速度建立目标函数，得到基质模量值；

判断所述基质模量值是否满足迭代终止条件；

若不满足，则按照模拟退火优化算法产生迭代后的基质模量和密度初值，并重新执行所述步骤：由所述岩层中流体的模量和密度、干岩石体积的剪切模量、以及基质模量和密度初值通过Gassmann方程计算岩层的体积和剪切模量；

若满足，则所述基质模量值为基质模量最优值。