[发明专利]一种致密砂岩储层的地震横波速度预测方法及装置

权利要求书

1.一种致密砂岩储层的地震横波速度预测方法，其特征在于，包括：

建立差分有效介质岩石物理模型和多孔隙岩石物理模型；

根据所述差分有效介质岩石物理模型、岩石基质的第一弹性模量获得含裂缝干岩石的第二弹性模量；

根据所述多孔隙岩石物理模型、所述第一弹性模量获得含凸边三角形干岩石的第三弹性模量和含凹边四边形干岩石的第四弹性模量；

获得流体的流体体积模量，并根据所述第二弹性模量、所述第三弹性模量、所述第四弹性模量和所述流体体积模量分别获得饱和流体岩石的第一地震纵波速度、第一地震横波速度、第二地震纵波速度、第二地震横波速度、第三地震纵波速度以及第三地震横波速度；

根据所述第一地震纵波速度、所述第二地震纵波速度和所述第三地震纵波速度确定第一孔隙体积比和第二孔隙体积比；

根据所述第一孔隙体积比、所述第二孔隙体积比、所述第一地震横波速度、所述第二地震横波速度以及所述第三地震横波速度对所述地震横波速度进行预测；

所述第一弹性模量包括第一体积模量和第一剪切模量，所述第二弹性模量包括第二体积模量和第二剪切模量；所述差分有效介质岩石物理模型为：

式中，K^*_DEM(y)为所述第二体积模量；G^*_DEM(y)为所述第二剪切模量；K_i为所述第一体积模量；G_i为所述第一剪切模量；y为干岩石中裂缝的体积分数；为与干岩石中裂缝相关的影响因子；α为第一等效形状因子；

所述第三弹性模量包括第三体积模量和第三剪切模量，所述第四弹性模量包括第四体积模量和第四剪切模量；所述多孔隙岩石物理模型为：

式中，g为第二等效形状因子；为孔隙度；K_dry为体积模量；G_dry为剪切模量；

所述根据所述第二弹性模量、所述第三弹性模量、所述第四弹性模量和所述流体体积模量分别获得饱和流体岩石的第一地震纵波速度、第一地震横波速度、第二地震纵波速度、第二地震横波速度、第三地震纵波速度以及第三地震横波速度包括：

根据所述第二弹性模量和所述流体体积模量获得饱和流体岩石的第五弹性模量；

根据所述第三弹性模量和所述流体体积模量获得饱和流体岩石的第六弹性模量；

根据所述第四弹性模量和所述流体体积模量获得饱和流体岩石的第七弹性模量；

根据所述第五弹性模量、所述第六弹性模量和所述第七弹性模量依次获得所述第一地震纵波速度、所述第一地震横波速度、所述第二地震纵波速度、所述第二地震横波速度、所述第三地震纵波速度以及所述第三地震横波速度；

所述第五弹性模量包括第五体积模量和第五剪切模量，所述第五弹性模量为：

式中，K_sat1为所述第五体积模量；G_sat1为所述第五剪切模量；

所述第六弹性模量包括第六体积模量和第六剪切模量，所述第七弹性模量包括第七体积模量和第七剪切模量；所述第六弹性模量为：

G_sat2＝G_dry1

所述第七弹性模量为：

G_sat3＝G_dry2

式中，K_sat2为所述第六体积模量；K_sat3为所述第七体积模量；G_sat2为所述第六剪切模量；G_sat3为所述第七剪切模量；K_dry1为所述第三体积模量；G_dry1为所述第三剪切模量；K_dry2为所述第四体积模量；G_dry2为所述第四剪切模量；

所述第一地震纵波速度、所述第一地震横波速度、所述第二地震纵波速度、所述第二地震横波速度、所述第三地震纵波速度以及所述第三地震横波速度分别为：

式中，V_{pc_crack}为所述第一地震纵波速度；V_{pc_convex}为所述第二地震纵波速度；V_{pc_concave}为所述第三地震纵波速度；V_{sc_crack}为所述第一地震横波速度；V_{sc_convex}为所述第二地震横波速度；V_{sc_concave}为所述第三地震横波速度；ρ_sat为饱和流体岩石的密度；

所述第一孔隙体积比为：

所述第二孔隙体积比为：

式中，w₁为所述第一孔隙体积比；w₂为所述第二孔隙体积比；V_pm为实测地震纵波速度；

所述地震横波速度为：

V_sc＝V_{sc_convex}-w₁(V_{sc_convex}-V_{sc_concave})

V_sc＝V_{sc_concave}-w₂(V_{sc_concave}-V_{sc_crack})

式中，V_sc为所述地震横波速度。