1.一种正交各向异性页岩岩石物理模型的建立方法，其特征在于，包括以下步骤：步骤S1、通过测井、录井以及室内测试，得到建模所需的页岩物性参数，所述页岩物性参数包括矿物组成及含量、地层流体组成及含量、孔隙度、饱和度；

步骤S2、利用Hashin-Shtrikman界限计算基质矿物混合物的等效弹性模量；

步骤S3、利用各向同性SCA模型和各向同性DEM模型，将无机质孔隙加入到基质矿物混合物中得到干基质，建立干基质页岩岩石物理模型并计算干基质的等效弹性模量；

步骤S4、利用各向异性SCA模型和各向异性DEM模型计算有机质矿物混合物的等效弹性刚度张量；

步骤S5、利用CL系数描述有机质矿物的成层分布，结合Bond变换得到取向分布的有机质矿物混合物等效弹性刚度张量；

步骤S6、利用各向异性DEM模型，将有机质孔隙加入到取向分布的有机质矿物混合物中得到干有机质，建立干有机质页岩岩石物理模型并计算干有机质的等效弹性刚度张量；

步骤S7、利用各向异性DEM模型，将干基质添加到干有机质中，得到包含空孔隙的干页岩，建立干VTI页岩石物理模型并计算干VTI页岩的刚度张量；

步骤S8、以干VTI页岩基质为背景，利用Schoenberg线性滑动模型的各向异性形式，将垂向分布的裂缝添加到其中，得到干ORT页岩岩石物理模型；

步骤S9、根据含气饱和度、含水饱和度，利用Wood公式计算孔隙内气-水混合物的体积模量，并转换成刚度张量；并结合干页岩呈ORT性质的情况，利用Brown-Korringa模型，将混合流体添加到干页岩中得到饱和流体页岩，建立饱和流体页岩岩石物理模型并计算饱和流体页岩的刚度张量。

2.根据权利要求1所述的一种正交各向异性页岩岩石物理模型的建立方法，其特征在于，所述步骤S2的具体过程是：利用Hashin-Shtrikman界限计算基质矿物混合物体积模量和剪切模量的上下限，并取上下限平均作为基质矿物混合物的体积模量和剪切模量。

3.根据权利要求2所述的一种正交各向异性页岩岩石物理模型的建立方法，其特征在于，所述Hashin-Shtrikman界限计算基质矿物混合物的公式如下：式中：KHS+为基质矿物混合物的等效体积模量上限，GPa；KHS-为基质矿物混合物的等效体积模量下限，GPa；μHS+为基质矿物混合物的等效剪切积模量上限，GPa；μHS-为基质矿物混合物的等效剪切模量下限，GPa。

4.根据权利要求1所述的一种正交各向异性页岩岩石物理模型的建立方法，其特征在于，所述步骤S3中各向同性SCA模型的计算公式如下：式中：vi为第i种材料的体积分数，无因次；P*i为第i种材料的几何因子1，无因次；Q*i为第i种材料的几何因子2； 为等效的体积模量，GPa； 为等效的剪切模量，GPa；

各向同性DEM模型的计算公式如下：

式中：K1、μ1分别为背景介质的体积模量和剪切模量，GPa；K2、μ2分别为包裹体的体积模量和剪切模量，GPa；v为包裹体的体积分数，无因次。

5.根据权利要求1所述的一种正交各向异性页岩岩石物理模型的建立方法，其特征在于，所述步骤S4中各向异性SCA模型的计算公式如下：式中： 为SCA模型等效的刚度张量，GPa； 为第n相物质的Eshelby刚度张量，无因次； 为四阶单位刚度张量，无因次； 为第n相物质的刚度张量，GPa；vn为第n相物质的体积分数，无因次；

各向异性DEM模型的计算公式如下：

式中： 为背景介质刚度张量，GPa； 为包裹体刚度张量，GPa； 为包裹体的Eshelby刚度张量，无因次； 为四阶单位刚度张量，无因次；v被添加项的体积，小数。

6.根据权利要求1所述的一种正交各向异性页岩岩石物理模型的建立方法，其特征在于，所述步骤S8中各向异性Schoenberg线性滑动模型计算公式如下：式中：C为添加裂缝后的干ORT页岩刚度矩阵，GPa；c11b、c12b、c13b、c33b、c44b、c66b为不含垂向裂缝的干VTI页岩刚度系数，GPa；ΔN、ΔV、ΔH描述裂缝特征的弱度，无因次。

7.根据权利要求1所述的一种正交各向异性页岩岩石物理模型的建立方法，其特征在于，所述步骤S9中Wood公式如下所示：式中：Kg为气体的体积模量，GPa；Kw为水的体积模量，GPa；Kf为混合流体的等效体积模量，GPa；sg为含气饱和度，无因次；μf为混合流体的等效体积模量，GPa；

Brown-Korringa模型的计算公式如下：

式中： 为饱和岩石的柔度张量，GPa-1； 为干岩石骨架的柔度张量，GPa-1； 为岩石基质的柔度张量，GPa-1；βfl、βgr分别为流体和岩石基质的可压缩系数，GPa-1；φ为孔隙度，无因次。