1.一种基于岩石物理建模的横观各向同性地层地应力预测方法,其特征在于,包括以下步骤:步骤S1、通过测井、录井以及室内测试得到页岩物性参数,所述页岩物性参数包括矿物组成及含量、地层流体组成及含量、孔隙度、饱和度;
步骤S2、建立横观各向同性岩石物理模型;
步骤S3、根据横观各向同性岩石物理模型计算横观各向同性的刚度系数;
步骤S4、根据刚度系数与岩石弹性参数的关系式,计算动态岩石弹性参数;
式中:E1为平行于层理面杨氏模量,GPa;E3为垂直于层理面的杨氏模量,GPa;ν12为平行于层理面,沿最小水平主应力方向的泊松比,无因次;ν13为平行于层理面,沿垂向方向的泊松比,无因次;ν31为垂直于层理面方向的泊松比,无因次;c11、c12、c13、c33为刚度矩阵中的刚度系数,GPa;
步骤S5、利用动态岩石弹性参数,并结合岩石弹性参数的动静态转换关系来计算静态岩石弹性参数;
步骤S6、根据下式计算横观各向同性岩石物理模型中各向异性的Biot系数;
式中: 为干页岩的刚度系数,GPa; 为背景介
质的刚度系数,GPa;α11、α33为平行于层理面和垂直于层理面方向的Biot系数,无因次;
步骤S7、基于Eaton法利用纵波时差来预测地层孔隙压力;
Pp=σv-(σv-pw)(ACn/AC)c
式中:σv为垂向地应力,MPa;pw为地层水静液柱压力,MPa;ACn为该深度点正常压实趋势线纵波时差,μs/ft;AC为实际的纵波时差,μs/ft;Pp为地层孔隙压力,MPa;
步骤S8、利用以下关系式,计算横观各向同性页岩地层最大、最小水平主应力;
式中:σH为最大水平主应力,MPa;σh为最小水平主应力,MPa;Pp为地层孔隙压力,MPa;E1为平行于层理面杨氏模量,GPa;E3为垂直于层理面的杨氏模量,GPa;ν12为平行于层理面,沿最小水平主应力方向的泊松比,无因次;ν13为平行于层理面,沿垂向方向的泊松比,无因次;
σv为垂向地应力,MPa;εH为最大水平主应力方向所对应的最大水平主应变,无因次;εh为最小水平主应力方向所对应的最小水平主应变,无因次。
2.根据权利要求1所述的一种基于岩石物理建模的横观各向同性地层地应力预测方法,其特征在于,所述横观各向同性岩石物理模型的具体建立过程为:步骤S21、利用Hashin-Shtrikman界限计算基质矿物混合物的等效弹性模量;
步骤S22、利用各向同性SCA模型和各向同性DEM模型,将无机质孔隙加入到基质矿物混合物中得到干基质,建立干基质页岩岩石物理模型并计算干基质的等效弹性模量;
步骤S23、利用各向异性SCA模型和各向异性DEM模型计算有机质矿物混合物的等效弹性刚度张量;
步骤S24、利用CL系数描述有机质矿物的成层分布,结合Bond变换得到取向分布的有机质矿物混合物等效弹性刚度张量;
步骤S25、利用各向异性DEM模型,将有机质孔隙加入到取向分布的有机质矿物混合物中得到干有机质,建立干有机质页岩岩石物理模型并计算干有机质的等效弹性刚度张量;
步骤S26、利用各向异性DEM模型,将干基质添加到干有机质中,得到包含空孔隙的干页岩,建立干页岩岩石物理模型并计算干页岩的刚度张量;
步骤S27、根据含气饱和度、含水饱和度,利用Wood公式计算有机质孔隙内气-水混合物的体积模量,并转换成刚度张量;并结合干页岩呈VTI性质的情况,利用Brown-Korringa模型,将混合流体添加到干页岩中得到饱和流体页岩,建立饱和流体页岩岩石物理模型并计算饱和流体页岩的刚度张量。
3.根据权利要求2所述的一种基于岩石物理建模的横观各向同性地层地应力预测方法,其特征在于,所述步骤S21的具体过程是:利用Hashin-Shtrikman界限计算基质矿物混合物体积模量和剪切模量的上下限,并取上下限平均作为基质矿物混合物的体积模量和剪切模量。
4.根据权利要求2所述的一种横观各向同性页岩岩石物理模型的建立方法,其特征在于,所述Hashin-Shtrikman界限计算基质矿物混合物的公式如下:式中:KHS+为基质矿物混合物的等效体积模量上限,GPa;KHS-为基质矿物混合物的等效体积模量下限,GPa;μHS+为基质矿物混合物的等效剪切积模量上限,GPa;μHS-为基质矿物混合物的等效剪切模量下限,GPa。
5.根据权利要求2所述的一种横观各向同性页岩岩石物理模型的建立方法,其特征在于,所述步骤S22中各向同性SCA模型的计算公式如下:式中:vi为第i种材料的体积分数,无因次;P*i为第i种材料的几何因子1,无因次;Q*i为第i种材料的几何因子2; 为等效的体积模量,GPa; 为等效的剪切模量,GPa;
各向同性DEM模型的计算公式如下:
式中:K1、μ1分别为背景介质的体积模量和剪切模量,GPa;K2、μ2分别为包裹体的体积模量和剪切模量,GPa;v为包裹体的体积分数,无因次。
6.根据权利要求2所述的一种横观各向同性页岩岩石物理模型的建立方法,其特征在于,所述步骤S23中各向异性SCA模型的计算公式如下:式中: 为SCA模型等效的刚度张量,GPa; 为第n相物质的Eshelby刚度张量,无因次; 为四阶单位刚度张量,无因; 为第n相物质的刚度张量,GPa;vn为第n相物质的体积分数,无因次;
各向异性DEM模型的计算公式如下:
式中: 为背景介质刚度张量,GPa; 为包裹体刚度张量,GPa; 为包裹体的Eshelby刚度张量,无因次; 为四阶单位刚度张量,无因次;v被添加项的体积,小数。
7.根据权利要求2所述的一种横观各向同性页岩岩石物理模型的建立方法,其特征在于,所述步骤S28中Wood公式如下所示:式中:Kg为气体的体积模量,GPa;Kw为水的体积模量,GPa;Kf为混合流体的等效体积模量,GPa;sg为含气饱和度,无因次;μf为混合流体的等效体积模量,GPa;
Brown-Korringa模型的计算公式如下:
式中: 为饱和岩石的柔度张量,GPa-1; 为干岩石骨架的柔度张量,GPa-1; 为岩石基质的柔度张量,GPa-1;βfl、βgr分别为流体和岩石基质的可压缩系数,GPa-1;φ为孔隙度,无因次。
8.根据权利要求1所述的一种横观各向同性页岩岩石物理模型的建立方法,其特征在于,所述步骤S5中岩石弹性参数的动静态转换关系通过测量的波速计算岩心的动态岩石弹性参数并和三轴压缩实验得到的静态岩石弹性参数建立得到。